Volumen 1.1
Volumen 1.1
Volumen 1.1
McGRAW-HILL
i^^S?liiftk^\i
ENGINEERING
SCIENCi!^^''^^'^^
Digitized by the Internet Archive
in 2010
http://www.archive.org/details/principlesofmech01yehh
civD aaiUBB^
V(^
t^^4
Principles of Mechanics of Solids and Fluids
Frankel •
Principles of the Properties of Materials
Hodge •
Plastic Analysis of Structures
\ O L U M E 1
•IMFACafiaa^lBRARYf
Hsuan Yeh
PROFESSOR OF MECHANIC AL KXGINKKRIXG
DIRECTOR OF TOWNE SCHOOL OF CIVIL
AND MECHANICAL ENGINEERING
UNIVERSITY OF PENNSYLVANIA
Joel I. Abrams
ASSISTANT PROFESSOR OF CIVIL KNGINEKRING
YALE UNIVERSITY
Guar
^Wi :bi=~ tk :
Preface
prepared, since 1952, by the senior author for use in a junior course at
The Johns Hopkins University. The central theme of the book, as
originally conceived, may be st^ited as the integration of the mathe-
matical subjects which are otherwise abstract for many students \\ith the
realities of a mechanics course and the demonstration of the \miversality
(^f the mathematical expressions which portniy the physical phenomena
encountered in the several phases of mechanics. It is the authors' hope
that the completed l-K>ok fulfills these aims.
We have made an attempt to present the fundamentals of mechanics
in a thoroughly mature manner to persons \Aith no (or little) previous
training in the subject, other than a first course in general physics and
the completion of two ye^irs of college-level mathematics. In doing this
we are reaffirming our strong belief that the professional engineer of
tomorrow should receive an early training, in w hich the strongest empha-
sis should l>e placed on mathematics, fimdamental science, and other
level.
and problems form an integral j^art of the book, and careful consideration
should be extended to tluMU. At this point the instructor has ample
opportunity to discuss soni(> of the rainilicat ions and practical aj^plica-
Preface vii
Hsuan Yeh
JoelT. Abrams
Contents
Preface V
3-1. Introduction 29
3-2. Scalar Product 29
3-3. Vector Product 33
3-4. Physical Significance of Scalar Product; Work; Line Integrals ... 36
3-5. Physical Significance of Vector Product; Moment 39
3-6. Moment of a Vector 43
3-7. Velocity Due to Rotation; Angular-velocity Vector 45
3-8. Triple Vector Products 46
Problems 49
ix
X Contents
Vector Algebra
It willbe shown later in the text that there are other physical quantities
of a still more complicated nature than vectors, and they therefore require
more specifications than those of a vector. However, we shall not con-
cern ourselves with such quantities until we reach Chap. 11.
Vectors are often denoted in print by a Gothic or a boldface letter, A.
In manuscript they are usually denoted by an arrow placed above, as, for
example, A . A plain letter A or the vector enclosed by the absolute-value
sign, |A|, shall mean the magnitude of the vector A.
4 Vector algebra
in Sec. 1-2 and Fig. 1-3 next ; we draw the vector B from the terminal of A.
The new vector (in this case vector C) that joins the origin or beginning of
A with the terminal of B is the sum of A and B. It can be seen that the
same resultant C is obtained if we first draw B, then draw A at the terminal
of B, and finally join the beginning of B with the terminal of A. In either
case the sum C (represented by oc) of the two vectors A (represented by oa
or he) and B (represented by oh or ac)
is the diagonal of a parallelogram
(oach)formed by having A and B as
the two neighboring sides in the man-
ner shown in Fig. 1-3. In many engi-
neering applications we term the sum
Fig. 1-3
two vectors A and B the result-
of the
two vectors.
ant of the
The magnitude C of the vector C may be related to A and B, the magni-
tudes of the vectors A and B, by use of the cosine law of trigonometry:
where 6 is the angle between the positive directions of the two vectors A
and B. By suitable trigonometric relations, the direction of C in Fig. 1-3
may be easily determined.
From the above graphical construction, it is clear that the addition of
vectors is commutative; that is,
(A -h B) -f C = A -f (B + C) (1-2)
hurifJcoordinate Hy.stem .such thai if we rotaU; t,he jixi.s iril,o iJie yy ./; axi.s, a
ri^hl-hiUKJ screw would advance; alon^ the posit/ive dinsciion of ifie ^; axin.
liy r(;ver,sin}^ tlie ruh; of Uie {^rii-pliicM,! a<J(htioii of vcjfitorM, it (;aii bo soon
tFial/ a vector A can b(; (Jocoinf)OHed into any ruimbor of component vectorH
Sec. 1-8] Add il ion of veclors by the method of components 7
along the three coordinate axes. The projected lengths A^, Ay, A^
represent the three components of A, or as they are sometimes called, the
three scalar components of A. From the geometry of Fig. 1-6, it is
evident that
A = VAJ + /!,2 + A,2 (1-6)
A = AA + .4,j + AM (1-7)
The above equation expresses any vector A as a sum of its three orthogo-
nal components, or vector components, and one of the most frequent is
Ax — A cos a
Ay = A cos & (1-8)
A^, = A cos 7
The three cosines cos a, cos ^, and cos 7 are called direction cosines of A.
They are related by the following equation:
cos- a + cos- /3 -f cos- 7 = 1 (1-9)
A = ^^i Ay} + AM
-I-
B = 5,1 By'} BM
-1- -f-
where this plane intersects the x axis, say at point c, determines the projection of A
along the x axis. This projected length Oc is usually denoted by the symbol Ax. It
should further be noted that Ax is an algebraic scalar quantity.
8 Vector algebra
C = C.i + C,j + CM
then Ci = Ax + x>x
Ly = Ay -\- By (1-11)
C. = A, + B,
We have thus proved a very important fact: The component along any
direction of the sum of vectors is equal to the sum of the components of the
vectors along that direction. Although Eqs. (1-11) are derived for the sum
of onlytwo vectors, they can be easily
extended to the sum of any number of
vectors.
The above fact can be conveniently
used to calculate the resultant (i.e.,
sum) of any number of vectors. The
procedure is decom-
as follows : First
pose each vector into its components
along the x, y, and z axes; next form
the algebraic sum of all the scalar com-
ponents along X and similarly for all
components along y and z; and then
finally form the resultant by vector
addition of the three resultant vec-
tor components. This is called the
method of components.
It was pointed out in Sec. 1-3 that
two vectors are equal if, and only if,
they have the same magnitude and
direction. Certainly, then, expressed
from the point of view of components,
Fig. 1-7
two vectors are equal if, and only if,
they have identical scalar components in any three noncoplanar direc-
tions.! The converse of this statement is obviously true also. From this
it follows that if a vector is equal to zero, each and every one of its scalar
components must necessarily be equal to zero. In general, then, a vector
equation A = is equivalent to the three scalar equations
A, = Ay = .4, =
Example 1-2. Given two forces A and B having magnitudes
A = 2 -\/lO lb and B = 15 lb, respectively, and directions as indicated
in Fig. 1-7. Find the resultant force R = (A -1- B).
A, = -^
Vio
X 2 VTO = 2 A, = -^
Vio
X 2 VTO = 6 .1, =
fix = By = 3.^ X 15 = 9 B, = +5 X 15 = 12
A = 2i + 6j + Ok lb
B = Oi 4- 9j -I- 12k lb
and the direction cosines of the resultant force are, from Eqs. (1-8),
2 15 12
cos a = cos
,
/3
= cos 7
\/373 -\/373 V'373
PROBLEMS
1-1. ^^'hat is the significanco of A/,l?
1-2. Show that |A| + 1B| > \\ + B\.
1-3. Determine the sum and difference of the two vectors A and B:
A = 6i - 4j - 6k
B = 4i - 2j - 8k
A = 4i + 5j + 3k
B = 4i - 5j + 3k
C = 3i + 5j - 4k
Express Fi and F2 in the form of Eq. (1-7). Also find the magnitude and the direction
cosines of
(a) Fi + F2.
(b) F, - F2.
(c) F2 - Fj.
Fig. P 1-6
Fui. P 1-11
Problems 11
of the piston exerts a horizontal force of 300 lb on the right end of tlie connecting rod.
Find the components of this force along the connecting rod and along a direction
perpendicular to the connecting rod.
1-12. An airplane climbs at a speed of 200 mph and at an angle of climb of 30°, both
relative to the air. If there is a vertically upwartl gust of 30 mph, what is the absolute
velocity of the airplane and its absolute angle of climb? (By absolute, we mean here
relative to the earth.)
200 mph
Fin. P 1-12
1-13. If n is a unit vector with direction cosines Ox, ay, and a., show that its scalar
components are simply equal to its direction cosines. Express n incomponent form.
1-14. A unit vector makes equal angles with each of the three coordinate axes.
Express this vector in terms of its vector components.
1-15. What is the locus of the extremity of the vector formed by A + A:B? Write
the equation of this locus if A and B are given in the form of Eq. (1-7) and A; is a real
number between — =o and + « (Assume that the origin of the above vector sum
.
F= (mv) (2-1)
J
where F = force
V = velocity
m = mass
mv = momentum
For an ordinary body, the mass is a constant, and we can write, alterna-
tively,
F = ma (2-2)
12
Sec. 2-2] Discussion of iS'ewlons laws 13
Although Newton's laws are familiar to readers of this text, there are
scA'eral aspects of these laws that need to be emphasized again. Past
experience shows that the following points are often not sufficiently clear.
1. A bod}' can be said to be in motion or at rest only when compared
internal forces and their influences do cancel out. On the other hand, if
we take only one of the two particles as our system, then only one of the
two forces is present.
The paragraph also illustrates the importance of clearly defining
last
one's system. Picking a system (or in the language of applied mechanics,
choosing afree body) is a prerequisite to the solution of almost all problems
in mechanics and sometimes requires rather careful thinking. Once a
free body (at present the particle) is chosen, all the forces that act on that
particle should be clearly indicated. The forces that need be considered
generally are the gravity force, or the weight of the body which is usually
assumed to act vertically downward (or more accurately, but still not
exactly, toward the mass center of the earth), all reactions and friction
forces due to the constraining of the particle, and finally all other forces
due to the mutual action of neighboring bodies and external agencies.
This concept will be illustrated by a number of examples in subsequent
discussions.
Clearly, what we choose as our free body depends to a large extent on
the idealizations we employ. The depend on the
idealizations themselves
type and accuracy of the information we desire. For example, if we are
interested in the motion of the earth around the sun, we might idealize
the earth as a particle traveling around another particle, the sun, fixed in
space. On the other hand, if we were interested in earthquake phenom-
ena, we might idealize the earth as an elastic sphere, fixed in space.
Neither of these models or idealizations exactly depicts the earth; the
validity of the idealization must be confirmed or rejected by experiment.
We are taking the view in this text that Newton's three laws are
fundamental postulates and thus serve as our basis for the study of
mechanics. Certainly this does not imply that these laws will remain in
their present or unaltered form forever. A little reflection on their
comparative newness in the historical sense would indicate that it
required many centuries for the enumeration of these principles in their
present form, and it is entirely possible that in centuries to come addi-
tional knowledge in the field of mechanics could lead to their modifica-
tion. Any future modification of such basic postulates is more than
likely to include the present laws as a particular case. This is so because
the present laws seem to apply so very well to most situations; therefore
one would expect any future discovery to be generalization rather than
contradiction.
Our discussion, then, will be centered around the use and interpretation
of these basic laws. Any discussion relative to the more philosophical
aspects of these laws will l)e omitted. Sufiico to say, then, that the
physical quantities, specifically force and mass, related by the laws exist
and that our knowledge of them is intuitive. Velocity and acceleration,
Sec. 2-3] Sialics and Ihe equations of equilibrium of a parlicle 15
of course, will he defined in terms of the spatial motion of a point, that is,
only if, the total, or resultant, force acting on the particle is ecjual to zero.
It follows, then, that a particle is in equilibrium when the resultant of the
forces acting on the particle is zero.
Thus, if Fi, F2, . . . ,
F„ are the forces acting on the particle, we may
express analytically the neces.sary and sufficient conditions for the etjui-
librium of a particle as
2 Fi = F, + F, + • • •
+ F„ = (2-3)
t = i
n
F, = 2F = (2-4)
where F, is the resultant force, or sum of all forces, acting on the particle.
The above vector equation can also be expressed in the form of three
scalar equations if we deal with the components of the vectors. Let us
write, for example,
Fi = F,A + F,J -I- F,,k
Then Eq. (2-3) or (2-4) is eciuivalent to the following three scalar ecjua-
tions, which, following Sec. 1-8, are both necessary and sufficient condi-
tions for the equilibrium of a particle:
sarily be equal in magnitude and opposite in direction and have the same
line of action, or be coUnear. On the other hand, if a particle in equilib-
rium is acted upon by three forces, the force polygon becomes a force
triangle, and the three forces must lie in one plane; i.e., they must be
coplanar.
Example 2-1. An electric lamp of weight W
is suspended at the end
of a string AB, the upper end of which is attached to the ceiling by three
Fig. 2-1
-ymmetry it is clear that the forces in the three strings AC, AD, and AE
are the same. Let us denote them by T2. Because of the ecjuiUbrium of
the joint A, we have, according to Eqs. (2-5),
where we have used the fact that the cosine of the angle between the
three strings and the vertical is
1 - —
V 3/2
3 Vl - a^SP
It should be noticed that in this problem the physical consideration
enabling us to obtain a solution was the fact that the entire system, in this
instance composed of the lamp and all the connecting strings, is in a state
of static equilibrium. Here we have tacitly assumed that
system in for a
equilibrium each and ever}' particle comprising that sj^stem must also be
in equilibrium. That this assumption is valid can be easily verified by a
consideration of the definition of equilibrium. The directions of the
forces in the strings are indicated in Fig. 2-lc. Remembering that Fig.
2-lc indicates the action of the strings on the particle .1, the particle -I in
turn reacts in an equal and opposite fashion on the strings. Therefore
the force acting on the string tends to elongate the string. We term this
type of force a tensile force; the string is in tension. Forces that tend to
contract or shorten a body are known as compressive forces; the body is in
compression. One characteristic of a flexible string is that it can sustain
only tensile forces.
can sustain both compressive and tensile forces. Rigid bars are therefore
essential elements of all engineering structures. In the analysis of
structures, it is usually permissible to assume that the bars are hinged
without friction at their two ends. In actual practice the bars may be
riveted or welded together, often with addition of gusset plates; but the
assumption of frictionless, hinged joints gives a sufficiently good answer as
a preliminary analysis, requiring only a fraction of the labor needed for a
more detailed analysis. Bars used in structures are usually slender in
cross section when compared Avith their lengths and are usually considered
iy
(c)
(6)
Fio. 2-3
\/2 IV sin 7
Fi= - — cos 7
sni 7
F, =
sm 7 — cos y
20 Newton's laws; concurrent force systems; equilibrium
The first and last equations can be solved simultaneously for F^ and F4.
Combined with the expression iorFo obtained before, we find the value of
Fz and Fi in terms of the weight W.
For example, if a = /?, then
W
2 cos 40° cos a (I — cot 7)
After we find F3 and Fi, the equation of equilibrium 2Fy = jdelds the
reaction F5.
Remarks. In the case of joint C, all the forces lie in one plane. This is
called a plane problem. In a plane problem only two of the three equa-
tions of equilibrium are nontrivial; hence we are able to solve for two
unknowns, magnitudes of the two forces (Fi and Fo). In
in this case the
the case of joint B the forces do not all lie in one plane. This is a space,
or three-dimensional, problem. In a space problem we have at our dis-
posal all three ecjuations of e(iuiHbrium; hence we are able to solve for
three unknowns, in this case the magnitudes of F3, F4, and Fa.
We also observe from the preceding solution that the initial choice of
sense of Fi and F2 is completely arbitrary. Whether this choice is correct
or not will be completely determined by the ensuing solution, a positive
rcsulf indicating the original selection was correct. However, once this
choice has been made, it is iinp(M-at i\(> that in all other free bodies that may
Sec. 2-5] Friclion 21
contain this force, the direction he consistent with the initial choice.
Remenil)ering that Fi. for example, acting at joint C in Fig. 2-3^, repre-
sents the action of member CA on -joint C and that in turn joint C,
according to Newton's third law, reacts on member CA in an efjual and
opposite fashion, we see that joint C "pushes" on CA, and thus CA is in
compression. It is usually not wise to attempt to guess beforehand
whether the force in the various members is tensile or compressive. A
better approach would be to solve the problem in a straightforward
manner, assuming that all the forces are tensile, and then subsequently
checking the results with physical intuition. In this fashion a sign
convention is established, positive indicating tension, negative indicating
compression.
It can be seen that if we had started by taking joint B instead of
joint C as our free body, we would have too many unknown magnitudes
{F i,P z,F i,F h) in our three equations of equilibrium. The proper sequence
of taking free bodies in a problem is therefore of importance if we wish
to solve a problem in the most convenient fashion. Occasionally a situ-
ation arises such that, regardless of the sequence of isolating joints, there
are more unknowns than there are available e<iuati()ns of e(iuilibrium.
Obviousl}^, then, a unique solution is not possible. We term these prob-
lems statically indeterminate (not de-
terminable by methods of statics
alone). They require for their solu-
tion additional information that
usually can be obtained by a con-
sideration of the deformation of
the structure caused by the various
forces.
2-5. Friction
point of contact or even a single line of contact, but rather may have a
surface of contact. The contact forces, then, may be distributed over
some finite area. In this event the normal force and friction force we are
considering here are simply the resultant of the distributed forces and
represent a convenient idealization that enables us to achieve a reasonable
solution for problems involving friction. More will be said of distributed
forces in later chapters (Chaps. 5 and 15).
F = iXsN (2-6)
The coefficient of static friction between the block and the inclined plane
l"iG. 2-5
is n,, and we assume that ju., > tan </>. (Unless Hs > tan 0. the block will
not stay at rest on the inclined plane.) Find (a) the horizontal force P
(applied to the center of the block in a plane parallel to the inclined plane)
which will just cause the block to begin to slip, and (h) the friction force F
if P than the value required to start the motion.
is less
Solution. We take the block as the free body and consider its cciuilib-
rium (Fig. 2-56). The forces acting on the block, referred to the coordi-
nate axes shown, are
W = IF sin i - IF cos j
N = iVj
P = Pk
F = FJ + /'\k
^F, = -W cos + A^ =
</) N = W cos <t>
ZF, = P + F, = F, = -P
Hence F = \/FJT~F? = VlF^ .sin^ + P^
24 Newton s laws; concurrent Jorce systems; equilibrium
(6) If P is less than the above value, then slipping is not imminent, and
F = VW^ sin2 + P2
where P can be any value smaller than given in (a) . It can be seen that
here F ?^ iXsN and that F < jisN.
one component part of the given system. A suitable free body is one
which is acted on by only two concurrent forces of known direction but
unknown magnitudes in a plane problem or by only three concurrent
forces of known direction but unknown magnitudes in a space problem.
Alternatively, in a plane problem the unknowns could consist of one
unknown magnitude and one unknown direction angle, and in a space
problem two unknown magnitudes and one unknown direction angle or
one unknown magnitude and two unknown direction angles, etc. The
essential feature is that there be a total of tivo unknowns in a plane
problem and a total of three unknowns in a space problem. (The reason
for such a choice was given in Sec. 2-4.)
2. Sketch the free body and show all the forces acting on it. The forces
will generally include all the known, or active, forces as well as all the
(such as the example of Sec. 2-4), repeat the above procedure for the next
free body system is solved.
until the entire
The problems end of this chapter are examples of a great variety
at the
of problems that can be completely sohed bj- following through this
procedure.
PROBLEMS
2-1. .\. "particlp" at .1 is .suspended by two strings, .1/^ and AC. If the particle
weighs 100 lb, find the force in each of the strings.
Fi
tlOOIb
Fu;. V 2-1 Fig. P 2-2
2-2. Three coplanar forces F,, Fa,and F3 are acting on a particle in equilil)rinm.
If Fiand a are known quantities, determine F2, F3, and such that F3 is a niininium.
2-3. A simple structure consists of a bar AB and an ideal string CB pin-connected
at their ends as shown. This structure carries a vertical force of 1,000 lb at end B.
Find the forces in each of the members.
1 1000 lb
2-4. A small bead of weight W can slide on a vertical circular wire hoop of radius a.
The bead is attached to point yl by a and spring con-
linear spring of initial length a
stant k. The coefficient of static friction between the bead and the wire is pi,. Find
the position of equilibrium. By spring con.stant A- we mean that the force in the spring
is given as F = k8, where 5 is the change in length of the spring.
26 Newton s laws; concurrent force systems; equilibrium
2-5. Verify the results of Example 2-1 by the use of all three of the equations of
equilibrium for the particle at A.
2-6. A weight W
is supported by the structure shown consisting of three hinged
bars AB, AC, and AD. Find the reactions in each of the bars.
Qw
2-7. A truss in the shape of a tetrahedron and consisting of six bars of equal length
is loaded by a vertical force P as shown, with the base triangle horizontal. Find the
reactions in all bars.
2-8. A
simple structure consists of three hinged bars AB, AC, and AD
oi equal
length 2a. Find the reaction in the bars due to a vertical load W
applied at the
joint A.
Fui. P 2-8
Problems 27
2-9. In the truss shown, A BCD and A'B'C'D' are two squares witli sides of lengths
(I and Ad, resiiectively. The vertical distance between the two squares is 2«. Find
the forces in the bars .1/^, H(\ HIV, C(", Cli', CD' due to tiie four vertical loads P.
Fic. V 2-9
2-10. Twelve bars are hinged togetlier to form a hexagon. One of the radial bars
isecjuipped with a turnbuckle so that the tension in the bar can be set at any value,
say T. Find the tension (or compression) in all bars.
Fui. r 2-10
2-11. Tworollers A and li are hinged at their centers to a rigid weightless bar .1/^
The rollerscan move without friction along the two inclined wails as shown. If the
weight of .1 is 10 lb and the weight of li is 40 lb, find the angle a between AB and the
iiorizoiital for the system to be in equilibrium.
28 Newton s laws; concurrent force systems; equilibrium
Fig. P 2-11
2-12. A circular cylinder A rests on top of two half circular cylinders B and C, all
having the same radius r. TF and that of B and C is W/2 each.
The weight of A is
Assume that the coefficient of static friction between the flat surfaces of the half
cylinders and the horizontal table top is /xs- Determine the maximum distance d
between the centers of the half cylinders to maintain equilibrium.
V ////////// ^/ 77/77
Fig. P 2-12
2-13. What horizontal force P must be applied to the weightless wedge in order to
raise the 1,000-lb weight? Assume that the coefficient of static friction between all
Vector Products
3-1. Introduction
hence also the term dot product) is a scalar algebraic number having the
t Not only are vector products of importance in mechanics, hut they have great
significance in other physical sciences such as optics, electricity, and magnetism.
29
30 Vector products
ij=jk = k.i =
Ji=jj =kk=l (3-i;
Fk!. 3-1
Since, as illustrated in Fig. 3-2, A cos 6 is
Fkj. 3-2
gcoiiicl ric.'illy as
A •
it .|(/) cos 0) . 1 (coiiipoiiciit (if It oil A)
or A •
It li{.\ COS 0) - -
/>(C0111|)011("llt ol' A (til It)
Plxprcsscd in words, this result may bestatiHl as follows: Tlu^ scalar prod-
uct of two \(>clors is \\\v pro(hict of the inngnitudc of I he lirst times the
magnifiidc of the second project(Ml on the lirst, or \i('(> \-ersa.
:
Fig. 3-;i
Commutative: A •B = B A •
A = .l,i + -4 J + -l.k
B = BA+ BJ + BM
Then. A- B = (.4,i -H-4 J
-H-4,k) -{BA -hBJ +B.-k) c
Since the dot product obeys the dis-
tributive law, we can expand the
right-hand side and obtain, with the
help of Eqs. (3-1).
Fin. .-^4
A B= • .4 ,B, -\- A ,B, -h .4 ,/^, (3-3)
Hence the scalar product of two rectors is equal to the sum of the products of
their respective orthogonal components. In particular, we see that
A-A = .4,24.^4^24.^4^2 = ^42,
Example 3-1. Let C = A B -|- (Fig. 3-4). Then
C C = (A
•
-h B) • (A + B) = A • A + 2A • B -j- B B
•
32 Vector products
A = 3i - 2j + 6k
B = 4i + Oj - 3k (Fig. 3-5)
Then A B =
• 3 •
4 + (-2) + •
6 •
(-3) =
Component of A on B = —^— = —^
Component of B on A = —— ^ "^ ~7~
cos = AB = -^
-6
e = 99°52'
A B 35
the plane determined by the three points. Let the directed line segments
be, ab, and ac be denoted by the vectors A, B, and C, respectively. Then
A - -3i - + j 4k
B - -i -F 2j + k
C = -4i +j + 5k
= B- 3 - 2 -f 4
be
'^
\/<) + 1 +16 V26
and finally, the perpendicular distance from point a to line t)c is
A X B (B X A)
Hence the vector product is not commutative, and the order of forming the
vector product is of great importance.
Also from the definition, it can be seen that the magnitude of A X JJ is
equal to the area of the parallelogram having A and B as its conterminous
sides (Fig. 3-7). This geometrical interpretation often proves to be
convenient.
t A more
convenient application of the right-hand-screw nile would be to let the
hand curl in the direction of the indicated rotation; then the thumb
fingers of the right
of the right hand points in the direction of the positive sen.se of the cross product.
34 Vector products
Fig. 3-8
Furthermore, the above three vector ijrodiicts lie in the i)lanc Oah, since
they are all perpendicular to C. Viewed in the plane Oab (Fig. o-Sb), we
find that the magnitude of the three vector products is proportional to the
lengths of the three sides of the triangle Oab and that the orientation of
these vector pi-odnds is iioniinl to Ihc corresponding sides of the triangle
Sec. 3-3] Vector product 35
Oab. The three vector products therefore form a closed triangle O'a'b'
(Fig. 3-8c) which is similar to the triangle Oah but rotated 90°. Hence
C X (A + B) = C X A + C: X B
and Eq. (3-4) is proved.
Applying the definition of cross product to the orthogonal miit vectors
i.j. k,
iXi=jX.j = kxk =
i X j = -j X i = k
(3-5)
j X k = -k X j = i
kXi= -iXk=j
Equations (3-5) lead us to the calculation of the cross product in terms
of their components:
1 i
k
A X B -3 4 -9i+j - 7k
-4 5
AB sin d A X B VT31
|A X B
Therefore d =B sin d =
A
= V 13L 26
The student should verify the fact that the same result would have beeu
obtained had we used the product A X C rather than A X B. This
result can be generalized if we let r represent a vector from a point P to
any arbitrary point along the line of action of a vector A. Then the
shortest distance from the point P to the line of action of A is given as
= X A[
=
d r sm
A
That r may have its terminal anywhere along the line of action of vector A
is amply demonstrated by Prob. 3-15.
d =
C
In physics avc learn that the work done by a force acting on a particle
is defined as the produc^t of the ('omi)onent of the force in the (.lirectiou
of the displacement of the particle and the displacement. In Fig. 3-10,
V is a varial)le forcef acting on a particle P which travels along a spatial
is defined as
d\V = Fidr) cos d
where 9 is the angle between F and dr. We see, however, that the Avork
d\V thus defined is just the scalar product (and conse(iuently a scalar
Fic. 3-10
fjuantitAM
dW = F . dr = Fx dx -f F„ dy -H F, dz
38 Vector products
W = r^ F 'dr =
jl^^'
F, dx +Fydy + F, dz (3-8)
Thus we note that the total work done by the force F while moving from
Pi to Po is simply the sum of the work done by each component of the
force
W = TF. + T7, + W, d
where Wx = \ F^dx etc.
J XI
^'F-rfr (3-9)
/.
is called the tine integral of F. The line integral is taken along a specified
path between the two points Pi and Po and is itself a scalar quantity.
Certainly, the line integral, in general, depends upon the path as well as
the end points of the path. A conservative force field is a particular type of
force field defined such that the work done does not depend on the path,
but rather depends onhj on the end points.
The determination of the total work done by a force F moving between
the points Pi andPo, or the evaluation of the associated line integral, may
be accomplished by the use of Eq. (3-8). Suppose that the path traveled
is expressed in parametric form, such that
''
=
^'(''
I + ^"('' t ^•(•^)
I ds (3-10)
I
so that the line integral has been converted into an ordinary type of
integral.
Example 3-7. Determine the work done by the force F = y-i — x^j
as it travels along the parabolic path y = x- from Pi(0,0) to P2(2,4).
Letting x = s, we have dx = ds, dy = 2s ds, and
F = .s^i - .s'^
This may he achieved since the path is fixed, tlius establishing a definite
relationship between the variables. Choosing, for example, .r as our
variable, we have
W =
f^'
{X' -2x^)dx = -«5
The student should verify tliis result by selecting y as the one variable of
integration.
FCQSd
: .
40 Vector products
the bolt, then Fig. 3-1 la is a view in a plane parallel to the xy plane.
M, = IF sin 6 (3-11)
where the subscript z denotes the axis about which the moment is taken.
It can be seen in Fig. 3-1 la that sin 6 represents the perpendicular
/
distance between the z axis and the F vector (which is in a plane normal
its extensions; this distance is called the moment arm of F.
to the z axis) or
Thus we may consider the moment of a force about an axis as the product
of the force and its moment arm. Furthermore, the value of M^ given
by Eq. (3-11) can also be considered as the sum of the moments of the
two components of F, namely, Fx and Fy. In other words, we can also
write
M, = xFy - ijFx (3-12)
where x and y, the components of the distance I, are the moment arms of
Fy and Fx about the z axis, respectively (Fig. 3-116). That the two
expressions for M^ as given by Eqs. (3-11) and (3-12) are identical can
readily be proved as follows
the X axis. Analogous to Eq. (3-12), we can form the moments about the
x and the y axes as follows:
Mx = yFz - zFy
My = zFx - xF, (3-13)
M, = xF„ - yFx
The last of Eqs. (3-13) is merely a repetition of E(|. (3-12). In these
expressions for niomenis, W(> ha\-e used llic lii^lil-lumd-scrcw rule: The
1
Fk;. 3-12
M - MA + M,j + MM
Substituting Eqs. (3-18) into the above expression, we find
i J k
M
42 Vector products
product of F and its moment arm. We see that this is in complete agree-
ment with our usual concept of "moment" as given in introductory
physics, t We also see that the axis of rotation due to the action of F
about the point is just the direction of the newly defined INI vector.
In particular, if we let the line L coincide, respectively, with the .r, y, and
z axes, we determine immediately that the components of projected in M
these directions arc the same as those given by Eqs. (3-13). For example,
Thus we have that the moments of a force about the three orthogonal
axes passing through a point are simply the scalar components of the
t See, for example, F. W. Sears and I\I. \\'. Zemansky, "College Physics," p. 52,
Acldison-We.sley PiiMishinf!; Company, Readiiifj;, Mass., l'.>52.
Sec-. 3-6] MonienI of a vector 43
moment vector at
point. In many practical applications the scalar
tliat
methods [Eqs. nuiy offer a more direct solution than the more
(.3-1^^)1
as
M = r X F
where r is the position vector of P relative to 0. The direction of M is the
axis al)out which there is a tendency for rotation to occur because of the
action of F, and the magnitude of M is a measure of the effectiveness of
turning about that axis.
2. Moment about an Axis. The component of the above-defined M
along any axis through is the moment of F about that axis passing
through [FCq. (8-1 3a)]. The sense of the moment is determined l)y tiie
right-hand-screw rule.
Example 3-8. A force vector F = i -|- 3j — k acts at a point
P{S, —I, 4). Find the moment of F about an axis passing through the
point Q(l,l,2) and having the direction of the \ector A = 4i -f- 3j -\- Ok.
Since the vector r from Q to P is 2i — 2j + 2k, the moment of F about
point (J is
i
J k
= X F - -2 = -4i
I
AI r 2 2 -f 4j + 8k
I
1 3-1
Hence the moment of F about the axis thiough (J along A is
M A ^ 44
A ~ V<iT
the property of a vector product. Any sliding of the vector A (Fig. 3-15)
would result in a change of the position vector from r to r + A;A. The
resulting moment is, by virtue of the fact that A X A = 0,
(r + AA) X A = r X A
Hence the moment of a vector is unchanged by sliding the vector along
its own line of action. When the vector in question is a force, this con-
cept sometimes referred to as the transmissihility of a force. It also
is
M= ri X A + r.> X (-A)
Because the moment is~ unchanged by a sliding of the vectors, this can
also be written as
M ri X A +.r, X (-A) = ri X (A A) =
Hence the vector sum of moments about an}' point of two equal, opposite,
and colinear vectors is equal to zero. This statement, which is almost
intuitive in nature, assumes importance when we consider the forces of
interaction between particles, which, because of Newton's third law, are
equal, opposite, and colinear.
Now suppose we have a system of vectors acting at a point P (Fig. 3-1 7).
Si'c. 3-7] Vehcily due lo rotation: angular-velocity vector 45
M = r X F, + r X F, + + r X F„
.AI = r X (Fi + F, + • • •
+ FJ
M = r X F,
Fic;. 3-1
Let us select any point on g> and let the position vector of P relative to
be represented by r. Then the velocity v of P due to rotation is given as
V = to X r (3-17)
In order to show that this definition agrees with our previously described
concept of velocity due to rotation, we observe that to X r is a vector having
direction perpendicular to r and w and a sense as determined by the right-
hand-screw rule and therefore has the same direction and sense as
previously described. Moreover, we have that
Thus the two vectors defining the velocity of point P have the same
magnitude and direction and by definition are equal. Equation (3-17)
isthereby verified.
According to Sec. 3-6, we can interpret the velocity v due to rotation
about a fixed line as the negative of the moment of to about the point P.
From this it follows that the location of point along w is immaterial, for
a new location 0' along o> merely changes r to r -1- ^w and therefore does
not modify v.
A • (B X C) = A
Sec. 3-8] Triple vector products 47
is as follows. As shown
magnitude of (IJ X C) is equal to
in Sec. o-o. the
the area of the parallelogram having B and C as its conterminous sides
and the direction of (B X C) is perpendicular to the plane of B and C (Fig.
3-19). Since A (B X C) has a value ecjual to the product of the
•
A X (B X C)
A X (B X C) = aB -f /3C
Hence B X C = 5.C^k
(A . C)B - (A B)C •
since A C =
• A^Cx + AyCy
and A B =
• A^B,
Hence A X (B X C) = (A • C)B - (A B)C (3-21)
The form of the above expressions shows clearly that the vector triple
product does not remain the same by any interchange of its elements.
Example 3-9. We have already seen evidence of the utility of triple
vector products where in Sec. 3-5 we used the scalar triple product to
b
Fig. 3-20
define the component of a moment vector along the direction of some line
L. As a finalexample of triple vector products we shall consider two
geometric problems.
First, let us determine the piercing point (point of intersection) of a
lineand a plane. Let the line ah intersect the plane at point P (Fig. 3-20).
Suppose we choose any three points, c, d, and e, lying in the plane.
Denote the directed line segments ab, aP, ac, cP, ce, and cd by the vectors
A, B, C, P, E, and D, respectively. Then, because the three vectors P,
E, and D are coplanar, we have
P . (E X D) =
Further, we can write that
P = B - C = ;//A - C
where m is some scalai' luiinlxM- (o be determined. Substituting this into
the previous expression,
(;/;A - C) . (E X n) =
Prohlenis 19
_ C • (E X D)
"" ~ A . (E' X D)
which can readily be evaluated according to Eq. (3-18). Once the value
of m is known. P can he readily determined and in turn the point of inter-
section P can be located.
As a second illustration, determine the vector projection of vector A on
some plane defined by its normal n (Prob. 3-10). The vector projection
of A on the normal n is, from Sec. 3-2,
(n . A)n
Now. considering the triangle formed
by the vector A, its vector projection on
the normal, and its vector projection B
on the plane itself (Fig. 3-21), we im-
mediately^ determine the vector projec-
tion on the plane as
B = A - (n- A)n
Using Eq. (3-21), this becomes
B = n X (A X n) Fig. 3-21
PROBLEMS
If d is the angle between two vectors A and B, derive an expression for ros 6
3-1.
in terms of the components of A and B.
3-2. Two line segments OP and OQ are formed by joining the origin O with the
points P(3, —5,2) and Q{ — 1,2,0), respectively. Find (a) the projection of OP on OQ,
(6) the angle between OP and OQ, (c) the unit vector perpendicular to the plane of
OP and OQ.
3-3. Given the force vector F whose components are
Fx = 3 lb Fy = -3 lb F^ = 1 lb
3-5. Prove that any two lines drawn from each extremity of the diameter of a circle
to a point on the circumference are mutually perpendicular.
3-6. Find the angle between the body diagonal and one edge of a unit cube, both
intersecting at a point.
3-7. Find the shortest distance from the line joining points a(3,l, — 1) and />(2,3,0)
to the origin O.
3-8. Show that A X A = by expanding the cross product in terms of its
components.
3-9. Using the method of Sec. 3-3, prove that the vector product is distributive for
the special case of three coplanar vectors.
50 Vector products
n = -4i + Bj + Ck
Hint: A normal to a plane is perpendicular to any two nonparallel lines in the plane.
3-11. IfC is normal to the vectors A and B, show that C is normal to A + B and
A — B. Whsbt does this result indicate about lines perpendicular to planes?
3-12. A force vector field is represented by
where F is the force acting on a particle located at point {x,y,z). Find the amount of
work done on the particle when the particle is moved along a straight line from the
origin to the point (2,1, —2).
3-13. A force vector F = 3i + j — k acts on the point (3,-2,3). Find the
moment of F about an axis through the origin and having the direction of the line
x-2,=2y + l=z.
3-14. A bar OP is 10 units long. a vertical wall through two guy
It is fastened to
wires at S and P and a ball joint at 0. The load W
is attached to the bar at T.
Assuming that the reactions in the two guy wires at S and P are Fi and F2, respec-
tively, find the moments due to the forces about the ball joint 0; i.e., find (a) moment
Ml due to Fi in terms of F,, {h) moment M2 due to Fo in terms of Ft, (c) moment M3
due to W.
Fig. P 3-14
Prohlems 51
3-13. Prove tliat the cross product of a vector A and any other vector that tcrini-
iiates on a hue jiarallel to A is constant; tliat is, show that A X H = A X C. This
ri'sult will lie extrenielv useful to us.
Fi(,. P :M5
3-16. A body rotates about an axis x = // = r with a constant angular speed
ri^id
of 1,000 rpm. Find the velocity and its components along x, y, and z axes of a parti-
cle at (1 in., —1 in., 3 in.) on the body.
3-17. Prove (A X B) (C X D) = (A C)(B D) - (A n)(B C).
• • • • •
3-18. Find the shortest distance between the two nonintersecting lines, one of
which is the body diagonal and tiie other of which is a face diagonal of a unit cube.
3-19. The force field of Prob. 3-12 acts on a i)article as it niak<'s one circuit around
a square in the xi/ plane of sides of length a. Find the work done on the particle.
Is the force field conservative?
CHAPTER 4
F, = 2F = (2-4)
F, + f, = (4-1)
n H n
Xow l)vXewton's third law of action and reaction, ail the internal forces
occur in e<iualand opposite pains, and thus we see that all the internal
forces must atld up to zero. Hence Ecj. (4-2) reduces to
2 F, =
i= 1
where F^ = I'F is the sum, or resultant, of all the external forces acting on
the system.
We now consider the moments of the forces. Let be any fixed point
in space and let the position \ector from to the ^th particle be r,. From
Eqs. (4-1) and (3-4), we obtain
r, X (F, + f,) - r, X F, + r, X f, =
We recognize that each of the two terms in the ai)o\c equation is a
moment (cf. Sec. 3-5). Adding n such equations for all /; particles, we
obtain
n n
Xow the internal forces f^ always occur in equal, opposite, and colinear
pairs. Let us, for call them fi2 and f^i,
example, pick one such pair and
so that fi2 is the force acting on particle
due to particle 2 and f-ii is the 1
The last relation is true because fi2 and £21 are colinear and act along the
line joining the two particles, although ri 7^ Vo (Sec. 3-G). The total
moment about (), of all the internal forces, is the sum of all such pairs of
moments and is therefore zero. Hence
y r, X f ,
= (4-5)
^ r. X F. = (4-6)
54 Statics of systems of particles and of rigid bodies
TVI,. = 2M = (4-7)
2F = 2IM = (4-8)
We can also deal "oith the components of F and M; the above two
vector equations are equivalent to the following six scalar equations:
//////^ ////yf
\y
TTT^^^TTTT TTTTTTT^T-
FiG. 4-1
still allows the body to rotate about any line through .4. The ball joint
at .4 is termed a reaction and accommodates anj' force necessarj- to
a frictionless hinged bar which offers a force onl^' along one direction,
namely, along the axis of the bar. If the axis of the bar at C does not
intersect AB. rotation about AB is prevented and the body is completely
56 Statics of systems of particles and of rigid bodies
SF = SM = (4-8)
A = A A + .1,J + .l,k
B = Byj+ BM
C= C,j
S(H'. 1-1] Procedures of solving problems in sialics of riqid bodies 57
:^¥ + 2F' =
i:M + 2M' =
Because the original system satisfies Eqs. (4-8),
vF.= o = A-|-B-fC
vM' = = Cy{AD)x - B,{AB)i -f B„(.4B)k
These vector ecjuations reduce to six scalar equations whose only solutions
are
Ax = Ay = A^ — By = B; — Cy =
We leave it Therefore the imposed con-
to the student to verify this.
straints offerno resistance and the body is in eciuilibrium without any
reactions at the supports A, B, and C. The supports can therefore be
removed without any effect on the state of the body. This proves that
the Eqs. (4-8), or their equivalent Eqs. (4-9), are both necessary and
sufficient for the equilibrium of a rigid body.
It now becomes apparent that the scope of statics is not unduly
restricted by the fact that Newton's laws are applicable, in a narrow
sense, to particles alone. Equations (4-8) or (4-9) indicate that a body of
finite size in equilibrium may be treated as a particle provided 2M = 0,
which in turn may be achieved if the forces are all concurrent. The tacit
assumptions made in Chap. 2 are thereby justified, and the scope of
particle mechanics is thus increased from a mathematical idealization to a
much more accurate physical reality.
The procedure outlined in Sec. 2-0 for solving problems of the statics of
a particle or a joint applies etiually well to statics of rigid bodies in general.
Only one modification is recjuired. In the case of a particle or a joint,
the lines of action of all the forces intersect in one point (i.e., they are
concurrent), which is the particle or the joint in (juestion. With respect
to this point, none of the forces contributes any moment. The moment
equation is identically zero, and hence does not yield any information.
58 Statics of systems of particles and of rigid bodies
TTTTt^TTT"
(a) (6)
Fig. 4-2
Following this reasoning, we see that a reaction that lies in a given plane,
but is otherwise unknown in magnitude and in direction, is equivalent to
two unknowns. Further, a reaction of unknown magnitude but known
direction constitutes one unknown.
If we have a problem in which all the acting forces lie in one plane, say
the xy plane, our useful equations are reduced from the six of Eqs. (4-9) to
only three:
^ZV., = :^F, = ^M, = (4-10)
Hence, for every plane free-body problem, we are generally able to solve
for three unknowns. The complete constraint of a rigid body in a plane
by means of three hinged bars is shown in Fig. 4-2a. The two hinged
Sec. 4-1] Procedures of solving problems in sialics of rigid bodies 59
(and unknowns) is reduced to five. Again, if all the forces are parallel to
some line, the number of useful independent eciuations is reduced to three,
whereas if all the forces intersect some line, there are only five independent
equations of equilibrium.
Although we have specified the form of the three scalar equations of
equilibrium [Eqs. (4-10)] for plane problems, it is not necessary that they
the same straight line. We leave it to the reader to verify that these
alternative sets of equilibrium efiuations do indeed satisfy Ecjs. (4-8) in
any plane problem.
In a similar fashion, it is possible to enumerate alternative sets of six
equations of etiuilibrium, satisfying Ecjs. (4-8), for three-dimensional
equilibrium problems. However, to explicitly describe these equations
and their restrictions more involved; the number of different sets is
is
native sets of equations of equilibrium lies in the fact that now we have
the opportunity to select our equations such that as few unknowns as
possible are involved in any one equation. Thereby we materially reduce
the effort required in the solution' of simultaneous algebraic equations.
As an extreme, if we are able to select six equations of equilibrium such
that each equation contained only one unknown, no simultaneous solution
would be necessary.
two forces must be equal in magnitude, but opposite in sense, and their
lines of action must be colinear.
Fig. 4-3
Examnle t-3. Tlie simple plane st riictiire shown in V\^. \-'m supports
a load of 3.000 lb. If the weijrht of the bars is neglected, find (a) the
reactions at the supports .1 and H and {h\ the reactions at the pin joint C.
Assume all joints are frictioidess.
12501b
62501b
Fic. 4-5
Solution, (a) First we take the entire structure as our free body. The
free-body diagram is shown in Fig. 4-56. The unknown reactions at the
supports .4 and B are
A = .l,i + .IJ at .4
B = BA at B
Applying the e(|uations of eriuilibrium.
i:F. =
SM,-' =
62 Statics of systems of particles and of rigid bodies
SF. = C. + HD =
ZFy = Cy + V5D - 3,000 =
2M/ = 9i%D) - 15 X 3,000 =
Solving these equations,
Fig. 4-6
(c) If we finally take bar AB as a free body, we can verify that all the
pulley B is vertical. Assume that the belt tension on the tight side of the
pulley is three times that on the loose side. Find the reactions in the two
supporting bearings, assuming that bearing 1) alone prevents motion in the
X directi(m (Fig. 4-0).
To transmit 50 hp at 800 rpm requires a tortjuc of
we find that
7'i = 1 .050 11) r. = 350 lb
1\ = 750 lb 7'4 = 250 lb
Let us now take the complete shaft-pulley as.sembly as our free body.
The forces acting on this free body are
Fig. 4-7
If the coefficient of static friction between the end B and the wall is /is,
find the angle a at which the bar will just start to slide down (Fig. 4-7).
When the bar just starts to slide down, the frictional force at B tends to
stop the motion and therefore tangent to the circular arc DBE. Taking
is
the bar AB as our free body, and considering the equilibrium of moments
about joint A, we have the following forces that contribute to the moment
about A :
M,r =
T»/i
r„ X
«-
W = W /t:.
- «"
., .
«
.
+
,
a IF
k
-2" V^" •'^m i .^^
S(v. 1-6] Couple 6.')
—F +
W sin
— rt =
— aF cos a + X \/l' — a- sin « =
Solving the first two cfpuitions (the thii'd eciuation ;ul(ls no new informa-
tion in this jiroblem),
^V
rJ? = — sni a
a\V cos a
A^ =
Ustt
tan a =
Vl' - a'
4-6. Couple
A couple is defined as a pair of ecjual and y^^^, 4.^
opposite forces whose lines of action do not
coincide. Let F and — F he the pair of forces passing through points .1
and B (Fig. 4-8). The moment of this couple, f al)out an arbitrary point
0, is
M = -r X F + (r + u) X F
or M = a X F (4-11)
where a is the vector connecting the two arbitrary j^oints on the lines of
action of the two forces B and A. Equation (4-11) shows that the result-
ant of a couple is independent of the location of the moment center. A
couple vector, therefore, has the same moment about all points in space
and thus is a free vector. Furthermore, since
|a X F| = /'' (projection of a ± F)
= /•' (perpendicular distance between F and — F)
= /''(/
we consider a couple as a force system, it has the peculiar property that its
Two force systems are said to be equipollent^ if they have (1) the same
total force and (2) the same total moment about the same moment center.
The significance of the concept of equipollence of force systems lies in the
following statement: Equipollent force systems acting on a rigid body
influence the equilibrium or motion of the body in exactly the same way.
To prove the above statement, let us assume that a system of forces S
consisting of Fi, F2, . . . , F„ acts on a rigid body (Fig. 4-9a). For each
of these forces, we can calculate its moment with respect to a common,
arbitrarily assigned moment center 0. We then form
2F = Fi + F2 -H • • •
+ F„
SM = Ml -1- M2 + • • •
+ M„
The graphical construction of SF and 2M is indicated in Fig. 4-9/>, in
which we define
^'- " ^^ (4-12)'
M, = SM ^
In general, F^ and Mr do not vanish, and hence the rigid body is not in
equilibrium.
t We use the term (MiuipoUont, rather tlian the term eeiuivnleiit or even mechanically
equivalent, because of its by its definition. If
precise technical meaninjr as supplied
we consider a defornuible bar acted upon by two and colinear tensile
equal, opposite,
forces, this force system would be e(juiiwllent to a force system composed of two com-
pressive forces, or even equipollent to a zero force system. However, the three force
systems would not be equivalent; certainly their effect on the deformable bar would be
different in each case.
Sec. 4-7] 'Eqnipolletice of force syslems in space 67
We now
add to the force system S a pair of ecjual and opposite forces
Fr and F„+i and a pair of equal and opposite moments Mr and M,.-i-i.
(To simplify our demonstration, we may let Fr and F„+i act at the point
so that neither force contributes any moment about 0. We also let Mr
and M„+i be two couples so that neither of them contributes any net
force.) Now the state of e(iuilil)rium or motion of the body is obviously
not changed by the addition of such e(iual and opposite pairs of forces and
pairs of moments. But now the sum of forces Fi + Fo + -f F„ • • •
(b)
Fi(i. 4-0
both zero. They actually satisfy the etiuations of equilibrium, and hence
can be removed from our consideration. Therefore the cfiuiiibrium or
motion of the body is due solely to the influence of a single force Fr acting
at and a couple Mr. Since the force system that consists of Fr at and
Mr is, by definition, equipollent to the original force system S, this proves
that ecjuipollent force systems exert the same influence on the equilibrium
or motion of a rigid body. In terms of force systems, it now becomes
apparent that the neces.sary and sufficient condition for the c(iuilibrium
of a rigid body is expressed by the fact that the total external force system
be equipollent to zero; i.e., Fr = Mr = 0.
It can be proved that if two force systems are equipollent with respect
to some moment center 0, they are etiuipoUent with respect to any other
moment center, say C/. The first part of the condition of equipollencc,
namely, that the resultant force is the same for both systems, is obviously
not influenced l)y any shift of the moment center. The second part of
68 Statics of systems of particles and of rigid bodies
2Mo = S(r X F)
is the same for both systems. If we now shift the moment center from
to 0', the new resultant moment will be
SF and M, = SM
This is called the reduction of the
force system to a force and a couple.
In particular, a single force F atP can
always be reduced to a force and
couple acting at some other point 0.
Here the equipollent system would
be
Fk;. 1-10
F, = F M, = r X F
where r is the position vector from to some point on the line of action
of F. As we shall see shortly, the converse is only valid when the force
A force system consisting of a force and a couple, with the same direc-
tion for both, is called a wrench or a screw. The original force .system is
F, = 2F and M; = ^M^.
where Ml represents the vector component of each individual moment in
the direction of Fr.
To conclude: A general force system can he reduced to (l) a force and a
couple in an infinite number of ways, depending on the location of the point
on which the force acts, or (2) a wrench in only one way.
For the special case when a force system consists of parallel forces
only, Mr is always perpendicular to Fr, and hence consists of cross com-
ponents only. From the previous discussion, it is .seen that such a force
system can be reduced to (1) a force and a couple in an infiiiile number of
ways; the force itself must be Fr = 2)F, but the value of the couple
depends on the position of the moment center; or (2) a single force in
one way, where again the single force must be Fr = 2F and its line; of
action is determined uniciuely by the perpendicular distance d fioni an
arbitrary point to the line of action of F^, d = Mr/Fr. One important
applicaticMi of such a force .system is when the parallel forces are the
weights of the particles of a system; then the location of the e(|uipollent
single force becomes the center of gravity of the .system and the single
force becomes the total weight (gravity force) of the .system (body).
This will be discussed in some detail in the next chapter.
()th(!r special forcesystems merit some comment. Concurrent force
systems (Mr =
about the point of concurrency) are clearly (!(|uipollent
to their resultant pas.sing through the point of application and may be
reduced to a force couple system in an infinite uMinbcr of ways and a
70 Statics of systems of particles and of rigid bodies
Fi=100lb
F2 = 80 lb
F, = 60 lb
Fig. 4-11
It should noted that the resultant moment of any coplanar force sys-
})e
tem about any point in the plane consists only of a component perpen-
dicular to the plane. The original force sj^stem is now reduced to an
equipollent system consisting of a single force Fr acting at and a single
couple Mr.
To further reduce this system consisting of a resultant force Fr and a
resultant couple Mr to an equipollent system consisting of a single force,
we note that M,. is jiorix^ulicular to F,. and thus Mrc = Mr. Therefore
Sec. 1-8] Reduclion of force systems; wrench 71
d = ^lr ^ 214. 2
=
Fr
~ 91.7
2:A in.
aiu r X F' = Mk = Mr
Letting r = .)vi + //,j + 2rkandsolv 1.16x,-3.57
yr = I.IG.JV — 3.57
Example 4-7. Given a general system of forces Fi, F2, and F3.
Reduce this system to an equipollent system of forces at the point and
further reduce the .system to a wrench.
Solution. The given forces and their respective position vectors are
Fi = 8i - 4j + 2k r, = Oj
F2 = 3i + 2j + 3k r, - Si + 2k
F3 = i + Oj - Ok rs = i — k
The given system is equipollent to a force couple system consisting of
Fr = 2F = 7i + 4j - i
Mr = 2M = 2r X F = I4i -f- 2k - Ok
Referring now to Fig. 4-10, we find that
_ _ Mr • F, Fr 50 ,_. ,
,
. , .
The equation of the line of action of the ecjuipollent wrench may now be
determined by letting
r X F' = MrC
72 Statics of systems of particles and of rigid bodies
d = Mr 2_ 12,507
33 66
PROBLEMS
4-1. A
"simple" beam is supported at one end by a fixed hinge and at the other
end by (The rollers are such that complete restraint is offered in the
"rollers."
vertical direction, both upward and downward, but no restraint is offered in the hori-
zontal direction.) The beam is loaded as shown. Find the reactions at each end.
20001b
10001b
12'
Fig. P 4-1
4-2. A beam, fixed in a wall so that no rotation or translation can occur, is loaded
as indicated. Determine the reactions necessary to maintain equilibrium.
20001b
12001b
|2000lb
2' I
4'
18'
Fig. P ^-2
A small truss-type highway bridge has a span length of 80 ft. The trusses,
4-3.
of which there are two, are assumed to consist of weightless bars connected together
with frictionless hinges. A loaded vehicle is positioned in panel he as shown. The
weight distril)ution of the vehicle is 2,800 lb on the rear axle and 2,200 lb on the front
axle. Assuming that half of the vehicle weiglit is distributed to each truss and that
the forces acting on tiie panel points h and c are carried through the floor system by
simple beam action, determine the reactions on a single truss. Find the forces in the
members /(; and U.i. Hint: Since llic truss is in (Mpiililiriuni, each portion is in equi-
Problems 73
librium. Hcncp, solcct a portion as a fnM> Ixxly (consitlcnnl also as a r\^,'u\ \mdy) that
will I'liabie a dctormination of the nMiuirt'd fort-cs.
Fig. P 4-3
4-4. A piano truss sliown is liiiipiod at .1 and siiiii)ly supported <ni rollers at D. The
bars CD and AB make G0° with the horizontal, and the bars (\\ and lilJ make 30°
with the horizontal. A total load of 3P is distributed at the two bottom joints as
shown. Find the reaetions in the bars AH, CD, and EF.
40001b
'.jX
1000 lb
4-5. Determine the reactions re(iuircd for efiuililtriuni shownfor the structure
Fig. P 4-7
4-8. A boom AB, 12 ft long, is held at one end A by a frictionless ball joint. It
carries a load of 1,000 lb attached to the middle. Besides the joint A, the boom is
also supported by two wires attached to the boom at B and C as shown. Find the
reactions in the wires and at the ball joint.
100 lb
Fic. P 4-10
i-11. A uiiifonn box lid of wfiight IF with hinges at .1 and li is lifted :}()° from the
horizontal by a rope attached at one lifted corner C. Find the reactions of the
hinges if (a) the rope is horizontal and parallel to the x axis, (b) the rope still lies in the
xij plane but makes 30° with x.
ABC lies in the xy plane and triangle DEF lies
4-12. In the truss shown, triangle
plane parallel to the xa plane. AB = AC = DE = DF = EH = DA =
in a vertical
FC = 1 ft. The truss is supported by a ball joint at A, a, ring guide at D, and a roller
support at C. A vertical load \V is api)lied at the joint F, and a horizontal load F is
applied at the joint E. Assume all joints are frictionless. (a) Find the reacitions at
the three supports (some reactions may be negative), (h) Find the reaction on the
l)ar AF. State whether it is tension or coinprcssion.
37501b
llOlb
4-13. The forces indicated are acting on an automobile. Find an equipollent sys-
tem of forces acting at 0. Further, reduce this system to a single equipollent force.
4-14. Determine the relation between the magnitude of the forces shown, so that
they may be reduced (i.e., are equipollent) to a single force.
76 Statics of systems of particles and of rigid bodies
4-15. A force system consists of two forces Fi and Fo and a couple C acting in the
plane of abed.
Fi = 5\h F. = V2 lb
Oa = 4 ft 06 = 3 ft 6c = 3 ft
4-17. What is the necessary relation betweenn, b, and c so that the three noninter-
secting forces (of equal magnitude P) acting as shown may lie nMluced to a single
force?
4-18. (live ail example of a body or system that satisfies I'^js. (4-S) and yet is not
in equilibrium.
4-19. tShow that Mrc and Fr are equipollent to a single force V, located at a position
r from O as defined by r = F, X Mrc/Fr. Show also that this results in the location
of Fr at a perpendi(Mdar distance d = Mrc/Fr-
CHAPTER 5
\
More Applications of Statics
j
5-1. Systems of Parallel Forces
j
It was discussed in Sec. 4-7 that a sj'stem of forces, in ^fncral, is
It is this last situation that applies when all the forc(^s in the system are
parallel, since then .M, is always perpendicular to Fr. A system of
78 More applications of statics
For this purpose we orientate our coordinate axes so that the y axis is
along the direction of the parallel forces
(Fig. 5-1). The forces in the given system
can therefore be represented by
Fi = Fij
F2 = F2J
F, = SF - (Fi + F2
+ •
)j = (SF)j (5-1:
F, = SF
The location of F^ must be such that
M, = SM (5-2)
for equipollence. But with all the forces having the y components only,
M, = -z,Fri + XrFM
and SM = i:{-zFi -\- xFk)
where x and z represent the coordinates of the lines of action of the forces.
Hence Eq. (5-2) is equivalent to the following two scalar equations:
XxF
I
SF
(5-3)
Zr = S2F
SF
Equations (5-3) give the x and z coordinates of the line of action of the
single resultant force equipollent to the given system of i)arallcl forces.
It is also clear that the total moment S!M of the gi\-ensystem of forces
vanishes about any center that lies on the liiu^ of action of the equipollent
single force. Thus a singl(>, force F^ = SF acting at coordinates x — Xr
Sec. 0-2] Center of gravily: center of mass: ceniroid 79
. and ij = ih on a rigid lH)dy would produce the same ellect as that caused
'
by the original system of parallel forces.
j
body is defined as the point (denoted by C) of action of the single, resul-
i tant force equipollent to the system of gravitational forces. Considering
the gravitational forces as comprising a system of parallel forces is valid
only when the distances between the various particles is sufficiently small
so that all particles lie in an essentially constant gravitational field.
Fortunately, the vast majority of bodies we deal with in engineering
applications satisfy this condition. Then Eqs. (5-3) can be written in
the following form for the calculation of the coordinates {.vcUc^c) of the
center of gravity:
xc
^VxiVi Sr.-w;.-
Sz,-ty»
zc =
where x„ r„ and iVi refer to the x coordinate, position vector, and weight
of the ith particle, respectively.
In the above we have extended E(is. (5-8) to obtain not only .re and Zc,
but also an equation for ijc. This can be done by merely rotating the
body as well as the coordinate axes so that either the z or x axis becomes
vertical, and hence parallel to the gravitational forces. The .same
reasoning that was used in deriving Eqs. (5-3) then gives a similar
equation for ijc.
f X dw I y dw I z dw
xc = . yc = ^ .dw Zc - r
- (.o-oj
\dw \ I dw
t Values of the centroid of some plane geometric shapes and the center of mass of
some homogeneous bodies are given in Appendixes .\ and B.
80 More applications of statics
where the integration is extended over the entire volume of the body.
Equations (5-5) may be written in the alternative form
/ r dw
re
/.-
dw
This concept of the center of gravity depends, for its definition, on the
existence of a gravitational field. Another concept not so dependent
may be defined. This is the idea of center of mass of a body, where we
define the center of mass as that point (denoted by C) having coordinates
/ X dm /
^ y dm j ^ z dm
XC = —r VC = —r 2c = —r
/ dm, / , dm / ^
dm
(5-6)
I
^ r dm
or re
dm
/.
the body. The numerators are often referred to as they^rs^ moment of the
the first moment of that portion of the mass to the other side of the plane.
In particular, then, the center of mass must lie in a plane or planes of
symmetry.
By analogy with Eqs. (5-5) or (5-6) , it can be defined that the centroid
of a volume V is located at the point (denoted by C) with coordinates
._/v xdV
X= ^, (5-7a)
/ x'dA
with simihir expressions for // and z. Finally, the centroid of a line seg-
ment L may be defined to be located at the point with the coordinates
! xdL
^ - 'r (5-7c)
has the shape of a slender wire, not necessarily in one plane), the center of
ma.ss of the body approximately coincides with the centroid of the line
segment which forms the spatial axis of the body [cf. Eqs. (5-6) and (5-7c)].
When this body possessing uniform material has a constant cross-section
area, the center of mass coincides exactly with tlu; centroid of the
spatial axis.
Although the previous discussion refers directly to continuous bodies,
these ideas may easily be (extended to include: cases of discont iimous bodies
of a finite number of j)ortions. For this event, relationships similar to
Eqs. (5-4) are obtained.
82 More applications of statics
- = ^' (5-8)
where Vi and pi represent the position vector and weight of each point,
respectively. If instead of having a discrete number of points in the set
we were to have a continuous distribution, the centroid would take the
form
_ jrdp .
^^-^^
""'-Jdp
where p now is a weighting function. Equations (5-8) and (5-9) can, as
before, be written in terms of the rectangular coordinates of the centroid.
It is noted that this more general definition of centroid is in perfect
agreement with our previous development, if we let the weight or weighting
function be the particular quantity of interest. In a certain sense, then,
the centroid represents the average position of the weighted set of points.
We also note from Eqs. (5-8) and (5-9) that the position of the centroid is
Sripi =
Suppose there is another point D which satisfies Eq. (5-8). Then
Sr^pi =
where r •
is the position vector measured from an origin at D. Subtracting
these two expressions,
Sp,(r. - rO = (Zp^i =
where d is the position vector from (' to D. The Inst st(>p follows since d
isa constant, independent of each individual point. T'rom this it follows,
provided Sp, 9^ 0, that d == 0. Thus the points (' and D coincide, and
there is a unique centroid as defined by Eq. (5-8).
The reader should refer to a book on integral calculus and review the
Sec Center of (fravily: center of muss: cenlroid 83
an axis in its plane is equal to the product of the length of the generating
curve and the distance of travel of its centroid during the rotation, if the
generating curve lies only on one side of the axis.
2. The volume of the solid generated by rotating any plane area about
an axis in its plane is equal to the product of the generating area and the
distance of travel of its centroid during the rotation, if the generating area
lies only on one side of the axis.
The proof of this theorem can be found
in any standard text on calculus. The
following examples show one applica-
tion of this theorem.
Example 5-1. Determine the /y
y = -
Fig. 5-2
The equations of the straight line and the parabola are given, respectively,
as
y 6-" //
= p.T-
dy
y =
-iri -lU dy
y = Hh
Example 5-2. Knowing that the volume V of a sphere of radius r is
V = %Trr^
:
and that the surface area is *S = Airr'^, use the theorem of Pappus to
4 _ Trr^
V - 7rr=* = 2'n-yA -^
4r
Hence Va
3^
S = 47rr2 = 2TryL
2r
and therefore
re (5-10)
J
Here the distributed ma>ss of each portion is considered to be equipollent
to having the total mass of the portion concentrated at a mass point,
namely, its own center of mass. We leave the proof of this statement as
an exercise (Prob. 5-25). The composite body can therefore be considered
as consisting merely of a finite number of mass points, and Eq. (5-10) can
be directly applied. A similar procedure can be used to find the centroid
of composite vokimes, of composite areas, or of composite line segments or
the center of gravity of a composite body.
Exaiiii>le 5-3. A crank-and-shaft assembly has the dimension and
shape as shown in Fig. 5-4. Assinning the assenU:»ly is made of uniform
material, find its center of gravity.
Since tlu; assembly is made of nnifoi'm niatcM'ial, the center of gravity
>t'C. -3] Pressure ht a sialic Ihiid 83
coincides with the centroid of the vohiine. Let us. for convenience, phice
our coorilinate system so that the indiviihial centroids of the five portions
all lie in tlie .r/y phme. Thus zc ^ 0, and the problem therefore is reduced
to a two-dimensional one. We also arbitrarily place the origin at the
left end of the assemblv and at the axis of the main shaft. The volumes
86 More applications of statics
Pv - Ps - ^
y dy
= „
t The proof, Imscd on tho p;onoral consifloration of stross tensors, w ill ho p;ivcn in
Chap. 15.
Sec. 5-3] Pressure in a sialic fluid 87
|P2dA-(p+ dp)dA
PidA p^dA dA
^ dA
dy
ydydA
along the x direction are the pressure forces on the two ends of the tube.
The equation of eciuilibrium 2Fi = gives
pi dA - P-,
dA =
or Pi = P2
p dA — (p + dp) dA = 7 dij dA
or dp = -y dy (5-11)
The negative sign indicates that the pressure decreases with increasing
altitude. Equation (o-ll) is referred to as the pressure-altitude relation
of a static fluid.
actual calculation of the pressure difference between two points
The
separated by a vertical distance in a compressil)le fluid re(|uires the
integration of Eq. (5-1 1). In order to do this, we nuist possess a relation
between 7 and p. For example, if the fluid under consideration is a
perfect gas under constant temperature, we have the well-known iso-
thermal relation
p = constant X 7
88 More applications of statics
This constant can be evaluated if we know the pressure and the specific
weight pi and 71 at a point of elevation yi. Integrating Eq. (5-11) then
gives the pressure p2 at some other point of elevation i/2 as
y2
- ^A
^— = exp (5-12)
Pi 7>i/7i /
As another example, consider a gas that follows a polytropic relation
p = constant X 7"
where the exponent n may have different values depending on the problem.
Integrating Eq. (5-11) gives
n pi /pA("-r)ln-
y2 - iji
— 1 (5-13)
n 1 7i
the manometer li(iuid with specific weight 72. The pressure at point A,
where the two liquids come into contact (the two licjuids must be immisci-
ble), is the same as the pressure at point B, since the two points have the
same altitude and are connected continuously by the same lic^uid 72 in
the tube ACB. Hence
Thus we obtain the j^ressure at I'roin the known ])ressure at 2 and the
1
of water or in inches of mercury (in. Hg). Finally, all the above vuiits
can be either a gauge pressure (e.g., psig) which is measured above the
atmospheric pressure or an absolute pressure (e.g., psia). The specific
weight 7 for water is 02.4 lb ft''. Hence
02 4
1 ft of water = ^^T ~ 0.488 psi
441
12
1 ft of water = ,., -. = 0.884 in. Hg
p\ = ZF = rCp A.l
We now let AA approach zero as a limit while at the same time the number
of forcesapproaches infinity. The summation in the above e(}uation
becomes an integration taken over the entire area of the plate; i.e.,
Fr = p dA
f^
I
xp dA
Similarly, Xr = ^ (•''>1G)
f^pdA
UP dA
Vr = f^
90 More applications of statics
It is recalled that only the x and y coordinates enter into the problem if the
direction of forces is along z.
Let us now apply Eqs. (5-16) to the case of a fiat surface immersed in a
liquid for Avhich the simple pressure-altitude relation Eq. (5-14) holds.
Figure 5-8 shows a flat surface AB submerged in a liquid and inclined at
an angle d with the horizontal. We take the intersection of the plane of
AB with the free surface as our x axis. The y axis is taken to lie in the
plane AB. This means that the plane AB becomes our xy plane.
Fig. 5-8
To compute the resultant pressure force Fr due to the liquid on one side
of AB and its location Xr and yr,
we proceed as follows. Take an element
of areadA formed by two horizontal lines on AB at a distance dy apart.
The pressure on dA due to the liquid is
p = yh = yy sin 6
Pc = yy yin 6 = yhc
Hence we can also write
Fr = pcA (5-18)
Thus the magnitude of the resultant pressure force on one side of aflat surface
is the product of the area of the s)irfa<-(' and the pressure at the centroid.
Sco. .i-1] Forces due lo Jhiid sialic pressure: plane surfaces 91
Since p = yy sin 6, Eq. (o-17) can also ho interpreted as the total volnme
under the pressure diagram acting on the plane area.
The location of the resultant pressure force Fr is described by the
coordinates x,, ijr of the piercing point of Fr with the given area. This
point is called the center of pressure. iMjuations (o-lOj yield
dA
Ja''^ III
Xr
yA
f.ydA
(5-19)
l.y^dA h
yr = -,
yA
l^udA
r 6r T
'''" M!
^"v - 12 12
bhH2a + b\
K 2A
If Ix and Ixy of a given area A are known with respect to a set of axes
X and y, we can compute the moment of inertia
Ix' and the product of
inertia Ix'y of A
with respect to a new; set of axes x' and ?/' which are
parallel to and displaced from the x and y axes (Fig. 5-9). The coordi-
nates of a small element of area dA are x, y, when referred to the original
t''
Sec. 5-1] Forces due lo fluid .sialic pressure: plane surfaces 93
and // =
f^
//'- dA =
j ^ (// + h)- dA
= 1^ if dA + 'Ih dA + Ir
j^ dA
f^,!,
= /. + 2byA + h-A (5-21)
Consider now the special case in which the original axes .c// are the
centroidal axes, meaning that their origin coincides with the centroid of
the given area. For this case x = y = 0, and Ecjs. (5-20) and (5-21) are
reduced to
x'y' = Izryr + dljA (5-22)
h' -^ I.r + h'A
where the subscript C is used to emphasize the fact that the original axes
can be only the centroidal axes if the above etjuations are to hold. Once
the original axes are designated by
xcyc, ^ve can afford to drop the prime
on the new axes x'y' and designate them
simply by xy (Fig. 5-10). Then, if x, y
are the coordinates of the centroid C
with respect to the new axes xy, we see
that X, y are just the (luantities a, h in
1 xu i xcyc I
XyA.
(5-23) Fic. 5-10
h = he + r-'^
Equations (5-23) are known as the parallel-axes theorems for the moment
and the product of inertia of an area. Sub.stituting Eqs. (5-23) into Eqs.
(5-19),
Xr
yA
+ x
(5-24)
yAf + ^
J zryr O and Xr = X
Example 5-4. Find the location of the resultant force on one side
of a rectangular plate submerged ver-
- - / tically in a liquid. The upper edge of
the plate is at the free surface (Fig.
,.2
r
Fx = y h — -^] r, a.t a vertical distance .,„,, y^ + ^ below B
{
\ 2/ 12(/i — r/2) 2
= 7r2 ( /i — r +
The force P required at A to hold the gate in place is
vy
G*I
B*
ib)
other hand, if the cone is rested on its side (along one of its generators),
the cone is said to be in a position of nrulral ((luilihrium.
By .stable, then, we imply that any slight departure from the neutral
position of the body willcause a net force, or moment, that tends to
restore the body to its neutral position. A floating body may or may not
be stable if its center of gravity is above the center of buoyancy. The
following analysis will show the reason why.
LetG be the center of gravity of a floating body of weight \V (Fig. 5-15)
and B its center of buoyancy. The body is in e(iuilibrium if undisturbed.
The buoyant force Fn is equal in magnitude to the weight of the body.
Now if the body is disturbed so that it rotates about its longitudinal axis z
t We .shall discuss this point further in See. Ki-2.
98 More applications of statics
by a small angle dd, the displaced volume will increase by the additional
submergence on one side of the z axis and decrease by the same amount on
the other side. This results in an increase of buoyant force on one side
and a decrease on the other. Hence a restoring couple is formed which
tends to counteract the overturning moment caused by the lateral dis-
placement of the center of gravity. To evaluate this couple, let us
consider a small element of volume dV in the two wedges that represent an
addition of submergence on one side and a reduction on the other. The
element of volume is
dV = (x dd) dA
M, = f 7 {dd)x-' dA = yidd)I
where I is the moment of inertia about the z axis of the entire cross-
section area of thebody at the free surface. This restoring moment M,
combined with the original buoyant force Fb forms an equipollent single
force F'g, which may be called the new buoyant force. The location of
the new buoyant force is determined by the distance BB' (Fig. 5-15),
where
, M. M. y{dd)I I
^^ ^ V ^^
-f;-~w- ~w~
where V is the displaced volume of the floating body. The intersection
of the new buoyant force with the y axis is the point AI, called the meta
center. The distance BM is
sin do
BM =
y (5-25)
Equation (5-25) gives the distance of the meta center from the center
of buoyancy. If BM
is larger than BG, the restoring moment is larger
than the overturning moment and the body is stable. On the other hand,
if BM is smaller than BG, the body is unstable. Finally, if coincides M
with G, the body is in neutral equilibrium. It should be emphasized that
/ in Eq. (5-25) is the moment of inertia of the cross-section area of the
floating body at the free surface. This cross-section area is, in general, not
the deck area.
Problems 99
FKOBLEMS
5-1. Determine the mass center of a ho niogenco us pyramid of base dimensions a
anil '' and height h.
5-2. Determine the eentroid of the area inrhided helween the i)arahola tj =
(h/b")!", the X axis,and the Hne x = 6.
5-3. A small cylinder of radius r and height is made of material of uniform {h'nsity.
/*
However, there is a large empty spherical cavity w hose radius is r/2 located as showu.
Locate the center of mass.
o-i. A beam is loaded by a distributed load, the intensity of load (force per unit
Fi(i. I' ^A
100 More applications of statics
radius of the earth is 4,000 miles. Where is the center of mass of the earth-moon sys-
tem relative to the surface of the earth?
5-6. Orange juice is held in an essentially weightless conical Dixie cup of height h
and base radius r. After some time there is some sedimentation, so that the density
of the orange juice is given by 7 = 1 + ky, where A; is a constant and y is the distance
from the top of the cup. Locate the center of gravity of the juice.
5-7. A uniform circular plate is supported in a horizontal position by three vertical
bars symmetrically located around the circumference of the plate. The plate itself
weighs Wi. A second load W2 is applied at the mid-point of the line joining the cen-
ter of the plate to one of the supports. Find the reactions in the supports.
Fig. P 5-7
5-8. Knowing that the surface area of a sphere of radius r is 47rr^ find the location
of the centroid of a semicircular arc by the theorem of Pappus.
5-9. Find the center of gravity of a slender wire OABCD consisting of three straight
portions OA, CB, and CD, each of length a, and a fourth portion in the shape of a
quarter circle of radius a. Assume the wire is of homogeneous material.
5-10. A somicylindrical shell of uniform wall thickness weighs wi Ih/lin ft along the
direction of the axis of the shell. It carries along its one edge a linearly distrihuted
load of Wi lb/ft. ?'ind the relation between the ratio «'i/h'j and the angle of inclina-
tion a when the shell rests on a horizontal table.
/////////7///////////
Fig. P 5-10
5-11. A semispherical shell of uniform wall thickness weighs H'l lb. It carries at
5-14. A pipeline contains natural gas, 7 = 0.045 lb/ft'. At one point a simple
water manometer shows the pressure to be 6 in. of water. What is the pressure in
inches of water at a iM)int in the line 500 ft higher than the first pf)int? Take 7 for
air as 0.07G5 Ib/ft^ and neglect the compres.sibility of air and gas for the relatively
small difference in elevation in this problem.
102 More applications of statics
5-15. Verify the magnitude and location of the pressure forces acting on the sluice
10'
8'
I
Fig. P 5-16 Fig. P 5-17
5-17. Locate the resultant force on each end of the cylindrical drum of 3-ft diameter
and 5-ft length immersed water as shown. Also find their respective magnitudes.
in
5-18. At one section of Grand Coulee Dam on the Columbia River, the dimensions
are approximately as follows: The area ABD is 13,750 ft 2, and its centroid Ci is 22.3 ft
to the right of B. The vertical distance DB is 430 ft. The area of the cross section
of the concrete dam is 98,000 ft^; its centroid, C2. (a) Compute the position where the
resultant force exerted by the ground on the base of the dam crosses the base, (b) If
the base BE of the dam is 400 ft, find the maximum and minimum stresses (intensity
of pressure) at the base (assume a linear stress distribution along the base).
430'
Problems 103
5-19. A plane gate woiphing; 500 lb/ft of length perpendicular to the pai)er is
hinged at one end O. The width of the gate is 10 ft, and its center of gravity is in
the middle. Find the (l(>pth li of water as a function of d for e(|uilii)riuni of the gate.
'
1 1 1
-
CHAPTER 6
Derivatives of Vectors;
Kinematics of a Particle
Fig. 6-1
The ratio AA/At is the average change in
A per unit of t during the increment At. As At approaches zero, we assume
that AA/At approaches a limit. (A is then said to be differentiable with
respect to t.) We use the notation
dA
-7-
,.AA
= hm — (6-1)
at Ai-»o At
we can form the derivative of dX/dt and obtain the second derivative of A,
which again is a vector function of the scalar variable t.
d (dX\ _ d'-X
dt\~dt) - w^ (^-2)
constant, and hence its derivative is zero, only when both its magnitude
and direction are constant.
At, then A, B, and C change by an amount AA, AB, and AC, respectively.
The new value of C is therefore
C + AC = (A + AA) + (B + AB)
Hence AC = AA + AB
Dividing each term by the scalar cjuantity At and passing to the limitf as
At^O yields
dC ^ d{X + B ^dX
) dB
dt dt ~ dt ^ dt ^ '
t It Ls asHumc'd here, and in what follows, that the vector fuiictioiis are of such a
form as to ensure that the hniit of a .sum is equal to the sum of the limits and that the
limit of a product i.s ecjual to the product of the limits.
:
dA _
~ dAx d\ dAy d^ dA^ dk
It 'dt^'^^^'dt^iir^
• \ A .
^ ^''dt^~dr^^ ^'H
_, a >
%^ _, *
The proof of the above depends on the distributive law for the dot and the
cross products of vectors. In particular, if B = A,
A-^
dt
=0
which implies either that dA/dt = or that A is perpendicular to dA/dt.
Soc. 6-3] Curves in space: princl/ml iiorniuls <in(i hinorntals 107
dt
ap sin pt i + hp cos p^ j
d'-r
-Tj-o = —ap- cos pt i — bp- sin pt j
dt
p-r
dt'-
Fic. ()-2
dr _ dx . dij . , dz .
(6-9)
ds ds ds ds
108 Derivatives of vectors; kinematics of a particle
We see that dx/ds, dy/ds, and dz/ds are the direction cosines of the tangent
to the curve at the point P; hence dr/ds is along the tangent to the curve in
the direction of increasing s (as was also pointed out in Sec. 6-1, where a
general scalar variable t was used instead of the particular variable s).
Furthermore, since
the magnitude of dr/ds, given by Eq. (6-9), is unity. Thus the derivative
of a position vector to a space curve, with respect to the distance s along the
curve, is a unit vector tangent to the space curve in the direction of increasing s.
det ^
ds
P = -en
ds p
(6-ll)t
where e„ is a unit vector along the direction of det/ds, called the direction
of the principal normal.
The direction of the principal normal can be visualized as follows : The
vector det/ds is zero only when e< is a constant, in which case the curve C
is a straight line. In general for nonzero det/ds, we can pick three points
on C and call themP, Pi, P2, where Pi, P2 are two points
at each side and in
the immediate neighborhood of P. These three points determine a plane.
When we let Pi and P2 approach P as a limit, the plane becomes the
osculating plane of C at P. The unit tangent e^ and its derivative det/ds
both lie in the osculating plane. The direction of the principal normal is
thus fixed as follows It is perpendicular to the tangent of the curve and
:
lies in the osculating plane. It should also be noted that the unit principal
normal e„ points to the center of curvature (Fig. 6-2).
t Equation (6-11) is the first of three equations called Frenet's formulas, which are
of great importance in differential geometry. The other two of Frenet's formulas will
be developed in Example 6-2.
Sec. 6-3] Curves in space: principal nornuils and hinorfnals 109
defines a unit vector perpendicular to both e^ and e„. The unit vector Cb
is called the unit binormal. The three mutually perpendicular vectors e^,
Cn. and Cfc form a moving trihedral, "moving" because all three continu-
ouslj^ change their directions along the curve.
Example 6-2. Since e5 is a unit vector, dcb/ds is a vector that is
perpendicular to Cb and thus lies in the plane of e^ and e„, the osculating
plane. Furthermore, because e^ • ej = 0,
dcb dci
-y-
ds
•
Ct + .
e6 • -y-
ds
=
dcb ,
1 dcb f.
ds p ds
^=
ds
re„ (6-13)
where t is the magnitude of dcb/ds and is called the torsion of the space
curve C.
An examination of Fig. 0-2 reveals that
e„ = Cft X et
de„
= ^ dni . dcb
X
-5— Cb X -. \-
-J- Ct
ds ds ds
de„
ds
-(^^e. + re.j (0-14)
dr , dt
The curvature 1/p, which represents the magnitude of det/ds, can now be
found to be
1 _ yx — xy
. dy dy dx dy
X = I y
.
1 d^y/dx'~
p [1 -1- idy/dxyf'
V = (6-15)
^;
the path of the particle. The acceleration oi the particle is defined as the
rate of change of velocity with time; hence
dy d-r
X2 — Xi = Vxdt
I
Vx, — Vx, = /
"
ttx dt
y = 2xx (a)
But x^ -{- y^ = v^
Hence
1 -\- 4x2
Differentiating Eq. (a) and substituting for x and x by the results just
obtained,
y =
= 'Zxx -{- 'Zx- = -j—
(i
At point P,
Hence
Sec. 6-6] kinematics of a parlicle: tangential, normal components 113
The magnitude of a
is 2v-/{b \/5) and the direction of a makes an angle
of tan~i _i..^ =
153°26' from the positive x direction. As an exercise,
the student should verify that the direction of v is tangent to the path
of the particle whereas the direction of a is not tangent.
rfr
_
ds
where cj is the unit vector along the tangent to the curve C. Hence the
velocity of the particle is
dr dr ds . ,,. ,„.
'' = 'dt
= lsJt = '^^ ^^-^^)
or the velociiij is a vector along the tatigeut with the magnitude of .<-.
The acceleration of the particle is obtained from Eq. (G-18) by the use
of Eqs. (6-16) and (6-4). Thus
a
114 Derivatives of vectors; kinematics of a particle
s — V
s =
For the given curve y = x^, dy/dx = 2x, and d^y/dx- = 2. The curva-
ture is, from Example 6-3,
1 _ d'^y/dx'^ 2
p
~ [1 + {dy/dxYf^
~ (1 + 4a;2)?
2v^
Hence ^"
(1 + 2vy2
4x2)^
2v'
5 V5
The direction of a is therefore along the principal normal. Since the
slope of the curve at P or tan 63°26', the direction of the principal
is 2,
use K(is. (6-10) and (()-ll) to ol)tain the derivatives of Cr and e^.f the
following direct method is perhaps more ilkistrative. From Fig. (i-o,
which shows the e^e^ plane, it can be s-een that l^ecause of a dilTerence of
A0. the unit vectors e^ and e.^ change by an amount
dcr de^
— Cr (0-21)
d^ ~d^
t Considering the unil circle of Fig. G-5 as the curve of Fig. G-2, we see that here
e^ e„ = — Cr = — r ds = d<t> p = 1
Hence -7-^
d<t>
=
dr
'
-r
ds
= c<t> hy Eq. (0-10)]
J^*
= ^= -c. [by Kq. (O-Il)]
Radial: f — r4>^
Tangential: r^ + 2r4> (6-26)
Axial z
r — constant
r(j)e4
Fig. 6-6 ^27r
where <j), the angular speed of the particle, is given by
V'r
a = —rcp-Cr = —
r2 + (p/27r)2
ured from the z axis. Positive directions of these coordinates are shown
in Fig. 6-7. It is easy to verify that any point in space is explicitly
described by this coordinate system. We shall use the unit vectors
Cr, e^, and e^ to denote the three orthogonal directions of the coordinates.
Sec. 6-8] hinenialics of a particle: spherical coordinales 117
Fig. 6-7
Fig. 0-8
.
Since the directions of the unit vectors are not fixed in space, it is
Ae^ = Ad ee
Ae0 = —Ader
Next, consider the change in unit
vectors due to an incremental change in
longitude A<^. As indicated in Fig. 6-8,
all three unit vectors suffer a change of
direction; Cr moves along a horizontal
circle of radius sin 6, ee moves along a
r = rcr (6-29)
PROBLEMS
6-1. Dorivo the compononts of acceleration in the cj'liiulrical coonnnate system
without using vectors. I'roceed as follows:
(a) PVoni the transformation equations, x = r cos 0, ?/ = sin find x and ij in /•
</>,
terms of f, r, and
ij), <i>.
(b) Find the components of x and y along Cr and c^. This gives Or and a,p, which
should check expressions (6-26).
6-2. The transformation equations from spherical to rectangular coordinates are
given as X = r sin d cos y — r sin d sin z = r cos d.
(j>, Proceeding as in Prob. 6-1,
<f>,
find the velocity and acceleration components in the spherical coordinate system.
6-3. Verify Eqs. (6-28) by first expressing the three unit vectors of the spherical
coordinate system in terms of three cartesian unit vectors and then determining their
rate of change.
6-4. If r = a cosh pt \ — b sinh pt j and a, b, and p are .scalar constants, find the
first and second derivatives of r with respect to the scalar varial)le I. What differen-
tial equation is satified by r?
6-5. Find a vector differential equation that is satisfied by
(fps^) 200
20 25 10
35 45 t (sec) 6 t (sec)
40
(a) (6)
250
Fig. P 6-9
6-10. A particle is traveling around the parabolic curve of Example 6-4 with a
constant x component of speed, I'x. Determine the velocity and acceleration of the
particle. Use both the rectangular and tangent normal coordinate systems.
6-11. A rocket-propelled sled travels along a straight track. It starts from rest at
t = 0, and at any subsequent time its speed is given as z; = at^ -\- bt. At time to,
"burnout" (propulsion ceases) occurs. Subsequently, brakes are applied to stop the
sled, the brakes producing a constant deceleration of a = — c. How far has the
rocket traveled during propulsion, what is its speed at burnout, and how much
farther must the rocket travel in order to come to rest?
6-12. An automobile is traveling along a highway whose profile may be approxi-
mated by the If the vehicle is traveling with constant speed, what
sine curve shown.
is its and acceleration? Suppose the vehicle has a constant x component of
velocity
speed; what is its velocity and acceleration? Assume that the vehicle remains in
contact with the surface.
2 I
2 I
6-13. A rod AB of length / slides down in the xy plane while keeping its two ends
A and B on the x and y axes, respectively. If the lower end A travels to the right
with a constant velocity va, find the acceleration of the upper end B and of the mid-
point C. What is the path of C?
6-14. A particle moves with constant velocity along a straight line parallel to the
X axis. Find (a) r and ^ in cylindrical coordinates; (6) the components of acceleration
in cylindrical coordinates.
Problems
Fi<;. P G-16
Fig. P 6-18
Fig. P 6-19
(c) Form the relation between the displacement of D and the crank angle B. (This
relation need not be
(d)
an
Find the velocity of
explicit one.)
Z) as a function of 6. For 6 = 30° and w = 120 rpm, what
I
is the velocity of D?
(e) Indicate the means of finding the acceleration of D.
CHAPIER 7
other points travel circular paths with centers on the axis and in planes
perpendicular to the axis of rotation. Finally, if there is one point in
the body, say point 0, that remains fixed in position, any displacement
of the body
is said to be a rotation about the point 0.
on the bisecting plane AOC, the angles AOC and A'OC are equal. Simi-
larly, since C also lies on the bisecting plane bOC, the angles BOC and
B'OC are equal. Thus the relative positions of the four points ABOC are
identical with the relative positions of the four points A' B'OC.
Hence,
with OC fixed in space, a rotation aliout OC by an appropriate angle will
displace .1 to A' and B to B'. The appropriate angle is, of course, the
angle AOA'. Therefore the original displacement can be effected by a
rotation about OC alone. This proves Euler's theorem.
See. 7-1] Displacemenis of a rhjid body 12'
Ar = An X r (7-1)
Ar = As + An X r (7-2)
We are now in a position to proceed with our original aim, the discus-
sion of the kinematics of a rigid body and, in particular, the basic theo-
rems to which wc previously referred.
Sec. 7-2] Kinenmlics of a rigid body 127
and. during some infinitesimal time interval At, moved to position /*'
Ar
= hm An
dr .. ,.
V = ^77 = lim —- —- X r
ai A/-.0 ^t :u-^o At
or V = 0) X r (7-3)
fore a common base point for the two rotations. Let P be another point
on the rigid body with a position vector r from (Fig. 7-3). The
128 Kinematics of a rigid body; relative motions
Vi = o>i X r V2 = a>2 X r
Vr = o)r X r
commutative; i.e.,
This result of course could have been anticipated since we have already
demonstrated that infinitesimal rotations are vector quantities.
Theorem The motion of a rigid body, due to a rotation of o about a
II.
base 'point and a translation, is identical with the motion due to a rotation
of the same to about any other base point 0' and a new, different translation.
The most general infinitesimal displacement of a rigid body occurring
in atime interval A^ is given by Eq. (7-2). Then the total velocity of a
point P, located by the position vector r from an arbitrarily chosen base
point 0, can be obtained by dividing Eq. (7-2) by and subsequently M
passing to the limit as A^ — > 0. Thus
dr d^ ,
It
= * + "'<
or vp = vo + to X r (7-6)
where vo is the total velocity of the base point 0, and hence is the velocity
of the translation part of the motion common to all points in the rigid
body (Pig. 7-4).
now suppose that the total velocity of P can also be considered as
Let us
due to a rotation of to' about a different base point 0' and a translation.
According to this consideration, the velocity of P is
where \o- is the total velocity of 0', and hence is also the velocity of the
translation part of the motion, the position vector oi P nieasureil
and r' is
vo' = vo + w X a (7-8)
v/' = Vo + w X (r' + a)
(7-9)
Vo' — Vo = to X a (7-10)
yoL = — ,—
CO"
<»>
and whose direction is parallel to gj. Thus the point 0' is oriented with
respect to so that a satisfies the equation
Vo • W w
CO''
^— o) = Vo +I
to X a
so that the location of a point on the screw axis is j^iven hy the expression
vo = — X a. Theorem l\ can
u) l)e stated for this special case:
Theorem Ha. The plane motion of a rigid body at ang instant ran alumgs
be eonsidered a.s o pure rotation about a partieular axis. (The intersection
of this axis with the plane of motion is often called the instantaneous center
of rotation, or the velocity pole.)
The instantaneous center, then, is a point of zero velocity. Its position
in space varies with time and depends on the time variation of u> and v.
Moreover, the instantaneous center need not lie within the physical
confines of the rigid body, hut merely associated with it. The locus of the
instantaneous center with respect to a fixed set of axes (space axes) is
usually referred to as the spcLce centrode, while the locus of the instan-
taneous center with respect to a set of axes fixed in the hody and mo\in}»;
with it (body axes) is called the body centrode. With respect to the
instantaneous center the velocity of any point is given as v/. = o> X b.
where b is the position vector of the point /^ relative to the instantaneous
center.
Many future misunderstandings can beavoidid il the reader takes note
that:
1. A base point "belongs" to the rigid body and tluMcforc nioxcs with
it. The velocity of the base point is the velocity of the translation |)art of
the motion.
2. Although a base point is a point of the body, it need not be in the
body. If a base point is not in the body, we can imagine that it lies in an
imaginary rigid extension of the rigid body.
3. The ba.se point can be any arbitrarily cho.sen point of the l)o<ly.
three points do not lie on the same line. This means that we have to
specify nine variables, three for each of the three points. However, since
the three mutual distances between these three points are fixed for a rigid
body, we have three equations that relate the nine variables. We there-
fore are left with only six independent variables to completely specify the
position of a rigid body.
In.stead of specifying the coordinates of three points, the position of a
rigidbody can more conveniently be described by the following method:
We specify the position of an arbitrary point
first of the body by its
coordinates Xo, yo, and Zq. The body is now still free to rotate about 0.
132 Kinematics of a rigid body; relative motions
Let OX, OY, OZ be a set of absolute coordinate axes in space (space axes)
and Ox, Oy, Oz a set of coordinate axes that are fixed in the rigid body and
therefore move with it (body axes) (Fig. 7-5). Let the two sets of axes be
cut by a sphere of unit radius and center 0. The first Eulerian angle 6 is
Fig. 7-5
the meridian plane ZOz and the plane and measured in the XY
ZOX is
plane. The third angle ^p is the angle between the meridian plane ZOz
and the plane zO.r and is measured in the .vi/ plane. It can be seen that
the first two angles 6 and determine the orientation of the moving axes
(j)
Oz; the third angle determines the amount of rotation about Oz. The
three Eulerian angles,! therefore, completely specify the orientation of
the rigid l)()dy.
t It is unrortuiialc (lial. tlic second and tliiid I'lnlcriaii antj;l('s arc defined din'crcntiy
by did'erenl anthors. Some British autliors use xp second anf^le and
for the for the <t>
third. Still others use tiiid \p that arc dilfcrent by 90° from the ones used here.
Sec. 7-1] Motion of a point of n riijid txtdy \\V.\
where i, j. k are unit Nectors aloufj; the connliiiate axes .r, //, z fixed on the
body. We now seek to find to or its coniponents when the motion of the
body is specified by its Eulerian angles 6, 0, and \^ as functions of time.
Let us extend the great circle /»/*, lying in the xij plane, to intersect the
meridian st at B. .Mso draw the radius ()V in the plane of the great circle
but perpendicular to the meridian. It is clear that the angular-velocity
component in the xij plane is <^ sin d, which can further he resolved into
two components — <^ sin B cos ^ and sin 6 sin ^ along tlie x and y axes,
respect i\ely. The angular-xelocity vector has a component in the ij;
We have seen that the two terms on the right, \o and m X r. represent the
velocities due to translation and due to rotation about (K respectively.
The acceleration of/'' is obtained by dilTercntiating Kq. (7-14) with
respect to time;
a/' = v/' = Vo + w X r -}- o> X r
where o> is the time rate of change <if t.> and is calleil the niKjiitar arrrfrrn-
134 Kinematics of a rigid body; relative motions
(c)
Fig. 7-
vo = w/?i
vo = cjRi
The position vector of any point /* on the wheel from the l)ase point is
The angular velocity and the angular acceleration of the wheel at any
instant are assumed to be
= — wk = — cjk
136 Kinematics of a rigid body; relative motions
vp = vo -H " X r -t- to X (o> X r) = (oiR -\- cbr sin (/> — wV cos (/>)i
Accordingly, then,
Notice that this concept cannot necessarily be extended for the deter-
mination of the acceleration since, in general, the instantaneous center
is a moving point in space.
Fig. 7-8
w^B and the angular acceleration w.ib of the driving crank AB, we can
find the motion of the other members of the linkage by our knowledge
of the motion of a rigid body (Sec. 7-4). In the following we let the
plane of motion be the xy plane. The linear velocities and accelerations
then have only i and j components. The angular velocities and accelera-
tions are along k. Considering the link BC of Fig. 7-8 as a rigid body,
the velocity of C is, by Eq. (7-1 -i),
For plane motions |a>/<c X r/,c| = omcBC, \i^ci) X rnc\ = c^cdDC, and
w.i/i X r,t/}| = u.niAB. Equating the two expressions for vc and sub-
138 Kinematics of a rigid body; relative motions
unknowns wbc and ojcd, either one of which will in turn give us vc.
The above discussion can be illustrated graphically by Fig. 7-11.
This method consists essentially of the graphical solution of the vector
equations just developed. At joint C let us draw v^, which we can easily
achieve once (o^b is specified. At the terminal of \b we draw a line
ubc >^ EC
Fig. 7-11
Vb Vc
BE CE
After determining the velocities of a four-bar linkage, its accelerations
are to be found in the following way. Considering the link BC (Fig. 7-8)
as a rigid body, the acceleration of C is, from Eq. (7-15),
tic/} X DC
Fic. 7-12
We are thus able to solve for the two unknowns umc and wcd, either one
of which will in turn give us ac.
The above consideration on accelerations can be illustrated graphically
by Fig. 7-12. At C we draw the vector an, which can be easily obtained
once (JiAB and w^ic are known. At the head of an we add the vector
<>>Bc X (<OBc X Tnc), which is in the direction from C to B (Sec. 7-4).
At the head of the la.st vector we draw a line perpendicular to CB. This
is line I, representing C^bc X rnc, which is yet of unknown magnitude.
Equations (7-16) therefore represent the total rate of change of the unit
vectors fixed in a vehicle due to the motion of the vehicle.
Let us next consider any vector quantity A. We can represent A, at
any instant, by its components along the moving coordinate axes x, y,
and z and write
A = .4,i + .4 J + .Lk
To an observer in the absolute, or Newtonian, frame of reference, the
rate of change of A is, from Eq. (6-4),
dA
-y- = A = i.i + A,j + i,k 4- A^ -f- .4j + .4,k (7-17)
dt
1
The first three terms of the ri<j;ht-hand side of the above ecjuation repre-
sent the rate ofchange of A relative to the vehicle (i.e., as if the vehicle
were not moving, or as detected by an observer at 0' traveling with the
vehicle). Let us use the svmbol
-^\^
\dt).,,
= AA + .l,j + AM (7-18)
f/A /f/A\ , .
/-, ,^x
It = \ii),., + "' ^ '^ <'-'»>
It
= (jt)r.^ + "^ "^
r = R-1- e
dR dg
^ dt~ ~dt ^t
Fig. 7-14
Since 9 is the position vector of
point P relative to a moving system of coordinates, we have, applying
Eq. (7-19),
%,. = *i +M+ *
is, by Eq. (7-18), the relative velocity of P
with respect to the vehicle.
We also see from Eq. (7-14) that the two terms R w/ X p in Eq. (7-20)
+
are the vehicle velocity at P. Equation (7-20) therefore tells us the
following: The absolute velocity is the sum of the relative velocity and the
vehicle velocity.
The absolute velocity v represents the velocity of point P that an
observer at the origin of the fixed XYZ coordinate system would
detect. The relative velocity pr is that velocity seen by an observer
at 0' traveling with the moving axes and is sometimes referred to as the
Sec. 7-7 Motion referred to a nioring coordumle svsleni 1 13
and p is the position vector of the point P relative to the base point and
can be expressed as
+ g^
(o>,) X 9 + .o, X j^ (5)
Again, applying Eq. (7-19), since we are forming derivatives with respect
to quantities referred to a moving coordinate system, this becomes
R =
146 Kinematics of a rigid body; relative motions
The total acceleration, being the sum of those terms, now becomes
4) = 0k ^ = (/)k
e = {a -\- L sin — L cos (/>)! (f) j
Qr = (t)L cos -f <^L
(/»sini
j <^
The various terms that make up the acceleration of the ball P are
obtained from E(i. (7-21) as follows:
: :
Problems 147
Vehicle acceleration
K =
w/ X p = — a)(a + L sin <^)k
to/ X (lo/ X o) = — a>-(a + /. sill </))i
(The same pr can also l)e obtained hy {liffeifntiating p^, assnniinji; con-
stant and i j.)
Coriolis acceleration
The reader is advised to check each term of the above expression and
to familiarize himself with the significance of the various terms.
Although in this example we have chosen the origin 0' of the moving
axes to coincide with the origin of the absolute axes, this need not be so.
The location of 0' is entirely arbitrary. A change of the location of 0'
merely redistributes the individual values of the three terms that make
up the vehicle acceleration. But the vehicle acceleration, which is
the sum of the three terms, is not affected. Nor do the relative accelera-
tion and Coriolis acceleration undergo any change by a shift of ()'.
PKOBLKMS
7-1. A unit cube undergoes an infinitesimal rotation of magnitude Ad about its
body diagonal OA. Find the displacement of the point A. Verify this result by
considering the rotation to be the sum of three infinitesimal rotations taken about the
three axes at O and also by considering the rotation to l)e the sum of three infinitesimal
rotations taken about the three face diagonals meeting at O.
7-2. Show that in a rigid body in motion, there is in general no point with zero
velocity but that there is, in general, one point with zero acceleration. (Under
special conditions there can be a ix)int or jjoints with zero velocit}' and no jjoint with
zero acceleration.)
7-3. Two points on a rigid body, .1 (1 in., 2 in., in.) and B (2 in., —',i in., in.),
have, at the instant under con.sideration, the following velocities: v^ = Gi + 4j in. /sec
and vb = ( — 12i 8j)c in. /sec.+ Locate the instantaneous center of rotation.
j
Devise a graphical procedure for locating the instantaneous center under similar (but
more general) circumstances. The constant c is some value to be determined.
7-4. Two points on a rigid body, .4(1,2,0) and 5(3,6,0), have, at the instant under
consideration, the following velocities: v^ has magnitude 4 and direction perpendicular
to hne AB, while vb has magnitude 10 and direction
perpendicular to line AB. All
coordinate distances are expressed in inches, and speed in inches per second. Locate
the instantaneous center of rotation. For this type of situation devise a graphicaj
procedure for locating the instantaneous center.
7-5. A wheel of radius a mounted on a shaft of radius b rolls on horizontal tracks
without slipping. If at a certain instant the wheel has an angular velocity of w and
an angular acceleration of w, find the velocity and acceleration of the following points
on the wheel: point A, or the point of contact; point B, or the highest point; point C
on the periphery of the wheel with OC inclined at angle with horizontal. (j)
y/^/////Ay/////////
Fig. P 7-5
7-6. Two points A and 5 of a rigid body are located 4 in. apart along the x axis.
At the instant under consideration.
VA = 3i — j + 2k in. /sec
Vb = 3i + 3j — k in. /sec
VA = i + 3k in. /sec 2
ti = i +2 radians /sec 2
7-7. Two plane wheels, rotating about perpoiidicidar axes through their centers,
are connected by a link C. Wheel A is rotating about a fixed center, while wheel B
Problems 149
is rolling without slipping along the periphery of wheel .1. Find a relationship
between the angular speeds of the two wheels and the link.
7-8. A wheel in plane motion is rolling without slipping inside a fixed eircular track
ol radius R. The radius of the wheel is r, and its angular speed and angular acceler-
ation are given as w and a. Determine the velocitj- and acceleration of tlic center of
the wheel.
Fig. P 7
7-9. A rocker is rolling without slipping along a plane surface. .\t the instant
shown, determine the velocity and acceleration of its center of mass. Tlie radius of
the rocker surface is R, and the center of mass of the rocker is located on the axis of
symmetry a distance a from the rocker surface.
•7-7-7y-r777-rz777'>y, y77-r7-r7-rT7-7-7T7777-
Fiti. P 7-i»
7-10. A wheel rolls without slipping on a track. Given the velocity of .1, find the
velocity of R by graphical construction.
TTTTT
Fio. P 7-10
150 Kinematics of a rigid body; relative motions
7-11. For the mechanism shown with a given velocity of A, find the velocity of P
by graphical construction.
(c) Locate the point, on the rigid body BC or its extensions, at which the acceler-
ation is zero.
Fig. P 7-12
liase point. Find the time rate of change of r if the body has an angular velocity w(,.
7-16. .\ wheel without slipping along a horizontal track with an angular
is rolling
-peed w. We define a fixed reference frame
i, j, and k and a moving reference frame
fixed in the wheel ei, e-i, and Cj. At the instant under consideration e^ makes an angle
e with j. The weight force of the wheel is W. Find the rate of change of relative W
to the moving reference system (a) using the results of Prob. 7-15 and (b) by express-
ing W in terms of the moving reference .sjstem.
Fk;. P 7-IG
7-17. .\n automobile is traveling along a meridian of the earth at a speed of GO mph
relative to the earth. What is the acceleration of the automobile at a location of
30''X latitude? If the automobile were traveling along a circle of constant latitude,
what would be its acceleration? Assume that the earth is a perfect sphere and rotates
with constant angular speed and that its mass center is fixed in space.
7-18. A bead slides down a circular hoop with a constant speed of 2 fps wliile tiie
hoop rotates about its vertical diameter with a constant angular speed of 1 radian/
sec. The radius of the hoop is 1 ft. Find the acceleration of the bead when the bead
is at the position shown.
Fig. P 7-18
152 Kinematics of a rigid body; relative motions
7-19. A helicopter blade oscillates in such a manner that the angle <j) varies with
time t according to
<j> = (t>o
-\-
<l>' sin pt
where <Ao is the neutral position of the blade and <j)' is the amplitude of oscillation.
The blade rotates about the vertical axis with the constant angular speed w. Find
the velocity and the acceleration of a point at the tip of the blade for any value of (j). I
7-20. A centrifugal pump runner rotates at 1,000 rpm. Just before leaving the
impeller, a particle of water P travels at the speed of 100 fps with the speed decreasing
at the rate of 5,000 ft/sec^, all relative to the impeller. Find the acceleration of P.
Fig. P 7-20
:
CHAPTER 8
Dynamics of a Particle
For a particle with constant mass m, Newton's second law gives the
relation between the motion of the particle and the force acting on it:
F = wa (8-1)
We can also write the equations of motion in terms of the normal and the
tangential components (Sec. 6-G)
F, = + =
m(r4> 2r<j>)
^ | (r^) (8-4)
F2 — m'z
Although the above equations of motion are, after Newton's second law,
specified to apply only to a particle, be proved subseciuently that
it will
they are sufficient to solve the motion of a body if the rotational elTect of
the body can be neglected. In many prol^lems, the body under considera-
tion is so small compared with its range of motion that the motion can be
con.sidered as a pure translation. The dynamics of the body then follows
15:5
154 Dynamics of a particle
that of a single particle having the mass of the body and located at the
center of mass of the body. Even if the rotational effect is considerable
and must be taken into account, the above equations, as will be shown
later, still define the motion of the mass center of the body.
In applying the above equations to the motion of a particle, if must be
F depends, in general, not only on time t but also on
realized that the force
position r and velocity v. Thus
F = F(r,v,0 (8-5)
Take, for example, the special case of a one-dimensional motion along the
X axis. We have, in general,
mx = F{x,x,t) (8-6)
X = Xo when t =
X = xo when t =
where Xo = initial displacement
±0 = speed
initial
In general, any two independent conditions are sufficient to enable one to
evaluate the undetermined constants of integration.
Let us consider some examples of one-dimensional motion. The
simplest case is when F, and consequently the acceleration x, is a constant.
After two successive integrations and corresponding evaluations of the
constants of integration according to the initial conditions, Eq. (8-6) then
yields
X = Xo + xot -^ I-
2 m
t'~
(8-7)
The reader easily recognizes that the above equation applies, for example,
to the displacement of a falling body under no other force except that due
to gravity (so-called "free" falling body) provided the gravity is assumed
to be a constant and air resistance is neglected.
The next simplest case is when the force is a function of time only, or
F = F{t). Proceeding again in a similar fashion, Eq. (8-6) gives
X = V andJ
..
n; = —
dv
= _ _=
dv dx
y
dv
dt dx dl dx
F{x) and or the type of solution desired. The examples and problems will
[
serve to illustrate this point.
I From the first substitution we obtain
j^ m dv
and, since v = dx/dt,
, _ mv dv
"^"^ ~ IW
Integrating these expressions and evaluating the constants of integration,
t - f{v) X = g{v)
mv dv
^•' = F{v)
X ^ g{v) or v = g'{x)
dx
Then dt =
(In the above solution the supcrscii])! piinic denotes anotiici- function.)
When the force is a function of displacement only, the ecjuation of
motion takes the form
mx = F(x) (8-6b)
: :
F(x) dx
V dv =
in
of .r or x. For example, the spring in the system may be such that its
deflection bej^ond a certain range is not proportional to the applied force;
or the resistance may be proportional to the square of velocity instead of
to velocity. For such cases the ditYerential equation becomes nonlinear
and no general analytic solution is available.
What has been said above about one-dimensional motion holds equally
well for throe-dimensional motion. The thro(^dimonsional motion, how-
ever, has this added complication: the force component along anij one
direction may now depend on the position and the velocity components in
all three directions and, of course, time. This means, in general, that
mx = F^{x,ij,z,x,y,z,t)
my = Fyix,y,z,x,y,z,t)
mz = F,{x,y,z,x,y,z,t)
t This led Laplace to make the following rcrnark.s (1820): "An intclliKPnt being . . .
could expre.ss in the same e(juation the motion.s of the largest bodies of the universe
and the motion of the smallest atoms. Nothing would be uncertain for him, and he
would see the future as well as the past in one glance." This concept is disclaimed
in quantum mechanics, in which only the probability of motions is considered.
sm = sm
coi CO
(-1)
[ t -\
F repeats itself, i.e., goes through a complete cycle, after a time interval of
27r/co. means that during a unit time, F goes through
It also co/27r cycles.
T = —
CO
= period, sec/cycle
27r
CO = 27r/ = "m" = circular frequency, radians/sec
1^
Fig. 8-2
mx = Fq sin wt (8-11)
speed i; = io at i = 0, we obtain
cohn
(co^ — sin wt) (8-13)
Constant-velocity motion
xo + wm
X = Xu + I i(j H
X = ^
m sin cot
X = — —
Fo
com
cos COf,
X = 5— Sm Oil = :,
OJ^W co-
X = cos oit
aj.f„, (8-14)
X = Xm sin oil
f
It can be seen from the first and third equations of (8-14) that
* + co^x = (8-15)
equation becomes
mx -\- kx = (8-15a)
V dv = — w'^x dx
Integrating, we find that
V = -y/vo^ — (xi^x^
dx
dt
t = - sin~^ —X -\- B
CO Vo
Xo ,
X = Xo cos
.
cof
,
i
.
sin wt
Sec. 8-3] Motion in a resist i tig ntediiint Id I
that
= Vo
—X ' .
X sill cot
03
This is in theform of Eqs. (8-14), noting that here the amplitude of the
f
resulting motion is Xm = yo/w.
P^ifi. 8-4
where </> represents the angle between the velocity vector and the hori-
zontal. Integrating the above e(iuations, using the substitutions
X = dx/dt and y = dy/dt, we obtain
log X = t -^ A
m
log B
h + ^) m
t -f
If xo and ijn are the initial speeds at / = 0, we have, after evaluating the
constants of integration in the above equations,
X = ioe-^*/""' (8-16)
= mg
y + (- + tO g-{klm)t (8-17)
:
Xom
(1 -e — (klm)t (8-18)
mot ni mg -
+ ^
,
?/o )
(1 e-(^^/-)0 (8-19)
F
F = -r
r
See. 8-1] Central forces: planetary and satellite niolions 163
r = —
at
x= —
mr
r (8-20)
d , s dv ,
d\
^^(r X v) = (8-21)
r X V = h = constant (8-22)
the latter two vectors are always constrained in one fixed plane. Thus
central-force motion is always planar, or two-dimensional. Second, we
examine the physical meaning of the magnitude of the constant vector
h. For this purpose let us consider a small time interval dt during which
the particle travels a distance dr = v dt (Fig. 8-5). The area dA "swept
out" by the position vector r during time dt (represented by the cross-
hatched area in Fig. 8-5) is numerically equal to
where the angle between the position and the velocity vectors.
\p is
The
integrated area A over a finite time interval is known as the sectorial area.
^^krv
dt 2
sin tA = ;^|rXv|=U
2 ^
[ironi Vai (8-22)]
The quantity h/2 called the sectorial, or areal, speed and. according to
is
10-3).
Tocontinue our discussion of central-force motion, we now consider a
specific problem, namely, the motion of a planet around the sun or that of
a satellite around the earth. We consider a system consisting of two
bodies, the sunand one of the planets, or the earth and an artificial satel-
lite. According to Newton's law of gravitation, the mutual force of
attraction between two particles acts along the line joining the two parti-
cles and may be expressed by the so-called inverse-square law as follows:
F=-«^e, (8-23)
m d\
GMm
dt r^
which simplifies to
^-"T^e. (8-24)
gj(vXh)=^^Xh + vX^=^Xh
which, from Eq. (8-2-4), becomes
_d (v
,
X ,
h)
,
= -
GM
— - (e. ,
X h) (8-25)
It is clear that Eq. (8-25) could be integrated if its right-hand side could
Using the expansion rule of a triple- vector product [Eq. (.'3-21)], this
becomes
e, X h = r[(er • v)er — (e^ • er)v] = r[(cr • v)er — v]
dr d , . dcr , dr
''^'dt^d^'^'^^'irt^-dt^'
and thus we have
_ ddr ,
dr _ dr
^
dt dt dt
the last step resulting since Cr isof constant magnitude of unity, and thus
er'(der/dt) =0 (Sec. 6-2). Substituting the expressions just found
into the expression for (cr X h),
/dr \ „ dcr
Cr X h =
,
'' = -T, e
'-{11'''-'')
r I
'-'-lit
This result can now be substituted into Ecj. (8-25), with the result that
V X h = GM(er + t) (8-26)
r • (v X h) = (r X v) •
h = h h = •
/i2
166 Dynamics of a particle
u^here (j) is the angle between the position vector r and the constant vector
E. The last expression may be written
= hyOM
r (8-27)
I -\- e cos
Orbit type
-r COS
:
^ Tg + r^ _ Jj
2 ~ I - e2 (8_2(j^
b = a{l - e-y- = L(l - 6^)-^^
Since the sectorial speed dA/dt = h/2 is constant, the period of an elliptic
orbit (i.e., the time required for a complete revolution) is
2^ _ 27raHl - e^)^^
h h
Substituting h = \/GML,
We now apply the results just obtained to the motion of artificial earth
satellites. To this end it is convenient to reformulate the basic e^iuations
in terms of the gravitational acceleration at earth's surface rather than
the universal gravitation constant and the mass of the earth. We note
that the weight of a particle at the earth's surface is given by
mge =
GMm
—y- I e
shut down so that the gravitational force is the only force acting. (Note :
1
£ = 5
/
(Vo X h) — CrO
h = rXv = roXvo
so that Vo X h = Vo X (ro X vo) = i'o^ro — (vo • ro)vo
where i/'o is the angle between ro and vo. Substituting into Eq. (8-32)
gives the relation we seek:
i/'o ^ sm-
• , ,
i^o + , 1
1 (8-33)
It will be useful now to consider the very special case of a circular orbit
at sea level, which of course cannot be attained because of the atmospheric
resistance but will give an estimate of the required speeds for a low-
altitude satellite. For such a case, we take
ro = Te "^0 = 2
^ =
and Eq. (8-33) gives
,2
=
\9ere }
where, in order to avoid confusion with the general case, the satellite
speed Wo in this special case has been denoted by Ve. Thus
""=*"+' <*-^""
''= fey fey ^""*"-'K^') fey
Given the initial conditions ro, Vo, ^o, at any instant after rocket burnout,
Eq. (8-36) enables us to compute the eccentricity e. From e and
h- _ ro-;'u- si n- i^o
7^1 ~
^ = GM = zr:n
g,re
= ^^ (
- L;) )
-^i"'- "Ao (8-37)
we may compute the apogee and perigee of the elliptic orbit by Eqs.
(8-28). The period of the satellite follows immediately from Eq. (8-30).
As a special case, consider the orbital speed for circular orbits at any
altitude. For such trajectories we take ^o = tt/'Z and c = 0, so that
Eq. (8-36) gives
Vo = V.
^P (8-38)
Thus the orbital speed for circular orbits is inversely proportional to the
square root of the altitude.
Equation (8-36) maj'^ also be used to solve for the so-called escape
velocity,which is really the stfiallcst speed recjuired to achieve an open
orbit. It is clear from Eq. (8-27) that this corresponds to a parabolic
orbit with e = 1, so that for the case of ^o = ir/2, Efj. (8-36) reduces to
fe"yfey-fe)fe"y-
This gives, letting ro/r,, = 1,
Hyperbolic orbit
Parabolic orbit
Elliptical orbit
Circular orbit
Fig. 8-8
F = mr
and form its dot product with the displacement dr:
F ' dr = mr • dr
From Eq. (6-8),
d .. .^ d , ^. .
dr
d&"'^==d&~^^^-'"='-''dt
Hence F • dr = -m-r (?;-) dt
traveling between /^ and P.. The term ^-imv-] is called the kimlic
cncryij of the particle. Hence E(i. (8-39) tells us that the tvork done on a
particle is equal to thechange of kinetic energy of the particle.
We have said before that F in general depends on time, position, and
velocity. It is clear that the evaluation of the
work integral for the
general case and we would have no advantage in deri\ing
is difficult,
F = F(x,!j,z)
then there is hope that the work integral can l)e evaluated l)y analytical
over a rough surface, the friction force between the two surfaces in
contact must necessarily be considered in the ecjuation of work and
energy.
Summarizing what we have accomplished, then We started with the :
''
F . dr
/
t The kinotio pnorgy f.-iri be oxplicitly expressed in terms of any convenient coordi-
nate system. See I'rob. 8-17.
172 Dynamics of a particle
which will in general depend not only upon the positions of the two end
points Pi, Pi, but also upon the path of travel of the particle. However,
for a particular type of force field, called the conservative force field, the
force ¥{x,y,z) such a function of position that the above line integral
is
-
V{Pi) =
P
-
F . dr ViP,) = j^' F . dr
such a conservative force can actually exist and, next, what is the
field
dV = T— dx -f T— dy -\- —- dz
dx dy dz
Hence the potential energy exists and the force field is conservative if
Eqs. (8-41) are satisfied.
From Eqs. (8-41) we can form the following partial derivatives:
dF^ ^ dW dFy _ dW
dy dy dx dx dx dy
For functions that are continuous and have continuous derivatives, the
order of differentiation is immaterial. Hence
^ - ^^ =
dx dij
dF BF
_?___«
Similarly. = (8-42)
dz dx
It can be shown that cither Eqs. (8-41) or Kqs. (8-42) are the necessary
dx dy dz
F = -TV (8-41a)
Ilcnce, in a conservative force field, the sum of the potential energy and the
kinetic energy is a constant. This is the principle of the conservation of
mechanical energy.
The above equations clearly satisfy Eqs. (8-42); hence the uniform
gravitational field is a conservative force field. By Eq. (8-40), we find
that the difference of potential energy between two points at elevations
y\ and iji is
V, - Fi = - !"'
i-W) dy = Wy, - Wy, (8-45)
J yi
the y distances refer to the center of the mass of the body (or system of
particles).
Another example of a conservative force field is the central force field
of Sec. 8-4. If we choose the central point as the origin, then the force
in a central field consists only of a radial component which is a function
of the radius only. This means that in the spherical coordinate system
(Sec. 6-8),
(«-*^)
I: : I'lt
It can be shown that this again satisfies Eqs. (8-42) ; hence this type of
force field is also conservative. The change of potential energy in this
case is
F2 - Fi = - jyr{r)dr (8-48)
One example of the central field is the force of attraction of two mass
particles mi and m-i, respectively, separated by a distance r (Sec. 8-4)
For this case, the indefinite integral on the right side of Eq. (8-48) can
be easily evahiated to obtain
,r ^ m\m->
= -G —^-' ,
The constant in Eq. (8-50) can be taken as zero for a given datum.
Without modification, Ecj. (8-49) gives the vector force of attraction
between any two particles. In order to find the attracting force between
two bodies of finite size it then becomes necessary to add vectorially the
interaction between each particle of each body. That this procedure
is a cumbersome one is (juite evident. However, if we make use of the
concept of potential energy, we can determine the attracting force
between two finite bodies by scalar methods alone. For convenience,
let us further define a poictitial which is equal to the potential energy for
a unit mass. For example, the potential at any point P in a gravitational
field due to a particle of mass mi at a distance r from P, by Eq. (8-50), is
<f>
= —Gnii/r. Subsequently, the force components acting on a unit
mass particle at P due to wi can be found by the use of Eqs. (8-41).
Suppose we desire to determine the potential at a distance r away from
a homogeneous spherical body of radius a, such that r > a. Using the
results of Prob. 8-19, the potential at a distance r from the center of a
spherical shell of radius a is = —Gnii/r, where nii is the total ma.ss of
<t)
situated outside a spherical shell is the same as if the spherical shell were
a particle concentrated at its center. Since the homogeneous sphere
may be regarded as a large number of concentric spherical shells, the
potential at a point outside a homogeneous sphere is <^ = —GM/r, where
now M
is the total mass of the sphere, and thus the force acting on a
unit mass at the point P is F/m = —{GM/r-)eT. This result has been
used previously in Sec. 8-4, where we idealized both the sun and its planet
as homogeneous spherical bodies and suV)sofiuently as particles.
On the point at which the potential is to be computed
the other hand, if
lies within the spherical shell, the results of Prob. 8-19 indicate that the
V = Ukr^ (8-51)
Fig. 8-9
to the length of the string. The
constraining force, then, is the force
in the string. The upper end
of the string is fixed at a point 0. The
pendulum with the amplitude a, which represents the maximum
oscillates
angle of swing, the motion taking place in a vertical plane.
Let us take the potential energy to be zero when the particle is at its
neutral position, i.e., when the string is vertical. For any angle 6
between the string and the vertical, the kinetic energy is 'jy^ma^d- and
the potential energy is mga(l — cos 6). At the extreme position, or
when 6 =^ a, the kinetic energy is zero and the potential energy is mga
(1 — cos a).
The principle of conservation of energy therefore gives
at ^ —
M tig
p .
V 2(co's 6 — cos a)
"^^
pt = r^ —
(8-52a)
Jo \/2(cos 6 cos a)
k = sm
.a- .6
sui = k
.
^m
.
:^ <t>
du = — 2k cos
,
d<l>
.
_, sin (6/2)
sm {a/2)
sin = sn (pt)
or sin .^
= sin - sn (pt) (8-55)
t Values of this function are tabulated in E. Jahnke and F. Enide. "Tal)leof Func-
tions," Dover Publications, New York, 1945.
178 Dynamics of a particle
where K{k) is called the complete elliptic integral of the first kind.t
Without going into details of elliptic functions, we can obtain a great
deal of information about the behavior of a simple pendulum in the
following manner:
1. If the amplitude a is so small that k^ can be neglected when com-
pf =
But sin ^ = sin (0/2)/sin (a/2) ~ d/a for small a, and hence small 6.
This means
e =- a sin 'pt (8-57)
keep the first-order small term. In other words, let us assume that
k « a/2 is a small enough quantity so that only the /c^ term is retained,
whereas higher powers of k {k^ and greater) can be neglected when com-
pared with unity. The integral on the right side of Eq. (8-53) can then
be expanded by the binomial theorem to become
f
*
(1 + y,k-- sin^ 0) d* = * + m- (I
- '^^
In order to find the period of the motion, we substitute (j) = 7r/2 (corre-
2 V ^ 4
This represents the amount of -pt for the pendulum to swing y^ cycle,
t Values of this function are tabulated in V.. Jahnke and F. Emde, "Table of Func-
tions," Dover Publications, New York, 1945.
Sec. 8-9] Example: the s'nnjAe petidnlum 179
since k- = sin^ (a/2) ~ a-/4 for small \alues of a. Thus the period
increaseswith the amplitude. I'or an amplitude of 0.4 radian, or
approximately 22.9°, the error in the period computed from the
elementary, or first-order, solution of Efj. (8-57) is shown by Ecj. (8-58)
to be per cent.
1 Moreover, the error involved here in the second-order
approximation compared with the exact solution of Eq. (8-56) is negligi-
bly small. On the other hand, if the amplitude is 1 radian or more, even
the second-order solution of E(|. (8-58) becomes inaccurate, and still
higher-order approximations, or the use of Eq. (8-50), may become
necessary.
So far we have solved the problem through the use of the conservation
of energy. It is also possible to solve the problem by the equations of
motion. We shall use Newton's law expressed in terms of the normal
and the tangential components [Eqs. (8-3)]. It is most natural to employ
these coordinates when dealing with constrained motion, and it is par-
ticularly advantageous when the path is circular. Referring to Fig. 8-9,
we notice that at any angle 6,
In Eq. (8.59) and for the remainder of this section, T stands for the
ad- = 2^ (cos 6 + C)
This is exactly Eci. (8-52). The solution from here on is therefore the
same as given before.
Let us consider once more the first-order approximation for small 6.
61 + ^0 =
a
Letting p^ = g/a,
d -{- p^6 =
:
The solution is, according to Sec. (8-2) and the prescribed initial
conditions,
d = a sin pt
Tmin ~ mg ll -
^j
Tmax « mgil -\- a2)
especially the force component normal to the motion, since this component
does not appear in the work-and-energy equation. In this respect it is
worthwhile noting that a combination method consisting of the determina-
tion of the speed by energy methods and the normal force by the equations
of motion is particularly efficient.
4. Instead of an exact solution, one can sometimes obtain approximate
Fi<;. 8-10
F = m[R + e,
-j- w/ X p + w/ X (w, X o) + 2o)/ X p,] (8-01)
R = o X (o> X R)
and Eq. (8-62) now becomes
niQr = F - m{o> X [(0 X (R + 9)] + 2a) X pr} (8-63)
= GMm R mw X (w X R) +N
or,writing this in terms of the coordinate
system of Fig. 8-10,
N
GMm «'• mw'^h
R^
Integrating the first and third of these and using the initial conditions
t = 0, x = ]j
= z = x = y = z = 0,
X = 2co/y sin X
i = 2oiy cos X — gt
Substituting these into the second of the e(iuations of motion and neglect-
ing terms in w-,
y = 2cogt cos X
Two integrations yield
y = ojg — COS X
a; =
y = —
2o}h cos X
5 J2h
184 Dynamics of a particle
PROBLEMS
8-1. A mass in rectilinear motion is acted upon by a force such that the displace-
ment is given asx = A cosh wt + B sinh ut. Describe the nature of the applied force.
8-2. Determine the motion of a mass resulting from the action of a repulsive force
proportional to its displacement. The initial conditions for the motion are, at t = 0,
X = xo, X = 0.
8-3. A projectile is shot vertically upward with an initial velocity of vo. If the air
resistance has the magnitude of kv^, find (a) the height reached by the projectile; (6)
the velocity of the projectile when it returns to the ground.
8-4. The resistance of the air to the fall of a raindrop can be considered propor-
tional to the velocity of fall. If the terminal velocity of the raindrop is found to be
20 fps, find the distance of the fall for the speed to reach 99.9 per cent of the terminal
velocity. (Assume that the initial velocity of the raindrop is zero.)
8-5. A particle weighing 100 lb travels horizontally in a viscous medium. The
drag force is proportional to the velocity. It is found that when the speed is 40 fps
the drag force is 60 lb.
(a) Find the distance of travel when the speed slows down from 40 to 10 fps because
of the viscous medium.
(b) Find the time elapsed for part a.
8-6. If a rocket for space travel does not possess any means of propulsion of its
own after it is projected from the ground, find the speed of projection required so that
the rocket can escape from the earth's gravitation. The gravitational force is
inversely proportional to the square of the distance between the object and the cen-
ter of the earth. The radius of the earth is 4,000 miles.
8-7. A particle falls
from rest in a medium that possesses a velocity type of resist-
ance of the form —kx^. Determine the terminal velocity of the particle.
8-8. A shot putter is capable of throwing the shot at an initial speed Vo. His hand
is at a height h above the ground upon release of the shot. Determine the angle <^
from the horizontal at which he should throw the shot in order to achieve maximum
range. Calculate this maximum range. (Neglect air resistance.)
8-9. A particle of weight W is traveling around a parabolic curve as specified in
Example 6-4. Find the reaction between the particle and the surface at its lowest
point of travel.
8-10. A particle traveling in the xy plane is acted upon by an attractive force (or
spring) proportional to the distance of the particle from the origin. Using the rec-
tangular coordinate system, determine the path of the particle. The initial condi-
tions are, al t = 0, x = Xo, y = 0, x = 0, and y = Vq.
8-11. A satellite is launched under the following initial conditions: ro, Vo, v^o = 7r/2.
Classify the various possible orbits as functions of the launching speed.
8-12. The apogee and perigee (measured from the earth's surface) of the following
three artificial satellites are estimated to be as follows:
Compute tlic period of cacli of the tlircc satellite Also coinjjute the orbital speed of
each at the perigee point.
Problems 185
8-13. An artificial satellite is designed to reach a circular orbit at an altitude of
550 miles. During launching the designed orbital speed is met at the re<iuired alti-
tude but the launching angle at rocket burnout is not exactly tangential to the earth's
surface as designed. Instead of a circular ()rbit, the satellite will follow an elliptical
orbit. Calculate the variation of a|K)gec and perigee with launching angle, using the
following values for the deviation of the launcliing angle: 1°, 2°, and '.i°.
8-11. ABCa smooth rail in the form of a vertical semicircle of l-ft radius.
is .\
Show that the work done moving the particle around the triangular
b^- this force in
path OabO is —y^Bab. system conservative?
Is this force
8-16. A force system is given by F = yzH + {xz"^ — 1 j + 'iixijz — l)k. Show )
geneous .spherical shell of radius a. The point is located outside the shell, i.e., r > a.
If r < a, what is the potential?
8-20. For a simple pendulum compute the period of oscillation according to the
methods of Eqs. (8-50) to (8-58) for amplitudi's of a = 30°, 75°, and !)0°. Compare
the results.
8-21. \ small bead of ma.ss in slideswithout friction along a vertical circular curve
of radius a. The l)ea<l starts at from an initial position, nuikiiig an angle of do
re.st
with the vertical axis. P'ind the reaction between the beatl and the wire. Fiiwl the
jM-riod of oscillation of the bead.
8-22. If the mass at the end of a pendulum travels around a circle in the hori-
zontal plane, the device ispendulum. Letting the radius of the circle
called a conical
of travel be r, the length of the string a, and the mass m, find (a) the tension in the
.string; (b) the angular speed 0; (r) the frecpiency and period of the motion.
186 Dynamics of a particle
Fig. P 8-22
Fig. P 8-2:^
8-24. A device called "centrifugal pendulum" has been used for the reduction of
torsional vibrations in aircraft engines. can be illustrated by a simple
Its principle
pendulum AB with one end A attached to an eccentric point that rotates with the
shaft and the other end B holding a mass m. Assume the distances OA = r and
AB = a.
(a) If the pendulum is deflected from the radial direction by a small angle it will (j),
oscillate back and forth about the radial direction. Set up the differential equation
of (f)
for small values of 4>.
8-25. A stretch of the Mississippi River at latitvide 40° during flood conditions is
1.5 miles wide. It is flowing to the south with a speed of 10 mph. Calculate the
difference in level between the two lianks of the river. Wliich bank is higher?
Prublenis 187
8-26. A heavy bead of mass m is free to move on a smooth c-ireular wire of radius a
w hich rotates with constant angular speed w about a fixed vertical diameter.
(a) Find the neutral position (i.e., the angle 6) of the bead for w' > g/a.
{b) If oj- > (j/a, find the period of small crscillations about the neutral position.
Fic. P 8-'_>()
8-27. .\ heavy particle moves with constant speed along a horizontal circle upon
tlicinner wall of an inverted right circular cone who.se axis is vertical. The half-cone
angle is a, and the circle of travel of the particle is at a height /( from the vertex of
the cone. Friction is negligible.
(a)Find the velocity of the particle.
(6) If through some slight disturbance the particle is shifted slightly in elevation
along the wall of the cone, show by differential equation that the particle will oscillate
up and down along the walls. Find the frequency of such small oscillations.
Fic. P 8-27
8-28. A particle at rest is dropped from a height of 100 ft. Find tiic deviation of
full due to the earth's rotation.
CHAPTER 9
Vibrating Systems
(ieiieraliz«Ml
190 Vibrating systems
Fig. 9-2
X = Xm sin o)t
X = Xm COS wt
is identical with that of Fig. 9-2, except for a phase shift of one-quarter
of a period along the abscissa.
motions in general, however, are not necessarilj^ simple
Periodic
harmonic. The motion shown in Fig. 9-1, for example, is periodic but
not simple harmonic. It is a great mathematical discovery that any
periodic motion can be considered as the sum of a series of simple har-
Sec. 9-3] One-degree-of-freedom vibrating systems 191
monic motions. Consider any periodic function /(/) with the period T
so that /(/ + 7')
/(/). =
Subject to certain restrictions, which are of no
concern in the present discussion, one can write
Vw.. \\-.\
wx -I- ci -I- kx = F{i) (9-4)
1. The restoring spring force acting on the mass is —kx, where k is the
spring constant, i.e., the number of units of force required for a unit
deflection of the spring. The negative sign of the spring force means
that the force always opposite to the displacement. This force is
is
damping, i.e., the number of force units of resistance per unit velocity
of motion (cf. Sec. 2-5). The negative sign of this resistance shows that
the resistance is always opposite to the velocity. This resistance is also
called the damping force. We often purposely introduce a damping
force in a vibrating system by means of, for example, a dashpot filled
with oil.
3. The external force F(t), also called the exciting force, is in general
some function of time. The most important kind of exciting force is a
sinusoidal function of time, given before as Eq. (8-10),
where Fo = amplitude
CO = circular frequency
An exciting force that is periodic but not sinusoidal can be written as a
Fourier series of an infinite number of terms, or harmonics (Sec. 9-2)
each term, or harmonic, is a sinusoidal function of time of the type of
Eq. (9-5). Since the differential equation (9-4) is linear, its solution
for a general periodic exciting force can be obtained by adding the solu-
tions for each harmonic of the exciting force.
For the type of exciting force represented by Eq. (9-5), the differential
equation (9-4) may be written in the following standard form
The reason for introducing the new parameters p and n will be made
clear very soon.
The problem of the idealized mechanical vibrating system described
above one of the more important problems of mechanics, partly because
is
the idealized system possesses all the essential features of more compli-
cated vibrating systems, and partly because the differential equation (9-4)
:
can also arise in other entirely ditTerent physical problems, for example,
in the theory of electrical circuits. Let us consider a series circuit
(Fig. 0-4) consisting of a variable volt^igc source E = K(t), a resistance R
(ohms), inductance L (hcnrys), and a capacitance C (farads). If i is the
instantaneous value of the current (amperes) in the circuit, then the
voltage drops (volts) across the three elements are as follows:
Across R : Ri
di
Across L :
J
Across C
This is the differential etiuation of an RLC circuit. It can lie seen that
Ecj. (9-8) can be recast into the form
of Eq. (9-f)) if we let
'iZ
p-
LC
ajid ^ = ^ (9-9) Eit)
If in the system there is neither exciting force nor damping, the differ-
194 Vibrating systems
X -^ pH = (9-11)
This is the same type of equation that we discussed in Sees. 8-2, 8-9,
and 9-2. The solution is a simple harmonic motion represented, in
general, by
X == A cos pt -\- B sin pi (9-12)
Substituting this trial function into Eq. (9-11), there results, after simpli-
fication, the following characteristic equation:
s2 -f p2 =
Thus the trial function constitutes a solution for values of s as given by
Si,2 = dzip, where i = \/— 1. Since we are dealing with a linear dif-
ferential equation, the total solution is composed of the sum of the two
where the two constants of integration A and B (or C and <f>) are to be
determined by the two initial conditions. If at ^= 0, a; = Xo and
X — ±0, we have then
A = Xo 5 =:= -
P
The solution is therefore
Consider the case of a mass at the end of a spring being displaced away
from its neutral position by an initial distance Xo and then released. The
initial conditions for this case are x = Xq and i = at the instant of
Sec. *-)-\] tree vibrations ivilfwul damping 1*).')
The initial displacement .ro in this case is also the maximum displacement
.r„,. i.e., the amplitude of the resulting vibration.
The period of the free vibration represented by Imj. (9-13) or (9-14) is
r = | = 2jf (9-,5)
= 1
= -L /I
/
•^
r f
27r 27r \m (9.10)
^
It can be seen from Eq. (9-13) or (9-14) that an undamped free xibration
lasts forever.
A different method, based on the energy eciuation, can also be applied
to find the period of a free undamped vibration. Assuming that the
motion is harmonic with amplitude ecjual to Xm, we can write that
where the circular frequency oj is as yet unknown. The net force acting
on the mass, in a free undamped vibration about an eciuilil)rium position,
is the spring force —kx, which is conservative. The potential energy
due to the spring force is, according to Eq. (8-51), V = kx-/2. The
kinetic energy of the mass is mx'^/2. The principle of the conservation
of mechanical energy then gives
Consider the following two positions of the mass: the neutral position at
which X = and x = x^o} and the extreme position at which x = Xm
and i = 0. According to Eq. (9-19),
j^niXmW = l-2kx„^
Hence o, = J- = p (9-20)
The frofiuency is therefore the same as that found by solving the differ-
ential equation (9-11).
196 Vibrating systems
It should be noted that the energy method described above gives the
exact frequency of vibration only if the vibration is truly simple harmonic.
Many practical problems are so complicated in nature that the solution
of the differential equation becomes impossible. For such cases, the
vibration can be assumed as simple harmonic, just as we did in
still
^777777777:^
Fig. 9-5
To +U 1 1 +
where k is the spring constant of each wire, Further, the horizontal
component of the force is
r. = To +n
z vi + ix/iy
T. = Tc
I
—X = -2Tof
9 I
or X -p IT/ 7
^
Comparing this result with Eq. (9-11), we see that the mass under-
Sec. 9-3] Free vibralions trilli (hunpiurj 197
X = Ce-" (9-21)
s- -\- 2/(.s + p- =
whose two roots of .s are
where Ci and C2 are two constants to be determined l)y the two initial
conditions.
The physical significance of the
above solution depends on the relative 1.0
magnitudes of n- and p'. We there-
fore distinguish the following three
cases.
Case I. Ifn- > p-, then both Si and
s-2 are real and negative and are given
as
Si,2 = —n±r
Fi«i. 9-G
where r = [x/n^ — p-| is a real, posi-
tive number such that n > r. The two functions e"'' and e'-' are therefore
exponentially decreasing functions of time (Fig. 9-G). The complete solu-
tion of Ecj. (9-28) represents, in this case, a creeping baeli to the neutral
position. The details of the motion can be easily worked out. Suppose
that the initial conditions are such that at = 0, x = .ro and x = 0. <
conditions,
Fig. 9-7
shown (Prob. 9-11), in general, that once the mass passes through the
equilibrium position (zero displacement), it can never pass through the
zero position again and can only approach zero as / °o ^
Case II. A system is critically damped if n^ = p^ or c = 2 \/mk. We
shall call this value of c the critical damping coefficient Cc', i.e..
The solution of a critically damped system is, from the theory of linear
differential equations for repeated roots of the characteristic equation,
conditions.
The nature of the motion in the critically damped case is somewhat
similar to the motion indicated in Fig. 9-7 for the overdamped case. In
other words, the mass tends to creep hack to the neutral position. This
case represents the transition between the overdamped creeping-back
behavior of Case I and the underdamped vibratory behavior of Case III.
Case III. If n^ < p'-, the system is said to be underdamped. We can,
from Eq. (9-22), write
Si,2 = —n ± i \/p- — n- = —n ± iq
1
where i = v—
q = \-\/p^ — n^\
X = e~"'(Cie'"' + C-iC-"")
Using the relations
e'" = cos y -{- i sin y
g-iy — gQg y — i sin y
We call
The natural logarithm of this ratio is called the logarithmic decrement; i.e.,
a = log £.
= 2™
Xb q
ri + a: = (9-32) t
where t = c/k. For reasons that will be made clear shortly, r is called
the relaxation time, or time constant. Thus
r^-^dt =
x
t The student should solve this differential equation by the method of trial
functions.
Sec. 9-6] Tfie case of negiujible mass; (he eqnalion tx }- x = 201
Integration gives
t = —T log X -\- constant
or X = Cie,-"^
The constant of integration (\ is evalualcHl hy tlie initial condition x- Xo
at = 0. The solution is
<
X = x^-'ir (9_33)
The motion of Eq. (9-3.'^) is shown by the curve in Fig. <)-<). We see that
the relaxation time, or time constant,
has the following significance: (1) it is
the time at the end of which the deflec-
tion reduces to \/e of its initial value;
and (2) it is the intercept on the ab-
scissa of the tangent to the curve
Sitt = 0.
applied to the mass and is kept there indefinitely. What is the resulting
motion of the mass? We see that the differential ecpiat ion goxcniing the
motion of the mass, for t > 0, is
F
TX -\- X = (9-:54)
X = Cie-"' +^ (9-35)
We notice that the first term on the right-hand side is the general solu-
tion, or the solution to the homoge-
neous equation tx -\- x = 0, and the
second term is the particular solution,
or that value of x such that tx -f x
equals F/k. The evaluation of the
constant of integration by the initial
condition x = at t = yields the
solution
Fk;. 9-10
The motion of K({. (9-36) is shown in
Fig. 9-10. The deflection is shown to approach its final value/*'/// asymp-
totically. The significance of the time constant t is again demonstrated.
202 Vibrating systems
e:=:
Li + Ri = E
The increase of current with time after
R *S is closed follows the curve of Fig. 9-10.
FiG. 9-11
The time constant in this case is L/R.
3. If in place of the inductance L in
Fig. 9-11 we have a capacitance (capacitor) C, the closing of the switch
represents charging of the capacitor. The charge Q of the capacitor is
Fo
X -^ p- = ^sm coi (9-37)
The
general solution of Eq. (9-37) consists of a homogeneous solution,
i.e., a solution of x +
p-x = 0, plus a particular solution of the non-
homogeneous eciuation. In order to oiitain the particular solution, let
us try the solution .r = x„, sin u)t. Upon substitution of the trial function
into Eq. (9-37), we find that
= ^'
ni{p'^ — 0}^)
F /k
X = A cos pt -\- H sin pt -f -,
1 - —
-. ;-—
(aj/p)2
;: sin mt (9-39)
It is recalled that the first two terms on the right side of the above
equation constitute the solution to the homogeneous differential eciuation
X -|- p'-x = (Sec. 9-4). The last term is the particular solution of the
complete differential equation (9-37). The combined motion represented
by Eq. (9-39) is therefore the sum of two simple-harmonic motions: an
undamped free rihration with fre(juency p and an undamped forced
vibration with frecjuency oj etjual to the forced or impressed frequency.
It might be noted that this combined motion is in general not periodic
unless (1) A and B are zero or (2) the ratio i)etween p and w is a rational
number.
Considering general starting conditions, such that at t = 0, x = Xo
and X = Xo, the motion of the ma.ss is given by
X = Xq cos pt -i sin pt + ,
—— —'wVp , -, ( sin ut sin pt ] (l)-39a)
p 1
\ P /
and consists of both a free and a forccui \ibration. i^ven it the ma.ss
started from rest so that xo = Xo = 0, the solution is
and is seen to be a function of the frequency ratio oo/p. Since Fo/k = 8st
is the static displacement of the mass under the action of a slowly applied
1
|AF
FJk l-(coVp^)
2 \
\
\ V.
1
-1
-2
-3
Fig. 9-12
x^ = 5.,(AF)
1
AF =
1 coVp^
One can from Fig. 9-12 that for w/p close to 0, x,„/(Fo/k) is approxi-
see
mately equal to unity. Physically, this means that for extremely slow
variations, F(t) can be considered to be a static force and the displacement
Sec. 9-7] Forced vihrulions ivilltoiil (hinij)liu/ 20.3
that is nearly the same as the natural, frequency of the system and the
response amplitude .r,„ becomes very large. The system is then said to
he in resonance. Resonance indicates physically that if the disturbing
force acts on the mass at the proper time and in the proper direction, the
displacements of the sj'stem ma}' become dangerou.sly large. Finally,
when co/p is very large, Xm becomes very small and approaches zero as a
limit. Physically, this means that when the exciting force changes
>ign too rapidl}', the mass refuses to follow, and therefore moves only
-lightly, although with the same frequency.
If u = p, the differential equation (9-37) becomes
X -{- p~ = — sm pt
m
The particular solution to this differential efjuation is not periodic as
is the case for o: ^ p. Instead, we have the following particular solution:
X = — ^
2mp
—
Fot
cos pt
Fot
X — A cos pt +B sin })t — ,^—^ cos pt
2mp
Evaluating the constants of integration for general starting conditions,
we find that the total solution is
X = Xo cos pt
, ,
-\
io
p
sm
.
pt
,
+ Imp
,
,^;
^0
— —/sin
I
\ p
/;/
— t cos pi
\
I
The same result would have been oljtained if we had evaluated Eq.
(9-39a) for the condition of = p. The (valuation leads to an indeter-
o>
mx + k(x — X)
where x is the absolute displacement of the mass from m its equilibrium
position. This equation can be rewritten in the form
Assuming now that the free vibrations are damped out, the forced vibra-
tion, from Eq. (9-38), becomes
X = D—sin ost
1 co^/p^
1X = Ds\n ut
»>>»»»: *
"7777777777777-
u = D sin wt
Yp2
Fig. 9-13
1 — coVp^
u ^ —D sin wt = —X
and the recording device indicates the motion of the vibrating body. It
should also be pointed out that under these conditions the vibrating mass
suffers only negligible absolute displacement; i.e., it remains essentially
motionless.
On the other hand, suppose p co. ^
Here the relative displacement is
approximately proportional to the negative acceleration of the vibrating
body; i.e.,
w^D sin wt X
p. V
Thus the recording device now presents a record of the acceleration of the
supporting body. The preceding condensed discussion serves as the
theoretical basis for many of the present-day vibration-measurement
devices and accelerometers.
:
After the preliniinary work done in the i)revious sections, we are ready
now to handle the complete problem of a forced vibration with \ iscous
(lamping. Let us consider the complete ecjuation (li-O)
Fo
X + 2nx + p2x = — sin w< (9-6)
?«
A particular solution to this differential eciuation can l)e found l)y letting
m
2noiN + (/;- - <^-)M =
Solving for M and A'^ and substituting into the trial expression for x,
Fo/m . , >\ , o .1
= + 4nw f^P" ~ ~
'•" -'^"
'^
ip-^
- c.r-
""'^ "' '''' '^'^
where S\ and S2by Eq. (9-22). The first two terms on the right
are given
side of this equation are thehomogeneous solution and represent a free
vibration with damping, as discussed in Sec. 9-5. With damping in the
system the free vibration eventually dies out. For this reason the free-
vibration part of the complete solution is called the transient. For an
underdamped or vibratory system the transient has a frequency of q
[Eq. (9-28)]. The last term of the complete solution, i.e., the particular
:
c/c^=0
Fig. 9-14
The response amplitude x„, (Eq. 9-41) of the steady-state vibration can
also be written in the following form
(9-44)
Fo/k {(1 - coVp2)2 + [2{n/p)(c./p)Yy^^
f
the amplitude factor decreases. Hence danii)inj!; reduces the response
ampHtude. Near the resonant fre(iuency the rethiction is particularly
pronounced, indicating that the severit-y of the resonance condition may
' be materially reduced by the introduction of some damnint:; The I 1-1
I
i .
j!
resonant frequency itself becomes snudler as the damping ratio c/cc
decreases. Differentiating Eq. (9-44) with respect to w and setting the
derivative to zero, we can find the resonant frequency:
}r = y/p'^ 2n'- = p. 1
- (9-45)
KO'
For small damping, as is common in the majority of problems, /('-
is much
smaller than p- and cor ~ p. For such cases, the maximum amplification
factor, from Eq. (9-44), is
Cc
(9-40)
2n 2c
It can also be seen from Fig. 9-14 that for very small or \-ery large
values of u/p, the respective amplification factor is nearly or 1. For
such values of co/p the influence of damping becomes insignificant.
Applying the above results to an electrical RLC circuit (Fig. 9-4), one
can write the transient and the steady-state responses to an a-c source,
E = Eo sin ot, by merely substituting Ecjs. (9-9) into Eqs. (9-43), (9-22),
(9-26), (9-41), and (9-42). Thus, for the underdamped case, we have
Transient
210 Vibrating systems
1/coC and coL are, respectively, the capacitive and the inductive reactances.
The resonant frequency is defined as the frequency at which the steady-
state current is largest. (Notice the difference of definitions of the
resonant frequency in mechanical and in electrical systems.) From
Eq. (9-52) we see that the frequency for maximum im, i.e., the resonant
frequency, is
^^'^^^
EC
at which frequency the capacitive and the inductive reactances are equal.
Thus the resonant current is E/R. It is limited only by the resistance of
the circuit.
Fq sin ut
X = Xm sin (cot —
Hence x = x^co cos (ut
Sec. 9-9] Vibralion isolalion 211
1 + (2 " "Y
(9-56)
i:-^-m\
The ratio Pm/Fo of the transmitted force to the applied force is called the
It can be seen that, although Fig. 9-16 is similar to Fig. 9-14 in many
major respects, there is the following important difference. For the
smallest possible transmissibility, damping is desirable only if oj/p is less
than y/2. For w/p larger than y/2, damping actually increases the
transmissibility. Since it is our desire to reduce the transmissibility to a
value much smaller than 1 (otherwise we might as well use a rigid support,
with transmissibility equal to 1), we want always to operate normally in
the range where u/p > \/2. Damping is therefore undesirable at nor-
mal operating speeds. Some damping, however, is provided for the pur-
pose of avoiding excessive vibration when the motor is near the resonance
speed, which might occur on starting or stopping, at which time damping
is then highly useful.
Example 9-3. Consider now a spring-and-dashpot-niountod motor,
rotating with a constant angular speed co and having an unbalanced weight
212 Vibrating systems
First let us consider the undamped motion of the system. The natural
Fig. 9-17
30 X 12 X 32.2
7.61 radians/sec
200
120 X 27r
= 12.57 radians/sec
60
12.57
= 1.65
7.61
P.n =
2 X 6 X (12.57)'-^
X 0.58 = 2.84 lb
32.2 X 12
Sec. 9-10] TorsioiKil rihnilion 213
c
^ _c_ ^ 3.94 X 3 2.2 X 12 ^ .
''
Cc 2Mp 2 X 200- X 7.61
Therefore, from Fig. 9-14 or Eq. (9-44). the amplitude factor r(Mluces to
AF = 0.42. while from Fig. 9-l() or Va\. (9-r)0), the transmissil)ility
increases to a value of 0.81. The resonant frefjuency, according to Eq.
(9-45), decreases to
and the amplitude factor at resonance, which may occur when starting or
stopping the motor, decreases to a value of .\F = 1.15.
There are many instances when the force applied to the system under
consideration is not sinusoidal or even periodic. Blast- and earthquake-
type forces are well-known examples. Since a general forcing function
can be expressed in terms of a trigonometric series, according to the dis-
cussion of Sec. 9-2, the response due a general forcing function may be
obtained by an extension of the previous methods together with use of the
principle of superposition. However, it is sometimes more convenient to
employ an alternative approach.
Consider a single-degree-of-freedom undamped system acted upon by
some general force. Referring to Eq. (9-6), the equation of motion for
this system is
jo m
Now, integrating the left side by parts, f and after substituting in the
proper limits,
dx
sm
. ,
pt — ,. ©.To
f + f
px cos pi = nF(r) sm .
pr dr
,
-TT
dt Jo
/
m
Performing a similar series of operations on the equation of motion, except
thatnow we start with the multiplying factor cos pr,
dx F(t)
-rj cos pt — Xo + px sin pt = f'
/ —'^— cos pr dr
dt Jo m
Solving these resulting expressions for x by multiplying each equation
by and
cos pt sin pt, respectively, subtracting and simplifying, the response
is given by
Xn 1 f'
X ^ Xo cos pt -\ sin pt H / F(t) sin p(t — r) dr (9-58)
p mp Jo
t Note that
I
Jo dr^
pT
-r-^ sin ^ dr = -7- sin
X^dr
vt \
jO
—piJO -r cos pr dr
dr
This result
is usually referred to as the convolution integral and repre-
sents the superposition of the free and forced vibration components of the
response of the system initially at rest and acted upon by some general
forcing function. easy to verify that the term (1/wp) sin p{t — t)
It is
in Eq. (9-59) represents the response of the system at time t to a unit
impulse acting at time r, for i > t (Prob. 9-21).
Example 9-4. In order to illustrate this method let us consider the
case of the forcing function having the form F(t) = Fo sin cot. Assuming
that the system starts from rest, the response as given by Eq. (9-59) is
X = —
mp
1 p
/
/o
Fo sm cor .sin p(t
^' — r) dr
[t
Fo
X = ;z 1 [cos (iOT -\- pr — pt) — COS (cor — pr -{- pt)] dr
Fo/k
X = (sin oit sin pt
1 — wVP" P /
j
This result, of course, conforms to that given in Eq. (9-396). which was
obtained by another method.
We have discussed in this chapter the theory of the free and forced
vibrations of a single-degree-of-freedom system with linear spring and
damping characteristics. Although such a sy.stem is only a somewhat
idealized model of actual situations, the theory developed can describe to
us the phenomena of resonance, the calculation of natural frequencies, the
reason for the creeping or the o.scillatory motions of transients, and the
concept of vibration isolation. Without any additional theory, many
practical problems can and have been solved in this manner. Other
problems, however, cannot be so idealized without losing the essential
character of the motion. Time does not permit us to delve into such more
complicated problems in this book, but in order to broaden our view, we
;
time only, the differential equation is still linear. A special class of prob-
lems involves a spring stiffness that varies sinusoidally with time. The
differential equation for this type of motion is expressed by
The above equation is equivalent to Eq. (9-10) if, in Eq. (9-10), we let
the spring contacts a curved guide which decreases the effective length of
the spring. Thus the stiffness of the spring increases with deflection,
i'he variation of the spring force with the deflection is shown in Fig.
Force
Fi(i. 0-21
9-206, which clearly indicates that the spring force is not proportional to
the deflection. Figure 9-2 la shows a mass m supported by one set of
springs for small deflections of m. When the deflection of m exceeds a
certain value, other sets of springs come into action, thereby increasing
the total spring stiffness. The variation of the spring force with deflec-
218 Vibrating systems
tion is shown in Fig. 9-216. Both cases shown here are examples of hard
springs, i.e., springs whose stiffness increases with the deflection. A
Belleville spring (Fig. 9-22), on the other hand, is an example of a soft
spring whose stiffness decreases with deflection. This device is used as a
spring washer in electronic equipment.
It hasbeen previously mentioned that at high velocities the resistance
of an object immersed in a fluid is roughly proportional to the square of
the velocity instead of to the velocity. On the other hand, the resistance
due to dry friction (Coulomb friction)
is independent of the velocity (Sec.
Fig. 9-22
but does not hold for the solu-
cients)
tions of nonlinear equations. This
principle states that if Xi is a solution of mx -{- ex -\- kx = Fi{t) and :r2 is
a solution of mx + ex -\- kx = F^^t), then Xi -\- X2 is a solution of
mx -\- ex -\- kx = Fi(t) + /^2(0- The validity of the principle of super-
position rests entirely on the fact that
kxi + kx2 = k{xi + X2) exi -\- ex2 = c(xi -\- ±2)
On
account of the principle of superposition, the complete motion of a
linear system under excitation is a combination of a free and a forced
vibration. This is not true for nonlinear systems. Furthermore, the
natural frequency (of free vibrations) of a linear system is independent
of the amplitude of vibration. With nonlinear systems, however, the
frequency of free vibrations changes with amplitude. We have seen this
phenomenon occurring in the case of a simple pendulum with large
amplitudes (Sec. 8-8). When a nonlinear system is under the action of a
sinusoidal exciting force F = Fo sin oot, the system may exhibit the follow-
ing typically nonlinear characteristics.
1. The amplitude of the response may suddenly jump up or down as
the exciting frequency 00 increases or decreases. The variation of the
response amplitude with co may appear as Fig. 9-23a for a hard spring or
Fig. 9-236 for a soft spring. Both curves are drawn for a constant Fo.
With increasing co, both curves follow the path 1-2-3-4. With decreas-
ing both curves follow the path 4-5-6-1. Thus, in the region 2-3-5-6,
CO,
the amplitude may have two different values for the same co (i.e., ampli-
Problems 219
of Figure 9-23a and b should be compared with Fig. 9-12 to show the
o).
Amplitude
Fig. 9-23
PROBLEMS
9-1. Show that the sum a„ sin iiwt + h„ co.s m^t can also be written A„ sin (nut +
0„). Find A„ and </>„ in terms of a„ and b,,-
9-2. A ma,ss is .supported by two springs of spring constants Ai and h> as shown.
If the spring.s are connected (a) in .series, (/;) in parallel, and (r) in parallel as shown,
determine the resultant spring constant and the freciuency of vibration.
2a
Fig. P 9-3
9-3. Two cylindrical rollers with distance 2a between the two centers roll
widi the
anRular speed w in opposite directions as .shown. .\ uniform bar of length / and weight
\y rests on top of the two rollers. The coefhcient of kinetic friction between tlie bar
and the roller is ma- If the bar is pulled slightly away from its middle position, it will
oscillate back and forth when relea.sed. Find the circular fretiuency of the resulting
oscillation.
220 Vibrating systems
9-4. AU tube is partially filled with water. The total length of the water column
is I. the water in one arm is depressed by a distance x, the water in the other arm
If
will rise by the same distance x. The system is then released. Write the equations
of motion of the water column and show that the motion is periodic. Find the period
of the motion.
/////////////////y///y/^^y/y///y///
'
\
^-^
Vo.Po ^o-Po
^V/////////////////^////A//////y//A
9-7. A cylinder contains gas under pressure po and a piston .4 of weight The W .
enclosed gas obeys Boyle's law; i.e., pF = constant. The initial volume of each half
of the cylinder is Fo; the piston moves freely within the cylinder without allowing any
leakage between compartments. Derive the equation of motion of the piston, assum-
ing it was given some initial displacement. What are the frequency and period of
small oscillations about the position of equilibrium?
9-8. A mass-spring-dashpot system is started in motion with the initial conditions
at i = 0, a; = Xo, =
the system is overdamped, find the resulting dis-
x •'Co- (a) If
constant is /c = 10 tons/in. and the damping coefficient is c = 8 ton sec /ft. A 1-ton
Problems 221
liummer falls on the niacliim- without rolioiiiul, aiul tlio two dfscoiul together after
the impact with an initial velocity of 2 fps. Fiiul
(a) The maximum displacement of the machine due to the impact.
(6) The displacement of the machine at the end of three complete cycles after the
maximum displacement.
(f) The minimum damping to prevent tlie machine from oscillatinK after impact.
9-10. Sliow that for the case of overdampeil free vibrations, a hody never pa.sses
tiiroush the equilibrium position if (a) it is released with initial velocity but no initial
displacement, or {h) it is released with initial displacement but no initial velocity.
9-11. Demonstrate that for the ca.se of overdamped free vibrations with general
starting conditions, once the mass passes through the equilibrium position it can never
pass through the equililirium position again, but can only approach it as < w. ^
9-12. A body weighing 20 lb is suspended by a spring and a dashpot that produces a
resistance of 0.02 lb at a speed of 1 in. /sec. In what ratio is the amplitude of vil)ra-
tion reduced after 10 cycles, if tlie spring constant is k = 5 lb/in.?
(6) The
frequency of oscillation.
(c) The
resonant frequency of the system.
(d) The
amplitude of the force at resonant freciuency necessary to cause the mass
to vibrate at 1 in. amplitude.
9-14. A simple undamped vibrating .system is initially at rest. If a constant force
F is applied to the ma.ss for a time interval t', find the motion of the mass after F is
removed.
9-15. A motor has an unbalanced weight of }4 lb at a distance of C) in. from the axis
of rotation. The rpm of the motor is 1,200. The motor weighs 500 lb and is mounted
at the middle of a platform. It is estinuited that a load of 10,000 lb acting at the mid-
dle of the platform will cau.se a deflection of 1 in. under the load. Assume a damping
ratio c/cc of 0.1.
(a) Find the amplitude of the motion of the motor.
(6) Find the amplitude of the force transmitted to the ground.
(c) How much should the stiffness of the platform be reduced so that the di.splace-
nient of the motor is limited to one-tenth of the value given in part a?
Fi(i. P 9-15
^C
R
I AAV-
Fic. P '.)-IG
222 Vibrating systems
and the charging begins. Write the differential equation for the charge, and find the
expression of the charge as a function of time.
9-17. The inductance of a coil is measured by adjusting a variable capacitor in
series with it The resistance of
to produce resonance at a frequency of 0.5 Mc/sec.
the circuit is 4 ohms, and the capacitance is 0.001 /if (microfarad). The amplitude
of the current at resonance is found by an ammeter to be 2.00 amp. The impressed
voltage and capacitance are now kept constant, while the frequency of the impressed
voltage is changed to 2 Mc/sec. Assuming that the resistance and the inductance do
not change with frequency, find the amplitude of the current.
9-18. Determine the complete solution for general starting conditions, for the case
of forced damped The
forcing function is of the form Fq sin ut.
vibrations.
9-19. A mass of weight W
suspended from a spring of spring constant k and is
is
acted upon by a vertical force Fo sin cot. Determine the range of values of w for which
the amplitude of steady-state vibrations exceeds twice the static deflection caused by
the force Fq.
9-20. A mass-spring vibrating system starts from rest and has applied to it forcing
functions of the form (a) F{t) = ior t = 0, Fit) = Foiort > and (fe) F(t) = F^t/ti
for < t < ti, F{t) = Fo for t > ti. Determine the responses of the system using
the method of Sec. 9-11.
9-21. Show that the response of a simple mass-spring system initially at rest and
acted upon by a unit impulse at time t is x(t) = (l/mp) sin p(t — t), iort > t. A unit
impulse is defined as the product of force F and increment of time dr such that its
resulting magnitude is unity. According to Newton's law of motion, F dr = m dv;
i.e., a unit impulse at t causes a unit change of momentum at t.
CHAPTER 10
We can write one such ecjuation for each of the particles in tlu; system.
Adding up all such equations, we obtain
t'--l'-i\t i>'i,r,]
According to Newton's third law, the internal forces occur in equal and
opposite pairs, thus:
2f, =
Furthermore, by Eqs. (5-6) (written in summation rather than in integral
form)
S(miri) = (2?ni)rc
where re is the position vector of the mass center or the center of gravity
for a constant gravitational field. We thus have, for the usual case of
time-independent mass,
2F,- = (2mi)rc
For brevity, let us use the notation F = SF^- for the resultant external
forceand ?n = linii for the total mass of the system. We therefore can
write
F = mrc (10-2)
This is the important equation that describes the motion of the mass
center of a system of particles and can be expressed in the following form
The motion mass center of a system of particles is identical with that
of the
of a particle ofmass m (the total mass of the system) concentrated at the
mass center and acted upon by a force F (the total external force on the
system). This explains why particle dynamics of Chaps. 8 and 9 is of
great practical importance, for it describes not only the motion of a
particle but also the motion of the mass center of a system, irrespective
of whether the system is a rigid body or not. It is important to under-
stand, however, that this result will not in general enable us to specify
the motion of each particle of the system.
An artillery shell may explode in mid-air into thousands of fragments
whose individual paths are difficult to trace, and yet the mass center of
all the fragments will follow a parabolic path (if air resistance is neglected
and the acceleration of gravity can be considered as constant) before as
well as after the explosion, until one of the fragments hits another system.
F = (mr)
I
can be integrated with respect to time and becomes
[The subscripts 1 and 2 on the right side of Eq. (10-3) denote the values
impulse (or linear impulse) and mir the mmnentum (or linear momentum).
Hence, for a particle, impulse is e(ptal to theThis change of momentum.
theorem, called the theorem of impulj^e and momentum, actually offers
no new information, .since E(|. (10-3) is no more and no less than Newton's
law in a somewhat different form, the tlilTerent form, of course, heiufj; a
first integral. Both impulse and momentum are vector ([uantities;
Eq. (10-3) thereby is a vector e^iuation that can i)e expressed i)y three
scalar equations.
Con.sidering now a sy.stem of particles, we can likewise integrate the
equation of motion [E(\. (10-2)]:
'
F dt = mrc-2 - mbci (10-4)
/,:
where /
'
F dt is the impulse of the resultant external force acting on
the system and mtc is the momentum of the total mass concentrated at
the mass center. Since
2(??i,r,) = mrc
we have, upon differentiating everj- term with respect to time,
particle having mass equal to the total mass of the sijstcm and moving with
the mass center of the system. Then we can express Ecj. (10-4) by the
following theorem: For a system of particles, the impulse due to the total
external force is equal to the change of the total tnomcntum of the entire
syste7n. Alternatively, we may expre.ss this theorem as follows: The
rate of change of the total momentum of a system of particles (or particle)
is equal to the resultant externallv applied force.
The scalar form of Ecj. (10-4J can be written in rectangular coordinates
as follows:
^
F^dt = m{xc2 — ici)
/:
Fydt ^ m{yc2 - yci)
_
/:
F, dt = m{zci — ici)
As a special case of Ecj. (10-4). if the total external force is zero along
any one direction, the linear momentum must be constant in that direc-
tion and the ma.ss center of the .system travels with constant velocity
along that same direction. When the resultant external force vanishes
in all three orthogonal directions, F = 0, we ol)tain the conservation
i.e.,
theorem of momentum, which states that with the absence of external forces,
226 Dynamics of systems of particles
the total momentum is conserved. In this event the mass center of the
system travels along a straight line with constant speed, and if the mass
center is initially remain so. The recoil of guns firing
at rest, it will
horizontally, for example, such that the mass center of the system
is
consisting of the gun and the projectile remains at rest immediately after
the firing. This means that the velocities of the gun and the projectile
must be of opposite direction and of magnitude inversely proportional
to their masses.
Impulsive forces are defined as those forces approaching infinite mag-
nitude and acting on the system (or particle) for an infinitesimal time,
such that their product approaches a definite limit. Thus an impulsive
force Fjr can be defined as
Solving,
mov sin a
R„ = mg ,
-\- m.og +
,
At
moVo cos a
Vc =
Sec. 10-3] Angular niunienhini of a syslem of particles 227
i
3-6).
According to the of Eq.
definition
,
(10-Oa), the origin may
be either fixed
I
or moving. However, to avoid future
confusion, we shall restrict Eq. (10-6a) to
the case where the origin is fixed (in the
F = ^ (mr)
r X F = r X jA^^)
Considering the vector identity [{Va{. 6-7)1,
d d
(r X mr) = r X mr + r X -
.. (mi)
dt
228 Dynamics of systems of particles
r X F = g (r X mr) = _
But r X F is the moment Mo. We thus have for a particle
Mo = ho (10-7)
Equation (10-7) states that the moment of the forces acting on a particle
about a fixed point is equal to the time rate of change of the angular momen-
tum of the particle about 0.
The same statement expressed in scalar form says that the moment
of the external forces acting on the particle and taken about some line
is equal to the time rate of change of the angular momentum of the
particle about the same line, the line in question passing through some
fixed point. The reader is reminded of the equation of motion of a par-
ticle along the tangential direction in cylindrical coordinates [Eqs. (8-4)],
F^r = in -r {r'^4>)
Moi = r, X F, + r,: X f .• = |, Ui
at
Summing up all such equations for all the particles in the given system,
Now 2(ri X fi) = (Sec. 4-1). Let us for brevity use the notation
Mo = S(ri X Fi) for the total external moment about and
Lo = 2Loi = S(r,: X ?n.r,:)
for the total angular moment about 0. We thus have for a system of
particles
Mo = Lo (10-8)
This important theorem states: For a system of particles, the total moment
of the external forces about a fixed origin is equal to the time rate of change
momentum about the same origin. As a special case,
of the total angular
we have the conservation theorem of angular momentum: with the absence
of external moments about a fixed origin 0, the total angular momentum
about of a system is conserved. Notice that the above theorems hold
•. 10-3] Angular momenliini of a sysleni of purlicles 229
lor any fixed origin O about which the moment and the angular momentum
are computed. The moment and the anguhir momentum of course
(hange in vakie when we shift the origin.
For the important special case of a particle, or system of particles,
acted upon by external forces directed toward a fixed origin -the central-
force problem —
the moment about the fixed origin vanishes and angular
momentum is conserved. Consecpiently, the path of the particle or
system is confined to one plane (Sec. 8-4).
It can be shown that the previous theo-
rem does not apply in general to a moving
origin. However, under special circum-
stances and, in particular, if the moving
origin is the mass center of the system, the
previous theorem again holds. To prove
this let r' and r, be the position vectors of
an arbitrary moving origin 0' and the /th
particle, respectively, both referred to a
fixed, or Newtonian, frame of reference
(Fig. 10-2). Let o, be the position vector
Fk:. 10-2
from the origin 0' to the ith. particle; p, is
therefore the position vector relative to 0'. Taking moments about 0'
of the forces acting on the ith particle and summing for all particles,
Mo' =-
y P. X |[m.(r'4-p.)]
Expanding this,
Recognizing that
of the relative angular momentum with respect to the arbitrary point 0'
and is denoted as Lo' [Eq. (10-6c)]. Hence we obtain, finally,
is such that its line of action passes through the center of mass, r' and qc
are parallel and the additional term vanishes. Finally, the third case
occurs when the arbitrary origin 0' coincides with the moving mass
center C, for then pc = 0. It is this last case that enables us to write
the important result
Mc = Lc (10-11)
center of the system. For these two points, the total moment of all
external forces is equal to the time rate of change of the total angular
momentum and relative angular momentum, respectively.
It is interesting to note that when the moving origin 0' is the mass
center of the system,
i.e., the angular momentum about th£ moving mass center due to the
relative velocities is the same as that due to the absolute velocities. We
leave the proof as an exercise (Prob. 10-7).
another interesting result that appears as a by-product of the
Still
But the last term can be rewritten according to Eci. (10-12). Further,
imagining the entire mass of the system concentrated at the mass center and
moving with it, and the angular momentum of the system about the mass
center due to the relative motion with respect to the mass center.
The vector equation (10-8) or (10-11) can also be written in terms of
its three components, in general form, as ff)llo\vs:
Along y :
y {zF. - xF.) = y m(zx — xz)
P' d
(mr-co) dt = mr^co
jto dt
= m(r rcoi — ro-coo)
ri-= (0 X r, = + wl sin 6 k
where r,-= + (Z sin i + cos / j). It follows, then, that the total
angular momentum becomes
ho = 2ri X mr-i = 2mri X r,- = 2??7paj sin d {
— cos 6 i -\- sin 6 j)
Mo = 2ml- sin 6 (
— 6) cos i + cj sin j -f co- cos d k)
I'p to this point, we have discussed the motion of particles and systems
of particles assuming that the mass of the system remained constant.
A' m K„
234 Dynamics of systems of particles
The bracketed term on the left represents the total linear momentum
first
at the time t -\- At; the second bracketed term on the left represents the
total linear momentum at the time t; their difference represents the
change of the total momentum during the time interval At. We have
assumed above formulation that after entering or before leaving,
in the
the additional mass in the system takes on the velocity of the system.
(This assumption is made merely for convenience. It is not a necessary
assumption.)
Expanding Eq. (10-14) and neglecting higher-order terms, we obtain,
after passing to the limit as A^ -^ 0, the following expression for the
motion of the system with variable mass:
7n ^ + KoVo -
at
K^Vi = F (10-15)
For those cases in which the velocity vectors in Eq. (10-15) are approxi-
mately parallel (i.e., rectilinear motion), the vector equation can be
reduced to a scalar equation in the direction of the motion:
where nio is the initial gross weight of rocket including full load of fuel
and is measured from the instant of firing. Further, if we neglect the
t
dv ,
T
dt= -^^m
or. substituting the expression for m and rearranging,
dv ^(-g-\- ^^ dt (10-20)
Assuming that n, T, and g are all con.stant during the period of rocket
operation, Eq. (10-20) can be easily integrated to yield the speed of the
rocket at any instant t after firing:
y = - log - gt
n mo — ,
fit
v, = I'log^'-^/. (10-20a)
Thus, in order to increase vt, we have to increase both ratios+ T/n and
mo/mfe, the latter being the ratio of the rocket weight with fuel to its
weight without fuel (i.e., the structure weight plus payload). With the
present chemical fuels and rocket-structure design, it has been found that
Vb does not quite reach the .satellite speed (Sec. 8-5) in a single-stage
rocket. Multistaging therefore becomes necessary for satellite launching
at present.
As a second example, consider the case of jet propulsion, say a jet plane
of mass m flying at a speed v. In this case air enters and leaves the jet
engine with Ko larger than Ki by a small percentage, the difference due to
the fuel added. If we neglect this small difference and again let A be the
discharge area of the jet nozzle, pe the discharge pressure, and Ve the dis-
charge jet speed relative to the engine, then
Ui = —V Uo = —Ve Ki = Ko = peAVe
where F' has the same significance as in Eq. (10-18). The quantity in
I
the brackets of Eq. (10-21) is commonly referred to as the thrust. Thus,
for jet propulsion, the thrust is given as
2
j' (F, + f,) . dvi = 2 y^m^Vi' I'
(10-23)
The left side of this equation represents the total work, and the right side
represents the change of the total kinetic energy.
Consider the total work on the left side of Eq. (10-23), written with
respect to the arbitrary moving point 0' (Fig. 10-2).
= r (I ''O
•
^^'
+ /i' ll ^0 •
^^' + A' I ^^'^ + ^^^ • ^'^
I
Because 2fi = 0, the second term vanishes. The total work therefore
becomes
W = f'F' dr + f' y^
(F, + fd .
dgi
I
where F = SF^ represents the resultant external force acting on the sys-
tem. In particular, if the point 0' coincides with the mass center of the
Equation (10-24) expresses the fact that the total work done by the
forces acting on a system of particles is etiual to work done by the external
forces due to the displacement of the mass ccMitcr, plus the work done by
all the forces (internal and external) due to the displacements relative to
the mass center and due to rotation about tiie mass center.
Consider now the total kinetic-energy term on the right side of Eq.
(10-23)
T = i:}2miv:- = 2i2m.(r' + p.) (r'
• -\-
p.)
= i.ym.r'^ -\- f' • 2w,p. + r>._,'".pr
Again, if the point 0'
is chosen coincident with the center of mass C, p, is
measured from the mass center, 2/«,p, and its derivative wm,p, are both
zero, and the second term therefore vanishes. Thus the total kinetic
energy can be written
Hence the total kinetic energy of a system of particles is the sum of the kinetic
energy of the mass center {as if the total mass of the syste77i were concentrated
there) and the total kinetic energy of the particles due to the relative motion
with respect to the mass center.
Equation (10-28) can now be written
This efiuation can be separated into two ecjuations. Taking the dot
product of E(i. (10-2) with drc and integrating (see Sec. 8-0),
We see that this is one part of Ec]. (10-2()). vSubtracting this equation
from Eq. (10-2G), we obtain the other part of Eq. (10-26):
^- (miA- •
i-i) = rriiV^ + r^ • (F,- + f,) (10-29)
Q = - \
Qdt (10-30)
Applying this definition to Eq. (10-29) and taking the time average
over a very long time (i.e., r ^^ oo), we find that the left side of Eq.
(10-29) averages to zero if the particle moves in a bounded region for
which r and r are finite at all times. Hence
TYliVi^ = -Xi' (Fi + fi)
Let us divide the above equation by 2 and sum up over all particles in the
system
HZ;^^ = -y^-Zvi' (F, + fO (10-31)
We now seek to apply the virial theorem [Eq. (10-31)] to a perfect gas.
Since the interaction between molecules that are not in direct contact is
neglected in a perfect gas, the virial of the system becomes that due to the
external forces alone. To be sure, there exist also internal forces of action
and reaction at the instant two molecules, but such a pair of
of collision of
forces, besidesbeing e(iual and opposite, act at the same position r antl
therefore contribute nothing to the virial. Hence
2r. • f, =
and E([. (10-31) becomes
1
2^m,rr = - ' 2i>, • F, (10-32)
and Xi = a on wall B,
pF = RT (10-36)
mv 2 = 3kT
Thus pV = NkT (10-38)
or alternatively, p = nkT (10-39)
where p = mN/V is the density of the gas. We thus are able to find the
root-mcan-square (rms) velocity of the molecules:
where /?, = R/mN is the gas constant per unit mass.f We shallhave
occasion to refer to this ecjuation when we discuss the velocity of sound
in a gas.
Fig. 10-7
for a gas on the basis of elementary kinetic theory of gases. The practical
aspect of the problem (for example, the influence of vi.scosity on fluid
motion) will be discussed in Volume 2.
^ = li.T
p
The 8ub.script s is often omitted in thermodynamics with tlic unthTstanding that the
gas constant in the c(|uation is that for a unit mas.s.
242 Dynamics of systems of particles
while the lower plate is held stationary, the speed of the fluid will vary
element ahcd of the fluid at one instant will assume the position ah'c'd
some time later. Let A be the area of one plate and y be the distance
between the two plates. It is found experimentally that for most com-
FiG. 10-8
mon fluids the shear stress FyA is proportional to the velocity gradient
U/y; i.e.,
F U (10-41)
A "^
y
equal tendency to travel in all six directions (i.e., ±.r, ±y, ±z), is nAvJQ>
per second. These molecules come on the average from a region which
is at a distance of L above A and at which the fluid velocity is greater than
F = eA^f
3 dy
(10-44)
'
pVsL
(10-46)
place, too many and too drastic assumptions are made in connection
with the structure of the particles and the forces of interaction between
them. Secondly, the elementary theory is not based on a rigorous
mathematical foundation. It is the purpose of statistical mechanics to
reduce the number of assumptions to as few as possible and to formulate
and apply the mathematical theory of probability in a precise manner.
For example, one of the fundamental problems in statistical mechanics is
the ergodic problem, i.e., the problem of determining the necessary and
sufficient conditions under which a time average of a single system can be
replaced by an ensemble average of a huge collection (or ensemble) of
replicas of the system. Another problem is the computation, based on
the so-called limit theorems, of the ultimate statistical properties of a
system possessing a large number of particles or degrees of freedom. It
can be appreciated that mathematical subtlety and precision are neces-
sities in this branch of science.
PROBLEMS
10-1. A .shell explodes at the hif^hest point of its path into two .splinters of equal
mass. It is observed that one .splinter falls down vertically. W hat is the path of the
other .splinter? Neglect air resistance in this prol)lem.
246 Dynamics of systems of particles
10-2. Neglecting friction and assuming massless pulleys and cord, find the acceler-
ation of the two weights Wi and W2 and the tension in the cord.
10-3. A 2-ton automobile traveling due east at
60 mph collides with a 5-ton truck traveling due
south at 40 mph. If the two vehicles lock together
after collision, find
(a) The velocity after collision.
(6) The velocity of the mass center before and
after collision.
(c) The kinetic energy of the two vehicles before
and after collision.
10-4. The impact of two bodies whose initial
velocities are along the same straight line is called
direct impact. the line of action of the impact
If
vi and Vo (see figure) are the speeds of the two bodies before the impact and v^ and v^
are the speeds after the impact, then the coefficient of restitution is
Vi — V2
Thus, for purely elastic bodies, e = 1, and for purely inelastic bodies, e = 0. If two
bodies of masses Wi and nio travel with speeds Vi, v-i along the same straight line as
shown, with Vi > v-i:
(a) Find the speeds Vj and v^ of the two
bodies after impact, using e as the coeffi-
cient of restitution.
(6) Find the loss of total kinetic energy
of the two bodies due to the impact.
(c) Evaluate the solution of cases a and Fig. P 10-4
h for nil = m-i = to and e = and e = 1.
{d) Evaluate the solution of cases a and h for TO2 = <» and e = and e = 1.
10-5. A pile 25 ft long, weighing 40 lb /ft, is being driven into the ground by a ham-
mer of 3,000 lb weight. With the last blow, the hammer is dropped freely from a
height of 16 ft and the pile penetrates 3 in. deeper into the ground. The coefficient
of restitution between the hammer and the pile is 0.5. Calculate the carrying capacity
of the pile. (Assume that the resistance from the ground to the pile is a constant dur-
ing the 3-in. movement after the last blow.)
10-6. A acted upon by a force which is always directed toward a fixed
particle is
/ —wvw
10-9. particles of masses m and 2m, respectively, are connected l)y a rigid bar
Two
of negligible mass and of length /. The system is initially at rest in the xy plane. A
force F of constant magnitude is applied to the mass m as shown; F lies in the xi/ plane
and is always normal to the bar. Find the angular velocitj' and the angular displace-
ment of the bar after a time interval t.
10-10. Two particles of masses m and 2m are moving with con.stant speeds v in
the xz plane and -\/5 v in the xy plane, respectively. At the instant under consider-
ation, m is located at point (2,0,2) and 2m is at point (1,1,0).
(a) Find the velocity of the mass center of the sj'stein.
(6) Find the angular momentum of the system about the origin.
(c) Find the angular momentum about the mass center (due to velocities relative
to the mass center).
(d) What is the physical significance of the dif'l'erence of {!>) and (c)?
C Jo
mi
a
248 Dynamics of systems of particles
10-12. Two particles of mass nii and m2, respectively, move about so that at a cer-
tain instant the relative velocity between the two is Vr and the velocity of their mass
center is Find the total kinetic energy of the system.
vc.
10-13. A pendulum (for determining the velocity of a bullet) consists of a
ballistic
bag of sand of weight W
suspended from a string of length I. The bullet, of weight w,
is fired horizontally into the sand. Find the velocity of the bullet if the pendulum
swings out to an angle a.
V >
W
Fig. P 10-13
10-14. Suppose, in Prob. 10-13, that the bullet is fired downward into the bag of
sand at an angle (3 with the horizontal. Find the speed of the bullet if the pendulum
swings out to angle a.
10-15. A box of sand slides down a plane of inclination a. with the horizontal. If
the mass of the box and its contents is equal to m at some time t and the sand spills
out of the box and continues down the inclined plane, find the equation of motion of
the center of mass of the box of sand. Assume that friction is negligible.
10-16. Suppose, in Prob. 10-15, that the sand is released in such a fashion that it
remains at rest after release. In this situation write the equation of motion of the
center of mass of the box of sand.
10-17. In Prob. 10-16, the mass of the box of sand and its contents at time I =
is mo. The sand escapes at rate k. If the box starts from rest, find its speed at some
later time t.
10-18. Given a rocket. Explain how you would conduct an experiment to measure
its thrust.
10-19. Under the conditions specified for Eq. (10-20), what physical conditions
should be striven for in order that a given rocket achieve maximum possible speed?
For this ideal condition, and assuming that the propellant has a specific impulse of
200 sec, determine the ratio of payload to weight of fuel necessary to achieve satellite
speed and escape velocity.
10-20. Estimate the rms velocity of hydrogen molecules at 15°C, considering the
density of hydrogen as one-sixteenth that of oxygen.
10-21. Suggest the dependence of the mean free path of an ideal gas on the pres-
sure and temperature. Estimate the mean free path of oxygen at 100°C and a pres-
sure of 0.001 atm.
CHAPTER 11
r = x-i 4- //j + 2k
The mass moments of inertia of the body with respect to the axes x, y,
and z are defined asf
= /^,
(//- + 2-) dm 3=
1^,
(r2 - .t2) dm
= 1^,
{z'- -f x'-) dm = j^,
(r^ - //2) dm (11-1)
= f (.t:2 -H //2) dm = f .
(r2 - z^-) dm
t Note the slight difTerence of nomenclature used here compared with that used in
Sec. 5-4. Here we use double subscripts to denote moment of inertia, whereas before
we used the single-subscript notation. The reason for the double-subscript notation
will become clear when we discuss rotation of coordinate axes in Sec. 1 1-4.
249
250 Moments and products of inertia
where the integration is carried through the entire volume of the body.
The quantity Ixx is termed the mass moment of inertia about the x axis,
lyy the mass moment of inertia about the y axis, etc. In the expressions
of Eqs. (11-1), we see that the integrands (y^ + z^), (z^ + x^), (x^ + y^)
are the distances of the mass element dm from the x, y, and z axes,
respectively. The mass moment of inertia of a body about an arbitrary
axis is therefore the summation (or integration) of the products of the
constituent masses of the body and the square of their distances to that
axis and thus can be thought of as a second moment of the mass about the
axis in question. According to this definition, we immediately notice
that the mass moment
of inertia must always be a positive quantity
greater than zero regardless of the shape of the body, its mass distribution,
or the location of the axis. It might be added that our definition of mass
moment of inertia is in agreement with what is usually given in elementary
physics, in which I alone often stands for /^^ here.f
The mass products of inertia of the body with respect to the coordinate
axes X, y, and z are defined as
Ixy = j^ xy dm
(11-2)
Again, the integrations are to be extended over the entire volume of the
body. The quantity Ixy is called the mass product of inertia of the body
with respect to the x and y axes (or more precisely to the two planes whose
normals are the x and y axes, respectively) and is simply the summation
(or integration) of the products of the constituent masses and their x and
y coordinates. Similar statements hold for ly^, etc. Unlike mass
moment of inertia, mass product of inertia is an algebraic quantity which
may be greater than, equal to, or less than zero.
can be seen from the definitions of Eqs. (11-2) that Ixy = lyx,
It
lyz =
Izy, Izx = I xz- The three mass moments of inertia and the six mass
products of inertia can be conveniently written in an array as follows:
Sec- 11-1] Introduction and dejiniliona 251
body, then Eqs. (11-1) and (11-2) can be written in the form
/xx
252 Moments and products of inertia
k) of the body about any axis is defined in such a way that mJc^ is equal to
the moment body about that axis. Thus the moment of
of inertia of the
inertia of a body remains the same if we replace the body by a particle
having the same total mass as the body and located at a distance k from
the axis.
Let us suppose that with respect to a given set of axes xyz, all the
moments and products of inertia hx, hy, etc., of a given body are known.
-, We wish to obtain the moments and
products of inertia I'^^, I'^y, etc., with
'^
respect to a new set of axes x'y'z' ob-
tained by a translation of xyz (Fig.
11-2). Let the origin 0' of the new
system be located at {a,b,c) in the xyz
system. We thus have the trans-
formation equations of the two coor-
dinate systems:
X = X —a
y' = y -b
z' = z — c
Fig. 11-2
By definition [Eqs. (11-1)], the
moment of inertia about the x' axis therefore becomes
where yc and zc are the coordinates of the mass center in the xyz system
and are given by Eq.s (5-6). Similarly,
where d is the distance between the two parallel axes x and x'. Equations
(11-6) and (11-7) are sometimes referred to as the parallel-axes theorem of
Sec. 11-3] Moments and producis of inerlia of some simple bodies 253
the moments and products of inertia. The special case of a plane area
has been treated by Eqs. (5-22) and (5-23) in Chap. 5. It must be
emphasized that Eqs. (11-6) and (ll-J) are valid only xijz axes
if the
are the center-of-mass axes, i.e., axes having their origin at the mass
center. Accordingly, moments of inertia taken with respect to axes
passing through the center of mass (or centroid) are indicated by an
asterisk in Eqs. (11-6) and (11-7), thus emphasizing this point. Ecjuation
(11-6) can be expressed by the following statement: The moment of
inertia of abody with respect to any axis is equal to the sum of (1) the
moment of inertia with respect to a parallel center-of-mass axis and (2)
the product of the mass of the body and the square of the distance between
the two axes. A similar statement holds for Eq. (11-7). Since each
term on the right side of Eq. (11-6) is positive, we can conclude that of
all parallel axes, the moment of inertia taken about an axis passing through
-b/2 Jr"
n.- pT"
^ J -c/2 Jr"
{<f+z^)dxd,dz '
-a/2^-^
^ pabc{b'^ -f c^)
12
mjb^ -H c")
Hence I* =
12
254 Moments and products of inertia
For the special case of a thin rectangular plate in the xy plane, we take
c = 0. Then
I* = mb^ ^^
/* =
ma"^ ,^ ni(a'^ + h^)
12 ~Y1 12
Circidar Cylinder {Fig. ll-4a). Let the radius of the cylinder be R and
the length be I. Again we choose the x, y, and z axes to be the centroidal
(b)
Fig. 11-4
(//* = (.1-2
+ ,^2) fifn ^ 2iTrHp dr
/* = 27r/p / r' dr =
^
Sphere. Again we choose the center of mass, in this case the center
of the sphere, as the origin of our coordinate axes. By definition, or
Eqs. (11-1), we have
/* 4. /* 4- I* = 2 [ r- dm
where r- = x- + //" + z-. In the case of a sphere, all three axes x, y, and z
are diameters. Hence
= = =
-^'AIr-'dm
/* /* /* 2,
Dividing the sphere into thin spherical shells of radius /• and thickness dr,
the mass of one such shell is dm = -iirr'-p dr. For the entire sphere of
radius R, we obtain
Thus all products of inertia are zero with respect to the particularly
chosen axes. Generalizing this, we can see that the products of inertia
1^
256 Moments and products of inertia
(mass, volume, or area) vanishwhen taken with respect to any two axes
whose symmetry, provided one of the two axes is
origin lies in a plane of
perpendicular to that plane. For a nonhomogeneous body, a plane of
geometrical (volume) symmetry is not necessarily a plane of symmetry
with respect to mass. For plane areas, these statements imply that when
one axis about which the product of inertia is to be found is an axis of
symmetry, the product of inertia vanishes.
Often in practice, we are faced with the problem of determining the
moments and products of inertia of irregular bodies, bodies composed of
several portions. In this event one simple integral will not suffice to find
the moments and products of inertia. However, a method of composite
bodies will enable us to solve the problem. According to this method,
the moments and products of inertia of a body composed of n portions
may be determined by expressions of the form
n n
-i XX / J-
-t XX;
xxi -l
-^ xy 7
/ -* xv: ' ' '
xyz). The unit vectors along these axes will be called ii, io, ij (instead of
the usual i, j, k). Similarly, the new rotated set of axes will be called
whose unit vectors are ij.
x'lXlr's, ij if,, (Fig. 11-5). In this new notation
the array of Eq. (11-3) is rewritten as
In
- (/) (11-8)
ky.xz
258 Moments and products of inertia
values of the complete set of Iki's, where the indices i, j, Jc, I can be any one
of the three numbers 1, 2, or 3. Before undertaking this task, we shall
first consider the simpler but more basic problem of describing the
components of a vector in terms of rotated coordinates.
In this section we shall consider the following problem : Given the three
rectangular components Xi, X2, Xz of a fixed vector r. What are its
components x'l, x'2, x'^ when the coordinate axes are rotated?
Let us use the symbol {x\,x^ for the angle between the positive direc-
tions of the axes x\ and Xj. We also adopt the notation = cos {x'i,x^.
a,-,-
Thus the nine a^-'s represent the nine direction cosines between the two
sets of axes x']_,x'2,x'^ and x^^x^.x^. They can be written
an =
or, in general,
Sec. 11-5] Orlltugonal IransJ'urnialions 259
»i = 2"^'»' iJ = 1. ^. 3 (11-17)
j
The
nine elements of the tran.sformation matrix of Eq. (11-18) are not
allindependent because the unit vectors in each set of coordinate systems
are mutually perpendicular and have unit magnitude. In other words,
we must have
i[ ' i'o — io • i'n = i'-, • i[ —
with similar relations for ii, U, 1,^. Substituting Va\. (11-15), or its
e(iuivalent, Eqs. (11-14), into the six equations (11-18), we obtain the
following six equations:
The first three equations of (11-19) are called orthogonality relations since
they come from the fact that i^, i2, is are mutually perpendicular. The
la.st three equations of (11-19) are called normalization conditions; they
come from the fact that ij, i2, is are each of unit magnitude. All six
equations can be represented by the following simple ecjuationit
t After the reader is accustomed to this type of notation, it is easy to derive Eq.
(11-20) without having to go through the process of writing out Eqs. (11-19). Let
us take Eq. (11-15):
"i = ) ««;«j
Similarly, \[ = y Okiii
Forming tlie dot ])n)d\ict of the two, rcinciiilx-riiig that the dot product is associative,
we obtain
i'i
•
i'k
= 5,t = ) ) (lijUiuij ' ij
but ij • i; = 8ji
J
which is the same as Eq. (11-20).
260 Moments and products of inertia
It is possible to show that Eqs. (11-20) and (11-21) are equivalent. The
nine elements ajj, therefore, are related by the six equations of (11-20) or
(11-21), and only three elements out
of the nine are independent. This,
with our previous knowledge that the orientation of a
of course, agrees
rigid body with one point fixed is determined by three independent
variables (Sec. 7-3).
Let us now return to the original problem we have proposed to solve.
Consider any vector r whose components are Xi, x^, Xz in the unprimed,
or original, coordinate system and x[, X2, x'^ in the primed, or rotated,
coordinate system. Forming the dot product of r with Eq. (11-15) and
noticing that r • i'- = x'^ and r = Xj, we obtain
• i_,
Xi = yajix'j (11-23)
the result is
But the coefficient of the last term can be reduced in the follow ing manner:
k I
(11 -24a)
(11-246)
This completes our task of finding the moments and products of inertia
with respect to a rotated coordinate system. I^(iuations (ll-24a) and
(11-246) are sometimes called tensor transformations; they provide us
with a means of computing the moments and products of inertia with
respect to a rotated set of axes in terms of the moments and products of
inertia with respect to an unrotated set of axes, or vice versa.
To Eq. (ll-24a), we shall write out the complete
illustrate the use of
erjuation for the following two elements of /'y:
Referring back to the usual xyz and x'y'z' notation and recalling that
J\i = I XX, 1 12 = —Jxy, etc., we can write out the moment of inertia /'„
:
and the product of inertia I'„j with respect to the rotated coordinates as
follows
Fig. 11-6
Equation (11-25) tells us that the moment of inertia 7^^ of the body
about any inclined axis x' passing through that point is dependent on
the three direction cosines of x', i.e., an, ai2, and au. Let us now lay off
on the inclined axis x' (which passes through the same origin 0) a distance
OP which is numerically equal to l/\/7^. If the three coordinates of P
with respect to the unprimed system are denoted by x, y, and z (Fig. 11-6),
then by this construction
x = aiiOP = an/^I'xx
y = av£>P = aii/Vl'xx
z = auOP = au/\/Pxx
Solving the above three equations for an, Oio, an and substituting into
Sec. 11-8J Principal monienis of inertia and principal axes 263
'
,
This is the eciuation of a (luadiic surlWe (Fig. Il-(i6). It is, in fact, the
equation of an ellipsoid about the center since the coefHeients of the
I
first three terms are positive here. The locus of the point P, therefore,
is an ellipsoid, called the inertia ellipsoid by Cauchy (1827). It is also
known as the momental The geometrical interpretation of the
ellipsoid.
ellipsoid is useful and enlightening. Thus the moment of inertia about
any line L passing through is numerically equal to the reciprocal of the
square of the distance from the origin to the point where the line L
intersects the ellipsoidal surface.
The concept
of the inertia ellipsoid is an important one. (As a matter
of the geometrical representation of a physical quantity by a
fact,
quadric surface is a useful device appearing frecjuently in various portions
of mechanics.) Our original unprimed axes xyz may, after all, be ori-
ented arbitrarily with respect to the body. With each different orienta-
tion of X, y, and z, we shall have different values of Ixx, Ix,,, etc. Thus
Eq. (11-27) will look different, since all its coefficients will be different.
Yet all such different-looking equations will describe the same ellipsoid
with respect to the body. In other words, the inertia ellipsoid, represent-
ing the inertial properties of the body at 0, is an invariant of the body and
is independent of the orientations of the coordinate axes, provided the origin
where /j^, /,„ and /j-j- are the moments of inertia with respe(;t to the
principal axes. They are called the principal moments of inertia. The
planes formed by the principal axes are termed principal planes. Since
the moment of inertia about any axis is numerically equal to the reciprocal
of the s(iuare of the radius along that axis, from the origin to its point
:
larly, the smallest principal moment of inertia is also the minimum among
all the possible moments of inertia of the body at 0. Every rigid body
therefore possesses three principal axes, corresponding to the three
principal diameters of the inertia ellipsoid. For the special cases where
the inertia ellipsoid becomes an ellipse of revolution or a sphere, there are
more than three principal axes (and principal moments of inertia), each
corresponding to one of the principal diameters of the inertia ellipsoid.
Comparing Eq. (11-28) with Eq. (11-27), we see that with respect to
the principal axes ^, rj, and f, all six products of inertia (/j„ etc.) vanish.
The moments and products of inertia with respect to the principal axes
therefore appear as
//,, \
/„ (11-29)
\0 /,,/
etc., where an, etc., are the direction cosines between the x'y'z' and the
^77f axes.
We see that Eqs. (11-30) and (11-31) are considerably simpler than
Eqs. (11-25) and (11-26). The computation of moments and products
of inertia about rotated axes can therefore be simplified if the unrotated
Sec. 11-8] Principal moments of inertia and principal axes 263
axes are the priiioipal axes. It is theicl'ore useful to locate the principal
axes and to compute the moments of inertia. To locate the
principal
principal axes of the ellipsoid such as the one defined l)y 1m i. (1 1-27) is a
problem of geometry. The derivation is somewhat lengthy, and we shall
merely state the result in the following theorem.!
Theorem. If a, 13, and y are the direction cosines of any one of the
three principal axes with respect to the given xyz coordinate system, the
following three equations must then all he satisfied:
(/xx -
- I, - /,,7 =
\)a ,13
^xx X ~I XI, — I IZ
^ i/X I III/
A 1 y, = (11-33)
i zi i zii i zz — A
general, three roots of X. It can be shown that the three roots of X Xi, —
X2, and X3 — areand are just the three principal moments of
all real
inertia Ijj, I,,,, and Iff. For each value of X, say Xj, we can solve for
the ratios ^i/at and 7,/ai by using any two of the three Eqs. (1 1-32). It
is noted here that onl}^ two of these three equations are independent by
virtue of Eq. (11-33). These two ratios together with the equation
tti^ + iSi^ + 7i~ = 1 are sufficient to give us one set of values of a„ /3,, 7,.
The direction cosines thus obtained are the ones belonging to that par-
ticular principal axis with respect to which the moment of inertia has the
assigned value of Xj. Thus, if we substitute Xi = /^^ in Eqs. (11-32),
the resulting ai, |8i, 71 are the direction cosines of the principal axis |.
Taking successively the three values of X, we can find all three sets of
direction cosines belonging to the three principal axes. As stated before,
the three principal axes are mutually perpendicular.
We thus have reached the following objective: Starting with a given
set of axes xyz and the array
/xx
(/)
=
t The discu.s.sioii in this .section, including thi.s theorem, turns out to be the con-
sequence of some more general theorems on symmetrical UiiHors. A rigorous proof of
this theorem will therefore be given in its proper place in Chap. 14 (Volume 2).
JSee Prob. 11-7.
;
I
S.r. 11-8] Principal ni(tnieiils of inerlia and princijnil axes 267
Using the results of Sec. 11-3, we find the moments of inertia about
the center-cjf-mass axes as
Substituting the value of Xi into the first two of Eqs. (11-32) and solving,
we obtain the following ratios of direction cosines:
^ = 1 :^ = i
the three direction cosines of the first principal axes ^ are found to be
„, = ^, = ^, = -L
Thus the principal axis ^ coincides with the body diagonal u of the cube.
Since the other two principal moments of inertia are equal, it is clear
that the inertia ellipsoid an ellipsoid of revolution whose axis is the
is
body diagonal. Therefore any axis at lying in the plane normal to the
^ axis is a principal axis. From symmetry considerations alone, the
diagonal plane passing through is a plane of symmetry, and thus two
principal axes lie in this plane; the third principal axis is perpendicular
to it.
Let the angle of rotation from xy to .1'//' be <^ (Fig- U-IO). Then
ttii = a22 = cos (f), ai2 = — ttoi = sin 0, and a.u = n-^-i = 0^ = 023 = 0.
Sec. 11-9] Plane rotation of uses: Mohr's-virde representation 269
With respect to the rotateil axes .c'//'. we have. Irmu Imjs. (11-2")) and
(11-20) or the more jj;eMeial Imis. (11-24).
/^„ = Ixs sin- -I- /„„ cos- -I- 2/,„ cos sin (1 1-84)
J't^ = i^xx - I„u) P"!^ <t> ^=ii> + (f'os- - sin- 0)/^„
/' = -^-^
' I'!'
I f_
/xx
" yy
ny =
270 Moments and products of inertia
If we now let the unrotated axes be the principal axes from which the ]
hi + hv I
hi —A
r =
hy-
— /
;incl products of inertia with respect to a pair of axes that are ()l)tained
l>y a phine rotation from the principal axes ^ and t;. By a reverse pro-
(cdure, the same method can also give us the orientation of the principal
;i\esand the principal moments of inertia from a given set of /„, I^y, /,„.
A little deliberation of the geometry of Fig. 11-11 will show that
Eqs. (11-3(3) and (11-37) come out as natural conse(iueMces in the method
of Mohr's circle.
A more formal proof of the Mohr's-cirde construction may be obtained
as follows. From Efjs. (11-35), we have
Squaring both sides of these ccjuations and adding the results, we obtain,
after some simplification,
(J,
/..
I XX ~r
^
' ii!i\'
)
j^
+
jfo
^'^ ^
V ——
I i xz ^ i yu\' I
J 2
This is the equation of a circle in the /^^, /^^ plane. The center of the
circle located on the abscissa a distance of (/„
is + lyy)/'^ from the origin,
and the radius is equal to \[{Ixx — Iyy)/2]- + Ixv^\^'-. It corresponds
exactly to the construction previously described.
Since the accuracy of such a graphical construction is usuall}' poor
unless great care is perhaps lies
exercised, the value of the Mohr's circle
in enabling us to visualize the over-all behavior of the problem and to
the moments and products of inertia for axes x' and y', rotated through
a positive or counterclockwise angle of 30° with the x, y axes at the
point 0. Clearly the plane rotation of axes applies in this case; we
therefore propose to use the Mohr's-circle construction.
We first compute the moments and products
of inertia with respect to
the unrotated axes at 0. This may
be accomplished by the use of the
method of composite areas (Sec. 11-3), the parallel-axes theorems, and
the tabular values given in Appendix A. Suppose we divide the figure
.
Sec. 11-9] Plane rolalion of axes; Mohr's-circlr rrf)resenlali()n 273
defined by
or 20 = 77°12'
and is measnred negatively from the x axis. The minor principal axis 77
4320 I
2376
Fi(i. ii-i:i
PROBLEMS
11-1. Show that III = I [r^ — (r • n)^] dm, where L is any line in space passing
through the origin and defined by the unit vector n, and r isthe position vector of an
element of mass relative to the origin. Verify Eq. (11-25) by use of this relationship.
11-2. direct integration find the mass moment of inertia of a homogeneous cone
By
of radius rand altitude h about the principal axes passing through its center of mass.
11-3. Determine the mass moment and products of inertia of a homogeneous
paraboloid of revolution about center-of-mass axes.
11-4. A composite body consists of a homogeneous hemisphere of radius r and mass
density pi atop a right circular cylinder of radius r and length I and mass density po.
Find the mass moment and product of inertia of the body about axes whose origin lies
at the base of the cylinder. Explain how you would obtain the mass moment and
products of inertia about center-of-mass axes.
11-5. From the result of Sec. 11-2, find the moment of inertia of a solid hemisphere
about a diametral axis A A. Then by means of the parallel-axes theorem, find the
moment of inertia about an axis BB tangent to the North Pole.
Fig. P 11-5
11-6. Write out the complete set of Eqs. (11-22) and (11-23) for a two-dimensional
orthogonal transformation assuming <p is the angle from x[ to Xi. Also write out the
orthogonality relations, Eqs. (11-20) and (11-21).
Fig. P 11-6
11-7. On the basis of Eq. (11-33), write out the cubic equation of X.
Problems 27.')
11-8. A uniform cylinder is of radius /?, lenRth /, and mass //(. Find tlir
solid
moment of inertia /^^and the produet of inertia /^^ about the inclined axes x'y'
shown. (The x' axis n)akes an enclosed angle of 'M° in the vertical plane witli the
base of the cylinder.)
Ill + lyv is invariant. Verify this result by the use of Ecis. (1 \-'.\\) and by the Molir's-
cirde construction.
11-13. If the inertia ellipsoid at the mass center of a homogeneous circular cylinder
is a sphere, find the length of the cylinder as a ratio to the radius.
11-14. Calcidatc the moment of inertia of a homogeneous right circular cone of
altitude h and radius of ba.se a with respect to a generator.
11-15. Derive Ecj. (11-24^).
1 1-16. Using \'a\. (1 1-24«), write out tiie complete exj)re.ssions for I .,., and Z,,. Al.so
Using these two equations, derive once more Kqs. ( 1 l-iiS) for the ca.se of a plane rota-
tion of axes without going through the general Eqs. (11-24) to (11-20).
11-18. For a plane rotation of axes, find the orientation of axes about which the
product of inertia is zero. Also determine the orientation of axes about which the
276 Moments and products of inertia
product of inertia is a maximum or minimum. How are these axes oriented with
respect to the principal axes?
11-19. Verify the results of Prob. 11-18 by the method of the Mohr's-circle con-
struction.
11-20. Find the area moment of inertia and the area product of inertia of the
Z-shaped area shown in the figure.
-h.
Fig. P 11-20
11-21. Determine the moment and product of inertia of the area included between
the X axis, the curve y = (fe/a'')x", and the line x = a with respect to the xy coordi-
nate system.
11-22. Find the area moment of inertia and the area product of inertia with
respect to the centroidal axes of a standard unequal angle section shown. Find the
principal axes and the principal moments of inertia. Draw Mohr's circle to check
the calculations.
-
CHAFIER 12
F = mrc (10-2)
The reader is reminded of the fact that the first e(|uation de.scribes the
translation of the ma.ss center, the second e(iuati()n describes the rotation
about an arbitrary fixed point, and the last equation describes the rota-
tion al)out the mass center. Whereas the first ecjuation is always needed
in finding the translatory part of the general motion of a rigid body,
either the second or the third equation may be used to obtain the rotatory
part of the motion of the body. The choice between these latter two is
of course be pointed out that Eq. (10-11) retains its general form for
moving points other than the mass center provided these moving points
meet the restrictions specified in Sec. \i)-'.i. In general, we shall restri(^t
ourselves to the u.se of K(\. (10-1 1) as applied to the moving ma.ss center.
f
p = .Ti + yj + zh (12-1)
and the angular velocity of the body is denoted by to. The velocity of a
point of the body relative to C is, by Eq. (7-10),
Lc = /^ (9 X e) dm (12-2)
The integrand of Eq. (12-2) can, with the help of Eq. (3-21), be written
(e X e)
= 9 X (w X e)
= (e •
e)<o - (p • w)?
= {x- -|- 7/2 -f- z-)io — (xcc^ + !jw,^ -t- zco,)q
= [{y- -f 2;-)cox — xyuy — xza),]i + — [ //.rco^ -|- (.r- -f z-)oo,i — yzu^,]j
-{- [
— zxw, — zyoiy + (.1-- + /y-)co,]k
Substituting this evaluation into Eq. (12-2) and applying the definitions
of moments and products of inertia [Eqs. (11-1) and (11-2)],
t In this soction and in the remaining portion of lliis cliaptor, for simplicity, wc
shall not \isc the asterisk notation to denote moments and ])roducts of inertia taken
with respe(rt to center-of-mass axes uidess such notation is warranted for clarity. The
context of the equation should indicate clearly which origin wc are using.
Sec. 1 2- 1] hauler's equal ions of molioii of a rigid (tody 279
Comparing tlir lij^lit side of this eciuatioii witli the array of i;(i. (1 l-:i), i.e.
r = xi + //j + zk
where .r. //, and z are the coordinates of any point relative to the fixed
point. The velocity of any point of the body is in this case due to a pure
rotation about the pivot. Hence
r = o> X r
Substituting these expressions int(j the prcxiously defined form of Lo. i.e.,
we find that Lo assumes exactly the .same form as K(\. (12-3). Hence,
tnomentutn of anij rigid body ahoid its ynass center and the
(hr atujular
angular momentum of a pivoted rigid body about its pivot can both he
expressed by Eq. (12-3), if in each case the point abont which the angular
momentum is taken is considered as the origin of the coordinates. In
addition, the angular momentum of any rigid body about .some arbitrary
moving point 0' can also be expre.s.sed by K(\. (12-3), provided 0' fulfills
the restrictions specified in Sec. 10-3. Furthermore, in both cases the
external moment taken about the proper point is etjual to the rate of
change of angular momentum. It is for this reason that we shall hence-
forth u.se the .same treatment for both cases. In particular, we shall
280 Dynamics of rigid bodies
This is the general vector equation of motion of a rigid body written for
axes fixed in the body. Its three scalar components are
Mx = Lx -f COyL,
(12-8)
Ly = lyyC^y (12-9)
Lz = IzzWz
Mr = Lx -{ (JifyLi — wfzLy
My = Ly -{ u3f,L, - w/xL, (12-12)
where Lx, Ly, L^ are given by Eris. (12-0) and are referred to the moving
triad. This method is justified only if we can simplify Eqs. (12-12) by
.selecting axes that are permanently principal axes although they may
not all be fixed in the body.Assuming that this may be achieved, and
it may be, provided that at least two of the three principal moments of
inertia are the same (i.e.. the inertia ellipsoid of the body under considera-
tion is an ellipsoid of revolution;, Eqs. (12-12j may now be rewritten
Once again, the origin of the moving triad may be cither at the fixed
point, the moving mass center, or any other moving point that satisfies
the conditions descril)ed in Sec. 10-3.
Example 12-1. Suppose we solve again Example 10-3, this time i)y
the u.^e of Euler's equations. Let us choo.se a set of principal axes fixed
in the body such that the x axis coincides with the axis of the bar, the
282 Dynamics of rigid bodies
y axis is perpendicular to the bar and lies in the plane of the bar and the
vertical axis of rotation, and the z
Let us treat first a rigid Ixxly in general ami then the special ease of a
pivoted body.
Rigid Body in General. Consider I^i. (12-1.")). The left side, repre-
senting the portion of work done hy the external forces due to the motion
relati\e to the mass center, is
= X
^ /,' ^ • Cp. P.) dt
2
/,"' F, • rfp. = j^ CO • M dl (12-lG)
^ ^2"'tPi' = / .
^
ip- dyn = '
2 / . (w X p) • (to X p) dm
= '-i / [('^)/2 — o),!j)- 4- (coj.r — u}j,z)- + (oixU — oj„.r)-] dm
= M<^-L (12-18)
The last step follows from the expression for L (p]cjs. 12-9). Substitut-
ing E(is. (12-10) and (12-18) into E(|. (12-15), we obtain for a rigid body
f"(.) •
M dt "._„o. iJ^ (12-19)
The total kineti(; energy (denoted li\- 7'), howcxcr, is given by Im|.
(10-25):
as the kinetic energy of rotation about the mass center. This means
that for a rigid body in general we can write the total work-and-energy
equation (12-14) as follows:
^ f^ F, . dn = t\1= ^mvc' \\
+ 3'^co • L
I'
(12-21)
This equation can be interpreted as the sum of the two equations (10-27)
and (12-19), namely,
2 2
F • drc = yimvc"^ (10-27)
/.^
/,
^
(0 • M d^ = Ho> L •
(12-19)
Of these two equations, the first is the work-and-energy equation for the
translatory part of the motion and the second is the Avork-and-energy
equation for the rotatory part of the motion about the mass center.
If the forces acting on the body are all conservative, the total work
'^jy.'dv, = Vr-V,
where V is the potential energy of the rigid body. Hence Eq. (12-21)
becomes, for the special case of conservative forces,
merely the product of the weight of the body and the elevation of the
mass center above some arbitrarily chosen datum.
Pivoted Body. We now consider the rotation of a pivoted body for
which the pivot is taken as the origin of the principal axes xiiz fixed in
the body. Since the origin is now permanently fixed with respect to an
absolute frame, the appropriate equation to be used in this case is
Eq. (12-14). Following the procedure used in the derivation of Eq.
Sec. 12-2] ^\ ork-and-enenjy equalion of (I rigid body 283
= j'^ a> • .M dt
which is seen to have the same form as Ef]. (12-19). The difference
between the two ecjuations (12-19) and (12-24) lies in the diiTerence of
the two origins. It can be seen that for a pivoted body [Eq. (12-24)]
using the pivot as the origin is more convenient, since IVI about the pivot
docs not involve the (usually unknown) external forces at the pivot.
For external forces that are conservative, Eq. (12-24) can be written
Since the equation of work and energy for a particle is a result of the
integration of the eciuations of motion, we would expect a similar result
to be forthcoming in the case of a rigid body. To verify this, let us
multiply each of the three of Euler's equations [E(is. (12-10)] by w^, u>y,
M • (O = -
^^
(/xxCOx- + I,„(^u- + Izz^^z"^) = :j ^ (to • I.)
12
/2 M -tor// = l.^to- L I
: :
The origin of this set of axes is the mass center C. For missiles with
cruciform fins, the y and z axes may be interchangeable. It might be
added that the x axis as defined
is in general not the direction of
flight.
6 = attitude
r/' = heading
<i>
= bank
The angle 6 is defined as the angle between the x axis and its own projec-
tion on the XY plane, and defined as the angle between the projection
\p is
of the X axis on the XF plane and X axis. Thus d and jointly determine rp
the orientation of the x axis. The third angle represents the amount
of roll or spin about the x axis, defined as the angle between the y axis
and the line of intersection (CD in Fig. 12-3) of the yz plane with the
XY plane. It is somewhat unfortunate that these three customarily
used angles in the analysis of air-borne vehicles are not the Eulorian
angles used in dynamics. For this reason we shall limit the usage of the
Sec. 12-3] l'^.r(uii/)l(': niolion of an (lir-hornc iTliicIf 287
symbols d, 4'. and in the inaiiiicr just (ie(iiu'i.l lor the reinainder of this
section only. Elsewhere in this chapter the same symbols represent
the Eulcrian anj^lcs introduced in Sec. 7-3.
The time rate of change of the thred angles 6. \p, and are represented </>
by the angular- velocity vectors 6, and whose directions, shown in ifr, <{>,
Fig. \'2-'-\. are perpendicular to the i)lancs in which each of the three
A' y plane
yz plane
where i, j. k are the unit vectors along .r, //, z. The svuii of 6. t[r. and is
<j>
6 — w„ cos — Wj sin
t In current practice the symbols p, </, r are often used instead of tux, w„, w,.
288 Dynamics of rigid bodies
For the special case of an airplane in level flight, 6 and are both (/>
The gravity force acts through the mass center C of the vehicle; hence
its moment about C is zero.
In addition, we also have the propulsive forces and aerodynamic
forces. The propulsive force, or thrust, is due to the action of a jet of
expelled air or gas behind a propeller, a jet engine, or a rocket (Sec. 10-4).
Although the main function of the propulsive force is to move the vehicle
forward, its moment about the mass center may also be used to adjust
[see Eq. (7-19)]. Thus the three scalar equations of Eq. (12-27) describ-
ing the motion of the mass center are
and their derivatives 5i, 82, ... Thus, in the most general form,
.
with analogous expressions for the other force and moment components.
Assuming such functional relationships are known, Eqs. (12-2(j), (12-28),
and (12-29) can then be used to calculate the orientation of the vehicle
at every instant of its flight. The same equations can also be uised to
predict the general behavior and the stability characteristics (i.e., the
response to a disturbance) of the vehicle. Furthermore, these equations
can be used to calculate the location of the vehicle at every instant.
To do this, the following auxiliary ecjuations are needed in relating
the velocity components Ui, U2, Us along the body axes x, y, and z to the
velocity components Ui, U2, U3 along the fixed axes X, )', and Z. Hy
the geometry of Fig. 12-3, it can be seen that
-|- ?/3(sin 4> i^iii 'A + cos </> sin 6 cos xj/)
('2 = Wi cos d sin \l/ -f ?f2(cos <(> cos ^ -j- sin sin d sin \f/) (12-30)
+ U3(— sin <f> cos ^ -j- cos sin sin ^)
The position vector re of the vehicle at any instant with respect to a fixed
/
reference, say its initial position, can then l;e obtained by integration.
rex
(12-31)
290 Dynamics of rigid bodies
t For a spinning missile, cox may be quite large and the quadratic terms involving
Wx may not be negligible in all six equations of motion.
Sec. 12-4] Slabilily of (he rolalional molion of o r'ujid f)ody 291
where i, j, and k are unit vectors oriented in the directions of the prin-
cipal inertia axes at the mass center. We also assume that ii is of the
form
il = CO + p (12-.S4)
hA - {lyu - hzH^h =
Iyytl„ - (/,, - /„)S2.12. = (12-35)
Piccalling that the initial velocity vectors also satisfied Euler's ecpiations
of the form of E(is. (12-35), E(is. (r2-3r)) then simplify to the following
perturbation e(| nations:
Vv — .
zz
j—f
J-iiy) 2^
^=^Vy = r>(\\
(12-39)
^ VV^ zz
As Case o, we see that if the quantity (/,,;;/ — ^z^ (J^x — L,,) is posi-
in
tive, the motion will be unstable; i.e., px will tend to "grow" without
bound. l^]xamining again the inequality of Eq. (12-32), we see that
(/,„ - /,,) < 0, (/.. - /,,) < 0. Therefore (/,„ - Lz){Ixx - I.J > 0,
Case
c. If the initial motion is assumed of the form oj, = a)„ = 0,
CO, =
constant, and proceeding in a fashion as in Cases a and b, we find
that pi — constant aiul that
It is readily seen from Eq. (TJ-iVJ) that (/„„ — /„)(/j. — /„) is negative.
Therefore px remains hounded and the motion to = oj^k is considered
stable.
We may therefore draw the following important conclusions concerning
rotational motions of a rigid body: Rotations about the axes correspond-
ing to the maximum and minimum principal moments of inertia are
stable motions, whereas rotation about the axis corresponding to the
intermediate principal moment f)f inertia is an unstable motion.
It should be noted that instability as defined here does not mean that
the disturbance would continue growing without limit. While the
exponential solution would tend to indi(;ate this, we must not forget that
the solution was obtained utilizing the assumption that /> 1^. As p «
tends to grow, this assumption would be violated, and the exact solution
should be determined from the iMiler eciuations without recourse to the
simplification. The important point we wish to make, however, is that
rotation about the axes of the intermediate principal moment of inertia
cannot be maintained without external forces (i.e., some form of control)
and will eventually lead to tumbling. The practical im])lications of this
characteristic are clear; for exami)le, any spinning i)rojectile must be
designed so that the spin axis corresponds to the maximum or minimum
principal moments of inertia and not the intermediate one. Numerous
other applications will be encountered in various fields of engineering.
One further remark should be made; the magnitudes of the dilTerences
of the moments of inertia (/„ — lyy), {I,,,, — Iii), {Izz — In) determine
the rate at which the disturbance will grow for unstable motions or the
fretpiency of the disturbance oscillation for stal)le motions. I^ow-
frequency oscillations will usually lead to relatively large excursions
away from the desired initial motion, (lenerally speaking, then, it is
desirable, for stable motions, to distribute the mass .so that the dilTer-
ences of the piiiicipal UKjments of inertia are relatively large.
Ill>J<^!l
="
U — ^z /xx)cO,COj (12-43)
i zz<j>z = \i XX J nil) '^^^ II
t A11lH)Uf>;li wo \v(M-c told that Cardan (1501-1575) was "unal)lo to assign a use
for it."
Sec. 12-6] Biyid lunly under no forces; Poinsol's repreaenialiou 20r>
to the motion of a rigid body which is pivoted at one point other than
the mass center provided there is no external force, or at least no external
(a)
Fir;. 12-4
moment about the pivot. This situation is, lunvever, rather diflicult
to realize in practice. Unless the pivot is the mass center, as in the cases
described in I'ig. 12-4, we would need at least one external force whose
moment about the pivot always counteracts that due to the gravitational
force.
Before we attempt to obtain an analytical solution to the problem
posed at the beginning of this section, we shall consider an elegant
geometrical representation of this problem due to Poinsot. From
Ec|s. (12-44) and (12-45). we .see that, with respect to the rigid body and
the principal axes x\fz, rigidly attached to the body, the angular-velocity
vector <o in a given body with initially specified motion is no longer
arbitrary. If we let the terminal of w be located at the point ha\ing
coordinates ix,ij,z) so that
X y z
But this line is just the line of action of the angular-momentum vector
L corresponding to the given o> [Eqs. (12-9) and (12-46)]. Hence the
statement made at the beginning of this paragraph is proved.
Since L is a constant both in magnitude and in direction when = 0, M
the direction of L can be represented by a fixed line in space, called the
dF_
{x - Xo) + T— I (2/ - ?/o) + -T- (z - zo) =
dx dy |o dz
X — xo y — yo z — zo
dF\ dF\ dF
dx |o dy \o dx
Sec. 12-6] Rigid body under no forces: Poinsol's represenlalion 297
invariable line (Fig. 12-5). At a distance of u> L/L from tlie mass
•
able line; this distance is constant in a given problem [cf. Eqs. (12-44)
and (12-45), from which the distance' is 2T/L] and represents the projec-
tion of o> on L. The plane so drawn is again fixed in space, and is called
the invariable plane.
By
this construction, the terminal of o> must lie on the invariable plane.
But we have just proved that the plane passing through the terminal
of o> and tangent to the Poinsot ellipsoid is normal to L, the invariable
line. Hence the invariable plane is always tangent to the Poinsot
ellipsoid. At the point of tangency, since it lies on the angular velocity
(0, the Poinsot ellipsoid has zero \elocity. The hodij therefore moves in
nvariable plane
Poinsot ellipsoid
Fk;. 12-5
such a manner that its rigidly attached Poinsot ellipsoid rolls, with C as
pivot, on the invariable plane. The angular-velocity vector at any
instant is represented by the vector from C to the point of contact of the
Poinsot ellipsoid and the invariable plane. The locus of the point of
contact on the invariable plane is called the herpolhode. We thus can
make the following statement: The polhode on the Poinsot ellipsoid rolls
without slipping on the herpolhode in the invariable plane. The so-called
space cone is obtained by joining the fixed point C to the herj)olhode,
while the body cone is formed by joining C to the polhode. f The angular-
velocity vector to therefore traces out the space cone in the ab.solute frame
of reference; it traces out at the same time a body cone with respect to
the body. Poinsot's representation also means that the body cone rolls
on the space cone.
Poinsot's representation of the motion of a force-free (or, more appro-
priately, moment-free) body clearly demonstrates the following impor-
t In general, both the body and space cones need not be riglit circuhir cones.
298 Dynamics of rigid bodies
about only one of the principal axes. Then o> and L are constant and
have the same direction and the motion will remain a pure rotation
about the same principal axis as if the axis were fixed. The body cone
and the space cone degenerate, in this special case, into a straight line
along (0 or L. Thus, only if motion is initiated by a rotation about a
principal axis will a steady rotation ensue.
stants,is indeed the motion that satisfies Euler's ecjuations (12-4.*}) for a
moment-free, sj'mmetrical body.
To show this, let us choose, in addition to the z axis already assi}»;netl
along the axis of synunetry of the l)ody. the xij axes such as to form a
right-hand system (Fig. 12-()). This coordinate system is fixed in the
body at C and rotates with it. Clearly this situation fulfills all the
conditions necessary to use Euler's ecjuations, for the .r// plane is per-
pendicular to the axis of .-^ymnu'try and thus the .r and ij axes are per-
Herpolhode
Polhode
Space cone
Fic. 12-6
manently principal axes. Kulcr's eciuations (12-43), for the case under
consideration, now become
I\6}x = — (/ — I \)0iy(jiz
/ici,, - (/ - /i)co,co. (12-49)
/w. =
From the last of Eqs. (12-49) we immediately determine that w, and
= /a;, are both con.stants. In view of this the first two of Eqs.
/.. (12-49)
become
/ - /.
=
ojjo;,/
/l
h- I OJjCJz =
300 Dynamics of rigid bodies
^y + n'-ooy = (12-50)
where w^ = {
—— j
co^
j
Since co^ and coo are constants, we see immediately that u is also con-
stant in magnitude. Moreover, from Eqs. (12-51) and (12-52),
Thus the extremity of o> travels around a circular path with center on the
z axis ofsymmetry, the polhode; the body cone is a right circular cone,
and the angular velocity o> precesses about the axis of symmetry with a
constant angular velocity of precession of magnitude equal to
Id' I oi.
This is the precessional velocity of the symmetry axis about the fixed
line in space repre.sented by L or the fixed axis CZ.
The magnitude of the angular-velocity vector o can be obtained from
lOqs. (12-53) and is eciual to
Since x and z are the principal axes, the angular momentum along these
two axes is
Lx = Iioix = /ico sin {d — a)
Lg = loiz = /co cos {d — a)
But the vector L is along the invariable line or the axis CZ. Hence
Lx =L sin 6
L^ =L cos d
/i
Thus tan d tan (6 — a) (12-57)
I
Body cone
Fig. 12-7
Fig. 12-6, in which the precession ^ and the spin about the axis of sym-
metry, t|, are both counterclockwise for an observer looking from the top
down. This motion is called direct, or progressive, precession.
Case 2. I > lx. This is the case of a flattened gyroscope such as
the one shown in Fig. 12-4a. Here 6 < d —
a, and thus to has to be on
the other side of L (Fig. 12-7). By Eq. (12-55) we find that ^ becomes
negative if 4> '•'^ positive. That t[r should be negative is also obvious
from a geometrical consideration since if and must combine vectorially t|r
to form 0). Viewed from the top down, the precession 4> is counterclock-
wise but the spin about the axis of symmetry is clockwise. This
t[r
See. 12-8] Symmelricul (jyroscope under exlernul moments 303
Fk;. 12-8
specified, applies to either of the two cases: (1) the gyroscope has one
point fixed (pivoted) in .space, and the origin is at that point; or (2) the
gyroscope has no fixed point, and the origin is the ma.ss center.
To describe the general motion of a gyroscope, let us assign a set of
ab.solute coordinates XYZ and a set of moving coordinates xijz (Fig. 12-8).
We let the z axis coincide with the axis of symmetry of the gyro.scope.
Furthermore we select the x axis to lie in the meridional plane ZOz,
perpendicular to the z axis. Thus the y axis is horizontal and lies in
We first compute the angular velocity of the movmg triad. With our
choice of axes, the third Eulerian angle \p (Sec. 7-3) is zero, and Eqs.
(7-11) to (7-13) become
co/x = —4> sin 6
o^fy = e (12-58)
CO/2 = ^ cos d
The total angular velocity of the gyroscope differs from this only by an
amount equal to the spin of the gyroscope about its own axis 4^. Thus
cox = — sin <^ 6
w, = d (12-59)
w^ = cos 6
4> -\- xj/
I
Substituting these values of angular speeds into Eqs. (12-13), we obtain
the following three scalar equations of motion
M, = I
j^
{cj> cose + rP)
The same set of equations of motion can also be derived if, instead of
adopting the moving frames of reference just described, we adopt a set of
body axes xyz which are fixed in the body and moving with it. As
before, the z axis is the axis of symmetry of the body. .4^ the instant
under consideration, we let the x axis lie in the meridional plane ZOz, per-
pendicular to the z axis. This appears to be the same as before, except
that here the x axis assumes such a position only for an instant. The
y axis is again horizontal and lies in the plane at the same instant. XY
With such a position of the body axes, = at the instant under con- i/'
sideration, and Eqs. (7-11) to (7-13) yield expressions for Wx, c^y, coz which
are identical to Eqs. (12-59). In order to apply Euler's equations [Eqs.
(12-10)], we also have to evaluate Wx, 03 y, and co,. They are obtained by
differentiating Eqs. (7-11) to (7-13) with respect to time before setting
1/' = 0. (Since xp = only instantaneously yp is not necessarily zero.)
Thus
u^ = d\p — — {(!> sii^ ^)
Now substituting Wx, co,,, w.-, Wx, co,;, and co, into Eqs. (12-10), the same set
of equations, Ecjs. obtained again. These are the general
(12-59), is
precession we mean that the velocity of precession 4> and the half-cone
angle ^ are bothcon.stant. With such a motion, Eqs. (r2-()0) are reduced to
il/. =
M^ = [I\P + (/ - Ii) (j) cos e\ 4> sin (iLMUj
il/. =
Thus a steady precession is possible if an external moment .1/,, ol the
magnitude given by the second equation above is supplied. If .1/,, = 0,
the .same etjuation shows a relation
between ^, <^, and 6. This relation
is just Eq. (12-55) of Sec. 12-7.
M = My = 1x1^4, (12-02)
to the shaft, the gyroscope will rise or fall to a new 6. This characteristic
of a gyroscope is the basis of many engineering devices, such as the gyro-
scopic compass, the gyroscopic stabilizer for ships and torpedoes, the
artificial horizon for airplanes, etc. It is also the cause of undesirable
shaft and bearing stresses in ship turbines, airplane jet turbines, airplane
propellers, etc., when the ship or the airplane rolls, turns, or pitches
(Probs. 12-6 and 12-7).
In the case of very large spin velocity i/- compared with <j) and 6, Eqs.
(12-61) can be approximated by
general along three different directions. To show that Eq. (12-63) can
sometimes be in gross error, we merely have to recall that in the free
motion of a gyroscope, 0, and d are all nonzero but the moment is zero,
i/-,
whereas Eq. (12-63) would erroneously show a nonzero My. For d = 90°,
Eq. (12-63) becomes the exact equation (12-62).
Special Case 2. Steady Precession Due to Gravity of a Top Pivoted at
One End (Fig. 12-10). Let the mass center C be located at distance a
along Oz from 0. The moment about due to the weight of the top is
M = ak X — mgK.) = mga
( sin 6 j
Comparing this equation with Eqs. (12-61), we find that a steady preces-
sion is possible under gravity alone if
For a given top, this equation relates the three quantities i/', (j), and d.
We see that the gyroscopic moment is in this case balanced by the moment
due to the weight of the top. Furthermore, we see from Eq. (12-64) that
the larger the spin ;/-, the smaller the precession For very large i/, and
<i>.
hence very small 4>, Eq. (12-64) can be approximated by Eq. (12-65)
even for 6 different from 90°.
Sec. 12-0] The gyrocompass 307
Fio. 12-10
consider the spinning projectile, the rolling hoop, the bicycle, etc. Many
of these phenomena may be deduced from the general efjuations of
motion (12-()0) used in conjunction with the work-energy principles.
the z axis and lies in the horizontal plane, and the y axis is perpendicular
Tangent (horizontal)
plane
Fig. 12-11
Since
The anguhir velocity of the gyroscope differs from this only because of
the spin about its own axis z, so that
xj/
CO = — oje cos X sin & i + (cog sin X + ^)j + (w^ cos X cos d -\- 4')k (12-07)
The only external moment due to the constraint supplied by the bearings
is about the x axis; thus M = Mi. Substituting these values of co/, co.
Assuming now that the spin is large compared with w,. and neglecting i/'
terms of the order of co^- when compared with terms of the order of co.i/',
This shows that for ip coc, the spin » ip is approximately constant. Equa-
tion (12-70) can now be rewritten
6 + /(•-
sin 6 = (12-71)
where n^ = — , — = constant
Mc = (12-72)
We can also write these two vector equations (10-2) and (12-72) in terms
of their cartesian components:
F^ = mxc Fy = mijc F^ = mzc n9 "q^
^^^"'"^^
Mc. = Mc, = Mc^ =
These six scalar equations completely determine the motion of the mass
center of the body and the external forces required for the prevention
of rotation of the body.
Since a rigid body in translation does not possess a pivot point, Eq.
(12-72) is valid in general only with respect to the mass center. Thus L
and M are not zero with respect to any fixed point. They are zero only with
respect to the mass center of the body.
It is of some interest (mainly academic) to examine the form of the
moment equations when they are taken with respect to points other than
the mass center. For some fixed point 0, and for some arbitrary moving
point 0', these equations take the following form:
Mo = j^
r X r dm (10-8)
Here we have iutroduced the iiitej^ration over the entire vohiine ot the
bodj' for a continuous mass cHstriWutiou.
Equation (10-8) becomes, upon performing the iniUcated (hlTcrentia-
tion. using Va\. (o-O) and recognizing that in translation the acceleration
of every point on the l)t)dy is the same.
CO:, = OJ,, =
But the z axis thus assigned may not be the j.^,, y^.vi
principal axis. We therefore cannot use the
Euler's ecjuations (12-10). In.stead, we have to start with the more
general treatment [Eqs. (12-6)], which becomes
Lz = — hzo:
L,j = —lyzO}
Lz = /«w
Substituting these components of the angular momentum into Ecjs. (12-8),
The third equation completely determines the rotation. The first two
equations determine the external moment required to keep the axis (i.e.,
the shaft) fixed. These external moments are usually supplied by bear-
ing reactions, i.e., forces exerted by the bearings to the shaft.
In the elementary treatment of rigid-body dynamics, it is often assumed
that the body is symmetrical with respect to either the xy plane or the
axis of rotation or both. When the body is so restricted,
I- yz -t xz O
(Sec. 11-8) and the angular-momentum vector L coincides with the direc-
tion of fa). Hence
M, =
Aly = (12-77)
rotation, since the products of inertia also vanish for such a case if the
origin chosen to lie in the plane of symmetry.
is
F = mrc
is, of course, also valid. This equation becomes, for the case of rotation
about a fixed axis,
find the four bearing reactions, i.e., the reactions along x and y for each of
the two bearings.
In some cases, depending on the body, it may prove to be advantageous
to select a fixed coordinate system with origin on the axis of rotation and
oriented in the directions of the principal axes. Under these conditions
Euler's equations may be used directly.
When the mass center of a rotor on the axis of rotation, the rotor
falls
is said to be statically balanced. For a statically balanced rotor, the
Sec. 12-11] Rolalion ahoul a fixed art's: dynamic halancing 3i:i
The work (if any) done by the gravitational forces belongs to the left
side of Ef). (12-79); this portion of the work is equal to mg{yc\ — ya)
[Eq. (12-23)].
Example 12-2. Given a thin rectangular plate rotating about a
diagonal axis with angular speed of w and angular acceleration of co.
_ 2m a^b^
L. " ~ T a2 -j-b^
, m , a"^ — b"^
3 a/ -j- b-
ly^ =
Substituting the.se values into Ecjs. (12-70), we immediately determine
the following recjuired couples:
t We are concerned here with dynamic loads only. The .static load (due to the
weight of the rotor, etc.) on the bearings is obtained by statics (Chap. 4) and is to be
added to the dynamic loads.
314 Dynamics of rigid bodies
6 a'- -\- h~
'\n a'
My = co^-ab
3 a'
^
+ , ,
0- .,
,, 2w . a-6-
3 a^ + O'^
Fig. 12-13
CO
(ai' -F 6k')
Va' + b^
and
^^, ^ m wab~
3 Va' + t'
?w a6(a- — 6^)
il/.' = - ^''" 2
fl^ + 6-^
a'b
M' =
3 V«' + ^'^
Fig. 12-14
These four equations enable us to solve for the four unknowns niAXA,
inAVA, tnuXB, and mBi/B. It can be seen that with the correction weights
and their locations so determined, the total assembly consisting of the
rotor and the correction weights is dynamically balanced. Furthermore,
since only the four products niAXA, etc., are significant, we can even
arbitrarily assign values for rriA and m^ or, alternatively, their distances
from the axis of the rotating shaft (Prob. 12-21).
Mx = Iyz<^^ — Jxz<J^
Fx = mxc
Fy = myc (12-81)
F, =
We find that the equations of motion (12-80) and (12-81) are of the
same form as Eqs. (12-76) and (12-78). The only difference is that the
moment equations (12-80), for plane motion, are taken with respect to
the mass center as origin, whereas the moment equations (12-76), for
rotation about a fixed axis, are taken with respect to a point on the axis
of rotation as origin.
The procedure in applying Eqs. (12-80) and (12-81) to solve a problem
of plane motion is not as clear-cut as the procedure for solving problems
of rotation about a fixed axis presented in the previous section. This is
because the unknown reactions do often contribute moment about the
.
mass center, and therefore the third of Ecjs. (12-80) alone is often not
sntlicient to solve for w. In general, we have to solve the simultaneous
efiuations consisting of all of Eqs. (12-80) and the first two of Ecjs. (12-81 ).
If there a kinematic relation in the problem relating .re, j/c, oiz< and w.-.
is
Fx = rnxc
Fy = myc (12-82)
The equation of work and energy for a body in plane motion is often
also of interest. The total kinetic energy of the body is, by K(\. (12-20),
This equation can be used in connection with Ya\. (12-21) or, for conserva-
tive external forces, with PJq. (12-22).
Example 12-4. A nonhomogeneous
circular cylinder with its mass center C at
a distance a from its geometrical center .1
(6) The forces acting on the body consist of the weight force mg
acting at C and of the normal force N
and the friction force F acting at B
(Fig. 12-16). This friction force is not
a dissipative type of friction force; i.e.,
it is not a sliding friction force, for point
w = ma cos
g + CO-/?
In
are given by the first two of Eqs. (12-84). The ratio F/N must be less
than the coefficient of static friction if our ori<j;inal as.sumption of rolling
without slipi)ing is to hold. If our original as.sumption proves to be
invalid, then a new solution must be obtained for a sliding cylinder. In
this event the friction force is a kinetic friction force; it has direction
opposing motion and magnitude ecjual to aia-.V, where ma- i^^ thecoeflicient
of kinetic friction; and the kinematic relationships pertaining to rolling
without slipping are no longer applicable.
PKOBLEMS
12-1. A
square plate is fixed at its center. At a certain in.stant, the anmilar-
velocity vector makes an angle of 30° with the normal to tiie plate. Find tiie angle
between the angular-momentum vector and the normal.
12-2. A circular disk is rigidly attached to a rotating shaft. The axis of .symmetry
of the disk makes an angle of '.M)° with the axis of the shaft. The disk is 1 ft in diame-
ter and ;J in. in thickness and weighs 100 lb. The shaft rotates at a uniform speed of
1,200 rpm in the direction shown.
(a) Find the orientation and magnitude of the angular-momentum vector.
(6) From the rate of change of the angular-momentum vector, find the moment
supplied bj' bearing react ioii.s.
(c) Find the kinetic energj- of the disk.
Fk;. V \1-1
l2-.'{. .\ body in the shape of a s(iu:irc prism of sides n, a, 'In is j)i voted at its center.
.\t a certain instant the body .spins with angular velocity u iibout ;i diagonal. There
are no external forces acting on the body.
(a) Find the time refpiired for the longitudinal axis of .symmetry of the body to
make one complete movement arfiund the surface of a cone.
(/>) Find the half-<'one angle.
(c) Find the kinetic energy of the body if its ma.ss is m.
320 Dynamics of rigid bodies
Fig. P 12-3
12-4. Express the kinetic energy of a symmetrical body rotating about a fixed
point in terms of the Eulerian angular velocities.
12-5. A homogeneous body having the shape of a solid of revolution moves about
a point. Its motion is represented by the uniform roUing of a cone of half-cone angle a
fixed in the body upon an equal cone fixed in space, the axis of the former being the axis
of the body. Show that the couple necessary to maintain the motion is of magnitude
where w is the resultant angular speed, and that the couple lies in the plane of the
axes of the cones.
12-6. A 2,000-lb auxiliary Laval turbine rotor mounted on a destroyer rotates at
10,000 rpm. The rotor withblades can be assumed, for the purpose of calculating
its
Fig. P 12-6
destroyer rolls inrough weather with an amplitude of 40° from the vertical, and it
takes a period of 10 sec to complete one cycle of roll. The distance between the two
bearings of the rotor shaft is 1 ft. Find the maximum reactions at the bearings.
12-7. The propeller of an airplane rotates at 1,200 rpm. It weighs 600 lb and has
a radius of gyration of 3 ft. If the airplane is traveling at 400 mph in a horizontal
circle of radius 3,000 ft, find the gyroscopic moment on the propeller shaft bearings.
12-8. Prove that a top cannot exhibit steady precession with its axis of symmetry
inclined at 6 with the vertical unless it has a sufliciently large angular velocity about
its axis of symmetry such that
12-9. A unifonn .s(iuaro plate rotates about its vertical diagonal IH) witli the
angular velocity w. The corner A is then suddenly fixed so that the plate now rotates
about O'O'. What is the new angular velocity?
as.sumed that the door is of \miform material and that there is no friction. Find the
value of F for which one of the wheels lifts off. (Which wheel?)
(^
HDk I I B
12-15. two uniform bars weighing 322 lb and 64.4 lb each. They
AB and CD are
are connected C and are hinged to the carriage at A and D. The weight
by a pin at
E of 32.2 lb is attached by a string to the bar AB at B. If the entire system moves to
the left with an acceleration of magnitude 20 ft/sec^, find the horizontal and vertical
components of the forces acting on each bar Sit A, C, and D.
12-16. A uniform stick is pivoted at one end and supported at a fairly small angle a.
To its end are attached two very Hght cups, the outer one of which contains a marble.
When the bar is allowed to fall to the horizontal position on a felt pad, it is found that
the marble is no longer in the outer cup but has jumped to the inner one. Explain
this experiment.
Fig. P 12-16
12-17.A body hinged about a horizontal axis is allowed to oscillate. The center
of mass C is at distance a from 0. If W
is the weight of the body and /o is the moment
of inertia of the body about 0, find the frequency of small oscillations. Find the
length of an equivalent simple pendulum of the same mass that has the same frequency
of small vibrations.
12-18. A connecting rod AB is suspended on a horizontal knife edge at A and made
to oscillate. The frequency found to be 60 cycles /min. Then the rod is suspended
is
vation of energy, find the subsequent angular speed as a function of 0. (b) What is
Problems 323
the amplitude of the motion? (r) Under wliat conditions will tlu- |>endulnm not
oscillate?
12-20. A 2- by 4-ft rectangular steel plate of 2 in. thickness is mounted in hearings
to rotate about a diagonal as shown. The weight of 2-in. -thick steel plates is SO psf.
The plate rotates at 1,200 rpm, and the distance between bearings is o ft. Calculate
the bearing reactions.
12-21. A horizontal shaft carries three unbalanced weights H'l = lb, Wi = 1 2 lb,
and H'3 = 3 lb. The distances of these weights from the shaft are ri = 8 in., r^ =
4 in., /-3 = 8 in. At each end of the shaft is a circular steel disk (,1 and li) of 1-in.
thickness. The entire as.sembly is by drilling one hole at a radius of
to be balanced
8 in. in each of the two disks. P'ind thediameters d of the two holes and the angles of
inclination of their radii. (The specific weight of steel is 0.283 lb/in.')
Fi(i. V 12-21
12-22. In order to retrieve a sjKml of thread by ])ulling at the thn>ad, what is the
limit on the angle a? Assume that there is sufficient friction between the floor and
the spool for rolling without slipping to occur.
324 Dynamics of rigid bodies
W
A homogeneous circular cylinder of weight has a thread wrapped around
12-23.
itsmiddle. One end of the thread is fixed at A. The cylinder is allowed to fall.
Find the tension in the thread and the acceleration of the center of the cylinder.
Fig. P 12-26
of the sphere. (6) At what time does rolling witliout .slij)ping occur? (c) What is
tion of motion in terms of the angle <j> that the center of the wluu'l makes witli the
vertical through the center of the circular track.
12-29. A bar AH of weight U'l lias at one end li a cylindrical roller, of weight \V>
and radius r, that rolls without slipping along a circular track of radius /{. Har and
roller are made of uniform niaterial. Initially AH makes an angle ^n with the verti-
cal. Assume no friction at the hinges .1 and H.
(a) Find the angular velocity of the har .1^ when it falls to the vertical position.
(b) Find the angidar acceleration of the har at the inilial posilion.
(c) What is the frequency of o.scill.itiiin if the initial angle <)>» is small?
Lagrange's Equations
13-1. Introduction
t Our consideration in this chapter will be limited to the so-called "holonomic con-
straints" for which this statement is true. Holonomic constraints are those con-
straints which impose certain definite functional relationships between the coordinates
of the particles in the system. (Such relationships are referred to as integrable
constraint conditions.) For a system with holonomic constraints the generalized
coordinates necessary for the complete specification of the configuration of the system
can be varied individually and arbitrarily without violating the constraints.
326
Sec. 13-1] Inlrodudion 327
For a system consisting of .V particles, the numlxT of generalized
coordinates re(iuired to desc-rihe all possible configurations of the system
would be 'SN none of the particles were constrained. (Juite often, and
if
F^ = mx Fy = my (8-2)
Equation (13-1) states that the total virtual work due to both active
and reactive forces is zero if the particle is in etiuilihrium. Now, if we
restrict ourselves to systems in which friction is not involved, then the
total reactive force at the constraint is perpendicular to the virtual dis-
placement. (For example, is constrained by its contact
if the particle
with a smooth surface, the reactive force at the contact is perpendicular
to the surface, whereas the virtual displacement satisfying the constraint
must be tangent to the surface.) Hence the reactive forces contribute
nothing to the virtual work. (Such constraints are known as workless
constraints.) Thus the principle of virtual work may be expressed as
follows: // a particle with workless constraints is in equilibrium, the total
virtualwork due to the applied forces must vanish under any arbitrary dis-
placement satisfying the constraints. The converse of this statement is
also true.
A more complete statement of the principle of virtual work con-
still
F = -^
dx
Fy =
so that F= — ^-«
d.C
— -T~J
Oij
',
dz
-
5IK = F • 6r = - ^
dx
dx - ^
dy
5ij - "^
dz
bz = (13-2)
|L = o
dx
^= dy
^' =
dz
(13-3)
330 Lagrange's equations
conservative force system if, and only if, the variation of the ^potential energy
about the equilibrium configuration vanishes. If the constraints imposed
y = V(qi,q2, . . . ,qn)
and bV = ^-
dqi
bqi
^
+ ^bq2+
dqi ^
• • •
-f- t-
dqn
5g„ =
f^f^...^f = Q
dqi dq2 dqn
(13-4)
t This conclusion could luivc also been reached by observing from Sec. 8-7 that
SV = —8W and tliat tlu'rcfore, for ('(luilibriuui. 8V = 0.
/
\<l.v)/ '^\d.v^),2\
m. "
I'pon neglecting terms that are of a higher order than dx'-, we deduce that
the change in potential energy is
/d-V\ 8x
AV = r - r., . . „.
\dx'-)o .,,
2!
Equilibrium; -r— = t— =
dx ay
Stable equihbnum; (
- ~- I — ^o ^r^ < ^
\dx dij dx^ dy-
dx- dy-
Unstable equilibrium :
f . . ) — ^V ^-r < ^'
dx' dy^
dx^ dy^
332 Lagrange's equations
Fig. 13-1
TV = -117(1 - cos 9)
The change of length of the spring at this position is rd; therefore the
potential energy of the spring is (Sec. 8-7)
^=
ad
= h-d - 117 sin e
and equilibrium can take place for angles 6, satisfying the transcendental
equation
e = T^ sm d
Icr-
S«'c. 13-3] Lagrange s equation for a jHirCirle 333
, 2
= /»'*- - ^yi t'os e
w= ''' - "''
Let Xi, .To, Xz be the rectangular coordinates of a particle and 71, 70,
where the index i may be any one of the three nunibeis I, 2, 3. 'i'he
scalarcomponents of the \elocities of the particle can therefore be written
in terms of the generalized cooidinates by differentiating the above
Thus any f, is a function of the \ariables dXi/d(/j and (jj for all values
of j from 1 to n. Considering the variation of i-, with respect to a par-
ticular ('/j while kccpitu/ (dl other ('//s and all ('ix,/(l{//s eonstont. we have,
t We restrict our {liscu.s.sion here to tho.sc .systems where flic nLitimis hefween tlie
^=P
dqj dqj
(13-9)
dqj
/ ^^i^Z
L^ dqj
(13-11)
i=l
3 3
dT
and -r-^ / mxi~ = / mxi^ (13-12)
dqj Zv dqj Zv dqj
i=l 1=1
where we have used Eq. (13-9). Since dT/dXi = mxi is the momentum
by analogy we
of the particle, call dT /dqj the generalized momentum.
Differentiating the preceding expression with respect to time, we obtain
3
dt
1=1
Here the index k is simply a dummy index used in connection with the
operation of summing from 1 to n. Furthermore, on account of Eq.
(13-8), this becomes
d /dx\ _ dXi
dt\dqjj dqj
Hence the last term on the right side of Eq. (13-13) reduces to
3 3
L( ^^' dt \dqjj
~ /^ ^-^
' ^q. Qq.
1=1 1=1
where we have used Eq. (13-1 1). Substituting this result into Eq. (13-13)
and taking into consideration that mXi = Fi by Newton's law,
3
d'^
dt \dqj/ dqj L^ dqj
Sec. 13-3] Lagrange's equation for a parlirle 335
The work done by the forces/-', for a siiuill (lisplaceineiit of (lie particle
is rSec. 3-4)
dW = y I\dxi (13-11)
dW
dqj
1=1
and call it Qj. Thus
dW
Qj =
^ (13-16)
F =^-E
'
dXi
which is obtained from Ya\. (13-14), we find the following analogy. The
work done upon a small displacement rf.r, is b\ dx„ while the work done
upon a small displacement dqj is Qj dqj. We therefore call f^/s the
generalized forces corresponding to the generalized coordinates qj. The
(luantity Qj may he thought of as the work-absorbing force along q,.
Qj = -^ dqj
(13-17)
336 Lagrange s equations
L = T -V (13-18)
qj may
be any one of the generalized coordinates in the problem, so that
we can write one Lagrangian equation for each of the g's. We therefore
always have available the same number of equations as the number of
unknowns, which are "picked beforehand to describe most conveniently the
configuration of the system. We thereby automatically obtain a neat
formulation of the problem and in many cases introduce considerable
simplification. It must be realized, however, that Lagrange's equation
helps us only to formulate the problem, i.e., it supplies us with the
differential equations, but that it does not tell us how to solve them.
Another important advantage of Lagrange's equations, so far as the
motion of the system is concerned, is that the forces exerted by the con-
straints do not appear in the equations since they do not contribute to
the work done, W. These forces are generally unknown, and their
presence in Newton's laws sometimes requires additional manipulations
(see, for instance, part b of Example
In addition, the kinematic
12-4).
relations necessary to deal with acceleration are not required;we need
only concern ourselves with the velocities involved in the formulation
of the kinetic energy.
We shall now illustrate the use of Lagrange's equations by two exam-
ples. The first example is completely trivial. The second one shows
the neatness of this new method.
Example 13-2. A particle is constrained to travel along a straight
line which may be considered as along the direction of x. The particle
has one degree of freedom; its position is specified by the single general-
ized coordinate q = x, measured from some fixed point along the line.
Therefore
T = — 2
r
dT dT .
(I /dT\
^ = 5^ =
dq dx
•
\, ^— = ^^ = wr
.
rn t^ = ""'
dqi dr (ft \dqj
( I
dT = dT ...
= m{rd' , ...
+,
np' sni^
.
, „,
^)
dq, dr
m -r {r-4> sin- 6) = Q^
t The small general displacement dr may Ix- fouiid liy forniiiiK tlic (liffcrciitial of
Fe
Sec. 13-1] Lagrange's eqnalion for o sysleni of /xirlirles 330
Also, the potential energy, using the initial position of the cylinder as a
datum, becomes
V = mgijc
^^ ^ 1 + hjmr-
Exaniple 13-5. For a symmetrical gyroscope pi\<)ted at a point
(_l"ig. 12-8), the kinetic energy is (Proh. 12-4)
dd
^ = /(^ -f 4, eo.s 0)
d\V = -.U, sin f/0 + -U„ r/9 -|- il/,(cos 6 d<i> -\- d^)
Hence Qe = -^ = Alu
dtp
With the above values for f/s. Lagrange's efpiations in this case are
identical with Kqs. (12-00) as derived in Sec. 12-8.
340 Lagrange's equations
In the preceding equation, the first term on the right, T^o, can be arbi-
trarily chosen to be zero if we agree to measure the potential energy of
the system from the neutral position. Hence
Fo =
FromEq. (13-17) ^=
oQi
-Qi
At qi = q^ = = r/„ =
the forces acting on the system are in
0,
equilibrium, and hence the component of the total force along any
Sec. 13-.')] Small frci' vihralions of conpled sy.slenis 311
= (13-22)
» >2
2 2^.^.
t=i y=i
bn =
i = \
where two subscripts are now necessary for ./, the first subscript k specify-
ing the particle under consideration. Substituting for its value as ./•
1 = 1 t = 1 j = \
Summing over all the .V particles of the system, we can write tiie kinetic
energy of the system in the following form:
II n
^ = ^^
I I "'>'?'''> (13-23)
.=1 >=i
;=1 1=1
A^ 3
L^ L^ \dqi/o\dqj/c
k=l 1=1
Since T is inherently a positive quantity, it also is a positive quadratic
form.
We see by Eqs. (13-22) and (13-23) that V and T are quadratic forms of
the generalized coordinates q and generalized velocities q, respectively.
Substituting these two equations into Lagrange's equation (13-19),
n n
52 = A2^\n{pl-V <i>)
^ ^ ^
Sec. 13-6] Small vlhrdlitins of lu'o-degree-of-freedoni syslenis 313
aiui A
can have nonzero
. (iiontrivial) solutions only it the determinant
vanishes, i.e.. onlv if
can he shown that roots of the frequency equation are all positive and
real provided the potential energy T is a positive real function. f We
shall denote the two positive roots as pi^, p-^. These two natural
circular frecjuencies are sometimes referred to as the normal frequencies.
In this discussion we
omit from further consideration the special
shall
cases of zero roots or repeated roots of the frefjuency eciuation. Zero
roots occur when it is possible to have displacements without changing
the potential energy of the system, and as such correspond physically
to a rigid-body motion of the system. Repeated roots indicate that the
system vibrates harmonically in two modes with the .same frequency but
with different amplitude ratios. Further di.scu.ssion will clarily this
point.
Corresponding to each of the normal frequencies pi and p^, we have a
solution of the type (take pi, for example)
We notice that the coefficients A,j appearing in Eqs. (l.V2()) are not
independent; they must .satisfy Iv|s. (1:^27) and are thereby related
by virtue of tho.se expressions. For example, substituting pr into the
t Since we have shown that T is a positive quadratic form, we indeed obtain posi-
tive real roots, the fn'(|uencies, from the characteri.stic c(|Uiiti()ii. W'c (h), however,
omit the proof of the above statement.
344 Lagrange s equations
O12 ai2Pi-
t Either of the two Eqs. (13-27) would he satisfactory, since only n — 1 of these
equations are independent by virtue of the frequency equations (13-28). It is easily
verified that the second of Eqs. (13-27) yields identical results.
} For purposes of better understanding of the theory, when solving problems the
student should follow the alternative approach ratlier than the one employed here.
Sec. 13-6) Sntall vlhrulions of I iro-deyrec-of-freedom syslenis 3 15
neutral positions. The.se coordinates can then l)e cho.sen as the general-
ized coordinates of the .system.
For the entire system, tiie kinetic and potential energies are
The four constants An, A12, 0i, 02 are determined by the initial con-
ditions, i.e., the values of Xi, ±1, X2, ±2 at t = 0. Although each mode of
vibration is sinusoidal with time, the superposition of the two modes is
not sinusoidal. In fact, it is not even periodic unless the ratio of the
two natural frequencies is a rational number. If the two natural fre-
quencies differ by only a slight amount, the combined motion will
exhibit the phenomenon of beats.
It should be noted that our solution of X21 = +1, X22 = — 1 in the
above example is due to the fact that we have chosen equal mass for the
two particles and the same stiffness ki for the top and the bottom springs.
For a more general two-degree-of-freedom system, the procedure out-
lined at the beginning of this section should be followed to determine
pi, P2 and the corresponding values of the two amplitude ratios. The
two amplitude ratios will not in general have values of +1 and —1.
n n
Here the constant C is introduced as a scale factor for rea.sons that will
he clarified .shortly. Suhstitutinji; E(|. (13-34) into Va\. (13-33). we obtain
the following /( algebraic eciuations for the n unknown amplitudes .{,:
^iH - ;/-V/i„
631 — p-'asi • •
•
= (13-3G)
h„\ — p-a„i • • • • • •
h„n
f>„„ — p-a„„
ha stated previou.sly, the roots are all real and positive as a roiise(nienr-e of posi-
t
tive values for V and T and a,/ = a,,, fc,, = /'/,.
1
All this discussion must sound familiar to the reader, because it is more
or less a repetition of precisely the same concepts already described in
Sec. 13-6 for a two-degree-of -freedom system. In the following, how-
ever, we shall bring forth yet another concept, the concept of normal
coordinates, which proves to be essential for a deeper understanding of
the problem.
Let us consider any two of the n normal modes, corresponding to two
normal frequencies pk and pi, respectively. The amplitudes Ajk and Aji
for these two modes satisfy the following equations, which are merely
Eq. (13-35):
n ft
fori = 1,2, . . . ,n. Now multiplying Eq. (13-38) by .4^7 and summing |
over all i, !
n n n n I
n n n n i
Y y hijAjiAik = X A praijAjiAik
i=l j = 1 i= 1 j =1
Since i and j are dummy indices, they can be interchanged in the last
expression, which becomes
n n n n
Remembering that a,^ = a^i and 6,j = ^ji (Sec. 13-5), we see that the left
sides of Eqs. (13-40) and (13-41) are identical. Taking the difference
of these two equations,
n n
- pr) a„.l,...l,, =
(Pa-
^ 2=
» = 1 i 1
Sec. 13-7] n-degree-of-freedoin syslenis; normal courdinules 349
a.jAaAjk = k (13-42)
2 ^ 9^ I
n n
,=1 y=i
is, from E{[S. (13-32) and (13-34), proportional to the magnitude of the
kinetic energy due to the mode pk. Since the magnitude of kinetic
energy cannot be zero, the last expression does not vanish. Now it has
been mentioned several times before that Eq. (13-35) can only determine
the ratios of Ajk and not their absolute values. Hence there will be no
loss of generality if we set
n n
a^jA^kAjk = (13-43)
^ 2
1=1 >=i
1
When the amplitudes are made to satisfy Eq. (13-43), they are said
to be normalized. It is now clear why we need a scale factor C in Ya\.
(13-34) or Ck in Etj. (13-37) (the subscript k is to identify C with the
particular p^-^). We need Ck to fix the magnitude of each normal mode
to fulfill the given initial if the amplitudes Aik are normalized.
conditions
Equations (13-42) and (13-43) can now be coniljined to form
ri )i
t = 1 j = l
where 8ki is the Kronecker delta. Equation (13-44) represents the rela-
tions that must exist between the n- normalized amplitudes, referred to
as orthogonality conditions.
The orthogonality conditions can be used to advantage in determining
the magnitude Ck and the phase angles </>*.- of the various normal modes
from the given 2n initial conditions. At / = 0, we have, from Eq.
(13-37),
n
9.(0) = Y,
^*'-^'*^^" '^'^ (l-*^--*^)
k=\
11
Assuming ^,(0) and r/,(0) are given, we can then multiply each side of
the above ecjuations by ttijAji and sum over both i and^. On account of
:
n n
Ci ^^^ ^i^ -
1 Vy V [qi{0)]aijAji (13-48)
/^
1=1 y=i
q, = 2 Ai,q', (13-49)
n n
lk= ^^
i^i y^i
a-u-AMi (13-50)
In terms of the new coordinates qi, the potential energy of Eq. (13-31)
becomes
n n
V = }i J^
X ^''^'^'
n n n n
^ ^^ 2 Z X 2 ^'i^'^^i'^k^'i
i=li=\ k=\ 1=1
n n n n
= H X X ( 2 2 ^^i^i^^i) ^'^<i'i
Substituting into this expression the following identity from Eq. (13-35)
n n
y hijAji = ) pi-aijAii
y=i y=i
9. = A,,q[
Yi
k=\
Thus we reach the remarkable result that, in terms of the new coordi-
nates ql-.
the expressions for potential enerjj;y and
kinetic ener}z;y are both
simplified to such an extent that all cross-product terms disappear.
These new coordinates which render this simplification possible are
known as normal coordinates of the system. The procedure involved in
such a simplification, i.e., in obtaining Fa\s. (13-o1) and (l3-.')2). is called
the reduction, or diagonalization, of quadratic forms.
In terms of normal coordinates, the Lagrange ecjuation
di \d<rj a,;
9', + Pr^A =
92 + P-r(/2 = (13-53)
7: + P«Vy'„ =
The solutions to these eciuations are, clearly,
Substituting these solutions into K(\. (13-49), we get the linal solution
of the generalized coordinates as
coordinates is
dW = ^ Fi dqi
By the use of Eq. (13-49), this work can also be expressed in terms of
the normal coordinates, becoming
dW = y y FiAij dq'j
y=ii=]
^'^ =
Wr 1(^'^''
1 = 1
(13-56)
Hence Lagrange's equations (13-15) yield, for the case of the forced
vibration of the n-degree-of-freedom system, the following set of n
independent, linear, nonhomogeneous differential equations:
n
Following now the solution given by Eq. (9-58), we see that solutions of
the n equations of Eq. (13-57) are given as
m, = 2.5 ksec/m.
..<5
;!!i
= 8,000 k/ in. *i
I
= 4
m^Q
/
/
*3 = 3,000 k/ in. /
/
/
/
X
I
t
/
•
}^//
1 kip=lk = l,0001b
(a) (6)
Ki<;. \\\-\\
and the total potential energy of the system (stored in the shear springs)
can be expressed as
Comparing these expressions for the kinetic and potential energies with
those given by Eqs. (13-31) and (13-32), we immediately find thatf
t The coefficients ai, and hij, when arranged in an array, arc usually referred to,
by analogy with the single-degree-of-freedom case, as the ?/;c.s,s, or inertia, matrix and
the stiffness matrix, respectively.
X Actually, there are more efficient methods for the calculation of the normal fre-
(|uencies and mode shapes than the straightforward apjilications of the general princi-
ples used here. These numerical methods are usually iterative in nature and often
produce the frequencies and mode shapes sinuiltaiieously. It can be appreciated
that these procedures constitute a major cITort in flie study of inulti-degree-of-free-
dom systems.
Sec. 13-7] n-degree-()f-freedom syslenis; nornuil cnnrdinales 355
.4,3= 1 .000 .4 = -
,,.H
'
1 .
7.-) 1 A r, = 0. 7.31
F = -
-=7
L
-:i t
I
\
..} --•I
"; "7
I
/ I
/
/ / /
I
/
Thu.s the general solution for free vibrations following E(j. (13-37) can
be written
For any given set of six initial conditions the six undetermined con-
stants C'l, C2, C's, 01, 02. and 03 can be easily determined by lOcjs. (13-47)
and (13-48). Suppo.se, for example, that Xi = ro; all other initial veloc-
C3 cos 03 = —7==
X
V90.55 102
j'o (2.5 X 0.2()8) C'3 sin 03 =
356 Lagrange's equations
a solution of which is
4>i = <l>2
= 4>3 =
Ci = 0.054 C2 = 0.0195 C3 = 0.007
.Ti = O.O2OV0 sin Pit - 0.0086^0 sin p^t + 0.00192^0 sin pst
X2 = 0.018^0 sin p^t + 0.00008vo sin p2t - 0.0033ro sin p-^t
and from the transformation equations of Eq. (13-50) the initial condi-
tions become
x'M = xm = x'M =
x[(0) = 2.5 X 0.3682;o
i:'2(0) = -2.5 X 0.443yo
^'3(0) = 2.5 X 0.268^0
</>l = </>2 = 03 =
C\ = 0.054yo C2 = 0.0195yo C3 = O.OOTi'o
The total solution for free vibrations can now be determined by the
use of Eq. (13-49). It is easily verified that exactly the same solution
as before will be obtained.
PROBLEMS
13-1. Suppose a particle were constrained to remain in contact with the surface
f{x,y,z) = constant. Derive the extremum conditions on the potential energy for
the equilibrium of the particle.
13-2. For the system shown, determine the force P necessary to maintain equi-
librium.
Problems 357
In.. V i;}-2
l.'}-3. In Example 13-1, suppose that kr^/Wl = 'Ki. What arc the possible ecpii-
libriiiiii eonfinurations and their corresponding classifications?
13-4. Two weightless bars of length / are connected together hy a frictionless pin.
The system issuj)ported hy a spring at the center and two frictionless rollers at the
ends. P'ind the largest value of F for stable equilibrium. The si)ring is undeformed
at the position shown.
^:r
FiG. P 13-1
13-5. Using Lagrange's efpiat ions, find the acceleration in I'rob. i'J-'Jl.
Fk;. P 13-'.»
358 Lagrange s equations
13-10. A
two-degree-of-freedom system consists of a mass ini suspended from a
A;i and a second mass mo coupled to mi by a spring of stiffness kz.
spring of stiffness
Find the natural frequencies and amplitude ratios of the two modes of vibrations.
m2
Fig. P 13-10
13-11. A uniform, slender bar of mass m and length I is supported by two springs
of equal spring constant k. The springs are situated symmetrically at a distance a
from the mid-point of the bar. Taking the vertical displacement of the mass center
of the bar and the angle of inclination of the bar as the two generalized coordinates,
write the differential equations by Lagrange's method and find the natural frequencies
of the vibration.
Fi<i. P 13-11
13-12. In the figure shown, body A, of mass nii, is free to slide in a horizontal
The other end B of the bar is attached to a second mass ?«>. A spring of spring con-
stant k is connected to bar AB at C, the other end of the spring being fixed. If
AB = I and AC = a, write the differential equations governing tlie oscillation.
13-13. A double iK'iululum is confined fo move in a plane. It consists of a rigid
body .1 hinged at <) and a rigid body B hinged to A at 0'. Tiie mass center of body
.1 is at ('a, and that of liody B at ('b. OCa = a, O'Cb = b, 00' = /. Tiu^ innineiit
of inertia of A with respect to is /i, and that of B with respect to O' is l-i. Both /'s
Problems 359
are about axes perpendicular to the plane of motion. Tiie masses of A and li are hi,
and mi, rospeetively. Derive tlie differential equations govertiinR small oscillations
of the system by means of LagranRian etiuations. (Hint: Use the two angles and ^
as the generalized coordinates.)
13-1 1. For the system shown, determine the e()Uatioiis of motion for small oscil-
lations. The mass of the compound pendulum is m>; its moment of inertia about the
center of ma.ss is 1 2- The car of mass itii is constrained to move horizontally.
13-15. In Prob. 13-10, suppose A-, = 100 lb/in., ki = 200 lb/in., w, = 8 lb secVft,
and mo = 12 Find the natural frecjuencies. the amplitude ratios, the
lb sec^/ft.
normalized mode shapes, the normal coordinates, and the respon.se of the system if
motion is initiated by giving ni-, an initial (^i^spla(•(ment of 1 in.
13-16. Verifj' that the mode shapes of Example 13-6 are orthogonal. What are
the normal coordinates in this ca.se?
13-17. Verify that the mode shapes and the normalized mode shapes of Example
13-7 are orthogonal. Suppose in this ca.se motion were initiated by imparting an
initial speed Vo to the coordinate x>. What would be the resjwn.se of the .system?
APPENDIX A
Al
A2 Properties of plane geometric shapes
Shape
ProjH'rhrs of phu\e (fconirlric shaf)es A3
*y
-I =
I = —
4r
y = z
Quarter circle
b ^\/i
nth-degree parabola
A4 Properties of plane geometric shapes
he =- sin 261
16 sin2 e)
^XCVC "
Circular segment
Circular sector
^xci/c ~ U
APPENDIX B
Izz, lyy, Jzz = mass moments of inertia tai^en with respect to x, y, z axes
fzy. I x:. I„2 = mass products of inertia taken with respect to x, y, z axes
^*x' ^tu^ i*z — mass moments of inertia taken with respect to axes pass-
ing through mass center and parallel to x, y, z axes
^ty, t*zz' Jlz = rnass products of inertia taken with respect to axes passing
through mass center and parallel to x, y, z axes
A5
A6 Properties of solid homogeneous bodies
i^y
2
?/c = 12
zc = /* = /* =0
rz Thin rod
2/c =
2C =
Ixy = /x. = /.. =
Thin circular rod
• c
'
'^s b
Rectangular prism
y
Properties of solid homogeneous bodies A7
V = ttH.
xc =
L
yc = - /* = mil
zc = I*
'xy = /*
'it = /*
'III
=0'-'
4xr'
V =
xc = Q
yc =0
zc =
* Sphere
27rr'
V =
xc
A8 Properties of solid homogeneous bodies
CHAPTER 1
cos 7 = 9/\/l53
(c) IF2 - F,| = \/T53 lb, cos a = -6/V153, cos = G/Vl-^^i. /3
cos 7 = -9/\/l53
1-14. n = (l/>/3)(i -j - k)
1 IX ^ ~ ^» _ y - ^v _ 2 - ^' _ K
CIIAI'TKK 2
CHAPTER 3
CHAPTER 4
CHAPTER 5
5-2. X =
n
—^
-\-2
b, y = ^
2«, +j—12^
r
5-21. (/ > 7 ft
CIIAPTKK 6
6-4. if - p*r =
6-5. (a) if - p=r = 0, (h) r - p=r =
6-6. 1/p = 30/(401)^2
6-7. 1-^=6+ 2ct, a, = 2f
6-8. y = a/bHl - e*') + (j'o + a/b)t + i/o
6-1 i. (a) f = (I'Vr) sin' (i>, 4> = {'2v^/r-) sin </> eo.s <^
(b) a =
6-1... (a) a, = 0. «,. = ^,)..^
^^, _^
O COS g \I
/; lA vp
6-16. = — aa> fl/i
sin ^ ( 1 H
V V ,
- a' sin' 0/
a r <* cos 9
= "
, 1 "I
H
I
Of a'co' sni' d
- "" ''''
"^
(/' a' sin' g)'^ (/'-a'sin'e)!'J
6-17. I'/. = — a)a[sin — {a/21) sin 26], ap = — a>-rt[c()s d + (a/t) cos 2S]
6-18. Vp = — aw sin g. ap = — ow' cos 5
6-19. (a) 01. = 19y°28', Or = 340°32'; (b) 0.045
rilVPTKK 7
7-3. r, = i%,r„ H,
7-4. r, = -}4,r, = -%
7-5. v^ = 0, a.4 = wVy, vfl = —uia + b)i, an = — w(o + b)i — w'oj
vc = — w(a sin + />)i + ua cos j </> </>
7-12. (a) wcD = 0, (b) o)Bc = 25 radians/sec', (r) 2 ft rij^lit of (' on the extension
of BC
7-13. (a) Vc = 34,850 in./min, (6) «Bf = 1,4.50 radians/niin
(c) ac = -28,000.000 in./min», (d) ubc = -20,000,000 radians/min»
(e) r^ = 18 in., r„ = ti \/3 in. (relative to li)
7-14. (d.\/d()„i = -w/ X A
7-1.1. t\r/(lt = U6 X r
7-16. (d\i'/dl)„i = irwCcos tf ei + sin d)
A 12 Answers to selected problems
4
7-18. as = -5/2i - 2 V3j + 2 V3 k ft/sec^, as = ^M ft/sec^
7-19. \ = —4>b sin <!> i + 4>b cos j — a)(a + b (^ cos 0)k
a = — co-(a
[ + 6 cos <^) — (p</>' cos pi)^^ cos (j> + p^^' sin pt b sin 0]i
— [ip4>' cos 7)<)^& sin <^ + p'^(t>' sin pi5 6 cos </)]j
CHAPTER 8
8-6. 7 miles/sec • ,
8-7. {mg/kyi- \
8-9. -y
W (2?;2 + ^)
i
8-14. W= 872 ft lb
I
8-15. No j
8-22. (a) ,
"^ , (6) (^- ^V^'\ (c) / = i- f ,
^
V^ |
r = 2.(^^^IHZZ«!)^^
8-23. cos <^ = %
8-24. (a) ^ + ^% = 0, (6) / = f \/^
8-25. West bank is 0.337 ft higher
8-26. (a) cos 6 = gf/w^a, (6) T = 27r/(co sin 6)
CHAPTER 9
9-3. -y/nkg/a
9-4. r = 27r V//2j7
9-5. (1 + 4a2x2).c + 4o2a:.t2 + 2ar/.r = 0, / = \/2^/2-n-
9-6. p = yf^^S
9.7. p = J^M^
V'olF
9-9. (a) 0.084 ft, (/>) O.OOOOOOO ft, (r) 17.67 ton sec/ft
9-13. (a) 0.01985 lb sec/in., (6) 13.9 radians/sec, (c) 13.88 radians/sec
(d) 0.2765 lb
—
= Fv — , . Fp
«>-! i. .t -— (\ cos pi ) sin pl\ + -p sin pi cos pli
F — —F ,
= r-
. .
X (1 t'os pt
,
) cos pl\ -\- sin pi sm />/i
CIIM'I IK 10
10-1. T = 2fo7/, y = y^ — y^gr- + /i (/i = lici^lit of shell iit hish point; io = hori-
zontal speed just before explosion: ?/o = vertical sjx'cd imparted to each
splinter by explosion)
10-2. a, = -2(/V2n-, - H'2)/(1F., + AWi), a-. = {/(2ir, - ir,)/(ir-. + AWx)
T = iWxWJiWo + 4H',)
10-3. (a) V = 25. li - 41.9j fps, {b) same as (a), (c) before impact, 7'e„r = 240
ton ft, 7',ruck = 208 ton ft; after impact, Tear = 74 ton ft, 7'truck = 185 ton ft
— — + nioVi 7»ie(t;i —
,
10-4. (a)
, ,
Vi
/
= niiVi -H m-fa
1
m^eivx Vi)
> i'2
'
= niiVi
——+ -[—
Vi)
,.> 1 — e* mim2 , ,,
2 7ni -|- 7/12
10-5. 81,000 lb
10-8. p = y/ka^/ml^ - ^/Z
10-9. <i>
= F//h(/, <^ = ^2Fiytul
10-10. (a) i.^i;(4i + 2j + k), (/*) 2)nv( - j - k), (r) 2.;„„.(i - 5j - lik)
H^ +UJ
10-1.3. V = -^^—^^-=;
\/2ff/(l - cos a)
10-17. — ^^
rr~ H ^m a
nio — kI
(4/^ - ZR^)
lo 16 V3
11-9. /max = 5, /min = 3
= -^^^^^^' - = - - -
/..
^"^
+^ (4^,2 66/ii + 3/^1^), /x,
^ (6 /iOCa h)
1) ^'' n + 3
11-22. /xx = 80.8 in.*, lyy = 38.8 in.«
/xy = -32.3 in. = -28°29' 4, cl>
CHAPTER 12
(wi — mi)g
12-26. -Cl =
Ttii + nii + I/R*
> . */6fl(l —
- cos
CO <t>o){mi -h 'im.)
12-29.
(/2 - r)(2»Hi + y/Hj)
=
— 3(; sin </>o (»ii + 2ot2)
(6) <i.
(72 - r)(2m, + Owj)
^'"'
' 2^ \(/e - r)(2m, + Oj/Ij)
12-.'il, /e^ = IL (li eos <^ - lyjiB = -r sin <^ (9 cos <^ - G)
4 4
CIIAPTKIt U
I.J-I.
_dfdv ^dVdf ^^^
_dfdV_^dVdf^^
dz i>x dz di '
dz dy dz dy
i:i-2. 02.5 U)
i.j-.i. e = 0, unstal)le; 6 = 73°3', stable
i:{-4. P = A-//2
W aQ sin 2a
3(»rx + H^b) - 2Pf^cos«a
l.{-6. 6} = m^a cos </>
/fl + 2mRa sin
7^^
— 3g sin <t>
(nil + 2?na)
13-7.
(R - r)(2ffi, + !)/«.)
13-9. = (2 _ ^/OA 9,
"
\-li/i.2 2 2»(j '-IV »ii »f2/ niimj
Judex
Index \21
Index A27
Moiiu-nts of inertia, of mass, of com- Natural fre(]uencv, of coupled systems,
posite bodies, 25(> 2 If)
units of measurement, 89
256
parallel-axes theorem for, 252, 253
Pressure-altitude relation, 87
relative, to principal axes, 2()4
Pressure forces, 89
to rotated axes, 257, 261
on curved surface, 94, 95
as tensor, 250
on flat plate, 89, 94
Mohr's circle for, 270
resultant, 90, 94
Progressive prec^ession, 'M)'2
on rocket, 234
Projectile with resistance, 184
Pressure gauge, 89
Projection of vector, on some line, 6
Pressure intensity {see Pressure)
on some plane, 49
Principal axes of inertia, 203 Propulsion shutdown, U)8
location, 2<')5
(.SV^ also Rocket burnout)
rotated moments of inertia referred to,
Propulsive force, 2S8
204 Pump, impeller, 152
by symmetry, 266, 267 mechanism, 150
Principal modes of vibration (see Normal
modes of vibration)
Principal luonients of inertia, 203, 2r)4 Quadratic forms, of kinetic energy, 342
Princii)al normal, 108 of potential energy, 341
Principal planes of inertia, 203 rethiction, 351
Principle, of angular momentum, for (Quadratic surface for inertia elii|)soid, 263
particle, 228 Quantum mechanics, 157//.
for syst<'m of particles, about fixed (^uick-return mechanism, 121
origin, 228
about moving mass center, 230
of Archimedes, 9(), 103 liadial components, of acct-leration, of
of conservation, of mechanical energy, particle, 115, 119
for particle, 173, 174 of rigiil body, 135
for rigid body, general, 284 of central force, 174
pivoted, 2S5 of velocity of particle, 115, 119
of momentum, angular, 228 Radius, of curvature, 108
linear, 225 of gyration, 251
of impuls*' and momentum, for Rate of change of vector (see Derivative
particrle, 225 of vector)
for system t)f particles, 225 IJayleigh's method, 196
of superi)osition, 214, 21S Reactions, of beam, 72
of virtual work, 329 of body in e(iuilibrium, 20, 55
of work and energy, for i)article, 171, of constraint, 22, 329
180 normal, 22, 329
A30 Index
Reactions, periodic, in bearing, 313 Resistance, of air, 184
redundant, 56 of series circuit, 193
for rotating rigid body, 313 Resisting medium, 161, 184
of rough contact, 22 Resonance, 205
of smooth contact, 22 210
in electric circuit,
workless, 329 subharmonic, 219
Reciprocating engine mechanism, 121, Resonant amplitude, 209
137 Resonant frequency, 205, 209
Recoil of guns, 226 Rest, 13
Rectangular components, of acceleration, Restitution, coefRcient, 246
111 Restoring force, 192, 340
of angular momentum, 280 Resultant, angular velocity, 126
of force, 15 external force, 53, 56, 236
of moment, 41 external moment, 54, 56
of position, 41, 111 of finite rotations, 126
of vector, 7 force, general, 68, 153
derivative, 106 parallel, 78
of velocity. 111 of infinitesimal rotations, 126
(See also Equations, of equilibrium, of Resultant moment, 68
motion) Resultant pressure force, 90, 94
Rectangular parallelepiped, moments of Retrograde precession, 303
inertia, 253 Revolution, body of, moment of inertia,
Rectilinear force, 157 266
Rectilinear motion, of particle, 154-157, motion, 298
184 Right-hand coordinate system, 6
constant velocity, 159 Right-hand-screw rule, 33, 40
free falling body, 154 Rigid bars, 18
simple harmonic, 159 colinearity of forces on, 60
of rigid body, 310 Rigid body, 54, 123
of system of particles, 234 angular momentum, about mass center,
Reduction, of angular velocities, 127 278, 279
of concurrent force systems, 69, 70 about pivot, 279
of general force system, to force and angular velocit}-, 45, 127
couple, 69 balancing, 312, 313
to wrench, 69 constraint, 54-56
of parallel force systems, 69, 78, 79 displacement, 54, 123
of quadratic forms, 351 Chasles' theorem for, 125
of rigid body motion, to rotation and Euler's theorem for, 124
translation, 130 equations of motion for, Euler's, 280,
to screw, 130 281
Redundant reaction, 56 general, 280
Reference frame (see Frames of reference) for moving coordinate axes, 281
Relative acceleration in rotating coordi- Eulerian angles for, 131-133
nate system, 145 freely moving, 277
Relative angular momentum, of rigid kinematics, point on, acceleration,
body, 278, 279 133-135
of system of particles, 230, 231 velocity, 45, 133
Relative kinetic energy, for rigid Ixxiy, theorems for, 127-131
283, 284 kinetic energy, general, 283, 284
for system of particles, 237 for pivoted, 284, 285
Relative velocity, of point on ri^iil body, plane motion, equations, 316-317
129 kinematics, 130, 131, 135
in rotating coordinate system, I 12 Poinsot's representation for motion,
Relaxation time, 200, 201" 29-4-298
ReiK'ated roots of frc(iuency e<iu:iti()n, potential energy, 284
343 rotation about fixed axis, 311, 313
Re])ulsive force, 162, 184 rotational stability, 290-293, 307
Index WW
Rigid liody, stable, oti Scalar product, gcoinelrical iulerprcla-
translation. 3ia-:ill tion, 30
typos of motion. J'.ll S<'.alar triple jiroduct, 46
work and oncrKy for, "JSJ 285 S<-alar variable. 104
(See also Air-borne vcliicle; Gyroscope) "Scotch crank," 121
Ring guide, 55 Screw, 130
/^Lr eireiiit, 193, 209 Second moment of m;is.<. 250
Roeker, 149 (See also Moments of inertia)
Rocket, engine, 234 Sections method. 73
escape velocity, 184 St'ctorial speed (areal speed), 163, 1()7
propulsion, 233-236 Semilatus rectum, 167
sf>ecific impulse, 235n. Series circuit, 19:}
thrust, 235 Series spring, 219
Rocket burnout, 168, 169 Shear force, 85
Rod motion, 302 Shear 244
stress, 85. 211.
Roller. 59 Shortest distance, from point to line. .32,
Rolling without slipping, 136, 317 36
Rolling contact. 13(5 between two nonintersecting lines, 36
lioiling cylinder. 317 Simple harmonic motion, 159, 160, 194
Rolling wheel, 135, 148 amplitude, 159
Hoot-mean-square velocity, 241, 244 fretiuency, circular, 160
Rotating frame of reference, 139, 140 period, 160
ilerivative of vector in, 140, 141 Simple pendulum, 176-180
motion of particle in, 181, 182 amplitude, 179
motion referred to, acceleration, 144, period, 179
145 as simple harmonic motion, 180
velocity, 142. 143 Simple structures, 17-21, 61
Rotation, of coordinate axes, 256-258 Simplifications, based, on order of mag-
of earth, 181 nitude, 183, 290, 309
about fixed axis, rigifl Ixxly, ((iiiations on symmetry, 26(), 290
of motion, 311, 312 Simultaneous dilTereiitial ecpiations, 157,
velocity, 45 342, 352
instantaneous axes, 127 Sine law of trigonometry, 35
instantaneous center, 131 Sliding friction (kinetic friction), 22. 176
plane, of axes, 208 Sliding vector, 2
Molir's circle for, 270, 271 Smooth contact, 329
Rotations, finite, 1, 126 Smooth surface, 22
infinitesimal, 125-127 Soft springs, 218
resultant, 127 Space axes, 132
Rough contact, 22 Space centrode, 131
Sj)ac(; cone, 297
Space structure, 74
235
Satellite, launching, Si)ecific gravity, 89
motion around earth, 164, 184. 185 Specific impulse, 2;{5«.
orbits, 1()7-170 Specific weight, 88
Scalar, examjjle, 1 of air, 101
Scalar components, of acceleration, of water, 89
111-119 Speed, 110
of vector, 7 S|)eed-time curve, 1 1 I