S. Widnall, J.

Peraire
16.07 Dynamics
Fall 2009
Version 2.0

Lecture L4 - Curvilinear Motion. Cartesian Coordinates

We will start by studying the motion of a particle. We think of a particle as a body which has mass,
but has negligible dimensions. Treating bodies as particles is, of course, an idealization which involves an
approximation. This approximation may be perfectly acceptable in some situations and not adequate in
some other cases. For instance, if we want to study the motion of planets, it is common to consider each
planet as a particle. This simplification is not adequate if we wish to study the precession of a gyroscope or
a spinning top.

Kinematics of curvilinear motion

In dynamics we study the motion and the forces that cause, or are generated as a result of, the motion.
Before we can explore these connections we will look first at the description of motion irrespective of the
forces that produce them. This is the domain of kinematics. On the other hand, the connection between
forces and motions is the domain of kinetics and will be the subject of the next lecture.

Position vector and Path

We consider the general situation of a particle moving in a three dimensional space. To locate the position of
a particle in space we need to set up an origin point, O, whose location is known. The position of a particle
A, at time t, can then be described in terms of the position vector, r, joining points O and A. In general,
this particle will not be still, but its position will change in time. Thus, the position vector will be a function
of time, i.e. r(t). The curve in space described by the particle is called the path, or trajectory.

We introduce the path or arc length coordinate, s, which measures the distance traveled by the particle along
the curved path. Note that for the particular case of rectilinear motion (considered in the review notes) the
arc length coordinate and the coordinate, s, are the same.

1
Using the path coordinate we can obtain an alternative representation of the motion of the particle. Consider
that we know r as a function of s, i.e. r(s), and that, in addition we know the value of the path coordinate
as a function of time t, i.e. s(t). We can then calculate the speed at which the particle moves on the path
simply as v = ṡ ≡ ds/dt. We also compute the rate of change of speed as at = s̈ = d2 s/dt2 .
We consider below some motion examples in which the position vector is referred to a fixed cartesian
coordinate system.

Example Motion along a straight line in 2D

Consider for illustration purposes two particles that move along a line defined by a point P and a unit vector
m. We further assume that at t = 0, both particles are at point P . The position vector of the first particle is
given by r 1 (t) = r P + mt = (rP x + mx t)i + (rP y + my t)j, whereas the position vector of the second particle
is given by r 2 (t) = r P + mt2 = (rP x + mx t2 )i + (rP y + my t2 )j.

Clearly the path for these two particles is the same, but the speed at which each particle moves along the
path is different. This is seen clearly if we parameterize the path with the path coordinate, s. That is,
we write r(s) = r P + ms = (rP x + mx s)i + (rP y + my s)j. It is straightforward to verify that s is indeed
the path coordinate i.e. the distance between two points r(s) and r(s + Δs) is equal to Δs. The two
motions introduced earlier simply correspond to two particles moving according to s1 (t) = t and s2 (t) = t2 ,
respectively. Thus, r 1 (t) = r(s1 (t)) and r 2 (t) = r(s2 (t)).

It turns out that, in many situations, we will not have an expression for the path as a function of s. It is
in fact possible to obtain the speed directly from r(t) without the need for an arc length parametrization of
the trajectory.

Velocity Vector

We consider the positions of the particle at two different times t and t + Δt, where Δt is a small increment
of time. Let Δr = r(r + Δt) − r(t), be the displacement vector as shown in the diagram.

2
The average velocity of the particle over this small increment of time is

Δr
v ave = ,
Δt

which is a vector whose direction is that of Δr and whose magnitude is the length of Δr divided by Δt. If
Δt is small, then Δr will become tangent to the path, and the modulus of Δr will be equal to the distance
the particle has moved on the curve Δs.
The instantaneous velocity vector is given by

Δr dr(t)
v = lim ≡ ≡ ṙ , (1)
Δt→0 Δt dt

and is always tangent to the path. The magnitude, or speed, is given by

Δs ds
v = |v | = lim ≡ ≡ s˙ .
Δt→0 Δt dt

Acceleration Vector

In an analogous manner, we can define the acceleration vector. Particle A at time t, occupies position
r(t), and has a velocity v(t), and, at time t + Δt, it has position r(t + Δt) = r(t) + Δr, and velocity
v(t + Δt) = v(t) + Δv. Considering an infinitesimal time increment, we define the acceleration vector as the
derivative of the velocity vector with respect to time,

Δv dv d2 r
a = lim ≡ = 2 . (2)
Δt→0 Δt dt dt

We note that the acceleration vector will reflect the changes of velocity in both magnitude and direction.
The acceleration vector will, in general, not be tangent to the trajectory (in fact it is only tangent when the
velocity vector does not change direction). A sometimes useful way to visualize the acceleration vector is to

3
translate the velocity vectors, at different times, such that they all have a common origin, say, O� . Then,
the heads of the velocity vector will change in time and describe a curve in space called the hodograph. We
then see that the acceleration vector is, in fact, tangent to the hodograph at every point.

Expressions (1) and (2) introduce the concept of derivative of a vector. Because a vector has both magnitude
and direction, the derivative will be non-zero when either of them changes (see the review notes on
vectors). In general, the derivative of a vector will have a component which is parallel to the vector itself,
and is due to the magnitude change; and a component which is orthogonal to it, and is due to the direction
change.

Note Unit tangent and arc-length parametrization

The unit tangent vector to the curve can be simply calculated as

et = v/v.

It is clear that the tangent vector depends solely on the geometry of the trajectory and not on the speed
at which the particle moves along the trajectory. That is, the geometry of the trajectory determines the
tangent vector, and hence the direction of the velocity vector. How fast the particle moves along the
trajectory determines the magnitude of the velocity vector. This is clearly seen if we consider the arc-length
parametrization of the trajectory r(s). Then, applying the chain rule for differentiation, we have that,
dr dr ds
v= = = et v ,
dt ds dt
where, ṡ = v, and we observe that dr/ds = et . The fact that the modulus of dr/ds is always unity indicates
that the distance traveled, along the path, by r(s), (recall that this distance is measured by the coordinate
s), per unit of s is, in fact, unity!. This is not surprising since by definition the distance between two
neighboring points is ds, i.e. |dr| = ds.

Cartesian Coordinates
When working with fixed cartesian coordinates, vector differentiation takes a particularly simple form. Since
the vectors i, j, and k do not change, the derivative of a vector A(t) = Ax (t)i + Ay (t)j + Az (t)k, is simply
Ȧ(t) = A˙x (t)i + A˙y (t)j + A˙z (t)k. That is, the components of the derivative vector are simply the derivatives
of the components.

4
Thus, if we refer the position, velocity, and acceleration vectors to a fixed cartesian coordinate system, we
have,

r(t) = x(t)i + y(t)j + z(t)k (3)

v(t) = vx (t)i + vy (t)j + vz (t)k = ẋ(t)i + ẏ(t)j + ż(t)k = ṙ(t) (4)

a(t) = ax (t)i + ay (t)j + az (t)k = v̇x (t)i + v̇y (t)j + v̇z (t)k = v̇(t) (5)

� �
Here, the speed is given by v = vx2 + vy2 + vz2 , and the magnitude of the acceleration is a = a2x + a2y + a2z .
The advantages of cartesian coordinate systems is that they are simple to use, and that if a is constant, or
a function of time only, we can integrate each component of the acceleration and velocity independently as
shown in the ballistic motion example.

Example Circular Motion

We consider motion of a particle along a circle of radius R at a constant speed v0 . The parametrization of
a circle in terms of the arc length is
s s
r(s) = R cos( )i + R sin( )j .
R R

Since we have a constant speed v0 , we have s = v0 t. Thus,

v0 t v0 t
r(t) = R cos( )i + R sin( )j .
R R
The velocity is
dr(t) v0 t v0 t
v(t) = = −v0 sin( )i + v0 cos( )j ,
dt R R

5
which, clearly, has a constant magnitude |v| = v0 . The acceleration is,
dr(t) v2 v0 t v2 v0 t
a(t) = = − 0 cos( )i − 0 sin( )j .
dt R R R R
Note that, the acceleration is perpendicular to the path (in this case it is parallel to r), since the velocity
vector changes direction, but not magnitude.
We can also verify that, from r(s), the unit tangent vector, et , could be computed directly as
dr(s) s s v0 t v0 t
et = = − sin( )i + cos( ) = − sin( )i + cos( )j .
ds R R R R

Example Motion along a helix

The equation r(t) = R cos ti + R sin tj + htk, defines the motion of a particle moving on a helix of radius R,
and pitch 2πh, at a constant speed. The velocity vector is given by
dr
v= = −R sin ti + R cos tj + hk ,
dt
and the acceleration vector is given by,
dv
a= = −R cos ti + −R sin tj .
dt
In order to determine the speed at which the particle moves we simply compute the modulus of the velocity
vector,
� �
v = |v| = R2 sin2 t + R2 cos2 t + h2 = R 2 + h2 .

If we want to obtain the equation of the path in terms of the arc-length coordinate we simply write,
�
ds = |dr| = vdt = R2 + h2 dt .
√
Integrating, we obtain s = s0 + R2 + h2 t, where s0 corresponds to the path coordinate of the particle
at time zero. Substituting t in terms of s, we obtain the expression for the position vector in terms of the
√ √ √
arc-length coordinate. In this case, r(s) = R cos(s/ R2 + h2 )i + R sin(s/ R2 + h2 )j + hs/ R2 + h2 k. The
figure below shows the particle trajectory for R = 1 and h = 0.1.

1.5

0.5

0
1
0.5 1
0 0.5
0.5 0
0.5
1 1

6
Example Ballistic Motion

Consider the free-flight motion of a projectile which is initially launched with a velocity v 0 = v0 cos φi +
v0 sin φj. If we neglect air resistance, the only force on the projectile is the weight, which causes the projectile
to have a constant acceleration a = −gj. In component form this equation can be written as dvx /dt = 0
and dvy /dt = −g. Integrating and imposing initial conditions, we get

vx = v0 cos φ, vy = v0 sin φ − gt ,

where we note that the horizontal velocity is constant. A further integration yields the trajectory

1
x = x0 + (v0 cos φ) t, y = y0 + (v0 sin φ) t − gt2 ,
2

which we recognize as the equation of a parabola.

The maximum height, ymh , occurs when vy (tmh ) = 0, which gives tmh = (v0 /g) sin φ, or,

v02 sin2 φ
ymh = y0 + .
2g

The range, xr , can be obtained by setting y = y0 , which gives tr = (2v0 /g) sin φ, or,

2v02 sin φ cos φ v 2 sin(2φ)

xr = x0 + = x0 + 0 .
g g

We see that if we want to maximize the range xr , for a given velocity v0 , then sin(2φ) = 1, or φ = 45o .
Finally, we note that if we want to model a more realistic situation and include aerodynamic drag forces of
the form, say, −κv 2 , then we would not be able to solve for x and y independently, and this would make the
problem considerably more complicated (usually requiring numerical integration).

ADDITIONAL READING
J.L. Meriam and L.G. Kraige, Engineering Mechanics, DYNAMICS, 5th Edition
2/1, 2/3, 2/4

7
MIT OpenCourseWare
http://ocw.mit.edu

16.07 Dynamics
Fall 2009

For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.