Chapter 20 Fixed-axis rotation in 3D

Download this chapter by section

All chapters and individual parts

The ideas of chapter 10 on 2D circular motion are extended to fixed axis rotation in 3D. The key difference here is the non-trivial use of the cross product for calculating velocities and accelerations. Fixed axis rotation is the simplest motion with which one can introduce the full moment of inertia matrix, where the diagonal terms are analogous to the scalar 2D moment of inertia and the off-diagonal terms have a “centripetal” interpretation. The main new application is dynamic balance.

Filename:tfigure8-rel-ang-vel
Figure 20.1: A car crankshaft is a complex three-dimensional object which is well approximated for many purposes as rotating about a fixed axis. The relative timing of the reciprocating pistons is controlled by this complex shape. “Connecting rods” are pinned to the cylinders at one end and to the short offset cylinders on the crankshaft at the other.

In chapter 10 we discussed the motions of particles and rigid bodies that lie in a plane and rotate about a fixed axis perpendicular to that plane. In this chapter we again are going to think about fixed axis rotation, but now in three dimensions. The axis of rotation might be in any skew direction and the rotating bodies might be arbitrarily complicated three-dimensional shapes. As a cartoon, imagine a rigid body skewered with a rigid rod and then turned by a motor that speeds up and slows down. More practically think of the crankshaft in a car engine (fig. 20.1). Other applications include accelerating or decelerating shafts of all kinds, gears, turbines, flywheels, pendula, and swinging doors.

To understand this motion we need to take a little more care with the kinematics because it now involves three dimensions, although in some sense the basic ideas are unchanged from the previous two-dimensional chapter. The three-dimensional mechanics naturally gets more involved.

This one special motion, rotation about a fixed axis, serves as our introduction to three-dimensional rigid body mechanics.

As for all motions of all systems, the momentum balance equations apply to any system or any part of a system that has fixed axis rotation. So our mechanics results will be based on these familiar equations:

Linear momentum balance:𝑭i=𝑳˙,
Angular momentum balance:𝑴i/O=𝑯˙/O.

and

Power balance:P=E˙K.

As always, we will evaluate the left hand sides of the momentum equations using the forces and moments in the free-body diagram. We evaluate the right hand sides of these equations using our knowledge of the velocities and accelerations of the various mass points.

The chapter starts with a discussion of kinematics. Then we consider the mechanics of systems with fixed axis rotation. The moment of inertia matrix is then introduced followed by a section where the moment of inertia is used as a shortcut in the evaluation of 𝑯˙/O. Finally we discuss dynamic balance, an important genuinely three-dimensional topic in machine design.

20.1 3-D description of circular motion

Let’s first assume each particle is going in circles around the z axis, as in the previous chapter. The figure below shows this two-dimensional situation first two-dimensionally (left) and then as a two-dimensional motion in a three-dimensional world (right).

Filename:tfigure4-2-const-rate3D
Filename:tfigure8-ang-vel-ex
Figure 20.2: When a point P is going in circles about the z axis, we define the unit vector 𝒆ˆR to be pointed from the axis to the point P. We define the unit vector 𝒆ˆθ to be tangent to the circle at P. Both of these vectors change in time as the point moves along its circular path.
Filename:tfigure8-ang-vel
Figure 20.3: The acceleration is the sum of two components. One directed towards the center of the circle in the 𝒆ˆR direction, and one tangent to the circle in the 𝒆ˆθ direction.

Either way, the velocity and acceleration are the same:

  • the velocity is tangent to the circle it is going around and is proportional in magnitude to the radius of the circle and also to its angular speed. That is, the direction of the velocity is in the direction 𝒆ˆθ and has magnitude ωR, where ω is the angular rate of rotation and R is the radius of the circle that the particle is going around.

  • the acceleration can be constructed as the sum of two vectors. One is pointed to the center of the circle and proportional in magnitude to both the square of the angular speed and to the radius. The other vector is tangent to the circle and equal in magnitude to the rate of increase of speed.

These two ideas are summarized by the following formulas:

𝒗=ωR𝒆ˆθand𝒂=ω2Rv2/R𝒆ˆR+Rθ¨v˙𝒆ˆθ (20.1)

with

v=ωR. (20.2)

The axis of rotation might not be the z-axis of a convenient xyz coordinate system. So the xy plane of circles might not be the xy plane of the coordinate system you might want to use for some other reasons. Fortunately, we can write the formulas 20.1 in a way that rids us of these problems.

Here are some formulas which are equivalent to the formulas 20.1 but which do not make use of the polar coordinate base vectors.

𝒗 = 𝝎×𝒓, (20.3)
𝒂 = 𝝎×(𝝎×𝒓)+𝝎˙×𝒓. (20.4)

To check that equations 20.3 and 20.4 are really equivalent to 20.1 we need to verify that the vector 𝝎×𝒓 is equal to ωR𝒆ˆθ, that the vector 𝝎×(𝝎×𝒓) is equal to the vector ω2R𝒆ˆR, and that 𝝎˙×𝒓 is equivalent to Rθ¨𝒆ˆθ.

First, define 𝒓 as the position of the point of interest relative to any point on the axis of rotation. If this point happens to be the center of the circle then 𝒓=𝑹. But, in general, 𝒓𝑹.

Now look at 𝒗=𝝎×𝒓 with respect to fig. 20.2. Using the right hand rule, it is clear that the direction of the cross product 𝝎×𝒓 is in fact the 𝒆ˆθ direction. What about the magnitude? The magnitude |𝝎×𝒓|=|𝝎||𝒓|sinϕ. But |𝒓|sinϕ=R. So the magnitude of 𝝎×𝒓 is ωR. That is,

𝝎×𝒓 = (|𝝎×𝒓|)(unit vector in the direction of 𝝎×𝒓) (20.5)
= (ωR)𝒆ˆθ (20.6)
= 𝒗 (20.7)

So 𝒗=𝝎×𝒓 is correct. The check of the second term in the acceleration formula follows the same reasoning. But the check of the first term involves the triple cross product.

Triple cross product

The formula for acceleration of a point on a rigid body includes the centripetal term 𝝎×(𝝎×𝒓). This expression is a special case of the general vector expression

𝑨×(𝑩×𝑪)

which is sometimes called the ‘vector triple product’ because its value is a vector (as opposed to the scalar value of the ‘scalar triple product’). The primary useful identity with vector triple products (see (17.45) ) is:

𝑨×(𝑩×𝑪)=(𝑨𝑪)𝑩(𝑨𝑩)𝑪. (20.8)

This formula may be remembered by the semi-mnemonic device ‘cab minus bac’ since 𝑨𝑪=𝑪𝑨 and 𝑨𝑩=𝑩𝑨. This formula is discussed in box 17.2 on page 17.2.

So now we can write

𝝎×(𝝎×𝒓)=𝝎×(ωR)𝒆ˆθ=ω2R𝒆ˆR𝑹=ω2𝑹=𝒂. (20.9)

Note that equations 20.3 and 20.4 are vector equations. They do not make use of any coordinate system. So, for example, we can use them even if 𝝎 is not in the z direction.

Angular velocity of a rigid body in 3D

Filename:tfigure8-ang-accel
Figure 20.4: A rigid body spinning about the z axis. Every point on the body, like point P at 𝒓, is going in circles. All of these circles have centers on the axis of rotation. All the points are going around at the same angular rate, θ˙=ω.

If a rigid body is constrained to rotate about an axis then all points on the body have the same angular rate about that axis. Hence one says that the body has an angular velocity. So the measure of rotation rate of a three-dimensional rigid body is the body’s angular velocity vector 𝝎. At any instant in time a given body has one and only one angular velocity 𝝎. Although we only discuss fixed axis rotation in this chapter, a given body has a unique angular velocity for general motion.

For rotations about the z-axis, 𝝎=ω𝒌ˆ. ω is the θ˙ shown in fig. 20.4. Since all points of a rigid body have the same θ˙, even if they have different θ’s, the definition is not ambiguous.

We would like to make this idea precise enough to be useful for calculations. Why, one may ask, do we talk about rotation rate ω or 𝝎 instead of just using the derivative of an angle θ, namely θ˙? The answer is that for a rigid body one would have trouble deciding what angle θ to measure.

First recall the situation for a two-dimensional rigid body. Consider all possible θ1,θ2,θ3,, the angles that all possible lines marked on a body could make with the positive x-axis, the positive y axis, or any other fixed line that does not rotate. As the body rotates all of these angles increment by the same amount. Therefore, each of these angles increases at the same rate. Because all these angular rates are the same, one need not define θ˙1=ω1, θ˙2=ω2, θ˙3=ω3, etc. for each of the lines. Every line attached to the body rotates at the same rate and we call this rate ω. So θ˙1=ω, θ˙=ω, θ˙3=ω, etc. Rather than say the lengthy phrase ‘the rate of rotation of every line attached to the rigid body is ω’, we instead say ‘the rigid body has angular velocity ω’. For use in vector equations, we define the angular velocity vector of a two-dimensional rigid body as 𝝎=ω𝒌ˆ and for a 3-D body rotating about an axis in the 𝝀ˆ direction as 𝝎=ω𝝀ˆ.

Filename:tfigure8-ang-accel-ex
Figure 20.5: The shadows of lines marked in a 3-D rigid body are shown on a plane perpendicular to the axis of rotation. The shadows rotate on the plane at the rate θ˙1=θ˙2=ω. The angular velocity vector is 𝝎=ω𝝀ˆ.

What do we mean by these angles θi for crooked lines in a three-dimensional body? We simply look at shadows of lines drawn in or on the body of interest onto a plane perpendicular to the axis of rotation; i.e., perpendicular to 𝝀ˆ. See fig. 20.5. The rate of change of their orientation (θ˙1=θ˙2=θ˙3) is ω, and 𝝎 is therefore ω𝝀ˆ. This intuitive geometric definition of ω in terms of the rotation of shadows has run its course. It gives you a picture but is not very convenient for developing formulas.

Example: What are the velocity and acceleration of one corner of a cube that is spinning about a diagonal?

Filename:tfigure8-ang-accel1
Figure 20.6: A spinning cube.

A one foot cube is spinning at 60rpm about the diagonal OC. What are the velocity and acceleration of point B?

First let’s find the velocity using 𝒗=𝝎×𝒓:

𝒗 = 𝝎×𝒓
= (60rpm𝝀OC)×𝒓OB
= (2πs1(ıˆ+ȷˆ+𝒌ˆ)3)×(1ft(ȷˆ+𝒌ˆ))
= (2π/3)(ȷˆ+𝒌ˆ)ft/s.

Now of course this equation could have been worked out with the first of equations 20.1 but it would have been quite tricky to find the vectors 𝒆ˆθ, 𝑹, and 𝒆ˆR! More simply, to find the acceleration we just plug in the formula 𝒂=𝝎×(𝝎×𝒓) as follows:

𝒂 = 𝝎×[𝝎×𝒓]
= (60rpm𝝀OC)×[(60rpm𝝀OC)×𝒓OB]
= (2πs1(ıˆ+ȷˆ+𝒌ˆ)3)
×[(2πs1(ıˆ+ȷˆ+𝒌ˆ)3)×(1ft(ȷˆ+𝒌ˆ)]
= (2π/3)2(2ıˆȷˆ𝒌ˆ)ft/s2.

The last line of calculation is eased by the calculation of velocity above where the term in square brackets, the velocity, was already calculated.

Relative motion of points on a rigid body

The relative velocity of two points A and B is defined to be

𝒗B/A𝒗B𝒗A

So, the relative velocity of two points glued to one rigid body, as observed from a Newtonian frame, is given by

𝒗B/A 𝒗B𝒗A (20.10)
= 𝝎×𝒓B/O𝝎×𝒓A/O (20.11)
= 𝝎×(𝒓B/O𝒓A/O) (20.12)
= 𝝎×𝒓B/A, (20.13)

where point O is a point in the Newtonian frame on the fixed axis of rotation. Clearly, since points A and B are fixed in the body their velocities and hence their relative velocity as observed in a reference frame fixed to is 𝟎. But, point A has some absolute velocity that is different from the absolute velocity of point B, as viewed from point O in the fixed frame. The relative velocity of points A and B, the difference in absolute velocity of the two points, is due to the difference in their positions relative to point O.

Similarly, the relative acceleration of two points glued to one rigid body spinning about a fixed axis is

𝒂B/A𝒂B𝒂A=𝝎×(𝝎×𝒓B/A)+𝝎˙×𝒓B/A. (20.14)

Again, the relative acceleration is due to the difference in the points’ positions relative to the point O fixed on the axis. Like their junior 2D cousins, these kinematics results, 20.13 and 20.14, are useful for calculating angular momentum relative to the center-of-mass as well as for the understanding of the motions of machines with moving connected parts.

Filename:sfig8-2-3
Figure 20.7: Two points, A and B on one body that has a fixed axis of rotation.

To repeat, for two points on one rigid body we have that

𝒓˙=B/A𝝎×𝒓B/A. (20.15)

Equation (20.15) is perhaps the most fundamental equation for those desiring an understanding of the motions of rigid bodies. Unless one desires to pursue matrix representations of rotation, equation (20.15) is the defining equation for 𝝎. There is always exactly one vector 𝝎 so that equation (20.15) is true for every pair of points on a rigid body.

Equation (20.15) is not so simple a defining equation as one would hope for such an intuitive concept as spinning. But, besides the pictorial definition with shadows, it’s the simplest definition we have.

Relative velocity and acceleration using rotating frames

If we glued a coordinate system xy to a rotating rigid body 𝒞, we would have what is called a rotating frame as shown in fig. 20.8.

Filename:sfig8-2-3a
Figure 20.8: A rotating rigid body 𝒞 with rotating frame xy attached.

The base vectors in this frame change in time the same way as did 𝒆ˆR and 𝒆ˆθ in section 15.1. That is

ddtıˆ=𝝎𝒞×ıˆandddtȷˆ=𝝎𝒞×ȷˆ.

If we now write the relative position of B to A in terms of ıˆ and ȷˆ, we have

𝒓B/A=xıˆ+yȷˆ.

Since the coordinates x and y rotate with the body to which A and B are attached, they are constant with respect to that body,

x˙=0 and y˙=0.

So

ddt(𝒓B/A) = ddt(xıˆ+yȷˆ)
= x˙0ıˆ+xddtıˆ+y˙0ȷˆ+yddtȷˆ
= x(𝝎𝒞×ıˆ)+y(𝝎𝒞×ȷˆ)
= 𝝎𝒞×(xıˆ+yȷˆ)𝒓B/A
= 𝝎𝒞×𝒓B/A.

If we now try to calculate the rate of change of 𝒗B/A,

ddt(𝒗B/A) = ddt(𝝎𝒞×𝒓B/A)
= d𝝎𝒞dt×𝒓B/A+𝝎𝒞×d𝒓B/Adt
𝒂B/A = 𝝎˙𝒞×𝒓B/A+𝝎𝒞×(𝝎𝒞×𝒓B/A).

Mechanics

Now that we know the velocity and acceleration of every point in the system we are ready, in principle, to find 𝑳˙ and 𝑯˙/O in terms of the angular velocity vector 𝝎, its rate of change 𝝎˙, and the position of all the mass in the system. This we do in the next section.

SAMPLE 20.1  Kinematics in 3-D—some basic questions: The following questions are about the velocity and acceleration formulae for the non-constant rate circular motion about a fixed axis:

𝒗 = 𝝎×𝒓
𝒂 = 𝝎˙×𝒓+𝝎×(𝝎×𝒓)
  1. 1.

    In the formulae above, what is r ? How is it different from 𝑹=R𝒆ˆR used in the formulae 𝒗=Rθ˙𝒆ˆθ?

  2. 2.

    What is the difference between θ˙ and 𝝎, and θ¨ and 𝝎˙?

  3. 3.

    Are the parentheses around the 𝝎×𝒓 term necessary in the acceleration formula?

  4. 4.

    Under what condition(s) can a particle have only tangential acceleration?

Solution

  1. 1.

    In the formulae for velocity and acceleration, 𝒓 refers to a vector from any point on the axis of rotation to the point of interest. Usually the origin of a coordinate system located on the axis of rotation is a convenient point to take as the base point for 𝒓. You can, however, choose any other point on the axis of rotation as the base point.

    The vector 𝒓 is different from 𝑹 in that 𝑹 is the position vector of the point of interest with respect to the center of the circular path that the point traces during its motion. See Fig. 5.2 of the text.

  2. 2.

    θ˙ and θ¨ are the magnitudes of angular velocity and angular acceleration, respectively, in planar motion, i.e., 𝛚=θ˙𝐤ˆ and 𝝎˙=θ¨𝒌ˆ. We have introduced these notations to highlight the simple nature of planar circular motion. Of course, you are free to use 𝝎=ω𝒌ˆ and 𝝎˙=α𝒌ˆ if you wish.

  3. 3.

    Yes, the parentheses around 𝝎×𝒓 in the acceleration formula are mandatory. The parentheses imply that this term has to be calculated before carrying out the cross product with 𝝎 in the formula. Since the term in the parentheses is the velocity, you may also write the acceleration formula as

    𝒂=𝝎˙×𝒓+𝝎×𝒗.

    Even if the formula is clear in your mind and you know which cross product to carry out first, it is a good idea to put the parentheses.

  4. 4.

    First of all let us identify the tangential and the normal (or radial) components of the acceleration:

    𝒂=𝝎˙×𝒓tangential+𝝎×(𝝎×𝒓).radial/normal

    Clearly, for a particle to have only tangential acceleration, the second term must be zero. For the second term to be zero we must have either 𝒓=𝟎 or 𝝎=𝟎. But if 𝒓=𝟎, then the tangential acceleration also becomes zero; the particle is on the axis of rotation and hence has no acceleration. Thus the condition that allows only tangential acceleration to survive is 𝝎=𝟎. Now remember that 𝝎˙ is not zero. Therefore, the condition we have found can be true only momentarily. This disappearance of the radial acceleration happens at start-up motions and in direction-reversing motions. margin: In all start-up motions, the velocity is zero but the acceleration is not zero at the start up (t=0). In direction-reversing motions, such as that of the washing machine drum during the wash-cycle, just at the moment when the direction of motion reverses, velocity becomes zero but the acceleration is non-zero.

SAMPLE 20.2  For a particle in circular motion, 𝝎=2rad/s𝒌ˆ,𝜶=4rad/s2𝒌ˆ and 𝒓G=3m(cos30ıˆ+sin30ȷˆ), where 𝒓G is the position vector of the particle. Find (i) 𝒗, (ii) the tangential acceleration 𝒂t, and (iii) the radial acceleration 𝒂r, and show the resulting vectors.

Solution

Filename:sfig8-2-3b
Figure 20.9: 𝒗=𝝎×𝒓G.
𝝎 = ω𝒌ˆ=2rad/s𝒌ˆ
𝜶 = α𝒌ˆ=4rad/s2𝒌ˆ
𝒓G = rGxıˆ+rGyȷˆ=3m(cos30ıˆ+sin30ȷˆ).

(i) From the given formulae, the linear velocity

𝒗 = 𝝎×𝒓G
= ω𝒌ˆ×(rGxıˆ+rGyȷˆ)=ωrGxȷˆ+ωrGy(ıˆ)
= 6m/s(cos30ȷˆsin30ıˆ)
= 3m/s(ıˆ+3ȷˆ).

The velocity vector 𝒗 is perpendicular to both 𝝎 and 𝒓G. These vectors are shown in Fig. 20.9. You should use the right hand rule to confirm the direction of 𝒗.

Answer: v=3m/s(ıˆ+3ȷˆ)

(ii) The tangential acceleration

Filename:sfig8-3-1
Figure 20.10: 𝒂t=𝜶×𝒓G.
𝒂t = 𝜶×𝒓G
= α𝒌ˆ×(rGxıˆ+rGyȷˆ)=αrGxȷˆαrGyıˆ
= 8m/s2(cos30ȷˆsin30ıˆ).

Since 𝝎 and 𝜶 are in the same direction, calculation of 𝒂t is similar to that of 𝒗 and 𝒂t has to be in the same direction as 𝒗. This vector is shown in Fig. 20.10. Once again, just as in the case of 𝒗 we could easily check that 𝒂t is perpendicular to both 𝜶 and 𝒓G.

Answer: at=4m/s2(3m/s(ıˆ+3ȷˆ)

(iii) Finally, the radial acceleration

Filename:sfig8-3-1a
Figure 20.11: 𝒂r=𝝎×(𝝎×𝒓G)
𝒂r = 𝝎×(𝝎×𝒓G)
= ω𝒌ˆ×(ωrGxȷˆωrGyıˆ)
= ω2rGx(ıˆ)ω2rGyȷˆ=ω2𝒓G
= 12m/s2(cos30ıˆsin30ȷˆ).

This cross product is illustrated in Fig. 20.11. Both from the illustration as well as the calculation you should be able to see that 𝒂r is in the direction of 𝒓G. In fact, you could show that 𝒂r=ω2𝒓G.

Answer: ar=6m/s2(3ıˆȷˆ)

SAMPLE 20.3

Filename:sfig8-5-wiper
Figure 20.12:
Filename:sfig8-5-wiper-a
Figure 20.13:

Simple 3-D circular motion. A system with two point masses A and B is mounted on a rod OC which makes an angle θ=45 with the horizontal. The entire assembly rotates about the y-axis with constant angular speed ω=3rad/s, maintaining the angle θ. Find the velocity of point A. What is the radius of the circular path that A describes? Assume that at the instant shown, AB is in the xy plane.

Solution The angular velocity of the system is

𝝎=ωȷˆ=3rad/sȷˆ.

Let 𝒓A be the position vector of point A. Then the velocity of point A is

𝒗 = 𝝎×𝒓A( and 𝒓A=𝒓C+𝒓A/C)
= ωȷˆ×(lcosθıˆ+lsinθȷˆ𝒓C+dcosθȷˆdsinθıˆ𝒓A/C)
= ωȷˆ×[(lcosθdsinθ)ıˆ+(lsinθ+dcosθ)ȷˆ]
= (ωlcosθωdsinθ)𝒌ˆ
= [3rad/s(1mcos450.5msin45)]
= 1.06m/s𝒌ˆ

Answer: v=1.06m/skˆ.

Filename:sfig8-5-2disks
Figure 20.14: Circular trajectory of point A as seen by looking down along the y-axis. At the instant shown, 𝒆ˆθ=𝒌ˆ.

We can find the radius of the circular path of A by geometry. However, we know that the velocity of A is also given by

𝒗=ωR𝒆ˆθ

where R is the radius of the circular path. At the instant of interest, 𝒆ˆθ=𝒌ˆ (see figure 20.14).

Thus𝒗=ωR𝒌ˆ.

Comparing with the answer obtained above, we get

1.06m/s𝒌ˆ = ωR𝒌ˆ
 R = 1.06m/s3rad/s
= 0.35m.

Answer: R=0.35m

SAMPLE 20.4

Filename:sfig8-4-1
Figure 20.15:

Velocity and acceleration in 3-D: The rod shown in the figure rotates about the y-axis at angular speed 10rad/s and accelerates at the rate of 2rad/s2. The dimensions of the rod are L=h=2m and r=1m. There is a small mass P glued to the rod at its free end. At the instant shown, the three segments of the rod are parallel to the three axes.

  1. 1.

    Find the velocity of point P at the instant shown.

  2. 2.

    Find the acceleration of point P at the instant shown.

Solution We are given:

𝝎=ωȷˆ=10rad/sȷˆ and 𝝎˙=ω˙ȷˆ=2rad/s2ȷˆ.
  1. 1.

    The velocity of point P is

    𝒗=𝝎×𝒓.

    At the instant shown, the position vector of point P (the vector 𝒓P/O) seems to be a good choice for 𝒓. margin: An even better choice, perhaps, is the vector 𝒓P/A. Remember, the only requirement on 𝒓 is that it must start at some point on the axis of rotation and must end at the point of interest. Thus,

    𝒓𝒓P/O=Lıˆ+hȷˆ+r𝒌ˆ.

    Therefore,

    Filename:sfig8-4-1a
    Figure 20.16: By drawing the velocity vector at point P (a vector tangent to the path) we see that 𝒗P must have a positive x-component and a negative z-component.
    𝒗 = ωȷˆ×(Lıˆ+hȷˆ+r𝒌ˆ)
    = ω(L𝒌ˆ+rıˆ)
    = 10rad/s(2m𝒌ˆ+1mıˆ)
    = (20ıˆ10𝒌ˆ)m/s.

    As a check, we look down the y-axis and draw a velocity vector at point P (tangent to the circular path at point P) without paying attention to the answer we got. From the top view in Fig. 20.16 we see that at least the signs of the components of 𝒗 seem to be correct.

    Answer: 𝒗=(20ıˆ10𝒌ˆ)m/s

  2. 2.

    The acceleration of point P is

    𝒂 = 𝝎˙×𝒓+𝝎×(𝝎×𝒓)
    = ω˙ȷˆ×(Lıˆ+hȷˆ+r𝒌ˆ)+ωȷˆ×ω(L𝒌ˆ+rıˆ)𝝎×𝒓
    = ω˙(L𝒌ˆ+rıˆ)+ω2(Lıˆr𝒌ˆ)
    = 2rad/s2(2m𝒌ˆ+1mıˆ)100(rad/s)2(2mıˆ+1m𝒌ˆ)
    = (98ıˆ+104𝒌ˆ)m/s2.

    We can check the sign of the components of 𝒂 also. Note that the tangential acceleration, 𝝎˙×𝒓, is much smaller than the centripetal acceleration, 𝝎×(𝝎×𝒓) . Therefore, the total acceleration is almost in the same direction as the centripetal acceleration, that is, directed from point P to A. If you draw a vector from P to A, you should be able to see that it has negative components along both the x- and z-axes. Thus the answer we have got seems to be correct, at least in direction.

    Answer: 𝒂=(98ıˆ+104𝒌ˆ)m/s2

Problems for 20.1 3-D description of circular motion

20.1.1  If 𝑨 and 𝑩 are perpendicular, what is |𝑨×(𝑨×𝑩)| in terms of |𝑨| and |𝑩|? (Hint: make a neat 3-D sketch.)

20.1.2  A circular disk of radius r=100mm rotates at constant speed ω about a fixed but unknown axis. At the instant when the disk is in the xy-plane, the angular velocity of each point on the diameter AB (which is parallel to the x-axis) is 𝒗=2m/s𝒌ˆ. At the same instant, the magnitude of the acceleration of the center of mass C is 50m/s2.

  1. (a)

    Find the angular velocity of the disk.

  2. (b)

    Find the location of the axis of rotation.

  3. (c)

    Find the acceleration of points A and B.

Filename:sfig8-6-2
Figure 20.17

20.1.3  In the expression for normal acceleration 𝒂n=𝝎×(𝝎×𝒓), the parenthesis around 𝝎×𝒓 is important because vector product is nonassociative. Showing 𝝎×(𝝎×𝒓)(𝝎×𝝎)×𝒓 is almost trivial. (why?) Take 𝝎1=ω1𝒌ˆ, 𝝎2=ω2ıˆ+ω1𝒌ˆ and 𝒓=aıˆbȷˆ and show that 𝝎1×(𝝎2×𝒓)(𝝎1×𝝎2)×𝒓.

20.1.4

For the given angular velocity 𝝎 and the position vector 𝒓 below, Compute the velocity 𝒗 in each of the three cases. Draw the circular orbit, and show the vectors 𝝎,𝒓, and 𝒗.

  1. (a)

    𝝎=2rad/s𝒌ˆ,𝒓=3mıˆ+4mȷˆ.

  2. (b)

    𝝎=4.33rad/sıˆ+2.5rad/sȷˆ,𝒓=0.5mıˆ0.87mȷˆ.

  3. (c)

    𝝎=2rad/sıˆ+2rad/sȷˆ+2.83rad/s𝒌ˆ,𝒓=3.5mıˆ+3.5mȷˆ4.95m𝒌ˆ.

20.1.5  A shell moves down a barrel. At the instant of interest the barrel is being raised at a constant angular speed ω1. Simultaneously the gun turret (the structure which holds the barrel) rotates about the vertical axis at a constant rate ω2. At the instant of interest, the barrel is tipped up an angle θ and the mass m is moving at speed v along the frictionless barrel. Neglecting gravity, find the force of the barrel on the shell.

Filename:sfig8-6-2a
Figure 20.18

A shell moves along a cannon barrel.

20.1.6  The following questions are about the velocity and acceleration formulae for the non-constant rate circular motion about a fixed axis:

𝒗 = 𝝎×𝒓
𝒂 = 𝝎˙×𝒓+𝝎×(𝝎×𝒓)
  1. (a)

    For planar circular motion about the z-axis, take

    𝝎˙=α𝒌ˆ,𝝎=ω𝒌ˆ, and 𝒓=R𝒆ˆR

    and show that the velocity and acceleration formulae reduce to the formulae from elementary physics: 𝒗=ωR𝒆ˆθ and 𝒂=ω2R𝒆ˆR+Rα𝒆ˆθ.

  2. (b)

    Draw a schematic picture of a particle (or a point on a rigid body) going in circles about a fixed axis in 3-D. Take two distinct points on the axis and draw vectors 𝒓1 and 𝒓2 from these two points to the particle of interest. Show, using the formulae for 𝒗 and 𝒂 above, that the acceleration and velocity of the particle are the same irrespective of which 𝒓 is used and thus conclude that the vector 𝒓 can be taken from any point on the axis of rotation.

20.1.7  A rectangular plate with width w=1m and length =2m is rigidly attached to a shaft along one of its edges. The shaft rotates about the x-axis. The center of mass of the plate goes around at a constant linear speed v=2.5m/s.

  1. (a)

    Find the angular velocity of the plate.

  2. (b)

    At the instant shown, find the velocity of point P.

Filename:efig1-2-28
Figure 20.19

20.1.8  The solid cylinder shown in the figure has an angular velocity 𝝎 of magnitude 40 rad/s. The vector 𝝎 lies in the yz-plane. The origin O of the coordinate frame is at the center of the cylinder and is at rest. Find the following quantities:

  1. (a)

    The vector 𝝎.

  2. (b)

    The velocity 𝒗A/B of point A relative to point B.

  3. (c)

    The acceleration 𝒂A of point A.

Filename:pfigure-blue-118-2
Figure 20.20

20.1.9  A rigid body rotates about a fixed axis with an angular acceleration 𝜶=3π2(ıˆ+1.5ȷˆ+1.5𝒌ˆ)rad/s2. At a certain instant, the x and y components of the velocity of a point A of the body are vx=3m/s and vy=4m/s, respectively. Find the magnitude v of the velocity of point A.

20.1.10  Particle attached to a shaft, no gravity. A shaft of negligible mass spins at a constant rate. The motor at A imposes this rate if needed. The bearing at A prevents translation of its end of the shaft in any direction but causes no torques other than that of the motor (which causes a torque in the 𝒌ˆ direction). The frictionless bearing at C holds that end of the shaft from moving in the ıˆ and ȷˆ directions but allows slip in the 𝒌ˆ direction. A bar of negligible mass is welded at right angles to the shaft. A small sphere, considered as a particle, of mass 1kg is attached to the free end of the bar. There is no gravity. The shaft is spinning counterclockwise (when viewed from the positive 𝒌ˆ direction) at 3rad/s. In the configuration of interest the rod BD is parallel to the ıˆ direction.

  1. (a)

    What is the position vector of the mass in the given coordinate system at the instant of interest?

  2. (b)

    What is angular velocity vector of the shaft?

  3. (c)

    What are the velocity and acceleration of the particle at the instant of interest?

Filename:pfigure-s95f3a
Figure 20.21

20.1.11  Consider the particle on a spinning shaft of problem 20.1 again.

  1. (a)

    What is the velocity of a point halfway down the rod BD relative to the velocity of the mass at D?

  2. (b)

    What is the acceleration of a point halfway down the rod BD relative to the acceleration of the mass at D?

20.1.12  Particle attached to a spinning and accelerating shaft, no gravity. Consider the particle on a spinning shaft of problem 20.1 again. At the instant of interest, the shaft is spinning counterclockwise (when viewed from the positive 𝒌ˆ direction) at ω=3rad/s and ω˙=5rad/s2. In the configuration of interest the rod BD is parallel to the ıˆ direction. In the coordinate system shown,

  1. (a)

    What is the position vector of the particle at the instant of interest?

  2. (b)

    What are angular velocity and angular acceleration vectors of the shaft?

  3. (c)

    What are the velocity and acceleration of the particle at the instant of interest?

Filename:Danef94s1q2
Figure 20.22

20.1.13  Particle attached to a spinning and accelerating shaft, still no gravity, relative motion. Reconsider the particle on the spinning and accelerating shaft of problem 20.1. What is the velocity and acceleration of a point halfway down the rod BD relative the velocity and acceleration of the particle?

20.1.14  A crooked rod spinning at a constant rate. A one meter long uniform rod is welded at its center to a shaft at an angle ϕ as shown. The shaft is supported by bearings at its ends and spins at a constant rate ω. At the instant of interest, the crooked rod lies in the xy-plane. In the coordinate system shown,

  1. (a)

    What is the position vector of the point P on the rod at the instant of interest? Answer: rP/O=0.5m[(cosϕ)ıˆ+sinϕȷˆ].

  2. (b)

    What are the velocity and acceleration of point P at the instant of interest?

    Answer: vP=(0.5ωsinϕkˆ)m/s, 𝒂P=(0.5ω2sinϕȷˆ)m/s2.

Filename:pfigure-blue-123-1
Figure 20.23

A crooked rod spins on a shaft.

20.1.15  The spinning crooked rod, again. Reconsider the crooked rod on a shaft in problem 20.23. What are the velocity and acceleration of point P relative to point Q?

Filename:pfigure-blue-119-2
Figure 20.24

A crooked rod spins on a shaft.

20.1.16  A crooked rod spinning at a variable rate. Consider the crooked rod shown in problem 20.23 again. The shaft is spinning counterclockwise (when viewed from the positive x direction) at a constant angular acceleration ω˙=5rad/s2. At the instant of interest, the crooked rod lies in the xy-plane and its angular speed is ω=10rad/s. In the coordinate system shown,

  1. (a)

    Find the position vector of the point P on the rod?

  2. (b)

    What are angular velocity and angular acceleration vectors of the rod and shaft?

  3. (c)

    What are the velocity and acceleration of point P?

  4. (d)

    What is the velocity and acceleration of point P relative to point Q?

20.1.17  A rigid rod OP of length =500mm is welded to a shaft AB at an angle θ=30. The shaft rotates about its longitudinal axis at a constant rate ω=120rpm. At the instant shown, the rod OP is in the the xz-plane.

  1. (a)

    Neatly draw the path of point P. At the instant shown, draw the basis vectors 𝒆ˆr and 𝒆ˆθ and express them in terms of the basis vectors ıˆ, ȷˆ, 𝒌ˆ.

  2. (b)

    Find the radius of the circular path of point P and calculate the velocity and acceleration of the point using planar circular motion formulae 𝒗=θ˙𝒆ˆθ and 𝒂=θ˙2R𝒆ˆr.

  3. (c)

    Find the position vector 𝒓𝒓P/O of point P. Use this position vector to compute the velocity and acceleration of point P using the general formulae 𝒗=𝝎×𝒓 and 𝒂=𝝎×𝝎×𝒓. Show that the answers obtained here are the same as those in (b).

Filename:pfigure-s95q14
Figure 20.25

20.1.18  A uniform rectangular plate ABCD of width and length 2 rotates about its diagonal AC at a constant angular speed ω=5rad/s. At the instant shown, the plate is in the xz-plane.

  1. (a)

    Find the velocity and acceleration of point B at the instant shown.

  2. (b)

    Find the velocity and acceleration of point B when the plate is in the yz-plane.

Filename:Danef94s3q2
Figure 20.26

20.1.19  A rectangular plate of width w=100mm and length =200mm rotates about one of its diagonals at a constant angular rate ω=3rad/s. There is a fixed coordinate system xyz with the origin at the center of the plate and the x-axis aligned with the axis of rotation. Another set of axes xyz, aligned with the principle axes of the plate, has the same origin but is attached to the plate and therefore, rotates with the plate. At the instant of interest, the plate is in the xy-plane and the two axes z and z coincide. Find the velocity of the corner point P as follows:

  1. (a)

    Write the angular velocity 𝝎 of the plate in terms of its components in the two coordinate systems.

  2. (b)

    Find the velocity of point P in both coordinate systems using the expressions for 𝝎 found in part (a).

  3. (c)

    Find the perpendicular distance d of point P from the axis of rotation, and show that the linear speed of point P calculated from velocities found in part (b) is ωd in each case.

Filename:bikefork1-ang-accel
Figure 20.27

20.1.20  A cone of angle β sits on its side on a horizontal disk of radius r=100mm as shown in the figure. The cone is attached to shaft AB that turns in its bearing. The disk rotates at a constant angular speed ω=120rpm. Assume that the cone rolls without slip on the disk. Thus, the disk drives the cone-and-shaft assembly causing it to turn in the bearing. Find the angular velocity of the cone as a function of the cone angle β. In particular, find the angular velocity of the cone for

  1. (a)

    β=30, and

  2. (b)

    β=90 (or 90 minus a hair if that makes you more comfortable).

Filename:bikefork-ang-accel
Figure 20.28

20.2 Dynamics of fixed-axis rotation

We now address mechanics questions concerning objects which are known a priori to spin about a fixed axis. We would like to calculate forces and moments if the motion is known. And we would like to determine the details of the motion, the angular acceleration in particular, if the applied forces and moments are known. Once the angular acceleration is known (as a function of some combination of time, angle and angular rate) the angular rate and angular position can be found by integration or solution of an ordinary differential equation.

The full content of the subject follows from the basic mechanics equations

linear momentum balance,𝑭i=𝑳˙
angular momentum balance,𝑴i/C=𝑯˙/C,
and power balance:P=E˙K+E˙P+E˙int.

The quantities 𝑳˙ and 𝑯˙/C are defined in terms of the position and acceleration of the system’s mass (see the second page of the inside cover). To evaluate 𝑳˙ and 𝑯˙/C for fixed-axis rotation we can use the kinematics relations from the previous chapter which determine velocity and acceleration of points on a body spinning about a fixed axis in terms of the position 𝒓 of the point of interest relative to any point on the axis.

𝒗 = 𝝎×𝒓,
𝒂 = 𝝎×(𝝎×𝒓)+𝝎˙×𝒓.

For fixed-axis rotation 𝝎=ω𝝀ˆ and 𝝎˙=ω˙𝝀ˆ with 𝝀ˆ a constant unit vector along the axis of rotation.

To solve problems we draw a free-body diagram, write the equations of linear and angular momentum balance, and evaluate the terms using the kinematics relations. In general this will lead to the evaluation of a sum or an integral. A short cut, the moment of inertia matrix, will be introduced in later sections.

Before proceeding to more difficult three-dimensional problems, let’s review a simple 2D problem.

Example: Spinning disk

The round flat uniform disk in fig. 20.29 is in the xy plane spinning at the constant rate 𝝎=ω𝒌ˆ about its center. It has mass mtot and radius R0. What force is required to cause this motion? What torque? What power?

From linear momentum balance we have:

𝑭i=𝑳˙=mtot𝒂cm=𝟎,

which we could also have calculated by evaluating the integral 𝑳˙𝒂𝑑m instead of using the general result that 𝑳˙=mtot𝒂cm. From angular momentum balance we have:

𝑴i/O = 𝑯˙/O
 𝑴 = 𝒓/O×𝒂𝑑m
= 0R002π(R𝒆ˆR)×(Rω2𝒆ˆR)mtotπRO2RdθdRdAdm
= 𝟎𝑑θ𝑑R
= 𝟎.
Filename:bikefork1-alt
Figure 20.29: A uniform disk turned by a motor at a constant rate.

Power balance is of limited use for constant rate circular motion. If all parts of a system move at constant angular rate at a constant radius then they all have constant speed. Thus the kinetic energy of the system is constant. So the power balance equation just says that the net power into the system is the amount dissipated inside (assuming no energy storage).

Example: Spinning disk with power balance

Consider the spinning disk from fig. 20.29 and the previous example. The power balance equation III gives

P=E˙K+E˙P0+E˙int0  P=𝒗𝒂𝑑m=0𝑑m= 0. (20.16)

In this example there is no force or torque acting on the disk so the power P must turn out to be zero. In other constant rate problems the force and moment will not turn out to be zero, but the kinetic energy of the system will still be constant and so, assuming no energy storage or dissipation, we will still have P=0.

A stick sweeps out a cone

Now we consider a genuinely three-dimensional problem involving fixed-axis, rigid body rotation. Consider a long narrow stick swinging in circles so that it sweeps out a cone (fig. 20.30). Each point on the stick is moving in circles around the z-axis at a constant rate ω. What is the relation between ω and the angle of the stick ϕ? The approach to this problem is, as usual, to draw a free-body diagram, write momentum balance equations, evaluate the left and right hand sides, and then solve for quantities of interest. The hard part of this problem is evaluating the right hand side of the angular momentum balance equations.

Filename:bikefork-alt
Figure 20.30: A spherical rigid body pendulum (uniform thin rod) going in circles at constant rate 𝝎.
Filename:pfigure-blue-128-1
Figure 20.31: Free-body diagram of the rod.

To simplify calculation, we look at the pendulum at the instant it passes through the yz-plane, assuming the xyz axes are fixed in space.

The free-body diagram shown in fig. 20.31 shows the gravity force at the center-of-mass, the reaction force at point O, and, consistent with the shown construction of the hinge, the moments at O perpendicular to the hinge.

Because we are interested in the relation between ϕ and ω and not the reaction force, at least for now, we look at angular momentum balance about point O.

𝑴O=𝑯˙/O

First, we show and discuss the results of evaluating the equation of angular momentum balance. Then, we will show the details of calculating 𝑴O and the details of several methods for calculating 𝑯˙/O.

Evaluation of 𝑴O

To find 𝑴O, we had to find the moment of the gravity force. The most direct method is to use the definition

𝑴O = 𝒓/O×𝑭
= 2[cosϕ(𝒌ˆ)+sinϕȷˆ]×(mg𝒌ˆ)
= mgsinϕ2ıˆ
Filename:tfigure8-syst-bods
Figure 20.32: The rod is shown on the yz plane. We use this figure to locate the bit of mass dm corresponding to the bit of length of rod ds.

One could also ‘slide’ the gravity force to the level of O (a force displaced in its direction of action is mechanically equivalent). Then you can see from the figure that the force is perpendicular to its position relative to O. Moment is then force (mg) times distance (2sinϕ) in the direction given by the right hand rule (ıˆ).

Evaluation of 𝑯˙/O

We now evaluate 𝑯˙/O by adding up the contribution to the sum from each bit of mass.

𝑯˙/O = 𝒓/O𝒓/O=scosϕ𝒌ˆ+ssinϕȷˆ×𝒂For constant rate circular motion, 
a = × ω ( × ω r / O )
= × ( ω ˆ k ) [ ( ω ˆ k ) ×
s ( - sin ϕ ˆ ȷ cos ϕ ˆ k ) ]
= - ω 2 s sin ϕ ˆ ȷ
 
dmdm=ρds, where  ρ = mass per unit length
= 0(scosϕ𝒌ˆ+ssinϕȷˆ)𝒓/O×(ω2ssinϕȷˆ)𝒂(ρds)
= 0s2cosϕsinϕω2ıˆρds (evaluating the cross product)
= cosϕsinϕω2ρıˆ0s2𝑑s (ϕ,ρ, and ω do not vary with s)
= cosϕsinϕω2(ρ33)ıˆ (evaluating the integral)
𝑯˙/O = cosϕsinϕω2(m23)ıˆ (because m=ρ).

We could have taken a short-cut in the calculation of acceleration 𝒂. Instead of using 𝒂=𝝎×(𝝎×𝒓), we could have used 𝒂=ω2𝑹 where 𝑹 is the radius of the circle each particle is traveling on. It is evident from the picture that the appropriate radius is 𝑹=ssinϕȷˆ, so 𝒂=ω2ssinϕȷˆ.

We will show two more methods for calculating 𝑯˙/O in section 20.4 on page 20.4 using the moment of inertia matrix.

The results for the conically swinging stick

We can now evaluate the terms in the angular momentum balance equation as

mgsinϕ2ıˆ+MOyȷˆ+MOz𝒌ˆ𝑴O=sinϕcosϕm23ω2ıˆ𝑯˙/O. (20.18)

We can get three scalar equations from eqn. 20.18 by dotting it with ȷˆ, 𝒌ˆ, and ıˆ to get

MOy=0 and MOz=0

and

ω2=3g2cosϕ.

Note that MOy=MOz=0. That is, for this special motion, the hinge joint at O could be replaced with a ball-and-socket joint.

Note that the solution for a point mass spherical pendulum is ω2=gcosϕ. That is, this stick would rotate at the same rate and angle as a point mass at the end of a rod of length 23. One could not easily anticipate this result. We point it out here to emphasize that the analysis of this rigid-body problem cannot be reduced a priori to any simple particle mechanics problem.

Filename:sfig8-7-2
Figure 20.33: Plot of non-dimensional rotational speed ω2g versus hang angle ϕ. For ω2/g<3/2 the only solution is ϕ=0 (hanging straight down). At or very close to ω2/g=3/2 a range of ϕ’s is possible. As ϕπ/2 and the rod becomes close to horizontal, the spin rate ω goes to infinity.

In fig. 20.33, non-dimensional rotational speed ω2g is plotted versus hang angle ϕ. As one might expect intuitively unless ω is high enough, (ω2>3g2), the only solution is hanging straight down (ϕ=0).

At the critical speed (ω2=3g2), the curve is nearly flat, implying that a range of hang angles ϕ is possible all with nearly the same angular velocity. As is also intuitively plausible, as the bar gets close to the horizontal (close to π2), the spin rate goes to infinity.

The scalar equations governing rotation about an axis

For two-dimensional motion of flat hinged objects we had the simple relation “M=Iα”. This formula captures our simple intuitions about angular momentum balance. When you apply torque to a body its rate of rotation increases. It turns out that, for three-dimensional motion of a rigid body about a fixed axis, the same result applies if we interpret the terms correctly.margin: Caution: For more general three-dimensional motion than rotation about a fixed axis the equation M=Iα does not apply. Trying to vectorize by underlining various terms gives the wrong answer.

If the axis of rotation goes through C and is in the direction 𝝀ˆ we can define M=𝝀ˆ𝑴C as the moment about the axis of rotation. We can similarly look at the 𝝀ˆ component of 𝑯˙/C (assume, for definiteness, that the system is continuous).

𝝀ˆ𝑯˙/C = 𝝀ˆ𝒓×𝒂𝑑m (20.19)
= 𝝀ˆ𝒓×[(ω𝝀ˆ)×((ω𝝀ˆ)×𝒓)+(ω˙𝝀ˆ)×𝒓]𝑑m.
= ω˙R2𝑑m,

where R is the distance of the mass points from the axis. The last line follows from the previous most simply by paying attention to directions and magnitudes when using the right-hand rule and the geometric definition of the cross product. We thus have derived the result that

M=Iα,

if by M we mean moment about the fixed axis and by I we mean R2𝑑m. Actually, the scalar we call I in the above equation is a manifestation of a more general matrix [𝑰] that we will explore in the next section.

SAMPLE 20.5

Filename:sfig8-7-2a
Figure 20.34: A carnival ride rotating at a constant speed

Going on a carnival ride at a constant rate. A carnival ride with roof AB and carriage BC is rotating about the vertical axis with constant angular velocity 𝝎=ωȷˆ. If the carriage with its occupants has mass m=100kg, find the tension in the inextensible and massless rod BC when θ=30o. What is the required angular speed ω (in revolutions/minute) to maintain this angle?

margin:

Solution

The free-body diagram of the carriage is shown in Fig. 20.35(a). The geometry of motion of the carriage is shown in Fig. 20.35(b). The carriage goes around a circle of radius r=OC with constant speed v=ωr. The only acceleration that the carriage has is the centripetal acceleration and at the instant of interest 𝒂=ω2rıˆ.

Filename:sfig8-7-2again
Figure 20.35:

The linear momentum balance (𝑭=𝑳˙) for the carriage gives:

T𝝀ˆCBmgȷˆ = m𝒂
or T(sinθıˆ+cosθȷˆ)mgȷˆ = mω2rıˆ (20.20)

Scalar equations from eqn. (20.20) are:

[eqn. (20.20)]ȷˆ  Tcosθmg = 0
 T = mgcosθ
= 100kg9.8m/s23/2
= 1133N.
[eqn. (20.20)]ıˆ  Tsinθ = mω2r
 ω2 = Tsinθmr
= Tsinθm(R+lsinθ)
= 1133N12100kg(4m+5m12)
= 0.871s2
 ω = 0.93rads
= 0.931s1rev2π60s1min
= 8.9rpm.

Answer: T=1133N,ω=8.9rpm

Alternatively,

we could also find the angular speed using angular momentum balance. The angular momentum balance about point B gives

𝑴/B=𝑯˙/B
𝑴/B = 𝒓C/B×(mgȷˆ)
= mglsinθ𝒌ˆ
𝑯˙/B = 𝒓C/B×(mω2rıˆ)
= mω2rlcosθ𝒌ˆ

Equating the two quantities, we get

mω2rlcosθ = mglsinθ
 ω2 = grtanθ
= gtanθR+lsinθ
= 9.8m/s20.5774m+5m12
= 0.871s2
ω = 0.93s1=8.9rpm

which is the same value as we found using the linear momentum balance .

SAMPLE 20.6

Filename:sfig8-7-2disks
Figure 20.36: A bar held by two strings rotates in 3-D.

A crooked bar rotating with a shaft in space. A uniform rod CD of mass m=2kg and length =1m is fastened to a shaft AB by means of two strings: AC of length R1=30cm, and BD of length R2=50cm. The shaft is rotating at a constant angular velocity 𝝎=5rad/s𝒌ˆ. There is no gravity. At the instant shown, find the tensions in the two strings.

margin:

Solution The free-body diagram of the rod is shown in Fig. 20.37. The linear momentum balance (𝑭=m𝒂 ) for the rod gives:

T1+T2=mω2rG. (20.21)
Filename:sfig8-4-4
Figure 20.37:

Usually, linear momentum balance gives us two scalar equations in 2-D and three scalar equations in 3-D. Unfortunately, in this case, it gives only one equation for two unknowns T1 and T2. Therefore, we need one more equation.

The angular momentum balance about point D gives:

𝑴/D = 𝑯˙/D,
where𝑴/D = 𝒓C/D×(T1ȷˆ)
= (sinθȷˆ+cosθ𝒌ˆ)×(T1ȷˆ)
= T1cosθıˆ,
and
𝑯˙/D = m𝒓P/D×(ω2RPȷˆ)𝑑m
= 0l(sinθȷˆ+cosθ𝒌ˆ)𝒓P/D×(ω2(R2lsinθ)RPȷˆ)mdldm
= mω2(R2cosθ0l𝑑lcosθsinθ0l2𝑑l)ıˆ
= mω2(R2cosθ22cosθsinθ33)ıˆ
= mω2cosθ(12R213sinθ)ıˆ.

Thus,

T1cosθ = mω2cosθ(12R213sinθ)
 T1 = mω2(12R213sinθ).

Substituting in (20.21) we get

T2=mω2(rG12R2+13sinθ).

Plugging in the given numerical values and noting that rG=(R1+R2)/2=40cm and sinθ=R2R1=20cm, we get

T1 = 2kg(51s)2(0.4m120.5m+130.2m)
= 9.17kgms2=9.17N
and T2 = 10.83N.

Answer: T1=9.1N,T2=10.9N

SAMPLE 20.7

Filename:sfig8-4-4a
Figure 20.38: A rectangular plate, mounted rigidly at an angle ϕ on a shaft, wobbles as the shaft rotates at a constant speed.

A crooked plate rotating with a shaft in space. A rectangular plate of mass m, length , and width b is welded to a shaft AB in the center. The long edge of the plate is parallel to the shaft axis but is tipped by an angle ϕ with respect to the shaft axis. The shaft rotates with a constant angular speed ω. The end B of the shaft is free to move in the z-direction. Assume there is no gravity. Find the reactions at the supports.

Solution A simple line sketch and the Free-Body Diagram of the system are shown in Fig. 20.39ab.

Filename:sfig8-4-4b
Figure 20.39:

The linear momentum balance equation for the shaft and the plate system is:

𝑭=mtotal𝒂cm.

Since the center-of-mass is on the axis of rotation, 𝒂cm=𝟎. Therefore,

(Ax+Bx)ıˆ+(Ay+By)ȷˆ+Az𝒌ˆ=𝟎
  Ax+Bx=0,Ay+By=0,Az=0. (20.22)

The angular momentum balance about the center-of-mass O is:

𝑴O=𝑯˙/O
  • Calculation of MO:

    𝑴O = 𝒓A/O×𝑭A+𝒓B/O×𝑭B+Mz𝒌ˆ (20.23)
    = 2𝒌ˆ×(Axıˆ+Ayȷˆ+Az𝒌ˆ)+2𝒌ˆ×(Bxıˆ+Byȷˆ)+Mz𝒌ˆ
    = 2(AyBy)ıˆ+2(BxAx)ȷˆ+Mz𝒌ˆ
  • Calculation of H˙/O: 𝑯˙/O can be computed in various ways. margin: 𝑯˙ could also be computed using the moment of inertia matrix of the body. See the next two text sections. Here, to compute 𝑯˙/O, we use

    𝑯˙/O=M𝒓dm/O×𝒂dm𝑑m,

    the formula which we have used so far. To carry out this integration for the plate, we take, as usual, an infinitesimal mass dm of the body, calculate its angular momentum about O, and then integrate over the entire mass of the body:

    𝑯˙/O=M𝒓dm/O×𝒂dm𝑑m

    We need to write carefully each term in the integrand. Let us define an axis w (don’t confuse this dummy variable w with ω) along the length of the plate (see Fig. 20.40(a)). We take an area element dA=dwdy on the plate as our infinitesimal mass. Fig. 20.40(b) shows this element and its coordinates.

    Filename:sfig8-4-4c
    Figure 20.40: Calculation of 𝑯˙/O: (a) mass element dm is shown on the plate, (b) the mass element as an area element and its geometry.
    dm = ρdA=mbdwdy(ρ=mass per unit area)
    𝒂dm = 𝝎×(𝝎×𝒓dm/O)
    𝒓dm/O = xıˆ+yȷˆ+z𝒌ˆ

    where

    x=wsinϕ,y=y,z=wcosϕ. (20.24)

    Therefore,

    𝒂cm = ω𝒌ˆ×(ω𝒌ˆ×(wsinϕıˆ+yȷˆ+wcosϕ𝒌ˆ))
    = ω𝒌ˆ×(ωwsinϕȷˆωyıˆ)=ω2(wsinϕıˆ+yȷˆ),
    𝒓dm/O×𝒂dm = (wsinϕıˆ+yȷˆ+wcosϕ𝒌ˆ)×[ω2(wsinϕıˆ+yȷˆ)]
    = ω2(w2sinϕcosϕȷˆ+wycosϕıˆ).

    Thus,

    𝑯˙/O = b/2b/2/2/2ω2(w2sinϕcosϕȷˆ+wycosϕıˆ)𝒓dm/O×𝒂dmmbdwdydm (20.25)
    = mbω2b/2b/2(/2/2(w2sinϕcosϕȷˆ+wycosϕıˆ)𝑑w)𝑑y
    = mbω2b/2b/2(sinϕcosϕw33|/2/2+ycosϕw22|/2/20)𝑑y
    = mω2212sinϕcosϕȷˆ.
  • Now, back to angular momentum balance: Now equating (20.23) and (20.25) and dotting both sides with ıˆ, ȷˆ, and 𝒌ˆ we get

    AyBy=0,BxAx=mω26sinϕcosϕ,Mz=0, (20.26)

    respectively. Solving (20.22) and (20.26) simultaneously we get

    Ay=By=0,Ax=mω212,sinϕcosϕBx=mω212sinϕcosϕ.

Answer: Ax=mω212sinϕcosϕ,Bx=mω212sinϕcosϕ,Ay=By=Az=Mz=0

SAMPLE 20.8

Filename:pfigure-s94h13p2
Figure 20.41: A short rod swings in 3-D.

A short rod as a 3-D pendulum. A uniform rod AB of mass m and length 2 is welded to a massless, inextensible thin rod OA at point A. Rod OA is attached to a ball and socket joint at point O. The rods are going around in a circle with constant speed maintaining a constant angle θ with the vertical axis. Assume that θ is small.

  1. 1.

    How many revolutions does the system make in one second?

Solution

  1. 1.

    The free-body diagram of the system (rod OA + rod AB) is shown in Fig. 20.42. Let ω be the angular speed of the system. Then, the number of revolutions in one second is n=ω/(2π). Therefore, to find the answer we need to calculate ω.

    Filename:pfigure-f93f5
    Figure 20.42:

    The angular momentum balance about point O gives 𝑴O=𝑯˙/O. Now,

    𝑴O = 𝒓G/O×(mg𝒌ˆ)
    = 4(sinθȷˆcosθ𝒌ˆ)×(mg𝒌ˆ)
    = 4mgsinθıˆ.
    and 𝑯˙/O = m𝒓dm/O×𝒂dm𝑑m
    = 35(sinθȷˆcosθ𝒌ˆ)𝒓dm/O×(ω2lsinθȷˆ)𝒂dmm2ddm
    = m2ω2sinθcosθıˆ352𝑑
    = 4932mω2sinθcosθıˆ.

    By equating the two quantities (𝑴O=𝑯˙/O ), we get

    4mgsinθıˆ = 4932mω2sinθcosθıˆ
     ω2 = 12g49cosθ.

    But for small θ, cosθ1. Therefore,

    ω=12g49=237g

    and the number of revolutions per unit time is n=2314πg.

    Answer: n=2314πg

  2. 2.

    Note, the natural frequency of this rod swinging back and forth as a simple pendulum turns out to be the same as the angular speed ω of the rotating system above for small θ.

Problems for 20.2 Dynamics of fixed axis rotation

20.2.1  Let P(x,y,z) be a point mass on an abstract massless body that rotates at a constant rate ω about the z-axis. Let C(0,0,z) be the center of the circular path that P describes during its motion.

  1. (a)

    Show that the rate of change of angular momentum of mass P about the origin O is given by

    𝑯˙/O=mω2Rz𝒆ˆθ

    where R=x2+y2 is the radius of the circular path of point P and m is its mass.

  2. (b)

    Is there any point on the axis of rotation (z-axis) about which 𝑯˙ is zero. Is this point unique?

Filename:pfigure-s94h13p3
Figure 20.43

20.2.2   System of Particles going in a circle. Assume for the problems below that you have two mass points going in circles about a common axis at a fixed common angular rate.

  1. (a)

    2-D (particles in the plane, axis of rotation normal to the plane). Show, by calculating the sum in 𝐇˙C that the result is the same when using the two particles or when using a single particle at the center of mass (with mass equal to the sum of the two particles). You may do the general problem or you may pick particular appropriate locations for the masses, particular masses for the masses, a specific point C, a specific axis and a particular rotation rate.

  2. (b)

    3-D (particles not in a common plane normal to the axis of rotation). Using particular mass locations that you choose, show that 2 particles give a different value for 𝐇˙C than a single particle at the center of mass. [Moral: In general, 𝑯˙/C𝒓cm/C×mtot𝒂cm.]

20.2.3  The rotating two-mass system shown in the figure has a constant angular velocity 𝝎=12rad/sıˆ. Assume that at the instant shown, the masses are in the xy-plane. Take m=1kg and =1m.

  1. (a)

    Find the linear momentum 𝑳 and its rate of change 𝑳˙ for this system.

  2. (b)

    Find the angular momentum 𝑯/O about point O and its rate of change 𝑯˙/O for this system.

Filename:p-s96-p3-3
Figure 20.44

20.2.4  Particle attached to a shaft, no gravity. A shaft of negligible mass spins at constant rate. The motor at A imposes this rate if needed. The bearing at A prevents translation of its end of the shaft in any direction but causes no torques other than that of the motor (which causes a torque in the 𝒌ˆ direction). The frictionless bearing at C holds that end of the shaft from moving in the ıˆ and ȷˆ directions but allows slip in the 𝒌ˆ direction. The shaft is welded orthogonally to a bar, also of negligible mass, which is attached to a small ball of (non-negligible) mass of 1kg. There is no gravity. The shaft is spinning counterclockwise (when viewed from the positive 𝒌ˆ direction) at 3rad/s. In the configuration of interest the rod BD is parallel with the ıˆ direction.

  1. (a)

    What is the tension in the rod BD? Answer: TDB=0.9N.

  2. (b)

    What are all the reaction forces at A and C? Answer: FA=0.3Nıˆ, 𝑭C=0.6Nıˆ

  3. (c)

    What is the torque that the motor at A causes on the shaft? Answer: Torquemotor=0.

  4. (d)

    There are two ways to do problems (b) and (c): one is by treating the shaft, bar and mass as one system and writing the momentum balance equations. The other is to use the result of (a) with action-reaction to do statics on the shaft BC. Repeat problems (b) and (c) doing it the way you did not do it the first time.

Filename:bikefork1-ang-mom
Figure 20.45

20.2.5   Particle attached to a shaft, with gravity. Consider the same situation as in problem 20.2 with the following changes: There is gravity (pointing in the 𝐢 direction). At the time of interest the rod BD is parallel to the +𝐣 direction.

  1. (a)

    What is the force on D from the rod BD? (Note that the rod BD is not a ‘two-force’ member. Why not?)

  2. (b)

    What are all the reaction forces at A and C?

  3. (c)

    What is the torque that the motor at A causes on the shaft?

  4. (d)

    Do problems (b) and (c) again a different way.

20.2.6  A turntable rotates at the constant rate of ω=2rad/s about the z axis. At the edge of the 2m radius turntable swings a pendulum which makes an angle ϕ with the negative z direction. The pendulum length is =10m. The pendulum consists of a massless rod with a point mass m=2kg at the end. The gravitational constant is g=10m/s2 pointing in the z direction. A perspective view and a side view are shown.

  1. (a)

    What possible angle(s) could the pendulum make in a steady motion. In other words, find ϕ for steady circular motion of the system. It turns out that this quest leads to a transcendental equation. You are not asked to solve this equation but you should have the equation clearly defined. Answer: Angles ϕ for steady circular motion are solutions to 10tanϕ=8+40sinϕ.

  2. (b)

    How many solutions are there? [hint, sketch the functions involved in your final equation.] Answer: ϕ=[1.3627rad,0.2735rad,1.2484rad,2,9813rad], four solutions, found numerically. Graphically, the solutions are the intersections of the curves f1(ϕ)=10tanϕ and f2(ϕ)=8+40sinϕ between π and π.

  3. (c)

    Sketch very approximately any equilibrium solution(s) not already sketched on the side view in the picture provided.

Filename:bikefork-ang-mom
Figure 20.46

20.2.7  Mr. X swings his daughter Y about the vertical axis at a constant rate ω=1rad/s. Y is merely one year old and weighs 12kg. Assuming Y to be a point mass located at the center of mass, estimate the force on each of her shoulder joints.

Filename:summer95f-5-a
Figure 20.47

20.2.8  A block of mass m=2kg supported by a light inextensible string at point O sits on the surface of a rotating cone as shown in the figure. The cone rotates about its vertical axis of symmetry at constant rate ω. Find the value of ω in rpm at which the block loses contact with the cone.

Filename:pfigure4-2-rp10
Figure 20.48

20.2.9  A small rock swings in circles at constant angular rate ω. It is hung by a string with length the other end of which is at the origin of a fixed coordinate system. The rock has mass m and the gravitational constant is g. The string makes an angle ϕ with the negative z-axis. The mass rotates so that if you follow it with the fingers of your right hand, your right thumb points up the z-axis. At the instant of interest the mass is passing through the yz-plane.

  1. (a)

    How long does the rock take to make one revolution? How does this compare with the period of oscillation of a simple pendulum. Answer: Period of revolution=2πω, Period of simple pendulum (for small oscillations)=2πLg, where is the length of the pendulum.

  2. (b)

    Find the following quantities in any order that pleases you. All solutions should be in terms of m, g, , ϕ, ω and the unit vectors ıˆ, ȷˆ, and 𝒌ˆ.

    • The tension in the string (T).

    • The velocity of the rock (𝒗).

    • The acceleration of the rock (𝒂).

    • The rate of change of the angular momentum about the +x axis (dHx/dt).

    Answer: T=mω2L (strange that ϕ drops out!), 𝒗rock=ωLsinϕıˆ, 𝒂rock=ω2Lsinϕȷˆ, dHxdt=mωL2sinϕcosϕ.

Filename:pfigure-blue-125-2
Figure 20.49

20.2.10  Consider the configuration in problem 20.2. What is the relation between the angle (measured from the vertical) at which the rock hangs and the speed at which it moves? Do an experiment with a string, weight, and stopwatch to check your calculation. Answer: ϕ=cos1(gω2L).

20.2.11  A small block of mass m=5kg sits on an inclined plane which is supported by a bearing at point O as shown in the figure. At the instant of interest, the block is at a distance =0.5m from where the surface of the incline intersects the axis of rotation. By experiments, the limiting rate of rotation of the turntable before the block starts sliding up the inclined plane is found to be ω=125rpm. Find the coefficient of friction between the block and the incline.

Filename:pfigure-blue-68-1
Figure 20.50

20.2.12  Particle sliding in circles in a parabolic bowl. As if in a James Bond adventure in a big slippery radar bowl, a particle-like human with mass m is sliding around in circles at speed v. The equation describing the bowl is z=CR2=C(x2+y2).

  1. (a)

    Find v in terms of any or all of R, g, and C.

  2. (b)

    Now say you are given ω,C and g. Find v and R  if you can. Explain any oddities.

Filename:pfigure-blue-58-1
Figure 20.51

20.2.13  Mass on a rotating inclined platform. A block of mass m=5kg is held on an inclined platform by a rod BC as shown in the figure. The platform rotates at a constant angular rate ω. The coefficient of friction between the block and the platform is μ. θ=30o, =1m, and g=9.81m/s2. (a) Assuming ω=5rad/s and μ=0 find the tension in the rod. (b) If μ=.3 what is the range of ω for which the block would stay on the ramp without slipping even if there were no rod? Answer: (a) T=69.2N, (b) 1.64rad/s<ω<3.47rad/s

Filename:pfigure-blue-157-1
Figure 20.52

20.2.14  A circular plate of radius r=10cm rotates in the horizontal plane about the vertical axis passing through its center O. Two identical point masses hang from two points on the plate diametrically opposite to each other as shown in the figure. At a constant angular speed ω, the two masses maintain a steady state angle θ from their vertical static equilibrium.

  1. (a)

    Find the linear speed and the magnitude of the centripetal acceleration of either of the two masses.

  2. (b)

    Find the steady state ω as a function of θ. Do the limiting values of θ, θ0 and θπ/2 give sensible values of ω?

Filename:summer95f-5
Figure 20.53

20.2.15  For the rod and mass system shown in the figure, assume that θ=60, ωθ˙=5rad/s, and =1m. Find the acceleration of mass B at the instant shown by

  1. (a)

    calculating 𝒂=𝝎×𝝎×𝒓 where 𝒓 is the position vector of point B with respect to point O, i.e., 𝒓=𝒓B/O, and by

  2. (b)

    calculating 𝒂=Rθ˙2𝒆ˆR where R is the radius of the circular path of mass B.

Filename:pfigure-blue-127-2
Figure 20.54

20.2.16  A spring of relaxed length 0 (same as the diameter of the plate) is attached to the mid point of the rods supporting the two hanging masses.

  1. (a)

    Find the steady angle θ as a function of the plate’s rotation rate ω and the spring constant k. (Your answer may include other given constants.)

  2. (b)

    Verify that, for k=0, your answer for θ reduces to that in problem 20.53.

  3. (c)

    Assume that the entire system sits inside a cylinder of radius 20cm and the axis of rotation is aligned with the longitudinal axis of the cylinder. The maximum operating speed of the rotating shaft is 200rpm (i.e., ωmax=200rpm). If the two masses are not to touch the walls of the cylinder, find the required spring stiffness k.

Filename:pfigure-s94h13p4
Figure 20.55

20.2.17  Flyball governors are used to control the flow of working fluid in steam engines, diesel engines, steam turbines, etc. A spring-loaded flyball governor is shown schematically in the figure. As the governor rotates about the vertical axis, the two masses tend to move outwards. The rigid and massless links that connect the outer arms to the collar, in turn, try to lift the central collar (mass m/4) against the spring. The central collar can move along the vertical axis but the top of the governor remains fixed. The mass of each ball is m=5 kg, and the length of each link is =0.25 m. There are frictionless hinges at points A, B, C, D, E, F where the links are connected. Find the amount the string is elongated, Δx, when the governor rotates at 100 rpm, given that the spring constant k=500 N/m, and that in this steady state, θ=90o. Answer: Δx=0.554m

Filename:pfigure-blue-90-2
Figure 20.56

20.2.18  Crooked rod. A 2 long uniform rod of mass m is welded at its center to a shaft at an angle ϕ as shown. The shaft is supported by bearings at its ends and spins at a constant rate ω.

  1. (a)

    What are 𝑯˙/O and 𝑳˙?

  2. (b)

    Calculate 𝐇˙/O using the integral definition of angular momentum. Answer: H˙/O=112mω22sinϕcosϕkˆ.

  3. (c)

    What are the six reaction components if there is no gravity? Answer: FA=(112mω22/dsinϕcosϕ)ȷˆ, 𝑭B=𝑭A.

  4. (d)

    Repeat the calculations, including gravity g. Which reactions change and by how much? Answer: FA=(112mω22/dsinϕcosϕ+12mg)ȷˆ, 𝑭B𝑭A.

Filename:s92f1p7
Figure 20.57

20.2.19  For the configuration in problem 20.2, find two positions along the rod where concentrating one-half the mass at each of these positions gives the same 𝑯˙/O.

20.2.20  Rod on a shaft. A uniform rod is welded to a shaft which is connected to frictionless bearings. The rod has mass m and length and is attached at an angle γ at its midpoint to the midpoint of the shaft. The length of the shaft is 2d. A motor applies a torque Mmotor along the axis of the shaft. At the instant shown in the figure, assume m, , d, ω and Mmotor are known and the rod lies in the yz-plane.

  1. (a)

    Find the angular acceleration of the shaft.

  2. (b)

    The moment that the shaft applies to the rod.

  3. (c)

    The reaction forces at the bearings.

Filename:twodisks-ang-mom
Figure 20.58

Rod on a shaft.

20.2.21  Rod spins on a shaft. A uniform rod with length R and mass m spins at constant rate ω about the z axis. It is held by a hinge at A and a string at B. Neglect gravity. Find the tension in the string in terms of ω, R and m. Answer: T=2mω2R6.

Filename:pfigure-blue-90-1
Figure 20.59

20.2.22  A uniform narrow shaft of length 2 and mass 2M is bent at an angle θ half way down its length. It is spun at constant rate ω.

  1. (a)

    What is the magnitude of the force that must be applied at O to keep the shaft spinning?

  2. (b)

    (harder) What is the magnitude of the moment that must be applied at O to keep the shaft spinning?

Filename:pfigure-blue-110-1
Figure 20.60

Spinning rod. No gravity.

20.2.23  For the configuration in problem 20.2, find two positions along the bent portion of the shaft where concentrating one-half its mass at each of these positions gives the same 𝑯˙/O and 𝑳˙.

20.2.24  Rod attached to a shaft. A uniform 1kg, 2m rod is attached to a shaft which spins at a constant rate of 1rad/s. One end of the rod is connected to the shaft by a ball-and-socket joint . The other end is held by two strings (or pin jointed rods) that are connected to a rigid cross bar. The cross bar is welded to the main shaft. Ignore gravity.

  1. (a)

    Draw a free-body diagram of the rod AB, including a bit of the strings.

  2. (b)

    Find the tension in the string DB. Two approaches:

    • (i)

      Write the six equations of linear and angular momentum balance (about any point of your choosing) and solve them (perhaps using a computer) for the tension TDB.

    • (ii)

      Write the equation of angular momentum balance about the axis AE. (That is, take the equation of angular momentum balance about point A and dot both sides with a vector in the direction of 𝒓E/A.) Use this calculation to find TDB and compare to method ii.

  3. (c)

    Find the reaction force (a vector) at A.

Filename:pfigure-blue-107-1
Figure 20.61

Uniform rod attached to a shaft.

20.2.25  3D, crooked bar on shaft falls down. A massless, rigid shaft is connected to frictionless hinges at A and B. A bar with mass Mb and length is welded to the shaft at 45o. At time t=0 the assembly is tipped a small angle θ0 from the yz-plane but has no angular velocity. It then falls under the action of gravity. What is θ at t=t1 (assuming t1 is small enough so that θ1 is much less than 1)?

Filename:pfigure-s94h14p4
Figure 20.62

A crooked rod welded to a shaft falls down.

20.3 Moment of inertia matrices [𝑰cm] and [𝑰O]

We now know how to find the velocity and acceleration of every bit of mass on a rigid body as it spins about a fixed axis. It is just a matter of doing integrals or sums to calculate the various motion quantities (momenta, energy) of interest. As the body moves and rotates the region of integration and the values of the integrands change. So, in principle, in order to analyze a rigid body one has to evaluate a different integral or sum at every different configuration. But there is a shortcut. A big sum (over all atoms, say), or a difficult integral is reduced to a simple multiplication using the moment of inertia. In three-dimensions this multiplication is a matrix multiplication.

[𝐈], the moment of inertia matrixmargin: In fact the moment of inertia matrix for a given object depends on what reference point is used. Most commonly when people say ‘the’ moment of inertia they mean to use the center-of-mass as the reference point. For clarity this moment of inertia matrix is often written as [𝑰cm] in this book. If a different reference point, say point O is used, the matrix is notated as [𝑰O]. , is defined for the purpose of simplifying the expressions for the angular momentum, the rate of change of angular momentum, and the energy of a system which moves like a rigid body.

First review the situation for flat objects in planar motion. A flat object spinning with 𝝎=ω𝒌ˆ in the xy plane has a mass distribution which gives a polar moment of inertia Izzcm or just ‘I’ so that:

𝑯cm = Iω𝒌ˆ (20.27)
𝑯˙cm = 𝟎 (20.28)
EK/cm = 12ω2I. (20.29)

Now, for a rigid body spinning in 3-D about a fixed axis with the angular velocity 𝝎 we need matrix multiplication, where the determination of the needed matrix is the central topic of this section.

𝑯cm = [Icm]𝝎 (20.30)
𝑯˙cm = 𝝎×[Icm]𝝎𝑯cm+[Icm]𝝎˙ (20.31)
EK/cm = 12𝝎([Icm]𝝎)=12𝝎𝑯cm. (20.32)

In detail, for example,

[Hx/cmHy/cmHz/cm]=[IxxcmIxycmIxzcmIxycmIyycmIyzcmIxzcmIyzcmIzzcm]xyz[ωxωyωz] (20.33)

where  𝑯cm=Hx/cmıˆ+Hy/cmȷˆ+Hz/cm𝒌ˆ. Note that the 2-D results are a special case of the 3-D results because, as you will soon see, for 2-D objects Ixz=Iyz=0. We postpone the use of these equations until section 20.4.

The moment of inertias in 3-D: [𝑰cm] and [𝑰O]

For the study of three-dimensional mechanics, including the simple case of constant rate rotation about a fixed axis, one often makes use of the moment of inertia matrix, defined below and motivated by the box 20.3 on page 20.3.

The distances x,y,z in the formulas below are the x,y,z components of the position of mass relative to a coordinate system which has either the center-of-mass (cm) or the point O as its origin. margin: Caution: While we have Ixy=xy𝑑m, some old books define Ixy=xy𝑑m. They then have minus signs in front of the off-diagonal terms in the moment of inertia matrix. They would say I12=Ixy. The numerical values in the matrix they write is the same as in the one we write. They just have a different sign convention in the definition of the components.

[𝑰] = [IxxIxyIxzIxyIyyIyzIxzIyzIzz] (20.34)
= [(y2+z2)𝑑mxy𝑑mxz𝑑mxy𝑑m(x2+z2)𝑑myz𝑑mxz𝑑myz𝑑m(x2+y2)𝑑m] (20.35)

If all mass is on the xy plane then it is clear that Ixz=Iyz=0 since z=0 for the whole xy plane. If rotation is also about the z-axis then 𝝎=ω𝒌ˆ.

Applying the formulas above we find that, if all of the mass is in the xy-plane and rotation is about the z-axis, the only relevant non-zero term in [𝑰] is Izz=(x2+y2)𝑑m . And, Ixx, Iyy, and Ixy don’t contribute to 𝑯, 𝑯˙, or EK. In this manner you can check that the three dimensional equations, when applied to two-dimensional bodies, give the same results that we found directly for two-dimensional bodies.

Example: Moment of inertia matrix for a uniform sphere

A sphere is a special shape which is, naturally enough, spherically symmetric. Therefore,

Ixxcm=Iyycm=Izzcm

and

Ixycm=Ixzcm=Iyzcm=0.

So, all we need is Ixxcm or Iyycm or Izzcm.

Filename:pfigure-blue-112-1
Figure 20.63:

Here is the trick:

Ixxcm = 13(Ixxcm+Iyycm+Izzcm)
= 13[(y2+z2)𝑑m+(x2+z2)𝑑m+(x2+y2)𝑑m]
= 23(x2+y2+z2)𝑑m
= 23r2𝑑m
= 230Rr2(4ρπr2dr)
= 83ρπ0Rr4𝑑r
= 815ρπR5
= 25mR2(m=43ρπR3).

So,

[𝑰cm]=25mR2[100010001].

The parallel axis theorem for rigid bodies in three dimensions

The 3-D parallel axis theorem is stated below and in the table on the inside back cover. It is derived in box 20.3 on page 20.3. The parallel axis theorem for rigid bodies in three dimensions is the equation

[𝑰O]=[𝑰cm]+m[ycm/o2+zcm/o2xcm/oycm/oxcm/ozcm/oxcm/oycm/oxcm/o2+zcm/o2ycm/ozcm/oxcm/ozcm/oycm/ozcm/oxcm/o2+ycm/o2] (20.36)

In this equation, xcm/o, ycm/o, and zcm/o are the x, y, and z coordinates, respectively, of the center-of-mass defined with respect to a coordinate system whose origin is located at some point O not at the center-of-mass cm. That is, if you know [𝑰cm], you can find [𝑰O] without doing any more integrals or sums. Like the 2-D parallel axis theorem, the primary utility of the 3-D parallel axis theorem is for the determination of [𝐈] for an object that is a composite of simpler objects. Such are not beyond the scope of this book in principle. But in fact, given the finite time available for calculation, we do not leave much time for practice of this tedious but routine calculation.

Matrices and tensors

We have just introduced the 3 by 3 moment of inertia matrix [𝑰]. We will find it in expressions having to do with angular momentum sitting next to either a vector 𝝎 or a vector 𝜶: [𝑰]𝝎or[𝑰]𝜶. What we mean by this expression is the three element column vector that comes from matrix multiplication of the matrix [𝑰] and the column vector for 𝝎, [ωxωyωz], an expression that only makes sense if everyone knows what bases are being used.

More formally, and usually only in more advanced treatments, people like to define a coordinate-free quantity called the tensor 𝑰¯¯. Then we would have

𝑰¯¯𝝎

by which we would mean the vector whose components would be found by [𝑰][ωxωyωz].

Eigenvectors and Eigenvalues

A square matrix [𝑨] when multiplied by a column vector [𝒗] yields a new vector [𝒘]. A given matrix has a few special vectors, somehow characteristic of that matrix, called eigenvectors. The vector 𝒗 is an eigenvector of [𝑨] if

[𝑨][𝒗] is parallel to [𝒗].

In other words, if

[𝑨][𝒗]=λ[𝒗]

for some λ.

The scalar λ is called the eigenvalue associated with the eigenvector [𝒗] of the matrix [𝑨]. The eigen-values and eigen-vectors of a matrix are found with a single command in many computer math programs. In statics you had little or no use for eigen-values and eigen-vectors. In dynamics, eigenvectors and eigenvalues are useful for understanding dynamic balance, 3-D rigid body rotations, and normal mode vibrations.

Eigenvectors of [𝑰]

The moment of inertia matrix is always a symmetric matrix. This symmetry means that [𝑰] always has a set of three mutually orthogonal eigenvectors. The importance of the eigenvectors of [𝑰] will be discussed in section 20.5 on dynamic balance. Sometimes a pair of the eigenvalues are equal to each other implying that any vector in the plane of the corresponding eigenvectors is also an eigenvector.

If the physical object has any natural symmetry directions these directions will usually manifest themselves in the dynamics of the body as being in the directions of the eigenvectors of the object’s moment of inertia matrix. For example, the dotted lines on fig. 2.34 are all in directions of eigenvectors for the objects shown. But even if an object is wildly asymmetric in shape, its moment of inertia matrix is always symmetric and thus all objects have moment of inertia matrices with at least three different eigenvectors at least three of which are mutually orthogonal.

Properties of [𝑰]

For those with experience with linear algebra various properties of the moment of inertia matrix [𝑰] are worth noting (although not worth proving here). Unless all mass is distributed on one straight line, the moment of inertia matrix is invertible (it is non-singular and has rank 3). Further, when invertible it is positive definite. In the special case that all the mass is on some straight line, the moment of inertia matrix is non-invertible and only positive semi-definite. The positive (semi) definiteness of the moment of inertia matrix is equivalent to the statement that the rotational kinetic energy of a body is always equal to or greater than zero. Finally, the eigenvalues of the moment of inertia matrix are all positive and have the property that no one can be greater than the sum of the other two (the same inequalities are satisfied by the lengths of the sides of a triangle, the “triangle inequality”).

Box 20.1 Discovering the moment of inertia matrix

Derivation 1

Here we present a direct derivation of the moment of inertia matrix; that is, a derivation in which the moment of inertia matrix arises as a convenient short hand. Assume a rigid body is moving in such a way that point O is fixed (i.e., it is either on the line of a hinge or a ball-and-socket joint).


Filename:tfigure4-inertia-fig0



The most basic kinematic relation for a rigid body is that

𝒗=𝝎×𝒓

where 𝒓/O=xıˆ+yȷˆ+z𝒌ˆ is the position of a point on the body relative to O and 𝝎, the angular velocity of the body.

Now, we tediously calculate and arrange the terms in the angular momentum about point O,

𝑯O = 𝒓/O×𝒗𝑑m
= 𝒓/O×(𝝎×𝒓/O)𝑑m
= (xıˆ+yȷˆ+z𝒌ˆ)×
[(ωxıˆ+ωyȷˆ+ωz𝒌ˆ)×(xıˆ+yȷˆ+z𝒌ˆ)]dm
= (xıˆ+yȷˆ+z𝒌ˆ)×
[(ωyzωzy)ıˆ+(ωzxωxz)ȷˆ+(ωxyωyx)𝒌ˆ]dm
= [(y(ωxyωyx)z(ωzxωxz))ıˆ
+(z(ωyzωzy)x(ωxyωyx))ȷˆ
+(x(ωzxωxz)y(ωyzωzy))𝒌ˆ]dm
= [((y2+z2)ωxxyωyxzωz)ıˆ
+(yxωx+(z2+x2)ωyyzωz)ȷˆ
+(zxωxxyωy+(x2+y2)ωz)𝒌ˆ]dm

Since the integral is over the mass and 𝝎 is constant over the body, we can pull 𝝎 out of the integral so that we may write the equation in matrix form. Writing 𝑯O as a column vector, we can rewrite the last equation as

[HOxHOyHOz]=[(y2+z2)𝑑mxy𝑑mxz𝑑mxy𝑑m(x2+z2)𝑑myz𝑑mxz𝑑myz𝑑m(x2+y2)𝑑m][𝑰O][ωxωyωz].


Finally, defining [𝑰O] by the matrix above, we can compactly write

𝑯/O=[𝑰O]𝝎

assuming O is a fixed point on the body where we represent 𝑯/O and 𝝎 in terms of x, y, and z components.

Center-of-mass inertia matrix

For any system moving, distorting, and rotating any crazy way, we have the general result that

𝑯/O=𝒓cm/O×mtot𝒗cmContribution of the system to 𝑯/O if treated as a particle at the system center-of-mass+(𝒓/cm×𝒗/cm=(𝒗𝒗cm)𝒗/cm)𝑑m𝑯cm

as you can verify by substituting 𝒗=𝒗cm+𝒗/cm and 𝒓=𝒓cm+𝒓/cm into the general definition of 𝑯/O=𝒓/O×𝒗𝑑m. For a rigid body, we have

𝒗/cm=𝝎×𝒓/cm.

So, by a derivation essentially identical to that for 𝑯/O, we get

𝑯cm=[𝑰cm]𝝎

with [𝑰cm] being defined using x, y, and z, as the distances from the center of mass rather than from point O. So, for a rigid body in general motion, we can find the angular momentum by

𝑯/O=𝒓cm/O×𝒗cmmtot+[𝑰cm]𝝎 (20.37)

Comment (aside)

In the special case that the body is rotating about point O, we also have

𝑯/O=[𝑰O]𝝎. (20.38)

You will see, if you look at the parallel axis theorem, that these two expressions 20.37 and 20.38 do in fact agree.

Derivation 2

Here, we present a less direct but perhaps more intuitive derivation of the moment of inertia matrix. We start with the special case of a 3-D rigid body spinning in circles at constant rate about a fixed axis.

To see from where the moment of inertia matrix comes, we will first calculate the angular momentum about point O of a general 3-D rigid body spinning about the z-axis with constant rate θ˙const.=ωz or 𝝎=ωz𝒌ˆ. We will refer to this case as (1).

Filename:tfigure4-inertia-fig2

Starting with the definition of angular momentum, we get

𝑯O = 𝒓/O×𝒗𝑑m
= (z𝒌ˆ+𝑹)×(θ˙R𝒆ˆθ)𝑑m.
Filename:tfigure4-inertia-fig1

Looking down the z-axis in the figure, we see that

𝑹 = R𝒆ˆR=Rcosθxıˆ+Rsinθyȷˆ,or
𝑹 = xıˆ+yȷˆ.

where R=x2+y2.

To compute the cross product in the integrand, we need

𝒌ˆ×𝒆ˆθ = 𝒌ˆ×(sinθıˆ+cosθȷˆ)=𝒆ˆRand
𝑹×𝒆ˆθ = R𝒆ˆR×𝒆ˆθ=R𝒌ˆ.

Therefore, we now have

𝑯O=zθ˙R𝒆ˆRdm+θ˙R2𝒌ˆ𝑑m=θ˙[z(Rcosθxıˆ+Rsinθyȷˆ)𝑑m+(x2+y2)𝒌ˆ𝑑m]=θ˙[(zx𝑑m)ıˆ+(zy𝑑m)ȷˆ+((x2+y2)𝑑m)𝒌ˆ].

To ‘un-clutter’ this expression, let’s define the following:

IxzO = xz𝑑m
IyzO = yz𝑑m
IzzO = (x2+y2)𝑑m

So, now, we have for case (1)

(𝑯O)1=IxzOωzıˆ+IyzOωzȷˆ+IzzOωz𝒌ˆ.

The substitutions we have defined form the elements of the third column of the inertia matrix, as we will see below in the general case. Let’s now move on to general 3-D rigid body motion and infer the first and second columns of the inertia matrix.

In general, the angular velocity of a rigid body is given by

𝝎=ωxıˆ+ωyȷˆ+ωz𝒌ˆ

So far, we have considered the special case above, 𝝎=ωz𝒌ˆ. But, we could have looked at 𝝎=ωxıˆ, case (2), and, similarly, would obtain instead the following angular momentum about point O

(𝑯O)2=IxxOωxıˆ+IyxOωxȷˆ+IzxOωx𝒌ˆ.

Likewise, for 𝝎=ωyȷˆ, case (3), we would obtain

(𝑯O)3=IxyOωyıˆ+IyyOωyȷˆ+IzyOωy𝒌ˆ.

Finally, for 𝝎=ωxıˆ+ωyȷˆ+ωz𝒌ˆ, we obtain

𝑯O=(𝑯O)1+(𝑯O)2+(𝑯O)3=IxzOωzıˆ+IyzOωzȷˆ+IzzOωz𝒌ˆ+IxxOωxıˆ+IyxOωxȷˆ+IzxOωx𝒌ˆ+IxyOωyıˆ+IyyOωyȷˆ+IzyOωy𝒌ˆ.

Collecting components, we get

𝑯O=HOxıˆ+HOyȷˆ+HOz𝒌ˆ

where

HOx=IxxOωx+IxyOωy+IxzOωz
HOy=IyxOωx+IyyOωy+IyzOωz
HOz=IzxOωx+IzyOωy+IzzOωz.

We can combine the above results into a matrix representation. Representing 𝑯O as a column vector, we can re-write the above set of three equations as a product of a matrix and the angular velocity written as a column vector.

[HOxHOyHOz]=[IxxOIxyOIxzOIyxOIyyOIyzOIzxOIzyOIzzO][ωxωyωz]

We define the coefficient matrix above to be the moment of inertia matrix about point O

[𝑰O]=[IxxOIxyOIxzOIyxOIyyOIyzOIzxOIzyOIzzO].

whose components are

IxxO=(y2+z2)𝑑mIxyO=xy𝑑mIxzO=xz𝑑mIyxO=yx𝑑mIyyO=(x2+z2)𝑑mIyzO=yz𝑑mIzxO=zx𝑑mIzyO=zy𝑑mIzzO=(x2+y2)𝑑m.

By inspection, one can see that IxyO=IyxO, IxzO=IzxO, and IyzO=IzyO. Thus, the inertia matrix is symmetric; i.e., [𝑰O]T=[𝑰O]. So, there are always at most only six, not nine, independent components in the inertia matrix to compute.

Box 20.2 3-D parallel axis theorem

In three dimensions, the two matrices [𝑰O] and [𝑰cm] are related to each other in a way similar to the two-dimensional case. Since the inertia matrix has six independent entries in it, the derivation involves six integrals. Let’s look at a typical term on the diagonal, say, IzzO and a typical off-diagonal term, say, IxyO. The calculation for the other terms is similar with a simple change of letters in the subscript notation.

First, IzzO=(x/O2+y/O2)𝑑m=Izzcm+m(xcm/O2+ycm/O2) by exactly the same reasoning used to derive the 2-D parallel axis theorem. We cannot do the last line in that derivation, however, since rcm/O2=xcm/O2+ycm/O2+zcm/O2xcm/O2+ycm/O2, because now, for three-dimensional objects, zcm/O0.

Now, let’s look at an off-diagonal term.

IxyO = x/Oy/O𝑑m
= (xcm/O+x/cmx/O)(ycm/O+y/cmy/O)dm
= xcm/Oycm/O𝑑mmycm/Ox/cm𝑑m0
xcm/Oy/cm𝑑m0x/cmy/cm𝑑mIxycm
= mxcm/Oycm/O+Ixycm

Similarly, we can calculate the other terms to get the whole 3-D parallel axis theorem.

[𝑰O]=[𝑰cm]+m[ycm/o2+zcm/o2xcm/oycm/oxcm/ozcm/oxcm/oycm/oxcm/o2+zcm/o2ycm/ozcm/oxcm/ozcm/oycm/ozcm/oxcm/o2+ycm/o2].

Again, one can think of this result as follows. The moment of inertia matrix about point O is the same as that for parallel axes through the center-of-mass plus the moment of inertia matrix for a point mass at the center-of-mass.

Relation between 2-D and 3-D parallel axis theorems

The (3,3) (lower right corner) element in the matrix of the 3-D parallel axis theorem is the 2-D parallel axis theorem.

SAMPLE 20.9

Filename:tfigure8-alt-app2c
Figure 20.64:

For the dumbbell shown in Figure 20.64, take m=0.5kg and =0.4m. Given that at the instant shown θ=30 and the dumbbell is in the yz-plane, find the moment of inertia matrix [𝑰O], where O is the midpoint of the dumbbell.

Solution The dumbbell is made up of two point masses. Therefore we can calculate [𝑰O] for each mass using the formula from the table on the inside back cover of the text and then adding the two matrices to get [𝑰O] for the dumbbell.
Now, from Table 4.9 of the text,

[𝑰O]=m[y2+z2xyxzxyx2+z2yzxzyzx2+y2]

For mass 1 (shown in Figure 20.65)

x=0,y=2cosθ,z=2sinθ.

Therefore,

Filename:sfig4-6-4a
Figure 20.65:
[𝑰O]mass1 = m[2400024sin2θ24cosθsinθ024cosθsinθ24cos2θ]
= 0.5kg[0.04m20000.04m2140.04m23400.04m2340.04m234]
= 0.02kgm2[1000143403414]

Similarly for mass 2,

x=0,y=2cosθ,z=2sinθ
  [𝑰O]mass2=0.02kgm2[1000143403414].

Therefore,

[𝑰O] = [𝑰O]mass1+[𝑰O]mass2
= 0.04kgm2[1000143403414].

SAMPLE 20.10

Filename:sfig4-6-3
Figure 20.66:

A uniform rod of mass m=2kg and length =12m is pivoted at one of its ends. At the instant shown, the rod is in the xy-plane and makes an angle θ=45 with the x-axis. Find the moment of inertia matrix [𝑰O] for the rod.

Solution The moment of inertia matrix [𝑰O] for a continuous system is given by

[𝑰O]=overallmass[y2+z2xyxzxyx2+z2yzxzyzx2+y2]𝑑m.
Filename:sfig4-6-3a
Figure 20.67:

Thus we need to carry out the integrals for the rod to find each component of the inertia matrix [𝑰O]. Let us consider an infinitesimal length element dl of the rod at distance l from O (see Fig 20.67). The mass of this element is dm=mmass/lengthdl, where m is the total mass.
The coordinates of this element are

x=lcosθ,y=lsinθ,z=0(since the rod is in the xy-plane).

Therefore,

Ixx = m(y2+z20)𝑑m=0l2sin2θy2mdldm
= msin2θ0l2𝑑l=msin2θ33=m23sin2θ.

Similarly,

Iyy = m(x2+z20)𝑑m=0l2cos2θmdl=m23cos2θ.
Izz = m(x2+y2)𝑑m=0l2m𝑑l=m23.
Ixy = mxy𝑑m=0l2cosθsinθmdl=m23cosθsinθ.
Ixz = mxz0dm=0.Iyz=myz0dm=0.

Thus,

[𝑰O]=m23[sin2θsinθcosθ0sinθcosθcos2θ0001].

Substituting the values of m, and θ, we get

[𝑰O]=0.083kgm2[110110001].

Problems for 20.3 Moment of inertia matrices: [𝑰cm] and [𝑰O]

20.3.1  Find the moment of inertia matrix [𝑰O] of the system shown in the figure, given that m=0.25kg and =0.5m. Which components in the matrix [𝑰O] will change if one mass is removed and the other one is doubled?

Filename:pfigure4-4-rp12
Figure 20.68

20.3.2  A dumbbell (two masses connected by a massless rod of length ) lies in the xy plane. Compute the moments and products of inertia of this body about the x, y, and z axes. What are the principle axes of this body?

Filename:pfigure-blue-38-2
Figure 20.69

20.3.3   Equivalent dumbbell. It is a great fact that, for the purposes of dynamics, every rigid body can be replaced by a massless rigid structure holding six masses. The structure looks like a child’s jack. It has three orthogonal bars with three pairs of point masses along the bars. The center of the jack is at the center of mass. The bars are in the directions of the eigenvectors of the moment of inertia matrix of the body (the principle directions). For flat objects one only needs four point masses. Where should one put four masses on a massless disk to make it equivalent to a uniform and ‘massful’ disk in problem 20.3(a). That is, what set of four masses has the same total mass, the same center of mass location, and the same moment of inertia matrix as a uniform disk?

20.3.4  Three identical solid spheres of mass m and radius r are situated at (b,0,0),(0,b,b), and (b,0,2b) relative to the origin O. Find all nine components of [𝑰O] of the system of three spheres.

Hint: [𝑰cm] for each sphere is

[25mr200025mr200025mr2].

Use the parallel axis theorem.

Filename:pfigure-blue-49-2
Figure 20.70

20.3.5  Find the moment of inertia matrix [𝑰O] for the ‘L’-shaped rod shown in the figure. The mass of the rod is m=1.5kg and each leg is =0.4m long. You may ignore the width of the rod.

Filename:pfigure4-4-rp13
Figure 20.71

20.3.6  Moment of Inertia Matrix. This problem concerns a uniform flat disk of radius R0 with mass Md=Mdisc. In the two parts of this problem, you are to calculate the moment of inertia matrix for the disk. In the first part, do it by direct integration and, in the second part, by using a change of basis.

  1. (a)

    Using a coordinate system that is lined up with the disk so that the z axis is normal to the disk and the origin of the coordinate system is at the center of mass (the center) of the disk, calculate the moment of inertia matrix for the disk. [Hint: use polar coordinates, and use dm=(M/πRO2)dA, and dA=RdRdθ.]

  2. (b)

    Use a coordinate system xyz aligned as follows: the origin of the coordinate system is at the center of the disk, the z axis makes an angle ϕ with the z axis which is normal to the plane of the disk (measured from the z axis towards the z axis), the disk intersects the xy plane along the y axis. The y axis and the y axis are coincident. Calculate the moment of inertia matrix for the disk. Alternate method: use the solution from (a) above with a change of basis.

Filename:f92h7p1
Figure 20.72

20.3.7  For each of the shapes below, mark the approximate location of the center of mass and the principle axes about the center of mass

Filename:pfigure-blue-52-2
Figure 20.73

20.3.8  Find the moment of inertia matrix [𝑰cm] for the uniform solid cylinder of length and radius r shown in the figure.

Filename:pfigure4-4-rp16
Figure 20.74

20.3.9  An ‘L’-shaped structure consists of two thin uniform rectangular plates each of width w=200mm. The vertical plate has height h=400mm and the horizontal plate has depth d=300mm. The plates are welded at right angles along the side measuring 200mm. The total mass of the structure is m=0.7kg. Find the moment of inertia matrix [𝑰O].

Filename:pfigure4-4-rp14
Figure 20.75

20.3.10  A solar panel consists of three rectangular plates each of mass m=1.2kg, length =1 and width w=0.5m The two side panels are inclined at 30 with the plane of the middle panel. Evaluate all components of the moment of inertia matrix [𝑰O] of the solar panel.

Filename:pfigure4-4-rp15
Figure 20.76

20.4 Mechanics using the moment of inertia matrix

Once one knows the velocity and acceleration of all points in a system one can find all of the motion quantities in the equations of motion by adding or integrating using the defining sums from the inside cover. This addition or integration is an impractical task for many motions of many objects where the required sums may involve billions and billions of atoms or a difficult integral. Linear momentum and the rate of change of linear momentum can be calculated by just keeping track of the center-of-mass of the system of interest. One would like something so simple for the calculation of angular momentum.

We are in luck if we are only interested in the two-dimensional motion of two-dimensional rigid bodies, the scalar moment of inertia from Chapter 7 is all we need. The luck is not so great for 3-D rigid bodies but still there is some simplificationmargin: For general motion of non-rigid bodies, a topic not covered in this book, the simplification for linear momentum still holds; linear momentum and its rate of change are given by the system mass times the velocity and acceleration of the center-of-mass. But there is no general simplification for the sums needed to evaluate angular momentum and energy or their rates of change. .The simplification is to use the moment of inertia matrix [𝑰cm] from the previous section. One may have to do a sum or integral to find IIzzcm or [𝑰cm] if an adequate table is not handy. But once you know [𝑰cm], say, you need not work with the integrals to evaluate angular momentum and its rate of change. Assuming that you are comfortable calculating and looking-up moments of inertia, we proceed to use it for the purposes of studying mechanics.

Let’s now consider again the conically swinging rod of fig. 20.30 on page 20.30. Method 1 for evaluating 𝑯˙/O was evaluating the sums directly. If we accept the formulae presented for rigid bodies in Table I at the back of the book, we can find all of the motion quantities by setting 𝝎=ω𝒌ˆ and 𝜶=𝟎.

Method 2 of evaluating 𝑯˙/O: using the xyz moment of inertia matrix about point O

In section 20.1 we examined the conical swinging of a straight uniform rod. Now let’s look at that example again using its moment of inertia matrix. For a rigid body in constant rate circular motion about an axis through O,

𝑯˙/O=𝝎×𝑯/O

because the 𝑯/O vector rotates with the body. For a rigid body rotating about point O,

𝑯/O=[𝑰O]𝝎.

We assumed at the outset that

𝝎=ω𝒌ˆ=[00ω].

So, the only trick is to find [𝑰O] for a rod in the configuration shown. We recall, look up, or believe for now that

[𝑰O]=[IxxOIxyOIxzOIyxOIyyOIyzOIzxOIzyOIzzO].
IxyO = x0y𝑑m=0 (x=0, all mass in the yz plane)
IxzO = x0z𝑑m=0 (x=0, all mass in the yz plane)
IyyO = (x20+z2)𝑑m (x=0, all mass in the yz plane)
= 0(scosϕ)2zρ𝑑s
= cos2ϕρ33=cos2ϕm23
IzzO = (x20+y2)𝑑m (x=0, all mass in the yz plane)
= (ssinϕ)2yρ𝑑s
= sin2ϕρ33=sin2ϕm23
IxxO = (y2+z2)𝑑m (both integrals have been evaluated above)
= IyyO+IzzO (perpendicular axis theorem)
= (cos2ϕ+sin2ϕ)m23
= m23
IyzO = yz𝑑m
= (ssinϕ)(scosϕ)ρ𝑑s
= sinϕcosϕρ33=sinϕcosϕm23.

Putting these terms all together in the matrix, we get

[𝑰O]=m23[1000cos2ϕcosϕsinϕ0cosϕsinϕsin2ϕ].

Now we can calculate 𝑯/O as

𝑯/O = [𝑰O]𝝎
= m23[1000cos2ϕcosϕsinϕ0cosϕsinϕsin2ϕ][𝑰O][00ω]𝝎
= m23[0ωcosϕsinϕωsin2ϕ]
= mω23sinϕ(cosϕȷˆ+sinϕ𝒌ˆ).

You may notice, by the way, that for this problem, 𝑯/O (in the direction of cosϕȷˆ+sinϕ𝒌ˆ) is perpendicular to the rod (in the direction of sinϕȷˆcosϕ𝒌ˆ). So 𝑯/O is not in the direction of 𝝎. margin: 𝑯/O is only parallel to 𝝎 if rotation is about one of the principal axes (eigen-vector directions) of the inertia matrix.

Now we calculate 𝑯˙/O as

𝑯˙/O = 𝝎×𝑯/O
= (ω𝒌ˆ)×[ωm23(cosϕsinϕȷˆ+sin2ϕ𝒌ˆ)]
= m23cosϕsinϕω2ıˆ,

the same result we got before.

Method 3 of calculating 𝑯˙/O: using [𝑰O] and a coordinate system lined up with the rod

Filename:tfigure4-spherical-rotaxis
Figure 20.77: The spherical pendulum using xyz axes aligned with the rod.

Here, as suggested above, we redo the problem using a rotated set of axes better aligned with the rod. This method makes calculation of the moment of inertia matrix quite a bit easier — we can even look it up in a table — but makes the determination of 𝝎 a little harder. We are stuck finding the components of 𝝎 in a rotated coordinate system.

Referring to the table of moment of inertias on the inside of the back cover and taking care because different coordinates are used, the moment of inertia matrix about point O for a thin uniform rod, in terms of the rotated coordinates is

[𝑰O]xyz=m23[100000001].

and the angular velocity in rotated coordinates is 𝝎=ωcosϕȷˆ+ωsinϕ𝒌ˆ, which can be written in component form as

[𝝎]xyz=ω[0cosϕsinϕ].

We calculate the angular momentum about point O as

𝑯/O=[𝑰O]𝝎vector equation, independent of coordinates
[H/OxH/OyH/Oz]=m23[100000001]xyzω[0cosϕsinϕ]xyzall components in xyz coordinate system
[H/OxH/OyH/Oz]=m23ω[00sinϕ]xyz.

Finally, we calculate 𝑯˙/O as

𝑯˙/O = 𝝎×𝑯/O
= [ω(cosϕȷˆ+sinϕ𝒌ˆ]×(ωm23sinϕ𝒌ˆ)
= sinϕcosϕmω223ıˆ.

This answer is again the same as what we got before because the x and x axis are coincident so ıˆ=ıˆ.

The multitude of ways to calculate 𝑯˙/O

We just showed three ways to calculate 𝑯˙/O for a conically swinging stick, but there are many more. You can get a sense of the possibilities by studying table I summarizing momenta and energy in the back of the book. Here are three basic choices.

(1) 𝑯˙/O = 𝒓/O×𝒂𝑑m,
(2) 𝑯˙/O = 𝝎×𝑯/O,or
(3) 𝑯˙/O = 𝒓cm/O×𝒂cmmtot+𝑯˙cm.

Choice (3) is the safest choice to make if you are in doubt, since it is the only one of the three choices that does not depend on point O being a fixed point (which it was for this example). For option (2) above, we can calculate 𝑯/O various ways as

(a) 𝑯/O = 𝒓/O×𝒗𝑑mor,
(b) 𝑯/O = [𝑰O]𝝎or,
(c) 𝑯/O = 𝒓cm/O×𝒗cmmtot+𝑯cm.

For option (3) above, we can calculate 𝑯˙cm as

(d) 𝑯˙cm = 𝒓/cm×𝒂/cm𝑑m or,
(e) 𝑯˙cm = 𝝎×𝑯cm.

For (c) and (e) above, we can calculate 𝑯cm as

(f) 𝑯cm = 𝒓/cm×𝒗/cm𝑑m or,
(g) 𝑯cm = [𝑰cm]𝝎.

For either of options (b) or (g), we can calculate [I] relative to the xyz axes (as we did for the second method on the previous pages) or relative to some rotated axes better aligned with the rod (as we did for the third method on the previous pages.)

As you can surmise from all the choices above, the list of options for the calculation of 𝑯˙/O is too long and boring to show here.

The pros and cons of these methods depend on the problem at hand. If you want to avoid integration you are pretty much stuck using either [𝑰O] or [𝑰cm] with (e) and (g) above. Avoiding the integrals depends on your having a table of moments of inertia (like table IV in the back of this book).

Energy of things going in circles at variable rate

The energy only depends on the speeds of the parts of a system, not their accelerations. So, as with constant rate motion,

EK = 12v2𝒗𝒗𝑑m
= 12(𝝎×𝒓)(𝝎×𝒓)𝑑m
= 12ω2RR is the distance from the axis to the mass2𝑑m
= 12𝝎[𝑰O]𝝎
= 12ω2IzzoFor rotation about the z axis the 9 terms in the matrix formula reduce to this one simple term.

The safest bet

The following are the most reliable (least prone to error) formulas for evaluating the motion quantities for a rigid body rotating about a fixed axes (they also apply to arbitrary motion of a rigid body).

𝑯/O = 𝒓cm/o×mtot𝒗cm+[𝑰O]𝝎 (20.47)
𝑯˙/O = 𝒓cm/o×mtot𝒂cm+𝝎×[𝑰cm]𝝎𝑯cm+[𝑰cm]𝝎˙ (20.48)
EK = 12mtotvcm2+𝝎([𝑰cm]𝝎) (20.49)
Filename:tfigure5-7
Figure 20.78: The two terms in the rate of change of angular momentum are shown.

Please survey table I at the back of the book which summarizes the ways of evaluating momenta and energy.

The geometry of 𝑯˙/O

It is possible to understand the formula for the change of angular momentum geometrically. Here is one way of looking at it. For this example it is most direct to look at 𝑯˙/O calculated using [𝑰O], but the same discussion would work with 𝑯˙cm calculated using [𝑰cm].

𝑯˙/O=𝝎×𝑯/OContribution from rotation of 𝑯/O+[𝑰O]𝝎˙Due to the changing length of 𝑯/O

The first term in the equation describes the rotation of the angular momentum vector. The second term describes its rate of change of length.

Since the body is spinning about a fixed axis, the orientation of the axis of rotation is not only fixed relative to the Newtonian frame, the room environment, but also relative to the body. At one instant of time we use a coordinate system to calculate 𝑯/O. After the body has rotated a little, if we now used a coordinate system that rotated the same amount as the body, a coordinate system ‘glued’ to the body, we could calculate 𝑯/O again.

This new calculation of 𝑯/O will be almost identical to the calculation before the small rotation, however, because the moment of inertia matrix does not change in time relative to a coordinate system that moves with the body. Also, the coordinates of 𝝎 will be unchanged except for possibly a multiplication by a constant because the direction of 𝝎 doesn’t change. Thus the only change of 𝑯/O as represented in this rotated coordinate system, is a possible change in its length due to a change in the spinning rate.

But this new coordinate system is a rotated coordinate system. To find the actual change in 𝑯/O we need to take this rotation into account. To picture this rotation consider the special case when the rotation rate is constant. Then the vector 𝑯/O is constant in the rotating coordinate system. That is, the vector 𝑯/O rotates with the body.

So the net change in 𝑯/O is a change due to rotation of the body added to a change due to the change in the rotation rate. For small angular changes, the direction of the first term is the same as the tangent to the circle that is traced by the tip of the angular momentum vector 𝑯/O as drawn on the body. To approximate the first term, consider the following reasoning. First, let’s denote the first term by (𝑯˙/O)rot=𝝎×𝑯/O, where

the subscript ‘rot’ indicates that this term is the contribution to 𝑯˙/O from the rotation of 𝑯/O. For small angular changes, (𝑯˙/O)rot=(d𝑯/O/dt)rot=𝝎×𝑯/O(Δ𝑯/O)rot/Δt. Thus, the term due to rotation is approximately, for small Δt, (Δ𝑯/O)rot=Δt(𝝎×𝑯0).

The contribution to 𝑯˙/O due to change of length of 𝑯/O is [𝑰O]𝝎˙. Similarly, this term is approximately, for small Δt, Δt([𝑰O]𝝎˙).

Rotation of a rigid body about a fixed axis: the general case

Consider a general rigid body of mass m and moment of inertia matrix with respect to the center-of-mass [𝑰cm] rotating about a fixed axis. Without loss of generality, let the axis of rotation be the 𝒌ˆ axis.

Filename:tfigure5-gen-rigid-body
Figure 20.79: Free-body diagram of a general rigid body rotating about the z-axis.

Linear momentum balance

Referring to the free-body diagram of the rigid body, linear momentum balance gives

𝑭 = 𝑳˙
𝑭net = m𝒂cm (20.50)
= m[ω𝒌ˆ×(ω𝒌ˆ×𝒓cm/O)+ω˙𝒌ˆ×𝒓cm/O].

The first term on the right hand side of equation 20.50, the centripetal term, is directed from the center-of-mass (

Filename:center-of-mass
) through the axis of rotation; that is, it lies in the xy-plane (or has no 𝒌ˆ component). The second term on the right hand side of equation 20.50, the tangential term, is normal to the plane determined by the axis and the center-of-mass. It is tangent to the circle that the center-of-mass travels on. It is zero if the center of mass is on the axis.

Angular momentum balance

Angular momentum balance about point O gives

𝑴O = 𝑯˙/O
𝑴net = 𝒓cm/O×m𝒂cmterm (i) (20.51)
+ω𝒌ˆ×{[𝑰cm]ω𝒌ˆ}term (ii)
+[𝑰cm](ω˙𝒌ˆ)term (iii).

Let’s look at each of the three terms on the right hand side in turn.

term (i)

The first term (i) on the right hand side of equation 20.51 is

term (i) = 𝒓cm/O×m𝒂cm𝒂cm=ω𝒌ˆ×(ω𝒌ˆ×𝒓cm/O)  +ω˙𝒌ˆ×𝒓cm/O
= m[𝒓cm/O×(ω𝒌ˆ×(ω𝒌ˆ×𝒓cm/O))
+𝒓cm/O×(ω˙𝒌ˆ×𝒓cm/O)]

Now, let’s consider the two parts of term (i) in turn.

Filename:tfigure5-term1-a
Figure 20.80: The first part of term (i)

The first part of term (i) is in the direction of 𝒌ˆ×𝒓cm/O. For example, if the center-of-mass is in the xz plane, this contribution to 𝑴net is in the ȷˆ direction and could be accommodated by reaction forces on the axis in the ıˆ and ıˆ directions.

Filename:tfigure5-term1-b
Figure 20.81: The second part of term (i)

Now, the second part of term (i) can be decomposed into a part along the 𝒌ˆ direction and a part perpendicular to 𝒌ˆ (along d). The part along 𝒌ˆ is

𝒓cm/O×(ω˙𝒌ˆ×𝒓cm/O)𝒌ˆ=md2ω˙𝒌ˆ.

The part of term (i) perpendicular to 𝒌ˆ has magnitude mω˙dw.

term (ii)

Now, let’s look at the second term in equation 20.51, term (ii), and expand it.

term (ii) = ω𝒌ˆ×[𝑰cm]ω𝒌ˆ
= ω𝒌ˆ×[IxxcmIxycmIxzcmIxycmIyycmIyzcmIxzcmIyzcmIzzcm][00ω]
= ω𝒌ˆ×[Ixzcmωıˆ+Iyzcmωȷˆ+Izzcmω𝒌ˆ]
= Ixzcmω2ȷˆIyzcmω2ıˆ

So, the second term in equation 20.51, term (ii), has no 𝒌ˆ component. The ıˆ and ȷˆ components are due to the off-diagonal terms in [𝑰cm]. The off-diagonal terms are responsible for dynamic imbalance.

term (iii)

Now, the third term in equation 20.51, term (iii), is

[𝑰cm]ω˙𝒌ˆ = [IxxcmIxycmIxzcmIxycmIyycmIyzcmIxzcmIyzcmIzzcm][00ω˙]
= Ixzcmω˙ıˆ+Iyzcmω˙ȷˆ+Izzcmω˙𝒌ˆ.

Putting terms (i), (ii), and (iii) back together

Now, let’s put the three terms back together into the equation of angular momentum balance about point O, equation 20.51. To help with the interpretation, assume the center-of-mass is in the xz plane. So, we start with

𝑴net=𝑯˙/O.

Breaking this vector equation into components, we get

Mxx=𝑴netıˆ = 𝑯˙/Oıˆ
Myy=𝑴netȷˆ = 𝑯˙/Oȷˆ
Mzz=𝑴net𝒌ˆ = 𝑯˙/O𝒌ˆ.

Now, using all the results so far, we find the torques about the x, y, and z axes. For the z-axis,

Mz=[md2d = distance of cm from axis of rotation+Izzcm]sometimes called ‘[𝑰O]’ ω˙. (20.52)

The only torque about the z-axis, the axis of rotation, is due to the acceleration about that axis.

Next, for the x-axis,

Mx=Iyzcmω2+Ixzcmω˙+mω˙dw. (20.53)

There is a torque about the x-axis due to off-diagonal terms in [𝑰cm]. One term is the dynamic imbalance term Iyzcm and the other is associated with angular acceleration. (Recall, the center-of-mass is on the xz plane.) There is also a term due to the center of mass tangential acceleration since the center-of-mass is on a plane that does not contain point O.

Finally, the torque about the y-axis is

My=Ixzcmω2+Iyzcmω˙+(mω2dw). (20.54)

Here, we have nearly the same types of terms as in equation 20.53. In this case, though, the centripetal acceleration causes the center-of-mass motion to contribute to the torque.

So, if

𝝎=ω𝒌ˆ and 𝝎˙=ω˙𝒌ˆ

then

𝑴O=𝑯˙/O = (Iyzcmω2+(Ixzcm+mdw)ω˙)ıˆ
+((Ixzcmmdw)ω2+Iyzcmω˙)ȷˆ
+(md2+Izzcm)ω˙𝒌ˆ.

SAMPLE 20.11  A scalar times a vector is a vector. A matrix times a vector is a vector. What is the difference? Find the vectors 𝑯1=IG𝝎 and 𝑯2=[IG]𝝎, if 𝝎=2rad/sıˆ+3rad/sȷˆ,IG=10kgm2 and [IG]=[5000520210]kgm2. Draw 𝝎,𝑯1 and 𝑯2.

Solution

Filename:sfig1-2-12
Figure 20.82:
𝑯1 = IG𝝎
= 10kgm2(2rad/sıˆ+3rad/sȷˆ)
= (20ıˆ+30ȷˆ)kgm2/s.
𝑯2 = [IG]𝝎
= [5000520210]kgm2{230}rad/s
= (10ıˆ+15ȷˆ6𝒌ˆ)kgm2/s.

These two vectors, 𝑯1 and 𝑯2, are shown along with 𝝎 in Fig. 20.82. Note that 𝑯1 has the same direction as 𝝎 but 𝑯2 does not. 𝑯1 and 𝝎 are both in the xy-plane but 𝑯2 is not; it is in 3-D.

Comments: Multiplying a vector by a scalar does not change the direction of the vector but multiplying by a matrix does change the direction, in general. Find the angles between 𝝎 and 𝑯1 and between 𝝎 and 𝑯2 to convince yourself.

SAMPLE 20.12

Filename:sfig4-6-5
Figure 20.83:

The composite rod OABCD shown in Figure 20.83 is made up of three identical rods OA, BC, CD of mass 0.5kg and length 20cm each, and the rod AB which is half of rod CD. The composite rod goes in circles about the y-axis at a constant rate ω=5rad/s. Find the angular momentum and the rate of change of angular momentum of the rod at the instant shown (i.e., when the rod is in the xy-plane).

  1. 1.

    Are all the components of [𝑰O] necessary to compute 𝑯/O and 𝑯˙/O? Find [𝑰O] or the necessary components of [𝑰O].

  2. 2.

    Find the angular momentum 𝑯/O.

  3. 3.

    Find the rate of change of angular momentum 𝑯˙/O.

  4. 4.

    If the rod were rotating about the z-axis instead, (i.e., if the motion were in the xy-plane) which components of [𝑰O] would be required to find 𝑯/O? What would be the value of 𝑯˙/O in that case?

Solution Since the rod rotates about the y-axis, and O is a fixed point on this axis,

𝝎=ωȷˆand𝑯/O=[𝑰O]𝝎.
  1. 1.

    Since

    𝑯/O = [IxxOIxyOIxzOIxyOIyyOIyzOIxzOIyzOIzzO]{0ω0} (20.62)
    = IxyOωıˆ+IyyOωȷˆ+IyzOω𝒌ˆ (20.63)
    = (IxyOıˆ+IyyOȷˆ+IyzO𝒌ˆ)ω,

    we only need to find three components of [𝑰O] to compute 𝑯/O, namely IxyO, IyyO and IyzO.
    We can compute these components by considering each rod individually. For any rod, the components can be calculated using the values of the components about the center-of-mass of the rod (see table IV in the back of the book) and then using the parallel axis theorem.

    Filename:sfig4-6-5a
    Figure 20.84:
    Rod OA:IyyO = Iyycm+ml24=112ml2+14ml2
    IxyO = Ixycm0+m(xcm/Oycm/O0)=0
    IyzO = Iyzcm0+m(ycm/O0zcm/O0)=0.
    Filename:sfig4-6-5b
    Figure 20.85:
    Rod AB:IyyO = Iyycm0+m2(xcm/O2+zcm/O20)=ml22
    IxyO = Ixycm0+m2(xcm/Olycm/Ol4)=ml28
    IyzO = IyzO0+m(ycm/Ozcm/O0)=0.

    Filename:sfig4-6-5c
    Figure 20.86:
    Rod BC:IyyO = Iyycm+m(xcm/O2+zcm/O20)=ml212+m9l24=73ml2
    IxyO = Ixycm0+m(xcm/Oycm/O)=m3l2l2=3ml24
    IyzO = 0.
    Filename:sfig4-6-5d
    Figure 20.87:
    Rod CD:IyyO = Iyycm0+m(xcm/O2+zcm/O20)=m(2l)2=4ml2
    IxyO = 0
    IyzO = 0.

    Thus for the entire rod,

    IyyO = 13ml2+12ml2+73ml2+4ml2=436ml2
    IxyO = 18ml234ml2=78ml2
    IyzO = 0.

    Answer: IyyO=436ml2,IxyO=78ml2,IyzO=0.

  2. 2.
    𝑯/O = (IyyOıˆ+IxyOȷˆ+IyzO𝒌ˆ)ω(from Eqn 20.63)
    = ml2ω(436ȷˆ78ıˆ)
    = 0.5kg(0.2m)25rad/s(7.17ȷˆ0.87)ıˆ)
    = (0.717ȷˆ0.087ıˆ)kgm2/s.

    Answer: 𝑯/O=(0.717ȷˆ0.087ıˆ)kgm2/s.

  3. 3.
    𝑯˙/O = 𝝎×𝑯/O=ωȷˆ×([𝑰O]𝝎)
    = 5rad/sȷˆ×(0.717ȷˆ0.087ıˆ)kgm2/s
    = 0.435Nm𝒌ˆ.

    Answer: 𝑯˙/O=0.435Nm𝒌ˆ.

  4. 4.

    If the rod were rotating in the xy-plane, it would be planar circular motion. The only component of [𝑰O] required for the calculation of 𝑯/O=IzzO𝝎 will be IzzO. Also,

    𝑯˙/O=𝝎×𝑯/O=𝝎×IzzO𝝎=0

    since 𝝎 is parallel to IzzO𝝎.

SAMPLE 20.13

Filename:sfig4-5-6
Figure 20.88:

A 0.5 kg uniform rectangular solar panel rotates about an off-centered axis zz with constant angular speed ω=1.5rad/s. Axis zz is parallel to the transverse z-axis of the plate and is 1.6m away from the center-of-mass of the plate. The in-plane moments of inertia Ixxcm and Iyycm of the plate are given: Ixxcm=4.0kgm2 and Iyycm=12.4kgm2. At the instant shown in Fig 20.88, calculate

  1. 1.

    the linear momentum of the panel,

  2. 2.

    the angular momentum of the panel, and

  3. 3.

    the kinetic energy of the panel.

Solution Since the plate rotates about the zz-axis which is parallel to the z-axis,

𝝎=ω𝒌ˆ=1.5rad/s𝒌ˆ.
  1. 1.

    Linear momentum:

    𝑳 = mtot𝒗cm=m(𝝎×𝒓cm/O)
    = 0.5kg(1.5rad/s𝒌ˆ×1.6mıˆ)
    = 1.2kgm/sȷˆ

    Answer: 𝑳=1.2kgm/sȷˆ

    Filename:sfig4-5-6a
    Figure 20.89:

    We can easily check the direction of 𝑳 since 𝑳=m𝒗cm, it has to be in the same direction as 𝒗cm. The center-of-mass goes in circles about O, therefore, 𝒗cm is tangential to the circular path, i.e. in the y-direction (see Fig 20.89).

  2. 2.

    Angular momentum:

    𝑯=IzzOω𝒌ˆ

    Thus to find 𝑯, we need to find IzzO. Since zzzz, we can use the parallel axis theorem to find IzzO if we know Izzcm. We are given the in-plane moments of inertia Ixxcm and Iyycm. Therefore, from the perpendicular axis theorem:

    Izzcm=Ixxcm+Iyycm=(4.0+12.4)kgm2=16.4kgm2.

    Now using the parallel axis theorem,

    IzzO = Izzcm+Mrcm/O2
    = 16.4kgm2+0.5kg(1.6m)2=17.68kgm2.
    Therefore,𝑯 = 17.68kgm2(1.5rad/s𝒌ˆ)
    = 26.52kgm2/s𝒌ˆ.

    Answer: 𝑯=26.52kgm2/s𝒌ˆ.

  3. 3.

    Kinetic energy:

    EK = 12IzzOω2=12(17.68kgm2)(1.5rad/s)2
    = 19.89kgm2s2==19.89Nm
    = 19.89J

    Answer: EK=19.89J

SAMPLE 20.14  The rotating crooked plate again. For the crooked plate considered in Sample 20.37 and shown again in Fig. 20.90,

Filename:sfig4-6-8
Figure 20.90: A rectangular plate of mass m rotates with shaft AB at a constant speed 𝝎=ω𝒌ˆ. A coordinate system xyz is aligned with the principle axes of the plate.
  1. 1.

    compute the moment of inertia matrix [𝑰] in the xyz coordinate system,

  2. 2.

    compute the rate of change of angular momentum 𝑯˙/O using the moment of inertia [𝑰]xyz computed above,

  3. 3.

    express 𝑯˙/O as a vector in the xyz coordinate system.


Solution A line sketch of the plate and the two coordinate systems attached to it are shown in Fig. 20.91. The set of basis vectors (ıˆ,ȷˆ,𝒌ˆ) and (ıˆ,ȷˆ,𝒌ˆ) associated with coordinate systems xyz and xyz respectively are shown separately for the sake of clarity. From the diagram of the two basis vector sets, we may write

Filename:sfig4-6-8a
Figure 20.91: A line sketch of the rotating plate along with the two coordinate systems xyz and xyz. The basis vectors associated with the two systems are also shown.
ıˆ=cosϕıˆsinϕ𝒌ˆȷˆ=ȷˆ𝒌ˆ=sinϕıˆ+cosϕ𝒌ˆ}. (20.64)
  1. 1.

    Calculation of [I]xyz:

    [𝑰]xyz=[IxxIxyIxzIxyIyyIyzIxzIyzIzz].

    Let us first consider the off diagonal terms of [𝑰]xyz.

    Since x=0 on the entire plate (the plate is in the yz plane and the origin is on the plate.),

    Ixy=mx0y𝑑m=0,Ixz=mx0z𝑑m=0.
    Filename:sfig4-6-8b
    Figure 20.92: Symmetry about the z axis implies that Iyz=0. Same result holds if we consider the symmetry about the y axis.

    How about Iyz? Well, you can calculate it two ways: (a) carry out the integration Iyz=myz𝑑m over the entire plate mass and find that Iyz=0, or (b) realize that for every mass element dm (=mbdA) with coordinates (+yz), there exists another element dm at (yz) such that the sum of their contributions to the integral is zero. Therefore, Iyz=myz𝑑m=0 margin: You can use a similar argument to find the off-diagonal terms in [𝑰] whenever there is such symmetry with respect to the coordinate axes. Now the other terms:

    Izz = m(x20+y2)𝑑m=/2/2b/2b/2y2mb𝑑y𝑑z
    = mb/2/2(y33|b/2b/2)𝑑z=mbb312/2/2𝑑z
    = mbb312=mb212,

    and similarly,

    Iyy = m(x20+z2)𝑑m=m212,
    Ixx = m(y2+z2)𝑑m=Iyy+Izz=m12(2+b2).

    Thus,

    [𝑰]xyz=m12[2+b20002000b2].
  2. 2.

    Calculation of H˙/O:

    𝑯˙/O=𝝎×𝑯/O and 𝑯/O=[𝑰]xyz{𝝎}xyz.

    The subscripts xyz in [𝑰] and 𝝎 have been used to denote that both [𝑰] and 𝝎 are expressed in xyz coordinate system. margin: We can express [𝑰] and 𝝎 in any coordinate system of our choice but both of them must be in the same system for their product to be valid. We have calculated [I]xyz above. Now we need to find {𝝎}xyz.

    Let 𝝎=ωxıˆ+ωyȷˆ+ωz𝒌ˆ. But, we can also write, 𝝎=ω𝒌ˆ. So,

    ω𝒌ˆ=ωxıˆ+ωyȷˆ+ωz𝒌ˆ (20.65)

    Dotting both sides of Eqn. (20.65) with ıˆ, ȷˆ, and 𝒌ˆ and using the relationships in (20.64) we get

    ωx=ω(𝒌ˆıˆsinϕ)=ωsinϕ,ωy=ω(𝒌ˆȷˆ0)=0,ωk=ω(𝒌ˆ𝒌ˆcosϕ)=ωcosϕ.

    Thus,

    {𝝎}xyz=ωsinϕıˆ+ωcosϕ𝒌ˆ

    So,

    𝑯/O = m12[2+b20002000b2]{ωsinϕ0ωcosϕ}
    = {m12(2+b2)ωsinϕ0m12b2ωcosϕ}
    or 𝑯/O = mω12[(2+b2)sinϕıˆ+b2cosϕ𝒌ˆ].

    Now we can easily compute 𝑯˙/O as follows:

    𝑯˙/O = 𝝎×𝑯/O
    = (ωsinϕıˆ+ωcosϕ𝒌ˆ)×mω12[(2+b2)sinϕıˆ+b2cosϕ𝒌ˆ]
    = mω212b2sinϕcosϕȷˆmω212(2+b2)sinϕcosϕȷˆ
    = mω2122sinϕcosϕȷˆ.
  3. 3.

    Now, back to the xyz coordinate system: Since ȷˆ=ȷˆ [Eqn (20.64)], we have

    𝑯˙/O=mω2122sinϕcosϕȷˆ

    which is the same result as obtained in Sample 20.37.

SAMPLE 20.15   The calculation of H˙ using the moment of inertia matrix. For the crooked plate considered in Sample 20.37, find the rate of change of angular momentum 𝑯˙/O, using the moment of inertia [𝑰O], calculated in the xyz coordinate system

Solution We calculate 𝑯˙/O using the following formula:

𝑯˙/O=[𝑰O]𝝎˙𝟎+𝝎×([𝑰O]𝝎)𝑯/O.

Note that point O is the center-of-mass of the plate. Let us write the expression for the angular momentum 𝑯/O in matrix form.

{Hx/OHy/OHz/O}=[IxxIxyIxzIxyIyyIyzIxzIyzIzz]{00ω}.

It is clear from the components of 𝝎 that we only need the last column of [𝑰] matrix since other elements of [𝑰] will multiply with zeros of 𝝎 vector. Carrying out the multiplication we get

{H0xH0yH0z}={IxzωIyzωIzzω}

or, in vector form

𝑯/O=Ixzωıˆ+Iyzωȷˆ+Izzω𝒌ˆ.

Therefore,

𝑯˙/O = ω𝒌ˆ×ω(Ixzıˆ+Iyzȷˆ+Izz𝒌ˆ)
= ω2Iyzıˆ+ω2Ixzȷˆ.

But,

Ixz=mxz𝑑m and Iyz=myz𝑑m

and x=wsinϕ,y=y,z=wcosϕ (see Fig. 20.40), hence

𝑯˙/O = ω2m(wycosϕıˆ+w2sinϕcosϕȷˆ)𝑑m
= b/2b/2/2/2ω2(w2sinϕcosϕȷˆ+wycosϕıˆ)mb𝑑w𝑑y

which is the same integral as obtained in Sample 20.37 for 𝑯˙/O. Therefore, the result is also the same:

𝑯˙/O=mω2212sinϕcosϕȷˆ.

Answer: H˙/O=mω2212sinϕcosϕȷˆ

SAMPLE 20.16  Direct application of formula: A rectangular plate is mounted on a massless shaft with the center-of-mass of the plate on the shaft axis. The shaft rotates about its axis with angular acceleration 𝝎˙=0.5rpm/s(ıˆ+ȷˆ). At the instant of interest, the angular velocity is 𝝎=100rpm(ıˆ+ȷˆ) and the components of the moment of inertia matrix of the plate are Ixx=2kgm2,Iyy=4kgm2,Izz=6kgm2 and Ixy=Iyz=Ixz=0.

  1. 1.

    Find the angular momentum of the plate about its mass-center and show that it is not in the same direction as the angular velocity.

  2. 2.

    Find the net moment acting on the plate.

Solution We are given the angular velocity, the angular acceleration, and the moment of inertia matrix of the plate:

𝝎 = 100rpm(ıˆ+ȷˆ)=10.47rad/s(ıˆ+ȷˆ)
𝝎˙ = 5rpm/s(ıˆ+ȷˆ)=0.52rad/s2(ıˆ+ȷˆ)
[𝑰cm] = [200040006]kgm2.
  1. 1.

    The angular momentum: The angular momentum of the plate is

    𝑯cm = [𝑰cm]𝝎
    = [200040006]kgm2{10.4710.470}rad/s
    = (20.94ıˆ+41.88ȷˆ)kgm2s1Nms
    = (20.94ıˆ+41.88ȷˆ)Nms.

    Answer: 𝑯cm=(20.94ıˆ+41.88ȷˆ)Nms

    There are many ways of showing that 𝑯cm is not parallel to 𝝎. We can simply draw the two vectors and show that they are not parallel

    We can, alternatively, take the cross product of the two vectors:

    𝑯cm×𝝎 = (20.94ıˆ+41.88ȷˆ)Nms×10.47rad/s(ıˆ+ȷˆ)
    = |ıˆȷˆ𝒌ˆ20.9441.88010.4710.470|Nm
    = 𝒌ˆ(219.24438.48)Nm=20.94Nm𝒌ˆ

    which is not zero, implying that the two vectors are not parallel.

  2. 2.

    The net moment: The net moment on the plate can be found by applying angular momentum balance:

    𝑴/cm = 𝑯˙cm=𝝎×[𝑰cm]𝝎𝑯cm+[𝑰cm]𝝎˙
    = 20.94Nm𝒌ˆ𝝎×𝑯cm+1.04Nmıˆ+2.08Nmȷˆ[𝑰cm]𝝎˙

    Answer: 𝑴/cm=(1.04ıˆ+2.08ȷˆ+20.94𝒌ˆ)Nm

SAMPLE 20.17

Filename:sfig5-5-2
Figure 20.93: A rod, welded to a tipped shaft AB, swings around the shaft axis if it is tipped slightly from the vertical plane yz.

A rod swings around a tipped axis in 3-D. A rigid shaft of negligible mass is connected to frictionless hinges at A and B. The shaft is tipped from the horizontal plane (xy-plane) such that the shaft axis makes an angle γ with the vertical axis. A uniform rod OC of mass m and length is welded to the shaft. At time t=0, the rod is tipped by a small angle ϕ from its position OC in the yz (or yz) plane. Find the equation of motion of the rod.


margin:

Solution Let ıˆ,ȷˆ,𝒌ˆ be the basis vectors associated with the xyz coordinate system. Since yz axes can be obtained by rotating yz axes counterclockwise about the x-axis by an angle 90oγ, we can relate the basis vectors of the two coordinate systems with the help of Fig. 20.94:

Filename:sfig5-5-2a
Figure 20.94: yz axes are obtained by rotating yz axes counterclockwise about the x-axis by an angle 90oγ.
ıˆ=ıˆ,ȷˆ=sinγȷˆ+cosγ𝒌ˆ,𝒌ˆ=cosγȷˆ+sinγ𝒌ˆ. (20.71)

The free-body diagram  of the shaft with rod OC is shown in Fig. 20.95. We can write angular momentum balance for this system about any point on the axis of rotation AB. However, rather than writing angular momentum balance about a point, let us write angular momentum balance about axis AB. This ‘trick’ will eliminate reactions 𝑹A and 𝑹B from our equations. Angular Momentum Balance about axis AB is:

𝝀ˆAB[𝑴O = 𝑯˙/O]
or ȷˆ𝑴O = ȷˆ𝑯˙/O.

Calculation of (ȷˆMO) : Since 𝑹A and 𝑹B pass through axis AB, they do not produce any moment about this axis. Therefore,

Filename:sfig5-5-2b
Figure 20.95: Free-body diagram of the bar and shaft system. In the inset, only the weight of the bar, mg, is shown for clarity of geometry.
ȷˆ𝑴O = ȷˆ[𝒓G/O×mg(𝒌ˆ)]
= ȷˆ[L2(sinϕıˆ+cosϕ𝒌ˆ)×mg(𝒌ˆ)]
= (sinγȷˆ+cosγ𝒌ˆȷˆ)[L2mg(sinϕȷˆ+cosϕcosγıˆusing Eqn.(20.71))]
= 2mgsinϕsinγ

Calculation of ȷˆH˙/O : Since the rod rotates about axis AB, we may write

𝝎 = ϕ˙ȷˆ,
𝝎˙ = ϕ¨ȷˆ.

Now,

𝑯˙/O=[𝑰O]{𝝎˙}+𝝎×[𝑰O]{𝝎}

where

[𝑰O]{𝝎˙} = [𝑰O]xyz{𝝎˙}xyz
= [IxxIxyIxzIxyIyyIyzIxzIyzIzz]{0ϕ¨0}
= {Ixyϕ¨Iyyϕ¨Iyzϕ¨}.
Similarly,
[𝑰O]{𝝎} = {Ixyϕ˙Iyyϕ˙Iyzϕ˙},
 𝝎×[𝑰O]{𝝎} = ϕ˙ȷˆ×[Ixyϕ˙ıˆ+Iyyϕ˙ȷˆ+Iyzϕ˙𝒌ˆ]
= Ixyϕ˙2𝒌ˆ+Iyzϕ˙2ıˆ.

Therefore,

ȷˆ𝑯˙/O = ȷˆ[(Ixyϕ¨ıˆ+Iyyϕ¨ȷˆ+Iyzϕ˙𝒌ˆ)+(Ixyϕ˙2𝒌ˆ+Iyzϕ˙2ıˆ)]
= Iyyϕ¨.

Now, setting ȷˆ𝑴O=ȷˆ𝑯˙/O, we get

L2mgsinϕsinγ = Iyyϕ¨
or ϕ¨ = mg(L/2)sinγIyysinϕ
or ϕ¨ = Csinϕ

where

C=mg(L/2)sinγIyy.

For rod OC,

Iyy=m(x2+z2)𝑑m=0Ll2mL𝑑l=13mL2.

Therefore, the equation of motion of the rod is

ϕ¨mg(L/2)sinγ(1/3)mL2sinϕ = 0
ϕ¨3gsinγ2Lsinϕ = 0.

Answer: ϕ¨3gsinγ2Lsinϕ=0

SAMPLE 20.18

Filename:sfig7-4-2
Figure 20.96:

Kinetic energy in 3-D rotation. A thin rod of mass 2kg and length L=12m is welded to a massless shaft at an angle θ=45. The shaft rotates about its longitudinal axis (y-axis) at 100rpm. The moment of inertia matrix of the rod about the weld point O is

[𝑰O]=0.08kgm2[110110002].

Find the kinetic energy of the rod.

Solution The rod rotates about the fixed point O with angular velocity

𝝎=ωȷˆwhereω=100rpm=10.47rad/s.

The kinetic energy of a rigid body rotating about a fixed point O is given by

EK=12𝝎[[𝑰O]𝝎]

For the given problem, let us write the moment of inertia matrix [𝑰O] of the rod as margin: This step is simply to facilitate computation. We carry out the multiplication and at the end substitute the value of K0.

[𝑰O] = K0[110110002]
where K0=0.08kgm2. Now
𝝎 = ωȷˆ={0ω0}T.
Therefore,
[𝑰O]𝝎 = K0[110110002]{0ω0}
= K0{ωω0}
= K0ωıˆ+K0ωȷˆ.

Substituting the expression in the formula for EK we get

EK = 12ωȷˆ(K0ωıˆ+K0ωȷˆ)
= 12K0ω2
= 120.08kgm2(10.47rad/s)2
= 4.38Nm=4.38J

Answer: EK=4.38J

Problems for 20.4 Mechanics using [𝑰cm] and [𝑰O]

20.4.1

Find the angle between 𝝎 and 𝑯=[𝑰cm]𝝎, if 𝝎=3rad/s𝒌ˆ and

[𝑰cm]=[200002000040]kgm2.

20.4.2  Find the projection of 𝑯=[𝑰cm]𝝎 in the direction of 𝝎, if 𝝎=(2ıˆ3ȷˆ)rad/s and [𝑰cm]=[5.40005.400010.8]lbmft2.

20.4.3  Compute 𝝎×([𝑰cm]𝝎) for the following two cases. In each case 𝝎=2.5rad/sȷˆ.

  1. (a)

    [𝑰cm]=[8000160008]kgm2, Answer: ω×([Icm]ω)=𝟎kgm2/s2.

  2. (b)

    [𝑰cm]=[8000164048]kgm2. Answer: ω×([Icm]ω)=25ıˆkgm2/s2.

20.4.4  Find the z component of 𝑯=[𝑰cm]𝝎 for 𝝎=2rad/s𝒌ˆ and

[𝑰cm]=[10.202.5602.5610.2000014.00]lbmft2.

20.4.5  For 𝝎=5rad/s𝝀ˆ, find 𝝀ˆ such that the angle between 𝝎 and 𝑯=[𝑰cm]𝝎 is zero if [𝑰cm]=[302060203]kgm2.

20.4.6  Calculation of H/O and H˙/O. You will consider a disc spinning about various axes through the center of the disc O. For the most part, you are concerned with calculating 𝑯˙/O. But, you should keep in mind that the purpose of the calculation might be to calculate reaction forces and moments. Every point on the disk moves in circles about the same axis at the same number of radians per second. These problems can be done at least three different ways: using the[I]’s from problem 20.3, direct integration of the definition of 𝑯˙/O, and using the equivalent dumbbell (with the definition of 𝑯˙/O). You should do them all three ways, first using the way you are most comfortable with.

  1. (a)

    Assume that the configuration is as in problem 20.3(a) and that the disk is spinning about the z axis at the constant rate of ω𝐤. What are 𝑯/O and 𝑯˙/O?

  2. (b)

    Assume that the configuration is as in problem 20.3(b) and that the disk is spinning about the z axis at the constant rate of ω𝐤. What are 𝑯/O and 𝑯˙/O? (In this problem there are two ways to use [I]: one is in the xyz coordinate system, the other is in a coordinate system aligned nicely with the disk. )

20.4.7  What is the angular momentum 𝑯 for:

  1. (a)

    A uniform sphere rotating with angular velocity 𝝎 (radius R, mass m)?

  2. (b)

    A rod shown in the figure rotating about its center of mass with these angular velocities: 𝝎=ωıˆ, 𝝎=ωȷˆ, and 𝝎=ω(ıˆ+ȷˆ)/2?

Filename:pfigure-blue-99-1
Figure 20.97

20.4.8  A hoop is welded to a long rigid massless shaft which lies on a diameter of the hoop. The shaft is parallel to the x-axis. At the instant of interest the hoop is in the xy-plane of a coordinate system which has its origin O at the center of the hoop. The hoop spins about the x-axis at 1rad/s. There is no gravity.

  1. (a)

    What is [𝑰cm] of the hoop?

  2. (b)

    What is 𝑯/O of the hoop?

Filename:pfigure-blue-47-2
Figure 20.98

20.4.9  Two thin rods, each of mass m and length d, are welded perpendicularly to an axle of mass M and length , which is supported by a ball-and-socket joint at A and a journal bearing at B. A force 𝑭=Fıˆ+G𝒌ˆ acts at the point P as shown in the figure. Determine the initial angular acceleration α and initial reactions at point B.

Filename:pfigure-blue-89-1
Figure 20.99

20.4.10  The thin square plate of mass M and side c is mounted vertically on a vertical shaft which rotates with angular speed ω0. In the fixed x,y,z coordinate system with origin at the center of mass (with which the plate is currently aligned) what are the vector components of 𝝎? What are the components of 𝑯cm?

Filename:pfigure-blue-35-2
Figure 20.100

20.4.11  If you take any rigid body (with weight ), skewer it with an axle, hold the axle with hinges that cause no moment about the axle, and hold it so that it cannot slide along the axle, it will swing like a pendulum. That is, it will obey the equations:

θ¨ =Csinθ,or (a)
ϕ¨ =+Csinϕ, (b)

where

θ is the amount of rotation about the axis measuring zero when the center of mass is directly below the z-axis; i. e., when the center of mass is in a plane containing the z-axis and the gravity vector,

ϕ is the amount of rotation about the axis measuring zero when the center of mass is directly above the z-axis,

C=Mtotalgdsinγ/Izz,

Mtotal is the total mass of the body,

g is the gravitational constant,

d is the perpendicular distance from the axis to the center of mass of the object,

γ is the angle of tip of the axis from the vertical,

Izz is the moment of inertia of the object about the skewer axis.

  1. (a)

    How many examples of such pendula can you think of and what would be the values of the constants? (e.g. a simple pendulum?, a broom upside down on your hand?, a door?, a box tipping on one edge?, a bicycler stopped at a red light and stuck in her new toe clips?, a person who got rigor mortis and whose ankles went totally limp at the same time?,…)

  2. (b)

    A door with frictionless hinges is misaligned by half an inch (the top hinge is half an inch away from the vertical line above the other hinge). Estimate the period with which it will swing back and forth. Use reasonable magnitudes of any quantities you need.

  3. (c)

    Can you verify the formula (a simple pendulum and a swinging stick are special cases)?

Filename:pfigure4-3Dpend
Figure 20.101

-3-D pendulum on a hinge.

20.4.12  A hanging rod going in circles. A uniform rod with length and mass m hangs from a well greased ball and socket joint at O. It swings in circles at constant rate 𝝎=ω𝒌 while the rod makes an angle ϕ with the vertical (so it sweeps a cone).

  1. (a)

    Calculate 𝑯/O and 𝑯˙/O.

  2. (b)

    Using 𝑯˙/O and angular momentum balance, find the rotation rate ω in terms of ϕ, g, , and m. How long does it take to make one revolution? Answer: ω2=3g2Lcosϕ, Trev=2πω.

  3. (c)

    If another stick rotates twice as fast while making the same angle ϕ, what must be its length compared to ?

  4. (d)

    Is there a simple way to describe what two uniform sticks of different length and different cone angle ϕ both have in common if their rate of rotation ω is the same? (If you accurately draw two such sticks, you should see the relation.)

  5. (e)

    Find the reaction force at the ball and socket joint. Is this force parallel to the stick?

Filename:pfigure-s94h6p3
Figure 20.102

20.4.13  From discrete to continuous. n equal beads of mass m/n each are glued to a rigid massless rod at equal distances l=/n from each other. For n=4, the system is shown in the figure, but you are to consider the general problem with n masses. The rod swings around the vertical axis maintaining angle ϕ.

  1. (a)

    Find the rate of rotation of the system as a function of g,,n, and ϕ for constant rate circular motion. [Hint: You may need these series sums: 1+2+3++n=n(n+1)2 and 12+22+32++n2=n(n+1)(2n+1)6.] Answer: ω2=3gnL(2n+1)cosϕ.

  2. (b)

    Set n=1 in the expression derived in (a) and check the value with the rate of rotation of a conical pendulum of mass m and length .

  3. (c)

    Now, put n in the expression obtained in part (a). Does this rate of rotation make sense to you? [Hint: How does it compare with the rate of rotation for the uniform stick of Problem 20.4?]

Filename:pfigure-s94h6p4
Figure 20.103

20.4.14  Spinning dumbbell. Two point masses are connected by a rigid rod of negligible mass and length 22. The rod is welded to a shaft that is along the y axis. The shaft is spinning at a constant rate ω driven by a motor at A (Not shown). At the time of interest the dumbbell is in the yz-plane. Neglect gravity.

  1. (a)

    What is the moment of inertia matrix about point O for the dumbbell at the instant shown? Answer:

    [𝑰O]=m2[400022022],

    in terms of the coordinate system shown.

  2. (b)

    What is the angular momentum of the dumbbell at the instant shown about point O. Answer: H/O=2m2ω(ȷˆ+kˆ).

  3. (c)

    What is the rate of change of angular momentum of the dumbbell at the instant shown? Calculate the rate of change of angular momentum two different ways:

    • (i)

      by adding up the contribution of both of the masses to the rate of change of angular momentum (i.e. 𝑯˙=𝒓i×(mi𝒂i)) and

    • (ii)

      by calculating Ixy and Iyz, calculating 𝑯, and then calculating 𝑯˙ from 𝝎×𝑯. Answer: H˙/O=2m2ω2ıˆ.

  4. (d)

    What is the total force and moment (at the CM of the dumbbell) required to keep this motion going? Do the calculation in the following ways:

    • (i)

      Using angular momentum balance (about any point of your choice) for the shaft dumbbell system.

    • (ii)

      By drawing free-body diagrams of the masses alone and calculating the forces on the masses. Then use action and reaction and draw a free-body diagram of the remaining (massless) part of the shaft-dumbbell system. Then use moment balance for this system. Answer: F=𝟎, 𝑴O=2mL2ω2ıˆ.

Filename:pg84-3
Figure 20.104

Crooked spinning dumbbell.

20.4.15  A uniform rectangular plate of mass m, height h, and width b spins about its diagonal at a constant angular speed ω, as shown in the figure. The plate is supported by frictionless bearings at O and C. The plate lies in the YZ-plane at the instant shown.

[Note: cos(β)=b/(h2+b2), and sin(β)=h/(h2+b2).]

  1. (a)

    Compute the angular momentum of the plate about G (with respect to the body axis Gxyz shown).

  2. (b)

    Ignoring gravity determine the net torque (moment) that the supports impose on the body at the instant shown. (Express your solution using base vectors and components associated with the inertial XYZ frame.)

  3. (c)

    Calculate the reactions at O or C for the configuration shown. Answer: RO=mω2bh(b2h2)12(b2+h2)32kˆ=RC.

  4. (d)

    In the case when h=b=1m what is the vertical reaction at C (still neglecting gravity).

Filename:pfigure-blue-36-1
Figure 20.105

20.4.16  A uniform rectangular plate spins about a fixed axis. The constant torque My is applied about the Y axis. The plate is held by bearings at the corners of the plate at A and B. At the instant shown the plate and the yz axis attached to the plate are in the YZ plane. Neglect gravity. At t=0 the plate is at rest.

  1. (a)

    At t=0+ (i.e., just after the start), what is the acceleration of point D? Give your answer in terms of any or all of c,d,m,My,ıˆ,ȷˆ, and 𝒌ˆ (the base vectors aligned with the fixed XYZ axes), and IYYcm [You may use IYYcm in your answer. Or you could use this result from the change of coordinates formulae in linear algebra: IYYcm=(cos2β)Iyycm+(sin2β)Izzcm to get an answer in terms of c, d, m, and My.] Answer: aD=MyIYYcm2dcc2+d2ıˆ.

  2. (b)

    At t=0+, what is the reaction at B? Answer: B=MyIYZcm2IYYcmc2+d2ıˆ.

Filename:pfigure-s95q8
Figure 20.106

20.4.17  A uniform semi-circular disk of radius r=0.2m and mass m=1.5kg is attached to a shaft along its straight edge. The shaft is tipped very slightly from a vertical upright position. Find the dynamic reactions on the bearings supporting the shaft when the disk is in the vertical plane (xz) and below the shaft.

Filename:pfigure4-rpi
Figure 20.107

20.4.18  A thin triangular plate of base b, height h, and mass m is attached to a massless shaft which passes through the bearings L and R. The plate is initially vertical and then, provoked by a teensy disturbance, falls, pivoting around the shaft. The bearings are frictionless. Given,

𝒓cm = (0,23b,13h),
IxxO = m2(b2+h2/3),
IyyO = mh2/6,
IzzO = mb2/2,
Iyycm = mh2/18,
IzyO = mbh/4,
  1. (a)

    what is the plate’s angular velocity as it passes through the downward vertical position?

  2. (b)

    Compute the bearing reactions when the plate swings through the downward vertical position.

Filename:pfigure-blue-86-1
Figure 20.108

20.4.19  A welded structural tubing framework is proposed to support a heavy emergency searchlight. To select properly the tubing dimensions it is necessary to know the forces and moments at the critical section AA, when the framework is moving. The frame moves as shown.

  1. (a)

    Determine the magnitude and direction of the dynamic moments (bending and twisting) and forces (axial and shear) at section AA of the frame if, at the instant shown, the frame has angular velocity ωo and angular acceleration αo. Assume that the mass per unit length of the pipe is ρ.

Filename:pfigure-blue-88-1
Figure 20.109

20.4.20  A uniform cube with side 0.25m is glued to a pipe along the diagonal of one of its faces ABCD. The pipe rotates about the y-axis at a constant rate ω=30rpm. At the instant of interest, the face ABCD is in the xy-plane.

  1. (a)

    Find the velocity and acceleration of point F on the cube.

  2. (b)

    Find the position of the center of mass of the cube and compute its acceleration 𝒂cm.

  3. (c)

    Suppose that the entire mass of the cube is concentrated at its center of mass. What is the angular momentum of the cube about point A?

Filename:pfigure4-2-rp3
Figure 20.110

20.4.21  A uniform square plate of mass m and side length b is welded to a narrow shaft at one corner. The shaft is perpendicular to the plate and has length . The upper end of the shaft is connected to a ball and socket joint at O. The assembly is rotating about the z axis through O at the fixed rate 𝝎=ω𝒌ˆ which happens to be just right to make the plate horizontal and the shaft vertical, as shown. The gravitational constant is g.

  1. (a)

    What is ω in terms of m, , b, and g?

  2. (b)

    Is the reaction at O on the rod-plate assembly in the direction of the line from G to O? (Why?)

Filename:pfigure-s94q6p1
Figure 20.111

20.4.22  The uniform cube shown is rotating at constant angular velocity about the fixed axis OA with angular speed 10rad/s counterclockwise looking towards O from A. Find:

  1. (a)

    𝝎

  2. (b)

    𝒗B

  3. (c)

    𝒂C

  4. (d)

    𝑯/O

  5. (e)

    𝑯˙/O

  6. (f)

    The total torque applied to the block in the configuration shown during this motion.

Filename:pfigure-spinningbrick
Figure 20.112

20.4.23  A spinning box. A solid box with dimensions a, a, and 2a (where a=.1m) and mass = 1kg. It is supported with ball and socket joints at opposite diagonal corners at A and B. It is spinning at constant rate of (50/π)rev/sec and at the instant of interest is aligned with the fixed axis ıˆ,ȷˆ, and 𝒌ˆ. Neglect gravity.

  1. (a)

    What is 𝝎?

  2. (b)

    What is 𝝎˙?

  3. (c)

    What is 𝑯/cm?

  4. (d)

    What is 𝑯˙/cm?

  5. (e)

    What is 𝒓B/A×𝑭B (where 𝑭B is the reaction force on at B)? A dimensional vector is desired.

  6. (f)

    If the hinges were suddenly cut off when was in the configuration shown, what would be 𝝎˙ immediately afterwards?

Filename:pfigure-s94f1p2
Figure 20.113

20.4.24  Uniform rod attached to a shaft. A uniform 1kg, 2m long rod is attached to a shaft which spins at a constant rate of 1rad/s. One end of the rod is connected to the shaft by a ball-and-socket joint. The other end is held by two strings (or pin jointed rods) that are connected to a rigid cross bar. The cross bar is welded to the main shaft. Ignore gravity. Find the reaction force at A.

Filename:s97f2
Figure 20.114

20.4.25  A round disk spins crookedly on a shaft. A thin homogeneous disk of mass m and radius r is welded to the horizontal axle AB. The normal to the plane of the disk forms an angle β with the axle. The shaft rotates at the constant rate ω.

  1. (a)

    What is the x component of the angular momentum of the disk about point G the center of mass of the disk?

  2. (b)

    What is the x component of the rate of change of angular momentum of the disk about G?

  3. (c)

    At the given instant, what is the total force and moment (about G) required to keep the disk spinning at constant rate?

Filename:pg92-2
Figure 20.115

A round disk spins crookedly on a shaft.

20.4.26  A uniform disk of mass m=2kg and radius r=0.2m is mounted rigidly on a massless axle. The normal to the disk makes an angle β=60o with the axle. The axle rotates at a constant rate ω=60rpm. At the instant shown in the figure,

  1. (a)

    Find the angular momentum 𝑯/O of the disk about point O. Answer: H/O=(0.108ıˆ+0.126k)kgm2/s

  2. (b)

    Draw the angular momentum vector indicating its magnitude and direction. Answer: H/O is in the xz plane at the instant of interest and at an angle α41, clockwise from the z-axis.

  3. (c)

    As the disk rotates (with the axle), the angular momentum vector changes. Does it change in magnitude or direction, or both? Answer: H/O does not change in magnitude (m, R, ω, β are constants). H/O rotates with the disk; i.,e., it changes direction.

  4. (d)

    Find the direction of the net torque on the system at the instant shown. Note, only the direction is asked for. You can answer this question based on parts (b) and (c). Answer: H˙/O= the rate of change of H/O. The tip of H/O goes in a circle. At the instant of interest, the change will be in the ȷˆ or ȷˆ direction. Thus, MO=H˙/O will be in the ȷˆ direction.

Filename:summer95p2-3
Figure 20.116

20.4.27  A uniform disk spins at constant angular velocity. The uniform disk of mass m and radius r spins about the z axis at the constant rate ω. Its normal 𝒏ˆ is inclined from the z axis by the angle β and is in the xz plane at the instant of interest. What is the total force and moment (relative to O) needed to keep this disk spinning at this rate? Answer in terms of m, r, ω, β and any base vectors you need.

Filename:p-f96-p3-3
Figure 20.117

20.4.28  Spinning crooked disk. A rigid uniform disk of mass m and radius R is mounted to a shaft with a hinge. An apparatus not shown keeps the shaft rotating around the z axis at a constant ω. At the instant shown the hinge axis is in the x direction. A massless wire, perpendicular to the shaft, from the shaft at A to the edge of the disk at B keeps the normal to the plate (the z direction) at an angle ϕ with the shaft.

  1. (a)

    What is the tension in the wire? Answer in terms of some or all of R,m,ϕ and ω.

  2. (b)

    What is the reaction force at the hinge at O?

Filename:s97p3-3
Figure 20.118

20.4.29  A disoriented disk rotating with a shaft. A uniform thin circular disk with mass m=2kg and radius R=0.25m is mounted on a massless shaft as shown in the figure. The shaft spins at a constant speed ω=3rad/s. At the instant shown, the disk is slanted at 30 with the xy-plane. Ignore gravity.

  1. (a)

    Find the rate of change of angular momentum of the disk about its center of mass. Answer: H˙cm=93128Nmȷˆ=0.122Nmȷˆ

  2. (b)

    Find the rate of change of angular momentum about point A. Answer: H˙A=26.89Nmȷˆ

  3. (c)

    Find the total force and moment required at the supports to keep the motion going. Answer: RA=4.56Nıˆ,RC=13.439Nıˆ,Mz=0.

Filename:pfigure-s94h8p1
Figure 20.119

20.5 Dynamic balance

Sometimes when something spins it can shake the structure that holds it. A familiar example is an ‘unbalanced’ car tire which makes a car shake. But rotor balance is important in all kinds of machines. Most often one seeks to eliminate or minimize the imbalance. For example the strange design of car crank-shafts is due in large part to an attempt to minimize its imbalance. But what does it mean, in the language of mechanics, to say a rotating body is ‘unbalanced’? It means that non-zero forces and/or torques are required to hold the axis still while the object spins at constant rate. Going back to the basic momentum balance equations:

Linear momentum balance:𝑭=𝑳˙

and

Angular momentum balance:𝑴C=𝑯˙/C,

we see that forces and torques are required if 𝑳˙𝟎 or if 𝑯˙/C𝟎.

So imbalance means 𝑳˙ and/or 𝑯˙/C is not zero. We break the concept of balance into the following two concepts:

  1. (1)

    An object is said to be statically balanced with respect to a given axis of rotation if 𝑳˙=𝟎 when it spins at constant rate about that axis.

  2. (2)

    An object is said to be dynamically balanced with respect to a given axis of rotation if both 𝑳˙=𝟎 and 𝑯˙=𝟎 for constant rate rotation about that axis.

The origin of the words ‘static’ balance and ‘dynamic’ balance is in how the imbalance can be measured. Static imbalance can be measured with a static test, dynamic imbalance requires a dynamic test.

Static balance

If the center-of-mass of a rigid body is on the fixed axis of rotation then it will not accelerate. Thus, 𝑳˙=𝒂cmmtot=𝟎 and the object is statically balanced. An equivalent definition of static balance is that the net force on the spinning body is zero. Whether or not this net force is so can be tested with a statics experiment. Put the axis of rotation on good bearings and see if the mass hangs down in any preferred direction. In a tire shop, this kind of balancing is sometimes called bubble balancing because of the bubble in the level measuring device.

Dynamic balance

The first condition for dynamic balance of an object with respect to spinning about an axis is that the object be statically balanced. The center-of-mass must lie on the axis of rotation. The second condition for dynamic balance, that 𝑯˙=𝟎, is a little more subtle. To make things more specific, let’s calculate the rate of change of angular momentum 𝑯˙/O with respect to a point O that is on the axis of rotation. So, for constant rate rotation about a fixed axis we have:

𝑯˙/O=𝝎×𝑯/O (20.78)

because 𝑯/O spins with the body. 𝑯˙/O is evidently zero if the angular momentum, 𝑯/O , is parallel to the angular velocity, 𝝎. However, 𝑯/O=[𝑰O]𝝎 which means 𝑯/O is parallel to 𝝎 only when 𝝎 is an eigenvector of [𝑰O]margin: When you go to ACME garage and get your car tires balanced you can politely ask the mechanic: “Ma’am, could you make sure that one of the eigenvectors of my wheel’s moment of inertia matrix is parallel to the car axle? Thanks.” The mechanic, if she is a decent person, will make sure that it is the appropriate eigenvector that is parallel to the axle. Otherwise the wheel would point straight sideways with the axle piercing the rubber. .

So an object which is spinning about a fixed axis is dynamically balanced if its center-of-mass is on the axis (static balance) and the angular velocity vector 𝛚 is an eigenvector of the moment of inertia matrix.

If we restrict margin: Caution: A common misperception about angular momentum is that it is always parallel to angular velocity. In general, angular momentum is not parallel to angular velocity. When is angular momentum about the center-of-mass, say, parallel to angular velocity? When 𝝎 is an eigenvector of [𝑰cm]. This is always the case for planar objects rotating about an axis perpendicular to the plane, but not generally in 3D. our attention to cases where the axis of rotation is the z-axis then this condition is easy to recognize. It is when Ixz=Iyz=0, that is, when the only non-zero element in the third column and in the third row of the [𝑰O] matrix is IzzO in the lower right corner.

Since these terms cause wobbling if one of the coordinate axis is an axis of rotation, the terms that are off the main diagonal in the moment of inertia matrix are sometimes called the ‘imbalance’ terms.

These terms lead to dynamic imbalance when the object is spun about the x, y, or z-axes. The off diagonal terms are also sometimes called the ‘centrifugal’ terms. This naming has an intuitive basis. One way to understand dynamic imbalance is to think of it being due to the unbalanced centrifugal pull of bits of mass that are spinning in circles.

Often, but not always, you can tell by inspection if something is dynamically balanced for rotation about a certain axis. If, for every bit of mass that is not on the axis there is another equal bit of mass that is exactly opposite, with respect to the axis, then the object is balanced. Some objects that do not meet this symmetry condition are also dynamically balanced, however.

Example: Some common objects

All of the figures in fig. 2.34 on page 2.34 are dynamically balanced about all the axes shown and about any axis perpendicular to any pair of these axes.

Example: A cube = a sphere

Both a cube and a sphere have moment of inertia matrices proportional to the identity matrix

[𝑰cm]cube = m26[100010001]
[𝑰cm]sphere = mD210[100010001].

So, for any 𝝎, we get that 𝑯cm is parallel to 𝝎 and, thus, both of these objects are dynamically balanced for rotation about any axis through their centers of mass.

Here is an example where the mathematics contradicts intuition; it does not seem like a cube should be balanced for rotation about a random skewed axis through its center-of-mass!

SAMPLE 20.19  Static and dynamic balance in 2-D motion.

Filename:sfig4-7-DH1
Figure 20.120:

A system with two equal masses at A and B rotates about point O at a constant rate ω=10rpm. The system is shown in the figure. The rod connecting the two masses has negligible mass.

  1. 1.

    Is the system statically balanced? If not, suggest a way to balance it.

  2. 2.

    Is the system dynamically balanced? If not, suggest a way to balance it.

Solution

  1. 1.
    Filename:sfig4-7-DH2
    Figure 20.121:

    Since the two masses are equal, the center-of-mass of the system is at the geometric center of rod AB, i.e., at a point O, halfway between A and B. Clearly, O is not on the axis of rotation which passes through O and is perpendicular to the plane of motion. Therefore, the system is not statically balanced.

    To balance the system, we must move the center-of-mass to O or pivot the system at O. Say we cannot move the pivot point. So, to move the center-of-mass O, we add mass m to m at A. From the definition of center-of-mass:

    (m+m)r1 = mr2
     m = m(r2r1)r1
    = 1kg(0.5m0.3m)0.3m
    = 0.67kg.

    Answer: Add 0.67 kg to mass at A.

  2. 2.
    Filename:sfig4-7-DH3
    Figure 20.122:

    The system, as given, is not dynamically balanced since it is not statically balanced. Let us check if it is dynamically balanced after adding m to A as suggested above.

    𝑴O=𝑯˙/O.

    Let us calculate 𝑯˙/O to see if it is zero:

    𝑯˙/O = 𝒓1×m1𝒂1+𝒓2×m2𝒂2
    = (r1𝒆ˆR)×m1(ω2r1𝒆ˆR)+(r2𝒆ˆR)×(ω2r2𝒆ˆR)m2
    = m1ω2r12(𝒆ˆR×𝒆ˆR𝟎)m2ω2r22(𝒆ˆR×𝒆ˆR𝟎)
    = 𝟎.

    The net torque on the system is zero, therefore it is dynamically balanced. Note that 𝑯˙/O=𝟎 irrespective of the values of m1,m2,r1, and r2. Thus 𝑯˙/O=𝟎 even for the system as given. But the system, as given, is not dynamically balanced because static balance is a necessary condition for dynamic balance.

SAMPLE 20.20

Filename:sfig4-7-1
Figure 20.123: Imbalance of a rotor due to point masses spinning at a constant rate.

Imbalance of a rotor due to rotating masses. Two masses m1 and m2 are attached to a massless shaft AB by massless rigid rods of length h each. The two masses are separated by distance b along the shaft and are in the same plane but on the opposite sides of the shaft. The shaft is rotating with a constant angular speed ω. The shaft is free to move along the z-axis at point B. Ignore gravity.

  1. 1.

    Is the system statically balanced?

  2. 2.

    What is the torque required to keep the motion going about the z-axis (i.e., Mz) ?

  3. 3.

    What are the reactions at the support points of the shaft?

Solution

  1. 1.

    A simple line sketch and the free-body diagram of the system are shown in Fig. 20.124. The linear momentum balance ( 𝑭=m𝒂 ) for the system gives:

    Filename:sfig4-7-1a
    Figure 20.124: (a) A simple line diagram of the system. (b) Free-Body Diagram of the system.
    Filename:sfig4-7-1b
    Figure 20.125: Accelerations of m1 and m2.
    (Ax+Bx)ıˆ+(Ay+By)ȷˆ+Az𝒌ˆ = m1𝒂1+m2𝒂2
    = m1ω2h(ıˆ)+m2ω2hıˆ
     Ax+Bx = ω2h(m2m1) (20.81)
    Ay+By = 0 (20.82)
    Az = 0

    Clearly, 𝑭𝟎 which means the system is not statically balanced.

  2. 2.

    Mz can be easily calculated by writing angular momentum balance about axis AB. In fact, we do not even need to write the equation in this case; since all the forces pass through this axis and the accelerations of m1 and m2 also pass through it,

    M/AB(due to reaction forces)=0 and H˙/AB=0.

    But

    M/AB = H˙/AB
    or Mz + M/AB(due to reaction forces)=H˙/AB
     Mz = 0.

    Answer: Mz=0

  3. 3.

    For the four unknown reactions at A and B (we have already found Az and Mz) we have two scalar equations so far (from Linear Momentum Balance ). Angular Momentum Balance about point A gives:

    𝑴/A=𝑯˙/A

    Now,

    𝑴/A = 𝒓B/A×𝑭B
    = (2+b)𝒌ˆ×(Bxıˆ+Byȷˆ)
    = Bx(2+b)ȷˆBy(2+b)ıˆ
    𝑯˙/A = 𝒓C/A×m1𝒂1+𝒓D/A×m2𝒂2
    = (𝒌ˆ+hıˆ)×m1(ω2hıˆ)+((+b)𝒌ˆhıˆ)×m2(ω2hıˆ)
    = ω2h(m2(+b)m1)ȷˆ

    Equating 𝑴/A and 𝑯˙/A we get

    Bx(2+b)ȷˆBy(2+b)ıˆ=ω2h(m2(+b)m1)ȷˆ

    Dotting both sides of the equation with ıˆ and ȷˆ, we get

    By = 0
    and Bx = ω2h(2+b)[m2(+b)m1].

    Substituting Bx and By in (20.81) and (20.82) we find

    Ax = ω2h(2+b)[m2m1(+b)]
    Ay = 0

    Answer: Ax=ω2h(2+b)[m2m1(+b)]Ay=Az=0Bx=ω2h(2+b)[m2(+b)m1]By=0Mz=0

SAMPLE 20.21

Filename:sfig4-7-2
Figure 20.126: Imbalance of a rotor due to point masses spinning at a constant rate.

Balancing a rotor with spinning masses. Consider the same system as in the previous sample problem (Sample 20.120). Assume that m1=m2=m.

  1. 1.

    Set m1=m2 in the reactions calculated. Is the system statically balanced now? Explain.

  2. 2.

    Is the system dynamically balanced? Explain.

  3. 3.

    Balance the system dynamically by (a) adjusting the geometry, (b) adding two masses to the system.

Solution

  1. Filename:sfig4-7-2a
    Figure 20.127: Reaction forces at A and B are required to keep the motion going.
  2. 1.

    Recall from the solution of Sample 20.120 that

    Ax = ω2h(2+b)[m2m1(+b)]
    Bx = ω2h(2+b)[m2(+b)m1]
    Ay = By=Az=Mz=0

    Setting m1=m2=m in the above expressions for the reactions, we get

    Ax=ω2h(2+b)(mb) and Bx=ω2h(2+b)(mb).

    Thus, Ax+Bx=0. Therefore, 𝑭=𝟎 which means the system is statically balanced.

    Static Balance:

    When the two masses are equal, the center-of-mass of the system is on the axis of rotation. Therefore, the static reactions are zero irrespective of the orientations of the two masses.

    The requirement for static balance is 𝑭=𝟎 in any static orientation of the system. As long as the center-of-mass of the system lies on the axis of rotation, the system will be in static balance.

    Dynamic Balance:

    When the masses are rotating at constant speed, reaction forces at A and B are required to keep the motion going. Clearly, 𝑴𝟎 now. Therefore, the system is not in dynamic balance.

    At the instant shown the reaction forces are Ax and Bx as shown in Fig. 20.127, but they change directions as the position of the masses changes. For example, if the masses are in the y–z plane the reaction forces will be Ay and By. The magnitudes of these forces will remain the same. Thus, we have rotating reaction forces at A and B. These rotary forces cause wear in the bearings and induce vibrations in the frame of the machine and the supporting structure. Therefore, these forces are undesirable.

    If we can somehow make these dynamic reactions zero, the bearings, the machine frame, the supporting structure, the machine operator and the company will all be very happy. So, how do we do it? Read on.

  3. 2.

    There are many ways in which the rotor can be dynamically balanced:

    • Trivial solution: Remove both the masses. Of course, there is nothing left to produce any non-zero 𝑯˙. Thus,

      𝑴/any point=𝟎. Also, 𝑭=𝟎

      .

    • Adjust geometry: Set b=0. The center-of-mass is still on the axis of rotation.

       𝑭 = 𝟎
          In addition,
      𝑯˙/A = 𝟎
       𝑴/A = 𝟎.
    • Add two masses: If b is required to be non-zero, we can balance the rotor by adding two masses in any two selected transverse (to the shaft-axis) planes to make 𝑯˙/A=𝟎. margin: Balancing a rotor by adding two masses in two selected transverse (to the shaft axis) planes is quite a general method. Even if the rotor has more than two spinning masses the net angular momentum due to all masses can be reduced to the angular momentum produced by an equivalent system with just two masses such as the one under consideration. For example:

      • We can take two equal masses m and m, and place them in the opposite directions of the masses already on the shaft as shown in Fig. 20.128(a). It should be clear that the net angular momentum about any point on the shaft will be zero now.

        Filename:sfig4-7-2b
        Figure 20.128: Dynamic balancing of the rotor by adding two masses. (a) The two added masses (shown by dotted circles) are the same as the original masses on the shaft. (b) The two added masses are different from the original masses (mm).
      • We can add two masses m and m different from the masses attached to the shaft and place them a distance c apart on the shaft. Let the length of the connecting rods of the new masses be h. For the net angular momentum to be zero we need

        mω2hc𝑯˙due to ms = mω2hb𝑯˙due to ms
         mhc = mhb

        Thus, we have the freedom margin: In practice, this freedom is usually restricted by geometric and space constraints. of selecting any combination of m, h and c to give the required product. Here are two examples:

        1. (a)

          Let m=m/2,c=b/2, then h=mhbmc=4h.

        2. (b)

          Let m=m/2,c=2b, then h=mhbmc=h.

        Filename:sfig4-7-2c
        Figure 20.129: Dynamic balancing of the rotor by adding two masses. (a) m=m/2, c=b/2 and h=4h. (b) m=m/2, c=2b and h=h.

Problems for 20.5 Dynamic balance

20.5.1  Is 𝝎=(2ıˆ+5𝒌ˆ)rad/s an eigenvector of [𝑰cm]=[8000164048]kgm2?

20.5.2  Find the eigenvalues and eigenvectors of [𝑰cm]=[5000530310]kgm2.

20.5.3  Three equal length (1m) massless bars are welded to the rigid shaft AB which spins at constant rate 𝝎=ωıˆ. Three equal masses (1kg) are attached to the ends of the bars. At the instant of interest the masses are all in the xy plane as shown. Where should what masses be added to make the system shown dynamically balanced? [The solution is not unique.]

Filename:pfigure-s94q7p1
Figure 20.130

20.5.4  Dynamic balance of a system of particles. Three masses are mounted on a massless shaft AC with the help of massless rigid rods. Mass m1=1kg and mass m2=1.5kg. The shaft rotates with a constant angular velocity 𝝎=5rad/s𝒌ˆ. At the instant shown, the three masses are in the xz-plane. Assume frictionless bearings and ignore gravity.

  1. (a)

    Find the net reaction force (the total of all the forces) on the shaft. Is the system statically balanced?

  2. (b)

    Find the rate of change of angular momentum 𝑯˙/A, using 𝑯˙/A=i𝒓i/A×mi𝒂i.

  3. (c)

    Find the angular momentum about point A, 𝑯/A. Use this intermediate result to find the rate of change of angular momentum 𝑯˙/A, using 𝑯˙/A=𝝎×𝑯/A and verify the result obtained in (b).

  4. (d)

    What is the value of the resultant unbalanced torque on the shaft? Propose a design to balance the system by adding an appropriate mass in an appropriate location.

Filename:pfigure-s94h7p2
Figure 20.131

20.5.5  For what orientation of the rod ϕ is the system shown dynamically balanced?

Filename:ch4-7-c
Figure 20.132

20.5.6  A uniform square plate with 1m sides and 1kg mass is mounted (as shown) with a shaft that is in the plane of the square but not parallel to the sides. What is the net force and moment required to maintain rotation at a constant rate?

Filename:pfigure-blue-87d-1
Figure 20.133

20.5.7  An equilateral uniform triangular plate with sides of length plate spins at constant rate ω about an axis through its center of mass. The axis is coplanar with the triangle. The angle ϕ is arbitrary. When ϕ=0, the base of the triangle is parallel to the axis of rotation. Find the net force and moment required to maintain rotation of the plate at a constant rate for any ϕ. For what angle ϕ is the plate dynamically balanced?

Filename:Mikef91p3
Figure 20.134

20.5.8  An abstract problem. This problem concerns a system rigidly attached to a massless rod spinning about a fixed axis at a fixed rate. That is, whether or not the system is a rigid body, it is moving as if it were a rigid body. Every point in the system is moving in circles about the same axis at the same number of radians per second. The position of a particle of the system relative to the axis is 𝑹. The position relative to some point C is 𝒓C. In all cases, the body rotates about the axis shown. The total mass of all the objects is Mtot. Point O is at the center of mass of the system and is at a distance d from the rotation axis and is on the x-axis at the moment of interest. Assume there is no gravity. At the point A, the rod is free to slide in the z-direction. But, at point B, the rod is constrained from moving in the z-direction. The FBD of the system is shown in the figure. Referring to the free-body diagram, there are six scalar unknown reactions in this 3 dimensional problem.

  1. (a)

    The general case: Imagine that you know 𝑯˙/O=H˙0xıˆ+H˙0yȷˆ+H˙0z𝒌ˆ and also 𝑳˙=L˙xıˆ+L˙yȷˆ+L˙z𝒌ˆ. How would you find the six unknown reactions Ax, Ay, MAz, Bx, By, Bz? [Reduce the problem to one of solving six equations in six unknowns.]

  2. (b)

    The General Result: Is MAz necessarily equal to zero? Prove your result or find a counter-example. Answer: MAz is always zero.

  3. (c)

    Static Balance:. If d=0 is it true that all reaction forces are equal to zero? Prove your result or find a counter example. Answer: No. See problems 20.4 and 20.2 for counter-examples.

Filename:pfigure-s94h7p3
Figure 20.135

20.5.9  A uniform block spins at constant angular velocity. The magnitude of the angular velocity is ω. The center of mass G is on the axis of rotation. The ıˆ, ȷˆ, and 𝒌ˆ directions are parallel to the sides with lengths b, c and d, respectively. Answer in terms of some or all of  ω, m, b, c, d, ıˆ, ȷˆ, and 𝒌ˆ.

  1. (a)

    Find 𝒗H.

  2. (b)

    Find 𝒂H.

  3. (c)

    What is the sum of all the moments about point G of all the reaction forces at A and B?

  4. (d)

    Are there any values of b, c, and d for which this block is dynamically balanced for rotation about an axis through a pair of corners? Why?

Filename:p-f96-p2-3
Figure 20.136

20.5.10  Dynamic Balance. Consider each of the objects shown in turn (so to speak), each spinning about the z axis which is horizontal. In each case the center of mass is half way between the bearings (not shown). Gravity is in the ȷˆ direction. State whether it is possible to spin the object at constant rate with constant bearing reactions of mgȷˆ/2 and no other applied torques (‘yes’) or might not be (‘no’). Give a brief explanation of your answer (enough to distinguish an informed answer from a guess).

  1. (a)

    Axis along centerline of a uniform rectangular plate.

  2. (b)

    Axis orthogonal to uniform rectangular plate and through its centroid.

  3. (c)

    Axis orthogonal to uniform rectangular plate and not through its centroid.

  4. (d)

    Axis orthogonal to uniform but totally irregular shaped plate and through its centroid.

  5. (e)

    Axis along the diagonal of a uniform rhombus.

  6. (f)

    Axis through the center of a uniform sphere.

  7. (g)

    Axis through the diagonal of a uniform cube.

  8. (h)

    Axis with two equal masses attached as shown.

  9. (i)

    Axis with three equal masses attached as shown, all in a common plane with the axis, with ϕ not equal to 30.

  10. (j)

    An arbitrary rigid body which has center of mass on the axis and two of its moment-of-inertia eigenvectors perpendicular to the axis. (no picture)

Filename:s97p3-2
Figure 20.137

20.5.11  Static and Dynamic Balance A series of bodies, each of uniform density and each with total mass m, rotate at a constant angular speed ω about a fixed horizontal axis. Ignore gravity. For each body state whether the body is (i) statically balanced and whether it is (ii) dynamically balanced. Give clear arguments using words or equations to support your claims. (iii, iv) For each body you must either (a) add one point mass m or (b) add two point masses each of mass m/2 (your choice) that maintain static and dynamic balance if they are balanced, or that make the bodies statically and dynamically balanced. Justify your placement with words and/or equations. The masses need not be added to the bodies, but could be attached off the bodies by structures with negligible mass. [Hint: none of the placements are unique. You may draw a side view if that helps clarify your placement.]

  1. (a)

    A rectangular plate (height h, length ) mounted with the axle perpendicular to the plate and through its center.

  2. (b)

    The same plate as in (a) above but mounted at an angle ϕπ/2 from the shaft.

  3. (c)

    The numerals ‘203’ cut out of a plate and connected by massless rods. Each letter has mass m/3 and the three center-of-mass points of the individual letters are colinear and equally spaced. The shaft goes through the center of the ‘0’ and is perpendicular to the plane of the letters.

  4. (d)

    A sphere with radius R where the shaft passes a distance d<R from the center.

Filename:pfigure-dynbalancers
Figure 20.138

20.5.12  Dynamic tire balancing. Where should you place what weights in order to dynamically balance the tire?

Assume that a car tire can be modeled as a uniform disk of radius one foot and mass 20lbm. Assume the tire has been mounted crooked by one degree (the plane of the disk makes an 89 angle with the car axle). This crooked wheel is statically balanced since its center of mass is on the axis of rotation. But it is dynamically unbalanced since it wobbles and will wobble the car. Car shops put little weights on unbalanced wheels to balance them. They figure out where to put the weights by using a little machine that measures the wheel wobble. But you don’t need such a machine for this problem because you have been told where all the mass is located.

There are many approaches to solving this problem and there are also many correct solutions. You should find any correct solution. [one approach: guess a good location to put two masses and then figure out how big they have to be to make the tire dynamically balanced.]

Filename:pfigure-s94h7p5
Figure 20.139