Chap 7
Chap 7
Chap 7
Classical Mechanics
From Thorton & Marions Classical Mechanics
Prepared by
Dr. Joseph M. Hahn
Saint Marys University
Department of Astronomy & Physics
October 17, 2004
Hamiltons Principle implies that the action J has a minimum along the
systems trajectory xi(t).
Consequently, each of the trajectories xi (t) obey the EulerLagrange eqns:
d L
L
= 0
xi dt x i
These equations are usually called the Lagrange eqns.
Note that Newtons Law can be recovered from the Lagrange eqns:
Consider the 1D motion of a particle moving in the potential U = U (x):
1
L(x, x)
= T U = mx 2 U (x)
2
L
U
so
=
=F
x
x
d
thus F =
mx = m
x as expected.
dt
Generalized coordinates
= 0 i=1,2,. . . ,s
qi dt qi
Example 7.5
Fig. 73.
What is T ?
Also need ms velocities:
a sin t + b cos
a cos t + b sin
2 + 2ab (
cos sin t + sin cos t)
x 2 + y 2 = (a)2 + (b)
sin( t)
1
1
2 + 2ab sin( t)]
so T = mv 2 = m[(a)2 + (b)
2
2
1
2 + 2ab sin( t)] mg(a sin t b cos )
and L = m[(a)2 + (b)
2
x
y
so v 2
but right ()
=
=
=
=
= 0
dt
L
= mab cos( t) mgb sin
where
L
and
= mb2 + mab sin( t)
d L
so
= mb2 + mab( ) cos( t)
dt
thus ab cos( t) gb sin b2 ab( ) cos( t) = 0
a
g
so + sin = 2 cos( t)
b
b
is the EOM.
This is the EOM for a pendulum that is driven
by an external torque (eg, the term on the right).
ie, the simple pendulum is recovered when = 0.
How would you solve the EOM?
Always keep in mind the distinction in the meaning of a partial derivative:
L
and a total derivative:
d
dt
(1)
(2)
Fig. 75
The Lagrangian is
L = T U
where U = mgz
1
T = mv 2
2
v = r
r + r + z
z = beads velocity in cylindrical coords
1
so L = m(r 2 + r22 + z 2) mgz
2
Is L written in terms of the systems generalized coordinates?
8
How do I simplify this further using the constraint imposed by the wire?
First note that = ,
and that z = cr 2 , so that z = 2crr and
L =
1
m(r 2 + r2 2 + 4c2 r2r 2) mgcr2
2
r dt
= 0
= mr 2 + 4mc2rr 2 2mgcr
= mr + 4mc2r2r
= m
r + 8mc2rr 2 + 4mc2r2r
+ 4c2rr 2 r 2 + 2gcr = 0
is the L EOM.
What is the condition for floating the bead
at some equilibrium height z = cr 2 > 0?
ie, how fast must the wire rotate for centrifugal force to balance gravity?
Since r = 0 and r = 0,
w2 = 2gc
is the angular at which the wire must spin in order to float the bead.
Example 7.9
Fig. 67
1
1
M y 2 + M R22
2
4
What is U ?
The Lagrangian is then
1
1
L = T U = M y 2 + M R22 + M gy sin
2
4
10
What does the noslip constraint tell us about the coordinates y and ?
What about the velocities?
Tip: put a dot on the disk, and use it to relate y arclength.
y = R and y = R
This allows us to write L in terms of a single generalized coordinate:
3
M R2 2 + M gR sin
4
The Lagrange equation for this system is
d L
L
= 0
dt
3
so M gR sin = M R2
2
2g
ie =
sin
3R
2gt
so (t)
=
sin assuming disk starts at rest
3R
gt2
and (t) =
sin
3R
L =
11
Problem Set #5
due Thursday November 10
at start of class
text problems 717, 727, 728, 733.
Exam #2
on Chapter 7 & Problem Sets 4 & 5
Thursday Nov. 17
The Hamiltonian H
Now lets derive another set of equations of motion from the Hamiltonian H.
This is usually obtained from the systems Lagrangian:
L
qj
L
d L
=
= p j
qj
dt qj
a simple example:
suppose L = L(x, y, x,
y)
L
L
py =
then px =
x
y
and H(px , py , x, y) = pxx + py y L(x, y, x,
y)
13
so by Chain Rule, dH =
X H
j
H
dpj +
dqj
pj
qj
H
dt
t
X
L
L
L
=
qj dpj + pj dqj
dt
dqj
dqj
q
t
j
j
j
L
L
= pj and
= p j ,
qj
qj
X
L
L
so RHS =
qj dpj
dqj
dt
q
t
j
j
Note that
Think of dH as the total change in H that results when you alter the
pj , qj , and t by small, arbitrary displacements dpj , dqj , dt.
Next bring RHSLHS:
X H
H
H L
qj dpj +
+ p j dqj +
+
dt = 0
p
q
t
t
j
j
j
Since the displacements dpj , dqj , dt are arbitrary,
what does that tell us about their coefficients?
Thus we get Hamiltons equations:
H
pj
H
pj =
qj
L
H
=
t
t
The systems H tells you how its ps and qq evolve over time.
qj =
14
15
The 7 Steps of H
L
qj
7. solve, if possible
16
Example 7.11
Fig. 79.
1. obtain L = T U :
1
1
U = kr2 = k(x2 + y 2 + z 2 )
2
2
How is U affected by the constraint?
x2 + y 2 = R2 so U = 12 k(R2 + z 2 ).
The particles KE in cylindrical coordinates is
1
T = m(r 2 + r2 2 + z 2)
2
How does the constraint alter T ?
r = R and r = 0 so T = 12 m(R2 2 + z 2 ).
17
1
1
and L = T U = m(R2 2 + z 2 ) k(R2 + z 2 )
2
2
2. get the generalized momenta:
L
= mR2 = angular momentum
L
pz =
= mz = vertical momentum
z
p =
k(R2 + z 2 )
+
2
2mR
2m 2
5. Construct H:
H(pi, qi , t) =
p2z
1
p2
+
+
k(R2 + z 2 )
=
2
2mR
2m 2
18
=
so =
p mR2
pz m
H
while p =
= 0 angular momentum is conserved
H
= kz
and pz =
z
qj =
19
H conservation
t
j
j
j
X d L L
=
qj
+
dt
t
j
j
L
= pj
qj
X d
d X
dH
L
(qj pj L) =
(qj pj L) =
=
dt
dt j
dt
t
j
but
so
N X
N
X
T (qi) =
aj,k qj qk where aj,k = ak,j
j=1 k=1
Examples include:
1
m(x 2 + y 2 + y 2) in Cartesian coords
2
1
or T = m(r 2 + r2 2 + r2 sin2 2) in spherical coords
2
T =
20
N
X
piqi L
i=1
N
X
i=1
N
X
i=1
L
qi T + U
qi
T
qi T + U
qi
Now calculate
N X
N
X
T
=
aj,k qj qk
qi
qi
j=1
=
=
consequently,
k=1
N
N X
X
j=1 k=1
N
X
ai,k qk +
k=1
N
X
i=1
N
X
aj,i qj = 2
j=1
N X
N
X
T
qi = 2
qi
N
X
k=1
ai,k qi qk = 2T
i=1 k=1
and thus H = 2T T + U = E
21
To summarize:
When (i.) L/t = 0 H = constant
But when (i.) L/t = 0
and (ii.) U = U (qi )
P
and (iii.) T = j,k aj,k qj qk H = E = constant
When conditions i., ii., and iii. hold, you can readily obtain the systems
Hamiltonian H by simply writing down the system energy E.
Just make sure that E is written as a function of the ps and qs rather than
the qs and qsyou
H
=0
qj
22
Problem 728
so
so
pr
m
p
=
mr2
r =
23
p2
k
p2r
+
H
pj
pr
H
=
pr m
p
H
=
which we already knew...
p mr2
H
qj
p2
k
H
=
r
mr3 r2
H
=0
24