Full Download Orbital Mechanics For Engineering Students Fourth Edition. Edition Howard D. Curtis PDF
Full Download Orbital Mechanics For Engineering Students Fourth Edition. Edition Howard D. Curtis PDF
com
https://ebookmass.com/product/orbital-mechanics-for-
engineering-students-fourth-edition-edition-howard-d-curtis/
OR CLICK HERE
DOWLOAD NOW
https://ebookmass.com/product/aerodynamics-for-engineering-
students-7th-edition-e-l-houghton/
ebookmass.com
https://ebookmass.com/product/statistical-mechanics-fourth-edition-r-
k-pathria/
ebookmass.com
https://ebookmass.com/product/chemistry-for-engineering-students-4th-
edition-lawrence-s-brown/
ebookmass.com
https://ebookmass.com/product/sourcing-and-procurement-with-
sap-s-4hana-ashlock/
ebookmass.com
https://ebookmass.com/product/978-0130387752-market-based-management/
ebookmass.com
https://ebookmass.com/product/advanced-health-assessment-and-
diagnostic-reasoning-3rd-edition-ebook-pdf/
ebookmass.com
https://ebookmass.com/product/pediatric-liver-transplantation-a-
clinical-guide-1st-edition-nedim-hadzic/
ebookmass.com
https://ebookmass.com/product/practical-applications-of-coaching-and-
mentoring-in-dentistry-janine-brooks/
ebookmass.com
Gilded Artist: A Golden Era Series Small Town Instalove
Romance Hayleigh Davis & Flirt Club
https://ebookmass.com/product/gilded-artist-a-golden-era-series-small-
town-instalove-romance-hayleigh-davis-flirt-club/
ebookmass.com
Orbital Mechanics for
Engineering Students
Orbital Mechanics for
Engineering Students
Fourth Edition
Howard D. Curtis
Professor Emeritus, Aerospace Engineering
Embry-Riddle Aeronautical University
Daytona Beach, Florida
Butterworth-Heinemann is an imprint of Elsevier
The Boulevard, Langford Lane, Kidlington, Oxford OX5 1GB, United Kingdom
50 Hampshire Street, 5th Floor, Cambridge, MA 02139, United States
# 2020 Elsevier Ltd. All rights reserved.
No part of this publication may be reproduced or transmitted in any form or by any means, electronic or
mechanical, including photocopying, recording, or any information storage and retrieval system, without
permission in writing from the publisher. Details on how to seek permission, further information about the
Publisher’s permissions policies and our arrangements with organizations such as the Copyright Clearance
Center and the Copyright Licensing Agency, can be found at our website: www.elsevier.com/permissions.
This book and the individual contributions contained in it are protected under copyright by the Publisher
(other than as may be noted herein).
MATLAB® is a trademark of The MathWorks, Inc. and is used with permission. The MathWorks does not
warrant the accuracy of the text or exercises in this book. This book’s use or discussion of MATLAB®
software or related products does not constitute endorsement or sponsorship by The MathWorks of a
particular pedagogical approach or particular use of the MATLAB® software.
Notices
Knowledge and best practice in this field are constantly changing. As new research and experience broaden
our understanding, changes in research methods, professional practices, or medical treatment may become
necessary.
Practitioners and researchers must always rely on their own experience and knowledge in evaluating and using
any information, methods, compounds, or experiments described herein. In using such information or methods
they should be mindful of their own safety and the safety of others, including parties for whom they have a
professional responsibility.
To the fullest extent of the law, neither the Publisher nor the authors, contributors, or editors, assume any liability
for any injury and/or damage to persons or property as a matter of products liability, negligence or otherwise,
or from any use or operation of any methods, products, instructions, or ideas contained in the material herein.
Library of Congress Cataloging-in-Publication Data
A catalog record for this book is available from the Library of Congress
British Library Cataloguing-in-Publication Data
A catalogue record for this book is available from the British Library
ISBN: 978-0-08-102133-0
The organization of the book remains the same as that of the third edition. Chapter 1 is a review of
vector kinematics in three dimensions and of Newton’s laws of motion and gravitation. It also focuses
on the issue of relative motion, crucial to the topics of rendezvous and satellite attitude dynamics. The
material on ordinary differential equation solvers will be useful for students who are expected to code
numerical simulations in MATLAB or other programming languages. Chapter 2 presents the vector-
based solution of the classical two-body problem, resulting in a host of practical formulas for the anal-
ysis of orbits and trajectories of elliptical, parabolic, and hyperbolic shape. The restricted three-body
problem is covered to introduce the notion of Lagrange points and to present the numerical solution of a
lunar trajectory problem. Chapter 3 derives Kepler’s equations, which relate position to time for the
different kinds of orbits. The universal variable formulation is also presented. Chapter 4 is devoted
to describing orbits in three dimensions. Coordinate transformations and the Euler elementary rotation
sequences are defined. Procedures for transforming back and forth between the state vector and the
classical orbital elements are addressed. The effect of the earth’s oblateness on the motion of an orbit’s
ascending node and eccentricity vector is described, pending a more detailed explanation in Chapter 10.
Chapter 5 is an introduction to preliminary orbit determination, including Gibbs’ and Gauss’ methods
and the solution of Lambert’s problem. Auxiliary topics include topocentric coordinate systems, Julian
xiii
xiv Preface
day numbering, and sidereal time. Chapter 6 presents the common means of transferring from one orbit
to another by impulsive delta-v maneuvers, including Hohmann transfers, phasing orbits, and plane
changes. Chapter 7 is a brief introduction to relative motion in general and to the two-impulse rendez-
vous problem in particular. The latter is analyzed using the Clohessy-Wiltshire equations, which are
derived in this chapter. Chapter 8 is an introduction to interplanetary mission design using patched
conics. Chapter 9 extends the patched conic method and the restricted three-body approach to lunar
trajectory analysis. Chapter 10 is an introduction to common orbital perturbations: drag, nonspherical
gravitational field, solar radiation pressure, and lunar and solar gravity. Chapter 11 presents those el-
ements of rigid body dynamics required to characterize the attitude of a space vehicle. Euler’s equa-
tions of rotational motion are derived and applied in a number of example problems. Euler angles, yaw-
pitch-roll angles, and quaternions are presented as ways to describe the attitude of rigid body.
Chapter 12 describes the methods of controlling, changing, and stabilizing the attitude of spacecraft
by means of thrusters, gyros, and other devices. Chapter 13 is a brief introduction to the characteristics
and design of multistage launch vehicles.
Chapters 1 through 4 form the core of a first orbital mechanics course. The time devoted to
Chapter 1 depends on the background of the student. It might be surveyed briefly and used thereafter
simply as a reference. What follows Chapter 4 depends on the objectives of the course.
Chapters 5 through 10 carry on with the subject of orbital mechanics. Chapter 6 on orbital maneu-
vers should be included in any case. Coverage of Chapters 5, 7, 8, and 9 is optional. However, if
Chapters 8 and 9 on interplanetary and lunar missions is to form a part of the course, then the solution
of Lambert’s problem (Section 5.3) must be studied beforehand.
Chapter 10 is appropriate for a course devoted exclusively to orbital mechanics with an introduction
to perturbations, which is a whole topic unto itself.
Chapters 11 and 12 must be covered if the course objectives include an introduction to spacecraft
dynamics. In that case Chapters 5, 7, 8, and 9 would probably not be studied in depth.
Chapter 13 is optional if the engineering curriculum requires a separate course in propulsion includ-
ing rocket dynamics.
The important topic of spacecraft control systems is omitted. However, the material in this book and
a course in control theory provide the basis for the study of spacecraft attitude control.
To understand the material and to solve problems requires using a lot of undergraduate mathemat-
ics. Mathematics, of course, is the language of engineering. Students must not forget that the English
mathematician and physicist Sir Isaac Newton (1642–1727) had to invent calculus so he could solve
orbital mechanics problems in more than just a heuristic way. Newton’s 1687 publication Mathemat-
ical Principles of Natural Philosophy (“the Principia”) is one of the most influential scientific works of
all time. It must be noted that his contemporary, the German mathematician Gottfried Wilhelm von
Leibnitz (1646–1716) is credited with inventing infinitesimal calculus independently of Newton in
the 1670s.
In addition to honing their math skills, students are urged to take advantage of computers (which,
incidentally, use the binary numeral system developed by Leibnitz). There are many commercially
available mathematics software packages for personal computers. Wherever possible they should be
used to relieve the burden of repetitive and tedious calculations. Computer-programming skills can
and should be put to good use in the study of orbital mechanics. The elementary MATLAB programs
referred to in Appendix D of this book illustrate how many of the procedures developed in the text can
Preface xv
be implemented in software. All the scripts were developed and tested using MATLAB version 9.2
(release 2017a). Information about MATLAB, which is a registered trademark of The MathWorks,
Inc., may be obtained from
The MathWorks, Inc.
3 Apple Hill Drive
Natick, MA 01760-2089, USA
www.mathworks.com
Appendix A presents some tables of physical data and conversion factors. Appendix B is a road map
through the first three chapters, showing how the most fundamental equations of orbital mechanics are
related. Appendix C shows how to set up the n-body equations of motion and program them in
MATLAB. Appendix D contains listings of all the MATLAB algorithms and example problems
presented in the text. Appendix E shows that the gravitational field of a spherically symmetric body
is the same as if the mass were concentrated at its center. Appendix F explains how to deal with a
computational issue that arises in some perturbation analyses.
ACKNOWLEDGEMENTS
Since the publication of the first three editions and during the preparation of this one, I have received
helpful criticism, suggestions, and advice from many sources locally and worldwide. I thank them all
and regret that time and space limitations prohibited the inclusion of some recommended additional
topics that would have enhanced the book.
It has been a pleasure to work with the people at Elsevier, in particular Joseph P. Hayton, Publisher;
Steve Merken, Senior Acquisitions Editor; and Nate McFadden, Senior Developmental Editor.
I appreciate their enthusiasm for the book, their confidence in me, and all the work they did to move
this project to completion.
Finally and most importantly, I must acknowledge the patience and support of my wife, Mary, who
was a continuous source of optimism and encouragement throughout the revision effort.
Howard D. Curtis
Daytona Beach, FL, United States
CHAPTER
1.2 VECTORS
A vector is an object that is specified by both a magnitude and a direction. We represent a vector graph-
ically by a directed line segment (i.e., an arrow pointing in the direction of the vector). The end opposite
the arrow is called the tail. The length of the arrow is proportional to the magnitude of the vector. Ve-
locity is a good example of a vector. We say that a car is traveling eastward at 80 km/h. The direction is
east and the magnitude, or speed, is 80 km/h. We will use boldface type to represent vector quantities
and plain type to denote scalars. Thus, whereas B is a scalar, B is a vector.
Orbital Mechanics for Engineering Students. https://doi.org/10.1016/B978-0-08-102133-0.00001-5
# 2020 Elsevier Ltd. All rights reserved.
1
2 CHAPTER 1 DYNAMICS OF POINT MASSES
FIG. 1.1
All of these vectors may be denoted A, since their magnitudes and directions are the same.
FIG. 1.2
Parallelogram rule of vector addition. A + B ¼ C.
Observe that a vector is specified solely by its magnitude and direction. If A is a vector, then all
vectors having the same physical dimensions, the same length, and pointing in the same direction as A
are denoted A, regardless of their line of action, as illustrated in Fig. 1.1. Shifting a vector parallel to
itself does not mathematically change the vector. However, the parallel shift of a vector might produce
a different physical effect. For example, an upward 5-kN load (force vector) applied to the tip of an
airplane wing gives rise to quite a different stress and deflection pattern in the wing than the same load
acting at the wing’s midspan.
The magnitude of a vector A is denoted kAk, or, simply A.
Multiplying a vector B by the reciprocal of its magnitude produces a vector that points in the di-
rection of B, but it is dimensionless and has a magnitude of one. Vectors having dimensionless mag-
nitude are called unit vectors. We put a hat (^) over the letter representing a unit vector. Then we can tell
simply by inspection that, for example, u ^ is a unit vector, as are B^ and ^e.
^A . As pointed out above,
It is convenient to denote the unit vector in the direction of the vector A as u
we obtain this vector from A as follows:
A
^A ¼
u (1.1)
A
^C ¼ C=C, u
Likewise, u ^F ¼ F=F, etc.
The sum or resultant of two vectors is defined by the parallelogram rule (Fig. 1.2). Let C be the sum
of the two vectors A and B. To form that sum using the parallelogram rule, the vectors A and B are
1.2 VECTORS 3
FIG. 1.3
Three-dimensional, right-handed Cartesian coordinate system.
shifted parallel to themselves (leaving them unaltered) until the tail of A touches the tail of B. Drawing
dotted lines through the head of each vector parallel to the other completes a parallelogram. The
diagonal from the tails of A and B to the opposite corner is the resultant C. By construction, vector
addition is commutative; that is,
A+B¼B+A (1.2)
A Cartesian coordinate system in three dimensions consists of three axes, labeled x, y, and z,
which intersect at the origin O. We will always use a right-handed Cartesian coordinate system,
which means if you wrap the fingers of your right hand around the z axis, with the thumb
pointing in the positive z direction, your fingers will be directed from the x axis toward the y axis.
Fig. 1.3 illustrates such a system. Note that the unit vectors along the x, y, and z axes are, respectively,
^i, ^j, and k.
^
In terms of its Cartesian components, and in accordance with the above summation rule, a vector A
is written in terms of its components Ax, Ay, and Az as
A ¼ Ax^i + Ay^j + Az k
^ (1.3)
According to the Pythagorean theorem, the magnitude of A in terms of its Cartesian components is
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
A¼ A2x + A2y + A2z (1.4)
From Eqs. (1.1) and (1.3), the unit vector in the direction of A is
^A ¼ cos θx^i + cos θy^j + cos θz k
u ^ (1.5)
4 CHAPTER 1 DYNAMICS OF POINT MASSES
FIG. 1.4
Direction angles in three dimensions.
where
Ax Ay Az
cos θx ¼ cos θy ¼ cos θz ¼ (1.6)
A A A
The direction angles θx, θy, and θz are illustrated in Fig. 1.4, and they are measured between the vector
and the positive coordinate axes. Note carefully that the sum of θx, θy, and θz is not in general known a
priori and cannot be assumed to be, say, 180 degrees.
EXAMPLE 1.1
Calculate the direction angles of the vector A ¼ ^i 4^j + 8k.
^
Solution
First, compute the magnitude of A by means of Eq. (1.4),
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
A ¼ 12 + ð4Þ2 + 82 ¼ 9
Then Eq. (1.6) yields
Ax 1
θx ¼ cos 1 ¼ cos 1 ) θx ¼ 83:62 degrees
A 9
Ay 4
θy ¼ cos 1 ¼ cos 1 ) θy ¼ 116:4 degrees
A 9
Az 8
θz ¼ cos 1 ¼ cos 1 ) θz ¼ 27:27 degrees
A 9
Observe that θx + θy + θz ¼ 227.3 degrees.
Multiplication and division of two vectors are undefined operations. There are no rules for com-
puting the product AB and the ratio A/B. However, there are two well-known binary operations on
1.2 VECTORS 5
vectors: the dot product and the cross product. The dot product of two vectors is a scalar defined as
follows:
A B ¼ AB cos θ (1.7)
where θ is the angle between the heads of the two vectors, as shown in Fig. 1.5. Clearly,
AB¼BA (1.8)
If two vectors are perpendicular to each other, then the angle between them is 90 degrees. It follows
from Eq. (1.7) that their dot product is zero. Since the unit vectors ^i, ^j, and k
^ of a Cartesian coordinate
system are mutually orthogonal and of magnitude 1, Eq. (1.7) implies that
^i ^i ¼ ^j ^j ¼ k
^k ^¼1
(1.9)
^i ^j ¼ ^i k^ ¼ ^j k
^¼0
Using these properties, it is easy to show that the dot product of the vectors A and B may be found in
terms of their Cartesian components as
A B ¼ Ax Bx + Ay By + Az Bz (1.10)
FIG. 1.5
The angle between two vectors brought tail to tail by parallel shift.
FIG. 1.6
Projecting the vector B onto the direction of A.
6 CHAPTER 1 DYNAMICS OF POINT MASSES
EXAMPLE 1.2
Let A ¼ ^i + 6^j + 18k
^ and B ¼ 42^i 69^j + 98k:
^ Calculate
(a) the angle between A and B;
(b) the projection of B in the direction of A;
(c) the projection of A in the direction of B.
Solution
First, we make the following individual calculations.
A B ¼ ð1Þð42Þ + ð6Þð69Þ + ð18Þð98Þ ¼ 1392 (a)
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
A ¼ ð1Þ2 + ð6Þ2 + ð18Þ2 ¼ 19 (b)
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
B ¼ ð42Þ2 + ð69Þ2 + ð98Þ2 ¼ 127 (c)
(a) According to Eq. (1.7), the angle between A and B is
AB
θ ¼ cos 1
AB
Substituting Eqs. (a), (b), and (c) yields
1392
θ ¼ cos 1 ¼ 54:77 degrees
19 127
(b) From Eq. (1.12), we find the projection of B onto A.
A AB
BA ¼ B ¼
A A
Substituting Eqs. (a) and (b) we get
1392
BA ¼ ¼ 73:26
19
(c) The projection of A onto B is
B AB
AB ¼ A ¼
B B
1.2 VECTORS 7
The cross product of two vectors yields another vector, which is computed as follows:
A B ¼ ðAB sinθÞ^
nAB (1.13)
where θ is the angle between the heads of A and B, and n^AB is the unit vector normal to the plane defined
by the two vectors. The direction of n^AB is determined by the right-hand rule. That is, curl the fingers of
the right hand from the first vector (A) toward the second vector (B), and the thumb shows the direction
^AB (Fig. 1.7). If we use Eq. (1.13) to compute B A, then n
of n ^AB points in the opposite direction,
which means
B A ¼ ðA BÞ (1.14)
Therefore, unlike the dot product, the cross product is not commutative.
The cross product is obtained analytically by resolving the vectors into Cartesian components.
A B ¼ Ax^i + Ay^j + Az k
^ Bx^i + By^j + Bz k
^ (1.15)
Expanding the right-hand side of Eq. (1.15), substituting Eq. (1.16), and making use of Eq. (1.14)
leads to
A B ¼ Ay Bz Az By ^i ðAx Bz Az Bx Þ^j + Ax By Ay Bx k
^ (1.17)
It may be seen that the right-hand side is the determinant of the matrix
2 3
^i ^j ^
k
4 Ax Ay Az 5
Bx By Bz
FIG. 1.7
^AB is normal to both A and B and defines the direction of the cross product A B.
n
8 CHAPTER 1 DYNAMICS OF POINT MASSES
where the two vertical bars stand for the determinant. Obviously, the rule for computing the cross prod-
uct, though straightforward, is a bit lengthier than that for the dot product. Remember that the dot prod-
uct yields a scalar whereas the cross product yields a vector.
The cross product provides an easy way to compute the normal to a plane. Let A and B be any two
vectors lying in the plane, or, let any two vectors be brought tail to tail to define a plane, as shown in
Fig. 1.7. The vector C ¼ A B is normal to the plane of A and B. Therefore, n ^ AB ¼ C=C, or
AB
^AB ¼
n (1.19)
kA Bk
EXAMPLE 1.3
Let A5 3^i + 7^j + 9k
^ and B56^i 5^j + 8k.
^ Find a unit vector that lies in the plane of A and B and is perpendicular to A.
Solution
The plane of vectors A and B is determined by parallel-shifting the vectors so that they meet tail to tail. Calculate the vector
D ¼ A B.
^i ^j k
^
D ¼ 3 7 9 ¼ 101^i + 78^j 27k
^
6 5 8
Note that A and B are both normal to D. We next calculate the vector C ¼ D A.
^i ^j ^
k
C ¼ 101 78 27 ¼ 891^i 828^j + 941k
^
3 7 9
C is normal to D as well as to A. A, B, and C are all perpendicular to D. Therefore, they are coplanar. Thus, C is not only
perpendicular to A, but it also lies in the plane of A and B. Therefore, the unit vector we are seeking is the unit vector in the
direction of C. That is
In the chapters to follow, we will often encounter the vector triple product, A (B C). By resolving
A, B, and C into their Cartesian components, it can easily be shown that the vector triple product can be
expressed in terms of just the dot products of these vectors as follows:
A ðB CÞ ¼ BðA CÞ CðA BÞ (1.20)
Because of the appearance of the letters on the right-hand side, this is often referred to as the “bac–cab
rule.”
1.3 KINEMATICS 9
EXAMPLE 1.4
If F ¼ E {D [A (B C)]}, use the bac–cab rule to reduce this expression to one involving only dot products.
Solution
First, we invoke the bac–cab rule to obtain
8 baccab rule
9
>
< zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}|fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{>
=
F ¼ E D ½BðA CÞ CðA BÞ
>
: >
;
Another useful vector identity is the “interchange of the dot and the cross”:
A ðB C Þ ¼ ðA B Þ C (1.21)
It is so-named because interchanging the operations in the expression A B C yields A B C.
The parentheses in Eq. (1.21) are required to show which operation must be carried out first, according
to the rules of vector algebra. (For example, (A B) C, the cross product of a scalar and a vector,
is undefined.) It is easy to verify Eq. (1.21) by substituting A ¼ Ax^i + Ay^j + Az k,
^ B ¼ Bx^i + By^j + Bz k,
^
^ ^ ^
and C ¼ Cx i + Cy j + Cz k and observing that both sides of the equal sign reduce to the same expression.
1.3 KINEMATICS
To track the motion of a particle P through Euclidean space, we need a frame of reference, consisting of
a clock and a nonrotating Cartesian coordinate system. The clock keeps track of time t, and the xyz axes
of the Cartesian coordinate system are used to locate the spatial position of the particle. In nonrelativ-
istic mechanics, a single “universal” clock serves for all possible Cartesian coordinate systems. So when
we refer to a frame of reference, we need to think only of the mutually orthogonal axes themselves.
The unit of time used throughout this book is the second (s). The unit of length is the meter (m), but
the kilometer (km) will be the length unit of choice when large distances and velocities are involved.
Conversion factors between kilometers, miles, and nautical miles are listed in Table A.3.
Given a frame of reference, the position of the particle P at a time t is defined by the position vector
r(t) extending from the origin O of the frame out to P itself, as illustrated in Fig. 1.8. The components of
r(t) are just the x, y, and z coordinates,
rðtÞ ¼ xðtÞ^i + yðtÞ^j + zðtÞk
^
10 CHAPTER 1 DYNAMICS OF POINT MASSES
FIG. 1.8
Position, velocity, and acceleration vectors.
The distance of P from the origin is the magnitude or length of r, denoted k r k or just r,
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
krk ¼ r ¼ x2 + y 2 + z2
As in Eq. (1.11), the magnitude of r can also be computed by means of the dot product operation,
pffiffiffiffiffiffiffiffi
r¼ rr
The velocity v and acceleration a of the particle are the first and second time derivatives of the position
vector,
dxðtÞ^ dyðtÞ^ dyðtÞ ^
vðtÞ ¼ i+ j+ k ¼ vx ðtÞ^i + vy ðtÞ^j + vz ðtÞk
^
dt dt dt
dvx ðtÞ^ dvy ðtÞ^ dvz ðtÞ ^
aðtÞ ¼ i+ j+ k ¼ ax ðtÞ^i + ay ðtÞ^j + az ðtÞk
^
dt dt dt
The derivatives of ^i, ^j, and k
^ are zero since axes of the Cartesian frame have fixed directions. It is con-
venient to represent the time derivative by means of an overhead dot. In this shorthand notation, if ( ) is
any quantity, then
dð Þ d2 ð Þ ⋯ d3 ð Þ
ð Þ¼ ð Þ¼ ð Þ¼ etc:
dt dt2 dt3
Thus, for example,
v ¼ r_
a ¼ v_ ¼ €r
vx ¼ x_ vy ¼ y_ vz ¼ z_
ax ¼ v_ x ¼ x€ ay ¼ v_ y ¼ y€ az ¼ v_ z ¼ z€
The locus of points that a particle occupies as it moves through space is called its path or trajectory.
If the path is a straight line, then the motion is rectilinear. Otherwise, the path is curved, and the motion
1.3 KINEMATICS 11
where the speed v is the magnitude of the velocity v. The distance ds that P travels along its path in the
time interval dt is obtained from the speed by
ds ¼ vdt
In other words,
v ¼ s_
The distance s, measured along the path from some starting point, is what the odometers in our auto-
mobiles record. Of course, s, _ our speed along the road, is indicated by the dial of the speedometer.
Note carefully that v 6¼ r_ (i.e., the magnitude of the derivative of r does not equal the derivative of
the magnitude of r).
EXAMPLE 1.5
The position vector in meters is given as a function of time in seconds as
r ¼ 8t2 + 7t + 6 ^i + 5t3 + 4 ^j + 0:3t4 + 2t2 + 1 k
^ ðmÞ (a)
At t ¼ 10 s, calculate (a) v (the magnitude of the derivative of r) and (b) r_ (the derivative of the magnitude of r).
Solution
(a) The velocity v is found by differentiating the given position vector with respect to time,
dr
v¼ ¼ ð16t + 7Þ^i + 15t2^j + 1:2t3 + 4t k
^
dt
The magnitude of this vector is the square root of the sum of the squares of its components,
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
v ¼ 1:44t6 + 234:6t4 + 272t2 + 224t + 49
Evaluating this at t ¼ 10 s, we get
v ¼ 1953:3m=s
(b) Calculating the magnitude of r in Eq. (a) leads to
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
r ¼ 0:09t8 + 26:2t6 + 68:6t4 + 152t3 + 149t2 + 84t + 53
The time derivative of this expression is
dr 0:36t7 + 78:6t5 + 137:2t3 + 228t2 + 149t + 42
r_ ¼ ¼ pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
dt 0:09t8 + 26:2t6 + 68:6t4 + 152t3 + 149t2 + 84t + 53
Substituting t ¼ 10 s yields
r_ ¼ 1935:5 m=s
FIG. 1.9
Orthogonal triad of unit vectors associated with the moving point P.
1.3 KINEMATICS 13
The center of curvature lies in the osculating plane. When the particle P moves an incremental distance
ds, the radial from the center of curvature to the path sweeps out a small angle, dϕ, measured in the
osculating plane. The relationship between this angle and ds is
ds ¼ ρdϕ
_ or
so that s_ ¼ ρϕ,
v
ϕ_ ¼ (1.30)
ρ
EXAMPLE 1.6
Relative to a Cartesian coordinate system, the position, velocity, and acceleration of a particle P at a given instant are
Substituting Eqs. (g) and (h) back into Eq. (e) finally yields the unit principal normal
^i ^j ^
k
^n ¼ 0:88409 0:0010104 0:46731 ¼ 0:25389^i + 0:8385^j 0:48214k
u ^ (i)
0:39233 0:5449 0:74106
The only unknown remaining in Eq. (d) is ρ, for which we appeal to Eq. (1.25),
v2
ρ¼ (j)
an
^n ,
The normal acceleration an is calculated by projecting the acceleration vector a onto the direction of the unit normal u
^n ¼ 16^i + 125^j + 30k
an ¼ a u ^ 0:25389^i + 0:8385^j 0:48214k ^ ¼ 86:287m=s2 (k)
Putting the values of v and an from Eqs. (f) and (k) into Eq. ( j) yields the radius of curvature,
229:42
ρ¼ ¼ 609:89m (l)
86:287
Upon substituting Eqs. (a), (i), and (l) into Eq. (d), we obtain the position vector of the center of curvature C,
rC ¼ 250^i + 630^j + 430k ^ + 609:89 0:25389^i + 0:8385^j 0:48214k ^
or
W ¼ mg (1.32)
where
GM
g¼ (1.33)
r2
g has units of acceleration (m/s2) and is called the acceleration of gravity. If planetary gravity is the only
force acting on a body, then the body is said to be in free fall. The force of gravity draws a freely falling
object toward the center of attraction (e.g., center of the earth) with an acceleration g. Under ordinary
conditions, we sense our own weight by feeling contact forces acting on us in opposition to the force of
gravity. In free fall, there are, by definition, no contact forces, so there can be no sense of weight. Even
though the weight is not zero, a person in free fall experiences weightlessness, or the absence of gravity.
Let us evaluate Eq. (1.33) at the surface of the earth, whose radius according to Table A.1 is
6378 km. Letting g0 represent the standard sea level value of g, we get
GM
g0 ¼ (1.34)
R2E
In SI units,
g0 ¼ 9:807m=s2 (1.35)
Substituting Eq. (1.34) into Eq. (1.33) and letting z represent the distance above the earth’s surface, so
that r ¼ RE + z, we obtain
R2E g0
g ¼ g0 ¼ (1.36)
ðRE + zÞ2 ð1 + z=RE Þ2
FIG. 1.10
Variation of the acceleration of gravity with altitude.
16 CHAPTER 1 DYNAMICS OF POINT MASSES
Commercial airliners cruise at altitudes on the order of 10 km (6 miles). At that height, Eq. (1.36)
reveals that g (and hence weight) is only three-tenths of a percent less than its sea level value.
Thus, under ordinary conditions, we ignore the variation of g with altitude. A plot of Eq. (1.36) out
to a height of 2000 km (the upper limit of low earth orbit operations) is shown in Fig. 1.10. The var-
iation of g over that range is significant. Even so, at space station altitude (400 km), weight is only
about 10% less than it is on the earth’s surface. The astronauts experience weightlessness, but they
clearly are not weightless.
EXAMPLE 1.7
Show that in the absence of an atmosphere, the shape of a low-altitude ballistic trajectory is a parabola. Assume the ac-
celeration of gravity g is constant and neglect the earth’s curvature.
Solution
Fig. 1.11 shows a projectile launched at t ¼ 0 s with a speed v0 at a flight path angle γ 0 from the point with coordinates
(x0, y0).
Since the projectile is in free fall after launch, its only acceleration is that of gravity in the negative y direction:
x€¼ 0
y€¼ g
Integrating with respect to time and applying the initial conditions leads to
x ¼ x0 + ðv0 cos γ 0 Þt (a)
1
y ¼ y0 + ðv0 sin γ 0 Þt gt2 (b)
2
Solving Eq. (a) for t and substituting the result into Eq. (b) yields
1 g
y ¼ y0 + ðx x0 Þtan γ 0 ðx x0 Þ2 (c)
2 v20 cos 2 γ 0
This is the equation of a second-degree curve, a parabola, as sketched in Fig. 1.11.
FIG. 1.11
Flight of a low-altitude projectile in free fall (no atmosphere).
1.5 NEWTON’S LAW OF MOTION 17
EXAMPLE 1.8
An airplane flies a parabolic trajectory like that in Fig. 1.11 so that the passengers will experience free fall (weightlessness).
What is the required variation of the flight path angle γ with speed v? Ignore the curvature of the earth.
Solution
Fig. 1.12 reveals that for a “flat” earth, dγ ¼ dϕ. That is,
γ_ ¼ ϕ_
It follows from Eq. (1.30) that
ρ_γ ¼ v (1.37)
The normal acceleration an is just the component of the gravitational acceleration g in the direction of the unit principal
normal to the curve (from P toward C). From Fig. 1.12, then,
an ¼ g cos γ (a)
Substituting the second of Eqs. (1.25) into Eq. (a) and solving for the radius of curvature yields
v2
ρ¼ (b)
gcos γ
Combining Eqs. (1.37) and (b), we find the time rate of change of the flight path angle,
g cos γ
γ_ ¼
v
FIG. 1.12
Relationship between dγ and dϕ for a “flat” earth.
In this equation, a is the absolute acceleration of the center of mass. The absolute acceleration is mea-
sured in a frame of reference that itself has neither translational nor rotational acceleration relative to
the fixed stars. Such a reference is called an absolute or inertial frame of reference.
Force is related to the primitive concepts of mass, length, and time by Newton’s second law. The
unit of force, appropriately, is the Newton, which is the force required to impart an acceleration of 1 m/
s2 to a mass of 1 kg. A mass of 1 kg therefore weighs 9.807 N at the earth’s surface. The kilogram is not
a unit of force.
Confusion can arise when mass is expressed in units of force, as frequently occurs in US engineer-
ing practice. In common parlance either the pound or the ton (2000 lb) is more likely to be used to
express the mass. The pound of mass is officially defined precisely in terms of the kilogram, as shown
in Table A.3. Since 1 lb of mass weighs 1 lb of force where the standard sea level acceleration of gravity
(Eq. 1.35) exists, we can use Newton’s second law to relate the pound of force to the Newton:
11b ðforceÞ ¼ 0:4536kg 9:807m=s2 ¼ 4:448N
The slug is the quantity of matter accelerated at 1 ft/s2 by a force of 1 lb. We can again use Newton’s
second law to relate the slug to the kilogram. Noting the relationship between feet and meters in
Table A.3, we find
1lb 4:448N kg m=s2
1 slug ¼ 2
¼ 2
¼ 14:59 ¼ 14:59kg
1ft=s 0:3048m=s m=s2
EXAMPLE 1.9
On a NASA mission, the space shuttle Atlantis orbiter was reported to weigh 239,255 lb just prior to liftoff. On orbit 18 at an
altitude of about 350 km, the orbiter’s weight was reported to be 236,900 lb. (a) What was the mass, in kilograms, of
Atlantis on the launchpad and in orbit? (b) If no mass was lost between launch and orbit 18, what would have been the
weight of Atlantis, in pounds?
Solution
(a) The given data illustrate the common use of weight in pounds as a measure of mass. The “weights” given are actually
the mass in pounds of mass. Therefore, prior to launch
0:4536kg
mlaunchpad ¼ 239,255lb ðmassÞ ¼ 108,500kg
1lbðmassÞ
In orbit,
0:4536kg
morbit 18 ¼ 236,900lb ðmassÞ ¼ 107,500kg
1lbðmassÞ
The decrease in mass is the propellant expended by the orbital maneuvering and reaction control rockets on the orbiter.
(b) Since the space shuttle launchpad at the Kennedy Space Center is essentially at sea level, the launchpad weight of
Atlantis in pounds (force) was numerically equal to its mass in pounds (mass). With no change in mass, the force of
gravity at 350 km would be, according to Eq. (1.36),
0 12
B 1 C
W ¼ 239,255lbðforceÞ @ A ¼ 215,000lbðforceÞ
350
1+
6378
1.5 NEWTON’S LAW OF MOTION 19
The integral of a force F over a time interval is called the impulse of the force,
ð t2
I¼ Fdt (1.39)
t1
Impulse is a vector quantity. From Eq. (1.38) it is apparent that if the mass is constant, then
ð t2
dv
I net ¼ m dt ¼ mv2 mv1 (1.40)
t1 dt
That is, the net impulse on a body yields a change mΔv in its linear momentum, so that
I net
Δv ¼ (1.41)
m
If Fnet is constant, then I net ¼ Fnet Δt, in which case Eq. (1.41) becomes
Fnet
Δv ¼ Δt ðif Fnet is constantÞ (1.42)
m
Let us conclude this section by introducing the concept of angular momentum. The moment of the
net force about O in Fig. 1.13 is
MO Þnet ¼ r Fnet
Substituting Eq. (1.38) yields
dv
MO Þnet ¼ r ma ¼ r m (1.43)
dt
But, keeping in mind that the mass is constant,
dv d dr d
rm ¼ ðr mvÞ mv ¼ ðr mvÞ ðv mvÞ
dt dt dt dt
Since v mv ¼ m(v v) ¼ 0, it follows that Eq. (1.43) can be written
dHO
MO Þnet ¼ (1.44)
dt
FIG. 1.13
The absolute acceleration of a particle is in the direction of the net force.
20 CHAPTER 1 DYNAMICS OF POINT MASSES
Thus, just as the net force on a particle changes its linear momentum mv, the moment of that force about
a fixed point changes the moment of its linear momentum about that point. Integrating Eq. (1.44) with
respect to time yields
ð t2
MO Þnet ¼ HO Þ2 HO Þ1 (1.46)
t1
The integral on the left is the net angular impulse. This angular impulse-momentum equation is the
rotational analog of the linear impulse-momentum relation given above in Eq. (1.40).
EXAMPLE 1.10
A particle of mass m is attached to point O by an inextensible string of length l, as illustrated in Fig. 1.14. Initially, the string
is slack when m is moving to the left with a speed v0 in the position shown. Calculate (a) the speed of m just after the string
becomes taut and (b) the average force in the string over the small time interval Δt required to change the direction of the
particle’s motion.
Solution
(a) Initially, the position and velocity of the particle are
FIG. 1.14
Particle attached to O by an inextensible string.
1.6 TIME DERIVATIVES OF MOVING VECTORS 21
FIG. 1.15
Displacement of a rigid body. (a) Change in orientation of an embedded vector A. (b) Differential rotation of A
about the instantaneous rotation axis.
the component of A normal to the axis of rotation. We can express the difference dA between A(t) and
A(t + dt) as
magnitude of dA
zfflfflfflfflfflfflfflfflfflfflfflffl}|fflfflfflfflfflfflfflfflfflfflfflffl{
dA ¼ ½ðkAk sinϕÞdθ n ^ (1.47)
where n^ is the unit normal to the plane defined by A and the axis of rotation, and it points in the direction
of the rotation. The angle ϕ is the inclination of A to the rotation axis. By definition,
dθ ¼ kωkdt (1.48)
where ω is the angular velocity vector, which points along the instantaneous axis of rotation, and its
direction is given by the right-hand rule. That is, wrapping the right hand around the axis of rotation,
with the fingers pointing in the direction of dθ, results in the thumb defining the direction of ω. This is
evident in Fig. 1.15(b). It should be pointed out that the time derivative of ω is the angular acceleration,
usually given the symbol α. Thus,
dω
α¼ (1.49)
dt
Substituting Eq. (1.48) into Eq. (1.47), we get
^ ¼ ðkωk kAk sinϕÞ^
dA ¼ kAk sinϕ kωkdt n n dt (1.50)
By definition of the cross product, ω A is the product of the magnitude of ω, the magnitude of A, the
sine of the angle between ω and A, and the unit vector normal to the plane of ω and A, in the rotation
direction. That is,
^
ω A ¼ kωk kAk sinϕ n (1.51)
Substituting Eq. (1.51) into Eq. (1.50) yields
dA ¼ ω Adt
1.6 TIME DERIVATIVES OF MOVING VECTORS 23
Eq. (1.52) is a formula we can use to compute the time derivative of any vector of constant magnitude.
EXAMPLE 1.11
Calculate the second time derivative of a vector A of constant magnitude, expressing the result in terms of ω and its de-
rivatives and A.
Solution
Differentiating Eq. (1.52) with respect to time, we get
d2 A d dA d dω dA
¼ ¼ ðω AÞ ¼ A+ω
dt2 dt dt dt dt dt
Using Eqs. (1.49) and (1.52), this can be written
d2 A
¼ α A + ω ðω AÞ (1.53)
dt2
EXAMPLE 1.12
Calculate the third derivative of a vector A of constant magnitude, expressing the result in terms of ω and its derivatives
and A.
Solution
d3 A d d2 A d
¼ ¼ ½α A + ω ðω AÞ
dt3 dt dt2 dt
d d
¼ ðα AÞ + ½ω ðω AÞ
dt
dt
dα dA dω d
¼ A+α + ðω AÞ + ω ðω AÞ
dt dt dt dt
dα dω dA
¼ A + α ðω AÞ + α ðω AÞ + ω A+ω
dt dt dt
dα
¼ A + α ðω AÞ + fα ðω AÞ + ω ½α A + ω ðω AÞg
dt
dα
¼ A + α ðω AÞ + α ðω AÞ + ω ðα AÞ + ω ½ω ðω AÞ
dt
dα
¼ A + 2α ðω AÞ + ω ðα AÞ + ω ½ω ðω AÞ
dt
d3 A dα
¼ A + 2α ðω AÞ + ω ½α A + ω ðω AÞ
dt3 dt
Let XYZ be a rigid inertial frame of reference and xyz a rigid moving frame of reference, as shown in
Fig. 1.16. The moving frame can be moving (translating and rotating) freely on its own accord, or it can
24 CHAPTER 1 DYNAMICS OF POINT MASSES
FIG. 1.16
Fixed (inertial) and moving rigid frames of reference.
be attached to a physical object, such as a car, an airplane, or a spacecraft. Kinematic quantities mea-
sured relative to the fixed inertial frame will be called absolute (e.g., absolute acceleration), and those
measured relative to the moving system will be called relative (e.g., relative acceleration). The unit
vectors along the inertial XYZ system are ^I, ^ ^ whereas those of the moving xyz system are
J, and K,
^i, ^j, and k.
^ The motion of the moving frame is arbitrary, and its absolute angular velocity is Ω. If, how-
ever, the moving frame is rigidly attached to an object, so that it not only translates but also rotates with
it, then the frame is called a body frame and the axes are referred to as body axes. A body frame clearly
has the same angular velocity as the body to which it is bound.
Let B be any time-dependent vector. Resolved into components along the inertial frame of refer-
ence, it is expressed analytically as
B ¼ BX^I + BY ^ ^
J + BZ K
where BX, BY, and BZ are functions of time. Since ^I, ^J, and K
^ are fixed, the time derivative of B is simply
dB dBX ^ dBY ^ dBZ ^
¼ I+ J+ K
dt dt dt dt
dBX/dt, dBY/dt, and dBZ/dt are the components of the absolute time derivative of B.
B may also be resolved into components along the moving xyz frame, so that, at any instant,
B ¼ Bx^i + By^j + Bz k
^ (1.54)
(unity) and, being attached to the xyz frame, they all have the angular velocity Ω. It follows from
Eq. (1.52) that
d^i d^j ^
dk
¼ Ω ^i ¼ Ω ^j ^
¼ Ωk
dt dt dt
Substituting these on the right-hand side of Eq. (1.55) yields
dB dBx ^ dBy ^ dBz ^
¼ i+ j+ k + Bx Ω ^i + By Ω ^j + Bz Ω k^
dt dt dt dt
dBx ^ dBy ^ dBz ^
¼ i+ j+ k + Ω Bx^i + Ω By^j + Ω Bz k ^
dt dt dt
dBx ^ dBy ^ dBz ^
¼ i+ j+ k + Ω Bx^i + By^j + Bz k
^
dt dt dt
In view of Eq. (1.54), this can be written as
dB dB
¼ + ΩB (1.56)
dt dt rel
where
dB dBx ^ dBy ^ dBz ^
¼ i+ j+ k (1.57)
dt rel dt dt dt
dB/dt)rel is the time derivative of B relative to the moving frame. Eq. (1.56) shows how the absolute
time derivative is obtained from the relative time derivative. Clearly, dB/dt ¼ dB/dt)rel only when the
moving frame is in pure translation (Ω ¼ 0).
Eq. (1.56) can be used recursively to compute higher order time derivatives. Thus, differentiating
Eq. (1.56) with respect to t, we get
d2 B d dB dΩ dB
¼ + B+Ω
dt2 dt dt rel dt dt
Using Eq. (1.56) in the last term yields
d2 B d dB dΩ dB
¼ + B+Ω + ΩB (1.58)
dt2 dt dt rel dt dt rel
Eq. (1.56) also implies that
d dB d2 B dB
¼ +Ω (1.59)
dt dt rel dt2 rel dt rel
where
d2 B d2 Bx ^ d2 By ^ d2 Bz ^
¼ i+ 2 j+ 2 k
dt2 rel dt2 dt dt
Substituting Eq. (1.59) into Eq. (1.58) yields
d2 B d2 B dB dΩ dB
¼ + Ω + B + Ω +ΩB
dt2 dt2 rel dt rel dt dt rel
26 CHAPTER 1 DYNAMICS OF POINT MASSES
FIG. 1.17
Absolute and relative position vectors.
1.7 RELATIVE MOTION 27
The cross product 2Ω vrel is called the Coriolis acceleration after Gustave Gaspard de Coriolis
(1792–1843), the French mathematician who introduced this term (Coriolis, 1835). Because of the
number of terms on the right, Eq. (1.70) is sometimes referred to as the five-term acceleration formula.
EXAMPLE 1.13
At a given instant, the absolute position, velocity, and acceleration of the origin O of a moving frame are
9
rO ¼ 100^I + 200^J + 300K
^ ðmÞ > >
=
vO ¼ 50^I + 30^J 10K^ ðm=sÞ ðgivenÞ (a)
>
>
aO ¼ 15^I + 40J^ + 25K^ m=s2 ;
Find (a) the velocity vrel and (b) the acceleration arel of P relative to the moving frame.
28 CHAPTER 1 DYNAMICS OF POINT MASSES
Solution
^ in terms of ^i, ^j, and k
Let us first use Eq. (c) to solve for ^I, ^J, and K ^ (three equations in three unknowns):
vrel ¼ v vO Ω rrel
^I J^ K^
¼ 70^I + 25^J 20K
^ 50^I + 30^J 10K
^ 1:0 0:4 0:6
200 300 150
¼ 70^I + 25J^ 20K
^ 50^I + 30J^ 10K
^ 240^I + 270^J 220K
^
or
so that
(b) To find the relative acceleration, we use the five-term acceleration formula, Eq. (1.70):
Fig. 1.18 shows the nonrotating inertial frame of reference XYZ with its origin at the center C of the
earth, which we shall assume to be a sphere. That assumption will be relaxed in Chapter 5. Embedded in
the earth and rotating with it is the orthogonal x0 y0 z0 frame, also centered at C, with the z0 axis parallel to
Z, the earth’s axis of rotation. The x0 axis intersects the equator at the prime meridian (0 degree
longitude), which passes through Greenwich in London, England. The angle between X and x0 is
θG, and the rate of increase of θG is just the angular velocity Ω of the earth. P is a particle (e.g., an
airplane or spacecraft), which is moving in an arbitrary fashion above the surface of the earth. rrel is
the position vector of P relative to C in the rotating x0 y0 z0 system. At a given instant, P is directly over
point O, which lies on the earth’s surface at longitude Λ and latitude ϕ. Point O coincides instanta-
neously with the origin of what is known as a topocentric-horizon coordinate system xyz. For our
purposes, x and y are measured positive eastward and northward along the local latitude and meridian,
FIG. 1.18
0 0 0
Earth-centered inertial frame (XYZ); earth-centered noninertial x y z frame embedded in and rotating with the
earth; and a noninertial, topocentric-horizon frame xyz attached to a point O on the earth’s surface.
30 CHAPTER 1 DYNAMICS OF POINT MASSES
respectively, through O. The tangent plane to the earth’s surface at O is the local horizon. The z axis is
the local vertical (straight up), and it is directed radially outward from the center of the earth. The unit
vectors of the xyz frame are ^i^jk,
^ as indicated in Fig. 1.18. Keep in mind that O remains directly below P,
so that as P moves, so do the xyz axes. Thus, the ^i^jk
^ triad, which comprises the unit vectors of a spherical
coordinate system, varies in direction as P changes location, thereby accounting for the curvature of
the earth.
Let us find the absolute velocity and acceleration of P. It is convenient to first obtain the velocity
and acceleration of P relative to the nonrotating earth, and then use Eqs. (1.66) and (1.70) to calculate
their inertial values.
The relative position vector can be written
^
rrel ¼ ðRE + zÞk (1.71)
where RE is the radius of the earth, and z is the height of P above the earth (i.e., its altitude). The time
derivative of rrel is the velocity vrel relative to the nonrotating earth,
^
vrel ¼
drrel ^ + ðRE + zÞ dk
¼ z_ k (1.72)
dt dt
^
To calculate dk=dt, we must use Eq. (1.52). The angular velocity ω of the xyz frame relative to the
nonrotating earth is found in terms of the rates of change of latitude ϕ and longitude Λ,
ω ¼ ϕ_ ^i + Λ_ cos ϕ^j + Λ_ sin ϕk
^ (1.73)
Thus,
^
dk ^ ¼ Λ_ cos ϕ^i + ϕ_ ^j
¼ ωk (1.74)
dt
Let us also record the following for future use:
d^j
¼ ω ^j ¼ Λ_ sin ϕ^i ϕ_ k
^ (1.75)
dt
d^i
¼ ω ^i ¼ Λ_ sin ϕ^j Λ_ cosϕk
^ (1.76)
dt
Substituting Eq. (1.74) into Eq. (1.72) yields the velocity in the nonrotating frame resolved along the
topocentric-horizon axes,
vrel ¼ x_^i + y_^j + z_ k
^ (1.77a)
where
x_ ¼ ðRE + zÞΛ_ cosϕ y_ ¼ ðRE + zÞϕ_ (1.77b)
It is convenient to use these results to express the rates of change of latitude and longitude in terms of
the components of relative velocity over the earth’s surface,
y_ x_
ϕ_ ¼ Λ_ ¼ (1.78)
RE + z ðRE + zÞ cos ϕ
1.7 RELATIVE MOTION 31
Substituting Eq. (1.74) through Eq. (1.76) together with Eqs. (1.78) and (1.79) into this expression
yields, upon simplification,
x_ ðz_ y_ tan ϕÞ ^ y_ z_ + x_ 2 tan ϕ ^ x_ 2 + y_ 2 ^
arel ¼ x€+ i + y€+ j + z€ k (1.80)
RE + z RE + z RE + z
Observe that the curvature of the earth’s surface is neglected by letting RE + z become infinitely large,
in which case
arel Þneglecting earth’s curvature ¼ x€^i + y€^j + z€k
^
That is, for a “flat earth,” the components of the relative acceleration vector are just the derivatives of
the components of the relative velocity vector.
For the absolute velocity we have, according to Eq. (1.66),
v ¼ vC + Ω rrel + vrel (1.81)
^ ¼ cos ϕ^j + sin ϕk,
From Fig. 1.18, it can be seen that K ^ which means the angular velocity of the earth is
^ ¼ Ω cos ϕ^j + Ω sin ϕk
Ω ¼ ΩK ^ (1.82)
Substituting this, together with Eqs. (1.71) and (1.77a) and the fact that vC ¼ 0, into Eq. (1.81) yields
v ¼ ½x_ + ΩðRE + zÞ cosϕ^i + y_^j + z_ k
^ (1.83)
From Eq. (1.70) the absolute acceleration of P is
_ rrel + Ω ðΩ rrel Þ + 2Ω vrel + arel
a ¼ aC + Ω
_ ¼ 0, we find, upon substituting Eqs. (1.71), (1.77a), (1.80), and (1.82), that
Since aC ¼ Ω
x_ ðz_ y_ tanϕÞ
a ¼ x€+ + 2Ωðz_ cos ϕ y_ sinϕÞ ^i
RE + z
y_ z_ + x_ tan ϕ
2
+ y€+ + Ω sinϕ½ΩðRE + zÞcos ϕ + 2x_ ^j (1.84)
RE + z
x_ 2 + y_ 2 ^
+ z€ Ωcosϕ½ΩðRE + zÞ cos ϕ + 2x_ k
RE + z
32 CHAPTER 1 DYNAMICS OF POINT MASSES
EXAMPLE 1.14
An airplane of mass 70,000 kg is traveling due north at a latitude 30°N, at an altitude of 10 km (32,800 ft), with a speed of
300 m/s (671 mph). Calculate (a) the components of the absolute velocity and acceleration along the axes of the
topocentric-horizon reference frame and (b) the net force on the airplane. Assume the winds aloft are zero.
Solution
(a) First, using the sidereal rotation period of the earth in Table A.1, we note that the earth’s angular velocity is
2π radians 2π radians 2π radians
Ω¼ ¼ ¼ ¼ 7:292 105 radians=s
sidereal day 23:93h 86,160 s
From Eq. (1.86a), the absolute velocity is
v ¼ ΩðRE + zÞcosϕ^i + y_^j ¼ 7:292 105 ð6378 + 10Þ 103 cos 30∘ ^i + 300^j
1.7 RELATIVE MOTION 33
FIG. 1.19
Components of the net force on the airplane.
or
3002
2
^
+ 7:292 105 ð6378 + 10Þ 103 cos 2 30∘ k
ð6378 + 10Þ 103
or
Fig. 1.19 shows the components of this relatively small force. The forward (y) and downward (negative z) forces are in
the directions of the airplane’s centripetal acceleration, caused by the earth’s rotation and, in the case of the downward
force, by the earth’s curvature as well. The westward force is in the direction of the Coriolis acceleration, which is due
to the combined effects of the earth’s rotation and the motion of the airplane. These net external forces must exist if the
airplane is to fly in the prescribed path.
In the vertical direction, the net force is that of the upward lift L of the wings plus the downward weight W of the
aircraft, so that
Fnet Þz ¼ L W ¼ 2769 ) L ¼ W 2769N
Other documents randomly have
different content
moon appears greater in the horizon, than when she is more
elevated. And because those reflected beams go transversely, there
will be generated, by the last article, redness. A possible cause
therefore is shown, why the moon as also the stars appear greater
and redder in the horizon, than in the midst of heaven. The same
also may be the cause, why the sun appears in the horizon greater
and of a colour more degenerating to yellow, than when he is higher
elevated. For the reflection from the little bodies between, and the
transverse motion of the medium, are still the same. But the light of
the sun is much stronger than that of the moon; and therefore, by
the last article, his splendour must needs by this perturbation
degenerate into yellowness.
But for the generation of these four colours, it is not necessary
that the figure of the glass be a prisma; for if it were spherical it
would do the same. For in a sphere the sunbeams are twice
refracted and twice reflected. And this being observed by Des
Cartes, and withal that a rainbow never appears but when it rains;
as also, that the drops of rain have their figures almost spherical; he
hath shown from thence the cause of the colours in the rainbow;
which therefore need not be repeated.
The cause of 15. Whiteness is light, but light perturbed by the
whiteness. reflections of many beams of light coming to the
eye together within a little space. For if glass or any other
diaphanous body be reduced to very small parts by contusion or
concussion, every one of those parts, if the beams of a lucid body be
from any one point of the same reflected to the eye, will represent
to the beholder an idea or image of the whole lucid body, that is to
say, a phantasm of white. For the strongest light is the most white;
and therefore many such parts will make many such images.
Wherefore, if those parts lie thick and close together, those many
images will appear confusedly, and will by reason of the confused
light represent a white colour. So that from hence may be deduced a
possible cause, why glass beaten, that is, reduced to powder, looks
white. Also why water and snow are white; they being nothing but a
heap of very small diaphanous bodies, namely, of little bubbles, from
whose several convex superficies there are by reflection made
several confused phantasms of the whole lucid body, that is to say,
whiteness. For the same reason, salt and nitre are white, as
consisting of small bubbles which contain within them water and air;
as is manifest in nitre, from this, that being thrown into the fire it
violently blows the same; which salt also doth, but with less
violence. But if a white body be exposed, not to the light of the day,
but to that of the fire or of a candle, it will not at the first sight be
easily judged whether it be white or yellow; the cause whereof may
be this, that the light of those things, which burn and flame, is
almost of a middle colour between whiteness and yellowness.
The cause of 16. As whiteness is light, so blackness is the
blackness. privation of light, or darkness. And, from hence it is,
first, that all holes, from which no light can be reflected to the eye,
appear black. Secondly, that when a body hath little eminent
particles erected straight up from the superficies, so that the beams
of light which fall upon them are reflected not to the eye but to the
body itself, that superficies appears black; in the same manner as
the sea appears black when ruffled by the wind. Thirdly, that any
combustible matter is by the fire made to look black before it shines.
For the endeavour of the fire being to dissipate the smallest parts of
such bodies as are thrown into it, it must first raise and erect those
parts before it can work their dissipation. If, therefore, the fire be
put out before the parts are totally dissipated, the coal will appear
black; for the parts being only erected, the beams of light falling
upon them will not be reflected to the eye, but to the coal itself.
Fourthly, that burning glasses do more easily burn black things than
white. For in a white superficies the eminent parts are convex, like
little bubbles; and therefore the beams of light, which fall upon
them, are reflected every way from the reflecting body. But in a
black superficies, where the eminent particles are more erected, the
beams of light falling upon them are all necessarily reflected towards
the body itself; and, therefore, bodies that are black are more easily
set on fire by the sun beams, than those that are white. Fifthly, that
all colours that are made of the mixture of white and black proceed
from the different position of the particles that rise above the
superficies, and their different forms of asperity. For, according to
these differences, more or fewer beams of light are reflected from
several bodies to the eye. But in regard those differences are
innumerable, and the bodies themselves so small that we cannot
perceive them; the explication and precise determination of the
causes of all colours is a thing of so great difficulty, that I dare not
undertake it.
Vol. I. Lat. & Eng.
C. XXVII.
Fig. 1-2]
Fig 1 Fig 2
CHAPTER XXVIII.
1.. The definition of sound, and the distinctions of sounds.—2. The cause of the
degrees of sounds.—3. The difference between sounds acute and grave.—4.
The difference between clear and hoarse sounds, whence.—5. The sound of
thunder and of a gun, whence it proceeds.—6. Whence it is that pipes, by
blowing into them, have a clear sound.—7. Of reflected sound.—8. From
whence it is that sound is uniform and lasting.—9. How sound may be helped
and hindered by the wind.—10. Not only air, but other bodies how hard
soever they be, convey sound.—11. The causes of grave and acute sounds,
and of concent.—12. Phenomena for smelling.—13. The first organ and the
generation of smelling.—14. How it is helped by heat and by wind.—15. Why
such bodies are least smelt, which have least intermixture of air in them.—16.
Why odorous things become more odorous by being bruised.—17. The first
organ of tasting; and why some savours cause nauseousness.—18. The first
organ of feeling; and how we come to the knowledge of such objects as are
common to the touch and other senses.
Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.
ebookmass.com