Démonstration des lois de Kepler

à l'aide du calcul différentiel et intégral

, . Serge Robert
Cegep de Samt-Jean-sur-Richelieu

Introduction Kepler découvrit quel' orbite de Mars était une ellipse

Johannes Kepler est un astronome allemand qui vécut et publia ses résultats en 1609 dans Astronomia Nova.
de 1571 à 1630. Né de parents pauvres, il obtint une Après plusieurs ennuis dans sa vie privée, où il perdit
bourse pour étudier le culte luthérien, mais il en fut ex- femme et enfants, il retourna en Allemagne et publia
clus un peu plus tard à cause de ses positions religieu- en 1619 Harmonices Mundi dans lequel il énonce sa
ses. Son refus d' adhérer au christianisme s'ajouta à ses troisième loi. 2 Voici les fameuses lois de Kepler.
problèmes et après un édit contre les protestants, il se
réfugia à Prague auprès de l' astronome danois Tycho Première loi: Les planètes décrivent des orbites ellip-
tiques dont/ 'un des foyers est le Soleil.
Brahe. Celui-ci le chargea de calculer le mouvement
de la planète Mars en utilisant les innombrables don-
nées qu' il avait déjà accumulées. Deuxième loi : La droite joignant le centre du Soleil
au centre d'une planète balaie des aires égales en des
intervalles de temps égaux.

Troisième loi : Le carré de la période de révolution

d'une planète est proportionnel au cube du demi
grand axe de l'ellipse.

Nt::wlon fut le premier à donner une véritable démons-

tration des lois de Kepler ; cela consacrait en quelque
sorte ses nouvelles théories sur le calcul différentiel et
intégral et sur sa loi de la gravitation universelle.

Dans les démonstrations qui vont suivre, nous utilise-

rons la seconde loi de Newton sous sa forme vecto-
rielle : la somme des forces appliquées sur un objet de
masse m est égale à sa masse multipliée par le vecteur
accélération :

Figure 11

10 - Bulletin AMQ, Vol. XLIII, n° 2, mai 2003

Le vecteur accélération se calcule en dérivant deux
P=r(t)cosO(t)•Ï +r(l)sinO(t),],
fois le vecteur position :
v = [r'(t) cos8- r(1)sin O·B'(t )]· i +[r'(t ) sinO+ r(t )cosO•B'(t ) J-J,
ii =[r' cosO - r' sin O·B' -r' sinO·B' -rcosB•(B')2 - rsinO·O"}ï
+[r'sinO+r' cos0 ·0' + r' cosO•B'-rsinO•(O')' + rcos B•B'} J
►
p ''
1
:
1
)' Le vecteur radial unitaire est : ü == cose •7 +sine·] et
r ' le vecteur normal unitaire, n= - sine· i = cose •]. On
peut alors réécrire le vecteur accélération ainsi :

2
â =[r• - r(B') } (cos8·i +sin0,})+[2r'O' +rO' ],(-sinO·i +cosO•])

Figure 2
( )) â = [r' -r(O')'}ü + [2r'O' +rO'J ·ii •

P = x(l)-T + y(t) · J ,
Comme la seule force en présence est radiale, la com-
v=x'(t)·T+ y '(t)·] , posante normale est nulle :
dr d0 d 18
ëi = x"(t)-T+y"(t)-]. 2- - ·- +r ·- 2 = O.
dt dt dt
Multiplions l'équation par r:
1. Dé monstration de la seconde loi de Kepler
dr d0 d 20
Nous allons commencer par démontrer la seconde loi 2r--·-+r2 •- =
2
O.
dt dt dt
de Kepler et ce comme corollaire du théorème sui-
vant: Cela devient une dérivée exacte :

Théorème 1. Si une masse se dé place dans un champ !!...[ri.d0] = O_

dt dt
de force central, en l 'absence de toute autre force, le
vecteur position décrit des aires égales en des temps On en déduit que :
égaux. d0
r 2 •- = k où k est une constante.
dt
Démonstration: On place l'origine du système de ré-
férence au centre de force et on utilise les coordonnées Désignons par A l' aire balayée par le vecteur position
polaires. pendant l' intervalle de temps

Figure 3 Figure 4

Bulletin AMQ, Vol. XLIII, n° 2, mai 2003 - 11

Un élément d'aire en coordonnées polaires est : la densité des lignes de force est inversement propor-
~ ~ r · 110 · & , tionnelle au carré de la distance. Entre autres, cela est
vrai pour la gravitation, pour les champs électriques et
pour l'intensité d.'une source lumineuse.

La première loi de Kepler sera une conséquence du

théorème suivant :

Théorème 2. Si un objet se déplace dans un champ de

Figure 5 force central inversementproportionnel au carré de la
De telle sorte que : distance entre le centre del 'objet et le centre deforce,
alors la trajectoire suivie par / 'objet est une conique.
A= fJrdrd0
~
Démonstration : Cette démonstration est la suite de la
démonstration précédente. La force de gravité est ra-
diale et dirigée vers le centre; à partir de la formule du
vecteur accélération obtenue en ( 1), nous obtenons :
Mais comme r 2d0 = kdt, il s'ensuit que
1
I ' k
A=-
2 f kdt = -(t
2 2 - t)
1 •
,,
Et, comme r i •-d() = k , 1·1 s'ensuit
. que :
Le vecteur position balaie des aires égales en des dt
temps égaux.
(2)
C.Q.F.D.

Nous voyons que la seule condition utilisée dans cette Comme le but est de trouver la trajectoire décrite par
démonstration est que le champ de force soit central et l'objet, il faut trouver r en fonction de 0 plutôt qu'.en
qu' il n' y ait pas d'autres forces en présence. fonction du temps. De plus, nous effectuerons le chan-
!
gement de variable u = - .
2. Démonstration de la première loi de Kepler r

Pour démontrer les deux autres lois Kepler nous au- du dr du du d0 . , .

rons besoin de la loi de la gravitation universelle de Comme - = - - 2 ·- et - = - - - ils ensuit
dt r dt dt d0 dt '
Newton :
que:

En fait, il suffit que la force soit inversement propor-

tionnelle au carré de la distance. Cette idée était déjà
dans l' air depuis un certain temps, ce n'est pas Newton
qui l'a inventée3. Elle s' explique faci lement pour un
champ de force qui rayonne dans l'espace à partir d' un
point central : comme la surface d'une sphère est 4nr 2 ,

12 - Bulletin AMQ, Vol. XLIII, n°2, mai 2003 - 17

En portant cette expression dans l'équation différen-
tielle (2), nous obtenons :

(3) (4) r(0)= K .

Ki cos0+ k2
C'est une équation linéaire d'ordre 2 non-homogène
de la forme ayw + by' + cy = c. Ce type d 'équation est Comme nous allons le voir, cela représente l' équation
étudié en Calcul III 4 • On commence par trouver la so- d 'une conique en coordonnées polaires.
d2u
lution générale del' équation homogène dei + u = 0 .
C.Q.F.D.

Pour ce, on calcule le discriminant du polynôme carac- 3. Équation des coniques e n coordonnées
téristique : t:i. = b 2 - 4ac = -4 . Comme sa valeur est polaires
négative, on aura comme solution Les coniques peuvent toutes être définies en fonction
d 'un point F (d it le foyer), d'une droite !),, (dite droite
directrice) et d' un nombre positif e (appelé excentri-
cité):
où
b d(P, F)-=:- e-d(P, t:i.)
a.=-- ,
2a
qui donne O dans ce cas-ci, et

✓4ac- b 2
~= --
2a

qui donne 1. On obtient donc

u(,0)= K 1sin0 + Ki cos0

►
pour l' équation homogène. Il suffit de trouver une so-
lution particufü:n:: à l'équation non-homogène et de les
additionner. Ainsi,

Figure 6

est la solution de l' équation (3).

Dans le cas d' une parabole, e = 1 ; dans le cas d' unV
Si nous plaçons l' axe des coordonnées polaires de telle
sorte que r(0) soit minimal lorsque 0 = 0 , u(0) de-
ellipse, e =:.a < 1 et dans le cas d ' une hyperbole,
vrait être alors maximal, ce qui donne K 1 = 0 et, puis-
que uw(O) = -K 2 <0, K2 > 0. Ainsi,
e =:.a > 1
5
. Le cercle est un cas limite où la droite direc-
trice est à l' infini et l' excentricité est nulle.

Bulle1in AMQ, Vol. XLIII, n° 2, mai 2003 - 13

r

Trouvons maintenant l'équation d' une conique en

coordonnées polaires en plaçant le pôle au foyer et en
prenant l'axe perpendiculaire à la droite directrice.

d(P, F) = e·d(P, t,.)

On voit que

d(P, F)= r et d(P,/',.)= p - x = p - rcos0 ;

d'où,

r = e ·( p-rcos0) = ep- ercos0

e·p
(5) r= .
l+e-cos0
Figure 7

L'équation (4) peut se ramener à cette forme; il s' agit

donc d'une conique. Par contre, si la vitesse de l'objet est comprise entre
certaines valeurs, il entrera en orbite sur une ellipse au-
1 tour de l' astre, sa vitesse déterminera les valeurs de a
(6) r(0) = K
et de b. Dans le cas d ' un satellite, cela déterminera
K 2cos0+,:i l' altitude à laquelle il tournera.
1
= - -, - - - - - ---,-
-K . ( 1+ -K2k2
- cos0)

•
k2 K
k2
=- --"K-=-----
K 2k 2
l + ~cos0
Figure 8

2
K k2 . Ce 1a comp1'ete 1a d'emonstrat1on
avec e = ~ . du Finalement, si sa vitesse n' est pas assez grande, il s'é-
crasera sur l' astre en décrivant une parabole. C' est ce
théorème 2.
qui est arrivé à deux reprises aux satellites de la NASA
sur Mars.
Pour achever la démonstration de la première loi de
Kepler, examinons les situations possibles. Si un objet
entre dans un champ gravitationnel quelconque, il
peut se produire trois cas : si sa vitesse linéaire est suf-
fisamment grande, il sera dévié par le champ de force
sur une hyperbole et poursuivra sa route sur une des
droites asymptotiques.

Figure 9

14 - Bulletin AMQ, Vol. XLIII, n°2, mai 2003

Comme la trajectoire des planètes autour du Soleil est Dans une ellipse, l'excentricité est égale à c/a, où c est
fermée, on peut conclure qu' il s'agit d'ellipses, ce qui la distance focale et c 2 = a 2 - b 2 • Ainsi,
termine la démonstratio n de la première loi de Kepler.

4. Dé monstration de la troisième loi de Kepler

Démontrons maintenant la troisième loi, en poursui-
vant la même démarche que dans les deux autres dé- d'où
monstrations.

Supposons que la trajectoire soit une ellipse; l' équa-

tion (6) peut alors s'écrire sous la forme:

k2
r(0)= - -- -
K(I + ecos0)

La distance minimale de r sera atteinte lorsque le co- Soit T la période d' une révolutio n complète; comme
sinus atteint son maximum à 0 = 0 ; d 'où il a été démontré précédemment, l'aire balayée par le
vecteur position entre le temps / 1 et le temps / 2 est
k2 donnée par la formule
r m,n - -- -
- K (I +e)

et la valeur maximale sera atteinte lorsque le cosinus

atteint son minimum à 0 = 1C :
L'aire totale de l'ellipse étant nab , nous obtenons:
k2
r - -- - 1
rn,n - K(l- e) · nab=-kT
2 '

Soient P le point le plus rapproché du foyer Fi et P' 2nab

le plus éloigné ; T =-
k '

P' p
ce qui démontre la troisième loi de Kepler.

On peut faire refaire à des étudiants de Calcul Ill toutes

Figure 10 ces démonstrations (en donnant évidemment beau-
coup d ' indications) comme travail de session dans le
cadre de l' activité d' intégration. C' est un bel exemple
de lien entre la physique et les mathématiques. En leur
Ym,n + r max k2
faisant faire une introduction historique plus détaillée,
a= 2 = K(l -e2 ) •
on touchera aussi à l' histoire des sciences. ■

Bulletin AMQ, Vol. XLlll, n°2, mai 2003 - 15

 Notes [2] Piskounov, N . ( 1980). Calcul différentiel et inté-
1
Photo tirée du site internet [Link] gral, Tomes let Il. Moscou, Éditions Mir.
.[Link]/[Link] ; elle provient du musée
autrichien der Sterwarte Kremsmünster. [3] Robert, S . ( I 992). Calcul différentiel et intégral.
Saint-Jean (CEGEP Saint-Jean-sur-Richelieu, Les
2
On peut consulter la référence[!] pour en savoir plus productions FRP.
sur Kepler.
(4] Robert, S. (1995). Algèbre, trigonométrie et géo-
3
On peut consulter [6], un magnifique livre consacré à métrie analytique. Saint-Jean (CEGEP Saint-Jean-
la vie de Newton. sur-Richelieu, Les productions FRP.

4
Pour plus de détails sur cette méthode, on peut con- (5] Tie rny, J. A. ( 1985). Differential equations. Boston
sulter la référence [3], pages 307 et suivantes. (MA), Allyn a nd Bacon inc.

5
Le lecteur aura deviné que a représente le demi-axe [6] Westfall, R. (1994). Newton. Paris, Flammarion.
horizontal et b, le demi-axe vertical, dans le cas d'une
ellipse. P our les équat io n s des coniques e n
coordonnées cartésiennes et pour une étude plus Serge Robert
détaillée de celles-ci, on peut consulter la référence Cégep Saint-J ean-sur-Richelieu
[4]. 30, boulevard du Séminaire, C.P. 1018
Saint-Jean-sur-Richelieu (Québec)
Références bibliographiques 13B 7B1
[1] Koyré, A. ( 1961 ). La révolution astronomique, Co-
pernic, Kepler, Borelli. Paris, Hermann. [Link]@ [Link]

Souscription au Fonds Maurice-L'Abbé pour les camps mathématiques

Nom: _ __ __ _ _ __
Oui! Je désire contribuer au financement des camos mathématiaues. Adresse : _ _ __ __ _ __

□ 20$ □ 30$ □ 50 $ □ 100 $ Code postal: _ __ _ _ _ __ _

AUTRES

□ PAR CHÈQU E À L'ORDRE DE L'AMQ Pour 20 $ ou plus, ou sur demande, vous

recevrez un reçu pour fin d'impôt.
□ VISA O MASTER CARO Date d'expiration : _ __ NE: 12 577 5858 RR 0001
NO. DE LA CARTE : _ __ _ _ _ _ _ _ _ _ __
D Je désire recevoir !.ln reçu pour fin d'impôt
SIGNATURE : _ _ _ _ __ _ _ _ _ _ _ __

7 400, boulevard Saint-Laurent, bureau 257, Montréal (Québec) H2R 2Y1 - 514-278-4263

16 - Bulletin AMQ, Vol. XLIIl, n° 2, mai 2003