Here we extend the idea of circular motion from particles to objects. A particle can go in circles. On a rigid object in 2D, all pairs of particles go in circles around each other. The key theoretical idea is of rigid-object rotation. The primary applications are pendulums, gear trains, and rotationally accelerating motors or brakes.

We covered the special case of circular motion of a particle in the previous chapter. Now we consider the kinematics and mechanics of rotating rigid objects in 2D. The key new conceptual idea is that of rotation of a rigid object.
For the systems in this chapter, for every system we show in a free-body diagram we have, as always,
Because you already know how to work with forces and moments (the left sides of the top two equations), the primary new skill in this chapter is the evaluation of , and for a rotating rigid object. That is, you need to understand the position, velocity and acceleration of points moving in circles. The rest of the skills used are universal, for example solving algebraic or differential equations, plotting, etc.
When two parts are glued together or attached by welding, gluing, several tight screws, bolts, rivets or the like we call the connection a ‘rigid attachment’. And, for the purposes of mechanics analysis, the two connected parts make up one bigger object. But most machines have various parts that are connected to each other, but not welded to each other. The most common such non-rigid attachment in engineering is a hinge. In 2D,
a hinge attachment between two objects keeps two points, one from each object, on top of each other while freely allowing relative rotation of the two objects about the hinge point.
In 3D a hinge keeps two lines, one on each body, coincident and allows relative rotation about that line. The common line, or in 2D the line orthogonal to the plane through the points, is called the hinge, the hinge axis, or the axis of relative rotation.
One example of a hinge is a car axle which allows rotation of a wheel relative to the car suspension. The hinge axis is the wheel axle
.
Physically, hinges are made various ways, sometimes by poking a cylindrical pin through the two objects and sometimes with ball bearings (see box 5.3 on page 5.3). So hinges are also called pin connections or bearings (fig. 16.3). In this chapter we limit our attention to a simple use of a hinge: one rigid part is hinged to a second fixed part that doesn’t move at all. Such a non-moving part can be thought of as connected to, or an extension of, the ground or ‘fixed frame’. In the simple case of interest in this chapter, one point on the moving part does not move and the rest of the part rotates about that point.
For definiteness and simplicity let’s call the hinge location 0 and the hinge axis through 0 the axis.
One function of the hinge is to make the part’s only possible motion to be rotation about O. Thus to understand the dynamics of a hinged part we need to understand the position, velocity and acceleration of points on a rigid object which rotates. This whole section is about the kinematics (the geometry of motion) for this rotation. We will measure the amount of rotation by the angle , and the rate of change of by the angular velocity (‘omega’), and the rate of change of this angular velocity by the angular acceleration (‘alpha’).
As simple as this topic seems at first glance, you should pay close attention to the meanings and uses of these quantities. The rest of the book completely depends on the material in this section.
Box 16.1 Rotation is uniquely defined for a rigid object (2D)
Most people will find it self-evident that, starting with a rigid object at a reference orientation, all lines marked on the object rotate by the same angle . Here, for the doubting, we demonstrate this fact.
A rigid object is defined this way:
For every pair of material points A and B on a rigid object the distance AB between them does not change as the object moves.
In particular, when a rigid object rotates all distances between all pairs of points are preserved. Thus, by the “side-side-side” similar triangle theorem of elementary geometry, all relative angles between marked line segments are preserved by the rotation. For example, for a triangle ABC the angle at B is constant as the object rotates. Now consider any pair of line segments on the object.
[By ‘on’ the object we don’t mean a projected image drawn with a light pen that can move around on the object relative to the atoms. Rather, by ‘on the object’ we mean something defined by a particular set of atoms that make up the object.]
If the segments do not cross we can extend them to a point of intersection B. Such a pair of intersecting lines is shown here before and after rotation.Initially BA makes an angle with a horizontal reference line. BC then makes an angle of . After rotation we measure the angle to the line BD. BA now makes an angle of . By the addition of angles in the rotated configuration line BC now makes an angle of which makes an increase by of the angle made by BC with the horizontal reference line. So both BA and BC rotate by the same angle .
We could use one of these two lines and compare it with an arbitrary third line through B and show that the third line also has equal rotation, and then a fourth, and so on. So all lines on the object through a point B rotate by the same angle . The demonstration for a pair of parallel lines, one of them through B, is easy, they stay parallel so always make a common angle with any reference line.
Because any line on the object either goes through B or is parallel to a line through B, all lines marked on a rigid object rotate by the same angle .
The rotation of a rigid object in 2D is thus unambiguously defined as the angle through which all lines on the object rotate.
We start by imagining the object in a distinguished configuration which we call the reference configuration, reference state or reference position. For example we could take the left figure in fig. 16.3 as the reference configuration. If possible it’s usually best to pick the reference state to be one in which a prominent feature of the object is aligned with the or axes. The reference state may or may not be the start of the motion of interest. Even if not, we measure an object’s rotation by the change, relative to the reference state, in the counterclockwise angle of a reference line marked in the object relative to a fixed line outside. Which reference line? Fortunately,
All real or imagined lines marked on a rotating rigid object rotate by the same angle, the rotation angle, . (See box 16.1).
In three dimensions things are more complicated. General rotation of a rigid object is then represented not with a single angle , but rather with 3 angles, or with a unit vector and an angle, or a 33 matrix. So we wait to discuss which of the 2D ideas here generalize to 3D and which do not.
We pick two orthogonal lines on the rotating object and give them distinguished status as object-fixed (or body-fixed) rotating coordinate axes and . Think of these axes as coordinate axes on a piece of graph paper that is glued to the object (see fig. 16.5a&b). It’s easiest if we start by assuming that the axes have the same origin 0 as the axes and are parallel with the fixed axes when the object is in the reference configuration (when ).
These rotating coordinate axes, and , have associated rotating base vectors and (fig. 16.6 and 16.7). So is always in the direction and always in the direction. We will use these rotating coordinates and base vectors to keep track of some particle of interest P that is ‘glued’ to the object. To start, note that particle P which is glued to the object has and coordinates that don’t change as the rotation progresses.
Example: A particle on the axis
If a particle P is fixed on the -axis at position , then we have.
for all time, even as the object rotates.
The position vector of a point P fixed to a rigid object hinged at O remains, as the rotation progresses,
| (16.1) |
with and both constant. These rotating coordinate system components, , are sometimes written as .
You will see that much of the math for rotating coordinates is reminiscent of that for polar coordinates. However, the spirit is a bit different. In polar coordinates the axis was picked to track a particular particle of interest. Here we pick axes that rotate with an extended object and use that one set of axes to track any and all particles of interest.
Note, even though neither nor change as changes, the point P they describe moves, in circles actually. How can the particle’s position change if its coordinates don’t change? Well, in eqn. (16.3) the change in position is represented by the base vectors changing as the object rotates. Thus we could write more explicitly that
| (16.3) |
Here we show more explicitly that the base vectors and depend on .
Just like for polar base vectors (see eqn. (15.1) on page 15.1) we can express the rotating base vectors in terms of the fixed base vectors and .
| (16.4) | |||||
Also we can express the fixed basis vectors in terms of the rotating vectors like this:
| (16.5) | |||||
Please review the section on dot products, 1.2, to see one derivation of these formulae.
We will use the phrase reference frame or just frame to mean “a coordinate system attached to a rigid object”. One can think of the coordinate grid as like an invisible metal framework (hence the word ‘frame’) that rotates with the object. We refer to a calculation based on the rotating coordinates in fig. 16.6 variously as “in the frame ” or “using the frame” or “in the frame
.
In computer calculations we usually manipulate lists and arrays of numbers and not geometric vectors. So on a computer we keep track of vectors by keeping track of their lists of components. Let’s look at a point fixed to the object and whose coordinates we know in the reference configuration:
Assuming the object axes and fixed axes coincide in the reference configuration, the object coordinates of a point are equal to the space fixed coordinates of the point in the reference configuration . We can think of the point as defined either way, so
Here is a question we often need to answer: What are the fixed basis coordinates of a point that has the rotating-frame coordinates ?
Here is one way to find the answer:
| (16.6) | |||||
so we can pull out the and coordinates compactly as,
| (16.7) |
But this can, in turn be written in matrix notation as
| (16.14) | |||||
| (16.15) | |||||
The matrix or is the rotation matrix for counterclockwise rotations by . As shown above, if you know the coordinates of a point fixed on an object before rotation, you can find its coordinates after rotation by multiplying the coordinate column vector by the matrix . You can remember what is by remembering its components or by remembering that
the first and second column of are the components of and , respectively, in the fixed coordinate system.
For example, the first column of consists of the and components of .
A feature of eqn. (16.14) is that the same matrix prescribes the coordinate change for every different point on the object. Thus for points called 1, 2 and 3 we have
A more compact way to write a matrix times a list of column vectors is to arrange the column vectors one next to the other in a matrix. By multiplying this matrix by we get a new matrix whose columns are the new coordinates of various points. For example,
| (16.16) |
Eqn. 16.16 is useful for computer animation of rotating things in video games (and in dynamics simulations too) where points 1,2, and 3 are points on an object.
Example: Rotate a picture
If a simple picture of a house is drawn by connecting the six points (fig. 16.5a) with the first point at , the second at , etc., and the sixth point on top of the first, we have,
After a counter-clockwise rotation about O, the coordinates of the house, in a coordinate system that rotates with the house, are unchanged (fig. 16.5b). But in the fixed (non-rotating, Newtonian) coordinate system the new coordinates of the rotated house points are,
as shown in fig. 16.5c.
SAMPLE 16.1
Computing rotated position of an object: A rigid object AOB is pinned at point O and is free to rotate about this point. Using the rotation matrix, find the coordinates of points A and B when and respectively.
††margin:
Solution Let be a set of axes glued to the object with origin at O which is also the origin of the space-fixed coordinate axes . We first write the coordinates of points A and B in the body-fixed axes.
The rotation matrix for the axes rotating counterclockwise along with the object is
We can find the coordinates of A and B in the rotated position by multiplying their coordinates with the rotation matrix. We first write the coordinates of both points A and B in a single matrix, one column for each, and multiply with to get the new coordinates:
Now, we calculate the new coordinates for the given values of .
For , we have,
Thus,
Answer:
Similarly, for , we get,
Answer:
SAMPLE 16.2
Computer program for object rotation: The object shown in the figure rotates counterclockwise with its angular position given by where . Write a computer program to animate the motion of the object by plotting its position at specified time instants. Use your program to plot the position of the object every two seconds for a total of 18 seconds.
Solution
In order to draw the object on the computer screen, we need to define it by taking a sufficient number of points on its boundary so that plotting all those points with connected line segments represents the object as closely as possible. In this case, we can select all the corner points on its boundary and define the object with the coordinates of these points. With the given geometry, it is fairly easy to find these coordinates. Denoting the coordinates of the th point by , we can define this object by the following set of coordinates:
Note that the last column is a repeat of the first column. This is essential so that when we plot the object using the and coordinates listed here, we get a closed boundary.
Animation of the angular motion of this object basically requires determining the rotated position of the object and plotting it at sufficiently small time intervals. To do this, we need to define the initial orientation (position), define the time increment, find the angle of rotation corresponding to the new time, compute the corresponding rotation matrix, multiply the object coordinates with the rotation matrix, separate out the and coordinates, and plot vs . We then repeat the whole sequence until we reach the final time.
The following pseudocode can be easily adapted to a computer program to do the required animation.
object = [ 0 2 4 4 2 0 0
-1 -1 -2 2 1 1 -1] % coordinates of points
t = 0Ψ % specify initial time
t_final = 36 % specify final time
delta_t = 2 % specify time increment
thetadot0 = 10*pi/180 % specify thetadot0 in rad/s
while t < t_final
t = t + delta_t % increment the time
theta = thetadot0*t % compute current theta
R = [ cos(theta) -sin(theta) % find the rotation matrix
sin(theta) cos(theta)]
new_position = R*object % find new coordinates
x = new_position(1st row) % extract all x-coordinates
y = new_position(2nd row) % extract all y-coordinates
plot y vs x % plot the new position
end
The plot obtained by implementing this code is shown in fig. 16.14.
16.1.1 Give the definition of each term below in words and, if possible, with an equation.
rotation angle
Answer: Rotating frame coordinates of a point.
Answer: Fixed basis coordinates of a point.
rotation matrix
Answer: A matrix that when operated on a column vector of fixed basis coordinates gives the coordinates in a rotated frame.
16.1.2 Assume that in the reference configuration, when , that the axes attached to an object are aligned with the axes. Consider a point P attached to the moving frame (object) that has coordinates . Find each of the quantities below in terms of some or all of and .
Answer:
The rotation matrix . Answer:
16.1.3 Assume that in the reference configuration, when , that the axes are aligned with the axes. Consider two points and attached to the moving frame. and have coordinates . Find each of the quantities below in terms of some or all of and .
Answer:
The rotation matrix . Answer:
16.1.4 The square plate OABC shown in the figure measures on a side. The coordinate axes are fixed to the plate and the axes are fixed in space. The plate is rotated by an angle as shown.
Find the coordinates and of the corner point B in the rotated frame.
Find the coordinates and of the same point in the fixed frame without using the rotation matrix.
Write the rotation matrix for the given rotation.
Find the coordinates of point B in the fixed frame using its coordinates and the rotation matrix .
Find the coordinates of point B in the rotated frame using its coordinates and the rotation matrix in its appropriate form.
Answer: ,
16.1.5 The rod shown in the figure rotates counterclockwise, starting from . Find the following quantities.
Write the position vector of points A and B as a list of numbers in an array, e.g., , at . Answer: ,
Write the rotation matrix for . Answer:
Find the coordinates of points A and B using the rotation matrix R for . Answer:
If the rod rotates at a constant angular speed , find the coordinates of points A and B after 2 seconds, assuming the rod is at when .
16.1.6 Find the angle of rotation corresponding to the following rotation matrices:
. Answer:
. Answer:
. Answer:
16.1.7 The rod shown in the figure rotates with constant angular speed . at .
Find the position vector and the velocity of point B at . Answer: and
Find the position vector of point B at . What is the net displacement of point B during the time interval ? Is this net displacement equal to ? Why not? Answer: Net displacement . The net displacement is not equal to because the motion is not along a straight line.
16.1.8 Write a computer program to animate the rotation of an object. Your input should be a set of and coordinates defining the object (such that plot y vs x draws the object on the screen) and the rotation angle . The output should be the rotated coordinates of the object.
From the geometric information given in the figure, generate coordinates of enough points to define the given object. Answer: (0, 0); (, 15); (, 0); (0, ); (15, ).
Using your program, plot the object at and .
Assume that the object rotates with constant angular speed . Find and plot the position of the object at .
16.1.9 The arrow shaped object shown in the figure is defined by the five points whose coordinates are given in some normalized units. The position shown is at . The time dependent angular position of the object is given by where and . Animate the motion of the object and show its position from to at every second.
Thus far we have talked about rotation, but not how it varies in time. Dynamics is about motion, velocities and accelerations, so we need to think about rotation rates and rotational accelerations.
A 2D rigid object’s net rotation is measured by the rotation angle . Thus, the simplest measure of rotation rate is . Because all marked lines rotate the same amount they all have the same rates of change. So etc. and, as for rotation, the concept of rotation rate of a rigid object transcends the concept of rotation rate of this or that particular line. We give this rotation rate of a rigid object a special name, angular velocity, and symbol, (omega).
Repeating, for all lines marked on a rigid object,
| (16.23) |
Often we think of angular velocity as a vector . Its direction is the axis of the rotation which, for objects in the plane is , pointed in the direction normal to the plane. The scalar part of is . So, the angular velocity vector is
| (16.24) |
with as defined in eqn. (16.23). Note is the angular velocity of the object (and of every line on it)
.
The usual sign convention for is the same as that for , positive is counter-clockwise (CCW) around the origin. That’s the direction your fingers wrap if you point your thumb in the direction and grab the axis.
Our first use of the angular velocity vector is to calculate the rate of change of the rotating unit base vectors and . We can find the rate of change of, say, , by taking the time derivative of the first eqn. (16.4), and using the chain rule while recognizing that . We can also make an analogy with polar coordinates (page 15.1), where we think of as like and as like . We found there that and . Either from differentiating or from the correlation with polar coordinates, we get
| (16.25) |
because and . Depending on the tastes of your lecturer, you may find Equations 16.25 are the most used equations from this point onward
.
Box 16.2 The fixed Newtonian reference frame
Now we can reconsider the concept of a Newtonian frame, a concept which we had to assume to write the equations of dynamics in the first place.
All of mechanics depends, of course, on the laws of mechanics. The laws of mechanics are equations which involve, in part, the positions of things as a function of time. But how position is perceived to change in time depends on your reference frame. And some reference frames are better than others. The best, from our point of view, are reference frames in which Newton’s laws are accurate. Such a reference frame is called a Newtonian frame. In engineering practice the frames we use as approximations of a Newtonian frame often seem, loosely speaking, somehow still. So we sometimes call such a frame the fixed frame and label it with a script capital . When we talk about velocity and acceleration of mass points, for use in the equations of mechanics, we are always talking about the velocity and acceleration relative to a ixed, or equivalently, Newtonian frame.
Assume and are the coordinates of a vector and is a fixed frame with fixed axis (with associated constant base vectors and ). When we write we mean . But we could be more explicit (and notationally ornate) and write the velocity of P in the Newtonian frame as
The in front of the time derivative (or in front of the dot) means that when we calculate a derivative we hold the base vectors of constant. This is no surprise, because for the base vectors are constant. In general, however, when taking a derivative in a given frame you
write vectors in terms of base vectors stuck to the frame, and
only differentiate the components.
But only for the ixed or ewtonian frame will accelerations calculated this way be directly applicable to Newton’s laws.
We will avoid the ornate notation of labeling frames when it is not needed. For example, if you don’t see any script capital letters floating around in front of derivatives, you can assume that we are taking derivatives relative to a fixed Newtonian frame.
Let’s call some rotating object (script capital B) to which is glued a coordinate system with base vectors and .
We now introduce the concept of derivative in a frame which we write, for the frame , as
which means, in words, the rate of change of something as viewed in the rotating frame . Now consider a point P at that is glued to the object. That is, the and coordinates of do not change in time.
Box 16.3 Plato’s discussion of spinning in circles as motion (or not)
Here is how a smart person, Plato, thought about these things two thousand
years ago.
Here he writes about an imagined discussion between Socrates and Glaucon
about how an object can maintain contradictory attributes
simultaneously.
“Socrates: Now let’s have a more precise agreement so that we won’t have any grounds for dispute as we proceed. If someone were to say of a human being standing still, but moving his hands and head, that the same man at the same time stands still and moves, I don’t suppose we’d claim that it should be said like that, but rather that one part of him stands still and another moves. Isn’t that so?
Glaucon: Yes it is.
Socrates: Then if the man who says this should become still more charming and make the subtle point that tops as wholes stand still and move at the same time when the peg is fixed in the same place and they spin, or that anything else going around in a circle on the same spot does this too, we wouldn’t accept it because it’s not with respect to the same part of themselves that such things are at the same time both at rest and in motion. But we’d say that they have in them both a straight and a circumference; and with respect to the straight they stand still since they don’t lean in any direction –while with respect to the circumference they move in a circle; and when the straight inclines to the right the left, forward, or backward at the same time that it’s spinning, then in no way does it stand still.
Glaucon: And we’d be right.”
This chapter is about things that are still with respect to their own parts (they do not distort) but in which the points do move in circles.
That is, relative to a moving frame, the velocity of a point glued to the frame is zero (no surprise).
We would like to know the velocity of such a point in the fixed frame. We just take the derivative, using the product rule and the differentiation rules we have developed for the rotating base vectors:
where is the simple way to write .
Box 16.4 Acceleration of a point on a rotating body, using
Leaving off the ornate pre-super-script for simplicity, we have
| (16.26) | |||||
To continue we need to use the product rule of differentiation for the cross product of two time dependent vectors like this:
| (16.27) | ||||
| (16.28) | ||||
| (16.29) | ||||
| (16.30) |
Substituting back into eqn. (16.26) we get
which derives eqn. (16.32).
Thus,
| (16.31) |
We can rewrite eqn. (16.31) in a minimalist or elaborate notation as
Both are correct. In the first case you have to use common sense to know what point you are talking about and that it is on an object rotating with absolute angular velocity . In the second case everything is laid out perfectly clearly (which is why it looks perfectly confusing). On the left side of the equation it says that we are interested in how point P moves relative to, not just any frame, but the fixed frame . On the right side we make clear that the rotation rate we are looking at is that of object relative to and not some other relative rotation. We further make clear that the formula only makes sense if the position of the point P is measured relative to a point which doesn’t move, namely 0.
What we have just found largely duplicates what we already learned in section 7.1 for points moving in circles. The slight generalization is that the same angular velocity can be used to calculate the velocities of multiple points on one rigid object. But the key idea remains: the velocity of a point going in circles is tangent to the circle it is going around and with magnitude proportional both to distance from the center and to the angular rate of rotation (fig. 16.20a).
Let’s again consider a point stuck on a rotating object and with position
Relative to the frame to which a point is attached, its acceleration is zero (again no surprise). But what is its acceleration in the fixed frame? We find this acceleration by writing the position vector and then differentiating twice, repeatedly using the product rule and eqn. (16.25) we get (see Box 16.2 for the details).
| (16.32) |
which is hardly intuitive at a glance
.
Recalling that in 2D we can use either the right hand rule or manipulation of unit vectors to rewrite eqn. (16.32) as
| (16.33) |
where and
.
Thus, as we found in section 15.1 for a particle going in circles, the acceleration can be written as the sum of two terms, a tangential acceleration due to increasing tangential speed, and a centrally directed (centripetal) acceleration due to the direction of the velocity continuously changing towards the center (see fig. 16.20b). The generalization we have made in this section is that the same can be used to calculate the acceleration for all the different points on one rotating object.
As you well know by now, the position of point B relative to point A is . Similarly the relative velocity and acceleration of two points A and B is defined to be
| (16.34) |
So, the relative velocity (as calculated relative to a fixed frame) of two points glued to one spinning rigid object is given by
where point is the point in the Newtonian frame on the fixed axis of rotation and is the angular velocity of . Repeating,
| (16.35) |
Because points A and B are fixed on their velocities and hence their relative velocity as observed in a reference frame fixed to are all . But, point A has an absolute velocity that is different from that of point B. So they have a relative velocity as seen in the fixed frame. And it is what you would get if B was just going in circles around A.
Similarly, the relative acceleration of two points glued to one rigid object spinning at constant rate is
| (16.36) |
Again, the relative acceleration is due to the difference in the points’ positions relative to the point fixed on the axis. These kinematics results, 16.2 and 16.36, are useful for calculating angular momentum relative to the center-of-mass. They are also sometimes useful for the understanding of the motions of machines with moving connected parts.
Equation 16.35 actually applies to any vector that has the property that it is fixed in the frame which is rotating at (fig. 16.22).
That is, all we needed in the derivation was that ’s components in the - reference frame be constant.
Examples of such would include more than just the relative position of two points fixed on the object, . Both and are also are fixed in the rotating frame. Also imagine a bug moving at a constant speed on a straight line marked on a body. That bug’s velocity relative to the rotating frame would also be constant as seen in the frame. For any whose representation as an arrow moves like a line segment drawn on a rigid object rotating at :
| (16.37) |
This is the most fundamental equation of rigid-object kinematics. Know it well. For future reference, equation (16.37) is also valid in three dimensions. Equation 16.37 is the generalization of the three equations
A coordinate system attached to a rotating rigid object , defines a reference frame (fig. 16.6 on page 16.6).
Recall, the base vectors in this frame change in time just like any other vector fixed in the rotating frame (eqn. (16.37))
Once we know how the base vectors and change with time, as per the above equations, we can find velocity and acceleration by differentiation of position.
For example, consider two points A and B fixed on a rotating object. The position of B relative to A is
where we are using and as a shorthand notation for and (see fig. 16.23).
The coordinates and are constant so
Now we just differentiate the relative position with respect to time,
We could similarly calculate by taking another derivative to get, after a calculation much like that above,
For the above calculation the points A and B were fixed on a rotating part. Most machines have many parts, at least some of which move in more complex ways than just circular motion. We will be able to understand such machines by considering points and parts that are moving relative to a part that is itself rotating. Such will be considered in Chapter 19.
Box 16.5 Angular velocity and the rotation matrix
The rotation matrix and the angular velocity vector are related. Because angular velocity is a rotation rate, in some sense it must be the derivative of the rotation. But the situation is a bit subtle because is a vector and is a matrix.
On the one hand we have the rotation equation:
This gives the coordinates of after rotation in terms of its coordinates before the rotation. Because the coordinates before rotation are fixed we can take the time derivative of both sides to get
On the other hand we have the vector equation
where . We can relate this vector equation to the matrix component equation above it by expressing the cross product with a matrix multiplication as described in box 1.3 on 1.3. Applying that result to this 2D case where
we get
Just keeping the 2D terms we have
But tells us the coordinates in terms of the reference coordinates. So
Comparing the equations above we see that
| (16.40) |
Using the notation from box 1.3 we can rewrite this as
| (16.41) |
Equation 16.41 is the way that angular velocity and the rate of change of the rotation matrix are related.
We can also solve for as
from which we can find in the off diagonal elements if we are given and .
We can understand eqn. (16.41) better, perhaps, by thinking of the rotation matrix as having two columns which are, respectively, the fixed-coordinate-system components of and :
Thus the columns of are the rates of change of these two unit vectors:
With a total abuse of notation (in general one does not allow cross products of a vectors with matrices) we can write the above equation more memorably as
Thus eqn. (16.41) is just a version of the fundamental kinematic equation 16.37 on page 16.37.
SAMPLE 16.3
A uniform bar AB of length rotates counterclockwise about point A with constant angular speed . At the instant shown in fig. 16.24 the linear speed of the center-of-mass C is .
What are the angular speed and angular velocity of the bar?
What is the linear velocity of point B?
By what angles do the angular positions of points C and B change in 2 seconds?
Solution
Let the angular velocity of the bar be where is the angular speed. We first need to find .
The linear speed of point C is given, . Now,
Therefore, the angular velocity of the bar is .
Answer:
Point B is at distance from the pivot point A. Thus it goes around a circle of radius (see fig. 16.26). Therefore,
Answer:
We can also write where
Let be the position of point C at some time and be the position at time . We want to find for .
The change in angular position of point B is the same as that of point C. In fact, all points on AB undergo the same change in angular position because AB is a rigid body.
Answer:
SAMPLE 16.4
A flywheel of diameter is made of cast iron. To avoid extremely high stresses and cracks it is recommended that the peripheral speed not exceed 6000 to 7000 ft/min. What is the corresponding rpm rating for the wheel?
Solution
| Diameter of the wheel | ||||
| radius of wheel |
Now,
Similarly, corresponding to
Thus the rpm rating of the wheel should read 955 – 1114 rpm.
Answer:
SAMPLE 16.5
Two gears A and B have the diameter ratio of 1:2. Gear A drives gear B. If the output at gear B is required to be 150 rpm, what should be the angular speed of the driving gear? Assume no slip at the contact point.
Solution Let C and C′ be the points of contact on gear A and B respectively at some instant . Since there is no relative slip between C and C′, both points must have the same linear velocity at instant . If the velocities are the same, then the linear speeds must also be the same. Thus
Answer:
SAMPLE 16.6
A uniform rigid rod AB of length is connected to two rigid links OA and OB. The assembly rotates at a constant rate about point O in the plane. At the instant shown, when rod AB is vertical, the velocities of points A and B are , and . Find the angular velocity of bar AB. What is the length of the links?
Solution
Let the angular velocity of the rod AB be . ††margin: We know that the rod rotates about the -axis but we do not know the sense of the rotation i.e., or . Here we have assumed that is in the positive direction, although just by sketching we can easily see that must be in the direction. Since we are given the velocities of two points on the rod we can use the relative velocity formula to find :
| (16.42) | |||||
| Thus, | |||||
| (16.43) | |||||
Answer:
Let be the angle between link OA and the horizontal axis. Now,
Dotting both sides of the equation with and we get
| (16.44) | |||||
| (16.45) |
Answer:
SAMPLE 16.7
Verify the relative velocity formula: The motor at O in Fig. 16.32 rotates the ‘L’ shaped bar OAB in counterclockwise direction at an angular speed which increases at . At the instant shown, the angular speed . Each arm of the bar is of length L = 2.
Find the velocity of point A.
Find the relative velocity and use the result to find the absolute velocity of point B .
Find the velocity of point B directly. Check the answer obtained in part (b) against the new answer.
Solution
As the bar rotates, every point on the bar goes in circles centered at point O. Therefore, we can easily find the velocity of any point on the bar using circular motion formula . Thus,
The velocity vector is shown in Fig. 16.33.
Answer:
Point B and A are on the same rigid body. Therefore, with respect to point A, point B goes in circles about A. Hence the relative velocity of B with respect to A is
These velocities are shown in Fig. 16.34.
Answer:
Since point B goes in circles of radius OB about point O, we can find its velocity directly using circular motion formula:
The velocity vector is shown in Fig. 16.35. Of course this velocity is the same velocity as obtained in part (b) above.
Answer:
Note: Nothing in this sample uses !
SAMPLE 16.8
Verify the relative acceleration formula: Consider the ‘L’ shaped bar of Sample 16.31 again. At the instant shown, the bar is rotating at and is slowing down at the rate of .
Find the acceleration of point A.
Find the relative acceleration of point B with respect to point A and use the result to find the absolute acceleration of point B ().
Find the acceleration of point B directly and verify the result obtained in (ii).
Solution
We are given:
Point A is going in circles of radius L. Hence,
Answer:
The relative acceleration of point B with respect to point A is found by considering the motion of B with respect to A. Since both the points are on the same rigid body, point B executes circular motion with respect to point A. Therefore,
and
Answer:
Since point B is going in circles of radius OB about point O, we can find the acceleration of B as follows.
This acceleration is, naturally again, the same acceleration as found in (ii) above.
Answer:
16.2.1 Give the definition of each term below in words and, but for the first, with an equation.
angular velocity Answer: The rate of change of any marked or imaginary line on a rigid object.
angular acceleration Answer: The rate of change of the angular velocity.
16.2.2 Assume that in the reference configuration, when , that the axes attached to an object are aligned with the axes. Consider a point P attached to the moving frame (object) that has coordinates . Find each of the quantities below in terms of some or all of and .
Answer:
Answer:
Answer:
16.2.3 Assume that in the reference configuration, when , the axes are aligned with the axes. Consider two points and attached to the moving frame. and have coordinates . Find each of the quantities below in terms of some or all of and .
Answer:
Answer:
Answer:
16.2.4 Find , if and . Answer:
16.2.5 A rod OB rotates with its end O fixed as shown in the figure with angular velocity and angular acceleration at the instant of interest. Find, draw, and label the tangential and normal acceleration of point B at . Answer: and .
16.2.6 A disc rotates at 15 rpm. How many seconds does it take to rotate by 180 degrees? What is the angular speed of the disc in rad/s? Answer: .
16.2.7 A motor turns a uniform disc of radius counter-clockwise about its mass center at a constant rate . The disc lies in the -plane and its angular displacement is measured (positive counter-clockwise) from the -axis. What is the angular displacement of the disc if it starts at and ? What are the velocity and acceleration of a point at position ? Answer: , , and .
16.2.8 Find the angular velocities of the second, minute, and hour hands of a clock. Answer: , , and .
16.2.9 A disc spins at a constant rate of two revolutions per second counter-clockwise about its geometric center, G, which is fixed. A point P is marked on the disk at a radius of one meter. At the instant of interest, point P is on the -axis of an -coordinate system centered at point G.
Draw a neat diagram showing the disk, the particle, and the coordinate axes.
What is the angular velocity of the disk, ? Answer: .
What is the angular acceleration of the disk, ? Answer: .
What is the velocity of point P? Answer: .
What is the acceleration of point P? Answer: .
16.2.10 A uniform rigid rod rotates at constant speed in the -plane about a peg at point . The center of mass of the rod may not exceed a specified acceleration . Find the maximum angular velocity of the rod. Answer:
16.2.11 A dumbbell AB is welded to a rigid arm OC such that OC is perpendicular to AB. Arm OC rotates about O at a constant angular velocity . At the instant when ,
Find the relative velocity of B with respect to A. Answer:
Find the acceleration of point B relative to point A. Answer:
16.2.12 Two discs and rotate at constant speeds about their centers. Disc rotates at and disc rotates at . Which is rotating faster? Answer: Disc is faster.
16.2.13 A motor turns a uniform disc of radius counter-clockwise about its mass center at a constant rate . The disc lies in the -plane and its angular displacement is measured (positive counter-clockwise) from the -axis. What are the velocity and acceleration of a point at position relative to the velocity and acceleration of a point at position ) on the disk? (.)
Answer: , and
16.2.14 A long rod has many holes along its length such that it can be pegged at any of the various locations. It rotates counter-clockwise at a constant angular speed about a peg whose location is not known. At some instant , the velocity of end is . After , the velocity of end is . If the rod has not completed one revolution during this period,
find the angular velocity of the rod, and Answer: .
find the location of the peg along the length of the rod. Answer: from .
16.2.15 A circular disc of radius rotates in the -plane about a point which is at a distance away from the center of the disk. At the instant of interest, the linear speed of the center is and the magnitude of its centripetal acceleration is .
Find the rotational speed of the disk. Answer: .
Is the given information enough to locate the center of rotation of the disk? Answer: No.
If the acceleration of the center has no component in the direction at the instant of interest, can you locate the center of rotation? If yes, is the point you locate unique? If not, what other information is required to make the point unique? Answer: If the acceleration of the center has no component in the direction the center of rotation can lie anywhere on the x-axis.
16.2.16 A uniform disc of radius is mounted eccentrically on a motor shaft at point . The motor rotates the disc at a constant angular speed. At the instant shown, the velocity of the center of mass is .
Find the angular velocity of the disc. Answer:
Find the point with the highest linear speed on the disc. What is its velocity? Answer: Speed is maximum at ;
16.2.17 The circular disc of radius rotates about its center . At a given instant, point on the disk has a velocity in the direction shown. At the same instant, the tangent of the angle made by the total acceleration vector of any point B with its radial line to is 0.6. Compute the angular acceleration of the disc. Answer:
16.2.18 Show that, for non-constant rate circular motion, the acceleration of all points in a given radial line are parallel.
16.2.19 A motor turns a uniform disc of radius about its mass center at a
variable angular rate with rate of change , counter-clockwise. The disc lies in the -plane
and its angular displacement is measured from the -axis, positive
counter-clockwise. What are the velocity and
acceleration of a point
at position relative to the velocity and acceleration of a
point at position ) on the disk? (.)
Answer: , and
16.2.20 The dumbbell AB shown in the figure rotates counterclockwise about point O with angular acceleration . Bar AB is perpendicular to bar OC. At the instant of interest, and the angular speed is 2 rad/s.
Find the velocity of point B relative to point A. Will this relative velocity be different if the dumbbell were rotating at a constant rate of 2 rad/s? Answer:
Without calculations, draw a vector approximately representing the acceleration of B relative to A.
Find the acceleration of point B relative to A. What can you say about the direction of this vector as the motion progresses in time? Answer:
16.2.21 Bit-stream kinematics of a CD. A Compact Disk (CD) has bits of data etched on concentric circular tracks. The data from a track is read by a beam of light from a head that is positioned under the track. The angular speed of the disk remains constant as long as the head is positioned over a particular track. As the head moves to the next track, the angular speed of the disk changes, so that the linear speed at any track is always the same. The data stream comes out at a constant rate bits/second. When the head is positioned on the outermost track, for which , the disk rotates at .
What is the number of bits of data on the outermost track? Answer: Megabits.
find the angular speed of the disk when the head is on the innermost track (), and Answer: 509.09 rpm.
find the number of bits on the innermost track. Answer: Megabits
16.2.22 2-D constant rate gear train. The angular velocity of the input shaft (driven by a motor not shown) is a constant, . What is the angular velocity of the output shaft and the speed of a point on the outer edge of disc , in terms of , , , and ? Answer: and
Gear B is welded to C and engages with A.
16.2.23 A horizontal disk of diameter is driven at a constant speed of . A small disk can be positioned anywhere between and on disk by sliding it along the overhead shaft and then fixing it at the desired position with a set screw (see the figure). Disk rolls without slip on disk . The overhead shaft rotates with disk and, therefore, its rotational speed can be varied by varying the position of disk . This gear system is called brush gearing. Find the maximum and minimum rotational speeds of the overhead shaft. Answer: and
16.2.24 Two points A and B are on the same machine part that is hinged at an as yet unknown location C. Assume you are given that points at positions and are supposed to move in given directions, indicated by unit vectors and . For each of the problem parts below, illustrate your results with two numerical examples (in consistent units): i) , , , and (thus ), and ii) a more complex example of your choosing.
Describe in detail what equations must be satisfied by the point .
Write a computer program that takes as input the 4 pairs of numbers , , and and gives as output the pair of numbers .
Find a formula of the form that explicitly gives the position vector for point C in terms of the 4 given vectors. Answer:
All of the basic equations of Dynamics concern the motion of mass. In particular, linear momentum, angular momentum, and energy are all expressed as sums over all the bits of mass in a system, with each bit of mass multiplied by some terms concerning position, velocity and acceleration.
From the earlier sections in this chapter we know how to find the velocity and acceleration of every bit of mass on a 2-D rigid object as it spins about a fixed axis. So it is just a matter of doing integrals or sums to calculate the various momentum and energy quantities of interest. As an object moves and rotates the region of integration and the values of the integrands change. So, in principle, in order to analyze a rigid object one has to evaluate a different integral or sum at every different configuration. But there is a shortcut: for a rotating rigid object a sum (over all atoms, say), or a difficult integral (for example, over the complex region representing a machine part) is reduced to simple multiplication.
The moment of inertia
††margin: The COM is the most common reference point. The moment of inertia of a given object depends on what reference point is used. Most often when people say ‘the’ moment of inertia they mean the moment of inertia with respect to the center-of-mass. For clarity, this moment of inertia is often notated as in this book. If a different reference point, say point is used, the moment of inertia is notated as .
simplifies the expressions for the angular momentum, the rate of change of angular momentum, and the energy of a rigid object. For more general motions the shortcuts need a matrix But for 2D mechanics only one component of the matrix is relevant, it is , called just or for short.
Here are the main results. A flat object spinning with in the plane has a mass distribution which gives, by means of a calculation which we will discuss shortly, a moment of inertia so that:
| (16.46) | |||||
| (16.47) | |||||
| (16.48) |
There are two main skills you need to develop associated with .
Using the formulas above properly in angular momentum and energy balance. This is covered in sec. 16.4.
Finding I of an object, as is covered in most freshman calculus texts and is reviewed in this section.
Some facts about moment of inertia are summarized on the inside back cover.
The definition of moment of inertia
is
where , and are the distances in the and directions of the bits of mass () from the center of mass and is the total distance. Similarly if all distances are measured relative to a point C or O (instead of relative to the center of mass) then similar formulas as above calculate or . See box 16.3 on page 16.3 for the calculation of the moments of inertia of some simple objects.
The term is sometimes called the polar moment of inertia, or polar mass moment of inertia to distinguish it from the and , terms which have little utility in planar dynamics ( and , as area inertias rather than mass inertias, are all-important when calculating the bending stiffness or the stress in elastic beams, however).
What, physically, is the polar moment of inertia? It is a measure of the extent to which mass is far from the given reference point. Every bit of mass contributes to in proportion to the square of its distance from the reference point. As we will show in sec. 16.4 on page 16.4 (see eqn. (16.51) on page 16.51) is just the quantity we need to do mechanics problems.
Another measure of the extent to which mass is spread from the reference point, besides the moment of inertia , is the radius of gyration,
.
The radius of gyration is defined as:
That is, the radius of gyration of an object is the radius of an equivalent ring of mass that has the same and the same mass as the given object. It is an equivalent single radius for the whole mass distribution
.
Example: Radius of gyration of a hoop
The radius of gyration about its center of a circular hoop of mass and radius about its center is .
For the most part it is which is of primary interest. Other reference points are useful
If the rigid object is hinged at a fixed point then can be used in a slight short cut in calculation of angular momentum and energy; and
If one wants to calculate the moment of inertia of a composite object about its center-of-mass then it is useful to first find the moment of inertia of each of its parts about that system’s center of mass (which is generally not the center-of-mass of any of the separate parts).
The planar parallel axis theorem is the equation
| (16.49) |
In this equation is the distance from the center-of-mass to a line parallel to the -axis which passes through point . See box 16.3 on page 16.3 for a derivation of the parallel axis theorem for planar objects.
Note that , always.
Moment of inertia of complex objects. One can calculate the moment of inertia of a composite object about its center of mass, in terms of the masses and moments of inertia of the separate parts. Say the position of the center of mass of is relative to a fixed origin, and the moment of inertia of that part about its center of mass is . We can then find the moment of inertia of the composite about its center-of-mass by the following sequence of calculations, in order:
You can reduce this recipe to one grand formula with lots of summation signs. But you would end up doing the calculations in about the order prescribed above in any case. This sequence of steps lends itself naturally to a computer spread sheet or to any program that deals easily with arrays of numbers. This is laid out for the similar center-of-mass calculation on page 2.2
.
Example: Inertia of a semicircle about its COM.
The moment of inertia of a semicircle with mass and radius about O (see fig. 16.52) is
From the example on page 2.2 the center of mass is at
By the parallel axis theorem
which we can solve for to get
The radius of gyration is that radius so that and is
The perpendicular axis theorem for planar objects is the equation
which is derived in box 16.3 on page 16.3. It gives the ‘polar’ inertia in terms of the inertias and . Unlike the parallel axis theorem, the perpendicular axis theorem does not have a three-dimensional counterpart. The theorem is of greatest utility when one wants to study the three-dimensional mechanics of a flat object and thus are in need of its full moment of inertia matrix.
Box 16.6 The 2-D parallel axis theorem and the perpendicular axis theorem
Sometimes one wants to know the moment of inertia relative to the center of mass. And, sometimes, if the object is held at a hinge joint at , relative to that hinge point There is a simple relation between these two moments of inertia known as the parallel axis theorem.
For the two-dimensional mechanics of two-dimensional objects, our only concern is and and not the full moment of inertia matrix. In this case, and . Now, let’s prove the theorem in two dimensions referring to the figure.
The cancellation comes from the definition of center of mass.
Sometimes, people write the parallel axis theorem more simply as
using the symbol to mean . One thing to note about the parallel axis theorem is that the moment of inertia about any point is always greater than the moment of inertia about the center of mass. For a given object, the minimum moment of inertia is about the center-of-mass.
Why the name parallel axis theorem? We use the name because the two ’s calculated are the moments of inertia about two parallel axes (both in the direction) through the two points and .
One way to think about the theorem is the following. The moment of inertia of an object about a point not at the center-of-mass is the same as that of the object about the plus that of a point mass located at the center-of-mass. If the distance from to the is larger than the outer radius of the object, then the term is larger than . The distance of equality of the two terms is the radius of gyration, .
For planar objects,
Similarly,
That is, the moment of inertia about the -axis is the sum of the inertias about the two perpendicular axes and . Note that the objects must be planar ( everywhere) or the theorem would not be true. For example, for a three-dimensional object.
Box 16.7 Some examples of 2-D Moment of Inertia
Here, we illustrate some simple moment of inertia calculations for two-dimensional objects. The needed formulas are summarized, in part, by the lower right corner components (that is, the elements in the third column and third row (3,3)) of the matrices in the table on the inside back cover.
If we assume that all mass is concentrated at one or more points, then the integral
reduces to the sum
which reduces to one term if there is only one mass,
So, if , , and , then . Note that, in this case, since the radius from the center-of-mass to the center-of-mass is zero.
In this case, the sum that defines reduces to two terms, so
Note that, if , then .
Consider a thin rod with uniform mass density, , per unit length, and length . We calculate as
If either or , then this expression reduces to . If , then is at the center-of-mass and
We can illustrate one last point. With a little bit of algebraic histrionics of the type that only hindsight can inspire, you can verify that the expression for can be arranged as follows:
That is, the moment of inertia about point is greater than that about the center of mass by an amount equal to the mass times the distance from the center-of-mass to point squared. This derivation of the parallel axis theorem is for one special case, that of a uniform thin rod.
For a hoop of uniform mass density, , per unit length, we might consider all of the points to have the same radius . So,
Or, a little more tediously,
This is the same as for a single point mass at a distance from the origin . It is also the same as for two point masses if they both are a distance from the origin. For the hoop, however, is at the center-of-mass so which is not the case for a single point mass.
Assume the disk has uniform mass density, , per unit area. For a uniform disk centered at the origin, the center-of-mass is at the origin so
For example, a plate of radius has the same moment of inertia as a hoop with a radius.
For the special case that the center of the plate is at point , the center-of-mass is also at and .
Note that for all planar objects (the perpendicular axis theorem). For a uniform rectangle, . But the integral is just the term often used for , the area moment of inertia, in strength of materials calculations for the stresses and stiffnesses of beams in bending. You may recall that for a rectangle. Similarly, . So, the polar moment of inertia can be recalled by remembering the area moment of inertia of a rectangle combined with the perpendicular axis theorem.
SAMPLE 16.9
Moment of inertia of point masses: A pendulum is made up of two unequal point masses and connected by a massless rigid rod of length . The pendulum is pivoted at distance along the rod from the small mass.
Find the moment of inertia of the pendulum.
Find the moment of inertia of the pendulum.
Find the radius of gyration of the pendulum.
Solution
First we need to find the center of mass of the system. Let the center of mass C be located at distance from the origin O. An easy way to find the location of C will be to consider the pendulum to lie along the x-axis with the origin at O (see fig. 16.54). Then
So, the moment of inertia about the center of mass is
Answer:
We can calculate in two different ways, one directly by summing the contributions of each mass about point O, and the other by using and the parallel axis theorem.
Answer:
The radius of gyration, by definition, is the distance that gives the desired moment of inertia if the entire mass of the system is concentrated there. Thus, if is the radius of gyration for the moment of inertia about point O, then
Thus the radius of gyration of the given pendulum is .
Answer:
SAMPLE 16.10
Moment of inertia of a rod: A uniform rigid rod AB of mass and length swings about the -axis passing through the pivot point O.
Find the moment of inertia of the bar using the fundamental definition .
Find using the parallel axis theorem given that where total mass, and total length of the rod. (You can find for many commonly encountered objects in the table on the inside backcover of the text).
Solution
Since we need to carry out the integral, , to find , let us consider an infinitesimal
length segment of the bar at distance from the pivot point
O. (see Figure 16.57). Let the mass of the infinitesimal
segment be .
Now the mass of the segment may be written as
| (mass per unit length of the bar) (length of the segment) | ||||
We also note that the distance of the segment from point O, . Substituting the values found above for and in the formula we get
Answer:
The parallel axis theorem states that
Since the rod is uniform, its center-of-mass is at its geometric center, i.e., at distance from either end. From the Fig 16.58 we can see that
Answer:
SAMPLE 16.11
Moment of inertia of a wheel with a cut-out: A uniform rigid wheel of radius is made eccentric by cutting out a portion of the wheel. The center-of-mass of the eccentric wheel is at C, a distance from the geometric center O. The mass of the wheel (after deducting the cut-out) is 3.2 lbm. The moment of inertia of the wheel about point O, , is 1.8 . We are interested in the moment of inertia of the wheel about points A and B on the perimeter.
Without any calculations, guess which point, A or B, gives a higher moment of inertia. Why?
Calculate , and and compare with the guess in (a).
Solution
The moment of inertia should be higher. Moment of inertia measures the geometric distribution of mass about the -axis. But the distance of the mass from the axis counts more than the mass itself (). The distance of the mass appears as a quadratic term in . The total mass is the same whether we take the moment of inertia about point A or about point B. However, the distribution of mass is not the same about the two points. Due to the cut-out being closer to point B there are more “’s” at greater distances from point B than from point A. So, we guess that
Answer:
If we know the moment of inertia (about the center-of-mass) of the wheel, we can use the parallel axis theorem to find and . In the problem, we are given . But,
Answer:
Clearly, , as guessed in (a).
SAMPLE 16.12
Moment of inertia: modeling a sphere as a point mass: A uniform solid sphere of mass and radius is attached to a massless rigid rod of length . The sphere swings in the plane. Find the error in calculating as a function of if the sphere is treated as a point mass concentrated at the center-of-mass of the sphere.
Solution
The exact moment of inertia of the sphere about point O can be calculated using parallel axis theorem:
If we treat the sphere as a point mass, the moment of inertia is
Therefore, the relative error in is
| error | ||||
From the above expression we see that for the error is very small. From the graph of error in Fig. 16.62 we see that even for , the error in due to approximating the sphere as a point mass is less than 2%.
16.3.1 Find an expression for the polar moment of inertia for each of the following systems about the specified points. Answer in terms of some or all of the variables given. If a position vector is given then you can assume the coordinates are given so that .
A particle with mass at about
Its center of mass
The origin
Answer: , .
Two particles with masses and a distance apart located at and , respectively, about
Mass 1
Mass 2
The center of mass of the system
The origin
Answer: , , , and
A collection of particles at locations about
The origin
The center of mass G of the system
Answer: and
A uniform line segment with length and mass with center of mass at location about
Its center of mass G
One end
The origin
Answer: , , and
.
A uniform hoop with mass and radius and center at about
Its center of mass G
A point on the hoop
The origin
Answer: , , and
.
A uniform disk with mass and radius and center at about
Its center of mass G
A point on the perimeter
The origin
Answer: , , and
.
A uniform rectangle with mass and sides and and with center at about
Its center of mass G
A corner
The origin
Answer: , , and
.
16.3.2 Two objects have masses and and polar moments of inertia and about their respective centers of mass G1 and G2 which are a distance apart. What is the moment of inertia of the system about its center of mass? Answer: .
16.3.3 A point mass is located at and in the -plane. Find the moment of inertia of the mass about the -axis. Answer: .
16.3.4 A small ball of mass is attached to a massless rod.
What is the value of of the ball about the other end of the rod? Answer: .
How much must you shorten the rod to reduce the moment of inertia of the ball by half? Answer: .
16.3.5 Two identical point masses are attached to the two ends of a rigid massless bar of length (one mass at each end). Locate a point along the length of the bar about which the polar moment of inertia of the system is more than that calculated about the mid point of the bar. Answer: At from either end
16.3.6 A dumbbell consists of a rigid massless bar of length and two identical point masses and , one at each end of the bar.
About which point on the dumbbell is its polar moment of inertia a minimum and what is this minimum value? Answer: , about the midpoint.
About which point on the dumbbell is its polar moment of inertia a maximum and what is this maximum value? Answer: , about either end
16.3.7 A light rigid rod of length has a point mass at end and a point mass at end . Point is the center of mass of the system. First, answer the following questions without any calculations and then do calculations to verify your guesses.
About which point , , or , is the polar moment of inertia of the system a minimum? Answer:
About which point is a maximum? Answer:
What is the ratio of and ? Answer:
Is the radius of gyration of the system greater, smaller, or equal to the length of the rod? Answer: smaller,
16.3.8 Three identical particles of mass are connected to three identical massless rods of length and welded together at point as shown in the figure.
Guess (no calculations) which of the three moment of inertia terms , , is the smallest and which is the biggest. Answer: Biggest: ; smallest: .
Calculate the three moments of inertia to check your guess. Answer: . .
If the orientation of the system is changed, so that one mass is along the -axis, will your answer to part (a) change? Answer: Answer does not change.
Find the radius of gyration of the system for the polar moment of inertia. Answer: .
16.3.9 Show that the polar moment of inertia of the uniform bar of length and mass , shown in the figure, is , in two different ways:
by using the basic definition of polar moment of inertia , and
by computing first and then using the parallel axis theorem.
16.3.10 Approximately locate the center of mass of the tapered rod shown in the figure and compute the polar moment of inertia . [Hint: use the variable thickness of the rod to define an approximate equivalent variable mass density per unit length.] Answer: and .
16.3.11 A short rod of mass and length hangs from an inextensible and massless rod of length .
Find the moment of inertia of the rod. Answer: .
Find the moment of inertia of the rod by considering it as a point mass located at its center of mass. Answer: .
Find the percent error in in treating the bar as a point mass by comparing the expressions in parts (a) and (b). Plot the percent error versus . For what values of is the percentage error less than ? Answer: Percent error in is . The percentage error less than if .
16.3.12 Parallel axis theorem? A massless square plate has four identical point masses located at its corners.
Find the polar moment of inertia . Answer:
Find a point on the plate about which the system’s moment of inertia is maximum. Answer: or , often, but not always, diagonally opposite the corner with the largest mass
Find the radius of gyration of the system about the rectangle center. Answer:
16.3.13 For the massless square plate with four point masses on the corners, the polar moment of inertia . Find of the system. Answer: .
16.3.14 A uniform square plate ( on edge) has a corner cut out. The total mass of the remaining plate is 3.
Where is the center of mass of the plate at the instant shown? Answer: .
What is the moment of inertia of the plate about point O? Answer: .
What is the moment of inertia of the plate about its center of mass? Answer: .
16.3.15 A uniform thin triangular plate of mass , height , and base lies in the -plane.
Set up the integral to find the polar moment of inertia of the plate. Answer: .
Show that by evaluating the integral in part (a).
Locate the center of mass of the plate and calculate . Answer: .
16.3.16 A uniform thin plate of mass is cast in the shape of a semi-circular disk of radius as shown in the figure.
Find the location of the center of mass of the plate Answer: .
Find the polar moment of inertia of the plate, . [Hint: It may be easier to set up and evaluate the integral for and then use the parallel axis theorem to calculate .] Answer: .
16.3.17 A uniform square plate of side has a circular cut-out of radius . The mass of the plate is .
Find the polar moment of inertia of the plate about its center. Answer: .
Plot versus .
Find the limiting values of for and .
Answer: and
.
16.3.18 A uniform thin circular disk of radius and mass has a rectangular slot of width cut into it as shown in the figure. You can approximate the right edge of the cutout as a straight line (thus your calculation would be exact for a system that has a very thin sliver of material closing the gap at the right).
Find the polar moment of inertia of the disk. Answer: .
Locate the center of mass of the disk and calculate . Answer: , , and .
estimate the error in your calculation of the mass, center of mass location and moment of inertia due to not subtracting out the sliver of material at the right edge of the slot. Answer: Not subtracting out the sliver of material at the right edge of the slot leads to an error of .
16.3.19 Calculate the polar moment of inertia of a uniform square plate with mass and side various ways. Choose numerical values for and if you want to address this as a numerical calculation rather than a theoretical one.
As a single uniform plate
As a composite of 4 plates, each with sides
As a composite of 100 plates, each with sides (These first three methods should give exactly the same answer).
For the composite of 100 plates now neglect the terms which concern the moments of inertia of each small square about its center of mass. What is the error in making this neglect? Answer: Exactly 1% .
We now want to find what we can about the forces on and motions of a single planar object that is hinged at one point. The most famous examples are pendula and gears. But many things move like pendula, like people or towers falling over or boats rocking in the water. And many machine parts besides gears rotate about a bearing including propellers, lawn mower blades, and printing-press rollers. None of these things are strictly planar objects, but, for many analyses, the planar approximation is more or less accurate
.
The basic method, as always, is to use a free-body diagram and kinematics results to evaluate terms in the linear momentum, angular momentum and energy balance equations. Then these equations are used to calculate unknown forces or to find differential equations of motion. The terms involving force and moment are found from a free-body diagram exactly as in statics. The terms involving motion, namely momentum, angular momentum and energy, are evaluated using the various tools developed earlier in this chapter.
We can evaluate all the momentum and energy terms in the equations of motion (inside cover), namely: and (for any reference point C of our choosing) if we can calculate the velocity and acceleration of every point in the system. For circular motion of a rigid object about point C, we know from sec. 16.2 that the velocities and accelerations of a point at are
where is the angular velocity of the object relative to a fixed frame. This situation is a little more complex than the special case of straight-line motion in chapter 14, where all points in a system had the same acceleration as each other, but is still quite manageable.
The most general case of 2D rigid-object circular motion is an arbitrarily shaped 2D rigid object with arbitrary and . The situation shown in fig. 16.75 shows the general case (but in that example exactly two forces are applied, as opposed to an arbitrary number).
For any system in any motion we have the general result that
For a rigid object, the center-of-mass is a particular point G that is fixed relative to the object. So the velocity and acceleration of that point can be expressed the same way as for any other point. So, for an object in planar rotational motion about 0
and
If the center-of-mass is at 0 the momentum and its rate of change are both zero. But if the center-of-mass is off the axis of rotation, there must be a net force on the object with a component parallel to (if ) and also a component orthogonal to (if ). This net force might be applied at 0 or G or any other place(s) on the object.
The angular momentum itself is easy enough to calculate, using the short hand notation that is the position vector of a point relative to hinge point O.
| (16.50) |
Here eqn. (16.50)c is a vector equation. But since both sides are in the direction we can dot both sides with to get the scalar moment equation eqn. (16.50)d, taking both and as positive when counterclockwise.
The moment of inertia shortcut. Note that an integral over all the mass of has been reduced to a simple multiplication by using that
This shortcut is especially useful because does not change as the object rotates about 0.
However, the moment of inertia is never needed for problem solving. One can just evaluate, for example, by evaluating the integral for the specific problem at hand.
To get the all important term for angular momentum balance, namely , for this system we can differentiate eqn. (16.50). We can also use the general expression for to write the angular momentum balance equation as follows.
| (16.51) |
We get from eqn. (16.51)e to eqn. (16.51)f by noting that is perpendicular to . Thus, using the right hand rule twice we get .
Eqn. 16.51f and eqn. (16.51)g are the vector and scalar versions for the rate of change of angular momentum with respect to point 0 for rotation of a planar object about 0. Repeating,
| (16.52) |
Again, can always be evaluated either of two ways: 1) as a sum or integral of over all the mass or 2) using the shortcut with .
Let’s assume that there are a set of point forces applied to an object
.
And, to be contrary, let’s assume the mass is continuously distributed (the derivation for rigidly connected point masses would be similar). The power balance equation for one rotating rigid object is (discussed below):
| (16.53) |
When not directly labeled, positions and moments are assumed to be relative to the hinge at 0. Derivation 16.53 is two derivations in one. The left side about power and the right side about kinetic energy. Let’s discuss one at a time.
Power. On the left side of eqn. (16.53) we note in (c) that the power of each force is the dot product of the force with the velocity of the point it touches. In (d) we use what we know about the velocities of points on rotating rigid bodies. In (e) we use the vector identity (see sec. 1.3, 1.3). In (f) we note that is common to all points so factors out of the sum. In (g) we note that is the moment of the force about pt O. And in (g) we sum the moments of the forces. So the power of a set of forces acting on a rigid object is the product of their net moment (about 0) and the object angular velocity,
| (16.54) |
Kinetic energy. On the right side of eqn. (16.53) we note in (c) that the kinetic energy is the sum of the kinetic energy of the mass increments. In (d) we use what we know about the velocities of these bits of mass, given that they are on a common rotating object. In (e) we use that the magnitude of the cross product of orthogonal vectors is the product of the magnitudes () and that the dot product of a vector with itself is its magnitude squared (). In (f) we factor out because it is common to all the mass increments and note that the remaining integral is constant in time for a rigid object. In (g) we carry out the derivative. In (h) we de-simplify the result from (g) in order to show a more general form that we will find later in 3D mechanics. Eqn. (h) follows from (g) because is parallel to for 2D rotations.
Note that we started here with the basic power balance equation from the front inside cover. Instead, we could have derived power balance from our angular momentum balance expression (see box 16.4 on 16.4).
Box 16.8 The relation between angular momentum balance and power balance
For this system, angular momentum balance can be derived from power balance and vice versa. Thus neither is essentially more fundamental than the other and both are reliable. First we can derive power balance from angular momentum balance as follows:
| (16.55) |
That is, when we dot both sides of the angular momentum equation with we get on the left side a term which we recognize as the power of the forces and on the right side a term which is the rate of change of kinetic energy.
The opposite derivation starts with the power balance fig. 16.53(g)
| (16.56) |
and, assuming , divide by to get the angular momentum equation for planar rotational motion.
Example: Spinning disk
The round flat uniform disk in fig. 16.76 is in the plane spinning at the constant rate about its center. It has mass and radius . What force is required to cause this motion? What torque? What power?
From linear momentum balance we have:
Which we could also have calculated by evaluating the integral instead of using the general result that . From angular momentum balance we have:
So the net force and moment needed are and . Like a particle that moves at constant velocity with no force, a uniform disk rotates at constant rate with no torque (at least in 2D).
Once one knows the velocity and acceleration of all points in a system one can find all of the motion quantities in the equations of motion by adding or integrating using the defining sums from earlier chapters. This addition or integration is an impractical task for many motions of many objects where the required sums may involve billions and billions of atoms or a difficult integral. As you recall from earlier chapters, the linear momentum and the rate of change of linear momentum can be calculated by just keeping track of the center-of-mass of the system of interest. One wishes for something so simple for the calculation of angular momentum.
It turns out that we are in luck if we are only interested in the two-dimensional motion of two-dimensional rigid bodies. The luck is not so great for 3-D rigid bodies but still there is some simplification. For general motion of non-rigid bodies there is no simplification to be had.
The simplification is to use the moment of inertia for the bodies rather than evaluating the momenta and energy quantities as integrals and sums. Of course one may have to do a sum or integral to evaluate or but once this calculation is done, one need not work with the integrals while worrying about the dynamics. At this point we will assume that you are comfortable calculating and looking-up moments of inertia. We proceed to use it for the purposes of studying mechanics.
For constant rate rotation, we can calculate the velocity and acceleration of various points on a rigid body using and . So we can calculate the various motion quantities of interest: linear momentum , rate of change of linear momentum , angular momentum , rate of change of angular momentum , and kinetic energy .
Consider a two-dimensional rigid body like that shown in fig. 16.77.
Now let us consider the various motion quantities in turn. First the linear momentum . The linear momentum of any system in any motion is . So, for a rigid body spinning at constant rate about point O (using ):
Similarly, for any system, we can calculate the rate of change of linear momentum as . So, for a rigid body spinning at constant rate,
That is, the linear momentum is correctly calculated for this special motion, as it is for all motions, by thinking of the body as a point mass at the center-of-mass.
Unlike the calculation of linear momentum, the angular momentum turns out to be something different than would be calculated by using a point mass at the center of mass. You can remember this important fact by looking at the case when the rotation is about the center-of-mass (point O coincides with the center-of-mass). In this case one can intuitively see that the angular momentum of a rigid body is not zero even though the center-of-mass is not moving. Here’s the calculation just to be sure:
We have defined the ‘polar’ moment of inertia as . In order to calculate for a specific body, assuming uniform mass distribution for example, one must convert the differential quantity of mass into a differential of geometric quantities. For a line or curve, ; for a plate or surface, , and for a 3-D region, . , , and are differential line, area, and volume elements, respectively. In each case, is the mass density per unit length, per unit area, or per unit volume, respectively. To avoid clutter, we do not define a different symbol for the density in each geometric case. The differential elements must be further defined depending on the coordinate systems chosen for the calculation; e.g., for rectangular coordinates, or, for polar coordinates, .
Since and always point in the direction for two dimensional problems people often just think of angular momentum as a scalar and write the equation above simply as ‘,’ the form usually seen in elementary physics courses.
The derivation above has a feature that one might not notice at first sight. The quantity called does not depend on the rotation of the body. That is, the value of the integral does not change with time, so is a constant. So, perhaps unsurprisingly, a two-dimensional body spinning about the -axis through has constant angular momentum about if it spins at a constant rate. ††margin: Note that the angular momentum about some other point than O will not be constant unless the center-of-mass does not accelerate (i.e., is at point O).
Now, of course we could find this result about constant rate motion of 2-D bodies somewhat more cumbersomely by plugging in the general formula for rate of change of angular momentum as follows:
| (16.57) |
Using moment of inertia about the center of mass. Often it is easier to think of the motion as composed of two parts, motion of the center of mass and motion relative to the center of mass, as explained in box 16.4 on page 16.4. Thus we have two terms for angular momentum and its rate of change :
| (16.58) | ||||
| (16.59) |
Kinetic energy. Finally, we can calculate the kinetic energy by adding up for all the bits of mass on a 2-D body spinning about the -axis:
| (16.60) |
The kinetic energy can also be written as a sum of contributions of motion of the center of mass and motion relative to the center of mass,
| (16.61) |
Example: Pendulum disk
For the disk shown in fig. 16.78, we can calculate the rate of change of angular momentum about point as
Alternatively, we could calculate directly
Note that we are using the planarity of the objects and of their motion for our calculations
.
The equation for linear momentum balance is the same as always, we just need to calculate the acceleration of the center-of-mass of the spinning body.
| (16.62) |
Finally, the kinetic energy for a planar rigid body rotating in the plane is:
Box 16.9 Simplifying using the center of mass
The definition of angular momentum relative to a point C is
If we rewrite as
and
then
So,
The reason is somewhat intuitive. It is what you would calculate if you were looking for the center-of-mass relative to the center of mass. More formally,
Similarly, because it is what you would calculate if you were looking for the velocity of the center-of-mass relative to the center of mass.
The central result of this box is that
angular momentum of any system is that due to motion of the center-of-mass plus motion relative to the center-of-mass.
SAMPLE 16.13
A rod in a constant rate circular motion: A uniform rod of mass and length is connected to a motor at end O. A ball of mass is attached to the rod at end B. The motor turns the rod in the counterclockwise direction at a constant angular speed . There is gravity pointing in the direction. Find the torque applied by the motor (i) at the instant shown and (ii) when . How does the torque change if the angular speed is doubled?
Solution
The FBD of the rod and ball system is shown in Fig. 16.80(a). Since the system is undergoing circular motion at a constant speed, the acceleration of the ball as well as every point on the rod is just radial (pointing towards the center of rotation O) and is given by where is the radial distance from the center O to the point of interest and is a unit vector along OB pointing away from O (Fig. 16.80(b)).
Angular Momentum Balance about point O gives
| (16.63) | |||||
| (16.64) | |||||
Answer:
(ii) Substituting the given values of in the above expression we get
Answer:
It is clear from the expression of the torque that it does not depend on the value of the angular speed ! Therefore, the torque will not change if the speed is doubled. In fact, as long as the speed remains constant at any value, the only torque required to maintain the motion is the torque to counteract the moments due to gravity at O.
SAMPLE 16.14
At the onset of circular motion: A rectangular plate of mass is pivoted at one of its corners as shown in the figure. The plate is released from rest in the position shown. Find the force on the support immediately after release.
Solution The free-body diagram of the plate is shown in Fig. 16.82. The force applied on the plate by the support is unknown.
The linear momentum balance for the plate gives
| (16.65) | |||||
Thus to find we need to find .
The angular momentum balance for the plate about the fixed support point O gives
| where | ||||
| and | ||||
| Thus, | ||||
From eqn. (16.65), the support force is now readily calculated:
Using the given numerical values of , and , , and , we get
Answer:
SAMPLE 16.15
A compound gear train. When the gear of an input shaft, often called the driver or the pinion, is directly meshed in with the gear of an output shaft, the motion of the output shaft is opposite to that of the input shaft. To get the output motion in the same direction as that of the input motion, an idler gear is used. If the idler shaft has more than one gear in mesh, then the gear train is called a compound gear train.
In the gear train shown in Fig. 16.83, the input shaft is rotating at 2000 rpm and the input torque is 200 N-m. The efficiency (defined as the ratio of output power to input power) of the train is 0.96 and the various radii of the gears are: and . Find
the input power and the output power ,
the output speed , and
the output torque.
Solution
The power:
Answer:
The angular speed of meshing gears can be easily calculated by realizing that the linear speed of the point of contact has to be the same irrespective of which gear’s speed and geometry is used to calculate it. Thus,
Answer:
The output torque,
Answer:
SAMPLE 16.16
An accelerating gear train. In the gear train shown in Fig. 16.85, the torque at the input shaft is and the angular acceleration is . The radii of the various gears are: and and the moments of inertia about the shaft axis passing through their respective centers are: . Find the output torque of the gear train.
Solution
Since the difference between the input power and the output power is used in accelerating the gears, we may write
Let be the output torque of the gear train. Then,
| (16.66) |
Now,
| (16.68) |
The velocity or acceleration of the point of contact between two meshing gears has to be the same irrespective of which meshing gear’s geometry and motion is used to compute them.
The different ’s and the ’s can be related by realizing that the linear speed or the tangential acceleration of the point of contact between any two meshing gears has to be the same irrespective of which gear’s speed and geometry is used to calculate it. Thus, using the linear speed and tangential acceleration calculations for points P and R in Fig. 16.86, we can find , and , . Considering the linear speed of point P, we find
and considering the tangential acceleration of point P, , we find
Similarly,
and
But
and
Now equating the two quantities, and , and canceling from both sides, we obtain
Answer:
SAMPLE 16.17
Drums used as pulleys. Two drums, A and B of radii and are welded together. The combined mass of the drums is and the combined moment of inertia about the -axis passing through their common center O is . A string attached to and wrapped around drum B supports a mass . The string wrapped around drum A is pulled with a force as shown in Fig. 16.87. Assume there is no slip between the strings and the drums. Find
the angular acceleration of the drums,
the tension in the string supporting mass , and
the acceleration of mass .
††margin:
Solution
The free-body diagram of the drums and the mass are shown in Fig. 16.88 separately where is the tension in the string supporting mass and and are the support reactions at O. Since the drums can only rotate about the -axis, let
Now, let us do angular momentum balance about the center of rotation O:
Since the motion is restricted to the -plane (i.e., 2-D motion), the rate of change of angular momentum may be computed as
Setting we get
| (16.69) |
Now, let us write linear momentum balance, , for mass :
Do we know anything about acceleration of the mass? Yes, we know its direction () and we also know that it has to be the same as the tangential acceleration of point D on drum B (why?). Thus,
| (16.70) | |||||
Therefore,
| (16.71) |
Calculation of : We now have two equations, (16.69) and (16.71), and two unknowns, and . Subtracting times Eqn.(16.71) from Eqn. (16.69) we get
Answer:
Calculation of tension : From equation (16.71):
Answer:
Calculation of acceleration of the mass: Since the acceleration of the mass is the same as the tangential acceleration of point D on the drum, we get (from eqn. (16.70))
Answer:
Comments: It is important to understand why the acceleration of the mass is the same as the tangential acceleration of point D on the drum. We have assumed (as is common practice) that the string is massless and inextensible. Therefore each point of the string supporting the mass must have the same linear displacement, velocity, and acceleration as the mass. Now think about the point on the string which is momentarily in contact with point D of the drum. Since there is no relative slip between the drum and the string, the two points must have the same vertical acceleration. This vertical acceleration for point D on the drum is the tangential acceleration .
SAMPLE 16.18 Energy method for pulley dynamics: Consider the pulley problem of Sample 16.86 again. Use energy method to
find the angular acceleration of the pulley, and
the acceleration of the mass.
Solution
In energy method we use speeds, not velocities. Therefore, we have to be careful in our thinking about the direction of motion. In the present problem, let us assume that the pulley rotates and accelerates clockwise. Consequently, the mass moves up against gravity.
The energy equation we want to use is
The power is given by where the sum is carried out over all external forces. For the mass and pulley system the external forces that do work are and .
The other external forces on the system—the reaction force of the support point O and the weight of the pulley—are acting at point O (see fig. 16.89). But, since point O is stationary, these forces do no work.
Therefore,
| (16.72) | |||||
Now we need to calculate the rate of change of kinetic energy . There are two objects here that have kinetic energy—the hanging mass and the pulley. Hence,
The hanging mass has pure translational motion and hence its kinetic energy is
where is the linear speed of the mass. If we assume the string to be inextensible, then the linear speed of the mass has to be the same as the tangential speed of point D of the pulley. Thus and
The pulley, on the other hand, has pure rotational motion about point O, and hence its kinetic energy is given by
Summing the two kinetic energies and differentiating with respect to time , we get
| (16.73) | |||||
| (16.74) |
Now equating the power and the rate of change of kinetic energy from eqns. eqn. (16.72) and eqn. (16.74), we get
From kinematics of circular motion,
Substituting these values in the power balance equation above, we get
Since the sign of is positive, our initial assumption of clockwise acceleration of the pulley is correct.
Answer:
From kinematics of circular motion,
Answer:
SAMPLE 16.19 Energy Accounting: Consider the pulley problem of Sample 16.86 again.
What percentage of the input energy (work done by the applied force ) is used in raising the mass by 1 ?
Where does the rest of the energy go? Provide an energy-balance sheet.
Solution
Let and be the input energy and the energy used in raising the mass by 1, respectively. Then the percentage of energy used in raising the mass is
Thus we need to calculate and to find the answer. is the work done by the force on the system during the interval in which the mass moves up by 1. Let be the displacement of the force during this interval. Since the displacement is in the same direction as the force (we know it is from Sample 16.86), the input-energy is
So to find we need to find .
For the mass to move up by 1 the inner drum B must rotate by an angle where
Since the two drums, A and B, are welded together, drum A must rotate by as well. Therefore the displacement of force is
and the energy input is
Now, the work done in raising the mass by is
Therefore, the percentage of input-energy used in raising the mass
The rest of the energy () goes in accelerating the mass and the pulley. Let us find out how much energy goes into each of these activities. Since the initial state of the system from which we begin energy accounting is not prescribed (that is, we are not given the height of the mass from which it is to be raised 1, nor do we know the velocities of the mass or the pulley at that initial height), let us assume that at the initial state, the angular speed of the pulley is and the linear speed of the mass is . At the end of raising the mass by 1 from this state, let the angular speed of the pulley be and the linear speed of the mass be . Then, the energy used in accelerating the pulley is
Similarly, the energy used in accelerating the mass is
We can calculate the percentage of input energy used in these activities to get a better idea of energy allocation. Here is the summary table:
| Activities | Energy Spent | |
|---|---|---|
| in Joule | as % of input energy | |
| In raising the mass by 1 | 19.62 | 49.05% |
| In accelerating the mass | 0.25 | 0.62 % |
| In accelerating the pulley | 20.13 | 50.33 % |
| Total | 40.00 | 100 % |
So, what would you change in the set-up so that more of the input energy is used in raising the mass? Think about what aspects of the motion would change due to your proposed design.
SAMPLE 16.20
Equation of motion of a swinging stick: A uniform bar of mass and length is pinned at one of its ends O. The bar is displaced from its vertical position by an angle and released (Fig. 16.90).
Find the equation of motion using momentum balance.
Find the reaction at O as a function of .
Solution
First we draw a simple sketch of the given problem showing relevant
geometry (Fig. 16.90(a)), and then a free-body diagram of the bar
(Fig. 16.90(b)).
(a) A line sketch of the swinging rod and (b) free-body diagram of the rod.
We should note for future reference that
Equation of motion using momentum balance: We can write angular momentum balance about point O as
Let us now calculate both sides of this equation:
| (16.75) | |||||
| (16.76) | |||||
| (16.77) |
Answer:
Reaction at O: Using linear momentum balance
Dotting both sides of with and and rearranging, we get
Now substituting the expression for from (16.77) in and , we get
| (16.78) | |||||
| (16.79) |
Answer:
Check: We can check the reaction force in the special case when the rod does not swing but just hangs from point O. The forces on the bar in this case have to satisfy static equilibrium. Therefore, the reaction at O must be equal to and directed vertically upwards. Plugging and (no motion) in Eqn. (16.78) and (16.79) we get and , the values we expect.
SAMPLE 16.21
Swinging stick dynamics using moment of inertia: A uniform bar of mass and length is pinned at one of its ends O. The bar is displaced from its vertical position by an angle and released (Fig. 16.93). Find the equation of motion of the stick.
††margin:
Solution
We repeat the problem solved in Sample 16.88 here with just one different step of finding the rate of change of angular momentum with the help of moment of inertia formula. As usual, we first draw a free-body diagram of the bar (Fig. 16.94). We assume, , and
We can write angular momentum balance about point O as
Let us now calculate both sides of this equation:
| (16.80) | |||||
| (16.81) | |||||
| (16.82) |
| (16.83) |
Answer:
SAMPLE 16.22 Numerical solution of the swinging stick motion: For the swinging stick considered in Samples 16.88 or 16.97, find the time that the rod takes to fall from to if it is released from rest at ?
Solution
The given initial angle is a big value of – big in that we cannot assume (obviously ). Therefore we may not use the linearized equation (16.86) to solve for explicitly. We have to solve the full nonlinear equation (16.83) to find the required time. Unfortunately, we cannot get a closed form solution of this equation using mathematical skills you have at this level. Therefore, we resort to numerical integration of this equation.
For numerical integration, we need to first write the given differential equation as a set of first order ordinary differential equations. To do so, we introduce as a new variable and rewrite eqn. (16.83) as
| (16.84) | |||||
| (16.85) |
Now we need to specify the initial conditions and the time duration for integration, and solve the equations using some ODE solver program. Here is a pseudo-code that lists the steps:
g = 9.81, L = 1 % define constants
ODES = { thetadot = omega
omegadot = -3*g/(2*L) * sin(theta) }
ICs = { theta_0 = pi/2
omega_0 = 0 }
solve ODES with ICs for t = 0 to 4 s
plot theta vs t and plot omega vs t
The results obtained from the numerical solution are shown in Fig. 16.96.
The problem of finding the time taken by the bar to fall from to numerically is nontrivial. It is called a boundary value problem. We have only illustrated how to solve initial value problems. However, we can get a fairly good estimate of the time from the solution obtained.
We first plot against time as shown in fig. 16.97. We find the values of and the corresponding values of that bracket . Now, we can use linear interpolation to find the value of at . Proceeding this way, we get (seconds), a little more than we get from the linear ODE in sample 16.97 ().
Comments: Additionally, we can also get the value of when using similar interpolation. In fact, from the vs plot, we find that at , . How does this result compare with the analytical value of from sample 16.97 (which did not depend on the small angle approximation)? Well, we found that
Thus, we get a fairly accurate value from numerical integration.
SAMPLE 16.23
The swinging stick dynamics with energy balance: Consider the same swinging stick as in Sample 16.88. The stick is, again, displaced from its vertical position by an angle and released (See Fig. 16.90).
Find the equation of motion using energy balance.
What is at if ?
Find the period of small oscillations about .
Solution
Equation of motion using energy balance: We use the power equation, , to derive the equation of motion of the bar. Now, the kinetic energy is given by
where is the speed of the infinitesimal mass element . Referring to fig. 16.99, we can write, , and . Thus,
and, therefore,
Calculation of power (): There are only two forces acting on the bar, the reaction force, and the force due to gravity, . Since the support point O does not move, no work is done by . Therefore,
Note that the negative sign stands for the work done against gravity.
Now,
Therefore,
Equating and we get
Answer:
This equation is, of course, the same as we obtained using balance of angular momentum in Sample 16.88.
Find at : We are given that at and (released from rest). This position is (1) shown in Fig. 16.100. In position (2) , i.e., the rod is vertical. Since there are no dissipative forces, the total energy of the system remains constant. Therefore, taking datum for potential energy as shown in Fig. 16.100, we may write
Answer:
Period of small oscillations: The equation of motion is
For small ,
| (16.86) | |||||
| or | |||||
| where |
Therefore,
| the circular frequency | ||||
| and the time period T |
Answer:
[Say for ]
SAMPLE 16.24 The swinging stick with a destabilizing torque. Consider the swinging stick of Sample 16.88 once again.
Find the equation of motion of the stick, if a torque is applied at end O and a force is applied at the other end A.
What happens if is big?
Solution
A free-body diagram of the bar is shown in Fig. 16.101. Once again, we can use to derive the equation of motion as in Sample 16.88. We calculated and in Sample 16.88. Calculation of remains the same in the present problem. We only need to recalculate .
| and | ||||
Therefore, from
Answer:
Now, setting and we get
| (16.87) |
Numerical Solution: We can numerically integrate (16.87) just as in the previous Sample to find . Here is the pseudo-code that can be used for this purpose.
g = 9.81, L = 1 % specify parameters
m = 1, C = 4
ODES = { thetadot = omega
omegadot = -(3*g/(2*L)) * sin(theta)
+ 3*C/(m*L^2) * theta }
ICs = { thetazero = pi/20
omegazero = 0 }
solve ODES with ICs until t = 10
Using this pseudo-code, we find the response of the pendulum. Figure 16.102 shows different responses for various values of . Note that for , it is the same case as unforced bar pendulum considered above.
From Fig. 16.102 it is clear that the bar has periodic motion for small , with the period of motion increasing with increasing values of . It makes sense if you look at Eqn. (16.87) carefully. Gravity acts as a restoring force while the applied torque acts as a destabilizing force. Thus, with the resistance of the applied torque, the stick swings more sluggishly making its period of oscillation bigger.
From Fig. 16.102, we see that at about the stability of the system changes completely. is not periodic anymore. It keeps on increasing at a faster and faster rate, that is, the bar makes complete loops about point O with ever increasing speed. Does it make physical sense? Yes, it does. As the value of is increased beyond a certain value (can you guess the value?), the applied torque overcomes any restoring torque due to gravity. Consequently, the bar is forced to rotate continuously in the direction of the applied force.
SAMPLE 16.25
A torsional pendulum with linear springs: A uniform rigid bar of mass and length is pinned at one end and connected to two springs, each with spring constant , at the other end. The bar is tweaked slightly from its vertical position. It then oscillates about its original position. The bar is timed for 20 full oscillations which take 12.5 seconds. Ignore gravity.
Find the equation of motion of the rod.
Find the spring constant .
What should be the spring constant of a torsional spring if the bar is attached to one at the bottom and has the same oscillating motion characteristics?
Solution
Refer to the free-body diagram in figure 16.104. Angular momentum balance for the rod about point O gives
Thus
However, for small and ,
| (16.88) |
Answer:
Comparing Eqn. (16.88) with the standard harmonic oscillator equation , we get
From the measured time for 20 oscillations, the time period (time for one oscillation) is
Now equating the measured with the derived expression for we get
Answer:
If the two linear springs are to be replaced by a torsional spring at the bottom, we can find the spring constant of the torsional spring by comparison. Let be the spring constant of the torsional spring. Then, as shown in the free-body diagram (see figure 16.105), the restoring torque applied by the spring at an angular displacement is . Now, writing the angular momentum balance about point O, we get
Comparing with the standard harmonic equation, we find the angular frequency
If this system has to have the same period of oscillation as the first system, the two angular frequencies must be equal, i.e.,
Answer:
SAMPLE 16.26
A spring loaded seesaw: A kid, modelled as a point mass with , is sitting at end B of a rigid rod AB of negligible mass. The rod is supported by a spring at end A and a pin at point O. The system is in static equilibrium when the rod is horizontal. Someone pushes the kid vertically downwards by a small distance and lets go. Given that ; find
the unstretched (relaxed) length of the spring,
the equation of motion (a differential equation relating the position of the mass to its acceleration) of the system, and
the natural frequency of the system.
If the rod is pinned at the midpoint instead of at O, what is the natural frequency of the system? How does the new natural frequency compare with that of a mass simply suspended by a spring with the same spring constant?
Solution
Static Equilibrium: The FBD of the (rod + mass) system is shown in Fig. 16.107. Let the stretch in the spring in this position be and the relaxed length of the spring be . The balance of angular momentum about point O gives:
Answer:
Equation of motion: As point B gets displaced downwards by a distance , point A moves up by a proportionate distance . From geometry, ††margin: Here, we are considering a very small so that we can ignore the arc the point mass B moves on and take its motion to be just vertical (i.e., for small ).
Therefore, the total stretch in the spring, in this position,
Now, Angular Momentum Balance about point O gives:
| (16.89) | |||||
| (16.90) | |||||
| (16.91) |
Answer:
The natural frequency of the system: We may also write the previous equation as
| (16.92) |
Substituting and in the expression for we get the natural frequency of the system
Answer:
Comparison with a simple spring mass system:
When , the equation of motion (16.92) becomes
and the natural frequency of the system is simply
which corresponds to the natural frequency of a simple spring mass system shown in Fig. 16.108.
In our system (with ) any vertical displacement of the mass at B induces an equal amount of stretch or compression in the spring which is exactly the case in the simple spring-mass system. Therefore, the two systems are mechanically equivalent. Such equivalences are widely used in modeling complex physical systems with simpler mechanical models.
16.4.1 An object consists of a massless bar with two attached masses and . The object is hinged at .
What is the moment of inertia of the object about point ()? Answer: .
Given , , and , what is , the angular momentum about point ? Answer: .
Given , , and , what is , the rate of change of angular momentum about point ? Answer: .
Given , , and , what is , the total kinetic energy? Answer: .
Assume that you don’t know , or but you do know that is applied to the rod, perpendicular to the rod at . What is ? (Neglect gravity.) Answer:
If were applied to instead of , would be bigger or smaller? Answer: will be bigger if is applied to instead of .
16.4.2 A uniform circular disc rotates at constant angular speed about the origin, which is also the center of the disc. It’s radius is . It’s total mass is .
What is the total force and moment required to hold it in place (use the origin as the reference point of angular momentum and torque). Answer: .
What is the total kinetic energy of the disk? Answer: .
16.4.3 The hinged disk of mass (uniformly distributed) is acted upon by a force shown in the figure. Determine the initial angular acceleration and the reaction forces at the pin . Answer: and .
16.4.4 A motor turns a bar. A uniform bar of length and mass is turned by a motor whose shaft is attached to the end of the bar at O. The angle that the bar makes (measured counter-clockwise) from the positive axis is . Neglect gravity.
Draw a free-body diagram of the bar.
Find the force acting on the bar from the motor and hinge at s. Answer: .
Find the torque applied to the bar from the motor at s. Answer: .
What is the power produced by the motor at s? Answer: .
16.4.5 A physical pendulum. A swinging stick is sometimes called a ‘physical’ pendulum. Take the ‘body’, the system of interest, to be the whole stick.
Draw a free-body diagram of the system.
Write the equation of angular momentum balance for this system about point .
Evaluate the left-hand-side as explicitly as possible in terms of the forces showing on your Free-Body Diagram.
Evaluate the right hand side as completely as possible. You may use the following facts:
where is the distance along the pendulum from the top, is the angle by which the pendulum is displaced counter-clockwise from the vertically down position, is vertically down, and is to the right. You will have to set up and evaluate an integral. Answer: .
16.4.6 A uniform one meter bar is hung from a hinge that is at the end. It is allowed to swing freely. .
What is the period of small oscillations for this pendulum? Answer: .
Suppose the rod is hung from one end. What is the period of small oscillations for this pendulum? Can you explain why it is longer or shorter than when it is hung by its end? Answer: The period of oscillations is longer if the rod is hung at because the torque due to weight of the portion above the hinge negates the torque due to the portion below the hinge. One way to see it is that if the bar is hinged at mid, net torque about the hinge becomes zero leading to infinitely large period of oscillation.
16.4.7 Motor turns a dumbbell. Two uniform bars of length and mass are welded at right angles. At the ends of the horizontal bar are two more masses . The bottom end of the vertical rod is attached to a hinge at O where a motor keeps the structure rotating at constant rate (counter-clockwise). What is the net force and moment that the motor and hinge cause on the structure at the instant shown? Answer: , .
16.4.8 The structure shown in the figure consists of two point masses connected by three rigid, massless rods such that the whole structure behaves like a rigid body. The structure rotates counterclockwise at a constant rate of . At the instant shown, find the force in each rod. Answer: , , and .
16.4.9 Balancing a system of rotating particles. A wire frame structure is made of four concentric loops of massless and rigid wires, connected to each other by four rigid wires presently coincident with the and axes. Three masses, , and , are glued to the structure as shown in the figure. The structure rotates counter-clockwise at a constant rate . There is no gravity.
Find the net force exerted by the structure on the support at the instant shown. Answer: (a) .
You are to put a mass at an appropriate location on the third loop so that the net force on the support is zero. Find the appropriate mass and the location on the loop. Answer: The net force on the support becomes zero if we put on radius loop at .
16.4.10 Motor turns a bent bar. Two uniform bars of length and uniform mass are welded at right angles. One end is attached to a hinge at O where a motor keeps the structure rotating at a constant rate (counterclockwise). What is the net force and moment that the motor and hinge cause on the structure at the instant shown.
neglecting gravity Answer: Net force: , Net moment: .
including gravity. Answer: Net force: , Net moment: .
A bent bar is rotated by a motor.
16.4.11 A uniform disk of mass and radius rotates about a hinge in the -plane. A point mass is fixed to the disk at a distance from the hinge. A motor at the hinge drives the disk/point mass assembly with constant angular acceleration . What torque at the hinge does the motor supply to the system? Answer: .
16.4.12 A rigid rod of length and total mass is held fixed at one end and whirled around in circular motion at a constant rate in the horizontal plane. Ignore gravity.
Find the tension in the rod as a function of , the radial distance from the center of rotation to any desired location on the rod. Answer:
Where does the maximum tension occur in the rod? Answer: at ; i.e., at the center of rotation.
At what distance from the center of rotation does the tension drop to half its maximum value? Answer: .
16.4.13 A thin uniform circular disc of mass and radius rotates in the plane about its center of mass point . Driven by a motor, it has rate of change of angular speed proportional to angular position, . The disc starts from rest at .
What is the rate of change of angular momentum about the origin at ? Answer: .
What is the torque of the motor at ? Answer: .
What is the total kinetic energy of the disk at ? Answer: .
16.4.14 Neglecting gravity, calculate at the instant shown for the system in the figure. Answer: .
16.4.15 The uniform square shown is released from rest at . What is immediately after release? Answer:
16.4.16 Acceleration of a trap door. A uniform bar AB of mass and a ball of the same mass are released from rest from the same horizontal position. The bar is hinged at end A. There is gravity.
Which point on the rod has the same acceleration as the ball, immediately after release. Answer: Point at from A
What is the reaction force on the bar at end A just after release? Answer: directed upwards.
16.4.17 A disk with radius has a string wrapped around it which is pulled with a force . The disk is free to rotate about the axis through normal to the page. The moment of inertia of the disk about is . A point is marked on the string. Given that (0) = 0 and that (0) = 0, what is ? Answer:
16.4.18 A uniform stick with length and mass is welded to a pulley hinged at the center . The pulley has negligible mass and radius . A string is wrapped many times around the pulley. At time , the pulley, stick, and string are at rest and a force is suddenly applied to the string. How long does it take for the pulley to make one full revolution? Ignore gravity. Answer: .
String wraps around a pulley with a stick glued to it.
16.4.19 Constant speed gear train. Gear A is connected to a motor (not shown) and gear B, which is welded to gear C, is connected to a taffy-pulling mechanism. Assume you know the torque and angular velocity of the input shaft. Assume the bearings and contacts are frictionless.
What is the input power? Answer: .
What is the output power? Answer: .
What is the output torque , the torque that gear C applies to its surroundings in the clockwise direction? Answer:
16.4.20 At the input to a gear box a force is applied to gear A. At the output, the machinery (not shown) applies a force of to the output gear. Gear rotates at constant angular rate , clockwise.
What is the angular speed of the right gear? Answer:
What is the velocity of point ? Answer:
What is ? Answer:
If the gear bearings had friction, would have to be larger or smaller in order to achieve the same constant velocity? Answer: will have to be smaller.
If instead of applying a to the left gear it is driven by a motor (not shown) at constant angular speed , what is the angular speed of the right gear? Answer: The left gear is constrained to rotate at .
Two gears with end loads.
16.4.21 A bevel gear system. A bevel-type gear system, shown in the figure, is used to transmit power between two shafts that are perpendicular to each other. The driving gear has a mean radius of and rotates at a constant speed . The mean radius of the driven gear is and the driven shaft is expected to deliver a torque of . Assuming no power loss, find the input torque supplied by the driving shaft. Answer: .
A bevel gear
16.4.22 Belt drives are used to transmit power between parallel shafts. Two parallel shafts, apart, are connected by a belt passing over the pulleys and fixed to the two shafts. The driver pulley rotates at a constant . The speed ratio between the pulleys and is 1:2.5. The input torque is . Assume no loss of power between the two shafts.
Find the input power. Answer:
Find the rotational speed of the driven pulley Answer: .
Find the output torque at . Answer:
16.4.23 In the belt drive system shown, assume that the driver pulley rotates at a constant angular speed . If the motor applies a constant torque on the driver pulley, show that the tensions in the two parts, and , of the belt must be different. Which part has a greater tension? Does your conclusion about unequal tension depend on whether the pulley is massless or not? Assume any dimensions you need. Answer:
16.4.24 Two racks connected by three constant rate gears. A force is applied to one rack. At the output, the machinery (not shown) applies a force of to the other rack.
Assume the gear-train is spinning at constant rate and is frictionless. What is ? Answer: .
If the gear bearings had friction would that increase or decrease to achieve the same constant rate?
If instead of applying a to the left rack it is driven by a motor (not shown) at constant speed , what is the speed of the right rack? Answer: .
If the angular velocity of the gear is increasing at rate does this increase or decrease at the given .
Two racks connected by three gears.
16.4.25 -3-D accelerating gear train. This is really a 2-D problem; each gear turns in a different parallel plane. Shaft B is rigidly connected to gears and . meshes with gear . Gears and are both rigidly attached to shaft AD. Gear meshes with which is welded to shaft A. Shaft A and shaft B spin independently. Assume you know the torque , angular velocity and the angular acceleration of the input shaft. Assume the bearings and contacts are frictionless.
What is the input power? Answer: .
What is the output power? Answer: .
What is the angular velocity of the output shaft? Answer: .
What is the output torque ? Answer: .
A 3-D gear train.
16.4.26 Gear A with radius is rigidly connected to a drum B with radius . The combined moment of inertia of the gear and the drum about the axis of rotation is . Gear A is driven by gear C which has radius . As the drum rotates, a mass is pulled up by a string wrapped around the drum. At the instant of interest, The angular speed and angular acceleration of the driving gear are and , respectively. Find the acceleration of the mass . Answer:
16.4.27 A spindle and pulley arrangement is used to hoist a mass as shown in the figure. Assume that the pulley is to be of negligible mass. When the motor is running at a constant ,
Find the velocity of the mass at . Answer: .
Find the tension in strings and . Answer: , .
16.4.28 Accelerating rack and pinion. The two gears shown are welded together and spin on a frictionless bearing. The inner gear has radius and negligible mass. The outer disk has radius and a uniformly distributed mass of . They are loaded as shown with the force on the massless rack which is held in place by massless frictionless rollers. At the time of interest the angular velocity is (though is not constant). The point P is on the disk a distance from the center. At the time of interest, point P is on the positive axis.
What is the speed of point P? Answer: .
What is the velocity of point P? Answer: .
What is the angular acceleration of the gear? Answer: .
What is the acceleration of point P? Answer: .
What is the magnitude of the acceleration of point P? Answer: .
What is the rate of increase of the speed of point P? Answer: .
Accelerating rack and pinion
16.4.29 Two racks connected by a gear. A force is applied to one rack. At the output the machinery (not shown) applies a constant force to the other rack.
Assume the gear is spinning at constant rate and is frictionless. What is ? Answer: .
If the gear bearing had friction, would that increase or decrease to achieve the same constant rate? Answer: If gear bearing has friction, has to increase.
If the angular velocity of the gear is increasing at rate , does this increase or decrease at the given . Answer: will decrease if the angular velocity of the gear is increasing.
If the output load is given then the motion of the machine can be found, assuming friction is negligible. Assume that the machine starts from rest with a given output load. So long as rack B moves down and the output force is positive, the output power is positive. The machine starts from rest and runs for a fixed amount of time . For what value of is the output work maximum? Once you find an answer, can you find an intuitive explanation for the answer? Answer: For work to be maximum, . Initially work increases with but too large leads to small thus giving smaller net work.
Two racks connected by a gear.
16.4.30 A rack and pinion constrained by a linear spring. Neglect gravity. The spring is relaxed when the angle . Assume the system is released from rest at . What is the acceleration of the point at the end of the stick when ? Answer in terms of any or all of , , , , , , and . [Hint: There are several steps of reasoning required. You might want to draw FBD(s), use angular momentum balance, set up a differential equation, solve it, plug values into this solution, and use the result to find the quantities of interest. Answer:
16.4.31 A small particle of mass is attached to the end of a thin rod of mass M (uniformly distributed), which is pinned at hinge O, as depicted in the figure.
Obtain the equation of motion governing the rotation of the rod. Answer: .
What is the natural frequency of the system for small oscillations ? Answer:
16.4.32 A rigid massless rod has two equal masses and () attached to it at distances and , respectively, measured along the rod from a frictionless hinge located at a point . The rod swings freely from the hinge. There is gravity. Let denote the angle of the rod measured from the vertical. Assume that and are known at the instant of interest.
What is ? Find in terms of , , , and . Answer: .
What is the force of the hinge on the rod? Solve in terms of , , , , , and any unit vectors you may need to define. Answer: .
Would you get the same answers if you put a mass at ? Why or why not? Answer: Answer will change because although total moment of force about the hinge stays same but moment of inertia reduces. In this case we get .
16.4.33 For the pendula in the figure :
Without doing any calculations, try to figure out the relative durations of the periods of oscillation for the five pendula (i.e. the order, slowest to fastest) Assume small angles of oscillation.
Calculate the period of small oscillations. [Hint: use balance of angular momentum about the point 0]. Answer:
(b) has a longer period than (e) does since in (b) the moment of inertia about the center of mass (located at the same position as the mass in (e)) is non-zero.
Rank the relative duration of oscillations and compare to your intuitive solution in part (a), and explain in words why things work the way they do.
16.4.34 A pegged compound pendulum. A uniform bar of mass and length hangs from a peg at point C and swings in the vertical plane about an axis passing through the peg. The distance from the center of mass of the rod to the peg can be changed by putting the peg at some other point along the length of the rod.
Find the angular momentum of the rod about point C. Answer: .
Find the rate of change of angular momentum of the rod about C. Answer: .
How does the period of the pendulum vary with ? Show the variation by plotting the period against . [Hint, you must first find the equations of motion, linearize for small , and then solve.] Answer:
Find the total energy of the rod (using point C as a datum for potential energy). Answer: .
Find when . Answer: .
Find the reaction force on the rod at C, as a function of , , , , and . Answer: .
For the given rod, what should be the value of (in terms of ) in order to have the fastest pendulum? Answer:
Test of Schuler’s pendulum. The pendulum with the value of obtained in (g) is called the Schuler’s pendulum. It is not only the fastest pendulum but also the “most accurate pendulum”. The claim is that even if changes slightly over time due to wear at the support point, the period of the pendulum does not change much. Verify this claim by calculating the percent error in the time period of a pendulum of length under the following three conditions: (i) initial and after some wear , (ii) initial and after some wear , and (iii) initial and after some wear . Which pendulum shows the least error in its time period? What is the connection between this result and the plot obtained in (c)? Answer: The error in the time period of the pendulum for the three cases is (i) , (ii) , and (iii) respectively.
16.4.35 A uniform stick of length and mass is a hair away from vertically up position when it is released with no angular velocity (a ‘hair’ is a technical word that means ‘very small amount, zero for some purposes’). It falls to the right. What is the force on the stick at point when the stick is horizontal. Solve in terms of , , , , and . Carefully define any coordinates, base vectors, or angles that you use. Answer: .
16.4.36 A massless 10 meter long bar is supported by a frictionless hinge at one end and has a kg point mass at the other end. It is released at from a tip angle of radians measured from vertically upright position (hinge at the bottom). Use .
Using a small angle approximation and the solution to the resulting linear differential equation, find the angle of tip at s and s. Use a calculator, not a numerical integrator. Answer: Angle of tip is at and at .
Using numerical integration of the non-linear differential equation for an inverted pendulum find at s and s.
Make a plot of the angle versus time for your numerical solution. Include on the same plot the angle versus time from the approximate linear solution from part (a).
Comment on the similarities and differences in your plots.
16.4.37 A zero length spring (relaxed length ) with stiffness supports the pendulum shown.
Find assuming , . Answer: .
Find as a function of and (and , , , and .) Answer: .
[Hint: use vectors (otherwise it’s hard)]
[Hint: For the special case, , the solution simplifies greatly.]
16.4.38 A spring-mass-damper system is depicted in the figure. The horizontal damping
force applied at B is given by
The dimensions and parameters are as follows:
For small , assume that and .
Determine the natural circular frequency of small oscillations about equilibrium for the pendulum shown. The static equilibrium position is (pendulum hanging vertically), so the spring is at its rest point in this position.
Idealize the pendulum as a point mass attached to a rigid massless rod of length , so . Also use the “small angle approximation” where appropriate. Answer: .
Sketch a graph of as a function of if the pendulum is released from rest at position when . Your graph should show the correct qualitative behavior, but calculations are not necessary.
16.4.39 The asymmetric dumbbell shown in the figure is pivoted in the center and also attached to a spring at one quarter of its length from the bigger mass. When the bar is horizontal, the compression in the spring is . At the instant of interest, the bar is at an angle from the horizontal; is small enough so that . If, at this position, the velocity of mass ‘’ is and that of mass is , evaluate the power term () in the energy balance equation. Answer: .
16.4.40 The dumbbell shown in the figure has a torsional spring with spring constant (torsional stiffness units are ). The dumbbell oscillates about the horizontal position with small amplitude . At an instant when the angular velocity of the bar is , the velocity of the left mass is and that of the right mass is . Find the expression for the power of the spring on the dumbbell at the instant of interest. Answer: .
16.4.41 A square plate with side and mass is hinged at one corner in a gravitational field . Find the period of small oscillation. Answer: .
16.4.42 A thin hoop of radius and mass is hung from a point on its edge and swings in its plane. Assuming it swings near to the position where its center of mass is below the hinge:
What is the period of its swinging oscillations? Answer: .
If, instead, the hoop was set to swinging in and out of the plane would the period of oscillations be greater or less? Answer: The period will reduce to .
What length of simple pendulum (a massless rod with a point mass at the end) has the same period of oscillation of this hoop, swinging in its plane? Answer:
16.4.43 Oscillating disk. A uniform disk with mass and radius pivots around a frictionless hinge at its center. It is attached to a massless spring which is horizontal and relaxed when the attachment point is directly above the center of the disk. Assume small rotations and the consequent geometrical simplifications. Assume the spring can carry compression. What is the period of oscillation of the disk if it is disturbed from its equilibrium configuration? [You may use the fact that, for the disk shown, , where is the angle of rotation of the disk.]. Answer: period = .
16.4.44
Robotics problem: Simplest balancing of an inverted pendulum. You are holding a stick upside down, one end is in your hand, the other end sticking up. To simplify things, think of the stick as massless but with a point mass at the upper end. Also, imagine that it is only a two-dimensional problem (either you can ignore one direction of falling for simplicity or imagine wire guides that keep the stick from moving in and out of the plane of the paper on which you draw the problem).
You note that if you model your holding the stick as just having a stationary hinge then you get . Assuming small angles, this hinge leads to exponentially growing solutions. Upside-down sticks fall over. How can you prevent this falling?
One way to do keep the stick from falling over is to firmly grab it with your hand, and if the stick tips, apply a torque in order to right it. This corrective torque is (roughly) how your ankles keep you balanced when you stand upright. Your task in this assignment is to design a robot that keeps an inverted pendulum balanced by applying appropriate torque.
Your model is: Inverted pendulum, length , point mass , and a hinge at the bottom with a motor that can apply a torque . The stick might be tipped an angle from the vertical. A horizontal disturbing force is applied to the mass (representing wind, annoying friends, etc).
Draw a picture and a FBD
Write the equation for angular momentum balance about the hinge point. Answer: .
Imagine that your robot can sense the angle of tip and its rate of change and can apply a torque in response to that sensing. That is you can make any function of and that you want. Can you find a function that will make the pendulum stay upright? Make a guess (you will test it below).
Test your guess the following way: plug it into the equation of motion from part (b), linearize the equation, assume the disturbing force is zero, and see if the solution of the differential equation has exponentially growing (i.e. unstable) solutions. Go back to (c) if it does and find a control strategy that works.
Pick numbers and model your system on a computer using the full non-linear equations. Use initial conditions both close to and far from the upright position and plot versus time.
If you are ambitious, pick a non-zero forcing function (say a sine wave of some frequency and amplitude) and see how that affects the stability of the solution in your simulations.
16.4.45 A thin rod of mass and length is hinged with a torsional spring of stiffness at A, and is connected to a thin disk of mass and radius at B. The spring is uncoiled when . Determine the natural frequency of the system for small oscillations , assuming that the disk is:
welded to the rod, and Answer: .
pinned frictionlessly to the rod. Answer: . Frequency higher than in (a)
16.4.46 Assume that the pulley shown in figure(a) rotates at a constant speed . Let the angle of contact between the belt and pulley surface be . Assume that the belt is massless and that the condition of impending slip exists between the pulley and the belt. The free-body diagram of an infinitesimal section of the belt is shown in figure(b).
Write the equations of linear momentum balance for section of the belt in the and directions.
Eliminate the normal force from the two equations in part (a) and get a differential equation for the tension in terms of the coefficient of friction and The contact angle . Answer: , where .
Show that the solution to the equation in part (b) satisfies , where and are the tensions in the lower and the upper segments of the belt, respectively.
16.4.47 A belt drive is required to transmit power from a diameter pulley rotating at a constant to a diameter pulley. The centers of the pulleys are located apart. The coefficient of friction between the belt and pulleys is .
(See problem 16.4.) Draw a neat diagram of the pulleys and the belt-drive system and find the angle of lap, the contact angle , of the belt on the driver pulley.
Find the rotational speed of the driven pulley. Answer: .
The belt in use has a rectangular cross-section. Find the maximum tensile stress in the belt. Answer: MPa.
16.4.48 Slippery money A round uniform flat horizontal platform with radius and mass is mounted on frictionless bearings with a vertical axis at 0. At the instant of interest it is rotating counter clockwise (looking down) with angular velocity . A force in the plane with magnitude is applied at the perimeter at an angle of from the radial direction. The force is applied at a location that is from the fixed positive axis. At the instant of interest a small coin sits on a radial line that is an angle from the fixed positive axis (with mass much much smaller than ). Gravity presses it down, the platform holds it up, and friction (coefficient=) keeps it from sliding.
Find the biggest value of for which the coin does not slide in terms of some or all of , and . Answer: .
16.4.49 Frequently parents will build a tower of blocks for their children. Just as frequently, kids knock them down. In falling (even when they start to topple aligned), these towers invariably break in two (or more) pieces at some point along their length. Why does this breaking occur? What condition is satisfied at the point of the break? Will the stack bend towards or away from the floor after the break? Answer: The stack will bend away from the floor after the break. This breaking occurs because gravity is not able to provide equal angular momentum to all the falling blocks.