An Introduction to Modern Variational Techniques in Mechanics and Engineering 1st Edition B. D. Vujanovic 2024 scribd download
An Introduction to Modern Variational Techniques in Mechanics and Engineering 1st Edition B. D. Vujanovic 2024 scribd download
https://ebookfinal.com
https://ebookfinal.com/download/an-introduction-
to-modern-variational-techniques-in-mechanics-and-
engineering-1st-edition-b-d-vujanovic/
https://ebookfinal.com/download/an-introduction-to-genetic-
engineering-3ed-edition-nicholl-d/
ebookfinal.com
https://ebookfinal.com/download/modern-trends-in-structural-and-solid-
mechanics-3-non-deterministic-mechanics-1st-edition-j-d-kaplunov/
ebookfinal.com
https://ebookfinal.com/download/materials-science-and-engineering-an-
introduction-10th-edition-edition-william-d-callister/
ebookfinal.com
https://ebookfinal.com/download/an-introduction-to-continuum-
mechanics-2ed-edition-reddy-j-n/
ebookfinal.com
An Introduction to Structural Mechanics 2001st Edition
Paul Smith
https://ebookfinal.com/download/an-introduction-to-structural-
mechanics-2001st-edition-paul-smith/
ebookfinal.com
https://ebookfinal.com/download/engineering-materials-2-third-edition-
an-introduction-to-microstructures-processing-and-design-d-r-h-jones/
ebookfinal.com
https://ebookfinal.com/download/an-introduction-to-modern-
astrophysics-2nd-edition-carroll/
ebookfinal.com
https://ebookfinal.com/download/an-introduction-to-rehabilitation-
engineering-1st-edition-rory-a-cooper/
ebookfinal.com
An Introduction to Modern Variational Techniques in
Mechanics and Engineering 1st Edition B. D. Vujanovic
Digital Instant Download
Author(s): B. D. Vujanovic, T. M. Atanackovic (auth.)
ISBN(s): 9781461264675, 1461264677
Edition: 1
File Details: PDF, 14.17 MB
Year: 2004
Language: english
B.D. Vujanovic
T.M. Atanackovic
An Introduction to
Modem Variational Techniques
in Mechanics and Engineering
QA808.V852003
53 1'.01'5 l-dc22 2003062781
AlI rights reserved. This work may noI be translated or copied in whole or in part without the written
permission of the publisher. Springer Science+Business Media, LLC.
except for brief excerpts in connection with reviews or scholarly
analysis. Use in connection with any form of information storage and retrieval, electronic adaptation,
computer software, or by similar or dissimilar methodology now known or hereafter developed is
forbidden.
The use in this publication of trade names, trademarks, service marks and similar terms, even if they
are not identified as such, is not to be taken as an expression of opinion as to whether or not they are
subject to proprietary rights.
www.birkhauser-science.com
Contents
Preface ix
Bibliography 333
Index 343
Preface
This book is devoted to the basic variational principles of mechani cs: the
Lagr ange-D 'Alemb ert different ial var iationa l principle and t he Hamil ton
integral variationa l principle. These two variational principles form t he
main subject of conte mpora ry ana lyt ical mechanics, and from th em the
whole colossal corpus of classical dynamics can be deductively derived as
a par t of physical theory.
In recent year s student s and resear chers of engineering and physics have
begun to realize the utility of variational principles and the vast possi-
bilities that they offer, and have applied th em as a powerful tool for the
st udy of linear and nonlinear problems in conservative and nonc onservative
dyn amic al systems.
The present book has evolved from a series of lectures to graduate stu-
dents and resear chers in engineering given by the aut hors at t he Depart-
ment of Mechani cs at the University of Novi Sad Serbi a , and numerous
foreign univ ersities.
The object ive of t he aut hors has been to acquaint t he reader with t he
wide possibilities to apply vari ati onal principles in num erou s problems of
cont emporary an alytical mechanics, for example, th e Noether theory for
findin g conservatio n laws of conservat ive and non conservative dynamical
syste ms, application of the Hamilton-Jacobi method and the field method
suitable for non conservative dynamical systems, the variational approach to
the mod ern op timal control theory, the application of variationa l methods
to sta bility and determining th e optima l sh ap e in th e elastic rod theory,
among others.
In order to reach a level of pr actical effectiveness , numerou s concrete ex-
amples are solved in orde r to clarify the vitality of t he t heory. It is hop ed
t ha t this book will be useful as a text in graduat e and senior und ergraduate
cours es wit h an emphas is on mechanics and/or applied ma thematics and
in graduate engineering courses. The exposition is intended to be sugges-
ti ve rather than (mathematically) rigorous . Thus, the mathematical level
has been kept as elementary as possible. Each chapter starts from widely
understood principles and brings the reader to th e forefront of the topic
in a logical way. An important part of the material present ed was already
publish ed by the a ut hors of th is book in th e numerous papers printed in
the current liter ature for th e last 10 or so years, and th e reader is dire cted
to t hese sources at the proper places in the text .
The book is divided int o two parts. The first part contains four chapters .
In t he first chapter we consider the basic form s of the Lagrange-D 'Alemb ert
x P reface
Differential Variational
Principles of Mechanics
Chapter 1
1.1 Introduction
The t ext mat erial of the pr esent chapte r is design ed to be a more or less self-
cont ained introduction to analyt ical mechanics expressed in an invar iant form
t hat is not connected to any privileged coordinate syste m . To accomplish this
goal we turn first to the Lagr ange-D'Alemb ert differential variational principle,
whose applica tions ar e very wide and encompa ss holonomic and nonh olonomi c
dyn amical syste ms and also conservative and purely nonconservative systems as
well. The elements of this part of contemporary analyt ical mechanics in fact,
constitute the content of this chapter.
f s(t , XI ,Yl ,ZI, ..., XN, YN,ZN) = 0, s = l, ...,k, wherek <3N. (1.2 .3)
The explicit dep endence on tim e in these relations means that physically the
const raints are in motion. Such constraints are usu ally referr ed to as rheonomic
or nonstationary, in cont rast to the cases when they are fixed in space or scle-
romic or stationary const ra ints, that is, they do not depend on time t explicitly,
namely 8fs/&t = O. It is to be not ed t ha t the case k = 3N is not of any int erest
IThe reader can find a rat her ex haust ive classificatio n of co nstraints as , for example, non -
holon om ic cons tmints. bilatera l, un ilateral, etc.• in t he monogr ap hs of P a rs [84], Sant illi [951.
a nd Papastavri dis [821. P a pastav rid is has also co nsi de red se rvo const raints.
1.2. Differential Equations of Motion 5
since we could solve the complete system (1.2.3) and find all 3N coordinates
Xi, Yi, Zi, (i = 1, ..., N) as functions of time t, which means that the motion of
the dynamical system is given in advance.
If the dynamical system is completely mobile without restraints, all 3N co-
ordinates Xl,Yl,Zl, .. . ,XN,YN,ZN can vary separately, and such a dynamical
system is said to have 3N degrees of freedom. Naturally, the existence of con-
straints reduces the number of independent coordinates. In fact, we can use
the equations of constraints to eliminate as many coordinates as there are con-
straints. This would bring the number of coordinates down to the number of
degrees of freedom. Namely, we can eliminate k of 3N coordinates from (1.2.3)
and express them as functions of 3N - k independent coordinates. Then, it is
said that the dynamical system has 3N - k degrees of freedom. Consequently,
the minimum number of the geometrical parameters that uniquely determine
the position of the dynamical system at each moment of time is known as the
number of degrees of freedom. Also, we tacitly assume that the constraints
(1.2.3) are independent, that is, that they have been reduced to the least pos-
sible number, which implies that the functions fI, ..., fk are not connected by
a relation B (fI, ..., fk) = O. In many practical situations the elimination of the
k redundant coordinates can be tedious or difficult, and there are advantages
in retaining more coordinates than the number of degrees of freedom. We will
pursue this possibility in the next paragraph.
Since the particles of the dynamical system are compelled to be in permanent
contact with the given constraints, we have to suppose that, as the result of
interaction between the particles and constraints, there are forces of constraints
R, = Rxiel + Ryie2 + Rzie3 , i = 1, ..., N , acting on the particles. The differential
equations of motion in the presence of holonomic constraints are
In contrast to the active forces Fi' which are fully specified, the forces of con-
straint are not furnished a priori. They are among the unknowns of the prob-
lem and must be obtained from the solution we seek. On the other hand, it
is easy to see that the problem posed by 3N differential equations (1.2.4) does
not constitute a sufficient set of equations for finding 3N unknown coordinates
Xl, YI, Zl, .. . , XN , YN, ZN as functions of time and 3N unknown orthogonal pro-
jections of the reaction forces R x i, Ryi,Rzi, i = 1, ..., N .
In order to establish a consistent problem we have to introduce some addi-
tional assumptions about the character offorces of constraints R;. It is sufficient
to require that the constraints are smooth, that is, that the reaction forces R,
are directed toward the normal of the hypersurfaces (1.2.3) and the magnitudes
of reaction forces are not limited. In addition, it is of vital importance to un-
derline some kinematical properties of the constraints.
Let us differentiate the expression (1.2.3) totally with respect to time
~(afs
L !l X,.. +!l
afs., afs .) afs_O
Y, +!l Z, + !l - , s = 1, ... . k, (1.2.5)
i=l UXi UYi UZi ut
6 Chapter 1. The Elements of Analyt ical Mechanics
T he velocity vector of the ith particle Vi = Xie l +Yie2+ ,iie 3 satisfies this relation
by all possible velocities t hat th e dyn ami cal system might have. However , t he
velocity vector that is compat ible sim ultaneously wit h (1.2.4) and (1.2.5) will be
referr ed to as th e actual velocity vector. Equivalentl y, t he actual displacement
vector of th e i th pa rticle dri = Vidt = r idt = dXie l + dYie2 + dZie3, which
satisfies at t he sam e time the different ial expression
N
~ (Of' o f. o f.
~ -0dx, + -,:;-dYi + -,:;-dzi
) of.
+ -;;-dt = 0, s = I , ..., k, (1.2.6)
i=1 Xi UYi UZi ut
(1.2.7)
This kind of displacement is int roduced in such a way t hat of smoot h constraints
and notwithst anding of scleronomic or rh eonomi c syst ems, t he relations
N
1. t
= 0, s = l , ..., k, (1.2.8)
constraints are smooth and the system is holonomic. Namely, the particles are
compelled to move on the constraints (surfaces, curves, et c.) and the reaction
forces are perpendicular to those surfaces, while the virtual displacement must
be tangent to them, hence the total virtual work vanishes . It should be stressed,
that for the rheonomic systems, the total work done by the forces of constraints
on the actual displacements dr, is not zero:
N
ER i · dr, =f O. (1.2.10)
i=l
Comparing conditions (1.2.8) and (1.2.9), we conclude that the reaction forces
can be expressed in terms of k multipliers >'5 (t) in th e following way:
k et, k et,
R Xi = 2:>'5-8. ' ti; = 2: >'5 a--:- ' i = 1, ..., N ,
5=1 X, 5=1 Y,
(1.2.11)
where the multipliers >'5 (t) are related to the magnitude of the forces of con-
straints. Therefore, the 3N differential equations (1.2.4),
(1.2.12)
become
mi X i
miZi (1.2.13)
(1.2.14)
text that follows, we will transform these differential equations into a different
form in which the Ai do not appear.
As a simple illustration of the pre ceding theory, consider the triangular prism
with inclination a th at moves with constant speed V on a horizontal plane Ox.
A particle of mass m slides down the smooth inclined face AB under the force of
gravity. The particle in motion is in contact with the moving constraint, which
represents the straight line AB moving parallel to the right, as shown in Figure
1.2.1, with constant velocity V.
a B
x
Figure 1.2.1
kx +y - k (Vt + a) = 0, (1.2.15)
.. .. \ aj \
mx = \ aj
A-
k
ax = >. ' my = -mg + A ay = A - mg, (1.2.16)
= M ' y·· -- - ~ - N, mg
i = -.!5JL
1 + k2 1 + k2 - A = 1 + k2 · (1.2.17)
x
Mt 2
= -2- + Vt, y=h-- .
Ne (1.2.19)
2
With (1.2.19) we can verify that the constraint (1.2.15) is identically satisfied.
According to (1.2.11), the projections of the reaction force Rare
k
Rx = k): = mg 1 + k 2 ' (1.2.20)
mii\ = F, (t, rI, ..., rN , VI, ..., VN) + R;, i = 1, ...,N, (1.3.1)
and
(1.3.2)
10 Chapter 1. The Elements of An alytical Mechanics
where the forces of constrai nts It; sa tisfy the condition of ort hogona lity (1.2.9):
N
LIt; · lir; =0. (1.3.3)
;= 1
that is, the velocity vector is the linear function of the generalized velocities.
Since the quantities ari/at and ari/aq5 depend only on the generalized co-
ordinates and time, it is easy to verify that the following two functional relations
are valid :
and
d ari aVi
-- = - , i = 1, ... ,N. (1.3 .13)
dt aq5 Bq,
The virtual displacement vector (or variation) of the ith particle, according to
(1.3.7), is of the form
n ari
8ri = L a8q5'
5=1 q5
i = 1, ...,N, (1.3 .14)
where the variations of the generalized coordinates are denoted by the symbol
8q5 '
As mentioned previously, the variation as a differential operator does not
produce any infinitesimal change upon time, that is, 8t = O. The variation of
velocity (1.3 .11) is found to be
(1.3 .15)
At the same time, the total time derivative of (1.3.14) (after simple interchange
of dummy indices) is of the form
d n " a2 . n a· d " a 2.
dt 8ri = LL a ~. 45 8qp + L
art dt (8q5) +L a 5~t8q5 ' (1.3 .16)
5=1 p=1 q5 qp 5=1 qs s=1 q
The difference between the last two expressions gives
that is, the variation of velocity is equal to the velocity of variation. However,
in the realm of the differential variational principles and especially in the for-
mulation of the Lagrange-D'Alembert differential principle, the commutative
rules (1.3 .18) are not obligatory. Readers can find a very exhaustive discussion
concerning the commutative and noncommutative rules in analytical mechanics
in [68], [751, and [82] .
1.3. An Invariant Form of Dynamics 13
n [N ..
~ ~m or o·- -~Fo .- 6q =0
ori N ori]
(1.3.19)
00"0 0 '0 s •
s =1 i= 1 q. i=1 q.
Transforming the first term in the following way and using relations (1.3.12)
and (1.3.13), one has
N
N •. ori
Lmiri ' - !:.- ~ m 'r o. ori
_ ~ nui» !!:.- ori
i= 1 oq. dt 0 " 0 q. 0i=1 " dto q.
d ~ . ori ~ . ori
dt 0 mir i ' q. a:- - 0 i=1
mir i '
q. a' (1.3.20)
Let us introduce the kinet ic energy of the dyn amica l system in terms of the
generalized coordinates
rN = r N(t,q,q)
wher e q = {ql , ..., qn} , q = {(1I , ..., qn} . Differentiating this function partially with
respect to generalized velocit ies q. and generalized coordinates q. we have
(1.3.22)
The second term in the brackets of equation (1.3.19) is usually term ed the
genemlized force of the corresponding genera lized coord inate q., s = 1, ..., n :
Q. =~ ori ~ (FO
0Fi . -,;- = 0
Xi 0Yi OZi)
xi-';- + Fyi-';- + F zi-';- , s = 1, ... , n .
i=1 vq . i=1 oq; vq. vq.
(1.3.24)
14 Chapter 1. The Elements of Analytical Mechanics
Since, according to (1.2.1), the active forces F, are functions of position ri,
velocities Vi, and time t, and in accordance with (1.3.7) and (1.3.11), these
vectors depend upon q's and q'e. Using the introduced notation, the virtual
work of the active forces and its transformation to the generalized coordinates
can be represented as
N n
(1.3.27)
In this case the generalized forces according to (1.3.7) and (1.3.24) can be written
which is actually the same expression for the partial derivative of the function
- I1 (t , Xl> Y1, Zl, ..., XN,YN, ZN) with respect to qs:
Qs
all
= --, (1.3.30)
oqs
(1.3.31)
1.3. An Invariant Form of Dynamics 15
t (!:...~ - or -
5=1 dt oqs oqs
Qs + OIT ) Sq,
8qs
= O. (1.3.32)
Since, according to (1.3.31), t he pote nt ial forces do not depend upon the gen-
era lized velocities qs th e last equation can be written in the form
(1.3.33)
L = L (t ,q], ..., qn,ql, ..., qn) = T (t ,q], ..., qn,ql , ..., qn) - IT (t ,ql , ..., qn) .
(1.3.34)
6L = i:
5=1
( oOL Sq, + o8~ 6qs + 8o~ (6qs)" -
qs qs qs
o8~qs (6qs)" ) • (1.3.35)
where we added and subtracted t he term (oL j oqs) (6qs)" . By using the iden-
ti ty (oL j oqs) (6qs)" = 1; [(oLj8qs) (6qs)] -1; (oLj8qs) Sq, th e last equa tion
becomes
6L = ~[(8L d 8qs
~ Oq5 - dt OL) Sq, + dtd (8qs
8L ) OL (d .)]
Sq, + 8qs dt Sq, - 6qs .
(1.3.36)
(1.3.37)
Note that this (covariant) vector will play a very important role in the subse-
quent text. Finally, with (1.3.38), the equation (1.3.37) becomes
d n n n
dt ~Ps8qs = 8L + ~ Qs8qs + ~Ps [(8qs)" - 8qs]. (1.3.39)
s=1 s=1 s=1
This is the second form of the central Lagrangian equation. According to our
best knowledge the form of this equation was first published by Lurie [68, p.
257] and Neimark and Fufaev [75, p. 133J.
It is of interest to note that the variational equation (1.3.39) can be success-
fully employed if the commutativity rule (1.3.18) is accepted. For this case we
have
(1.3.40)
(1.4.2)
where the potential function IT and generalized forces Qs have the structure
given by (1.3.21) and (1.3.24), respectively.
Note that we have used the generally accepted summation convention, which
means that whenever an index occurs two times in a term, it is implied that the
terms are to be summed over all possible values of the index. Thus, in (1.4.2)
1.4. Euler-Lagrangian Equations 17
(1.4.3)
It is easy to see that if we int egrat e th e Eul er-Lagrangian equat ions of motion
and find the generalized coord inates as functions of time ql (t ) , ..., qn (t ), then
the moti on of the dyn ami cal syste m is complete ly determined. Ind eed , ent ering
with qs (t) into (1.3 .5) we det ermine th e motion of each particle r i = r, (t) , i =
1, ..., N , and the reaction forces follow from (1.3.1) as R , = m ifi (t) --:" F i. 2 From
the pr evious analysis it follows tha t th e Euler-Lagrangian equat ions play the
cent ra l role in the study of moti on of holonomic dynamical syste ms, since they do
not contain redundant coord inates, th ey are free of equat ions of constra ints, and
they are ind ependent of unkn own Lagrangi an multipliers and reacti on forces. In
gene ra l, t he Euler-Lagrangian equations can be considered as t he corne rst one
of t he whole of ana lytical dynamics, and according t o Pars [84, p. 76] t hey
are "right ly regarded as one of t he outstanding intellectu al achievements of
manki nd ."
For the sake of completeness , note that t he Eu ler-Lagrangian equat ions
(1.4 .1) can be written in t he form
d 8L 8L
dt 8qs - 8qs = Qs, s = 1, ...,11, (1.4.4)
which follows dir ectl y from (1.3.33). Here th e Lagrangian function L is defined
as t he difference between th e kineti c energy and potential funct ion, as indicated
by (1.3.34). It is of special interest to not e that in th e case for which t he genera l-
ized nonconservative forces are equal to zero, Qs = 0, that is, the case in which
all active forces acting on t he dynamical system are of t he potent ial charac ter ,
the Euler-Lagran gian equations of motion depend solely on one function L :
.!!:.- 8L _ 8L _ 0
dt 8qs Bq, - ,
s = 1, ..., n . (1.4.5)
We ca ll this typ e of dyn ami cal system th e Lagrangian dyn amical system. The
fact t hat the Lagr an gian function depends explicitly on th e tim e as indic ated
by (1.3.34) signifies that the dynamic al syste m is subjected to th e rheonomic
constraints given by equations (1.2.3). However , if t he Lagrangian func tion is
form ed as the difference betwe en t he kineti c energy and potenti al function which
does not depend explicit ly upon ti me t, t hat is, if
(1.4.6)
2U nder li ned indices shou ld n ot be summed .
18 Ch ap ter 1. The Element s of An alytical Mechanics
which we call the potential eneryy function, and also if the kinet ic energy does
not depend on time t , t hen we call th e dynamical syste m characterized by t he
Lagrangi an
(1.4.7)
the conservative dynam ical system. Not e that the Lagr angi an dyn amic al system
whose Lagrangian function is given by L(qll " " qn, (h , ...,iln, t) Le.BL / Ot f. 0 is
also nonconservative since the tot al energy is not conserved . This fact will be
discu ssed in the proc eeding paragraphs.
T = To + T] + T2, (1.4.10)
where
where in t he last expression t he coefficients ask (t , q], ..., qn) are given by
Br; Br;
ask = aks = L m ;-Bq,
N
;= ]
.- .
8 qk
(1.4 .14)
1.4. Eul er-Lagr angi an Equ ations 19
From this expression it is seen t hat th e kineti c energy for th e scleronomic dy-
nam ical systems repr esent s a homogeneous quadrat ic form with respect to th e
genera lized velocities qs .
Tak ing into account the form of the kinetic ener gy (1.4.11)-(1.4.13) let us
write th e Euler-Lagr angian equations. First, we calculate the following expres-
sion:
(1.4.16 )
where t he symbol
Q _ an _ of _ oKs _ Daks qk
s Bq, Elqs Elt at '
k, s,m 1, ... ,n, (1.4 .19)
which are t he explicit form of the Euler-Lagrangian equat ions we have been
seeking.
As mentioned pr eviously, t he generalized accelerati ons ente r into (1.4.19)
linearl y, an d equat ions (1.4.19) can be solved with respect to lis since according
to the eq uation (1.4.3) t he matrix A with elements ams is not singular; that is,
Since the matrix A = lams] is a square nonsingular matrix, let us multiply equa-
tion (1.4.19) by inverse matrix A-I = [a:;r1] and, taking into account that
for k i= r,
-1
aksa sr = Dkr = {O1 t:
lor k -_ r, (1.4.21)
..
qr + a,r
-1 [k ] ..
rn, s qkqm = a,r
-1 (Q s -
all of oKs / oak,.)
oqs - oqs - -rit - -at qk .
(1422)
. .
..
q; + rrkmqkqm
" _ -1
- asr
(Q _0 (IlBq,+ F) _ oKat s _ oaks
s
. )
at qk· (2 )
lA. 4
It is to be noted that the Kristoffel symbols of the first and second kind , (1.4.18)
and (1.4.23), play an important role in Riemannian geometry and geometrical
interpretation of classical mechanics. For an elaborate discussion of the appli-
cations of Riemannian geometry in classical mechanics, see, for example, [106]
and [71].
In the case of scleronomic dyn amical systems, the equations of motion (1.4.19)
and (1.4.20) are reduced to the simpler forms
..
aksqk+ [krn,s ] qkqm=
. . all'
Qs - 8 s= 1, ... ,n, (1.4.25)
qs
or
.. rr"
qr + kmqkqm
- 1
= a,r (Q s -
0 (Iloq,+ F)) ' r = 1, ... ,n, (1.4.26)
(1.4.28)
1.4. Euler-Lagrangian Equations 21
Note that the function f (t, ql , ..., qn) figuring in equation (1.4.27) is frequently
referred to as the gauge function.
Moreover, in many instances it is possible to find a Lagrangian function L *
satisfying the relation (1.4.28) , where L* is not formed by the rule (1.4.27). As
an example, consider the simple harmonic oscillator whose differential equation
of motion is
(1.4.30)
L = ~ (q + qwtanwt)2 (1.4.31)
will generate the corr ect differential equation (1.4.29) . For w = 1 it is also
demonstrated in [92] and [96] that the following two Lagrangian functions,
and
q
1
L 2 = 2" (qcost - qsint) In .
(qcost - q sin
.
q sm t + q cos t
t) , (1.4.33)
will also generate the sam e differential equation (1.4.29) and all three La-
grangians (1.4.31)-(1.4.33) are not formed by the rule L = T - IT or the pre-
scription (1.4.27) .
The fact that the form of the Lagrangian functions in dynamics is not unique
raises the question of finding functions L for a given holonomic dynamical sys-
tem whose differential equat ions of motion are given in advance. This important
problem, usually referred to as th e inverse Lagrangian problem , was first con-
sidered by Helmholtz in 1887 and later studied by numerous aut hors (see, for
example, [95] , [122], and [27]).
Finally, it is of interest to note that the Euler-Lagrangian equations are form
invariant with respect to any one-to-one (i.e., punctual) transformation of two
syste ms of gener alized coordinates; the "old" generalized coordinates q}, ... , qn
are tr ansformed to the "new" generalized coordinates QI, ..., Qn by
Let us consider the given Lagrangian function L (t, ql , ..., qn,ql , ..., qn) . By using
(1.4.34) one has
(1.4.35)
22 Chapter 1. The Elements of Analytical Mechanics
For the analysis that follows we will need the following two identities, which are
similar to (1.3.12), (1.3.13):
aqi d aqi
(1.4.36)
aQj - dt aQ)'
Now we have
aL' aL aqj aL aqj
(1.4.37)
aQi = aqj aQi = aqj aQi ·
Hence ,
(1.4.38)
d 8L 8L
----=0, s=l , ...,n , (1.4.39)
dt aqs Bq;
have a very suitable form for finding conservation laws (or first int egrals) of
th e dynamical systems whose behavior can be completely described by the La-
grangian function L(t, qi , ... , qn, qI, ..., qn). By the term conservation laws or first
integrals we underst and some specific functional relations between physical and
geometrical parameters figuring in dynamical systems, which are satisfied iden-
tic ally du e to the differential equations of motion of th e dynamical system in
question. The existe nce of conservat ion laws can considerably simplify the inte-
gration of the differential equat ions of motion. Before discussing two important
1.4. Euler-Lagrangian Equations 23
8L
Ps = Bq,' s = 1, ..., n . (1.4.40)
It can happen very often that there are generalized coordinates that do not
occur in the Lagrangian function L, although their time derivatives (generalized
velocities) do . Such coordinates are usually referred to as ignorable or cyclic
coordinates.
Let qj be an ignorable coordinate, where j is a fixed particular integer. In
this case 8L/8qj = 0, and from the equation (1.4.39) for s = i , it follows that
the momentum or cyclic integral
8L
Pj = -8' = const. (1.4.41)
qj
(1.4.43)
d (8L , ) 8L (1.4.44)
dt 8qs qs - L = -fit'
If the Lagrangian function L does not depend explicitly upon time, that is,
8L/8t = 0, we arrive at the Jacobi conservation law (first integral)
8L.qs - L
-8' = E = const. (1.4.45)
qs
If the dynamical system can be represented by means of the kinetic and
potential energy, then the Lagrangian function is of the form
L =T-fI. (1.4.46)
24 Chapter 1. The Elements of Analytical Mechanics
Supposing that the kinetic energy is of the form given by equation (1.4.15),
T = (1/2) aijqiqj, and the potential energy according to (1.4.6) depends only
on qs, we have
aT.
?i7qi
..
= aijqiqj = 2T = ?i7qi,
aL. (1.4.47)
vqi v~
T + IT = E = const., (1.4.48)
which is a familiar expression for the conservation of the total mechanical energy:
the sum of the kinetic and potential energy is constant during the motion of the
system if the Lagrangian function is given in the form (1.4.46) and if it does not
depend explicitly upon time, namely aLlot = O. It can also be stated that the
existence of the energy integral is a consequence of the time invariance of the
Lagrangian function L with respect to the time translation
ll-r
m1
Figure 1.4.1
At the initial moment t = 0, r = ro, and r (0) = O. The kinetic and potential
energy of the system are
It is obvious that the generalized coordinate 'P is ignorable and the momentum
(cyclic) integral reads
P<p = et. 2· C
o<p = mlr 'P = = const. (1.4.52)
Since the Lagrangian function does not depend explicitly on time t, we have
also the energy integral T + IT = E, which in our case becomes
(1.4.53)
Since our dynamical system has two degrees of freedom, we can base our analysis
on two conservation laws (1.4.52) and (1.4.53), ignoring the Euler-Lagrangian
equations of motion.
At t = 0 the constant C is found to be C = mlrOvO, and from (1.4.52) it
follows that
. rovo
'P = -;:2 ' (1.4.54)
(1.4.55)
Entering with (1.4.54) and (1.4.55) into (1.4.53) we find after simple calculation
the following relation depending only upon variable r :
(1.4.56)
The roots of the quadratic expression with respect to r in the square brackets
are found to be
(1.4.57)
Therefore, rl > 0, r2 < 0, and the particle ml will oscillate between the bound-
ary circles "o and rl. The equation (1.4.56) can be written in the form
2m2g 1
f2 = 2 (ro - r) (r - rl) (r - r2) = <I> (r). (1.4.58)
ml +m2r
26 Ch apter 1. The Elements of Analytical Mechanics
From th is equation it follows that <I> (r) must be positive during moti on . Since
the root r 2 is negative, th e expression
(1.4.59)
is always positiv e. The right-hand side of (1.4.58) wiIl be positive in the following
two cases:
The trajectories of th e particles ml for the cases (a) and (b) are depicted in
Figure 1.4.2a and b.
.. . G
......•.........•.•.....
<~:~- - -
A
o
\.. \\ r
" .....
.................................
Figure 1.4.2
It is clear that the case (a) is the motion of the particle ml for which , at the
initial moment t = 0, the cent rifugal force of ml is greater than the weight of
the particle m 2, t ha t is,
(1.4.61)
(1.4 .62)
(1.4.63)
and t he distan ce of the particle m l will st art to decrease from the initi al dist ance
corr espond ing to point A . Finally it may happen t hat the cent rifuga l force of
lA . Euler-Lagrangian Equations 27
the mass m1 is equal to the weight m29 of the particle m2, and we have the
motion of m1 along the circle of the radius ro for which '
(1.4.64)
Let us briefly describe the motion for the case (a). During the initial increasing
period from A to B the angular velocity <p given by (1.4.54) will decrease;
the velocity and centrifugal force are going to decrease also until the point m1
reaches the point B, where the weight m19 is equal to the centrifugal force. In
the subsequent motion from B to C the regime of motion is going to be the
opposite.
From equation (1.4 .64) it follows that
dr
dt=±--· (1.4.65)
J<I> (r)'
since <I> (r) > 0 we have to take the plus sign when r is increasing and the minus
sign when r is decreasing. For the case (a) the radius is increasing from ro to
r1 in accordance with the equation
dr
L
T
t- (1.4.66)
- TO J<I> (r)'
Therefore, the period of this oscillatory motion from the point A to B and back
to C is finally
(104.69)
a = r1 - rO. (1.4.70)
3The circular motion of the p article ml and stability of the motion for which (1.4.64) is
satisfied is considered in [681, pp . 635-638.
28 Chapter 1. The Elements of Analytical Mechanics
The integral figuring in (1.4.69) is the elliptic type, and the integration can be
accomplished for given ro, r}, and r2.
It is of interest to note that for this type of oscillatory motion we can use
two equivalent formulas for the approximate period proposed by Pars ([84, p.
10]):
or
where
with
A = 2m2g (1.4.75)
m1+ m2
Finally, the second generalized coordinate <p can be found from the conservation
law (1.4.54)
(1.4.76)
(1.4.77)
lA. Euler-Lagrangian Equations 29
which is a basic invariant of Vn . Note that the quantity aij (q1, ..., qj) represents
the covariant metric tensor in the configuration space Vn . Note also that the
quantity ds introduced by (1.4.77) is referred to as the basic line element in
Vn and configuration space defined in this way is also the Riemannian space.
Together with the covariant metric tensor aij we can also introduce the con-
travariant metric tensor aij (which is denoted by aij1 in equation (1.4.26)), and
we have that aika kj = 5{ , where 5{ denotes the Kronecker delta symbol and
is equal to unity or zero when i and j are equal or different, respectively. Let
us consider an arbitrary scleronomic dynamical system subject to generalized
forces Qi (ql, ..., qn) . Using the tensorial notation we can write the differential
equations of motion of this system in two different forms (already indicated
by (1.4.25) and (1.4.26)) . In the covariant form the differential equations of
motion are (nonconservative forces are equal to zero , Qs = 0, in (1.4.25) and
Qi = -8IT/8qi )
aiih + [jm,i]tfti'" = Qi; S = 1, ... ,n. (1.4.78)
As demonstrated earlier in this section (see (1.4.19)), these equations can be
derived from the corresponding Euler-Lagrangian equations.
The equations of motion can also be written in the contravariant form
;l+r;kljjljk=Qi=aijQj, r=l, ..., n . (1.4.79)
Note that the symbols rim, j] and r;k
denote the Kristoffel symbols of the first
and second kind, respectively, introduced earlier in this section.
After this introduction in the tensor notations, we now consider two mo-
tions: the undisturbed motion qi (t) and the infinitesimal disturbed motion
xi (t). These two motions take place along the neighboring curves in the con-
figuration space Vn . We suppose that the undisturbed motion qi (t) satisfies
identically the differential equations (1.4.79). We call these trajectories, accord-
ing to Synge, the natural trajectories C. Let qr, r = 1, ..., n, be a point P of
C and qr + z" the coordinates of the corresponding (simultaneous) point P* of
disturbed natural trajectory C., z" being infinitesimally small. We shall call the
vector with contravariant components x T the disturbance vector. If the distur-
bance vector between simultaneous configurations remains permanently small,
we say that the undisturbed motion is stable in the kinematical sense .
In order to obtain the differential equations of the disturbed motion, we
substitute the generalized coordinates qi by qi + xi in (1.4.79) and calculate the
corresponding Kristoffel symbols and generalized forces in the disturbed region
as
(104.80)
where the asterisk denotes the quantities that should be calculated at p •. There-
fore, expanding r.A and Q. i and retaining the first powers of small quantities
only, we have
(a)
o xi . i r i . j k
6t =x + jkq X (1.4.83)
(a) ,
d 0 x'
dt Ot
X + riJ'kqooj Xk + riJ'kqii :Xk + -
ooi 8 ri'k'l;j q'm Xk
8q-
m J
Substituting ijj from (1.4.79) in this expression we obtain the second absolute
derivative of the disturbance vector in the form [68, p. 624]
Substituting Xi from this equation into the equation of disturbed motion (1.4 .82) ,
we have after simple manipulations the following equation of disturbed motion
in the contravariant form [106, p. 79J, [68, p. 625J:
(1.4.86)
4Concerning the geometrical interpretation of (1.4 .82) which follows, Synge wrote [106 ,
p. 78]. "The use of the tensorial notation is of greatest importance. The appearance of the
Riemannian curvature tensor makes it difficult to believe that similar results could be obtained
without the use of this method ."
lA. Euler-Lagrangian Equations 31
8 ri
R kjm j -- 8qj 8 ri rs ri rims rsjk· (1.4.87)
mk - 8qm jk + mk js -
where
i 8Qi i k
\JjQ = -8. + rjkQ
qJ
(1.4.89)
where Rkmjr denotes the covariant form of the Riemann-Kristoffel tensor Rprmn
= a ps R rmn", The explicit form of the Riemann-Kristoffel tensor is
1 8 2a 8 2a 8 2a 8 2a
R rsmn - --
28 8
rn
-+ - sm
8qr aqn
- rm
- -aqSaqn
--- -sn-
aqr aqm
qs qm
Note that for the case of the potential forces, that is, Qi = -aITj8qi, where
IT = IT (qI ,..., c") is the potential energy, the equations (1.4.90) become
(a)
air
2 i
sot + R kmjrq·kq·j
2
X
X
m
= -X m ( a
2
IT
8q m8qr -
r srm 8qS
aIT)
= -X m'<7v m (8IT)
8qr .
(1.4.92)
It is fair to say that the Dutch had several frigates and brigs
abreast the intervals in their line, which did good service, raking the
English ships as they came through and luffed up to leeward of the
enemy’s line.
As it was, Admiral Duncan met and fought the Dutch fleet before a
98- and two 74-gun ships which De Winter had expected could join
him.
Admiral De Winter, in his official report of the action, attributed his
failure to four causes: first, the superiority of the British in large
ships; secondly, their having been together at sea for many weeks,
and hence well accustomed to work together; thirdly, the advantage
of the attack, and fourthly, the early retreat of some of his ships,
and the bad sailing of some of the others. He also expressed his
belief that, if his signals had been obeyed as promptly as Admiral
Duncan’s were, some of the English fleet would have been brought
into the Texel, instead of the Dutch ships going to England. His
statement about the English ships being so long together was not
altogether correct. Captain Williamson, of the English ship Agincourt,
64, was court-martialed for his conduct in this action. He was
accused of disobedience of signals and failure to go into action; and
also, on a second charge, of cowardice or disaffection. The first
charge was found proved, but not the second, and Williamson
received a very severe sentence. It was proved on this trial that
some of Admiral Duncan’s fleet did not know other ships in the same
fleet. In the great fleet actions of those days, between ships-of-the-
line, it was not customary for frigates and smaller ships to fire, or to
be fired at, unless they provoked it; and the Dutch frigates,
corvettes and brigs formed in this action a second line, and fought
well. The Dutch were, indeed, an enemy not to be despised, and
Admiral Duncan did full justice to the determined way in which most
of them fought.
Scarcely was the British fleet, with its prizes, pointed to the
westward, when a gale of wind came on, which scattered and
endangered the whole of them. The injured masts fell, and the
vessels leaked through shot-holes which in any ordinary weather
would have been above the reach of the water.
On the 13th, the Delft, 50, a prize, exhibited a board with the
words chalked on it “The ship is sinking.” Assistance was sent, and
most of the men removed; but several of the prize crew and many of
the prisoners perished in her, so quickly did she founder.
The Monnikendam frigate was wrecked on a shoal; but all on
board were saved; and the Ambuscade frigate, being driven on the
Dutch coast, was recaptured. One by one the rest of the scattered
fleet and the prizes reached English ports.
Admiral Duncan was made a peer, and Vice-Admiral Onslow a
baronet, for this action. Gold medals were presented to the flag-
officers and captains; and the thanks of Parliament were voted to
the fleet.
We often realize more of the real spirit of a fight from private
accounts and comments than from the official reports; and we,
therefore, add a few remarks and anecdotes from such sources. In
the first place, the promptitude and decision of Admiral Duncan on
meeting the Dutch fleet is especially to be noticed. “The British
Admiral soon perceived that if he waited to form his line (the enemy
drawing fast in with the land) there would be no action.” He,
therefore, hoisted the signal to make all sail, break the line, and
engage the enemy to leeward; and for close action, which last signal
flew until it was shot away. This signal could not be mistaken, and,
coupled with the gallant Admiral’s example, superseded all former
ones.
If further proof of the superior efficacy of such a mode of attack
be wanting, it is to be found not only in the declaration of the brave
Dutch Admiral, but also in the testimony of Lord Nelson, who,
although not acquainted with Lord Duncan, wrote to him, after the
battle of the Nile, to tell him how “he had profited by his example.”
The Dutch Admiral De Winter said, “Your not waiting to form a line
ruined me; if I had got nearer to the shore, and you had attacked, I
should probably have drawn both fleets on, and it would have been
a victory to me, being on my own coast.”
It is a fact that many of the vessels of Admiral Duncan’s fleet were
intended for Indiamen, and not so stoutly built as men-of-war
usually are; and many of his ships were in bad condition, and had
not had time to complete their stores when called away from
Yarmouth Roads to encounter the enemy.
Among other incidents of this action, it is recorded that, when the
main-top-gallant mast of the Venerable was shot away, a seaman
named Crawford went aloft with another flag, and hammer and
nails, and nailed the flag to the topmast-head.
Had Duncan’s fleet been of as good material as that of Lord St.
Vincent, it is probable that every Dutch ship would have been taken.
When the action ceased the English fleet were in only nine fathoms
of water, and a severe gale was nearly upon them; and the wonder
is that they saved themselves and so many of their prizes, in their
battered condition.
Captain Inglis, of the Belliqueux, of 64 guns, owing either to a
long absence from active service, or an inaptitude to the subject,
sometimes apparent in sea officers, had neglected to make himself a
competent master of the signal-book, and on the morning of the day
of the battle, when it became necessary to act with promptitude in
obedience to signals, found himself more puzzled than enlightened
by it, and, throwing it with contempt upon the deck, exclaimed, in
broad Scotch: “D—n me, up wi’ the hellum, and gang intil the middle
o’t!”
In this manner he bravely anticipated the remedy in such cases
provided by Nelson, who, in his celebrated “Memorandum,” observes
that, “when a captain should be at a loss he cannot do very wrong if
he lay his ship alongside of the enemy.”
In strict conformity with this doctrine the Belliqueux got herself
very roughly treated by the van of the Dutch fleet.
Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.
ebookfinal.com