Voir la source de Problème à deux corps

[[Image:Eclipsing binary star animation 2.gif|thumb|é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 [[point matériel|points matériels]]<ref>Ceci est valable dans le cas où les dimensions des deux corps sont très petites devant leur distance au cours du mouvement.</ref> en interaction mutuelle ([[Force_conservative|conservative]]), le système global étant considéré comme isolé<ref>Cette approximation revient à négliger l'influence d'autres corps, à considération de l'importance relative de leurs actions sur chacun des deux corps. Par exemple pour le mouvement d'une planète autour du soleil, l'interaction dominante est bien sûr celle de l'étoile sur la planète, on peut négliger au moins en première approximation les effets des interactions des autres corps du système solaire tant sur le Soleil que sur la planète considérée. Toutefois il conviendra d'en tenir compte de façon perturbative pour une description plus complète.</ref>. Dans cet article, seul sera abordé le problème à deux corps en [[mécanique newtonienne|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 [[gravitation|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 à N corps|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é<ref>En fait, à deux problèmes à un corps indépendants, mais le mouvement du centre d'inertie est trivial.</ref>. 

Par ailleurs, les résultats obtenus permettent de rendre compte des trajectoires des [[planète|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 {{s-|XVII|e}}. 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 [[Isaac Newton|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 [[orbite|éléments d'orbite]] ainsi que les équations de [[Équation de Kepler|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 [[Orbite|éléments d'orbite]]).

== Situation envisagée et notations ==
Le problème à deux corps est celui de deux corps de masse ''m''{{ind|1}} et ''m''{{ind|2}}, assimilés à des [[point matériel|points matériels]] ''M''{{ind|1}} et ''M''{{ind|2}}, respectivement, en interaction mutuelle. La force exercée par ''M''{{ind|1}} sur ''M''{{ind|2}} dérive d'un potentiel ''attractif'' ''V''(''r'') et est notée <math>\overrightarrow{F_{12}}</math>: du fait de la troisième loi de Newton (ou principe des actions réciproques) il est évident que <math>\overrightarrow{F_{21}} = - \overrightarrow{F_{12}}</math>.

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

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

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

:<math>\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}</math>, (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 :

:<math>m_1\ddot{\overrightarrow{r_1}} + m_2\ddot{\overrightarrow{r_2}} = \left( m_{1} + m_{2}\right) \ddot{\overrightarrow{R_{C}}} = \overrightarrow{0}</math>, avec ''C'' centre de masse du système, de vecteur position <math>\overrightarrow{R_{C}} = \overrightarrow{OC} = \frac{m_1\overrightarrow{r_1} + m_2\overrightarrow{r_2}}{m_{1} + m_{2}}</math>.

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<sub>c</sub>'') (qui sera galiléen, du fait du mouvement rectiligne et uniforme du corps auquel il est lié, (''R'') étant supposé galiléen), dit [[barycentre_(physique)|''barycentrique'']], pour réécrire les équations du mouvement précédentes.

=== Introduction de la notion de particule fictive ===

En posant <math>\overrightarrow{r} \stackrel{\text{déf}}{=} \overrightarrow{r_2} - \overrightarrow{r_1}</math> il est possible d'écrire :
:<math>\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}</math>, (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 <math>\ddot{\vec{R}_c}=\vec{0}</math>, pour obtenir :
:<math>\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}}</math>.
Cette équation est en fait celle du mouvement d'un seul corps, à trois degrés de liberté :
:<math>\mu \ddot{\overrightarrow{r}} = \overrightarrow{F}</math>, 
avec <math>\mu \stackrel{\text{déf}}{=} \frac{m_{1}m_{2}}{m_1 + m_2}</math>, ''masse réduite'' du système, et <math>\overrightarrow{F} =\overrightarrow{F_{12}}</math>.

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 <math>\overrightarrow{r}</math>, les trajectoires des corps ''M''{{ind|1}} et ''M''{{ind|2}} se déduisant par [[homothétie]], d'après les formules précédentes sur <math>\vec{r}_{1,2}</math> .

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 <math>m_1 \gg m_2</math>, 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, <math>\mu \simeq m_2</math>. À 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 ===
{{Article détaillé|Moment cinétique spécifique}}

Dans le cas particulier très important d'une force centrale, on a <math>\overrightarrow{F} = F(r)\overrightarrow{e_r}</math> avec <math>\overrightarrow{e_r} \stackrel{\text{déf}}{=} \dfrac{\overrightarrow{r}}{r}</math>, le théorème du [[moment cinétique]] au centre de force noté ''O'' s'écrit :
:<math>\dot{\overrightarrow{L}} = \overrightarrow{r}\times \overrightarrow{F} = \overrightarrow{0}</math>,
ce qui implique <math>\overrightarrow{L} = \overrightarrow{r}\times\left(\mu\dot{\overrightarrow{r}}\right) = \overrightarrow{\mathrm{cte}}</math>.

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''<ref>Par homothétie, il en est bien sûr de même de celles des particules réelles ''M1'' et ''M2''.</ref>{{,}}<ref>Toutefois, si <math>\overrightarrow{L} =\overrightarrow{0}</math>, le mouvement est dit ''dégénéré'' et se réduit à une droite, la notion de planéité de la trajectoire n'a aucun sens</ref>, le problème est donc à deux degrés de liberté. 

Dans le plan de la trajectoire, défini comme celui généré par <math>\overrightarrow{r}\left(t = 0\right) = \overrightarrow{r}_0</math> et <math>\dot{\overrightarrow{r}}\left(t = 0\right) = \overrightarrow{v}_0</math>, il est judicieux de se placer en coordonnées cylindro-polaires d'axe la direction ''Oz'' de <math>\overrightarrow{L} = L\overrightarrow{e_z}</math>, avec ''θ'' angle entre <math>\overrightarrow{r}_0</math> et <math>\overrightarrow{r}</math>, il vient :
:<math>\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}</math>,
par suite:
:<math>\mu r^2\dot{\theta} = L = \mathrm{cte}</math>(3).

Or l'aire élémentaire balayée par le rayon-vecteur <math>\overrightarrow{r}</math> pendant d''t'' est donnée par :
: <math>\mathrm{d}\mathcal{A} = \dfrac{1}{2}r^2\mathrm{d}\theta</math>.
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) :

:<math>\dot{\mathcal{A}} = \frac{\mathrm{d}\mathcal{A}}{\mathrm{d}t} = \frac{1}{2}C = \mathrm{cte}</math>,
avec <math>C \stackrel{\text{déf}}{=} \frac{L}{\mu}</math>, 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 <math>\dot{\theta}</math> 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 ===
{{article détaillé|Énergie orbitale spécifique}}

Le mouvement étant conservatif puisque la force dérive d'une énergie potentielle ''V''(''r''), l'énergie totale est une [[intégrale première|intégrale première du mouvement]] :
:<math>H = \frac{1}{2}\mu\dot{\overrightarrow{r}}^2 +V(r) = \mathrm{cte}</math>, or <math>\dot{\overrightarrow{r}}^2 = \dot{r}^2 + r^2\dot{\theta}^2</math>, soit en utilisant l'expression de la valeur du moment cinétique ''L'' :
:<math>H = \frac{1}{2}\mu\dot{r}^2+U_{\text{eff}}\left(r\right)</math>,                 (4),
avec <math>U_{\text{eff}}\left(r\right) \stackrel{\text{déf}}{=} \frac{L^2}{2\mu r^2} +V(r)</math> ''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 <math>\dot{r}</math>, il vient : <math>\dot{r}=\frac{\mathrm{d}r}{\mathrm{d}t}= \sqrt{\frac{2}{\mu}\left(H-V(r)\right)-\frac{L^2}{\mu^2r^2}}</math>.

Il est alors possible de séparer les variables, et d'intégrer entre deux instants ''t''{{ind|0}} et ''t'', auxquels correspondent respectivement les positions radiales ''r''{{ind|0}} et ''r'', pour obtenir :
:<math>t-t_0 = 	\int_{r_0}^{r} \frac{\mathrm{d}r'}{\sqrt{\frac{2}{\mu}\left(H-V(r')\right)-\frac{L^2}{\mu^2r'^2}}}</math>,                         (4bis).

Ceci correspond, de façon implicite, à l'équation horaire <math> r(t) </math>.

Compte tenu de (4), il est alors possible d'obtenir une expression similaire pour ''θ'' :
:<math>\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}}}</math>,               (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 <math>U_{\text{eff}}\left(r\right)</math>, le terme <math>\frac{L^2}{2\mu r^2}</math>, 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 <math>\lim_{r\to 0^+}r^2\left|V(r)\right| = 0</math>, et donc le terme en {{frac|1|''r''{{2}}}} domine quand <math>r\to 0^+</math>. 

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, <math>\lim_{r \to \infty} V(r) = 0^+</math>. Par convention ''V''(''r'') < 0<ref>Cette dernière condition n'est pas absolument nécessaire, on obtient également une cuvette de potentiel, certes ''infinie'', avec un potentiel spatial harmonique de la forme <math>V(r)=\dfrac{1}{2}kr^2</math> avec ''k''>0, mais cet exemple ne sera pas considéré ensuite.</ref>.

Avec ces conditions, valables pour les potentiels physiques usuels, le potentiel effectif présente un minimum absolu unique, noté <math>-U_{\min}<0</math> pour <math>r = r_m</math> tel que <math>\left.\frac{\mathrm{d}U_{\text{eff}}}{\mathrm{d}r}\right|_{r=r_m} = 0</math>, et possède donc une ''cuvette de potentiel'' (cf. figure ci-contre, avec un potentiel newtonien).
[[Fichier:Courbe_potentiel_effectif.jpg|thumb|alt=alternative à l'image|Allure de la courbe du potentiel effectif ''U''{{ind|eff}}(''r'') avec les contributions des différents termes]]

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

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

* Énergie totale ''H''>0 : puisque <math>U_{\text{eff}}\left(r\right) \longrightarrow 0^-</math> quand <math>r \to \infty</math>, cette condition est vérifiée pour toute valeur de <math> r \geqslant r_{\min}</math> avec ''r''{{ind|min}} distance minimale d'approche<ref>Cette distance minimale d'approche correspond à une vitesse radiale nulle à distance finie.</ref> telle que <math>U_{\text{eff}}\left(r_{\min}\right) = H</math>. La particule fictive peut donc aller à l'infini avec une [[vitesse radiale]] positive, égale à <math>\dot{r}_{\max} = \sqrt{\dfrac{2H}{\mu}}</math><ref>En l'infini, la vitesse est purement radiale : en effet, le terme orthoradial est celui en {{frac|1|''r''{{2}}}}, et donc tend vers 0 à grande distance.</ref>. 

* É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 <math>U_{\text{eff}}\left(r_{\min}\right) = 0</math>, soit <math>V\left(r_{\min}\right)= -\frac{L^2}{2\mu r_{\min} ^2} </math>.

* Énergie totale <math>0 > H > -U_{\min}</math> : 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 <math>U_{\text{eff}}\left(r_{\min,\max}\right) = H <0</math> (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''{{ind|min}} à ''r''{{ind|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''{{ind|max}} à ''r''{{ind|min}}, le mouvement radial est donc périodique de période ''T'' donnée par :
:<math>T = 2\int_{r_{\min}}^{r_{\max}} \frac{\mathrm{d}r}{\sqrt{\frac{2}{\mu}\left(H-V(r)\right)-\frac{L^2}{\mu^2r^2}}}</math> 
:, et en une période ''radiale'' où l'angle ''θ'' varie de la quantité Δ''θ'', donnée par l'expression intégrale de ''θ'' :

:<math>\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}}}</math>.

[[Fichier:Etude_potentiel_effectif.png|thumb|left|alt=alternative à l'image|Mouvements permis dans le potentiel effectif ''U''{{ind|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 <math>\Delta\theta = \frac{2\pi m}{n}</math> avec ''m'' et ''n'' entiers. Dans ce cas, et seulement dans ce cas, le rayon vecteur <math>\overrightarrow{r}</math> retrouve sa valeur initiale après ''n'' périodes « radiales » ''T'', puisque qu'alors il aura « tourné » de 2''mp'' : ''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''{{2}} (ce dernier cas ne sera pas considéré plus avant) : ce résultat constitue le [[théorème de Bertrand]].

* Énergie totale <math>H = - U_{\min} < 0</math> : il s'agit du cas limite où la seule valeur permise de ''r'' est ''r{{ind|m}}'', et dans ce cas la trajectoire est un ''cercle'' de ce rayon.

* Énergie totale <math>H < - U_{\min}</math> : 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 <math>L \ne 0</math>. Si ''L'' = 0, on a tout simplement <math>\dot{\theta}=0</math> à 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 à <math>V(r) \leqslant H</math>, vérifiée dans tous les cas si <math>H \geqslant 0</math>. 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 ==
{{Article détaillé|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 <math>V(r)= -\frac{Gm_1 m_2}{r}</math>, et donc <math> \overrightarrow{F} = - G \frac{m_{1}m_{2}}{r^2}\overrightarrow{e_r} =- G \frac{\mu M_T}{r^2}\overrightarrow{e_r}</math>, avec <math>M_T \stackrel{\text{déf}}{=} m_1 + m_2</math> 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 <math>1/r</math> se caractérise par l'existence d'un invariant supplémentaire particulier, l'[[vecteur de Runge-Lenz|invariant de Runge-Lenz]], donné par:
: <math>\overrightarrow{e} = \frac{\left(\dot{\overrightarrow{r}}\times\overrightarrow{L}\right)}{GM_T\mu}-\overrightarrow{e_r}</math>,   (5)<ref>En toute rigueur le vecteur de Runge-Lenz est classiquement défini par <math>\overrightarrow{p}\times\overrightarrow{L}-GM_T\mu^2\overrightarrow{e_r}</math>.</ref>
{{démonstration|
:<math>\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]</math>, où il a été tenu compte de <math>\overrightarrow{L}=\overrightarrow{\mathrm{cte}}</math>, ce qui donne :

:<math>\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}}</math>, où il a été tenu compte de la planéité de la trajectoire, or compte tenu de <math>\dot{\theta}\overrightarrow{e_{\theta}} = \dot{\overrightarrow{e_r}}</math> il vient aussitôt :
:<math>\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}</math>, d'où le résultat.}}

==== Équation de la trajectoire de la particule fictive ====
Il est évident que <math>\overrightarrow{L}\cdot\overrightarrow{e} = 0</math> donc <math>\overrightarrow{e}</math> est contenu dans le plan du mouvement. Par suite, il est possible de prendre comme angle polaire l'angle ''w'' entre <math>\overrightarrow{e}</math> et <math>\overrightarrow{r}</math>, avec bien sûr <math>\dot{\theta}=\dot{w}</math><ref>Cela revient à un changement d'origine, qui ne remet pas en cause les résultats précédents, notamment le fait que ''L'' = cte, puisque ceci ne dépend que de la valeur de la vitesse angulaire <math>\dot{\theta}=\dot{w}</math>.</ref>, et en notant ''e'' la norme de <math>\overrightarrow{e}</math> il est facile de vérifier :
:<math>\overrightarrow{r}\cdot\overrightarrow{e} = re \cos{w} =  \frac{\overrightarrow{r}\cdot\left(\dot{\overrightarrow{r}}\times\overrightarrow{L}\right)}{GM_T\mu} - r</math>, soit compte tenu de l'identité 
:<math>\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}</math> :

:<math>r\left(1+e \cos{w}\right)= p</math>, avec <math>p \stackrel{\text{déf}}{=} \frac{L^2}{GM_T\mu^2}</math>. 

Physiquement, ''p'' = ''r{{ind|m}}'', valeur de ''r'' pour laquelle ''U''{{ind|eff}}(''r'') est minimal. De fait puisque <math>U_{\text{eff}}(r) = \frac{L^2}{2\mu r^2}-\frac{GM_T \mu}{r}</math>, il vient facilement par dérivation à <math>\frac{L^2}{\mu r_m}=GM_T \mu</math> soit en effet ''r{{ind|m}}'' = ''p''.

L'équation de la trajectoire obtenu est donc celle d'une '''conique''' d'excentricité ''e'' et de paramètre ''p'' :
: <math>r = \frac{p}{\left(1+e \cos{w}\right)}</math>,       (6), 
dont le centre de force occupe l'un des foyers. Le vecteur <math>\overrightarrow{e}</math> est donc dirigé vers le point de distance minimale (ou ''périastre''), distance notée ''q'' avec <math>q=p/\left(1+e\right)</math>, correspondant à ''w'' = 0. En ce point, la vitesse angulaire orbitale <math>\dot{w_0} = \dot{\theta_0}</math> 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 [[Lois de Kepler|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 [[Loi de conservation|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''{{exp|2}}, 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. [[Vecteur de Runge-Lenz|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 [[Niveau d'énergie|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''{{exp|2}} 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{{ind|c}}''), l'axe ''Ox'' de la conique, orienté vers le point de distance minimale au foyer de vecteur unitaire <math>\overrightarrow{e_x}</math>, ''Oz'' la direction du moment cinétique <math>\overrightarrow{L}</math>, et comme au périastre <math>\dot{r} =0</math> et ''r'' = ''q'', l'intégrale première de l'énergie peut être exprimée en fonction de cette distance minimale, il vient :
:<math>H =  \frac{L^2}{2\mu q^2}-\frac{GM_T\mu}{q}</math>,
ce qui implique <math>L^2=2\mu q^2H+2GM_T\mu^2q</math>.

De même, il est possible d'exprimer le vecteur excentricité de façon simple :
: <math>\overrightarrow{e}=\left(\frac{L^2}{GM_T\mu^2 q}-1\right)\overrightarrow{e_x}</math>.
{{démonstration|<math>\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}</math>, or puisque <math>q\dot{\theta_0} = \frac{L}{\mu q}</math>, d'où le résultat.}}

En remplaçant dans cette expression celle obtenue pour ''L''{{2}}, l'expression de l'excentricité de la conique se met sous la forme :
:<math>e=1+\frac{2qH}{GM_T\mu}</math>.

Il est alors possible d'éliminer la distance minimale ''q'' dans cette expression compte tenu de <math>q = \frac{p}{1+e} = \frac{L^2}{G^2M_T ^2 \mu^3 (1+e)}</math>, 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é :
:<math>(e+1)(e-1)=\frac{2L^2H}{G^2M_T ^2\mu^3}</math>, soit finalement : <math>e=\sqrt{1+\frac{2L^2 H }{G^2M_T ^2\mu^3}}</math>.

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

* Cas de la trajectoire hyperbolique : ''e'' > 1 ce qui implique ''H'' > 0<ref>On pourrait aussi avoir ''L''=0 pour ''H'' quelconque, toutefois le paramètre ''p'' est alors nul et on n'obtient plus une parabole, mais simplement une "conique" dégénérée en droite : voir remarque plus haut sur la dégénérescence du mouvement. Ce cas trivial ne sera plus envisagé par la suite.</ref>, comme il a été dit le mobile fictif peut aller en l'infini, avec une vitesse purement radiale non nulle, donnée par <math>\dot{r}_{\max} = \sqrt{\dfrac{2H}{\mu}}</math> ;

* 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 <math>-\frac{GM_T\mu}{2q} < H <0</math> .

* Cas de la trajectoire circulaire : ''e'' = 0 donc <math>H= -\frac{GM_T\mu}{2q} = - \frac{G^2M_T ^2 \mu^3}{2L^2}</math>.

'''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 "[[orbite|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{{ind|1}}'' et ''F{{ind|2}}'' est constante : <math>d(F_1,M)+d(F_2,M) = 2a</math>, ''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). 
[[Fichier:Paramètres d'une ellipse.png|thumb|alt=alternative à l'image|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 à <math>w=\pi</math>, il vient : <math>Q = p/\left(1-e\right)</math>, et on en déduit <math>a = p/\left(1-e^2\right)</math>. 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 <math>c=a-q = e\frac{p}{1-e^2}=ea</math>. On déduit le ''demi-petit axe'', noté ''b'' : <math>b=\sqrt{a^2-c^2}=a\sqrt{1-e^2}</math>. On peut éliminer ''e'' entre les équations donnant ''a'' et ''b'', il vient l'expression du paramètre de l'ellipse : <math>p=b^2/a</math>.

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

==== Aspects énergétiques - Équation des forces vives ====
Pour obtenir l'équation qui tire son nom, « [[Force vive (physique)|forces vives]] », du concept [[Leibniz|leibnizien]] de la force, on remplace l'expression de la distance au périastre <math>q = a\left(1-e\right)</math> dans la relation (9) excentricité-énergie, il vient :
:<math>1-e = -\frac{2qH}{GM_T\mu} = -\frac{2Ha\left(1-e\right)}{GM_T\mu}</math>, d'où l'on tire l'expression de l'énergie mécanique totale en fonction de la valeur du demi-grand axe ''a'' :
:<math>H =- \frac{GM_T\mu}{2a}</math>,     (10bis). 

Par suite, l'énergie totale est égale à la moitié de l'énergie potentielle de gravitation pour ''r'' = ''a''. Comme ''H'' = cte, il vient :
:<math>H = - \frac{GM_T\mu}{2a} = \frac{1}{2}\mu v^2-\frac{GM_T \mu}{r}</math>, 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'' :
:<math>v=\sqrt{GM_T \left(\frac{2}{r}-\frac{1}{a}\right)}</math>,               (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 <math>\left|H\right|\to 0</math> on a <math>a\to\infty</math>, 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 <math>\dot{\mathcal{A}} = \frac{L}{2\mu}</math>, on obtient en intégrant sur la période ''T'' du mouvement l'aire totale de l'ellipse, égale à <math>\pi ab</math>, d'où l'identité <math>\pi ab = \frac{LT}{2\mu}</math>, qui donne <math>L^2=\frac{4\pi^2\mu^2a^2b^2}{T^2} = 4\pi^2\mu^2p\frac{a^3}{T^2}</math>, où il a été fait usage de <math>p=b^2/a</math>. 

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

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 <math>\frac{a^3}{T^2} = \frac{k}{4\pi^2}</math>, avec <math>k\stackrel{\text{déf}}{=} GM_S</math> ''constante de Gauss''.

En faisant cette approximation, on peut écrire pour deux planètes quelconques du système solaire, avec des notations évidentes :
:<math>\frac{a_1 ^3}{T_1 ^2}=\frac{a_2 ^3}{T_2 ^2}</math>.
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étrie|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 : <math>n \stackrel{\text{déf}}{=} \frac{2\pi}{T}</math>, on a donc d'après la troisième loi de Kepler :
:<math>n = \sqrt{\frac{GM_T}{a^3}}</math>.

==== 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 <math>\dot{\theta} = L/\mu p^2</math> 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 [[#Relation entre l'excentricité et les intégrales premières du mouvement|3.1-4]], cela correspond au minimum physique de l'énergie totale ''H'' avec  <math>H= -\frac{GM_T\mu}{2R} = - \frac{G^2M_T ^2 \mu^3}{2L^2}</math>, donc une énergie cinétique radiale <math>\dfrac{1}{2}\mu\dot{r}^2</math> nulle. L'énergie potentielle effective <math>U_{\text{eff}}(r)</math> est donc minimale, ce qui correspond physiquement au point où les termes de barrière centrifuge en 1/''r''{{2}} 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 :

:<math>v_1 = \sqrt{\frac{GM_T}{R}}</math>,              (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'' = {{unité|6400|km}} (surface de la Terre), ce qui donne environ ''v''{{ind|1}} ≈ {{unité|7,9|km||s|-1}}.

=== 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''{{ind|1}}, l'autre foyer ''F''{{ind|2}} 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.

[[File:Mouvement keplerien parabolique.png|thumb|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 : <math>r = \frac{p}{1+ \cos w}</math>. Le périastre ''P'' correspond à ''w'' = 0 et est situé à la distance ''q'' = ''p''/2 du foyer ''F'', et on a <math>r \to +\infty</math> pour <math>w \to \pi</math>. 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é : <math>\frac{1}{2}\mu v^2=\frac{GM_T \mu}{r}</math>, 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 :
<math>v=\sqrt{\frac{2GM_T}{r}}</math>,   (12).

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

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 <math>v_2=\sqrt{\frac{2GM_T}{R}}=\sqrt{2}v_1</math>.

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''{{ind|2}} ≈ {{unité|11,2|km||s|-1}}.

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)<ref>On confond ici centre de l'astre et centre de masse du système {sonde - astre}, vu les rapports entre les masses des deux objets.</ref> 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 <math>w_\infty = \arccos{\left(-\frac{1}{e}\right)}</math>, les deux directions <math>w = w_\infty</math> et <math>w=2\pi-w_\infty</math>, symétriques par rapport au grand axe (''FP''), définissent les asymptotes à la courbe de la trajectoire. Les valeurs de ''w'' telles que <math>w_\infty < w<2\pi-w_\infty</math> 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''.
[[File:Trajectoire hyperbolique.png|thumb|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 <math>q=\frac{p}{1+e}</math> 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 <math>Q=\frac{p}{\left|1-e\right|}=\frac{p}{e-1}</math> du foyer, la distance entre ''P'' et ''A'' correspondant étant <math>2a\stackrel{\text{déf}}{=} Q-q=2\frac{p}{e^2-1}</math>. Cette valeur du "demi-grand axe" ''a'' permet d'écrire <math>p=a(e^2-1)</math>, <math>q=a(e-1)</math> et <math>Q=a(e+1)</math>, 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 :
<math>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)}</math>, or on a <math>L^2=GM_T \mu^2 p=GM_T \mu^2 a(e^2-1)</math>, ce qui donne par substitution dans l'équation précédente :
:<math>H=\frac{GM_T \mu}{a(e-1)}\left(\frac{e-1)}{2}-1\right)</math>, soit finalement la relation : <math>H=\frac{GM_T\mu}{2a}</math>,    (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 :
:<math>v=\sqrt{GM_T \left(\frac{2}{r}+\frac{1}{a}\right)}</math>,   (14).

== Illustrations ==
Quelques animations représentant les orbites de deux corps (disques blancs) autour du barycentre (croix rouge).
{|
| valign="top" | [[Image:Orbit1.gif|framed|Deux corps de même masse]]
| valign="top" | [[Image:Orbit2.gif|framed|Un corps de masse légèrement plus importante que l'autre]]
| valign="top" | [[Image:Orbit3.gif|framed|Un corps de masse plus importante que l'autre]]
| valign="top" | [[Image:Orbit4.gif|framed|Un corps de masse bien plus importante que l'autre]]
|-
| valign="top" colspan="2" | [[Image:Orbit5.gif|framed|Deux corps de même masse, orbites elliptiques]]
| valign="top" colspan="4" | [[Image:Wobbling binary.gif|framed|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<ref>{{article | langue=en | prénom=Ethan | nom=Siegel | titre=This is why scientists will never exactly solve general relativity | périodique=[[Forbes (magazine)|Forbes]] | date=4/12/2019 | lire en ligne= https://www.forbes.com/sites/startswithabang/2019/12/04/this-is-why-scientists-will-never-exactly-solve-general-relativity/ | accès url=libre}}</ref>{{,}}<ref>{{ouvrage | langue=en | prénom1=Michal | nom1=Křížek | prénom2=Lawrence | nom2=Somer | titre=Mathematical aspects of paradoxes in cosmology | sous-titre=can mathematics explain the contemporary cosmological crisis ? | lieu=Cham | éditeur=[[Springer Science+Business Media|Springer]], hors {{coll.}} | date=7/2023 | numéro d'édition=1 | pages totales={{XVI}}-264 | format livre={{dunité|15,4|23,6|cm}} | isbn10=3-031-31767-X | isbn1=978-3-031-31767-5 | isbn2=978-3-031-31770-5 | ean=9783031317675 | oclc=1374095021 | doi=10.1007/978-3-031-31768-2 | présentation en ligne=https://link.springer.com/book/10.1007/978-3-031-31768-2 | lire en ligne={{Google Livres|id=tPbJEAAAQBAJ}} | consulté le=22 octobre 2023 | partie={{chap.|5}}, {{sec.|5.3}} | passage=102}}.</ref>.

== Références ==
{{Références|colonnes=2}}
'''Ouvrages utiles :'''

* Dumoulin et Parisot, ''Astronomie pratique et informatique'', Masson, Paris, 1987.
* Perez, ''Cours de physique : mécanique'' - {{4e|é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|Cinématique > Mouvement elliptique]]
* [[Lois de Kepler]]
{{Palette|Orbites}}
{{Portail|physique|astronomie}}

[[Catégorie:Physique théorique]]
[[Catégorie:Physique mathématique]]
[[Catégorie:Mécanique céleste]]
[[Catégorie:Orbite]]