Etape 11 : Correction : Mouvement d’un satellite géostationnaire (version professeur niveau avancé)

Télécharger le pdf | Télécharger le notebook | Lancer le notebook sur binder (lent)

Document :

GOES-17 est le deuxième satellite de la génération actuelle de satellites météorologiques exploités par l’Administration nationale des océans et de l’atmosphère (NOAA). Il s’agit d’un satellite géostationnaire qui vise à fournir des images haute résolution visibles et infrarouges et des observations de la foudre sur plus de la moitié du globe. Le satellite a été lancé dans l’espace le 1er mars 2018 par un véhicule Atlas V (541) depuis la base aérienne de Cape Canaveral, en Floride. Il avait une masse de lancement de 5 192 kg (sa masse sèche (sans le carburant (ergols)) est de 2 857 kg). Le 12 mars, GOES-17 a rejoint GOES-16 (lancé en 2016) sur une orbite géosynchrone à 35 786 km au-dessus de la Terre (soit un rayon orbital de 42 164 km). GOES-17 est devenu opérationnel le 12 février 2019 sous le nom de GOES-West. Sa durée de vie utile prévue est de 15 ans.

L’orbite géosynchrone est une orbite géocentrique sur laquelle un satellite dit géostationnaire se déplace dans le même sens que la Terre (d’ouest en est) et dont la période orbitale est égale à la période de rotation sidérale de la Terre (soit environ 23 h 56 min 4 s). Un satellite géostationnaire reste donc toujours à la verticale d’un même lieu sur Terre.

source : Wikipédia (texte remanié)

satellite.jpg

orbites.png

Rappels mathématiques :

cercle.jpg

Les coordonnées cartésiennes du point M décrivant un cercle de rayon R centré sur l’origine O du repère sont :

\((x_M=R\times \cos{\theta}\;;\;y_M=R\times \sin{\theta})\)

Le vecteur unitaire \(\overrightarrow{u}=\frac{\overrightarrow{OM}}{OM}=\frac{\overrightarrow{OM}}{R}\) a pour coordonnées : \(\overrightarrow{u} \begin{pmatrix} \frac{x_M}{R}=\cos{\theta} \\ \frac{y_M}{R}=\sin{\theta} \end{pmatrix}\)

Problématique : Comment la force d’attraction gravitationnelle exercée par la Terre sur le satellite GOES-17 influence la variation de son vecteur vitesse moyenne? Quel paramètre est-il important de prendre en compte lors de cette étude?

L’étude du mouvement du satellite GOES-17 aura lieu dans le référentiel géocentrique supposé galiléen auquel on associe un repère cartésien orthonormé fixe dont l’origine est au centre de la Terre.

  1. Quelle est la nature du mouvement du satellite GOES-17 ? (Répondre dans la cellule ci-dessous en double-cliquant dessus si besoin).

  2. En vous basant sur vos connaissances issues de la classe de seconde (en physique et en programmation), réfléchir aux différentes parties que comportera le programme permettant de répondre à la problématique. (Répondre dans la cellule ci-dessous en double-cliquant dessus si besoin).

  • codage de la trajectoire du satellite (cellule 2)

  • codage des coordonnées du vecteur vitesse (cellule 3)

  • codage des coordonnées du vecteur variation de vitesse (cellule 4)

  • codage des coordonnées du vecteur force d’attraction gravitationnelle (cellule 5)

  • affichage d’un graphique représentant la trajectoire, le vecteur variation de vitesse du satellite et le vecteur force d’attraction gravitationnelle (cellule 6)

[1]:

# cellule 1 : import des bibliothèques

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

  1. Quelle valeur en radian prend l’angle \(\theta\) lorsqu’il est parcouru par le segment OM en une période T? (Répondre dans la cellule ci-dessous en double-cliquant dessus si besoin).

note codage LaTeX : double-cliquer sur cette cellule pour voir comment coder l’écriture des lettres grecques :

  • théta : \(\theta\)

  • pi: \(\pi\)

  1. En déduire l’expression de l’angle \(\theta\) en fonction du temps t et de la période T (ligne 8 de la cellule 2).

  2. En déduire les coordonnées x et y de la position du satellite en vous aidant des rappels mathématiques (lignes 10 et 11 de la cellule 2).

Note codage : la constante pi ainsi que les fonctions cos et sin sont fournies par la bibliothèque numpy : - np.pi - np.cos() - np.sin()

[2]:

# cellule 2 : coordonnées de la position du satellite

# Fonction permettant de calculer les coordonnées
# de la position du satellite

def coordposition (t, T) :

    theta=2*np.pi/T*t

    x=R*np.cos(2*np.pi/T*t)
    y=R*np.sin(2*np.pi/T*t)

    return x,y

Afin de calculer les coordonnées du vecteur vitesse, notées vx et vy on crée une fonction coordvit(u) :

  • u est une liste et représente une des coordonnées (x ou y) du vecteur position.

  • vu est une liste et représente une des coordonnées (vx ou vy) du vecteur vitesse.

  • vui est la valeur à l’instant \(t_i\) de la coordonnée étudiée du vecteur vitesse. La liste vu contiendra ces valeurs vui.

  1. Compléter la ligne 6 de la cellule 3 permettant de calculer les valeurs vui prises par la coordonnée vu du vecteur vitesse à chaque instant du mouvement.

Note codage : la variable t est déclarée dans le programme principal et est donc globale. Elle est ainsi reconnue au sein de toute fonction…

On rappelle que la :math:`i^{ème}` valeur d’une liste u est codée u[i]

  1. Ecrire les deux lignes de code permettant de calculer les valeurs des coordonnées vx et vy (lignes 10 et 11 de la cellule 3).

[3]:

# cellule 3 : coordonnées du vecteur vitesse du satellite

def coordvit(t,u):
    vu=[]
    for i in range (len(u)-1):
        vui=(u[i+1]-u[i])/(t[i+1]-t[i])
        vu.append(vui)
    return vu

On appelle vecteur variation de vitesse au temps \(t_i\), le vecteur : \(\overrightarrow{\Delta v}(t_i)=\overrightarrow{v}(t_{i+1})-\overrightarrow{v}(t_{i})\).

On désire calculer les coordonnées notées dvx et dvy du vecteur \(\overrightarrow{\Delta v}\)

  1. En vous basant sur le modèle de la fonction précédente, créer une fonction coordvarvit(vu) permettant de calculer les valeurs d’une coordonnée notée dvu du vecteur variation de vitesse (lignes 3 à 8 de la cellule 4).

  2. Ecrire les deux lignes de code permettant de calculer les valeurs des coordonnées dvx et dvy (lignes 10 et 11 de la cellule 4).

[4]:

# cellule 4 : coordonnées du vecteur variation de vitesse

def coordvarvit(vu):
    dvu=[]
    for j in range (len(vu)-1):
        dvuj=vu[j+1]-vu[j]
        dvu.append(dvuj)
    return dvu

  1. A l’aide du document, déterminer la valeur de la masse m du satellite après avoir consommé 1000 kg de carburant (ligne 4 de la cellule 5).

  2. Donner l’expression vectorielle de la force d’attraction gravitationnelle exercée par la Terre sur le satellite \(\overrightarrow{F_{T/S}}\) en fonction de G, \(M_T\), m, R et \(\overrightarrow{u}\). (Répondre dans la cellule ci-dessous en double-cliquant dessus si besoin).

\(\overrightarrow{F_{T/S}}=-G\times\frac{M_T\times m}{R^2}\times \overrightarrow{u}\)

note codage LaTeX : double-cliquer sur cette cellule pour voir comment coder l’écriture :

  • d’un vecteur : \(\overrightarrow{vecteur}\) ou \(\vec{vecteur}\)

  • d’une fraction : \(\frac{numérateur}{dénominateur}\)

  • d’un indice : \(x_{indice}\)

  • d’un exposant : \(y^{exposant}\)

  • d’un signe x : \(\times\)

  1. Déterminer les expressions des coordonnées Fx et Fy de la force d’attraction gravitationnelle exercée par la Terre sur le satellite en vous aidant de vos connaissances et des rappels mathématiques (lignes 5 et 6 de la cellule 5).

[5]:

# cellule 5 : coordonnées du vecteur
# force d'attraction gravitationnelle

def coordforce (x,y,MT,m,R):
    Fx=(G*MT*m/(R**2))*(-x/R)
    Fy=(G*MT*m/(R**2))*(-y/R)
    return Fx,Fy

  1. Compléter la ligne de code 7 permttant d’afficher la trajectoire du satellite.

  2. Compléter les lignes de code 10, 11 afin d’afficher les légendes des axes.

  3. Sur le modèle des lignes de code 19, 20 et 21 , compléter les lignes de code 22, 23 et 24 permettant de tracer le vecteur force d’attraction gravitationnelle.

[6]:

# cellule 6 : Tracé du graphique permettant de visualiser
# la trajectoire du satellite ainsi que  les
# vecteurs variation de vitesse et force d'attraction
# gravitationnelle en fonction de la durée Te entre
# deux positions successives du satellite.

def trace (x,y,dvx,dvy,Fx,Fy,Te,n) :
    plt.subplot(1,3,n)
    plt.plot(x,y,'k+')
    plt.plot(0,0,'go')
    plt.text(1000000,1000000,"Terre")
    plt.xlabel('x(m)')
    plt.ylabel('y(m)')
    plt.xlim(-50000000,50000000)
    plt.ylim(-50000000,50000000)
    plt.text(20000000,40000000,"Te="+ str(Te) +" jours \n")

    for k in range(len(dvx)):
        plt.arrow(x[k],y[k],30000*dvx[k],30000*dvy[k],
                  facecolor='b',edgecolor='b',width=200000,
                  head_width=1000000,length_includes_head=True)
        plt.arrow(x[k],y[k],10000*Fx[k],10000*Fy[k],
                  facecolor='r',edgecolor='r',width=200000,
                  head_width=1000000,length_includes_head=True)

  1. A l’aide du document, déterminer les valeurs du rayon R de la trajectoire et la période de révolution T du satellite (lignes 3 et 4 de la cellule 7).

[7]:

# cellule 7 : Programme principal permettant d'afficher
# 3 graphiques avec des durées entre deux positions successives
# du satellite différentes

R=42164000   # rayon en m
T=84164      # période de révolution en s

MT=5.972*10**24          # masse de la Terre en kg
m=4192                   # masse du satellite en kg
G=6.67408*10**(-11)      # constante de gravitation universelle
                         # en m^3.kg^-1.s^-2

plt.figure (figsize=(15,5))


# Graphique 1
Te1=3000      # durée entre deux positions successives en jours
t1=np.arange(0,T,Te1)
x1,y1=coordposition (t1,T)
vx1=coordvit(t1,x1)
vy1=coordvit(t1,y1)
dvx1=coordvarvit(vx1)
dvy1=coordvarvit(vy1)
Fx1,Fy1=coordforce(x1,y1,MT,m,R)
trace(x1,y1,dvx1,dvy1,Fx1,Fy1,Te1,1)

# Graphique 2
Te2=2000      # durée entre deux positions successives en jours
t2=np.arange(0,T,Te2)
x2,y2=coordposition (t2,T)
vx2=coordvit(t2,x2)
vy2=coordvit(t2,y2)
dvx2=coordvarvit(vx2)
dvy2=coordvarvit(vy2)
Fx2,Fy2=coordforce(x2,y2,MT,m,R)
trace(x2,y2,dvx2,dvy2,Fx2,Fy2,Te2,2)

# Graphique 3
Te3=1000      # durée entre deux positions successives en jours
t3=np.arange(0,T,Te3)
x3,y3=coordposition (t3,T)
vx3=coordvit(t3,x3)
vy3=coordvit(t3,y3)
dvx3=coordvarvit(vx3)
dvy3=coordvarvit(vy3)
Fx3,Fy3=coordforce(x3,y3,MT,m,R)
trace(x3,y3,dvx3,dvy3,Fx3,Fy3,Te3,3)

plt.title ("Comparaison du vecteur variation de vitesse "
               "et de la force d'attraction gravitationnelle \n"
               "Cas du mouvement d'un satellite géostationnaire "
               "autour de la Terre \n",horizontalalignment='right')

plt.show()
../_images/etape_11_mouvement_satellite_geostationnaire_prof_niveau_avance_29_0.png

  1. Conclusion : Répondre à la problématique dans la cellule suivante en double-cliquant dessus si besoin.