19.01.2015 Views

My title - Departamento de Matemática da Universidade do Minho

My title - Departamento de Matemática da Universidade do Minho

My title - Departamento de Matemática da Universidade do Minho

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

7 OSCILLATIONS AND CYCLES 54<br />

Observe that if one of the bodies is much bigger than the other (like the Sun and the Earth), say<br />

m 1 ≫ m 2 , then the center of mass nearly coinci<strong>de</strong>s with the position q 1 of the bigger body, while<br />

the reduced mass m is essentially the mass m 2 of the smaller one (hence it looks like the Earth<br />

moving around the Sun, as Galileo had suggested).<br />

Central forces.<br />

Consi<strong>de</strong>r the Newton equation<br />

m¨r = F (|r|) ˆr<br />

<strong>de</strong>scribing the motion of a particle (planet) of mass m in a central force field F . Conservation<br />

of angular momentum implies that the motion is planar, hence we may take r ∈ R 2 . In polar<br />

coordinates r = ρe iθ , the equations reed<br />

¨ρ − ρ ˙θ 2 = F (ρ)/m<br />

ρ¨θ + 2 ˙ρ ˙θ = 0 .<br />

The second equation says that the “areal velocity” (“veloci<strong>da</strong><strong>de</strong> areal”) l = ρ 2 ˙θ is a constant of<br />

the motion (Kepler’s second law).<br />

Taking Newton’s gravitational force F (ρ) = − GmM<br />

ρ<br />

(where M is the mass of the Sun and G is<br />

2<br />

the gravitational constant), the first equation may be written as<br />

m¨ρ = − ∂ ∂ρ V l (ρ) ,<br />

where we <strong>de</strong>fined the ”effective potential energy” as V l (ρ) = 1 2 m l2<br />

ρ<br />

− G mM 2 ρ<br />

is<br />

E = 1 2 m ˙ρ2 + 1 2 m l2<br />

ρ 2 − GmM ρ .<br />

. The conserved energy<br />

Kepler’s effective potential and some energy level sets.<br />

Now we set ρ = 1/x and look for a differential equation for x as a function of θ. Computation<br />

shows that dx/dθ = − ˙ρ/l, and, using conservation of l, that d 2 x/dθ 2 = −ρ 2 ¨ρ/l 2 . There follows<br />

that the first Newton equation reads<br />

we get<br />

d 2 x<br />

dθ 2 + x = − 1<br />

l 2 x 2 m F (1/x) .<br />

d 2 x<br />

dθ 2 + x = −GM l 2 .

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

Saved successfully!

Ooh no, something went wrong!