Démonstration des lois de Kepler

Les lois de Kepler ont été découvertes à partir des observations de Tycho Brahe à la fin du XVI siècle et de leur analyse poussée par Johannes Kepler dans les décennies qui ont suivi.

En 1687 dans les Philosophiae Naturalis Principia Mathematica, Isaac Newton introduit la force de gravitation, qui se veut être à la fois une explication aux mouvements des planètes, et à la pesanteur sur Terre. Le problème exposé ici est de démontrer que la seule expression de la loi universelle de la gravitation, combinée au principe fondamental de la dynamique, fournit une justification des lois empiriques de Kepler. Il fut résolu par Newton.

Une autre démonstration géométrique a été donnée par Richard Feynman, durant ses cours. Il ne l'a pas publiée, mais Brian Beckman l'a fait en 2006 dans "The Journal of Symbolic Geometry", volume 1.

Pour simplifier, prenons le Soleil comme origine du référentiel, l'axe z perpendiculaire à la droite passant par le Soleil et la planète et perpendiculaire à la direction de la vitesse de la planète au temps t = 0. L'axe x dans la direction correspondante à la distance la plus petite entre le Soleil et la planète. La distance la plus grande entre le Soleil et la planète sera dans la direction -x.

• $M$ la masse du Soleil.
• $m$ la masse de la planète.
• $O$ la position du Soleil. C'est l'origine du référentiel.
• $P$ la position de la planète.
• $t=0$ l'instant où la planète se trouve le plus proche du Soleil.
• $P_0$ la position de la planète au temps $t=0$.

Elle s'appelle le périhélie. C'est le point de la trajectoire près de Hélios; soit le plus proche du Soleil.

• $\overrightarrow{r}=\overrightarrow{OP}$ le vecteur allant du Soleil à la planète.
• $r=|\overrightarrow{r}|=|\overrightarrow{OP}|$ = la distance du Soleil à la planète.
• $\overrightarrow{r_0}=\overrightarrow{O{P}_0}$ le vecteur allant du Soleil à la planète au temps $t=0$.
• $\overrightarrow{V}$ la vitesse de la planète.
• $\overrightarrow{V_0}$ la vitesse de la planète au temps $t=0$.
• $\overrightarrow{F}$ la force d'attraction de la planète par le Soleil.
• $\overrightarrow{L}=m\cdot \overrightarrow{r}\wedge \overrightarrow{V}$ le moment cinétique relativement à l'origine où se trouve le Soleil.

C.f. le produit vectoriel "$\wedge$" pour des informations sur cet opérateur.

Dans un référentiel immobile par rapport au Soleil, la trajectoire d'une planète se trouve dans un plan.

• Démonstration :

Cela résulte du fait que le Soleil attire la planète selon une force centrale. C'est-à-dire une force qui est toujours dirigée de la planète vers le Soleil.

En effet, étant données une position $\overrightarrow{r_0}$ et une vitesse $\overrightarrow{V_0}$ initiales, cela définit un plan. Selon nos conventions ci-dessus, c'est le plan passant par l'origine O, contenant les axes x et y.

Puisque la force est centrale, elle et l'accélération sont dans une direction se trouvant dans ce même plan. Donc les variations de vitesses et les variations de positions resteront dans ce même plan. En conclusion toute la trajectoire restera dans ce plan.

La deuxième loi de Kepler donnera une deuxième démonstration de cette partie.

Soit A(t) l'aire de la surface balayée par le rayon vecteur $\vec{r}$ durant le mouvement, alors cette seconde loi stipule que des aires égales sont balayées dans des temps égaux.

• Démonstration :

En dérivant le moment cinétique $\overrightarrow{L}=m\cdot \overrightarrow{r}\wedge \overrightarrow{V}$ par rapport au temps, on obtient :

$\frac{d\overrightarrow{L}}{dt}=m\cdot \overrightarrow{V}\wedge \overrightarrow{V}+m\cdot \overrightarrow{r}\wedge \overrightarrow{F}/m$ qui vaut 0.

Donc le moment cinétique est un vecteur constant. Cela résulte du fait que la force est centrale (cf. mouvement à force centrale)

D'autre part : $\frac{|\vec{L}|}{m}=r\cdot v\cdot sin(\theta )=2\cdot \frac{dA(t)}{dt}$

En conséquence : $A(t)=\frac{|\vec{L}|}{2\cdot m}\cdot t$

Dans un référentiel immobile par rapport au Soleil, la trajectoire d'une planète est elliptique, un foyer étant le Soleil. Le Soleil n'est un des foyers qu'approximativement, du fait que sa masse M est très supérieure à celle de la masse m de la planète. Pour être exact, il faudrait se placer au centre de gravité du système Soleil - planète.

• Démonstration :

Ce fait est directement lié au fait que l'hodographe est un cercle dans le cas de l'attraction universelle de Newton.

Le principe fondamental de la dynamique s'écrit :

$\frac{d\vec{V}}{dt}=-\frac{G\cdot M}{r^2}\cdot \hat{r}$     où     $\hat{r}=\frac{\vec{r}}{r}$

Par changement de variable, on a :

(Attention : L'angle noté θ dans ce paragraphe est différent de celui du paragraphe précédent)

Rappelons l'expression de la loi des aires obtenue au paragraphe précédent : $A(t)=\frac{|\vec{L}|}{2\cdot m}\cdot t$

ainsi, l'aire balayée pendant un temps $dt$ égale $\frac{|\vec{L}|}{2\cdot m}\cdot dt$. Or cet aire vaut $\frac{1}{2}\cdot r^2\cdot d\theta$.

Donc $\frac{d\theta }{dt}=\frac{|\vec{L}|}{m\cdot r^2}$

En combinant ces deux relations, on obtient :$\frac{d\vec{V}}{d\theta }=-\frac{G\cdot M\cdot m}{|\vec{L}|}\cdot \hat{r}$

D'où en intégrant par rapport à θ : $\vec{V}=\frac{G\cdot M\cdot m}{|\vec{L}|}\cdot \hat{\theta }+{\vec{V}}_0$    car    $\frac{d\hat{\theta }}{d\theta }=-\hat{r}$

$\hat{\theta }$ = le vecteur unité perpendiculaire à $\hat{r}$ dans le plan de la trajectoire, dirigé dans le sens le plus proche de celui de la vitesse.

Ce résultat s'appelle : théorème de Hermann, Laplace, Runge, Lenz, Hamilton : cf. vecteur de Runge-Lenz.

L'hodographe est un cercle excentré par rapport à l'origine des vitesses. Il en résulte que la trajectoire est une conique : si l'origine des vitesses est à l'intérieur du cercle, la conique est une ellipse ; si elle est à l'extérieur, c'est une hyperbole ; cas limite : une parabole. Ce théorème de cinématique est très vieux, mais on l'attribue à tort à Hamilton ; il était encore enseigné dans le cours de cosmographie de « math-élem » (terminale S actuelle) : cf. par exemple Lebossé, cours de mathématiques élémentaires. On va en donner la démonstration due à Landau, très algébrique :

Le théorème précédent se réécrit après multiplication vectorielle par $\frac{\vec{L}}{G\cdot M\cdot m}$ et simplifications :

$\vec{e}$ s'appelle le vecteur d'exentricité. On montre qu'il est constant. (C.f. lien externe n°2)

puis en effectuant le produit scalaire avec $\vec{r}$ et en simplifiant :

$e\cdot r\cdot cos(\theta )=p-r$, avec $p=\frac{L^2}{G\cdot M\cdot m^2}$

soit $r=\frac{p}{1+e\cdot cos(\theta )}$

Ce qui est la définition focale d'une conique en coordonnées polaires, d'excentricité $e=\frac{c}{a}$ et de paramètre $p=\frac{b^2}{a}$, d'angle polaire ayant pour origine le périhélie, comme il est usuel de le faire.

Si on prend l'aphélie comme origine : $r=\frac{p}{1-e\cdot cos(\theta )}$

Le carré de la période $T$ varie comme le cube du demi-grand axe $a$ : $\frac{a^3}{T^2}=\frac{G\cdot M}{4\cdot {\pi }^2}$

$G\cdot M$ s'appelle la constante de Gauss : elle est connue avec une extraordinaire précision, dix chiffres significatifs et vaut 0.01720209895 (alors que $G$ n'est connue qu'avec 5 chiffres significatifs).

Démonstration :

• En utilisant : $A(T)=\pi \cdot a\cdot b=\frac{L}{2\cdot m}\cdot T$   et   $p=\frac{b^2}{a}=\frac{L^2}{G\cdot M\cdot m^2}$,

puis en éliminant $L$ entre ces égalités, on obtient le résultat annoncé.

Par conséquent, toutes les ellipses de même grand axe, quelle que soit leur excentricité $e$, ont la même période de révolution jusqu'à la circulaire où $e=0$.