Formule de MacCullagh - Définition
La liste des auteurs de cet article est disponible ici.
Introduction
La formule de MacCullagh est une importante formule pour décrire le champ de gravité d'un corps non sphérique à une distance suffisamment grande mais pas très grande par rapport aux dimensions linéaires du corps. Elle s'obtient en partant de la définition de base du potentiel gravifique, sans effectuer un développement général de Laplace dont elle représente un cas particulier limité à une approximation d'ordre 2. Elle est très utilisée en géodésie théorique et en géophysique.
Expression générale du potentiel gravifique
En coordonnées sphériques, désignons la composante radiale du vecteur-position
avec
La constante G dénote la constante de gravitation de Newton.
Comme (1 + x) − 1 / 2 = 1 − x / 2 + 3x2 / 8 + o(x3) lorsque
![\frac{1}{\sqrt{1 - 2 \frac{r'}{r} \cos \psi + \left(\frac{r'}{r}\right)^2}} = 1 + \frac{r'}{r} \cos \psi - \frac{1}{2} \left(\frac{r'}{r}\right)^2 + \frac{3}{2} \left(\frac{r'}{r}\right)^2 \cos^2\psi + o\left[\left(\frac{r'}{r}\right)^3\right]](https://static.techno-science.net/illustrationWebp/Definitions/autres/f/f9f0c3575aadd392776974e4ac962116_ea1d88bde00087b77d75fa7e06af98b5.png)
Par conséquent, on a
Le potentiel gravifique en P est donc
Formule de MacCullagh pour les corps à symétrie de révolution
Le moment d'inertie I autour d'un axe passant par le centre de masse O d'un corps peut s'écrire en termes des moments principaux d'inertie A, B, C obtenus par diagonalisation préalable de la matrice 3x3 représentant le tenseur d'inertie au moyen de la relation
Ici, nx, ny, nz désignent les cosinus directeurs de l'axe OP par rapport aux axes principaux d'inertie Ox, Oy, Oz, respectivement.
En faisant l'hypothèse que le corps est de révolution autour de l'axe Oz, c'est-à-dire qu'il possède une symétrie axiale d'axe Oz, on a A = B. Désignons par φ l'angle que fait la direction OP avec le plan Oxy. Alors nz = sinφ et, à cause de la relation liant
Avec A = B, la formule de MacCullagh devient
Dans le cas d'une symétrie axiale, le moment d'inertie I prend lui-même la forme
La formule de MacCullagh, pour un corps de révolution, prend donc finalement la forme suivante :
Théorème de MacCullagh
On voit ainsi qu'à l'approximation d'ordre 2, le potentiel gravifique V(P) peut se décomposer en quatre termes distincts V0(P), V1(P), V2(P), V3(P), soit
Le premier terme,
correspond au potentiel d'une distribution sphérique, lequel est le même que celui d'une masse ponctuelle égale à la masse comprise dans toute la sphère. A grande distance, c'est évidemment le terme dominant, puisque sa décroissance en 1/r est plus lente que celle des autres termes. C'est le potentiel newtonien classique utilisé pour établir les trajectoires keplériennes des planètes. On dit de lui qu'il correspond à un monopôle, ou à une distribution de masse monopolaire.
Le deuxième terme,
correspond à une distribution de masse dipolaire, autrement dit à un dipôle. En choisissant l'origine des coordonnées au centre de masse du corps, on annule ce terme. En effet, en se référant à la figure ci-contre, on peut écrire
où x0 désigne la composante le long d'un axe Ox de la position du centre de masse.
Les troisième et quatrième termes, V2 et V3, se rapportent à une distribution de masse quadrupolaire, autrement dit à un quadrupôle. Pour V2 on trouve successivement
![V_2 = G r^{-3} \int_M r'^2 \, \mathrm dM = G r^{-3} \int_M (x^2+y^2+z^2) \mathrm dM = \frac12 G r^{-3} \int_M [(y^2+z^2)+(x^2+z^2)+(x^2+y^2)] \, \mathrm dM = \frac12 G r^{-3} (A+B+C)](https://static.techno-science.net/illustrationWebp/Definitions/autres/a/a9f9b8d97aa7003db89def483a3651d3_05da227d0b6d711785514a1cf0dff980.png)
où A, B, C désignent les trois moments d'inertie par rapport aux axes Ox, Oy, Oz, respectivement. En définissant le moment d'inertie moyen
on a donc
Pour V3, on trouve
où I est le moment d'inertie autour de la direction OP. On aboutit ainsi à la formule (parfois aussi appelée théorème) de MacCullagh :
Celle-ci s'applique très utilement à des corps presque sphériques comme la Terre et les planètes, et c'est une approximation encore valable à distance suffisante pour un corps de symétrie quelconque. Il existe évidemment une limitation implicite à sa validité par le fait que la série a été tronquée à l'ordre 2 : des harmoniques supérieurs au degré 2 sont bien sûr nécessaires pour représenter de manière plus précise le champ de gravité extérieur en toute généralité.
