# 25. Rigid Body Moving Freely

Michael Fowler

## Angular Momentum and Angular Velocity

In contrast to angular velocity, the angular momentum of a body depends on the point with respect to which it is defined.  For now, we take it (following Landau, of course) as relative to the center of mass, but we denote it by $\stackrel{\to }{L},$ following modern usage.  This “intrinsic” angular momentum is like the Earth’s angular momentum from its diurnal rotation, as distinct from its orbital angular momentum in going around the Sun.

That is,

$\stackrel{\to }{L}=\sum _{n}{\stackrel{\to }{r}}_{n}×{m}_{n}{\stackrel{\to }{v}}_{n}=\sum _{n}{\stackrel{\to }{r}}_{n}×{m}_{n}\left(\stackrel{\to }{\Omega }×{\stackrel{\to }{r}}_{n}\right)=\sum _{n}{m}_{n}\left[{r}_{n}^{2}\stackrel{\to }{\Omega }-{\stackrel{\to }{r}}_{n}\left({\stackrel{\to }{r}}_{n}\cdot \stackrel{\to }{\Omega }\right)\right]=I\stackrel{\to }{\Omega }$ ,

where $I$ is the inertia tensor: this just means ${L}_{i}={I}_{ik}{\Omega }_{k}$.

Explicitly, taking the principal axes as the $\left({x}_{1},{x}_{2},{x}_{3}\right)$ axes,

${L}_{1}={I}_{1}{\Omega }_{1},\text{ }\text{ }{L}_{2}={I}_{2}{\Omega }_{2},\text{ }\text{ }{L}_{3}={I}_{3}{\Omega }_{3}.\text{ }\text{ }$  For anything with spherical inertial symmetry (such as a cube or a tetrahedron!) $\stackrel{\to }{L}=I\stackrel{\to }{\Omega }.$

Landau defines a rotator as a collection of massive particles all on a line.  (I guess that includes diatomic molecules, and, for example, CO2, neglecting electrons and nuclear size).  We know there are only two physical rotational degrees of freedom for these molecular rotators (thanks to quantum mechanics) and obviously the two principal axes are perpendicular to the line of masses, and degenerate.  Again, then, $\stackrel{\to }{L}=I\stackrel{\to }{\Omega }.$

## Precession of a Symmetrical Top

A more interesting case is the free rotation (zero external torque) of a symmetrical top, meaning ${I}_{1}={I}_{2}\ne {I}_{3}.$ We can take any pair of orthogonal axes, perpendicular to the body’s symmetry axis, as the ${x}_{1},{x}_{2}$ axes.  We’ll choose ${x}_{2}$, following Landau, as perpendicular to the plane containing $\stackrel{\to }{L}$ and the momentary position of the ${x}_{3}$ axis, so in the diagram here ${x}_{2}$ is perpendicularly out from the paper/screen, towards the viewer.

This means the angular momentum component ${L}_{2}=0$  and therefore ${\Omega }_{2}=0$.  Hence $\stackrel{\to }{\Omega }$ is in the same plane as $\stackrel{\to }{L},{x}_{3}$, and so the velocity $\stackrel{\to }{v}=\stackrel{\to }{\Omega }×\stackrel{\to }{r}$ of every point on the axis of the top is perpendicular to this plane (into the paper/screen).  The axis of the top $O{x}_{3}$ must be rotating uniformly about the direction of $\stackrel{\to }{L}.$

The spin rate of the top around its own axis is

${\Omega }_{3}={L}_{3}/{I}_{3}=\left(L/{I}_{3}\right)\mathrm{cos}\theta .$ The angular velocity vector $\stackrel{\to }{\Omega }$ can be written as a sum of two components, one along the body’s axis $O{x}_{3}$ and one parallel to the angular momentum $\stackrel{\to }{L}$ (these components are shown dashed in the figure)

$\stackrel{\to }{\Omega }={\stackrel{\to }{\Omega }}_{\text{precession}}+{\stackrel{\to }{\Omega }}_{3}$  The component along the body’s axis $O{x}_{3}$ does not contribute to the precession, which all comes from the component along the (fixed in space) angular momentum vector.

The speed of precession follows from ${\Omega }_{\text{precession}}\mathrm{sin}\theta ={\Omega }_{1}$ and ${\Omega }_{1}={L}_{1}/{I}_{1}=\left(L/{I}_{1}\right)\mathrm{sin}\theta ,$ so ${\Omega }_{\text{precession}}=L/{I}_{1}.$Note also the ratio of precession rate to spin around axis is

${\Omega }_{\text{precession}}/{\Omega }_{3}=\left({I}_{3}/{I}_{1}\right)\mathrm{sec}\theta .$ This means the precession rate and the spin are very comparable, except when $\theta$ is near $\pi /2,$ when the precession becomes much faster. Remember this is the body’s precession with no external torque, and is clearly completely different$—$much faster precession$—$than the familiar case of a fast spinning top under gravity.

## Throwing a Football…

If you throw a football and manage to give it only spin about the long axis, it will stay pointing that way (apart from drag effects, which will tend to line up spin direction with velocity).  If when you throw it you also add some angular velocity along a shorter axis, it will precess (wobble). Given the angle, the ratio of precession to spin is fixed.

If you take a ball that is already spinning fast about its long axis, have the two ends of the long axis (its points) in your palms, then, as you throw it, give it a quick twist by moving one hand downwards and the other up as you throw, to give it significant angular velocity about a short axis, at the same time keeping the fast spin about the long axis, once the ball leaves your hands, the angle between the spin and the total angular momentum, the angle of wobble, is completely determined by the ratio of the two angular velocities.

## Equations of Motion for Rigid Body with External Forces

### Translation

A free rigid body has six degrees of freedom (for instance, the coordinates of the center of mass and the orientation of the body).  Therefore, there are six equations of motion, three for the rate of change of spatial position of the center of mass, in other words for the components of the velocity $\stackrel{\to }{V},$ and three for the rate of change of orientation, the angular velocity $\stackrel{\to }{\Omega }$.

These equations are of course nothing but Newton’s laws, easily derived by summing over the set of equations ${\stackrel{\to }{f}}_{i}=d\left({m}_{i}{\stackrel{\to }{v}}_{i}\right)/dt$ for each particle.

Denoting the total momentum of the body by $\stackrel{\to }{P},$

$\sum _{n}\frac{d}{dt}\left({m}_{n}{\stackrel{\to }{v}}_{n}\right)=\frac{d\stackrel{\to }{P}}{dt}=\sum _{n}{\stackrel{\to }{f}}_{n}=\stackrel{\to }{F}$and $\stackrel{\to }{P}=M\stackrel{\to }{V},$ where $\stackrel{\to }{V}=d\stackrel{\to }{R}/dt$ is the velocity of the center of mass. (This can be established by differentiating with respect to time the definition of the center of mass, $M\stackrel{\to }{R}=\sum _{n}{m}_{n}{\stackrel{\to }{r}}_{n}.$ )

The total force on all the particles is a sum of the total external force on the body and the sum of internal forces between particles$—$but these internal forces come in equal and opposite pairs, from Newton’s Third Law, and therefore add to zero.

The bottom line, then, is that the rate of change of momentum of a rigid body equals the total external force on the body.  If this force is from a time-independent potential, then

$\stackrel{\to }{F}=-\text{ }\text{ }\partial V/\partial \stackrel{\to }{R}$ because if the body is moved through $\delta \stackrel{\to }{R}$ (without rotation, hence the partial derivative), each individual particle moves through the same $\delta \stackrel{\to }{R}$, the work done by the external potential on the ${n}^{\text{th}}$ particle is ${\stackrel{\to }{f}}_{n}^{\text{ext}}\cdot \delta \stackrel{\to }{R}=-\delta {V}_{n}$ , and summing over all the particles gives $\stackrel{\to }{F}\cdot \delta \stackrel{\to }{R}=-\delta {V}_{\text{tot}},$ giving the above equation as $\delta \stackrel{\to }{R}\to 0.$

### Rotation

To derive the equation of motion for rotation of a rigid body, we choose the inertial frame in which the center of mass is momentarily at rest, and take the center of mass as the origin.

The rate of change of angular momentum about the center of mass (origin),

$\stackrel{˙}{\stackrel{\to }{L}}=\left(d/dt\right)\sum _{n}{\stackrel{\to }{r}}_{n}×{\stackrel{\to }{p}}_{n}=\sum _{n}\left[\left({\stackrel{˙}{\stackrel{\to }{r}}}_{n}×{\stackrel{\to }{p}}_{n}\right)+\left({\stackrel{\to }{r}}_{n}×{\stackrel{˙}{\stackrel{\to }{p}}}_{n}\right)\right]=\sum _{n}{\stackrel{\to }{r}}_{n}×{\stackrel{\to }{f}}_{n}=\stackrel{\to }{K}$ where we dropped the ${\stackrel{˙}{\stackrel{\to }{r}}}_{n}×{\stackrel{\to }{p}}_{n}$ term because ${\stackrel{˙}{\stackrel{\to }{r}}}_{n}={\stackrel{\to }{v}}_{n}$ is parallel to ${\stackrel{\to }{p}}_{n}=m{\stackrel{\to }{v}}_{n},$ then we used ${\stackrel{\to }{f}}_{n}={\stackrel{˙}{\stackrel{\to }{p}}}_{n}$ to get the total moment of the external forces about the center of mass, the torque.

The angular momentum about the center of mass is the same in any inertial frame, since the extra term on adding a velocity ${\stackrel{\to }{v}}_{0}$ to each mass is $\sum {\stackrel{\to }{r}}_{n}×{m}_{n}{\stackrel{\to }{v}}_{0}=\text{ }-{\stackrel{\to }{v}}_{0}×\sum {m}_{n}{\stackrel{\to }{r}}_{n}=0$ from the definition of the center of mass.

If the center of mass is not at the origin, denote the particle coordinates by ${\stackrel{\to }{\rho }}_{n}=\stackrel{\to }{R}+{\stackrel{\to }{r}}_{n}$ in the usual notation, so

a sum of an intrinsic (“spin”) angular momentum and an extrinsic (“orbital”)angular momentum.

Similarly, if the torque of external forces relative to the center of mass is $\stackrel{\to }{K}=\sum _{n}{\stackrel{\to }{r}}_{n}×{\stackrel{\to }{f}}_{n}$ as defined above, then relative to the new origin the torque is

that is, the torque about the new origin is the torque about the center of mass plus the torque about the new origin of the total external force acting at the center of mass.

An important special case is that of a couple: a pair of equal but oppositely directed forces, acting along parallel but separated lines (like two hands oppositely placed turning a steering wheel).  The forces add to zero, so from the above equation a couple exerts the same torque about any origin.

More generally, the term couple is often used (including by Landau) to refer to any set of forces that add to zero, but give a nonzero torque because of their lines of action, and such a set give the same torque about any origin.

Exercise: prove that for a rigid body freely falling in a uniform gravitational field, the angular momentum about the center of mass remains constant, but about another point it will in general be changing.   What about a charged rigid body moving in space (no gravity) through a uniform electric field?