MISSION GALILEO Trajectoire C3 : De la Terre à la Terre en passant par l'astéroïde GASPRA |
Téléchargement possible de tout un groupe d'exercices
III VOL TERRE - TERRE( Ellipse C3)
(
Voir documentation )Les USA ont en octobre 1989 lancé une sonde spatiale nommée Galiléo, destinée après un voyage très tourmenté et complexe, à atteindre la planète Jupiter pour l'étudier ainsi que ses satellites.
Galiléo emprunte une orbite héliocentrique C1 qui, partant de la Terre, survole la planète Vénus en y réalisant un premier tremplin gravifique le 9/02/1990 la sonde et renvoyée vers la Terre où elle effectue un second tremplin gravitationnel le 8 décembre 1990.
Nous nous intéressons à l'orbite C3 après le tremplin sur la Terre spécialement étudié pour amener Galileo à survoler à nouveau la Terre 2 ans après, afin d'y réaliser l'ultime tremplin et gain de vitesse, qui projettera la sonde vers sa planète de destination, Jupiter
NB: Les dates non calculées seront toujours prises à 0 H, car nous n'avons pas la prétention de calculer de telles missions à l'heure près.
NB 2 : Ce problème vous donnera l'occasion d'utiliser les routines du BDL( qui ne fonctionnent qu'en intranet ou si vous les téléchargez pour les exécuter chez vous)
Téléchargement éventuel des routines
EPHEMBDL.ZIPDonnées utiles de ce vol
:Constante de gravitation du Soleil: mS = 13.271244 1010 km3s-2
Constante de gravitation de la Terre: mTerre = 39.86 104 km3s-2
Rayon de la Terre : RT = 6378 km
Arrivée et survol de la Terre : 8 décembre 1990
Survol de Gaspra le 29 novembre 1991
Hypothèses de calcul : Nous adoptons obligatoirement l'hypothèse des sphères d'influence. On rappelle que vue de la planète la sphère d'influence a un rayon infini et que vue du Soleil, elle a un rayon négligeable donc nul.
QUESTIONS
( NB : Solution dans cette même page )
1°)
En intranet exécutez la routine planeph.exe en allant la chercher dans le répertoire MECASPA/EPHEMERISur Internet, vous devez télécharger la routine et ses ressources(voir plus haut) et l'exécuter ensuite.
Ceci, pour récupérer les coordonnées rectangulaires et la vitesse de la Terre le 8/12/1990 .
NB : Vous trouvez toutes explications sur l'utilisation de la routine planeph.exe dans
ephemeri.htm.2°) A PROPOS DE GASPRA( Astéroïde N°951) :
a) Quel est le demi grand axe de l'orbite de Gaspra?
NB : Les dates juliennes s'obtiennent avec
DATEJULI.EXE et les dates calendaires avec DATE_CAL.EXE.b) Connaissant les paramètres orbitaux de GASPRA, utilisez la routine
PAR_RV_S.EXE d'abord pour déterminer un passage au périgée, puis pour trouver ses coordonnées héliocentriques, sa vitesse et sa distance au Soleil le 29/10/1991. Vous pourrez vérifier que GASPRA se trouve au périgée de son orbite le 7 avril 1989 à 7 h 11 mn.3°) PROBLEME DE LAMBERT :
On souhaite déterminer C3 en utilisant la résolution du
problème dit de Lambert:"Il n'existe qu'une seule conique képlérienne, joignant 2 points donnés du système solaire avec une durée de voyage imposée"
La solution numérique de ce problème, dans le cas elliptique est fournie par 2 routines au choix sur le site
DEUX_PTS.EXE ou LAMBERT1.EXE que vous allez mettre en œuvre. ( Eventuellement vous pouvez comparer les 2 routines )a ) Donner pour Terre-Gaspra, la durée du voyage en jours entre les 2 dates( Utilisez 2 fois la routine
datejuli.exe )b) Avec les positions vectorielles de la Terre et de Gaspra, exécutez
LAMBERT1.EXE pour déterminer les caractéristiques orbitales de C3. Les valeurs des anomalies excentriques au départ et à l'arrivée vous renseigneront sur le parcours emprunté sur C3.NB : CALCULEZ LA PERIODE SUR C3 ET COMPAREZ A CELLE DE LA TERRE ( 365.25 jours)?
c) Pour calculer la vitesse de Galileo au départ près de la Terre, ou au passage près de Gaspra, vous utilisez la routine
PAR_RV_S.EXE.Quelle est la vitesse à l'infini de Galileo à la sortie de la sphère d'influence de la Terre ?
Au passage près de Gaspra ?
4°) SURVOL HYPERBOLIQUE DE LA TERRE LE 8/12/1990 :
a) NB : La vitesse à l'infini d'entrée dans la sphère d'influence de la Terre a été calculée dans un autre exercice ( Galileo de Vénus à la terre ) et valait 8.909 km/s de composantes héliocentriques ci-contre : |
Vx = -0.62 km/s Vy = - 8.727 km/s Vz = 1.68 km/s |
Comparer les normes des vitesses à l'infini d'entrée et de sortie de la sphère d'influence de la Terre.
La documentation donne une altitude de survol au périgée de 1000 km. Calculer alors
a) Le demi angle d'ouverture des asymptotes de l'hyperbole.
b) L'effet
DV de tremplin gravitationnel exercé par la planète Terre. Commentez. Comparez au gain réel sur la norme de la vitesse héliocentrique.Comparez avec le résultat de
Galileo de Vénus à la terre5°) SURVOL HYPERBOLIQUE DE LA TERRE LE 8/12/1992 :
Vous avez maintenant en main tous les outils de calcul et pourrez donc, le 8/012/1992, 2 ans après le précédent tremplin sur la Terre, calculer :
La vitesse héliocentrique de Galileo, celle de la Terre, la vitesse à l'infini d'arrivée dans la sphère d'influence de la Terre.
SOLUTION
1°) Coordonnées rectangulaires et vitesse de la Terre le 8/12/1990.
Pour la Terre le 8/12/1990 on a :
Position X= 36.269106 km X= 142.835106 km X= 3391 km |
Vitesse Vx= - 2535905 km/j = - 29.351 km/s Vy= 624591 km/j = 7.229 km/s Vz= 93 km/j = 0.00107 km/s |
2°) A PROPOS DE GASPRA( Astéroïde N°951) :
a) Quel est le demi grand axe de l'orbite de Gaspra?
Le moyen mouvement n permet de calculer la période et donc le demi grand axe de Gaspra.
Le BDL donne le moyen mouvement de Gaspra, n = 0.3001555054°/jour ====> Période T=1999.378 jours et un demi grand axe a(Gaspra) = 330.494 106 km
b) Passage au périgée. Coordonnées héliocentriques le 29/10/1991 et distance au Soleil.
Une anomalie moyenne de 173°.099744 correspond à une anomalie excentrique pour la position
j = 174°.12084Le 5/11/90; le temps écoulé depuis le périgée est donné par la routine: 576.70021767 jours solaires moyens, ce qui place le périgée à la date julienne JJ = - 3921.20021767 soit 7 avril 1989 à 7 h 12 mn.
Depuis le passage au périgée jusqu'au 29/10/1991 il s'est écoulé 934.7002 jours. Ou bien on peut aussi calculer l'anomalie moyenne au 29/10/1991 alors qu'il s'est écoulé 358 jours depuis le .5/11/1990
M(29/10/1991) = 173°.099744+0.3001555054*358 = 280°.55541 Une nouvelle exécution de PAR_RV_S.EXE donne pour le 29/10/1990:
Position de Gaspra X= 70.519 106 km X= - 322.067 106 km X= 11.708 106 km |
Vitesse de Gaspra Vx= 18.518 km/s Vy= 7.672 km/s Vz= 1.103 km/s |
La distance au Soleil est donc de 329.904 106 km
.3°) DETERMINATION DE C3 :
a) Durée du voyage entre les 2 dates.
Avec datejuli.exe on trouve les dates juliennes:
JJ(8/12/1990)= - 3311.5 et JJ(29/10/1991)= -2986.5, ce qui donne une durée de voyage de 325 jours solaires moyens.
3°) Caractéristiques et paramètres orbitaux de C3.
a) L'exécution de LAMBERT1.EXE donne les résultats suivants :
Paramètres orbitaux :
a = 237.7498 106 km
e = 0.43053
i = 4°.538
w
= 223°.098W
= 255°.769j
d = 332°.007= 5.7946153 rd ( anomalie au point de départ = Terre )j
a = 154°.201( anomalie au point d'arrivée = Gaspra )Ces résultats correspondent à un trajet comme ci-dessous, avec passage au périgée, arrivée sur Gaspra avant l'apogée dans la trajectoire montante:
REMARQUE CAPITALE
: Les experts du JPL ont " parfaitement calculé leur coup " en effectuant un survol de la Terre diaboliquement précis pour imprimer à C3 une période de 731.8 jours, double de celle de la Terre, de manière à retrouver celle-ci exactement 2 ans après, pour un troisième et ultime tremplin.b) Calculer la vitesse de Galileo soit au départ près de la Terre, soit à l'arrivée près de Vénus
On exécute 2 fois PAR_RV_S.EXE ce qui fournit confirmation des positions et apporte en plus les vitesses héliocentriques de Galileo.
Vitesse au départ de la Terre Vx = - 35.136 km/s Vy = - 0.9751 km/s Vz = - 2.722 km/s V = 35.255 km/s |
Vitesse au passage GASPRA Vx = 15.656 km/s Vy = 0.2037 km/s Vz = 1.2004 km/s V = 15.7 km/s |
4°) SURVOL HYPERBOLIQUE DE LA TERRE LE 8/12/1990 :
Vitesse à l'infini de Galileo à la sortie de la sphère d'influence de la Terre ?
Au passage près de Gaspra ?
On calcule par différence la vitesse relative de Galileo, par rapport à la planète Terre et on trouve :
Vitesse à l'infini de sortie de la sphère d'influence de la Terre le 08/12/1990 Vx = -5.785 km/s Vy = - 6.254 km/s Vz = - 2.723 km/s Norme de Vinfini2 = 8.944 km/s donc C3V2 = 79.994 km²/s2 |
Vitesse héliocentrique au passage près de Gaspra le 29/10/1990 Vx = 15.656 km/s Vy = 0.2037 km/s Vz = 1.2004 km/s |
a) Le demi angle
b d'ouverture des asymptotes de l'hyperbole.Il peut se calculer de 2 manières:
1) Par les vitesses à l'infini :
2) Par le périgée de l'hyperbole :
La vitesse à l'infini d'entrée, permet de calculer, par conservation de l'énergie, le demi grand axe de l'hyperbole :
a =
mT/C3T2 = 4983. km or Rp = a(e-1) =1000+6378 donne e= 2.4690685 = 1/cosb donc b = 66°.226NB : La petite différence tient à notre ignorance des corrections de trajectoire et aux simplifications d'usage et aussi à notre choix des dates à 0 h( une plus grande précision eut été ridicule face à l'hypothèse des sphères d'influence )
Le problème précédent avait fourni
b = 66°.108 ce qui est satisfaisant, compte tenu que nous utilisons l'hypothèse simplificatrice des sphères d'influence et que nous ignorons les petits ajustements inévitables sur la vitesse. D'ailleurs nous aurions dû trouver la même norme de vitesse à l'infini en entrée et en sortie ( notre erreur est de 0.39% par rapport à la moyenne des vitesses)b) Gain réel de vitesse héliocentrique. La vitesse de Galileo avant tremplin était de 30.055 km/s, et atteint 35.255 après tremplin. Le
gain réel de vitesse est la différence des normes donc vaut 5.2 km/s, ce qui est conséquent. Rappelons que l'effet total de tremplin ( différence des vecteurs vitesses ) est de 7.224 km/s, calculé ci-après.b) L'effet
DV de tremplin gravitationnel exercé par la planète Terre.Nous savons d'après le cours que
DV = 2 Vinfini cosb = 7.224 km/sL'effet de tremplin remplace au moins 2 moteurs qu'il aurait fallu emporter pour réaliser la modification du vecteur vitesse héliocentrique entre l'arrivée et la sortie de la sphère d'influence de la Terre. Un tel tremplin est déjà relativement important. Une étude détaillée montrerait que le gain de vitesse en norme est sensible et que le tremplin a également servi à "faire tourner" le vecteur vitesse.
5°) SURVOL HYPERBOLIQUE DE LA TERRE LE 8/12/1992:
Coordonnées rectangulaires et vitesse de la Terre le 8/12/1992.
Position X= 35.017 106 km X= 143.13 106 km X= 2179 km Norme = R = 147.351 106 |
Vitesse Vx= - 29.41 km/s Vy= 6.962 km/s Vz= 0.0001 km/s |
Naturellement ces valeurs sont très proches de celles du 8/12/1990
Vitesse de Galileo le 8/12/1992.
Du 8/12/1990 au 8/12/1992 il s'est écoulé 731 jours, l'anomalie moyenne sur l'orbite de Galileo vaut en radians:
M(8/12/1992)= M(8/12/1990)+360/T*731 = 5.7946153- 0.43053*sin(332.007) + 2
p*731/731.797 =12.273033M(8/12/1992) = 703°.19298 mod 360 = 343°.19298
PAR_RV_S.EXE donne alors une vitesse héliocentrique et une vitesse à l'infini d'arrivée:
Vitesse héliocentrique Vx= - 35.03 km/s Vy= 0.014 km/s Vz= - 2.722 km/s De norme 35.136 km/s |
Vitesse à l'infini par rapport à la Terre Vx= - 5.62 km/s Vy=- 6.948 km/s Vz= -2.722 km/s De norme 9.342 km/s |
NB : La vitesse à l'infini d'entrée dans la sphère d'influence de la Terre est de 9.342 km/s, un peu différente de celle que l'on trouvera en sortie, après le tremplin Terre et qui vaut 8.855 km/s. Normalement ces 2 valeurs sont très voisines. Alors où est le problème?
Réponse probable : le détour par Gaspra est certes intéressant, mais ne donne pas exactement une trajectoire qui rejoint parfaitement la Terre. Certes c'est très près, mais une petite correction de trajectoire a dû être nécessaire;
Si on recherche la trajectoire C3* qui joint Gaspra le 29/10/1991 à la Terre le 8/12/1992, on trouve avec DEUX_PTS.EXE :
Paramètres orbitaux de C3* :
a = 237.84 106 au lieu 237.7498 106 km pour C3
e = 0.42948 au lieu de 0.43053 pour C3
i = 4°.6169 au lieu de 4°.538 pour C3
w
= 222°.54 au lieu de 223°.098 pour C3W
= 256°.26 au lieu de 255°.769 pour C3j
d = 154°.33Le calcul de la vitesse à l'infini est alors beaucoup plus satisfaisante, 8.904 km:s en norme à comparer à 8.855 km/s; soit 0.55 % d'écart. Voir la trajectoire suivante C4 de la Terre vers Jupiter
Guiziou Robert février 2002