Problème à deux corps

étoile binaire

Le problème à deux corps est un modèle théorique important en mécanique, qu'elle soit classique ou quantique, dans lequel sont étudiés les mouvements de deux corps assimilés à des points matériels en interaction mutuelle (conservative), le système global étant considéré comme isolé. Dans cet article, seul sera abordé le problème à deux corps en mécanique classique (voir par exemple l'article atome d'hydrogène pour un exemple en mécanique quantique), d'abord dans le cas général d'un potentiel attractif, puis dans le cas particulier très important où les deux corps sont en interaction gravitationnelle, ou mouvement képlérien, lequel est un sujet important de la mécanique céleste.

L'importance de ce problème vient en premier lieu de son caractère exactement intégrable, contrairement au problème à trois corps et plus. En effet le problème à deux corps, qui possède a priori six degrés de liberté, peut se ramener en fait à la résolution d'un problème à un corps à un seul degré de liberté.

Par ailleurs, les résultats obtenus permettent de rendre compte des trajectoires des planètes dans le système solaire (dans le référentiel héliocentrique) ainsi que celle de leurs satellites naturels ou artificiels, au moins en première approximation. On retrouve alors les lois de Kepler, mis en évidence par l'analyse des observations astronomiques dès le XVII^e siècle. Ainsi, la situation envisagée est loin d'être purement académique. La première résolution de ce problème a été exposée par Newton, qui a énoncé la loi fondamentale de la mécanique classique : le résultat est annoncé dans les propositions 57 à 65 de ses Principia.

Cet article a pour objet l'exposé et le traitement général du problème à deux corps, avec la démonstration des lois de Kepler et l'étude détaillée des différents types de trajectoires envisageables. La question de la détermination des éléments d'orbite ainsi que les équations de Kepler et de Barker et leurs applications font l'objet d'articles séparés (cf les articles mouvement képlérien, équation de Kepler et éléments d'orbite).

Situation envisagée et notations

Le problème à deux corps est celui de deux corps de masse m₁ et m₂, assimilés à des points matériels M₁ et M₂, respectivement, en interaction mutuelle. La force exercée par M₁ sur M₂ dérive d'un potentiel attractif V(r) et est notée ${\overrightarrow {F_{12}}}$: du fait de la troisième loi de Newton (ou principe des actions réciproques) il est évident que ${\overrightarrow {F_{21}}}=-{\overrightarrow {F_{12}}}$.

Le système global étant considéré comme isolé, il s'agit d'étudier le mouvement de M₁ et M₂ par rapport à un référentiel (R) supposé galiléen, dont le repère d'espace associé a pour origine O.

Les notations suivantes sont adoptées par la suite : ${\overrightarrow {r_{1}}}={\overrightarrow {OM_{1}}}$, ${\overrightarrow {r_{2}}}={\overrightarrow {OM_{2}}}$, et ${\overrightarrow {r}}{\stackrel {\text{déf}}{=}}{\overrightarrow {r_{12}}}={\overrightarrow {M_{1}M_{2}}}={\overrightarrow {r_{2}}}-{\overrightarrow {r_{1}}}$.

Les équations du mouvement dans (R) de chacun des corps s'écrivent alors, en utilisant la relation fondamentale de la dynamique :

${\begin{cases}m_{1}{\ddot {\overrightarrow {r_{1}}}}={\overrightarrow {F_{21}}}=-{\overrightarrow {F_{12}}}\\m_{2}{\ddot {\overrightarrow {r_{2}}}}={\overrightarrow {F_{12}}}\end{cases}}$, (1).

La stratégie de résolution du problème est alors la suivante : tout d'abord se ramener à l'étude du mouvement d'un seul corps, en introduisant la notion de particule fictive ; ensuite, se ramener à un problème à une seule dimension, facilement résoluble.

Réduction à un problème à un corps

Le fait que le système soit isolé permet de séparer le mouvement trivial de son centre d'inertie de celui d'un corps par rapport à l'autre et en fait de se ramener à l'étude du mouvement d'une seule particule, dite fictive.

Conservation de la quantité de mouvement - référentiel barycentrique

L'addition des deux équations du mouvement donne aussitôt :

$m_{1}{\ddot {\overrightarrow {r_{1}}}}+m_{2}{\ddot {\overrightarrow {r_{2}}}}=\left(m_{1}+m_{2}\right){\ddot {\overrightarrow {R_{C}}}}={\overrightarrow {0}}$, avec C centre de masse du système, de vecteur position ${\overrightarrow {R_{C}}}={\overrightarrow {OC}}={\frac {m_{1}{\overrightarrow {r_{1}}}+m_{2}{\overrightarrow {r_{2}}}}{m_{1}+m_{2}}}$.

Par suite, et comme cela était attendu pour un système isolé, le mouvement du centre de masse C dans (R) est rectiligne et uniforme (ou à la limite, au repos), et il est possible de se placer dans le référentiel du centre de masse (R_c) (qui sera galiléen, du fait du mouvement rectiligne et uniforme du corps auquel il est lié, (R) étant supposé galiléen), dit barycentrique, pour réécrire les équations du mouvement précédentes.

Introduction de la notion de particule fictive

En posant ${\overrightarrow {r}}{\stackrel {\text{déf}}{=}}{\overrightarrow {r_{2}}}-{\overrightarrow {r_{1}}}$ il est possible d'écrire :

${\begin{cases}{\overrightarrow {r_{1}}}={\overrightarrow {R}}_{C}-{\frac {m_{2}}{m_{1}+m_{2}}}{\overrightarrow {r}}\\{\overrightarrow {r_{2}}}={\overrightarrow {R}}_{C}+{\frac {m_{1}}{m_{1}+m_{2}}}{\overrightarrow {r}}\end{cases}}$, (2).

Il suffit de prendre la différence entre les deux équations du mouvement (1), et de tenir compte du caractère galiléen du référentiel barycentrique qui implique ${\ddot {{\vec {R}}_{c}}}={\vec {0}}$, pour obtenir :

${\ddot {\vec {r}}}=\left({\frac {1}{m_{1}}}+{\frac {1}{m_{2}}}\right){\vec {F}}_{12}=\left({\frac {m_{1}+m_{2}}{m_{1}m_{2}}}\right){\overrightarrow {F_{12}}}$.

Cette équation est en fait celle du mouvement d'un seul corps, à trois degrés de liberté :

$\mu {\ddot {\overrightarrow {r}}}={\overrightarrow {F}}$,

avec $\mu {\stackrel {\text{déf}}{=}}{\frac {m_{1}m_{2}}{m_{1}+m_{2}}}$, masse réduite du système, et ${\overrightarrow {F}}={\overrightarrow {F_{12}}}$.

Dans le référentiel barycentrique le problème se ramène donc au mouvement d'une particule dite fictive, de masse μ et de rayon-vecteur ${\overrightarrow {r}}$, les trajectoires des corps M₁ et M₂ se déduisant par homothétie, d'après les formules précédentes sur ${\vec {r}}_{1,2}$ .

Il convient de remarquer que dans le cas particulier important où l'un des corps a une masse beaucoup plus importante que le second (corps central, généralement une étoile, ou une « grosse » planète), par exemple si $m_{1}\gg m_{2}$, le centre de masse du système est pratiquement confondu avec ce corps central, et la masse réduite est pratiquement égale à celle de l'autre corps, $\mu \simeq m_{2}$. À noter toutefois que pour le mouvement de la Lune, qui dans le système solaire possède la plus forte masse relative d'un satellite par rapport à sa planète (1/81 Mt), cette approximation est relativement peu précise.

Intégrale première du moment cinétique - Planéité de la trajectoire - Loi des aires

Dans le cas particulier très important d'une force centrale, on a ${\overrightarrow {F}}=F(r){\overrightarrow {e_{r}}}$ avec ${\overrightarrow {e_{r}}}{\stackrel {\text{déf}}{=}}{\dfrac {\overrightarrow {r}}{r}}$, le théorème du moment cinétique au centre de force noté O s'écrit :

${\dot {\overrightarrow {L}}}={\overrightarrow {r}}\times {\overrightarrow {F}}={\overrightarrow {0}}$,

ce qui implique ${\overrightarrow {L}}={\overrightarrow {r}}\times \left(\mu {\dot {\overrightarrow {r}}}\right)={\overrightarrow {\mathrm {cte} }}$.

Physiquement, cela impose que le vecteur position et le vecteur vitesse de la particule fictive soient à tout instant perpendiculaires à un vecteur constant: la trajectoire de M est donc plane, le problème est donc à deux degrés de liberté.

Dans le plan de la trajectoire, défini comme celui généré par ${\overrightarrow {r}}\left(t=0\right)={\overrightarrow {r}}_{0}$ et ${\dot {\overrightarrow {r}}}\left(t=0\right)={\overrightarrow {v}}_{0}$, il est judicieux de se placer en coordonnées cylindro-polaires d'axe la direction Oz de ${\overrightarrow {L}}=L{\overrightarrow {e_{z}}}$, avec θ angle entre ${\overrightarrow {r}}_{0}$ et ${\overrightarrow {r}}$, il vient :

${\overrightarrow {L}}=\left(r{\overrightarrow {e_{r}}}\right)\times \mu \left({\dot {r}}{\overrightarrow {e_{r}}}+\left(r{\dot {\theta }}\right){\overrightarrow {e_{\theta }}}\right)=\mu r^{2}{\dot {\theta }}{\overrightarrow {e_{z}}}$,

par suite:

$\mu r^{2}{\dot {\theta }}=L=\mathrm {cte}$(3).

Or l'aire élémentaire balayée par le rayon-vecteur ${\overrightarrow {r}}$ pendant dt est donnée par :

$\mathrm {d} {\mathcal {A}}={\dfrac {1}{2}}r^{2}\mathrm {d} \theta$.

La vitesse aréolaire est donc constante pour la particule fictive (il en est de même par homothétie pour les corps réels) :

${\dot {\mathcal {A}}}={\frac {\mathrm {d} {\mathcal {A}}}{\mathrm {d} t}}={\frac {1}{2}}C=\mathrm {cte}$,

avec $C{\stackrel {\text{déf}}{=}}{\frac {L}{\mu }}$, la constante des aires.

Par suite, le rayon vecteur de chaque particule balaie des aires égales en des temps égaux. Cette propriété est en fait valable pour tout mouvement à force centrale. De par l'expression de C, il est facile de voir que la vitesse angulaire de la particule fictive est inversement proportionnelle à la distance r, et est donc maximale quand celle-ci est minimale, soit au périastre dans le mouvement képlérien - cf. infra.

Remarques :

• Un mouvement circulaire uniforme est un cas simple de mouvement obéissant à la loi des aires ;
• L'expression de L implique que la vitesse angulaire instantanée ${\dot {\theta }}$ ne change jamais de signe le long de la trajectoire : en quelque sorte, la particule fictive, et donc les deux corps "réels" "tournent" toujours dans le "même sens" sur leurs trajectoires ;
• Si L=0, le mouvement est rectiligne : ce cas particulier (dit dégénéré) sera exclu par la suite.

Intégrale première de l'énergie - potentiel effectif

Le mouvement étant conservatif puisque la force dérive d'une énergie potentielle V(r), l'énergie totale est une intégrale première du mouvement :

$H={\frac {1}{2}}\mu {\dot {\overrightarrow {r}}}^{2}+V(r)=\mathrm {cte}$, or ${\dot {\overrightarrow {r}}}^{2}={\dot {r}}^{2}+r^{2}{\dot {\theta }}^{2}$, soit en utilisant l'expression de la valeur du moment cinétique L :

$H={\frac {1}{2}}\mu {\dot {r}}^{2}+U_{\text{eff}}\left(r\right)$, (4),

avec $U_{\text{eff}}\left(r\right){\stackrel {\text{déf}}{=}}{\frac {L^{2}}{2\mu r^{2}}}+V(r)$ potentiel effectif.

Finalement le problème se ramène à l'étude du mouvement d'un seul corps à un seul degré de liberté r. Ceci est toujours valable dans le cas du problème à deux corps, quelle que soit la nature du potentiel d'interaction.

Résolution analytique

D'après (4), il est possible d'exprimer la vitesse radiale ${\dot {r}}$, il vient : ${\dot {r}}={\frac {\mathrm {d} r}{\mathrm {d} t}}={\sqrt {{\frac {2}{\mu }}\left(H-V(r)\right)-{\frac {L^{2}}{\mu ^{2}r^{2}}}}}$.

Il est alors possible de séparer les variables, et d'intégrer entre deux instants t₀ et t, auxquels correspondent respectivement les positions radiales r₀ et r, pour obtenir :

$t-t_{0}=\int _{r_{0}}^{r}{\frac {\mathrm {d} r'}{\sqrt {{\frac {2}{\mu }}\left(H-V(r')\right)-{\frac {L^{2}}{\mu ^{2}r'^{2}}}}}}$, (4bis).

Ceci correspond, de façon implicite, à l'équation horaire $r(t)$.

Compte tenu de (4), il est alors possible d'obtenir une expression similaire pour θ :

$\theta -\theta _{0}=\int _{r_{0}}^{r}{\frac {{\frac {L}{r'^{2}}}\mathrm {d} r'}{\sqrt {2\mu \left(H-V(r')\right)-{\frac {L^{2}}{r'^{2}}}}}}$, (4ter).

Ces deux expressions sont, en pratique, difficiles à utiliser. Elles permettent toutefois de discuter qualitativement de la nature des mouvements possibles.

Étude qualitative des mouvements possibles

Dans l'expression de $U_{\text{eff}}\left(r\right)$, le terme ${\frac {L^{2}}{2\mu r^{2}}}$, toujours positif, correspond à une barrière centrifuge. Le potentiel V(r) étant supposé :

• attractif : on a donc V(r) monotone croissant pour tout r ;
• régulier en l'origine : c'est-à-dire $\lim _{r\to 0^{+}}r^{2}\left|V(r)\right|=0$, et donc le terme en ¹⁄_r² domine quand $r\to 0^{+}$.

Bien que cela ne soit pas indispensable, dans la suite V(r) sera également supposé borné à l'infini, soit avec un choix judicieux de l'origine des potentiels, $\lim _{r\to \infty }V(r)=0^{+}$. Par convention V(r) < 0.

Avec ces conditions, valables pour les potentiels physiques usuels, le potentiel effectif présente un minimum absolu unique, noté $-U_{\min }<0$ pour $r=r_{m}$ tel que $\left.{\frac {\mathrm {d} U_{\text{eff}}}{\mathrm {d} r}}\right|_{r=r_{m}}=0$, et possède donc une cuvette de potentiel (cf. figure ci-contre, avec un potentiel newtonien).

Allure de la courbe du potentiel effectif U_eff(r) avec les contributions des différents termes

Par ailleurs, on a que, d'après l'expression de H : $H-U_{\text{eff}}(r)={\frac {1}{2}}\mu {\dot {r}}^{2}\geqslant 0$, les valeurs de r permises doivent être telles que $U_{\text{eff}}(r)\leqslant H$.

Il est dès lors possible d'envisager qualitativement les cas suivants (L étant non nul) :

• Énergie totale H>0 : puisque $U_{\text{eff}}\left(r\right)\longrightarrow 0^{-}$ quand $r\to \infty$, cette condition est vérifiée pour toute valeur de $r\geqslant r_{\min }$ avec r_min distance minimale d'approche telle que $U_{\text{eff}}\left(r_{\min }\right)=H$. La particule fictive peut donc aller à l'infini avec une vitesse radiale positive, égale à ${\dot {r}}_{\max }={\sqrt {\dfrac {2H}{\mu }}}$.

• Énergie totale H = 0 : cas limite où la particule peut aller à l'infini, mais avec une vitesse radiale nulle, d'après la formule précédente, la distance minimale d'approche étant telle que $U_{\text{eff}}\left(r_{\min }\right)=0$, soit $V\left(r_{\min }\right)=-{\frac {L^{2}}{2\mu r_{\min }^{2}}}$.

• Énergie totale $0>H>-U_{\min }$ : dans ce cas, la particule est confinée dans une région précise de l'espace, comprise entre les deux valeurs de r telles que $U_{\text{eff}}\left(r_{\min ,\max }\right)=H<0$ (points d'arrêt). La vitesse radiale s'annule à chacun de ces points, toutefois ce n'est pas le cas de la vitesse angulaire instantanée (cf. expression de H et remarques à son sujet).

Du fait du confinement du mouvement, la durée Δt pour que r varie de r_min à r_max peut être obtenue aisément en utilisant l'expression intégrale donnant t. Comme l'invariance dans le temps du hamiltonien implique que cette durée sera la même pour aller de r_max à r_min, le mouvement radial est donc périodique de période T donnée par :

$T=2\int _{r_{\min }}^{r_{\max }}{\frac {\mathrm {d} r}{\sqrt {{\frac {2}{\mu }}\left(H-V(r)\right)-{\frac {L^{2}}{\mu ^{2}r^{2}}}}}}$

, et en une période radiale où l'angle θ varie de la quantité Δθ, donnée par l'expression intégrale de θ :

$\Delta \theta =2\int _{r_{\text{min}}}^{r_{\text{max}}}{\frac {{\frac {L}{r^{2}}}\mathrm {d} r}{\sqrt {2\mu \left(H-V(r)\right)-{\frac {L^{2}}{r^{2}}}}}}$.

Mouvements permis dans le potentiel effectif U_eff(r) en fonction des valeurs de l'énergie totale H, cf. discussion dans le texte.

Toutefois, il est important de souligner que le confinement du mouvement n'implique nullement que la trajectoire du mobile soit une courbe fermée. Pour cela, il faudrait en effet que $\Delta \theta ={\frac {2\pi m}{n}}$ avec m et n entiers. Dans ce cas, et seulement dans ce cas, le rayon vecteur ${\overrightarrow {r}}$ retrouve sa valeur initiale après n périodes « radiales » T, puisque qu'alors il aura « tourné » de 2mp : nT sera en fait la période du mouvement et de la fonction θ(t). Une telle situation correspond à des périodes radiale et angulaire qui sont commensurables, et il s'agit d'une condition nécessaire et suffisante pour que le mouvement borné s'effectue selon une courbe fermée. Cette condition n'est remplie que par le potentiel newtonien en 1/r et le potentiel harmonique spatial isotrope V(r)=kr² (ce dernier cas ne sera pas considéré plus avant) : ce résultat constitue le théorème de Bertrand.

• Énergie totale $H=-U_{\min }<0$ : il s'agit du cas limite où la seule valeur permise de r est r_m, et dans ce cas la trajectoire est un cercle de ce rayon.

• Énergie totale $H<-U_{\min }$ : ces valeurs de r sont interdites et ne correspondent à rien de physique.

Il sera montré plus bas que chacun de ces cas correspond aux formes particulières de trajectoires pour le mouvement képlérien, à savoir une hyperbole, une parabole, une ellipse et un cercle. La discussion précédente peut être résumée graphiquement sur la figure ci-contre.

Dégénérescence du mouvement

L'étude précédente a été faite en supposant $L\neq 0$. Si L = 0, on a tout simplement ${\dot {\theta }}=0$ à tout moment, et le mouvement est purement radial : il est dit dégénéré. La discussion précédente se simplifie, la condition (4bis) précédente se résume à $V(r)\leqslant H$, vérifiée dans tous les cas si $H\geqslant 0$. Dans le cas contraire, il est facile de vérifier que la particule "tombe" sur le centre de force.

Cas du mouvement képlérien

Le mouvement képlérien correspond au cas où les deux corps sont en interaction gravitationnelle, c'est-à-dire avec le potentiel d'interaction $V(r)=-{\frac {Gm_{1}m_{2}}{r}}$, et donc ${\overrightarrow {F}}=-G{\frac {m_{1}m_{2}}{r^{2}}}{\overrightarrow {e_{r}}}=-G{\frac {\mu M_{T}}{r^{2}}}{\overrightarrow {e_{r}}}$, avec $M_{T}{\stackrel {\text{déf}}{=}}m_{1}+m_{2}$ masse totale du système. Tout se passe alors comme si la particule fictive M subissait l'interaction gravitationnelle d'un corps affecté de la masse totale du système placé à l'origine O du rayon-vecteur. Si les résultats généraux précédents, valables d'ailleurs pour tout mouvement dans un potentiel central conservatif V(r), permettraient d'ores et déjà de déterminer l'équation horaire de r = r(t), le potentiel newtonien possède une intégrale première particulière, le vecteur de Runge-Lenz, qui permet d'obtenir de façon simple l'équation de la trajectoire.

Invariant de Runge-Lenz - équations de la trajectoire

Existence d'une intégrale première additionnelle

Le potentiel newtonien en $1/r$ se caractérise par l'existence d'un invariant supplémentaire particulier, l'invariant de Runge-Lenz, donné par:

${\overrightarrow {e}}={\frac {\left({\dot {\overrightarrow {r}}}\times {\overrightarrow {L}}\right)}{GM_{T}\mu }}-{\overrightarrow {e_{r}}}$, (5)

Démonstration

${\ddot {\overrightarrow {r}}}\times {\overrightarrow {L}}={\frac {\mathrm {d} \left({\dot {\overrightarrow {r}}}\times {\overrightarrow {L}}\right)}{\mathrm {d} t}}={\ddot {\overrightarrow {r}}}\times \left[{\overrightarrow {r}}\times \left(\mu {\dot {\overrightarrow {r}}}\right)\right]$, où il a été tenu compte de ${\overrightarrow {L}}={\overrightarrow {\mathrm {cte} }}$, ce qui donne :

${\frac {\mathrm {d} \left({\dot {\overrightarrow {r}}}\times {\overrightarrow {L}}\right)}{\mathrm {d} t}}=-{\frac {GM_{T}{\overrightarrow {e_{r}}}}{r^{2}}}\times \left(\mu r^{2}{\dot {\theta }}{\overrightarrow {e_{z}}}\right)=GM_{T}\mu {\dot {\theta }}{\overrightarrow {e_{\theta }}}$, où il a été tenu compte de la planéité de la trajectoire, or compte tenu de ${\dot {\theta }}{\overrightarrow {e_{\theta }}}={\dot {\overrightarrow {e_{r}}}}$ il vient aussitôt :

${\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {\left({\dot {\overrightarrow {r}}}\times {\overrightarrow {L}}\right)}{GM_{T}\mu }}-{\overrightarrow {e_{r}}}\right)={\overrightarrow {0}}$, d'où le résultat.

Équation de la trajectoire de la particule fictive

Il est évident que ${\overrightarrow {L}}\cdot {\overrightarrow {e}}=0$ donc ${\overrightarrow {e}}$ est contenu dans le plan du mouvement. Par suite, il est possible de prendre comme angle polaire l'angle w entre ${\overrightarrow {e}}$ et ${\overrightarrow {r}}$, avec bien sûr ${\dot {\theta }}={\dot {w}}$, et en notant e la norme de ${\overrightarrow {e}}$ il est facile de vérifier :

${\overrightarrow {r}}\cdot {\overrightarrow {e}}=re\cos {w}={\frac {{\overrightarrow {r}}\cdot \left({\dot {\overrightarrow {r}}}\times {\overrightarrow {L}}\right)}{GM_{T}\mu }}-r$, soit compte tenu de l'identité

${\overrightarrow {r}}\cdot \left({\dot {\overrightarrow {r}}}\times {\overrightarrow {L}}\right)={\overrightarrow {L}}\cdot \left({\overrightarrow {r}}\times {\dot {\overrightarrow {r}}}\right)={\frac {L^{2}}{\mu }}$ :

$r\left(1+e\cos {w}\right)=p$, avec $p{\stackrel {\text{déf}}{=}}{\frac {L^{2}}{GM_{T}\mu ^{2}}}$.

Physiquement, p = r_m, valeur de r pour laquelle U_eff(r) est minimal. De fait puisque $U_{\text{eff}}(r)={\frac {L^{2}}{2\mu r^{2}}}-{\frac {GM_{T}\mu }{r}}$, il vient facilement par dérivation à ${\frac {L^{2}}{\mu r_{m}}}=GM_{T}\mu$ soit en effet r_m = p.

L'équation de la trajectoire obtenu est donc celle d'une conique d'excentricité e et de paramètre p :

$r={\frac {p}{\left(1+e\cos {w}\right)}}$, (6),

dont le centre de force occupe l'un des foyers. Le vecteur ${\overrightarrow {e}}$ est donc dirigé vers le point de distance minimale (ou périastre), distance notée q avec $q=p/\left(1+e\right)$, correspondant à w = 0. En ce point, la vitesse angulaire orbitale ${\dot {w_{0}}}={\dot {\theta _{0}}}$ est maximale. L'angle w est appelé anomalie vraie en astronomie.

Par homothétie, chacun des corps réels décrit dans le référentiel barycentrique une conique dont le centre de masse occupe l'un des foyers. La valeur de e détermine la nature de la conique :

• Si e > 1, la trajectoire est une hyperbole : dans ce cas les corps célestes ne sont pas liés et peuvent aller à l'infini ;

• Si e = 1, la trajectoire est une parabole ;

• Si 0 < e < 1, la trajectoire est une ellipse. C'est le cas des planètes et de la plupart des autres corps du système solaire, où par surcroît le centre de masse est pratiquement confondu avec la position du Soleil, ce qui permet de retrouver la première loi de Kepler (1609) : « Au cours de leur mouvement autour du Soleil, les planètes décrivent des ellipses dont le Soleil occupe l'un des foyers ». La seconde loi (1609 également) découle directement de la constance de la vitesse aréolaire, et porte d'ailleurs le nom de loi des aires : « Le rayon-vecteur reliant le Soleil à une planète balaie des aires égales en des temps égaux ».

• Si e=0, la trajectoire est circulaire.

Remarque sur le caractère particulier du champ newtonien

L'obtention de trajectoires fermées pour les champs newtoniens est remarquable, et découle en fait d'une symétrie particulière. En effet, d'après le théorème de Noether, les lois de conservation sont liées à l'existence de symétrie particulière du problème.

Ainsi l'invariance par translation du système global des deux corps (liée au caractère supposé isolé du système) conduit à la conservation de la quantité de mouvement du système global, tandis que l'invariance par rotation (isotropie) du champ central conduit à celle du moment cinétique et l'invariance par translation dans le temps (qui suppose l'absence de « frottements ») à la conservation de l'énergie totale. L'existence de ces intégrales premières permet de passer successivement de 6 à 3 puis 2 et enfin un degré de liberté. Toutefois, même pour le mouvement fini, les deux corps n'ont aucune raison de décrire des courbes fermées, et seule une symétrie additionnelle y conduit, et se traduit par l'existence de l'intégrale première particulière du champ en 1/r², le vecteur de Runge-Lenz (ou, de façon équivalente, le vecteur excentricité).

Il est clair d'après ce qui précède que l'équation de la trajectoire s'obtient simplement par projection du rayon-vecteur sur ce vecteur constant particulier qui définit une direction particulière de l'espace, en pratique l'axe de la conique. Ce type de symétrie ne peut cependant s'interpréter correctement qu'à quatre dimensions, cf. article particulier.

En mécanique quantique, cet observable supplémentaire se retrouve dans l'étude de l'atome d'hydrogène, qui correspond à un problème à deux corps quantique, et ceci se traduit par la dégénérescence « accidentelle » des niveaux d'énergie de l'électron, qui ne dépendent pas du nombre quantique orbital l lié au moment cinétique. Là encore, une symétrie additionnelle spécifique au champ coulombien, et interprétable seulement à 4 dimensions, explique ce phénomène.

En creux, ces considérations montrent que dans un cas où il est tenu compte des perturbations dues aux autres corps célestes, le potentiel subi ne sera plus en 1/r² et le vecteur de Runge-Lenz ne sera plus strictement une intégrale première du mouvement. Les trajectoires ne seront plus rigoureusement des coniques, courbes fermées (pour 0 < e < 1). C'est de fait ce qui est observé, avec des orbites elliptiques qui « tournent » lentement dans l'espace, phénomène d'avance du périhélie, qui en toute rigueur s'interprète dans le cadre de la relativité générale.

Relation entre l'excentricité et les intégrales premières du mouvement

Dans un repère d'espace Oxyz fixe lié à (R_c), l'axe Ox de la conique, orienté vers le point de distance minimale au foyer de vecteur unitaire ${\overrightarrow {e_{x}}}$, Oz la direction du moment cinétique ${\overrightarrow {L}}$, et comme au périastre ${\dot {r}}=0$ et r = q, l'intégrale première de l'énergie peut être exprimée en fonction de cette distance minimale, il vient :

$H={\frac {L^{2}}{2\mu q^{2}}}-{\frac {GM_{T}\mu }{q}}$,

ce qui implique $L^{2}=2\mu q^{2}H+2GM_{T}\mu ^{2}q$.

De même, il est possible d'exprimer le vecteur excentricité de façon simple :

${\overrightarrow {e}}=\left({\frac {L^{2}}{GM_{T}\mu ^{2}q}}-1\right){\overrightarrow {e_{x}}}$.

Démonstration

${\overrightarrow {e}}={\frac {\left(q{\dot {\theta _{0}}}{\overrightarrow {e_{y}}}\right)\times \left(L{\overrightarrow {e_{z}}}\right)}{GM_{T}\mu }}-{\overrightarrow {e_{x}}}=\left({\frac {r_{\min }{\dot {\theta _{0}}}L}{GM_{T}\mu }}-1\right){\overrightarrow {e_{x}}}$, or puisque $q{\dot {\theta _{0}}}={\frac {L}{\mu q}}$, d'où le résultat.

En remplaçant dans cette expression celle obtenue pour L², l'expression de l'excentricité de la conique se met sous la forme :

$e=1+{\frac {2qH}{GM_{T}\mu }}$.

Il est alors possible d'éliminer la distance minimale q dans cette expression compte tenu de $q={\frac {p}{1+e}}={\frac {L^{2}}{G^{2}M_{T}^{2}\mu ^{3}(1+e)}}$, et d'obtenir une relation entre l'excentricité e et les deux intégrales premières H et L. Il vient par substitution dans l'expression précédente de l'excentricité :

$(e+1)(e-1)={\frac {2L^{2}H}{G^{2}M_{T}^{2}\mu ^{3}}}$, soit finalement : $e={\sqrt {1+{\frac {2L^{2}H}{G^{2}M_{T}^{2}\mu ^{3}}}}}$.

Ces deux dernières expressions n'ont bien sûr de sens physique que si $e\geqslant 0$, et permettent de retrouver les différents cas vus plus haut :

• Cas de la trajectoire hyperbolique : e > 1 ce qui implique H > 0, comme il a été dit le mobile fictif peut aller en l'infini, avec une vitesse purement radiale non nulle, donnée par ${\dot {r}}_{\max }={\sqrt {\dfrac {2H}{\mu }}}$ ;

• Cas de la trajectoire parabolique : e = 1 implique H = 0, là-aussi le mobile peut aller en l'infini avec une vitesse nulle.

• Cas de la trajectoire elliptique : 0 < e < 1 implique que l'on a $-{\frac {GM_{T}\mu }{2q}}<H<0$ .

• Cas de la trajectoire circulaire : e = 0 donc $H=-{\frac {GM_{T}\mu }{2q}}=-{\frac {G^{2}M_{T}^{2}\mu ^{3}}{2L^{2}}}$.

Remarque : dans le cas d'un engin spatial en orbite, capable de modifier son énergie cinétique et donc H et L, l'expression précédente de e montre qu'il est possible, en « choisissant » correctement l'impulsion correctrice, de modifier la valeur de l'excentricité e de la trajectoire : ceci est couramment employé pour amener les satellites vers les orbites désirées en utilisant des "orbites de transfert" appropriées. Il est même possible le cas échéant d'échapper à l'attraction terrestre, pour des sondes interplanétaires par exemple.

Cas du mouvement elliptique

Le mouvement képlérien elliptique est très important pour l'astronomie (orbites des planètes dans le système solaire, satellites artificiels, par exemple). Il sert en tout cas de point de départ à des calculs plus avancés, la prise en compte de l'influence d'autres corps célestes ou d'autres facteurs se faisant le plus souvent comme perturbation de ce mouvement « idéal ».

Principaux paramètres de la trajectoire

Une ellipse peut se définir comme le lieu des points M tels que la somme des distances à deux points fixes appelés foyers notés F₁ et F₂ est constante : $d(F_{1},M)+d(F_{2},M)=2a$, a étant le demi-grand axe de l'ellipse, c'est-à-dire la demi-distance entre les deux sommets que sont pour une trajectoire le périastre P et l'apoastre A, point le plus éloigné du centre de force (qui occupe l'un des foyers) de l'orbite (cf. figure).

Principaux paramètres de la trajectoire elliptique, cf. discussion dans le texte.

En notant Q la distance au foyer de ce dernier point, qui correspond à $w=\pi$, il vient : $Q=p/\left(1-e\right)$, et on en déduit $a=p/\left(1-e^{2}\right)$. L'ellipse possède un centre de symétrie O à mi-distance des foyers, par lequel passent ses deux axes de symétrie perpendiculaires, dits grand et petit axe.

Le paramètre p de l'ellipse correspond lui à la valeur de r pour w = p/2.

La demi-distance au foyer est notée c, on a de façon évidente $c=a-q=e{\frac {p}{1-e^{2}}}=ea$. On déduit le demi-petit axe, noté b : $b={\sqrt {a^{2}-c^{2}}}=a{\sqrt {1-e^{2}}}$. On peut éliminer e entre les équations donnant a et b, il vient l'expression du paramètre de l'ellipse : $p=b^{2}/a$.

En combinant ces formules, on obtient respectivement pour le périastre q et l'apoastre Q : $q=a(1-e)$ et $Q=a(1+e)$.

Aspects énergétiques - Équation des forces vives

Pour obtenir l'équation qui tire son nom, « forces vives », du concept leibnizien de la force, on remplace l'expression de la distance au périastre $q=a\left(1-e\right)$ dans la relation (9) excentricité-énergie, il vient :

$1-e=-{\frac {2qH}{GM_{T}\mu }}=-{\frac {2Ha\left(1-e\right)}{GM_{T}\mu }}$, d'où l'on tire l'expression de l'énergie mécanique totale en fonction de la valeur du demi-grand axe a :

$H=-{\frac {GM_{T}\mu }{2a}}$, (10bis).

Par suite, l'énergie totale est égale à la moitié de l'énergie potentielle de gravitation pour r = a. Comme H = cte, il vient :

$H=-{\frac {GM_{T}\mu }{2a}}={\frac {1}{2}}\mu v^{2}-{\frac {GM_{T}\mu }{r}}$, d'où l'on tire l'expression de la valeur de la vitesse v de la particule sur la trajectoire en fonction de r et a, dite équation des forces vives :

$v={\sqrt {GM_{T}\left({\frac {2}{r}}-{\frac {1}{a}}\right)}}$, (10ter).

Cette équation est d'usage fréquent en astronautique. On constate donc que la mesure du demi-grand axe est directement liée à l'énergie totale de la particule fictive, et que si $\left|H\right|\to 0$ on a $a\to \infty$, comme on doit s'y attendre puisque la trajectoire elliptique tend alors vers une trajectoire parabolique.

Troisième loi de Kepler

Le mouvement elliptique étant fini, il est périodique de période T, dite période de révolution ou période orbitale, les lois physiques étant invariantes par translation dans le temps. Or cette période est facile à mesurer pour un corps céleste donné par des observations astronomiques, comme l'est également la valeur du demi-grand axe a. Or il existe une relation simple entre période de révolution et demi-grand axe, mise en évidence expérimentalement pour la première fois par Kepler en 1618.

La vitesse aréolaire étant constante avec ${\dot {\mathcal {A}}}={\frac {L}{2\mu }}$, on obtient en intégrant sur la période T du mouvement l'aire totale de l'ellipse, égale à $\pi ab$, d'où l'identité $\pi ab={\frac {LT}{2\mu }}$, qui donne $L^{2}={\frac {4\pi ^{2}\mu ^{2}a^{2}b^{2}}{T^{2}}}=4\pi ^{2}\mu ^{2}p{\frac {a^{3}}{T^{2}}}$, où il a été fait usage de $p=b^{2}/a$.

Or par définition, $p={\frac {L^{2}}{GM_{T}\mu ^{2}}}$ ce qui donne finalement en éliminant le moment cinétique la relation fondamentale : ${\frac {a^{3}}{T^{2}}}={\frac {GM_{T}}{4\pi ^{2}}}$.

Autrement dit : les cubes des demi-grands axes des orbites sont proportionnels aux carrés des périodes de révolution.

Dans le cas d'une planète du système solaire, la masse du Soleil est pratiquement égale à la masse totale du système, et on écrit ${\frac {a^{3}}{T^{2}}}={\frac {k}{4\pi ^{2}}}$, avec $k{\stackrel {\text{déf}}{=}}GM_{S}$ constante de Gauss.

En faisant cette approximation, on peut écrire pour deux planètes quelconques du système solaire, avec des notations évidentes :

${\frac {a_{1}^{3}}{T_{1}^{2}}}={\frac {a_{2}^{3}}{T_{2}^{2}}}$.

Ainsi la connaissance du demi-grand axe d'une planète (par exemple la Terre, que l'on peut mesurer très précisément par des méthodes astrométriques) et des périodes orbitales (par des observations) permet de déterminer les demi-grands axes de toutes les planètes du système solaire. On peut procéder de même pour tous les autres "systèmes" (e.g. une étoile et ses planètes, une planète et ses satellites...) où l'on peut négliger la masse d'un corps donné devant celle du corps "central" (étoile, planète).

On appelle moyen mouvement, noté n, la vitesse angulaire moyenne du corps céleste sur son orbite : $n{\stackrel {\text{déf}}{=}}{\frac {2\pi }{T}}$, on a donc d'après la troisième loi de Kepler :

$n={\sqrt {\frac {GM_{T}}{a^{3}}}}$.

Cas limite de la trajectoire circulaire - vitesse de satellisation minimale

Quand e = 0, l'ellipse est dégénérée et se ramène à un cercle de rayon R = p et de centre O, confondu avec les deux foyers. Tout axe passant par le centre est axe de symétrie de la trajectoire. D'après (3) la vitesse angulaire ${\dot {\theta }}=L/\mu p^{2}$ est constante (on retrouve donc bien une « loi des aires » triviale).

Sur le plan énergétique, d'après la formule (10) et comme il a été dit précédemment (cf. partie 3.1-4, cela correspond au minimum physique de l'énergie totale H avec $H=-{\frac {GM_{T}\mu }{2R}}=-{\frac {G^{2}M_{T}^{2}\mu ^{3}}{2L^{2}}}$, donc une énergie cinétique radiale ${\dfrac {1}{2}}\mu {\dot {r}}^{2}$ nulle. L'énergie potentielle effective $U_{\text{eff}}(r)$ est donc minimale, ce qui correspond physiquement au point où les termes de barrière centrifuge en 1/r² et du potentiel attractif "centripète" en 1/r s'équilibrent.

Ce n'est que pour cette énergie mécanique minimale que l'on peut avoir satellisation du mobile fictif autour du centre de force à la distance R (et donc également d'un corps "réel" autour de l'autre), une valeur plus basse conduit à la « retombée » de la particule fictive sur le centre de force. À cette énergie mécanique minimale correspond une valeur particulière de la vitesse v dite vitesse de satellisation minimale ou « première vitesse cosmique ». En effet, d'après l'équation des forces vives (10ter), on pour la trajectoire circulaire de rayon R = a une vitesse purement orthoradiale et de valeur constante donnée par :

$v_{1}={\sqrt {\frac {GM_{T}}{R}}}$, (11).

Dans les applications (astronautique), cette vitesse est en fait indépendante de la masse du corps « satellite », puisque la masse totale équivaut alors à celle de l'astre « central ».

Exemple : Pour la Terre, on a au minimum R = 6 400 km (surface de la Terre), ce qui donne environ v₁ ≈ 7,9 km s⁻¹.

Cas du mouvement parabolique

Le mouvement parabolique est un cas limite du mouvement elliptique quand l'excentricité e tend vers 1. Intuitivement, cela correspond à une ellipse de plus en plus allongée, le périastre P se rapprochant du foyer F₁, l'autre foyer F₂ se retrouvant lui « projeté » de plus en plus loin. À la limite il est rejeté à l'infini tout comme l'apoastre A, et l'ellipse « s'ouvre » au niveau du point A pour donner une parabole.

Principaux paramètres de la trajectoire parabolique, cf. discussion dans le texte

Principaux paramètres de la trajectoire

Ce cas correspond à e = 1, et l'équation polaire au foyer de la trajectoire s'écrit alors : $r={\frac {p}{1+\cos w}}$. Le périastre P correspond à w = 0 et est situé à la distance q = p/2 du foyer F, et on a $r\to +\infty$ pour $w\to \pi$. La direction (FP) est axe de symétrie de la courbe, et il n'existe par d'apocentre à distance finie.

La figure ci-contre récapitule les principales caractéristiques de la trajectoire.

Aspect énergétiques - vitesse de libération

D'après la relation (10) et comme indiqué précédemment, ce cas limite correspond à H=0. Dans ce cas, l'énergie cinétique totale est toujours égale à l'énergie potentielle de gravitation, soit à r donné : ${\frac {1}{2}}\mu v^{2}={\frac {GM_{T}\mu }{r}}$, il vient aussitôt une expression simple de la vitesse sur la trajectoire à r quelconque, qui correspond à l'équation des forces vives pour le mouvement képlerien parabolique : $v={\sqrt {\frac {2GM_{T}}{r}}}$, (12).

Cette relation correspond à l'équation dite des forces vives trouvée pour le mouvement elliptique, relation (10ter), avec $a\to \infty$.

Cette valeur de la vitesse est maximale au périastre situé à la distance r = q, où la vitesse est purement orthoradiale. On définit dès lors pour une distance R donnée la vitesse de libération ou "seconde vitesse cosmique" comme $v_{2}={\sqrt {\frac {2GM_{T}}{R}}}={\sqrt {2}}v_{1}$.

Il s'agit alors de la vitesse minimale qu'il faut imprimer (de façon orthoradiale) à la particule fictive placée à la distance R donnée du centre de force pour que celle-ci puisse « échapper » à l'attraction gravitationnelle exercée par celui-ci, en lui permettant d'aller à l'infini, suivant une trajectoire parabolique.

Exemple : Pour la Terre depuis sa surface, on a v₂ ≈ 11,2 km s⁻¹.

Concrètement, il est possible de donner à un objet tel qu'une sonde spatiale une trajectoire parabolique de foyer le centre C d'un astre donné (comme la Terre) et de sommet un point donné M de l'espace tel que R=CM en communiquant à l'engin une vitesse de valeur égale à la vitesse de libération pour R, et dirigée perpendiculairement à la direction radiale (CM). Ceci peut parfaitement être fait depuis une orbite elliptique initiale, en partant du périastre : en effet, en ce point la vitesse de l'engin est orthoradiale et a une valeur maximale, bien que la vitesse de libération au périastre soit plus élevée que celle à l'apoastre.

Cas du mouvement hyperbolique

Principaux paramètres de la trajectoire

Si e > 1, la valeur de r tend vers l'infini pour $w_{\infty }=\arccos {\left(-{\frac {1}{e}}\right)}$, les deux directions $w=w_{\infty }$ et $w=2\pi -w_{\infty }$, symétriques par rapport au grand axe (FP), définissent les asymptotes à la courbe de la trajectoire. Les valeurs de w telles que $w_{\infty }<w<2\pi -w_{\infty }$ correspondent à des valeurs de r négatives : il s'agit en fait de l'autre branche de l'hyperbole, qui serait parcourue dans le cas d'un champ répulsif (cf. diffusion Rutherford). Physiquement, la seule branche parcourue est celle la plus proche du foyer F.

Principaux paramètres du mouvement képlerien hyperbolique, pour e=2 et e=5.

Les deux asymptotes de l'hyperbole se coupent en un point O du grand axe, centre de symétrie de la courbe mathématique complète, à deux branches. On peut définir un point F’ symétrique de F par rapport à O qui est le foyer de la deuxième branche de l'hyperbole. Comme pour toutes les coniques, le périastre P est situé à la distance $q={\frac {p}{1+e}}$ du foyer et constitue le sommet de la trajectoire. On peut définir un "apoastre" A correspondant à w = p, et correspondant au sommet de la deuxième branche. Ce point est situé à la distance $Q={\frac {p}{\left|1-e\right|}}={\frac {p}{e-1}}$ du foyer, la distance entre P et A correspondant étant $2a{\stackrel {\text{déf}}{=}}Q-q=2{\frac {p}{e^{2}-1}}$. Cette valeur du "demi-grand axe" a permet d'écrire $p=a(e^{2}-1)$, $q=a(e-1)$ et $Q=a(e+1)$, la distance du centre de symétrie O à un foyer étant ae.

La courbe ci-contre résume les caractéristiques de la trajectoire en montrant deux exemples d'hyperboles, de même paramètre, une avec e = 2 et l'autre avec e=5.

Aspects énergétiques - équation des forces vives

Il est possible de mettre en évidence comme pour le cas de la trajectoire elliptique une relation entre H et le demi-grand axe a de l'hyperbole défini plus haut. En effet au périastre P où r = q = a(e–-1) la vitesse est purement orthoradiale et l'énergie mécanique H s'exprime sous la forme : $H={\frac {1}{2}}\mu v^{2}-{\frac {GM_{T}\mu }{q}}={\frac {L^{2}}{2\mu a^{2}(e-1)^{2}}}-{\frac {GM_{T}\mu }{a(e-1)}}$, or on a $L^{2}=GM_{T}\mu ^{2}p=GM_{T}\mu ^{2}a(e^{2}-1)$, ce qui donne par substitution dans l'équation précédente :

$H={\frac {GM_{T}\mu }{a(e-1)}}\left({\frac {e-1)}{2}}-1\right)$, soit finalement la relation : $H={\frac {GM_{T}\mu }{2a}}$, (13).

Cette relation est identique à celle obtenue pour le mouvement elliptique, en changeant a en –a. On obtient alors, en procédant de la même façon que pour le cas elliptique, l'équation des forces vives pour le mouvement hyperbolique :

$v={\sqrt {GM_{T}\left({\frac {2}{r}}+{\frac {1}{a}}\right)}}$, (14).

Illustrations

Quelques animations représentant les orbites de deux corps (disques blancs) autour du barycentre (croix rouge).

Deux corps de même masse	Un corps de masse légèrement plus importante que l'autre	Un corps de masse plus importante que l'autre	Un corps de masse bien plus importante que l'autre
Deux corps de même masse, orbites elliptiques		Un corps de masse plus importante que l'autre en mouvement linéaire

Relativité générale

En relativité générale, il n'existe pas de solution analytique et exacte de l'équation d'Einstein pour deux corps ou plus.

Références

Ouvrages utiles :

• Dumoulin et Parisot, Astronomie pratique et informatique, Masson, Paris, 1987.
• Perez, Cours de physique : mécanique - 4^e édition, Masson, Paris, 2001.
• Landau et Lifchitz, Cours de physique - Tome I : Mécanique, Ellipses - Marketing, Paris, 1994.

Sur Internet :

• Luc Duriez, Cours de mécanique céleste classique, université Lille-I / IMCCE, 2002, très complet, disponible à l'adresse suivante en pdf : http://www.imcce.fr/fr/formations/cours/CoursMC_Duriez/mc/CoursMCecr_Duriez.pdf .

Articles connexes

• Problème à N corps
• Mouvement à force centrale
• Cinématique > Mouvement elliptique
• Lois de Kepler