04.08.2013 Views

Nonlinear Mechanics - Physics at Oregon State University

Nonlinear Mechanics - Physics at Oregon State University

Nonlinear Mechanics - Physics at Oregon State University

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

4.1. ONE-DIMENSIONAL SYSTEMS 49<br />

and<br />

∂<br />

∂ψ F1(ψ, J) = 1 [<br />

H1 − H1(ψ, J)<br />

ω0(J)<br />

] ≡ − ˜ H1(ψ, J)<br />

ω0<br />

(4.16)<br />

˜H1 is the periodic part of H1. We are left with a differential equ<strong>at</strong>ion th<strong>at</strong><br />

is easy to integr<strong>at</strong>e.<br />

F1(ψ, J) = − 1<br />

∫<br />

dψ<br />

ω0(J)<br />

˜ H1(ψ, J) (4.17)<br />

4.1.1 Summary<br />

I will summarize all these technical details in the form of an algorithm for<br />

doing first order perturb<strong>at</strong>ion theory. Remember th<strong>at</strong> the object is to find<br />

equ<strong>at</strong>ions of motion in the form q = q(t) and p = p(t). We do this in three<br />

steps: (1) Find q = q(I, ψ) and p = p(I, ψ). (2) Find I = I(J, φ) and<br />

ψ = ψ(J, φ). (3) J and ˙φ are constant, so φ = ˙φt + φ0.<br />

1. Identify the H0 part of the Hamiltonian. Find the transform<strong>at</strong>ion<br />

equ<strong>at</strong>ions q = q(I, ψ) and p = p(I, ψ) using the Hamiltonian-Jacobi<br />

equ<strong>at</strong>ion as described in the previous section. Use (??) to get ω0.<br />

2. Equ<strong>at</strong>ion (4.8) can be used to find J in terms of the total energy<br />

E. The integral presents no difficulties in principle, especially if the<br />

Hamiltonian is separable. In fact, textbooks never bother to do this.<br />

It seems sufficient to display the results in terms of J, the assumption<br />

being th<strong>at</strong> we could find J = J(E) if we really had to.<br />

3. The first order correction to the energy is obtained from the integral in<br />

(4.15). Get the first order correction to the frequency by differenti<strong>at</strong>ing<br />

it with respect to J.<br />

4. The gener<strong>at</strong>ing function F1 is calcul<strong>at</strong>ed from (4.17). It is then substituted<br />

into (4.4) and (4.5). These give implicit equ<strong>at</strong>ions for ψ =<br />

ψ(J, φ) and I = I(J, φ). Unfortun<strong>at</strong>ely, it is usually impossible to<br />

invert them to obtain these formula explicitly.<br />

4.1.2 The simple pendulum<br />

The pendulum makes a nice example<br />

H = l2<br />

+ mgR(1 − cos θ)<br />

2mR2

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!