27.05.2013 Views

Etude de la dynamique autour des points de Lagrange

Etude de la dynamique autour des points de Lagrange

Etude de la dynamique autour des points de Lagrange

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

tel-00422422, version 1 - 6 Oct 2009<br />

Figure 2.5 – Symétrie du problème restreint <strong>de</strong>s trois corps.<br />

ce p<strong>la</strong>n intersecte à nouveau le p<strong>la</strong>n y = 0 orthogonalement, alors cette trajectoire et son image<br />

miroir se raccor<strong>de</strong>ront tant en position qu’en vitesse au niveau du p<strong>la</strong>n y = 0 et formeront donc<br />

une seule et unique orbite périodique.<br />

Considérons une condition initiale X0 appartenant au p<strong>la</strong>n y = 0 et avec un vecteur vitesse<br />

orthogonal au p<strong>la</strong>n y = 0,<br />

X0 =<br />

<br />

<br />

x0, 0, z0, 0, ˙y0, 0 ·<br />

Pour que l’orbite partant <strong>de</strong> X0 soit périodique, il faut et il suffit qu’au temps t1 où l’on recoupe<br />

le p<strong>la</strong>n y = 0 on obtienne un point <strong>de</strong> coordonnées<br />

X(t1) =<br />

<br />

<br />

x, 0, z, 0, ˙y, 0<br />

(vecteur vitesse orthogonal au p<strong>la</strong>n (y = 0)). Si cette condition est remplie, l’orbite construite<br />

sera alors périodique <strong>de</strong> pério<strong>de</strong> T = 2t1.<br />

Pour calculer numériquement une orbite périodique, on se donne une condition initiale X0<br />

<strong>de</strong> <strong>la</strong> forme précé<strong>de</strong>nte. On intègre numériquement les équations du mouvement à partir <strong>de</strong><br />

X0, et on corrige cette condition initiale jusqu’à ce que les composantes en x et z <strong>de</strong> <strong>la</strong> vitesse<br />

soient nulles au moment où <strong>la</strong> trajectoire intersecte le p<strong>la</strong>n y = 0. Numériquement, on résout ce<br />

problème en utilisant une métho<strong>de</strong> <strong>de</strong> Newton.<br />

Correction <strong>de</strong>s conditions initiales par une métho<strong>de</strong> <strong>de</strong> Newton.<br />

On cherche X0 = (x0, 0, z0, 0, ˙y0, 0) et t1 tels que<br />

⎛<br />

⎜<br />

⎝<br />

y(t1)<br />

˙x(t1)<br />

˙z(t1)<br />

⎞<br />

⎟<br />

⎠ =<br />

⎛<br />

⎜<br />

⎝<br />

φ2(X0, t1)<br />

φ4(X0, t1)<br />

φ6(X0, t1)<br />

24<br />

⎞<br />

⎟<br />

⎠ =<br />

⎛<br />

⎜<br />

⎝<br />

0<br />

0<br />

0<br />

⎞<br />

⎟<br />

⎠ ·

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

Saved successfully!

Ooh no, something went wrong!