The
general equation governing rotation of a rigid body:
Having
dealt with situations where components of are
constant, we now ask what happens when is
also changed. For this let me look at the expression
for the angular momentum in the principal axis frame
again. It is
![1](22_3_clip_image004.gif)
I now give a slightly different derivation for the rate
of change of .
In doing this derivation I keep in mind that as a rigid
body rotates, the unit vectors along its principal axes
also rotate and their rate of change is (see previous
lecture)
![1](22_3_clip_image008.gif)
Now I differentiate to
get
![1](22_3_clip_image010.gif)
Here the first term is due to the change in the components
of along
the principal axis and the second term is the
change in due
to its rotation. Notice that we recover the formula
derived earlier if the components of do
not change with time, i.e. .
Let me repeat the interpretation of the equation:
at any instant we take the body rotating in the
principal axes frame at that time, i.e. the frame
is frozen at its position at that time and the
body is taken to be rotating in it. To see this
geometrically, let me take a two-dimensional
case. Shown in figure 10 are the principal axes
1 and 2 of a rigid body at times t and (t+ Δt) .
In time interval Δt the
body and the frame attached to it rotate by an
angle ,
and ω1 and ω2 change
to ω1 + Δω1 and ω2 + Δω2 .
With these changes let me calculate changes in
the components L1 and L2 in
the frame frozen at time t .
![1](22_3_clip_image015.gif)
Looking at the figure, where I have shown all the changes
that have taken place during the time interval Δt , we get in the frame at time t
![1](22_3_clip_image017.gif)
and
![1](22_3_clip_image019.gif)
So the total change in the angular momentum is
![1](22_3_clip_image021.gif)
Dividing both sides by Δt and taking proper
limit gives
![1](22_3_clip_image023.gif)
This gives you some idea about where this equation comes
from. Of course in a more accurate treatment, rotations
about the other axes also have to be taken into account.
For infinitesimal rotations, they can all be added up
and give the general equation
![1](22_3_clip_image025.gif)
This gives
![1](22_3_clip_image027.gif)
Each one of these rates of change should be equal to the
component of the torque in that direction .Thus
![1](22_3_clip_image029.gif)
These are the most general equations governing the dynamics
of a rigid body and are known as Euler's equations. I
now use it to explain the third experiment I had suggested
in the beginning of Lecture 21.
Example 4: Hold a rectangular box
at a height with one of its faces perpendicular to
the vertical, give it a spin and let it drop (see
figure 11). Describe its subsequent rotational motion.
![1](22_3_clip_image030.gif)
This is an example of torque-free ( )
motion because there is no torque on the box
about its centre of mass. Thus its rotational motion
is governed by the equations
![1](22_3_clip_image034.gif)
For a box similar to the one shown in figure 11 we
would generally have I3 > I2 > I1
.
Let me first consider the case when the box is given
a spin about its principal axis 1. Let me also assume
that in the process I also disturb it and give it
very small angular velocities ω2 and ω3 about
its axes 2 and 3, respectively. Since both ω2 and ω3 are very small, their product is second-order
in smallness and will be ignored. The Euler equations
and there are then as given below.
![1](22_3_clip_image036.gif)
The first equation implies that ω1 is
a constant. Let me call it the spin rate ω0 .
Using this fact the other two equations are dealt
with as follows. Differentiate equation (II) with
respect to time to get
![1](22_3_clip_image038.gif)
and substitute for from
equation (III) to obtain
![1](22_3_clip_image042.gif)
Since I3 > I2 > I1 , the equation above is of
the form
![1](22_3_clip_image044.gif)
Its solution is of the form
![1](22_3_clip_image046.gif)
One can similarly get equation for ω3 also
and see that it also has similar oscillatory
solution. This implies that as the box falls down
it spins about axis 1 and oscillates about axes 2
and 3. Since magnitudes of ω2 and ω3 are
small, you see the box fall essentially spinning
only. The same thing will happen if we give initial
spin about axis 3. However something different happens
when the initial spin is about axis 2. Assuming ω1 and ω3 to
be small, in this case the Euler equations take
the following form right after the release of the
box.
![1](22_3_clip_image048.gif)
The second equation above implies that ω2 is
a constant and with I3 > I2 > I1 , the
other two equations take the form
![1](22_3_clip_image050.gif)
Solution of these equations is of the form
![1](22_3_clip_image052.gif)
which indicates that right after the release, the angular
velocities about axes 1 and 3 will grow very fast and
take on a large value. Thus the box will start rotating
about all three axes and that is what you observe. Thus
we see that a rigid body is stable when it is given a
spin about the axes having the smallest or the largest
moment of inertia. However, if given a spin about the
axis with intermediate moment of inertia, it will be
unstable. Next I take up the case of precessing top that
I had not solved by employing Euler's equations earlier.
This is an example where a torque is also being applied
on the system
Example 5: Apply Euler's equations
to a precessing top and get its precession frequency Ω.
The top has a mass m and is spinning
at a rate of ωS (see
figure 12). Its centre of gravity is at a distance l from
the pivot point.
![1](22_3_clip_image054.gif)
I have already discussed about the principal axes
of the top in example 2 above. With the
Euler's equations for the top are
![1](22_3_clip_image058.gif)
Now in applying Euler's equations you have to keep
in mind that the top is spinning. As such its principle
axes 2 and 3 also rotate about axis 1 with angular
frequency ωS . So the components
of angular frequency and torque in the direction
of these axes also change with time. Taking time
at which the position of the top is shown in figure
12 to be t = 0, I draw in figure 13
the position of axes 2 and 3 at time t .
In this figure, I have neglected the angle W
t through which the top and therefore the torque
vector itself has rotated. In other words I have
assumed that .
Thus the angular velocity and torque are shown
where they were at t = 0 .
![1](22_3_clip_image063.gif)
Looking at figure 13, it is clear that the components
of the angular velocity and the torque are
![1](22_3_clip_image065.gif)
Substituting these in the Euler's equation for the top
gives
![1](22_3_clip_image067.gif)
The first of these equations gives ω1 =
constant = ωS.
The other two equations give the same answer
which is
![1](22_3_clip_image069.gif)
This is the answer that we have seen earlier. In solving
the Euler's equations for the top, we made the
assumption of .
Further we assumed that the top only precesses
about the vertical. However, there is no reason why
it cannot posses a horizontal angular velocity ΩH also.
Assuming the existence of Ω and ΩH and
then solving the Euler's equations will give
a more complete solution for the motion of a
spinning top. It in fact gives the nutating motion
also. You may want to try getting this general
solution.
With this lecture I end of the topic of rigid-body rotation.
|