2.1.22Analytical Mechanics

Inertia tensor — principal axes, principal moments

1,943 words9 min readdifficulty · medium4 backlinks

1. WHAT is the inertia tensor — derived from scratch

We want L\vec{L} for a rigid body rotating about a fixed point with angular velocity ω\vec{\omega}.

Start from one particle. Particle aa at position ra\vec{r}_a, mass mam_a. In a rigid body its velocity is va=ω×ra\vec{v}_a = \vec{\omega}\times\vec{r}_a. Its angular momentum:

La=mara×va=mara×(ω×ra)\vec{L}_a = m_a\,\vec{r}_a\times\vec{v}_a = m_a\,\vec{r}_a\times(\vec{\omega}\times\vec{r}_a)

Why this step? Angular momentum is r×mv\vec{r}\times m\vec{v} by definition, and rigidity forces v=ω×r\vec{v}=\vec\omega\times\vec r.

Use the BAC–CAB identity A×(B×C)=B(AC)C(AB)\vec{A}\times(\vec{B}\times\vec{C}) = \vec{B}(\vec{A}\cdot\vec{C}) - \vec{C}(\vec{A}\cdot\vec{B}):

La=ma[ω(ra2)ra(raω)]\vec{L}_a = m_a\big[\vec{\omega}\,(r_a^2) - \vec{r}_a(\vec{r}_a\cdot\vec{\omega})\big]

Why this step? This separates ω\vec\omega out so we can read off a linear map acting on it.

Sum over particles (or integrate dm\int dm). Look at one component, say LxL_x. With r=(x,y,z)\vec r=(x,y,z) and ω=(ωx,ωy,ωz)\vec\omega=(\omega_x,\omega_y,\omega_z):

Lx=ama[ωx(x2+y2+z2)xa(xaωx+yaωy+zaωz)]L_x = \sum_a m_a\big[\omega_x(x^2+y^2+z^2) - x_a(x_a\omega_x + y_a\omega_y + z_a\omega_z)\big] =ωx ⁣ma(y2+z2)Ixxωy ⁣maxyIxyωz ⁣maxzIxz= \omega_x\!\underbrace{\sum m_a(y^2+z^2)}_{I_{xx}} - \omega_y\!\underbrace{\sum m_a\,x y}_{-I_{xy}} - \omega_z\!\underbrace{\sum m_a\,x z}_{-I_{xz}}

So L=Iω\vec{L} = \mathbf{I}\,\vec{\omega} where:


2. WHY L\vec L isn't parallel to ω\vec\omega — and the fix

Because I\mathbf I has off-diagonal terms, Iω\mathbf I\vec\omega generally points elsewhere. A symmetric real matrix can always be diagonalized by an orthogonal (rotation) matrix. That means there exists a special rotated coordinate frame where:

I=(I1000I2000I3)\mathbf I' = \begin{pmatrix} I_1&0&0\\0&I_2&0\\0&0&I_3\end{pmatrix}

WHY these are exactly the eigenvectors: "L\vec L parallel to ω\vec\omega" means Iω=λω\mathbf I\vec\omega = \lambda\vec\omega. That is the eigenvalue equation. So:

HOW to find them (recipe):

  1. Write I\mathbf I in any convenient frame.
  2. Solve the characteristic cubic det(Iλ1)=0\det(\mathbf I - \lambda\mathbb 1)=0 → three roots I1,I2,I3I_1,I_2,I_3.
  3. For each λ=Ik\lambda=I_k, solve (IIk1)e^k=0(\mathbf I-I_k\mathbb1)\hat e_k=0 → axis direction; normalize.
  4. Shortcut: any symmetry axis of the body is automatically a principal axis.
Figure — Inertia tensor — principal axes, principal moments

3. Rotational KE & a key invariant

T=12mava2=12ωL=12ωTIωT = \tfrac12\sum m_a v_a^2 = \tfrac12\,\vec\omega\cdot\vec L = \tfrac12\,\vec\omega^{\mathsf T}\mathbf I\,\vec\omega

In principal axes this becomes the clean diagonal form:

T=12(I1ω12+I2ω22+I3ω32)\boxed{\,T=\tfrac12(I_1\omega_1^2+I_2\omega_2^2+I_3\omega_3^2)\,}


4. Worked examples


5. Common mistakes (Steel-manned)


6. Flashcards

What does the inertia tensor relate?
It maps angular velocity to angular momentum: L=Iω\vec L=\mathbf I\vec\omega.
Why is L\vec L generally not parallel to ω\vec\omega?
Because I\mathbf I has nonzero off-diagonal (product of inertia) terms, so it rotates/scales ω\vec\omega unequally.
Define principal axes.
Directions in which I\mathbf I is diagonal; equivalently eigenvectors of I\mathbf I where Lω\vec L\parallel\vec\omega.
Define principal moments of inertia.
The eigenvalues I1,I2,I3I_1,I_2,I_3 of I\mathbf I — diagonal entries in the principal frame.
Formula for a diagonal moment, e.g. IxxI_{xx}.
Ixx=ma(ya2+za2)I_{xx}=\sum m_a(y_a^2+z_a^2) (distance² from the x-axis).
Formula for a product of inertia IxyI_{xy}.
Ixy=maxayaI_{xy}=-\sum m_a x_a y_a (note the minus sign).
Why is I\mathbf I symmetric?
Ixy=mxy=IyxI_{xy}=-\sum mxy=I_{yx}; products are symmetric in the two indices.
Equation that finds principal axes.
Ie^k=Ike^k\mathbf I\hat e_k=I_k\hat e_k with det(IIk1)=0\det(\mathbf I-I_k\mathbb1)=0.
Rotational KE in principal axes.
T=12(I1ω12+I2ω22+I3ω32)T=\tfrac12(I_1\omega_1^2+I_2\omega_2^2+I_3\omega_3^2).
Frame-independent check on principal moments.
I1+I2+I3=TrI=2mara2I_1+I_2+I_3=\operatorname{Tr}\mathbf I=2\sum m_a r_a^2.
Shortcut for finding a principal axis.
Any axis of symmetry of the body is automatically a principal axis.
What is a spherical top?
A body with I1=I2=I3I_1=I_2=I_3 (e.g. cube/sphere); every axis through CM is principal.

Recall Feynman: explain it to a 12-year-old

Imagine spinning a wonky lump of clay. If the lump is heavier on one side, when you twist it, it doesn't just spin straight — it wobbles and tries to flop sideways. The "wobble recipe" of the lump is a little grid of numbers (the inertia tensor). For most ways of spinning, the lump pushes its "spin-effort" sideways. But there are exactly three special straight lines through the lump where, if you spin around them, it spins clean and steady with no flop. Those three lines are the principal axes, and how hard it is to spin around each is the principal moment. Finding them is just asking: "which spin directions keep behaving themselves?"

Connections

Concept Map

has

rigidity gives

plug into

BAC-CAB identity

defines linear map

gives

diagonal entries

off-diagonal

nonzero means

symmetric so

yields

make I diagonal

restore

Rigid body spinning

Angular velocity omega

v equals omega cross r

L equals sum m r cross v

Separate omega out

Inertia tensor I

L equals I omega

Moments of inertia

Products of inertia

L not parallel to omega

Diagonalize by rotation

Principal axes

Principal moments

L parallel to omega

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, jab koi extended (phaila hua) body ghoomta hai, to har chhota tukda alag-alag tarike se rotation ko "resist" karta hai. Isi distribution ko capture karne ke liye hum inertia tensor I\mathbf{I} banate hain — ek 3×33\times3 matrix. Single particle me L\vec L aur ω\vec\omega parallel hote hain, par body me aksar nahi! Kyunki I\mathbf I ke andar off-diagonal (product of inertia) terms hote hain, jo ω\vec\omega ko tilt kar dete hain. Yahi reason hai ki wonky object ghumate waqt wobble karta hai.

Ab kamaal ki baat: I\mathbf I ek symmetric matrix hai, isliye hamesha kuch special directions milti hain jahan ghoomne par L\vec L wapas ω\vec\omega ke parallel ho jaata hai. In directions ko principal axes bolte hain, aur unke corresponding numbers ko principal moments I1,I2,I3I_1, I_2, I_3. Mathematically ye bilkul eigenvalue problem hai: Ie^k=Ike^k\mathbf I\hat e_k = I_k\hat e_k. Jis direction me Lω\vec L\parallel\vec\omega, wahi eigenvector, aur scaling factor wahi eigenvalue.

Practical shortcut: agar body me koi symmetry axis hai, to woh apne aap principal axis ban jaata hai — alag se diagonalize karne ki zaroorat nahi. Aur ek check rakho: $I_1+I_2+I_3 =

Go deeper — visual, from zero

Test yourself — Analytical Mechanics

Connections