Les points de Lagrange

Point de Lagrange

Pour les articles homonymes, voir Lagrange.

Un point de Lagrange (noté L₁ à L₅), ou, plus rarement, point de libration, est une position de l'espace où les champs de gravité de deux corps en orbite l'un autour de l'autre, et de masses substantielles, se combinent pour compenser exactement la force centrifuge. Dans le cas où les deux corps sont en orbite circulaire, ces points représentent les endroits où un troisième corps de masse négligeable resterait immobile par rapport aux deux autres, au sens où il accompagnerait à la même vitesse angulaire leur rotation autour de leur centre de gravité commun sans que sa position par rapport à eux n'évolue. Au nombre de cinq, ces points se scindent en deux points stables dénommés L₄ et L₅, et en trois points instables notés L₁ à L₃. Ils sont nommés en l'honneur du mathématicien français Joseph-Louis Lagrange. Ils interviennent dans l'étude de certaines configurations d'objets du Système solaire (principalement pour les points stables) et dans le placement de divers satellites artificiels (principalement pour les points instables).

Historique

En mécanique céleste, il est un sujet qui a passionné de nombreux mathématiciens : c'est le problème dit des trois corps. Newton, après avoir énoncé sa loi qui exprime que « les corps s'attirent avec une force proportionnelle au produit de leur masse et inversement proportionnelle au carré de la distance de leurs centres », a cherché à décrire le comportement de trois corps sans y parvenir. Il faut attendre le mathématicien Joseph-Louis Lagrange qui, en 1772, étudia le cas d'un petit corps, de masse négligeable, soumis à l'attraction de deux plus gros : le soleil et, par exemple, une planète. Il découvrit qu'il existait des positions d'équilibre pour le petit corps, des endroits où toutes les forces se compensent.

Définition

Un objet de faible masse situé exactement en ces points n'en bouge plus relativement aux deux autres corps, et tourne de concert avec eux (par exemple une planète et le Soleil).

Trois des points de Lagrange sont situés sur l'axe reliant les deux corps. Dans le cas d'un grande dissymétrie de masse entre ceux-ci, deux points sont situés proches et de part et d'autre du corps peu massif, alors que le troisième est quasiment situé à l'opposé du corps peu massif par rapport au corps massif.

Si on donne en exemple les points de Lagrange du système Soleil-Terre, ces cinq points sont notés et définis comme suit (échelle non respectée) :

L₁ : sur la ligne définie par les deux masses, entre celles-ci, la position exacte dépendant du rapport de masse entre les deux corps ; dans le cas où l'un des deux corps a une masse beaucoup plus faible que l'autre, le point L₁ est situé nettement plus près du corps peu massif que du corps massif.
L₂ : sur la ligne définie par les deux masses, au-delà de la plus petite. Dans le cas où l'un des deux corps a une masse beaucoup plus faible, la distance de L₂ à ce corps est comparable à celle entre L₁ et ce corps.
L₃ : sur la ligne définie par les deux masses, au-delà de la plus grande. Dans le cas où l'un des deux corps est notablement moins massif que l'autre, la distance entre L₃ et le corps massif est comparable avec celle entre les deux corps.

Les deux derniers points de Lagrange forment avec les deux corps des triangles équilatéraux.

L₄ et L₅ : sur les sommets des deux triangles équilatéraux dont la base est formée par les deux masses. Sans qu'il y ait de consensus précis, L₄ est celui des deux points en avance sur l'orbite de plus petite des masses, dans son orbite autour de la grande, et L₅ est en retard. Ces points sont parfois appelés points de Lagrange triangulaires ou points Troyens, du fait que c'est le lieu où se trouvent les astéroïdes troyens du système Soleil-Jupiter. Contrairement aux trois premiers points, ces deux derniers ne dépendent pas des masses relatives des deux autres corps.

Calcul de la position des points de Lagrange

Le calcul de la position des points de Lagrange se fait en considérant l'équilibre d'un corps de masse négligeable entre le potentiel gravitationnel créé par deux corps en orbite et la force centrifuge. La position des points L₄ et L₅ peut être obtenue analytiquement. Celle des trois autres points L₁ à L₃ s'obtient en résolvant numériquement ou éventuellement à l'aide d'un développement limité une équation algébrique. La position de ces trois points est donnée dans le tableau ci-dessous dans le cas où la masse d'un des deux corps (en l'occurrence le numéro 2) est négligeable devant l'autre, situé à une distance R du précédent. Les positions sont données le long de l'axe reliant les deux corps, dont l'origine est identifiée au centre de gravité du système, et dont l'orientation va du corps 1 au corps 2. Les quantités r₂ et q dénotent respectivement la position du corps 2 sur l'axe et le rapport de la masse du corps le plus léger à la masse totale des deux corps. Enfin, on utilise la quantité ε définie par ε = (q / 3)^1/3.

Point	Position par rapport au centre de gravité du système
L₁	$r_2 + R \left(- \epsilon + \frac{1}{3} \epsilon^2 + \frac{1}{9} \epsilon^3 + \mathcal{O}(\epsilon^4)\right)$
L₂	$r_2 + R \left(\epsilon + \frac{1}{3} \epsilon^2 - \frac{1}{9} \epsilon^3 + \mathcal{O}(\epsilon^4)\right)$
L₃	$- R + R\left(- \frac{5}{12} q + \frac{1127}{20736} q^3 + \mathcal{O}(q^4)\right)$

Dans la littérature, on trouve parfois des expression quelque peu différentes, du fait que l'origine de l'axe est prise ailleurs que sur le centre de gravité, et que l'on utilise comme terme à la base du développement limité le rapport entre les deux masses plutôt que le rapport de la plus petite à la masse totale, c'est-à-dire que l'on utilise parfois la quantité q' définie par $q' = \frac{m_2}{M_1} = \frac{q}{1 - q}$ .

Détail du calcul — Introduction

Préliminaires

On note M₁ et m₂ la masse deux deux corps, la masse du premier étant supposée supérieure ou égale à celle du second. Les deux corps sont supposés être en orbite circulaire, leur séparation étant R. Les deux corps orbitent autour de leur centre de gravité commun. On note r₁ et r₂ les distances algébriques des deux corps par rapport à un axe orienté du corps 1 au corps 2 (c'est-à-dire que r₁ va être négatif et r₂ positif). Le centre de gravité est défini par l'équation

M 1 r 1 + m 2 r 2 = 0

avec par définition, de la distance R,

r 2 - r 1 = R

Ces deux équations ont pour solution

$r_1 = - \frac{m_2}{M} R$ ,

$r_2 = \frac{M_1}{M} R$ ,

où on a noté M = M₁ + m₂ la masse totale du système.

Les deux corps orbitent l'un autour de l'autre à une vitesse angulaire ω, dont la valeur est donnée par la troisième loi de Kepler :

$\omega^2 = \frac{G M}{R^3}$ ,

G étant la constante de gravitation.

Si l'on se place dans le référentiel tournant avec les deux corps, c'est-à-dire à la vitesse angulaire ω, un corps immobile sera soumis, outre aux forces gravitationnelles des deux corps, à la force centrifuge. Si on notre r le rayon vecteur de ce corps, la force centrifuge par unité de masse f_c à laquelle il sera soumis s'écrit

${\boldsymbol f}_{\rm c} = {\boldsymbol r} \; \omega^2$ .

Équation fondamentale

La définition d'un point de Lagrange est que la somme des forces gravitationnelles et inertielles s'annulent en ces points. En notant r le rayon vecteur du ou des points en question, on a ainsi

$- \frac{G M_1 ({\boldsymbol r} - {\boldsymbol r}_1)}{||{\boldsymbol r} - {\boldsymbol r}_1||^3} - \frac{G m_2 ({\boldsymbol r} - {\boldsymbol r}_2)}{||{\boldsymbol r} - {\boldsymbol r}_2||^3} + {\boldsymbol r} \; \omega^2 = 0$ ,

les doubles barres indiquant que l'on prend la norme des vecteurs considérés. On remplace ensuite la vitesse angulaire ω par sa valeur issue de la troisième loi de Kepler, ce qui donne

que l'on simplifie immédiatement par la constante de gravitation

$- M_1 \frac{{\boldsymbol r} - {\boldsymbol r}_1}{||{\boldsymbol r} - {\boldsymbol r}_1||^3} - m_2 \frac{{\boldsymbol r} - {\boldsymbol r}_2}{||{\boldsymbol r} - {\boldsymbol r}_2||^3} + M \frac{{\boldsymbol r}}{R^3} = 0$ .

C'est la résolution de cette équation qui donne les différents points de Lagrange.

Les deux cas à considérer

La projection de cette équation perpendiculairement au plan de l'orbite, dont la normale est donnée par un vecteur noté $\hat{\boldsymbol z}$ donne immédiatement

${\boldsymbol r} \cdot \hat{\boldsymbol z} = 0$ ,

ce qui implique que l'ensemble des points de Lagrange est situé dans le plan de l'orbite. La résolution de l'équation se fait donc dans le plan orbital. Deux cas sont à considérer :

celui où l'on cherche un point le long de l'axe formé par les deux corps,
celui où l'on cherche un point en dehors de cet axe.

Le second cas s'avère être le plus simple à étudier.

Détail du calcul — Les points L₄ et L₅

Cas des points L₄ et L₅

On suppose que le rayon vecteur r n'est pas parallèle à l'axe passant par les deux corps. On projette donc l'équation fondamentale perpendiculairement à cet axe, direction que l'on suppose définie par un vecteur noté $\hat{\boldsymbol y}$ . Par définition, cette direction étant perpendiculaire à l'axe reliant les deux corps, on a

$\hat {\boldsymbol y} \cdot {\boldsymbol r}_1 = \hat {\boldsymbol y} \cdot {\boldsymbol r}_2 = 0.$ .

L'équation fondamentale se réécrit donc

$- M_1 \frac{\hat {\boldsymbol y} \cdot {\boldsymbol r}}{||{\boldsymbol r} - {\boldsymbol r}_1||^3} - m_2 \frac{\hat {\boldsymbol y} \cdot {\boldsymbol r}}{||{\boldsymbol r} - {\boldsymbol r}_2||^3} + M \frac{\hat {\boldsymbol y} \cdot {\boldsymbol r}}{R^3} = 0$ .

Les termes en $\hat {\boldsymbol y} \cdot {\boldsymbol r}$ se simplifient, ce qui donne

$- \frac{M_1}{||{\boldsymbol r} - {\boldsymbol r}_1||^3} - \frac{m_2}{||{\boldsymbol r} - {\boldsymbol r}_2||^3} + \frac{M}{R^3} = 0$ .

On définit maintenant la direction $\hat {\boldsymbol s}$ comme la perpendiculaire à r. Comme r n'est pas colinéaire à r₁ et r₂, les quantités $\hat {\boldsymbol s} \cdot {\boldsymbol r}_{1,2}$ ne sont pas nulles. En projetant l'équation fondamentale le long de s, on obtient

$M_1 \frac{\hat {\boldsymbol s} \cdot {\boldsymbol r}_1}{||{\boldsymbol r} - {\boldsymbol r}_1||^3} + m_2 \frac{\hat {\boldsymbol s} \cdot {\boldsymbol r}_2}{||{\boldsymbol r} - {\boldsymbol r}_2||^3} = 0$ .

Or, d'après le théorème de Thalès, les projections de r₁ et r₂ le long de $\hat {\boldsymbol s}$ sont dans le même rapport que les projections de ces vecteurs le long de l'axe reliant les deux corps. Il s'ensuit que l'équation précédente peut se réécrire

$M_1 \frac{r_1}{||{\boldsymbol r} - {\boldsymbol r}_1||^3} + m_2 \frac{r_2}{||{\boldsymbol r} - {\boldsymbol r}_2||^3} = 0$ .

Le barycentre des deux corps implique, comme vu précédemment, que

M 1 r 1 + m 2 r 2 = 0

La combinaison de cette équation et celle qui précède implique donc que les deux distances $||{\boldsymbol r} - {\boldsymbol r}_1||$ et $||{\boldsymbol r} - {\boldsymbol r}_2||$ sont identiques, leur valeur étant notée R' :

$||{\boldsymbol r} - {\boldsymbol r}_1|| = ||{\boldsymbol r} - {\boldsymbol r}_2|| \equiv R'$ .

En injectant ce résultat sur la projection le long de r, il vient alors

$- \frac{M_1}{R'^3} - \frac{m_2}{R'^3} + \frac{M}{R^3} = 0$ .

En multipliant le tout par R'³ et en se souvenant que M est la somme des deux masses, on obtient finalement

$M \left(\frac{R'^3}{R^3} - 1 \right) = 0$ ,

ce qui donne au final

$||{\boldsymbol r} - {\boldsymbol r}_1|| = ||{\boldsymbol r} - {\boldsymbol r}_2|| = R$ ,

c'est-à-dire que les points cherchés forment un triangle équilatéral avec les deux corps du système. Ces triangles sont de plus inclus dans le plan orbital, ce qui donne deux points possibles, notés comme annoncé L₄ et L₅, étant situé de part et d'autre de l'axe reliant les deux corps.

En utilisant le théorème de Pythagore, la distance D de ces deux points de Lagrange du centre de gravité du système s'écrit

$D^2 = (-|r_1| + R / 2)^2 + (\sqrt{3}R / 2)^2$ ,

ce qui donne

$D^2 = r_1^2 + \frac{1}{4} R^2 - R |r_1| + \frac{3}{4} R^2$ ,

ce qui donne

$D^2 = R^2 + r_1^2 - R |r_1|$ .

En utilisant le fait que $R = | r 1 | + r 2$ , il vient

$D^2 = r_2^2 + |r_1| R = r_1^2 + r_2 R$ .

La distance est donc supérieure aux distances de chacun des deux corps au centre de gravité du système. Ces points de Lagrange sont donc au-delà de l'orbite du corps le moins massif et ne sont pas strictement situés sur celle-ci, quoique que ce soit quasiment le cas dans la limite où la masse du corps le plus léger devient négligeable par rapport à celle de son compagnon.

Détail du calcul — Les points L₁ à L₃

Cas des points L₁ à L₃

Dans le cas où l'on considère des points de Lagrange situés sur l'axe reliant les deux corps, trois sous-cas sont à considérer :

Le cas où le ou les points sont entre les corps 1 et 2 ;
Le cas où le ou les points sont à l'opposé du corps 2 par rapport au corps 1 ;
Le cas où le ou les points sont à l'opposé du corps 1 par rapport au corps 2 .

Dans ces trois cas, l'équation fondamentale se réécrit de la façon suivante :

Cas 1 : la projection des deux forces sur l'axe a des signes opposés (la projection de la force exercée par le corps 1 est négative, celle de la force exercée par le corps 2 est positive), ce qui donne

$- \frac{M_1}{(r - r_1)^2} + \frac{m_2}{(r - r_2)^2} + M \frac{r}{R^3} = 0$ ,

avec

r 1 < r < r 2

Cas 2 : la projection des deux forces sur l'axe est négative, ce qui donne

$- \frac{M_1}{(r - r_1)^2} - \frac{m_2}{(r - r_2)^2} + M \frac{r}{R^3} = 0$ ,

avec

r 1 < r 2 < r

Cas 3 : la projection des deux forces sur l'axe est positive, ce qui donne

$\frac{M_1}{(r - r_1)^2} + \frac{m_2}{(r - r_2)^2} + M \frac{r}{R^3} = 0$ ,

avec

r < r 1 < r 2

Chacune de ces trois équations peut se ramener à une équation polynomiale du cinquième degré, pour laquelle il n'existe pas de solution analytique exacte, sauf cas particulier (comme celui des deux masses identiques par exemple).

L'unicité des solutions dans chacun des trois cas se déduit du fait que l'équation à résoudre sur l'équilibre des forces dérive d'un potentiel U, donné par

$U = -\frac{M_1}{||{\boldsymbol r} - {\boldsymbol r}_1||} - \frac{m_2}{||{\boldsymbol r} - {\boldsymbol r}_2||} - \frac{1}{2} \frac{M r^2}{R^3}$ .

Ce potentiel représente des pôles en r₁ et r₂, et correspond en dehors de ces valeurs à la somme de trois termes concaves et est donc localement concave. Il ne possède donc qu'un seul extréme local dans chacun des domaines où il est défini, c'est-à-dire dans chacun des trois cas cités plus haut.

Solutions pour L₁ à L₃ dans le cas où le rapport entre les masses est faible

Forme réduite et solution dans le cas où le rapport entre les masses est faible

Quand le rapport entre m₂ et M₁ (ou entre m₂ et M) est faible, on peut trouver une solution approchée pour la position de chacun des points en effectuant un développement limité à partir d'une solution approchée facile à trouver. Pour simplifier les notations, on effectue un changement d'échelle afin d'exprimer toutes les longueurs en unité de la séparation R et les masse en unité de la masse totale M. On pose ainsi

u = r / R

u 1,2 = r 1,2 / R

et on définit le petit paramètre q par

$q \equiv m_2 / M$ ,

à partir de quoi on peut exprimer

M 1 / M = 1 - q

u 1 = q

u 2 = 1 - q

Dans ce cas là, les trois équations écrites ci-dessus prennent la forme plus simple

Cas 1 :

$- \frac{1 - q}{(u + q)^2} + \frac{q}{(u - 1 + q)^2} + u = 0$ ,

avec

- q < u < 1 - q

Cas 2 : la projection des deux forces sur l'axe est négative, ce qui donne

$- \frac{1 - q}{(u + q)^2} - \frac{q}{(u - 1 + q)^2} + u = 0$ ,

avec

1 - q < u

Cas 3 : la projection des deux forces sur l'axe est positive, ce qui donne

$\frac{1 - q}{(u + q)^2} + \frac{q}{(u - 1 + q)^2} + u = 0$ ,

avec

u < - q

Le point L₁

Quand la masse du corps 2 est négligeable, l'attraction de celui-ci est négligeable sauf si la particule d'épreuve est très proche. Or, quand l'attraction du corps 2 est négligeable, l'équilibre entre l'attraction du corps 1 et la force centrifuge est tel que la distance du point d'équilibre est de l'ordre de R. Quand le point d'équilibre est situé à l'opposé du corps 2, on est dans le cas du point de Lagrange L₃, qui est donc, en gros, situé à l'opposé du corps 2 par rapport au corps 1. Dans le cas contraire, on va donc supposer que le point d'équilibre est plutôt proche du corps 2 (et donc à nouveau situé à la distance R du corps 1), mais néanmoins suffisamment éloigné pour que l'attraction du corps 2 exercée sur la particule d'épreuve reste petite par rapport à celle du corps 1. On pose donc à partir de la forme réduite

u = u 2 + ε' = 1 - q + ε

où ici ε' est une quantité petite et négative (on suppose ici que le point est entre les deux corps). L'équation réduite se transforme alors en

$- \frac{1 - q}{(1 + \epsilon')^2} + \frac{q}{\epsilon'^2} + 1 - q + \epsilon' = 0$ .

On effectue un développement limité au premier ordre de l'attraction produite par le corps 1 :

$- (1 - q) (1 - 2 \epsilon') + \frac{q}{\epsilon'^2} + 1 - q + \epsilon' = 0$ .

Les termes en 1 - q se simplifient, et il reste

$(3 - 2 q) \epsilon' + \frac{q}{\epsilon'^2} = 0$ .

Toujours en ne gardant que les termes d'ordre le plus bas en q, il vient

$\epsilon' = - \left(\frac{q}{3}\right)^\frac{1}{3}$ .

On peut par la suite continuer le calcul, en développement l'écart du point au corps 2 en puissances de ε'. On pose ainsi

$u = 1 - q + \epsilon' + x \epsilon'^2 + \mathcal{O} (\epsilon'^3)$ .

L'équation fondamentale réduite donne alors

$- \frac{1 - q}{(1 + \epsilon' + x \epsilon'^2)^2} + \frac{q}{(\epsilon' + x \epsilon'^2)^2} + 1 - q + \epsilon' + x \epsilon'^2 = 0$ .

On peut factoriser le second terme avec q / ε'², que l'on peut remplacer par sa valeur, soit -3 ε'. On obtient alors

- (1 - q)(1 + ε' + x ε' 2) - 2 - 3ε'(1 + x ε') - 2 + 1 - q + ε' + x ε' 2 = 0

On effectue ensuite un développement limité des deux premiers termes, au second ordre pour le premier et au premier ordre pour le suivant, ce qui donne

- (1 - q)(1 - 2ε' - 2 x ε' 2 + 3ε' 2) - 3ε'(1 - 2 x ε') + 1 - q + ε' + x ε' 2 = 0

d'où on déduit que x vaut un tiers, ce qui donne

$u = 1 - q + \epsilon' + \frac{1}{3} \epsilon'^2 + \mathcal{O}(\epsilon^3)$ .

Le développement peut ensuite être continué suivant la même procédure. À l'ordre suivant, on a ainsi

$u = 1 - q + \epsilon' + \frac{1}{3} \epsilon''^2 - \frac{1}{9} \epsilon'^3 + \mathcal{O}(\epsilon'^4)$ .

Le point L₂

Le cas du point L₂ se résout exactement comme dans la section précédente, si ce n'est que le signe du second terme de l'équation fondamentale est négatif. On pose donc

u = 1 - q + ε

ε étant cette fois-ci supposé petit et positif, et on a ainsi

$- \frac{1 - q}{(1 + \epsilon)^2} - \frac{q}{\epsilon^2} + 1 - q + \epsilon = 0$ .

La résolution à l'ordre le plus bas donne

$- (1 - q)(1 - 2 \epsilon) - \frac{q}{\epsilon^2} + 1 - q + \epsilon = 0$ ,

qui après annulation des termes donne

$3 \epsilon = \frac{q}{\epsilon^2}$ ,

soit

$\epsilon = \left(\frac{q}{3}\right)^\frac{1}{3} = \epsilon'$ .

Cela correspond au signe près au même résultat que précédemment. La suite du développement de la solution se fait de même que précédemment. On part de

u = 1 - q + ε + x ε 2

et on injecte ce résultat dans l'équation fondamentale

$- \frac{1 - q}{(1 + \epsilon + x \epsilon^2)^2} - \frac{q}{(\epsilon + x \epsilon^2)^2} + 1 - q + \epsilon + x \epsilon^2 = 0$ .

Comme précédemment, on transforme cette expression selon

- (1 - q)(1 - 2ε - 2 x ε 2 + 3ε 2) - 3ε(1 - 2 x ε) + 1 - q + ε + x ε 2 = 0

ce que l'on résout en

$x = \frac{1}{3}$ ,

soit

$u = 1 - q + \epsilon + \frac{1}{3} \epsilon^2 + \mathcal{O} (\epsilon^3)$ .

Cette expression est identique à celle du premier point de Lagrange en remplaçant ε' par ε, mais ces deux points sont dissymétriques : comme le signe de ε, ε' change entre le point L₁ et le point L₂, la correction du second ordre, toujours positive, rapproche le point L₁ du corps 2 alors qu'elle éloigne le point L₂ : les deux points ne sont plus à égale distance du corps 2. Pour la Terre, le rapport de masse est de 1 / 300 000, et ε est de l'ordre 0,01, ce qui place les deux points par rapport à la Terre à une distance d'environ un centième de la distance Terre-Soleil, soit dans les 1 500 kilomètres. Le terme de second ordre est de l'ordre d'un trente-millième de la distance Terre-Soleil, soit dans les 5 000 km. Le point L₁ est donc environ 10 000 km plus près de la Terre que ne l'est L₂.

Enfin, on peut poursuivre le développement à l'ordre supérieur, ce qui donne, tous calculs faits

$u = 1 - q + \epsilon + \frac{1}{3} \epsilon^2 - \frac{1}{9} \epsilon^3 + \mathcal{O} (\epsilon^4)$ .

Le point L₃

Dans le cas 3, qui va correspondre au point L₃, l'équation fondamentale s'écrit

$\frac{1 - q}{(u + q)^2} + \frac{q}{(u - 1 + q)^2} + u = 0$ .

Comme le point est suppose au-delà du corps 1 par rapport au corps 2, il est plus proche du corps le plus massif, dont l'attraction va être prépondérante par rapport à l'autre corps. Dans la situation où l'on se place, le point recherché a donc a position approximée par

$\frac{1 - q}{(u + q)^2} + u = 0$ .

La solution approchée de cette équation est bien sûr

$u \simeq - 1$ .

Pour trouver les écarts à cette valeurs, on écrit dans l'équation fondamentale

u = - 1 + v

et on résout l'équation en prenant en compte les premiers termes en q. On obtient ainsi

$\frac{1 - q}{(-1 + v + q)^2} + \frac{q}{(-1 - 1 + v + q)^2} - 1 + v = 0$ .

Les quantités v et q étant petites devant R, le premier terme s'écrit

$\frac{1 - q}{(-1 + v + q)^2} = (1 - q) (1 - v - q)^{-2} \simeq (1 - q) + 2 (v + q)$ .

Le second terme étant négligeable par rapport au précédent (il est proportionnel à q), il peut s'approximer en

$\frac{q}{(-1 - 1 + v + q)^2} \simeq \frac{q}{4}$ .

En combinant l'ensemble de ces termes, on obtient

$1 - q + 2 (v + q) + \frac{q}{4} - 1 + v = 0$ ,

ce qui donne

$\frac{5}{4} q + 3 v = 0$ ,

c'est-à-dire

$v = - \frac{5}{12} \frac{m_2}{M} + \mathcal{O}\left(\frac{m_2^2}{M^2}\right)$ .

On peut sans difficulté continuer ce calcul en posant désormais

$u = -1 - \frac{5}{12} q + w$ ,

w étant cette fois proportionnel à q². L'équation fondamentale devient alors

$\frac{1 - q}{\left(-1 - \frac{5}{12} q + w + q\right)^2} + \frac{q}{\left(-1 - 1 - \frac{5}{12} q + w + q\right)^2} - 1 - \frac{5}{12} q + w = 0$ ,

c'est-à-dire

$\frac{1 - q}{\left(1 - \frac{7}{12} q - w\right)^2} + \frac{q}{\left(2 - \frac{7}{12} q - w\right)^2} - 1 - \frac{5}{12} q + w = 0$ .

En développant cette expression au second ordre en q, on trouve

$w = 0 + \mathcal{O}(q^3)$ ,

c'est-à-dire que w est au plus en q³. En refaisant le calcul dans ce cadre là, on trouve finalement

$w = \frac{1127}{20736} q^3 + \mathcal{O}(q^4)$ .

Il est rarement utile de pousser le calcul jusque là : dans une configuration Soleil-Planète, le dernier terme correctif est au mieux de l'ordre !)^-9, puisque le plus gros rapport de masse Planète-Soleil, dans le cas de Jupiter est de l'ordre d'un millième. Le terme q³ est donc, pour Jupiter, de l'ordre d'un milliardième, ce qui, au vu de la taille de son orbite, correspond à une correction d'une cinquantaine de mètre étant donné que la fraction en facteur de q³ est de l'ordre d'un vingtième. Pour le système Terre-Soleil (distance de l'ordre de 150 millions de kilomètres, rapport de masse de l'ordre de 1 / 300 000), la dernière correction est une fraction de micron !

Stabilité

Le calcul ci-dessus n'indique en rien si les points de Lagrange sont stables. La stabilité ou non de ces points est du reste peu intuitive. Dans le référentiel tournant avec les deux corps, une particule d'épreuve peut être vue comme soumise à un potentiel incluant la contribution gravitationnelle et celle de la force centrifuge. Ce potentiel, noté Ω, s'écrit ainsi

$\Omega = - \frac{G M_1}{|\boldsymbol{r} - \boldsymbol{r}_1|} - \frac{G m_2}{|\boldsymbol{r} - \boldsymbol{r}_2|} - \frac{1}{2} |\boldsymbol{r}|^2 \omega^2$ .

Tous les termes de ce potentiel sont négatifs et décroissent à mesure que l'on s'éloigne des masses (pour les deux premiers termes) ou du centre de gravité du système (pour le second). On peut ainsi montrer que les points de Lagrange L₄ et L₅ sont des maxima locaux du potentiel Ω (voir ci-dessous) et que les trois autres points sont des points selles. D'ordinaire, une position d'équilibre (déterminée par l'annulation des dérivées du potentiel) est stable uniquement si on se situe dans des minima locaux du potentiel. Cependant, étant donné que l'on est dans un référentiel tournant, le référentiel est non inertiel. Un objet se déplaçant dans ce référentiel, par exemple au voisinage d'une position d'équilibre, va être soumis à la force de Coriolis, et son mouvement ne dépend pas uniquement de la forme du potentiel. Pour étudier la stabilité des points de Lagrange, il faut donc tenir compte de la force de Coriolis.

Pour calculer la stabilité des points de Lagrange, il faut ainsi étudier l'équation du mouvement d'un objet situé au voisinage d'un de ces points. En notant δR le vecteur de coordonnées δX et δY donnant l'écart d'un tel objet à un des points de Lagrange (que l'on suppose confiné au plan orbital, l'équation du mouvement s'écrit

$\ddot{\boldsymbol{\delta R}} = 2 \boldsymbol{\omega} \wedge \dot {\boldsymbol {\delta R}} + \boldsymbol {\delta f}$ ,

où δf représente la force par unité de masse exercée sur l'objet. Cette force est petite du fait qu'au point de Lagrange la force (constituée d'une composante gravitationnelle et de la force centrifuge) est nulle et que l'on se place à proximité d'un tel point. Cette force peut se calculer en termes d'un développement limité. Par exemple, pour la composante X, on a

$\delta f_X = f_X (0) + \frac{\partial f_X}{\partial X} \delta X + \frac{\partial f_X}{\partial Y} \delta Y$ .

Le premier terme correspond à la force s'exerçant au point de Lagrange, force qui est nulle par construction. Par ailleurs, la force dérivant d'un potentiel, on peut exprimer les dérivées de la force en termes de dérivées secondes du potentiel :

$\delta f_X = - \frac{\partial^2 \Omega}{\partial X^2} \delta X - \frac{\partial^2 \Omega}{\partial X \partial Y} \delta Y$ .

On peut ainsi exprimer l'équation du mouvement en termes des composantes selon

$\ddot{\delta X} = 2 \omega \dot{\delta Y} - \frac{\partial^2 \Omega}{\partial X^2} \delta X - \frac{\partial^2 \Omega}{\partial X \partial Y} \delta Y$ ,

$\ddot{\delta Y} = - 2 \omega \dot{\delta X}- \frac{\partial^2 \Omega}{\partial X \partial Y} \delta X - \frac{\partial^2 \Omega}{\partial Y^2} \delta Y$ .

Ce groupe d'équation peut être mis sous la forme d'un système de quatre équations différentielles du premier ordre :

$\frac{\partial}{\partial t}\left(\begin{array}{c} \delta X \\ \delta Y \\ \dot{\delta X} \\ \dot{\delta Y} \end{array} \right) = \left(\begin{array}{cccc} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ -\Omega_{,xx} & - \Omega_{,xy} & 0 & -2 \omega \\ - \Omega_{,xy} & -\Omega_{,yy} & 2 \omega & 0 \end{array}\right) \left(\begin{array}{c} \delta X \\ \delta Y \\ \dot{\delta X} \\ \dot{\delta Y} \end{array} \right)$ ,

où les dérivées partielles du potentiel Ω ont été notée en indice précédé d'un virgule (par exemple, Ω_,xx correspond à $\partial^2 \Omega / \partial x^2$ ).

La stabilité du point de Lagrange considéré est obtenu en recherchant les solutions de cette équation. Pour cela, il suffit de trouver des solutions de type exponentielle, en $e Γ t$ . On va ainsi procéder à la diagonalisation de la matrice ci-dessus, que l'on notera A. Les valeurs propres trouvées vont correspondre aux quantités Γ ci-dessus, les écarts à la position d'équilibre étant alors une certaine combinaison d'au plus quatre exponentielles. La stabilité du système est assuré par le fait que les exponentielles ne croissent pas au cours du temps, c'est-à-dire que les quantités Γ sont soit négatives, soit complexes à partie réelle négatives. En fait, il n'est pas nécessaire de diagonaliser complètement la matrice, il sufft d'en trouver les valeurs propres, c'est-à-dire les solutions de l'équation

det(A - λ I) = 0

Ce déterminant s'écrit

$\left|\begin{array}{cccc} - \lambda & 0 & 1 & 0 \\ 0 & -\lambda & 0 & 1 \\ -\Omega_{,xx} & - \Omega_{,xy}, & - \lambda & - 2 \omega \\ - \Omega_{,xy} & - \Omega_{yy} & 2 \omega & -\lambda \end{array} \right|$ ,

et il vaut

$\lambda^4 + (\Omega_{,xx} + \Omega_{,yy} + 4 \omega^2) \lambda^2 + \Omega_{,xx} \Omega_{,yy} - \Omega_{,xy}^2$ .

Cette équation peut se ramener à une équation polynomiale du second ordre en λ². Les solutions de l'équation de départ sont donc deux couples de nombres opposés deux à deux. Pour que deux nombres opposés soient négatifs ou nuls ou alors de partie réelle négtive ou nulle, il faut obligatoirement qu'ils soient des nombres imaginaires purs, donc que les solutions de l'équation en λ² soient réelles et négatives. Pour que ces solutions soient réelles, il faut donc que le discriminant soit positif, soit ici

$(\Omega_{,xx} + \Omega_{,yy} + 4 \omega^2)^2 - 4 (\Omega_{,xx} \Omega_{,yy} - \Omega_{,xy}^2) > 0$ .

Une fois ceci obtenu, il faut que les deux solutions réelles soient négative, ce qui implique que simultanément leur somme noit négative et leur produit positif, ce qui implique

Ω, x x + Ω, y y + 4ω 2 > 0

$\Omega_{,xx} \Omega_{,yy} - \Omega_{,xy}^2 > 0$ .

La stabilité d'un point de Lagrange est soumise à la réalisation de ces trois contraintes. Parmi ces constraintes, la dernière a une interprétation simple : le signe de la quantité $\Omega_{,xx} \Omega_{,yy} - \Omega_{,xy}^2$ détermine si la position considérée est un extrémum local ou un point selle. En l'occurence, la positivité de cette quantité implique qu'elle doive être un extrémum local, condition nécessaire mais non suffisante à la stabilité du point de Lagrange. Quand cette quantité est négative, on a un point selle et le point de Lagrange est instable. Par contre, de façon plus surprenante, un point de Lagrange peut être stable s'il correspond à un maximum local du potentiel, c'est-à-dire que Ω_,xx + Ω_,yy peut être négatif, pourvu que cette quantité ne dépasse pas la valeur critique de -4 ω². En pratique, c'est ce qui se produit dans certains cas pour les points de Lagrange L₄ et L₅. L'interprétation physique de cete situation est que la stabilité est alors assurée par la force de Coriolis. Un objet légèrement décalé d'un tel point va s'en éloigner dans un premier temps en façon radiale, avant de voir sa trajectoire incurvée par la force de Coriolis. Si le potentiel est partout décroissant autour du point, alors il est possible que la force de Coriolis force l'objet à tourner autour du point de Lagrange, à l'instar des nuages dans une dépression qui ne se dirigent pas vers le cœur de la dépression, mais sont contraints à une trajectoire circulaire autour de celui-ci.

Suite du calcul

Préliminaire

Pour étudier la stabilité des points de Lagrange, on va être amené à calculer les dérivées successives du potentiel. Ce potentiel fait intervenir la distance |r - r₁|. Il faut donc connaître les dérivées des différentes puissance d'une telle quantité. En coordonnées cartésiennes, cette quantité s'écrit

$|{\boldsymbol{r}} - {\boldsymbol{r}}_1| = \sqrt{(x-x_1)^2 + (y - y_1)^2 + (z - z_1)^2}$ .

Sa dérivée par rapport à l'une des coordonnées x, y, z, collectivement notées x_i s'écrit donc

$\partial_i |{\boldsymbol{r}} - {\boldsymbol{r}}_1| = \frac{1}{2} \frac{2 (x_i - x_{1\;i})}{\sqrt{(x-x_1)^2 + (y - y_1)^2 + (z - z_1)^2}} = \frac{x_i - x_{1\;i}}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|}$ .

La dérivée d'une puissance p quelconque de cette grandeur est donc

$\partial_i (|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^p) = p (x_i - x_{1\;i}) |{\boldsymbol{r}} - {\boldsymbol{r}}_1|^{p - 2}$ .

En adaptant ce résultat aux dérivées secondes des quantités intervenant dans le potentiel, on a

$\partial_{ij} \left(\frac{1}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|}\right) = - \frac{\delta_{ij}}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^3} + 3 \frac{(x_i - x_{1\;i}) (x_j - x_{1\;j})}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^5}$ ,

ce qui pour le potentiel complet donne

$\partial_{ij} \Omega = \delta_{ij}\frac{G M_1}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^3} - 3 \frac{(x_i - x_{1\;i}) (x_j - x_{1\;j}) G M_1}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^5} + \delta_{ij}\frac{G m_2}{|{\boldsymbol{r}} - {\boldsymbol{r}}_2|^3} - 3 \frac{(x_i - x_{2\;i}) (x_j - x_{2\;j}) G m_2}{|{\boldsymbol{r}} - {\boldsymbol{r}}_2|^5} - \omega^2 \delta_{ij}$ ,

où δ_ij représente le symbole de Kronecker. C'est la valeur de ces dérivées partielles qu'il faut calculer pour déterminer la stabilité des différents points de Lagrange. C'est pour les points de Lagrange L₄ et L₅ que ce calcul est le plus simple.

Cas des points de Lagrange L₄ et L₅

Ces points se caractérisent par le fait que leur distance aux deux corps est identique et égale à R :

$|{\boldsymbol{r}} - {\boldsymbol{r}}_1| = |{\boldsymbol{r}} - {\boldsymbol{r}}_2| = R$ .

Par ailleurs, on peut utiliser la troisième loi de Kepler pour passer des quantités du type G M / R³ à ω, et on connaît les coordonnées exactes des points de Lagrange. En évaluant les dérivées du potentiels au points de Lagrange L₄ ou L₅, on a

$x_i - x_{1\;i} = \left(\frac{1}{2} R, \pm \frac{\sqrt{3}}{2} R \right)$ ,

$x_i - x_{2\;i} = \left(- \frac{1}{2} R, \pm \frac{\sqrt{3}}{2} R \right)$ ,

le signe + s'appliquant pour L₅ et le signe - pour L₄. Au final, la matrice recherchée a pour composantes

$\partial_{ij} \Omega = \omega^2 \left(\begin{array}{cc} -\frac{3}{4} & \mp \frac{3\sqrt{3}}{4}(1-2q) \\ \mp \frac{3\sqrt{3}}{4}(1-2q) & -\frac{9}{4}\end{array} \right)$ .

Le déterminant de cette matrice vaut

$\det(\partial_{ij} \Omega) = \Omega_{,xx} \Omega_{,xy} - \Omega_{,xy}^2 = \frac{27}{16} \left(1 - (1 - 2 q)^2 \right) \omega^2 = \frac{27}{4} \omega^2 q (1 - q)$ ,

qui est toujours positif puisque q est confiné entre 0 et 1. Cette première condition de stabilité est établie. La seconde condition de stabilité s'écrit

$\Omega_{,xx} + \Omega_{,yy} + 4 \omega^2 = \omega^2\left(-\frac{3}{4} - \frac{9}{4} + 4 \right) = + \omega^2$ ,

quantité là encore positive. Enfin, le discriminant donne

$\left(\Omega_{,xx} + \Omega_{,yy} + 4 \omega^2\right)^2 - 4 \left(\Omega_{,xx} \Omega_{,xy} - \Omega_{,xy}^2\right) = \omega^2 \left(1 - 27 q (1 - q) \right) = \omega^2(27 q^2 - 27 q + 1)$ .

La stabilité des deux points est, au final, déterminée par la positivité de la quantité $27 q 2 - 27 q + 1$ . Les zéros q_a, q_b de ce polynôme sont donnés par la formule usuelle, qui ici indique

$q_{a,b} = \frac{1}{2} \pm \frac{1}{2} \sqrt{\frac{23}{27}}$ .

Ce polynôme a ainsi des valeurs négatives sur la plage $q \in \left]\frac{1}{2} - \frac{1}{2} \sqrt{\frac{23}{27}}, \frac{1}{2} + \frac{1}{2} \sqrt{\frac{23}{27}} \right[ \simeq ]0,\!03852, 0,\!96148[$ . Ainsi, la stabilité de ces deux points de Lagrange n'est assurée que si la plus petite masse n'excède pas 3,852% de la masse totale, ou, de façon équivalente, que le rapport des deux masses n'excède pas 4,006%^[1].

Cette condition est vérifiée pour toutes les configurations de type Soleil-Planète (où q n'excède pas environ un millième pour Jupiter), ou pour le système Terre-Lune (où q est de l'ordre de 1/80, soit 1,25%).

Cas des points de Lagrange L₁ à L₃

Les trois points de Lagrange L₁ à L₃ sont situés sur l'axe reliant les deux corps. Dans la formule donne les dérivées secondes, les quantités y_i - y_{1 i} sont nulles, alors que leurs analogues en x s'identifient aux distances entre un des corps et le point de Lagrange considéré. En conséquence, la matrice des dérivées secondes s'écrit

$\partial_{ij} \Omega = \left(\begin{array}{cc} - 2\frac{G M_1}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^3} - 2\frac{G m_2}{|{\boldsymbol{r}} - {\boldsymbol{r}}_2|^3} - \omega^2 & 0 \\ 0 & \frac{G M_1}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^3} + \frac{G m_2}{|{\boldsymbol{r}} - {\boldsymbol{r}}_2|^3} - \omega^2 \end{array} \right)$ .

Le terme Ω_,xx est manifestement négatif. Le signe du déterminant de la matrice est déterminé par celui de Ω_,yy : si ce dernier est positif, alors le point de Lagrange est un point selle et il est instable. On peut réécrire ce terme en utilisant la troisième loi de Kepler :

$\Omega_{,yy} = \omega^2 \left((1 - q)\frac{R^3}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^3} + q \frac{R^3}{|{\boldsymbol{r}} - {\boldsymbol{r}}_2|^3} - 1 \right)$ .

Le cas de L₁

Le point de Lagrange L₁ est situé entre les deux corps. Sa distance à ceux-ci, |r - r₁| et |r - r₂| est donc, à chaque fois, strictement inférieure à R. On a ainsi

$\left.\Omega_{,yy}\right|_{L_2} = \omega^2 \left((1 - q)\frac{R^3}{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|^3} + q \frac{R^3}{|{\boldsymbol{r}} - {\boldsymbol{r}}_2|^3} - 1 \right) > \omega^2 \left( (1 - q) + q - 1 \right) = 0$ .

Cette quantité est donc strictement positive, ce qui assure que le déterminant est négatif, c'est-à-dire que L₁ est un point selle, ce qui en fait un point instable.

Le cas de L₂ et L₃

On pose, pour simplifier les notations,

$u_1 = \frac{|{\boldsymbol{r}} - {\boldsymbol{r}}_1|}{R}$ ,

$u_2 = \frac{|{\boldsymbol{r}} - {\boldsymbol{r}}_2|}{R}$ .

On s'intéreresse donc au signe de la quantité

$\frac{\Omega_{,yy}}{\omega^2} = (1 - q) \frac{1}{u_1^3} + q \frac{1}{u_2^3} - 1$ ,

soit

$\frac{\Omega_{,yy}}{\omega^2} = (1 - q) \left(\frac{1}{u_1^3} - 1\right) + q \left(\frac{1}{u_2^3} - 1\right)$ ,

sachant que u₁ et u₂ sont reliés entre eux par le fait que leur différence est égale à 1 et qu'ils définissent un point de Lagrange, soit la relation

$- (1 - q) \frac{1}{u_1^2} - q \frac{1}{u_2^2} + \frac{|{\boldsymbol{r}}|}{R} = 0$ .

La distance du point de Lagrange au centre de gravité du système peut s'écrire, pour le point L₂,

$\frac{|{\boldsymbol{r}}|}{R} = u_2 + (1 - q) = u_1 - q$ ,

relations que l'on peut combiner en

$\frac{|{\boldsymbol{r}}|}{R} = q u_2 + (1 - q) u_1$ .

La position du point L₂ est donc donnée par

$- (1 - q) \left(\frac{1}{u_1^2} - u_1\right) - q \left(\frac{1}{u_2^2} - u_2 \right) = 0$ .

On pose alors

$A = (1 - q) \left(\frac{1}{u_1^2} - u_1\right)$ ,

$B = q \left(\frac{1}{u_2^2} - u_2 \right)$ .

On a donc, d'une part

A + B = 0

et, d'autre part

$\frac{\Omega_{,yy}}{\omega^2} = \frac{A}{u_1} + \frac{B}{u_2}$ .

Autrement dit,

$\frac{\Omega_{,yy}}{\omega^2} = \frac{A}{u_1} + \frac{B}{u_1} + B \left(\frac{1}{u_2} - \frac{1}{u_1} \right)$ .

Le premier terme du membre de droite est nul en vertu de la relation A + B = 0. Il reste donc

$\frac{\Omega_{,yy}}{\omega^2} = B \left(\frac{1}{u_2} - \frac{1}{u_1} \right)$ .

Or, pour le point L₂, on est situé plus près du corps 2 que du corps 1. Par conséquent, u₂ est plus petit que u₁, et, par conséquent, $\left(\frac{1}{u_2} - \frac{1}{u_1} \right)$ est positif. Le signe de la dérivée seconde correspond donc à celui de B, qui lui-même est déterminé par la valeur de u₂ : si cette quantité est supérieure à 1, alors B est négatif, alors que, dans le cas contraire, B est positif, ce qui implique que le point est instable. Le point de Lagrange L₂ est situé au-delà du corps 2. La force totale (gravitationnelle plus centrifuge) s'exerçant dans cette région est d'abord tournée vers le corps 2 quand on est proche de celui-ci, puis s'annule en L₂ et est après dirigée à l'opposé de L₂. Au point tel que u₂ vaut 1, la composante de cette force, selon l'axe reliant les deux corps, est donnée (à une constante multiplicative positivite près) par

$\delta f \propto - (1 - q) \left(\frac{1}{u_1^2} - u_1\right) - q \left(\frac{1}{u_2^2} - u_2 \right)$ ,

avec, ici,

u 2 = 1

u 1 = 2

soit

$\left.\delta f \right|_{u_2 = 1} \propto \frac{7}{4} (1 - q)$ .

Cette quantité étant strictement positive, le point u₂ = 1 de l'axe est situé au-delà du point L₂. En conséquence, au point L₂, u₂ est inférieur à 1, donc B est positif, donc le point est bien un point selle, ce qui assure de son instabilité. Une démonstration strictement analogue peut être faite pour le point L₃, ce qui achève la démonstration de leur instabilité du fait de leur caractère de point selle.

Temps caractéristiques en L₁ et L₂ pour les systèmes à grande hétérogénéité de masse

Une des applications les plus importantes de l'instabilité des points de Lagrange, L₁ et L₂, réside dans le fait que des satellites artificiels peuvent être envoyés en ces points du système Terre-Soleil (voir ci-dessous). Pour de tels satellites, des corrections de trajectoires régulières doivent être appliquées afin de conserver le satellite au voisinage du point. Ce temps caractéristique peut être évalué dans le cas où le rapport de masse des deux corps du système est élevé. Dans ce cas, le temps caractéristique γ^-1 d'instabilité est donné par

$\frac{1}{\gamma} = \frac{T}{2 \pi \sqrt{1 + 2\sqrt{7}}}$ ,

où T est la période orbitale du système. Dans le cas du système Terre-Soleil, où T est légèrement supérieur à 365 jours, le temps caractéristique d'instabilité est alors de 23 jours et 4 heures.

Par ailleurs, la composante stable de la trajectoire se fait à la pulsation

$\omega_{\rm Traj} = \omega \sqrt{2 \sqrt{7} - 1}$ ,

soit, de façon équivalente, avec la période

$T_{\rm Traj} = \frac{T}{\sqrt{2 \sqrt{7} - 1}}$ ,

ce qui, dans le même cas de figure que ci-dessus, donne une période de 176 jours.

Démonstration

L'équation donnant les valeurs propres du système est toujours

$\lambda^4 + (\Omega_{,xx} + \Omega_{,yy} + 4 \omega^2) \lambda^2 + \Omega_{,xx} \Omega_{,yy} - \Omega_{,xy}^2 = 0$ ,

avec, pour les points L₁ et L₂,

Ω, x y = 0

Ω, x x = - 2Ω, y y - 3ω 2

$\Omega_{,yy} = \omega^2 \left(\frac{q}{u_2^3} + \frac{1 - q}{u_1^3} - 1\right)$ .

en se restreignant aux termes d'ordre le plus bas en q, u₁ vaut 1, et u₂ est déterminé par la relation donnée par le preimier tableau de cette page. On a ainsi

$\Omega_{,yy} = \omega^2 \left(3 + 1 - 1\right) = 3 \omega^2$ ,

Ω, x x = - 9ω 2

L'équation polynomiale devient alors

λ 4 - 2ω 2 λ 2 - 27ω 4 = 0

dont les solutions sont

$\lambda_\pm^2 = \left(1 \pm 2 \sqrt{7}\right) \omega^2$ .

La solution positive à cette équation indique que les écart au point d'équilibre croissent exponentiellement au cours du temps selon la relation

$\delta X \propto e^{\lambda_+ t}$ ,

avec

$\lambda_+ = \omega \sqrt{1 + 2 \sqrt{7}} = \frac{2 \pi}{T} \sqrt{1 + 2 \sqrt{7}}$ .

Le temps caractéristique associé est donc

$\frac{1}{\gamma} = \frac{1}{\lambda_+} = \frac{T }{2 \pi \sqrt{1 + 2 \sqrt{7}}}$ .

De la même façon, il existe des trajectoires périodiques dont la pulsation est donnée par les racines complexes de l'équation, soit

$\omega_{\rm Traj} = \omega \sqrt{2 \sqrt{7} - 1}$ ,

c'est-à-dire une période de

$T_{\rm Traj} = \frac{T}{\sqrt{2 \sqrt{7} - 1}}$ .

La structure des orbites en présence d'instabilité

Une fois les valeurs propres d'un point instable connues, une trajectoire au voisinage d'un point de Lagrange va être une combinaison linéaire des vecteurs propres associés aux valeurs propres. En notant λ'_i l'une de ces valeurs propres, le vecteur propre associé a pour composantes

$V_i = \left(\begin{array}{c} 1 \\ \mu_i \\ \lambda_i \\ \lambda_i \mu_i \end{array}\right)$ ,

avec

$\mu_i = - \frac{1}{2 \omega} \left(\frac{\Omega_{,xx}}{\lambda} + \lambda \right)$ ,

et une trajectoire est de la forme

$\left(\begin{array}{c} \delta X \\ \delta Y \\ \dot{\delta X} \\ \dot{\delta Y} \end{array}\right) = \sum_i \alpha_i V_i$ ,

où les quantités α_i sont des nombres quelconques déterminés par la valeur des δX, δY et de leur dérivée à un instant donné. Dans le cas des trois points de Lagrange instables, le déterminant de la matrice des dérivées seconde est négatif, ce qui implique que le discriminant de l'équation du second degré en λ² possède des racines réelles de signe opposé, et, qu'au final, les valeurs propres recherchées sont deux nombres imaginaires purs opposés et deux nombres réels opposés. Une trajectoire générique comprend donc, dans le plan orbital, une composante périodique (liée aux racines imaginaires pures), une composante amortie (liée à la racine réelle positive), et une composante instable. Pour une position δX, δY donnée, il est toujours possible de choisir une vitesse telle que les deux vecteurs propres aux racines réelles ne contribuent pas à la solution correspondante. La trajectoire obtenue est alors périodique, la période étant donnée par la racine complexe. Une telle solution n'est cependant pas stable. Un écart de trajectoire infime va en réalité rajouter à la trajectoire une composante instable, qui va, peu à peu, éloigner la trajectoire de sa composante périodique. On dit que la trajectoire obtenue n'est pas dynamiquement stable. Ceci est une généralisation du fait qu'un objet situé exactement sur un point de Lagrange instable est dans une situation instable : un petit écart à cette position d'équilibre, inéluctablement généré par les perturbations causées par les autres corps du système, finira par éloigner l'objet de sa position initiale. La même chose se produit pour des trajectoires situées autour du point d'équilibre instable.

Pertinence du concept

Le calcul ci-dessus fait référence à une configuration où les deux corps du système sont en orbite circulaire. Néanmoins, le concept de point de Lagrange prévaut pour tout type d'orbite, y compris elliptique. On peut donc définir ces points dans tout système à deux corps liés gravitationnellement. Par contre, les trajectoires, stables ou instables, autour des différents points de Lagrange dépendent explicitement de la circularité ou non de l'orbite des deux corps du système.

Applications

Étant donné les questions de stabilité évoquées plus haut, on ne trouve pas d'objet naturel autour des points L₁, L₂ et L₃ dans le système solaire. Cependant, ils représentent tout de même un intérêt pour les réalisations scientifiques, car ils permettent des économies de combustible pour le contrôle d'orbite et d'attitude. Ceci n'est pas valable pour le point L₃, du fait de son éloignement de la Terre dont la seule application était que les auteurs de science-fiction et de bande dessinée aimaient y placer une Anti-Terre d'autant plus utopique que la masse de cette planète-jumelle y était bien trop élevée par rapport à la théorie énoncée plus haut. En revanche, des missions spatiales utilisent L₁ et L₂ : c'est le cas de la sonde SoHO (Solar and Heliospheric Observatory) une station d'observation du Soleil située, depuis 1995, au point L₁ à 1,5 million de km de la Terre. Le point L₂ est particulièrement bien adapté pour observer le cosmos. Depuis 2001, son voisinage est occupé par le satellite WMAP, chargé d'étudier le fond diffus cosmologique. Il a été rejoint par les satellites Planck et Herschel en 2009, et le sera par Gaia en 2011 et par le James Webb Space Telescope en 2013.

L₄ et L₅ étant stables, on y trouve de nombreux corps naturels. Dans le système Jupiter-Soleil, plusieurs centaines d'astéroïdes, appelés astéroïdes Troyens, s'y agglutinent (près de 1800 en avril 2005). On en compte quelques-uns dans les systèmes Neptune-Soleil et Mars-Soleil. Curieusement, il semblerait que le système Saturne-Soleil ne soit pas en mesure d'en accumuler, à cause des perturbations joviennes. On trouve également des objets à ces points dans le système Saturne-satellites de Saturne : Saturne-Téthys avec Télesto et Calypso aux points L₄ et L₅, et Saturne-Dioné avec Hélène au point L₄ et Pollux au point L₅. Dans le système Terre-Soleil, il n'y a pas d'objet connu de grande taille aux points Troyens, mais on y a découvert une légère surabondance de poussières en 1950. De légers nuages de poussières sont également présents pour le système Terre-Lune ; cela a fait renoncer à y placer un télescope spatial comme il avait été envisagé.

En science fiction, de par leur stabilité, les points L₄ et L₅ du système Terre-Lune abritent souvent de gigantesques colonies spatiales.

Référence

(en) Richard Fitzpatrick, Analytical Classical Dynamics, cours donné à l'Université du Texas à Austin Voir en ligne (versions pdf et html).

Voir aussi

Liste d'objets situés à un point de Lagrange

Note

↑ Si l'on définit q comme le rapport de la plus petite masse à la totale, seules les valeurs de q inférieures à 0,5 ont un sens, puisque les valeurs plus grandes correspondent au rapport de la plus grande masse à la masse totale.

Portail de l’astronomie
Portail de l’astronautique

Ce document provient de « Point de Lagrange ».

Catégories : Mécanique céleste | Joseph-Louis Lagrange

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Les points de Lagrange de Wikipédia en français (auteurs)

Игры ⚽ Поможем решить контрольную работу

Regardez d'autres dictionnaires:

Points de Lagrange — Point de Lagrange Pour les articles homonymes, voir Lagrange. Un point de Lagrange (noté L1 à L5), ou, plus rarement, point de libration, est une position de l espace où les champs de gravité de deux corps en orbite l un autour de l autre, et de… … Wikipédia en Français
Les Villes de l'espace — Les Villes de l espace : vers le peuplement, l industrialisation et la production d énergie dans l espace Auteur Gerard K. O Neill Genre Colonisation de l espace Version originale Titre original The High Frontier: Human Colonies in Space… … Wikipédia en Français
Les phases de la lune — Phase lunaire En astronomie, une phase lunaire désigne une portion de Lune illuminée par le Soleil et vue à partir de la Terre. La Lune tournant en orbite autour de la Terre, les positions relatives du Soleil, de la Terre et de la Lune changent… … Wikipédia en Français
LES MONTS DE VERDUN 1917 — Bataille de Verdun Pour les articles homonymes, voir Verdun … Wikipédia en Français
Point de Lagrange — Pour les articles homonymes, voir Lagrange. Un point de Lagrange (noté L1 à L5), ou, plus rarement, point de libration, est une position de l espace où les champs de gravité de deux corps en orbite l un autour de l autre, et de masses… … Wikipédia en Français
Point de lagrange — Pour les articles homonymes, voir Lagrange. Un point de Lagrange (noté L1 à L5), ou, plus rarement, point de libration, est une position de l espace où les champs de gravité de deux corps en orbite l un autour de l autre, et de masses… … Wikipédia en Français
Lagrange inversion theorem — In mathematical analysis, the Lagrange inversion theorem, also known as the Lagrange Bürmann formula, gives the Taylor series expansion of the inverse function of an analytic function. Theorem statementSuppose the dependence between the variables … Wikipedia
Dualité dans les programmes d'optimisation — Multiplicateur de Lagrange Pour les articles homonymes, voir Théorème de Lagrange. La méthode des multiplicateurs de Lagrange permet de trouver un optimum, sur la figure le point le … Wikipédia en Français
Multiplicateur De Lagrange — Pour les articles homonymes, voir Théorème de Lagrange. La méthode des multiplicateurs de Lagrange permet de trouver un optimum, sur la figure le point le … Wikipédia en Français
Multiplicateur de Lagrange — Pour les articles homonymes, voir Théorème de Lagrange. La méthode des multiplicateurs de Lagrange permet de trouver un optimum, sur la figure le point le plus élevé possible, tout en satisfaisant une contrainte, sur la figure un … Wikipédia en Français

Dictionnaires et Encyclopédies sur 'Academic'

Les points de Lagrange

Point de Lagrange

Sommaire

Historique

Définition

Calcul de la position des points de Lagrange

Préliminaires

Équation fondamentale

Les deux cas à considérer

Cas des points L₄ et L₅

Cas des points L₁ à L₃

Forme réduite et solution dans le cas où le rapport entre les masses est faible

Le point L₁

Le point L₂

Le point L₃

Stabilité

Préliminaire

Cas des points de Lagrange L₄ et L₅

Cas des points de Lagrange L₁ à L₃

Le cas de L₁

Le cas de L₂ et L₃

Temps caractéristiques en L₁ et L₂ pour les systèmes à grande hétérogénéité de masse

La structure des orbites en présence d'instabilité

Pertinence du concept

Applications

Référence

Voir aussi

Note

Regardez d'autres dictionnaires:

Share the article and excerpts

Dictionnaires et Encyclopédies sur 'Academic'

Wikipédia en Français

Les points de Lagrange

Point de Lagrange

Sommaire

Historique

Définition

Calcul de la position des points de Lagrange

Préliminaires

Équation fondamentale

Les deux cas à considérer

Cas des points L4 et L5

Cas des points L1 à L3

Forme réduite et solution dans le cas où le rapport entre les masses est faible

Le point L1

Le point L2

Le point L3

Stabilité

Préliminaire

Cas des points de Lagrange L4 et L5

Cas des points de Lagrange L1 à L3

Le cas de L1

Le cas de L2 et L3

Temps caractéristiques en L1 et L2 pour les systèmes à grande hétérogénéité de masse

La structure des orbites en présence d'instabilité

Pertinence du concept

Applications

Référence

Voir aussi

Note

Regardez d'autres dictionnaires:

Share the article and excerpts

Direct link

Cas des points L₄ et L₅

Cas des points L₁ à L₃

Le point L₁

Le point L₂

Le point L₃

Cas des points de Lagrange L₄ et L₅

Cas des points de Lagrange L₁ à L₃

Le cas de L₁

Le cas de L₂ et L₃

Temps caractéristiques en L₁ et L₂ pour les systèmes à grande hétérogénéité de masse