The restrictions of rigid body motion require that distance between any two particles within the body remain constant. For example
will be constant during the motion. Let's take the time derivative of this quantity,
For this to vanish either must vanish or the dot product must vanish. The first corresponds to a translation of the object.
For the dot product to vanish, the velocity difference must be perpendicular to the position difference. Specifically,
At first it might appear that you can have a different for each pair of particles. It turns out that this is not the case. Let's calculate
On the other hand we could have calculated it as follows,
For these expressions to be equal for arbitrary values of , and , the three values of must be equal to each other so we will drop the subscripts and call this quantity the angular velocity of the rotation.
To summarize, a rigid body has two types of motion:
Translation of the entire body, and
Uniform rotation about a particular axis (the direction of and an angular speed given by the magnitude of .
We can use some point within or outside the body as the origin and denote its position at any time by . The velocity of a particular part of the body is
where is the location of the part of the body and is the velocity of the origin.
Let's move the origin to . The velocity of the particle does not change so we have
so
The vector is constant with respect to changes in the origin. A useful example of this is a wheel that rolls without slipping.
Rigid bodies are such an important part of classical mechanics that we have developed special techniques to calculate their kinetic energies and angular momenta. Let's calculate the kinetic energy of a mass in rigid body motion. Its velocity is
where we have written the coordinates relative to the origin ().
Because the quantities and are the same for every particle, we can sum up the kinetic energy and angular momentum for all of the particles to get
where is the total mass, is the location of the centre of mass and is the moment of inertia of the body about the origin. The total angular momentum is
so the angular momentum is the generalized momentum conjugate to the angular velocity; if the potential does not depend on angle, the angular momentum is a first integral of the motion.
The moment of inertia is a positive definite matrix. We know this because when we multiply it by any non-zero vector , we get the kinetic energy of the body that is necessarily positive. Such a matrix can always be diagonalized. From a physical point of view, that means that the eigenvector of the matrix are known as the principal axes of the body. If the angular velocity points along one of the principal axes, then the angular momentum is parallel to the angular velocity. Furthermore, if one chooses the principal axes of the body to be one's coordinate axes, the diagonal element of the moment of inertia matrix vanish, so for example,
Sometimes it is easiest to calculate the moment of inertia about a particular point but you are interested in the moment of inertia about another point. There is a straightforward prescription to achieve this. First, let's use the center of mass of the system as the origin, so we have
Let's write this same equation with respect to a new origin so we have
If we compare this result to the earlier result we find
On a component by component basis we have
To change from one origin to another you first use the equation above to get the moment of inertia relative to the centre of mass and then use it a second time to get the moment of inertia about the new origin.
Let's calculate the moment of inertia of an ellipsoid. Let's assume that the axes of the ellipsoid are line up along the coordinate axes. Along the axis the ellipsoid spans from to . Along the and axes, the bounds are and respectively.
Because the coordinate axes are the principal axes of the ellipsoid, the off-diagonal components of the matrix will vanish. We will calculate the sum of over the ellipsoid first. We have
The region of integration is rather complicate in the Cartesian coordinates but this gives us a place to start. We would like to
simplify things a bit by defining
In these new coordinates the ellipsoid becomes a sphere of unit radius. We also have to include the Jacobian of the coordinate transformation
Now our integral looks like
Let's calculate the density of the ellipsoid in terms of its and volume. The volume of the ellipsoid is
so we can rewrite the density in the moment of inertia integral
The easiest way to integrate over a unit sphere is to use spherical coordinates. Let's do it:
By symmetry we can calculate the results for the other two coordinates. We have
and the moment of inertia matrix about the centre of mass have the following non-zero components:
This is an example of reducing a general shape to a symmetric one through a transformation of variables. One important result of this technique is rule.
We saw earlier that if the velocity of the origin vanished, then we could write
. Now if there are no forces on the body, the kinetic energy of the body is conserved. If there are no forces on the body, there are no torques, so the angular momentum is conserved as well. Does this mean that the angular velocity is constant?
The answer is of course no because the equation only says that the component of the angular velocity along the direction of the angular momentum is constant. The other components can change. For example, if the angular velocity lies along one of the principal axes (eigenvectors of the moment of inertia matrix), then the angular momentum points in the same direction and the angular velocity must stay constant in direction and magnitude.
A usual picture to understand the force-free motion of a rigid body is the ellipsoid of inertia. Let's use the principal axes of the body as a coordinate frame to write the kinetic energy
The ellipsoid of inertia is the locus of values of that give a kinetic energy of the body.
Let's calculate the shape of the ellipsoid of inertia of an ellipsoid. We have
so along the , and directions, the axes of the ellipsoid of inertia are
The size of the ellipsoid of inertia increases as and decreases with the mass and size of the actual ellipsoidal body. The shape of the ellipsoid of inertia reflects the shape of the physical ellipsoid. The largest axis of the ellipsoid is also the largest axis of the corresponding ellipsoid of inertia. If two axes of the ellipsoid are the same size, the corresponding axes of the ellipsoid of inertia will be equal as well.
To conserve energy the angular velocity must remain on the ellipsoid of inertia and the orientation of ellipsoid of inertia determines the orientation of the body. We can use the conservation of angular momentum to orient the ellipsoid. If we remember that we see that the normal to the ellipsoid of inertia is the angular momentum; therefore, as the ellipsoid of inertia rotates about the angular velocity, it must remain tangent to a plane perpendicular to the angular momentum, the invariable plane. Furthermore, the angular velocity is restricted to move along the ellipsoid of inertia along curves of constant . This curve is called the polhode, and it traces the path of the angular velocity through the body. The angular velocity also traces a path in the invariable plane called the herpolhode. The herpolhode is the path of the angular velocity through space.
We can also write that so the centre of the ellipsoid of inertia must remain a constant distance above the invariable plane. The height of the centre of the ellipsoid is , so when the ellipsoid is down low the angular momentum is high relative to the kinetic energy.
Keeping the ellipsoid of inertia tangent to the invariable plane preserves the direction of the angular momentum. The magnitude of the angular momentum is also conserved. This defines a second ellipsoid through the equation
so the polhodes are the intersections of the ellipsoid of inertia with the angular momentum ellipsoid. The angular momentum ellipsoid deviates more from a sphere than the ellipsoid of inertia; this gives an idea of where the polhodes lie for small and large angular momenta.
If the body has two principal axes with equal moments of inertia, the polhodes are circles centered on the axis with the unique moment of inertia, and the herpolhodes are circles in the invariable plane. If the unique moment of inertia is larger than the others, the body is called oblate. Otherwise it is prolate -- in analogy to the spheroids with the similar properties. In this case, the free rotation of the body consists of a constant precession of the angular velocity in a circle about the angular momentum in the space frame (the space cone) and about the unique principal axis in the body frame (the body cone). In the oblate case the body cone rolls within the space cone, and in the prolate case, the body cone rolls on the outside of the space cone.
If no two moments of inertia are equal to each other, the polhodes are much more complicated. The polhodes near the axes with the largest and smallest moment of inertia are closed, but those near the intermediate axis are not; consequently, rotation about the intermediate axis is not stable as we shall see when we analyze the free rotation using Euler's equations.
So far we have tried to find a graphical description of the motion of an object without any torques. Although "the polhode rolls on the herpolhode without slipping" might paint a nice Victorian picture, it doesn't exploit the mathematical language of physics to which we are now accustomed. We know that the change in the angular momentum is equal to the torques on the body
Let's write the angular momentum in terms of the principal axes of the body. Here we will be looking at how the principal axes move so we will use the unit vectors and to denote the current direction of the principal axes of the body. We can write
Let's define the angular momentum relative to the body axes as
to get
where the torque and angular velocity are written relative to the body axes. One can imagine using Euler's equations would be quite cumbersome if one had to include the torques because one would always have to transform the torques from the inertia frame to the body frame. Let's write out these equations component by component
These equations are most powerful when there are no torques, so the right-hand sides are zero. In this case we can see immediately that if any two moments of inertia differ than the angular velocity must lie along one of the principal axes to remain constant.
If two of the moments of inertia are equal to each other (let's take ), we can solve the torque-free Euler's equations exactly. We have
We can solve this with
with
yielding the solution
This result coincides with the graphical picture that the polhodes of a symmetric body are circles centered on the axis with the unique moment of inertia.
We can calculate the expected rate for the Earth to precess due to its oblatness. Let's estimate the ratio of the moments of inertia in terms of the radii of the Earth at the equator and at the poles, we have
where in the last step we have assume that the difference between the polar and equatorial radius is small. We have
, so the ratio is about . What is for the Earth? It is , so the period of the free precession of the Earth should be about 290 days. The period is somewhat longer about 433 days w:Chandler_wobble. You can check out a plot of the wobble at the Paris Observatory. Why the Earth wobbles is a bit of a puzzle since because the Earth is not a solid body, without a driving force it would have damped long ago. Since the period of the wobble is similar to a year, seasonal changes are thought to be to blame.
I first learned about the wobble while spending a summer at the Pulkovo Astronomical Observatory near St. Petersburg. I saw a plot of the location of the Earth's pole on the surface of the Earth -- it moves in a irregular circular shape with a radius of 3 to 15 meters.
The wobble can be measured by accurate measurements of the locations of stars relative to landmarks on the Earth. I asked my hosts why did they need to keep track the wobble so accurately. "To aim our missiles at your cities -- of course" came the reply.
Unfortunately, we can't find such a straightforward to the solution of the motion of a triaxial body. In the symmetric case, one of the components of the angular velocity was constant with time. This is not the case for a triaxial body. We have the following Euler's equations
Let's take and define the following three positive ratios of moments of inertia
Let's assume that the rotation is nearly about the major axis, so . We will only include the latter components of the angular velocity to first order, so we have
If we define we have the following solution for and ,
so the angular velocity travels in an ellipse centered on the major axis.
Let's assume that the rotation is nearly about the minor axis, so . We will only include the latter components of the angular velocity to first order, so we have
If we define we have the following solution for and ,
so the angular velocity travels in an ellipse centered on the minor axis.
Let's assume that the rotation is nearly about the intermediate axis, so . We will only include the latter components of the angular velocity to first order, so we have
If we define we have the following solution for and ,
so the angular velocity diverges exponetially away from the intermediate axis.
Although it is possible to include the effects of torques on the motion of a rigid body within the framework of Euler's equations, it is rather cumbersome. Here we will develop the Lagrangian treatment of a top -- a symmetric rigid body with a torque from the Earth's gravity.
In a Lagrangian treatment it is necessary to find a set of independent coordinates to describe the position of the system at any time. For rigid bodies a convenient set of such coordinates are called the Euler angles. There are several different conventions for the Euler angles. The convention used here makes the analysis of the top simpler.
To get from one set of coordinates to the other, we have to do a series of three rotations. We imagine that the primed coordinates are the body coordinates and the unprimed coordinates are the inertia coordinates; therefore, it is straightforward to write the kinetic energy in the primed coordinates and the potential energy in the unprimed coordinates. Here are the three steps in words and in terms of a rotation matrix
Rotate by an angle about the initial axis. The final position of the axis is called the line of nodes.
Rotate about the line of nodes by an angle . This changes the direction of the axis and pitches the axis out of the initial plane.
Finally rotate about the new axis by an angle . This takes the axis out of the initial plane.
Euler proved that any rotation can be decomposed into these three rotations. One can see this mathematically by constructing the rotation matrix for these three successive transformations () and verifying that it is the most general rotation matrix.
To calculate the kinetic energy of the body, we have to find the angular velocity in terms of the coordinates, and and their time derivatives. In any coordinate system, the angular velocity is directed along the axis of rotation. We can construct the total angular velocity vector using the following information:
Changes in produce rotations about the axis (inertial frame).
Changes in produce rotations about the axis (body frame).
Changes in produce rotations about the line of nodes.
We can write the angular velocity relative to the inertial frame to get
or relative to the body frame
The first form is useful to visualize what the body is doing, while the second form is good for calculating the kinetic energy because we have used the primed coordinate system to characterize the location of the principal axes of the body.
For simplicity let's specialize to symmetric bodies with , so we have
Because the primed frame is fixed to the body, the moments of inertia in the kinetic energy are constant with time. Let's first calculate
This simplication only obtains for a symmetric body. The component is
The Lagrangian is the difference between the kinetic and the potential energy, so we have
We notice immediately that and don't appear in the Lagrangian. The system looks the same as in turns around the and axes. This means that we have two conserved quantities, reducing the three dimensional problem to a single dimension (we will choose this to be ). Because the Lagrangian does not depend on time, the Hamiltonian is conserved as well so we can use the techniques outlined in Linear Motion to understand the top.
To calculate the conserved momenta we have to calculate the partial derivatives of the Lagrangian with respect to the generalized velocities. We have
where the last step defines the conserved quantity that has units of angular velocity. Similarly for the coordinate , we have
Finally we have the Hamiltonian that is also conserved,
As we mentioned earlier, these three integrals are sufficient to solve for the motion. We can use the two conserved momenta to solve for the generalized velocities, and , in terms of the conserved quantities and and the value of the third coordinate . We have
The equation for is cubic so it will generally have three roots (designated at -- it must have at least one root because for large positive values of is positive and
for large negative values of is negative. As the top moves from straight up to straight down, ranges from 1 to -1, so only the roots within this range are important. Furthermore, the equation for has a root at .
The figure shows a typical curve for as a function of . The motion is bounded between and . Depending on the value of , the motion in the -direction can be monotonic or not. If , the motion will have a cusp at the upper end of the range.
The figures to the right show the motion of the axis as a function of time for several values of with the curve for above. The bounds of the motion in all three cases is from to . How can we set up a top to execute these three classes of motion? If we release the axis of the top from a stationary position, the axis will first drop (nutation) because of gravity. As it drops the angular momentum about the vertical axis must be conserved so the tip of the top starts moving sideways as well (precession) -- yielding the motion in the middle figure. On the other hand if we release the tip of the top with some motion in the direction of the precession, we will get the lower figure. Finally to get the rather freaky looking motion in the top figure, we release the top with some motion opposing the direction of the precession.
Specifically let's examine the case where we release the top from a stationary position, so where is the middle root of the curve. Initially we have
Unfortunately, the motion of the top in general cannot be solved in closed form, but we can make some progress in understanding the motion of the so-called fast top. If we look at the expression above, we have the ratio
The numerator is the potential energy of the top while the denominator is the related to the kinetic energy of the top about its spin axis. The condition for a fast top is
In a fast top, the kinetic energy is much larger than the potential energy; consequently, the extent of the nutation () is small, so
so we can write
The extent of the nutation is
and we can rewrite the equation for as
Let's take the time derivative of both sides to yield
A sleeping top is a top that is set in motion with its spin axis vertical and the axis remains vertical. When the top is vertical it is impossible to distinguish changes in from changes in so the relationship between and and isn't very useful. Let's look at the equation for . We have
When the top points upwards, . We let go of the top with so we have