Chapter 19 Mechanics of constrained particles and rigid objects

Download this chapter by section

All chapters and individual parts

The dynamics of particles and rigid objects is studied using the relative-motion kinematics ideas from chapter chapter 18. This is the capstone chapter for a two-dimensional dynamics course. After this chapter a good student should be able to navigate through and use most of the skills in the concept map inside the back cover.

We have studied the mechanics of particles and rigid bodies with constraints that require progressively more involved kinematics. We now proceed to study the mechanics of more complex systems: particles with constrained paths, particles moving relative to moving frames, and mechanisms with several parts.

The basic strategy throughout is to use, in combination, the following skills:

  • 1.

    Basic modeling. Describe a system in an appropriate way using the language of particle and rigid body mechanics. As described in Chapter 2, the force modeling and kinematic modeling are coupled. Where relative motion is freely allowed there is no force. And where motion is caused or prevented there is a force. Here is where you decide the constitutive (force) laws you are using for springs, contact, gravity, etc.

  • 2.

    Draw free-body diagrams of the system of interest and of its parts. These diagrams should show what you do and do not know about the constraint forces (e.g., at a pin connection cut free in a free-body diagram the FBD should show an arbitrary force and no moment). These are exactly the same free-body diagrams that one would draw for statics.

  • 3.

    Kinematics calculations. Pick appropriate configuration variables, as many as there are degrees of freedom. Then write the velocities, accelerations, angular velocities and angular accelerations of interest in terms of the configuration variables and their first and second time derivatives, possibly using methods from Chapter 10. Often this is the hardest part of the analysis.

  • 4.

    Use appropriate balance equations: linear momentum, angular momentum, or power balance equations.

  • 5.

    Solve the balance equations for unknown forces or accelerations of interest. Sometimes this can be done by hand by writing out components and solving simultaneous equations or by using appropriate dot products. And sometimes it is best done by setting up a matrix equation and solving on the computer.

  • 6.

    Solve the differential equations to find how the basic configuration variables change with time. For some special problems this can be done by hand, but most often involves computer solution.

  • 7.

    Plug the ODE solution from (6) above into the equations from kinematics (3 above) and the balance laws (4 above). This is not a different skill from (3) or (4), it is just applied at a different time in the work.

  • 8.

    Make plots of how forces, positions and velocities change with time, or of trajectories. Animations are also often nice.

These skills are used to solve dynamics problems which often fall into one of these 4 categories.

  • a.

    Kinematics. These are problems where only geometry is used, where the kinematics constraints determine what you are interested in, independent of the forces or time history. A classic example is determining the path of a point on a given four-bar linkage. More basic examples include finding position or acceleration from a given velocity history.

  • b.

    Instantaneous dynamics. These are problems where the positions and velocities of all points are given and you need to find forces or accelerations. Often these are “first-motion” problems: what are accelerations and forces immediately after something is released from rest?

  • c.

    “Inverse dynamics.” These problems are called “inverse” because they are backwards of the original hard dynamics problems ((d) below). In these problems the motion is given as a function of time, and you have to calculate the forces. These problems are easier than non “inverse” problems because the differential equations from the balance laws don’t need to be solved. A classic example is a slider-crank where the motion of the crank is known a priori to be at constant rate and you need to find the torque required to keep that motion. Usually in science “inverse” problems are harder. In dynamics this kind of “inverse” problem is easier than the non “inverse” problems.

  • d.

    Dynamics analysis. You are given some information about forces and constraints and you have to find the motion and more about the forces. These are the capstone problems that require use of all the skills.

A flow chart showing how these problem types are solved using the basic skill components ideas is shown in the chart in the inside back cover.

As you solve a problem, at any instant in time you should be able to place your work on this chart.

In the sense of putting all the basic ideas together, this chapter completes the book. But in this chapter we only consider two-dimensional models and motions. Three-dimensional models and motions involve the same assembly of basic ideas, but more difficult kinematics, so are postponed.

19.1 Mechanics of a constrained particle

The kinematics of time-varying base vectors help us deal with some more difficult particle motion mechanics problems. For one point mass it is easy to write balance of linear momentum. It is:

𝑭=m𝒂.

The mass of the particle m times its vector acceleration 𝒂 is equal to the total force on the particle 𝑭. No problem.

Now, however, we can write this equation in five somewhat distinct ways.

  1. 1.

    In general abstract vector form: 𝑭=m𝒂.

  2. 2.

    In cartesian coordinates: Fxıˆ+Fyȷˆ+Fz𝒌ˆ=m[x¨ıˆ+y¨ȷˆ+z¨𝒌ˆ].

  3. 3.

    In polar coordinates:  FR𝒆ˆR+Fθ𝒆ˆθ+Fz𝒌ˆ=m[(R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨)𝒆ˆθ+z¨𝒌ˆ].

  4. 4.

    In path coordinates:

    Ft𝒆ˆt+Fn𝒆ˆn=m[v˙𝒆ˆt+(v2/ρ)𝒆ˆn]

All of these equations are always right. Additionally, for a given particle moving under the action of a given force there are many more correct equations that can be found by shifting the origin and orientation of the coordinate systems. For example for a moving frame with origin 0, rotation rate 𝝎 and angular acceleration 𝜶:

  1. 5.

    𝑭=m{𝒂0ω2𝒓+𝜶×𝒓+𝒂/+2𝝎×𝒗/}

where, to simplify the notation, all motions are relative to and positions relative to 0 unless explicitly indicated by a / or /0.

This is quite a collection of kinematic tools. In general we want to choose the best tools for the job. But to get a sense let’s first look at a simple problem using each of these kinematic approaches, some of which are rather inappropriate.

A particle that moves with no net force

In the special case that a particle has no force on it we know intuitively, or from the verbal statement of Newton’s First Law, that the particle travels in a straight line at constant speed. As a first example, let’s try to find this result using the vector equations of motion five different ways: in the general abstract form, in cartesian coordinates, in polar coordinates, in path coordinates, and relative to a moving frame (see fig. 19.1).

Filename:tfigure8-rel-ang-vel
Figure 19.1: A particle P moves. One can track its motion using the general vector form 𝒓, Cartesian coordinates in the fixed frame =0ıˆȷˆ, Polar coordinates using 𝒆ˆR&𝒆ˆθ, path coordinates 𝒆ˆt&𝒆ˆn and cartesian coordinates in a rotating frame =0ıˆȷˆ.

General abstract form. The equation of linear momentum balance is 𝑭=m𝒂 or, if there is no force, 𝒂=𝟎, which means that d𝒗/dt=𝟎. So 𝒗 is a constant. We can call this constant 𝒗0. So after some time the particle is where it was at t=0, say, 𝒓0, plus its velocity 𝒗0 times time. That is:

𝒓=𝒓0+𝒗0t. (19.1)

This vector relation is a parametric equation for a straight line. The particle moves in a straight line, as expected.

Cartesian coordinates. If instead we break the linear momentum balance equation into cartesian coordinates we get

Fxıˆ+Fyȷˆ+Fz𝒌ˆ=m(x¨ıˆ+y¨ȷˆ+z¨𝒌ˆ).

Because the net force is zero and the net mass is not negligible,

x¨=0,y¨=0, and z¨=0.

These equations imply that x˙, y˙, and z˙ are all constants, let’s call them vx0, vy0, vz0. So x, y, and z are given by

x=x0+vx0t,y=y0+vy0t,&z=z0+vz0t.

We can put these components into their place in vector form to get:

𝒓=xıˆ+yȷˆ+z𝒌ˆ=(x0+vx0t)ıˆ+(y0+vy0t)ȷˆ+(z0+vz0t)𝒌ˆ. (19.2)

Note that there are six free constants in this equation representing the initial position and velocity. Equation 19.2 is a cartesian representation of equation 19.1; it describes a straight line being traversed at constant rate.

Polar/cylindrical coordinates. 

When there is no force, in polar coordinates we have:

FR0𝒆ˆR+Fθ0𝒆ˆθ+Fz0𝒌ˆ=m[(R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨)𝒆ˆθ+z¨𝒌ˆ].

This vector equation leads to the following three scalar differential equations, the first two of which are coupled non-linear equations (neither can be solved without the other).

R¨Rθ˙2 = 0
2R˙θ˙+Rθ¨ = 0
z¨ = 0

A tedious calculation will show that these equations are solved by the following functions of time:

R = d2+[v0(tt0)]2 (19.3)
θ = θ0+tan1[v0(tt0)/d]
z = z0+vz0t,

where θ0, d, t0, v0, z0, and vz0 are constants. Note that, though eqn. (19.1) looks different than eqn. (19.2), there are still 6 free constants. From the physical interpretation you know that eqn. (19.1) must be the parametric equation of a straight line. And, indeed, you can verify that picking arbitrary constants and using a computer to make a polar plot of eqn. (19.1) does in fact show a straight line. From eqn. (19.1) it seems that polar coordinates’ main function is to obfuscate rather than clarify. For the simple case that a particle moves with no force at all, we have to solve non-linear differential equations whereas using cartesian coordinates we get linear equations which are easy to solve and where the solution is easy to interpret.

But, if we add a central force, a force like earth’s gravity acting on an orbiting satellite (the force on the satellite is directed towards the center of the earth), the equations become almost intolerable in cartesian coordinates. But, in polar coordinates, the solution is almost as easy (which is not all that easy for most of us) as the solution 19.1. So the classic analytic solutions of celestial mechanics are usually expressed in terms of polar coordinates.

Path coordinates. 

When there is no force, 𝑭=m𝒂 is expressed in path coordinates as

Ft0𝒆ˆt+Fn0𝒆ˆn=m(v˙𝒆ˆt+(v2/ρ)𝒆ˆnv2𝜿).

That is,

v˙=0andv2/ρ=0.

So the speed v must be constant and the radius of curvature ρ of the path infinite. That is, the particle moves at constant speed in a straight line.

Relative to a rotating reference frame Let’s look at the equations using a frame that shares an origin with but is rotating at a constant rate 𝝎=ω𝒌ˆ relative to . Thus

𝜶=𝟎 and 𝒂0/0=𝟎

and we have that 𝑭=m𝒂 is written as

𝑭 = m𝒂
𝟎 = m{𝒂0𝟎ω2𝒓+𝜶𝟎×𝒓+𝒂/+2𝝎×𝒗/} (19.4)
= ω2𝒓+𝒂/+2ω𝒌ˆ×𝒗/.

Now, using the rotating base vectors ıˆ and ȷˆ, we have that Fxıˆ+Fyȷˆ=0ıˆ+0ȷˆ and

𝒓=xıˆ+yȷˆ,𝒗=x˙ıˆ+y˙ȷˆ, and 𝒂=x¨ıˆ+y¨ȷˆ

so eqn. (19.4) can be rewritten as

𝟎 = ω2(xıˆ+yȷˆ)+(x¨ıˆ+y¨ȷˆ)+2ω𝒌ˆ×(x˙ıˆ+y˙ȷˆ)

which in turn can be broken into components and written as:

x¨ = ω2x+2ωy˙ (19.5)
y¨ = ω2y2ωx˙

which makes up a pair of second order linear differential equations. With some work, someone good at ODEs can solve this with pencil and paper. But most of us would use a computer for such a system. If, in some consistent units, we had

y(0)=0,y˙(0)=0,x(0)=0,x˙(0)=1,ω=1

then the solution turns out to be

x = tcos(ωt)
y = tsin(ωt)

as you can check by substituting into eqn. (19.1). That is, a particle which goes in a straight line away from the origin goes, in spirals as seen in the rotating framemargin: And the straight line motion could have been more complicated yet, as seen in the rotating frame, with variable rate rotation and acceleration.This is why we are lucky the earth is not spinning fast or at highly variable rates. Otherwise the motion we would see, us humans here on the rotating earth, would be particles with no force on them spiraling around all over the place. And that’s what Isaac Newton would have seen. And he would have written “A particle in motion tends to go in crazy spirals, and a particle that is initially at rest also goes bonkers” and the equations we use through out this book would have been much harder to discover. .

Constrained motion

A particle in a plane has 2 degrees of freedom. There is basically only one kind of constraint — to a path. When constrained to a path the particle has one remaining degree of freedom so its configuration can be described with one variable. For a given problem you must think about

  • What force(s) constrain the motion to the path?

  • What do you want to use for a configuration variable?

You use these ideas to

  • draw an appropriate free-body diagram, and then

  • calculate the velocity and acceleration in terms of the configuration variable and its derivatives.

After these key first steps you plug into equations of motion and solve for what you are interested in. Of course the needed math could be difficult or impossible, but the work is somewhat routine from a mechanics point of view.

Filename:tfigure8-ang-vel-ex
Figure 19.2: A point-mass bead slides on a rigid immobile frictionless wire. The free-body diagram shows that the only force on the bead is in the direction normal to the wire.

Example: Bead on frictionless wire

A bead slides on a frictionless wire with a crazy but smooth shape. No forces are applied to the bead besides the constraint force (see fig. 19.2).

Fig. 19.2b shows a free-body diagram where 𝒏ˆ is the normal to the wire at the point of interest. It doesn’t matter if you use for 𝒏ˆ=𝒆ˆn, or 𝒏ˆ= a vector always, say, to the left, just so long as you know what you mean by 𝒏ˆ. The free-body diagram shows that you know that the constraint force is normal to the path (the frictionless wire) but that you don’t know how big it is (F is an unknown scalar).

For some purposes, especially general problems like this where no specific path is given, the most appropriate configuration variable is s, the arc length along the path. If the path is given we assume we know the position at any given arc length by the functions

x(s) and y(s).

So

𝒓=𝒓(s),𝒗=s˙𝒆ˆt, and 𝒂=(s˙2/ρ)𝒆ˆn+s¨𝒆ˆt.

Now we can write linear momentum balance

𝑭 = m𝒂 (19.6)
{F𝒏ˆ = m((s˙2/ρ)𝒆ˆn+s¨𝒆ˆt)}
{}𝒆ˆt  v˙ = 0 and
{}𝒆ˆn  F = mv2/ρ.

Eqn. 19.6 tells us that the bead moves at constant speed, no matter what the shape of the wire. It also tells us that the more curved the wire, the bigger the constraint force needed to keep the bead on the wire.

Because this is a 1 DOF system, any one equation of motion should give us the result. Instead of linear momentum balance we could have used power balance to get the same result like this:

P = EK˙
 𝑭tot𝒗 = ddt{12mv2}
 F𝒆ˆnv𝒆ˆt = mvv˙
 0 = v˙.

This is natural enough. The only force on the particle is perpendicular to its motion, so does no work. So the particle must have constant kinetic energy and its speed must be constant.

The example above doesn’t seem that applicable; how often does one see beads on frictionless wires? It is a bit more useful than it appears. For example, if a car is coasting on an open road and tire resistance and sideslip can be ignored, the reasoning above shows that the car is neither speeded nor slowed by turning. Similarly, an idealization of an airplane wing is as something which only causes force perpendicular to motion. So in quick maneuvers where the gravity force is relatively small, the plane maintains its speed.

Filename:tfigure8-ang-vel
Figure 19.3: A bead slides on a frictionless wire on the curve implicitly defined by y=cs2/2, where s is arc-length along the curve measured from the origin.

Example: The hard brachistochrone problem.

Filename:tfigure8-ang-accel
Figure 19.4: The roller coaster that gets from A to B the fastest is the one with a track in the shape of the brachistochrone y=cs2/2.

Here is a puzzle proposed by Johann Bernoulli in June 1696 margin: The original brachistochrone (least time) puzzle:
“I, Johann Bernoulli, greet the most clever mathematicians in the world. Nothing is more attractive to intelligent people than an honest, challenging problem whose possible solution will bestow fame and remain as lasting monument. Following the example set by Pascal, Fermat, etc., I hope to earn the gratitude of the entire scientific community by placing before the finest mathematicians of our time a problem which will test their methods and the strength of their intellect. If someone communicates to me the solution of the proposed problem, I shall then publicly declare him worthy of praise.

Let two points A and B be given in a vertical plane. Find the curve that a point M, moving on a path AMB must follow such that, starting from A, it reaches B in the shortest time under its own gravity. Newton. Besides Bernoulli’s brother Jakob, one of the people to solve the puzzle was 55 year old Isaac Newton. He attempted to keep his solution, said to have been worked out in one evening while he also invented the calculus of variations, anonymous. But Bernoulli supposedly saw through this deception, commenting “I recognize the lion by his paw print”, which was presumably not a comment about Newton’s handwriting.
:

given that a roller coaster has to coast from rest at one place to another place that is no higher, what shape should the track be to make the trip as quick as possible.

The solution. Finding the solution, or even verifying it, is a problem in the calculus of variations, i.e., too advanced for this book. The solution turns out to be the brachistochrone curve that obeys the following relationship between arc-length s from the origin and vertical position y (drawn accurately in fig. 19.3):

y=12cs2(for|s|1/c). (19.7)

Starting at the origin this curve is close to y=cx2/2 but gets a bit higher (bigger y) because for a given value of x, s is greater than x. The curve terminates at vertical tangents at s=±1/c where y=1/2c. (see fig. 19.3). To solve the puzzle this curve is scaled (by choosing a value of c) and displaced so that it has a vertical tangent at A and also so the curve goes through B. The idea that the hard math seems to be expressing, is that the particle should first build up as much speed as it can (by going straight down) and then head off in the right direction (see fig. 19.4).

On the other hand, here is an easier problem that is a virtual setup for the techniques now at hand.

Example: The easy brachistochrone problem

How long does it take for a particle to slide back and forth on a frictionless wire with y=cs2/2 as driven by gravity? (see fig. 19.3)

Let’s use s as our configuration variable. The power balance equation is:

P = EK˙
(F𝒆ˆn)(v𝒆ˆt)=0  mgȷˆ(x˙ıˆ+y˙ȷˆ) = ddt{12mv2}
 mgy˙ = mddts˙22
 mgcss˙ = ms˙s¨
Assuming s˙0  s¨+gcs=0.

This, remarkably, is the simple harmonic oscillator equation with general solution

s=Acos(gct)+Bsin(gct).

Thus the period of oscillation (T such that gcT=2π) is

T=2πgc

which is independent of the amplitude of oscillation

margin: Actually, the amplitude can’t be arbitrarily large. The solution to the defining eqn. (19.7) only makes sense for |s|<1/c. For |s|>1/c there is no curve satisfying eqn. (19.7).

.

The key to this quick solution was using a configuration variable that made the expression for the velocity simple, and using an equation of motion that didn’t involve the unknown reaction force F which we also didn’t care about. We could have got the same equation of motion by writing 𝑭=m𝒂 and eliminated F𝒆ˆn by dotting both sides with a convenient vector orthogonal to F𝒆ˆn, say 𝒗.

The brachistochrone is a famous curve that has various interesting properties (e.g., Box 19.1 on page 19.1).

Box 19.1 Some brachistochrone curiosities

The brachistochrone is a cycloid. There is no straightforward way to draw the curve y=cs2/2 because the formula doesn’t tell you the x coordinates of the points. You could find them by integrating dx=ds2dy2 numerically or with calculus tricks. But it turns out (see below) that the curve with y=cs2/2 is described by the parametric equations

  x=r(ϕ+sinϕ)  with r=14c.
  y=r(1cosϕ).

Filename:tfigure-brachwheel


This is the path of a particle on the perimeter of a wheel that rolls against a horizontal ceiling a distance 2r above the origin, as you can verify by adding up distances in the picture above (see page 18.1). We will show below that the upside down cycloid and the curve y=cs2/2 are one and the same.

Note that the osculating circle of this cycloid at its lowest point has radius 4r=1/c, just the length of a simple pendulum that, for small oscillations, has the same frequency of oscillation as the bead on the brachistochrone. A point mass swinging on a string is like a bead on a frictionless circular wire and this, in turn, is close to the motion of a bead on a brachistochrone wire for small oscillations.

Galileo (1564-1642). Well before Bernoulli’s challenge, Galileo was interested in things rolling and sliding on ramps. He knew that the shortest distance between two points is a straight line, and had noted that a ball rolling down an appropriately curved ramp gets to its destination faster than a ball traveling the shortest route. A ball going on a straight ramp just doesn’t pick up much speed, and when it finally has its greatest speed the trip is over. Imagine sliding straight sideways; it takes forever on a straight-line route. Better, he must have reasoned, to get the ball rolling fast at the start and then go fast for most of its journey, possibly slowing at the end. Galileo thought the best shape was the bottom of a circle (or fraction thereof), which isn’t far off either in shape or concept, but isn’t quite right. Galileo was apparently obsessed with cycloids for other reasons but didn’t see their connection to this problem.

Filename:tfigure-brachpen

A constant period pendulum. For clock time keeping, a pendulum is better than a bead on a wire because the friction of sliding is avoided. Unfortunately, a simple pendulum has a period which is longer if the amplitude is bigger. Not much longer, 18% if the swinging is ±90 and only 1.7% longer if the amplitude is ±30, but enough to annoy clock designers. A bead sliding frictionlessly on the path y=s2/2 has the nice property that the period does not depend on the amplitude. But any real bead sliding on any real wire has substantial friction. So, at first blush the brachistochrone curve, despite its nice constant-period property, cannot be used to keep time.

But Huygens, one of the smart old timers, looked for a curve that, when a string wraps around it, makes the end follow the brachistochrone curve. To this day you can see fancy old clocks with this wrapping device, a solid piece with two cuspoidal shapes at the hinge of the swinging (“isochronous” or “tautochrone”) pendulum which wraps around it.

Geometry. The two key features discussed above, that the curve y=cs2/2 is a cycloid, and that a cycloid can be generated by wrapping a string around another cycloid, can be found from the geometric construction below. Two cycloids are shown, one from wheel 1 rolling under line L1 and another from wheel 2 rolling under line L2 a distance 2r below. Both wheels have radius r. Imagine that the cycloids A1M1B and A2M2B are drawn by wheels always arranged with vertically aligned rolling contact points C1 and C2 and with points M1 and M2 initially aligned vertically a distance 4r apart at A1 and A2.

The two cycloids are thus the same shape but are displaced with one being 2r below and πr sideways from the other.

Because both wheels have rolled the same distance (A1C1) they have rotated the same amount and M2 is as far forward of C1C2D as M1 is behind. Similarly M1 is as far above L2 as M2 is below. So the line M1M2 is bisected by the point C2.

Because C1C2 is the diameter of a circle with M1 on the perimeter, angle C1M1C2 is a right angle. Because material point C1 on the wheel has zero velocity the velocity of M1 (and thus the tangent to the curve) is orthogonal to C1M1. Thus the line M1M2 is tangent to the upper cycloid.

The rolling of wheel 2 instantaneously rotating about C2 makes the tangent to the lower cycloid orthogonal to M1M2, the condition for the motion of M2 to be from the wrapping of an inextensible line around the curve A1M1B. This shows that cycloid A2M2B is generated by the wrapping of a line anchored at A1 about the upper cycloid. And this is Huygen’s wrapping mechanism for making a pendulum bob follow a cycloid. Because of this wrapping generation, the arc-length s+s of A1M1B must be 4r and the arc length s of M1B is M1M2 so the length M1C2 is s/2. By the similarity of the two right triangles that share the length s/2 of M1C2:

ys/2=s/22r  y=14rs22=cs22

Filename:tfigure-brach2wheels

which shows that the upper cycloid is the curve y=cs2/2 if c=1/4r, where s is measured from B. This was the equation used to show the constant period nature of the sliding motion of a bead on a frictionless cycloidal curve using power balance.

.

Example: A collar on two rotating rods

Filename:tfigure8-ang-accel-ex
Figure 19.5: A point-mass collar slides simultaneously on 2 rods.

Consider a pair of collars hinged together as a point mass m at P. Each slides frictionlessly on a rod about whose rotation everything is known (see fig. 19.5. What is the force of rod 1 on the mass? For this 2 degree of freedom system let’s use configuration variables θ1 and θ2, and two sets of rotating base vectors: 𝝀ˆ1𝒏ˆ1 and 𝝀ˆ2𝒏ˆ2. These rotating base vectors can be written in terms of the θs, ıˆ and ȷˆ in the standard manner. Assume we know 1 and 2 in this configuration. First find ˙1 and ˙2 by thinking of the velocity of the collar two different ways:

𝒗 = 𝒗 (19.8)
{˙1𝝀ˆ1+θ˙11𝒏ˆ1 = ˙2𝝀ˆ2+θ˙22𝒏ˆ2}
{}𝒏ˆ2 ˙1=θ˙22θ˙11𝒏ˆ1𝒏ˆ2𝝀ˆ1𝒏ˆ2
{}𝒏ˆ1 ˙2=θ˙11θ˙22𝒏ˆ2𝒏ˆ1𝝀ˆ2𝒏ˆ1.

Having found ˙1 and ˙2 we can find the velocity 𝒗 by evaluating either side of eqn. (19.8). Now we apply identical reasoning with the acceleration. The result looks messy, but the approach is straightforward:

𝒂 = 𝒂 (19.9)
{¨1𝝀ˆ1+θ¨11𝒏ˆ1+2˙1θ˙1𝒏ˆ1 = ¨2𝝀ˆ2+θ¨22𝒏ˆ2+2˙2θ˙2𝒏ˆ2}
{}𝒏ˆ2  ¨1=θ¨22+2˙2θ˙2θ¨11𝒏ˆ1𝒏ˆ22˙1θ˙1𝒏ˆ1𝒏ˆ2𝝀ˆ1𝒏ˆ2
{}𝒏ˆ1  ¨2=θ¨11+2˙1θ˙1θ¨22𝒏ˆ2𝒏ˆ12˙2θ˙2𝒏ˆ2𝒏ˆ1𝝀ˆ2𝒏ˆ1

We use the results from eqn. (19.8) for ˙1 and ˙2 to evaluate the right hand sides of the expressions for ¨1 and ¨2 in eqn. (19.9). So now either the left hand side or the right hand side of the second of Eqns. 19.9 can be used to evaluate the acceleration 𝒂, all the terms in both expressions have been found.

To find the forces we use linear momentum balance and the free-body diagram

𝑭tot = m𝒂 (19.10)
{F1𝒏ˆ1+F2𝒏ˆ2 = m{¨1𝝀ˆ1+θ¨11𝒏ˆ1+2˙1θ˙1𝒏ˆ1}}
{}𝝀ˆ2  F2={¨1𝝀ˆ1+θ¨11𝒏ˆ1+2˙1θ˙1𝒏ˆ1}𝝀ˆ2𝒏ˆ1𝝀ˆ2
{}𝝀ˆ1  F1={¨1𝝀ˆ1+θ¨11𝒏ˆ1+2˙1θ˙1𝒏ˆ1}𝝀ˆ1𝒏ˆ2𝝀ˆ1

When actually evaluating the expressions above one can write the base vectors in terms of ıˆ and ȷˆ or use geometry.

Often when working out a problem it is best to not substitute numbers until the end of a problem. This example shows the opposite. If we left the expressions for ˙1 and ˙2 with letters and substituted that into the expressions for ¨1 and ¨2 and left those expressions intact while substituting for the acceleration 𝒂 we would have large expressions for the force components F1 and F2. On the other hand, by using numbers as the calculation progresses the formulas do not grow so much in complexity.

As is the case with most mechanism-mechanics problems, the hard work in getting the dynamics equations is in the kinematics. Generally there are no great short-cuts. There are alternative methods. In this case the location of the base points and the two angles determine the base and two angles of a triangle. This triangle can be solved for the location of the point P. Once that position is known in terms of θ1 and θ2 the velocity and acceleration can be found by differentiation.

As a robot manipulator, this design has the advantage that no motors need to be displaced. It has the disadvantage of requiring good sliding joints.

An alternative solution of the kinematics of this problem would be to use trigonometry to find the position of point P in terms of the angles θ1 and θ2. Then the acceleration of point P is found by taking two time derivatives. The result is approximately equal in the complexity of its appearance to the results used above. That method requires more cleverness at the start (solving an angle-side-angle triangle) and then just brute force differentiation using the chain rule and the product rule.

Filename:tfigure8-ang-accel1
Figure 19.6: The base of a pendulum is vertically vibrated.

Example: Inverted pendulum with a vibrating base

Assume that the base 0 of an inverted point-mass pendulum of length is vibrated according to (see fig. 19.6)

𝒓0=dsinωtıˆ.

The point P thus has acceleration

𝒂P = 𝒂0+𝒂P/0
= dω2sinωtıˆ+θ¨𝒆ˆθθ˙2𝒆ˆR

Now apply linear momentum balance as

𝑭tot = m𝒂P
{mgıˆ+T𝒆ˆR = m{dω2sinωtıˆ+θ¨𝒆ˆθθ˙2𝒆ˆR}}
{}𝒆ˆθ  gıˆ𝒆ˆθ = dω2sinωtıˆ𝒆ˆθ+θ¨
 gsinθ = dω2sinωtsinθ+θ¨

which you write as

0=θ¨+(dω2sinωtg)sinθ/orθ¨=(gdω2sinωt)sinθ/

depending on whether you are analytically or numerically inclined. This is a second order non-linear ordinary differential equation. If ω=0 or d=0 then this is the classic inverted pendulum equation and has solutions that show that the pendulum doesn’t stay near upright. But, you can find by analytic cleverness or numerical integration that for some values of d and ω that the pendulum does not fall down! Just shaking the base keeps the pendulum up (ω2d>g for all cases where this is possible). This isn’t just academic nonsense, the device can be built and the balancing demonstrated.

One alternative to using linear momentum balance in the equations above would be to use angular momentum balance about the point 0’. The resulting vector equation

𝒆ˆR×(mgıˆ)=𝒓P/0×(m𝒂)

yields the same second order scalar ODE.

The vibrating mechanism shown is a “Scotch yoke”. An eccentric disk is mounted to the shaft of a constant angular velocity motor. The rectangular slot moves up and down sinusoidally as the disk wobbles.

SAMPLE 19.1

Filename:sfig8-2-3
Figure 19.7:

A bead on a straight wire. A straight wire is hung between points A and B in the xy plane as shown in the figure. A bead slides down the wire from point A. Write the geometric constraint equation for the bead’s motion and derive the conditions on velocity and acceleration components of the bead due to the constraint.

Solution The constraint on the bead’s motion is that its path must be along the wire, i.e., a straight line between points A and B. Thus the geometric constraint on the motion is expressed by the equation of the path which is

y=hhx.

Since the bead is constrained to move on this path, its velocity and acceleration vectors are also constrained to be directed along AB. This imposes conditions on their x and y components that are easily derived by differentiating the geometric constraint equation with respect to t. Thus,

y˙ = hx˙,
y¨ = hx¨.

Answer: y=hhx,y˙=hx˙,y¨=hx¨

SAMPLE 19.2

Filename:sfig8-2-3a
Figure 19.8:

A particle sliding on a parabolic path. A particle slides on a parabolic trough given by y=ax2 where a is a constant. Write the geometric constraints of motion (on the path, velocity, and acceleration) of the particle. Write the velocity and acceleration of the particle at a generic location (x,y) on its path.

Solution The geometric constraint on the path of the particle is already given, y=ax2. Differentiating the path constraint with respect to time, we get the constraint on velocity and acceleration components.

y˙ = 2axx˙,
y¨ = 2axx¨+2ax˙2.

Now, at a point (x,y), we can write the velocity and acceleration of the particle as

𝒗 = x˙ıˆ+y˙ȷˆ=x˙ıˆ+2axx˙ȷˆ,
𝒂 = x¨ıˆ+y¨ȷˆ=x¨ıˆ+(2axx¨+2ax˙2)ȷˆ.

Answer: v=x˙ıˆ+2axx˙ȷˆ,a=x¨ıˆ+(2axx¨+2ax˙2)ȷˆ

SAMPLE 19.3

Filename:sfig8-2-3b
Figure 19.9:

Circular motion of a particle. A particle is constrained to move on a frictionless circular path of radius R0 with constant angular speed θ˙. There is no gravity. Find the equation of motion of the particle in the x-direction and show that this motion is simple harmonic.

Solution This is a simple problem that you have solved before, probably a few times. Here, we do this problem again just to show how it works out with the constraint machinery in evidence. The geometric constraint on the path of the particle is R=R0 (in polar coordinates). This constraint gives us R˙=0 and R¨=0. Then the acceleration of the particle (in polar coordinates), 𝒂=(R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨)𝒆ˆθ reduces to 𝒂=Rθ˙2𝒆ˆR, (of course).

Filename:sfig8-3-1
Figure 19.10:

The free-body diagram of the particle shows that there is only one force acting on the particle, the normal reaction N of the path acting in the 𝒆ˆR direction. Therefore, the linear momentum balance gives,

N𝒆ˆR=m𝒂=m(Rθ˙2𝒆ˆR)  N=mRθ˙2.

But, to write the equation of motion in the x-direction, we need to write the linear momentum balance in the x-direction. We can write 𝑭=m𝒂 using mixed basis vectors as N𝒆ˆR=m(x¨ıˆ+y¨ȷˆ). Dotting this equation with ıˆ, we get

x¨=Nm(𝒆ˆRıˆ)=mRθ˙2mcosθ=θ˙2(Rcosθ)=θ˙2x

or, x¨+θ˙2x=0, which is the equation of simple harmonic motion in x. You can easily show that the motion in the y-direction is also simple harmonic (y¨+θ˙2y=0).

Answer:

SAMPLE 19.4

A bead slides on a straight wire. Consider the problem of the bead sliding on a straight, inclined, frictionless wire of Sample 19.6 again. Find the position of the bead x(t) and y(t) assuming it slides under gravity starting from rest at A.

Solution To find the position of the bead, we need to write the equation of motion and solve it. This is a single DOF system and, therefore, one scalar equation of motion should suffice.

Filename:sfig8-3-1a
Figure 19.11:

The free-body diagram of the bead is shown in fig. 19.11. Using basis vectors (𝝀ˆ,𝒏ˆ) and (ıˆ,ȷˆ) we write the LMB for the bead as

mgȷˆ+N𝒏ˆ=m𝒂=m(x¨ıˆ+y¨ȷˆ).

We can easily eliminate the constraint force N from this equation by dotting this equation with 𝝀ˆ, which gives

mg(ȷˆ𝝀ˆ) = m[x¨(ıˆ𝝀ˆ)+y¨(ȷˆ𝝀ˆ)]
 gsinθ = x¨cosθy¨sinθ.

But, from the geometric constraint y=h(h/)x, we have y¨=(h/)x¨=(tanθ)x¨. Therefore,

gsinθ=x¨cosθ+x¨tanθsinθ  x¨=gsinθcosθ.

Since gsinθcosθ is constant, we integrate the equation of motion easily to find x(t)=12gsinθcosθt2 since x(0)=0,x˙(0)=0. And, since y=hxtanθ, we have y(t)=h12gsin2θt2.

Answer: x(t)=12ght2sinθcosθ,y(t)=h12gt2sin2θ

SAMPLE 19.5

Filename:sfig8-5-wiper
Figure 19.12:

A bead sliding down a parabolic trough. Consider the problem of Sample 19.7 again. Find the equation of motion of the bead.

Solution This is, again, a one DOF system. Therefore, we will get a single scalar equation of motion. The free-body diagram shown in fig. 19.13 shows two forces acting on the bead. The constraint force N acts normal to the path. Let 𝒆ˆt and 𝒆ˆn be unit vectors tangential and normal to the path, respectively. Then the linear momentum balance gives

Filename:sfig8-5-wiper-a
Figure 19.13:
mgȷˆ+N𝒆ˆn=m𝒂=m(x¨ıˆ+y¨ȷˆ).

To eliminate the unknown constraint force N from this equation, we can take a dot product of this equation with 𝒆ˆt. However, we must first find 𝒆ˆt. Now 𝒆ˆt is the unit tangent vector. So, we can find it by finding a tangent vector to the path (remember gradient of a function f?) and then dividing it by the length of the vector. That is doable but a little complicated. All we need here is the dot product with a vector normal to 𝒆ˆn. Why not use the velocity vector 𝒗=x˙ıˆ+y˙ȷˆ? The velocity vector is always tangential to the path. Furthermore, we know that from the geometric constraint (y=ax2), y˙=2axx˙ and y¨=2axx¨+2ax˙2. Therefore, 𝒗=x˙ıˆ+2axx˙ȷˆ. Now, dotting the LMB equation with 𝒗 we get,

mg(ȷˆ𝒗)2axx˙ = m[x¨(ıˆ𝒗)x˙+y¨(ȷˆ𝒗)2axx˙]
2gaxx˙ = x¨x˙+2axy¨x˙
= x¨+2ax(2axx¨+2ax˙2)y¨
= x¨(1+4a2x2)+4a2xx˙2.

Rearranging the terms above, we get the required equation of motion:

x¨+4a2x1+4a2x2x˙2+2gax1+4a2x2=0.

As you can see, this is a nonlinear ODE. Analytical solution of this equation is rather difficult. We can, however, always solve it numerically. Note that a solution of this equation only gives you x(t), i.e., the x coordinate of the position of the bead. But, you can always find the y coordinate since y=ax2.

Answer: x¨+4a2x1+4a2x2x˙2+2gax1+4a2x2=0

Comment: Note that if we consider x and x˙ to be very small so that we can ignore the x˙2 terms completely, and take 1+4a2x21, then the equation of motion becomes

x¨+(2ga)x=0

which is the equation of simple harmonic motion with frequency 2ga. Thus, if we consider a shallow parabola, and release the bead close to the origin, it executes simple harmonic motion, much like a simple pendulum. This is an intuitively realizable motion.

SAMPLE 19.6

Filename:sfig8-5-2disks
Figure 19.14: A pin is constrained to move in a groove and a slotted arm.

Constrained motion of a pin. During a small interval of its motion, a pin of 100 grams is constrained to move in a groove described by the equation R=R0+kθ where R0=0.3m and k=0.05m. The pin is driven by a slotted arm AB and is free to slide along the arm in the slot. The arm rotates at a constant speed ω=6rad/s. Find the magnitude of the force on the pin at θ=60.

Solution Let 𝑭 denote the net force on the pin. Then from the linear momentum balance

𝑭=m𝒂

where 𝒂 is the acceleration of the pin. Therefore, to find the force at θ=60 we need to find the acceleration at that position.

From the given figure, we assume that the pin is in the groove at θ=60. Since the equation of the groove (and hence the path of the pin) is given in polar coordinates, it seems natural to use polar coordinate formula for the acceleration. For planar motion, the acceleration is

a=(R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨)𝒆ˆθ.

We are given that θ˙ω=6rad/s and the radial position of the pin R=R0+kθ. Therefore, margin: Note that R is a function of θ and θ is a function of time, therefore R is a function of time. Although we are interested in finding R˙ and R¨ at θ=60, we cannot first substitute θ=60 in the expression for R and then take its time derivatives (which will be zero).

θ¨ = dθ˙dt=0( since θ˙= constant)
R˙ = ddt(R0+kθ)=kθ˙ and
R¨ = kθ¨=0.

Substituting these expressions in the acceleration formula and then substituting the numerical values at θ=60, (remember, θ must be in radians!), we get

𝒂 = (R0+kθ)θ˙2Rθ˙2𝒆ˆR+2kθ˙22R˙θ˙𝒆ˆθ
= (0.3m+0.05mπ3)(6rad/s)2𝒆ˆR+20.05m(6rad/s)2𝒆ˆθ
= 13.63m/s2𝒆ˆR+3.60m/s2𝒆ˆθ.

Therefore the net force on the pin is

𝑭 = m𝒂
= 0.1kg(13.63𝒆ˆR+3.60𝒆ˆθ)m/s2
= (1.36𝒆ˆR+0.36𝒆ˆθ)N

and the magnitude of the net force is

F=|𝑭|=(1.36N)2+(0.36N)2=1.41N.

Answer: F=1.41N

SAMPLE 19.7

Filename:sfig8-4-1
Figure 19.15:

A puck sliding on a frictional rotating table. A horizontal turntable rotates with constant angular speed ω=100rpm. A puck of mass m=0.1kg is gently placed on the rotating turntable. The puck begins to slide. The coefficient of friction between the puck and the turntable is 0.25. Find the equation of motion of the puck using

  1. 1.

    A fixed reference frame with cartesian coordinates

  2. 2.

    A rotating reference frame with cartesian coordinates

Solution

  1. 1.

    Equation of motion using a fixed reference frame: The puck has two DOF on the turntable. So, we will need two configuration variables, say x and y, and we will have to find equation of motion for each variable.

    Let us use a fixed cartesian coordinate system with the origin at the center of the turntable. Let 𝒓P=xıˆ+yȷˆ be the position of the puck at some instant t, so that its velocity is 𝒗P=x˙ıˆ+y˙ȷˆ and acceleration is 𝒂P=x¨ıˆ+y¨ȷˆ.

    Filename:sfig8-4-1a
    Figure 19.16: A partial free-body diagram of the puck. For linear momentum balance we need to consider only the forces acting in the plane of motion.

    The free-body diagram of the puck should show three forces — the force of gravity (in 𝒌ˆ direction), the normal reaction of the turntable (in 𝒌ˆ direction) and the friction force 𝑭. Since, there is no motion in the vertical direction, we know that N=mg and that F=|𝑭|=μN=μmg. But, what is the direction of the friction force? Well, we know that it acts in the opposite direction of the relative slip, so that

    𝑭=μN𝒗rel|𝒗rel|.

    So, we need to find 𝒗rel. Now 𝒗rel is the velocity of the puck relative to the turntable, or more precisely, relative to the point on the turntable just underneath the puck. Let us denote that point by P. Clearly, P goes in circles with constant speed, so that its velocity is

    𝒗P=𝝎×𝒓P=θ˙𝒌ˆ×(xıˆ+yȷˆ)=θ˙xȷˆθ˙yıˆ.

    Therefore, the relative velocity, 𝒗rel is

    𝒗rel=𝒗P𝒗P=(x˙+θ˙y)ıˆ+(y˙θ˙x)ȷˆ

    Now the linear momentum balance for the puck in the xy plane gives

    μmg𝒗rel|𝒗rel| = m(x¨ıˆ+y¨ȷˆ)
     x¨ = μg|𝒗rel|(𝒗relıˆ)=μg(x˙+θ˙y)(x˙+θ˙y)2+(y˙θ˙x)2
    and y¨ = μg|𝒗rel|(𝒗relȷˆ)=μg(y˙θ˙x)(x˙+θ˙y)2+(y˙θ˙x)2.

    These are coupled nonlinear ODEs that represent the equations of motion of the puck.

    Answer: x¨=μg(x˙+θ˙y)(x˙+θ˙y)2+(y˙θ˙x)2,y¨=μg(y˙θ˙x)(x˙+θ˙y)2+(y˙θ˙x)2

    Note that these equations are valid only as long as there is relative slip between the puck and the turntable. If the puck stops sliding due to friction, it simply goes in circles with the turntable and, therefore, its equations of motion then are x¨=θ˙2x,y¨=θ˙2y.

  2. 2.

    Equation of motion using a rotating reference frame: Now we derive the equations of motion using a rotating reference frame, , with (x,y) coordinate axes, fixed to the rotating turntable. Let the position of the puck in the rotating frame be 𝒓P/O=xıˆ+yȷˆ. Note that the velocity of the puck in the rotating reference frame is 𝒗P/=x˙ıˆ+y˙ȷˆ, and acceleration is 𝒂P/=x¨ıˆ+y¨ȷˆ.

    Now, from the linear momentum balance (𝑭=m𝒂) for the puck, we get

    Filename:sfig8-6-2
    Figure 19.17: Axes (x,y) represent the rotating frame fixed to the rotating turntable; ie (x,y) rotate with angular velocity 𝝎=θ˙𝒌ˆ.
    μmg𝒗rel|𝒗rel|=m(𝒂P+𝒂P/+2𝝎×𝒗P/)

    where we have used the three term acceleration formula for 𝒂P. Here,

    𝒂P = θ˙2(xıˆ+yȷˆ)
    𝒂P/ = x¨ıˆ+y¨ȷˆ
    2𝝎×𝒗P/ = 2θ˙y˙ıˆ+2θ˙x˙ȷˆ

    Note that the point P, coincident with P and fixed on the turntable, is stationary with respect to the rotating frame. Therefore, the relative velocity of P as observed in the rotating frame is 𝒗rel=𝒗P/=x˙ıˆ+y˙ȷˆ. Substituting these terms in the LMB equation above, we have

    μgx˙ıˆ+y˙ȷˆx˙2+y˙2=(θ˙2x+x¨2θ˙y˙)ıˆ+(θ˙2y+y¨+2θ˙x˙)ȷˆ

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

    x¨ = θ˙2xμgx˙x˙2+y˙2+2θ˙y˙
    y¨ = θ˙2yμgy˙x˙2+y˙22θ˙x˙

    These are the required equations of motion for the puck in the rotating frame.

    Answer: x¨=θ˙2xμgx˙x˙2+y˙2+2θ˙y˙,y¨=θ˙2yμgy˙x˙2+y˙22θ˙x˙


    Once we find a solution x(t) and y(t) of these equations, we can find the solution in the fixed frame by transforming (x,y) to (x,y) through

    {x(t)y(t)}=[cos(θ˙t)sin(θ˙t)sin(θ˙t)cos(θ˙t)]{x(t)y(t)}

    Also, note that when the solution of the equations of motion in the rotating reference frame brings the puck to halt, the puck stops with respect to the rotating turntable. To an observer in the fixed frame, the puck will be going in circles with a constant θ˙.

SAMPLE 19.8

Filename:sfig8-6-2a
Figure 19.18: A collar slides on a frictional bar and finally shoots off the end. The bar rotates with constant angular speed.

A collar sliding on a frictional rod. A collar of mass m=0.5lbm slides on a massless rigid rod OA of length =8ft. The rod rotates counterclockwise with a constant angular speed θ˙=5rad/s. The coefficient of friction between the rod and the collar is μ=0.3. At time t=0s, the bar is horizontal, the collar is at R0=1ft and has radial speed R˙0=μθ˙R0 towards the pivot O. Ignore gravity.

  1. 1.

    How does the position of the collar change with time (i.e., what is the equation of motion of the collar)?

  2. 2.

    Plot the path of the collar starting from t=0s till the collar shoots off the end of the bar.

  3. 3.

    How long does it take for the collar to leave the bar?

Solution

  1. 1.
    Filename:efig1-2-28
    Figure 19.19: Free-Body Diagram of the collar. The only forces on the collar are from the bar: the normal force N and the friction force Fs=μN.

    A Free-Body Diagram of the collar at a general position (R,θ) is shown in Fig. 19.19. The geometry of the position vector and basis vectors is shown in Fig. 19.20. In vector notation, the forces on the collar are 𝑵=N𝒆ˆθ acting normal to the rod and the force of friction 𝑭s=μN𝒆ˆR acting along the rod. Thus linear momentum balance   for the collar is:

    Filename:pfigure-blue-118-2
    Figure 19.20: Geometry of the collar position at an arbitrary time during its slide on the rod.
    𝑭 = m𝒂
    μN𝒆ˆR+N𝒆ˆθ = m[(R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨0)𝒆ˆθ (19.11)

    Note that θ¨=0 because the rod is rotating at a constant rate. Dotting both sides of eqn. (19.11) with 𝒆ˆR and 𝒆ˆθ we get

    [Eqn. (19.11)]𝒆ˆR  μN = m(R¨Rθ˙2)
    or R¨Rθ˙2 = μNm
    [Eqn. (19.11)]𝒆ˆθ  N = 2mR˙θ˙.

    Eliminating N from the last two equations we get

    R¨+2μθ˙R˙θ˙2R=0.

    Since θ˙=ω is constant, the above equation is of the form

    R¨+CR˙ω2R=0 (19.12)

    where C=2μω and ω=θ˙.

    Note that we could derive eqn. (19.12) from eqn. (19.11) in a single step by taking a dot product of eqn. (19.11) with 𝒆ˆR+μ𝒆ˆθ. Why? Look at the left hand side of eqn. (19.11). The net reaction force is N(μ𝒆ˆR+𝒆ˆθ) acting in the direction (μ𝒆ˆR+𝒆ˆθ). If we dot the reaction force with a vector normal to its direction, we get rid of the reaction force. A vector normal to (μ𝒆ˆR+𝒆ˆθ) is 𝒆ˆR+μ𝒆ˆθ.

    [Eqn. (19.11)](𝒆ˆR+μ𝒆ˆθ)  0 = m[(R¨Rθ˙2)𝒆ˆR+2R˙θ˙𝒆ˆθ](𝒆ˆR+μ𝒆ˆθ)
     0 = R¨+2μθ˙R˙θ˙2R

    which is the same equation as eqn. (19.12).

    Solution of equation (19.12): The characteristic equationmargin: If we assume a solution of the form R(t)=eλt with some unknown λ and plug back into the given differential equation, we get an algebraic equation in λ which is called the characteristic equation. For more information, see box 11.1 on page 11.1, eqn. (11.9). associated with Eqn. (19.12) is

    λ2+Cλω2 = 0
     λ = C±C2+4ω22
    = ω(μ±μ2+1).

    Therefore, the solution of Eqn. (19.12) is

    R(t) = Aeλ1t+Beλ2t
    = Ae(μ+μ2+1)ωt+Be(μμ2+1)ωt.

    Substituting the given initial conditions: R(0)=R0=1ft and R˙(0)=μθ˙R0, we getmargin: From the general solution and the given initial conditions, we have: R0 = A+B μωR0 = Aλ1+Bλ2 Solving these two equations simultaneously, we get A=λ2+μωλ2λ1R0, and B=λ1+μωλ2λ1R0. Substituting the values of λ1 and λ2, we get A=B=R02=1ft2.

    R(t)=1ft2[e(μ+μ2+1)ωt+e(μμ2+1)ωt]. (19.13)

    Answer: R(t)=1ft2[e(μ+μ2+1)ωt+e(μμ2+1)ωt].

  2. 2.

    To draw the path of the collar we need both R and θ. Because θ˙=5rad/s= constant,

    θ=θ˙t=(5rad/s)t.
    Filename:pfigure-s95f3a
    Figure 19.21: Finding the final time tf from the graph of R(t) such that R(tf)=Rf=8ft.

    Now we can take various values of t from 0s to, say, 1s, and calculate values of θ and R. Plotting all these values of R and θ, however, does not give us the correct path after the collar reaches the end of the bar, at R=length of the bar=8ft. So we need to find the time tf when R(tf)=8ft. Equation (19.13) is a nonlinear algebraic equation for t which we can solve iteratively on a computer, or with some patience, even on a calculator. One way to find tf would be to simply plot R(t) and find the intersection with R=8 (see fig. 19.21) and read the corresponding value of t. Following either of these methods we find that tf0.74518s here. Now, we can plot the path of the collar by computing R and θ from t=0 to t=tf and making a polar plot on a computer as follows (pseudocode).

    tf = 0.74518                % final value of t
    t = 0:tf/100:tf;            % take 101 points in [0 tf]
    R0 = 1; w = 5; mu = .3;     % initialize variables
    f1 = -mu + sqrt(mu^2 +1);   % first partial exponent
    f2 = -mu - sqrt(mu^2 +1);   % second partial exponent
    R = 0.5*R0*(exp(f1*w*t) + exp(f2*w*t));  % calculate R
    theta = w*t;                % calculate theta
    polarplot(theta,r)
    
    Filename:Danef94s1q2
    Figure 19.22: Plot of the path of the collar until it leaves the rod.

    The plot produced thus is shown in Fig. 19.22.


  3. 3.

    The time tf computed above was

    tf=0.7452s.

    By plugging this value in the expression for R(t) (Eqn. (19.13) we get, indeed,

    R=8ft.

    Answer: tf=0.7452s

Note, although we have not checked it explicitly, the increase in kinetic energy of the particle comes from the work of the force from the rod on the particle:

ΔEK=𝑭𝒗𝑑t

SAMPLE 19.9

Filename:pfigure-blue-123-1
Figure 19.23: A collar of mass m slides on a massless and frictionless bar, bent at an angle ϕ with the vertical axis. The bar is rotating at a constant angular speed ω about the z-axis.

A collar sliding on a rotating rod in 3-D. A massless and frictionless rod AB is rotating about the vertical axis through point A. The rod is bent at an angle ϕ from the vertical and is rotating with a constant angular speed ω about the vertical axis. A small collar of mass m slides on the rod. Assume that at time t=0 the collar is released from a rest position with respect to the rod at a distance R0 from the axis of rotation. There is no gravity.

  1. 1.

    Find the equation of motion for the collar (a differential equation for the position of the collar).

  2. 2.

    How does the distance of the collar from the vertical axis change with time?

  3. 3.

    For ϕ=π/2, show that the solution obtained in (ii) above is the same as that obtained in Sample 19.15 for μ=0.


Solution Since the bar is bent and it rotates about the vertical axis, it sweeps a conical surface about the axis of rotation. As the collar slides on the rod, it traces a path on this surface. Let (R,θ,z) be the cylindrical coordinates of the collar at any general time t (Fig. 19.24(a)). Since the rod is frictionless and there is no gravity, the only force acting on the collar is the normal reaction from the rod. This force is shown in the free-body diagram of the collar in Fig. 19.24(b). Note that there are a lot of possible directions for a normal to the rod at the collar. In fact, any vector in the plane perpendicular to the rod at the collar is normal to the rod. So, at this point let us write

𝑵=N𝒏ˆ

where 𝒏ˆ is a unit normal to the rod. Thus, if 𝝀ˆ is a unit vector along the rod, then

𝝀ˆ𝒏ˆ=0 (19.14)

where

𝝀ˆ=𝒓AB|𝒓AB|=sinϕ𝒆ˆR+cosϕ𝒌ˆ.
Filename:pfigure-blue-119-2
Figure 19.24: (a) Cylindrical coordinates R and z of the collar (θ is not shown) and the orientation of the cylindrical basis vectors. (b) Free-body diagram of the collar (there is no gravity). (c) Instantaneous R-theta plane of the collar.
  1. 1.

    Equation of Motion: We now write the linear momentum balance   for the collar:

    𝑭=m𝒂

    where

    𝑭 = N𝒏ˆ
    𝒂 = (R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨0)𝒆ˆθ+z¨𝒌ˆ

    Dotting both sides of the equation 𝑭=m𝒂 by 𝝀ˆ we get

    N𝒏ˆ𝝀ˆ0 = m[(R¨Rθ˙2)𝒆ˆR𝝀ˆsinϕ+2R˙θ˙𝒆ˆθ𝝀ˆ0+z¨𝒌ˆ𝝀ˆcosϕ]
    (R¨Rθ˙2)sinϕ+z¨cosϕ=0.

    This expression is an equation of motion of the collar. However, it is only a single equation in terms of derivatives of two variables R and z. Now, from geometry

    R=ztanϕ  z¨=R¨/tanϕ.

    Substituting this relationship in the equation of motion we get

    R¨(sinϕ+cosϕtanϕ)Rθ˙2=0.

    Noting that θ˙=ω= a constant, we may write the above equation, with some trigonometric simplifications, as

    R¨+(ωsinϕ)2R=0 (19.15)

    which is an equation of motion of the collar in terms of its distance R from the vertical axis.

  2. 2.

    Solution of the equation of motion: The solution of Eqn. (19.15) is given by margin: You may find the solution by solving the corresponding characteristic equation (see the solution of Eqn. (19.15) in Sample 19.15, for example), or by looking up the table of “The Simplest ODEs and Their Solutions” on page 22 of the text, or by guessing a solution yourself if you have some experience with ODEs.

    R(t)=C1e(ωsinϕ)t+C2e(ωsinϕ)t

    where C1 and C2 are arbitrary constants to be determined from the initial conditions. Substituting the given initial conditions: R(0)=R0 and R˙(0)=0 we get

    C1=C2=R02.

    Therefore, the solution may be written as

    R(t)=R02[e(ωsinϕ)t+e(ωsinϕ)t]. (19.16)
  3. 3.

    Special case, ϕ=π/2: Substituting ϕ=π/2 in Eqn. (19.16) we get

    R(t)=R02[eωt+eωt]

    which is the same solution as obtained in Eqn. (19.16) in Sample 19.15 for μ=0 and R0=1ft.

SAMPLE 19.10

Filename:pfigure-s95q14
Figure 19.25: A motor driven bead slides down a helical wire-frame described by the equations R=R0θ,z=2R0θ.

Osculating circle in 3-D. A small bead is driven down a a wire-frame bent in the shape of a conical helix by a tiny motor imbedded in the bead. The combined mass of the bead and the motor is m=200 gm. The shape of the helix is given: R=R0θ,z=2R0θ and where R0=0.3m. At the instant when θ=2radians, the angular speed and the angular acceleration of the bead are θ˙=1rad/s and θ¨=2rad/s2. Find

  1. 1.

    the net normal force on the bead and

  2. 2.

    the radius of the osculating circle at the instant given.


margin:

Solution

  1. 1.

    We can find the net normal force on the bead from the linear momentum balance of the bead (see the free-body diagram of the bead):

    Filename:Danef94s3q2
    Figure 19.26: Free-body diagram of the bead. Note that the net normal force is Fn=𝑭𝒆ˆn where 𝑭=N𝒆ˆn+T𝒆ˆtmg𝒌ˆ.
    𝑭 = m𝒂
    N𝒆ˆn+T𝒆ˆtmg𝒌ˆ = m(at𝒆ˆt+an𝒆ˆn)
     Fn 𝑭𝒆ˆn=man.

    Thus we need to find the normal acceleration of the bead. We can write the position of the bead as

    𝒓 = RR0θ𝒆ˆR+z2R0θ𝒌ˆ
     𝒗 R˙𝒆ˆR+Rθ˙𝒆ˆθ+z˙𝒌ˆ
    = R0θ˙𝒆ˆR+R0θθ˙𝒆ˆθ+2R0θ˙𝒌ˆ
    = R0θ˙(𝒆ˆR+θ𝒆ˆθ+2𝒌ˆ),
    and 𝒂 (R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨)𝒆ˆθ+z¨𝒌ˆ
    = (R0θ¨R0θθ˙2)𝒆ˆR+(2R0θ˙2+R0θθ¨)𝒆ˆθ+2R0θ¨𝒌ˆ
    = R0[(θ¨θθ˙2)𝒆ˆR+(θ˙2+θθ¨)𝒆ˆθ+2θ¨𝒌ˆ]

    Substituting the given numerical values for θ,θ˙,R0 and θ¨ in the above expressions for 𝒗 and 𝒂 we get the velocity and the acceleration of the bead at the instant of interest:

    𝒗 = 0.3m/s(𝒆ˆR+2𝒆ˆθ+2𝒌ˆ)
    𝒂 = 1.2m/s2(2𝒆ˆθ+𝒌ˆ).

    In path coordinates,

    𝒂=𝒂t+𝒂n=at𝒆ˆt+an𝒆ˆn

    where

    𝒆ˆt = 𝒗|𝒗| (19.17)
    = 0.3m/s(𝒆ˆR+2𝒆ˆθ+2𝒌ˆ)0.31+4+4m/s (19.18)
    = 13(𝒆ˆR+2𝒆ˆθ+2𝒌ˆ). (19.19)

    Therefore,

    𝒂t = (𝒂𝒆ˆt)𝒆ˆt
    = (4.8+2.43m/s2)𝒆ˆt
    = 0.8m/s2(𝒆ˆR+2𝒆ˆθ+2𝒌ˆ)
    𝒂n = 𝒂𝒂t
    = [2.4𝒆ˆθ+1.2𝒌ˆ0.8m/s2(𝒆ˆR+2𝒆ˆθ+2𝒌ˆ)]
    = (0.8𝒆ˆR+0.8𝒆ˆθ0.4𝒌ˆ)m/s2
    an = |𝒂n|=1.2m/s2
    𝒆ˆn = 𝒂nan
    = 23𝒆ˆR+23𝒆ˆθ13𝒌ˆ.

    Hence, the net normal force on the bead

    𝑭n = man𝒆ˆn
    = 0.2kg1.2m/s213(2𝒆ˆR+2𝒆ˆθ𝒌ˆ)
    = (0.16𝒆ˆR+0.16𝒆ˆθ0.08𝒌ˆ)N.

    Answer: Fn=0.08N(2eˆR+2eˆθ+kˆ)

  2. 2.

    For calculating the radius ρ of the osculating circle, we note that

    an = v2ρ
     ρ = v2an=(0.9m/s)21.2m/s2=0.675m.

    Answer: ρ=0.675m

Problems for 19.1 Mechanics of a constrained particle

19.1.1

A bead slides on a frictionless circular hoop . The mass of the bead mbead=2grams,

mass of hoop mhoop=1kg and

radius of hoop Rhoop=3m.

Neglect gravity. The center of the circular hoop is the origin O of a fixed (Newtonian) coordinate system Oxyz. The hoop is on the xy plane. The hoop is kept from moving by little angels who let the bead slide by unimpeded. At t=0, the speed of the bead is 4m/s, it is traveling counter-clockwise (looking down the z-axis), and it is on the +x-axis. There are no other external forces applied.

  1. (a)

    At t=0 what is the bead’s kinetic energy ?

  2. (b)

    At t=0 what is the bead’s linear momentum ?

  3. (c)

    At t=0 what is the bead’s angular momentum about the origin?

  4. (d)

    At t=0 what is the bead’s acceleration?

  5. (e)

    At t=0 what is the radius of the osculating circle of the bead’s path ρ?

  6. (f)

    At t=0 what is the force of the hoop on the bead?

  7. (g)

    At t=0 what is the net force of the angel’s hands on the hoop?

  8. (h)

    At t=27s what is the x-component of the bead’s linear momentum?

19.1.2  A warehouse operator wants to move a crate of weight W=100lb from a 6ft high platform to ground level by means of a roller conveyer as shown. The rollers on the inclined plane are well lubricated and thus assumed frictionless; the rollers on the horizontal conveyer are frictional, thus providing an effective friction coefficient μ. Assume the rollers are massless.

  1. (a)

    What are the kinetic and potential energies(pick a suitable datum) of the crate at the elevated platform, point A?

  2. (b)

    What are the kinetic and potential energies of the crate at the end of the inclined plane just before it moves onto the horizontal conveyor, point B?

  3. (c)

    Using conservation of energy, calculate the maximum speed attained by the crate at point B, assuming it starts moving from rest at point A.

  4. (d)

    If the crate slides a distance x, say, on the horizontal conveyor, what is the energy lost to sliding in terms of μ, W, and x?

  5. (e)

    Calculate the value of the coefficient of friction, μ, such that the crate comes to rest at point C.

Filename:bikefork1-ang-accel
Figure 19.27

19.1.3  On a wintry evening, a student of mass m starts down the steepest street in town, Steep Street, which is of height h and slope θ. At the top of the slope she starts to slide (coefficient of dynamic friction μ).

  1. (a)

    What is her initial kinetic and potential energy?

  2. (b)

    What is the energy lost to friction as she slides down the hill?

  3. (c)

    What is her velocity on reaching the bottom of the hill? Ignore air resistance and all cross streets (i.e. assume the hill is of constant slope).

  4. (d)

    If, upon reaching the bottom of the street, she collides with another student of mass M and they embrace, what is their instantaneous mutual velocity just after the embrace?

  5. (e)

    Assuming that friction still acts on the level flats on the bottom, how much time will it take before they come to rest?

Filename:bikefork-ang-accel
Figure 19.28

19.1.4   Reconsider the system of blocks in problem 2.3, this time with equal mass, m1=m2=m. Also, now, both blocks are frictional and sitting on a frictional surface. Assume that both blocks are sliding to the right with the top block moving faster. The coefficient of friction between the two blocks is μ1. The coefficient of friction between the lower block and the floor is μ2. It is known that μ1>2μ2.

  1. (a)

    Draw free-body diagrams of the blocks together and separately.

  2. (b)

    Write the equations of linear momentum balance for each block.

  3. (c)

    Find the acceleration of each block for the case of frictionless blocks. For the frictional blocks, find the acceleration of each block. What happens if μ1>>μ2?

19.1.5  An initially motionless roller-coaster car of mass 20 kg is given a horizontal impulse 𝑭𝑑t=Pıˆ at position A, causing it to move along the track, as shown below.

  1. (a)

    Assuming that the track is perfectly frictionless from point A to C and that the car never leaves the track, determine the magnitude of the impulse I so that the car “just makes it” over the hill at B.

  2. (b)

    In the ensuing motion, assume that the horizontal track C-D is frictional and determine the coefficient of friction required to bring the car to a stop at D.

Filename:bikefork1-alt
Figure 19.29

19.1.6  Masses m1=1kg and m2kg move on the frictionless varied terrain shown. Initially, m1 has speed v1=12m/s and m2 is at rest. The two masses collide on the flat section. The coefficient of restitution in the collision is e=0.5. Find the speed of m2 at the top of the first hill of elevation h1=1m. Does m2 make it over the second hill of elevation h2=4m?

Filename:bikefork-alt
Figure 19.30

19.1.7   A Scotch yoke is a device for converting rotary motion into linear motion. In this case it is used to make a horizontal platform go up and down. The motion of the platform is y=Asin(ωt) with constant ω. A dead bug with mass m is standing on the platform. Her feet have no glue on them. There is gravity.

  1. (a)

    Draw a free-body diagram of the bug.

  2. (b)

    What is the acceleration of the platform and, hence, the bug?

  3. (c)

    Write the equation of linear momentum balance for the bug.

  4. (d)

    What condition must be met so that the bug does not bounce the platform?

  5. (e)

    What is the maximum spin rate ω that can be used if the bug is not to bounce off the platform?

Filename:pfigure-blue-128-1
Figure 19.31

19.1.8  A bug of mass m walks straight-forwards with speed vA and rate of change of speed v˙A on a straight light (assumed to be massless) stick. The stick is hinged at the origin so that it is always horizontal but is free to rotate about the z axis. Assume that the distance the bug is from the origin, , the angle the stick makes with the x axis, ϕ, and its rate of change, ϕ˙, are known at the instant of interest. Ignore gravity. Answer the following questions in terms of ϕ, ϕ˙, , m, vA, and v˙A.

  1. (a)

    What is ϕ¨?

  2. (b)

    What is the force exerted by the rod on the support?

  3. (c)

    What is the acceleration of the bug?

Filename:tfigure8-syst-bods
Figure 19.32

19.1.9  A new kind of gun. Assume ω=ω0 is a constant for the rod in the figure. Assume the mass is free to slide. At t=0, the rod is aligned with the x-axis and the bead is one foot from the origin and has no radial velocity (dR/dt=0).

  1. (a)

    Find a differential equation for R(t). Answer: R¨(t)ωo2R(t)=0.

  2. (b)

    Turn this equation into a differential equation for R(θ). Answer: d2R(θ)dθ2R(θ)=0.

  3. (c)

    How far will the bead have moved after one revolution of the rod? How far after two? Answer: R(θ=2π)=267.7ft, R(θ=4π)=1.43×105ft.

  4. (d)

    What is the speed of the bead after one revolution of the rod (use ω0=2πrad/s=1rev/sec)? Answer: v=1682.3ft/s

  5. (e)

    How much kinetic energy does the bead have after one revolution and where did it come from? Answer: EK=(2.83×106)×(mass)(ft/s)2.

Filename:sfig8-7-2
Figure 19.33

19.1.10  The new gun gets old and rusty.

Reconsider the bead on a rod in problem 19.1. This time, friction cannot be neglected. The friction coefficient is μ. At the instant of interest the bead with mass m has radius R0 with rate of change R˙0. The angle θ is zero and ω is a constant. Neglect gravity.

  1. (a)

    What is R¨ at this instant? Give your answer in terms of any or all of R0, R˙0, ω, m, μ, 𝒆ˆθ, 𝒆ˆR, ıˆ, and ȷˆ. Answer: R¨=R0ω22μR˙0ω.

  2. (b)

    After a very long time it is observed that the angle ϕ between the path of the bead and the rod/trough is nearly constant. What, in terms of μ, is this ϕAnswer: ϕ=sin1[1(1+μ2μ)2+1].

19.1.11  A newer kind of gun. As an attempt to make an improvement on the ‘new gun’ demonstrated in problem 19.1, a person adds a length to the shaft on which the bead slides. Assume there is no friction between the bead (mass m) and the wire. Assume the bead starts at s=0 with s˙=v0. The rigid rod, on which the bead slides, rotates at a constant rate ω=ω0. Find s(t) in terms of ,m,v0, and ω0.

Filename:sfig8-7-2a
Figure 19.34

19.1.12   Slippery bead on straight rotating stick. A long stick with mass ms rotates at a constant rate and a bead, modeled as a point mass, slides on the stick. The stick rotates at constant rate 𝝎=ω𝒌ˆ. Neglect gravity. The bead has mass mb. Initially the mass is at a distance R0 from the hinge point on the stick and has no radial velocity (R˙=0). The initial angle of the stick is θ=0, measured counterclockwise from the positive x axis.

  1. (a)

    What is the torque, as a function of the net angle the stick has rotated, θ, required in order to keep the stick rotating at a constant rate ?

  2. (b)

    What is the path of the bead in the xy-plane? (Draw an accurate picture showing about one half of one revolution.)

  3. (c)

    How long should the stick be if the bead is to fly in the negative x-direction when it gets to the end of the stick?

  4. (d)

    Add friction. How does the speed of the bead over one revolution depend on μ, the coefficient of friction between the bead and the wire? Make a plot of |𝒗one revolution| versus μ. [Note, you have been working with μ=0 in the problems above.]

19.1.13  Mass on a lightly greased slotted turntable or spinning uniform rod. Assume that the rod/turntable in the figure is massless and also free to rotate. Assume that at t=0, the angular velocity of the rod/turntable is 1rad/s, that the radius of the bead is one meter, and that the radial velocity of the bead, dR/dt, is zero. The bead is free to slide on the rod. Where is the bead at t=5sec?

Filename:sfig8-7-2again
Figure 19.35

19.1.14   A bead slides in a frictionless slot in a turntable. The turntable spins a constant rate ω. The slot is straight and goes through the center of the turntable. If the bead is at radius Ro with R˙=0 at t=0, what are the components of the acceleration vector in the directions normal and tangential to the path of the bead after one revolution? Neglect gravity.

19.1.15  A bead of mass m=1gm is constrained to slide in a straight frictionless slot in a disk which is spinning counterclockwise at constant rate ω=3rad/s. At time t=0, the slot is parallel to the x-axis and the bead is in the center of the disk moving out (in the plus x direction) at a rate R˙o=.5m/s. After a net rotation θ of one and one eighth (1.125) revolutions, what is the force 𝑭 of the disk on the bead? Express this answer in terms of ıˆ and ȷˆ. Make the unreasonable assumption that the slot is long enough to contain the bead for this motion.

Filename:sfig8-7-2disks
Figure 19.36

19.1.16  Bead on springy leash in a slot on a turntable. The bead in the figure is held by a spring that is relaxed when the bead is at the origin. The constant of the spring is k. The turntable speed is controlled by a strong stiff motor.

  1. (a)

    Assume ω=0 for all time. What are possible motions of the bead?

  2. (b)

    Assume ω=ω0 is a constant. What are possible motions of the bead? Notice there are two cases depending on the value of ω0. What is going on here?

Filename:sfig8-4-4
Figure 19.37

Mass in a slot.

19.1.17  A small bead with mass m slides without friction on a rigid rod which rotates about the z axis with constant ω (maintained by a stiff motor not shown in the figure). The bead is also attached to a spring with constant K, the other end of which is attached to the rod. The spring is relaxed when the bead is at the center position. Assume the bead is pulled to a distance d from the center of the rod and then released with an initial R˙=0. If needed, you may assume that K>mω2.

  1. (a)

    Derive the equation of motion for the position of the bead R(t).

  2. (b)

    How is the motion affected by large versus small values of K?

  3. (c)

    What is the magnitude of the force of the rod on the bead as the bead passes through the center position? Neglect gravity. Answer in terms of m, ω, d, K.

  4. (d)

    Write an expression for the Coriolis acceleration. Give an example of a situation in which this acceleration is important and explain why it arises.

Filename:sfig8-4-4a
Figure 19.38

19.1.18  A small ball of mass mb=500grams may slide in a slender tube of length l=1.2m. and of mass mt=1.5kg. The tube rotates freely about a vertical axis passing through its center C. (Hint: Treat the tube as a slender rod.)

  1. (a)

    If the angular velocity of the tube is ω=8rad/s as the ball passes through C with speed relative to the tube, v0, calculate the angular velocity of the tube just before the ball leaves the tube. Answer: ω=4rad/s, where the minus sign ‘-’ means ‘just before leaving’.

  2. (b)

    Calculate the angular velocity of the tube just after the ball leaves the tube. Answer: ω+=4rad/s, where the plus sign ‘+’ means ‘just after leaving’.

  3. (c)

    If the speed of the ball as it passes Cis v=1.8m/s, determine the transverse and radial components of velocity of the ball as it leaves the tube. Answer: vr=3.841m/s, vt=2.4m/s.

  4. (d)

    After the ball leaves the tube, what constant torque must be applied to the tube (about its axis of rotation) to bring it to rest in 10s? Answer: torque=0.072Nm

Filename:sfig8-4-4b
Figure 19.39

19.1.19  Toy train car on a turn-around. The 0.1kg toy train’s speedometer reads a constant 1m/s when, heading west, it passes due north of point O. The train is on a level straight track which is mounted on a spinning turntable whose center is ‘pinned’ to the ground. The turntable spins at the constant rate of 2rad/s. What is the force of the turntable on the train? (Don’t worry about the z-component of the force.) Answer: F=0.6Nȷˆ.

Filename:sfig8-4-4c
Figure 19.40

19.1.20  Due to forces not shown, the cart moves to the right with constant acceleration ax. The ball B has mass mB. At time t=0, the string AB is cut. Find

  1. (a)

    the tension in string BC before cutting, Answer: TBC=m3413(2ax+g).

  2. (b)

    the absolute acceleration of the mass at the instant of cutting, Answer: a=134[(25ax+15g)ıˆ+(15ax25g)ȷˆ]

  3. (c)

    the tension in string BC at the instant of cutting. Answer: TBC=m134(5ax+3g).

Filename:pfigure-s94h13p2
Figure 19.41

19.1.21  The forked arm mechanism pushes the bead of mass 1kg along a frictionless hyperbolic spiral track given by r=0.5θ(m/rad). The arm rotates about its pivot point at O with constant angular acceleration θ¨=1rad/s2 driven by a motor (not shown). The arm starts from rest at θ=0.

  1. (a)

    Determine the radial and transverse components of the acceleration of the bead after 2s have elapsed from the start of its motion.

  2. (b)

    Determine the magnitude of the net force on the mass at the same instant in time.

Filename:pfigure-f93f5
Figure 19.42

19.1.22  A rod is on the palm of your hand at point A. Its length is . Its mass m is assumed to be concentrated at its end at C. Assume that you know θ and θ˙ at the instant of interest. Also assume that your hand is accelerating both vertically and horizontally with 𝒂hand=ahxıˆ+ahyȷˆ. Coordinates and directions are as marked in the figure.

  1. (a)

    Draw a Free-Body Diagram of the rod.

  2. (b)

    Assume that the hand is stationary. Solve for θ¨ in terms of g, , m, θ, and θ˙.

  3. (c)

    If the hand is not stationary but θ¨ has been determined somehow, find the vertical force of the hand on the rod in terms of θ¨, θ˙, θ, g, , m, ahx, and ahy.

Filename:pfigure-s94h13p3
Figure 19.43

Rod on a Moving Support(2-D).

19.1.23  Balancing a broom. Assume the hand is accelerating to the right with acceleration 𝒂=aıˆ. What is the force of the hand on the broom in terms of m,,θ,θ˙,a,ıˆ,ȷˆ, and g? (You may not have any 𝐞ˆR or 𝐞ˆθ in your answer.)

Filename:p-s96-p3-3
Figure 19.44

19.2 Mechanics of one-degree-of-freedom 2-D mechanisms

A one-degree-of-freedom mechanism is a collection of parts linked together so that they can move in only one waymargin: No physical system can move in only one way. The idea that a machine has one degree of freedom is an idealization that takes the rigid-object description of the parts and the ideal nature of the connections literally. In fact, all parts deform at least a little, and no connections are so precise as to be exact geometric constraints. It may be reasonable to respect the standard rigid-object idealizations and consider a machine as one-degree-of-freedom for basic analysis, as we do in this section. When, for example, trying to figure out why a machine vibrates in an undesirable way it may also be reasonable to relax some of the assumptions we make here and consider more degrees of freedom. .

The word “freedom” has to be taken lightly here because in practice even the one “freedom” is often controlled or restricted. Frankly, most machine designers don’t trust the laws of mechanics to enforce motions that they want. Instead they choose kinematic restrictions that enforce the desired motions and then use a motor, a computer controlled actuator or big flywheel to keep that motion moving at a prescribed rate.

We consider here machines that can move in just one way, whether or not that one motion is free. So in this sense, “one”-degree-of-freedom machines include machines with no freedom at all, just so long as they move in only one way.

Some familiar examples of one-degree-of-freedom mechanisms are a 1-D spring and mass, a pendulum, a slider-crank, a grounded 4-bar linkage, and a gear train.

Most ideal constraints are workless constraints

A fruitful equation for studying one-degree-of-freedom mechanisms is power balance, or for conservative systems, energy balance. The reason these equations are so useful is because most ideal connections are workless. That is:

the net work of the interaction forces (and moments) of a pair of parts that are connected with the standard ideal connections is zero. This includes welds, frictionless hinges, frictionless sliding contact, rolling contact, or parts connected by a massless inextensible link.

Example: A frictionless hinge is a workless constraint

Filename:bikefork1-ang-mom
Figure 19.45: A frictionless hinge is a workless constraint. The net work of the interaction force on the two contacting bodies is zero.

Body 𝒜 is connected to body by a frictionless hinge at C (see fig. 19.45). The force on body at C from 𝒜 is 𝑭C and the force on 𝒜 from body at C is 𝑭C. The power of the interaction force on body is P𝒜on=𝑭C𝒗C. This power contributes to the increase in the kinetic energy of . The power of the interaction force on 𝒜 is Pon𝒜=𝑭C𝒗C=P𝒜on. So the contribution to the increase in the kinetic energy of body 𝒜 is minus the contribution to body and the net power on the system of two bodies is zero. Writing this out,

PtotalThe net power of the pair of interaction forces on the pair of bodies = Pon𝒜+P𝒜on
= 𝑭C𝒗C+(𝑭C)𝒗C
= (𝑭C𝑭C)𝒗C
= 0.

Basically the same situation holds for all the standard ideal connections as explained in the box on page 19.2.

If one of two interacting bodies is known to be stationary, like the ground, then the work of the constraint forces is zero on both of the bodies. Thus the work of the hinge force on a pendulum, and the ground reaction forces on a frictionlessly sliding body or the ground force on a perfectly rolling body is zero. But be careful with the words “workless constraint forces”, however.

The workless constraint connecting moving bodies 𝒜 and is likely to do positive work on one of the bodies and negative work on the other.

It is just the net work on the two bodies which is zero.

Box 19.2 Ideal constraints and workless constraints

All of the ideal constraints we consider are interactions between two bodies 𝒜 and . One of these could be the ground. Let’s take the interaction force 𝑭 and moment 𝑴 to be the force and moment of 𝒜 on . The point of interaction is A on 𝒜 and B on . By the principle of action and reaction, the net power of the interaction force on the two bodies is

P = 𝑭𝒗B+𝑴𝝎+(𝑭)𝒗A+(𝑴)𝝎𝒜
= 𝑭𝒗B/A+𝑴𝝎/𝒜.

All of our ideal constraints are designed to exactly make these dot products zero. The ideal hinge is considered in the text. Another example is perfect rolling. In that case the interaction moment is assumed to be zero. The no-slip condition means that 𝒗B/A=𝟎. On the other hand for frictionless sliding there 𝒗B/A can have a component tangent to the surfaces. But that is exactly the direction where the friction force is assumed to be zero.

And so it is for all of the ideal “workless” constraints.

Examples of non-workless constraints, that is, interactions that contribute to the energy equations are: sliding with non-zero friction, joints with non-zero friction torques, joints with motors, or interactions mediated by springs, dampers or actuators.

Energy method: single degree of freedom systems

Although linear and angular momentum balance apply to a single degree of freedom system and all of its parts, often one finds what one wants with a single scalar equation, namely energy or power balance.

Imagine a complex machine that only has one degree of freedom, meaning the position of the whole machine is determined by a single configuration variable, call it q. Further assume that the machine has no motion when q˙=0. The variable q could be, for example, the angle of one of the linked-together machine parts. Also, assume that the machine has no dissipative parts: no friction, no collisions, no inelastic deformation. Because q characterizes the position of all of the parts of the system we can, in principle, calculate the potential energy of the system as a function of q,

EP=EP(q).

We find this function by adding up the potential energies of all the springs in the machine and the gravitational potential energies of the parts.

Similarly we can write the system’s kinetic energy in terms of q and its rate of change q˙. Because at any configuration the velocity of every point in the system is proportional to q˙ we can write the kinetic energy as:

EK=M(q)q˙2/2

where M(q) is a function that one can determine by calculating the machine’s total kinetic energy in terms of q and q˙ and then factoring q˙2 out of the resulting expression.

Now, if we accept the equation of mechanical energy conservation we have

constant =ET conservation of energy,
0 =ddtET taking one time derivative,
=ddt[EP+EK] total energy is potential plus kinetic
=ddt[EP(q)+12M(q)q˙2] substituting from paragraphs above
 0 =ddq[EP(q)]q˙+12ddq[M(q)]q˙q˙2+M(q)q˙q¨
0 =ddq[EP(q)]+(12ddq[M(q)])q˙2+M(q)q¨ cancelling q˙
0 =f1(q)+f2(q)q˙2+f3(q)q¨ (19.20)
withf1(q)ddq[EP(q)],
f2(q)12ddq[M(q)], and
f3(q)M(q).

The cancellation of q˙ above lacks mathematical rigor, but doesn’t cause problemsmargin: The cancellation of the factor q˙ from equation 19.20 depends on q˙ being other than zero. While moving, q˙ is not zero. Strictly we cannot cancel the q˙ term from the equation at the instants when q˙=0. However, to say that a differential equation is true except for certain instants in time is, in practice to say that it is always true, at least if we make reasonable assumptions about the smoothness of the motions. .The equation of motion is complicated because when we take the time derivative of a function of M(q) and EP(q) we have to use the chain rule. Also, because we have products of terms, we had to use the product rule.

Eqn. 19.20 is the general equation of motion of a conservative one-degree-of-freedom system. It is really just a special case of the equation of motion for one-degree-of-freedom systems found from power balance. Rather than memorizing eqn. (19.20) it is probably best to look at its derivation as an algorithm to be reproduced on a problem by problem basis.

Example: Spring and mass

Filename:bikefork-ang-mom
Figure 19.46: The familiar one degree of freedom spring and mass system.

Although the motion of a spring and mass system can be found easily enough from linear momentum balance, it is also a good example for energy balance (see fig. 19.46). Using conservation of energy for the spring and mass system:

ET = constant
0 = ddtET
= EK˙+EP˙
= ddt(mv2/2)+ddt(kx2/2)
= mvv˙+kxx˙
v=x˙  0 = mx¨+kx.

Similarly power balance could have been used to get the same result, looking at just the mass

P = ddtEK
(kx)(x˙) = ddt(mv2/2)=mvv˙
 0 = kx+mx¨

as before.

Example: Pendulum

Filename:summer95f-5-a
Figure 19.47: A rigid body suspended from a frictionless hinge is an energy conserving one-degree-of-freedom mechanism.

Consider a rigid body with mass m and moment of inertia Io about a hinge which is a distance from the center-of-mass (see fig. 19.47). The familiar simple pendulum is another single degree of freedom system for which the equation of motion can be found from conservation of energy.

ET = constant
0 = ddtET=EK˙+EP˙
= ddt(Ioω2/2)+ddt(mgcosθ)
= Ioωω˙+mg(sinθ)θ˙
ω=θ˙  0 = θ¨+mgIosinθ.

the pendulum equation that we have derived before by this and other means (angular momentum balance about point o).

The above examples are old friends which are handled easily with other techniques. Here is a problem which is much more difficult without the energy method.

Example: Three bars act like a simple pendulum

Filename:pfigure4-2-rp10
Figure 19.48: Three identical uniform bars are pinned with frictionless hinges and swing, obviously, something like a simple pendulum.

Assume all three bars in the structure shown in fig. 19.48are of equal length and have mass m uniformly distributed along their length. It is intuitively obvious that this device swings back and forth something like a simple pendulum. But how can we get the laws of mechanics to tell us this? One approach, which will work in the end, is to draw free-body diagrams of all the parts, write linear and angular momentum balance for each, and then add and subtract equations to eliminate the unknown constraint forces at the various hinges.

The more direct approach is to write the energy equation, adding up the potential and kinetic energies of the parts, all evaluated in terms of the single configuration variable θ. Taking the potential energy to be zero at θ=π/2 (when all centers of mass are at hinge height) we have

ET = constant
0 = ddtET=EP˙+EK˙
= ddt((Ioω2/2)+(Ioω2/2)+(m(ω)2/2))
+ddt((gm(/2)cosθ
gm(/2)cosθ)gmcosθ)
Io=m2/3  0 = ddt(5m2ω2/6)+ddt(2gmcosθ)
= (5m2ωω˙/3+2gmωsinθ)
 0 = 5ω˙/3+2gsinθ
 0 = θ¨+6g5sinθ

which is the same governing equation as for a point-mass pendulum with length 5/6. This is just half way between the following two cases. If the side links had no mass the equation would have been the same as for a point mass pendulum with length

0=θ¨+gsinθ

and if the bottom link had no mass the equation would be the same as a stick hanging from one end which goes back and forth like a point mass pendulum with length 2/3 according to

0=θ¨+3g2sinθ.
Filename:pfigure-blue-125-2
Figure 19.49: A person rides a bike. The pedaling leg is idealized as a pair of equal and opposite forces acting on the seat and pedal.

Example: One on the rim is like two on the frame.

A bicycle transmission is such that the speed of the bike relative to the ground is n times the speed of the pedal relative to the frame :

vbike=nvpedal/bike

Assume the kinetic energy of the relative motion of a rider’s legs can be neglected, as can be the weight of the rider’s leg. At the moment in question the velocity of the pedal is parallel to the direction from the seat to the leg. Thus the free-body diagram of the bike/person system, leaving out the pedaling leg is as shown in fig. 19.49. Let’s assume the bike and rider have mass M and that the wheels have mass mr and mf concentrated on the rim (the hubs are considered part of the frame and the spokes are neglected). Neglecting air resistance etc.the power balance equation is:

P = EK˙ (19.21)

Let’s do some side calculations for evaluating the terms in eqn. (19.21). First, the only forces that do work on the system as drawn are the force on the pedal and the force on the seat.

P = 𝑭P𝒗seat+𝑭P𝒗pedal
= 𝑭P𝒗seat+𝑭P(𝒗bike+𝒗pedal/bike)
= 𝑭P𝒗seat+𝑭P𝒗bike0+𝑭P𝒗pedal/bikeFPvpedal/bike
= FPvbike/n

The net power of the leg is expressed by the compression it carries times its extension rate. The kinetic energy of the wheel comes from both its rotation and its translation. The moment of inertia of a hoop about its center is I=mr2. For rolling contact |ωR|=v so, for one wheel:

EKwheel = mvbike2/2+Iω2/2
= mvbike2/2+(mR2)(v/R)2/2
= mvbike2.

The kinetic energy of a rolling hoop is twice that of a point mass moving at the same speed. Putting these results back in to eqn. (19.21) we have

P = EK˙ (19.22)
FPvbike/n = ddt(Mvbike2/2+(mr+mf)vbike2) (19.23)
= (Mv˙bike+2(mr+mf)v˙bike)vbike (19.24)
 FP/n = (M+2(mr+mf))v˙bike (19.25)
 v˙bike = FPn(M+2(mr+mf)). (19.26)

The bigger the pedal force, the bigger the acceleration, obviously. The higher the gear ratio, the less the acceleration; the faster gears let you pedal slower for a given bike speed, but demand more pedal force for a given acceleration. A heavier bike accelerates less. But the contribution to slowing a bike is twice as much for mass added to the rim as for mass added to the frame or body.

Some comments. n typically ranges from about 1.7 to 8 for a new 21 speed bike and is about 5 for an adult European, Indian or Chinese 1-speed. For a given speed of bicycle riding your feet go n times slower relative to your body than for walking or running at that speed. This calculation is for accelerating a bike on level ground with no wind and rolling resistance. The net speed of a bike in a bike race is not so dependent on weight, because the main enemy is wind resistance. To the extent that weight is a problem it is for steady uphill travel. In this case the mass on the rim makes the same contribution as mass on the frame.

Vibrations

The preponderance of systems where vibrations occur is not due to the fact that so many systems look like a spring connected to a mass, a simple pendulum, or a torsional oscillator. Instead there is a general class of systems which can be expected to vibrate sinusoidally near some equilibrium position. These systems are one-degree -of-freedom (one DOF) near an energy minimum.

In detail why this works out is explained in Box 19.2.

Box 19.3 One degree of freedom systems near a potential energy minimum are harmonic oscillators

In order to specialize to the case of oscillations, we want to look at a one degree of freedom system near a stable equilibrium point, a potential energy minimum.

At a potential energy minimum we have, as you will recall from ‘max-min’ problems in calculus, that dEP(q)/dq=0. To keep our notation simple, let’s assume that we have defined q so that q=0 at this minimum. Physically this means that q measures how far the system is from its equilibrium position. That means that if we take a Taylor series approximation of the potential energy the expression for potential energy can be expressed as follows:

EP const+dEPdq0q+12d2EPdq2Kequivq2+
 dEPdq Kequivq (19.27)

Applying this result to equation( 19.20) we get:

0=Kequivq+12ddq[M(q)]q˙2+M(q)q¨. (19.28)

We now write M(q) in terms of its Taylor series. We have

M(q)=M(0)+dM/dq|0q+ (19.29)

and substitute this result into equation 19.28. We have not finished using our assumption that we are only going to look at motions that are close to the equilibrium position q=0 where q is small. The nature of motion close to an equilibrium is that when the deflections are small, the rates and accelerations are also small. Thus, to be consistent in our approximation we should neglect any terms that involve products of q,q˙,orq¨. Thus the middle term involving q˙2 is negligibly smaller than other terms. Similarly, using the Taylor series for M(q), the last term is well approximated by M(0)q¨, where M(0) is a constant which we will call Mequiv. Now we have for the equation of motion:

0=ddtET  0=Kequivq+Mequivq¨, (19.30)

which you should recognize as the harmonic oscillator equation. So we have found that for any energy conserving one degree of freedom system near a position of stable equilibrium, the equation governing small motions is the harmonic oscillator equation. The effective stiffness is found from the potential energy by Kequiv=d2EP/dq2 and the effective mass is the coefficient of q˙2/2 in the expansion for the kinetic energy EK. The displacement of any part of the system from equilibrium will thus be given by

Asin(λt)+Bcos(λt) (19.31)

with λ2=Kequiv/Mequiv, and A and B determined by the initial conditions. So we have found that all stable non-dissipative one-degree-of-freedom systems oscillate when disturbed slightly from equilibrium and we have found how to calculate the frequency of vibration.

Examples of 1 DOF harmonic oscillators

In the previous section, we have shown that any non-dissipative one-degree-of-freedom system that is near a potential energy minimum can be expected to have simple harmonic motion. Besides the three examples we have given so far, namely,

  • a spring and mass,

  • a simple pendulum, and

  • a rigid body and a torsional spring,

there are examples that are somewhat more complex, such as

  • a cylinder rolling near the bottom of a valley,

  • a cart rolling near the bottom of a valley, and a

  • a four bar linkage swinging freely near its energy minimum.

The restriction of this theory to systems with only one-degree-of-freedom is not so bad as it seems at first sight. First of all, it turns out that simple harmonic motion is important for systems with multiple-degrees-of-freedom. We will discuss this generalization in more detail later with regard to normal modes. Secondly, one can also get a good understanding of a vibrating system with multiple-degrees-of-freedom by modeling it as if it has only one-degree-of-freedom.

Example: Cylinder rolling in a valley

Consider the uniform cylinder with radius r rolling without slip in an cylindrical ‘ideal’ valley of radius R.

For this problem we can calculate EK and EP in terms of θ. Briefly,

EP = mg(Rr)cosθ
EK = 12(32mr2)(θ˙(Rr)r)2
= 34m(Rr)2θ˙2

So we can derive the equation of motion using the fact of constant total energy.

Filename:pfigure-blue-68-1
Figure 19.50: Cylinder rolling without slip in a cylinder. This gives the pendulum equation, but the effective pendulum length does not tend to R as r0 (see text).
0 = ddt(ET)=ddt(EK+EP)
= ddt(mg(Rr)cosθEP+34m(Rr)2θ˙2EK)
= (mg(Rr)sinθ)θ˙+32(Rr)2θ˙θ¨
 0 = mg(Rr)sinθ+32(Rr)2mθ¨

Now, assuming small angles, so θsinθ, we get

g(Rr)θ+32(Rr)2θ¨ = 0 (19.32)
θ¨+(23g(Rr))λ2θ = 0 (19.33)

This is, naturally, our old friend the harmonic oscillator equation. The period is a funny combination of terms. If rR it looks like a point mass pendulum with length 3R/2, more than R. That is, the rolling effect doesn’t go away and make the roller act like a point mass even when the radius goes to zero. See page 17.4 for the angular momentum approach to this problem.

Although, in some abstract way the energy approach always works, practically speaking it has limitations for systems where the configuration is not easily found from one configuration variable.

Example: A four bar linkage analyzed with energy methods

Although probably not usually the best approach, energy methods can be used to find the motions of a 4-bar linkage. Take a four-bar linkage with one bar grounded. Assume the bars all have different lengths. This is a one-DOF system which can use the angle θ of one of the links as a configuration variable. But finding the potential energy as a single formula in terms of all of the links in terms of θ is more trigonometry than most of us like. And then finding the kinetic energy in terms of θ and θ˙ is close enough to impossible that people don’t do it.

So, though it is true that there are functions EP(θ) and EK(θ,θ˙) and that the equations of motion could be written in terms of them, it is really not practical to do so.

How do you find the motions of a 4-bar linkage in practice? It’s more tricky. One approach is to solve the kinematics by integrating kinematic differential equations, as in Sample 18.5 on page 18.5. Then you set up and solve the balance equations of the separate parts as in Sample 19.95 on page 19.95

SAMPLE 19.11

Filename:pfigure-blue-58-1
Figure 19.51:

A plate pendulum. A 2a×2b rectangular plate of mass m hangs from two parallel, massless links EA and FD of length each. The links are hinged at both ends so that when the plate swings, its edges AD and BC remain horizontal at all times. The only driving force present is gravity. Find the equation of motion of the plate.


Solution

Filename:pfigure-blue-157-1
Figure 19.52:

The given system is a single DOF system. So, we need just one configuration variable and the equation of motion will be just one scalar equation in this variable. Let us take angle θ (fig. 19.52) as our configuration variable.

Filename:summer95f-5
Figure 19.53:

The free-body diagram of the plate is shown in fig. 19.53. Note that the link forces 𝑭1 and 𝑭2 act along the links because massless links are two-force bodies. Let (x,y) be the coordinates of the center-of-mass. Then the linear momentum balance for the plate gives

𝑭1+𝑭2mgȷˆ=m𝒂=m(x¨ıˆ+y¨ȷˆ).

Now, we can eliminate both the unknown link forces from this equation by dotting the equation with 𝒏ˆ=cosθıˆ+sinθȷˆ, a unit vector normal to the links. Then, we have

mg(ȷˆ𝒏ˆ) = m[x¨(ıˆ𝒏ˆ)+y¨(ȷˆ𝒏ˆ)]
gsinθ = x¨cosθ+y¨sinθ. (19.34)

Now we need to find a relationship between x and θ, and y and θ, so that we can write x¨ and y¨ in terms of our configuration variable θ and its derivatives. From fig. 19.52, we have

𝒓G = (sinθ+a)xıˆ+(cosθb)yȷˆ
 x¨ = (cosθθ¨sinθθ˙2)
y¨ = (sinθθ¨+cosθθ˙2)

Substituting these expressions for x¨ and y¨ in eqn. (19.34) , we get

gsinθ = θ¨
 θ¨+gsinθ = 0.

Answer: θ¨+gsinθ=0

This is the equation of a simple pendulum! Well, the plate does behave just like a simple pendulum in the given mechanism. From the expressions for the x and y coordinates of the center-of-mass, we have

xa = sinθ
y+b = cosθ
 (xa)2 + (y+b)2=2

that is, the center-of-mass follows a circle of radius centered at (a,b). Since the orientation of the plate never changes (AD and BC always remain horizontal), the plate has no angular velocity. Thus the motion of the plate is equivalent to the motion of a particle of mass m going in a circle centered at (a,b) and driven by gravity. That is the simple pendulum.

SAMPLE 19.12

Filename:pfigure-blue-127-2
Figure 19.54:

Equation of motion from power balance. A slider crank mechanism is shown in fig. 19.54 where the crank is a uniform wheel of mass m and radius R anchored at the center and the connecting rod AB is a massless rod of length . The rod is driven by a piston at B with a known force F(t)=F0cosΩt. There is no gravity. Find the equation of motion of the wheel.


Solution The given mechanism is a one DOF system. So, let us choose a single configuration variable, θ for specifying the configuration of the system and derive an equation that determines θ. Since the applied force is given and the point of application of the force has a simple motion (vertical), it will be easy to calculate power of this force. Also, the connecting rod is massless, so it does not enter into dynamic calculations. The wheel rotates about its center and, therefore, it is easy to calculate its kinetic energy. So, we use the power balance, E˙K=P, here to find the equation of motion of the wheel. Since EK=12Izzcmθ˙2 and P=𝑭𝒗B, we have

Izzcmθ¨θ˙=F(t)ȷˆvBȷˆ=F(t)vB (19.35)
Filename:pfigure-s94h13p4
Figure 19.55:

Now, we need to find vB and express it using the configuration variable θ and its derivatives. There are several ways we could find vB. Vectorially, we could write, 𝒗BvBȷˆ=𝒗A+𝝎AB×𝒓B/A where 𝝎AB=ϕ˙𝒌ˆ. Dotting both sides of this equation with ıˆ and ȷˆ we can find ϕ in terms of θ and vB in terms of θ and θ˙. But, for a change, let us use geometry here.

See fig. 19.55. From triangle ABO, we have

Rsin(90oϕ) = sin(90o+θ)
 cosϕ = Rcosθ
 sinϕϕ˙ = Rsinθθ˙
ϕ˙ = Rsinθsinϕθ˙
Now,
yB = Rsinθsinϕ
 vB y˙B=Rcosθθ˙cosϕϕ˙
= Rθ˙cosθRcosθRsinθsinϕθ˙
= Rθ˙(cosθRsinθcosθ2R2cos2θ)

Substituting this expression for vB in power balance eqn. (19.35), we get

Izzcmθ¨θ˙ = F(t)Rθ˙(cosθsin2θ2(/R)2cos2θ)
θ¨ = RF0cosΩt12mR2(cosθsin2θ2(/R)2cos2θ).

This is the required equation of motion. As is evident, it is a nonlinear ODE which requires numerical solution on a computer if we would like to plot θ(t).

Answer: θ¨=2F0cosΩtmR(cosθsin2θ2(/R)2cos2θ)

SAMPLE 19.13

Filename:pfigure-blue-90-2
Figure 19.56: End A of bar AB is free to slide on the frictionless horizontal surface while end B is going in circles with a disk rotating at a constant rate.

Instantaneous dynamics of slider crank. A uniform rigid rod AB of mass m and length =4R has one of its ends pinned to the rim of a disk of radius R. The other end of the bar is free to slide on a frictionless horizontal surface. A motor, connected to the center of the disk at O, keeps the disk rotating at a constant angular speed ωD. At the instant shown, end B of the rod is directly above the center of the disk making θ to be 30o.

  1. 1.

    Find all the forces acting on the rod.

  2. 2.

    Is there a value of ωD which makes end A of the rod lift off the horizontal surface when θ=30o?


Solution The disk is rotating at constant speed. Since end B of the rod is pinned to the disk, end B is going in circles at constant rate. The motion of end B of the rod is completely prescribed. Since end A can only move horizontally (assuming it has not lifted off yet), the orientation (and hence the position of each point) of the rod is completely determined at any instant during the motion. Therefore, the rod represents a zero degree of freedom system.

  1. 1.

    Forces on the rod: The free-body diagram of the rod is shown in Fig. 19.57. The pin at B exerts two forces Bx and By while the surface in contact at A exerts only a normal force N because there is no friction. Now, we can write the momentum balance equations for the rod. The linear momentum balance (𝑭=m𝒂) for the rod gives

    Filename:s92f1p7
    Figure 19.57: Free-body diagram of the bar.
    Bxıˆ+(By+Nmg)ȷˆ=m𝒂G. (19.36)

    The angular momentum balance about the center-of-mass G (𝑴/G=𝑯˙/G) of the rod gives

    𝒓A/G×Nȷˆ+𝒓B/G×(Bxıˆ+Byȷˆ)=Izz/Gαrod𝒌ˆ. (19.37)

    From these two vector equations we can get three scalar equations (the Angular Momentum Balance gives only one scalar equation in 2-D since the quantities on both sides of the equation are only in the 𝒌ˆ direction), but we have six unknowns — Bx,By,N,𝒂G (counts as two unknowns), and αrod. Therefore, we need more equations. We have already used the momentum balance equations, hence, the extra equations have to come from kinematics.

    𝒗A = 𝒗B+𝝎rod×𝒓A/B𝒗A/B
    or vAıˆ = ωDRıˆ+ωrod𝒌ˆ×(cosθıˆsinθȷˆ)
    = (ωDR+ωrodsinθ)ıˆωrodcosθȷˆ

    Dotting both sides of the equation with ȷˆ we get

    0=ωrodcosθ  ωrod=0.

    Also,

    𝒂A = 𝒂B+𝝎˙×𝒓A/B+𝝎rod𝟎×(𝝎rod×𝒓A/B)𝒂A/B
    or aAıˆ = ωD2Rȷˆ+ω˙rod𝒌ˆ×(cosθıˆsinθȷˆ)
    = (ωD2R+ω˙rodcosθ)ȷˆ+ω˙rodsinθıˆ.

    Dotting both sides of this equation by ȷˆ we get

    ω˙rod=ωD2Rcosθ. (19.38)

    Now, we can find the acceleration of the center-of-mass:

    𝒂G = 𝒂B+𝝎˙×𝒓G/B+𝝎rod𝟎×(𝝎rod×𝒓G/B)𝒂G/B
    = ωD2Rȷˆ+ω˙rod𝒌ˆ×12(cosθıˆsinθȷˆ)
    = (ωD2R+12ω˙rodcosθ)ȷˆ+12ω˙rodsinθıˆ.

    Substituting for ω˙rod from eqn. (19.38) and 30o for θ above, we obtain

    𝒂G=12ωD2R(13ıˆ+ȷˆ).

    Substituting this expression for 𝒂G in eqn. (19.36) and dotting both sides by ıˆ and then by ȷˆ we get

    Bx = 123mωD2R,
    By+N = 12mωD2R+mg (19.39)

    From eqn. (19.37)

    12[(ByN)cosθBxsinθ]𝒌ˆ = 112m2(ωD2Rcosθ)𝒌ˆ
    or ByN = 16mωD2Rcos2θ+Bxtanθ (19.40)
    = 29mωD2R16mωD2R
    = 718mωD2R

    From eqns. (19.39) and (19.40)

    By=12(mg89mωD2R)

    and

    N=12(mg19mωD2R).
  2. 2.

    Lift off of end A: End A of the rod loses contact with the ground when normal force N becomes zero. From the expression for N from above, this condition is satisfied when

    29mωD2R = mg
     ωD = 3gR.

SAMPLE 19.14

Filename:twodisks-ang-mom
Figure 19.58:

A bar sliding on a sliding wedge. A bar AB of mass m and length is hinged at end A and rests on a wedge of mass M at the other end B. The contact at B is frictionless. The wedge is free to slide horizontally without any friction. The motion of the system is driven only by gravity. Find the equation of motion of the bar using

  1. 1.

    momentum balance

  2. 2.

    energy (or power) balance.


Solution

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

(a) Momentum Balance: The bar and the wedge make up a single DOF system. To derive the equations of motion of the bar, let us choose θ as the configuration variable. The free-body diagram of the bar and the wedge are shown in fig. 19.59. Note that the normal reaction at B is normal to the wedge surface, i.e., 𝐍=N𝐧ˆ. Now, the angular momentum balance for the bar about point A gives,

𝑯˙A = 𝑴A
IzzAθ¨𝒌ˆ = 2𝒆ˆR×(mgȷˆ)+𝒆ˆR×N𝒏ˆ (19.41)
= 12mgcosθ𝒌ˆ+Ncos(αθ)𝒌ˆ

where the last line follows from the fact that 𝒆ˆR=cosθıˆ+sinθȷˆ, the unit normal 𝒏ˆ=sinαıˆ+cosαȷˆ, and so, 𝒆ˆR×ȷˆ=cosθ𝒌ˆ and 𝒆ˆR×𝒏ˆ=cos(αθ)𝒌ˆ. We now need to eliminate the unknown normal reaction N from the above equation. Since the wedge is constrained to move only horizontally, we can write the linear momentum balance for the wedge as

N𝒏ˆıˆ=Mx¨  N=Mx¨𝒏ˆıˆ=Mx¨sinα (19.42)
Filename:pfigure-blue-110-1
Figure 19.60:

Thus, we have found N in terms of x¨ that we can use in eqn. (19.41) to get rid of N. But, we now need to express x¨ in terms of our configuration variable θ and its derivatives. Consider the triangle ABC formed by the bar and the slanted edge of the wedge. Let x= AC denote the horizontal position of the wedge. Then, from the law of sines, we have xsin(αθ)=sinα, so that

x = sinαsin(αθ)
 x˙ = sinαcos(αθ)(θ˙) (19.43)
 x¨ = sinα[θ¨cos(αθ)+θ˙2sin(αθ)]. (19.44)

Now, substituting for N in eqn. (19.41) from eqn. (19.42), using the expression for x¨ from above, and dotting the resulting equation with 𝒌ˆ, we get

13m2θ¨ = mg2cosθ+M2cos(αθ)sin2α[θ¨cos(αθ)+θ˙2sin(αθ)]
 θ¨ = 3mgsin2αcosθ+3Mθ˙2sin2(αθ)2msin2α+6Mcos2(αθ) (19.45)
= 3[(g/)sin2αcosθ+(M/m)θ˙2sin2(αθ)]2[sin2α+3(M/m)cos2(αθ)]

Answer: θ¨=3[(g/)sin2αcosθ+(M/m)θ˙2sin2(αθ)]2[sin2α+3(M/m)cos2(αθ)]

(b) Power balance: Now we derive the equation of motion for the bar using power balance E˙K=P. For power balance, we have to consider both the bar and the wedge. The bar rotates about the fixed point A, therefore, its kinetic energy is (1/2)IzAzθ˙2. The wedge moves with horizontal speed x˙, therefore, its kinetic energy is (1/2)Mx˙2. Thus, EK=1/2)IAzzθ˙2+(1/2)Mx˙2. The only force that contributes to power is the force of gravity on the rod because 𝒗G(mgȷˆ) is non-zero. The sliding contact at B is frictionless and hence the net power due to the contact force there is zero. Now, 𝒗G=12θ˙𝒆ˆθ. So, P=12mgθ˙(𝒆ˆθȷˆ)=12mgθ˙cosθ. Thus the power balance for the system gives

IzzAθ¨θ˙+Mx¨x˙=12mgθ˙cosθ. (19.46)

We can simplify this equation further. Note that, θ¨θ˙=ddt(12θ˙2) and x¨x˙=ddt(12x˙2). So that we can write the above equation as

ddt(12IzzAθ˙2)+ddt(12Mx˙2) = 12mgθ˙cosθ
or d(IzzAθ˙2+Mx˙2) = mgcosθdθ
 IzzAθ˙2+Mx˙2 = Cmgsinθ

where C is a constant of integration to be determined from initial conditions. For example, when the bar begins to slide from rest, we have, at t=0, x˙=0 and θ˙=0. At that instant, if θ(0)=θ0, then C=mgsinθ0. So, we can write

IzzAθ˙2+Mx˙2=mg(sinθ0sinθ).

Now, replacing x˙ in this equation with the expression we obtained in eqn. (19.43), we have

IzzAθ˙2+M2θ˙2cos2(αθ)sin2α = mg(sinθ0sinθ)
 θ˙2 = mg(sinθ0sinθ)IzzA+M2cos2(αθ)sin2α
= mgsin2α(sinθ0sinθ)(1/3)m2sin2α+M2cos2(αθ)
 θ˙ = 3gsinαsinθ0sinθsin2α+3(M/m)cos2(αθ).

This is a first order, nonlinear, ODE compared to the second order equation we got in eqn. (19.45). However, again, we need to resort to numerical solution if we wish to solve for θ(t), in which case, this reduction to the first order equation does not save much work.

Answer: θ˙=3gsinαsinθ0sinθsin2α+3(M/m)cos2(αθ)

Problems for 19.2 One-degree-of-freedom 2-D mechanisms

19.2.1  A conservative vibratory system has the following equation of conservation of energy.

m(θ˙)2mg(1cosθ)+K(aθ)2=E0

where E0 is a constant.

  1. (a)

    Obtain the differential equation of motion of this system by differentiating this energy equation with respect to t.

  2. (b)

    Determine the circular frequency of small oscillations of the system in part (a). HINT: (Let sinθ and cosθ1θ2/2).

19.2.2  A motor at O turns at rate ωo whose rate of change is ω˙o. At the end of a stick connected to this motor is a frictionless hinge attached to a second massless stick. Both sticks have length L. At the end of the second stick is a mass m. For the configuration shown, what is θ¨? Answer in terms of ωo, ω˙o, L, m, θ, and θ˙. Ignore gravity.

Filename:pfigure-blue-107-1
Figure 19.61

19.2.3

In problem 18.3, find

  1. (a)

    the angular momentum about point O, Answer: H/O=[mL2ω1+12mr2(ω1+ω2)]kˆ.

  2. (b)

    the rate of change of angular momentum of the disk about point O, Answer: H˙/O=𝟎.

  3. (c)

    the angular momentum about point C, Answer: H/C=12mr2(ω1+ω2)kˆ.

  4. (d)

    the rate of change of angular momentum of the disk about point C. Answer: H˙/C=𝟎.

Assume the rod is massless and the disk has mass m.

Filename:pfigure-s94h14p4
Figure 19.62

19.2.4  Robot arm, 2-D . The robot arm AB is rotating about point A with ωAB=5rad/s and ω˙AB=2rad/s2. Meanwhile the forearm BC is rotating at a constant angular speed with respect to AB of ωBC/AB=3rad/s. Gravity cannot be neglected. At the instant shown, find the net force acting on the object P which has mass m=1kg. Answer: F=ma=109.3Nıˆ19.54Nȷˆ.

Filename:pfigure-blue-112-1
Figure 19.63

19.2.5  A crude model for a column, shown in the figure, consists of two identical rods of mass m and length with hinge connections, a linear torsional spring of stiffness K attached to the center hinge (the spring is relaxed when θ=0), and a load P applied at the top end.

  1. (a)

    Obtain the exact nonlinear equation of motion.

  2. (b)

    Obtain the squared natural frequency for small motions θ.

  3. (c)

    Check the stability of the straight equilibrium state θ=0 for all P0 via minimum potential energy. How do these results compare with those from part(b)?

Filename:tfigure8-alt-app2c
Figure 19.64

19.2.6  A thin uniform rod of mass m rests against a frictionless wall and on a frictionless floor. There is gravity.

  1. (a)

    Draw a free-body diagram of the rod.

  2. (b)

    The rod is released from rest at θ=θ00. Write the equation of motion of the rod. Answer: θ¨=3g2Lcosθ.

  3. (c)

    Using the equation of motion, find the initial angular acceleration, 𝝎˙AB, and the acceleration of the center of mass, 𝒂G, of the rod. Answer: From the diagram, we see 𝝎˙AB=θ¨(t=0)𝒌ˆ=3g2Lcosθ0,

    𝒂G=34gcosθ0[sinθ0ıˆcosθ0ȷˆ].

  4. (d)

    Find the reactions on the rod at points A and B. Answer: RA=34mgsinθ0cosθ0ıˆ, 𝑹B=mg[134gcos2θ0]ȷˆ.

  5. (e)

    Find the acceleration of point B. Answer: aB=32gsinθ0cosθ0ıˆ.

  6. (f)

    When θ=θ02, find 𝝎AB and the acceleration of point A. Answer: At θ=θ02, 𝝎AB=3gL(sinθ0sinθ02)𝒌ˆ, 𝒂A=3g[12cos2θ02+sinθ02(sinθ0sinθ02)]ȷˆ.

Filename:sfig4-6-4a
Figure 19.65

19.2.7  Bar leans on a crooked wall. A uniform 3lbm bar leans on a wall and floor and is let go from rest. Gravity pulls it down.

  1. (a)

    Draw a Free-Body Diagram of the bar.

  2. (b)

    Kinematics: find the velocity and acceleration of point B in terms of the velocity and acceleration of point A

  3. (c)

    Using equations of motion, find the acceleration of point A.

  4. (d)

    Do the reaction forces at A and B add up to the weight of the bar? Why or why not? (You do not need to solve for the reaction forces in order to answer this part.)

Filename:sfig4-6-3
Figure 19.66

19.2.8  Two blocks, each with mass m, slide without friction on the wall and floor shown. They are attached with a rigid massless rod of length that is pinned at both ends. The system is released from rest when the rod makes an angle of 45. What is the acceleration of the block B immediately after the system is released?

Filename:sfig4-6-3a
Figure 19.67

19.2.9  Slider Crank. 2-D . No gravity. Refer to the figure in problem 18.78. What is the tension in the massless rod AB (with length L) when the slider crank is in the position with θ=0 (piston is at maximum extent)? Assume the crankshaft has constant angular velocity ω, that the connecting rod AB is massless, that the cylinder walls are frictionless, that there is no gas pressure in the cylinder, that the piston has mass M and the crank has radius R.

Filename:pfigure4-4-rp12
Figure 19.68

19.2.10  In problem 18.84, find the force on the wheel at point D due to the rod.

Filename:pfigure-blue-38-2
Figure 19.69

19.2.11  In problem 18.4, find the force on the rod at point P.

Filename:pfigure-blue-49-2
Figure 19.70

19.2.12  Slider-Crank mechanism. The slider crank mechanism shown is used to push a 2lbm block P. Arm AB and BC are each 0.5ft long. Given that arm AB rotates counterclockwise at a constant 2 revolutions per second, what is the force on P at the instant shown? Answer: There are many solution methods. 8π2 poundals to the left. (Note: 1 poundal =1ftlbms2.)

Filename:pfigure4-4-rp13
Figure 19.71

19.2.13  The large masses m at C and A were supported by the light triangular plate ABC, whose corners follow the guide. B enters the curved guide with velocity V. Neglecting gravity, find the vertical reactions (in the y direction) at A and B. (Hint: for the rigid planar body ABC, find y¨A and y¨B in terms of aycm and θ¨, assuming θ˙=0 initially.)

Filename:f92h7p1
Figure 19.72

19.2.14  An idealized model for a car comprises a rigid chassis of mass Mc and four identical rigid disks (wheels) of mass Mw and radius R, as shown in the figure. Initially in motion with speed v0, the car is momentarily brought to rest by compressing an initially uncompressed spring of stiffness k. Assume no frictional losses.

  1. (a)

    Assuming no wheel slip, determine the compression Δ of the spring required to stop the car. Answer: Δ=(Mc+6Mw)v02k.

  2. (b)

    While the car is in contact with the spring but still moving forward, in which direction is the tangential force on any one of the wheels (due to contact with the ground) pointing? Why? (Illustrate with a sketch). Answer: Tangential forces point toward the wall!

  3. (c)

    Repeat part (a) assuming now that the ground is perfectly frictionless from points A to B shown in the figure. Answer: Δ=(Mc+4Mw)v02k.

Filename:pfigure-blue-52-2
Figure 19.73

19.2.15  Assume Greg Lemond’s riding “tuck” was so good that you can neglect air resistance when you think about him and his bike. Further, you can regard his and the bike’s combined mass (all 70 kg)) as concentrated at a point in his stomach somewhere. Greg’s left foot has just fallen off the pedal so he is only pedaling with his right foot, which at the moment in question is at its lowest point in the motion. You note that, relative to the ground the right foot is only going 3/4 as fast as the bike (since it is going backwards relative to the bike), though you can’t make out all the radii of his frictionless gears and rigid round wheels. Greg, ever in touch with his body, tells you he is pushing back on the pedal with a force of 70 N. You would like to know Greg’s acceleration.

  1. (a)

    In your first misconceived experiment you set up a 70 kg bicycle in your laboratory and balance it with strings that cause no fore or aft forces. You tie a string to the vertically down right pedal and pull back with a force of 70 N. What acceleration do you measure?

  2. (b)

    What is Greg’s actual acceleration? [Hint: Greg’s massless leg is pushing forward on his body with a force of 70 N] You can neglect the mass of the wheels and other transmission parts (chain, crank, etc).

Filename:pfigure4-4-rp16
Figure 19.74

19.2.16  Which way does the bike accelerate? A bicycle with all frictionless bearings is standing still on level ground. A horizontal force F is applied on one of the pedals as shown. There is no slip between the wheels and the ground. The bicycle is gently balanced from falling over sideways. It is heavy enough so that both wheels stay on the ground. Does the bicycle accelerate forward, backward, or not at all? Make any reasonable assumptions about the dimensions and mass. Justify your answer as clearly as you can, clearly enough to convince a person similar to yourself but who has not seen the experiment performed.

Filename:pfigure4-4-rp14
Figure 19.75

19.3 Dynamics of multi-degree-of-freedom 2-D mechanisms

A typical machine has many parts. They may work in concert as a one degree-of-freedom system or they may be designed to move independently.

One degree-of-freedom machines. Some mechanisms are designed so the parts work together to do something, one thing, well. A car going straight down a road has dozens of moving parts: cylinders, connecting rods, crank shaft, cam shaft, transmission gears, drive shaft, differential, wheel axles, wheels and the car itself. All of these parts move, each in its own different way, approximately like a rigid object. But in the simplest description

margin: A car has one degree-of-freedom? For studying the acceleration of a car going straight one might ignore the flexibility of the drive shaft and the bounce of the suspension, etc.thus leading to a 1 DOF model. But if you want to see how acceleration causes the car to tip back, or how the car oscillates when the engine lugs, then a 2, or more, DOF model would be appropriate.

they move as a one-degree-of freedom mechanism towards the end of fighting friction and moving the car down the road. If the crankshaft rotates a small amount, each part moves a corresponding amount; the amount of rotation of the crank shaft determines the position of all of the parts. And all of the velocities and accelerations of all of the bits of mass in the car can be determined by the rotation θcrank and its derivatives ( dotθcrank and ddotθcrank). Thus a car, even with lots of pieces moving this way and that, might be well-modeled for some purposes as a one degree-of-freedom mechanism. For one-degree-of-freedom mechanisms the methods of Section  19.2 may be appropriate.

Multi DOF systems In some modeling of machines it is important to keep track of multiple degrees of freedom (See fig. 18.65 on page 18.65 for examples of simple mechanisms). There are various reasons that a multi-DOF analysis is relevant:

  • Some machines intentionally have various ways of moving that are controlled by various motors. Classic examples include:

    • Robots,

    • Animal bodies.

  • Some systems intentionally have various degrees of freedom so as to respond smoothly to disturbances, for example the suspensions on the bottoms of cars and washing machines.

  • Some systems have undesirable degrees of freedom due to the parts which are nominally rigid actually being elastic. Thus a machine which is intended to have one degree of freedom might have various undesirable vibration modes.

Mechanical analysis of multi-DOF systems. To solve problems with multiple degrees of freedom the basic strategy is, as described at the start of the chapter:

  1. 1.

    Draw FBDs of each object,

  2. 2.

    Pick configuration variables,

  3. 3.

    Write linear and angular momentum balance equations

  4. 4.

    Solve the equations for variables of interest (usually forces and second derivatives of the configuration variables).

  5. 5.

    Set up and solve the resulting differential equations (if you are trying to find the motion).

There are two basic approaches to these multi-object problems which, for lack of better language we call “brute-force” and “clever”.

  • A.

    In the brute force approach you write three times as many scalar balance equations as you have objects (limiting attention to 2D where each object has 3 DOFs and 3 independent momentum balance equations). That is, for example, for each free-body diagram you write linear momentum balance and angular momentum balance about the center of mass. Then you take this set of 3n equations and add and subtract them to solve for variables of interest.

    This approach is quite suitable for computers so most commercial general purpose dynamic simulators use a variant of this approach. For individual use with packaged software, the brute-force approach is generally both more reliable and more time consuming.

  • B.

    In the clever approach you write as many scalar momentum balance equations as you have unknowns. For example, if you have 2 degrees of freedom and you are concerned with motions and not reaction forces, you write 2 equations. You do this by finding momentum balance equations that do not include the variables you are not interested in. Usually this involves using angular momentum balance about hinge points, or linear momentum balance orthogonal to sliding contacts.

    The clever approach does not always work; the four-bar linkage is the classic problem case. However the desire to find such minimal sets of equations of motion it is historically important

    margin: Attempts to automate the clever approach, that is to quickly find minimal equations of motion, led to Lagrange equations which led to Hamilton’s equations which led to quantum mechanics (but we won’t be that clever here).

    .

At this point in the subject there are no quick simple problems. All problems are involved, especially if taken from start all the way to plotting solutions to the differential equations. The examples that follow emphasize getting to the equations of motion. The skills for numerically solving the differential equations and plotting the solutions are the same as from the start of dynamics so are not discussed until the sample problems which show all the work from setup to solution.

Example: Block sliding on sliding block: clever approach

Block 1 with mass m1 rolls without friction on ideal massless wheels at A and B (see fig. 19.76). Block 2 with mass m2 rolls down the tipped top of block 1 on ideal massless rollers at C and D. The locations of G1 relative to A and B, and of G2 relative to C and D are known. How do blocks 1 and 2 move?

First look at the free-body diagram of the system and note that there are no unknown forces in the ıˆ direction. So, for the system

{𝑭i = 𝑳˙}ıˆ
 0 = m1x¨+m2(x¨ıˆ+y¨ȷˆ)ıˆ (19.47)
= (m1+m2)x¨y¨m2cosθ.

Looking at the free-body diagram of mass 2 note that there are no unknown forces in the ȷˆ direction, so

Filename:pfigure4-4-rp15
Figure 19.76: Block 2 rolls on block 1 which rolls on the ground. All rollers are ideal (frictionless and massless). The moving ıˆȷˆ frame moves with the lower block and is oriented with the slope.
{𝑭i = 𝑳˙}ȷˆ
 m2gsinθ = m2(x¨ıˆ+y¨ȷˆ)ȷˆ
m2gsinθ = m2(cosθx¨+y¨). (19.48)

Eqns. 19.47 and 19.48 are a system of two equations in the two unknowns x¨ and y¨

(m1+m2)x¨m2cosθy¨=0m2cosθx¨+m2y¨=m2gsinθ

which can be solved for x¨ and y¨ by hand or on the computer. Finding x(t) and y(t) is then easy because both x¨ and y¨ are constants.

Now we look at the same example, but proceed in a more naive manner.

Example: Block sliding on sliding block: brute force approach

Now we look at the free-body diagrams of the two separate blocks. We will use 3 balance equations from each free-body diagram, taking account of the kinematic constraints.

For the lower block we have

AMBG1 𝑴/G1=𝑯˙/G1
where𝑴/G1=𝒓D/G1×(FDıˆ)
+𝒓C/G1×(FCıˆ)
+𝒓B/G1×FBȷˆ
+𝒓A/G1×FAȷˆ
and𝑯˙/G1=𝟎 ( 𝝎1=𝟎)
andLMB 𝑭i=𝑳˙
where𝑭i=FDıˆFCıˆ+FBȷˆ+FAȷˆm1gȷˆ
and𝑳˙=m1x¨ıˆ

Similarly for the upper block:

AMBG2 𝑴/G2=𝑯˙/G2
where𝑴/G1=𝒓D/G2×FDıˆ
+𝒓C/G2×FCıˆ
and𝑯˙/G2=𝟎 ( 𝝎1=𝟎)
andLMB 𝑭i=𝑳˙
where𝑭i=FDıˆ+FCıˆm2gȷˆ
and𝑳˙=m2(x¨ıˆ+y¨ȷˆ)

Eqns. 19.3-19.3 can be written as scalar equations by dotting the LMB equations with ıˆ and ȷˆ and the AMB equations with 𝒌ˆ. All is known in these 6 equations but the six scalars: FA,FB,FC,FD,x¨, and y¨. These could be set up as a matrix equation and solved on the computer, or you could try to find your way through by adding and subtracting equations. In any case you could solve for x¨ and y¨ and thus have differential equations to solve to find the motions.

One quick inference one can make is from looking at the equations: no term in the coefficients of the unknowns depends on x,x˙,y, or y˙. So all of the reactions FA,FB,FC, and FD as well as the accelerations x¨ and y¨ are constants in time (until the upper mass hits the ground).

Example: Block sliding on sliding block: even more brute force

This “multi-body” problem can be solved in an even more naive and more brute force manner. The method is the same as shown in Section 14.1 on page 14.1 for a one dimensional problem.

We would use 6 configuration variables: the x and y coordinates of G1 and G2 and the rotations of the two bodies:

x1,y1,θ1,x2,y2, and θ2

That the two bodies don’t rotate would be expressed indirectly by noting that the accelerations of points A and B on the lower mass must be in the ıˆ direction. These two equations would be added to the 6 linear and angular momentum balance equations. Similar constraint equations would be written for the interactions at C and D. Altogether there are now 6 configuration variables and 4 constraint forces. But there are 6 differential equations of motion and 4 constraint equations. Thus at one instant in time a set of 10 simultaneous equations needs to be solved. Then these are used to evaluate the right hand sides in the differential equations.

These 10 equations are an impractical mess for solving one simple problem like this. But these equations lend themselves to easy automation and this is closest to the approach used by general purpose dynamics simulators.

The example above is particularly simple because block 1 moves in a straight line without rotating and block 2 moves in a straight line without rotating relative to block 1. Even for a system of just two objects the situation could be much more complex if the first object had a complex motion and the second a complex motion relative to the first. But because of the preponderance of hinges in the world, circular motion, and motion relative to circular motion, is the most complex motion that need be considered by many engineers. Here is a version of the most common example of that class.

Example: A two link robot arm: finesse the finding of reactions

Filename:tfigure4-spherical-rotaxis
Figure 19.77: A two link robot arm. Free-body diagrams are shown of the whole system (including the motor torque at the shoulder Ms) and of the fore-arm (including the motor torque at the elbow Me) .

A robot arm has two links. There are motors that apply known torque Ms at the shoulder (reacted by the base and a torque Me at the elbow (reacted by the upper arm). Dimensions are as marked. This system has 2 degrees of freedom. So we need 2 configuration variables and 2 independent balance equations to find the motion.

The natural configuration variables are the angles of the upper arm relative to a fixed reference and the angle of the lower arm relative to a fixed reference. It would also be natural to instead use the angle of the lower arm relative to the upper arm. This leads to simpler equations in the end, but more work in set up.

Angular momentum balance for the system about the shoulder contains no unknown reaction forces, nor does angular momentum balance of the fore-arm about the hinge. So we base our work on these two equations:

System AMB/O  𝑴/O = 𝑯˙/O (19.53)
Forearm AMB/E  𝑴/E = 𝑯˙/E. (19.54)

The goal, equations of motion, is reached by evaluating the left and right sides of these equations in terms of known geometric and mass quantities as well as the configuration variables When we write 𝑴/O and 𝑯˙/O we implicitly mean for the whole system. And 𝑴/E and 𝑯˙/E apply just to the forearm.

At each step in the calculations below imagine the results can be substituted into the later steps. We don’t do that here because the expressions grow in size. Further, if the angles and their rates of change are known, as they are when doing most dynamics problems, the intermediate calculations will result in numbers, and these do not become more numerous as the calculation proceeds.

𝝀ˆ1=cosθ1ıˆ+sinθ1ȷˆ and 𝝀ˆ2=cosθ2ıˆ+sinθ2ȷˆ
𝒓G1/O=1𝝀ˆ1,𝒓E/O=3𝝀ˆ1 and 𝒓G2/E=2𝝀ˆ2
and 𝒓G2/O=𝒓E/O+𝒓G2/E.
𝒂G1/O=θ˙12𝒓G1/O+θ¨1𝒌ˆ×𝒓G1/O ,
𝒂E/O=θ˙12𝒓E/O+θ¨1𝒌ˆ×𝒓E/O ,
𝒂G2/E=θ˙22𝒓G2/E+θ¨2𝒌ˆ×𝒓G2/E and 𝒂G2/O=𝒂E/O+𝒂G2/E.

These terms are all we need to evaluate the 4 terms in Eqns. 19.53 and 19.54.

𝑴/O=𝒓G1/O×(m1gıˆ)+𝒓G2/O×(m2gȷˆ)+Ms𝒌ˆ,
𝑴/E=𝒓G2/E×(m2gıˆ)+Me𝒌ˆ,
𝑯˙/O=m1𝒓G1/O×𝒂G1/O+θ¨1I1𝒌ˆ+m2𝒓G2/O×𝒂G2/O+θ¨2I2𝒌ˆ,
and 𝑯˙/E=m2𝒓G2/E×𝒂G2/O+θ¨2I2𝒌ˆ.

Once these are substituted into Eqns. 19.53 and 19.54 one has 2 vector equations with only 𝒌ˆ components. In other words we have two scalar equations in the two unknowns θ¨1 and θ¨2. One can go through the algebra and solve for them explicitly, but the expressions are quite complex, even when simplified. At given values of θ1,θ˙1,θ2, and θ˙2 however these are just two linear equations in two unknowns.

Closed kinematic chains

When a series of mechanical links is open you can not go from one link to the next successively and get back to your starting point. Such chains include a pendulum (1 link), a double pendulum (2 links), a 100 link pendulum, and a model of the human body (so long as only one foot is on the ground). A closed chain has at least one loop in it. You can go from link to next and get back to where you started. A slider-crank, a 4-bar linkage, and a person with two feet on the ground are closed chains.

Closed chains are kinematically difficult because they have fewer degrees of freedom than they have joints. So some of the joint angles depend on the others. The values of any minimal set of configuration variables, say some of the joint angles, determines all of the joint angles, but by geometry that is difficult or impossible to express with formulas.

Example: Four bar linkage.

It is impractically difficult to write the positions velocities and accelerations of a 4-bar linkage in terms of θ, θ˙ and θ¨ of any one of its joints. Why? Because finding all of the bar angles from one angle is a trigonometric mess. And differentiating that mess once (for velocities) and then once again (for accelerations) makes a huge mess.

SAMPLE 19.15

Filename:tfigure5-7
Figure 19.78:

Dynamics of sliding wedges. A wedge shaped body of mass m2 sits on a frictionless ground. Another wedge shaped body of mass m1 is gently placed on the inclined face of the stationary wedge. The top wedge starts to slide down. The coefficient of friction between the two wedges is μ. Find the sliding acceleration of the top wedge along the incline (i.e., the relative acceleration of m1 with respect to m2).


Solution

Filename:tfigure5-gen-rigid-body
Figure 19.79:

The free-body diagrams of the two wedges are shown in fig. 19.79. Note that the friction force is μN since the wedges are sliding with respect to each other (if they were not sliding already then the friction force is an unknown force FμN). Let the absolute acceleration of m2 be 𝒂2=a2ıˆ. Then, the absolute acceleration of m1 is 𝒂1=𝒂2+𝒂1/2=a2ıˆ+arel𝝀ˆ. Now, we can write the linear momentum balance for m1 and m2 as follows.

N𝒏ˆm1gȷˆμN𝝀ˆ = m1(a2ıˆ+arel𝝀ˆ) (19.55)
(Rm2g)ȷˆN𝒏ˆ+μN𝝀ˆ = m2a2ıˆ (19.56)

where 𝝀ˆ=cosαıˆsinαȷˆ and 𝒏ˆ=sinαıˆ+cosαȷˆ. Here, we have 4 independent scalar equations (from the two 2-D vector equations) in four unknowns N,R,a2, and arel. Thus, we can certainly solve for them. We are, however, only interested in arel. So, we should try to find the answer with fewer calculations. Dotting eqn. (19.55) with 𝝀ˆ, we have

m1arel=m1a2cosαμN+m1gsinα (19.57)

So, to find arel, we need a2 and N. Dotting eqn. (19.55) with 𝒏ˆ, we have

m1a2sinα=Nm1gcosα, (19.58)

and dotting eqn. (19.56) with ıˆ, we have

m2a2=Nsinα+μNcosα. (19.59)

Solving eqn. (19.58) and (19.59) simultaneously, and using new variables (for convenience) M=m1/m2,C=cosα, and S=sinα, we get

a2=MC(μCS)1MS(μCS)g,N=m1gC1MS(μCS)

Substituting these expression in eqn. (19.57), we get

arel=gSgC[MC(μCS)μ]1MS(μCS)

Answer: arel=gsinαgcosα[m1m2cosα(μcosαsinα)μ]1m1m2sinα(μcosαsinα)

Note that when there is no friction (μ=0), the expression for arel reduces to

arel=gS+gMC2S1+MS2

and if we let m2 (i.e., m2 represents fixed ramp) so that M0, then arel=gsinα which is the acceleration of a point mass down a frictionless ramp of slope tanα.

SAMPLE 19.16

Filename:tfigure5-term1-a
Figure 19.80:

Dynamics of a new gun. A new gun consists of a uniform rod AB of mass m1 and a small collar C of mass m2 that slides freely on the rod. A motor at A rotates the rod with constant torque T.

  1. 1.

    Find the equations of motion of the collar.

  2. 2.

    Show that if T=0 then the equations of motion imply conservation of angular momentum about point A.


Solution

Filename:tfigure5-term1-b
Figure 19.81:
  1. 1.

    Let us denote the configuration of the collar with R, the radial distance from the fixed point A along the rod, and θ, the angular displacement of the rod. We need to find differential equations that determine R and θ as functions of time. The free-body diagram of the whole system (rod and collar together) and that of the collar is shown in fig. 19.81. We can write angular momentum balance for the whole system about point A so that the unknown reaction force F at A does not enter the equations. Noting that the acceleration of the collar is 𝒂C=(R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨)𝒆ˆθ, and letting I1IzAz be the moment of inertia of the rod about A, we have,

    𝑴A = 𝑯˙A
    T𝒌ˆ = I1θ¨𝒌ˆ+R𝒆ˆR×m2[(R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨)𝒆ˆθ]
    = I1θ¨𝒌ˆ+m2R2θ¨𝒌ˆ+2m2RR˙θ˙𝒌ˆ

    Dotting this equation with 𝒌ˆ, we have

    θ¨=TI1+m2R22m2RI1+m2R2R˙θ˙. (19.60)

    Thus we have obtained the equation of motion for θ. Now we consider the free-body diagram of the collar alone and write the linear momentum balance for it in the 𝒆ˆR direction, i.e., 𝐞ˆR(𝐅=m𝐚), so that we do not have to care about the unknown normal reaction 𝑵. So, we have,

    0 = 𝒆ˆRm2[(R¨Rθ˙2)𝒆ˆR+(2R˙θ˙+Rθ¨)𝒆ˆθ]
    = m2(R¨Rθ˙2)
     R¨ = Rθ˙2. (19.61)

    Thus we have the required equations of motion. Note that eqn. (19.60) and (19.61) are coupled nonlinear differential equations. So, to find θ(t) and R(t) we need to solve them numerically.

    Answer: θ¨=TI1+m2R22m2RI1+m2R2R˙θ˙,R¨=Rθ˙2

  2. 2.

    Now we set T=0 in our equations of motion. Note that the equation for R is independent of T. The equation for θ becomes

    θ¨=2m2RI1+m2R2R˙θ˙.  (I1+m2R2)θ¨+2m2RR˙θ˙=0

    But the last expression is simply H˙A for the system. Thus we have H˙A=0 which implies that HA= constant. That is conservation of angular momentum about point A.

SAMPLE 19.17  Numerical solutions of new gun equations. Consider Sample 19.79 again. Set up the equations of motion for numerical solution. Take T=1Nm,=1m,m2=1kg, and m1=m2/3. Carry out numerical solutions for the following cases.

  1. 1.

    Let the system start from rest at θ=0 and R=0.1m. Find the solution from t=0 to t=1s. Plot R(t),θ(t) and R(θ) (in polar coordinates).

  2. 2.

    Find the solution till the collar leaves the rod. What is the speed of the collar at this instant?

  3. 3.

    Compute and plot the total energy of the system as a function of time. Also, plot the work done by the torque as a function of time and show that the work done is equal to the total energy of the system at each instant.

  4. 4.

    Vary torque T and carry out solutions for several values of T. Find the terminal value of θf (when the collar leaves the rod) for each T. Justify your observation about θf by plotting R˙/θ˙ as a function of T.

Solution

We first need to write the equations of motion, eqn. (19.60) and (19.61), as a set of first order ODEs. We can easily do so by introducing new variables ωθ˙ and vRR˙, so that we have,

(θ˙ω˙R˙v˙R)=(ωTI1+m2R22m2RI1+m2R2vRωvRRω2)

Given the values of all constants, we only need to specify the initial conditions for θ,ω,R, and vR for solving these equations numerically.

  1. 1.

    We use the following pseudocode to carry out the numerical solution.

         Set T = 1, L = 1, m2 = 1, m1 = m2/3
         Let I1 = m1*L^2/3, I2 = m2*R^2,
         ODEs = {thetadot = w,
                 wdot = (T-2*m2*R*vR*w)/(I1+I2),
                 Rdot = vR,
                 vRdot = R*w^2}
         IC = {theta = 0, w = 0, R = 0.1, vR = 0}ΨΨ
         Solve ODEs with IC for t=0 to t=1
         Plot t vs R, Plot t vs theta,
         Polarplot theta vs R
    
    Filename:sfig1-2-12
    Figure 19.82:
    Filename:sfig4-6-5
    Figure 19.83:

    The R(t) and θ(t) plots obtained from the numerical solution are shown in fig. 19.83 and the polar plot of R(θ) is shown in fig. 19.82.

  2. 2.

    We do not know a priori the value of t at which the collar leaves the rod. So, we have to carry out the solution for some assumed tf which gives us R(tf)> so that we know the collar has gone past the end of the rod. We then plot R(t), including the unreal value of R(tf)>, and find the time t at which R(t)=, either by zooming into the graph or by interpolation (although, there are various sophisticated algorithms to find this t). Following the method of zooming into the graph

    (see fig. 19.84)

    we find the terminal value of t to be 1.147s. We carry out the numerical solution again from t=0 to tf=0.147s and find that R(tf)=1m,vR(tf)=2.13m/s, and ω(tf)=1.03rad/s, so that vf=R˙2+(Rθ˙)2=2.37m/s. This is the terminal speed of the collar.

    Filename:sfig4-6-5a
    Figure 19.84: Finding the time t at which the collar leaves the rod from the graph of R(t).
    Filename:sfig4-6-5b
    Figure 19.85: Angle θ at which the collar leaves the rod vs torque T.
  3. 3.

    The work done by the torque is W=Tθ at any instant. The system only possesses kinetic energy. So the energy of the system at any instant is E=E1+E2 where E1=12I1ω2 and E2=12m2(R˙2+R2ω2). Computing these quantities for the solution obtained above, we plot W and E vs t as shown in fig. 19.86. Clearly, W=E.

    Filename:sfig4-6-5c
    Figure 19.86: Work done by the torque and the kinetic energy of the system.
    Filename:sfig4-6-5d
    Figure 19.87: The ratio of terminal R˙/θ˙ vs torque T.
  4. 4.

    Now we take several values of T (0.1, 0.5, 1, 1.5, 2, 2.5, and 3) and carry out the numerical solutions for each T. We note the terminal values of θ,ω(=θ˙), and vR(=R˙). By plotting the terminal value of θ against T (fig. 19.85), we see that the collar leaves the rod at exactly the same θ=2.86rad for each T! But this is possible only if R˙ and θ˙ both change in the same proportion for each T. So, plot the ratio R˙/θ˙ just for the terminal values against T and find that the ratio is indeed constant (see fig. 19.87).

SAMPLE 19.18

Filename:sfig4-5-6
Figure 19.88:

Dynamics of a sliding-base pendulum. A cart of mass M slides down a frictionless inclined plane as shown in the figure. A simple pendulum of mass m and length hangs from the center-of-mass of the cart. Find the equation of motion of the pendulum.


Solution

Filename:sfig4-5-6a
Figure 19.89:

Let us measure the angular displacement of the pendulum with respect to the cart with angle θ measured anticlockwise from the normal to the inclined plane. Let s be the position of the cart along the inclined plane from some reference point. Then, the acceleration of the cart can be written as 𝒂C=s¨𝝀ˆ and the acceleration of the pendulum mass as 𝒂=𝒂C+𝒂rel=s¨𝝀ˆ+θ¨𝒆ˆθθ˙2𝒆ˆR.

The free-body diagram of the cart and the pendulum system is shown in fig. 19.89. Writing angular momentum balance of the system about point C, we get

𝑴C = 𝑯˙C
𝒆ˆR×(mgȷˆ) = 𝒆ˆR×m(s¨𝝀ˆ+θ¨𝒆ˆθθ˙2𝒆ˆR)
mglsin(θα)𝒌ˆ = ms¨cosθ𝒌ˆ+m2θ¨𝒌ˆ
 θ¨ = gsin(θα)s¨cosθ. (19.62)

To find s¨, we write the linear momentum balance for the whole system in the 𝝀ˆ (so that we do not involve the unknown normal reaction N) direction.

𝝀ˆ(Mgȷˆmgȷˆ) = 𝝀ˆ[Ms¨𝝀ˆ+m(s¨𝝀ˆ+θ¨𝒆ˆθθ˙2𝒆ˆR)]
 (M+m)gsinsinα = (M+m)s¨+mθ¨cosθmθ˙2sinθ
 s¨ = gsinαmM+m(θ¨cosθ+θ˙2sinθ)

Substituting eqn. (LABEL:eq:cartpend.s) in eqn. (19.62) and rearranging terms, we get

θ¨=gsinθ(1+mM)cosα1+mMsin2θ+mMθ˙2sinθcosθ1+mMsin2θ

Answer: θ¨=gsinθ(1+mM)cosα1+mMsin2θ+mMθ˙2sinθcosθ1+mMsin2θ

Note that if we set α=0 and let M so that the cart behaves like a fixed ground, then we recover the equation of simple pendulum, θ¨=gsinθ, from the equation of motion above. It is a good practice to carry out such simple checks wherever possible.

Remarks: We could write the equations of motion, eqn. (19.62) and eqn. (LABEL:eq:cartpend.s) in the coupled form as

[cosθmcosθM+m](θ¨s¨)=(gsin(θα)(M+m)gsinαmθ˙2sinθ)

and leave it at that, since for most computational purposes, it is enough. It is not so hard to find expressions for θ¨ and s¨ from here by solving the matrix equation, even by hand.

SAMPLE 19.19

Filename:sfig4-6-8
Figure 19.90:

Resonant capture. A slightly unbalanced motor mounted on an elastic machine part is modeled as a spring mass system with a simple pendulum of mass m and length ϵ driven by a constant torque T as shown in the figure. The spring has stiffness k and the motor has mass M. There is no friction between mass M and the horizontal surface.

  1. 1.

    Find the equation of motion of the system.

  2. 2.

    Take M=m=1kg,k=1N/m,T=15×103Nm. Solve (numerically) the equations of motion with zero initial conditions and plot x(t) and θ(t) for t=0 to 100s.

Solution

Filename:sfig4-6-8a
Figure 19.91:
  1. 1.

    The free-body diagram of the system is shown in fig. 19.91. Let the angular displacement of the eccentric mass m at some instant t be θ. At the same instant, let the displacement of the motor be x from the relaxed state of the spring. Then we can write the acceleration of the motor as x¨ıˆ and that of the eccentric mass as 𝒂P=x¨ıˆ+ϵθ¨𝒆ˆθϵθ˙2𝒆ˆR. Now, we can write the angular momentum balance for the system about point C (fixed in the stationary frame of reference but instantly coincident with the center-of-mass of motor M) as

    T𝒌ˆ = ϵ𝒆ˆR×m(x¨ıˆ+ϵθ¨𝒆ˆθϵθ˙2𝒆ˆR)
    = mϵx¨(𝒆ˆR×ıˆ)+mϵ2θ¨(𝒆ˆR×𝒆ˆθ)
    = mϵ(x¨sinθ+ϵθ¨)𝒌ˆ
     ϵθ¨sinθx¨ = Tmϵ (19.64)

    This is just one scalar equation in θ¨ and x¨. We need one more independent equation θ¨ and x¨ without involving any other unknowns. So, we write the linear momentum balance for the system in the x-direction:

    Filename:sfig4-6-8b
    Figure 19.92:
    Filename:sfig5-5-2
    Figure 19.93:
    kx = Mx¨+m(x¨ıˆ+ϵθ¨𝒆ˆθϵθ˙2𝒆ˆR)ıˆ
    = (M+m)x¨mϵθ¨sinθmϵθ˙2cosθ
     ϵsinθθ¨(M+mm)x¨ = kmxϵθ˙2cosθ. (19.65)

    Thus we have the required equations of motion. We can write eqn. (19.64) and (19.65) compactly as

    [ϵsinθϵsinθm+MM](θ¨x¨)=(Tmϵkmxϵθ˙2cosθ.)
  2. 2.

    We use the following pseudocode to solve the equations of motion. Note that we first convert the two second order ODEs into four first order ODEs by introducing new variables ω=θ˙ and u=x˙.

        Set T = 0.015, m = 1, M = 1, k = 1, e = 1
        A=[e -sin(theta); e*sin(theta) -(1+M/m)];
        b = [T/(m*e); k/m*x-e*omega^2*cos(theta)];
        solve A*acln = b for acln  % acln = accelerations
        ODEs = {omega = thetadot, u  = xdot,
                omegadot = acln(1), udot = acln(2)}
        IC   = {theta = 0, x = 0, omega = 0, u = 0}
        Solve ODEs with IC for t=0 to t=100
    

    The plots of x(t) and θ(t) obtained from the numerical solution are shown in fig. 19.92 and fig. 19.93, respectively. Note the resonance of M for the given values of the system.

SAMPLE 19.20

Filename:sfig5-5-2a
Figure 19.94:

Dynamics using a rotating and translating coordinate system. Consider the rotating wheel of Sample 18.52 which is shown here again in Figure 19.94. At the instant shown in the figure find

  1. 1.

    the linear momentum of the mass P and

  2. 2.

    the net force on the mass P.

For calculations, use a frame attached to the rod and a coordinate system in with origin at point A of the rod OA.


Solution

Filename:sfig5-5-2b
Figure 19.95: The velocity of point P is the sum of two terms: the velocity of O and the velocity of P relative to O.

We attach a frame to the rod. We choose a coordinate system xyz in this frame with its origin O at point A. We also choose the orientation of the primed coordinate system to be parallel to the fixed coordinate system xyz (see Fig. 19.95), i.e., ıˆ=ıˆ,ȷˆ=ȷˆ, and 𝐤ˆ=𝐤ˆ.

  1. 1.

    Linear momentum of P: The linear momentum of the mass P is given by

    𝑳=m𝒗P.

    Clearly, we need to calculate the velocity of point P to find 𝑳. Now,

    𝒗P = 𝒗P+𝒗rel=𝒗O+𝒗P/O𝒗P+𝒗rel.

    Note that O and P are two points on the same (imaginary) rigid body OAP. Therefore, we can find 𝒗P as follows:

    𝒗P = 𝝎×𝒓O/O𝒗O+𝝎×𝒓P/O𝒗P/O
    = ω1𝒌ˆ×(cosθıˆ+sinθȷˆ)+ω1𝒌ˆ×r(cosθıˆsinθȷˆ)
    = ω1[(+r)cosθȷˆ(r)sinθıˆ]
    = 3rad/s[2.5mcos30ȷˆ1.5msin30ıˆ]
    = (6.50ȷˆ2.25ıˆ)m/s(same as in Sample 18.52.),
    𝒗rel = 𝒗P/
    = ω2𝒌ˆ×r(cosθıˆsinθȷˆ)
    = ω2r(cosθȷˆ+sinθıˆ)
    = (2.16ȷˆ+1.25ıˆ)m/s
    = (2.16ȷˆ+1.25ıˆ)m/s.

    Therefore,

    𝒗P = 𝒗P+𝒗rel
    = (4.34ȷˆ3.50ıˆ)m/sand
    𝑳 = m𝒗P
    = 0.5kg(4.34ȷˆ3.50ıˆ)m/s
    = (1.75ıˆ+2.17ȷˆ)kgm/s.

    Answer: 𝑳=(1.75ıˆ+2.17ȷˆ)kgm/s

  2. 2.

    Net force on P: From the

    𝑭=m𝒂

    for the mass P we get 𝑭=m𝒂P. Thus to find the net force 𝑭 we need to find 𝒂P.

    The calculation of 𝒂P is the same as in Sample 18.52 except that 𝒂P is now calculated from

    𝒂P=𝒂O+𝒂P/O

    where

    Filename:sfig7-4-2
    Figure 19.96:
    𝒂O = 𝝎×(𝝎×𝒓O/O)
    = ω12𝒓O/O
    = ω12(cosθıˆ+sinθȷˆ)
    = (3rad/s)22m(cos30ıˆ+sin30ȷˆ)
    = (15.59ıˆ+9.00ȷˆ)m/s2,
    𝒂P/O = 𝝎×(𝝎×𝒓P/O)
    = ω12𝒓P/O
    = ω12r(cosθıˆsinθȷˆ)
    = (3rad/s)20.5m(cos30ıˆsin30ȷˆ)
    = (3.90ıˆ2.25ȷˆ)m/s2.

    Thus,

    𝒂P=(19.49ıˆ+6.75ȷˆ)m/s2

    which, of course, is the same as calculated in Sample 18.52. The other two terms, 𝒂cor and 𝒂rel, are exactly the same as in Sample 18.52. Therefore, we get the same value for 𝒂P by adding the three terms:

    𝒂P=(17.83ıˆ+3.63ȷˆ)m/s2.

    The net force on P is

    𝑭 = m𝒂P
    = 0.5kg(17.83ıˆ3.63ȷˆ)m/s2
    = (8.92ıˆ+1.81ȷˆ)N.

    Answer: 𝑭=(8.92ıˆ+1.81ȷˆ)N

SAMPLE 19.21

Filename:pfigure-blue-99-1
Figure 19.97:

Inverse dynamics of a four bar mechanism. A four bar mechanism ABCD consists of three uniform bars AB, BC, and CD of length 1,2,3, and mass m1,m2,m3, respectively. The mechanism is driven by a torque T applied at A such that bar AB rotates at constant angular speed. Write equations to find the torque T at some instant t.


Solution This is an inverse dynamics problem, that is, we are given the motion and we are supposed to find the forces (torque T in this case) that cause that motion. We are given that rod AB rotates at constant angular speed, say θ˙. From kinematics, we can find out angular velocities and angular accelerations of the other two bars as well as the accelerations of center-of-mass of each rod. Then we can write the momentum balance equations and compute the forces and moments required to generate this motion. So, in contrast to what we usually do, let us do the kinematics first. Please see Sample 18.5 on page 18.5. We found the angular velocities, β˙ (eqn. (18.110)) and ϕ˙ (eqn. (18.109)), of rods BC and CD, respectively, in terms of θ˙. We can rewrite those equations as

[2sinβ3sinϕ2cosβ3cosϕ](β˙ϕ˙)=θ˙(1sinθ1cosθ) (19.66)

We wrote this equation in matrix form to make it easier for us to find the angular accelerations which we do by simply differentiating this equation once:

[2cosββ˙3cosϕϕ˙2sinββ˙3sinϕϕ˙](β˙ϕ˙) + [2sinβ3sinϕ2cosβ3cosϕ](β¨ϕ¨)
= 1(θ¨sinθ+θ˙2cosθθ¨cosθθ˙2sinθ)

Rearranging terms we get

[2sinβ3sinϕ2cosβ3cosϕ](β¨ϕ¨) = [2cosβ3cosϕ2sinβ3sinϕ](β˙2ϕ˙2) (19.82)
+1[sinθcosθcosθsinθ](θ¨θ˙2)

Thus, we can find the angular accelerations of BC and CD, β¨ and ϕ¨, because the quantities on the right hand side are known (θ¨(=0) and θ˙ are given, and β˙ and ϕ˙ are determined by eqn. (19.66)). Now, we can find the accelerations of center-of-mass of each rod as follows.

𝒂G1 = 12θ˙2𝝀ˆ1 (19.83)
𝒂G2 = 𝒂B+𝒂G2/B=1θ˙2𝝀ˆ122β˙2𝝀ˆ2+2β¨𝒏ˆ2 (19.84)
𝒂G3 = 32ϕ˙2𝝀ˆ3 (19.85)

We are now ready to write momentum balance equations. Since we are only interested in finding the torque T, we should try to write equations involving minimum number of unknown forces. So, we draw free-body diagrams of the whole mechanism, of part BCD, and of bar CD alone; and write angular momentum balance equations about appropriate points so that we involve only the unknown torque T and the unknown reaction 𝑹D at D. Thus, we will have only three scalar unknowns T, RDx and RDy (since 𝑹D=RDxıˆ+RDyȷˆ). So, we will need only three independent equations.

Consider the free-body diagram of the whole mechanism. We can write angular momentum balance about point A for the whole mechanism as

Filename:pfigure-blue-47-2
Figure 19.98:
T𝒌ˆ+𝒓D/A×𝑹D=𝑯˙A=𝑯˙1/A+𝑯˙2/A+𝑯˙3/A (19.86)

where

𝑯˙1/A = I1θ¨𝒌ˆ+𝒓G1×m1𝒂G1
𝑯˙2/A = I2β¨𝒌ˆ+𝒓G2×m2𝒂G2=I2β¨𝒌ˆ+(𝒓B+𝒓G2/B)×m2𝒂G2
𝑯˙3/A = I3ϕ¨𝒌ˆ+𝒓G3×m3𝒂G3=I3ϕ¨𝒌ˆ+(𝒓D+𝒓G3/D)×m3𝒂G3.
Filename:pfigure-blue-89-1
Figure 19.99:

Similarly, the angular momentum balance about point B for BCD gives

𝒓D/B×𝑹D=𝑯˙2/B+𝑯˙3/B (19.87)

where

𝑯˙2/B = I2β¨𝒌ˆ+𝒓G2/B×m2𝒂G2
𝑯˙3/B = I3ϕ¨𝒌ˆ+𝒓G3/B×m3𝒂G3

and angular momentum balance of bar CD about point C gives

Filename:pfigure-blue-35-2
Figure 19.100:
𝒓D/C×𝑹D=𝑯˙3/C=I3ϕ¨𝒌ˆ+𝒓G3/C×m3𝒂G3. (19.88)

Note that we can easily write the position vectors in terms of 1,2,3 and the unit vectors (𝝀ˆ1,𝒏ˆ1), (𝝀ˆ2,𝒏ˆ2) and (𝝀ˆ3,𝒏ˆ3) where

𝝀ˆ1=cosθıˆ+sinθȷˆ, 𝒏ˆ1=sinθıˆ+cosθȷˆ
𝝀ˆ2=cosβıˆ+sinβȷˆ, 𝒏ˆ2=sinβıˆ+cosβȷˆ
𝝀ˆ1=cosϕıˆ+sinϕȷˆ, 𝒏ˆ1=sinϕıˆ+cosϕȷˆ.

We can put all the three angular momentum balance equations, (19.86), (19.87), and (19.88), in one matrix equation by dotting both sides of the equations with 𝒌ˆ and assembling them as follows.

[1040𝒌ˆ(2𝝀ˆ23𝝀ˆ3)×ıˆ𝒌ˆ(2𝝀ˆ23𝝀ˆ3)×ȷˆ0𝒌ˆ(3𝝀ˆ3×ıˆ)𝒌ˆ(3𝝀ˆ3×ȷˆ)](TRDxRDy)=(H˙123/AH˙23/BH˙3/C) (19.98)
(19.99)

where H˙123/A=𝒌ˆ(𝑯˙1/A+𝑯˙2/A+𝑯˙3/A), H˙23/B=𝒌ˆ(𝑯˙2/B+𝑯˙3/B), and H˙3/C=𝒌ˆ𝑯˙3/C.

Note that we know the 𝑯˙’s on the right hand side and the matrix on the left side can be evaluated for any given (θ,β,ϕ). Thus we can solve for T,RDx, and RDy.

SAMPLE 19.22

Numerical solution of the inverse dynamics problem. Consider Sample 19.95 again. Using numerical solutions on a computer, find and plot torque T against θ for one complete cycle of the drive arm AB. Take m1=m2=m3=1kg and 1=400mm,2=4002mm,3=8002mm, and 4=1200mm.


Solution

Since we have to plot T against θ for one complete revolution, we need to find angular velocities, angular accelerations and center-of-mass accelerations for several values of θ and then solve for T for each of those θ’s. We can do this several ways. One way would be to first solve kinematic equations to find θ(t),β(t), and ϕ(t) at discrete times over one complete cycle and then compute all other quantities at each (θ(ti),β(ti),ϕ(ti)) where ti represents a discrete time. So, let us follow this method step by step with pseudocodes. Here, we assume that we have vector functions called dot and cross that compute the dot product and the cross product of two vectors that are given as input arguments.

Step-1: solve for angular positions.

Specify the given geometry

L1=0.4, L2=0.4*sqrt(2), L3=0.8*sqrt(2), L4=1.2

and use the pseudocode of Sample 18.5 to find θ(ti),β(ti),ϕ(ti) for, say, 100 values of ti between 0 and 1 sec. Now, for each triad of (θ(ti),β(ti),ϕ(ti)), follow all the steps below.

Step-2: solve for angular velocities.

Since θ˙=2πrad/s is given, we only need to solve for β˙ and ϕ˙. We use eqn. (18.110) and eqn. (18.109) to compute β˙ and ϕ˙ as follows (or modify the pseudocode of Sample 18.5 to save β˙ and ϕ˙ along with the values for β and ϕ.

  define thdot=thetadot, bdot=betadot, pdot=phidot
  thdot = 2*pi   % this is given
  set th = theta(ti), b = beta(ti), p = phi(ti)
  set unit vectors
    l1=[cos(th) sin(th) 0]’, n1=[-sin(th) cos(th) 0]’
    l2=[cos(b) sin(b) 0]’, n2=[-sin(b) cos(b) 0]’
    l3=[cos(p) sin(p) 0]’, n3=[-sin(p) cos(p) 0]’
  bdot = -(L1/L2)*(cross(n1,l3)/cross(n2,l3))*thdot
  pdot = (L1/L3)*(cross(n1,l2)/cross(n3,l2))*thdot
Step-3: solve for angular accelerations.

Now that we have (θ,β,ϕ) and the corresponding values of (θ˙,β˙,ϕ˙), we can use eqn. (19.82) to calculate β¨ and ϕ¨ (we are given θ¨=0).

  define thddot=thetaddot, bddot=betaddot,
         pddot=phiddot
  thddot = 0   % this is given
  B = [-L2*sin(b) L3*sin(p); -L2*cos(b) L3*cos(p)]
  C = L1*[sin(th) cos(th); cos(th) -sin(th)]
  D = [-cos(b) cos(p); sin(b) -sin(p)]
  c = [thddot thdot^2]’,  d = [L2*bdot^2 L3*pdot^2]’
  assume w = [bddot pddot]’
  solve B*w = C*c + D*d for w

So, now we know θ¨,β¨,ϕ¨ also. We are now ready to compute 𝑯˙’s required for dynamic calculations.

Step-4: set up equations and solve for unknown forces.

We need to set up and solve eqn. (19.99). Note that we need to compute several quantities for this equation but the vector computations are more or less straightforward.

   % set mass and inertia properties
   m1 = 1, m2 = 1, m3 = 1
   I1 = m1*L1^2/12, I2 = m2*L2^2/12, I3 = m3*L3^2/12
   % set fixed unit vectors
   i = [1 0 0]’, j = [0 1 0]’, k = [0 0 1]’
   % compute position vectors
   rA = [0;0;0],  rB = rA+L1*l1
   rC = rB+L2*l2, rD = L4*l4
   rG1 = L1/2*l1
   rG2 = rB+L2/2*l2
   rG3 = rD+L3/2*l3
   % compute center-of-mass accelerations
   aG1 = 0.5*L1*(tddot*n1-tdot^2*l1)
   aG2 = 2*aG1+0.5*L2*(bddot*n2-bdot^2*l2)    % aB = 2*aG1
   aG3 = 0.5*L3*(pddot*n3-pdot^2*l3)
Filename:pfigure4-3Dpend
Figure 19.101: Torque T as a function of θ over one complete cycle of motion (θ(0)=π/2,θ(1)=5π/2).
   % compute Hdot_cms

   Hdot_cm1 = I1*tddot*uk
   Hdot_cm2 = I2*bddot*uk
   Hdot_cm3 = I3*pddot*uk
   % compute Hdots
   Hdot_123_A = Hdot_cm1 + cross(rG1, m1*aG1)
                 + Hdot_cm2 + cross(rG2, m2*aG2)
                   + Hdot_cm3 + cross(rG3, m3*aG3)
   Hdot_23_B = Hdot_cm2 + cross(rG2-rB, m2*aG2)
                 + Hdot_cm3 + cross(rG3-rB, m3*aG3)
   Hdot_3_C = Hdot_cm3 + cross(rG3-rC, m3*aG3)
   % set up the linear eqns for torque and RD
   b = [dot(Hdot_123_A,k) dot(Hdot_23_B,k) dot(Hdot_3_C,k)]
   A = [1  dot(k,cross(rD,i))     dot(k,cross(rD,j))
        0  dot(k,cross(rD-rB,i))  dot(k,cross(rD-rB,j))
Ψ        0  dot(k,cross(rD-rC,i))  dot(k,cross(rD-rC,j))]
   % let forces = [T  RDx  RDy]’
   solve A*forces = b for forces
Step-5, repeat calculations.

Now repeat Step-2 – Step-4 for each triad (θ,β,ϕ) obtained in Step-1 and save the corresponding values of T in a vector. Finally,

   plot T vs theta
Filename:pfigure-s94h6p3
Figure 19.102: Torque T as a function of time over one complete cycle of motion.

The plot thus obtained is shown in fig. 19.101. We can also plot T vs time (as shown in fig. 19.102), and, of course, expect to see the same graph of T since θ is just a linear function of t. Note that the area under the graph of T over one complete cycle must equal zero since the net impulse must be zero over one cycle.

Problems for 19.3 Multi-degree-of-freedom 2-D mechanisms

19.3.1  Particle on a springy leash. A particle with mass m slides on a rigid horizontal frictionless plane. It is held by a string which is in turn connected to a linear elastic spring with constant k. The string length is such that the spring is relaxed when the mass is on top of the hole in the plane. The position of the particle is 𝒓=xıˆ+yȷˆ. For each of the statements below, state the circumstances in which the statement is true (assuming the particle stays on the plane). Justify your answer with convincing explanation and/or calculation.

  1. (a)

    The force of the plane on the particle is mg𝒌ˆ.

  2. (b)

    x¨+kmx=0

  3. (c)

    y¨+kmy=0

  4. (d)

    r¨+kmr=0,where r=|𝒓|

  5. (e)

    r= constant

  6. (f)

    θ˙=constant

  7. (g)

    r2θ˙= constant.

  8. (h)

    m(x˙2+y˙2)+kr2= constant

  9. (i)

    The trajectory is a straight line segment.

  10. (j)

    The trajectory is a circle.

  11. (k)

    The trajectory is not a closed curve.

Filename:pfigure-s94h6p4
Figure 19.103

Particle on a springy leash.

19.3.2  “Yo-yo” mechanism of satellite de-spinning. A satellite, modeled as a uniform disk, is “de-spun” by the following mechanism. Before launch two equal length long strings are attached to the satellite at diametrically opposite points and then wound around the satellite with the same sense of rotation. At the end of the strings are placed 2 equal masses. At the start of de-spinning the two masses are released from their position wrapped against the satellite. Find the motion of the satellite as a function of time for the two cases: (a) the strings are wrapped in the same direction as the initial spin, and (b) the strings are wrapped in the opposite direction as the initial spin.

19.3.3  A particle with mass m is held by two long springs each with stiffness k so that the springs are relaxed when the mass is at the origin. Assume the motion is planar. Assume that the particle displacement is much smaller than the lengths of the springs.

  1. (a)

    Write the equations of motion in cartesian components. Answer: k(xıˆ+yȷˆ)=m(x¨ıˆ+y¨ȷˆ).

  2. (b)

    Write the equations of motion in polar coordinates. Answer: kreˆr=m(r¨rθ˙2)eˆr+m(rθ¨+2r˙θ˙)eˆt.

  3. (c)

    Express the conservation of angular momentum in cartesian coordinates. Answer: ddt(xy˙yx˙)=0.

  4. (d)

    Express the conservation of angular momentum in polar coordinates. Answer: ddt(r2θ˙)=0.

  5. (e)

    Show that (a) implies (c) and (b) implies (d) even if you didn’t note them a-priori. Answer: dot equation in (a) with (xȷˆyıˆ) to get (xy¨yx¨)=0 which can be rewritten as (c); dot equation in (b) with 𝒆ˆt to get (rθ¨+2r˙θ˙)=0 which can be rewritten as (d).

  6. (f)

    Express the conservation of energy in cartesian coordinates. Answer: 12m(x˙2+y˙2)+12k(x2+y2)=const.

  7. (g)

    Express the conservation of energy in polar coordinates. Answer: 12m(r˙2+(rθ˙)2)+12kr2=const.

  8. (h)

    Show that (a) implies (f) and (b) implies (g) even if you didn’t note them a-priori. Answer: dot equation in (a) with 𝒗=x˙ıˆ+y˙ȷˆ to get m(x¨x˙+y¨y˙)+k(x˙x+y˙y)=0 which can be rewritten as (f).

  9. (i)

    Find the general motion by solving the equations in (a). Describe all possible paths of the mass. Answer: x=A1sin(ωt+B1), y=A2sin(ωt+B2), where ω=km. The general motion is an ellipse.

  10. (j)

    Can the mass move back and forth on a line which is not the x or y axis? Answer: Yes. Consider B1=B2=0, A1=1, and A2=2.

Filename:pg84-3
Figure 19.104

19.3.4  Cart and pendulum A mass mB=6kg hangs by two strings from a cart with mass mC=12kg. Before string BC is cut string AB is horizontal. The length of string AB is r=1m. At time t=0 all masses are stationary and the string BC is cleanly and quietly cut. After some unknown time tvert the string AB is vertical.

  1. (a)

    What is the net displacement of the cart xC at t=tvert? Answer: 0.33mıˆ

  2. (b)

    What is the velocity of the cart vC at t=tvert? Answer: vC=1.82m/sıˆ.

  3. (c)

    What is the tension in the string at t=tvert? Answer: T240N.

  4. (d)

    (Optional) What is t=tvert? [You will either have to leave your answer in the form of an integral you cannot evaluate analytically or you will have to get part of your solution from a computer.]

Filename:pfigure-blue-36-1
Figure 19.105

19.3.5  A dumbbell slides on a floor. Two point masses m at A and B are connected by a massless rigid rod with length . Mass B slides on a frictionless floor so that it only moves horizontally. Assume this dumbbell is released from rest in the configuration shown. [Hint: What is the acceleration of A relative to B?]

  1. (a)

    Find the acceleration of point B just after the dumbbell is released.

  2. (b)

    Find the velocity of point A just before it hits the floor

Filename:pfigure-s95q8
Figure 19.106

19.3.6   After a winning goal one second before the clock ran out a psychologically stunned hockey player (modeled as a uniform rod) stands nearly vertical, stationary and rigid. The player’s perfectly slippery skates start to pop out from under her as she falls. Her height is , her mass m, her tip from the vertical θ, and the gravitational constant is g.

  1. (a)

    What is the path of her center of mass as she falls? (Show clearly with equations, sketches or words.)

  2. (b)

    What is her angular velocity just before she hits the ice, a millisecond before she sticks out her hands and breaks her fall (first assume her skates remain in contact the whole time and then check the assumption)?

  3. (c)

    Find a differential equation that only involves θ, its time derivatives, m, g, and . (This equation could be solved to find θ as a function of time. It is a non-linear equation and you are not being asked to solve it numerically or otherwise.)

19.3.7  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:pfigure4-rpi
Figure 19.107

19.3.8  A model for a yo-yo consists of a thin disk of mass M and radius R and a light drum of radius r, rigidly attached to the disk, around which a light inextensible cable is wound. Assuming that the cable unravels without slipping on the drum, determine the acceleration aG of the center of mass.

Filename:pfigure-blue-86-1
Figure 19.108

19.3.9  A uniform rod with mass mR pivots without friction about point A in the xy-plane. A collar with mass mC slides without friction on the rod after the string connecting it to point A is cut. There is no gravity. Before the string is cut, the rod has angular velocity ω1.

  1. (a)

    What is the speed of the collar after it flies off the end of the rod? Use the following values for the constants and initial conditions: mR=1kg, mC=3kg, a=1m, =3m, and ω1=1rad/s

  2. (b)

    Consider the special case mR=0. Sketch (approximately) the path of the motion of the collar from the time the string is cut until some time after it leaves the end of the rod.

Filename:pfigure-blue-88-1
Figure 19.109

19.3.10  Assume the rod in the figure for problem 18.18 has polar moment of inertia Iozz. Assume it is free to rotate. The bead is free to slide on the rod. Assume that at t=0 the angular velocity of the rod is 1rad/s, that the radius of the bead is one meter and that the radial velocity of the bead, dR/dt, is zero.

  1. (a)

    Draw separate free-body diagrams of the bead and rod.

  2. (b)

    Write equations of motion for the system. Answer:

    R¨Rθ˙2 = 0
    (Izz+mR2)θ¨+2mRR˙θ˙ = 0
  3. (c)

    Use the equations of motion to show that angular momentum is conserved. Answer: The second equation in part (b) can be rewritten in the form ddt[Izz+mR2)θ˙]=0. The quantity inside the derivative is angular momentum; thus, it is conserved and equal to a constant, say, (H/o)0, which can be found in terms of the initial conditions.

  4. (d)

    Find one equation of motion for the system using: (1) the equations of motion for the bead and rod and (2) conservation of angular momentum.

    Answer: R¨R[(Ho)02(IzzO+mR2)2]=0.

  5. (e)

    Write an expression for conservation of energy. Let the initial total energy of the system be, say, E0.

    Answer: E0=12mR˙2+12ω(H/o)0.

  6. (f)

    As t goes to infinity does the bead’s distance go to infinity? Its speed? The angular velocity of the turntable? The net angle of twist of the turntable? Answer: The bead’s distance goes to infinity and its speed approaches a constant. The turntable’s angular velocity goes to zero and its net angle of twist goes to a constant.

Filename:pfigure4-2-rp3
Figure 19.110

Coupled motion of bead and rod and turntable.

19.3.11  A primitive gun rides on a cart (mass M) and carries a cannon ball (of mass m) on a platform at a height H above the ground. The cannon ball is dropped through a frictionless tube shaped like a quarter circle of radius R.

  1. (a)

    If the system starts from rest, compute the horizontal speed (relative to the ground) that the cannon ball has as it leaves the bottom of the tube. Also find the cart’s speed at the same instant.

  2. (b)

    Compute the velocity of the cannon ball relative to the cart.

  3. (c)

    If two balls are dropped simultaneously through the tube, what speed does the cart have when the balls reach the bottom? Is the same final speed also achieved if one ball is allowed to depart the system entirely before the second ball is released? Why/why not?

Filename:pfigure-s94q6p1
Figure 19.111

19.3.12   Numerically simulate the coupled system in problem 19.3. Use the simulation to show that the net angle of the turntable or rod is finite.

  1. (a)

    Write the equations of motion for the system from part (b) in problem 19.3 as a set of first order differential equations. Answer:

    θ˙ = ω
    ω˙ = 2mRvωIzz+mR2
    R˙ = v
    v˙ = ω2R
  2. (b)

    Numerically integrate the equations of motion.

19.3.13  Two frictionless prisms of similar right triangular sections are placed on a frictionless horizontal plane. The top prism weighs W and the lower one nW. The prisms are held in the initial position shown and then released, so that the upper prism slides down along the lower one until it just touches the horizontal plane. The center of mass of a triangle is located at one-third of its height from the base. Compute the velocities of the two prisms at the moment just before the upper one reaches the bottom. Answer: v2x=2g(ab)tanϕ(1+1n+(n+1n)2tan2ϕ), v2y=v2xn+1ntanϕ, v1x=v2xn, v1y=0, where 2 refers to the top wedge and 1 refers to the lower wedge.

Filename:pfigure-spinningbrick
Figure 19.112

19.3.14  Mass slides on an accelerating cart. 2D. A cart is driven by a powerful motor to move along the 30 sloped ramp according to the formula: d=d0+vot+a0t2/2 where d0, v0, and a0 are given constants. The cart is held from tipping over. The cart itself has a 30 sloped upper surface on which rests a mass (given mass m). The surface on which the mass rests is frictionless. Initially the mass is at rest with regard to the cart.

  1. (a)

    What is the force of the cart on the mass? [in terms of g,d0,v0,t,a0,m,g,ıˆ, and ȷˆ.] Answer: F=m(g+a0)[34ıˆ+34ȷˆ].

  2. (b)

    For what values of d0, v0, t and a0 is the acceleration of the mass exactly vertical (i.e., in the ȷˆ direction)? Answer: For a0=g, the acceleration of the mass is exactly vertical; d0, v0, and t could be anything.

Filename:pfigure-s94f1p2
Figure 19.113

19.3.15  A thin rod AB of mass WAB=10lbm and length LAB=2ft is pinned to a cart C of mass WC=10lbm, the latter of which is free to move along a frictionless horizontal surface, as shown in the figure. The system is released from rest with the rod in the horizontal position.

  1. (a)

    Determine the angular speed of the rod as it passes through the vertical position (at some later time). Answer: angular speed = 8.8rad/s.

  2. (b)

    Determine the displacement x of the cart at the same instant. Answer: displacement = 0.5ft.

  3. (c)

    After the rod passes through vertical, it is momentarily horizontal but on the left side of the cart. How far has the cart moved when this configuration is reached?

Filename:s97f2
Figure 19.114

19.3.16  As shown in the figure, a block of mass m rolls without friction on a rigid surface and is at position x (measured from a fixed point). Attached to the block is a uniform rod of length which pivots about one end which is at the center of mass of the block. The rod and block have equal mass. The rod makes an angle θ with the vertical. Use the numbers below for the values of the constants and variables at the time of interest:

= 1m
m = 2kg
θ = π/2
dθ/dt = 1rad/s
d2θ/dt2 = 2rad/s2
x = 1m
dx/dt = 2m/s
d2x/dt2 = 3m/s2
  1. (a)

    What is the kinetic energy of the system?

  2. (b)

    What is the linear momentum of the system (momentum is a vector)?

Filename:pg92-2
Figure 19.115

19.3.17  A pendulum of length hangs from a cart. The pendulum is massless except for a point mass of mass mp at the end. The cart rolls without friction and has mass mc. The cart is initially stationary and the pendulum is released from rest at an angle θ. What is the acceleration of the cart just after the mass is released? [Hints: 𝐚P=𝒂C+𝒂P/C. The answer is 𝒂C=(g/3)𝐢 in the special case when mp=mc and θ=π/4. ]

Filename:summer95p2-3
Figure 19.116

19.3.18  Due to the application of some unknown force F, the base of the pendulum A is accelerating with 𝒂A=aAı. There is a frictionless hinge at A. The angle of the pendulum θ with the x axis and its rate of change θ˙ are assumed to be known. The length of the massless pendulum rod is . The mass of the pendulum bob M. There is no gravity. What is θ¨? (Answer in terms of aA, , M, θ and θ˙.)

Filename:p-f96-p3-3
Figure 19.117

19.3.19  Pumping a Swing Can a swing be pumped by moving the support point up and down?

For simplicity, neglect gravity and consider the problem of swinging a rock in circles on a string. Let the rock be mass m attached to a string of fixed length . Can you speed it up by moving your hand up and down? How? Can you make a quantitative prediction? Let xS be a function of time that you can specify to try to make the mass swing progressively faster.

Filename:s97p3-3
Figure 19.118

19.3.20  Using free-body diagrams and appropriate momentum balance equations, find differential equations that govern the angle θ and the vertical deflection y of the system shown. Be clear about your datum for y. Your equations should be in terms of θ,y and their time-derivatives, as well as M,m,, and g.

Filename:pfigure-s94h8p1
Figure 19.119

19.3.21  The two blocks shown are released from rest at t=0. There is no friction and the cable is initially taut. (a) What is the tension in the cable immediately after release? (Use any reasonable value for the gravitational constant). (b) What is the tension after 5 s?

Filename:sfig4-7-DH1
Figure 19.120

19.3.22  Carts A and B are free to move along a frictionless horizontal surface, and bob C is connected to cart B by a massless, inextensible cord of length , as shown in the figure. Cart A moves to the left at a constant speed v0=1m/s and makes a perfectly plastic collision (e=0) with cart B which, together with bob C, is stationary prior to impact. Find the maximum vertical position of bob C, hmax, after impact. The masses of the carts and pendulum bob are mA=mB=mC=10kg.

Filename:sfig4-7-DH2
Figure 19.121

19.3.23  A double pendulum is made of two uniform rigid rods, each of length . The first rod is massless. Find equations of motion for the second rod. Define any variables you use in your solution.

Filename:sfig4-7-DH3
Figure 19.122

19.3.24  A model of walking involves two straight legs. During the part of the motion when one foot is on the ground, the system looks like the picture in the figure, confined to motion in the xy plane. Write two equations from which one could find θ¨1, and θ¨2 given θ1,θ2,θ˙1,θ˙2 and all mass and length quantities.

[Hint: 𝑴/A=𝑯˙/Awhole system𝑴/B=𝑯˙/Bfor bar BC]
Filename:sfig4-7-1
Figure 19.123

Advanced problems in 2D motion

19.3.25  A pendulum is hanging from a moving support in the xy-plane. The support moves in a known way given by 𝒓(t)=r(t)ıˆ. For the following cases, find a differential equation whose solution would determine θ(t), measured clockwise from vertical; then find an expression for d2θ/dt2 in terms of θ and r(t):

  1. (a)

    with no gravity and assuming the pendulum rod is massless with a point mass of mass m at the end,

  2. (b)

    as above but with gravity,

  3. (c)

    assuming the pendulum is a uniform rod of mass m and length .

19.3.26  Robotics problem: balancing a broom stick by sideways motion. Try to balance a broom stick by moving your hand horizontally. Model your hand contact with the broom as a hinge. You can model the broom as a uniform stick or as a point mass at the end of a stick — your choice.

  1. (a)

    Equation of motion. Given the acceleration of your hand (horizontal only), the current tip, and the rate of tip of the broom, find the angular acceleration of the broom. Answer: For point mass: ϕ¨=(g/L)sinϕ(ahand/L)cosϕ.

  2. (b)

    Control? Can you find a hand acceleration in terms of the tip and the tip rate that will make the broom balance upright? Answer: There are many correct solutions. Test your solution with a computer simulation. Show your result with appropriate plots.

Filename:sfig4-7-1a
Figure 19.124

19.3.27  Balancing the broom again: vertical shaking works too. You can balance a broom by holding it at the bottom and applying appropriate torques, as in problem 16.147 or by moving your hand back and forth in an appropriate manner, as in problem 19.124. In this problem, you will try to balance the broom differently. The lesson to be learned here is more subtle, and you should probably just wonder at it rather than try to understand it in detail.

In the previous balancing schemes, you used knowledge of the state of the broom (θ and θ˙) to determine what corrective action to apply. Now balance the broom by moving it in a way that ignores what the broom is doing. In the language of robotics, what you have been doing before is ‘closed-loop feedback’ control. The new strategy, which is simultaneously more simple minded and more subtle, is an ‘open loop’ control. Imagine that your hand connection to the broom is a hinge.

  1. (a)

    Picture and model. Assume your hand oscillates sinusoidally up and down with some frequency and some amplitude. The broom is instantaneously at some angle from vertical. Draw a picture which defines all the variables you will use. Use any mass distribution that you like.

  2. (b)

    FBD. Draw a FBD of the broom.

  3. (c)

    Momentum balance. Write the equation of angular momentum balance about the point instantaneously coinciding with the hinge.

  4. (d)

    Kinematics. Use any geometry and kinematics that you need to evaluate the terms in the angular momentum balance equation in terms of the tip angle and its time derivatives and other known quantities (take the vertical motion of your hand as ‘given’). [Hint: There are many approaches to this problem.]

  5. (e)

    Equation of motion. Using the angular momentum balance equation, write a governing differential equation for the tip angle. Answer: ϕ¨[y¨+gL]sinϕ=0, where y(t) is the vertical displacement of your hand and ϕ(t) is the angle of the broom from the vertical.

  6. (f)

    Simulation. Taking the hand motion as given, simulate on the computer the system you have found.

  7. (g)

    Stability? Can you find an amplitude and frequency of shaking so that the broom stays upright if started from a near upright position? You probably cannot find linear equations to solve that will give you a control strategy. So, this problem might best be solved by guessing on the computer. Successful strategies require the hand acceleration to be quite a bit bigger than g, the gravitational constant. The stability obtained is like the stability of an undamped uninverted pendulum — oscillations persist. You improve the stability a little by including a little friction in the hinge.

  8. (h)

    Dinner table experiment for nerdy eaters. If you put a table knife on a table and put your finger down on the tip of the blade, you can see that this experiment might work. Rapidly shake your hand back and forth, keeping the knife from sliding out from your finger but with the knife sliding rapidly on the table. Note that the knife aligns with the direction of shaking (use scratch-resistant surface). The knife is different from the broom in some important ways: there is no gravity trying to ‘unalign’ it and there is friction between the knife and table that is much more significant than the broom interaction with the air. Nonetheless, the experiment should convince you of the plausibility of the balancing mechanism. Because of the large accelerations required, you cannot do this experiment with a broom.

Filename:sfig4-7-1b
Figure 19.125

The sketch of the knife on the table goes with part (h).

19.3.28  Double pendulum. The double pendulum shown is made up of two uniform bars, each of length and mass m. The pendulum is released from rest at ϕ1=0 and ϕ2=π/2. Just after release what are the values of ϕ¨1 and ϕ¨2? Answer in terms of other quantities. Answer: ϕ¨1=0,ϕ¨2=3g/2.

Filename:sfig4-7-2
Figure 19.126

19.3.29  A rocker. A standing dummy is modeled as having massless rigid circular feet of radius R rigidly attached to their uniform rigid body of length L and mass m. The feet do not slip on the floor.

  1. (a)

    Given the tip angle ϕ, the tip rate ϕ˙ and the values of the various parameters (m, R, L, g) find ϕ¨. [You may assume ϕ and ϕ˙ are small.]

    Answer: ϕ¨+[3(RL2)gL2]ϕ=0.

  2. (b)

    Using the result of (a) or any other clear reasoning find the conditions on the parameters (m, R, L, g) that make vertical passive dynamic standing stable. [Stable means that if the person is slightly perturbed from vertically up that their resulting motion will be such that they remain nearly vertically up for all future time.]

    Answer: This first-order approximation show that stable standing requires R>L2. R=L2 is neutrally stable (like a wheel). (Someone questions this answer, so don’t trust it too much.)

Filename:sfig4-7-2a
Figure 19.127

19.3.30  Consider a rigid spoked wheel with no rim. Assume that when it rolls a spoke hits the ground and doesn’t bounce. The body just swings around the contact point until the next spoke hits the ground. The uniform spokes have length R. Assume that the mass of the wheel is m, and that the polar moment of inertia about its center is I (use I=mR2/2 if you want to get a better sense of the solution). Assume that just before collision number n, the angular velocity of the wheel is ωn, the kinetic energy is Tn, the potential energy (you must clearly define your datum) is Un. Just after collision n the angular velocity of the wheel is ωn+. The Kinetic Energy is TN+, the potential energy (you must clearly define your datum) is Un+. The wheel has k spokes (pick k=4 if you have trouble with abstraction). This problem is not easy. It can be answered at a variety of levels. The deeper you get into it the more you will learn.

  1. (a)

    What is the relation between ωn and Tn?

  2. (b)

    What is the relation between ωn and ωn+?

  3. (c)

    Assume ‘rolling’ on level ground. What is the relation between ωn+ and ωn+1+ (the angular velocities just after two successive collisions?

  4. (d)

    Assume rolling down hill at slope θ. What is the relation between ωn+ and ωn+1?

  5. (e)

    Can it be true that ωn+ = ω(n+1)+? About how fast is the wheel going in this situation?

  6. (f)

    As the number of spokes m goes to infinity, in what senses does this wheel become like an ordinary wheel?

Filename:sfig4-7-2b
Figure 19.128