Key Theorems for computing \(\left[I\right]^{S/O}\)¶
There are three useful theorems to compute moments of inertia. They are:
Rotation theorem
Parallel axis theorem (or translation theorem)
Composite theorem
Let us examine them in further detail.
Rotation theorem¶
\(X-Y-Z\) make up a cartesian coordinate system.
\(O\) is the origin.
\(\hat{n}_x,\;\hat{n}_y,\;\hat{n}_z\) are unit vectors directed along \(X-Y-Z\) respectively. The unit vectos represent a reference frame \(N\).
Also, the figure shows a rocket which has a reference frame \(R\) attached to it. You are given the inertia matrix of the rocket \(R\) about \(O\) along the unit vectors of frame \(N\). Naturally, this reference frame also has \(3\) mutually orthogonal unit vectors:
Now, we can define another inertia matrix for \(R\) about \(O\) along the newly introduced rotating reference frame’s unit vectors \(\hat{r}_x, \; \hat{r}_y, \; \hat{r}_z\). This matrix and its elements are:
Question: How are the two matrices related?
We begin by assuming that \(R\) is a particle of mass, \(m_R\).
Let us consider a product of inertia from each matrix.
\(P^*\) is a position vector from \(O\) to the origin of the frame \(R\).
We know from our discussion on direction \(\text{cosine}\) matrices that the unit vectors of \(R\) can be related to the unit vectors of \(N\) as:
\({}^{R}C^{N}\) is the direction \(\text{cosine}\) matrix of \(R\) in \(N\). It is a \(3\times3\) matrix, which we shall assume has the following elements:
So, we can now easily derive:
And then, we can rewrite the product of inertia \(J^{R/O}_{xy}\) as:
This leads to
or in terms of matrix multiplication.
More generally, it can be shown that
Parallel axes theorem¶
In the previoussection we presumed that the rocket was a point mass. What happens if this assumption is relaxed and rockets are better approximated as a system of particles?
\(P^*\) is a position vector from \(O\) to \(R^*\), the origin of the frame \(R\).
\(R^*\) is he mass center of the rocket.
\(\vec{P}_i\) is a position vector from \(O\) to \(\vec{P}_i\), the \(i^{\text{th}}\) particle on the rocket body.
\(r_i\) is a position vector from \(R^*\) to \(\vec{P}_i\). So,
Now, we know that the product of inertia for this system of particles along the \(\hat{n}_x,\;\hat{n}_y,\;\hat{n}_z\) directions given by:
The terms highlighted in green go to zero because \(\sum_im_ir_i = 0\) by the definition of mass centers.
\(I^{R/R^*}_{xy}\) is the product of inertia of \(R\) about the mass center \(R^*\).
\(I^{R^*/O}_{xy}\)is the product of inertia of \(R^*\) about \(O\).
This result extrapolates to the moment of inertia scalars. For example,
This result extrapolates to the inertia matrix,
This rule is only valid when going through the mass center.