26. Euler’s Angles

Michael Fowler

Introduction

So far, our analysis of rotational motion has been of essentially one dimensional, or more precisely one angular parameter, motion: rotating about an axis, rolling, precessing and so on. But this leaves out many interesting phenomena, for example the wobbling of a slowing down top, nutation, and so on. We need a well-defined set of parameters for the orientation of a rigid body in space to make further progress in analyzing the dynamics.

The standard set is Euler’s Angles. What you see as you watch a child’s top beginning to wobble as it slows down is the direction of the axis this is given by the first two of Euler’s angles: θ,ϕ  the usual spherical coordinates, the angle θ  from the vertical direction and the azimuthal angle ϕ  about that vertical axis.  Euler’s third angle, ψ,  specifies the orientation of the top about its own axis, completing the description of the precise positioning of the top. To describe the motion of the wobbling top as we see it, we evidently need to cast the equations of motion in terms of these angles.

Definition

The rotational motion of a rigid body is completely defined by tracking the set of principal axes x 1 , x 2 , x 3 , with origin at the center of mass, as they turn relative to a set of fixed axes X,Y,Z .  The principal axes can be completely defined relative to the fixed set by three angles: the two angles θ,ϕ  fix the direction of x 3 ,  but that leaves the pair x 1 , x 2  free to turn in the plane perpendicular to x 3 ,  the angle ψ  fixes their orientation. 

To see these angles, start with the fixed axes, draw a circle centered at the origin in the horizontal X,Y  plane. Now draw a circle of the same size, also centered at the same origin, but in the principal axes x 1 , x 2  plane.  Landau calls the line of intersection of these circles (or discs) the line of nodes. It goes through the common origin, and is a diameter of both circles.

The angle between these two planes, which is also the angle between Z, x 3  (since they’re the perpendiculars to the planes) is labeled θ.  

The angle between this line of nodes and the X  axis is ϕ.  It should be clear that θ,ϕ  together fix the direction of x 3 ,  then the other axes are fixed by giving ψ,  the angle between x 1  and the line of nodes ON.  The direction of measurement of ϕ,ψ  around Z, x 3  are given by the right-hand or corkscrew rule.

Angular Velocity and Energy in Terms of Euler’s Angles

Since the position is uniquely defined by Euler’s angles, angular velocity is expressible in terms of these angles and their derivatives.  

The strategy here is to find the angular velocity components along the body axes x 1 , x 2 , x 3  of θ ˙ , ϕ ˙ , ψ ˙  in turn.  Once we have the angular velocity components along the principal axes, the kinetic energy is easy.

(You might be thinking: wait a minute, aren't the axes embedded in the body? Don't they turn with it?  How can you talk about rotation about these axes? Good point: what we're doing here is finding the components of angular velocity about a set of axes fixed in space, not the body, but momentarily coinciding with the principal axes of the body.)

From the diagram, θ ˙  is along the line ON,  and therefore in the x 1 , x 2  plane: notice it is at an angle ψ  with respect to x 1 .  Its components are therefore θ ˙ = θ ˙ cosψ, θ ˙ sinψ,0 .  

Now ϕ ˙  is about the Z  axis. The principal axis x 3  is at angle θ  to the Z  axis, so ϕ ˙  has component ϕ ˙ cosθ  about x 3 ,  and ϕ ˙ sinθ  in the x 1 , x 2  plane, that latter component along a line perpendicular to ON , and therefore at angle ψ  from the x 2  axis. Hence ϕ ˙ = ϕ ˙ sinθsinψ, ϕ ˙ sinθcosψ, ϕ ˙ cosθ .  

The angular velocity ψ ˙  is already along a principal axis, x 3 .  

To summarize, the Euler angle angular velocities (components along the body’s principal axes) are:

θ ˙ = θ ˙ cosψ, θ ˙ sinψ,0 , ϕ ˙ = ϕ ˙ sinθsinψ, ϕ ˙ sinθcosψ, ϕ ˙ cosθ , ψ ˙ = 0, 0, ψ ˙ ,  from which, the angular velocity components along those in-body axes x 1 , x 2 , x 3  are:

Ω 1 = ϕ ˙ sinθsinψ+ θ ˙ cosψ, Ω 2 = ϕ ˙ sinθcosψ θ ˙ sinψ, Ω 3 = ϕ ˙ cosθ+ ψ ˙ .  For a symmetric top, meaning I 1 = I 2 I 3 ,  the rotational kinetic energy is therefore

T rot = 1 2 I 1 Ω 1 2 + Ω 2 2 + 1 2 I 3 Ω 3 2 = 1 2 I 1 ϕ ˙ 2 sin 2 θ+ θ ˙ 2 + 1 2 I 3 ϕ ˙ cosθ+ ψ ˙ 2 .  For this symmetrical case, as Landau points out, we could have taken the x 1  axis momentarily along the line of nodes ON,  giving

Ω = θ ˙ , ϕ ˙ sinθ, ϕ ˙ cosθ+ ψ ˙ .

Free Motion of a Symmetrical Top

As a warm up in using Euler’s angles, we’ll redo the free symmetric top covered in the last lecture. With no external torques acting the top will have constant angular momentum L .  

We’ll take L  in the fixed Z  direction. The axis of the top is along x 3 .  

 Taking the x 1  axis along the line of nodes ON  (see figure below) at the instant considered, the constant angular momentum L = I 1 Ω 1 , I 1 Ω 2 , I 3 Ω 3 = I 1 θ ˙ , I 1 ϕ ˙ sinθ, I 3 ϕ ˙ cosθ+ ψ ˙ .

Remember, this new x 1  axis (see diagram!) is perpendicular to the Z  axis we’ve taken L  along, so L 1 = I 1 θ ˙ =0,  and θ  is constant, meaning that the principal axis x 3  describes a cone around the constant angular momentum vector L .  The rate of precession follows from the constancy of L 2 = I 1 ϕ ˙ sinθ.   Writing the absolute magnitude of the angular momentum as L,   L 2 =Lsinθ,  (remember L  is in the Z  direction, and x 1  is momentarily along ON  )  so the rate of precession ϕ ˙ =L/ I 1 .  Finally, the component of L  along the x 3  axis of symmetry of the top is Lcosθ= I 3 Ω 3 ,  so the top’s spin along its own axis is Ω 3 = L/ I 3 cosθ.  

Motion of Symmetrical Top around a Fixed Base with Gravity:  Nutation

Denoting the distance of the center of mass from the fixed bottom point P  as  (along the axis) the moment of inertia about a line perpendicular to the axis at the base point is

I 1 = I 1 +M 2 .

( I 1  being usual center of mass moment.)

The Lagrangian is ( P  being the origin, I 3  in direction θ,ϕ  ) L= 1 2 I 1 ϕ ˙ 2 sin 2 θ+ θ ˙ 2 + 1 2 I 3 ϕ ˙ cosθ+ ψ ˙ 2 Mgcosθ.  Notice that the coordinates ψ,ϕ  do not appear explicitly, so there are two constants of motion:

p ψ =L/ ψ ˙ = I 3 ϕ ˙ cosθ+ ψ ˙ = L 3 , p ϕ =L/ ϕ ˙ = I 1 sin 2 θ+ I 3 cos 2 θ ϕ ˙ + I 3 ψ ˙ cosθ= L Z .  That is, the angular momentum about x 3  is conserved, because the two forces acting on the top, the gravitational pull at the center of mass and the floor reaction at the bottom point, both act along lines intersecting the axis, so never have torque about x 3 .  The angular momentum about Z  is conserved because the gravitational torque acts perpendicular to this line.

We have two linear equations in ψ ˙ , ϕ ˙  with coefficients depending on θ  and the two constants of motion L 3 , L Z .  The solution is straightforward, giving

ϕ ˙ = L Z L 3 cosθ I 1 sin 2 θ , ψ ˙ = L 3 I 3 cosθ L Z L 3 cosθ I 1 sin 2 θ .  The (conserved) energy

E= 1 2 I 1 Ω 1 2 + Ω 2 2 + 1 2 I 3 Ω 3 2 +Mgcosθ = 1 2 I 1 ϕ ˙ 2 sin 2 θ+ θ ˙ 2 + 1 2 I 3 ϕ ˙ cosθ+ ψ ˙ 2 +Mgcosθ.  Using the constants of motion to express ψ ˙ , ϕ ˙  in terms of θ  and the constants L Z , L 3  , then subtracting a θ  independent term to reduce clutter, E =EMg L 3 2 /2 I 3 ,  we have

E = 1 2 I 1 θ ˙ 2 + V eff θ , V eff θ = L Z L 3 cosθ 2 2 I 1 sin 2 θ Mg 1cosθ .  The range of motion in θ  is given by E > V eff θ .  For L 3 L Z ,   V eff θ  goes to infinity at θ=0,π . It has a single minimum between these points.  (This isn’t completely obvious one way to see it is to change variable to u=cosθ , following Goldstein. Multiplying throughout by sin 2 θ , and writing θ ˙ 2 sin 2 θ= u ˙ 2  gives a one dimensional particle in a potential problem, and the potential is a cubic in u.  Of course some roots of  E = V eff θ  could be in the unphysical region u >1.  In any case, there are at most three roots, so since the potential is positive and infinite at θ=0,π,  it has at most two roots in the physical range.)

From the one-dimensional particle in a potential analogy, it’s clear that θ  oscillates between these two points θ 1  and θ 2 .  This oscillation is called nutation.  Now ϕ ˙ = L Z L 3 cosθ / I 1 sin 2 θ  could change sign during this oscillation, depending on whether or not the angle cos 1 L Z / L 3  is in the range. Visualizing the path of the top center point on a spherical surface centered at the fixed point, as it goes around it oscillates up and down, but if there is this sign change, it will “loop the loop”, going backwards on the top part of the loop.

Steady Precession

Under what conditions will a top, spinning under gravity, precess at a steady rate?  The constancy of L 3 , L Z  mean that Ω 3 = ϕ ˙ cosθ+ ψ ˙ , and  Ω pr = ϕ ˙  are constants.

The θ  Lagrange equation is

I 1 θ ¨ = I 1 ϕ ˙ 2 sinθcosθ I 3 ϕ ˙ cosθ+ ψ ˙ ϕ ˙ sinθ+Mgsinθ

For constant θ,   θ ¨ =0,  so, with Ω 3 = ϕ ˙ cosθ+ ψ ˙ , and  Ω pr = ϕ ˙ ,   

I 1 Ω pr 2 cosθ I 3 Ω 3 Ω pr +Mg=0.

Since this is a quadratic equation for the precession rate, there are two solutions in general: on staring at a precessing top, this is a bit surprising!  We know that for the top, when it’s precessing nicely, the spin rate Ω 3  far exceeds the precession rate Ω pr .   Assuming I 1 , I 3  to be of similar size, this means the first term in the quadratic is much smaller than the second. If we just drop the first term, we get the precession rate

Ω precess (slow) = Mg I 3 Ω 3 , Ω 3 Ω precess .  Note that this is independent of angle the torque varies as sinθ,  but so does the horizontal component of the angular momentum, which is what’s changing.

This is the familiar solution for a child’s fast-spinning top precessing slowly.  But this is a quadratic equation, there’s another possibility: in this large Ω 3  limit, this other possibility is that Ω pr  is itself of order Ω 3 , so now in the equation the last term, the gravitational one, is negligible, and

Ω precess (fast) I 3 Ω 3 / I 1 cosθ.

This is just the nutation of a free top!  In fact, of course, both of these are approximate solutions, only exact in the limit of infinite spin (where one goes to zero, the other to infinity), and a more precise treatment will give corrections to each arising from the other.  Landau indicates the leading order gravitational correction to the free body nutation mode.

Stability of Top Spinning about Vertical Axis

(from Landau) For θ= θ ˙ =0, L 3 = L Z , E =0.  Near θ=0,  

V effective θ = L Z L 3 cosθ 2 2 I 1 sin 2 θ Mg 1cosθ L 3 2 1 2 θ 2 2 2 I 1 θ 2 1 2 Mg θ 2 = L 3 2 /8 I 1 1 2 Mg θ 2  The vertical position is stable against small oscillations provided   L 3 2 >4 I 1 Mg,  or Ω 3 2 >4 I 1 Mg/ I 3 2 .

Exercise: suppose you set the top vertical, but spinning at less than Ω 3 crit , the value at which it is just stable.  It will fall away, but bounce back, and so on. Show the maximum angle it reaches is given by cos θ/2 = Ω 3 / Ω 3 crit .