newton euler equations

%23%20%2F%2F%2F%20script%0A%23%20requires-python%20%3D%20%22%3E%3D3.11%22%0A%23%20dependencies%20%3D%20%5B%22marimo%3E%3D0.23.9%22%5D%0A%23%20%0A%23%20%5Btool.marimo.export%5D%0A%23%20lock_kind%20%3D%20%22observed%22%0A%23%20%2F%2F%2F%0A%0Aimport%20marimo%0A%0A__generated_with%20%3D%20%220.25.2%22%0Aapp%20%3D%20marimo.App()%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20Newton-Euler%20equations%20for%20rigid%20bodies%0A%0A%20%20%20%20%3E%20Renato%20Naville%20Watanabe%2C%20Marcos%20Duarte%3Cbr%3E%0A%20%20%20%20%3E%20%5BLaboratory%20of%20Biomechanics%20and%20Motor%20Control%5D(https%3A%2F%2Fbmclab.pesquisa.ufabc.edu.br%2F)%3Cbr%3E%0A%20%20%20%20%3E%20Federal%20University%20of%20ABC%2C%20Brazil%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20How%20to%20use%20this%20guide%0A%0A%20%20%20%20Start%20from%20linear%20and%20angular%20momentum%2C%20then%20combine%20forces%20into%20equivalent%20force%E2%80%93couple%20systems.%20We%20will%20use%20moments%20about%20the%20center%20of%20mass%20and%20identify%20the%20extra%20care%20needed%20for%20moving%20reference%20points.%20A%20free-body%20diagram%20should%20support%20every%20term%20in%20the%20equations.%0A%0A%20%20%20%20Try%20the%20distributed%20**Exercises**%20before%20opening%20their%20collapsed%20worked%20solutions.%20Then%20use%20the%20**Checkpoint%20questions**%20and%20attempt%20the%20**Problems**%20independently%3B%20their%20panels%20contain%20final%20answers%20only.%20Predict%20the%20effect%20of%20a%20parameter%20change%20before%20running%20it%2C%20and%20check%20units%20and%20limiting%20cases.%0A%0A%20%20%20%20**Opening%20challenge.**%20Can%20two%20force%20systems%20have%20the%20same%20resultant%20force%20but%20produce%20different%20rotational%20motions%20of%20a%20rigid%20body%3F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Python%20setup%0A%0A%20%20%20%20The%20dependency%20declaration%20at%20the%20top%20of%20this%20file%20lists%20the%20packages%20used%20by%20this%20guide.%20Marimo%20runs%20cells%20in%20dependency%20order.%20Run%20the%20notebook%20in%20an%20editor%20to%20use%20its%20displays%20and%20controls%3B%20the%20computational%20examples%20below%20explain%20their%20inputs%20before%20using%20them.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Mechanics%0A%0A%20%20%20%20In%20Mechanics%20we%20are%20interested%20in%20the%20study%20of%20motion%20(including%20deformation)%20and%20forces%20(and%20the%20relation%20between%20them)%20of%20anything%20in%20nature.%0A%0A%20%20%20%20As%20a%20good%20rule%20of%20thumb%2C%20we%20model%20the%20phenomenon%20of%20interest%20as%20simple%20as%20possible%2C%20with%20just%20enough%20complexity%20to%20understand%20the%20phenomenon.%0A%0A%20%20%20%20For%20example%2C%20we%20could%20model%20a%20person%20jumping%20as%20a%20particle%20(the%20center%20of%20gravity%2C%20with%20no%20size)%20moving%20in%20one%20direction%20(the%20vertical)%20if%20all%20we%20want%20is%20to%20estimate%20the%20jump%20height%20and%20relate%20that%20to%20the%20external%20forces%20to%20the%20human%20body.%20So%2C%20mechanics%20of%20a%20particle%20might%20be%20all%20we%20need.%0A%0A%20%20%20%20However%2C%20if%20the%20person%20jumps%20and%20performs%20a%20somersault%2C%20to%20understand%20this%20last%20part%20of%20the%20motion%20we%20have%20to%20model%20the%20human%20body%20as%20one%20of%20more%20objects%20which%20displaces%20and%20rotates%20in%20two%20or%20three%20dimensions.%20In%20this%20case%2C%20we%20would%20need%20what%20is%20called%20mechanics%20of%20rigid%20bodies.%0A%0A%20%20%20%20If%2C%20besides%20the%20gross%20motions%20of%20the%20segments%20of%20the%20body%2C%20we%20are%20interested%20in%20understanding%20the%20deformation%20in%20the%20the%20human%20body%20segments%20and%20tissues%2C%20now%20we%20would%20have%20to%20describe%20the%20mechanical%20behavior%20of%20the%20body%20(e.g.%2C%20how%20it%20deforms)%20under%20the%20action%20of%20forces.%20In%20this%20case%20we%20would%20have%20to%20include%20some%20constitutive%20laws%20describing%20the%20mechanical%20properties%20of%20the%20body.%0A%0A%20%20%20%20In%20the%20chapter%20mechanics%20of%20rigid%20bodies%2C%20the%20body%20deformation%20is%20neglected%2C%20i.e.%2C%20the%20distance%20between%20every%20pair%20of%20points%20in%20the%20body%20is%20considered%20constant.%20Consequently%2C%20the%20position%20and%20orientation%20of%20a%20rigid%20body%20can%20be%20completely%20described%20by%20a%20corresponding%20coordinate%20system%20attached%20to%20it.%0A%0A%20%20%20%20Let's%20review%20some%20Newton's%20laws%20of%20motion%20for%20a%20particle%20and%20then%20extend%20these%20equations%20to%20motion%20of%20rigid%20bodies.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Recapitulation%0A%0A%20%20%20%20%23%23%23%20Newton's%20laws%20of%20motion%0A%0A%20%20%20%20The%20Newton's%20laws%20of%20motion%20describe%20the%20relationship%20between%20the%20forces%20acting%20on%20a%20body%20and%20the%20resultant%20linear%20motion%20due%20to%20those%20forces%3A%0A%0A%20%20%20%20-%20**First%20law**%3A%20An%20object%20will%20remain%20at%20rest%20or%20in%20uniform%20motion%20in%20a%20straight%20line%20unless%20an%20external%20force%20acts%20on%20the%20body.%0A%20%20%20%20-%20**Second%20law**%3A%20The%20acceleration%20of%20an%20object%20is%20directly%20proportional%20to%20the%20net%20force%20acting%20on%20the%20object%20and%20inversely%20proportional%20to%20the%20mass%20of%20the%20object%3A%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BF%7D%7D%20%3D%20m%5Cmathbf%7B%5Cvec%7Ba%7D%7D.%24%3C%2Fspan%3E%0A%20%20%20%20-%20**Third%20law**%3A%20Whenever%20an%20object%20exerts%20a%20force%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BF%7D%7D_1%24%3C%2Fspan%3E%20(action)%20on%20a%20second%20object%2C%20this%20second%20object%20simultaneously%20exerts%20a%20force%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BF%7D%7D_2%24%3C%2Fspan%3E%20on%20the%20first%20object%20with%20the%20same%20magnitude%20but%20opposite%20direction%20(reaction)%3A%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BF%7D%7D_2%20%3D%20%E2%88%92%5Cmathbf%7B%5Cvec%7BF%7D%7D_1.%24%3C%2Fspan%3E%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%20Linear%20momentum%0A%0A%20%20%20%20The%20linear%20momentum%2C%20or%20quantity%20of%20motion%2C%20is%20defined%20as%20the%20product%20between%20mass%20and%20velocity%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BL%7D%7D%20%3D%20m%5Cmathbf%7B%5Cvec%7Bv%7D%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20%23%23%23%20Angular%20momentum%0A%0A%20%20%20%20In%20analogy%20to%20the%20linear%20momentum%2C%20the%20angular%20momentum%20is%20the%20quantity%20of%20movement%20of%20a%20particle%20rotating%20around%20an%20axis%20passing%20through%20any%20point%20O%20at%20a%20distance%20%24%5Cmathbf%7B%5Cvec%7Br%7D%7D%24%20to%20the%20particle%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BH_O%7D%7D%20%3D%20%5Cmathbf%7B%5Cvec%7Br_%7BO%7D%7D%7D%20%5Ctimes%20%5Cmathbf%7B%5Cvec%7BL%7D%7D%0A%20%20%20%20%24%24%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%20Torque%20(moment%20of%20force)%0A%0A%20%20%20%20In%20analogy%20to%20the%20second%20Newton's%20law%20for%20the%20linear%20case%2C%20torque%20or%20moment%20of%20force%20(or%20simply%20moment)%20is%20the%20time%20derivative%20of%20angular%20momentum%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20%5Cfrac%7Bd%5Cmathbf%7B%5Cvec%7BH_O%7D%7D%7D%7Bdt%7D%20%3D%20%5Cfrac%7Bd%7D%7Bdt%7D(%5Cmathbf%7B%5Cmathbf%7B%5Cvec%7Br%7D%7D%20%5Ctimes%20%5Cmathbf%7B%5Cvec%7BL%7D%7D%7D)%20%3D%20%5Cfrac%7Bd%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%7D%7Bdt%7D%20%5Ctimes%20%5Cmathbf%7B%5Cvec%7BL%7D%7D%20%2B%20%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%20%5Ctimes%20%5Cfrac%7Bd%5Cmathbf%7B%5Cvec%7BL%7D%7D%7D%7Bdt%7D%20%3D%20%5Cfrac%7Bd%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%7D%7Bdt%7D%20%5Ctimes%20(m%5Cmathbf%7B%5Cmathbf%7B%5Cvec%7Bv%7D%7D%7D)%20%2B%20%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%20%5Ctimes%20%5Cfrac%7Bd(m%5Cmathbf%7B%5Cvec%7Bv%7D%7D)%7D%7Bdt%7D%20%3D%20%5Cmathbf%7B%5Cvec%7Bv%7D%7D%20%5Ctimes%20(m%5Cmathbf%7B%5Cmathbf%7B%5Cvec%7Bv%7D%7D%7D)%20%2B%20%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%20%5Ctimes%20%5Cfrac%7Bd(m%5Cmathbf%7B%5Cvec%7Bv%7D%7D)%7D%7Bdt%7D%20%3D%200%20%2B%20%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%20%5Ctimes%20%5Cmathbf%7B%5Cvec%7BF%7D%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%20%5Ctimes%20%5Cmathbf%7B%5Cvec%7BF%7D%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20(r_%7BO_x%7D%5C%3A%5Cmathbf%7B%5Chat%7Bi%7D%7D%2Br_%7BO_y%7D%5C%3A%5Cmathbf%7B%5Chat%7Bj%7D%7D%2Br_%7BO_z%7D%5C%3A%5Cmathbf%7B%5Chat%7Bk%7D%7D)%20%5Ctimes%20%20(F_x%5C%3A%5Cmathbf%7B%5Chat%7Bi%7D%7D%2BF_y%5C%3A%5Cmathbf%7B%5Chat%7Bj%7D%7D%2BF_z%5C%3A%5Cmathbf%7B%5Chat%7Bk%7D%7D)%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Where%20the%20symbol%20%24%5Ctimes%24%20stands%20for%20the%20%5Bcross%20product%5D(https%3A%2F%2Fmolab.marimo.io%2Fgithub%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fnotebooks%2Fscalar_vector.py)%20mathematical%20function.%3Cbr%3E%0A%20%20%20%20The%20the%20moment%20of%20force%20can%20be%20calculated%20as%20the%20determinant%20of%20the%20following%20matrix%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20%5Cdet%20%5Cbegin%7Bbmatrix%7D%0A%20%20%20%20%5Cmathbf%7B%5Chat%7Bi%7D%7D%20%26%20%5Cmathbf%7B%5Chat%7Bj%7D%7D%20%26%20%5Cmathbf%7B%5Chat%7Bk%7D%7D%20%5C%5C%0A%20%20%20%20r_%7BO_x%7D%20%26%20r_%7BO_y%7D%20%26%20r_%7BO_z%7D%20%5C%5C%0A%20%20%20%20F_x%20%26%20F_y%20%26%20F_z%0A%20%20%20%20%5Cend%7Bbmatrix%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20(r_%7BO_y%7DF_z-r_%7BO_z%7DF_y)%5Cmathbf%7B%5Chat%7Bi%7D%7D%2B(r_%7BO_z%7DF_x-r_%7BO_x%7DF_z)%5Cmathbf%7B%5Chat%7Bj%7D%7D%2B(r_%7BO_x%7DF_y-r_%7BO_y%7DF_x)%5Cmathbf%7B%5Chat%7Bk%7D%7D%0A%20%20%20%20%24%24%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20The%20magnitude%20of%20moment%20of%20force%20can%20also%20be%20calculated%20by%20the%20geometric%20equivalent%20formula%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%7C%7C%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%7C%7C%20%3D%20%7C%7C%5Cmathbf%7Br_O%7D%20%5Ctimes%20%5Cmathbf%7B%5Cvec%7BF%7D%7D%7C%7C%20%3D%20%7C%7C%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%7C%7C%5C%3A%7C%7C%5Cmathbf%7B%5Cvec%7BF%7D%7D%7C%7C%5Csin(%5Ctheta)%24%3C%2Fspan%3E%0A%0A%20%20%20%20Where%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Ctheta%24%3C%2Fspan%3E%20is%20the%20angle%20between%20the%20vectors%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%24%3C%2Fspan%3E%20and%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BF%7D%7D%24%3Cspan%3E.%0A%0A%20%20%20%20The%20animation%20below%20illustrates%20the%20relationship%20between%20force%2C%20torque%2C%20and%20momentum%20vectors%3A%0A%0A%20%20%20%20%3Ccenter%3E%3Cfigure%20aria-describedby%3D%22figure-caption-1%22%3E%3Cimg%20src%3D%22https%3A%2F%2Fgithub.com%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fimages%2FTorqueAnim.gif%3Fraw%3D1%22%20alt%3D%22Torque%20animation%22%20width%3D%22300%22%2F%3E%3C%2Ffigure%3E%3C%2Fcenter%3E%0A%0A%20%20%20%20%3Cspan%20id%3D%22figure-caption-1%22%20class%3D%22figure-caption%22%20style%3D%22display%3Ablock%3Btext-align%3Acenter%3Bfont-style%3Aitalic%22%3EFigure.%20Relationship%20between%20force%20(%24%5Cmathbf%7B%5Cvec%7BF%7D%7D%24)%2C%20torque%20(%24%5Cmathbf%7B%5Cvec%7BM%7D%7D%24)%2C%20linear%20momentum%20(%24%5Cmathbf%7B%5Cvec%7BL%7D%7D%24)%20and%20angular%20momentum%20(%24%5Cmathbf%7B%5Cvec%7BH%7D%7D%24).%20Adapted%20from%20%3Ca%20href%3D%22https%3A%2F%2Fen.wikipedia.org%2Fwiki%2FFile%3ATorque_animation.gif%22%3EWikipedia%3C%2Fa%3E.%3C%2Fspan%3E%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%20Moment%20of%20inertia%0A%0A%20%20%20%20Let's%20use%20the%20example%20above%2C%20where%20%24%5Cmathbf%7B%5Cvec%7Br_O%7D%7D%24%20and%20%24%5Cmathbf%7B%5Cvec%7BF%7D%7D%24%20are%20orthogonal%20and%20derive%20an%20expression%20for%20the%20magnitude%20of%20these%20quantities%20as%20the%20equivalent%20of%20Newton's%20second%20law%20for%20angular%20motion%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20M_O%20%3D%20r_OF%20%3D%20r_Oma%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Replacing%20the%20linear%20acceleration%20%24a%24%20by%20the%20angular%20acceleration%20%24%5Calpha%24%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20M_O%20%3D%20r_Omr_O%5Calpha%20%3D%20mr_O%5E2%20%5Calpha%0A%20%20%20%20%24%24%0A%0A%20%20%20%20In%20analogy%20to%20Newton's%20second%20law%2C%20where%20the%20constant%20of%20proportionality%20between%20%24a%24%20and%20%24F%24%20is%20called%20inertia%20(mass)%2C%20the%20constant%20of%20proportionality%20between%20%24M_O%24%20and%20%24%5Calpha%24%20is%20called%20rotational%20inertia%20or%20moment%20of%20inertia%2C%20%24I_O%3Dmr_O%5E2%24%20for%20a%20particle%20with%20mass%20%24m%24%20rotating%20at%20a%20distance%20%24r%24%20from%20the%20center%20of%20rotation%20O.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**Exercise%201%20%E2%80%94%20Compute%20a%20vector%20moment.**%20A%20force%20%24(0%2C10%2C0)%24%20N%20acts%20at%20%24(0.2%2C0%2C0)%24%20m%20from%20the%20origin.%20Find%20its%20moment%20and%20identify%20the%20rotation%20axis.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Solution%20to%20Exercise%201%0A%0A%20%20%20%20%24%5Cmathbf%20M%3D%5Cmathbf%20r%5Ctimes%5Cmathbf%20F%3D(0%2C0%2C2)%24%20N%20m.%20The%20positive%20z%20component%20follows%20the%20right-hand%20rule.%20Reversing%20the%20force%20reverses%20the%20moment.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Principle%20of%20transmissibility%20and%20Principle%20of%20moments%0A%0A%20%20%20%20On%20the%20effects%20of%20forces%2C%20there%20are%20two%20important%20principles%3A%0A%0A%20%20%20%20%23%23%23%20Principle%20of%20transmissibility%0A%0A%20%20%20%20%3E%20*For%20rigid%20bodies%20with%20no%20deformation%2C%20an%20external%20force%20can%20be%20applied%20at%20any%20point%20on%20its%20line%20of%20action%20without%20changing%20the%20resultant%20effect%20of%20the%20force.*%0A%0A%20%20%20%20%23%23%23%20Varignon's%20Theorem%20(Principle%20of%20Moments)%0A%0A%20%20%20%20%3E%20*The%20moment%20of%20a%20force%20about%20a%20point%20is%20equal%20to%20the%20sum%20of%20moments%20of%20the%20components%20of%20the%20force%20about%20the%20same%20point.*%3Cbr%3E%0A%20%20%20%20Note%20that%20the%20components%20of%20the%20force%20don't%20need%20to%20be%20orthogonal.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**Exercise%202%20%E2%80%94%20Replace%20a%20displaced%20force.**%20Move%20the%20preceding%20force%20to%20the%20origin.%20What%20must%20be%20added%20to%20preserve%20its%20external%20effect%20on%20a%20rigid%20body%3F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Solution%20to%20Exercise%202%0A%0A%20%20%20%20Add%20a%20free%20couple%20%24(0%2C0%2C2)%24%20N%20m%20at%20the%20same%20time.%20The%20resultant%20force%20alone%20is%20unchanged%2C%20but%20without%20the%20couple%20its%20moment%20about%20the%20origin%20would%20become%20zero.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Equivalent%20systems%0A%0A%0A%20%20%20%20A%20set%20of%20forces%20and%20moments%20is%20considered%20equivalent%20if%20its%20resultant%20force%20and%20sum%20of%20the%20moments%20computed%20relative%20to%20a%20given%20point%20are%20the%20same.%20Normally%2C%20we%20want%20to%20reduce%20all%20the%20forces%20and%20moments%20being%20applied%20to%20a%20body%20into%20a%20single%20force%20and%20a%20single%20moment.%0A%0A%20%20%20%20We%20have%20done%20this%20with%20particles%20for%20the%20resultant%20force.%20The%20resultant%20force%20is%20simply%20the%20sum%20of%20all%20the%20forces%20being%20applied%20to%20the%20body.%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cvec%7B%5Cbf%7BF%7D%7D%20%3D%20%5Csum%5Climits_%7Bi%3D1%7D%5En%20%5Cvec%7B%5Cbf%7BF_i%7D%7D%24%3C%2Fspan%3E%0A%0A%0A%20%20%20%20where%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cvec%7B%5Cbf%7BF_i%7D%7D%24%3C%2Fspan%3E%20is%20each%20force%20applied%20to%20the%20body.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20Similarly%2C%20the%20total%20moment%20applied%20to%20the%20body%20relative%20to%20a%20point%20O%20is%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cvec%7B%5Cbf%7BM_O%7D%7D%20%3D%20%5Csum%5Climits_%7Bi%7D%5Cvec%7B%5Cbf%7Br_%7Bi%2FO%7D%7D%7D%20%5Ctimes%20%5Cvec%7B%5Cbf%7BF_i%7D%7D%20%2B%20%5Csum%5Climits_%7Bi%7D%5Cvec%7B%5Cbf%7BM_%7BF_i%7D%7D%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20where%20%24%5Cvec%7B%5Cbf%7Br_%7Bi%2FO%7D%7D%7D%24%20is%20the%20vector%20from%20the%20point%20O%20to%20the%20point%20where%20the%20force%20%24%5Cvec%7B%5Cbf%7BF_i%7D%7D%24%20is%20being%20applied%20and%20%24%5Cvec%7B%5Cbf%7BM_%7BF_i%7D%7D%7D%24%20are%20the%20free%20moments.%0A%0A%20%20%20%20Free%20moments%20are%20moments%20created%20by%20a%20pair%20of%20forces.%20They%20do%20not%20sum%20um%20to%20the%20resultant%20force%2C%20since%20they%20have%20have%20the%20same%20magnitude%20and%20opposite%20directions%2C%20but%20create%20a%20moment%20because%20there%20is%20a%20distance%20between%20these%20forces.%0A%0A%20%20%20%20%3Ccenter%3E%3Cfigure%3E%3Cimg%20src%3D%22https%3A%2F%2Fgithub.com%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fimages%2FfreeMoment.png%3Fraw%3D1%22%20width%3D500%2F%3E%3C%2Ffigure%3E%3C%2Fcenter%3E%0A%0A%20%20%20%20But%20where%20the%20resultant%20force%20should%20be%20applied%20in%20the%20body%3F%20If%20the%20resultant%20force%20were%20applied%20to%20any%20point%20different%20from%20the%20point%20O%2C%20it%20would%20produce%20an%20additional%20%20moment%20to%20the%20body%20relative%20to%20point%20O.%20So%2C%20the%20resultant%20force%20must%20be%20applied%20to%20the%20point%20O.%0A%0A%20%20%20%20So%2C%20any%20set%20of%20forces%20can%20be%20reduced%20to%20a%20moment%20relative%20to%20a%20chosen%20point%20O%20and%20a%20resultant%20force%20applied%20to%20the%20point%20O.%0A%0A%20%20%20%20To%20compute%20the%20resultant%20force%20and%20moment%20relative%20to%20another%20point%20O'%2C%20the%20new%20moment%20is%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cvec%7B%5Cbf%7BM_%7BO'%7D%7D%7D%20%3D%20%5Cvec%7B%5Cbf%7BM_O%7D%7D%20%2B%20%5Cvec%7B%5Cbf%7Br_%7BO'%2FO%7D%7D%7D%20%5Ctimes%20%5Cvec%7B%5Cbf%7BF%7D%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20And%20the%20resultant%20force%20is%20the%20same.%0A%0A%20%20%20%20It%20is%20worth%20to%20note%20that%20if%20the%20resultant%20force%20%20%24%5Cvec%7B%5Cbf%7BF%7D%7D%24%20is%20zero%2C%20than%20the%20moment%20is%20the%20same%20relative%20to%20any%20point.%0A%0A%20%20%20%20%3Ccenter%3E%3Cfigure%3E%3Cimg%20src%3D%22https%3A%2F%2Fgithub.com%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fimages%2FequivalentSystem.png%3Fraw%3D1%22%20width%3D850%2F%3E%3C%2Ffigure%3E%3C%2Fcenter%3E%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Mechanics%20(dynamics)%20of%20rigid%20bodies%0A%0A%20%20%20%20A%20%5Brigid%20body%5D(https%3A%2F%2Fen.wikipedia.org%2Fwiki%2FRigid_body)%20is%20a%20model%20(an%20idealization)%20for%20a%20body%20in%20which%20deformation%20is%20neglected%2C%20i.e.%2C%20the%20distance%20between%20every%20pair%20of%20points%20in%20the%20body%20is%20considered%20constant.%20This%20definition%20also%20also%20implies%20that%20the%20total%20mass%20of%20a%20rigid%20body%20is%20constant.%0A%0A%20%20%20%20Consequently%2C%20the%20motion%20of%20a%20rigid%20body%20can%20be%20completely%20described%20by%20its%20pose%20(position%20and%20orientation)%20in%20space.%20In%20a%20three-dimensional%20space%2C%20at%20least%20three%20coordinates%20and%20three%20angles%20are%20necessary%20to%20describe%20the%20pose%20of%20the%20rigid%20body%2C%20totalizing%20six%20degrees%20of%20freedom%20for%20a%20rigid%20body.%20This%20also%20implies%20that%20we%20will%20need%20six%20equations%20of%20motion%20for%20these%20components%20to%20describe%20the%20dynamics%20of%20a%20rigid%20body.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**Exercise%203%20%E2%80%94%20Use%20both%20balances.**%20A%203%20kg%20body%20has%20resultant%20force%20%24(6%2C0)%24%20N%20and%20center%20moment%200.4%20N%20m%2C%20with%20center%20inertia%200.2%20kg%20m%C2%B2.%20Find%20its%20planar%20accelerations.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Solution%20to%20Exercise%203%0A%0A%20%20%20%20%24%5Cmathbf%20a_C%3D%5Cmathbf%20F%2Fm%3D(2%2C0)%24%20m%2Fs%C2%B2%20and%20%24%5Calpha%3DM_C%2FI_C%3D2%24%20rad%2Fs%C2%B2.%20Translation%20and%20rotation%20are%20distinct%20parts%20of%20the%20rigid-body%20motion.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Euler's%20laws%20of%20motion%20(for%20a%20rigid%20body)%0A%0A%20%20%20%20Euler's%20laws%20of%20motion%20extend%20Newton's%20laws%20of%20motion%20for%20particles%20for%20the%20motion%20of%20a%20rigid%20body.%0A%0A%20%20%20%20**First%20law**%3A%20The%20linear%20momentum%20of%20a%20body%20is%20equal%20to%20the%20product%20of%20the%20mass%20of%20the%20body%20and%20the%20velocity%20of%20its%20center%20of%20mass%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BL%7D%7D%20%3D%20m%5Cmathbf%7B%5Cvec%7Bv%7D%7D_%7Bcm%7D%24%3C%2Fspan%3E%0A%0A%20%20%20%20And%20calculating%20the%20time%20derivative%20of%20this%20equation%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BF%7D%7D%20%3D%20m%5Cmathbf%7B%5Cvec%7Ba%7D%7D_%7Bcm%7D%24%3C%2Fspan%3E%0A%0A%20%20%20%20**Second%20law**%3A%20The%20rate%20of%20change%20of%20angular%20momentum%20about%20a%20point%20that%20is%20fixed%20in%20an%20inertial%20reference%20frame%20is%20equal%20to%20the%20resultant%20external%20moment%20of%20force%20about%20that%20point%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20%5Cfrac%7Bd%5Cmathbf%7B%5Cvec%7BH_O%7D%7D%7D%7Bdt%7D%24%3C%2Fspan%3E%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%20Derivation%20of%20the%20Euler's%20laws%20of%20motion%0A%0A%20%20%20%20**First%20law**%3A%0A%0A%20%20%20%20The%20sum%20of%20the%20linear%20momentum%20of%20all%20the%20particles%20of%20a%20rigid%20body%20(considering%20the%20body%20as%20a%20discrete%20sum%20of%20elements%2C%20but%20this%20also%20holds%20for%20the%20continuous%20case)%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BL%7D%7D%20%3D%20%5Csum%20m_i%5Cmathbf%7B%5Cvec%7Bv%7D%7D_i%24%3C%2Fspan%3E%0A%0A%20%20%20%20Looking%20at%20the%20definition%20of%20center%20of%20mass%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7Br%7D%7D_%7Bcm%7D%20%3D%20%5Cfrac%7B1%7D%7Bm_%7BB%7D%7D%5Csum%20m_%7Bi%7D%5Cmathbf%7B%5Cvec%7Br%7D%7D_i%20%5Cquad%20%5Ctext%7Bwhere%7D%20%5Cquad%20m_%7BB%7D%20%3D%20%5Csum%20m_%7Bi%7D%24%3Cspan%3E%0A%0A%20%20%20%20By%20differentiation%2C%20the%20velocity%20of%20the%20center%20of%20mass%20is%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7Bv%7D%7D_%7Bcm%7D%20%3D%20%5Cfrac%7B1%7D%7Bm_%7BB%7D%7D%5Csum%20m_%7Bi%7D%5Cmathbf%7B%5Cvec%7Bv%7D%7D_i%24%3C%2Fspan%3E%0A%0A%20%20%20%20And%20finally%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BL%7D%7D%20%3D%20m_%7BB%7D%20%5Cmathbf%7B%5Cvec%7Bv%7D%7D_%7Bcm%7D%20%3D%20m_B%20%5Cmathbf%7B%5Cvec%7Bv%7D%7D_%7Bcm%7D%24%3C%2Fspan%3E%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20We%20can%20get%20the%20second%20equation%20of%20the%20first%20law%20calculating%20the%20time%20derivative%20of%20the%20equation%20above.%3Cbr%3E%0A%20%20%20%20Another%20way%20to%20derive%20this%20second%20equation%20is%20considering%20the%20effects%20of%20all%20forces%20acting%20on%20each%20particle%20of%20the%20rigid%20body%20and%20apply%20Newton's%20second%20law%20to%20them%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Csum%20%5Cmathbf%7B%5Cvec%7BF%7D%7D_i%20%3D%20%5Csum%20m_i%5Cmathbf%7B%5Cvec%7Ba%7D%7D_i%24%3C%2Fspan%3E%0A%0A%20%20%20%20With%20respect%20to%20the%20origin%20of%20these%20forces%2C%20they%20can%20be%20divided%20in%20two%20types%3A%20external%20and%20internal%20forces%20to%20the%20rigid%20body.%20Internal%20forces%20are%20interaction%20forces%20between%20particles%20inside%20the%20body%20and%20because%20of%20Newton's%20third%20law%20(action%20and%20reaction)%20they%20cancel%20each%20other.%20So%2C%20the%20equation%20above%20becomes%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Csum%20%5Cmathbf%7B%5Cvec%7BF%7D%7D_%7Bexternal%7D%20%3D%20%5Csum%20m_i%5Cmathbf%7B%5Cvec%7Ba%7D%7D_i%24%3C%2Fspan%3E%0A%0A%20%20%20%20But%20the%20acceleration%20of%20the%20center%20of%20mass%20is%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7Ba%7D%7D_%7Bcm%7D%20%3D%20%5Cfrac%7B1%7D%7Bm_B%7D%5Csum%20m_%7Bi%7D%5Cmathbf%7B%5Cvec%7Ba%7D%7D_i%24%3C%2Fspan%3E%0A%0A%20%20%20%20And%20finally%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24%5Cmathbf%7B%5Cvec%7BF%7D%7D%20%3D%20%5Csum%20%5Cmathbf%7B%5Cvec%7BF%7D%7D_%7Bexternal%7D%20%3D%20m_B%5Cmathbf%7B%5Cvec%7Ba%7D%7D_%7Bcm%7D%24%3C%2Fspan%3E%0A%0A%20%20%20%20This%20means%20that%20for%20a%20rigid%20body%20the%20internal%20forces%20between%20the%20particles%20of%20the%20body%20do%20not%20contribute%20to%20changing%20the%20total%20momentum%20nor%20changing%20the%20acceleration%20of%20the%20center%20of%20mass.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**Second%20law**%3A%0A%0A%20%20%20%20For%20a%20complete%20derivation%20of%20the%20second%20Euler's%20law%20of%20motion%2C%20see%20for%20example%20Taylor%20(2005)%20or%20%5Bhttp%3A%2F%2Femweb.unl.edu%2FNEGAHBAN%2FEM373%2Fnote19%2Fnote19.htm%5D(http%3A%2F%2Femweb.unl.edu%2FNEGAHBAN%2FEM373%2Fnote19%2Fnote19.htm).%0A%0A%20%20%20%20Let's%20derive%20the%20second%20Euler's%20law%20of%20motion%20for%20a%20simpler%20case%20of%20a%20rigid%20body%20rotating%20in%20a%20plane.%0A%0A%20%20%20%20First%2C%20a%20general%20consideration%20about%20the%20total%20angular%20momentum%20of%20a%20rotting%20rigid%20body%3A%3Cbr%3E%0A%20%20%20%20The%20total%20angular%20momentum%20of%20a%20rigid%20body%20rotating%20around%20a%20point%20%24O%24%20can%20be%20expressed%20as%20the%20angular%20momentum%20of%20the%20body%20center%20of%20mass%20around%20the%20point%20%24O%24%20plus%20the%20sum%20of%20the%20angular%20momentum%20of%20each%20particle%20around%20the%20body%20center%20of%20mass%20(for%20a%20proof%20see%20page%20368%20of%20Taylor%2C%202005)%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BH_O%7D%7D%20%3D%20%5Cmathbf%7B%5Cvec%7Br%7D%7D_%7Bcm%2FO%7D%20%5Ctimes%20m%5Cmathbf%7B%5Cvec%7Bv%7D%7D_%7Bcm%2FO%7D%20%2B%20%5Csum%20%5Cmathbf%7B%5Cvec%7Br%7D%7D_%7Bi%2Fcm%7D%20%5Ctimes%20m_i%5Cmathbf%7B%5Cvec%7Bv%7D%7D_%7Bi%2Fcm%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20For%20a%20two-dimensional%20case%2C%20where%20the%20rigid%20body%20rotates%20around%20its%20own%20center%20of%20mass%20and%20also%20rotates%20around%20another%20parallel%20axis%20(fixed)%2C%20the%20second%20term%20of%20the%20right%20side%20of%20the%20equation%20above%20can%20be%20simplified%20to%20%24%5Csum%20(m_i%5Cmathbf%7Br%7D%5E2_%7Bi%2Fcm%7D)%20%5Cmathbf%7B%5Cvec%7B%5Comega%7D%7D%24%20and%20calculating%20the%20time%20derivative%20of%20the%20whole%20equation%2C%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20%5Cfrac%7Bd%5Cmathbf%7B%5Cvec%7BH_O%7D%7D%7D%7Bdt%7D%20%5Crightarrow%20%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20%5Cfrac%7Bd%5Cleft(%20%5Cmathbf%7B%5Cvec%7Br%7D%7D_%7Bcm%2FO%7D%20%5Ctimes%20m%5Cmathbf%7B%5Cvec%7Bv%7D%7D_%7Bcm%2FO%7D%20%2B%20%5Csum%20%5Cmathbf%7B%5Cvec%7Br%7D%7D_%7Bi%2Fcm%7D%20%5Ctimes%20m_i%5Cmathbf%7B%5Cvec%7Bv%7D%7D_%7Bi%2Fcm%7D%20%5Cright)%7D%7Bdt%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20the%20second%20Euler's%20law%20of%20motion%20simplifies%20to%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BM_O%7D%7D%20%3D%20%5Cmathbf%7B%5Cvec%7Br%7D%7D_%7Bcm%2FO%7D%20%5Ctimes%20m%5Cmathbf%7B%5Cvec%7Ba%7D%7D_%7Bcm%7D%20%2B%20I_%7Bcm%7D%20%5Cmathbf%7B%5Cvec%7B%5Calpha%7D%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20where%20%24%5Cmathbf%7B%5Cvec%7Br%7D%7D_%7Bcm%7D%24%20is%20the%20position%20vector%20of%20the%20center%20of%20mass%20with%20respect%20to%20the%20point%20%24O%24%20about%20which%20moments%20are%20summed%2C%20%24%5Cmathbf%7B%5Cvec%7B%5Calpha%7D%7D%24%20is%20the%20angular%20acceleration%20of%20the%20body%20about%20its%20center%20of%20mass%2C%20and%20%24I_%7Bcm%7D%24%20is%20the%20moment%20of%20inertia%20of%20the%20body%20about%20its%20center%20of%20mass.%0A%0A%20%20%20%20If%20%24d%24%20is%20the%20(shortest)%20distance%20between%20the%20point%20%24O%24%20and%20the%20line%20of%20the%20acceleration%20vector%2C%20then%20the%20equation%20above%20becomes%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7BM%7D%20%3D%20ma_%7Bcm%7Dd%20%2B%20I_%7Bcm%7D%20%5Cmathbf%7B%5Calpha%7D%0A%20%20%20%20%24%24%0A%0A%20%20%20%20Note%20that%20if%20the%20acceleration%20of%20the%20center%20of%20mass%20is%20zero%20or%20the%20sum%20of%20moments%20of%20force%20is%20calculated%20around%20the%20center%20of%20mass%20(then%20%24%5Cmathbf%7Br%7D_%7Bcm%7D%3D0%24)%2C%20this%20case%20of%20rotation%20in%20a%20plane%20simplifies%20to%20the%20well-known%20simple%20equation%3A%0A%0A%20%20%20%20%24%24%0A%20%20%20%20%5Cmathbf%7B%5Cvec%7BM_%7Bcm%7D%7D%7D%20%3D%20I_%7Bcm%7D%20%5Cmathbf%7B%5Cvec%7B%5Calpha%7D%7D%0A%20%20%20%20%24%24%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20*Three-dimensional%20case*%0A%0A%20%20%20%20In%20the%20three-dimensional%20space%2C%20if%20we%20describe%20the%20rotation%20of%20a%20rigid%20body%20using%20a%20rotating%20reference%20frame%20with%20axes%20parallel%20to%20the%20principal%20axes%20of%20inertia%20(referred%20by%20the%20subscripts%201%2C2%2C3)%20of%20the%20body%2C%20the%20Euler's%20second%20law%20becomes%3A%0A%0A%20%20%20%20%3Cspan%20class%3D%22notranslate%22%3E%24M_1%20%3D%20I_1%5Cdot%7B%5Comega_1%7D%20%2B%20(I_3-I_2)%5Comega_2%5Comega_3%24%24M_2%20%3D%20I_2%5Cdot%7B%5Comega_2%7D%20%2B%20(I_1-I_3)%5Comega_3%5Comega_1%24%24M_3%20%3D%20I_3%5Cdot%7B%5Comega_3%7D%20%2B%20(I_2-I_1)%5Comega_1%5Comega_2%24%3C%2Fspan%3E%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**Exercise%204%20%E2%80%94%20Know%20the%203D%20limitation.**%20Why%20is%20%24%5Cmathbf%20M%3DI%5Cboldsymbol%5Calpha%24%20generally%20insufficient%20in%20a%20rotating%20body%20frame%3F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Solution%20to%20Exercise%204%0A%0A%20%20%20%20The%20inertia%20is%20a%20tensor%20and%20the%20rotating%20basis%20contributes%20%24%5Cboldsymbol%5Comega%5Ctimes(%5Cmathbf%20I%5Cboldsymbol%5Comega)%24.%20In%20principal%20body%20axes%2C%20%24%5Cmathbf%20M%3D%5Cmathbf%20I%5Cdot%7B%5Cboldsymbol%5Comega%7D%2B%5Cboldsymbol%5Comega%5Ctimes(%5Cmathbf%20I%5Cboldsymbol%5Comega)%24.%20The%20simpler%20scalar%20form%20is%20appropriate%20for%20fixed-axis%20planar%20rotation.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%0A%20%20%20%20return%20(mo%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Checkpoint%20questions%0A%0A%20%20%20%201.%20What%20must%20match%20for%20two%20force%20systems%20to%20be%20externally%20equivalent%3F%0A%20%20%20%202.%20When%20is%20a%20force%20movable%20without%20adding%20a%20couple%3F%0A%20%20%20%203.%20Why%20is%20the%20center%20of%20mass%20a%20useful%20point%20for%20taking%20moments%3F%0A%20%20%20%204.%20Does%20zero%20resultant%20force%20ensure%20static%20equilibrium%3F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Checkpoint%20answers%0A%0A%20%20%20%201.%20Resultant%20force%20and%20resultant%20moment%20about%20one%20common%20point.%0A%20%20%20%202.%20Along%20its%20line%20of%20action%20for%20an%20ideal%20rigid%20body.%0A%20%20%20%203.%20The%20translational%20contribution%20simplifies%2C%20giving%20the%20center%20angular-momentum%20balance.%0A%20%20%20%204.%20No%3B%20a%20couple%20can%20still%20cause%20angular%20acceleration.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Problems%0A%0A%20%20%20%201.%20(Sample%20problem%202%2F5%2C%20Taylor's%20book)Calculate%20the%20magnitude%20of%20the%20moment%20about%20the%20base%20point%20*O*%20of%20the%20600-N%20force%20in%20five%20different%20ways%20for%20the%20structure%20shown%20below%20(hint%3A%20use%20the%20equation%20for%20torque%20in%20different%20ways%2C%20and%20also%20the%20principles%20of%20moments%20and%20of%20transmissibility)%0A%20%20%20%20%3Ccenter%3E%3Cfigure%3E%3Cimg%20src%3D%22https%3A%2F%2Fgithub.com%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fimages%2Ftorque2.png%3Fraw%3D1%22%20alt%3D%22Torque%22%20width%3D%22250%22%2F%3E%3C%2Ffigure%3E%3C%2Fcenter%3E%0A%0A%20%20%20%202.%20(Example%2017.7%2C%20Hibbeler's%20book)%20A%20uniform%20%2450kg%24%20crate%20rests%20on%20a%20horizontal%20surface%20for%20which%20the%20coefficient%20of%20kinetic%20friction%20is%20%24%5Cmu_k%20%3D%200.2%24.%20Determine%20the%20acceleration%20if%20a%20force%20of%20%24P%20%3D%20600%20N%24%20is%20applied%20to%20the%20crate%20as%20shown%20in%20the%20figure.%3Cbr%3E%0A%20%20%20%20%3Ccenter%3E%3Cfigure%3E%3Cimg%20src%3D%22https%3A%2F%2Fgithub.com%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fimages%2Fhibbeler_17_7.png%3Fraw%3D1%22%20alt%3D%22crate%22%20width%3D%22350%22%2F%3E%3C%2Ffigure%3E%3C%2Fcenter%3E%0A%0A%20%20%20%203.%20(Example%2017.10%2C%20Hibbeler's%20book)%20At%20the%20instant%20shown%20in%20the%20figure%2C%20the%20%2420kg%24%20slender%20rod%20has%20an%20angular%20velocity%20of%20%24%5Comega%20%3D%205%20rad%2Fs%24.%20Determine%20the%20angular%20acceleration%20and%20the%20horizontal%20and%20vertical%20components%20of%20reaction%20of%20the%20pin%20on%20the%20rod%20at%20this%20instant.%3Cbr%3E%0A%20%20%20%20%3Ccenter%3E%3Cfigure%3E%3Cimg%20src%3D%22https%3A%2F%2Fgithub.com%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fimages%2Fhibbeler_17_10.png%3Fraw%3D1%22%20alt%3D%22rod%22%20width%3D%22350%22%2F%3E%3C%2Ffigure%3E%3C%2Fcenter%3E%0A%0A%20%20%20%204.%20(Example%2017.12%2C%20Hibbeler's%20book)%20The%20slender%20rod%20shown%20in%20the%20figure%20has%20a%20mass%20%24m%24%20and%20length%20%24l%24%20and%20is%20released%20from%20rest%20when%20%24%5Ctheta%20%3D%200%5Eo%24.%20Determine%20the%20horizontal%20and%20vertical%20components%20of%20force%20which%20the%20pin%20at%20%24A%24%20exerts%20on%20the%20rod%20at%20the%20instant%20%24%5Ctheta%20%3D%2090%5Eo%24.%3Cbr%3E%0A%20%20%20%20%3Ccenter%3E%3Cfigure%3E%3Cimg%20src%3D%22https%3A%2F%2Fgithub.com%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fimages%2Fhibbeler_17_12.png%3Fraw%3D1%22%20alt%3D%22pendulum%22%20width%3D%22300%22%2F%3E%3C%2Ffigure%3E%3C%2Fcenter%3E%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%2F%2F%2F%20details%20%7C%20Final%20answers%20to%20original%20Problems%0A%0A%20%20%20%201.%20Moment%20magnitude%3A%202609.85%20N%20m%2C%20clockwise.%0A%0A%20%20%20%202.%20Acceleration%3A%2010.038%20m%2Fs%C2%B2%20right%20(%24g%3D9.81%24%20m%2Fs%C2%B2).%20The%20normal%20resultant%20lies%200.46697%20m%20right%20of%20the%20center%2C%20still%20inside%20the%200.5%20m%20half-width%2C%20so%20the%20assumed%20translating%20contact%20is%20consistent.%0A%0A%20%20%20%203.%20The%20figure%20specifies%20%24%5Comega%3D5%24%20rad%2Fs%20(the%20original%20%E2%80%9CS%E2%80%9D%20is%20a%20typographical%20error).%20With%20counterclockwise%20positive%3A%20%24%5Calpha%3D-5.905%24%20rad%2Fs%C2%B2%2C%20%24O_x%3D-750%24%20N%20and%20%24O_y%3D19.05%24%20N.%0A%0A%20%20%20%204.%20At%20the%20downward%20vertical%20position%3A%20%24A_x%3D0%24%2C%20%24A_y%3D5mg%2F2%24%20upward%3B%20the%20rod%20is%20uniform%20and%20released%20from%20the%20horizontal.%0A%20%20%20%20%2F%2F%2F%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20Go%20deeper%0A%0A%20%20%20%20Use%20%5BKineticsFundamentalConcepts%5D(https%3A%2F%2Fmolab.marimo.io%2Fgithub%2FBMClab%2FBMC%2Fblob%2Fmaster%2Fnotebooks%2Fkinetics_fundamental_concepts.py)%20to%20connect%20this%20guide%20with%20the%20surrounding%20mechanics%20and%20modeling%20material.%20Return%20to%20the%20assumptions%20behind%20each%20equation%20when%20applying%20it%20to%20a%20new%20experiment.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%23%20Video%20lectures%20on%20the%20Internet%0A%0A%20%20%20%20%5BMIT%20OpenCourseWare%5D(https%3A%2F%2Focw.mit.edu%2Fcourses%2F2-003sc-engineering-dynamics-fall-2011%2Fpages%2Ffree-body-diagrams-and-fictitious-forces%2F)%20%E2%80%94%20Free-body%20diagrams%3A%20practice%20selecting%20a%20body%20and%20frame%20before%20writing%20force%20and%20moment%20equations.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%23%20References%0A%0A%20%20%20%20-%20Hibbeler%20RC%20(2010)%20%5BEngineering%20Mechanics%20Dynamics%5D(https%3A%2F%2Fdrive.google.com%2Ffile%2Fd%2F1sDLluWCiBCog2C11_Iu1fjv-BtfVUxBU%2Fview).%2012th%20Edition.%20Pearson%20Prentice%20Hall.%0A%0A%20%20%20%20-%20Ruina%20A%2C%20Rudra%20P%20(2019)%20%5BIntroduction%20to%20Statics%20and%20Dynamics%5D(https%3A%2F%2Fruina.tam.cornell.edu%2FBook%2Findex.html).%20Oxford%20University%20Press.%0A%0A%20%20%20%20-%20Taylor%20JR%20(2005)%20%5BClassical%20Mechanics%5D(https%3A%2F%2Fbooks.google.com.br%2Fbooks%3Fid%3DP1kCtNr-pJsC).%20University%20Science%20Books.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A