Mouvement d'un solide ponctuel en microgravité ( mécanique )

référentiel géocentrique : solide formé par le centre de la terre et par les centres de 3 étoiles lointaines ; l'origine du repère est le centre de la Terre ; les trois axes pointent vers des étoiles lointaines qui paraissent fixes.

Le référentiel géocentrique est pratiquement galiléen au voisinage de la Terre si on ne tient pas compte des attractions des autres astres.

Placer un véhicule spatial ( masse m) ur une orbite circulaire autour de la Terre à l'altitude h = 300 km :

L'étude est faîte dans le référentiel géocentrique

1^ère méthode : ce véhicule est lancé verticalement depuis la surface de la Terre avec une vitesse v₀ afin qu'il atteigne l'altitude h ; sa vitesse devient nulle ; on lui communique alors une vitesse horizontale v₁.

Calculs des vitesses en écrivant la conservation de l'énergie mécanique lors de la montée :

énergie potentielle à la distance r du centre de la Terre : -mMG/r ( l'origine de cette énergie potentielle est prise à l'infini)

à la surface de la Terre, l'énergie mécanique vaut : ½mv₀² - mMG/R

à l'altitude h où la vitesse est nulle, l'énergie mécanique vaut : - mMG/(R+h)

L'énergie mécanique se conserve : ½mv₀² - mMG/R = - mMG/(R+h) ; v₀² =2MG[ 1/R - 1/(R+h]= 2MGh / (R(R+h))

De plus en assimilant le poids à l'attraction terrestre ( à la surface de la Terre) : mg = GMm/R² soit GM= g₀R²

d'où : v₀² =2g₀Rh / (R+h) = 2*9,8*6,37 10⁶*3 10⁵ / 6,67 10⁶ =5,61 10⁶ ; v₀ =2370 m/s.

Ecrire la 2^è loi de Newton sur l'orbite circulaire : le satellite n'est soumis qu'à la force de gravitation de la Terre ;

l'accélération est centripète de valeur a_N= v₁² / (R+h) ; d'où GMm/(R+h)² = mv₁² / (R+h)

v₁² =GM/(R+h) = g₀R²/(R+h) ; v₁ =R[g₀/(R+h)]^½= 6,37 10⁶[9,8/6,67 10⁶]^½= 7721 m/s.

vitesse angulaire du satellite : w= v₁/(R+h) =1,16 10^-3 rad/s.

En réalité il est plus économique de :

- Traverser verticalement l'atmosphère afin de limiter le freinage ; puis en modifiant l'inclinaison de la poussée du moteur on gagne l'orbite elliptique de Holmann dont l'apogée est l'altitude h. Cette apogée étant atteinte, on met à nouveau en marche le moteur pour obtenir la vitesse v₁.

- Dans le cas de satellite géostationnaire, dont l'orbite est dans le plan équatorial, il vaut mieux que la base de lancement soit proche de l'équareur.

Le référentiel d'étude (r) est lié au satellite. Le centre de masse du satellite, noté A, est le centre de ce référentiel. L’objet est assimilé à un point matériel P de masse m repéré, dans (r), par les cordonnées du point P ( x, y, z) ; la position initiale P₀ a les coordonnées (x₀, y₀, z₀).

Forces s'exerçant, dans (r), sur l’objet P :

attraction terrestre :

force d'inertie :

HP projection de CP sur la plan de la trajectoire du satellite ( les grandeurs écrites en bleu et gras sont des vecteurs)

force de Coriolis : La force de Coriolis, perpendiculaire à la vitesse, ne travaille pas.

L’attraction terrestre dérive de l’énergie potentielle : -GMm/ CP= -m w² (R+h)³/ CP.

La force d'inertie centrifuge dérive de l’énergie potentielle : -½m w² HP².

Exprimons cette énergie potentielle au voisinage de A : x, y et z petits deavnt R+h.

on note, pour simplifier l'écriture r= R+h ; avec : HP² = (r+x)²+y²= r²+2rx+x²+y² et CP²=HP²+z²

d'où 1/CP= ( r²+2rx+x²+y² +z²)^-½ .

mettre r² en facteur commun : 1/CP= r^-1 ( 1+2x/r +(x²+y² +z²)/r²)^-½ .

Or : (1+e)^-½ voisin de 1-½e +3/8 e ² avec e = 2x/r +(x²+y² +z²)/r²

d'où : 1/CP= r^-1[1-x/r -(x²+y² +z²)/ (2r²) + 3x²/(2r²)]

soit E_P= -mw² r²[1-x/r -(x²+y² +z²)/ (2r²) + 3x²/(2r²)] -½mw²[ r²+2rx+x²+y²] = -mw²[1,5 r²-½(z²-3x²)].

Ecrire la seconde loi de Newton :

mx" = -dE_p/dx+F_{ic x} = 3mw²x + 2 mwy' (1)

my" = -dE_p/dy+F_{ic y} = -2 mwx' (2)

mz"= -dE_p/dz+F_{ic z} = -mw²z ou bien z" +w²z =0 (3)

La solution de (3) est du type : z= z₀ cos(wt+j) avec z(0)= z₀ et z'(0) = 0 ; or z'= -z₀ w sin (wt+j) d'où 0 = -z₀ w sin j--> j=0 ou p.

de plus z(0) = 0 = z₀ cos j donne j=0. z= z₀ cos(wt).

Par intégration (2) devient : y'= -2w(x-x₀) (4)

Repport dans (1) : x"=3w²x - 4 w²(x-x₀) soit x"+ w²x = 4 w²x₀.

la solution est du type : x= A cos(wt) + B sin(wt)+ 4x₀.

or x(0)=x₀ soit x₀ = A + 4x₀ --> A= -3x₀ ; de plus x'(0)=0 soit -Aw sin(0)+ Bw cos(0) =0 --> B=0 ; x= x₀ (4-3cos(wt)).

Repport dans (4) : y'= 6wx₀[ cos(wt)-1] ; par intégration y = 6x₀[ sin(wt)-wt] + y₀.

La trajectoire est une cycloïde décalée par rapport à l'axe des x. Initialement la force d'inertie centrifuge l'emporte ; puis la vitesse augmente et la force de Coriolis courbe la trajectoire.