next up previous
Next: Eulerian Angles Up: Rigid Body Rotation Previous: Principal Axes of Rotation

Euler's Equations

The fundamental equation of motion of a rotating body [see Equation (456)],
\begin{displaymath}
{\bf T} = \frac{d{\bf L}}{dt},
\end{displaymath} (501)

is only valid in an inertial frame. However, we have seen that ${\bf L}$ is most simply expressed in a frame of reference whose axes are aligned along the principal axes of rotation of the body. Such a frame of reference rotates with the body, and is, therefore, non-inertial. Thus, it is helpful to define two Cartesian coordinate systems, with the same origins. The first, with coordinates $x$, $y$, $z$, is a fixed inertial frame--let us denote this the fixed frame. The second, with coordinates $x'$, $y'$, $z'$, co-rotates with the body in such a manner that the $x'$-, $y'$-, and $z'$-axes are always pointing along its principal axes of rotation--we shall refer to this as the body frame. Since the body frame co-rotates with the body, its instantaneous angular velocity is the same as that of the body. Hence, it follows from the analysis in Section 7.2 that
\begin{displaymath}
\frac{d{\bf L}}{dt} = \frac{d{\bf L}}{dt'} + \mbox{\boldmath$\omega$}\times{\bf L}.
\end{displaymath} (502)

Here, $d/dt$ is the time derivative in the fixed frame, and $d/dt'$ the time derivative in the body frame. Combining Equations (501) and (502), we obtain
\begin{displaymath}
{\bf T} = \frac{d{\bf L}}{dt'} + \mbox{\boldmath$\omega$}\times{\bf L}.
\end{displaymath} (503)

Now, in the body frame let ${\bf T}= (T_{x'},\,T_{y'},T_{z'})$ and $\mbox{\boldmath$\omega$}= (\omega_{x'},\,\omega_{y'},\,\omega_{z'})$. It follows that ${\bf L} = (I_{x'x'}\,\omega_{x'},\,I_{y'y'}\,\omega_{y'},\,I_{z'z'}\,\omega_{z'})$, where $I_{x'x'}$, $I_{y'y'}$ and $I_{z'z'}$ are the principal moments of inertia. Hence, in the body frame, the components of Equation (503) yield
$\displaystyle T_{x'}$ $\textstyle =$ $\displaystyle I_{x'x'}\,\dot{\omega}_{x'} - (I_{y'y'}-I_{z'z'})\,\omega_{y'}\,\omega_{z'},$ (504)
$\displaystyle T_{y'}$ $\textstyle =$ $\displaystyle I_{y'y'}\,\dot{\omega}_{y'} - (I_{z'z'}-I_{x'x'})\,\omega_{z'}\,\omega_{x'},$ (505)
$\displaystyle T_{z'}$ $\textstyle =$ $\displaystyle I_{z'z'}\,\dot{\omega}_{z'} - (I_{x'x'}-I_{y'y'})\,\omega_{x'}\,\omega_{y'},$ (506)

where $\dot{~}=d/dt'$. Here, we have made use of the fact that the moments of inertia of a rigid body are constant in time in the co-rotating body frame. The above equations are known as Euler's equations.

Consider a rigid body which is constrained to rotate about a fixed axis with constant angular velocity. It follows that $\dot{\omega}_{x'}=\dot{\omega}_{y'} = \dot{\omega}_{z'}=0$. Hence, Euler's equations, (504)-(506), reduce to

$\displaystyle T_{x'}$ $\textstyle =$ $\displaystyle - (I_{y'y'}-I_{z'z'})\,\omega_{y'}\,\omega_{z'},$ (507)
$\displaystyle T_{y'}$ $\textstyle =$ $\displaystyle - (I_{z'z'}-I_{x'x'})\,\omega_{z'}\,\omega_{x'},$ (508)
$\displaystyle T_{z'}$ $\textstyle =$ $\displaystyle - (I_{x'x'}-I_{y'y'})\,\omega_{x'}\,\omega_{y'}.$ (509)

These equations specify the components of the steady (in the body frame) torque exerted on the body by the constraining supports. The steady (in the body frame) angular momentum is written
\begin{displaymath}
{\bf L} = (I_{x'x'}\,\omega_{x'},\,I_{y'y'}\,\omega_{y'},\,
I_{z'z'}\,\omega_{z'}).
\end{displaymath} (510)

It is easily demonstrated that ${\bf T}= \mbox{\boldmath$\omega$}\times{\bf L}$. Hence, the torque is perpendicular to both the angular velocity and the angular momentum vectors. Note that if the axis of rotation is a principal axis then two of the three components of $\omega$ are zero (in the body frame). It follows from Equations (507)-(509) that all three components of the torque are zero. In other words, zero external torque is required to make the body rotate steadily about a principal axis.

Suppose that the body is freely rotating: i.e., there are no external torques. Furthermore, let the body be rotationally symmetric about the $z'$-axis. It follows that $I_{x'x'} = I_{y'y'} = I_\perp$. Likewise, we can write $I_{z'z'} = I_\parallel$. In general, however, $I_\perp\neq I_\parallel$. Thus, Euler's equations yield

$\displaystyle I_\perp\,\frac{d\omega_{x'}}{dt'} + (I_{\parallel}-I_{\perp})\,\omega_{z'}\,\omega_{y'}$ $\textstyle =$ $\displaystyle 0,$ (511)
$\displaystyle I_\perp\,\frac{d\omega_{y'}}{dt'} - (I_{\parallel}-I_{\perp})\,
\omega_{z'}\,\omega_{x'}$ $\textstyle =$ $\displaystyle 0,$ (512)
$\displaystyle \frac{d\omega_{z'}}{dt'}$ $\textstyle =$ $\displaystyle 0.$ (513)

Clearly, $\omega_{z'}$ is a constant of the motion. Equation (511) and (512) can be written
$\displaystyle \frac{d\omega_{x'}}{dt'} + {\mit\Omega}\,\omega_{y'}$ $\textstyle =$ $\displaystyle 0,$ (514)
$\displaystyle \frac{d\omega_{y'}}{dt'} - {\mit\Omega}\,\omega_{x'}$ $\textstyle =$ $\displaystyle 0,$ (515)

where ${\mit\Omega}= (I_{\parallel}/I_\perp-1)\,\omega_{z'}$. As is easily demonstrated, the solution to the above equations is
$\displaystyle \omega_{x'}$ $\textstyle =$ $\displaystyle \omega_\perp\,\cos({\mit\Omega}\,t'),$ (516)
$\displaystyle \omega_{y'}$ $\textstyle =$ $\displaystyle \omega_\perp\,\sin({\mit\Omega}\,t'),$ (517)

where $\omega_\perp$ is a constant. Thus, the projection of the angular velocity vector onto the $x'$-$y'$ plane has the fixed length $\omega_\perp$, and rotates steadily about the $z'$-axis with angular velocity ${\mit\Omega}$. It follows that the length of the angular velocity vector, $\omega=(\omega_{x'}^2+\omega_{y'}^2+\omega_{z'}^2)^{1/2}$, is a constant of the motion. Clearly, the angular velocity vector makes some constant angle, $\alpha$, with the $z'$-axis, which implies that $\omega_{z'} = \omega\,\cos\alpha$ and $\omega_\perp = \omega\,\sin\alpha$. Hence, the components of the angular velocity vector are
$\displaystyle \omega_{x'}$ $\textstyle =$ $\displaystyle \omega\,\sin\alpha\,\cos({\mit\Omega}\,t'),$ (518)
$\displaystyle \omega_{y'}$ $\textstyle =$ $\displaystyle \omega\,\sin\alpha\,\sin({\mit\Omega}\,t'),$ (519)
$\displaystyle \omega_{z'}$ $\textstyle =$ $\displaystyle \omega\,\cos\alpha,$ (520)

where
\begin{displaymath}
{\mit\Omega} =\omega\,\cos\alpha \left(\frac{I_\parallel}{I_\perp}-1\right).
\end{displaymath} (521)

We conclude that, in the body frame, the angular velocity vector precesses about the symmetry axis (i.e., the $z'$-axis) with the angular frequency ${\mit\Omega}$. Now, the components of the angular momentum vector are
$\displaystyle L_{x'}$ $\textstyle =$ $\displaystyle I_\perp\,\omega\,\sin\alpha\,\cos({\mit\Omega}\,t'),$ (522)
$\displaystyle L_{y'}$ $\textstyle =$ $\displaystyle I_\perp\,\omega\,\sin\alpha\,\sin({\mit\Omega}\,t'),$ (523)
$\displaystyle L_{z'}$ $\textstyle =$ $\displaystyle I_\parallel\,\omega\,\cos\alpha.$ (524)

Thus, in the body frame, the angular momentum vector is also of constant length, and precesses about the symmetry axis with the angular frequency ${\mit\Omega}$. Furthermore, the angular momentum vector makes a constant angle $\theta $ with the symmetry axis, where
\begin{displaymath}
\tan\theta = \frac{I_\perp}{I_\parallel}\,\tan\alpha.
\end{displaymath} (525)

Note that the angular momentum vector, the angular velocity vector, and the symmetry axis all lie in the same plane: i.e., ${\bf e}_{z'}\cdot{\bf L}\times \mbox{\boldmath$\omega$}=0$, as can easily be verified. Moreover, the angular momentum vector lies between the angular velocity vector and the symmetry axis (i.e., $\theta<\alpha$) for a flattened (or oblate) body (i.e., $I_\perp< I_\parallel$), whereas the angular velocity vector lies between the angular momentum vector and the symmetry axis (i.e., $\theta>\alpha$) for an elongated (or prolate) body (i.e., $I_\perp>I_\parallel$).

Let us now consider the most general motion of a freely rotating asymmetric rigid body, as seen in the body frame. Since a freely rotating body experiences no external torques, its angular momentum vector ${\bf L}$ is a constant of the motion in the inertial fixed frame. In general, the direction of this vector varies with time in the non-inertial body frame, but its length remains fixed. This can be seen from Equation (502): if $d{\bf L}/dt=0$ then the scalar product of this equation with ${\bf L}$ implies that $d L^2/dt'=0$. It follows from Equation (510) that

\begin{displaymath}
L^2 = I_{x'x'}^{\,2}\,\omega_{x'}^{\,2} + I_{y'y'}^{\,2}\,\o...
...'}^{\,2}
+ I_{z'z'}^{\,2}\,\omega_{z'}^{\,2} ={\rm constant}.
\end{displaymath} (526)

The above constraint can also be derived directly from Euler's equations, (504)-(506), by setting $T_{x'}=T_{y'}=T_{z'}=0$. A freely rotating body subject to no external torques clearly has a constant rotational kinetic energy. Hence, from Equation (469),
\begin{displaymath}
\mbox{\boldmath$\omega$}\cdot {\bf L} = I_{x'x'}\,\omega_{x'...
...ega_{y'}^{\,2}
+ I_{z'z'}\,\omega_{z'}^{\,2} ={\rm constant}.
\end{displaymath} (527)

This constraint can also be derived directly from Euler's equations. We conclude that, in the body frame, the components of $\mbox{\boldmath$\omega$}$ must simultaneously satisfy the two constraints (526) and (527). These constraints are the equations of two ellipsoids whose principal axes coincide with the principal axes of the body, and whose principal radii are in the ratio $I_{x'x'}^{-1}:I_{y'y'}^{-1}:I_{z'z'}^{-1}$ and $I_{x'x'}^{-1/2}:I_{y'y'}^{-1/2}:I_{z'z'}^{-1/2}$, respectively. In general, the intersection of these two ellipsoids is a closed curve. Hence, we conclude that the most general motion of a freely rotating asymmetric body, as seen in the body frame, is a form of irregular precession in which the tip of the angular velocity vector $\mbox{\boldmath$\omega$}$ periodically traces out the aforementioned closed curve. It is easily demonstrated that the tip of the angular momentum vector ${\bf L}$ periodically traces out a different closed curve.


next up previous
Next: Eulerian Angles Up: Rigid Body Rotation Previous: Principal Axes of Rotation
Richard Fitzpatrick 2011-03-31