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 x⁢y 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 x⁢y 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 x⁢y 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 x⁢y plane. This problem does fit in with the methods of this chapter. The symmetry of the sphere with respect to the x⁢y 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 chapter††margin: 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 x⁢y 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]x⁢y=[𝒓⇀0′/0]x⁢y⏟displacement+[R⁢(θ)]⁢[𝒓⇀P/0′]x′⁢y′⏟rotation (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/0∗yP/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 θ=π/3⁢rad=60∘ relative to that configuration, then its new position is:

[𝒓⇀P/0]x⁢y = [𝒓⇀0′/0]x⁢y+[R⁢(θ)]⁢[𝒓⇀P/0′]x′⁢y′
= [x0′/0y0′/0]+[cos⁡θ−sin⁡θsin⁡θcos⁡θ]⁢[xP/0∗yP/0∗]
= {[34]+[cos⁡π/3−sin⁡π/3sin⁡π/3cos⁡π/3]⁢[12]}⁢m
= [3.5−35+3/2]⁢m
⇒𝒓⇀P/0 = ((3.5−3)⁢ıˆ+(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∗]x⁢y = [𝒓⇀P/0′∗]x′⁢y′=[x′y′] (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 x⁢y 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 x⁢y 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 x⁢y 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 x⁢y plane .

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

𝒗⇀B/A = dd⁢t⁢(xB/A′⁢ıˆ′+yB/A′⁢ȷˆ′)
= xB/A′⁢dd⁢t⁢ıˆ′+yB/A′⁢dd⁢t⁢ȷˆ′
= 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]x⁢y = dd⁢t⁢[xB/AyB/A]
= dd⁢t⁢{[cos⁡θ−sin⁡θsin⁡θcos⁡θ]⁢[xB/A∗yB/A∗]}
= [−θ˙⁢sin⁡θ−θ˙⁢cos⁡θθ˙⁢cos⁡θ−θ˙⁢sin⁡θ]⁢[xB/A∗yB/A∗]
= [0−ωω0]⁢[cos⁡θ−sin⁡θsin⁡θcos⁡θ]⁢[xB/A∗yB/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 = dd⁢t⁢𝒗⇀B/A (17.17)
= dd⁢t⁢{𝝎⇀ℬ×𝒓⇀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]x⁢y=[R⁢(θ)]⁢[𝒓⇀B/A∗]x′⁢y′

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=ω0⁢A+θ˙.]

Computer graphics

Given one point given by the x⁢y 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 45⁢deg rotation to

[xy]=[R]⁢[x0y0] = [cos⁡π4−sin⁡π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]≡[x⁢y 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 know††margin: 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 5⁢rad/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.3⁢m/s⁢ıˆ+0.6⁢m/s⁢ȷˆ and the plate has angular velocity 𝝎⇀=2⁢rad/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.2⁢m⁢ıˆ)=−0.2⁢m⁢ıˆ

    Thus,

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

    Answer: v⇀A=(0.3⁢ıˆ+0.2⁢ȷˆ)⁢m/s

  2. 2.

    Let x⁢y axes with origin at A be fixed in space with x-axis along AB at the instant of interest. Let x′⁢y′ 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 x⁢y are fixed in space and the axes x′⁢y′ are fixed to the rotating plate in the orientation shown. Thus in the x′⁢y′ axes, the coordinates of point C will always be xC′=0.2⁢m,yC′=0.
    [xC′yC′]=[0.2⁢m0]

    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]x⁢y=[𝒗⇀A]x⁢y+[0−ωω0]⁢[cos⁡θ−sin⁡θsin⁡θcos⁡θ]⁢[xC/A′yC/A′].

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

    [𝒗⇀C]x⁢y = [0.3⁢m/s0.2⁢m/s]+[0−2⁢rad/s2⁢rad/s0]⁢[12−323212]⁢[0.2⁢m0]
    = [0.3⁢m/s0.2⁢m/s]+[−0.35⁢m/s0.2⁢m/s]=[−0.05⁢m/s0.4⁢m/s].

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

    Answer: 𝒗⇀C=−0.05⁢m/s⁢ıˆ+0.4⁢m/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 𝝎⇀=5⁢rad/s⁢𝒌ˆ and the linear velocity of a point C on the body is 𝒗⇀C=2⁢m/s⁢ıˆ−5⁢m/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|=(2⁢m/s⁢ıˆ−5⁢m/s⁢ȷˆ)/(22+52⁢m/s)=0.37⁢ıˆ−0.93⁢ȷˆ, we can immediately write 𝒏ˆ=0.93⁢ıˆ+0.37⁢ȷˆ. So, now,

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

Thus

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

Answer: r⇀O/C=1⁢m⁢ıˆ+0.4⁢m⁢ȷˆ

Alternatively,

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

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

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

2−5⁢y = 0
−5+5⁢x = 0.

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

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

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 x⁢y-plane. At an instant when the baton axis is at θ=60∘ from the horizontal, the velocity of end A of the baton is 𝒗⇀A=2⁢m/s⁢ıˆ+3⁢m/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 ℓ=12⁢m 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 x⁢y-plane, let 𝝎⇀=ω⁢𝒌ˆ. Then,

Filename:sfig8-6-2
Figure 17.17:
𝒗⇀B = 𝒗⇀A+𝝎⇀×𝒓⇀B/A=𝒗⇀A+ω⁢𝒌ˆ×ℓ⁢(cos⁡θ⁢ıˆ−sin⁡θ⁢ȷˆ)⏟𝒓⇀B/A
or−vB⁢ıˆ = (2⁢ıˆ+3⁢ȷˆ)⁢m/s+ω⁢ℓ⁢(cos⁡θ⁢ȷˆ+sin⁡θ⁢ıˆ)
= (2⁢ıˆ+3⁢ȷˆ)⁢m/s+ω⁢⋅⁢12⁢m⁢⋅⁢(12⁢ȷˆ+32⁢ıˆ)

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

0 = 3⁢m/s+ω2⁢m⁢⋅⁢12
⇒ ⁢ω = −3⁢ms⁢⋅⁢41⁢m=−4⁢3⁢rad/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/s−3⁢m/s⁢⋅⁢(12⁢ȷˆ+32⁢ıˆ)
= (2−32)⁢m/s⁢ıˆ+(3−32)⁢m/s⁢ȷˆ
= 12⁢m/s⁢ıˆ+32⁢m/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 10⁢ft 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 10⁢ft/s and 12⁢ft/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 =60⁢mph⁢ıˆ=88⁢ft/s⁢ıˆ
𝒂⇀A = acceleration of the truck =3⁢mph/s=4.4⁢ft/s2⁢ıˆ
|𝒗⇀B/A| = vB/A= magnitude of relative velocity of B =10⁢ft/s
|𝒂⇀B/A| = aB/A= magnitude of relative acceleration of B =12⁢ft/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ℓ=10⁢ft/s10⁢ft=1⁢rad/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/A2ℓ2−ω4=(12⁢ft/s210⁢ft)2−(1⁢rad/s)4
    = ±0.663⁢rad/s2.

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

    Answer: ω⇀=1⁢rad/s⁢kˆ,ω⇀˙=0.663⁢rad/s2⁢kˆ

  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⁡θ⁢ıˆ)
    = 88⁢ft/s⁢ıˆ+10⁢ft/s⁢(12⁢ȷˆ−12⁢ıˆ)
    = 80.93⁢ft/s⁢ıˆ+7.07⁢ft/s⁢ȷˆ.
    𝒂⇀B = 𝒂⇀A+𝒂⇀B/A
    = 𝒂⇀A+𝝎⇀˙×𝒓⇀B/A−ω2⁢𝒓⇀B/A
    = aA⁢ıˆ+ω˙⁢𝒌ˆ×ℓ⁢(cos⁡θ⁢ıˆ+sin⁡θ⁢ȷˆ)−ω2⁢ℓ⁢(cos⁡θ⁢ıˆ+sin⁡θ⁢ȷˆ)
    = (aA−ω˙⁢ℓ⁢sin⁡θ−ω2⁢ℓ⁢cos⁡θ)⁢ıˆ+(ω˙⁢ℓ⁢cos⁡θ−ω2⁢ℓ⁢sin⁡θ)⁢ȷˆ
    = (4.4⁢ft/s2−0.66s2⋅10⁢ft⋅12−1s2⋅10⁢ft⋅12)⁢ıˆ
    +(0.66s2⋅10⁢ft⋅12−1s2⋅10⁢ft⋅12)⁢ȷˆ
    = −7.34⁢ft/s2⁢ıˆ−2.40⁢ft/s2⁢ȷˆ.

    Answer: v⇀B=(80.93⁢ıˆ+7.07⁢ȷˆ)⁢ft/s,a⇀B=−(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⁡θ⁢ȷˆ)−ω2⁢ℓ2⁢(cos⁡θ⁢ıˆ+sin⁡θ⁢ȷˆ)
    = aA⁢ıˆ+ω˙⁢ℓ2⁢(cos⁡θ⁢ȷˆ−sin⁡θ⁢ıˆ)−ω2⁢ℓ2⁢(cos⁡θ⁢ıˆ+sin⁡θ⁢ȷˆ)
    = 4.4⁢ft/s2⁢ıˆ+0.663⁢rad/s2⋅10⁢ft2⋅(12⁢ȷˆ−12⁢ıˆ)
    −(1⁢rad/s)2⋅10⁢ft2⋅(12⁢ıˆ+12⁢ȷˆ)
    = −1.48⁢ft/s2⁢ıˆ−1.19⁢ft/s2⁢ȷˆ.

    Answer: a⇀G=−(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 ℓ=1⁢m 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)=θ0⁢sin⁡(λ⁢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
    θ⁢𝒌ˆ = θ0⁢sin⁡(λ⁢t)⁢𝒌ˆ
    ⇒ ⁢𝝎⇀ ≡ θ˙⁢𝒌ˆ=θ0⁢λ⁢cos⁡(λ⁢t)⁢𝒌ˆ (17.28)
    ⇒ ⁢𝝎⇀˙ ≡ θ¨⁢𝒌ˆ=−θ0⁢λ2⁢sin⁡(λ⁢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⁡θ⁢ıˆ−ℓ⁢θ˙2⁢sin⁡θ⁢ıˆ+ℓ⁢θ˙2⁢cos⁡θ⁢ȷˆ
    = ℓ⁢(θ¨⁢cos⁡θ−θ˙2⁢sin⁡θ)⁢ıˆ+[t/3+ℓ⁢(θ¨⁢sin⁡θ+θ˙2⁢cos⁡θ)]⁢ȷˆ

    where θ=θ0⁢sin⁡(λ⁢t),θ˙=θ0⁢λ⁢cos⁡(λ⁢t), and θ¨=−θ0⁢λ2⁢sin⁡(λ⁢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/18−ℓ⁢cos⁡θ)⁢ȷˆ.

    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 x′⁢y′ 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 x′⁢y′ coordinates in terms of given quantities.

  2. (b)

    Write the position vector of point P relative to O in the x⁢y 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=30⁢cm, r=25⁢cm, θ=60∘, and ϕ=45∘, find [𝒓⇀C/0]x⁢y, and [𝒓⇀P/0]x⁢y at the instant shown.

  6. (f)

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

Filename:Danef94s1q2
Figure 17.22:

17.1.2  A uniform rigid rod AB of length ℓ=1⁢m 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=−1⁢m/s⁢ıˆ and 𝒗⇀B=1⁢m/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=−2⁢ft/s⁢(ıˆ+ȷˆ) and 𝒗⇀D=−(2⁢ft/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, ρ=25⁢mi, and its speed increases uniformly from 50⁢mph to 60⁢mph. 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.2⁢m/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: ω⇀=−v0⁢sin⁡ϕℓ⁢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=60⁢rpm. The radius of the disk is R=0.25⁢m.

  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 ℓ=2⁢m 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=2⁢s. At the instant shown, θ=60∘, and rod AC has an angular velocity 𝝎⇀=0.5⁢rad/s⁢𝒌ˆ and angular acceleration 𝜶⇀=1⁢rad/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 ℓ=4⁢ft 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)=θ0⁢sin⁡ω⁢t where θ0=π/6 and ω=π⁢rad/s. The cart moves at a constant speed v0=2⁢ft/s on a sinusoidal track represented by y⁢(x)=y0⁢sin⁡λ⁢x where y0=0.5⁢ft 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=10⁢m/s⁢ıˆ when the center of mass is at its highest point h=6⁢m. 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 measures††margin: 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+Iz⁢zcm⁢𝝎⇀⏞Angular momentum due to motion relative to the center-of-mass (17.33)
and ⁢𝑯⇀˙/O = 𝒓⇀cm/O×(mtot⁢𝒂⇀cm)+Iz⁢zcm⁢𝝎⇀˙. (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=Iz⁢zcm⁢𝝎⇀ and 𝑯⇀˙cm=Iz⁢zcm⁢𝝎⇀˙. (17.35)

The kinetic energy and its rate of change are given by

EK = 12⁢mtot⁢vc⁢m2⏟𝒗⇀cm⋅𝒗⇀cm⏞kinetic energy from center-of-mass motion +12⁢Iz⁢zcm⁢ω2⏞kinetic energy relative to the center-of-mass (17.36)
and ⁢E˙K = mtot⁢vc⁢m⁢v˙c⁢m⏟𝒗⇀cm⋅𝒂⇀cm+Iz⁢zcm⁢ω⁢ω˙ (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)+Iz⁢zcm⁢𝝎⇀˙(b)or∑𝑴⇀cm=Iz⁢zcm⁢𝝎⇀˙,andPower:𝑭⇀tot⋅𝒗⇀cm+𝝎⇀⋅𝑴⇀cm=mtot⁢v⁢v˙+Iz⁢zcm⁢ω⁢ω˙.(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=mtot⁢acmx,

    • –

      {LMB}⋅ȷˆ or, equivalently, ∑Fy=mtot⁢acmy, and

    • –

      {AMB}⋅𝒌ˆ or, equivalently, ∑Mcm=Iz⁢zcm⁢ω˙.

  • •

    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⇒𝟎⇀ = Iz⁢zcm⁢ω˙⁢𝒌ˆ

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

∑𝑴⇀cm⁢i=𝑯⇀˙cm ⇒ {𝒓⇀A/G×F⁢ıˆ=Iz⁢zcm⁢θ¨⁢𝒌ˆ}
{}⋅𝒌ˆ ⇒ θ¨+F⁢ℓIz⁢zcm⁢sin⁡θ=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×𝑭⇀ = Iz⁢zcmω˙𝒌ˆ}
{}⋅𝒌ˆ ⇒ Iz⁢zcm⁢ω˙ = 𝒓⇀C/G×𝑭⇀⋅𝒌ˆ.

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

ω˙=c⁢ℓIz⁢zcm⁢(𝝀ˆ×𝒗⇀⋅𝒌ˆ−θ˙⁢ℓ)

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/Iz⁢zcm,

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=∑mi⁢aix,(a){∑𝑭⇀i=∑mi⁢𝒂⇀i}⋅ȷˆ⇒∑Fiy=∑mi⁢aiy,and(b){∑𝒓⇀i×𝑭⇀i=∑𝒓⇀i×mi⁢𝒂⇀i}⋅𝒌ˆ⇒∑rix⁢Fiy−riy⁢Fix=∑mi⁢(rix⁢aiy−riy⁢aix).(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 x⁢y 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=Iz⁢zcm⁢ω.

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

H˙cmz=Iz⁢zcm⁢ω˙.

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))⋅𝒌ˆ+Iz⁢zcm⁢ω˙

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 x⁢y plane, or its mass is symmetrically distributed about the x⁢y plane, and it is in planar motion in the x⁢y plane then

∑mi⁢aiz=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=∫dd⁢t⁢𝒓⇀⁢𝑑m=dd⁢t⁢∫𝒓⇀⁢𝑑m=dd⁢t⁢(mtot⁢𝒓⇀cm)
= mtot⁢dd⁢t⁢𝒓⇀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×𝒗⇀cm⁢mtot.

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×𝒗⇀cm⁢mtot.

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×𝒗⇀cm⁢mtot.

Then defining Iz⁢zcm≡∫r/cm2⁢𝑑m we get the desired result:

𝑯⇀/O = 𝒓⇀cm×𝒗⇀cm⁢mtot+Iz⁢zcm⁢ω⁢𝒌ˆ.

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

𝑯⇀˙/O = 𝒓⇀cm×𝒂⇀cm⁢mtot+Iz⁢zcm⁢ω˙⁢𝒌ˆ.

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⏟𝟎⇀×𝒗⇀cm⁢mtot+Iz⁢zcm⁢ω⁢𝒌ˆ=Iz⁢zcm⁢ω⁢𝒌ˆ.

and similarly

𝑯⇀˙cm = Iz⁢zcm⁢ω˙⁢𝒌ˆ.

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
= 12⁢mtot⁢vc⁢m2+12⁢Iz⁢zcm⁢ω2

and

E˙K ≡ dd⁢t⁢EK
= mtot⁢v⁢v˙+Iz⁢zcm⁢ω⁢ω˙.

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𝑭⇀⋅𝒗⇀⁢d⁢t⏟d⁢𝒓⇀
= ∫𝒓⇀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:

W12≡∫t1t2P⁢𝑑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

𝑭⇀hand⁢on⁢train.

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

Pforce⁢on⁢train=𝑭⇀hand⁢on⁢train⋅𝒗⇀train≠0.

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:

𝑭⇀train⁢on⁢hand=−𝑭⇀hand⁢on⁢train

while the velocity of the hand is zero so

Pforce⁢on⁢hand=𝑭⇀train⁢on⁢hand⋅𝒗⇀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 B⁢C. 𝑨⇀×(𝑩⇀×𝑪⇀) is perpendicular to both 𝑨⇀ and 𝑩⇀×𝑪⇀, so it is perpendicular to the normal to the plane of B⁢C. 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 ℓ=1⁢m and mass mr=1⁢kg, and a rigid square plate of side 1⁢m and mass mp=10⁢kg are launched in motion on a frictionless plane (e.g., an ice hockey rink) with exactly the same initial velocities 𝒗⇀cm⁢(0)=10⁢m/s⁢ıˆ+1⁢m/s⁢ȷˆ and 𝝎⇀⁢(0)=1⁢rad/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=1⁢sec and at t=3⁢sec.


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 x⁢y-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 x⁢y-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+𝒗⇀cm0⁢t (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|=|𝒗⇀cm0⁢t|) at time t. Similarly, the angular momentum balance about the center-of-mass of the rod gives

    ∑𝑴⇀cm = 𝑯⇀˙cm
    𝟎⇀ = Iz⁢zcm⁢𝝎⇀˙
    ⇒ ⁢𝝎⇀ = ∫𝟎⇀⁢𝑑t=constant=𝝎⇀0=θ˙0⁢𝒌ˆ
    ⇒ ⁢θ = ∫θ˙0⁢𝑑t=θ0+θ˙0⁢t (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=θ˙0⁢t) during any given time.

  2. 2.

    We can find the position of the rod at t=1⁢s and t=3⁢s 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=1⁢s) = 𝒗⇀cm0⋅(1⁢s)=(10⁢m/s⁢ıˆ+1⁢m/s⁢ȷˆ)⋅(1⁢s)⁢10⁢m⁢ıˆ+1⁢m⁢ȷˆ
    𝒓⇀cm⁢(t=3⁢s) = 𝒗⇀cm0⋅(3⁢s)=30⁢m⁢ıˆ+3⁢m⁢ȷˆ.
    θ⁢(t=1⁢s) = θ˙0⋅(1⁢s)=(1⁢rad/s)⋅(1⁢s)=1⁢rad
    θ⁢(t=3⁢s) = θ˙0⋅(3⁢s)=3⁢rad.

    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 3⁢m 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 1⁢m above the water level). If the initial vertical velocity of the diver’s center-of-mass is 3⁢m/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
−m⁢g⁢ȷˆ = m⁢y¨⁢ȷˆ
⇒ ⁢y¨ = −g
∑𝑴⇀cm = 𝑯⇀˙cm
𝟎⇀ = Iz⁢zcm⁢θ¨⁢𝒌ˆ
⇒ ⁢θ¨ = 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˙0⁢t−12⁢g⁢t2
θ⁢(t) = θ0+θ˙0⁢t.

We need to find the initial angular speed θ˙0 such that θ=π when y=1⁢m (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=1⁢m at that instant, and that y0=3⁢m and y˙0=3⁢m/s. We have,

y = y0+y˙0−12⁢g⁢t2
⇒ ⁢t = y˙0±y˙02+2⁢g⁢(y0−y)g
= 3⁢m/s±(3⁢m/s)2+2⋅9.8⁢m/s2⋅(3⁢m−1⁢m)9.8⁢m/s2
= 1.15⁢ or −0.53⁢s.

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

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

Answer: θ˙0=2.73⁢rad/s

SAMPLE 17.8

Filename:sfig8-4-4c
Figure 17.40:

A plate tumbling in space. A rectangular plate of mass m=0.5⁢kg,Iz⁢zcm=2.08×10−3⁢kg⋅m2, and dimensions a=200⁢mm and b=100⁢mm is pushed by a force 𝑭⇀=0.5⁢N⁢ıˆ, acting at d=30⁢mm 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.5⁢N⁢ıˆ0.5⁢kg=1⁢m/s2⁢ıˆ.

    which is the initial acceleration of the mass-center.

    Answer: a⇀cm=1⁢m/s2⁢ıˆ

  2. 2.

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

    𝑴⇀cm = 𝑯⇀˙cm
    F⁢d⁢𝒌ˆ = Iz⁢zcm⁢𝝎⇀˙
    ⇒ ⁢𝝎⇀˙ = F⁢dIz⁢zcm⁢𝒌ˆ=0.5⁢N⋅0.03⁢m2.08⁢kg⋅m2=7.2⁢rad/s2⁢𝒌ˆ

    which is the initial angular acceleration of the plate.

    Answer: ω⇀˙=7.2⁢rad/s2⁢kˆ

  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, 𝝎⇀˙=F⁢d/Iz⁢zcm⁢𝒌ˆ, 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/c⁢m which makes a fixed angle α(=tan−1⁡db/2) with the transverse axis of the plate. This angle is shown here merely for ease of calculation.
    𝑴⇀cm = 𝑯⇀˙cm
    𝒓⇀P/cm×𝑭⇀ = Iz⁢zcm⁢θ¨⁢𝒌ˆ.
      Now,
    𝒓⇀P/cm = −r⁢[cos⁡(θ+α)⁢ıˆ+sin⁡(θ+α)⁢ȷˆ]
    𝑭⇀ = F⁢ıˆ
    ⇒ ⁢𝒓⇀P/cm×𝑭⇀ = F⁢r⁢sin⁡(θ+α)⁢𝒌ˆ.
      Thus,
    θ¨ = F⁢rIz⁢zcm⁢sin⁡(θ+α)

    where r=d2+(b/2)2 and α=tan−1⁡(2⁢d/b).

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

    Answer: x¨=Fm,y¨=0,θ¨=F⁢rIz⁢zcm⁢sin⁡(θ+α)

  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˙0⁢t+12⁢Fm⁢t2
    y⁢(t) = y0+y˙0⁢t.

    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)=F2⁢m⁢t2,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=𝒓⇀cm−r⁢cos⁡(θ+α)⁢ıˆ−r⁢sin⁡(θ+α)⁢ȷˆ. 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∑𝑭⇀⁢d⁢t
    m⁢𝒗⇀cm⁢(t)−m⁢𝒗⇀cm⁢(0) = ∫0t𝑭⇀⁢𝑑t
    ⇒ ⁢𝒗⇀cm⁢(t) = 𝒗⇀cm⁢(0)+1m⁢∫02𝑭⇀⁢𝑑t
    = 𝟎⇀+10.5⁢kg⁢∫02(0.5⁢N⁢ıˆ)⁢𝑑t
    = 2⁢m/s⁢ıˆ.

    Answer: v⇀cm⁢(2⁢s)=2⁢m/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∑𝑴⇀cm⁢d⁢t
    Iz⁢zcm⁢𝝎⇀⁢(t)−Iz⁢zcm⁢𝝎⇀⁢(0) = ∫0t∑𝑴⇀cm⁢d⁢t
    ⇒ ⁢𝝎⇀⁢(t) = 𝝎⇀⁢(0)+1Iz⁢zcm⁢∫0t∑𝑴⇀cm⁢d⁢t
    = 𝝎⇀⁢(0)+1Iz⁢zcm⁢∫0t(𝒓⇀P/cm×𝑭⇀)⁢𝑑t
    = 𝟎⇀+1Iz⁢zcm⁢∫0t(F⁢r⁢sin⁡(θ+α)⁢𝒌ˆ)⁢𝑑t
    = F⁢rIz⁢zcm⁢(∫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=1⁢m, length b=2⁢m, and mass m=1⁢kg in the figure is sliding on the x⁢y-plane with no friction. At the moment in question, point C is at xC=3⁢m and yC=2⁢m. The linear momentum is 𝑳⇀=4⁢ıˆ+3⁢ȷˆ⁢(kg⁢⋅⁢m/s) and the angular momentum about the center of mass is 𝑯⇀c⁢m=5⁢𝒌ˆ⁢(kg⁢⋅⁢m/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¨c⁢m? What is θ¨A⁢B? 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: a⇀c⁢m=Fm⁢ıˆ.

  2. (b)

    What is the angular acceleration of the stick? Answer: α⇀=−6⁢Fm⁢ℓ⁢kˆ.

  3. (c)

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

  4. (d)

    (relatively harder) What additional force would have to be applied to point B to make point B’s acceleration zero? Answer: F⇀B=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 x⁢y plane. A massless peg is attached to a point on its perimeter. This disk has radius of 1⁢m and mass of 10⁢kg. A constant force of F=1000⁢N⁢ı is applied to the peg for .0001⁢s (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: α⇀=F⁢dm⁢R22⁢kˆ.

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 𝑭⇀=1⁢N⁢ıˆ 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=1⁢s. 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.0⁢ft 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.0⁢ft/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 x⁢y-plane. What is the total kinetic energy of the disk?

17.2.13

Calculate the energy stored in a spring using the expression EP=12⁢k⁢δ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=50⁢N and has speed v=2⁢m/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 Mi⁢n. If the driven gear rotates at a constant speed ωo⁢u⁢t, 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 Vi⁢n=c⁢o⁢n⁢s⁢t. Assuming that Vo⁢u⁢t=7⁢Vi⁢n and the force between the left belt and roller is Fi⁢n=3⁢lb:

  1. (a)

    What is Fo⁢u⁢t (draw a picture defining the signs of Fi⁢n and Fo⁢u⁢t)?. Answer: Fo⁢u⁢t=37⁢lb.

  2. (b)

    Is Fo⁢u⁢t greater or less than the Fi⁢n? (Assume Fi⁢n>0.) Why? Answer: Fo⁢u⁢t is always less than the Fi⁢n.

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ϕ)ıˆ+(1−cosϕ)ȷˆ)]
⇒𝒗⇀D = R⁢[(ϕ˙⁢(1−cos⁡ϕ)⁢ıˆ+ϕ˙⁢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[ϕ˙(1−cos(2nπ⏞1)⏟0]ıˆ+ϕ˙sin⁡(2⁢n⁢π)⏟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.

𝒗⇀A⏞Velocity of point on the wheel touching the ground=𝒗⇀B⏞Velocity 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
⇒𝒗⇀C⏟vC⁢ıˆ+𝝎⏟⇀ω⁢𝒌ˆ×𝒓⇀A/C⏟−R⁢ȷˆ = 𝟎⇀
⇒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=𝒗⇀circle⁢center+𝝎⇀×𝒓⇀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,

𝒗⇀C⏞A point on the spider = 𝒗⇀C⏞A 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
ω𝒫⁢(2⁢RP)⁢𝒆ˆθ = ωℛ⁢rR⁢𝒆ˆθ
⇒ωℛ = 2⁢RPrR⁢ω𝒫⏟ω𝒫=rCRP⁢ω𝒮
ωℛ = 2⁢RPrR⁢rC⏞rC=rS+RPRP⁢ω𝒮
= 2⁢(rS+RP)rR⏟rR=rS+2⁢RP⁢ω𝒮
⇒ωℛω𝒮 = 2⁢1+RPrS1+2⁢RPrS= 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=3⁢m slip along the frictionless wall and floor shown in Figure 17.59. At the instant shown, when θ=60∘, the angular speed θ˙=1.15⁢rad/s and the angular acceleration θ¨=2.5⁢rad/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.15⁢rad/s⁢𝒌ˆ
𝝎⇀˙ = θ¨⁢𝒌ˆ=−2.5⁢rad/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,𝒂⇀A⁢and⁢𝒂⇀B.

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

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

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

vA⁢ȷˆ⁢⋅⁢ıˆ⏟0 = (vB+θ˙⁢L⁢sin⁡θ)⁢ıˆ⁢⋅⁢ıˆ⏟1+θ˙⁢L⁢cos⁡θ⁢ȷˆ⁢⋅⁢ıˆ⏟0
⇒ ⁢0 = vB+θ˙⁢L⁢sin⁡θ
⇒ ⁢vB = −θ˙⁢L⁢sin⁡θ=−(−1.15⁢rad/s)⁢⋅⁢3⁢m⁢⋅⁢32
= 2.99⁢m/s.

Answer: 𝒗⇀B=2.9⁢m/s⁢ıˆ

Similarly,

𝒂⇀A = 𝒂⇀B+𝝎⇀˙×𝒓⇀A/B+𝝎⇀×(𝝎⇀×𝒓⇀A/B)⏞−ω2⁢𝒓⇀A/B
aA⁢ȷˆ = aB⁢ıˆ+θ¨⁢𝒌ˆ×L⁢(−cos⁡θ⁢ıˆ−sin⁡θ⁢ȷˆ)−θ˙2⁢L⁢(−cos⁡θ⁢ıˆ−sin⁡θ⁢ȷˆ)
= (aB+θ¨⁢L⁢sin⁡θ+θ˙2⁢L⁢cos⁡θ)⁢ıˆ+(−θ¨⁢L⁢cos⁡θ+θ˙2⁢L⁢sin⁡θ)⁢ȷˆ.

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

0 = aB+θ¨⁢L⁢sin⁡θ+θ˙2⁢L⁢cos⁡θ
⇒ ⁢aB = −θ¨⁢L⁢sin⁡θ−θ˙2⁢L⁢cos⁡θ
= −(−2.5⁢rad/s2⁢⋅⁢3⁢m⁢⋅⁢32)−(−1.15⁢rad/s)2⁢⋅⁢3⁢m⁢⋅⁢12
= 4.51⁢m/s2.

Answer: 𝒂⇀B=4.51⁢m/s2

SAMPLE 17.11

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

A cylinder of diameter 500⁢mm rolls down an inclined plane with uniform acceleration (of the center-of-mass) a=0.1⁢m/s2. At an instant t0, the mass-center has speed v0=0.5⁢m/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.25⁢m.

  1. 1.

    From the kinematics of pure rolling,

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

    Answer: ω=2⁢rad/s,ω˙=0.4⁢rad/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)=ω=2⁢rad/s and θ⁢(t=0)=0, to get

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

    Answer: θ=0.76⁢rev

  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.25⁢m)⋅(4.8)=1.2⁢m.

    Answer: x=1.2⁢m

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

SAMPLE 17.12  Condition of pure rolling. A cylinder of radius R=20⁢cm rolls on a flat surface with absolute angular speed ω=12⁢rad/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⁢ıˆ=(12⁢rad/s)⋅(0.2⁢m)⁢ıˆ=2.4⁢m/s⁢ıˆ.
Case(ii)

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

  1. 1.

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

  2. 2.

    Velocity of the center:

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

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

  1. 1.

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

  2. 2.

    Velocity of the center:

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

Answer: (a):(i)⁢𝒗⇀P=𝟎⇀,(ii)⁢𝒗⇀P=1⁢m/s⁢ıˆ,(iii)⁢𝒗⇀P=−1⁢m/s⁢ıˆ,(b):(i)⁢𝒗⇀C=2.4⁢m/s⁢ıˆ,(ii)⁢𝒗⇀C=3.4⁢m/s⁢ıˆ,(iii)⁢𝒗⇀C=1.4⁢m/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=0⁢s 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=0⁢s. 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=(R−r)⁢𝝀ˆO⁢C+ℓ⁢𝝀ˆC⁢P

    where

    𝝀ˆO⁢C = a unit vector along OC =−sin⁡ϕ⁢ıˆ−cos⁡ϕ⁢ȷˆ,
    𝝀ˆC⁢P = a unit vector along CP =−sin⁡θ⁢ıˆ+cos⁡θ⁢ȷˆ.

    Thus,

    𝒓⇀P=[−(R−r)⁢sin⁡ϕ−ℓ⁢sin⁡θ]⁢ıˆ+[−(R−r)⁢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⁢(−𝝀ˆO⁢C) = −ϕ˙⁢𝒌ˆ×(R−r)⁢𝝀ˆO⁢C
    or −ω⁢r⁢(𝒌ˆ×𝝀ˆO⁢C) = −ϕ˙⁢(R−r)⁢(𝒌ˆ×𝝀ˆO⁢C)
    ⇒ ⁢rR−r⁢ω = ϕ˙.

    Integrating the last expression with respect to time, we obtain

    ϕ=rR−r⁢ω⁢t.

    Let

    q=rR−r,

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

    𝒓⇀P=[−(R−r)⁢sin⁡(q⁢ω⁢t)−ℓ⁢sin⁡(ω⁢t)]⁢ıˆ+[−(R−r)⁢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=−ω⁢[(R−r)⁢q⁢cos⁡(q⁢ω⁢t)+ℓ⁢cos⁡(ω⁢t)]⁢ıˆ+ω⁢[(R−r)⁢q⁢sin⁡(q⁢ω⁢t)−ℓ⁢sin⁡(ω⁢t)]⁢ȷˆ.

    Substituting (R−r)⁢q=r in 𝒗⇀P we get

    𝒗⇀P=−ω⁢r⁢[{cos⁡(q⁢ω⁢t)+ℓr⁢cos⁡(ω⁢t)}⁢ıˆ−{sin⁡(q⁢ω⁢t)−ℓr⁢sin⁡(ω⁢t)}⁢ȷˆ]. (17.54)
  3. 3.

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

    𝒂⇀P=−ω2⁢r⁢[−{q⁢sin⁡(q⁢ω⁢t)+ℓr⁢sin⁡(ω⁢t)}⁢ıˆ−{q⁢cos⁡(q⁢ω⁢t)−ℓr⁢cos⁡(ω⁢t)}⁢ȷˆ]. (17.55)

Answer: r⇀P=[−(R−r)⁢sin⁡(q⁢ω⁢t)−ℓ⁢sin⁡(ω⁢t)]⁢ıˆ+[−(R−r)⁢cos⁡(q⁢ω⁢t)+ℓ⁢cos⁡(ω⁢t)]⁢ȷˆv⇀P=−ω⁢r⁢[{cos⁡(q⁢ω⁢t)+ℓr⁢cos⁡(ω⁢t)}⁢ıˆ−{sin⁡(q⁢ω⁢t)−ℓr⁢sin⁡(ω⁢t)}⁢ȷˆ]a⇀P=−ω2⁢r⁢[−{q⁢sin⁡(q⁢ω⁢t)+ℓr⁢sin⁡(ω⁢t)}⁢ıˆ−{q⁢cos⁡(q⁢ω⁢t)−ℓr⁢cos⁡(ω⁢t)}⁢ȷˆ]

SAMPLE 17.14  The rolling disk: instantaneous kinematics. For the rolling disk in Sample 17.63, let R=4⁢ft,r=1⁢ft 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=1⁢s, 3⁢s, and 5.25⁢s,

  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=4⁢r,

q=rR−r=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 = −3⁢r⁢[sin⁡(ω⁢t3)⁢ıˆ+cos⁡(ω⁢t3)⁢ȷˆ],
    𝒓⇀P = 𝒓⇀C+r⁢[−s⁢i⁢n⁢(ω⁢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 = ω2⁢r3⁢[sin⁡(ω⁢t3)⁢ıˆ+cos⁡(ω⁢t3)⁢ȷˆ],
    𝒂⇀P = ω2⁢r⁢[{13⁢sin⁡(ω⁢t3)+sin⁡(ω⁢t)}⁢ıˆ+{13⁢cos⁡(ω⁢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=1⁢ft, ω=π⁢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=1⁢s, 3⁢s, and 5.25⁢s.

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=8⁢r,

  2. 2.

    R=4⁢r, and

  3. 3.

    R=2⁢r.


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=8⁢r. 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=8⁢r.
  2. 2.
    Filename:pfigure4-4-rp12
    Figure 17.68: Path of point P and the center C of the disk for R=4⁢r.

    Similarly, for R=4⁢r 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=2⁢r.

    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)=ω⁢t⁢Rw−Rp⁢sin⁡(ω⁢t)
    y⁢(t)=Rw−Rp⁢cos⁡(ω⁢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 1⁤14 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=1⁢m/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=2⁢b 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⁢ω˙⁢𝒌ˆ (where⁢I=Iz⁢zcm) (17.57)
EK =mtot⁢vc⁢m2/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=Pd⁢i⁢s⁢s=Ffriction⋅vslip (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+ω˙⁢Iz⁢zcm⁢𝒌ˆ⇒𝟎⇀=R⁢ȷˆ×(−m⁢ω˙⁢R⁢ıˆ)+ω˙⁢Iz⁢zcm⁢𝒌ˆ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 𝑭⇀t⁢o⁢t=F⁢ıˆ+n⁢ȷˆ=m⁢g⁢ȷˆ. 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⁢ȷˆ−m⁢g⁢ȷˆ=m⁢a⁢ıˆ}{}⋅ȷˆ⇒N=m⁢g{}⋅ıˆ⇒a=−μ⁢gAMB/G:⇒−R⁢μ⁢m⁢g=Iz⁢zcm⁢ω˙⇒v=v0−μ⁢g⁢t and ω=−μ⁢R⁢m⁢g⁢t/Iz⁢zcm

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=−ω⁢R⇒v0−μ⁢g⁢t=−(−μ⁢R⁢m⁢g⁢t/Iz⁢zcm)⁢R⇒t=v0μ⁢g⁢(1+m⁢R2Iz⁢zcm).

Note that the energy lost during sliding is less than μ⁢m⁢g 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×(−m⁢g⁢ȷˆ)=𝒓⇀G/C×m⁢𝒂⇀G+Iz⁢zcm⁢ω˙⁢𝒌ˆ{−R⁢sin⁡ϕ⁢m⁢g⁢𝒌ˆ=(R⁢𝒏ˆ)×(m⁢aG⁢𝝀ˆ)+Iz⁢zcm⁢ω˙⁢𝒌ˆ}{}⋅𝒌ˆ⇒−R⁢m⁢g⁢sin⁡ϕ=−R⁢m⁢aG−Iz⁢zcm⁢aG/R⇒aG=g⁢sin⁡ϕ1+Iz⁢zcm/(m⁢R2).

Which is less than the acceleration of a block sliding on a ramp without friction: a=g⁢sin⁡ϕ (unless the mass of the rolling ball is concentrated at the center with Iz⁢zcm=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=dd⁢t⁢ET⇒0=EK˙+EP˙=dd⁢t⁢(m⁢v2/2+Iz⁢zcm⁢ω2/2)+dd⁢t⁢(m⁢g⁢y)=m⁢v⁢v˙+Iz⁢zcm⁢ω⁢ω˙+m⁢g⁢y˙=m⁢v⁢v˙+Iz⁢zcm⁢(v/R)⁢v˙/R−m⁢g⁢(sin⁡ϕ)⁢vassuming v≠0⇒0=(m+Iz⁢zcm/R2)⁢v˙−m⁢g⁢sin⁡ϕ⇒v˙=g⁢sin⁡ϕ1+Iz⁢zcm/(m⁢R2)

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⁢𝝀ˆ−m⁢g⁢ȷˆ=m⁢aG⁢𝝀ˆ}{}⋅𝒏ˆ⇒N=m⁢g⁢cos⁡ϕ{}⋅𝝀ˆ⇒F+m⁢g⁢sin⁡ϕ=m⁢g⁢sin⁡ϕ1+Iz⁢zcm/(m⁢R2)F=−m⁢g⁢sin⁡ϕ1+m⁢R2/Iz⁢zcmCritical condition:⇒μ=|F|N=tan⁡ϕ1+m⁢R2/Iz⁢zcm

If Iz⁢zcm 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 Iz⁢zcm≈2⁢m⁢R2/5) the steepest slope for rolling without slip is a steep ϕ=tan−1⁡(7/2)≈74∘. A metal hoop on the other hand (with μ≈.3 and Iz⁢zcm≈m⁢R2) will only roll without slip for slopes less than about ϕ=tan−1⁡(.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×(−m⁢g⁢ȷˆ)=𝒓⇀G/C×m⁢𝒂⇀G+Iz⁢zcm⁢ω˙⁢𝒌ˆr⁢m⁢g⁢sin⁡θ⁢𝒌ˆ=−r⁢𝒆ˆr×(m⁢((R−r)⁢θ¨⁢𝒆ˆθ−(R−r)⁢θ˙2⁢𝒆ˆr))+Iz⁢zcm⁢ω˙⁢𝒌ˆ.

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

(R−r)⁢θ¨=−g⁢sin⁡θ1+Iz⁢zcm/m⁢r2,

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(R−r)⁢(1+Iz⁢zcm/m⁢r2))⁢θ=0.

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

T=2⁢π⁢(R−r)⁢(1+Iz⁢zcm/m⁢r2)g.

For a marble or ball bearing in a sideways glass (with R−r≈4⁢cm=.04⁢m, Iz⁢zcm/m⁢r2≈2/5 and g≈10⁢m/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, −m⁢g⁢ȷˆ, 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 = 𝒓⇀c⁢m/C⏞R⁢ȷˆ×(F⁢𝝀ˆ)
= R⁢ȷˆ×F⁢(−cos⁡ϕ⁢ıˆ−sin⁡ϕ⁢ȷˆ)
= F⁢R⁢cos⁡ϕ⁢𝒌ˆ
and
𝑯⇀˙/C = 𝒓⇀cm/C×m⁢𝒂⇀cm+Iz⁢zcm⁢𝝎⇀˙
= R⁢ȷˆ×m⁢x⏟¨−ω˙⁢R⁢ıˆ+Iz⁢zcm⁢ω˙⁢𝒌ˆ
= m⁢ω˙⁢R2⁢𝒌ˆ+Iz⁢zcm⁢ω˙⁢𝒌ˆ
= (Iz⁢zcm+m⁢R2)⁢ω˙⁢𝒌ˆ.
Thus,
F⁢R⁢cos⁡ϕ⁢𝒌ˆ = (Iz⁢zcm+m⁢R2)⁢ω˙⁢𝒌ˆ
⇒ ⁢ω˙≡θ¨ = F⁢R⁢cos⁡ϕIz⁢zcm+m⁢R2

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¨=−F⁢R2⁢cos⁡ϕIz⁢zcm+m⁢R2.

Answer: θ¨=F⁢R⁢cos⁡ϕIz⁢zcm+m⁢R2

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 2⁢ft and mass 20⁢lbm rolls without slipping on a horizontal surface. The kinetic energy of the wheel is 1700⁢ft⁢⋅⁢lbf. 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 = 12⁢m⁢vcm2+12⁢Icm⁢ω2
    = 12⁢m⁢ω2⁢r2+12⁢Icm⁢ω2
    = 12⁢(m⁢r2+Icm⏟12⁢m⁢r2)⁢ω2
    = 34⁢m⁢r2⁢ω2
    ⇒ ⁢ω2 = 4⁢EK3⁢m⁢r2
    = 4×1700⁢ft⁢⋅⁢lbf3×20⁢lbm⁢⋅⁢1⁢ft2
    = 4×1700×32.2⁢lbm⁢⋅⁢ft/s23×20⁢lbm⁢⋅⁢ft
    = 3649.33⁢1s2
    ⇒ ⁢ω = 60.4⁢rad/s.

    Answer: ω=60.4⁢rad/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 = EK1−EK2Δ⁢t
    = 1700⁢ft⁢⋅⁢lbf−05⁢s=340⁢ft⁢⋅⁢lbf/s
    = 340⁢ft⁢⋅⁢lbf/s⋅1⁢ hp550⁢ft⁢⋅⁢lbf/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⏞𝟎⇀−m⁢g⁢ȷˆ⋅𝒗⇀c⁢m⏞v⁢ıˆ+F⁢𝝀ˆ⋅𝒗⇀c⁢m⏞v⁢ıˆ
= −m⁢g⁢v⁢(ıˆ⋅ȷˆ⏟0)+F⁢v⁢(𝝀ˆ⋅ıˆ⏟−cos⁡ϕ)
= −F⁢v⁢cos⁡ϕ
and
E˙K = dd⁢t⁢(12⁢m⁢x˙2+12⁢Iz⁢zcm⁢ω2⏞(x˙/R)2)
= 12⁢dd⁢t⁢[(m+Iz⁢zcmR2)⁢x˙2]
= (m+Iz⁢zcmR2)⁢x˙⁢x¨.
Thus,
−F⁢v⁢cos⁡ϕ = (m+Iz⁢zcmR2)⁢x˙⁢x¨
or −F⁢x˙⁢cos⁡ϕ = (m+Iz⁢zcmR2)⁢x˙⁢x¨
⇒ ⁢x¨ = −F⁢cos⁡ϕm+Iz⁢zcmR2.

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

θ¨=(F/R)⁢cos⁡ϕm+Iz⁢zcmR2.

Answer: x¨=−F⁢cos⁡ϕm+Iz⁢zcm/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=1⁢kg and radius R=0.4⁢m 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⁢𝒏ˆ×(−m⁢g⁢ȷˆ)
= −R⁢m⁢g⁢sin⁡α⁢𝒌ˆ
and
𝑯⇀˙/C = −Iz⁢zcm⁢ω˙⁢𝒌ˆ+𝒓⇀O/C⏞R⁢𝒏ˆ×m⁢𝒂⇀cm⏞R⁢ω˙⁢𝝀ˆ
= −Iz⁢zcm⁢ω˙⁢𝒌ˆ+m⁢R2⁢ω˙⁢(𝒏ˆ×𝝀ˆ)
= −(Iz⁢zcm+m⁢R2)⁢ω˙⁢𝒌ˆ.
Thus,
−R⁢m⁢g⁢sin⁡α⁢𝒌ˆ = −(Iz⁢zcm+m⁢R2)⁢ω˙⁢𝒌ˆ
⇒ ⁢ω˙ = g⁢sin⁡αR⁢[1+Iz⁢zcm/(m⁢R2)].

Answer: ω˙=g⁢sin⁡αR⁢[1+Iz⁢zcm/(m⁢R2)]

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 β=g⁢sin⁡α/R⁢(1+Iz⁢zcm/m⁢R2). Then,

ω ≡ θ˙=β⁢t+C1
θ = 12⁢β⁢t2+C1⁢t+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=R⋅12⁢g⁢sin⁡αR⁢(1+Iz⁢zcm/m⁢R2)⁢t2
= 12⁢g⁢sin⁡α1+12⁢m⁢R2m⁢R2⁢t2=13⁢(g⁢sin⁡α)⁢t2
x⁢(2⁢s) = 13⋅9.8⁢m/s2⋅sin⁡(30∘)⋅(2⁢s)2=6.53⁢m.

Answer: x⁢(2⁢s)=6.53⁢m

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=2⁢m) angular speeds of the disk, respectively. We know that in rolling, the kinetic energy is given by

EK=12⁢m⁢vcm2⏞(ω⁢R)2+12⁢Iz⁢zcm⁢ω2=12⁢(m⁢R2+Iz⁢zcm)⁢ω2.

Therefore,

Δ⁢EK=EK2−EK1=12⁢(m⁢R2+Iz⁢zcm)⁢(ω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=(−m⁢g⁢ȷˆ)⋅(d⁢𝝀ˆ)=−m⁢g⁢d⁢(ȷˆ⋅𝝀ˆ⏞−sin⁡α)=m⁢g⁢d⁢sin⁡α. (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

m⁢g⁢d⁢sin⁡α = 12⁢(m⁢R2+Iz⁢zcm)⁢(ω22−ω12)
⇒ ⁢ω22 = ω12+2⁢m⁢g⁢d⁢sin⁡αm⁢R2+Iz⁢zcm=ω12+2⁢g⁢d⁢sin⁡αR2⁢(1+Iz⁢zcmm⁢R2)
= ω12+4⁢g⁢d⁢sin⁡α3⁢R2.
Substituting the values of g,d,α,R, etc., and setting ω1=0, we get
ω22 = 4⋅(9.8m/s2)⋅(2m)⋅(sin(30∘)3⋅(0.4⁢m)2=81.67/s2
⇒ ⁢ω2 = 9.04⁢rad/s.

The corresponding speed of the center-of-mass is

vcm=ω2⁢R=9.04⁢rad/s⋅0.4⁢m=3.61⁢m/s.

Answer: vcm=3.61⁢m/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⁢ıˆ−m⁢g⁢ȷˆ)⁢Δ⁢t = m⁢(v2−v1)⁢𝝀ˆ. (17.62)
Dotting eqn. (17.62) with 𝝀ˆ gives
(−F−m⁢g⁢(ȷˆ⋅𝝀ˆ⏟−sin⁡α))⁢Δ⁢t = m⁢(v2−v1)
(−F+m⁢g⁢sin⁡α)⁢Δ⁢t = m⁢R⁢(ω2−ω1). (17.63)

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

(−F⁢R⁢𝒌ˆ)⁢Δ⁢t = −Iz⁢zcm⁢(ω2−ω1)⁢𝒌ˆ
⇒ ⁢F⁢R⁢Δ⁢t = Iz⁢zcm⁢(ω2−ω1). (17.64)

Note that the other forces (N and m⁢g) 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+m⁢g⁢sin⁡α)⁢Δ⁢t⋅R+F⁢R⁢Δ⁢t = m⁢R⁢(ω2−ω1)⋅R+Iz⁢zcm⁢(ω2−ω1)
or ⁢m⁢g⁢R⁢sin⁡α⁢Δ⁢t = (Iz⁢zcm+m⁢R2)⁢(ω2−ω1)
or ⁢g⁢sin⁡α⁢Δ⁢t = R⁢(1+Iz⁢zcmm⁢R2)⁢(ω2−ω1)
⇒ ⁢ω2 = ω1+g⁢sin⁡αR⁢(1+Iz⁢zcmm⁢R2)⁢Δ⁢t.

Answer: ω2=ω1+g⁢sin⁡αR⁢(1+Iz⁢zcmm⁢R2)⁢Δ⁢t


Note: From the answer obtained, we clearly see that ω2=ω1+ω˙⁢Δ⁢t where ω˙ is the angular acceleration given by the expression ω˙=g⁢sin⁡αR⁢(1+Iz⁢zcm/m⁢R2). 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=ℓ2⁢cos⁡θ⁢ and ⁢yG=ℓ2⁢sin⁡θ,

    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⁢ıˆ+(N2−m⁢g)⁢ȷˆ=m⁢𝒂⇀G=m⁢𝒓⇀¨G.

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

    𝒓⇀¨G=ℓ2⁢[(−θ¨⁢sin⁡θ−θ˙2⁢cos⁡θ)⁢ıˆ+(θ¨⁢cos⁡θ−θ˙2⁢sin⁡θ)⁢ȷˆ].

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

    N1 = −12⁢m⁢ℓ⁢(θ¨⁢sin⁡θ+θ˙2⁢cos⁡θ)
    N2 = 12⁢m⁢ℓ⁢(θ¨⁢cos⁡θ−θ˙2⁢sin⁡θ)+m⁢g.

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

    ∑𝑴⇀/G = (−N1⁢ℓ2⁢sin⁡θ+N2⁢ℓ2⁢cos⁡θ)⁢𝒌ˆ
    = 12⁢m⁢ℓ⁢(θ¨⁢sin⁡θ+θ˙2⁢cos⁡θ)⁢ℓ2⁢sin⁡θ⁢𝒌ˆ
    +[12⁢m⁢ℓ⁢(θ¨⁢cos⁡θ−θ˙2⁢sin⁡θ)+m⁢g]⁢ℓ2⁢cos⁡θ⁢𝒌ˆ
    = (14⁢m⁢ℓ2⁢θ¨+12⁢m⁢g⁢ℓ⁢cos⁡θ)⁢𝒌ˆ
    and
    𝑯⇀˙/G = Iz⁢z/G⁢𝝎⇀˙=112⁢m⁢ℓ2⁢θ¨⁢(−𝒌ˆ),

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

    14⁢m⁢ℓ2⁢θ¨+12⁢m⁢g⁢ℓ⁢cos⁡θ = −112⁢m⁢ℓ2⁢θ¨
    or ⁢(112+14)⁢ℓ2⁢θ¨ = −12⁢g⁢ℓ⁢cos⁡θ
    or ⁢θ¨ = −3⁢g2⁢ℓ⁢cos⁡θ (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⁢ωd⁢t=d⁢ωd⁢θ⋅d⁢θd⁢t=ω⁢d⁢ωd⁢θ.

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

    ∫0ωω⁢𝑑ω = −∫θ0θ3⁢g2⁢ℓ⁢cos⁡θ⁢d⁢θ
    ⇒ 12⁢ω2 = −3⁢g2⁢ℓ⁢(sin⁡θ−sin⁡θ0)
    ⇒ ω = ±3⁢gℓ⁢(sin⁡θ0−sin⁡θ).

    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.,

    𝝎⇀=θ˙⁢(−𝒌ˆ)=−3⁢gℓ⁢(sin⁡θ0−sin⁡θ)⁢𝒌ˆ.

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,and⁢m⁢g 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 = EK1⏟K.E.+EP1⏟P.E.
    = 0+m⁢g⁢h1
    = m⁢g⁢ℓ2⁢sin⁡θ0
    andE2 = EK2+EP2=12⁢m⁢vG2+12⁢Iz⁢zG⁢ω2⏟EK2+m⁢g⁢h2.

    Equating E1 and E2 we get

    m⁢g⁢ℓ2⁢(sin⁡θ0−sin⁡θf) = 12⁢(m⁢vG2+112⁢m⁢ℓ2⏟Iz⁢zG⁢ω2)
    org⁢ℓ⁢(sin⁡θ0−sin⁡θf) = vG2+112⁢ℓ2⁢ω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⁢ω
    ⇒ ⁢ω = 2⁢vGℓ.

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

    g⁢ℓ⁢(sin⁡θ0−sin⁡θf) = vG2+112⁢ℓ2⁢⋅⁢4⁢vG2ℓ2
    = 43⁢vG2
    ⇒ ⁢vG = 3⁢g⁢ℓ4⁢(sin⁡θ0−sin⁡θf)
    = 0.46⁢g⁢ℓ.

    Answer: vG=0.46⁢g⁢ℓ

  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 m⁢g. This force acts at the center-of-mass G. Therefore,

    P = 𝑭⇀⁢⋅⁢𝒗⇀=−m⁢g⁢ȷˆ⁢⋅⁢𝒗⇀G
    = −m⁢g⁢ȷˆ⁢⋅⁢[ℓ2⁢(−sin⁡θ⁢ıˆ+cos⁡θ⁢ȷˆ)⁢θ˙]
    = −m⁢g⁢ℓ2⁢θ˙⁢cos⁡θ.
      Now, the rate of change of kinetic energy is
    E˙K = dd⁢t⁢(12⁢m⁢vG2+12⁢Iz⁢zG⁢ω2)
    = dd⁢t⁢(12⁢m⁢ℓ2⁢ω24+12⁢m⁢ℓ212⁢ω2)
    = m⁢ℓ24⁢ω⁢ω˙+m⁢ℓ212⁢ω⁢ω˙
    = m⁢ℓ23⁢ω⁢ω˙≡m⁢ℓ23⁢θ˙⁢θ¨(since ω=θ˙ and ω˙=θ¨).

    Now equating P and E˙K we get

    m⁢ℓ23⁢θ˙⁢θ¨ = −m⁢g⁢ℓ2⁢θ˙⁢cos⁡θ
    ⇒ ⁢θ¨ = −3⁢g2⁢ℓ⁢cos⁡θ

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

    Answer: θ¨=−3⁢g2⁢ℓ⁢cos⁡θ

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=1⁢kg, radius r=10⁢cm and two identical outer disks of mass M=2⁢kg 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+2⁢M)⁢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+2⁢M)⁢g⁢ȷˆ+N⁢𝒏ˆ−Ff⁢𝝀ˆ = (m+2⁢M)⁢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+2⁢M)⁢g⁢(−ȷˆ)
= r⁢𝒏ˆ×(m+2⁢M)⁢g⁢(−ȷˆ)
= −(m+2⁢M)⁢g⁢r⁢sin⁡θ⁢𝒌ˆ(see Fig. 17.94)
and
𝑯⇀˙/C = Iz⁢zG⁢𝝎⇀˙+𝒓⇀G/C×mtotal⁢𝒂⇀G
= Iz⁢zG⁢(−ω˙⁢𝒌ˆ)−mtotal⁢ω˙⁢r2⁢𝒌ˆ
= (Iz⁢zG+mtotal⁢r2)⁢(−ω˙⁢𝒌ˆ)
= [(12⁢m⁢r2+2⋅12⁢M⁢R2)+(m+2⁢M)⏞mtotal⁢r2]⁢(−ω˙⁢𝒌ˆ)
= −[32⁢m⁢r2+M⁢(R2+2⁢r2)]⁢ω˙⁢𝒌ˆ.

Thus,

−(m+2⁢M)⁢g⁢r⁢sin⁡θ⁢𝒌ˆ = −[32⁢m⁢r2+M⁢(R2+2⁢r2)]⁢ω˙⁢𝒌ˆ
⇒ ⁢ω˙ = (m+2⁢M)⁢g⁢r⁢sin⁡θ32⁢m⁢r2+M⁢(R2+2⁢r2). (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+2⁢M)⁢g⁢r2⁢sin⁡θ32⁢m⁢r2+M⁢(R2+2⁢r2).

Now we can solve for R in terms of aG:

32⁢m⁢r2+M⁢(R2+2⁢r2) = (m+2⁢M)⁢g⁢r2⁢sin⁡θaG
⇒ ⁢M⁢(R2+2⁢r2) = (m+2⁢M)⁢gaG⁢r2⁢sin⁡θ−32⁢m⁢r2
⇒ ⁢R2 = (m+2⁢M)⁢gM⁢aG⁢r2⁢sin⁡θ−3⁢m2⁢M⁢r2−2⁢r2.

Since we require aG≤0.2⁢g we get

R2 ≥ ((m+2⁢M)⁢gM⋅0.2⁢g⁢sin⁡θ−3⁢m2⁢M−2)⁢r2
≥ (5⁢kg0.4⁢kg⋅12−3⁢kg4⁢kg−2)⁢(0.1⁢m)2
≥ 0.035⁢m2
⇒ ⁢R ≥ 0.187⁢m.

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

Answer: R≥18.7⁢cm

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

m⁢v˙⁢ıˆ = −μ⁢N⁢ıˆ+(N−m⁢g)⁢ȷˆ (17.69)
eqn. (17.69)⋅ȷˆ⁢ ⇒ ⁢N = m⁢g
eqn. (17.69)⋅ıˆ⁢ ⇒ ⁢m⁢v˙ = −μ⁢N=−μ⁢m⁢g
⇒ ⁢v˙ = −μ⁢g
⇒ ⁢v = v0−μ⁢g⁢t. (17.70)

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

−Iz⁢zcm⁢ω˙⁢𝒌ˆ = −μ⁢N⁢r⁢𝒌ˆ=−μ⁢m⁢g⁢r⁢𝒌ˆ
⇒ ⁢ω˙ = μ⁢m⁢g⁢rIz⁢zcm
⇒ ⁢ω = ω0⏟0+μ⁢m⁢g⁢rIz⁢zcm⁢t. (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−μ⁢g⁢t∗ = μ⁢m⁢g⁢r2Iz⁢zcm⁢t∗
⇒ ⁢v0 = μ⁢g⁢t∗⁢(1+m⁢r2Iz⁢zcm)
⇒ ⁢t∗ = v0μ⁢g⁢(1+m⁢r2Iz⁢zcm).

Now, for a sphere, Iz⁢zcm=25⁢m⁢r2. Therefore,

t∗=v0μ⁢g⁢(1+m⁢r225⁢m⁢r2)=2⁢v07⁢μ⁢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∗=2⁢v07⁢μ⁢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=10⁢cm and mass m=1⁢kg is launched horizontally with initial velocity v0=5⁢m/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=5⁢m/s, and ω0=0. We also know that at t=t∗, vcm=vt∗=ωt∗⁢r, 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∗ = m⁢vt∗⁢ıˆ−m⁢v0⁢ıˆ
(−μ⁢N⁢ıˆ+(N−m⁢g)⁢ȷˆ)⁢t∗ = m⁢(vt∗−v0)⁢ıˆ.
Dotting the above equation with ȷˆ and ıˆ, respectively, we get
N = m⁢g
−μ⁢N⏟m⁢g⁢t∗ = m⁢(vt∗−v0)
⇒ −μ⁢g⁢t∗ = vt∗−v0. (17.72)

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

∑𝑴⇀cm⁢t∗ = (𝑯⇀cm)t∗−(𝑯⇀cm)0
(−μ⁢N⁢r⁢𝒌ˆ)⁢t∗ = (Iz⁢zcm⁢ωt∗−Iz⁢zcm⁢ω0⏟0)⁢(−𝒌ˆ)
or −μ⁢m⁢g⁢r⁢t∗ = −Iz⁢zcm⁢ωt∗
⇒ ⁢ωt∗ = μ⁢m⁢g⁢r⁢t∗/Iz⁢zcm
⇒ ⁢vt∗ ≡ ωt∗⁢r=μ⁢m⁢g⁢r2⁢t∗/Iz⁢zcm.

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

−μ⁢g⁢t∗ = μ⁢m⁢g⁢r2⁢t∗/Iz⁢zcm−v0
⇒ ⁢t∗ = v0μ⁢g⁢(1+m⁢r2Iz⁢zcm)

which is, of course, the same expression we obtained for t∗ in Sample 17.94. Again, noting that Iz⁢zcm=25⁢m⁢r2 for a sphere, we calculate the time of transition as

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

Answer: t∗=0.73⁢s

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, 𝜶⇀D⁢i⁢s⁢k=−θ¨⁢𝒌ˆ? Make any reasonable assumptions you need that are consistent with the figure information and the laws of mechanics. State your assumptions. Answer: α⇀D⁢i⁢s⁢k=−43⁢rad/s2⁢kˆ.

  2. (b)

    Find the acceleration of the point A in the figure. Answer: a⇀A=83⁢m/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)⁢(1−Ri/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 Ef⁢r⁢i⁢c=m⁢R2⁢θ˙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 200⁢kg is carried on a conveyer belt which moves at V0=3⁢m/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=2⁢m/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=2⁢2⁢g⁢R3.

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 Iz⁢zc⁢m=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=12⁢a⁢t2”]

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) 𝒂⇀c⁢m=g⁢sin⁡θ⁢ıˆ, (b) 𝒂⇀c⁢m=23⁢g⁢sin⁡θ⁢ıˆ, (c) 𝒂⇀c⁢m=12⁢g⁢sin⁡θ⁢ıˆ, (d) 𝒂⇀c⁢m=23⁢g⁢sin⁡θ⁢ıˆ. 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=2⁢R0).

  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=2⁢M0).

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×(−𝑷⇀)=Iz⁢zcm1⁢(ω1+−ω1−)⁢𝒌ˆ𝒓⇀C/G2×𝑷⇀=Iz⁢zcm2⁢(ω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

𝒗⇀C⁢1+=𝒗⇀C⁢2+

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−+Iz⁢zcm⁢ω−⁢𝒌ˆ = 𝒓⇀G/A×m⁢𝒗⇀G++Iz⁢zcm⁢ω+⁢𝒌ˆ

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+=ω+⁢Iz⁢zA⁢𝒌ˆ. So one way of expressing the post-collision state is as

ω+=(𝒓⇀G/A×m⁢𝒗⇀G−+Iz⁢zcm⁢ω−⁢𝒌ˆ)⋅𝒌ˆIz⁢zA⁢ 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+−m⁢v⁢ıˆ
P⁢𝒏ˆ = m⁢𝒗⇀G2+−𝟎⇀

gives 4 scalar equations which are supplemented by the restitution equation (using e=1)

(Δ⁢𝒗⇀+)⋅𝒏ˆ=−e⁢(Δ⁢𝒗⇀−)⋅𝒏ˆ
⇒−v⁢cos⁡θ=𝒗⇀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+ = v⁢sin⁡θ⁢(sin⁡θ⁢ıˆ+cos⁡θ⁢ȷˆ),
𝒗⇀G2+ = v⁢cos⁡θ⁢(cos⁡θ⁢ıˆ−sin⁡θ⁢ȷˆ),and
P = m⁢v⁢cos⁡θ.

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=−v0⁢sin⁡θ. 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⁢𝒗⇀1⏟0+m2⁢𝒗⇀2 = m1⁢𝒗⇀1++m2⁢𝒗⇀2+
    or −m2⁢v0⁢ȷˆ = m1⁢v1+⁢ıˆ+m2⁢(v2t+⁢𝒆ˆt+v2n+⁢𝒆ˆn). (17.74)

    Dotting both sides of eqn. (17.74) with 𝒆ˆn gives

    −m2⁢v0⁢(𝒆ˆn⋅ȷˆ⏞cos⁡θ)=m1⁢v1+⁢(𝒆ˆn⋅ıˆ⏞−sin⁡θ)+m2⁢v2t+⁢(𝒆ˆn⋅𝒆ˆt⏞0)+m2⁢v2n+⁢(𝒆ˆn⋅𝒆ˆn⏞1)
    or −m2⁢v0⁢cos⁡θ=−m1⁢v1+⁢sin⁡θ+m2⁢v2n+. (17.75)

    Dotting both sides of eqn. (17.74) with 𝒆ˆt gives

    −m2⁢v0⁢(𝒆ˆt⋅ȷˆ⏞sin⁡θ)=m1⁢v1+⁢(𝒆ˆt⋅ıˆ⏞cos⁡θ)+m2⁢v2t+⁢(𝒆ˆt⋅𝒆ˆt⏞1)+m2⁢v2n+⁢(𝒆ˆt⋅𝒆ˆn⏞0)
    or −m2⁢v0⁢sin⁡θ=m1⁢v1+⁢cos⁡θ+m2⁢v2t+. (17.76)

    Answer: −m2⁢v0⁢cos⁡θ=−m1⁢v1+⁢sin⁡θ+m2⁢v2n+,−m2⁢v0⁢sin⁡θ=m1⁢v1+⁢cos⁡θ+m2⁢v2t+

  2. 2.

    Now, we rearrange eqn. (17.75) and 17.76 along with the third given equation, v+2t=−v0⁢sin⁡θ, 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.

    [−m1⁢sin⁡θ0m2−m1⁢cos⁡θm20010]⁢{v1+v2t+v2n+}={−m2⁢v0⁢cos⁡θ−m2⁢v0⁢sin⁡θ−v0⁢sin⁡θ}.

    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=10⁢mm above the center of the ball. The ball has mass m=0.2⁢kg and radius r=30⁢mm. Immediately after the strike, the center-of-mass of the ball moves with linear speed v=1⁢m/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 m⁢g,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 (m⁢g,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⁢ıˆ)⁢d⁢t=d⁢𝑳⇀⁢ ⇒ ⁢∫(F⁢ıˆ)⁢d⁢t=𝑳⇀2−𝑳⇀1

where L2−L1=Δ⁢𝑳⇀ 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⁢𝒗⇀1⏟0=𝟎⇀. Immediately after the strike, 𝒗⇀=v⁢ıˆ=1⁢m/s.

Thus𝑳⇀2=m⁢𝒗⇀=0.2⁢kg⁢⋅⁢1⁢m/s⁢ıˆ=0.2⁢N⁢⋅⁢s⁢ıˆ.
Hence∫(F⁢ıˆ)⁢d⁢t=0.2⁢N⁢⋅⁢s⁢ıˆor∫F⁢d⁢t=0.2⁢N⁢⋅⁢s. (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⁢ ⇒ ⁢∫∑𝑴⇀cm⁢d⁢t=∫d⁢𝑯⇀cm=(𝑯⇀cm)2−(𝑯⇀cm)1.
Because𝑯⇀cm = Ic⁢mz⁢z⁢𝝎⇀ and just before the strike, ⁢𝝎⇀=𝟎⇀,
(𝑯⇀cm)1 ≡ angular momentum just before the strike =𝟎⇀
(𝑯⇀cm)2 ≡ angular momentum just after the strike =Ic⁢mz⁢z⁢ω⁢𝒌ˆ,
∫∑𝑴⇀cm⁢d⁢t = Ic⁢mz⁢z⁢ω⁢𝒌ˆ=25⁢m⁢r2⁢ω⁢𝒌ˆ⁢ (since for a sphere, Ic⁢mz⁢z=25⁢m⁢r2).
But ∑𝑴⇀cm=−F⁢d⁢𝒌ˆ,
therefore −∫(F⁢d)⁢d⁢t⁢𝒌ˆ=25⁢m⁢r2⁢ω⁢𝒌ˆ
or −d⏟constant⁢∫F⁢d⁢t=25⁢m⁢r2⁢ω⁢ ⇒ ⁢ω=−5⁢d2⁢m⁢r2⁢∫F⁢d⁢t.

Substituting the given values and ∫F⁢d⁢t=0.2⁢N⁢⋅⁢s from equation 17.77 we get

ω=−5⁢(0.01⁢m)2⁢⋅⁢0.2⁢kg⁢⋅⁢(0.03⁢m)2⁢⋅⁢0.2⁢N⁢⋅⁢s=−27.78⁢rad/s.

The negative value makes sense because the ball will spin clockwise after the strike, but we assumed that ω was anticlockwise.

Answer: ω=−27.78⁢rad/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+
    Iz⁢zcm⁢𝝎⇀−+𝒓⇀G/C×m⁢𝒗⇀G− = Iz⁢zcm⁢𝝎⇀++𝒓⇀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+=Iz⁢zC⁢𝝎⇀=13⁢m⁢ℓ2⏟Iz⁢zC⁢ω⁢(−𝒌ˆ).
    𝑯⇀C− = 𝒓⇀G/C×m⁢𝒗⇀G−=(ℓ/2)⁢𝝀ˆ×m⁢v⁢(−ȷˆ)
    = −m⁢v⁢ℓ2⁢cos⁡θ⁢𝒌ˆ(since ⁢𝝀ˆ=cos⁡θ⁢ıˆ+sin⁡θ⁢ȷˆ)
    𝑯⇀C+ = Iz⁢zcm⁢𝝎⇀++𝒓⇀G/C×m⁢(𝝎⇀+×𝒓⇀G/C)
    = −Iz⁢zcm⁢ω⁢𝒌ˆ+(ℓ/2)⁢𝝀ˆ×m⁢(−ω⁢𝒌ˆ×ℓ/2⁢𝝀ˆ)⏟ω⁢ℓ/2⁢(−𝒏ˆ)
    = −112⁢m⁢ℓ2⁢ω⁢𝒌ˆ−14⁢m⁢ℓ2⁢ω⁢𝒌ˆ=−13⁢m⁢ℓ2⁢ω⁢𝒌ˆ.

    Now, equating 𝑯⇀C− and 𝑯⇀C+ we get

    ω=3⁢v2⁢ℓ⁢cos⁡θ, ⇒ ⁢𝝎⇀=−3⁢v2⁢ℓ⁢cos⁡θ⁢𝒌ˆ.

    Answer: ω⇀=−3⁢v2⁢ℓ⁢cos⁡θ⁢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=−ω⁢𝒌ˆ×ℓ⁢𝝀ˆ
    = −ω⁢ℓ⁢𝒏ˆ=−3⁢v2⁢cos⁡θ⁢(−sin⁡θ⁢ıˆ+cos⁡θ⁢ȷˆ)

    but, for small θ, cos⁡θ≈1, and sin⁡θ≈0. Therefore, 𝒗⇀B/C=−3⁢v2⁢ȷˆ.

    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: v⇀B/C=−(3/2)⁢v⁢ȷˆ

SAMPLE 17.30

Filename:sfig4-7-2a
Figure 17.127:

Tipping box. A box of mass m and dimensions 2⁢a and 2⁢b 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⁢ȷˆ)×m⁢v⁢ıˆ=−m⁢a⁢v⁢𝒌ˆ

Let the box have angular velocity 𝝎⇀+=ω⁢𝒌ˆ just after impact. Then,

𝑯⇀A+ = Iz⁢zcm⁢𝝎⇀++𝒓⇀G/A×m⁢𝒗⇀G+=Iz⁢zcm⁢ω⁢𝒌ˆ+r⁢𝝀ˆ×m⁢(ω⁢𝒌ˆ×r⁢𝝀ˆ)
= Iz⁢zcm⁢ω⁢𝒌ˆ+m⁢r2⁢ω⁢𝒌ˆ=112⁢(4⁢a2+4⁢b2)⁢m⁢ω⁢𝒌ˆ+m⁢(a2+b2)⁢ω⁢𝒌ˆ
= 43⁢(a2+b2)⁢m⁢ω⁢𝒌ˆ.

Now equating the two momenta, we get

ω+=−3⁢a4⁢(a2+b2)⁢v⁢ ⇒ ⁢𝝎⇀+=−3⁢a4⁢(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 Iz⁢zA=43⁢m⁢(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
EPt−EP+ =−(EKt−EK+)
m⁢g⁢((b2+a2)−b) =−(0−12⁢(ω+)2⁢Iz⁢zA)
(ω+)2 =32⁢g⁢(1−b(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: v≤2a⁢2⁢g⁢(b2+a2)3/2⁢(1−b(b2+a2))/3.

SAMPLE 17.31

Filename:pfigure-s94q7p1
Figure 17.130:

Ball hits the bat. A uniform bar of mass m2=1⁢kg and length 2⁢ℓ=1⁢m hangs vertically from a hinge at A. A ball of mass m1=0.25⁢kg comes and hits the bar horizontally at point D with speed v=5⁢m/s. The point of impact D is located at d=0.75⁢m 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⁢𝒗⇀−+Iz⁢zA⁢𝝎⇀−
    = d⁢ȷˆ×m1⁢v⁢(−ıˆ)+𝟎⇀=m1⁢d⁢v⁢𝒌ˆ.
      Similarly,
    𝑯⇀A+ = 𝒓⇀D/A×m1⁢𝒗⇀++Iz⁢zA⁢𝝎⇀+
      but, 𝒗⇀+=𝝎⇀+×𝒓⇀D/A=−ω⁢d⁢ıˆ, where 𝝎⇀+=ω⁢𝒌ˆ (let). Hence,
    𝑯⇀A+ = d⁢ȷˆ×m1⁢(−ω⁢d⁢ıˆ)+13⁢m2⁢(2⁢ℓ)2⁢ω⁢𝒌ˆ
    = (m1⁢d2+43⁢m2⁢ℓ2)⁢ω⁢𝒌ˆ.

    Equating the two momenta, we get

    ω = m1⁢d⁢vm1⁢d2+(4/3)⁢m2⁢ℓ2
    = vd⁢(1+43⁢m2m1⁢(ℓd)2)
    ⇒ ⁢𝒗⇀D = 𝝎⇀+×𝒓⇀D/A=ω⁢d⁢(−ıˆ)
    = −v1+43⁢m2m1⁢(ℓd)2⁢ıˆ.

    Now, substituting the given numerical values, v=5⁢m/s, m1=0.25⁢kg, m2=1⁢kg, ℓ=0.5⁢m, and d=0.75⁢m, we get 𝒗⇀D=−2.08⁢m/s⁢ıˆ

    Answer: v⇀D=−2.08⁢m/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=∫𝑭⇀D⁢d⁢t = 𝑳⇀+−𝑳⇀−=m1⁢(𝒗⇀+−𝒗⇀−)
    = m1⁢(−v1+43⁢m2m1⁢(ℓd)2⁢ıˆ+v⁢ıˆ)
    = m1⁢v⁢(1−11+43⁢m2m1⁢(ℓ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.73⁢kg⁢⋅⁢m/s⁢ıˆ. The impulse on the bar is equal and opposite. Therefore, the impulse on the bar is −𝑷⇀D=−0.73⁢kg⁢⋅⁢m/s⁢ıˆ.

    Answer: Impulse at D=−0.73⁢kg⁢⋅⁢m/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)⁢d⁢t=𝑳⇀+−𝑳⇀−=𝑳⇀+=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 ≡ ∫𝑭⇀A⁢d⁢t=∫𝑭⇀D⁢d⁢t+m2⁢(−ω⁢ℓ⁢ıˆ)
    = m1⁢v⁢(1−11+43⁢mr⁢q2)⁢ıˆ−m2⁢ℓ⁢vd⁢(1+43⁢mr⁢q2)⁢ıˆ
    = m1⁢v⁢(43⁢mr⁢q21+43⁢mr⁢q2)⁢ıˆ−m2⁢v⁢(q1+43⁢mr⁢q2)⁢ıˆ
    = (4/3)⁢m2⁢q2−m2⁢q1+43⁢mr⁢q2⁢v⁢ıˆ=q⁢(4⁢q−3)3⁢(1+43⁢mr⁢q2)⁢m2⁢v⁢ıˆ.

    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/(m1⁢v), that is, the nondimensional impulse at A, normalized with respect to the initial linear momentum m1⁢v 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 ⁢P⇀A=𝟎⇀

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 2⁢a=2⁢m and 2⁢b=1⁢m and mass mP=2⁢kg drifts in space at a uniform speed vp=10⁢m/s (relative to some local Newtonian reference frame) in the direction shown. A circular disk of radius R=0.25⁢m and mass mD=1⁢kg is heading towards the plate at a linear speed v0=vD=1⁢m/s directed normal to the facing edge of the plate. In addition, the disk is spinning at ωD=5⁢rad/s in the clockwise direction. The plate and the disk collide at point A of the plate, located at d=0.8⁢m 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 x⁢y 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⁢ıˆ−ωD⁢R⁢ȷˆ
𝒗⇀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=e⏟1⁢(−vP−vD)=−vP−vD. (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)
mD⁢vDy+ = 0. (17.81)
Dotting LMB for the plate with ıˆ and ȷˆ, respectively, we get
mP⁢(vPx++vP) = P (17.82)
mP⁢vPy+ = 0. (17.83)
Dotting AMB for the disk and the plate with 𝒌ˆ, we get
IDcm⁢(ωD++ωD) = 0 (17.84)
IPcm⁢ωP+ = −P⁢d. (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

[mD0000010mD0000000mP000−1000mP0000000IDcm0000000IPcmd10−100d0]⁢{vDx+vDy+vPx+vPy+ωD+ωP+P}={mD⁢vD0−mP⁢vP0−IDcm⁢ωD0−vP−vD}.

Substituting the given numerical values for the masses and the pre-collision velocities, and the moments of inertia, IDcm=(1/2)⁢mD⁢R2 and IPcm=(1/12)⁢mP⁢(4⁢a2+4⁢b2), and then solving the matrix equation on a computer, we get,

𝒗⇀D+ =−8.70⁢m/s⁢ıˆ,𝒗⇀P+=−5.15⁢m/s⁢ıˆ
𝝎⇀D+ =−5⁢rad/s⁢𝒌ˆ,𝝎⇀P+=−9.31⁢rad/s⁢𝒌ˆ
P =9.70⁢kg⁢⋅⁢m/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=−5⁢rad/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 10⁢m rope tied to it. A frictionless ice skater of mass 50⁢kg and speed 3⁢m/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=2⁢L/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 vh⁢i⁢t 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 vh⁢i⁢t 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).