## 2.10  Axisymmetric Tops

We have all played with a top at one time or another. For the purposes of analysis we will consider an idealized top that does not wander around. Thus, an ideal top is a rotating rigid body, one point of which is fixed in space. Furthermore, the center of mass of the top is not at the fixed point, which is the center of rotation, and there is a uniform gravitational acceleration.

For our top we can take the Lagrangian to be the difference of the kinetic energy and the potential energy. We already know how to write the kinetic energy -- what is new here is that we must express the potential energy in terms of the configuration. In the case of a body in a uniform gravitational field this is easy. The potential energy is the sum of ``mgh'' for all the constituent particles:

where g is the gravitational acceleration, h = · , and the unit vector indicates which way is up. Rewriting the vector to the constituents in terms of the vector to the center of mass, the potential energy is

where the last sum is zero because the center of mass is the origin of . So the potential energy of a body in a gravitational field with uniform acceleration is very simple: it is just Mgh, where M is the total mass and h = · is the height of the center of mass.

Here we consider an axisymmetric top (see figure 2.4). Such a top has an axis of symmetry of the mass distribution, so the center of mass is on the symmetry axis and the fixed point is also on the axis of symmetry.

In order to write the Lagrangian we need to choose a set of generalized coordinates. If we choose them well we can take advantage of the symmetries of the problem. If the Lagrangian does not depend on a particular coordinate, the conjugate momentum is conserved, and the complexity of the system is reduced.

The axisymmetric top has two apparent symmetries. The fact that the mass distribution is axisymmetric implies that neither the kinetic nor potential energy is sensitive to the orientation of the top about that symmetry axis. Additionally, the kinetic and potential energy are insensitive to a rotation of the physical system about the vertical axis, because the gravitational field is uniform.

We can take advantage of these symmetries by choosing appropriate coordinates, and we already have a coordinate system that does the job -- the Euler angles.12 We choose the reference orientation so that the symmetry axis is vertical. The first Euler angle, , expresses a rotation about the symmetry axis. The next Euler angle, , is the tilt of the symmetry axis of the top from the vertical. The third Euler angle, , expresses a rotation of the top about the z axis. The symmetries of the problem imply that the first and third Euler angles do not appear in the Lagrangian. As a consequence the momenta conjugate to these angles are conserved quantities. Let's work out the details.

First, we develop the Lagrangian explicitly. The general form of the kinetic energy about a point is given by equation 2.30. The top is constrained so that it pivots about a fixed point that is not at the center of mass. So the moments of inertia that enter the kinetic energy are the moments of inertia of the top with respect to the pivot point, not with respect to the center of mass. If we know the moments of inertia about the center of mass we can write the moments of inertia about the pivot in terms of them (see exercise 2.2). So let's assume the principal moments of inertia of the top about the pivot are A, B, and C, and A = B because of the symmetry.13 We can use the computer to help us figure out the Lagrangian for this special case:

(show-expression
((T-rigid-body 'A 'A 'C)
(up 't
(up 'theta 'phi 'psi)
(up 'thetadot 'phidot 'psidot))))

We can rearrange this a bit to get

In terms of Euler angles, the potential energy is

where R is the distance of the center of mass from the pivot. The Lagrangian is L = T - V. We see that the Lagrangian is indeed independent of and , as expected.

There is no particular reason to look at the Lagrange equations. We can assign that job to the computer when needed. However, we have already seen that it may be useful to examine the conserved quantities associated with the symmetries.

The energy is conserved, because the Lagrangian has no explicit time dependence. Also, the energy is the sum of the kinetic and potential energy E = T + V, because the kinetic energy is a homogeneous quadratic form in the generalized velocities. The energy is

Two of the generalized coordinates do not appear in the Lagrangian, so there are two conserved momenta. The momentum conjugate to is

The momentum conjugate to is

The state of the system at a moment is specified by the tuple ( t; , , ; , , ). The two coordinates and that do not appear in the Lagrangian do not appear in the Lagrange equations or the conserved momenta. So the evolution of the remaining four state variables, , , , and , depends only on those remaining state variables. This subsystem for the top has a four-dimensional state space. The variables that did not appear in the Lagrangian can be determined by integrating the derivatives of these variables, which are determined separately by solving the independent subsystem.

The evolution of the top is described by a four-dimensional subsystem and two auxiliary quadratures.14 This subdivision is a consequence of choosing generalized coordinates that incorporate the symmetries. However, the choice of generalized coordinates that incorporate the symmetries also gives conserved momenta. We can make use of these momenta to simplify the formulation of the problem further. Each integral can be used to locally eliminate one dimension of the subsystem. In this case the subsystem has four dimensions and there are three integrals, so the system can be completely reduced to quadratures. For the top, this can be done analytically, but we think it is a waste of time to do so. Rather, we are interested in extracting interesting features of the motion. We concentrate on the energy integral and use the two conserved momenta to eliminate and . After a bit of algebra we find:

Along a path , where D(t) is substituted for , this is an ordinary differential equation for . This differential equation involves various constants, some of which are set by the initial conditions of the other state variables. The solution of the differential equation for involves no more than ordinary integrals. So the top is essentially solved. We could continue this argument to obtain the qualitative behavior of : Using the energy (2.66), we can plot the trajectories in the plane of versus and see that the motion of is simply periodic. However, we will defer continuing along this path until chapter 3, when we have developed more tools for analysis.

Let's get real. Let's make a top out of an aluminum disk with a steel rod through the center to make the pivot. Measuring the top very carefully, we find that the moment of inertia of the top about the symmetry axis is about 6.60 × 10-5 kg m2, and the moment of inertia about the pivot point is about 3.28 × 10-4 kg m2. The combination gMR is about 0.0456 kg m2 s-2. We spin the top up with an initial angular velocity of = 140 rad s-1 (about 1337 rpm). The top initially has = = = 0 and is initially tilted with = 0.1 rad. We then kick it so that = - 15 rad s-1. Figures 2.5-2.8 display aspects of the evolution of the top for 2 seconds. The tilt of the top (measured by ) varies in a periodic manner. The orientation about the vertical is measured by : we see that the top also precesses, and the rate of precession varies with . We also see that as the top bobs up and down the rate of rotation of the top oscillates -- the top spins faster when it is more vertical. The plot of tilt versus precession angle shows that in this case the top executes a looping motion. If we do not kick it but just let it drop, then the loop disappears, leaving just a cusp. If we kick it in the other direction, then there is no cusp nor any looping motion.

Exercise 2.11.  Kinetic energy of the top
The rotational kinetic energy of the top can be written in terms of the principal moments of inertia with respect to the pivot point and the angular velocity vector of rotation with respect to the pivot point. Show that this formulation of the kinetic energy yields the same value that one would obtain by computing the sum of the rotational kinetic energy about its center of mass and the kinetic energy of the motion of the center of mass.

Exercise 2.12.  Nutation of the top

a.  Carry out the algebra to obtain the energy (2.66) in terms of and .

b.  Numerically integrate the Lagrange equations for the top to obtain figure 2.5, versus time.

c.  Note that the energy is a differential equation for in terms of , with conserved quantities p, p, and E determined by initial conditions. Can we use this differential equation to obtain as a function of time? Explain.

Exercise 2.13.  Precession of the top
Consider a top that is rotating so that is constant.

a.  Using the angular momentum integrals, compute the rate of precession as a function of the conserved angular momenta and the equilibrium value of .

b.  For to be at an equilibrium the acceleration D2 must be zero. Use the Lagrange equation for to find the rate of precession at the equilibrium in terms of the equilibrium and .

c.  Find an approximate expression for the precession rate in the limit that is large.

d.  The Newtonian rule is that the rate of change of the angular momentum is the torque. Assume the top is spinning so fast that the angular momentum is nearly the same as the angular momentum of rotation about the symmetry axis. By equating the rate of change of this vector angular momentum to the gravitational torque on the center of mass develop an approximate formula for the precession rate.

e.  Numerically integrate the top to check your deductions.

12 That the axisymmetric top can be solved in Euler angles is, no doubt, the reason for the traditional choice of the definition of these. For other problems, the Euler angles may offer no particular advantage.

13 Here, we do not require that C be larger than A = B, because they are not measured with respect to the center of mass.

14 Traditionally, evaluating a definite integral is known as performing a quadrature.