Mechanics Notes
Mechanics Notes
net/publication/342510225
Mechanics II
CITATIONS READS
0 1,723
1 author:
Babar Ahmad
COMSATS Institute of Information Technology Islamabad, Pakistan
23 PUBLICATIONS 30 CITATIONS
SEE PROFILE
All content following this page was uploaded by Babar Ahmad on 28 June 2020.
ANALYTIC DYNAMICS I
1
MECHANICS II
ANALYTIC DYNAMICS I
All rights are reserved. No part of this book may be reproduced, or transmitted, in any form
or by any means – including, but not limited to electronic, mechanical, photocopying, recording,
or, otherwise or used for any commercial purpose what so ever without the prior written permission
of the writer and, if author considers necessary, formal license agreement with author may be
executed.
Edition: 2020
Price
Disclaimer: The author has used his best efforts for this publication to meet the quality
standards, but does not assume, and hereby disclaims for any loss or damage caused by the errors
or omissions in this publication, where such errors or emissions result from negligence, accident,
or any other cause.
To My Parents
Mr. & Mrs. Rana Muhammad Hanif (late)
(May ALLAH bless them)
PREFACE
Mechanics is one of the most important course in maximum disciplines of science and engineering.
No matter what your interest in science or engineering, mechanics will be important for you.
Mechanics is a branch of physics which deals with the bodies at rest and in motion. During
the early modern period, scientists such as Galileo, Kepler, and Newton laid the foundation for
what is now known as classical mechanics. Hence there is an extensive use of mathematics in its
foundation.
Mechanics is core course for undergraduate Mathematics, Physics and many engineering
disciplines. It appears under different names as Analytical/Classical Mechanics, Theoretical
Mechanics, Mechanics I, Mechanics II, Mechanics III, Analytical Dynamics.
This textbook is designed to support teaching activities in Theoretical Mechanics specially
Dynamics. It covers the contents of “Mechanics” for many undergraduate science and engineering
programs. It presents simply and clearly the main theoretical aspects of mechanics.
It is assumed that the students have completed their courses in Calculus, Linear Algebra and
Differential Equations. This book also lay the foundations for further studies in physics, physical
sciences, and engineering.
For each concept a number books, documents and lecture notes are consulted. I wish to express
my gratitude to the authors of such works.
In Chapter 1, preliminary theory of dynamics of rigid body is given. Newton’s laws and some
types of motion are also given in this chapter. Chapter 2 is about one dimensional motion. In this
chapter, the kinematics of body is discussed by using graphical method, differentiation and
integration. In chapter 3, the motion is also one dimensional, but the particle is restricted to move
under gravity. This motion is discussed without and with air resistance. In chapter 4, the kinematics
of a body in two and three dimensional coordinate systems is presented.
From Chapter 5, we will learn angular kinematics. In this chapter we will learn kinematics of a
particle in polar coordinates. In chapter 6, Simple Harmonic Motion is discussed under natural,
forced and damping aspects. Two dimensional projectile motion (without and with air resistance)
is presented in chapter 7. Chapter 8 starts with Kepler’s postulates to discuss the motion under
Central Force. Some useful information about planetary motion is also available. Chapter 9 is
about Small Oscillation. The stability of horizontally and vertically modulated pendulum under
various forces is discussed. Chapter 10 is about rotational dynamics. Chapter 10 contains Rotation
of rigid bodies in two space and three space. The concepts of rotational matrix, angle of rotation
and axis of rotation are also discussed in this chapter. Euler’s theorem and Chasle’s theorem are
also given in this chapter. This chapter also contains kinematics of a body in cylindrical and
spherical coordinate systems.
In a book of this concept, level and size, there may be a possibility that some misprint might have
remained uncorrected. If you find such misprints or want to give some suggestions for its
improvement, please write me at: babar.sms@gmail.com
Dr. Babar Ahmad
Islamabad, Pakistan
June, 2020
vi
Table of Contents
2 Rectilinear Motion 9
2.1 Rectilinear Motion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1.1 Position Vector . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.1.2 Displacement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.1.3 Instantaneous Velocity . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.1.4 Instantaneous Acceleration . . . . . . . . . . . . . . . . . . . . . . . 11
2.2 Constant and Variable Quantities . . . . . . . . . . . . . . . . . . . . . . . 12
2.3 Graphical Methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.3.1 Time-Displacement Curve . . . . . . . . . . . . . . . . . . . . . . . 12
2.3.2 Time-Velocity Curve . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.3.3 Time-Acceleration Curve . . . . . . . . . . . . . . . . . . . . . . . . 15
vii
viii CONTENTS
5 Angular Motion 81
5.1 Angular Kinematics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
5.1.1 Angular Displacement . . . . . . . . . . . . . . . . . . . . . . . . . . 81
5.1.2 Angular Velocity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
5.1.3 Angular Acceleration . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
5.2 Motion in Polar Coordinates or
Circular Motion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
5.2.1 Uniform Circular Motion . . . . . . . . . . . . . . . . . . . . . . . . 84
5.2.2 Position Vector . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
5.2.3 Velocity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86
5.2.4 Acceleration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
5.2.5 Centripetal Force . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
5.2.6 Relation Between Linear and Angular Acceleration (in magnitude) . 89
5.3 Motion in Radial and Transverse Plane . . . . . . . . . . . . . . . . . . . . 91
5.3.1 Position Vector . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
5.3.2 Radial and Transverse Components of Velocity . . . . . . . . . . . . 92
CONTENTS ix
Bibliography 322
xii CONTENTS
Chapter 1
In this chapter some basic concepts of Mechanics are presented for the revision of elementary
knowledge of the students. These concepts are essential part of mechanics. Some types of
Motion are defined, their detail will be in next chapters. Newton’s three laws of motion
and Newton’s gravitational law are also mentioned.
1. Kinematics is the study of motion without reference to the forces which cause motion.
It is the geometric theory of motion.
2. Kinetics is the study of motion with reference to the forces which cause motion.
1
2 1 Basic Review of Mechanics
A man standing on ground is observing the flight of an aeroplane moving away from him.
At first instant he can see its shape very clear. At next instant, the shape becomes smaller
and next becoming smaller and smaller. At last instant he sees a dot and then nothing.
This dot can be regarded as particle (see Fig. 1.1).
Body: An aggregate of particles is known as body.
Kinds of bodies: There are various types of bodies e.g rigid, elastic, plastic, elasto plastic,
etc
The bodies are of two categories:
a) Rigid bodies
b) Deformable bodies
Rigid body A body is considered rigid when the relative movement between its parts is
negligible.
In Mechanics we will assume the bodies to be perfectly rigid, no deformation. This is never
true in the real world, everything deforms a little when a load is applied. These deforma-
tions are small and will not significantly affect the conditions of equilibrium or motion, so
we will neglect the deformations.
~r(t) = r(t)r̂
If position function of a particle is given, then differentiation is used to define the notions
of instantaneous velocity and acceleration, and if acceleration function is given, then inte-
gration is used to define the notions of instantaneous velocity and position. First we define
the relations of instantaneous velocity and acceleration for a particle in rectilinear motion
in the form of derivatives.
d 0
~v (t) = ~r(t) = ~r (1.3.2)
dt
The velocity vector ~v always points in the direction of motion. The magnitude of the ve-
locity, v = |~v | is known as the speed.
Units commonly used to express velocity and speed are m/s and f t/sec.
Thus, if a particle in rectilinear motion has velocity function ~v (t), then its instantaneous
acceleration at time (acceleration function) is defined as
∆~v
~a(t) = lim
∆t→0 ∆t
d
= ~v (t) (1.3.3)
dt
If the velocity changes in either magnitude or direction (or both), the particle must have
an acceleration.
Using (4.1.1), the acceleration function in terms of the position function is
d2 00
~a(t) = 2
~r(t) = ~r (1.3.4)
dt
The acceleration vector ~a can point anywhere. Units commonly used to express acceleration
are m/s2 or f t/sec2
1.3.4 Displacement
If the velocity vector ~v (t) of a moving particle over the time interval [t1 , t2 ] is given, then
the displacement ∆~r = ~r(t2 ) − ~r(t1 ) or simply ~r(t) is
Zt2
~r(t) = ~v (t)dt (1.3.5)
t1
1.4.1 Translation
A motion by which a body shifts from one point in space to another in such a way that
the straight lines joining the initial and final positions of each of the points of the body
are parallel. i.e. every point of the body moves an equal distance. If the path of motion
is straight line, the motion is said a rectilinear motion. The line might be an x − axis,
a y − axis, or a coordinate line inclined at some angle. If the paths are curved lines, the
motion is a curvilinear motion. The motion of a bullet fired from a gun is the example
of such motion.
1.4.2 Rotation
A motion by which an extended body changes orientation, with respect to other bodies in
space, without changing position (e.g., the motion of a spinning top). Rotation is of two
types.
a) Rotation about a line: Rotation of a rigid body about a line l is such a motion in
which the position in space of every point of the body that lies on l in unchanged.
The line l is called the axis of rotation.
b) Rotation about a point: Rotation of a rigid body about some point P of itself is
such a motion in which the position in space of the point P is unchanged.
1.4.3 Oscillation
A motion, which continually repeats in time with a fixed period (e.g., the motion of a
pendulum in a grandfather clock).
above mentioned types of motion are not entirely distinct: e.g., circular motion contains
elements of both rotational and oscillatory motion.
r2
G = g (1.6.3)
M
substituting (1.6.3) into (1.5.2)
F~ = m~a
comparing (1.6.4) with (1.6.1), the acceleration at the surface of the Earth is
~a = ~g
g = 9.81m/s2
g = 32.2f t/s2
Chapter 2
Rectilinear Motion
The position vector of a point P with reference to some point is discussed in chapter 1. It
can be represented as r(t). Usually rectilinear motion is considered along x − axis, so the
position vector can be represented as x(t).
In rectilinear motion, if position function of a particle is given, then differentiation is used to
define the notions of instantaneous velocity and acceleration and if acceleration function is
given, then integration is used to define the notions of instantaneous velocity and position.
If velocity function is given, then integration is used to define the notions of position and
differentiation is used to define the notions of instantaneous acceleration.
First we define the relations of instantaneous velocity and acceleration for a particle in
rectilinear motion in the form of derivatives.
9
10 2 Rectilinear Motion
2.1.2 Displacement
If the coordinate of a particle at time t1 is x(t1 ) and x1 be its position vector with respect
to some reference point usually is the origin , next its coordinate at a later time t2 is x(t2 )
and the position vector be x2 , then the displacement describes the change in position of the
particle and is given as
∆x = x2 − x1
gives the distance travelled by the particle in time interval (t1 , t2 ) Units commonly used to
express displacement and distance are m and f t.
d
v(t) = x(t) = ẋ (2.1.3)
dt
The velocity of a moving particle always points in the direction of motion.
Speed: The magnitude of the velocity, |v| is known as the speed. Units commonly used to
2.1 Rectilinear Motion 11
∆v
a(t) = lim
∆t→0 ∆t
dv
a = (2.1.4)
dt
Also we can write
dv dx
a =
dx dt
dv
a(v) = v (2.1.5)
dx
Equation (2.1.5) represents acceleration as a function of velocity. From equation (2.1.5) we
can write
adx = vdv (2.1.6)
Equation (2.1.6) is known as differential relation.
Using equation (2.1.3) in equation (2.1.4) we can express the acceleration function in terms
of the position function as
d2
d dx
a(t) = = 2 x(t)
dt dt dt
= ẍ (2.1.7)
The acceleration vector a can point anywhere. Units commonly used to express acceleration
are m/s2 or f t/s2 .
12 2 Rectilinear Motion
dx
v =
dt
represents the velocity of the particle at (x, t). If the velocity is constant, the slope of the
time-displacement curve remains the same. That is, the curve is a straight line.
2.3 Graphical Methods 13
v = v(t)
Here t is independent variable and v is dependent variable. For time-velocity curve, we take
t along horizontal axis and v along vertical axis. The plane so formed is known as tv−plane.
The slope of the tangent to the curve at (x, t)
dv
a =
dt
represents the acceleration of the particle at (x, t). If the acceleration is constant, the slope
of the time-velocity curve remains the same. That is, the curve is a straight line. Here we
have an important result:
Remark 2.3.1. The slope time-velocity curve of a particle moving in a straight line gives its
acceleration and area under the curve gives the distance travelled by the particle.
Example 2.3.1. A particle starts to move from origin O. It moves with constant accelera-
tion a for time interval t1 . At t1 it acquires a velocity v, then it starts to move with uniform
velocity v and continued for time interval t2 . After that its velocity starts to decrease with
retardation 2a and comes to rest after time interval t3 . If t is the time taken by the particle
from rest to rest, find the distance travelled by the particle in terms of v, t and a.
Solution The time velocity graph of this motion is shown in Fig. 2.3. According to
remark (2.3.1), area under the velocity curve gives the distance travelled by the particle.
14 2 Rectilinear Motion
The motion of a particle has three phases. In first phase it moves with velocity v and
acceleration a for time interval t1 , forming a triangle AOM . Its area will give the distance
travelled by the particle along line OA. According to Fig. 2.3
v2
A1 = (2.3.1)
2a
In second phase it moves with constant velocity v for time interval t2 , forming a rectangle
ABM N . Its area will give the distance travelled by the particle along line AB. According
to Fig.
|M N | = t2 and |AM | = v
A2 = |M N ||AM |
= vt2
In third phase it moves with velocity −v and retardation 2a for time interval t3 , forming
a triangle BN C. Its area will give the distance travelled by the particle along line BC.
According to Fig.
|N C| = t3 and |BN | = v
Remark 2.3.2. The slope time-acceleration curve of a particle moving in a straight line gives
its jerk and area under the curve gives the speed of the particle for defined interval.
Example 2.3.2. Let x(t) = t3 − 6t2 be the position function of a particle moving along
x−axis. Find the velocity, speed and acceleration functions, and show the graphs of position,
velocity, speed and acceleration versus time. From these graphs, determine when the particle
x(t) = t3 − 6t2
Solution Using (2.1.1), the displacement of the particle for time interval 0 ≤ t ≤ 6 is
= 0m
Using remark 2.3.1, the distance may be calculated by finding the area under the curve
6
Z
3t2 − 12t dt
|x| =
0
4 6
Z Z
3t2 − 12t dt + 3t2 − 12t dt
= −
0 4
4 6
= − t3 − 6t2 0 + t3 − 6t2 4
= 36 m/s
The speed is
|v(t)| = |3t2 − 12t|
2.4 Kinematics by Using Integration 19
Using remark 2.3.2, the distance may be calculated by finding the area under the curve
6
Z
|v| = (6t − 12) dt
0
2 4
Z Z
= − (6t − 12) dt + (6t − 12) dt
0 2
2 6
= − 3t2 − 12t0 + 3t2 − 12t2
at any time t0 , then we have sufficient information to find the constant of integration and
determine a unique velocity function. If acceleration function is given, the following two
cases arise:
a) Acceleration function is constant.
consider (2.1.4)
d
a(t) = a = v(t)
dt
dv = adt
(2.4.1) gives the velocity of the particle moving with constant acceleration at any time t.
ẍ = a(t)
d
(ẋ) = a(t)
dt
integrating with respect to t
Z
v = ẋ = a(t)dt + C1 (2.4.2)
dx
v= = v0 + at
dt
dx = (v0 + at)dt
1 2
x = at + v0 t (2.5.2)
2
If the velocity function of particle is constant over the time interval [0, t] and an initial
condition x(0) = x0 , then (2.1.3) can be written as
dx
= v
dt
It is first order separable differential equation, separating variables
dx = vdt
2.5 Displacement by Using Integration 23
Integrating
Zx Zt
dx = v dt
x0 0
x − x0 = vt
or
x = x0 + vt (2.5.3)
(2.5.3) gives the position of the particle moving with constant velocity at any time t.
b) Velocity function is variable.
If the velocity function of particle is constant over the time interval [0, t] and an initial
condition x(0) = x0 , then (2.1.3) can be written as
dx
= v(t)
dt
It is first order separable differential equation, separating variables
dx = v(t)dt
Integrating
Zx Zt
dx = v(t)dt
x0 0
or
Zt
x = x0 + v(t)dt (2.5.4)
0
(2.5.4) gives the position of the particle moving with variable velocity at any time t.
If the acceleration function of particle is constant over the time interval [0, t] and boundary
condition x(0) = 0 and x(t1 ) = x1 , then (4.1.3) can be written as
ẍ = a
1 2
x = at + C1 t + C2 (2.5.7)
2
Where C1 and C2 are constants of integration. Using x(0) = 0, we have C2 = 0, then (2.5.7)
becomes
1 2
x = at + C1 t (2.5.8)
2
1 2
x1 = at + C1 t1
2 1
or
x1 a
C1 = − t1 (2.5.9)
t1 2
(2.5.10) gives the position of the particle moving with constant acceleration at any time
t from origin
O. A comparison of equations (2.5.2) and (2.5.10) shows that the term
x1 a
t1 − 2 t1 in (2.5.10) is the velocity of the particle at t = 0.
If the acceleration function of particle is time dependent, then (4.1.3) can be written as
ẍ = a(t)
Equation (2.5.14) gives the velocity of the particle moving with constant acceleration a and
initial conditions v(0) = v0 and x(0) = x0 , after it has moved a displacement x. If the
particle starts from origin, then equation (2.5.14) becomes
Equation (2.5.15) gives the velocity of the particle starts to move from origin with velocity
v0 and constant acceleration a at any time t after it has moved a displacement x. Equation
(2.5.14) can also be written as
Equation (2.5.15) gives the displacement of the particle moving with constant acceleration
a and initial conditions v(0) = v0 and x(0) = x0 , after it has gained velocity v during time
t. If the particle starts from origin, then equation (2.5.15) becomes
Equation (2.5.17) gives the displacement of the particle starts to move from origin with
velocity v0 and constant acceleration a at any time t after it has gained a velocity v.
Equation (2.5.17) also gives the acceleration of the particle that would change its velocity
from v0 to v in a given displacement x.
dv
a(v) = v
dx
is first order separable differential equation. On separating variables we get
vdv
dx =
a(v)
On integration we have
Z
vdv
x = +A (2.5.18)
a(v)
Where A is constant of integration and can be determined if we know the velocity of the
particle for some value of x
When the acceleration depends only on the velocity v, then equation (2.1.4) may be written
as
dv
a(v) =
dt
is first order separable differential equation. On separating variables we get
dv
dt =
a(v)
On integration we have
Z
dv
t = +B (2.5.19)
a(v)
Where B is constant of integration and can be determined if we know the velocity of the
particle for some value of t
2.5 Displacement by Using Integration 27
Example 2.5.1. A particle starts to move from origin with acceleration a = kv 3 along
a straight line. If the initial velocity is v0 , find the velocity and the time spend when the
Solution Here the acceleration depends only on the velocity v, then equation (2.1.5)
can be written as
dv
v = kv 3 (2.5.20)
dx
The particle starts from origin, it means at x = 0, v = v0 . This initial condition can be
written as
v(0) = v0 (2.5.21)
Where B is constant of integration and can be determined from the following initial con-
dition. Since the particle starts from origin, it means at t = 0, v = v0 . Hence the initial
condition is
v(0) = v0 (2.5.24)
Using equation (2.5.24) in equation(2.5.23), we get
1
0 = − +B
2kv02
1
or B =
2kv02
Hence the time of motion is
1 1
t = −
2kv02 2kv 2
1 1 1
= − (2.5.25)
2k v02 v 2
Equation (2.5.25) gives the time of motion in terms of velocities. It may be calculated in
terms of initial velocity and displacement by using equation (2.5.22) in equation (2.5.25)
!
1 1 (1 − kxv0 )2
t = −
2k v02 v02
1
1 − 1 − 2kxv0 + k 2 x2 v02
= 2
2kv0
x
= (2 − kxv0 ) (2.5.26)
2v0
The displacement can be calculated by using equation (2.5.18)
Z
vdv
x = +C
kv 3
Z
dv
= +C
kv 2
1
= − +C (2.5.27)
kv
Using initial condition given by equation (2.5.21), the constant of integration is
1
C =
kv0
Then equation (2.5.27) becomes
1 1 1
x = − (2.5.28)
k v0 v
Equation (2.5.28) gives the displacement moved by the particle in terms of velocities.
2.5 Displacement by Using Integration 29
dv
2.5.7 Geometric Interpretation of a = v dx
Consider a particle moves in displacement velocity plane. At any time t the particle is at
dv
P (x, v). The slope of the curve at any point is given by dx . Thus if θ is the angle which
the tangent at P of the displacement velocity curve makes with the x axis as shown in Fig.
2.9, then
dv
tan θ = (2.5.29)
dx
Let P N be the normal to the curve at P and P M be the perpendicular to the x axis as
∠M P N = θ
MN = P M tan θ (2.5.30)
PM = v (2.5.31)
dv
v = a(x)
dx
is first order separable differential equation. On separating variables we get
vdv = a(x)dx
On integration we have
v2
Z
= a(x)dx + A (2.5.33)
2
Where A is constant of integration and can be determined if we know the position of the
particle for some value of x. The velocity may be written as
s Z
v = ± 2 a(x)dx + A (2.5.34)
s Z
dx
or = ± 2 a(x)dx + A (2.5.35)
dt
Then equation (2.5.35) is first order separable differential equation. On separating variables
we get
dx
dt = ± q R
2 a(x)dx + A
On integration we have
Z
dx
t = ± q R +B (2.5.36)
2 a(x)dx + A
Where B is constant of integration and can be determined if we know the position of the
particle for some value of t
Example 2.5.2. A particle moves along a straight line with acceleration a = x3 , where x
is the displacement of the particle from a fixed point O. If at t = 0, its distance from origin
c2
is c and velocity is √
2
find the velocity and the time spend when the particle has moved a
displacement x.
2.5 Displacement by Using Integration 31
Solution Here acceleration depends only on the velocity x, then consider equation
(2.1.5)
dv
v = x3
dx
Following equation (2.5.33), its solution is
v2
Z
= x3 dx + A
2
x4
= +A (2.5.37)
4
Where A is constant of integration and can be determined from the given initial condition.
2
That is, at x = c, v = √c 2 . Then equation (2.5.37) becomes
c4 c4
= +A
4 4
⇒ A = 0
Equation (2.5.37) becomes
v2 x4
=
2 4
x2
or v = √ (2.5.38)
2
Equation (2.5.38) gives the velocity of the particle when it has moved a displacement x
from fixed point O. Then equation (2.5.38) can be written as
dx x2
= √ (2.5.39)
dt 2
Equation (2.5.39) is first order separable differential equation. On separating variables we
get
√ dx
dt = 2 2
x
On integration we have
√
2
t = − +B (2.5.40)
x
Where B is constant of integration and can be determined from given initial condition.
That is at t = 0, x = c. Then equation (2.5.40) becomes
√
2
0 = − +B
√ c
2
or B =
c
32 2 Rectilinear Motion
Example 2.5.3. A particle moves along x − axis with acceleration a(t) = 3 sin 3t; v(0) =
Example 2.5.4. A particle moves along an r − axis with velocity v(t) = t2 − 2t m/s; (see
Fig. 2.12)
(a) Find the displacement of the particle during the time interval 0 ≤ t ≤ 3.
(b) Find the distance of the particle during the time interval 0 ≤ t ≤ 3.
v(t) = t2 − 2t
Using (1.3.5), the displacement of the particle during the time interval 0 ≤ t ≤ 3 is
Zt2 Z3
t2 − 2t dt
v(t)dt =
t1 0
3
1 3 2
= t −t
3 0
1 3
= 3 − 32 = 9 − 9
3
= 0 (2.5.47)
34 2 Rectilinear Motion
Also from graph 2.17, we see that the displacement or position vector is a zero vector. (b)
For distance, see Fig. 2.12, v(t) ≤ 0 for 0 ≤ t ≤ 2 and v(t) ≥ 0 for 2 ≤ t ≤ 3 Thus, it
follows from (2.5.5) that the distance traveled is
Zt2 Z2 Z3
2
t2 − 2t dt
|v(t)|dt = − t − 2t dt +
t1 0 2
2 3
1 3 1
= − t − t2 + t3 − t2
3 0 3 2
1 3 2 1 3 2 1 3 2
= − 2 −2 + 3 −3 − 2 +2
3 3 3
4 4
= +
3 3
8
= m (2.5.48)
3
Example 2.5.5. A particle moves starts to move along a straight line from origin, so at
ẍ = t2 + sin t + et (2.5.49)
v(0) = 0 (2.5.50)
x(0) = 0 (2.5.51)
1 3
v(t) = ẋ = t − cos t + et + A (2.5.52)
3
Where A is constant of integration and can be calculated by using (2.5.50) in (2.5.52)
0 = 0−1+1+A
⇒ A = 0
1 3
v(t) = t − cos t + et (2.5.53)
3
(b) The peed at any time t is
36 2 Rectilinear Motion
1 3 t
|v(t)| = t − cos t + e (2.5.54)
3
1 4
x = t − sin t + et + B (2.5.55)
12
Where B is constant of integration and can be calculated by using (2.5.51) in (2.5.55)
0 = 0−0+1+B
⇒ B = −1
1 4
x = t − sin t + et − 1 (2.5.56)
12
(d) Distance at any time t is
1 4 t
|x| = t − sin t + e − 1
(2.5.57)
12
2.5 Displacement by Using Integration 37
Example 2.5.6. A car starts from rest and with constant acceleration achieves a velocity
of 15 m/s when it travels a distance of 200 m. Determine the acceleration of the car and
v 2 = v02 + 2ac (x − x0 )
(15)2 = 0 + 2ac (200 − 0)
ac = 0.562m/s2
v = v0 + ac t
15 = 0 + 0.562 t
t = 26.66s
Example 2.5.7. A car moves in a straight line such that for a short time its velocity is
defined by v = (3t2 + 2t) f t/sec. Determine its position and acceleration when t = 3 s.
When t = 0, x = 0.
38 2 Rectilinear Motion
When t = 3 s
a = 6(3) + 2
a = 20 f t/s2
Example 2.5.8. A particle starts to move from origin O in a straight line with uniform
v0 = 0
t = 4s
v4 = 60 miles/hour
60 × 1760 × 3
=
60 × 60
= 88 f t/s
40 2 Rectilinear Motion
(b) Let x1 represents the displacement travelled by the particle in first 3 seconds. It can be
determined by using equation (2.5.2)
1
x1 = v0 t + at2
2
1
= 0(3) + 22(3)2
2
= 99 f t
(c) Let x2 represents the displacement travelled by the particle in last 3 seconds. It can be
determined by subtracting displacement travelled in first second from displacement travelled
in 4 seconds. Both displacements will be determined by using equation (2.5.2)
1 1
x2 = 0(4) + 22(4)2 − 0(1) + 22(1)2
2 2
= 176 − 11
= 165 f t
(d) Let x3 represents the displacement travelled by the particle in last 4th second. It can
be determined by subtracting displacement travelled in first 3 second from displacement
2.5 Displacement by Using Integration 41
Corollary 2.5.1. A particle moves in a straight line with uniform acceleration a. At time
1
x = v0 + a(2n − 1) (2.5.58)
2
Proof: Let x1 and x2 be the displacement travelled by the particle in the first n and
n − 1 units of time respectively. Then it follows from equation (2.5.2) that
1
x1 = v0 n + an2
2
1
and x2 = v0 (n − 1) + a(n − 1)2
2
Then the distance travelled by the particle in the nth unit of time is
1 2 1 2
x = x1 − x2 = v0 n + an − v0 (n − 1) + a(n − 1)
2 2
1 2 1 2
= v0 n + an − v0 n − v0 + a(n − 2n + 1)
2 2
1
= v0 + a(2n − 1)
2
Example 2.5.9. In example 2.5.8 the displacement travelled in 4th seconds can be calculated
1
x4 = v0 + a(2n − 1)
2
1
= 0 + 22[2(4) − 1]
2
= 77 f t
42 2 Rectilinear Motion
Exercise
1. Find velocity and acceleration of the particle for t ≥ 0 for the following position
curves.
2. Let x(t) = 2t3 − 21t2 + 60t + 3 be the position function of a particle moving along
x − axis, Find the velocity, speed and acceleration functions. Also plot the graphs of
position, velocity, speed and acceleration versus time. From these graphs, determine
when the particle is speeding up and slowing down.
3. Let a particle moves with constant acceleration a along straight line. Obtain the
following equations of motion by graphical method.
(a) vf = vi + at
(b) x = vi t + 12 at2
(c) 2ax = vf2 − vi2
4. Let r1 = 15t2 + 10t + 20 and r2 = 5t2 + 40t, t ≥ 0, be the position functions of cars
A and B that are moving along parallel straight lanes of a highway.
5. Two particles A and B are moving in a straight line in the same direction in such a
way that A is accelerated and B is retard. At point O, A is accelerated at the rate of
2 f t/s2 having velocity 45 miles/hour, while B is retard at the rate of 8 f t/s2 having
velocity 90 miles/hour.
(a) At what time both particles have same velocity, find this velocity.
(b) At what time both particles have same displacement from O, find this displace-
ment.
6. Two particles A and B are moving in a straight line in the same direction in such a
way that A is accelerated and B is retarded. At point O, A is accelerated at the rate
of 1.2 m/s2 having velocity 60 km/hour, while B is retard at the rate of 7.1 m/s2
having velocity 120 km/hour.
2.5 Displacement by Using Integration 43
(a) At what time both particles have same velocity, find this velocity.
(b) At what time both particles have same displacement from O, find this displace-
ment.
7. A particle starts to move from rest along straight line from origin. At any time t, its
acceleration is given as following. Find velocity, speed, displacement and distance of
the particle at ant time t.
(a) a(t) = tn
(b) a(t) = a cos t + b sin t
(c) a(t) = 3t2 − 4t
1
(d) a(t) = t
(e) a(t) = 6et + 2t
8. A particle moves along a straight line with acceleration a = −n2 x, where x is the
displacement of the particle from a fixed point O. Let the particle starts from rest
from origin. Find the velocity and the time spend when the particle has moved a
displacement x. Also find this displacement as a function of time.
9. A particle moves along a straight line with acceleration a = µx, where x is the
displacement of the particle from a fixed point O. Let the particle starts from rest
at a displacement x0 from origin. Find the velocity, displacement and the time spend
when the particle has moved a displacement x.
10. A particle moves along a straight line with acceleration a = x, where x is the displace-
ment of the particle from a fixed point O. If at t = 0, its displacement from origin
and velocity are c. Find the velocity and the time spend when the particle has moved
a displacement x.
11. A particle moves along a straight line with acceleration a = x2 , where x is the dis-
placement of the q
particle from a fixed point O. If at t = 0, its distance from origin is
3
c and velocity is 2c3 . Find the velocity and the time spend when the particle has
moved a displacement x.
12. A particle starts to move from origin along a straight line with velocity vi . If it has
acceleration v 3 , where v is the velocity of the particle at any time t, find the velocity
and the time spend when the particle has moved a displacement x. Also find this
displacement as a function of time.
13. A particle starts to move from origin along a straight line with velocity vi . If it suffers
a retardation equal to the square of the velocity of the particle at any time t, find the
velocity and the time spend when the particle has moved a displacement x. Also find
this displacement as a function of time.
44 2 Rectilinear Motion
14. A particle starts to move from rest from origin O along a straight line. It moves with
uniform acceleration a till it attains a velocity v. The motion is then retarded and
the particle comes to rest after travelling a total distance x. Find the retardation and
the total time taken by the particle from rest to rest.
15. Two particles A and B are moving in a straight line and are accelerated uniformly
such that if acceleration of A is a then acceleration of B is 12 a. Both particles starts
from origin at the same time. The motion is such that when a particle attains the
maximum velocity v, the motion is retarded uniformly in a way that retardation of A
is 12 a and retardation of B is a. Then the two particles comes to rest simultaneously
at a distance x from the starting point. Find the distance between the points where
the two particles attain their maximum velocities.
16. At time t = 0, a gunner detects a plane approaching him with a velocity v, the
horizontal and vertical displacements of the plane being h and k respectively. His gun
can fire a shell vertically upwards with an initial velocity u . Find the time when he
should fire the gun and the condition on u so that he may be able to hit the plane if
it continues its flight in the same horizontal line.
Chapter 3
A motion under the force of gravity is called projectile motion and a body executing such a
motion is called projectile. This motion is in the vertically upward or downward direction,
hence also called vertical motion. Sometimes air resistance also takes place in this motion.
So we will discuss it as free projectile motion and resisted projectile motion. For example,
an object dropped from a height, thrown vertically upwards or thrown at an angle (oblique).
We can subdivide this motion in two categories namely
a) One dimensional Projectile Motion.
45
46 3 Vertical Motion Under Gravity
the center of the earth with a constat acceleration. Therefore, all of the equations we have
derived for constant acceleration apply to an object in free fall, neglecting air resistance). All
objects fall with a constant acceleration of about 9.80 m/s2 which we call g, the acceleration
due to gravity.
In this chapter, we will consider vertically free and resisted motions. Assume that both
gravity and mass remain constant and, for convenience, choose the downward direction as
the positive direction. Since it is one dimensional motion, so the reference point may not
be with respect to coordinate system. First of all consider free vertical motion.
vf = vi + gt (3.1.1)
1
x = vi t + gt2 (3.1.2)
2
2gx = vf2 − vi2 (3.1.3)
Example 3.1.1. A body is dropped (at rest) from a height of h meters. If the motion is
free fall, then at what speed will it hit the ground? Also find the time required for it.
Solution Since it is one dimensional motion, the reference axis may be z − axis only.
As the body starts from rest, so the initial data is
t0 = 0
dz
v0 = (0) = 0
dt
z0 = 0
At time t the body is at P . At P , the only force acting on the particle is the gravitational
force W = mg. Hence by Newton’s second law of motion
F = mg
d2 z
m = mg
dt2
d2 z
= g
dt2
dv
= g
dt
is first order differential equation in variable v and can be solved as (separating variables)
v(t) = gt + v0
3.1 One Dimensional Projectile Motion or Vertical Motion 47
Example 3.1.2. A body of mass 5 slugs is dropped from a height of 100 f t with zero
velocity. Assuming no air resistance, find
(b) an expression for the position of the body at any time t, and
Solution: Since it is one dimensional motion, the reference axis may be x − axis only.
As the body starts from rest, so the initial data is
t0 = 0
v0 = 0
x0 = 0
Choose the coordinate system as in Fig. (3.2). Since there is no air resistance, so the only
F = mg
ma = mg (3.1.4)
Since
dv
a =
dt
then (3.1.4) becomes
dv
= g
dt
3.1 One Dimensional Projectile Motion or Vertical Motion 49
This differential equation is linear or, in differential form, separable; its solution is
v = gt + C1 (3.1.5)
initially the body has zero velocity i.e. at t = 0, v = 0 or v(0) = 0 using in (3.1.5), we have
0 = g(0) + C1 (3.1.6)
or C1 = 0. Thus,
v = gt
v = 32t (3.1.7)
Since
dx
v = ,
dt
then (3.1.7) becomes
dx
= 32t (3.1.8)
dt
This differential equation is also both linear and separable; its solution is
x = 16t2 + C2 (3.1.9)
0 = (16)(0)2 + C2
x = 16t2 (3.1.10)
v = 32(2.5)
= 80 f t/s
Example 3.1.3. A body is dropped from a height of h meters with speed v0 . Assuming no
air resistance, find
50 3 Vertical Motion Under Gravity
(b) an expression for the position of the body at any time t, and
Solution Since it is one dimensional motion, the reference axis may be z − axis only.
As the body starts from rest, so the initial data is
t0 = 0
dz
v0 = (0) = v0
dt
z0 = 0
At time t the body is at P . At P , the only force acting on the particle is the gravitational
F = mg
d2 z
m = mg
dt2
d2 z
= g
dt2
dv
= g
dt
3.1 One Dimensional Projectile Motion or Vertical Motion 51
is first order differential equation in variable v and can be solved as (separating variables)
v(t) = gt + v0
dz
v(t) = = gt
dt
now it is first order differential equation in variable z and can be solved as (separating
variables)
1 2
z(t) = gt + z0
2
Using initial condition z(0) = 0, the above relation is
1 2
z(t) = gt
2
Let it hits the ground after time t1 then z(t1 ) = h, then we have
1 2
h = gt
2 1
or time is
s
2h
t1 =
g
v(t1 ) = gt1
p
= 2hg.
Since it is one dimensional motion, the reference axis may be z − axis only.
√
Example 3.1.4. A particle is projected vertically upward with a velocity 2gh and at
the same time, an other is dropped from a height h with zero velocity. Assuming no air
resistance, find the height where they meet each other.
Solution Since it is one dimensional motion, the reference axis may be x − axis only.
Choose the coordinate system as in Fig. (3.4). Let both the particles meet at time t at P ,
first particle has distance x from A and the second particle has distance h − x from B. Let
x(t) be the position variable for the first particle and y(t) be the position variable for the
second particle Then the first particle has coordinates at P (x, v1 , t) and the second particle
has coordinates at P (h − x, v2 , t) = P (y, v2 , t). Since there is no air resistance, so the only
force acting on the body is w = mg.
52 3 Vertical Motion Under Gravity
t0 = 0
dx p
v1 (0) = (0) = 2gh
dt
x0 = 0
By Newton’s second law of motion, the equation of motion of the first particle is:
F = −mg
a = −g
dv1
= −g (3.1.11)
dt
Integrating (3.1.11) with respect to t
or
p
2gh = C1
3.1 One Dimensional Projectile Motion or Vertical Motion 53
1 p
x(t) = − gt2 + 2gh t + C2 (3.1.14)
2
Using initial data for first particle, at t = 0, x = 0 i.e. ,
x(0) = 0
C2 = 0
t0 = 0
dy
v2 (0) = (0) = 0
dt
y0 = 0
By Newton’s second law of motion, the equation of motion of the second particle is:
F = mg
a = g
dv2
= g (3.1.16)
dt
Integrating (3.1.16) with respect to t
v2 (t) = gt + C1 (3.1.17)
v2 (0) = 0
54 3 Vertical Motion Under Gravity
0 = g(0) + C1
or
C1 = 0
v2 (t) = gt
dy
= gt (3.1.18)
dt
(3.1.18) gives the velocity of second ball at any time t. Integrating (3.1.18) with respect to
t
1 2
y(t) = gt + C2 (3.1.19)
2
Using initial data for first particle, at t = 0, y = 0 i.e. ,
y(0) = 0
C2 = 0
Example 3.1.5. A body of mass m is thrown vertically into the air with an initial velocity
v0 . If the body encounters an air resistance proportional to its velocity, find
(b) an expression for the velocity of the body at any time t, and
(c) the time at which the body reaches its maximum height.
Solution
Choose the coordinate system as in Fig. (3.5). Consider a body of mass m is thrown
vertically into the air from A. At A, the initial data is
t0 = 0
dx
v(0) = (0) = v0
dt
x0 = 0
After time t the body is at P having coordinates P (x, v, t). Then by Newton’s second law
of motion, at time t, the net force acting on a body is
dv
F = m (3.1.25)
dt
where F is the net force on the body and v is the velocity of the body. At P there are two
forces acting on the body:
56 3 Vertical Motion Under Gravity
(1) the force due to gravity given by the weight W = −mg of the body, and
(2) the force due to air resistance given by −kv, where k > 0 is a constant of proportionality.
The minus sign is required because this force opposes the velocity; that is, it acts in the
downward (see Fig. 3.5 ). The net force F on the body is, therefore,
F = −mg − kv (3.1.26)
dv
m = −mg − kv
dt
or
dv k
+ v = −g (3.1.27)
dt m
(3.1.27) is the equation of motion in the coordinate system.
(3.1.27) is a linear differential equation. Its solution will give the expression for velocity at
any time t. The integrating factor is
k
I.F = emt
d k t k
ve m = −ge m t (3.1.28)
dt
3.1 One Dimensional Projectile Motion or Vertical Motion 57
The body reaches its maximum height when v = 0. Hence to calculate time t required to
reach maximum height, we put v = 0 in (3.1.31)
mg mg − k t
0 = − + v0 + e m
k k
mg mg − k t
= v0 + e m
k k
mg k t mg
em = v0 +
k k
k
t k mg
e m = v0 +
mg k
k kv0
emt = +1
mg
and time is
m kv0
t = ln +1 (3.1.33)
k mg
Example 3.1.6. A steel ball weighing 4.9 N is dropped from a height of 100 m with no
velocity. As it falls, the ball encounters air resistance numerically equal to 0.2v (in Newton),
where v is the velocity of the ball (in m/s). Find
(b) the time required for the ball to hit the ground.
Solution Choose the coordinate system as in Fig. (3.6). The weight of the ball is
w = 4.9 N , its mass is 0.5 kg. Let the ball is dropped from A. At A, the initial data is
t0 = 0
v0 = 0
x0 = 0
After time t the body is at P having coordinates P (x, v, t). Then by Newton’s second law
dv
F = m (3.1.34)
dt
where F is the net force on the body and v is the velocity of the body. At P there are two
forces acting on the body:
(1) the force due to gravity given by the weight W = 4.9 N of the body, and
(2) the force due to air resistance is −0.2 v. The net force F on the body is,
F = mg − kv (3.1.35)
3.1 One Dimensional Projectile Motion or Vertical Motion 59
I.F = e0.4t
C1 = 24.5 (3.1.39)
C2 = 61.25 (3.1.44)
The ball hits the ground when x(t) = 100, then (3.1.45) has the form
Although (3.1.46) cannot be solved explicitly for t, we can approximate the solution by trial
and error, substituting different values of t into (3.1.45) until we locate a solution to the
degree of accuracy we need. Such approximation is illustrated in table 3.1
From table 3.1, we see that the ball hits the ground at time t = 2.5 s
Time Distance
0 0
1 44.69289718
2 82.72860095
2.1 86.25773044
2.2 89.74454666
2.3 93.19070873
2.4 96.59781073
2.5 99.96738423 ∼
= 100
Alternatively, we note that for any large value of t, the negative exponential term will be
negligible. A good approximation is obtained by setting exponential is essentially zero, then
24.5t = 38.75
t = 1.58 s (3.1.47)
Example 3.1.7. (a) A small stone of mass m is thrown vertically upwards with initial
speed V . If the air resistance at speed v is mkv 2 , where k > 0 constant, show that the stone
2
1
2
returns to its starting point with speed V 1 − k Vg
kv 2
(b) If kv 2 << g i.e. g << 1 in the above case, find the velocity with which stone strikes
the ground.
3.1 One Dimensional Projectile Motion or Vertical Motion 61
Solution Here we consider two way motion, firstly upward and next downward motion.
Upward Motion:
The upward motion of the small stone is illustrated in the Fig 3.7.
0
Let OO = H be the maximum height of attained by it. Let at time t, its position be
at P with distance x ≤ H from the ground. Let its upward velocity at P be ẋ, and the
acceleration in the direction of increasing x is
dv
a =
dt
dv dx dv
= =v
dx dt dx
d 1 2
= v (3.1.48)
dx 2
For the problem at hand, there are two forces acting on the body at P :
(1) the force due to gravity acting in downward direction given by the weight w of the body,
which equals mg, and
(2) the force due to air resistance acting in downward direction given by mkv 2 , where k > 0
is a constant of proportionality.
The minus sign is required because the forces oppose the velocity; that is, they act in the
downward, or negative, direction (see Fig. 3.7). The net force F on the body is, therefore,
the equation of motion is:
F = −mg − mkv 2
ma = −m(g + kv 2 )
using (3.1.48), we have
d 2
v + 2kv 2 = −2g
dx
d
+ 2k v 2 = −2g (3.1.49)
dx
62 3 Vertical Motion Under Gravity
(3.1.49) is first order linear nonhomogeneous differential equation and can be solved by the
method of undetermined coefficient. The homogeneous part of (3.1.49) is:
d
+ 2k v 2 = 0 (3.1.50)
dx
The characteristic equation is
D + 2k = 0
D = −2k
vp2 = C (3.1.52)
then
d 2
v = 0 (3.1.53)
dx p
using above result in (3.1.49), we have
2kC = −2g
g
C = − (3.1.54)
k
Hence the particular solution is:
g
vp2 = − (3.1.55)
k
Hence the general solution is:
g
v 2 (x) = Ae−2kx − (3.1.56)
k
Since the stone is thrown with initial speed V with distance x = 0, i.e.
v(0) = V (3.1.57)
After attaining the maximum height H, the stone rests for a while and then returns
back. Let at time t, its position be at P with distance y ≤ H from the top.
Let its downward velocity at P be ẏ, and the acceleration in the direction of increasing y is
d 1 2
a = v (3.1.63)
dy 2
64 3 Vertical Motion Under Gravity
F = mg − mkv 2
ma = m(g − kv 2 )
v(0) = 0 (3.1.66)
Exercises
1. A stone of mass 1 kg is dropped from the top of a tower with zero velocity. Assuming
no air resistance, it hits the ground after 3 seconds. Is the mass is helpful to increase
the velocity? Also find
(a) the height of the tower and
(b) the velocity with witch the stone hits the ground.
2. A body of mass 5 kg is dropped from a height of 20 m with zero velocity. Assuming
no air resistance, find
(a) an expression for the velocity of the body at any time t.
(b) an expression for the position of the body at any time t, and
(c) the time required to reach the ground.
3. A small marble ball is thrown vertically upward with a velocity. Assuming no air
resistance, it returns back and hit the ground after 8 seconds. Find
(a) the maximum height reached by the ball and
(b) the velocity with witch the ball is thrown up.
4. Assuming no air resistance, a particle is projected vertically upward with a velocity
v0 from origin. It passes through a point at a height h from origin at time t1 . After
reaching its maximum height it returns back and passes through the same point at
time t2 . Show that
2v0
t1 + t2 =
g
2h
t1 t2 =
g
5. Two particles
√ are projected
√ simultaneously in the vertically upward direction with
velocities 2gh and 2gk, (k > h), where h and k are the maximum heights attained
by the particles. After a time t, when the two particles are still in flight, another
particle is projected upward with a velocity u. Find the condition so that the third
particle may meet the first two particles during their upward flights.
6. A particle is projected vertically upward. After a time t, another particle is sent up
from the same point with the same velocity and meets the first at height h during the
downward flight of the first. Find the velocity of projection.
7. The acceleration of a particle falling freely under the gravitational pull is equal to
k
x2
, where x is the distance of the particle from the center of the earth and k is some
constant. Find the velocity of the particle if it is let fall from an altitude R, on striking
the surface of the earth if the radius of earth is r and the air offers no resistance to
motion.
3.1 One Dimensional Projectile Motion or Vertical Motion 67
8. A ball of weight 10 N is dropped from a height of 50 m with zero velocity. If the body
encounters an air resistance proportional to half its velocity, find
10. A small stone of mass m is thrown vertically upwards with initial speed V . If the air
resistance at speed v is mkv, where k > 0 constant, show that the stone returns to its
starting point with speed U given by the relation
− kg (V +U )
g − kU = (g + kV ) e
68 3 Vertical Motion Under Gravity
Chapter 4
In this chapter we will discuss the motion of a particle in 2 and 3 dimensional cartesian
coordinate system.
tem
The concepts of position vector, displacement, velocity and acceleration will be discussed
in these systems.
4.1.2 Displacement
(a) Displacement in 2 space
69
70 4 Motion in Two and Three Dimensional Cartesian Coordinate Systems
Consider a particle moves along a curve C in 2 space as shown in Fig. 4.1 . Let at time t1
it is at point A = (x1 , y1 ) whose position vector relative to O is ~r1 = x1 î + y1 ĵ. At a latter
time t2 it is at point B = (x2 , y2 ) whose position vector relative to O is ~r2 = x2 î + y2 ĵ.
Then displacement describes the change in position of the particle and is given as
Where ∆x = x2 −x1 and ∆y = y2 −y1 are displacements along x axis and y axis respectively.
Consider a particle moves along a curve C in 3 space as shown in Fig. 4.2. Let at time t1
it is at point A = (x1 , y1 , z1 ) whose position vector relative to O is ~r1 = x1 î + y1 ĵ + z1 k̂.
At a latter time t2 it is at point B = (x2 , y2 , z2 ) whose position vector relative to O is
~r2 = x2 î + y2 ĵ + z2 k̂. Then the displacement describes the change in position of the particle
and is given as
∆t = t2 − t1
displacement
average velocity =
time internal
∆~r
~vavg =
∆t
And instantaneous velocity is defined as
∆~r
~v = lim
∆t→0 ∆t
d~r(t) 0
~v (t) = = ~r (t) (4.1.1)
dt
(a) Average Velocity in 2 space is
∆xî + ∆y ĵ
~vavg =
∆t
∆x ∆y
= î + ĵ
∆t ∆t
(b) Similarly Average Velocity in 3 space is
72 4 Motion in Two and Three Dimensional Cartesian Coordinate Systems
∆x ∆y ∆z
~vavg = î + ĵ + k̂
∆t ∆t ∆t
(c) Instantaneous Velocity in 2 space is
d~r
~v =
dt
d
= xî + y ĵ
dt
dx dy
= î + ĵ
dt dt
= vx î + vy ĵ
dy
Where vx = dx
dt is the horizontal scalar component of velocity and vy = dt is the vertical
scalar component of velocity as shown in Fig. 4.3.
~v = vx î + vy ĵ + vz k̂
∆~v
~a(t) = lim
∆t→0 ∆t
d
= ~v (t) (4.1.2)
dt
Using (4.1.1), the acceleration function in terms of the position function is
d2
d dx
~a(t) = = 2 ~x(t)
dt dt dt
00
= ~x (t) (4.1.3)
∆vx î + ∆vy ĵ
~aavg =
∆t
∆vx ∆vy
= î + ĵ
∆t ∆t
(b) Similarly Average Acceleration in 3 space is
d~v
~a =
dt
d
= vx î + vy ĵ
dt
dvx dvy
= î + ĵ
dt dt
= ax î + ay ĵ
dv
Where ax = dv y
dt is the horizontal scalar component of acceleration and ay = dt is the
x
vertical scalar component of acceleration as shown in Fig. 4.4. the acceleration function in
74 4 Motion in Two and Three Dimensional Cartesian Coordinate Systems
~a = ax î + ay ĵ + az k̂
Example 4.1.1. A car is moving along a path C as shown in Fig. 4.5 Let at time t1 = 2s
the car is at point A having position vector ~r1 = (1m)î + (4m)ĵ. After a later time t2 = 5s
the car is at point B having position vector ~r2 = (6m)î + (6m)ĵ. Find its displacement and
average velocity.
Solution The given data is
t1 = 2s
t2 = 5s
~r1 = (1m)î + (4m)ĵ
~r2 = (6m)î + (6m)ĵ
Time internal to move from A to B is
∆t = t2 − t1
= 5 − 2 = 3s
4.1 Two and Three Dimensional Cartesian Coordinate System 75
∆t = t2 − t1
= 5 − 2 = 3s
Example 4.1.3. A bus is moving along a path C as shown in Fig. 4.7 Let at time t1 = 2s
the bus is at point A having position vector ~r1 = (−2m)î + (4m)ĵ + (5m)k̂. After a later
time t2 = 5s the bus is at point B having position vector ~r2 = (3m)î + (2m)ĵ + (5m)k̂. Find
t1 = 2s
t2 = 5s
~r1 = (−2m)î + (4m)ĵ + (5m)k̂
~r2 = (3m)î + (2m)ĵ + (5m)k̂
∆t = t2 − t1
= 5 − 2 = 3s
∆x = x2 − x1
= 3 − (−2) = 5m
∆y = y2 − y1
= 2 − 4 = −2m
∆z = z2 − z1
= 5 − 5 = 0m
The displacement is
∆~r = ∆xî + ∆y ĵ + ∆z k̂
= 5î − 2ĵ
78 4 Motion in Two and Three Dimensional Cartesian Coordinate Systems
Example 4.1.4. A particle moves along a path. Let ~r = î + 4t2 ĵ be its position vector.
d~r
~v =
dt
d
= î + 4t2 ĵ
dt
d d
4t2 ĵ
= (1)î +
dt dt
= 0î + 8tĵ
= (8t m/s)ĵ
d~r
~a =
dt
d
= 8tĵ
dt
= (8 m/s2 )ĵ
~r(2) = î + 4(2)2 ĵ
= (1m)î + (16m)ĵ
Velocity of a particle at t = 2s is
~v (2) = 8(2)ĵ
= (16m/s)ĵ
4.1 Two and Three Dimensional Cartesian Coordinate System 79
At time t = 2s acceleration is
~a(2) = 8ĵ
= (8m/s2 )ĵ
Example 4.1.5. A particle moves along a path. Let the components of its position vector
are
x = t
y = 4t2 − 3
z = 1
~r = tî + 4t2 − 3 ĵ + k̂
d~r
~v =
dt
d
tî + 4t2 − 3 ĵ + k̂
=
dt
d d d
4t2 − 3 ĵ + (1)k̂
= (t)î +
dt dt dt
= î + 8tĵ + 0k̂
= (1m)î + (8t m/s)ĵ
d~r
~a =
dt
d
= î + 8tĵ
dt
= 0î + 8ĵ
= (8 m/s2 )ĵ
80 4 Motion in Two and Three Dimensional Cartesian Coordinate Systems
Exercises
1. A particle is moving along a path. Let at time t1 = 2s the car is at point A having
position vector ~r1 = (2m)î + (2m)ĵ. After a later time t2 = 4s the car is at point
B having position vector ~r2 = (3m)î + (5m)ĵ. Find its displacement and average
velocity.
2. A car starts to move from point A having position vector ~r1 = (3m)î + (2m)ĵ + (4m)k̂.
After a later time t = 4s the car is at point B having position vector ~r2 = (3m)î +
(5m)ĵ. Find its displacement and average velocity.
3. A particle is moving along a path. Let at time t1 = 2s the car is at point A having
velocity ~v1 = (3m)î + (2m)ĵ. After a later time t2 = 4s the car is at point B having
velocity ~v2 = (3m)î + (5m)ĵ. Find its average acceleration.
4. A car starts to move from point A having velocity ~r1 = (3m)î + (2m)ĵ + (4m)k̂. After
a later time t = 4s the car is at point B having velocity ~v2 = (1m)î + (3m)ĵ + (3m)k̂.
Find its average acceleration.
5. A particle moves along a path. Let ~r = 2t3 î − 21t2 ĵ be its position vector.
6. A particle moves along a path. Let ~r = 2t2 − 4 î + 4t2 ĵ + et k̂ be its position vector.
7. Analyze the motion of the particle for t ≥ 0 for the following position vectors.
Angular Motion
In this chapter we will discuss the motion of a particle with rotational effects. First of all
consider some basic concepts of rotational kinematics.
s = rθ
s = 2πr
s
θ =
r
81
82 5 Angular Motion
Angular velocity is defined as the rate of change of angular coordinate with respect to t
(time) and is denoted by ω~.
If a body is at point A and after a short interval of time ∆t it reaches at point B (see Fig.
5.2, then the change in its angular displacement is ∆θ and average angular velocity is
∆θ
ω
~ avg = â
∆t
5.1 Angular Kinematics 83
where â a unit vector in the direction of angular velocity. The instantaneous angular velocity
is
∆θ
ω
~ = lim â
∆t→0 ∆t
dθ
= â
dt
= θ̇â (5.1.1)
ω = θ̇ (5.1.2)
The directions are measured by right hand rule. The angular displacement (in magnitude)
is
θ = ωt (5.1.3)
∆ω
α
~ avg = b̂
∆t
where b̂ a unit vector in the direction of angular acceleration. The instantaneous angular
acceleration is
∆ω
α
~ = lim b̂
∆t→0 ∆t
dω
= b̂
dt
d2 θ
= b̂
dt2
or
d d2 θ~
α
~ = ω
~ = 2 (5.1.4)
dt dt
The magnitude of angular acceleration is
dω
α = (5.1.5)
dt
84 5 Angular Motion
Circular Motion
A particle executes circular motion if it travels around a circle or a circular arc. The velocity
is always directed tangent to the circle in the direction of the motion and the acceleration
is always directed radially inward (see Fig. 5.3).
dv
= 0
dt
~
~r = OP
= xî + y ĵ
5.2 Motion in Polar Coordinates or
Circular Motion 85
Its magnitude is
~ |
r = |OP
Let the position vector ~r makes an angle θ with x axis, then completing right angle triangle
OAP , we have
x = r cos θ (5.2.1)
y = r sin θ (5.2.2)
If P is a point in polar plane coordinates system then its coordinates are P (r, θ). The point
O is known as pole, x axis as initial line and y axis as terminal line. The number r is
called the radial coordinate of P and the number θ the angular coordinate (or polar angle)
of P . The number r is given as
~ |
r = |OP
p
= x2 + y 2 (5.2.3)
The number θ is given as
θ = ∠AOP
y
= arctan (5.2.4)
x
Using (5.2.1) and (5.2.2), the position vector of a particle executing circular motion is
~r = r cos θî + r sin θĵ
= rr̂ (5.2.5)
where
r̂ = cos θî + sin θĵ
is a unit vector in the direction of position vector. Its detail will be in next section.
86 5 Angular Motion
5.2.3 Velocity
Since the velocity is the time rate of change of position vector
d~r d
~v = = r cos θî + r sin θĵ
dt dt
= r −θ̇ sin θî + θ̇ cos θĵ
= rθ̇ − sin θî + cos θĵ
= rω − sin θî + cos θĵ (5.2.6)
= rω θ̂ (5.2.7)
θ̂ = − sin θî + cos θĵ
is a unit vector perpendicular to position vector. Its detail will be in next section.
The magnitude of velocity is
v = rω (5.2.8)
d~r dx dy
~v = = î + ĵ
dt dt dt
= vx î + vy ĵ
π
Since ~v makes an angle θ + 2 with x axis (see Fig. 5.7), so it can be written in its
rectangular components as
π π
~v = v cos θ + î + v sin θ + ĵ
2 2
= −v sin θî + v cos θĵ
= v − sin θî + cos θĵ
= v θ̂ = rω θ̂ (5.2.10)
5.2 Motion in Polar Coordinates or
Circular Motion 87
5.2.4 Acceleration
The acceleration of a particle at P is the time derivative of its velocity
d~v d
~a = = rω − sin θî + cos θĵ
dt dt
d
= r ω − sin θî + cos θĵ + ω̇ − sin θî + cos θĵ
dt
h i
= r ω −θ̇ cos θî − θ̇ sin θĵ + ω̇ − sin θî + cos θĵ
h i
= r −ω 2 cos θî + sin θĵ + ω̇ − sin θî + cos θĵ
drω
= −rω 2 r̂ + θ̂
dt
dv
= −rω 2 r̂ + θ̂
dt
88 5 Angular Motion
Acceleration in circular motion has both the radial (ar ) and the tangential(at ) components.
When the particle is moving with uniform speed then above relation becomes
~a = −rω 2 r̂
F = ma
For circular motion, the acceleration is given by (5.2.15), then the centripetal force is
v2
Fc = m (5.2.16)
r
It acts in the direction of centripetal acceleration.
a = rα (5.2.17)
Example 5.2.1. A particle of mass 2 kg moves in a circle of radius 0.5 m with a linear
speed of 15 m/s. Find the angular speed and the force required to keep it in circular path.
m = 2 kg
r = 0.5 m
v = 15 m/s
The force required to keep it in circular path is the centripetal force and is given by (5.2.16)
v2
Fc = m
r
(15)2
= 2
0.5
= 900 N
Example 5.2.2. A particle moves in a circle of radius 1 m. If its speed uniformly increases
r = 1m
vi = 5 m/s
vf = 10 m/s
dt = 2 s
The particle is executing uniform circular motion The tangential acceleration is average
acceleration and given by
dv
a =
dt
vf − vi 10 − 5
= =
dt 2
2
= 2.5 m/s
5.3 Motion in Radial and Transverse Plane 91
a
α =
r
2.5
=
1
= 2.5 rad/s2
In some problems a particle P may not moving along a circular path but may be located
in a better way using its polar coordinates. In such cases it becomes convenient to resolve
the velocity and acceleration into components parallel and perpendicular to the tip of the
vector OP~ . These components are called the radial and transverse components. In order
to do this, attach 2 unit vectors to the tip of P . The vector r̂ is directed along OP and the
vector θ̂ is obtained by rotating r̂ about π2 radians counterclockwise. The vector r̂ defines
the radial direction, that is the direction in which P would move if r were increased and θ
kept constant. And the vector θ̂ defines the radial direction, that is the direction P would
move if θ were increased and r kept constant. Its graphical representation is given in Fig.
5.8 Notice that in this coordinate system, unlike the xy coordinate system, the unit vectors
constantly change direction, thus they have time derivatives. These two unit vectors are
92 5 Angular Motion
given by
~r = rr̂ (5.3.3)
The above relation is already given by (5.2.5). The position vector in this system is same
as in polar plane system. Also it has only radial component as shown in Fig. 5.9.
d~r
~v =
dt
d(rr̂)
=
dt
dr̂
= ṙr̂ + r (5.3.4)
dt
5.3 Motion in Radial and Transverse Plane 93
dr̂ dr̂ dθ
=
dt dθ dt
= (− sin θ, cos θ)θ̇
= θ̇(− sin θ, cos θ)
Using (5.3.2)
dr̂
= θ̇θ̂ (5.3.5)
dt
is the velocity of the particle at any time t. The radial component of velocity is
vr = ṙ
vt = rθ̇
Example 5.3.1. A particle is constrained to move along the equiangular spiral r = aebθ , so
that the radius vector moves with constant angular velocity ω. Determine the velocity and
acceleration components.
Solution : From (5.1.3), angular distance at any time t is
θ(t) = ωt
θ̇ = ω
θ̈ = 0
5.3 Motion in Radial and Transverse Plane 95
r = aebθ
r(t) = aebωt
ṙ = ωabebωt
r̈ = ω 2 ab2 ebωt
~v = ṙr̂ + rθ̇θ̂
bωt
= ωabe
D r̂ + aωebωtE θ̂
bωt bωt
= abωe , ωae
vr = abωebωt
vt = aωebωt
at = 2abω 2 ebωt
Example 5.3.2. Suppose we have an object that travels in a circle of constant radius 3
meters with a constant angular velocity of 2 radians per second. What are the expressions
Solution :
The given data is
r = 3
ω = 2
ṙ = 0
r̈ = 0
θ̇ = ω = 2 rad/s
θ = ωt = 2t rad
θ̈ = 0
r(t) = 3
The velocity would have a radial speed of 0m/s and a tangential speed of 6m/s, both of
which were constant; i.e. the object would be traveling in a circle with a constant speed of
6m/s. In Cartesian coordinates, the expression for velocity would be
so that
which is simply
a = −12r̂ + 0θ̂
The acceleration has a constant magnitude of 12m/s2 and is directed towards the origin.
(Remember that positive r̂ points away from the origin.) This is the centripetal acceleration
as
v2
ac =
r
62
= = 12 m/s2
3
There is no tangential component of the acceleration.
In Cartesian coordinates, the expression for acceleration would be:
Let P be the position of the particle at any time t. Let Q be the neighboring position along
the curve at time t + ∆t. Then the small distance from P to Q is
P
d Q = ∆s
Let ~v and ~v + ∆~v be the velocities along the tangents to the curve at P and Q respectively.
Then from the Fig.
Which is the displacement during the time ∆t. Then the average velocity is
∆~r
~vavg = (5.4.3)
∆t
and proceeding to the limit ∆t → 0, we obtain the actual velocity of the particle at any
time t
∆~r
~v = lim
∆t→0 ∆t
d~r
=
dt
Here ∆r and ∆t are very very small quantities and are differentiable functions of time t.
Again consider Eq.(5.4.3)
∆~r ∆~r ∆s
~vavg = =
∆t ∆s ∆t
Now proceeding limit as Q → P , we have
∆~r ∆~r ∆s
~v = lim = lim
Q→P ∆t Q→P ∆s ∆t
dr ∆~r ∆s
= = lim lim
dt ∆s→0 ∆s ∆t→0 ∆t
d~r ds
=
ds dt
d~r
= ṡ (5.4.4)
ds
Here d~
r
ds is a vector in the direction of P~Q. Next we show that it is a unit vector as
∆r
= |∆r|
∆s ∆s
chord P Q |P Q|
= =
arc P Q P
d Q
5.4 Motion in Tangent and Normal Plane 99
~v = ṡt̂ (5.4.6)
|~v | = v = ṡ (5.4.7)
From (5.4.28), the path or arc length moved by the particle in time interval T is
ZT
√
s = vdt (5.4.8)
0
dt̂
= κ
ds
then
dt̂
= κn̂
ds
100 5 Angular Motion
If v exceeds vm then the string will break, and the weight will cease to be subject to a
centripetal force, so it will fly off with velocity vm along the straight-line which is tangential
to its executed circular path.
There is another force, known as centrifugal force, equal in magnitude but acts in the
opposite direction of centripetal force.
Example 5.4.1. A particle of mass 0.50 kg is attached with a massless string of length one
meter. Let it is moving in a circle with speed 10 m/s. Find the tension in the string.
Solution The given data is
m = 0.50 kg
r = 1m
v = 5 m/s
v2
Fc = m
r
(10)2
T = 0.5
1
= 50 N
2π
√
Show that the distance moved by the particle in one full turn of the helix is ω a2 ω 2 + b2
The particle starts at t = 0, and rotates in a circle of radius a about z axis with angular
velocity ω. So the time required to complete one full turn of the helical path is T = 2π
ω
At any time t(o < t < T ), the position vector is
~r = ~r(t)
= a cos ωtî + a sin ωtĵ + btk̂
= ha cos ωt, a sin ωt, bti (5.4.16)
the velocity is
d~r
~v =
dt
= h−aω sin ωt, aω cos ωt, bi (5.4.17)
Example 5.4.3. A particle is constrained to move along a circular helix so that its coordi-
π
nates at any time t are (a cos θ, a sin θ, aθ tan α), where a > 0, 0 < α < 2 are constants. The
speed increases linearly with time t from zero at t = 0 to V at t = T . Find the acceleration
at any time t < T , the motion taking place in the sense of θ increasing and starting from
~r = ~r(θ)
= ahcos θ, sin θ, θ tan αi (5.4.19)
the velocity is
d~r d~r dθ
~v = =
dt dθ dt
= aθ̇h− sin θ, cos θ, tan αi (5.4.20)
But
~v = v t̂ (5.4.23)
v = at (5.4.25)
The speed increases linearly with time t from zero at t = 0 to V at t = T , then acceleration
can be found as
V
a = (5.4.26)
T
Using Eq. (5.4.26), Eq. (5.4.25) becomes
V
ṡ = v = t (5.4.27)
T
Using Eq. (5.4.27). Eq. (5.4.21) becomes
V
~v = t t̂ (5.4.28)
T
104 5 Angular Motion
Which is the velocity of the particle at any time t, with t̂ is given by Eq.(5.4.24)
Next the acceleration of the particle is given by
F = −mg sin θ
ma = −mg sin θ
a = −g sin θ (5.4.42)
Here g is the acceleration due to gravity near the surface of the earth. The negative sign
on right hand side of (5.4.42) implies that θ and a always acts in opposite directions. Since
r = l is fixed, then using trigonometric relation
s = rθ = lθ
du
a = = lθ̈ (5.4.43)
dt
Using (5.4.43) in (5.4.42)
lθ̈ = −g sin θ
g
θ̈ − sin θ = 0
l
θ̈ + ω 2 sin θ = 0 (5.4.44)
r
g
with ω = is the frequency of oscillation.
l
Exercises
2. A radar fixed on the ground, tracks the circular motion of a rocket under the gravi-
tational force only. At an instant, the rocket is at a distance of r = 75 km away from
the radar with an inclination of θ = π3 with the ground. At that instant the rocket
has linear speed of ṙ = 1700 m/s and angular speed as ω = 0.8 degree/s. Find
3. A particle is constrained to move along a circle helix so that its coordinates at any
time t are (a cos θ, a sin θ), where a > 0, is constants. The speed increases linearly
with time t from zero at t = 0 to V at t = T . Find the acceleration at any time
t < T , the motion taking place in the sense of θ increasing and starting from the
point (a, 0).
4. A particle is constrained to move along a circular helix so that its coordinates at any
time t are (a cos θ, a sin θ, bθ), where a > 0, b 6= 0 are constants. The speed increases
linearly with time t from zero at t = 0 to V at t = T . Find the acceleration at any
time t < T , the motion taking place in the sense of θ increasing and starting from
the point (a, 0, 0).
108 5 Angular Motion
Chapter 6
A periodic motion in which the displacement is symmetrical about a point is called harmonic
motion. It may be discussed as simple harmonic motion, damped harmonic motion, forced
harmonic motion and forced and damping harmonic motion. First consider simple harmonic
motion.
Fr = −kx (6.1.1)
where k is a constant of proportionality called the spring constant. The spring is essentially
characterized by the number k. This force is called restoring force as it is always directed
toward the equilibrium position, opposite the displacement. That is, when the block is
109
110 6 Simple Harmonic Motion
displaced to the right of O in Figure 6.2, then the displacement is positive and the restoring
force is directed to the left. When the block is displaced to the left of O, then the dis-
placement is negative and the restoring force is directed to the right. This is the only force
acting on the block. Then by Newtons second law of motion, its equation of motion is
F = Fr
ma = −kx
k
a = − x (6.1.2)
m
Take
r
k
ω = (6.1.3)
m
6.1 Simple Harmonic Motion 111
Fr = −ks
And the mass attains an equilibrium position at which its weight W = mg is balanced by
the restoring force ks (see Fig. 6.3 (b)).
mg = ks
mg − ks = 0 (6.1.6)
112 6 Simple Harmonic Motion
If the mass is pulled down a distance x from its equilibrium position, the new restoring
force of the spring is
Fr = −k(x + s)
Assuming that there are no retarding forces acting on the system, then by Newtons second
law of motion, its equation of motion is
F = Fr + W
ma = −kx − ks + mg (6.1.7)
k
a = − x
m
The same expression for acceleration as given in (6.1.2) and consequently we have the same
equation of motion. Hence we have same equation of motion for both horizontal and vertical
oscillatory motions.
Motion
: Let the block of mass m is pulled a distance x = x1 from its equilibrium position (from
O at x = 0) and is released from rest at time t = 0. Then its initial data (at P ) is
t = t0 = 0
v = v0 = 0
x = x0 = x1
From (6.1.4), we can say that the motion is taking place in such a way that when the
particle is moving towards O (equilibrium position), the acceleration is acting along it so
that as the time progresses, the velocity becomes higher and higher and when the particle
is moving away from O, the acceleration is acting against it so that as the time progresses,
the velocity becomes lesser and lesser. Let at time t = t1 , the particle is at O with velocity
v = v1 . Then the final data (at O) is
t = t1
v = v1
x = x1
6.2 Expression of velocity and Position in Simple Harmonic Motion 113
Integrating
x
sin−1 = ωt + B (6.2.7)
x1
The constant of integration B can be evaluated by using initial data
x1
sin−1 = ω(0) + B
x1
B = sin−1 (1)
π
= (6.2.8)
2
114 6 Simple Harmonic Motion
Since
cos (ωt) ≤ 1
it follows that the numerical value of x cannot be greater than x1 . Thus the particle is
bound to stay within a distance x1 from the fixed point O. The length x1 is called the
amplitude of the motion. The point O is referred to as the centre of motion. Following
equation (6.2.5), the velocity of the particle at amplitude of the motion is zero and at the
centre of motion
Equation (6.2.10) gives the maximum velocity of the particle executing simple harmonic
motion. If we consider the motion starts from O, then the initial data (at O) is
Figure 6.4: Plot of displacement vs. time for Simple harmonic motion.
t0 = 0
v0 = v
x0 = 0
6.2 Expression of velocity and Position in Simple Harmonic Motion 115
0
sin−1 = ω(0) + B
x1
B = sin−1 (0)
= 0 (6.2.11)
x
sin−1 = ωt
x1
x
= sin (ωt)
x1
x = x1 sin (ωt) (6.2.12)
(6.2.9) gives the displacement of the object from equilibrium position O, at any time t,
Figure 6.5: Plot of displacement vs. time for Simple harmonic motion.
when we consider its motion starts from P and (6.2.12) gives the displacement of the object
from equilibrium position P , at any time t, when we consider its motion starts from O.
From (6.2.9) and (6.2.12), the displacement has different expression due to different bound-
ary conditions. Similarly other boundary conditions give still different forms for displace-
ment.
For a general solution of (6.1.4), the acceleration can be written as
d2 x
a = (6.2.13)
dt2
116 6 Simple Harmonic Motion
d2 x
= −ω 2 x
dt2
d2 x
+ ω2x = 0 (6.2.14)
dt2
(6.2.14) is second order homogeneous differential equation. Its characteristic equation is
m2 + ω 2 = 0 (6.2.15)
m = ± iω
If we suppose the constants as C1 = A cos φ and C2 = −A sin φ, then (6.2.16) can be written
as
Similarly
(6.2.17) and (6.2.18) are the general expressions of displacement for simple harmonic motion.
From (6.2.20), (6.2.21) and (6.2.22), we observe that the displacement, velocity and accel-
eration of the object executing simple harmonic motion are the same after an addition of
2π 2π
ω in time t. Therefore simple harmonic motion is periodic of period ω . This time period
is denoted by T .
2π
T = (6.2.23)
ω
Equation (6.2.23) indicates that time period depends only on frequency of oscillation and
is independent of amplitude.
vmax = x1 ω (6.2.24)
The maximum speed occurs at equilibrium position of the oscillation (at x = 0).
From (6.2.15), we can find the magnitude of maximum acceleration as
2π
amax = maximum of −x1 ω 2 cos ω t +
ω
118 6 Simple Harmonic Motion
amax = x1 ω 2 (6.2.25)
The magnitude of acceleration is greatest at the ends of the oscillation (at x = ±x1 ). The
Figure 6.6: Displacement velocity and acceleration of Simple harmonic motion at mean and end
points.
ẍ + ω 2 (x − b) = 0
Before going next, we present some definitions, considering the above spring mass example.
Oscillation or Vibration: The distance covered by the particle in time period T is called
oscillation
In spring mass example, the particle is at rest at P having distance x = x1 from O. Next
it is released from P to move towards O, at time t = 0. As the time progresses, it moves
towards O and its velocity increases to its maximum value v = x1 ω at O. The particle then
moves to the left of O and its velocity decreases and becomes zero at x = −x1 . Next the
particle retrace its motion backwards and comes to rest at x = x1 , i.e. it reaches at P and
completes one round known as one oscillation or vibration. The motion is then repeated.
6.3 Simple Harmonic Motion with Centre other than Origin 119
Frequency: The number of oscillation per unit time is called frequency. It is denoted by
f . Mathematically can be written as
1
f = (6.3.1)
T
The frequency of spring mass system is
ω
f = (6.3.2)
2π r
1 k
or =
2π m
Equation (6.3.2) indicates that the angular frequency ω and the frequency f are closely
related with a factor 2π. We normally express ω in rad/s and f in cycle/s or Hz (Hertz).
However, the real dimensions of both are 1/s in SI system.
Amplitude: The maximum distance covered by the particle on either side of equilibrium
position (i.e. from O) is called amplitude of simple harmonic motion. In above motion,
sin ωt ≤ 1 so x = x1 is the amplitude. It is also equal to one-half of the total range of
motion 12 (xmax − xmin ). In spring mass example, xmax = x1 and xmin = −x1 , and the
amplitude is
1
A = (xmax − xmin )
2
1
= (x1 + x1 )
2
= x1 (6.3.3)
Phase of the Motion: The time varying quantity (ωt + φ) is called the phase of the
motion.
Phase Constant or Phase Angle: The constant φ is called the phase constant or phase
angle. Its value depends on the displacement and velocity of the particle at t = 0. The
general expression for displacement is
x(t) = A cos(ωt + φ)
at t = 0, the displacement is
at t = 0, the displacement is
v(0)
= −ω tan φ
x(0)
or
−1 v(0)
φ = tan − (6.3.6)
ωx(0)
F = −kx
6.4 Kinetic and Potential Energies of Spring Mass System 121
Since the force is conservative, there exist a potential function U such that
dU
F = −
dx
Or we can write
Zb Zb
U = F = kx
a a
Let the limits of integration are from 0 → x, then the potential energy of the system is
Zx
U = kx
0
1 2
= kx (6.4.3)
2
The graph of the potential energy of spring mass system is illustrated in Fig. 6.8. The
potential energy is maximum at its endpoints (x = ±x1 ) and is zero (minimum) at its mean
position (at x = 0)
1 dx2 1 dẋ2
− k = m
2 dt 2 dt
d 1 2 d 1 2
− kx = mẋ
dt 2 dt 2
d 1 2 1 2
mẋ + kx = 0
dt 2 2
1 1
mẋ2 + kx2 = C (constant)
2 2
E =K +U = C (6.4.6)
Hence the total energy of spring mass system is constant as shown in Fig. 6.9
Example 6.4.1. A block of mass 680 g is fastened to a spring whose spring constant is
e) What is the maximum speed vmax of the oscillating block, and where is the block when
f ) What is the magnitude of maximum acceleration amax of the motion, and where is the
m = 680 g = 0.68 kg
k = 65 N/m
x(0) = 11 cm = 0.11 m
v(0) = ẋ(0) = 0m/s
r
k
ω =
m
r
65
=
0.68
= 9.87 rad/s
2π
T =
ω
2π
=
9.87
= 0.636594 ∼
= 0.64 s
1
f =
T
1
= = 1.57086 ∼
= 1.6 Hz
0.636594
d) Using (6.3.3), the amplitude is:
A = x1
= 0.11 m
vmax = x1 ω
= (0.11)(9.87) = 1.0857
∼
= 1.1 m/s
amax = x1 ω 2
= (0.11)(9.87)2 = 10.715859
∼
= 10.7 m/s2
x(0) = 11 cm = 0.11 m
v(0) = 0 m/s
h) Using (6.2.17), the displacement function x(t) for the moving block is:
x(t) = A cos(ωt + φ)
= 0.11 cos(9.87t + 0)
= 0.11 cos(9.87t) m
i) Using (6.4.2), the kinetic energy function for the moving block is:
1
k x21 − x2
K =
2
1
(65) (0.11)2 − x2
=
2
= 0.3933 − 32.5x2
126 6 Simple Harmonic Motion
j) Using (6.4.3), the potential energy function for the moving block is:
1 2
U = kx
2
1
= (65)x2
2
= 32.5x2
k) The extreme values of kinetic and potential energies exist at mean and end positions.
The kinetic energy is minimum at end points (at x = ±0.11) and maximum at mean
position (at x = 0). The values are as under:
Kmin = 0
Kmax = 0.3933
The potential energy is maximum at end points (at x = ±0.11) and minimum at mean
position (at x = 0). The values are as under:
Umax = 0.3933
Umin = 0
l) The total energy of the system is
E = K +U
= (0.3933 − 32.5x2 ) + 32.5x2
= 0.3933J
The energy of the system is constant.
6.5 Relation Between Uniform Circular Motion and Simple Harmonic Motion 127
Harmonic Motion
There is a correspondence between simple harmonic motion and uniform circular motion.
Consider a particle P moves at constant speed v around a circle of radius r. At the instant
~ = ~r makes an angle θ with x axis.
t, P has coordinates (x, y). Then the radius vector OP
Also the x coordinate in polar form is
x = r cos θ
The velocity component along x axis is
dx dθ
vx = = −r sin θ
dt dt
= −r sin θ ω
The acceleration component along x axis is
dvx dθ
ax = = −r cos θ ω
dt dt
= −xω 2
Hence the system executes simple harmonic motion.
Example 6.5.1. The wheel of a car has a radius of 0.30 m and it being rotated at 15
revolutions per second on a tire-balancing machine. Determine the angular speed and the
r = 0.30 m
f = 15 Hz
And the speed at which the outer edge of the wheel moving is
2πr 2π(0.3)
v = =
T 0.0667
= 28.274 m/s
θ̈ = −ω 2 sin θ
with
r
g
ω = (6.5.2)
l
is the frequency of oscillation. The angle θ is defined with respect to the equilibrium
position. When θ > 0 , the bob has moved to the right, and when θ < 0 , the bob has
moved to the left. The object will move in a circular arc centered at the pivot point in
the absence of any dissipation due to air resistance or frictional forces acting at the pivot.
When there is small oscillation, i.e. θ is very very small, we can write sin θ ∼
= θ, then above
equation can be written as
θ̈ = −ω 2 θ
2π 2π
ω = =
T 2.5
= 2.5133 rad/s
b) The length of the pendulum is
Since the angular frequency is
r
g
ω =
l
or the length is
g 9.8
l = 2
=
ω (2.5133)2
= 1.5515 m
130 6 Simple Harmonic Motion
If we take this pendulum on the Moon, it will oscillate under Moon’s acceleration of gravity
with same length, so its angular frequency is
r s
g 1.67
ω = =
l (1.5515)
= 1.0375 rad/s
2π 2π
T = =
ω 1.0375
= 6.056 ∼
= 6.06 s
F = −bẋ
Where b > 0 is damping constant and negative sign indicates that the damping force acts
in a direction opposite to the direction of motion. Then by Newton’s second law of motion,
the equation of motion is
F = −kx − bẋ
mẍ + bẋ + kx = 0 (6.5.3)
(6.5.3) is second order homogenous linear differential equation with constant coefficients.
Its linear standard form is
b k
ẍ + ẋ + x = 0
m m
ẍ + β ẋ + ω 2 x = 0 (6.5.4)
b
with β = m is damping constant. Let
x = emt (6.5.5)
6.5 Relation Between Uniform Circular Motion and Simple Harmonic Motion 131
d2 + βd + ω 2 = 0 (6.5.6)
β
d=−
2
and the solution is
Thus x → 0 as t → ∞. Hence the motion decays with time as shown in Fig. 6.12.
1
d = (−β ± iγ)
2
with −γ 2 = β 2 − 4ω 2 and the solution is
A = Be−(β/2)t (6.5.12)
d = −a, −b
with
p
−β + β 2 − 4ω 2
−a = d1 =
p2
−β − β 2 − 4ω 2
−b = d2 =
2
6.5 Relation Between Uniform Circular Motion and Simple Harmonic Motion 133
Thus x → 0 as t → ∞. Hence the motion decays with time as shown in Fig. 6.14. The
three damping cases, collectively are shown in Fig. 6.15.
Example 6.5.4. A block of mass 500 g is fastened to a spring whose spring constant is
2 N/m. The block is pulled a distance one meter from its equilibrium position on a surface
and is released with velocity 1 m/s. Assuming 5 is the damping constant of the surface.
Find
g) the amplitude.
Solution In this problem, we have
m = 500 g = 0.5 kg
k = 2 N/m
β = 5
And the initial data is (at t = 0 )
x(0) = 1 m
v(0) = ẋ(0) = 1 m/s
a) the angular frequency.
First of all we will find its angular frequency, for this use (6.1.3)
r
k
ω =
m
r
2
=
0.5
= 2 rad/s
6.5 Relation Between Uniform Circular Motion and Simple Harmonic Motion 135
Fr = kx = 2x
b = βm
= 5(0.5) = 2.5
Fd = 2.5v
β 2 − 4ω 2 = 25 − 4(4)
= 9>0
ẍ + β ẋ + ω 2 x = 0
ẍ + 5ẋ + 4x = 0 (6.5.16)
136 6 Simple Harmonic Motion
d2 + βd + ω 2 = 0
d2 + 5d + 4 = 0 (6.5.17)
d = −1 and − 4 (6.5.18)
C1 + C2 = 0 (6.5.20)
or
8
e3t =
5
1 8
t = ln
3 5
= 0.157s (6.5.25)
g) The amplitude
ẍ(0.157) < 0
5 −t 32 −4t
ẍ(t) = e − e
3 3
5 −0.157 32 −4(0.157)
ẍ(0.157) = e − e
3 3
= −4.267 < 0
5 −0.157 2 −4(0.157)
x(0.157) = e − e
3 3
xmax = 1.069 m (6.5.26)
138 6 Simple Harmonic Motion
Example 6.5.5. A block of mass 250 g is fastened to a spring whose spring constant is
4 N/m and can move on a surface. The block is released from its equilibrium position to
move towards left with a velocity 3 m/s . Assuming the damping force numerically equals
to two times the instantaneous velocity acts on the system. Determine the extremum of its
motion.
Solution In this problem, we have
m = 250 g = 0.25 kg
k = 4 N/m
Fd = 2v = 2ẋ
And the initial data is (at t = 0 )
x(0) = 0 m
v(0) = ẋ(0) = −3 m/s
The restoring force is
Fr = kx = 4x
Then by Newton’s second law of motion, the equation of motion is
F = −Fd − Fr
mẍ = −2ẋ − 4x
0.25ẍ = −2ẋ − 4x
ẍ + 8ẋ + 16x = 0 (6.5.27)
(6.5.27) is the equation of motion. From it, we have
β = 8
ω = 4
The nature of the damping is determined as
β 2 − 4ω 2 = 64 − 4(16)
= 0
The system is under critical damping.
Considering (6.5.27), the characteristic equation is
d2 + βd + ω 2 = 0
d2 + 8d + 16 = 0
(d + 4)2 = 0 (6.5.28)
6.5 Relation Between Uniform Circular Motion and Simple Harmonic Motion 139
d = −4 and − 4
A1 = 0
A2 = −3 (6.5.32)
−3(1 − 4t)e−4t = 0
or
1
t = = 0.25s (6.5.35)
4
The position function has an extremum at t = 0.25. For maxima, it has to follow
ẍ(0.25) < 0
140 6 Simple Harmonic Motion
ẍ(0.25) > 0
x(0.25) = −3te−4(0.25)
xmin = −0.276 m
The particle is towards left, at a distance 0.276m, from equilibrium position, and the am-
plitude of the motion is
xmax = 0.276 m
Example 6.5.6. A block of mass 500 g is fastened to a spring whose spring constant is
5 N/m and can move on a surface. The block is pulled a distance 2 cm towards left from
its equilibrium position on a surface and released from rest. Assuming the damping force
6.5 Relation Between Uniform Circular Motion and Simple Harmonic Motion 141
numerically equals to the instantaneous velocity acts on the system. Determine its path of
motion.
m = 500 g = 0.5 kg
k = 5 N/m
Fd = v = ẋ
x(0) = −2 m
v(0) = ẋ(0) = 0m/s
Fr = kx = 5x
F = −Fd − Fr
mẍ = −ẋ − 5x
0.5ẍ = −ẋ − 5x
ẍ + 2ẋ + 10x = 0 (6.5.36)
β = 2
√
ω = 10 rad/s
β 2 − 4ω 2 = 4 − 4(10)
= −36 < 0
d2 + βd + ω 2 = 0
d2 + 2d + 10 = 0
d = −1 ± 3i
142 6 Simple Harmonic Motion
Consider spring mass system with f (t), a driving force causing an oscillatory motion of
the support of the spring. Let there is no damping force, then by Newton’s second law of
motion, its equation of motion is
F = −Fr + f (t)
mẍ + kx = f (t)
k f (t)
ẍ + x =
m m
2
ẍ + ω x = F (t) (6.5.41)
6.5 Relation Between Uniform Circular Motion and Simple Harmonic Motion 143
k
With ω 2 = m and F (t) = fm(t)
with frequency γ 6= ω.
Initially the system is resting at equilibrium position. Next an external force F (t) = F0 sin γt
is applied and it executes oscillatory motion. Then (6.5.41) takes the form
ẍ + ω 2 x = F0 sin γt (6.5.42)
x(0) = 0
ẋ(0) = 0
(6.5.42) is second order linear nonhomogeneous differential equation and has solution
A1 ω 2 − γ 2 = 0
(6.5.47)
A2 ω 2 − γ 2 = F0
(6.5.48)
A1 = 0
C1 = 0
Example 6.5.7. A 128 lb weight is attached to a spring having a spring constant of 64lb/f t.
The weight is pulled a distance 6 inches towards left, from its equilibrium position. From
there, it is started in motion with no initial velocity by applying an external force F (t) =
8sin4t to the weight. Assuming no air resistance, find the subsequent motion of the weight.
Solution In this problem, we have
W = 128 lb
k = 64lb/f t
F (t) = 8 sin 4t
x(0) = −6 in = −0.5 f t
v(0) = ẋ(0) = 0f t/s
Then
8
F (t) = sin 4t = 2 sin 4t
4
And the angular frequency frequency is calculated by using (6.1.3)
r
k
ω =
m
64 √
r
= = 16
4
= 4 rad/s
ẍ + ω 2 x = F0 sin γt
ẍ + 16x = 2 sin 4t (6.5.53)
146 6 Simple Harmonic Motion
(6.5.53) is second order linear nonhomogeneous differential equation and has solution
Using initial conditions, a particular solution is calculated. Using initial condition x(0) =
−0.5, (6.5.60) implies that
C1 = −0.5
1
ẋ(t) = 0.5 sin 4t + 4C2 cos 4t + t sin 4t − cos 4t (6.5.61)
4
1
C2 =
16
Using these coefficients, a particular solution is
1 1 1
x(t) = − cos 4t + sin 4t − t cos 4t (6.5.62)
2 16 4
Its graphical representation is given in Fig. 6.19.
During damped oscillatory motion, the oscillations eventually die away due to frictional
energy losses. To maintain the motion of a damped oscillator, an external force f (t) is
148 6 Simple Harmonic Motion
applied to the system. Let the damping force is proportional to instantaneous velocity,
then by Newton’s second law of motion, its equation of motion is
F = −Fr − Fd + f (t)
mẍ + bẋ + kx = f (t)
b k f (t)
ẍ + ẋ + x =
m m m
ẍ + β ẋ + ω 2 x = F (t)
Example 6.5.8. A body 200 g mass is attached to a spring having a spring constant of
2N/m. The weight is pulled a distance 50 cm towards right, from its equilibrium posi-
tion. From there, it is started in motion with no initial velocity by applying an external
force f (t) = 5cos4t to the mass. Assuming there is 1.2 v m/s a resistance force, find the
1
m = 200 g = kg
5
k = 2N/m
f (t) = 5 cos 4t
x(0) = 50 cm = 0.5 m
v(0) = ẋ(0) = 0m/s
Fr = kx = 2x m
Fr = k ẋ = 1.2ẋ m/s
f (t) = 5 cos 4t
F = −Fr − Fd + f (t)
mẍ + bẋ + kx = f (t)
1
ẍ + 1.2ẋ + 2x = 5 cos 4t
5
ẍ + 6ẋ + 10x = 25 cos 4t (6.5.63)
6.5 Relation Between Uniform Circular Motion and Simple Harmonic Motion 149
(6.5.63) is second order linear nonhomogeneous differential equation and has solution
d2 + 6d + 10 = 0 (6.5.65)
d = −3 ± i
50
A2 =
51
150 6 Simple Harmonic Motion
Exercises
1. A block of weight 490 N is fastened to a spring whose spring constant is 1.2 N/m.
The block is pulled a distance one meter from its equilibrium position on a surface and
is released with velocity 3 m/s. Assuming 2 is the damping constant of the surface.
Find
a) the angular frequency.
b) the restoring force.
c) the damping force function.
d) the nature of the damping.
e) the path of motion.
f) the time for which the function has a maximum.
g) the amplitude.
2. A particle describing simple harmonic motion has speeds 6 m/s and 4 m/s when its
distances from the center are 4 m and 4.5 m respectively. Find
a) the angular frequency.
b) its acceleration at these distances.
c) the time period of motion.
d) the amplitude of motion.
3. A particle describing simple harmonic motion has velocities 5f t/sec and 4f t/sec.When
its distances from the centre are 12f t and 13f t respectively find the time period of
the motion
152 6 Simple Harmonic Motion
4. The maximum velocity that a particle executing simple harmonic motion of amplitude
a attains.is v if it is disturbed in such a way its maximum velocity becomes nv find
the change in the amplitude and the time period of the motion.
5. A point describes simple harmonic motion in such a way that its velocity and acceler-
ation at point P are u and f respectively and the corresponding quantities at another
point Q are v and g find the distance P Q.
v 2 = n2 ax2 + 2bx + c
show that P executes a simple harmonic motion. Find the centre, the amplitude and
the time period of the motion.
7. A particle describes simple harmonic motion with frequency N . If the greatest velocity
is V , find the amplitude and the maximum value of the acceleration of the particle.
Also show√ that the velocity v at a distance x from the centre of motion is given by
v = 2π a2 − x2 , where a is the amplitude.
Chapter 7
Two dimensional projectile motion is namely horizontal motion and vertical motion. These
two motions are independent of each other. In cricket a bowler bowls a ball, a batsman hit
a ball by bat and fielder throw a ball. In all activities we observe two dimensional projectile
motion. This motion is simply known as projectile motion. The ball (object) executing this
motion is called projectile. Here we will discuss this motion as:
a) Angle of projection
b) Projection velocity
153
154 7 Two Dimensional Projectile Motion
an angle α with the horizontal. The point O is named as point of projection, the velocity ~v0
is the velocity of projection and the angle α is called angle of projection. In the absence of
air resistance, the particle executes projectile motion. The initial velocity ~v0 can be written
as
~r(0) = ~0
x(0)î + y(0)ĵ = 0î + 0ĵ
implies that
x(0) = 0 (7.2.2)
and
y(0) = 0 (7.2.3)
at t = 0
At P , the only force acting on the particle is force of gravity acting in the downward
direction. Then by Newton’s second law of motion its equation of motion is
F~ ~
= W
m~a = −m~g
ẍî + ÿ ĵ = −g ĵ (7.2.7)
ẍ(t) = 0 (7.2.8)
ÿ(t) = −g (7.2.9)
ẋ(t) = A1 (7.2.10)
At t = 0, (7.2.10) becomes
ẋ(0) = A1 (7.2.11)
A1 = v0 cos α (7.2.12)
(7.2.13) gives the horizontal scalar component of velocity of the particle at any time t.
Integrating it with respect to t
At t = 0, (7.2.14) becomes
B1 = 0 (7.2.16)
(7.2.17) gives the horizontal component of position of the particle at any time t.
Next for vertical component integrate (7.2.9) with respect to t
At t = 0, (7.2.18) becomes
ẏ(0) = A2 (7.2.19)
A2 = v0 sin α (7.2.20)
(7.2.21) gives the vertical scalar component of velocity of the particle at any time t.
Integrating (7.2.21) with respect to t
1
y(t) = − gt2 + (v0 sin α)t + B2 (7.2.22)
2
At t = 0, (7.2.22) becomes
1
y(0) = − g(0)2 + (v0 sin α)(0) + B2 (7.2.23)
2
Using (7.2.3), (7.2.23) becomes
B2 = 0 (7.2.24)
158 7 Two Dimensional Projectile Motion
1
y(t) = − gt2 + (v0 sin α)t (7.2.25)
2
(7.2.25) gives the vertical component of position of the particle at any time t.
From (7.2.17), we can find the time required to reach the particle at P as
x(t)
t = (7.2.26)
v0 cos α
(7.2.27) gives the path of the projectile of the particle at any time t.
(x − h)2 = −4p(y − k)
with (h, k) is the vertex, p is the distance from the vertex to the focus and the vertex to
the directrix. The axis of the parabola is
x = h
y = k
7.2 Projectile Motion of a Projectile Projected from Origin 159
2v02 cos2 α
and the length of the latus rectum is g . The focus is a point F (x1 , y1 ) on the axis of
the parabola with
and
1
y1 = ordinate of the vertex − (length of the latus rectum)
4
2 2 2 2
v0 sin α 1 2v0 cos α
= −
2g 4 g
v 2
= − 0 cos2 α − sin2 α
2g
v2
= − 0 cos 2α
2g
Hence the focus is
v02 v2
F = sin 2α, − 0 cos 2α
2g 2g
Height of Directrix The directrix of a parabola is a line perpendicular to the axis of the
parabola and is given by
1
y = height of the vertex + (length of the latus rectum)
4
2 2 2 2
v0 sin α 1 2v0 cos α
= +
2g 4 g
2
v0
sin2 α + cos2 α
=
2g
or
v02
y2 = (7.2.30)
2g
(7.2.30) gives the height of directrix of parabola.
Time of Flight As the particle is moving under gravity, so it will strike horizontal axis (or
plane) after time t. At that time we have
y = 0
1 2
− gt + (v0 sin α)t = 0
2
1
t − gt + (v0 sin α) = 0
2
Since t 6= 0, then
1
− gt + (v0 sin α) = 0
2
or
2v0 sin α
t = tr = (7.2.31)
g
7.2 Projectile Motion of a Projectile Projected from Origin 161
2v0 sin α
x(t) = (v0 cos α)
g
v02
xR = sin 2α (7.2.32)
g
v02
xR = sin (π − 2α)
g
v02 h π i
= sin 2 −α
g 2
v02
= sin (π − 2α)
g
v02
= sin 2α (7.2.33)
g
(7.2.32) and (7.2.33) implies that the horizontal range is the same for both the angles of
projection α and π2 − α as shown in Fig. 7.5.
Maximum Horizontal Range The maximum horizontal range is given by the maximum
of right hand side of (7.2.32). This maximum is accompanied with the maximum of sin 2α
162 7 Two Dimensional Projectile Motion
v02
xR−max = (7.2.34)
g
Hight of the Projectile The vertex of the parabola is the highest point of the trajectory,
so ordinate of the vertex gives the hight attained by the projectile.
v02 sin2 α
yH = (7.2.35)
2g
π
When the projectile has maximum horizontal range, angle of projection is 4. In this case
height of the projectile is
v02 π
y3 = sin2
2g 4
2
v0 1
=
2g 2
v02
= (7.2.36)
4g
Maximum Hight of the Projectile The projectile can attain maximum hight if
sin2 α = 1
which implies that α = π2 , then the projectile will be projected vertically upward and the
motion will be one dimensional motion. Hence the height will be
v02
ymax = (7.2.37)
2g
(7.2.37) gives the maximum hight of the projectile.
Example 7.2.1. A golfer hit a ball and it is projected with a speed of 10 m/s at an incli-
π
nation 6 with the ground.
Solution
1
y(t) = − gt2 + (v0 sin α)t
2
1 π
= − 9.8 t2 + 10(sin )t
2 6
2
= (−4.9 t + 5t) m
v02
xR = sin 2α
g
(10)2 π
= sin 2
9.8 6
= 8.837 m
v02
xR−max =
g
(10)2
=
9.8
= 10.2041 m
164 7 Two Dimensional Projectile Motion
2v0 sin α
t1 =
g
π
2(10) sin 6
=
9.8
= 1.0204 s
v02 sin2 α
y1 =
2g
π
(10)2 sin2
6
=
2(9.8)
= 1.2755 m
Example 7.2.2. A built is fired from a cannon having muzzle velocity 1 mile/s.
Solution Muzzle velocity The velocity with which a bullet or shell leaves the muzzle of
a gun is called muzzle velocity.
Here the initial speed is in mile/s. It should be in mile/h or f t/s. Let it be in f t/s. As
v0 = 5280 f t/s
v02
xR−max =
g
(5280)2
=
32
= 871200 f t
7.2 Projectile Motion of a Projectile Projected from Origin 165
v02
yr =
4g
(5280)2
=
4(32)
= 217800 f t
Corollary 7.2.1. A projectile projected from origin with initial speed v0 with angle of
projection α (moving under gravity only), has horizontal range xR and maximum height
Solution: Horizontal range xR of a projectile projected from origin with initial speed v0
with angle of projection α (moving under gravity only), is given by (7.2.32)
v02
xR = sin 2α
g
and maximum height yH is given by (7.2.35)
v02 sin2 α
yH =
2g
Here the goal is to express v0 in terms of xR , yH and g. From (7.2.32) we can write
v02 v2
x2R = 4 sin2 α 0 cos2 α
g g
Using (7.2.35), we can write
v02
x2R = 8yH cos2 α
g
166 7 Two Dimensional Projectile Motion
or
x2R
v02 cos2 α = g (7.2.41)
8yH
Also from (7.2.35), we can write
v02 sin2 α = 2gyH (7.2.42)
Adding (7.2.41) and (7.2.43)
x2R
v02 cos2 α + sin2 α
= g + 2gyH (7.2.43)
8yH
x2 + 16yH 2
v02 = g R (7.2.44)
8yH
Taking square root we have
s
2
x2R + 16yH
v0 = g
8yH
Using (7.2.38) in (7.2.41) we can write
x2
cos2 α =
x2R + 16yH
2
or
xR
cos α = q
x2R + 16yH
2
i.e.
xR
α = arccos q
2 2
xR + 16yH
or
4yH
sin α = q
x2R + 16yH
2
i.e.
4yH
α = arcsin q
x2R + 16yH
2
7.2 Projectile Motion of a Projectile Projected from Origin 167
Corollary 7.2.2. A projectile projected from origin with initial speed v0 (moving under
gravity only) has same horizontal range xR for two angles of projections, namely α and
π
2 − α. If y1 is the maximum height with angle α and y2 is the maximum height with angle
π
2 − α, then show that
√
x R = 4 y1 y2
v02
xR = sin 2α
g
v02 sin2 α
y1 =
2g
π
And for angle of projection 2 − α, the horizontal range xR is same given by (7.2.32)
v02
xR = sin 2α
g
v02 π
y2 = sin2 −α
2g 2
2 2
v0 cos α
=
2g
Corollary 7.2.3. A projectile projected from origin with initial speed v0 (moving under
gravity only) has maximum horizontal range xR−max . If y is its maximum height, then
show that
xR−max
Its maximum height is y = (7.2.45)
4
√
Its initial speed is v0 = xR−max g (7.2.46)
r
2xR−max
Its time of flight is tr = (7.2.47)
g
Solution: When a projectile has maximum horizontal range its angle of projection is
π
4. In this case maximum horizontal range xR−max is given by (7.2.34)
v02
xR−max =
g
and maximum hight attained by the projectile is given by (7.2.36)
v02
y =
4g
Using (7.2.34), we have
xR−max
y =
4
From (7.2.34), initial speed of a projectile can be written as
√
v0 = xR−max g
Corollary 7.2.4. A particle of mass m is projected from origin with initial speed v0 and
angle of projection α. At any time t the particle is at A(x, y). Another particle of same
7.2 Projectile Motion of a Projectile Projected from Origin 169
mass is dropped from B, a point on the directrix of the trajectory, vertically above A as
shown in Fig. 7.6 . Neglecting air resistance in both motions, show that both particles have
same speed at A.
Solution: This can be shown by the law of conservation of energy. First consider the
particle executing two dimensional projectile motion.
At O its speed is v0 , so its kinetic energy is
1
T1 = mv 2
2 0
And its height is y = 0, so its potential energy is
U1 = mgh = 0
Total energy at O is
1
E = T1 + U1 = mv 2 + 0
2 0
1
= mv 2
2 0
Next at A its speed is v1 , so its kinetic energy is
1
T2 = mv 2
2 1
And its height is y, so its potential energy is
U2 = mgh = mgy
170 7 Two Dimensional Projectile Motion
Total energy at A is
1
E = T2 + U2 = mv 2 + mgy
2 1
By law of conservation of energy, we have
Total energy at O = Total energy at A
1 1
mv02 = mv 2 + mgy
2 2 1
v02 = v12 + 2gy (7.2.48)
Next we consider the particle executing one dimensional projectile motion.
At B its speed is vb = 0, so its kinetic energy is
1
T3 = mv 2 = 0
2 b
v02
Since B lies on the directrix of the parabolic trajectory, so its height is y = 2g , and its
potential energy is
v02
U3 = mgh = mg
2g
v02
= m
2
Total energy at B is
1
E = T3 + U3 = 0 + mv02
2
1
= mv 2
2 0
As the particle is dropped, at A it gains a speed v2 , so its kinetic energy is
1
T4 = mv 2
2 2
And its height is y, so its potential energy is
U2 = mgh = mgy
Total energy at P is
1
E = T4 + U4 = mv 2 + mgy
2 2
By law of conservation of energy, we have
Total energy at B = Total energy at A
1 1
mv02 = mv 2 + mgy
2 2 2
v02 = v22 + 2gy (7.2.49)
7.2 Projectile Motion of a Projectile Projected from Origin 171
v12 = v22
or v1 = v2
Hence proved.
Corollary 7.2.5. A particle of mass m is projected from origin with initial speed v0 and
angle of projection α. Neglecting air resistance, find the least speed v0 , so that it passes
Solution: A particle of mass m is projected from origin with initial speed v0 and angle
of projection α. At any time t the particle is at A(x, y). Let ~r be its position vector
~r = xî + y ĵ
d~r dx dy
~v = = î + ĵ
dt dt dt
= v0 cos αî + (v0 sin α − gt) ĵ
172 7 Two Dimensional Projectile Motion
Since
1
y = v0 sin αt − gt2
2
so the speed of the particle at A is
q
v = v02 − 2gy
s
v02
= 2g −y
2g
p
= 2g (height of the directrix − ordinate of A) (7.2.50)
If B is a point on the directrix of the trajectory, vertically above A as shown in Fig. 7.7,
then speed of the particle at A is
p
v = 2g |AB| (7.2.51)
If the particle passes through P with speed v1 , then by (7.2.50), its speed is
p
v1 = 2g (height of the directrix − ordinate of P)
s
v02
= 2g − y1 (7.2.52)
2g
If M is a point on the directrix of the trajectory, vertically above P as shown in Fig. 7.7,
then using (7.2.51) speed of the particle at P is
p
v1 = 2g |M P | (7.2.53)
v02
2g |M P | = 2g − y1
2g
7.3 Parabola of Safety 173
or
If the particle passes through Q with speed v2 , then by (7.2.50), its speed is
p
v2 = 2g (height of the directrix − ordinate of Q)
s
v02
= 2g − y2 (7.2.55)
2g
If N is a point on the directrix of the trajectory, vertically above Q as shown in Fig. 7.7,
then using (7.2.51) speed of the particle at Q is
p
v2 = 2g |N Q| (7.2.56)
v02
2g |N Q| = 2g − y2
2g
or
2v02 = 2g (|M P | + y1 + |N Q| + y2 )
or v02 = g (y1 + y2 + |M P | + |N Q|) (7.2.58)
If S is the focus of parabolic trajectory, then using the focus-directrix property of the
parabola, |M P | = |P S| and |N Q| = |SQ| as shown in Fig. 7.7. Then (7.2.58) can be
written as
Now v02 is least when |P S| + |SQ| is least, which is least when S lies on P Q, i.e., when
|P S| + |SQ| = |P Q|
or |M P | + |N Q| = |P Q|
Hence the least speed v0 , so that it passes through two points P and Q at heights y1 and
y2 respectively is
p
(v0 )min = g (y1 + y2 + |P Q|) (7.2.59)
174 7 Two Dimensional Projectile Motion
g sec2 α 2
y(x) = − x + x tan α
2v02
sec2 α = 1 + tan2 α
gx2 g
2 tan2 α − x tan α + 2 x2 + y = 0 (7.3.1)
2v0 2v0
7.3 Parabola of Safety 175
or we can write
2g g 2
x +y = 1
v02 2v02
v02
g 2
x + y =
2v02 2g
or
2v 2 v02
2
x = − 0 y− (7.3.2)
g 2g
(7.3.2) represents a parabola, known as parabola of safety as shown in Fig. 7.8. Its vertex
is
v2
V = 0, 0 ,
2g
v02
y1 = ,
g
2v02
the length of the latus rectum is g and origin is the focus
F = O = (0, 0)
Example 7.3.1. A shell bursts on contact with the ground and pieces from it fly in all
directions with all speeds up-to 80 f t/s. Prove that a man standing 100 f t away is in
Solution The point of contact with the ground may be taken as origin O. From O all
pieces fly in different directions with all speeds up-to 80 f t/s. We will discuss projectile
motion for one of the pieces with initial speed v0 = 80 f t/s. We will discuss projectile
motion for one of the pieces that has 100 f t horizontal range. Here we have two cases:
This horizontal range may have two angles of projections. If first angle of projection is α,
then the other angle of projection is π2 − α. These angles may be calculated using (7.2.32),
describing the horizontal range of a projectile.
v02
xR = sin 2α
g
and α is
1 gxR
α = arcsin
2 v2
0
1 32(100)
= arcsin
2 (80)2
1 1
= arcsin
2 2
1 π
=
2 6
π
= = 15◦
12
π
the other angle 2 − α is
π π π
−α = −
2 2 12
5π
= = 75◦
12
Time of flight of the projectile is given by (7.2.31)
2v0 sin α
t =
g
π
Let t1 be the time of flight with elevation α = 12
π
2(100) sin 12
t1 =
32
7.4 Projectile Motion of a Projectile Projected from a Relative Height from Origin 177
5π
t2 be the time of flight with elevation α = 12
5π
2(100) sin 12
t2 =
32
The man is danger for the time t1 − t2
π
2(100) sin 4
t3 =
32
5
= √
2
x(0) = 0 (7.4.1)
178 7 Two Dimensional Projectile Motion
and
y(0) = y0 (7.4.2)
B 3 = y0 (7.4.5)
Example 7.4.1. A flying squirrel launches itself from the top of a 10 m high tree. The
squirrel leaves the tree with a velocity of 5 m/s making an inclination 30◦ relative to a level
above the ground and approaches a shorter tree of height of 5 m located 8.78 m away from
Solution The height of the shorter tree can be set as zero level above the ground and
the horizontal line from the mid of longer tree to the top of small tree can be considered
as x axis and y axis along the longer tree. Its center can be considered as the origin. The
system is illustrated in Fig. 7.12. Here C is the top of longer tree from where the squirrel
jumps and approaches to D, executing projectile motion. The given data is
v0 = 5 m/s
xR = 8.78 m
Y0 = 10 − 5 = 5 m
α = 30◦
180 7 Two Dimensional Projectile Motion
vx = v0 cos α
= 5 cos 30◦
= 4.33 m/s
1
y(t) = y0 − gt2 + (v0 sin α)t
2
1
= 5 − (9.8)t2 + (5 sin 30◦ )t
2
= −4.9t2 + 2.5t + 5
g sec2 α 2
y(x) = y0 − x + x tan α
2v02
9.81 sec2 (30) 2
= 5− x + x tan(30)
252
= 5 − 0.1307x2 + 0.5773x
= −0.13x2 + 0.58x + 5
~r(0) = ~r0
x(0)î + y(0)ĵ = 0î + y0 ĵ
implies that
x(0) = 0 (7.4.9)
and
y(0) = y0 (7.4.10)
at t = 0
ẋ(0) = v0 (7.4.12)
ẏ(0) = 0 (7.4.13)
At P , the only force acting on the particle is force of gravity acting in the downward
direction. Then by Newton’s second law of motion its equation of motion is
F~ ~
= W
m~a = −m~g
ẍî + ÿ ĵ = −g ĵ (7.4.14)
ẍ(t) = 0 (7.4.15)
ÿ(t) = −g (7.4.16)
ẋ(t) = A1 (7.4.17)
At t = 0, (7.4.17) becomes
ẋ(0) = A1 (7.4.18)
A1 = v 0 (7.4.19)
ẋ(t) = v0 (7.4.20)
7.4 Projectile Motion of a Projectile Projected from a Relative Height from Origin 183
(7.4.20) gives the horizontal scalar component of velocity of the particle at any time t.
Integrating it with respect to t
x(t) = v0 t + B1 (7.4.21)
At t = 0, (7.4.21) becomes
B1 = 0 (7.4.23)
x(t) = v0 t (7.4.24)
(7.4.24) gives the horizontal component of position of the particle at any time t.
Next for vertical component integrate (7.4.16) with respect to t
At t = 0, (7.4.25) becomes
ẏ(0) = A2 (7.4.26)
A2 = 0 (7.4.27)
(7.4.28) gives the vertical scalar component of velocity of the particle at any time t.
Hence the velocity of the particle at any time t is
1
y(t) = − gt2 + B2 (7.4.30)
2
At t = 0, (7.4.30) becomes
1
y(0) = − g(0)2 + B2 (7.4.31)
2
184 7 Two Dimensional Projectile Motion
B 2 = y0 (7.4.32)
Using (7.4.34), (7.4.29) will give the velocity with which the projectile will hit the ground.
r
2y0
~v (t) = v0 î − g ĵ (7.4.35)
g
From (7.4.24), we can find the time required to reach the particle at P as
x(t)
t = (7.4.36)
v0
Using (7.4.36), (7.4.33) becomes
2
1 x(t)
y(t) = y0 − g (7.4.37)
2 v0
(7.4.37) gives the path of the projectile of the particle at any time t.
Example 7.4.2. A cannonball is fixed at the top of a cliff. The height of the cliff is 20 m
from ground level. A ball is fired from the cannon horizontally with a speed 100 m/s.
x(t) = v0 t
= 100t m (7.4.38)
1
y(t) = y0 − gt2
2
= 20 − 4.9t2 m
r
2y0
t =
g
r
2(10)
=
9.8
= 1.4 s
x(t) = 100(1.4)
= 140 m (7.4.39)
186 7 Two Dimensional Projectile Motion
y
tan β =
x
with P (x, y) is a point common to both planes (see Fig. 7.13). Clearly this line passes
through the origin, so its equation is
y = x tan β (7.4.40)
g sec2 α 2
y = − x + x tan α
2v02
g sec2 α 2
x tan β = − x + x tan α
2v02
7.4 Projectile Motion of a Projectile Projected from a Relative Height from Origin 187
x
|OP | = r =
cos β
= x sec β
188 7 Two Dimensional Projectile Motion
v02 sec2 β
xRu = x sec β = (sin (2α − β) − sin β) (7.4.43)
g
If the plane p is inclined downward to the horizontal, then the angle β is replaced with −β
and its intersection with the vertical plane is the line
y = −x tan β (7.4.44)
The x coordinate of the point where the projectile hits the plane p is given by replacing β
by −β in (7.4.42)
v02 sec(−β)
x = (sin (2α − (−β)) − sin(−β))
g
v02 sec β
xRd = = (sin (2α + β) + sin β) (7.4.45)
g
Time of Flight on an Inclined Plane The time of flight up the inclined plane p is defined
as the time taken by the projectile to achieve the full range. This time can be calculated as
distance
time =
speed
The distance is the horizontal distance x is given by (7.4.41) and speed is horizontal speed
v0 cos α
1 v02 sec β
t = 2 cos α sin (α − β)
v0 cos α g
2v0
= sec β sin (α − β) (7.4.46)
g
Maximum Range on an Inclined Plane Since β is a fixed angle, then for a given value
of v0 , the range up the inclined plane is maximum if
Hence maximum range up the inclined plane is obtained by using (7.4.47) in (7.4.43)
v02
xRu−max = (1 − sin β)
g cos2 β
v02
= (1 − sin β)
g 1 − sin2 β
v02
= (7.4.48)
g (1 + sin β)
v02
xRd−max = (7.4.49)
g (1 − sin β)
First consider the particle is projected in the backward direction, then both velocities have
opposite directions. Consider a cartesian plane as the vertical plane with x axis along
horizontal and y axis along vertical. The point O is the point of projection. Since ~v0 is
~ will be considered negative. The initial velocity ~v0
considered positive in section 7.2.1 so V
can be written as
~ = v0 cos αî + v0 sin αĵ − V î
~v0 = ~v (0) − V (7.5.1)
~r(0) = ~0
x(0)î + y(0)ĵ = 0î + 0ĵ
implies that
x(0) = 0 (7.5.2)
190 7 Two Dimensional Projectile Motion
and
y(0) = 0 (7.5.3)
at t = 0
At P , the only force acting on the particle is force of gravity acting in the downward
direction. Then by Newton’s second law of motion its equation of motion is
F~ ~
= W
m~a = −m~g
ẍî + ÿ ĵ = −g ĵ (7.5.7)
7.5 Projectile Motion with Horizontal Relative Motion 191
ẍ(t) = 0 (7.5.8)
ÿ(t) = −g (7.5.9)
ẋ(t) = A1 (7.5.10)
At t = 0, (7.5.10) becomes
ẋ(0) = A1 (7.5.11)
A1 = v0 cos α − V (7.5.12)
(7.5.13) gives the horizontal scalar component of velocity of the particle at any time t.
Integrating it with respect to t
At t = 0, (7.5.14) becomes
B1 = 0 (7.5.16)
(7.5.17) gives the horizontal component of position of the particle at any time t.
From section 7.2.1 the vertical component of position of the particle at any time t is
1
y(t) = − gt2 + (v0 sin α)t
2
From (7.5.17), we can find the time required to reach the particle at P as
x(t)
t = (7.5.18)
v0 cos α − V
192 7 Two Dimensional Projectile Motion
d2 xR 2v0
= sin α (V − v0 cos α)
dα2 g
For maximum horizontal range, we must have
d2 xR
<0
dα2
7.6 Projectile Motion of a Projectile Projected with a Relative speed from Origin 193
~r(0) = ~0
x(0)î + y(0)ĵ = 0î + 0ĵ
implies that
x(0) = 0 (7.6.2)
and
y(0) = 0 (7.6.3)
at t = 0
At P , the only force acting on the particle is force of gravity acting in the downward
direction. Then by Newton’s second law of motion its equation of motion is
F~ ~
= W
m~a = −m~g
ẍî + ÿ ĵ = −g ĵ (7.6.7)
ẍ(t) = 0 (7.6.8)
ẋ(t) = A1 (7.6.9)
7.6 Projectile Motion of a Projectile Projected with a Relative speed from Origin 195
At t = 0, (7.7.4) becomes
ẋ(0) = A1 (7.6.10)
Using (7.6.5), (7.7.5) implies that
A1 = V + v0 cos α (7.6.11)
Using (7.6.11), (7.7.4) becomes
ẋ(t) = V + v0 cos α (7.6.12)
(7.7.6) gives the horizontal scalar component of velocity of the particle at any time t.
Integrating it with respect to t
x(t) = (v0 cos α + V ) t + B1 (7.6.13)
At t = 0, (7.7.7) becomes
x(0) = (v0 cos α + V ) (0) + B1 (7.6.14)
Using (7.6.2), (7.7.8) becomes
B1 = 0 (7.6.15)
Using (7.7.9), (7.7.7) becomes
x(t) = (V + v0 cos α) t (7.6.16)
(7.7.10) gives the horizontal component of position of the particle at any time t.
Since the motion of the cart is along horizontal only, so at ant time t, vertical component
speed is given by (7.2.21) and vertical component of position is given by (7.2.25). From
(7.7.10), we can find the time required to reach the particle at P as
x(t)
t = (7.6.17)
V + v0 cos α
Time of Flight As the particle is moving under gravity, so it will strike horizontal axis
(or plane) after time t. This time is obtained by taking vertical component of speed equal
to zero. Since there is no change in it so same relation is for time of flight that is given by
(7.2.31)
2v0 sin α
t = tr =
g
Horizontal Range Let the particle hits the x axis at A after projected from O. Then the
distance |OA| is known as horizontal range and is calculated by using (7.2.31) in (7.7.10)
2v0 sin α
x(t) = (V + v0 cos α)
g
2V v0 sin α v02
xr = + sin 2α (7.6.18)
g g
196 7 Two Dimensional Projectile Motion
∆x = xr − xR
2V v0 sin α
= (7.6.19)
g
ẋ(0) = v0 cos α
ẏ(0) = v0 sin α
7.7 Projectile Motion with Air Resistance 197
and
x(0) = 0
y(0) = 0
then at P , two force are acting on the particle, one is force of gravity, acting in the downward
direction and the other is air resistance, opposing the velocity. Then by Newton’s second
law of motion its equation of motion is
F~ ~ + F~r
= W
m~a = −m~g − k~v
ẍî + ÿ ĵ = −g ĵ − k ẋî − k ẏ ĵ
= −k ẋî − (g + k ẏ)ĵ (7.7.1)
ẍ(t) = −k ẋ
ẍ(t) + k ẋ = 0 (7.7.2)
ẇ(t) + kw = 0
w = A1 e−kt
ẋ(t) = A1 e−kt (7.7.4)
At t = 0, (7.7.4) becomes
ẋ(0) = A1 (7.7.5)
A1 = v0 cos α
(7.7.6) gives the horizontal velocity of the particle at any time t. Integrating it with respect
to t
1
x(t) = (v0 cos α)e−kt + B1 (7.7.7)
−k
198 7 Two Dimensional Projectile Motion
At t = 0, (7.7.7) becomes
1
x(0) = − (v0 cos α)(1) + B1 (7.7.8)
k
Using initial condition x(0) = 0, (7.7.8) implies that
1
B1 = (v0 cos α) (7.7.9)
k
Using (7.7.9), (7.7.7) becomes
1 1
x(t) = − (v0 cos α)e−kt + (v0 cos α)
k k
v0 cos α −kt
= 1−e (7.7.10)
k
(7.7.10) gives the horizontal scalar component of position of the particle at any time t. This
time can be calculated as
kx
1 − e−kt = (7.7.11)
v cos α
0
kx
e−kt = 1−
v0 cos α
kx
−kt = ln 1 −
v0 cos α
1 kx
−t = ln 1 − (7.7.12)
k v0 cos α
Next for vertical component consider (7.7.3), replace ẏ by z and ÿ by ż, (7.7.2) becomes
linear first order differential equation
ż(t) + kz = −g
Exercises
1. A footballer running on ground level kicks a football and the ball is projected with a
speed of 10 m/s at an inclination 35◦ with the ground.
3. The maximum range of Shaheen missile from ground to ground mark is 1500 km.
4. A shell bursts on contact with the ground and pieces from it, fly in all directions. If
a piece is flying making an angle π6 with the ground with a speed of 500 m/s.
5. A particle of mass m is projected from the ground with a velocity ~v0 , making an angle
α with the ground. If the air resistance is proportional to the square of the velocity
of the body.
If a particle is moving under the action of a force which is always directed towards or away
from a fixed point and depends on its distance from that point, then the force is called a
central force, the motion is known as central force motion and the fixed point is called the
center of the force. The fixed point is usually taken as the origin. The orbital movement
of planets and satellites are such motions. The laws which govern this motion were first
postulated by Kepler (1571-1630). His interest was in describing the motion of planets
around the sun. He postulated the following laws:
K1 The orbits of the planets are ellipses with the Sun at one focus
K2 The line joining a planet to the Sun sweeps out equal areas in equal intervals of time
K3 The square of the period of a planet is proportional to the cube of the major axis of
its elliptical orbit
F~ (r) = F r̂
where r̂ is a unit vector in the radial direction. Particularly, if the force is inversely pro-
portional to the square of the distance between the particle and the origin, such as the
gravitational force given as,
µm
F~ (r) = − 2 r̂ (8.1.1)
r
203
204 8 Motion Under Central Force
~ = m1~r1 + m2~r2
R (8.1.2)
m1 + m2
Also
Then
and
~+ ~r
~r2 = R (8.1.6)
m1 + m2
Similarly
~− ~r
~r1 = R (8.1.7)
m1 + m2
The attractive force on m1 due to m2 is
m1 m2
F~ = G r̂
r2
where G is the gravitational constant And by Newton’s second law of motion, this force is
F~ = m1~r̈1 (8.1.8)
where p is the linear momentum of the particle. Therefore, the angular momentum L ~ is
always perpendicular to the plane defined by the particle’s position vector ~r and velocity ~v .
The rate of change of the angular momentum L ~ equals the net torque
~
dL
~τ = = ~r × F~
dt
Using (8.2.1)
~
dL
= ~ṙ × m~v + ~r × m~v̇
dt
= ~ṙ × m~ṙ + ~r × m~r̈
= 0 + ~r × F~
~r × F~ = rr̂ × F (r)r̂ = 0
Hence, the angular momentum L ~ is constant. Consequently, the particle’s position ~r and
~ as shown in Fig. 8.3.
velocity ~v always lie in a single plane perpendicular to L,
F~ (r) = m~r̈
208 8 Motion Under Central Force
r2 θ̇ = h (constant) (8.2.5)
h
r2 = (8.2.6)
θ̇
and
h
θ̇ = (8.2.7)
r2
8.2 Motion under Central Force 209
0
F (r) = −U (r)
K2 The line joining a planet to the Sun sweeps out equal areas in equal intervals of time.
That is
dA
= constant (8.2.8)
dt
Proof Consider a planet P moves around the sun S. The trajectory is in the xy plane.
At any time t, P has coordinates P (r, θ) relative to S. At time t + dt, the planet is at Q,
having position vector r + dr as shown in Fig. 8.4. Since dr is very very small, then
dr ∼
= ds = rdθ
Let dA be the area swept by ~r, which is the area of the small region SP Q, that is
1 2
dA = r dθ
2
The areal velocity of the planet is
dA 1 2 dθ
= r
dt 2 dt
1 2
= r θ̇ (8.2.9)
2
210 8 Motion Under Central Force
d2
1 1 µ
+ = (8.2.16)
dθ2 r r h2
1
(8.2.16) is a linear second order nonhomogeneous differential equation for r as a function
of θ. Its general solution is
1 1 1
= +
r r c r p
m2 + 1 = 0
Where A1 and A2 are constants of integration, let A1 = e hµ2 cos ψ and A2 = −e hµ2 sin ψ,
with e and ψ are constants. Then
1 µ
= e 2 cos(θ + ψ)
r c h
1 µ
= (1 + e cos(θ + ψ))
r h2
The constant ψ = 0 and e > 0 if we rotate the base line θ = 0, then the equation describing
the trajectory will be
h2
r = (8.2.17)
µ (1 + e cos θ)
212 8 Motion Under Central Force
h2
l = (8.2.19)
µ
Hence the reduced mass is
h2
µ = (8.2.20)
l
The constant e ≥ 0 is called the eccentricity, defining the conic surface as follows
When e < 1, the trajectory given by (8.2.17) is an ellipse, thus proving Keplers first law.
The point in the trajectory which is closest to the focus is called the periapsis and is denoted
8.2 Motion under Central Force 213
by π. For elliptical orbits, the point in the trajectory which is farthest away from the focus
is called the apoapsis and is denoted by α. When considering orbits around the earth, these
points are called the perigee and apogee, whereas for orbits around the sun, these points are
called the perihelion and aphelion, respectively.
At perihelion and aphelion, the radial velocity is zero, so ~v is orthogonal to the radius vector
~r . At these points the areal velocity of the planet is
dA 1
= rv
dt 2
from the focus S, so is the aphelion. If a is the semi major axis and b is the semi minor axis
of the ellipse, then
a2 = b2 + c2 (8.2.22)
Also the distance between the focus S and the center of the ellipse is
c = ae = a − rπ (8.2.23)
The eccentricity is
c
e = (8.2.24)
a
214 8 Motion Under Central Force
a2 = b2 + (ae)2 (8.2.25)
or
b2 = a2 (1 − e2 ) (8.2.26)
b2
l = (8.2.28)
a
Using (8.2.28), (8.2.26) becomes
l = a(1 − e2 ) (8.2.29)
or
l
a = (8.2.30)
(1 − e2 )
2a = rπ + rα (8.2.31)
2l
2a = rπ + rα = (8.2.32)
(1 − e2 )
Apply partial fraction technique on right hand side of (8.2.32), then we have
l l
2a = rπ + rα = + (8.2.33)
(1 + e) (1 − e)
l
rπ = = a(1 − e) (8.2.34)
(1 + e)
and
l
rα = = a(1 + e) (8.2.35)
(1 − e)
8.2 Motion under Central Force 215
A = πab
p
= πa2 (1 − e2 ) (8.2.36)
Consider (8.2.10)
dA h
=
dt 2
which is first order differential equation and has solution
h
A = T (8.2.37)
2
where T is constant of integration, known as period of the orbit. Equating (8.2.38) and
(8.2.39)
h p
T = πa2 (1 − e2 )
2
2π 2 p
h = a (1 − e2 ) (8.2.38)
T
Using (8.2.29), (8.2.38) can be written as
2
2 2π
h = a3 l (8.2.39)
T
h2
µ =
l
2
2π
= a3
T
or
4π 2
2
T = a3 (8.2.40)
µ
Example 8.2.1. The planet Mercury orbits the Sun in 87.97 days, in an elliptic orbit with
semi major axis is 57.91 × 106 km and semi minor axis is 56.67 × 106 km. Find
2. the eccentricity.
6. its speed
Figure 8.7: The planet Mercury moves in elliptic orbit around the Sun
1. To find the center of the elliptic path, we use (8.2.23), to calculate the distance of the
focus S(0, 0) from the center C(−c, 0) as
8.2 Motion under Central Force 217
c = a − rπ
= 57.91 × 106 − 46.00 × 106
= 11.91 × 106 km (8.2.41)
c
e =
a
11.91 × 106
=
57.91 × 106
≈ 0.2057 (8.2.42)
b2
l =
a
2
56.67 × 106
=
57.91 × 106
= 36.54 × 106 km (8.2.43)
A = πab
= π 56.67 × 106 57.91 × 106
(8.2.44)
15 2
= 8.37 × 10 km (8.2.45)
Consider (8.2.10)
dA A
=
dt T
8.37 × 1015
=
7.6 × 105
= 1.1 × 109 km2 /s (8.2.46)
dA 1
= rv
dt 2
218 8 Motion Under Central Force
6. So its speed
a) at its perihelion is
2 dA
v =
r dt
At perihelion, rπ = 46.00 × 106 km, then the speed is
2 1.1 × 109
v =
46.00 × 106
= 47.8 km/s (8.2.47)
2 1.1 × 109
v =
69.82 × 106
= 31.5 km/s (8.2.48)
8.2 Motion under Central Force 219
Exercise A satellite is launched in a direction parallel to the surface of the earth with
a speed 36, 900 km/h in an elliptic orbit from perigee, 500 km away from the earth. Let
63700 be the radius of the earth, then find
2. the eccentricity.
Small Oscillation
dU
F = − (9.1.1)
dx
Let the particle initially be at rest at a local minimum of U (x). Let this minimum is zero
and it exists at x = x0 . Then
U (x0 ) = 0 (9.1.2)
and
dU
= 0 (9.1.3)
dx0
Let it be given a small kick, from its equilibrium position, so that it moves back and forth
around it, executing simple harmonic motion of small amplitude. To see this, expand U (x)
in a Taylor series around the equilibrium point, x0 .
1 00
U (x) = U (x0 ) + U 0 (x0 )(x − x0 ) + U (x0 )(x − x0 )2
2!
1 000
+ U (x0 )(x − x0 )3 + · · · (9.1.4)
3!
Using (9.1.2) and (9.1.3) in (9.1.4), we are left with the U 00 (x0 ) and higher-order terms.
But for sufficiently small displacements, these higher-order terms are negligible compared
to the U 00 (x0 ) term, and we are left with
1 00
U (x) ∼
= U (x0 )(x − x0 )2 (9.1.5)
2
221
222 9 Small Oscillation
a) Equilibrium position.
b) Angular frequency.
Solution Using (9.1.3), the equilibrium point, x0 is
U 0 (x) = 0 (9.1.10)
A B
− 3+ 2 = 0 (9.1.11)
x x
A
x = 2 (9.1.12)
B
or
A
x0 = 2 (9.1.13)
B
The spring constant is given by using Using (9.1.7). First the second derivative of U (x) is
A B
U 00 (x) = 4
− 3 (9.1.14)
x x
Using (9.1.13),(9.1.14) becomes,
r
U 00 (x0 )
ω =
r m
B4
= (9.1.15)
8mA3
(9.1.15) gives the frequency of oscillation.
9.2 Equilibrium and Stability 223
F = 0
0
−U (x) = 0
or
U 0 (x) = 0
U 0 (x) = 0
U 00 (x) > 0
U 0 (x) = 0
U 00 (x) < 0
Example 9.2.1. In spring mass system with k > 0, the potential energy function is
1 2
U (x) = kx
2
and
U 0 (x) = kx
then
U 0 (x) = 0
gives
x = 0
224 9 Small Oscillation
Next
U 00 (x) = k
U 00 (0) = k > 0
Force
In this section, the application of small but fast oscillation is presented accompanied with
the stability of the system. This stabilization was initialized by Stephenson in 1908. Later,
in 1950 Kapitza renewed it. In 1961 Landau introduced oscillation under the harmonic force.
For stabilization, the potential energy function is obtained by using averaging procedure.
For this consider one dimensional motion of a particle of mass m under the action of a time
independent potential field force
−dU
F1 = (9.3.1)
dx
and a small but fast oscillating force
The system has time period TU due to force F1 , so may oscillate due to this force with
frequency ωU = T1U and T is the time period due the force F2 , and oscillation due to this
force is ω = T1 . The fast oscillation means the frequency ω ωU . Also the coefficients f1 , f2
are functions of co-ordinates only. The magnitude of F2 is not assumed small in comparison
with the force F1 but we shall assume that the oscillation of the particle (denoted below by
ξ) due to this force is small.
Since the system has one dimensional motion it means it depends only on the space co-
ordinate x. Then by Newton’second law of motion, equation of motion of the particle
is
mẍ = F1 + F2 (9.3.3)
Since the particle transverse a slow path and at the same time execute fast but small
oscillations of frequency ω about the path, the function x(t) (path) is a sum of both parts
as shown in Fig. 9.5.17. If we represent X(t) the slow path and ξ(t) the fast but small
oscillations, the path is
x(t) = X(t) + ξ(t) (9.3.4)
By (9.3.4) the functions f and U will also be transformed. Using Taylor’s expansion in
9.3 Motion in a Rapidly Oscillating Field under Harmonic Force 225
dU d d ∂U dU d2 U
= [U (x) = U (X + ξ)] = [U (X) + ξ ]= +ξ (9.3.5)
dx dX dX ∂X dX dX 2
substituting Eq.(9.3.5) in Eq.(9.3.3) we have
dU d2 U df
mẌ + mξ¨ = − −ξ 2
+ f (X, t) + ξ (9.3.6)
dX dX dX
This equation involves both fast and slow motions, which must be separately equal. For
the fast oscillating term we can put simply
The term C1 = ξ˙0 is the initial speed of the small oscillation. Let at t = 0, ξ˙0 = 0, then
above relation is
1
mξ˙ = (−f1 sin ωt + f2 cos ωt)
ω
Again integrating, we have
1
mξ = − (f1 cos ωt + f2 sin ωt) + C2
ω2
The term C2 = ξ0 is the initial position of the small oscillation. Let at t = 0, ξ0 = 0, then
above relation is
f
mξ = − 2
ω
or
f
ξ=− (9.3.9)
mω 2
(9.3.9) gives the small but fast oscillation. Next the mean value of a T periodic function is
denoted by bar and is given by
1 T
Z
¯
f = f (t)dt (9.3.10)
T 0
Using Eq. (9.3.10), the mean value of f (t) over its period T = 2π ω is
Z T
1
f¯ = (f1 cos ωt + f2 sin ωt) dt
T 0
T
1 f1 f2
= − sin ωt + cos ωt
T ω ω 0
1 h i
= − f1 [sin(2π) − sin(0)] + f2 [cos(2π) − cos(0)]
2π
= 0
Hence the mean value of f (t) over its period T is zero. Since ξ(t) contains f , so its mean
value over same period is also zero. Also
¯
ξ¨ = 0
During this time averaging, the function X(t) remains invariant. We therefore have x =
X(t), i.e. X(t) describes the slow motion of the particle averaged over the rapid oscillations.
Using this time averaging technique, we shall derive a relation which will be the function
X(t) only.
Next time averaging (9.3.6) over its period (0 → T ). As calculated above, the mean values
of the first powers of f and ξ over this period are zero, so the mean value of (9.3.6) over
this period is
dU df
mẌ = − −ξ
dX dX
dU 1 df
= − − 2
f (9.3.11)
dX mω dX
9.3 Motion in a Rapidly Oscillating Field under Harmonic Force 227
This relation represent smooth motion of the particle averaged over small oscillation. Con-
sider
1 df 2 df
= f (9.3.12)
2 dX dX
Using Eq. (9.3.12), Eq. (9.3.11) becomes
dU 1 1 df 2
mẌ = − − (9.3.13)
dX mω 2 2 dX
(9.3.13) is a function of X only so can be written as
2
!
d f
mẌ = − U+ (9.3.14)
dX 2mω 2
2
f
In (9.3.14) the quantity U + 2mω 2
can be regarded as potential energy function. We call
it effective potential energy function and is given by
2
f
Uef f = U+ (9.3.15)
2mω 2
Using (9.3.15), (9.3.14) can be written as
dUef f
mẌ = − (9.3.16)
dX
2
(9.3.16) represents that after averaging the net force is conservative. To calculate f we
first calculate the quantity f 2 from (9.3.2) as
To simplify this calculation we make use of the following trigonometric orthogonal relations.
1 T
Z
sin ωt cos ωt = sin ωt cos ωtdt = 0
T 0
1 T 1 T 1
Z Z
2 2 1
cos ωt = cos ωtdt = (1 + cos 2ωt)dt =
T 0 T 0 2 2
Z T Z T
1 1 1 1
sin kω 2 t = sin2 ωtdt = (1 − cos 2ωt)dt =
T 0 T 0 2 2
228 9 Small Oscillation
f12 + f22
Uef f = U + (9.3.20)
4mω 2
Thus the motion of the particle averaged over the small but fast oscillations is the same
as if the constant potential U were augmented by a constant quantity proportional to the
squared amplitude of the variable field.
The system moving under the action small but fast oscillation will be stabilized by mini-
mizing (9.3.20).
An external force sin ωt of high frequency with amplitude a is applied at the pivot point to
have a horizontal modulation, as shown in Fig. (9.5). This high frequency means ω ω0 .
Then at any time the position of the particle is
x = a sin ωt + l sin φ
y = l cos φ
We can write down the equation of motion of the above system by finding Euler-Lagrange
equations. The time independent potential energy function is given by (9.4.1), so we need
9.4 Stabilization of Modulated Pendulum with Harmonic Force 229
to find its kinetic energy only. For this the velocity components are
L = K −U
mφ̇2 l2 ma2 ω 2
= + cos2 ωt + maω φ̇l cos ωt cos φ + mgl cos φ
2 2
1
cos2 ωt = (1 + cos2ωt)
2
d t sin2ωt
= +
dt 2 4ω
230 9 Small Oscillation
Also consider
d
(sin φ cos ωt) = −ω sin ωt sin φ + φ̇ cos ωt cos φ
dt
then
d
φ̇ cos ωtcosφ = (sin φ cos ωt) + ω sin ωt sin φ
dt
then Lagrangian becomes
mφ̇2 l2 ma2 ω 2 d t
sin2ωt
L = + + + mgl cos φ
2 2 dt 2 4ω
d
+ mlaω (sin φ cos ωt) + ω sin ωt sin φ
dt
mφ̇2 l2
= + mgl cos φ + mlaω 2 sin ωt sin φ
2
ma2 ω 2 d t
sin 2ωt d
+ + + mlaω (sin φ cos ωt)
2 dt 2 4ω dt
Which is of the form
d
L0 = L + f (φ, t)
dt
d
Since Lagrangian remains invariant if we add or subtract dt f (q, t). So omitting total deriva-
tives we have
mφ̇2 l2
L = + mgl cos φ + mlaω 2 sin ωt sin φ
2
Here q = lφ is the only generalized coordinates. So equation of motion is
!
∂L d ∂L
− =0 (9.4.2)
∂(lφ) dt ∂(lφ) ˙
∂L 1 ∂L
=
∂(lφ) l ∂(φ)
1
= [mgl(− sin φ) + mlaω 2 sin ωt cos φ]
l
= m[g(− sin φ) + aω 2 sin ωt cos φ]
! !
∂L 1 ∂L
=
˙
∂(lφ) l ∂(φ) ˙
1
= [ml2 φ̇]
! l
d ∂L
= mlφ̈
dt ˙
∂(lφ)
9.4 Stabilization of Modulated Pendulum with Harmonic Force 231
a2 ω 2
Uef f 2
Uef f =
e = − cos φ + cos φ (9.4.6)
mgl 4gl
Stability
The positions of stable equilibrium correspond to the minima of Uef f in [−π, π]. The first
step is to differentiate (9.4.5) with respect to φ.
2a2 ω 2
dUef f
= mgl sin φ − sin φ cos φ
dφ 4gl
a2 ω 2
= mgl sin φ 1 − cos φ (9.4.7)
2gl
dUef f
To find the critical points we have dφ = 0, which means
a2 ω 2
mgl sin φ 1 − cos φ = 0
2gl
Since m 6= 0, g 6= 0 and l 6= 0. That means
sin φ = 0
or
a2 ω 2
1− cos φ = 0
2gl
232 9 Small Oscillation
First consider
sin φ = 0
⇒ φ = 0, π
Next consider
a2 ω 2
1− cos φ = 0
2gl
a2 ω 2
cos φ = 1
2gl
2gl
cos φ =
a2 ω 2
2gl
φ = ± arccos
a2 ω 2
(2) Stability at φ = π
The system at φ = π is stable if Uef f is minimum at φ = π. This minimum exist if
d2 Uef f
|φ=π > 0
dφ2
Since cos(π) = −1
a2 ω 2
⇒ mgl −1 − > 0
2gl
a2 ω 2
− 1+ > 0
2gl
Which is impossible, hence the system is unstable φ = π.
(3) Stability at φ = ± arccos a2gl2 ω2
−2gl a2 ω 2
mgl + >0
a2 ω 2 2gl
or
−2gl a2 ω 2
+ >0
a2 ω 2 2gl
a2 ω 2 2gl
> 2 2
2gl a ω
It means that
a2 ω 2 > 2gl
2gl
Hence if a2 ω 2 > 2gl the position given by cos φ = a2 ω 2
is stable. It is also conditional stable
position.
All these stable points are shown in Fig. 9.3.
We can write the summery as:
234 9 Small Oscillation
L = K −U
m 22 m 2 2d t sin2ωt d
= φ̇ l + a ω + + mlaω ω cos φ sin ωt + (cos ωt cos φ)
2 2 dt 2 4ω dt
mga d
+ mglcosφ − (sinωt)
ω dt
Omitting total derivative
m 22
L = φ̇ l + mlaω 2 cos φ sin ωt + mgl cos φ
2
Here lφ is the only generalized coordinate
∂L 1 ∂L
=
∂lφ l ∂φ
1
−mlaω 2 sin φ sin ωt − mgl sin φ
=
l
−maω 2 sin φ sin ωt − mg sin φ
=
∂L ∂L
=
∂lφ̇ l∂ φ̇
1
= mφ̇l2 = mlφ̇
l
d ∂L
= mlφ̈
dt ∂lφ̇
Hence the equation of motion is
comparing it with Eq. (9.3.3) the external force acting on the bob is
a2 ω 2
mglsinφ 1 + cosφ = 0
2gl
Since m 6= 0, g 6= 0 and l 6= 0. That means
a2 ω 2
sinφ 1 + cosφ = 0
2gl
so we have
sinφ = 0
or
a2 ω 2
1+ cosφ = 0
2gl
First consider
sinφ = 0
⇒ φ = 0, π
Next consider
a2 ω 2
1+ cosφ = 0
2gl
a2 ω 2
cosφ = −1
2gl
2gl
cosφ = −
a2 ω 2
9.4 Stabilization of Modulated Pendulum with Harmonic Force 237
and
2gl
φ = ± arccos − 2 2
a ω
d2 Uef f a2 ω 2
= mgl cos φ + cos 2φ (9.4.13)
dφ2 2gl
(1) Stability at φ = 0
d2 Uef f
|φ=0 > 0
dφ2
Since cos(0) = 1
a2 ω 2
⇒ mgl 1 + > 0 always > 0
2gl
(2) Stability at φ = π
d2 Uef f
|φ=π > 0
dφ2
Since cos(π) = −1
a2 ω 2
⇒ mgl −1 + > 0
2gl
a2 ω 2
−1 + > 0
2gl
It means that
a2 ω 2 > 2gl
Hence if a2 ω 2 > 2gl the position given by φ = π is stable. It is conditional stable position.
238 9 Small Oscillation
2gl
(3) Stability at φ = ± arccos a2 ω 2
d2 Uef f −2gl a2 ω 2 4g 2 l2
= mgl + 2· 4 4 −1
dφ2 a2 ω 2 2gl a ω
2gl
= mgl −1
a2 ω 2
2gl
Now | cosφ |= a2 ω 2
<1
d2 Uef f
then we have dφ2
<0
so the position given by cosφ = a−2gl
2 ω 2 is not stable.
All these stable points are shown in Fig. 9.6.
We can write the summery as:
(1) the position φ = 0(vertically downward position) is always stable. Graphically min-
imization of Uef f at φ = 0 is shown in Fig. 9.7 (a) .
(2) If a2 ω 2 > 2gl then φ = π (vertically upward position) is also stable. Graphically
minimization of Uef f at φ = π is shown in Fig. 9.7 (b)
−2gl
(3) the position given by cos φ = a2 ω 2
is unstable.
9.5 Motion in a Rapidly Oscillating Field under Arbitrary Periodic Force with Zero Mean
239
−dU
F1 = (9.5.1)
dx
and a fast oscillating periodic force with zero mean. The Fourier expansion of this periodic
force is
∞ h
X i
f (x, t) = ak (x) cos(kωt) + bk (x) sin(kωt) , (9.5.2)
k=1
The system has time period TU due to force F1 , so may oscillate due to this force with
frequency ωU = T1U and T is the time period due the force F2 , and oscillation due to this
force is ω = T1 . The fast oscillation means the frequency ω ωU . The coefficients ak and bk
are functions of co-ordinates only given by (9.5.3) and (9.5.4) respectively.
2 T
Z
ak (x) = f (x, t) cos kωt dt ; (9.5.3)
T 0
2 T
Z
bk (x) = f (x, t) sin kωt dt . (9.5.4)
T 0
dU (x)
mẍ = − + f (x, t) . (9.5.5)
dx
We present the movement as a slow path and at the same time execute fast but small
oscillations of frequency ω about the path: x(t) = X(t) + ξ(t), where ξ(t) corresponds to
these small oscillations. The mean value of the function ξ(t) over its period T is zero, and
the function X(t) changes only slightly in that time. Denoting this average by a bar, we
therefore have x̄ = X(t). Now Taylor’s expansion in powers of ξ up to the first order term
provides us:
dU dU d2 U
= +ξ . (9.5.6)
dx dX dX 2
Substituting (9.5.6) in (9.5.5) we have:
dU d2 U df
mẌ + mξ¨ = − −ξ 2
+ f (X, t) + ξ . (9.5.7)
dX dX dX
240 9 Small Oscillation
This equation involves both fast and slow terms, which must evidently be separately equal.
For the fast term we can put simply
dU d2 U df
mẌ = − −ξ 2
+ξ .
dX dX dX
Integrating (9.5.8) with the function f given by (9.5.2) and regarding X as a constant, we
get
∞
1 X 1
ξ=− a k cos kωt + bk sin kωt
mω 2 k2
k=1
Next we average equation (9.5.7) with respect to the time interval [0, T ]: Since ξ¯ = 0 and
f¯ = 0,
dU df
mẌ = − +ξ (9.5.9)
dX dX
and
∞
df X
= a´k cos kωt + b´k sin kωt ,
dX
k=1
where a´k = dak /dX and b´k = dbk /dX. Then we apply the time averaging:
∞
df 1 X ak a´j
ξ· =− · cos kωt cos jωt+
dX mω 2 k2
k,j=1
bk a´j ak b´j
+ · sin kωt cos jωt + · cos kωt sin jωt +
k2 k2 #
bk b´j
+ 2 · sin kωt sin jωt .
k
Since
and
0 k 6= j
cos kωt cos jωt = sin kωt sin jωt =
1 k=j
2
(9.5.11)
9.5 Motion in a Rapidly Oscillating Field under Arbitrary Periodic Force with Zero Mean
241
then we have
∞
1 X 1 da2k db2k
df
ξ· =− + . (9.5.12)
dX 4mω 2 k 2 dX dX
k=1
dUef f
mẌ = − ,
dX
where the effective potential energy is defined as
1
∞
X a2k + b2k
Uef f = U + . (9.5.14)
4mω 2 k2
k=1
and shown in Fig. 9.8. The Fourier expansion of (9.5.15) is given by (9.5.16)
∞ h
X i
2
f = mω cos φ ak (x) cos(kωt) + bk (x) sin(kωt) (9.5.16)
k=1
This force is applied at the pivot point to have a horizontal modulation, as shown in Fig.
(9.5.1). Since the first Fourier coefficient is
242 9 Small Oscillation
Z T
2
a0 (x) = f (x, t) dt (9.5.17)
T 0
Hence the zero meaning f = 0 is equivalent to show a0 = 0. The first Fourier coefficient of
(9.5.15) is
2 T
Z
a0 = sin ωt dt
T 0
Using (9.5.10) we have
a0 = 0 (9.5.18)
Hence the mean value of f is zero. The other Fourier coefficients are given by (9.3.3)
2 T
Z
ak = sin ωt cos kωt dt
T 0
Using (9.5.11) we have
ak = 0
The other fourier coefficient is
Z T
2
bk = sin ωt sin kωt dt
T 0
Again using (9.5.11) we have the result for k = 1 only
2 T
Z
b1 = sin2 ωt dt
T 0
= 1 (9.5.19)
Hence with Fourier coefficient
bk = mω 2 cos φ (9.5.20)
the force acting on the particle is
f = mω 2 cos φ sin(ωt) (9.5.21)
Using (9.5.14) the effective potential energy is
ω2
2
Uef f = mgl − cos φ + cos φ .
4gl
Which is same as given by (9.4.5) with a = 1. Hence the stability of extremum positions
φ = 0 , π , ± arccos 2gl/ω 2 is as following.
(1) If ω 2 < 2gl , the downward position φ = 0 is stable.
(2) The Vertically upward position φ = π is unstable.
2gl
(3) If ω 2 > 2gl the position given by cos φ = ω2
is stable.
The stable points with horizontal modulation are shown in Fig. 9.10. The same we will
observe for vertical modulation.
9.5 Motion in a Rapidly Oscillating Field under Arbitrary Periodic Force with Zero Mean
243
where the function R is T -periodic R(t + T, n) ≡ R(t, n). Let’s introduce the triangular
shape linear force:Rs (t) given by (10.3.2)
4t T
T if 0 ≤ t < 4 ;
4 T T 3T
RS (t) = T 2 −t if 4 ≤t< 4 ; (9.5.22)
4(t−T )
3T
if ≤t<T .
T 4
Whic is equivalent to show that R̄s = 0, and we can proceed further. The other Fourier
coefficients are given by Eq.(9.3.3)
Z T
2
ak = RS (t) cos kωt dt
T 0
244 9 Small Oscillation
2 4 t sin kωt cos kωt T4 sin kωt 3T4 4 t sin kωt cos kωt 3T4
= + 2 2 +2 − + 2 2
T T kω k ω kω T4 T kω k ω
T
0 4
2 4 t sin kωt cos kωt T sin kωt T
+ + 2 2 3T − 4
T T kω k ω kω 3T4
4
" ! #
T sin kω T4 cos kω T4
2 4 2 3T T
= + +sin kω − sin kω
T T 4 kω k2 ω2 kω 4 4
" ! !#
2 4 3T sin kω 3T
4 T sin kω T4 cos kω 3T
4 − cos kω 4
T
− − +
TT 4 kω 4 kω k2 ω2
" ! !#
2 4 sin kωT 3T sin kω 3T4 cos kωT − cos kω 3T 4
+ T − +
TT kω 4 kω k2 ω2
" #
2 sin kωT − sin kω 3T4
− 4
T kω
" #
2 8 cos kω T4 4 8 cos kω 3T
4 4 cos kωT
= − − +
T T k2 ω2 T k2 ω2 T k2 ω2 T k2 ω2
2π
using T = ω , we have
ak = 0 (9.5.24)
2 4 −t cos kωt sin kωt T4 − cos kωt 3T4 4 −t cos kωt sin kωt 3T4
= + 2 2 +2 T − + 2 2
T T kω k ω kω T kω k ω
T
0 4 4
2 4 −t cos kωt sin kωt T cos kωt T
+ + 2 2 3T + 4
T T kω k ω kω
3T
4 4
9.5 Motion in a Rapidly Oscillating Field under Arbitrary Periodic Force with Zero Mean
245
" ! #
T cos kω T4 sin kω T4
2 4 2 3T T
= − + + − cos kω + cos kω
T T 4 kω k2 ω2 kω 4 4
" ! !#
2 4 3T cos kω 3T
4 T cos kω T4 sin kω 3T
4 − sin kω 4
T
+ − −
T T 4 kω 4 kω k2 ω2
" ! !#
2 4 cos kωT 3T cos kω 3T4 sin kωT − sin kω 3T
4
+ −T − +
T T kω 4 kω k2 ω2
" #
8 cos kωT − cos kω 3T 4
−
T kω
2π
using T = ω , we have
8 h πi
bk = sin k
k2 π2 2
when k is even, bk = 0
for odd bk can be written as;
8 (−1)j
b2j+1 =
π 2 (2j + 1)2
8 (−1)j
bk = mω 2 cos φ (9.5.25)
π 2 (2j + 1)2
Then the force acting on the particle is
∞
8 X (−1)j
2 2π(2j + 1)t
RS (t) = mω cos φ 2 sin .
π (2j + 1)2 T
j=0
This force also satisfy the property R̄l = 0 as the Fourier coefficient a0 = 0. For horizontal
modulation, the force acting on the particle is
ak = 0 (9.5.27)
4
bk = mω 2 cos φ (9.5.28)
(2k − 1)π
For horizontal modulation the the stability results at nontrivial positions are summarized in
table 9.1. At nontrivial positions, we see that when we apply external triangular type force,
whose area under the curve is less as compared to area under sine curve, we stabilize the
pendulum with relatively high frequency, as compared to external sine force. And when we
apply external rectangular linear force, whose area under the curve is more as compared to
area under sine curve. Hence we can stabilize the pendulum with relatively low frequency
as compared to external sine force.
9.5 Motion in a Rapidly Oscillating Field under Arbitrary Periodic Force with Zero Mean
247
Here the vertical downward position φ = 0 is always stable, the inverse point φ = π is
stable under the condition ω 2 > 2gl .
In place of external harmonic force, we use an arbitrary periodic force. Then the force
acting on the particle is
f = mω 2 sin φ · R(t, n) ,
where R(t, n) represents an arbitrary periodic force.
For vertical modulation the point φ = 0 is always stable while for upper point φ = π we
reproduce conditions as:
Exercise
1. Discuss the stability of a pendulum in the following cases using 9.4
(a) Modulation with a high frequency horizontal small oscillating cos ωt external
force.
(b) Modulation with a high frequency vertical small oscillating cos ωt external force.
2. Discuss the stability of a modulated pendulum with a high frequency horizontal small
oscillating following external force using 9.5.14
(a) f (t) = cos ωt
(b)
4t
T if 0 ≤ t < T4 ;
4 T
t if T4 ≤ t < 3T
f (t) = T 2 − 4 ;
4(t−T ) 3T
if ≤t<T
T 4
(c)
1
2 if 0 ≤ t < 61 T ;
1
≤ t < 13 T ;
1 if 6T
1 1
≤ t < 12 T ;
2 if 3T
f (t) =
− 12 if 1
2T ≤ t < 23 T ;
2
−1 ≤ t < 56 T ;
if 3T
−1 if 5
≤t<T
6T
2
(d)
8t T
T if 0 ≤ t < 8 ;
T 3T
1 if 8 ≤t< 8 ;
8 T 3T 5T
f (t) = T 2 − t if 8 ≤t< 8 ;
5T 7T
−1 if ≤t< ;
8 8
8(t−T )
7T
T if 8 ≤t<T
(e)
8t 3T
T if 0 ≤ t < 8 ;
8 T 3T
< 5T
f (t) = T 2 − t if 8 ≤t 8 ;
8(t−T ) 5T
T if 8 ≤t<T
9.5 Motion in a Rapidly Oscillating Field under Arbitrary Periodic Force with Zero Mean
249
(f)
1 0≤t<τ
f (t, n) =
−(n − 1) τ ≤ t < T , n > 1 ∈ Z.
3. Discuss the stability of a modulated pendulum with a high frequency vertical small
oscillating following external force using 9.5.14
(c)
1
2 if 0 ≤ t < 61 T ;
1
≤ t < 13 T ;
1 if 6T
1 1
≤ t < 12 T ;
2 if 3T
f (t) =
− 12 if 1
2T ≤ t < 23 T ;
2
−1 ≤ t < 56 T ;
if 3T
−1 if 5
≤t<T
6T
2
(d)
8t T
T if 0 ≤ t < 8 ;
T 3T
1 if 8 ≤t< 8 ;
8 T 3T 5T
f (t) = T 2 − t if 8 ≤t< 8 ;
5T 7T
−1 if ≤t< ;
8 8
8(t−T )
7T
T if 8 ≤t<T
(e)
8t 3T
T if 0 ≤ t < 8 ;
8 T 3T
< 5T
f (t) = T 2 − t if 8 ≤t 8 ;
8(t−T ) 5T
T if 8 ≤t<T
250 9 Small Oscillation
(f)
1 0≤t<τ
f (t, n) =
−(n − 1) τ ≤ t < T , n > 1 ∈ Z.
9.5 Motion in a Rapidly Oscillating Field under Arbitrary Periodic Force with Zero Mean
251
Figure 9.4: (a): Stable point φ = 0 of Kapitza Pendulum with Horizontal Oscillation under
the condition ω 2 < 2gl
(b): Stable points given by cosφ = 2glω2
of Kapitza Pendulum with Horizontal Oscillation
2
under the condition ω > 2gl
252 9 Small Oscillation
Figure 9.7: (a): Stable point φ = 0 of Kapitza Pendulum with vertical Oscillation
(b): Stable points given by φ = π of Kapitza Pendulum with Vertical Oscillation under the
condition ω 2 > 2gl
254 9 Small Oscillation
Figure 9.14: (a): Stable point φ = 0 of Kapitza Pendulum with Horizontal Oscillation under
the condition ω 2 < 1.2159gl
(b): Stable points given by cosφ = 1.2159gl
ω2
of Kapitza Pendulum with Horizontal Oscillation
under the condition ω 2 > Sn gl
258 9 Small Oscillation
Figure 9.15: (a): Stable point φ = 0 of Kapitza Pendulum with vertical Oscillation
(b): Stable points given by φ = π of Kapitza Pendulum with Vertical Oscillation under the
condition ω 2 > Sn gl
Chapter 10
Rotation is very common in our daily life. We can see it in fan, clock (with needles), swing,
merry go round, speedometer and dvd in use. In rotation the orientation of a point or line
is fixed, if a point is fixed, we say rotation is about a point and if a line is fixed, we say
rotation is about a line. The line is known as the axis of rotation. In rotation the angular
displacement is called angle of rotation. For example, if we move our arm, we have rotation
about our shoulder (may be consider as a point) and if we close or open the door, there is
rotation about fixed side, can be regarded as the axis of rotation. In chapter 5, we have
already discussed angular motion. In this chapter, we will discuss some particular aspects
of rotation in 2 − space (Rotation around origin) and rotation in 3 − space (Rotation about
an axis passing through origin). First consider rotation in 2 − space.
x = r cos φ (10.1.1)
y = r sin φ (10.1.2)
259
260 10 Rotation and Rigid Bodies
(10.1.7) and (10.1.8) give the coordinates of a point after a rotation θ about origin. Above
equations in matrix form can be written as
0
x cos θ − sin θ x
= (10.1.9)
0
y sin θ cos θ y
cos θ − sin θ
Here R = is known as rotational matrix. Its action on the coordinates
sin θ cos θ
is:
0
x l11 l12 x
= (10.1.10)
0
y l21 l22 y
~
~x = RX (10.1.11)
10.1 Rotations in 2 − space (Rotation about point) 261
(10.1.11) is known as vector equation. It gives the relation how the rotational matrix R
~ The inverse transformation of (10.1.9) is
connects the vectors ~x and X.
0
x cos φ sin φ x
= (10.1.12)
0
y − sin φ cos φ y
If the coordinates of a point are given before rotation, we will use (10.1.9) to calculate its
coordinates after rotation, and if we are given coordinates of a point after rotation, we will
use (10.1.12) to calculate its coordinates before rotation.
From above discussion a general result can be stated as:
If the rotation about origin is anticlockwise, then coordinates of a point can be calculated by
using (10.1.9) and if the rotation about origin is clockwise, then coordinates of a point can
be calculated by using (10.1.12)
Example 10.1.1. Let a point (2, 4) is rotated about the origin through an angle of θ = 30◦ .
Solution Here the rotation is anticlockwise, so we will use (10.1.7) and (10.1.8) to find
new coordinates of A. Given A = (2, 4) and θ = 30◦ so x = 2 and y = 4, then
√
3
cos θ = cos 30 = (10.1.13)
2
and
1
sin θ = sin 30 = (10.1.14)
2
we obtain
√ !
0 3 1
x = 2 −4
2 2
√ !
0 1 3
y = 2 +4
2 2
Here we observe that there is same coordinate transformation or rotational matrix when a
point has anticlockwise rotation about origin and coordinate axes have clockwise rotation
about origin.
Example 10.1.2. If the coordinate axes are rotated about the origin through an angle of
and
1
sin θ = sin 30 =
2
we obtain
√ !
0 3 1
x = 2 +4
2 2
√ !
0 1 3
y = −2 +4
2 2
Example 10.1.3. Suppose that the axes of an x́ý−coordinate system are rotated through
in xy−coordinates.
Solution Here we have θ = 45◦ , then
1
cos θ = cos 45 = √
2
and
1
sin θ = sin 45 = √
2
Using the rotation equations in (10.1.7) and (10.1.8), we obtain
0 x y
x = √ − √
2 2
0 x y
y = √ + √
2 2
Substituting these into equation (10.1.28) yields
y 2 y 2
x x y x y x
√ −√ − √ −√ √ +√ + √ +√ −6 = 0
2 2 2 2 2 2 2 2
x2 y 2
+ = 1 (10.1.28)
12 4
Example 10.1.4. Suppose that the axes of an xy−coordinate system are rotated through
x2 + y 2 = a2 (10.1.29)
in x́ý−coordinates.
Solution Since we have given equation of a circle before rotation, so will will use inverse
transformation
x = x́ cos θ + ý sin θ
y = −x́ sin θ + ý cos θ
in given equation.
x2 + y 2 = a2
(x́ cos θ + ý sin θ)2 + (−x́ sin θ + ý cos θ)2 = a2
x́2 cos2 θ + ý 2 sin2 θ + 2x́ý cos θ sin θ + x́2 sin2 θ + ý 2 cos2 θ − 2x́ý cos θ sin θ = a2
x́2 + ý 2 = a2
It means that x is transformed into x́ and y into ý and the equation of a circle remains
invariant.
266 10 Rotation and Rigid Bodies
Remark 10.1.1. A circle centred at origin is same after any rotation about the origin.
Rotational matrices are orthogonal (i.e. they preserve lengths and angles), hence if R−1 is
inverse and RT is transpose of R then the criteria for a matrix to be rotational is
a) detR = |R| = +1
b) RRT = RT R = I or R−1 = RT
If any of above conditions violate, the matrix is not rotational matrix.
direction about the origin through an angle θ relative to xy−coordinate system. This
10.2 Rotational Matrix in 2 − space 267
rotation is shown in Fig. 10.5. In (10.2.1), lij are cosines of the angles between the axes.
Then rotation matrix can be written as
0
0
cos Ox, Ox cos Ox, Oy
R = 0
0
cos Oy, Ox cos Oy, Oy
In Fig.
10.5, the
angles between
0 the axes are as following.
0
∠ Ox , Ox = θ ∠ Ox , Oy = 90 − θ
0 0
∠ Oy , Ox = 90 + θ ∠ Oy , Oy = θ
Then the rotational matrix is
cos (θ) cos (90 + θ)
R =
cos (90 − θ) cos (θ)
cos θ − sin θ
=
sin θ cos θ
Example 10.2.1. Determine whether the following matrices are rotational matrices.
√
3 1
2 −2
1) A = √
1 3
2 2
1 −2
2) B =
1
2 1
268 10 Rotation and Rigid Bodies
3 1
3) C =
5 2
a) |R| = 1
b) RRT = RT R = I or R−1 = RT
1) The determinant of A is
√
3
− 12
2 3 1
|A| = √ =
4 + 4 =1
1 3
2 2
√
3 1
2 √2
The transpose of A is AT = Next AAT is
3
− 21 2
√ √
3 3
2 − 21 2
1
AAT = √ √2
1 3 3
2 2 − 12 2
√ √
3 1 3 3
+ − 1 0
= √4 4
√ 4 4 =
3 3 3 1
4 − 4 4 + 4 0 1
= I
2) The determinant of B is
1 −2
|B| =
=1+1=2
21 1
3) The determinant of C is
3 1
|C| = =6−5=1
5 2
10.3 Rotations in 3 − space 269
3 5
The transpose of C is C T = Next CC T is
1 2
3 1 3 5
CC T =
5 2 1 2
9 + 1 15 + 2 10 17
= =
15 + 2 25 + 4 17 29
6= I
As one of the conditions is not satisfied, the matrix C is not rotational matrix.
l11 l12 l13
R=
l21 l22 l23
in other
words
0 0 0
cos OX1 , O x1 cos OX1 , O x2 cos OX1 , O x3
0 0 0
R = cos OX2 , O x1 cos OX2 , O x2 cos OX2 , O x3
0 0 0
cos OX3 , O x1 cos OX3 , O x2 cos OX3 , O x3
Determinant of the matrix R denoted by |R| is the volume of the rectangular parallelepiped
0
whose sides have unit length and lie along the axes of the trihedral O x1 x2 x3 . Since both
the trihedral are right-handed, so this volume is +1
|R| = +1 (10.3.3)
The matrix R is called rotational matrix. If RT denotes the transpose of R, the orthogonality
conditions (10.3.1) show that
RRT = RT R = I (10.3.4)
10.3 Rotations in 3 − space 271
RT = R−1 (10.3.5)
(10.3.4) and (10.3.5) confirm that the rotational matrix R is a proper orthogonal matrix.
0
For next discussion, we consider O and O coincide, then rotational matrix is
cos (∠X1 Ox1 ) cos (∠X1 Ox2 ) cos (∠X1 Ox3 )
R=
cos (∠X2 Ox1 ) cos (∠X2 Ox2 ) cos (∠X2 Ox3 )
cos (∠X3 Ox1 ) cos (∠X3 Ox2 ) cos (∠X3 Ox3 )
a) |R| = 1
b) RRT = RT R = I or R−1 = RT
a) The determinant of A is
272 10 Rotation and Rigid Bodies
√
1
2 − 34 43
√
|A| = 34 5 3
8 8
√3 √3 7
−
4 8 8
√ √ √ !
1 35 3 3 21 3 3 3 3 5 3
= − + + + +
2 64 64 4 32 32 4 32 32
1 9 3
= + +
4 16 16
= 1
The transpose of A is
√
1 3 3
2 4 −
√4
AT
= 3 5 3
−
√4 √8 8
3 3 7
4 8 8
Next AAT is
√ √
1 3 3
2 − 43 1 3
−
√4 2 4 √4
AAT
= 3 5 3 −3 5 3
4√ √8 8 √4 √8 8
3 3 7 3 3 7
− 4 8 8 4 8 8
√ √ √
1 9 3 3 15 3 3 3 3 7 3
4 + 16 + 16 8 − 32 + 32 −8
√
− 32
√
+ 32
√
= 3 15 3 9 25 3 3 3 5 3 7 3
− 32√+ 32 √ 16 + 64√+ 64 √ − 16 + 64 + 64
√8 √
− 83 − 3 3 7 3
32 + 32 − 3163 + 5643 + 7643 3 3 49
16 + 64 + 64
1 0 0
=
0 1 0
0 0 1
= I
As both conditions are satisfied, the matrix A is rotational matrix.
cos θ − sin θ 0
R =
sin θ cos θ 0
(10.4.1)
0 0 1
cos θ sin θ 0
R =
− sin θ cos θ 0
0 0 1
274 10 Rotation and Rigid Bodies
0 sin θ cos θ
Solution Here the axis of rotation is x3 axis, to find rotational matrix we can use
(10.4.1) with θ = π6
cos( π6 ) − sin( π6 ) 0
R = π π
sin( 6 ) cos( 6 ) 0
0 0 1
√
3
2 − 12 0
√
= 1 3
2 2 0
0 0 1
P (X1 , X2 , X3 ) = P (x1 , x2 , x3 )
Solving (10.5.4), we have the vector (10.5.2), known as the axis of rotation.
Let P = (x, y, z), be a point on the axis of rotation, then its position vector relative to
origin is: ~x = (x, y, z)T then the vector equation (10.5.4), becomes
√
1 3
x 2 − 34 x
√4
y = 3 5 3 y
4√ √8 8
3 3 7
z − 4 8 8 z
x = 0
√
3y − 3z = 0
The above system is in reduced Echelon form with x = 0 and z as free variable. Set z = 1,
we have y = √13 .
T
Hence the axis of rotation is x = 0, √13 , 1 and the unit vector along the axis of rotation
D √ E
is x̂ = 0, 12 , 23 .
Normal vector indicates that axis of rotation lies in yz plane. More information about its
direction can be calculated as
√
z 3/2 √
tan θ = = = 3
y 1/2
and
π
θ =
3
Hence axis of rotation passing through origin lies in the yz− plane and is inclined at an
angle of θ = π3 to the positive y − axis.
10.6 Euler’s Angles (φ, θ and ψ) 277
x~0 = Rφ ~x (10.6.1)
0
The next rotation is counterclockwise through an angle θ radian about x1 axis. This rotation
0 0
takes place in x2 x3 plane,then the rotational matrix is
1 0 0
Rθ =
0 cos θ − sin θ
0 sin θ cos θ
0 0 0 00 00 00
After second rotation the coordinate system Ox1 x2 x3 transform to Ox1 x2 x3 and the vector
278 10 Rotation and Rigid Bodies
equation is
R = Rψ Rθ Rφ (10.6.5)
10.7 Rotation about arbitrary axis in 3 − space 279
cos ψ − sin ψ 0 1 0 0 cos φ − sin φ 0
R =
sin ψ cos ψ 0
0 cos θ − sin θ sin φ
cos φ 0
0 0 1 0 sin θ cos θ 0 0 1
cos ψ cos φ − cos θ sin φ sin ψ − cos ψ sin φ − cos θ cos φ sin ψ sin θ sin ψ
=
sin ψ cos φ − cos θ cos ψ sin φ − sin φ sin ψ + cos θ cos φ cos ψ − sin θ cos ψ
sin θ sin φ sin θ cos φ cos θ
000 000
The lines common to the planes x1 x2 plane and x1 x2 planes is called the line of nodes.
Let a rigid body rotates about an arbitrary axis OC passing through origin O. Let n̂ be the unit
vector in this direction. Let P (initial position) be a fixed point in the body having position vector
~r. Let it makes angle φ with OC axis.
~
OP = ~r (10.7.1)
280 10 Rotation and Rigid Bodies
Let the point P has anticlockwise rotation of angle θ, then Q be its new (final) position with position
vector ~r1 and is shown in Fig. 10.12 (a)
~
OQ = ~r1 (10.7.2)
By vector addition, the position vector of this new position of the particle is
~ = ~r1
OQ = ~ + N~Q
ON (10.7.3)
~ is along the axis of rotation. To calculate N~Q, consider right angle triangle N M Q as shown in
ON
Fig. 10.13 (a)
ON = r cos φ = 1 r cos φ
= k n̂ kk ~r k cos φ
= n̂ · ~r
~ can be written as
ON
~
ON = (n̂ · ~r) n̂ (10.7.5)
For N~M and M~Q, consider right angle triangle N M Q in the circle, we have
NM = N Q cos θ
MQ = N Q sin θ
282 10 Rotation and Rigid Bodies
Next N Q = N P (radius of the circle with centre at N ) hence above equations can be written as
NM = N P cos θ (10.7.6)
MQ = N P sin θ (10.7.7)
N P = OP sin φ = r sin φ
= 1 r sin φ
= k n̂ kk ~r k sin φ
= k n̂ × ~r k (10.7.9)
MQ = k n̂ × ~r k sin θ (10.7.10)
It can be viewed from the Fig. (10.17), M Q is in the direction of n̂ × ~r. Let ê be the unit vector in
this direction.
M~Q = k n̂ × ~r k ê sin θ
~
OP = ~ + N~P
ON
then N~P is
N~P = ~ − ON
OP ~ (10.7.12)
~ = ~r1
OQ = (n̂ · ~r) n̂ + [~r − (n̂ · ~r) n̂] cos θ + (n̂ × ~r) sin θ
In this case rotational matrix can be calculated by (10.7.15) into vector equation
r~1 = R~r
where r~1 is vector after rotation, ~r is vector before rotation and R is rotational matrix. Let
~r = hx, y, zi
We will convert all terms on right side of (10.7.15) in matrix form. First consider the term ~r cos θ
a
(1 − cos θ) (n̂ · ~r) n̂ = (1 − cos θ) (ax + by + cz) b
c
2
a x + aby + caz
= (1 − cos θ) abx + b y + bcz
2
acx + bcy + c2 z
2
a ab ac x
= (1 − cos θ) ab b
2
bc y
ac bc c2 x
2
a (1 − cos θ) ab (1 − cos θ) ac (1 − cos θ)
= ab (1 − cos θ) b2 (1 − cos θ) bc (1 − cos θ) ~r
ac (1 − cos θ) bc (1 − cos θ) c2 (1 − cos θ)
î ĵ k̂
(n̂ × ~r) sin θ = a b c sin θ
x y z
Here
2
a (1 − cos θ) + cos θ ab (1 − cos θ) − c sin θ ac (1 − cos θ) + b sin θ
R =
ab (1 − cos θ) + c sin θ b2 (1 − cos θ) + cos θ bc (1 − cos θ) − a sin θ
ac (1 − cos θ) − b sin θ bc (1 − cos θ) + a sin θ c2 (1 − cos θ) + cos θ
(10.7.18)
Example 10.7.1. Find rotational matrix if a rotation of angle θ radian is made about h0, 0, 1i axis
in counter clockwise direction.
Solution Here the axis of rotation is x3 axis with a = 0, b = 0 and c = 1. To find rotational
matrix we can use (10.7.18)
(0)2 (1 − cos θ) + cos θ (0)(0) (1 − cos θ) − (1) sin θ (0)(1) (1 − cos θ) + (0) sin θ
R =
(0)(0) (1 − cos θ) + (1) sin θ (0)2 (1 − cos θ) + cos θ (0)(1) (1 − cos θ) − (0) sin θ
2
(0)(1) (1 − cos θ) − (0) sin θ (0)(1) (1 − cos θ) + (0) sin θ (1) (1 − cos θ) + cos θ
10.7 Rotation about arbitrary axis in 3 − space 287
then we get
cos θ − sin θ 0
R = sin θ
cos θ 0
0 0 1
π
Example 10.7.2. A point (2, 1, 1) has rotation of angle θ = 3 about an axis passing through origin
π
lying in the yz− plane and is inclined at an angle of 3 to the positive y − axis. Determine:
Solution
Consider cartesian coordinate system with O as origin. The position vector of point P (2, 1, 1) is
~r = h2, 1, 1i
π
Let OC be the axis of rotation making an angle 3 to the positive y − axis. Accordingly the axis
π
of rotation lies in yz plane with inclination 3 with y − axis. All this information is shown in Fig.
10.18 The unit vector along axis of rotation is
D π π E
n̂ = 0, cos , sin
3 3
* √ +
1 3
= 0, ,
2 2
π
θ =
3
π
~r cos θ = h2, 1, 1i cos
3
1 1 1
= h2, 1, 1i = 1, ,
2 2 2
√ !* √ +
1 1+ 3 1 3
(1 − cos θ) (n̂ · ~r) n̂ = 0, ,
2 2 2 2
* √ √ +
1+ 3 3+ 3
= 0, ,
8 8
10.7 Rotation about arbitrary axis in 3 − space 289
As
290 10 Rotation and Rigid Bodies
a) RRT = RT R = I or R−1 = RT
b) detR = |R| = 1
Next we study Euler’s theorem, was given by Euler in 1776. It has a fundamental importance in the
theory of rotations of a rigid body.
Theorem 10.8.1. The general displacement of a rigid body with one point fixed is a rotation about
some axis.
Proof : Let us take a regular trihedral OX1 X2 X3 , fixed in space. Introduce another regular
trihedral Ox1 x2 x3 , fixed in the rigid body. Initially both systems coincides with O as same origin.
The position of the body is completely determined by giving the coordinates P (X1 , X2 , X3 ) with
respect to origin O relative to the axes OX1 X2 X3 , while P (x1 , x2 , x3 ) with respect to origin O
relative to the axes Ox1 x2 x3 . If the body undergoes some displacement, then the representation of
~x ~
= RX (10.8.1)
10.8 Euler’s Theorem 291
~x ~
= X (10.8.2)
~x = R~x (10.8.3)
Hence Euler’s theorem will be true if it can be shown that thee exist a vector ~x having the same
components in both the systems. Consider the rotational matrix with property |R| = 1, then for
the eigen vector ~x having eigen value λ, the vector equation is:
R~x = λ~x
Equations (10.8.3) and (10.8.4), with orthogonal matrix R specify the physical motion of a rigid
body with one point fixed always has the eigen value λ = 1.
Alternate Approach: Then the characteristic equation with I3×3 unit matrix
|R − λI| = 0
292 10 Rotation and Rigid Bodies
|R − I| = 0
0
This can be shown as follows: Let R is the transpose of R. Consider
0 0 0
R (R − I) = R R−R
0
Since R is orthogonal matrix, then R R = I, so we can write
0 0
R (R − I) = I −R
0
= −(R − I)
0
= −(R − I)
0
Since R is rotational matrix, so we have |R| = |R | = 1
0
also (R − I) = |R − I|, then above relation may be written as:
0 0
R |(R − I)| = (−1)3 (R − I)
0
|(R − I)| = (−1) (R − I)
|(R − I)| = 0
~
AX ~
= λX (10.8.4)
~
X ~
= RX
10.8 Euler’s Theorem 293
so we have:
1 l11
x l12 l13 x1
x = l l22 l23
x2
2 21
x2 l31 l32 l33 x3
The above system is same as given by (10.5.4) to calculate the axis of rotation.
Example 10.8.1. Determine whether the following matrix is rotational matrix. If yes, find angle
of rotation and axis of rotation.
2
3 − 13 2
3
2 2
R=
3 3 − 31
− 13 2
3
2
3
a) |R| = 1
b) RRT = RT R = I or R−1 = RT
a) The determinant of R is
2
− 13 2
3 3
|R| = 2
3
2
3 − 13
1 2 2
−3
3 3
294 10 Rotation and Rigid Bodies
1
= [2(4 + 2) + (4 − 1) + 2(4 + 2)]
27
= 1
b) The transpose of R is
2 2
3 3 − 31
T
R = 1 2 2
−3 3 3
2
3 − 13 2
3
Next RRT is
2
3 − 13 2
3
2
3
2
3 − 13
T
2 2
RR =
3 3 − 13 −1
3
2
3
2
3
− 13 2
3
2
3
2
3 − 13 2
3
4 4 1
+ + 9 9 9 − 92 + 4
9 − 2
9
4
9 − 2
9 − 2
9
2 4 2 1 4 4
=
−9 + 9 − 9 9 + 9 + 9 − 29 − 2
9 + 4
9
4 2 2
9 − 9 − 9 − 92 − 2
9 + 4
9
4
9 + 1
9 + 4
9
1 0 0
=
0 1 0
0 0 1
= I
2 2 2
T r(R) = + + =2
3 3 3
10.8 Euler’s Theorem 295
Axis of Rotation: Let P = (x1 , x2 , x3 ), be the position of the particle on the axis of rotation,
T
Then its position vector is: ~x = (x1 , x2 , x3 ) then the vector equation (10.8.3), becomes
2
x1 3 − 13 2
3 x1
x = 2 2
− 13 x
2 3 3 2
x2 − 13 2
3
2
3 x3
2 1 2
x1 = x1 − x2 + x3
3 3 3
2 2 1
x2 = x1 + x2 − x3
3 3 3
1 2 2
x3 = − x1 + x2 + x3
3 3 3
x1 + x2 − 2x3 = 0
2x1 − x2 − x3 = 0
x1 − 2x2 + x3 = 0
x1 + x2 − 2x3 = 0
−x2 + x3 = 0
The above system is in reduced Echelon form with x3 as free variable. Set x3 = 1, we have x2 = 1
and x1 = 1.
296 10 Rotation and Rigid Bodies
T
Hence the axis of rotation is x = (1, 1, 1) and the unit vector along the axis of rotation is
1 1 1
x̂ = √ ,√ ,√
3 3 3
= hcos α, cos β, cos γi
Where α, β and γ are angles that the axis of rotation makes with coordinate axes. Hence
1
cos α = √
3
and
1
α = arccos √ =β=γ
3
= 0.955 radian = 54.7◦
The axis of rotation makes an angle of 54.7◦ with all the coordinate axes.
Consider two coordinate systems with common origan O. Let OX1 X2 X3 be a fixed system and
let Ox1 x2 x3 is rotating with angular velocity ω about some arbitrary axis passing through origin.
Initially both systems coincides with O as same origin. Let iˆ0 , jˆ0 , kˆ0 be unit vectors in fixed system
while î, ĵ, k̂ are unit vectors in rotating system The position of the body is completely determined by
giving the coordinates P (X1 , X2 , X3 ) with respect to origin O relative to the axis OX1 X2 X3 , while
P (x1 , x2 , x3 ) with respect to origin O relative to the axis Ox1 x2 x3 . Let r be the position vector of
point P then
and
~r = x1 î + x2 ĵ + x3 k̂ (10.9.2)
10.9 Dynamics of Rotating Coordinate System 297
The transformation equations from moving system to fixed system can be obtained by taking the
dot product of (10.9.1) and (10.9.2) with iˆ0 , jˆ0 , kˆ0 respectively
~r · iˆ0 = X1
~r · jˆ0 = X2 (10.9.3)
~r · kˆ0 = X3
(The equations for inverse transformation can be obtained by taking dot product of (10.9.1) and
(10.9.2) with î, ĵ, k̂)
298 10 Rotation and Rigid Bodies
The dot products on right hand side of (10.9.5) are simply the cosines of the angles between the
corresponding axes. If the rotation of angle θ is made about x3 axis in anticlockwise direction, the
system will transform as
X3 = (0)x1 + (0)x2 + x3
The linear velocity ~v of a particle having position vector ~r and rotating with angular velocity ω
~
d~r
~v = ~ × ~r
=ω (10.9.8)
dt
10.9 Dynamics of Rotating Coordinate System 299
d~
r d~
r
Here dt f ix is the velocity of fixed system relative to fixed system and dt rot is velocity of rotating
system relative to fixed system. The time derivative of ~r, however will be different in two systems.
In fixed system
d~r
= Ẋ1 iˆ0 + Ẋ2 jˆ0 + Ẋ3 kˆ0 (10.9.10)
dt f ix
However, in rotating system the unit vectors are changing in directions, hence their time derivatives
will appear with respect to fixed axis.
d~r dî dĵ dk̂
= ẋ1 î + ẋ2 ĵ + ẋ3 k̂ + x1 + x2 + x3 (10.9.11)
dt f ix dt dt dt
In rotating frame the unit vectors î, ĵ, k̂ are treated as constant unit vectors, then the velocity in
the rotating system is
d~r
= ẋ1 î + ẋ2 ĵ + ẋ3 k̂ (10.9.12)
dt rot
Since î, ĵ, k̂ are the unit vectors in a system rotating with angular velocity ω, then (10.9.8) can be
applied as a special case as
dî
= ~ × î
ω
dt
dĵ
= ~ × ĵ
ω (10.9.13)
dt
dk̂
= ~ × k̂
ω
dt
(10.9.14) can be treated as an operator equation which gives the relation between the time derivative
in the fixed and the rotating coordinate systems.
d d
= ~×
+ω (10.9.15)
dt f ix dt rot
300 10 Rotation and Rigid Bodies
d~
r ∂~
r
The operator equation (10.9.15) can be operated on any vector. If we denote dt rot = ∂t (derivative
relative to moving frame) in (10.9.14), then we have
d~r ∂~r
~v = = ~ × ~r
+ω (10.9.16)
dt ∂t
d ∂
= ~×
+ω (10.9.17)
dt ∂t
dF~ ∂ F~
= ~ × F~
+ω (10.9.18)
dt ∂t
10.9.1 Acceleration
d~v ∂~v
~a = = ~ × ~v
+ω (10.9.19)
dt ∂t
Corollary 10.9.1. : The angular acceleration is the same in fixed and rotating systems.
Proof : Consider the operator equation (10.9.15) with angular velocity vector ω
~.
d~
ω d~
ω
= +ω~ ×ω~
dt f ix dt rot
~ ×ω
Since ω ~ = 0, so we have
d~
ω d~
ω
= (10.9.21)
dt f ix dt rot
Hence the angular acceleration is the same in fixed and rotating systems.
10.9 Dynamics of Rotating Coordinate System 301
Example 10.9.1. The points (a, 2a, −a), (−a, −a, a) and (a, a, a) of a rigid body have instantaneous
D√ √ E D E D E
velocities 23 v, 0, 23 v , − √13 v, 0, − √13 v and 0, − √13 v, √13 v . Show that the body has the line
D E
through the origin having direction cosines √13 , − √13 , − √13 as instantaneous axis of rotation and
1
that the magnitude of the angular velocity is 2a v.
Solution:
Set O as origin in the body. Let A = (a, 2a, −a), B = (−a, −a, a) and C = (a, a, a) be the points of
~rC = ha, a, ai
D√ √ E D E D E
3 3
Also let ~vA = 2 v, 0, 2 v , ~vB = − √13 v, 0, − √13 v and ~vC = 0, − √13 v, √13 v be the linear
velocities of A, B and C respectively. Let OP be the instantaneous axis of rotation with angular
velocity ω, then
ω
~ = hω1 , ω2 , ω3 i
~vA = ~ × ~rA
ω
302 10 Rotation and Rigid Bodies
î ĵ k̂
= ω1 ω2 ω3
a 2a −a
But
*√ √ +
3 3
vA = v, 0, v (10.9.23)
2 2
ω2 = −ω1 (10.9.27)
~vB = ~ × ~rB
ω
î ĵ k̂
= ω1 ω2 ω3
−a −a a
But
1 1
vB = − √ v, 0, − √ v (10.9.29)
3 3
10.9 Dynamics of Rotating Coordinate System 303
1
aω2 + aω3 = −√ v (10.9.30)
3
−aω1 − aω3 = 0 (10.9.31)
1
−aω1 + aω2 = −√ v (10.9.32)
3
ω3 = −ω1 (10.9.33)
1 1 v
ω1 = √ v=√ (10.9.34)
2a 3 3 2a
1 v
ω2 = −√ (10.9.35)
3 2a
1 v
ω3 = −√ (10.9.36)
3 2a
Hence
1 v 1 v 1 v
ω
~ = √ , −√ , −√ (10.9.37)
3 2a 3 2a 3 2a
v 1 1 1
= √ , −√ , −√ (10.9.38)
2a 3 3 3
Since
1
√ , − √1 , − √1
=
3 1
3 3
so we can write
v
ω = (10.9.39)
2a
Polar coordinates can be extended to three dimensions in a very straightforward manner by adding
the z coordinate. Every point in space is determined by the r and θ coordinates of its projection in
the xy plane, and its z coordinates. Consider Oxyz cartesian coordinate system. Let P (x, y, z) be
a point and |P B| is perpendicular from P on xOy plane. Then
∠AOB = θ
~
|OB| = r
|P~B| = z
Then (r, θ, z) are called the cylindrical polar coordinates of P . Let P varies with time, and
~ = ~a = OB
OP ~ + BP
~
= rr̂ + z ẑ
Where r̂, θ̂, ẑ are unit vectors in the directions of r, θ, z. They constitute a right hand frame for
which
r̂ × θ̂ = ẑ
θ̂ × ẑ = r̂
ẑ × r̂ = θ̂
Since P is in a frame specified by these unit vectors and moving with angular velocity ω along Oz
axis
ω = θ̇ẑ
The velocity of P can be given by using operator equation (10.9.15) with position vector ~a.
d~a ∂~a
~v = = ~ × ~a
+ω (10.10.3)
dt ∂t
∂
~v = rr̂ + 0θ̂ + z ẑ + 0r̂ + 0θ̂ + θ̇ẑ × rr̂ + 0θ̂ + z ẑ
∂t
= ṙr̂ + 0θ̂ + ż ẑ + 0r̂ + rθ̇θ̂ + 0ẑ
= ṙr̂ + rθ̇θ̂ + ż ẑ (10.10.4)
The acceleration of P can be calculated by using operator equation (10.9.15) with velocity vector ~v .
∂
~a = ṙr̂ + rθ̇θ̂ + ż ẑ + 0r̂ + 0θ̂ + θ̇ẑ × ṙr̂ + rθ̇θ̂ + ż ẑ
h∂t i
= r̈r̂ + rθ̈ + ṙθ̇ θ̂ + z̈ ẑ + −rθ̇2 r̂ + ṙθ̇θ̂ + 0ẑ
h i
= r̈ − rθ̇2 r̂ + rθ̈ + 2ṙθ̇ θ̂ + z̈ ẑ (10.10.5)
In spherical coordinates, we utilize two angles and a distance to specify the position of a particle.
Consider Oxyz cartesian coordinate system. Let P (x, y, z) be a point. The circle in the P OZ plane
306 10 Rotation and Rigid Bodies
∠AOB = φ
∠P OZ = θ
OP = r
Then (r, θ, φ) are called the spherical polar coordinates of P . Let P varies with time, and
~
OP = ~r
Where r̂, θ̂, φ̂ are unit vectors in the directions of r, θ, φ. They constitute a right hand frame for
which
r̂ × θ̂ = φ̂
θ̂ × φ̂ = r̂
φ̂ × r̂ = θ̂
In this motion two angles are involved, hence the frame has angular velocity as
The units vectors r̂, θ̂ and k̂ lies in the same plane constituted by r̂ and k̂. The unit vectors r̂ and
π
θ̂ are at right angles and k̂ makes angles θ with r̂ and 2 + θ with θ̂ as shown in Fig. 10.27. Then k̂
can be written as
π
k̂ = cos θr̂ + cos + θ θ̂
2
= cos θr̂ − sin θθ̂ (10.11.3)
The velocity of P can be given by using operator equation (10.9.15) with position vector ~r given by
308 10 Rotation and Rigid Bodies
(10.33).
d~r ∂~r
~v = = ~ × ~r
+ω (10.11.5)
dt ∂t
∂
~v = rr̂ + 0θ̂ + 0φ̂ + φ̇ cos θr̂ − φ̇ sin θθ̂ + θ̇φ̂ × rr̂ + 0θ̂ + 0ẑ
∂t
= ṙr̂ + 0θ̂ + 0φ̂ + 0r̂ + rθ̇θ̂ + rφ̇ sin θφ̂
= ṙr̂ + rθ̇θ̂ + rφ̇ sin θφ̂ (10.11.6)
The acceleration of P can be calculated by using operator equation (10.9.15) with velocity vector ~v .
∂
~a = ṙr̂ + rθ̇θ̂ + rφ̇ sin θφ̂ + φ̇ cos θr̂ − φ̇ sin θθ̂ + θ̇φ̂ × ṙr̂ + rθ̇θ̂ + rφ̇ sin θφ̂
h∂t i
= r̈r̂ + rθ̈ + ṙθ̇ θ̂ + rφ̈ sin θ + ṙφ̇ sin θ + rθ̇φ̇ cos θ φ̂
h i
+ −rφ̇2 sin2 θ − rθ̇2 r̂ + ṙθ̇ − rφ̇2 sin θ cos θ θ̂ + rθ̇φ̇ cos θ + ṙφ̇ sin θ φ̂
= r̈ − rθ̇2 − rφ̇2 sin2 θ r̂ + rθ̈ + 2ṙθ̇ − rφ̇2 sin θ cos θ θ̂
+ 2rθ̇φ̇ cos θ + rφ̈ sin θ + 2ṙφ̇ sin θ φ̂ (10.11.7)
10.12 Screw Motion 309
A motion combination of rotation and translation motion is termed as screw motion as illustrated
in Fig. 10.28.
Screw motion consist a of pairs of vectors, such as forces and moments and angular and linear
velocity, that arise in the kinematics and dynamics of rigid bodies. The mathematical framework was
developed by Sir Robert Stawell Ball in 1876 for application in kinematics and statics of mechanisms
(rigid body mechanics).
A spatial displacement of a rigid body can be defined by a rotation about a line and a translation
along the same line, called a screw displacement. This is known as Chasles’ theorem.
The most general rigid body displacement can be produced by a translation along a line followed
(or preceded) by a rotation about that line.
To establish this result, discovered by Chasles in 1830, we note that the positions in space of any
310 10 Rotation and Rigid Bodies
three non-collinear points of a rigid body determine the position of the rigid body. Let A, B, C be
any three non-collinear points of the body which are displaced to other positions A1 , B2 , C2 . Thus
the points A, B, C determine the initial position of the body and A1 , B2 , C2 determine the final
position. This screw displacement takes place in two steps. In first step, the motion is translation,
in which A moves A1 , B moves B1 and C moves C1 as shown in Fig. 10.29. In second step, the
motion is rotation about A1 in which B1 moves B2 and C1 moves C2 . The point A is called base
point. We can choose any other point as base point. Then the translation will be altered but rotation
will be independent of the choice of base point.
Let a rigid body rotates about a fixed point O. Let ω be its angular velocity about the instantaneous
axis OC at time t. Let P be a fixed point in the body such that
~
OP = ~r (10.14.1)
10.14 Fundamental Properties of Screw Motion 311
making an angle θ with OC axis. Draw a perpendicular P M from P on OC. Then in triangle OP M
PM = r sin θ (10.14.2)
v = ωr (10.14.3)
Where r is the radius of the circle. From Fig 10.30 we can see that the circular path has radius
r sin θ, so the speed of P is
v = ωr sin θ (10.14.4)
~v = ~ × ~r
ω (10.14.5)
If in addition, the body as a whole has a constant translation velocity ~v0 , then the total linear
velocity of the particle at P is
~v ~ × ~r
= ~v0 + ω (10.14.6)
Corollary 10.14.1. For a rotating rigid body about a fixed point, prove that
curl(~v ) = 2~
ω
312 10 Rotation and Rigid Bodies
Figure 10.31: Rigid body with OC as axis of rotation having translation motion
Solution Consider
ω
~ = hω1 , ω2 , ω3 i
and
~r = hx1 , x2 , x3 i
î ĵ k̂
~ × ~r = ω1
~v = ω ω2 ω3
x1 x2 x3
= hx3 ω2 − x2 ω3 , x1 ω3 − x3 ω1 , x2 ω1 − x1 ω2 i (10.14.7)
Next
curl(~v ) = ∇ × ~v (10.14.8)
Where
∂ ∂ ∂
∇ = , , (10.14.9)
∂x1 ∂x2 ∂x3
10.14 Fundamental Properties of Screw Motion 313
î ĵ k̂
curl (~v ) = ∂ ∂ ∂
∂x1 ∂x2 ∂x3
x3 ω2 − x2 ω3 x1 ω3 − x3 ω1 x2 ω1 − x1 ω2
∂ ∂
= (x2 ω1 − x1 ω2 ) − (x1 ω3 − x3 ω1 ) î
∂x2 ∂x3
∂ ∂
+ (x3 ω2 − x2 ω1 ) − (x2 ω1 − x1 ω2 ) ĵ
∂x3 ∂x1
∂ ∂
+ (x1 ω3 − x3 ω1 ) − (x3 ω2 − x2 ω3 ) k̂
∂x1 ∂x2
1 i=j
∂xi
∂xj =
0 i 6= j
so we have
= 2 hω1 , ω2 , ω3 i
= 2~
ω
In screw motion, the axis of rotation lie in the direction of the translation motion i.e. the linear
velocity v and angular velocity ω has the same direction. Consider a rigid body. Let O be a fixed
point in the rigid body, having velocity v and ω be the instantaneous angular velocity of the body.
314 10 Rotation and Rigid Bodies
Let OC be the axis of rotation. Let P be another point of the body such that
~
OP = ~r
v~1 ~ × ~r
= ~v + ω (10.14.10)
~ × v~1
ω = 0 (10.14.11)
0 = ~ × ~v + ω
ω ~ × (~
ω × ~r)
= ~ × ~v + (~
ω ~ − ω 2~r
ω · ~r) ω
or we can write
ω 2~r = ~ × ~v + (~
ω ω · ~r) ω
~ (10.14.12)
~ × ~v
ω
~r = + λ~
ω (10.14.13)
ω2
ω ·~
(~ r)
where λ = ω2 As λ varies, (10.14.13) represents the vector equation of straight line passing
~ ×~
ω v
through the point represented by ~r0 = ω2 and having the same direction as ω
~ . At the particular
10.14 Fundamental Properties of Screw Motion 315
instant considered every point on the line (10.14.13), having the velocity parallel to ω
~ . Hence the
instantaneous motion of the body is a screw motion about the line. This line is called external axis
or axis of screw. Every point on the axis moves along it and the body turns about the axis.
Since the velocity of P is v~1 given by (10.14.10). Taking its dot product with ω
~
~ · v~1
ω = ~ · (~v + ω
ω ~ × ~r)
= ~ · ~v + ω
ω ~ · (~
ω × ~r)
= ~ · ~v + 0
ω
and we have
~ · v~1
ω = ~ · ~v
ω (10.14.14)
~ ·v~1
ω
Next consider the quantity ω2 , with the velocity v~1 is parallel to ω
~ . Then
~ · v~1
ω ωv1 cos(0)
=
ω2 ω2
v1
=
ω
distance per second
=
angle per second
or
~ · v~1
ω distance
= = distance per unit angle (10.14.15)
ω2 angle
316 10 Rotation and Rigid Bodies
The ratio (10.14.15) is called pitch of the screw, i.e. distance moved by the particle per unit angle.
The rotations and other rotational quantities derived from it may also be treated as vector quanti-
ties. The addition of such two vector quantities, should however commutative, but the reality is not
so. Consider two regular trihedral systems with common origin O, one fixed and the other is rotat-
π
able. Let OX1 X2 X3 , fixed in space and Ox1 x2 x3 is rotatable. First we made a rotation of θ1 = 2
about X3 − axis. In this rotation, x3 will remain unchanged, while x1 and x2 will interchange their
π
positions . Considering the new positions, we made an other rotation of θ2 = 2 about X2 − axis.
In this rotation, x1 (along X2 − axis) will remain unchanged, while x2 and x3 will interchange their
positions (see Fig. 10.34 left side). In this way we complete the operation θ1 + θ2 . Next we reverse
π
the order. First we made a rotation of θ2 = 2 about X2 − axis. In this rotation, x2 will remain
unchanged, while x1 and x3 will interchange their positions. Considering the new positions, we made
π
an other rotation of θ1 = 2 about X3 − axis. In this rotation, x1 (along X3 − axis) will remain
unchanged, while x2 and x3 will interchange their positions (see Fig. 10.34 right side). In this way
we complete the operation θ2 + θ1 . Above Fig. 10.34, shows that addition of two finite rotations
π π
θ2 = 2 & θ2 = 2 is not commutative. i.e., θ1 + θ2 does not yield the position of the body as by
θ2 + θ1
The body which is subjected to these two finite rotations is not found to be in the same final state.
Hence finite rotations can not be regarded as vector quantities.
Consider two regular trihedral systems with common origin O, one fixed and the other is rotatable.
Let OX1 X2 X3 , fixed in space and Ox1 x2 x3 is rotatable. We made a rotation of a very small angle
about X3 − axis. The axes X1 − axis and X2 − axis moves but maintaining the same direction.
Here θ ≈ 0 so sin θ ≈ θ. Then we shall see that, for very small rotations ∆θ1 &∆θ2 , the addition will
be commutative. To represent infinitesimal rotation as a vector, we draw a straight line along the
axis of rotation with the properties.
10.14 Fundamental Properties of Screw Motion 317
~
OP = ~r (10.14.16)
Let the rigid body turns through a small angle ∆θ1 in the positive sense about an axis through O
specified by unit vector â1 , then
ω
~ = ∆θ1 â1 (10.14.17)
318 10 Rotation and Rigid Bodies
Now suppose the body rotates through another small angle ∆θ2 in the opposite sense about an axis
through O specified by a unit vector â2 . and the movement of P is
Since ∆θ1 and ∆θ2 are very small, consider only first order term, we can write
If the operations were performed in the reverse order, the position of P relative to O, considering
only first order term,is
If ∆t be the time interval in which these operations take place, the velocity of P relative to O is
~r12 − ~r ∆θ1 ∆θ2
lim = lim â1 + â2 × ~r
∆t→0 ∆t ∆t→0 ∆t ∆t
= θ̇1 â1 + θ̇2 â2 × ~r
= (~ ~ 2 ) × ~r
ω1 + ω
where ω
~ 1 = θ̇1 â1 vector angular velocity about â1
and ω
~ 2 = θ̇2 â2 vector angular velocity about â2
This result confirms that vector quantities about a point of a rigid body are commutative with
respect to addition, as the other quantities are.
320 10 Rotation and Rigid Bodies
Exercises
1. consider the point (x, y) = (1, 1) in a frame rotated by 45o counter-clockwise. Find its new
position.
(a) Find the xy−coordinates of the point whose x́ý−coordinates are (−2, 6).
√
(b) Find an equation of the curve 3x́ý + ý 2 = 6 in xy−coordinates.
√
(a) Find the xy−coordinates of the point whose x́ý−coordinates are 1, − 3 .
√
(b) Find an equation of the curve 2x́2 + 2 3x́ý = 3 in xy−coordinates.
4. Find the angle and axis of rotation corresponding to the following rotational matrices.
5
2 − 12 3
4
2 1
(a)
3 2 − 34
− 23 4
5 − 52
1 2
5 5 − 15
1 3
(b)
−5 5 − 25
2
5 − 15 − 45
10.14 Fundamental Properties of Screw Motion 321
1 3
4 4 − 14
3
(c)
4 − 14 1
4
− 14 3
4 − 34
5. The points (a, 0, 0), 0, √a3 , 0 and (0, 0, 2a) of a rigid body have instantaneous velocities
√
hu, 0, 0i, hu, 0, vi and u + v, 3v, v2 respectively, referred to a regular trihedral system. Find
the magnitude and direction of spin of the body and the point at which the certain axis cuts
the xz-plane.
322 BIBLIOGRAPHY
Bibliography
[1] Q.K.Ghori, Introduction to Mechanics, West Pak Publishing Co. Ltd., Lahore 1971.
[2] D. Halliday, R. Resnick and J. Walker, Fundamentals of Physics,John Wiley & Sons, lnc., 8th
ed. extended, 2008.
[3] D. G. Zill and M. R. Cullen, Differential Equations with Boundary-Value Problems, 7th Edi-
tion, 2009.
[4] Richard Bronson and G. B. Costa, Schaum’s Outline of Differential Equations, McGraw-Hill
Companies, Inc., 3rd Edition, 2006.
[5] H. Anton, I. Bivens and S. Davis, Calculus, John Wiley & Sons, lnc., Early Transcendentals
10th edition, 2012.
[6] L. D. Landau, E. M. Lifshitz, Mecanics, Pergamon Press, Oxford: butterworth, 3rd edition,
15 2005.
[7] B. Ahmad, S. Borisenok, Control of effective potential minima for Kapitza oscillator by peri-
odical kicking pulses, Physics Letters A, 373, 701-707, 2009.
[8] B. Ahmad, Stabilization of Kapitza Oscillator by Periodic Symmetric Pulses, Nonlinear Dy-
namics, 62, 499-506, 2010.
[9] K. Singh, Mechanics, Department of Mathematics, Directorate of Distance Education, Guru
Jambheshwar University of Science & Technology Hisar-125001.
[10] S.T. Thornton and J.B. Marion, Classical Dynamics of Particles and Systems, Cengage Learn-
ing India Private Ltd. Third Indian Reprint, 2009.
[11] S.S. Bhavikatti, A textbook of Classical Mechanics, New Age International Publishers, Second
Edition Reprint, 2009.
[12] R.G. Takwale and P.S. Puranik, Introduction to Classical Mechanics, Tata McGraw-Hill Pub-
lishing Company Ltd. 39th reprint 2009.
[13] Lecture notes of Mechanics MIT, USA.
[14] Wikipedia.
Index
323
324 INDEX
time-acceleration curve, 15
time-displacement curve, 12
time-velocity curve, 13
trajectory of a projectile, 154, 177
translation, 6
transverse component, 208
two body problem, 204