|
|
| (One intermediate revision by one other user not shown) |
| Line 1: |
Line 1: |
| {{About|the mass moment of inertia of a rotating object|area moment of inertia in beam bending|second moment of area}}
| | A reverse mortgage will give you the answer this retirement dilemma. This option sells your property a piece at a time, in place of all at once. Also, you get to reside in your home. You should use a mortgage calc... <br><br>If you are similar to retired adults, you own a home but have almost no else for retirement. However, if you sell your home, you"ll not have a location to call home! So here"s your problem: you require cash to live on, however the only thing that you own of value could be the area you live. <br><br>A reverse mortgage can give you the answer this pension problem. This program carries your property an item at a time, in place of all at once. Also, you get to reside in your house. You can use a mortgage calculator to find out the monthly cost of home equity loans or refinancing. Also, you may use this [http://browse.Deviantart.com/?q=mortgage+calculator mortgage calculator] to find out just how much your mortgage would cost you altogether. <br><br>First, call a genuine estate agent. They"ll be more than pleased to tell you how much your house could sell for, and how to boost its value. According to your amount of savvy and the time you could commit to it, this could pay off handsomely. The reason is that the total amount that a reverse mortgage can pay you is based on your home"s value. So, when there is an easy way to improve the value of the house, do-it before applying for a reverse mortgage. <br><br>You can use a mortgage calculator to discover if you should obtain a home equity mortgage before you get your reverse mortgage. The mortgage calculator will tell you just how much, in total, a equity loan would charge you for the short time between the repairs and the reverse mortgage. But be careful. Don"t save money remodeling than it will increase your home"s value. Also, if you love anything about your home, do not change it. All things considered, you still get to live in it. <br><br>Ok, since you know just how much your property would sell for, it is time and energy to look into a reverse mortgage loan. You can use a particular mortgage calculator to discover how much each different loan would give you. This mortgage calculator bases its effects on four things: your age, your house"s worth, your house"s area and your bank. More than one company offers a calculator, so it is best to talk with AARP to find out if it is an appropriate program. The mortgage calculator on the web site really is easy, however it is a great place to start. <br><br>But how come it called financing? Because, if you are completed with the house, money is wanted by the lender, not the house. Of course, if the home sells for significantly more than you had been settled, your heirs could get some of it. This is a depth you must workout when you get the loan. Be taught more on this affiliated use with by visiting [http://journals.fotki.com/calmcan110/Tooth-Steven/ Tooth Steven | Fotki.com, photo and video sharing made easy.]. Again, you"ll find mortgage calculator programs to assist you figure this out. If you still have financing on your home, you will have to cover it off before you get your hard earned money. <br><br>Once you"ve done your own research, it"s time to keep in touch with a professional. The real estate agent that you talked to before must be happy to provide you a list of good lenders and mortgage brokers. They will walk you through the process. Read every file. Ask questions about something that you may not understand. And soon, instead of paying a mortgage each month, you will be able to receive a check instead..<br><br>If you are you looking for more on independent health insurance, [http://www.kiwibox.com/berserkfer774/blog just click the following website], have a look at our own web-page. |
| | |
| {{Classical mechanics|cTopic=Fundamental concepts}}
| |
| [[File:Samuel Dixon Niagara.jpg|thumb|right|A tightrope walker uses the moment of inertia of the long rod to help maintain balance. This is Samuel Dixon crossing the Niagara river in 1890.]]
| |
| '''Moment of inertia''' is the mass property of a rigid [[Physical body|body]] that defines the torque needed for a desired change in [[angular velocity]] about an axis of rotation. Moment of inertia depends on the shape of the body and may be different around different axes of rotation. A larger moment of inertia around a given axis requires more torque to increase the rotation, or to stop the rotation, of a body about that axis. Moment of inertia depends on the amount and distribution of its mass, and can be found through the sum of moments of inertia of the masses making up the whole object, under the same conditions. For example, if <math>m_a+m_b=m_c</math>, then <math>I_a+I_b=I_c</math>. In [[classical mechanics]], '''moment of inertia''' may also be called '''mass moment of inertia''', '''rotational inertia''', '''polar moment of inertia''', or the '''angular mass'''.
| |
| | |
| For planar movement of a body, the trajectories of all of its points lie in parallel planes, and the rotation occurs only about an axis perpendicular to this plane. In this case, the body has a single moment of inertia, which is measured around this axis.
| |
| | |
| For spatial movement of a body, the moment of inertia is defined by its symmetric 3x3 inertia matrix. The inertia matrix is often described as a symmetric [[rank (linear algebra)|rank]] two [[tensor]], having six independent components. The inertia matrix includes off-diagonal terms called products of inertia that couple torque around one axis to acceleration about another axis. Each body has a set of mutually perpendicular axes, called [[principal axis (mechanics)|principal axes]], for which the off-diagonal terms of the inertia matrix are zero, and a torque around a principal axis only affects the acceleration about that axis.
| |
| | |
| ==Introduction==
| |
| | |
| [[File:Маховик.jpg|thumb|upright|A [[flywheel]] is a wheel with a large moment of inertia used to smooth out motion in machines. This example is in a Russian museum.]]
| |
| When a body is rotating around an axis, a [[torque]] must be applied to change its [[angular momentum]]. The amount of torque needed for any given change in angular momentum is proportional to the size of that change. Moment of inertia has units of kg·m{{sup|2}} in [[SI]] units and lb{{sub|m}} ft{{sup|2}} in [[United States customary units|US]] units.
| |
| | |
| In 1673 [[Christiaan Huygens]] introduced this parameter in his study of the oscillation of a body hanging from a pivot, known as a [[compound pendulum]].<ref name="mach">{{cite book |last=Mach |first=Ernst |title=The Science of Mechanics |year=1919 |pages=173–187. |url=http://archive.org/details/scienceofmechani005860mbp |accessdate=02/04/2013}}</ref> The term ''moment of inertia'' was introduced by [[Leonhard Euler]] in his book ''Theoria motus corporum solidorum seu rigidorum'' in 1765,<ref name="mach"/><ref name="Euler1730">{{Cite book |last=Euler |first=Leonhard |title=Theoria motus corporum solidorum seu rigidorum: Ex primis nostrae cognitionis principiis stabilita et ad omnes motus, qui in huiusmodi corpora cadere possunt, accommodata |publisher=Cornell University Library |date=1 January 1765 |language=Latin |isbn=978-1-4297-4281-8}}</ref> and it is incorporated into [[Euler's laws#Euler's second law|Euler's second law]].
| |
| | |
| The natural frequency of oscillation of a compound pendulum is obtained from the ratio of the torque imposed by gravity on the mass of the pendulum to the resistance to acceleration defined by the moment of inertia. Comparison of this natural frequency to that of a simple pendulum consisting of a single point of mass provides a mathematical formulation for moment of inertia of an extended body.<ref name="Marion 1995">{{cite book |last=Marion |first=JB |last2=Thornton |first2=ST |year=1995 |title=Classical dynamics of particles & systems |edition=4th |publisher=Thomson |isbn=0-03-097302-3}}</ref><ref name="Symon 1971">{{cite book |last=Symon |first=KR |year=1971 |title=Mechanics |edition=3rd |publisher=Addison-Wesley |isbn=0-201-07392-7}}</ref>
| |
| | |
| Moment of inertia also appears in [[angular momentum|momentum]], [[kinetic energy]], and in [[rigid body dynamics|Newton's laws of motion]] for a rigid body as a physical parameter that combines its shape and mass. There is an interesting difference in the way moment of inertia appears in planar and spatial movement. Planar movement has a single scalar that defines the moment of inertia, while for spatial movement the same calculations yield a 3x3 matrix of moments of inertia, called the inertia matrix or inertia tensor.<ref name="Tenenbaum 2004">{{cite book |last=Tenenbaum |first=RA |year=2004 |title=Fundamentals of Applied Dynamics |publisher=Springer |isbn=0-387-00887-X}}</ref><ref name="Kane">
| |
| {{cite book
| |
| |first=T. R.
| |
| |last=Kane
| |
| |first2=D. A.
| |
| |last2=Levinson
| |
| |url=
| |
| |title=Dynamics, Theory and Applications
| |
| |publisher=McGraw-Hill
| |
| |location=New York
| |
| |year=1985
| |
| |isbn=
| |
| }}</ref>
| |
| | |
| The moment of inertia of a rotating flywheel is used in a machine to resist variations in applied torque in order to smooth its rotational output. The moment of inertia of an airplane about its longitudinal, horizontal and vertical axes determines how steering forces on the control surfaces of its wings, elevators and tail affect the plane in roll, pitch and yaw.
| |
| | |
| === Simple pendulum ===
| |
| A simple pendulum is a point mass suspended by a string so that its movement is constrained to a circle around a pivot point. The mass of a simple pendulum supported by a light string accelerates due to the force of gravity.
| |
| | |
| The moment of inertia of the pendulum about the pivot point is its resistance to movement due to the torque due to gravity. Mathematically, it is the ratio of the torque due to gravity about the pivot of a pendulum to its angular acceleration about that pivot point. For a simple pendulum this is found to be the product of the mass of the particle with the square of its distance to the pivot. This is shown as follows:
| |
| | |
| The force of gravity on the mass of a simple pendulum generates a torque <math>\boldsymbol \tau = \mathbf{r}\times \mathbf{F}</math> around the axis perpendicular to the plane of the pendulum movement. Here '''F''' is the tangential component of the net force on the mass. Associated with this torque is an [[angular acceleration]] <math>\boldsymbol\alpha</math> of the string and mass around this axis. Since the mass is constrained to a circle the tangential acceleration of the mass is <math>\bold a = \boldsymbol\alpha \times \bold r </math>. Since <math>F = ma</math> the torque equation becomes:
| |
| :<math>\boldsymbol\tau = \mathbf{r}\times \mathbf{F} = \mathbf{r}\times (m \boldsymbol\alpha \times \bold r) = (mr^2)\boldsymbol\alpha = I\alpha \bold e,</math>
| |
| where '''e''' is a unit vector perpendicular to the plane of the pendulum. (The second to the last step occurs because of the [[BAC-CAB rule]] using the fact that <math>\boldsymbol\alpha</math> is always perpendicular to '''r'''.) The quantity {{math|''I'' {{=}} ''mr''<sup>2</sup>}} is the ''moment of inertia'' of this single mass around the pivot point.
| |
| | |
| The quantity {{math|''I'' {{=}} ''mr''<sup>2</sup>}} also appears in the [[angular momentum]] of a simple pendulum, which is calculated from the velocity {{math|''v'' {{=}} '''ω'''×'''r'''}} of the pendulum mass around the pivot, where {{math|'''ω'''}} is the [[angular velocity]] of the mass about the pivot point. This angular momentum is given by
| |
| :<math>\mathbf{L} = \mathbf{r}\times(m\mathbf{v}) = (mr^2)\boldsymbol\omega = I\omega\bold e,</math>
| |
| using math similar to that used to derive the previous equation.
| |
| | |
| Similarly, the kinetic energy of the pendulum mass is defined by the velocity of the pendulum around the pivot to yield
| |
| :<math>E_\text{K} = \frac{1}{2}m\mathbf{v}\cdot\mathbf{v} = \frac{1}{2}(mr^2)\omega^2 = \frac{1}{2}I\omega^2.</math>
| |
| | |
| This shows that the quantity {{math|''I'' {{=}} ''mr''<sup>2</sup>}} is how mass combines with the shape of a body to define rotational inertia. The moment of inertia of an arbitrarily shaped body is the sum of the values {{mvar|mr}}{{sup|2}} for all of the elements of mass in the body.
| |
| | |
| ===Compound pendulum===
| |
| [[File:Mendenhall gravimeter pendulums.jpg|thumb|left|Pendulums used in Mendenhall [[gravimeter]] apparatus, from 1897 scientific journal. The portable gravimeter developed in 1890 by Thomas C. Mendenhall provided the most accurate relative measurements of the local gravitational field of the Earth.]]
| |
| A [[compound pendulum]] is a body formed from an assembly of particles or continuous shapes that rotates rigidly around a pivot. Its moments of inertia is the sum the moments of inertia of each of the particles that is composed of.<ref name="B-Paul">{{cite book
| |
| | last =Paul
| |
| | first =Burton
| |
| | title =Kinematics and Dynamics of Planar Machinery
| |
| | publisher =Prentice Hall
| |
| | date =June 1979
| |
| | isbn =978-0135160626
| |
| }}</ref><ref name=Resnick>{{cite book|last=Walker|first=David Halliday, Robert Resnick, Jearl|title=Fundamentals of physics|year=2005|publisher=Wiley|location=Hoboken, NJ|isbn=9780471216438|edition=7th}}</ref>{{rp|395–396}}<ref>{{cite book|last=French|first=A.P.|title=Vibrations and waves|year=1971|publisher=CRC Press|location=Boca Raton, FL|isbn=9780748744473}}</ref>{{rp|51–53}}
| |
| The [[Resonance|natural]] [[Angular frequency|frequency]] (<math> \omega_n</math>) of a compound pendulum depends on its moment of inertia, <math>I_P</math>,
| |
| :<math> \omega_n = \sqrt{\frac{mgr}{I_P}},</math>
| |
| where <math>m</math> is the mass of the object, <math>g</math> is local acceleration of gravity, and <math>r</math> is the distance from the pivot point to the centre of mass of the object. Measuring this frequency of oscillation over small angular displacements provides an effective way of measuring moment of inertia of a body.<ref name ="Uicker"/>{{rp|516–517}}
| |
| | |
| Thus, to determine the moment of inertia of the body, simply suspend it from a convenient pivot point <math>_P</math> so that it swings freely in a plane perpendicular to the direction of the desired moment of inertia, then measure its natural frequency or period of oscillation (<math>t</math>), to obtain
| |
| :<math> I_P = \frac{mgr}{\omega_n^2} = \frac{mgrt^2}{4\pi^2},</math>
| |
| where <math>t</math> is the period (duration) of oscillation (usually averaged over multiple periods).
| |
| | |
| The moment of inertia of the body about its center of mass, <math>I_C</math>, is then calculated using the [[parallel axis theorem]] to be
| |
| :<math> I_C = I_P - mr^2,</math>
| |
| where <math>m</math> is the mass of the body and <math>r</math> is the distance from the pivot point <math>_P</math> to the center of mass <math>_C</math>.
| |
| | |
| Moment of inertia of a body is often defined in terms of its ''[[radius of gyration]]'', which is the radius of a ring of equal mass around the center of mass of a body that has the same moment of inertia. The radius of gyration <math>K</math> is calculated from the body's moment of inertia <math>I_C</math> and mass <math>m</math> as the length,<ref name="Beer"/>{{rp|1296–1297}}
| |
| :<math>K=\sqrt{\frac{I_C}{m}}.</math>
| |
| | |
| ===Center of oscillation===
| |
| A simple pendulum that has the same natural frequency as a compound pendulum defines the length <math>L</math> from the pivot to a point called the [[center of oscillation]] of the compound pendulum. This point also corresponds to the [[center of percussion]]. The length <math>L</math> is determined from the formula,
| |
| :<math> \omega_n = \sqrt{\frac{g}{L}}=\sqrt{\frac{mgr}{I_P}},</math>
| |
| or
| |
| :<math>L=\frac{g}{\omega_n^2} = \frac{I_P}{mr}.</math>
| |
| | |
| The [[seconds pendulum]], which provides the "tick" and "tock" of a grandfather clock, takes one second to swing from side-to-side. This is a period of two seconds, or a natural frequency of π radians/second for the pendulum. In this case, the length <math>L</math> is given by,
| |
| :<math> L = \frac{g}{\omega_n^2}=\frac{9.81 \ \mathrm{m/s^2}}{(3.14 \ \mathrm{rad/s})^2}=0.99 \ \mathrm{m}.</math>
| |
| | |
| Notice that the centre of oscillation of the seconds pendulum must be adjusted to accommodate use in locations with different values for the local acceleration of gravity. [[Kater's pendulum]] is an example of a compound pendulum that is used to measure gravity called a [[gravimeter]].
| |
| | |
| ==Measuring moment of inertia==
| |
| The moment of inertia of complex systems such as a vehicle or airplane around its vertical axis can be measured by suspending the system from three points to form a trifilar pendulum. A trifilar pendulum is a platform supported by three wires designed to oscillate in torsion around its vertical centroidal axis.<ref>H. Williams, [http://www.ima.org.uk/_db/_documents/maths07_williams_huw.pdf Measuring the inertia tensor], presented at the IMA Mathematics 2007 Conference.</ref> The period of oscillation of the trifilar pendulum yields the moment of inertia of the system.<ref>Gracey, William, The experimental determination of the moments of inertia of airplanes by a simplified compound-pendulum method, [http://naca.central.cranfield.ac.uk/reports/1948/naca-tn-1629.pdf NACA Technical Note No. 1629], 1948</ref>
| |
| | |
| == Definition ==
| |
| The moment of inertia, I, is defined as the ratio of an applied torque to the angular acceleration along a principal axis of the object, where then τ, <math>\alpha</math> and I are [[scalar (physics)|scalars]], that is
| |
| | |
| :<math>\tau = I \alpha.</math>
| |
| | |
| An equivalent definition of I uses the angular momentum ''L'' as follows,
| |
| | |
| :<math>L = I \omega,</math>
| |
| | |
| where <math>\omega</math> is the [[angular velocity]] of the object.
| |
| | |
| Apply this definition to a simple pendulum to see that the moment of inertia of the mass ''m'' about the pivot point at a distance ''r'' is
| |
| :<math> I=mr^2.</math>
| |
| This generalizes to define the moment of inertia of a body about an axis '''S''' as the sum of all elemental point masses ''dm'' each multiplied by the square of its perpendicular distance ''r'' to the axis.
| |
| | |
| == Calculating moment of inertia about an axis ==
| |
| [[File:Rolling Racers - Moment of inertia.ogv|thumb|300px|right|Four objects racing down a plane while rolling without slipping. From back to front: spherical shell (red), solid sphere (orange), cylindrical ring (green) and solid cylinder (blue). The time for each object to reach the finishing line depends on their moment of inertia. ([[:File:Rolling Racers - Moment of inertia.ogv|Details]], [[:File:Rolling Racers - Moment of inertia.gif|Animated GIF version]])]]
| |
| | |
| The moment of inertia about an axis of a body is calculated by summing mr<sup>2</sup> for every particle in the body, where ''r'' is the perpendicular distance to the specified axis. In order to see how moment of inertia arises in the study of the movement of an extended body, it is convenient to consider a rigid assembly of point masses. (This equation can be used for axes that are not principal axes provided that it is understood that this does not fully describe the moment of inertia.<ref>In that situation this moment of inertia only describes how a torque applied along that axis causes a rotation about that axis. But, torques not aligned along a principal axis will also cause rotations about other axes.</ref>)
| |
| | |
| Consider the kinetic energy of an assembly of {{mvar|N}} masses {{mvar|m<sub>i</sub>}} that lie at the distances {{mvar|r{{sub|i}} }} from the pivot point ''P'', which is the nearest point on the axis of rotation. It is the sum of the kinetic energy of the individual masses,<ref name="Uicker">
| |
| {{cite book
| |
| |first=John J.
| |
| |last=Uicker
| |
| |first2=Gordon R.
| |
| |last2=Pennock
| |
| |first3=Joseph E.
| |
| |last3=Shigley
| |
| |title=Theory of Machines and Mechanisms
| |
| |edition=4th
| |
| |publisher=Oxford University Press
| |
| |year=2010
| |
| |isbn=978-0195371239
| |
| }}</ref>{{rp|516–517}}<ref name="Beer"/>{{rp|1084–1085}} <ref name="Beer">{{cite book|last=Beer|first=Ferdinand P.|coauthors=E. Russell Johnston, Jr., Phillip J. Cornwell| title=Vector mechanics for engineers: Dynamics|year=2010|publisher=McGraw-Hill|location=Boston|isbn=978-0077295493|edition=9th ed.}}</ref>{{rp|1296–1300}}
| |
| | |
| : <math>E_\text{K} = \sum_{i=1}^N \frac12\,m_i \mathbf{v}_i\cdot\mathbf{v}_i = \sum_{i=1}^N \frac12\,m_i (\omega r_i)^2 = \frac12\, \omega^2 \sum_{i=1}^N m_i r_i^2.</math>
| |
| | |
| This shows that the moment of inertia of the body is the sum of each of the mr<sup>2</sup> terms, that is
| |
| : <math> I_P = \sum_{i=1}^N m_i r_i^2.</math>
| |
| Thus, moment of inertia is a physical property that combines the mass and distribution of the particles around the rotation axis. Notice that rotation about different axes of the same body yield different moments of inertia.
| |
| | |
| The moment of inertia of a continuous body rotating about a specified axis is calculated in the same way, with the summation replaced by the [[multiple integral|integral]],
| |
| : <math> I_P = \int_V \rho(\mathbf{r})\,\mathbf{r}^2 \, dV.</math>
| |
| Again '''r''' is the [[radius vector]] to a point in the body from the specified axis through the pivot ''P'', and {{mvar|ρ}}('''r''') is the mass density at each point '''r'''. The integration is evaluated over the volume {{mvar|V}} of the body. The moment of inertia of a flat surface is similar with the mass density being replaced by its areal mass density with the integral evaluated over its area.
| |
| | |
| '''''Note on second moment of area''''': The moment of inertia of a body moving in a plane and the [[second moment of area]] of a beam's cross-section are often confused. The moment of inertia of body with the shape of the cross-section is the second moment of this area about the z-axis perpendicular to the cross-section, weighted by its density. This is also called the ''polar moment of the area'', and is the sum of the second moments about the x and y axes.<ref>Walter D. Pilkey, [http://books.google.com/books?id=4hEsqvplmFMC&pg=PA437&dq=polar+moment+of+inertia&hl=en&sa=X&ei=1vxkUbj1JIr-rQH-5oC4Bg&ved=0CF4Q6AEwCTgK#v=onepage&q=%22polar%20moment%20of%20inertia%22&f=false Analysis and Design of Elastic Beams: Computational Methods], John Wiley, 2002.</ref> The stresses in a [[beam (structure)|beam]] are calculated using the second moment of the cross-sectional area around either the x-axis or y-axis depending on the load.
| |
| | |
| === Example calculation of moment of inertia ===
| |
| {{main|List of moments of inertia}}
| |
| <!--this does not seem helpful
| |
| [[File:Sodium-hydride-monomer-CRC-MW-3D-balls.png|right|thumb|100px]]
| |
| The moment of inertia of a '''diatomic molecule''' formed by two identical atoms around its center of mass is given by,
| |
| :<math> I = 2mr^2,</math>
| |
| where ''m'' is the mass of each atom and ''r'' is the distance from the center of mass to each atom.
| |
| -->
| |
| <!---
| |
| '''Diatomic molecule''': The moment of inertia of a [[diatomic molecule]] about a perpendicular axis through its center of mass is calculated by first locating the center of mass. For convenience let the x-axis be aligned with the axis of the molecule and place ''m''{{sub|1}} at the origin and ''m''{{sub|2}} at the distance ''d'', then the distance ''r'' to the center of mass is
| |
| :<math> r=\frac{m_2d}{m_1+m_2}.</math>
| |
| | |
| The moment of inertia through the center of mass perpendicular to the molecule is given by
| |
| :<math> I = m_1r^2 + m_2(d-r)^2 = m_1(\frac{m_2d}{m_1+m_2})^2+m_2(\frac{m_1d}{m_1+m_2})^2=\frac{m_1m_2}{m_1+m_2}d^2.</math>
| |
| -->
| |
| [[File:Moment of inertia rod center.svg|thumb|100px|right]]
| |
| | |
| The moment of inertia of a '''compound pendulum''' constructed from a thin disc mounted at the end of a thin rod that oscillates around a pivot at the other end of the rod, begins with the calculation of the moment of inertia of the thin rod and thin disc about their respective centers of mass.<ref name="Beer"/>
| |
| | |
| *The moment of inertia of a '''thin rod''' with constant cross-section {{mvar|s}} and density ρ and with length {{mvar|l}} about a perpendicular axis through its center of mass is determined by integration.<ref name="Beer"/>{{rp|1301}} Align the x-axis with the rod and locate the origin its center of mass at the center of the rod, then
| |
| : <math> I_{C, \text{rod}} = \int \rho\,x^2 dV = \int_{-\ell/2}^{\ell/2} \rho\,x^2 sdx= \rho s\frac{x^3}{3}\bigg|_{-\ell/2}^{\ell/2} =\frac{\rho s}{3} (\ell^3/8 + \ell^3/8) = \frac{1}{12}\, m\ell^2,</math>
| |
| where {{math|''m'' {{=}} ''ρsℓ''}} is the mass of the rod.
| |
| | |
| *The moment of an inertia of a '''thin disc''' of constant thickness {{mvar|s}}, radius {{mvar|R}}, and density {{mvar|ρ}} about an axis through its center and perpendicular to its face (parallel to its axis of [[rotational symmetry]]) is determined by integration.<ref name="Beer"/>{{rp|1301}} Align the z-axis with the axis of the disc and define a volume element as {{math|''dV'' {{=}} ''sr drdθ''}}, then
| |
| :<math> I_{C, \text{disc}}=\int \rho r^2 dV =\int_0^{2\pi} \int_0^R \rho r^2 (s r dr d\theta) = 2\pi \rho s \frac{R^4}{4} = \frac{1}{2}mR^2,</math>
| |
| where {{math|''m'' {{=}} ''πR''{{sup|2}}''ρs''}} is its mass.
| |
| | |
| *The moment of inertia of the compound pendulum is now obtained by adding the moment of inertia of the rod and the disc around the pivot point ''P'' as,
| |
| :<math> I_P = I_{C, \text{rod}} + M_\text{rod}(L/2)^2 + I_{C, \text{disc}} + M_\text{disc}(L+R)^2,</math>
| |
| where ''L'' is the length of the pendulum. Notice that the parallel axis theorem is used to shift the moment of inertia from the center of mass to the pivot point of the pendulum.
| |
| | |
| A [[list of moments of inertia]] formulas for standard body shapes provides a way to obtain the moment of inertial of a complex body as an assembly of simpler shaped bodies. The [[parallel axis theorem]] is used to shift the reference point of the individual bodies to the reference point of the assembly.
| |
| | |
| [[File:Moment of inertia solid sphere.svg|100px|right|thumb]]
| |
| As one more example, consider the moment of inertia of a solid sphere of constant density about an axis through its center of mass. This is determined by summing the moment of inertias of the thin discs that form the sphere. If the surface of the ball is defined by the equation<ref name="Beer"/>{{rp|1301}}
| |
| :<math> x^2+y^2+z^2 =R^2,</math>
| |
| then the radius ''r'' of the disc at the cross-section z along the z-axis is
| |
| :<math>r(z)^2 = x^2+y^2 = R^2 - z^2.</math>
| |
| Therefore, the moment of inertia of the ball is the sum of the moment of inertias of the discs along the z-axis,
| |
| :<math> I_{C, \text{ball}} = \int_{-R}^{R} \frac{\pi \rho}{2} r(z)^4 dz= \int_{-R}^{R} \frac{\pi \rho}{2} (R^2 - z^2)^2 dz
| |
| =\frac{\pi \rho}{2}(R^4z - 2R^2z^3/3+z^5/5)\bigg|_{-R}^{R}
| |
| </math>
| |
| :<math>=\pi \rho(1-2/3+1/5)R^5= \frac{2}{5}mR^2,
| |
| </math>
| |
| where {{math|''m'' {{=}} (4/3)π''R''{{sup|3}}''ρ''}} is the mass of the ball.
| |
| | |
| == Moment of inertia in planar movement of a rigid body ==
| |
| If a [[mechanical system]] is constrained to move parallel to a fixed plane, then the rotation of a body in the system occurs around an axis {{math|''k''}} perpendicular to this plane. In this case, the moment of inertia of the mass in this system is a scalar known as the ''polar moment of inertia''. The definition of the polar moment of inertia can be obtained by considering momentum, kinetic energy and Newton's laws for the planar movement of a rigid system of particles.<ref name="B-Paul"/><ref name="Uicker"/><ref name="Goldstein">{{cite book |last=Goldstein |first=H. |year=1980 |title=Classical Mechanics |edition=2nd |publisher=Addison-Wesley |isbn=0-201-02918-9}}</ref><ref>L. D. Landau and E. M. Lifshitz, [http://archive.org/details/Mechanics_541 Mechanics], Vol 1. 2nd Ed., Pergamon Press, 1969.</ref>
| |
| | |
| If a system of {{mvar|n}} particles, {{math|''P{{sub|i}}'', ''i'' {{=}} 1,...,''n''}}, are assembled into a rigid body, then the momentum of the system can be written in terms of positions relative to a reference point '''R''', and absolute velocities {{math|'''v'''{{sub|i}} }}
| |
| :<math>\Delta \mathbf{r}_i = \mathbf{r}_i - \mathbf{R}, \quad \mathbf{v}_i = \boldsymbol\omega\times(\mathbf{r}_i - \mathbf{R}) + \mathbf{V},</math>
| |
| where ω is the angular velocity of the system and {{math|'''V'''}} is the velocity of {{math|'''R'''}}.
| |
| | |
| For planar movement the angular velocity vector is directed along the unit vector {{mvar|k}} which is perpendicular to the plane of movement. Introduce the unit vectors {{math|'''e'''{{sub|i}} }}from the reference point {{math|'''R'''}} to a point {{math|'''r'''{{sub|i}} }}, and the unit vector {{math|'''t'''{{sub|i}} {{=}} ''k'' × '''e'''{{sub|i}} }} so
| |
| :<math> \Delta r_i\mathbf{e}_i = \mathbf{r}_i-\mathbf{R}, \quad \mathbf{v}_i = \omega \Delta r_i\mathbf{t}_i + \mathbf{V},\quad i=1,\dots, n.</math>
| |
| This defines the relative position vector and the velocity vector for the rigid system of the particles moving in a plane.
| |
| | |
| '''''Note on the cross product''''': When a body moves parallel to a ground plane, the trajectories of all the points in the body lie in planes parallel to this ground plane. This means that any rotation that the body undergoes must be around an axis perpendicular to this plane. Planar movement is often presented as projected onto this ground plane so that the axis of rotation appears as a point. In this case, the angular velocity and angular acceleration of the body are scalars and the fact that they are vectors along the rotation axis is ignored. This is usually preferred for introductions to the topic. But in the case of moment of inertia, the combination of mass and geometry benefits from the geometric properties of the cross product. For this reason, in this section on planar movement the angular velocity and accelerations of the body are vectors perpendicular to the ground plane, and the cross product operations are the same as used for the study of spatial rigid body movement.
| |
| | |
| ===Angular momentum in planar movement===
| |
| [[File:Cup of Russia 2010 - Yuko Kawaguti (2).jpg|thumb|left|upright|A figure skater can reduce her moment of inertia by pulling in her arms, allowing her to spin faster due to [[conservation of angular momentum]].]]
| |
| The angular momentum vector for the planar movement of a rigid system of particles is given by<ref name="B-Paul"/><ref name="Uicker"/>
| |
| :<math>\begin{align} \mathbf{L} & = \sum_{i=1}^n m_i (\mathbf{r}_i-\mathbf{R})\times\mathbf{v}_i \\
| |
| & = \sum_{i=1}^n m_i \Delta r_i\mathbf{e}_i \times(\omega \Delta r_i\mathbf{t}_i + \mathbf{V}) \\
| |
| & = (\sum_{i=1}^n m_i \Delta r_i^2)\omega \vec{k} + (\sum_{i=1}^n m_i \Delta r_i\mathbf{e}_i)\times\mathbf{V}. \\
| |
| \end{align}</math>
| |
| | |
| Use the [[center of mass]] '''C''' as the reference point so
| |
| :<math>\Delta r_i\mathbf{e}_i = \mathbf{r}_i-\mathbf{C},\quad \sum_{i=1}^n m_i\Delta r_i\mathbf{e}_i=0,</math>
| |
| and define the moment of inertia relative to the center of mass I{{sub|C}} as
| |
| :<math> I_C= \sum_{i=1}^n m_i\Delta r_i^2,</math>
| |
| then the equation for angular momentum simplifies to<ref name="Beer"/>{{rp|1028}}
| |
| :<math> \mathbf{L} = I_C \omega \vec{k}.</math>
| |
| The moment of inertia {{mvar|I{{sub|C}} }} about an axis perpendicular to the movement of the rigid system and through the center of mass is known as the ''polar moment of inertia''.
| |
| | |
| For a given amount of angular momentum, a decrease in the moment of inertia results in an increase in the angular velocity. Figure skaters can change their moment of inertia by pulling in their arms. Thus, the angular velocity achieved by a skater with outstretched arms results in a greater angular velocity when the arms are pulled in, because of the reduced moment of inertia.
| |
| | |
| ===Kinetic energy in planar movement===
| |
| [[File:Lever shear flywheel.jpg|thumb|This 1906 rotary shear uses the moment of inertia of two flywheels to store kinetic energy which when released is used to cut metal stock (International Library of Technology, 1906).]]
| |
| The kinetic energy of a rigid system of particles moving in the plane is given by<ref name="B-Paul"/><ref name="Uicker"/>
| |
| :<math>E_\text{K} = \frac{1}{2}\sum_{i=1}^n m_i \mathbf{v}_i\cdot\mathbf{v}_i = \frac{1}{2}\sum_{i=1}^n m_i (\omega \Delta r_i\mathbf{t}_i + \mathbf{V})\cdot(\omega \Delta r_i\mathbf{t}_i + \mathbf{V}).</math>
| |
| This equation expands to yield three terms
| |
| :<math>E_\text{K} = \frac{1}{2}\omega^2\sum_{i=1}^n m_i \Delta r_i^2(\mathbf{t}_i\cdot\mathbf{t}_i) + \omega\mathbf{V}\cdot(\sum_{i=1}^n m_i \Delta r_i\mathbf{t}_i) + \frac{1}{2}(\sum_{i=1}^n m_i) \mathbf{V}\cdot\mathbf{V}.</math>
| |
| Let the reference point be the center of mass '''C''' of the system so the second term becomes zero, and introduce the moment of inertia I{{sub|C}} so the kinetic energy is given by<ref name="Beer"/>{{rp|1084}}
| |
| :<math>E_\text{K} = \frac{1}{2}I_C \omega^2 + \frac{1}{2}M\mathbf{V}\cdot\mathbf{V}.</math>
| |
| The moment of inertia I{{sub|C}} is the ''polar moment of inertia'' of the body.
| |
| | |
| ===Newton's laws for planar movement===
| |
| [[File:Johndeered.jpg|thumb|A 1920's John Deere tractor with the spoked [[flywheel]] on the engine. The large moment of inertia of the flywheel smooths the operation of the tractor]]
| |
| Newton's laws for a rigid system of N particles, {{math|''P{{sub|i}}'', ''i'' {{=}} 1,..., ''N''}}, can be written in terms of a [[resultant force]] and torque at a reference point {{math|'''R'''}}, to yield<ref name="B-Paul"/><ref name="Uicker"/>
| |
| :<math> \mathbf{F} = \sum_{i=1}^N m_i\mathbf{A}_i,\quad \mathbf{T} = \sum_{i=1}^N (\mathbf{r}_i-\mathbf{R})\times (m_i\mathbf{A}_i), </math>
| |
| where {{math|'''r'''{{sub|i}} }} denotes the trajectory of each particle.
| |
| | |
| The [[kinematics]] of a rigid body yields the formula for the acceleration of the particle {{mvar|P{{sub|i}} }} in terms of the position {{math|'''R'''}} and acceleration {{math|'''A'''}} of the reference particle as well as the angular velocity vector {{mvar|ω}} and angular acceleration vector {{mvar|α}} of the rigid system of particles as,
| |
| :<math> \mathbf{A}_i = \boldsymbol\alpha\times(\mathbf{r}_i-\mathbf{R}) + \boldsymbol\omega\times\boldsymbol\omega\times(\mathbf{r}_i-\mathbf{R}) + \mathbf{A}.</math>
| |
| For systems that are constrained to planar movement, the angular velocity and angular acceleration vectors are directed along {{math|''k''}} perpendicular to the plane of movement, which simplifies this acceleration equation. In this case, the acceleration vectors can be simplified by introducing the unit vectors {{math|'''e'''{{sub|i}} }} from the reference point {{math|'''R'''}} to a point {{math|'''r'''{{sub|i}} }} and the unit vectors {{math|'''t'''{{sub|i}} {{=}} ''k'' × '''e'''{{sub|i}} }}, so
| |
| :<math> \mathbf{A}_i = \alpha(\Delta r_i\mathbf{t}_{i}) - \omega^2(\Delta r_i\mathbf{e}_{i}) + \mathbf{A}.</math>
| |
| | |
| This yields the resultant torque on the system as
| |
| :<math>\boldsymbol\tau = \sum_{i=1}^N (m_i\Delta r_i\mathbf{e}_i)\times (\alpha(\Delta r_i\mathbf{t}_{i}) - \omega^2(\Delta r_i\mathbf{e}_{i}) + \mathbf{A}) = (\sum_{i=1}^N m_i\Delta r_i^2)\alpha\vec{k} + (\sum_{i=1}^N m_i\Delta r_i\mathbf{e}_i)\times\mathbf{A},</math>
| |
| where {{math|'''e'''{{sub|i}} × '''e'''{{sub|i}} {{=}} 0}}, and {{math|'''e'''{{sub|i}} × '''t'''{{sub|i}} {{=}} ''k''}} is the unit vector perpendicular to the plane for all of the particles {{mvar|P{{sub|i}} }}.
| |
| | |
| Use the [[center of mass]] {{math|'''C'''}} as the reference point and define the moment of inertia relative to the center of mass {{mvar|I{{sub|C}} }}, then the equation for the resultant torque simplifies to<ref name="Beer"/>{{rp|1029}}
| |
| :<math>\boldsymbol\tau = I_C\alpha\vec{k}.</math>
| |
| The parameter {{mvar|I{{sub|C}} }} is the ''polar moment of inertia'' of the moving body.
| |
| | |
| ==The inertia matrix for spatial movement of a rigid body==
| |
| The scalar moments of inertia appear as elements in a matrix when a system of particles is assembled into a rigid body that moves in three dimensional space. This inertia matrix appears in the calculation of the angular momentum, kinetic energy and resultant torque of the rigid system of particles.<ref name="Marion 1995"/><ref name="Symon 1971"/><ref name="Tenenbaum 2004"/><ref name="Kane"/><ref name="Tsai">L. W. Tsai, Robot Analysis: The mechanics of serial and parallel manipulators, John-Wiley, NY, 1999.</ref>
| |
| | |
| An important application of the inertia matrix and Newton's laws of motion is the analysis of a spinning top. This is discussed in the article on [[gyroscopic precession]]. A more detailed presentation can be found in the article on [[Euler's equations (rigid body dynamics)|Euler's equations of motion]].
| |
| | |
| Let the system of particles {{math|''P{{sub|i}}'', ''i'' {{=}} 1,..., ''n''}} be located at the coordinates {{math|'''r'''{{sub|i}} }} with velocities {{math|'''v'''{{sub|i}} }} relative to a fixed reference frame. For a (possibly moving) reference point {{math|'''R'''}}, the relative positions are
| |
| :<math>\Delta\mathbf{r}_i=\mathbf{r}_i - \mathbf{R}</math>
| |
| and the (absolute) velocities are
| |
| :<math>\mathbf{v}_i = \boldsymbol\omega\times\Delta\mathbf{r}_i + \mathbf{V}_R</math>
| |
| where {{mvar|ω}} is the angular velocity of the system, and {{math|'''V'''}}<sub>''R''</sub> is the velocity of {{math|'''R'''}}.
| |
| | |
| ===Angular momentum===
| |
| If the reference point ''R'' in the assembly, or body, is chosen as the center of mass ''C'', then its angular momentum takes the form,<ref name="Marion 1995"/><ref name="Kane"/>
| |
| :<math> \mathbf{L} = \sum_{i=1}^n m_i\Delta\mathbf{r}_i\times \mathbf{v}_i = \sum_{i=1}^n m_i \Delta\mathbf{r}_i\times(\boldsymbol\omega\times\Delta\mathbf{r}_i),</math>
| |
| where the terms containing {{math|'''V'''}}<sub>''R''</sub> sum to zero by definition of the [[center of mass]].
| |
| | |
| In order to define the inertia matrix, introduce the skew-symmetric matrix [{{mvar|B}}] constructed from a vector {{math|'''b'''}} that performs the cross product operation, such that
| |
| :<math> [B]\mathbf{y} =\mathbf{b}\times\mathbf{y}.</math>
| |
| This matrix [{{mvar|B}}] has the components of {{math|'''b''' {{=}} }} {{mvar|(b{{sub|x}}, b{{sub|y}},b{{sub|z}})}} as its elements, in the form
| |
| :<math> [B] = \begin{bmatrix} 0 & -b_z & b_y \\ b_z & 0 & -b_x \\ -b_y & b_x & 0 \end{bmatrix}.</math>
| |
| | |
| Now construct the skew-symmetric matrix [Δr<sub>i</sub>]= [{{mvar|r{{sub|i}}-C}}] obtained from the relative position vector Δ'''r'''<sub>i</sub>={{math|'''r'''{{sub|i}} - '''C'''}}, and use this skew-symmetric matrix to define,
| |
| :<math> \mathbf{L} = (-\sum_{i=1}^n m_i [\Delta r_i]^2)\boldsymbol\omega = [I_C]\boldsymbol\omega,</math>
| |
| where [{{mvar|I{{sub|C}} }}] defined by
| |
| :<math> [I_C] = -\sum_{i=1}^n m_i[\Delta r_i]^2,</math>
| |
| is the inertia matrix of the rigid system of particles measured relative to the center of mass ''C''.
| |
| | |
| ===Kinetic energy===
| |
| The kinetic energy of a rigid system of particles can be formulated in terms of the [[center of mass]] and a matrix of mass moments of inertia of the system. Let the system of particles {{math|''P{{sub|i}}'', ''i'' {{=}} 1,...,''n''}} be located at the coordinates '''r'''{{sub|i}} with velocities '''v'''{{sub|i}}, then the kinetic energy is<ref name="Marion 1995"/><ref name="Kane"/>
| |
| :<math>E_\text{K} = \frac{1}{2}\sum_{i=1}^n m_i \mathbf{v}_i\cdot\mathbf{v}_i = \frac{1}{2}\sum_{i=1}^n m_i (\boldsymbol\omega\times\Delta\mathbf{r}_i + \mathbf{V}_C)\cdot(\boldsymbol\omega\times\Delta\mathbf{r}_i + \mathbf{V}_C),</math>
| |
| where Δ'''r'''<sub>i</sub>= '''r'''<sub>i</sub>-'''C''' is the position vector of a particle relative to the center of mass.
| |
| | |
| This equation expands to yield three terms
| |
| :<math>E_\text{K} = \frac{1}{2}\sum_{i=1}^n m_i (\boldsymbol\omega\times\Delta\mathbf{r}_i) \cdot(\boldsymbol\omega\times\Delta\mathbf{r}_i)) + \sum_{i=1}^n m_i \mathbf{V}_C\cdot(\boldsymbol\omega\times\Delta\mathbf{r}_i)) + \frac{1}{2}\sum_{i=1}^n m_i \mathbf{V}_C\cdot\mathbf{V}_C.</math>
| |
| The second term in this equation is zero because '''C''' is the center of mass. Introduce the skew-symmetric matrix [Δr{{sub|i}}] so the kinetic energy becomes
| |
| :<math>E_\text{K} = \frac{1}{2}\sum_{i=1}^n m_i ([\Delta r_i]\boldsymbol\omega)\cdot([\Delta r_i]\boldsymbol\omega) + \frac{1}{2}(\sum_{i=1}^n m_i) \mathbf{V}_C\cdot\mathbf{V}_C.</math>
| |
| :<math>E_\text{K} = \frac{1}{2}\sum_{i=1}^n m_i (\boldsymbol\omega^T[\Delta r_i]^T[\Delta r_i]\boldsymbol\omega) + \frac{1}{2}(\sum_{i=1}^n m_i) \mathbf{V}_C\cdot\mathbf{V}_C.</math>
| |
| :<math>E_\text{K} = \frac{1}{2}\boldsymbol\omega\cdot(-\sum_{i=1}^n m_i [\Delta r_i ]^2) \boldsymbol\omega + \frac{1}{2}(\sum_{i=1}^n m_i) \mathbf{V}_C\cdot\mathbf{V}_C.</math>
| |
| | |
| Thus, the kinetic energy of the rigid system of particles is given by
| |
| :<math>E_\text{K} = \frac{1}{2}\boldsymbol\omega\cdot[I_C]\boldsymbol\omega + \frac{1}{2}M\mathbf{V}_C^2.</math>
| |
| where [I{{sub|C}}] is the inertia matrix relative to the center of mass and M is the total mass.
| |
| | |
| ===Resultant torque===
| |
| The inertia matrix appears in the application of Newton's second law to a rigid assembly of particles. The resultant torque on this system is,<ref name="Marion 1995"/><ref name="Kane"/>
| |
| :<math>\boldsymbol\tau = \sum_{i=1}^n (\mathbf{r_i}-\mathbf{R})\times (m_i\mathbf{a}_i),</math>
| |
| where '''a'''{{sub|i}} is the acceleration of the particle P{{sub|i}}. The [[kinematics]] of a rigid body yields the formula for the acceleration of the particle P{{sub|i}} in terms of the position '''R''' and acceleration '''A''' of the reference point, as well as the angular velocity vector ω and angular acceleration vector α of the rigid system as,
| |
| :<math> \mathbf{a}_i = \boldsymbol\alpha\times(\mathbf{r}_i-\mathbf{R}) + \boldsymbol\omega\times\boldsymbol\omega\times(\mathbf{r}_i-\mathbf{R}) + \mathbf{A}_R.</math>
| |
| | |
| Use the center of mass '''C''' as the reference point, and introduce the skew-symmetric matrix [Δr<sub>i</sub>]=[r{{sub|i}}-C] to represent the cross product ('''r'''{{sub|i}} - '''C''')x, in order to obtain
| |
| :<math>\boldsymbol\tau = (-\sum_{i=1}^n m_i[\Delta r_i]^2)\boldsymbol\alpha + \boldsymbol\omega\times(-\sum_{i=1}^n m_i [\Delta r_i]^2)\boldsymbol\omega</math>
| |
| This calculation uses the identity
| |
| :<math> \Delta\mathbf{r}_i\times(\boldsymbol\omega\times(\boldsymbol\omega\times\Delta\mathbf{r}_i )) + \boldsymbol\omega\times(\Delta\mathbf{r}_i\times(\Delta\mathbf{r}_i\times\boldsymbol\omega))=0,</math>
| |
| obtained from the Jacobi identity for the triple [[cross product]].
| |
| | |
| Thus, the resultant torque on the rigid system of particles is given by
| |
| :<math>\boldsymbol\tau = [I_C]\boldsymbol\alpha + \boldsymbol\omega\times[I_C]\boldsymbol\omega,</math>
| |
| where [I{{sub|C}}] is the inertia matrix relative to the center of mass.
| |
| | |
| ===Parallel axis theorem===
| |
| {{Main|Parallel axis theorem}}
| |
| The inertia matrix of a body depends on the choice of the reference point. There is a useful relationship between the inertia matrix relative to the center of mass '''C''' and the inertia matrix relative to another point '''R'''. This relationship is called the parallel axis theorem.<ref name="Marion 1995"/><ref name="Kane"/>
| |
| | |
| Consider the inertia matrix [I{{sub|R}}] obtained for a rigid system of particles measured relative to a reference point '''R''', given by
| |
| :<math> [I_R] = -\sum_{i=1}^n m_i[r_i-R]^2.</math>
| |
| Let '''C''' be the center of mass of the rigid system, then
| |
| :<math> \mathbf{R} = (\mathbf{R}-\mathbf{C}) + \mathbf{C} = \mathbf{d} + \mathbf{C},</math>
| |
| where '''d''' is the vector from the center of mass '''C''' to the reference point '''R'''. Use this equation to compute the inertia matrix,
| |
| :<math> [I_R] = -\sum_{i=1}^n m_i[r_i - C - d]^2.</math>
| |
| | |
| Expand this equation to obtain
| |
| :<math> [I_R] = (-\sum_{i=1}^n m_i [r_i - C]^2) + (\sum_{i=1}^n m_i[r_i - C])[d] + [d](\sum_{i=1}^n m_i[r_i - C]) + (-\sum_{i=1}^n m_i)[d][d].</math>
| |
| The first term is the inertia matrix [I{{sub|C}}] relative to the center of mass. The second and third terms are zero by definition of the center of mass '''C'''. And the last term is the total mass of the system multiplied by the square of the skew-symmetric matrix [d] constructed from '''d'''.
| |
| | |
| The result is the parallel axis theorem,
| |
| :<math> [I_R] = [I_C] - M[d]^2,</math>
| |
| where '''d''' is the vector from the center of mass '''C''' to the reference point '''R'''.
| |
| | |
| '''Note on the minus sign''': By using the skew symmetric matrix of position vectors relative to the reference point, the inertia matrix of each particle has the form ''-m''[r]<sup>2</sup>, which is similar to the ''mr''<sup>2</sup> that appears in planar movement. However, to make this to work out correctly a minus sign is needed. This minus sign can be absorbed into the term ''m''[r]<sup>T</sup>[r], if desired, by using the skew-symmetry property of [r].
| |
| | |
| ==The inertia matrix and the scalar moment of inertia around an arbitrary axis==
| |
| The relationship between the inertia matrix of a rigid body and the scalar moment of inertia of the same body about a specified axis is important and rarely presented in detail. The following calculation expands the derivation presented by Kane and Levinson.<ref name="Kane"/>
| |
| | |
| Let a rigid assembly of rigid system of {{mvar|N}} particles, {{math|''P{{sub|i}}'', ''i'' {{=}} 1,...,''N''}}, have coordinates '''r'''{{sub|i}}. Choose '''R''' as a reference point and compute the moment of inertia around an axis L defined by the unit vector '''S''' through the reference point '''R'''. The moment of inertia of the system around this line L='''R'''+t'''S''' is computed by determining the perpendicular vector from this axis to the particle P{{sub|i}} given by
| |
| :<math> \Delta\mathbf{r}_i ^\perp= (\mathbf{r}_i-\mathbf{R}) - (\mathbf{S}\cdot(\mathbf{r}_i-\mathbf{R}))\mathbf{S} = [[I]-[\mathbf{S}\mathbf{S}^T]](\Delta\mathbf{r}_i),</math>
| |
| where [I] is the identity matrix and ['''S''' '''S'''{{sup|T}}] is the outer product matrix formed from the unit vector '''S''' along the line L.
| |
| | |
| In order to relate this scalar moment of inertia to the inertia matrix of the body, introduce the skew-symmetric matrix [S] such that [S]'''y'''='''S''' x '''y''', then we have the identity
| |
| :<math> -[S]^2 = [I]-[\mathbf{S}\mathbf{S}^T],</math>
| |
| which relies on the fact that '''S''' is a unit vector.
| |
| | |
| The magnitude squared of the perpendicular vector is
| |
| :<math> |\Delta\mathbf{r}_i^\perp|^2 = (-[S]^2(\Delta\mathbf{r}_i)) \cdot (-[S]^2(\Delta\mathbf{r}_i)) = -\mathbf{S}\cdot[\Delta r_i][\Delta r_i]\mathbf{S}.</math>
| |
| The simplification of this equation uses the identity
| |
| :<math> (\mathbf{S}\times(\mathbf{S}\times(\Delta\mathbf{r}_i))) \cdot \mathbf{S}\times(\mathbf{S}\times(\Delta\mathbf{r}_i)) = (\mathbf{S}\times(\mathbf{S}\times(\Delta\mathbf{r}_i)))\times\mathbf{S}\cdot (\mathbf{S}\times(\Delta\mathbf{r}_i)),</math>
| |
| where the dot and the cross products have been interchanged. Expand the cross products to compute
| |
| :<math> -(\Delta\mathbf{r}_i)\times\mathbf{S}\cdot(\mathbf{S}\times(\Delta\mathbf{r}_i)=-\mathbf{S}\cdot[\Delta r_i][\Delta r_i]\mathbf{S},</math>
| |
| where [Δr{{sub|i}}] is the skew symmetric matrix obtained from the vector Δ'''r'''='''r'''{{sub|i}}-'''R'''.
| |
| | |
| Thus, the moment of inertia around the line L through '''R''' in the direction '''S''' is obtained from the calculation
| |
| :<math> I_L = \sum_{i=1}^N m_i |\Delta\mathbf{r}_i^\perp|^2= -\sum_{i=1}^N m_i \mathbf{S}\cdot[\Delta r_i]^2\mathbf{S},</math>
| |
| or
| |
| :<math>I_L = \mathbf{S}\cdot(-\sum_{i=1}^N m_i [\Delta r_i]^2)\mathbf{S}=\mathbf{S}\cdot[I_R]\mathbf{S}=\mathbf{S}^T[I_R]\mathbf{S},</math>
| |
| where [I{{sub|R}}] is the moment of inertia matrix of the system relative to the reference point '''R'''.
| |
| | |
| This shows that the inertia matrix can be used to calculate the moment of inertia of a body around any specified rotation axis in the body.
| |
| | |
| == The inertia tensor ==
| |
| The inertia matrix is often described as the inertia tensor, which consists of the same moments of inertia and products of inertia about the three coordinate axes.<ref name="Kane"/><ref name="Goldstein"/> The inertia tensor is constructed from the nine component tensors, (the symbol <math>\otimes</math> is the [[tensor product]])
| |
| :<math> \mathbf{e}_i\otimes\mathbf{e}_j,\quad i,j=1,2,3,</math>
| |
| where '''e'''{{sub|i}}, i=1,2,3 are the three orthogonal [[unit vector]]s defining the inertial frame in which the body moves. Using this basis the inertia tensor is given by
| |
| :<math>\mathbf{I}=\sum_{i=1}^3\sum_{j=1}^3 I_{ij}\mathbf{e}_i\otimes\mathbf{e}_j.</math>
| |
| This tensor is of degree two because the component tensors are each constructed from two basis vectors. In this form the inertia tensor is also called the ''inertia binor''.
| |
| | |
| For a rigid system of particles {{math|''P{{sub|k}}'', ''k'' {{=}} 1,...,''N''}} each of mass m{{sub|k}} with position coordinates '''r'''{{sub|k}}=(x{{sub|k}}, y{{sub|k}}, z{{sub|k}}), the inertia tensor is given by
| |
| :<math>\mathbf{I} =\sum_{k=1}^Nm_k((\mathbf{r}_k\cdot\mathbf{r}_k)\mathbf{E}-\mathbf{r}_k\otimes\mathbf{r}_k),</math>
| |
| where '''E''' is the identity tensor
| |
| :<math>\mathbf{E}= \mathbf{e}_1\otimes\mathbf{e}_1+ \mathbf{e}_2\otimes\mathbf{e}_2+\mathbf{e}_3\otimes\mathbf{e}_3.</math>
| |
| <!--In this case, the components of the inertia tensor are given by
| |
| : <math>I_{11} = I_{xx} = \sum_{k=1}^{N} m_{k} (y_{k}^{2}+z_{k}^{2}),\,\! </math>
| |
| : <math>I_{22} = I_{yy} = \sum_{k=1}^{N} m_{k} (x_{k}^{2}+z_{k}^{2}),\,\!</math>
| |
| : <math>I_{33} = I_{zz} =\sum_{k=1}^{N} m_{k} (x_{k}^{2}+y_{k}^{2}),\,\!</math>
| |
| : <math>I_{12} = I_{21} = I_{xy} = -\sum_{k=1}^{N} m_{k} x_{k} y_{k},\,\!</math>
| |
| : <math>I_{13} = I_{31} = I_{xz} = -\sum_{k=1}^{N} m_{k} x_{k} z_{k},\,\!</math>
| |
| : <math>I_{23} = I_{32} = I_{yz} = -\sum_{k=1}^{N} m_{k} y_{k} z_{k}.\,\!</math>-->
| |
| | |
| The inertia tensor for a continuous body is given by
| |
| : <math>\mathbf{I}=\int_V \rho(\mathbf{r}) \left( \left( \mathbf{r} \cdot \mathbf{r} \right) \mathbf{E} - \mathbf{r}\otimes \mathbf{r}\right)\, dV,</math>
| |
| where '''r''' defines the coordinates of a point in the body and ρ('''r''') is the mass density at that point. The integral is taken over the volume ''V'' of the body. The inertia tensor is symmetric because I{{sub|ij}}= I{{sub|ji}}.
| |
| | |
| Alternatively it can also be written in terms of the [[Hat operator|hat]] operator as:
| |
| : <math>\mathbf{I}=\int_V \rho(\mathbf{r}) (\hat r)^2 \, dV,</math>
| |
| | |
| The inertia tensor can be used in the same way as the inertia matrix to compute the scalar moment of inertia about an arbitrary axis in the direction '''n''',
| |
| :<math> I_n = \mathbf{n}\cdot\mathbf{I}\cdot\mathbf{n},</math>
| |
| where the dot product is taken with the corresponding elements in the component tensors. A product of inertia term such as I{{sub|12}} is obtained by the computation
| |
| :<math>I_{12} = \mathbf{e}_1\cdot\mathbf{I}\cdot\mathbf{e}_2,</math>
| |
| and can be interpreted as the moment of inertia around the x-axis when the object rotates around the y-axis.
| |
| | |
| The components of tensors of degree two can be assembled into a matrix. For the inertia tensor this matrix is given by,
| |
| : <math>[I] = \begin{bmatrix}
| |
| I_{11} & I_{12} & I_{13} \\
| |
| I_{21} & I_{22} & I_{23} \\
| |
| I_{31} & I_{32} & I_{33}
| |
| \end{bmatrix}=\begin{bmatrix}
| |
| I_{xx} & I_{xy} & I_{xz} \\
| |
| I_{xy} & I_{yy} & I_{yz} \\
| |
| I_{xz} & I_{yz} & I_{zz}
| |
| \end{bmatrix}.
| |
| </math>
| |
| It is common in rigid body mechanics to use notation that explicitly identifies the x, y, and z axes, such as I{{sub|xx}} and I{{sub|xy}}, for the components of the inertia tensor.
| |
| | |
| ===Identities for a skew-symmetric matrix===
| |
| In order to compute moment of inertia of a mass around an axis, the perpendicular vector from the mass to the axis is needed. If the axis L is defined by the unit vector '''S''' through the reference point '''R''', then the perpendicular vector from the line L to the point '''r''' is given by
| |
| :<math> \Delta\mathbf{r}_i^\perp = (\mathbf{r}_i-\mathbf{R}) - (\mathbf{S}\cdot(\mathbf{r}_i-\mathbf{R}))\mathbf{S} = [[I]-[\mathbf{S}\mathbf{S}^T]](\Delta\mathbf{r}_i),</math>
| |
| where [I] is the identity matrix and ['''S''' '''S'''{{sup|T}}] is the outer product matrix formed from the unit vector '''S''' along the line L. Recall that skew-symmetric matrix [S] is constructed so that [S]'''y'''='''S''' x '''y'''. The matrix [I-SS<sup>T</sup>] in this equation subtracts the component of Δ'''r'''='''r'''-'''R''' that is parallel to '''S'''.
| |
| | |
| The previous sections show that in computing the moment of inertia matrix this operator yields a similar operator using the components of the vector Δ'''r''' that is
| |
| :<math> [I|\Delta\mathbf{r}|^2-\Delta\mathbf{r}\Delta\mathbf{r}^T].</math>
| |
| | |
| It is helpful to keep the following identities in mind In order to compare the equations that define the inertia tensor and the inertia matrix.
| |
| | |
| Let [R] be the skew symmetric matrix associated with the position vector '''R'''=(x, y, z), then the product in the inertia matrix becomes
| |
| :<math> -[R]^2= -\begin{bmatrix} 0 & -z & y \\ z & 0 & -x \\ -y & x & 0 \end{bmatrix}^2 = \begin{bmatrix}
| |
| y^2+z^2 & -xy & -xz \\ -y x & x^2+z^2 & -yz \\ -zx & -zy & x^2+y^2 \end{bmatrix}.</math>
| |
| | |
| This can be viewed as another way of computing the perpendicular distance from an axis to a point, because the matrix formed by the outer product ['''R''' '''R'''{{sup|T}}] yields the identify
| |
| :<math> -[R]^2 = |\mathbf{R}|^2[I] -[\mathbf{R}\mathbf{R}^T]=
| |
| \begin{bmatrix} x^2+y^2+z^2 & 0 & 0 \\ 0& x^2+y^2+z^2 & 0 \\0& 0& x^2+y^2+z^2 \end{bmatrix}- \begin{bmatrix}x^2 & xy & xz \\ yx & y^2 & yz \\ zx & zy & z^2\end{bmatrix},</math>
| |
| where [I] is the 3x3 identity matrix.
| |
| | |
| Also notice, that
| |
| :<math> |\mathbf{R}|^2 = \mathbf{R}\cdot\mathbf{R} =\operatorname{tr}[\mathbf{R}\mathbf{R}^T],</math>
| |
| where ''tr'' denotes the sum of the diagonal elements of the outer product matrix, known as its trace.
| |
| | |
| == The inertia matrix in different reference frames ==
| |
| The use of the inertia matrix in Newton's second law assumes its components are computed relative to axes parallel to the inertial frame and not relative to a body-fixed reference frame.<ref name="Kane"/><ref name="Goldstein"/> This means that as the body moves the components of the inertia matrix change with time. In contrast, the components of the inertia matrix measured in a body-fixed frame are constant.
| |
| | |
| ===Body frame inertia matrix===
| |
| Let the body frame inertia matrix relative to the center of mass be denoted [I{{sub|C}}{{sup|B}}], and define the orientation of the body frame relative to the inertial frame by the rotation matrix [A], such that,
| |
| :<math>\mathbf{x}=[A]\mathbf{y},</math>
| |
| where vectors '''y''' in the body fixed coordinate frame have coordinates '''x''' in the inertial frame. Then, the inertia matrix of the body measured in the inertial frame is given by
| |
| :<math> [I_C]=[A][I_C^B][A^T].</math>
| |
| Notice that [A] changes as the body moves, while [I{{sub|C}}{{sup|B}}] remains constant.
| |
| | |
| ===Principal axes===
| |
| Measured in the body frame the inertia matrix is a constant real symmetric matrix. A real symmetric matrix has the [[eigendecomposition of a matrix|eigendecomposition]] into the product of a rotation matrix [Q] and a diagonal matrix [Λ], given by
| |
| :<math>[I_C^B]=[Q][\Lambda][Q^T],</math>
| |
| where
| |
| :<math>[\Lambda]= \begin{bmatrix}
| |
| I_{1} & 0 & 0 \\
| |
| 0 & I_{2} & 0 \\
| |
| 0 & 0 & I_{3}
| |
| \end{bmatrix}.
| |
| </math>
| |
| The columns of the rotation matrix [Q] define the directions of the principal axes of the body, and the constants I{{sub|1}}, I{{sub|2}} and I{{sub|3}} are called the '''principal moments of inertia'''. This result was first shown by [[James Joseph Sylvester|J. J. Sylvester (1852)]], and is a form of [[Sylvester's law of inertia]].<ref name=syl852>
| |
| {{cite journal|author=Sylvester, J J | title=A demonstration of the theorem that every homogeneous quadratic polynomial is reducible by real orthogonal substitutions to the form of a sum of positive and negative squares | journal=Philosophical Magazine (Ser. 4)| volume=4 | issue=23 | pages=138–142 | year=1852 | url=http://www.maths.ed.ac.uk/~aar/sylv/inertia.pdf | doi= 10.1080/14786445208647087 | accessdate=2008-06-27}}
| |
| </ref><ref name=norm>
| |
| {{cite book|author=Norman, C.W.| title=Undergraduate algebra | publisher=[[Oxford University Press]] | pages=360–361 | year=1986 | isbn=0-19-853248-2 }}
| |
| </ref>
| |
| | |
| For bodies with constant density an axis of rotational symmetry is a principal axis.
| |
| | |
| ===Inertia ellipsoid===
| |
| [[File:Triaxial Ellipsoid.jpg|thumb|right|An ellipsoid with the semi-principal diameters labeled a, b, and c.]]
| |
| The moment of inertia matrix in body-frame coordinates is a quadratic form that defines a surface in the body called [[Poinsot's ellipsoid]].<ref>
| |
| {{cite book
| |
| |first=Matthew T.
| |
| |last=Mason
| |
| |url=http://books.google.com/books?id=Ngdeu3go014C
| |
| |title=Mechanics of Robotics Manipulation
| |
| |publisher=MIT Press
| |
| |year=2001
| |
| |isbn=978-0-262-13396-8
| |
| |accessdate=2012-11-16
| |
| }}</ref> Let [Λ] be the inertia matrix relative to the center of mass aligned with the principal axes, then the surface
| |
| :<math> \mathbf{x}^T[\Lambda]\mathbf{x}=1,</math>
| |
| or
| |
| :<math> I_1x^2 + I_2y^2 + I_3z^2 =1,</math>
| |
| defines an [[ellipsoid]] in the body frame. Write this equation in the form,
| |
| :<math> \frac{x^2}{(1/\sqrt{I_1})^2} + \frac{y^2}{(1/\sqrt{I_2})^2} + \frac{z^2}{(1/\sqrt{I_3})^2} = 1,</math>
| |
| to see that the semi-principal diameters of this ellipsoid are given by
| |
| :<math> a=\frac{1}{\sqrt{I_1}}, \quad b=\frac{1}{\sqrt{I_2}}, \quad c=\frac{1}{\sqrt{I_3}}.</math>
| |
| | |
| Let a point '''x''' on this ellipsoid be defined in terms of its magnitude and direction, '''x'''=|'''x'''|'''n''', where '''n''' is a unit vector. Then the relationship presented above, between the inertia matrix and the scalar moment of inertia I{{sub|n}} around an axis in the direction '''n''', yields
| |
| :<math> \mathbf{x}^T[\Lambda]\mathbf{x}=|\mathbf{x}|^2\mathbf{n}^T[\Lambda]\mathbf{n} = |\mathbf{x}|^2I_n = 1. </math>
| |
| Thus, the magnitude of a point '''x''' in the direction '''n''' on the inertia ellipsoid is
| |
| :<math> |\mathbf{x}| = \frac{1}{\sqrt{I_n}}.</math>
| |
| | |
| <!---duplicated above
| |
| ===Parallel axis theorem for the inertia matrix===
| |
| {{Main|Parallel axis theorem}}
| |
| It is useful to note here that if the moment of inertia matrix or tensor is relative to the [[center of mass]], then it can be determined relative to any other reference point in the body using the parallel axis theorem. If [I{{sub|C}}{{sup|B}}] is the moment of inertia matrix in the body frame relative to the center of mass '''C''', then the moment of inertia matrix [I{{sub|R}}{{sup|B}}] in the same frame but relative to a different point '''R''' is given by
| |
| :<math>[I_R^B] = [I_C^B] - M[d]^2, </math>
| |
| where M is the mass of the body, and [d] is the skew-symmetric matrix obtained from the vector '''d'''= '''R'''-'''C'''.
| |
| | |
| The tensor form of the parallel axis theorem is given by
| |
| : <math> \mathbf{I}_R^B = \mathbf{I}_C^B + M((\mathbf{d} \cdot \mathbf{d}) \mathbf{E} - \mathbf{d} \otimes \mathbf{d}).
| |
| </math>
| |
| -->
| |
| | |
| == See also ==
| |
| * [[Central moment]]
| |
| * [[Instant center of rotation]]
| |
| * [[List of moments of inertia]]
| |
| * [[List of moment of inertia tensors]]
| |
| * [[Rotational energy]]
| |
| * [[Stretch rule]]
| |
| * [[Tire balance]]
| |
| | |
| == References ==
| |
| {{Reflist}}
| |
| <!--
| |
| * {{cite book|last=Goldstein|first=H.|year=1980|title=Classical Mechanics|edition=2nd|publisher=Addison-Wesley|isbn=0-201-02918-9}}.
| |
| * {{cite book|last1=Landau|first1=LD|last2=Lifshitz|first=EM|year=1976|title=Mechanics|edition=3rd|publisher=Pergamon Press|isbn=0-08-021022-8 (hardcover)}}; ISBN 0-08-029141-4 (softcover).}}
| |
| * {{cite book|last1=Marion|first1=JB|last2=Thornton|first2=ST. |year=1995|title=Classical Dynamics of Systems and Particles|edition=4th|publisher=Thomson|isbn=0-03-097302-3}}.
| |
| * {{citation|last=Sylvester|first=J J |authorlink=James Joseph Sylvester| title=A demonstration of the theorem that every homogeneous quadratic polynomial is reducible by real orthogonal substitutions to the form of a sum of positive and negative squares | journal=Philosophical Magazine | volume=IV | pages=138–142 | year=1852 | url=http://www.maths.ed.ac.uk/~aar/sylv/inertia.pdf | accessdate=2008-06-27}}
| |
| * {{cite book|last=Symon|first=KR|year=1971|title=Mechanics|edition=3rd|publisher=Addison-Wesley|isbn=0-201-07392-7}}.
| |
| * {{cite book|last=Tenenbaum|first=RA|year=2004|title=Fundamentals of Applied Dynamics|publisher=Springer|isbn=0-387-00887-X}}.
| |
| -->
| |
| | |
| == External links ==
| |
| {{Commons category|Moments of inertia}}
| |
| * [http://www.lightandmatter.com/html_books/0sn/ch04/ch04.html Angular momentum and rigid-body rotation in two and three dimensions]
| |
| * [http://hyperphysics.phy-astr.gsu.edu/hbase/mi.html Lecture notes on rigid-body rotation and moments of inertia]
| |
| * [http://kwon3d.com/theory/moi/iten.html The moment of inertia tensor]
| |
| * [http://www.phy.hk/wiki/englishhtm/Balance.htm An introductory lesson on moment of inertia: keeping a vertical pole not falling down (Java simulation)]
| |
| * [http://hypertextbook.com/physics/mechanics/rotational-inertia/ Tutorial on finding moments of inertia, with problems and solutions on various basic shapes]
| |
| * [http://www.cs.cmu.edu/afs/cs/academic/class/16741-s07/www/ Notes on mechanics of manipulation: the angular inertia tensor]
| |
| | |
| {{tensors}}
| |
| | |
| [[Category:Concepts in physics]]
| |
| [[Category:Mechanics]]
| |
| [[Category:Physical quantities]]
| |
| [[Category:Rigid bodies]]
| |
| [[Category:Quadratic forms]]
| |
| [[Category:Rotation]]
| |
| [[Category:Rotation in three dimensions]]
| |
| | |
| {{Link GA|de}}
| |