Chapter 16 Circular motion of a rigid object

Download this chapter by section

All chapters and individual parts

Here we extend the idea of circular motion from particles to objects. A particle can go in circles. On a rigid object in 2D, all pairs of particles go in circles around each other. The key theoretical idea is of rigid-object rotation. The primary applications are pendulums, gear trains, and rotationally accelerating motors or brakes.

Filename:tfigure8-rel-ang-vel
Figure 16.1: All the points on a gear move in circles, assuming the axle is not translating.
Refer to caption
Filename:tfigure8-ang-vel-ex
Figure 16.2: All points on a flywheel move in circles. The points further from the center move faster, so they move further while the camera shutter is open, so they show a bigger motion blur (the constant-thickness spokes look fatter near the rim).

We covered the special case of circular motion of a particle in the previous chapter. Now we consider the kinematics and mechanics of rotating rigid objects in 2D. The key new conceptual idea is that of rotation of a rigid object.

Mechanics of circular motion

For the systems in this chapter, for every system we show in a free-body diagram we have, as always,

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

Because you already know how to work with forces and moments (the left sides of the top two equations), the primary new skill in this chapter is the evaluation of 𝑳˙,𝑯˙/C, and E˙K for a rotating rigid object. That is, you need to understand the position, velocity and acceleration of points moving in circles. The rest of the skills used are universal, for example solving algebraic or differential equations, plotting, etc.

16.1 Rotation of a rigid object in planar circular motion

Filename:tfigure8-ang-vel
Figure 16.3: a) An object, b) rotated counterclockwise an angle θ about 0.
Filename:tfigure8-ang-accel
Figure 16.4: Rotation of lines on a rotating rigid object. Some real or imagined lines marked on the rigid object are shown. They make the angles θ1, θ2, θ3, with respect to various fixed lines which do not rotate. As the object rotates, each of these angles increases by the same amount.

When two parts are glued together or attached by welding, gluing, several tight screws, bolts, rivets or the like we call the connection a ‘rigid attachment’. And, for the purposes of mechanics analysis, the two connected parts make up one bigger object. But most machines have various parts that are connected to each other, but not welded to each other. The most common such non-rigid attachment in engineering is a hinge. In 2D,

a hinge attachment between two objects keeps two points, one from each object, on top of each other while freely allowing relative rotation of the two objects about the hinge point.

In 3D a hinge keeps two lines, one on each body, coincident and allows relative rotation about that line. The common line, or in 2D the line orthogonal to the plane through the points, is called the hinge, the hinge axis, or the axis of relative rotation.

One example of a hinge is a car axle which allows rotation of a wheel relative to the car suspension. The hinge axis is the wheel axle

margin: The word axis is obviously related to the word axle. More generally the word axis means ’line’. For example the x and y axes are generally not axles about which anything rotates, they are just lines.

.

Physically, hinges are made various ways, sometimes by poking a cylindrical pin through the two objects and sometimes with ball bearings (see box 5.3 on page 5.3). So hinges are also called pin connections or bearings (fig. 16.3). In this chapter we limit our attention to a simple use of a hinge: one rigid part is hinged to a second fixed part that doesn’t move at all. Such a non-moving part can be thought of as connected to, or an extension of, the ground or ‘fixed frame’. In the simple case of interest in this chapter, one point on the moving part does not move and the rest of the part rotates about that point.

For definiteness and simplicity let’s call the hinge location 0 and the hinge axis through 0 the z axis.

One function of the hinge is to make the part’s only possible motion to be rotation about O. Thus to understand the dynamics of a hinged part we need to understand the position, velocity and acceleration of points on a rigid object which rotates. This whole section is about the kinematics (the geometry of motion) for this rotation. We will measure the amount of rotation by the angle θ, and the rate of change of θ by the angular velocity 𝝎 (‘omega’), and the rate of change of this angular velocity 𝝎 by the angular acceleration 𝜶 (‘alpha’).

As simple as this topic seems at first glance, you should pay close attention to the meanings and uses of these quantities. The rest of the book completely depends on the material in this section.

Box 16.1 Rotation is uniquely defined for a rigid object (2D)

Most people will find it self-evident that, starting with a rigid object at a reference orientation, all lines marked on the object rotate by the same angle θ. Here, for the doubting, we demonstrate this fact.

A rigid object is defined this way:

For every pair of material points A and B on a rigid object the distance |AB| between them does not change as the object moves.

In particular, when a rigid object rotates all distances between all pairs of points are preserved. Thus, by the “side-side-side” similar triangle theorem of elementary geometry, all relative angles between marked line segments are preserved by the rotation. For example, for a triangle ABC the angle at B is constant as the object rotates. Now consider any pair of line segments on the object.

[By ‘on’ the object we don’t mean a projected image drawn with a light pen that can move around on the object relative to the atoms. Rather, by ‘on the object’ we mean something defined by a particular set of atoms that make up the object.]

If the segments do not cross we can extend them to a point of intersection B. Such a pair of intersecting lines is shown here before and after rotation.Initially BA makes an angle θ0 with a horizontal reference line. BC then makes an angle of θABC+θ0. After rotation we measure the angle to the line BD. BA now makes an angle of θ0+θ. By the addition of angles in the rotated configuration line BC now makes an angle of θABC+θ0+θ which makes an increase by θ of the angle made by BC with the horizontal reference line. So both BA and BC rotate by the same angle θ.

We could use one of these two lines and compare it with an arbitrary third line through B and show that the third line also has equal rotation, and then a fourth, and so on. So all lines on the object through a point B rotate by the same angle θ. The demonstration for a pair of parallel lines, one of them through B, is easy, they stay parallel so always make a common angle with any reference line.

Because any line on the object either goes through B or is parallel to a line through B, all lines marked on a rigid object rotate by the same angle θ.

The rotation of a rigid object in 2D is thus unambiguously defined as the angle through which all lines on the object rotate.

Filename:tfigure-rotationisunique

Rotation of an object counterclockwise by θ

We start by imagining the object in a distinguished configuration which we call the reference configuration, reference state or reference position. For example we could take the left figure in fig. 16.3 as the reference configuration. If possible it’s usually best to pick the reference state to be one in which a prominent feature of the object is aligned with the x or y axes. The reference state may or may not be the start of the motion of interest. Even if not, we measure an object’s rotation by the change, relative to the reference state, in the counterclockwise angle θ of a reference line marked in the object relative to a fixed line outside. Which reference line? Fortunately,

All real or imagined lines marked on a rotating rigid object rotate by the same angle, the rotation angle, θ. (See box 16.1).

In three dimensions things are more complicated. General rotation of a rigid object is then represented not with a single angle θ, but rather with 3 angles, or with a unit vector and an angle, or a 3×3 matrix. So we wait to discuss which of the 2D ideas here generalize to 3D and which do not.

Rotated coordinates and base vectors ıˆ and ȷˆ

Filename:tfigure8-ang-accel-ex
Figure 16.5: a) A house is drawn by connecting lines between 6 points, b) the house and coordinate system are rotated, thus its coordinates in the rotating system do not change c) But the coordinates in the original system do change
Filename:tfigure8-ang-accel1
Figure 16.6: A rotating rigid object 𝒞 with rotating coordinates xy rigidly attached.

We pick two orthogonal lines on the rotating object and give them distinguished status as object-fixed (or body-fixed) rotating coordinate axes x and y. Think of these axes as xy coordinate axes on a piece of graph paper that is glued to the object (see fig. 16.5a&b). It’s easiest if we start by assuming that the xy axes have the same origin 0 as the xy axes and are parallel with the fixed xy axes when the object is in the reference configuration (when θ=0).

These rotating coordinate axes, x and y, have associated rotating base vectors ıˆ and ȷˆ (fig. 16.6 and 16.7). So ıˆ is always in the x direction and ȷˆ always in the y direction. We will use these rotating coordinates and base vectors to keep track of some particle of interest P that is ‘glued’ to the object. To start, note that particle P which is glued to the object has x and y coordinates that don’t change as the rotation progresses.

Example: A particle on the x axis

If a particle P is fixed on the x-axis at position x=3cm, then we have.

𝒓P = 3cmıˆ

for all time, even as the object rotates.

The position vector of a point P fixed to a rigid object hinged at O remains, as the rotation progresses,

𝒓P = xıˆ+yȷˆ, (16.1)

with x and y both constant. These rotating coordinate system components, [𝒓]xy=[x,y], are sometimes written as [𝒓]xy=[xy].

You will see that much of the math for rotating xy coordinates is reminiscent of that for polar coordinates. However, the spirit is a bit different. In polar coordinates the 𝒆ˆR axis was picked to track a particular particle of interest. Here we pick axes that rotate with an extended object and use that one set of axes to track any and all particles of interest.

Filename:sfig8-2-3
Figure 16.7: Fixed coordinate axes and rotating coordinate axes.

Note, even though neither x nor y change as θ changes, the point P they describe moves, in circles actually. How can the particle’s position change if its coordinates don’t change? Well, in eqn. (16.3) the change in position is represented by the base vectors changing as the object rotates. Thus we could write more explicitly that

𝒓P = xıˆ(θ)+yȷˆ(θ). (16.3)

Here we show more explicitly that the base vectors ıˆ and ȷˆ depend on θ.

Filename:sfig8-2-3a
Figure 16.8: The x and y coordinates of a point fixed on a rotating object stay constant while the base vectors ıˆ and ȷˆ rotate with the object.

Just like for polar base vectors (see eqn. (15.1) on page 15.1) we can express the rotating base vectors in terms of the fixed base vectors and θ.

ıˆ = cosθıˆ+sinθȷˆ, (16.4)
ȷˆ = sinθıˆ+cosθȷˆ.

Also we can express the fixed basis vectors in terms of the rotating vectors like this:

ıˆ = cosθıˆsinθȷˆ (16.5)
ȷˆ = sinθıˆ+cosθȷˆ.

Please review the section on dot products, 1.2, to see one derivation of these formulae.

We will use the phrase reference frame or just frame to mean “a coordinate system attached to a rigid object”. One can think of the coordinate grid as like an invisible metal framework (hence the word ‘frame’) that rotates with the object. We refer to a calculation based on the rotating coordinates in fig. 16.6 variously as “in the frame 𝒞” or “using the xy frame” or “in the ıˆȷˆ frame

margin: Advanced aside. Sometimes a reference frame is defined as the set of all coordinate systems that could be attached to a rigid object. Two coordinate systems, even if rotated with respect to each other, then represent the same frame so long as they rotate together. Some of the results we will develop only depend on this more slack definition of frame, that the coordinates are glued to the object with no mind of their orientation in the reference configuration.

.

In computer calculations we usually manipulate lists and arrays of numbers and not geometric vectors. So on a computer we keep track of vectors by keeping track of their lists of components. Let’s look at a point fixed to the object and whose coordinates we know in the reference configuration:

[𝒓Pref]xy=[xrefyref].

Assuming the object axes and fixed axes coincide in the reference configuration, the object coordinates of a point [𝒓P]xy are equal to the space fixed coordinates of the point in the reference configuration [𝒓Pref]xy. We can think of the point as defined either way, so

[𝒓P]xy=[𝒓Pref]xy.

The rotation matrix [R]

Here is a question we often need to answer: What are the fixed basis coordinates of a point that has the rotating-frame coordinates [𝒓]xy=[xy]?

Here is one way to find the answer:

𝒓P = xıˆ+yȷˆ (16.6)
= x(cosθıˆ+sinθȷˆ)+y(sinθıˆ+cosθȷˆ)
= ((cosθ)x(sinθ)y)xıˆ+((sinθ)x+(cosθ)y)yȷˆ

so we can pull out the x and y coordinates compactly as,

[𝒓P]xy=[xy]=[cosθx+sinθ(y)sinθx+cosθ(y)]. (16.7)

But this can, in turn be written in matrix notation as

[xy] = [cosθsinθsinθcosθ][R][xy],or (16.14)
[𝒓P]xy = [R][𝒓P]xy,or (16.15)
[𝒓P]xy = [R][𝒓Pref]xy.

The matrix [R] or [R(θ)] is the rotation matrix for counterclockwise rotations by θ. As shown above, if you know the coordinates of a point fixed on an object before rotation, you can find its coordinates after rotation by multiplying the coordinate column vector by the matrix [R]. You can remember what [R] is by remembering its components or by remembering that

the first and second column of [R] are the components of ıˆ and ȷˆ, respectively, in the fixed coordinate system.

For example, the first column of [R] consists of the x and y components of ıˆ.

A feature of eqn. (16.14) is that the same matrix [R] prescribes the coordinate change for every different point on the object. Thus for points called 1, 2 and 3 we have

[x1y1] =[R][x1y1],[x2y2]=[R][x2y2] &
[x3y3] =[R][x3y3].

A more compact way to write a matrix times a list of column vectors is to arrange the column vectors one next to the other in a matrix. By multiplying this matrix by [R] we get a new matrix whose columns are the new coordinates of various points. For example,

[x1x2x3y1y2y3]=[R][x1x2x3y1y2y3]. (16.16)

Eqn. 16.16 is useful for computer animation of rotating things in video games (and in dynamics simulations too) where points 1,2, and 3 are points on an object.

Example: Rotate a picture

If a simple picture of a house is drawn by connecting the six points (fig. 16.5a) with the first point at (x,y)=(1,2), the second at (x,y)=(3,2), etc., and the sixth point on top of the first, we have,

[xy points BEFORE][133211224542].

After a 30 counter-clockwise rotation about O, the coordinates of the house, in a coordinate system that rotates with the house, are unchanged (fig. 16.5b). But in the fixed (non-rotating, Newtonian) coordinate system the new coordinates of the rotated house points are,

[xy points AFTER] = [R][xy points BEFORE]=[R][xy points]
= [3/2.5.53/2][133211224542]
[0.11.60.60.81.10.12.23.25.05.34.02.2]

as shown in fig. 16.5c.

SAMPLE 16.1

Filename:sfig8-2-3b
Figure 16.9:

Computing rotated position of an object: A rigid object AOB is pinned at point O and is free to rotate about this point. Using the rotation matrix, find the coordinates of points A and B when θ=30 and 110 respectively.


margin:

Solution Let xy be a set of axes glued to the object with origin at O which is also the origin of the space-fixed coordinate axes xy. We first write the coordinates of points A and B in the body-fixed axes.

𝒓A|xy=[1m0],𝒓B|xy=[1m1m].

The rotation matrix for the xy axes rotating counterclockwise along with the object is

R=[cosθsinθsinθcosθ].

We can find the coordinates of A and B in the rotated position by multiplying their xy coordinates with the rotation matrix. We first write the coordinates of both points A and B in a single matrix, one column for each, and multiply with R to get the new coordinates:

[𝒓A𝒓B]xy=[cosθsinθsinθcosθ][xAxByAyB].

Now, we calculate the new coordinates for the given values of θ.

  • For θ=30, we have,

    [𝒓A𝒓B]xy = [0.8660.50.50.866][1m1m01m]
    = [0.866m0.366m0.5m1.366m].
    Filename:sfig8-3-1
    Figure 16.10: Coordinates of points A and B after a rotation of 30.

    Thus,

    Answer: [𝒓A]xy=[0.866m0.5m],[𝒓B]xy=[0.366m1.366m].

  • Similarly, for θ=110, we get,

    [𝒓A𝒓B]xy = [0.3420.940.940.342][1m1m01m]
    = [0.342m1.282m0.94m0.598m].
    Filename:sfig8-3-1a
    Figure 16.11: Coordinates of points A and B after a rotation of 110.

    Answer: [𝒓A]xy=[0.342m0.94m],[𝒓B]xy=[1.282m0.598m].

SAMPLE 16.2

Filename:sfig8-5-wiper
Figure 16.12:

Computer program for object rotation: The object shown in the figure rotates counterclockwise with its angular position θ given by θ(t)=θ˙0t where θ˙0=10/s. Write a computer program to animate the motion of the object by plotting its position at specified time instants. Use your program to plot the position of the object every two seconds for a total of 18 seconds.

Solution

In order to draw the object on the computer screen, we need to define it by taking a sufficient number of points on its boundary so that plotting all those points with connected line segments represents the object as closely as possible. In this case, we can select all the corner points on its boundary and define the object with the coordinates of these points. With the given geometry, it is fairly easy to find these coordinates. Denoting the coordinates of the kth point by (xk,yk), we can define this object by the following set of coordinates:

object=[x1x2x6x1y1y2y6y1]=[02442001122111].
Filename:sfig8-5-wiper-a
Figure 16.13: Coordinates of the six corner points are sufficient to define this object.

Note that the last column is a repeat of the first column. This is essential so that when we plot the object using the x and y coordinates listed here, we get a closed boundary.

Animation of the angular motion of this object basically requires determining the rotated position of the object and plotting it at sufficiently small time intervals. To do this, we need to define the initial orientation (position), define the time increment, find the angle of rotation θ corresponding to the new time, compute the corresponding rotation matrix, multiply the object coordinates with the rotation matrix, separate out the x and y coordinates, and plot x vs y. We then repeat the whole sequence until we reach the final time.

The following pseudocode can be easily adapted to a computer program to do the required animation.

object = [ 0  2  4  4  2  0  0
          -1 -1 -2  2  1  1 -1]     % coordinates of points
t = 0Ψ                              % specify initial time
t_final = 36                        % specify final time
delta_t = 2                         % specify time increment
thetadot0 = 10*pi/180               % specify thetadot0 in rad/s

while t < t_final
      t = t + delta_t               % increment the time
      theta = thetadot0*t           % compute current theta
      R = [ cos(theta) -sin(theta)  % find the rotation matrix
            sin(theta)  cos(theta)]
      new_position = R*object       % find new coordinates
      x = new_position(1st row)     % extract all x-coordinates
      y = new_position(2nd row)     % extract all y-coordinates
      plot y vs x                   % plot the new position
end

The plot obtained by implementing this code is shown in fig. 16.14.

Filename:sfig8-5-2disks
Figure 16.14: Rotated position of the object plotted at 2 second intervals.

Problems for 16.1 Rotation of a rigid body

Preparatory Problems

16.1.1   Give the definition of each term below in words and, if possible, with an equation.

  1. (a)

    rotation angle

  2. (b)

    (x,y) Answer: Rotating frame coordinates of a point.

  3. (c)

    (x,y) Answer: Fixed basis coordinates of a point.

  4. (d)

    rotation matrix Answer: A matrix that when operated on a column vector of fixed basis coordinates gives the coordinates in a rotated frame.
    [cosθsinθsinθcosθ]

16.1.2  Assume that in the reference configuration, when θ=0, that the xy axes attached to an object are aligned with the xy axes. Consider a point P attached to the moving frame (object) that has coordinates (x,y). Find each of the quantities below in terms of some or all of x,y,ıˆ,ȷˆ,θ,θ˙ and θ¨.

  1. (a)

    𝒓P Answer: [cosθsinθsinθcosθ][xy]

  2. (b)

    The rotation matrix [R]. Answer: [R]=[cosθsinθsinθcosθ]

16.1.3  Assume that in the reference configuration, when θ=0, that the xy axes are aligned with the xy axes. Consider two points P1 and P2 attached to the moving frame. P1 and P2 have coordinates (x1,y1) (x2,y2). Find each of the quantities below in terms of some or all of x1,y1,x2,y2,ıˆ,ȷˆ,θ,θ˙ and θ¨.

  1. (a)

    𝒓P1/P2 Answer: rP1/P2=[cosθsinθsinθcosθ][(x1x2)(y1y2)]

  2. (b)

    The rotation matrix [R]. Answer: [R]=[cosθsinθsinθcosθ]

16.1.4  The square plate OABC shown in the figure measures =1m on a side. The coordinate axes xy are fixed to the plate and the axes xy are fixed in space. The plate is rotated by an angle θ=60 as shown.

  1. (a)

    Find the coordinates x and y of the corner point B in the rotated frame.

  2. (b)

    Find the coordinates x and y of the same point in the fixed frame without using the rotation matrix.

  3. (c)

    Write the rotation matrix R for the given rotation.

  4. (d)

    Find the coordinates of point B in the fixed frame using its (x,y) coordinates and the rotation matrix R.

  5. (e)

    Find the coordinates of point B in the rotated frame using its (x,y) coordinates and the rotation matrix R in its appropriate form.

Answer: [rB]xy=[1m1m],[rB]xy=[0.37m1.37m], R=[0.50.870.870.5]

Filename:sfig8-4-1
Figure 16.15:

16.1.5  The rod shown in the figure rotates counterclockwise, starting from θ=0. Find the following quantities.

  1. (a)

    Write the position vector of points A and B as a list of numbers in an array, e.g., [𝐫A]xy=[xAyA], at θ=0. Answer: [rA]xy=[1m0m], [𝒓B]xy=[2m0m]

  2. (b)

    Write the rotation matrix R for θ=45. Answer: R=[12121212]

  3. (c)

    Find the coordinates of points A and B using the rotation matrix R for θ=45. Answer: [rArB]xy=[122122]m

  4. (d)

    If the rod rotates at a constant angular speed ω=π/3rad/s, find the coordinates of points A and B after 2 seconds, assuming the rod is at θ=0 when t=0.

Filename:sfig8-4-1a
Figure 16.16:

16.1.6  Find the angle of rotation corresponding to the following rotation matrices:

  1. (a)

    [0.70710.70710.70710.7071]. Answer: θ=45

  2. (b)

    [0.70710.70710.70710.7071]. Answer: θ=135

  3. (c)

    [0.70710.70710.70710.7071]. Answer: θ=315

More-Involved Problems

16.1.7  The rod shown in the figure rotates with constant angular speed ω=10rad/s. θ=0 at t=0s.

  1. (a)

    Find the position vector 𝒓B|t1 and the velocity 𝒗B|t1 of point B at t1=1s. Answer: rB|t1=1.68ıˆ1.09ȷˆ and 𝒗B|t1=20m/s[0.840.54]

  2. (b)

    Find the position vector 𝒓B|t2 of point B at t2=1.2s. What is the net displacement of point B during the time interval Δt=t2t1? Is this net displacement equal to 𝒗B|t1Δt? Why not? Answer: Net displacement =(3.366ıˆ+0.016ȷˆ)m. The net displacement is not equal to 𝒗B|t1Δt because the motion is not along a straight line.

Filename:sfig8-6-2
Figure 16.17:

16.1.8  Write a computer program to animate the rotation of an object. Your input should be a set of x and y coordinates defining the object (such that plot y vs x draws the object on the screen) and the rotation angle θ. The output should be the rotated coordinates of the object.

  1. (a)

    From the geometric information given in the figure, generate coordinates of enough points to define the given object. Answer: (0, 0); (153, 15); (303, 0); (0, 303); (15, 153).

  2. (b)

    Using your program, plot the object at θ=20, 60, 100, 160, and 270.

  3. (c)

    Assume that the object rotates with constant angular speed ω=2rad/s. Find and plot the position of the object at t=1s, 2s,and3s.

Filename:sfig8-6-2a
Figure 16.18:

16.1.9  The arrow shaped object shown in the figure is defined by the five points whose coordinates are given in some normalized units. The position shown is at t=0. The time dependent angular position of the object is given by θ(t)=C1t2+C2t where C1=0.25rad/s2 and C2=0.1rad/s. Animate the motion of the object and show its position from t=0 to t=5s at every second.

Filename:efig1-2-28
Figure 16.19:

16.2 Angular velocity and acceleration of a rigid object in planar circular motion

Angular velocity of a rigid object: 𝝎

Thus far we have talked about rotation, but not how it varies in time. Dynamics is about motion, velocities and accelerations, so we need to think about rotation rates and rotational accelerations.

A 2D rigid object’s net rotation is measured by the rotation angle θ. Thus, the simplest measure of rotation rate is θ˙dθdt. Because all marked lines rotate the same amount θ they all have the same rates of change. So θ˙1=θ˙2=θ˙3= etc.  and, as for rotation, the concept of rotation rate of a rigid object transcends the concept of rotation rate of this or that particular line. We give this rotation rate of a rigid object a special name, angular velocity, and symbol, ω (omega).

Repeating, for all lines marked on a rigid object,

ωθ˙1=θ˙2=θ˙3==θ˙. (16.23)

Often we think of angular velocity as a vector 𝝎. Its direction is the axis of the rotation which, for objects in the xy plane is 𝒌ˆ, pointed in the +z direction normal to the xy plane. The scalar part of 𝝎 is ω. So, the angular velocity vector is

𝝎ω𝒌ˆ (16.24)

with ω as defined in eqn. (16.23). Note 𝝎 is the angular velocity of the object (and of every line on it)

margin: We could also think of rotation rate as a derivative of the rotation matrix [R]. This point of view and its relation to the vector point of view here, is described in Box 16.2 on page 16.2.

.

The usual sign convention for ω is the same as that for θ, positive is counter-clockwise (CCW) around the origin. That’s the direction your fingers wrap if you point your thumb in the +z direction and grab the z axis.

Rate of change of ıˆ, ȷˆ

Our first use of the angular velocity vector 𝝎 is to calculate the rate of change of the rotating unit base vectors ıˆ and ȷˆ. We can find the rate of change of, say, ıˆ, by taking the time derivative of the first eqn. (16.4), and using the chain rule while recognizing that θ=θ(t). We can also make an analogy with polar coordinates (page 15.1), where we think of 𝒆ˆR as like ıˆ and 𝒆ˆθ as like ȷˆ. We found there that 𝒆ˆ˙R=θ˙𝒆ˆθ and 𝒆ˆ˙θ=θ˙𝒆ˆR. Either from differentiating ıˆ or from the correlation with polar coordinates, we get

ıˆ˙=θ˙ȷˆ or ıˆ˙=𝝎×ıˆandȷˆ˙=θ˙ıˆ or ȷˆ˙=𝝎×ȷˆ (16.25)

because ȷˆ=𝒌ˆ×ıˆ and ıˆ=𝒌ˆ×ȷˆ. Depending on the tastes of your lecturer, you may find Equations 16.25 are the most used equations from this point onward

margin: Eqn. 16.25 is sometimes considered the definition of 𝝎. In this view, popularized by Tom Kane, 𝝎 is that vector which determines ıˆ˙ and ȷˆ˙ by the formulas ıˆ˙=𝝎×ıˆ and ȷˆ˙=𝝎×ȷˆ. Then one needs to show that such a vector exists and that it is 𝝎=ıˆ×ıˆ˙. Luckily this reasoning leads to the same 𝝎 as our 𝝎=θ˙𝒌ˆ.

.

Box 16.2 The fixed Newtonian reference frame

Now we can reconsider the concept of a Newtonian frame, a concept which we had to assume to write the equations of dynamics in the first place.

All of mechanics depends, of course, on the laws of mechanics. The laws of mechanics are equations which involve, in part, the positions of things as a function of time. But how position is perceived to change in time depends on your reference frame. And some reference frames are better than others. The best, from our point of view, are reference frames in which Newton’s laws are accurate. Such a reference frame is called a Newtonian frame. In engineering practice the frames we use as approximations of a Newtonian frame often seem, loosely speaking, somehow still. So we sometimes call such a frame the fixed frame and label it with a script capital . When we talk about velocity and acceleration of mass points, for use in the equations of mechanics, we are always talking about the velocity and acceleration relative to a ixed, or equivalently, Newtonian frame.

Assume x and y are the coordinates of a vector 𝒓P and is a fixed frame with fixed axis (with associated constant base vectors ıˆ and ȷˆ). When we write 𝒓˙P we mean x˙ıˆ+y˙ȷˆ. But we could be more explicit (and notationally ornate) and write the velocity of P in the Newtonian frame as

d𝒓Pdt𝒓˙Pby which we meanx˙ıˆ+y˙ȷˆ.

The in front of the time derivative (or in front of the dot) means that when we calculate a derivative we hold the base vectors of constant. This is no surprise, because for the base vectors are constant. In general, however, when taking a derivative in a given frame you

  • write vectors in terms of base vectors stuck to the frame, and

  • only differentiate the components.

But only for the ixed or 𝒩ewtonian frame will accelerations calculated this way be directly applicable to Newton’s laws.

We will avoid the ornate notation of labeling frames when it is not needed. For example, if you don’t see any script capital letters floating around in front of derivatives, you can assume that we are taking derivatives relative to a fixed Newtonian frame.

Velocity of a point fixed on a rigid object

Let’s call some rotating object (script capital B) to which is glued a coordinate system xy with base vectors ıˆ and ȷˆ.

We now introduce the concept of derivative in a frame which we write, for the frame , as

ddt

which means, in words, the rate of change of something as viewed in the rotating frame . Now consider a point P at 𝒓P that is glued to the object. That is, the x and y coordinates of 𝒓P do not change in time.

Box 16.3 Plato’s discussion of spinning in circles as motion (or not)

Here is how a smart person, Plato, thought about these things two thousand years ago. Here he writes about an imagined discussion between Socrates and Glaucon about how an object can maintain contradictory attributes simultaneously.

“Socrates: Now let’s have a more precise agreement so that we won’t have any grounds for dispute as we proceed. If someone were to say of a human being standing still, but moving his hands and head, that the same man at the same time stands still and moves, I don’t suppose we’d claim that it should be said like that, but rather that one part of him stands still and another moves. Isn’t that so?

Glaucon: Yes it is.

Socrates: Then if the man who says this should become still more charming and make the subtle point that tops as wholes stand still and move at the same time when the peg is fixed in the same place and they spin, or that anything else going around in a circle on the same spot does this too, we wouldn’t accept it because it’s not with respect to the same part of themselves that such things are at the same time both at rest and in motion. But we’d say that they have in them both a straight and a circumference; and with respect to the straight they stand still since they don’t lean in any direction –while with respect to the circumference they move in a circle; and when the straight inclines to the right the left, forward, or backward at the same time that it’s spinning, then in no way does it stand still.

Glaucon: And we’d be right.”

This chapter is about things that are still with respect to their own parts (they do not distort) but in which the points do move in circles.

d𝒓Pdt𝒓˙P=x˙ıˆ+y˙ȷˆ=𝟎.

That is, relative to a moving frame, the velocity of a point glued to the frame is zero (no surprise).

We would like to know the velocity of such a point in the fixed frame. We just take the derivative, using the product rule and the differentiation rules we have developed for the rotating base vectors:

Filename:pfigure-blue-118-2
Figure 16.20: Velocity and acceleration of two points on a rigid object rotating about 0.
𝒓P = xıˆ+yȷˆ
 𝒗P=𝒓˙P = ddt(xıˆ+yȷˆ)=xıˆ˙+yȷˆ˙=x(𝝎×ıˆ)+y(𝝎×ȷˆ)
= 𝝎×(xıˆ+yȷˆ)

where 𝒓˙P is the simple way to write d𝒓Pdt.

Box 16.4 Acceleration of a point on a rotating body, using 𝝎

Leaving off the ornate pre-super-script for simplicity, we have

𝒂P=𝒗˙P = ddt(ddt(xıˆ+yȷˆ)) (16.26)
= ddt(x(𝝎×ıˆ)+y(𝝎×ȷˆ)).

To continue we need to use the product rule of differentiation for the cross product of two time dependent vectors like this:

ddt(𝝎×ıˆ) =𝝎˙×ıˆ+𝝎×ıˆ˙ (16.27)
=𝝎˙×ıˆ+𝝎×(𝝎×ıˆ), (16.28)
ddt(𝝎×ȷˆ) =𝝎˙×ȷˆ+𝝎×ȷˆ˙ (16.29)
=𝝎˙×ȷˆ+𝝎×(𝝎×ȷˆ). (16.30)

Substituting back into eqn. (16.26) we get

𝒂P =x(𝝎˙×ıˆ+𝝎×(𝝎×ıˆ))
+y(𝝎˙×ȷˆ+𝝎×(𝝎×ȷˆ))
=𝝎˙×(xıˆ+yȷˆ)+𝝎×(𝝎×(xıˆ+yȷˆ)))
=𝝎˙×𝒓P+𝝎×(𝝎×𝒓P)

which derives eqn. (16.32).

Thus,

𝒗P=𝝎×𝒓P (16.31)

We can rewrite eqn. (16.31) in a minimalist or elaborate notation as

𝒗 = 𝝎×𝒓or
d𝒓Pdt = 𝝎/×𝒓P/O.

Both are correct. In the first case you have to use common sense to know what point you are talking about and that it is on an object rotating with absolute angular velocity 𝝎. In the second case everything is laid out perfectly clearly (which is why it looks perfectly confusing). On the left side of the equation it says that we are interested in how point P moves relative to, not just any frame, but the fixed frame . On the right side we make clear that the rotation rate we are looking at is that of object relative to and not some other relative rotation. We further make clear that the formula only makes sense if the position of the point P is measured relative to a point which doesn’t move, namely 0.

What we have just found largely duplicates what we already learned in section 7.1 for points moving in circles. The slight generalization is that the same angular velocity 𝝎 can be used to calculate the velocities of multiple points on one rigid object. But the key idea remains: the velocity of a point going in circles is tangent to the circle it is going around and with magnitude proportional both to distance from the center and to the angular rate of rotation (fig. 16.20a).

Acceleration of a point on a rotating rigid object

Let’s again consider a point stuck on a rotating object and with position

𝒓P=xıˆ+yȷˆ.

Relative to the frame to which a point is attached, its acceleration is zero (again no surprise). But what is its acceleration in the fixed frame? We find this acceleration by writing the position vector and then differentiating twice, repeatedly using the product rule and eqn. (16.25) we get (see Box 16.2 for the details).

𝒂P=𝝎˙×𝒓P+𝝎×(𝝎×𝒓P) (16.32)

which is hardly intuitive at a glance

margin: Although the form eqn. (16.32) is not of much immediate use, if you are going to continue on to the mechanics of mechanisms or three-dimensional mechanics, you should follow the derivation of eqn. (16.32) carefully.

.

Recalling that in 2D 𝝎=ω𝒌ˆ we can use either the right hand rule or manipulation of unit vectors to rewrite eqn. (16.32) as

𝒂P=ω˙𝒌ˆ×𝒓Pω2𝒓P (16.33)

where ω=θ˙ and ω˙=θ¨

margin: The derivation of eqn. (16.33) can be written more informally and briefly like this (using minimalist notation): 𝒂 = 𝒗˙=ddt(𝝎×𝒓) = 𝝎˙×𝒓+𝝎×𝒓˙ = 𝝎˙×𝒓+𝝎×(𝝎×𝒓) = 𝝎˙×𝒓ω2𝒓.

.

Thus, as we found in section 15.1 for a particle going in circles, the acceleration can be written as the sum of two terms, a tangential acceleration ω˙𝒌ˆ×𝒓P due to increasing tangential speed, and a centrally directed (centripetal) acceleration ω2𝒓P due to the direction of the velocity continuously changing towards the center (see fig. 16.20b). The generalization we have made in this section is that the same 𝝎 can be used to calculate the acceleration for all the different points on one rotating object.

Relative motion of points on a rigid object

As you well know by now, the position of point B relative to point A is 𝒓B/A𝒓B𝒓A. Similarly the relative velocity and acceleration of two points A and B is defined to be

𝒗B/A𝒗B𝒗A and 𝒂B/A𝒂B𝒂A (16.34)

So, the relative velocity (as calculated relative to a fixed frame) of two points glued to one spinning rigid object is given by

𝒗B/A 𝒗B𝒗A
= 𝝎×𝒓B/O𝝎×𝒓A/O
= 𝝎×(𝒓B/O𝒓A/O)
= 𝝎×𝒓B/A,

where point O is the point in the Newtonian frame on the fixed axis of rotation and 𝝎=𝝎 is the angular velocity of . Repeating,

𝒗B/A=𝝎×𝒓B/A. (16.35)

Because points A and B are fixed on their velocities and hence their relative velocity as observed in a reference frame fixed to 𝒞 are all 𝟎. But, point A has an absolute velocity that is different from that of point B. So they have a relative velocity as seen in the fixed frame. And it is what you would get if B was just going in circles around A.

Filename:pfigure-s95f3a
Figure 16.21: The acceleration of B relative to A if they are both on the same rotating rigid object.

Similarly, the relative acceleration of two points glued to one rigid object spinning at constant rate is

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

Again, the relative acceleration is due to the difference in the points’ positions relative to the point O fixed on the axis. These kinematics results, 16.2 and 16.36, are useful for calculating angular momentum relative to the center-of-mass. They are also sometimes useful for the understanding of the motions of machines with moving connected parts.

The fundamental 𝝎 equation

Equation 16.35 actually applies to any vector 𝑸 that has the property that it is fixed in the frame which is rotating at 𝝎 (fig. 16.22).

Filename:Danef94s1q2
Figure 16.22: A vector 𝑸 is fixed to a rotating object. So its rate of change is the velocity of its tip. The magnitude is θ˙|𝑸| and the direction is orthogonal to 𝑸.

That is, all we needed in the derivation was that 𝑸’s components (Qx,Qy) in the ıˆ-ȷˆ reference frame be constant.

Examples of such 𝑸 would include more than just the relative position of two points fixed on the object, 𝒓B/A. Both ıˆ and ȷˆ are also are fixed in the rotating frame. Also imagine a bug moving at a constant speed on a straight line marked on a body. That bug’s velocity relative to the rotating frame would also be constant as seen in the frame. For any 𝑸 whose representation as an arrow moves like a line segment drawn on a rigid object rotating at 𝝎:

𝑸˙=𝝎×𝑸. (16.37)

This is the most fundamental equation of rigid-object kinematics. Know it well. For future reference, equation (16.37) is also valid in three dimensions. Equation 16.37 is the generalization of the three equations

𝒓˙=𝝎×𝒓,ıˆ˙=𝝎×ıˆandȷˆ˙=𝝎×ȷˆ.

Calculating relative velocity directly, using rotating frames

A coordinate system xy attached to a rotating rigid object 𝒞, defines a reference frame 𝒞 (fig. 16.6 on page 16.6).

Recall, the base vectors in this frame change in time just like any other vector fixed in the rotating frame (eqn. (16.37))

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

Once we know how the base vectors ıˆ and ȷˆ change with time, as per the above equations, we can find velocity and acceleration by differentiation of position.

For example, consider two points A and B fixed on a rotating object. The position of B relative to A is

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

where we are using x and y as a shorthand notation for xB/A and yB/A (see fig. 16.23).

Filename:pfigure-blue-123-1
Figure 16.23: Coordinates in a rotating frame.

The coordinates x and y are constant so

x˙=0 and y˙=0.

Now we just differentiate the relative position 𝒓B/A with respect to time,

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

We could similarly calculate 𝒂B/A by taking another derivative to get, after a calculation much like that above,

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

For the above calculation the points A and B were fixed on a rotating part. Most machines have many parts, at least some of which move in more complex ways than just circular motion. We will be able to understand such machines by considering points and parts that are moving relative to a part that is itself rotating. Such will be considered in Chapter 19.

Box 16.5 Angular velocity 𝝎 and the rotation matrix [R]

The rotation matrix [R] and the angular velocity vector 𝝎=ω𝒌ˆ are related. Because angular velocity is a rotation rate, in some sense it must be the derivative of the rotation. But the situation is a bit subtle because 𝝎 is a vector and [R] is a matrix.

On the one hand we have the rotation equation:

[xy][𝒓]xy=[cosθsinθsinθcosθ][R][xy][𝒓ref]xy.

This gives the coordinates of 𝒓 after rotation in terms of its coordinates before the rotation. Because the coordinates [𝒓ref]xy before rotation are fixed we can take the time derivative of both sides to get

[x˙y˙][𝒓˙]xy=θ˙[sinθcosθcosθsinθ][R˙][xy][𝒓ref]xy.

On the other hand we have the vector equation

𝒓˙=𝝎×𝒓.

where 𝝎=ω𝒌ˆ. We can relate this vector equation to the matrix component equation above it by expressing the cross product with a matrix multiplication as described in box 1.3 on 1.3. Applying that result to this 2D case where

[𝝎]xyz=[00ω]and[𝒓]xyz=[xy0]

we get

[ωyωx0][𝒓˙]xy=[0ω0ω00000][𝒮(𝝎)][xy0][𝒓ref]xy.

Just keeping the 2D terms we have

[𝝎×𝒓]xy = [𝒮(𝝎)][xy]
= [0ωω0][xy].

But [R] tells us the coordinates in terms of the reference coordinates. So

[𝝎×𝒓]xy=[0ωω0][R][xy][xy].

Comparing the equations above we see that

[R˙]=[0ωω0][R]. (16.40)

Using the notation from box 1.3 we can rewrite this as

[R˙]=[𝒮(𝝎)][R]. (16.41)

Equation 16.41 is the way that angular velocity 𝝎 and the rate of change [R˙] of the rotation matrix are related.

We can also solve for [𝒮(𝝎)] as

[𝒮(𝝎)]=[R˙][R]1

from which we can find ±ω in the off diagonal elements if we are given [R] and [R˙].


We can understand eqn. (16.41) better, perhaps, by thinking of the rotation matrix [R] as having two columns which are, respectively, the fixed-coordinate-system components of ıˆ and ȷˆ:

[R]=[[ıˆ]xy|[ȷˆ]xy]

Thus the columns of [R˙] are the rates of change of these two unit vectors:

[R˙]=[[ıˆ˙]xy|[ȷˆ˙]xy]=[[𝝎×ıˆ]xy|[𝝎×ȷˆ]xy].

With a total abuse of notation (in general one does not allow cross products of a vectors with matrices) we can write the above equation more memorably as

[R˙]=𝝎×[R].

Thus eqn. (16.41) is just a version of the fundamental kinematic equation 16.37 on page 16.37.

SAMPLE 16.3

Filename:pfigure-blue-119-2
Figure 16.24:

A uniform bar AB of length =50cm rotates counterclockwise about point A with constant angular speed ω. At the instant shown in fig. 16.24 the linear speed vC of the center-of-mass C is 7.5cm/s.

  1. 1.

    What are the angular speed and angular velocity of the bar?

  2. 2.

    What is the linear velocity of point B?

  3. 3.

    By what angles do the angular positions of points C and B change in 2 seconds?

Solution

Filename:pfigure-s95q14
Figure 16.25:

Let the angular velocity of the bar be 𝝎=θ˙𝒌ˆ where θ˙ is the angular speed. We first need to find θ˙.

  1. 1.

    The linear speed of point C is given, vC=7.5cm/s. Now,

    vC = θ˙rC
     θ˙ = vCrC=7.5cm/s25cm=0.3rad/s.

    Therefore, the angular velocity of the bar is 𝝎=θ˙𝒌ˆ=0.3rad/s𝒌ˆ.

    Answer: θ˙=0.3rad/s, and 𝝎=0.3rad/s𝒌ˆ

  2. 2.

    Point B is at distance from the pivot point A. Thus it goes around a circle of radius (see fig. 16.26). Therefore,

    Filename:Danef94s3q2
    Figure 16.26: From the given geometry, 𝒆ˆr=cosθıˆ+sinθȷˆ, 𝒆ˆθ=sinθıˆ+cosθȷˆ, and 𝒗B=|𝒗B|𝒆ˆθ.
    𝒗B = 𝝎×𝒓B=θ˙𝒌ˆ×(cosθıˆ+sinθȷˆ)
    = θ˙(cosθȷˆsinθıˆ)
    = 0.3rad/s50cm(32ȷˆ12ıˆ)
    = 15cm/s(32ȷˆ12ıˆ).

    Answer: 𝒗B=15cm/s(32ȷˆ12ıˆ)

    We can also write 𝒗B=15cm/s𝒆ˆθ where 𝒆ˆθ=32ȷˆ12ıˆ.

  3. 3.

    Let θ1 be the position of point C at some time t1 and θ2 be the position at time t2. We want to find Δθ=θ2θ1 for t2t1=2s.

    dθdt = θ˙=constant=0.3rad/s.
     dθ = (0.3rad/s)dt.
     θ1θ2𝑑θ = t1t2(0.3rad/s)𝑑t.
     θ2θ1 = 0.3rad/s(t2t1)
    orΔθ = 0.3rads2s=0.6rad.

    The change in angular position of point B is the same as that of point C. In fact, all points on AB undergo the same change in angular position because AB is a rigid body.

    Answer: ΔθC=ΔθB=0.6rad

SAMPLE 16.4

Filename:bikefork1-ang-accel
Figure 16.27:

A flywheel of diameter 2ft is made of cast iron. To avoid extremely high stresses and cracks it is recommended that the peripheral speed not exceed 6000 to 7000 ft/min. What is the corresponding rpm rating for the wheel?


Solution

Diameter of the wheel = 2ft.
radius of wheel = 1ft.

Now,

v = ωr
 ω = vr=6000ft/min1ft
= 6000radmin1rev2πrad
= 955rpm.

Similarly, corresponding to v=7000ft/min

ω = 7000ft/min1ft
= 7000radmin1rev2πrad
= 1114rpm.

Thus the rpm rating of the wheel should read 955 – 1114 rpm.

Answer: ω=955to 1114rpm.

SAMPLE 16.5

Filename:bikefork-ang-accel
Figure 16.28:

Two gears A and B have the diameter ratio of 1:2. Gear A drives gear B. If the output at gear B is required to be 150 rpm, what should be the angular speed of the driving gear? Assume no slip at the contact point.


Solution Let C and C be the points of contact on gear A and B respectively at some instant t. Since there is no relative slip between C and C, both points must have the same linear velocity at instant t. If the velocities are the same, then the linear speeds must also be the same. Thus

Filename:bikefork1-alt
Figure 16.29:
vC = vC
 ωArA = ωBrB
 ωA = ωBrBrA
= ωB2rr=2ωB
= (2)(150rpm)
= 300rpm.

Answer: ωA=300rpm

SAMPLE 16.6

Filename:bikefork-alt
Figure 16.30:

A uniform rigid rod AB of length =0.6m is connected to two rigid links OA and OB. The assembly rotates at a constant rate about point O in the xy plane. At the instant shown, when rod AB is vertical, the velocities of points A and B are 𝒗A=4.64m/sȷˆ1.87m/sıˆ, and 𝒗B=1.87m/sıˆ4.64m/sȷˆ. Find the angular velocity of bar AB. What is the length R of the links?


Solution

Let the angular velocity of the rod AB be 𝝎=ω𝒌ˆ. margin: We know that the rod rotates about the z-axis but we do not know the sense of the rotation i.e., +𝐤ˆ or 𝒌ˆ. Here we have assumed that 𝝎 is in the positive 𝒌ˆ direction, although just by sketching 𝒗A we can easily see that 𝝎 must be in the 𝒌ˆ direction.   Since we are given the velocities of two points on the rod we can use the relative velocity formula to find 𝝎:

𝒗B/A=𝝎×𝒓B/A = 𝒗B𝒗A
orω𝒌ˆ𝝎×ȷˆ𝒓B/A = (1.87ıˆ4.64ȷˆ)m/s(4.64ȷˆ1.87ıˆ)m/s
orω(ıˆ) = (1.87ıˆ+1.87ıˆ)m/s(4.64ȷˆ4.64ȷˆ)m/s
= 3.74ıˆm/s
 ω = 3.74m/s (16.42)
= 3.740.6rad/s
= 6.233rad/s
Thus,
𝝎 = 6.233rad/s𝒌ˆ. (16.43)

Answer: 𝝎=6.23rad/s𝒌ˆ

Let θ be the angle between link OA and the horizontal axis. Now,

Filename:pfigure-blue-128-1
Figure 16.31: The given velocity vector 𝒗A. Since 𝒗A=𝝎×𝒓A where 𝒓A=Rcosθıˆ+Rsinθȷˆ, we can find R from the given 𝒗A.
𝒗A=𝝎×𝒓A = ω𝒌ˆ×R(cosθıˆsinθȷˆ)𝒓A
or(4.64ȷˆ1.87ıˆ)m/s = ωR(cosθȷˆ+sinθıˆ)

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

1.87m/s = ωRsinθ (16.44)
4.64m/s = ωRcosθ (16.45)

Squaring and adding Eqns (16.44) and (16.45) together we get

ω2R2 = (4.64m/s)2+(1.187m/s)2
= 25.026m2/s2
 R2 = 25.026m2/s2(6.23rad/s)2
= 0.645m2
 R = 0.8m

Answer: R=0.8m

SAMPLE 16.7

Filename:tfigure8-syst-bods
Figure 16.32: An ‘L’ shaped bar rotates at speed ω about point O.

Verify the relative velocity formula: The motor at O in Fig. 16.32 rotates the ‘L’ shaped bar OAB in counterclockwise direction at an angular speed which increases at ω˙=2.5rad/s2. At the instant shown, the angular speed ω=4.5rad/s. Each arm of the bar is of length L = 2ft.

  1. 1.

    Find the velocity of point A.

  2. 2.

    Find the relative velocity 𝒗B/A(=𝝎×𝒓B/A) and use the result to find the absolute velocity of point B (𝒗B=𝒗A+𝒗B/A).

  3. 3.

    Find the velocity of point B directly. Check the answer obtained in part (b) against the new answer.


Solution

  1. 1.
    Filename:sfig8-7-2
    Figure 16.33: 𝒗A=𝝎×𝒓A is tangential to the circular path of point A.

    As the bar rotates, every point on the bar goes in circles centered at point O. Therefore, we can easily find the velocity of any point on the bar using circular motion formula 𝒗=𝝎×𝒓. Thus,

    𝒗A = 𝝎×𝒓A=ω𝒌ˆ×Lıˆ=ωLȷˆ
    = 4.5rad/s2ftȷˆ=9ft/sȷˆ.

    The velocity vector 𝒗A is shown in Fig. 16.33.

    Answer: 𝒗A=9ft/sȷˆ

  2. 2.
    Filename:sfig8-7-2a
    Figure 16.34: 𝒗B/A=𝝎×𝒓B/A and 𝒗B=𝒗A+𝒗B/A.

    Point B and A are on the same rigid body. Therefore, with respect to point A, point B goes in circles about A. Hence the relative velocity of B with respect to A is

    𝒗B/A = 𝝎×𝒓B/A
    = ω𝒌ˆ×Lȷˆ=ωLıˆ
    = 4.5rad/s2ftıˆ=9ft/sıˆ.
    and 𝒗B = 𝒗A+𝒗B/A
    = 9ft/s(ıˆ+ȷˆ).

    These velocities are shown in Fig. 16.34.

    Answer: 𝒗B/A=9ft/sıˆ,𝒗B=9ft/s(ıˆ+ȷˆ)

  3. 3.
    Filename:sfig8-7-2again
    Figure 16.35: 𝒗B=𝝎×𝒓B.

    Since point B goes in circles of radius OB about point O, we can find its velocity directly using circular motion formula:

    𝒗B = 𝝎×𝒓B
    = ω𝒌ˆ×(Lıˆ+Lȷˆ)=ωL(ȷˆıˆ)
    = 9ft/s(ıˆ+ȷˆ).

    The velocity vector is shown in Fig. 16.35. Of course this velocity is the same velocity as obtained in part (b) above.

    Answer: 𝒗B=9ft/s(ıˆ+ȷˆ)

Note: Nothing in this sample uses ω˙!

SAMPLE 16.8

Filename:sfig8-7-2disks
Figure 16.36: The ‘L’ shaped bar rotates counterclockwise while slowing down.

Verify the relative acceleration formula: Consider the ‘L’ shaped bar of Sample 16.31 again. At the instant shown, the bar is rotating at 4rad/s and is slowing down at the rate of 2rad/s2.

  1. (i)

    Find the acceleration of point A.

  2. (ii)

    Find the relative acceleration 𝒂B/A of point B with respect to point A and use the result to find the absolute acceleration of point B (𝒂B=𝒂A+𝒂B/A).

  3. (iii)

    Find the acceleration of point B directly and verify the result obtained in (ii).

Solution

We are given:

𝝎=ω𝒌ˆ=4rad/s𝒌ˆ, and 𝝎˙=ω˙𝒌ˆ=2rad/s2𝒌ˆ.
  1. Filename:sfig8-4-4
    Figure 16.37:
  2. (i)

    Point A is going in circles of radius L. Hence,

    𝒂A = 𝝎˙×𝒓A+𝝎×(𝝎×𝒓A)=𝝎˙×𝒓Aω2𝒓A
    = ω˙𝒌ˆ×Lıˆω2Lıˆ=ω˙Lȷˆω2Lıˆ
    = 2rad/s2ftȷˆ(4rad/s)22ftıˆ
    = (4ȷˆ+32ıˆ)ft/s2.

    Answer: 𝒂A=(4ȷˆ+32ıˆ)ft/s2

  3. (ii)
    Filename:sfig8-4-4a
    Figure 16.38:

    The relative acceleration of point B with respect to point A is found by considering the motion of B with respect to A. Since both the points are on the same rigid body, point B executes circular motion with respect to point A. Therefore,

    𝒂B/A = 𝝎˙×𝒓B/A+𝝎×(𝝎×𝒓B/A)=𝝎˙×𝒓B/Aω2Lȷˆ
    = ω˙𝒌ˆ×Lȷˆω2Lȷˆ
    = ω˙Lıˆω2Lȷˆ=2rad/s22ftıˆ(4rad/s)22ftȷˆ
    = (4ıˆ32ȷˆ)ft/s2,

    and

    𝒂B = 𝒂A+𝒂B/A=(28ıˆ36ȷˆ)ft/s2.

    Answer: 𝒂B=(28ıˆ+36ȷˆ)ft/s2

  4. (iii)
    Filename:sfig8-4-4b
    Figure 16.39:

    Since point B is going in circles of radius OB about point O, we can find the acceleration of B as follows.

    𝒂B = 𝝎˙×𝒓B+𝝎×(𝝎×𝒓B)
    = 𝝎˙×𝒓Bω2𝒓B
    = ω˙𝒌ˆ×(Lıˆ+Lȷˆ)ω2(Lıˆ+Lȷˆ)
    = (ω˙Lω2L)ȷˆ+(ω˙Lω2L)ıˆ
    = (432)ft/s2ȷˆ+(432)ft/s2ıˆ
    = (36ȷˆ28ıˆ)ft/s2.

    This acceleration is, naturally again, the same acceleration as found in (ii) above.

    Answer: 𝒂B=(28ıˆ+36ȷˆ)ft/s2

Problems for 16.2 Angular velocity

Preparatory Problems

16.2.1   Give the definition of each term below in words and, but for the first, with an equation.

  1. (a)

    angular velocity Answer: The rate of change of any marked or imaginary line on a rigid object.

  2. (b)

    angular acceleration Answer: The rate of change of the angular velocity.

16.2.2  Assume that in the reference configuration, when θ=0, that the xy axes attached to an object are aligned with the xy axes. Consider a point P attached to the moving frame (object) that has coordinates (x,y). Find each of the quantities below in terms of some or all of x,y,ıˆ,ȷˆ,θ,θ˙ and θ¨.

  1. (a)

    𝒓P Answer: [cosθsinθsinθcosθ][xy]

  2. (b)

    𝒗P Answer: vP=ω×rP

  3. (c)

    𝒂P Answer: aP=ω˙×rP+ω×(ω×rP)

16.2.3  Assume that in the reference configuration, when θ=0, the xy axes are aligned with the xy axes. Consider two points P1 and P2 attached to the moving frame. P1 and P2 have coordinates (x1,y1) (x2,y2). Find each of the quantities below in terms of some or all of x1,y1,x2,y2,ıˆ,ȷˆ,θ,θ˙ and θ¨.

  1. (a)

    𝒓P1/P2 Answer: rP1/P2=(x2x1)ıˆ+(y2y1)ȷˆ

  2. (b)

    𝒗P1/P2 Answer: vP1/P2=θ˙((x2x1)ȷˆ(y2y1)ıˆ)

  3. (c)

    𝒂P1/P2 Answer: aP1/P2=θ¨((x2x1)ȷˆ(y2y1)ıˆ)θ˙2((x2x1)ıˆ+(y2y1)ȷˆ)

16.2.4   Find 𝒗=𝝎×𝒓, if 𝝎=1.5rad/s𝒌ˆ and 𝒓=2mıˆ3mȷˆ. Answer: (5.4ıˆ+3.0ȷˆ)m/s

16.2.5  A rod OB rotates with its end O fixed as shown in the figure with angular velocity 𝝎=5rad/s𝒌ˆ and angular acceleration 𝜶=2rad/s2𝒌ˆ at the instant of interest. Find, draw, and label the tangential and normal acceleration of point B at θ=60. Answer: at=20ȷˆin/s2 and 𝒂n=250ıˆin/s2.

Filename:sfig8-4-4c
Figure 16.40:

16.2.6  A disc rotates at 15  rpm. How many seconds does it take to rotate by 180 degrees? What is the angular speed of the disc in  rad/s? Answer: rad/s=π2rad/s.

16.2.7   A motor turns a uniform disc of radius R counter-clockwise about its mass center at a constant rate ω. The disc lies in the xy-plane and its angular displacement θ is measured (positive counter-clockwise) from the x-axis. What is the angular displacement θ(t) of the disc if it starts at θ(0)=θ0 and θ˙(0)=ω? What are the velocity and acceleration of a point P at position 𝒓=xıˆ+yȷˆ? Answer: θ(t)=θ0+ωt, 𝒗=ω(xȷˆyıˆ), and 𝒂=ω(xıˆ+yȷˆ).

16.2.8  Find the angular velocities of the second, minute, and hour hands of a clock. Answer: ωsec=π30rad/s, ωmin=2π3600rad/s, and ωhour=2π43200rad/s.

16.2.9   A disc 𝒞 spins at a constant rate of two revolutions per second counter-clockwise about its geometric center, G, which is fixed. A point P is marked on the disk at a radius of one meter. At the instant of interest, point P is on the x-axis of an xy-coordinate system centered at point G.

  1. (a)

    Draw a neat diagram showing the disk, the particle, and the coordinate axes.

  2. (b)

    What is the angular velocity of the disk, 𝝎𝒞? Answer: ω𝒞=4πrad/s.

  3. (c)

    What is the angular acceleration of the disk, 𝝎˙C? Answer: ω˙C=0rad/s2.

  4. (d)

    What is the velocity 𝒗P of point P? Answer: vp=4πȷˆm/s.

  5. (e)

    What is the acceleration 𝒂P of point P? Answer: ap=16π2ıˆm/s2.

16.2.10  A uniform rigid rod rotates at constant speed in the xy-plane about a peg at point O. The center of mass of the rod may not exceed a specified acceleration amax=0.5m/s2. Find the maximum angular velocity of the rod. Answer: ωmax=1rad/s

Filename:pfigure-s94h13p2
Figure 16.41:

16.2.11  A dumbbell AB is welded to a rigid arm OC such that OC is perpendicular to AB. Arm OC rotates about O at a constant angular velocity 𝝎=10rad/s𝒌ˆ. At the instant when θ=60,

  1. (a)

    Find the relative velocity of B with respect to A. Answer: vB/A=4m/s(ıˆ+3ȷˆ)

  2. (b)

    Find the acceleration of point B relative to point A. Answer: aB/A=40m/s2(3ıˆȷˆ)

Filename:pfigure-f93f5
Figure 16.42:

More-Involved Problems

16.2.12  Two discs 𝒜 and rotate at constant speeds about their centers. Disc 𝒜 rotates at 100rpm and disc rotates at 10rad/s. Which is rotating faster? Answer: Disc 𝒜 is faster.

16.2.13  A motor turns a uniform disc of radius R counter-clockwise about its mass center at a constant rate ω. The disc lies in the xy-plane and its angular displacement θ is measured (positive counter-clockwise) from the x-axis. What are the velocity and acceleration of a point P at position 𝒓P=cıˆ+dȷˆ relative to the velocity and acceleration of a point Q at position 𝒓Q=0.5(dıˆ+cȷˆ) on the disk? (c2+d2<R2.)

Answer: vP/Q=ω((c+d2)ȷˆ(dc2)ıˆ), and 𝒂P/Q=ω2((c+d2)ıˆ+(dc2)ȷˆ)

16.2.14  A 0.4m long rod AB has many holes along its length such that it can be pegged at any of the various locations. It rotates counter-clockwise at a constant angular speed about a peg whose location is not known. At some instant t, the velocity of end B is 𝒗B=3m/sȷˆ. After π20s, the velocity of end B is 𝒗B=3m/sıˆ. If the rod has not completed one revolution during this period,

  1. (a)

    find the angular velocity of the rod, and Answer: ω=30rad/s.

  2. (b)

    find the location of the peg along the length of the rod. Answer: 110m from B.

Filename:pfigure-s94h13p3
Figure 16.43:

16.2.15  A circular disc of radius r=250mm rotates in the xy-plane about a point which is at a distance d=2r away from the center of the disk. At the instant of interest, the linear speed of the center C is 0.60m/s and the magnitude of its centripetal acceleration is 0.72m/s2.

  1. (a)

    Find the rotational speed of the disk. Answer: θ˙=1.2rad/s.

  2. (b)

    Is the given information enough to locate the center of rotation of the disk? Answer: No.

  3. (c)

    If the acceleration of the center has no component in the ȷˆ direction at the instant of interest, can you locate the center of rotation? If yes, is the point you locate unique? If not, what other information is required to make the point unique? Answer: If the acceleration of the center has no component in the ȷˆ direction the center of rotation can lie anywhere on the x-axis.

Filename:p-s96-p3-3
Figure 16.44:

16.2.16  A uniform disc of radius r=200mm is mounted eccentrically on a motor shaft at point O. The motor rotates the disc at a constant angular speed. At the instant shown, the velocity of the center of mass is 𝒗G=1.5m/sȷˆ.

  1. (a)

    Find the angular velocity of the disc. Answer: ω=11.25rad/skˆ

  2. (b)

    Find the point with the highest linear speed on the disc. What is its velocity? Answer: Speed is maximum at 𝒓=5r3ıˆ; 𝒗max=3.75m/sȷˆ

Filename:bikefork1-ang-mom
Figure 16.45:

16.2.17  The circular disc of radius R=100mm rotates about its center O. At a given instant, point A on the disk has a velocity vA=0.8m/s in the direction shown. At the same instant, the tangent of the angle θ made by the total acceleration vector of any point B with its radial line to O is 0.6. Compute the angular acceleration α of the disc. Answer: α=38.4rad/s2

Filename:bikefork-ang-mom
Figure 16.46:

16.2.18  Show that, for non-constant rate circular motion, the acceleration of all points in a given radial line are parallel.

Filename:summer95f-5-a
Figure 16.47:

16.2.19  A motor turns a uniform disc of radius R about its mass center at a variable angular rate ω with rate of change ω˙, counter-clockwise. The disc lies in the xy-plane and its angular displacement θ is measured from the x-axis, positive counter-clockwise. What are the velocity and acceleration of a point P at position 𝒓P=cıˆ+dȷˆ relative to the velocity and acceleration of a point Q at position 𝒓Q=0.5(dıˆ+cȷˆ) on the disk? (c2+d2<R2.) Answer: vP/Q=ω((c+d2)ȷˆ(dc2)ıˆ), and
𝒂P/Q=(ω˙(dc2)+ω2(c+d2))ıˆ+(ω˙(c+d2)ω2(dc2))ȷˆ

16.2.20  The dumbbell AB shown in the figure rotates counterclockwise about point O with angular acceleration 3rad/s2. Bar AB is perpendicular to bar OC. At the instant of interest, θ=45 and the angular speed is 2  rad/s.

  1. (a)

    Find the velocity of point B relative to point A. Will this relative velocity be different if the dumbbell were rotating at a constant rate of 2  rad/s? Answer: 12m/s(ıˆ+ȷˆ)

  2. (b)

    Without calculations, draw a vector approximately representing the acceleration of B relative to A.

  3. (c)

    Find the acceleration of point B relative to A. What can you say about the direction of this vector as the motion progresses in time? Answer: aB/A=(0.353ıˆ+2.474ȷˆ)m/s2

Filename:pfigure4-2-rp10
Figure 16.48:

16.2.21  Bit-stream kinematics of a CD. A Compact Disk (CD) has bits of data etched on concentric circular tracks. The data from a track is read by a beam of light from a head that is positioned under the track. The angular speed of the disk remains constant as long as the head is positioned over a particular track. As the head moves to the next track, the angular speed of the disk changes, so that the linear speed at any track is always the same. The data stream comes out at a constant rate 4.32×106 bits/second. When the head is positioned on the outermost track, for which r=56mm, the disk rotates at 200rpm.

  1. (a)

    What is the number of bits of data on the outermost track? Answer: 1.296 Megabits.

  2. (b)

    find the angular speed of the disk when the head is on the innermost track (r=22mm), and Answer: 509.09 rpm.

  3. (c)

    find the number of bits on the innermost track. Answer: 0.509 Megabits

16.2.22  2-D constant rate gear train. The angular velocity of the input shaft (driven by a motor not shown) is a constant, ωinput=ωA. What is the angular velocity ωoutput=ωC of the output shaft and the speed of a point on the outer edge of disc C, in terms of RA, RB, RC, and ωA? Answer: ωC=ωARARB and vC=ωARARCRB

Filename:pfigure-blue-125-2
Figure 16.49:

Gear B is welded to C and engages with A.

16.2.23  A horizontal disk 𝒟 of diameter d=500mm is driven at a constant speed of 100rpm. A small disk 𝒞 can be positioned anywhere between r=10mm and r=240mm on disk 𝒟 by sliding it along the overhead shaft and then fixing it at the desired position with a set screw (see the figure). Disk 𝒞 rolls without slip on disk 𝒟. The overhead shaft rotates with disk 𝒞 and, therefore, its rotational speed can be varied by varying the position of disk 𝒞. This gear system is called brush gearing. Find the maximum and minimum rotational speeds of the overhead shaft. Answer: ωmin=10rpm and ωmax=240rpm

Filename:pfigure-blue-68-1
Figure 16.50:

16.2.24   Two points A and B are on the same machine part that is hinged at an as yet unknown location C. Assume you are given that points at positions 𝒓A and 𝒓B are supposed to move in given directions, indicated by unit vectors 𝝀ˆA and 𝝀ˆB. For each of the problem parts below, illustrate your results with two numerical examples (in consistent units): i) 𝒓A=1ıˆ, 𝒓B=1ȷˆ, 𝝀ˆA=1ȷˆ, and 𝝀ˆB=1ıˆ (thus 𝒓C=𝟎), and ii) a more complex example of your choosing.

  1. (a)

    Describe in detail what equations must be satisfied by the point 𝒓C.

  2. (b)

    Write a computer program that takes as input the 4 pairs of numbers [𝒓A], [𝒓B], [𝝀ˆA] and [𝝀ˆB] and gives as output the pair of numbers [𝒓C].

  3. (c)

    Find a formula of the form 𝒓C= that explicitly gives the position vector for point C in terms of the 4 given vectors. Answer: rC=rB+(rArB)λˆA((rC×λˆB)×λˆB)λˆA((rC×λˆB)×λˆB)

16.3 Polar moment of inertia

Filename:pfigure-blue-58-1
Figure 16.51: A bit of mass dm on a general planar object.

All of the basic equations of Dynamics concern the motion of mass. In particular, linear momentum, angular momentum, and energy are all expressed as sums over all the bits of mass in a system, with each bit of mass multiplied by some terms concerning position, velocity and acceleration.

From the earlier sections in this chapter we know how to find the velocity and acceleration of every bit of mass on a 2-D rigid object as it spins about a fixed axis. So it is just a matter of doing integrals or sums to calculate the various momentum and energy quantities of interest. As an object moves and rotates the region of integration and the values of the integrands change. So, in principle, in order to analyze a rigid object one has to evaluate a different integral or sum at every different configuration. But there is a shortcut: for a rotating rigid object a sum (over all atoms, say), or a difficult integral (for example, over the complex region representing a machine part) is reduced to simple multiplication.

The moment of inertia

Icmmargin: The COM is the most common reference point. The moment of inertia of a given object depends on what reference point is used. Most often when people say ‘the’ moment of inertia they mean the moment of inertia with respect to the center-of-mass. For clarity, this moment of inertia is often notated as Icm in this book. If a different reference point, say point O is used, the moment of inertia is notated as I0.

simplifies the expressions for the angular momentum, the rate of change of angular momentum, and the energy of a rigid object. For more general motions the shortcuts need a 3×3 matrix [𝑰cm] But for 2D mechanics only one component of the matrix [𝑰cm] is relevant, it is Izzcm, called just I or J for short.

Here are the main results. A flat object spinning with 𝝎=ω𝒌ˆ in the xy plane has a mass distribution which gives, by means of a calculation which we will discuss shortly, a moment of inertia Icm so that:

𝑯cm = Icmω𝒌ˆ (16.46)
𝑯˙cm = Icmω˙𝒌ˆ (16.47)
EK/cm = 12ω2Icm. (16.48)

There are two main skills you need to develop associated with I.

  • Using the formulas above properly in angular momentum and energy balance. This is covered in sec. 16.4.

  • Finding I of an object, as is covered in most freshman calculus texts and is reviewed in this section.

Some facts about moment of inertia are summarized on the inside back cover.

The moments of inertia in 2-D : [𝑰cm] and [𝑰O].

The definition of moment of inertia

margin: What do the words ‘moment of inertia’ mean? The word inertia, in this case, can be translated to mean mass. So it’s ‘moment of mass’. The word ‘moment’ is similar to the word ‘moment’ as in torque, or (force)(distance), but slightly more general. The more general word moment is ‘the product of a quantity and its distance, raised to a power, from a reference point,’ the n’th moment of a quantity rel to 0 =xn(quantity at x)𝑑x. So, what we used to call ‘moment’ (meaning torque) we can think of as short for ‘the first moment of perpendicular force’. And the moment of inertia is ‘the second moment of mass’. The words ‘first’ and ‘second’ refer to the exponent n.

Icm is

Icm (x2+y2r2)𝑑m
=r2(mtotA)The mass per unit area.dAdmfor a uniform planar object

where x=x/cm, y=y/cm and r=r/cm are the distances in the x and y directions of the bits of mass (dm) from the center of mass and r=r/cm is the total distance. Similarly if all distances are measured relative to a point C or O (instead of relative to the center of mass) then similar formulas as above calculate IC or IO. See box 16.3 on page 16.3 for the calculation of the moments of inertia of some simple objects.

The term Icm=Izzcm is sometimes called the polar moment of inertia, or polar mass moment of inertia to distinguish it from the Ixx and Iyy, terms which have little utility in planar dynamics (Ixx and Iyy, as area inertias rather than mass inertias, are all-important when calculating the bending stiffness or the stress in elastic beams, however).

What, physically, is the polar moment of inertia? It is a measure of the extent to which mass is far from the given reference point. Every bit of mass contributes to I in proportion to the square of its distance from the reference point. As we will show in sec. 16.4 on page 16.4 (see eqn. (16.51) on page 16.51) Icm is just the quantity we need to do mechanics problems.

Radius of gyration

Another measure of the extent to which mass is spread from the reference point, besides the moment of inertia Icm, is the radius of gyration, rgyr

margin: k for stiffness. The radius of gyration rgyr is sometimes called k but we save k for stiffness in this book.

.

The radius of gyration is defined as:

rgyrI/m  rgyr2m=I.

That is, the radius of gyration of an object is the radius of an equivalent ring of mass that has the same I and the same mass as the given object. It is an equivalent single radius for the whole mass distribution

margin: Radius of gyration and average distance of mass. The radius of gyration is an equivalent radius of the collection of mass. But the radius of gyration is not exactly the average radius of the distributed mass. Rather the radius of gyration is the square root of the mean square radius: rgyrr2𝑑mmtot. For the special cases where all the mass is at the same radius R (as for a circular hoop or for masses distributed on the vertices of, say, a regular pentagon) the radius of gyration is R. For all other cases the radius of gyration is larger than the average radius of the mass.

.

Example: Radius of gyration of a hoop

The radius of gyration about its center of a circular hoop of mass m and radius R about its center is rgyr=R.

Other reference points

For the most part it is Icm which is of primary interest. Other reference points are useful

margin: Notation: simple vs precise. When just the symbol I is used one assumes I=Izzcm, and if just Icm is used, one assumes Icm=Izzcm. Similarly I0 and IC are assumed to mean Izz0 and IzzC, respectively.
  1. 1.

    If the rigid object is hinged at a fixed point O then I0 can be used in a slight short cut in calculation of angular momentum and energy; and

  2. 2.

    If one wants to calculate the moment of inertia of a composite object about its center-of-mass then it is useful to first find the moment of inertia of each of its parts about that system’s center of mass (which is generally not the center-of-mass of any of the separate parts).

  3. 3.

    Sometimes it is easiest to set up the integral for moment of inertia about a special point C that is not the center of mass and then use the parallel axis theorem (eqn. (16.49)) to find the moment of inertia about the center of mass (see, for example the semi-circle example on 16.3).

The parallel axis theorem

The planar parallel axis theorem is the equation

IzzC=Izzcm+mtotrcm/C2d2. (16.49)

In this equation d=rcm/C is the distance from the center-of-mass to a line parallel to the z-axis which passes through point C. See box 16.3 on page  16.3 for a derivation of the parallel axis theorem for planar objects.

Note that IzzCIzzcm, always.

Moment of inertia of complex objects.  One can calculate the moment of inertia of a composite object about its center of mass, in terms of the masses and moments of inertia of the separate parts. Say the position of the center of mass of mi is (xi,yi) relative to a fixed origin, and the moment of inertia of that part about its center of mass is Ii. We can then find the moment of inertia of the composite Itot about its center-of-mass (xcm,ycm) by the following sequence of calculations, in order:

(1)mtot=mi(2)xcm=[ximi]/mtotycm=[yimi]/mtot(3)di2=(xixcm)2+(yiycm)2(4)Itot=[Iicm+midi2].

You can reduce this recipe to one grand formula with lots of summation signs. But you would end up doing the calculations in about the order prescribed above in any case. This sequence of steps lends itself naturally to a computer spread sheet or to any program that deals easily with arrays of numbers. This is laid out for the similar center-of-mass calculation on page 2.2

margin: The tidy recipe just presented is actually more commonly used, with slight modification, in strength of materials than in dynamics. The need for finding area moments of inertia of strange beam cross sections arises more frequently than the need to find the polar mass moment of inertia of a strange cutout shape.

.

Example: Inertia of a semicircle about its COM.

The moment of inertia of a semicircle with mass m and radius r about O (see fig. 16.52) is

Filename:pfigure-blue-157-1
Figure 16.52: The moment of inertia of a semi-circle about the COM is calculated by calculating about point O and then using the parallel axis theorem.
I0=r/02𝑑m=mr2.

From the example on page 2.2 the center of mass is at

yG=2πr.

By the parallel axis theorem

I0=Icm+myG2

which we can solve for Icm to get

Icm=I0myG2=mr24π2mr2=(14π2)mr20.6mr2.

The radius of gyration is that radius ρgyr so that mρgyr2=Icm and is

ρgyr=14π2r.

The perpendicular axis theorem for planar rigid bodies

The perpendicular axis theorem for planar objects is the equation

Izz=Ixx+Iyy

which is derived in box 16.3 on page 16.3. It gives the ‘polar’ inertia Izz in terms of the inertias Ixx and Iyy. Unlike the parallel axis theorem, the perpendicular axis theorem does not have a three-dimensional counterpart. The theorem is of greatest utility when one wants to study the three-dimensional mechanics of a flat object and thus are in need of its full moment of inertia matrix.

Box 16.6 The 2-D parallel axis theorem and the perpendicular axis theorem

Sometimes one wants to know the moment of inertia relative to the center of mass. And, sometimes, if the object is held at a hinge joint at O, relative to that hinge point O There is a simple relation between these two moments of inertia known as the parallel axis theorem.

2-D parallel axis theorem

For the two-dimensional mechanics of two-dimensional objects, our only concern is Izzo and Izzcm and not the full moment of inertia matrix. In this case, Izzo=r/o2𝑑m and Izzcm=r/cm2𝑑m. Now, let’s prove the theorem in two dimensions referring to the figure.

Filename:tfigure4-para-axis

IzzO = r/O2𝑑m
= (x/O2+y/O2)𝑑m
= [(xcm/O+x/cm)x/O2+(ycm/O+y/cm)y/O2]𝑑m
= [(xcm/O2+2xcm/Ox/cm+x/cm2)+
(ycm/O2+2ycm/Oy/cm+y/cm2)]dm
= (xcm/O2+ycm/O2)𝑑mm+2xcm/Ox/cm𝑑m0+
2ycm/Oy/cm𝑑m0+(x/cm2+y/cm2)𝑑m
= rcm/O2m+(x/cm2+y/cm2)𝑑mIzzcm
= Izzcm+rcm/O2d2m

The cancellation y/cm𝑑m=x/cm𝑑m=0 comes from the definition of center of mass.

Sometimes, people write the parallel axis theorem more simply as

I0=Icm+md2orJO=Jcm+md2

using the symbol J to mean Izz. One thing to note about the parallel axis theorem is that the moment of inertia about any point O is always greater than the moment of inertia about the center of mass. For a given object, the minimum moment of inertia is about the center-of-mass.

Why the name parallel axis theorem? We use the name because the two I’s calculated are the moments of inertia about two parallel axes (both in the z direction) through the two points cm and O.

One way to think about the theorem is the following. The moment of inertia of an object about a point O not at the center-of-mass is the same as that of the object about the cm plus that of a point mass located at the center-of-mass. If the distance from O to the cm is larger than the outer radius of the object, then the d2m term is larger than Izzcm. The distance of equality of the two terms is the radius of gyration, rgyr.

Perpendicular axis theorem (applies to planar objects only)

For planar objects,

IzzO = |𝒓|2𝑑m
= (x/O2+y/O2)𝑑m
= x/O2𝑑m+y/O2𝑑m
= IyyO+IxxO

Similarly,

Izzcm=Ixxcm+Iyycm.

That is, the moment of inertia about the z-axis is the sum of the inertias about the two perpendicular axes x and y. Note that the objects must be planar (z=0 everywhere) or the theorem would not be true. For example, Ixxo=(y/O2+z/O2)𝑑my/O2𝑑m for a three-dimensional object.

Filename:tfigure4-perp-axis-thm

Box 16.7 Some examples of 2-D Moment of Inertia

Here, we illustrate some simple moment of inertia calculations for two-dimensional objects. The needed formulas are summarized, in part, by the lower right corner components (that is, the elements in the third column and third row (3,3)) of the matrices in the table on the inside back cover.

One point mass

Filename:tfigure4-pointmassIo

If we assume that all mass is concentrated at one or more points, then the integral

Izzo=r/o2𝑑m

reduces to the sum

Izzo=ri/o2mi

which reduces to one term if there is only one mass,

Izzo=r2m=(x2+y2)m.

So, if x=3, y=4, and m=0.1lbm, then Izzo=2.5lbm2. Note that, in this case, Izzcm=0 since the radius from the center-of-mass to the center-of-mass is zero.

Two point masses

Filename:tfigure4-inertia-fig0

In this case, the sum that defines Izzo reduces to two terms, so

Izzo=ri/o2mi=m1r12+m2r22.

Note that, if r1=r2=r, then Izzo=mtotr2.

A thin uniform rod

Filename:tfigure4-2DrodIo

Consider a thin rod with uniform mass density, ρ, per unit length, and length . We calculate Izzo as

Izzo = r2dmρds
= 12s2ρ𝑑s(s=r)
= 13ρs3|12(sinceρconst.)
= 13ρ(13+23).

If either 1=0 or 2=0, then this expression reduces to Izzo=13m2. If 1=2, then O is at the center-of-mass and

Izzo=Izzcm=13ρ((2)3+(2)3)=ml212.

We can illustrate one last point. With a little bit of algebraic histrionics of the type that only hindsight can inspire, you can verify that the expression for Izz0 can be arranged as follows:

Izz0 = 13ρ(13+23)
= ρ(1+2)m(212d)2+ρ(1+2)312m2/12
= md2+m212
= md2+Izzcm

That is, the moment of inertia about point O is greater than that about the center of mass by an amount equal to the mass times the distance from the center-of-mass to point O squared. This derivation of the parallel axis theorem is for one special case, that of a uniform thin rod.

A uniform hoop

Filename:tfigure4-2DhoopIo

For a hoop of uniform mass density, ρ, per unit length, we might consider all of the points to have the same radius R. So,

Izzo=r2𝑑m=R2𝑑m=R2𝑑m=R2m.

Or, a little more tediously,

Izzo = r2𝑑m
= 02πR2ρR𝑑θ
= ρR302π𝑑θ
= 2πρR3=(2πρR)mR2=mR2.

This Izzo is the same as for a single point mass m at a distance R from the origin O. It is also the same as for two point masses if they both are a distance R from the origin. For the hoop, however, O is at the center-of-mass so Izzo=Izzcm which is not the case for a single point mass.

A uniform disk

Filename:tfigure4-2DdiskIo

Assume the disk has uniform mass density, ρ, per unit area. For a uniform disk centered at the origin, the center-of-mass is at the origin so

Izzo=Izzcm = r2𝑑m
= 0R02πr2ρr𝑑θ𝑑r
= 0R2πρr3𝑑r
= 2πρr44|0R=πρR42=(πρR2)R22
= mR22.

For example, a 1kg plate of 1m radius has the same moment of inertia as a 1kg hoop with a 70.7cm radius.

Uniform rectangular plate

Filename:tfigure4-2DplateIo

For the special case that the center of the plate is at point O, the center-of-mass is also at O and Izzo=Izzcm.

Izzo=Izzcm = r2𝑑m
= b2b2a2a2(x2+y2)ρdxdydm
= b2b2ρ(x33+xy2)|x=a2x=a2dy
= ρ(x3y3+xy33)|x=a2x=a2|y=b2y=b2
= ρ(a3b12+ab312)
= m12(a2+b2).

Note that r2𝑑m=x2𝑑m+y2𝑑m for all planar objects (the perpendicular axis theorem). For a uniform rectangle, y2𝑑m=ρy2𝑑A. But the integral y2dA is just the term often used for I, the area moment of inertia, in strength of materials calculations for the stresses and stiffnesses of beams in bending. You may recall that y2𝑑A=ab312=Ab212 for a rectangle. Similarly, x2𝑑A=Aa212. So, the polar moment of inertia J=Izzo=m112(a2+b2) can be recalled by remembering the area moment of inertia of a rectangle combined with the perpendicular axis theorem.

SAMPLE 16.9

Filename:summer95f-5
Figure 16.53:

Moment of inertia of point masses: A pendulum is made up of two unequal point masses m and 2m connected by a massless rigid rod of length 4r. The pendulum is pivoted at distance r along the rod from the small mass.

  1. 1.

    Find the moment of inertia Izzcm of the pendulum.

  2. 2.

    Find the moment of inertia IzzO of the pendulum.

  3. 3.

    Find the radius of gyration of the pendulum.


Solution

  1. 1.

    First we need to find the center of mass of the system. Let the center of mass C be located at distance rcm from the origin O. An easy way to find the location of C will be to consider the pendulum to lie along the x-axis with the origin at O (see fig. 16.54). Then

    Filename:pfigure-blue-127-2
    Figure 16.54:
    xcm = miximi
     rcm = m(r)+2m(3r)m+2m=53r.

    So, the moment of inertia about the center of mass is

    Izzcm=miri/cm2=m(r+53r)2+2m(3r53r)2=323mr2.

    Answer: Izzcm=10.67mr2

  2. 2.

    We can calculate IzzO in two different ways, one directly by summing the contributions of each mass about point O, and the other by using Izzcm and the parallel axis theorem.

    IzzO = mr2+2m(3r)2=19mr2
    or IzzO = Izzcm+mtotrcm/O2
    = 323mr2+3m(53r)2=19mr2.

    Answer: IzzO=19mr2

  3. 3.

    The radius of gyration, by definition, is the distance that gives the desired moment of inertia if the entire mass of the system is concentrated there. Thus, if rgyr is the radius of gyration for the moment of inertia about point O, then

    Filename:pfigure-s94h13p4
    Figure 16.55:
    IzzO = (3m)rgyr2
     19mr2 = 3mrgyr2
     rgyr = 193r=2.52r.

    Thus the radius of gyration rgyr of the given pendulum is rgyr=2.52r.

    Answer: rgyr=2.52r

SAMPLE 16.10

Filename:pfigure-blue-90-2
Figure 16.56:

Moment of inertia of a rod: A uniform rigid rod AB of mass M=2kg and length 3=1.5m swings about the z-axis passing through the pivot point O.

  1. 1.

    Find the moment of inertia IzzO of the bar using the fundamental definition IzzO=mr/O2𝑑m.

  2. 2.

    Find IzzO using the parallel axis theorem given that Izzcm=112m2 where m= total mass, and = total length of the rod. (You can find Izzcm for many commonly encountered objects in the table on the inside backcover of the text).

Solution

Filename:s92f1p7
Figure 16.57:
  1. 1.

    Since we need to carry out the integral, IzzO=mr/O2𝑑m, to find IzzO, let us consider an infinitesimal length segment d of the bar at distance from the pivot point O. (see Figure 16.57). Let the mass of the infinitesimal segment be dm.
    Now the mass of the segment may be written as

    dm = (mass per unit length of the bar) (length of the segment)
    = M3d(Note:massunitlength=totalmasstotallength).

    We also note that the distance of the segment from point O, r/O=. Substituting the values found above for r/O and dm in the formula we get

    IzzO = 2()2r/O2M3ddm
    = M32()2𝑑=M3[33]2
    = M3[833(33)]=M2
    = 2kg(0.5m)2=0.5kgm2.

    Answer: IzzO=0.5kgm2

  2. 2.

    The parallel axis theorem states that

    Filename:twodisks-ang-mom
    Figure 16.58:
    IzzO=Izzcm+MrO/cm2.

    Since the rod is uniform, its center-of-mass is at its geometric center, i.e., at distance 32 from either end. From the Fig 16.58 we can see that

    rO/cm=AGAO=32=2
    Therefore,IzzO = 112M(3)2Izzcm+M(2)2
    = 912M2+M24=M2
    = 0.5kgm2(same as in (a), of course)

    Answer: IzzO=0.5kgm2

SAMPLE 16.11

Filename:pfigure-blue-90-1
Figure 16.59:

Moment of inertia of a wheel with a cut-out: A uniform rigid wheel of radius r=1ft is made eccentric by cutting out a portion of the wheel. The center-of-mass of the eccentric wheel is at C, a distance e=r3 from the geometric center O. The mass of the wheel (after deducting the cut-out) is 3.2 lbm. The moment of inertia of the wheel about point O, IzzO, is 1.8 lbmft2. We are interested in the moment of inertia Izz of the wheel about points A and B on the perimeter.

  1. 1.

    Without any calculations, guess which point, A or B, gives a higher moment of inertia. Why?

  2. 2.

    Calculate IzzC, IzzA and IzzB and compare with the guess in (a).


Solution

  1. 1.

    The moment of inertia IzzB should be higher. Moment of inertia Izz measures the geometric distribution of mass about the z-axis. But the distance of the mass from the axis counts more than the mass itself (IzzO=mr/O2𝑑m). The distance r/O of the mass appears as a quadratic term in IzzO. The total mass is the same whether we take the moment of inertia about point A or about point B. However, the distribution of mass is not the same about the two points. Due to the cut-out being closer to point B there are more “dm’s” at greater distances from point B than from point A. So, we guess that

    Answer: IzzB>IzzA

  2. 2.

    If we know the moment of inertia IzzC (about the center-of-mass) of the wheel, we can use the parallel axis theorem to find IzzA and IzzB. In the problem, we are given IzzO. But,

    Filename:pfigure-blue-110-1
    Figure 16.60:
    IzzO = IzzC+MrO/C2(parallel axis theorem)
     IzzC = IzzOMrO/C2
    = 1.8lbmft23.2lbm(1ft3)2rO/C=e=r3
    = 1.44lbmft2
    Now,IzzA = IzzC+MrA/C2=IzzC+M(2r3)2
    = 1.44lbmft2+3.2lbm(2ft3)2
    = 2.86lbmft2
    and IzzB = IzzC+MrB/C2=IzzC+M(r+r3)2
    = 1.44lbmft2+3.2lbm(1ft+1ft3)2
    = 7.13lbmft2

    Answer: IzzC=1.44lbmft2,IzzA=2.86lbmft2,IzzB=7.13lbmft2

    Clearly, IzzB>IzzA, as guessed in (a).

SAMPLE 16.12

Filename:pfigure-blue-107-1
Figure 16.61:

Moment of inertia: modeling a sphere as a point mass: A uniform solid sphere of mass m and radius r is attached to a massless rigid rod of length . The sphere swings in the xy plane. Find the error in calculating IzzO as a function of r/ if the sphere is treated as a point mass concentrated at the center-of-mass of the sphere.


Solution

The exact moment of inertia of the sphere about point O can be calculated using parallel axis theorem:

IzzO = Izzcm+m2
= 25mr2+m2.(See Table IV on inside cover)

If we treat the sphere as a point mass, the moment of inertia IzzO is

I~zzO=m2.

Therefore, the relative error in IzzO is

error = IzzOI~zzOIzzO
= 25mr2+m2m225mr2+m2
= 25r2225r22+1

From the above expression we see that for r the error is very small. From the graph of error in Fig. 16.62 we see that even for r=/5, the error in IzzO due to approximating the sphere as a point mass is less than 2%.

Filename:pfigure-s94h14p4
Figure 16.62: Relative error in IzzO of the sphere as a function of r/.

Problems for 16.3 Polar moment of inertia

Preparatory Problems

16.3.1  Find an expression for the polar moment of inertia for each of the following systems about the specified points. Answer in terms of some or all of the variables given. If a position vector 𝒓 is given then you can assume the coordinates (x,y) are given so that 𝒓=xıˆ+yȷˆ.

  1. (a)

    A particle with mass m at 𝒓 about

    1. (a)

      Its center of mass

    2. (b)

      The origin

    Answer: Izzcm=0, IzzO=mr2.

  2. (b)

    Two particles with masses m1 and m2 a distance d apart located at 𝒓1 and 𝒓2, respectively, about

    1. i)

      Mass 1

    2. ii)

      Mass 2

    3. i)

      The center of mass of the system

    4. iv)

      The origin

    Answer: Izzm1=m2|r212|, Izzm2=m1|r122|, Izzcm=(m1m2m1+m2)d2, and Izzcm=m1|𝒓1|2+m2|𝒓2|2

  3. (c)

    A collection of n particles mi at locations 𝒓i/0 about

    1. i)

      The origin

    2. ii)

      The center of mass G of the system

    Answer: IzzO=i=1nmi|ri|2 and
    Izzcm=i=1nmi|𝒓i𝒓cm|2

  4. (d)

    A uniform line segment with length and mass m with center of mass at location 𝒓G/0 about

    1. i)

      Its center of mass G

    2. ii)

      One end

    3. iii)

      The origin

    Answer: IzzG=ml212, IzzA=ml23, and
    IzzO=ml212+m|𝒓G/O|2.

  5. (e)

    A uniform hoop with mass m and radius R and center at 𝒓G/0 about

    1. i)

      Its center of mass G

    2. ii)

      A point on the hoop

    3. iiii)

      The origin

    Answer: IzzG=mR2, IzzP=2mR2, and
    IzzO=IzzG+m|𝒓G/O|2.

  6. (f)

    A uniform disk with mass m and radius R and center at 𝒓G/0 about

    1. i)

      Its center of mass G

    2. ii)

      A point on the perimeter

    3. iii)

      The origin

    Answer: IzzG=mR22, IzzP=32mR2, and
    IzzO=IzzG+m|𝒓G/O|2.

  7. (g)

    A uniform rectangle with mass m and sides and w and with center at 𝒓G/0 about

    1. i)

      Its center of mass G

    2. ii)

      A corner

    3. iii)

      The origin

    Answer: IzzG=m12(l2+w2), IzzP=m3(l2+w2), and
    IzzO=IzzG+m|𝒓G/O|2.

16.3.2  Two objects have masses m1 and m2 and polar moments of inertia I1 and I2 about their respective centers of mass G1 and G2 which are a distance d apart. What is the moment of inertia of the system about its center of mass? Answer: IzzG=I1+I2+(m1m2m1+m2)d2.

16.3.3  A point mass m=0.5kg is located at x=0.3m and y=0.4m in the xy-plane. Find the moment of inertia of the mass about the z-axis. Answer: Izz=0.125kgm2.

16.3.4  A small ball of mass 0.2kg is attached to a 1m massless rod.

  1. (a)

    What is the value of Izz of the ball about the other end of the rod? Answer: 0.2kgm2.

  2. (b)

    How much must you shorten the rod to reduce the moment of inertia of the ball by half? Answer: 0.29m.

More-Involved Problems

16.3.5  Two identical point masses are attached to the two ends of a rigid massless bar of length (one mass at each end). Locate a point along the length of the bar about which the polar moment of inertia of the system is 20% more than that calculated about the mid point of the bar. Answer: At 0.72 from either end

16.3.6  A dumbbell consists of a rigid massless bar of length and two identical point masses m and m, one at each end of the bar.

  1. (a)

    About which point on the dumbbell is its polar moment of inertia Izz a minimum and what is this minimum value? Answer: (Izz)min=m2/2, about the midpoint.

  2. (b)

    About which point on the dumbbell is its polar moment of inertia Izz a maximum and what is this maximum value? Answer: (Izz)max=m2, about either end

16.3.7  A light rigid rod AB of length 3 has a point mass m at end A and a point mass 2m at end B. Point C is the center of mass of the system. First, answer the following questions without any calculations and then do calculations to verify your guesses.

  1. (a)

    About which point A, B, or C, is the polar moment of inertia Izz of the system a minimum? Answer: C

  2. (b)

    About which point is Izz a maximum? Answer: A

  3. (c)

    What is the ratio of IzzA and IzzB? Answer: IzzA/IzzB=2

  4. (d)

    Is the radius of gyration of the system greater, smaller, or equal to the length of the rod? Answer: smaller, rgyr=IzzC/(3m)=2

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

16.3.8  Three identical particles of mass m are connected to three identical massless rods of length and welded together at point O as shown in the figure.

  1. (a)

    Guess (no calculations) which of the three moment of inertia terms IxxO, IyyO, IzzO is the smallest and which is the biggest. Answer: Biggest: IzzO; smallest: IyyO=IxxO.

  2. (b)

    Calculate the three moments of inertia to check your guess. Answer: IxxO=32m2=IyyO. IzzO=3m2.

  3. (c)

    If the orientation of the system is changed, so that one mass is along the x-axis, will your answer to part (a) change? Answer: Answer does not change.

  4. (d)

    Find the radius of gyration of the system for the polar moment of inertia. Answer: rgyr=.

Filename:tfigure8-alt-app2c
Figure 16.64:

16.3.9  Show that the polar moment of inertia IzzO of the uniform bar of length and mass m, shown in the figure, is 13m2, in two different ways:

  1. (a)

    by using the basic definition of polar moment of inertia IzzO=r2𝑑m, and

  2. (b)

    by computing Izzcm first and then using the parallel axis theorem.

Filename:sfig4-6-4a
Figure 16.65:

16.3.10  Approximately locate the center of mass of the tapered rod shown in the figure and compute the polar moment of inertia Izzcm. [Hint: use the variable thickness of the rod to define an approximate equivalent variable mass density per unit length.] Answer: xcm=40mm and Izzcm=733.3kgmm2.

Filename:sfig4-6-3
Figure 16.66:

16.3.11  A short rod of mass m and length h hangs from an inextensible and massless rod of length .

  1. (a)

    Find the moment of inertia IzzO of the rod. Answer: IzzO=mh212+m(+h2)2.

  2. (b)

    Find the moment of inertia of the rod IzzO by considering it as a point mass located at its center of mass. Answer: IzzO=m(+h2)2.

  3. (c)

    Find the percent error in IzzO in treating the bar as a point mass by comparing the expressions in parts (a) and (b). Plot the percent error versus h/. For what values of h/ is the percentage error less than 5%? Answer: Percent error in IzzO is 100(1+12(h+0.5)2)1. The percentage error less than 5% if h/<1.32.

Filename:sfig4-6-3a
Figure 16.67:

16.3.12  Parallel axis theorem? A massless square plate ABCD has four identical point masses located at its corners.

  1. (a)

    Find the polar moment of inertia Izzcm. Answer: Izzcm=2m2

  2. (b)

    Find a point P on the plate about which the system’s moment of inertia Izz is maximum. Answer: PA,B,C, or D, often, but not always, diagonally opposite the corner with the largest mass

  3. (c)

    Find the radius of gyration of the system about the rectangle center. Answer: rgyr=/2

Filename:pfigure4-4-rp12
Figure 16.68:

16.3.13  For the massless square plate with four point masses on the corners, the polar moment of inertia Izzcm=0.6kgm2. Find Ixxcm of the system. Answer: IxxO=IyyO=0.3kgm2.

Filename:pfigure-blue-38-2
Figure 16.69:

16.3.14  A uniform square plate (2m on edge) has a corner cut out. The total mass of the remaining plate is 3kg.

  1. (a)

    Where is the center of mass of the plate at the instant shown? Answer: xcm=1.67m,ycm=0.83m.

  2. (b)

    What is the moment of inertia of the plate about point O? Answer: IzzO=8kgm2.

  3. (c)

    What is the moment of inertia of the plate about its center of mass? Answer: Izzcm=1.83kgm2.

Filename:pfigure-blue-49-2
Figure 16.70:

16.3.15  A uniform thin triangular plate of mass m, height h , and base b lies in the xy-plane.

  1. (a)

    Set up the integral to find the polar moment of inertia IzzO of the plate. Answer: IzzO=2mbh0b0hx/b(x2+y2)𝑑y𝑑x.

  2. (b)

    Show that IzzO=m6(h2+3b2) by evaluating the integral in part (a).

  3. (c)

    Locate the center of mass of the plate and calculate Izzcm. Answer: xcm=2b3,ycm=h3.

Filename:pfigure4-4-rp13
Figure 16.71:

16.3.16  A uniform thin plate of mass m is cast in the shape of a semi-circular disk of radius R as shown in the figure.

  1. (a)

    Find the location of the center of mass of the plate Answer: xcm=0,ycm=4R3π.

  2. (b)

    Find the polar moment of inertia of the plate, Izzcm. [Hint: It may be easier to set up and evaluate the integral for IzzO and then use the parallel axis theorem to calculate Izzcm.] Answer: Izzcm=(9π231)mR218π2.

Filename:f92h7p1
Figure 16.72:

16.3.17  A uniform square plate of side =250mm has a circular cut-out of radius r=50mm. The mass of the plate is m=12kg.

  1. (a)

    Find the polar moment of inertia of the plate about its center. Answer: Izzcm=m26(13π(r/)4)1π(r/)2).

  2. (b)

    Plot Izzcm versus r/.

  3. (c)

    Find the limiting values of Izzcm for r=0 and r=. Answer: limr0Izzcm=m26 and
    limrIzzcm=m26(3π1π1).

Filename:pfigure-blue-52-2
Figure 16.73:

16.3.18  A uniform thin circular disk of radius r=100mm and mass m=2kg has a rectangular slot of width w=10mm cut into it as shown in the figure. You can approximate the right edge of the cutout as a straight line (thus your calculation would be exact for a system that has a very thin sliver of material closing the gap at the right).

  1. (a)

    Find the polar moment of inertia IzzO of the disk. Answer: IzzO=0.01010904kgm2.

  2. (b)

    Locate the center of mass of the disk and calculate Izzcm. Answer: xcm=1.64mm, ycm=0mm, and Izzcm=0.01010364kgm2.

  3. (c)

    estimate the error in your calculation of the mass, center of mass location and moment of inertia due to not subtracting out the sliver of material at the right edge of the slot. Answer: Not subtracting out the sliver of material at the right edge of the slot leads to an error of 0.053%.

Filename:pfigure4-4-rp16
Figure 16.74:

16.3.19  Calculate the polar moment of inertia of a uniform square plate with mass m and side various ways. Choose numerical values for m and if you want to address this as a numerical calculation rather than a theoretical one.

  1. (a)

    As a single uniform plate

  2. (b)

    As a composite of 4 plates, each with sides /2

  3. (c)

    As a composite of 100 plates, each with sides /10 (These first three methods should give exactly the same answer).

  4. (d)

    For the composite of 100 plates now neglect the terms which concern the moments of inertia of each small square about its center of mass. What is the error in making this neglect? Answer: Exactly 1% .

16.4 Dynamics of a rigid object in planar circular motion

We now want to find what we can about the forces on and motions of a single planar object that is hinged at one point. The most famous examples are pendula and gears. But many things move like pendula, like people or towers falling over or boats rocking in the water. And many machine parts besides gears rotate about a bearing including propellers, lawn mower blades, and printing-press rollers. None of these things are strictly planar objects, but, for many analyses, the planar approximation is more or less accurate

margin: Accuracy of planar analysis. The planar analysis might be accurate for one of two different reasons: 1) the objects are nearly planar, or 2) the projection of the 3D equations to 2D gives the planar equations.

.

The basic method, as always, is to use a free-body diagram and kinematics results to evaluate terms in the linear momentum, angular momentum and energy balance equations. Then these equations are used to calculate unknown forces or to find differential equations of motion. The terms involving force and moment are found from a free-body diagram exactly as in statics. The terms involving motion, namely momentum, angular momentum and energy, are evaluated using the various tools developed earlier in this chapter.

Mechanics and the motion quantities

We can evaluate all the momentum and energy terms in the equations of motion (inside cover), namely: 𝑳,𝑳˙,𝑯/C,𝑯˙/C,EK and E˙K (for any reference point C of our choosing) if we can calculate the velocity and acceleration of every point in the system. For circular motion of a rigid object about point C, we know from sec. 16.2 that the velocities and accelerations of a point at 𝒓=𝒓/C are

𝒗 = 𝝎×𝒓,
𝒂 = 𝝎˙×𝒓+𝝎×(𝝎×𝒓),
= 𝝎˙×𝒓ω2𝒓

where 𝝎 is the angular velocity of the object relative to a fixed frame. This situation is a little more complex than the special case of straight-line motion in chapter 14, where all points in a system had the same acceleration as each other, but is still quite manageable.

The most general case of 2D rigid-object circular motion is an arbitrarily shaped 2D rigid object with arbitrary ω and ω˙. The situation shown in fig. 16.75 shows the general case (but in that example exactly two forces are applied, as opposed to an arbitrary number).

Filename:pfigure4-4-rp14
Figure 16.75: a) A rigid object rotates about 0. It has center of mass at G. Forces 𝑭1 and 𝑭2 are applied. b) The free-body diagram shows all the forces acting on the object, including the reaction force (𝑭0=F0xıˆ+F0yȷˆ) at 0 and the gravity force mg at G.

Linear momentum: 𝑳 and 𝑳˙

For any system in any motion we have the general result that

𝑳=mtot𝒗cm and 𝑳˙=mtot𝒂cm.

For a rigid object, the center-of-mass is a particular point G that is fixed relative to the object. So the velocity and acceleration of that point can be expressed the same way as for any other point. So, for an object in planar rotational motion about 0

𝑳=mtot𝝎×𝒓G/0

and

𝑳˙=mtot(𝝎˙×𝒓G/0ω2𝒓G/0)𝒂G.

If the center-of-mass is at 0 the momentum and its rate of change are both zero. But if the center-of-mass is off the axis of rotation, there must be a net force on the object with a component parallel to 𝒓0/G (if ω0) and also a component orthogonal to 𝒓0/G (if ω˙0). This net force might be applied at 0 or G or any other place(s) on the object.

Angular momentum: 𝑯/O and 𝑯˙/O

The angular momentum itself is easy enough to calculate, using the short hand notation that 𝒓 is the position vector 𝒓/0 of a point relative to hinge point O.

𝑯/O=allmass𝒓×𝒗𝑑m(a)=𝒓×(𝝎×𝒓)𝑑m(b)=ω𝒌ˆr2𝑑m(c)=IzzOIzzO is the ‘polar’ moment of inertia.ω𝒌ˆ  H0=ωIzzO(d) (16.50)

Here eqn. (16.50)c is a vector equation. But since both sides are in the 𝒌ˆ direction we can dot both sides with 𝒌ˆ to get the scalar moment equation eqn. (16.50)d, taking both Mnet and ω as positive when counterclockwise.

The moment of inertia shortcut. Note that an integral over all the mass of 𝒓×𝒗 has been reduced to a simple multiplication by using that

r2𝑑m=Izz0.

This shortcut is especially useful because r2𝑑m does not change as the object rotates about 0.

However, the moment of inertia is never needed for problem solving. One can just evaluate, for example, 𝑯/O by evaluating the integral 𝒓×𝒗𝑑m for the specific problem at hand.

To get the all important term for angular momentum balance, namely 𝑯˙/0, for this system we can differentiate eqn. (16.50). We can also use the general expression for 𝑯˙/O to write the angular momentum balance equation as follows.

𝑯˙/O=rate of change of angular momentum/0(a)=allmass𝒓×𝒂𝑑m(b)=𝒓×(ω2𝒓+ω˙𝒌ˆ×𝒓𝒂)𝑑m(c)=𝒓×(ω˙𝒌ˆ×𝒓)𝑑m(because 𝒓×𝒓=𝟎)(d)=𝒓×(ω˙𝒌ˆ×𝒓)𝑑m(e)𝑯˙/O=ω˙𝒌ˆr2𝑑m=ω˙Izz0𝒌ˆ(f)  H˙/0=ω˙r2𝑑m=ω˙Izz0(g) (16.51)

We get from eqn. (16.51)e to eqn. (16.51)f by noting that 𝒓 is perpendicular to 𝒌ˆ. Thus, using the right hand rule twice we get 𝒓×(𝒌ˆ×𝒓)=r2𝒌ˆ.

Eqn. 16.51f and eqn. (16.51)g are the vector and scalar versions for the rate of change of angular momentum with respect to point 0 for rotation of a planar object about 0. Repeating,

𝑯˙/O = ω˙𝒌ˆr2𝑑m=ω˙Izz0𝒌ˆ
andH˙/0 = ω˙r2𝑑m=ω˙Izz0. (16.52)

Again, 𝑯˙/O can always be evaluated either of two ways: 1) as a sum or integral of 𝒓/0×𝒂 over all the mass or 2) using the shortcut with Izz0.

Power and Energy

Let’s assume that there are a set of point forces applied to an object

margin: Forces distributed on a curve, surface or throughout a volume could be treated similarly with sums replaced by line, surface, or volume integrals.

.

And, to be contrary, let’s assume the mass is continuously distributed (the derivation for rigidly connected point masses would be similar). The power balance equation for one rotating rigid object is (discussed below):

Net power in=rate of change of kinetic energy(a)P=E˙K(b)allappliedforces𝑭i𝒗i=ddtallmass12v2𝑑m(c)𝑭i(𝝎×𝒓i)=ddt12(𝝎×𝒓)(𝝎×𝒓)𝑑m(d)𝝎(𝒓i×𝑭i)=ddt12ω2r2𝑑m(e)𝝎(𝒓i×𝑭i)=ddt(12ω2)r2𝑑m(f)𝝎𝑴i=ω˙ωr2𝑑m(g)𝝎𝑴tot=𝝎˙(𝝎r2𝑑m)𝑯/0(h) (16.53)

When not directly labeled, positions and moments are assumed to be relative to the hinge at 0. Derivation 16.53 is two derivations in one. The left side about power and the right side about kinetic energy. Let’s discuss one at a time.

Power. On the left side of eqn. (16.53) we note in (c) that the power of each force is the dot product of the force with the velocity of the point it touches. In (d) we use what we know about the velocities of points on rotating rigid bodies. In (e) we use the vector identity 𝑨𝑩×𝑪=𝑩𝑪×𝑨 (see sec. 1.3, 1.3). In (f) we note that 𝝎 is common to all points so factors out of the sum. In (g) we note that 𝒓×𝑭i is the moment of the force about pt O. And in (g) we sum the moments of the forces. So the power of a set of forces acting on a rigid object is the product of their net moment (about 0) and the object angular velocity,

P=𝝎𝑴tot. (16.54)

Kinetic energy. On the right side of eqn. (16.53) we note in (c) that the kinetic energy is the sum of the kinetic energy of the mass increments. In (d) we use what we know about the velocities of these bits of mass, given that they are on a common rotating object. In (e) we use that the magnitude of the cross product of orthogonal vectors is the product of the magnitudes (|𝑨×𝑩|=AB) and that the dot product of a vector with itself is its magnitude squared (𝑨𝑨=A2). In (f) we factor out ω2 because it is common to all the mass increments and note that the remaining integral is constant in time for a rigid object. In (g) we carry out the derivative. In (h) we de-simplify the result from (g) in order to show a more general form that we will find later in 3D mechanics. Eqn. (h) follows from (g) because 𝝎 is parallel to 𝝎˙ for 2D rotations.

Note that we started here with the basic power balance equation from the front inside cover. Instead, we could have derived power balance from our angular momentum balance expression (see box 16.4 on 16.4).

Box 16.8 The relation between angular momentum balance and power balance

For this system, angular momentum balance can be derived from power balance and vice versa. Thus neither is essentially more fundamental than the other and both are reliable. First we can derive power balance from angular momentum balance as follows:

𝑴net=ω˙𝒌ˆr2𝑑m𝝎𝑴net=𝝎(ω˙𝒌ˆr2𝑑m).P=E˙K (16.55)

That is, when we dot both sides of the angular momentum equation with 𝝎 we get on the left side a term which we recognize as the power of the forces and on the right side a term which is the rate of change of kinetic energy.

The opposite derivation starts with the power balance fig. 16.53(g)

𝝎𝑴i=ω˙ωr2𝑑m(g)  ω(𝒌ˆ𝑴i)=ω˙ωr2𝑑m  (𝒌ˆ𝑴i)=ω˙r2𝑑m (16.56)

and, assuming ω0, divide by ω to get the angular momentum equation for planar rotational motion.

Example: Spinning disk

The round flat uniform disk in fig. 16.76 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:

Filename:pfigure4-4-rp15
Figure 16.76: A uniform disk turned by a motor at a constant rate. The free-body diagram shows a force and moment equivalent to the force system that acts including gravity, bearing forces, etc. A bit of mass dm occupies the space between the radial lines at θ and θ+dθ and between the circles at radii R and R+dR.
𝑴i/O = 𝑯˙/O
 𝑴 = 𝒓/O×𝒂𝑑m
= 0R002π(R𝒆ˆR)×(Rω2𝒆ˆR)(mtotπRO2)ρRdθdRdAdm
= 𝟎𝑑θ𝑑R
= 𝟎.

So the net force and moment needed are 𝑭=𝟎 and 𝑴=𝟎. Like a particle that moves at constant velocity with no force, a uniform disk rotates at constant rate with no torque (at least in 2D).

Using moment-of-inertia in 2-D circular motion dynamics

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 earlier chapters. 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. As you recall from earlier chapters, the 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 wishes for something so simple for the calculation of angular momentum.

It turns out that we are in luck if we are only interested in the two-dimensional motion of two-dimensional rigid bodies. The luck is not so great for 3-D rigid bodies but still there is some simplification. For general motion of non-rigid bodies there is no simplification to be had.

The simplification is to use the moment of inertia for the bodies rather than evaluating the momenta and energy quantities as integrals and sums. Of course one may have to do a sum or integral to evaluate IIzzcm or [𝑰cm] but once this calculation is done, one need not work with the integrals while worrying about the dynamics. At this point we will assume that you are comfortable calculating and looking-up moments of inertia. We proceed to use it for the purposes of studying mechanics.

For constant rate rotation, we can calculate the velocity and acceleration of various points on a rigid body using 𝒗=𝝎×𝒓 and 𝒂=𝝎×(𝝎×𝒓). So we can calculate the various motion quantities of interest: linear momentum 𝑳, rate of change of linear momentum 𝑳˙, angular momentum 𝑯, rate of change of angular momentum 𝑯˙, and kinetic energy EK.

Consider a two-dimensional rigid body like that shown in fig. 16.77.

Filename:tfigure4-spherical-rotaxis
Figure 16.77: A two-dimensional body is rotating around the point O at constant rate ω. A differential bit of mass dm is shown. The center-of-mass is also shown.

Now let us consider the various motion quantities in turn. First the linear momentum 𝑳. The linear momentum of any system in any motion is 𝑳=𝒗cmmtot. So, for a rigid body spinning at constant rate ω about point O (using 𝝎=ω𝒌ˆ):

𝑳=𝒗cmmtot=𝝎×𝒓cm/omtot.

Similarly, for any system, we can calculate the rate of change of linear momentum 𝑳˙ as 𝑳˙=𝒂cmmtot. So, for a rigid body spinning at constant rate,

𝑳˙=𝒂cmmtot=𝝎×(𝝎×𝒓cm/o)mtot.

That is, the linear momentum is correctly calculated for this special motion, as it is for all motions, by thinking of the body as a point mass at the center-of-mass.

Unlike the calculation of linear momentum, the angular momentum turns out to be something different than would be calculated by using a point mass at the center of mass. You can remember this important fact by looking at the case when the rotation is about the center-of-mass (point O coincides with the center-of-mass). In this case one can intuitively see that the angular momentum of a rigid body is not zero even though the center-of-mass is not moving. Here’s the calculation just to be sure:

𝑯/O=𝒓/O×𝒗𝑑m(by definition of 𝑯/O)=𝒓/O×(𝝎×𝒓/O)𝑑m(using 𝒗=𝝎×𝒓)=(x/Oıˆ+y/Oȷˆ)×[(ω𝒌ˆ)×(x/Oıˆ+y/Oȷˆ)]𝑑m(substituting 𝒓/O and 𝝎)={(x/O2+y/O2)𝑑m}ω𝒌ˆ(doing cross products)={r/O2𝑑m}ω𝒌ˆ=IzzOIzzO is the ‘polar’ moment of inertia.ω𝒌ˆ

We have defined the ‘polar’ moment of inertia as Izzo=r/o2𝑑m. In order to calculate Izzo for a specific body, assuming uniform mass distribution for example, one must convert the differential quantity of mass dm into a differential of geometric quantities. For a line or curve, dm=ρd; for a plate or surface, dm=ρdA, and for a 3-D region, dm=ρdV. d, dA, and dV are differential line, area, and volume elements, respectively. In each case, ρ is the mass density per unit length, per unit area, or per unit volume, respectively. To avoid clutter, we do not define a different symbol for the density in each geometric case. The differential elements must be further defined depending on the coordinate systems chosen for the calculation; e.g., for rectangular coordinates, dA=dxdy or, for polar coordinates, dA=rdrdθ.

Since 𝑯 and 𝝎 always point in the 𝒌ˆ direction for two dimensional problems people often just think of angular momentum as a scalar and write the equation above simply as ‘H=Iω,’ the form usually seen in elementary physics courses.

The derivation above has a feature that one might not notice at first sight. The quantity called IzzO does not depend on the rotation of the body. That is, the value of the integral does not change with time, so IzzO is a constant. So, perhaps unsurprisingly, a two-dimensional body spinning about the z-axis through O has constant angular momentum about O if it spins at a constant rate. margin: Note that the angular momentum about some other point than O will not be constant unless the center-of-mass does not accelerate (i.e., is at point O).

𝑯˙/O=𝟎.

Now, of course we could find this result about constant rate motion of 2-D bodies somewhat more cumbersomely by plugging in the general formula for rate of change of angular momentum as follows:

𝑯˙/O=𝒓/O×𝒂𝑑m=𝒓/O×(𝝎×(𝝎×𝒓/O))𝑑m=(x/Oıˆ+y/Oȷˆ)×[ω𝒌ˆ×(ω𝒌ˆ×(x/Oıˆ+y/Oȷˆ))]𝑑m=𝟎. (16.57)

Using moment of inertia about the center of mass. Often it is easier to think of the motion as composed of two parts, motion of the center of mass and motion relative to the center of mass, as explained in box 16.4 on page 16.4. Thus we have two terms for angular momentum 𝑯/O and its rate of change 𝑯˙/O:

𝑯/O =𝒓G/O×m𝒗cm+Izzcmω𝒌ˆ (16.58)
𝑯˙/O =𝒓G/O×m𝒂cm+Izzcmα𝒌ˆ (16.59)

Kinetic energy. Finally, we can calculate the kinetic energy by adding up 12mivi2 for all the bits of mass on a 2-D body spinning about the z-axis:

EK=12v2𝑑m=12(ωr)2𝑑m=12ω2r2𝑑m=12Izzoω2. (16.60)

The kinetic energy can also be written as a sum of contributions of motion of the center of mass and motion relative to the center of mass,

EK=12mvcm2+12Izzcmω2. (16.61)

Example: Pendulum disk

Filename:tfigure5-7
Figure 16.78:

For the disk shown in fig. 16.78, we can calculate the rate of change of angular momentum about point O as

𝑯˙/O = 𝒓G/O×m𝒂cm+Izzcmα𝒌ˆ
= R2mθ¨𝒌ˆ+Izzcmθ¨𝒌ˆ
= (Izzcm+R2m)θ¨𝒌ˆ.

Alternatively, we could calculate directly

𝑯˙/O = IzzOα𝒌ˆ
= (Izzcm+R2m)by the parallel axis theoremθ¨𝒌ˆ.

Note that we are using the planarity of the objects and of their motion for our calculations

margin: 2D vs 3D. Beware of falling into the common misconception that the formula M=Iα applies in three dimensions by just thinking of the scalars as vectors and matrices. In 3D the formula 𝑯˙/O=[𝑰O]𝝎˙𝜶 is only correct when 𝝎 is zero or when 𝝎 is an eigen vector of [I/O]. That is, the vector equation Moments about O=[𝑰O]𝜶 is generally wrong in 3D.

.

The equation for linear momentum balance is the same as always, we just need to calculate the acceleration of the center-of-mass of the spinning body.

𝑳˙=mtot𝒂cm=mtot[𝝎×(𝝎×𝒓cm/O)+𝝎˙×𝒓cm/O] (16.62)

Finally, the kinetic energy for a planar rigid body rotating in the plane is:

EK=12𝝎([𝑰cm]𝝎)+12mvcm2𝒗cm=𝝎×𝒓cm/O.

Box 16.9 Simplifying 𝑯/C using the center of mass

The definition of angular momentum relative to a point C is

𝑯/C=𝒓i/C×mi𝒗i.

If we rewrite 𝒗i as

𝒗i=(𝒗i𝒗cm)+𝒗cm=𝒗i/cm+𝒗cm

and

𝒓i=(𝒓i𝒓cm)+𝒓cm=𝒓i/cm+𝒓cm

then

𝑯/C = (𝒓cm+𝒓i/cm)×[𝒗cm+𝒗i/cm]mi.
= 𝒓cm×𝒗cmmi+𝒓i/cm×𝒗i/cmmi
+𝒓cm×𝒗i/cmmi+𝒓i/cm×𝒗cmmi
= 𝒓cm×𝒗cmmtot+𝒓i/cm×𝒗i/cmmi
+𝒓cm×[𝒗i/cmmi]𝟎+[𝒓i/cmmi]𝟎×𝒗cm

So,

𝑯/C=𝒓cm×𝒗cmmtotcontribution of center of mass motion+𝒓i/cm×𝒗i/cmmicontribution of motion relative to center-of-mass.

The reason 𝒓i/cmmi=𝟎 is somewhat intuitive. It is what you would calculate if you were looking for the center-of-mass relative to the center of mass. More formally,

𝒓i/cmmi = (𝒓i𝒓cm)mi
= 𝒓imimtot𝒓cmmtot𝒓cm
= 𝟎.

Similarly, 𝒗i/cmmi=𝟎 because it is what you would calculate if you were looking for the velocity of the center-of-mass relative to the center of mass.

The central result of this box is that

angular momentum of any system is that due to motion of the center-of-mass plus motion relative to the center-of-mass.

SAMPLE 16.13

Filename:tfigure5-gen-rigid-body
Figure 16.79: A rod goes in circles at a constant rate.

A rod in a constant rate circular motion: A uniform rod of mass m and length is connected to a motor at end O. A ball of mass m is attached to the rod at end B. The motor turns the rod in the counterclockwise direction at a constant angular speed ω. There is gravity pointing in the ȷˆ direction. Find the torque applied by the motor (i) at the instant shown and (ii) when θ=0, 90, 180. How does the torque change if the angular speed is doubled?


Solution

The FBD of the rod and ball system is shown in Fig. 16.80(a). Since the system is undergoing circular motion at a constant speed, the acceleration of the ball as well as every point on the rod is just radial (pointing towards the center of rotation O) and is given by 𝒂=ω2r𝝀ˆ where r is the radial distance from the center O to the point of interest and 𝝀ˆ is a unit vector along OB pointing away from O (Fig. 16.80(b)).

Filename:tfigure5-term1-a
Figure 16.80: A rod goes in circles at a constant rate.

Angular Momentum Balance about point O gives

𝑴O=𝑯˙/O
𝑴O = 𝒓G/O×(mgȷˆ)+𝒓B/O×(mgȷˆ)+M𝒌ˆ (16.63)
= 2cosθmg𝒌ˆcosθmg𝒌ˆ+M𝒌ˆ
= (M32mgcosθ)𝒌ˆ
𝑯˙/O = 𝒓B/O×m𝒂B(𝑯˙/O)ball+m𝒓dm/O×𝒂dm𝑑m(𝑯˙/O)rod (16.64)
= 𝝀ˆ×(mω2𝝀ˆ)+ms𝝀ˆ𝒓dm/O×(ω2s𝝀ˆ)𝒂dm𝑑m
= 𝟎(since 𝝀ˆ×𝝀ˆ=𝟎)

(i) Equating (16.63) and (16.64) we get

M=32mgcosθ.

Answer: M=32mgcosθ

(ii) Substituting the given values of θ in the above expression we get

M(θ=0)=32mg,M(θ=90)=0M(θ=180)=32mg

Answer: M(0)=32mg,M(90)=0M(180)=32mg

margin: The values obtained above make sense (at least qualitatively). To make the rod and the ball go up from the 0 position, the motor has to apply some torque in the counterclockwise direction. In the 90 position no torque is required for the dynamic balance. In the 180 position the system is accelerating downwards under gravity; therefore, the motor has to apply a clockwise torque to make the system maintain a uniform speed.

It is clear from the expression of the torque that it does not depend on the value of the angular speed ω! Therefore, the torque will not change if the speed is doubled. In fact, as long as the speed remains constant at any value, the only torque required to maintain the motion is the torque to counteract the moments due to gravity at O.

SAMPLE 16.14

Filename:tfigure5-term1-b
Figure 16.81: A rectangular plate is released from rest from the position shown.

At the onset of circular motion: A 2×4 rectangular plate of mass 20lbm is pivoted at one of its corners as shown in the figure. The plate is released from rest in the position shown. Find the force on the support immediately after release.

Solution The free-body diagram of the plate is shown in Fig. 16.82. The force 𝑭 applied on the plate by the support is unknown.

Filename:sfig1-2-12
Figure 16.82: (a) The free-body diagram of the plate. (b) Computation of the integral in 𝑯˙/O=ω˙𝒌ˆmr2𝑑m. (c) The geometry of motion. From the given dimensions,   𝒆ˆR=aıˆbȷˆ(a2+b2), 𝒆ˆθ=bıˆ+aȷˆ(a2+b2), and rG/O=a2+b22.

The linear momentum balance for the plate gives

𝑭 = m𝒂G
𝑭mgȷˆ = m(θ¨rG/O𝒆ˆθθ˙2rG/O𝒆ˆR) (16.65)
= mθ¨rG/O𝒆ˆθ(since θ˙=0 at t=0).

Thus to find 𝑭 we need to find θ¨.

The angular momentum balance for the plate about the fixed support point O gives

𝑴O = 𝑯˙O
where
𝑴O = 𝒓G/O×mg(ȷˆ)
= (a2ıˆb2ȷˆ𝒓G/O)×mg(ȷˆ)=mga2𝒌ˆ,
and
𝑯˙/O = ω˙𝒌ˆmr2𝑑m=θ¨𝒌ˆ0b0a(x2+y2)r2mabdxdydm
= m(a2+b2)3θ¨𝒌ˆ.
Thus,
mga2𝒌ˆ = m(a2+b2)3θ¨𝒌ˆ
 θ¨ = 3ga2(a2+b2)
= 332.2ft/s24ft2(16+4)ft2=9.66rad/s2.

From eqn. (16.65), the support force is now readily calculated:

𝑭 = mgȷˆ+mθ¨rG/O𝒆ˆθ
= mgȷˆ+mθ¨a2+b22bıˆ+aȷˆ(a2+b2)
= 12mθ¨bıˆ+(mg+12mθ¨a)ȷˆ

Using the given numerical values of m,a, and b, θ¨=9.66rad/s2, and g=32.2ft/s2, we get

𝑭=(6ıˆ+8ȷˆ)lbf.

Answer: 𝑭=(6ıˆ+8ȷˆ)lbf

SAMPLE 16.15

Filename:sfig4-6-5
Figure 16.83: A compound gear train.

A compound gear train. When the gear of an input shaft, often called the driver or the pinion, is directly meshed in with the gear of an output shaft, the motion of the output shaft is opposite to that of the input shaft. To get the output motion in the same direction as that of the input motion, an idler gear is used. If the idler shaft has more than one gear in mesh, then the gear train is called a compound gear train.

In the gear train shown in Fig. 16.83, the input shaft is rotating at 2000 rpm and the input torque is 200 N-m. The efficiency (defined as the ratio of output power to input power) of the train is 0.96 and the various radii of the gears are: RA=5cm,RB=8cm,RC=4cm, and RD=10cm. Find

  1. 1.

    the input power Pin and the output power Pout,

  2. 2.

    the output speed ωout, and

  3. 3.

    the output torque.

Solution

  1. 1.

    The power:

    Pin = Minωin=200Nm2000rpm
    = 400000Nmrevmin2π1rev1min60s
    = 41887.9Nm/s42 kW.
     pout = efficiencyPin=0.9642 kW40 kW

    Answer: Pin=42 kW,pout=40 kW

  2. 2.
    Filename:sfig4-6-5a
    Figure 16.84:

    The angular speed of meshing gears can be easily calculated by realizing that the linear speed of the point of contact has to be the same irrespective of which gear’s speed and geometry is used to calculate it. Thus,

    vP = ωinRA=ωBRB
     ωB = ωinRARB
    and vR = ωCRC=ωoutRD
     ωout = ωCRCRD
    ButωC = ωB
     ωout = ωinRARBRCRD
    = 2000rpm58410=500rpm.

    Answer: ωout=500rpm

  3. 3.

    The output torque,

    Mout=Poutωout = 40 kW500rpm=405001000Nmsminrev1rev2π60s1min
    = 764Nm.

    Answer: Mout=764Nm

SAMPLE 16.16

Filename:sfig4-6-5b
Figure 16.85: An accelerating compound gear train.

An accelerating gear train. In the gear train shown in Fig. 16.85, the torque at the input shaft is Min=200Nm and the angular acceleration is αin=50rad/s2. The radii of the various gears are: RA=5cm,RB=8cm,RC=4cm, and RD=10cm and the moments of inertia about the shaft axis passing through their respective centers are: IA=0.1kgm2,IBC=5IA,ID=4IA. Find the output torque Mout of the gear train.


Solution

Since the difference between the input power and the output power is used in accelerating the gears, we may write

PinPout=EK˙

Let Mout be the output torque of the gear train. Then,

PinPout=MinωinMoutωout. (16.66)

Now,

EK˙ = ddt(EK)
= ddt(12IAωin2+12IBCωBC2+12IDωout2)
= IAωinω˙in+IBCωBCω˙BC+IDωoutω˙out
= IAωinαin+5IAωBCαBC+4IAωoutαout. (16.68)
Filename:sfig4-6-5c
Figure 16.86:

The velocity or acceleration of the point of contact between two meshing gears has to be the same irrespective of which meshing gear’s geometry and motion is used to compute them.

The different ω’s and the α’s can be related by realizing that the linear speed or the tangential acceleration of the point of contact between any two meshing gears has to be the same irrespective of which gear’s speed and geometry is used to calculate it. Thus, using the linear speed and tangential acceleration calculations for points P and R in Fig. 16.86, we can find ωB, αB and ωout, αout. Considering the linear speed of point P, we find

vP = ωinRA=ωBRB
 ωB = ωinRARB,

and considering the tangential acceleration of point P, (aP)θ, we find

(aP)θ = αinRA=αBRB
 αB = αinRARB.

Similarly,

vR = ωCRC=ωoutRD
 ωout = ωCRCRD,

and

(aR)θ = αCRC=αoutRD
 αout = αCRCRD.

But

ωC = ωB=ωBC
 ωout = ωinRARBRCRD

and

αC = αB=αBC
 αout = αinRARBRCRD.

Substituting these expressions for ωout,αout,ωBC and αBC in equations (16.66) and (16.68), we get

PinPout = MinωinMoutωinRARBRCRD
= ωin(MinMoutRARBRCRD).
EK˙ = IA[ωinαin+5ωinαin(RARB)2+4ωinαin(RARBRCRD)2]
= IAωin[αin+5αin(RARB)2+4αin(RARBRCRD)2].

Now equating the two quantities, PinPout and EK˙, and canceling ωin from both sides, we obtain

MoutRARBRCRD = MinIAαin[1+5(RARB)2+4(RARBRCRD)2]
Mout58410 = 200Nm5kgm2rad/s2[1+5(58)2+4(58410)2]
Mout = 735.94Nm
736Nm.

Answer: Mout=736Nm

SAMPLE 16.17

Filename:sfig4-6-5d
Figure 16.87: Two drums with strings wrapped around are used to pull up a mass m.

Drums used as pulleys. Two drums, A and B of radii Ro=200mm and Ri=100mm are welded together. The combined mass of the drums is mD=20kg and the combined moment of inertia about the z-axis passing through their common center O is Izz/O=1.6kgm2. A string attached to and wrapped around drum B supports a mass m=2kg. The string wrapped around drum A is pulled with a force F=20N as shown in Fig. 16.87. Assume there is no slip between the strings and the drums. Find

  1. 1.

    the angular acceleration of the drums,

  2. 2.

    the tension in the string supporting mass m, and

  3. 3.

    the acceleration of mass m.

margin:

Solution

The free-body diagram  of the drums and the mass are shown in Fig. 16.88 separately where T is the tension in the string supporting mass m and Ox and Oy are the support reactions at O. Since the drums can only rotate about the z-axis, let

𝝎=ω𝒌ˆ and 𝝎˙=ω˙𝒌ˆ.
Filename:sfig4-5-6
Figure 16.88: Free-body diagram of the drums and the mass m. T is the tension in the string supporting mass m and Ox and Oy are the reactions of the support at O.

Now, let us do angular momentum balance about the center of rotation O:

𝑴O=𝑯˙/O
𝑴O = TRi𝒌ˆFRo𝒌ˆ
= (TRiFRo)𝒌ˆ.

Since the motion is restricted to the xy-plane (i.e., 2-D motion), the rate of change of angular momentum 𝑯˙/O may be computed as

𝑯˙/O = Izz/cmω˙𝒌ˆ+𝒓cm/O×𝒂cmmD
= Izz/Oω˙𝒌ˆ+𝒓O/O0×𝒂cm0mD
= Izz/Oω˙𝒌ˆ.

Setting 𝑴O=𝑯˙/O we get

TRiFRo=Izz/Oω˙. (16.69)

Now, let us write linear momentum balance, 𝑭=m𝒂, for mass m:

(Tmg)ȷˆ𝑭=m𝒂.

Do we know anything about acceleration 𝒂 of the mass? Yes, we know its direction (±ȷˆ) and we also know that it has to be the same as the tangential acceleration (𝒂D)θ of point D on drum B (why?). Thus,

𝒂 = (𝒂D)θ (16.70)
= ω˙𝒌ˆ×(Riıˆ)
= ω˙Riȷˆ.

Therefore,

Tmg=mω˙Ri. (16.71)
  1. 1.

    Calculation of ω˙: We now have two equations, (16.69) and (16.71), and two unknowns, ω˙ and T. Subtracting Ri times Eqn.(16.71) from Eqn. (16.69) we get

    FRo+mgRi = (Izz/O+mRi2)ω˙
     ω˙ = FRo+mgRi(Izz/O+mRi2)
    = 20N0.2m+2kg9.81m/s20.1m1.6kgm2+2kg(0.1m)2
    = 2.038kgm2/s21.62kgm2
    = 1.2581s2

    Answer: ω˙=1.26rad/s2kˆ

  2. 2.

    Calculation of tension 𝐓: From equation (16.71):

    T = mgmω˙Ri
    = 2kg9.81m/s22kg(1.26s2)0.1m
    = 19.87N

    Answer: T=19.87N

  3. 3.

    Calculation of acceleration of the mass: Since the acceleration of the mass is the same as the tangential acceleration of point D on the drum, we get (from eqn. (16.70))

    𝒂 = (𝒂D)θ=ω˙Riȷˆ
    = (1.26s2)0.1m
    = 0.126m/s2ȷˆ

    Answer: a=0.13m/s2ȷˆ

Comments: It is important to understand why the acceleration of the mass is the same as the tangential acceleration of point D on the drum. We have assumed (as is common practice) that the string is massless and inextensible. Therefore each point of the string supporting the mass must have the same linear displacement, velocity, and acceleration as the mass. Now think about the point on the string which is momentarily in contact with point D of the drum. Since there is no relative slip between the drum and the string, the two points must have the same vertical acceleration. This vertical acceleration for point D on the drum is the tangential acceleration (𝒂D)θ.

SAMPLE 16.18   Energy method for pulley dynamics: Consider the pulley problem of Sample 16.86 again. Use energy method to

  1. 1.

    find the angular acceleration of the pulley, and

  2. 2.

    the acceleration of the mass.


Solution

In energy method we use speeds, not velocities. Therefore, we have to be careful in our thinking about the direction of motion. In the present problem, let us assume that the pulley rotates and accelerates clockwise. Consequently, the mass moves up against gravity.

  1. 1.

    The energy equation we want to use is

    P=E˙K.

    The power P is given by P=𝑭i𝒗i where the sum is carried out over all external forces. For the mass and pulley system the external forces that do work are F and mg.

    The other external forces on the system—the reaction force of the support point O and the weight of the pulley—are acting at point O (see fig. 16.89). But, since point O is stationary, these forces do no work.

    Filename:sfig4-5-6a
    Figure 16.89: Free-body-diagram of the pulley-mass system together. Note that the forces acting at the support point O do no work since point O is stationary.

    Therefore,

    P = 𝑭𝒗A+m𝒈𝒗m (16.72)
    = FıˆvAıˆ+(mgȷˆ)vDȷˆ𝒗m
    = FvAmgvD.

    Now we need to calculate the rate of change of kinetic energy E˙K. There are two objects here that have kinetic energy—the hanging mass and the pulley. Hence,

    E˙K=(E˙K)m+(E˙K)pulley.

    The hanging mass has pure translational motion and hence its kinetic energy is

    (EK)m=12mv2

    where v is the linear speed of the mass. If we assume the string to be inextensible, then the linear speed v of the mass has to be the same as the tangential speed of point D of the pulley. Thus v=vD and

    (EK)m=12mvD2.

    The pulley, on the other hand, has pure rotational motion about point O, and hence its kinetic energy is given by

    (EK)pulley=12IzzOω2.

    Summing the two kinetic energies and differentiating with respect to time t, we get

    E˙K = ddt(12mvD2+12IzzOω2) (16.73)
    = mvDv˙D+IzzOωω˙. (16.74)

    Now equating the power and the rate of change of kinetic energy from eqns. eqn. (16.72) and eqn. (16.74), we get

    FvAmgvD=mvDv˙D+IzzOωω˙.

    From kinematics of circular motion,

    vA = ωRo,
    vD = ωRi
    and v˙D(aD)θ = ω˙Ri.

    Substituting these values in the power balance equation above, we get

    ω(FRomgRi) = ωω˙(mRi2+IzzO)
     ω˙ = FRomgRi(IzzO+mRi2)
    = 20N0.2m2kg9.81m/s20.1m1.6kgm2+2kg(0.1m)2
    = 1.2581s2.(same as the answer before.)

    Since the sign of ω˙ is positive, our initial assumption of clockwise acceleration of the pulley is correct.

    Answer: ω˙=1.26rad/s2

  2. 2.

    From kinematics of circular motion,

    am = (aD)θ
    = ω˙Ri
    = 0.126m/s2.

    Answer: am=0.13m/s2

SAMPLE 16.19  Energy Accounting: Consider the pulley problem of Sample 16.86 again.

  1. 1.

    What percentage of the input energy (work done by the applied force F) is used in raising the mass by 1 m?

  2. 2.

    Where does the rest of the energy go? Provide an energy-balance sheet.


Solution

  1. 1.

    Let Wi and Wh be the input energy and the energy used in raising the mass by 1m, respectively. Then the percentage of energy used in raising the mass is

    % of input energy used=WhWi×100.

    Thus we need to calculate Wi and Wh to find the answer. Wi is the work done by the force F on the system during the interval in which the mass moves up by 1m. Let s be the displacement of the force F during this interval. Since the displacement is in the same direction as the force (we know it is from Sample 16.86), the input-energy is

    Wi=Fs.

    So to find Wi we need to find s.

    For the mass to move up by 1 m  the inner drum B must rotate by an angle θ where

    1m=θRi  θ=1m0.1m=10rad.

    Since the two drums, A and B, are welded together, drum A must rotate by θ as well. Therefore the displacement of force F is

    s=θRo=10rad0.2m=2m,

    and the energy input is

    Wi=Fs=20N2m=40J.

    Now, the work done in raising the mass by 1m is

    Wh=mgh=2kg9.81m/s21m=19.62J.

    Therefore, the percentage of input-energy used in raising the mass

    =19.62Nm40×100=49.05%49%.
  2. 2.

    The rest of the energy (=51%) goes in accelerating the mass and the pulley. Let us find out how much energy goes into each of these activities. Since the initial state of the system from which we begin energy accounting is not prescribed (that is, we are not given the height of the mass from which it is to be raised 1m, nor do we know the velocities of the mass or the pulley at that initial height), let us assume that at the initial state, the angular speed of the pulley is ωo and the linear speed of the mass is vo. At the end of raising the mass by 1m  from this state, let the angular speed of the pulley be ωf and the linear speed of the mass be vf. Then, the energy used in accelerating the pulley is

    (ΔEK)pulley =final kinetic energyinitial kinetic energy
    =12Iωf212Iωo2
    =12I(ωf2ωo2)assuming constant acceleration, ωf2=ωo2+2αθ, or ωf2ωo2=2αθ.
    =Iαθ(from Sample 16.88α=1.258rad/s2. )
    =1.6kgm21.258rad/s210rad
    =20.13Nm=20.13J.

    Similarly, the energy used in accelerating the mass is

    (ΔEK)mass = final kinetic energyinitial kinetic energy
    = 12mvf212mvo2
    = 12m(vf2vo22ah)
    = mah
    = 2kg0.126m/s21m
    = 0.25J.

    We can calculate the percentage of input energy used in these activities to get a better idea of energy allocation. Here is the summary table:

    Activities Energy Spent     
    in Joule as % of input energy
    In raising the mass by 1m 19.62 49.05%
    In accelerating the mass 0.25 0.62 %
    In accelerating the pulley 20.13 50.33 %
    Total 40.00 100 %

So, what would you change in the set-up so that more of the input energy is used in raising the mass? Think about what aspects of the motion would change due to your proposed design.

SAMPLE 16.20

Filename:sfig4-6-8
Figure 16.90: A uniform rod swings in the plane about its pinned end O.

Equation of motion of a swinging stick: A uniform bar of mass m and length is pinned at one of its ends O. The bar is displaced from its vertical position by an angle θ and released (Fig. 16.90).

  1. 1.

    Find the equation of motion using momentum balance.

  2. 2.

    Find the reaction at O as a function of (θ,θ˙,g,m,).


Solution First we draw a simple sketch of the given problem showing relevant geometry (Fig. 16.90(a)), and then a free-body diagram of the bar (Fig. 16.90(b)).

Filename:sfig4-6-8a
Figure 16.91:

(a) A line sketch of the swinging rod and (b) free-body diagram of the rod.

We should note for future reference that

𝝎 = ω𝒌ˆθ˙𝒌ˆ
𝝎˙ = ω˙𝒌ˆθ¨𝒌ˆ
  1. 1.

    Equation of motion using momentum balance: We can write angular momentum balance about point O as

    𝑴O=𝑯˙/O.

    Let us now calculate both sides of this equation:

    Filename:sfig4-6-8b
    Figure 16.92: Computation of 𝑯˙/O by integration over the rod.
    𝑴O = 𝒓G/O×mg(ȷˆ) (16.75)
    = 2(sinθıˆcosθȷˆ)×mg(ȷˆ)
    = 2mgsinθ𝒌ˆ.
    𝑯˙/O = ω˙𝒌ˆmr2𝑑m (16.76)
    = θ¨𝒌ˆ0s2m𝑑s(dm=m/ds)
    = mθ¨[s33|0]=m23θ¨𝒌ˆ

    Equating (16.75) and (16.76) we get

    2mgsinθ = m23θ¨
    or ω˙+3g2sinθ = 0
    or θ¨+3g2sinθ = 0. (16.77)

    Answer: θ¨+3g2sinθ=0

  2. 2.

    Reaction at O: Using linear momentum balance

    𝑭 = m𝒂G,
    where 𝑭 = Rxıˆ+(Rymg)ȷˆ,
    and 𝒂G = 2ω˙(cosθıˆ+sinθȷˆ)+2ω2(sinθıˆ+cosθȷˆ)
    = 2[(ω˙cosθω2sinθ)ıˆ+(ω˙sinθ+ω2cosθ)ȷˆ].

    Dotting both sides of 𝑭=m𝒂G with ıˆ and ȷˆ and rearranging, we get

    Rx = m2(ω˙cosθω2sinθ)
    m2(θ¨cosθθ˙2sinθ),
    Ry = mg+m2(ω˙sinθ+ω2cosθ)
    mg+m2(θ¨sinθ+θ˙2cosθ).

    Now substituting the expression for θ¨ from (16.77) in Rx and Ry, we get

    Rx = msinθ(34gcosθ+2θ˙2), (16.78)
    Ry = mg(134sin2θ)+m2θ˙2cosθ. (16.79)

    Answer: R=m(34gcosθ+2θ˙2)sinθıˆ+[mg(134sin2θ)+m2θ˙2cosθ]ȷˆ

    Check: We can check the reaction force in the special case when the rod does not swing but just hangs from point O. The forces on the bar in this case have to satisfy static equilibrium. Therefore, the reaction at O must be equal to mg and directed vertically upwards. Plugging θ=0 and θ˙=0 (no motion) in Eqn. (16.78) and (16.79) we get Rx=0 and Ry=mg, the values we expect.

SAMPLE 16.21

Filename:sfig5-5-2
Figure 16.93: A uniform rod swings in the plane about its pinned end O.

Swinging stick dynamics using moment of inertia: A uniform bar of mass m and length is pinned at one of its ends O. The bar is displaced from its vertical position by an angle θ and released (Fig. 16.93). Find the equation of motion of the stick.


margin:

Solution

We repeat the problem solved in Sample 16.88 here with just one different step of finding the rate of change of angular momentum with the help of moment of inertia formula. As usual, we first draw a free-body diagram of the bar (Fig. 16.94). We assume, 𝝎=ω𝒌ˆθ˙𝒌ˆ, and 𝝎˙=ω˙𝒌ˆθ¨𝒌ˆ

Filename:sfig5-5-2a
Figure 16.94: The free-body diagram of the rod.

We can write angular momentum balance about point O as

𝑴O=𝑯˙/O.

Let us now calculate both sides of this equation:

Filename:sfig5-5-2b
Figure 16.95: Radial and tangential components of 𝒂G. Since the radial component is parallel to 𝒓G, 𝒓G×𝒂G=24ω˙𝒌ˆ.
𝑴O = 𝒓G/O×mg(ȷˆ) (16.80)
= 2(sinθıˆcosθȷˆ)×mg(ȷˆ)
= 2mgsinθ𝒌ˆ.
𝑯˙/O = Izz/G𝝎˙+𝒓G×m𝒂G (16.81)
= m212ω˙𝒌ˆ+𝒓G×m(ω˙𝒌ˆ×𝒓Gω2𝒓G𝒂G)
= m212ω˙𝒌ˆ+m24ω˙𝒌ˆ
= m23ω˙𝒌ˆ (16.82)

where the last step, 𝒓G×m𝒂G=m24ω˙𝒌ˆ, should be clear from Fig. 16.95. Equating (16.75) and (16.76) we get

2mgsinθ = m23ω˙
or ω˙+3g2sinθ = 0
or θ¨+3g2sinθ = 0. (16.83)

Answer: θ¨+3g2sinθ=0

SAMPLE 16.22   Numerical solution of the swinging stick motion: For the swinging stick considered in Samples 16.88 or 16.97, find the time that the rod takes to fall from θ=π/2 to θ=0 if it is released from rest at θ=π/2?

Solution

The given initial angle π/2 is a big value of θ – big in that we cannot assume sinθθ (obviously 11.5708). Therefore we may not use the linearized equation (16.86) to solve for t explicitly. We have to solve the full nonlinear equation (16.83) to find the required time. Unfortunately, we cannot get a closed form solution of this equation using mathematical skills you have at this level. Therefore, we resort to numerical integration of this equation.

For numerical integration, we need to first write the given differential equation as a set of first order ordinary differential equations. To do so, we introduce ω as a new variable and rewrite eqn. (16.83) as

θ˙ = ω (16.84)
ω˙ = 3g2sinθ (16.85)

Now we need to specify the initial conditions and the time duration for integration, and solve the equations using some ODE solver program. Here is a pseudo-code that lists the steps:

   g = 9.81,   L = 1      % define constants
   ODES = { thetadot = omega
            omegadot = -3*g/(2*L) * sin(theta) }
    ICs = { theta_0 = pi/2
            omega_0 = 0 }
   solve ODES with ICs for t = 0 to 4 s
   plot theta vs t and plot omega vs t

The results obtained from the numerical solution are shown in Fig. 16.96.

Filename:sfig7-4-2
Figure 16.96: Numerical solution is shown by plotting θ and ω against time.

The problem of finding the time taken by the bar to fall from θ=π/2 to θ=0 numerically is nontrivial. It is called a boundary value problem. We have only illustrated how to solve initial value problems. However, we can get a fairly good estimate of the time from the solution obtained.

We first plot θ against time t as shown in fig. 16.97. We find the values of t and the corresponding values of θ that bracket θ=0. Now, we can use linear interpolation to find the value of t at θ=0. Proceeding this way, we get t=0.48 (seconds), a little more than we get from the linear ODE in sample 16.97 (t=0.41).

Comments: Additionally, we can also get the value of θ˙ω when θ=0 using similar interpolation. In fact, from the ω vs t plot, we find that at t=0.48s, ω=5.42rad/s. How does this result compare with the analytical value of ω from sample 16.97 (which did not depend on the small angle approximation)? Well, we found that

ω=3g=39.81m/s21m=5.4249s1.

Thus, we get a fairly accurate value from numerical integration.

Filename:pfigure-blue-99-1
Figure 16.97: Time t corresponding to θ=0.

SAMPLE 16.23

Filename:pfigure-blue-47-2
Figure 16.98: Work done by the force of gravity in moving from G to G 𝑭𝑑𝒓=mgȷˆhȷˆ=mgh.
Filename:pfigure-blue-89-1
Figure 16.99: The infinitesimal mass dm considered in the calculation of EK.

The swinging stick dynamics with energy balance: Consider the same swinging stick as in Sample 16.88. The stick is, again, displaced from its vertical position by an angle θ and released (See Fig. 16.90).

  1. 1.

    Find the equation of motion using energy balance.

  2. 2.

    What is θ˙ at θ=0 if θ(t=0)=π/2?

  3. 3.

    Find the period of small oscillations about θ=0.

Solution

  1. 1.

    Equation of motion using energy balance: We use the power equation, EK˙=P, to derive the equation of motion of the bar. Now, the kinetic energy is given by

    EK=12mv2𝑑m

    where v is the speed of the infinitesimal mass element dm. Referring to fig. 16.99, we can write, dm=(m/)ds, and v=ωsθ˙s. Thus,

    EK = 120θ˙2s2m𝑑s
    = mθ˙220s2𝑑s
    = 16m2θ˙2

    and, therefore,

    EK˙=ddt(16m2ω2)=13m2ωω˙=13m2θ˙θ¨.

    Calculation of power (P): There are only two forces acting on the bar, the reaction force, 𝑹(=Rxıˆ+Ryȷˆ) and the force due to gravity, mgȷˆ. Since the support point O does not move, no work is done by 𝑹. Therefore,

    W = Work done by gravity force in moving from G to G.=mgh

    Note that the negative sign stands for the work done against gravity.

    Now,

    h=OGOG′′=22cosθ=2(1cosθ).

    Therefore,

    W = mg2(1cosθ)
    and P = W˙=dWdt=mg2sinθθ˙.

    Equating EK˙ and P we get

    mg2sinθθ˙ = 13m2θ˙θ¨
    or θ¨+3g2sinθ = 0.

    Answer: θ¨+3g2sinθ=0

    This equation is, of course, the same as we obtained using balance of angular momentum in Sample 16.88.

  2. 2.

    Find ω at θ=0: We are given that at t=0,θ=π/2 and θ˙ω=0 (released from rest). This position is (1) shown in Fig. 16.100. In position (2) θ=0, i.e., the rod is vertical. Since there are no dissipative forces, the total energy of the system remains constant. Therefore, taking datum for potential energy as shown in Fig. 16.100, we may write

    EK10+V1 = EK2+V20
    or mg2 = 12mv2𝑑m
    = 16m2ω2(see part (a))
     ω = ±3g
    Filename:pfigure-blue-35-2
    Figure 16.100: The total energy between positions (1) and (2) is constant.

    Answer: ω=±3g

  3. 3.

    Period of small oscillations: The equation of motion is

    θ¨+3g2sinθ=0.

    For small θ, sinθθ

    θ¨+3g2θ=0 (16.86)
    or θ¨+λ2θ=0
    where λ2=3g2.

    Therefore,

    the circular frequency = λ=3g2,
    and the time period T = 2πλ=2π23g.

    Answer: T=2π23g

    [Say for g=9.81m/s2,=1mwe getT4=π22319.81s=0.4097s]

SAMPLE 16.24   The swinging stick with a destabilizing torque. Consider the swinging stick of Sample 16.88 once again.

  1. 1.

    Find the equation of motion of the stick, if a torque 𝑴=M𝒌ˆ is applied at end O and a force 𝑭=Fıˆ is applied at the other end A.

  2. 2.

    Take F=0 and M=Cθ. For C=0 you get the equation of free oscillations obtained in Sample 16.88 or 16.97. For small C, does the period of the pendulum increase or decrease?

  3. 3.

    What happens if C is big?


Solution

  1. 1.

    A free-body diagram of the bar is shown in Fig. 16.101. Once again, we can use 𝑴O=𝑯˙/O to derive the equation of motion as in Sample 16.88. We calculated 𝑴O and 𝑯˙/O in Sample 16.88. Calculation of 𝑯˙/O remains the same in the present problem. We only need to recalculate 𝑴O.

    𝑴O = M𝒌ˆ+𝒓G/O×mg(ȷˆ)+𝒓A/O×𝑭
    = M𝒌ˆ2mgsinθ𝒌ˆ+Fcosθ𝒌ˆ
    = (M+Fcosθ2mgsinθ)𝒌ˆ
    and
    𝑯˙/O = mθ¨23𝒌ˆ(see Sample 16.88)
    Filename:pfigure4-3Dpend
    Figure 16.101: Free-body diagram of the bar with applied torque 𝑴 and force 𝑭

    Therefore, from 𝑴O=𝑯˙/O

    M+Fcosθ2mgsinθ = mθ¨23
     θ¨+3g2sinθ3Fmcosθ3Mm2 = 0.

    Answer: θ¨+3g2sinθ3Fmcosθ3Mm2=0

  2. 2.

    Now, setting F=0 and M=Cθ we get

    θ¨+3g2sinθ3Cθm2=0 (16.87)

    Numerical Solution: We can numerically integrate (16.87) just as in the previous Sample to find θ(t). Here is the pseudo-code that can be used for this purpose.

       g = 9.81, L = 1      % specify parameters
       m = 1, C = 4
       ODES = { thetadot = omega
                omegadot = -(3*g/(2*L)) * sin(theta)
                          + 3*C/(m*L^2) * theta      }
       ICs =  { thetazero = pi/20
                omegazero = 0  }
       solve ODES with ICs until t = 10
    

    Using this pseudo-code, we find the response of the pendulum. Figure 16.102 shows different responses for various values of C. Note that for C=0, it is the same case as unforced bar pendulum considered above.

    From Fig. 16.102 it is clear that the bar has periodic motion for small C, with the period of motion increasing with increasing values of C. It makes sense if you look at Eqn. (16.87) carefully. Gravity acts as a restoring force while the applied torque acts as a destabilizing force. Thus, with the resistance of the applied torque, the stick swings more sluggishly making its period of oscillation bigger.

    Filename:pfigure-s94h6p3
    Figure 16.102: θ(t) with applied torque M=Cθ for C=0, 1, 2, 4, 4.905, 5. Note that for small C the motion is periodic but for large C (C4.4) the motion becomes aperiodic.
  3. 3.

    From Fig. 16.102, we see that at about C4.9 the stability of the system changes completely. θ(t) is not periodic anymore. It keeps on increasing at a faster and faster rate, that is, the bar makes complete loops about point O with ever increasing speed. Does it make physical sense? Yes, it does. As the value of C is increased beyond a certain value (can you guess the value?), the applied torque overcomes any restoring torque due to gravity. Consequently, the bar is forced to rotate continuously in the direction of the applied force.

SAMPLE 16.25

Filename:pfigure-s94h6p4
Figure 16.103:

A torsional pendulum with linear springs: A uniform rigid bar of mass m=2kg and length =1m is pinned at one end and connected to two springs, each with spring constant k, at the other end. The bar is tweaked slightly from its vertical position. It then oscillates about its original position. The bar is timed for 20 full oscillations which take 12.5 seconds. Ignore gravity.

  1. 1.

    Find the equation of motion of the rod.

  2. 2.

    Find the spring constant k.

  3. 3.

    What should be the spring constant of a torsional spring if the bar is attached to one at the bottom and has the same oscillating motion characteristics?


Solution

  1. Filename:pg84-3
    Figure 16.104:
  2. 1.

    Refer to the free-body diagram in figure 16.104. Angular momentum balance for the rod about point O gives

    𝑴O=𝑯˙/O
    where𝑴O = 2kxsinθcosθ𝒌ˆ
    = 2k2sinθcosθ𝒌ˆ,
    and𝑯˙/O = IzzOθ¨𝒌ˆ=13m2IzzOθ¨𝒌ˆ.

    Thus

    13m2θ¨=2k2sinθcosθ.

    However, for small θ,cosθ1 and sinθθ,

      θ¨+6k2m2θ=0. (16.88)

    Answer: θ¨+6kmθ=0

  3. 2.

    Comparing Eqn. (16.88) with the standard harmonic oscillator equation x¨+λ2x=0, we get

    angular frequencyλ = 6km,
    and the time periodT = 2πλ
    = 2πm6k.

    From the measured time for 20 oscillations, the time period (time for one oscillation) is

    T=12.520s=0.625s

    Now equating the measured T with the derived expression for T we get

    2πm6k = 0.625s
     k = 4π2m6(0.625s)2
    = 4π22kg6(0.625s)2
    = 33.7N/m.

    Answer: k=33.7N/m

  4. 3.

    If the two linear springs are to be replaced by a torsional spring at the bottom, we can find the spring constant of the torsional spring by comparison. Let ktor be the spring constant of the torsional spring. Then, as shown in the free-body diagram (see figure 16.105), the restoring torque applied by the spring at an angular displacement θ is ktorθ. Now, writing the angular momentum balance about point O, we get

    Filename:pfigure-blue-36-1
    Figure 16.105:
    𝑴O = 𝑯˙/O
    ktorθ𝒌ˆ = IzzO(θ¨𝒌ˆ)
     θ¨+ktorIzzOθ = 0.

    Comparing with the standard harmonic equation, we find the angular frequency

    λ=ktorIzzO=ktor13m2.

    If this system has to have the same period of oscillation as the first system, the two angular frequencies must be equal, i.e.,

    ktor13m2 = 6km
     ktor = 6k132=2k2
    = 2(33.7N/m)(1m)2
    = 67.4Nm.

    Answer: ktor=67.4Nm

SAMPLE 16.26

Filename:pfigure-s95q8
Figure 16.106:

A spring loaded seesaw: A kid, modelled as a point mass with m=10kg, is sitting at end B of a rigid rod AB of negligible mass. The rod is supported by a spring at end A and a pin at point O. The system is in static equilibrium when the rod is horizontal. Someone pushes the kid vertically downwards by a small distance y and lets go. Given that AB=3m,AC=0.5m,k=1kN/m; find

  1. 1.

    the unstretched (relaxed) length of the spring,

  2. 2.

    the equation of motion (a differential equation relating the position of the mass to its acceleration) of the system, and

  3. 3.

    the natural frequency of the system.

If the rod is pinned at the midpoint instead of at O, what is the natural frequency of the system? How does the new natural frequency compare with that of a mass m simply suspended by a spring with the same spring constant?


Solution

  1. 1.
    Filename:pfigure4-rpi
    Figure 16.107: Free-body diagrams

    Static Equilibrium: The FBD of the (rod + mass) system is shown in Fig. 16.107. Let the stretch in the spring in this position be yst and the relaxed length of the spring be 0. The balance of angular momentum about point O gives:

    𝑴/o=𝑯˙/o = 𝟎(no motion)
     (kyst)d1(mg)d2 = 0
     yst = mgkd2d1
    = 10kg9.8m/s221000N/m=0.196m
    Therefore, 0 = ACyst
    = 0.5m0.196m=0.304m.

    Answer: 0=30.4cm

  2. 2.

    Equation of motion: As point B gets displaced downwards by a distance y, point A moves up by a proportionate distance ya. From geometry, margin: Here, we are considering a very small y so that we can ignore the arc the point mass B moves on and take its motion to be just vertical (i.e., sinθθ for small θ).

    y d2θ  θ=yd2
    ya d1θ=d1d2y

    Therefore, the total stretch in the spring, in this position,

    Δy=ya+yst=d1d2y+d2d1mgk

    Now, Angular Momentum Balance about point O gives:

    𝑴/o = 𝑯˙/o
    𝑴/o = 𝒓B×mgȷˆ+𝒓A×kΔyȷˆ (16.89)
    = (d2mgd1kΔy)𝒌ˆ
    𝑯˙/o = 𝒓B×m𝒂=𝒓B×my¨ȷˆ (16.90)
    = d2my¨𝒌ˆ (16.91)

    Equating (16.89) and (16.91) we get

    d2mgd1kΔy = d2my¨
    or d2mgd1k(d1d2y+d2mgd1k) = d2my¨
    or d2mgkd12d2yd2mg = d2my¨
    or y¨+km(d1d2)2y = 0

    Answer: y¨+km(d1d2)2y=0

  3. 3.

    The natural frequency of the system: We may also write the previous equation as

    y¨+λ2y=0whereλ2=kmd12d22. (16.92)

    Substituting d1= and d2=2 in the expression for λ we get the natural frequency of the system

    λ=12km=121000N/m10kg=5s1

    Answer: λ=5s1

  4. 4.

    Comparison with a simple spring mass system:

    Filename:pfigure-blue-86-1
    Figure 16.108:

    When d1=d2, the equation of motion (16.92) becomes

    y¨+kmy=0

    and the natural frequency of the system is simply

    λ=km

    which corresponds to the natural frequency of a simple spring mass system shown in Fig. 16.108.

    In our system (with d1=d2 ) any vertical displacement of the mass at B induces an equal amount of stretch or compression in the spring which is exactly the case in the simple spring-mass system. Therefore, the two systems are mechanically equivalent. Such equivalences are widely used in modeling complex physical systems with simpler mechanical models.

Problems for 16.4 Dynamics

Preparatory Problems

16.4.1  An object consists of a massless bar with two attached masses m1 and m2. The object is hinged at O.

  1. (a)

    What is the moment of inertia of the object about point O (IzzO)? Answer: Izzcm=m112+m22.

  2. (b)

    Given θ, θ˙, and θ¨, what is 𝑯/O, the angular momentum about point O? Answer: H/O=θ˙(m112+m22)kˆ.

  3. (c)

    Given θ, θ˙, and θ¨, what is 𝑯˙/O, the rate of change of angular momentum about point O? Answer: H˙/O=θ¨(m112+m22)kˆ.

  4. (d)

    Given θ, θ˙, and θ¨, what is T, the total kinetic energy? Answer: EK=θ˙22(m112+m22)kˆ.

  5. (e)

    Assume that you don’t know θ, θ˙ or θ¨ but you do know that F1 is applied to the rod, perpendicular to the rod at m1. What is θ¨? (Neglect gravity.) Answer: θ¨=F11m112+m22

  6. (f)

    If F1 were applied to m2 instead of m1, would θ¨ be bigger or smaller? Answer: θ¨ will be bigger if F1 is applied to m2 instead of m1.

Filename:pfigure-blue-88-1
Figure 16.109:

16.4.2  A uniform circular disc rotates at constant angular speed ω about the origin, which is also the center of the disc. It’s radius is R. It’s total mass is M.

  1. (a)

    What is the total force and moment required to hold it in place (use the origin as the reference point of angular momentum and torque). Answer: F=0N.

  2. (b)

    What is the total kinetic energy of the disk? Answer: EK=14MR2ω2.

Filename:pfigure4-2-rp3
Figure 16.110:

16.4.3  The hinged disk of mass m (uniformly distributed) is acted upon by a force P shown in the figure. Determine the initial angular acceleration and the reaction forces at the pin O. Answer: θ¨0=PmR and 𝑹1=P(sin30ıˆcos30ȷˆ).

Filename:pfigure-s94q6p1
Figure 16.111:

16.4.4  A motor turns a bar. A uniform bar of length and mass m is turned by a motor whose shaft is attached to the end of the bar at O. The angle that the bar makes (measured counter-clockwise) from the positive x axis is θ=2πt2/s2. Neglect gravity.

  1. (a)

    Draw a free-body diagram of the bar.

  2. (b)

    Find the force acting on the bar from the motor and hinge at t=1s. Answer: RBM=4πm2(4πıˆȷˆ).

  3. (c)

    Find the torque applied to the bar from the motor at t=1s. Answer: M=4πrad/s23m2.

  4. (d)

    What is the power produced by the motor at t=1s? Answer: P=16π2s33m2.

Filename:pfigure-spinningbrick
Figure 16.112:

16.4.5  A physical pendulum. A swinging stick is sometimes called a ‘physical’ pendulum. Take the ‘body’, the system of interest, to be the whole stick.

  1. (a)

    Draw a free-body diagram of the system.

  2. (b)

    Write the equation of angular momentum balance for this system about point O.

  3. (c)

    Evaluate the left-hand-side as explicitly as possible in terms of the forces showing on your Free-Body Diagram.

  4. (d)

    Evaluate the right hand side as completely as possible. You may use the following facts:

    𝒗 = lθ˙cosθȷˆ+lθ˙sinθıˆ
    𝒂 = lθ˙2[cosθıˆ+sinθȷˆ]
    +lθ¨[cosθȷˆsinθıˆ]

    where is the distance along the pendulum from the top, θ is the angle by which the pendulum is displaced counter-clockwise from the vertically down position, ıˆ is vertically down, and ȷˆ is to the right. You will have to set up and evaluate an integral. Answer: MO=2mgsinθkˆ=H˙/O=m23θ¨kˆ.

Filename:pfigure-s94f1p2
Figure 16.113:

16.4.6   A uniform one meter bar is hung from a hinge that is at the end. It is allowed to swing freely. g=10m/s2.

  1. (a)

    What is the period of small oscillations for this pendulum? Answer: T=1.62s.

  2. (b)

    Suppose the rod is hung 0.4m from one end. What is the period of small oscillations for this pendulum? Can you explain why it is longer or shorter than when it is hung by its end? Answer: The period of oscillations is longer if the rod is hung at 0.4m because the torque due to weight of the portion above the hinge negates the torque due to the portion below the hinge. One way to see it is that if the bar is hinged at mid, net torque about the hinge becomes zero leading to infinitely large period of oscillation.

More-Involved Problems

16.4.7  Motor turns a dumbbell. Two uniform bars of length and mass m are welded at right angles. At the ends of the horizontal bar are two more masses m. The bottom end of the vertical rod is attached to a hinge at O where a motor keeps the structure rotating at constant rate ω (counter-clockwise). What is the net force and moment that the motor and hinge cause on the structure at the instant shown? Answer: Fx=0,Fy=7mω22, 𝑴O=0𝒌ˆ.

Filename:s97f2
Figure 16.114:

16.4.8  The structure shown in the figure consists of two point masses connected by three rigid, massless rods such that the whole structure behaves like a rigid body. The structure rotates counterclockwise at a constant rate of 60rpm. At the instant shown, find the force in each rod. Answer: F1=0N, F2=2π2N, and F3=22π2N.

Filename:pg92-2
Figure 16.115:

16.4.9  Balancing a system of rotating particles. A wire frame structure is made of four concentric loops of massless and rigid wires, connected to each other by four rigid wires presently coincident with the x and y axes. Three masses, m1=200grams, m2=150grams and m3=100grams, are glued to the structure as shown in the figure. The structure rotates counter-clockwise at a constant rate θ˙=5rad/s. There is no gravity.

  1. (a)

    Find the net force exerted by the structure on the support at the instant shown. Answer: (a) 𝑭=0.33Nıˆ0.54Nȷˆ.

  2. (b)

    You are to put a mass m at an appropriate location on the third loop so that the net force on the support is zero. Find the appropriate mass and the location on the loop. Answer: The net force on the support becomes zero if we put m=0.84kg on 0.3m radius loop at θ=302.33.

Filename:summer95p2-3
Figure 16.116:

16.4.10  Motor turns a bent bar. Two uniform bars of length and uniform mass m are welded at right angles. One end is attached to a hinge at O where a motor keeps the structure rotating at a constant rate ω (counterclockwise). What is the net force and moment that the motor and hinge cause on the structure at the instant shown.

  1. (a)

    neglecting gravity Answer: Net force: 𝑭net=(3mω2L2)ıˆ(mω2L2)ȷˆ, Net moment: 𝑴net=𝟎.

  2. (b)

    including gravity. Answer: Net force: 𝑭net=(3mω2L2)ıˆ+(2mgmω2L2)ȷˆ, Net moment: 𝑴net=3mgL2𝒌ˆ.

Filename:p-f96-p3-3
Figure 16.117:

A bent bar is rotated by a motor.

16.4.11   A uniform disk of mass M and radius R rotates about a hinge O in the xy-plane. A point mass m is fixed to the disk at a distance R/2 from the hinge. A motor at the hinge drives the disk/point mass assembly with constant angular acceleration α. What torque at the hinge does the motor supply to the system? Answer: T=αR24(2M+m).

16.4.12  A rigid rod of length and total mass m is held fixed at one end and whirled around in circular motion at a constant rate ω in the horizontal plane. Ignore gravity.

  1. (a)

    Find the tension in the rod as a function of r, the radial distance from the center of rotation to any desired location on the rod. Answer: T(r)=mω22L(L2r2)

  2. (b)

    Where does the maximum tension occur in the rod? Answer: at r=0; i.e., at the center of rotation.

  3. (c)

    At what distance from the center of rotation does the tension drop to half its maximum value? Answer: r=L/2.

16.4.13  A thin uniform circular disc of mass M and radius R rotates in the xy plane about its center of mass point O. Driven by a motor, it has rate of change of angular speed proportional to angular position, α=kθ5/2. The disc starts from rest at θ=0.

  1. (a)

    What is the rate of change of angular momentum about the origin at θ=π3rad? Answer: H˙/O=(kπ3)2.5(MR22)kˆ.

  2. (b)

    What is the torque of the motor at θ=π3rad? Answer: T=k(π3)2.5(MR22).

  3. (c)

    What is the total kinetic energy of the disk at θ=π3rad? Answer: EK=k2(π3)3.5(MR22).

16.4.14  Neglecting gravity, calculate α=ω˙=θ¨ at the instant shown for the system in the figure. Answer: θ¨=100rev/sec.

Filename:s97p3-3
Figure 16.118:

16.4.15  The uniform square shown is released from rest at t=0. What is α=ω˙=θ¨ immediately after release? Answer: θ¨=3g2Lsinθ

Filename:pfigure-s94h8p1
Figure 16.119:

16.4.16  Acceleration of a trap door. A uniform bar AB of mass m and a ball of the same mass are released from rest from the same horizontal position. The bar is hinged at end A. There is gravity.

  1. (a)

    Which point on the rod has the same acceleration as the ball, immediately after release. Answer: Point at 2L/3 from A

  2. (b)

    What is the reaction force on the bar at end A just after release? Answer: mg/4 directed upwards.

Filename:sfig4-7-DH1
Figure 16.120:

16.4.17  A disk with radius R has a string wrapped around it which is pulled with a force F. The disk is free to rotate about the axis through O normal to the page. The moment of inertia of the disk about O is Io. A point A is marked on the string. Given that xA(0) = 0 and that x˙A(0) = 0, what is xA(t)? Answer: xA(t)=Ft22(m+IR2)

Filename:sfig4-7-DH2
Figure 16.121:

16.4.18   A uniform stick with length and mass Mo is welded to a pulley hinged at the center O. The pulley has negligible mass and radius Rp. A string is wrapped many times around the pulley. At time t=0, the pulley, stick, and string are at rest and a force F is suddenly applied to the string. How long does it take for the pulley to make one full revolution? Ignore gravity. Answer: Trev=2MO2π3FRp.

Filename:sfig4-7-DH3
Figure 16.122:

String wraps around a pulley with a stick glued to it.

Gears, Belt Drives, and Gear Trains

16.4.19  Constant speed gear train. Gear A is connected to a motor (not shown) and gear B, which is welded to gear C, is connected to a taffy-pulling mechanism. Assume you know the torque Minput=MA and angular velocity ωinput=ωA of the input shaft. Assume the bearings and contacts are frictionless.

  1. (a)

    What is the input power? Answer: Pin=MAωA.

  2. (b)

    What is the output power? Answer: Pout=MAωA.

  3. (c)

    What is the output torque Moutput=MC, the torque that gear C applies to its surroundings in the clockwise direction? Answer: MC=MARBRA

Filename:sfig4-7-1
Figure 16.123:

16.4.20   At the input to a gear box a 100lbf force is applied to gear A. At the output, the machinery (not shown) applies a force of FB to the output gear. Gear A rotates at constant angular rate ω=2rad/s, clockwise.

  1. (a)

    What is the angular speed of the right gear? Answer: ωB=ωRARC

  2. (b)

    What is the velocity of point P? Answer: VP=ωRARBRC

  3. (c)

    What is FB? Answer: FB=FARCRB

  4. (d)

    If the gear bearings had friction, would FB have to be larger or smaller in order to achieve the same constant velocity? Answer: FB will have to be smaller.

  5. (e)

    If instead of applying a 100lbf to the left gear it is driven by a motor (not shown) at constant angular speed ω, what is the angular speed of the right gear? Answer: The left gear is constrained to rotate at ωB=ωRARC.

Filename:sfig4-7-1a
Figure 16.124:

Two gears with end loads.

16.4.21  A bevel gear system. A bevel-type gear system, shown in the figure, is used to transmit power between two shafts that are perpendicular to each other. The driving gear has a mean radius of 50mm and rotates at a constant speed ω=150rpm. The mean radius of the driven gear is 80mm and the driven shaft is expected to deliver a torque of Mout=25Nm. Assuming no power loss, find the input torque supplied by the driving shaft. Answer: Min=15.625Nm.

Filename:sfig4-7-1b
Figure 16.125:

A bevel gear

16.4.22  Belt drives are used to transmit power between parallel shafts. Two parallel shafts, 3m apart, are connected by a belt passing over the pulleys A and B fixed to the two shafts. The driver pulley A rotates at a constant 200rpm. The speed ratio between the pulleys A and B is 1:2.5. The input torque is 350Nm. Assume no loss of power between the two shafts.

  1. (a)

    Find the input power. Answer: Pin=7.33kilo-watts

  2. (b)

    Find the rotational speed of the driven pulley B Answer: 500rpm.

  3. (c)

    Find the output torque at B. Answer: Mout=140Nm

Filename:sfig4-7-2
Figure 16.126:

16.4.23  In the belt drive system shown, assume that the driver pulley rotates at a constant angular speed ω. If the motor applies a constant torque MO on the driver pulley, show that the tensions in the two parts, AB and CD, of the belt must be different. Which part has a greater tension? Does your conclusion about unequal tension depend on whether the pulley is massless or not? Assume any dimensions you need. Answer: Mout=140Nm

Filename:sfig4-7-2a
Figure 16.127:

16.4.24  Two racks connected by three constant rate gears. A 100lbf force is applied to one rack. At the output, the machinery (not shown) applies a force of FB to the other rack.

  1. (a)

    Assume the gear-train is spinning at constant rate and is frictionless. What is FB? Answer: FB=100lbf.

  2. (b)

    If the gear bearings had friction would that increase or decrease FB to achieve the same constant rate?

  3. (c)

    If instead of applying a 100lbf to the left rack it is driven by a motor (not shown) at constant speed v, what is the speed of the right rack? Answer: vright=v.

  4. (d)

    If the angular velocity of the gear is increasing at rate α does this increase or decrease FB at the given ω.

Filename:sfig4-7-2b
Figure 16.128:

Two racks connected by three gears.

16.4.25  -3-D accelerating gear train. This is really a 2-D problem; each gear turns in a different parallel plane. Shaft B is rigidly connected to gears G4 and G5. G3 meshes with gear G6. Gears G6 and G5 are both rigidly attached to shaft AD. Gear G5 meshes with G2 which is welded to shaft A. Shaft A and shaft B spin independently. Assume you know the torque Minput, angular velocity ωinput and the angular acceleration αinput of the input shaft. Assume the bearings and contacts are frictionless.

  1. (a)

    What is the input power? Answer: Pin=2500πNm/s.

  2. (b)

    What is the output power? Answer: Pout=2500πNm/s.

  3. (c)

    What is the angular velocity ωoutput of the output shaft? Answer: ωoutput=25πrad/s.

  4. (d)

    What is the output torque Moutput? Answer: Moutput=100Nm.

Filename:sfig4-7-2c
Figure 16.129:

A 3-D gear train.

16.4.26  Gear A with radius RA=400mm is rigidly connected to a drum B with radius RB=200mm. The combined moment of inertia of the gear and the drum about the axis of rotation is Izz=0.5kgm2. Gear A is driven by gear C which has radius RC=300mm. As the drum rotates, a 5kg mass m is pulled up by a string wrapped around the drum. At the instant of interest, The angular speed and angular acceleration of the driving gear are 60rpm and 12rpm/s, respectively. Find the acceleration of the mass m. Answer: am=0.188m/s2

Filename:pfigure-s94q7p1
Figure 16.130:

16.4.27  A spindle and pulley arrangement is used to hoist a 50kg mass as shown in the figure. Assume that the pulley is to be of negligible mass. When the motor is running at a constant 100rpm,

  1. (a)

    Find the velocity of the mass at B. Answer: VB=3.77m/s.

  2. (b)

    Find the tension in strings AB and CD. Answer: TAB=500N, TCD=1000N.

Filename:pfigure-s94h7p2
Figure 16.131:

16.4.28  Accelerating rack and pinion. The two gears shown are welded together and spin on a frictionless bearing. The inner gear has radius 0.5m and negligible mass. The outer disk has 1m radius and a uniformly distributed mass of 0.2kg. They are loaded as shown with the force F=20N on the massless rack which is held in place by massless frictionless rollers. At the time of interest the angular velocity is ω=2rad/s (though ω is not constant). The point P is on the disk a distance 1m from the center. At the time of interest, point P is on the positive y axis.

  1. (a)

    What is the speed of point P? Answer: |vp|=2m/s.

  2. (b)

    What is the velocity of point P? Answer: vp=2ıˆm/s.

  3. (c)

    What is the angular acceleration α of the gear? Answer: α=100rad/s2.

  4. (d)

    What is the acceleration of point P? Answer: ap=(100ıˆ4ȷˆ)m/s2.

  5. (e)

    What is the magnitude of the acceleration of point P? Answer: |ap|=100.08m/s2.

  6. (f)

    What is the rate of increase of the speed of point P? Answer: d|va|dt=100m/s2.

Filename:ch4-7-c
Figure 16.132:

Accelerating rack and pinion

16.4.29  Two racks connected by a gear. A 100lbf force is applied to one rack. At the output the machinery (not shown) applies a constant force FB to the other rack.

  1. (a)

    Assume the gear is spinning at constant rate and is frictionless. What is FB? Answer: FB=100lbf.

  2. (b)

    If the gear bearing had friction, would that increase or decrease FB to achieve the same constant rate? Answer: If gear bearing has friction, FB has to increase.

  3. (c)

    If the angular velocity of the gear is increasing at rate α, does this increase or decrease FB at the given ω. Answer: FB will decrease if the angular velocity of the gear is increasing.

  4. (d)

    If the output load FB is given then the motion of the machine can be found, assuming friction is negligible. Assume that the machine starts from rest with a given output load. So long as rack B moves down and the output force FB is positive, the output power is positive. The machine starts from rest and runs for a fixed amount of time T. For what value of FB is the output work maximum? Once you find an answer, can you find an intuitive explanation for the answer? Answer: For work to be maximum, FB=50lbf. Initially work increases with FB but too large FB leads to small ω0 thus giving smaller net work.

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

Two racks connected by a gear.

16.4.30  A rack and pinion constrained by a linear spring. Neglect gravity. The spring is relaxed when the angle θ=0. Assume the system is released from rest at θ=θ0. What is the acceleration of the point P at the end of the stick when θ=0? Answer in terms of any or all of m, R, , θ0, k, ıˆ, and ȷˆ. [Hint: There are several steps of reasoning required. You might want to draw FBD(s), use angular momentum balance, set up a differential equation, solve it, plug values into this solution, and use the result to find the quantities of interest. Answer: aP=3θ02R2kmȷˆ

Filename:Mikef91p3
Figure 16.134:

Pendulum Problems

16.4.31  A small particle of mass m is attached to the end of a thin rod of mass M (uniformly distributed), which is pinned at hinge O, as depicted in the figure.

  1. (a)

    Obtain the equation of motion governing the rotation θ of the rod. Answer: ((M+4m)29)θ¨+((M+4m)g6)sinθ=0.

  2. (b)

    What is the natural frequency of the system for small oscillations θ? Answer: f=12π3g(M+4m)2(M+4m)

Filename:pfigure-s94h7p3
Figure 16.135:

16.4.32  A rigid massless rod has two equal masses mB and mC (mB=mC=m) attached to it at distances 2 and 3, respectively, measured along the rod from a frictionless hinge located at a point A. The rod swings freely from the hinge. There is gravity. Let ϕ denote the angle of the rod measured from the vertical. Assume that ϕ and ϕ˙ are known at the instant of interest.

  1. (a)

    What is ϕ¨? Find ϕ¨ in terms of m, , g, ϕ and ϕ˙. Answer: ϕ¨=5g13sinϕ.

  2. (b)

    What is the force of the hinge on the rod? Solve in terms of m, , g, ϕ, ϕ˙, ϕ¨ and any unit vectors you may need to define. Answer: F=5m(ϕ˙2cosϕ+ϕ¨sinϕ)ıˆ2mgıˆ+5ml(ϕ¨cosϕϕ˙2sinϕ)ȷˆ.

  3. (c)

    Would you get the same answers if you put a mass 2m at 2.5? Why or why not? Answer: Answer will change because although total moment of force about the hinge stays same but moment of inertia reduces. In this case we get ϕ¨=5g12.5sinϕ.

Filename:p-f96-p2-3
Figure 16.136:

16.4.33  For the pendula in the figure :

  1. (a)

    Without doing any calculations, try to figure out the relative durations of the periods of oscillation for the five pendula (i.e. the order, slowest to fastest) Assume small angles of oscillation.

  2. (b)

    Calculate the period of small oscillations. [Hint: use balance of angular momentum about the point 0]. Answer:

    1. (b)

      T=2.29s

    2. (e)

      T=1.99s

    (b) has a longer period than (e) does since in (b) the moment of inertia about the center of mass (located at the same position as the mass in (e)) is non-zero.

  3. (c)

    Rank the relative duration of oscillations and compare to your intuitive solution in part (a), and explain in words why things work the way they do.

Filename:s97p3-2
Figure 16.137:

16.4.34  A pegged compound pendulum. A uniform bar of mass m and length hangs from a peg at point C and swings in the vertical plane about an axis passing through the peg. The distance d from the center of mass of the rod to the peg can be changed by putting the peg at some other point along the length of the rod.

  1. (a)

    Find the angular momentum of the rod about point C. Answer: 𝑯c=(m212+md2)θ˙𝒌ˆ.

  2. (b)

    Find the rate of change of angular momentum of the rod about C. Answer: 𝑯˙c=(m212+md2)θ¨𝒌ˆ.

  3. (c)

    How does the period of the pendulum vary with d? Show the variation by plotting the period against d. [Hint, you must first find the equations of motion, linearize for small θ, and then solve.] Answer: T=2πg/112(d/)+d

  4. (d)

    Find the total energy of the rod (using point C as a datum for potential energy). Answer: ET=12(m212+md2)θ˙2mgdcosθ.

  5. (e)

    Find θ¨ when θ=π/6. Answer: θ¨=6gd12d2+2.

  6. (f)

    Find the reaction force on the rod at C, as a function of m, d, , θ, and θ˙. Answer: 𝑹c=m(12gd12d2+2cosθsinθ+θ˙2dsinθ)ıˆ+(mg12mgd12d2+2sinθ2+mθ˙2dcosθ)ȷˆ.

  7. (g)

    For the given rod, what should be the value of d (in terms of ) in order to have the fastest pendulum? Answer: d=0.29

  8. (h)

    Test of Schuler’s pendulum. The pendulum with the value of d obtained in (g) is called the Schuler’s pendulum. It is not only the fastest pendulum but also the “most accurate pendulum”. The claim is that even if d changes slightly over time due to wear at the support point, the period of the pendulum does not change much. Verify this claim by calculating the percent error in the time period of a pendulum of length =1m under the following three conditions: (i) initial d=0.15m and after some wear d=0.16m, (ii) initial d=0.29m and after some wear d=0.30m, and (iii) initial d=0.45m and after some wear d=0.46m. Which pendulum shows the least error in its time period? What is the connection between this result and the plot obtained in (c)? Answer: The error in the time period of the pendulum for the three cases is (i) 1.767%, (ii) 0.034%, and (iii) 0.470% respectively.

Filename:pfigure-dynbalancers
Figure 16.138:

16.4.35  A uniform stick of length and mass m is a hair away from vertically up position when it is released with no angular velocity (a ‘hair’ is a technical word that means ‘very small amount, zero for some purposes’). It falls to the right. What is the force on the stick at point O when the stick is horizontal. Solve in terms of , m, g, ıˆ, and ȷˆ. Carefully define any coordinates, base vectors, or angles that you use. Answer: R=m2θ˙2ıˆ+7mg4ȷˆ.

Filename:pfigure-s94h7p5
Figure 16.139:

16.4.36   A massless 10 meter long bar is supported by a frictionless hinge at one end and has a 3.759kg point mass at the other end. It is released at t=0 from a tip angle of ϕ=.02radians measured from vertically upright position (hinge at the bottom). Use g=10ms2.

  1. (a)

    Using a small angle approximation and the solution to the resulting linear differential equation, find the angle of tip at t=1s and t=7s. Use a calculator, not a numerical integrator. Answer: Angle of tip is 0.031rad at t=1s and 10.967rad at t=7s.

  2. (b)

    Using numerical integration of the non-linear differential equation for an inverted pendulum find ϕ at t=1s and t=7s.

  3. (c)

    Make a plot of the angle versus time for your numerical solution. Include on the same plot the angle versus time from the approximate linear solution from part (a).

  4. (d)

    Comment on the similarities and differences in your plots.

16.4.37  A zero length spring (relaxed length 0=0) with stiffness k=5N/m supports the pendulum shown.

  1. (a)

    Find θ¨ assuming θ˙=2rad/s, θ=π/2. Answer: θ¨=0rad/s2.

  2. (b)

    Find θ¨ as a function of θ˙ and θ (and k, , m, and g.) Answer: θ¨=sinθm(Dkmg).

[Hint: use vectors (otherwise it’s hard)]

[Hint: For the special case, kD=mg, the solution simplifies greatly.]

Filename:pfigure-blue-81-1
Figure 16.140:

16.4.38  A spring-mass-damper system is depicted in the figure. The horizontal damping force applied at B is given by FD=cy˙B
The dimensions and parameters are as follows:

rB/0 = 2ft
rA/0 = =3ft
k = 2lbf/ft
c = 0.3lbfs/ft

For small θ, assume that sin(θ)θ and cos(θ)1.

  1. (a)

    Determine the natural circular frequency of small oscillations about equilibrium for the pendulum shown. The static equilibrium position is θ=0 (pendulum hanging vertically), so the spring is at its rest point in this position.

    Idealize the pendulum as a point mass attached to a rigid massless rod of length 1, so IO=m2. Also use the “small angle approximation” where appropriate. Answer: ωn=kB2+AmgmA2.

  2. (b)

    Sketch a graph of θ as a function of t(t0) if the pendulum is released from rest at position θ=0.2rad when t=0. Your graph should show the correct qualitative behavior, but calculations are not necessary.

Filename:pfigure-blue-137-2
Figure 16.141:

16.4.39  The asymmetric dumbbell shown in the figure is pivoted in the center and also attached to a spring at one quarter of its length from the bigger mass. When the bar is horizontal, the compression in the spring is ys. At the instant of interest, the bar is at an angle θ from the horizontal; θ is small enough so that yL2θ. If, at this position, the velocity of mass ‘m’ is vȷˆ and that of mass 3m is vȷˆ, evaluate the power term (𝑭𝒗) in the energy balance equation. Answer: P=2mgkv2(Lθ4+ys).

Filename:pfigure-blue-78-2
Figure 16.142:

16.4.40  The dumbbell shown in the figure has a torsional spring with spring constant k (torsional stiffness units are kNmrad). The dumbbell oscillates about the horizontal position with small amplitude θ. At an instant when the angular velocity of the bar is θ˙𝒌ˆ, the velocity of the left mass is Lθ˙ȷˆ and that of the right mass is Lθ˙ȷˆ. Find the expression for the power P of the spring on the dumbbell at the instant of interest. Answer: P=mgLθ˙kθθ˙.

Filename:pfigure-blue-130-2
Figure 16.143:

16.4.41  A square plate with side and mass m is hinged at one corner in a gravitational field g. Find the period of small oscillation. Answer: T=2π223g.

Filename:pfigure-s94h10p4
Figure 16.144:

16.4.42  A thin hoop of radius R and mass M is hung from a point on its edge and swings in its plane. Assuming it swings near to the position where its center of mass G is below the hinge:

  1. (a)

    What is the period of its swinging oscillations? Answer: T=2π2Rg.

  2. (b)

    If, instead, the hoop was set to swinging in and out of the plane would the period of oscillations be greater or less? Answer: The period will reduce to T=2π3R2g.

  3. (c)

    What length of simple pendulum (a massless rod with a point mass at the end) has the same period of oscillation of this hoop, swinging in its plane? Answer: =2R

Filename:pfigure-blue-81-2
Figure 16.145:

16.4.43  Oscillating disk. A uniform disk with mass m and radius R pivots around a frictionless hinge at its center. It is attached to a massless spring which is horizontal and relaxed when the attachment point is directly above the center of the disk. Assume small rotations and the consequent geometrical simplifications. Assume the spring can carry compression. What is the period of oscillation of the disk if it is disturbed from its equilibrium configuration? [You may use the fact that, for the disk shown, 𝑯˙/O=12mR2θ¨𝒌ˆ, where θ is the angle of rotation of the disk.]. Answer: period = π2mk.

Filename:pg147-2
Figure 16.146:

16.4.44

Robotics problem: Simplest balancing of an inverted pendulum. You are holding a stick upside down, one end is in your hand, the other end sticking up. To simplify things, think of the stick as massless but with a point mass at the upper end. Also, imagine that it is only a two-dimensional problem (either you can ignore one direction of falling for simplicity or imagine wire guides that keep the stick from moving in and out of the plane of the paper on which you draw the problem).

You note that if you model your holding the stick as just having a stationary hinge then you get ϕ¨=glsinϕ. Assuming small angles, this hinge leads to exponentially growing solutions. Upside-down sticks fall over. How can you prevent this falling?

One way to do keep the stick from falling over is to firmly grab it with your hand, and if the stick tips, apply a torque in order to right it. This corrective torque is (roughly) how your ankles keep you balanced when you stand upright. Your task in this assignment is to design a robot that keeps an inverted pendulum balanced by applying appropriate torque.

Your model is: Inverted pendulum, length , point mass m, and a hinge at the bottom with a motor that can apply a torque Tm. The stick might be tipped an angle ϕ from the vertical. A horizontal disturbing force F(t) is applied to the mass (representing wind, annoying friends, etc).

  1. (a)

    Draw a picture and a FBD

  2. (b)

    Write the equation for angular momentum balance about the hinge point. Answer: F(t)cosϕmgsinϕ+Tm=m2ϕ¨.

  3. (c)

    Imagine that your robot can sense the angle of tip ϕ and its rate of change ϕ˙ and can apply a torque in response to that sensing. That is you can make Tm any function of ϕ and ϕ˙ that you want. Can you find a function that will make the pendulum stay upright? Make a guess (you will test it below).

  4. (d)

    Test your guess the following way: plug it into the equation of motion from part (b), linearize the equation, assume the disturbing force is zero, and see if the solution of the differential equation has exponentially growing (i.e. unstable) solutions. Go back to (c) if it does and find a control strategy that works.

  5. (e)

    Pick numbers and model your system on a computer using the full non-linear equations. Use initial conditions both close to and far from the upright position and plot ϕ versus time.

  6. (f)

    If you are ambitious, pick a non-zero forcing function F(t) (say a sine wave of some frequency and amplitude) and see how that affects the stability of the solution in your simulations.

Filename:pfigure-s94h9p1
Figure 16.147:

16.4.45  A thin rod of mass m and length is hinged with a torsional spring of stiffness K at A, and is connected to a thin disk of mass M and radius R at B. The spring is uncoiled when θ=0. Determine the natural frequency ωn of the system for small oscillations θ, assuming that the disk is:

  1. (a)

    welded to the rod, and Answer: ωn=gL(M+m2)+K(M+m3)L2+MR22.

  2. (b)

    pinned frictionlessly to the rod. Answer: ωn=gL(M+m2)+K(M+m3)L2. Frequency higher than in (a)

Filename:pfigure-blue-151-4
Figure 16.148:

Mixed Problems

16.4.46  Assume that the pulley shown in figure(a) rotates at a constant speed ω. Let the angle of contact between the belt and pulley surface be θ. Assume that the belt is massless and that the condition of impending slip exists between the pulley and the belt. The free-body diagram of an infinitesimal section ab of the belt is shown in figure(b).

  1. (a)

    Write the equations of linear momentum balance for section ab of the belt in the ıˆ and ȷˆ directions.

  2. (b)

    Eliminate the normal force N from the two equations in part (a) and get a differential equation for the tension T in terms of the coefficient of friction μ and The contact angle θ. Answer: T1T2=Ro+RisinθRoRisinθ, where θ=tan1μ.

  3. (c)

    Show that the solution to the equation in part (b) satisfies T1T2=eμθ, where T1 and T2 are the tensions in the lower and the upper segments of the belt, respectively.

Filename:pfigure4-rpf
Figure 16.149:

16.4.47  A belt drive is required to transmit 15 kW power from a 750mm diameter pulley rotating at a constant 300rpm to a 500mm diameter pulley. The centers of the pulleys are located 2.5m apart. The coefficient of friction between the belt and pulleys is μ=0.2.

  1. (a)

    (See problem 16.4.) Draw a neat diagram of the pulleys and the belt-drive system and find the angle of lap, the contact angle θ, of the belt on the driver pulley.

  2. (b)

    Find the rotational speed of the driven pulley. Answer: N2=450rpm.

  3. (c)

    (See the figure in problem 16.4.) The power transmitted by the belt is given by power = net tension × belt speed, i.e., P=(T1T2)v, where v is the linear speed of the belt. Find the maximum tension in the belt. [Hint: T1T2=eμθ (see problem 16.4).] Answer: T1max=2793.8N.

  4. (d)

    The belt in use has a 15mm×5mm rectangular cross-section. Find the maximum tensile stress in the belt. Answer: σmax=37.25 MPa.

16.4.48  Slippery money A round uniform flat horizontal platform with radius R and mass m is mounted on frictionless bearings with a vertical axis at 0. At the instant of interest it is rotating counter clockwise (looking down) with angular velocity 𝝎=ω𝒌ˆ. A force in the xy plane with magnitude F is applied at the perimeter at an angle of 30 from the radial direction. The force is applied at a location that is ϕ from the fixed positive x axis. At the instant of interest a small coin sits on a radial line that is an angle θ from the fixed positive x axis (with mass much much smaller than m). Gravity presses it down, the platform holds it up, and friction (coefficient=μ) keeps it from sliding.

Find the biggest value of d for which the coin does not slide in terms of some or all of F,m,g,R,ω,θ,ϕ, and μ. Answer: dmax=μmgRF2+(mRω)2.

Filename:pfigure-slipperymoney
Figure 16.150:

16.4.49  Frequently parents will build a tower of blocks for their children. Just as frequently, kids knock them down. In falling (even when they start to topple aligned), these towers invariably break in two (or more) pieces at some point along their length. Why does this breaking occur? What condition is satisfied at the point of the break? Will the stack bend towards or away from the floor after the break? Answer: The stack will bend away from the floor after the break. This breaking occurs because gravity is not able to provide equal angular momentum to all the falling blocks.

Filename:pfigure-blue-129-1
Figure 16.151: