1 Graded Problems: PHY 5246: Theoretical Dynamics, Fall 2015 Assignment # 2, Solutions
1 Graded Problems: PHY 5246: Theoretical Dynamics, Fall 2015 Assignment # 2, Solutions
1 Graded Problems: PHY 5246: Theoretical Dynamics, Fall 2015 Assignment # 2, Solutions
Assignment # 2, Solutions
1 Graded Problems
Problem 1
(1.a)
The coordinates can be written using spherical coordi-
nates as:
(1.b)
The equilibrium points are defined as the points where the bead, if placed there (i.e. with zero
velocity), does not move. They are given therefore by the condition
2 g
θ̈ = 0 −→ sin θ ω cos θ − =0 , (5)
R
and are g g
θ0 = 0, π , and θ0 = arccos if <1 . (6)
Rω 2 Rω 2
Let us notice that we could get to the same result by interpreting the r.h.s of Eq. (4) as a force,
F (θ) and defining an effective potential, Veff (θ) such that
∂Veff (θ)
F (θ) = −
∂θ
1
Veff (θ) = − mR2 sin2 θω 2 − mgR cos θ .
2
Then, the equilibrium condition is equivalent to the condition that gives the extrema of Veff (θ),
i.e.
∂Veff (θ)
θ̈ = 0 −→ F (θ) = 0 −→ =0 . (7)
∂θ
The nature of the equilibrium points (stable vs unstable) can then be determined by looking at
the second derivative of Veff (θ),
∂ 2 Veff (θ)
2
= −mR2 ω 2 cos2 θ + mR2 ω 2 sin2 θ + mgR cos θ , (8)
∂θ
at each equilibrium point. We find that θ0 = π is always an unstable equilibrium point, since:
∂ 2 Veff (θ)
= −mR2 ω 2 − mgR < 0 , (9)
∂θ2 θ0 =π
while the nature of θ0 = 0 and θ0 = arccos Rωg 2 depends on the ratio g/(Rω 2). Indeed,
∂ 2 Veff (θ)
g
2 2
= mR ω − 1 (10)
∂θ2 θ0 =0 Rω 2
and θ0 = 0 is a stable equilibrium point if g/(Rω 2) > 1 while it is unstable if g/(Rω 2) < 1.
Viceversa,
∂ 2 Veff (θ)
g 2
2 2
= mR ω 1 − , (11)
∂θ2 θ0 =arccos( g ) Rω 2
Rω 2
and θ0 = arccos(g/(Rω 2)) is a stable equilibrium point if g/(Rω 2) < 1, while it is unstable if
g/(Rω 2) > 1. The shape of Veff (θ) therefore changes in going from g/(Rω 2) > 1 to g/(Rω 2) < 1.
In the first case Veff (θ) has a minimum at θ0 = 0 and no other extrema except a maximum at
θ0 = π, while in the second case, θ0 = 0 becomes a maximum (as well as θ0 = π) and a new
minimum develops at θ0 = arccos(g/(Rω 2)). Of course the potential is symmetric with respect
to θ, so an analogous discussion holds for θ between π and 2π. We can actually say that for
g/(Rω 2) < 1 there are two stable equilibrium points, one to the right and one to the left of the
θ = 0 position, both at an angle θ0 = arccos(g/(Rω 2)) from the vertical.
(1.c)
To find the frequency of small oscillations about the stable equilibrium positions (θ0 ) we consider
a small displacement from θ0 (let’s call it δ), i.e. we write
θ = θ0 + δ , (12)
and plug it into the equation of motion. We then expand the equation of motion in δ and its
derivatives, keeping only up to linear terms, i.e. we write
sin(θ0 + δ) = sin θ0 cos δ + cos θ0 sin δ = sin θ0 (1 + O(δ 2 )) + cos θ0 (δ + O(δ 3 ) , (13)
cos(θ0 + δ) = cos θ0 cos δ − sin θ0 sin δ = cos θ0 (1 + O(δ 2)) − sin θ0 (δ + O(δ 3 ) ,
and we get the equation of motion for small oscillations about θ0 in the form
g h g i
δ̈ = ω 2 sin θ0 cos θ0 − sin θ0 + δ ω 2 cos2 θ0 − sin2 θ0 − cos θ0 . (14)
R R
The terms in the first parenthesis cancel because θ0 is a solution of θ̈ = 0 (or equivalently a
minimum of Veff and therefore satisfies the extreme condition in Eq. (7)), and the equation for δ
reduces to
δ̈ + Ω2 δ = 0 , (15)
which is the equation of a one-dimensional harmonic oscillator with frequency
g
Ω2 = cos θ0 − ω 2 cos2 θ0 − sin2 θ0 .
(16)
R
We then have two cases:
• g/(Rω 2) > 1, the stable equilibrium point is θ0 = 0 and the frequency of small oscillation
about θ0 is
g 1/2
Ω= − ω2 , (17)
R
• g/(Rω 2) < 1, the stable equilibrium point is θ0 = arccos(g/(Rω 2)) and the frequency of
small oscillation about θ0 is
1/2
g2
Ω = ω sin θ0 = ω 1 − 2 4 . (18)
Rω
At the critical value of ω 2 = g/R, the system is transitioning between the two cases, and the
critical values all occur at θ = 0.
Problem 2
The system has two degrees of freedom and given the symmetry of the prob-
lem, it is natural to choose polar coordinates, r and θ. Indeed, we will denote
by r the elongation of the spring from its unscratched length. They are the z
generalized coordinates of the problem, in terms of which we can write the
kinetic and potential energy of m as x
1
T = m(ṙ 2 + (l + r)2 θ̇2 ) , (19) l
2
1 θ
V = −mg(l + r) cos θ + kr 2 ,
2 r
m
and the Lagrangian of the system as
1 1
L = T − V = m(ṙ 2 + (l + r)2 θ̇2 ) + mg(l + r) cos θ − kr 2 . (20)
2 2
The equation of motion for r is then
m(r̈ − (r + l)θ̇2 ) = mg cos θ − kr , (21)
while the equation of motion for θ is
m[(l + r)θ̈ + 2ṙ θ̇] = −mg sin θ . (22)
Notice that Eq. (21) and (22) are of the form mar = Fr and maθ = Fθ , where (ar , Fr ) and (aθ , Fθ )
are the components of the acceleration and of the total applied force in the radial and tangent
direction respectively.
y CM
θ
a
α y
x
x
φ
l−r
the origin is conserved,
∂L pθ m2
pθ = = m1 r 2 θ̇ = constant ⇒ θ̇ = .
∂ θ̇ m1 r 2
Thus we can rewrite the r equation of motion as
p2θ
(m1 + m2 )r̈ − + m2 g = 0,
m1 r 3
which is now just a 1d problem.
2 Non-graded Problems
Problem 5 (1.15 of Goldstein’s book)
U(r, v) = V (r) + σ · L
(a)
(b)
Putting σ in the ẑ direction it is easy (but rather lengthy) to show the result
X ∂ri
Qj = Fi · .
i
∂qj
(c)
1
T = m(ṙ 2 + r 2 θ̇2 + r 2 sin2 θφ̇2 ),
2
where the equations of motion are
d ∂T ∂T
− = Qj .
dt ∂ q̇j ∂qj
Problem 6 (1.18 of Goldstein’s book)
m k
L= (aẋ2 + 2bẋẏ + cẏ 2) − (ax2 + 2bxy + cy 2 ).
2 2
Equations of motion:
• Case 1, a = c = 0.
ÿ + ω 2 y = 0
ẍ + ω 2 x = 0
p
where ω = k/m. So we have two decoupled 2d harmonic oscillators with the same fre-
quency.
• Case 2, b = 0, c = −a.
ẍ + ω 2x = 0
ÿ + ω 2 y = 0.
the system decouples at the level of the Lagrangian. The b2 − ac 6= 0 condition is the condition
for the transformation matrix to not be singular.
θ1
m1
θ2
m2