Here is a second approach to the kinematics of particle motion. Now, instead of using constant base vectors, we use time-varying base vectors. The discussion of polar coordinates started in Chapter 15 is completed here. Path coordinates, where one base vector is parallel to the velocity and the others orthogonal to that, are introduced. The challenging kinematics topic of relative motion is introduced in two stages: first using rotating base vectors connected to a moving rigid object and then using the more abstract notation associated with frame-dependent differentiation and the famous “five term acceleration formula.”
Many parts of practical machines and structures move in ways that can be idealized as straight-line motion (Chapter 14) or circular motion (Chapter 15). But often an engineer must analyze parts with more general motions as we began to study in Chapter 17.
In principle one can study all motions of all things using one fixed, say , coordinate system. If one knows the and coordinates of all points at all times then one can evaluate the linear momentum, the angular momentum, and their rates of change. In this way one can do all of mechanics. But when a machine has various parts, each moving relative to the other, it turns out it is helpful to make use of additional base vectors besides those fixed to a Newtonian (“fixed”) reference frame. That is, the formulas for velocity and acceleration are in some senses simplified (or clarified) by using moving base vectors. Most often these moving base vectors move with the moving parts.
You have seen the time-varying base vectors and , the polar coordinate base vectors used to describe circular motion. These are the ideas on which we build here. Altogether we discuss 4 approaches that use time varying base vectors:
Polar coordinates are generalized to include more than just circular motion;
Path coordinates and path base vectors;
General rotating base vectors and coordinate systems with an origin that moves are introduced; and
Formulas are presented for differentiation ‘in’ moving frames that don’t depend on any particular choice of base vectors.
The basic idea is to try to use coordinate systems that most simply describe the motions of interest, even if these coordinate systems are superficially confusing because they rotate and move. Time-varying base vectors are difficult at first. Like many shortcuts, they have a cost. In the end, however, they can aid intuition and simplify calculations.
As you learned in Chapter 15, when a particle moves in a plane while going in circles around the origin its position, velocity, and acceleration can be described like this:
These three equations say that the position is the distance from the origin times a unit vector towards the point; that the velocity is tangent to the circle of motion; and that the acceleration has a centripetal component proportional to the speed squared, and a tangential component, tangent to the circle of motion with magnitude equal to the rate of change of speed.
We will now generalize these results two different ways.
First we will use polar base vectors for non-constant .
Next we will use path base vectors to show that, in a sense to be explained, the 2nd and 3rd formulas above (for and ) apply to any wild motion in 2D or 3D.
As mentioned, in principle these new methods are not needed. We could just use one fixed coordinate system with base vectors , and and write the velocity and acceleration of a point at position as
as in Chapters 9-12. But, as for the circular motion of Chapter 15, rotating base vectors are helpful for simplifying some kinematics and mechanics problems.
The extension of polar coordinates to 3 dimensions as cylindrical coordinates is shown in fig. 18.1.
Rather than identifying the location of a point by its , and coordinates, a point is located by its cylindrical coordinates
the distance to the point from the axis,
the angle that the most direct line from the axis to the point makes with the positive direction,
the conventional coordinate of the particle,
and base vectors:
a unit vector normal to the axis that points from the axis
to the particle
(in 2-D ,
in 3-D
a unit vector in the direction of the shadow of in the plane)
a vector in the plane normal to (Formally ),
the conventional base vector.
The position vector of a particle is
The component of position, velocity and acceleration is the same in cylindrical coordinates as in Cartesian coordinates. If only two-dimensional problems are being considered then
As the particle moves, the values of its coordinates , and change as do the base vectors and .
Example: Oblong path
A particle that moves on the oblong path with is shown in fig. 18.2. The position vector is
Note that, unlike for circular motion, is not tangent to the particle’s path in general. The base vector is only tangent to the path at those points where the path is closest to, or furthest from, the origin (which, in this example, are also the points where the path crosses the and axes).
The velocity and acceleration are found by differentiating the position , taking account that the base vectors and also change with time just as they did for circular motion:
We find the velocity by taking the time derivative of the position, using the product rule of differentiation:
| (18.1) | |||||
This formula is intuitive. The velocity is the sum of three vectors: one, , due to moving towards or away from the axis; one, , having to do with the angle being swept; and in 3-D, one, , for motion perpendicular to the plane. In 2-D this is shown in fig. 18.3.
But, as has been emphasized before, this isn’t a new vector but just a new way of representing the same vector:
The vector can be represented in different base vector systems, in this case cartesian and polar.
To find the acceleration, we differentiate once again. The resulting formula has new terms generated by the product rule of differentiation.
| (18.2) | |||||
The acceleration for an arbitrary planar path is shown in fig. 18.4.
Four of the five terms comprising the polar coordinate formula for acceleration are easy to understand.
is just the acceleration due to the distance from the origin changing with time.
is the familiar centripetal acceleration.
is the acceleration due to rotation proceeding at a faster and faster rate. And
is the same as for Cartesian coordinates.
The Coriolis term.
The difficult term in the polar coordinate expression for the acceleration is the
term, called the Coriolis acceleration, after the civil engineer Gustave-Gaspard Coriolis who first wrote about it in 1835 (in a slightly more general context).
The presence of the ‘2’ in this term is due to the two effects from which it derives: 1 from the change of the term in the velocity and 1 from the change of the term in the velocity ().
The Coriolis acceleration occurs even if both and are constant. That is, a particle that moves at constant speed ( = constant) on a straight line that is itself rotating at constant rate ( constant) does not have a straight-line path, and thus has some acceleration. The Coriolis term catches this acceleration. Here is a situation in which the Coriolis term is the only non-zero term in the general polar-coordinate acceleration formula (eqn. (18.2)).
Example: The simplest Coriolis example
Consider, like above, a particle moving at constant speed along a line which is itself rotating at constant about the origin. Let’s look at a particle as it passes through the origin at time and a small amount of time later (see fig. 18.5). At time the direction of the scribed line has changed by an angle . So that, even if at that later time the direction of has changed so has changed by an amount . But the rotation of the line does continue, so the velocity includes a part in the direction with magnitude . That is is changed by both and by so
as predicted by the general polar coordinate acceleration formula.
But one need not get confounded by a desire to understand every term intuitively. Equation (18.2) is a way of describing the same acceleration we have described with cartesian coordinates. Namely,
Example: and are given functions.
Say and vary in time according to
with and as given constants. Then at any the position, velocity, and acceleration are (see fig. 18.6)
with polar base vectors
If we substituted these expressions for the polar base vectors into the expressions for , and we would get the same cartesian representation (a giant mess that we don’t show) that we would get from using and with , , and . That is , , and even if the representation is different.
Another way still to describe the velocity and acceleration is to use base vectors which are defined by the motion. In particular, the path base vectors used are:
the unit tangent to the path , and
the unit normal to the path .
Somewhat surprisingly at first glance, only two base vectors are needed to define the velocity and acceleration, even in three dimensions.
The base vectors can be described geometrically and analytically. Let’s start out with a geometric description.
As a particle moves through space it traces a path . At the instant of interest the path has a unique tangent line. The unit tangent is along this line in the direction of motion, as shown in fig. 18.7.
Less clear is that the path has a unique ‘kissing’ plane. One line on this plane is the tangent line. The other line needed to define this plane is determined by the position of the particle just before and just after the time of interest. Just before and just after the time of interest the particle is a little off the tangent line (unless the motion happens to be a straight line and the tangent plane is not uniquely determined). Three points, the position of the particle just before, just at, and just after the instant of interest determine the tangent plane.
Another way to picture the tangent plane is to find the circle in space that is tangent to the path and which turns at the same rate and in the same direction as the path turns. This circle, which touches the path so intimately, is called the osculating or ‘kissing’ circle. The tangent plane is the plane of this circle. (See fig. 18.7).
The unit normal is the unit vector which is perpendicular to the unit tangent and is in the tangent plane. It is pointed in the direction from the edge of the osculating circle towards the center of the circle as shown in fig. 18.7. For 2-D motion in the plane the osculating plane is the plane and the osculating circle is in the plane. The path base vectors are unit vectors that vary along the path, always tangent and normal to the path (see fig. 18.8).
The path of a particle can also be parameterized by arc length along the path, as explained in any introductory calculus text. So the path in space is also
, where the arc length is the path “coordinate”. The unit tangent is:
Using the chain rule with this is also
To define the unit normal let’s first define the curvature of the path as the rate of change of the tangent (rate in terms of arc length).
The unit normal is the unit vector in the direction of the curvature
Finally, the binormal is the unit vector perpendicular to and :
For 2-D motion the binormal is always in the direction. The radius of the osculating circle is
Note that, in general, the polar and path coordinate basis vectors are not parallel; i.e., is not parallel to and is not parallel to . For example, consider a particle moving on an elliptical path in the plane shown in fig. 18.9. In this case, the polar coordinate and path coordinate basis vectors are only parallel where the major and minor axes intersect the path.
Although it is not necessarily easy to compute the path basis vectors and , they lead to simple expressions for the velocity and acceleration:
| (18.3) | ||||
| (18.4) |
This formula for velocity is intuitive: velocity is speed times a unit vector in the direction of motion. The formula for acceleration is more interesting. It says that the acceleration of any particle at any time is given by the same formula as the formula for acceleration of a particle going around in circles at non-constant rate.
The first term is tangent to the path (also tangent to the osculating circle), as shown in fig. 18.10. If we draw the osculating circle at the point of interest and fix it in space, then this first term is in a polar coordinate system centered at the center of that fixed circle.
The term is directed towards the center of the osculating circle. This term is associated with change of direction and does not vanish even if the speed is constant. This normal acceleration is perpendicular to the path.
That is, the two terms in the path-coordinate acceleration formula correspond exactly to the two terms for acceleration for a particle going in circles.
Example: Estimating the direction of acceleration
By looking at the path of a particle (e.g., see fig. 18.11) and just knowing whether it is speeding up or slowing down one can estimate the direction of the acceleration.
If the particle is known to be speeding up at A then so the tangential acceleration is in the direction of the velocity. Without thinking about the normal acceleration you know that the acceleration vector is pointed in the half plane of directions shown.
If at point B nothing is known about the rate of change of speed, you still know that the acceleration must be in the half plane shown because that is the direction of and .
If at C you know that the particle is slowing down then you know that . But you also can see the curve is to the right so is to the right. So the acceleration must be in the quadrant shown.
One can use the information about the curvature to further restrict the possible accelerations at point A also. At point B there is nothing more to know unless you know how the speed is changing with time.
Earlier we found the curvature by assuming the particle’s path was parameterized by arc length . A second way of calculating the curvature (and then the unit normal ) is to calculate the normal part of the acceleration. First calculate the acceleration. Then subtract from the acceleration that part which is parallel to the velocity.
The normal acceleration is so
Assume that you know the position as a function of time in either cartesian or polar coordinates. Then, say, at a particular time of interest when the particle is at , you can calculate the velocity of the particle using:
and the acceleration using
From these expressions we can calculate all the quantities used in the path coordinate description. So we repeat what we have already said but in an algorithmic form. Here is one set of steps one can follow. This recipe is of little practical use, but does show that the motion explicitly determines the path base vectors as well as the osculating circle.
Calculate .
Calculate .
Calculate .
Calculate
Calculate the radius of curvature as .
Write a parametric equation for the osculating circle as
where is the parameter used to parameterize the points on the circle of the point on the curve . As ranges from 0 to the point goes from around the circle and back. The plane of the osculating circle is determined by and . For planar curves, the osculating circle is in the plane of the curve.
Example: Particle on the rim of a tire
A particle P on the rim of a tire whose center is moving at constant speed has position given by
where the origin is at the ground contact time . When the particle is at its highest point and
At that midpoint
as shown in fig. 18.12. The osculating circle has 4 times the radius of the tire. Note the intimacy of the osculating circle’s kiss with the cycloidal path.
See also the inside back cover, table II, row 3 for future reference.
For 2-D problems just set , and in these equations.
See the inside back cover table II, row 4 and the text under the table for future reference:
Both polar coordinates and path coordinates define base vectors in terms of the motion of a particle of interest relative to a fixed coordinate system.
SAMPLE 18.1
Acceleration in polar coordinates. A bug walks along the spiral section of a natural shell. The path of the bug is described by the equation where and . The bug’s radial distance from the center of the spiral is seen to be increasing at a constant rate of . Find the and components of the acceleration of the bug at .
Solution In polar coordinates, the acceleration of a particle in planar motion is
Since we know the position of the bug,
Since the radial distance R of the bug is increasing at a constant rate , , that is,
Therefore,
Now substituting , and in the above expression, we get
But, at
therefore,
Answer:
SAMPLE 18.2 Going back and forth between and . Given the position of a particle in polar coordinates and its radial and angular velocity and radial and angular acceleration , find and . Also, find the inverse relationship.
Solution
Polar to Cartesian: In polar coordinates, we are given , and . We need to find , and . Let us consider the velocity first. The velocity of a point is in cartesian coordinates and in polar coordinates. Thus,
where and . Dotting this equation with and , respectively, we get
| or, | |||||
| or, | |||||
| (18.13) | |||||
Thus given and at , we can find and . Similarly, from the acceleration formula, , we derive
| (18.24) |
It is not necessary to split the terms on the right hand side. We could have kept them together as and but we split them to keep the radial acceleration term and angular acceleration in evidence.
Cartesian to Polar: Given , and at , we can now find , and easily by inverting eqn. (18.13) and eqn. (LABEL:eq:10.1.aptoac):
| (18.34) |
| (18.45) |
Note that in eqn. (LABEL:eq:10.1.actoap) we need and in order to compute and . This, however, is no problem since we have and from eqn. (18.34). Of course, and are required too, which are easily computed as and .
SAMPLE 18.3 Velocity in path coordinates. The path of a particle, stuck at the edge of a disk rolling on a level ground with constant speed, is called a cycloid. The parametric equations of a cycloid described by a particle is , where is a dimensionless time. Find the velocity of the particle at
,
, and
and express the velocity in terms of path basis vectors .
Solution The position of the particle is given:
| (18.49) |
In terms of path basis vectors, the velocity is given by
Here,
| (18.50) |
Substituting the values of in equations 18.49 and 18.50 we get
at :
Answer:
at :
Answer:
at :
Answer:
SAMPLE 18.4 Path coordinates in 2-D. A particle traverses a limacon with constant angular speed .
Find the normal and tangential accelerations ( and ) of the particle at .
Find the radius of the osculating circle and draw the circle at .
Solution
The equation of the path is
The path is shown in Fig. 18.14. Since the equation of the path is given in polar coordinates, we can calculate the velocity and acceleration using the polar coordinate formulae:
| (18.51) | |||||
| (18.52) |
So, we need to find for computing and . From the given equation for
which give velocity and acceleration at any . Now substituting we get the velocity and acceleration at the desired point:
Thus we know the velocity and the acceleration of the particle in polar coordinates. Now we proceed to find the tangential and the normal components of acceleration (acceleration in path coordinates). In path coordinates
where and are unit vectors in the directions of the tangent and the principle normal of the path. We compute these unit vectors as follows.
| So, | ||||
| and | ||||
| Therefore, | ||||
Thus,
Answer:
In path coordinates the acceleration is also expressed as
where is the radius of the osculating circle. Since we already know the speed and the normal component of acceleration we can easily compute the radius of the osculating circle.
Answer:
18.1.1 A particle moves along the two paths (1) and (2) as shown.
In each case, determine the velocity of the particle in terms of , , and . Answer: plot(2): and .
Find the and coordinates of the path as functions of and or and . Answer: plot(2): and or and
18.1.2 For the particle path (1) in problem 18.17, find the acceleration of the particle in terms of , , and .
18.1.3 For the particle path (2) in problem 18.17, find the acceleration of the particle in terms of , , and .
18.1.4 A body moves with constant velocity in a straight line parallel to and at a distance from the -axis.
Calculate in terms of , and .
Calculate the component of acceleration.
Carefully define, using a sketch and/or words, any variables, coordinate systems, and reference frames you use. Express your answers using any convenient coordinate system (just make sure its orientation has been clearly defined).
18.1.5 Picking apart the polar coordinate formula for velocity. This problem concerns a small mass that sits in a slot in a turntable. Alternatively you can think of a small bead that slides on a rod. The mass always stays in the slot (or on the rod). Assume the mass is a little bug that can walk as it pleases on the rod (or in the slot) and you control how the turntable/rod rotates. Name two situations in which one of the terms is zero but the other is not in the two term polar coordinate formula for velocity, . You should thus gain some insight into the meaning of each of the two terms in that formula. Answer: One situation: ;the case of a particle constrained to move in a circle.
18.1.6 Picking apart the polar coordinate formula for acceleration. Reconsider the configurations in problem 18.18. This time, name four situations in which all of the terms, but one, in the four term polar coordinate formula for acceleration, , are zero. Each situation should pick out a different term. You should thus gain some insight into the meaning of each of the four terms in that formula. Answer: One situation: ; the case of no rotation ().
18.1.7 The two differential equations:
have the general solution
where , , , and are arbitrary constants. This solution could be checked by plugging back into the differential equations — you need not do this (tedious) substitution. The solution describes a curve in the plane. That is, if for a range of values of , the values of and were calculated and then plotted using polar coordinates, a curve would be drawn. What can you say about the shape of this curve?
[Hints:
Actually make a plot using some random values of the constants and see what the plot looks like.
Write the equation for a particle in polar coordinates and think of a force that would be relevant to this problem.
The answer is something simple. ]
18.1.8 A car driver on a very boring highway is carefully monitoring her speed. Over a one hour period, the car travels on a curve with constant radius of curvature, , and its speed increases uniformly from to . What is the acceleration of the car half-way through this one hour period, in path coordinates?
18.1.9 Find expressions for , , , , and the radius of curvature , at any position (or time) on the given particle paths for
18.1.10 A particle travels at non-constant speed on an elliptical path given by . Carefully sketch the ellipse for particular values of and . For various positions of the particle on the path, sketch the position vector ; the polar coordinate basis vectors and ; and the path coordinate basis vectors and . At what points on the path are and parallel(or and parallel)?
In this section you will learn about rotating reference frames, how to take the derivative of a vector ‘in’ a rotating frame, and how to use that derivative to find the derivative in a Newtonian or fixed frame. We start by showing the alternative, just using one frame with one set of fixed base vectors.
To motivate the sections that follow we first show the “fixed base vector” method.
Consider the task of determining the acceleration of a bug walking at constant speed as it walks on a straight line marked on the surface of a tire rolling at constant rate. Artificial as this problem seems, it is similar to the sort of calculation needed in the kinematics of mechanisms. For now, imagine you really care how strong the bug’s legs need to be to hold on (unreasonably neglecting air friction). So knowing the bug’s acceleration determines the net force on it by . Now we try to find by taking two time derivatives of position.
If we want to avoid using rotating base vectors we have to write an expression for the position of the bug in terms of and components. Choosing a suitable origin of the coordinate system we have
| (18.53) |
To find the velocity we take the time derivative, taking account that both and are functions of . Thus, for example looking at the term both the product rule and chain rule need be applied. Proceeding we get
| (18.54) | |||||
To get the acceleration of the bug we differentiate one more time. This time we use the product rule and chain rule again, but get to use the simplification for this problem that the rolling and bug walking are at constant rate so and :
| (18.55) | |||||
which is a bit of a mess. We could regroup the terms, but there would still be 6 of them.
The moving-reference-frame methods that follow don’t change this answer. But they give a somewhat simpler derivation. And they also group the terms in a physically meaningful way. One would be hard pressed to make sense of all the terms in eqn. (18.55). With the time-varying base vector methods below we can interpret the terms.
A reference frame is a coordinate system††margin: A fine point for experts: reference frame vs coordinate system. There is a semantic debate about the degree to which the phrases “coordinate system” and “ reference frame” are synonymous. For simplicity we take them to mean the same thing. An alternative definition distinguishes reference frames from coordinate systems. It turns out that coordinate systems which are rotated, but not rotating, with respect to each other both calculate the same time-derivative of a given vector. Because these coordinate systems are equivalent in this regard they are sometimes called the same reference frame. That is, some people consider one reference frame to be the set of all coordinate systems that are glued to each other, no matter what their position or orientation. In this way of thinking, a frame is made manifest by the use of one of its coordinate systems, but no particular coordinate system is unique to the frame. . It has an origin and a set of preferred mutually orthogonal directions represented by base vectors. You can think of a reference frame as a giant piece of graph paper, or in 3-D as a giant jungle gym, that permeates space. It has the look of a wire frame. Because we will use various frames, we name them. We always have one frame that we think of as fixed for the purposes of Newtonian mechanics. We call this frame (or sometimes ). Most often we choose a frame that is ‘glued’ to the ground with an origin at a convenient point and with at least one base vector lined up with something convenient (e.g., up, sideways, along a slope, along the edge of an important part, etc.). is a frame in which the mechanics laws we use are accurate. We define it by its origin and the direction of its coordinate axes, thus we would write
is or is .
where we would generally have a picture showing the position of the origin and the orientation of the coordinate axes (see fig. 18.21).
When we write casually ‘position ’ of a point we mean . When we write ‘velocity ’ we mean as calculated in . That is, if then we define the derivative of with respect to in as
The script shows explicitly that when we take the time derivative of the vector we take the time derivative of its components, using the components associated with and holding constant the base vectors associated with . That is
are just fancy ways of writing what we have been calling .
The elaborate notation just makes explicit how is defined. The only need for this elaborate notation is if there is ambiguity. There is only ambiguity if more than one reference frame is used in a given problem.
Let’s add a second reference frame called glued to and oriented with the roof of the building. We will always use script capital letters ( or ) to name reference frames. We define by writing
is or is
and by drawing a picture (see fig. 18.22).
This new frame, as we have drawn it, is also a good Newtonian or fixed frame. So we could write all positions using the coordinates and base vectors and then proceed with all of our mechanics equations. The only confusion being that gravity doesn’t point in the direction, but in some crooked direction relative to and (which we would have to work out from the angle of the roof).
Although one hardly notices when using just a single fixed frame, we actually use frames for three somewhat distinct purposes:
To define a vector. For example if we were tracking the motion of a cannon ball at P we could define its position vector as , using frame to define . Or we could define as using frame to define .
To assign coordinate values to a given vector. For example, the vector could be written as
Alternatively, if we just want to look at the components of a given vector we use to indicate the components of the vector in the square brackets, , using the base vectors of . Thus
where we have used to put the components in their standard column form (although this is a picky detail). Note that although that
That is, the components of a vector are different when expressed in different frames. Two different lists of numbers represent the same vector. The magnitude, however, is independent of frame (, to two decimal places).
To find the rate of change of a given vector. The position of P relative to A changes with time. We can calculate this rate of change two different ways. First using frame
if we are clear in our minds that and are the coordinates of P relative to A. But we can also calculate the rate of change of the same vector using frame as
For the two frames and because the two frames are not rotating relative to each other. Specifically, for and the formula for finding and from and does not involve time. Similarly, the formulas for finding and from and do not involve time. For frames that are rotated with respect to each other but not rotating, the two time derivatives of a given vector are related the same way the vector itself is related to itself in the two frames. The vectors are the same but their coordinates are different. That is, for rotated but not relatively rotating frames
| and | ||||
| and |
Going back and forth between these three uses of frames with ease is one of the advanced skills of a person who can analyze the dynamics of complex systems (And being confused about the distinctions is an almost universal part of learning advanced dynamics).
Example: Two fixed frames and
Consider and both to be fixed to the ground. Let’s look at where P is moving up at constant rate (see fig. 18.23). First look at the position using both frames:
Now look at the rate of change of position using both frames. First :
Then the rate of change of as calculated in :
You can quickly verify that by noting that and .
So long as is not rotating with respect to then the rate of change of a given vector is the same in both reference frames.
Now look at a third reference frame that is glued to the roof of the car as it starts up hill (see fig. 18.24). We define by the origin of its coordinate system 0′′ and its time-varying base vectors and . The issues with defining a vector with and with writing components using are the same as for . However taking the time derivative of a given vector in is different than taking the time derivative in or because is rotating relative to them.
Because dynamics involves the time derivatives of so many different vectors (e.g. , , , , and ) it is easier to think about the derivative of some arbitrary or general vector, call it , and then apply what we learn to these other vectors.
Recalling our three uses of frames:
To define a vector.
To express the coordinates of a given vector.
To take the time derivative of a vector.
we see that items [III.] and [I.] can be combined. That is, once a vector is defined clearly by some means then we can define a new vector as the derivative of that vector in, say, moving frame . Once this new vector is defined it can be expressed in terms of the coordinates of any convenient frame.
Example: Derivative in a moving frame of a constant vector
Consider as the relative position vector of the points A and P that do not move in the fixed frame . That is, the points A and P don’t move in the ordinary sense of the words (see fig. 18.25). Now also look at the frame that is rotating with respect to at the rate . We have
So we can now calculate the derivative in each frame by holding the corresponding base vectors as constant. So
That is, the stationary vector and the rotating frame define a new vector, the derivative of in . This is also , the difference between the velocity of P and the velocity of A in the frame . This new vector can be expressed in any coordinate system of choice for example the system. So we wrote above
which looks mixed up but isn’t. The frame is used to help define a vector which is then expressed in the coordinates of .
Given a new vector , the derivative of as calculated in a rotating frame , one calculation of common use is the determination of the derivative of the same vector in the fixed frame .
First think of a line segment that is marked between two points that are glued to a moving frame . We know (at least in 2-D and for fixed axis rotation) that
Likewise for any vector which is fixed in . It is especially useful to apply this formula to unit base vectors, so
| (18.57) | ||||
In some minds, Eqns. 18.2 are the core of rigid body kinematics. Box 18.2 on 18.2 shows how these relations give ‘the Q dot’ formula: For any time dependent vector
| (18.58) |
or more simply, but less explicitly,
where is the time derivative of relative to the moving frame of interest (in this case ). The ‘Q dot’ formula says that
The derivative of a vector with respect to a Newtonian frame (or ‘absolute derivative’) can be calculated as the derivative of the vector with respect to a moving frame , plus a term that corrects for the rotation of frame relative to frame .
Note that if is a constant in the frame , like the relative position vector of two points glued to , then and the formula reduces to
The formula 18.58 is useful for the derivation of a variety of formulas and is also useful in the solution of problems.
While we have shown how to use this formula to calculate the rate of change of a vector with respect to a Newtonian frame, the formula can be used to calculate its rate of change with respect to a non-Newtonian frame. Letting and be two possibly non-Newtonian frames, the formula for the rate of change of with respect to frame is
| (18.59) |
Both and could be non-Fixed (non-Newtonian).
For a vector fixed in ,
For any time dependent vector ,
Some examples of applying the formula are:
Absolute velocity of a point relative to :
Rate of change of a rotating unit vector which is fixed in :
One way to calculate velocity, acceleration is to express the position of a particle in terms of a combination of base vectors, some of which change in time. Velocity and acceleration are then determined by directly differentiating the expression for position, taking account that the base vectors themselves are changing. This method is sometimes convenient for bodies connected in series, one body to the next, etc.
The overall approach is as follows:
Glue a coordinate system to every moving body. If needed, also create moving frames that move independently of any particular body.
Call the basis vectors associated with these frames , , for the fixed frame ; , , for the moving frame ; and , , for the moving frame , etc.
Evaluate all of the relative angular velocities; , , etc. in terms of the scalar angular rates , , etc. and the base vectors glued to the frames.
Express all of the absolute angular velocities in terms of the relative angular velocities.
Differentiate to get the angular accelerations using, for example,
Write the position of all points of interest in terms of the various base vectors.
Differentiate the position to get the velocities (again using , etc.)
Differentiate again to get acceleration.
First, reconsider the bug crawling on the tire in fig. 18.27.
Example: Absolute velocity of a point moving relative to a moving frame: Bug crawling on a tire
We write the position of the bug in terms of the various basis vectors as
To get the absolute velocity of the bug at the instant shown, we differentiate the position of the bug once, using the product rule and the rates of change of the rotating basis vectors with respect to the fixed frame, to get
| (18.60) | |||||
Example: Absolute acceleration of a point moving relative to a moving frame (2-D): Bug crawling on a tire, again
Differentiating equation 18.60 from the example above again, we get the absolute acceleration of the bug at the instant shown,
In the varying base vector method, we calculate the velocity of a point by looking at the position as the sum of two position vectors, one of which is expressed in the moving base vectors. We then differentiate the position, taking account that the base vectors of the moving frame change with time. In general
We could calculate similarly using a combination of the product rule of differentiation and the facts that , , and ,
and would get a formula with 15 non-zero terms.
Box 18.1 The formula
We think about some vector as a quantity that could be represented by an arrow. We can write using the coordinates of the fixed Newtonian frame with base vectors : . Similarly we could write in terms of the coordinates of some moving and rotating frame with base vectors : . Now of course
Similarly, so long as what we mean by is its derivative in a fixed frame. That is, we use as an informal notation for . We can calculate the same way we have from the start of the book, namely,
We didn’t have to use the product rule of differentiation because the unit vectors , , and , associated with a fixed frame, are constant in time.
What if we wanted to use the coordinate information that was given to us by a person who was moving and rotating with the moving frame ? Now we calculate taking account that the base vectors change in time.
The first term in the product rule is just the derivative of in the moving frame . That is, is calculated by differentiating the components in holding the base vectors in fixed. The second term depends on evaluating , , and . We know (at least for 2-D and for fixed-axis rotation in 3-D) that
| (18.62) | ||||
Eqns. 18.2 are the core of rigid body kinematics.
Now we can go back to the second group of terms in Eqn. 18.2.
| ? | ||||
Going back to Eqn. 18.2 we get the desired result:
| (18.63) |
or more simply, but less explicitly,
| (18.64) |
Here is a geometrical ‘derivation’ of the formula in two dimensions. Referring to the figure at right, we look at a vector at two successive times. We then look at how seems to change in a frame that rotates slightly as changes. The picture shows how to account for the difference between the change of as perceived by the two different frames.
In detail the parts (a) to (e) of the picture show the following.
Part (a) shows a vector at time .
Part (b) shows at time and the change in , .
Part (c) is like (a) but shows a moving body or frame .
Part (d) shows the change in , , that would occur if were fixed (constant) in .
Part (e) shows the change in that would be observed in the moving frame .
Part (f) shows the net change in , , that is the same as that in (b) above; here, it is shown as the sum of the two contributions from (d) and (e).
Thus, using , for small is composed of two parts: (1) the observed in , and (2) the change in which would occur if were constant in and thus rotating with it. Dividing by gives the ‘ formula’, .
Two different looks at the change in the vector , , over a time interval .
SAMPLE 18.5
Acceleration of a point moving in a rotating frame. It is given that the arm OAB rotates with counterclockwise angular acceleration and at the instant shown the angular speed . Also, at the same instant, the particle P is sliding (roughly) down with speed and acceleration . Find the absolute acceleration of the particle at the given instant. Take in the figure.
Solution Let us attach a body frame to the rigid arm OAB. For calculations we fix a coordinate system in this frame such that the origin of the coordinate system coincides with O, and at the given instant, the axes are aligned with the inertial coordinate axes . Since is fixed in the frame and rotates with the rigid arm with and , the basis vectors , and rotate with the same and .
In the rotating (primed) coordinate system,
Now, we use the formula to evaluate and , i.e.,
Also, note that is constant since in frame , the motion of the particle is always along the tube, i.e., along the negative axis (see Fig. 18.30). Thus, , , , and . Substituting these quantities in , we get:
| (18.65) | |||||
Now substituting , , , and noting that at the given instant, we get:
We can find by differentiating Eq. (18.65) and noting again that at the given instant
:
Answer:
SAMPLE 18.6
Rate of change of unit vectors. A circular disk is welded to a rigid rod AB. The rod rotates about point A with angular velocity . A frame is attached to the disk and therefore rotates with the same . Two coordinate systems, and are fixed in frame as shown in the figure.
Find the rate of change of unit vectors using the formula.
Express the and vectors in terms of and and verify the results obtained above for and by direct differentiation.
Solution Since the disk is welded to the rod and frame is fixed in the disk, the frame rotates with .
To find the rate of change of the unit vectors using the formula, we substitute the desired unit vector in place of in the formula (eqn. (18.58)). For example,
It should be clear that , since does not change with respect to an observer sitting in frame . Therefore,
Similarly,
Answer:
Since and , we get their rates of change by direct differentiation as
Here we have used the fact that , the angle between the unit vectors and , remains constant during the motion. The results obtained are the same as in part (a).
SAMPLE 18.7
Rate of change of a position vector. A rigid rod OAB rotates counterclockwise about point O with constant angular speed . A collar C slides out on the bent arm AB with constant speed with respect to the arm. Find the velocity of the collar using the formula.
††margin:
Solution Let be the position vector of the collar. Then the velocity of the collar is . Let the rod OAB be the rotating frame . Now we can find using the formula:
To compute , let us first find , the rate of change of as seen in frame (this term represents the velocity of the collar you see if you sit on the rod and watch the collar; also called ).
Note that the vector does not change in frame since both its magnitude, , and direction, , remain fixed in . Therefore,
because does not change in and = speed of the collar with respect to the arm (see Figure 18.34). Thus,
Hence,
where we have used the fact that at the given instant, and
Answer:
18.2.1 Express the basis vectors () associated with axes and in terms of the standard basis () for .
18.2.2 Body frames are frames of reference attached to a body in motion. The orientation of a coordinate system attached to a body frame is shown in the figure at some instant of interest. For , express the basis vectors in terms of the standard basis ().
18.2.3 Find the components of (a) in the rotated basis () and (b) in the standard basis (). ( and are in the same direction and and are in the plane.)
18.2.4 A particle travels in a straight line in the -plane parallel to the -axis at a distance in the positive direction. The position of the particle is denoted by . The angle of measured positive counter-clockwise from the axis is decreasing at a constant rate with magnitude . If the particle starts on the axis at , what is in cartesian coordinates?
18.2.5 Given that and that , find
in terms of and ,
in terms of and .
18.2.6 A bug walks on a turntable. In polar coordinates, the position of the bug is given by , where the origin of this coordinate system is at the center of the turntable. The coordinate system is attached to the turntable and, hence, rotates with the turntable. The kinematical quantities describing the bug’s motion are , , , and . A fixed coordinate system has origin at the center of the turntable. As the bug walks through the center of the turntable:
What is its speed?
What is its acceleration?
What is the radius of the osculating circle (i.e., what is the radius of curvature of the bug’s path?).
Now that we have some comfort with moving frames we can develop formulas that are not so strongly attached to base vectors. That is, we take account that the base vectors rotate with the frame, but develop formulas that don’t use the base vectors explicitly. Thus the formulas we develop here work equally for any frame that is glued to the rotating frame of choice, independent of its orientation.
Imagine that you know the absolute velocity of some point on an object , say the center of a car tire and the angular velocity of the tire, . Finally, imagine you also know the relative velocity of point P, , say of a bug crawling on the tire.
If the frame is translating or rotating, the velocity of particle P relative to the frame is not the absolute velocity (the velocity relative to a Newtonian frame). The absolute velocity in this case is , or more simply , or more simply still, just . We want to know the relationship between the relative velocity and the absolute velocity (otherwise known as just ).
Let’s start by looking at the position. The position of a point P that is moving is:
where is the origin of a coordinate system which is glued to the rigid object, as shown in fig. 18.39.
To find the absolute velocity of point P we will use the formula, equation 18.58, for computing the rate of change of a vector. The velocity of P is the rate of change of its position. Here, we use
The formula 18.58 was used in the calculation to compute
| (18.66) |
Thus, the ‘three term velocity formula’.
| (18.67) |
Another way to write the formula for absolute velocity is as
where P’ is a point glued to which is instantaneously coincident with P, so the absolute velocity of is
| (18.68) |
Reconsider the bug crawling on the tire, object , in fig. 18.40. To find the absolute velocity of the bug, we need to be concerned with how the bug moves relative to the tire and how the tire moves relative to the ground.
Example: Absolute velocity of a point moving relative to a moving frame (2-D): Bug crawling on a tire, again
At the instant of interest, the direction of the bug’s absolute velocity depends upon the relative magnitudes of and as well as the orientation of and .
As we noted earlier, another way to write the formula for absolute velocity is
where, in the example above, and . At the instant of concern, we can think of the absolute velocity of the bug as the velocity of the mark labeled under the bug plus the velocity of the bug relative to the tire.
We would like to find acceleration of a point using information about its motion relative to a moving frame. The result, the ‘five term acceleration formula’ is the most complicated formula in this book. (For reference, it is in Table II, 5c).
The acceleration of a point relative to an object or frame is the acceleration you would calculate if you were looking at the particle while you translated and rotated with the frame and took no account of the outside world. That is, if the position of a particle P relative to the origin of a coordinate system in a moving frame is given by:
then the acceleration of the particle P relative to the frame is:
That is, the acceleration relative to the frame takes no account of (a) the motion of the frame or of (b) the rotation of the base vectors with the frame to which they are fixed.
Reconsider the bug labeled point P crawling on the tire, object B, in fig. 18.42.
Example: Acceleration relative to a frame (2-D): Bug crawling on a tire, again
If we are sitting on the tire, all that we see is the bug crawling in a straight line at non-constant rate relative to us. Thus, its acceleration relative to the tire is
So, at the instant of interest, the bug has an acceleration relative to the tire frame parallel to the -axis.
Imagine that you know the absolute acceleration of some point at the center of a frame , say the center of a car tire. Imagine you also know the angular velocity of the tire, , and the angular acceleration, . Then, you can find the absolute acceleration of a piece of gum labeled point stuck to the sidewall (see fig. 18.44).
If we start with the equation 18.68 for the absolute velocity of a point glued to a moving frame on page 18.68 and differentiate with respect to time, we get the absolute acceleration of a point fixed in a moving frame as follows:
| (18.69) | |||||
Example: Absolute acceleration of a point glued to a moving frame (2-D): Bug crawling on a tire, again
Here, the acceleration of point P′ glued to the tire, relative to the tire is zero, (see fig. 18.44). The angular velocity of the wheel with respect to the ground is . The angular speed is increasing at a rate . Thus, . The position of relative to is .
In this example, the absolute acceleration is due to:
the increase in the translational speed of the tire relative to the ground (acceleration of origin of moving frame),
its going in circles at a non-constant rate about point relative to the ground (‘tangential term’), and
‘centripetal term’ towards the origin of the moving frame. (In three-dimensional problems, this term is directed towards an axis through that goes through O’).
If we start with the equation for absolute velocity 18.3 on page 18.3 and differentiate with respect to time we get the absolute acceleration of a point P using a moving frame . To do this calculation we need to use the product rule of differentiation. Refer to the formula, eqn. (18.63) on page 18.63. Here is the calculation:
The collection of terms is the acceleration of a point which is glued to body and is instantaneously coincident with P. It is the same as using in equation 18.69. To repeat, the result is
The ‘five term’ acceleration formula
| (18.70) | ||||
| (18.71) |
This ‘five-term-acceleration’ formula is both famous and infamous. It’s famous because it is given a lot of emphasis by some instructors
,
and infamous because it takes some getting used to. Eqn. 18.71 is the ‘three term acceleration formula’. It combines the first three terms in the 5-term formula and interprets them as the acceleration of the point P′ on B that instantaneously coincides with P. The best way to get used to the five term acceleration formula is to find situations where some of the terms drop out.
Reconsider the bug labeled point P crawling on the tire, object , in fig. 18.44. To find the absolute acceleration of the bug we need to think about how the bug moves relative to the tire and how the tire moves relative to the ground.
Example: Absolute acceleration of a point moving relative to a moving frame (2-D): Bug crawling on a tire, again
Referring to the five term acceleration formula, equation 18.70 on page 18.70, the absolute acceleration of the bug is
So, at the instant of interest, the bug’s absolute acceleration is due to:
the translational acceleration of the tire,
the centripetal acceleration of going in circles of radius about the center of the tire as it rolls, , pointing at the center of the tire,
the tangential acceleration of going in circles about the center of the tire as the tire rolls at non-constant rate, ,
the acceleration of the bug relative to the tire as it crawls on the line, , and
the Coriolis acceleration caused, in part, by the change in direction, relative to the ground, of the velocity of the bug relative to the tire, .
Items 1, 2 and 3 sum to be the acceleration of point on the tire but instantaneously coinciding with moving point P.
We can now give a different interpretation of the expressions we have been using and .
Rather than thinking of as the difference between and we can think of as the where is a frame with origin that moves with point A and which has no rotation rate relative to . That is
Similarly,
Box 18.2 Relation between moving frame formulae and polar coordinate formulae
A similarity exists between the polar coordinate velocity formula
and the second two terms in the ‘three-term’ velocity formula
In fact, we have tried to build your understanding of moving frames by means of that connection.
Similarly, the polar coordinate formula for acceleration
is somehow closely linked to the last four terms of the ‘5-term’ acceleration formula
Let’s make these connections explicit. Imagine a particle P moving around on the -plane.
Let’s create a moving frame with rotating coordinate system attached to it whose origin is coincident with origin of a coordinate system attached to a fixed frame . Let this frame rotate in exactly such a way so that the particle is always on the -axis. So, in this frame, , , and . Also, the frame motion is characterized by , , and . So, if we plug in the three-term velocity formula, we get
which is the polar coordinate velocity formula.
Similarly, if we plug into the five-term acceleration formula, we get
Again, we recover the appropriate polar coordinate formula.
We have just shown how the polar coordinate formulae are special cases of the relative motion formulae.
In problems where we want the rotating frame to be a
rotating object on which a particle moves, the polar coordinate formulae only
correspond term by term with the relative motion formulae if the particle path
is a straight radial line fixed on a object, as in the example
of a bug walking on a straight line scribed on the
surface of a rotating CD or a bead sliding in a tube rotating about an axis perpendicular to the tube.
SAMPLE 18.8
A ‘T’ shaped tube is welded to a massless rigid arm OAB which rotates about O at a constant rate . See also Sample 18.29. At the instant shown a particle P is falling down in the vertical section of the tube with speed . Find the absolute velocity of the particle. Take in the figure.
††margin:
Solution Let us attach a frame to arm OAB. Thus rotates with OAB with angular velocity where . To do calculations in we attach a coordinate system to at point O. At the instant of interest the rotating coordinate system coincides with the fixed coordinate system . (Since the entire motion is in the -plane, the -axis is not shown in the figure). Let be a point coincident with but fixed in . Now,
where
and
| Velocity relative to the frame | ||||
Thus,
Answer:
Comments: The kinematics calculation is equivalent to the vector addition shown in Figure 18.48. The velocity of P is the sum of and .
SAMPLE 18.9
Acceleration of a point in a rotating frame. Consider the rotating tube of Sample 18.3 again. The arm OAB rotates with counterclockwise angular acceleration and, at the instant shown, its angular speed . Also, at the same instant, the particle P falls down with speed and acceleration . Find the absolute acceleration of the particle at the given instant. Take in the figure.
Solution
We consider a frame , with coordinate axes , fixed to the arm OAB and thus rotating with and . The acceleration of point P is given by
where
| acceleration of a point P′ that is fixed in | ||||
| and at the moment coincides with P, | ||||
| Coriolis acceleration, and | ||||
| acceleration of P relative to frame . |
Now we calculate each of these terms separately. For calculating , imagine a rigid rod from point O to point P, rotating with the frame . Mark the far end of the rod as P′ (same as point P). The acceleration of this end of the rod is . To find the relative terms and , freeze the motion of the frame at the given moment and watch the motion of point P. The non-intuitive term has no such simple physical interpretation but has a simple formula. Thus,
Adding the three terms together, we get
Answer:
Note that the single term encompasses three terms of the five term acceleration formula.
SAMPLE 18.10
A small collar P is pinned to a rigid rod AB at length along the rod. The collar is free to slide in a straight track on a disk of radius . The disk rotates about its center O at a constant . At the instant shown, when and the collar is at a distance in the track from the center O, find
the angular velocity of the rod AB and
the velocity of point P relative to the disk.
Solution
We will think of P in two ways: one as attached to the rod and the other as sliding in the slot. First, let us attach a frame to the disk. Thus rotates with the disk with angular velocity . We attach a coordinate system to at point O. At the instant of interest, the rotating coordinate system coincides with the fixed coordinate system . Now let us consider point which is fixed on the disk (and hence in ) and coincides with point P at the instant of interest. We can write the velocity of P as:
where
In the last expression, , we do not know the magnitude of and hence have left it as an unknown , but its direction is known because has to be along the track and the track at the given instant is along the -axis. Thus,
| (18.72) |
Now let us consider the motion of rod AB. Let be the angular velocity of AB at the instant of interest where is unknown. Since P is pinned to the rod, it executes circular motion about A with radius . Therefore,
| (18.73) |
But, and this trivial formula is the key, . Therefore, from Eqn. (18.72) and (18.73),
| (18.74) |
Taking dot product of both sides of the above equation with we get
Again taking the dot product of both sides of Eqn. (18.74) with we get
Answer:
SAMPLE 18.11
Spinning wheel on a rotating rod in 2-D . A rigid body OA is attached to a wheel that is massless except for three point masses P, Q, and R, placed symmetrically on the wheel. Each of the three masses is . The rod OA rotates about point O in the counterclockwise direction at a constant rate . The wheel rotates with respect to the arm about point A with angular acceleration and at the instant shown it has angular speed . Note that both and are given with respect to the arm.
Using a rotating frame attached to the rod and a coordinate system attached to the frame with origin at O, find
the velocity of the mass P and
the acceleration of the mass P.
Solution
Frame is attached to the rod. We choose a coordinate system in frame with its origin at O and, at the instant, aligned with the fixed coordinate system . We consider a point P′ momentarily coincident with point P but fixed in frame . Since P′ is fixed in , it rotates with with . To visualize the motion of P′ imagine a rigid rod from O to P′ (see Fig. 18.54). Now we can calculate the velocity and acceleration of point P′ as follows.
Velocity of point P:
Now we calculate the two terms separately:
Since the wheel rotates with angular speed with respect to the rod, an observer sitting in frame would see a circular motion of point P about point A. Therefore,
But at the instant of interest, So,
Therefore,
Answer:
Acceleration of point P: We can similarly find the acceleration of point P:
where
| acceleration of point P′ | ||||
| Coriolis acceleration | ||||
| acceleration of P relative to frame | ||||
The term encompasses three terms of the five term acceleration formula. The last line in the calculation of follows from the fact that at the instant of interest and .
Now adding the three parts of we get
Answer:
18.3.1 A bug walks on a straight line engraved on a rotating turntable (the bug’s path in the room is not a straight line). The line passes through the center of the turntable. The bug’s speed on this line is 1 (the bug’s absolute speed is not 1). The turntable rotates at a constant rate of 2 revolutions every in a positive sense about the -axis. Its surface is always in the -plane. At time , the engraved line is aligned with the -axis, the bug is at the origin and headed towards the positive direction.
What is the component of the bug’s velocity at ?
What is the component of the bug’s velocity at seconds?
What is the component of the bug’s acceleration at ?
What is the radius of curvature of the bug’s path at ?
18.3.2 Actual path of bug trying to walk a straight line. A straight line is inscribed on a horizontal turntable. The line goes through the center. Let be the angle of rotation of the turntable which spins at constant rate . A bug starts on the outside edge of the turntable of radius and walks towards the center, passes through it, and continues to the opposite edge of the turntable. The bug walks at a constant speed , as measured by how far her feet move per step, on the line inscribed on the table. Ignore gravity.
Picture. Make an accurate drawing of the bug’s path as seen in the room (which is not rotating with the turntable). In order to make this plot, you will need to assume values of and and initial values of and . You will need to write a parametric equation for the path in terms of variables that you can plot (probably and coordinates). You will also need to pick a range of times. Your plot should include the instant at which the bug walks through the origin. Make sure your and - axes are drawn to the same scale. A computer plot would be nice.
Calculate the radius of curvature of the bug’s path as it goes through the origin.
Accurately draw (say, on the computer) the osculating circle when the bug is at the origin on the picture you drew for (a) above.
Force. What is the force on the bug’s feet from the turntable when she starts her trip? Draw this force as an arrow on your picture of the bug’s path.
Force. What is the force on the bug’s feet when she is in the middle of the turntable? Draw this force as an arrow on your picture of the bug’s path.
18.3.3 A small bug is crawling on a straight line scratched on an old record. The scratch is a distance from the center of the turntable. The turntable is turning clockwise at a constant angular rate . The bug is walking, relative to the turntable, at a constant rate , straight along the scratch in the -direction. At the instant of interest, everything is aligned as shown in the figure. The bug has a mass gram.
What is the bug’s velocity?
What is the bug’s acceleration?
What is the sum of all forces acting on the bug?
Sketch the path of the bug in the neighborhood of its location at the time of interest (indicate the direction the bug is moving on this path).
18.3.4 Arm OC rotates with constant rate . Disc D of radius rotates about point C at constant rate measured with respect to the arm OC. What are the absolute velocity and acceleration of point P on the disc, and ? (To do this problem will require defining a moving frame of reference. More than one choice is possible.)
Answer: , .
Pick a suitable moving frame and do the problem.
Pick another suitable moving frame and redo the problem. Make sure the answers are the same.
18.3.5 For the configuration in problem 18.3 what is the absolute angular velocity of the disk, ? Answer: .
18.3.6 For the configuration in problem 18.3, taking to be the angular velocity of disk relative to the rod, what is the absolute angular acceleration of the disk, ? What is the absolute angular acceleration of the disk if and are not constant?
18.3.7 A turntable oscillates with displacement . The disc of the turntable rotates with angular speed and acceleration and . A small bug walks along line with velocity relative to the turntable. At the instant shown, the turntable is at its maximum amplitude , the line is currently aligned with the -axis, and the bug is passing through point on line . Point is a distance from the center of the turntable, point . Find the absolute acceleration of the bug, . Answer: .
18.3.8 A small kg toy train engine is going clockwise at a constant rate (relative to the track) of on a circular track of radius . The track itself is on a turntable that is rotating counter- clockwise at a constant rate of . The dimensions are as shown. At the instant of interest the train is pointing due south () and is at the center of the turntable.
What is the velocity of the train relative to the turntable ?
What is the absolute velocity of the train ?
What is the acceleration of the train relative to the turntable ?
What is the absolute acceleration of the train ?
What is the total force acting on the train?
Sketch the path of the train for one revolution of the turntable (surprise)?
18.3.9 A giant bug walks on a horizontal disk. An frame is attached to the disk.
The disk is rotating about the
axis (out of the paper) and simultaneously translating with respect to an inertial
frame plane. In each of the cases shown in the figure, determine the total
force acting on the bug. In each case, the dotted line is scratched on the disk and
is the path the bug follows walking in the direction of the arrow. The
location of the bug is marked with a dot. At the instants shown, the
coordinate system
shown is aligned with the inertial frame (which is not shown in
the pictures).
In this problem:
= the position of the center of the disk,
= the velocity of the center of the disk,
= the acceleration of the center of the disk,
the angular velocity of the disk (i.e., ),
= the angular acceleration of the disk,
= the speed of the bug traversing the dotted line (arc length on the
disk per unit time),
= , and
= mass of the bug.
18.3.10 Repeat questions (a)-(f) for the toy train in problem 18.3 going counter-clockwise at constant rate (relative to the track) of on the circular track.
18.3.11 A nostalgic honeybee alights on a phonograph turntable that is being carried by a carnival goer who is riding on a carousel. The situation is sketched below. The carousel has angular velocity of , which is increasing (accelerating) at ; the phonograph rotates at a constant 33 1/3 . The honeybee is at the outer edge of the phonograph record in the position shown in the figure; the radius of the record is 7 inches. Calculate the magnitude of the acceleration of the honeybee.
18.3.12 Consider a turntable on the back of a pick-up truck. A bug walks on a line on the turntable. The line may or may not be drawn through the center of the turntable. The truck may or may not be going at constant rate. The turntable may or may not be spinning, and, if it spins, it may or may not go at a constant rate. The bug may be anywhere on the line and may or may not be walking at a constant rate.
Draw a picture of the situation. Clearly define all variables you are going to use. What is the moving frame, and what is the reference point on that frame?
One term at a time. For each term in the five term acceleration formula, find a situation where all but one term is zero. Use the turntable as the moving frame, the bug as the particle of interest. Answer: One case: , where refers to the bug; the bug stands still on the turntable, the truck accelerates, and the turntable does not rotate.
Two terms at a time. (harder) How many situations can you find where a pair of terms is not zero but all other terms are zero? (Don’t try to do all 10 cases unless you really think this infamous formula is fun. Try at least one or two.) Answer: One case: , where is the turntable; the bug accelerates along a straight line on the turntable, the truck accelerates, and the turntable does not rotate.
Picking apart the five term acceleration formula.
An ideal mechanism or linkage is a collection of rigid objects constrained to move relative to the ground and each other by hinges but which still has some possible motion(s). People also use the word mechanism or linkage more loosely to include any collection of machine parts connected by any means.
The analysis of the kinematics of mechanisms is an important part of machine design. Mechanisms synthesis, coming up with a mechanism design which has desired motions, is obviously key in creative design, and nowadays in computer aided design.
Finally, the determination of the dynamics of a mechanism, how it will move and with what forces, is completely dependent on understanding the kinematics of the mechanism. The whole subject of mechanism kinematic analysis, although in some sense a subset of dynamics, is actually a huge and infinitely complex subject in itself
and also a useful subject in itself. Often kinematics is the central interest in machine design, and mechanics (force and acceleration) analysis is only carried out if something that shouldn’t do so, breaks or shakes. This section presents some of the basic ideas in kinematic analysis. The overall question in mechanism kinematic analysis is this:
Given a collection of parts and a description of how they are connected, in what ways can they move?
Without getting into the details of the motions yet, the first question to answer is simpler than finding the motions, but just counting them: In how many ways can the mechanism move?
The number of degrees of freedom (DOF) of a mechanism is the number of different ways it can move. More precisely
The number of degrees of freedom of a mechanism is the minimum number of configuration variables needed to describe all possible configurations of the mechanism.
The minimum number of configuration variables is a property of the mechanism. The choice of what these variables are, however is not unique to a given mechanism.
Example: A particle in a plane has 2 degrees of freedom.
The set of ‘configurations’ of a particle in a plane is the set of positions of the particle. This is fully described by its and coordinates. Thus . But the configuration is also determined by the particle’s polar coordinates and . And there are an infinite number of other pairs of numbers that could be used to describe the configurations (e.g., the and coordinates, the and coordinates with and , etc). The minimum number of configuration variables, 2, is unique, but the choice of variables is not.
For planar mechanisms one can often determine the number of degrees of freedom by the following formula
| (18.75) |
The formula starts with the number of ways one rigid object can move (2 translations and a rotation makes 3) and one particle can move (just 2 translations) and then subtracts the restrictions to the motion. In eqn. (18.75) the number of constraints is counted as the number of degrees of freedom restricted by the connections.
See fig. 18.64 for some standard idealized connections and their number of constraints (assuming they are already constrained to a plane).
2 for a pin joint: a pin joint restricts relative motion in 2 directions but still allows relative rotation. If three objects are connected at one pin then it counts as two pin joints and thus reductions in the number of degrees of freedom. There are 6 reductions for 4 objects connected at one pin, etc.
3 for a welded connection: a weld restricts relative translation in two directions as well as relative rotation (). So two parts that are welded together have degrees of freedom. That is, any collection of rigid objects welded together is the same as one rigid object. The word ‘weld’ is meant to include any collection of bolts, glue, string, rivets or bailing wire that prevents any relative motion.
1 for a sliding contact: the sliding contact restricts relative translation normal to the contact surfaces and allows translation tangent to the surfaces. Relative rotation is also allowed.
2 for a keyed sliding contact: allows relative translation in one direction but disallows translation in one direction as well as rotation.
1 for a massless link hinged at its ends to two objects: this keeps the distance between two points fixed which is one restriction (alternatively the bar adds 3 degrees of freedom and each hinge subtracts 2 ( degree of freedom).
2 for a rolling contact: relative slip is not allowed nor is interpenetration.
Figure 18.65 shows some examples of simple mechanisms and the number of degrees of freedom. In each case we look at eqn. (18.75): .
An object connected to ground by a hinge has 1 degree of freedom; the set of all possible configurations can be described by the angle of the object: .
An unconstrained object has 3 degrees of freedom; the set of all configurations can be described by the and coordinates of a reference point and by the rotation : .
A bead on a wire has 1 degree of freedom; its configuration is fully determined by the distance the bead has advanced along the wire relative to a reference mark: .
A statically determinate truss has 0 degrees of freedom; a statically determinate truss has no way to deform: . Note that the number of pin restrictions in our count here is more than twice the number of joints in the study of trusses in statics. In statics, we focussed on joints as restricted by bars. Here we look at bars as restricted by joints. A given joint counts 2, 4 or 6 pin restrictions depending on whether it connects 2, 3 or 4 bars. Thus the 11 bar truss shown has 7 joints (by truss-analysis counting) but, by the counting here it has pin restrictions plus 3 ground restrictions for 33 restrictions in all. Without restrictions the 11 bars had 33 degrees of freedom. So, because 33-33 = 0, and none of the restrictions is redundant, the truss has zero degrees of freedom.
A rolling wheel has 1 degree of freedom; its configuration is fully determined either by the net angle it has rolled or by the coordinate of its center: , .
A double pendulum or two-link robot arm has 2 degrees of freedom; its configuration is determined by the net rotations of its two links (or by the rotation of the first link and the relative rotation of the second link): .
A cart with two rolling wheels has 1 degree of freedom; given the position of the cart, say , the rotation of the wheels is determined. Similarly for a bicycle. For this model of a bicycle we are neglecting the steering degree of freedom and also the freedoms of the crank and pedals. For a vehicle with 2 wheels in 2D we have: (2 hinges and 2 rolling contacts) .
A “four” bar linkage has 1 degree of freedom; the angle of any one of the bars determines the angles of the others: (there are 4 pin joints between the bars, one pin joint to ground and one roller connection to the ground).
A slider crank has 1 degree of freedom; the rotation of the crank determines the configuration of the system (there are 3 pins and one keyed connection) .
An ideal gear train (with all gears pinned to ground) has 1 degree of freedom; the amount of rotation of any one gear determines the rotation of all of the gears: In this case the counting formula is wrong. Say there are 2 gears, then (two pins and one rolling contact) . The rolling constraint prevents interpenetration, but this was already prevented by the hinges at the center of the gears. The constraints are redundant and the system has more degrees of freedom than eqn. (18.75) indicates.
A redundant swing with one horizontal bar suspended by 3 parallel struts has 1 degree of freedom; the angle of one upright link determines the full configuration of the mechanism. The counting formula is again wrong: underestimates the number of degrees of freedom because the constraints are redundant.
A 2-D 10-link model of a person with one foot on the ground has 10 degrees of freedom; the angles of the 10 links (4 arm links, 4 leg links, a body and a head) determine the full configuration of the mechanism. There are no redundant constraints so the counting formula also works using :
| (18.76) | ||||
| (18.77) |
What are the 10 hinges? One at the ground, 1 at each knee and elbow, 2 at the hip and 3 at the shoulders (one hinge for each object linked to the trunk).
Once we know the number of degrees of freedom of a system it is often useful to settle on one set of configuration variables. In this book will be 1, 2 or at most 3. Thus we pick 1,2 or 3 variables.
Example: Straight line motion
Chapter 14 on straight line motion was mostly about one-degree-of-freedom systems (). These systems could all be characterized by the single configuration variable , the displacement along the line of a reference point on the object relative to a reference point on the ground. All the positions,velocities and accelerations of all points in the system could be found in terms of , and (in fact all points had and ).
Example: Circular motion about a fixed axis
In chapters 15 and 16 we were almost entirely focussed on systems with one degree of freedom well characterized by the one configuration variable, the rotation angle . For such motions positions, velocities, and accelerations of all points were determined by the initial positions of the points and by equations which you know well by now.
For more general motions we almost always take inspiration from the two examples above. We use the translation of a conspicuous point, or we use the rotation of a conspicuous object for a configuration variable. And more of the same if the system has more than 1 degree of freedom. The natural choice of configuration variables for some simple mechanisms is given in the text discussing fig. 18.65.
Often our main kinematic task is to express the full configuration of the system as well as all the velocities and accelerations of all its parts in terms of the positions of the parts, the configuration variables, and their first and second time derivatives.
One last simple kinematic fact is needed before we can plug and chug with the theory we have so far and apply it to general kinematic mechanisms. It concerns the addition of rotations and rotation rates. The following example basically tells the whole story
Example: Double pendulum and the addition of rotation rates
The commonly used configuration variables for the double pendulum shown in fig. 18.66 are and . To actually know the configuration of the system obviously we need to know which is given by
[Aside: One reason for choosing instead of as a configuration variable is that if one was measuring or controlling the second link, say as a robotic arm, the angle can be measured more easily than . Also, it turns out (in hindsight) that the differential equations of motion are slightly simpler using instead of .]
Looking at the bars as being glued to reference frames ( for the fixed frame, for bar AB, and for bar BC), the above example shows that
| (18.78) | |||||
| which is often written with the simple notation | |||||
Which can only be given strict meaning by the more elaborate eqn. (18.78) above it.
Example: Double pendulum (see previous example)
Take , and and eqn. (18.78) is self evident from the addition of angles.
One approach to mechanisms is to do what one can with high-school geometry and trigonometry, the laws of sines (see page 1.5), and so on.
Example: Rod on step using geometry and trigonometry
One end of a rod slides on the ground. The other end slides on a corner at A (see fig. 18.67). Given that we can find as follows:
As the above calculation shows, this problem doesn’t need the heavy machinery of our moving-frame vector methods. But it provides an instructive example.
Example: Rod on step using moving-frame methods (see previous example)
We look at point A and note that we can think of it as a fixed point in the fixed frame and also as a point that is moving relative to the translating and rotating frame . We evaluate its velocity both ways.
| (eqn. (18.67)) | |||||
| ( ) | |||||
| ( ) | |||||
as we had before. The key equation was the ‘three term velocity formula’ eqn. (18.67) on page 18.67 and the observation that relative to frame point A slides along the rod. Note that we never had to explicitly use the rotating coordinates associated with frame to do this calculation.
You should understand the examples above, and the needed background material, before going on to the following examples.
Example: Slider Crank using geometry and trigonometry
The slider crank mechanism (fig. 18.68) was briefly introduced in the context of statics where its forces could be analyzed assuming inertial terms were negligible (see 7.3). But it is a commonly used mechanism (e.g., in every car) and its motions are of central interest. The angle is the most natural configuration variable for this system. One would like to know the position, velocity and acceleration of the slider ( and ) in terms of and .
| (18.80) | |||||
The positive corresponds to C being to the right of 0. The negative corresponds to point C being to the left of 0. The mechanism just doesn’t work for a full revolution of the link 0A if as you can see from the picture or from that the above giving imaginary values for near -1, near 1, and near .
To get the velocity of point C we just take the derivative of eqn. (18.80) above.
To get the acceleration we differentiate once again. For simplicity let’s assume the crank rotates at constant rate, so is a constant and . Cranking out the derivative of eqn. (18.4), so to speak, we get
So we now know the position, velocity and acceleration of point C in terms of and . You should commit the solution eqn. (18.4) to memory. Just kidding.
Plots of , and from these equations are shown in fig. 18.69ab for two different extremes of slider crank design: one with a very long connecting rod that gives sinusoidal motion, and one with a connecting rod just barely long enough to prevent locking that gives intermittent motion.
Unlike some more complex mechanisms, the slider crank is solvable in that one can write a formula for the position of any point of interest in terms of the single configuration variable . For more complex mechanisms this may not be possible. Further, even if possible the above example shows that the differentiation required to find velocity and acceleration can lead to a bit of a mess.
A different approach is to assume that at some value of the configuration variable ( for the slider crank) that the full configuration of the system is known. That is, that the locations of all points are known. Then we can use our vector methods to find velocities and accelerations of all points of interest.
Example: Slider crank using vector methods (see previous example)
Take the slider crank of fig. 18.68 to be in some known configuration. We now try to find the velocity and acceleration of point C in terms of the positions of the points 0, B, and C as well as and .
The basic approach is to write true things, and then solve for unknowns. First work on velocities. The basic idea is to look at the closure condition. That is, the velocity of point C as calculated by working down the linkage from 0 to A to C has to be consistent with the velocity of C as calculated in the fixed frame.
| (18.83) |
eqn. (18.83) is a 2-D vector equation in the 2 unknown scalars and . It could be solved as a pair of equations, or solved directly by first dotting both sides with to find and dotting both sides with to find . These yield
where everything on the right of these equations is assumed known. Without grinding out the vector products in terms, say, of components, we can just use that we know and .
We proceed to find the accelerations by similar means, assuming is a constant so :
| (18.84) |
eqn. (18.84) is a 2-D vector equation in the two scalar unknowns and . We can set this up as two equations in two unknowns. Or we can solve for directly by dotting both sides with and we can solve for directly by crossing both sides with or by dotting with a vector perpendicular to (e.g., )
Although we have presented an algorithm rather than a formula, we have found the velocity and acceleration of C without writing any large equations of the type needed in the previous example. The shortcoming is that this method depends on knowing the full configuration at the time of interest.
Example: Four bar linkage using geometry and trigonometry
fig. 18.70 shows a “four-bar linkage”. Please see 7.3 for an introduction to 4-bar linkages in the context of statics. Four bar linkages are solvable in the sense that one can write equations for the positions of any point of interest in terms of the single configuration variable marked in fig. 18.70. But the formulas are really a mess. And the first and second time derivatives are an unbelievable mess.
The four-bar linkage is about as complex a system as can be solved in this sense, and it is probably too-complex for this solution to be useful in the kinematic analysis of accelerations.
For complex mechanisms one is often stuck using vector methods, like we are for practical purposes stuck with the 4-bar linkage. But the vector methods based on the current configuration are not crippled by complexity.
Example: Four-bar linkage using relative velocities and accelerations
Assuming the configuration is known (i.e., that , , and are known), we can proceed with the 4-bar linkage just as we did for the slider crank. We enforce closure by picking a point and thinking about its velocity two different ways
We could pick any point, say C. From the fixed frame we know that the velocity of C is zero. Working around the linkage, link by link, we know it is the sum of relative velocities as
which is equivalent to two scalar equations in the two unknowns and . This equation can be solved directly for by taking the dot product of both sides with a vector perpendicular to (such as or ) and for by taking the dot product of both sides with for a vector perpendicular to (such as or ) to get
The dot product with is used to get a scalar on the top and bottom of the fraction, both vectors are already only in the direction. Now that and are known the velocity of any point on the mechanism is known. For example
The angular accelerations of the two links are found by the same method. For simplicity let’s assume that the driving crank 0A spins at constant rate so . Looking at the acceleration of point C two ways we have
Because and are already known, this is one equation in the two unknowns and . They can be solved for by taking the dot product of both sides with and for by taking the dot product of both sides with .
At this point you know , , , , , , , and and can thus calculate the position, velocity and acceleration of any point in the mechanism.
SAMPLE 18.12
Two rods connected in a one-DOF mechanism. A mechanism consists of two rods AB and CD connected together at P with a collar pinned to AB but free to slide on CD. Rod AB is driven with and . At the instant shown, and . The length of rod AB is . At the instant shown,
Find the angular velocity and angular acceleration of rod CD.
Find the velocity and acceleration of the collar with respect to rod CD.
Solution Here, we are interested in instantaneous kinematics of this mechanism. Since point P is on rod AB as well as on rod CD, its velocity and acceleration can be expressed in terms of the angular motion of rod AB or that of rod CD. Let us consider rod AB first. Let and be basis vectors attached to rod AB that rotate with the rod. Since P is fixed on rod AB, it executes simple circular motion about A with and where and , respectively. Then
| (18.85) | |||||
| (18.86) |
Now let us consider rod CD and express the velocity and acceleration of point P in terms of motion of rod CD. Let the angular velocity and angular acceleration of rod CD be and , respectively. Let and be basis vectors attached to rod CD. Let the instantaneous position of point P on CD be . Since the collar can slide along CD, we can write the velocity and acceleration of point P as
| (18.87) | |||||
| (18.88) |
Thus, to find all kinematic quantities of interest, all we need now is to figure out
a few dot products between the two sets of basis vectors. This is easily done by
writing out , and .
††margin:
We have,
,
,
,
.
Therefore,
Substituting the dot products in the expressions for , and
we get
Substituting the given values of , and , we get
Answer:
SAMPLE 18.13
Kinematics of a Link rod in a one-DOF mechanism. In machines we often encounter mechanisms and links in which the ends of a link or a rod are constrained to move on a specified geometric path. A simplified typical link AB is shown in Fig. 18.74.
Link AB is a uniform rigid rod of length . End A of the rod is attached to a collar which slides on a horizontal track. End B of the rod is attached to a uniform disk of radius which rotates about its center O. At the instant shown, when and , end A is observed to move at to the left.
Find the angular velocity of the rod.
Find the angular velocity of the disk.
Find the velocity of the center-of-mass of the rod.
Solution Let the angular velocities of the rod and the disk be and respectively, where and are unknowns. We are given where .
Point B is on the rod as well as the disk. Hence, the velocity of point B can be found by considering either the motion of the rod or the disk. Considering the motion of the rod we write,
| (18.89) | |||||
Now considering the motion of the disk (and noting that ), we write,
By equating the and components of the above equation we get
| (18.91) | |||||
| (18.92) |
Answer:
From eqn. (18.92)
Thus
Answer:
[At this point, it is a good idea to check our algebra by substituting the values of and in equations (18.89) and (18.90) to calculate .] ††margin: Substituting in Eqn. (18.89) and plugging in the given values of other variables we get Similarly, substituting in Eqn. (18.89) and plugging in the other given values we get which checks with the found above.
Now we can calculate the velocity of the center-of-mass of the rod by considering either point A or point B as a reference:
Answer:
We could easily check our calculation by taking point B as a reference and writing
By plugging in appropriate values we get, of course, the same value as above.
Comment: We used the standard basis vectors and for all our vector calculations in this sample. We can shorten these calculations by choosing other appropriate basis vectors as we show in the following samples.
SAMPLE 18.14
A two-DOF mechanism. A two degree-of-freedom mechanism made of three rods and two sliders is shown in the figure. At the instant shown, the crank AB is rotating with angular velocity and angular acceleration . At the same instant, the collar at end C of the link rod CD is sliding on the vertical rod with velocity and acceleration . Find the angular velocity and angular acceleration of the link rod CD.
Solution Once again, we are interested in instantaneous kinematics — we wish to find the angular velocity and acceleration of rod CD at the given instant. This problem is just like the previous sample problem except that end C of the link rod CD is not fixed but free to slide on the vertical bar. But the velocity and acceleration of point C is given; so it is exactly like the previous sample (there, the velocity and acceleration of point C was identically zero). So, we adopt the same line of attack. We figure out the velocity and acceleration of point B using the kinematics of rod AB. We then write the velocity and acceleration of the same point using the kinematics of rod CD (this will involve the unknown angular velocity and acceleration of CD that we are interested in). Equate the two and solve for the unknowns we are interested in.
Let the angular velocity and acceleration of rod CD be and , respectively. Let and be base vectors rotating with rod AB, and and be the base vectors rotating with rod CD (see fig. 18.76). Considering rod AB, we have
| (18.93) | |||||
| (18.94) |
Considering rod CD, we have
| (18.95) | |||||
| (18.96) |
Now equating eqn. (18.93) and (18.95), and dotting both sides with and , we get
| (18.97) | |||||
| (18.98) |
where the dot products among the basis vectors are easily found from either their geometry (see fig. 18.77) or from their component representation (see previous sample). Following exactly the same procedure, we get, from eqn. (18.94) and 18.96,
| (18.99) | |||||
| (18.100) |
Now, note that although we are only interested in finding and . So, we only need eqn. (18.98) and eqn. (18.100). But, eqn. (18.100) requires on the right hand side and, therefore, we do need eqn. (18.97). We can, however, happily ignore eqn. (18.99).
Now, to find the numerical values of and , we need to find and in addition to all other given values. Consider triangle ABC in fig. 18.76. Using the law of sines (), we get and . Now, substituting all known numerical values in eqns. (18.97), (18.98), and (18.100), we get
Answer:
18.4.1 Slider crank kinematics (No FBD required!). 2-D . Assume are given. The crank mechanism parts move on the plane with the direction being along the piston. Vectors should be expressed in terms of and components.
What is the angular velocity of the crank OA? Answer: .
What is the angular acceleration of the crank OA? Answer: .
What is the velocity of point A? Answer: .
What is the acceleration of point A?
Answer: .
What is the angular velocity of the connecting rod AB? [Geometry fact: ] Answer: .
For what values of is the angular velocity of the connecting rod AB equal to zero (assume )? (you need not answer part (e) correctly to answer this question correctly.) Answer: , .
18.4.2 Slider-Crank. Consider a slider-crank mechanism. Given , , , , and , can you find the velocity and acceleration of B? There are many ways to do this problem.
18.4.3 The crank AB with length in the crank mechanism shown rotates at a constant rate counter-clockwise. The initial angle of rotation is at . The connecting rod BC has a length of .
What is the velocity of point B at the end of the crank when ?
What is the velocity of point C at the end of the connecting rod when ?
What is the angular velocity of the connecting rod BC when ?
What is the angular velocity of the connecting rod BC as a function of time?
18.4.4 The two rods AB and DE, connected together through a collar C, rotate in the vertical plane. The collar C is pinned to the rod AB but is free to slide on the frictionless rod DE. At the instant shown, rod AB is rotating clockwise with angular speed and angular acceleration . Find the angular velocity of rod DE. Answer:
18.4.5 Reconsider problem 18.81. The two rods AB and DE, connected together through a collar C, rotate in the vertical plane. The collar C is pinned to the rod AB but is free to slide on the frictionless rod DE. At the instant shown, rod AB is rotating clockwise with angular speed and angular acceleration . Find the angular acceleration of rod DE.
Answer:
18.4.6 Collar is constrained to slide on a horizontal rod to the right at constant speed . It is connected by a pin joint to one end of a rigid bar with length which makes an angle with the horizontal at the instant of interest. A second collar connected by a pin joint to the other end of the rigid bar slides on a vertical rod.
Find the velocity of point . Answer in terms of , , , and .
Find the angular speed of the rod? Answer in terms of , , , and .
18.4.7 A bar of length , body , connects sliders A and B on an L-shaped frame, body , which itself is rotating at constant speed about an axle perpendicular to the plane of the figure through the point and relative to a fixed frame , . At the instant shown, body is aligned with the axes, slider A is from point , slider B is from .
The speed and acceleration of slider B relative to the frame are and , respectively.
Determine:
The absolute velocity of the slider A, , and
The absolute angular velocity of bar AB, body , .
18.4.8 The link is supported by a wheel at D and its end A is constrained so that it only has horizontal velocity. No slipping occurs between the wheel and the link. The wheel has an angular velocity and radius . The distance .
Given: , , and . .
Determine the angular velocity of the link AB and velocity of the point A.
18.4.9 A solid cylinder of radius and mass rolls without slip along the ground. A thin rod of mass and length is attached by a frictionless pin P to the cylinder’s rim and its right end is dragged at a constant speed along the (frictionless) horizontal ground.
For the position shown (where P is directly above the contact point), find P’s velocity and the rod’s angular velocity .
Find P’s acceleration and the angular acceleration of the rod.
18.4.10 The slotted link CB is driven in an oscillatory motion by the link ED which rotates about D with constant angular velocity . The pin P is attached to ED at fixed radius and engages the slot on CB as shown. Find the angular velocity and acceleration and of CB when .
18.4.11 The rod of radius shown has a constant angular velocity of counterclockwise. Knowing that rod AD is long and distance , determine the acceleration of collar D when .
In this section we consider three types of problems where the kinematics involves solution of differential equations. In most cases this means computer solution is involved for this type of problem. Here are the three problem types:
I. Closed kinematic chains. The main simple example is a 4-bar linkage with one bar grounded. This system has one degree of freedom, but it is difficult to directly calculate the positions, velocities and accelerations of all points in terms of one variable. Instead, the constraint that the linkage is closed is sometimes most-easily expressed as a differential equation.
II. Rolling contact with not-round objects. For non-round rollers and cams solving for configuration, velocity and acceleration can sometimes be best done with integration. A side benefit from studying this topic is the observation that all non-translational motions are equivalent to rolling of some kind.
III. Contact with ideal wheels and skates, looking down. Cars, tricycles, trailers, grocery carts and sleighs have wheels and have dynamics that is sometimes well characterized by planar analysis, where the plane is the horizontal plane. In this view the simple model of a wheel is as something that prevents sideways motion but allows motion in the direction of travel (like for some of the trike and car problems in 1-D constrained motion). Such problems are called non-holonomic (see box 18.5 on page 18.5).
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 the 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 do 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: configuration variables
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.
However, given a configuration, the constraint on the rates and accelerations is relatively easy to express, always yielding linear equations.
Box 18.3 Skates, wheels and non-holonomic constraints
Of the words in this book “non-holonomic” is probably the most obscure. This is because the subject of mechanics was mostly stolen from engineers by physicists about 100 years ago. And physicists, the authors of most introductions to mechanics, had no use for non-holonomic mechanics as it wasn’t useful for the development of quantum mechanics. So many people are unaware of the word, the subject or its utility.
In two dimensions the word non-holonomic in effect means the mechanics of objects constrained by ideal skates or massless ideal wheels. Often these non-holonomic constraints are described as “non integrable”. Literally, the word non-holonomic means “not whole”. But in what sense is a rolling ideal wheel “non-integrable” or less “whole” than anything else?
A constrained rigid body. Consider a rigid body that is free to slide on a plane. It has three degrees of freedom described by and , all measured relative to a fixed reference frame . Point C on the body has relative position where and are constants. Now let’s constrain the body at point C one of these two different ways:
a) Pin the body to the ground with an ideal hinge at point C. This keeps point C from moving but allows the body to rotate (holonomic).
b) Put an ideal wheel or skate under the body at C that prevents sliding sideways to the skate but allows point C to move parallel to the skate and also allows rotation about the skate (non-holonomic).
Pin Constraint. In the first case, for the pin, we could describe the constraint with the phrase ‘point C on the body can have no velocity’ and then write and calculate:
| (18.101) | |||||
The last two equations are two differential equations in the three variables and . They are “integrable” in the sense that they are equivalent to
| (18.102) |
where and are integration constants that need to be set by the starting configuration. Solving for and in terms of :
Here we have derived the obvious, that a pinned body has one independent configuration variable , but we did so starting with a vector expression of constraint in terms of velocities (eqn. (18.101)). Then we wrote the constraint as two scalar constraints on the derivatives of configuration variables and then “integrated” them to write constraints on the configuration variables, finally eliminating two of the configuration variables.
Skate constraint. Now consider the same body constrained by a skate or ideal wheel at C instead of a pin. The skate is aligned with the so point C can only move in the direction. The body is still free to rotate about the point C (to steer). Thus,
| (18.103) | |||||
As for the hinge where we found 2 constant functions, we might want to find the function that satisfies the differential equation above, namely
| (18.104) |
Another math nightmare. How do we find this ? You can’t find one. This is the crux of the matter. Neither your calculus professor nor Ramanujan could find one either. No computer can find one, or even a numerical approximation of a solution. There is no function that solves eqn. (18.104). The solution fundamentally does not exist. That is why we say the skate/wheel constraint eqn. (18.103) is “non-integrable”.
Parallel parking We can use physical reasoning to show that no function can solve eqn. (18.104). If such an did exist it would mean that only the set of configurations with position and angles consistent with would be allowed by the skate constraint (assuming depends nontrivially on at least one of the variables). This means there would be some angles and positions that the body couldn’t get to. Remember, we are not doing mechanics, just kinematics. So we can see what configurations are geometrically allowed while still respecting the constraint. The simple observation that motivates the answer is this:
Even though the skate constrains to not have a sideways component, point C can get to a point that is straight sideways.
How? Like a car can move sideways into a parking space without skidding sideways; by parallel parking. More generally, the body can get to any position and any orientation by the following moves. First rotate the body so the skate aims to its new goal. Then slide the skate to its new goal. And finally rotate the body to its new desired orientation.
Thus, the skate constraint does not disallow any configurations! Yet the constraint does disallow some velocities (the skate can’t go sideways). In this way, the skate constraint is not “whole”. It constrains velocities without constraining configurations.
Counting degrees of freedom. How many degrees of freedom does a body with a skate constraint have? There are two different answers. By counting possible configurations there are three degrees of freedom (it takes three variables to describe all possible configurations). But at any configuration the velocity can be described by 2 numbers ( and ). Whenever the number of configuration degrees of freedom is greater than the number of velocity degrees of freedom (for example, 3¿2) there are non-holonomic constraints.
One might like more examples. But besides artificial mathematical ones, there are none. The only smooth non-holonomic constraint in 2D mechanics is the ideal skate or wheel.
Example: Four bar linkage: configuration rates
If you write the relative velocities of the ends of the bars in terms of configuration rates and and then write the chain closure equation you get a linear equation in the rates. Likewise if you write the closure condition in terms of acceleration. The coefficients in these equations are likely to be complex functions of the configuration, so integrating these equations requires numerics. But the constraint is linear.
Thus, as shown in the last sample of Sect. 10.4, one way to calculate the evolving configurations of a closed chain is to integrate the velocity relations numerically.
When two objects roll on each other they maintain contact and do not slip relative to each other. That is to say rolling of one rigid curve on another means:
The instantaneous relative motion of with respect to is a rotation about the contact point at the common tangent C, and
The sequence of points C moves the same distance on both curves.
For simplicity let’s take to be a curve fixed in space on which rigid curve rolls. Take a reference point of interest fixed on body to be O’. So,
| (18.105) |
and is, say, the rotation of a axis fixed in relative to a axis fixed in . If we use the rotation of body as our configuration variable, we now know how to find the velocity of all points in terms of their positions and the rotation rate. Thus we can find the rate of change of the configuration. To proceed as time progresses we also need to know how the position of point C evolves. Not the material point C on either body, but the location of mutual contact.
If we assume that both curves are parameterized by arc-length going counter-clock wise, if we take curvature as positive if directed towards the interior of each curve’s body, then the condition of maintaining contact requires that
and is the advance along curve . To maintain tangency, the angles must be maintained so
Altogether this gives
To find the acceleration of material point O’ on we differentiate eqn. (18.105) with respect to time:
where is normal to the curves and directed towards the interior of . Thus the acceleration of all points on is the same as if the body were pinned at C plus an acceleration due to rolling. This rolling acceleration is small if either of the bodies is sharp (has very large ) and large if the bodies are nearly conformal.
As one rigid body moves arbitrarily on a plane with some non-zero rotation rate we can find a point at relative position where . That is a place where the velocity due to rotation about O’ exactly cancels the velocity of O’.
Crossing both sides with and using that we get
as the point “on” the body that has no velocity. This point does not literally have to be on the body, rather it is fixed to the reference frame defined by the body.
As motion progresses a sequence of such points C is traced on the ground. Similarly a sequence of points is traced on the body. These two sequences are called the space curve and the body curve (or “polohodie” and “herpolhodie” in older books). The motion of body is thus a rolling of the body curve on the space curve.
As a machine designer this means you can generate any desired motion by rolling of appropriate shapes. Move the object in the desired manner, draw the space curve and body curve, make parts with those shapes, and the desired motion occurs by a rolling of those shapes.
If we look down on an ideal skate or wheel at point C on a rigid body and assume that the skate is oriented with the positive axis at point C on the body then we know that
and hence the velocity of any point G on the body of interest is
The acceleration is found by differentiating this expression as
It is interesting to note that the Coriolis-like term does not have the usual factor of 2 one encounters in holonomic problems. To find the trajectory of the point C, say, one needs to integrate the velocity like this:
SAMPLE 18.15
Kinematics of a four bar linkage. A four bar linkage ABCD is shown in the figure (fourth bar is the ground AD) at some instant . The driving bar AB rotates with angular velocity . Find the angular velocities of rods BC and CD as a function of . How can you solve for the positions of the bars at any if the initial configuration is as shown in the figure?
Solution Let the angles that rods AB, BC, and CD make with the horizontal (-axis) be , and , respectively. Then, we can write and . We have to find and .
Note that the motion of point B is a simple circular motion about point A with given angular velocity . Thus, the velocity of point B is known. Now, we can find the velocity of point C two ways: (i) by considering rod CB: , and (ii) by considering rod CD: Either way the velocity must be the same. Thus, we have a 2-D vector equation with two unknowns and . We can get two independent scalar equations from the vector equation and thus we can solve for the desired unknowns.
Let us use the rotating base vectors with rods AB, BC, and CD, respectively. Note that these base vectors are basically the pairs; we use just for the sake of easy subscripting. Now,
| (18.106) | |||||
| (18.107) |
| (18.108) |
Dotting eqn. (18.108) with (to eliminate term), we get
| (18.109) |
Similarly, dotting eqn. (18.108) with (to eliminate term), we get
| (18.110) |
We are practically done at this point with the kinematics — we have found and as functions of . The various dot products are just geometry and vector algebra. To write them explicitly, we note that , , etc. Thus,
Answer:
Note that the expressions for and are coupled, nonlinear, first order ordinary differential equations. To be able to find and , we need to integrate , and . Here, we set up these differential equations for numerical integration. Although, we can use any given (e.g., or or whatever), for definiteness in our numerical integration, let us take a constant , that is, let (say). So, our equations are,
and the initial conditions are
Here is a pseudocode that we use to integrate these equations numerically for a period of seconds (one complete revolution of AB).
ODEs = {thetadot = C,
betadot = (l1/l2)*sin(phi-theta)/sin(beta-phi)*C,
phidot = (l1/l3)*sin(beta-theta)/sin(beta-phi)*C}
IC = {theta(0) = pi/2, beta(0) = pi/4, phi(0) = 3*pi/4}
Set C=10, l1=.4, l2=.4*sqrt(2), l3=.8*sqrt(2)
Solve ODEs with IC for t=0 to t=pi/5
After we get the angles, we can compute the coordinates of points B and C at each instant and plot the mechanism at those instants. Plots thus obtained from our numerical solution are shown in fig. 18.91 where the configuration of the mechanism is shown at 9 equally spaced times between to .
18.5.1 Double pendulum. The double pendulum shown is made up of two uniform bars, each of length and mass . At the instant shown, , , , and are known. For the instant shown, answer the following questions in terms of , , , , and .
What is the absolute velocity of point ?
What is the velocity of point relative to point ?
What is the absolute velocity of point ?
18.5.2 Double pendulum, Again. For the double pendulum in problem 18.5, and are also known at the instant shown. For the instant shown, answer the following questions in terms of , , , , , , and .
What is the absolute acceleration of point ?
What is the acceleration of point relative to point ?
What is the absolute acceleration of point ?