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:
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.
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.
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.
Use appropriate balance equations: linear momentum, angular momentum, or power balance equations.
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.
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.
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.
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.
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.
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?
“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.
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.
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:
The mass of the particle 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.
In general abstract vector form: .
In cartesian coordinates: .
In polar coordinates: .
In path coordinates:
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 , rotation rate and angular acceleration :
where, to simplify the notation, all motions are relative to and positions relative to 0 unless explicitly indicated by a or .
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.
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).
General abstract form. The equation of linear momentum balance is or, if there is no force, , which means that . So is a constant. We can call this constant . So after some time the particle is where it was at , say, , plus its velocity times time. That is:
| (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
Because the net force is zero and the net mass is not negligible,
These equations imply that , , and are all constants, let’s call them , , . So , , and are given by
We can put these components into their place in vector form to get:
| (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:
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).
A tedious calculation will show that these equations are solved by the following functions of time:
| (19.3) | |||||
where , , , , , and are constants. Note that, though eqn. (19.1) looks different than eqn. (19.2), there are still 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, is expressed in path coordinates as
That is,
So the speed 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 we have that is written as
| (19.4) | |||||
Now, using the rotating base vectors and , we have that and
so eqn. (19.4) can be rewritten as
which in turn can be broken into components and written as:
| (19.5) | |||||
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
then the solution turns out to be
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 frame††margin: 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. .
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.
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 , 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 ( is an unknown scalar).
For some purposes, especially general problems like this where no specific path is given, the most appropriate configuration variable is , 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
So
Now we can write linear momentum balance
| (19.6) | |||||
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:
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.
Example: The hard brachistochrone problem.
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 from the origin and vertical position (drawn accurately in fig. 19.3):
| (19.7) |
Starting at the origin this curve is close to but gets a bit higher (bigger ) because for a given value of , is greater than . The curve terminates at vertical tangents at where . (see fig. 19.3). To solve the puzzle this curve is scaled (by choosing a value of ) 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 as driven by gravity? (see fig. 19.3)
Let’s use as our configuration variable. The power balance equation is:
This, remarkably, is the simple harmonic oscillator equation with general solution
Thus the period of oscillation ( such that ) is
which is independent of the amplitude of oscillation
.
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 which we also didn’t care about. We could have got the same equation of motion by writing and eliminated by dotting both sides with a convenient vector orthogonal to , 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 because the formula doesn’t tell you the coordinates
of the points. You could find them by integrating
numerically or with calculus tricks. But it turns out (see below) that the
curve with is described by the parametric equations
with .
.
Note that the osculating circle of this cycloid at its lowest point has radius , 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.
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 and only 1.7% longer if the amplitude is , but enough to annoy clock designers. A bead sliding frictionlessly on the path 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 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 and another from wheel 2 rolling under line a distance below. Both wheels have radius . Imagine that the cycloids and are drawn by wheels always arranged with vertically aligned rolling contact points and and with points and initially aligned vertically a distance apart at and .
The two cycloids are thus the same shape but are displaced with one being below and sideways from the other.
Because both wheels have rolled the same distance () they have rotated the same amount and is as far forward of as is behind. Similarly is as far above as is below. So the line is bisected by the point .
Because is the diameter of a circle with on the perimeter, angle is a right angle. Because material point on the wheel has zero velocity the velocity of (and thus the tangent to the curve) is orthogonal to . Thus the line is tangent to the upper cycloid.
The rolling of wheel 2 instantaneously rotating about makes the tangent to the lower cycloid orthogonal to , the condition for the motion of to be from the wrapping of an inextensible line around the curve . This shows that cycloid is generated by the wrapping of a line anchored at 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 of must be and the arc length of is so the length is . By the similarity of the two right triangles that share the length of :
which shows that the upper cycloid is the curve if , where 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
Consider a pair of collars hinged together as a point mass 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 and , and two sets of rotating base vectors: and . These rotating base vectors can be written in terms of the s, and in the standard manner. Assume we know and in this configuration. First find and by thinking of the velocity of the collar two different ways:
| (19.8) | |||||
Having found and 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) | |||||
We use the results from eqn. (19.8) for and to evaluate the right hand sides of the expressions for and 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
| (19.10) | |||||
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 and with letters and substituted that into the expressions for and and left those expressions intact while substituting for the acceleration we would have large expressions for the force components and . 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 and 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 and . 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.
Example: Inverted pendulum with a vibrating base
Assume that the base of an inverted point-mass pendulum of length is vibrated according to (see fig. 19.6)
The point P thus has acceleration
Now apply linear momentum balance as
which you write as
depending on whether you are analytically or numerically inclined. This is a second order non-linear ordinary differential equation. If or 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 and that the pendulum does not fall down! Just shaking the base keeps the pendulum up ( 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
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
A bead on a straight wire. A straight wire is hung between points A and B in the 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
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 and components that are easily derived by differentiating the geometric constraint equation with respect to . Thus,
Answer:
SAMPLE 19.2
A particle sliding on a parabolic path. A particle slides on a parabolic trough given by where 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 on its path.
Solution The geometric constraint on the path of the particle is already given, . Differentiating the path constraint with respect to time, we get the constraint on velocity and acceleration components.
Now, at a point , we can write the velocity and acceleration of the particle as
Answer:
SAMPLE 19.3
Circular motion of a particle. A particle is constrained to move on a frictionless circular path of radius with constant angular speed . There is no gravity. Find the equation of motion of the particle in the -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 (in polar coordinates). This constraint gives us and . Then the acceleration of the particle (in polar coordinates), reduces to , (of course).
The free-body diagram of the particle shows that there is only one force acting on the particle, the normal reaction of the path acting in the direction. Therefore, the linear momentum balance gives,
But, to write the equation of motion in the -direction, we need to write the linear momentum balance in the -direction. We can write using mixed basis vectors as . Dotting this equation with , we get
or, , which is the equation of simple harmonic motion in . You can easily show that the motion in the -direction is also simple harmonic ().
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 and 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.
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
We can easily eliminate the constraint force from this equation by dotting this equation with , which gives
But, from the geometric constraint , we have . Therefore,
Since is constant, we integrate the equation of motion easily to find since . And, since , we have
Answer:
SAMPLE 19.5
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 acts normal to the path. Let and be unit vectors tangential and normal to the path, respectively. Then the linear momentum balance gives
To eliminate the unknown constraint force from this equation, we can take a dot product of this equation with . However, we must first find . Now is the unit tangent vector. So, we can find it by finding a tangent vector to the path (remember gradient of a function ?) 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 . Why not use the velocity vector ? The velocity vector is always tangential to the path. Furthermore, we know that from the geometric constraint (), and . Therefore, . Now, dotting the LMB equation with we get,
Rearranging the terms above, we get the required equation of motion:
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 , i.e., the coordinate of the position of the bead. But, you can always find the coordinate since .
Answer:
Comment: Note that if we consider and to be very small so that we can ignore the terms completely, and take , then the equation of motion becomes
which is the equation of simple harmonic motion with frequency . 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
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 where and . 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 . Find the magnitude of the force on the pin at .
Solution Let denote the net force on the pin. Then from the linear momentum balance
where is the acceleration of the pin. Therefore, to find the force at we need to find the acceleration at that position.
From the given figure, we assume that the pin is in the groove at . 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
We are given that and the radial position of the pin . Therefore, ††margin: Note that is a function of and is a function of time, therefore is a function of time. Although we are interested in finding and at , we cannot first substitute in the expression for and then take its time derivatives (which will be zero).
Substituting these expressions in the acceleration formula and then substituting the numerical values at , (remember, must be in radians!), we get
Therefore the net force on the pin is
and the magnitude of the net force is
Answer:
SAMPLE 19.7
A puck sliding on a frictional rotating table. A horizontal turntable rotates with constant angular speed . A puck of mass 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
A fixed reference frame with cartesian coordinates
A rotating reference frame with cartesian coordinates
Solution
Equation of motion using a fixed reference frame: The puck has two DOF on the turntable. So, we will need two configuration variables, say and , 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 be the position of the puck at some instant , so that its velocity is and acceleration is .
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 and that . 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
So, we need to find . Now 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 . Clearly, goes in circles with constant speed, so that its velocity is
Therefore, the relative velocity, is
Now the linear momentum balance for the puck in the plane gives
These are coupled nonlinear ODEs that represent the equations of motion of the puck.
Answer:
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 .
Equation of motion using a rotating reference frame: Now we derive the equations of motion using a rotating reference frame, , with coordinate axes, fixed to the rotating turntable. Let the position of the puck in the rotating frame be . Note that the velocity of the puck in the rotating reference frame is , and acceleration is .
Now, from the linear momentum balance () for the puck, we get
where we have used the three term acceleration formula for . Here,
Note that the point , 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 . Substituting these terms in the LMB equation above, we have
Dotting this equation with and , respectively, we get
These are the required equations of motion for the puck in the rotating frame.
Answer:
Once we find a solution and of these equations, we can find the
solution in the fixed frame by transforming () to () through
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
A collar sliding on a frictional rod. A collar of mass slides on a massless rigid rod OA of length . The rod rotates counterclockwise with a constant angular speed . The coefficient of friction between the rod and the collar is . At time , the bar is horizontal, the collar is at and has radial speed towards the pivot O. Ignore gravity.
How does the position of the collar change with time (i.e., what is the equation of motion of the collar)?
Plot the path of the collar starting from till the collar shoots off the end of the bar.
How long does it take for the collar to leave the bar?
Solution
A Free-Body Diagram of the collar at a general position 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 acting normal to the rod and the force of friction acting along the rod. Thus linear momentum balance for the collar is:
| (19.11) |
Note that because the rod is rotating at a constant rate. Dotting both sides of eqn. (19.11) with and we get
Eliminating from the last two equations we get
Since is constant, the above equation is of the form
| (19.12) |
where 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 . Why? Look at the left hand side of eqn. (19.11). The net reaction force is acting in the direction . If we dot the reaction force with a vector normal to its direction, we get rid of the reaction force. A vector normal to is .
which is the same equation as eqn. (19.12).
Solution of equation (19.12): The characteristic equation††margin: If we assume a solution of the form 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
Therefore, the solution of Eqn. (19.12) is
Substituting the given initial conditions: and , we get††margin: From the general solution and the given initial conditions, we have: Solving these two equations simultaneously, we get , and . Substituting the values of and , we get
| (19.13) |
Answer:
To draw the path of the collar we need both and . Because constant,
Now we can take various values of from to, say, , and calculate values of and . Plotting all these values of and , however, does not give us the correct path after the collar reaches the end of the bar, at . So we need to find the time when . Equation (19.13) is a nonlinear algebraic equation for which we can solve iteratively on a computer, or with some patience, even on a calculator. One way to find would be to simply plot and find the intersection with (see fig. 19.21) and read the corresponding value of . Following either of these methods we find that here. Now, we can plot the path of the collar by computing and from to 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)
The plot produced thus is shown in Fig. 19.22.
The time computed above was
By plugging this value in the expression for (Eqn. (19.13) we get, indeed,
Answer:
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:
SAMPLE 19.9
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 slides on the rod. Assume that at time the collar is released from a rest position with respect to the rod at a distance from the axis of rotation. There is no gravity.
Find the equation of motion for the collar (a differential equation for the position of the collar).
How does the distance of the collar from the vertical axis change with time?
For , show that the solution obtained in (ii) above is the same as that obtained in Sample 19.15 for .
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 be the cylindrical coordinates of the collar at any general time (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
where is a unit normal to the rod. Thus, if is a unit vector along the rod, then
| (19.14) |
where
Equation of Motion: We now write the linear momentum balance for the collar:
where
Dotting both sides of the equation by we get
This expression is an equation of motion of the collar. However, it is only a single equation in terms of derivatives of two variables and . Now, from geometry
Substituting this relationship in the equation of motion we get
Noting that a constant, we may write the above equation, with some trigonometric simplifications, as
| (19.15) |
which is an equation of motion of the collar in terms of its distance from the vertical axis.
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.
where and are arbitrary constants to be determined from the initial conditions. Substituting the given initial conditions: and we get
Therefore, the solution may be written as
| (19.16) |
SAMPLE 19.10
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 gm. The shape of the helix is given: and where . At the instant when , the angular speed and the angular acceleration of the bead are and . Find
the net normal force on the bead and
the radius of the osculating circle at the instant given.
††margin:
Solution
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):
Thus we need to find the normal acceleration of the bead. We can write the position of the bead as
Substituting the given numerical values for and in the above expressions for and we get the velocity and the acceleration of the bead at the instant of interest:
In path coordinates,
where
| (19.17) | |||||
| (19.18) | |||||
| (19.19) |
Therefore,
Hence, the net normal force on the bead
Answer:
For calculating the radius of the osculating circle, we note that
Answer:
19.1.1
A bead slides on a frictionless circular hoop . The mass of the bead ,
mass of hoop and
radius of hoop .
Neglect gravity. The center of the circular hoop is the origin of a fixed (Newtonian) coordinate system . The hoop is on the plane. The hoop is kept from moving by little angels who let the bead slide by unimpeded. At , the speed of the bead is , it is traveling counter-clockwise (looking down the -axis), and it is on the -axis. There are no other external forces applied.
At what is the bead’s kinetic energy ?
At what is the bead’s linear momentum ?
At what is the bead’s angular momentum about the origin?
At what is the bead’s acceleration?
At what is the radius of the osculating circle of the bead’s path ?
At what is the force of the hoop on the bead?
At what is the net force of the angel’s hands on the hoop?
At what is the -component of the bead’s linear momentum?
19.1.2 A warehouse operator wants to move a crate of weight from a 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.
What are the kinetic and potential energies(pick a suitable datum) of the crate at the elevated platform, point A?
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?
Using conservation of energy, calculate the maximum speed attained by the crate at point B, assuming it starts moving from rest at point A.
If the crate slides a distance , say, on the horizontal conveyor, what is the energy lost to sliding in terms of , , and ?
Calculate the value of the coefficient of friction, , such that the crate comes to rest at point C.
19.1.3 On a wintry evening, a student of mass starts down the steepest street in town, Steep Street, which is of height and slope . At the top of the slope she starts to slide (coefficient of dynamic friction ).
What is her initial kinetic and potential energy?
What is the energy lost to friction as she slides down the hill?
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).
If, upon reaching the bottom of the street, she collides with another student of mass and they embrace, what is their instantaneous mutual velocity just after the embrace?
Assuming that friction still acts on the level flats on the bottom, how much time will it take before they come to rest?
19.1.4 Reconsider the system of blocks in problem 2.3, this time with equal mass, . 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 . The coefficient of friction between the lower block and the floor is . It is known that .
Draw free-body diagrams of the blocks together and separately.
Write the equations of linear momentum balance for each block.
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 ?
19.1.5 An initially motionless roller-coaster car of mass 20 is given a horizontal impulse at position A, causing it to move along the track, as shown below.
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 so that the car “just makes it” over the hill at 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.
19.1.6 Masses and move on the frictionless varied terrain shown. Initially, has speed and is at rest. The two masses collide on the flat section. The coefficient of restitution in the collision is . Find the speed of at the top of the first hill of elevation . Does make it over the second hill of elevation ?
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 with constant . A dead bug with mass m is standing on the platform. Her feet have no glue on them. There is gravity.
Draw a free-body diagram of the bug.
What is the acceleration of the platform and, hence, the bug?
Write the equation of linear momentum balance for the bug.
What condition must be met so that the bug does not bounce the platform?
What is the maximum spin rate that can be used if the bug is not to bounce off the platform?
19.1.8 A bug of mass walks straight-forwards with speed and rate of change of speed 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 axis. Assume that the distance the bug is from the origin, , the angle the stick makes with the axis, , and its rate of change, , are known at the instant of interest. Ignore gravity. Answer the following questions in terms of , , , , , and .
What is ?
What is the force exerted by the rod on the support?
What is the acceleration of the bug?
19.1.9 A new kind of gun. Assume is a constant for the rod in the figure. Assume the mass is free to slide. At , the rod is aligned with the -axis and the bead is one foot from the origin and has no radial velocity .
Find a differential equation for . Answer: .
Turn this equation into a differential equation for . Answer: .
How far will the bead have moved after one revolution of the rod? How far after two? Answer: , .
What is the speed of the bead after one revolution of the rod (use )? Answer:
How much kinetic energy does the bead have after one revolution and where did it come from? Answer: .
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 has radius with rate of change . The angle is zero and is a constant. Neglect gravity.
What is at this instant? Give your answer in terms of any or all of , , , , , , , , and . Answer: .
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: .
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 with . The rigid rod, on which the bead slides, rotates at a constant rate . Find in terms of , and .
19.1.12 Slippery bead on straight rotating stick. A long stick with mass 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 . Initially the mass is at a distance from the hinge point on the stick and has no radial velocity (). The initial angle of the stick is , measured counterclockwise from the positive axis.
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 ?
What is the path of the bead in the -plane? (Draw an accurate picture showing about one half of one revolution.)
How long should the stick be if the bead is to fly in the negative -direction when it gets to the end of the stick?
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 versus . Note, you have been working with 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 , the angular velocity of the rod/turntable is , that the radius of the bead is one meter, and that the radial velocity of the bead, , is zero. The bead is free to slide on the rod. Where is the bead at ?
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 with at , 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 gm is constrained to slide in a straight frictionless slot in a disk which is spinning counterclockwise at constant rate . At time , the slot is parallel to the -axis and the bead is in the center of the disk moving out (in the plus direction) at a rate . After a net rotation of one and one eighth () 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.
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 . The turntable speed is controlled by a strong stiff motor.
Assume for all time. What are possible motions of the bead?
Assume is a constant. What are possible motions of the bead? Notice there are two cases depending on the value of . What is going on here?
Mass in a slot.
19.1.17 A small bead with mass slides without friction on a rigid rod which rotates about the axis with constant (maintained by a stiff motor not shown in the figure). The bead is also attached to a spring with constant , 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 from the center of the rod and then released with an initial . If needed, you may assume that .
Derive the equation of motion for the position of the bead .
How is the motion affected by large versus small values of ?
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 , , , .
Write an expression for the Coriolis acceleration. Give an example of a situation in which this acceleration is important and explain why it arises.
19.1.18 A small ball of mass may slide in a slender tube of length . and of mass . The tube rotates freely about a vertical axis passing through its center C. (Hint: Treat the tube as a slender rod.)
If the angular velocity of the tube is as the ball passes through with speed relative to the tube, , calculate the angular velocity of the tube just before the ball leaves the tube. Answer: , where the minus sign ‘-’ means ‘just before leaving’.
Calculate the angular velocity of the tube just after the ball leaves the tube. Answer: , where the plus sign ‘+’ means ‘just after leaving’.
If the speed of the ball as it passes is , determine the transverse and radial components of velocity of the ball as it leaves the tube. Answer: , .
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 ? Answer: torque=
19.1.19 Toy train car on a turn-around. The toy train’s speedometer reads a constant when, heading west, it passes due north of point . 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 . What is the force of the turntable on the train? (Don’t worry about the -component of the force.) Answer: .
19.1.20 Due to forces not shown, the cart moves to the right with constant acceleration . The ball has mass . At time , the string is cut. Find
the tension in string before cutting, Answer: .
the absolute acceleration of the mass at the instant of cutting, Answer:
the tension in string at the instant of cutting. Answer: .
19.1.21 The forked arm mechanism pushes the bead of mass along a frictionless hyperbolic spiral track given by . The arm rotates about its pivot point at with constant angular acceleration driven by a motor (not shown). The arm starts from rest at .
Determine the radial and transverse components of the acceleration of the bead after have elapsed from the start of its motion.
Determine the magnitude of the net force on the mass at the same instant in time.
19.1.22 A rod is on the palm of your hand at point . Its length is . Its mass is assumed to be concentrated at its end at . Assume that you know and at the instant of interest. Also assume that your hand is accelerating both vertically and horizontally with . Coordinates and directions are as marked in the figure.
Draw a Free-Body Diagram of the rod.
Assume that the hand is stationary. Solve for in terms of , , , , and .
If the hand is not stationary but has been determined somehow, find the vertical force of the hand on the rod in terms of , , , , , , , and .
Rod on a Moving Support(2-D).
19.1.23 Balancing a broom. Assume the hand is accelerating to the right with acceleration . What is the force of the hand on the broom in terms of , and ? (You may not have any or in your answer.)
A one-degree-of-freedom mechanism is a collection of parts linked together so that they can move in only one way††margin: 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.
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
Body is connected to body by a frictionless hinge at C (see fig. 19.45). The force on body at C from is and the force on from body at C is . The power of the interaction force on body is . This power contributes to the increase in the kinetic energy of . The power of the interaction force on is . 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,
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
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 . On the other hand for frictionless sliding there 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.
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 . Further assume that the machine has no motion when . The variable 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 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 ,
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 and its rate of change . Because at any configuration the velocity of every point in the system is proportional to we can write the kinetic energy as:
where is a function that one can determine by calculating the machine’s total kinetic energy in terms of and and then factoring out of the resulting expression.
Now, if we accept the equation of mechanical energy conservation we have
| (19.20) | ||||
The cancellation of above lacks mathematical rigor, but doesn’t cause problems††margin: The cancellation of the factor from equation 19.20 depends on being other than zero. While moving, is not zero. Strictly we cannot cancel the term from the equation at the instants when . 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 and 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
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:
Similarly power balance could have been used to get the same result, looking at just the mass
as before.
Example: Pendulum
Consider a rigid body with mass and moment of inertia 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.
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
Assume all three bars in the structure shown in fig. 19.48are of equal length and have mass 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 (when all centers of mass are at hinge height) we have
which is the same governing equation as for a point-mass pendulum with length . 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
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 according to
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 times the speed of the pedal relative to the frame :
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 and that the wheels have mass and 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:
| (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.
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 . For rolling contact so, for one wheel:
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
| (19.22) | |||||
| (19.23) | |||||
| (19.24) | |||||
| (19.25) | |||||
| (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. 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 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.
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 . To keep our notation simple, let’s assume that we have defined so that at this minimum. Physically this means that 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:
| (19.27) |
Applying this result to equation( 19.20) we get:
| (19.28) |
We now write in terms of its Taylor series. We have
| (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 where 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 . Thus the middle term involving is negligibly smaller than other terms. Similarly, using the Taylor series for , the last term is well approximated by , where is a constant which we will call . Now we have for the equation of motion:
| (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 and the effective mass is the coefficient of in the expansion for the kinetic energy . The displacement of any part of the system from equilibrium will thus be given by
| (19.31) |
with , and and 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.
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 rolling without slip in an cylindrical ‘ideal’ valley of radius .
For this problem we can calculate and in terms of . Briefly,
So we can derive the equation of motion using the fact of constant total energy.
Now, assuming small angles, so , we get
| (19.32) | |||||
| (19.33) |
This is, naturally, our old friend the harmonic oscillator equation. The period is a funny combination of terms. If it looks like a point mass pendulum with length , more than . 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 and 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
A plate pendulum. A rectangular plate of mass 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
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.
The free-body diagram of the plate is shown in fig. 19.53. Note that the link forces and act along the links because massless links are two-force bodies. Let be the coordinates of the center-of-mass. Then the linear momentum balance for the plate gives
Now, we can eliminate both the unknown link forces from this equation by dotting the equation with , a unit vector normal to the links. Then, we have
| (19.34) |
Now we need to find a relationship between and , and and , so that we can write and in terms of our configuration variable and its derivatives. From fig. 19.52, we have
Substituting these expressions for and in eqn. (19.34) , we get
Answer:
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 and coordinates of the center-of-mass, we have
that is, the center-of-mass follows a circle of radius centered at . 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 going in a circle centered at and driven by gravity. That is the simple pendulum.
SAMPLE 19.12
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 and radius 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 . 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, , here to find the equation of motion of the wheel. Since and , we have
| (19.35) |
Now, we need to find and express it using the configuration variable and its derivatives. There are several ways we could find . Vectorially, we could write, where . Dotting both sides of this equation with and we can find in terms of and in terms of and . But, for a change, let us use geometry here.
See fig. 19.55. From triangle ABO, we have
| Now, | ||||
Substituting this expression for in power balance eqn. (19.35), we get
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 .
Answer:
SAMPLE 19.13
Instantaneous dynamics of slider crank. A uniform rigid rod AB of mass and length has one of its ends pinned to the rim of a disk of radius . 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 . At the instant shown, end B of the rod is directly above the center of the disk making to be .
Find all the forces acting on the rod.
Is there a value of which makes end A of the rod lift off the horizontal surface when ?
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.
Forces on the rod: The free-body diagram of the rod is shown in Fig. 19.57. The pin at B exerts two forces and while the surface in contact at A exerts only a normal force because there is no friction. Now, we can write the momentum balance equations for the rod. The linear momentum balance () for the rod gives
| (19.36) |
The angular momentum balance about the center-of-mass G () of the rod gives
| (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 — (counts as two unknowns), and . Therefore, we need more equations. We have already used the momentum balance equations, hence, the extra equations have to come from kinematics.
Dotting both sides of the equation with we get
Also,
Dotting both sides of this equation by we get
| (19.38) |
Now, we can find the acceleration of the center-of-mass:
Substituting for from eqn. (19.38) and for above, we obtain
Substituting this expression for in eqn. (19.36) and dotting both sides by and then by we get
| (19.39) |
From eqn. (19.37)
Lift off of end A: End A of the rod loses contact with the ground when normal force becomes zero. From the expression for from above, this condition is satisfied when
SAMPLE 19.14
A bar sliding on a sliding wedge. A bar AB of mass and length is hinged at end A and rests on a wedge of mass 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
momentum balance
energy (or power) balance.
Solution
(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., . Now, the angular momentum balance for the bar about point A gives,
| (19.41) | |||||
where the last line follows from the fact that , the unit normal , and so, and . We now need to eliminate the unknown normal reaction from the above equation. Since the wedge is constrained to move only horizontally, we can write the linear momentum balance for the wedge as
| (19.42) |
Thus, we have found in terms of that we can use in eqn. (19.41) to get rid of . But, we now need to express 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 AC denote the horizontal position of the wedge. Then, from the law of sines, we have , so that
| (19.43) | |||||
| (19.44) |
Now, substituting for in eqn. (19.41) from eqn. (19.42), using the expression for from above, and dotting the resulting equation with , we get
| (19.45) | |||||
Answer:
(b) Power balance: Now we derive the equation of motion for the bar using power balance . 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 . The wedge moves with horizontal speed , therefore, its kinetic energy is . Thus, . The only force that contributes to power is the force of gravity on the rod because is non-zero. The sliding contact at B is frictionless and hence the net power due to the contact force there is zero. Now, . So, . Thus the power balance for the system gives
| (19.46) |
We can simplify this equation further. Note that, and . So that we can write the above equation as
where is a constant of integration to be determined from initial conditions. For example, when the bar begins to slide from rest, we have, at , and . At that instant, if , then . So, we can write
Now, replacing in this equation with the expression we obtained in eqn. (19.43), we have
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 , in which case, this reduction to the first order equation does not save much work.
Answer:
19.2.1 A conservative vibratory system has the following equation of conservation of energy.
where is a constant.
Obtain the differential equation of motion of this system by differentiating this energy equation with respect to .
Determine the circular frequency of small oscillations of the system in part (a). HINT: (Let and .
19.2.2 A motor at O turns at rate whose rate of change is . At the end of a stick connected to this motor is a frictionless hinge attached to a second massless stick. Both sticks have length . At the end of the second stick is a mass . For the configuration shown, what is ? Answer in terms of , , , , , and . Ignore gravity.
19.2.3
In problem 18.3, find
the angular momentum about point , Answer: .
the rate of change of angular momentum of the disk about point , Answer: .
the angular momentum about point , Answer: .
the rate of change of angular momentum of the disk about point . Answer: .
Assume the rod is massless and the disk has mass .
19.2.4 Robot arm, 2-D . The robot arm AB is rotating about point A with and . Meanwhile the forearm BC is rotating at a constant angular speed with respect to AB of . Gravity cannot be neglected. At the instant shown, find the net force acting on the object P which has mass . Answer: .
19.2.5 A crude model for a column, shown in the figure, consists of two identical rods of mass and length with hinge connections, a linear torsional spring of stiffness attached to the center hinge (the spring is relaxed when ), and a load applied at the top end.
Obtain the exact nonlinear equation of motion.
Obtain the squared natural frequency for small motions .
Check the stability of the straight equilibrium state for all via minimum potential energy. How do these results compare with those from part(b)?
19.2.6 A thin uniform rod of mass rests against a frictionless wall and on a frictionless floor. There is gravity.
Draw a free-body diagram of the rod.
The rod is released from rest at . Write the equation of motion of the rod. Answer: .
Using the equation of motion, find the initial angular acceleration, , and the acceleration of the center of mass, , of the rod. Answer: From the diagram, we see ,
.
Find the reactions on the rod at points and . Answer: , .
Find the acceleration of point . Answer: .
When , find and the acceleration of point . Answer: At , , .
19.2.7 Bar leans on a crooked wall. A uniform bar leans on a wall and floor and is let go from rest. Gravity pulls it down.
Draw a Free-Body Diagram of the bar.
Kinematics: find the velocity and acceleration of point B in terms of the velocity and acceleration of point A
Using equations of motion, find the acceleration of point A.
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.)
19.2.8 Two blocks, each with mass , 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?
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 (with length ) when the slider crank is in the position with (piston is at maximum extent)? Assume the crankshaft has constant angular velocity , that the connecting rod is massless, that the cylinder walls are frictionless, that there is no gas pressure in the cylinder, that the piston has mass and the crank has radius .
19.2.10 In problem 18.84, find the force on the wheel at point due to the rod.
19.2.11 In problem 18.4, find the force on the rod at point .
19.2.12 Slider-Crank mechanism. The slider crank mechanism shown is used to push a block P. Arm AB and BC are each 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. poundals to the left. (Note: poundal .)
19.2.13 The large masses at and were supported by the light triangular plate ABC, whose corners follow the guide. B enters the curved guide with velocity . Neglecting gravity, find the vertical reactions (in the direction) at A and B. (Hint: for the rigid planar body , find and in terms of and , assuming initially.)
19.2.14 An idealized model for a car comprises a rigid chassis of mass and four identical rigid disks (wheels) of mass and radius , as shown in the figure. Initially in motion with speed , the car is momentarily brought to rest by compressing an initially uncompressed spring of stiffness . Assume no frictional losses.
Assuming no wheel slip, determine the compression of the spring required to stop the car. Answer: .
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!
Repeat part (a) assuming now that the ground is perfectly frictionless from points A to B shown in the figure. Answer: .
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 ) 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.
In your first misconceived experiment you set up a 70 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?
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).
19.2.16 Which way does the bike accelerate? A bicycle with all frictionless bearings is standing still on level ground. A horizontal force 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.
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
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 and its derivatives ( and ). 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:
Draw FBDs of each object,
Pick configuration variables,
Write linear and angular momentum balance equations
Solve the equations for variables of interest (usually forces and second derivatives of the configuration variables).
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”.
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 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.
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
.
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 rolls without friction on ideal massless wheels at A and B (see fig. 19.76). Block 2 with mass 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
| (19.47) | |||||
Looking at the free-body diagram of mass 2 note that there are no unknown forces in the direction, so
| (19.48) |
Eqns. 19.47 and 19.48 are a system of two equations in the two unknowns and
which can be solved for and by hand or on the computer. Finding and is then easy because both and 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
Similarly for the upper block:
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: and . 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 and 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 or . So all of the reactions and as well as the accelerations and 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 and coordinates of and and the rotations of the two bodies:
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
A robot arm has two links. There are motors that apply known torque at the shoulder (reacted by the base and a torque 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:
| (19.53) | |||||
| (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 and we implicitly mean for the whole system. And and 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.
| and | ||||
| and | ||||
| and | ||||
| , | ||||
| , | ||||
| and |
| and |
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 and . One can go through the algebra and solve for them explicitly, but the expressions are quite complex, even when simplified. At given values of and however these are just two linear equations in two unknowns.
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
Dynamics of sliding wedges. A wedge shaped body of mass sits on a frictionless ground. Another wedge shaped body of mass 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 with respect to ).
Solution
The free-body diagrams of the two wedges are shown in fig. 19.79. Note that the friction force is since the wedges are sliding with respect to each other (if they were not sliding already then the friction force is an unknown force ). Let the absolute acceleration of be . Then, the absolute acceleration of is . Now, we can write the linear momentum balance for and as follows.
| (19.55) | |||||
| (19.56) |
where and . Here, we have 4 independent scalar equations (from the two 2-D vector equations) in four unknowns , and . Thus, we can certainly solve for them. We are, however, only interested in . So, we should try to find the answer with fewer calculations. Dotting eqn. (19.55) with , we have
| (19.57) |
So, to find , we need and . Dotting eqn. (19.55) with , we have
| (19.58) |
and dotting eqn. (19.56) with , we have
| (19.59) |
Solving eqn. (19.58) and (19.59) simultaneously, and using new variables (for convenience) , and , we get
Substituting these expression in eqn. (19.57), we get
Answer:
Note that when there is no friction , the expression for reduces to
and if we let (i.e., represents fixed ramp) so that , then which is the acceleration of a point mass down a frictionless ramp of slope .
SAMPLE 19.16
Dynamics of a new gun. A new gun consists of a uniform rod AB of mass and a small collar C of mass that slides freely on the rod. A motor at A rotates the rod with constant torque .
Find the equations of motion of the collar.
Show that if then the equations of motion imply conservation of angular momentum about point A.
Solution
Let us denote the configuration of the collar with , 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 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 , and letting be the moment of inertia of the rod about A, we have,
Dotting this equation with , we have
| (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 direction, i.e., , so that we do not have to care about the unknown normal reaction . So, we have,
| (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 and we need to solve them numerically.
Answer:
Now we set in our equations of motion. Note that the equation for is independent of . The equation for becomes
But the last expression is simply for the system. Thus we have which implies that 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 , and . Carry out numerical solutions for the following cases.
Let the system start from rest at and . Find the solution from to . Plot and (in polar coordinates).
Find the solution till the collar leaves the rod. What is the speed of the collar at this instant?
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.
Vary torque and carry out solutions for several values of . Find the terminal value of (when the collar leaves the rod) for each . Justify your observation about by plotting as a function of .
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 , so that we have,
Given the values of all constants, we only need to specify the initial conditions for , and for solving these equations numerically.
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
We do not know a priori the value of at which the collar leaves the rod. So, we have to carry out the solution for some assumed which gives us so that we know the collar has gone past the end of the rod. We then plot , including the unreal value of , and find the time at which , either by zooming into the graph or by interpolation (although, there are various sophisticated algorithms to find this ). Following the method of zooming into the graph
(see fig. 19.84)
we find the terminal value of to be 1.147. We carry out the numerical solution again from to and find that , and , so that This is the terminal speed of the collar.
The work done by the torque is at any instant. The system only possesses kinetic energy. So the energy of the system at any instant is where and . Computing these quantities for the solution obtained above, we plot and vs as shown in fig. 19.86. Clearly, .
Now we take several values of (0.1, 0.5, 1, 1.5, 2, 2.5, and 3) and carry out the numerical solutions for each . We note the terminal values of , and . By plotting the terminal value of against (fig. 19.85), we see that the collar leaves the rod at exactly the same for each ! But this is possible only if and both change in the same proportion for each . So, plot the ratio just for the terminal values against and find that the ratio is indeed constant (see fig. 19.87).
SAMPLE 19.18
Dynamics of a sliding-base pendulum. A cart of mass slides down a frictionless inclined plane as shown in the figure. A simple pendulum of mass and length hangs from the center-of-mass of the cart. Find the equation of motion of the pendulum.
Solution
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 be the position of the cart along the inclined plane from some reference point. Then, the acceleration of the cart can be written as and the acceleration of the pendulum mass as .
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
| (19.62) |
To find , we write the linear momentum balance for the whole system in the (so that we do not involve the unknown normal reaction ) direction.
Substituting eqn. (LABEL:eq:cartpend.s) in eqn. (19.62) and rearranging terms, we get
Answer:
Note that if we set and let so that the cart behaves like a fixed ground, then we recover the equation of simple pendulum, , 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
and leave it at that, since for most computational purposes, it is enough. It is not so hard to find expressions for and from here by solving the matrix equation, even by hand.
SAMPLE 19.19
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 and length driven by a constant torque as shown in the figure. The spring has stiffness and the motor has mass . There is no friction between mass and the horizontal surface.
Find the equation of motion of the system.
Take . Solve (numerically) the equations of motion with zero initial conditions and plot and for to 100.
Solution
The free-body diagram of the system is shown in fig. 19.91. Let the angular displacement of the eccentric mass at some instant be . At the same instant, let the displacement of the motor be from the relaxed state of the spring. Then we can write the acceleration of the motor as and that of the eccentric mass as . 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
| (19.64) |
This is just one scalar equation in and . We need one more independent equation and without involving any other unknowns. So, we write the linear momentum balance for the system in the -direction:
| (19.65) |
Thus we have the required equations of motion. We can write eqn. (19.64) and (19.65) compactly as
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 .
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 and obtained from the numerical solution are shown in fig. 19.92 and fig. 19.93, respectively. Note the resonance of for the given values of the system.
SAMPLE 19.20
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
the linear momentum of the mass P and
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
We attach a frame to the rod. We choose a coordinate system 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 (see Fig. 19.95), i.e.,
Linear momentum of P: The linear momentum of the mass P is given by
Clearly, we need to calculate the velocity of point P to find . Now,
Note that O′ and P′ are two points on the same (imaginary) rigid body OAP′. Therefore, we can find as follows:
Therefore,
Answer:
Net force on P: From the
for the mass P we get . Thus to find the net force we need to find .
The calculation of is the same as in Sample 18.52 except that is now calculated from
where
Thus,
which, of course, is the same as calculated in Sample 18.52. The other two terms, and , are exactly the same as in Sample 18.52. Therefore, we get the same value for by adding the three terms:
The net force on P is
Answer:
SAMPLE 19.21
Inverse dynamics of a four bar mechanism. A four bar mechanism ABCD consists of three uniform bars AB, BC, and CD of length , and mass , respectively. The mechanism is driven by a torque applied at A such that bar AB rotates at constant angular speed. Write equations to find the torque at some instant .
Solution This is an inverse dynamics problem, that is, we are given the motion and we are supposed to find the forces (torque 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
| (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:
Rearranging terms we get
| (19.82) | |||||
Thus, we can find the angular accelerations of BC and CD, and , because the quantities on the right hand side are known ( 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.
| (19.83) | |||||
| (19.84) | |||||
| (19.85) |
We are now ready to write momentum balance equations. Since we are only interested in finding the torque , 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 and the unknown reaction at D. Thus, we will have only three scalar unknowns , and (since ). 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
| (19.86) |
where
Similarly, the angular momentum balance about point B for BCD gives
| (19.87) |
where
and angular momentum balance of bar CD about point C gives
| (19.88) |
Note that we can easily write the position vectors in terms of and the unit vectors , and where
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.
| (19.98) | |||
| (19.99) |
where , , and .
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 , and .
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 against for one complete cycle of the drive arm AB. Take and , and .
Solution
Since we have to plot 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 for each of those ’s. We can do this several ways. One way would be to first solve kinematic equations to find , and at discrete times over one complete cycle and then compute all other quantities at each where 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.
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 for, say, 100 values of between 0 and 1 sec. Now,
for each triad of , follow all the steps
below.
Since 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
Now that we have and the corresponding values of , we can use eqn. (19.82) to calculate and (we are given ).
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.
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)
% 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
Now repeat Step-2 – Step-4 for each triad obtained in Step-1 and save the corresponding values of in a vector. Finally,
plot T vs theta
The plot thus obtained is shown in fig. 19.101. We can also plot vs time (as shown in fig. 19.102), and, of course, expect to see the same graph of since is just a linear function of . Note that the area under the graph of over one complete cycle must equal zero since the net impulse must be zero over one cycle.
19.3.1 Particle on a springy leash. A particle with mass 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 . 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 . 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.
The force of the plane on the particle is .
constant
constant
constant.
constant
The trajectory is a straight line segment.
The trajectory is a circle.
The trajectory is not a closed curve.
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 is held by two long springs each with stiffness 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.
Write the equations of motion in cartesian components. Answer: .
Write the equations of motion in polar coordinates. Answer: .
Express the conservation of angular momentum in cartesian coordinates. Answer: .
Express the conservation of angular momentum in polar coordinates. Answer: .
Show that (a) implies (c) and (b) implies (d) even if you didn’t note them a-priori. Answer: dot equation in (a) with to get which can be rewritten as (c); dot equation in (b) with to get which can be rewritten as (d).
Express the conservation of energy in cartesian coordinates. Answer: .
Express the conservation of energy in polar coordinates. Answer: =const.
Show that (a) implies (f) and (b) implies (g) even if you didn’t note them a-priori. Answer: dot equation in (a) with to get which can be rewritten as (f).
Find the general motion by solving the equations in (a). Describe all possible paths of the mass. Answer: , , where . The general motion is an ellipse.
Can the mass move back and forth on a line which is not the or axis? Answer: Yes. Consider , , and .
19.3.4 Cart and pendulum A mass hangs by two strings from a cart with mass . Before string BC is cut string AB is horizontal. The length of string AB is . At time all masses are stationary and the string BC is cleanly and quietly cut. After some unknown time the string AB is vertical.
What is the net displacement of the cart at ? Answer:
What is the velocity of the cart at ? Answer: .
What is the tension in the string at ? Answer: .
(Optional) What is ? [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.]
19.3.5 A dumbbell slides on a floor. Two point masses 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?]
Find the acceleration of point B just after the dumbbell is released.
Find the velocity of point A just before it hits the floor
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 , her tip from the vertical , and the gravitational constant is .
What is the path of her center of mass as she falls? (Show clearly with equations, sketches or words.)
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)?
Find a differential equation that only involves , its time derivatives, , , 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 and radius . 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 .
Find at a later time in terms of any or all of , , , and .
Does move sideways as the hoop falls and unrolls?
19.3.8 A model for a yo-yo consists of a thin disk of mass and radius 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 of the center of mass.
19.3.9 A uniform rod with mass pivots without friction about point in the -plane. A collar with mass slides without friction on the rod after the string connecting it to point is cut. There is no gravity. Before the string is cut, the rod has angular velocity .
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: , , , , and
Consider the special case . 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.
19.3.10 Assume the rod in the figure for problem 18.18 has polar moment of inertia . Assume it is free to rotate. The bead is free to slide on the rod. Assume that at the angular velocity of the rod is , that the radius of the bead is one meter and that the radial velocity of the bead, , is zero.
Draw separate free-body diagrams of the bead and rod.
Write equations of motion for the system. Answer:
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 . The quantity inside the derivative is angular momentum; thus, it is conserved and equal to a constant, say, , which can be found in terms of the initial conditions.
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: .
Write an expression for conservation of energy. Let the initial total energy of the system be, say, .
Answer: .
As 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.
Coupled motion of bead and rod and turntable.
19.3.11 A primitive gun rides on a cart (mass ) and carries a cannon ball (of mass ) on a platform at a height above the ground. The cannon ball is dropped through a frictionless tube shaped like a quarter circle of radius .
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.
Compute the velocity of the cannon ball relative to the cart.
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?
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.
Write the equations of motion for the system from part (b) in problem 19.3 as a set of first order differential equations. Answer:
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 and the lower one . 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: , , , , where 2 refers to the top wedge and 1 refers to the lower wedge.
19.3.14 Mass slides on an accelerating cart. 2D. A cart is driven by a powerful motor to move along the sloped ramp according to the formula: where , , and are given constants. The cart is held from tipping over. The cart itself has a sloped upper surface on which rests a mass (given mass ). The surface on which the mass rests is frictionless. Initially the mass is at rest with regard to the cart.
What is the force of the cart on the mass? [in terms of , and .] Answer: .
For what values of , , and is the acceleration of the mass exactly vertical (i.e., in the direction)? Answer: For , the acceleration of the mass is exactly vertical; , , and could be anything.
19.3.15 A thin rod of mass and length is pinned to a cart of mass , 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.
Determine the angular speed of the rod as it passes through the vertical position (at some later time). Answer: angular speed = .
Determine the displacement of the cart at the same instant. Answer: displacement = .
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?
19.3.16 As shown in the figure, a block of mass rolls without friction on a rigid surface and is at position (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:
What is the kinetic energy of the system?
What is the linear momentum of the system (momentum is a vector)?
19.3.17 A pendulum of length hangs from a cart. The pendulum is massless except for a point mass of mass at the end. The cart rolls without friction and has mass . 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: . The answer is in the special case when and . ]
19.3.18 Due to the application of some unknown force , the base of the pendulum A is accelerating with There is a frictionless hinge at A. The angle of the pendulum with the 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 . There is no gravity. What is ? (Answer in terms of , , , and .)
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 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 be a function of time that you can specify to try to make the mass swing progressively faster.
19.3.20 Using free-body diagrams and appropriate momentum balance equations, find differential equations that govern the angle and the vertical deflection of the system shown. Be clear about your datum for . Your equations should be in terms of and their time-derivatives, as well as , and .
19.3.21 The two blocks shown are released from rest at . 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 ?
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 and makes a perfectly plastic collision () with cart B which, together with bob C, is stationary prior to impact. Find the maximum vertical position of bob C, , after impact. The masses of the carts and pendulum bob are .
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.
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 plane. Write two equations from which one could find given and all mass and length quantities.
19.3.25 A pendulum is hanging from a moving support in the -plane. The support moves in a known way given by . For the following cases, find a differential equation whose solution would determine , measured clockwise from vertical; then find an expression for in terms of and :
with no gravity and assuming the pendulum rod is massless with a point mass of mass at the end,
as above but with gravity,
assuming the pendulum is a uniform rod of mass 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.
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: .
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.
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.
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.
FBD. Draw a FBD of the broom.
Momentum balance. Write the equation of angular momentum balance about the point instantaneously coinciding with the hinge.
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.]
Equation of motion. Using the angular momentum balance equation, write a governing differential equation for the tip angle. Answer: , where is the vertical displacement of your hand and is the angle of the broom from the vertical.
Simulation. Taking the hand motion as given, simulate on the computer the system you have found.
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 , 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.
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.
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 . The pendulum is released from rest at and . Just after release what are the values of and ? Answer in terms of other quantities. Answer: .
19.3.29 A rocker. A standing dummy is modeled as having massless rigid circular feet of radius rigidly attached to their uniform rigid body of length and mass . The feet do not slip on the floor.
Given the tip angle , the tip rate and the values of the various parameters (, , , ) find . [You may assume and are small.]
Answer: .
Using the result of (a) or any other clear reasoning find the conditions on the parameters (, , , ) 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 . is neutrally stable (like a wheel). (Someone questions this answer, so don’t trust it too much.)
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 . Assume that the mass of the wheel is , and that the polar moment of inertia about its center is (use if you want to get a better sense of the solution). Assume that just before collision number , the angular velocity of the wheel is , the kinetic energy is , the potential energy (you must clearly define your datum) is . Just after collision the angular velocity of the wheel is . The Kinetic Energy is , the potential energy (you must clearly define your datum) is . The wheel has spokes (pick 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.
What is the relation between and ?
What is the relation between and ?
Assume ‘rolling’ on level ground. What is the relation between and (the angular velocities just after two successive collisions?
Assume rolling down hill at slope . What is the relation between and ?
Can it be true that = ? About how fast is the wheel going in this situation?
As the number of spokes goes to infinity, in what senses does this wheel become like an ordinary wheel?