Moment d'inertie

Moment d'inertie

En serrant ses bras le long du corps, cette patineuse diminue son moment d'inertie, augmentant sa vitesse de rotation, puisque son moment cinétique est conservé.

Données clés

Unités SI	kg m2
Dimension	M·L 2
Nature	Grandeur tensorielle extensive
Symbole usuel	I, JΔ
Lien à d'autres grandeurs	$I=\sum _{i}m_{i}r_{i}^{2}$

Le moment d'inertie d'un système physique est une grandeur qui caractérise son inertie vis-à-vis des mouvements de rotation. Sa masse caractérise son inertie vis-à-vis des mouvements de translation. Ce moment d'inertie dépend de la valeur et de la répartition des masses au sein du système. Il a pour dimension ${\textstyle M\cdot L^{2}}$ (produit d'une masse par le carré d'une longueur) et s'exprime donc en kg·m² selon le système international d'unités (S.I.).

Dans le cas de la rotation d'une masse autour d'un axe fixe, le moment d'inertie par rapport à cet axe est une grandeur scalaire qui apparaît dans les expressions du moment cinétique et de l'énergie cinétique de rotation de ce corps. Toutefois, dans le cas général d'une rotation autour d'un axe dont la direction varie au cours du temps, il est nécessaire d'introduire un tenseur symétrique du second ordre, le tenseur d'inertie. Il est toujours possible de choisir un système d'axes, dits axes principaux d'inertie, tel que la matrice représentative de ce tenseur prenne une forme diagonale. Les trois moments correspondants sont les moments principaux d'inertie. Dans le cas particulier d'un solide homogène, ils ne dépendent que de la forme géométrique de celui-ci.

En mécanique des matériaux, l'appellation de « moment d'inertie » est parfois utilisée pour déterminer la contrainte dans une poutre soumise à une flexion. Il s'agit alors d'une notion physique différente, encore appelée moment quadratique, qui a pour dimension L⁴ (exprimée en m⁴ dans le S.I.).

Approche empirique

Résistance au mouvement

Lequel de ces mouvements demande le plus d'énergie ?

Le moment d'inertie est au mouvement de rotation l'analogue de la masse pour le mouvement de translation : il reflète la « résistance » qu'oppose un corps à sa mise en mouvement.

Cette difficulté est d'autant plus grande, dans le cas de la rotation d'un solide, que les masses en son sein se trouvent éloignées de l'axe de rotation. Par exemple, dans le cas d'un balai pris en main au milieu du manche (cf. figure ci-contre), il est plus aisé de le faire tourner autour de l'axe du manche (cas 1), qu'autour de l'axe transversal indiqué (cas 2). En effet, dans ce dernier cas, la brosse, dont la masse relative par rapport au manche est importante, est située plus loin de l'axe de rotation que dans le premier cas (il y a également une asymétrie dans la répartition des masses autour de l'axe de rotation). Comme pour un solide en rotation, la vitesse linéaire d'un point croît en proportion avec cet éloignement, il est nécessaire, à vitesse angulaire égale, de communiquer une plus grande énergie cinétique aux points éloignés. D'où la plus grande « résistance » du balai à tourner autour d'un axe transversal qu'autour de l'axe du manche.

Énergie nécessaire au mouvement

Dans cette analogie, l'énergie nécessaire à la mise en mouvement se transforme (d'une manière ou d'une autre) en énergie cinétique. Dans le cas d'un mouvement de translation, l'énergie cinétique d'un point de masse m avec une vitesse v est donnée par la formule $E_{k}={\frac {1}{2}}mv^{2}$. Dans le cas d'une rotation d'un solide autour d'un axe fixe à la vitesse angulaire ω, il est possible de montrer que l'énergie cinétique totale s'écrit sous la forme $E_{k}={\frac {1}{2}}J_{\Delta }\,\omega ^{2}$, où le terme $J_{\Delta }$ est précisément le moment d'inertie par rapport à l'axe de rotation, c'est-à-dire une forme analogue à celle de l'énergie cinétique du mouvement de translation. Le moment d'inertie apparaît ainsi comme l'analogue de la masse (inerte) du solide dans le cas du mouvement de translation, et est donc en relation avec la « résistance » de celui-ci à sa mise en rotation. Dans l'exemple du balai, le moment d'inertie est plus faible dans le cas d'une rotation autour de l'axe du manche que dans le cas d'une rotation autour de l'axe transversal, d'où une plus grande « résistance à tourner » dans le second cas.

Mise en évidence et définition du moment d'inertie

Soit un solide, considéré comme composé de plusieurs points matériels $M_{i}$ de masse $m_{i}$, dont les distances mutuelles sont fixes. Ce système est en mouvement de rotation autour d'un axe $\Delta$, fixe dans le référentiel d'étude, à la vitesse angulaire $\omega$, qui est la même (à un instant donné) pour tous les points du système. Le moment d'inertie autour de l'axe $\Delta$ apparaît alors naturellement dans les expressions de l'énergie cinétique et du moment cinétique du solide considéré.

Énergie cinétique de rotation

L'énergie cinétique étant une grandeur extensive, c'est-à-dire que sa valeur sur un système complexe est la somme des valeurs sur les parties élémentaires, l'énergie cinétique totale du solide considéré s'exprime donc sous la forme :

$E_{k}=\sum _{i}{\tfrac {1}{2}}m_{i}v_{i}^{2}=\sum _{i}{\tfrac {1}{2}}m_{i}(\omega r_{i})^{2}={\tfrac {1}{2}}\omega ^{2}\sum _{i}m_{i}r_{i}^{2}\ {\overset {\mathrm {def} }{=}}\ {\tfrac {1}{2}}J_{\Delta }.\omega ^{2}$,

avec :

• $r_{i}=\mathrm {H} _{i}\mathrm {M} _{i}$, distance du point $\mathrm {M} _{i}$ à l'axe de rotation $\Delta$ ($\mathrm {H} _{i}$ désignant la projection de $\mathrm {M} _{i}$ sur $\Delta$) ;
• $v_{i}={\frac {dr_{i}}{dt}}$, la vitesse du point $M_{i}$ ;
• $J_{\Delta }=\sum _{i}m_{i}r_{i}^{2}$.

Cette expression possède une analogie évidente avec celle de l'énergie cinétique pour un point matériel (ou un solide en translation), où celle-ci s'écrit $E_{k}={\tfrac {1}{2}}mv^{2}$. Dans le cas de la rotation, la vitesse angulaire ω est l'homologue de la vitesse linéaire v, et le moment d'inertie $J_{\Delta }$ comme l'homologue de la masse m. Toutefois, le moment d'inertie dépend également de la répartition des masses autour de l'axe $\Delta$ et sa valeur dépend donc du choix de celui-ci.

En tant qu'homologue de la masse (inerte) pour la rotation, ce terme reflète la «résistance» du solide à sa mise en rotation, comme il avait été indiqué dans l'approche empirique précédente.

Moment cinétique de rotation du solide

De la même façon que pour l'énergie cinétique, et du fait du caractère extensif de cette grandeur, le moment cinétique du solide par rapport à un point quelconque $\mathrm {O}$ de l'axe s'exprime alors sous la forme, en utilisant les mêmes notations que précédemment :

${\vec {L}}_{\mathrm {O} }=\sum _{i}m_{i}{\overrightarrow {\mathrm {OM} _{i}}}\wedge {\vec {v}}_{i}=\sum _{i}m_{i}({\overrightarrow {\mathrm {OH} _{i}}}+{\vec {r}}_{i})\wedge {\vec {v}}_{i}$

Comme la vitesse d'un point $\mathrm {M} _{i}$ est donnée par ${\vec {v}}_{i}={\vec {\omega }}\wedge {\vec {r}}_{i}$ où ${\vec {\omega }}$ désigne le vecteur vitesse angulaire, la composante du moment cinétique colinéaire à l'axe de rotation $\Delta$ correspond au second terme de la somme, et se met sous la forme :

${\vec {L}}_{\Delta }=\sum _{i}m_{i}{\vec {r}}_{i}\wedge {\vec {v}}_{i}=\sum _{i}m_{i}{\vec {r}}_{i}\wedge \left({\vec {\omega }}\wedge {\vec {r}}_{i}\right)=\sum _{i}m_{i}\left(r_{i}^{2}{\vec {\omega }}-({\vec {\omega }}\cdot {\vec {r}}_{i}){\vec {r}}_{i}\right)=\left(\sum _{i}m_{i}r_{i}^{2}\right){\vec {\omega }}=J_{\Delta }{\vec {\omega }}$, puisque ${\vec {\omega }}\cdot {\vec {r}}_{i}=0$ par hypothèse.

Dans les deux cas apparaît la grandeur caractéristique, qui ne dépend que de la géométrie des masses du solide, appelé moment d'inertie $J_{\Delta }$ par rapport à l'axe $\Delta$ :

$J_{\Delta }=\sum _{i}m_{i}r_{i}^{2}$.

Remarque : il convient de souligner que ${\vec {L}}_{\mathrm {O} }$, moment cinétique du solide au point $\mathrm {O}$ de l'axe de rotation n'est pas en général colinéaire à cet axe. Ce n'est que si l'axe de rotation coïncide avec un axe principal d'inertie du solide, ce qui est notamment le cas s'il est un axe de symétrie matérielle du solide, c'est-à-dire à la fois axe de symétrie géométrique et que pour tout couple de points $(\mathrm {M} _{i},\mathrm {M} _{j})$ symétriques l'un par rapport à l'autre on ait $m_{i}=m_{j}$, que ${\vec {L}}_{\mathrm {O} }$ sera colinéaire à l'axe de rotation. En fait on aura alors ${\vec {L}}_{\mathrm {O} }=J_{\Delta }{\vec {\omega }}$ par rapport à tout point $\mathrm {O}$ situé sur l'axe de rotation.

Unités du moment d'inertie

Du fait de sa définition, le moment d'inertie a les dimensions d'une masse par le carré d'une longueur, soit M·L². Son unité dans le système international d'unités pourra donc naturellement être exprimée en kg⋅m², unité qui n'a pas de nom propre.

On peut cependant remarquer que dans la relation $E={\frac {1}{2}}J_{\Delta }.\omega ^{2}$, la vitesse de rotation ω n'est pas exprimée en s⁻¹, mais en rad⋅s⁻¹. Toutefois, étant défini comme le rapport de deux longueurs, le radian, considéré comme unité dérivée du système international depuis la 20^e conférence générale du BIPM, est sans dimension, par suite ceci n'affecte pas l'homogénéité de l'équation aux dimensions de l'expression de l'énergie cinétique. Il ne serait cependant pas faux que le moment d'inertie soit exprimé en kg⋅m²⋅rad⁻², mais ce choix est rarement retenu en pratique. Son seul intérêt serait de rappeler qu'il s'agit d'une unité spécifiquement rattachée au mouvement de rotation, comme pour toutes les unités où le radian apparaît.

Généralisation à un solide quelconque

Par extension dans un solide considéré comme un ensemble continu de points matériels $x$ affectés d'une masse volumique $\rho$, le moment d'inertie s'écrit :

$J_{\Delta }\ {\overset {\mathrm {def} }{=}}\ \iiint d(x,\Delta )^{2}\;\mathrm {d} m=\iiint d(x,\Delta )^{2}\,\rho \;\mathrm {d} V$

où

• $d(x,\Delta )$ est la distance entre le point $x$ et l'axe Δ ;
• $\mathrm {d} V$ est un volume élémentaire autour de $x$ ;
• $\mathrm {d} m$ est la masse de ce volume élémentaire ;
• $\rho$ est la masse volumique autour de $x$.

Cette définition peut également prendre une forme vectorielle :

$J_{\Delta }\ {\overset {\mathrm {def} }{=}}\ \iiint {\bigl \|}{\vec {e}}_{\Delta }\wedge {\overrightarrow {\mathrm {OM} }}{\bigr \|}^{2}\mathrm {d} m=\iiint {\bigl \|}{\vec {e}}_{\Delta }\wedge {\overrightarrow {\mathrm {OM} }}{\bigr \|}^{2}\,\rho \;\mathrm {d} V$

où

• $\mathrm {O}$ est un point sur l'axe $\Delta$ ;
• $\mathrm {M}$ est un point du solide considéré ;
• ${\vec {e}}_{\Delta }$ est un vecteur unitaire de l'axe $\Delta$ ;
• $\wedge$ désigne le produit vectoriel et $\|\,\|$ la norme euclidienne.

Généralisation à un système quelconque

En toute rigueur, la notion de moment d'inertie n'est définie que si la quantité $\sum _{i}m_{i}r_{i}^{2}$ peut être isolée des expressions de l'énergie cinétique ou du moment cinétique de rotation, c'est-à-dire que dans la mesure où la vitesse angulaire est la même pour les points du système à un instant donné. Ceci n'est valable strictement que dans le cas du modèle du solide indéformable.

Toutefois, la définition précédente peut s'étendre à un système déformable, dès lors qu'il ne présente pas de rotation différentielle, ou que l'on peut négliger l'effet de celle-ci, de façon qu'il soit possible de considérer que tous les points du système ont à un instant donné la même vitesse angulaire. Par exemple un système articulé, constitué de plusieurs solides reliés entre eux par des liaisons, en rotation autour d'un axe fixe, si les vitesses angulaires de rotation entre les diverses parties sont petites devant la vitesse angulaire de rotation « globale ». Toutefois, pour un système déformable, le moment d'inertie n'est plus constant dans le temps.

Moment d'inertie et moment cinétique

Dans la comparaison entre le mouvement de rotation et le mouvement de translation, le moment cinétique est l'homologue de la quantité de mouvement, et pour un système isolé est une grandeur conservative.

Comme indiqué précédemment, il découle de la définition du moment d'inertie que plus les masses constituant un solide sont réparties loin de l'axe de rotation, plus son moment d'inertie par rapport à cet axe est important. Comme la composante axiale du moment cinétique est égale au moment d'inertie multiplié par sa vitesse angulaire de rotation, et du fait de son caractère conservatif, si le moment d'inertie d'un système diminue du fait d'une variation de sa géométrie interne, sa vitesse angulaire de rotation doit augmenter (et inversement).

Ainsi, le patineur sur glace rapproche les bras de son corps lors d'une pirouette. Cela a pour effet de diminuer son moment d'inertie, ce qui, par conservation du moment cinétique, implique une plus grande vitesse de rotation.

De même, les enfants jouant à mettre en rotation un tourniquet en courant à ses côtés atteignent une vitesse de rotation limitée par leur course. Ils peuvent cependant ensuite sauter sur le tourniquet en mouvement, puis se positionner en son centre, augmentant ainsi significativement la vitesse de rotation initialement obtenue.

Cas général : tenseur d'inertie

La notion de moment d'inertie a été mise en évidence à partir du mouvement de rotation autour d'un axe fixe d'un solide. Or le mouvement général d'un solide par rapport à un référentiel (R) quelconque peut se décomposer en celui de son centre d'inertie C (affecté de la masse totale du système) et un mouvement de rotation propre autour de C dans le référentiel lié à ce point, en translation par rapport à (R), appelé référentiel barycentrique (noté (R^*))

Il est alors possible d'exprimer comme précédemment le moment cinétique et l'énergie cinétique propres du système, c'est-à-dire évalué dans (R^*), notés respectivement ${\vec {L}}^{*}$ et $E_{c}^{*}$, ce qui permet de mettre également en évidence une grandeur ne dépendant que de la géométrie des masses du solide, et généralisant la notion précédente, qui ne se réduit cependant plus à une grandeur scalaire, mais sera représenté par un tenseur, le tenseur d'inertie (appelé aussi opérateur ou matrice d'inertie).

Mise en évidence et définition du tenseur d'inertie

Le moment cinétique propre ${\vec {L}}^{*}$ d'un solide de vecteur rotation instantané ${\vec {\omega }}(t)$ s'écrit :

${\vec {L}}^{*}=\sum _{i}m_{i}{\vec {r}}_{i}\wedge {\vec {v}}_{i}^{*}=\sum _{i}m_{i}{\vec {r}}_{i}\wedge \left({\vec {\omega }}\wedge {\vec {r}}_{i}\right)=\sum _{i}m_{i}\left(r_{i}^{2}{\vec {\omega }}-\left({\vec {\omega }}\cdot {\vec {r}}_{i}\right){\vec {r}}_{i}\right)$,

Il est possible d'exprimer par exemple la composante suivant x en coordonnées cartésiennes, ce qui donne :

$L_{x}^{*}=\sum _{i}m_{i}r_{i}^{2}\omega _{x}-\sum _{i}m_{i}\left(x_{i}\omega _{x}+y_{i}\omega _{y}+z_{i}\omega _{z}\right)x_{i}=\left(\sum _{i}m_{i}\left(y_{i}^{2}+z_{i}^{2}\right)\right)\omega _{x}-\left(\sum _{i}m_{i}x_{i}y_{i}\right)\omega _{y}-\left(\sum _{i}m_{i}x_{i}z_{i}\right)\omega _{z}$.

Dans cette expression, les facteurs entre parenthèses représentent respectivement le moment d'inertie du solide par rapport à l'axe Ox, noté $I_{\mathrm {O} x}$, et deux termes homogènes à un moment d'inertie, appelés produits d'inertie, notés $I_{xy}$ et $I_{xz}$. Il est possible d'écrire $L_{x}$ sous la forme :

$L_{x}={\begin{bmatrix}I_{\mathrm {O} x}&-I_{xy}&-I_{xz}\end{bmatrix}}{\vec {\omega }}$,

et en procédant de même pour les autres composantes, il vient finalement l'expression du moment cinétique propre sous la forme :

${\vec {L}}^{*}={\bar {\bar {I}}}{\vec {\omega }}$,

avec ${\bar {\bar {I}}}$ tenseur (ou opérateur) d'inertie, qui est défini par :

${\bar {\bar {I}}}={\begin{bmatrix}I_{\mathrm {O} x}&-I_{xy}&-I_{xz}\\-I_{xy}&I_{\mathrm {O} y}&-I_{yz}\\-I_{xz}&-I_{yz}&I_{\mathrm {O} z}\end{bmatrix}}={\begin{bmatrix}\sum _{i}m_{i}(y_{i}^{2}+z_{i}^{2})&-\sum _{i}m_{i}x_{i}y_{i}&-\sum _{i}m_{i}x_{i}z_{i}\\-\sum _{i}m_{i}x_{i}y_{i}&\sum _{i}m_{i}(x_{i}^{2}+z_{i}^{2})&-\sum _{i}m_{i}y_{i}z_{i}\\-\sum _{i}m_{i}x_{i}z_{i}&-\sum _{i}m_{i}y_{i}z_{i}&\sum _{i}m_{i}(x_{i}^{2}+y_{i}^{2})\end{bmatrix}}$,

expression dans laquelle les éléments diagonaux sont les moments d'inertie du solide par rapport aux divers axes, et les éléments non diagonaux sont les produits d'inertie. De même, l'énergie cinétique propre s'écrit : $E_{c}^{*}={\tfrac {1}{2}}{\vec {\omega }}\cdot ({\bar {\bar {I}}}{\vec {\omega }})$.

Il est clair que ${\bar {\bar {I}}}$ ne dépend que de la géométrie des masses du solide. Il découle des relations précédentes que dans le cas général le moment cinétique propre du système n'est pas colinéaire à l'axe instantané de rotation, les relations précédentes généralisant celles obtenues dans le cas de la rotation autour d'un axe fixe.

Le moment d'inertie par rapport à un axe quelconque, de direction donnée par le vecteur unitaire ${\vec {n}}$ est alors donné par $I={\vec {n}}\cdot ({\bar {\bar {I}}}{\vec {n}})$.

Caractère tensoriel de ${\bar {\bar {I}}}$ - Axes principaux d'inertie

Expression générale du tenseur ${\bar {\bar {I}}}$

Il est facile de montrer que ${\bar {\bar {I}}}$ est effectivement un tenseur. En effet, en adoptant la notation $x_{i,1}=x_{i},x_{i,2}=y_{i},x_{i,3}=z_{i}$ il est possible de remarquer que la composante $I_{\alpha \beta }$ de ${\bar {\bar {I}}}$ se met sous la forme suivante :

$I_{\alpha \beta }=\left(\sum _{i}m_{i}r_{i}^{2}\right)\delta _{\alpha \beta }-\sum _{i}m_{i}x_{i,\alpha }x_{i,\beta }$.

Le premier terme est le produit d'un scalaire (le moment d'inertie par rapport au point O, $\sum _{i}m_{i}r_{i}^{2}$) par un tenseur (le tenseur de Kronecker $\delta _{ij}$). Le second terme correspond à une somme dans lequel chaque terme correspond au produit d'un scalaire (la masse m_i) par $T_{i,\alpha \beta }=x_{i,\alpha }x_{i,\beta }$. Or il s'agit là des composantes du tenseur résultant du produit tensoriel du vecteur ${\vec {r}}_{i}$ par lui-même, donc celles d'un tenseur. Par suite ${\bar {\bar {I}}}$ est donc bien un tenseur d'ordre deux : ceci était nécessaire de façon à assurer le caractère invariant par changement du système de coordonnées des expressions précédentes de ${\vec {L}}^{*}$ et $E_{c}^{*}$. Ce tenseur est évidemment symétrique.

Axes principaux d'inertie et éléments de symétrie matérielle du solide

Du fait de son caractère symétrique il est toujours possible de choisir un système d'axe tel que la matrice représentant ${\bar {\bar {I}}}$ soit diagonale. De tels axes sont dits axes principaux d'inertie. Les moments d'inertie correspondants sont appelés moments principaux d'inertie, et sont notés $I_{1},I_{2},I_{3}$. Leurs valeurs dépendent de la forme géométrique du solide et de la distribution de la masse en son sein, donc de l'expression que prend sa masse volumique $\rho ({\vec {r}})$ en chaque point du solide. Pour un solide homogène ρ est constante, et les moments principaux d'inertie ne dépendent alors que de la forme géométrique du solide.

En général un solide quelconque possède trois moments principaux d'inertie différents, il est appelé toupie asymétrique. Si deux moments principaux d'inertie sont égaux, par exemple $I_{1}=I_{2}\neq I_{3}$, le corps est qualifié de toupie symétrique, et si tous les moments principaux sont égaux, de toupie sphérique. Par exemple un parallélépipède homogène quelconque sera une toupie asymétrique, un cône ou un cylindre homogène, une toupie symétrique et une sphère homogène une toupie sphérique. La Terre, du fait de son aplatissement aux pôles, est également en général considérée comme une toupie symétrique.

La présence d'éléments de symétrie matérielle simplifie grandement la recherche des axes principaux d'inertie. En effet en présence de tels éléments certains produits d'inertie, par nature impairs par réflexions, s'annulent, ce qui permet de diagonaliser facilement la matrice représentant ${\bar {\bar {I}}}$.

Un élément de symétrie (point, axe, plan) matérielle est non seulement un élément par rapport auquel le solide est géométriquement symétrique, mais aussi pour lequel sa masse volumique présente la même symétrie. Ainsi un cylindre homogène comporte-t-il un axe de symétrie matérielle (son axe), par lequel passe une infinité de plans de symétrie matérielle, ainsi qu'un autre plan de symétrie matérielle qui est celui perpendiculaire à son axe passant par le milieu du cylindre. En revanche, si le cylindre est constitué de deux demi-cylindres tous deux homogènes, mais de masses volumiques différentes, accolés dans un plan contenant leurs axes, le cylindre possède toujours le plan de symétrie matérielle précédent, mais son axe n'est plus axe de symétrie matérielle. En revanche le plan de symétrie commun des deux demi-cylindres est toujours un plan de symétrie matérielle du système.

Il est possible de montrer, compte tenu des expressions précédentes des produits d'inertie, les propriétés suivantes :

• tout axe de symétrie matérielle est axe principal d'inertie ;
• tout axe perpendiculaire à un plan de symétrie matérielle est axe principal d'inertie.

Ainsi un cylindre homogène a pour axes principaux d'inertie son axe ainsi que tout axe qui y est perpendiculaire, passant par son centre. Par suite, deux des moments principaux d'inertie sont égaux et le tenseur d'inertie prend dans cette base la forme suivante :

${\bar {\bar {I}}}={\begin{bmatrix}I_{1}&0&0\\0&I_{1}&0\\0&0&I_{3}\end{bmatrix}}$. Il s'agit donc d'une toupie symétrique.

Relation avec le moment quadrupolaire d'une distribution de masses

Le potentiel de gravitation créé par une distribution quelconque de points matériels $M_{i}$ de masses $m_{i}$ ne se ramène pas en général à la forme obtenue dans le cas d'une distribution de masse à symétrie sphérique. Toutefois, la plupart des corps célestes (étoiles, planètes,…) possèdent approximativement cette symétrie et les écarts à la sphéricité demeurent faibles. Ces écarts sont évidemment liés à la répartition de la matière au sein de la distribution de masse, et donc doivent être (au moins pour les premières corrections) en relation avec le tenseur d'inertie de la distribution (assimilée à un solide parfait dans la suite) : de fait il est possible de montrer facilement que la première correction non nulle au potentiel sphérique fait intervenir une grandeur tensorielle, le tenseur de moment quadrupolaire de la distribution de masse, dont les composantes s'expriment de façon simple en fonction de celles du tenseur d'inertie.

De façon générale, à grande distance d'une distribution de masse le potentiel créé peut se mettre sous la forme d'un développement multipolaire : chaque point matériel constituant la distribution (de masse totale $\sum _{i}m_{i}=M$), repéré par le vecteur ${\vec {r}}_{i}={\overrightarrow {OM_{i}}}$ par rapport à une origine $O$ génère en ${\vec {r}}={\overrightarrow {OM}}$ le potentiel

$\phi _{i}({\vec {r}})=-{\frac {Gm_{i}}{r}}\sum _{k=0}^{\infty }\left({\frac {r_{i}}{r}}\right)^{k}P_{k}(\cos \theta _{i})$,

dans lequel $\theta _{i}$ est l'angle entre ${\vec {r}}$ et ${\vec {r}}_{i}$ et $P_{k}(z)$ est le polynôme de Legendre d'ordre $k$.

Comme $r>>r_{i}$, il est possible de se limiter aux 3 premiers termes, ce qui donne :

$\phi _{i}({\vec {r}})\approx -{\frac {Gm_{i}}{r}}\left(1+{\frac {r_{i}\cos \theta _{i}}{r}}+{\frac {r_{i}^{2}(3\cos ^{2}\theta _{i}-1)}{2r^{2}}}\right)$,

ce qui, compte tenu de $\cos \theta _{i}={\frac {{\vec {r}}_{i}\cdot {\vec {r}}}{rr_{i}}}$, devient :

$\phi _{i}({\vec {r}})\approx -{\frac {Gm_{i}}{r}}\left(1+{\frac {{\vec {r}}\cdot {\vec {r}}_{i}}{r^{2}}}+{\frac {(3({\vec {r}}\cdot {\vec {r}}_{i})^{2}-r^{2}r_{i}^{2})}{2r^{4}}}\right)$.

Le potentiel créé par la distribution en $M$ est égal à la somme des potentiels $\phi _{i}({\vec {r}})$, et chacun des trois termes se met alors sous la forme :

• terme d'ordre 0 ou polaire : $\Phi ^{(0)}({\vec {r}})=-{\frac {GM}{r}}$. C'est le potentiel à symétrie sphérique créé par une masse ponctuelle $M$ à l'origine $O$ ;
• terme d'ordre 1 ou dipolaire : $\Phi ^{(1)}({\vec {r}})=-{\frac {G{\vec {r}}}{r^{3}}}\cdot (\sum _{i}m_{i}{\vec {r}}_{i})$. Or, en prenant pour origine le centre de masse du système, $\sum _{i}m_{i}{\vec {r}}_{i}={\vec {0}}$, et le terme dipolaire est nul ;
• terme d'ordre 2 ou quadrupolaire : $\Phi ^{(2)}({\vec {r}})=-{\frac {G}{2r^{5}}}\sum _{i}m_{i}\left[3\left({\vec {r}}\cdot {\vec {r}}_{i}\right)^{2}-r^{2}r_{i}^{2}\right]$.

Ce dernier terme est donc la première correction, a priori non nulle, traduisant la non sphéricité à grande distance du potentiel de gravitation créé par une distribution de masse. La somme sur $i$ qui apparaît dans celui-ci peut être réécrite sous la forme :

$\sum _{i}m_{i}\left[3\left({\vec {r}}\cdot {\vec {r}}_{i}\right)^{2}-r^{2}r_{i}^{2}\right]={\vec {r}}\cdot ({\bar {\bar {Q}}}{\vec {r}})$,

où ${\bar {\bar {Q}}}$ est le moment quadrupolaire de la distribution de masse :

${\bar {\bar {Q}}}={\begin{bmatrix}3\sum _{i}m_{i}x_{i}^{2}-\sum _{i}m_{i}r_{i}^{2}&3\sum _{i}m_{i}x_{i}y_{i}&3\sum _{i}m_{i}x_{i}z_{i}\\3\sum _{i}m_{i}x_{i}y_{i}&3\sum _{i}m_{i}y_{i}^{2}-\sum _{i}m_{i}r_{i}^{2}&3\sum _{i}m_{i}y_{i}z_{i}\\3\sum _{i}m_{i}x_{i}z_{i}&3\sum _{i}m_{i}y_{i}z_{i}&3\sum _{i}m_{i}z_{i}^{2}-\sum _{i}m_{i}r_{i}^{2}\end{bmatrix}}$,

Chaque composante de ce tenseur se met donc sous la forme : $Q_{\alpha \beta }=3\sum _{i}m_{i}x_{i,\alpha }x_{i,\beta }-\delta _{\alpha \beta }\sum _{i}m_{i}r_{i}^{2}$, et compte tenu de l'expression générale des composantes du tenseur d'inertie d'un solide, $I_{\alpha \beta }=\left(\sum _{i}m_{i}r_{i}^{2}\right)\delta _{\alpha \beta }-\sum _{i}m_{i}x_{i,\alpha }x_{i,\beta }$, il vient finalement l'expression générale : $Q_{\alpha \beta }=\mathrm {Tr} ({\bar {\bar {I}}})\delta _{\alpha \beta }-3I_{\alpha \beta }$,

où $\mathrm {Tr} ({\bar {\bar {I}}})$ désigne la trace de la matrice représentant le tenseur d'inertie, c'est-à-dire la somme de ses termes diagonaux.

Les axes principaux d'inertie, pour lesquels la matrice représentant ${\bar {\bar {I}}}$ est diagonale, constituent également une base dans lequel ${\bar {\bar {Q}}}$ est diagonale, il vient pour les composantes de ce dernier l'expression $Q_{\alpha }=(I_{1}+I_{2}+I_{3})-3I_{\alpha }$.

Dans le cas où la distribution de masse est à symétrie sphérique, tous les moments principaux d'inertie sont égaux, et alors ${\bar {\bar {Q}}}=0$. Ce résultat est physiquement évident, et en fait dans ce cas tous les termes d'ordre supérieur du développement multipolaire précédent sont également nuls.

Une planète comme la Terre se comporte comme une toupie symétrique pour laquelle $I_{1}=I_{2}<I_{3}$, l'axe principal d'inertie selon $Oz$ correspondant pratiquement à son axe de rotation. Dans ce cas, d'après les formules précédentes, le potentiel gravitationnel à grande distance se met sous la forme :

$\Phi ({\vec {r}})\approx {\frac {-Gm}{r}}+{\frac {GM_{T}(I_{3}-I_{1})}{r^{3}}}\left[3\cos ^{2}\theta -1\right]$, où $\theta$ est l'angle entre la direction de ${\vec {r}}$ et celle de l'axe principal d'inertie $Oz$.

Ellipsoïde d'inertie

Le moment d'inertie d'un solide quelconque de tenseur d'inertie ${\bar {\bar {I}}}$ par rapport à un axe quelconque dont la direction est donnée par le vecteur unitaire ${\vec {n}}$ est donné par :

$I={\vec {n}}\cdot ({\bar {\bar {I}}}{\vec {n}})$,

en posant ${\vec {\rho }}={\frac {\vec {n}}{\sqrt {I}}}$ cette relation peut se mettre sous la forme :

${\vec {\rho }}\cdot ({\bar {\bar {I}}}{\vec {\rho }})=1$,

en explicitant avec les composantes $(\rho _{x},\rho _{y},\rho _{z})$ de ce vecteur, il vient l'équation d'un ellipsoïde :

$I_{\mathrm {O} x}\rho _{1}^{2}+I_{\mathrm {O} y}\rho _{2}^{2}+I_{\mathrm {O} z}\rho _{3}^{2}+2I_{xy}\rho _{1}\rho _{2}+2I_{xz}\rho _{1}\rho _{3}+2I_{yz}\rho _{2}\rho _{3}=1$,

laquelle prend dans les axes principaux d'inertie une forme particulièrement simple :

$I_{1}\rho _{1}^{2}+I_{2}\rho _{2}^{2}+I_{3}\rho _{3}^{2}=1$.

Cet ellipsoïde est appelé ellipsoïde d'inertie. Dans le cas d'une toupie symétrique, il s'agit d'un ellipsoïde de révolution, et dans le cas d'une toupie sphérique, d'une sphère. Cette notion n'a plus en général aujourd'hui qu'un intérêt historique, toutefois il est intéressant de remarquer que les axes principaux d'inertie sont les axes principaux de l'ellipsoïde d'inertie.

Moments d'inertie particuliers

Pour les exemples suivants, nous considérerons des solides homogènes ($\rho$ constant) et de masse $M$.

Boule

Pour une boule homogène de rayon $R$ et de centre $O$, les moments d'inertie au centre de la boule par rapport aux trois axes sont égaux. On peut donc écrire :

$J_{\Delta }={\frac {1}{3}}(J_{\mathrm {O} x}+J_{\mathrm {O} y}+J_{\mathrm {O} z})={\frac {1}{3}}\int _{V}[(y^{2}+z^{2})+(z^{2}+x^{2})+(x^{2}+y^{2})]\,\mathrm {d} m={\frac {2}{3}}\int _{V}(x^{2}+y^{2}+z^{2})\,\mathrm {d} m={\frac {2}{3}}\int _{V}r^{2}\,\mathrm {d} m$

En désignant par $\rho$ la masse volumique, $\mathrm {d} m=\rho \,4\pi r^{2}\mathrm {d} r$ donc :

$J_{\Delta }={\frac {2}{3}}\int _{0}^{R}r^{2}(\rho \,4\pi r^{2}\mathrm {d} r)={\frac {8}{15}}\pi R^{5}\rho$

Comme la masse de la boule est $M={\frac {4}{3}}\pi R^{3}\rho$, on obtient :

$J_{\Delta }={\frac {2}{5}}M\,R^{2}$.

Sphère (creuse)

Pour une sphère (par parenthèses une sphère est une surface, donc creuse), comme pour une boule, les moments d'inertie passant par son centre sont égaux. En l'occurrence, si son rayon est $R$, on a :

$J_{\Delta }={\frac {2}{3}}MR^{2}$

Barre

Dans le cas d'une barre de section négligeable et de longueur $L$, le moment d'inertie selon un axe perpendiculaire à la barre est, en son centre :

$J_{\Delta }=\int _{-L/2}^{+L/2}x^{2}\rho \,\mathrm {d} x={\frac {1}{12}}\rho L^{3}={\frac {1}{12}}ML^{2}$, (avec $M=\rho L$)

Ici, $\rho$ exprime une masse linéique (masse par unité de longueur).

Cela est vrai si l'axe de rotation passe par le milieu de la barre. La formule est différente si l'axe de rotation passe par son extrémité.

Carré

Dans le cas d'un carré de côté $a$, le moment d'inertie selon un axe perpendiculaire au plan du carré est, en son centre :

$J_{\Delta }=\int _{-a/2}^{+a/2}\int _{-a/2}^{+a/2}(x^{2}+y^{2})\rho \,\mathrm {d} x\mathrm {d} y={\frac {\rho }{6}}a^{4}={\frac {1}{6}}Ma^{2}$, (avec $M=\rho a^{2}$)

Ici, $\rho$ exprime une masse surfacique (masse par unité de surface).

Rectangle

Dans le cas d'un rectangle de grand côté $b$ et de petit côté $c$, le moment d'inertie selon un axe perpendiculaire au plan du rectangle (ici l'axe Oz) est, en son centre :

$J_{\Delta }=\int _{-b/2}^{+b/2}\int _{-c/2}^{+c/2}(x^{2}+y^{2})\rho (x,y)\,{\mathrm {d} }x{\mathrm {d} }y={\frac {\rho }{12}}bc^{3}+{\frac {\rho }{12}}cb^{3}={\frac {1}{12}}M(b^{2}+c^{2})$, (avec $M=\rho bc$)

Ici, $\rho$ exprime une masse surfacique (masse par unité de surface) pour une surface homogène, elle ne dépend donc pas de x et y. Remarquons que si $b=c$, on se ramène au cas du carré.

Parallélépipède

Dans le cas d'un parallélépipède de hauteur $h$, de grand côté $b$ et de petit côté $c$, le moment d'inertie selon l'axe le long de sa hauteur et en son centre est le même que celui pour un rectangle. En d'autres termes, la hauteur du parallélépipède ne joue aucun rôle :

$J_{\Delta }={\frac {1}{12}}M(b^{2}+c^{2})$

Cylindre plein

On utilisera les coordonnées cylindriques pour simplifier les calculs. Dans le cas d'un cylindre de rayon $R$ et de hauteur $h$, le moment d'inertie selon l'axe de révolution Oz du cylindre est :

${J_{\Delta }}_{z}=\int _{0}^{R}r^{2}\rho 2\pi rh\,\mathrm {d} r={\frac {1}{2}}\rho \pi R^{4}h={\frac {1}{2}}MR^{2}$, (avec $M=\rho \pi R^{2}h$)

Ici, $\rho$ exprime une masse volumique (masse par unité de volume).

On peut également déterminer le moment d'inertie selon tout axe Ox perpendiculaire à l'axe de révolution Oz du cylindre. Il vaut :

${J_{\Delta }}_{x}={\frac {MR^{2}}{4}}+{\frac {Mh^{2}}{12}}$

Cylindre creux

Dans le cas d'un cylindre creux de rayons intérieur $R_{1}$ et extérieur $R_{2}$, et de hauteur $h$, le moment d'inertie selon l'axe du cylindre est :

$J_{\Delta }=\int _{R_{1}}^{R_{2}}r^{2}(\rho 2\pi rh\,\mathrm {d} r)=2\rho \pi h\left[{\frac {r^{4}}{4}}\right]_{R_{1}}^{R_{2}}={\frac {1}{2}}\rho \pi h(R_{2}^{4}-R_{1}^{4})={\frac {1}{2}}M(R_{2}^{2}+R_{1}^{2})$, (où la masse du cylindre $M=\rho \pi h(R_{2}^{2}-R_{1}^{2})$)

Ici, $\rho$ exprime une masse volumique (masse par unité de volume).

Cône

Pour un cône (plein) dont la base a un rayon $R$, son moment d'inertie le long de sa hauteur est :

$J_{\Delta }={\frac {3}{10}}MR^{2}$, avec $M=\rho \,{\frac {\pi R^{2}h}{3}}$.

A noter que l'expression du moment d'inertie en fonction de sa masse et du rayon de sa base ne dépend pas de la hauteur du cône.

Tore

Pour un tore de paramètres $R$ et $r$, son moment d'inertie selon son axe de révolution est:

$J_{\Delta }=M\left(R^{2}+{\frac {3}{4}}r^{2}\right)$

Anneau fin

Pour un anneau fin (d'épaisseur négligeable) de rayon $R$ et de densité linéique (masse par unité de longueur) $\mu$, tous les éléments sont à la même distance de l'axe donc :

$J_{\Delta }=M\,R^{2}=(\mu \,2\pi R)R^{2}=2\pi R^{3}\mu$

Théorème de transport (ou théorème de Huygens-Steiner)

Soit l'axe $\Delta$ passant par le centre de masse d'un objet de masse $M$, et un axe $\Delta '$ parallèle à $\Delta$ et distant de $d$. En calculant comme précédemment le moment d'inertie, on retrouve la relation établie par Christian Huygens connue sous le nom de théorème de transport, théorème de Huygens, théorème de Steiner ou théorème des axes parallèles, qui donne le moment d'inertie $J_{\Delta '}$ en fonction de $J_{\Delta }$ :

$J_{\Delta '}=J_{\Delta }+M\,d^{2}$

À l'énergie cinétique de rotation propre d'un corps, s'ajoute celle de « translation » circulaire du centre de masse auquel on a affecté la masse totale du solide.

Une conséquence immédiate du théorème de Huygens est qu'il est moins coûteux (en énergie) de faire tourner un corps autour d'un axe passant par le centre de masse.

Voir aussi

Bibliographie

• Perez, Cours de physique : mécanique - 6^e édition, Masson, Paris, 2001, chapitre 17;
• Herbert Goldstein, Charles P. Poole Jr. et John L. Safko, Classical Mechanics [détail des éditions], chapitre 5;
• K. Gieck et R. Gieck (trad. G. Bendit), Formulaire technique, Gieck Verlag, 1997;
• Lev Landau et Evgueni Lifchits, Physique théorique, t. 1 : Mécanique [détail des éditions], chapitre 6.

Articles connexes

• Moment cinétique
• Moment d'inertie normalisé
• Moment quadratique

Liens externes

• Notices dans des dictionnaires ou encyclopédies généralistes :

  • Britannica
  • Den Store Danske Encyklopædi
  • Store norske leksikon
• Notices d'autorité :

  • LCCN
  • GND
  • Israël
  • Tchéquie