Guide de Ion Des Astres
Guide de Ion Des Astres
Guide de
localisation
des astres
Guide
de localisation
des astres
Christian Gentili
ISBN : 978-2-7598-0059-9
Imprimé en France
Avant-propos .......................................................................................... 15
Chapitre 1
Mécanique céleste
1 Structure cosmique ...................................................................... 20
1.1 Notre système solaire ......................................................................... 20
Soleil ..................................................................................................... 20
Planètes ................................................................................................. 20
Astéroïdes et comètes ............................................................................ 21
Autres corps........................................................................................... 22
Satellites naturels ................................................................................ 22
Satellites artificiels .............................................................................. 22
1.2 Galaxies ................................................................................................ 23
Notre Galaxie ........................................................................................ 23
Amas galactiques .................................................................................. 24
2 Durées ................................................................................................ 25
2.1 Jours et années .................................................................................... 25
Les différents jours ................................................................................ 25
Jour sidéral ......................................................................................... 25
Jour solaire ......................................................................................... 25
Les différentes années ........................................................................... 27
Année sidérale .................................................................................... 27
Année tropique ................................................................................... 27
Année julienne, calendriers julien et grégorien ..................................... 28
5
Guide de localisation des astres
3 Distances .......................................................................................... 30
3.1 Unité astronomique ............................................................................ 31
3.2 Année-lumière et parsec .................................................................... 33
Exercices .............................................................................................. 36
E1-1 Jour Julien ........................................................................................ 36
E1-2 Périodes des planètes ....................................................................... 37
Chapitre 2
2 Définition bifocale........................................................................ 49
2.1 Cas de l’ellipse .................................................................................... 49
2.2 Cas de l’hyperbole .............................................................................. 49
Exercices .............................................................................................. 52
E2-1 Géométrie de l’ellipse en fonction du demi-grand axe
et de l’excentricité............................................................................... 52
E2-2 Équation commune aux coniques en fonction du paramètre
et de l’excentricité............................................................................... 54
6
Sommaire
Chapitre 3
exercices ............................................................................................... 68
E3-1 Vitesses en certains points de l’ellipse ; relations déduites .............. 68
E3-2 Vitesses des planètes du système solaire........................................... 70
Chapitre 4
7
Guide de localisation des astres
Exercices .............................................................................................. 82
LI KI
E4-1 Accélération centrale de la forme Γ = − 3 I .............................. 82
ρ LI C I
2
E4-2 Cas particulier : accélération centrale de la forme Γ = − 3 I ... 84
ρ
Chapitre 5
Anomalie et orbite
1 Anomalies ....................................................................................... 88
1.1 Temps du passage ............................................................................... 88
Définitions ............................................................................................. 88
Équation de Kepler ............................................................................... 89
1.2 Relations entre anomalies vraie et excentrique .............................. 89
Rayons ρ et ρ' ........................................................................................ 89
Relations trigonométriques liant θ et ϕ ................................................ 89
Application à l’expression de la vitesse et de l’accélération ................ 90
8
Sommaire
Exercices .............................................................................................. 97
E5-1 Temps de l’intersection d’une comète avec l’écliptique ................. 97
E5-2 Durée et date des saisons................................................................. 102
E5-3 Réflexion sur l’ellipse ....................................................................... 104
Chapitre 6
3e loi de Kepler
et référentiels
1 Lois de Kepler pour un centre attractif fixe .................... 110
1.1 Rappel des hypothèses....................................................................... 110
e
1.2 3 loi de kepler..................................................................................... 111
Formulation ........................................................................................... 111
Expressions numériques ........................................................................ 112
9
Guide de localisation des astres
Chapitre 7
Énergie et vitesse
selon la trajectoire
1 Caractérisation énergétique des trajectoires .................. 128
1.1 Énergie potentielle, cinétique et mécanique ................................... 128
Énergie potentielle EP ........................................................................... 128
Énergie cinétique EC ............................................................................. 129
Énergie mécanique E ............................................................................ 130
Énergie fournie EF ................................................................................. 131
1.2 Trajectoire déduite de la loi de conservation de l’énergie ............. 132
Chapitre 8
10
Sommaire
Chapitre 9
11
Guide de localisation des astres
Chapitre 10
Coordonnées équatoriales
et écliptiques
1 Voûte céleste .................................................................................. 220
1.1 Constellations ..................................................................................... 220
1.2 Magnitude ........................................................................................... 222
1.3 Nom des étoiles .................................................................................. 222
1.4 Catalogues des étoiles et objets non stellaires............................... 223
12
Sommaire
Chapitre 11
13
Guide de localisation des astres
14
Avant-propos
Notre univers ne cesse de nous fasciner et nombreux sont ceux qui consacrent
leur vie professionnelle ou leurs loisirs à son observation, afin d’en percer les
mystères…
Il y a de par le monde des milliers de clubs et de sociétés d’astronomie.
Parallèlement, les programmes spatiaux d’exploration de notre système
solaire sont toujours actifs et sont passés du stade de compétition à celui, très
fructueux, de coopérations internationales orchestrées par des organisations
prestigieuses telles que la NASA, l’Institut russe de la recherche spatiale,
l’Agence spatiale européenne…
Pointer un télescope sur un astre ou une antenne sur un satellite, déduire
sa position de celle des étoiles et connaître l’instant du crépuscule, définir
le mouvement d’une sonde dans le milieu interplanétaire, déterminer à un
instant donné la position d’une comète ou d’un satellite artificiel sur l’or-
bite assignée par sa mission sont autant de problèmes qui, nous l’espérons,
trouveront leur solution dans le présent ouvrage. Dans cet objectif, et afin
de s’adresser au plus grand nombre, nous avons adopté une démarche très
progressive, guidée également par le souci d’éviter au lecteur d’avoir à
rechercher dans d’autres sources les supports mathématiques et physiques
nécessaires à la compréhension.
Le milieu au sein duquel nous évoluerons, le système solaire, est présenté
dans un chapitre introductif où nous précisons un certain nombre de gran-
deurs astronomiques fondamentales. À ce stade, il convient de saluer le travail
des astronomes, en particulier dans le domaine de la mesure et des unités ;
15
Guide de localisation des astres
16
Avant-propos
17
This page intentionally left blank
Chapitre 1
Mécanique céleste
1 Structure cosmique
Le Soleil est plutôt une petite étoile ; certaines sont quelques dizaines de fois
plus lourdes. En dessous d’une masse critique, une étoile ne peut pas naître
car son cœur ne peut atteindre les températures nécessaires à l’amorçage de
réactions nucléaires.
Planètes
Par rapport à une étoile, une planète est un astre mobile dont la brillance est
due à la réflexion de la lumière de l’étoile voisine, le coefficient de réflexion
étant l’albédo.
20
1. Mécanique céleste
Les planètes ont des masses bien inférieures à celles des étoiles. Dans notre
système solaire, la plus grosse planète, Jupiter, est 318 fois plus massive que
notre Terre pour laquelle :
MT = 5,973 64 1024 kg
Jupiter est mille fois plus léger que notre Soleil.
Par rapport à la masse solaire, le Bureau des longitudes (voir bibliographie)
donne :
MT = 3,003 489 6 10-6 MS MJ = 9,545 942 9 10-4 MS
Les planètes, du mot grec signifiant « vagabond », orbitent autour du Soleil
selon des trajectoires faiblement elliptiques, sauf pour la plus proche, Mercure,
et la plus lointaine, Pluton. Les ellipses sont décrites à peu près dans un même
plan, sauf pour Pluton !
Le plan dans lequel orbite la Terre est appelé écliptique car c’est dans ce plan
qu’ont lieu les éclipses en relation avec notre planète.
On distingue les planètes inférieures, qui sont plus proches du Soleil que la
Terre, et les planètes dites supérieures, qui en sont plus éloignées.
Le classement suivant est plus intéressant :
– en s’éloignant du Soleil, Mercure, Vénus, Terre et Mars forment un groupe
d’astres constitués d’un noyau métallique entouré d’un manteau rocheux ; ce
sont les planètes telluriques ;
– au-delà, Jupiter, Saturne, Uranus et Neptune constituent le groupe des
planètes joviennes (nom dérivé de Jupiter). Elles sont beaucoup plus volu-
mineuses que les planètes telluriques, mais caractérisées par une plus faible
densité due à une fraction importante d’éléments légers, tels que le dihydro-
gène et l’hélium ; ainsi la surface de ces astres est fluide, et leur limite avec
l’atmosphère mal définie.
Ces quatre planètes possèdent des anneaux, mais seul celui de Saturne est
spectaculaire. La planète Neptune quant à elle se caractérise par son axe de
rotation incliné de 98° par rapport à la perpendiculaire au plan de son orbite,
d’où son surnom de « planète couchée ».
Pluton n’est pas une planète jovienne ; précisons qu’une partie de sa trajec-
toire est plus près du Soleil que celle de Neptune, mais il ne peut y avoir de
risque d’interférence vu la forte inclinaison de l’orbite par rapport à l’éclip-
tique (17,2°). Notons que la période de révolution de Pluton est 3/2 de celle de
Neptune (respectivement 247,7 et 164,8 années).
Astéroïdes et comètes
En dessous d’une certaine taille, on parle d’astéroïdes et non de planètes :
Pluton, d’un diamètre à peine supérieur à 2200 km est une petite planète ;
21
Guide de localisation des astres
Cérès, d’un diamètre de 940 km, est le plus gros astéroïde. Actuellement, on
connaît les orbites de quelques milliers d’astéroïdes ; la grande majorité de
ceux-ci ont des orbites quasiment circulaires confinées entre Mars et Jupiter,
qui constituent la ceinture principale. Toutefois certaines orbites ont une
excentricité notable, faisant que les périhélies peuvent être proches de Mars
voire de la Terre… À la différence des planètes parfaitement sphériques, les
astéroïdes de petite taille ont tendance à être d’une forme irrégulière.
Bien au delà de l’orbite de Pluton existe un nuage de petits astéroïdes consti-
tuant le nuage de Oort, du nom de l’astronome qui en fit l’hypothèse.
Certains corps quittent ce nuage et adoptent des orbites héliocentriques très
elliptiques : ce sont les comètes. Elles sont constituées d’un mélange de pous-
sières et de glace. Si la comète s’approche près du soleil, c’est-à-dire dès que
l’ordre de grandeur de sa distance est compris entre Mars et Jupiter, son noyau
commence à s’évaporer, donnant naissance à une queue spectaculaire dirigée
vers l’extérieur de la trajectoire, comme l’avaient noté les Chinois plus de
mille ans avant J.-C.
Autres corps
On considère ici les corps orbitant autour de planètes dont la masse leur est
bien supérieure ; on les dénomme satellites. On distingue les satellites natu-
rels et les satellites artificiels.
Satellites naturels
Les planètes inférieures n’ont pas de satellites. La Terre possède un seul satel-
lite, la Lune, dont la masse représente environ 1,2 % de celle de notre astre.
Les deux satellites de Mars sont minuscules et s’apparentent à des astéroïdes.
Sur les 16 satellites de Jupiter, 12 sont petits et ressemblent à des astéroïdes
également, mais les quatre autres, découverts par Galilée (1564-1642), ont
une dimension et une masse de l’ordre de grandeur de celle de la Lune. Ce
sont, en s’éloignant de Jupiter : Io, Europe, Ganymède et Callisto.
Saturne a de nombreux satellites assez petits, sauf Titan qui est plus gros que
la Lune. Le satellite Triton a la particularité de tourner autour de Neptune dans
le sens rétrograde.
Notre système solaire compte ainsi 16 objets planétaires de diamètre supérieur
à 2500 km.
Satellites artificiels
C’est au cours des années 1957-1958 que les premiers satellites artificiels
de la Terre commencèrent à peupler notre espace : le 4 octobre 1957, le
premier spoutnik fut mis en orbite par les Soviétiques ; entre 1958 et 1959,
22
1. Mécanique céleste
1.2 Galaxies
Notre Galaxie
Les galaxies sont principalement constituées d’une concentration d’étoiles
dont l’origine est la force gravitationnelle. Le nombre d’astres ainsi rassem-
blés peut atteindre plusieurs centaines de milliards.
On distingue trois principaux groupes de galaxies : les galaxies irrégulières,
les galaxies elliptiques et les galaxies spirales. Ces dernières sont les plus
nombreuses et notre Galaxie, qui compte environ cent milliards d’étoiles, est
de ce type.
Vue de face, la Galaxie a la forme d’un disque où les étoiles sont concentrées
essentiellement au centre ainsi que dans quatre bras spiraux. Le Soleil est dans
le bras d’Orion, considéré comme une ramification, ou plutôt une protubé-
rance du bras du Sagittaire.
Le disque a une certaine épaisseur, et, observée par la tranche, la Galaxie
montre un renflement symétrique au centre appelé bulbe ; d’autre part, elle
est vue sur la voûte céleste comme une bande ou traînée laiteuse, d’où elle tire
son nom de Voie lactée.
23
Guide de localisation des astres
≈ 26 000 al
≈ 100 000 al
Amas galactiques
Dans l’univers, les galaxies sont en général groupées en amas, dits amas galac-
tiques, de plusieurs milliers d’éléments. La Voie lactée appartient à un petit
amas galactique appelé Groupe local. Ses plus proches voisines sont deux
petites galaxies satellites, appelées le Petit Nuage de Magellan et le Grand
Nuage de Magellan. Un peu plus loin, à deux millions d’années-lumière (voir
paragraphe 3), on rencontre la galaxie d’Andromède (M31) très semblable à
notre Voie lactée, ainsi qu’une autre galaxie spirale, M33. Généralement, les
galaxies sont uniquement identifiées par un numéro dans le catalogue d’astro-
nomie où elles sont répertoriées.
24
1. Mécanique céleste
2 Durées
25
Guide de localisation des astres
m
m
26
1. Mécanique céleste
Année tropique
L’axe de rotation de la Terre est incliné de 23° par rapport à la perpendiculaire à
l’écliptique ; il n’est pas pour autant fixe par rapport aux étoiles : il décrit, dans
le sens rétrograde, un cône de demi-angle 23°. La période de cette précession
est légèrement inférieure à 26 000 ans (on retient plutôt 25 800 ans). Solstices
et équinoxes vont donc dériver avec la même période : c’est le phénomène dit
de précession des équinoxes.
Si on veut que l’année sidérale représente le cycle des saisons, il faut la corriger
de la précession ; on obtient alors l’année tropique. La révolution et la préces-
sion étant de sens opposés, la correction annuelle de 365,2564/26 000 jours
est négative :
27
Guide de localisation des astres
2.2 Périodes
Périodes synodique et sidérale
• La période synodique correspond au temps de réapparition d’une même
configuration de deux planètes par rapport au Soleil.
Les configurations caractéristiques sont :
28
1. Mécanique céleste
– la conjonction, lorsque les deux planètes sont alignées avec le Soleil d’un
même côté ;
– l’opposition, lorsque les deux planètes sont de part et d’autre alignées avec
le Soleil ;
– la quadrature, lorsque les directions des planètes avec le Soleil forment un
angle droit.
• La période sidérale est le temps TP mis par une planète quelconque pour
accomplir une révolution complète autour du Soleil par rapport à un obser-
vateur lié aux étoiles fixes. Elle correspond à l’année sidérale pour la planète
Terre.
29
Guide de localisation des astres
360° + θ
S T PS
Remarque
Un autre raisonnement consiste à considérer les vitesses de chacune des
planètes :
- La vitesse de la planète supérieure est de 360/TPs degrés par jour ; pendant
la période synodique SPs le nombre de degrés décrit est
360° θ S
θ = SPs → = Ps
TPs 360° TPs
- La vitesse de la Terre est de 360/TT degrés par jour ; pendant la même
durée SPs le nombre de degrés décrit est
360° SPs S
θ + 360° = SPs → + 1 = Ps
TT TPs TT
Divisant par la période synodique, il vient :
1 1 1
+ =
TPs SPs TT
On retrouve la relation valable pour les planètes supérieures.
3 Distances
30
1. Mécanique céleste
31
Guide de localisation des astres
orbite lunaire :
entrée pénombre - dist. 357 694 km à 3 h 40
1 h 06 - diam. ombre 1°,548
- diam. pénombre 2°,624
entrée ombre
2 h 03
début totalité
3 h 14
maximum éclipse Lune Terre Soleil
3 h 40
fin totalité
4 h 06
sortie ombre
5 h 17
sortie pénombre
6 h 15
32
1. Mécanique céleste
33
Guide de localisation des astres
2θ
Étoile
Terre D
position 1
Soleil a T
Terre
position 2
34
1. Mécanique céleste
35
Guide de localisation des astres
exercices
E1-1 Jour Julien
Le jour julien (JJ) est le nombre de jours écoulés depuis le 1er janvier
4712 avant J.-C. à 12h TU. On peut calculer le JJ correspondant à une date
donnée à l’aide de la formule suivante :
JJ = P.E.[365,25 (A + 4716)] + P.E.[30,6001 (M + 1)] + J + B + T − 1524,5
P.E. = Partie Entière
J,M,A = date (J jour du mois, M numéro du mois, A année)
– Si M = 1 ou 2, on remplace A par A -1 et M par M + 12
– Si J,M,A est une date du calendrier julien, B=0
– Si J,M,A est une date du calendrier grégorien, B = 2 − C + P.E.[C/4],
où C = P.E.[A/100]
L’heure UTC (ou TU) est exprimée sous la forme HH MM SS, et
T = HH/24 + MM/1440 + SS/86400
On définit également le jour julien modifié (JJM ou MJD pour Modified
Julian Day) qui se déduit du JJ en retranchant 2 400 000,5 j.
1. Calculer le JJ correspondant aux dates suivantes :
1er janvier an 0 12h TU ; 15 octobre 1582 0h TU ; 4 septembre 1983 0h TU
2. Calculer le JJ correspondant aux 1er janvier 1001 et 2001 à 0h TU, justifier
la différence par rapport à l’année de 365,25 j.
3. À quelle date correspond 2 400 000,5 JJ ?
4. Calculer le JJ et le MJD à la date du 28 février 2003 à 23 h 59 min
59,9 s TU.
solutions
1. 1er janvier an 0, 12 h TU → 1 721 058 JJ
15 octobre 1582, 0 h TU → 2 299 160,5 JJ
4 septembre 1983, 0 h TU → 2 445 581,5 JJ
2. 1er janvier 1001, 0 h TU → 2 086 673,5 JJ
1er janvier 2001, 0 h TU → 2 451 910,5 JJ
exercices
36
1. Mécanique céleste
Remarque
– L’époque J2000,0 , c’est-à-dire la date référence du 1er janvier 2000
à 0 h TU correspond à 2 451 544,5 JJ.
– L’époque J2000 , c’est-à-dire la date référence du 1er janvier 2000
à 12 h TU correspond à 2 451 545,0 JJ.
37
Guide de localisation des astres
solutions
1.
a)
période sidérale pér. synodique
Planète années jours jours
Mercure 0,2409 88 115,9
Vénus 0,6152 224,7 583,9
Terre 1,0000 365,2564 --------
Mars 1,8809 687 779,9
Jupiter 11,8631 4 333 398,88
Saturne 29,4593 10 760 378,09
Uranus 84,0219 30 689 369,66
Neptune 164,7693 60 182 367,49
b)
– Pour Mercure, le nombre de jours séparant les deux passages est donné par
la différence des dates exprimées en jours juliens :
2 454 048,403 JJ – 2 452 766,828 JJ soit 1281,575 jours .
Compte tenu de la période synodique de 115,9 jours, cela correspond à la
onzième conjonction, et à 14,5 périodes de rotation de Mercure, et 3,5 de
la Terre.
Mercure a bien effectué 11 rotations de plus que la Terre.
– On calcule de même pour Vénus :
2 456 084,555 JJ – 2 453 164,847 JJ soit 2919,708 jours .
Compte tenu de la période synodique de 583,9 jours, cela correspond à la
cinquième conjonction et à 13 périodes de rotation de Vénus, et 8 de la Terre.
Vénus a bien effectué 5 rotations de plus que la Terre.
2.
exercices
Par rapport à Jupiter, Mars est une planète inférieure. La formule se note :
1 1 1
= − SM = 816,45 j = 2,235 années
SM TM TJ
Par rapport à Jupiter, Saturne est une planète supérieure. La formule se note :
1 1 1
= − SS = 7254,3 j = 19,861 années
SS TJ TS
38
1. Mécanique céleste
3.
• Supposons, pour illustrer, qu’au départ on ait l’alignement Soleil-Terre-
Lune ; après un mois sidéral TL il faudra encore attendre quelques jours pour
arriver à la pleine Lune car la Terre aura tourné autour du Soleil de :
TL
360° = 26°, 9285
TT
La Lune a parcouru 360° en TL jours ; les 26°,9285 seront parcourus en :
TL T T2
× 360° L = L = 2, 0437 j
360° TT TT
• Effectuons la deuxième itération :
Pendant ces 2,0437 j la Terre s’est déplacée de :
360° TL2 T2
× = 360° L2 = 2°, 0143
TT TT TT
Pour parcourir ces 2°,0143 , la Lune mettra le temps supplémentaire :
TL T2 T3
× 360° L2 = L2 = 0,1529 j
360° TT TT
Correction à ajouter : ΔTL = 2,1966 j
• Ainsi après deux itérations, la valeur du mois synodique est :
TL2 T3
SL2 = TL + + L2 ; soit SL2 = 29,5183 j
TT TT
⎡ T T2 T3 ⎤
SL = TL ⎢1 + L + L2 + L3 + ..... ⎥
⎣ TT TT TT ⎦
La raison de cette progression géométrique est inférieure à l’unité, d’où :
1 1
SL = TL =
TL 1 1
1− −
TT TL TT
SL = 29,5306 j
Remarques
1– Phases de la Lune et lunaisons
L’alignement Soleil-Terre-Lune correspond à la pleine lune (PL) ; la nou-
velle lune (NL) correspond à l’alignement Soleil-Lune-Terre. À partir de
cette dernière, les phases intermédiaires sont le premier croissant, le pre-
mier quartier (PQ), la lune gibbeuse croissante ; puis après la pleine lune,
39
Guide de localisation des astres
On déduit que sur l’année 2004 la valeur moyenne des 12 lunaisons comp-
tées à partir de la pleine lune est de 29,4980 j. Dans le tableau ci-dessous,
les lunaisons sont cette fois calculées entre deux nouvelles lunes consécu-
tives et la valeur moyenne est de 29,5621 j !
Nouvelle lune Mois lunaire Nouvelle lune Mois lunaire
Date heure TU j, h, min jours Date heure TU j, h, min jours
21/01 21 h 05 min 17/07 11 h 24 min
29 j 12 h 13 min 29,5090 29 j 14 h 00 min 29,5833
20/02 09 h 18 min 16/08 01 h 24 min
29 j 13 h 23 min 29,5576 29 j 13 h 05 min 29,5451
20/03 22 h 41 min 14/09 14 h 29 min
29 j 14 h 40 min 29,6111 29 j 12 h 19 min 29,5132
19/04 13 h 21 min 14/10 02 h 48 min
exercices
40
1. Mécanique céleste
exercices
41
This page intentionally left blank
Chapitre 2
Rappels sur les coniques
Une conique est l’ensemble C des points M du plan tels que le rapport des
distances à un point fixe, le foyer F, et à une droite fixe, la directrice (D), est
constant. H étant la projection de M sur (D), on a :
MF
= e (e est l’excentricité).
MH
y
H
M(x,y)
(D) ax + by + c = 0
F(x0 ,y0 )
O x
M ∈ C ⇔ MF2 – e2 MH2 = 0
(ax + by + c)2
⇔ (x – x0)2 + (y – y0)2 – e2 = 0
a 2 + b2
Équation de la forme :
Ax2 + 2Bxy + Cy2 + 2Dx + 2Ey + F = 0
44
2. Rappels sur les coniques
y
(D)
M H
F i D x
k
45
Guide de localisation des astres
Cas de l’ellipse
La figure 2.3 fait apparaître en plus de F et (D), le foyer F’ et la directrice (D’)
symétriques par rapport au centre O caractérisé comme suit :
Le foyer F est l’origine ;
A et A' ont pour coordonnées A (ρ0 , 0) et A'(−ρπ , 0) ;
ρ0 − ρπ
les coordonnées du milieu O de AA' sont O ( , 0).
2
On déduit :
A ' A = ρ0 + ρπ F ' F = ρπ − ρ0
Ox est l’axe focal ; Oy l’axe non focal.
y
B M
H’ H
ρ
θ
D’ A’ F’ O F A D x
B’
(D’)
(D)
46
2. Rappels sur les coniques
Détermination de a, c, b
ρ0 + ρπ p
Classiquement, on pose a = OA a = =
2 1 − e2
ρ − ρ0 ep
c = OF c = π =
2 1 − e2
On déduit c/a = e ρ 0 = a – c ; ρπ = a + c
p
On introduit b tel que b= a 2 − c2 b = ρ0ρπ =
1 − e2
On déduit b/a = 1 − e 2 , c’est le coefficient d’aplatissement.
La distance 2a représente le grand axe ; la distance 2b le petit axe.
La distance 2c est la distance focale.
a apparaît comme la moyenne arithmétique des rayons au péricentre et à
l’apocentre.
b apparaît comme la moyenne géométrique de ces mêmes rayons.
Le rapport moyenne géométrique sur moyenne arithmétique vaut 1 − e 2 .
Dans la littérature, la distance moyenne au centre est tout simplement le demi-
grand axe de l’ellipse. On pourrait également retenir comme rayon moyen ρm
la racine carrée du produit du demi-grand axe par le demi-petit axe :
ρm = (ab)1/2
Cela se justifierait par le fait que le cercle de rayon ρm a la même surface que
celle de l’ellipse.
Détermination de p, k, d, e
Considérons le paramètre p et la distance foyer directrice k :
p = a(1−e2) ⇒ p = b2/a
k = p/e ⇒ k = b2/c
Soit d la distance de O à (D), d = OD = OF + FD = c + k, d’où
d = a2/c
Rappel e = c /a
Précisions sur le vocabulaire
Nous avons utilisé les noms péricentre et apocentre pour désigner la position
respectivement la plus proche et la plus éloignée du foyer considéré de l’el-
lipse. Si le foyer correspond à la position d’un astre quelconque, on parlera de
périastre et d’apoastre. Si l’astre est la Terre on retient les termes périgée et
apogée ; et pour le Soleil, périhélie et aphélie.
47
Guide de localisation des astres
Cas de l’hyperbole
p
On part de ρ0 = a − c =
1+ e
c
et on pose = e
a
p 1+ e a 2 − c2
On vérifie alors que : k = = (a − c) =
e e c
Introduisant b tel que : b2 = c2 − a2
b2 b2
on obtient k = − , p = −
c a
Afin que les paramètres k et p soient positifs, il convient que a et c soient
négatifs.
On retient finalement les deux jeux de formules :
p ep p
a = c = b =
1− e 2
1− e 2
e −1
2
b2 b2 a2
p = − k = − d =
a c c
48
2. Rappels sur les coniques
Cas où e = 1
L’origine O est prise au milieu du segment FD.
L’équation polaire de la parabole est
k
ρ=
1 + cos θ
On déduit ρcos θ = k − ρ
Le système : | x = ρcos θ − k/2
| y = ρsin θ
s’écrit : | (x + k/2)2 + y2 = ρ2
| x = k/2 − ρ
et compte tenu de l’expression de ρ dans la deuxième égalité, l’équation carté-
sienne de la parabole est :
y2 = −2 kx
où k représente la distance algébrique FD .
2 Définition bifocale
49
Guide de localisation des astres
M H H’
F A D O D’ A’ F’ x
e = 0,5
e = 0,2
e→0
•
F D
F
(D)
50
2. Rappels sur les coniques
3.1 Cas e → 0
– La conique, ellipse, se réduit normalement au point F ;
MF
cela se déduit de la définition = e
MH
– Toutefois, considérant l’équation polaire, on voit qu’il existe en fait deux
possibilités dépendant de la distance foyer directrice :
• k = constante,
Alors p → 0 et ρ → 0 : on obtient le point F ; c et a sont nuls.
• k → ∞,
Alors p = ek → Ctte = Rp ; et comme ρ → p on obtient un cercle de rayon Rp.
On note que c = 0,
a = b = Rp,
d = ∞.
Ainsi l’ellipse devient un cercle d’équation :
ρ = Rp ou x 2 + y 2 = R p2
3.2 Cas e → ∞
On déduit :
p p
ρ ≈ soit ρ cos θ ≈ ≈ k
e cos θ e
La conique, branche d’hyperbole, devient la droite d’équation :
ρ cos θ = k ou x = k
Cette droite est la directrice (D).
• À la figure 2.5, la distance foyer directrice était une constante. On peut
également tracer les courbes fonction du paramètre excentricité en considé-
rant que la distance ρ0 foyer-extremum est une donnée, ρ0 = Rp, comme dans
le cas de la mise en orbite d’un satellite d’une planète P.
La représentation sera fournie à la figure 7.1.
Dans cette hypothèse,
ek 1+ e
Rp = ; et k = R p
1+ e e
Pour e → 0 k → ∞ la courbe tend vers le cercle ;
e = 1 k = 2Rp la courbe est la parabole ;
e → ∞ k → Rp la courbe tend vers une droite.
51
Guide de localisation des astres
exercices
E2-1 Géométrie de l’ellipse en fonction du demi-grand axe
et de l’excentricité
1. Reporter sur une figure les différentes dimensions de l’ellipse en fonction
du demi-grand axe a et de l’excentricité e :
On fera ainsi apparaître les expressions de OB, OF, OD, FA', FA, FD.
Positionner le paramètre p.
2. Calculer pour la planète Mercure :
– OD, OF, OP, OB, ainsi que FD, FA, FP, FB, FP', FA' ;
– les angles (FA, FB) et (FA, FP').
Données : a = 0,387 098 ua ; e = 0,205 632.
3. Calculer pour la Terre la différence entre les rayons vecteurs à l’aphélie et
au périhélie ; l’exprimer en km et commenter.
Calculer également le demi-petit axe, ainsi que la différence en km avec le
demi-grand axe.
Données : a = 1,000 001 ua ; e = 0,016 710.
solutions
1.
B
P’ P
a(1-e )1/2
a(1-e )
• •
A’ F’ O F A D
B’
a ae a(1-e )/e
exercices
a(1+e) a(1-e)
a/e
52
2. Rappels sur les coniques
2. Planète Mercure
– Calcul de OD, OF, OP, OB :
OD = d = a/e → OD = 1,882 479 ua
OF = c = ae → OF = 0,079 600 ua
OA = a (donnée) → OA = 0,387 098 ua
OP = (c2 + p2)1/2 = a(e4 – e2 + 1)1/2 → OP = 0,379 179 ua
OB = b = a(1−e2)1/2 → OB = 0,378 825 ua
– Calcul de FD, FA, FP, FB, FP', FA' :
FD = k = a(1−e2)/e → FD = 1,802 880 ua
FA = ρ0 = a(1−e) → FA = 0,307 498 ua
FP = p = a(1−e2) → FP = 0,370 730 ua
FB = ρθb = a → FB = 0,387 098 ua
FP' = ρθp’ = (4c2+p2)1/2 = a(1+e2) → FP' = 0,403 466 ua
FA' = ρπ = a(1+e) → FA' = 0,466 698 ua
– Calcul de (FA, FB) et (FA, FP') :
e
(FA, FB) = 90° + arctg(c/b) = 90° + arctg → (FA, FB) = 101°,866
1 − e2
2e
(FA, FP') = 90° + arctg(2c/p) = 90° + arctg → (FA, FP') = 113°,240
1 − e2
3. Planète Terre
FA' – FA = ρπ − ρ0 = 2ae (c’est la distance focale) → FA' – FA = 0,033 420 ua ;
or, voir chapitre précédent, 1 ua = 149,597 871 106 km, d’où
FA' – FA = 4 999 566 km
La différence des rayons vecteurs à l’aphélie et au périhélie, c’est-à-dire la
distance focale, est de pratiquement cinq millions de kilomètres, soit 13 fois
la distance Terre-Lune, ou encore 125 périmètres terrestres.
OB = b = a(1 – e2)1/2 → OB = 0,999 861 ua
OA – OB = a – b → OA – OB = 0,000 140 ua
OA – OB = 20 885 km
Ainsi la différence entre le grand axe et le petit axe de notre planète est prati-
exercices
Remarque
Les rayons au péricentre et à l’apocentre sont parfois notés respectivement
q et Q :
q = FA = ρ0 = a(1−e)
Q = FA' = ρπ = a(1+e)
53
Guide de localisation des astres
solutions
Le paramètre est désigné par p ; l’excentricité par e.
• Cercle (e = 0)
Y y Y
A’ O(F) A
(-p,0) (p,0) X, x
Y y Y
A’ F’ O F A
(-a,0) (a,0) X, x
exercices
54
2. Rappels sur les coniques
Y=y Y=y
F A A’ F’
X=x X=x
Sommet A :
On a établi l’équation de la parabole, y2 = − 2kx
Or p = ek = k
On note donc Y2 = −2pX
Sommet A' :
Y2 = 2pX
Introduisant ε tel que défini précédemment, l’équation s’écrit :
Y2 = −2εpX
• Hyperbole (e > 1)
Y y Y
F A O A’ F’ x, X
exercices
55
Guide de localisation des astres
Remarque
Avant de clore ce chapitre, précisons que les termes conique, ellipse, pa-
rabole et hyperbole trouvent leur origine en particulier dans les travaux
de deux mathématiciens célèbres de la Grèce antique, à savoir Menaech-
mus (environ 380-320 avant J.-C.) et Appolonius de Perge (environ 262-
190 avant J.-C.).
C’est Menaechmus qui le premier définit les coniques comme intersection
d’un cône et d’un plan ; peut-être même aurait-il inventé les mots parabole
et hyperbole pour désigner les courbes obtenues respectivement lorsque le
plan est parallèle à une génératrice du cône ou lorsqu’il coupe les deux
nappes de celui-ci.
Environ un siècle plus tard, Appolonius de Perge écrivit un important
ouvrage en huit tomes qu’il intitula Coniques, toujours par référence à
l’intersection d’un cône avec un plan ; les termes ellipse, parabole, hyper-
bole y étaient introduits. Ce livre allait bien au-delà de tout ce qui avait été
étudié au paravant sur le sujet, y compris par le célèbre Euclide (environ
325-265 avant J.-C.).
exercices
56
Chapitre 3
Mouvements
à accélération centrale
Vz
z
Vθ
P
Vρ
z Y
K
J
F
y
θ
ρ
x
P
X
• Coordonnées de P
ρ est le rayon polaire, θ l’angle polaire et z la cote ; il vient immédiatement :
LLI I I
FP = ρ I + zk
• VitesseLLIde P I I I
LI dFP I dI I dI dI I
V = = ρ I + ρ + zk
; or = θ = θ J
dt dt dt dθ
on obtient par conséquent :
LI I I I
V = ρ I + ρθ J + zk
Les coordonnées cylindriques du vecteur vitesse sont donc :
Vρ = ρ
Vθ = ρθ
Vz = z
58
3. Mouvements à accélération centrale
LI
Pour un mouvement plan, les composantes de V sont la vitesse radiale Vρ et
la vitesse transversale Vθ.
• Accélération de P
LI
LI dV LI I I I
Γ = − ρθ 2 ) I + (2ρθ
, soit : Γ = (ρ + ρθ
) J +
zk
dt
Les coordonnées du vecteur accélération sont par conséquent :
Γρ = (ρ − ρθ 2 )
Γθ = (2ρθ + ρθ )
Γz =
z
V
A
P
Γ
F
59
Guide de localisation des astres
dθ dS V
F
ρ P
Γ θ
P0
Figure 3.3 Coordonnées polaires de P.
Le
LI vecteur
I vitesse
I a pour expression :
V = ρ I + ρθ J
Sa norme au carré est :
I 2
V = ρ 2 + ρ2θ 2 , on notera simplement V .
2
60
3. Mouvements à accélération centrale
⎡1 ⎛ d ⎛ 1⎞ ⎞ ⎤
2
C2 ⎡1 d2 ⎛ 1 ⎞ ⎤
V = C ⎢ 2 + ⎜ ⎜ ⎟⎟ ⎥
2 2
et Γ = − ⎢ + 2 ⎜ ⎟⎥
⎢⎣ ρ ⎝ dθ ⎝ ρ ⎠ ⎠ ⎥⎦ ρ2 ⎣ ρ dθ ⎝ ρ ⎠ ⎦
I μ I
Application dans le cas où Γ = − 2 I
ρ
La signification de la constante μ sera donnée plus tard, au chapitre 6.
• La formule de Binet relative à l’accélération conduit à :
μ C2 ⎡ 1 d2 ⎛ 1 ⎞ ⎤
= 2 ⎢ + ⎥
ρ 2
ρ ⎣ ρ dθ2 ⎜⎝ ρ ⎟⎠ ⎦
d’où en isolant les constantes C et μ
1 d2 ⎛ 1 ⎞ μ
+ = 2
ρ dθ2 ⎜⎝ ρ ⎟⎠ C
La solution générale de l’équation sans second membre est :
1
= Bcos(θ − θ0 )
ρ
B et θ0 étant les deux constantes d’intégration.
La solution particulière de l’équation avec second membre est :
1 μ
= 2
ρ C
On déduit la solution générale de l’équation avec second membre :
1 μ
= 2 + Bcos(θ − θ0 )
ρ C
61
Guide de localisation des astres
62
3. Mouvements à accélération centrale
p2 ρ0 + ρ π p
On sait également que : ρ0ρπ = et =
1 − e2 2 1 − e2
Ainsi :
μ μ
V0 Vπ = =
ρ0 + ρ π a
2
Autre résultat intéressant :
V0 Vπ ρ 0ρ π
= = 1 − e2
V0 + Vπ ρ0 + ρ π
2 2
Enfin, les rapports entre ρπ et ρ0 d’une part, Vπ et V0 d’autre part s’expriment
simplement en fonction de e :
1+ e 1− e
ρ π = ρ0 et Vπ = V0
1− e 1+ e
ce qui conduit à :
ρπVπ = ρ0V0
Parabole (e = 1)
L’expression de V en fonction de e devient :
μ μ
V = 2 1 + cos θ = 2 1 + cos θ
C p
et Vθ→±π ⇒ 0
63
Guide de localisation des astres
Hyperbole (e > 1)
On a noté que −θM < θ < θM avec cos(±θM) = −1/e.
Les expressions des valeurs particulières sont :
μ μ μ
V0 = (1 + e ) = (1 + e ) = 1+ e
C p ρ0
μ 2 μ 2
Vθ → ± θM ⇒ e −1 = e −1
C p
Ellipse
μ μ V02 1 V02
Γ0 = − (1 + e ) 2
= − = − = −
p2 ρ20 ρ0 1 + e p
μ μ V2 1 V2
Γπ = − (1 − e )2 = − 2 = − π = − π
p 2
ρπ ρπ 1 − e p
Parabole
μ μ V02 V02
Γ 0 = −4 = − = − = −
p2 ρ20 2ρ0 p
Γθ→±π ⇒ 0
64
3. Mouvements à accélération centrale
Hyperbole
μ μ V02 1 V02
Γ0 = − (1 + e ) 2
= − = − = −
p2 ρ20 ρ0 1 + e p
Γθ→±θm ⇒ 0
V
J
α
I
M
N
ρ
ρ’ α
θ
F x
65
Guide de localisation des astres
On retiendra :
1 + e cos θ
tg α =
e sin θ
Donnons quelques valeurs particulières du couple (α,θ) :
– Pour θ = 0 → α = π/2
– Pour θ = π/2 → tg α = 1/e
– L’angle θ tel que θ + α = π, c’est-à-dire tel que le vecteur vitesse soit hori-
zontal, satisfait à :
1 + e cos θ
tg α = tg (π−θ) = − tg θ =
e sin θ
La relation est : cos θ = − e, pourvu que e ≤ 1,
1 − e2
cos α = cos (π − θ) = − cos θ = + e, soit encore tg α =
e
module : V = ρ2 + ρ '2 θ = NM θ
composante radiale : Vρ = ρ ' θ = FN θ
composante orthoradiale : Vθ = ρθ = FM θ
C
Si le mouvement obéit à la loi des aires, θ = 2 , et
ρ
1 ρ' 2 C 1
V = C 2 + 4 = ρ sin α (c’est la première formule de Binet)
ρ ρ
ρ' C cos α
Vρ = C 2 =
ρ ρ sin α
1
Vθ = C
ρ
66
3. Mouvements à accélération centrale
Si en plus l’accélération est en 1/ρ2, alors ρ = p/(1 + e cos θ), et les expressions
s’écrivent :
C
V = 1 + e 2 + 2e cos θ
p
C
Vρ = e sin θ
p
C
Vθ = (1 + e cos θ)
p
Remarques : LLLI LI
1– Le moment du vecteur vitesse par rapport à F étant FM ∧ V , on en
tire que
LLLI LI
FM ∧ V = ρV sin α = C (voir figure 3.4)
d’où : C = ρV sinα = ρVθ
C V 1 + e cos θ
et : sin α = = θ =
ρV V 1 + e 2 + 2e cos θ
Vρ e sin θ
cos α = =
V 1 + e 2 + 2e cos θ
Pour θ = 0, ou θ = π, α = π/2, ce qui donne :
C = ρ0V0 = ρπVπ
L’intérêt de la relation C = ρV sin α apparaîtra par exemple aux exercices
E5-1, E7-4 et E8-1.
2– L’expression de l’accélération centrale en fonction de la première
constante dynamique C et du paramètre p est :
C2
Γ = − 3
(1 + e cos θ)2
p
67
Guide de localisation des astres
exercices
E3-1 Vitesses en certains points de l’ellipse ; relations déduites
On considère l’ellipse de l’exercice E2-1 sur laquelle figurent certains points
particuliers.
1. Calculer le module de la vitesse et ses composantes en A, P, B, P', et A'.
2. Calculer les produits VAVA' , VPVP' , et VB2 .
3. On donne la relation suivante entre anomalie vraie θ et excentrique ϕ :
cos ϕ − e
cos θ = (voir chapitre 5)
1 − e cos ϕ
En déduire l’expression du module de la vitesse en fonction de ϕ.
Calculer le produit Vϕ1Vϕ2 pour deux positions symétriques par rapport au
petit axe de l’ellipse ; conclure.
4. Calculer les produits ρVθ en A, P, B, P', et A' ; remarques ?
solutions
1. Module V de la vitesse, composantes Vρ et Vθ
μ
Afin de simplifier l’écriture, on pose : m =
a
C μ C m
On a établi que = ; or p = a(1− e2) , d’où : =
p p 1 − e2
p
On en déduit les expressions du module et des composantes de la vitesse ;
1 + e 2 + 2e cos θ e sin θ 1 + e cos θ
V = m Vρ = m Vθ = m
1 − e2 1− e 2
1 − e2
– en A, θ = 0,
1+ e
VA = V0 = m VρA = 0 VθA = VA
1− e
– en P, θ = π/2,
1 + e2 e 1
VP = m VρP = m VθP = m
1 − e2 1− e 2
1 − e2
– en B, cos(π−θ) = FO/FB = e, d’où : cos θ = − e et sin θ = (1− e2)1/2, il vient
exercices
donc :
μ
VB = m = VρB = me VθB = m 1 − e 2
a
– en P', cos(π−θ) = FF'/FP' = 2e/(1 + e2), d’où : cos θ = − 2e/(1 + e2)
et sin θ = (1−e2)/(1 + e2),
et par conséquent :
68
3. Mouvements à accélération centrale
1 − e2 1 − e2 1 − e2
VP ' = m VρP ' = me VθP ' = m
1 + e2 1 + e2 1 + e2
– en A', θ = π,
1− e
VA ' = Vπ = m VρA' = 0 VθA' = VA'
1+ e
2. Produits VAVA', VPVP', et VB2
On a immédiatement :
VAVA' (= V0Vπ) = VPVP' = VB2 = m2 = μ/a
On déduit :
ρAVθA (= ρ0V0) = ρPVθP = ρBVθB = ρP'VθP' = ρA'VθA' (= ρπVπ) =
μa 1 − e 2 ≡ C
69
Guide de localisation des astres
En effet :
a (1 − e 2 ) μ 1 + e cos θ
C = ρVθ = = μa 1 − e 2
1 + e cos θ a 1 − e2
D’autre part :
1 ua = 149,597 870 61 10 6 km (voir chapitre 1) ;
la valeur de la constante héliocentrique de la gravitation est :
μS = 1,327 124 4 10 20 m3. s –2 (voir chapitre 6) ;
1 jour = 86 400 s .
solutions
• Formules utilisées :
μ S 1 + e 2 + 2e cos θ
Module de la vitesse V =
a 1 − e2
μS 1 + e
Vitesse au périhélie VA = V0 =
a 1− e
exercices
μS 1 − e
Vitesse à l’aphélie VA ' = Vπ =
a 1+ e
μS
Vitesse moyenne VB = V0 Vπ =
a
Les vitesses sont exprimées en km.s -1 , μS = 13,271 244 1010 km3. s –2.
70
3. Mouvements à accélération centrale
• Tableau :
Remarques
1– Le tableau précédent donne les vitesses des planètes en certains points de
leur orbite. La vraie vitesse moyenne sur la trajectoire elliptique passe par
la connaissance du périmètre et sera donnée au chapitre 8 (paragraphe 2.1).
Ci-dessous, nous avons calculé les vitesses d’entraînement définies comme
suit :
2πR
VE =
J
où R est le rayon équatorial de la planète et J sa période de rotation, mesurée
avec certaines réserves pour les planètes joviennes qui n’ont pas de sol.
71
Guide de localisation des astres
72
Chapitre 4
Mouvement dans
le repère de Frenet
τ J
V
α I
n
ρ
j
R
ϕ
θ
F i
W
ρ’
N α
74
4. Mouvement dans le repère de Frenet
I I
I dτ dτ ds
n = =
dϕ ds dϕ
ds
Posant R = il vient :
dϕ
I LLLI
I dτ I d 2 FM
n = R soit : n = R
ds ds 2
2 Rayon de courbure
2.1 Définition
ds
La grandeur R =
dϕ
introduite ci-dessus est par définition le rayon de courbure en M.
M’
M
dθ ϕ + dϕ
ϕ
θ
Ω
75
Guide de localisation des astres
2.2 Expressions
Calcul de R(ρ, ρ', ρ'')
ds ds dθ dθ
R = = = ρ 2 + ρ '2
dϕ dθ dϕ dϕ
dϕ dα
or ϕ = θ+α → = 1+
dθ dθ
ρ dα ρ '2 − ρρ "
et tgα = → (1 + tg2α ) =
ρ' dθ ρ '2
(on a dérivé la tangente par rapport à θ)
En conséquence,
dα ρ '2 − ρρ " 1 ρ '2 − ρρ "
= = 2
dθ ρ' 2
1 + ρ / ρ'
2 2
ρ + ρ '2
et
dϕ dα ρ2 + 2ρ '2 − ρρ "
= 1+ =
dθ dθ ρ 2 + ρ '2
Le rayon de courbure a donc pour expression :
(ρ2 + ρ '2 ) 3 / 2
R =
ρ2 + 2ρ '2 − ρρ "
76
4. Mouvement dans le repère de Frenet
On exprime le numérateur :
p2 p2 e 2 sin 2 θ
ρ 2 + ρ '2 = +
(1 + e cos θ)2 (1 + e cos θ)2 (1 + e cos θ)2
p2
= (1 + e 2 + 2e cos θ)
(1 + e cos θ)4
p3
(ρ2 + ρ '2 )3 / 2 = (1 + e 2 + 2e cos θ)3 / 2
(1 + e cos θ)6
On obtient ainsi l’expression du rayon de courbure de la conique en fonction
de l’angle θ :
(1 + e 2 + 2e cos θ)3 / 2
R = p
(1 + e cos θ)3
Remarques
C
1– On a établi au chapitre précédant que V = 1 + e 2 + 2e cos θ
p
On en déduit :
ρ3 p3 3 avec C2
R = V p =
p 2 C3 μ
On retient :
p 3 3 p 1 3 3
R = ρV = = ρV
C 3
sin α
3
μC
77
Guide de localisation des astres
Par conséquent :
d 2s ρρ '+ ρ ' ρ '' 2
Γτ = = θ + ρ2 + ρ '2
θ
dt 2 ρ 2 + ρ '2
et donc :
d 2s ρρ '+ ρ ' ρ '' C2 ρ2 + ρ '2 −2C2ρ '
= +
dt 2 ρ 2 + ρ '2 ρ ρ 2 + ρ '2 ρ
4 5
C2
= (ρ2ρ '+ ρρ ' ρ ''− 2ρ2ρ '− 2ρ '3 )
ρ 5
ρ + ρ'
2 2
78
4. Mouvement dans le repère de Frenet
79
Guide de localisation des astres
soit :
(ρ2 + 2ρ '2 − ρρ '')
Γ = C2
ρ5
Cas où la trajectoire est une conique
2
p2 (1 + e cos θ)5 C2 ⎛ 1 + e cos θ⎞⎞
Γ = C2 =
(1 + e cos θ)3 p5 p ⎜⎝ p ⎟⎠
C2 μ
Or, = μ , d’où : Γ = 2
p ρ
Cela veut dire que si le mouvement obéit à la loi des aires, et que la trajectoire
est une conique, alors l’accélération est en 1/ρ2.
Expressions vectorielles
Le mouvement satisfaisant toujours à la loi des aires, on peut aussi écrire,
LI I I
compte tenu de Γ = Γ τ τ + Γ n n :
LI C2 (ρ2 + 2ρ '2 − ρρ '') I I
Γ = 5 ( −ρ ' τ + ρ n )
ρ ρ + ρ'
2 2
80
4. Mouvement dans le repère de Frenet
ρ
Or, on a vu que = tgα
ρ'
ce qui permet de déduire :
LI C2 (ρ2 + 2ρ '2 − ρρ '') 1 I I
Γ = 4 (− cos ατ + siin α n)
ρ ρ + ρ'
2 2 sin α
avec
1 1 + tg2α ρ '2 + ρ 2
= =
sin α tgα ρ
D’où finalement :
LI C2 I I
Γ = − 5 (ρ2 + 2ρ '2 − ρρ '')(cos ατ − sin α n)
ρ
LI I I
Écrivons maintenant Γ dans le repère {M, I , J }.
I I I I
La matrice de passage de {M, I , J } à {M, τ,n } est telle que
I I
τ cos α sin α I
I = I
n − sin α cos α J
On obtient donc :
LI C2 I
Γ = − 5 (ρ2 + 2ρ '2 − ρρ '') I
ρ
On démontre ainsi le théorème réciproque suivant :
Si un mouvement satisfait à la loi des aires, il est à accélération centrale ;
et, comme on vient de le voir :
Si en plus la trajectoire est une conique, alors le module de l’accélération est
proportionnel à 1/ρ2, le coefficient étant μ = C2/p.
On retrouve également la deuxième formule de Binet :
C2 2 C2 ⎛ 1 2ρ '2 − ρρ '' ⎞
Γ = − (ρ + 2ρ ' 2
− ρρ '') = − +
ρ5 ρ2 ⎝⎜ ρ ρ3 ⎟⎠
C2 ⎛ 1 d2 ⎛ 1 ⎞ ⎞
= − +
ρ2 ⎜⎝ ρ dθ2 ⎜⎝ ρ ⎟⎠ ⎟⎠
81
Guide de localisation des astres
exercices
LI KI
E4-1 Accélération centrale de la forme Γ = − 3 I
ρ
1. Déterminer la trajectoire.
2. Que peut-on dire de l’angle du vecteur vitesse avec le rayon vecteur ?
Donner les expressions de la vitesse et de ses composantes.
3. Calculer le rayon de courbure.
4. Établir les composantes de l’accélération dans le repère de Frenet.
solutions
1. Trajectoire
La formule de Binet relative à l’accélération permet d’écrire :
C2 ⎛ 1 d2 ⎛ 1 ⎞ ⎞ K
+ = 3 ; soit en isolant les constantes :
ρ2 ⎜⎝ ρ dθ2 ⎜⎝ ρ ⎟⎠ ⎟⎠ ρ
⎛1 d2 ⎛ 1 ⎞ ⎞ K d2 ⎛ 1 ⎞ K
ρ⎜ + 2 ⎜ ⎟⎟
= 2 , puis,
ρ 2 ⎜ ⎟
= 2 −1
⎝ ρ dθ ⎝ ρ ⎠ ⎠ C dθ ⎝ ρ ⎠ C
K 1
On pose K' = − 1 , et u = ;
C2 ρ
d’où l’équation différentielle du second ordre :
u'' – K'u = 0 ,
dont la solution est du type u(θ) = Ae α ( θ − θ0 )
avec par identification α = ± K ' ; K > C2
Moyennant un choix convenable de l’axe polaire, la trajectoire a pour équation :
K − C2
± θ
ρ = ρ0 e C
ρ 1 ε
tg α = = = → α = Ctte
ρ' εk k
82
4. Mouvement dans le repère de Frenet
On dit que « les rayons vecteurs sont coupés sous un angle constant ».
ε = + 1 → 0 < α < π/2 ε = − 1 → π/2 < α < π
• Module et composantes :
LI I I I
V = V τ = Vρ I + Vθ J
C
Le mouvement obéit à la loi des aires, donc θ = 2 ; ainsi :
ρ
2
1 ⎛ ρ' ⎞ 1 k2 C 1
V = C + ⎜ 2⎟ = C + = 1 + K' → V = K
ρ 2
⎝ρ ⎠ ρ 2
ρ 2
ρ ρ
C ρ' 1 1 1
Vρ = = εkC = ε K ' C → Vρ = ε K − C2
ρ ρ ρ ρ ρ
C
Vθ =
ρ
3. Rayon de courbure
K 2
ρ2 = ρ20 e 2 ε k θ ρ '2 = k 2 ρ 2 → ρ2 + ρ '2 = ρ2 (1 + k 2 ) = ρ
C2
K 2
ρ2 + 2ρ '2 − ρρ '' = ρ2 (1 + 2 k 2 − k 2 ) = ρ2 (1 + k 2 ) = ρ
C2
K
On déduit : R = ρ
C
4. Composantes de l’accélération dans le repère de Frenet
LI I I
Γ = Γττ + Γnn
• Accélération tangentielle
d 2s 2 1 ρ ' ρ + 2ρ ' − ρρ ''
2 2
C2 K / C2
Γτ = = − C = − ε k ρ
dt 2 ρ4 ρ ρ 2 + ρ '2 ρ4 K / C
1
Γ τ = − ε K K − C2
ρ3
• Accélération normale
V2 1 ρ2 + 2ρ '2 − ρρ '' C2 K
Γn = = C2 4 = ρ
exercices
R ρ ρ 2 + ρ '2 ρ4 C
C K
Γn =
ρ3
LI
On vérifie que Γ = K
ρ3
83
Guide de localisation des astres
LI C I
2
E4-2 Cas particulier : accélération centrale de la forme Γ = − 3 I
ρ
1. Déterminer la trajectoire, la tracer ; on précisera les asymptotes ainsi que
la sous-normale et la sous-tangente.
2. Donner les expressions de la vitesse et de ses composantes.
3. Calculer le rayon de courbure et établir les composantes de l’accélération
dans le repère de Frenet.
solutions
1. Trajectoire
Dans ce cas particulier où K = C2, la formule de Binet relative à l’accélération
conduit à :
d2 ⎛ 1 ⎞
= 0
dθ2 ⎜⎝ ρ ⎟⎠
1 θ A
On retient la solution = , soit ρ =
ρ A θ
Il s’agit d’une spirale hyperbolique (figure 4.3) ;
quand θ → ∞ ρ → 0 (point asymptotique) ;
quand θ → 0 ρ → ∞ et ρ sin θ → A (ordonnée de l’asymptote horizontale).
X
y
I
A
α
V M
Y
ρ
T A
θ
F x
α
ρ’
N
exercices
• Sous-normale n, sous-tangente t :
M est un point de la spirale, le rayon FM est porté par l’axe FX. La perpendi-
culaire à FX est FY ; la tangente en M coupe cet axe en T, la perpendiculaire
en M le coupe en N.
84
4. Mouvement dans le repère de Frenet
– Calcul de n = FN
LLLLI I I I 1I I I
MN = − ρ I + nJ , et un vecteur tangent est T = V = ρ ' I + ρJ
θ
On déduit :
I LLLLI A ρ ρ2
T MN = 0 → −ρρ ' + nρ = 0 → n = ρ' , ρ ' = − = − = −
θ2 θ A
– Calcul de t = FT
Considérant les triangles FTM et FMN, il vient :
ρ2
t = − = A
ρ'
La sous-tangente à la spirale hyperbolique est égale à la constante d’intégration.
2. Vitesse
1 ρ '2 1 1 C ρ2
V = C + = C + = 1 +
ρ2 ρ4 ρ2 A2 ρ A2
ρ' C
Vρ = C = − → La vitesse radiale est constante.
ρ 2
A
C
Vθ =
ρ
3. Rayon de courbure R
ρ2 + ρ '2 = ρ2 (1 + ρ2 / A 2 ) et ρ2 + 2ρ '2 − ρρ '' = ρ2
d’où :
3/ 2
⎛ ρ2 ⎞
R = ρ ⎜1 + 2 ⎟
⎝ A ⎠
• Accélération
C2 1 1 1 1
Γτ = , Γ n = C2
A ρ2 1 + ρ2 / A 2 ρ3 1 + ρ 2 / A
Remarque
Nous avons noté au chapitre introductif que de nombreuses galaxies possé-
exercices
85
This page intentionally left blank
Chapitre 5
Anomalie et orbite
1 Anomalies
Définitions
Kepler en est à l’origine. On considère une planète P d’orbite une ellipse dont
le foyer F est par exemple le Soleil.
On se réfère à la figure 5.1 ci-dessous.
ϕ
θ
F’ O H F A
88
5. Anomalie et orbite
Équation de Kepler
Sachant que l’aire de l’ellipse est πab (voir chapitre 8, paragraphe 2.2), la loi
des aires permet d’écrire :
Aire ( AFP) πab
=
t T
Or,
b b
Aire(AFP) = Aire(AFQ) = [Aire(AOQ) – Aire (FOQ)]
a a
b1 2
= (a ϕ − ca sin ϕ )
a2
ab
= (ϕ − e sin ϕ )
2
On obtient l’équation de Kepler :
ϕ − e sinϕ = ωmt = M
conduisant à :
ϕ − e sin ϕ
t = T
2π
89
Guide de localisation des astres
d’où :
cos ϕ − e 1 − e 2 sin ϕ
cos θ = sin θ =
1 − e cos ϕ 1 − e cos ϕ
e + cos θ 1 − e 2 sin θ
cos ϕ = sin ϕ =
1 + e cos θ 1 + e cos θ
Donnons également des relations entre les angles moitiés.
Partant des égalités :
θ (1 − e )(1 + cos ϕ ) 1− e ϕ
2 cos 2 = 1 + cos θ = = 2 cos 2
2 1 − e cos ϕ 1 − e cos ϕ 2
θ (1 + e )(1 − cos ϕ ) 1+ e ϕ
2 sin 2 = 1 − cos θ = = 2 sin 2
2 1 − e cos ϕ 1 − e cos ϕ 2
on écrit :
ϕ ϕ
θ 1 − e cos θ 1 + e sin
cos = 2 sin = 2
2 1 − e cos ϕ 2 1 − e cos ϕ
Les tangentes des demi-anomalies sont ainsi liées par un coefficient de propor-
tionnalité, simple fonction de l’excentricité :
θ 1+ e ϕ ϕ 1− e θ
tg = tg tg = tg
2 1− e 2 2 1+ e 2
En conclusion, les formules permettent de traiter deux types de problèmes :
soit, connaissant l’anomalie vraie, en déduire le temps écoulé depuis le passage
au péricentre ; soit, connaissant le temps écoulé, en déduire l’anomalie vraie.
La grandeur intermédiaire est dans les deux cas l’anomalie excentrique.
90
5. Anomalie et orbite
C μ e sin ϕ
Vρ = e sin θ → Vρ =
p a 1 − e cos ϕ
C μ 1 − e2
Vθ = (1 + e cos θ) → Vθ =
p a 1 − e cos ϕ
1 + e cos θ 1 − e2
tgα = → tgα =
e sin θ e sin ϕ
• l’expression de l’accélération centrale en fonction de ϕ :
μ μ 1
Γ = − (1 + e cos θ)2 → Γ = −
p2 a (1 − e cos ϕ )2
2
I
Vecteur tangent unitaire τ et équation de la normale en P
LLLI LLLI
I dOP dOP dϕ
τ = =
ds dϕ ds
I I
α étant l’angle ( i , τ ), les coordonnées sont
x' y'
cos α = sin α =
x' + y'
2 2
x ' + y '2
2
91
Guide de localisation des astres
On retient :
− a sin ϕ
I a sin ϕ + b2 cos 2 ϕ
2 2
τ
b cos ϕ
a sin ϕ + b2 cos 2 ϕ
2 2
LLLI I
PM τ = 0 → ( x − x p ) cos α + ( y − y p )sin α = 0
− a sin ϕ
→ ( x − a cos ϕ )
a sin ϕ + b2 cos 2 ϕ
2 2
b cos ϕ
+ ( y − b sin ϕ ) = 0
a sin ϕ + b2 cos 2 ϕ
2 2
I
Vecteur normal unitaire n et équation de la tangente en P
n
− a sin ϕ
sin(α + π / 2) = cos α =
a sin ϕ + b2 cos 2 ϕ
2 2
92
5. Anomalie et orbite
− b cos ϕ
→ ( x − a cos ϕ )
a sin 2 ϕ + b2 cos 2 ϕ
2
−a sin ϕ
+ ( y − b sin ϕ ) = 0
a sin ϕ + b2 cos 2 ϕ
2 2
F’ F
93
Guide de localisation des astres
Remarques
1– La construction géométrique montre que l’ellipse est la courbe enve-
loppe des médiatrices de F'G, où F' est à l’intérieur du cercle, obtenue
lorsque le point G a décrit le cercle. Si F' était à l’extérieur du cercle, la
courbe enveloppe serait l’hyperbole.
2– Considérons maintenant la normale en P. La courbe enveloppée par
celle-ci lorsque G décrit le cercle est la développée de l’ellipse ; c’est une
astroïde. On dit aussi que l’ellipse est la développante !
Il y a ainsi deux définitions pour la développée d’une courbe ; à savoir
lieu des centres de courbure (voir paragraphe 2.2), ou bien enveloppe des
normales.
Rayon de courbure
I I
Rappelons que α est ici l’angle ( i , τ ).
ds
Par définition le rayon de courbure est R =
dα
On introduit la différentielle dϕ de l’anomalie excentrique ; ainsi :
ds ds dα
=
dϕ dα dϕ
et en dérivant tg α :
y' dα x ' y ''− y ' x '' dα x ' y ''− y ' x ''
tg α = → (1 + tg2 α) = → =
x' dϕ x' 2
dϕ x '2 + y '2
Or
ds dϕ ds
R = avec = x '2 + y '2
dϕ dα dϕ
Le rayon de courbure a donc pour expression :
( x '2 + y '2 ) 3 / 2
R =
x ' y ''− y ' x ''
ce qui pour l’ellipse conduit à :
(a 2 sin 2 ϕ + b2 cos 2 ϕ ) 3 / 2
R =
ab
soit encore, en fonction uniquement du cosinus de l’anomalie excentrique :
a2
R = (1 − e 2 cos 2 ϕ )3 / 2
b
Au paragraphe suivant nous allons déterminer le lieu du centre de courbure Ω
en fonction du paramètre ϕ.
94
5. Anomalie et orbite
d’où
c2
xd = + cos3 ϕ
LLLI a
OΩ(ϕ ) =
c2 3
yd = − sin ϕ
b
C’est l’équation paramétrique d’une astroïde, représentée à la figure 5.3.
Donnons deux valeurs particulières :
LLLI c2
•ϕ=0
OΩ(0) = ( , 0) , le point représentatif est Ad ;
a
LLLI c 2
• ϕ = π/2 OΩ( π / 2) = (0, − ) , le point représentatif est Bd .
b
b2
Pour ϕ = 0, le rayon de courbure est : R (0) = OA − OA d =
a
a2
Pour ϕ = π/2, le rayon de courbure est : R ( π / 2) = Bd O + OB =
b
La moyenne géométrique des rayons de courbure en A et en B est égale à ab ,
c’est-à-dire ρm (voir chapitre 2, paragraphe 1.2).
Pour la Terre :
R(0) = 0,999 722 ua
R(π/2) = 1,000 141 ua
ρm = 0,999 931 ua
Les points Ad et Bd sont positionnés sur la figure à l’aide d’une construction
géométrique simple :
Q est le premier sommet du rectangle tangent à l’ellipse en A et B ; de ce
point, on trace la perpendiculaire au segment AB qui coupe l’axe des x en Ad
et l’axe des y en Bd .
95
Guide de localisation des astres
En effet,
Ad A AQ b2
– les triangles AdAQ et AQB étant semblables → = → Ad A =
AQ QB a
Bd B BQ a2
– les triangles BdBQ et BQA étant semblables → = → Bd B =
BQ QA b
La figure 5.3 qui suit précise toutes les notations utilisées :
B Q
P
O Ad F A
x
Bd
96
5. Anomalie et orbite
exercices
E5-1 Temps de l’intersection d’une comète avec l’écliptique
La trajectoire elliptique de la Comète de Halley coupe l’écliptique – plan de
l’orbite terrestre – au point Nd (appelé nœud descendant).
1. Quel est le temps mis par la comète pour atteindre ce point après son
passage au périhélie ? La date au périhélie étant le 9 février 1986 à 11 h
00 min 25 s, quelle est celle au nœud descendant ?
2. Calculer la vitesse au périhélie ainsi que le rayon et la vitesse au nœud
descendant.
3. Calculer la constante des aires et les produits rayon-vitesse au périhélie et
au nœud descendant ; conclure.
4. Déduire l’angle αNd du vecteur vitesse avec le rayon vecteur au point Nd ,
et calculer les composantes radiale et orthoradiale de la vitesse.
5. Calculer la date, le rayon vecteur et la vitesse de la comète au nœud ascen-
dant Na avant son passage au périhélie ; on précise que le Soleil appartient à
la ligne des nœuds.
Données numériques :
– angle des rayons vecteurs au périhélie et à l’intersection de l’écliptique
θNd = 68°,153 48
– excentricité e = 0,967 277 8
– période sidérale TH = 76 années
– rayon au périhélie ρ0 = 87,829 116 10 6 km
– constante héliocentrique de la gravitation μS = 13,271 244 10 10 km3s -2.
solutions
1. Temps tNd
On connaît l’anomalie vraie θNd comptée à partir de la position au périhélie ;
on en déduit l’anomalie excentrique ϕNd par la relation (chapitre 5, para-
graphe 1.2) :
exercices
1 − e 2 sin θ Nd
sin ϕ Nd = → sin ϕNd = 0,173 168
1 + e cos θ Nd
ϕNd = 9°,972 046 = 0,174 045 rd
L’équation de Kepler conduit à :
T
t Nd = (ϕ Nd − e sin ϕ Nd ) H ; avec TH = 27 759 jours (l’année est de 365,25 jours)
2π
97
Guide de localisation des astres
• ρNd et VNd
p
ρ Nd =
1 + e cos θ Nd
avec
p = ρ0 (1 + e ) p = 172,784 270 106 km = 1,154 992 ua
l’application numérique donne ρNd = 127,052 401 106 km = 0,849 293 ua
μS
VNd = 1 + e 2 + 2e cos θ Nd
p
on calcule : VNd = 45,163 km/s
98
5. Anomalie et orbite
• VρNd et VθNd
VρNd = VNd cos αNd → VρNd = 24,882 km/s
VθNd = VNd sin αNd = CH / ρNd VθNd = 37,690 km/s
Remarques
1 1 2
1– On pourra vérifier que + = , et plus généralement
ρ Nd ρ Na p
1 1 2
+ = .
ρ(θ) ρ(θ + π ) p
exercices
99
Guide de localisation des astres
3–
Trajectoire de la comète anomalie
rayon
entre les nœuds Na et Nd
ua °
2,5 100
69,72
rayon anomalie
2,0 1,843 50
0
1,5 0
0,587
0,0 -150
-110 -100 -90 -80 -70 -60 -50 -40 -30 -20 -10 0 10 20 30
100
5. Anomalie et orbite
Halley est une comète de courte période. On distingue en effet deux caté-
gories principales de comètes : celles à longue période, dont l’aphélie se
situe souvent à une distance de 50 000 ua correspondant au nuage de Oort
(voir chapitre 1), et celles dont la période est inférieure à 200 années.
Les comètes de cette deuxième catégorie pourraient être des comètes à
longue période qui auraient subi les forces gravitationnelles des grosses
planètes ; l’autre hypothèse est qu’elles proviendraient de la ceinture de
Kuiper située au voisinage de l’écliptique et qui s’étend de l’orbite de
Neptune jusqu’à environ 100 ua.
Les comètes de courte période sont généralement regroupées en deux
familles : celles de la famille de Jupiter, qui ne vont guère au-delà de cette
planète et ont des orbites proches de l’écliptique, et celles de la famille
de Halley.
Depuis 1995 la nomenclature des comètes, recommandée pour les dési-
gner, est établie comme suit :
– la première lettre fait référence à la période (P pour courte et C pour
longue ou non périodique), et si la comète a disparu on utilise la lettre D ;
– l’année de découverte vient ensuite ;
– elle est suivie d’une lettre représentant la quinzaine du mois de cette
découverte (lettre A pour la première quinzaine de janvier, F pour la
deuxième quinzaine de mars), sachant que le I n’est pas utilisé ;
– puis un nombre indique le rang de la découverte dans la quinzaine ;
– enfin, le nom du premier découvreur et éventuellement du deuxième
peut être rajouté.
On a répertorié à ce jour pas moins de 1000 comètes, dont près de 200
sont de courte période et numérotées par ordre chronologique de décou-
verte. Voici quelques exemples :
1P/Halley, la première de la liste !
exercices
101
Guide de localisation des astres
Équinoxe
de printemps
Aphélie
Équinoxe
d’automne
1. Calculer
a) la durée correspondant à θ0 ;
b) la durée des saisons.
exercices
102
5. Anomalie et orbite
solutions
La résolution s’effectue à partir de :
– l’équation de Kepler donnant le temps en fonction de l’anomalie excen-
trique
T
t = (ϕ − e sin ϕ )
2π
– la relation entre anomalie excentrique et anomalie vraie
1 − e 2 sin θ
sin ϕ =
1 + e cos θ
On s’aide d’un tableur, en utilisant la bonne détermination de la fonction
Arc sin (voir exercice E5-1) selon les valeurs θB = Arc cos(−e) = 90°,9574 et
θB' = 360° − θB = 269°,0426.
1.
a) Durée t0 entre solstice d’hiver et périhélie
θ = 360° − 12°,9375 θ = 347°,0625
d’où ϕ = 347°,275 124 → ta = 352,559 682 j
durée cherchée t0 = 365,2564 – ta → t0 = 12,696 718 j
b) Durée des saisons
• L’hiver ; il correspond au temps entre solstice d’hiver et équinoxe de printemps.
La durée du périhélie à l’équinoxe de printemps correspond à une anomalie vraie
θ = 90° − θ0 , soit θ = 77°,0625 ;
on déduit :
ϕ = 76°,131 145 th = 76,299 728 j
Il faut rajouter t0 pour avoir la durée DH de l’hiver : DH = 88,996 j
• Le printemps ; il correspond au temps entre équinoxe de printemps et sols-
tice d’été.
exercices
103
Guide de localisation des astres
104
5. Anomalie et orbite
solutions
1. Points et vecteurs
P = (a cos ϕ , b sin ϕ )
F = F1 = (c, 0), F' = F-1 = (– c, 0) ; on introduit ε = ±1 , ce qui permet de
regrouper : Fε = (εc, 0)
LLI
PF ε = (εc − a cos ϕ , − b sin ϕ )
LLI
PF ε = ε 2 c 2 + a 2 cos 2 ϕ − 2εac cos ϕ + b2 − b2 cos 2 ϕ
on déduit :
PFε . N b
cos(PFε , N) = =
PFε N (a + c cos ϕ )(a − c cos ϕ )
• Sinus
LLI LI
On calcule le produit vectoriel de PF ε et N :
LLI LI I
PF ε ∧ N = (− εac sin ϕ + a 2 sin ϕ cos ϕ − b2 sin ϕ cos ϕ ) k
I
= (− εac sin ϕ + ε 2 c 2 sin ϕ cos ϕ ) k
d’où :
LLI LI I
PF ε ∧ N = − εc sin ϕ (a − εc cos ϕ ) k
exercices
et
LLI LI
PF ε ∧ N = c sin ϕ (a − εc cos ϕ )
on déduit
LLI LI
LLI LI PF ε ∧ N c sin ϕ
sin(PF ε , N) = LLI LI =
PF ε N (a + c cos ϕ )(a − c cos ϕ )
105
Guide de localisation des astres
LLI LI
Calcul de la tangente de l’angle (PF ε , N)
• En fonction de l’anomalie excentrique ϕ, on obtient immédiatement
LLI LI c
tg (PF ε , N) = sin ϕ
b
On vérifie que pour ϕ = π/2 , la tangente est égale à c/b (voir figure 2.6.).
• En fonction de l’anomalie vraie θ
Se référant au chapitre 2, on a :
ep 1 − e2 p c e
c = ; b = → =
1 − e2 1 − e2 b 1 − e2
Au début du présent chapitre on a établi que :
1 − e 2 sin θ
sin ϕ =
1 + e cos θ
D’où l’expression en fonction de e et θ :
LLI LI e sin θ 1
tg (PF ε , N) = = (voir chapitre 3, paragraphe 4.1).
1 + e cos θ tgα
2. Conclusion
Nous pouvons en déduire, en employant le langage de l’optique, que le rayon
incident provenant d’un foyer se réfléchit à la surface de l’ellipse de sorte qu’il
passe par l’autre foyer et que les angles d’incidence et de réflexion sont égaux
et de signe contraire. Les opticiens disent que l’ellipse est « stigmatique »
pour son couple de foyers.
Remarque
u
P u’
F’ F
N
exercices
106
5. Anomalie et orbite
Posant :
LLI LLLI
I FP I F 'P
u = LLLI u' = LLLLLI
FP 2 F ' P2
LLI LLLI
dFP dF ' P
et notant que le vecteur tangent =
dϕ dϕ
LLI
I I dFP
il vient : ( u + u ') = 0
I dϕ I I
Le vecteur N est donc colinéaire à u + u ' et par conséquent porté par
la bissectrice.
exercices
107
This page intentionally left blank
Chapitre 6
3e loi de Kepler
et référentiels
μS = GMS
Le Soleil ayant pour masse MS = 1,988 9 1030 kg
μS = 1,327 1 1020 N.m2/kg (m3.s-2)
On a vu au chapitre 1 que la masse du Soleil est une constante dérivée. La
valeur plus précise de la constante d’accélération ou constante héliocentrique
de la gravitation est :
110
6. 3e loi de Kepler et référentiels
111
Guide de localisation des astres
Expressions numériques
Donnons une expression simple de
a3
T2 = (4π2/G). où 4π2/G = 5,916 51 1011 SI
M
dérivée des unités suivantes, précisées au chapitre introductif :
• l’année julienne (aj) = 365,25 24 3600 = 365,25 86 400
= 31 557 600 s
• l’unité astronomique (ua) = 149,597 870 61 109 m
• la masse du Soleil (MS) = 1,988 9 1030 kg
T(2aj) ≈ a 3( ua )
Remarque
Il n’est pas surprenant que le coefficient KS soit si proche de l’unité ; mais
pour obtenir rigoureusement KS = 1, il aurait par exemple fallu choi-
sir comme unités liées à la Terre :
• l’année sidérale TSid (unité notée sid) ;
• la distance moyenne Terre-Soleil aT = (TSid/2π)2/3 .(GMS)1/3
(unité notée at).
La 3e loi de Kepler, par rapport à ces unités et à la masse du Soleil,
s’écrit :
M( Ms ) MS
a 3( at )
et on obtient en toute rigueur : T(2Sid ) =
M( Ms )
112
6. 3e loi de Kepler et référentiels
Comme nous venons de le voir, les planètes du Soleil sont étudiées dans un
référentiel héliocentrique. Les satellites de la Terre, la Lune par exemple, sont
étudiés dans un référentiel géocentrique. Le mouvement d’un satellite d’une
planète quelconque sera étudié dans un référentiel lié à celle-ci.
Le problème est que ces référentiels masquent le caractère réciproque de la
gravitation universelle…
γA γB
A B
mA FAB FBA mB
113
Guide de localisation des astres
On déduit :
I I I I I
• m A Γ A + m BΓ B + m C Γ C + … + m K Γ K = 0 ,
ce qui signifie que la quantité d’accélération est nulle.
Effectuant une première intégration, il vient :
I I I I I I
• m A VA + m B VB + m C VC + … + m K VK = C tte = P , où P est une constante
d’intégration ; cette égalité traduit le fait que la quantité de mouvement est
constante.
Enfin, une dernière intégration conduit à :
LLLI LLLI LLLI LLLI I I
• m A OA +I m B OB + m C OC + … + m K OK = Pt + Q ,
le vecteur Q étant la deuxième constante d’intégration.
114
6. 3e loi de Kepler et référentiels
I I I I
On choisit P = 0 ; Q = 0 ;
l’origine O introduite lors de la deuxième intégration correspond alors au
barycentre G, d’où l’écriture finale :
LLLI LLLI LLLI LLLI I
• m A GA + m B GB + m C GC + …….+ m K GK = 0
On définit ainsi un référentiel barycentrique, c’est-à-dire dont l’origine
du repère d’espace est le centre d’inertie des corps A, B , C, … , K, et dont
l’orientation est rapportée à trois étoiles fixes. Un tel référentiel est bien sûr
galiléen. Un exemple est le référentiel de Copernic dont l’origine du repère
est le centre d’inertie du système solaire.
Le référentiel héliocentrique ou de Kepler est quant à lui constitué d’un
repère ayant pour origine le centre d’inertie du Soleil et en translation par
rapport à celui de Copernic. Considérer le référentiel de Kepler comme
galiléen revient donc à confondre le centre d’inertie du Soleil avec celui du
système solaire ; il est vrai qu’à lui seul, le Soleil représente 99,87 % de la
masse du système solaire !
115
Guide de localisation des astres
116
6. 3e loi de Kepler et référentiels
1 + m sat / m p m sat ⎛ M ⎞
= Λ donne = ⎜1 + S ⎟ Λ − 1
1 + MS / m p mp ⎝ mp ⎠
m sat ⎛ mp ⎞ mp
= ⎜1 + ⎟ Λ−
MS ⎝ MS ⎠ MS
117
Guide de localisation des astres
On notera que la cohérence avec la valeur des masses des satellites impose de
calculer le demi-grand axe à dix mètres près !
Prenons comme autre exemple la Lune, dont les données sont :
aL = 383 398 km TL = 27,3217 j mL = 1,230 % de mT
On peut vérifier que la valeur ci-dessus de la période de révolution lunaire ne
permet pas de retrouver le rapport mL / mT ; en revanche, la valeur TL = 27,1780 j
calculée à l’exercice E6-2 (hors perturbations) convient parfaitement.
118
6. 3e loi de Kepler et référentiels
exercices
E6-1 Système de deux corps dans le référentiel du centre
de masse
On considère deux corps de masses respectives m1 et m2 , qui dans le réfé-
LLLLI RIdu laboratoire,
rentiel LLLLI I muni d’un repère de centre O, ont pour position :
OM 1 = ρ1 et OM 2 = ρ2 . I LLLLLLI I I
La position relative des deux corps est : ρ = M 1 M 2 = ρ2 − ρ1
m1 m2
On posera : M = m1 + m2 , et m =
m1 + m2
I
Ces deux corps sont soumis
I Ià un champ
I de pesanteur g et à la force d’attrac-
tion gravitationnelle F = F21 = − F12 .
1. On se place dans le référentiel R .
a) Déterminer la position G du centre de masse (CM), sa vitesse et son accé-
lération.
b) Le champ de pesanteur est maintenant supposé négligeable : les deux corps
forment un système isolé. Montrer qu’il y a conservation de la quantité de
mouvement et l’exprimer en fonction de la vitesse du CM.
c) Écrire l’expression du moment cinétique et de l’énergie cinétique.
d) Montrer que le mouvement relatif des deuxI corps se ramène à celui d’un
corps unique de masse m soumis à la force F .
2. On se place dans le référentiel R G du CM.
a) Déterminer les positions des points M1 et M2 par rapport à G, et en déduire
les quantités de mouvement de chaque corps et du système.
LI
b) Donner l’expression du moment cinétique σ G .
c) Donner également l’expression de l’énergie cinétique EC .
3. Montrer que les grandeurs cinétiques obtenues dans le référentiel R du
laboratoire sont la somme d’un premier terme lié au mouvement du CM et
d’un deuxième lié au mouvement relatif des deux corps.
solutions
1. Référentiel R I LI
exercices
119
Guide de localisation des astres
LLLI I I
I d 2 OG 1 ⎛ d 2ρ1 d 2ρ 2 ⎞
ΓG = = m + m
m1 + m 2 ⎜⎝ dt 2 ⎟⎠
1 2
dt 2 dt 2
I I
d 2ρ1 I I d 2ρ 2 I I
Or m1 2 = m1 g − F et m 2 2 = m 2 g + F
dt dt
I I
Ce qui donne : Γ G = g → le mouvement du centre de masse est recti-
ligne, uniformément varié.
I
b) Quantité de mouvement p
LLLI I I
I d 2 OG d 2ρ1 d 2ρ 2
g ≈ 0 → = 0 → m1 2 + m 2 2 = 0
dt 2 dt dt
Après intégration, et en introduisant les quantités de mouvement de chaque
corps, on obtient :
I I
dρ dρ I I I
m1 1 + m 2 2 = p1 + p2 = p = C tte
dt dt
I I I
Or m1V1 + m 2 V2 = (m1 + m 2 )VG
LLLI
I I dOG
D’où : p = MVG = M = C tte
dt
LI
c) Moment cinétique σ , énergie cinétique ec
I I
I I dρ1 I dρ 2 I I I I
σ = ρ1 ∧ m1 + ρ2 ∧ m 2 = ρ1 ∧ p1 + ρ2 ∧ p2
dt dt
I 2 I 2
1 ⎛ dρ ⎞ 1 ⎛ dρ ⎞ 1 I 1 I
e c = m1 ⎜ 1 ⎟ + m 2 ⎜ 2 ⎟ = m1V12 + m 2 V22
2 ⎝ dt ⎠ 2 ⎝ dt ⎠ 2 2
d) Mouvement relatif I
Chacun des deux corps est soumis à la force d’attraction mutuelle F telle
que :
I I I I
d 2ρ d 2ρ
F = m 2 22 → m1 F = m1m 2 22
dt dt
I I I I
d 2ρ1 d 2ρ1
F = − m1 2 → m 2 F = − m 1m 2 2
dt dt
I I I
I m 1 m 2 ⎛ d 2ρ 2 d 2ρ1 ⎞ I d 2ρ
D’où, après addition, F = − → F = m 2
exercices
m1 + m 2 ⎜⎝ dt 2 dt 2 ⎟⎠ dt
2. Référentiel R G du CM
I I I
a) Positions de M1 et M2, quantités de mouvement P1 , P2 et P
LLLLI LLLLI LLLI I I I I
I m ρ + m 2 ρ2 m (ρ − ρ1 )
GM1 = OM1 − OG = ρ1 − 1 1 = − 2 2
m1 + m 2 m1 + m 2
120
6. 3e loi de Kepler et référentiels
LLLI I
⎛ LLLI m1 I ⎞ ⎛ dOG dρ ⎞
+ OG + ρ ∧ m2 +m ⎟
⎝ M ⎠ ⎜⎝ dt dt ⎠
Développant et simplifiant on obtient :
LLLI I
I LLLI dOG I I I dρ
σ = OG ∧ M + σ G avec σG = ρ ∧ m
dt dt
121
Guide de localisation des astres
• Énergie cinétique :
On développe ec sachant que :
I 2 LLLI 2 I 2 LLLI I
⎛ dρ1 ⎞ ⎛ dOG ⎞ ⎛ m 2 ⎞ ⎛ dρ ⎞
2
m2 dOG dρ
⎜⎝ ⎟⎠ = ⎜ ⎟ + ⎜⎝ ⎟⎠ − 2
dt ⎝ dt ⎠ ⎝ M ⎠ dt M dt dt
I 2 LLLI 2 I LLLI I
⎛ dρ2 ⎞ ⎛ dOG ⎞ ⎛ m1 ⎞ ⎛ dρ ⎞
2 2
m dOG dρ
⎜⎝ ⎟⎠ = ⎜ ⎟ + ⎜ ⎟ +2 1
dt ⎝ dt ⎠ ⎝ M ⎠ ⎝ dt ⎠ M dt dt
D’où l’expression :
LLLI 2 I 2
1 ⎛ dOG ⎞ 1 ⎛ dρ ⎞
e c = M⎜ + E C avec E C = m ⎜ ⎟
2 ⎝ dt ⎟⎠ 2 ⎝ dt ⎠
solutions
exercices
1. Mouvement et repère
Le seul mouvement possible est un mouvement de rotation autour du centre
de masse G, les centres M1 et M2 de la Planète et de son satellite étant alignés
avec celui-ci.
Le repère d’étude est le repère barycentrique.
I
Sur la figure 6.2, I est le vecteur unitaire orienté de M1 vers M2 .
122
6. 3e loi de Kepler et référentiels
I
a M2
F
G
r2
r1
M1
2. Rayons r1 et r2 , période T
m2 m1
G étant le barycentre, il vient : r1 = a r2 = a
m1 + m 2 m1 + m 2
Considérant par exemple le satellite : I mm I
• Il est soumis à la force d’attraction gravitationnelle F = −G 12 2 I
a
• la composante orthoradiale de son accélération est nulle
Γ = rθ ≡ 0 → θ = C tte = ω
θ 2
• la composante radiale a donc pour expression Γ ρ = − r2θ 2 = − r2 ω 2
(voir chapitre 3, paragraphe 1).
Le principe fondamental de la dynamique permet d’écrire :
mm
G 1 2 2 = m 2 r2 ω 2
a
Compte tenu de l’expression de r2 , on retiendra finalement
mm m1m 2
G 12 2 = aω2 = m aω 2
a m1 + m 2
exercices
123
Guide de localisation des astres
mT 1
rL = a = a → rL = 378 739,5 km
mT + mL 1+ k
Le CM du système Terre-Lune se trouve donc à 4 658 km du centre de la Terre
dont le rayon est de 6 378 km ; ainsi, on peut dire en première approximation
que la Lune tourne autour de notre planète qui elle reste immobile.
a3 a3
T = 2π = 2π → T = 27,1780 j (2 348 179 s)
Gm T (1 + k ) μ T (1 + k )
Dans l’approximation géocentrique, on obtient T = 27,3446 j ; compte tenu
des différentes perturbations, la valeur à retenir est en réalité T = 27,3217 j.
4. Couple Pluton-Charon, mP et mC
• La connaissance de la distance Charon-CM permet d’en déduire le rapport
de la masse du satellite à celle de la planète :
1 mC a
rC = a d’où k = = − 1 → k = 0,156 61
1+ k mP rC
• D’après la 3e loi de Kepler corrigée, posant M = mP + mC , on déduit :
a3
GM = 4 π 2 → GM = 955,222 109 m3.s-2
T2
et M = 143,156 1020 kg ; d’où
M
mP = → mP = 123,8 1020 kg
1+ k
kM
m C = km P = → mC = 19,4 1020 kg
1+ k
exercices
Remarques
1– Pour le système Terre-Lune on donne à la figure 6.3 la trajectoire autour
du Soleil du centre de masse G, laquelle est une ellipse, ainsi que l’allure
des trajectoires de notre planète et de son satellite, sachant que la Lune
n’orbite pas dans le plan de l’écliptique !
124
6. 3e loi de Kepler et référentiels
Rayons du Soleil
NL NL
PQ DQ
PL
Trajectoires
- du CM
- de la Lune
- de la Terre
125
This page intentionally left blank
Chapitre 7
Énergie et vitesse selon
la trajectoire
Nous introduisons différentes gran- Les deux premiers exercices concer-
deurs énergétiques et démontrons nent les satellites géostationnaires.
que pour un corps soumis à une force Nous déterminons les caractéristi-
d’origine gravitationnelle, l’énergie sur ques d’une orbite géosynchrone,
la trajectoire est une constante inver- puis comment faire passer le satellite
sement proportionnelle au grand axe. de l’orbite basse de lancement à son
Nous verrons que le signe de cette orbite opérationnelle : c’est le trans-
deuxième constante dynamique déter- fert d’orbite de Hohmann.
mine la nature de la conique, et que Les deux exercices suivants nous
la connaissance des deux constantes feront voyager dans notre système
dynamiques permet d’accéder à l’ex- solaire. Nous apprendrons à déter-
centricité. Il sera enfin montré que miner la trajectoire à énergie minimale
la trajectoire est parfaitement déter- d’une sonde vers Mars et la durée d’un
minée par la vitesse de lancement du tel voyage, qui ne peut se reproduire
satellite comparée aux vitesses mini- que tous les 780 jours ! Nous précise-
males de satellisation et de libération, rons également les positions de notre
respectivement appelées première et planète et de l’astre rouge au départ
deuxième vitesse cosmique. et à l’arrivée de la sonde. Nous calcu-
lerons enfin la vitesse à communiquer
pour un poser en douceur. Le dernier
exercice propose toutefois deux autres
issues à ce voyage : soit la sonde
passe à une distance donnée de Mars
et, profitant de cet effet catapulte, elle
continue son voyage ; soit on la satel-
lise en diminuant sa vitesse…
Guide de localisation des astres
128
7. Énergie et vitesse selon la trajectoire
Énergie cinétique EC
On calcule l’énergie liée au mouvement, en utilisant la formule de Binet rela-
tive à la vitesse ;
1
EC = mV 2
2
soit :
⎡ 2
⎛ d ⎛ 1⎞ ⎞ ⎤
2
1 2 ⎛ 1⎞
E C = mC ⎢⎜ ⎟ + ⎜ ⎜ ⎟ ⎟ ⎥
2 ⎢⎣⎝ ρ ⎠ ⎝ dθ ⎝ ρ ⎠ ⎠ ⎥⎦
1 1 + e cos θ
Or, = et C 2 = pμ
ρ p
d’où le résultat :
1 1 + 2e cos θ + e 2
EC = mμ
2 p
• Donnons d’autres expressions de l’énergie cinétique :
1 1 + e cos θ + e (e + cos θ)
EC = mμ
2 p
d’où on déduit :
1 1 e (e + cos θ) 1
E C = − E P + mμ soit en mettant − E P en facteur,
2 2 p 2
e (e + cos θ) ⎤
E C = − E P ⎡⎢1 +
1
2 ⎣ 1 + e cos θ ⎥⎦
Le deuxième terme du crochet a été exprimé en fonction de l’anomalie excen-
trique ϕ, ce qui permet l’écriture condensée :
1
E C = − E P (1 + e cos ϕ )
2
Notons qu’en dérivant par rapport au temps l’équation de Kepler établie au
chapitre 5, on a :
2π
= ϕ (1 − e cos ϕ )
T
De ces deux dernières relations faisant intervenir l’anomalie excentrique, on
déduit :
T dϕ EP
=
π dt EP + EC
129
Guide de localisation des astres
Énergie mécanique E
La variation de l’énergie cinétique s’obtient en introduisant
I la vitesse dans
l’expression du travail élémentaire de la somme F des forces auxquelles le
corps est soumis :
I
I LI
L dV I ⎛1 I
δW = F dl = m Vdt = d mV 2 ⎞ = dE C
dt ⎝2 ⎠
Si on considère d’autre part que les forces agissantes dérivent d’un poten-
tiel, c’est-à-dire qu’en particulier il n’y a pas de force de frottement, il vient
alors :
I LI
L
F dl = − dE P = dE C
Or l’énergie mécanique E du corps assimilé à un point matériel est la somme
des énergies cinétique et potentielle :
E = EC + EP
d’où : dE ≡ 0
Par conséquent l’énergie mécanique est une constante du mouvement ; c’est
la deuxième constante dynamique. La première constante est celle de la loi
des aires, comme vu au chapitre 3. Compte tenu des dernières expressions
de EC, on obtient immédiatement :
1 e2 − 1
E = mμ
2 p
Rappelons que : - lorsque e < 1 on a établi p = a(1 − e2) = b2/a
- lorsque e > 1 on a posé b2 = c2 − a2
d’où p = − a(e2 − 1) = − b2/a et on a vu que dans ce cas,
a est devenu négatif, donc p reste positif.
On en déduit la deuxième expression de l’énergie totale, valable quelque soit
le type de conique :
1 μ
E = − m
2 a
130
7. Énergie et vitesse selon la trajectoire
Énergie
Excentricité Courbe Demi-grand axe Remarques
mécanique
e=0 cercle a = p = RP E<0 EC = Cte
0<e<1 ellipse e→1 ⇒a→+∞
–
E<0 EC < |EP|
e=1 parabole |a| = ∞ E=0 θ → ± π ⇒ EC → 0
e>1 hyperbole e → 1+ ⇒ a → − ∞ E>0 EC > |EP|
Énergie fournie EF
Cette notion d’énergie fournie trouve son intérêt dans l’étude des satellites
artificiels. Au repos, au sol, l’énergie potentielle est :
μ
E P (R T ) = − m où RT est le rayon terrestre.
RT
Si on néglige la vitesse d’entraînement de la Terre, l’énergie fournie a pour
expression :
1 μ μ
E F = E − E P (R T ) = − m + m
2 a RT
soit : GM T ⎛ R ⎞
EF = m 1 − T⎟
R T ⎜⎝ 2a ⎠
MT étant la masse de la Terre.
Dans le cas d’une orbite basse circulaire, a ≈ RT et bien sûr
1 GM T
EF ≈ m = EC
2 RT
Pour une orbite circulaire de rayon Rsat,
GM T ⎛ RT ⎞
EF = m ⎜ 1−
RT ⎝ 2R sat ⎟⎠
131
Guide de localisation des astres
soit, en explicitant V2 :
⎛ dρ ⎞ ⎛ dθ ⎞ μ
2 2
2E
+ ρ2 −2 =
⎝ dt ⎠ ⎝ dt ⎠ ρ m
• La première constante dynamique provient du fait qu’un mouvement à accé-
lération centrale satisfait à la loi des aires :
dθ dθ C dρ dρ C d ⎛ 1⎞
ρ2 = C ; ainsi, = 2 et = = −C ⎜ ⎟
dt dt ρ dt dθ ρ 2
dθ ⎝ ρ ⎠
2
1 ⎛ du ⎞ 2E
Posant u = , on obtient la relation : C2 + C 2 u 2 − 2 μu =
ρ ⎝ dθ ⎠ m
μ μ2 μ2
2
⎛ du ⎞ 2E
que l’on transforme comme suit : + u2 − 2 2 u + 4 = 4 +
⎝ dθ ⎠ C C C mC2
d’où l’équation différentielle finale :
μ μ2
2 2
⎛ du ⎞ ⎛
+ u − 2⎞ = 4 +
2E
⎝ dθ ⎠ ⎝ C ⎠ C mC2
Une solution réelle est possible pour
1 μ2
E ≥ − m 2
2 C
Elle s’écrit :
μ μ2 2E
u = + + cos(θ − θ0 )
C2 C 4
mC2
soit finalement, en réintroduisant le rayon ρ
1 μ ⎡ 2E C 2 ⎤
= 2 ⎢1 + 1 + cos(θ − θ0 ) ⎥
ρ C ⎣⎢ m μ 2
⎥⎦
La solution est effectivement une conique dont un foyer est l’origine des
rayons vecteurs, et pour laquelle la constante d’intégration θ0 est l’angle du
grand axe avec l’axe origine Ox.
1 1 + e cos(θ − θ0 )
Par identification avec = on déduit :
ρ p
132
7. Énergie et vitesse selon la trajectoire
C2
paramètre p =
μ
2E C 2
excentricité e = 1+
m μ2
2 Vitesses et trajectoires
133
Guide de localisation des astres
134
7. Énergie et vitesse selon la trajectoire
V0 →
V∞ = V 0 e − 1 vitesse
e +1
de lancement
V∞ = 0
V0 < Vmin →
trajectoire
balistique
V0
V0 = Vmin →
orbite circulaire
rasante
Rp
Vmin < V0 < Vlib →
F orbite elliptique
Vπ = V 0 1 − e
1+ e
V0 = Vlib →
trajectoire
parabolique
V0 > Vlib →
trajectoire
hyperbolique
135
Guide de localisation des astres
Il vient :
ρ0 → ρh = R P + h (donnée)
1
ρ π = (R P + h) 2
(ellipse)
⎛ Vlibh ⎞
⎜⎝ V ⎟⎠ − 1
h
V0 → Vh (donnée)
⎡⎛ Vlibh ⎞ 2 ⎤ R + h (ellipse)
Vπ = Vh ⎢⎜ ⎟ − 1⎥ = Vh P
⎢⎣⎝ Vh ⎠ ⎥⎦ ρπ
2
⎛V ⎞
Vθ → θm = V∞ = Vh 1 − ⎜ libh ⎟ (hyperbole)
⎝ Vh ⎠
136
7. Énergie et vitesse selon la trajectoire
1 2 ⎡ ⎛ Vh ⎞ ⎤
2
E
= − Vmin h⎢ 2 − ⎜⎝ V ⎟⎠ ⎥
m 2 ⎢⎣ min h ⎥⎦
On retiendra comme expression de l’énergie massique :
1 2 ⎡⎛ Vlibh ⎞ ⎤
2
E
= − Vh ⎢⎜ ⎟ − 1⎥
m 2 ⎢⎣⎝ Vh ⎠ ⎥⎦
et on vérifie que dans le cas de l’ellipse, c’est-à-dire, pour Vh < Vlib h :
E 1
= − Vh Vπ
m 2
Période
1μ
3/ 2
⎛ a⎞ E
e
La 3 loi de Kepler conduit à T = 2 πμ ⎜ ⎟ Or = −
⎝ μ⎠ m 2a
Ainsi, dans le cas d’une trajectoire elliptique, on a :
3/ 2
⎛
3/ 2
⎛ m ⎞ 1 ⎞
T = 2 πμ ⎜ = 2 π GM P ⎜ 2 2⎟
⎝ 2 E ⎟⎠ ⎝ Vlibh − Vh ⎠
On retiendra comme expression de la période :
2
Vmin
T = 2 π (R P + h) h
(V )
3/ 2
2
lib h − Vh2
• Avant de clore ce paragraphe, remarquons que pour certains problèmes les
données sont l’altitude de lancement h et le point le plus éloigné du foyer ρπ .
On détermine alors Vh par la relation déduite de l’expression de ρπ :
1
Vh = Vlibh
R +h
1+ P
ρπ
137
Guide de localisation des astres
Remarques
1– Les exercices E7-1 et E9-3 concernent des satellites en orbite géos-
tationnaire, c’est-à-dire sur une orbite géosynchrone située dans le plan
équatorial et dont l’excentricité est nulle. Avec l’exercice E11-1 on traitera
d’orbites semi-synchrones...
2– C’est le britannique Sir Arthur Charles Clarke, né en 1917, qui, dans
un article publié dans le Wireless World en 1945, proposa d’utiliser des
satellites en orbite géostationnaire pour les télécommunications.
Le calcul d’une telle orbite, souvent baptisée orbite de Clarke, avait été
effectué dans un ouvrage, Le Problème du Voyage Spatial, publié à Berlin
en 1928. L’auteur en était Herman Potocnik (1892-1929) également connu
sous son pseudonyme de Noordung. Le russe Konstantin Tsiolkovsky (1857-
1935) avait déjà eu l’idée de ce type d’orbite.
Le premier satellite de télécommunication vraiment stationnaire, c’est-à-
dire dont l’inclinaison de l’orbite par rapport au plan équatorial est négli-
geable, est le Syncom 3, lancé par les États-Unis en 1964 pour les jeux
olympiques de Tokyo.
Actuellement, avec plus de 180 satellites en exploitation, on assiste à l’en-
combrement de l’orbite géostationnaire.
138
7. Énergie et vitesse selon la trajectoire
exercices
E7-1 Satellite géostationnaire
1.
a) Calculer l’altitude et la vitesse d’un satellite géostationnaire en utilisant
les formules générales dans lesquelles e → 0.
b) Établir directement les formules dans le cas particulier.
2.
a) Calculer les énergies massiques (potentielle, cinétique et mécanique).
b) Dresser un tableau de ces trois énergies en fonction de l’état du satellite :
repos au sol, en orbite géostationnaire ; en déduire les variations de ces éner-
gies.
3. Donner l’expression littérale de l’énergie fournie en fonction des rayons
terrestre et orbital d’une part et des vitesses cosmique (1re) et orbitale d’autre
part.
Données numériques :
– constante géocentrique de la gravitation μT = 3,986 004 418 1014 m3.s-2
– rayon équatorial RT = 6,378 14 10 6 m
– jour sidéral terrestre T = 86 164,095 6 s.
solutions
1.
a) Rayon Rsat , altitude Hsat et vitesse Vsat
• Considérations préliminaires :
μT = 398 600,4418 km3.s-2 RT = 6 378,14 km
La période du satellite géostationnaire est par définition égale à la période de
rotation de la Terre : le satellite est synchrone de la Terre, c’est-à-dire géosyn-
chrone et T = 86 164,0956 s.
L’orbite, dans le plan équatorial, tend vers un cercle : dans les formules géné-
rales e → 0 .
• Calcul de Rsat
a3
La 3e loi de Kepler s’écrit : μ T = 4 π 2
exercices
T2
p
Or quand e → 0 , a = → p = Rsat
1 − e2
⎛ μT 2 ⎞
1/ 3
139
Guide de localisation des astres
• Calcul de Vsat
Les formules de Binet ont conduit à
μ μ 1 + e 2 + 2e cos θ
V 2 = (1 + e 2 + 2e cos θ) =
p a 1 − e2
1/ 3
⎛ μ T2 ⎞
Or a = ⎜ T 2 ⎟ on déduit donc :
⎝ 4π ⎠
⎛ μ 1 + e 2 + 2e cos θ
2/3
V 2 = 2π T ⎞
⎝ T⎠ 1 − e2
⎛ μ ⎞
1/ 3
GM T
Il faut au préalable calculer E P (R T ) = − m , qui représente l’énergie
RT
potentielle du satellite au sol.
GM T
AN par unité de masse : − = − 62,495 MJ/kg
RT
140
7. Énergie et vitesse selon la trajectoire
Énergie EP GMT EC 1 2 E EP EC
= – = V = +
État du satellite m R m 2 m m m
Repos au sol (*)
R = RT V supposée nulle − 62,495 0 − 62,495
Orbite géostationnaire
R = Rsat V = Vsat − 9,454 + 4,727 − 4,727
Variation d’énergie + 53,041 + 4,727 + 57,768
(*) Voir plus loin exercice E9-2.
3. Énergie fournie EF
EF 1 GM T GM T
E F = E − E P (R T ) → = − +
m 2 R sat RT
EF GM T ⎛ RT ⎞
→ = ⎜ 1−
m RT ⎝ 2R sat ⎟⎠
GM T
La première vitesse cosmique est : Vmin =
RT
GM T
d’autre part : Vsat =
R sat
Vsat2 R EF 1
On en tire : 2
= T et = Vmin
2
− Vsat2
Vmin R sat m 2
On remarquera que lorsque le rayon de l’orbite circulaire passe de RT à l’infini :
→ Vsat décroît de Vmin à 0
1 2
→ EF/m croît de
2
Vmin à Vmin
2
141
Guide de localisation des astres
2.
a) Quelles sont les énergies massiques sur les orbites basse, de transfert
et géostationnaire ? En déduire les accroissements d’énergie massique à
communiquer au périgée et à l’apogée.
b) Quelle est la durée du transfert de Hohmann ?
3. La vitesse atteinte à l’apogée est inférieure de ΔVt2 à la vitesse initiale-
ment calculée.
a) Expliquer pourquoi l’origine de cet écart est dû à une altitude de lancement
insuffisante ; calculer Δh1.
b) Calculer les accroissements de vitesse effectivement communiqués au satel-
lite afin qu’il soit quand même géostationnaire.
c) Calculer les nouveaux accroissements d’énergie massique nécessaires et
les comparer aux valeurs nominales.
Données numériques :
– μT = 3,986 004 418 105 km3.s-2 ; rayon terrestre RT = 6 378,14 km
– altitude de l’orbite basse h1 = 200 km
– altitude de l’orbite géostationnaire h2 = 35 786,03 km
– ΔVt2 = 5 m.s-1.
solutions
1. Vitesses et orbites
a) Vitesses V1 et V2
Orbite circulaire basse :
μT 1
V1 = = Vmin → V1 = 7,7843 km.s-1
R T + h1 h
1+ 1
RT
Orbite géostationnaire :
μT 1
V2 = = Vmin → V2 = 3,0747 km.s-1
R T + h2 h
1+ 2
RT
b) Orbite de transfert
Demi-grand axe :
(R T + h1 ) + (R T + h 2 ) h + h2
exercices
at = = RT + 1 → at = 24 371 km
2 2
La trajectoire étant elliptique,
⎛2 1⎞ μT ⎛ R + h⎞
1/ 2
Vt2 = μ T ⎜ − ⎟ avec ρ = RT + h → Vt = 2− T
⎝ ρ at ⎠ ⎜
RT + h ⎝ a t ⎟⎠
142
7. Énergie et vitesse selon la trajectoire
ΔV2
Δ V1
RT h1
exercices
h2
143
Guide de localisation des astres
⎛ 1 1 ⎞
(Vt 2 − ΔVt 2 )2 = 2μ T ⎜ − ⎟ < Vt22
⎝ R T + h 2 2a t ⎠
'
Par conséquent : 2a 't < 2a t → le satellite n’a été lancé qu’à l’altitude h1 − Δh1
L’accroissement de vitesse au périgée a toutefois été corrigée pour que le
satellite atteigne quand même l’altitude h2.
144
7. Énergie et vitesse selon la trajectoire
E 't 1
= − Vt'1Vt'2 = − 8,1857 MJ/kg
m 2
Les nouveaux accroissements d’énergie massique à communiquer sont :
ΔE1' E' E'
– au périgée : = t − 1 = + 22,3318 MJ/kg
m m m
soit + 212 kJ/kg par rapport au nominal ;
145
Guide de localisation des astres
ΔE '2 E E'
– à l’apogée : = 2 − t = + 3,4588 MJ/kg
m m m
soit + 8 kJ/kg par rapport au nominal.
146
7. Énergie et vitesse selon la trajectoire
solutions
1. Caractérisation de l’ellipse et vitesses d’entraînement
a) Demi-grand axe at , excentricité e, paramètre p
ρT + ρ M
at = → at = 1,888 108 km
2
p p ρ 1− e
Au périhélie ρT = ; à l’aphélie ρ M = , d’où T = ;
1+ e 1− e ρM 1+ e
on déduit :
ρ M − ρT
e = → e = 0,208
ρ M + ρT
2ρ M
p = ρT (1 + e ) , or 1 + e = ; on déduit l’expression du paramètre :
ρ M + ρT
ρT ρ M
p = 2 → p = 1,208 ua = 1,807 108 km
ρT + ρ M
b) Vitesses VP au périhélie et VA à l’aphélie, constante C des aires
⎛ 2 1⎞
Les vitesses se déduisent de la formule V 2 = μ ⎜ − ⎟
⎝ ρ a⎠
⎛ 2 1⎞ μ S 2ρ M
VP = μS ⎜ − ⎟ = → VP = 32,73 km.s-1
⎝ ρT at ⎠ ρT ρT + ρ M
⎛ 2 1⎞ μ S 2ρT
VA = μS ⎜ − ⎟ = → VA= 21,48 km.s-1
⎝ ρM at ⎠ ρ M ρT + ρ M
ρT ρ M
C = μS p = 2μ S → C = 4,896 109 km2.s-1
ρT + ρ M
c) Vitesses d’entraînement VT et VM
• La Terre tourne autour du Soleil en un an (vitesse ωT, période TT) et dans
le référentiel héliocentrique il faut tenir compte de la vitesse d’entraînement
communiquée à la sonde :
exercices
2π
VT = ρT ω T = ρT = 1,496 108 2π/(365,26 86400) km.s-1
TT
→ VT = 29,78 km.s –1
On vérifie, voir exercice E3-2, que 29,29 km.s-1 < VT < 30,29 km.s-1 ; par ailleurs,
μS
on rappelle que VT = : l’année sidérale est une donnée redondante !
ρT
147
Guide de localisation des astres
TM2 TT2
Or, la troisième loi de Kepler permet d’écrire = , d’où :
ρ3M ρ3T
3/ 2
⎛ρ ⎞
TM = ⎜ M ⎟ TT → TM = 1,881 TT
⎝ ρT ⎠
AN : VM = (1,524/1,881) VT → VM = 24,13 km.s –1
On note que 21,97 km.s-1 < VM < 26,50 km.s-1 et on rappelle que
μS ρT
VM = = VT
ρM ρM
1
Vlib. h = 2 Vmin → Vlib.h = 10,93 km.s –1
1 + h / RT
148
7. Énergie et vitesse selon la trajectoire
Remarque
VI = Vlib2 . h + V∞2
π − θM D D ⎛ρ ⎞ 1 ⎛ ρ + ρM ⎞ ⎛ ρT ⎞ 1 ⎛ ρ + ρM ⎞
= M = M⎜ T⎟ = ⎜ T = ⎜ T
2π TM TT ⎝ ρM ⎠ 2 ⎝ 2ρT ⎟⎠ ⎜⎝ ρ ⎟⎠
M 2 ⎝ 2ρ M ⎟⎠
soit :
⎡ ⎛ ρT + ρ M ⎞ ⎤
3/ 2
θM = π ⎢1 − ⎜ ⎥ → θM = 44°,36
⎢⎣ ⎝ 2ρ M ⎟⎠ ⎥⎦
149
Guide de localisation des astres
Sonde
Mars
position au départ
θT θM
Mars Terre
position à l’arrivée Soleil position au départ
Terre
position à l’arrivée
Remarques
1– Ce problème donne un aperçu de la méthode de calcul de la trajectoire
à énergie minimale ; il devra être repris en tenant compte du fait que les
orbites réelles de la Terre et de Mars sont elliptiques et qu’elles ne sont
pas dans un même plan. Il faudra alors connaître la date du lancement
correspondant à une configuration initiale des deux planètes par rapport au
Soleil, qui ne se produit que tous les 780 jours (voir exercice E1-2).
exercices
150
7. Énergie et vitesse selon la trajectoire
151
Guide de localisation des astres
solutions
1. Trajectoire relative
a) Les vitesses de la sonde et de Mars sont parallèles :
V∞ = VM − VA → V∞ = 2,65 km/s
Par définition r = μ M / V∞2 → r = 6 099 km
b) Par rapport à Mars, la sonde en P est considérée à l’infini et son énergie
mécanique est simplement égale à son énergie cinétique :
E 1
= V∞2 > 0
m 2
axe
polaire
V’∞
D
θ T
Planète Mars
Paramètre d’impact d
exercices
V∞
P
152
7. Énergie et vitesse selon la trajectoire
On retiendra :
ρmin = − r + d2 + r 2 → ρmin = 3 961 km
Notons que l’altitude hmin de la sonde est :
hmin = ρmin − RM → hmin = 564 km
153
Guide de localisation des astres
L’angle des deux vecteurs vitesse à l’infini est égal à l’angle D des deux
asymptotes. Soit θ∞ l’angle polaire rendant le rayon infini, il satisfait à
1 + e cos θ∞ = 0 ; d’où :
1
cos θ∞ = − → θ∞ = ± 127°,32
exercices
e
D’autre part :
D = π − 2( π − θ∞) = 2θ∞ − π ; compte tenu que e2−1 = p/r ,
on déduit les égalités :
2 p−r d2 − r 2
cos D = − cos 2θ∞ = 1 − 2 cos 2 θ∞ = 1 − = =
e2 p+r d2 + r 2
154
7. Énergie et vitesse selon la trajectoire
2 pr 2dr
sin D = = 2
p+r d + r2
D e2 − 1 D 1
Pour l’angle moitié, on obtient facilement cos = , sin = , et
2 e 2 e
D 1 r r μ
tg = = = = M2
2 e −1
2
p d dV∞
D’une façon ou d’une autre, on calcule D = 74°,64
I I
b) Avant l’interaction, les vecteurs VA et VM étaient colinéaires ; après l’in-
teraction
I' I ' on aI la relation :
VA = V∞ + VM
V∞’
α
’
VA
VM
sin α sin D
= ; exprimant le sinus on en déduit :
V∞ VA'
V 2dr 2drV∞
sin α = ∞ 2 =
VA d + r
' 2
(V∞ + VM )(d + r 2 )2 − 2V∞VM (d 4 − r 4 )
2 2 2
→ α = 6°,225
155
Guide de localisation des astres
I
La connaissance de VA' à l’aphélie permet de déterminer la nouvelle trajec-
toire de la sonde (voir exercice E8-3) .
4. Satellisation
a) Sur l’orbite circulaire,
μM
V 'T = → V 'T = 3, 29 km / s or VT = 5,35 km/s
ρmin
Il faut donc freiner la sonde par une poussée dirigée vers l’avant, afin qu’elle
perde la vitesse :
Remarques
1– Pour un observateur Martien, l’énergie mécanique de la sonde n’a pas
varié : la sonde s’éloigne à l’infini avec la même vitesse sur trajectoire
qu’elle avait en venant de l’infini.
Pour un observateur lié au Soleil, Mars, en interagissant avec la sonde, a
modifié la vitesse de celle-ci en grandeur et direction. Dans ce problème,
la planète rattrape la sonde, qui voit en conséquence sa vitesse s’accroître
de 2,09 km/s et dévier de plus de 6° ; cela se fait au détriment de l’énergie
de la planète : cette perte infinitésimale d’énergie provenant du travail des
forces gravitationnelles non nul dans le repère de Copernic.
exercices
156
7. Énergie et vitesse selon la trajectoire
exercices
157
This page intentionally left blank
Chapitre 8
Formules pour
mouvement elliptique
Dans les problèmes liés aux trajec- deuxième, la caractérisation se fait en
toires et aux positions à un instant fonction de la période et des vitesses
donné des planètes ou satellites sur au péricentre et à l’apocentre. Le
celles-ci, il convient de disposer rapi- troisième exercice nous propose un
dement des formules adaptées afin retour sur la comète de Halley, déjà
de s’éviter de fastidieux calculs. Tout rencontrée à l’exercice E5-1, et dont
au long des chapitres précédents, on précise la trajectoire à partir de la
nous nous sommes efforcés de fournir simple connaissance du rayon et du
de nombreuses relations que nous vecteur vitesse en un point donné.
complétons ici par un formulaire en Il sera suivi d’un dernier exercice qui
première partie du présent chapitre. La traitera de la trajectoire héliocen-
deuxième partie permet d’accéder au trique de la sonde Giotto à partir de
périmètre de la trajectoire elliptique en la connaissance des rayons vecteurs
fonction de l’excentricité et à l’aire d’un au départ de la Terre et au passage
secteur balayé par le rayon polaire. au voisinage de la comète. Il s’agit en
fait du problème de Lambert qui sera
Si les deux premiers exercices viennent résolu à l’aide d’un tableur.
compléter ce formulaire, les suivants
vont au-delà en nous permettant de Ce chapitre termine une première
déterminer des trajectoires à partir partie de l’ouvrage où nous avons
d’un minimum de données. étudié en coordonnées polaires les
Au premier exercice, nous établissons mouvements plans découlant d’une
un ensemble de formules permet- accélération centrale. Il convient main-
tant de caractériser un mouvement tenant de situer ces mouvements dans
elliptique en fonction des rayons au d’autres repères appropriés : ce sera
péricentre et à l’apocentre, ainsi que principalement l’objet des trois chapi-
de la vitesse au péricentre. Avec le tres qui suivront.
Guide de localisation des astres
→ a3
μ = 4π2
T2
• Énergie
1 μ μ a2
E = − m avec = 4π2 2
2 a a T
1 ⎛ 2 πa ⎞
2
→ 1
E = − m = − mV0 Vπ
2 ⎝ T ⎠ 2
Remarque
ρ0 + ρ π
T = 2π 2
V0 Vπ
160
8. Formules pour mouvement elliptique
C = μp avec p = a (1 − e 2 ) → C = μa 1 − e 2
• Trajectoire
a (1 − e 2 )
(voir paragraphe 1.1) → ρ =
1 + e cos θ
On déduit :
ρ0 + ρ π
ρ0 = a (1 − e ) ρ π = a (1 + e ) → = a
2
• Vitesses
μ
V = 1 + e 2 + 2e cos θ
p
μ 1 + e 2 + 2e cos θ
→ V =
a 1 − e2
On déduit :
μ 1+ e μ 1− e μ
V0 = Vπ = → V0 Vπ =
a 1− e a 1+ e a
• Période
a3 a
T2 = 4 π 2 → T = 2πa
μ μ
• Énergie
1 μ
→ E = − m
2 a
161
Guide de localisation des astres
• Trajectoire
a (1 − e 2 ) m 1 − e2
ρ = → ρ = μ
1 + e cos θ 2 E 1 + e cos θ
On déduit :
m m
ρ0 = μ (1 − e ) ρπ = μ (1 + e )
2E 2E
→ ρ0 + ρ π m
= μ
2 2E
• Vitesses
μ 1 + e 2 + 2e cos θ
V =
a 1 − e2
2E 1 + e 2 + 2e cos θ
→ V =
m 1 − e2
On déduit :
2E 1+ e 2E 1− e
V0 = Vπ =
m 1− e m 1+ e
2E
→ V0 Vπ =
m
• Période
a ⎛ m ⎞
3/ 2
T = 2πa → T = 2 πμ ⎜
μ ⎝ 2 E ⎟⎠
162
8. Formules pour mouvement elliptique
On déduit :
2/3 2/3
⎛ T⎞ ⎛ T⎞
ρ0 = μ1 / 3 (1 − e ) ρπ = μ1 / 3 (1 + e )
⎝ 2π ⎠ ⎝ 2π ⎠
ρ0 + ρ π
2/3
⎛ T⎞
→ = μ 1/ 3
2 ⎝ 2π ⎠
• Vitesses
2E 1 + e 2 + 2e cos θ
V =
m 1 − e2
−1 / 3
⎛ T ⎞ 1 + e 2 + 2e cos θ
→ V = ⎜
⎝ 2 πμ ⎟⎠ 1 − e2
On déduit :
−1 / 3 −1 / 3
⎛ T ⎞ 1+ e ⎛ T ⎞ 1− e
V0 = ⎜ Vπ = ⎜
⎝ 2 πμ ⎟⎠ 1− e ⎝ 2 πμ ⎟⎠ 1+ e
−1 / 3
⎛ T ⎞
→ V0 Vπ = ⎜
⎝ 2 πμ ⎟⎠
• Énergie
1/ 3 −2 / 3
m ⎛ T ⎞ E 1⎛ T ⎞
= ⎜ → = − ⎜
2E ⎝ 2 πμ ⎟⎠ m 2 ⎝ 2 πμ ⎟⎠
E
1.5. Relations en fonction de la période T, de l’énergie massique m ,
de l’excentricité e
• Constante des aires
1/ 3 1/ 3
⎛ T⎞ m ⎛ T ⎞
C = μ2/3 1 − e2 avec = ⎜
⎝ 2π ⎠ 2E ⎝ 2 πμ ⎟⎠
T 2E
→ C = 1 − e2
2π m
• Trajectoire
m 1 − e2 T 2 E 1 − e2
ρ = μ → ρ =
2 E 1 + e cos θ 2 π m 1 + e cos θ
On déduit :
T 2E T 2E
ρ0 = (1 − e ) ; ρ π = (1 + e )
2π m 2π m
163
Guide de localisation des astres
ρ0 + ρ π T 2E
→ =
2 2π m
• Vitesses
(voir 1.3)
2E 1 + e 2 + 2e cos θ
→ V =
m 1 − e2
2E 1+ e 2E 1− e
V0 = Vπ =
m 1− e m 1+ e
2E
→ V0 Vπ =
m
• Coefficient de l’accélération centrale
3/ 2
⎛ m ⎞ T ⎛2E ⎞
3/ 2
T = 2 πμ ⎜ → μ =
⎝ 2 E ⎟⎠ 2 π ⎜⎝ m ⎟⎠
B
M’
θ
dθ
θ A
164
8. Formules pour mouvement elliptique
Application à l’ellipse
Expression en coordonnées polaires
p
Le rayon vecteur a pour expression ρ(θ) =
1 + e cos θ
pe sin θ
de dérivée ρ '(θ) =
(1 + e cos θ)2
La longueur de l’arc d’ellipse est donc :
θB 1 + e 2 + 2e cos θ
l AB = p ∫ dθ
θA (1 + e cos θ)2
Or la formule de Binet relative à la vitesse nous a permis d’établir que :
C
V(θ) = 1 + e 2 + 2e cos θ C étant la constante des aires.
p
Ainsi, faisant apparaître également ρ(θ), il vient :
1 θB 2
C ∫θA
l AB = ρ (θ)V(θ) dθ
165
Guide de localisation des astres
ϕB
d’où l’on tire : l AB = ∫
ϕA
a 2 sin 2 ϕ + b2 cos 2 ϕ dϕ
Exprimons cette intégrale en fonction du paramètre p et de l’excentricité e.
Notant que :
p
a 2 sin 2 ϕ + b2 cos 2 ϕ = a 1 − e 2 cos 2 ϕ , avec a =
1 − e2
il en découle les égalités :
ϕB p ϕB
l AB = a ∫ 1 − e 2 cos 2 ϕ dϕ = ∫ 1 − e 2 cos 2 ϕ dϕ
ϕA 1 − e2 ϕA
Approche du périmètre
Les fonctions élémentaires existantes ne permettent pas d’obtenir une formule
exacte donnant le périmètre de l’ellipse.
Obtention par développement en série d’une intégrale elliptique
• Partant de l’expression de la longueur d’un arc d’ellipse en fonction de l’ano-
malie excentrique ϕ, et introduisant le paramètre t tel que :
166
8. Formules pour mouvement elliptique
π
t = ϕ−
2
on a bien sûr sin ϕ = cos t et cos ϕ = − sin t, et on déduit :
tB tB
l AB = ∫ tA
a 2 cos 2 t + b2 sin 2 t dt = a ∫
tA
1 − e 2 sin 2 t dt
π⎡
2 2 2
⎛ 1 ⎞ 2 1 ⎛ 1.3 ⎞ 4 1 ⎛ 1.3.5 ⎞ 6 1 ⎛ 1.3.5.7 ⎞ 8
2
⎤
E (e ) = ⎢1 − 1⎝ ⎠ e − ⎝ e − e − e − ...⎥
2⎣ 2 3 2.4 ⎠ ⎝
5 2.4.6 ⎠ ⎝
7 2.4.6.8 ⎠ ⎦
2
1 ⎛ (2n)! ⎞ 2 n
et le terme général du crochet a pour expression : − e
2n − 1 ⎜⎝ (2n n!)2 ⎟⎠
Écrivant sous forme de fractions rationnelles irréductibles les coefficients des
puissances de e, le périmètre pourra se calculer comme suit :
⎡ ⎛1
2
1 ⎛ 1.3 ⎞ 4
2
P = π(a + b) ⎢1 + 1 ⎞ u 2 + 2 u
⎣ ⎝ 2 ⎠ 3 ⎝ 2.4 ⎠
2
1 ⎛ 1.3.5 ⎞ 6
2
1 ⎛ 1.3.5.7 ⎞ 8 ⎤
+ u + u − ⋅ ⋅ ⋅⎥
5 ⎝ 2.4.6 ⎠
2
7 ⎝ 2.4.6.8 ⎠
2
⎦
167
Guide de localisation des astres
soit :
P 1 1 1
= 1 + 2 u 2 + 2 u 4 + 2 u6
π(a + b) 2 8 16
25 8 49 10 441 12
+ 2
u + 2
u + u −⋅⋅⋅
128 256 1024 2
• D’autre part Ramanujan, mathématicien indien du siècle dernier, donna des
formules approximatives du périmètre en fonction de a, b et éventuellement
(a − b)2
de h = u 2 =
(a + b)2
– Première formule de Ramanujan :
P ≈ π ⎡⎣3(a + b) − (a + 3b)(3a + b) ⎤⎦ = π (a + b) ⎡⎣3 − 4 − h ⎤⎦
– Deuxième formule de Ramanujan :
⎡ 3h ⎤
P ≈ π (a + b) ⎢1 +
⎣ 10 + 4 − 3h ⎥⎦
On déduit que la Terre parcourt les 939,886 465 106 km de son orbite à la
vitesse moyenne de 29,7827 km/s (à comparer à Vb = 29,7847 km/s : voir
exercice E3-2).
On vérifiera que, vu les faibles excentricités des planètes, on peut avec une
très bonne approximation considérer que P ≈ π(a + b) ; ainsi pour Mercure :
P a + b
≈ = 0, 382962
2π 2
168
8. Formules pour mouvement elliptique
• Par rapport aux planètes, les comètes ont des excentricités élevées. Il faut
donc s’attendre d’une part à ce que la convergence des séries permettant d’at-
teindre le périmètre soit plus lente et d’autre part à ce que les formules de
Ramanujan soient moins précises.
Pour la comète de Halley, les données sont :
a = 17,941 990 ua e = 0,967 277 8 T = 27 759 j
Le tableau ci-après donne une idée de la convergence de la série exponentielle,
et de la précision des deux formules de Ramanujan :
série en e2n P/2π P V
(nombre de termes) ua ua km/s
1 17,942 112,733 7,032
2 13,745 86,364 5,387
--------------- ----------- ----------- ----------
6 12,489 78,468 4,894
7 12,433 78,118 4,873
8 12,395 77,879 4,858
9 12,368 77,709 4,847
10 12,348 77,584 4,839
11 12,333 77,489 4,833
12 12,321 77,416 4,829
--------------- ----------- ----------- ----------
→∞ 12,265 77,064 4,807
Ramanujan
1re formule 12,266 77,073 4,807
2e formule 12,268 77,081 4,808
169
Guide de localisation des astres
ρ dρ
ρ dθ
θ
A
θB
θA
O
Figure 8.2 Aire d’un secteur en coordonnées
polaires.
Application à l’ellipse
L’aire d’un secteur d’ellipse d’origine le foyer F et délimité par les rayons
ρ(θ0) et ρ(θ1) a donc pour expression :
p2 θ1 dθ
S =
2 ∫θ0 (1 + e cos θ)2
e<1
On déduit :
− e sin θ 2 1− e θ
F (θ) = + Arctg tg
(1 − e )(1 + e cos θ) (1 − e )
2 2 3/ 2
1+ e 2
ou encore,
1 ⎡ − e 1 − e 2 sin θ 1 − e θ⎤
F (θ) = ⎢ + 2 Arctg tg ⎥
(1 − e 2 )3 / 2 ⎢⎣ 1 + e cos θ 1 + e 2 ⎥⎦
π
On calcule : F(0) = 0, F (π) =
(1 − e 2 )3 / 2
p2
L’aire du secteur d’ellipse s’écrit : S =
2
[ F (θ1 ) − F (θ0 )]
p2
Pour θ0 = 0, on notera simplement : S(θ) = F (θ)
2
170
8. Formules pour mouvement elliptique
On retiendra finalement :
ab ⎡ − e 1 − e 2 sin θ 1 − e θ⎤
S(θ) = ⎢ + 2 Arctg tg ⎥
2 ⎢⎣ 1 + e cos θ 1 + e 2 ⎦⎥
Angle médian
La figure 8.3 visualise celui-ci.
θμ
A’ A
F
P’
L’angle médian est l’angle θμ tel que l’astre passe autant de temps d’un côté
de l’axe P'P que de l’autre.
Compte tenu de la loi des aires, il s’en suit que :
Aire (AFP) = Aire (A'FP)
D’où :
F(θμ) – F(0) = F(π) – F(θμ) → 2 F(θμ) = F(π)
L’angle médian satisfait donc la relation :
− e 1 − e 2 sin θμ 1 − e θμ π
+ 2 Arctg tg =
1 + e cos θμ 1+ e 2 2
On note que : e → 0 ⇒ θμ → π/2
e → 1 ⇒ θμ → π
La relation peut s’exprimer plus simplement en fonction de l’anomalie excen-
trique ϕμ correspondant à θμ ; en effet, on a établi au chapitre 5 que :
171
Guide de localisation des astres
1 − e 2 sin θ ϕ 1− e θ
sin ϕ = et tg = tg
1 + e cos θ 2 1+ e 2
π
Il vient donc : ϕ μ − e sin ϕ μ =
2
T
Pour cet angle ϕμ on doit vérifier que le temps de parcours t AP = ; c’est
4
bien ce que fournit l’équation de Kepler.
Enfin, l’angle θk tel que l’Astre est situé k pour cent du temps à droite de P'P
(côté périastre) satisfait à :
− e 1 − e 2 sin θ k 1 − e θk
+ 2 Arctg tg = kπ
1 + e cos θ k 1+ e 2
172
8. Formules pour mouvement elliptique
exercices
E8-1 Caractérisation du mouvement elliptique par rayons
et vitesse
Les données du mouvement sont les rayons ρ0 au péricentre et ρπ à l’apo-
centre, ainsi que la vitesse V0 au péricentre.
Calculer :
1. l’excentricité ;
2. le demi-grand axe, le demi-petit axe, la demi-distance focale ;
3. le paramètre, la distance foyer-directrice, la directrice ;
4. la vitesse à l’apocentre et la vitesse en fonction de l’angle polaire ;
5. la période, l’énergie massique ;
6. les constantes de l’accélération et des aires.
solutions
1. Excentricité e
ρπ 1+ e ρ π − ρ0
De = on tire → e =
ρ0 1− e ρ π + ρ0
2. Demi-grand axe a, demi-petit axe b, demi-distance focale c
(Voir chapitre 2, paragraphe 1.2)
ρ0 + ρ π ρ π − ρ0
a = b = ρ0 ρ π c =
2 2
3. Paramètre p, distance foyer-directrice k, directrice d
(Voir chapitre 2, paragraphe 1.2)
p = b2/a , k = b2/c , d = a2/c
ρ0 ρ π ρ0 ρ π 1 (ρ0 + ρ π )2
p = 2 k = 2 d =
ρ0 + ρ π ρ π − ρ0 2 ρ π − ρ0
4. Vitesses, Vπ à l’apocentre , V en fonction de l’angle θ
ρ0
•ρ π Vπ = ρ0 V0 (chapitre 3, paragraphe 3.1) → Vπ = V0
ρπ
1 + e 2 + 2e cos θ
exercices
• V = V0 Vπ
1 − e2
Or :
ρ0 2 ρ2π + ρ20 ρ π ρ0
V0 Vπ = V0 ; 1 + e2 = 2 ; 1 − e2 = 4
ρπ (ρπ + ρ0 )2 (ρπ + ρ0 )2
173
Guide de localisation des astres
V0
d’où → V = 1 + ρ20 / ρ2π + (1 − ρ20 / ρ2π ) cos θ
2
5. Période T, énergie massique E/m
(Voir chapitre 8, paragraphe 1)
ρ0 + ρ π ρ0 + ρ π ρπ
• T = π → T = π
V0 Vπ V0 ρ0
1 E 1 ρ0 2
• E = − mV0 Vπ → = − V0
2 m 2 ρπ
6. Constantes μ de l’accélération et C des aires
μ 1 ρ
• De V0 Vπ = on déduit : → μ = V02 (ρ0 + ρ π ) 0
a 2 ρπ
ab
• C = 2π
T
Compte tenu des expressions de a, b, et T obtenues ci-dessus, on a :
→ C = ρ0 V0
Remarque
À partir de l’équation polaire, on montre facilement que :
1 1 2
+ =
ρ(θ) ρ(θ + π ) p
1 1 2
et que + 2 = 2 (1 + e 2 cos 2 θ)
ρ (θ) ρ (θ + π )
2
p
ce qui pour θ = 0 conduit à :
1 1 2 1 1 1 + e2
+ = et + 2 = 2
ρ0 ρπ p ρ0
2
ρπ p2
d’où on en déduit p, puis e.
et Vπ à l’apocentre.
Calculer :
1. l’excentricité ;
2. le demi-grand axe, le demi-petit axe, la demi-distance focale ;
3. le paramètre, la distance foyer-directrice, la directrice ;
174
8. Formules pour mouvement elliptique
solutions
Préliminaires
L’aire balayée par le rayon polaire durant une période est égale à la surface de
l’ellipse, la loi des aires permet d’écrire :
1
πab = CT
2
C est la première constante dynamique, elle satisfait aux égalités :
2π
C = ab = ρV sin α = ρ0 V0 = ρπ Vπ
T
avec :
p p
ρ0 = , ρπ =
1+ e 1− e
1. Excentricité e
V0 ρ 1+ e
De = π =
Vπ ρ0 1− e
On tire :
V0 − Vπ 2Vπ 2V0
e = 1− e = 1+ e =
V0 + Vπ V0 + Vπ V0 + Vπ
On tire :
2π 2
ρ0 V0 = a (1 − e )V0 = a 1 − e2
T
d’où
exercices
T 1− e T
a = V0 , b = (1 − e )V0 , c = ae
2π 1 + e 2π
Il vient finalement :
T T V0 Vπ T V0 − Vπ
a = V0 Vπ b = c = V0 Vπ
2π π V0 + Vπ 2 π V0 + Vπ
175
Guide de localisation des astres
1 + e 2 + 2e cos θ
• On a vu au chapitre 3 que : V = V0 Vπ , et on connaît
1 − e2
l’expression de l’excentricité en fonction des vitesses V0 et Vπ ; toutefois,
ρ V
notant que 0 = π , on reprendra directement la formule de la vitesse
ρπ V0
établie à l’exercice précédent :
V0
V = 1 + Vπ2 / V02 + (1 − Vπ2 / V02 ) cos θ
2
1
soit encore : V = V02 + Vπ2 + (V02 − Vπ2 ) cos θ
2
5. Constante C, énergie massique E/m et constante μ
T (V0 Vπ )3 / 2
• Partant de C = ρ0V0 , on a immédiatement : C =
π V0 + Vπ
exercices
E 1
• On a vu au chapitre 8 que : = − V0 Vπ
m 2
a3 e T
• μ = 4π2 (3 loi de Kepler) → μ = (V0 Vπ )3 / 2
T2 2π
176
8. Formules pour mouvement elliptique
α V
I
P
ρ
Soleil
Figure 8.4
177
Guide de localisation des astres
solutions
1. Énergie massique E/m
L’énergie mécanique de la comète de masse m est la somme des énergies ciné-
tique et potentielle, d’où :
E 1 μ
= V2 − s → E/m = −24,722 115 106 J.kg-1
m 2 ρ
La condition pour que la trajectoire soit elliptique(*) est que l’énergie soit
négative :
1 2 μs 1
V − < 0 , soit ρV 2 < μ s ; nous sommes donc dans ce cas.
2 ρ 2
ρV 2
On pose K = et on en déduit que 0 < K < 2 ;
μs
1
On notera qu’il faut satisfaire à sin 2 α ≤ , ce qui est toujours
K (2 − K )
vérifié dans l’intervalle {0, 2}.
(*) la trajectoire ne peut être circulaire dans la mesure où α est différent de π/2.
178
8. Formules pour mouvement elliptique
C ρV sin α
VP = = → VP = 54,521 702 km.s-1
ρP ρP
5. Période de rotation T
⎛ m ⎞ 32 1
T = 2 πμ s ⎜ = 2 πμ s → T = 27 759 j (de 86 400 s)
⎝2E⎠⎟ 3
(VP VA ) 2
179
Guide de localisation des astres
180
8. Formules pour mouvement elliptique
solutions
1. Éléments de la trajectoire elliptique
a) Expression de l’excentricité e et du paramètre p
p
L’équation polaire en H est : ρH =
1 + e cos θH
p
L’équation polaire en T est : ρT = avec
exercices
1 + e cos θT
θT = θH + Δθ (1) où Δθ = 44°,544 737
On en déduit :
ρT − ρH cos θH − cos θT
e = (2) p = ρHρT (3)
ρH cos θH − ρT cos θT ρH cos θH − ρT cos θT
181
Guide de localisation des astres
b) Expression de la période TG
a3 p
Partant de TG2 = 4 π 2 avec a = on écrit :
μS 1 − e2
2π ⎛ p ⎞
3/ 2
TG = (4)
μS ⎝ 1 − e2 ⎠
2. Équation de Kepler et durée des trajets
a) Durée Δt du trajet de l’aphélie terrestre à la comète
La date en T est le 4 juillet 1985 à 17 h 48 min 53 s soit 0,742 275 jour.
La date en H est le 10 mars 1986 à 8 h 50 min 55 s soit 0,368 692 jour.
On retiendra comme durée du trajet : Δt = 248,626 j
b) Durée tT – tH du trajet de la comète à l’aphélie terrestre
L’équation de Kepler donne la durée t du trajet à partir du péricentre :
T
t = (ϕ − e sin ϕ ) , ϕ étant l’anomalie excentrique.
2π
La sonde Giotto a une trajectoire elliptique de période TG , et pour aller de H
à T elle mettrait le temps tT – tH tel que :
TG
t T − t H = [ (ϕ T − e sin ϕ T ) − (ϕ H − e sin ϕ H ) ] (5)
2π
La durée Δt du trajet est liée au temps tT – tH par la relation :
Δt = TG − ( t T − t H ) (6)
Afin d’exprimer la relation (5) uniquement en fonction de l’anomalie vraie θH,
il nous reste à donner les expressions convenables de l’anomalie excentrique.
c) Anomalies excentriques ϕH et ϕT
• Anomalies excentriques et anomalies vraies en H et T, sont liées par les
relations :
1 − e 2 sin θH 1 − e 2 sin θT
sin ϕ H = (7) sin ϕ T = (7’)
1 + e cos θH 1 + e cos θT
• Comme vu à l’exercice E5-1, la détermination de la fonction trigonomé-
trique inverse va dépendre du couple :
exercices
182
8. Formules pour mouvement elliptique
183
Guide de localisation des astres
184
8. Formules pour mouvement elliptique
Remarques
1– La vitesse relative obtenue est sensiblement supérieure à celle corres-
pondant à la mission réelle où le plan de la trajectoire était légèrement
incliné sous l’écliptique, conduisant à un rendez-vous le 14 mars 1986,
c’est-à-dire 4 jours plus tard.
Soit Nδ ce nouveau point ; se reportant à l’exercice E5-1 on en déduit
qu’à partir du périhélie de la comète :
tNδ = 32,91 jours
d’où
ϕNδ = 11°,031 954 θNδ = 73°,650 501 soit
θNδ − θNd = 5°,497 021
On déduit alors :
ρNδ = 0,907 809 ua
VNδ = 43,646 km/s avec αNδ = 53°,888 VρNδ = 25,723 km/s
VθNδ = 35,260 km/s
Le problème peut être repris connaissant l’inclinaison de la trajectoire de
Giotto et après avoir déterminé le nouvel angle Δθ...
2– Anticipant sur les prochains chapitres, et en particulier sur le chapitre 11,
nous avons représenté sur la figure 8.5 la trajectoire de la sonde Giotto en y
faisant figurer quelques points remarquables et les vitesses en ceux-ci.
Voici quelques éléments du cheminement :
Détermination de la date du passage de notre planète à l’apogée
en 1985.
Pour dégrossir, le 4 juillet à 0 h correspond à 5294,5 jours avant le
1er janvier 2000 à 12 h TU.
En fait, on atteint le temps τ pour une anomalie moyenne M telle que :
M = − (360 14 + 180) = − 5220°
On déduit ainsi le nombre de jours avant la référence J2000, soit
5293,757 725 264 jours !
D’où la date du : 4 juillet 1985 à 17 h 48 min 53 s.
À cette date, pour la Terre :
a = 1,000 001 018 ua e = 0,016 714 724
Le rayon à l’aphélie est : ρA = 1,016 716 ua
L’argument du périhélie ωT = 102°,688 147 448 = 102° 41’ 17’’,33
exercices
185
Guide de localisation des astres
I
Nœud
descendant Halley
VH
H(Nd)
Aphélie
Terre T
VP
Δθ
VT
θH
Périhélie
Aphélie S Giotto
Giotto A P
VA ωT
Périhélie
ΩH Terre T’
γ
Nœud
ascendant Halley
Na
186
Chapitre 9
Repérage des astres
par rapport à la Terre
Nous nous proposons maintenant, à Un premier exercice nous fait calculer
partir de notre planète Terre, de situer la plus courte distance entre deux
un astre : autre planète du système points de coordonnées géographi-
solaire, satellite, étoile… Après avoir ques connues. Le deuxième exercice
défini les coordonnées terrestres du montre l’intérêt d’une base de lance-
point d’observation, nous préciserons ment équatoriale : on bénéficie
le repère local dans lequel par exemple pleinement de la vitesse d’entraîne-
une monture alt-azimutale peut être ment de la Terre lors d’un tir vers l’Est,
positionnée sur un objet céleste donné. et il n’est plus besoin de faire pivoter le
Un autre repère lié à la Terre sera plan de l’orbite de transfert d’un satel-
ensuite défini par rapport à l’équateur lite géostationnaire. Nous complétons
et au méridien de Greenwich, l’angle ainsi les exercices E7-1 et E7-2. Le troi-
horaire d’un astre y sera affecté de sième exercice concerne également les
la rotation de notre planète que les satellites géostationnaires ; il s’agit
montures équatoriales des instru- maintenant d’orienter une antenne
ments d’observation astronomique parabolique dans la direction d’un
savent compenser. À partir du triangle satellite de télécommunication grâce à
de position, à l’aide des relations de la connaissance de ses angles d’azimut
Gauss, nous établirons les formules de et de site.
passage entre coordonnées horaires Nous quitterons enfin les satellites
et locales. artificiels pour nous intéresser aux posi-
tions du Soleil à différents moments
de la journée, en précisant l’heure de
passage au méridien et les instants des
crépuscules civil et nautique.
Guide de localisation des astres
1 Coordonnées géographiques
188
9. Repérage des astres par rapport à la Terre
méridien du lieu
PN
90°
parallèle du lieu
60°
P L
méridien origine
(Greenwich)
C 30°
lat.
équateur
90° long.
-30°
60° G O
30°
-30°
-60°
PS -90°
189
Guide de localisation des astres
Remarques
1– Le périmètre de l’équateur étant de 40 075 km, on en déduit que :
• 1° correspond à 111,32 km à l’équateur,
• 1’, soit 1/60e de degré, correspond à 1 855,3 mètres.
Le mille marin ou mille nautique ou simplement mille, est une unité de
distance correspondant à un arc de méridien de une minute d’angle. Sur la
Terre, le méridien n’est pas rigoureusement un cercle ; on a choisi comme
valeur de cette unité celle correspondant à la mesure à une latitude de 45°,
soit 1 852 m. On retient donc :
1Nm = 1 852 m
2– L’angle complémentaire de la latitude est la colatitude qui correspond
à l’arc PNP. Elle a donc pour valeur :
(CPN, CP) = 90° − L
2 Coordonnées locales
Les coordonnées locales sont aussi dénommées coordonnées horizontales.
On se propose de préciser la position d’un astre dans un repère terrestre local.
190
9. Repérage des astres par rapport à la Terre
Nord
direction de Zénith
l’astre
PN
h horizon du lieu
P
Z
Ouest C Est
projection de
l’astre
méridien du lieu
PS
191
Guide de localisation des astres
direction de
l’astre Zénith
90°-h
horizon du
projection de h lieu
l’astre CA P
Pied de l’astre PA
c dz
P’
C
vertical de
l’astre
192
9. Repérage des astres par rapport à la Terre
direction de
Z
l’astre
axe vertical
axe horizontal
horizon du lieu
193
Guide de localisation des astres
3 Coordonnées horaires
3.2 Déclinaison
La droite menée du centre de la Terre en direction de l’astre coupe la sphère
en PA.
– L’angle entre cette droite et sa projection sur le plan équatorial est la décli-
naison D. Tout comme la latitude, la déclinaison est comptée positivement du
plan équatorial vers le pôle Nord. Elle est aussi notée δ.
– L’angle complémentaire est l’angle entre les directions du pôle Nord et de
l’astre, considérées à partir de C. Il est dans le plan contenant l’axe des pôles et
le point PA qui est perpendiculaire au plan équatorial, l’intersection étant ainsi
la projection de la droite CPA.
Par conséquent, cet angle qui correspond à l’arc PNPA vaut 90° − D ; il porte le
nom de distance polaire. On note :
(CPN , CPA) = 90° − D
194
9. Repérage des astres par rapport à la Terre
N Soleil
tropique du cancer
PG
D
équateur
tropique du capricorne
On constate que pour L > + 66° 33’ 29’’ c’est le jour permanent ;
L < − 66° 33’ 29’’ c’est la nuit permanente.
195
Guide de localisation des astres
On notera cet angle AHG (en anglais Greenwich Hour Angle) comme le font
certains auteurs, toutefois la notation en vigueur dans les éphémérides est
souvent :
– AHao pour l’angle horaire des planètes,
– AHvo pour l’angle horaire du Soleil.
– L’angle entre le méridien du lieu et le cercle horaire de l’astre est l’angle
horaire local compté du méridien de l’observateur vers celui de l’astre avec
les mêmes conventions que pour AHG.
On note cet angle AHL (en anglais Local Hour Angle), sachant que dans les
éphémérides,
– AHag correspond à l’angle horaire local des planètes,
– AHvg correspond à l’angle horaire local du Soleil.
On constate que l’angle entre le méridien origine et le méridien local correspond
à la longitude de l’observateur, elle aussi comptée positivement vers l’Ouest.
La relation liant angles horaires et longitude est par conséquent :
AHG = AHL + G
AHG se déduisant des éphémérides, on utilise plutôt :
AHL = AHG – G
Déclinaison et angles horaires sont représentés à la figure 9.6 :
direction de PN
l’astre
PA
C
D
AHG O
AHL
PS
Figure 9.6 Coordonnées horaires d’un astre : D ou δ est la déclinaison ; AHG et AHL
sont les angles horaires origine et local.
196
9. Repérage des astres par rapport à la Terre
197
Guide de localisation des astres
axe horaire
(polaire)
direction de PN
L
l’astre
axe des
déclinaisons
PN
90°- L
90°- D
P
90°- h
PA
198
9. Repérage des astres par rapport à la Terre
z
A
b
c
C
a
B O y
C’
199
Guide de localisation des astres
ˆ
cos a = cos b cos c + sin b sin c cos A
ˆ
ˆ = sin b sin A
sin a sin B
ˆ = cos b sin c − sin b cos c cos A
sin a cos B ˆ
(*) Carl Friedrich Gauss est né à Brunswick, en Allemagne, le 30 avril 1777. Très tôt, ses talents
de mathématicien furent connus et grâce à Ferdinand, duc de Brunswick, il put poursuivre des
études et intégrera l’Université de Göttingen en 1795 où il restera comme professeur jusqu’à sa
mort le 23 février 1855.
Il eut à effectuer la topographie de plusieurs régions d’Allemagne ; à cette occasion, il donna
un support théorique aux erreurs aléatoires et fit progresser la géométrie non euclidienne. Celle
des espaces elliptiques ou géométrie sphérique trouve une application en cartographie grâce
aux relations entre les éléments du triangle sphérique.
Le Français Adrien-Marie Legendre (1752-1833) s’est également intéressé à la géodésie et à la
trigonométrie sphérique.
200
9. Repérage des astres par rapport à la Terre
Détermination de la hauteur h
Les arcs PNPA et PNP étant connus, de même que l’angle P̂N entre ceux-ci, la
substitution donne :
cos(90° − h) = cos(90° − L )cos(90° − D)
+ sin(90° − L )sin(90° − D)cos Pˆ N
soit :
sin h = sin L sin D + cos L cos D cos Pˆ N
avec :
P̂N = AHL (astre à l’Ouest) ou 360° − AHL (astre à l’Est)
mais dans les deux cas :
cos Pˆ N = cos AHL
Ainsi, la formule donnant la hauteur est :
sin h = sin L sin D + cos L cos D cos AHL
La hauteur de l’astre dépend directement de sa déclinaison et de la latitude
d’observation ; et elle dépend de la longitude par l’intermédiaire de l’angle
horaire local, puisque AHL = AHG − G.
Astre
N h D-L
D L
201
Guide de localisation des astres
Détermination de l’azimut Z
• En utilisant la relation donnant l’angle en fonction des trois côtés :
On considère l’angle P̂ compris entre les arcs PPN et PPA ; appliquant la formule
donnant l’angle  du triangle sphérique, on écrit :
cos(90° − D) − cos(90° − L)cos(90° − h)
cos Pˆ =
sin(90° − L )sin(90° − h)
soit :
sin D − sin L sin h
cos Pˆ =
cos L cos h
si 0° ≤ AHL < 180° Z = 360° − Pˆ (astre à l’Ouest)
si 180° ≤ AHL < 360° Z = Pˆ (astre à l’Est)
202
9. Repérage des astres par rapport à la Terre
Point astronomique
La formule de la hauteur dans laquelle on explicite le temps, s’écrit :
sin h = sin L sin D(t) + cos L cos D(t) cos [AHG(t) – G]
Pour un astre donné, à une heure donnée, les éphémérides permettent de calculer
la déclinaison et l’angle horaire origine. On en déduit une méthode, basée sur
la mesure de la hauteur d’un astre, pour faire le point astronomique ; c’est-à-
dire déterminer la latitude et la longitude du lieu d’observation P :
À l’instant donné t1, on mesure la hauteur h1 de l’astre sélectionné ; puis à
l’instant t2 on mesure la hauteur h2 du même astre ou d’un deuxième. Sur une
plateforme mobile (bateau ou aéronef) on utilise à cette fin un sextant ; sur
terre on utilise un théodolite.
203
Guide de localisation des astres
Quelles que soient les positions relatives des méridiens de l’astre et du lieu,
on obtient finalement :
sin D = sin h sin L + cos h cos L cos Z
cos D sin AHL = − cos h sin Z
cos D cos AHL = sin h cos L − cos h sin L cos Z
et la tangente de l’angle horaire local a pour expression :
− cos h sin Z
tg AHL =
sin h cos L − cos h sin L cos Z
204
9. Repérage des astres par rapport à la Terre
exercices
E9-1 Orthodromie
On se propose de calculer la plus courte distance, ou géodésique, entre deux
points P1 et P2 sur notre planète de centre C.
1. Exprimer les coordonnées cartésiennes x, y, z d’un point P dont les coor-
données sphériques sont (R, L, G). On positionnera sur une figure ce point
dans le référentiel choisi d’origine C.
2. Donner l’expression de l’angle θ entre les vecteurs CP1(R1, L1, G1) et
CP2(R2, L2, G2).
3. Quelle est la longueur L de l’arc de grand cercle reliant le centre Guya-
nais de Kourou à Toulouse Blagnac ?
Données : Kourou, L2 = 5° 14’ N G2 =52° 45’ W
Toulouse, L1 = 43° 38’ N G1 = 1° 22’ E
Rayon terrestre = 6378 km
solutions
1. Coordonnées x, y, z de P
C
L
j
i
x
G
p y
205
Guide de localisation des astres
b) Dans quel plan doit se situer la vitesse d’injection V si l’on veut que le
point M corresponde au périgée ?
LLLLI
c) Quelles sont les coordonnées du vecteur position CM de module R et du
I
vecteur vitesse V d’angle a avec la direction de l’Est ?
d) Préciser le plan de l’orbite de transfert. À quelle relation, fonction des
angles a et L, satisfait son inclinaison i sur l’équateur ? Commenter.
206
9. Repérage des astres par rapport à la Terre
solutions
1. Plan de l’orbite de transfert
a) Voir figure 9.12a.
b) Si l’on veut que le point d’injection
I LLLI corresponde au périgée de l’orbite de
transfert, il faut satisfaire à : V ⊥I CM
→ V estLLLdans
I le plan horizontal contenant M.
exercices
207
Guide de localisation des astres
z
V
M a
Vsina
Est
L
C M
y L
G
x C
i≥L
L’inclinaison de l’orbite de transfert ne peut donc être inférieure à la latitude
du lieu de lancement. Le cas limite où l’inclinaison est égale à la latitude est
obtenu pour un tir vers l’Est (a = 0) donc parallèle au plan équatorial.
Pour a = ± π/2, i = ± π/2 → on obtient une orbite polaire.
208
9. Repérage des astres par rapport à la Terre
En conclusion, seule une base équatoriale permet d’obtenir toutes les inclinai-
sons désirées.
2. Gain de vitesse
a) La vitesse tangentielle de rotation de la Terre à l’équateur est :
2 πR T
VT = R T ω T = → VT = 465 m.s-1
JSid
La vitesse tangentielle à un petit cercle de latitude L est : VL = 465 cos L m.s-1
Kourou : L = 5°,23 → VL = 463 m.s-1
Cap Canaveral : L = 28°,45 → VL = 409 m.s-1
Baïkonour : L = 46°,25 → VL = 322 m.s-1
b) L’énergie cinétique massique au sol à la latitude L est :
EC 1
= VL2 = 0,108 cos 2 L MJ / kg
m 2
E GM T 1 E
= − + VL2 → = (− 62, 495 + 0,108 cos 2 L ) MJ / kg
m RT 2 m
3. Modification de l’orbite
a) Il faut que l’orbite de transfert puisse pivoter autour de l’axe des nœuds,
intersection du plan équatorial avec cette orbite (voir chapitre 11). Donc le
point P correspondant au périgée doit être dans le plan équatorial ; en consé-
quence, la ligne des apsides est confondue avec la ligne des nœuds.
b) La vitesse VA à l’apogée étant plus faible qu’au périgée, l’impulsion de
vitesse ΔV à appliquer sera également plus faible et impliquera en consé-
quence une moindre consommation de propergol.
c)
δi
DV VA
i
A
V2
exercices
209
Guide de localisation des astres
Plan de l’orbite
de transfert
VA A
DV
VP
V2 P
GTO
GEO
Plan de l’orbite
exercices
géostationnaire
210
9. Repérage des astres par rapport à la Terre
Remarque
Le lancement d’un satellite de la base de Kourou en Guyane vers l’orbite
géostationnaire comprend :
– la phase de montée verticale hors de l’atmosphère qui dure 3 minutes et
est assurée par le premier étage de la fusée ;
– la phase intermédiaire qui est assurée par le deuxième étage et au cours
de laquelle montée et mise en vitesse se poursuivent pendant 2 minutes,
mais le rôle de cette phase est également d’infléchir la trajectoire qui doit
être horizontale en fin de combustion ;
– la phase d’accélération horizontale qui dure environ 12 minutes de façon
à ce que le satellite atteigne la vitesse de 10,239 km/s au passage à la verti-
cale de l’équateur ; ajoutons qu’avant sa séparation à l’issue de cette phase,
le satellite est orienté et stabilisé par mise en rotation sur lui-même.
À l’apogée, c’est-à-dire après 5 h 15 min 32 s à compter du passage au
périgée, le satellite ayantLLL
été
I correctement orienté, son moteur d’apogée
fournira l’accroissement ΔV de vitesse requis.
211
Guide de localisation des astres
Nord
F P0 P2 P4 satellite S
W P1
P3
solutions
1. Identification des angles D, L, G, A et S
D̂ = (FW, FS)
exercices
LI les angles θ, D, G et L
2. Relation entre
(PP4) ⊥ V
LI → FP4 cos θ = R
(PP3) ⊥ V → FP3 cos L = R
212
9. Repérage des astres par rapport à la Terre
nom position
exercices
azimut A site S θ
du satellite D
EUTELSAT II F4 10°,0 E 10°,176 ou 10° 11’ E 33°,727 ou 33° 44’ N 49°,046 ou 49° 03’ N
ASTRA 1 19°,2 E 22°,011 ou 22° 01’ E 31°,886 ou 31° 53’ N 50°,735 ou 50° 44’ N
TÜRKSAT 1C 42°,0 E 47°,871 ou 47° 52’ E 22°,574 ou 22° 34’ N 59°,397 ou 59° 24’ N
INTELSAT 804 64°,0 E 67°,981 ou 67° 59’ E 09°,719 ou 09° 43’ N 71°,707 ou 71° 42’ N
EXPRESS 6A 80°,0 E 80°,688 ou 80° 41’ E -0°,579 ou 00° 35’ S 81°,880 ou 81° 53’ N
213
Guide de localisation des astres
Remarques
Le satellite russe EXPRESS 6A ne pourra être vu puisqu’il se trouve juste
sous l’horizon local. Il faut en effet remplir la condition :
R 1 R
cos θ ≥ → cos(D − G) ≥ = 0, 228 708
R+H cos L R + H
soit :
D − G < 76°,779
d’où
D < 79°,111 E
On notera que si le satellite est trop bas sur l’horizon, la qualité du signal
reçu sera affectée par le bruit rayonné par la Terre.
solutions
1. Variations horaires V de l’angle horaire et d de la déclinaison
On transforme les données sexagésimales en données décimales :
– le 1er juillet : AHG = 179°,0917 D = 23°,1
– le 2 juillet : AHG = 179° D = 23°,0283
214
9. Repérage des astres par rapport à la Terre
215
Guide de localisation des astres
Remarques
1– Les heures ont été calculées en faisant abstraction de la réfraction
d’angle r0 = 34’ ; tenant compte de celle-ci, le cosinus de l’angle horaire
local a comme nouvelle expression :
− sin r0
cos AHL = − tg L tgD
cos L cos D
On calcule alors :
exercices
lever coucher
terme correctif = − 0,016 335 − 0,016 329
cos AHL = − 0,503 905 − 0,502 767
AHL = 239°,7413 120°,1832
t= 3 h,8885 19 h,9220
D’où : tl = 3 h 53 min 18 s tc = 19 h 55 min 19 s
216
9. Repérage des astres par rapport à la Terre
Remarque
Si on tient compte des deux corrections opérées ci-dessus, h n’est plus nul :
la hauteur est en effet légèrement négative puisqu’elle vaut −(r0 + 0,5s).
Utilisant la relation donnant tg Z pour calculer les azimuts au lever et au
exercices
coucher, il vient :
Au lever tg Zl = + 1,290 658 → Zl = + 52°,2315
par conséquent : Zl = 52° 13’ 53’’ et al = 232° 13’ 53’’ (− 127° 46’ 7’’)
Au coucher tg Zc = − 1,294 615 → Zc = − 52°,3164 = 307°,6836
par conséquent : Zc = 307° 41’ 1’’ et ac = 127° 41’ 1’’
217
Guide de localisation des astres
h tl Zl tc Zc
− 50’ 3 h 51 min 16 s 52° 13’ 51” 19 h 57 min 22 s 307° 41’ 03”
− 6° 3 h 09 min 03 s 44° 00’ 03” 20 h 39 min 29 s 315° 53’ 37”
− 12° 2 h 09 min 28 s 31° 32’ 01” 21 h 38 min 50 s 328° 18’ 27”
− 18° 0 h 07 min 52 s 3° 17’ 01” 23 h 35 min 00 s 355° 14’ 52”
Les précisions des résultats fournis au cours de cet exercice sont un peu
illusoires ; elles sont données dans le but de mieux suivre les étapes des
calculs effectués sans approximation à l’aide d’un tableur.
exercices
218
Chapitre 10
Coordonnées équatoriales
et écliptiques
1 Voûte céleste
1.1 Constellations
Lorsque de l’hémisphère Nord on observe le ciel nocturne au voisinage de
l’Étoile polaire, les constellations apparaissent tourner autour de celle-ci d’Est
en Ouest, comme illustré sur la figure 10.1 :
Nord
Petite Ourse
Ouest Est
Polaire
PN
P
Observateur
220
10. Coordonnées équatoriales et écliptiques
Pôle Nord
céleste Solstice
d’été
Soleil
Équinoxe Gem
d’automne Cnc
Tau
Leo
Vir Ari
Lib Psc
Équateur
Sco Aqr céleste
Cap Équinoxe
de printemps
Sgr
221
Guide de localisation des astres
1.2 Magnitude
Une constellation représente donc une configuration apparente d’étoiles, telles
qu’elles nous semblent positionnées sur la voûte céleste, donc abstraction faite
du paramètre distance. Une autre grandeur a toutefois été introduite : c’est la
magnitude apparente m, nombre caractérisant l’éclat, c’est-à-dire l’inten-
sité, d’une étoile.
Plus un astre est brillant, plus sa magnitude est faible. Dans l’échelle de
Pogson, pour un rapport d’intensité de 1 à 100, les magnitudes vont de 0 à 5.
Le rapport q entre deux magnitudes consécutives étant constant, il vient :
q5 = 100 d’où q = 100,4 ≈ 2,512
Les étoiles les plus brillantes peuvent avoir des magnitudes négatives, comme
on peut le constater sur le tableau en 2.5. La notion de magnitude est égale-
ment appliquée aux planètes : ainsi, celle de Mars varie entre – 2,6 (position la
plus proche de la Terre) et + 2,0 ; celle d’Uranus entre + 5,65 et + 6,07.
Notons que pour bien comparer les luminosités des étoiles, les magni-
tudes doivent être recalculées en considérant des éloignements identiques :
la magnitude absolue M est ainsi calculée en supposant toutes les étoiles à
10 pc (parsecs) de la Terre.
222
10. Coordonnées équatoriales et écliptiques
2 Coordonnées équatoriales
223
Guide de localisation des astres
Le point gamma, aussi dénommé point vernal, est de fait dans la constellation
des Poissons.
Deux coordonnées angulaires permettent de positionner tout objet céleste dans
un repère lié au référentiel de la sphère céleste. Le premier angle sera défini
par rapport au plan de l’équateur céleste (confondu avec l’équateur terrestre) ;
le deuxième angle sera mesuré dans le plan précédant à partir de l’axe des
équinoxes.
La figure 10.3 donne les coordonnées équatoriales de deux étoiles dans un tel
repère :
Pôle Nord
Méridien céleste de Méridien céleste du
céleste
l’étoile E1 point vernal
Hémisphère E2
boréal
γ’
E1
C δ
Équateur D
céleste
γ
e1 α e2
AV
Hémisphère
austral
2.2 Déclinaison
Elle correspond à la latitude d’un point P de la Terre en coordonnées géogra-
phiques.
C’est l’angle D de la ligne de visée de l’étoile, menée du centre C de la sphère,
avec l’équateur céleste : la définition est identique à celle donnée pour les
coordonnées horaires au chapitre 9. La déclinaison est positive, de 0° à + 90°,
vers le pôle Nord céleste ; négative, de 0° à – 90°, vers le pôle Sud céleste.
La déclinaison est aussi notée δ par les astronomes.
224
10. Coordonnées équatoriales et écliptiques
225
Guide de localisation des astres
226
10. Coordonnées équatoriales et écliptiques
L’angle horaire origine d’un astre est aussi noté AHao, et celui d’une étoile
AH*o. Pour l’angle horaire local, les notations sont respectivement AHag
et AH*g.
• Des formules ci-dessus, on déduit par exemple l’heure de passage tm d’une
étoile au méridien.
Au méridien, AHLe(tm) = 0
D’où on déduit : AHGγ(tm) = AHGγ(0) + Vγ tm = G – AVe
G − AVe − AHGγ (0)
Ainsi : tm =
Vγ
Par exemple, le 1er juillet 2004 à Paris, pour Arcturus, les données sont :
G = – 2° 20’ 15” α = 14 h 15 min 52 s T0 = 18 h 37 min 32 s
et on trouve tm = 19 h 25 min 47 s
(Résultat en cohérence avec ceux des exercices E10-3 et E10-4).
227
Guide de localisation des astres
Ces étoiles sont reportées sur la carte du ciel de la figure 10.4. Le Soleil y
est également représenté en certains points de l’écliptique. Ses positions sont
données à l’équinoxe de printemps (EP), au solstice d’été (SE), à l’équinoxe
d’automine (EA) et au solstice d’hiver (SH) ; ainsi que le 21 de chaque mois
intermédiaire pour l’années 2004.
12h
Spica
10h 14h
EA
Régulus
Arcturus
16h
8h
Antarès
Procyon
Pollux
Sirius
Castor
4h 20h
EP 22h
2h
Fomalhaut
24h
Figure 10.4 Position sur la carte céleste des étoiles les plus brillantes et du Soleil
évoluant sur l’écliptique.
3 Coordonnées écliptiques
228
10. Coordonnées équatoriales et écliptiques
Pôle Nord
écliptique céleste
PN
QN
P
b
γ’
δ
C
céleste
Équateur l
écliptique α
γ
PS QS
céleste écliptique
Pôle Sud
229
Guide de localisation des astres
PN
ε 90°- δ
QN 90°- b
Angles au sommet :
On notera en préliminaire que l’axe γ ' γ des équinoxes est perpendiculaire aux
axes des pôles, d’où :
π
P̂N = + α (angle entre le méridien céleste passant par QN
2
et celui passant par P).
π
Q̂ N = − l (angle entre les demi-cercles QNPNQS et QNPQS).
2
230
10. Coordonnées équatoriales et écliptiques
Cas du Soleil
Sa latitude écliptique bS étant nulle, sa déclinaison et son ascension droite
s’expriment simplement en fonction de sa longitude écliptique lS :
sin δ S = sin ε sin l S
tg α S = cos ε tg l S
D’autre part, puisque sin bS = 0, on obtient directement la relation entre
coordonnées équatoriales :
tg δ S = tg ε sin α S
4 Temps
231
Guide de localisation des astres
Équation du temps
La correction à apporter au temps vrai pour atteindre le temps moyen est
appelée équation du temps E. Elle satisfait aux égalités :
E = tm – tv = Hm – H = α− αm
Notons que pour d’autres pays, l’équation du temps peut correspondre à la
différence tv – tm .
Nombres d’expressions plus ou moins précises de cette équation font appa-
raître explicitement deux termes principaux :
• l’équation du centre C de période un an, due à l’excentricité e de l’orbite
terrestre ;
• la réduction à l’équateur R de période six mois, due à l’obliquité ε de l’éclip-
tique sur l’équateur.
Fondamentalement, l’équation du temps correspond à la partie non linéaire
de l’ascension droite α du Soleil. Pour cet astre, on a établi au paragraphe
précédent que :
tg α = cos ε tg l
l étant la longitude écliptique du Soleil égale, à une constante près, à l’ano-
malie vraie θ de son mouvement elliptique apparent autour de la Terre :
l = ϖ + 180° + θ
Nous verrons au chapitre suivant que la constante ϖ est la longitude du péri-
centre de la Terre ; ainsi, au premier janvier 2000 :
ϖ = 102° 56’ 14”,45
D’autre part, nous savons exprimer θ en fonction de l’anomalie moyenne M
de notre planète ; cette anomalie a, au premier ordre en t, pour expression :
232
10. Coordonnées équatoriales et écliptiques
Temps civil
Il est midi vrai lorsque le Soleil passe au méridien du lieu. À cet instant :
tv = 0 et tm = E.
Il est midi moyen lorsque tm = 0, soit : tv = −E ; par exemple, en février il
faudra attendre plusieurs minutes avant que le Soleil ne soit au méridien.
Il fut décidé de décaler de 12 heures le temps moyen afin d’éviter le désa-
grément du changement de date en milieu de journée induit par la périodicité
de 24 heures. Le temps civil tc est donc défini comme suit :
tc = tm + 12
Le temps civil correspondant au passage physique du Soleil au méridien de
Greenwich est souvent appelé T.Pass dans les éphémérides destinées aux
navigateurs.
233
Guide de localisation des astres
Temps légal
Des nations appliquent certains décalages en nombres entiers d’heures par
rapport au temps universel. En France, l’heure légale tl satisfait à :
tl = UT + 1 pour l’heure d’hiver ;
tl = UT + 2 pour l’heure d’été.
En 2004, l’heure d’été débute le dimanche 28 mars à 2 h du matin et s’achève
le dimanche 31 octobre à 3 h du matin. En 2005 les jours du changement
tombent respectivement le 27 mars et le 30 octobre (dernier dimanche du mois
à la même heure).
234
10. Coordonnées équatoriales et écliptiques
235
Guide de localisation des astres
L’heure UTC est déduite des équations ci-dessus depuis le 1er janvier 1972, et
le temps universel coordonné est l’échelle de temps légale depuis 1978. Il est
diffusé par le réseau des satellites GPS et est accessible en permanence avec
une précision de 1 μs. L’UTC, aux décalages horaires près, est également
transmis par horloges parlantes et stations de radiodiffusion, et des émetteurs
peuvent automatiquement recaler les horloges munies de récepteurs appro-
priés.
Si on ne cherche pas une précision meilleure que la seconde, il y a alors équi-
valence entre UTC et UT1.
D’autre part, le lien avec le temps des éphémérides s’exprime par l’égalité :
TE = UTC + n + 32,184 s
Pour l’année 2000, n = 32.
236
10. Coordonnées équatoriales et écliptiques
ΔT = UT + UT (235,9 / 86164,1)
soit, plus rigoureusement :
ΔT = UT 1,002 737 91
• On déduit le temps sidéral à Greenwich à l’instant donné,
soit T0 + ΔT.
• Le temps sidéral à l’instant donné tl et au lieu considéré de longitude G
est finalement :
T = T0 + ΔT − G
237
Guide de localisation des astres
(*) Voir par exemple : Graham Woan, The Cambridge Handbook of Physics Formulas,
Cambridge University Press 2000.
238
10. Coordonnées équatoriales et écliptiques
exercices
E10-1 Matrices rotations appliquées au changement
de systèmes de coordonnées
Le but de cet exercice est d’obtenir les formules de passage d’un système de
coordonnées à l’autre à partir de la rotation qui les lie.
1. On considère
I I M dans un premier temps, le repère orthonormé {O ; x, y, z} et
sa base { i, j,k } représentés à la figure suivante :
k r
O Φ j y
i
Ψ A’
x
239
Guide de localisation des astres
solutions
1. Matrices rotations
a) Coordonnées de A
Coordonnées ⏐Ψ = (Ox, OA') Coordonnées ⏐x = r cos Φ cos Ψ
polaires : ⏐Φ = (OA', OA) cartésiennes : ⏐y = r cos Φ sin Ψ
⏐ r = OA ⏐z = r sin Φ
b) Matrice R(θx)
Z z
θx
Y
K
J
θx
I O y
X,x
1 0 0 1 0 0
R (θ x ) = 0 cos θ x sin θ x R i (θ x ) = R (− θ x ) = 0 cos θ x − sin θ x
0 − sin θ x cos θ x 0 sin θ x cos θ x
240
10. Coordonnées équatoriales et écliptiques
c) Matrice R(θy)
Z z
θy
O Y,y
x θy I
I I I
Ii = cos
I
θ y I + sin θ y K
Ij = J I I
k = − sin θ y I + cos θ y K
On déduit :
cos θ y 0 − sin θ y
R (θ y ) = 0 1 0
sin θ y 0 cos θ y
d) Matrice R(θz)
Z,z
K
Y
J
θz
O y
x I
θz
exercices
I I I
Ii = cos θzII − sin θz IJ
Ij = sin
I z
θ I + cos θz J
k = K
241
Guide de localisation des astres
On déduit :
cos θz sin θz 0
R (θz ) = − sin θz cos θz 0
0 0 1
2. Formules de passage entre coordonnées horaires et horizontales
a) Repère et coordonnées horaires
• On se réfère à la figure 9.6, dans laquelle p est l’intersection du méridien
local avec l’équateur ; il vient :
Ox = Cp Oz = CPN Oy = Cy (trièdre direct)
• Coordonnées ⏐Ψ = 360° − AHL Coordonnées ⏐ cos D cos AHL
horaires ⏐Φ = D horaires ⏐ − cos D sin AHL
polaires ⏐ r = 1 cartésiennes ⏐ sin D
b) Repère et coordonnées horizontales
LLI
• On se réfère
LLI à la figure 9.2 ; PS y représente la direction Sud dans le plan
horizontal, PZ est la verticale locale, ainsi :
Ox = PS Oz = PZ Oy = Py (trièdre direct)
• Coordonnées ⏐Ψ = 360° − a Coordonnées ⏐ cos h cos a
horizontales ⏐Φ = h horizontales ⏐ − cos h sin a
polaires ⏐ r = 1 cartésiennes ⏐ sin h
c) Angle de rotation et matrice
PN
Z (zénith)
HN
L
équateur
C θy
HS (horizon)
PS
Compte tenu de la distance de l’astre, on confond les centres des deux repères :
P est en C, et dans la figure 10.9, on a représenté les traces de l’équateur et de
l’horizon local dans le plan du méridien du lieu.
On en déduit que le passage du repère horaire au repère horizontal s’effectue
par une rotation d’axe Cy et d’angle :
θy = 90° − L
242
10. Coordonnées équatoriales et écliptiques
Par conséquent :
sin L 0 − cos L
R (θ y ) = R (90° − L ) = 0 1 0
cos L 0 sin L
d) Formules de changement de coordonnées
• Horaires à horizontales :
cos h cos a = cos D cos AHL sin L − sin D cos L
cos h sin a = cos D sin AHL
sin h = cos D cos AHL cos L + sin D sin L
• Horizontales à horaires :
cos D cos AHL = cos h cos a sin L + sin h cos L
cos D sin AHL = cos h sin a
sin D = − cos h cos a cos L + sin h sin L
243
Guide de localisation des astres
solutions
1. Formules de passage entre coordonnées équatoriales et écliptiques
a) Axes de références et coordonnées polaires et cartésiennes
On se reporte à la figure 10.5.
• Coordonnées équatoriales
Axes : Ox = Cγ Oz = CPN Oy = Cy (trièdre direct)
Coordonnées ⏐α ⏐cos δ cos α
polaires : ⏐δ cartésiennes : ⏐cos δ sin α
⏐1 ⏐sin δ
• Coordonnées écliptiques
Axes : Ox = Cγ Oz = CQN Oy = Cy (trièdre direct)
Coordonnées ⏐λ ⏐cos β cos λ
polaires : ⏐β cartésiennes : ⏐cos β sin λ
⏐1 ⏐sin β
b) Matrice R(ε)
La rotation d’axe Ox = Cγ et d’angle θx = ε permet de passer des coordonnées
équatoriales à écliptiques ; la matrice rotation est :
1 0 0
R(ε ) = 0 cos ε sin ε
0 − sin ε cos ε
c) Formules de changement de coordonnées
• Coordonnées équatoriales à écliptiques
cos β cos λ = cos δ cos α (1)
cos β sin λ = cos δ sin α cos ε + sin δ sin ε (2)
sin β = − cos δ sin α sin ε + sin δ cos ε (3)
exercices
244
10. Coordonnées équatoriales et écliptiques
2. Recul du point γ
a) Expressions de dα/dt et dδ/dt en fonction de dλ/dt
• Dérivée des formules (4), (5) et (6) :
dδ dα dλ
sin δ cos α + cos δ sin α = cos β sin λ (7)
dt dt dt
dδ dα dλ
− sin δ sin α + cos δ cos α = cos β cos ε cos λ (8)
dt dt dt
dδ dλ
cos δ = cos β sin ε cos λ (9)
dt dt
Combinant (7) et (8) il vient :
dα dλ
cos δ = (sin α cos β sin λ + cos α cos ε cos β cos λ)
dt dt
dδ dλ
avec : cos δ = sin ε cos β cos λ
dt dt
• Tenant compte des relations (1) et (2), on obtient finalement :
dα dλ
= (cos ε + sin ε sin α tg δ) (10)
dt dt
dδ dλ
= sin ε cos α (11)
dt dt
b) Application numérique
• ε = 23°,439 278 ; pour dt correspondant à un an, on calcule :
dα = (46 ",15 + 20 ", 01sin α tg δ) / an
dα = (3, 077 + 1, 334 sin α tg δ) s / an dδ = 20 ", 01 cos α / an
• Au 1/07/04,
αSirius = 6 h 45 min 20,8 s = 101°,336 67 δSirius = −16° 43’ 21’’ = −16°,7225
– on calcule
dαSirius = 2,7 s dδSirius = − 4’’
αVéga = 18 h 37 min 5,5 s = 279°,272 91 δVéga = 38° 47’ 17’’ = 38°,788 056
– on calcule
dαVéga = 2,0 s dδVéga = + 3’’
Au 1/07/05,
exercices
245
Guide de localisation des astres
1
– 1er cas, α = 90° → tg δ = − = tg(ε − 90°)
tg ε
1
– 2e cas, α = 270° → tg δ = = tg (90° − ε )
tg ε
Les couples invariables sont donc :
(α = 6 h, δ = ε − 90°) qui correspond au pôle Sud écliptique
(β = − 90°, λ indéterminé).
(α = 18 h, δ = 90° − ε) qui correspond au pôle Nord écliptique
(β = + 90°, λ indéterminé) .
d) Quelles sont les heures légales tl de lever et tc de coucher, ce 1er juillet 2004
à Paris ?
3. On veut maintenant connaître les azimuts Zl de l’étoile à son lever et Zc à
son coucher ; quelle relation les lie ?
Montrer que cos Z s’exprime simplement en fonction de la latitude d’observa-
tion et de la déclinaison de l’étoile.
246
10. Coordonnées équatoriales et écliptiques
solutions
1. Angle horaire local AHLe (H) de l’étoile
Pour h = 0 les équations de passage des coordonnées horizontales aux coor-
données horaires s’écrivent :
sin D = cos L cos Z
cos D sin AHL e = − sin Z
cos D cos AHL e = − sin L cos Z
Considérant la troisième et la première équation, on déduit :
sin L sin D
cos AHL e = −
cos D cos L
soit :
cos H = − tgL tgD
La latitude d’observation est L = + 48° 50’ 11’’ = 48°,8364
La déclinaison de Sirius est D (ou δ) = −16° 43’ 21’’ = −16°,7225
• Sans correction de réfraction on calcule
cos H = 0,343 633 H = ± 69°,9016
→ H = ± 4 h 54 min 06 s
• Tenant compte de la réfraction r0 à l’horizon, le cosinus de l’angle horaire
s’écrit :
sin r0 sin L sin D
cos H = − −
cos L cos D cos L cos D
sin r0
r0 = 34’ = 0°,5667 → = 0,015 689
exercices
cos D cos L
Cette expression se déduit de la relation :
sin h = sin L sin D + cos L cos D cos H dans laquelle h = − r0
On obtient donc : cosH = 0,327 944 → H = ± 70°,8560
Au lever : Hl = − 70°,8560 = 289°,1440 → Hl = 19 h 16 min 35 s
Au coucher : Hc = + 70°,8560 → Hc = 4 h 43 min 25 s .
247
Guide de localisation des astres
sin L sin D
On en tire : cos Z = sin D cos L + sin L
cos L
sin D
d’où après simplification : cos Z =
cos L
248
10. Coordonnées équatoriales et écliptiques
249
Guide de localisation des astres
solutions
1. Conditions de visibilité
N Étoile
H
δ L
exercices
H’
S
250
10. Coordonnées équatoriales et écliptiques
Castor 35 32 28 5 40 11 73 02 30 14 45 54 324 27 32 23 51 37
251
Guide de localisation des astres
10
5
0
-5
-10
-15
-20
0 30 60 90 120 150 180 210 240 270 300 330 360
date (jours)
252
Chapitre 11
Repérage des orbites
et éphémérides
Dans un premier temps, nous position- Une première application illustrant ce
nons la trajectoire elliptique de l’astre dernier chapitre porte sur les satellites
par rapport à un plan et une direction semi-synchrones largement utilisés par
de référence : à cette fin, trois angles les soviétiques pour leurs télécommuni-
sont définis. Ayant repéré à un instant cations. L’exercice nous donne le cahier
donné l’astre sur l’ellipse, nous indi- des charges permettant de préciser
quons ensuite comment, grâce à la la trajectoire, puis de déterminer en
connaissance des six variables ellipti- fonction du temps, les coordonnées
ques, nous passons par exemple des géographiques de la trace au sol du
coordonnées polaires de l’astre ou du satellite, que l’on retranscrira sur une
satellite artificiel à ses coordonnées carte du monde.
équatoriales terrestres. L’objet des deux exercices suivants sera
d’obtenir les coordonnées équatoriales
terrestres de certains points particu-
liers de la trajectoire de Vénus. Dans le
premier exercice, nous déterminerons
les dates de passage en ces points et les
coordonnées polaires de ceux-ci. Pour
ces dates, les grandeurs intermédiaires
à calculer seront alors les coordonnées
écliptiques héliocentriques de Vénus et
de la Terre. Dans le dernier exercice,
toujours pour ces mêmes dates, nous
déduirons les coordonnées équato-
riales terrestres de Vénus à partir des
éphémérides afin d’estimer la précision
des résultats précédents.
Guide de localisation des astres
1 Éléments d’orbites
Plan de référence
On le choisit selon le couple d’astres considéré ; par exemple :
– Pour une planète orbitant autour du Soleil, le plan de référence est le plan
écliptique.
– Pour un satellite de la Terre, le plan de référence estI le plan équatorial.
Le plan de référence contient le foyer F, et on notera K le vecteur unitaire qui
lui est normal en F et dirigé vers le Nord (le haut !).
254
11. Repérage des orbites et éphémérides
Ligne de référence
Afin de compléter la référence spatiale, on considère une ligne de référence
située dans le plan de référence :
– si ce plan est le plan équatorial ou écliptique, alors la ligne de référence est
l’intersection de ces deux plans correspondant à la ligne des équinoxes ;
– si le plan de référence est autre, on choisira dans ce plan un axe de référence
passant également
I par F.
On notera I le vecteur unitaire de l’axe de référence s’appuyant sur F ; sur
la ligne des
I Iéquinoxes
I il est dirigé vers le point γ. On définit ainsi une base
directe { I , J , K } du repère de référence.
Ascension droite
L’angle Ω dénommé ascension droite du nœud ascendant NA est mesuré
dans le plan de référence ; c’est l’angle orienté entre la direction de référence
et la ligne des nœuds. Sur l’écliptique, Ω est la longitude du nœud ascen-
dant. I I
Ω est l’angle, compté dans le sens direct, des deux vecteurs I et L asc :
I I
Ω = ( I , L asc )
Inclinaison
L’inclinaison i est l’angle entre le plan de référence et le plan orbital.
I
Soit N est la normale à la ligne des nœuds au nœud ascendant, contenue dans
le plan de référence et orientée vers l’Est.
I
Soit N dep la normale à la ligne des nœuds au nœud ascendant contenue dans
le demi-plan orbital supérieur ; elle est orientée dans le sens de déplacement
du satellite.
L’angle d’inclinaison est alors défini par :
I I
i = ( N est , N dep )
255
Guide de localisation des astres
Argument du péricentre
L’angle ω, argument du péricentre, est mesuré dans le plan orbital. C’est
l’angle orienté entre la ligne des nœuds et la ligne joignant l’apocentre et le
péricentre, appelée ligne des apsides.
I
À partir du centre F, on considère un vecteur L per dirigé vers le péricentre ;
I I
ainsi ω est l’angle, compté dans le sens direct, des vecteurs L asc et L per :
I I
ω = ( L asc , L per )
Si l’inclinaison est faible, la ligne des nœuds est mal définie entraînant
une dégradation de la précision sur l’ascension droite et l’argument du
péricentre. On a ainsi été amené à définir une nouvelle grandeur angu-
laire, la longitude du péricentre ϖ telle que :
ϖ=Ω+ω
D’autre part, l’angle entre le nœud ascendant et la position de l’astre est tel que :
I I I I I I
( L asc , L sat ) = ( L asc , L per ) + ( L per , L sat ) = ω + θ
I I
ω + θ = ( L asc , L sat )
256
11. Repérage des orbites et éphémérides
Date au périhélie
Disposant de la formulation en fonction du temps de la longitude moyenne et
de la longitude du péricentre ; on en déduit, lorsque ces grandeurs diffèrent
d’un multiple m de 360°, la date de passage tp,m des planètes au périhélie.
L’égalité à satisfaire s’écrit :
λ(tp,m − t0) − ϖ (tp,m − t0) = 360 m (m = (0), 1, 2…)
La date t0 à partir de laquelle sont calculés les passages au périhélie correspond
à la référence J2000 (voir paragraphe suivant). Le tableau suivant donne, pour
les planètes Mercure à Neptune, les instants des deux premiers passages au
périhélie exprimés en nombre de jours après t0 et convertis (*) en date calen-
daire (J, M, A, h, min).
257
Guide de localisation des astres
258
11. Repérage des orbites et éphémérides
foyer F ω
i
NA
apocentre Ω
ligne des noeuds
Figure 11.1 Localisation d’un astre par rapport à un plan et une direction de référence.
259
Guide de localisation des astres
260
11. Repérage des orbites et éphémérides
Historique de sa résolution
Dès l’antiquité, l’Almageste de Ptolémée, gigantesque ouvrage en treize
volumes, permettait déjà de calculer pour une date quelconque, les positions
du Soleil, de la Lune et des Planètes.
Du XIIIe au XVIe siècle, trois tables célèbres vont se succéder ; ce sont :
– les Tables Alphonsines, du nom du roi Alphonse de Castille (1221-1284).
– les Tables Pruténiques, dénommées ainsi en l’honneur du roi de Prusse ;
établies par l’astronome allemand Érasme Reinhold (1511-1553) à partir du
modèle copernicien, elles ont détrôné les tables précédentes par la qualité des
éphémérides déduites ;
– les Tables Rodolphines, en hommage à l’empereur Rodolphe II de Hongrie,
éditées à Ulm en 1627.
Une meilleure précision est encore atteinte par ces dernières tables, Kepler
appliquant à la cosmologie copernicienne, délivrée des épicycles, les mesures
rigoureuses des astres effectuées par le danois Tycho Brahé (1546-1601).
261
Guide de localisation des astres
Des éphémérides vont ensuite être régulièrement publiées, par exemple dans :
– la Connaissance des temps en France, dès 1679, grâce à Joachim d’Alencé
(1640-1707) ;
– le Nautical Almanach and Astronomical Ephemeris en Grande-Bretagne,
dès 1766, grâce à l’astronome Nevil Maskelyne (1732-1811).
Les théories planétaires, à partir desquelles ont été construites les éphémérides
récentes, doivent beaucoup aux travaux des Français Laplace (1749-1827) et Le
Verrier (1811-1877) et de l’Américain Newcomb (1835-1909) ; en particulier :
– Laplace apporta une contribution théorique déterminante avec, par exemple,
le théorème sur l’invariabilité des grands axes et des moyens mouvements.
– Le Verrier découvrit Neptune en cherchant les perturbations qu’apporterait à
Uranus, découverte par Herschel en 1781, une planète située sur l’écliptique à
une distance donnée par la loi empirique de Titius-Bode (*).
Sa théorie planétaire a longtemps été la référence.
– La théorie des planètes de Newcomb, utilisée dans les éphémérides améri-
caines jusque très récemment, ne diffère fondamentalement pas de celle de
Le Verrier ; mais les résultats y sont beaucoup plus précis pour les planètes
telluriques.
Les précisions globales données par les théories planétaires actuelles sont
toujours meilleures que la seconde d’angle, la théorie donnant 0″,10 pour les
planètes telluriques.
(*) Titius montra en 1766 que la distance a des planètes au Soleil, exprimée commodément
en ua, suivait une progression mathématique simple. Sa loi fut publiée en 1772 par son compa-
triote allemand Johann Elert Bode (1747-1826). On l’exprime par la relation :
a = 0,4 + 0,3x2n −1 , n ≥ 1
(n est un entier correspondant au rang de la planète. Pour Mercure, n = 0 et a = 0,4 ; pour Jupiter,
n = 5 d’où a = 5,2).
262
11. Repérage des orbites et éphémérides
I I
ω + θ = ( L asc , L sat )
Choisissons le trièdre orthogonal direct {O, x0, y0, z0} du référentiel orbital
tel que :
Ox0 = FNA
Oy0 dans le demi-plan orbital supérieur au plan de référence.
On en déduit les coordonnées orbitales cartésiennes :
x 0 = ρ cos(ω + θ)
y 0 = ρ sin(ω + θ)
z0 = 0
• Une première rotation autour de Ox0 va faire passer au trièdre intermédiaire
{O, xi, yi, zi}
tel que :
Oxi = FNA
θx = (Oy0, Oyi) = − i (i, angle d’inclinaison entre plans équatorial et
orbital).
La matrice rotation correspondante est :
1 0 0
R (− i ) = 0 cos i − sin i
0 sin i cos i
d’où les coordonnées cartésiennes intermédiaires :
x i = ρ cos(ω + θ)
y i = ρ sin(ω + θ) cos i
z i = ρ sin(ω + θ)sin i
• Une deuxième rotation d’axe Ozi , perpendiculaire au plan de référence, fait
passer au trièdre {O, xΩ, yΩ, zΩ} caractérisé par :
OzΩ = Ozi LLLLI LLI
θz = (Oxi, OxΩ) = ( FN A , Fγ ) = − Ω
La matrice rotation est :
cos Ω − sin Ω 0
R(− Ω) = sin Ω cos Ω 0
0 0 1
On en tire les coordonnées cartésiennes :
x Ω / ρ = cos(ω + θ) cos Ω − sin(ω + θ) cos i sin Ω
y Ω / ρ = cos(ω + θ)sin Ω + sin(ω + θ) cos i cos Ω
z Ω / ρ = sin(ω + θ)sin i
263
Guide de localisation des astres
Compte tenu du choix du trièdre {O, xΩ, yΩ, zΩ}, il s’agit des coordonnées
cartésiennes équatoriales.
Les coordonnées équatoriales polaires α et δ satisfont à :
cos(ω + θ)sin Ω + sin(ω + θ) cos i cos Ω
tg α =
cos(ω + θ) cos Ω − sin(ω + θ) cos i sin Ω
sin δ = sin(ω + θ)sin i
On a donc établi les relations permettant de passer de l’anomalie vraie aux
coordonnées équatoriales du satellite.
264
11. Repérage des orbites et éphémérides
planète
S
z
repère écliptique
héliocentrique
T y
repère écliptique
point γ x
géocentrique
x = xP − xT
Elles sont donc : y = yP − yT
z = zP − zT
Les coordonnées (xp , yp , zp) et (xT , yT , zT) correspondent à (xΩ, yΩ, zΩ).
Des coordonnées cartésiennes écliptiques (x, y, z) ainsi calculées, on en déduit
les nouvelles longitude et latitude qui satisfont à :
yP − yT
tg λ =
xP − xT
zP − zT
sin β =
( x P − x T ) + ( y P − y T )2 + (z P − z T )2
2
265
Guide de localisation des astres
exercices
E11-1 Satellites semi-synchrones MOLNYA
Les satellites géostationnaires ne permettent pas de couvrir les régions de
latitude élevée ; en outre, la qualité de la transmission est dégradée s’ils
apparaissent trop bas sur l’horizon. Aussi, dès 1967, les soviétiques ont mis
au point un système de télécommunications par satellites apte à assurer
une couverture permanente du pays. Quatre satellites sont répartis sur une
orbite nettement inclinée sur l’équateur et fortement elliptique. Ces satel-
lites, dénommés MOLNYA (MOLNIYA dans la littérature anglo-saxonne),
« éclair» en russe, ont tous un phasage 2.1 : ils sont semi-synchrones, c’est-
à-dire de période TMol moitié de celle de la Terre.
L’inclinaison i du plan de l’orbite satisfait en fait à la relation 5 cos²i – 1 = 0 .
Ceci permet de stopper la dérive de l’axe des apsides dont sont affectées
les orbites des satellites artificiels du fait de l’aplatissement de la Terre.
D’autre part, au nœud ascendant NA , le satellite est à 65° Est du méridien de
Greenwich.
1. On considère un point M de la surface terrestre de latitude L et longitude
G. Par ailleurs on note T la trace du satellite, c’est-à-dire l’intersection du
rayon FS avec le sol (F est le centre de la Terre et S la position du satellite à
l’instant considéré).
a) Calculer la vitesse angulaire de rotation Ω T de la Terre autour de son
axe des pôles et en déduire, en fonction de L , la vitesse tangentielle VM du
point M ; la calculer pour L = i.
b) On veut que la trace T apparaisse confondue le plus longtemps possible
avec le point M : I
Quelles sont les implications sur le vecteur vitesse VA à l’apogée ?
En déduire la position du plan d’orbite dans le référentiel géocentrique, et
faire un schéma précisant l’orbite autour de la Terre.
2. On se propose maintenant de calculer les caractéristiques de l’orbite d’un
satellite satisfaisant aux contraintes ci-dessus.
a) Connaissant le rayon Rsat d’une orbite géostationnaire, en déduire le grand
axe a de l’orbite MOLNYA.
exercices
266
11. Repérage des orbites et éphémérides
solutions
1. Étude de l’orbite
I
a) VM
2π
La vitesse de rotation de la Terre est : Ω T = = 7,2921 10-5 rd/s
T
exercices
267
Guide de localisation des astres
I
b) Conditions sur VA et plan d’orbite
C’est à l’apogée que la vitesse du satellite est la plus faible ; donc, si possible,
I
la ligne des apsides PA devra passer par M. D’autre part, VA devra être coli-
I I
néaire à VM afin que la vitesse VT de la trace le soit également. En consé-
I
quence, VA sera orthogonal au plan contenant l’axe polaire et la ligne des
apsides. Il n’y a pas d’impossibilité, car l’on sait que pour une trajectoire
elliptique, l’angle α du vecteur vitesse avec le rayon vecteur est précisément
égal à 90° au périgée P et à l’apogée A (voir chapitre 3, paragraphe 4.1).
LLLI
Rayon vecteur FA et vecteur vitesse en A déterminent le plan d’orbite dont
l’intersection avec l’équateur donne la ligne des nœuds. La ligne des nœuds
située dans le plan équatorial est perpendiculaire à l’axe des pôles tout comme
I
VA ; elle est donc parallèle à ce vecteur puisqu’ils sont dans un même plan à
savoir le plan d’orbite.
I
Le plan d’orbite contenant VA , il en découle qu’il est perpendiculaire au plan
défini par l’axe des pôles et la ligne des apsides, tout comme le plan équa-
torial. L’angle i de ces deux plans, qui correspond à l’angle de la ligne des
apsides avec sa projection sur l’équateur, est donc égal à L :
i = 63°,435
La figure 11.3 illustre ces considérations géométriques.
2. Caractéristiques de l’orbite MOLNYA
a) Grand axe a
a3 μ
Pour un satellite de la Terre la 3e loi de Kepler permet d’écrire : = T2
T2 4π
Pour un satellite géostationnaire : a = Rsat et T correspond à la période sidérale ;
pour un satellite MOLNYA : a est l’inconnue, sa période TMol est la moitié
de la période sidérale.
On déduit :
a3 4a 3 3
R sat μ R
a = 3 sat
= = = T2 d’où :
2
TMol T 2
T 2
4π 4
L’application numérique donne : a = 26 562 km
exercices
b) Expression de l’excentricité e
La relation donnant la vitesse VA en fonction du grand axe a et de la période
TMol (voir chapitre 8, paragraphe 1) est :
2 πa 1 − e 4 πa 1 − e
VA = =
TMol 1 + e T 1+ e
268
11. Repérage des orbites et éphémérides
VA
A
ΩT
i
ND
T
NA
P
GNa
Ω
γ
Figure 11.3 Position de l’orbite d’un satellite MOLNYA dans le référentiel géocen-
exercices
trique.
269
Guide de localisation des astres
VA
Le rayon à l’apogée est ρ A = a (1 + e ) et le rapport représente la
ρA
vitesse angulaire θ A par rapport au foyer F. La vitesse angulaire du point T est
bien sûr identique, d’où :
VA 4 πR 1 − e
VT = R =
ρA T (1 + e )3
2π
On veut que VT = VM = R cos i
T
On en tire par conséquent la relation permettant d’accéder à l’excentricité :
1− e cos 2 i
= = 0, 05
(1 + e )3 4
et on calcule que e = 0,737 66
VA
On aurait pu directement écrire que le rapport , qui correspond
ρ A cos i
à la vitesse angulaire de A par rapport à l’axe des pôles, devait être égal
à la vitesse de rotation 2π de la Terre.
T
270
11. Repérage des orbites et éphémérides
1 − e 2 sin θ TMol
a) sin ϕ = b) t = (ϕ − e sin ϕ )
1 + e cos θ 2π
a (1 − e 2 ) 2 πa 1 + e + 2e cos θ
2
c) H = −R d) V =
1 + e cos θ TMol 1 − e2
On en déduit le tableau demandé :
θ ϕ t t H V
degrés degrés secondes h min s km km/s
(P) 0 0,000 0,0 0 h 00 min 00 s,0 590 9,9699
15 5,857 184,8 0 h 03 min 04 s,8 692 9,8865
30 11,888 380,7 0 h 06 min 20 s,7 1 010 9,6380
45 18,286 601,4 0 h 10 min 01 s,4 1 579 9,2289
60 25,288 865,7 0 h 14 min 25 s,7 2 468 8,6669
75 33,204 1 203,8 0 h 20 min 03 s,8 3 789 7,9625
(NA) 90 42,468 1 667,2 0 h 27 min 47 s,2 5 730 7,1296
105 53,713 2 350,9 0 h 39 min 10 s,9 8 587 6,1856
120 67,881 3 437,8 0 h 57 min 17 s,8 12 806 5,1525
135 86,338 5 284,7 1 h 28 min 04 s,7 18 932 4,0608
150 110,819 8 534,3 2 h 22 min 14 s,3 27 147 2,9618
165 142,564 13 986,5 3 h 53 min 06 s,5 35 741 1,9800
(A) 180 180,000 21 541,0 5 h 59 min 01 s,0 39 777 1,5052
Au périgée P : θ = 0° α = Ω − 90°
au nœud ascendant NA : θ = 90° α = Ω
À l’apogée A : θ = 180° α = Ω + 90°
au nœud descendant ND : θ = 270° α = Ω + 180°
271
Guide de localisation des astres
θ° α° θT° = Vγ .t t G° L°
(P) 0 -18,034 0,000 0 h 00 min 00 s,0 18,034 W -63,435
30 34,205 1,591 0 h 06 min 20 s,7 32,614 E -50,768
60 57,488 3,617 0 h 14 min 25 s,7 53,871 E -26,565
(NA) 90 71,966 6,966 0 h 27 min 47 s,2 65,000 E 0,000
272
11. Repérage des orbites et éphémérides
Moscou
60
30
-30
-60
équatoriales de Vénus
• Préambule : Dans ce problème, certaines données sont à rechercher dans
des ouvrages : les éléments moyens de Vénus et de la Terre correspondant aux
trois variables métriques (demi-grand axe a, excentricité e, incidence i) et aux
trois variables angulaires (longitude moyenne λ, longitude du péricentre ϖ ,
273
Guide de localisation des astres
274
11. Repérage des orbites et éphémérides
275
Guide de localisation des astres
q a ρ ρ
e
point ° ′ ′′ ua ua km
P 0° 0,723 329 820 0,006 769 901 0,718 432 949 107 476 039
Nd 125° 05’ 38” 0,723 329 820 0,006 769 799 0,726 122 805 108 626 425
A 180° 0,723 329 820 0,006 769 754 0,728 226 585 108 941 146
exercices
Na 305° 05’ 33” 0,723 329 820 0,006 769 651 0,720 492 618 107 784 161
276
11. Repérage des orbites et éphémérides
277
Guide de localisation des astres
q a ρ ρ
e
point ° ′ ′′ ua ua km
P 78° 38’ 34” 1,000 001 018 0,016 706 860 0,996 443 566 149 065 836
Nd 154° 10’ 40” 1,000 001 018 0,016 706 771 1,014 985 850 151 839 722
A 187° 15’ 29” 1,000 001 018 0,016 706 731 1,016 569 367 152 076 613
Na 262° 56’ 39” 1,000 001 018 0,016 706 641 1,001 777 767 149 863 821
Nd 103° 00’ 49” 153° 20’ 03” 153° 45’ 27” − 0,225 020 182 − 0,989 728 343
− 33 662 540 − 148 061 253
A 103° 00’ 55” 187° 30’ 11” 187° 22’ 48” 0,352 239 775 − 0,953 593 424
52 694 320 − 142 655 546
Na 103° 01’ 08” 264° 50’ 49” 263° 53’ 42” 0,996 357 395 0,104 070 344
149 052 945 15 568 702
Point xV – xT yV – yT zV – zT Distance à la
(date) ua / km ua / km ua / km Terre (ua / km)
P 0,519 832 045 0,565 551 994 0,034 807 310 0,768 951 209
(21/03/04) 77 765 767 84 605 374 5 207 100 115 033 463
Nd 0,058 220 814 0,283 023 101 0 0,288 949 371
(07/06/04) 8 709 710 42 339 653 0 43 226 211
exercices
A 0,130 489 721 0,409 496 615 – 0,035 282 500 0,431 230 680
(12/07/04) 19 520 984 61 259 822 – 5 278 187 64 511 191
Na – 0,830 885 511 0,597 163 332 0 1,023 217 855
(28/09/04) – 124 298 703 89 334 363 0 153 071 212
278
11. Repérage des orbites et éphémérides
• Coordonnées équatoriales α et δ
cos ε sin λ − sin ε tg β
L’ascension droite α se déduit de : tg α =
cos λ
La déclinaison δ se déduit de : sin δ = cos ε sin β + sin ε cos β sin λ
L’obliquité ε dépend légèrement du temps :
ε ≈ 23° 26′ 21″,412 - 468″,096 τ
À la précision affichée, elle vaut 23° 26′ 19″
Le tableau suivant donne la longitude, la latitude, l’ascension droite ainsi
que la déclinaison en « °,′ ,″ » ; l’ascension droite y est également convertie
en « h, min, s ».
279
Guide de localisation des astres
solutions
1. Étapes des calculs et tableau des résultats
a) Étapes
On donne les détails pour la première date :
– conversion en temps des éphémérides → 21/03/2004 à 22 h 31 min 11 s
(ajout de 65 s)
– le tableau qui convient va du 0 mars 0 h au 2 avril 0 h et DT = 33 jours
– jour de l’année du 0 mars (29 février) = 60 ; du 21 mars = 81 → ΔJ = 21
– t – t0 = 21 + 22/24 + 31/1440 + 11/86400 = 21,938 322 jours
– le jour julien correspondant au 29 février à 0 h est 2 453 064,5
– t – t0 = 2 453 086,437 567 + 65/86400 – 2 453 064,5 = 21,938 319 jours (à
retenir)
– x = – 1 + 2 (t – t0)/DT = 0,329 595 110 → θ = Arccos x = 1,234 921 637
b) Tableau
exercices
280
11. Repérage des orbites et éphémérides
– Résultats :
Remarque
Les exercices de ce dernier chapitre nous ont montré comment obtenir en
fonction du temps la trace au sol d’un satellite, ou les coordonnées équa-
toriales terrestres d’une planète à partir de leur position connue sur l’or-
bite. Un des objectifs de cet ouvrage était en effet de nous rendre aptes à
traiter ce genre de problème, mais il nous permettra d’en appréhender bien
d’autres, comme nous avons pu le constater à travers la résolution d’un
certain nombre d’exercices proposés dans les chapitres précédents.
Toutefois, au-delà de cet aspect didactique, nous avons également voulu,
lorsque l’occasion se présentait, préciser notre environnement spatial et
exercices
évoquer les noms de certains scientifiques qui ont permis de percer les
secrets de la mécanique qui l’anime.
Puisse cet ouvrage inciter le lecteur à voyager dans l’espace en observant
par exemple le mouvement des astres sous la vaste coupole d’un planéta-
rium, ou encore en suivant les périples de véhicules de croisière interpla-
nétaires et autres sondes spatiales… Il pourra alors choisir sa destination
281
Guide de localisation des astres
282
Index
A Astroïde 95
Albédo 20 Austral 223
Altitude 190 Axe
Amas des équinoxes 223
galactique 24 focal 46
globulaire 24 non focal 46
Angle Axe du monde 223
médian 171 Azimut 190
Angle horaire
local 196 B
origine 195 Boréal 223
sidéral 225 Bulbe 23
Année Bureau des longitudes 21
julienne 28
sidérale 27
tropique 27
C
Année-lumière 34 Calendrier
Anomalie grégorien 28
excentrique 88 julien 28
moyenne 88 Catalogue FK4 223
vraie 88 Ceinture principale 22
Apex 24 Centre
Aphélie 47 attractif 110
Apoastre 47 Centre de courbure 75
Apocentre 46 Cercle
Apogée 47 d’égales hauteurs 193
Argument du péricentre 256 équatorial 188
Ascension horaire de l’astre 195
verse 225 osculateur 96
Ascension droite 225 principal 88
du nœud ascendant 255 Circularisation 157, 207
Astéroïde 21 Coefficient d’aplatissement 47
Astre 20 Colatitude 190
283
Guide de localisation des astres
Comète 22 Formules
Conique 44 de Binet 61
Conjonction 29 Foyer 44
Constante Fuseaux horaires 234
de gravitation universelle 110
dynamique
G
première 60
géocentrique de la gravitation 111 Galaxie 23
héliocentrique de la gravitation Géodésique 205
110 Géosynchrone 139
Constante de Gauss 26 Grand axe 47
Constante dynamique Groupe local 24
deuxième 130
Constellation 220 H
Coordonnées
horizontales 190 Hauteur 190
terrestres 188 Heure
Crépuscule 218 légale 234
Horizon du lieu 190
D
Déclinaison 194 I
Déférent 72 Inclinaison 255
Développée 94 Intégrale elliptique 167
Déviation gravitationnelle 151 Intensité 222
Directrice 44
Distance
focale 47
J
polaire 194 Jour
zénithale 192 sidéral 25
solaire 25
E Jour julien 36
Écliptique 21, 27 modifié 36
Effet catapulte 156
Éléments moyens 258 L
Enveloppe 94
Latitude 188
Épicycle 72
à la méridienne 202
Équant 72
Équateur céleste 27, 223 écliptique 229
Équation Ligne
de Kepler 89 des apsides 256
du temps 232 des équinoxes 255
Étoile 20 des nœuds 254
double 125 Loi des aires 60
polaire 28 Longitude 188
Excentricité 44 du péricentre 256
écliptique 229
F moyenne 257
Formule Luminosité 222
de Ramanujan 168 Lunaison 39
284
Index
M Périhélie 47
Magnitude Période
absolue 222 sidérale 29
apparente 222 synodique 28
Méridien 188 Petit axe 47
Phase 39
céleste 223
Pied de l’astre 192
du lieu 191
Planète 20, 110
origine 188
jovienne 21
Mille marin 190
tellurique 21
Mois
Point
anomalistique 41
astronomique 203
draconitique 41
Point gamma 27
sidéral 37 Point vernal 27
synodique 38 Pôle écliptique 229
tropique 41 Potentiel gravitationnel 128
Monture Précession des équinoxes 27
alt-azimutale 193 Problème de Lambert 180
équatoriale 197
Moyen mouvement 88 Q
Quadrature 29
N
Nadir 190 R
Nœud Rayon de courbure 75
ascendant 254 Référentiel
descendant 254 barycentrique 115
Nomenclature des comètes 101 de Copernic 115
Nuage de Oort 22 Galiléen 113
Nuit civile 251 héliocentrique 110
héliocentrique ou de Kepler 115
O relatif 115
Obliquité 27, 229 Réfraction 216
Rétrograde 273
Opposition 29
Révolution 27
Orbite
Rotation 25
de Clarke 138
géostationnaire 138
S
géosynchrone 138
Orbite d’insertion 157 Satellite 22
constellation de 23
semi-synchrone 266
P synchrone 139
Parallaxe annuelle 35 Sens rétrograde 25
Parallèle 188 Sextant 192
Paramètre 45 Sonde 146
d’impact 151 Sphère céleste 27, 223
Parsec 34 Sphère d’influence 152
Périastre 47 Spirale
Péricentre 46 de Descartes 82
Périgée 47 hyperbolique 84
285
Guide de localisation des astres
T V
Temps Variable
Atomique International 235 angulaire 258
civil 233 elliptique 258
des éphémérides 235 métrique 258
moyen 232 Variation horaire 197
sidéral 226 Vertical de l’astre 191
local 226 Vitesse
sidéral local 236 aréolaire 60
cosmique
Sidéral Moyen de Greenwich 238
deuxième 134
Universel 233
première 134
universel coordonné 235
d’entraînement 71, 146, 207
vrai 231 d’injection 146
Termes séculaires 260 de lancement 134
Titius-Bode (loi de) 262 de libération 134
Trace 266 minimale de satellisation 133
Trajectoire à énergie minimale 150 Voie lactée 23
Triangle de position 199
Z
U Zénith 190
Unité astronomique 32 Zodiaque 222
286