# Advanced Classical Mechanics/Rigid Bodies

## Definition[edit | edit source]

- In a rigid body, no part of the body moves relative to another part of the body regardless of the force applied.
- A rigid body is equivalent to a system of particles restricted to have fixed distances apart.

There are two methods for solving for the motion of rigid bodies.

- Write out the kinetic energy and use the Lagrangian to derive the equations of motion
- Use Euler's equations to get the motion of the center of mass and the rotation.

## The Motion of a Rigid Body[edit | edit source]

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 .

### Angular Velocity[edit | edit source]

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.

## Moment of Inertia[edit | edit source]

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 ().

### Kinetic Energy[edit | edit source]

Its kinetic energy is

and multiplying out the terms we get

We can simplify the result using the following fact

to get

### Angular Momentum[edit | edit source]

Let's look at the angular momentum of the particle, we have

Here we can use the vector triple product to simplify the result

which for us yields

If then

In this case, let's write out component by component, we have

defining the moment of inertia matrix for that particular particle.

### Putting It All Together[edit | edit source]

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

Using the result for the scalar triple product we can show that

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.

### Manipulating the Moment of Inertia[edit | edit source]

#### Principal Axes[edit | edit source]

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,

and the kinetic energy is

#### Change of Origin[edit | edit source]

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.

### Calculating a Moment of Inertia[edit | edit source]

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.

## Ellipsoid of Inertia[edit | edit source]

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.

### Ellipsoid of Inertia of an Ellipsoid[edit | edit source]

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.

### Motion of the Ellipsoid of Inertia[edit | edit source]

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.

## Euler's Equations[edit | edit source]

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.

### Chandler Wobble[edit | edit source]

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.

### Triaxial Body[edit | edit source]

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

giving the following equations

#### Major Axis[edit | edit source]

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.

#### Minor Axis[edit | edit source]

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.

#### Intermediate Axis[edit | edit source]

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.

## Lagrangian Treatment of Rigid-Body Motion[edit | edit source]

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.

### Euler Angles[edit | edit source]

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.

### The Kinetic Energy[edit | edit source]

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

Combining these results yields the kinetic energy

### The Potential Energy[edit | edit source]

The potential energy is pretty simple using the Euler angles. If is the distance between the pivot point and the center of mass of the top, we have

### The Lagrangian[edit | edit source]

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.

#### Integrals of the Motion[edit | edit source]

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

#### Conservation of Energy[edit | edit source]

Since both the Hamiltonian and the angular velocity along the axis are conserved, let's define the conserved quantity

which we can rewrite as

where

We can make further progress by looking at which will simplify the potential. Let's write

so we can write in terms of and its time derivative

What remains is to solve for to get the boundaries of the motion

### The Motion of a Symmetric Top[edit | edit source]

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

so we can write

#### The Fast Top[edit | edit source]

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

Let's define so we have

with the solution

and

The frequency of the nutation is

and

and

#### Sleeping Top[edit | edit source]

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

so

and

Using the definition of we find

and

so there is a double root at and another root at

If this root is at , then the top will sleep, so

is the condition for a sleeping top.