# 4.6: Force-free Motion of a Rigid Asymmetric Top

- Page ID
- 6951

By “asymmetric top” I mean a body whose three principal moments of inertia are unequal. While we often think of a “top” as a symmetric body spinning on a table, in this section the “top” will not necessarily be symmetric, and it will not be in contact with any table, nor indeed subjected to any external forces or torques.

A complete description of the motion of an asymmetric top is quite complicated, and therefore all that we shall attempt in this chapter is a qualitative description of certain aspects of the motion. That our description is going to be “qualitative” does not by any means imply that this section is not going to be replete with equations or that we can give our poor brains a rest.

The first point that we can make is that, *provided that no external torques act on the body*, its angular momentum \(\bf L \) is constant in magnitude and direction. A second point is that, *provided the body is rigid and has no internal degrees of freedom*, the rotational kinetic energy \(T\) is constant. I deal briefly with nonrigid bodies in Section 4.7. Although the angular velocity vector \( \boldsymbol\omega \) is by no means fixed in either magnitude and direction, and the body can tumble over and over, these two conditions impose some constraints of the magnitude and direction of \( \boldsymbol\omega \).

We are going to examine these two conditions to see what constraints are imposed on \( \boldsymbol\omega \). One of the things we shall find is that rotation of a body about a principal axis of greatest or of least moment of inertia is stable against small displacements, whereas rotation about the principal axis of intermediate moment of inertia is unstable.

Absence of an external torque means that the angular momentum is constant:

\[ L^{2} = L^{2}_{1} + L^{2}_{2} + L^{2}_{3} = constant, \tag{4.6.1}\label{eq:4.6.1} \]

so that, at all times,

\[ I^{2}_{1}\omega_{1}^{2}+I^{2}_{2}\omega_{2}^{2}+I^{2}_{3}\omega_{3}^{2} = L^{2} \tag{4.6.2}\label{eq:4.6.2} \]

Thus, for a given \( L\). the angular velocity components always satisfy

\[ \frac{\omega^{2}_{1}}{(L/I_{1})^{2}} + \frac{\omega^{2}_{2}}{(L/I_{2})^{2}} + \frac{\omega^{2}_{3}}{(L/I_{3})^{2}}= 1. \tag{4.6.3}\label{eq:4.6.3} \]

That is to say, the angular velocity vector is constrained such that the tip of the vector \( \boldsymbol\omega\)** **is always on the surface of an ellipsoid of semi axes \( \frac{L}{I_{1}}\), \( \frac{L}{I_{2}}\), \( \frac{L}{I_{3}}\),

In addition to the constancy of angular momentum, the kinetic energy is also constant:

\[ \frac{1}{2}I_{1}\omega_{1}^{2} + \frac{1}{2}I_{2}\omega_{2}^{2} + \frac{1}{2}I_{3}\omega_{3}^{2}= T \tag{4.6.4}\label{eq:4.6.4} \]

Thus the tip of the angular velocity vector must also be on the surface of the ellipsoid

\[ \frac{\omega^{2}_{1}}{(\sqrt{2T/I_{1}})^{2}}+\frac{\omega^{2}_{2}}{(\sqrt{2T/I_{2}})^{2}} +\frac{\omega^{2}_{3}}{(\sqrt{2T/I_{3}})^{2}} = 1. \tag{4.6.5}\label{eq:4.6.5} \]

This ellipsoid (which is similar in shape to the momental ellipsoid) has semi axes \( \sqrt{2T/I_{1}}, \sqrt{2T/I_{2}}, \sqrt{2T/I_{3}}. \)

Thus, how ever the body tumbles over and over, \( \boldsymbol\omega \) is constrained in magnitude and direction so that its tip is on the curve where these two ellipses intersect.

Suppose, that we have a rigid body with

- \( I_{1}\) = 0.2 kgm
^{2} - \( I_{2}\) = 0.3 kgm
^{2} - \( I_{3}\) = 0.5 kgm
^{2}

and that we set it in motion such that the angular momentum and kinetic energy are L = 4 J s and T = 20 J.

(The angular momentum and kinetic energy will be determined by the magnitude and direction of the initial velocity vector by which it is set in motion.)

The tip of \( \boldsymbol\omega \) is constrained to be on the curve of intersection of the two ellipsoids

\[ \frac{\omega_{1}^{2}}{20^{2}}+\frac{\omega_{2}^{2}}{13.3^{2}}+ \frac{\omega_{3}^{2}}{8^{2}} = 1 \tag{4.6.6}\label{eq:4.6.6} \]

and

\[ \frac{\omega_{1}^{2}}{14.14^{2}}+\frac{\omega_{2}^{2}}{11.55^{2}}+ \frac{\omega_{3}^{2}}{8.94^{2}} = 1 \tag{4.6.7}\label{eq:4.6.7} \]

It is not easy (or I do not find it so) to imagine what this curve of intersection looks like in three-dimensional space, but one of my students, Leif Petersen, prepared the drawing below, and I am grateful to him for permission to reproduce it here. You can see that the curve of intersection is not a plane curve.

In case it’s of any help, you might want to note that equations \( \ref{eq:4.6.6}\) and \( \ref{eq:4.6.7}\) can be written

\[ 4 \omega^{2}_{1} + 9 \omega^{2}_{2} + 25 \omega^{2}_{3} = 1600 \tag{4.6.8}\label{eq:4.6.8} \]

and

\[ 2 \omega^{2}_{1} + 3 \omega^{2}_{2} + 5 \omega^{2}_{3} = 400 \tag{4.6.9}\label{eq:4.6.9} \]

but I’m going to leave the equations in the form \( \ref{eq:4.6.6}\) and \( \ref{eq:4.6.7}\), and in figure IV.7, I’ll sketch one octant of the two ellipsoidal surfaces.

The continuous blue curve shows an octant of the ellipsoid \( L\) = constant, and the dashed black curve shows an octant of the ellipsoid \( T\) = constant. The angular momentum vector can end only on the curve (not drawn) where the two ellipsoids intersect. Two points on the curve are indicated in Figure IV.7. If, for example, \( \boldsymbol\omega \) is oriented so that \( \omega_{1} \) = 0, the other two components must be \( \omega_{2} \) = 8.16 and \( \omega_{3} \) = 6.32. If it is oriented so that \( \omega_{3} \) = 0, the other two components must be \( \omega_{3} \) = 7.30 and \( \omega_{1} \) = 8.16. If \( \omega_{3} \) = 0, there are no real solutions for \( \omega_{1} \) and \( \omega_{2} \). This means that, for the given values of \( L\) and \( T\), \( \omega_{3} \) cannot be zero.

Now I’m going to address myself to the stability of rotation when a symmetric top is initially set to spin about one of its principal axes, which I’ll take to be the \( z\)-axis. We’ll suppose that initially \( \omega_{1} \) = \( \omega_{2} \) = 0, and \( \omega_{3} \) = \( \Omega\). In that case the angular momentum and the kinetic energy are \( L = I_{3}\Omega \) In any subsequent motion, the tip of \( T = I_{3}\Omega^{2} \) is restricted to move along the curve of intersection of the ellipsoids given by equations \( \ref{eq:4.6.3}\) and \( \ref{eq:4.6.5}\). That is to say, along the curve of intersection of the ellipsoids

\[ \frac{\omega_{1}^{2}}{(\frac{I_{3}}{I_{1}}\Omega)^2} + \frac{\omega_{2}^{2}}{(\frac{I_{3}}{I_{2}}\Omega)^2} +

\frac{\omega_{3}^{2}}{\Omega^{2}} = 1 \tag{4.6.10}\label{eq:4.6.10} \]

and

\[ \frac{\omega_{1}^{2}}{(\sqrt{\frac{I_{3}}{I_{1}}}\Omega)^2} + \frac{\omega_{2}^{2}}{(\sqrt{\frac{I_{3}}{I_{2}}}\Omega)^2} +

\frac{\omega_{3}^{2}}{\Omega^{2}} = 1 \tag{4.6.11}\label{eq:4.6.11} \]

For a specific example, I’ll suppose that the moments of inertia are in the ratio \( 2:3:5\), and we’ll consider three cases in turn.

Case I. Rotation about the axis of *least* moment of inertia. That is, we’ll take \( I_{3}=2\), \( I_{1}=3\), \( I_{2}=5\). Since \( I_{3}\) is the smallest moment of inertia, each of the ratios \( \frac{I_{3}}{I_{1}}\) and \( \frac{I_{3}}{I_{2}}\)_{ }are less than 1, and \( \sqrt{\frac{I_{3}}{I_{1}}} > \frac{I_{3}}{I_{1}} \quad and \sqrt{\frac{I_{3}}{I_{2}}} > \frac{I_{3}}{I_{2}} \) The two ellipsoids are

\[ \frac{\omega^{2}_{1}}{(0.667 \Omega)^2} + \frac{\omega^{2}_{2}}{(0.400 \Omega)^2} + \frac{\omega^{2}_{3}}{\Omega^{2}} = 1 \tag{4.6.12}\label{eq:4.6.12} \]

and

\[ \frac{\omega^{2}_{1}}{(0.816\Omega)^2} + \frac{\omega^{2}_{2}}{(0.632 \Omega)^2} + \frac{\omega^{2}_{3}}{\Omega^{2}} = 1 \tag{4.6.13}\label{eq:4.6.13} \]

I’ll try and sketch these:

Initially, we suppose, the body was set in motion rotating about the \( z\)-axis with angular speed \( \Omega \), which determines the values of \( L\) and \( T\), which will remain constant. The tip of the vector \( \omega \) is constrained to remain on the surface of the ellipsoid \( L\) = 0 and on the ellipsoid \( T\) = 0, and hence on the intersection of these two surfaces. But these two surface touch only at one point, namely (\( \omega_{1}, \omega_{2}, \omega_{3} \)) = (0 , 0 , \( \Omega \)). Thus there the vector \( \omega \) remains, and the rotation is stable.

Case II. Rotation about the axis of *greatest* moment of inertia. That is, we’ll take \( I_{3}=5\), \( I_{1}=2\), \( I_{2}=3\). Since \( I_{3}\) is the greatest moment of inertia, each of the ratios \( \frac{I_{3}}{I_{1}}\) and \( \frac{I_{3}}{I_{2}}\) are greater than 1, and \( \sqrt{\frac{I_{3}}{I_{1}}} < \frac{I_{3}}{I_{1}}\) and \( \sqrt{\frac{I_{3}}{I_{2}}} < \frac{I_{3}}{I_{2}}\) The two ellipsoids are

\[ \frac{\omega_{1}^{2}}{(2.50\Omega)^{2}} + \frac{\omega_{2}^{2}}{(1.67\Omega)^{2}} + \frac{\omega_{3}^{2}}{\Omega^{2}}=1 \tag{4.6.14}\label{eq:4.6.14} \]

and

\[ \frac{\omega_{1}^{2}}{(1.58\Omega)^{2}} + \frac{\omega_{2}^{2}}{(1.29\Omega)^{2}} + \frac{\omega_{3}^{2}}{\Omega^{2}}=1 \tag{4.6.15}\label{eq:4.6.15} \]

I’ll try and sketch these:

Again, and for the same reason as for Case I, we see that this motion is stable.

Case III. Rotation about the *intermediate* axis. That is, we’ll take \( I_{3}=3\), \( I_{1}=5\), \( I_{2}=2\). This time \( \frac{I_{3}}{I_{1}}\) is less than 1 and \( \frac{I_{3}}{I_{2}}\) is less than 1, and \( \sqrt{\frac{I_{3}}{I_{1}}} < \frac{I_{3}}{I_{1}}\) and \( \sqrt{\frac{I_{3}}{I_{2}}} < \frac{I_{3}}{I_{2}} \) The two ellipsoids are

\[ \frac{\omega^{2}_{1}}{(0.60 \Omega)} + \frac{\omega^{2}_{2}}{(1.50 \Omega)^2} + \frac{\omega^{2}_{3}}{\Omega^{2}} = 1 \tag{4.6.16}\label{eq:4.6.16} \]

and

\[ \frac{\omega^{2}_{1}}{(0.77\Omega)^2} + \frac{\omega^{2}_{2}}{(1.22 \Omega)^2} + \frac{\omega^{2}_{3}}{\Omega^{2}} = 1 \tag{4.6.17}\label{eq:4.6.17} \]

I'll try and sketch these:

Unlike the situation for Cases I and II, in which the two ellipsoids touch at only a single point, the two ellipses for Case III intersect in the curve shown as a dotted line in figure IV.10. Thus \( \omega\) is not restricted to lying along the \( z\)-axis, but it can move anywhere along the dotted line. The motion, therefore, is not stable.

You should experiment by throwing a body in the air in such a manner as to let it spin around one of its principal axes. A rectangular block will do, though the effect is particularly noticeable with something like a table-tennis bat.

Here is another approach to reach the same result. We imagine an asymmetric top spinning about one of its principal axes with angular velocity \( \boldsymbol\omega \ = \omega \hat{\bf z}). It is then given a small perturbation, so that its angular velocity is now.

\[ \omega = \epsilon \hat{ \bf x} + \eta\hat{\bf y} + \omega_{z}\hat{ \bf z} \tag{4.6.18}\label{eq:4.6.18} \]

Here the “hatted” quantities are the unit orthogonal vectors; \( \epsilon \) and \( \eta \) are supposed small compared with \( \omega_{z} \). Euler’s equations are :

\[ I_{1} \dot{\epsilon} = \eta \omega_{Z} (I_{2} - I_{3}), \tag{4.6.19}\label{eq:4.6.19} \]

\[ I_{1} \dot{\eta} = \omega_{Z} \epsilon (I_{3} - I_{1}), \tag{4.6.20}\label{eq:4.6.20} \]

\[ I_{1} \dot{\omega_{Z}} = \epsilon \eta (I_{1} - I_{2}). \tag{4.6.21}\label{eq:4.6.21} \]

If \( \epsilon \eta << \dot{\omega_{Z}} , \) then \( \omega_{Z}\) is approximately constant. Elimination of \( \eta \) from the first two equations yields

\[ \ddot{\epsilon} = - [\frac{(I_{2} - I_{3})(I_{1}-I_{3})\omega_{z}^{2}}{I_{1}I_{2}}]\epsilon \tag{4.6.22}\label{eq:4.6.22} \]

Elimination of \( \epsilon \) instead results in a similar equation in \( \eta \).

If \( I_{3}\)_{ }is either the largest or the smallest of the three moments of inertia, the two parentheses in the denominator have the same sign, so the expression in the brackets is positive. Equation \( \ref{eq:4.6.22}\) is then the equation for simple harmonic motion, and the motion is stable. If, however, \( I_{3}\) is intermediate between the other two, the two parentheses have opposite sign, and the expression in brackets in negative. In that case \( \epsilon\) and \( \eta\) increase exponentially, and the motion is unstable.

Mr Neil Honkanen of the University of Victoria conducted an experiment to illustrate the stability of rotation about the three principal axes. The body in question was a small “brick” of mild steel (density 7.83 g/cm^{3}) of dimensions 3/8 inch × 3/4 inch × 1 1/2 inch, mass 54.1 g. In round figures, this corresponds to principal moments of inertia \( A_{0}\) = 2 ×10^{-6 }kg m^{2}, \( B_{0}\) = 7 × 10^{-6 }kg m^{2}, \( C_{0}\) = 8 × 10^{-6 }kg m^{2}. He suspended it from an electromagnet, which he set in rotation at about 25 revolutions per second, and then let it fall, while photographing it stroboscopically. He did three experiments rotation respectively about the three principle axes. You can see from the photographs below that the rotation is stable when the rotation is about the axes of greatest or least moment of inertia, but is unstable when the rotation is about the axis of intermediate moment of inertia.