Understanding the Hamiltonian of a Pendulum on a Riemannian Manifold, Exercises of Classical and Relativistic Mechanics

The mechanics of a pendulum on a riemannian manifold, t ∗s1, using the cotangent space t ∗ and the state of the pendulum described by a point (q(t), p(t)). The document derives the kinetic and potential energy of the pendulum, and then calculates the hamiltonian of the pendulum, h(q, p). The document also discusses hamilton's equations and their relation to the motion of the pendulum. Additionally, it includes digressions on the approximation of the pendulum as a harmonic oscillator and the significance of elliptic integrals.

Typology: Exercises

2011/2012

Uploaded on 07/19/2012

davdas900
davdas900 🇮🇳

4.6

(5)

110 documents

1 / 3

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
the cotangent space T
q(t)S1, and the state of the pendulum is described by a point (q(t), p(t)) in
TS1. However, we can treat the time derivative of an angle as a real number using the isomorphism
S1
=R/2πZ, so we get an isomorphism Tq(t)S1
=R, and thus an isomorphism T
q(t)S1
=R. We
thus have
TS1
=S1×R
and using this we can think of (q(t), p(t)) as a point in S1×R. Of course, most physicists do all this
without making such a fuss about it!
Now here’s where you come in....
1. Using what you already know about the mechanics of point particles, show that the kinetic energy
of the pendulum is
K( ˙q) = 1
2mr2˙q2
2. Assuming the force of gravity is a vector pointing down with magnitude mG, show that we can
assume the potential energy of the pendulum to be
V(q) = mG cos q.
3. Using what you know about classical mechanics on a Riemannian manifold, show that the
Hamiltonian of the pendulum,
H:TS1R,
is given by
H(q, p) = p2
2mr2mG cos q.
4. Work out Hamilton’s equations for the pendulum and show that
˙q=p
mr2,
˙p=mG sin q,
and thus
¨q=G
r2sin q.
Digression 1: Note that ˙qis really the angular velocity of the pendulum, while p=mr2˙qis really
its angular momentum.
Digression 2: If the angle qstays smal l, we can use the approximation sin q'qto approximate the
pendulum by a harmonic oscillator with
¨q=G
r2q.
However, when the angle becomes large the pendulum becomes very different from the harmonic
oscillator. For example, if the pendulum starts out at q=π, p = 0 at time zero, it will stay there for
all times, balanced upside down! This is an unstable equilibrium.
5. Plot the level curves of Has a function of (q , p)(R/2πZ)×R. Since energy is conserved, the
state (q(t), p(t)) must stay on one of these level curves as it evolves in time. Use this to qualitatively
2
docsity.com
pf3

Partial preview of the text

Download Understanding the Hamiltonian of a Pendulum on a Riemannian Manifold and more Exercises Classical and Relativistic Mechanics in PDF only on Docsity!

the cotangent space T (^) q∗(t)S^1 , and the state of the pendulum is described by a point (q(t), p(t)) in

T ∗S^1. However, we can treat the time derivative of an angle as a real number using the isomorphism S^1 ∼= R/ 2 πZ, so we get an isomorphism Tq(t)S^1 ∼= R, and thus an isomorphism T (^) q∗(t)S^1 ∼= R. We thus have T ∗S^1 ∼= S^1 × R

and using this we can think of (q(t), p(t)) as a point in S^1 × R. Of course, most physicists do all this without making such a fuss about it!

Now here’s where you come in....

  1. Using what you already know about the mechanics of point particles, show that the kinetic energy of the pendulum is

K( ˙q) =

mr^2 q˙^2

  1. Assuming the force of gravity is a vector pointing down with magnitude mG, show that we can assume the potential energy of the pendulum to be

V (q) = −mG cos q.

  1. Using what you know about classical mechanics on a Riemannian manifold, show that the Hamiltonian of the pendulum, H: T ∗S^1 → R,

is given by

H(q, p) =

p^2 2 mr^2

− mG cos q.

  1. Work out Hamilton’s equations for the pendulum and show that

q˙ = p mr^2

p ˙ = −mG sin q,

and thus

¨q = −

G

r^2

sin q.

Digression 1: Note that q˙ is really the angular velocity of the pendulum, while p = mr^2 q˙ is really its angular momentum.

Digression 2: If the angle q stays small, we can use the approximation sin q ' q to approximate the pendulum by a harmonic oscillator with

q ¨ = −

G

r^2 q.

However, when the angle becomes large the pendulum becomes very different from the harmonic oscillator. For example, if the pendulum starts out at q = π, p = 0 at time zero, it will stay there for all times, balanced upside down! This is an unstable equilibrium.

  1. Plot the level curves of H as a function of (q, p) ∈ (R/ 2 πZ) × R. Since energy is conserved, the state (q(t), p(t)) must stay on one of these level curves as it evolves in time. Use this to qualitatively

describe the behavior of the pendulum for various different values of the energy. In particular, find stable and unstable equilibria.

  1. Supposing that the pendulum’s energy equals E ∈ R, show that

q˙ = ±

mr^2

(E + mG cos q). (1)

  1. To reduce the clutter and focus on essentials, switch to units where mr^2 = mG = 1. Using equation (1), and taking the positive square root, show that

t =

dq √ 2(E + cos q)

  1. If we could do the integral in equation (2), we’d know t as a function of q. Then we could solve for q as a function of t and we’d be done! Unfortunately, this integral cannot be done using elementary functions — it’s a so-called elliptic integral. To bring it into Jacobi’s favorite form, let’s work not with q but with

x =

E + 1

sin(q/2)

Show that x˙ =

(1 − x^2 )(1 − k^2 x^2 )

where k, the so-called modulus, is given by

k =

E + 1

Conclude that

t =

dx √ (1 − x^2 )(1 − k^2 x^2 )

This is Jacobi’s elliptic integral of the first kind. When we solve for x as a function of t, we get an elliptic function.

Digression: To do integrals involving the square root of a quadratic function of x, you need inverse trig functions. However, for integrals involving the square root of a cubic or quartic function of x, you need inverse elliptic functions — or in other words, elliptic integrals. Why are they called ‘elliptic’? Well, if you work out the circumference of the ellipse

x^2 a^2

y^2 b^2

you get 4 a times this: ∫ (^1)

0

1 − k^2 x^2 1 − x^2

dx

where k^2 = 1 − b^2 /a^2. The stuff under the square root here is not a quartic in x, but the integral is closely related to the one we’ve been discussing: it’s called Jacobi’s elliptic integral of the second kind.

  1. I’ve been a bit sloppy about the limits of integration in equation (3). Show that if we start our clock so that t = 0 when our pendulum happens to be pointing straight down, we have

t =

∫ (^) x

0

dy √ (1 − y^2 )(1 − k^2 y^2 )