Talk:PlanetPhysics/Inertia Tensor

Original TeX Content from PlanetPhysics Archive

%%% This file is part of PlanetPhysics snapshot of 2011-09-01 %%% Primary Title: inertia tensor %%% Primary Category Code: 45.40.-f %%% Filename: InertiaTensor.tex %%% Version: 7 %%% Owner: bloftin %%% Author(s): bloftin %%% PlanetPhysics is released under the GNU Free Documentation License. %%% You should have received a file called fdl.txt along with this file. %%% If not, please write to gnu@gnu.org. \documentclass[12pt]{article} \pagestyle{empty} \setlength{\paperwidth}{8.5in} \setlength{\paperheight}{11in}

\setlength{\topmargin}{0.00in} \setlength{\headsep}{0.00in} \setlength{\headheight}{0.00in} \setlength{\evensidemargin}{0.00in} \setlength{\oddsidemargin}{0.00in} \setlength{\textwidth}{6.5in} \setlength{\textheight}{9.00in} \setlength{\voffset}{0.00in} \setlength{\hoffset}{0.00in} \setlength{\marginparwidth}{0.00in} \setlength{\marginparsep}{0.00in} \setlength{\parindent}{0.00in} \setlength{\parskip}{0.15in}

\usepackage{html}

% this is the default PlanetMath preamble. as your knowledge % of TeX increases, you will probably want to edit this, but % it should be fine as is for beginners.

% almost certainly you want these \usepackage{amssymb} \usepackage{amsmath} \usepackage{amsfonts}

% used for TeXing text within eps files %\usepackage{psfrag} % need this for including graphics (\includegraphics) %\usepackage{graphicx} % for neatly defining theorems and propositions %\usepackage{amsthm} % making logically defined graphics %\usepackage{xypic}

% there are many more packages, add them here as you need them

% define commands here

\begin{document}

The inertia tensor is straight forward to calculate in theory from equation (1). However, in practice finding the inertia tensor of an \htmladdnormallink{object}{http://planetphysics.us/encyclopedia/TrivialGroupoid.html} is quite involved and high precision is needed in applications such as spacecraft design.

\begin{equation} I = \left[ \begin{array}{ccc} \displaystyle\sum_i m_i (r^2_i - x^2_i) & -\displaystyle\sum_i m_i x_i y_i & -\displaystyle\sum_i m_i x_i z_i \\ -\displaystyle\sum_i m_i y_i x_i & \displaystyle\sum_i m_i (r^2_i - y^2_i) & -\displaystyle\sum_i m_i y_i z_i \\ -\displaystyle\sum_i m_i z_i x_i & -\displaystyle\sum_i m_i z_i y_i & \displaystyle\sum_i m_i (r^2_i - z^2_i) \end{array} \right] \end{equation}

Units for the inertia tensor are $[kg/m^2]$. Terminology to reference the elements of the inertia tensor are the \htmladdnormallink{moments of inertia}{http://planetphysics.us/encyclopedia/MomentOfInertia.html}, which are the diagonal elements, $I_{11}$, $I_{22}$, and $I_{33}$.

and the products of inertia, which are the off-diagonal elements. A quick inspection of (1) reveals that the inertia tensor is symmetric, so the the products of inertia are

$$ I_{12} = I_{21} $$ $$ I_{13} = I_{31} $$ $$ I_{23} = I_{32} $$

To help build intuition on working with the inertia tensor, a few simple cases are looked at. For the example of a \htmladdnormallink{solid}{http://planetphysics.us/encyclopedia/CoIntersections.html} sphere, we found that the \htmladdnormallink{moment of inertia of a Solid Sphere}{http://planetphysics.us/encyclopedia/MomentOfInertiaOfASolidSphere.html} was

$$I = \frac{2}{5} M R^2 $$

This value was just the rotational inertia about a single axis that was a diameter of the sphere. It is easy to see that the general case about any of the sphere's body axes, x,y or z is taken care of by the inertia tensor with the moments of inertia being

$$ I_{11} = I_{22} = I_{33} = I $$

and the products of inertia equal to 0. For this case of the body axes being diameters of the sphere, then the body axes are also the principle axes and this says a lot about the symmetry of the chosen axes. For an axis of rotation that is not a diameter of the sphere, then there will be non-zero products of inertia and no more symmetry, so this axis is not a principle axis. If you visualize moving the origin of the body axes anywhere other than the orgin, then the set of coordinate body axes are not principle axes of rotation.

Another simple case is for the \htmladdnormallink{moment of inertia of a solid cylinder}{http://planetphysics.us/encyclopedia/MomentOfInertiaOfASolidCylinder.html}. Here we find that the diagonal elements are different based on their symmetry. For the moment of inertia about the z axis, we found

$$ I_{33} = \frac{1}{2} M R^2 $$

and for the symmetric axes of x and y

$$ I_{11} = I_{22} = \frac{1}{4} M R^2 + \frac{1}{12} M L^2 $$

If you imagine a very dense pencil, where R is small and L is large, you might have an inertia tensor that looks like

$$ I = \left[ \begin{array}{ccc} 100 & 0 & 0 \\ 0 & 100 & 0\\ 0 & 0 & 10 \end{array} \right] $$

Notice how this gives us a feel for how the \htmladdnormallink{mass}{http://planetphysics.us/encyclopedia/Mass.html} of the body is distributed. For a stubby cylinder, like a disk the inertia tensor would be something like

$$ I = \left[ \begin{array}{ccc} 10 & 0 & 0 \\ 0 & 10 & 0\\ 0 & 0 & 100 \end{array} \right] $$

\end{document}