Chapter 17 General planar motion of a single rigid object

Download this chapter by section

All chapters and individual parts

The main goal here is to generate equations of motion for general planar motion of a (planar) rigid object that may roll, slide or be in free flight. Multi-object systems are also considered so long as they do not involve other kinematic constraints between the bodies. Features of the solution that can be obtained from analysis are discussed, as are numerical solutions.

Many machine and structural parts move in straight-lines (Chapter 14) or circles (Chapters 13). But other things have more general motions, like a plane in unsteady flight or a connecting rod in a car engine. Keeping track of such motion is a bit more difficult.

In this chapter we will use these two modeling approximations:

  • The objects are planar, or symmetric with respect to a plane; and

  • They have planar motions in that plane.

Filename:tfigure8-rel-ang-vel
Figure 17.1: Planar motion of a 3D car. If the car is symmetrical it can be studied by the means of this chapter.

A planar object is one where the whole object is flat and all its matter is confined to one plane, say the xy plane. This is a palatable approximation for a piece cut out of flat sheet metal. For more substantial real objects, like a full car, the approximation seems at a glance to be terrible. But it turns out that so long as the motion is planar and the car is reasonably idealized as symmetrical (left to right) that treating the car as equivalent to it being squished into a plane does not introduce any more approximation. Thus, even in this 3-D world we live in with 3-D objects, it is fruitful to do 2-D analysis of the type you will learn in this chapter.

A planar motion is one where the velocities of all points are in the same constant plane, say a fixed xy plane, at all times and where points with, say, the same z coordinate have the same velocity. The positions of the points do not have to be in the same plane for a planar motion. Each point stays in a plane, but different points can be in different planes, with each plane parallel to the others.

Example: A car going over a hill

Assume the road is straight in map view in, say, the x direction. Assume the whole width of the road has the same hump. Although the car is clearly not planar, the car motion is probably close to planar, with the velocities of all points in the car in the xy plane (see fig. 17.1)

Filename:tfigure8-ang-vel-ex
Figure 17.2: Planar motion of a skewered sphere. This can be studied by the means in this chapter.

Example: Skewered sphere

A sphere skewered and rotating about a fixed axis in the 𝒌ˆ direction has a planar motion (see fig. 17.2). The points on the object do not all lie in a common plane. But all of the velocities are orthogonal to 𝒌ˆ and thus in the xy plane. This problem does fit in with the methods of this chapter. The symmetry of the sphere with respect to the xy plane makes it so that the three-dimensional mass distribution does not invalidate the two-dimensional analysis.

Filename:tfigure8-ang-vel
Figure 17.3: Planar motion of a planar object. But the plane of the motion is not the plane of the object. This is not a natural topic for this chapter.

Example: Skewered plate

A flat rectangular plate with normal 𝒏ˆ has a fixed axis of rotation in the direction 𝝀ˆ that makes a 45 angle to 𝒏ˆ (see fig. 17.3). This is a planar object (a plane normal to 𝒏ˆ) in planar motion (all velocities are in the plane normal to 𝝀ˆ). But the plane of motion is not the plane of the mass distribution, the object is not symmetric with respect to a motion plane, so this example does not fit into the discussion of this chaptermargin: Actually, a two-dimensional analysis of the plate in this example would be legitimate in this sense. Project all the plate’s mass into the plane normal to the 𝝀ˆ direction. The projections of the forces on this plane would be correctly predicted, but three-dimensional effects, like those associated with dynamic imbalance, would be lost in this projection. .

No real object is exactly planar and no real motion is exactly a planar motion. But many objects are relatively flat and thin or symmetrical and many motions are approximately planar motions. Thus many, if not most, simple engineering analyzes assume planar motion. For bodies that are approximately symmetric about the xy plane of motion (such as a car, if the asymmetrically placed driver’s mass etc.  is neglected), there is no loss in doing a two-dimensional planar rather than full three dimensional analysis.

The plan of this chapter. We start with planar kinematics. Then we evaluate and use expressions for the rates of change of linear and angular momentum for planar bodies. Finally we discuss rolling, sliding and collisions.

17.1 Description of motion: planar rigid-object kinematics

As a rigid object moves, how do the points on it move? There are two reasons to ask this question. First, velocities and accelerations of mass points are needed to apply the momentum-balance equations. Second, formulas for positions, velocities and accelerations of points are useful to understand mechanisms, machines where various parts (each one usually idealized as a rigid object) are connected to each other with hinges and bearings of one type or another.

The central observation in all rigid-object kinematics is the definition of a rigid object:

All pairs of points on a single rigid object keep constant distance from each other as the object moves.

In this section you will learn how to use rigidity to calculate positions, velocities and accelerations of all points (millions and billions of them) on a rigid object given only a few numbers (9 of them). This goal is achieved by putting together the ideas from Chapter 11 (arbitrary motion of one particle), Chapter 14 (straight-line motion), and Chapter chapter 15 (circular motion of a rigid object in a plane).

Displacement and rotation

When a planar object (say a machine part) moves from one configuration in the plane to another it has a displacement and a rotation. It starts in some reference configuration *. We mark a reference point on the body that, in the reference configuration, coincides with a fixed reference point, say 0

margin: The body never has to pass through this reference position, however. For example, the position of an airplane flying from New York to Mumbai is measured relative to a point in the Gulf of Guinea 1000 miles west of Gabon ( the location of O longitude and O latitude) even though the airplane never goes there (nor does anyone want it to).

.

We also mark a (directed) line on the body that, in the reference configuration, coincides with a fixed reference line, say the positive x axis.

Filename:tfigure8-ang-accel
Figure 17.4: Rotation of object is measured by the rotation of real or imagined lines marked on the object. The lines make different angles: θ1θ2, θ2θ3 etc, but θ˙1=θ˙2=θ˙3=. Angular velocity is defined as 𝝎=ω𝒌ˆ with ωθ˙1=θ˙2=θ˙3=.
Filename:tfigure8-ang-accel-ex
Figure 17.5: The displacement and rotation of a planar object relative to a reference configuration.

We could measure rotation by measuring the rotation of any line segment connecting any pair of points fixed to the object. For each line we keep track of the angle that line makes with a line fixed in space, say the positive x or y axis. It’s simplest to stick to the convention that counter-clockwise rotations are positive (fig. 17.4).

The angles θ1, θ2, all change with time and are all different from each other. But all the angles change the same amount, just like in section 16.1. We can pick any one line we like for definiteness and measure the object rotation by the rotation of that line. So

The net motion of a rigid planar object is described by translation, the vector displacement of a reference point from a reference position 𝒓o/o=𝒓oo, and a scalar rotation θ of the reference line from its reference orientation.

The position of a point on a moving rigid object.

Given that P on the object is at 𝒓P/0 in the reference (*) configuration, where is it (What is 𝒓P/0?) after the object has been displaced by 𝒓0/0 and rotated an angle θ?

The base-independent or direct vector notation. An easy way to treat this is to write the new position of P as (see fig. 17.5), in the base-independent or direct vector representation of the position of P,

𝒓P/0=𝒓0/0+𝒓P/0. (17.1)

The vector 𝒓0/0 describes translation, that’s half the story. The other term 𝒓P/0 we find by rotating 𝒓P/0 as we did in Section 16.2. This ‘base-independent’ formula is correct no matter what base vectors are used to represent the vectors in the formula (e.g., ıˆ and ȷˆ or ıˆ and ȷˆ).

Fixed basis or component representation. For a given basis, say ıˆ and ȷˆ associated with x and y, we can write the coordinates of a point as,

[𝒓P/0]xy=[𝒓0/0]xydisplacement+[R(θ)][𝒓P/0]xyrotation (17.2)

or, writing out all the components of the vectors and matrices, We get the fixed basis or component representation of the motion:

[xPyP]=[x0/0y0/0]+[cosθsinθsinθcosθ][xP/0yP/0]. (17.3)

As the motion progresses the displacement [x0/0y0/0] changes with time as does the rotation angle θ. In this fixed basis or component representation of the motion of eqn. (17.3) we get the components of the position in terms of base vectors that are fixed in space.

Example:

If in the reference position a particle on a rigid object is at 𝒓P/0=(1ıˆ+2ȷˆ)m and the object displaces by 𝒓0/0=(3ıˆ+4ȷˆ)m and rotates by θ=π/3rad=60 relative to that configuration, then its new position is:

[𝒓P/0]xy = [𝒓0/0]xy+[R(θ)][𝒓P/0]xy
= [x0/0y0/0]+[cosθsinθsinθcosθ][xP/0yP/0]
= {[34]+[cosπ/3sinπ/3sinπ/3cosπ/3][12]}m
= [3.535+3/2]m
𝒓P/0 = ((3.53)ıˆ+(5+3/2)ȷˆ)m

Changing-base representation. Finally, the changing base representation uses base vectors ıˆ,ȷˆ that are aligned with ıˆ,ȷˆ in the reference configuration but which are glued to the rotating object. If we define x and y as the x and y components of P in the reference (*) configuration we have that

[𝒓P/0]xy = [𝒓P/0]xy=[xy] (17.9)
so𝒓P/0 = (x0/0ıˆ+y0/0ȷˆ)+(xıˆ+yȷˆ). (17.10)

The second equation above is the changing base representation.

Although all three notations are in some sense equivalent, they have their best uses depending on context. The base-independent or direct-vector notation eqn. (17.1) is the most direct, least ornate and shortest to write. The component or fixed base representation eqn. (17.3) is the best for computer calculations. And the changing-base notation eqn. (17.10) is the most explicit for pencil and paper work.

Angular velocity

Filename:tfigure8-ang-accel1
Figure 17.6: It is generally best to take positive ω to be counterclockwise when viewed from the positive z axis.

Because all lines on object rotate at the same rate (at a given instant) ’s rotation rate is the single number we call ω (‘omega b’). As for the case of pure rotation, we define a vector angular velocity with magnitude ω which is perpendicular to the xy plane:

𝝎=ωθ˙𝒌ˆ

where θ˙ is the rate of change of the angle of any line marked on object .

So long as you are careful to define angular velocity by the rotation of line segments and not by the motion of individual particles, the concept of angular velocity in general motion is defined exactly as for an object rotating about a fixed axis. The key kinematic fact is, in words:

On a rigid object in 2D any given point is moving in circles about any other given point (relative to that point).

So you can think of the general planar motion of a rigid object as the general motion of a point plus uniform circular motion about that point.

When a rigid object moves it always has an angular velocity (possibly zero). If we call the object (script B for body), we then call the object’s angular velocity 𝝎. Again, it is generally best to use the CCW sign convention that when ω>0 the object is rotating counterclockwise when viewed looking in from a point outside the object on the positive z axis (that is, ‘looking down’, see fig. 17.6).

The angular velocity vector 𝝎 of an object describes its rate and direction of rotation. For planar motions 𝝎=ω𝒌ˆ.

Relative velocity of two points on a rigid object

For any two points A and B glued to a rigid object the relative velocity 𝒗B/A of the points (‘the velocity of B relative to A’) is given by the cross product of the angular velocity of the object with the relative position of the two points (fig. 17.7):

𝒗B/A𝒗B𝒗A=𝝎×𝒓B/A. (17.11)
Filename:sfig8-2-3
Figure 17.7: The relative velocity of points A and B is in the xy plane and perpendicular to the line segment AB.

That is, the relative velocity of two points on a rigid object is the same as would be predicted for one of the points if the other were stationary. The derivation of this formula is the same as for planar circular motion.

Note that we use three-dimensional cross products even though we are doing planar 2D kinematics. Generally the plane of motion is the xy plane and 𝝎 will be in the z direction. Because 𝝎×𝒓 must be perpendicular to 𝝎 it is perpendicular to the z axis. So the three-dimensional cross product 𝝎×𝒓 gives a vector that is perpendicular to 𝒓 and is in the xy plane .

We can also represent the relative velocity in the changing base notation as

𝒗B/A = ddt(xB/Aıˆ+yB/Aȷˆ)
= xB/Addtıˆ+yB/Addtȷˆ
= xB/A(𝝎×ıˆ)+yB/A(𝝎×ȷˆ).

That is, the relative velocity is the same as the relative position, but with the base vectors substituted with their rates of change.

Finally, we can use the fixed-base or component notation for position and differentiate the terms (see box 16.2 on page 16.2):

[𝒗B/A]xy = ddt[xB/AyB/A]
= ddt{[cosθsinθsinθcosθ][xB/AyB/A]}
= [θ˙sinθθ˙cosθθ˙cosθθ˙sinθ][xB/AyB/A]
= [0ωω0][cosθsinθsinθcosθ][xB/AyB/A]

where xB/A and yB/A are the components of the position of B with respect to A in the reference configuration and hence do not change with time.

Absolute velocity of a point on a rigid object

Filename:sfig8-2-3a
Figure 17.8: Absolute velocity of a point on a rigid object: 𝒗B=𝒗A+ω×𝒓B/A

If one knows the velocity of one point on a rigid object and one also knows the angular velocity of the object, then one can find the velocity of any other point. How? By addition. Say we know the velocity of point A, the angular velocity of the object, and the relative position of A and B, then

𝒗B = 𝒗A+(𝒗B𝒗A) (17.16)
= 𝒗A+𝒗B/A
= 𝒗A+𝝎×𝒓B/A𝒓B𝒓A.

That is, the absolute velocity of the point B is the absolute velocity of the point A plus the velocity of the point B relative to the point A (fig. 17.8).

Equation (17.16) is a 2D vector equation, and thus equivalent to 2 scalar equations, in 4 vectors which have 7 independent components (𝝎 only has one independent component). Generally one could use eqn. (17.16) to solve for any 2 components in terms of the other 5.

Filename:sfig8-2-3b
Figure 17.9: The center of rotation C is that point ‘on’ the object that has zero velocity. So long as ω0 there always is such a point.

Instantaneous Center of Rotation. Interestingly, for any planar object with non-zero ω there is some point C ‘on’ the object that has exactly zero velocity (fig. 17.9). By ‘on’ the object we mean a point that moves with the object in a rigid way. Imagine a giant sheet of clear plexiglass glued to the object and extending way beyond the object. By ‘on’ the object we mean a point on the plexiglass. Why must there be such a point? Consider a point A on the object with some velocity. Relative to A the other points are going in circles. Each radial line extending from A has a relative velocity orthogonal to that line. Thus, considering all radial lines, every velocity direction is represented. And the relative velocities have magnitudes from 0 to infinity as points on a radial line are scanned from A outwards. Thus each possible relative velocity is had by some point ‘on’ the body. And one of those exactly cancels the velocity of A. Thus some point C has no velocity and all other points on the body rotate around it. This point is called the center of rotation. For other than purely circular motion, as an object moves, its center of rotation changes with time. See Sample 17.1 on page 17.1 for more details.

When trying to understand a mechanism it is often useful to think about the locus of locations of the center of rotation of each of the parts.

Angular acceleration

The angular acceleration 𝜶 (‘alpha’) of a rigid object is the rate of change of angular velocity, 𝜶=𝝎˙. The angular acceleration of a object is 𝜶.

As for angular velocity, in 2D angular acceleration is perpendicular to the plane of motion:

𝝎=ω𝒌ˆand𝜶=α𝒌ˆ=ω˙𝒌ˆ.

In 2D the angular acceleration is only due to the speeding up or slowing down of the rotation rate; i.e., α=ω˙=θ¨.

Relative acceleration of two points on a rigid object

For any two points A and B glued to a rigid object , the acceleration of B relative to A is

𝒂B/A = ddt𝒗B/A (17.17)
= ddt{𝝎×𝒓B/A}
= 𝝎˙×𝒓B/A+𝝎×(𝒗B/A),
= 𝝎˙×𝒓B/A+𝝎×(𝝎×𝒓B/A),
= α𝒌ˆ×𝒓B/A+(ω2𝒓B/A),

This is the base-independent or direct-vector expression for relative acceleration. If point A has no acceleration, this formula for relative acceleration reduces to the total acceleration of a point undergoing circular motion from Chapter 15.

Equation (17.17) could also be derived, with some algebra, by taking two time derivatives of the relative position coordinate expression

[𝒓B/A]xy=[R(θ)][𝒓B/A]xy

or by taking two time derivatives of the changing base vector expression

𝒓B/A=xB/Aıˆ+yB/Aȷˆ.

Absolute acceleration of a point on a rigid object

If one knows the acceleration of one point on a rigid body and the angular velocity and acceleration of the body, then one can find the acceleration of any other point. How? Add the absolute acceleration of one point to the acceleration of the second relative to the first.

𝒂B = 𝒂A+(𝒂B𝒂A)=𝒂A+𝒂B/A (17.18)
= 𝒂A+𝝎×(𝝎×𝒓B/A)+𝝎˙×𝒓B/A
= 𝒂Aω2𝒓B/A+α𝒌ˆ×𝒓B/A

This is the base-independent or direct-vector expression for acceleration. The fixed-base (component) and changing-base notations are somewhat more complex.

Equation 17.18 is often called ‘the three term acceleration formula’ because the acceleration of B is the sum of three terms. First is 𝒂A, the acceleration of some point A on the object. Second is 𝝎×(𝝎×𝒓B/A), the centripetal acceleration of B going in circles relative to A. It is directed from B towards A. Third is 𝝎˙×𝒓B/A, due to the change in the magnitude of the angular velocity. The third term is in the direction normal to the line from A to B (tangent to the circle of B going around A).

Example: Robot arm

Filename:sfig8-3-1
Figure 17.10: A two link robot arm.

Given the configuration shown in fig. 17.10 the acceleration of point B can be found by thinking of link AB as the object in eqn. (17.18) and using what you know about circular motion to find the acceleration of A:

𝒂B = 𝒂Aω2𝒓B/A+α𝒌ˆ×𝒓B/A
= (ωOA2ȷˆω˙OAıˆ)(ωAB2ıˆ)+(ω˙AB𝒌ˆ×(ıˆ))
= (ω˙OA+ωAB2)ıˆ+(ωOA2+ω˙AB)ȷˆ

[Note that ωABθ˙ where θ is the angle between the links. Rather ωAB=ω0A+θ˙.]

Computer graphics

Given one point given by the xy pair [x0y0] we can find out what happens to it by rotation [R] as

[xy]=[R][x0y0].

For example the point [02] gets changed by a 45deg rotation to

[xy]=[R][x0y0] = [cosπ4sinπ4sinπ4cosπ4][02]
[.7.7.7.7][02][1.41.4].

A translation is just a vector addition. For example the point [1.41.4] gets translated a distance 2 in the y direction by the addition of [xtyt]=[02] like this

[xy]translated=[xy]+[xtyt]=[1.41.4]+[02]=[1.43.4].

Putting these together the point [x0y0] gets rotated and translated by first multiplying by the rotation matrix and then adding the translation:

[xy]=[R][x0y0]+[xtyt] [.7.7.7.7][02]+[02][1.43.4].

A collection of points all rotated the same amount and then all translated the same amount keep their relative distances.

A picture is a set of points on a plane. If all the points are rotated and translated the same amount the picture is rotated and translated. Thus a picture of a rigid object described by points is rigidly rotated and translated. On a computer line drawings are often represented as a connect-the-dots picture. The picture is represented by the x and y coordinates of the reference dots at the corners. These can be stored in an array with the first row being the x coordinates and the second row the y coordinates as explained on page 16.1. Each column of this matrix represents one point of the connect-the-dots picture. Thus a primitive picture of a house at the origin is given by the array

[P0][xy points originally]=[022100002321]

with the lower left corner of the house at the origin.

To rotate this picture we rotate each of the columns of the matrix P0]. But this is exactly what is accomplished by the matrix multiplication [R][P0]. To translate the points you add the translation vector to each of the columns of the resulting matrix. Thus the whole picture rotated by 45 and translated up by 1 is given by

Filename:sfig8-3-1a
Figure 17.11: (a) A house drawn as 6 dots connected by line segments. The first and last point are the same. (b) The same house but rigidly rotated and translated.
[Pnew]=[R][P0]+[xtyt] [.7.7.7.7][P0]+[01]

which gives a new array of points that, when connected give the picture shown. We have allowed the informal notation of adding a column matrix to a rectangular matrix, by which we mean adding to each column of the rectangular matrix.

To animate the motion of a house flying in, say, the Wizard of Oz, you would first define the house as the set of points [P0]. Then define, maybe by means of numerical solution of differential equations, a set of rotations and translations. Then for each rotation and translation the picture of the house should be drawn, one after the other. The sequence of such pictures is an animation of a flying and spinning house.

Summary of the kinematics of one rigid object in general 2D motion

You can use the position of one reference point and the rotation of the object as simple kinematic measures of the entire motion of the object. If you know the position, velocity, and acceleration of one point on a rigid object (represented by 6 scalars, say) , and you know the angular rate and angular acceleration (3 scalars) then you can find the position, velocity and acceleration of any point on the object (also given its initial position). In 2D, just 9 numbers tell you the position, velocity, and acceleration of any of the billions of points whose initial positions you knowmargin: In 1D it takes just 3 numbers and in 3D just 18. The unusual pattern (3,9,18) comes from rotation being characterized by 0, 1, and 3 numbers in 1, 2, and 3 dimensions, respectively. .

SAMPLE 17.1

Filename:sfig8-5-wiper
Figure 17.12: The rigid triangular plate ABC rotates counterclockwise at 5rad/s. Velocity of point B is given. We have to find the velocity of points A and C.

Velocity of a point on a rigid body in planar motion. An equilateral triangular plate ABC is in motion in the x-y plane. At the instant shown in the figure, point B has velocity 𝒗B=0.3m/sıˆ+0.6m/sȷˆ and the plate has angular velocity 𝝎=2rad/s𝒌ˆ.

  1. 1.

    Find the velocity of point A using vector calculations.

  2. 2.

    Find the velocity of point C using matrix calculations.

Solution

  1. 1.

    We are given 𝒗B and ω, and we need to find 𝒗A, the velocity of point A on the same rigid body. We know that,

    𝒗A=𝒗B+𝝎×𝒓A/B

    Thus, to find 𝒗A, we need to find 𝒓A/B. Let us take an x-y coordinate system whose origin coincides with point A of the plate at the instant of interest and the x-axis is along AB. Then,

    𝒓A/B=𝒓A𝒓B=𝟎(0.2mıˆ)=0.2mıˆ

    Thus,

    𝒗A = 𝒗B+𝝎×𝒓A/B
    = (0.3ıˆ+0.6ȷˆ)m/s+2rad/s𝒌ˆ×(0.2ıˆ)m
    = (0.3ıˆ+0.2ȷˆ)m/s.

    Answer: vA=(0.3ıˆ+0.2ȷˆ)m/s

  2. 2.

    Let xy axes with origin at A be fixed in space with x-axis along AB at the instant of interest. Let xy axes be fixed to the plate with origin at O and the x-axis along AC as shown in fig. 17.13. Thus the coordinates of point C in the rotating axes are

    Filename:sfig8-5-wiper-a
    Figure 17.13: Velocity of point A is known. The axes xy are fixed in space and the axes xy are fixed to the rotating plate in the orientation shown. Thus in the xy axes, the coordinates of point C will always be xC=0.2m,yC=0.
    [xCyC]=[0.2m0]

    Clearly, at the instant of interest, the rotating axes make an angle of θ=π/3 with the fixed axes. Now, we can calculate the velocity of point C as follows.

    [𝒗C]xy=[𝒗A]xy+[0ωω0][cosθsinθsinθcosθ][xC/AyC/A].

    Using θ=π/3, the angle that the rotating axes make with the fixed axes at the instant of interest, and ω=2rad/s, we get

    [𝒗C]xy = [0.3m/s0.2m/s]+[02rad/s2rad/s0][12323212][0.2m0]
    = [0.3m/s0.2m/s]+[0.35m/s0.2m/s]=[0.05m/s0.4m/s].

    Thus, 𝒗C=0.05m/sıˆ+0.4m/sȷˆ. We can also verify this answer by computing 𝒗C=𝒗B+𝝎×𝒓C/B.

    Answer: 𝒗C=0.05m/sıˆ+0.4m/sȷˆ

SAMPLE 17.2

The instantaneous center of rotation. A rigid body is in planar motion. At some instant t, the angular velocity of the body is 𝝎=5rad/s𝒌ˆ and the linear velocity of a point C on the body is 𝒗C=2m/sıˆ5m/sȷˆ. Find the instantaneous center of rotation.


Solution

Let a point O be the instantaneous center of rotation. We need to find the location of O using the given information. Whether O lies inside the object or outside is irrelevant here since the object boundary is not specified. So, let us consider an abstract rigid object as shown in fig. 17.14. What we know about O is that it has zero velocity (since it is the center of rotation).

Filename:sfig8-5-2disks
Figure 17.14:

Let point O be a distance r away from C in some unknown direction 𝒏ˆ (to be found). Thus, 𝒓O/C=r𝒏ˆ is the position vector of point O with respect to point C. Now,

𝒗C = 𝒗O𝟎+𝝎×𝒓C/O𝒓O/C
= ω𝒌ˆ×r𝒏ˆ=ωr(𝒌ˆ×𝒏ˆ).

Thus, 𝒌ˆ×𝒏ˆ must be parallel to and in the opposite direction of 𝒗C. Let 𝒗C=vC𝝀ˆ. Then, we must have 𝒌ˆ×𝒏ˆ=𝝀ˆ, which says that 𝒏ˆ must be normal to 𝝀ˆ. Since 𝝀ˆ=𝒓C/|𝒓C|=(2m/sıˆ5m/sȷˆ)/(22+52m/s)=0.37ıˆ0.93ȷˆ, we can immediately write 𝒏ˆ=0.93ıˆ+0.37ȷˆ. So, now,

vC𝝀ˆ = ωr(𝒌ˆ×𝒏ˆ𝝀ˆ)
 vC = ωr
 r = vC/ω
= 29m/s/(5rad/s)=1.08m.

Thus

𝒓O/C=r𝒏ˆ=1.08m(0.93ıˆ+0.37ȷˆ)=1mıˆ+0.4mȷˆ.

Answer: rO/C=1mıˆ+0.4mȷˆ

Alternatively,

We can find 𝒓O/C purely by vector algebra. Let 𝒓O/C=(xıˆ+yȷˆ)m. So,

𝒗O = 𝒗C+𝝎×𝒓O/C
or 𝟎 = 2m/sıˆ5m/sȷˆ+(5rad/s)𝒌ˆ×(xıˆ+yȷˆ)m
= (25y)m/sıˆ+(5+5x)m/sȷˆ.

Dotting this equation with ıˆ and ȷˆ respectively, we get

25y = 0
5+5x = 0.

Solving these two equations simultaneously, we get, x=1 and y=0.4. Thus,

𝒓O/C=1mıˆ+0.4mȷˆ

as obtained above.

SAMPLE 17.3

Filename:sfig8-4-1
Figure 17.15:

A cheerleader throws her baton up in the air in the vertical xy-plane. At an instant when the baton axis is at θ=60 from the horizontal, the velocity of end A of the baton is 𝒗A=2m/sıˆ+3m/sȷˆ. At the same instant, end B of the baton has velocity in the negative x-direction (but |𝒗B| is not known). If the length of the baton is =12m and the center-of-mass is in the middle of the baton, find the velocity of the center-of-mass.


Solution

We are given:𝒗A = (2ıˆ+3ȷˆ)m/s
and𝒗B = vBıˆ

where vB=|𝒗B| is unknown. We need to find 𝒗G. Using the relative velocity formula for two points on a rigid body, we can write:

𝒗G=𝒗A+𝝎×𝒓G/A (17.25)
Filename:sfig8-4-1a
Figure 17.16:

Here, 𝒗A and 𝒓G/A are known. Thus, to find 𝒗G, we need to find 𝝎, the angular velocity of the baton. Since the motion is in the vertical xy-plane, let 𝝎=ω𝒌ˆ. Then,

Filename:sfig8-6-2
Figure 17.17:
𝒗B = 𝒗A+𝝎×𝒓B/A=𝒗A+ω𝒌ˆ×(cosθıˆsinθȷˆ)𝒓B/A
orvBıˆ = (2ıˆ+3ȷˆ)m/s+ω(cosθȷˆ+sinθıˆ)
= (2ıˆ+3ȷˆ)m/s+ω12m(12ȷˆ+32ıˆ)

Dotting both sides of this equation with ȷˆ we get:

0 = 3m/s+ω2m12
 ω = 3ms41m=43rad/s.

Now substituting the appropriate values in Eqn (17.25) we get:

𝒗G = 𝒗A+ω𝒌ˆ×2(cosθıˆsinθȷˆ)𝒓G/A
= 𝒗A+ω2(cosθȷˆ+sinθıˆ)
= (2ıˆ+3ȷˆ)m/s3m/s(12ȷˆ+32ıˆ)
= (232)m/sıˆ+(332)m/sȷˆ
= 12m/sıˆ+32m/sȷˆ

Answer: 𝒗G=(0.5ıˆ+0.87ȷˆ)m/s

SAMPLE 17.4

Filename:sfig8-6-2a
Figure 17.18: A board in the back of an accelerating truck.

A board in the back of an accelerating truck. A 10ft long uniform board AB rests in the back of a flat-bed truck as shown in Fig. 17.18. End A of the board is hinged to the bed of the truck. The truck is going on a level road at 55 mph. In preparation for overtaking a vehicle in the front, the trucker accelerates at a rate of 3 mph/s. At the instant when the speed of the truck is 60 mph, the magnitude of the relative velocity and relative acceleration of end B with respect to the bed of the truck are 10ft/s and 12ft/s2, respectively. There is wind and at this instant, the board has lost contact with point C. If the angle θ between the board and the bed is 45 at the instant mentioned, find

  1. 1.

    the angular velocity and angular acceleration of the board,

  2. 2.

    the absolute velocity and absolute acceleration of point B, and

  3. 3.

    the acceleration of the center-of-mass G of the board.

Solution At the instant of interest, the quantities given are:

𝒗A = velocity of the truck =60mphıˆ=88ft/sıˆ
𝒂A = acceleration of the truck =3mph/s=4.4ft/s2ıˆ
|𝒗B/A| = vB/A= magnitude of relative velocity of B =10ft/s
|𝒂B/A| = aB/A= magnitude of relative acceleration of B =12ft/s2.

Let 𝝎=ω𝒌ˆ be the angular velocity and 𝝎˙=ω˙𝒌ˆ be the angular acceleration of the board.

  1. 1.

    The relative velocity of end B of the board with respect to end A is

    𝒗B/A = 𝝎×𝒓B/A=ω𝒌ˆ×(cosθıˆ+sinθȷˆ)
    = ω(cosθȷˆsinθıˆ)
     ω = |𝒗B/A|=vB/A=10ft/s10ft=1rad/s.

    Note that we have taken the positive value for ω because the board is rotating counterclockwise at the instant of interest (since it is given that the board has lost contact with point C).

    Similarly, we can compute the angular acceleration:

    𝒂B/A = 𝝎˙×𝒓B/Aω2𝒓B/A
    = ω˙𝒌ˆ×(cosθıˆ+sinθȷˆ)ω2(cosθıˆ+sinθȷˆ)
    = ω˙(cosθȷˆsinθıˆ)ω2(cosθıˆ+sinθȷˆ)
     |𝒂B/A| = (ω˙)2+(ω2)2=aB/A (given)
     aB/A2 = (ω˙)2+(ω2)2
     ω˙ = aB/A22ω4=(12ft/s210ft)2(1rad/s)4
    = ±0.663rad/s2.

    Once again, we select the positive value for ω˙ since we assume that the board accelerates counterclockwise.

    Answer: ω=1rad/skˆ,ω˙=0.663rad/s2kˆ

  2. 2.

    The absolute velocity and the absolute acceleration of the end point B can be found as follows.

    𝒗B = 𝒗A+𝒗B/A
    = vAıˆ+vB/A(cosθȷˆsinθıˆ)
    = 88ft/sıˆ+10ft/s(12ȷˆ12ıˆ)
    = 80.93ft/sıˆ+7.07ft/sȷˆ.
    𝒂B = 𝒂A+𝒂B/A
    = 𝒂A+𝝎˙×𝒓B/Aω2𝒓B/A
    = aAıˆ+ω˙𝒌ˆ×(cosθıˆ+sinθȷˆ)ω2(cosθıˆ+sinθȷˆ)
    = (aAω˙sinθω2cosθ)ıˆ+(ω˙cosθω2sinθ)ȷˆ
    = (4.4ft/s20.66s210ft121s210ft12)ıˆ
    +(0.66s210ft121s210ft12)ȷˆ
    = 7.34ft/s2ıˆ2.40ft/s2ȷˆ.

    Answer: vB=(80.93ıˆ+7.07ȷˆ)ft/s,aB=(7.34ıˆ+2.40ȷˆ)ft/s2.

  3. 3.

    Now, we can easily calculate the acceleration of the center-of-mass as follows.

    𝒂G = 𝒂A+𝒂G/A
    = aAıˆ+𝝎˙×𝒓G/Aω2𝒓G/A
    = aAıˆ+ω˙𝒌ˆ×2(cosθıˆ+sinθȷˆ)ω22(cosθıˆ+sinθȷˆ)
    = aAıˆ+ω˙2(cosθȷˆsinθıˆ)ω22(cosθıˆ+sinθȷˆ)
    = 4.4ft/s2ıˆ+0.663rad/s210ft2(12ȷˆ12ıˆ)
    (1rad/s)210ft2(12ıˆ+12ȷˆ)
    = 1.48ft/s2ıˆ1.19ft/s2ȷˆ.

    Answer: aG=(1.48ıˆ+1.19ȷˆ)ft/s2

Comments: This problem is admittedly artificial. We, however, solve this problem to show kinematic calculations.

SAMPLE 17.5

Filename:efig1-2-28
Figure 17.19:

Tracking motion. A cart moves along a suspended curved path. A rod AB of length =1m hangs from the cart. End A of the rod is attached to a motor on the cart. The other end, B, hangs freely. The motor rotates the rod such that θ(t)=θ0sin(λt) (θ is measured with respect to a fixed vertical line) while the cart moves along the path such that 𝒓A=tıˆ+t318ȷˆ, where all variables (r, t, etc.) are nondimensional.

  1. 1.

    Find the velocity and acceleration of point B as a function of nondimensional time t.

  2. 2.

    Take θ0=π/3 and λ=6. Find and plot the position of the bar at t=0,0.1,0.3,0.9,1,1.1,1.2, and 1.5. Find and draw 𝒗B and 𝒂B at the specified t.


Solution

  1. 1.

    The velocity and acceleration of point B are given by

    𝒗B = 𝒗A+𝒗B/A=𝒗A+𝝎×𝒓B/A
    𝒂B = 𝒂A+𝝎˙×𝒓B/Aω2𝒓B/A.

    Thus, in order to find the velocity and acceleration of point B, we need to find the velocity and acceleration of point A and the angular velocity and angular acceleration of the bar. We are given the position of point A and the angular position of the bar as functions of t. We can, therefore, find 𝒗A, 𝒂A, 𝝎, and 𝝎˙ by differentiating the given functions with respect to t.

    𝒓A = tıˆ+t318ȷˆ
     𝒗A 𝒓˙A=ıˆ+(t2/6)ȷˆ (17.26)
     𝒂A 𝒗˙A=(t/3)ȷˆ (17.27)
      and
    θ𝒌ˆ = θ0sin(λt)𝒌ˆ
     𝝎 θ˙𝒌ˆ=θ0λcos(λt)𝒌ˆ (17.28)
     𝝎˙ θ¨𝒌ˆ=θ0λ2sin(λt)𝒌ˆ. (17.29)

    So,

    𝒗B = 𝒗A+𝝎×(sinθıˆcosθȷˆ) (17.30)
    = ıˆ+(t2/6)ȷˆ+θ˙(sinθȷˆ+cosθıˆ)
    = (1+θ˙cosθ)ıˆ+(t2/6+θ˙sinθ)ȷˆ
    𝒂B = 𝒂A+θ¨𝒌ˆ×(sinθıˆcosθȷˆ)θ˙2(sinθıˆcosθȷˆ)
    = (t/3)ȷˆ+θ¨sinθȷˆ+θ¨cosθıˆθ˙2sinθıˆ+θ˙2cosθȷˆ
    = (θ¨cosθθ˙2sinθ)ıˆ+[t/3+(θ¨sinθ+θ˙2cosθ)]ȷˆ

    where θ=θ0sin(λt),θ˙=θ0λcos(λt), and θ¨=θ0λ2sin(λt)=λ2θ. Thus 𝒗B and 𝒂B are functions of t.

  2. 2.

    The position of the rod at any time t is specified by the position of the two end points A and B (or alternatively, the position of A and the angle of the rod θ). The position of point A is easily determined by substituting the value of t in the given expression for 𝒓A. The position of end B is given by

    𝒓B = 𝒓A+𝒓B/A=tıˆ+(t3/18)ȷˆ+(sinθıˆcosθȷˆ)
    = (t+sinθ)ıˆ+(t3/18cosθ)ȷˆ.

    To compute the positions, velocities, and accelerations of end points A and B at the given instants, we first compute θ, θ˙, and θ¨, and then substitute them in the expressions for 𝒓A,𝒓B,𝒗A,𝒗B,𝒂A, and 𝒂B. A pseudocode for computer calculation is given below.

    t = [0 0.1 0.3 0.9 1.0 1.1 1.2 1.5]
    theta0=pi/3, L=.5, lam=6
    for each t, compute
      theta = theta0*sin(lam*t)
      w = lam*theta0*cos(lam*t)
      wdot = -lam^2*theta
           % Position of A and B
      xA=t, yA=t^3/18
      xB = xA + L*sin(theta)
      yB = yA - L*cos(theta)
           % Velocity of A and B
      uA = 1,  vA = t^2/6
      uB = uA + L*w.*cos(theta)
      vB = vA + L*w.*sin(theta)
           % Acceleration of A and B
      axA = 0,  ayA = t/3
      axB = L*wdot*cos(theta) - L*w^2*sin(theta)
      ayB = ayA + L*wdot*sin(theta) + L*w^2*cos(theta)
    

    From the above calculation, we get the desired quantities at each t. For example, at t=0 we get,

    xA = 0, yA = 0, xB = 0, yB = -0.5
    uA = 1, vA = 0, uB = 4.14, vB = 0, axB = 0, ayB = 19.74
    

    which means,

    𝒓A=𝟎,𝒓B=0.5ȷˆ,𝒗A=ıˆ,𝒗B=4.14ıˆ,𝒂B=19.74ȷˆ.

    The position of the bar, the velocity vectors at points A and B, and the acceleration vector at B, thus obtained, are shown in fig. 17.20 graphically. margin: We can take several values of t, say 400 equally spaced values between t=0 and t=4, and draw the bar at each t to see its motion and the trajectory of its end points. fig. 17.21 shows such a graph.

Filename:pfigure-blue-118-2
Figure 17.20: Position, velocity of the end points A and B, and acceleration of point B at various time instants.
Filename:pfigure-s95f3a
Figure 17.21: Graph of closely spaced configuration of the bar between t=0 to t=4.

Problems for 17.1 Rigid object kinematics

Preparatory Problems

17.1.1  A disk of radius R is hinged at point O at the edge of the disk, approximately as shown. It rotates counterclockwise with angular velocity θ˙=𝝎. A bolt is fixed on the disk at point P at a distance r from the center of the disk. A frame xy is fixed to the disk with its origin at the center C of the disk. The bolt position P makes an angle ϕ with the x-axis. At the instant of interest, the disk has rotated by an angle θ.

  1. (a)

    Write the position vector of point P relative to C in the xy coordinates in terms of given quantities.

  2. (b)

    Write the position vector of point P relative to O in the xy coordinates in terms of given quantities.

  3. (c)

    Write the expressions for the rotation matrix R(θ) and the angular velocity matrix 𝒮(𝝎).

  4. (d)

    Find the velocity of point P relative to C using R(θ) and the angular velocity matrix 𝒮(𝝎).

  5. (e)

    Using R=30cm, r=25cm, θ=60, and ϕ=45, find [𝒓C/0]xy, and [𝒓P/0]xy at the instant shown.

  6. (f)

    Assuming that the angular speed is ω=10rad/s at the instant shown, find [𝒗C/0]xy and [𝒗P/0]xy taking other quantities as specified above.

Filename:Danef94s1q2
Figure 17.22:

17.1.2  A uniform rigid rod AB of length =1m rotates at a constant angular speed ω about an unknown fixed point. At the instant shown, the velocities of the two ends of the rod are 𝒗A=1m/sıˆ and 𝒗B=1m/sȷˆ.

  1. (a)

    Find the angular velocity of the rod.

  2. (b)

    Find the center of rotation of the rod.

Filename:pfigure-blue-123-1
Figure 17.23:

17.1.3  A square plate ABCD rotates at a constant angular speed about an unknown point in its plane. At the instant shown, the velocities of the two corner points A and D are 𝒗A=2ft/s(ıˆ+ȷˆ) and 𝒗D=(2ft/s)ıˆ, respectively.

  1. (a)

    Find the center of rotation of the plate.

  2. (b)

    Find the acceleration of the center of mass of the plate.

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

More-Involved Problems

17.1.4  A ten foot ladder is leaning between a floor and a wall. The top of the ladder is sliding down the wall at one foot per second. (The foot is simultaneously sliding out on the floor). When the ladder makes a 45 degree angle with the vertical what is the speed of the midpoint of the ladder?

17.1.5  The slender rod AB rests against the step of height h, while end “A” is moved along the ground at a constant velocity vo. Find ϕ˙ and ϕ¨ in terms of x, h, and vo. Is ϕ˙ positive or negative? Is ϕ¨ positive or negative?

Filename:pfigure-s95q14
Figure 17.25:

17.1.6  Consider the motion of a rigid ladder which can slide on a wall and on the floor as shown in the figure. The point A on the ladder moves parallel to the wall. The point B moves parallel to the floor. Yet, at a given instant, both have velocities that are consistent with the ladder rotating about some special point, the center of rotation (COR). Define appropriate dimensions for the problem.

  1. (a)

    Find the COR for the ladder when it is at some given position (and moving, of course). [Hint: if a point A is ‘going in circles’ about another point C, that other point C must be in the direction perpendicular to the motion of A.

  2. (b)

    As the ladder moves, the COR changes with time. What is the set of points on the plane that are the COR’s for the ladder as it falls from straight up to lying on the floor?

Filename:Danef94s3q2
Figure 17.26:

17.1.7   A car driver on a very boring highway is carefully monitoring her speed. Over a one hour period, the car travels on a curve with constant radius of curvature, ρ=25mi, and its speed increases uniformly from 50mph to 60mph. What is the acceleration of the car half-way through this one hour period, in path coordinates?

17.1.8  The bar AB shown in the figure is 1 m long. The end A of the bar is dragged to the left on the horizontal floor at a constant speed v0=0.2m/s. The other end of the bar drags along the inclined plane that makes an angle ϕ=60 with the horizontal.

  1. (a)

    Does the end B of the rod also have a constant speed? Does it imply that the rod has a constant angular velocity? Guess your answer and write down your reasoning. No calculations.

  2. (b)

    Find an expression for the angular velocity of the rod in terms of θ, ϕ, , and v0. Answer: ω=v0sinϕcos(ϕθ)kˆ

  3. (c)

    Find the velocity of end B when θ=30.

Filename:bikefork1-ang-accel
Figure 17.27:

17.1.9  The disk shown in the figure drives the slider AB. The disk rotates clockwise at a constant angular speed ωD=60rpm. The radius of the disk is R=0.25m.

  1. (a)

    Find the velocity of collar A when θ=30.

  2. (b)

    Find the acceleration of collar A when θ=30.

  3. (c)

    Find the maximum and minimum velocity of the collar and the corresponding angle θ of the slider.

Filename:bikefork-ang-accel
Figure 17.28:

17.1.10  A rod AC of length =2m is mounted on a cart at point C. The cart slides down under gravity on a frictionless surface inclined at ϕ=45 from the horizontal. Assume that bar has negligible mass compared to that of the cart. The cart is released from rest and it reaches the position shown in t=2s. At the instant shown, θ=60, and rod AC has an angular velocity 𝝎=0.5rad/s𝒌ˆ and angular acceleration 𝜶=1rad/s2𝒌ˆ. Find the velocity and acceleration of the center of mass of the rod at the instant shown.

Filename:bikefork1-alt
Figure 17.29:

17.1.11  A vertical bar pendulum AB of length =4ft is mounted on a cart at point A. The cart has a motor that drives the bar pendulum back and forth about the vertical position such that θ(t)=θ0sinωt where θ0=π/6 and ω=πrad/s. The cart moves at a constant speed v0=2ft/s on a sinusoidal track represented by y(x)=y0sinλx where y0=0.5ft and λ=π/4/ft.

  1. (a)

    Find the velocity of point B when θ=0.

  2. (b)

    Find and plot the position of point B for one full oscillation of the bar, assuming the bar starts from its leftmost angular position (that is, θ=θ0) when the cart is on one of the crests of the track.

Filename:bikefork-alt
Figure 17.30:

17.1.12  The center of mass of a javelin travels on a more or less parabolic path while the javelin rotates during its flight. In a particular throw, the velocity of the center of mass of a javelin is measured to be 𝒗C=10m/sıˆ when the center of mass is at its highest point h=6m. As the javelin lands on the ground, its nose hits the ground at G such that the javelin is almost tangent to the path of the center of mass at G. Neglect the air drag and lift on the javelin.

  1. (a)

    Given that the javelin is at an angle θ=45 at the highest point, find the angular velocity of the javelin. Assume the angular velocity is constant during the flight and that the javelin makes less than a full revolution.

Filename:pfigure-blue-128-1
Figure 17.31:

17.2 General planar mechanics of a rigid-object

We now apply the kinematics ideas of the last section to the general mechanics principles in Table I in the inside cover. The goal is to understand the relation between forces and motion for a planar object in general 2-D motion. The simple measuresmargin: Advanced aside: What we call “simple measures” are examples of “generalized coordinates” in more advanced books. The idea sounds intimidating, but is simply this: If something can only move in a few ways, you should only keep track of the motion with that many variables. The kinematics of a rigid object (Sect. 17.1) allow us to “evaluate” the motion quantities, namely linear momentum, angular momentum, kinetic energy, and their rates of change in terms of these “simple measures”. By “evaluate” we mean express the motion quantities in terms of these measures. The alternative is as sums over Avogadro’s number of particles (There are on the order of 1023 atoms in a typical engineering part.). Even neglecting atoms and viewing matter as continuous we would still be stuck with integrals over complicated regions if we did not describe the motion with as few variables as possible. In the case of 2-D rigid object motion, the position of a reference point (x and y) with the rotation θ is called a set of minimal coordinates. These, and their time derivatives are the minimal information needed to describe all important mechanics motion quantities. of motion will be the displacement of one reference point 0 on the object and the rotation of the object. We also need the first and second time derivatives of the displacement and rotation, altogether we use 𝒓0,𝒗0, 𝒂0,θ,𝝎 and 𝜶.

We will treat all bodies as if they are squished into the plane (planar) and moving in the plane (planar motion). But the analysis is sensible for an object that is symmetric with respect to the plane containing the velocities (see Box 17.2 on page 17.2).

The balance laws for a rigid object

As always, once you have defined the system and the forces acting on it by drawing a free-body diagram, the basic momentum balance equations are applicable (and exact for engineering purposes). Namely,

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

The same point 0, any point, is used on both sides of the angular momentum balance equation.

We also have power balance which, for a system with no internal energy or dissipation, is

Power balance:P=E˙K.

The left hand sides of the momentum balance equations are evaluated the same way, whether the system is composed of one object or many, whether the bodies are deformable or not, and whether the points move in straight lines, circles, hither and thither, or not at all. It is the right hand sides of the momentum equations that involve the motion.

Similarly, in the energy balance equations the applied power P only depends on the position of the forces and the motions of the material points at those positions. But the kinetic energy EK and its rate of change depend on the motion of the whole system.

You already know how to evaluate the momenta and energy, and their rates of change, for a variety of special cases, namely

  • Systems composed of particles where all the positions and accelerations are known (chapter 13);

  • Systems where all points have the same acceleration. That is, the system moves like a rigid object that does not rotate (chapter 14); and

  • Systems where all points move in circles about the same fixed axis in 2D (chapter 15).

Now we go on to consider the general 2-D motions of a planar rigid object. It’s now a little harder to evaluate 𝑳,𝑳˙,𝑯/O,𝑯˙/O,EK and E˙K. But not much.

Summary of the motion quantities

Table I in the back of the book describes the motion quantities for various special cases, including the planar motions we consider in this chapter. Most relevant is row (7).

The basic idea is that the momenta for general motion, which involves translation and rotation, are the sum of the momenta (both linear and angular, and their rates of change too) from those two effects.

Namely, the linear momentum is described, as for any system with any motion, by the motion of the center-of-mass

𝑳=mtot𝒗cm and 𝑳˙=mtot𝒂cm, (17.32)

and the angular momentum has two contributions, one from the motion of the center-of-mass and one from rotation of the object about the center of mass,

𝑯/O = 𝒓cm/O×(mtot𝒗cm)Angular momentum due to motion of the center-of-mass+Izzcm𝝎Angular momentum due to motion relative to the center-of-mass (17.33)
and 𝑯˙/O = 𝒓cm/O×(mtot𝒂cm)+Izzcm𝝎˙. (17.34)

These simplifications, for 2D objects moving in the plane, are discussed in box 17.2 on page  17.2.

An important special case for the angular momentum evaluation is when the reference point is coincident with the center-of-mass. Then the angular momentum and its rate of change simplify to

𝑯cm=Izzcm𝝎 and 𝑯˙cm=Izzcm𝝎˙. (17.35)

The kinetic energy and its rate of change are given by

EK = 12mtotvcm2𝒗cm𝒗cmkinetic energy from center-of-mass motion +12Izzcmω2kinetic energy relative to the center-of-mass (17.36)
and E˙K = mtotvcmv˙cm𝒗cm𝒂cm+Izzcmωω˙ (17.37)

The relations above are easily derived from the general center of mass theorems (see box 17.2 on page 17.2 for some of these derivations).

Equations of motion

Putting together the general balance equations and the expressions for the motion quantities we can now write linear momentum balance, angular momentum balance and power balance as:

LMB:𝑭i=mtot𝒂cm,(a)AMB:𝑴O=𝒓cm/O×(mtot𝒂cm)+Izzcm𝝎˙(b)or𝑴cm=Izzcm𝝎˙,andPower:𝑭tot𝒗cm+𝝎𝑴cm=mtotvv˙+Izzcmωω˙.(c) (17.38)

Independent equations?

Equations are only independent if no one of them can be derived from the others. When counting equations and unknowns one needs to make sure one is writing independent equations. How many independent equations are in the set eqns. (17.38)abc applied to one free object diagram? The short answer is 3.

The linear momentum balance equation 17.38a yields two independent equations by dotting with any two non-parallel vectors (say, ıˆ and ȷˆ). Dotting with a third vector yields a dependent equation.

For any one reference point the angular momentum equation 17.38b yields one scalar equation. It is a vector equation but always has zero components in the ıˆ and ȷˆ directions. But angular momentum equation can yield up to three independent equations by being applied to any set of three non-collinear points.

The power balance equation is one scalar equation.

In total, however, the full set of equations above only makes up a set of three independent equations.

To avoid thinking about what is or is not an independent set of equations some people prefer to stick with one of the canonical sets of independent equations:

  • The coordinate based set (“old standard”)

    • {LMB}ıˆ or, equivalently, Fx=mtotacmx,

    • {LMB}ȷˆ or, equivalently, Fy=mtotacmy, and

    • {AMB}𝒌ˆ or, equivalently, Mcm=Izzcmω˙.

  • Moment only (good for eliminating unknown reaction forces)

    • {AMB about pt A}𝒌ˆ   (A is any point, on or off the object)

    • {AMB about pt B}𝒌ˆ   (B is any other point)

    • {AMB about pt C}𝒌ˆ   (C is a third point not on the line AB)

  • Two moments and a force component

    • {AMB about pt A}𝒌ˆ   (A is any point, on or off the object)

    • {AMB about pt B}𝒌ˆ   (B is any other point)

    • {LMB}𝝀ˆ   (where 𝝀ˆ is not perpendicular to the line AB)

  • Two force components and a moment (also good for eliminating unknown forces)

    • {LMB}𝝀ˆ1   (where 𝝀ˆ1 is any unit vector)

    • {LMB}𝝀ˆ2   (where 𝝀ˆ2 is any other unit vector)

    • {AMB about pt A}𝒌ˆ   (A is any point, on or off the object)

Any of these will always do the job. The power balance equation is often used as a consistency check rather than an independent equation.

From a theoretical point of view one might ask the related question of which of the equations of motion can be derived from the others. There are many answers. Here are some of them:

  • Power balance follows from LMB and AMB,

  • AMB about three non-collinear points implies LMB, and

  • LMB and power balance yield AMB

Interestingly, there is no way to derive angular momentum balance from linear momentum balance without the questionable microscopic assumptions discussed in box LABEL:threeten on page LABEL:threeten.

Some simple examples

Here we consider some simple examples of unconstrained motion of a rigid object.

Example: The simplest case: no force and no moment.

If the net force and moment applied to an object are zero we have:

LMB𝟎 = mtot𝒂cmand
AMB𝟎 = Izzcmω˙𝒌ˆ

so 𝒂cm=𝟎 and ω˙=0 and the object moves at constant speed in a constant direction with a constant rate of rotation, all determined by the initial conditions. Throw an object in space and its center-of-mass goes in a straight line and it spins at constant rate (subject to the 2-D restrictions of this chapter).

Filename:tfigure8-syst-bods
Figure 17.32: The X marked at the center-of-mass of a thrown spinning clipboard follows a parabolic trajectory.

Example: Constant force applied to the center-of-mass.

In this case angular momentum balance about the center-of-mass again yields that the rotation rate is constant. Linear momentum balance is now the same as for a particle at the center-of-mass, i.e., the center-of-mass has a parabolic trajectory.

Near-earth (constant) gravity provides a simple example. An ‘X’ marked at the center-of-mass of a clipboard tossed across a room follows a parabolic trajectory (see fig. 17.32).

Example: Constant force not at the center-of-mass.

Assume the only force applied to an object is a constant force 𝑭=Fıˆ at A (see fig. 17.33). Then linear momentum balance gives us that

𝑭i=𝑳˙Fıˆ=m𝒂G𝒂G=F/mıˆ=constant.

So if the object starts at rest, the point G will move in a straight line in the ıˆ direction (The common intuition that point G will be pulled up is incorrect). Angular momentum about the center-of-mass gives

𝑴cmi=𝑯˙cm {𝒓A/G×Fıˆ=Izzcmθ¨𝒌ˆ}
{}𝒌ˆ θ¨+FIzzcmsinθ=0,

with =|𝒓A/G|, which is the pendulum equation. That is, the object can swing back and forth about θ=0 just like a pendulum, approximately sinusoidally if the angle θ starts small and with θ˙ initially also small. [One might wonder how to do this experiment. One way would be with a jet on a space craft that keeps re-orienting itself to keep in a constant spatial direction as the object changes orientation. Another would be with a string tied to A and pulled from a great distance.]

Filename:sfig8-7-2
Figure 17.33: The only force applied to the object is the constant force 𝑭=Fıˆ applied at point A. The resulting motion is a constant acceleration of the center-of-mass 𝒂G=(F/m)ıˆ and an oscillatory motion of θ identical to that of a pendulum hinged at G. Note that the center of mass goes in a straight line.
Filename:sfig8-7-2a
Figure 17.34: A rocket is pointed in the direction 𝝀ˆ which makes an angle θ with the positive x axis. The position and velocity of the center-of-mass at G are called 𝒓 and 𝒗. The velocity of the tail is 𝒗C

Example: The flight of an arrow or rocket.

As a primitive model of an arrow or rocket assume that the only force is from drag on the fins at C and that this force opposes motion according to

𝑭=c𝒗C

where c is a drag coefficient (see fig. 17.34). From linear momentum balance we have:

𝑭i=𝑳˙  𝑭 = m𝒂
c𝒗C = m𝒗˙
m𝒗˙ = c(𝒗+𝝎×𝒓C/G)
= c(𝒗+θ˙𝒌ˆ×(𝝀ˆ))
(𝒌ˆ×𝝀ˆ=𝒏ˆ)  𝒗˙ = cm(θ˙𝒏ˆ𝒗).

So if 𝒗,θ and θ˙ are known the acceleration 𝒗˙ is calculated by the formula above.

Similarly angular momentum balance about G gives

𝑴G=𝑯˙G  {𝒓C/G×𝑭 = Izzcmω˙𝒌ˆ}
{}𝒌ˆ  Izzcmω˙ = 𝒓C/G×𝑭𝒌ˆ.

Then, making the same substitutions as before for 𝒓C/G and 𝑭 we get

ω˙=cIzzcm(𝝀ˆ×𝒗𝒌ˆθ˙)

which determines the rate of change of ω if the present values of 𝒗,θ and θ˙ are known.

Setting up differential equations for solution

If one knows the forces and torques on an object in terms of its position, velocity, orientation and angular velocity one then has a ‘closed’ set of differential equations. That is, one has enough information to define the equations for a mathematician or a computer to solve them.

The full set of differential equations is gathered from linear and angular momentum balance and also from simple kinematics. Namely, one has the following set of 6 first order differential equations:

x˙=vx,v˙x=Fx/m,y˙=vy,v˙y=Fy/m,θ˙=ω, andω˙=Mcm/Izzcm,

where the positions and velocities are the positions and velocities of the center-of-mass. The expressions for Fx,Fy, and Mcm may well be complicated, as in the rocket example above.

Box 17.1 2-D mechanics makes sense in a 3-D world

The math for two-dimensional mechanics analysis is simpler than the math for three-dimensional analysis. And thus easier to learn first. But we do actually live in a three-dimensional world. You might wonder at the utility of learning something that is not right. There are three answers.

  1. 1.

    Two-dimensional analysis can give partial information about the three-dimensional system that is exactly the same as the three-dimensional analysis would give by projection, no matter what the motion;

  2. 2.

    if the motion is planar the 2-D kinematics can be used; and

  3. 3.

    if the object is planar or symmetric about the motion plane, and any constraints that hold the object are also symmetric about the motion plane, the 2-D analysis is not only correct, but complete.

Of course no machine is exactly planar or exactly symmetric, but if the approximation seems reasonable most engineers will accept a small loss in accuracy for great gain in simplicity.

a) Projection

First let’s relax our assumption of 2-D motion. Consider arbitrary 3-D motions of an arbitrarily complex system. If we take the dot product of the linear momentum equations with ıˆ and ȷˆ and the angular momentum balance equation with 𝒌ˆ we get

{𝑭i=mi𝒂i}ıˆFix=miaix,(a){𝑭i=mi𝒂i}ȷˆFiy=miaiy,and(b){𝒓i×𝑭i=𝒓i×mi𝒂i}𝒌ˆrixFiyriyFix=mi(rixaiyriyaix).(c) (17.39)

These are exactly the equations of 2-D mechanics. That is, if we only consider the planar components of the forces, the planar components of the positions, and the planar components of the motions, we get a correct but partial set of the 3-D equations. In this sense 2-D analysis is correct but incomplete.

b) Planar motion

If all the velocities of the parts of a 3-D system have no z component the motion is planar (in the xy plane). Thus the right-hand sides of eqns. (17.39) are not just projections, but the whole story. Further, in the case of rigid-object motion, the 2-D kinematics equation

𝒗P=𝒗G+ω𝒌ˆ×𝒓P/G=𝒗G+ω𝒌ˆ×(rP/Gxıˆ+rP/Gyȷˆ) (17.40)

also applies (the z component of the position drops out of the cross product) and the expression for, say, the z component of the angular momentum of an object about its center-of-mass is

Hcmz=Izzcmω.

Differentiating, or adding up the mi𝒂i terms we get,

H˙cmz=Izzcmω˙.

Similarly, the z component of the full angular momentum balance equation for a 3-D rigid object in planar motion is the same as the z component of eqn. (17.38)b.

𝑴O𝒌ˆ=(𝒓cm/O×(mtot𝒂cm))𝒌ˆ+Izzcmω˙

So for planar motion of 3-D rigid bodies one can do a correct 2-D analysis with the full ease of analyzing a planar object.

But this result is deceptively simple. The free object diagram in 3-D most likely shows forces in the z direction, pairs of forces in the x or y directions that are applied at points with the same x and y coordinates but different z values, or moments with components in the x or y directions. Full information about these force and moment components can’t be found from 2-D analysis. That is,

the nature of the forces that it takes to keep a system in planar motion can’t be found from a planar analysis.

For example, a system rotating about the z axis which is statically balanced but is dynamically imbalanced

has no net x or y reaction force, as a planar analysis would reveal, yet the bearing reaction forces are not zero.

Another example would be a plan view of a car in a turn (assuming a stiff suspension). A 2-D analysis could be accurate, but would not be complete enough to describe the tire reaction forces needed to keep the car flat.

c) Symmetric bodies and planar bodies

If the rigid object has all its mass in the xy plane, or its mass is symmetrically distributed about the xy plane, and it is in planar motion in the xy plane then

miaiz=0 and {𝒓i×mi𝒂i}ıˆ=0 and {𝒓i×mi𝒂i}ȷˆ=0

where 𝒓 is measured relative to any point in the plane. Thus, by linear and angular momentum balance,

Fz=0 and {𝒓i×𝑭i}ıˆ=0 and {𝒓i×𝑭i}ȷˆ=0

so

A planar object or a symmetric object in planar motion needs no force in the z direction and no moment in the x or y direction to keep it in the plane.

Systems that are symmetric or flat and moving in an approximately planar manner, are thus both accurately and completely modeled with a 2-D analysis. A slight generalization of the result is to any object or collection of objects whose centers of mass are on the plane and each of which is dynamically balanced for rotation about a 𝒌ˆ axis through its center-of-mass.

Box 17.2 The center-of-mass theorems for 2-D rigid bodies

That all the particles in a system are part of one planar object in planar motion (in that plane) allows highly useful simplification of the expressions for the motion quantities, namely Eqns. 17.32 to 17.36. We can derive these expressions from the center-of-mass theorems from chapter LABEL:sec:KonigSchmonig. For completeness, we repeat some of those derivations here. To save space, we only use the integral () forms for the general expressions; the derivations with sums () are similar. In all cases position, velocity, and acceleration are relative to a fixed point in space (that is 𝒓,𝒗, and 𝒂 mean 𝒓/0,𝒗/0, and 𝒂/0 respectively).

Linear momentum.

Here we show that to evaluate linear momentum and its rate of change you only need to know the motion of the center of mass.

𝑳 𝒗𝑑m=ddt𝒓𝑑m=ddt𝒓𝑑m=ddt(mtot𝒓cm)
= mtotddt𝒓cm=mtot𝒗cm

By identical reasoning, or by differentiating the expression above with respect to time,

𝑳˙=mtot𝒂cm

Thus for linear momentum balance one need not pay heed to rotation, only the center-of-mass motion matters.

Angular momentum.

Here we attempt a derivation like the one above but get slightly more complicated results. For simplicity we evaluate angular momentum and its rate of change relative to the origin, but a very similar derivation would hold relative to any fixed point C.

𝑯/O 𝒓×𝒗𝑑m
= (𝒓𝒓cm+𝒓cm)×(𝒗𝒗cm+𝒗cm)𝑑m
= (𝒓/cm+𝒓cm)×(𝒗/cm+𝒗cm)𝑑m
= 𝒓/cm×𝒗/cm𝑑m+𝒓cm×𝒗cm𝑑m
+𝒓/cm×𝒗cm𝑑m+𝒓cm×𝒗/cm𝑑m
= 𝒓/cm×𝒗/cm𝑑m+𝒓cm×𝒗cm𝑑m
+(𝒓/cm𝑑m)𝟎×𝒗cm+𝒓cm×(𝒗/cm𝑑m)𝟎
= 𝒓/cm×𝒗/cm𝑑m+𝒓cm×𝒗cmmtot.

This much is true for any system in any motion. For a rigid object we know about the motions of the parts. Using the center-of-mass as a reference point we know that for all points on the object 𝒗/cm=𝝎×𝒓/cm. Thus we can continue the derivation above, following the same reasoning as was used in chapter 7 for circular motion of rigid bodies:

𝑯/O = 𝒓/cm×(𝝎×𝒓/cm)𝑑m+𝒓cm×𝒗cmmtot.

Using the identity for the triple cross product (see box 17.2 on page 17.2) or using the geometry of the cross product directly with 𝝎=ω𝒌ˆ as in chapters 7 and 8 we get

𝑯/O = ω𝒌ˆr/cm2𝑑m+𝒓cm×𝒗cmmtot.

Then defining Izzcmr/cm2𝑑m we get the desired result:

𝑯/O = 𝒓cm×𝒗cmmtot+Izzcmω𝒌ˆ.

A similar derivation, or differentiation of the result above (and using that (ddt𝒓)×𝒗=𝒗×𝒗=𝟎) gives

𝑯˙/O = 𝒓cm×𝒂cmmtot+Izzcmω˙𝒌ˆ.

The results above hold for any reference point, not just the origin of the fixed coordinate system. Thus, relative to a point instantaneously coinciding with the center-of-mass

𝑯cm = 𝒓cm/cm𝟎×𝒗cmmtot+Izzcmω𝒌ˆ=Izzcmω𝒌ˆ.

and similarly

𝑯˙cm = Izzcmω˙𝒌ˆ.

Kinetic energy.

Unsurprisingly the expression for kinetic energy and its rate of change are also simplified by derivations very similar to those above. Skipping the details (or leaving them as an exercise for the peppy reader):

EK 12𝒗𝒗𝑑m
= 12mtotvcm2+12Izzcmω2

and

E˙K ddtEK
= mtotvv˙+Izzcmωω˙.

Box 17.3 The work of a moving force and of a couple

The work of a force acting on an object from state one to state two is

W12=t1t2P𝑑t.

But sometimes we like to think not of the time integral of the power, but of the path integral of the moving force. So we rearrange this integral as follows.

W12 = t1t2P𝑑t (17.41)
= t1t2𝑭𝒗dtd𝒓
= 𝒓1𝒓2𝑭𝑑𝒓

The validity of equation 17.41 depends on the force acting on the same material point of the moving object as it moves from position 1 to position 2; i.e., the force moves with the object. If the material point of force application changes with time, eqn. (17.41) is senseless and should be replaced with the following more generally applicable equation:

W12t1t2P𝑑t=t1t2𝑭𝒗𝑑t (17.42)

where 𝒗 is the velocity of the material point at the instantaneous location of the applied force.

Hand drags on a passing train: a subtlety

There is a subtle distinction between eqn. (17.41) and eqn. (17.42). As an example think of standing still and dragging your hand on a passing train. Your hand slows down the train with the force

𝑭handontrain.

It might seem that the work of the hand on the train is zero because your hand doesn’t move; work is force times distance and the distance is zero and eqn. (17.41) superficially evaluates to

𝒓1𝒓2𝑭𝑑𝒓=0.

But we have violated the condition for the validity of eqn. (17.41): the force be applied to a fixed material point as time progresses. Whereas on the train your hand smears a whole line of material points.

Clearly your hand does slow the train, so it must do (negative) work on the train, as eqn. (17.42) correctly shows because

Pforceontrain=𝑭handontrain𝒗train0.

The power of the hand force on the train is the force on the train dotted with the velocity of the train (not with the velocity of your hand). Thus, your hand does negative work on the train. eqn. (17.42) applies to the train and eqn. (17.41) does not.

On the other hand (so to speak) if one looks at the power of the force on the hand we have:

𝑭trainonhand=𝑭handontrain

while the velocity of the hand is zero so

Pforceonhand=𝑭trainonhand𝒗hand=0.

So the train does no work on your hand while your hand does (negative) work on the train. The difference, of course, is mechanical energy lost to heat.

Work of an applied torque

By thinking of an applied torque as really a distribution of forces, the work of an applied torque is the sum of the contributions of the applied forces. If a collection of forces equivalent to a torque is applied to one rigid object the power of these forces turns out to be 𝑴𝝎. At a given angular velocity a bigger torque applies more power. And a given torque applies more power to a faster spinning object.

Here’s a quick derivation for a collection of forces 𝑭i that add to zero acting at points with positions 𝒓i relative to a reference point on the object o.

P = 𝑭i𝒗i (17.43)
= 𝑭i(𝒗o+𝝎×𝒓i/o)
= 𝒗o𝑭i𝟎+𝑭i(𝝎×𝒓i/o)
= 𝝎(𝒓i/o×𝑭i)
= 𝝎𝒓i/o×𝑭i
= 𝝎𝑴o

Work of a general force distribution

A general force distribution has, by reasoning close to that above, a power of:

P = 𝑭tot𝒗o+𝝎𝑴o. (17.44)

For a given force system applied to a given object in a given motion any point o can be used. The terms in the formula above will depend on o, but the sum does not. Besides the center-of-mass, another convenient locations for o is a fixed hinge, in which case the applied force has no power.

Box 17.4 The vector triple product 𝑨×(𝑩×𝑪)

The formula

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

can be verified by writing each of the vectors in terms of its orthogonal components (e.g., 𝐀=Axıˆ+Ayȷˆ+Az𝐤ˆ) and checking equality of the 27 terms on the two sides of the equations (only 12 are non-zero). If this 20 minute proof seems tedious it can be replaced by a more abstract geometric argument partly presented below that surely takes more than 20 minutes to grasp.

Geometry of the vector triple product


Filename:tfigure1-14



Because 𝑩×𝑪 is perpendicular to both 𝑩 and 𝑪 it is perpendicular to the plane of 𝑩 and 𝑪, that is, it is ‘normal’ to the plane BC. 𝑨×(𝑩×𝑪) is perpendicular to both 𝑨 and 𝑩×𝑪, so it is perpendicular to the normal to the plane of BC. That is, it must be in the plane of 𝑩 and 𝑪. But any vector in the plane of 𝑩 and 𝑪 must be a combination of 𝑩 and 𝑪. Also, the vector triple product must be proportional in magnitude to each of 𝑨, 𝑩 and 𝑪. Finally, the triple cross product of 𝑨×(𝑩×𝑪) must be the negative of 𝑨×(𝑪×𝑩) because 𝑩×𝑪=𝑪×𝑩. So the identity

𝑨×(𝑩×𝑪)=(𝑨𝑪)𝑩(𝑨𝑩)𝑪


is almost natural: The expression above is almost the only expression that is a linear combination of 𝑨 and 𝑩 that is linear in both, also linear in 𝑪 and switches sign if 𝑩 and 𝑪 are interchanged. These properties would be true if the whole expression were multiplied by any constant scalar. But a test of the equation with three unit vectors shows that such a multiplicative constant must be one. This reasoning constitutes an informal derivation of the identity 17.45.

Using the triple cross product in dynamics equations



We will use identity 17.45 for two purposes in the development of dynamics equations:

  1. 1.

    In the 2D expression for acceleration, the centripetal acceleration is given by 𝝎×(𝝎×𝑹) simplifies to ω2𝑹if𝝎𝑹. This equation follows by setting 𝑨=𝝎, 𝑩=𝝎 and 𝑪=𝑹 in equation 17.45 and using 𝑹𝝎=0 if 𝝎𝑹. In 3D 𝝎×(𝝎×𝒓) gives the vector shown in the lower figure below.

    Filename:tfigure1-f
    Filename:tfigure1-15


  2. 2.

    The term 𝒓×(𝝎×𝒓) will appear in the calculation of the angular momentum of a rigid body. By setting 𝑨=𝒓, 𝑩=𝝎 and 𝑪=𝒓, in equation 17.45 and use 𝒓𝒓=r2 because 𝒓𝒓 we get the useful result that 𝒓×(𝝎×𝒓)=r2𝝎(𝒓𝝎)𝒓.

SAMPLE 17.6

Filename:sfig8-7-2again
Figure 17.35:

Free planar motion. A rigid rod of length =1m and mass mr=1kg, and a rigid square plate of side 1m and mass mp=10kg are launched in motion on a frictionless plane (e.g., an ice hockey rink) with exactly the same initial velocities 𝒗cm(0)=10m/sıˆ+1m/sȷˆ and 𝝎(0)=1rad/s𝒌ˆ. Both the rod and the plate have their center-of-mass at the baseline at t=0.

  1. 1.

    Which of the two is farther from the base line in 3 seconds and which one has undergone a greater number of revolutions?

  2. 2.

    Find and draw the position of the bar at t=1sec and at t=3sec.


Solution

  1. 1.

    The free-body diagram of the rod is shown in fig. 17.36. There are no forces acting on the rod in the xy-plane. Although there is force of gravity and the normal reaction of the surface acting on the rod, these forces are inconsequential since they act normal to the xy-plane. Therefore, we do not include these forces in our free-body diagram .

    Filename:sfig8-7-2disks
    Figure 17.36:

    The linear momentum balance for the rod gives

    𝑭 = mr𝒂cm
    𝟎 = mr𝒗˙cm
     𝒗cm = 𝟎𝑑t=constant=𝒗cm0
     𝒓cm = 𝒗cm0𝑑t=𝒓cm0+𝒗cm0t (17.46)

    It is clear from the analysis above that in the absence of any applied forces, the position of the body depends only on the initial position and the initial velocity. Since both the plate and the rod start from the same base line with the same initial velocity, they travel the same distance from the base line during any given time period; mass or its geometric distribution plays no role in the motion. Thus the center-of-mass of the rod and the plate will be exactly the same distance (|𝒓cm(t)𝒓cm0|=|𝒗cm0t|) at time t. Similarly, the angular momentum balance about the center-of-mass of the rod gives

    𝑴cm = 𝑯˙cm
    𝟎 = Izzcm𝝎˙
     𝝎 = 𝟎𝑑t=constant=𝝎0=θ˙0𝒌ˆ
     θ = θ˙0𝑑t=θ0+θ˙0t (17.47)

    Thus the angular position of the body is also, as expected, independent of the mass and mass distribution of the body, and depends entirely on the initial position and the initial angular velocity. Therefore, both the rod and the plate undergo exactly the same amount of rotation (θ(t)θ0=θ˙0t) during any given time.

  2. 2.

    We can find the position of the rod at t=1s and t=3s by substituting these values of t in eqns. (17.46) and 17.47. For convenience, let us assume that 𝒓cm0=𝟎. From the initial configuration of the rod, we also know that θ0=0.

    𝒓cm(t=1s) = 𝒗cm0(1s)=(10m/sıˆ+1m/sȷˆ)(1s)10mıˆ+1mȷˆ
    𝒓cm(t=3s) = 𝒗cm0(3s)=30mıˆ+3mȷˆ.
    θ(t=1s) = θ˙0(1s)=(1rad/s)(1s)=1rad
    θ(t=3s) = θ˙0(3s)=3rad.

    Accordingly, we show the position of the rod in fig. 17.37.

    Filename:sfig8-4-4
    Figure 17.37:

SAMPLE 17.7

Filename:sfig8-4-4a
Figure 17.38:

A passive rigid diver. An experimental model of a rigid diver is to be launched from a diving board that is 3m above the water level. Say that the initial velocity of the center-of-mass and the initial angular velocity of the diver can be controlled at launch. The diver is launched into the dive in almost vertical position, and it is required to be as vertical as possible at the very end of the dive (which is marked by the position of the diver’s center-of-mass at 1m above the water level). If the initial vertical velocity of the diver’s center-of-mass is 3m/s, find the required initial angular velocity for the vertical entry of the diver into the water.


Solution

Filename:sfig8-4-4b
Figure 17.39:

Once the diver leaves the diving board, it is in free flight under gravity, i.e., the only force acting on it is the force due to gravity. The free-body diagram of the diver is shown in fig. 17.39. The linear momentum balance for the diver gives

𝑭 = m𝒂cm
mgȷˆ = my¨ȷˆ
 y¨ = g
𝑴cm = 𝑯˙cm
𝟎 = Izzcmθ¨𝒌ˆ
 θ¨ = 0.

From these equations of motion, it is clear that the linear and the angular motions of the diver are uncoupled. We can easily solve the equations of motion to get

y(t) = y0+y˙0t12gt2
θ(t) = θ0+θ˙0t.

We need to find the initial angular speed θ˙0 such that θ=π when y=1m (the center-of-mass position at the water entry). From the expression for θ(t), we get, θ˙0=π/t. Thus we need to find the value of t at the instant of water entry. We can find t from the expression for y(t) since we know that y=1m at that instant, and that y0=3m and y˙0=3m/s. We have,

y = y0+y˙012gt2
 t = y˙0±y˙02+2g(y0y)g
= 3m/s±(3m/s)2+29.8m/s2(3m1m)9.8m/s2
= 1.15 or 0.53s.

We reject the negative value of time as meaningless in this context. Thus it takes the diver 1.15s to complete the dive. Since, the diver must rotate by π during this time, we have

θ˙0=π/t=π/(1.15s)=2.73rad/s.

Answer: θ˙0=2.73rad/s

SAMPLE 17.8

Filename:sfig8-4-4c
Figure 17.40:

A plate tumbling in space. A rectangular plate of mass m=0.5kg,Izzcm=2.08×103kgm2, and dimensions a=200mm and b=100mm is pushed by a force 𝑭=0.5Nıˆ, acting at d=30mm away from the mass-center, as shown in the figure. Assume that the force remains constant in magnitude and direction but remains attached to the material point P of the plate. There is no gravity.

  1. 1.

    Find the initial acceleration of the mass-center.

  2. 2.

    Find the initial angular acceleration of the plate.

  3. 3.

    Write the equations of motion of the plate (for both linear and angular motion).


Solution

The only force acting on the plate is the applied force 𝑭. Thus, fig. 17.40 is also the free-body diagram of the plate at the start of motion.

  1. 1.

    From the linear momentum balance we get,

    𝑭 = m𝒂cm
     𝒂cm = 𝑭m=0.5Nıˆ0.5kg=1m/s2ıˆ.

    which is the initial acceleration of the mass-center.

    Answer: acm=1m/s2ıˆ

  2. 2.

    From the angular momentum balance about the mass-center, we get

    𝑴cm = 𝑯˙cm
    Fd𝒌ˆ = Izzcm𝝎˙
     𝝎˙ = FdIzzcm𝒌ˆ=0.5N0.03m2.08kgm2=7.2rad/s2𝒌ˆ

    which is the initial angular acceleration of the plate.

    Answer: ω˙=7.2rad/s2kˆ

  3. 3.

    To find the equations of motion, we can use the linear momentum balance and the angular momentum balance as we have done above. So, why aren’t the equations obtained above for the linear acceleration, 𝒂cm=F/mıˆ, and the angular acceleration, 𝝎˙=Fd/Izzcm𝒌ˆ, qualified to be called equations of motion? Because, they are not valid for a general configuration of the plate during its motion. The expressions for the accelerations are valid only in the initial configuration (and hence those are initial accelerations).

    Let us first draw a free-body diagram of the plate in a general configuration during its motion (see fig. 17.41). Assume the center-of-mass to be displaced by xıˆ and yȷˆ, and the longitudinal axis of the plate to be rotated by θ𝒌ˆ with respect to the vertical (inertial y-axis). The applied force remains horizontal and attached to the material point P, as stated in the problem.

    Filename:pfigure-s94h13p2
    Figure 17.41: Free-body diagram of the plate at some instant t when the longitudinal axis of the plate makes an angle θ(t) with the fixed vertical axis.

    The linear momentum balance gives

    𝑭=m𝒂cm 𝒂cm=𝑭m
    or x¨ıˆ+y¨ȷˆ=Fmıˆ x¨=Fm,y¨=0.

    Since F/m is constant, the equations of motion of the center-of-mass indicate that the acceleration is constant and that the mass-center moves in the x-direction.

    Similarly, we now use angular momentum balance to determine the rotation (angular motion) of the plate. The angular momentum balance about the mass-center give

    Filename:pfigure-f93f5
    Figure 17.42: Geometry of the plate at the instant t when the longitudinal axis of the plate makes an angle θ(t) with the fixed y-axis. The position of point P is 𝒓P/cm which makes a fixed angle α(=tan1db/2) with the transverse axis of the plate. This angle is shown here merely for ease of calculation.
    𝑴cm = 𝑯˙cm
    𝒓P/cm×𝑭 = Izzcmθ¨𝒌ˆ.
      Now,
    𝒓P/cm = r[cos(θ+α)ıˆ+sin(θ+α)ȷˆ]
    𝑭 = Fıˆ
     𝒓P/cm×𝑭 = Frsin(θ+α)𝒌ˆ.
      Thus,
    θ¨ = FrIzzcmsin(θ+α)

    where r=d2+(b/2)2 and α=tan1(2d/b).

    Thus, we have got the equations of motion for both the linear and the angular motion.

    Answer: x¨=Fm,y¨=0,θ¨=FrIzzcmsin(θ+α)

  4. 4.

    The equations of linear motion of the plate are very simple and we can solve them at once to get

    x(t) = x0+x˙0t+12Fmt2
    y(t) = y0+y˙0t.

    If the plate starts from rest (x˙0=0,y˙0=0) with the center-of-mass at the origin (x0=0,y0=0), then we have

    x(t)=F2mt2,andy(t)=0.

    Thus the center-of-mass moves along the x-axis with acceleration F/m.

    The equation of angular motion of the plate is not so simple. In fact, it is a nonlinear ODE. It is difficult to get an analytical solution of this equation. However, we can solve it numerically using, say, a Runge-Kutta ODE solver:

    Filename:pfigure-s94h13p3
    Figure 17.43: Variation of plate orientation θ with time.
    ODEs = {thetadot = w, wdot = (F*r/Icm)*sin(theta+a)}
    IC   = {theta(0) = 0, w(0) = 0}
    Set  F=.5, d=0.03; b=0.1; Icm=2.08e-03
    compute r = sqrt(d^2+.25*b^2),  a = atan(2*d/b)
    Solve ODEs with IC for t=0 to t=10
    Plot theta(t)
    

    The plot obtained from this calculation is shown in fig. 17.43 which shows that the plate oscillates in its plane as the center of mass races along the x-axis.

The trajectory of point P can now be obtained from 𝒓P=𝒓cmrcos(θ+α)ıˆrsin(θ+α)ȷˆ. Using this relationship, the trajectory of point P, along with the trajectory of the center of mass, is shown in fig. 17.44 for the first 3 seconds.

Filename:p-s96-p3-3
Figure 17.44: Motion of the plate during the first 3 seconds: Position of the plate at various time instants along with the trajectories of the center of mass and point P.

SAMPLE 17.9  Impulse-momentum. Consider the plate problem of Sample 17.39 (page 17.39) again. Assume that the plate is at rest at t=0 in the vertical upright position and that the force acts on the plate for 2 seconds.

  1. 1.

    Find the velocity of the center-of-mass of the plate at the end of 2 seconds.

  2. 2.

    Can you also find the angular velocity of the plate at the end of 2 seconds?


Solution

  1. 1.

    Since we are interested in finding the velocity at a particular instant t, given the velocity at another instant t=0, we can use the impulse-momentum equations to find the desired velocity.

    𝑳2𝑳1 = t1t2𝑭dt
    m𝒗cm(t)m𝒗cm(0) = 0t𝑭𝑑t
     𝒗cm(t) = 𝒗cm(0)+1m02𝑭𝑑t
    = 𝟎+10.5kg02(0.5Nıˆ)𝑑t
    = 2m/sıˆ.

    Answer: vcm(2s)=2m/sıˆ

  2. 2.

    Now, let us try to find the angular velocity the same way, using angular impulse-momentum relation. We have,

    (𝑯cm)2(𝑯cm)1 = t1t2𝑴cmdt
    Izzcm𝝎(t)Izzcm𝝎(0) = 0t𝑴cmdt
     𝝎(t) = 𝝎(0)+1Izzcm0t𝑴cmdt
    = 𝝎(0)+1Izzcm0t(𝒓P/cm×𝑭)𝑑t
    = 𝟎+1Izzcm0t(Frsin(θ+α)𝒌ˆ)𝑑t
    = FrIzzcm(0tsin(θ+α)𝑑t)𝒌ˆ.

    Now, we are in trouble; how do we evaluate the integral? In the integrand, we have θ which is an implicit function of t. Unless we know how θ depends on t we cannot evaluate the integral. To find θ(t) we have to solve the equation of angular motion we derived in the previous sample. However, we were not able to solve for θ(t) analytically; we had to resort to a numerical solution. Thus, it is not possible to evaluate the integral above and, therefore, we cannot find the angular velocity of the plate at the end of 2 seconds using impulse-momentum equations. We could, however, find the desired velocity easily from the numerical solution.

Problems for 17.2 Dynamics of a rigid object

17.2.1  The uniform rectangle of width a=1m, length b=2m, and mass m=1kg in the figure is sliding on the xy-plane with no friction. At the moment in question, point C is at xC=3m and yC=2m. The linear momentum is 𝑳=4ıˆ+3ȷˆ(kgm/s) and the angular momentum about the center of mass is 𝑯cm=5𝒌ˆ(kgm/s2). Find the acceleration of any point on the body that you choose. (Mark it.) [Hint: You have been given some redundant information.]

Filename:bikefork1-ang-mom
Figure 17.45

17.2.2  The vertical pole AB of mass m and length is initially, at rest on a frictionless surface. A tension T is suddenly applied at A. What is x¨cm? What is θ¨AB? What is x¨B? Gravity may be ignored.

Filename:bikefork-ang-mom
Figure 17.46

17.2.3  Force on a stick in space. 2-D . No gravity. A uniform thin stick with length and mass m is, at the instant of interest, parallel to the y axis and has no velocity and no angular velocity. The force 𝑭=Fıˆ with F>0 is suddenly applied at point A. The questions below concern the instant after the force 𝑭 is applied.

  1. (a)

    What is the acceleration of point C, the center of mass? Answer: acm=Fmıˆ.

  2. (b)

    What is the angular acceleration of the stick? Answer: α=6Fmkˆ.

  3. (c)

    What is the acceleration of the point A? Answer: aA=4Fmıˆ.

  4. (d)

    (relatively harder) What additional force would have to be applied to point B to make point B’s acceleration zero? Answer: FB=F2ıˆ

Filename:summer95f-5-a
Figure 17.47

17.2.4  A uniform thin rod of length and mass m stands vertically, with one end resting on a frictionless surface and the other held by someone’s hand. The rod is released from rest, displaced slightly from the vertical. No forces are applied during the release. There is gravity.

  1. (a)

    Find the path of the center of mass.

  2. (b)

    Find the force of the floor on the end of the rod just before the rod is horizontal.

17.2.5   A uniform disk, with mass center labeled as point G, is sitting motionless on the frictionless xy plane. A massless peg is attached to a point on its perimeter. This disk has radius of 1m and mass of 10kg. A constant force of F=1000Nı is applied to the peg for .0001s (one ten-thousandth of a second).

  1. (a)

    What is the velocity of the center of mass of the center of the disk after the force is applied?

  2. (b)

    Assuming that the idealizations named in the problem statement are exact is your answer to (a) exact or approximate?

  3. (c)

    What is the angular velocity of the disk after the force is applied?

  4. (d)

    Assuming that the idealizations named in the problem statement are exact is your answer to (c) exact or approximate?

17.2.6  A uniform thin flat disc is floating in space. It has radius R and mass m. A force F is applied to it at a distance d from the center in the y direction. Treat this problem as two-dimensional.

  1. (a)

    What is the acceleration of the center of the disc? Answer: a=Fmȷˆ.

  2. (b)

    What is the angular acceleration of the disk? Answer: α=FdmR22kˆ.

Filename:pfigure4-2-rp10
Figure 17.48

17.2.7  A uniform 1kg plate that is one meter on a side is initially at rest in the position shown. A constant force 𝑭=1Nıˆ is applied at t=0 and maintained henceforth. If you need to calculate any quantity that you don’t know, but can’t do the calculation to find it, assume that the value is given.

  1. (a)

    Find the position of G as a function of time (the answer should have numbers and units).

  2. (b)

    Find a differential equation, and initial conditions, that when solved would give θ as a function of time. θ is the counterclockwise rotation of the plate from the configuration shown.

  3. (c)

    Write computer commands that would generate a drawing of the outline of the plate at t=1s. You can use hand calculations or the computer for as many of the intermediate commands as you like. Hand work and sketches should be provided as needed to justify or explain the computer work.

  4. (d)

    Run your code and show clear output with labeled plots. Mark output by hand to clarify any points.

Filename:pfigure-blue-125-2
Figure 17.49

17.2.8  A uniform rectangular metal beam of mass m hangs symmetrically by two strings as shown in the figure.

  1. (a)

    Draw a free-body diagram of the beam and evaluate 𝑭.

  2. (b)

    Repeat (a) immediately after the left string is cut.

Filename:pfigure-blue-68-1
Figure 17.50

17.2.9  A uniform slender bar AB of mass m is suspended from two springs (each of spring constant K) as shown. Immediately after spring 2 breaks, determine

  1. (a)

    the angular acceleration of the bar,

  2. (b)

    the acceleration of point A, and

  3. (c)

    the acceleration of point B.

Filename:pfigure-blue-58-1
Figure 17.51

17.2.10  Two small spheres A and B are connected by a rigid rod of length =1.0ft and negligible mass. The assembly is hung from a hook, as shown. Sphere A is struck, suddenly breaking its contact with the hook and giving it a horizontal velocity v0=3.0ft/s which sends the assembly into free fall. Determine the angular momentum of the assembly about its mass center at point G immediately after A is hit. After the center of mass has fallen two feet, determine:

  1. (a)

    the angle θ through which the rod has rotated,

  2. (b)

    the velocity of sphere A,

  3. (c)

    the total kinetic energy of the assembly of spheres A and B and the rod, and

  4. (d)

    the acceleration of sphere A.

Filename:pfigure-blue-157-1
Figure 17.52

The next several problems concern Work, power and energy

17.2.11  Verify that the expressions for work done by a force F, W=FΔS, and by a moment M, W=MΔθ, are dimensionally correct if ΔS and Δθ are linear and angular displacements respectively.

17.2.12  A uniform disc of mass m and radius r rotates with angular velocity ω𝒌ˆ. Its center of mass translates with velocity 𝒗=x˙ıˆ+y˙ȷˆ in the xy-plane. What is the total kinetic energy of the disk?

17.2.13

Calculate the energy stored in a spring using the expression EP=12kδ2 if the spring is compressed by 100 mm and the spring constant is 100 N/m.

17.2.14  In a rack and pinion system, the rack is acted upon by a constant force F=50N and has speed v=2m/s in the direction of the force. Find the input power to the system.

Filename:summer95f-5
Figure 17.53

17.2.15  The driving gear in a compound gear train rotates at constant speed ω0. The driving torque is Min. If the driven gear rotates at a constant speed ωout, find:

  1. (a)

    the input power to the system, and

  2. (b)

    the output torque of the system assuming there is no power loss in the system; i.e., power in = power out.

Filename:pfigure-blue-127-2
Figure 17.54

17.2.16  An elaborate frictionless gear box has an input and output roller with Vin=const. Assuming that Vout=7Vin and the force between the left belt and roller is Fin=3lb:

  1. (a)

    What is Fout (draw a picture defining the signs of Fin and Fout)?. Answer: Fout=37lb.

  2. (b)

    Is Fout greater or less than the Fin? (Assume Fin>0.) Why? Answer: Fout is always less than the Fin.

Filename:pfigure-s94h13p4
Figure 17.55

17.3 Kinematics of rolling and sliding

Pure rolling in 2-D

In this section, we would like to add to the vocabulary of special motions by considering pure rolling. Most commonly, one discusses pure rolling of round objects on flat ground, like wheels and balls, and rolling of round things on other round things like gears and cams.

2-D rolling of a round wheel on level ground

The simplest case, the no-slip rolling of a round wheel, is an instructive starting point. First, we define the geometric and kinematic variables as shown in fig. 17.56.

Filename:pfigure-blue-90-2
Figure 17.56: Pure rolling of a round wheel on a level support.

For convenience, we pick a point D which was at xD=0 at the start of rolling, when xC=0. The key to the kinematics is that:

The arc length traversed on the wheel is the distance traveled by the wheel center.

That is,

xC = sD
= Rϕ
vC=x˙C = Rϕ˙
aC=v˙C=x¨C = Rϕ¨

So the rolling condition amounts to the following set of restrictions on the position of C, 𝒓C, and the rotations of the wheel ϕ:

𝒓C=Rϕıˆ+Rȷˆ,𝒗C=Rϕ˙ıˆ,𝒂C=Rϕ¨ıˆ,𝝎=ϕ˙𝒌ˆ, and𝜶=𝝎˙=ϕ¨𝒌ˆ.

If we want to track the motion of a particular point, say D, we could do so by using the following parametric formula:

𝒓D = 𝒓C+𝒓D/C
= R(ϕıˆ+ȷˆ)+R(sinϕıˆcosϕȷˆ)
= R[(ϕsinϕ)ıˆ+(1cosϕ)ȷˆ)]
𝒗D = R[(ϕ˙(1cosϕ)ıˆ+ϕ˙sinϕȷˆ)] (17.48)
𝒂D =assuming ϕ˙= constant Rϕ˙2(sinϕıˆ+cosϕȷˆ).

Note that if ϕ=0 or 2π or 4π, etc., then the point D is on the ground and eqn. (17.48) correctly gives that

𝒗D=R[ϕ˙(1cos(2nπ1)0]ıˆ+ϕ˙sin(2nπ)0ȷˆ=𝟎.

Instantaneous kinematics

Instead of tracking the wheel from its start, we could analyze the kinematics at the instant of interest. Here, we make the observation that the wheel rolls without slip. Therefore, the point on the wheel touching the ground has no velocity relative to the ground.

𝒗AVelocity of point on the wheel touching the ground=𝒗BVelocity of ground =  0 (17.49)

Now, we know how to calculate the velocity of points on a rigid body. So,

𝒗A=𝒗C+𝒗A/C,

where, since A and C are on the same rigid body (fig. 17.56), we have from eqn. (16.2) that

𝒗A/C=𝝎×𝒓A/C.

Putting this equation together with eqn. (17.49), we get

𝒗A = 𝒗B
𝒗CvCıˆ+𝝎ω𝒌ˆ×𝒓A/CRȷˆ = 𝟎
vCıˆ+ωRıˆ = 𝟎
vC=ωR. (17.50)

We use 𝒗C=vıˆ since the center of the wheel goes neither up nor down.

Note that if you measure the angle by ϕ, like we did before, then 𝝎=ϕ˙𝒌ˆ so that positive rotation rate is in the counter-clockwise direction. Thus, vC=ωR=(ϕ˙)R=ϕ˙R.

Since there is always some point of the wheel touching the ground, we know that vC=ωR for all time. Therefore,

𝒂C=v˙Cıˆ=ω˙Rıˆ.

Rolling of round objects on round surfaces

For round objects rolling on or in another round object, the analysis is similar to that for rolling on a flat surface. Common applications are the so-called epicyclic, hypo-cyclic, or planetary gears (See Box 17.3 on planetary gears on page 17.3).

Filename:s92f1p7
Figure 17.57: One round object rolling on another round object.

Referring to fig. 17.57, we can calculate the velocity of C with respect to a fixed frame two ways and compare:

𝒗C = 𝒗B+𝒗C/B
𝒗C = 𝒗A𝟎+𝒗B/A𝟎𝒗B+𝒗C/B.
θ˙(R1+R2)𝒆ˆθ = ωR2𝒆ˆθ
ω = θ˙(R1+R2)R2=θ˙(1+R1R2).

Example: Two quarters.

The formula above can be tested in the case of R1=R2 by using two quarters or two dimes on a table. Roll one quarter, call it , around another quarter pressed fast to the table. You will see that as the rolling quarter travels around the stationary quarter one time, it makes two full revolutions. That is, the orientation of changes twice as fast as the angle of the line from the center of the stationary quarter to its center. Or, in the language of the calculation above, ω=2θ˙.

Sliding

Although wheels and balls are known for rolling, they do sometimes slide such as when a car screeches at fast acceleration or sudden braking or when a bowling ball is released on the lane.

The sliding velocity is the velocity of the material point on the wheel (or ball) relative to its contacting substrate. In the case of pure rolling, the sliding velocity is zero. In the case of a ball or wheel moving against a stationary support surface, whether flat or curved, the sliding velocity is

𝒗sliding=𝒗circlecenter+𝝎×𝒓contact/center (17.51)

Example: Bowling ball

Filename:twodisks-ang-mom
Figure 17.58: The bowling ball is sliding so long as vGωR

The velocity of that point on the bowling ball which is instantaneously in contact with the alley (ground) is 𝒗C=vGıˆ+ω𝒌ˆ×𝒓C/G=(vG+ωR)ıˆ. So unless ω=vG/R the ball is sliding.

Note that, if sliding, the friction force on the ball opposes the slip of the ball and tends to accelerate the ball’s rotation towards rolling. That is, for example, if the ball is not rotating the sliding velocity is vGıˆ, the friction force is in the ıˆ direction and angular momentum balance about the center-of-mass implies ω˙<0 and a clockwise rotational acceleration. No matter what the initial velocity or rotational rate the ball will eventually roll.

Box 17.5 The Sturmey-Archer hub

In 1903, the year the Wright Brothers first flew powered airplanes, the Sturmey-Archer company patented the internal-hub three-speed bicycle transmission. This marvel of engineering was sold on the best bikes until finicky but fast racing bicycles using derailleurs started to push them out of the market in the 1960’s. Now, a hundred years later, internal bicycle hubs (now made by Shimano and Sachs) are having something of a revival, particularly in Europe. These internal-hub transmissions utilize a system called planetary gears, gears which roll around other gears. See the figure below.

Filename:tfigure7-2D-planetary

In order to understand this gear system, we need to understand its kinematics—the motion of its parts. The central ‘sun’ gear is stationary, at least we treat it as stationary in this discussion since it is fixed to the bike frame, so it is fixed in body . The ‘planet’ gears roll around the sun gear. Let’s call one of these planets 𝒫. The spider 𝒮 connects the centers of the rolling planets. Finally, the ring gear rotates around the sun.

Filename:tfigure7-2D-planetary-b

The gear transmission steps up the angular velocity when the spider 𝒮 is driven and ring , which moves faster, is connected to the wheel. The transmission steps down the angular velocity when the ring gear is driven and the slower spider is connected to the wheel. The third ‘speed’ in the three-speed gear transmission is direct drive (the wheel is driven directly).

What are the ‘gear ratios’ in the planetary gear system? The ‘trick’ is to recognize that for rolling contact that the contacting points have the same velocity, 𝒗A=𝒗B and 𝒗D=𝒗E. Let’s define some terms.

𝝎𝒮 = ω𝒮𝒌ˆ angular velocity of the spider
𝝎𝒫 = ω𝒫𝒌ˆ angular velocity of the planet
𝝎 = ω𝒌ˆ angular velocity of the ring

Now, we can find the relation of these angular velocities as follows. Look at the velocity of point C in two ways. First,

𝒗CA point on the spider = 𝒗CA point on the planetary gear
𝝎𝒮×𝒓C = 𝒗B𝟎+𝝎𝒫×𝒓C/B
ω𝒮rC = ω𝒫RP
ω𝒫 = rCRPω𝒮 (17.52)

Next, let’s look at point D and E:

𝒗D = 𝒗E
𝒗A+𝒗D/A = 𝝎×𝒓R
𝟎+𝝎𝒫×𝒓D/A = ω𝒌ˆ×𝒓R
ω𝒫(2RP)𝒆ˆθ = ωrR𝒆ˆθ
ω = 2RPrRω𝒫ω𝒫=rCRPω𝒮
ω = 2RPrRrCrC=rS+RPRPω𝒮
= 2(rS+RP)rRrR=rS+2RPω𝒮
ωω𝒮 = 21+RPrS1+2RPrS= angular velocity step-up.
Filename:tfigure7-2D-planetary-a

Typically, the gears have radius ratio of RPrS=32 which gives a gear ratio of 54. Thus, the ratio of the highest gear to the lowest gear on a Sturmey-Archer hub is 54/45=2516=1.5625. You might compare this ratio to that of a modern mountain bike, with eighteen or twenty-one gears, where the ratio of the highest gear to the lowest is about 4:1.

SAMPLE 17.10

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

Falling ladder: The ends of a ladder of length L=3m slip along the frictionless wall and floor shown in Figure 17.59. At the instant shown, when θ=60, the angular speed θ˙=1.15rad/s and the angular acceleration θ¨=2.5rad/s2. Find the absolute velocity and acceleration of end B of the ladder.


Solution

Since the ladder is falling, it is rotating clockwise. From the given information:

𝝎 = θ˙𝒌ˆ=1.15rad/s𝒌ˆ
𝝎˙ = θ¨𝒌ˆ=2.5rad/s2𝒌ˆ.

We need to find 𝒗B, the absolute velocity of end B, and 𝒂B, the absolute acceleration of end B.

Since the end A slides along the wall and end B slides along the floor, we know the directions of 𝒗A,𝒗B,𝒂Aand𝒂B.

Let 𝒗A=vAȷˆ,𝒂A=aAȷˆ,𝒗B=vBıˆand𝒂B=aBıˆ where the scalar quantities vA,aA,vBandaB are unknown.

Now,𝒗A = 𝒗B+𝒗A/B=𝒗B+𝝎×𝒓A/B
orvAȷˆ = vBıˆ+θ˙𝒌ˆ×L(cosθıˆsinθȷˆ)𝒓A/B
= (vB+θ˙Lsinθ)ıˆθ˙Lcosθȷˆ.
Filename:pfigure-blue-110-1
Figure 17.60:

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

vAȷˆıˆ0 = (vB+θ˙Lsinθ)ıˆıˆ1+θ˙Lcosθȷˆıˆ0
 0 = vB+θ˙Lsinθ
 vB = θ˙Lsinθ=(1.15rad/s)3m32
= 2.99m/s.

Answer: 𝒗B=2.9m/sıˆ

Similarly,

𝒂A = 𝒂B+𝝎˙×𝒓A/B+𝝎×(𝝎×𝒓A/B)ω2𝒓A/B
aAȷˆ = aBıˆ+θ¨𝒌ˆ×L(cosθıˆsinθȷˆ)θ˙2L(cosθıˆsinθȷˆ)
= (aB+θ¨Lsinθ+θ˙2Lcosθ)ıˆ+(θ¨Lcosθ+θ˙2Lsinθ)ȷˆ.

Dotting both sides of this equation with ıˆ (as we did for velocity) we get:

0 = aB+θ¨Lsinθ+θ˙2Lcosθ
 aB = θ¨Lsinθθ˙2Lcosθ
= (2.5rad/s23m32)(1.15rad/s)23m12
= 4.51m/s2.

Answer: 𝒂B=4.51m/s2

SAMPLE 17.11

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

A cylinder of diameter 500mm rolls down an inclined plane with uniform acceleration (of the center-of-mass) a=0.1m/s2. At an instant t0, the mass-center has speed v0=0.5m/s.

  1. 1.

    Find the angular speed ω and the angular acceleration ω˙ at t0.

  2. 2.

    How many revolutions does the cylinder make in the next 2 seconds?

  3. 3.

    What is the distance travelled by the center-of-mass in those 2 seconds?


Solution

This problem is about simple kinematic calculations. We are given the velocity, x˙, and the acceleration, x¨, of the center-of-mass. We are supposed to find angular velocity ω, angular acceleration ω˙, angular displacement θ in 2 seconds, and the corresponding linear distance x along the incline. The radius of the cylinder R=diameter/2=0.25m.

  1. 1.

    From the kinematics of pure rolling,

    ω = x˙R=0.5m/s0.25m=2rad/s,
    ω˙ = x¨R=0.1m/s20.25m=0.4rad/s2.

    Answer: ω=2rad/s,ω˙=0.4rad/s2

  2. 2.

    We can find the number of revolutions the cylinder makes in 2 seconds by solving for the angular displacement θ in this time period. Since,

    θ¨ω˙=constant,

    we integrate this equation twice and substitute the initial conditions, θ˙(t=0)=ω=2rad/s and θ(t=0)=0, to get

    θ(t) = ωt+12ω˙t2
     θ(t=2s) = (2rad/s)(2s)+12(0.4rad/s)(4s2)
    = 4.8rad=4.82πrev=0.76rev.

    Answer: θ=0.76rev

  3. 3.

    Now that we know the angular displacement θ, the distance travelled by the mass-center is the arc-length corresponding to θ, i.e.,

    x=Rθ=(0.25m)(4.8)=1.2m.

    Answer: x=1.2m

    Note that we could have found the distance travelled by the mass-center by integrating the equation x¨=0.1m/s2 twice.

SAMPLE 17.12  Condition of pure rolling. A cylinder of radius R=20cm rolls on a flat surface with absolute angular speed ω=12rad/s under the conditions shown in the figure (In cases (ii) and (iii), you may think of the ‘flat surface’ as a conveyor belt). In each case,

  1. 1.

    Write the condition for pure rolling.

  2. 2.

    Find the velocity of the center C of the cylinder.

Filename:pfigure-s94h14p4
Figure 17.62:

Solution

Filename:pfigure-blue-112-1
Figure 17.63: The cylinder rolls on the flat surface. Instantaneously, point P on the cylinder is in contact with point Q on the flat surface. For pure rolling, points P and Q must have the same velocity.

At any instant during rolling, the cylinder makes a point-contact with the flat surface. Let the point of instantaneous contact on the cylinder be P, and let the corresponding point on the flat surface be Q. The condition of pure rolling, in each case, is 𝒗P=𝒗Q, that is, there is no relative motion between the two contacting points (a relative motion will imply slip). Now, we analyze each case.

Case(i)

In this case, the bottom surface is fixed. Therefore,

  1. 1.

    The condition of pure rolling is: 𝒗P=𝒗Q=𝟎.

  2. 2.

    Velocity of the center:

    𝒗C = 𝒗P+𝝎×𝒓C/P=𝟎+(ω𝒌ˆ)×Rȷˆ
    = ωRıˆ=(12rad/s)(0.2m)ıˆ=2.4m/sıˆ.
Case(ii)

In this case, the bottom surface moves with velocity 𝒗=1m/sıˆ. Therefore, 𝒗Q=1m/sıˆ. Thus,

  1. 1.

    The condition of pure rolling is: 𝒗P=𝒗Q=1m/sıˆ.

  2. 2.

    Velocity of the center:

    𝒗C = 𝒗P+𝝎×𝒓C/P=v0ıˆ+ωRıˆ
    = 1m/sıˆ+2.4m/sıˆ=3.4m/sıˆ.
Case(iii)

In this case, the bottom surface moves with velocity 𝒗=1m/sıˆ. Therefore, 𝒗Q=1m/sıˆ. Thus,

  1. 1.

    The condition of pure rolling is: 𝒗P=𝒗Q=1m/sıˆ.

  2. 2.

    Velocity of the center:

    𝒗C = 𝒗P+𝝎×𝒓C/P=v0ıˆ+ωRıˆ
    = 1m/sıˆ+2.4m/sıˆ=1.4m/sıˆ.

Answer: (a):(i)𝒗P=𝟎,(ii)𝒗P=1m/sıˆ,(iii)𝒗P=1m/sıˆ,(b):(i)𝒗C=2.4m/sıˆ,(ii)𝒗C=3.4m/sıˆ,(iii)𝒗C=1.4m/sıˆ

SAMPLE 17.13

Filename:tfigure8-alt-app2c
Figure 17.64: A uniform disk of radius r rolls without slipping inside a fixed cylinder.

Motion of a point on a disk rolling inside a cylinder. A uniform disk of radius r rolls without slipping with constant angular speed ω inside a fixed cylinder of radius R. A point P is marked on the disk at a distance (<r) from the center of the disk. At a general time t during rolling, find

  1. 1.

    the position of point P,

  2. 2.

    the velocity of point P, and

  3. 3.

    the acceleration of point P


Solution

Let the disk be vertically below the center of the cylinder at t=0s such that point P is vertically above the center of the disk (Fig. 17.65). At this instant, Q is the point of contact between the disk and the cylinder. Let the disk roll for time t such that at instant t the line joining the two centers (line OC) makes an angle ϕ with its vertical position at t=0s. Since the disk has rolled for time t at a constant angular speed ω, point P has rotated counter-clockwise by an angle θ=ωt from its original vertical position P’.

Filename:sfig4-6-4a
Figure 17.65: Geometry of motion: keeping track of point P while the disk rolls for time t, rotating by angle θ=ωt inside the cylinder.
  1. 1.

    Position of point P: From Fig. 17.65(b) we can write

    𝒓P=𝒓C+𝒓P/C=(Rr)𝝀ˆOC+𝝀ˆCP

    where

    𝝀ˆOC = a unit vector along OC =sinϕıˆcosϕȷˆ,
    𝝀ˆCP = a unit vector along CP =sinθıˆ+cosθȷˆ.

    Thus,

    𝒓P=[(Rr)sinϕsinθ]ıˆ+[(Rr)cosϕ+cosθ]ȷˆ.

    We have thus obtained an expression for the position vector of point P as a function of ϕ and θ. Since we also want to find velocity and acceleration of point P, it will be nice to express 𝒓P as a function of t. As noted above, θ=ωt; but how do we find ϕ as a function of t? Note that the center of the disk C is going around point O in circles with angular velocity ϕ˙𝒌ˆ. The disk, however, is rotating with angular velocity 𝝎=ω𝒌ˆ about the instantaneous center of rotation, point D. Therefore, we can calculate the velocity of point C in two ways:

    𝒗C = 𝒗C
    or 𝝎×𝒓C/D = ϕ˙𝒌ˆ×𝒓C/O
    or ω𝒌ˆ×r(𝝀ˆOC) = ϕ˙𝒌ˆ×(Rr)𝝀ˆOC
    or ωr(𝒌ˆ×𝝀ˆOC) = ϕ˙(Rr)(𝒌ˆ×𝝀ˆOC)
     rRrω = ϕ˙.

    Integrating the last expression with respect to time, we obtain

    ϕ=rRrωt.

    Let

    q=rRr,

    then, the position vector of point P may now be written as

    𝒓P=[(Rr)sin(qωt)sin(ωt)]ıˆ+[(Rr)cos(qωt)+cos(ωt)]ȷˆ. (17.53)
  2. 2.

    Velocity of point P: Differentiating Eqn. (17.53) once with respect to time we get

    𝒗P=ω[(Rr)qcos(qωt)+cos(ωt)]ıˆ+ω[(Rr)qsin(qωt)sin(ωt)]ȷˆ.

    Substituting (Rr)q=r in 𝒗P we get

    𝒗P=ωr[{cos(qωt)+rcos(ωt)}ıˆ{sin(qωt)rsin(ωt)}ȷˆ]. (17.54)
  3. 3.

    Acceleration of point P: Differentiating Eqn. (17.54) once with respect to time we get

    𝒂P=ω2r[{qsin(qωt)+rsin(ωt)}ıˆ{qcos(qωt)rcos(ωt)}ȷˆ]. (17.55)

Answer: rP=[(Rr)sin(qωt)sin(ωt)]ıˆ+[(Rr)cos(qωt)+cos(ωt)]ȷˆvP=ωr[{cos(qωt)+rcos(ωt)}ıˆ{sin(qωt)rsin(ωt)}ȷˆ]aP=ω2r[{qsin(qωt)+rsin(ωt)}ıˆ{qcos(qωt)rcos(ωt)}ȷˆ]

SAMPLE 17.14  The rolling disk: instantaneous kinematics. For the rolling disk in Sample 17.63, let R=4ft,r=1ft and point P be on the rim of the disk. Assume that at t=0, the center of the disk is vertically below the center of the cylinder and point P is on the vertical line joining the two centers. If the disk is rolling at a constant speed ω=πrad/s, find

  1. 1.

    the position of point P and center C at t=1s, 3s, and 5.25s,

  2. 2.

    the velocity of point P and center C at those instants, and

  3. 3.

    the acceleration of point P and center C at the same instants as above.

Draw the position of the disk at the three instants and show the velocities and accelerations found above.


Solution

The general expressions for position, velocity, and acceleration of point P obtained in Sample 17.63 can be used to find the position, velocity, and acceleration of any point on the disk by substituting an appropriate value of in equations (17.53), (17.54), and (17.55). Since R=4r,

q=rRr=13.

Now, point P is on the rim of the disk and point C is the center of the disk. Therefore,

for point P:  = r,
for point C:  = 0.

Substituting these values for , and q=1/3 in equations (17.53), (17.54), and (17.55) we get the following.

  1. 1.

    Position:

    𝒓C = 3r[sin(ωt3)ıˆ+cos(ωt3)ȷˆ],
    𝒓P = 𝒓C+r[sin(ωt)ıˆ+cos(ωt)ȷˆ].
  2. 2.

    Velocity:

    𝒗C = ωr[cos(ωt3)ıˆsin(ωt3)ȷˆ],
    𝒗P = ωr[{cos(ωt3)+cos(ωt)}ıˆ{sin(ωt3)sin(ωt)}ȷˆ].
  3. 3.

    Acceleration:

    𝒂C = ω2r3[sin(ωt3)ıˆ+cos(ωt3)ȷˆ],
    𝒂P = ω2r[{13sin(ωt3)+sin(ωt)}ıˆ+{13cos(ωt3)cos(ωt)}ȷˆ].

    We can now use these expressions to find the position, velocity, and acceleration of the two points at the instants of interest by substituting r=1ft, ω=πrad/s, and appropriate values of t. These values are shown in Table 17.1.

    t 1s 3 s 5.25 s
    𝒓C (ft) 3(32ıˆ12ȷˆ) 3ȷˆ 3(12ıˆ12ȷˆ)
    𝒓P (ft) 𝒓Cȷˆ 𝒓Cȷˆ 4(12ıˆ12ȷˆ)
    𝒗C ( ft/s) π(12ıˆ+32ȷˆ) πıˆ π(12ıˆ12ȷˆ)
    𝒗P ( ft/s) π(12ıˆ+32ȷˆ) 2πıˆ 0
    𝒂C(ft/s2) π23(32ıˆ+12ȷˆ) π23ȷˆ π23(12ıˆ+12ȷˆ)
    𝒗P(ft/s2) 11.86(.24ıˆ+.97ȷˆ) 2π23ȷˆ 13.16(12ıˆ+12ȷˆ)
    Table 17.1: Position, velocity, and acceleration of point P and point C

The velocity and acceleration of the two points are shown in Figures 17.66(a) and (b) respectively.

Filename:sfig4-6-3
Figure 17.66: (a) Velocity and (b) Acceleration of points P and C at t=1s, 3s, and 5.25s.

It is worthwhile to check the directions of velocities and the accelerations by thinking about the velocity and acceleration of point P as a vector sum of the velocity (same for acceleration) of the center of the disk and the velocity (same for acceleration) of point P with respect to the center of the disk. Since the motions involved are circular motions at constant rate, a visual inspection of the velocities and the accelerations is not very difficult. Try it.

SAMPLE 17.15  The rolling disk: path of a point on the disk. For the rolling disk in Sample 17.63, take ω=πrad/s. Draw the path of a point on the rim of the disk for one complete revolution of the center of the disk around the cylinder for the following conditions:

  1. 1.

    R=8r,

  2. 2.

    R=4r, and

  3. 3.

    R=2r.


Solution In Sample 17.63, we obtained a general expression for the position of a point on the disk as a function of time. By computing the position of the point for various values of time t up to the time required to go around the cylinder for one complete cycle, we can draw the path of the point. For the various given conditions, the variable that changes in Eqn. (17.53) is q. We can write a computer program to generate the path of any point on the disk for a given set of R and r. Here is a pseudocode to generate the required path on a computer according to Eqn. (17.53).

A pseudocode to plot the path of a point on the disk:
(pseudo-code) program rollingdisk
%-------------------------------------------------------------
% This code plots the path of any point on a disk of radius
% ’r’ rolling with speed ’w’ inside a cylinder of radius ’R’.
% The point of interest is distance ’l’ away from the center of
% the disk. The coordinates x and y of the specified point P are
% calculated according to the relation mentioned above.
%--------------------------------------------------------------
phi = pi/50*[1,2,3,...,100]   % make a vector phi from 0 to 2*pi
x = R*cos(phi)                % create points on the outer cylinder
y = R*sin(phi)
plot y vs x                   % plot the outer cylinder
hold this plot                % hold to overlay plots of paths
q = r/(R-r)                   % calculate q.
T = 2*pi/(q*w)                % calculate time T for going around-
                              % the cylinder once at speed ’w’.
t = T/100*[1,2,3, ..., 100]   % make a time vector t from 0 to T-
                              % taking 101 points.

rcx = -(R-r)*(sin(q*w*t))     % find the x coordinates of pt. C.
rcy = -(R-r)*(cos(q*w*t))     % find the y coordinates of pt. C.
rpx = rcx-l*sin(w*t)          % find the x coordinates of pt. P.
rpy = rcy + l*cos(q*t)        % find the y coordinates of pt. P.
plot rpy vs rpx               % plot the path of P and the path
plot rcy vs rcx               % of C. For path of C

Once coded, we can use this program to plot the paths of both the center and the point P on the rim of the disk for the three given situations. Note that for any point on the rim of the disk l=r (see Fig 17.65).

  1. 1.

    Let R=4 units. Then r=0.5 for R=8r. To plot the required path, we run our program rollingdisk with desired input,

    R = 4
    r = 0.5
    w = pi
    l = 0.5
    execute rollingdisk

    The plot generated is shown in Fig.17.67 with a few graphic elements added for illustrative purposes.

    Filename:sfig4-6-3a
    Figure 17.67: Path of point P and the center C of the disk for R=8r.
  2. 2.
    Filename:pfigure4-4-rp12
    Figure 17.68: Path of point P and the center C of the disk for R=4r.

    Similarly, for R=4r we type:

    R = 4
    r = 1
    w = pi
    l = 1
    execute rollingdisk

    to plot the desired paths. The plot generated in this case is shown in Fig.17.68

  3. 3.
    Filename:pfigure-blue-38-2
    Figure 17.69: Path of point P and the center C of the disk for R=2r.

    The last one is the most interesting case. The plot obtained in this case by typing:

    R = 4
    r = 2
    w = pi
    l = 2
    execute rollingdisk

    is shown in Fig.17.69. Point P just travels on a straight line! In fact, every point on the rim of the disk goes back and forth on a straight line. Most people find this motion odd at first sight. You can roughly verify the result by cutting a hole twice the diameter of a coin (say a US quarter or dime) in a piece of cardboard and rolling the coin around inside while watching a marked point on the perimeter.

A curiosity. We just discovered something simple about the path of a point on the edge of a circle rolling in another circle that is twice as big. The edge point moves in a straight line. In contrast one might think about the motion of the center G of a straight line segment that slides against two straight walls as in sample 17.89. A problem couldn’t be more different. Naturally the path of point G is a circle (as you can check physically by looking at the middle of a ruler as you hold it sliding against a wall-floor corner).

Problems for 17.3 Kinematics of rolling and sliding

17.3.1   A stone in a wheel. A round wheel rolls to the right. At time t=0 it picks up a stone from the road. The stone is stuck in the edge of the wheel. You want to know the direction of the rock’s motion just before and after it next hits the ground. Here are some candidate answers:

  • When the stone approaches the ground its motion is tangent to the ground.

  • The stone approaches the ground at angle x (you name it).

  • When the stone approaches the ground its motion is perpendicular to the ground.

  • The stone approaches the ground at various angles depending on the following conditions(…you list the conditions.)

Although you could address this question analytically, you are to try to get a clear answer by looking at computer generated plots. In particular, you are to plot the pebble’s path for a small interval of time near when the stone next touches the ground. You should pick the parameters that make your case for an answer the strongest. You may make more than one plot.

Here are some steps to follow:

  1. (a)

    Assuming the wheel has radius Rw and the pebble is a distance Rp from the center (not necessarily equal to Rw). The pebble is directly below the center of the wheel at time t=0. The wheel spins at constant clockwise rate ω. The x-axis is on the ground and x(t=0)=0. The wheel rolls without slipping. Using a clear well labeled drawing (use a compass and ruler or a computer drawing program), show that

    x(t)=ωtRwRpsin(ωt)
    y(t)=RwRpcos(ωt)
  2. (b)

    Using this relation, write a program to make a plot of the path of the pebble as the wheel makes a little more than one revolution. Also show the outline of the wheel and the pebble itself at some intermediate time of interest. [Use any software and computer that pleases you.]

  3. (c)

    Change whatever you need to change to make a good plot of the pebble’s path for a small amount of time as the pebble approaches and leaves the road. Also show the wheel and the pebble at some time in this interval.

  4. (d)

    In this configuration the pebble moves a very small distance in a small time so your axes need to be scaled. But make sure your x- and y- axes have the same scale so that the path of the pebble and the outline of the wheel will not be distorted.

  5. (e)

    How does your computer output buttress your claim that the pebble approaches and leaves the ground at the angles you claim?

  6. (f)

    Think of something about the pebble in the wheel that was not explicitly asked in this problem and explain it using the computer, and/or hand calculation and/or a drawing.

17.3.2  A uniform disk of radius r rolls at a constant rate without slip. A small ball of mass m is attached to the outside edge of the disk. What is the force required to hold the disk in place when the mass is above the center of the disk?

17.3.3  Rolling at constant rate. A round disk rolls on the ground at constant rate. It rolls 114 revolutions over the time of interest.

  1. (a)

    Particle paths. Accurately plot the paths of three points: the center of the disk C, a point on the outer edge that is initially on the ground, and a point that is initially half way between the former two points. [Hint: Write a parametric equation for the position of the points. First find a relation between ω and vC. Then note that the position of a point is the position of the center plus the position of the point relative to the center.] Draw the paths on the computer, make sure x and y scales are the same.

  2. (b)

    Velocity of points. Find the velocity of the points at a few instants in the motion: after 14, 12, 34, and 1 revolution. Draw the velocity vector (by hand) on your plot. Draw the direction accurately and draw the lengths of the vectors in proportion to their magnitude. You can find the velocity by differentiating the position vector or by using relative motion formulas appropriately. Draw the disk at its position after one quarter revolution. Note that the velocity of the points is perpendicular to the line connecting the points to the ground contact.

  3. (c)

    Acceleration of points. Do the same as above but for acceleration. Note that the acceleration of the points is parallel to the line connecting the points to the center of the disk.

Filename:pfigure-blue-49-2
Figure 17.70

17.3.4  The concentric wheels are welded to each other and roll without slip on the rack and stationary support. The rack moves to the right at vr=1m/s. What is the velocity of point A in the middle of the wheels shown?

Filename:pfigure4-4-rp13
Figure 17.71

17.3.5  Questions (a) - (e) refer to the cylinders in the configuration shown in the figure. Question (f) is closely related. Answer the questions in terms of the given quantities (and any other quantities you define if needed).

  1. (a)

    What is the speed (magnitude of velocity) of point c?

  2. (b)

    What is the speed of point P?

  3. (c)

    What is the magnitude of the acceleration of point c?

  4. (d)

    What is the magnitude of the acceleration of point P?

  5. (e)

    What is the radius of curvature of the path of the particle P at the point of interest?

  6. (f)

    In the special case of A=2b what is the curve which particle P traces (for all time)? Sketch the path.

Filename:f92h7p1
Figure 17.72

A little cylinder (with outer radius b and center at point c) rolls without slipping inside a bigger fixed cylinder (with inner radius A and center at point O). The absolute angular velocity of the little cylinder ω is constant. P is attached to the outside edge of the little cylinder. At the instant of interest, P is on the line between O and c.

17.4 Mechanics of contacting bodies: rolling and sliding

A typical machine part has forces that come from contact with other parts. In fact, with the major exceptions of

  • Gravity,

  • Electromagnetic forces inside motors, and

  • Magnetic attraction/repulsion

most of the forces that act on bodies of engineering interest come from contact. Many of the forces you have drawn in free-body diagrams have been contact forces: The force of the ground on a wheel, of an axle on a bearing, of any part on any other part it touches.

We’d now like to consider some mechanics problems that involve sliding or rolling contact. Once you understand the kinematics from the previous section, there is nothing new in the mechanics. As always, the mechanics is linear momentum balance, angular momentum balance and energy balance. Because we are considering single rigid objects in 2D the expressions for the motion quantities are especially simple (as you can look up in Table I at the back of the book):

𝑳˙ =mtot𝒂cm, (17.56)
𝑯˙/C =𝒓cm/C×(mtot𝒂cm)+Iω˙𝒌ˆ (whereI=Izzcm) (17.57)
EK =mtotvcm2/2+Iω2/2 (17.58)

The key to success, as usual, is the drawing of appropriate free-body diagrams (see Chapter 3). For rolling and sliding the two cases one needs to consider as possible are, self-evidently, rolling, where the contact point has no relative velocity and the tangential reaction force is unknown but less than μN, and sliding where the relative velocity could be anything and the tangential reaction force is usually assumed to have a magnitude of μN but oppose the relative motion.

Work-energy relations and impulse-momentum relations are useful to solve some problems both with and without slip. In pure rolling contact the contact force does no work because the material point of contact has no velocity. However, when there is sliding mechanical energy is dissipated. The rate of loss of kinetic and potential energy is

Rate of frictional dissipation=Pdiss=Ffrictionvslip (17.59)

where vslip is the relative velocity of the contacting slipping points. If either the friction force (ideal lubrication) or sliding velocity (no slip) is zero there is no dissipation.

As for various problems throughout the text, it is often a savings of calculation to use angular momentum balance (or moment balance in statics) relative to a point where there are unknown reaction forces. For rolling and slipping problems this often means making use of contact points.

Filename:pfigure-blue-52-2
Figure 17.73: A ball rolls or slides on level ground.

Example: Pure rolling on level ground

A ball or wheel rolling on level ground, with no air friction etc, rolls at constant speed (see fig. 17.73). This is most directly deduced from angular momentum balance about the contact point C:

𝑴C=𝑯˙C𝒓G/C×mgȷˆ=𝒓G/C×m𝒂G+ω˙Izzcm𝒌ˆ𝟎=Rȷˆ×(mω˙Rıˆ)+ω˙Izzcm𝒌ˆdotting with 𝒌ˆω˙=0ω=constant.

Because for rolling vG=ωR we thus have that vG is a constant. [The result can also be obtained by combining angular momentum balance about the center-of-mass with linear momentum balance.]

Finally, linear momentum balance gives the reaction force at C to be 𝑭tot=Fıˆ+nȷˆ=mgȷˆ. So,

assuming point contact, there is no rolling resistance.

Example: Bowling ball with initial sliding

A bowling ball is released with an initial speed of v0 and no rotation rate. What is its subsequent motion? To start with, the motion is incompatible with rolling, the bottom of the ball is sliding to the right. So there is a frictional force which opposes motion and F=μN (see fig. 17.73). Linear and angular momentum balance give:

LMB:{Fıˆ+Nȷˆmgȷˆ=maıˆ}{}ȷˆN=mg{}ıˆa=μgAMB/G:Rμmg=Izzcmω˙v=v0μgt and ω=μRmgt/Izzcm

Thus the forward speed of the ball decreases linearly with time while the counter-clockwise angular velocity decreases linearly with time.

This solution is only appropriate so long as there is rightward slip, vG>ωR. Just like for a sliding block, there is no impetus for reversal, and the ball switches to pure rolling when

v=ωRv0μgt=(μRmgt/Izzcm)Rt=v0μg(1+mR2Izzcm).

Note that the energy lost during sliding is less than μmg times the distance the center of the ball moves during slip.

Filename:pfigure4-4-rp16
Figure 17.74: A ball rolls or slides down a slope.

Example: Ball rolling down hill.

Assuming rolling we can find the acceleration of a ball as it rolls downhill (see fig. 17.74). We start out with the kinematic observations that 𝒂G=aG𝝀ˆ, that Rω=vG and that Rω˙=aG. Angular momentum balance about the stationary point on the ground instantaneously coinciding with the contact point gives

AMB/C𝒓G/C×(mgȷˆ)=𝒓G/C×m𝒂G+Izzcmω˙𝒌ˆ{Rsinϕmg𝒌ˆ=(R𝒏ˆ)×(maG𝝀ˆ)+Izzcmω˙𝒌ˆ}{}𝒌ˆRmgsinϕ=RmaGIzzcmaG/RaG=gsinϕ1+Izzcm/(mR2).

Which is less than the acceleration of a block sliding on a ramp without friction: a=gsinϕ (unless the mass of the rolling ball is concentrated at the center with Izzcm=0). Note that a very small ball rolls just as slowly. In the limit as the ball radius goes to zero the behavior does not approach that of a point mass that slides; the rolling remains significant.

Example: Ball rolling down hill: energy approach

We can find the acceleration of the rolling ball using power balance or conservation of energy. For example

0=ddtET0=EK˙+EP˙=ddt(mv2/2+Izzcmω2/2)+ddt(mgy)=mvv˙+Izzcmωω˙+mgy˙=mvv˙+Izzcm(v/R)v˙/Rmg(sinϕ)vassuming v00=(m+Izzcm/R2)v˙mgsinϕv˙=gsinϕ1+Izzcm/(mR2)

as before.

Example: Does the ball slide?

How big is the coefficient of friction μ needed to prevent slip for a ball rolling down a hill? Use angular momentum balance to find the normal and frictional components of the contact force, using the rolling example above.

LMB (𝑭tot=m𝒂G){N𝒏ˆ+F𝝀ˆmgȷˆ=maG𝝀ˆ}{}𝒏ˆN=mgcosϕ{}𝝀ˆF+mgsinϕ=mgsinϕ1+Izzcm/(mR2)F=mgsinϕ1+mR2/IzzcmCritical condition:μ=|F|N=tanϕ1+mR2/Izzcm

If Izzcm is very small (the mass concentrated near the center of the ball) then small friction is needed to prevent rolling. For a uniform rubber ball on pavement (with μ1 and Izzcm2mR2/5) the steepest slope for rolling without slip is a steep ϕ=tan1(7/2)74. A metal hoop on the other hand (with μ.3 and IzzcmmR2) will only roll without slip for slopes less than about ϕ=tan1(.6)31.

Filename:pfigure4-4-rp14
Figure 17.75: A ball rolls in a round cross-section bowl.

Example: Oscillations of a ball in a bowl.

A round ball can oscillate back and forth in the bottom of a circular cross section bowl or pipe (see fig. 17.75). Similarly, a cylindrical object can roll inside a pipe. What is the period of oscillation? Start with angular momentum balance about the contact point

𝒓G/C×(mgȷˆ)=𝒓G/C×m𝒂G+Izzcmω˙𝒌ˆrmgsinθ𝒌ˆ=r𝒆ˆr×(m((Rr)θ¨𝒆ˆθ(Rr)θ˙2𝒆ˆr))+Izzcmω˙𝒌ˆ.

Evaluating the cross products (using that 𝒆ˆr×𝒆ˆθ=𝒌ˆ) and using the kinematics from the previous section (that (Rr)θ˙=rω) and dotting the left and right sides with 𝒌ˆ gives

(Rr)θ¨=gsinθ1+Izzcm/mr2,

the tangential acceleration is the same as would have been predicted by putting the ball on a constant slope of θ. Using the small angle approximation that sinθ=θ the equation can be rearranged as a standard harmonic oscillator equation

θ¨+(g(Rr)(1+Izzcm/mr2))θ=0.

If all the ball’s mass were concentrated in its middle, keeping r a fixed non-zero size, (so Izzcm=0, like a lead pellet inside a styrofoam ball) this is the same as for a simple pendulum with length Rr. For any parameter values the period of small oscillation is

T=2π(Rr)(1+Izzcm/mr2)g.

For a marble or ball bearing in a sideways glass (with Rr4cm=.04m, Izzcm/mr22/5 and g10m/s2) this gives about one oscillation every half second. See page 19.2 for the energy approach to this problem.

SAMPLE 17.16

Filename:pfigure4-4-rp15
Figure 17.76

Equation of motion of a driven wheel Consider a wheel of mass m and radius R driven by an axle force 𝑭 as shown in the figure. The wheel rolls to the left without slipping. Write the equation of motion of the wheel.


Solution

Filename:tfigure4-spherical-rotaxis
Figure 17.77: FBD of a wheel with mass m. Force F is applied by the axle.

The free-body diagram of the wheel is shown in figure 17.77. We can write the equation of motion of the wheel in terms of either the center-of-mass position x or the angular displacement of the wheel θ. Since in pure rolling, these two variables share a simple relationship (x=Rθ), we can easily get the equation of motion in terms of x if we have the equation in terms of θ and vice versa. Let 𝝎=ω𝒌ˆ and 𝝎˙=ω˙𝒌ˆ.

Since all the forces are shown in the free-body diagram, we can readily write the angular momentum balance for the wheel. We choose the point of contact C as our reference point for the angular momentum balance (because the gravity force, mgȷˆ, the friction force Ffrictionıˆ, and the normal reaction of the ground Nȷˆ, all pass through the contact point C and therefore, produce no moment about this point). We have

𝑴C = 𝑯˙/C
where
𝑴C = 𝒓cm/CRȷˆ×(F𝝀ˆ)
= Rȷˆ×F(cosϕıˆsinϕȷˆ)
= FRcosϕ𝒌ˆ
and
𝑯˙/C = 𝒓cm/C×m𝒂cm+Izzcm𝝎˙
= Rȷˆ×mx¨ω˙Rıˆ+Izzcmω˙𝒌ˆ
= mω˙R2𝒌ˆ+Izzcmω˙𝒌ˆ
= (Izzcm+mR2)ω˙𝒌ˆ.
Thus,
FRcosϕ𝒌ˆ = (Izzcm+mR2)ω˙𝒌ˆ
 ω˙θ¨ = FRcosϕIzzcm+mR2

which is the equation of motion we are looking for.

Note that we can easily substitute θ¨=x¨/R in the equation of motion above to get the equation of motion in terms of the center-of-mass displacement x as

x¨=FR2cosϕIzzcm+mR2.

Answer: θ¨=FRcosϕIzzcm+mR2

Comments: We could have, of course, used linear momentum balance with angular momentum balance about the center-of-mass to derive the equation of motion. Note, however, that the linear momentum balance will essentially give two scalar equations in the x and y directions involving all forces shown in the free-body diagram . The angular momentum balance , on the other hand, gets rid of some of them. Depending on which forces are known, we may or may not need to use all the three scalar equations. In the final equation of motion, we must have only one unknown.

SAMPLE 17.17

Filename:tfigure5-7
Figure 17.78

Energy and power of a rolling wheel. A wheel of diameter 2ft and mass 20lbm rolls without slipping on a horizontal surface. The kinetic energy of the wheel is 1700ftlbf. Assume the wheel to be a thin, uniform disk.

  1. 1.

    Find the rate of rotation of the wheel.

  2. 2.

    Find the average power required to bring the wheel to a complete stop in 5 s.


Solution

  1. 1.

    Let ω be the rate of rotation of the wheel. Since the wheel rotates without slip, its center-of-mass moves with speed vcm=ωr. The wheel has both translational and rotational kinetic energy. The total kinetic energy is

    EK = 12mvcm2+12Icmω2
    = 12mω2r2+12Icmω2
    = 12(mr2+Icm12mr2)ω2
    = 34mr2ω2
     ω2 = 4EK3mr2
    = 4×1700ftlbf3×20lbm1ft2
    = 4×1700×32.2lbmft/s23×20lbmft
    = 3649.331s2
     ω = 60.4rad/s.

    Answer: ω=60.4rad/s

    Note: This rotational speed, by the way, is extremely high. At this speed the center-of-mass moves at 60.4  ft/s!

  2. 2.

    Power is the rate of work done on a body or the rate of change of kinetic energy. Here we are given the initial kinetic energy, the final kinetic energy (zero) and the time to achieve the final state. Therefore, the average power is,

    P = EK1EK2Δt
    = 1700ftlbf05s=340ftlbf/s
    = 340ftlbf/s1 hp550ftlbf/s
    = 0.62 hp

    Answer: P=0.62 hp

SAMPLE 17.18

Filename:tfigure5-gen-rigid-body
Figure 17.79: FBD of a rolling wheel.

Equation of motion of a rolling wheel from energy balance. Consider the wheel with mass m from figure 17.79. The free-body diagram of the wheel is shown here again. Derive the equation of motion of the wheel using energy balance. Assume CCW (counter clockwise) rotation is positive.


Solution

From energy balance, we have

P = E˙K
where
P = 𝑭i𝒗i
= Ffrictionıˆ𝒗C𝟎+Nȷˆ𝒗C𝟎mgȷˆ𝒗cmvıˆ+F𝝀ˆ𝒗cmvıˆ
= mgv(ıˆȷˆ0)+Fv(𝝀ˆıˆcosϕ)
= Fvcosϕ
and
E˙K = ddt(12mx˙2+12Izzcmω2(x˙/R)2)
= 12ddt[(m+IzzcmR2)x˙2]
= (m+IzzcmR2)x˙x¨.
Thus,
Fvcosϕ = (m+IzzcmR2)x˙x¨
or Fx˙cosϕ = (m+IzzcmR2)x˙x¨
 x¨ = Fcosϕm+IzzcmR2.

We can also write the equation of motion in terms of θ by replacing x¨ with θ¨R giving,

θ¨=(F/R)cosϕm+IzzcmR2.

Answer: x¨=Fcosϕm+Izzcm/R2

Comments: In the equations above (for calculating P), we have set 𝒗C=𝟎 because in pure rolling, the instantaneous velocity of the contact point is zero. Note that the force due to gravity is normal to the direction of the velocity of the center-of-mass. So, the only power supplied to the wheel is due to the force F𝝀ˆ acting at the center-of-mass.

SAMPLE 17.19

Filename:tfigure5-term1-a
Figure 17.80

Equation of motion of a rolling disk on an incline. A uniform circular disk of mass m=1kg and radius R=0.4m rolls down an inclined shown in the figure. Write the equation of motion of the disk assuming pure rolling, and find the distance travelled by the center-of-mass in 2s.


Solution

The free-body diagram of the disk is shown in fig. 17.81. In addition to the base unit vectors ıˆ and ȷˆ, let us use unit vectors 𝝀ˆ and 𝒏ˆ along the plane and perpendicular to the plane, respectively, to express various vectors.

Filename:tfigure5-term1-b
Figure 17.81

We can write the equation of motion using linear momentum balance or angular momentum balance. However, note that if we use linear momentum balance we have two unknown forces in the equation. On the other hand, if we use angular momentum balance about the contact point C, these forces do not show up in the equation. So, let us use angular momentum balance about point C:

𝑴C = 𝑯˙/C
where
𝑴C = 𝒓O/C×m𝒈=R𝒏ˆ×(mgȷˆ)
= Rmgsinα𝒌ˆ
and
𝑯˙/C = Izzcmω˙𝒌ˆ+𝒓O/CR𝒏ˆ×m𝒂cmRω˙𝝀ˆ
= Izzcmω˙𝒌ˆ+mR2ω˙(𝒏ˆ×𝝀ˆ)
= (Izzcm+mR2)ω˙𝒌ˆ.
Thus,
Rmgsinα𝒌ˆ = (Izzcm+mR2)ω˙𝒌ˆ
 ω˙ = gsinαR[1+Izzcm/(mR2)].

Answer: ω˙=gsinαR[1+Izzcm/(mR2)]

Note that in the above equation of motion, the right hand side is constant. So, we can solve the equation for ω and θ by simply integrating this equation and substituting the initial conditions ω(t=0)=0 and θ(t=0)=0. Let us write the equation of motion as ω˙=β where β=gsinα/R(1+Izzcm/mR2). Then,

ω θ˙=βt+C1
θ = 12βt2+C1t+C2.

Substituting the given initial conditions θ˙(0)=0 and θ(0)=0, we get C1=0 and C2=0, which implies that θ=12βt2. Now, in pure rolling, x=Rθ. Therefore,

x(t) = Rθ(t)=12βt2=R12gsinαR(1+Izzcm/mR2)t2
= 12gsinα1+12mR2mR2t2=13(gsinα)t2
x(2s) = 139.8m/s2sin(30)(2s)2=6.53m.

Answer: x(2s)=6.53m

SAMPLE 17.20

Filename:sfig1-2-12
Figure 17.82

Using Work and energy in pure rolling. Consider the disk of Sample 17.79 rolling down the incline again. Suppose the disk starts rolling from rest. Find the speed of the center-of-mass when the disk is 2m down the inclined plane.


Solution We are given that the disk rolls down, starting with zero initial velocity. We are to find the speed of the center-of-mass after it has travelled 2malong the incline. We can, of course, solve this problem using equation of motion, by first solving for the time t the disk takes to travel the given distance and then evaluating the expression for speed ω(t) or x(t) at that t. However, it is usually easier to use work energy principle whenever positions are specified at two instants, speed is specified at one of those instants, and speed is to be found at the other instant. This is because we can, presumably, compute the work done on the system in travelling the specified distance and relate it to the change in kinetic energy of the system between the two instants.

Filename:sfig4-6-5
Figure 17.83

In the problem given here, let ω1 and ω2 be the initial and final (after rolling down by d=2m) angular speeds of the disk, respectively. We know that in rolling, the kinetic energy is given by

EK=12mvcm2(ωR)2+12Izzcmω2=12(mR2+Izzcm)ω2.

Therefore,

ΔEK=EK2EK1=12(mR2+Izzcm)(ω22ω12). (17.60)

Now, let us calculate the work done by all the forces acting on the disk during the displacement of the mass-center by d along the plane. Note that in ideal rolling, the contact forces do no work. Therefore, the work done on the disk is only due to the gravitational force:

W=(mgȷˆ)(d𝝀ˆ)=mgd(ȷˆ𝝀ˆsinα)=mgdsinα. (17.61)

From work-energy principle (integral form of power balance, P=E˙K ), we know that W=ΔEK. Therefore, from eqn. (17.60) and eqn. (17.61), we get

mgdsinα = 12(mR2+Izzcm)(ω22ω12)
 ω22 = ω12+2mgdsinαmR2+Izzcm=ω12+2gdsinαR2(1+IzzcmmR2)
= ω12+4gdsinα3R2.
Substituting the values of g,d,α,R, etc., and setting ω1=0, we get
ω22 = 4(9.8m/s2)(2m)(sin(30)3(0.4m)2=81.67/s2
 ω2 = 9.04rad/s.

The corresponding speed of the center-of-mass is

vcm=ω2R=9.04rad/s0.4m=3.61m/s.

Answer: vcm=3.61m/s

SAMPLE 17.21

Filename:sfig4-6-5a
Figure 17.84

Impulse and momentum calculations in pure rolling. Consider the disk of Sample 17.79 rolling down the incline again. Find an expression for the rolling speed (ω) of the disk after a finite time Δt, given the rolling speed ω1 at some instant t1.


Solution

Once again, this problem can be solved by integrating the equation of motion (as done in Sample 17.79). However, we will solve this problem here using impulse-momentum relationship. Let the desired angular speed at t2=t1+Δt be ω2. We need to find ω2, given the ω1 at t=t1. Since the forces acting on the disk do not change during this time (assuming pure rolling), it is easy to calculate impulse and then relate it to the change in the momenta of the disk between the two instants.

Filename:sfig4-6-5b
Figure 17.85

Now, from the linear impulse momentum relationship, 𝑭Δt=𝑳2𝑳1, we have

(F𝝀ˆ+Nıˆmgȷˆ)Δt = m(v2v1)𝝀ˆ. (17.62)
Dotting eqn. (17.62) with 𝝀ˆ gives
(Fmg(ȷˆ𝝀ˆsinα))Δt = m(v2v1)
(F+mgsinα)Δt = mR(ω2ω1). (17.63)

Similarly, the angular impulse-momentum relationship about the mass-center, 𝑴OΔt=(𝑯/O)2(𝑯/O)1, gives

(FR𝒌ˆ)Δt = Izzcm(ω2ω1)𝒌ˆ
 FRΔt = Izzcm(ω2ω1). (17.64)

Note that the other forces (N and mg) do not produce any moment about the mass-center as they pass through this point. We can now eliminate the unknown force F from eqn. (17.63) and eqn. (17.64) by multiplying eqn. (17.63) with R and adding to eqn. (17.64):

(F+mgsinα)ΔtR+FRΔt = mR(ω2ω1)R+Izzcm(ω2ω1)
or mgRsinαΔt = (Izzcm+mR2)(ω2ω1)
or gsinαΔt = R(1+IzzcmmR2)(ω2ω1)
 ω2 = ω1+gsinαR(1+IzzcmmR2)Δt.

Answer: ω2=ω1+gsinαR(1+IzzcmmR2)Δt


Note: From the answer obtained, we clearly see that ω2=ω1+ω˙Δt where ω˙ is the angular acceleration given by the expression ω˙=gsinαR(1+Izzcm/mR2). This is the same expression we obtained for the angular acceleration in Sample 17.79.

SAMPLE 17.22

Filename:sfig4-6-5c
Figure 17.86: A ladder, modeled as a uniform rod of mass m and length , falls from a rest position at θ=θo(<π/2) such that its ends slide along frictionless vertical and horizontal surfaces.

Falling ladder. A ladder AB, modeled as a uniform rigid rod of mass m and length , rests against frictionless horizontal and vertical surfaces. The ladder is released from rest at θ=θo(θo<π/2). Assume the motion to be planar (in the vertical plane).

  1. 1.

    As the ladder falls, what is the path of the center-of-mass of the ladder?

  2. 2.

    Find the equation of motion (e.g., a differential equation in terms of θ and its time derivatives) for the ladder.

  3. 3.

    How does the angular speed ω(=θ˙) depend on θ?


Solution Since the ladder is modeled by a uniform rod AB, its center-of-mass is at G, half way between the two ends. As the ladder slides down, the end A moves down along the vertical wall and the end B moves out along the floor. Note that it is a single degree of freedom system as angle θ (a single variable) is sufficient to determine the position of every point on the ladder at any instant of time.

  1. 1.

    Path of the center-of-mass:

    Filename:sfig4-6-5d
    Figure 17.87

    Let the origin of our x-y coordinate system be the intersection of the two surfaces on which the ends of the ladder slide (see Fig. 17.87). The position vector of the center-of-mass G may be written as

    𝒓G = 𝒓B+𝒓G/B (17.65)
    = cosθıˆ+2(cosθıˆ+sinθȷˆ)
    = 2(cosθıˆ+sinθȷˆ).

    Thus the coordinates of the center-of-mass are

    xG=2cosθ and yG=2sinθ,

    from which we get

    xG2+yG2=24

    which is the equation of a circle of radius 2. Therefore, the center-of-mass of the ladder follows a circular path of radius 2 centered at the origin. Of course, the center-of-mass traverses only that part of the circle which lies between its initial position at θ=θo and the final position at θ=0.

  2. 2.

    Equation of motion:

    Filename:sfig4-5-6
    Figure 17.88: The free-body diagram of the ladder.

    The free-body diagram of the ladder is shown in Fig. 17.88. Since there is no friction, the only forces acting at the end points A and B are the normal reactions from the contacting surfaces.

    Now, writing the linear momentum balance (𝑭=m𝒂) for the ladder we get

    N1ıˆ+(N2mg)ȷˆ=m𝒂G=m𝒓¨G.

    Differentiating eqn. (17.65) twice we get 𝒓¨G as

    𝒓¨G=2[(θ¨sinθθ˙2cosθ)ıˆ+(θ¨cosθθ˙2sinθ)ȷˆ].

    Substituting this expression in the linear momentum balance equation above and dotting both sides of the equation by ıˆ and then by ȷˆ we get

    N1 = 12m(θ¨sinθ+θ˙2cosθ)
    N2 = 12m(θ¨cosθθ˙2sinθ)+mg.

    Next, we write the angular momentum balance for the ladder about its center-of-mass, 𝑴/G=𝑯˙/G, where

    𝑴/G = (N12sinθ+N22cosθ)𝒌ˆ
    = 12m(θ¨sinθ+θ˙2cosθ)2sinθ𝒌ˆ
    +[12m(θ¨cosθθ˙2sinθ)+mg]2cosθ𝒌ˆ
    = (14m2θ¨+12mgcosθ)𝒌ˆ
    and
    𝑯˙/G = Izz/G𝝎˙=112m2θ¨(𝒌ˆ),

    where 𝝎˙=θ¨(𝒌ˆ) because θ is measured positive in the clockwise direction (𝒌ˆ). Now, equating the two quantities 𝑴/G=𝑯˙/G and dotting both sides with 𝒌ˆ we get

    14m2θ¨+12mgcosθ = 112m2θ¨
    or (112+14)2θ¨ = 12gcosθ
    or θ¨ = 3g2cosθ (17.66)

    which is the required equation of motion. Unfortunately, it is a nonlinear equation which does not have a nice closed form solution for θ(t).

  3. 3.

    Angular Speed of the ladder: To solve for the angular speed ω(=θ˙) as a function of θ we need to express eqn. (17.66) in terms of ω,θ, and derivatives of ω with respect to θ. Now,

    θ¨=ω˙=dωdt=dωdθdθdt=ωdωdθ.

    Substituting in eqn. (17.66) and integrating both sides from the initial rest position to an arbitrary position θ we get

    0ωω𝑑ω = θ0θ3g2cosθdθ
     12ω2 = 3g2(sinθsinθ0)
     ω = ±3g(sinθ0sinθ).

    Since end B is sliding to the right, θ is decreasing; hence it is the negative sign in front of the square root which gives the correct answer, i.e.,

    𝝎=θ˙(𝒌ˆ)=3g(sinθ0sinθ)𝒌ˆ.

Note: To do this problem we have assumed that the upper end of the ladder stays in contact with the wall as it slides down. If a real ladder were sliding against a slippery wall and floor, would it not lose contact? The answer is yes, it would. One way of finding when the contact would be lost is to calculate the normal reaction N1 and finding out at what value of θ it passes through zero. The answer depends on θ0. For example, if θ0=80, then N1 is zero at about θ=41. You can verify this result by substituting the expressions for θ¨ and θ˙ in the expression for N1 and solving for θ when N1=0.

SAMPLE 17.23

Filename:sfig4-5-6a
Figure 17.89

The falling ladder again. Consider the falling ladder of Sample 17.58 again. The mass of the ladder is m and the length is . The ladder is released from rest at θ=80.

  1. 1.

    At the instant when θ=45, find the speed of the center-of-mass of the ladder using energy.

  2. 2.

    Derive the equation of motion of the ladder using work-energy balance.


Solution

Filename:sfig4-6-8
Figure 17.90
  1. 1.

    Since there is no friction, there is no loss of energy between the two states: θ0=80 and θf=45. The only external forces on the ladder are N1,N2,andmg as shown in the free-body diagram. Since the displacements of points A and B are perpendicular to the normal reactions of the walls, N1 and N2, respectively, no work is done by these forces on the ladder. The only force that does work is the force due to gravity. But this force is conservative. Therefore, the conservation of energy holds between any two states of the ladder during its fall.

    Filename:sfig4-6-8a
    Figure 17.91

    Let E1 and E2 be the total energy of the ladder at θ0 and θf, respectively. Then

    E1=E2(conservation of energy).
    NowE1 = EK1K.E.+EP1P.E.
    = 0+mgh1
    = mg2sinθ0
    andE2 = EK2+EP2=12mvG2+12IzzGω2EK2+mgh2.

    Equating E1 and E2 we get

    mg2(sinθ0sinθf) = 12(mvG2+112m2IzzGω2)
    org(sinθ0sinθf) = vG2+1122ω2. (17.67)

    Clearly, we cannot find vG from this equation alone because the equation contains another unknown, ω. So we need to find another equation which relates vG and ω. To find this equation we turn to kinematics. Note that

    𝒓G = 2(cosθıˆ+sinθȷˆ)
     𝒗G = 𝒓˙G=2(sinθθ˙ıˆ+cosθθ˙ȷˆ)
     vG = |𝒗G|=24(cos2θ+sin2θ)θ˙2
    = 2θ˙=2ω
     ω = 2vG.

    Substituting the expression for ω in eqn. (17.67) we get

    g(sinθ0sinθf) = vG2+11224vG22
    = 43vG2
     vG = 3g4(sinθ0sinθf)
    = 0.46g.

    Answer: vG=0.46g

  2. 2.

    Equation of motion: Since the ladder is a single degree of freedom system, we can use the power equation to derive the equation of motion:

    P=E˙K.

    For the ladder, the only force that does work is mg. This force acts at the center-of-mass G. Therefore,

    P = 𝑭𝒗=mgȷˆ𝒗G
    = mgȷˆ[2(sinθıˆ+cosθȷˆ)θ˙]
    = mg2θ˙cosθ.
      Now, the rate of change of kinetic energy is
    E˙K = ddt(12mvG2+12IzzGω2)
    = ddt(12m2ω24+12m212ω2)
    = m24ωω˙+m212ωω˙
    = m23ωω˙m23θ˙θ¨(since ω=θ˙ and ω˙=θ¨).

    Now equating P and E˙K we get

    m23θ˙θ¨ = mg2θ˙cosθ
     θ¨ = 3g2cosθ

    which is the same expression as obtained in Sample 17.85 (b).

    Answer: θ¨=3g2cosθ

Note: Just like in the last problem, we have assumed that the upper end of the ladder stays in contact with the vertical wall. This assumption is used implicitly in the kinematics of the ladder. In energy calculations, the normal reactions of the two walls make no contribution as the work done by them is zero and hence the assumption of wall contact does not figure in the calculations above.

SAMPLE 17.24

Filename:sfig4-6-8b
Figure 17.92: A composite wheel made of three uniform disks rolls down an inclined wedge without slipping.

Rolling on an inclined plane. A wheel is made up of three uniform disks— the center disk of mass m=1kg, radius r=10cm and two identical outer disks of mass M=2kg each and radius R. The wheel rolls down an inclined wedge without slipping. The angle of inclination of the wedge with horizontal is θ=30. The radius of the bigger disks is to be selected such that the linear acceleration of the wheel center does not exceed 0.2g. Find the radius R of the bigger disks.


Solution

Since a bound is prescribed on the linear acceleration of the wheel and the radius of the bigger disks is to be selected to satisfy this bound, we need to find an expression for the acceleration of the wheel (hopefully) in terms of the radius R.

The free-body diagram of the wheel is shown in Fig. 17.93. In addition to the weight (m+2M)g of the wheel and the normal reaction N of the wedge surface there is an unknown force of friction Ff acting on the wheel at point C. This friction force is necessary for the condition of rolling motion. You must realize, however, that FfμN because there is neither slipping nor a condition of impending slipping. Thus the magnitude of Ff is not known yet.

Filename:sfig5-5-2
Figure 17.93: Free-body diagram of the wheel.

Let the acceleration of the center-of-mass of the wheel be

𝒂G=aG𝝀ˆ

and the angular acceleration of the wheel be

𝝎˙=ω˙𝒌ˆ.

We assumed 𝝎˙ to be in the negative 𝒌ˆ direction. But, if this assumption is wrong, we will get a negative value for ω˙.

Now we write the equation of linear momentum balance for the wheel:

𝑭 = mtotal𝒂cm
(m+2M)gȷˆ+N𝒏ˆFf𝝀ˆ = (m+2M)aG𝝀ˆ

This 2-D vector equation gives (at the most) two independent scalar equations. But we have three unknowns: N,Ff, and aG. Thus we do not have enough equations to solve for the unknowns including the quantity of interest, aG. So, we now write the equation of angular momentum balance for the wheel about the point of contact C (using 𝒓G/C=r𝒏ˆ):

Filename:sfig5-5-2a
Figure 17.94: Geometry of unit vectors. This diagram can be used to find various dot and cross products between any two unit vectors. For example, 𝒏ˆ×ȷˆ=sinθ𝒌ˆ.
𝑴C = 𝑯˙/C
where
𝑴C = 𝒓G/C×(m+2M)g(ȷˆ)
= r𝒏ˆ×(m+2M)g(ȷˆ)
= (m+2M)grsinθ𝒌ˆ(see Fig. 17.94)
and
𝑯˙/C = IzzG𝝎˙+𝒓G/C×mtotal𝒂G
= IzzG(ω˙𝒌ˆ)mtotalω˙r2𝒌ˆ
= (IzzG+mtotalr2)(ω˙𝒌ˆ)
= [(12mr2+212MR2)+(m+2M)mtotalr2](ω˙𝒌ˆ)
= [32mr2+M(R2+2r2)]ω˙𝒌ˆ.

Thus,

(m+2M)grsinθ𝒌ˆ = [32mr2+M(R2+2r2)]ω˙𝒌ˆ
 ω˙ = (m+2M)grsinθ32mr2+M(R2+2r2). (17.68)

Now we need to relate ω˙ to aG. From the kinematics of rolling,

aG=ω˙r.

Therefore, from Eqn. (17.68) we get

aG = (m+2M)gr2sinθ32mr2+M(R2+2r2).

Now we can solve for R in terms of aG:

32mr2+M(R2+2r2) = (m+2M)gr2sinθaG
 M(R2+2r2) = (m+2M)gaGr2sinθ32mr2
 R2 = (m+2M)gMaGr2sinθ3m2Mr22r2.

Since we require aG0.2g we get

R2 ((m+2M)gM0.2gsinθ3m2M2)r2
(5kg0.4kg123kg4kg2)(0.1m)2
0.035m2
 R 0.187m.

Thus the outer disks of radius 20cm will do the job.

Answer: R18.7cm

SAMPLE 17.25

Filename:sfig5-5-2b
Figure 17.95

Which one starts rolling first — a marble or a bowling ball? A marble and a bowling ball, made of the same material, are launched on a horizontal platform with the same initial velocity, say v0. The initial velocity is large enough so that both start out sliding. Towards the end of their motion, both have pure rolling motion. If the radius of the bowling ball is 16 times that of the marble, find the instant, for each ball, when the sliding motion changes to rolling motion.


Solution

Let us consider one ball, say the bowling ball, first. Let the radius of the ball be r and mass m. The ball starts with center-of-mass velocity 𝒗o=v0ıˆ. The ball starts out sliding. During the sliding motion, the force of friction acting on the ball must equal μN (see the FBD). The friction force creates a torque about the mass-center which, in turn, starts the rolling motion of the ball. However, rolling and sliding coexist for a while, till the speed of the mass-center slows down enough to satisfy the pure rolling condition, v=ωr. Let the instant of transition from the mixed motion to pure rolling be t.

Filename:sfig7-4-2
Figure 17.96: Free-body diagram of the ball during sliding.

From linear momentum balance , we have

mv˙ıˆ = μNıˆ+(Nmg)ȷˆ (17.69)
eqn. (17.69)ȷˆ  N = mg
eqn. (17.69)ıˆ  mv˙ = μN=μmg
 v˙ = μg
 v = v0μgt. (17.70)

Similarly, from angular momentum balance about the mass-center, we get

Izzcmω˙𝒌ˆ = μNr𝒌ˆ=μmgr𝒌ˆ
 ω˙ = μmgrIzzcm
 ω = ω00+μmgrIzzcmt. (17.71)

At the instant of transition from mixed rolling and sliding to pure rolling, i.e., at t=t, v=ωr. Therefore, from eqn. (17.70) and eqn. (17.71), we get

v0μgt = μmgr2Izzcmt
 v0 = μgt(1+mr2Izzcm)
 t = v0μg(1+mr2Izzcm).

Now, for a sphere, Izzcm=25mr2. Therefore,

t=v0μg(1+mr225mr2)=2v07μg.

Note that the expression for t is independent of mass and radius of the ball! Therefore, the bowling ball and the marble are going to change their mixed motion to pure rolling at exactly the same instant. This is not an intuitive result.

Answer: t=2v07μg for both.

SAMPLE 17.26

Filename:pfigure-blue-99-1
Figure 17.97

Transition from a mix of sliding and rolling to pure rolling, using impulse-momentum. Consider the problem in Sample 17.94 again: A ball of radius r=10cm and mass m=1kg is launched horizontally with initial velocity v0=5m/s on a surface with coefficient of friction μ=0.2. The ball starts sliding, rolls and slides simultaneously for a while, and then rolls without sliding. Find the time it takes to start pure rolling.


Solution

Let us denote the time of transition from mixed motion (rolling and sliding) to pure rolling by t. At t=0, we know that vcm=v0=5m/s, and ω0=0. We also know that at t=t, vcm=vt=ωtr, where r is the radius of the ball. We do not know t and vt. However, we are considering a finite time event (during t) and the forces acting on the ball during this duration are known. Recall that impulse momentum equations involve the net force on the body, the time of impulse, and momenta of the body at the two instants. Momenta calculations involve velocities. Therefore, we should be able to use impulse-momentum equations here and find the desired unknowns.

Filename:pfigure-blue-47-2
Figure 17.98: Free-body diagram of the ball during sliding.

From linear impulse-momentum, we have

(𝑭)t = mvtıˆmv0ıˆ
(μNıˆ+(Nmg)ȷˆ)t = m(vtv0)ıˆ.
Dotting the above equation with ȷˆ and ıˆ, respectively, we get
N = mg
μNmgt = m(vtv0)
 μgt = vtv0. (17.72)

Similarly, from angular impulse-momentum relation about the mass-center, we get

𝑴cmt = (𝑯cm)t(𝑯cm)0
(μNr𝒌ˆ)t = (IzzcmωtIzzcmω00)(𝒌ˆ)
or μmgrt = Izzcmωt
 ωt = μmgrt/Izzcm
 vt ωtr=μmgr2t/Izzcm.

Substituting this expression for vt in eqn. (17.72), we get

μgt = μmgr2t/Izzcmv0
 t = v0μg(1+mr2Izzcm)

which is, of course, the same expression we obtained for t in Sample 17.94. Again, noting that Izzcm=25mr2 for a sphere, we calculate the time of transition as

t=2v07μg=2(5m/s)7(0.2)(9.8m/s2)=0.73s.

Answer: t=0.73s

Problems for 17.4 Mechanics of contact

17.4.1  A uniform disc of mass m and radius r rolls without slip at constant rate. What is the total kinetic energy of the disk?

17.4.2  A non-uniform disc of mass m and radius r rolls without slip at constant rate. The center of mass is located at a distance r2 from the center of the disc. What is the total kinetic energy of the disc when the center of mass is directly above the center of the disc?

17.4.3  Falling hoop. A bicycle rim (no spokes, tube, tire, or hub) is idealized as a hoop with mass m and radius R. G is at the center of the hoop. An inextensible string is wrapped around the hoop and attached to the ceiling. The hoop is released from rest at the position shown at t=0.

  1. (a)

    Find yG at a later time t in terms of any or all of m, R, g, and t.

  2. (b)

    Does G move sideways as the hoop falls and unrolls?

Filename:pfigure-blue-89-1
Figure 17.99

17.4.4  A uniform disk with radius R and mass m has a string wrapped around it. The string is pulled with a force F. The disk rolls without slipping.

  1. (a)

    What is the angular acceleration of the disk, 𝜶Disk=θ¨𝒌ˆ? Make any reasonable assumptions you need that are consistent with the figure information and the laws of mechanics. State your assumptions. Answer: αDisk=43rad/s2kˆ.

  2. (b)

    Find the acceleration of the point A in the figure. Answer: aA=83m/s2ıˆ.

Filename:pfigure-blue-35-2
Figure 17.100

17.4.5  If a pebble is stuck to the edge of the wheel in problem 17.3, what is the maximum speed of the pebble during the motion? When is the force on the pebble from the wheel maximum? Draw a good FBD including the force due to gravity.

Filename:pfigure4-3Dpend
Figure 17.101

17.4.6  Spool Rolling without Slip and Pulled by a Cord. The light-weight spool is nearly empty but a lead ball with mass m has been placed at its center. A force F is applied in the horizontal direction to the cord wound around the wheel. Dimensions are as marked. Coordinate directions are as marked.

  1. (a)

    What is the acceleration of the center of the spool? Answer: aC=(F/m)(1Ri/RO)ıˆ.

  2. (b)

    What is the horizontal force of the ground on the spool?

    Answer: (Ri/RO)Fıˆ.

Filename:pfigure-s94h6p3
Figure 17.102

17.4.7  A film spool is placed on a very slippery table. Assume that the film and reel (together) have mass distributed the same as a uniform disk of radius Ri. What, in terms of Ri,Ro,m,g,ıˆ,ȷˆ, and F are the accelerations of points C and B at the instant shown (the start of motion)?

Filename:pfigure-s94h6p4
Figure 17.103

17.4.8  Again, Spool Rolling without Slip and Pulled by a Cord. Reconsider the spool from problem 17.4. This time, a force F is applied vertically to the cord wound around the wheel. In this case, what is the acceleration of the center of the spool? Is it possible to pull the cord at some angle between horizontal and vertical so that the angular acceleration of the spool or the acceleration of the center of mass is zero? If so, find the angle in terms of Ri, R0, m, and F.

17.4.9  A napkin ring lies on a black velvet tablecloth. The thin ring (of mass m, radius r) rolls without slip as a mischievous child pulls the tablecloth (mass M) out with acceleration 𝒂A. The ring starts at the right end (x=d). You can make a reasonable physical model of this situation with an empty soda can and a piece of paper on a flat table.

  1. (a)

    What is the ring’s acceleration as the tablecloth is being withdrawn?

  2. (b)

    How far has the tablecloth moved to the right from its starting point x=0 when the ring rolls off its left-hand end?

  3. (c)

    Clearly describe the subsequent motion of the ring. Which way does it end up rolling at what speed?

  4. (d)

    Would your answer to the previous question be different if the ring slipped on the cloth as the cloth was being pulled out?

Filename:pg84-3
Figure 17.104

17.4.10  A block of mass M is supported by two rollers (uniform cylinders) each of mass m and radius r. They roll without slip on the block and the ground. A force F is applied in the horizontal direction to the right, as shown in the figure. Given F, m, r, and M, find:

  1. (a)

    the acceleration of the block,

  2. (b)

    the acceleration of the center of mass of this block/roller system,

  3. (c)

    the reaction at the wheel bases,

  4. (d)

    the force of the right wheel on the block,

  5. (e)

    the acceleration of the wheel centers, and

  6. (f)

    the angular acceleration of the wheels.

Filename:pfigure-blue-36-1
Figure 17.105

17.4.11  Dropped spinning disk. 2-D . A uniform disk of radius R and mass m is gently dropped onto a surface and doesn’t bounce. When it is released it is spinning clockwise at the rate θ˙0. The disk skids for a while and then is eventually rolling.

  1. (a)

    What is the speed of the center of the disk when the disk is eventually rolling (answer in terms of g, μ, R, θ˙0, and m)? Answer: speed v=Rθ˙3.

  2. (b)

    In the transition from slipping to rolling, energy is lost to friction. How does the amount lost depend on the coefficient of friction (and other parameters)? How does this loss make or not make sense in the limit as μ0 and the dissipation rate zero? Answer: The energy lost to friction is Efric=mR2θ˙026. The energy lost to friction is independent of μ for μ>0. Thus, the energy lost to friction is constant for given m, R, and θ˙0. As μ0, the transition time to rolling . It is not true, however, that the energy lost to friction 0 as μ0. Since the energy lost is constant for any μ>0, the disk will slip for longer and longer times so that the distance of slip goes to infinity. The dissipation rate 0 since the constant energy is divided by increasing transition time. The energy lost is zero only for μ=0.

Filename:pfigure-s95q8
Figure 17.106

17.4.12  Disk on a conveyor belt. A uniform metal cylinder with mass of 200kg is carried on a conveyer belt which moves at V0=3m/s. The disk is not rotating when on the belt. The disk is delivered to a flat hard platform where it slides for a while and ends up rolling. How fast is it moving (i.e. what is the speed of the center of mass) when it eventually rolls? Answer: V=2m/s.

Filename:pfigure4-rpi
Figure 17.107

17.4.13  A rigid hoop with radius R and mass m is rolling without slip so that its center has translational speed vo. It then hits a narrow bar with height R/2. When the hoop hits the bar suddenly it sticks and doesn’t slide. It does hinge freely about the bar, however. The gravitational constant is g. How big is vo if the hoop just barely rolls over the bar? Answer: vo=22gR3.

Filename:pfigure-blue-86-1
Figure 17.108

17.4.14  2-D rolling of an unbalanced wheel. A wheel, modeled as massless, has a point mass (mass =m) at its perimeter. The wheel is released from rest at the position shown. The wheel makes contact with coefficient of friction μ.

  1. (a)

    What is the acceleration of the point P just after the wheel is released if μ=0?

  2. (b)

    What is the acceleration of the point P just after the wheel is released if μ=2?

  3. (c)

    Assuming the wheel rolls without slip (no-slip requires, by the way, that the friction be high: μ=) what is the velocity of the point P just before it touches the ground?

Filename:pfigure-blue-88-1
Figure 17.109

17.4.15  Spool and mass. A reel of mass M and moment of inertia Izzcm=I rolls without slipping upwards on an incline with slope-angle α. It is pulled up by a string attached to mass m as shown. Find the acceleration of point G in terms of some or all of M,m,I,R,r,α,g and any base vectors you clearly define.

Filename:pfigure4-2-rp3
Figure 17.110

17.4.16  Two objects are released on two identical ramps. One is a sliding block (no friction), the other a rolling hoop (no slip). Both have the same mass, m, are in the same gravity field and have the same distance to travel. It takes the sliding mass 1 s to reach the bottom of the ramp. How long does it take the hoop? [Useful formula: “s=12at2”]

Filename:pfigure-s94q6p1
Figure 17.111

17.4.17  The hoop is rolled down an incline that is 30 from horizontal. It does not slip. It does not fall over sideways. It is let go from rest at t = 0.

  1. (a)

    At t = 0+ what is the acceleration of the hoop center of mass?

  2. (b)

    At t = 0+ what is the acceleration of the point on the hoop that is on the incline?

  3. (c)

    At t = 0+ what is the acceleration of the point on the hoop that is furthest from the incline?

  4. (d)

    After the hoop has descended 2 vertical meters (and traveled an appropriate distance down the incline) what is the acceleration of the point on the hoop that is (at that instant) furthest from the incline?

Filename:pfigure-spinningbrick
Figure 17.112

17.4.18  A uniform cylinder of mass m and radius r rolls down an incline without slip, as shown below. Determine: (a) the angular acceleration α of the disk; (b) the minimum value of the coefficient of friction μ that will insure no slip.

Filename:pfigure-s94f1p2
Figure 17.113

17.4.19  Race of rollers. A uniform disk with mass M0 and radius R0 is allowed to roll down the quite slip-resistant (μ=1) 30 ramp shown. It is raced against four other objects (A, B, C and D), one at a time. Who wins the races, or are there ties? First try to construct any plausible reasoning. Good answers will be based, at least in part, on careful use of equations of mechanics. Answer: Accelerations of the center of mass, where ıˆ is parallel to the slope and pointing down: (a) 𝒂cm=gsinθıˆ, (b) 𝒂cm=23gsinθıˆ, (c) 𝒂cm=12gsinθıˆ, (d) 𝒂cm=23gsinθıˆ. So, the block is fastest, all uniform disks are second, and the hollow pipe is third.

  1. (a)

    Block A has the same mass and has center of mass a distance R0 from the ramp. It rolls on massless wheels with frictionless bearings.

  2. (b)

    Uniform disk B has the same mass (MB=M0) but twice the radius (RB=2R0).

  3. (c)

    Hollow pipe C has the same mass (MC=M0) and the same radius (RC=R0).

  4. (d)

    Uniform disk D has the same radius (RD=R0) but twice the mass (MD=2M0).

Can you find a round object which will roll as fast as the block slides? How about a massless cylinder with a point mass in its center? Can you find an object which will go slower than the slowest or faster than the fastest of these objects? What would they be and why? (This problem is harder.)

Filename:s97f2
Figure 17.114

17.4.20  A roller of mass M and polar moment of inertia about the center of mass IG is connected to a spring of stiffness K by a frictionless hinge as shown in the figure. Consider two kinds of friction between the roller and the surface it moves on:

  1. 1.

    Perfect slipping (no friction), and

  2. 2.

    Perfect rolling (infinite friction).

  1. (a)

    What is the period of oscillation in the first case?

  2. (b)

    What is the period of oscillation in the second case?

Filename:pg92-2
Figure 17.115

17.4.21  A uniform cylinder of mass m and radius R rolls back and forth without slipping through small amplitudes (i.e., the springs attached at point A on the rim act linearly and the vertical change in the height of point A is negligible). The springs, which act both in compression and tension, are unextended when A is directly over C.

  1. (a)

    Determine the differential equation of motion for the cylinder’s center.

  2. (b)

    Calculate the natural frequency of the system for small oscillations.

Filename:summer95p2-3
Figure 17.116

17.4.22  Hanging disk, 2-D. A uniform thin disk of radius R and mass m hangs in a gravity field g from a pair of massless springs each with constant k. In the static equilibrium configuration the springs are vertical and attached to points A and B on the right and left edges of the disk. In the equilibrium configuration the springs carry the weight, the disk counter-clockwise rotation is ϕ=0, and the downwards vertical deflection is y=0. Assume throughout that the center of the disk only moves up and down, and that ϕ is small so that the springs may be regarded as vertical when calculating their stretch (sinϕϕ and cosϕ1).

  1. (a)

    Find ϕ¨ and y¨ in terms of some or all of ϕ,ϕ˙,y,y˙,k,m,R, and g.

  2. (b)

    Find the natural frequencies of vibration in terms of some or all of k,m,R, and g.

Filename:p-f96-p3-3
Figure 17.117

17.4.23  A disk rolls in a cylinder. For all of the problems below, the disk rolls without slip and rocks back and forth due to gravity.

  1. (a)

    Sketch. Draw a neat sketch of the disk in the cylinder. The sketch should show all variables, coordinates and dimension used in the problem.

  2. (b)

    FBD. Draw a free-body diagram of the disk.

  3. (c)

    Momentum balance. Write the equations of linear and angular momentum balance for the disk. Use the point on the cylinder which touches the disk for the angular momentum balance equation. Leave as unknown in these equations variables which you do not know.

  4. (d)

    Kinematics. The disk rolling in the cylinder is a one-degree-of-freedom system. That is, the values of only one coordinate and its derivatives are enough to determine the positions, velocities and accelerations of all points. The angle that the line from the center of the cylinder to the center of the disk makes from the vertical can be used as such a variable. Find all of the velocities and accelerations needed in the momentum balance equation in terms of this variable and its derivative. [Hint: you’ll need to think about the rolling contact in order to do this part.]

  5. (e)

    Equation of motion. Write the angular momentum balance equation as a single second order differential equation.

  6. (f)

    Simple pendulum? Does this equation reduce to the equation for a pendulum with a point mass and length equal to the radius of the cylinder, when the disk radius gets arbitrarily small? Why, or why not?

Filename:s97p3-3
Figure 17.118

A disk rolls without slip inside a bigger cylinder.

17.4.24  A uniform hoop of radius R1 and mass m rolls from rest down a semi-circular track of radius R2 as shown. Assume that no slipping occurs. At what angle θ does the hoop leave the track and what is its angular velocity ω and the linear velocity 𝒗 of its center of mass at that instant?

If the hoop slides down the track without friction, so that it does not rotate, will it leave at a smaller or larger angle θ than if it rolls without slip (as above)? Give a qualitative argument to justify your answer.

HINT: Here is a geometric relationship between angle ϕ hoop turns through and angle θ subtended by its center when no slipping occurs: ϕ=[(R1+R2)/R2]θ. (You may or may not need to use this hint.)

Filename:pfigure-s94h8p1
Figure 17.119

17.5 Rigid-object collision mechanics

This section extends the particle collisions discussion in Section 13.2 that starts on page 13.2.

2D collisions


For collisions between rigid bodies with more general motions before and after the collisions we depend on the three key ideas

  • I.

    Collision forces are big,

  • II.

    Collisions are quick, and

  • III.

    The laws of mechanics apply during the collision.

There are two extra assumptions that are needed in simple analysis:

  • IV.

    Collision forces are few. For a given rigid body there is one, or at most two non-negligible collision forces. This is the real import of idea (I) above. Because collision forces are big most other forces can be neglected.

  • V.

    The collision force(s) act at a well defined point which does not move during the collision.

Based on these assumptions one then uses linear and angular momentum balance in their time-integrated form.

Filename:sfig4-7-DH1
Figure 17.120: Two bodies collide at point C. The only non-negligible collision impulse is 𝑷 acting on body 2 (and 𝑷 on body 1) at point C. The material points on the contacting bodies are C1 and C2. The outward normal to body 1 at C1 is 𝒏ˆ.

Example: Two bodies in space

Two bodies collide at point C. The impulse acting on body 2 is 𝑷=𝑭coll𝑑t. If the mass and inertia properties of both bodies is known, as are the velocities and rotation rates before the collision we have the following linear and angular momentum balance equations for the two bodies:

𝑷=m1(𝒗G1+𝒗G1)𝑷=m2(𝒗G2+𝒗G2)𝒓C/G1×(𝑷)=Izzcm1(ω1+ω1)𝒌ˆ𝒓C/G2×𝑷=Izzcm2(ω2+ω2)𝒌ˆ. (17.73)

These make up 6 scalar equations (2 for each momentum equation, 1 for each angular momentum equation). There are 8 scalar unknowns: 𝒗G1+ (2), 𝒗G2+ (2), ω1+ (1), ω2+ (1), and 𝑷 (2). Thus the motion after the collision cannot be determined.

[Note that linear and angular momentum balance for the system would give equations which could be obtained by adding and subtracting combinations of the equations above. So adding system momentum balance equations does not add information (i.e., adds linearly dependent equations).]

So, as for 1-D collisions, momentum balance is not enough to determine the outcome of the collision. Eqns. 17.73 aren’t enough. A thousand different models and assumptions could be added to make the system solvable. But there are only two cases that are non-controversial and also relatively simple: 1) sticking collisions, and 2) frictionless collisions.

Sticking collisions

A ‘perfectly-plastic’ sticking collision is one where the relative velocities of the two contacting points are assumed to go suddenly to zero. That is

𝒗C1+=𝒗C2+

Writing 𝒗C1+=𝒗G1++(ω1+𝒌ˆ)×𝒓C/G1 and similarly for 𝒗C2 thus adds a vector equation (2 scalar equations) to the equation set 17.73. This gives 8 equations in 8 unknowns.

A little cleverness can reduce the problem to one of solving only 4 equations in 4 unknowns. Linear momentum balance for the system, angular momentum balance for the system and angular momentum balance for object 2 make up 4 scalar equations. None of these equations includes the impulse 𝑷. Because the system moves as if hinged at C1 after the collision, the state of motion after the system is fully characterized by 𝒗G1+, ω1+, and ω2+. Thus we have 4 equations in 4 unknowns.

Filename:sfig4-7-DH2
Figure 17.121: Sticking collision with an immovable object. The box sticks at A and then rotates about A. Angular momentum about point A is conserved in the collision.

Example: One body is hugely massive: collision with an immovable object

If body 2, say, is huge compared to body 1 then it can be taken to be immovable and collision problems can be solved by only considering body 1 (see fig. 17.121). In the case of a sticking collision the full state of the system after the collision is determined by ω1+. This can be found from the single scalar equation obtained from angular momentum balance about the collision point.

𝑯A = 𝑯A+
𝒓G/A×m𝒗G+Izzcmω𝒌ˆ = 𝒓G/A×m𝒗G++Izzcmω+𝒌ˆ

Because the state of the system before the collision is assumed known (the left “-” side of the equation, and because the post-collision (+) state is a rotation about A, this equation is one scalar equation in the one unknown ω+. Note that 𝑯A+ could also be evaluated as 𝑯A+=ω+IzzA𝒌ˆ. So one way of expressing the post-collision state is as

ω+=(𝒓G/A×m𝒗G+Izzcmω𝒌ˆ)𝒌ˆIzzA and 𝒗G+=ω+𝒌ˆ×𝒓G/A.

Note also that the same 𝒓G/A is used on the right and left sides of the equation because only the velocity and not the position is assumed to jump during the collision.

The collision impulse 𝑷 can then be found from linear momentum balance as

𝑷=m(𝒗G+𝒗G).

Sticking collisions are used as models of projectiles hitting targets, of robot and animal limbs making contact with the ground, of monkeys and acrobats grabbing hand holds, and of some particularly dead and frictional collisions between solids (such as when a car trips on a curb).

Frictionless collisions

The second special case is that of a frictionless collision. Here we add two assumptions:

  1. 1.

    There is no friction so 𝑷=P𝒏ˆ. The number of unknowns is thus reduced from 8 to 7.

  2. 2.

    There is a coefficient of (normal) restitution e.

The normal restitution coefficient is taken as a property of the colliding bodies. It is a given number with 0<e<1 with this defining equation:

(𝒗C2+𝒗C1+)𝒏ˆ=e(𝒗C2𝒗C1)𝒏ˆ.

This says that the normal part of the relative velocity of the contacting points reverses sign and its magnitude is attenuated by e. This adds a scalar equation to the set Eqns. 17.73 thus giving 7 scalar equations (4 momentum, 2 angular momentum, 1 restitution) for 7 unknowns (4 velocity components, 2 angular velocities and the normal impulse).

The most popular application of the frictionless collision model is for billiard or pool balls, or carom pucks. These things have relatively small coefficients of friction.

We state without proof that a frictionless collision with e=1 conserves energy.

Filename:sfig4-7-DH3
Figure 17.122: Frictionless collision between two identical round objects. Ball one is initially moving to the right, ball 2 is initially stationary. The impulse of ball 1 on ball 2 is 𝑷.

Example: Pool balls

Assume one ball approaches the other with initial velocity 𝒗G1+=vıˆ and has an elastic frictionless collision with the other ball at a collision angle of θ as shown in fig. 17.122. Defining 𝒏ˆcosθıˆsinθȷˆ we have that 𝑷=P𝒏ˆ. To determine the outcome of the equation we have the angular momentum balance equations (about the center-of-mass) which trivially tell us that

ω1+=ω2+=0

because the balls start with no spin and the frictionless collision impulses 𝑷=P𝒏ˆ and 𝑷=P𝒏ˆ have no moment about the center-of-mass. Linear momentum balance for each of the balls

P𝒏ˆ = m𝒗G1+mvıˆ
P𝒏ˆ = m𝒗G2+𝟎

gives 4 scalar equations which are supplemented by the restitution equation (using e=1)

(Δ𝒗+)𝒏ˆ=e(Δ𝒗)𝒏ˆ
vcosθ=𝒗G2+𝒏ˆ𝒗G1+𝒏ˆ

which together make 5 scalar equations in the 5 scalar unknowns 𝒗G1+, 𝒗G2+, and P (each vector has 2 unknown components). These have the solution

𝒗G1+ = vsinθ(sinθıˆ+cosθȷˆ),
𝒗G2+ = vcosθ(cosθıˆsinθȷˆ),and
P = mvcosθ.

The solution can be checked by plugging back into the momentum and restitution equations. Also, as promised, this e=1 solution conserves kinetic energy. The solution has the interesting property that the outgoing trajectories of the two balls are orthogonal for all θ but θ=0 in which case ball 1 comes to rest in the collision. [The solution can be found graphically by looking for two outgoing vectors which add to the original velocity of mass 1, where the sum of the squares of the outgoing speeds must add to the square of the incoming speed.]

Frictional collisions

For a collision with friction, but not so much that total sticking is accurate, the modeling is complex and subtle. As of this writing there are no standard acceptable ways of dealing with such situations. Commercial simulation packages should be used for such with skeptical caution. They are generally defective in that they can predict only a limited range of phenomena and/or they can create energy even with innocent input parameters.

Why is it hard to find a good collision law

Ideally one would like a rule to determine how bodies move after a collision from how they move before the collision. Such a rule would be called a collision law or a constitutive relation for collisions. That accurate collision laws are rare at best might be surmised from the basic problem that the phrase rigid body collisions is in some sense a contradiction in terms, an oxymoron. The force generated in the contact comes from material deformation, and deformation is just what we generally try to neglect when doing rigid body mechanics.

There is a temptation to say that one wants to continue to neglect deformation during the collision, but for an infinitesimal contact region. And some collision laws are formulated with this approach. Even then, there are no reliable models for the deformation in that small region, and such laws are doomed to inaccuracy in situations where the deformation is not so limited.

For complex shaped bodies touching at various points that are generally not known a priori, no collision law is reliably accurate.

SAMPLE 17.27   For two masses conservation of momentum is expressed by the vector equation m1𝒗1+m2𝒗2=m1𝒗1++m2𝒗2+. Suppose 𝒗1=𝟎,𝒗2=v0ȷˆ,𝒗1+=v1+ıˆ and 𝒗2+=v2t+𝒆ˆt+v2n+𝒆ˆn, where 𝒆ˆt=cosθıˆ+sinθȷˆ and 𝒆ˆn=sinθıˆ+cosθȷˆ.

  1. 1.

    Obtain two independent scalar equations from the momentum equation corresponding to projections in the 𝒆ˆn and 𝒆ˆt directions.

  2. 2.

    Assume that you are given another equation v+2t=v0sinθ. Set up a matrix equation to solve for v1+,v2t+, and v2n+ from the three equations.


Solution

  1. 1.

    The given equation of conservation of linear momentum is

    m1𝒗10+m2𝒗2 = m1𝒗1++m2𝒗2+
    or m2v0ȷˆ = m1v1+ıˆ+m2(v2t+𝒆ˆt+v2n+𝒆ˆn). (17.74)

    Dotting both sides of eqn. (17.74) with 𝒆ˆn gives

    m2v0(𝒆ˆnȷˆcosθ)=m1v1+(𝒆ˆnıˆsinθ)+m2v2t+(𝒆ˆn𝒆ˆt0)+m2v2n+(𝒆ˆn𝒆ˆn1)
    or m2v0cosθ=m1v1+sinθ+m2v2n+. (17.75)

    Dotting both sides of eqn. (17.74) with 𝒆ˆt gives

    m2v0(𝒆ˆtȷˆsinθ)=m1v1+(𝒆ˆtıˆcosθ)+m2v2t+(𝒆ˆt𝒆ˆt1)+m2v2n+(𝒆ˆt𝒆ˆn0)
    or m2v0sinθ=m1v1+cosθ+m2v2t+. (17.76)

    Answer: m2v0cosθ=m1v1+sinθ+m2v2n+,m2v0sinθ=m1v1+cosθ+m2v2t+

  2. 2.

    Now, we rearrange eqn. (17.75) and 17.76 along with the third given equation, v+2t=v0sinθ, so that all unknowns are on the left hand side and the known quantities are on the right hand side of the equal sign. These equations, in matrix form, are as follows.

    [m1sinθ0m2m1cosθm20010]{v1+v2t+v2n+}={m2v0cosθm2v0sinθv0sinθ}.

    This equation can be easily solved on a computer for the unknowns.

SAMPLE 17.28  Cueing a billiard ball.

Filename:sfig4-7-1
Figure 17.123:

A billiard ball is cued by striking it horizontally at a distance d=10mm above the center of the ball. The ball has mass m=0.2kg and radius r=30mm. Immediately after the strike, the center-of-mass of the ball moves with linear speed v=1m/s. Find the angular speed of the ball immediately after the strike. Ignore friction between the ball and the table during the strike.


Solution

Filename:sfig4-7-1a
Figure 17.124: FBD of the ball during the strike. The nonimpulsive forces mg,N, and f can be ignored in comparison to the strike force F.

Let the force imparted during the strike be F. Since the ball is cued by giving a blow with the cue, F is an impulsive force. Impulsive forces, such as F, are in general so large that all non-impulsive forces are negligible in comparison during the time such forces act. Therefore, we can ignore all other forces (mg,N,f) acting on the ball from its free-body diagram during the strike.

Now, from the linear momentum balance of the ball we get

Fıˆ=𝑳˙or(Fıˆ)dt=d𝑳  (Fıˆ)dt=𝑳2𝑳1

where L2L1=Δ𝑳 is the net change in the linear momentum of the ball during the strike. Since the ball is at rest before the strike, 𝑳1=m𝒗10=𝟎. Immediately after the strike, 𝒗=vıˆ=1m/s.

Thus𝑳2=m𝒗=0.2kg1m/sıˆ=0.2Nsıˆ.
Hence(Fıˆ)dt=0.2NsıˆorFdt=0.2Ns. (17.77)

To find the angular speed we apply the angular momentum balance. Let ω be the angular speed immediately after the strike and 𝝎=ω𝒌ˆ. Now,

𝑴cm=𝑯˙cm  𝑴cmdt=d𝑯cm=(𝑯cm)2(𝑯cm)1.
Because𝑯cm = Icmzz𝝎 and just before the strike, 𝝎=𝟎,
(𝑯cm)1 angular momentum just before the strike =𝟎
(𝑯cm)2 angular momentum just after the strike =Icmzzω𝒌ˆ,
𝑴cmdt = Icmzzω𝒌ˆ=25mr2ω𝒌ˆ (since for a sphere, Icmzz=25mr2).
But 𝑴cm=Fd𝒌ˆ,
therefore (Fd)dt𝒌ˆ=25mr2ω𝒌ˆ
or dconstantFdt=25mr2ω  ω=5d2mr2Fdt.

Substituting the given values and Fdt=0.2Ns from equation 17.77 we get

ω=5(0.01m)20.2kg(0.03m)20.2Ns=27.78rad/s.

The negative value makes sense because the ball will spin clockwise after the strike, but we assumed that ω was anticlockwise.

Answer: ω=27.78rad/s.

SAMPLE 17.29

Filename:sfig4-7-1b
Figure 17.125:

Falling stick. A uniform bar of length and mass m falls on the ground at an angle θ as shown in the figure. Just before impact at point C, the entire bar has the same velocity v directed vertically downwards. Assume that the collision at C is plastic, i.e., end C of the bar gets stuck to the ground upon impact.

  1. 1.

    Find the angular velocity of the bar just after impact.

  2. 2.

    Assuming θ to be small, find the velocity of end B of the bar just after impact.

Solution We are given that the impact at point C is plastic. That is, end C of the bar has zero velocity after impact. Thus end C gets stuck to the ground. Then we expect the rod to rotate about point C as the rest of the bar moves (perhaps faster) to touch the ground. The free-body diagram of the bar is shown in fig. 17.126 during the impact at point C. Note that we can ignore the force of gravity in comparison to the large impulsive force Fc due to impact at C.

Filename:sfig4-7-2
Figure 17.126: The free-body diagram of the bar during collision. The impulsive force at the point of impact C is so large that the force of gravity can be completely ignored in comparison.
  1. 1.

    Now, if we carry out angular momentum balance about point C, there will be no net moment acting on the bar, and therefore, angular momentum about the impact point C is conserved. Distinguishing the kinematic quantities before and after impact with superscripts ‘-’ and ‘+’, respectively, we get from the conservation of angular momentum about point C,

    𝑯C = 𝑯C+
    Izzcm𝝎+𝒓G/C×m𝒗G = Izzcm𝝎++𝒓G/C×m𝒗G+.

    Now, we know that 𝝎=𝟎 since every point on the bar has the same vertical velocity 𝒗=vȷˆ, and that just after impact, 𝒗G+=𝝎+×𝒓G/C where we can take 𝝎+=ω(𝒌ˆ). Thus,

    margin: Since C is a fixed point for the motion of the bar after impact, we could calculate 𝑯C+ as follows. 𝑯C+=IzzC𝝎=13m2IzzCω(𝒌ˆ).
    𝑯C = 𝒓G/C×m𝒗G=(/2)𝝀ˆ×mv(ȷˆ)
    = mv2cosθ𝒌ˆ(since 𝝀ˆ=cosθıˆ+sinθȷˆ)
    𝑯C+ = Izzcm𝝎++𝒓G/C×m(𝝎+×𝒓G/C)
    = Izzcmω𝒌ˆ+(/2)𝝀ˆ×m(ω𝒌ˆ×/2𝝀ˆ)ω/2(𝒏ˆ)
    = 112m2ω𝒌ˆ14m2ω𝒌ˆ=13m2ω𝒌ˆ.

    Now, equating 𝑯C and 𝑯C+ we get

    ω=3v2cosθ,  𝝎=3v2cosθ𝒌ˆ.

    Answer: ω=3v2cosθkˆ

  2. 2.

    The velocity of the end B is now easily found using 𝒗B=𝒗C+𝒗B/C=𝒗B/C and 𝒗B/C=𝝎×𝒓B/C. Thus,

    𝒗B/C = 𝝎×𝒓B/C=ω𝒌ˆ×𝝀ˆ
    = ω𝒏ˆ=3v2cosθ(sinθıˆ+cosθȷˆ)

    but, for small θ, cosθ1, and sinθ0. Therefore, 𝒗B/C=3v2ȷˆ.

    Thus, end B of the bar speeds up by one and a half times its original speed due to the plastic impact at C.

    Answer: vB/C=(3/2)vȷˆ

SAMPLE 17.30

Filename:sfig4-7-2a
Figure 17.127:

Tipping box. A box of mass m and dimensions 2a and 2b moves along a horizontal surface with uniform speed v. Suddenly, it bumps into an obstacle at A. Assume that the impact is plastic and point A is at the lowest level of the box. What is the maximum v the box can have so that it does not tip over after the impact.

Solution Whether the box can tip or not depends on whether it gets sufficient initial angular speed just after collision to overcome the restoring moment due to gravity about the point of rotation A. So, first we need to find the angular velocity of the box immediately following the collision. The free-body diagram of the box during collision is shown in fig. 17.128. There is an impulse 𝑷 acting at the point of impact. If we carry out the angular momentum balance about point A, we see that the impulse at A produces no moment impulse about A, and therefore, the angular momentum about point A has to be conserved. That is, 𝑯A+=𝑯A.

Filename:sfig4-7-2b
Figure 17.128: The free-body diagram of the box during collision.

Now,

𝑯A=𝒓G/A×m𝒗G=(bıˆ+aȷˆ)×mvıˆ=mav𝒌ˆ

Let the box have angular velocity 𝝎+=ω𝒌ˆ just after impact. Then,

𝑯A+ = Izzcm𝝎++𝒓G/A×m𝒗G+=Izzcmω𝒌ˆ+r𝝀ˆ×m(ω𝒌ˆ×r𝝀ˆ)
= Izzcmω𝒌ˆ+mr2ω𝒌ˆ=112(4a2+4b2)mω𝒌ˆ+m(a2+b2)ω𝒌ˆ
= 43(a2+b2)mω𝒌ˆ.

Now equating the two momenta, we get

ω+=3a4(a2+b2)v  𝝎+=3a4(a2+b2)v𝒌ˆ. (17.78)

Thus we know the angular velocity immediately after impact. Now let us find out if it is enough to get over the hill, so to speak. Once the impact is over (in a few milliseconds), the usual forces show up on the free-body diagram (see fig. 17.129).

Filename:sfig4-7-2c
Figure 17.129: The free-body diagram of the box just after the collision is over.

In this post-collision part of the motion energy is conserved. Calling just after collision ‘+’ and the top position ‘t’ (where G is directly above A), using that IzzA=43m(a2+b2) and solving for the critical condition, that the block has just enough energy to get to the top position, we have,

EP++EK+ =EKt+EPt
EPtEP+ =(EKtEK+)
mg((b2+a2)b) =(012(ω+)2IzzA)
(ω+)2 =32g(1b(b2+a2))(b2+a2).

This tells us how big ω+ has to be. Substituting this ω+ into eqn. (17.78) we find how big v has to be. We find that the box does not tip over so long as v is small enough, as given by

Answer: v2a2g(b2+a2)3/2(1b(b2+a2))/3.

SAMPLE 17.31

Filename:pfigure-s94q7p1
Figure 17.130:

Ball hits the bat. A uniform bar of mass m2=1kg and length 2=1m hangs vertically from a hinge at A. A ball of mass m1=0.25kg comes and hits the bar horizontally at point D with speed v=5m/s. The point of impact D is located at d=0.75m from the hinge point A. Assume that the collision between the ball and the bar is plastic.

  1. 1.

    Find the velocity of point D on the bar immediately after impact.

  2. 2.

    Find the impulse on the bar at D due to the impact.

  3. 3.

    Find and plot the impulsive reaction at the hinge point A as a function of d, the distance of the point of impact from the hinge point. What is the value of d which makes the impulse at A zero?

Solution The free-body diagram of the ball and the bar as a single system is shown in fig. 17.131 during impact. There is only one external impulsive force 𝑭A acting at the hinge point A. We take the ball and the bar together here so that the impulsive force acting between the ball and the bar becomes internal to the system and we are left with only one external force at A. Then, the angular momentum balance about point A yields 𝑯˙A=𝟎 since there is no net moment about A. Thus the angular momentum about A is conserved during the impact.

  1. 1.

    Let us distinguish the kinematic quantities just before impact and immediately after impact with superscripts ‘-’ and ‘+’, respectively. Then, from the conservation of angular momentum about point A, we get 𝑯A=𝑯A+. Now,

    Filename:pfigure-s94h7p2
    Figure 17.131: The free-body diagram of the ball and the bar together during collision. The impulsive force at the point of impact is internal to the system and hence, does not show on the free-body diagram.
    𝑯A = (𝑯A)ball+(𝑯A)bar
    = 𝒓D/A×m1𝒗+IzzA𝝎
    = dȷˆ×m1v(ıˆ)+𝟎=m1dv𝒌ˆ.
      Similarly,
    𝑯A+ = 𝒓D/A×m1𝒗++IzzA𝝎+
      but, 𝒗+=𝝎+×𝒓D/A=ωdıˆ, where 𝝎+=ω𝒌ˆ (let). Hence,
    𝑯A+ = dȷˆ×m1(ωdıˆ)+13m2(2)2ω𝒌ˆ
    = (m1d2+43m22)ω𝒌ˆ.

    Equating the two momenta, we get

    ω = m1dvm1d2+(4/3)m22
    = vd(1+43m2m1(d)2)
     𝒗D = 𝝎+×𝒓D/A=ωd(ıˆ)
    = v1+43m2m1(d)2ıˆ.

    Now, substituting the given numerical values, v=5m/s, m1=0.25kg, m2=1kg, =0.5m, and d=0.75m, we get 𝒗D=2.08m/sıˆ

    Answer: vD=2.08m/sıˆ

  2. 2.

    To find the impulse at D due to the impact, we can consider either the ball or the bar separately, and find the impulse by evaluating the change in the linear momentum of the body. Let us consider the ball since it has only one impulse acting on it. The free-body diagram of the ball during impact is shown in fig. 17.132. From the linear impulse-momentum relationship we get,

    𝑷D=𝑭Ddt = 𝑳+𝑳=m1(𝒗+𝒗)
    = m1(v1+43m2m1(d)2ıˆ+vıˆ)
    = m1v(111+43m2m1(d)2)ıˆ.
    Filename:ch4-7-c
    Figure 17.132: Separate free-body diagrams of the ball and the bar during collision.

    Substituting the given numerical values, we get 𝑷D=0.73kgm/sıˆ. The impulse on the bar is equal and opposite. Therefore, the impulse on the bar is 𝑷D=0.73kgm/sıˆ.

    Answer: Impulse at D=0.73kgm/sıˆ

  3. 3.

    Now that we know the impulse at D, we can easily find the impulse at A by applying impulse-momentum relationship to the bar. Since, the bar is stationary just before impact, its initial momentum is zero. Thus, for the bar,

    (𝑭A𝑭D)dt=𝑳+𝑳=𝑳+=m2𝒗cm+.

    Denoting the impulse at A with 𝑷A, the mass ratio m2/m1 by mr, and the length ratio /d by q, and noting that 𝒗cm+=ω𝒌ˆ×ȷˆ=ωıˆ, we get

    𝑷A 𝑭Adt=𝑭Ddt+m2(ωıˆ)
    = m1v(111+43mrq2)ıˆm2vd(1+43mrq2)ıˆ
    = m1v(43mrq21+43mrq2)ıˆm2v(q1+43mrq2)ıˆ
    = (4/3)m2q2m2q1+43mrq2vıˆ=q(4q3)3(1+43mrq2)m2vıˆ.

    Now, we are ready to graph the impulse at A as a function of q/d. However, note that a better quantity to graph will be PA/(m1v), that is, the nondimensional impulse at A, normalized with respect to the initial linear momentum m1v of the ball. The plot is shown in fig. 17.133. It is clear from the plot, as well as from the expression for 𝑷A, that the impulse at A is zero when q=3/4 or d=4/3=2/3(2), that is, when the ball strikes at two thirds the length of the bar. Note that this location of the impact point is independent of the mass ratio mr.

    Filename:pfigure-blue-87d-1
    Figure 17.133: Plot of normalized impulse at A as a function of q=/d.

    Answer: d=2/3(2) for PA=𝟎

Comment: This particular point of impact D (when d=2/3(2)) which induces no impulse at the support point A is called the center of percussion. If you imagine the bar to be a bat or a racquet and point A to be the location of your grip, then hitting a ball at D gives you an impulse-free shot. In sports, point D is called a sweet spot.

SAMPLE 17.32

Filename:Mikef91p3
Figure 17.134:

Flying dish collides with solar panel. A uniform rectangular plate of dimensions 2a=2m and 2b=1m and mass mP=2kg drifts in space at a uniform speed vp=10m/s (relative to some local Newtonian reference frame) in the direction shown. A circular disk of radius R=0.25m and mass mD=1kg is heading towards the plate at a linear speed v0=vD=1m/s directed normal to the facing edge of the plate. In addition, the disk is spinning at ωD=5rad/s in the clockwise direction. The plate and the disk collide at point A of the plate, located at d=0.8m from the center of the long edge. Assume that the collision is frictionless and purely elastic. Find the linear and angular velocities of the plate and the disk immediately after the collision.

Solution We are given the linear and angular velocities of the two bodies before their collision. We have to find their linear and angular velocities after the collision. Let the linear velocities of the plate and the disk after the collision be 𝒗P+ and 𝒗D+, respectively. We assume that the collision takes place in the plane of the two bodies (that is, all motions before and after collision are planar). Then the unknown angular velocities of the plate and the disk after the collision can be assumed as 𝝎P+=ωP+𝒌ˆ and 𝝎D+=ωD+𝒌ˆ. Thus, in total, we have 6 scalar unknowns here — 4 for linear velocities of the disk and the plate (each velocity has two components) and 2 for the two angular velocities.

Since we are trying to find velocities after a collision, given the velocities before the collision, we can use linear and angular momentum-impulse relations. Let 𝑷 be the impulse acting during the collision. Note that the collision is frictionless and hence 𝑷 must act normal to the two colliding surfaces. Therefore, the direction of the impulse is known, only its magnitude is not known. Thus, in total, we have 7 unknowns now. To solve for them, we need 7 independent equations.

We have 6 independent equations from the linear and angular impulse-momentum balance for the two bodies (3 each). We need one more equation. That equation is the relationship between the normal components of the relative velocities of approach and departure with the coefficient of restitution e (=1 for elastic collision). Thus we have 7 equations for 7 unknowns. First we set them up, then we solve them with a calculator or computer.

Filename:pfigure-s94h7p3
Figure 17.135: Separate free-body diagrams of the disk and the plate during collision. The impulse 𝑷 acts normal to the colliding surfaces and hence we can assume that on the plate, 𝑷=Pıˆ. Therefore, the impulse on the disk is Pıˆ.
Filename:p-f96-p2-3
Figure 17.136: The free-body diagram of the disk and the plate together during collision. The impulsive force at the point of impact is internal to the system and hence, does not show on the free-body diagram .

The free-body diagrams of the disk and the plate together and the two separately are shown in fig. 17.136 and 17.135, respectively. Using an xy coordinate system oriented as shown in fig. 17.136, we can write

LMB for disk:mD(𝒗D+𝒗D)=PıˆLMB for plate:mP(𝒗P+𝒗P)=PıˆAMB for disk:IDcm(𝝎D+𝝎D)=𝟎AMB for plate:IPcm(𝝎P+𝝎P)=𝒓A/P×Pıˆkinematics:𝒗AD+𝒗AP+=e(𝒗AP𝒗AD)

where, in the last equation 𝒗AD and 𝒗AP refer to the velocities of the material points located at A on the disk and on the plate, respectively. Other linear velocities in the equations above refer to the velocities at the center-of-mass of the corresponding bodies. We are given that 𝒗D=vDıˆ,𝒗P=vPıˆ,𝝎D=ωD𝒌ˆ, and 𝝎P=𝟎. Let us assume that 𝝎D+=ωD+𝒌ˆ,𝝎P+=ωP+𝒌ˆ,𝒗D+=vDx+ıˆ+vDy+ȷˆ, and similarly, 𝒗P+=vPx+ıˆ+vPy+ȷˆ. Then,

𝒗AD = 𝒗D+𝝎D×𝒓A/O=vDıˆωDRȷˆ
𝒗AD+ = 𝒗D++𝝎D+×𝒓A/O=vDx+ıˆ+(vDy++ωD+R)ȷˆ
𝒗AP = 𝒗P=vPıˆ
𝒗AP+ = 𝒗P++𝝎P+×𝒓A/P=(vPx+ωP+d)ıˆ+(vPy+ωP+b)ȷˆ.

Substituting these quantities in the kinematics equation above and dotting with ıˆ, the unit normal at A, we get

vDx+vPx++ωP+d=e1(vPvD)=vPvD. (17.79)

Now, let us extract the scalar equations from the impulse-momentum equations for the disk and the plate by dotting with appropriate unit vectors. Dotting LMB for the disk with ıˆ and ȷˆ, respectively, we get

mD(vDx+vD) = P (17.80)
mDvDy+ = 0. (17.81)
Dotting LMB for the plate with ıˆ and ȷˆ, respectively, we get
mP(vPx++vP) = P (17.82)
mPvPy+ = 0. (17.83)
Dotting AMB for the disk and the plate with 𝒌ˆ, we get
IDcm(ωD++ωD) = 0 (17.84)
IPcmωP+ = Pd. (17.85)

We have all the equations we need. Let us rearrange these equations in a matrix form, taking the known quantities to the right and putting all unknowns to the left side. We then, write eqns. (17.80)–(17.85), and then eqn. (17.79) as

[mD0000010mD0000000mP0001000mP0000000IDcm0000000IPcmd10100d0]{vDx+vDy+vPx+vPy+ωD+ωP+P}={mDvD0mPvP0IDcmωD0vPvD}.

Substituting the given numerical values for the masses and the pre-collision velocities, and the moments of inertia, IDcm=(1/2)mDR2 and IPcm=(1/12)mP(4a2+4b2), and then solving the matrix equation on a computer, we get,

𝒗D+ =8.70m/sıˆ,𝒗P+=5.15m/sıˆ
𝝎D+ =5rad/s𝒌ˆ,𝝎P+=9.31rad/s𝒌ˆ
P =9.70kgm/s.

You can easily check that the results obtained satisfy the conservation of linear momentum for the plate and the disk taken together as one system.


Comments: In this particular problem, the equations are simple enough to be solved by hand. For example, eqns. (17.81), (17.83), and (17.84) are trivial to solve and immediately give, vDy+=0,vPy+=0, and ωD+=ωD=5rad/s. The rest of the equations can be solved by usual eliminations and substitutions, etc. The key is that you have 7 linearly independent equations for the 7 unknowns.

Problems for 17.5 Collisions

17.5.1  The two blocks shown in the figure are identical except that one rests on two springs while the other one sits on two massless wheels. Draw free-body diagrams of each mass as each is struck by a hammer. Here we are interested in the free-body diagrams only during collision. Therefore, ignore all forces that are much smaller than the impulsive forces. State in words why the forces you choose to show should not be ignored during the collision.

Filename:s97p3-2
Figure 17.137

17.5.2  These problems concern two colliding masses. In (a) the smaller mass hits the hanging mass from above at an angle 45 with the vertical. In (b) the second case the smaller mass hits the hanging mass from below at the same angle. Assuming perfectly elastic impact between the balls, find the velocity of the hanging mass just after the collision. [Note, these problems are not well posed and can only be solved if you make additional modeling assumptions.]

Filename:pfigure-dynbalancers
Figure 17.138

17.5.3  A narrow pole is in the middle of a pond with a 10m rope tied to it. A frictionless ice skater of mass 50kg and speed 3m/s grabs the rope. The rope slowly wraps around the pole. What is the speed of the skater when the rope is 5m long? (A tricky question.)

Filename:pfigure-s94h7p5
Figure 17.139

17.5.4  The masses m and 3m are joined by a light-weight bar of length 4. If point A in the center of the bar strikes fixed point B vertically with velocity V0, and is not permitted to rebound, find θ˙ of the system immediately after impact.

Filename:pfigure-blue-81-1
Figure 17.140

Neglect gravity!

17.5.5  Two equal masses each of mass m are joined by a massless rigid rod of length . The assembly strikes the edge of a table as shown in the figure, when the center of mass is moving downward with a linear velocity v and the system is rotating with angular velocity ϕ˙ in the counter-clockwise sense. The impact is ’elastic’. Find the immediate subsequent motion of the system, assuming that no energy is lost during the impact and that there is no gravity. Show that there is an interchange of translational and rotational kinetic energy.

Filename:pfigure-blue-137-2
Figure 17.141

17.5.6  In the absence of gravity, a thin rod of mass m and length is initially tumbling with constant angular speed ωo, in the counter-clockwise direction, while its mass center has constant speed vo, directed as shown below. The end A then makes a perfectly plastic collision with a rigid peg O (via a hook). The velocity of the mass center happens to be perpendicular to the rod just before impact.

  1. (a)

    What is the angular speed ωf immediately after impact?

  2. (b)

    What is the angular speed 10 seconds after impact? Why?

  3. (c)

    What is the loss in energy in the plastic collision?

Filename:pfigure-blue-78-2
Figure 17.142

17.5.7  A gymnast of mass m and extended height L is performing on the uneven parallel bars. About the x, y, z axes which pass through her center of mass, her radii of gyration are kx, ky, and kz, respectively. Just before she grasps the top bar, her fully extended body is horizontal and rotating with angular rate ω; her center of mass is then stationary. Neglect any friction between the bar and her hands and assume that she remains rigid throughout the entire stunt.

  1. (a)

    What is the gymnast’s rotation rate just after she grasps the bar? State clearly any approximations/assumptions that you make.

  2. (b)

    Calculate the linear speed with which her hips (CM) strike the lower bar. State all assumptions/approximations.

  3. (c)

    Describe (in words, no equations please) her motion immediately after her hips strike the lower bar if she releases her hands just prior to this impact.

Filename:pfigure-blue-130-2
Figure 17.143

17.5.8  An acrobat modeled as a rigid body with uniform rigid mass m of length l. She falls without rotation in the position shown from height h where she was stationary. She then grabs a bar with a firm but slippery grip. What is h so that after the subsequent motion the acrobat ends up in a stationary handstand? [ Hint: What quantities are preserved in what parts of the motion?] Answer: h=2L/3

Filename:pfigure-s94h10p4
Figure 17.144

17.5.9  A crude see-saw is built with two supports separated by distance d about which a rigid plank (mass m, length L) can pivot smoothly. The plank is placed symmetrically, so that its center of mass is midway between the supports when the plank is at rest.

  1. (a)

    While the left end is resting on the left support, the right end of the plank is lifted to an angle θ and released. At what angular velocity ω1 will the plank strike the right hand support?

  2. (b)

    Following the impact, the left end of the plank can pivot purely about the right end if d/L is properly chosen and the right end does not bounce. Find ω2 under these circumstances.

Filename:pfigure-blue-81-2
Figure 17.145

17.5.10  Baseball bat. In order to convey the ideas without making the calculation too complicated, some of the simplifying assumptions here are highly approximate. Assume that a bat is a uniform rigid stick with length L and mass ms. The motion of the bat is a pivoting about one end held firmly in place with hands that rotate but do not move. The swinging of the bat occurs by the application of a constant torque Ms at the hands over an angle of θ=π/2 until the point of impact with the ball. The ball has mass mb and arrives perpendicular to the bat at an absolute speed vb at a point a distance from the hands. The collision between the bat and the ball is completely elastic.

  1. (a)

    To maximize the speed vhit of the hit ball How heavy should a baseball bat be? Where should the ball hit the bat? Here are some hints for one way to approach the problem.

    • Find the angular velocity of the bat just before collision by drawing a FBD of the bat etc.

    • Find the total energy of the ball and bat system just before the collision.

    • Draw a FBD of the ball and of the bat during the collision (with this model there is an impulse at the hands on the bat). Call the magnitude of the impulse of the ball on the bat (and vice versa) F𝑑t.

    • Use various momentum equations to find the angular velocity of the bat and velocity of the ball just after the collision in terms of F𝑑t and other quantities above. Use these to find the energy of the system just after collision.

    • Solve for the value of F𝑑t that conserves energy. As a check you should see if this also predicts that the relative separation speed of the ball and bat (at the impact point) is the same as the relative approach speed (it should be).

    • You now can calculate vhit in terms of mb,ms,Ms,L,, and vb.

    • Find the maximum of the above expression by varying ms and . Pick numbers for the fixed quantities if you like.

  2. (b)

    Can you explain in words what is wrong with a bat that is too light or too heavy?

  3. (c)

    Which aspects of the model above do you think lead to the biggest errors in predicting what a real ball player should pick for a bat and place on the bat to hit the ball?

  4. (d)

    Describe as clearly as possible a different model of a baseball swing that you think would give a more accurate prediction. (You need not do the calculation).