Leçons sur les hypothèses cosmogoniques (Poincaré, 1911)/Chapitre 07

Libraire scientifique A. Hermann et fils (p. 131-189).

CHAPITRE vii.

THÉORIE DE Sir G. H. DARWIN.


I. — Généralités.

94.Dans l’histoire, tant passée que future, des astres du système solaire, Sir G. H. Darwin attribue un rôle essentiel à l’influence des marées[1].

Considérons les marées produites sur une planète T par un astre troublant L (fig. 25) : la planète sera, par exemple, la Terre et l’astre Figure 25
fig.25.
troublant la Lune. Rendons-nous compte de l’action du frottement de la marée sur la rotation de la planète et de la réaction qui en résulte sur le mouvement de l’astre producteur de la marée. Si la théorie statique des marées pouvait être appliquée, la surface des océans terrestres serait un ellipsoïde allongé vers la Lune L. Mais les frottements dus à la viscosité ont pour effet de produire un décalage, un retard de la marée sur l’instant du passage de la Lune au méridien. Le grand axe de l’ellipsoïde terrestre fera donc un certain angle avec la ligne TL qui joint les centres de gravité des deux astres[2]. Par conséquent, la résultante des attractions de la Lune sur les molécules terrestres ne passe pas par le centre de gravité T de la Terre, mais est appliquée à un certain point A : cette résultante a donc un moment par rapport au point T et tend, par suite, à ralentir le mouvement de rotation de la Terre.

De même, l’attraction de la Terre sur la Lune n’est pas dirigée suivant le rayon vecteur LT ; ce n’est plus rigoureusement une force centrale : elle a une petite composante tangentielle qui va troubler le mouvement orbital de la Lune. Cette composante tangentielle, dirigée dans le sens même du mouvement de la Lune sur son orbite produira évidemment un effet inverse de celui d’une résistance de milieu, cette dernière étant une force tangentielle dirigée à l’opposé du mouvement orbital de la Lune. Il y aura donc augmentation du grand axe de l’orbite lunaire et, conséquemment, diminution de la vitesse angulaire de révolution.

Cette augmentation du grand axe de l’orbite lunaire pouvait encore se prévoir autrement. L’effet de freinage des marées est de diminuer la vitesse angulaire de la Terre et par suite son moment de rotation. D’autre part, le moment de rotation total du système Terre-Lune doit rester constant. Donc le moment de rotation dû au mouvement orbital de la Lune doit augmenter, c’est-à-dire que la distance de la Lune à la Terre doit s’accroître.

Nous prévoyons donc dès maintenant et sans calculs que les deux principaux effets du frottement des marées sont la diminution de la rotation terrestre et l’augmentation corrélative de la distance de la Lune. En d’autres termes, la durée du jour augmentera, ainsi que celle du mois.

95.Pour expliquer l’accélération séculaire du moyen mouvement de la Lune dont la gravitation, comme on le sait, ne paraît pas rendre entièrement compte, Delaunay avait déjà proposé d’admettre une augmentation de la durée du jour sidéral, due précisément à l’action du frottement des marées sur le fond des océans. Dans cette hypothèse, l’accélération séculaire de la Lune ne serait pas réelle ; ce ne serait qu’une apparence due à ce que notre unité de temps, le jour sidéral, augmenterait peu à peu. Bien plus, l’effet des marées, loin de produire réellement une accélération du moyen mouvement de la Lune, augmenterait au contraire la distance de cet astre à la Terre et sa période de révolution, comme nous venons de le remarquer. N’est-il pas à craindre alors que ce retard effectif ne l’emporte sur l’accélération apparente ? Pour répondre à cette objection, Delaunay soutenait que le coefficient du retard imprimé à la Lune est moindre que celui du retard de la rotation de la Terre, ce qui se trouve être exact, semble-t-il. Mais, si l’on ne tient compte que du frottement des océans actuels sur leur fond, l’effet total est tout à fait minime et insignifiant. Il faudrait supposer que l’intérieur du globe terrestre n’est pas rigide, qu’il est plus ou moins visqueux, et que, par conséquent, les marées internes y produisent des frottements intenses.

II. — L’excentricité et l’inclinaison de l’orbite lunaire sont supposées nulles.

96.Quoi qu’il en soit, et pour en revenir au point de vue cosmogonique, on doit supposer, avec Sir G. H. Darwin, que dans ses états antérieurs la Terre était fluide et visqueuse. Elle a alors subi des marées dans toute sa masse (bodily tides), et les frottements dus à ces marées internes étaient incomparablement plus énergiques que ceux qu’on peut attribuer actuellement aux marées océaniques.

97.Nous nous proposons d’étudier de plus près les variations que subissent, du fait de ces marées, la distance de la Lune et la rotation terrestre.

Pour simplifier, nous supposerons tout d’abord que le plan de l’orbite lunaire coïncide avec celui de l’équateur terrestre et que cette orbite est circulaire. Il est bien clair que, par raison de symétrie, l’orbite restera alors indéfiniment dans le plan équatorial. Mais restera-t-elle toujours circulaire ? Montrons qu’il en sera bien ainsi et que, si l’excentricité est initialement nulle, elle le restera toujours. Soient la force perturbatrice perpendiculaire au rayon vecteur[3] et le chemin parcouru par la Lune pendant le temps  : le travail de la force perturbatrice est

désigne l’angle de avec . Soit le moment de rotation ; sa dérivée est égale au moment de la force perturbatrice ; on a donc l’égalité

désigne le rayon vecteur et l’angle de avec la perpendiculaire à ce rayon vecteur.

Dans le cas d’une orbite circulaire, et sont tous deux égaux à 1, et l’on a

étant le moyen mouvement. Nous tirons alors des formules précédentes

(1)

Or, est la différentielle de la constante des forces vives

et le moment de rotation a pour valeur

désigne la masse de l’ensemble Terre-Lune, et sont le grand axe et l’excentricité de l’orbite lunaire. L’équation (1) s’écrit donc

L’orbite étant supposée circulaire, nous faisons  ; il vient

(2)

Or, on a

les termes en se détruisent donc dans l’équation (2) qui donne alors

L’excentricité ne subit donc pas de variation : nulle au début, elle restera nulle.

On pourrait faire ici une objection. Il n’est pas étonnant, dira-t-on, que soit nul, et d’ailleurs cela ne prouve rien : en effet

et, pour une orbite circulaire, est nul. Mais il est facile de reconnaître que, dans le calcul ci-dessus de (ou pour être plus précis, de sa partie séculaire), nous n’avons négligé que des termes en . Donc est de l’ordre de et est de l’ordre de  : il est par suite nul pour une orbite circulaire, et l’objection n’a pas de portée.

98.Considérons donc la Lune décrivant autour de la Terre une orbite circulaire dans le plan de l’équateur. Désignons par la vitesse angulaire de rotation de la Terre autour de son axe, par la vitesse de révolution de la Lune autour de la Terre, et posons

Nous allons écrire que le moment de rotation du système est constant (principe des aires), et que l’énergie mécanique diminue par suite du frottement (principe de dégradation de l’énergie).

Soit le moment d’inertie de la Terre autour de son axe. Le moment de rotation de la Terre est et sa demi-force vive est .

Tenons compte à présent du mouvement de révolution de la Lune et de la Terre autour de leur centre de gravité commun. Le moment de rotation dû à ce mouvement est proportionnel à , car il a pour valeur

il est donc proportionnel à qui lui-même est proportionnel à d’après l’équation

et la constante des forces vives est proportionnelle à c’est-à-dire à .

Nous choisirons les unités de façon à simplifier les coefficients de proportionnalité. Nous prendrons tout d’abord les unités de masse et de longueur de manière que Ensuite nous choisirons l’unité de temps de telle sorte que, pour le moment de rotation du système Terre-Lune dans sa révolution autour de son centre de gravité soit égal à , c’est-à-dire à 1. Alors, le moment de rotation dû au mouvement orbital est non seulement proportionnel, mais égal à . D’ailleurs, l’équation du viriel (no 74, p. 91)

qui, puisque le mouvement est circulaire, s’écrit ici

nous apprend que l’énergie totale est égale à c’est-à-dire à la demi-force vive changée de signe. On en conclut immédiatement qu’avec les unités choisies, la constante des forces vives, proportionnelle à est égale à

Pour l’ensemble du système Terre-Lune, nous avons donc le moment de rotation total

et l’énergie totale

Le moment de rotation reste toujours constant : nous écrivons donc

(3)

Quant à l’énergie, elle va constamment en diminuant, absorbée qu’elle est par le frottement qui la transforme en chaleur : si donc nous posons

ira toujours en décroissant.

Si nous remplaçons par sa valeur tirée de l’équation (3), devient fonction de . Nous obtiendrons ses maxima et ses minima en annulant  :

Annuler revient donc à écrire

ou encore, se rappelant la définition de ,

Ainsi, lorsque l’énergie est maximum ou minimum, la vitesse angulaire de rotation de la Terre est égale à la vitesse angulaire de révolution de la Lune.

Figure 26
fig.26.

99.Prenons, avec Sir G. H. Darwin, pour abscisse et pour ordonnée ; et traçons (fig. 26) la droite

(3)

et la courbe

(4)

Cette courbe, Sir G. H. Darwin la nomme courbe de rigidité parce que si le point représentatif est sur cette courbe, la durée de rotation de la Terre égale la durée de révolution de la Lune, et l’ensemble Terre-Lune tourne d’un seul bloc à la façon d’un corps solide.

Deux cas sont à distinguer suivant la valeur de la constante . Ou bien la droite (3) coupe la courbe de rigidité (4) en deux points C et D ; ou bien la droite (3) ne coupe pas la courbe (4).

Prenons d’abord le premier cas, et étudions en fonction de Nous avons

Si nous prenons pour abscisse et pour ordonnée, cette équation représente une courbe telle que celle de la figure 27 : les points C′ et Figure 27
fig.27.
D’ où passe par un maximum et par un minimum correspondent aux mêmes abscisses que les points C et D de la figure 26 situés sur la courbe de rigidité.

Dans le second cas, où la droite (3) ne coupe pas la courbe de rigidité (4), la fonction de ne présente plus ni maximum ni minimum, et la courbe de la figure 27 est remplacée par celle de la figure 28.

Nous pouvons maintenant suivre les changements subis par le système formé par la planète et son satellite.

Supposons que l’état initial soit représenté par un point de la droite ABCDE (fig. 26) : alors le point représentatif restera toujours sur cette droite, mais de telle façon que aille toujours en décroissant.

Si le point représentatif initial est situé entre A et B (fig. 26), le point (fig. 27) est situé sur la branche A′B′ et, puisque diminue toujours, il décrit cette branche en allant de A′ vers B′ . Le point (fig. 26) décrit donc la portion de droite AB en avançant toujours à droite, jusqu’au point B qui représente l’état final. En ce point on a c’est-à-dire autrement dit, le satellite finira par tomber sur la planète. Remarquons que, dans ce cas, le satellite circule toujours autour de sa planète dans un sens inverse de celui de la rotation de celle-ci Figure 28
fig.28.

Nous arriverions à une conclusion analogue — le satellite tomberait finalement sur la planète — si le point représentatif initial était situé entre B et C, car alors le point (fig. 27) décrirait la branche de courbe C′B′′ en allant toujours vers B′′. Cette fois, les deux mouvements de rotation de la planète et de circulation du satellite sont de même sens  ; mais puisqu’on a

c’est-à-dire

le mouvement angulaire du satellite est plus rapide que celui de la planète : le mois est plus court que le jour.

Supposons maintenant que le point représentatif soit entre C et D (fig. 26). Les deux mouvements sont alors directs et le mois est plus long que le jour  : c’est le cas offert par le système Terre-Lune. Alors le point (fig. 27) décrira l’arc C′D′ en marchant vers D′ puisque décroît. L’état final est représenté par le point D de la courbe de rigidité (fig. 26), pour lequel les deux durées de rotation et de révolution sont égales.

Enfin, si le point représentatif était situé entre D et E, l’état final serait encore le point D, puisque le point (fig. 27) décrirait alors la branche E′D′. Dans ce cas, le mois est toujours plus court que le jour le mouvement du satellite est toujours direct mais celui de la planète a pu commencer par être rétrograde

Dans le cas où la droite

ne couperait pas la courbe de rigidité, la figure 27 serait remplacée par la figure 28. L’état final serait toujours le point B′ ou B′′ c’est-à-dire que le satellite finirait toujours par tomber sur la planète.

La plupart des satellites connus correspondent au cas où les deux mouvements sont de même sens, mais avec le mois plus long que le jour : le point représentatif est alors (fig. 26) entre C et D.

100.Étudions spécialement le système Terre-Lune. Pour ce système, on a les valeurs numériques suivantes

et le point représentatif se trouve dans une position telle que P. Ce point se déplace lentement vers le point D qui représente l’état final. Si, au lieu de chercher à prévoir l’avenir, nous remontons dans le passé, nous pouvons dire que le système Terre-Lune est parti de l’état initial représenté par C. Dans cet état initial, le mois égalait le jour et leur durée commune était de 5h,36. La durée du mois et celle du jour se sont mises ensuite à croître, la première plus vite que la seconde, en même temps que la Lune s’éloignait de la Terre. Lorsque l’état final D sera atteint, le mois sera redevenu égal au jour, leur durée commune étant d’environ 55 jours actuels. Cet état final serait définitif si la Terre et la Lune existaient seules. Mais le Soleil produit aussi sur la terre des marées qui continueront à retarder sa rotation : le jour et le mois continueront donc à croître, quoique lentement, et la Lune s’éloignera de plus en plus de la Terre, qui pourra finir par perdre son satellite.

Dans l’état initial, la durée de révolution de la Lune, 5h,36, était courte et par suite la distance de cet astre à la Terre était faible : le calcul montre que cette distance était de 21/2 rayons terrestres. Donc, d’après les idées de Sir G. H. Darwin, la Lune aurait pris naissance tout près de la Terre, et son orbite se serait peu à peu élargie et dilatée. Nous remarquons là un contraste complet avec la théorie de Faye selon laquelle, nous l’avons vu, le satellite se formerait très loin de l’astre central, son orbite se rétrécissant peu à peu à mesure que la masse de l’astre central augmenterait par suite de la condensation.

101.Dans tout ce qui précède, nous avons regardé , moment d’inertie de la Terre, comme constant. Or, la Terre, en se contractant par suite du refroidissement, diminue de volume et décroît. Nous étudierons un peu plus loin l’effet qui en résulte. Cependant disons dès maintenant que, pour la Terre, cet effet paraît peu important et ne change pas beaucoup l’allure générale des phénomènes.

102.Pour le système Terre-Lune, le point représentatif actuel P (fig. 26) est assez voisin du point D, puisque est sensiblement plus grand que Par conséquent, ce système est relativement assez près de son état final.

Si nous nous occupons maintenant des systèmes formés par les autres planètes et leurs satellites, nous constatons, au contraire, que le rapport se trouve être extrêmement petit. Par conséquent, leur point représentatif, au lieu d’être voisin du point D, est voisin du point C. Ces systèmes sont donc beaucoup plus près de leur état initial que de leur état final. Cela tient à ce que, pour ces systèmes, la masse des satellites est extrêmement faible par rapport à la masse de la planète ; les marées sont donc peu importantes et n’ont pas encore eu le temps de produire un effet considérable.

Pour le système formé par le Soleil et l’ensemble des planètes, la masse des planètes est sans doute fort petite par rapport à celle du Soleil, mais les rayons de gyration sont très grands, et l’on constate que le rapport est aussi très grand. Il semble donc que ce système doive être très près de son état final. Mais cette affirmation ne supporte pas l’examen. Les planètes, en effet, produisent sur le Soleil des marées qui retardent sa rotation ; et c’est de ces marées qu’il faut tenir compte si l’on cherche les variations de la rotation du Soleil. Mais, inversement, le Soleil produit des marées sur les planètes, et ces dernières marées sont beaucoup plus fortes que les premières : leur effet devrait donc s’être produit depuis longtemps, c’est-à-dire que les planètes devraient toujours tourner une même face vers le Soleil, ce qui n’a pas lieu.

III. — Cas général.

103.Dans l’exposé qui précède, nous avons, pour simplifier, supposé nulles l’excentricité de l’orbite lunaire et l’inclinaison de cet orbite sur le plan de l’équateur terrestre. Nous allons à présent nous affranchir de ces hypothèses et entrer dans le détail de la théorie de Sir G. H. Darwin.

Auparavant, il est nécessaire de rappeler quelques points de la théorie statique des marées[4].

Dans cette théorie statique, on suppose que la mer prend à chaque instant sa forme d’équilibre : la surface libre de l’océan est donc une surface de niveau relativement à la somme des potentiels de la gravité et de l’astre perturbateur qui produit la marée.

Le potentiel de la gravité peut se représenter par

(5)

désignant la valeur constante de ce potentiel au niveau moyen des mers et la dénivellation de la particule liquide superficielle envisagée[5].

Le potentiel dû à l’astre perturbateur, nous le désignons par . Soient la masse de l’astre perturbateur L (qui sera, par exemple, la Lune), sa distance à un point A de la surface de la mer (fig. 29), Nous avons

(6)

Comme nous rapportons les points de la mer à des axes de coordonnées invariablement liés à la Terre, il faut, pour pouvoir regarder ces axes comme fixes, appliquer à chaque point A les forces apparentes dues à leur mouvement. Mais, puisqu’il ne s’agit ici que de l’équilibre, la force centrifuge composée n’intervient pas ; il ne reste que la force d’inertie dans le mouvement d’entraînement du point A Figure 29
fig.29.
avec les axes. La force provenant de la rotation diurne a déjà été prise en considération dans Il suffit donc, aux forces réelles, d’ajouter la force d’inertie due à la translation des axes, c’est-à-dire, puisque l’origine est au centre de la Terre, une force accélératrice égale et contraire à l’accélération que l’astre perturbateur tend à imprimer à ce point.

Soient les composantes rectangulaires de . À chaque point A, on devra appliquer une force de composantes

comme ces composantes ne dépendent que du temps, et non des coordonnées du point A, elles peuvent être considérées comme les dérivées partielles de la fonction

(7)

en appelant le rayon moyen des mers, égal sensiblement à la distance du point A au centre T de la Terre, et en désignant par l’angle de avec le rayon TA. Finalement, en écrivant que la somme des trois expressions (5), (6) et (7) est égale à une constante, nous obtiendrons l’équation de la surface libre des océans rapportée à des axes invariablement liés à la Terre :

(8)

Telle est l’équation qui, dans la théorie statique, donne la dénivellation Nous allons la transformer.

104.L’accélération que la Lune L imprime au centre de la Terre T a pour valeur et est dirigée suivant la droite TL Donc l’angle n’est autre que l’angle désigné par sur la figure 29 (c’est la distance zénithale géocentrique de la Lune), et l’on a

Comme on a

le potentiel dû à la Lune peut s’écrire

en négligeant les termes de l’ordre de . Nous écrirons

en posant

Par suite, nous pouvons écrire l’équation (8) ainsi

et étant des constantes ne dépendant pas du point A, cette équation nous donne finalement, pour la dénivellation

La constante du second membre se déterminerait en écrivant que le volume total de la mer est constant, c’est-à-dire que l’intégrale

étendue à tous les éléments de la surface des mers, est nulle. On reconnaît ainsi que, si la mer recouvre toute la surface de la Terre, cette constante est nulle, et l’on a alors la formule définitive de la dénivellation statique

(9)

105.Dans la théorie statique, on ne tient pas compte de l’inertie des eaux de la mer. Dans quelle mesure cela est-il légitime ? L’inertie, proportionnelle à l’accélération, joue un rôle dans les marées à courte période ; au contraire, pour les marées à longue période, l’effet de l’accélération, c’est-à-dire de l’inertie, est négligeable. Mais que doit-on entendre par longue et courte période ? La période d’une marée pourra être qualifiée de longue ou de courte suivant que cette période sera ou ne sera pas très considérable par rapport à la période d’oscillation propre de la mer.

Par conséquent, si la période d’oscillation propre est très courte, toutes les marées, même les marées semi-diurnes, pourront être considérées comme étant à longue période, et la théorie statique s’appliquera. Or, si, avec Sir G. H. Darwin, nous considérons la Terre aux époques reculées où elle était entièrement liquide, nous trouvons que sa période propre d’oscillation est effectivement très petite par rapport à la demi-journée : il est donc, légitime de se contenter de la théorie statique.

Au contraire, l’inertie jouerait un rôle important, et la théorie statique deviendrait tout à fait insuffisante, si l’on voulait étudier le mouvement des océans actuels, dont la période propre d’oscillation est de l’ordre du demi-jour.

106.Mais si, dans le calcul de l’action des marées produites par la Lune sur la Terre, nous pouvons négliger l’inertie, nous devons en revanche tenir compte de la viscosité, puisque ce sont précisément ses effets que nous voulons étudier.

Reprenons la formule de la dénivellation statique

(9)

Le mouvement de l’astre perturbateur L étant connu en fonction du temps , le troisième membre de cette formule peut être développé sous forme trigonométrique

les étant des fonctions des coordonnées du lieu géographique et les étant des constantes. Nous envisagerons séparément chacun des termes de la somme et, pour tenir compte de la viscosité, nous ajouterons au second membre de l’équation qui donne un terme négatif proportionnel L’équation en prend ainsi la forme

ou encore

(10)

étant une constante proportionnelle à la viscosité du liquide.

Pour intégrer cette équation linéaire à second membre, nous considérerons l’équation suivante :

(11)

désigne l’imaginaire et nous ne conserverons que la partie réelle de sa solution. Cette solution est de la forme

en la substituant dans l’équation (11), il vient pour déterminer l’équation

Posant

(12)

on aura

par suite

Et l’intégrale de l’équation (11) est

dont la partie réelle

est l’intégrale cherchée de l’équation (10).

Nous voyons donc que l’effet de la viscosité est d’abord de réduire l’amplitude de la marée dans le rapport de à 1, ensuite d’amener pour la marée un retard de phase, un décalage, égal à

Si l’astre sur lequel se produit la marée, au lieu d’être liquide et visqueux, était solide et parfaitement élastique, nous aurions bien une réduction de l’amplitude, mais nous n’aurions pas de décalage. S’il était solide et visqueux, nous aurions à la fois une réduction de l’amplitude plus grande que et un décalage.

107.Reprenons la formule (12) qui définit

(12)

Si la viscosité est très faible, c’est-à-dire si est très petit, l’angle est petit, il se confond sensiblement avec sa tangente et il est proportionnel à c’est-à-dire à la vitesse de la marée ; de même et seront proportionnels à quant à il sera presque égal à . Au contraire, si la viscosité était très forte, serait très grand, aussi, serait voisin de et l’on aurait sensiblement

dans ce cas, serait très sensiblement égal à  ; , égal alors à , serait proportionnel à et serait aussi proportionnel à .

Dans le cas actuel devons-nous considérer la viscosité comme faible ou comme forte ? Nous devons la considérer comme faible, car une substance comme la poix nous donnerait, si nous faisions le calcul, un angle très petit et non voisin de . La poix doit donc être considérée, au point de vue qui nous occupe, comme un corps à viscosité très faible, et il y a lieu d’admettre que la Terre, à l’époque où elle était liquide, n’offrait pas une viscosité incomparablement supérieure à celle de la poix.

108.Revenons à la formule de la dénivellation statique

(9)

et indiquons maintenant comment on pourrait développer effectivement son troisième membre sous forme trigonométrique.

Représentons sur une sphère le pôle boréal P, le lieu géographique A et l’astre perturbateur L, de façon à former un triangle sphérique PAL (triangle Pôle-Zénith-Astre) dont les trois côtés sont

PA  colatitude du lieu,
PL complément de la déclinaison de l’astre,
LA distance zénithale géocentrique de l’astre.

L’angle APL est l’angle horaire de l’astre relativement au point A où on étudie la marée. La formule fondamentale de la trigonométrie sphérique, appliquée au triangle sphérique PAL, nous donne l’équation suivante :

d’où nous tirons

Si nous portons cette valeur dans la formule (9), nous voyons que contient trois sortes de termes[6] :

1o Des termes indépendants de l’angle horaire  : ils donneront les marées à longue période ;

2o Des termes en dont la période est voisine d’un demi-jour : ils donneront les marées semi-diurnes ;

3o Des termes en dont la période est voisine d’un jour : ils correspondent aux marées diurnes.

Introduisons les éléments osculateurs de l’orbite de l’astre L. Soient, sur la sphère céleste (fig. 30), E le plan de l’équateur terrestre, O le plan de l’orbite de l’astre troublant L, N le nœud de l’orbite sur Figure 30
fig.30.
l’équateur, A0 l’origine des longitudes géographiques (intersection de l’équateur avec le méridien de Paris), A′ l’intersection de l’équateur avec le méridien du lieu géographique A, le périgée de l’astre. Nous avons

A′A0  longitude du lieu,
A0N angle horaire du nœud à Paris.

Nous désignons par l’arc N égal à la longitude du périgée moins la longitude du nœud ; et par la longitude moyenne de l’astre dans son orbite, comptée à partir du nœud N.

Alors, si nous développons le troisième membre de la formule (9), autrement dit en série de fonctions trigonométriques dont l’argument est une fonction linéaire du temps , nous obtiendrons des termes tels que

(13)

étant une fonction des coordonnées du lieu A, et étant des entiers. C’est en effet des quantités

que dépendent les coordonnées horaires de l’astre et

Puisque dans l’expression (13), ne dépend comme que des coordonnées du lieu A, et non de l’astre, nous l’isolerons en développant le cosinus sous la forme

(14)

Et finalement les différents termes du développement trigonométrique de seront de la forme

(15)

où la constante vaut ou , suivant que l’on prend le cosinus ou le sinus qui figure dans la formule (14), et où nous désignons par le coefficient qui dépend des coordonnées du lieu géographique.

De quelle nature seront ces coefficients en tant que fonctions des coordonnées du lieu ? Ils seront évidemment des fonctions sphériques du second ordre comme l’est lui-même le troisième membre de la formule (9). Nous avons désigné chacun de ces coefficients par un produit  : la lettre désigne une constante numérique, et la lettre une fonction sphérique du second ordre multipliée par un nombre convenable de telle façon que l’intégrale

étendue à tous les éléments de surface de la sphère, ait une même valeur constante pour toutes les fonctions sphériques que nous aurons à envisager.

Soient les coordonnées du lieu par rapport à trois axes rectangulaires invariablement liés à la Terre et passant par son centre, l’axe des étant l’axe de rotation, le plan des étant le méridien de Paris. Nous aurons (en prenant pour unité le rayon terrestre )

Les fonctions sphériques qu’introduit alors le développement trigonométrique du troisième membre de la formule (9) sont les suivantes[7] :

Faisons maintenant quelques remarques sur les trois entiers Tout d’abord on aura

0, 1 ou 2,

suivant que le terme considéré correspond à une marée à longue période, à une marée diurne, ou à une marée semi-diurne ; car, dans l’expression (15), et ne varient que très lentement et peuvent être regardés comme sensiblement constants. Ensuite, si l’inclinaison est nulle, on aura

car alors les deux plans E et O (fig. 30) coïncident, le point N devient indéterminé, et les seuls angles qui interviennent sont et . Par suite, dans les termes indépendants de l’inclinaison, on aura

Dans les termes qui contiennent l’inclinaison à la puissance , on verrait facilement que

doit être au plus égal à et de même parité que . Enfin, dans les termes indépendants de l’excentricité, n’intervient pas, on a donc

et dans les termes qui contiennent l’excentricité à la puissance on verrait que est au plus égal à et de même parité que .

Comme nous ne conserverons dans la suite que des termes du premier ordre au plus par rapport à l’excentricité et à l’inclinaison , nous aurons, dans les termes indépendants de et de

dans les termes en ,

dans les termes en ,

109.Bref, nous avons prouvé que le troisième membre de la formule (9), qui représente (au facteur près) le potentiel générateur de la marée, peut se développer en une série de termes de la forme (15) :

ce que nous écrivons simplement

en posant pour abréger

Si nous voulons maintenant tenir compte de la viscosité comme nous l’avons fait au no 106 (p. 146), nous trouverons que ce potentiel perturbateur produit une dénivellation

(16)

est donné par la formule

représente ici la vitesse de la marée que nous appelions au no 106)[8].

Nous nous proposons de chercher l’action de la Terre, ainsi déformée par la marée (16), sur un corps extérieur. Pour fixer les idées, nous supposerons que la marée (16) est produite par le Soleil et nous chercherons les perturbations que cette marée solaire (16) fait subir au mouvement de la Lune.

Les quantités sont donc relatives au Soleil. Nous appellerons les mêmes quantités relatives à la Lune. Comme nous cherchons l’action, sur l’orbite de la Lune, du bourrelet soulevé par la marée solaire à la surface de la Terre, nous introduirons une fonction perturbatrice qui sera le potentiel dû à l’attraction de ce bourrelet. Figure 31
fig.31.
Soit la distance de la Lune L à l’élément du bourrelet (fig. 31). Nous aurons alors

l’intégrale étant étendue à toute la surface de la sphère terrestre.

Nous pouvons développer comme nous avons développé au no 104 (p. 144) et écrire

étant respectivement des fonctions sphériques d’ordre 0, 1, 2, par rapport aux coordonnées du lieu géographique A. Comme est lui-même une fonction sphérique du second ordre par rapport aux coordonnées du même point A, on a

Nous avons donc simplement pour notre fonction perturbatrice

(17)

Développons sous forme trigonométrique de la même manière que nous avons développé un peu plus haut : nous aurons

étant trois entiers, étant une constante égale à 0 ou à étant un coefficient numérique, et une fonction sphérique du second ordre telle que l’intégrale

étendue à toute la sphère ait une valeur constante donnée la même pour toutes les fonctions sphériques Nous écrirons simplement

(18)

en désignant, pour abréger, par

l’argument du cosinus. Alors, d’après (16) et (18), l’expression (17) de peut s’écrire, en faisant sortir du signe tout ce qui ne dépend pas des coordonnées du lieu A,

(19)

Telle est l’expression de la fonction perturbatrice dont nous avons à chercher l’action sur l’orbite de la Lune.

110.Nous appliquerons la méthode de la variation des constantes. Nous commencerons par définir la position de la Lune par un système de six éléments canoniques[9]. Aux trois quantités nous adjoindrons, pour achever de déterminer la position de la Lune, les trois suivantes :

représentant le demi-grand axe, l’excentricité et l’inclinaison sur l’équateur de l’orbite lunaire.

Le vecteur des aires a alors pour valeur (à un facteur constant près dépendant des masses et dont nous faisons abstraction)

et la projection de ce vecteur sur la perpendiculaire au plan de l’équateur a pour valeur

Les six éléments

forment un système de variables canoniques[10] entre lesquelles existent les équations de Hamilton. Nous aurons en particulier

(20)

Telles sont les équations auxquelles nous conduit l’application de la méthode de la variation des constantes, pour calculer les perturbations des éléments lunaires sous l’action de la fonction perturbatrice

111.Le calcul qui précède est relatif à l’action, sur l’orbite lunaire, du bourrelet liquide soulevé sur les océans terrestres par la marée solaire. Ne pourrait-on pas appliquer le même calcul à l’action, sur l’orbite lunaire, du bourrelet liquide soulevé par la marée lunaire elle-même ? On le peut certainement, mais à condition de prendre quelques précautions : devenant alors égaux à , la fonction perturbatrice se trouve dépendre de la variable de deux manières différentes ; elle en dépend par et ensuite par . Dans le calcul de il faut donc supposer d’abord puis dériver par rapport à et enfin faire Les mêmes précautions doivent être prises dans le calcul des autres dérivées de En d’autres termes, il faut distinguer la Lune en tant qu’astre troublant producteur des marées et en tant qu’astre troublé par ces marées. C’est ainsi que Sir G. H. Darwin appelle notre satellite Diane quand il est troublant et Lune quand il est troublé : alors sont les coordonnées de Diane, , sont les coordonnées de la Lune, et l’on a

Moyennant cette précaution, notre analyse s’applique à l’action des marées lunaires sur la Lune elle-même,

112.Nous reprendrons donc les formules (20), et comme nous ne nous occupons que des effets séculaires, nous ne conserverons aux seconds membres que les termes constants indépendants du temps Nous avons d’après la formule (19)

(20bis)

Dans ces expressions, le temps ne figure que dans le produit

produit qui, développé, fournit des termes en

et des termes en

Nous aurons donc, aux seconds membres des expressions (20), un terme indépendant de chaque fois que ou se réduira à une constante. Comme n’est jamais une constante, les termes séculaires correspondent à

c’est-à-dire à

(21)

Si l’astre troublant et l’astre troublé sont différents, cette condition n’est jamais remplie[11], car le rapport des moyens mouvements de ces deux astres n’est pas commensurable. Nous pouvons donc dire que la marée solaire ne produit pas d’effet séculaire sur l’orbite lunaire.

Mais si l’astre troublant et l’astre troublé sont identiques, la condition (21) est remplie chaque fois que et n’est remplie que dans ce cas : on a alors

qui donne le terme constant

Les seconds membres des formules (20bis), dans lesquels on ne conserve que les termes constants indépendants de , deviennent donc

Nous pouvons maintenant effacer les accents « prime » devenus inutiles, puisque tout se rapporte au même astre, la Lune. Nous remarquons que le facteur

est le même pour tous les termes, puisque nous avons eu soin de choisir les fonctions sphériques de façon que cette intégrale ait toujours la même valeur constante (note de la p. 151). Nous n’écrivons donc pas ce facteur, qui est le même partout (quitte à le rétablir plus tard), et les formules (20bis) prennent la forme

(20ter)

Rappelons-nous que l’angle est très petit et que par suite est proportionnel à , c’est-à-dire à la vitesse de la marée . Or,

et l’on a

 vitesse angulaire de rotation de la Terre.
vitesse angulaire de rotation de la Lune.
à très peu près.

Par suite, dans les formules (20ter), nous pouvons remplacer par

qui lui est proportionnel.

Finalement les formules (20), dans lesquelles on ne conserve que ce qui est relatif aux termes séculaires, sont devenues

(21)

Ces dernières formules, par un calcul immédiat, nous fournissent les suivantes :

(22)

Ces trois formules vont nous permettre de calculer les variations séculaires du grand axe, de l’excentricité et de l’inclinaison de l’orbite lunaire, provoquées par le bourrelet liquide de la marée lunaire elle-même. Dans la troisième de ces formules (22), qui nous servira à calculer la variation de l’inclinaison nous négligerons les termes en car

s’annule avec l’excentricité, et, dans le calcul relatif à l’inclinaison, nous pouvons, en première approximation, supposer l’excentricité nulle.

113.Les seconds membres de (22) sont des sommes de termes dont chacun provient d’un des termes

(15)

en lesquels on a pu décomposer le potentiel producteur de la marée lunaire. Parmi ces termes, nous ne conserverons ici que ceux qui sont du premier ordre au plus par rapport à l’excentricité et à l’inclinaison . La théorie des marées[12] nous enseigne que ces termes sont les suivants :

Termes semi-diurnes :
(principal semi-diurne),
(elliptique majeur),
(elliptique mineur)
Termes diurnes :
(principal diurne),
(sidéral) ;
Terme à longue période :
(elliptique mensuel),

Ces six termes sont ceux que nous conserverons dans , (seul le premier est indépendant de et de ),

114.Tels qu’ils sont, ces six termes ne sont pas tout à fait de la forme (15) (puisque figure dans l’argument) ; or, il convient de les rendre de cette forme.

1o Les trois termes semi-diurnes sont de la forme

(sd)

qui peut s écrire

or,

et

sont justement des fonctions sphériques qui donnent à l’intégrale

la valeur (note de la p. 151). Donc chaque terme semi-diurne (sd) fournit deux termes de la forme (15) ayant respectivement pour coefficient et , et pour argument et . Par suite, chaque terme semi-diurne fournit aux seconds membres des formules (22) deux termes égaux, pour lesquels on a

2o Les deux termes diurnes sont de la forme

(d)

ce qui peut s’écrire

or,

et

sont justement des fonctions sphériques qui donnent à l’intégrale

la valeur . Donc chaque terme diurne (d) fournit deux termes de la forme (15) ayant respectivement pour coefficient et , et pour argument et . Par suite, chaque terme diurne fournit aux seconds membres des formules (22) deux termes égaux, pour lesquels on a

3o Le terme à longue période est de la forme

et, comme

est une fonction sphérique donnant à l’intégrale

la valeur , ce terme est de la forme (15), ayant pour coefficient . Il donnera donc aux seconds membres des formules (22) un seul terme pour lequel on aura

Comme chaque terme diurne ou semi-diurne en fournit deux égaux dans les seconds membres des formules (22), tandis que le terme longue période n’en fournit qu’un seul, il y a lieu de multiplier par 2 la quantité relative à chaque terme diurne ou semi-diurne. Au lieu de cela, nous diviserons par 2 la quantité relative au terme à longue période[13], pour lequel on devra prendre, par conséquent,

115.Bref, si nous prenons les six termes

nous pouvons former le tableau suivant à double entrée :


Pour le calcul des quantités , a été remplacé par et par . Dans les deux dernières colonnes on n’a conservé les termes en et en que pour le terme principal (). Enfin, dans l’avant dernière colonne, on ne s’est pas occupé des termes et qui contiennent en facteur ; dans la dernière colonne, on ne s’est pas occupé des termes et qui contiennent en facteur ; car, pour le calcul des variations de l’excentricité, nous supposons l’inclinaison nulle et, pour le calcul des variations de l’inclinaison, nous supposons l’excentricité nulle.

116.Le Tableau précédent fournit tous les éléments nécessaires pour calculer

au moyen des formules (22). Nous trouvons ainsi

Ces trois formules s’écrivent (en négligeant au second membre de la première les termes en et à côté du premier terme qui est fini)

(23)

Rappelons-nous que nous avons fait abstraction, aux seconds membres des formules (20) ou (22), d’un même facteur constant. Ce facteur, il est facile de le voir, est négatif. Désignons-le par et rétablissons-le maintenant dans les formules (23) qui deviennent ainsi

(24)

Au lieu de conserver les variables et , introduisons l’excentricité et l’inclinaison  : nous avons, à des termes près d’ordre supérieur,

et par suite

Les équations (24) s’écrivent donc

(25)

Ces dernières équations donneront les variations du grand axe, de l’excentricité et de l’inclinaison de l’orbite lunaire par suite de l’effet des marées.

117.Nous voulons aussi calculer la variation de la vitesse angulaire de rotation de la Terre. Nous nous servirons pour cela du théorème des aires. Dans la Section 'II' de ce Chapitre, nous avons posé

et l’équation des aires nous a donné (no 98)

(3)

En vertu de la troisième loi de Képler, est proportionnel à la racine carrée du grand axe, c’est-à-dire à ce que nous appelons ici Donc, en adoptant des unités convenables, cette équation (3) peut s’écrire

Mais cette équation n’a été établie qu’en supposant nulle l’inclinaison de l’orbite lunaire sur l’équateur, hypothèse que nous abandonnons ici. Représentons (fig. 32) par O le plan de l’orbite lunaire, Figure 32
fig.32.
par E celui de l’équateur terrestre, et par Π le plan invariable (plan du maximum des aires).

Modifiant un peu les notations précédentes, nous appellerons l’inclinaison de l’orbite sur le plan invariable et l’inclinaison de celui-ci sur l’équateur.

Le théorème des aires donne alors les deux équations suivantes :

Les angles et étant supposés petits, ces équations s’écrivent approximativement

(26)

Ces équations nous serviront à calculer les variations des éléments terrestres et puisque les variations des éléments lunaires et sont données par les équations (25). Nous aurons d’abord

(27)

Ensuite, puisque nous appelons maintenant ce qui, dans les équations (25), est désigné simplement par nous devons, dans la troisième équation (25), remplacer par ( n’étant pas changé) ; cette équation devient

(28)

Alors la seconde équation (26) donne

d’où nous tirons

(29)

118.Réunissons dans un tableau les deux premières équations (25) et les équations (27), (28) et (29) :

(30)

C’est de ces équations que nous allons tirer les variations des éléments.

119.Remarquons d’abord que les deux premières ne contiennent ni ni . Elles donnent les variations de la rotation terrestre et de la distance moyenne de la Lune. Ce sont ces variations que nous avons discutées en détail dans la Section II. Rappelons les résultats de cette Figure 33
fig.33.
discussion. Prenant pour abscisse et pour ordonnée (fig. 33), nous traçons la droite

et la courbe de rigidité

qui coupe la droite en deux points C et D. Le point représentatif du système Terre-Lune est parti de C, état dans lequel le jour et le mois étaient tous deux égaux à 5h36m ; actuellement, ce point représentatif est en P et l’état final sera représenté par D, où le jour redeviendra égal au mois, leur durée commune étant de 55 jours actuels.

120.Passons maintenant à la variation de l’excentricité donnée par la troisième équation (30). Il s’agit de savoir si est positif ou négatif, pour reconnaître si l’excentricité croit ou décroît. Or, le signe de dépend du signe de l’expression

Si

l’excentricité croît. Si, au contraire,

l’excentricité décroît. Traçons sur la figure 33 la courbe

c’est-à-dire la courbe

Cette courbe (représentée en trait ponctué) coupe la droite AE en deux points C′ et D′.

Lorsque le point représentatif () parcourait le segment de droite CC′, on avait

c’est-à-dire

l’excentricité a donc commencé par décroître.

En P (état actuel), on est entre C′ et D′, par suite

et l’excentricité est en train de croître.

Enfin, lorsque le point représentatif parcourra D′D, l’excentricité recommencera à décroître.

Si l’on trace la courbe (fig. 34) qui représente les variations de l’excentricité en fonction de (ce qui est possible puisqu’on connaît l’état actuel P), on constate que cette courbe présente une asymptote verticale correspondant à l’abscisse du point C et que l’excentricité passe par un minimum en C′, par un maximum en D′, puis décroît ensuite jusqu’au point final D où elle s’annule.

Il ne faudrait pas croire que l’asymptote verticale signifie que l’excentricité a été initialement très grande. Les équations (30), en effet, supposent essentiellement très petit et cessent d’être applicables dès que devient grand.

D’ailleurs, ce n’est pas ainsi que le problème se pose. Supposons que l’excentricité ait été initialement nulle, la troisième équation (30) donnerait

donc

est une solution : elle correspond à une orbite circulaire restant indéfiniment circulaire[14]. L’important est de savoir si cette solution est Figure 34
fig.34.
stable ou instable. Supposons que ait subi une petite variation : la troisième équation (30) étant de la forme

la solution

sera stable ou instable, suivant que , c’est-à-dire

sera négatif ou positif.

Donc, tant que le point représentatif (fig. 33) est situé entre C et C′, la solution

est stable ; l’orbite reste circulaire. Mais à partir du point C′ cette solution devient instable, et l’excentricité, ayant cessé d’être nulle par suite d’une petite perturbation quelconque, croît jusqu’en D′ ; puis, à partir de là, devenant négatif, elle décroît jusqu’en D où elle s’annule.

La courbe représentative de l’excentricité (fig. 34) se composerait ainsi de la portion de droite cC′ et du morceau de courbe C′D′D.

Par suite, le fait qu’actuellement l’orbite lunaire est excentrique n’implique pas forcément qu’à l’origine l’excentricité était différente de zéro : le frottement des marées a pu, d’après Sir G. H. Darwin, faire naître une excentricité qui n’existait pas initialement.

121.Étudionsmaintenant les variations des inclinaisons et , données par les deux dernières équations (30). Comme ce qui nous intéresse c’est l’angle que fait l’orbite avec l’équateur, nous ajoutons ces deux équations : il vient

(31)

Donc croîtra ou décroîtra suivant le signe de la quantité

Remplaçant par cette quantité s’écrit (au facteur positif près)

ce polynôme en présente deux variations de signe ; il a donc au plus deux racines positives. Dans le cas de la Lune, il a effectivement deux racines positives qui correspondent aux abscisses de deux certains points C′′ et D′′ situés entre C et D (fig. 33).

Nous pouvons donc faire pour l’inclinaison la même discussion que pour l’excentricité. Nous avons pour l’équation (31) la solution

Cette solution est stable lorsque le point représentatif est entre C et C′′ ; elle devient instable entre C′′ et D′′. Si donc nous partons d’un état initial où l’inclinaison est nulle, l’inclinaison restera nulle au début ; puis, lorsque le point représentatif sera arrivé en C′′, si elle cesse d’être nulle par suite d’une petite perturbation quelconque, elle augmentera jusqu’en D′′ ; ensuite elle diminuera jusqu’en D où elle s’annulera de nouveau.

Nous voyons donc que le frottement des marées a pu faire naître une inclinaison du plan de l’orbite sur l’équateur, qui n’existait pas initialement.

122.Nous venons d’exposer, d’après Sir G. H. Darwin, l’évolution passée du système Terre-Lune. Mais quelle a pu être la durée de cette évolution ? Nos formules ne nous l’apprennent pas, car il y entre le coefficient de viscosité inconnu qu’avait la Terre quand elle était encore pâteuse. Néanmoins nous pouvons reconnaître facilement qu’il y a un maximum de l’action perturbatrice, c’est-à-dire un minimum de temps nécessaire.

Les seconds membres des formules (20ter) sont de la forme

(32)

Les et les sont connus par la théorie des marées. Mais les , définis par les équations

sont inconnus puisque dépend de la viscosité. Or, un sinus est toujours inférieur à l’unité en valeur absolue. Nous exagérerons donc l’expression (32), par suite nous diminuerons le temps nécessaire à la variation de , si nous remplaçons les par l’unité et si nous prenons tous les termes avec le même signe.

Époques Durée du jour  Durée du mois (en jours actuels)  Incli- naison  Inverse de l’aplatis- sement  Distance (en rayons terrestres actuels)  ξ/n Chaleur dégagée
0 23h56m 27,32 23°28' 232 60,4 '4,01' 0
46 300 000 15°30° 18,62 20°40' 96 46,8 '2,28' 225
56 600 000 9°55° 8,17 17°20' 40 27,0 '1,11' 760
56 800 000 7°50° 3,59 15°30' 25 15,6 '0,67' 1300
56 810 000 6°45° 1,58 14°25' 18 9,0 '0,44' 1760

C’est de cette façon qu’a procédé Sir G. H. Darwin pour introduire le temps dans son analyse. Il a dressé le Tableau ci-dessus, dans lequel la première colonne indique le nombre d’années écoulées à partir de l’époque actuelle, et en arrière (ces époques sont de plus en plus rapprochées à mesure qu’on remonte dans le passé, parce que, la Lune étant alors plus proche de la Terre, les effets des marées étaient plus considérables).

L’énergie mécanique du système Terre-Lune a toujours été en diminuant, le frottement la transformant en chaleur. Si cette chaleur avait été employée à échauffer la Terre, elle aurait élevé sa température d’un certain nombre de degrés (Fahrenheit) : c’est ce nombre de degrés qui figure à la dernière colonne du Tableau. Sir G. H. Darwin fait remarquer que l’on pourrait peut-être invoquer cette cause pour expliquer l’origine de la chaleur interne du globe.

123.Sir G. H. Darwin a aussi essayé de calculer une valeur du coefficient de viscosité inconnu de la Terre, en partant de l’accélération séculaire de la Lune. Des observations d’éclipses dans l’antiquité ont permis d’évaluer cette accélération à 10″. Or, le calcul indique une accélération théorique de 6″ seulement. Il y a donc une accélération de 4″ que la gravitation n’explique pas[15]. Sir G. H. Darwin cherche, comme le fit autrefois Delaunay (no 95), à en rendre compte par l’augmentation du jour sidéral, due au frottement des marées. La Lune s’éloignant de la Terre subit réellement, non pas une accélération, mais un retard, et son accélération apparente ne serait due qu’à la différence entre le retard réel de la rotation terrestre et le retard réel de la révolution lunaire, le premier de ces deux retards étant plus grand que l’autre. C’est en égalant à 4″ la différence de ces deux retards, telle que la lui donnent ses formules, que Sir G. H. Darwin a calculé le coefficient de viscosité de l’intérieur de la Terre. Si l’on adoptait le coefficient ainsi obtenu, on trouverait pour la durée de l’évolution plusieurs milliards d’années ; mais il importe d’observer que la Terre a pu être autrefois beaucoup plus liquide qu’aujourd’hui.

124.Jusqu’ici, nous avons toujours adopté des unités particulières destinées à simplifier les formules. Cela était légitime, parce que nous étudiions l’action des marées produites par un même astre, la Lune, sur un même astre, la Terre. Mais si nous voulons maintenant comparer l’action des marées produites par divers astres sur un même astre, ou par un même astre sur divers astres (par exemple si nous envisageons le système formé par une planète et plusieurs satellites, ou bien le système formé par le Soleil et l’ensemble des planètes), il faudra rétablir les coefficients de proportionnalité. C’est ce que nous allons faire.

125.De même que la marée lunaire a une influence sur la longueur du jour et du mois, la marée solaire a une influence sur la longueur du jour et de l’année. Désignons par les masses du Soleil, de la Lune, de la Terre. Nous appelons toujours

le moment d’inertie de la Terre,
sa vitesse angulaire de rotation,
la vitesse de révolution de la Lune,
le demi-grand axe de son orbite.

L’équation des aires, appliquée au système Terre-Lune, supposé seul, s’écrit

(33)

représente le moment de rotation dû à la rotation terrestre, le moment de rotation dû au mouvement orbital de la Lune. Le moment de rotation dû au mouvement orbital de la Terre autour du centre de gravité du système Terre Lune étant très petit, nous le négligeons à côté de . Quant à celui qui est dû à la rotation de la Lune sur elle-même, nous en faisons abstraction, envisageant ici la Lune comme un simple point matériel sans dimensions.

Écrivons maintenant l’équation des aires pour le système Terre-Soleil, supposé seul. Appelant

la vitesse de révolution de la Terre autour du Soleil,
le demi-grand axe de l’orbite terrestre,

nous aurons

(34)

représente toujours le moment de rotation dû à la rotation terrestre ; représente le moment de rotation dû au mouvement orbital de la Terre autour du Soleil, à côté duquel celui qui est dû au mouvement orbital du Soleil, autour du centre de gravité du système Terre-Soleil, est négligeable. Nous faisons aussi abstraction du moment de rotation dû à la rotation du Soleil sur lui-même, envisageant ici le Soleil comme un simple point matériel, ou, si l’on préfère, comme une sphère absolument rigide dont le moment de rotation est constant.

Transformons les équations (33) et (34), en introduisant les quantités

et en nous servant de la troisième loi de Képler qui donne

Les équations des aires (33) et (34) s’écrivent alors

Considérons d’abord l’action de la marée lunaire. Cette marée donne, pendant un certain temps , à la rotation de la Terre une variation , et il en résulte pour une variation donnée par

(35)

De même, la marée solaire fait, dans le même temps , subir à une variation et il en résulte pour une variation donnée par

(36)

Comparons à  : il s’agit de l’action de deux marées différentes sur un même astre, la Terre. Le rapport de à sera donc approximativement égal au rapport des carrés des coefficients des deux marées. La marée solaire est environ trois fois moindre que la marée lunaire. On a donc

L’augmentation de la durée du jour provenant de la marée solaire est donc environ 9 fois moindre que celle qui provient de la marée lunaire.

Les équations (35) et (36) donnent ensuite

or, le second membre est extrêmement petit ; est donc incomparablement plus petit que  ; ce qui signifie que la marée solaire n’a qu’une influence insignifiante sur la distance moyenne de la Terre au Soleil, c’est-à-dire sur la longueur de l’année.

126.Les marées que la Lune soulève sur la Terre ralentissent la rotation terrestre. De même les marées que la Terre pourrait soulever sur la Lune exerceraient une influence sur la rotation de la Lune. Actuellement, la Terre ne soulève pas de marées sur la Lune, puisque la Lune nous présente toujours le même hémisphère. Mais aux époques reculées, la Terre a dû soulever des marées sur son satellite, et c’est précisément, ainsi que l’avait déjà affirmé Laplace, à l’action de ces marées qu’est due l’égalité actuelle de la durée de rotation et de la durée de révolution de la Lune.

Soit la masse de l’astre perturbateur qui produit la marée sur un astre sphérique T de rayon (fig. 35). Appelons la distance Figure 35
fig.35.
TM du centre de la sphère T à l’astre perturbateur et la distance MA de l’astre perturbateur à un point A de la surface de l’astre T. Le potentiel dû à l’astre perturbateur est

Le potentiel générateur de la marée est

il est proportionnel à

La dénivellation statique est

elle est proportionnelle à

Le bourrelet liquide, dû à cette dénivellation, produit sur l’astre M un potentiel perturbateur

désigne la densité du bourrelet liquide soulevé sur l’astre T. Or, nous avons, le signe indiquant la proportionnalité,

et peut, sous le signe être remplacé par [16] qui est proportionnel à Nous pouvons donc écrire

Le couple qui fait varier la rotation de l’astre T est proportionnel à

puisque

représente le travail virtuel produit dans une petite rotation de la sphère T.

Comme est lui-même proportionnel à [17] le couple est proportionnel à

Si, maintenant, nous tenons compte de la viscosité de l’astre T, nous trouvons

couple perturbateur

La dérivée de la vitesse de rotation de l’astre T s’obtiendra en divisant le couple perturbateur par le moment d’inertie de la sphère T qui est proportionnel à

Bref, nous aurons


L’angle est défini par l’équation

étant une constante proportionnelle à la viscosité et étant proportionnel à la vitesse de la marée. Pour la marée principale semi-diurne, on trouve, en faisant le calcul,

étant le coefficient de viscosité de l’astre T. Comme l’angle est très petit (en prenant pour le coefficient de viscosité de la poix, ou trouve que cet angle un dépasse pas 10″ pour le système Terre-Lune), nous pouvons écrire

Or, nous avons pour la gravite , à la surface de l’astre T

et par suite

Finalement, il vient

Si la sphère T est la Terre et l’astre M la Lune, ces proportionnalités représentent la diminution de la rotation terrestre due à la marée lunaire.

Si, maintenant, nous voulions avoir la diminution de la rotation de la Lune, due aux marées que la Terre produit sur la Lune, nous aurions

désigne cette fois la masse de la Terre, et la densité et le rayon de la Lune, sa vitesse angulaire de rotation.

Actuellement, les durées de rotation et de révolution de la Lune étant les mêmes, est nul et il en est de même de l’action retardatrice de la Terre sur la Lune est épuisée. Mais il n’en a pas toujours été ainsi : à l’origine était comparable à , et alors on avait

ce qui, puisque

s’écrit

Comme les densités et de la Terre et de la Lune sont du même ordre, ce rapport est comparable à

32 000 environ.

Donc l’action retardatrice de la Terre sur la Lune a dû être au début environ 32 000 fois plus forte que l’action retardatrice de la Lune sur la Terre. On conçoit donc que cette action ait pu être assez puissante pour avoir arrêté le mouvement de rotation de la Lune relativement à la Terre, et pour avoir forcé notre satellite à nous tourner toujours la même face.

127.Si, maintenant, on étudie le système formé par le Soleil et les différentes planètes, ou constate d’abord que l’action produite sur la rotation du Soleil par les marées que les planètes soulèvent sur le Soleil est tout à fait insignifiante. Quant à l’action retardatrice de la marée solaire sur la rotation des différentes planètes, son coefficient de proportionnalité est

ce qui peut s’écrire

en négligeant la vitesse angulaire de révolution de la planète, à côté de sa vitesse angulaire de rotation.

Or, , masse du Soleil, est le même pour toutes les planètes, et l’on a

étant la masse de la planète envisagée. Le coefficient de proportionnalité auquel est proportionnel est donc pour chaque planète

Sir G. H. Darwin donne le Tableau suivant (où représente le moment d’inertie de la planète) :

Planètes
Mercure 
9,10000 1 000,00002
Vénus 
8.1 11,00002
La Terre 
1 1,00002
Mars 
0,026 0,89002
Jupiter 
2,3 0,00005
Saturne 
0,11 0,00002

Le calcul, en ce qui concerne Mercure et Vénus, a été fait en supposant à ces deux planètes une durée de rotation de 24 heures. On voit que, pour ces planètes intérieures, l’action retardatrice de la marée solaire est forte. S’il est vrai, comme le prétendent plusieurs observateurs, que ces deux planètes tournent toujours une même face vers le soleil[18], ce fait pourrait s’expliquer par l’action retardatrice de la marée solaire. La théorie expliquerait aussi pourquoi les planètes extérieures Uranus et Neptune, n’ayant subi qu’une marée solaire très faible, ont conservé leur rotation rétrograde.

IV. Influence accélératrice du refroidissement.

128.Le refroidissement séculaire de la Terre la contracte et diminue son moment d’inertie ; par suite, en vertu de la loi des aires, sa vitesse de rotation doit s’accroître de ce chef. Il y a donc lieu, à côté de l’influence retardatrice des marées, étudiée dans les pages précédentes, de tenir compte de l’influence accélératrice du refroidissement.

Nous avons vu plus haut (no 126) que le couple retardateur dû à la marée est proportionnel à

c’est-à-dire à

ou encore à

Or, (masse de la Lune) et (masse de la Terre) sont des constantes, et nous regarderons aussi (distance de la Lune à la Terre) comme constant. Par suite, le couple est proportionnel à

Le moment de rotation d