13.7: Diagonalize the Inertia Tensor
( \newcommand{\kernel}{\mathrm{null}\,}\)
Finding the three principal axes involves diagonalizing the inertia tensor, which is the classic eigenvalue problem discussed in appendix 19.1. Solution of the eigenvalue problem for rigid-body motion corresponds to a rotation of the coordinate frame to the principal axes resulting in the matrix
{I}⋅ω=Iω
where I comprises the three-valued eigenvalues, while the corresponding vector ω is the eigenvector. Appendix 19.1 gives the solution of the matrix relation
{I}⋅ω=I{I}ω
where I are three-valued eigen values for the principal axis moments of inertia, and {I} is the unity tensor, equation (A.2.4).
{I}≡{100010001}
Rewriting ??? gives
({I}−I{I})⋅ω=0
This is a matrix equation of the form A⋅ω=0 where A is a 3×3 matrix and ω is a vector with values ωx,ωy,ωz. The matrix equation A⋅ω=0 really corresponds to three simultaneous equations for the three numbers ωx,ωy,ωz. It is a well-known property of equations like ??? that they have a non-zero solution if, and only if, the determinant det(A) is zero, that is
det(I−II)=0
This is called the characteristic equation, or secular equation for the matrix I. The determinant involved is a cubic equation in the value of I that gives the three principal moments of inertia. Inserting one of the three values of I into equation (13.4.6) gives the corresponding eigenvector ω. Applying the above eigenvalue problem to rigid-body rotation corresponds to requiring that some arbitrary set of body-fixed axes be the principal axes of inertia. This is obtained by rotating the body-fixed axis system such that
L1=I11ω1+I12ω2+I13ω3=Iω1L2=I21ω1+I22ω2+I23ω3=Iω2L3=I31ω1+I32ω2+I33ω3=Iω3
or
(I11−I)ω1+I12ω2+I13ω3=0I21ω1+(I22−I)ω2+I23ω3=0I31ω1+I32ω2+(I33−I)ω3=0
These equations have a non-trivial solution for the ratios ω1:ω2:ω3 since the determinant vanishes, that is
|(I11−I)I12I13I21(I22−I)I23I31I32(I33−I)|=0
The expansion of this determinant leads to a cubic equation with three roots for I. This is the secular equation for I whose eigenvalues are the principal moments of inertia.
The directions of the principal axes, that is the eigenvectors, can be found by substituting the corresponding solution for I into the prior equation. Thus for eigensolution I1 the eigenvector is given by solving
(I11−I1)ω11+I12ω21+I13ω31=0I21ω11+(I22−I1)ω21+I23ω31=0I31ω11+I32ω21+(I33−I1)ω31=0
These equations are solved for the ratios ω11:ω21:ω31 which are the direction numbers of the principle axis system corresponding to solution I1. This principal axis system is defined relative to the original coordinate system. This procedure is repeated to find the orientation of the other two mutually perpendicular principal axes.