Mathematics For Mechanical Engineers
Mathematics For Mechanical Engineers
MECHANICAL ENGINEERS
NGUYEN VIET KHOA
I
PREFACE
Science and mathematics are integral parts of engineering. Science
provides the laws of the natural world and describes them as
relationships between different quantities using mathematical
formulas. Mechanical engineering is all about applying science and
mathematics to solve problems we face in our daily lives such as
mechanics of solids, mechanics of fluids, mechanics of structures and
so on. This text book is written to provide mathematical knowledge to
meet the needs of engineering students in the early years of mechnical
engineering degrees provided they have learned calculus I and II.
However, the book will also be valuable to a much wider readership
including those on other courses, e.g. engineering diplomas, and those
who already have some basic knowledge of mechanics, but wish to
explore more advanced topics, such as works done by variable forces,
circulations of fluid in vector fields, vibrations of structures.
The book contains ten chapters. Scalar fields and vector fields are the
subjects of Chapter 1, as they constitute the knowledge of the vector
and scalar fields that is very important tool for describing many
physical concepts such as gravitation, velocity fields of wind and fluid
which affect the behavior of objects over a large region of space. This
chapter presents the basic knowledge of vectors, vector fields and
their properties. The vector fields will be applied to form engineering
problems, e.g. to define work done by a force field in moving an
object along a curve as line integrals, to define the rate of fluid flow
across as surface integrals.
Students can seek more materials beyond the book from a list of
references provided at the end of the book.
Acknowledgments
I would like to thank Thao Thi Bich Dao for her careful reading the
text, checking the examples, and drawing and modifying some of the
figures to make it readable.
CONTENTS
INTRODUCTION ............................................................................... I
CHAPTER 1. SCALAR FIELDS, VECTOR FIELDS ................ 1
1.1 Scalar Fields............................................................................... 2
1.1.1 Definition of Scalar Field ................................................ 2
1.1.2 Representations of Scalar Fields...................................... 3
1.2 Vectors ....................................................................................... 6
1.2.1 Adding Vectors ................................................................ 6
1.2.2 Triangle Law.................................................................... 7
1.2.3 Parallelogram Law ........................................................... 7
1.2.4 Scalar ............................................................................... 8
1.2.5 Subtracting Vectors ......................................................... 8
1.2.6 Components of Vectors in 2- or 3-Dimensional.............. 8
1.2.7 Length of Vector ............................................................ 10
1.2.8 Adding Algebraic Vectors ............................................. 10
1.2.9 Multiplying Algebraic Vectors ...................................... 11
1.2.10 3-D Algebraic Vectors ................................................. 11
1.2.11 Components of n-Dimensional Vectors ...................... 12
1.2.12 Properties of Vectors ................................................... 12
1.2.13 Unit Vector .................................................................. 12
1.2.14 Dot Product .................................................................. 13
1.2.15 Direction Angles and Cosines ..................................... 15
1.2.16 Projections ................................................................... 16
1.2.17 Cross Product............................................................... 18
1.2.18 Cross Product in Physics ............................................. 20
1.2.19 Equations of Lines and Planes ..................................... 21
1.2.20 Scalar Equations of a Line ........................................... 22
1.2.21 Equation of Line Segments.......................................... 25
1.2.22 Equation of Planes ....................................................... 26
1.2.23 Parallel planes .............................................................. 28
1.2.24 Nonparallel planes ....................................................... 28
1.2.25 Distance from a Point to a Plane ................................. 30
1.3 Vector Fields ............................................................................ 32
1.3.1 Vector Functions............................................................ 32
1.3.2 Limit of a Vector Function ............................................ 32
1.3.3 Continuous of Vector Function and Space Curve ......... 33
1.3.4 Derivatives of Vector Functions .................................... 34
1.3.5 Integrals of Vector Functions ........................................ 35
1.3.6 Vector Fields.................................................................. 36
1.3.7 Vector Fields on 2 ...................................................... 36
1.3.8 Vector Fields on 3 ....................................................... 38
1.4 Velocity Field .......................................................................... 38
1.5 Gravitational Field ................................................................... 40
1.6 Conservative Vector Fields...................................................... 42
CHAPTER 2. GRADIENT, DIVERGENCE, ROTATIONAL
AND LAPLACIAN OPERATORS ..................... 53
2.1 Gradient Operator .................................................................... 53
2.1.1 Partial Derivative ........................................................... 53
2.1.2 Gradient Operator .......................................................... 56
2.1.3 Directional Derivative ................................................... 57
2.1.4 Tangent Planes to a Level Surfaces ............................... 60
2.1.5 Normal Line ................................................................... 61
2.2 Curl Operator ........................................................................... 62
2.2.1 Conservative Vector Fields ........................................... 64
2.2.1 Physical Meaning of Curl .............................................. 66
2.3 Divergence Operator ................................................................ 67
2.4 Laplace Operator ..................................................................... 68
CHAPTER 3. MULTIPLE INTEGRALS ................................ 73
3.1 Double Integrals over Rectangles ............................................ 73
3.2 Iterated Integrals ...................................................................... 78
3.3 Double Integrals over General Regions ................................... 79
3.3.1 Type I Region ................................................................ 80
3.3.2 Type II region ................................................................ 82
3.4 Double Integrals in Polar Coordinates ..................................... 89
3.5 Triple Integrals......................................................................... 97
3.5.1 Definition of Triple Integrals ......................................... 97
3.5.2 Integral Over a General Bounded Region ................... 100
3.5.3 Applications of Triple Integrals ................................... 107
3.6 Triple Integrals in Cylindrical Coordinates ........................... 110
3.7 Triple Integrals in Spherical Coordinates .............................. 114
3.8 Change of Variables in Multiple Integrals ............................. 118
3.8.1 Change of Variables in Double Integrals .................... 118
3.8.2 Change of Variables in Triple Integrals ...................... 130
CHAPTER 4. DIFFERENTIAL FORMS ...............................139
4.1 1-Forms .................................................................................. 139
4.2 Exactness in R2 ...................................................................... 140
4.3 2-Forms .................................................................................. 148
4.4 Exactness in R3 ...................................................................... 150
4.4.1 Convert a Vector Field to a 2-Form............................. 151
4.4.2 Derivative of 1-Forms ................................................. 152
4.4.3 Curl and Derivative of 1-Forms................................... 152
4.4.4 Derivative of 2-Forms ................................................. 153
4.4.5 Surface Integrals .......................................................... 154
CHAPTER 5. LINE INTEGRALS AND CIRCULATION OF A
FIELD .............................................................165
5.1 Line Integrals ......................................................................... 165
5.1.1 Line Integrals in 2-D Coordinates ............................... 165
5.1.2 Interpretation of a Line Integral................................... 167
5.1.3 Piecewise-smooth Curve ............................................. 168
5.1.4 Mass ............................................................................. 169
5.1.5 Curve Orientation ........................................................ 173
5.1.6 Line Integrals in 3-D Coordinates ............................... 175
5.1.7 Line Integrals of Vector Fields .................................... 178
5.1.8 Line Integral in Differential Forms and Vector Fields 165
5.2 Fundamental Theorem for Line Integrals .............................. 183
5.3 Line Integrals as Circulation.................................................. 186
5.4 Physical Interpretation of Line Integrals ............................... 192
CHAPTER 6. FLOW OF A FIELD THROUGH A SURFACE 203
6.1 Parametric Surfaces ............................................................... 203
6.1.1 Tangent Planes ............................................................. 211
6.1.2 Smooth Surface............................................................ 211
6.1.3 Surface Area ................................................................ 213
6.2 Surface Integral ...................................................................... 218
6.3 Flow of a Field Through a Surface ........................................ 227
CHAPTER 7. GREEN, STOKE AND OSTROGRADSKI
THEOREMS....................................................239
7.1 Green Theorem ...................................................................... 239
7.1.1 Reverse Direction ........................................................ 244
7.1.2 Regions with Holes ...................................................... 247
7.2 Stokes Theorem ..................................................................... 249
7.3 Ostrogradski Theorem ........................................................... 257
CHAPTER 8. FOURIER SERIES ..........................................271
8.1 Fourier Series Definition ....................................................... 271
8.2 Lemma ................................................................................... 271
8.3 Fourier Convergence Theorem .............................................. 272
8.4 Functions with Period 2L....................................................... 279
8.5 Complex form of Fourier series ............................................. 284
8.6 Complex form of Fourier series ............................................. 286
CHAPTER 9. FOURIER TRANSFORM ................................299
9.1 From Fourier Series to Fourier Transform ............................ 273
9.2 Theorems of Fourier Transform ............................................ 301
9.3 The Sampling Theorem ......................................................... 313
9.4 Discrete Fourier Transform ................................................... 315
9.4.1 Forward Transform ...................................................... 317
9.4.2 Inverse Transform........................................................ 318
9.6 Spectrum Analysis ................................................................. 319
CHAPTER 10. LAPLACE TRANSFROM .............................327
10.1 Definition, Basic Properties ................................................. 327
10.2 Further Laplace Transform Theorems ................................. 334
10.3 The Inverse Laplace Transform ........................................... 338
10.4 Laplace Transform for Solving Ordinary Differential
Equations .............................................................................. 344
10.4.1 First Order Differential Equations ............................. 345
10.4.2 Second Order Differential Equations ........................ 348
10.4.3 Some Applications of Laplace Transform to Mechanical
Engineering........................................................................... 349
REFERENCES .....................................................................371
Scalar and vector fields 1
This chapter will give the definitions and examples of scalar and
vector quantities. As you proceed through this text book, give careful
attention to the scalar and vector nature of each quantity. This is very
important for engineers to understand the nature of engineering
problems. The concepts of scalar fields are presented first since it is
familiar with the students. The vector fields will be presented later in
details because it will provide a new mathematical presentation of
engineering problems.
1 2
f x, y , z (1.1)
x2 y d z 2 x2 y d z2
2 2
The first way to a presentation is to use the contour map. In this way
we need to fix one of independent variables (z, for example) and then
show a contour map for the two remaining dimensions. In contour
map, each curve represents a line of constant values of the function f.
We show such a contour map in the xy-plane at z = 0 for Equation 1.1.
1 2
f x, y,0
x2 y d x2 y d
2 2
Scalar and vector fields 5
The Figure 1.3 shows a contour map consisting of 100 contour levels
with d =0.5. A series of these maps for various (fixed) values of z will
give a feel for the properties of the scalar function.
1.1.2.2 Color-Coding
1.2 Vectors
In engineering problems, vector is used to indicate a quantity that has
both magnitude and direction (such as acceleration, velocity,
displacement or force). A vector is in general represented by an arrow
or a directed line segment. The length of the arrow represents the
magnitude of the vector, while the arrow points in the direction of the
vector. In this text book, we denote a vector by either: printing a letter
in boldface, for example v, or putting an arrow above the letter, for
example v .
B
Figure 1.6 presents a vector which
describes the velocity vector v of an
object moving from initial point A
(the tail) to terminal point B (the tip).
A
We can write this vector as follows:
Figure 1.6. A vector
v AB
If the vector u CD has the same length and the same direction as v
even though it is in a different position then we say u and v are equal
and write u v as shown in Figure 1.7.
B D
A C
The
resulting displacementvector
AB is called the sum of AM and
MB . We write:
M
A
AB AM MB
Figure 1.8. Adding vectors
1.2.2 Triangle Law
The above law for adding two vectors
is sometime called the Triangle Law. B
Figure 1.9 will explain the reason why
they call this name. As can be seen
from this picture vector AB
corresponding to the side AB of the A
triangle ABC has the same initial and M
terminal points with the combination of
two vectors AM MB which
Figure 1.9. Triangle Law
correspond to two other sides of the
triangle AM and MB.
u v vu
1.2.4 Scalar
In this context, we call the real
number c a scalar to distinguish it
from a vector. For example, we
want to describe a vector 2v to be
2v
the same as vector v v , which has
the same direction as v but is twice v
as long. In general, we multiply a 0.5v
vector by a scalar as follows. If c is
a scalar and v is a vector, the scalar
-v -1.5v
multiple cv is the vector whose
length is |c| times the length of v
and whose direction is the same as
v if c 0 and is opposite to v if
c 0 . Look at Figure 1.11, we see
that real numbers work like scaling Figure 1.11. Scalar multiplies of v
factors here. That’s why we call
them scalars.
u
1.2.5 Subtracting Vectors v
u-v
The difference u – v of two vectors
implies that:
-v
u – v = u + (–v)
y
z
(a1, a2)
a2
O a1 x
x
y
a = 〈a1, a2〉 a = 〈a1, a2, a3〉
a a1 , a2 or a a1 , a2 , a3
x1 a1 x2
y1 a2 y2
z1 a3 z2
Thus,
a1 x2 x1
a2 y2 y1
a3 z2 z1
10 Mathematics for mechanical engineers
a x2 x1 , y2 y1 , z2 z1
a a12 a22
a+b a1 b1 , a2 b2 a+b b b2
b1
at least for the case where the a
components are positive. This a2 a2
means that, to add algebraic O a1 b1 x
vectors, we add their
components.
Figure 1.14. Adding Algebraic Vectors
Similarly, to subtract vectors, we subtract components.
a-b a1 b1 , a2 b2
Scalar and vector fields 11
ca ca1 , ca2
a1 , a2 , a3 b1 , b2 , b3 a1 b1 , a2 b2 , a3 b3
a1 , a2 , a3 b1 , b2 , b3 a1 b1 , a2 b2 , a3 b3
c a1 , a2 , a3 ca1 , ca2 , ca3
Solution: We have:
a 12 22 32 = 14
a + b = 1, 2, 3 + 2, 0, 2
= 1 + ( 2), 2 + 0, 3 + 2
= 1, 2, 5
a b = 1, 2, 3 2, 0, 2
= 1 ( 2), 2 0, 3 2
= 3, 2,1
2a = 21, 2, 3
= 2.(1), 2(2), 2(3)
= 2, 4, 6
12 Mathematics for mechanical engineers
2a 3b = 21, 2, 3 32, 0, 2
= 2, 4, 6 6, 0, 6
= 8, 4, 0
From now on we denote V2 as the set of all 2-D vectors, V3 as the set
of all 3-D vectors. More generally, we will later need to consider the
set Vn of all n-dimensional vectors. An n-dimensional vector is
defined as:
a a1 , a2 , , an
where a1 , a2 ,, an are real numbers that are called the components of
a.
1. a + b = b + a 2. a + (b + c) = (a + b) + c
3. a + 0 = a 4. a + ( a) = 0
5. c(a + b) = ca + cb 6. (c + d )a = ca + da
7. (cd )a = c (da) 8. 1a = a
1 a
| u | a
|a| |a|
1
| u || ca || c || a | | a | 1
|a|
Example 1.2: Find the unit vector in the direction of the vector
i 2 j 2k
| i 2 j 2k | 12 (2) 2 ( 2) 2 9 3
1 1 2 2
(i 2 j 2k ) i j k
3 3 3 3
a b a1 , a2 b1 , b2 a1b1 a1b2
a b 1,3 3, 1 1 3 3 1 0
c d 1,3, 2 2, 2,1 1 2 3 2 2 1 10
e f i j 2k i 2 j k 1 1 1 2 2 1 1
1. a a=|a|2
2. a b b a
3. a (b c) a b a c
4. (ca) b c(a b) a (cb)
5. 0 a 0
B
a-b
b
θ A
O a
x
y
a b a b cos
2 2 2 2 2
AB OA OB 2 OA OB cos a b 2 a b cos
In addition:
AB a-b a b a b a a a b b a b b
2 2
2 2
= a 2a b b
Thus,
2a b 2 a b cos
or
a b a b cos
a b
cos
| a || b |
ai a
cos 1 (1.2)
| a || i | | a |
a2 a
cos , cos 3 (1.3)
|a| |a|
16 Mathematics for mechanical engineers
a3
a
γ
β a2
a1
α y
x
a a1 , a2 , a3
a cos , a cos , a cos
a cos , cos , cos
Thus,
1
a cos , cos , cos
|a|
This states that the direction cosines of a are the components of the
unit vector in the direction of a.
1.2.16 Projections
Scalar and vector fields 17
Figure 17 depicts representations PQ and PR of two vectors a and b
with the same initial point P. Let S be the foot of the perpendicular
from R to the line containing PQ . Then, the vector with
representation PS is called the vector projection of b onto a and is
denoted by proja b. Then,
PS b cos
a b a b cos a b cos
shows that the dot product of a and b can be interpreted as the length
of a multiply by the scalar projection of b onto a, i.e.
a b a
| b | cos b
|a| |a|
θ S a
P
Q
|b|cos θ = compab
Example 1.4: Find the work done by the force of 100 N applied to pull
a wagon a distance of 150 m along a horizontal path, provided that the
handle of the wagon is held at an angle of 30° above the horizontal.
18 Mathematics for mechanical engineers
W F D F D cos 30o
100 150 cos 30 12990 Nm 12990 J
30° F
30°
D
Now we learn that unlike the dot product, the cross product a b of
two vectors a and b is a vector. For this reason, it is also called the
vector product. Note that a b is defined only when a and b are three-
dimensional (3-D) vectors.
The reason for the particular form of this definition is that the cross
product defined in this way has many useful properties. In order to
make this definition easier to remember, we use the notation of
determinants. A determinant of order 2 is defined by:
a b
ad bc
c d
a1 a2 a3
b b3 b b b b
b1 b2 b3 a1 2 a2 1 3 a3 1 2 (1.4)
c2 c3 c1 c3 c1 c2
c1 c2 c3
a2 a3 a1 a3 a1 a2
ab i j k (1.5)
b2 b3 b1 b3 b1 b2
In view of the similarity between the Equations 1.4 and 1.5, we can
write:
i j k
a b a1 a2 a3
b1 b2 b3
a2 a3 a a3 a a2
(a b ) a a1 1 a2 1 a3
b2 b3 b1 b3 b1 b2
a1 (a2b3 a3b2 ) a2 (a1b3 a3b1 ) a3 (a1b2 a2b1 )
a1a2b3 a1b2 a3 a1a2b3 b1a2 a3 a1b2 a3 b1a2 a3 0
a b a b sin
Corollary 1.2: Two nonzero vectors a and b are parallel if and only if:
ab 0
If we apply Theorems 1.3 and 1.4 to the standard basis vectors i, j, and
k using / 2 , we obtain:
i j k, j k i, k i j, j i k, k j i, i k j
1. a b b a
2. ca b c a b a cb
3. a b c a b a c
4. a b c a c b c
5. a b c a b c
6. a b c a c b a b c
τ r F r F sin
where θ is the angle between the position vector and force vector. As
can be seen from Figure 1.20, vector F can have two components: one
is in the direction of r and another is perpendicular to r. We can
observe that the only component of F that can cause a rotation is the
one perpendicular to r, that is, F sin . The magnitude of the torque
is equal to the area of the parallelogram determined by r and F.
Example 1.5: Find the magnitude of the torque about the center of a
bolt given that the bolt is tightened by applying a 60 N force to a 0.2
m wrench, as shown in Figure 1.21.
τ r F r F sin 70o
0.2 60 sin 70o Figure 1.21. A bolt is
12sin 70 11.28 N.m
o
tightened by a wrench
r = r0 + a
a = tv
Thus,
r = r0 + t v (1.6)
x x0 at
y y0 bt (1.7)
z z0 ct
these equations for t, equate the results, then we can derive the
following equations:
x x0 y y0 z z0
(1.8)
a b c
a 1, b –5, c 4
x 2 t , y 4 – 5t , z –3 4t
x2 y 4 z 3
1 5 4
We can see from Example 1.6 that direction numbers of the line L
through the points P0 x0 , y0 , z0 and P1 x1 , y1 , z1 are:
x1 – x0, y1 – y0, z1 – z0
x x0 y y0 z z0
x1 x0 y1 y0 z1 z0
Example 1.7:
Solution:
a 4, b 3, c 5
x3 y 5 2
4 3 5
Figure 1.23. Vector AB
We know from Equation 1.6 that the vector equation of a line through
the tip of the vector r0 in the direction of a vector v is r r0 tv . If the
line also passes through the tip of r1, then we can take v r1 – r0 . So,
its vector equation is:
r r0 t r1 – r0 1 – t r0 tr1
26 Mathematics for mechanical engineers
n r r0 0
a, b, c x x0 , y y0 , z z0 0
a x x0 b y y0 c z z0 0 (1.9)
Example 1.8: Find an equation of the plane through the point 2,3,1
with normal vector n 3, 4, 6 . Find the intercepts and sketch the
plane.
Solution: Putting
a 3, b 4, c 6, x0 2, y0 3, z0 1
3 x – 2 4 y – 3 6 z – 1 0
or
z
3x 4 y 6 z 24 (0,0,4)
ax by cz d 0 (1.10)
Conversely, it can be shown that, if a, b, and c are not all 0, then the
linear Equation 1.10 represents a plane with normal vector a, b, c .
3x 4 y 6 z 24 and –6 x – 8 y –12 z 12
1.2.24 Nonparallel
planes
n1
n2
If two planes are not parallel,
then they intersect in a θ
straight line. The angle
between the two planes is
defined as the acute angle θ
between their normal vectors
as illustrated in Figure 1.25
Figure 1.25. Nonparallel planes
Example 1.9: Find:
Solution
a) From the given plane equations, the normal vectors of these planes
are:
154
b) From the symmetric equations for the line we see that the line L is
determined if we know a point on it and a vector parallel to it. So we
first need to find a point on L then find vector v parallel to L:
One way is to find the point where the line intersects the xy-
plane by setting z 0 in the equations of both planes.
This provides the equations x 2 y 1 and 3x – y 3
whose solution is x 1, y 0 .
So, the point 1, 0, 0 lies on L.
i j k
v n1 n 2 1 2 3 5 i 8 j 7 k
3 1 1
x 1 y z
5 8 7
30 Mathematics for mechanical engineers
b x1 x0 , y1 y0 , z1 z0
Since P0 lies in the plane, its coordinates satisfy the equation of the
plane, hence:
Thus,
L1 : x –1 2t , y 1 t , z 3 – t
L2 : x 2 – s, y 1 3s, z –3 2 s
Solution: Since the two lines L1 and L2 are not intersected, they can
be viewed as lying on two parallel planes P1 and P2 as shown in
Figure 1.27. Thus, the distance between L1 and L2 is equal to the
distance between P1 and P2 . Thus, to solve the problem we will find a
point on P1 and an equation of P2 then apply Equation 1.11. We have
the direction of L1 and L2 are:
v1 2,1, –1 (direction of L1 )
n
v 2 –1,3, 2 (direction of L2 )
v2
L2
So the normal vector of the two lines
is:
n v1 v 2 v1
L1
i j k Figure 1.27. Parallel planes
2 1 1 5 i 3 j 7k
1 3 2
5 x – 2 – 3 y –1 7 z 3 0
or
5 x – 3 y 7 z 14 0
r t f t , g t , h t f t i g t j h t k
ln(1 t )
r (t ) (1 t )i sin t j k
t
ln 1 t
lim t 0 r t lim t 0 1 t i lim t 0 sin t j lim t 0 k
t
ik
x f t , y g t , z h t
dr r t h r t
r t lim h 0
dt h
Theorem 1.6: If r t f t , g t , h t f t i g t j h t k ,
where f, g, and h are differentiable functions, then:
r t f t , g t , h t f t i g t j h t k
Proof: We have
r t
1
lim t 0 [r t t r t ]
t
1
lim t 0 f t t , g t t , h t t f t , g t , h t
t
f t t f t g t t g t h t t h t
lim t 0 , ,
t t t
f t , g t , h t
Example 1.12:
1
r t 1 t i sin t j e 2 t k
2
Solution:
Scalar and vector fields 35
r t i cos t j e 2 t k
1
b) Since r 0 i k and r 0 i j k , the unit tangent vector
2
at the point 1, 0,1 is:
r 0 i j k 1 1 1
T(0) i j k
| r 0 | 111 3 3 3
d
1. u t v t u t v t
dt
d
2. cu t cu t
dt
d
3. f t u t f t u t f t u t
dt
d
4. u t v t u t v t u t v t
dt
d
5. u t v t u t v t u t v t
dt
d
6. u f t f t u f t (Chain Rule)
dt
b
r t dt
a
n
lim n r ti* t
i 1
n n n
lim n f ti* t i g ti* t j h ti* t k
i 1 i 1 i 1
Thus
a
b
r (t ) dt
a
b
f (t ) dt i
b
a h(t ) dt k
g (t ) dt j
b
F x, y P x , y i Q x, y j P x, y , Q x, y
F Pi Qj
Here P and Q are scalar functions of two variables. They are called
scalar fields to distinguish them from vector fields.
x F x xi yj – yi xj – xy yx 0
| F ( x, y ) | ( y ) 2 x 2 x 2 y 2 | x |
38 Mathematics for mechanical engineers
F x, y , z P x, y , z i Q x, y , z j R x , y , z k
F x, y , z zk
Another example is the wind field which is very important for weather
forecasting. Figure 1.33 presents the wind field in North Vietnam on
Monday 16 August 2021 at 4:50 PM. They indicate the wind speed
40 Mathematics for mechanical engineers
mMG
| F |
r2
r x
So,
2
r2 x
The gravitational force exerted on the airplane acts toward the origin.
The unit vector in this direction is:
x
|x|
mMG
F (x ) x
| x |3
Physicists often use the notation r instead of x for the position vector.
So, we can write:
Scalar and vector fields 41
mMG
F r (1.12)
r3
x xi yj zk and x x 2 y 2 z 2
F ( x, y , z )
mMGx mMGy mMGz
i 2 j 2 k
(x y z )
2 2 2 3/2
(x y z )
2 2 3/ 2
( x y 2 z 2 )3/ 2
0.5
-0.5
0.5
0.5
0
0
y -0.5 -0.5 x
mMG
f ( x, y , z )
x2 y2 z2
Then,
f ( x, y, z )
f f f
i j k
x y z
mMGx mMGy mMGz
i 2 j 2 k
(x y z )
2 2 2 3/2
(x y z )
2 2 3/ 2
( x y 2 z 2 )3/ 2
F ( x, y, z )
F ( x, y , z )
mMGx mMGy mMGz
i 2 j 2 k
(x y z )
2 2 2 3/2
(x y z )
2 2 3/ 2
( x y 2 z 2 )3/ 2
mMG
f ( x, y , z )
x2 y2 z2
so that
Scalar and vector fields 43
f ( x, y, z ) F( x, y, z ) .
EXCERCISES
I. Find the sum of the given vector and illustrate geometrically.
1. 0, 4 , 3, 0
2. 1, 3 , 3, 2
3. 1, 0, 2 , 2, 2, 0
4. 0, 0,3 , 1, 0, 0 .
1. a 2, 4 , b 4, 2
2. a i 3 j, b 2i 5 j
3. a i 3j 2k , b 2i 2 j k
4. a 3i 2k , b 2 j 3k
5. a i 3j k , b 2i k .
III. Find a unit vector that has the same direction as the given vector
1. 3, 4
2. 5i 3j
3. 3, 4, 5
4. 2i 3j 4k .
1. a 2, 4 , b 4, 2
2. a 1, 3 , b 3, 2
3. a i 3 j, b 2i 5 j
4. a 1, 0, 2 , b 2, 2, 0
5. a 0, 0,3 , b 1, 0, 0
6. a i 3j 2k , b 2i 2 j k
7. a 3i 2k , b 2 j 3k
8. a i 3j k , b 2i k .
1. a 2, 4 , b 4, 2
2. a 1, 3 , b 3, 2
3. a i 3 j, b 2i 5 j
4. a 1, 0, 2 , b 2, 2, 0
5. 0, 0,3 , 1, 0, 0
6. a i 3j 2k , b 2i 2 j k
7. a 3i 2k , b 2 j 3k
8. a i 3j k , b 2i k
1. a 0, 0,3 , b 1, 0, 0
2. a 1,3, 9 , b= 1, 3, 9
4. a 0, 0, 3 , b 2,1,5
5. a i 3j 2k , b 2i 2 j k
6. a 3i 2k , b 2 j 3k
7. a i 3j k , b 2i k .
1. The plane through the point 1, 2,3 and with normal vector
2, 0, 3 .
2. The plane through the point 3, 2,3 and with normal vector
1,3, 2 .
4. The plane through the point 3, 2,3 and perpendicular to the line
x t , y 3 2t , z 5 3t .
VIII. Find the distance from the point to the given plane
1. 3, 1, 2 , x 2 y 3z 0
2. 0, 3, 5 , 2x y 4z 2
Scalar and vector fields 47
3. 2, 1, 4 , 4 x y 3z 9
4. 3,1,9 , x 3 y 5z 9
1. x 2 y 3z 2, 2 x 4 y 6 z 5
2. P : 2 x y 5 z 1, Q :4 x 2 y 10 z 5 .
ANSWERS
I.
1. 0, 4 3, 0 3, 4
2. 1, 3 3, 2 2, 1
3. 1, 0, 2 2, 2, 0 1, 2, 2
4. 0, 0, 3 1, 0, 0 1, 0, 3
48 Mathematics for mechanical engineers
II.
2.
a i 3j, b 2i 5 j
a b 1, 3 2, 5 1, 2
2a 2 1, 3 2, 6 22 6 40 2 10
2
2a b 2 1, 3 2,5 2, 6 2,5
4.
a 3i 2k , b 2 j 3k
a b 3, 0, 2 0, 2,3 3, 2,1
2a 6,0, 4 62 02 4 52
2
62 2 7 89
2 2
III. Find a unit vector that has the same direction as the given vector
5 3
2. v 5i 3 j v 52 32 34 u ,
34 34
2
2
4. v 2i 3j 4k v 32 42 29
2 3 4
u i j k
29 29 29
4. a 1, 0, 2 , b 2, 2, 0 a.b 1 2 0 2 2 0 2
6. a i 3 j 2k , b 2i 2 j k
V.
a 10, b 13
50 Mathematics for mechanical engineers
a.b 9 9
cos a, b a, b 1420
a b 10 13 130
a 5, b 8
a.b 2 2
cos a, b a, b 1080 26'
a b 5 8 40
6. a i 3 j 2k , b 2i 2 j k
a.b 10
cos a, b a, b 10
a b 3 11
8. a i 3j k , b 2i k
a.b 1 1
cos a, b a, b 82015'
a b 11 5 55
VI.
6. a 3, 0, 2 , b 2, 0,1 a b 0, 7, 0
VII.
Scalar and vector fields 51
2 x 1 0 y 2 3 z 3 0 or 2 x 3 z 11
1 x 3 2 y 2 3 z 3 0
or x 2 y 3z 16
VIII.
1
The point 0, 0, is in the plane Q : 4 x 2 y 10 z 5 .
2
The distance between P and Q is the same as the distance from
1
0, 0, to P. This distance is
2
1
5. 1
2 7
D .
2 12 5 2 30
2 2
52 Mathematics for mechanical engineers
Gradient, divergence, curl and laplacian operators 53
f ( x0 h, y0 ) f ( x0 , y0 )
f x ( x0 , y0 ) lim
h0 h
f ( x0 , y0 h) f ( x0 , y0 )
f y ( x0 , y0 ) lim
h 0 h
f z
f x ( x, y ) f x f ( x, y ) f1 D1 f Dx f
x x x
f z
f y ( x, y ) f y f ( x, y ) f 2 D2 f Dy f
y y y
54 Mathematics for mechanical engineers
f x x, y 3 x 2 9 x 2 y 2
Thus,
f y x, y 6 x3 y 8 y 3
Then,
A x x0 B y y0 C z z0 0
z z0 a x x0 b y y0 (2.1)
Suppose that Equation 2.1 represents the tangent plane at P, then this
tangent plane will intersect the plane y y0 by the tangent line T1 .
Setting y y0 in Equation 2.1 gives:
z z0 a x x0
or
z ax ax0 z0
z z0 f x x0 , y0 x x0 f y x0 , y0 y y0 (2.2)
dz f dx f dy
ds x ds y ds
dz z dx z dy
(2 xy 3 y 4 )(2 cos 2 s) ( x 2 12 xy 3 )( sin s )
ds x ds y ds
dz
(0 3)(2.1) (0 0)( 0) 6
ds s 0
i j k
x y z
f f
f ( x, y) f x ( x, y ), f y ( x, y) i j
x x
f ( x, y ) f x , f y 2 xy , x 2 2 y cos y 2
f (1, 0) 0,1
f ( x0 ha, y0 hb) f ( x0 , y0 )
Du f ( x0 , y0 ) limh0
h
Theorem 2.2: With the notation for the gradient vector, we can express
the directional derivative as:
Du f ( x, y ) f x ( x, y )a f y ( x, y )b
f x ( x, y ), f y ( x, y ) a, b
f ( x , y ) u
g (h) g (0)
g (0) lim
h 0 h
f ( x0 ha, y0 hb) f ( x0 , y0 )
lim (2.3)
h 0 h
Du f ( x0 , y0 )
f x f y
g (h) f x ( x, y )a f y ( x, y )b (2.4)
x h y h
If we put h 0 into Equation 2.4 and compare with Equation 2.3, then
Du f ( x0 , y0 ) f x ( x0 , y0 ) a f y ( x0 , y0 )b
Du f ( x, y ) f x ( x, y ) a f y ( x, y )b
f ( x, y ) y 3 8 x i 3xy 2 j
f (1, 2) 16 i 12 j
Note that v is not a unit vector so we need to find the unit vector in the
direction of v. Since | v | 10 , the unit vector in the direction of v is:
v 1 3
u i j
|v| 10 10
Therefore,
1 3
Du f (1, 2) f (1, 2) u (16i 12 j) i j
10 10
16 1 12 3 52
10 10
Gradient, divergence, curl and laplacian operators 59
f ( x, y , z ) f x ( x, y , z ), f y ( x, y , z , ), f z ( x, y , z )
or,
f f f
f f x , f y , f z i j k
x y z
a) The gradient of f
Solution:
f ( x, y , z ) f x ( x, y , z ), f y ( x, y, z ), f z ( x, y, z )
ye xz xyze xz , xe xz , x 2 ye xz
b) At 1, 2, 0 , we have:
f (1, 2, 0) 2,1, 2
3 2 1
u i j k
14 14 14
60 Mathematics for mechanical engineers
Du f (1, 2, 0) f (1, 2, 0) u
3 2 1 6
2i j 2k i j k
14 14 14 14
r t x t , y t , z t
r t0 x0 , y0 , z0
F x t , y t , z t k
F dx F dy F dz
0 (2.5)
x dt y dt x dt
Gradient, divergence, curl and laplacian operators 61
F Fx , Fy , Fz
r(t ) x t , y t , z t
F r t 0
r t0 x0 , y0 , z0
So,
F x0 , y0 , z0 r t0 0
Fx ( x0 , y0 , z0 )( x x0 ) Fy ( x0 , y0 , z0 )( y y0 ) Fz ( x0 , y0 , z0 )( z z0 ) 0
x x0 y y0 z z0
Fx ( x0 , y0 , z0 ) Fy ( x0 , y0 , z0 ) Fz ( x0 , y0 , z0 )
R Q P R Q P
curl F i j k
y z z x x y
i j k
x y z
f f f f f f
f i j k i j k
x y z x y z
i j k
F
x y z
P Q R
R Q P R Q P
i j k
y z z x x y
curl F
Gradient, divergence, curl and laplacian operators 63
R Q P R Q P
i j k
y z z x x y
curl F
curl F F (2.6)
curl F F
i j k
x y z
x2 xyz 2 y2 z
y 2 z xyz 2 i y 2 z x 2 j
y z x z
xyz 2 x 2 k
x y
2 yz 2 xyz i 0 0 j yz 2 0 k
2 yz 1 x i yz 2 k
curl f 0
curl f f
i j k
x y z
f f f
x y z
2 f 2 f 2 f 2 f 2 f 2 f
i j k
y z z y z x x z x y y x
0i 0 j 0k 0
F x, y , z x 2 i xyz 2 j zy 2k
is not conservative.
curl F 2 yz 1 x i yz 2 k
Example 2.7:
Solution:
a) We have:
curl F F
i j k
x y z
2 xyz x2 z x2 y
x 2 x 2 i 2 xy 2 xy j 2 xz 2 xz k 0
f x x, y , z 2 xyz
f y x, y , z x 2 z (2.7)
f z x, y , z x 2 y
f x , y , z x 2 yz g y , z
66 Mathematics for mechanical engineers
f y x, y , z x 2 z g y y , z ,
Thus, comparing this equation with the second equation in (2.7) gives:
g y y, z 0
g y, z h z
Thus:
f x , y , z x 2 yz h z
f z x, y , z x 2 y
h z 0
Thus,
hz C
Therefore:
f x , y , z x 2 yz C (C – constant).
P Q R
div F
x y z
i j k
x y z
div F F
2
div F F
x
x xyz 2 zy 2 2 x xz 2 y 2
y z
div curl F 0
div curl F F
R Q P R Q P
x y z y z x z x y
2 R 2Q 2 P 2 R 2Q 2 P
0
x y x z y z y x z x z y
a a b 0
equation for heat and fluid flow, wave propagation, and quantum
mechanics.
2 f 2 f 2 f
div f f
x 2 y 2 z 2
2 f 2 f 2 f
2 f 0
x 2 y 2 z 2
F Pi Qj Rk
2 F 2 P i 2Q j 2 R k
EXERCISES
1) F y 2 i x 2 j z 2k
2) F x y 2 i xyj zk
3) F x 2 yzi xy 2 zj xyz 2k
4) F ln yi e xy j yxk
5) F e x i e y j e z k
1) f x, y , z 4 x 4 y 2 z
2) f x, y , z 3 x 2 2 y 2 2 z 2
3) f x, y, z x ln y 2 z 2
4) f x , y , z x e y z
5) f x, y , z 5 x 2 5 y 2 x y 4 z
ANSWERS
I.
2. F x y 2 i xyj zk
div F=.F=
x
x y 2 xy z 2 x
y z
Gradient, divergence, curl and laplacian operators 71
div F 1, 0, 0
i j k
curl F= F= 3 yk
x y z
x y2 xy z
2 F= 0
4. F ln yi e xy j yxk
div F xe xy
div F 1 xy e xy i x 2 e xy j
or div F e xy 1 xy , x 2 e xy , 0
1
curl F x, y, ye xy
y
2 F= x 2e xy
II.
2. f x, y , z 3 x 2 2 y 2 2 z 2
2 f 2 f 2 f
f 2 2 2 3 2 2 7 0
2
x y z
4. f x , y , z x e y z
72 Mathematics for mechanical engineers
2 f 2 f 2 f
2 f 0 e y z e y z 2e y z 0 .
x 2 y 2 z 2
5. f x, y , z 5 x 2 5 y 2 x y 4 z
2 f 2 f 2 f
2 f 10 10 0 0
x 2 y 2 z 2
R a , b c, d
x, y | a x b, c y d
Let S be the solid that lies above R and Figure 3.1. Surface
under the graph of f, that is, above a solid
S x, y , z 3 | 0 z f x, y , x, y
Rij xi –1 , xi y j –1 , y j x, y | xi –1 x xi , y j –1 y y j
A xy
In each rectangle Rij we choose a sample point xij* , yij* . Then, we can
approximate the part of S that lies above each Rij by a thin rectangular
column with base Rij and height f xij* , yij* . The volume of this
Multiple integrals 75
column is the height of the column times the area of the base
rectangle:
f xij* , yij* A
This double sum is called a double Riemann sum. We can see that as
the rectangular becomes smaller, i.e as m and n become larger the
approximation given in this equation becomes better. So, we would
expect that:
m n
V lim
m , n
f ( x , y ) A
i 1 j 1
*
ij
*
ij (3.1)
m n
f ( x, y ) dA lim
m , n
f ( x , y
i 1 j 1
*
ij
*
ij ) A
R
m n
f ( x, y ) dA lim
m , n
( x , y ) A
i 1 j 1
i i (3.2)
R
From Equations 3.1 and 3.2, the volume V of the solid that lies above
the rectangle R and below the surface z f x, y is:
V f ( x, y ) dA
R
4 4
V f ( xi , y j ) A
i 1 j 1
Figure 3.4 shows the approximating rectangular boxes and the graph
of the elliptic paraboloid.
3.1. [ f ( x, y ) g ( x, y )] dA f ( x, y ) dA g ( x, y ) dA
R R R
3.2. cf ( x, y) dA c f ( x, y) dA
R R
3.3. f ( x, y) dA g ( x, y) dA
R R
78 Mathematics for mechanical engineers
R = {(x, y) |a ≤ x ≤ b, c ≤ y ≤ d
then
b d d b
f ( x, y ) dA
R
a c
f ( x, y ) dy dx
c
a
f ( x, y ) dx dy
f is bounded on R.
f is discontinuous only on a finite number
of smooth curves.
The iterated integrals exist.
x 2 y dA , where
2
Example 3.2: Evaluate the double integral
R
R x, y | 0 x 1, 0 y 2
1 2
( x 2 y ) dA ( x 2 2 y ) dy dx
2
0 0
R
1 y 2
x 2 y y 2 dx
0 y 0
1
1 2 x3 10
(2 x 4) dx 2
4x
0
3 0 3
2 1
( x 2 y ) dA (x 2 y ) dx dy
2 2
0 0
R
x 1
2 x3
3 2 xy dy
0
x0
2
2 1 1 2 10
2 y dy y y
1
3 3 1 3
The results are the same whether we first integrate with respect to x or
y. It is noted that the result in this example is negative. However,
nothing is wrong with that because we do not calculate a volume, we
just calculate a given integral with non-positive function f.
f ( x, y ) if ( x, y ) is in D
F ( x, y )
0 if ( x, y ) is in R but not in D
f ( x, y) dA F ( x, y) dA
D R
y
The type I region is defined as a region which lies between the graphs
of two continuous functions of x, that is,
Multiple integrals 81
D x, y | a x b, g1 x y g 2 x
To compute f ( x, y ) dA
D
over a region D of type I, we choose a
b d
f ( x, y) dA F ( x, y) dA
D R
a c
F ( x, y ) dy dx
The explanation for this formula is addressed here. From Figure 3.8
we see that F x, y 0 if y g1 x or y g 2 x because x, y
then lies outside D. Thus,
d g2 ( x ) g2 ( x )
c
F ( x, y ) dy
g1 ( x )
F ( x, y ) dy
g1 ( x )
f ( x, y ) dy
because F f when
g1 x y g 2 x
then
82 Mathematics for mechanical engineers
b g2 ( x )
f ( x, y) dA
D
a g1 ( x )
f ( x, y ) dy dx
D x, y | c y d , h1 y x h2 y
y y
d d
𝑥 = ℎ (𝑥) D
D 𝑥 = ℎ (𝑥)
𝑥 = ℎ (𝑥) 𝑥 = ℎ (𝑥)
c c
O x O x
d h2 ( y )
D
f ( x, y ) dA
c
h1 ( y )
f ( x, y ) dx dy
D x, y | –2 x 3, x 2 y x 6
( x 2 y) dA
D
3 x6
( x 2 y ) dy dx
2 x 2
1
[ xy y 2 ] yy xx2 6 dx
1
1
( x 4 x3 6 x 36) dx
1
3
x5 x 4
3x 2 36 x
5 4 2
625
Figure 3.10. Type I region
4
Example 3.4: Find the volume of the solid that lies under the plane
x 2 y z 0 and above the region D in the xy–plane bounded by the
line y 2 x and the parabola y x 2 .
D x, y | 0 x 2, x 2 y 2 x
V ( x 2 y ) dA
D
2 2x
( x 2 y ) dy dx
0 x2
2 y 2 x
xy y 2 2 dx
0 yx
2
x5 x 4 28
x x 6 x dx 2 x3
2
4 3 2
0 5 4 0 5
1
D x, y | 0 y 4, y x y
2
Multiple integrals 85
V ( x 2 y ) dA
D
4 y
1
( x 2 y ) dx dy
0 2
y
x y
4 x2
2 2 xy dy
0
x 12 y
4 y y2
2 y 3/2 y 2 dy
0
2 8
4
3 4 1 28
y 3 y 5/2 y 2
8 5 4 0 5
D x, y | 0 x 1.5, x / 2 y 3 – x / 2 .
So,
V (3 x 2 y ) dA
D
3/ 2 3 x / 2
(3 x 2 y ) dy dx
0 x /2
3/ 2 y 3 x / 2
3 y xy y 2 y x / 2 dx
0
1 3/ 2
4 0
4 x 2 12 x 9 dx
3/ 2
1 4 x3 9
6 x 2 9 x
4 3 0 8
sin x 2 dx dy sin x 2 dA
2 2
0 y
D
where D x, y | 0 y 2, y x 2
D x, y | 0 x 2, 0 y x
sin x dA sin x 2 dy dx
2 x
2
0 0
D
y sin x 2 y 0 dx
2 yx
x sin x 2 dx
2
0
2
1 1
cos x 2 (1 cos 4)
2 0 2
3.4. f x, y g x, y dA f x, y dA g x, y dA
D D D
3.5 cf x, y dA c f x, y dA
D D
88 Mathematics for mechanical engineers
If f x, y g x, y for all x, y in y D
D, then:
3.6. f ( x, y) dA g ( x, y) dA
D D
D1
D2
If D D1 D2 , where D1 and D2 x
don’t overlap except perhaps on their
boundaries as shown in Figure 3.14, Figure 14. Union regions
then:
3.7. f x, y dA f x, y dA f x, y dA
D D1 D2
3.8. 1 dA A D
D
3.9. mA( D ) f x, y dA MA D
D
Property 3.9 can be easily proved from Properties 3.5, 3.6, and 8.
r 2 x 2 y 2 , x r cos , y r sin
The regions in the first figure are special cases of a polar rectangle as
illustrated in Figure 3.18.
90 Mathematics for mechanical engineers
R r , | 0 r a , 0 2 R r , | 0 , a r b
a0 ba0
Figure 3.17. Regions that can be described in polar coordinates
R r , | a r b,
Figure 3.18. Polar rectangles
By this way, the circles r ri and the rays i divide the polar
rectangle R into the small polar rectangles shown in Figure 3.19. Here,
the “center” of the polar subrectangle:
Rij ( r , ) | ri –1 r ri , j –1 i
ri* 1
2 ri –1 ri
*j 12 j –1 j
Ai 12 ri 2 12 ri 21
12 (ri 2 ri 21 )
12 (ri ri 1 )(ri ri 1 )
ri* r
The rectangular coordinates of the center of Rij are ri*cos *j , ri* sin *j .
So, a typical Riemann sum is:
m n
f (r
i 1 j 1
i
*
cos *j , ri* sin *j ) Ai
m n
f (ri* cos *j , ri* sin *j ) ri* r
i 1 j 1
92 Mathematics for mechanical engineers
m n
g (r ,
i 1 j 1
i
* *
j ) r
b
a
g (r , ) dr d
Therefore, we derive:
m n
f ( x, y ) dA lim
m , n
f (r
i 1 j 1
i
*
cos *j , ri* sin *j ) Ai
R
m n
lim
m , n
g (r ,
i 1 j 1
i
* *
j ) r
b
g ( r , ) dr d
a
b
f ( r cos , r sin ) r dr d
a
0 a r b,
where 0 – 2 , then
b
f ( x, y ) dA
R
a
f (r cos , r sin )r dr d (3.3)
R x, y | y 0, 1 x 2 y 2 4
1 r 2, 0
R r , |1 r 2, 0
( x 2 y ) dA
2
Figure 3.20. Half ring
R
2
(r 2 cos 2 2r sin ) r dr d
0 1
2
(r 3 cos2 2r 2 sin ) dr d
0 1
2
r4
2
cos 2 r 3 sin d
0
4 3 1
15 14
cos 2 sin d
0
4 3
15 14
1 cos 2 sin d
0
8 3
15 28
8 3
Example 3.8: Evaluate the volume of the solid bounded by: the plane
z 0 and the paraboloid z 4 – x 2 – y 2 as shown in Figure 3.21.
94 Mathematics for mechanical engineers
R r , | 0 r 2, 0 2
4 – x2 – y 2 4 – r 2
V (4 x 2 y 2 ) dA
D
2 1
(4 r
2
) r dr d
0 0
1
2 1 r4 7
d (4r r ) dr 2 2r 2
3
0 0
4 0 2
2 4 x2
V (4 x 2 y 2 ) dA (4 x 2 y 2 ) dy dx
2 4 x 2
D
D ( r , ) | , h1 ( ) r h2 ( )
then
h2 ( )
f ( x, y) dA
D
h1 ( )
f (r cos , r sin ) r dr d
Example 3.9: Use a double integral to find the area enclosed by a half
of the cardioid r 1– cos .
D r , | 0 , 0 r 1 – cos
Figure 3.23.
A half of the
cardioid
96 Mathematics for mechanical engineers
A( D) dA
D
1 cos
r dr d
0 0
[ 12 r 2 ]10cos d
0
1 cos
2
1
2 d
0
1
4 0
(3 cos cos 2 ) d 43
Figure 3.24 shows that the solid lies above the paraboloid z x 2 y 2
and under the paraboloid z 4 3x 2 3 y 2 . The projection of solid
onto xy-plane is the disk D x 2 y 2 1 . In polar coordinates, we have:
x2 y2 r 2
4 3x 2 3 y 2 4 3r 2
D {(r , ) | 0 2 , 0 r 1}
So, we have:
V 4 (1 x 2 y 2 ) dA
D
1 r r dr d
2 1
4 2
0 0
1
r2 r4
8 2
2 4 0
B x, y , z a x b, c y d , r z s
Bijk xi 1 , xi y j 1 , y j zk 1 , z k
l m n
f x
i 1 j 1 k 1
*
ijk
*
, yijk *
, zijk V
Bij
z
x x
y
l m n
f x, y, z dV lim
l , m, n
f x
i 1 j 1 k 1
*
ijk
*
, yijk *
, zijk V
B
We can choose the sample point to be any point in the sub-box. For
simplicity reason, we choose it to be the point xi , y j , zk that we get a
simpler-looking expression:
l m n
f x, y, z dV lim
l , m , n
f x , y , z V
i 1 j 1 k 1
i j k
B
s d b
f x, y, z dV f x, y, z dx dy dz
B
r c a
b s d
f x, y, z dV f x, y, z dy dz dx
B
a r c
100 Mathematics for mechanical engineers
B x, y , z 0 x 1, 1 y 2, 0 z 1
1 2 1
xyz dV xy 1 z dx dy dz
2
0 1 0
B
x 1
1 x 2 y 1 z
2
dy dz
0 1
2 x 0
1 1 2
y 1 z dy dz
2 0 1
y 2
1 1 3 1
y 2 1 z dz 1 z dz
4 0 y 1 4 0
1
3 z 2 3
z
4 2 0 8
f x, y, z dV F x, y, z dV
E B
The solid region of type I it the region which lies between the graphs
of two continuous functions of x and y. That is,
E x, y, z x, y D, u x, y z u
1 2 x, y
We notice from Figure 3.27 that the upper boundary of the solid E is
the surface with equation z u2 x, y and the lower boundary is the
surface z u1 x, y . Using the same procedures to derive the formula
for double integrals in region I, it can be shown that, if E is a type I
region, then:
u2 x , y
f x, y, z dz dA
f x, y, z dV
E D
u1 x , y
The meaning of the inner integral on the right side of Equation 3.6 is
that x and y are held fixed. Therefore, u1 x, y and u2 x, y are
regarded as constants. f x, y , z is integrated with respect to z. In
particular, if the projection D of E onto the xy-plane is a type I plane
region, then
102 Mathematics for mechanical engineers
E x, y , z a x b, g1 ( x ) y g 2 ( x ), u1 ( x, y ) z u2 ( x, y )
b g2 ( x ) u2 ( x , y )
f x, y, z dV
E
a g1 ( x )
u1 ( x , y )
f x, y , z dz dy dx
D
x
E x, y , z c y d , h1 ( y ) x h2 ( y ), u1 ( x, y ) z u2 ( x, y )
So,
d h2 ( y ) u2 ( x , y )
f x, y, z dV
E
c h1 ( y )
u1 ( x , y )
f x, y , z dz dx dy
E x, y, z | 0 x 3,
3 x
0 y ,
2
0 z 3 x 2 y
3 x
3 3 x 2 y
x dV
E
0 0
2
0
x dz dy dx
3 x
xz
3 z 3 x 2 y
2
z 0
dy dx
0 0
3 x
3
2
(3 x x 2 2 xy )dy dx
0 0
3 x
3 y
3 xy x 2 y xy 2
2
dx
0 y 0
Figure 3.29. The given
1 3 3
4 0
x 6 x 2 9 x dx solid tetrahedron
3
1 x4 9 27
2 x3 x 2
4 4 2 0 16
E x, y, z y, z D, u ( y, z ) x u ( y, z )
1 2
f x, y, z dV f x, y, z dx dA
u2 ( y , z )
(3.4)
u1 ( y , z )
E D
x E
y
E x, y, z x, z D, u ( x, z) y u
1 2 x, z
where:
D
y
f x, y, z dV f x, y, z dy dA
u 2( x , z )
(3.5)
u1 ( x , z )
E D
Notice that there might be two possible expressions for the integral
depending on whether D is a type I or type II plane region using each
of Equations 3.4 and 3.5.
E x, y, z 0 y 2, y x y, y 2 x2 z y2 x2
Therefore, we have:
x z 2 dz dx dy
2 y y 2 x2
x 2 y 2 dV
2
0 y y 2 x2
E
y 2 dV x 2 z 2 dy dA
2
x
2
E
D2
x2 z 2
D2
2 x 2 z 2 x z dA
2 2
2 4 x2
2 4 x 2 2 x2 z 2 x z dz dx
2 2
x r cos , z r sin
That gives:
x
E
2
D2
z 2 dV 2 x 2 z 2 x z dA
2 2
2 2
2 r r r dr d
2
0 0
d 2r r dr
2 2
3 4
0 0
2
r 4 r 5 16
2
2 5 0 5
V E dV
E
1 dV dz dA u 2 ( x, y ) u1 ( x, y ) dA
u2 ( x , y )
u1 ( x , y )
E D D
V T dV
T
3 3 x / 2 3 x 2 y
dz dy dx
0 0 0
3 3 x / 2 9
3 x 2 y dy dx
0 0 4
m x, y, z dV
E
The moments of the object about the three coordinate planes can be
expressed as:
110 Mathematics for mechanical engineers
M yz x x, y, z dV
E
M xz y x, y, z dV
E
M xy z x, y, z dV
E
M yz M xz M
x , y , z xy
m m m
If the density is constant, the center of mass of the solid is called the
centroid of E. The moments of inertia about the three coordinate axes
are:
I x y 2 z 2 x, y , z dV
E
I y x 2 z 2 x, y , z dV
E
I z x 2 y 2 x, y , z dV
E
x r cos z
y r sin
zz
z
To convert from rectangular to
cylindrical coordinates the
y
following equation are used: r
r 2 x2 y2 x
tan y / x
Figure 3.34. Cylindrical coordinates
zz
D
x
E x, y , z | x, y D, u1 z u2 x, y
D r , | , h1 r h2
f ( x, y, z ) dV f x, y, z dz dA
u2 ( x , y )
E D
u1 ( x , y )
f x, y, z dV
E
h2 ( ) u2 r cos , r sin
f r cos , r sin , z r dz dr d
h1 ( ) u1 r cos , r sin
Example 3.15: Figure 3.36 presents a solid which lies within: The
cylinder x 2 y 2 1, below the plane z 2 , above the cone
z x 2 y 2 . The density at any point is proportional to its distance
from the axis of the cylinder with a proportionality constant K. Find
the mass of E.
E r , , z | 0 2 , 0 r 1, r z 2
f x, y , z x 2 y 2 r
f x, y, z Kr
m x, y, z dV
E
Thus,
m K x 2 y 2 dV
E Figure 36. The given solid
2 1 2
( Kr ) r dz dr d
0 0 r
2 1
Kr 2 r dr d
2
0 0
K 2 2r r dr
1
2 3
0
1
2 r 4 5 K
2 K r 3
3 4 0 6
2 4 x2 4
Example 3.16: Evaluate
2 4 x 2
x2 y 2
x 2 dz dy dx
E
{ x, y, z | 2 x 2, 4 x 2 y 4 x 2 , x 2 y 2 z 4}
This describes the solid region E in which the projection of E onto the
xy-plane is the disk x 2 y 2 4 , the upper surface of E is the
paraboloid z x 2 y 2 and the lower surface is the plane z 4 . Using
the cylindrical coordinates, we have:
E (r , , z ) | 0 2 , 0 r 2, r 2 z 4
2 4 x2 4
2 4 x 2 x2 y 2
x 2 dz dy dx x 2 dV
E
2 2 4
r cos 2 r dz dr d
0 0 r2
2 2 4 2 2 4
2
cos .r 2 dz dr d 2 r dz dr d
0 0 r 0 0 r2
2
2 r4
d r 4 r dr 4 2r 2 16
2
2 2
0 0
4 0
0
0 z
x sin cos
Figure 3.38. Spherical coordinate
y sin sin
z cos
2 x2 y2 z 2
f x, y, z dV
E
d b
f sin cos , sin sin , cos sin d d d
2
c a
E , , a b, , c d
Figure 3.39.
B , , |1 2, 0 , 0
2 2
x2 y 2 z 2 2
2
y dV 2 sin 2 sin 2 2 sin d d d
2 2 2
0 0 1
B
2
2 sin 3 d 2 sin 2 d 4 d
0 0 1
Multiple integrals 117
2
1 sin 2 2 5
cos cos
3 2
0 2
2 0 5 1
2 31 31
34 5 30
This gives: sin cos or .
4
E , , 0 2 , , 0 2
4 2
V ( E ) dV
E
2 /2 2
2 sin d d d
0 /4 0
118 Mathematics for mechanical engineers
2
2 /2 3
d sin d
0 /4
3 0
16 / 2
3 /4
sin d
16
cos /4
/2
3
8 2
3
b d
a
f ( x) dx f ( g (u)) g '(u) du
c
or
b d dx
a
f ( x) dx f ( x(u ))
c du
du
where x g u and a g c , b g d .
x r cos , y r sin
Multiple integrals 119
f ( x, y) dA f (r cos , r sin ) r dr d
R S
T u , v x, y
x x u, v , y y u, v
v y
T (x1, y1)
S
(u1, v1)
R
T-1
u x
u G x, y , v H x, y
u v 2u v
x , y
3 3
v 3, 0 u 4
u 3 2u 3
x , y .
3 3
7 5
which is the line segment from (1, –1) to , in the xy-plane.
3 3
u 4, 3 v 3
4v 8v
x , y
3 3
Eliminating v, we obtain:
y x 4, 1/ 3 x 7 / 3
7 5 1 11
which is the line segment from , to , in the xy-plane.
3 3 3 3
v 3, 0 u 4
u 0, –3 v 3
x0 , y0 T u0 , v0
ru gu (u0 , v0 ) i hu (u0 , v0 ) j
x y
i j
u u Figure 3.45. Image R of S
rv gv (u0 , v0 ) i hv (u0 , v0 ) j
x y
i j
v v
a r (u0 u , v0 ) r (u0 , v0 )
b r (u0 , v0 v) r (u0 , v0 )
Remember that,
r (u0 u , v0 ) r (u0 , v0 )
ru lim
u 0 u
So,
i j k x y x x
x y u u u v
ru rv 0 k k
u u x y y y
x y v u u v
0
v u
x x
( x, y ) u v x y x y
(u, v) y y u v v u
u v
( x, y )
A u v (3.7)
(u , v )
( x, y )
Rij Sij
(u , v)
f ( x, y ) dA
R
m n
f ( xi , y j ) A
i 1 j 1
m n
( x, y )
f ( g (ui , v j ), h(ui , v j )) u v
i 1 j 1 (u, v)
( x, y )
f ( g (u, v), h(u, v)) (u, v) du dv
S
Then,
( x, y )
f ( x, y ) dA f ( x(u, v), y(u, v)) (u, v) du dv
R S
( x, y )
dA du dv
(u , v)
b d dx
a
f ( x) dx f ( x(u ))
c du
du
We notice the similarity between the single and the double integrals:
dx
Instead of the derivative , we have the absolute value of the
du
( x, y )
Jacobian .
(u , v)
x g r , r cos
y h r , r sin
Multiple integrals 127
x x
( x, y ) r
(r , ) y x
r
cos r sin
r cos 2 r sin 2 r 0
sin r cos
( x, y )
f ( x, y) dx dy f (r cos , r sin ) (u, v) dr d
R S
b
f (r cos , r sin ) r dr d
a
This is the same as the formula for double integral in polar coordinate.
The x-axis.
The parabolic y 2 4 x 8 and the line segment
y 3x 27 , y 0 .
Solution:
Figure 3.48 shows that R is not an easy region for evaluating the
integral. Thus, we need to change variables to evaluate the given
integral. We knew that T S R where S is the rectangular
0, 3 0, 2 . So, we see that S is a much simpler region than R. To
derive the formula of the double integral uv-plane we compute
Jacobian of the variable change:
128 Mathematics for mechanical engineers
x x
( x, y ) u v 2u 1
2u 2 v 0, v 0
(u , v ) y y v u
u v
( x, y )
y dA uv (u, v) dA
R S
2 3
(uv)(2u 2 v) du dv
0 0
2 3
(2u 3v uv 2 ) du dv
0 0
u 3
1
2 1
u 4 v u 2v 2 dv
0 2 2
u 0
2 81 9
v v 2 dv
0
2 2
2
81 3
v 2 v 3 93
4 2 0
This example show that the derived integral was not very difficult to
solve as we were given a suitable change of variables. However, if we
are not supplied with a transformation, we need to think of an
appropriate change of variables.
Multiple integrals 129
u x y, v – 2 x y
uv 2u v
x , y
3 3
x x 1 1
( x, y ) u v 3 3 1
(u , v ) y y 2 1 3
u v 3 3
S u , v | –3 v 3, 0 u 4
( x, y )
x y dA u
R S
(u, v)
du dv
1 3 4
3 3 0
u du dv
1 3 4
dv udu
3 3 0
4
1 u2
3 3 16
3 2 0
x g u , v , w , y h u , v, w , z k u , v , w
x x x
u v w
( x, y , z ) y y y
(u, v, w) u v w
z z z
u v w
f ( x, y, z ) dV
R
( x, y , z )
f ( x (u , v, w), y (u, v, w), z (u, v, w)) du dv dw
S
(u , v , w )
x sin cos
y sin sin
z cos
( x, y, z )
2 sin 2 sin
( , , )
f ( x, y, z) dV
R
Obviously, we can see that this is the same with the formula for triple
integrals in spherical coordinates as presented in previous.
Multiple integrals 133
EXERCISE
I. Calculate the iterated integral:
2 1
1. 1 2xy dxdy
1 0
2x xy dxdy
1 3
3
2.
1 0
x y 2 dxdy
2 2
2
3.
0 0
1
4.
0 0
x sin ydxdy
ln 2 ln 5
5.
0 0
e 2 x y dxdy
5 x y 2 2 y 3 dA , R x, y | 0 x 2, 0 y 2
2
1.
R
2. sin x y dA , R x, y | 0 x , 0 y
R
xy 2
3. R x2 1 dA , R x, y | 0 x 1, 0 y 5
1 2x2
4. R 2 y 2 dA , R x, y | 0 x 2, 0 y 2
xye xy dA , R x, y | 0 x 2, 0 y 1
2
5.
R
6. x 2 y 3 dA , R x, y | 0 x 2, 0 y x
R
2y
7. dA , R x, y | 0 x 2 y , 0 y 1
R x 2 3
134 Mathematics for mechanical engineers
8) e 2 x / y dA , R x, y | y x y 3 , 0 y 1
R
4 x
1. f x, y dydx
0 0
1 x
2. f x, y dydx
0 0
1 x
3. f x, y dydx
0 x2
1 x2
4. f x, y dydx
0 x4
1. E
xdV , E x, y, z |0 y 2, 0 x
1 y2 , 0 z y
2. xydV , E x, y , z |0 x 1, 0 y x, x z 3 x
E
3. yz cos x 3 dV , E x, y , z |0 x 1, 0 y x, 0 z y 2
E
Multiple integrals 135
4. E
xy 2 dV , where E lies under the plane z 2 x y and above
the region in the xy-plane bounded by the curves y x , y 0 and
x 1.
VI. Use a triple integral to find the volume of the given solid E:
1. Find the volume of the solid E that lies above the cone
z x 2 y 2 and below the sphere x 2 y 2 z 2 4 .
2. Evaluate E
zdV , where E lies above the paraboloid z x 2 y 2
and below the plane z 2 y .
2 4 y2 4 x2 y 2
1.
2 4 y 2 0
z x 2 y 2 z 2 dzdxdy
x y 2 z 2 dzdxdy
1 9 y 2 16 x 2 y 2
2
2.
0 9 y2 x2 y 2
ANSWERS
I.
2. 81
136 Mathematics for mechanical engineers
1
4. 2 sin 2
2
II.
2. sin x y dA 0
R
1 2 x2
4.
11 2
R 2 y 2 dA 3 arctan 2
32
6. R
x 2 y 3 dA
7
1
8. R
e2 x / y dA
8
1 e2
III.
2. V 8
4. V 8
IV.
1 1
2. f x, y dxdy
0 y
1 4 y
4. f x, y dxdy
0 y
V.
1
2. E
xydV
2
1 1 2 x y 989
4. E
xy 2 dV
0
y2 0
xy 2 dzdxdy
3024
Multiple integrals 137
2y
zdV 2 zdzdA , where D is the disk x 2 y 1 1
2
2.
E x y2
D
E x, y, z |
2 y 2, 4 y 2 x 4 y 2 , 0 z 4 x 2 y 2
This describes the upper half of a sphere with radius of 2. Thus the
domain in the spherical coordinates will be:
E , , |0 2, 0 , 0 2 .
2
138 Mathematics for mechanical engineers
Differential forms 139
4.1 1-forms
F x, y, z dx G x, y, z dy
f f
df dx dy
x y
F x, y , z dx G x, y , z dy H x, y , z dz
f f f
df dx dy dz (4.1)
x y z
Where i, j, k are the standard unit vectors along the x, y, z axes. Under
this set up, the gradient f corresponds to df. The differential notation
is very suggestive and ultimately quite powerful. We can see that we
can write the previous concepts using by using differential notations.
For example, we can see how to derive Equation 4.1 just by using
differential notations. Suppose that that variables x, y, z depend on
some parameter t, and f depends on x, y, z, then the chain rule says:
f f x f y f z
t x t y t z t
Thus, the formula for df can be obtained by canceling the notation dt:
f f f
df dx dy dz
x y z
4.2 Exactness in 2
However, most differential forms are not exact. To see why, we notice
that Equation 4.2 is equivalent to:
f f
F ,G
x y
F 2 f 2 f G
(4.3)
y yx xy x
Differential forms 141
But Equation 4.3 would fail for most examples such as x3 ydx e y dy .
Indeed, we can see that:
F x3 y , G e y
F G
x3 , 0
y x
Thus:
F G
y x
F G
Definition 4.3: We call a differential Fdx Gdy closed if and
y x
are equal.
dy F x, y
dx G x, y
Fdx Gdy 0
f x df C
df Fdx Gdy
Line Integrals
b dx dy
Fdx Gdy F x t , y t G x t , y t dt
C a
dt dt
Differential forms 143
Lemma 4.1:
C
Fdx Gdy Fdx Gdy
C
C
Fdx Gdy Hdz
dx dy
b
a
F x t , y t , z t G x t , y t , z t
dt dt
dz
H x t , y t , z t dt
dt
C1
C2
This can be easily proved using the definition of the line integral.
Indeed, if df and C1 : a, b 2 then:
b df
C1
df
a dt
dt
df
dt f C1 b f C1 a
b
a dt
A similar setting up for the line integral over C2 gives same answer:
df
a
b
df dt f C 2 b f C 2
C2 a dt
If the C is closed as shown in Figure 4.1, which means that the starting
point coincides the endpoint, then we have:
b
C1
df df
C2 C1
f C 2 b f C 2 a
f C1 b f C 1 a 0 a C2
f x, y Fdx Gdy
C
x y
f x, y F t , y0 dt G x, t dt
x0 y0
Differential forms 145
df Fdx Gdy
or
f f
F, G
x y
f y
G x, y dy G x, y
y y 0
f x y
F t , y0 dt G x, t dt
x x 0 x x y0
y G x, t
F x , y0 dt
y0 x
G x, y F x, y F x, t
x y t
Thus,
f y F x, t
F x, y0 dt
x y0 t
F x, y0 F x, y F x, y0 F x, y
So,
f f
df dx dy F x, y dx G x, y dy
x y
146 Mathematics for mechanical engineers
F d F3z
However, consider the case that a ball is kicked into the sky. Under
the gravitation field, the ball will finally fall toward the ground. If the
trajectory of the ball is not straight, but it is a parabol C as shown in
the right picture of Figure 4.3. The gravitation field acting on the ball
presented by a vector F F1 , F2 , F3 is not constant, then the work
done by the gravitational field would be the integral line:
C
F1dx F2 dy F3 dz
or
Differential forms 147
C
F ds
From this example, we can see that the differences between the
differential forms and vectors: To understand the nature of physics,
equations are expressed by coordinates, while the general principles
behind are often differential forms. The differential form represents a
coordinate free way of dealing with differentiation.
mMGt mMG 1 t 2
2t dt
1
C Fx dx Fz dz 1 2
t 1 t t 1 t
3/ 2 3/ 2
2 2 2
mMGt 2tmMG 1 t 2 dt
1
1
1
mMGt mMG 2t 3 dt
1
1 1
mMGt 2 mMG 2t 4
0
2 1
4 1
The work done by the gravitational field on the ball is zero can be
explained by as follows. We can see that the trajectory is symmetric
with respect to the z-axis. The ball moves up in the first half trajectory
so the work done by the gravitational field is negative. While, the ball
moves down in the second half trajectory so the work done is positive.
Thus, the works done in the first half of trajectory is equal to that of
the second trajectory but with opposite sign, so that the total works in
both half of trajectory is zero.
148 Mathematics for mechanical engineers
4.3 2-form
1 T
u v
2
u v vT u
In 3 we have:
0 ae bd af cd
1
a, b, c d , e, f ae bd 0 bf ce
2
af cd bf ce 0
We can see that the nonzero entries are basically the same as for the
cross product:
a, b, c d , e, f bf ce i af cd j ae bd k
1 T 1
u v u v vT u vT uTT uT vTT
T T
2 2
1 T 1
v u uT v uT v vT u u v
2 2
u v v u
uu 0
c u v cu v u cv
u v w u v u w
These rules are very useful to work with wedge products rather than
the Definition 1.
Now we will use the wedge product to define the 2-form which gives
a convenient way to deal with essential geometric features of the cross
product and surface integrals. Suppose F1dx G1dy H1dz and
F2 dx G2 dy H 2 dz are two 1-forms. We express the wedge product
of these two 1-forms with the use of wedge product properties as
follows:
F1G2 dx dy G1 H 2 dy dz H1 F2 dz dx u u 0
G1 F2 dx dy H1G2 dy dz F1 H 2 dz dx u v v u
F1G2 G1 F dx dy G1 H 2 H1G2 dy dz
H1 F2 F1 H 2 dz dx
F x, y, z dx dy G x, y, z dy dz H x, y, z dz dx
F x, y , z dx dy G x, y , z dy dz H x, y , z dz dx
x dx ydy dx 2dy
3
Fy Gx ; Fz H x ; Gz H y
F fx ; G f y ; H fz
any piecewise C 1
curve connecting 0, 0, 0 to x0 , y0 , z0 , then
df .
F1dy dz F2 dz dx F3dx dy
F1 j k F2k i F3i j=F1i F2 j F3k
F x, y , z dx G x, y , z dy H x, y , z dz
d d d
d f df fd
d dx d dy d dz 0
Recall that:
df f x dx f y dy f z dz
So
G x Fy k H y Gz i Fz H x j= V
u v w v u w v w u ...
Notice that:
u vw 0
u1 u 2 v w u1 v w u 2 v w
d dx dy d dy dz d dz dx 0
dx dy dz dy dx dz dy dz dx ...
Fx dx dy dz Fy dy dy dz Fz dz dy dz
Gx dx dz dx Gy dy dz dx Gz dz dz dx
H x dx dx dy H y dy dx dy H z dz dx dy
Fx dx dy dz Gy dx dy dz H z dx dy dz
Fx Gy H z dx dy dz
x f u, v
y g u, v
z h u, v
u, v D
x x
dx u v
u v
y y
dy u v
u v
Thus,
x x y y
dx dy u v u v
u v u v
x y x y
du dv
u v v u
x, y
du dv
u, v
Differential forms 155
x, y
Notice that is Jacobian of the variable change from the xy-
u, v
coordinates to uv-coordinates. Similarly, we have the following
equations:
y z y z
dy dz du dv
u v v u
y, z
du dv
u, v
z x z x
dz dx du dv
u v v u
z, x
du dv
u, v
f x, y, z dS f x, y, z T T
S D u v dudv
x y z
Tu i j k
u u u
x y z
Tv i j k
v v v
Thus,
x y z x y z
Tu Tv i j k i j k
u u u v v v
i j k
x y z
u u u
x y z
v v v
y z z y z x x z x y y x
i j k
u v u v u v u v u v u v
y , z z , x x, y
= i j k
u, v u, v u, v
So
2 2 2
y , z z , x x, y
Tu Tv
u, v u, v u, v
f x, y, z dS
S
2 2
y , z z , x x, y
2 (4.5)
f x, y, z
dudv
u, v u, v u, v
D
Differential forms 157
S is the surface z x y 2 , 0 x 1,
0 y 2 as shown in Figure 4.5.
Solution: We regard x, y as
parameters, then:
y y
y, z x y z y Figure 4.5. Graph of S
1
x, y z z x y
x y
z z
z , x x y x z
2 y
x, y x x x y
x y
x x
x, y x y x y
1
x, y y y x y
x y
y dS y 1 2 y 1
2 2 2
dxdy
S D
1 2
y 1 1 4 y 2 dy dx
0 0
1 2
dx 2 y 1 2 y 2 dy
0 0
2
2 14 23 (1 2 y 2 )3/ 2
0
13 2
3
158 Mathematics for mechanical engineers
F dS
S S
F1dy dz F2 dz dx F3 dx dy
F dS F ndS
S S
Tu Tv
n
Tu Tv
F x, y , z F1 x, y , z i F2 x, y, z j F3 x, y , z k
then,
Differential forms 159
x, y y, z z, x
S F dS D u, v u, v u, v dudv
F1 F2 F3
x, y y , z z , x
F dudv
u, v u , v u, v
D
F Tu Tv dudv
D
Tu Tv
F Tu Tv dudv
D Tu Tv
F ndS
S
For a closed oriented surface S, its orientation is defined such that the
normal vector points outward. That is, the positive orientation of S
coincides the direction of the outward normal vector.
x y 2 z 2 , y y, z z
F xi zk xdy dz zdx dy
y, z x, y
S1 F d S S1 y, z y, z dydz
x z (4.4)
We have:
y y
y, z y z y z
1
y, z z z y z
y z
z z
z , x y z z x
2 y
y, z x x z z
y z
x x
x, y y z x y
2 z
y, z y y z y
y z
y , z z , x x, y
Tu Tv = i j k
u, v u, v u, v
= 1, 2 y, 2 z
F dS F n dS
S1 S1
y S1
2
z 2 2 z 2 dydz
y 2
3 z 2 dydz
S1
cos 3sin r dr
2 1
2 2 3
0 0
2 1
0
2 cos 2 d 0 r 3dr
y, z x, y
1, 0
y, z y, z
Therefore,
S2
F dS xdydz 0dydz 0
S2 S2
Finally,
F dS
S S1
F dS F dS
S2
162 Mathematics for mechanical engineers
EXERCISE
1. 3 ydx 2 xdy
2. 3 ydx 3xdy
3. ln ydx x 2 dy
4. e y dx e x dy
III Calculate:
1. C
dx dy , where C is a circle r sin , 0
2. C
ydx xdy , where C is a circle 0 r R, 0 2
IV. Let xdx ydy zdz, xdx zdy ydz, xydz . Calculate:
1.
2.
3.
4. d
Differential forms 163
5.
D
where D is the square 0 x 1, 0 y 1, z 1 oriented
upward normal.
6.
D
where D is the square 0 x 1, 0 y 1, z 1 oriented
upward normal.
ANSWERS
I.
3 y 3x
2. df 3 ydx 3xdy is exact 1-form since
y x
II.
C
xdy xdy xdy xdy xdy
S1 S2 S3 S4
2 2
xdy r cos d r sin r 2 cos 2 d
C 0 0
4
2 r 2 1 cos 2 r 2 1 sin u
d 2 r2
0 4 4 0
164 Mathematics for mechanical engineers
IV.
1.
xz yx dx dy y 2 z 2 dy dz xy zx dz dx
D D
x, y y, z z, x
xz yx y2 z2 xy zx dxdy
D
x, y x, y x, y
where:
x x y y
x, y x y x y y, z x y
1, 0
x, y y y x y x, y z z
x y x y
z z
z , x x y
0
x, y x x
x y
Thus,
1 1
xz yx dxdy x yx dydx
D D 0 0
1
1 xy 2 1x 1
x dx dx
0
2 0 0 2 4
Line integrals and circulation of a field 165
x x t , y y t , a t b (5.1)
y
P0
Pi
P1 C
Pi-1 Pn
P2
x
O
a ti-1 ti b t
r t x t i y t j
xi x ti , yi y ti
n
Notice that the sum f xi* , yi* si is similar to a Riemann sum. We
i 1
have known that if f is a continuous function, then the limit in
Definition 5.1 always exists. We have also known that the length of C
is:
2 2
L ds
b dx dy dt
C a dt dt
2 2
dx dy
f x, y ds f x t , y t dt
b
C a dt dt
(5.3)
Line integrals and circulation of a field 167
y
Example 5.1: Evaluate
C
4 ds
x
where C is the arc of the parabola
x2 1
y from 1, to (2, 2) as
2 2
shown in Figure 5.3.
t2
x t, y , 1 t 2
2
2
y 2 1t
C x 1 t 2 1 t dt 1 2t 1 t dt 1 1 t d 1 t
2 2
2 1/ 2
4 ds 4 2 2 2
2 2
3
3/ 2 2
1 3
1 t 2 5 5 2 2
168 Mathematics for mechanical engineers
C
f x, y ds f x, y ds f x, y ds ... f x, y ds
C1 C2 Cn
y
Example 5.2: Evaluate C
4 ds , where C consists of the arc
x
x2
C1 of the parabola y followed
2
by the vertical line segment C2 from
2, 0 to 2, 2 .
y 2
C1
4 ds 5 5 2 2
x 3
On C2 , x is fixed at 2 so we choose y as the parameter. Thus, the
equations of C2 are:
x 2, y t , 0 t 2
Line integrals and circulation of a field 169
and
2 2
y 2 t dx dt 2
C2 4 x ds 0 4 2 dt dt dt 0 2 tdt 4
Thus,
y y y 2
C
4 ds 4 ds 4 ds 5 5 2 2 4
x C1 x C 2 x 3
5.1.4 Mass
Consider a metal wire with shape like
a curve C that has non-uniform cross
section as shown in Figure 5.5. That is,
the linear density x, y of the wire
varies over the length of the ring.
Then, the mass of the segment of the
wire from Pi 1 to Pi in this figure can
be approximated by xi* , yi* si . So,
the total mass of the wire is
approximately xi* , yi* si . By Figure 5.5. Metal wire
increasing the number of segments of
the wire, we obtain the mass m of the metal wire as the limiting value
of these approximations:
n
m lim xi* , yi* si x, y ds
n C
i 1
y
Recall function f x, y 4
in Example 5.2. If this function
x
describes the density of a wire, then the integral in Example 5.2 would
represent the mass of the wire. We now define that the center of mass
of the wire with density function ρ is the point x , y , where:
170 Mathematics for mechanical engineers
1
m C
x x x, y ds
(5.4)
1
y y x, y ds
m C
Example 5.3: Consider a wire that it has the shape of a line segment
y x, and its cross section increases from the base 0, 0 to the top
2, 2 as presented in Figure 5.6. Suppose that the linear density at
any point is proportional to its distance from the line y 0 by the
proportional factor k. Find the center of mass of the wire.
x t , y t , ds 2dt , 0 t 2
x, y ky
m ky ds
C
2
2ktdt
0
2
t2 Figure 5.6. Line segment with
2k 2 2k
2 0 non-uniform cross section
1 1
x
m C
x x, y ds xkyds
m C
2
1 2 1 t3 4
2 2k
0 ktt 2 dt
2 3 0 3
Line integrals and circulation of a field 171
and,
1
m C
y y x, y ds
1
2 2k C
y ky ds
2
1 2 1 t3 4
2t dt
2
2 2 0 2 3 0 3
4 4
Thus, the center of mass is: ,
Figure 5.7. Centre of mass
3 3
as shown in Figure 5.7. In special case, where si are replaced by
either:
xi xi xi 1
yi yi yi 1
Then, we have two integrals called the line integrals of f along C with
respect to x and y:
C
P x, y dx Q x, y dy P x, y dx Q x, y dy
C C
(5.6)
This formula is exactly the same with the line integral using
differential form as presented in Definition 4.4 in Chapter 4. To
evaluate this line integral, we can parameterize the curve whose
geometric description is given as presented in previous. However, if
we need to evaluate the line integral over a line segment, it is useful to
172 Mathematics for mechanical engineers
r t 1 t r0 tr1 0 t 1 (5.7)
1
r t 1 t 1, t 2, 2 , 0 t 1
2
Then,
x t 1 t 1 2t 3t 1
1 3 1
y t 1 t 2t t
2 2 2
Line integrals and circulation of a field 173
3
and dx 3dt , dy dt . Thus, the parametric representation for the
2
line segment is:
3 1
x 3t 1, y t , 0 t 1
2 2
1 3 1 3
C1 ydx x dy 0 2 t 2 3 dt 2 3t 1 2 dt
3 1
9t 1 dt
2 0
1
3 9t 2 21
t
2 2 0 4
x2
x x, y , 1 x 2
2
x2
2
C2
ydx 2 x dy
1 2
dx 2 x.xdx
1 2
5 x 2 dx
2 1
2
1 5 3 15
x
2 3 1 2
x x t
y y t
at b
A C
a b t
B
A
–C
We denote –C as the curve which has the same end points as C but
with the opposite orientation (from initial point B to terminal point A),
we have:
C
f x, y dx f x, y dx
C
C
f x, y dy f x, y dy
C
C
f x, y ds f x, y ds
C
Line integrals and circulation of a field 175
This is because si is always positive, whereas xi and yi change
sign when we reverse the orientation of C.
r t x t i y t j z t k
x x t , y y t , z z t , a t b
C
f x, y, z ds
2 2 2 (5.8)
dx dy dz
f x t , y t , z t
b
a dt dt dt
b
a
f r t r t dt
b
C
ds r t dt L
a
176 Mathematics for mechanical engineers
C
P x, y, z dx Q x, y, z dy R x, y, z dz (5.9)
t
x cos t , y sin t , z , 0 t 2
2
C
x 2 ds
2 2 2
dx dy dz
cos 2 t dt
2
0 dt dt dt
1 2 1
2 0
1 cos 2t sin 2 t cos 2 t dt
2
Figure 5.10. Graph of C
6 2 6
4 0
1 cos 2t dt
2
Line integrals and circulation of a field 177
r t 1 – t 2, 2, 0 t 4,6,5
2 2t , 2 4t ,5t
x 2 2t , y 2 4t , z 5t
0 t 1
Thus,
1
C1
2 y dx z dy x dz 2 2 4t 2dt 5t 4 dt 2 2t 5 dt
0
1 1
18 6t dt 18t 3t 2 0 21
0
r t 1 – t 4, 6, 5 t 4, 6, 0 4, 6,5 – 5t
or
x 4, y 6, z 5 – 5t , 0 t 1
Then, dx 0 dy . Thus,
1
C2
2 y dx z dy x dz 4 5 dt 20
0
C
2 y dx z dy x dz 21 20 1
F x, y , z P x, y , z i Q x, y , z j R x , y , z k
P0
Pi
P1
Pi-1 Pn
P2
O x
a ti-1 ti b t
Pi-1 Pi
O Pn
x P0 y
The first graph in Figure 5.12 shows the two-dimensional case. The
second shows the three-dimensional one. Suppose a point
Pi * xi* , yi* , zi* on the ith subarc corresponding to ti* . If si from Pi 1 to
Pi is small, we can approximate si as a vector with direction the
same with the unit tangent vector T ti* at Pi* . Therefore, the work
done by the force F to move the object from Pi 1 to Pi is
approximately
180 Mathematics for mechanical engineers
F xi* , yi* , zi* si T ti* F xi* , yi* , zi* T ti* si
F( x
i 1
i
*
, yi* , zi* ) T( xi* , yi* , zi* ) si (5.10)
where T( xi* , yi* , zi* ) is the unit tangent vector at the point x, y , z on
C. We can see that Obviously, when n becomes larger these
approximations will become better. Thus, the work W done by the
force field F is computed as the limit of the Riemann sums in Formula
5.9:
W F x, y, z T x, y, z ds F Tds (5.11)
C C
r t
T t
r t
Remember that:
ds x t y t z t dt r t dt
2 2 2
So, using Equation 5.8, the work done expressed in Equation 5.11 is
now rewritten:
b r t b
W F r t r t dt F r t r t dt
a
r t a
Line integrals and circulation of a field 181
b
C
F dr F r t r t dt F Tds
a C
F r t t 4i t 3 j
y
and
r t i 2tj
Figure 5.13 shows the force field and the curve in the example. The
work done is negative because the field impedes movement along the
curve.
Notice that C
F dr F dr because the unit tangent vector T in
C
t
is the circular helix given by x cos t , y sin t , z , 0 t 2 .
2
Solution: We have:
t
r t cos t i sin t j k
2
1
r ' t sin t i cos t j k
z
2
F r t sin t i cos t j 2t k
Thus,
C
F dr
2 Figure 5.14. The circular helix
F r t r t dt
0
2
2 t2
tdt 2 2
0 2 0
Figure 5.14 shows the circular helix in this example and vectors acting
at points on C.
F Pi Qj Rk
b
C
F dr F r t r t dt
a
b
P i Q j R k x t i y t j z t k dt
a
P x t , y t , z t x t
b
Q x t , y t , z t y t dt
a
R x t , y t , z t z t
Pdx Qdy Rdz
C
C
F dr P dx Q dy R dz
C
be expressed as F dr where F x, y , z – yi xj 4 zk .
C
b
a
F ( x) dx F (b) F (a) (5.12)
where F is continuous on a, b .
Equation 5.12 is also called the Net Change Theorem: The integral of
a rate of change is the net change. Similar to the Fundamental
Theorem of Calculus, the following theorem can be regarded as a
version of the Fundamental Theorem for line integrals.
184 Mathematics for mechanical engineers
f dr f r b f r a
C
f dr f r t r t dt
b
C a
b f dx f dy f dz
dt
a
x dt y dt z dt
b d
f r t dt f r b f r a
a dt
z
y B(x2, y2, z2)
C
B(x2, C
y)
A(x1, y1) O y
x
O x A(x1, y1, z1)
The last step follows from the FTC. This theorem is still true for
piecewise-smooth curves. You can prove it by subdividing C into a
finite number of smooth curves and adding the resulting integrals.
f dr f x2 , y2 f x1 , y1
C
f dr f x2 , y2 , z2 f x1 , y1 , z1
C
mMG
f x, y , z
x2 y2 z2
W F dr f dr
C C
f 2, 2,1 f 5,10,10
mMG mMG
22 22 12 52 102 102
1 1 4
mMG mMG
3 15 15
186 Mathematics for mechanical engineers
C1
F dr F dr
C2
C1
F dr F dr
C2
This says that the line integral of a conservative vector field depends
on only end points of the path, or we can say that the line integral is
independent of path.
The notation
C
vdr or C
vdr is sometimes used to indicate that the
line integral is calculated using the positive orientation of the closed
curve C. Another notation for the positively oriented boundary curve
of D is D . So the circulation can be presented by the notation
vdr . The circulation can be positive or negative, depending on the
D
orientation of C compared to the flow of F as demonstrated in Figure
5.16.
Line integrals and circulation of a field 187
C C
T
T v
a)
C
vd r 0 b)
C
vd r 0
F x, y – y , x
x2 y 2 9
Since
F dr F c t c t dt
2
C 0
2
F 3cos t ,3sin t 3sin t ,3cos t dt
0
2
0
3sin t ,3cos t 3sin t ,3cos t dt
9sin t 9 cos 2 t dt 9dt 18
2 2
2
0 0
y2
x2 1
9
F dr F c t c t dt
2
C 0
2
F cos t ,sin t sin t ,3cos t dt
0
2
0
0, 2 cos t sin t ,3cos t dt
2 2
6 cos 2 tdt 3
0 0
1 cos 2t dt 6
Line integrals and circulation of a field 189
C
F dr F dr F dr F d r
C1 C2 C1 C2
F dr 0
C
F dr 0 F dr
C1 C2
F dr F dr F dr
C1 C2
Hence,
C1
F dr F dr
C2
Theorem 5.2: C
F dr is independent of path in D if and only if:
C
F dr 0 for every closed path C in D.
We love to know if there are other vector fields that are conservative?
The following theorem says that the only vector fields that are
190 Mathematics for mechanical engineers
f ( x, y ) 0 F dr Pdx Qdy
x x C2 x C2
Line integrals and circulation of a field 191
x
f ( x, y ) P dx Q dy P t , y dt P x, y
x x C2 x x1
y
f ( x, y ) P dx Q dy Q x, t dt Q x, y
y y C2 y y1
Thus,
f f
F Pi Q j i j f
x y
f f
P and Q
x y
P 2 f 2 f Q
y yx xy x
P Q
y x
The converse of Theorem 5.4 is true only for a special type of region
as you can see in the following section.
Figure 5.21 shows some examples for different types of curves and
regions.
P Q
0, 1
y x
P Q
As , F is not
y x
conservative by Theorem 5.5.
If we choose a curve C as
shown in Figure 5.22. The
vectors in the figure that start
on the closed curve C all
appear to point in roughly the
same direction as C. Thus, it Figure 5.22. Non-conservative field
looks as if F d r 0 and so
C
F x, y y 2 2 x i 2 xy – y j
is conservative.
Solution: Let
P x, y y 2 2 x
Q x, y – y 3 xy
then:
P Q
2y
y x
Example 5.14:
Solution:
f x x, y y3 2 x (5.13)
f y x, y y 3 xy 2 (5.14)
f x, y xy3 x 2 g y (5.15)
f y x, y 3 xy 2 g y (5.16)
g y y
196 Mathematics for mechanical engineers
1 2
g y y K
2
1 2
f x, y xy 3 x 2 y K
2
r 0 2,0
r 0,3
17
F dr f dr f 0,3 f 2, 0
C C 2
F r t m r t .
W F dr
C
F r t r t dt
b
a
b
mr t r t dt
a
m bd
r t r t dt
2 a dt
m bd m 2 b
r t dt r t
2
2 a dt 2 a
m
2
r ' b r ' a
2 2
Therefore,
2 2
W 12 m v ( b ) 12 m v ( a ) (5.17)
2
The quantity 1
2 m v (t ) , which is a half the mass times the square of
the speed, is called the kinetic energy of the object. So, Equation 5.17
can be expressed as:
W K B – K A (5.18)
This says that the work done by the force field along C is equal to the
change in kinetic energy at the endpoints of C.
198 Mathematics for mechanical engineers
W F dr P dr
C C
P r b P r a P A P B
P A K A P B K B
EXERCISES
I. Evaluate the line integral, where C is the given curve:
1.
C
yds , C is described by x 2t 2 , y t , 0 t 1 .
y
2.
C x
ds , C is described by x t , y t 2 , 0 t 1 .
3.
C
x 3 yds , C is the half of the circle x 2 y 2 4, y 0 .
1. F x, y x 2 y 3i yx 3/ 2 j , r t t 2 i 2t 3 j , 0 t 1
2. F x, y xyi yx 2 j+xzk , r t t 2 i 2t 2 j t 3k , 0 t 2
4. F x, y sin xi yj cos zk , r t t 2 i tj t 2k , 0 t 1
III. Find the work done by the force field F that move an object on C.
IV. Show that the line integral is independent of path and evaluate the
integral:
ANSWERS
I.
2
y 1t 1
2.
C x
ds
0 t
4t 2 1dt
12
5 5 1
5 49 2
xy 3 ds 2 x x 1 dx
3
4.
C 0 20
6.
C
1
xe y ds e 2 y 1 e 2 y dy
0
1
3
1 e2 2 2 .
3/2
II.
352
Fdr 5t 7 4t 5 dt
2
2.
C 0 3
1 3
4.
C
Fdr 3tdt
0 2
.
Line integrals and circulation of a field 201
III.
F dr 2t 3 3t 2 t dt 2t 3 48t 2 256t dt
5 8
C 0 5
1071 1971
450
2 2
/2
C
F dr 3
0
sin10 t cos 2 t sin 2 t cos10 t dt
63
1024
IV.
2.
y
1 ye x e x , e x e x
x
1 ye dx e
x x
therefore, dy is independent of path.
C
F x, y 1 ye x i e x j f x, y x ye x is a potential
function for F.
1
1 ye dx e dy F.dr f 1,1 f 0, 0 1 .
x x
C C e
202 Mathematics for mechanical engineers
Flow of a field through a surface 203
Imagine that water is flowing through a net in the rivers. If the net
surface is perpendicular to the water flow, a lot of water will flow
through the net surface and the flux will be large. However, if the net
is parallel to the flow, water will not flow through the net surface, and
the flux will be zero. If we know the velocity field F of water, we can
calculate the total amount of water flowing through the net surface by
adding up the component of the vector F that is perpendicular to the
surface. This establishes the flow of a field through a surface.
r u, v x u, v i y u, v j z u, v k (6.1)
204 Mathematics for mechanical engineers
x x u, v , y y u, v , z z u, v (6.2)
z
u
r S
r(u, v)
D (u, v)
y
v x
Example 6.1: Identify and plot the surface with vector equation:
r u , v vi 3cos uj 3sin uk
x v, y 3 cos u , z 3sin u
y 2 z 2 9 cos 2 u 9 sin 2 u
domain as:
0 u / 2, 0 v 5
then
y 0, z 0, 0 x 5
Figure 6.3. Quarter-cylinder
So, we obtained only the quarter-cylinder
with length 5 as shown in Figure 6.3.
Grid Curve
v z
S
r
D
(𝑢 , 𝑣 )
𝐶1 y
x 𝐶
u
0 u 6
0 v 2
x 3 sin v cos u
y u cos v
z 3 sin v sin u Figure 6.3. Grid curve of a spiral tube
If v is constant, then sin v and cos v are constant. So, the grid curves
with v constant are the spiral curves. If u is constant the curves must
look like circles. The reason is that, if u is kept constant, for example
u 0 , then the equations
x 3 sin v
y cos v
z0
Flow of a field through a surface 207
r u , v xi yj zk
where
x x u, v , y y u, v , z z u, v
r u , v r0 ua vb
r x, y , z
r0 x0 , y0 , z0
a a1 , a2 , a3
b b1 , b2 , b3
x a sin cos
y a sin sin
z a cos
Figure 6.7. The sphere
as the parametric equations of the sphere.
Flow of a field through a surface 209
y 2 cos , z 2 sin , x x
where:
x 2 y 2 z 2 , y y, z z
r x , y 2 y 2 z 2 i yj z k
x f y, z , y y, z z
x x, y y
r ( x, y ) x i y j x 2 y 2 k
x r cos
y r sin
z x2 y2 r
r r , r cos i r sin j rk
r 0, 0 2
Flow of a field through a surface 211
x y z
rv (u0 , v0 ) i (u0 , v0 ) j (u0 , v0 ) k (6.3)
v v v
v z
r P0
D
(𝑢 , 𝑣 )
x C1 C2 y
u
x y z
ru (u0 , v0 ) i (u0 , v0 ) j (u0 , v0 ) k (6.4)
u u u
Example 6.8: Find the tangent plane to the surface with parametric
equations x x, y y, z x 2 y 2 at the point (1, 1, 2).
a x x0 b y y0 c z z0 0
r x i y j x2 y2 k
x y z
ru i j k i 2 xk
u u u
x y z
rv i j k j 2yk
v v v
i j k
ru rv 1 0 2 x 2 x i 2 y j k
0 1 2y
2 x 1 2 y 1 z 2 0
or
2x 2 y z 2 0
Rij D
Sij
S
z
r
Pij
Let ru* ru ui* , v*j and rv* rv ui* , v*j be the tangent vectors at Pij as
given by Equations 6.3 and 6.4. As presented in Section 3.8 of
214 Mathematics for mechanical engineers
Chapter 3, the two edges of the patch that start from Pij can be
approximated by two vectors as shown. These vectors, in turn, can be
approximated by the vectors uru* and vrv* . So, Sij is approximated
by the parallelogram determined by the vectors uru* and vrv* as
illustrated in Figure 6.13. Thus, the area of this parallelogram is:
m n
| r
i 1 j 1
*
u rv* | u v
We can see intuitively that this approximation will get closer to the
area of S as we increase the number of subrectangles. Therefore,
taking the limit of this Riemann sum will lead to the double integral:
| r r | du dv
D
u v
r u, v x u, v i y u, v j z u, v k , u, v D
then, the surface area of S is:
Flow of a field through a surface 215
A( S ) | ru rv | dA
D
where:
x y z
ru i j k
u u u
x y z
rv i j k
v v v Figure 6.14. A sphere
D , 0 , 0 2
i j k
x y z
r r
x y z
i j k
a cos cos a cos sin a sin
a sin sin a sin cos 0
Then,
a 2 sin 2 a 2 sin
2 2
A | r r | dA a 2 sin d d a 2 d sin d
0 0 0 0
D
a 2 (2 )2 4 a 2
x x , y y , z f x, y
Thus,
f f
rx i k , ry j k
x y
and
ki j
f f f
rx ry 1 0 i jk (6.5)
x x y
f
0 1
y
Then, we have:
Flow of a field through a surface 217
2 2 2 2
f f z z
| rx ry | 1 1 (6.6)
x y x y
2 2
z z
A( S ) 1 dA (6.7)
D x y
2 2
z z
A 1 dA
D x y
1 2 x 2 y dA
2 2
1 4( x 2 y 2 ) dA
D
218 Mathematics for mechanical engineers
2 2
A 1 4r 2 r dr d
0 0
2 2
d r 1 4r 2 dr
0 0
2 18 23 (1 4r 2 )3/ 2
2
(17 17 1)
6
r u, v x u, v i y u, v j z u, v k u, v D
Rij D
Sij
S
z
m n
f ( P ) S
i 1 j 1
*
ij ij
Next, we take the limit of the Riemann sum as the number of patches
increases and finally we define the surface integral of f over the
surface S as:
m n
f ( x, y, z ) dS lim
m , n
f ( P ) S
i 1 j 1
*
ij ij (6.8)
S
We can see the similarity of the definition of the surface integral with
the definition of a line integral and the definition of a double integral.
f ( x, y, z) dS f (r(u, v)) | r
S D
u rv | dA (6.9)
This equation is also true when D is not a rectangle. Formula 6.9 gives
us the way to compute a surface integral by converting it into a double
220 Mathematics for mechanical engineers
x x u, v , y y u, v , z z u, v
1 dS | r r
S D
u v | dA A( S )
y
2
Example 6.11: Compute the surface integral dS , where S is the
S
unit sphere x2 y 2 z 2 1 .
x a sin cos
y a sin sin
z a cos
0 , 0 2
| r r | a 2 sin sin
y dS (sin sin ) | r r | dA
2 2
S D
2
(sin 2 sin 2 sin d d
0 0
Flow of a field through a surface 221
2
sin 2 d sin 3 d
0 0
2 1
(1 cos 2 ) d (sin sin cos 2 ) d
0 2 0
2
1 1 1 4
sin 2 cos cos3
2 2 0 3 0 3
x x , y y , z g x, y
So, we get:
g g
rx i k , ry j k
x y
Thus,
g g
rx ry i jk (6.10)
x x
and
2 2
z z
| rx ry | 1 (6.11)
x y
f ( x, y, z ) dS
S
2 2 (6.12)
z z
f ( x, y, g ( x, y )) 1 dA
D x y
222 Mathematics for mechanical engineers
2 2
y y
f ( x, y, z ) dS
S D
f ( x, h( x, z ), z ) 1 dA
x z
z z
2 x, 2
x y
2 2
z z
S x dS D x 1 x y dA
1 1
x 1 4 x 2 4 dx dy
0 0
1 1
dy x 5 4 x 2 dx
0 0
1
1 2
. (5 4 x 2 )3/ 2
8 3 0
Figure 6.17. Surface S
1
12
27 5 5
If S is a piecewise-smooth surface, which consists of a finite sum of
smooth surfaces S1, S2,..., Sn that intersect only along their boundaries,
then the surface integral of f over S is evaluated as:
f ( x, y, z ) dS f ( x, y, z ) dS f ( x, y, z ) dS
S S1 Sn
Flow of a field through a surface 223
For S1, since this surface is cylinder so we use the cylinder parameters
r, θ and z to parametrize the surface. In this case, r=1, then, as
presented in Example 6.5 we write its parametric equations as:
x cos
y sin
zz
0 2
0 z 1 x 1 cos
Thus,
i j k
r rz sin cos 0
0 0 1
cos i sin j
Figure 6.18. Graph of S
and Jacobian is:
| r rz | cos 2 sin 2 1
2 2
z z
S z dS D (1 x) 1 x y dA
3
2 1
(1 r cos ) 1 1 0 r dr d
0 0
2 1
2 (r r
2
cos ) dr d
0 0
2
2
0
12 13 cos d
2
sin
2
2 3 0
2
Therefore,
z dS z dS z dS z dS
S S1 S2 S3
3
2
0 2 3
2 2
Oriented Surfaces
To define surface integrals of vector fields, we need to distinguish the
orientable surfaces and nonorientable surfaces. Consider the Möbius
strip as shown in Figure 6.19. If we draw a line continuously along the
Möbius trip starting from point P on one side of the trip, then it would
end up at point P, but on the other side of the trip. Then, if we
continue to draw in the same direction, it would end up at the same
point P without ever lifting having crossed an edge. This says that a
Möbius strip really has only one side, or in the other words, the
Möbius strip is nonorientable surface.
The nonorientable surfaces are outside of the scope of this text book.
We consider only orientable (two-sided) surfaces. We choose a
surface S that has a tangent plane at every point x, y, z on S (except
at any boundary point) as shown in Figure 6.20 to begin with.
Obviously, in this case the surface has two unit normal vectors n1 and
Flow of a field through a surface 225
We now want to find the formula for the normal vector of a surface. If
a surface z g x, y is given as the graph of g, we know that the
normal vector of surface z at any point can be determined by the cross
product of two unit tangent vector n rx ry . Thus, we use Equation
6.10 to build up the formula for the unit normal vector:
g g
i jk
x y
n rx ry (6.13)
2 2
g g
1
x y
ru rv
n (6.14)
| ru rv |
and
| r r | a 2 sin
r r 1
n sin cos i sin sin j cos k r ( , )
| r r | a
n r r
where ρ, v, and n are evaluated at some point on Sij. Look at this term,
we know that v n is the component of the vector ρv in the direction
of the unit vector n. Adding up these terms and taking the limit, we
derive the surface integral of the function ρv · n over S using
Definition 6.1:
v n dS ( x, y, z ) v( x, y, z ) n( x, y, z) dS
S S
(6.15)
This equation says that the rate of flow of a fluid through a surface S
is the surface integral of the function ρv · n over S. If we denote
F v , then F is also a vector field on 3 . Then, the integral in
Equation 6.15 is written in the form of abbreviation:
228 Mathematics for mechanical engineers
F n dS
S
F dS F n dS
S S
This integral is also called the flow (or flux) of F over S. We can see
that the term F n is the component of F in the normal direction of S.
Therefore, this equation interprets that the surface integral of a vector
field over S is the surface integral of its normal component over S.
ru rv r r
F dS F r r
S S u v
dS F(r(u, v)) u v ru rv dA
D
ru rv
F dS F (r r ) dA
S D
u v (6.16)
Solution: Since the surface is the unit sphere, we use the parametric
equation as:
Flow of a field through a surface 229
0 , 0 2
Therefore,
F r , r r
sin 3 cos 2 sin 3 sin 2 sin cos 2
sin 3 cos 2 sin cos 2
F dS
S
F (r r ) dA
D
0.5
-0.5
-1
1 1
0 0
-1 -1
y x
g g
F (rx ry ) ( P i Q j R k ) i jk
x y
g g
F dS P x Q y R dA
S D
(6.17)
For S1, we can see that the xy-domain D is the disk x 2 y 2 4 which
is the projection of S1 on the xy-plane as shown in Figure 6.24. Denote
that:
P x , y , z y , Q x, y , z x, R x , y , z z x 2 y 2
then,
g g
2 x, 2y
x y
F dS
S1
g g
P Q R dA
D
x y
[ y (2 x) x (2 y ) x 2 y 2 ] dA
D
( x 2 4 xy y 2 ) dA
D
2 1
(r 4r cos sin ) r dr d
2 2
0 0
2 1
r (1 2sin 2 ) dr d
3
0 0
2 1 Figure 6.24. Graph of S
0
1 2sin 2 d 0 r 3dr
1
(2 )
4 2
For the disk S2 at z 4 , its orientation is upward. So, its unit normal
vector is n k and we have:
232 Mathematics for mechanical engineers
F dS F (k ) dS z dA 16
S2 S2 D
33
F dS F dS F dS 2 16
S S1 S2
2
Q 0 E dS (6.17)
S
F K u
F dS K u dS
S S
u x, y , z C x 2 y 2 z 2
F x, y , z K u KC 2 xi 2 yj 2 zk
1
n xi yj zk
a
Thus,
2 KC 2
F n (x y2 z2 )
a
However, on S, we have:
x2 y 2 z 2 a2
Thus,
F n 2aKC
F dS F n dS
S S
2aKC dS
S
2aKCA( S )
2aKC (4 a 2 )
8KC a3
Flow of a field through a surface 235
EXERCISES
I. Evaluate the surface integral.
1. S
x 2 y 2 zdS , S is a part of the plane z 1 2x 3 y that lies above
the rectangle 0,1 0,1 .
and 0, 0,3 .
3. S
xydS , S is a part of the plane x y z 2 that lies in the first
octant.
4.
S
xy 2 dS , S is the surface z x 2 y 2 , 0 x 1 , 0 y 1 .
5.
S
y 2 z 2 dS , S is the surface z x 3 y 3 , 0 x 1 , 0 y 1 .
6.
S
yz 2 dS , S is the surface x y 2 z , 0 y 1 , 0 z 1 .
II. Find the flux of vector field F across the oriented surface S
(Evaluate the surface integral F dS for the given vector field F
S
and the oriented surface S). For closed surfaces, use the positive
orientation.
3. F x, y , z xi yj zk , S is the hemisphere x 2 y 2 z 2 9 , y 0
oriented in the direction of the positive y-axis.
ANSWERS
I.
1 3 3 x 9 11
xyzdS xy 3 3x y 11dydx
S 0 0 40
1 1 2
4. S
xy 2 dS
0 0
xy 2 2dxdy
6
1 1 3 6 2
6. S
yz 2 dS
0 0
yz 2 2 4 y 2 dydz
18
.
II.
F dS x ye 4 y 2 x 2 y 3e 4 y x 3e 2 y dxdy
1 1
2
2.
S 0 0
1
192
31 24e2 29e4
4. F dS
S S1
F dS F dS
S2
On 𝑆 :
F dS y 2 z 2 2 y 2 2 z 2 dA
S1
y 2 z 2 1
3
3 y z 2 dA
2
y 2 z 2 1
2
On 𝑆 : S2
F dS 1 dA
y z 2 1
2
Hence F dS 2
S
238 Mathematics for mechanical engineers
Green, Ostrogradski and Stoke theorems 239
D D
Q P
C
P dx Q dy
D
dA
x y
We usually use the notation for the positively oriented boundary curve
of D by ∂D. Then, the equation in Green’s Theorem can be written as:
240 Mathematics for mechanical engineers
Q P
x y dA
D
D
P dx Q dy (7.1)
Figure 7.2. A simply region can be divided into type III regions
P
C
P dx
D
y
dA (7.2)
and
Q
C
Q dy
D
x
dA (7.3)
are true, then Green’s Theorem follows immediately for the region D.
We prove Equation 7.2 by expressing D as a type I region as shown in
Figure 7.3:
D x, y | a x b, g1 x y g 2 x
P b g 2 ( x ) P
D y dA a g1 ( x) y ( x, y) dy dx
(7.4)
b
[ P ( x, g 2 ( x)) P( x, g1 ( x))] dx
a
where the last step follows from the FTC. Now we compute the line
integral in Equation 7.2. C can be rewritten as the union of four curves
C1, C2, C3, and C4.
𝑦 = 𝑔 (𝑥)
On C1 we regard x as the parameter and
get the parametric equations as: C3
x x, y g1 x , a x b D
C2
C4
So,
C1
b 𝑦 = 𝑔 (𝑥)
C1 P ( x , y ) dx a P ( x , g1 ( x )) dx
a b
Figure 7.3. Type I region
242 Mathematics for mechanical engineers
x x, y g 2 x , a x b
Then,
b
C3
P ( x, y ) dx
C3
P ( x, y ) dx P ( x, g 2 ( x )) dx
a
C2
P ( x, y ) dx 0 P ( x, y ) dx
C4
Therefore,
C
P ( x, y ) dx
P( x, y ) dx P( x, y ) dx P( x, y ) dx P ( x, y ) dx
C1 C2 C3 C4
b b
P( x, g1 ( x)) dx P( x, g 2 ( x)) dx
a a
Comparing this expression with the one in Equation 7.4, we see that:
P
C
P ( x, y ) dx
D
y
dA
A similar treatment yields Equation 7.3 for regions of type II. Then,
by adding Equations 7.2 and 7.3, we obtain Green’s Theorem. Since
Green’s Theorem gives us two choices to evaluate an integral, either a
line integral or a double integral, so we can choose the convenient
integral to work with. The following examples show that while line
integrals are complicated to evaluate, they are easier to be evaluate by
the double integrals using Green’s Theorem.
from 0, 0 to 2, 0
from 2, 0 to 1,1
from 1,1 to 0, 0
C
x 2 dx 2 xy dy
Q P
dA
D
x y
1 2 y
(2 y 0) dx dy
0 y
1
[2 xy ]xx 2y y dy
0
1
4 y y 2 dy
0 Figure 7.4. The region D
1
y 2 y3 2
4
2 3 0 3
C
(3 y x 2 3) dx 7 x y cos y dy
7 x y cos y (3 y x 2 3) dA
D
x y
2 2
(7 3) r dr d
0 0
2 2
4 d r dr 16
0 0
In Examples 7.1 and 7.2, we found that the double integral was easier
to evaluate than the line integral. Sometimes, though, it’s easier to
evaluate the line integral, and Green’s Theorem is used in the reverse
direction. For example, if P x, y Q x, y 0 on the curve C, then
the double integral over D bounded by C is
Q P
x y dA
D
C
P dx Q dy 0
P x , y 0 , Q x, y x
P x, y y , Q x, y 0
y x
P x, y , Q x, y
2 2
so that:
Q P
1
x y
Green, Ostrogradski and Stoke theorems 245
Then, from Green’s Theorem we can see that the following formulas
can be applied to find the area of D:
A
x dy y dx 1
2 x dy y dx (7.5)
C C c
x2 y2
Example 7.3: Find the area enclosed by the ellipse 2 2 1 .
a b
x a cos t , y b sin t , 0 t 2
A 1
2 x dy y dx
C
2
1
2 (a cos t )(b cos t ) dt (b sin t )(a sin t ) dt
0
ab 2
2 0
dt ab
C1
D1 C
C3
– C3
D2
C2
Q P
C1 C2
P dx Q dy
D1
dA
x y
Q P
C2 ( C3 )
P dx Q dy
D2
dA
x y
Summing up these two equations, the line integrals along C3 and –C3
cancel. Thus, we have:
Q P
C1 C2
P dx Q dy
D
dA
x y
C
e x dx 3 xy 2 dy (3xy 2 ) (e x ) dA
D
x y
3 y 2 dA
D
2
3 ( r 2 sin 2 ) r dr d
0 1
3 2
1 cos 2 d r 3dr
2 0 1
2
3 1 1 4 45
[ sin 2 ]0 r
4
2 2 1 8
Q P Q P Q P
x y dA x y dA x y dA
D D' D"
P dx Q dy P dx Q dy
D ' D
Observe that in
D '
P dx Q dy and D
P dx Q dy the line integrals
along the common boundary lines are in opposite directions, they
vanish. Then, we obtain:
Q P
x y dA
D
D '
P dx Q dy
D
P dx Q dy
P dx Q dy P dx Q dy
C1 C2
P dx Q dy
C
yi xj
Example 7.5: Suppose F x, y
, find F dr where
x2 y2 C
C
P dx Q dy
C '
P dx Q dy
Q P y 2 x2 y 2 x2
dA ( x 2 y 2 ) 2 ( x 2 y 2 ) 2 dA 0
D x y D
Therefore,
C
P dx Q dy P dx Q dy
C
That is,
C
F dr F dr
C
Thus,
C
F dr F dr
C
2
F(r (t )) r(t ) dt
0
2 ( a sin t )( a sin t ) (a cos t )(a cos t )
dt
0 a 2 cos 2 t a 2 sin 2 t
2
dt 2
0
Then,
C
F dr curl F dS
S
The Stoke’s Theorem states that “the surface integral of the curl of a
function over a surface bounded by a closed surface is equal to the
line integral of the particular vector function around that surface.”
F dr F T ds
C C
curl F dS curl F n dS
S S
C
F T ds curl F n dS
S
curl F dS
S
S
F dr (7.6)
Consider the special case where the surface S is flat and lies in the xy-
plane with upward orientation. The unit normal is n 0i 0 j k . If
F Pi Qj Rk is a vector field, then Stokes’ Theorem becomes:
Q P
C
F dr curl F dS
S S x
y
k dA
z g x, y , x, y D
curlF dS
S
(7.7)
R Q z P R z Q P
dA
D y z x z x y x y
If we parametrize C1 by equations:
x x t , y y t , a t b
x x t , y y t , z g x t , y t , a t b
b dx dy dz
C F d r a P dt Q dt R dt dt
b dx dy z dx z dy
P Q R dt
a
dt dt x dt y dt
b z dx z dy
P R Q R dt
a
x dt y dt
z z
P R dx Q R dy
C1
x y
z z
Q R P R dA
D
x y y x
Notice that step 3 is derived by using Chain Rule and the last step is
obtained by using Green’s Theorem. Since P, Q, and R are functions
of x, y, and z where z is itself a function of x and y. Thus, using Chain
Rule and rearrange, we obtain:
z z
C
F dr Q R P R dA
D
x y y x
Q Q z R z R z z 2z
R
D
x z x x y z x y xy
P P z R z R z z 2 z
R dA
y z y y x z y x yx
R Q g P R g Q P
dA
D
y z x z x y x y
Observe that the last integral coincides the double integral in Equation
6.17 in Chapter 6:
g g
F dS P x Q y R dA
S D
254 Mathematics for mechanical engineers
R Q P R Q P
We also observe that , , are three
y z z x x y
components of vector curlF . Hence, the last integral becomes:
R Q g P R g Q P
y z x z x y x y dA
D
curl FdS
S
Thus,
C
F dr curl F dS
S
F x, y, z y 2i xj z 2k
C is the curve of intersection of a half of the sphere
x 2 y 2 z 2 4, z 0 and the cylinder x 2 y 2 2 y as
shown in Figure 7.11. (Orient C to be counterclockwise
when viewed from above.)
x2 y 2 2 y
x2 y 2 z 2 4
Subtracting, we get z 4 2 y .
C
F dr could be evaluated directly as usual. However, we can see
that it can be solve easier by using Stokes’ Theorem. We have:
i j k
curl F 1 2 y k
x y z
y2 x z2
As can be seen from the figure, there are some surfaces with boundary
C, but the simplest surface is the region S in the sphere
x 2 y 2 z 2 4 that is bounded by C. Here, S is oriented upward and
the positive orientation of C is counterclockwise.
F dr curl F dS 1 2 y dA
C
S D
2 1
1 2 1 r sin r dr d
0 0
1
2 3r 2 r3
2 sin d
0
2 3 0
2 3 2
sin d 3
0
2 3
x2 y 2 z 2 4
x2 y 2 2
r t 2 cos t i 2 sin t j 2 k , 0 t 2
Thus,
r t 2 sin t i 2 cos t j
curl F dS
S
C
F dr
2
F (r (t )) r '(t ) dt
0
2
2
2 cos t sin t 2 2 sin t cos t dt
0
2
0 dt 0
0
curl F dS
S1
C
F dr curl F dS
S2
Green, Ostrogradski and Stoke theorems 257
Then,
This theorem says that the surface integral of a vector field F over a
closed surface S, which is called the flux through the surface, is equal
to the triple integral of the divergence over the region E bounded by S.
P Q R
div F
x y z
Hence,
P Q R
div F dV x dV y dV z dV
E E E E
258 Mathematics for mechanical engineers
F dS F n dS
S S
P i Q j R k n dS
S
P i n dS Q j n dS R k n dS
S S S
P
P i n dS x dV
S E
(7.8)
Q
Q j n dS y dV
S E
(7.9)
R
R k n dS z dV
S E
(7.10)
E x, y, z x, y D, u x, y z u
1 2 x, y
R u2 x , y R
dV x, y, z dz dA
E
z D
u1 x , y z
R
z dV R x, y, u x, y R x, y, u x, y dA
E D
2 1 (7.11)
Green, Ostrogradski and Stoke theorems 259
R k n dS 0 dS 0
S3 S3
Figure 13. E is type I
Then, we have: region
R k n dS R k n dS R k n dS
S S1 S2
(7.12)
R k n dS R x, y, u x, y dA
S2 D
2
R k n dS R x, y, u x, y dA
S1 D
1
R k n dS R x, y, u x, y R x, y, u x, y dA
S D
2 1
R
R k n dS z dV
S E
Similar proofs can be obtained for Equations 7.8 and 7.9 when E is
regarded as a type 2 or type 3 region, respectively. We can notice the
analogue of the methods for proving the Divergence Theorem and
Green’s Theorem.
div F yz 3 y x 2 3
x y z
The unit sphere S is the boundary of the unit ball B given by:
x2 y2 z 2 1
4
3V B 3 1 4
3
F x, y , z cos yz x 2 i 3 xyj e xy k .
div F
x
cos yz x 2 3xy e xy
y z
2 x 3x x
E x, y, z 0 x 2, 2 y 2, 0 z 2 y 2
Thus, from Ostrogradski's Theorem we have:
F dS div F dV
S E
x dV
E
2 2 2 y 2
x dz dy dx
0 2 0
2 2
xdx 2 y 2 dy
0 2
8 2 16 2
2 Figure 7.14. E is
3 3 type III region
The Ostrogradski's Theorem is also true for regions that are finite
unions of simple solid regions. The proof is similar to the one we have
done to extend Green’s Theorem.
For example, Figure 7.15 describes the solid region E bounded by the
closed surface S2. E consists of the solid region E1 that is bounded by
the closed surface S1 and the solid region E2 that is bounded by S1 and
S2, where S1 lies inside S2. Let n1 and n2 be outward normal vectors of
S1 and S2. Then, the boundary surface of E2 is: S S1 S2 . The
262 Mathematics for mechanical engineers
div F dV F dS F n
E1 S1 S1
1 dS (7.12)
div F dV F dS
E2 S
F n dS F n dS
S1 S2
(13)
F n1 dS F n 2 dS
S1 S2
F n1 dS F n 2 dS
S1 S2
F n1 dS F n1 dS F n 2 dS
S1 S1 S2
F n 2 dS
S2
Let’s consider the electric field. Suppose that we want to find the flux
Q
of an electric field E x 3 x over S2. Then, we choose S1 as a
x
small sphere with radius a and center the origin. Observe that:
Green, Ostrogradski and Stoke theorems 263
Qx Qy Qz
E x i j k
x y z
2 3/2
x y z
2 3/ 2
x y2 z2
2 2 2 2 2 3/2
Then,
div E
Qx Qy Qz
x x 2 y 2 z 2 3/ 2
y x 2 y 2 z 2 3/2
z x 2 y 2 z 2 3/ 2
x y2 z2 3x 2 x 2 y 2 z 2
2 3/ 2 1/2
Q
x y z
3
2 2 2
x z 3y x y z2
3/2 1/ 2
2
y2 2 2 2 2
Q
x y z 2 2 2 3
x z 3z x y
2 3/2
z2
1/ 2
2
y2 2 2 2
Q 0
x y z
3
2 2 2
E dS E dS div E dV E dS E n dS
S2 S1 E2 S1 S1
The result says that the integral of the field E over surface S2 is equal
to the surface integral over S1. Because S1 is a sphere, the normal
x
vector at x is . Therefore,
x
Q x Q Q Q
En x 4 x x 2 2
3
x x x x a
Q Q
E dS E n dS 2 dS A S1
S2 S1
a S1
a2
Q
2
4 a 2 4 Q
a
This result says that the electric flux of E through any closed surface
S2 that contains the origin is 4πεQ.
have:
1
div F P0 lim
a 0 V B
F dS (7.15)
a Sa
Equation 7.15 states that div F P0 is the net rate of outward flux per
unit volume at P0. If div F P 0 , P is called a source. In this case
the net flow is outward near P. If div F P 0 , P is called a sink.
The net flow is inward near P.
incoming vector are shorter than the outgoing arrows. Thus, the net
flow is inward. So, div F P2 0 and P2 is a source.
EXERCISES
I. Evaluate the line integral by using Green’s Theorem.
1.
C
ydx xdy , C is the triangle with vertices (0, 0), (3, 0), (3, 3), and
(0, 3).
2. C
xy 2 dx x 2dy , C is the rectangle with vertices (0, 0), (3, 0) and
(3, 3).
3.
C
xydx x 2dy , C is the circle with the center the origin and radius
1.
4. C
ydx xdy , C consists of segments (-3, 4) to (1, 4) and parabol
y x 2 2 x 1 from (-3, 4) to (1, 4).
II. Use the Green’s Theorem to find the work done on an object:
2. By the force F xi x3 3xy 2 j in moving the object from (-2, 0)
1/2
along the x-axis to (2, 0) then along the semicircle y 4 x 2 to
the starting point.
2. F x, y, z xe yz i ye xz j ze xy k , S is the hemisphere
x 2 y 2 z 2 4, z 0 , oriented upward.
Green, Ostrogradski and Stoke theorems 267
IV. Verify that Stokes’ Theorem is true for the given vector field F
and surface S.
V. Verify that the Frean- Ostrogradski Theorem is true for the vector
field F on the region E.
4. F x, y , z xi yj zk , E is the ball x 2 y 2 z 2 4 .
ANSWERS
I.
268 Mathematics for mechanical engineers
9
3 x
2.
C
xy 2 dx x 2dy
0 2 x 2 xy dydx 4
0
64
1 4
4. ydx xdy 2 dydx .
C 2
3 x 2 x 1 3
II.
2
W 3 r 2 .rdrd 12
0 0
III.
Thus:
S2 x2 y 2 4 x2 y2 4
IV.
Green, Ostrogradski and Stoke theorems 269
C1 : r1 t 1 t i tj, 0 t 1
F r1 0, 0, 0 , r1 1,1, 0
C2 : r2 t 1 t j tk , 0 t 1
F r2 0, t 1 t , t , r2 0, 1,1
C3 : r3 t ti 1 t k , 0 t 1
F r3 0, 0,1 t , r3 1, 0, 1
Thus:
1 1 1 1
F.dr 0 0 dt 0
t 1 t t
dt 0 t 1 dt 6
C
g g
1, 1
x y
g g
curl F.dS P x Q y R dA
S D
y xy x 1 x y dA
D
1
x 2 xy y x dydx
1 1 x
2
0 0 6
V.
2. We should show that FdS div FdV , in which the surface S
S E
FdS xy , 2 x y, z . 2 x, 2 y,1 dA
2 2
S1 D
6 x 2 y 2 1 x 2 y 2 dA
D
FdS xy , 2 x y, z . 0, 0, 1 dA x 2 y 2 1 dA
2 2
S2 D D
1 x 2 y 2
E D 0
Oxy – plane x y 1 .
2 2
4. We should show that FdS div FdV , in which the surface S
S E
is the sphere x y z 4 2 2 2
and the solid E is the ball
x2 y2 z 2 4 .
Fourier series 271
f x a0 an cos nx bn sin nx
n 1
8.2 Lemma
To study the convergence of Fourier series we need to present the
following lemma:
Lemma 8.1:
sin mx dx 0
cos mx dx 0,
m0
272 Mathematics for mechanical engineers
cos mx sin nx 0
0 mn
cos mx cos nx dx
mn
0 mn
sin mx sin nx dx mn
2
am
0
f ( x) cos mxdx, bm 0, m N
2
f ( x)sin mx dx,
bm am 0, m N
0
f x a0 an cos nx bn sin nx , x (8.2)
n 1
series f x cn x a , the coefficients is computed in terms of
n
n 1
f
n
a
derivatives cn . For Fourier series, we use integrals instead
n!
of differentials. Let's integrate both sides of Equation 8.2 and assume
that it’s permissible to integrate the series term-by-term, we get:
f ( x)dx a0 dx
a
n 1
n cos nx bn sin nx dx
2 a0 an cos nx dx bn sin nx dx
n 1 n 1
cos nx dx 0
and:
sin nx dx 0
So,
f ( x)dx 2 a
0
1
a0
2 f ( x)dx
(8.3)
f ( x) cos mx dx
a0 an cos nx bn sin nx dx
n 1
(8.4)
a0 cos mx dx an cos nx cos mx dx
n 1
bn sin nx cos mx dx
n 1
sin nx cos mx dx 0
m, n
0 nm
sin nx cos mx dx
nm
So,
f x cos mx dx am cos mx cos mx dx am
1
f x cos nx dx
an n 1, 2, 3,... (8.5)
1
bn
f x sin nx dx
n 1, 2, 3,... (8.6)
1
a0
2 f x dx
1
f x cos nx dx
an (8.7)
1
f x sin nx dx
bn
Example 8.1: Find the Fourier coefficients and Fourier series of the
square-wave function f defined by:
0 if x
f x
1 if 0 x
f x 2 f x
Using Equation 8.7 the formulas for the Fourier coefficients, are
obtained as:
1 1 0 1 1
a0 f x dx 0dx 1dx
2 2 2 0 2
1
an
f x cos nx dx
1 0 1 0
0dx
cos nx dx
1 sin nx
0
n 0
1
sin n sin 0 0
1
bn
f x sin nx dx
1 0 1
0dx cos nx dx
0
1 cos nx
n 0
1
cos n cos 0
n
0 n is even
2
n n is odd
1 2
sin 2k 1 x
2 k 1 2k 1
The graphs in Figure 8.2 present some of the partial sums of this
Fourier series. Notice that, as n increases the graph of S n x becomes
closer to the graph of square-wave function. Observe that the graph of
S n x is approaching the graph of f x , except where x 0 or x is
equal to m, where m is an integer. The graphs suggest that f x is
equal to the sum of its Fourier series except at the points where is
discontinuous.
y
y
278 Mathematics for mechanical engineers
y
y
y
f a lim x a f x
f a lim x a f x
Theorem 8.1: Let f be a function with period 2 and f and its first
derivative f are piecewise continuous on , that has only finite
jump discontinuities, then:
Fourier series 279
1
f x f x
2
f 0 lim x 0 f x 1
f 0 lim x 0 f x 0
and similarly for the other discontinuities. The average of these left
and right limits is 1/2, so for any integer n the Fourier Convergence
Theorem gives:
1 2 f x n n
sin 2k 1 x
2 k 1 2k 1 0 n n
a0 an cos nt bn sin nt (8.8)
n 1
280 Mathematics for mechanical engineers
where:
1
a0
2 g t dt
1
an
g t cos nt dt
(8.9)
1
bn
g t sin nt dt
x
Now, we substitute t into Equation 8.9 then dt
dx .
L L
Therefore, if f is a piecewise continuous function on L, L , its
Fourier series is:
n x n x
a0 an cos bn sin
n 1 L L
Here
1 L
f x dx
2 L L
a0
1 L n x
an f x cos dx
L L L
1 L n x
bn f x sin dx
L L L
f x x , 1 x 1
f x 2 f x , x
1 1
2 1
a0 x dx
1 0 1 1
x dx xdx
2 1 2 0
1 0 1 1 1
x2 x2
4 1 4 0 2
For n 1 , we compute
1 1
an x cos n x dx 2 x cos n x dx
1 0
1
x 2 1
an 2 sin n x sin n x dx
n 0 n 0
cos n x
1
2 2
0 2 2 cos n 1
n n 0 n
1
bn x sin n x dx 0
1
282 Mathematics for mechanical engineers
1 2 cos n 1
cos n x
2 n 1 n 2 2
0 n is even
2
an 2 2 cos n 1 4
n n 2 2 n is odd
1 4 4 4
2 cos x 2 cos 3 x cos 5 x
2 9 25 2
1 4
cos 2k 1 x
2 n 1 2k 12 2
1 4
f x cos 2k 1 x
2 n 1 2k 12 2
1 x 0
f x
1 0 x
Solution:
a0
1
2 f x dx 2 dx
1 0
0
dx 0
1 1
f x cos nx dx cos nx dx cos nx dx 0
0
an
0
1
f x sin nx dx
bn
1 0
sin nx dx sin nx dx
0
2
sin nx dx
0 n is even
2
1 cos nx 4
n n is odd
4
f x sin 2k 1 x
k 0 2k 1
S4 S6
S2
a0
f ( x) an cos nx bn sin nx
2 n 1
Then,
einx e inx
cos nx
2
e e inx
inx
sin nx
2i
Therefore,
Fourier series 285
and:
Now, let:
a0
c0
2
a ibn
cn n n 1
2
a b
c n n n n 1
2
then
f x ce
n
n
inx
an ibn
cn
2
1
f ( x ) cos nx dx i f ( x) sin nx dx
2
1
2
f ( x)e inx dx
cn c n
1
f ( x) cos nx dx i f ( x) sin nx dx
2
1
2
f ( x)e inx dx
1
f ( x )e
inx
cn dx
2
nx
ce
i
f x n
L
n
nx
1 L i
2 L L
cn f x e L
dx
2 y x,t 2 y x,t
(8.10)
x 2 F t 2
Fourier series 287
f(x)
rL L
Where is the mass per unit length of the string and F is the tension
in the string. The wave equation for the vibrating string is solved by
separation of variables. The solution y x, t of the wave equation is
assumed in the form:
y x,t X x T t
X x T t X x T t (8.11)
F
1
F 2
X x T t
2
X x T t
288 Mathematics for mechanical engineers
It should be noted that the left size of this equation depends only on
variable x while the right size depends only on variable t. Thus, both
sizes of this function should be a constant:
X x T t
2
X x T t
X x X x 0 (8.12)
T t T t 0 (8.13)
Let's solve Equation 8.12 first. This is the second order differential
equation. The characteristic equation of Equation 8.12 is:
k2 0
k
x
C1e C2 e x 0
X x C1 x C2 0
C1 sin x C2 cos x 0
y 0 ,t 0
y L,t 0
Fourier series 289
y x,0 f x, 0
y x,0
y x,0 0
X L C1e L
C2 e L
0
X 0 0 C1.0 C2 0 C2 0
X L 0 C1.L 0 C1 0
X 0 0 C1.0 C2 cos 0 C2 0
X L 0 C1 sin x 0
L n , n 1, 2,...
Then, we have
290 Mathematics for mechanical engineers
n 2 2
n
L2
n x
X x C1 sin (8.14)
L
k 2 2 0
k
Thus,
t
C1e C2e t 0
T t C1t C2 0
C1 sin t C2 cos
t 0
T t C1 sin
t C2 cos t
Applying the initial condition, gives:
T 0 C1 0 C1 0
Thus
n t
Tn t C2 cos (8.15)
L
n x n t n x n t
X n x Tn t C1C2 sin cos bn sin cos
L L L L
This describes the vibration of nth mode of the string. The general
solution of the wave equation y x, t is a linear superposition of the
modes,
n x n t
y x, t bn X n x Tn t bn sin cos (8.16)
n 1 n 1 L L
F
Now, if we substitute into Equation 8.16 then we have:
n x F n t
y x, t bn sin cos (8.17)
n 1 L L
n x
y x, 0 bn sin
n 1 L
2 L n x
bn f x sin dx (8.18)
L 0 L
n F
with angular frequency n . When a guitar string is plucked,
L
F
the fundamental frequency, 1 will dominates and it is
L
associated with the pitch of the sound made by the string.
X2(x)
X1(x)
X3(x)
X4(x)
x
0 x rL
f x rL
1 x / L rL x L
1 r
Using Equation 8.18, the Fourier coefficients are given by the sum of
two integrals:
2 rL x n x 2 L 1 x / L n x
bn
L 0 rL
sin
L
dx
L rL 1 r
sin
L
dx
(8.19)
2sin r n
2 2 , n0
n r 1 r
2sin r n F n t n x
y x, t cos sin (8.20)
n 1 n r 1 r L
2 2
L
bn
294 Mathematics for mechanical engineers
bn
bn
bn
As can be seen from Formula 8.20, the parameter r affects the Fourier
coefficients. Figure 8.7 shows the spectrum of Fourier coefficients for
different values of r. As observed from this figure, the harmonics are
affected significantly by r. The harmonics are weakest and the sound
is thinnest when the string is plucked in the middle, r 12 . As the
string is plucked closer to the bridge, r 110 , r 210 , the harmonics
strengthen and the sound becomes richer.
Fourier series 295
EXERCISES
Find the Fourier coefficients and Fourier series of a function f if it is
defined on the interval , and it is periodic with period 2 :
2 if x 0
1. f x
2 if 0 x
1 if x 0
2. f x
x if 0 x
3. f x 2 x
4. f x x 2 1
sin x if x 0
5. f x
0 if 0 x
1 if x 0
6. f x
1 if 0 x
1 if x 1
7. f x f x 4 f x
1 if 1 x 2
0 if 2 x 1
8. f x 1 if 1 x 1 f x 4 f x
0 if 1 x 2
1 if 2 x 1
9. f x 0 if 1 x 1 f x 4 f x
1 if 1 x 2
10. f x sin x x 1
296 Mathematics for mechanical engineers
ANSWERS
1.
1 1 0 1
a0
2 f x dx 2 2dx 2
0
2dx 0
1 1 0 1
f x cos nx dx 2 cos nx dx 2 cos nx dx
an
0
2sin nx 2sin nx
0
n n
1 1 0 1
bn f x sin nx dx 2sin nx dx 2sin nx dx
0
4 cos nx 4
n
3.
1 1
a0
2 f x dx 2 2 xdx 0
1 1
f x cos nx dx 2 x cos nx dx 0
an
1 1
f x sin nx dx 2 x sin nx dx
bn
5.
Fourier series 297
1 1 0 1 1
a0
2 f x dx 2 sin xdx 2
0
0dx
1 1 0 1
f x cos nx dx sin x cos nx dx
an 0dx
0
cos nx 1
n 2 1
1 1 0 1
f x sin nx dx sin x sin nx dx
bn 0dx
0
sin nx
n 2 1
7.
1 2 1 1 1 1 1 2
a0 f x dx 1dx 1dx 1dx 0
4 2 4 2 4 1 4 1
1 2
f x cos n x dx
2 2
an
1 1 1 1 1 2
1cos n x dx 1cos n x dx 1cos n x dx
2 2 2 1` 2 1
sin 2 x
n
1 2
f x sin n x dx
2 2
bn
1 1 1 1 1 2
1sin n x dx 1sin n x dx 1sin n x dx 0
2 2 2 1` 2 1
9.
298 Mathematics for mechanical engineers
1 2 1 1 1 1 1 2 1
a0 f x dx 1dx 0 dx 1dx
4 2 4 2 4 1 4 1 2
1 2
an f x cos n x dx
2 2
1 1 1 1 1 2
1cos n x dx 1cos n x dx 1cos n x dx
2 2 2 1` 2 1
sin 2 x
n
1 2
bn f x sin n x dx
2 2
1 1 1 1 1 2
1sin n x dx 0sin n x dx 1sin n x dx 0
2 2 2 1` 2 1
Fourier transform 299
2 nx
i
f x cn e L
n
where
2 nx
1 L /2 i
cn f x e L dx
L L/ 2
We can see that the function f is the sum of components with angular
2 n n
frequency n , or with the ordinary frequency n n . Here,
L 2 L
1
the spacing of the Fourier components is n . Substituting these
L
into the above expressions, we have:
300 Mathematics for mechanical engineers
f x F ei 2 x d (9.1)
F f x e i 2 x dx (9.2)
1
f t F eit d
2
F f t e it dt
Fourier transform 301
Where, is called angular frequency whose units are radians per second.
af a f
f g f g
Thus, we have the following theorem.
Theorem 9.1 (Linearity): Let f and F, g and G are two pairs of Fourier
integrals: f F and g G , and let a, b be constants. Then,
af bg aF bG
1 x 1 x
e xU x e e sgn x
2 2
302 Mathematics for mechanical engineers
1 x0
1
U x x0
2
0 x0
1 x0
sgn x 0 x0
1 x0
1 x 1 x
Find the Fourier transforms of components e and e sgn x .
2 2
F e xU x e x e i 2 x dx
0
1
e
1 i 2 x
1 i 2 0
1
1 i 2
1 i 2
1 2 1 2
2 2
Separating the even and odd parts of f and the real and imaginary parts of
F, gives:
1 x 1 x 1 i 2
e e sgn x
1 2 1 2
2 2
2 2
2
e x
1 2
2
i 4
sgn x
x
e
1 2
2
f x F
f * x F *
f * x F *
Proof: Here we will give the proof for the first relationship. By making
the
f x limt f x e i 2 x dx
t
t
limt f ei 2 d
t
t
limt f ei 2 d
t
t
f e i 2 d
F
F x f v
or
1 f x f x
Proof: The forward Fourier transform of a function F(x) is:
304 Mathematics for mechanical engineers
F x e i 2 x dx F x e
i 2 x
dx
F e
i 2 x
d
f x
Notice that this formula has the form of an inverse transform, with x and
exchanged, and also with a minus sign on , that is, f .
f x a e i 2 a F
ei 2 bx f x F b
Proof: We give the proof for the first equation. The forward transform
integral is:
f x a
f x a ei 2 x dx
f x a e i 2 x dx f ei 2 a d
e 2 a f ei 2 d
e 2 a F
The second equation can be also proved by using the forward transform:
ei 2 bx f x e i 2 x dx f x e i 2 b x dx F b
H H ei
f t F H ei v F f t
here f t is the output which has been amplified by the amplifier. Now
consider the special case of an amplifier with linear phase, by which we
mean the phase response is linearly proportional to frequency:
2 a
Then,
H H 0 e i 2 a
Applying the amplifying procedure for a pure tone input, cos 2 0t , the
output is:
Thus, a cosine input cos 2 0t will be passed to the output as:
H 0 cos 2 0 t a
1 1
f x cos 2 0 x F 0 F 0
2 2
1 1
f x e 0 dx f x e 0 dx
i 2 x i 2 x
2 2
1 1
F 0 F 0
2 2
1
f ax F
a a
Proof: Using the forward Fourier transform, with making the change of
variable ax , we have:
d
f ax e i 2 x dx f e i 2 / a
a
1 1
f e i 2 / a d F
a a a
f x i 2 F
Proof: Writing down the integral, and using directly the definition of the
derivative, we obtain:
f f x e i 2 x dx
f x x f x i 2 x
lim x 0 e dx
x
f x x i 2 x f x
lim x 0 e dx e i 2 x dx
x x
Shift ei 2x F F
lim x0 x x
ei 2x 1 L ' hopital
F lim x 0 F lim x 0 i 2 ei 2x
x
i 2 F
308 Mathematics for mechanical engineers
Theorem 9.9 (Integral): Suppose f is integrable with f x dx 0 and
f F , then:
x F
f x dx
i 2
F
If is integrable, then the inverse relationship also holds.
i 2
d x
f x dx f x
dx
d x
dx
f x dx f x dx i 2
x
f x dx F
d x
f x dx and f x , this relationship becomes:
dx
Removing the terms
i 2
x
f x dx F
Therefore,
x F
f x dx
i 2
f x dx F 0
f 0 F dx
Fourier transform 309
This theorem says that the value of f at 0 is equal to the area bounded by
F and the frequency-axis, and the value of F at 0 is equal to the area
bounded by f and the time-axis.
f x e i 2 x dx F
If we set 0 then,
f x dx F 0
f 0 F ei 2 0 dx
F dx
f h f h x d
t
f h x d lim t f h x d
t
t x t
limt f h x d limt f x h d
t x t
x t
limt f x h d
x t
310 Mathematics for mechanical engineers
f x h d
f x h d
b
f h f h x d
a
f g f * t g t dt f * t g t dt
f g f * t g t dt
f * t g t dt f * g
f h FH
Proof:
f h f h x d e i 2 x dx
f h x e i 2 x d dx
f h x e i 2 x dx d
Fubini's Theorem
Fourier transform 311
f h x e i 2 x dx d
f e i 2 x H d Shift Theorem
f e i 2 x d H
F H
fh F H
Correlation: f g F *G
Autocorelation: f f F *G F
2
f f f * g
f g f * g F *G
f f f * f F *F F
2
312 Mathematics for mechanical engineers
Orthogonality property
From the definition of the Fourier transform, we have:
F f t e i 2 t dt
F ei 2 t d e i 2 t dt
F ei 2 2 t d dt
0 t0
t t dt 1
0 t 0
h t t dt h 0
F d F
we have:
ei 2 dt
2 t
f x x dx f x x L dx
Definition 9.4 (Comb Function): The comb function is defined as a
infinite unit impulse train:
III x x n
n
1 x
f s x =f x III
x x
1 x
f x x x n f x x nx
n n
1 x
f x III f x nx
x x n
1 x
f s x f x III F III x
x x
1 x
f x nx f x x III x F III x
n
s
x n x
f x f nx sinc x
n
1
When the continuous analog signal is sampled at a frequency s , the
x
resulting discrete signal has more frequency components than did the
analog signal. To be precise, the frequency components of the analog
signal are repeated at the sample rate. How many samples are necessary
to ensure we are preserving the information contained in the signal? If
the signal contains high frequency components, we will need to sample at
a higher rate to avoid losing information that is in the signal. In general,
to preserve the full information in the signal, it is necessary to sample at
twice the maximum frequency of the signal. This is known as the
Nyquist rate. The Sampling Theorem states that a signal can be exactly
reproduced if it is sampled at a frequency s where s is greater than
twice the maximum frequency in the signal.
N 1
f x dx f xn xn 1 xn
n 0
F f x e i 2 x dx
X /2
f x e i 2 x dx
X /2
N
f xn e i 2 xn x
n0
X X
where the interval , is chosen such that f is acceptably close to
2 2
X X
zero for x , and x is the sampling interval.
N N
Fourier transform 317
F
f x e i 2 x dx f nx e
n
i 2 nx
x
M M
m , m ,...,0,..., 1
2 2
where
1 1
2x 2x 1
s
M M x M
N
2
F m f nx e
N
i 2 m nx
x
n
2
N
2 mn
2 i M M
N
f nx e M
x, m
2
,..., 0,..., 1
2
n
2
Notice that, the very important step in making the sum in a DFT is setting
M N so the vectors ( f nx ) and ( F m ) have the same number
N
of samples. For convenience, we can shift the vector f nx 2
N into
n
2
vector f n n 0 as follows:
N 1
N
f n N x , n
2
,..., 1
f c n
f nx , N
n 0 ,..., 1
2
Thus,
N 1 i 2 mn
F m f c n e N
x
n 0
N
1 2 mn
2
N N
F m e
i
f nx N
, m ,...,0,..., 1
m
N 2 2
2
N
into vector Fc m m0 as follows:
N 1
We can shift the vector F m 2
N
m
2
Fourier transform 319
N
F m N , m
2
,..., 1
Fc m
F m , N
m 0 ,..., 1
2
Therefore,
N 2 mn
f nx Fc m e
i
N
m 0
frequencies are contained in the wave, and we don't know how the
loudness of each frequency component is.
Recall that, any signal that is periodic in the time domain can be derived
from the sum of sine and cosine signals of different frequency and
amplitude. Such a sum is referred to as a Fourier series. In complex form
the signal with period of T can be expanded as:
2 nt
i
f t cn e T
n
T
F f t e i 2 t dt f t e i 2 t dt
0
2 nt
cne ce
T i T
T
e i 2 t dt n
i 2 n
e i 2 t dt
0 0
n n
ce
T i 2 n t
n dt
0
n
n
here: n
T
0 n
ce
T
F
i 2 n
dt T
0 cn dt cnT n
0 n
n
This says that the Fourier transform of a periodic signal gives coefficients
cn of the Fourier series of the signal at corresponding frequency 2T n .
Thus, the coefficients cn can be calculated from the Fourier transform of
the signal as:
1
cn F
T
2 mn
1 1 1 N 1 i
cn F
T N t
F m fc n e N t
N t n 0
N 1 2 mn N 1 2 mn
1 i 1 i
N
f n e
n0
c
N
N
f nt e
n 0
N
As can be seen from the spectrum presented in Figure 9.5, the harmonic
with frequency of 40 Hz has amplitude of 1, while the harmonic with
frequency of 80 Hz has amplitude of 5 as described in the given equation.
Since the spectrum is symmetric over the y-axis, we just present only one
side of the spectrum. The amplitude of one side is therefore doubled.
N=1024;
dt=1/N; t=-1:dt:1;
omega1=40*2*pi;
omega2=80*2*pi;
f=sin(omega1*t)+5*sin(omega2*t);
Fourier transform 323
NFFT=2^nextpow2(N);
F=abs(fft(f,NFFT)/N);
F=2*F;
plot(F(1:round(NFFT/2));
xlabel('Frequency','color','black','fontsize',10);
ylabel('Amplitude','color','black','fontsize',10);
324 Mathematics for mechanical engineers
EXERCISES
I. Calculate following Fourier coefficients:
1 1
1. fˆ xe ix dx
2 1
1
1 x e
1
2. fˆ ix
dx
2 1
1 /2
3. fˆ cos xe ix dx
2 /2
1 /2
4. fˆ sin xe ix dx
2 /2
1. f x t
2. f x 1
3. f x e at H t
4. f x e at H t
5) f x e
a t
6. f x te at H t
7. f x rect x 1 ei x
8. f x e x
2
Fourier transform 325
ANSWERS
I.
1.
1 1 1 1
fˆ xe i x dx x cos x i sin x dx
2 1
2 1
sin x x cos x
1
i 1 2i 1
2 1 x sin x dx 2 1 2 dx
0
i 2 sin cos
2
2.
1 1 1 1
fˆ 1 x e i x dx 1 x cos x dx
1
2 2 1
2 1
2 1 x cos x dx
0
sin x
1
2 1 2 1
0 sin x dx
2 0
1 x dx
2 0
2 1 cos
2
3.
1 /2 2 /2
fˆ cos xe i x dx cos x cos x dx
2 /2
2 0
2 /2
cos x 1 x cos x 1 x dx
2 0
/2
2 sin x 1 x sin x 1 x
2 x 1 x 1 0
326 Mathematics for mechanical engineers
cos
2 2
1
12
2 /2 2
fˆ 1 cos x cos xdx
2 0 4
II.
x x e i x dx e 0 1
1 1 0 1
1 ei x d
e
2 2 2
1
1 1 2
2
e at H t e at H t e it dt e at e it dt
0
u
1 1 1
lim u e i a t lim u e i a u
i a 0 a i i a
2 2
a
1 1 u 1
lim u e a i
a i i a a i
Laplace transform 327
f s F t e st dt
0
Proof: We have,
aF1 bF2 aF1 bF2 e st dt
0
aF e bF2 e st dt
1
st
0
a F1e st dt b F2 e st dt
0 0
a F1 b F2
Theorem 10.2 (First Shift): If there exists constants M and α such that
F t Me at , that is, F t is of exponential order, then:
e bt F t f s b
e bt F limT e bt Fe st dt
T
0
e bt Fe st dt
0
since the limit exists
s b t
Fe dt
0
f s b
1 1
1 e st dt e st
0 s 0 s
T
t limT te st dt
0
T
T t T t
te dt e st e st dt
st
0 s 0
0 s
T
T sT t st
e 2e
s s 0
T sT 1 sT 1
e 2e 2
s s s
T 1 1
limT e sT 2 e sT 2
s s s
T 1 1
limT e sT limT 2 e sT limT 2
s s s
T 1 1
limT sT limT 2 sT limT 2
se se s
1
00 2
s
Thus, we have:
1
t
s2
n!
Corollary 10.1: t n , where n is a positive integer.
s n 1
tn nt n 1 st n
t e t n 1
n
t e dt e st
n st
0 s 0
0 s s
2 2 2!
t 2 t 3 21
s s s
Assume that:
n!
t n
s n 1
t n 1
n 1 !
s n2
Put m = n + 1 then:
m ! n 1 ! n 1 !
t m n2
s m 1 s n 1 1 s
F t eat f s a
1
So, if F t t then f s .
s2
Thus,
1
F t e at
s a
2
n!
Now, if F(t) = tn then f s . Therefore:
s n 1
n!
t n e at
s a
n 1
e
is t
1
lim t
is is
1 s 2 t
e i s t e i s t e i s t e
limt limt limt 0
i s i s e i s t
i s e i s t
Thus,
1
eit
is
s 1
2 i 2
s 1 s 1
cos t i sin t
cos t i sin t
Equating real and imaginary parts, the two results are obtained:
s
cos t
s 1 2
1
sin t 2
s 1
d
tF t f s
ds
and in general:
dn
t n F t 1 f s .
n
ds n
F t e st F t dt
d d
f s e st F t dt te st F t dt tF t
ds ds 0 0
d
tF t f s
ds
dn
t n F t 1 f s
n
ds n
dn
t n 1e st F t dt 1 f s
n
0 ds n
334 Mathematics for mechanical engineers
or
d n 1
t n 1e st F t dt 1 f s
n 1
0 ds n 1
Thus,
d n 1
t n 1 F t 1 f s
n 1
ds n 1
1
Solution: If F t sin t then f s . Using Theorem 10.3 we
s 1
2
have:
d d 1
t sin f s
ds ds s 2 1
2s
s 1
2 2
1 s
F at f
a a
F at e st F at dt
0
Laplace transform 335
du s
Changing the variable u at then dt , st u and replacing
a a
these into the integral, we have:
1 as u 1 s
F at e F u du f
a 0 a a
1
sin t
s 12
1 1
sin t
s 2
s 2
2
1
F t e st F t dt F 0 sf s
F t e st F t se st F t dt
0 0
F 0 sf s
F t s 2 f s sF 0 F 0
F t e st F dt
0
e st F se st F dt
0 0
F 0 se st F s 2 e st Fdt
0 0
F 0 sF 0 s f s
2
This establish the result. The general result which can be proved by
induction, is
F
n
t s n f s s n1 F 0 s n 2 F 0 ... F n 1 0
1
cos t sin t
1
s 2 sin 0
s 2
s
2
s 2
0 t0
H t
1 t0
Thus
1
H t e st H t dt e st .1dt
0 0 s
H t t0 He st dt e st dt
0 t0
e st e st0
s t0 s
H t t0 F t t0 e st0 f s
H t t0 F t t0 e st H t t0 F t t0 dt
0
e st F t t0 dt
t0
F u du put u t t0
s u t0
e
0
e st0 e su F u du
0
e st0
f s
F t f s
1 f s F t
aF1 t bF2 s a F1 s b F2 s
af1 s bf 2 s
F1 t 1 f1 s
F2 t 1 f 2 s
and,
aF1 t a1 f1 s
bF2 t b1 f 2 s
Thus,
Finally,
a
Example 10.8: Find 1 2 2
s a
340 Mathematics for mechanical engineers
a 1 1 1
s a
2 2
2 s a s a
Notice that,
1
L 1
s
1
L e at .1 f s a
sa
1
L e at .1 f s a
sa
1 1 1 1 at at
L1 e e sinh at
2 s a s a 2
s 2
Example 10.9: Find 1 3
s 3
s2 1 6 9
s 3 s 3 s 3 s 3 3
3 2
s 2 1 1 6 1 9
1 3
1 2
3
s 3 s 3 s 3
s 3
Laplace transform 341
9
e3t 6te 3t t 2 e 3t
2
1
s 3
s s 1 s 2
s 3 3 4 1
s s 1 s 2 2 s 3 s 1 6 s 2
1
s 3 1 3 1 4 1 1
s s 1 s 2 2s 3 s 1 6 s 2
3 1 4 1 1 1
1 1 1
2 s 3 s 1 6 s 2
3 4 1
e t e 2 t
2 3 6
0 t0
t
0 t 0
h t t dt h 0
t dt 1
t 1
0
h t t dt h t t dt h t t dt h 0
0
If h t e st , we have
e st t dt e st t dt
0 0
lim 0 e st t dt
0
e st t dt
0
e0 1
Thus,
t 1
s2
Example 10.11: Determine 1 2
s 1
Solution: We have:
s2 1
1 2 1 1 2
s 1 s 1
1
1 1 1 2
s 1
t sin t
Laplace transform 343
1 f s h s f t h t
t
f h f h t d
0
Proof: We have:
f t h t e st f h t d dt
t
0 0
t
e
st
f h t d dt
0 0
t
e st f h t d dt e st f h t dtd
0 0 0
f e st h t dtd
0
e st h t dt e su h u du
0
t
e st e su h u du
0
e h s
st
Thus
f h e st h t dt
f e s h s d
h s f s
f s h s
F t sf s F 0
and
F t s 2 f s sF 0 F 0
y t 3 y t 0
where y 0 1 .
y t 3 y t 0
sy s y 0 3 y s 0
1
y s
s3
Whence,
1
y t 1 e 3t
s3
346 Mathematics for mechanical engineers
using the standard form. That this is indeed the solution is easily
checked. Let us look at the same equation, but with a different
boundary condition.
y t 3 y t 0
where y 1 1 .
Solution: Taking the same procedure as before, but we don't know the
value of y (0) , so we have the solution:
y 0
y t 1 y 0 e 3t
s3
y 1 y 0 e 3 1
which gives
y 0 e3
y 3 y cos 3t
given y 0 0 .
s
sy s y 0 3 y s
s 9
2
s
y s
s 3 s 2 9
1
1 e et
s 3
s
1 cos 3t
s 92
s
y t 1
s 3 s 2 9
t
e3 t cos 3 d
0
t
e 3t e3 cos 3 d
0
1 1 1
e 3t e3t cos 3t e3t sin 3t
6 6 6
1 1 1
cos 3t sin 3t e 3t
6 6 6
348 Mathematics for mechanical engineers
y t y t 0
with y 0 1, y 0 0 .
s 2 y s sy 0 y 0 y s 0
s
y s
s 1
2
s
y t 1 cos t
s 1
2
y t y t t
Laplace transform 349
with y 0 1, y 0 0 .
1
s 2 y s sy 0 y 0 y s
s2
s 1
y s 2 2
s 1 s s 1
2
s 1 1
y t 1 1 2 2
s 1
2
s s 1
cos t
0
t sin d
cos t sin t t
The first two terms are the complementary functions and the third is
the particular integral. The whole is easily checked to be the correct
solution. The Laplace transform method provides the solution of the
differential equation with a general right hand side in a simple and
straightforward manner.
dv
m f
dt
Here: m, v and f are the mass, velocity of the object. f is the total
forces acting on the parachutist.
dv
m mg cv
dt
or
dv c
m vg
dt m
c g
s v t v 0 v t
m s
or
c c g
s v s v0 v s
m m s
g v0
v s
c c
ss s
m m
A B C
s s c s c
m m
g 1 1 v
0
mc s s c s c
m m
g t
c c
t
v t 1 e v0 e
m m
mc
This is the velocity of the falling parachutist with the initial velocity
v 0 v0 as a function of time.
my t cy t ky t f t
352 Mathematics for mechanical engineers
m s 2 y s sy 0 y 0 c sy s y 0 ky s f s
f s ms c y 0 my 0
y s
ms 2 cs k
f s
y s
ms 2 cs k
y t 5 y t 6 y t 2e t
2
s 2 y s sy 0 y 0 5 sy s y 0 6 y s
s 1
2
s 2
5s 6 y s
s 1
s 5
2 s5
y s
s 1 s 2 s 3 s 2 s 3
Using the partial fraction method, we have:
354 Mathematics for mechanical engineers
1 1 1
y s
s 1 s 2 s 3
Finally, applying the inverse Laplace transform, the solution of the
second order differential equation is obtained:
y t e t e 2 t e 3 t
s5 3 2
y s
s 2 s 3 s 2 s 3
Thus, the solution of the homogeneous equation is found by taking the
inverse Laplace transform:
y t 3e 2t 2e3t
y t 6 y t 9 y t sin t
Laplace transform 355
1
s 2 y s sy 0 y 0 6 sy s 6 y 0 9 y s
s 1
2
1 1
y s
s 1 s 3 2
2
1
sin t
s 1
2
1
te 3t
s 3
2
1 1 t
1 2 2
e 3 sin t d
s 1 s 3
0
e 3t 3 2
5t 3 cos t sin t
50 50 25
The first term is the particular solution and is called the transient
response by engineers since it dies away for large times. The final two
terms are the complementary function and are called the steady state
response by engineers since it persists as long as the force remains.
After a “long time” the influence of the first term vanishes, the
356 Mathematics for mechanical engineers
3 2 1
cos t sin t sin t
50 25 10
with amplitude and phase form. Where, 1/10 is the amplitude of the
vibrating response and is the phase difference between input and
output. Here, cos 4 / 5 , sin 3 / 5 . It is worth noting that the
response frequency is the same as the forcing frequency that is
important in practical applications.
d4y
k 4 f x (10.1)
dx
y
f(x)
k s 4 y s s 3 y 0 s 2 y 0 sy 0 y 0 f s
Thus,
f s y 0 y 0 y 0 y 0
y s 4
2
3
(10.2)
ks s s s s4
However, the length of beam is finite, so we need to use the
Heaviside’s step function to describe the domain of the beam. Here we
358 Mathematics for mechanical engineers
g x x3
3!
1 g s f s 1 4 f s x f d
x 3
s 0
f s
1 y s 1 4
ks
y 0 1 y 0 1 y 0 1 y 0
1 2 3 4
s s s s
or
1 x
y x 1 H x L x f d
3
6k 0
(10.3)
1 1
y 0 xy 0 x 2 y 0 x 3 y 0
2 6
Laplace transform 359
y 0 0, y L 0
y 0 0 y L 0
1 x
y x 1 H x L x f d
3
6k 0
(10.4)
1
xy 0 x 3 y 0
6
The application of boundary conditions at x L is less easy. One
method is to put u x y x and Equation 10.1 becomes:
d 2u
k f x
dx 2
Applying Laplace transforms as before, the solution of this equation is
found in the form:
1 x
u x 1 H x L x f d u 0 xu 0
k 0
1 x
y x 1 H x L x f d
k 0
y 0 xy 0
y L 0
we obtain:
1 L
y 0 L f d
kL 0
(10.5)
x 6kx 0 6
Applying the boundary conditions at x = L again for this equation,
yields:
1 L 1
y 0 L f d L2 y 0
3
6kL 0 6
Substituting Equation 10.5 into this equation we obtain:
1 L L L
y 0 L f d L f d
3
(10.6)
6kL 0 6k 0
From Equations 10.4, 10.5 and 10.6, the general solution to the
problem in terms of integrals is derived:
1 x
y x 1 H x L x f d
3
6k 0
(10.7)
1 L
x L 2 L x f d
6kL 0
2 2
1 x 1 L
y x x d x L 2 2 L x 2 d
3
6 0 6 0
4 3
x x 9x
, x 0, L
24 4 8
This deflection of the beam is depicted in Figure 5.
f s F t F(t)
1 1
s
1
t
s2
1 t n 1
(n=1, 2, ...)
sn n 1!
1 1
(n=1, 2,...) t n 1e at
s a
n
n 1!
1 1
(k>0) t k 1e at
s a
k
1 1
s t
1 2
s3/2 t /
1 t x 1
, x0
sx x
1
eat
sa
s
cos at
s a2
2
a
sin at
s a2
2
s
cosh at
s a2
2
a
sinh at
s a2
2
Laplace transform 363
1 ebt e at
, ab
s a s b ba
s bebt ae at
, ab
s a s b ba
1 sin at at cos at
s 2
a 2 2
2a 3
s t sin at
s 2
a 2 2
2a
sin at at cos at
2
s
s 2
a 2 2
2a
s3 1
cos at at sin at
s 2
a 2 2
2
1 e bt e at
sa sb s b a t 3
e a / s cos 2 at
s s t
a/ s n/ 2
e
s n 1
1
a
J n 2 at
1 t
sn n t
e a/ s
H t a
s
e a / s
a
erfc
s 2 t
2
a / s
ae a / 4t
e
2 t3
1 at
erf
s sa a
364 Mathematics for mechanical engineers
1 e at erf at
s s a a
1
sa b
1
eat
t
2
beb t erfc b t
1
J 0 at
s2 a2
a sin at
tan 1
s t
t
a 0t a
1 as F t
2
tanh 2 t a t 2a
as 2 a
F t F t 2a
1 0t a
1 as F t
tanh 1 a t 2a
s 2
F t F t 2a
n/ 2
e a/ s
s n 1
t
a
J n 2 at
a as t
coth sin
a s2 2
2 a
t
sin 0ta
a F t a
a s 2 1 eas
2 2 0 a t 2a
F t F t 2a
t
1 e as F t sin 0t a
a
as 2 s 1 e as
F t F t a
1 e at / 2 3at 3at
3 sin cos e 3at / 2
s a3
3
3a 2 2 2
Laplace transform 365
s e at / 2 3at 3at
3 sin cos e 3at / 2
s a3
3
3a 2 2 2
s2 1 at 3at
e 2e cos
at / 2
s3 a3 3 2
s e at / 2 3at 3at
2
3 sin cos e3at / 2
s a3
3
3a 2 2
s2 1 at at / 2 3at
e 2e cos
s3 a3 3 2
1
1
4a 3
sin at cosh at cos at sinh at
s 4a 4
4
s sin at sinh at
s 4a 4
4
2a 2
366 Mathematics for mechanical engineers
EXERCISES
I. Find the Laplace transform:
1. t sin 2t
2. t 2 2t 1
3. t 2 cos 3t
4. t 2 e at
1. 1
s 1
s 2s 8
2
3s 7
2. 1
s 2s 5
2
e 7 s
3. 1
s 3
3
dx 2 dx
3. 2
3 2 x sin t , x 0 x 0 0
dt dt
2
dx dx
4. 2 4 2 x 5, x 0 x 0 1
dt dt
2
dx
5. 2 2 y 5sin 3t , y 0 y 0 1
dt
ANSWERS
I.
1
1. If F t sin t then f s . Using Dilation Theorem we
s 1
2
have:
1 s 1 1 a
sin at f 2 .
a a a s 2
s a2
1
a
Using Theorem 10.3 we have:
d d 2 4s
t sin 2t f s
ds ds s 2 s 2 2 2
2
II.
1. Put t x , then:
x t d dx , and:
0 xt
t x0
Thus:
t 0
f g f g t d f t x g x dx
0 t
t t
g x g t x dx g g t d g f
0 0
2. We have:
f g h 0 0 f g x d h t x dx
t x
368 Mathematics for mechanical engineers
f gh
t t
f g x h t x dxd
0
Let u x we have:
f gh
f g u h t u du d
t t
0 0
f g h
s 1
1
5
s 2s 8
2
6 s 2 6 s 4
1
s 1 1
1
1
5
s 2s 8
2
6 s 2 6 s 4
1 1 5 1 1 5
1 1 e 2t e 4t
6 s 2 6 s 4 6 6
3s 7
2. The denominator of the rational function does not
s 2s 5 2
3s 7 3 s 1 10
s 2 s 5 s 12 4
2
So
3s 7 3 s 1 10
1 1
s 2s 5 s 1 4
2 2
3 s 1 10
1 1
s 1 4 s 1 4
2 2
31
s 1 51 2
s 1 4 s 1 4
2 2
in the numerator. Thus, we can use the Second Shift Theorem since:
e 7 s
e 7 s f s
s 3
3
1
where f s . We first find F t as the inverse Laplace
s 3
3
1 1
transform of f s : 1 t 2 e3t using the First Shift Theorem.
s 3
3
2
1
This says that if F t t 2e3t . We recall the Second Shift Theorem:
2
H t t0 F t t0 e st0 f s
370 Mathematics for mechanical engineers
1 1
Now, put t0 7 , F t t 2e3t , f s to this equation, we
s 3
3
2
1 1
obtain: H t 7 t 7 e3t 7 e 7 s
2
. So
s 3
3
2
e 7 s 1
1 H t 7 t 7 e3t 7
2
s 3
3
2
REFERENCES