Lois de Kepler, démonstration - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs est disponible ici.

Les lois de Kepler ont été découvertes à partir des observations de Tycho Brahé (1546 - 1601) et de leur analyse poussée par Johannes Kepler (1571 - 1630). Elles ont été démontrées par Isaac Newton (1642 - 1727) en 1687 dans les Principia, œuvre considérée comme un monument de la pensée humaine. Les démonstrations qui suivent ne sont pas celles de Newton, pour des raisons diverses (cf. Principia et Calculus).

Des explications plus abordables se trouvent sous lois de Kepler.
Le lien externe : Démonstration des trois lois de Kepler et propriétés d'une ellipse est plus complet.

Conventions et notations pour les démonstrations qui vont suivre

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{OP_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.

Première partie de la première loi (1609)

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.

Deuxième loi, loi des aires (1609)

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.

A(t) = \frac{\| \vec L \|}{2 \cdot m} \cdot t
  • 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. force centrale, mouvement)
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

Deuxième partie de la première loi (1609)

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 :
\frac{d \vec V}{d\theta} \cdot \frac{d \theta}{dt} = - \frac{G \cdot M}{r^2} \cdot \hat r
(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. invariant 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 vitesse est à l'intérieur du cercle, la conique est une ellipse; si elle est à l'extérieur, c'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 :

\frac{\vec L \wedge \vec V}{G \cdot M \cdot m} =  \frac{1}{\| \vec L \|} \cdot \vec L \wedge \hat \theta + \frac{\vec L \wedge \vec V_0}{G \cdot M \cdot m}  = \hat r + \vec e
\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) = r - p, 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)}

Troisième loi (1618)

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'a que 5 misérables 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 supprimant L \, de ces égalités, on obtient le résultat annoncé.

On remarquera ce fait extraordinaire, toutes les ellipses de même grand axe, quelle que soit leur excentricité e \,, ont la même période jusqu'à la circulaire où e = 0 \,.

La disposition des planètes est assez remarquable (loi de Titius-Bode), tant par leur faible excentricité que par leur faible inclinaison sur le plan écliptique moyen. Le scénario cosmogonique de Levinson-Morbidelli le justifie, écartant ainsi Pluton et (90377) Sedna des planètes ordinaires et les faisant entrer dans la catégorie des " petits " objets de la ceinture de Kuiper.

Page générée en 2.808 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales | Partenaire: HD-Numérique
Version anglaise | Version allemande | Version espagnole | Version portugaise