Engineering Mechanics Solved Problems PDF
Engineering Mechanics Solved Problems PDF
Engineering Mechanics Solved Problems PDF
-
Engineering Mechanics I
Phongsaen PITAKWATCHARA
Phongsaen.P@chula.ac.th
Contents
Preface xiv
1 Introduction to Statics 1
1.1 Basic Concepts . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Scalars and Vectors . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2.1 Representation of Vectors . . . . . . . . . . . . . . . . . . 3
1.2.2 Vector Manipulation . . . . . . . . . . . . . . . . . . . . . 3
1.2.3 Coordinate Systems . . . . . . . . . . . . . . . . . . . . . . 5
1.2.4 2-D Rectangular Coordinate System . . . . . . . . . . . . 6
1.2.5 3-D Rectangular Coordinate System . . . . . . . . . . . . 6
1.3 Newton’s Laws . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4 Gravitational Law . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 Force Systems 9
2.1 Overview of Forces . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.2 2-D Rectangular Coordinate Systems . . . . . . . . . . . . . . . . 11
2.3 2-D Force, Moment, and Couple . . . . . . . . . . . . . . . . . . . 17
2.4 2-D Resultants . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.5 3-D Rectangular Coordinate Systems . . . . . . . . . . . . . . . . 35
2.6 3-D Force, Moment, and Couple . . . . . . . . . . . . . . . . . . . 42
2.6.1 Moment about an axis . . . . . . . . . . . . . . . . . . . . 43
2.6.2 3-D couple . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
2.6.3 3-D equivalent force-couple system . . . . . . . . . . . . . 45
2.7 3-D Resultants . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
2.7.1 Resultants of special force systems . . . . . . . . . . . . . 50
2.7.2 Wrench resultant – Force-Couple equivalence . . . . . . . . 51
3 Equilibrium 57
3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.2 Mechanical System Isolation (FBD) . . . . . . . . . . . . . . . . . 58
3.3 2-D Equilibrium Conditions . . . . . . . . . . . . . . . . . . . . . 62
3.4 3-D Equilibrium Conditions . . . . . . . . . . . . . . . . . . . . . 87
CONTENTS ii
4 Friction 111
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
4.2 Types of Friction . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
4.3 Dry Friction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
4.3.1 Mechanism of Friction . . . . . . . . . . . . . . . . . . . . 112
4.3.2 Friction Cone . . . . . . . . . . . . . . . . . . . . . . . . . 115
4.3.3 Solving the Dry Friction Problems . . . . . . . . . . . . . . 115
CONTENTS iii
Bibliography 379
List of Figures
2.1 For rigid body, the force can be moved along the line of action . .
and the effect is unchanged . . . . . . . . . . . . . . . . . . . . . . 11
2.2 Components vs. orthogonal projection of a vector onto . . . . . .
the same coordinate system . . . . . . . . . . . . . . . . . . . . . 12
2.3 Parallel force addition . . . . . . . . . . . . . . . . . . . . . . . . 12
2.4 Planar force vector and its rectangular components . . . . . . . . 12
2.5 Example 2.1 ([2], pp. 29) . . . . . . . . . . . . . . . . . . . . . . . 13
2.6 Solution to example 2.1 . . . . . . . . . . . . . . . . . . . . . . . . 13
2.7 Example 2.2 ([2], pp. 30) . . . . . . . . . . . . . . . . . . . . . . . 14
2.8 Solution to example 2.2 . . . . . . . . . . . . . . . . . . . . . . . . 14
2.9 Example 2.3 ([2], pp. 32) . . . . . . . . . . . . . . . . . . . . . . 15
2.10 Solution to example 2.3 . . . . . . . . . . . . . . . . . . . . . . . . 15
2.11 Example and solution of 2.4 ([1], pp. 36) . . . . . . . . . . . . . . 16
2.12 Example 2.5 ([1], pp. 37) . . . . . . . . . . . . . . . . . . . . . . . 17
2.13 Solution to example 2.5 . . . . . . . . . . . . . . . . . . . . . . . . 17
2.14 Moment MA of F about point A . . . . . . . . . . . . . . . . . . 18
2.15 Principle of moment . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.16 Two dimensional couple . . . . . . . . . . . . . . . . . . . . . . . 19
2.17 Couple generated by equal and opposite pair of forces . . . . . . . 20
2.18 Example and solution of 2.6 ([1], pp. 44) . . . . . . . . . . . . . . 20
2.19 Example and solution of 2.7 . . . . . . . . . . . . . . . . . . . . . 21
2.20 Example 2.8 ([1], pp. 47) . . . . . . . . . . . . . . . . . . . . . . . 22
LIST OF FIGURES v
LIST OF FIGURES vi
3.1 Common action of forces in two dimensional analysis ([1], pp. 111) 59
3.2 Common action of forces in two dimensional analysis (continued)
([1], pp. 112) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
3.3 Examples of drawing the free body diagram ([1], pp. 115) . . . . . 63
3.4 Example 3.1 ([1], pp. 118) . . . . . . . . . . . . . . . . . . . . . . 64
3.5 Example 3.1 (continued) ([1], pp. 119) . . . . . . . . . . . . . . . 65
3.6 Solution to example 3.1 . . . . . . . . . . . . . . . . . . . . . . . . 66
3.7 Solution to example 3.1 (continued) . . . . . . . . . . . . . . . . . 67
3.8 Special cases of equilibrium in two dimension ([1], pp. 122) . . . . 69
3.9 Two-force member ([1], pp. 122) . . . . . . . . . . . . . . . . . . . 70
3.10 Three-force member ([1], pp. 123) . . . . . . . . . . . . . . . . . . 70
3.11 An alternative equilibrium conditions ([1], pp. 123) . . . . . . . . 71
3.12 Another alternative equilibrium conditions ([1], pp. 123) . . . . . 71
3.13 Two Dimensional Constraints and Statical Determinacy . . . . . .
([1], pp. 125) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
3.14 Example 3.2 ([1], pp. 137) . . . . . . . . . . . . . . . . . . . . . . 73
3.15 Solution to example 3.2 . . . . . . . . . . . . . . . . . . . . . . . . 74
3.16 Example 3.3 ([1], pp. 139) . . . . . . . . . . . . . . . . . . . . . . 74
3.17 Solution to example 3.3 . . . . . . . . . . . . . . . . . . . . . . . . 75
3.18 Solution to example 3.3 (continued) . . . . . . . . . . . . . . . . . 75
3.19 Example 3.4 ([2], pp. 131) . . . . . . . . . . . . . . . . . . . . . . 76
3.20 Solution to example 3.4 . . . . . . . . . . . . . . . . . . . . . . . . 77
3.21 Example 3.5 ([1], pp. 140) . . . . . . . . . . . . . . . . . . . . . . 77
3.22 Solution to example 3.5 . . . . . . . . . . . . . . . . . . . . . . . . 77
3.23 Solution to example 3.5 (continued) . . . . . . . . . . . . . . . . . 78
3.24 Example 3.6 ([1]) . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
3.25 Solution to example 3.6 . . . . . . . . . . . . . . . . . . . . . . . . 79
3.26 Solution to example 3.6 (continued) . . . . . . . . . . . . . . . . . 79
3.27 Example 3.7 ([2], pp. 136) . . . . . . . . . . . . . . . . . . . . . . 81
3.28 Solution to example 3.7 . . . . . . . . . . . . . . . . . . . . . . . . 81
3.29 Example 3.8 ([1], pp. 144) . . . . . . . . . . . . . . . . . . . . . . 83
3.30 Solution to example 3.8 . . . . . . . . . . . . . . . . . . . . . . . . 83
3.31 Solution to example 3.8 (continued) . . . . . . . . . . . . . . . . . 83
3.32 Example 3.9 ([2], pp. 137) . . . . . . . . . . . . . . . . . . . . . . 84
3.33 Solution to example 3.9 . . . . . . . . . . . . . . . . . . . . . . . . 84
3.34 Example 3.10 ([2], pp. 141) . . . . . . . . . . . . . . . . . . . . . 85
3.35 Solution to example 3.10 . . . . . . . . . . . . . . . . . . . . . . . 86
3.36 Solution to example 3.10 (continued) . . . . . . . . . . . . . . . . 86
5.1 Some examples of the distributed forces ([1], pp. 237) . . . . . . . 136
5.2 Examples of symmetric objects ([1], pp. 239) . . . . . . . . . . . . 138
5.3 Determination of the line centroid ([1], pp. 240) . . . . . . . . . . 139
5.4 Determination of the area centroid ([1], pp. 240) . . . . . . . . . . 139
5.5 Select the differential element, of which the volume is known, . . .
to be as large as possible ([1], pp. 241) . . . . . . . . . . . . . . . 141
5.6 If possible, divide the object in the manner that can be integrated
in one continuous operation to cover the whole body ([1], pp. 241) 142
5.7 The higher order terms can safely be omitted from the integrand .
([1], pp. 241) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 143
5.8 Select the coordinates that suit to the shape and boundaries . . .
of the object ([1], pp. 242) . . . . . . . . . . . . . . . . . . . . . . 143
5.9 Divide the object in the manner that the centroid be expressed . .
easily in terms of the selected coordinates ([1], pp. 242) . . . . . . 144
5.10 Example 5.1 ([1], pp. 245) . . . . . . . . . . . . . . . . . . . . . . 144
5.11 Solution to example 5.1 ([1], pp. 245) . . . . . . . . . . . . . . . . 145
5.12 Example 5.2 ([2], pp. 258) . . . . . . . . . . . . . . . . . . . . . . 146
5.13 Solution to example 5.2 . . . . . . . . . . . . . . . . . . . . . . . . 147
5.14 Example 5.3 ([2], pp. 261) . . . . . . . . . . . . . . . . . . . . . . 148
5.15 Solution to example 5.3 . . . . . . . . . . . . . . . . . . . . . . . . 149
5.16 Example 5.4 ([2], pp. 262) . . . . . . . . . . . . . . . . . . . . . . 150
5.17 Solution to example 5.4 . . . . . . . . . . . . . . . . . . . . . . . . 150
5.18 The centroid of a complex body may be determined discretely . .
([1], pp. 256) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 151
5.19 Example 5.5 ([2], pp. 269) . . . . . . . . . . . . . . . . . . . . . . 152
5.20 Solution to example 5.5 . . . . . . . . . . . . . . . . . . . . . . . . 153
5.21 Example 5.6 ([2], pp. 274) . . . . . . . . . . . . . . . . . . . . . . 154
5.22 Solution to example 5.6 . . . . . . . . . . . . . . . . . . . . . . . . 155
5.23 Creation of the revolved surface object ([1], pp. 266) . . . . . . . 155
5.24 Creation of the revolved volume object ([1], pp. 266) . . . . . . . 156
5.25 Example 5.7 ([2], pp. 280) . . . . . . . . . . . . . . . . . . . . . . 157
5.26 Example 5.8 ([1], pp. 273) . . . . . . . . . . . . . . . . . . . . . . 159
LIST OF FIGURES ix
5.27 Free body diagram of an infinitesimal fluid ([1], pp. 308) . . . . . 159
5.28 A flat surface immersed in the fluid medium with the . . . . . . .
pertinent coordinates ([1], pp. 312) . . . . . . . . . . . . . . . . . 161
5.29 Pressure distribution overlaying the flat surface ([1], pp. 312) . . . 162
5.30 Side view profile of the pressure distribution ([1], pp. 310) . . . . 163
5.31 A cylindrical surface with constant width immersed in the . . . .
fluid medium ([1], pp. 311) . . . . . . . . . . . . . . . . . . . . . . 164
5.32 Side view profile of the pressure distribution on the . . . . . . . .
cylindrical surface ([1], pp. 311) . . . . . . . . . . . . . . . . . . . 165
5.33 Free body diagram of a block of liquid column over the . . . . . .
cylindrical surface ([1], pp. 311) . . . . . . . . . . . . . . . . . . . 167
5.34 Concept of buoyancy force ([1], pp. 314) . . . . . . . . . . . . . . 168
5.35 Example 5.9 ([1], pp. 320) . . . . . . . . . . . . . . . . . . . . . . 169
5.36 Solution to example 5.9 . . . . . . . . . . . . . . . . . . . . . . . . 169
5.37 Example 5.10 ([2], pp. 328) . . . . . . . . . . . . . . . . . . . . . 170
5.38 Solution to example 5.10 . . . . . . . . . . . . . . . . . . . . . . . 171
5.39 Example 5.11 ([1]) . . . . . . . . . . . . . . . . . . . . . . . . . . 171
5.40 Solution to example 5.11 . . . . . . . . . . . . . . . . . . . . . . . 172
5.41 Example 5.12 ([2], pp. 329) . . . . . . . . . . . . . . . . . . . . . 173
5.42 Solution to example 5.12 . . . . . . . . . . . . . . . . . . . . . . . 174
5.43 Example 5.13 ([1], pp. 322) . . . . . . . . . . . . . . . . . . . . . 175
5.44 Solution to example 5.13 . . . . . . . . . . . . . . . . . . . . . . . 175
5.45 Example 5.14 ([1], pp. 324) . . . . . . . . . . . . . . . . . . . . . 176
5.46 Solution to example 5.14 . . . . . . . . . . . . . . . . . . . . . . . 177
5.47 Example 5.15 ([1], pp. 325) . . . . . . . . . . . . . . . . . . . . . 178
5.48 Solution to example 5.15 . . . . . . . . . . . . . . . . . . . . . . . 179
5.49 Example 5.16 ([1], pp. 333) . . . . . . . . . . . . . . . . . . . . . 179
5.50 Solution to example 5.16 . . . . . . . . . . . . . . . . . . . . . . . 180
7.1 Several descriptions of the particle position ([3], pp. 22) . . . . . . 186
7.2 Rectilinear motion showing the distance s and its change ∆s . . .
measured along the straight line path ([3], pp. 22) . . . . . . . . . 188
7.3 Plots of the displacement, velocity, and acceleration with respect .
to time and their related quantities ([3], pp. 24) . . . . . . . . . . 190
7.4 Example 7.1 ([3], pp. 33) . . . . . . . . . . . . . . . . . . . . . . . 193
7.5 Solution to example 7.1 . . . . . . . . . . . . . . . . . . . . . . . . 194
7.6 Example 7.2 ([4], pp. 30) . . . . . . . . . . . . . . . . . . . . . . . 195
7.7 Solution to example 7.2 . . . . . . . . . . . . . . . . . . . . . . . . 195
7.8 Example 7.3 ([3], pp. 34) . . . . . . . . . . . . . . . . . . . . . . . 197
7.9 Solution to example 7.3 . . . . . . . . . . . . . . . . . . . . . . . . 197
7.10 Example 7.4 ([4], pp. 32) . . . . . . . . . . . . . . . . . . . . . . . 199
7.11 Example 7.5 ([3], pp. 37) . . . . . . . . . . . . . . . . . . . . . . . 200
7.12 Example 7.6 ([4], pp. 35) . . . . . . . . . . . . . . . . . . . . . . . 201
LIST OF FIGURES x
LIST OF FIGURES xi
9.1 Common plane motion of rigid bodies ([3], pp. 333) . . . . . . . . 307
9.2 Rotation motion of a rigid body ([3], pp. 333) . . . . . . . . . . . 309
9.3 Rotation about a fixed axis viewing from top ([3], pp. 335) . . . . 310
9.4 Rotation about a fixed axis in perspective view ([3], pp. 335) . . . 310
9.5 Example 9.1 ([3], pp. 342) . . . . . . . . . . . . . . . . . . . . . . 311
9.6 Solution to example 9.1 . . . . . . . . . . . . . . . . . . . . . . . . 312
9.7 Example 9.2 ([3], pp. 343) . . . . . . . . . . . . . . . . . . . . . . 313
9.8 Example 9.3 ([3], pp. 345) . . . . . . . . . . . . . . . . . . . . . . 315
9.9 Applications of rolling without slipping ([3], pp. 345) . . . . . . . 315
9.10 Example 9.4 ([3], pp. 346) . . . . . . . . . . . . . . . . . . . . . . 316
9.11 Solution to example 9.4 ([3], pp. 346) . . . . . . . . . . . . . . . . 316
9.12 Example 9.5 ([3], pp. 350) . . . . . . . . . . . . . . . . . . . . . . 317
9.13 Solution to example 9.5 . . . . . . . . . . . . . . . . . . . . . . . . 317
9.14 Example 9.6 ([4], pp. 340) . . . . . . . . . . . . . . . . . . . . . . 318
9.15 Solution to example 9.6 . . . . . . . . . . . . . . . . . . . . . . . . 319
9.16 Example 9.7 ([4], pp. 342) . . . . . . . . . . . . . . . . . . . . . . 320
9.17 Example 9.8 ([3], pp. 354) . . . . . . . . . . . . . . . . . . . . . . 321
9.18 Solution to example 9.8 . . . . . . . . . . . . . . . . . . . . . . . . 322
9.19 Example 9.9 ([3], pp. 355) . . . . . . . . . . . . . . . . . . . . . . 323
9.20 Example 9.10 ([4], pp. 345) . . . . . . . . . . . . . . . . . . . . . 324
List of Tables
5.1 Table of areas, their moments, and centroids of fig. 5.19 . . . . . . 152
5.2 Table of areas, their moments, and centroids of fig. 5.21 . . . . . . 153
Preface
Chapter 1
Introduction to Statics
• Time is the measure of the succession of events. Often, we are more inter-
ested in the change of physical quantities with respect to time, e.g. v = dr
dt
,
instead of time variable itself.
• Rigid body is a body whose relative movement between its parts are
negligible relative to the gross motion of the body. For example the motion
of an ingot can be analyzed by assuming the object being rigid.
• Nonrigid body is a body whose relative movement between its parts are
significant relative to the gross motion of the body. Knowledge of the
mechanics of the deformable material must be used along with Dynamics
in order to determine the absolute motion of the nonrigid bodies.
Let us consider some examples to see the difference of each term. If we have
an object and consider the very small substance of the body. For differential
element analysis of the body, the small substance can be treated as a particle.
However, the substance must be handled as connecting objects had the molecular
effects in the body are of concern. Or think of an airplane. Even of its huge size,
the whole airplane may be modeled as a point in flight speed analysis along the
route. But if the rotational motion, such as yawing or pitching, of the airplane
body is important, its size does matter.
The next two examples are to show whether an object is considered rigid or
nonrigid depends on how much detailed of the problem we would like to analyze.
Truss can just be looked as a rigid body for the preliminary design of truss
structure. But we must think of the truss elasticity if we were to choose the
material for that truss. A stiff linkage of the robot may be considered a rigid
body. However, the n-connecting linkages, treated as a whole, to form the robot
arm is an example of nonrigid body. Note the body-fixed inertia of the nonrigid
body is not constant.
• Scalars are quantities for which only the magnitude can describe com-
pletely. Time, volume, density, speed, energy, and mass are some examples.
• Vectors are quantities for which both the magnitude and the direction
are needed to completely describe. Examples are displacement, velocity,
acceleration, force, moment, and momentum. Vectors cand be classified
into 3 types: free vector, sliding vector, and fixed vector.
• Free Vector is a vector whose action is not confined with a unique line in
space. That is, only its magnitude and direction do matter. Some examples
are the displacement vector of a pure translational rigid object, or the couple
vector of a rigid body. Free vector is free to slide and translate as long as its
direction and magnitude are maintained. In other words, its line of action
and point of application do not matter.
Figure 1.2: Parallelogram law and Head to tail for vector addition ([1], pp. 6)
points coincide. Then parallelogram law is applied to find the resultant vector,
v. Vectors can also be added up using the head to tail, as done in the last figure.
However, in contrary to the parallelogram method, the head to tail method does
not guarantee correct line of action of the resultant vector. It may be obtained
separately with the help of the principle of moment. Neither method gives the
correct point of application. In summary, only the magnitude and direction of
the resultant vector are ensured. If one choose the parallelogram law, the correct
line of action can be obtained as well.
Some familiar algebraic laws also hold for the vector addition operation. They
are
In the algebraic approach, cosine law and sine law are used in determining the
magnitude and direction of the resultant vector from the addition of two vectors,
as shown in Fig. 1.3. The resultant vector, drawn in red, is one side of the triangle
Figure 1.3: Determination of the magnitude and direction of the resultant vector
Figure 1.4: Vector components along different coordinate systems ([1], pp. 6)
forming from the other two vectors, v1 and v2 . The magnitude, i.e. the length,
of the resultant vector v is found by the cosine law:
v 2 = v12 + v22 − 2v1 v2 cos β (1.1)
And the direction can be determined from the sine law and pruning the impossible
solution:
v v2
= (1.2)
sin β sin θ
It can be concluded from the equation that the three angles are dependent. Only
two of them are enough to specify the orientation of the vector.
F = attraction force
gravitational constant value = 6.673 × 10−11 m3 / kg · s2
G =
m = mass of the involving bodies
r = distance between the bodies
Hence there is always the attraction force between the earth and the object.
This gravitational force is called the weight of the body.
Gme
W =m = mg (1.11)
r2
where
Chapter 2
Force Systems
Figure 2.1: For rigid body, the force can be moved along the line of action and
the effect is unchanged
Figure 2.2: Components vs. orthogonal projection of a vector onto the same
coordinate system
F = Fx + Fy = Fx i + Fy j (2.1)
The components are the orthogonal projection of the vector onto the respective
axes which are determined by the dot product of the vector and the unit vector
along the axes.
Fx = F · i = F cos θ
(2.2)
Fy = F · j = F sin θ
The magnitude and direction of the force vector F follow immediately as
q
F = Fx2 + Fy2 (2.3)
θ = arctan2 (Fy , Fx ) (2.4)
Example 2.1 ([2], Prob. 2/9) If the two equal tension T in the pulley cable
together produce a force of 5 kN on the pulley bearing, calculate T.
Solution: Use the parallelogram law and the cosine law to determine the
non-orthogonal components. From the force vector addition and by the cosine
law,
52 = T 2 + T 2 + 2T × T cos 60◦
T = 2.89 kN
Example 2.2 ([2], Prob. 2/11) While steadily pushing the machine up an
incline, a person exerts a 180 N force P as shown. Determine the components of
P which are parallel and perpendicular to the incline.
Solution: Draw the coordinate axes and the force vector. Carefully indi-
cate the angles. Then project the force to the respective axes.
Pt = 180 cos (10 + 15) = 163.1 N
Pn = −180 sin (10 + 15) = −76.1 N
Example 2.3 ([2], Prob. 2/19) Determine the resultant R of the two forces
applied to the bracket. Write R in terms of unit vectors along the x- and y- axes
shown.
Solution: Draw the coordinate axes and the force vectors. Carefully indi-
cate the angles. After that, project the forces to the respective axes. Then
algebraically add the components that are on the same axis to obtain the answer.
The problem is extended to determine the components along the non-orthogonal
coordinate system.
Force components along the x-y coordinate system are
Rx = 200 cos (15 + 20) − 150 sin (10 + 20) = 88.8 N
Ry = 200 sin (15 + 20) + 150 cos (10 + 20) = 244.6 N
Example 2.4 ([1], Prob. 2/20) It is desired to remove the spike from the timber
by applying force along its horizontal axis. An obstruction A prevents direct
access, so that two forces, one 1.6 kN and the other P, are applied by cables as
Solution: To remove the spike, the efforted force must point in the direc-
tion along the spike axis. This means there is no net force in the perpendicular
direction. P and 1.6 kN must add up to give the resultant force in the horizontal
direction.
Because there is no net force in the y- direction, the summation of force is
100 150
Ry = P sin atan − 1.6 sin atan =0
200 200
Therefore,
P = 2.15 kN
The axial tension is the total force along the x- direction.
100 150
T = Rx = P cos atan + 1.6 cos atan = 3.20 kN
200 200
Example 2.5 ([1], Prob. 2/26) As it inserts the small cylindrical part into
a close fitting circular hole, the robot arm exerts a 90 N force P on the part
parallel to the axis of the hole as shown. Determine the components of the force
which the part exerts on the robot along axes (a) parallel and perpendicular to
the arm AB, and (b) parallel and perpendicular to the arm BC.
Solution: The indicated force P is the force done by the robot on the
cylindrical part. Therefore the force exerted by the part on the robot is −P. We
set up the coordinate frame n1 t1 and n2 t2 where their axes are perpendicular
and parallel to the link AB and BC, respectively. Then −P is projected onto
these rectangular coordinate system.
Force done by the part on the robot written in n1 t1 frame is
−P = −90 cos 45n1 + 90 sin 45t1 = −63.6n1 + 63.6t1 N
MA = r × F (2.5)
where r is the position vector from A to any point along the line of action of F.
If only the magnitude is considered, the formula can be written as
MA = F · r (sin α) = F · d (2.6)
“The moment of a force about any point is equal to the sum of the
moments of the components of the force about the same point”
This theorem is helpful in determining the moment of the force from its com-
ponents. In many cases, moments of some components may be trivial to calculate.
M = rA × F + rB × (−F) = r × F (2.8)
M = F · r (sin α) = F d (2.9)
The couple vector’s direction is perpendicular to the plane established by those
two lines of action of the forces. It is a free vector and so no moment center.
Only the magnitude and direction are enough to describe the couple.
For rigid body, several pairs of equal and opposite forces can give the same
couple. Therefore it is unique to calculate the couple from a given pair of forces
but it is non-unique to determine the pair of forces which will produce that value
of couple. See fig. 2.17. This property can be applied in determining the effect
of the couple from the equivalent pair of forces. Effect from some specific pair of
forces may be trivial to calculate.
Example 2.6 ([1], Prob. 2/36) Calculate the moment of the 250 N force on the
handle of the monkey wrench about the center of the bolt.
Solution: There are many ways in determining the answer. The most
suitable way, however, is to factor the force into components where the moments
can be found easily. Varignon’s theorem is then used in calculating the desired
moment from the sum of moments of their components.
Project 250 N force into x-y coordinate frame and sum the moments induced
by these components about the center of the bolt.
Example 2.7 Calculate the moment of the 240 N force on the handle of the
prong about the instantaneous supporting point O.
and F are described as vectors according to the coordinate system. The moment
can then be calculated directly as MO = r × F.
r = 0.03i + 0.35j m
F = 240 cos 10i − 240 sin 10j N
MO = r × F = −84.0k Nm
Example 2.8 ([1], Prob. 2/46) The force exerted by the plunger of cylinder
AB on the door is 40 N directed along the line AB, and this force tends to
keep the door closed. Compute the moment of this force about the hinge
O. What force FC normal to the plane of the door must the door stop at C
exert on the door so that the combined moment about O of the two forces is zero?
FC = MO /0.825 = 8.53 N
Example 2.9 ([1], Prob. 2/52) While inserting a cylindrical part into the
circular hole, the robot exerts the 90 N force on the part as shown. Determine
the moment about point A, B, and C of the force which the part exerts on the
robot.
From fig. 2.23, the position vector rAC written in x-y coordinate system is
rAC = (0.55 cos 60 + 0.45 cos 45) i+(0.55 sin 60 − 0.45 sin 45) j = 0.593i+0.158j m
M−P at C about A
= rAC × F = 55.23 Nm CCW
Therefore moment of the force which the part exerts on the robot about A is
MA = MC + M−P at C about A
= 68.7 Nm CCW
Example 2.10 ([2], Prob. 2/55) As part of a test, the two aircraft engines are
revved up and the propeller pitches are adjusted so as to result in the fore and aft
thrusts shown. What force F must be exerted by the ground on each of the main
braked wheels at A and B to counteract the turning effect of the two propeller
thrusts? Neglect any effect of the nose wheel C, which is turned 90◦ and unbraked.
Solution: Calculate the couple generated by the thrust forces and equate
it to the one produced by the braked forces. The resultant couple is zero, that is,
MC = 2 × 5 − F × 3 = 0
F = 3.33 kN
Example 2.11 ([2], Prob. 2/59) A lug wrench is used to tighten a square-head
bolt. If 250 N forces are applied to the wrench as shown, determine the
magnitude F of the equal forces exerted on the four contact points on the 25
mm bolt head so that their external effect on the bolt is equivalent to that of
the two 250 N forces. Assume that the forces are perpendicular to the flats of
the bolt head.
Solution: Calculate the couple generated by the two 250 N forces. Four
forces exerted at the sides of the square-head bolt are such that they give the
negative value to balance the applied couple. The equivalent couple system at
the bolt head is shown in fig. 2.26.
250 × 0.7 = 2 (F × 0.025)
F = 3500 N
of three figures are equivalent in the sense that they make the rigid body moves
in the same manner. To repeat, force F acting at point A is equivalent to the
force F acting at point B plus the compensated couple M.
Resultant is the simplest force combination which can replace the original
system of forces, moments, and couples without altering the external effect of the
system on the rigid body. Several tools are used in determining the resultant.
accuracy. The magnitude, direction, and line of action obtained by this method
are correct. Figure 2.29 shows an example of determining the resultant force by
the principle of transmissiblity and the parallelogram law.
2. Move all forces so the new lines of action pass through point O. The
new lines are parallel to the original ones. By the force-couple equivalent
method, the couples are generated to preserve the motion effect of the
original force system.
4. Locate the correct line of action of the resultant force R. This is done by
applying the inverse of force-couple equivalent method. Force R must be
moved to a new location so that the generated couple cancels the couple
MO . Only the force is left and hence is the resultant force.
Figure 2.30: Step 1 and 2 of the force-couple equivalent method ([1], pp. 59)
Figure 2.31: Step 3 and 4 of the force-couple equivalent method ([1], pp. 59)
Example 2.13 ([2], Prob. 2/61) Calculate the moment of the 1200 N force
about pin A of the bracket. Begin by replacing the 1200 N force by a force-couple
system at point C. Calculate the moment of the 1200 N force about the pin at B.
Example 2.14 ([1], Prob. 2/70) The combined drive wheels of a front-wheel-
drive automobile are acted on by a 7000 N normal reaction force and a friction
force F, both of which are exerted by the road surface. If it is known that the
resultant of these two forces makes a 15◦ angle with the vertical, determine the
equivalent force-couple system at the car mass center G. Treat this as a 2D
problem.
Solution: First, determine the resultant force of the normal and friction
forces. Since the resultant force makes a 15◦ angle with the vertical,
Then move the resultant force acting at the front wheel to the car mass center
G. Consequently, the couple must be provided to compensate for the induced
moment of the original force system.
MG = 7000 × 1 + 7246.9 sin 15 × 0.5 = 7937.8 Nm CW
Example 2.15 ([2], Prob. 2/75) Determine and locate the resultant R of the
two forces and one couple acting on the I-beam.
Solution: We first select an arbitrary point O to which all forces and cou-
ples will be moved to. Then the equivalent force-couple at point O is
determined.
R = 8 − 5 = 3 kN downward
MO = 25 − 5 × 2 − 8 × 2 = 1 kNm CW
Finally, the resultant force is found by locating the correct line of action. This
step is essentially done by applying the principle of moment.
1
3d = 1, d = 1/3 m & x = 4 m
3
Example 2.16 ([1], Prob. 2/83) If the resultant of the two forces and couple
M passes through point O, determine M.
Solution: The key to this problem is that the resultant passes through
point O means there is no moment at point O. By summing all moment from
the two forces and one unknown couple M at point O to zero, we can find the
solution.
MO = M − 400 × 0.15 cos 30 − 320 × 0.3 = 0
M = 148 Nm CCW
Example 2.17 ([1], Prob. 2/88) The directions of the two thrust vectors of
an experimental aircraft can be independently changed from the conventional
forward direction within limits. For the thrust configuration shown, determine
the equivalent force-couple system at point O. Then replace this force-couple
system by a single force and specify the point on the x-axis through which the
line of action of this resultant passes.
Example 2.18 ([1], Prob. 2/93) Two integral pulleys are subjected to the belt
tensions shown. If the resultant R of these forces passes through the center O, de-
termine T and the magnitude of R and the CCW angle θ it makes with the x-axis.
Solution: The key to this problem is that the resultant passes through
R = (200 + 150 − 160 cos 30 − 60 cos 30) i + (160 sin 30 + 60 sin 30) j
Example 2.19 ([1], Prob. 2/97) A rear-wheel-drive car is stuck in the snow
between other park cars as shown. In an attempt to free the car, three students
exert forces on the car at point A, B, and C while the driver’s actions result
in a forward thrust of 200 N acting parallel to the plane of rotation of each
rear wheel. Treating the problem as 2D, determine the equivalent force-couple
system at the car center of mass G and locate the position x of the point on the
car centerline through which the resultant passes. Neglect all forces not shown.
Solution: The resultant force is simply the addition of all forces acting on
the car. Using the given x-y coordinate system,
R = (200 + 400 + 200 + 250 sin 30) i + (250 cos 30 + 350) j = 925i + 566.5j N
It is obvious that the 400 N force in x-direction and the y-component of 250
N cause no moment about G. Therefore the moment at point G is
Note that the x-component of the resultant force does not contribute the
moment about G. Hence only Ry will be used in the moment calculation, from
which the position on the centerline of the car where the resultant passes can be
deduced.
566.5 × x = 690, x = 1.218 m
Example 2.20 ([1], Prob. 2/98) An exhaust system for a pickup truck is
shown in the figure. The weights Wh , Wm , and Wt of the headpipe, muffler, and
tailpipe are 10, 100, and 50 N, respectively, and act at the indicated points. If
the exhaust pipe hanger at point A is adjusted so that its tension FA is 50 N,
determine the required forces in the hangers at point B, C, and D so that the
force-couple system at point O is zero. Why is a zero force-couple system at O
desirable.
Solution: The key is that for static, if the force-couple at point O is zero,
the force-couple at any point is zero too! In particular, the moment at point E
is zero. That is
Wh ×(0.2 + 1.3 + 0.9)+Wm ×(0.65 + 0.9)+Wt ×0.4−FA ×(1.3 + 0.9)−FB ×0.9 = 0
FB = 98.9 N
We use the knowledge that the force components in horizontal and vertical
direction are zero to determine the remaining unknown forces, FC and FD .
FA + FB + FC cos 30 + FD cos 30 − Wh − Wm − Wt = 0
FD sin 30 − FC sin 30 = 0
FC = FD = 6.415 N
Therefore the pipe is in equilibrium without external reaction force at support
O. The stress at O, which is at the joint, is zero. This effectively prevents the
breakage of the exhaust system.
F = Fx + Fy + Fz = Fx i + Fy j + Fz k (2.14)
The components are the orthogonal projection of the vector onto the respective
axes which are determined by the dot product of the vector and the unit vector
along the axes.
Fx = F · i = F cos θx
Fy = F · j = F cos θy (2.15)
Fz = F · k = F cos θz
The magnitude and direction of the force vector F follow immediately as
q
F = Fx2 + Fy2 + Fz2 (2.16)
θx = arccos (Fx , F )
θy = arccos (Fy , F ) (2.17)
θz = arccos (Fz , F )
And the following relationship can be seen directly from the figure;
Fxz = F sin θy
Fxy = F sin θz (2.18)
Fyz = F sin θx
Figure 2.43: Force vector and its rectangular components ([1], pp. 66)
We define the directional unit vector, nF to be the unit vector of which its
direction is the same as the vector F:
F = F nF (2.20)
Figure 2.44: Direction of force vector by two points ([1], pp. 66)
Orthogonal projection
The magnitude of the orthogonal projection of F in the n-direction is deter-
mined by the dot product of F with n. Mathematically,
Fn = F · n = F nF · n = F cos θ (2.23)
Fn = Fn n = (F · n) n (2.25)
Example 2.21 ([2], Prob. 2/99) In opening a door which is equipped with a
heavy duty return mechanism, a person exerts a force P of magnitude 32 N as
shown. Force P and the normal n to the face of the door lie in a vertical plane.
Express P as a vector and determine the angles θx , θy , and θz which the line of
action P makes with the positive x-, y-, and z-axes.
Figure 2.45: Direction of force vector by two angles ([1], pp. 67)
Figure 2.46: Orthogonal projection of F onto the n-direction ([1], pp. 67)
onto the z-axis and the x-y plane. Then the force in the x-y plane is further
decomposed into the x- and y- axes, respectively.
Example 2.22 ([1], Prob. 2/112) The rectangular plate is supported by hinges
along its side BC and by the cable AE. If the cable tension is 300 N, determine
the projection onto line BC of the force exerted on the plate by the cable.
Note that E is the midpoint of the horizontal upper edge of the structural support.
A = (−0.4, 0, 1.2 sin 25) B = (0, 0, 1.2 sin 25) C = (0, 1.2 cos 25, 0)
Example 2.23 ([2], Prob. 2/104) The power line is strung from the power-pole
arm at A to point B on the same horizontal plane. Because of the sag of the
cable in the vertical plane, the cable makes an angle of 15◦ with the horizontal
where it attaches to A. If the cable tension at A is 800 N, write T as a vector
and determine the magnitude of its projection onto the x-z plane.
angle θ so the new orientation is along AB. Finally the line is vertically rotated
downward by 15◦ and hence its direction is the same as that of T. With the
known directional unit vector, we can write T in x-y-z frame. The magnitude of
its projection onto the x-z plane is followed immediately.
The angle θ is
θ = atan (1.5/10) = 8.53◦
Therefore the tension force T can be described as
It can also be determined by first calculating the angle that T made with the
y-axis:
T·j
θy = acos = 81.76◦
T
Referring to fig. 2.43 and eq. 2.18,
MO = r × F (2.26)
MO = r × F = (rx i + ry j + rz k) × (Fx i + Fy j + Fz k)
= (ry Fz − rz Fy ) i + (rz Fx − rx Fz ) j + (rx Fy − ry Fx ) k
(2.27)
This formula can be remembered easily by recognizing the expression as the
following determinant;
i j k
r × F ⊲ rx ry rz (2.28)
Fx Fy Fz
Figure 2.52: Moment determination by the principle of moment ([1], pp. 75)
We can proofcheck this vector cross product and gain understanding in the
process of determining the moment by applying the principle of moment. Ac-
cording to fig 2.52,
Mx = ry Fz − rz Fy
My = rz Fx − rx Fz (2.29)
Mz = rx Fy − ry Fx
which confirms the mathematical operation.
Mλ = (MO · n) n = (r × F · n) n (2.30)
There is an interesting point of the moment about an axis. Point O can be any
point on the axis λ and it still gives equal value of Mλ even, of course, MO is
distinct for each of the chosen point O.
Figure 2.55: Three dimensional equivalent force-couple system ([1], pp. 77)
2. Move the force so the new line of action pass through B. The new line
are parallel to the original one. By the force-couple equivalent method, the
couple must be generated to preserve the motion effect of the original force
system. If r is the position vector from B to any point on the line of action
of F, the compensated couple is M = r × F.
Figure 2.55 illustrates the procedure.
Example 2.24 ([1], Prob. 2/123) The helicopter is drawn here with certain
3-D geometry given. During a ground test, a 400 N aerodynamic force is applied
to the tail rotor at P as shown. Determine the moment of this force about point
O of the airframe.
Solution: For this simple force P, we can determine the moment component-wise.
It is obvious that the force is in y-direction. Hence it does not cause moment in
the y-direction. By observation,
MO = (400 × 1.2) i + (400 × 6) k = 480i + 2400k N
Example 2.25 ([2], Prob. 2/118) In picking up a load from position A, a cable
tension T of magnitude 21 kN is developed. Calculate the moment that T
produces about the base O of the construction crane.
Example 2.26 ([2], Prob. 2/122) The special-purpose milling cutter is sub-
jected to the force of 1200 N and a couple of 240 Nm as shown. Determine the
moment of this system about point O.
r = 0.2i + 0.25k m
The moment of this system about O is the summation of the moment induced
by the force and the free vector couple;
Example 2.27 ([1], Prob. 2/141) A 5 N vertical force is applied to the knob of
the window-opener mechanism when the crank BC is horizontal. Determine the
moment of the force about point A and about the line AB.
Solution: As usual, describe the related quantities using the same coordi-
nate frame.
r = 75 cos 30i + 75j + 75 sin 30k mm
F = −5k N
The moment about the line AB is then the component of MA about AB, which
can be found by the dot product.
X
M= (r × F) (2.33)
4. Locate the correct line of action of the resultant R. This step is the reverse
of the force-couple equivalence. The goal is to solve for the piercing point
of the resultant with the object. In 3-D problem, it involves determining 2
unknowns in the moment equation, r × R = M. The matrix equation is of
rank 2. That is one of the equation is degenerated.
and X
M= (r × F) = 0 (2.39)
Parallel forces If the lines of action of all forces are parallel, the resultant force
can be obtained algebraically. The magnitude of the resultant force, R, is the
magnitude of the algebraic sum of the given forces.
X
|R| = F (2.40)
Wrench resultant If the resultant force and the resultant moment, as shown
in fig. 2.61, the resultant system is called wrench resultant. An example is the
reaction force and moment of the screwdriver. If the force and moment are in the
same direction, it is called positive wrench. Otherwise it is the negative wrench.
Wrench resultant is the simplest form to visualize the effect of general force
system on to the object. That is the object is simultaneously translating and
rotating about the unique axis: the screw axis.
Figure 2.61: A special force system: wrench system ([1], pp. 89)
3. Transform the couple M2 into the equivalent pair of R and −R, with −R
applied at O to cancel R, as shown in fig. 2.62c.
Example 2.28 ([1], Prob. 2/148) The pulley and gear are subjected to the loads
shown. For these forces, determine the equivalent force-couple system at point O.
R = (800 + 200 − 1200 sin 10) i + 1200 cos 10j = 792i + 1182j N
Therefore the resultant moment is the sum of the moments of three forces:
Example 2.29 ([1], Prob. 2/152) Two upward loads are exerted on the small
3D truss. Reduce these two loads to a single force-couple system at point O.
Show that R is perpendicular to MO . Then determine the point in the x-z plane
through which the resultant passes.
Solution: Since the moment involves the cross product of the position and
force vector, R must be perpendicular to MO . However we have to show it
explicitly for this problem.
By inspection,
Consequently, the resultant force R passes through (2.4, −0.6) m in the x-z plane.
Example 2.30 ([2], Prob. 2/141) Replace the two forces acting on the block
by a wrench. Write the moment M associated with the wrench as a vector and
specify the coordinates of the point P in the x-y plane through which the line of
action of the wrench passes.
R = Fi − Fk
MO = F (a + c) j − F bk
Next, we project MO onto the direction parallel and perpendicular to nR , for
which its value is
R 1 1
nR = = √ i− √ k
|R| 2 2
Consequently, the components of MO are
Fb Fb
Mk = (MO · nR ) nR = i− k
2 2
Fb Fb
M⊥ = MO − Mk = − i + F (a + c) j − k
2 2
Thirdly, we transform the couple M⊥ into pair of forces R and −R. Suppose
r is the position vector of the desired piercing point in the x-y plane which can
be described as
r = xi + yj
The moment equation is formed to transform between the pair of forces and the
couple;
M⊥ = r × R
in which the values of M⊥ and R are substituted to solve for r. The coordinates
are
x = a + c, y = b/2
Finally, the wrench consists of R and Mk which pass through the x-y plane
at x = a + c and y = b/2.
Example 2.31 ([1], Prob. 2/159) The resultant of the two forces and couple
may be represented by a wrench. Determine the vector expression for the
moment M of the wrench and find the coordinates of the point P in the x-z
plane through which the resultant force of the wrench passes.
R = 100i + 100j N
Another way to determine the wrench is to assume the point where the wrench
passes. Let point P in the x-z plane, where the wrench passes, has the coordinate
(x, 0, z). Consequently, the moment of the force system about P is
Note that this moment at P must be equal to the couple of the wrench pass-
ing through P , which is parallel to the resultant force. That is MP k R. By
comparing the ratio of their components, we have
100 100
=
100z 20 − 100z
and
10 − 100x = 0
As the result,
x = 0.1 m, z = 0.1 m
and substituting back into the moment equation, we have
MP = 10i + 10j Nm
Chapter 3
Equilibrium
3.1 Introduction 58
3.1 Introduction
Newton’s first law states about the condition of a particle to be at rest or to be
moving with a constant velocity. This is called the equilibrium condition. Later,
physicists have extended the equilibrium condition to a body.
When a body is in equilibrium, the resultant on the body is zero. And if the
resultant on a body is zero, the body is in equilibrium. Hence the necessary and
sufficient conditions for equilibrium of a body are
X
F=O (3.1)
and X
M=0 (3.2)
These equations reveal that it is the prerequisite to determine the resultant
force and moment acting on the body before the equilibrium condition can be
applied. The method of free body diagram, explained in the next section, is a
tool which helps in determining the correct resultant.
Figure 3.1 and 3.2 show common examples of modeling the action of forces
in two dimensional problems. Beware that these are not the FBDs. Only the
specific action forces, after the surroundings have been removed, are shown. Here
are some explanation of each particular case.
1. Force by a flexible cable is always the tension force. Weight of the cable
may be significant and hence make the cable sags.
Figure 3.1: Common action of forces in two dimensional analysis ([1], pp. 111)
Figure 3.2: Common action of forces in two dimensional analysis (continued) ([1],
pp. 112)
3. Ideally, roller, rocker, smooth guide, or slider eliminate the frictional force.
That is the supports cannot provide the resistance to motion in the tan-
gential direction.
4. Pin connection provides support force in any direction normal to the pin
axis. If the joint is not free to turn, a resisting couple may also be supported.
5. The built-in or fixed support of the beam is capable of supporting the axial
force, the shear force, and the bending moment.
7. Remote action force has the same overall effects on a rigid body as direct
contact force of equal magnitude and direction.
8. On the FBD, the force exerted on the body to be isolated by the body to
be removed is indicated.
9. The sense of the force exerted on the FBD by the removed bodies opposes
the movement which would occur if those bodies were removed.
10. If the correct sense cannot be known at first place, the sense of the scalar
component is arbitrarily assigned. Upon computation, a negative algebraic
sign indicates that the correct sense is opposite to that assigned.
3. All forces that act on the isolated body by the removed contacting and
attracting bodies are represented on the isolated body diagram. Forces
should be indicated by vector arrows, each with its magnitude, direction,
and sense. Consistency of the unknowns must be carried throughout the
calculation.
Example 3.1 ([1], Prob. 3/A,B) Complete the FBDs given in fig. 3.4 and 3.5.
Figure 3.3: Examples of drawing the free body diagram ([1], pp. 115)
are the resultant of the force system acting on the body of interest. Since the
object is in equilibrium, the resultant must be nulled and hence the zero value on
the right hand side. There are some things needed to be mentioned in applying
the equations.
1. The x-y coordinate system and the moment point O can be chosen arbi-
trarily.
2. Complete equilibrium in 2-D motion must satisfy all three equations. How-
ever, they are independent to each other. That is, equilibrium may only be
satisfied in some generalized coordinates.
At this point, the reader is armed with enough information to start practicing
the equilibrium problems. However it is helpful to know the characteristics of
some special cases of equilibrium problems. If we classify the problem according
to the force system, the following cases of equilibrium in two dimension, fig. 3.8,
may result. P
If the force system is collinear, only one equilibrium equation Fx = 0 is
effective. The other equations (perpendicular force and moment equations) are
satisfied automatically and hence contribute nothing in solving the problems.
For the concurrent force system, the moment at the point of concurrency is
always satisfied. Therefore only the force equilibrium equations are usable.
In the case of parallel forces, the force in the perpendicular direciton is null.
Hence only the force equilibrium equation in the parallel direction and the
moment equation are active. Note that in general case of force system, all three
equilibrium equations are valid.
Two-force member is a body under the action of two force only. The
equilibrium condition of the two-force member requires the acting forces be
equal, opposite, and collinear. Only one equation of force along its direction is
then effective. One common assumption of the two-force member is that the
weight of the member is negligible. Sometimes it is the trade-off between the
simplification in solving the equilibrium problem and the accuracy of the answer.
Figure 3.9 depicts some two-force members.
Figure 3.8: Special cases of equilibrium in two dimension ([1], pp. 122)
Equation 3.3, 3.4, and 3.5 are not the only set of equations for equilibrium.
Here we represent another two sets of equations which also guarantee the equi-
librium of an object. The first set of equilibrium conditions are
X
Fx = 0 (3.6)
X
MA = 0 (3.7)
X
MB = 0 (3.8)
where AB must P not be perpendicular to the x-direction. See fig. 3.11 for the
equation setup.
P Fx = 0 forces the reaction to be only in the perpendicular
P
direction. MA = 0 requires that reaction to pass through A, and so do MB =
0 constraints the reaction to pass through B. However all three requirements will
be satisfied only if the resultant is zero, otherwise the moment equations will
make the nonzero resultant break the force equation.
Another set of equilibrium conditions are
X
MA = 0 (3.9)
X
MB = 0 (3.10)
X
MC = 0 (3.11)
where A, B, andPC are not on P the same straight line. See fig. 3.12 for the
equation setup. MA = P 0 and MB = 0 force the reaction to pass through
A and B. The equation MC = 0 does the same thing. However all three
requirements will be satisfied only if the resultant is zero, otherwise the last
moment equation will make the nonzero resultant break the first two equations.
Before giving some examples, let us review the guideline in solving the equi-
Figure 3.13: Two Dimensional Constraints and Statical Determinacy ([1], pp.
125)
librium problems.
1. List known / unknown quantities. Then count the number of unknowns and
the number of available independent equations. If the number of unknowns
is greater than the number of equations, the problem cannot be solved solely
by the equilibrium conditions.
Example 3.2 ([1], Prob. 3/32) In a procedure to evaluate the strength of the
triceps muscle, a person pushes down on a load cell with the palm of his hand as
indicated in the figure. If the load-cell reading is 160 N, determine the vertical
tensile force F generated by the triceps muscle. The mass of the lower arm is 1.5
kg with mass center at G. State any assumptions.
the forces acting on it are the pushing force from the load-cell, the mass of
the arm, the tensile force by the triceps muscle, and the compressive force by
the humerus bone. The FBD of the system is shown in fig. 3.15. From the
figure, there are two unknowns. Hence two equations are required, which are
one moment and one force equations. First let us take the moment about O to
eliminate the compressive force.
P
[ MO = 0] −T × 25 − 1.5g × 150 + 160 × 300 = 0
T = 1832 N
C = 1977 N
Example 3.3 ([1], Prob. 3/39) The exercise machine consists of a lightweight
cart which is mounted on small rollers so that it is free to move along the inclined
ramp. Two cables are attached to the cart – one for each hand. If the hands
are together so that the cables are parallel and if each cable lies essentially in a
vertical plane, determine the force P which each hand must exert on its cable
in order to maintain an equilibrium position. The mass of the person is 70 kg,
the ramp angle is 15◦ , and the angle β is 18◦ . In addition, calculate the force R
Solution: Assume the cart is lightweight and the rail friction is negligible.
Then we select the machine and the man as our system and draw the FBD as
shown in fig. 3.17. It is seen that there are two unknowns which are solvable
using two force equations.
P
[ Fx′ = 0] 70g sin 15 − T cos 9 = 0 T = 179.9 N
P
[ Fy′ = 0] R − 70g cos 15 − T sin 9 = 0 R = 691 N
To determine the force exerted by each hand, we choose the pulley and cable
as our system of interest and draw its FBD, as depicted in fig. 3.18. Only the
force exerted by each hand, P, is the unknown. The force along the line of action
of T must be null which means
P
[ Fx′ = 0] −T + 4P cos 9 = 0 P = 45.5 N
Example 3.4 ([2], Prob. 3/29) A uniform ring of mass m and radius r carries
Solution: Write the FBD of the system which consists of the ring and the
eccentric mass. The force acting on the system are the weights, the friction
force, and the normal force. Here the friction is not negligible and indeed makes
the system be in equilibrium. We have three unknowns, which are F , N, and θ.
However, the problem only asks us for the angle θ.
First we should take the moment equation about O to eliminate the unknown
N.
mo gb sin θ
P
[ MO = 0] F r − mo gb sin θ = 0 F = r
Then we use the force equilibrium equation in the x′ -direction and substi-
tute the expression of F . Note that the normal force does not appear in the
equation.
h i
F − (mo + m) g sin α = 0 θ = asin rb 1 + mmo sin α
P
[ Fx′ = 0]
Example 3.5 ([1], Prob. 3/43) The hook wrench or pin spanner is used to
turn shafts and collars. If a moment of 80 Nm is required to turn the 200 mm
diameter collar about its center O under the action of the applied force P,
determine the contact force R on the smooth surface at A. Engagement of the
pin at B may be considered to occur at the periphery of the collar.
Solution: First select the system of the wrench and the collar as shown in
fig. 3.22. From the FBD, the applied force P must be
P
[ MO = 0] 80 − P × 0.375 = 0 P = 213.3 N
Now we must select only the spanner as our system which will reveal the
reaction force at A as shown in fig. 3.23. Note that this is a three-force member.
The moment equilibrium at B suppress the reaction R to appear in the equation.
P
[ MB = 0] NA × 0.1 sin 60 − P × (0.375 + 0.1 cos 60) = 0 NA = 1047 N
Example 3.6 ([1]) The small crane is mounted on one side of the bed of a
pickup truck. For the position θ = 40◦ , determine the magnitude of the force
supported by the pin at O and the oil pressure p against the 50 mm-diameter
piston of the hydraulic cylinder BC.
Solution: To determine such forces, we select the linkage ACO as our sys-
tem and draw its FBD. Note that this is a three-force member. Hence three
reaction forces must meet at a common point, which helps in drawing the correct
forces shown in fig. 3.25. However, before setting up the equilibrium equations,
let us consider the geometry at BCDO to determine the moment arm of the
hydraulic force. Referring to fig. 3.26,
360 + 340 sin 40 − 110 cos 40
α = atan = 56.2◦
340 cos 40 + 110 sin 40
Now we are ready to apply the equilibrium conditions to the system. Start
with the moment equation at O to eliminate unknown reaction at O.
P
[ MO = 0] 120g × (785 + 340) cos 40 − C × d = 0 C = 5063 N
Example 3.7 ([2], Prob. 3/48) The rubber-tired tractor shown has a mass
of 13.5 Mg with the center of mass at G and is used for pushing or pulling
heavy loads. Determine the load P which the tractor can pull at a constant
speed of 5 km/h up the 15-percent grade if the driving force exerted by the
ground on each of its four wheels is 80 percent of the normal force under that
wheel. Also find the total normal reaction NB under the rear pair of wheels at B.
100
NA + NB − 13500g × =0
P
[ Fy′ = 0] √
152 +1002
We see that there are three unknowns in two equations. Another equation
from the moment equilibrium must be solved simultaneously. Let us choose to
take the moment about A.
15
−13500g × √
152 +1002
× 0.825 = 0
Example 3.8 ([1], Prob. 3/59) Pulley A delivers a steady torque (moment) of
100 Nm to a pump through its shaft at C. The tension in the lower side of the
belt is 600 N. The driving motor B has a mass of 100 kg and rotates clockwise.
Determine the magnitude R of the force on the supporting pin at O.
of the pulley A. Since the pulley A delivers 100 Nm torque to a pump, there
must be the moment reaction resisting the clockwise rotation of it. See fig. 3.30.
The moment equilibrium is then used to determine the tension force, T .
P
[ MC = 0] (600 − T ) × 0.225 − 100 = 0
T = 155.6 N
Next we draw the FBD of the motor, fig. 3.31, of which the supporting
forces O and P are revealed. Determining three unknowns are relatively
straightforward.
P
[ MD = 0] Oy × 0.25 − 600 × (0.2 − 0.075) − 100g × 0.125
Oy = 906 N
P
[ Fx = 0] T cos 30 + 600 − Ox = 0 Ox = 734.7 N
p
O= Ox2 + Oy2 = 1.17 kN
P
[ Fy = 0] T sin 30 − 100g − P + Oy = 0 P = 2.8 N
Note the tension in the rod at the left leg results from the compressed
spring to resist the rotation of the motor stator.
Example 3.9 ([2], Prob. 3/52) When setting the anchor so that it will dig into
the sandy bottom, the engine of the 40 Mg cruiser with the center of gravity at
G is run in reverse to produce a horizontal thrust T of 2 kN. If the anchor chain
makes an angle of 60◦ with the horizontal, determine the forward shift b of the
center of buoyancy from its position when the boat is floating free. The center of
buoyancy is the point through which the resultant of the buoyant force passes.
Solution: When the boat is free floating (no thrust or tension), the buoy-
ancy force is equal to the weight and acting at the C.G. At the times the boat
moves backward, the anchor chain is pulled against the thrust force. This
changes the buoyancy force, both in magnitude and point of application, to
maintain the equilibrium. Figure 3.33 is the FBD of this system. Apply the
equilibrium equations to solve for three unknowns; A, B, and b.
P
[ Fx = 0] A cos 60 − 2000 = 0 A = 4 kN
P
[ Fy = 0] B − 40000g − A sin 60 = 0 B = 395, 864 N
P
[ MA = 0] 40000g × 8 − 2000 × 3 − Bx = 0 x = 7.915 m
b = 8 − x = 85.2 mm
Example 3.10 ([2], Prob. 3/63) A special jig for turning large concrete pipe
sections (shown dotted) consists of an 80 Mg sector mounted on a line of rollers
at A and a line of rollers at B. One of the rollers at B is a gear which meshes
with a ring of gear teeth on the sector so as to turn the sector about its geometric
center O. When α = 0, a counterclockwise torque of 2460 Nm must be applied
to the gear at B to keep the assembly from rotating. When α = 30◦ , a clockwise
torque of 4680 Nm is required to prevent rotation. Locate the mass center G of
the jig by calculating r and θ.
Solution: The problem gives two different postures that are in equilibrium
and ask for the locating the C.G. by the parameters r and θ. Therefore two
equations coming from each posture of equilibrium must be constituted and
solved. For the case when α = 0, a CCW torque is applied to the gear B.
Therefore the reaction force at the meshing teeths with the jig must be as shown
in fig. 3.35. For the case of α = 30◦ , a CW torque is applied instead. Hence
the reaction force at the meshing teeths is as shown in fig. 3.36. In each case,
the tangential component of the reaction force is determined by the moment
equilibrium equation:
α = 0◦ : 2460 − F1 × 0.24 = 0
P
[ MB = 0] F1 = 10250 N
The reaction is imparted to the jig, shown in fig. 3.37, along with the one
from the rollers A. To avoid determining the unnecessary unknowns, NA and
NB , we should use the moment equation taken at the center point O. For each
case, we have
r = 367 mm θ = 79.8◦
1. The x-y-z coordinate system and the moment point O can be chosen arbi-
trarily.
2. Complete equilibrium in 3-D motion must satisfy all six equations. How-
ever, they are independent to each other. That is, equilibrium may only be
satisfied in some generalized coordinates.
Figure 3.38 shows common examples of modeling the action of forces in three
dimensional problems. Beware these are not the FBDs. Only the specific action
forces, after the surroundings have been removed, are shown. Here are some
explanation of each particular case.
1. Reaction force of the member in contact with the smooth surface, or ball-
supported member, is normal to the surface and directed toward the mem-
ber. It is usually denoted by N.
Figure 3.38: Common action of forces in three dimensional analysis ([1], pp. 147)
2. Reaction force of the member in contact with the rough surface has both
the normal and tangential components. To determine these forces, we first
seek for the contact plane which is tangent to both contact surfaces. Normal
force, N, acts in the direction normal to the plane while tangential force,
F, lies in the plane. Tangential force is commonly recognized as the friction
force.
4. The ball and socket joint constraints the point on each linkage to always
be together. This requires the joint to support a force R.
6. The thrust bearing support prevents the shaft from moving and rotating in
all direction except the rotational motion about the shaft axis. Hence three
dimensional reaction forces, Rx , Ry , and Rz , must be supplied. In addition,
the resisting moment in x and z- direction must be provided as well.
At this point, the reader is armed with enough information to start practicing
the equilibrium problems. However it is helpful to know the characteristics of
some special cases of equilibrium problems. If we classify the problem according
to the force system, the following cases of equilibrium in three dimension, fig. 3.39,
may result.
For the concurrent force system, the moment at the point of concurrency
is always satisfied. Therefore only the force equilibrium equations are usable.
If the force system has all forces being concurrent with a line, i.e. all forces
intersecting a common line, then the moment of the forces about that line is
automatically null. Therefore only five equilibrium equations are effective. In
the case of parallel forces, the force in the perpendicular direction is null. Hence
only the force equilibrium equation in the parallel direction and the moment
equations about the lines in the plane perpendicular to the parallel forces are
active. Note that in general case of force system, all six equilibrium equations
are valid.
Figure 3.39: Special cases of equilibrium in three dimension ([1], pp. 148)
then we cannot solve it. This is because the system has more constraints than
necessary to maintain the equilibrium. This is called statically indeterminate
system. Extra equations obtained from force-deformation material properties
must also be applied to solve for the redundant constraints. If the number of
unknown is equal or less than the number of equilibrium equations, the system
is statically determinate. With some wrong installation of the supports, the
number of active constraints may be reduced unintentionally. See case (b) and
(c) of fig. 3.40.
Example 3.11 ([2], Prob. 3/71) The light right angle boom which supports
the 400 kg cylinder is supported by three cables and a ball-and-socket joint at O
attached to the vertical x-y surface. Determine the reactions at O and the cable
tensions.
Solution: The FBD of this problem is shown in fig. 3.42. There are six
unknowns; three reaction force components at O and three cable tensions.
We can simply set up three force equilibrium equations and three moment
equilibrium equations about any point. However solving the resulting equations
may be intractable. A little thought before writing down the solution can be
helpful. For this problem, the moment equilibrium about line OB suppresses
five unknowns to show up. Only TAC will appear in the equation, which is easy
to solve. Examining the FBD in this manner helps avoid solving simultaneous
equations unnecessarily.
First, let us determine the pertaining unit vectors:
nAC = −0.408i + 0.408j − 0.816k, nBD = 0.707j − 0.707k
nBE = −k, nOE = i
nOB = 0.6i + 0.8k, nOD = 0.6i + 0.8j
TAC = 4808.8 N
[(2k) × TAC nAC + (0.75i + 2k) × (−400gj) + (1.5i) × TBE nBE ] · nOD = 0
TBE = 654 N
Figure 3.40: Three Dimensional Constraints and Statical Determinacy ([1], pp.
149)
The last cable tension TBD is found from taking the moment about line OE.
P
[ MOE = 0]
[(2j) × TBD nBD + (0.75i + 2k) × (−400gj) + (2k) × TAC nAC ] · nOE = 0
TBD = 2775.1 N
Ox = 1962 N, Oy = 0 N, Oz = 6540 N
P
[ F = 0]
Example 3.12 ([1], Prob. 3/67) The 600 kg industrial door is a uniform
rectangular panel which rolls along the fixed rail D on its hanger-mounted
wheels A and B. The door is maintained in a vertical plane by the floor-mounted
guide roller C, which bears against the bottom edge. For the position shown
compute the horizontal side thrust on each of the wheels A and B, which must
be accounted for in the design of the brackets.
Solution: The FBD of the door including the hanger A and B is shown in
fig. 3.44. Due to the offset of the door’s weight and the normal forces at the
wheels, the normal force at the guide rollor C points in the direction shown to
counter for the induced moment. Hence the side thrusts on the wheels must
point in the opposite direction to balance the force.
P
[ MAB = 0] 600g × 0.15 − NC × 3 = 0, NC = 294.3 N
P
[ MAz = 0] NC × 0.6 − Bx × 3 = 0, Bx = 58.86 N
P
[ Fx = 0] Ax + Bx − NC = 0, Ax = 235.44 N
Example 3.13 ([1], Prob. 3/73) The smooth homogeneous sphere rests in the
120◦ groove and bears against the end plate which is normal to the direction
of the groove. Determine the angle θ, measured from the horizontal, for which
the reaction on each side of the groove equals the force supported by the end plate.
Solution: First the FBD will be drawn according to the given condition
that the reaction on each side of the groove must equal the force supported by
the end plate. However, it is rather difficult to draw and visualize the FBD of the
sphere in three dimensions. Therefore the FBD will be drawn in two orthogonal
views; along and perpendicular to the V-groove. Forces acting on the sphere will
then be the projected components onto that view. Also the coordinate frame
x-y-z is set at the center of the sphere for convenience. See fig. 3.46.
Now the equilibrium conditions are ready to be applied under the condition
that N1 = N2 = Nr = N.
P
[ Fy = 0] N1 = N2 = N
P
[ Fz = 0] mg cos θ = 2N cos 30
P
[ Fx = 0] Nr = mg sin θ
Example 3.14 ([1], Prob. 3/75) The mass center of the 30 kg door is in the
center of the panel. If the weight of the door is supported entirely by the lower
hinge A, calculate the magnitude of the total force supported by the hinge at B.
P P
[ Fx = 0] MAy = 0 30g × 0.36 − Bx × 1.5 = 0, Bx = Ax = 70.6 N
P P
[ Fy = 0] [ MAx = 0] By × 1.5 − 30g × 0.9 = 0, By = Ay = 176.6 N
q
B = Bx2 + By2 = 190.2 N
Example 3.15 ([1], Prob. 3/82) One of the three landing pads for the Mars
Viking lander is shown in the figure with its approximate dimensions. The mass
of the lander is 600 kg. Compute the force in each leg when the lander is resting
on a horizontal surface on Mars. Assume equal support by the pads and consult
table D/2 in Appendix D as needed.
Solution: From the problem statement, it can be deduced that each land-
ing pad supports 200 kg force. However the gravitational constant on Mars
is g = 3.73m/s2 . The leg with the landing pad is sectioned out and its FBD
is drawn. Forces in each strut are then revealed. Note that the forces are all
concurrent at the center of ball and socket. Hence the available equilibrium
conditions are reduced to three only.
First the relating unit vectors are determined:
FDC = 1049.1 N
Example 3.16 ([2], Prob. 3/88) The uniform 15 kg plate is welded to the
vertical shaft, which is supported by bearings A and B. Calculate the magnitude
of the force supported by bearing B during application of the 120 Nm couple
to the shaft. The cable from C to D prevents the plate and shaft from turning,
and the weight of the assembly is carried entirely by bearing A.
Solution: The plate and the attached vertical shaft are selected as the
system and its FBD is drawn. By the statement that the weight of the assembly
is carried entirely by bearing A, bearing B will experience only the radial forces.
Since the only vertical force is the assembly’s weight, the vertical component
reaction at the bearing A must be equal to the weight. The convenient coordinate
frame is defined and now it is ready to set up the equilibrium conditions.
First, take moment about the z-axis so all unknowns except the tension T
are eliminated.
P
[ MOz = 0] 120 + 0.6i × T nDC · k = 0, T = 632.9 N
Then, take moment about the line in y-direction passing A so the reac-
tion forces at A and By are eliminated.
P
MAy = 0 −Bx × 0.2 + 15g × 0.3 + Tx × 0.68 = 0, Bx = 2265 N
Example 3.17 ([2], Prob. 3/95) The uniform 900 × 1200 mm trap door has
a mass of 200 kg and is propped open by the light strut AB at the angle
θ = atan 43 . Calculate the compression FB in the strut and the force supported
by the hinge D normal to the hinge axis. Assume that the hinges act at the
extreme ends of the lower edge.
Solution: The trap door is isolated and its FBD is drawn. There are the
reactions at C and D, the strut force, and its weight. Reactions at D are
decomposed along the coordinate frame axes. It is unnecessary to do so for the
reactions at C, however.
Since the strut lies along the obligue line, it is helpful to determine the direc-
tional unit vector.
Example 3.18 ([1], Prob. 3/95) The uniform rectangular panel ABCD has a
mass of 40 kg and is hinged at its corner A and B to the fixed vertical surface.
A wire from E to D keeps edges BC and AD horizontal. Hinge A can support
thrust along the hinge axis AB, whereas hinge B supports force normal to the
hinge axis only. Find the tension T in the wire and the magnitude B of the force
supported by hinge B.
Solution: The FBD of the panel ABCD is shown in fig. 3.56. Note that
there is no reaction component along the hinge axis at B. Assigning the
coordinate frame x-y-z as shown,
Take moment about the x-axis so only the tension force remains:
P
[ MAx = 0] [0.6j × 40g (− cos 30k − sin 30i)] · i + (1.2j × TDE nDE ) · i = 0
TDE = 278.55 N
Next, take moment about point A along the y-axis. Only Bz then ap-
pears in the equation.
P
MAy = 0 [1.2i × 40g (− cos 30k − sin 30i)] · j − 2.4Bz = 0
Bz = 169.9 N
Finally, take moment about the vertical line AE. Reaction forces at A,
the tension force, the panel’s weight, and Bz do not contribute any moment.
Therefore
By = 0 N ⇒ Bn = 169.9 N
P
[ MAE = 0]
Example 3.19 ([2], Prob. 3/98) Under the action of the 40 Nm torque (couple)
applied to the vertical shaft, the restraining cable AC limits the rotation of
the arm OA and attached shaft to an angle of 60◦ measured from the y-axis.
The collar D fastened to the shaft prevents downward motion of the shaft in
its bearing. Calculate the bending moment M, the compression P , and the
shear force V in the shaft at section B. (note: Bending moment, expressed as a
vector, is normal to the shaft axis, and shear force is also normal to the shaft axis.)
Solution: The shaft is sectioned at B (fig. 3.58) so that the reaction forces and
moments are exposed. Reactions are the compressive force, the shear force, and
the bending moment. However the shaft is free to rotate about its axis. Hence
no supporting moment in this direction.
Take moment about the z-axis, so the tension force can be determined.
P
[ Mz = 0] 40 + (0.18j × TAC nAC ) · k = 0 ⇒ TAC = 419.3 N
The shear force is the remaining of the tension force after subtracting the
compressive force.
q q
2
V = Vx2 + Vy2 = TAC − P 2 = 273.5 N
P
[ MB = 0] MBx i + MBy j + 40 + (−0.09k + 0.18j) × TAC nAC = 0
q
MBx = 42.87 Nm, MBy = 20.0 Nm Mb = MB2 x + MB2 y = 47.3 Nm
Example 3.20 ([1]) Determine the turning force at the knob DE and the
normal forces at each leg.
To determine the normal forces at each leg, consider the FBD of the whole
assembly shown on the right of fig. 3.60. Note of the friction forces which are
not negligible. Taking moment at B along the x-axis,
P 200 × 0.2 − 100 sin 15 × 0.35 − P sin 30 × 0.075 − NC × 0.5 =
[ MBx = 0]
0, NC = 58.13 N
Taking moment at B along the y-axis,
P
MBy = 0
Example 3.21 ([1], Prob. 3/114) The drum and shaft are welded together
and have a mass of 50 kg with mass center at G. The shaft is subjected to a
torque (couple) of 120 Nm, and the drum is prevented from rotating by the cord
wrapped securely around it and attached to point C. Calculate the magnitudes
of the forces supported by bearings A and B.
Solution: Select the drum and shaft unit and draw the FBD. Each bear-
ing supports the shaft by the radial force which is projected into the x- and z-
direction as shown in fig. 3.62.
First, take moment about B along the y-axis.
P
MBy = 0 T × 0.15 − 120 = 0, T = 800 N
Then, take moment about B along the z- and x-axes to obtain the bear-
ing forces at A.
P
[ MBz = 0] T cos 66.87 × 0.36 − Ax × 0.7 = 0, Ax = 161.6 N
P
[ MBx = 0] 50g × 0.3 + T sin 66.87 × 0.36 − Az × 0.7 = 0, Az = 588.6 N
Finally, use the force equilibrium conditions to determine the bearing forces at B.
P
[ Fx = 0] Ax + Bx − T cos 66.87 = 0, Bx = 152.6 N
P
[ Fz = 0] Az + Bz − 50g − T sin 66.87 = 0, Bz = 637.6 N
Therefore,
p p
A= A2x + A2z = 610.4 N, B = Bx2 + Bz2 = 655.6 N
Example 3.22 ([1], Prob. 3/115) Determine the reaction force and moment at
the double U-joint support O.
P
[ MAB = 0] Oz = 0 N
To determine the cable forces, take the moment about O normal to the y-axis
because the unknown moment M always lies along the y-axis.
P
[ Mz = 0] (1.8i × TBC nBC ) · k + (2.1j × TAD nAD ) · k = 0
P
[ Mx = 0] (2.1j × TBC nBC ) · i + (2.1j × TAD nAD ) · i − 50g × 2.1 = 0
Now the moment equilibrium about y-axis will give the supporting moment
at the double U-joint.
P
[ My = 0] M + 50g × x + (1.5i × TBC nBC ) · j = 0
M = 365.66 − 490.5x
Finally, the force equilibrium along the x- and y- axis will help determine the
remaining supporting forces.
P
[ Fy = 0] Oy + TBC nBC · j + TAD nAD · j = 0, Oy = 1429 N
P
[ Fx = 0] Ox + TBC nBC · i + TAD nAD · i = 0, Ox = 410 N
q
O = Ox2 + Oy2 + Oz2 = 1487 N
Chapter 4
Friction
4.1 Introduction
When two objects are in contact, the forces of action and reaction between con-
tacting surfaces are developed. To the mutual contacting surface, these forces
have their components both in the tangential and normal directions. Force com-
ponent in the tangential direction is known as the friction force. Whenever a
tendency exists for one contacting surface to slide along another sur-
face, the developed friction force is always in the direction opposing
this tendency.
In some systems, friction is undesirable because it just plainly changes the
system characteristics from the required behavior. In particular, where the slid-
ing motion between parts occurs, the developed friction force results in a loss
of energy. However, in many cases, friction instead functions the systems. For
example, many mechanisms employ friction as an extra force to retain their equi-
librium states.
F = N tan α (4.1)
α, called the friction angle, is the angle the resulting reaction force deviated
from the normal direction. This angle can be determined approximately by the
experiment. The value depends largely on the materials of two contacting surface.
Fig. 4.1(c) explains why the parameter is at best obtained as an estimated value.
Local geometry of the contact point determines each infinitesimal reaction force
which adds up vectorially to give the gross reaction force R. A multitude of
parameters dictate the value of the macro friction angle α, which therefore is the
source of uncertainties.
Relationship of the friction force F and the applied force P is plotted in
fig. ??(d). There are three regimes of the friction force development. In the
beginning of the application, the friction force equally catch up the applied force
as shown by the 45◦ straight line dictating F = P . During this phase, the
block does not move at all. This can be understood by determining the resulting
force acting along the tangential direction. The force is zero which indicates the
equilibrium status and hence no motion of the block.
The first regime lasts until the applied force P reaches
Fmax = µs N (4.2)
for which the graph enters the second regime. µs is the static coefficient of friction.
By comparing eq. 4.1 and 4.2,
µs = tan αs (4.3)
Hence the coefficient of friction depends on the material of the mating surfaces.
At this point, the block is about to move, called the impending motion, be-
cause the current applied force is the largest force resistable by the physically
largest generated friction force Fmax , governed by eq. 4.2. From the equation,
the maximum friction force depends on the normal force and the static coeffi-
cient of friction. The implication of the impending motion is the validity of the
equilibrium condition and the friction force reaches the maximum value.
If the applied force is further increased, the block is no longer in the equilib-
rium. That is the frictional force is now not enough to hold the block at rest.
From the unbalanced resulting force, the block will then start moving in the di-
rection of the applied force. Moreover, the magnitude of the friction force itself
is known to be dropped off from the maximum value Fmax . The new value of the
friction in the third regime is
Fk = µk N (4.4)
where µk is the kinetic coefficient of friction determined experimentally. Its value
is usually less than µs and also depends on the material types of the contact
surfaces. Practically, the friction force decreases as the relative velocity of the
mating surfaces increases.
In summary, there are three regions of the friction force development in tran-
sitioning of the object from rest to motion. These are
1. No Motion is the region up to the point of slippage or impending motion.
Friction force is determined by the equations of equilibrium because the
system is in equilibrium. When the motion is not impending, F < Fmax .
2. Impending Motion is the moment where the body is on the verge of slipping.
Static friction force reaches the maximum value. For a given pair of mating
surfaces, F = Fmax = µs N.
3. Motion The body starts moving in the direction of the applied force. Here,
friction force drops to a lower value called kinetic friction F = µk N. It will
drop further with higher velocity.
Unknown status of the problem This is the difficult case. The prob-
lem is tackled by first assuming the system is in static equilibrium and using the
equilibrium condition to determine the required friction force F . The result is
then investigated to conclude the validity of the equilibrium assumption:
1. F < µs N indicates that the friction force for the assumed equilibrium
condition can be provided and so the body is in static equilibrium as firstly
assumed.
3. F > µs N This is impossible because the surface cannot support more fric-
tion than µs N. Consequently, the equilibrium assumption is invalid and
motion occurs instead. This makes the friction force be directed by the
kinetic friction force; F = µk N. And the equilibrium conditions are no
longer held. The motion is accelerated.
Using the above guidelines, the following representative problems are now
presented and solved.
Example 4.1 ([1], SP. 6/1) Determine the maximum angle θ which the ad-
justable incline may have with the horizontal before the block of mass m begins
to slip. The coefficient of static friction between the block and the inclined
surface is µs .
Solution: From the given statement, the block is on the verge of slipping.
Therefore the friction force is the maximum value F = µs N acting upward the
incline. According to the free body diagram in fig. 4.4, the following equilibrium
equations can be formulated.
P
[ Fy = 0] N − mg cos θ = 0
P
[ Fx = 0] µs N − mg sin θ = 0
µs = tan θ or θ = tan−1 µs
tan−1 µs to the normal line when the motion impends, and noting that from the
geometry α = θ, it then can be concluded that
θ = tan−1 µs
Example 4.2 ([1], SP. 6/2) Determine the range of values which the mass mo
may have so that the 100 kg block shown in the figure will neither start moving
up the plane nor slip down the plane. The coefficient of static friction for the
contact surface is 0.30.
N = 922 N
mo = 62.4 kg
mo = 6.0 kg
6.0 ≤ mo ≤ 62.4 kg
Example 4.3 ([1], SP. 6/3) Determine the magnitude and direction of the
friction force acting on the 100 kg block shown if, first, P = 500 N and, second,
P = 100 N. The coefficient of static friction is 0.20, and the coefficient of kinetic
friction is 0.17. The force are applied with the block initially at rest.
Solution: After reading the problem statement, the motion status of the
block cannot be deduced. The problem falls in the unknown status case and the
N = 1092.85 N
P
[ Fx = 0] 500 cos 20 − F − 100g sin 20 = 0
F = 134.3 N
N = 956.04 N
P
[ Fx = 0] F + 100 cos 20 − 100g sin 20 = 0
F = 241.55 N
The friction force then becomes the kinetic friction force in the upward direction
of which its magnitude is
Example 4.4 ([1], SP. 6/4) The homogeneous rectangular block of mass m,
width b, and height H is placed on the horizontal surface and subject to a
horizontal force P which moves the block along the surface with a constant
velocity. The coefficient of kinetic friction between the block and the surface is
µk . Determine
(a) the greatest value that h may have so that the block will slide
without tipping over and
(b) the location of a point C on the bottom face of the block through which
the resultant of the friction and normal force acts if h = H/2.
Solution: In the first question, as the point of application of the applied
force P increases, the block will be more likely about to tip over. When the
block is about to tip over, the resultant supporting force from the ground must
be acting through the last contact point at the far corner of the block, labeled as
point A in fig. 4.10. For this problem, the block is being moved with a constant
velocity. Therefore the developed friction is the kinetic friction force. From the
friction cone viewpoint, the resultant force R must be such that its direction is
forming the angle
θ = tan−1 µk
with the normal direction. Since the block is in equilibrium (moving with constant
velocity), all three forces acting on the block must meet at the common point
denoted as B in the figure. The location of this point will dictate the height h of
the applied force P . Because the ground reaction force must pass through point
A and B, the following geometric relationship can be written:
b/2h = tan θ = µk
h = b/2µk
For the second question, the height of the point of application is fixed to be
h = H/2. Nevertheless, the system is still in equilibrium and hence three forces
acting on the block must also meet at the common point. In this case the point is
known; it is the center of gravity point G. Since the ground reaction force must
pass through this point and its direction must be such that
θ = tan−1 µk
the location of the point C on the bottom face through which the resultant
reaction force acts must satisfy the following relationship:
x/ (H/2) = tan θ = µk
In other words, the location of the point C from the middle vertical line of the
block is
x = µk H/2
Example 4.5 ([1], SP. 6/5) The three flat blocks are positioned on the 30◦
incline as shown, and a force P parallel to the incline is applied to the middle
block. The upper block is prevented from moving by a wire which attaches it
to the fixed support. The coefficient of static friction for each of the three pairs
of mating surfaces is shown. Determine the maximum value which P may have
before any slipping takes place.
The maximum friction forces possible developed at each mating surface there-
fore are
F1max = µs N1 = 76.46 N
F2max = µs N2 = 271.86 N
F3max = µs N3 = 458.8 N
The normal force will next be used to determine the friction forces right before
any slippage as follow.
which is not possible. Therefore the 40 kg-block cannot stay still and the
requirement of no slippage cannot be satisfied. Hence this case will not happen.
which is possible. Therefore the 50 and 40 kg-block will not slip relative
to each other and be at rest. To determine the maximum value of the applied
force P , apply the equilibrium condition at the middle block.
P
[ Fx = 0] P − F1 − F2 + 50g sin 30 = 0
P = 93.8 N
Example 4.6 ([1], Prob. 6/6) The light bar is used to support the 50 kg block
in its vertical guides. If the coefficient of static friction is 0.30 at the upper and
0.40 at the lower end of the bar, find the friction force acting at each end for
x = 75 mm. Also find the maximum value of x for which the bar will not slip.
Solution: The most distinct aspect of the system is that the light bar is a
two-force member. Therefore, if it is in equilibrium, the reaction forces acted
at both ends must be equal, opposite, and acting along the bar’s longtitudinal
direction, depicted in fig. 4.14. Decomposing the reaction force into the normal
and friction force, the normal force must be such that it supports the weight of
the 50 kg-block.
P
[ Fy = 0] N − 50g = 0, N = 490.5 N
Because the required situation of this problem is static equilibrium, the de-
veloped friction must not exceed the maximum value in addition to the two-force
member equilibrium condition. The alternative view to this constraint is that
the direction of the reaction force must reside in the static friction cone. Hence
the angle that the reaction force on each end made with the normal (vertical)
direction cannot exceed
φsA = tan−1 µsA = 21.8◦
φsB = tan−1 µsB = 16.7◦
Consider the case where the lower end of the bar is displaced by x = 75 mm.
Direction of the reaction force is along the propped bar. By geometry, the angle
the reaction force made with the normal direction is
The reaction force R is inside the static friction cone and therefore the system
is in equilibrium. The friction force can then readily be determined as
F = N tan θ = 126.6 N
To find the maximum value of x for which the bar will not yet slip, one of the
reaction forces must be on the static friction cone surface. That reaction force
is at the upper end where the coefficient of friction is lower. The lower reaction
force will still be inside the static friction cone. At B, the direction of the reaction
force, and so the bar, must then be such that
Example 4.7 ([1]) Find the tension in the cable and force P that makes the 15
kg lower block
F1max + 8g sin 20 − T = 0, T = 49 N
b) P pushing up There are two possible situations. Either the 15 & 8 kg-
blocks impend to slide upward together or the 15 kg-block impends to slide
upward alone. If both blocks are about to go together, the cable will get slack.
This implies T = 0 N. Using the equilibrium condition on the 8 kg-block, this
requires the resisting friction force of
P
[ Ft = 0] 8g sin 20 − F1 = 0, F1 = 26.84 N
to slide upward alone will be realized instead. Using the equilibrium conditions
along the tangential direction on each block, the applied force P and the cable
tension T can be determined.
P
[ Ft = 0] −P + F1max + F2max + 15g sin 20 = 0, P = 157.3 N
Example 4.8 ([2], Prob. 6/29) The uniform slender rod of mass m and length
L is initially at rest in a centered horizontal position on the fixed circular surface
of radius R = 0.6L. If a force P normal to the bar is gradually applied to its
end until the bar begins to slip at the angle θ = 20◦ , determine the coefficient of
static friction.
Solution: The status of the bar is impending motion. From the no slip
condition, the distance between the current and the initial contact points of the
bar must be equal to the length of the curve traced on the surface, which is
π
[a = θr] a = 20 × 180
R = πR/9
After drawing the free body diagram and recognizing the angle α the reaction
force made with the normal direction is also the angle of the static friction cone,
as shown in fig. 4.18, the coefficient of static friction can be determined as
L/2 − πR/9
µs = tan α = = 0.211
L/ (2 tan 20)
Example 4.9 ([2], Prob. 6/32) The three identical rollers are stacked on a
horizontal surface as shown. If the coefficient of static friction µs is the same
for all pairs of contacting surfaces, find the minimum value of µs for which the
rollers will not slip.
Solution: After careful consideration of the sketch in fig. 4.19, if the rollers start
slipping, the lower rollers tend to roll out at the upper contacts while they tend
to slide out at the lower contacts. The top most roller will then start falling
down. This described motion is shown in fig. 4.20. The slipping will occur if
one or more contacts of the lower rollers impend to slip. To understand the
phenomenon, the free body diagram of the lower left roller is drawn as depicted
in fig. 4.20. Applying the moment equilibrium at point O,
P
[ MO = 0] FA = FB
From the figure, it can be concluded that NA < NB . Therefore FAmax <
FBmax . Consequently, FA will reach the limit value before FB . This implies
slipping will occur first at A and so FA = FAmax . FB , which may be less than
FBmax , will be determined by the equilibrium equation.
The roller is recognized as a three-force member. Therefore the line of action
of these three forces must intersect at a unique point for the system to be in
equilibrium. The only possible meeting point is B. In addition, since the reaction
force at A reaches the limit, RA must be making the angle tan−1 µs with NA . From
the geometry of the triangle OAB sketched in fig. ??, it can be stated that
tan−1 µs = 15◦
Example 4.10 ([2], Prob. 6/41) The industrial truck is used to move the solid
1200 kg roll of paper up the 30◦ incline. If the coefficients of static and kinetic
friction between the roll and the vertical barrier of the truck and between the
roll and the incline are both 0.40, compute the required tractive force P between
the tires of the truck and the horizontal surface.
Solution: The moving paper roll is making contacts with the vertical bar-
rier of the truck and the incline. This problem is difficult because it is not clear
from the statement how the roll is moving. After drawing the free body diagram
of the roll shown in fig. 4.22, three possibilities are listed.
1. A and B both slip
2. Only B slips
3. Only A slips
Each case has been investigated with the assumption of dynamical equilib-
rium. After the calculation, only the last case — i.e. the roller slips at A —
is viable. For this case, because the roll slips at A, the kinetic friction force is
determined as
FA = 0.4NA
Applying the equilibrium conditions to the roller, other unknown forces can
be determined as follow:
P
[ MO = 0] FB = FA = 0.4NA
P
[ Fx = 0] NA − FB cos 30 − NB sin 30 = 0, NB = 1.307NA
P
[ Fy = 0] −0.4NA − 1200g − FB sin 30 + NB cos 30 = 0
P = 22.1 kN →
Chapter 5
Distributed Forces
5.1 Introduction
In reality, forces are applied over a region of which often its area is negligible.
When forces are applied over a region whose dimensions are not negligible com-
pared with other pertinent dimensions, we must account for the actual manner
in which the force is distributed by summing up the effects of the distributed
force over the region. For this purpose, we need to know the intensity of the force
at any location and we will use the integration to determine their total effect.
Figure 5.1 are common engineering examples of the distributed forces.
Figure 5.1: Some examples of the distributed forces ([1], pp. 237)
Most of the objects subject to the earth gravity can be safely assumed that the
intensity of the earth’s force field over the body is uniform. The other reasonable
assumption is that the field of force due to the gravitational attraction is parallel.
Applying these assumptions to the principle of moment simplifies the derivation
of the C.G. because the cross product degenerates to a simple multiplication.
Line We assume the cross sectional area is constant over the length. The
centroid of the line is determined as
R
rc dL
rC = R (5.4)
dL
Note that the centroid may not always be on the line. See fig. 5.3.
Area We assume the thickness is constant over the entire area. The cen-
troid of the area is determined as
R
rc dA
rC = R (5.5)
dA
R
rc dA is called the first moment of area. Note that the centroid may not always
be on the surface. See fig. 5.4.
or ZZ
A= dxdy
or ZZZ
V = dxdydz
Figure 5.5: Select the differential element, of which the volume is known, to be
as large as possible ([1], pp. 241)
On the other hand, if the boundary is explained using the polar coordinates,
as in fig. 5.8b, the area should be calculated as
Z θo
1 2
A= r dθ
0 2
Figure 5.6: If possible, divide the object in the manner that can be integrated in
one continuous operation to cover the whole body ([1], pp. 241)
Example 5.1 ([1], SP. 5/3) Locate the centroid of the area of a circular sector
with respect to its vertex.
Figure 5.7: The higher order terms can safely be omitted from the integrand ([1],
pp. 241)
Figure 5.8: Select the coordinates that suit to the shape and boundaries of the
object ([1], pp. 242)
R Rr Rr ro sin α
XA = xc dA X 0
2αro dro = 0 α
2αro dro
2r sin α
X=
3α
Infinitesimal sector Here the area is divided into infinitesimal sectors as shown
in the right of fig. 5.11. Again, from the observation, Y = 0. Only the x-
coordinate of the centroid, X, will be determined. The infinitesimal area of a
sector at the angle θ, as highlighted in red, is
1
dA = r (rdθ)
2
The centroid of this ring can be determined as a separted problem of finding the
centroid of the infinitesimal sector of radius r located at the angle θ. As a result,
Figure 5.9: Divide the object in the manner that the centroid be expressed easily
in terms of the selected coordinates ([1], pp. 242)
Example 5.2 ([2], Prob. 5/23) Locate the centroid of the area shown in the
figure by direct integration.
Solution: It can be observed that the given area is symmetric about the
line y = a − x. Therefore the centroid will lie on this line. This implies
Y =a−X
The shaded area can be viewed as the summation of the horizontal strips.
Each strip is built from subtracting the horizontal strip of the length a with
the differential sector located at the corresponding angle θ. These strips will be
evaluated through the whole range of θ varying from 0 to π/2.
Each infinitesimal area is then
1
dA = ady − a (adθ)
2
Transform the y- to θ- coordinate to perform the integration under the same
variable. We have the following relationship:
y = a − a cos θ
dy = a sin θdθ
The x-coordinate centroid of these differential elements can be observed as
xc1 = a/2
2
xc2 = a sin θ
3
for the horizontal and the sector strips, respectively. Now it is ready to determine
the x-coordinate centroid of this area.
R
XA = xc dA
Z π/2 Z π/2 Z π/2
2 1 a 2 2 1
X a sin θ − dθ = × a sin θdθ − a sin θ × a2 dθ
0 2 0 2 0 3 2
2a
X=
3 (4 − π)
Recall the centroid location on the line y = a − x. Hence, the y-coordinate
centroid of the area is
(10 − 3π) a
Y =a−X =
3 (4 − π)
Example 5.3 ([2], Prob. 5/35) Determine the coordinates of the centroid of the
volume obtained by revolving the shaded area about the z-axis through 90◦ angle.
Solution: This revolved volume is symmetric about the vertical plane making
45◦ with the x-axis. Therefore X = Y .
The given volume may be dissected as shown in fig. 5.15. The y-coordinate
may be chosen as the integration variable varying from 0 to a. The circular edge
in the yz plane has the following relation:
y 2 + (z − a)2 = a2
R
ZV = zc dV
a Z a p
π a − a2 − y 2 π
Z p p
Z 2 2
a − a − y ydy = × 2 2
a − a − y ydy
0 2 0 2 2
a
Z=
4
Example 5.4 ([2], Prob. 5/40) Locate the center of mass G of the steel half ring.
Solution: It is readily seen that the planes of symmetry of the steel half
ring is the middle vertical plane and the zero horizontal plane. Therefore the
centroid lies on the intersection line of these two planes. Consequently, only the
distance r is needed to be determined. To do so, the differential volume shaded
in red, as shown in fig. 5.17, is selected. If the x- coordinate varies from −a to
a, the whole steel half ring is addressed. The other coordinates then need to be
expressed as the functions of x variable. From the circular-shape cross section
along the ring, we can write
x2 + y 2 = a2
Or √
y= a2 − x2
for the upper half of the circle. The differential volume can now be expressed as
√
dV = 2 a2 − x2 dx × π (R − x)
Using the result of Ex. 5.1, the centroid of the differential volume is located
at
(R − x) sin (π/2) 2
rc = = (R − x)
π/2 π
The center of mass G of the steel half ring can now be determined by
substituting the above quantities into the moment equation:
R
rV = rc dV
Z a √ Z a
2 √
r 2 2
2π (R − x) a − x dx = (R − x) × 2π (R − x) a2 − x2 dx
−a −a π
a2 + 4R2
r=
2πR
Figure 5.18: The centroid of a complex body may be determined discretely ([1],
pp. 256)
object with arbitrary shape, much like the numerical integration techniques per-
formed by mathematical programs to solve the integration that has no analytical
solution. Error magnitude of the result may be specified and achieved by modi-
fying the shape and size of the elements. Formulae to determine the coordinates
of the center of mass discretely can be expressed as follow:
P
mi xi
X = Pi
m
Pi i
mi y i
Y = Pi (5.7)
i mi
P
mi z i
Z = Pi
i mi
If the density is uniform throughout the object, the mass m in the above
equations may be replaced by the length, l, the area, A, or the volume, V ,
depending on its geometry. The equations will now determine the coordinates of
the centroid discretely.
For the object that has holes or hollows inside its body, they will be accounted
by subtracting their moment contribution from the fictitious body with no such
features. That is the mass or volume will become negative. Also, it is advised
to tabulate the related arguments for the systematic summation purpose. Such
arguments are, for example, mi , rci, and mi rci .
Example 5.5 ([2], Prob. 5/54) Determine the coordinates of the centroid of
the shaded area.
Solution: We divide the area into three parts: the rectangle, the triangle,
and the negative area circle. Each of these can be determined readily their area
and the centroid as shown in Table 5.1 with respect to the given coordinate
frame in fig. 5.20.
Table 5.1: Table of areas, their moments, and centroids of fig. 5.19
After the table is created, the centroid of the shaded area can be determined
directly from eq. (5.7) as
h P i
i Ai x i
X= P X = 0.02+8.4375E-3−2.262E-3
0.1+0.01875−0.0113
= 243.6 mm
i Ai
h P i
Ay
Y = Pi Ai i i Y = 0.0125+1.5625E-3−1.4137E-3
0.1+0.01875−0.0113
= 117.7 mm
i
Example 5.6 ([2], Prob. 5/72) A cylindrical container with an extended rect-
angular back and semicircular ends is all fabricated from the same sheet-metal
stock. Calculate the angle α made by the back with the vertical when the
container rests in an equilibrium position on a horizontal surface.
Solution: Since the container rests in an equilibrium position, its free body
diagram will be drawn and the equilibrium condition will be applied. This will
allow us to solve for the required angle α. The force acting on the body would
be the ground reaction force and the gravitational force. Hence it is necessary
to locate its center of gravity. Since the container is fabricated from the same
material, the center of gravity would be the same as the centroid. In addition,
the metal sheet used has uniformly ‘thin’ thickness. Therefore, the centroid of
the volume would be the same as the centroid of the area of the sheet.
We divide the container into four parts: the rectangular back, two semi-
circular ends, and the half-cylindrical shell. Each of these can be determined
readily their area and the centroid as shown in Table 5.2 with respect to the
given coordinate frame in fig. 5.22.
Table 5.2: Table of areas, their moments, and centroids of fig. 5.21
After the table is created, the centroid can be determined directly from
eq. (5.7) as
h P i
Ax 16E6+12πE6+4.5πE6
X = Pi Ai i i X = 8E4+6πE4+1.125πE4+1.125πE4 = 200 mm
i
h P i
Ay -9πE6−1.6875πE6−1.6875πE6
Y = Pi Ai i i Y = 8E4+6πE4+1.125πE4+1.125πE4 = −114.62 mm
i
h P i
Az 8E6−18E6−2.25E6−2.25E6
Z = Pi Ai i i Z = 8E4+6πE4+1.125πE4+1.125πE4 = −42.75 mm
i
These are the x, y, and z coordinates of the center of gravity as well. Looking
perpendicular to the yz-plane, the free body diagram along with its coordinate
system are displayed in fig. 5.22. There will be the equal and opposite ground
reaction force to the gravity force to make the container be in equilibrium. Both
forces are in vertical direction and pass through the center of gravity. Since a
radius is always perpendicular to the tangent line of the circle, the radius will
be the vertical line, coincident with the direction of forces, perpendicular to the
ground. If α is the angle made by the back with the vertical (see fig. 5.22), the
following geometric relationship can be stated:
−1 150 − |Y |
α = tan = 39.6◦
|Z|
Figure 5.23: Creation of the revolved surface object ([1], pp. 266)
Figure 5.24: Creation of the revolved volume object ([1], pp. 266)
Revolved Surface From fig. 5.23, the infinitesimal area of each strip
is
dA = θydL
using the fact that the surface is constructed by revolving the planar curve L
about the x-axis by the angle θ. The infinitesimal area of the strip then would
be the length, θy, times the width, dL. Since the other strips are constructed by
revolving the segment of the line dL with the same angle θ, the total area of the
revolved surface would be Z
A = θ ydL = θyL (5.8)
Revolved Object From fig. 5.24, the infinitesimal volume of each strip
is
dV = θydA
using the fact that the object is constructed by revolving the planar area A about
the x-axis by the angle θ. The infinitesimal volume of the strip then would be
the extruded length, θy, times the cross sectional area, dA. Since the other strips
are constructed by revolving the differential area dA with the same angle θ, the
total volume of the revolved object would be
Z
V = θ ydA = θyA (5.9)
volume may be viewed as the extruded volume of the cross sectional area A along
the circular path of the centroid (with the length θy) about the revolving axis.
The usage of this theorem may be for directly determining the area or volume
of the revolving object. The other usage comes from viewing the equation in
the opposite way. Had we know the surface area or volume of the object by
some means (possibly from the experiment), we may apply eq. (5.8) or (5.9) to
determine the corresponding centroid of the planar curve or the planar area.
Example 5.7 ([2], Prob. 5/89) A hand-operated control wheel made of
aluminum has the proportions shown in the cross-sectional view. The area of the
total section shown is 15, 200 mm2 , and the wheel has a mass of 10 kg. Calculate
the distance r to the centroid of the half-section. The aluminum has a density
of 2.69 Mg/m3 .
Solution: The mass and density of the control wheel are given. Hence its
volume can be determined as
Since the volume is created from the full-round revolution of half of the sec-
tion depicted in fig. 5.25, the planar area to revolve is 15200
2
×1E-6 = 7600E-6 m2 .
From the Pappus theorem, the distance r to the centroid can be calculated as
Example 5.8 ([1], Prob. 5/93) Calculate the mass m of concrete required to
construct the arched dam shown. Concrete has a density of 2.40 Mg/m3 .
Solution: The concrete mass can be determined from the knowledge of the
density and the volume of the dam. From fig. 5.26, this particular dam is
constructed as if we rotate the cross section A-A about the pole of the 200 m
radius circle for 60◦ . Therefore, we may use Pappus theorem in determining the
volume.
First, to avoid direct integration of the differential volume, the centroidal
coordinate with respect to the axis of revolution is necessary. To avoid direct
integration of the differential area, we view the cross sectional area as obtained
by subtracting the quarter area from the square area. Hence the centroidal
coordinate can be determined from the known centroids of the composite areas
in a discrete manner as
80×80×(40+120)− π ×702 ×( 4×70 +120)
h i
y A
P
y = Pi ciAi i y = 80×80−
4
π
×702
3π
i 4
y = 175.52 m
Now the Pappus theorem may be used to calculate the revolved volume as
π π
× 702
[V = θyA] V = 3
× 175.52 × 80 × 80 − 4
V = 468, 985.2 m3
Consequently, the concrete mass required may be determined readily as
Figure 5.27: Free body diagram of an infinitesimal fluid ([1], pp. 308)
we then can apply it to determine the force acting upon the infinitesimal area dA
on the surface of the volume.
Note that all pressure forces in every direction except the vertical one cancels.
From the free body diagram and the equilibrium condition, the pressure force on
the lower face must balance the summation of the force acting on the upper face
and the fluid’s weight. This can be reduced to the following relation:
dp = ρgdh (5.11)
where ρ is the density of the fluid, g is the gravitational acceleration, and dh is
the infinitesimal height measured along the gravitational acceleration direction.
It is observed that the pressure is the function of the dimension in the vertical
direction solely. Additionally, the pressure increases with the depth. In fact, with
a constant ρ, we may integrate the above relation to determine the pressure as
follow:
p = p0 + ρgh = absolute pressure (5.12)
p0 is the pressure at h = 0. Normally p0 is selected to be the pressure at
the sea water surface level, which is called the atmospheric pressure. Mostly,
however, we offset p0 to be 0 (zero). Consequently the pressure p is no longer be
the absolute (or actual) pressure. Rather, it will be the relative pressure to the
atmospheric pressure. This is the same value as that read off the pressure sensor.
Therefore this pressure p is commonly called the gage pressure.
p = ρgh = gage pressure (5.13)
The unit of pressure is the unit of force over the square of the unit of dimen-
sion. For the SI unit system, it turns out to be N/mm2 . To pay tribute to Pascal,
a French mathematician who introduce the notion of the pressure, 1 N/mm2 is
equal to 1 Pa.
Now we are in the position ready to determine the resultant force due to the
pressure distribution on the surface. To simplify the development, we divide the
surface into two categories: the flat surface and the curved surface. Furthermore,
for the curved surface, we confine the scope to the study of the cylindrical
surface with constant width only. Let’s begin with the flat surface case.
Figure 5.28: A flat surface immersed in the fluid medium with the pertinent
coordinates ([1], pp. 312)
Figure 5.29: Pressure distribution overlaying the flat surface ([1], pp. 312)
where y is the distance measured along the projection of the vertical line onto
the surface. x is the width of a specific strip at the y position. Lastly, h is the
vertical distance of the strip measured from the surface level of the fluid. Note
that the pressure at this strip is ρgh (y). Hence the force acting on the strip can
be viewed as the volume of the extruded strip with the extruded length be equal
to the pressure at that strip.
Since the pressure force acts perpendicular to the surface and the surface is
flat, we can obtain the total force over the whole surface by summing the forces
acting on all strips. In continuous domain, this will become the integrating
operation: Z Z
R = dR = ρgh (y) x (y) dy (5.14)
Figure 5.30: Side view profile of the pressure distribution ([1], pp. 310)
due to the linear variation of pressure along the depth, p = ρgh. Therefore we
can interpret the determination of the magnitude of the resultant force
to be equivalent to calculating the volume of the pressure distribution.
This fact allows one to avoid evaluating the direct integration had the volume
of the pressure distribution been readily recognized. The trapezoidal sectional
shape may assist the volume calculation, and furthermore determination of the
line of action of the resultant force.
Revisiting the above derivation and making use of the centroidal coordinate
definition, we can write
Z
R = ρg hdA = ρghA = pav A (5.15)
Figure 5.31: A cylindrical surface with constant width immersed in the fluid
medium ([1], pp. 311)
R
ypxdy
Y = R (5.16)
R
RY = ydR pxdy
The implication of this equation can be obtained from the principle. Line
of action of the resultant force will pass through the centroid of the
volume of pressure distribution (denoted as point C in fig. 5.29). The above
formula yields the y-coordinate of the centroid. Additionally, the piercing point
of the line of action to the flat plate is called the center of pressure (denoted as
point P in fig. 5.29). This point is the projection of the centroid of the volume
of pressure distribution along the pressured force direction onto the plate. Note
that it is usually not the same point as the centroid of the plate area.
A special case is when the width of the plate is constant, i.e. the rectangular
flat plate. In other words, x is constant. This simplifies the above expression for
determining the line of action to be
ydA′
R
Y = R (5.17)
dA′
where dA′ = pdy is the infinitesimal area of the trapezoid along the y-coordinate.
See fig. 5.30. For this case, line of action of the resultant force will pass
through the centroidal y-coordinate of the area A′ of the side view
profile of the pressure distribution. The location of this trapezoidal centroid
may be evaluated using technique of composite bodies of triangle and rectangle,
as depicted in fig. 5.30.
Next let’s analyze a specific case of the curved surface, the cylindrical surface
with constant width.
Figure 5.32: Side view profile of the pressure distribution on the cylindrical sur-
face ([1], pp. 311)
surface of the constant width b immersed in the fluid medium. We would like to
determine the resultant force R acting on it.
Similar approach as the flat surface case may be applied here where the plate
is first divided along the side view profile into infinitesimal strips (see fig. 5.32).
Then the corresponding force magnitude acting on a strip becomes
dR = pdA = ρgh (L) bdL
where L is the curved distance measured along the side view profile of the cylin-
drical surface. h is the vertical distance of the strip measured from the surface
level of the fluid. Again, ρgh (L) is the pressure at the strip of location L. Hence
the magnitude of force acting on the strip can be viewed as the volume of the
extruded strip with the extruded length be equal to the pressure at that strip.
The pressure forces act perperdicular to the strips. However, they are not
parallel as in the case of flat surface which mean that the magnitude of the
resultant force cannot be obtained simply by summing the force magnitude of
every strip altogether. Instead, we have to project the force acting on each strip
into two common components and then sum the components algebraically. The
resultant force may in turn be determined from summing vectorially these two
components back. If x- and y- components have been chosen as denoted in
fig. 5.32, then the magnitude of the resultant force along each direction may be
written as
Z Z
Rx = b (pdL)x = b pdy (5.18)
Z Z
Ry = b (pdL)y = b pdx (5.19)
Note that we cannot evaluate the magnitude of the resultant force by calcu-
lating the volume of pressure distribution, or by using eq. 5.15. This is because
the pressure forces acting perpendicular to the strips are not parallel to each
other and so cannot be summed algebraically. As a result, besides the magni-
tude, we need to specify its direction and line of action. The direction can be
easily determined, for example, by the angle the vector made with the x-direction;
−1 Ry
θ = tan (5.20)
Rx
The line of action may be determined by applying the definition of centroid
to evaluate the x- and y-centroidal coordinates separately. From fig. 5.32, the
x-centroidal coordinate of Ry is determined from the principal of moment
R R
xdRy xpdx
X= R = R (5.21)
R
Ry X = xdRy dRy pdx
Similarly, the y-centroidal coordinate of Rx may be determined as
R R
ydRx ypdy
Y = R = R (5.22)
R
Rx Y = ydRx dRx pdy
The intersection of the line of action of Rx and Ry will be a point through which
the resultant force passes.
An alternative method of determining the resultant force R on the cylindrical
surface, which is easier for most cases, is to apply the equilibrium condition on
block of liquid. Casting this condition is sensible from the intrinsic assumption
of the fluid static, where the liquid is not flowing and incompressible. One can
then think casually the block of liquid as a rigid body object.
A typical free body diagram of a block of liquid column over the cylindrical
surface viewing from the side is shown in fig. 5.33. In the figure, Px and Py are
the pressure resultant forces acting on the top and side faces of the liquid block.
W is the weight of the liquid block passing through the centroid of the area
ABC (the plate has constant width). Note that all of them can be determined
readily compared to directly evaluating the resultant force. The reaction force
R from the cylindrical surface acting on liquid block may be determined from
the equilibrium conditions. Force equilibrium condition will reveal the mag-
nitude while the moment equilibrium condition will help answer the line of action.
Figure 5.33: Free body diagram of a block of liquid column over the cylindrical
surface ([1], pp. 311)
fluid block surface. By the equilibrium condition, this force would have the same
magnitude as the weight of the fluid block and act vertically upward through the
C.M. of the fluid lump. See fig. 5.34c). Hence, this force is called the buoyancy
force. This force is not to be changed for if the fluid block be replaced
with any immersed object of the same shape. Therefore, regarding to the
buoyancy force, the following can be concluded
1. The buoyancy force is the pressure resultant forces exerting by the sur-
rounding fluid on the surface of the immersed object.
3. Its line of action passes through the center of mass of the displaced fluid.
B = ρf Vf g (5.23)
where ρf and Vf are the density and volume of the displaced fluid, respectively.
Suppose the material of the object has the density ρo , and if the object is fully
immersed, its weight will be equal to W = ρo Vf g. When ρo < ρf , the object’s
weight will be less than the buoyancy force. This makes the object to buoy up to
the surface level and possibly expose to the air, which can be think of as another
fluid medium. It will rise up to the equilibrium position where the summation of
the buoyancy forces be equal to the object’s weight.
Substituting the depth of the centroid of the circular area, h, of which its
value is 3000 m into the above equation and assuming the constant compressive
Example 5.10 ([2], Prob. 5/169) The cross section of a fresh-water tank with
a slanted bottom is shown. A rectangular door 1.6 × 0.8 m in the bottom of the
tank is hinged at A and is open against the pressure of the water by the cable
under a tension P as shown. Calculate P .
Solution: Assume that the door has negligible weight and volume. There-
fore its weight and buoyancy force may be omitted. To determine P , we need
to calculate the pressure force acting on the gate and apply the equilibrium
condition to it. Fig. 5.38 shows the free body diagram of the gate and the
corresponding pressure distribution and forces acting on it. Since the door
surface is flat and has constant width, it is appropriate to calculate the pressure
force magnitude from the volume of pressure distribution. Furthermore, we may
divide the volume into two parts – rectangular and triangular volume. Hence
the indicated force magnitudes are
2
−P × 1.6 cos 30◦ + F1 × 0.8 + F2 × 1.6 ×
P
[ MA = 0] 3
=0
P = 12566 N
Example 5.11 ([1]) A rectangular wood block of density 800 kg/m3 floats
stably in the salted water, of which its surface is topped with the spilled oil.
The salt water and the oil have density of 1030 and 900 kg/m3 , respectively.
Determine the height, h, of the block that is not sunk.
Solution: If all of the pressure forces are referred to the atmospheric pres-
sure forces, the free body diagram of the block may be drawn as shown in
fig. 5.40. Recall that the buoyancy force is equal to the weight of the displaced
fluid, we separate it into two parts; one from the oil pressure and the other from
the salt water pressure. Establishing the equilibrium condition in the vertical
h = 70.4 mm
Solution: For the first part, it is assumed that the buoy is weighted at its
base, and so it remains vertical. After drawing the free body diagram of the
buoy in fig. 5.42a), we may apply the equilibrium force condition along the
vertical direction,
P
[ Fy = 0] B − Mg − T = 0
For the case when h = 0.6 m, the buoyancy force would be determined as
π × 0.32
B = ρV g = 1030 × × 1.8 × 9.81 N
4
Substituting the values into the above equilibrium equation, we can solve for the
cable tension
T = 402.7 N
When the cable goes slack, T = 0 and the buoyancy force must decrease to
counter the buoy’s weight solely. Consequently, the new height h must increase.
The new equilibrium equation would become
P
[ Fy = 0] B − Mg = 0
π × 0.32
1030 × × (2.4 − h) × g − 90 × g = 0
4
h = 1.164 m
In the second part, the buoy’s weight passes through its centroid, showing in
the new free body diagram of fig. 5.42b). Therefore, the buoy might be tilted
instead of just upright only. Let the tilted angle be θ. From the free body
diagram, the weight acts through the buoy’s centroid and the buoyancy force
passes through the centroid of the immersed portion. We then can set up the
equilibrium in moment about the base point of the buoy as follow
P
[ Mo = 0] B × 0.75 tan θ − 90 × g × 1.2 sin θ = 0
The buoyancy force magnitude is determined from the weight of the dis-
placed fluid:
π × 0.32 1.5
B = ρV g = 1030 × × ×g
4 cos θ
Substitute this value into the above equilibrium condition, it requires
Therefore,
θ = 0◦ or ± 29.44◦
Example 5.13 ([1], Prob. 5/186) A fresh water channel 3 m wide is blocked at
its end by a rectangular barrier, shown in section ABD. Supporting struts BC
are spaced every 0.6 m along the 3 m width. Determine the compression C in
each strut. Neglect the weights of the members.
Solution: The free body diagram of the barrier is shown in fig. 5.44 show-
ing the pressure distribution, the corresponding resultant force, the compressive
strut forces and the supporting force at the hinge. The magnitude of the
resultant force may be determined from the volume of the pressure distribution
which is
Substituting this value into the moment equilibrium condition about the hinge
point,
C = 2.354 kN
Example 5.14 ([1], Prob. 5/193) The barge crane of rectangular proportions
has a 6 × 15 m cross section over its entire length of 40 m. If the maximum
permissible submergence and list in sea water are represented by the position
shown, determine the corresponding maximum safe mass mo that the barge can
handle at 10 m extended position of the boom. Also find the total displacement
m in the metric tons of the unloaded barge. The distribution of machinery and
ballast places the CG of the barges, minus the mass mo , at the center of the
hull.
Solution: Free body diagram of the crane reveals the pertinent forces and
make clear the equilibrium posture of it. The buoyancy force B is equal to the
weight of the displaced fluid,
40
B = ρV g = 1030 × 15 × 6 × ×g
2
and acts through its C.M. The submerged portion is of prismatic shape and
hence the easiness of determining the location of its centroid, as indicated by
the intersection of the 2 m horizontal line and the 5 m vertical line shown in
fig. 5.46. The maximum permissible safe mass may be determined from setting
the moment equilibrium condition of the crane critical posture about the C.G.
of the barges.
P
[ MG = 0] B cos θ × 2.5 − B sin θ × 1 − mo g × (10 cos θ + 23 sin θ) = 0
mo = 203 Mg
The unloaded weight of the barge itself may be determined using the force
equilibrium condition along the vertical direction;
P
[ Fy = 0] B − mg − mo g = 0
m = 1651 Mg
Example 5.15 ([1], Prob. 5/197) The fresh water side of a concrete dam has
the shape of a vertical parabola with vertex at A. Determine the position b of
the base point B through which acts the resultant force of the water against the
dam face C.
Solution: This problem falls to the case of cylindrical surface with con-
stant width. Therefore, we will use the alternative method of free body diagram
of the block of fluid since it is more convenient. The free body diagram of the
fluid column on the dam face C is shown in fig. 5.48. First, we will determine
the fluid weight. The shape of the dam surface can be described by a vertical
parabola with the vertex at A. If we set the x-y coordinate frame with the porigin
2
at A, the surface may be written in equation of the form y = ax or x = y/a
in the first quadrant.
From fig. 5.47, the dam surface passes through (27, 36). Substitute this point
into the equation, we can solve for the coefficient
y 36
a= 2
= 2
x 27
To determine the fluid weight, we use the density definition which requires
the knowledge of the fluid volume. The fluid volume may be calculated from the
cros-sectional area times the depth of the fluid block. See fig. 5.48.
R
36
V = 0 xdy h = 32 xy|y=36 h = 648h m2
[V = Ah] y=0
Its line of action passes through the centroid of the volume. This point
projected onto the cross-sectional view will be coincident with the centroid of
the area. Hence the distance X shown may be determined as
R R 36 x
AX = xc dA 648 X = 0 2
xdy
X = 10.125 m
There is the pressure force from the sided fluid acting on the vertical left sur-
face. We may determine its magnitude from the volume of pressure distribution,
1
F = × 36 × (1000 × 9.81 × 36) × h
2
The line of action passes through the centroid of the volume of pressure distribu-
tion which, for this simple geometry, is 12 m above the vertex point A.
We decompose the resultant force from the dam acting on the parabolic
surface passing through B into the horizontal and vertical components. It is
obvious that the magnitude of Ry must be equal to the fluid weight. From the
figure, Rx passes through point D which is located at the distance b to the left
of point B. Moment of all forces about D must be balanced. This condition is
used to determine the distance b.
P
[ MD = 0] −F × 18 − W × 10.125 + W × b = 0
b = 28.125 m
Example 5.16 ([1], Prob. 5/220) A flat plate seals a triangular opening in the
vertical wall of a tank of liquid of density ρ. The plate is hinged about the upper
edge O of the triangle. Determine the force P required to hold the gate in a
closed position against the pressure of the liquid.
Solution: This flat plate has varying width along the depth. Hence the
explicit formula of the volume of pressure distribution is not available. One
might resort to determining the volume from integrating the volume of the
infinitesimal strip over the surface and formulating the moment equilibrium
condition about O. However, an extra step of evaluating the line of action of the
pressure force must be performed.
Rather, since the problem does not ask for the pressure force. We may deter-
mine the holding force P by formulating the moment equilibrium condition about
O and dealing with the infinitesimal distributed force directly. See fig. 5.50. A
particular infinitesimal force dR acting at the depth of y below the hinge edge
may be expressed as
dR = pdA = ρg (h + y) xdy
The width x must be expressed in terms of the changing variable y. From the
geometric relationship of similar triangle,
x b y
= → x=b 1−
a−y a a
This leads to
y
dR = ρg (h + y) × b 1 − dy
a
These infinitesimal forces induce the moments around O, for which their
summation must balance the moment produced by the holding force P . Therefore
R y=a
−P a + ydR = 0
P
[ Mo = 0] y=0
Substitute the expression for dR into the equation and evaluate the inte-
gral,
ρgba a
P = h+
6 2
Chapter 6
Introduction to Dynamics
• Time is the measure of the succession of events. Often, we are more interested
in the change of physical quantities with respect to time, e.g. v = dr dt
,
instead of time variable itself.
• Mass is the measure of the inertia of a body. The inertia indicates the resis-
tance to a change in motion.
• Rigid body is a body whose relative movement between its parts are negligi-
ble relative to the gross motion of the body. For example the motion of an
ingot can be analyzed by assuming the object being rigid.
• Nonrigid body is a body whose relative movement between its parts are sig-
nificant relative to the gross motion of the body. Knowledge of the mechan-
ics of the deformable material must be used along with Dynamics in order
to determine the absolute motion of the nonrigid bodies.
Let us consider some examples to see the difference of each term. If we have
an object and consider the very small substance of the body. For differential
element analysis of the body, the small substance can be treated as a particle.
However, the substance must be handled as connecting objects had the molecular
effects in the body are of concern. Or think of an airplane. Even of its huge size,
the whole airplane may be modeled as a point in flight speed analysis along the
route. But if the rotational motion, such as yawing or pitching, of the airplane
body is important, its size does matter.
The next two examples are to show whether an object is considered rigid or
nonrigid depends on how much detailed of the problem we would like to analyze.
Truss can just be looked as a rigid body for the preliminary design of truss
structure. Buf we must think of the truss elasticity if we were to choose the
material for that truss. A stiff linkage of the robot may be considered a rigid
body. However, the n-connecting linkages, treated as a whole, to form the robot
arm is an example of nonrigid body. Note the body-fixed inertia of the nonrigid
body is not constant.
ΣF = 0 ⇔ a = 0
ΣF = ma (6.1)
This fact is used very often in drawing the free body diagram (FBD).
F = attraction force
gravitational constant value = 6.673 × 10−11 m3 / kg · s2
G =
m = mass of the involving bodies
r = distance between the bodies
Hence there is always the attraction force between the earth and the object.
This gravitational force is called the weight of the body.
Gme
W =m = mg (6.3)
r2
where
Chapter 7
Kinematics of Particles
7.1 Introduction
Kinematics is the study of the motion of bodies with no consideration to the forces
that accompany the motion. Kinematics analysis is a prerequisite to kinetics,
which is the study of the relationship between the motion and the corresponding
forces that causes the motion or are generated as a result of the motion.
In this chapter, we are interested in kinematics of particles. A particle is
a body whose physical dimensions are so small compared with the radius of
curvature of its path. Therefore an airplane may be considered as a particle if
the problem of interested is related to the flying route between two cities, even
the size of the airplane is huge compared to a human. In other words, the body
rotation effect is insignificant compared to the translation.
Figure 7.1: Several descriptions of the particle position ([3], pp. 22)
to the projected point is denoted as r. The angle measured from the reference
to the radial direction is called θ. And the distance measured along the line
perpendicular to the plane from the projected point to the particle is called z.
{r, θ, z} forms the cylindrical coordinate system that might be used to describe
the position.
Yet another way to specify where the particle be is to employ the spherical
coordinate system. This is possible by first determining the origin, the reference
plane, and the reference direction. The distance measured from the origin to the
point P is called R. Next we project this distance vector onto the reference plane.
The angle measured from the reference line to the projected radial vector is noted
as θ, or the azimuth angle. Finally, we measure another angle from the radial
vector to the distance vector. This angle is called the elevation angle, denoted as
φ.
Each description has its own advantage and hence is suitable to some par-
ticular tasks. For example, the spherical coordinate system is largely employed
in navigation system, such as rocket, radar, or satellite system. The cylindrical
coordinate system might be appropriate for some robots that have this structure
for pick and place operation. The rectangular coordinate system is suitable for
general tasks since it is easiest to understand the position of an object in three
dimensions by telling the width, height, and depth.
Motion of P We are not satisfied only just being able to describe where
the particle is because it is moving. We must be able to describe the motion
of the particle. In other words, we will also be interested in the change of its
position. This is achieved in general by differentiating the position with respect
to time. Before the analysis, first we must decide which kind of the coordinate
system we will be using. Other than intrinsically different types of coordinate
systems described above, they may also be categorized according to whether the
frames used are fixed (not moving) or moving.
Fixed reference coordinate frame In this case, the coordinates are measured
with respect to the fixed reference frame. For example, we are standing on the
ground and seeing the airplane flying. That is we are observing the motion of
the airplane from the fixed reference frame. What we are measuring is really the
absolute motion of the particle. Hence this is called absolute motion analysis.
Unfortunately, often the motion is convolved. It is quite inconvenient to
describe the position of the particle using the fixed reference frame solely. In
such situation, it is more advantage to employ the moving reference coordinate
frame.
Figure 7.2: Rectilinear motion showing the distance s and its change ∆s measured
along the straight line path ([3], pp. 22)
fly on another plane following the target plane in a manner that the relative
location between them is fixed. Consequently, we will see the plane not moving
in this moving reference frame! What we are measuring is just the motion of the
particle relative to the moving frame. When we add this relative motion to the
motion of the reference frame, we will recover the absolute motion.
The velocity or, to be precise, the instantaneous velocity is the average ve-
locity as the time interval ∆t approaches zero. Recalling the definition of the
differentiation, this may be written as
ds
v = lim ∆s/∆t = = ṡ (7.2)
∆t→0 dt
In words, the velocity is the time rate of change of the position coordinate.
We can do what we just did with the position coordinate to the velocity to
obtain the acceleration. Let the change in the velocity occurred during the time
dv d2 s
a = lim ∆v/∆t = = v̇ = 2 = s̈ (7.4)
∆t→0 dt dt
In words, the acceleration is the time rate of change of the velocity. Multiplying
eq. 7.2 and 7.4 together, and integrating with respect to time, we obtain the
relationship:
vdv = ads (7.5)
which is sometimes useful as we shall see.
3. Both velocity and acceleration are free vectors. Therefore their changes
happen from two sources; change in magnitude and change in direction.
In solving the problems, one may need to perform the integration to the
definitions and the relationship of kinematics parameters in eq. 7.2, 7.4, and 7.5.
Some of their physical meanings are, for example,
Z s2 Z t2
ds = vdt → s2 − s1 = area under v−t curve
s1 t1
Z v2 Z t2
dv = adt → v2 − v1 = area under a−t curve
v1 t1
v2 s2
v22 − v12
Z Z
vdv = ads → = area under a−s curve
v1 s1 2
ds
v= → v = slope of s−t curve
dt
dv
a= → a = slope of v−t curve
dt
Figure 7.3: Plots of the displacement, velocity, and acceleration with respect to
time and their related quantities ([3], pp. 24)
These relationships are depicted in fig. 7.3. The interpretation suggests an alter-
native way of numerical methods in evaluating the kinematic parameters.
Typically the acceleration (in general could be a function of displacement,
velocity, or time explicitly) of the object will be given. Then the problem will
ask for other kinematical parameters. Conceptually, these can be determined
by integrating the acceleration and solve for the quantities, such as velocity,
position, or time. To make the reader become acquainted with such problems, in
the following, we will consider several case studies where the given acceleration
be functions of other kinematic quantities.
b) a = f (t) In this case, the acceleration would be the result of the synthetic
forces applied for specific purposes. For example, the explosion force in the piston
cylinder from appropriate control of the valve timing and fuel amount causes the
automobile to run at a desired motion. Compared to the constant acceleration,
we cannot pull out the acceleration of the integral sign. The pertinent kinematic
quantities may then be determined as follow:
Z v Z t Z t
dv = f (t) dt → v = vo + f (t) dt
vo 0 0
Z s Z t Z t Z tZ θ
ds = vdt → s = so + vdt = so + vo t + f (ζ) dζdθ
so 0 0 0 0
s = h (t)
On the other hand, if the information on time is not provided, we may derive
the displacement relationship from
Rv v
Rs
[vdv = ads] vo f (v)
dv = so
ds
for which we may solve for the displacement as a function of the velocity
Z v
v
s = so + dv = g (v)
vo f (v)
d) a = f (s) Many devices possess the spring behavior to some degrees. That is,
they provide the resistive forces that are function of the deformation. Another
type of forces in this category is the field force, such as the attraction force. From
the gravitational law, the attraction force will be inversely proportional to the
square of the separated distance. These forces induce the acceleration that is a
function of the displacement. Hence we may solve for the velocity
Z v Z s Z s
2 2
vdv = f (s) ds → v = vo + 2 f (s) ds
vo so so
as a function of displacement
v = g (s)
The displacement may be determined explicitly by
Rs 1
ds
v= dt
t= so g(s)
ds
s = h (t)
Example 7.1 ([3], Prob. 2/23) Small steel balls fall from rest through the open-
ing at A at the steady rate of 2 per second. Find the vertical separation h of two
consecutive balls when the lower one has dropped 3 meters. Neglect air resistance.
Solution: The accceleration for the free-falling object must be equal to the
constant of gravity. Also, we can integrate for the velocity and the displacement
explicitly. In other words,
a=g
v = vo + gt
s = so + vo t + gt2 /2
Since the ball is dropped from rest at the reference level of the opening, we have
the initial conditions
vo = 0, so = 0
Hence
s = gt2 /2
Apply this relation to the ball already dropped by 3 m, the time spent would
be
3 = gt2l /2, tl = 0.782 s
Since the balls are released with the rate of 2 per second, the time the consecutive
ball spent so far would then be
tu = tl − 0.5 = 0.282 s
Use the above relation again to solve for the distance travelled, which is related
to the vertical separation h:
su = 3 − h = gt2u /2 → h = 2.61 m
and
R sC R t/3600 R t/3600
sD
ds = 0
vdt = 0
60dt, 3 − sC = t/60
The distance traveled during the braking may be calculated indirectly from
the area under v-t curve as (fig. 7.7)
1 4
sC − sB = × × (100 + 60) = 4/45
2 3600
These three equations may be used to solve for the time spent
t = 65.5 sec
s = sB = 1.819 km
Example 7.3 ([3], Prob. 2/27) The 350-mm spring is compressed to a 200-mm
length, where it is released from rest and accelerates the sliding block A. The
acceleration has an initial value of 130 m/s2 and then decreases linearly with
the x-movement of the block, reaching zero when the spring regains its original
350-mm length. Calculate the time t for the block to go a) 75 mm and b) 150 mm.
Solution: We may draw the change in the acceleration with respect to the
block motion from the system description as shown in fig. 7.9. Accordingly, the
acceleration-displacement relationship may be written as
130
a=− x = −866.7x
0.15
This graph inspires us to apply the velocity-displacement relationship
through the acceleration:
Rv Rx Rx
[vdv = ads] 0
vdv = −0.15
ads = −0.15
−866.7sds
√
v 2 = −866.7x2 + 19.5, v = 29.44 0.0225 − x2
assuming that the block move positively to the right. The time spent for the
block to travel to the distance x would then be
Rt Rx 1
Rx
ds √ 1
v= dt 0
dt = −0.15 v
ds = −0.15 29.44 0.0225−s2
ds
h x πi
t = 0.034 sin−1 +
0.15 2
By substituting the position of the box x = 75 mm, the time spent would be
71.14 ms. For the position at x = 150 mm, the time spent is 106.7 ms.
Alternatively, one may use the definition of the acceleration
[a = s̈] ẍ + 866.7x = 0
We may deduce the initial conditions from the problem statements as follow.
Substituting these conditions into the above solution form, we may be able to
solve for the coefficients:
A = 0, B = −0.15
Therefore, the motion of this block behaves according to
for which we may inversely solve for the time spent. The answers agree with the
other approach.
Example 7.4 ([4], Prob. 2/31) A train that is traveling at 130 km/h applies
its brakes as it reaches point A and slows down with a constant deceleration. Its
decreased velocity is observed to be 96 km/h as it passes a point 0.8 km beyond
A. A car moving at 80 km/h passes point B at the same instant that the train
reaches point A. In an unwise effort to beat the train to the crossing, the driver
’steps on the gas’. Calculate the constant acceleration a that the car must have
in order to beat the train to the crossing by 4 s and find the velocity v of the car
as it reaches the crossing.
Solution: Since all objects in the problem move with constant accelera-
tion, we may refer to the instant formula applied to the constant acceleration
case, namely,
v = vo + at
v 2 = vo2 + 2a (s − so )
s = so + vo t + at2 /2
For the train, the driver applies the brake causing the constant deceleration.
This makes the velocity being reduced from 130 to 96 km/h for the distance
traveled 0.8 km. Therefore the deceleration may be determined:
By this deceleration, the train would reach the intersection at the time counted
from the time it reaches point A
For the car to beat the train to the crossing by 4 s, the time it must spent
counting from point B must be
Example 7.5 ([3], Prob. 2/44) The horizontal motion of the plunger and shaft
is arrested by the resistance of the attached disk that moves through the oil
bath. If the velocity of the plunger is vo in the position A where x = 0 and t = 0,
dv
Rv 1
Rt
dt, v = vo e−kt
a= dt vo −kv
dv = 0
ds
Rx Rt vo
v e−kt dt, x = 1 − e−kt
v= dt 0
ds = 0 o k
To express the velocity in terms of the displacement, we may use the fol-
lowing differential relationship:
Rv 1
Rx
[vdv = ads] vo −kv
vdv = 0
ds, v = vo − kx
Example 7.6 ([4], Prob. 2/44) The electronic throttle control of a model train
is programmed so that the train speed varies with position as shown in the plot.
Determine the time t required for the train to complete one lap.
Solution: From the velocity-displacement graph, the train travels with the
constant velocity of 0.25 m/s in the straight line segments. Hence, the time used
for each segment may be simply calculated by
Next, for a quarter of the circle, i.e. from 2 to 2 + π2 m, the velocity decreases
dv
v dv dv
vdv = ads, a = dt ds
= dt
Example 7.7 ([4], Prob. 2/52) A bumper, consisting of a nest of three springs,
is used to arrest the horizontal motion of a large mass that is traveling at 40 m/s
as it contacts the bumper. The two outer springs cause a deceleration propor-
tional to the spring deformation. The center spring increases the deceleration
rate when the compression exceeds 0.5 m as shown on the graph. Determine the
maximum compression x of the outer spring.
Solution: The associated deceleration of the mass with respect to the dis-
placement is given. This leads us to apply the following relationship
R0
[vdv = ads] 40
vdv = area under a−s curve
In particular,
0 − 402 1 1
= − × 0.5 × 1000 − × (x − 0.5) × (1000 + 1000 + 4000 (x − 0.5))
2 2 2
for which we recognize the negative area due to the deceleration. Also, we employ
the equation of the second line segment in determining the height of the end
point. Solving the equation for the maximum compression x corresponding to
null motion, we have
x = 0.831 m
Figure 7.14: Position vector describing the location of the particle ([3], pp. 40)
the curved path. Hence it is always greater than or equal to the magnitude of
the displacement vector.
Figure 7.15: Velocity vector always tangent to the path of the particle ([3], pp.
40)
vector that has the direction tangent to the path. Had we know the trajectory of
the particle, we will know the direction of the velocity at every instant as drawn
in fig. 7.15.
If we consider the magnitude of the displacement alone and treating the time
interval ∆t → 0, the magnitude value of the average velocity would approach
|∆r| ∆s
|v| = lim = lim
∆t→0 ∆t ∆t→0 ∆t
using the fact that |∆r| will approach ∆s as ∆t → 0. This leads us to the
definition of the speed,
ds
v= = ṡ = |v| (7.10)
dt
Consequently, the speed may be determined from the magnitude of the velocity.
Figure 7.16: Acceleration vector always pointing toward the region containing
the center of curvature ([3], pp. 40)
along the wide and long directions. For this coordinate system, we have to set up
the coordinate frame which consists of the origin O, and two mutually perpendic-
ular directions denoted by x and y symbols for x- and y-directions. Associated
with these directions are the unit vectors i and j, respectively. From the elemen-
tary properties of the vectors, we may write the position vector describing the
location of the particle specifically as
That is the position vector is decomposed into two subvectors along the x- and
y-directions. If these directions are unchanged, we may determine the velocity
from its definition:
v = ṙ = ẋi + ẏj = vx i + vy j (7.14)
di dj
with dt = 0 and dt = 0. The acceleration may be determined by differentiating
the velocity, of which its straightforward relationship is
Figure 7.17: Position, velocity, and acceleration described in x-y coordinate sys-
tem ([3], pp. 43)
From the above kinematic equations, if the x- and y-components of the acceler-
ation are independently generated or determined, i.e. ax and ay at a specific point
are related by the same instant of time only, we may integrate each individual
expression with respect to time to obtain the velocity and position components.
One may imagine the curvilinear motion be generated from the superposition of
two perpendicular rectilinear motions simultaneously. In fact this kind of motion
is found in machines like x-y plotter or the gantry crane in the factory. Therefore
the rectangular coordinate system is suitable for the complex motion which is na-
tively generated from the mutually orthogonal rectilinear independent motions
simultaneously.
If one has solved for the parametric equations in time x = f1 (t) and y = f2 (t),
the equation of the curved path, written in general form as y = f (x), may
be determined by algebraically eliminating the time parameter t in those two
parametric equations. Consulting the diagram shown in fig. 7.17, some of the
common quantities may be calculated from the kinematical vector components
as shown below.
vy dy
direction of the velocity : tan θ = = (7.16)
vx dx
q
speed : v = vx2 + vy2 (7.17)
q
magnitude of the acceleration : a = a2x + a2y (7.18)
Projectile Motion A typical motion that is suitably represented using the rect-
angular coordinate system is the projectile motion. It is the motion of the thrown
object. For the first-run analysis, we shall neglect the aerodynamic drag force.
Also it is assumed that the effects of the curvature of the rounded earth and its
rotation are small enough. The range of the altitude of the whole motion is in
the order that the gravitational force variation may be omitted.
Figure 7.18: Motion profile of the projectile motion showing variations in the
velocity ([3], pp. 44)
With these simplification, the resulting trajectory affixed with the velocity at
certain points are shown in fig. 7.18. We set up the coordinate system such that
its origin coincides with the initial point. The axes directions are chosen to point
in the horizontal and vertical directions, lining up the y-direction opposite to the
gravitational acceleration direction. Therefore, we may write
ax = 0, ay = −g
vy2 = vy o 2 − 2g (y − yo )
Example 7.8 ([4], Prob. 2/81) A particle is ejected from the tube at A with a
velocity v at an angle θ with the vertical y-axis. A strong horizontal wind gives
the particle a constant horizontal acceleration a in the x-direction. If the particle
strikes the ground at a point directly under its released position, determine the
height h of point A. The downward y-acceleration may be taken as the constant g.
vx o = −v sin θ, xo = 0
vy o = v cos θ, yo = 0
We then integrate the acceleration relationship to obtain the velocity and position
as a function of time:
1
ax = a, vx = −v sin θ + at, x = −vt sin θ + at2
2
1
ay = g, vy = v cos θ + gt, y = vt cos θ + gt2
2
At the point where the particle hits the ground, x = 0 and y = h. Using the
x (t) relationship, we may solve for the hitting time
at 2v sin θ
0=t − v sin θ , t =
2 a
Substituting the time t into y (t) relationship, the height h of point A may then
be determined:
2v 2 g
h= sin θ cos θ + sin θ
a a
Example 7.9 ([3], Prob. 2/75) Electrons are emitted at A with a velocity u at
the angle θ into the space between two charged plates. The electric field between
the plates is in the direction E and repels the electrons approaching the upper
plate. The field produces an acceleration of the electrons in the E-direction of
eE/m, where e is the electron charge and m is its mass. Determine the field
strength E that will permit the electrons to cross one-half of the gap between
the plates. Also find the distance s.
Solution: The electron traveling between two plates is subject to the elec-
tric field force of which its acceleration along the horizontal and vertical
components may be written as
eE
ax = 0, ay = −
m
See fig. 7.22. The discharged electron will have the zero velocity component in
y-direction when it reaches one-half of the gap. Therefore,
0 = (u sin θ)2 − 2 eE
m
× b
2
The required field strength that permits the electron to cross one-half of
the gap would follow the relationship that
mu2
E= eb
sin2 θ
When the electron falls back to the emitting plate, its coordinate values would
be (s, 0). Substitute these values into the x-coordinate relationship, we may solve
for the time
s
x = s = ut cos θ, t =
u cos θ
Use this expression to eliminate the time in the y-coordinate relationship, the
horizontal distance s may be determined.
s 2b
0 = s tan θ 1 − tan θ , s = 0 or
2b tan θ
Example 7.10 ([4], Prob. 2/87) Water is ejected from the nozzle with a speed
vo = 14 m/s. For what value of the angle θ will the water land closest to the wall
after clearing the top? Neglect the effects of wall thickness and air resistance.
Where does the water land?
The corresponding nozzle angle may be solved for after substituting the time
back into either one of the positional relationship.
θ = 50.64◦ , 41.43◦
Trajectory sketch of these two solutions is drawn in fig. 7.23. It is seen that
path (1) with the nozzle angle of 50.64◦ will make the water jet land closest to the
wall. For the coordinate frame employed, coordinate values of the place where
the water lands are (x, 0) where x is the horizontal distance measured from the
nozzle A. The y-coordinate relation will allow us to solve for the spent time
1
0 = 0.3 + 14t sin 50.64 − gt2 , t = 2.234 sec
2
Then, the associated x-coordinate value would be
Hence the water lands at 19.835 − 19 = 0.835 m to the right of the wall.
Solution: Two forces act on the bullet, namely the fluid drag force and
the gravitational force. The drag force produces the deceleration that is
proportional to the velocity which, for the given coordinate frame, we may write
g g
e−kt −
vy = vo sin θ + k k
dy
Rt g g
Ry
e−kt −
vy = dt 0
vo sin θ + k k
dt = 0
dy
g
1
1 − e−kt − kg t
y= k
vo sin θ + k
Example 7.12 ([3], Prob. 2/100) A projectile is launched with speed vo from
point A. Determine the launch angle θ that results in the maximum range R up
the incline of angle α (where 0 ≤ α ≤ 90◦ ). Evaluate your results for α = 0, 30,
and 45◦ .
Solution: For the projectile motion under the gravity field, ax = 0 and
ay = −g. Integrating them with respect to time and applying the initial velocity
conditions, the velocity components would be
vx = vo cos θ, vy = vo sin θ − gt
Then, we integrate the velocity expressions to determine the displacement along
x- and y-directions as
1
x = vo t cos θ, y = vo t sin θ − gt2
2
Consider the particle meeting the incline at B, we may use the x-coordinate
expression to solve for the corresponding time as
R cos α = vo t cos θ
R cos α
t=
vo cos θ
Substitute the time into the y-coordinate expression to construct the relation
between R and θ;
2
R cos α g R cos α
R sin α = vo sin θ −
vo cos θ 2 vo cos θ
which, after some manipulation, may be reduced to
Since the change in the projectile launch angle θ causes the change in the
inclined distance R, the maximum distance happens when dR dθ
= 0. Therefore we
may take the total differentiation on the previous relation,
−2vo2 tan α (2 cos θ sin θdθ) = vo2 (2 cos 2θdθ) − g cos αdR
Example 7.13 ([4], Prob. 2/95) Determine the equation for the envelope a
of the parabolic trajectories of a projectile fired at any angle but with a fixed
muzzle velocity u. (Hint: Substitute m = tan θ, where θ is the firing angle with
the horizontal, into the equation of the trajectory. The two roots m1 and m2 of
the equation written as a quadratic in m give the two firing angles for the two
trajectories shown such that the shells pass through the same point A. Point A
will approach the envelope a as the two roots approach equality.) Neglect air
resistance and assuming g is constant.
ax = 0, vx = u cos θ, x = ut cos θ
1
ay = −g, vy = u sin θ − gt, y = ut sin θ − gt2
2
for a particular firing angle θ and velocity u. Next, we should determine the
projectile trajectory by eliminating the common time parameter t. The trajectory
equation is then
x2
g
y = x tan θ −
2 u2 cos2 θ
which confirm the parabolic shape of the path.
Suppose we would like the trajectory to pass through a point A that has
the coordinate (x, y) by varying the firing angle solely. The required angle may
be determined from the above trajectory equation. The angle θ appears as the
variable of two transcendental functions, namely tan θ and cos θ. Both are related
by the following trigonometric identity:
1
1 + tan2 θ = sec2 θ =
cos2 θ
Let m = tan θ. Hence cos12 θ = 1 + m2 . Make use of these new expression with
the trajectory equation, we may rearrange and rewrite it as
a quadratic function of the variable m. Hence there are two roots according
to the fact that point A may be reached from two distinct paths as seen in
fig. 7.26. Interestingly, as those two paths approach each other, the point (x, y)
will become a point on the envelope of the family of trajectories. Consequently,
the two distinct roots must become two repeated roots. This calls for the zero
discriminant. Specifically,
2
2xu2 − 4gx2 2yu2 + gx2 = 0
In other words, the following relation must be hold for any point (x, y) on the
envelope
u2 gx2
y= −
2g 2u2
That is, it is the envelope equation.
Figure 7.27: Successive n-t coordinate frames showing the instantaneous center
of curvature ([3], pp. 55)
0en + 0et . For the velocity, one may be tempted to take the derivative of the
position vector with respect to time. This would give us the null velocity, which
is quite surprising since no matter how the particle actually moves, the obtained
velocity would always be zero. What went wrong is that the position is the relative
position. In other words, we measure the position of the particle with respect to
the moving frame. The resulting velocity would then be the relative velocity. An
observer moving with the n-t frame would see the particle not moving.
Nevertheless, we will resort to considering fig. 7.27 to determine the absolute
′
velocity. At a moment, the particle is at A. After time elapsed, it moves to A
′
along the curved path. As the elapsed time converges to zero, point A and A
becomes coincident. The displacement becomes smaller and approach ds = ρdβ,
where ρ is the radius of curvature of the path and dβ is the sweeping angle the
radius of curvature traveled during that small interval. Taking the limit of the
displacement as time interval converges to zero is equivalent to differentiating
the displacement with respect to time, which is, according to the definition, the
speed. Therefore,
ds ρdβ
v= = = ρβ̇
dt dt
Since the n-t coordinate frame has the n-axis pointing toward the center of
curvature and perpendicular to the trajectory, it is implied that the t-axis, which
is tangent to the path, will be in the direction of the velocity vector. Therefore,
the velocity vector may be expressed as
That is, the velocity points along the t-axis. As the particle travels along +t-
direction, β̇ > 0.
Differentiating the (absolute) velocity to determine the acceleration, we have
dv
a= = v ėt + v̇et
dt
Chulalongkorn University Phongsaen PITAKWATCHARA
which can be further simplified into fundamental terms. Consider the change in
et , which is a free unit vector. Therefore, the change will get involved with the
change in its direction only. Figure 7.28 illustrates et being rotated by the small
′
angle dβ to a new vector et . Difference between these two vectors is the vector
det :
det = |et | dβen = dβen
Taking the limit as time duration approaches zero, we have
v2
an = v β̇ = ρβ̇ 2 =
ρ
at = v̇ = s̈ = ρβ̈ + ρ̇β̇
From the last relationship, one may observe that the tangential acceleration be-
haves as if it is the rectilinear acceleration. That is it is the time derivative of
the velocity, or the second time derivative of the displacement. We can extend
our perspective on rectilinear motion to the curvilinear motion as long as the re-
lated components are all in tangential direction. The normal direction quantities
are there to make the path deviate from the straight line path. Therefore the
following relation holds as well:
vdv = at ds
Consider fig. 7.29 which shows the change in velocity and its components.
As time duration approaches zero, an indicates the change in the direction of v
while at indicates the change in the magnitude. Note that the component an will
always direct toward the center of curvature. The component at , however, will
be in +t-direction if the speed v is increasing. It may be in −t-direction had the
speed is decreasing. Consequently though the acceleration direction does not lie
in either one pricipal direction solely as for the velocity case (in the tangential
direction), it will tend to the side (divided by the curvature tangent) where the
center of curvature lives.
A final comment is that the kinematics quantities are intrinsic. Even they
may be represented in some specific coordinate frame, the inherent property such
Figure 7.29: Infinitesimal change in the velocity and the relationship with the
acceleration represented in n-t coordinate system ([3], pp. 55)
as the magnitude and the direction must not be changed. For the problem using
the x-y and the n-t coordinate frames, the following may be stated:
v = vx i + vy j = vet
a = ax i + ay j = an en + at et
Example 7.14 ([3], Prob. 2/122) The camshaft drive system of a four-cylinder
automobile engine is shown. As the engine is revved up, the belt speed v changes
uniformly from 3 m/s to 6 m/s over a 2 second interval. Calculate the mag-
nitudes of the accelerations of point P1 and P2 half way through this time interval.
Solution: Since the timing belt velocity is increased uniformly, the (tangential)
acceleration is constant. Therefore
dv ∆v 6−3
at = = = = 1.5 m/s2
dt ∆t 2
This constant acceleration also implies that, at the time half of the interval, the
belt speed is
3+6
v= = 4.5 m/s
2
Point P 1 is moving in circular path around the sprocket. Therefore it will
have both the tangential and normal components of the acceleration. The
For point P 2, it is moving along the tangent between the camshaft and the
intermediate sprocket. Hence the path is straight line and so the acceleration is
just the tangential acceleration of the belt.
aP 2 = at = 1.5 m/s2
Example 7.15 ([3], Prob. 2/126) A baseball player releases a ball with the
initial conditions shown in the figure. Determine the radius of curvature of the
trajectory a) just after release and b) at the apex. For each case, compute the
time rate of change of the speed.
Solution: We must use the n-t coordinate system since the problem asks
for the radius of curvature ρ. Refer to the velocity and acceleration sketches
shown in fig. 7.32. At just after release, the normal acceleration component is
an = g cos 30
because the throwing angle is 30◦ up with respect to the horizontal line. From
this, we may determine the radius of curvature
302
[an = v 2 /ρ] g cos 30 = ρ
, ρ = 105.9 m
The time rate of change of the speed is the tangential acceleration, which
is obvious from the sketch as
Same quantities for the baseball at the apex may be determined in the same
manner. At the apex,
Example 7.16 ([4], Prob. 2/112) Pin P in the crank P O engages the horizontal
slot in the guide C and controls its motion on the fixed vertical rod. Determine
the velocity ẏ and the acceleration ÿ of guide C for a given value of the angle θ
if a) θ̇ = ω and θ̈ = 0 and b) if θ̇ = 0 and θ̈ = α.
[v = rω] v = rω
[an = v 2 /r] an = rω 2
[at = rα] at = 0
v=0
an = 0
at = rα
ẏ = rω sin θ, ÿ = rω 2 cos θ
Example 7.17 ([4], Prob. 2/129) The pin P is constrained to move in the slot-
ted guides that move at right angles to one another. At the instant represented,
A has a velocity to the right of 0.2 m/s which is decreasing at the rate of 0.75
m/s each second. At the same time, B is moving down with a velocity of 0.15
m/s which is decreasing at the rate of 0.5 m/s each second. For this instant,
determine the radius of curvature ρ of the path followed by P . Is it possible to
determine also the time rate of change of ρ.
v = 0.2i − 0.15j
a = −0.75i + 0.5j
Recalling the basic fact that velocity vector lies along the t-axis, we are able to
determine its direction as shown in fig. 7.36. Then, we will be able to determine
the direction of the n-axis from the direction of the t-axis and the acceleration.
That is, it will be making right angle to the t-axis, where the positive direction
will be pointing to the side where the direction of the acceleration resides. As
a result, we will be able to completely indicate the n-t coordinate frame. See
fig. 7.36.
Knowing the orientation of the frame n-t with respect to frame x-y, the unit
vector et may be expressed as
et = 0.8i − 0.6j
Hence the acceleration components in the normal and tangential directions may
be determined as follow:
at = (a · et ) et = −0.72i + 0.54j
an = a − at = −0.03i − 0.04j
Now we may determine the radius of curvature by
0.252
[an = v 2 /ρ] ρ= 0.05
= 1.25 m
For the question of determining ρ̇, consider the following kinematic relation
which contains this term;
v
at = v̇ = ρβ̈ + ρ̇β̇ = ρβ̈ + ρ̇
ρ
Currently, this equation has two unknowns, namely ρ̇ (the one we need) and β̈.
Therefore, ρ̇ cannot be determined until β̈ is known.
Example 7.18 ([3], Prob. 2/134) As a handling test, a car is driven through
the slalom course shown. It is assumed that the car path is sinusoidal and that
the maximum lateral acceleration is 0.7g. If the testers wish to design a slalom
through which the maximum speed is 80 km/h, what cone spacing L should be
used?
Solution: From the geometry of the path in fig. 7.37, we may write down
its equation as π
y = 3 sin ωx = 3 sin x
L
It is specified that the maximum lateral acceleration during the course is
an = 0.7g. Since an = v 2 /ρ, the value will be at peak if the velocity is maximum
and the radius of curvature becomes minimum. From the sketched sinusiodal
trajectory, the peak point of the waveform will have minimum radius of curvature.
Associated with the maximum speed of 80 km/h, the minimum value would be
2
vmax (80 × 10/36)2
ρmin = = = 71.9 m
(an )max (0.7g)
Now we need to think of some relations which will link the radius of curvature
ρ to the unknown cone spacing L. Since L is contained in the trajectory equation,
we may found the following relation from calculus to be helpful.
3/2
1 + y′2
ρ=
y ′′
y ′ = 0, y ′′ = ∓3ω 2
Substitute these derivatives into the x-y-frame radius of curvature formula, and
solve for the desired L:
1
∓71.9 =
∓3ω 2
π
ω = 0.0681 = → L = 46.14 m
L
Example 7.19 A particle starts from rest at the origin and moves along the
positive branch of the curve y = 2x3/2 so that the distance s measured from the
origin along the curve varies with the time t according to s = 2t3 , where x, y,
and s are in millimeters and t is in seconds. Find the magnitude of the total
acceleration of the particle when t = 1 s.
Solution: Recall that the curvilinear motion is the extension of the straight
line motion with the normal acceleration added. Therefore along the tangential
direction, we have
y = 2x3/2
we have
√ 3
y ′ = 3 x & y ′′ = √
2 x
Therefore, the coordinate x, i.e. the location where we would like to evaluate the
radius of curvature, must be known first.
However the elapsed time is given instead. Consequently, the relationship
between the traveling distance, coordinates, and time must be established. Rela-
tion between distance traveled and time are specified to be s = 2t3 . When t = 1
second, the traveled distance is 2 mm.
The link between the distance traveled and its coordinates may be determined
from the differential geometry as
Z Zq Zq
2 2
ds = (dx) + (dy) = 1 + (y ′)2 dx
Substitute y ′ (x) for this particular path into the above relation and solve for x,
we then have Z x
√
2= 1 + 9xdx
0
x = 0.913 mm, y = 1.746 mm
Now we are ready to determine the total acceleration at t = 1 second. From
the above equations, at that time, v = 6 mm/s and at = 12 mm/s2 . From the
path integration, x = 0.913 mm, y ′ = 2.8665, y ′′ = 1.57, and so ρ = 17.8 mm.
The normal acceleration now can be determined as
62
an = = 2.02 mm/s2
17.8
Consequently the total acceleration becomes
q √
a = a2n + a2t = 2.022 + 122 = 12.17 mm/s2
r = ret (7.23)
Figure 7.38: r-θ coordinate frame attached to the particle and the x-y reference
coordinate frame ([3], pp. 68)
Since this is the absolute position, the velocity may be determined by differ-
entiating the above expression to obtain
v = ṙ = ṙer + r ėr
which can be further resolved into fundamental terms. Consider the change in
er , which is a free unit vector. Therefore, the change of the vector will come from
the change in its direction only. Figure 7.39 illustrates er (and eθ ) being rotated
′
by the small angle dθ to a new vector er . Difference between these two vectors
is the vector der :
der = |er | dθeθ = dθeθ
Taking the limit as time duration approaches zero, we have
Figure 7.39: Change in er and eθ by the infinitesimal rotation ([3], pp. 68)
distance. Hence it shows how fast the position vector r contracts or stretches.
About the θ-component, vθ = r θ̇ will be non-zero if the position vector is ro-
tating. Therefore vθ indicates the change of the position vector’s direction. The
magnitude of the velocity is readily calculated due to the orthogonality of the
axes: q
v = vr2 + vθ2
Differentiating the velocity to determine the acceleration and applying the
expressions for the time rate of change of the unit vectors, we have
Figure 7.40: Infinitesimal change in the velocity and the relationship with the
acceleration represented in r-θ coordinate system ([3], pp. 69)
must neither be null for this Coriolis acceleration to be non-zero. The magnitude
of the acceleration is readily calculated due to the orthogonality of the axes:
q
a = a2r + a2θ
One final note is that ar 6= v̇r and aθ 6= v̇θ . This is because we must also
account for change in the direction of both vr and vθ , which are ṙ θ̇ along eθ and
−r θ̇2 along er , respectively.
Example 7.20 ([3], Prob. 2/147) The rocket is fired vertically and tracked by
the radar shown. When θ reaches 60◦ , other corresponding measurements give
the values r = 9 km, r̈ = 21 m/s2 , and θ̇ = 0.02 rad/s. Calculate the magnitudes
of the velocity and acceleration of the rocket at this position.
Solution: At the position where θ = 60◦ , the r-θ frame may be set up as
shown in fig. 7.42. From the problem statement, the velocity and acceleration
vectors point vertically upward. Therefore we may relate their geometric
projections to the kinematic relationships as
Example 7.21 ([3], Prob. 2/151) Link AB rotates through a limited range of
the angle β, and its end A causes the slotted link AC to rotate also. For the
instant represented where β = 60◦ and β̇ = 0.6 rad/s constant, determine the
corresponding values of ṙ, r̈, θ̇, and θ̈.
Hence,
r̈ = −0.0135 m/s2 , θ̈ = 2.31 × 10−4 rad/s2
Example 7.22 ([4], Prob. 2/150) The slotted arm OA forces the small pin to
move in the fixed spiral guide defined by r = Kθ. Arm OA starts from rest
at θ = π/4 and has a constant counterclockwise angular acceleration θ̈ = α.
Determine the magnitude of the acceleration of the pin when θ = 3π/4.
Solution: First of all, the pin P is constrained to move in the fixed spiral
guide. This implies that the path of the pin described in r-θ coordinates shall be
r = Kθ
ṙ = K θ̇, r̈ = K θ̈
Some information regarding to θ and its derivatives has been given. The
motion lasts from θ = π/4 to θ = 3π/4, with the constant angular acceleration
θ̈ = α. With the motion start from rest, θo = 0. Therefore we may be able to
compute the ending angular velocity from
h i
θ̇2 = θ̇o2 + 2θ̈ (θ − θo ) θ̇2 = 2α 3π − π4 = πα
4
The r-coordinate and its derivatives at the ending may be computed from the
path equations shown above. The results are
√
r = 3Kπ/4, ṙ = K πα, r̈ = Kα
Substituting the parameter values and calculate for the magnitude, we have
a = 10.753Kα
Example 7.23 ([3], Prob. 2/160) The circular disc rotates about its center
O with a constant angular velocity ω = θ̇ and carries the two spring-loaded
plungers shown. The distance b that each plunger protudes from the rim
of the disc varies according to b = bo sin 2πnt, where bo is the maximum
protusion, n is the constant frequency of oscillation of the plungers in the radial
slots, and t is the time. Determine the maximum magnitudes of the r- and θ-
components of the acceleration of the ends A of the plungers during their motion.
r = ro + b = ro + bo sin 2πnt
For θ-coordinate, the disc is rotating with constant angular velocity. There-
fore,
θ̇ = ω, θ̈ = 0
Substitute the above parameters into each acceleration component’s expres-
sion:
ar = r̈ − r θ̇2 = −ro ω 2 − 4π 2 n2 + ω 2 bo sin 2πnt
|ar |max = ro ω 2 + 4π 2 n2 + ω 2 bo
Example 7.24 ([3], Prob. 2/164) The small block P starts from rest at
time t = 0 at point A and moves up the incline with constant acceleration a.
Determine ṙ and θ̇ as a function of time.
For determining θ̇, one may start from the following relation
y
tan θ =
x
Differentiating the equation,
xẏ − y ẋ
θ̇ sec2 θ =
x2
Substituting the following relation
1 r2
sec2 θ = =
cos2 θ x2
Chulalongkorn University Phongsaen PITAKWATCHARA
Example 7.25 ([4], Prob. 2/159) The block P slides on the surface shown with
constant speed v = 0.6 m/s and passes point O at time t = 0. If R = 1.2 m,
determine the following quantities at time t = 2 (1 + π/3): r, θ, ṙ, θ̇, r̈, and θ̈.
Solution: The block slides with constant speed. Hence the distance trav-
eled after t = 2 (1 + π/3) sec is
Since the length of the horizontal floor is 1.2 m, the block moves up the
quarter guide by 0.4π m, for which the corresponding position of the block may
be located with the angle
0.4π π
θ= = = 60◦
1.2 3
measured from the vertical downward line in counter-clockwise direction.
To determine kinematical parameters of (r, θ) coordinates, we need to start
from the problem’s geometry. Referring to fig. 7.48, at this position,
With these (x, y) intermediate coordinates, the polar coordinate description may
be calculated with ease as
p y
r = x2 + y 2 = 2.318 m, θ = tan−1 = 15◦
x
Chulalongkorn University Phongsaen PITAKWATCHARA
For the velocity, we may employ the n-t description due to the information
given and the circular path with known geometry. At θ = 60◦ , the constant
velocity of 0.6 m/s has the direction as depicted in fig. 7.49. With the overlaying
of the n-t and r-θ coordinate frames, we may project the velocity onto r- and
θ-directions, and recall the formula to solve for the derivative parameters as
follow.
This may be decomposed into r- and θ-direction for which r̈ and θ̈ may
be resolved.
h i
ar = r̈ − r θ̇ 2
−0.3 cos 45 = r̈ − 1.2 × 0.1832 , r̈ = −0.134 m/s2
h i
aθ = r θ̈ + 2ṙ θ̇ 0.3 sin 45 = 1.2θ̈ + 2 × 0.424 × 0.183, θ̈ = 0.0245 rad/s2
Example 7.26 ([4], Prob. 2/160) The slotted arm OA oscillates about O
within the limits shown and drives the crank CP through the pin P . For an
interval of the motion, θ̇ = K, a constant. Determine the magnitude of the
corresponding total acceleration of P for any value of θ within the range for
which θ̇ = K. Use polar coordinates r and θ. Show that the magnitudes of the
velocity and acceleration of P in its circular path are constant.
Solution: The analysis will be start from the underlying geometry. The r-
and θ-coordinates must always form the isosceles triangle OCP as drawn in
fig. 7.51. Therefore β = 2θ. It is given that the slotted arm OA rotates with
the constant angular velocity θ̇ = K. Consequently, the change of β may be
determined as
β̇ = 2K, β̈ = 0
From the mechanism constraint, pin P moves along the circular path centering
at C with the radius b and angle β. Therefore its velocity and acceleration may
be determined as
v = ρβ̇ = 2bK = constant
an = v 2 /ρ = 4bK 2 , at = v̇ = 0
a = 4bK 2 = constant
Example 7.27 ([3], Prob. 2/169) The earth satellite has a velocity v = 17, 970
km/h as it passes the end of the semiminor axis at A. Gravitational attraction
produces an acceleration a = ar = −1.556m/s2 as calculated from the gravi-
tational law. For this position, calculate the rate v̇ at which the speed of the
satellite is changing and the quantity r̈.
Solution: The velocity must be tangent to the trajectory. With the satel-
lite position at one end of the vertical semiminor axis, the velocity is in
horizontal direction. Also, the n-t coordinate frame may be set up as depicted
in fig. 7.53. Because the only acceleration of the satellite comes from the
gravitational attraction, its direction points toward the earth along the radial
line.
As shown in fig. 7.53, we may expressed both the velocity and the acceleration
using r-θ coordinate system. Therefore, the rate v̇ may be determined from the
tangential acceleration:
One may conjecture that the rate r̈ may be determined using the expression of
the radial acceleration. However, the rate θ̇ in the formula is not known yet. It
may be found from the velocity relationship:
Substituting its value into the following relation, r̈ may now be solved.
rotating rate of the crank AC, the acceleration points along the normal axis
toward the center of rotation C.
Since the velocity and acceleration analysis calls for the position informa-
tion as well, we will first employ the geometric constraint that the piston-crank
mechanism forms the triangular loop OAC. Using the cosine law, we have
Example 7.29 ([4], Prob. 2/166) If the slotted arm is revolving CCW at the
constant rate of 40 rev/min and the cam is revolving clockwise at the constant
rate of 30 rev/min, determine the magnitude of the acceleration of the center of
the roller A when the cam and arm are in the relative position for which θ = 30◦ .
The limacon has the dimensions b = 100 mm and c = 75 mm.
Figure 7.57: Absolute and relative motion of A with respect to B using the pure
translating moving coordinate frame {xy} ([3], pp. 91)
Similar arguments, that shall be skipped, may be used to derive for the relative
acceleration equation:
aA = aB + aA/B (7.32)
and
aA/B = ẍi + ÿj (7.33)
It is trivial to show for the relative motion with the pure translating coordinate
frame that
r B/A = −r A/B , vB/A = −v A/B , aB/A = −aA/B (7.34)
Lastly, if the pure translating moving frame has constant velocity, aB = 0. From
the relative acceleration equation, we can conclude that aA = aA/B . Implication
of this result is that the determination of the absolute acceleration can also be
made on the inertial frame, a pure translating frame that has no acceleration.
Consequently, Newton’s 2nd law of motion holds in the inertial as well as in the
fixed reference frame.
Example 7.30 ([3], Prob. 2/195) The car A has a forward speed of 18 km/h
and is accelerating at 3m/s2 . Determine the velocity and acceleration of the car
relative to observer B, who rides in a nonrotating chair on the ferris wheel. The
angular rate ω = 3 rev/min of the ferris wheel is constant.
and
aB = ρβ̇ 2 (− cos 45i − sin 45j) = −0.628i − 0.628j m/s2
on the fact that the tangential component is null due to constant angular
velocity. Consequently, the velocity and acceleration of A relative to B may be
determined directly as
vA/B = v A − v B v A/B = 5i − (2i − 2j) = 3i + 2j m/s
Example 7.31 ([3], Prob. 2/197) Hockey player A carries the puck on his stick
and moves in the direction shown with a speed vA = 4 m/s. In passing the puck
to his stationary teammate B, by what shot angle α should the direction of his
shot trail the line of sight if he launches the puck with a speed of 7 m/s relative
to himself?
Solution: The puck must be passed from the player A to B with the di-
rection being 45◦ relative to his movement. Because the player A is moving, he
will see the puck moving in the direction of (45 + α)◦ instead. See the relevant
velocity diagram corresponding to the relative velocity relationship
v P = vA + v P/A
in fig. 7.60. Using the law of sine with this triangle, the shot angle α may be
determined:
7 4
= , α = 23.8◦
sin 45 sin α
Example 7.32 ([3], Prob. 2/210) The aircraft A with radar detection equip-
ment is flying horizontally at 12 km and is increasing its speed at the rate of 1.2
m/s each second. Its radar locks onto an aircraft flying in the same direction
and in the same vertical plane at an altitude of 18 km. If A has a speed of 1000
km/h at the instant that θ = 30◦ , determine the values of r̈ and θ̈ at this same
instant if B has a constant speed of 1500 km/h.
10 10
vB/A = v B − v A v B/A = 1500 × 36
i − 1000 × 36
i = 138.89i m/s
Example 7.33 ([3], Prob. 2/211) A batter hits the baseball A with an initial
velocity of vo = 30 m/s directly toward fielder B at an angle of 30◦ to the
horizontal; the initial position of the ball is 0.9 m above the ground level. Fielder
B requires 41 sec to judge where the ball should be caught and begins moving to
that position with constant speed. Because of great experience, fielder B choose
his running speed so that he arrives at the “catch position” simultaneously with
the baseball. The catch position is the field location at which the ball altitude
is 2.1 m. Determine the velocity of the ball relative to the fielder at the instant
the catch is made.
Solution: The ball was hit and then moved freely in the air subject to the
constant gravitational acceleration
a = −gj
This implies the zero acceleration component in the horizontal direction and
v y = −14.19j m/s
t = 2.976 s
This implies, from fig. 7.62, that the catcher B must move by the dis-
tance 77.32 − 65 = 12.32 m to the right, to meet the ball, in 2.976 − 0.25 = 2.726
second; due to the idling thinking time. Hence the constant running velocity of
the catcher is
12.32
[v = s/t] vB = 2.726
= 4.52 i m/s
Consequently, the velocity of the ball relative to the catcher at the time
of catching is
vA/B = v A − v B v A/B = (25.98i − 14.19j) − 4.52i = 21.46i − 14.19j m/s
Chapter 8
Kinetics of Particles
8.1 Introduction
This chapter is a continuation of the previous one where now the kinetics of
particles will be investigated. Kinetics is the study of the relations between the
forces and the motion. As for the introductory material, we will not seriously
concern whether the forces cause the motion or the motion generates the forces.
This is the issue called causality problem.
The object of interested in this chapter will be confined to the particles.
Relatively speaking, the body whose physical dimensions are so small compared
with the radius of curvature of its trajectory may safely be treated as the particle.
There are at least three approaches for solving the kinetics problems. They
are (a) Newton-Euler’s method (b) work and energy method and (c) momentum
method. Each has its own pros and cons. Nevertheless, we will be employing
only Newton-Euler’s method. Roughly, Newton’s law governs the force and trans-
lational motion while Euler’s law relates the moment to the rotational motion.
Since the concept of the rotation does not apply for the particle object, only
Newton’s second law of motion will be mentioned in this chapter.
For most engineering problems on earth, the acceleration measured with re-
spect to the reference frame fixed to the earth’s surface may be treated as abso-
lute. Hence it may be used in Newton’s law. Sometimes, the absolute acceleration
may be determined via the use of relative motion. If the pure translating moving
reference frame is used, the (absolute) acceleration of A might be calculated from
aA = aB + aA/B
Another caveat about this empirical law is that it does not hold when the velocity
of the order of the speed of light is involved. In that case, reader is urged to refer
to advanced topic on theory of relativity.
2. Forward dynamic Conversely, the applied force are provided. Rather, the
resulting motion will be determined, whether at the instant or as a func-
tion of time.It is far more difficult to solve compared to the inverse dynamic
problems, because the type of equations are a system of differential equa-
tions. For simple form of forcing function, we may be able to determine the
closed form solution, as for the rectilinear motion problems. Unfortunately,
the force may be described as some function of mixing parameters of time,
displacement, velocity, and acceleration. Then, only the numerical solution
may be obtained.
Unconstrained motion
As the name suggested, unconstrained motion is the motion free of constraints.
Behavior of the motion is determined by the initial conditions and the explicit
external forces. A well known example would be the projectile motion. For
particle under free motion in general has three degrees of freedom. Therefore
three scalar equations of motion may be applied and integrated to obtain the
motion.
Constrained motion
In this case, motion of the particle is partially or totally constrained (and
determined) by restraining guides, for example. Its initial motion and explicit
external forces still influence the motion, of course. Hence all forces, i.e. both
applied and reactive forces (or external force and interaction force with the
constraints), that act on the particle must be accounted for in Newton’s law of
motion. The number of degrees of freedom and associated differential equations
of motion may be reduced, depending on the type of constraints.
Additionally, in several cases, other coordinate systems, such as n-t or r-θ, might
be more appropriate.
Example 8.1 ([3], Prob. 3/19) The coefficient of static friction between the
flat bed of the truck and the crate it carries is 0.30. Determine the minimum
stopping distance s that the truck can have from a speed of 70 km/h with
constant deceleration if the crate is not to slip forward.
Solution: If the crate is not to slip forward, motion of it and the truck
must be the same at all time. This implies they must have the same acceleration.
Additionally, it is required that the truck decelerates with a maximum constant
deceleration that the crate still not slipping forward. This induces the friction
force between the truck bed and the crate to reach the maximum value. That is
the crate is on the impending status.
To determine the deceleration, we draw the free body diagram of the crate
and assign the x-axis, as shown in fig. 8.2. By the Coulomb’s model of friction,
and the specified coefficient of static friction, µs , of 0.3, the static friction force
acting on the crate is
[Fs = µs N] Fs = 0.3 × mg
[v 2 = vo2 + 2a(s − so )] 10 2
0 = 70 × 36
+ 2 (−0.3g) s, s = 64.2 m
Example 8.2 ([3], Prob. 3/20) If the truck of previous problem comes to stop
from an initial forward speed of 70 km/h in a distance of 50 m with uniform
deceleration, determine whether or not the crate strikes the wall at the forward
end of the flat bed. If the crate does strike the wall, calculate its speed relative
to the truck as the impact occurs. Use the friction coefficients µs = 0.3 and
µk = 0.25.
Solution: From the result of the previous problem, we see that the stop-
ping distance of 50 m is less than the minimum value of 64.2 m for non-slipping
condition. Therefore the crate will slip. As a result, the motion of the truck and
the crate will not be the same. For the truck, its uniform deceleration is
10 2
[v 2 = vo2 + 2a(s − so )] 0 = 70 × + 2atruck × 50, atruck = −3.781 m/s2
36
10
[v = vo + a(t − to )] 0 = 70 × 36
− 3.781tstop , tstop = 5.14 s
Consider the crate. With the material properties of the mating surface, the
static and kinetic friction force may be calculated as
Assume the crate and the truck go together as before, then atruck = acrate .
From the free body diagram of the crate (fig. 8.4), applying Newton’s law to
determine the required friction force:
which is greater than the static friction. This indicates the assumed situ-
ation is not possible. The crate must then slip and the friction drops to kinetic
friction. Re-applying Newton’s law, the acceleration of the crate may then be
determined.
To determine whether the crate strikes the wall, the analysis must be based
upon the relative motion because both the truck and the crate move. Their
relative acceleration is
The crate is moving with forward acceleration relative to the truck. Hence it
will slip forward. So we need to determine whether it will strike the wall.
To do so, we will calculate the time that the crate need for traveling 3 m
relative to the truck bed. Since the relative acceleration is constant, relationship
between the displacement and the time spent is simplified to
1
s = so + vo (t − to ) + a(t2 − t2o )
2
Note that for this problem, the kinematical quantities in the formula refer to
their relative values. Therefore,
1
3= × 1.331 × t2 , tstrike = 2.123 s
2
is the time spent for the crate to hit the wall, which is less than the time for the
truck to come to stop. Consequently, the crate will strike the wall before the
truck stops. The relative speed as the impact occurs may be calculated from
Example 8.3 ([3], Prob. 3/23) If the coefficients of static and kinetic friction
between the 20-kg block A and the 100-kg cart B are both essentially the same
value of 0.50, determine the acceleration of each part for (a) P = 60 N and (b)
P = 40 N.
Solution: From our intuition, if we apply too strong pulling force, block A
will move relatively to the right of block B. However if the applied force is small
enough, both blocks will move together. In either case, the free body diagram
will be as shown in fig. 8.6. The normal force between the blocks, NA , is to
balance the weight of block A. Therefore NA = 20g N. In turn, the maximum
supplied friction force (static) would be
Fmax = 0.5NA = 98.1 N
(a) Due to the cable-pulley arrangement, the effective pulling force is 120 N
which is greater than maximum resisting friction force Fmax . Hence block A
will slip forward relative to B. The corresponding friction force is the kinetic
frition, which happens to be equal to the static value for this problem. Their
accelerations are governed by Newton’s law:
(b) Now the effective pulling force is 80 N which is less than Fmax . Hence
block A will not slip relative to B. Since they are moving together, we may treat
them as a unified block of 120 kg subject to the pulling force of 80 N. Therefore,
the acceleration becomes
Developing friction force may be determined from either the isolated free
body diagram of block A or B. If we use block A, then
Solution: When the pendulum reaches the steady state motion, it will
possess the same acceleration as of the plate, namely a upward the incline. From
the free body diagram in fig. 8.8, there are two forces, i.e. tension and gravity
force, acting on the concentrated swinging mass, with their lines of action shown.
Choosing the coordinate frame x-y aligning with its natural motion description,
Newton’s law may then be set up conveniently;
Example 8.5 ([3], Prob. 3/26) For the friction coefficients µs = 0.25 and
µk = 0.20, calculate the acceleration of each body and the tension T in the cable.
sA + 2sB + c = l
where sA and sB are the displacement of A and B measured from the center of the
pulley positively along the incline (for A) and downward (for B). The constant
c is to take care for the portions that wrap around the pulleys, for example.
Differentiating the constraint equation to come at its acceleration form;
aA + 2aB = 0
Hence the downward motion of A will be twice as fast as the upward motion of
B, as one may understand.
For kinetic analysis, we start by drawing the free body diagram of both
objects. See fig. 8.10. Since block A involves with the friction force, one may
prepare for the calculation of the normal force and maximum supplied friction
force:
aA + 2aB = 0
Example 8.6 ([3], Prob. 3/35) A bar of length l and negligible mass connects
the cart of mass M and the particle of mass m. If the cart is subject to a
constant acceleration a to the right, what is the resulting steady-state angle θ
that the freely pivoting bar makes with the vertical? Determine the net force
P (not shown) that must be applied to the cart to cause the specified acceleration.
Solution: Free body diagram of the cart and the pendulum are shown in
fig. 8.12. Because the system is to be considered when the pendulum attains the
steady angle, both travel with the same motion, i.e. the constant acceleration a
to the right. The suitable coordinate system is then the x-y coordinate frame
shown.
The applied force P may be determined from the equation of motion of the
cart:
Solution: Free body diagram of A and B are drawn in fig. 8.14 where the
tension in the cable is denoted T . Let the positive displacement of each body be
in the direction measured from the referencing bump to the body. See fig. 8.14.
Applying Newton’s law to each of them in succession, we have
300 − 3T = 35aB
independently, but are related by the cable-pulley mechanism. With the ideal
assumption of inelasticity and rolling without slipping, the cable length, l, may
be expressed as
2sA + 3sB + c = l
where c takes care of all constant lengths along the cable route. Differentiating
the relation, we obtain the acceleration relationship
2aA + 3aB = 0
that shall be used in solving for the unknowns. The answers would be
Example 8.8 ([3], Prob. 3/45) The sliders A and B are connected by a light
rigid bar and move with negligible friction in the slots, both of which lie in a
horizontal plane. For the position shown, the velocity of A is 0.4 m/s to the
right. Determine the acceleration of each slider and the force in the bar at this
instant.
sA = sB = 0.2588 m
Change in sA and sB are governed by the cosine law as
It is given that vA = 0.4 m/s. Substituting the value into the velocity constraint
equation, we have
vB = −0.4 m/s
agreeing with our intuition.
Differentiating the general velocity relationship to obtain the acceleration re-
lation which shall be incorporated with the result from the kinetic analysis:
Substitute all known kinematical parameters into the equation, we have the re-
lationship between aA and aB at this instant as
Kinetic analysis starts with the free body diagrams of both carts, which are
simple because of the assumption of light rigid bar and negligible friction. Also,
the system is oriented to lie in a horizontal plane which makes the gravity force
pointing perpendicularly to the paper. Hence it will not affect the planar motion
of the carts. Applying Newton’s law onto each cart along its traveling direction,
we have
−T cos 15 = 3aB
Example 8.9 ([4], Prob. 3/46) With the blocks initially at rest, the force P is
increased slowly from zero to 260 N. Plot the accelerations of both masses as
functions of P .
Solution: Free body diagram of both blocks are sketched in fig. 8.18. First, we
calculate the normal forces, the static, and the kinetic frictions;
P − FA − FBk = 42a
We may use these equations in determining the range of P that makes two
blocks go together. When P increases to the value of FBmax = 113.3 N, FB
will drop from FBmax = 113.3 N to FBk = 75.54 N. Block B will then start
to move, for which we may apply the above equations of motion:
FA = 35a
68.67 = 35a
Between these extremum values, we may solve for the (common) accelera-
tion directly by viewing both blocks as a new single block of 77 kg. Hence,
the internal friction force FA vanishes and the acceleration
P −75.54
[ΣF = ma] P − FBk = 77a, a = 77
Example 8.10 ([4], Prob. 3/47) The system is released from rest in the position
shown. Calculate the tension T in the cord and the acceleration a of the 30-kg
block. The small pulley attached to the block has negligible mass and friction.
(Suggestion: First establish the kinematic relationship between the accelerations
of the two bodies.)
Solution: We see that motion of both blocks are related by the cable-pulley
mechanism. Definitions of kinematic parameters are illustrated in fig. 8.21. The
sliding block displacement x causes the change in the hypotenuse b of the right
triangle. In turn, this will make the hanging block to move, denoted as the
displacement y, according to the unchanged cable length. Mathematically,
b2 = c2 + x2
b+y =l
where c and l are some positive constants. Differentiating the relationship to
obtain the acceleration equations;
b̈ 4
=
ẍ 5
Next, the kinetics aspect will be analyzed. Since the system involves friction,
of which its direction depends on the sliding direction, we assume that the cylinder
block moves downward. Hence the sliding block will move to the left due to the
constraint which causes the friction force to point rightward. Based on this
assumption, free body diagram of the objects may be drawn, as depicted in
fig. 8.21.
4 4
[ΣFx = max ] F −T × 5
= 0.25N − T × 5
= 30ẍ
Now we have two effective equations for three unknowns: T , ẍ, and b̈. One
more equation comes from the kinematic constraint. Substituting the acceleration
of b and x obtained from the equations of motion into it, we may solve for the
tension;
4 b̈ (T /15 − g) × 30
= = , T = 137.9 N
5 ẍ 7.5g − 0.7T
And the acceleration of the block becomes
7.5g − 0.7T
ẍ = = −0.766 m/s2
30
analysis must be carried out even its related quantities have not been asked for.
Note that the complete results might not yet be obtained.
After that, one may draw the helpful free body diagram and set up Newton’s
law along the axes of the selected coordinate frame. For typical ones, they are
ΣFx = mẍ
ΣFy = mÿ (8.2)
As the last word, the positive sense of the force and accceleration must be
consistent.
Example 8.11 ([3], Prob. 3/58) The member OArotates about a horizontal
axis through O with a constant counterclockwise velocity ω = 3 rad/s. As it
passes the position θ = 0, a small block of mass m is placed on it at a radial
distance r = 450 mm. If the block is observed to slip at θ = 50◦ , determine the
coefficent of static friction µs between the block and the member.
Solution: While either r-θ or n-t coordinate system is equally suit for this
problem, the latter shall be adopted. Hence, the kinematic parameters right at
the time of slipping in n-t frame notation may be expressed as
ρ = 0.45 m, ρ̇ = 0
β = 50◦ , β̇ = 3 rad/s, β̈ = 0
Note that ρ̇ = 0 because the block is not yet slipping. However ρ̈ 6= 0.
Free body diagram of the block is shown in fig. 8.23. Here it is assumed
that the block slides down relative to the member and hence the friction
points upward. Setting up the equations of motion and substituting the known
kinematic parameters, we have
Example 8.12 ([3], Prob. 3/69) A 2-kg sphere S is being moved in a vertical
plane by a robotic arm. When the arm angle θ is 30◦ , its angular velocity about
a horizontal axis through O is 50 deg/s CW and its angular acceleration is
200 deg/s2 CCW. In addition, the hydraulic element is being shortened at the
constant rate of 500 mm/s. Determine the necessary minimum gripping force P
if the coefficient of static friction between the sphere and the gripping surfaces is
0.5. Compare P to the minimum gripping force Ps required to hold the sphere
in static equilibrium in the 30◦ position.
Solution: We will use r-θ coordinates in describing the motion of the sphere. The
frame is attached to the ball in the free body diagram as shown in fig. 8.25. The
unknowns are components of the friction force, Fr and Fθ . From the specified
kinematics, we may write
r = 1 m, ṙ = −0.5 m/s, r̈ = 0
π π
θ = 30◦ , θ̇ = −50 × = −0.873 rad/s, θ̈ = 200 × = 3.49 rad/s2
180 180
Setting up the equations of motion in r-θ frame,
If the robot is not moving at the same posture, the required minimum gripping
force will be generated purely just to cancel the gravity force, as depicted in
fig. 8.25. The friction force in that case may be determined as
Example 8.13 ([4], Prob. 3/78) A flatbed truck going 100 km/h rounds a
horizontal curve of 300 m radius inwardly banked at 10◦ . The coefficient of static
friction between the truck bed and the 200 kg crate it carries is 0.70. Calculate
[ΣFn = mv 2 /ρ] m 10 2
N sin 10 + F cos 10 = 300
100 × 36
Other direction that we know its motion is the vertical direction. Since it
is not moving,
N = 2021.52 N, F = 165.9 N
Hence the required friction force of 165.9 N will be provided. Therefore the crate
tends to slide up due to high speed circular motion. However, it is still to far
from sliding up. We can increase the truck speed yet the crate does not move
relative to the truck bed.
Example 8.14 ([3], Prob. 3/79) The flatbed truck starts from rest on a road
whose constant radius of curvature is 30 m and whose bank angle is 10◦ . If the
constant forward acceleration of the truck is 2 m/s2 , determine the time t after
the start of motion at which the crate on the bed begins to slide. The coefficient
of static friction between the crate and the truck bed is µs = 0.3, and the truck
motion occurs in a horizontal plane.
Solution: First, free body diagrams of the crate for the static and dy-
namic case are drawn in fig. 8.29. After applying the equilibrium condition for
the static case, we may solve for the normal and friction force:
Fs = 0.3Ns = 579.66 N
[at = v̇] v = at t = 2t
Next, applying Newton’s law along the y-, n-, and t-direction to arrive at the
equations of motion.
There are three unknowns, namely N, Fsn , and t, with two effective
equations. The third equation may be formulated by recognizing that the magni-
tude of the friction must be equal to the static value when the crate begins sliding.
2
[Fsn + Fst2 = Fs2 ] 2
Fsn + 4002 = (0.3N)2
√
Fsn = 0.09N 2 − 160000
Substitute the expression into the equilibrium equation, the normal force
N = 2076.47, 1919.24 N
When the truck enters the curve, the normal acceleration is developed. From
the free body diagram, it is necessary that N > Ns to match the positive com-
ponent of the normal inertial force perpendicular to the incline. Therefore
N = 2076.47 N > Ns
Fst = 400 N
4t2
477.55 cos 10 + 2076.47 sin 10 = 200 × , t = 5.58 sec
30
Example 8.15 ([3], Prob. 3/77) The small object is placed on the inner surface
of the conical dish at the radius shown. If the coefficient of static friction between
the object and the conical surface is 0.30, for what range of angular velocities
ω about the vertical axis will the block remain on the dish without slipping?
Assume that speed changes are made slowly so that any angular acceleration
may be neglected.
Solution: Figure 8.31 shows free body diagrams of the problem for two
cases, i.e. on the verge of slipping downward and upward for minimum and
maximum allowable angular velocity. Because the angular acceleration may be
neglected, the total static friction force (0.3N) points along the conical dish
as denoted by Fs . Their directions are oriented in a way to match with the
magnitude of the centerfugal acceleration.
With the specified kinematic parameters, we have
ρ = 0.2 m, ρ̇ = 0, ρ̈ = 0
For the minimum angular velocity, Newton’s law allows us to write
[ΣFn = mrω 2 ] 2
N sin 30 − 0.3N cos 30 = m(0.2ωmin )
[ΣFn = mrω 2 ] 2
N sin 30 + 0.3N cos 30 = m(0.2ωmax )
Therefore the angular velocity value between 3.405 < ω < 7.214 rad/s will not
make the block sliding.
Example 8.16 ([4], Prob. 3/82) The 2-kg slider fits loosely in the smooth slot
of the disk, which rotates about a vertical axis through point O. The slider is
free to move slightly along the slot before one of the wires becomes taut. If
the disk starts from rest at time t = 0 and has a constant clockwise angular
acceleration of 0.5 rad/s2 , plot the tensions in wires 1 and 2 and the magnitude N
of the force normal to the slot as functions of time t for the interval 0 ≤ t ≤ 5 sec.
with the zero initial conditions, the velocity and angle-time function become
θ̇ = 0.5t, θ = 0.25t2
h i
ΣFθ = m(r θ̈ + 2ṙ θ̇) N sin 45 − T sin 45 = 2 (0.1 × 0.5) = 0.1
Figure 8.34: Plots of tension wires and normal force acting on the slider
A, determine its speed v when it reaches B and the magnitude N of the force
exerted by the guide on the wheels just prior to reaching B. Neglect any friction
and any loss of mass of the rocket.
Solution: Free body diagram of the vehicle showing the gravity force, the
normal force, and the propulsion force in action is depicted in fig. 8.36. Because
the vehicle is moving along the constrained circular path, we decide to use n-t
coordinate frame. Newton’s law then gives the following governing equations:
The unknowns are the normal force and velocity, which may be determined
from the velocity-acceleration relationship. With the vehicle starts from rest at
A, we have
Rθ
[vdv = at ds] v 2 /2 = a (rdθ)
0 t
N = 3mg sin θ + 2T θ
Nθ= π2 = 3mg + T π
Example 8.18 ([3], Prob. 3/100) A hollow tube rotates about the horizontal
axis through point O with constant angular velocity ωo . A particle of mass m is
θ̇ = ωo , θ̈ = 0
Substitute the angle function into it, the differential equation of r(t) is
obtained:
r̈ − ωo2r = g sin ωo t
General solution r(t) is the sum of the particular and homogeneous solution.
r(t) = rp + rh
Since the forcing function is sinusoidal function, it is known that the forced re-
sponse of such linear system will also be sinusoidal function of the same frequency.
It is obvious that a family of the particular solution has the form
rp (t) = Cg sin ωo t
Substitute the solution into the differential equation, parameter C may be deter-
mined.
1
−Cωo2 g sin ωo t − Cωo2 g sin ωo t = g sin ωo t, C = − 2
2ωo
Homogeneous solution is the solution of the equation when the forcing function
is null. For this problem, we look for the solution, rh (t), of
r̈ − ωo2 r = 0
rh (t) = Aest
where the coefficient A and B are determined from the initial conditions:
r(0) = 0 = A + B
and
g
ṙ(0) = 0 = Aωo − Bωo −
2ωo
Solving these simultaneous equations for A and B:
g g
A= , B = −
4ωo2 4ωo2
Substitute back into the general solution, we may determine the radial parameter
as a function of time;
g
r(t) = (sinh θ − sin θ)
2ωo2
Example 8.19 ([4], Prob. 3/98) The small pendulum of mass m is suspended
from a trolley that runs on a horizontal rail. The trolley and pendulum are
initially at rest with θ = 0. If the trolley is given a constant acceleration a = g,
determine the maximum angle θmax through which the pendulum swings. Also
find the tension T in the cord in terms of θ.
When the angle reaches the maximum or minimum value, its angular velocity
will be zero. Applying this fact to the above equation, the necessary condition
for being at the apex is
sin θ + cos θ = 1
For the maximum value, the equation yields θmax = π2 .
Tension T in the cord may then be determined directly from solving the
equation of motion in the normal direction:
Substitute the velocity expression into the equation, the tension as a func-
tion of the swinging angle is
Example 8.20 ([4], Prob. 3/100) A small object is released from rest at A
and slides with friction down the circular path. If the coefficient of friction
is 0.2, determine the velocity of the object as it passes B. (Hint: Write the
equations of motion in the n- and t- directions, eliminate N, and substitute
vdv = at rdθ. The resulting equation is a linear nonhomogeneous differential
equation of the form dy/dx+f (x) y = g (x), the solution of which is well known.)
Solution: Newton’s law may be used to solve for the object’s motion, known as
the inverse dynamic problem. Free body diagram ofthe sliding mass is drawn
in fig. 8.42. Equations of motion along the normal and tangential direction are
then
d(θ̇2 ) 2
+ 0.4(θ̇2 ) = g (cos θ − 0.2 sin θ)
dθ 3
up (θ) is the forced response of a sinusoidal function: 32 g(cos θ − 0.2 sin θ). Since
the system is linear, the response will also be sinusoidal function of the same
frequency. Therefore a family of the particular solution has the form
Substitute the solution into the differential equation, parameters A and B may
be determined.
2
−A sin θ + B cos θ + 0.4(A cos θ + B sin θ) = g(cos θ − 0.2 sin θ)
3
1.2 2 0.48
A= g, B = − g
3.48 3 3.48
Homogeneous solution is the solution of the equation when the forcing function
is null. For this problem, we look for the solution, uh (θ), of
du
+ 0.4u = 0
dθ
which has the well known form of the exponential function:
uh (θ) = Cesθ
uh (θ) = Ce−0.4θ
Example 8.21 ([3], Prob. 3/101) A small collar of mass m is given an initial
velocity of magnitude vo on the horizontal circular track fabricated from a
slender rod. If the coefficient of kinetic friction is µk , determine the distance
traveled before the collar comes to rest. (Hint: Recognize that the friction force
depends on the net normal force.)
Solution: Free body diagram of the collar, depicted in fig. 8.44, reveals
the forces acting on it. Note that the normal force is decomposed to Nv and Nh
according to the natural coordinate frame n-t used. The friction force then lies
along the t-axis. Setting up the equations of motion of the moving collar along
y-, n-, and t-axes:
[ΣFy = 0] Nv = mg
[ΣFn = man ] Nh = mv 2 /r
p
[ΣFt = mat ] −F = −µk Nv2 + Nh2 = mat
µk
p
[vdv = at ds] vdv = − mr r 2 m2 g 2 + m2 v 4 ds
Integrating the differential relation during the motion of the collar and
assuming the distance traveled be s, we have
Z 0 Z s
−r 2
p d(v ) = ds
vo 2µk r 2 g 2 + (v 2 )2 0
p !
r vo2 + vo4 + r 2 g 2
s= ln
2µk rg
Example 8.22 ([4], Prob. 3/101) The slotted arm OB rotates in a horizontal
plane about point O of the fixed circular cam with constant angular velocity
θ̇ = 15 rad/s. The spring has a stiffness of 5 kN/m and is uncompressed when
θ = 0. The smooth roller A has a mass of 0.5 kg. Determine the normal force N
that the cam exerts on A and also the force R exerted on A by the sides of the
slot when θ = 45◦ . All surfaces are smooth. Neglect the small diameter of the
roller.
have
0.12 + r 2 + 0.2r cos θ = 0.22
Differentiate the equation to derive the velocity and acceleration constraints:
2ṙ 2 + 2rr̈ + 0.2r̈ cos θ − 0.4ṙ θ̇ sin θ − 0.2r θ̈ sin θ − 0.2r θ̇2 cos θ = 0
With the given motion of the circular cam, θ = π4 , θ̇ = 15 rad/s, and θ̈ = 0.
By the constraints, we may solve for the remaining parameters.
Kinetics analysis starts from the free body diagram showing all forces acting
on the roller. The angle the normal force made with the r-axis, β, may be
determined by applying the law of sine to the hidden triangle:
0.2 0.1
= , β = 20.7◦
sin 135 sin β
Spring force may be determined alternatively from the compressed length. If
the distance from O to the roller is r, at any instant the spring is compressed by
r − 0.1. The compressive force F then becomes
F = 5000 × (r − 0.1)
Example 8.23 ([4], Prob. 3/102) The small cart is nudged with negligible
velocity from its horizontal position at A onto the parabolic path that lies in a
vertical plane. Neglect friction and show that the cart maintains contact with
the path for all values of k.
Solution: To show that the cart maintains contact with the path is equivalent to
show that the reaction force (normal force solely in this case) must be greater
than zero. This inspires us to use the n-t coordinate frame since the normal force
N will align with the n-axis at all time. Drawing the free body diagram of the
cart is shown in fig. 8.48 where the only forces acting upon it is the gravity and
the normal force. Formulation of its equations of motion follws straightforwardly;
2
[ΣFn = man ] −N + mg cos θ = m vρ
dy d2 y
= tan θ = 2kx, = 2k
dx dx2
Hence by the differential calculus,
3/2
[1+y′ 2 ] (1+4k2 x2 )
3/2
ρ= y ′′ ρ= 2k
According to the n-t frame notation, t-axis is coincident with the tangent of
the path. Hence, the angle θ that it made with the horizontal reference direction
may be determined as
−1 dy
θ = tan
dx
Nevertheless, cos θ may be evaluated directly if one recognize the trigonometry
identity:
√
[1 + tan2 θ = sec2 θ] cos θ = 1/ 1 + 4k 2 x2
v 2 = 2gy = 2kgx2
Chapter 9
9.1 Introduction
This chapter studies the plane kinematics of rigid bodies. As usual, we will
not consider how the motion is achieved in kinematics. This is the subject of
kinetics. However rigid bodies motion differ from the particles motion in that
the angular motion must be taken care of properly. Typical plane motion of
rigid bodies are shown in fig. 9.1. Knowledge of the rigid bodies motion are
important in designing a mechanism to perform the desired motion. Kinematics
is the prerequisite to Kinetics. Together, they will help us determine the motion
resulting from the applied force. Now some important notions will be reviewed
before the concrete analysis.
Figure 9.1: Common plane motion of rigid bodies ([3], pp. 333)
Rigid body is a system of particles for which the distance between the particles
remain unchanged. Thus there will be no change in the position vector of any
particle measured from the body-fixed coordinate system.
Translation is the motion in which every line in the body remains paral-
lel to its original position at all time. In other words, there is no rotation of any
line in the body. The motion of the body is completely specified by the motion
of any point in the body, since all points have the same motion.
Rotation is the motion in which all particles move in circular paths about
the axis of rotation. All lines in the body which are perpendicular to the axis
of rotation rotate through the same angle in the same time. Therefore circular
motion of a point helps describe the rotating motion.
There are two main approaches for the kinematics analysis. The first approach
employs the geometry of the problem at hand for direct calculation of the absolute
displacement, velocity, and acceleration. The second method, which is easier for
most of the time, use the principle of relative motion. Later sections contain
details and examples of both methods.
9.2 Rotation
Rotation of a rigid body is described by its angular motion, which is dictated
by the change in the angular position (specified by the angle θ measured from
any fixed line) of any line attached to the body. Figure 9.2 shows the rotation
motion of a rigid body. θ1 and θ2 are the angles of line 1 and 2 measured from
any fixed reference line.
Their difference is the angle β of which its value does not matter. Its conse-
quences can be written as
θ2 = θ1 + β
∆θ2 = ∆θ1
θ˙2 = θ˙1
θ¨2 = θ¨1
ω = θ̇ (9.1)
α = ω̇ = θ̈ (9.2)
ω dω = α dθ or θ̇ dθ̇ = θ̈ dθ (9.3)
Note the analogies between the linear and angular motion.
Figure 9.3: Rotation about a fixed axis viewing from top ([3], pp. 335)
Figure 9.4: Rotation about a fixed axis in perspective view ([3], pp. 335)
at = v̇ = rα (9.5)
v2
an = rω 2 = (9.6)
r
In vector form (fig. 9.4), if ω is the angular velocity of the position vector r
(constant magnitude) which is also the angular velocity of the rigid body, motion
of A can be written as
v = ṙ = ω × r (9.7)
a = v̇ = ω × (ω × r) + α × r (9.8)
an = ω × (ω × r) (9.9)
at = α × r (9.10)
Example 9.1 ([3], Prob. 5/22) The two V-belt pulleys form an integral unit
and rotate about the fixed axis at O. At a certain instant, point A on the belt
of the smaller pulley has a velocity vA = 1.5 m/s, and point B on the belt of
the larger pulley has an acceleration aB = 45 m/s2 as shown. For this instant,
determine the magnitude of the acceleration aC of point C and sketch the vector
in your solution.
Example 9.2 ([3], Prob. 5/26) A V-belt speed-reduction drive is shown where
pulley A drives the two integral pulleys B which in turn drive pulley C. If
dω
α= dt
ωA = α1 t
Assume the pulleys roll without slipping with the belt and the belt’s elasticity
is negligible. By these kinematic constraints, the angular velocity and acceleration
of pulley B and C are
2
r1 r1
ωB = α1 t, ωC = α1 t
r2 r2
2
r1 r1
αB = α1 , αC = α1
r2 r2
Point P moves in circular path with the center of the curvature at the
center point of pulley C. Therefore its acceleration comprises of the normal and
tangential acceleration, of which the formulas are well familiar.
4
[an = rω 2 ] aPn = r2 rr21 (α1 t)2
2
[at = rα] aPt = r2 rr12 α1
Example 9.3 ([3], SP. 5/4) A wheel of radius r rolls on a flat surface without
slipping. Determine the angular motion of the wheel in terms of the linear
motion of its center O. Also determine the acceleration of a point on the rim of
the wheel as the point comes into contact with the surface on which the wheel
rolls.
Solution: The problem states the rolling without slipping condition of the
wheel motion. This is in fact important which helps identifying the kinematic
relationship. In the absolute motion approach, the kinematic relationship must
first be determined. For this specific problem with the given condition, it can be
concluded that the displacement of the center O must be equal to the arc length
along the rim of the wheel that rolls over the flat surface. With the depicted
′
figure, fig. 9.8, the distance s is equal to the arc length C A. Mathematically,
s = rθ
Differentiating the above relationship, we have the velocity and the accelera-
tion relations;
vO = rω
aO = rα
where
ẋ = r θ̇ (1 − cos θ) = vO (1 − cos θ)
ẏ = r θ̇ sin θ = vO sin θ
ẍ = v̇O (1 − cos θ) + vO θ̇ sin θ = aO (1 − cos θ) + rω 2 sin θ
ÿ = v̇O sin θ + vO θ̇ cos θ = aO sin θ + rω 2 cos θ
When point C comes to contact, θ = 0. Substitute this value into the above
equations.
ẋ = 0 ẏ = 0
ẍ = 0 ÿ = rω 2
Point C has zero velocity as expected from the constraint of rolling without
slipping. However it has the acceleration of rω 2 pointing normal off the surface.
This is in fact the tangential acceleration along the cycloidal path to take it from
rest.
There are frequent situations for which the rolling without slipping motion can
be assumed. See fig. 9.9. There we can apply these basic relations immediately.
Example 9.4 ([3], SP. 5/5) The load L is being hoisted by the pulley and cable
arrangement shown. Each cable is wrapped securely around its respective pulley
so it does not slip. The two pulleys to which L is attached are fastened together
to form a single rigid body. Calculate the velocity and acceleration of the load
L and the corresponding angular velocity ω and angular acceleration α of the
Example 9.5 ([3], Prob. 5/35) The telephone-cable reel rolls without slipping
on the horizontal surface. If point A on the cable has a velocity vA = 0.8m/s to
the right, compute the velocity of the center O and the angular velocity ω of the
reel. (Be careful not to make the mistake of assuming that the reel rolls to the
left.)
Solution: From the previous analysis, rolling without slipping implies the
velocity at the contact point is zero. Also there is no slippage of the cable at
the inner hub, which implies the velocity of the contact rim is the same as
the velocity of the wrapped cable. Observe the motion of any line on the reel,
here line OC, to determine the angular motion. From the proportional velocity
profile on the disc, fig. 9.13,
vO 0.8
= , vO = 1.2 m/s →
0.9 0.6
1.2
[vO = ωr] ω= 0.9
= 1.333 rad/s CW
Example 9.6 ([4], Prob. 5/37) The cable from the drum A turns the double
wheel B, which rolls on its hubs without slipping. Determine the angular velocity
ω and angular acceleration α of drum C for the instant when the angular velocity
and angular acceleration of A are 4rad/s and 3rad/s2 , respectively, both in the
CCW direction.
Solution: Drum A imparts its motion to the double wheel B through the
wrapped cable. In the same manner, drum C rolls CCW by the CW rotation of
the outer wheel B through the wrapped cable. Assume there is no slip of the
wrapped cable, which makes its motion the same as the tangential component of
the departed point on the drum. Therefore the motion of the upper cable is
The lower cable has the same motion as the tangential component of the
departed point from wheel B as
Figure 9.15 displays the velocity profile distribution of the double wheel B.
Finally, drum C rolls by the motion of the tensioning lower cable:
Example 9.8 ([3], Prob. 5/54) Show that the expressions v = rω and at = rα
hold for the motion of the center O of the wheel which rolls on the concave
or convex circular arc, where ω and α are the absolute angular velocity and
acceleration, respectively, of the wheel. (Hint: Follow the sample problem and
allow the wheel to roll a small distance. Be very careful to identify the correct
absolute angle through wihch the wheel turns in each case in determining its
angular velocity and angular acceleration.)
Solution: Let θ be the angular motion of the wheel while β be the angu-
lar motion of the line connecting the wheel center and the center of the curvature
arc as shown in fig. 9.18. Consider each case separately.
Concave arc:
The rolling distance can be calculated from either the wheel or the constraining
arc, which must be equal. This is the underlying geometrical relationship.
rolling distance = Rβ = r (θ + β)
s = (R − r) β
s = (R − r) β = rθ
Differentiate with respect to time to get the velocity and the tangential ac-
celeration of O:
v = ṡ = (R − r) β̇ = r θ̇ = rω
at = v̇ = (R − r) β̈ = r ω̇ = rα
Convex arc:
The rolling distance can be calculated from either the wheel or the constraining
arc, which must be equal. This is the underlying geometrical relationship.
rolling distance = Rβ = r (θ − β)
s = (R + r) β
s = (R + r) β = rθ
Differentiate with respect to time to get the velocity and the tangential ac-
celeration of O:
v = ṡ = (R + r) β̇ = r θ̇ = rω
at = v̇ = (R + r) β̈ = r ω̇ = rα
Therefore, it can be conclude that the velocity and the tangential acceleration
of the center of the wheel, O,
v = rω
at = rα
hold independent of the curvature, R, of the terrain!
Example 9.9 ([3], Prob. 5/56) The Geneva wheel is a mechanism for producing
intermittent rotation. Pin P in the integral unit of wheel A and locking plate
B engages the radial slots in wheel C thus turning wheel C one-fourth of a
revolution for each revolution of the pin. At the engagement postion shown,
θ = 45◦ . For a constant CW angular velocity ω1 = 2rad/s of wheel A, determine
the corresponding CCW angular velocity ω2 of the wheel C for θ = 20◦ . (Note
that the motion during engagement is governed by the geometry of triangle
O1 O2 P with changing θ.)
Example 9.10 ([4], Prob. 5/54) The rod AB slides through the pivoted collar
as end A moves along the slot. If A starts from rest at x = 0 and moves to the
right with a constant acceleration of 0.1 m/s2 , calculate the angular acceleration
α of AB at the instant when x = 150mm.
Example 9.11 ([3], Prob. 5/57) The punch is operated by a simple harmonic
oscillation of the pivoted sector given by θ = θ0 sin 2πt where the amplitude is
π
θ0 = 12 rad (15◦ ) and the time for one complete oscillation is 1 second. Determine
π
the acceleration of the punch when (a) θ = 0 and (b) θ = 12 .
0 = 2ẏ 2 +2y ÿ −0.28ÿ cos θ +0.28ẏ θ̇ sin θ +0.28ẏ θ̇ sin θ +0.28y θ̈ sin θ +0.28y θ̇2 cos θ
Motion of the pivoted sector is a simple harmonic oscillation described by
θ = θo sin 2πt
θ̇ = 2πθo , θ̈ = 0
Substitute these parameters into the triangular relation, motion of the punch is
y = 0.24, 0.04 m
ẏ = 0 m/s
ÿ = −0.909 m/s2
π 1
When θ = 12
, t= 4
s. Consequently,
θ̇ = 0, θ̈ = − (2π)2 θo
Substitute these parameters into the triangular relation, motion of the punch is
y = 0.2284, 0.042 m
ẏ = 0 m/s
ÿ = 0.918 m/s2
Example 9.12 ([3], Prob. 5/58) One of the most common mechanisms is the
slider-crank. Express the angular velocity ωAB and angular acceleration αAB of
the connecting rod AB in terms of the crank angle θ for a given constant crank
speed ω0 . Take ωAB and αAB to be positive counterclockwise.
l sin β = r sin θ
Figure 9.23: General motion of the rigid body: translation and rotation ([3], pp.
356)
Figure 9.24: General motion of the rigid body: translational and rotational ve-
locity ([3], pp. 357)
Figure 9.25: Velocity propagation among rigid bodies ([3], pp. 366)
Motion of the rigid body can be divided into two parts: translation and
rotation. In fig. 9.23, after the translation of the rigid body expressed by the
motion of point B, the body appears to undergo fixed axis rotation about B with
A executing circular motion as shown in fig. 9.23 (b). Hence the relationship for
circular motion describes the relative portion of A’s motion.
With B as the reference point, the total displacement of A is
frame shown in eq. 9.11. This time, point A and B are coincident points on
different rigid bodies for the instant. See fig. 9.25. Two points are at the position
of the pin in the slot. In this case, the distance between two points are not
constrained to be fixed even it happens to be zero at the moment. Therefore the
relative term vA/B can no more be determined by ω × rA/B .
3. Vector/Graphic approach
Vector algebra approach For this approach, we write each vetor in terms of i-
and j- components. Since the problem of interest is planar, there will be two
scalar equations for each vector equation. This implies at most two unknowns
can be solved per equation.
Vector/Graphic approach Here, both the vector algebra and graphical techniques
are used altogether to determine the solutions in the most convenient way. Sketch
of the vector polygon representing the vector equation is helpful. With the
sketch, we may see the convenient directions along which the projection of the
vectors yields the simple scalar component equation. Ultimately, simultaneous
equations may be avoided by a careful choice of the projecting directions.
Example 9.13 ([3], SP. 5/7) The wheel of radius r = 300 mm rolls to the right
without slipping and has a velocity v0 = 3 m/s of its center O. Calculate the
velocity of point A on the wheel for the instant represented.
Example 9.14 ([3], SP. 5/10) The power screw turns at a speed that gives the
threaded collar C a velocity of 0.25 m/s vertically down. Determine the angular
velocity of the slotted arm when θ = 30◦ .
Solution: Imagine there are two points: point A on the slotted arm and
point B on the collar. A and B coincide at θ = 30◦ . Because of the sliding
contact constraint from the slot, vA/B has the direction along the slot (away
from O). Also vA is always perpendicular to the link AO. See fig. 9.29.
vA = vB + vA/B vA = 0.25 cos 30 = 0.217 m/s
Example 9.15 ([3], Prob. 5/78) The rotation of the gear is controlled by the
horizontal motion of end A of the rack AB. If the piston rod has a constant
velocity 300 mm/s during a short interval of motion, determine the angular ve-
locity of the gear and the angular velocity of AB at the instant when x = 800 mm.
Setting up the velocity equation and draw the relevant velocity polygon
(fig. 9.31), we have
vC = vA + vC/A vC = 0.3 cos θ = ωO × 0.2, ωO = 1.45 rad/s CW
√
vC/A = 0.3 sin θ = ωAB × 0.82 − 0.22 , ωAB = 0.0968 rad/s CCW
Example 9.16 ([3], Prob. 5/83) The flywheel turns CW with a constant speed
of 600 rev/min, and the connecting rod AB slides through the pivoted collar
at C. For the position of θ = 45◦ , determine the angular velocity of AB using
the relative velocity relations. Choose a point D on AB coincident with C as a
reference point whose direction of velocity is known.
vA can be determined from either the motion of the flywheel or the motion
of link AB. From this relationship and the velocity polygon in fig. 9.33, ωAB can
be found.
vA = vD + vA/D vA = rω = 0.2 × (600 × 2π/60) = 12.566 m/s
Example 9.17 ([4], Prob. 5/85) The Geneva mechanism is shown again here.
By relative motion principle, determine the angular velocity of wheel C for
θ = 20◦ . Wheel A has a constant CW angular velocity ω1 = 2 rad/s.
Solution: Let P be the point on the wheel A at the knob, while Q be the
point on the wheel C that is coincident with P when θ = 20◦ . Their relative
velocity is constrained to be along the slot on the wheel C. Therefore
√
vP = vQ + vP/Q vP = ω × r = 2 × 0.2/ 2 = 0.283 m/s
Example 9.18 ([3], Prob. 5/91) At the instant represented, a = 150 mm and
b = 125 mm, and the distance a + b between A and C is decreasing at the rate of
0.2 m/s. Determine the common velocity v of points B and D for this instant.
Solution: Because B and D are on the same rigid body translating along
the slot, vB = vD . The velocity relationship is obvious:
vB = vA + vB/A vD = vC + vD/C
From the statement, block A and C become closer which means ȧ + ċ = −0.2.
And from the imposing mechanism, A moves to the right and B to the left.
Therefore
vA − vC = 0.2i m/s
Referring to the velocity diagram in fig. 9.37,
vB vB
+ = 0.2, vB = 0.0536j m/s
tan 22.62 tan 36.87
Example 9.19 ([3], Prob. 5/92) The wheel rolls without slipping. For the in-
stant portrayed, when O is directly under point C, link OA has a velocity v = 1.5
m/s to the right and θ = 30◦ . Determine the angular velocity ω of the slotted link.
Let P be the point at the pin on the disk and Q be the coincident point on
the slotted arm. According to the figure, their relative velocity directs along the
slot:
vP = vQ + vP/Q vQ = vP cos (15 + 23.78) = 2.26 m/s
Example 9.20 ([4], Prob. 5/88) Ends A and C of the connected links are
controlled by the vertical motion of the piston rods of the hydraulic cylinders.
For a short interval of motion, A has an upward velocity of 3 m/s, and C has a
downward velocity of 2 m/s. Determine the velocity of B for the instant when
y = 150 mm.
vB = vC + vB/C = vA + vB/A
The related velocity diagram is shown in fig. 9.40, from which, by the law of
In case the body motion is just the translation motion, the ICZV is located
at infinity along the perpendicular line to the velocity.
Example 9.22 ([3], Prob. 5/121) Determine the angular velocity ω of the ram
head AE of the rock crusher in the position for which θ = 60◦ . The crank OB
has an angular speed of 60 rev/min. When B is at the bottom of its circle, D
and E are on a horizontal line through F , and lines BD and AE are vertical.
The dimensions are OB = 100 mm, BD = 750 mm, and AE = ED = DF = 375
mm. Carefully construct the configuration graphically, and use the method of
ICZV.
Solution: Graphical method using ICZV is the best way for this problem,
which is capture in fig. 9.43. Here we explain another approach using the
vector/graphic method. From the initial and current posture, write the distance
between point O and F :
Example 9.23 ([3], Prob. 5/122) The shaft at O drives the arm OA at a
clockwise speed of 90 rev/min about the fixed bearing at O. Use the method of
ICZV to determine the rotational speed of gear B (gear teeth not shown) if
(a) ring gear D is fixed and
(b) ring gear D rotates CCW about O with a speed of 80 rev/min.
Solution: It is easily observed that the motion of gear A is influenced
by
vA = 90 a ↓
a
vC = 2vA = 2 × 90 a = ωB × , ωB = 360 rev/min CW
2
When the ring gear D rotates at 80 rev/min CCW (see fig. 9.46),
3a
vA = 90 a ↓, vD = 80 × = 120 a ↑
2
10/7 a
vC = (90 a) = ωB × , ωB = 600 rev/min CW
3/7 2
Example 9.24 ([4], Prob. 5/118) The large roller bearing rolls to the left on
its outer race with a velocity of its center O of 0.9 m/s. At the same time, the
central shaft and inner race rotate CCW with an angular speed of 240 rev/min.
Determine the angular velocity of each of the rollers.
Solution: The vertical velocity profile is imposed on fig. 9.47. From the
problem, it is specified that vO = 0.9 m/s and ωi = 240 rev/min = 8π rad/s
CCW.
From the velocity profile, the ICZV of the inner race is 0.9
8π
m lower to point
O. Therefore, velocity of the point on the roller contacting with the inner race is
0.9
v = 8π × 0.05 − = 0.3566 m/s →
8π
If we imagine the bearing as one rigid disk rolling with the same angular
velocity. That value is
0.9
ωO = = 7.2 rad/s CCW
0.125
which is also the angular velocity of the outer ring. Hence, velocity of the point
on the roller contacting with the outer race is
Consider the velocity profile across the roller. Using the similar triangle rela-
tionship,
0.18 x
= , x = 16.77 mm
0.3566 0.05 − x
Therefore, its angular velocity is
0.18
ωroller = = 10.732 rad/s CW
x
Figure 9.48: Relative acceleration between point A and B ([3], pp. 381)
where
= rω 2 = vA/B
2
aA/B n
= ω × (ω × r) , aA/B n
/r (9.19)
aA/B t
= α × r, aA/B t
= rα = v̇A/B (9.20)
3. Vector/Graphic approach
Vector algebra approach Each term in the equation is represented in the selected
coordinate system, which may be x − y, n − t, or r − θ coordinates. For
two dimensional problems, each vector equation has two independent scalar
equations. Therefore at most two unknowns can be solved per equation. This
method is appropriate for writing the computer codes to solve the problems or
when the exact analytical expression is needed.
Example 9.25 ([4], Prob. 5/138) If the wheel in each case rolls on the circular
surface without slipping, determine the acceleration of the point C on the wheel
momentarily in contact with the circular surface. The wheel has an angular
velocity ω and an angular acceleration α.
Applying these expressions to the acceleration equation for each case, we have
(a) aC = aO + aC/O
aC = (aO )t + (aO )n + aC/O t
+ aC/O n
2 2
r ω R
= rαi + j − rαi + rω 2j = rω 2 j
R−r R−r
Example 9.26 ([4], Prob. 5/144) The simplified clam-shell bucket is shown.
With the block at O considered fixed and with the constant velocity of the
control cable at C equal to 0.5 m/s, determine the angular acceleration α of the
right-hand bucket jaw when θ = 45◦ as the bucket jaws are closing.
Velocity analysis Here it may be simpler to apply the ICZV method. First,
locate the ICZV of the right jaw. Since the velocity direction of B and C are
known, ICZV is the intersection of the radial lines as shown in fig. 9.52. Let that
ICZV point be called point D. Then, for the given configuration θ = 45◦ ,
On the bucket,
where
2
= 0.2604 m/s2 2
= 0.0782 m/s2
aB/C n
= BCωBC (aB )n = BOωBO
Substitute these values in the above relationship. From the acceleration dia-
gram in fig. 9.52 along the horizontal direction, we have
0.2604 cos 22.5 + aB/C t sin 22.5 − (aB )t cos 50.345 − 0.0782 sin 50.345 = 0
Hence
Therefore
Example 9.27 ([3], Prob. 5/152) The mechanism where the flexible band F
attached to the sector at E is given a constant velocity of 4 m/s as shown. For
the instant when BD is perpendicular to OA, determine the angular acceleration
of BD.
Solution: Motion of the flexible band causes the rotation of the sector,
which then causes the motion of the link AD and BD, respectively. The sector
rotates about O by the pulling force of the band. Consequently, point E on the
sector moves in circular path about O with the velocity of 4 m/s; the velocity of
the band. Hence, at this instant,
and
vD/A /vA = 5/4, vD/A = 3.125 m/s
Hence,
Acceleration analysis Since the flexible band pulls the sector with a con-
stant velocity, (aE )t = 0. Therefore the angular acceleration of the sector is zero.
This implies the total acceleration of A, aA , is in the normal direction pointing
towards O. Consequently, aD can be determined as
[an = rω 2 ]
aD = aA + aD/A
of which their magnitudes can be determined from the lengths of the polygon as
3 4
= 11.72 m/s2
39.0625 ×
− aD/A t × − 14.0625 = 0, aD/A t
5 5
4 3
(aD )t = 11.72 m/s2
39.0625 × + aD/A t × + (aD )t − 50 = 0,
5 5
Hence,
Example 9.28 ([3], Prob. 5/156) Elements of the switching device are shown.
If the velocity v of the control rod is 0.9 m/s and is slowing down at the rate of
Velocity analysis The configuration here suggests the use of ICZV for the
analysis. Since the motion of point A and B on the linkage ABC are known,
the ICZV point of ABC is the intersection of the perpendicular lines to their
velocity as depicted in fig. 9.55.
0.9
[ω = v/r] ωAB = 0.075 cos 30
= 13.856 rad/s CCW
With the vectors in the equation projected onto the vertical direction, the
following equation can be written:
6 + 14.4 sin 30 − aA/B t sin 60 = 0
With the vectors in the equation projected onto the horizontal direction, the
following equation can be written:
aA = 14.4 cos 30 + aA/B t cos 60
Solving these two simultaneous equations,
aA/B t = 15.242 m/s2 , aA = 20.09 m/s2
Then, aC can be directly determined from the following relationship. See the
acceleration diagram in fig. 9.56.
aC = (14.4 cos 30 + 15.242 cos 60) i + (6 − 14.4 sin 30 + 15.242 sin 60) j
= 20.09i + 12.0j, aC = 23.4 m/s2
Example 9.29 ([4], Prob. 5/149) An oil pumping rig is shown in the figure.
The flexible pump rod D is fastened to the sector at E and is always vertical as
it enters the fitting below D. The link AB causes the beam BCE to oscillate as
the weighted crank OA revolves. If OA has a constant CW speed of 1 rev every
3 s, determine the acceleration of the pump rod D when the beam and the crank
OA are both in the horizontal position shown.
Solution: The flexible pump rod D and the sector E are connected such
that the motion of the sector in the tangential direction is imparted to the
rod entering and exiting vertically the hole tube. According to the motion
transmitted by the linkage OA, AB, and BC forming the four bar linkage
mechanism, the sector itself moves back and forth along the specific circular arc.
To determine the acceleration of the pumping rod, kinematic analysis along the
motion transmission from the crank OA should be perform starting from the
velocity-level analysis as follow.
Velocity analysis First, let us convert the unit of the angular velocity of
the crank OA.
1 2π
ωOA = rev/s = rad/s CW constant
3 3
At the position shown in fig. 9.57, the beam BCE and the crank OA are both
in the horizontal position, for which the velocity at A can be simply calculated
as
vA = ωOA × OA = 0.4π m/s ↑
In this case, ωAB and ωCE can be conveniently determined using the method
of ICZV. Since the direction of vA and vB are known, the location of the ICZV
point of the connecting rod AB is the intersection of the perpendicular lines to
those velocity as depicted in fig. 9.58. Of course the ICZV point of the beam
BCE is the fixed rotating point C. After calculating the exact coordinates of
the ICZV point,
The acceleration polygon is sketched in fig. 9.58. Projecting all vectors onto
the horizontal direction, the following relation holds:
2.632 − 0.046 cos 78.1 − aB/A t cos 11.9 − (aB )t sin 16.7 − 0.496 cos 16.7 = 0
Similarly, projecting all vectors onto the vertical direction, the following relation
holds:
−0.046 sin 78.1 + aB/A t sin 11.9 − (aB )t cos 16.7 + 0.496 sin 16.7 = 0
Consequently, the acceleration of the pumping rod can now be readily determined:
(aB )t
[α = at /r] αCE = √
0.92 +32
= 0.1724 rad/s2 CW
Hence,
aA/B
αAB = t
= 24.467 rad/s2 CW
AB
Using the acceleration relationship between B and D to determine the accel-
eration at D. Since the link has no rotation at this moment,
aD/B n = 0
aD = aB + aD/B aD = (−7.896 + 7.34 cos 14.4775) i + (7.34 sin 14.4775) j
aD = −0.789i + 1.835j m/s2 , aD = 1.997 m/s2
Figure 9.61: Motion relative to the rotating axes frame ([3], pp. 395)
It is evident that i and j are not constant since their directions change. Their
rate of changes with respect to time can be determined by investigating the
drawing of the unit vector i and j rotating by the angle of dθ for the elapsed time
of dt, shown in fig. 9.61. With the definition of the differentiation, the following
equations can be derived:
di
= ω× i = ωj (9.22)
dt
dj
= ω× j = −ωi (9.23)
dt
Figure 9.62: Visualization of relative velocity equation using the virtual coinci-
dent point P ([3], pp. 396)
Differentiating eq. 9.21 with respect to time and making use of the above
relations, the following relative velocity relationship results:
Another way to comprehend the relative velocity equation using the rotating
frame is to visualize the situation shown in fig. 9.62. There are two particles A
and B which may be on different rigid bodies. Imagine there is a rotating plate
on which the particle B is situated. On this plate, imagine the virtual point P
currently coincident with A. As a result, the relationship between vP and vB is
indicated by
vP = vB + ω× rP/B = vB + ω× rA/B
Since point A is on the different body, vA 6= vP . From the above relation,
it can be concluded that the velocity of A as seen from P is vrel . In fact, the
observer at any location fixed to the rotating plane will see A to be moving
with the velocity vrel . For this particular illustration, the direction of vrel is
tangent to the path (slot) fixed in the rotating plate. That is, the virtual slot
is constructed to be coincided with the trajectory seen by the observer fixed to
the rotating plate. Note that it is not the absolute path of A (which must be
measured by the fixed observer).
From the illustration in fig. 9.62, the magnitude of vrel , or the relative speed,
is ṡ. Since P is the point fixed to the moving plate that is instantaneously
coincident with A, it can be concluded that vrel = vA/P . That is, the velocity of
A seen by the rotating observer B or P is the same, and it is equal to the velocity
of A relative to P . Be caution that it is not the same as the velocity of A relative
to B.
To summarize, observers moving with different velocities (different vB ) on the
same rotating x-y frame see the target moving with the same velocity of vrel .
Observers moving with the same velocity, but are on different rotating x-y frames
(different ω), see the target moving with different velocity vrel . Non-rotating
observers will see the resultant of circular motion plus the relative velocity. The
following relative velocity equations summarize the relative motion of point A
and B from different point of views.
vA = vB + ω × r + vrel (9.25)
vA = vB + vA/B (9.26)
vA = vB + vP/B + vA/P (9.27)
vA = vP + vA/P (9.28)
This expression means the time derivative of V measured in the fixed frame
(the total time derivative) is equal to the time derivative of V as measured in the
rotating frame plus the compensation due to rotation of the reference frame.
More insight can be seen from the vector diagram in fig. 9.63. Vector V
′
changes in both direction and magnitude to V . In the figure, x-y frame changes
by rotating with the angular velocity ω while X-Y frame is fixed. During the
time interval dt, the observer in rotating x-y frame see the change in magnitude
of V, dV , plus the change in direction, V dβ, due to the relative rotation of V to
x-y. The change that the rotating observer recognized is called (dV)xy . What it
does not notice is the rotation of V induced by the rotation of x-y, V dθ. Imagine
that V is fixed to x-y. Therefore, its direction changes by the rotation of x-y
frame, which is not known to the observer rotating together.
This fundamental relationship can be used to derive the relative acceleration
equation from the relative velocity equation.
With the use of eq. 9.30, the differentiation terms above can be further exemplified
as
The last equation implies the angular acceleration seen in the rotating frame
is the absolute angular acceleration because the vector of the angular velocity
aligns with the angular velocity of the observing frame. Substituting these terms
into the acceleration relationship, the relative acceleration equation employing
the rotating frame results:
aA = aB + ω̇× rA/B + ω× ω× rA/B + 2ω× vrel + arel (9.34)
Figure 9.64: Visualization of relative acceleration equation using the virtual co-
incident point P ([3], pp. 398)
In summary, arel is the change of vrel observed from any location in the
rotating frame. This is the acceleration
of A seen by any rotating observer. It is
not equal to v̇rel = ω× vrel + arel which is just the change of vrel observed
in
the nonrotating frame. Also, it is not equal to aA/P = 2ω× vrel + arel which is
the acceleration of A seen by the coincident point P in the nonrotating frame. If
the observer is not coincident to A, he will see the resultant of normal, tangential
(portion of the relative motion which makes up the relative circular motion),
Coriolis, and relative acceleration (the rest of the relative motion by virtual slot
on the rotating frame). The following relative acceleration equations summarize
the relative motion of point A and B from different point of views.
aA = aB + ω̇× rA/B + ω× ω× rA/B + 2ω× vrel + arel (9.35)
aA = aB + aA/B (9.36)
aA = aB + aP/B + aA/P (9.37)
aA = aP + aA/P (9.38)
Example 9.31 ([3], Prob. 5/183) The crank OA revolves clockwise with a
constant angular velocity of 10 rad/s within a limited arc of its motion. For the
position θ = 30◦ , determine the angular velocity of the slotted link CB and the
acceleration of A as measured relative to the slot in CB.
2ω CB × v = 10 m/s2
rel
aA = vA2 /OA = 20 m/s2
Construct the acceleration diagram as depicted in fig. 9.66 and perform the geo-
metrical analysis, the remaining acceleration can be determined.
velocity of 2 rad/s.
Geometry analysis First, the geometric parameters at this instant are de-
termined. Consider the triangle O1 P O2 when θ = 20◦ shown in fig. 9.68.
200
√
2
sin 20
tan α = 200 , α = 35.78◦
200 − √
2
cos 20
2 2
2 200 200
P O2 = √ sin 20 + 200 − √ cos 20 , P O2 = 82.7 mm
2 2
Velocity analysis Next, the velocity relationship of the pin-coincident points A
on the wheel A and P on the wheel C is formulated.
vA = vP + vrel
Draw the velocity diagram and perform the geometrical analysis to unveil vP and
vrel .
Since the path of A seen on the rotating wheel C is along the slot, the
direction of arel is parallel to the slot. Considering the unknowns in the above
equation, they are the magnitudes of ω̇ 2 × rA/O and arel for the two scalar
equations to solve. See the acceleration diagram in fig. 9.68. The magnitudes of
the known terms are
√
aO2 = 2 × 2 × 0.2/ 2 = 0.5656 m/s2
2ω 2 × v = 0.9 m/s2
rel
recognizing the fact that point A is moving on the circular path at constant
angular velocity.
With the acceleration diagram, projecting the sides of the vector polygon onto
the ω̇ 2 × rA/O direction, the following relationship can be written:
ω̇ 2 × rA/O = 0.9 + 0.5656 cos 34.22
From the appendix and assuming the axis of the rotation of the earth is normal
to the equilateral plane, ωe = 0.7292 × 10−4 rad/s constant CCW. Therefore,
To determine the velocity of the space shuttle, recall the circular-orbit shuttle
has the normal acceleration towards the center of the earth of which its value is
2
vA2
R
(aA )n = g =
R+h R+h
Consulting the appendix for the accurate value of the gravitational constant,
g = 9.814 m/s2 . Thus
vA = 7766.79 m/s ←
Applying the relative velocity equation above, the velocity seen by the ob-
server B, vrel , can now be determined.
Let us first count the unknowns in the equation. From the velocity analysis,
(aA )n has been determined. However, the magnitude of (aA )t is the unknown.
The other unknown is the acceleration seen by the observer, arel . Hence there
are three scalar unknowns in this two-dimensional vector equation for which not
all of them may be solved.
Assume that the shuttle orbits with the constant velocity magnitude, (aA )t =
0 and so aA = (aA )n . As a result, only arel provides two scalar unknowns which
can be solved from the equation. Computing the values of each term in the
acceleration equation,
2
6378 2
−9.81 j = − 0.7292 × 10−4 × 6378 × 103 j
6378 + 240
2
− 0.7292 × 10−4 × 240 × 103 j − 2 × 0.7292 × 10−4 × 7284.205j + arel
Here the angular velocity of the earth is presumed to be constant as well which
makes aB = (aB )n . Solving the equation for the acceleration of the space shuttle
seen by the observer B,
arel = 8.018 m/s2 ↓
Example 9.34 ([3], Prob. 5/179) The figure shows the vanes of a centrifugal
pump impeller which turns with a constant clockwise speed of 200 rev/min. The
fluid particles are observed to have an absolute velocity whose component in the
r-direction is 3 m/s at discharge from the vane. Furthermore, the magnitude of
the velocity of the particles measured relative to the vane is increasing at the
rate of 24 m/s2 just before they leave the vane. Determine the magnitude of the
total acceleration of a fluid particle an instant before it leaves the impeller. The
radius of curvature ρ of the vane at its end is 200 mm.
Velocity analysis
vA = vP + ω× rA/P + vrel
With the help of the velocity diagram shown in fig. 9.71 and expressing the
2
√
20π 20π
aA = −0.15 × × 3 2 (cos 45er + sin 45eθ )
er + 2 ×
3 3
√ 2
3 2
+24 (cos 45er − sin 45eθ ) + (− cos 45er − sin 45eθ )
0.2
= 13.187er + 45.04eθ
aA = 46.93 m/s2
Solution: Rotation of the screw causes the relative translation of the col-
lar A. With the given thread specification, three turns advances the collar by
1 cm. Let A be the point on the collar and P be the coincident point on the
screw. Since the screw rotates at a constant speed 100 rev/min,
100 1
cm/min =
vA/P = m/s
3 180
in the direction along the screw towards C.
Velocity analysis Using the velocity relationship and the velocity diagram
in fig. 9.73, the absolute velocity of the collar and the links’ angular velocities
can be determined.
vA = vP + vrel
1/180
vA = = ωOA × 0.2, ωOA = 0.032 rad/s CCW
cos 29.745
vP = vrel tan 29.745 = 0.0032 m/s, ωP C = vP /CP = 0.00788 rad/s CCW
rA/P = 0,
aA = aP + 2ω P C × vrel + arel
The corresponding acceleration diagram is drawn in fig. 9.73. For this case,
arel = 0 because the screw turns at a constant rate. The centrifugal acceleration
component of aA and aP can be calculated as
[an = rω 2 ] 2
(aA )n = 0.2 × ωOA = 204.8 × 10−6 m/s2
(aA )t = ω̇OA ×0.2 = −25.41×10−6 cos 29.745◦ +(87.556 + 133.8)×10−6 sin 29.745◦
Example 9.36 ([3], Prob. 5/185) Determine the angular acceleration of link
EC in the position shown, where ω = β̇ = 2 rad/s and β̈ = 6rad/s2 when
θ = β = 60◦ . Pin A is fixed to link EC. The circular slot in link DO has a
radius of curvature of 150 mm. In the position shown, the tangent to the slot at
the point of contact is parallel to AO.
rA/P = 0
Velocity analysis With the current geometry of the mechanism, the velocity
relationship and its diagram (fig. 9.75) of point A and P can be determined as
vA = vP + vrel
Bibliography
BIBLIOGRAPHY 380