Mouvement à la Poinsot

En mécanique du solide, on appelle mouvement à la Poinsot, le mouvement d'un solide autour de son centre de gravité G, le moment des forces extérieures par rapport à G étant nul. Ce mouvement est caractérisé par la conservation du moment cinétique $\overrightarrow{L_G}$ et de l'énergie cinétique de rotation $E_c=\frac{1}{2}⟨\overrightarrow{L_G},\overrightarrow{\mathrm{\Omega }}⟩$, demi-produit scalaire du moment cinétique et du vecteur de rotation instantanée. Il existe 3 cas :

• le solide est à symétrie sphérique. Ses moments principaux d'inertie sont égaux : A = B = C. Alors le mouvement se réduit à une simple rotation uniforme d'axe le moment cinétique.
• le solide est à symétrie de révolution : A = B et C est différent. On parle de mouvement d'Euler-Poinsot de la toupie. Ce mouvement est équivalent à celui d'un cône roulant sans glissement sur un autre cône fixe. Le mouvement de la Terre en est un exemple.
• le solide est quelconque : C > B > A. Le mouvement est intégrable grâce aux fonctions elliptiques de Jacobi. Il se répète régulièrement après une certaine variation de la précession.

Solide à symétrie sphérique

Fichier:CubePoinsot.ogv Pour un tel solide, tous les moments principaux d'inertie sont égaux. Notons Modèle:Math leur valeur commune. Si $\overrightarrow{\mathrm{\Omega }}$ est le vecteur de rotation instantanée et $\overrightarrow{L_G}$ le moment cinétique du solide, on a $\overrightarrow{L_G}=I\overrightarrow{\mathrm{\Omega }}$. $\overrightarrow{L_G}$ étant constant, il en est de même de $\overrightarrow{\mathrm{\Omega }}$. Le mouvement est un simple mouvement de rotation uniforme d'axe le moment cinétique, à la vitesse angulaire $\frac{L_G}{I}$. Son énergie cinétique de rotation est $\frac{L_G^2}{2I}$.

Ce cas concerne la sphère, la boule, mais aussi le cube.

Solide à symétrie de révolution

Fichier:ToupiePoinsot.ogv Fichier:ToupiePoinsotCone.ogv

Soit un solide de révolution, d'axes principaux d'inertie $\overrightarrow{I}$, $\overrightarrow{J}$, $\overrightarrow{K}$, de moments d'inertie respectifs (A, A, C). Le référentiel (G, $\overrightarrow{I}$, $\overrightarrow{J}$, $\overrightarrow{K}$) sera dit référentiel lié au solide. Le référentiel d'origine G et dont les axes sont parallèles à un référentiel galiléen sera dit référentiel fixe. On rappelle que le moment cinétique $\overrightarrow{L_G}$ est constant. On peut le supposer colinéaire au vecteur $\overrightarrow{k}$ du référentiel fixe. On dispose des résultats suivants :

Théorème d'EulerModèle:,Modèle:, :

1. Les vecteurs $\overrightarrow{L_G}$, $\overrightarrow{K}$ et $\overrightarrow{\mathrm{\Omega }}$ sont coplanaires.
2. l'angle de nutation $\theta$ entre $\overrightarrow{L_G}$ et $\overrightarrow{K}$ est constant.
3. Le vecteur de rotation instantanée $\overrightarrow{\mathrm{\Omega }}$ a un module constant.
4. Le vecteur $\overrightarrow{\mathrm{\Omega }}$ décrit un cône d'axe $\overrightarrow{K}$ dans le référentiel lié au solide, dit cône du solide.
5. Le vecteur $\overrightarrow{\mathrm{\Omega }}$ décrit un cône d'axe $\overrightarrow{L_G}$ dans le référentiel fixe, dit cône de base.
6. Le cône du solide roule sans glissement le long du cône de base. À chaque instant, la génératrice commune aux deux cônes est dirigée par $\overrightarrow{\mathrm{\Omega }}$.
7. La vitesse angulaire de précession $\dot{\psi }$ est constante : $\dot{\psi }=\frac{L_G}{A}$.
8. La vitesse angulaire de rotation propre $\dot{\phi }$ est constante : $\dot{\phi }=L_G\cos (\theta )(\frac{1}{C}-\frac{1}{A})$.

En général, le mouvement n'est pas périodique dans le référentiel fixe car les vitesses angulaires de la précession et de la rotation propre ne sont pas commensurables.

Exemple de la Terre

En première approximation, la Terre a un moment cinétique constant et sa forme est essentiellement celle d'un solide à symétrie de révolution de forme ellipsoïdale, avec $\frac{A}{C-A}=305$. S'il s'agissait d'un ellipsoïde homogène, on aurait $\frac{A}{C-A}=\frac{a^2+c^2}{a^2-c^2}=297$, où a est le rayon équatorial et c le rayon d'un méridien. Mais la Terre est plus dense au centre. Le mouvement de précession libre calculé précédemment est de période 305 jours sidéraux. Le mouvement de rotation propre est d'un jour sidéral, soit Modèle:Unité. L'angle de nutation est très faible. Au pôle Nord, la polhodie est un cercle d'environ 10 m : le vecteur rotation et le moment cinétique sont presque alignés. Ce mouvement est actuellement suivi par l'IERS, qui observe des écarts entre la réalité expérimentale et le théorème d'Euler. Les raisons en sont multiples :

• En premier lieu, la Terre n'est pas un solide, mais c'est un corps doté d'une certaine élasticité, et dont la rigidité est voisine de celle de l'acier. Newcomb et Chandler reprirent la théorie d'Euler en tenant compte de ce fait, et trouva que le mouvement de précession d'Euler devait être changé en un mouvement de précession de Modèle:Unité.
• Par ailleurs, la Terre se comporte comme un matériau doué de viscosité ; le rebond glaciaire relève le bouclier canadien, entraînant une dérive de la polhodie.
• De plus, le mouvement des océans (courants et marées) et le mouvement de l'atmosphère perturbe la polhodie (voir rotation de la Terre).
• Il existe des tremblements de Terre (voir effet Sumatra), des mouvements de convection profonde correspondant à la tectonique des plaques, un noyau dont une partie est liquide et agité de mouvements de convection importants (voir géomagnétisme terrestre). À moment cinétique quasi constant, cela correspond à une rotation légèrement irrégulière, parfaitement perceptible compte tenu de la précision actuelle (10Modèle:Exp).

La modélisation précédente ne tient pas compte de la précession des équinoxes, décrite par Hipparque en 200 av. J.-C. : à long terme, le moment cinétique de la Terre n'est pas constant en raison de l'action de la Lune et du Soleil sur le bourrelet équatorial. Il subit donc par effet gyroscopique une précession, bien plus lente que la précession libre du mouvement à la Poinsot.

Solide quelconque

Ce cas, dit de la toupie asymétrique, concerne le mouvement d'un objet quelconque, sans symétrie particulière, et est bien plus compliqué. En effet, dans ce mouvement, les trois angles d'Euler (nutation, précession et rotation propre) varient tous les trois.

Mise en équation

On se place dans le référentiel lié au solide (G, $\overrightarrow{I}$, $\overrightarrow{J}$, $\overrightarrow{K}$), les axes étant des axes principaux d'inertie, les moments d'inertie étant respectivement A, B et C. On supposera par exemple que C > B > A. Notons $({\omega }_1,{\omega }_2,{\omega }_3)$ les composantes du vecteur de rotation instantanée $\overrightarrow{\mathrm{\Omega }}$. Les composantes du moment cinétique $\overrightarrow{L_G}$ sont $(A{\omega }_1,B{\omega }_2,C{\omega }_3)$. Selon les règles de dérivation dans une composition des mouvements, la conservation du moment cinétique dans le référentiel fixe s'exprime à partir de ses composantes dans le référentiel lié au solide sous la forme :

$0={\frac{\mathrm{d}\overrightarrow{L_G}}{\mathrm{d}t}}_{(\mathrm{f}\mathrm{i}\mathrm{x}\mathrm{e})}={\frac{\mathrm{d}\overrightarrow{L_G}}{\mathrm{d}t}}_{(\mathrm{s}\mathrm{o}\mathrm{l}\mathrm{i}\mathrm{d}\mathrm{e})}+\overrightarrow{\mathrm{\Omega }}\wedge \overrightarrow{L_G}$

ce qui donne les équations d'Euler suivantes :

$\{\begin{array}{c}A{\dot{\omega }}_1+(C-B){\omega }_2{\omega }_3=0\\ B{\dot{\omega }}_2+(A-C){\omega }_3{\omega }_1=0\\ C{\dot{\omega }}_3+(B-A){\omega }_1{\omega }_2=0\end{array}$

Résolution dans le référentiel lié au solide

En utilisant la constance du module du moment cinétique et de l'énergie cinétique de rotation, on a également :

$L_G^2=A^2{\omega }_1^2+B^2{\omega }_2^2+C^2{\omega }_3^2$

$2E_c=A{\omega }_1^2+B{\omega }_2^2+C{\omega }_3^2$

On détermine alors l'expression de $({\omega }_1,{\omega }_2,{\omega }_3)$ en fonction du temps Modèle:Math de la façon suivante. Les deux équations ci-dessus permettent d'exprimer ${\omega }_1$ et ${\omega }_3$ en fonction de ${\omega }_2$ et des constantes $L_G$ et $E_c$. On reporte ces expressions de ${\omega }_1$ et ${\omega }_3$ dans la deuxième équation d'Euler, permettant d'obtenir une équation différentielle non linéaire en ${\omega }_2$, à variables séparables. Celle-ci s'exprime simplement si on utilise les variables auxiliaires suivantes :

$\tau =t\sqrt{\frac{(C-B)(L_G^2-2E_{c}A)}{ABC}}$

$s={\omega }_2\sqrt{\frac{B(C-B)}{2E_{c}C-L_G^2}}$

$k^2=\frac{(B-A)(2E_{c}C-L_G^2)}{(C-B)(L_G^2-2E_{c}A)}$

On obtient alors $\frac{\mathrm{d}s}{\mathrm{d}\tau }=\sqrt{1-s^2}\sqrt{1-k^2{s}^2}$. En choisissant l'origine des temps lorsque ${\omega }_2$ prend une valeur nulle, cette équation différentielle s'intègre sous la forme :

$\tau =\underset{0}{\overset{s}{\int }}\frac{1}{\sqrt{1-x^2}\sqrt{1-k^2{x}^2}}\mathrm{d}x$

qui est l'intégrale elliptique de première espèce $F(s;k)$. On exprime Modèle:Math en fonction de τ, et donc ${\omega }_2$ en fonction du temps Modèle:Math en prenant la réciproque de cette fonction, ce qui conduit à $s=\mathrm{s}\mathrm{n}(\tau )$, où Modèle:Math est l'une des fonctions elliptiques de Jacobi (les autres étant Modèle:Math et Modèle:Math). On en déduit ensuite ${\omega }_1$ et ${\omega }_3$. Les expressions finales sont :

${\omega }_1=\sqrt{\frac{2E_{c}C-L_G^2}{A(C-A)}}\mathrm{c}\mathrm{n}(\tau )$

${\omega }_2=\sqrt{\frac{2E_{c}C-L_G^2}{B(C-B)}}\mathrm{s}\mathrm{n}(\tau )$

${\omega }_3=\sqrt{\frac{L_G^2-2E_{c}A}{C(C-A)}}\mathrm{d}\mathrm{n}(\tau )$

Ces fonctions sont périodiques de période $T=4K\sqrt{\frac{ABC}{(C-B)(L_G^2-2E_{c}A)}}$ avec $K=\underset{0}{\overset{1}{\int }}\frac{1}{\sqrt{1-x^2}\sqrt{1-k^2{x}^2}}$. Ainsi, dans le référentiel lié au solide, le vecteur de rotation instantanée a une variation périodique.

Résolution dans le référentiel fixe

On remonte ensuite à l'expression des angles d'Euler en fonction du temps, en prenant la direction $\overrightarrow{k}$ du référentiel fixe dans la direction du vecteur $\overrightarrow{L_G}$, puis en exprimant $\overrightarrow{k}$ dans le référentiel $(\overrightarrow{I},\overrightarrow{J},\overrightarrow{K})$ lié au solide en fonction des angles de nutation θ et de rotation propre φ. On obtient ainsi :

$\overrightarrow{L_G}=L_G\overrightarrow{k}=L_G(\cos (\theta )\overrightarrow{K}+\sin (\theta )\sin (\phi )\overrightarrow{I}+\sin (\theta )\cos (\phi )\overrightarrow{J})=A{\omega }_1\overrightarrow{I}+B{\omega }_2\overrightarrow{J}+C{\omega }_3\overrightarrow{K}$

On en déduit :

$\cos (\theta )=\frac{C{\omega }_3}{L_G}=\sqrt{\frac{C(L_G^2-2E_{c}A)}{L_G^2(C-A)}}\mathrm{d}\mathrm{n}(\tau )$

$\tan (\phi )=\frac{A{\omega }_1}{B{\omega }_2}=\sqrt{\frac{A(C-B)}{B(C-A)}}\frac{\mathrm{c}\mathrm{n}(\tau )}{\mathrm{s}\mathrm{n}(\tau )}$

Quant à la précession $\psi$, on l'obtient par l'expression des composantes du vecteur de rotation instantanée $\overrightarrow{\mathrm{\Omega }}$ selon les vecteurs $\overrightarrow{I}$ et $\overrightarrow{J}$ en fonction des angles d'Euler, et qui valent respectivement ${\omega }_1=\dot{\psi }\sin (\theta )\sin (\phi )+\dot{\theta }\cos (\phi )$ et ${\omega }_2=\dot{\psi }\sin (\theta )\cos (\phi )-\dot{\theta }\sin (\phi )$. On en déduit :

$\dot{\psi }=\frac{{\omega }_1\sin (\phi )+{\omega }_2\cos (\phi )}{\sin (\theta )}=L_G\frac{A{\omega }_1^2+B{\omega }_2^2}{A^2{\omega }_1^2+B^2{\omega }_2^2}$

$\psi$ est une primitive de la fonction ainsi trouvée. Au bout d'une période Modèle:Math donnée au paragraphe précédent, la précession a tourné d'un certain angle et le mouvement se répète relativement à cet angle, mais il n'y a pas de raison que cet angle soit commensurable à un tour complet, de sorte que le mouvement dans le référentiel fixe n'est pas en général périodique.

Polhodie et herpolhodie

Exemple d'ellipsoïdes (E1) en gris et (E2) en jaune. À un facteur constant près, le vecteur de rotation instantanée $\overrightarrow{\mathrm{\Omega }}$ décrit la polhodie, courbe d'intersection des deux ellipsoïdes.

Poinsot a décrit le mouvement géométriquement comme suit :

Dans le référentiel lié au solide, soit (E1) l'ellipsoïde d'inertie, d'équation $A{X}^2+B{Y}^2+C{Z}^2=1$, et soit (E2) l'ellipsoïde d'équation $A^2{X}^2+B^2{Y}^2+C^2{Z}^2=\frac{L_G^2}{2E_c}$.

En raison des deux équations $2E_c=A{\omega }_1^2+B{\omega }_2^2+C{\omega }_3^2$ (conservation de l'énergie cinétique de rotation) et $L_G^2=A^2{\omega }_1^2+B^2{\omega }_2^2+C^2{\omega }_3^2$ (conservation du module du moment cinétique), le vecteur $\frac{\overrightarrow{\mathrm{\Omega }}}{\sqrt{2E_c}}$ décrit la courbe d'intersection de ces deux ellipsoïdes. Cette courbe s'appelle la polhodie du mouvement.

Si l'on trace sur l'ellipsoïde d'inertie les diverses polhodies pour des valeurs différentes de Modèle:Math, on constate qu'elles se répartissent en quatre zones. Toujours avec la convention Modèle:Math > Modèle:Math > Modèle:Math, deux zones sont constituées des polhodies (en bleu sur la figure ci-contre) qui entourent l'axe Modèle:Math correspondant au moment cinétique Modèle:Math le plus élevé. Deux autres zones sont constituées des polhodies (en rouge sur la figure ci-contre) qui entourent l'axe Modèle:Math correspondant au moment cinétique Modèle:Math le plus faible. Ces polhodies sont fermées, et l'on retrouve ainsi le fait que le vecteur de rotation instantanée $\overrightarrow{\mathrm{\Omega }}$ se déplace périodiquement dans le référentiel lié au solide. Les quatre zones précédentes sont limitées par deux courbes séparatrices se coupant selon l'axe intermédiaire Modèle:Math. Cette configuration s'interprète mécaniquement par le fait qu'un mouvement de rotation proche des axes Modèle:Math et Modèle:Math sera stable (et plus stable selon Modèle:Math que selon Modèle:Math), alors qu'il sera instable autour de l'axe intermédiaire Modèle:Math. À titre d'exemple, il est possible de faire tourner une boîte d'allumettes autour de l'axe perpendiculaire aux deux grandes faces de la boîte, qui est le plus stable (ce que chacun sait également, quand il effectue des ricochets sur l'eau avec des pierres plates) ; il l'est également (mais avec plus de difficulté) autour de son grand axe. Mais il n'y a aucun espoir de lancer la boîte d'allumettes selon l'axe perpendiculaire aux grattoirs : c'est bien une position du mouvement stationnaire, mais elle est instable (effet Djanibekov).

Roulement sans glissement de l'ellipsoïde d'inertie d'un solide sur un plan (P) fixe. La polhodie est en rouge dans l'ellipsoïde, et l'herpolhodie en vert dans le plan (P). Le vecteur de rotation instantanée est colinéaire au vecteur noir qui relie le centre d'inertie du solide au point de contact entre l'ellipsoïde et le plan. Le moment cinétique est colinéaire au vecteur vertical noir.

Considérons maintenant le plan (P) tangent à l'ellipsoïde d'inertie au point $\frac{\overrightarrow{\mathrm{\Omega }}}{\sqrt{2E_c}}$. La normale à (P) est dirigée par le gradient de la fonction $A{X}^2+B{Y}^2+C{Z}^2$ en ce point, à savoir $\frac{1}{\sqrt{2E_c}}(2A{\omega }_1,2B{\omega }_2,2C{\omega }_3)$, qui, à un facteur près, n'est autre que le moment cinétique $\overrightarrow{L_G}$. L'équation du plan tangent est donc :

$\overrightarrow{u}\in \mathrm{(}\mathrm{P}\mathrm{)}⟺⟨\overrightarrow{L_G},\overrightarrow{u}⟩=⟨\overrightarrow{L_G},\frac{\overrightarrow{\mathrm{\Omega }}}{\sqrt{2E_c}}⟩$

Comme $⟨\overrightarrow{L_G},\overrightarrow{\mathrm{\Omega }}⟩=2E_c$, cette équation se simplifie en :

$⟨\overrightarrow{L_G},\overrightarrow{u}⟩=\sqrt{2E_c}$

$E_c$ et $\overrightarrow{L_G}$ étant constants, il s'agit d'un plan qui reste fixe pendant le mouvement. Le point de contact $\frac{\overrightarrow{\mathrm{\Omega }}}{\sqrt{2E_c}}$ entre ce plan et l'ellipsoïde d'inertie appartenant à l'axe de rotation instantanée, sa vitesse instantanée est nulle, ce qui signifie que l'ellipsoïde d'inertie roule et pivote sans glissement sur le plan fixe (P). La trace du point de contact dans le plan (P) s'appelle herpolhodie (on peut démontrer par ailleurs que cette courbe est sans point d'inflexion).

Voir aussi

Bibliographie

• Modèle:Ouvrage
• Modèle:Ouvrage
• Modèle:Ouvrage

Articles connexes

• Effet Djanibekov
• Gyroscope
• Mouvement de Lagrange de la toupie
• Pendule simple

Liens externes

• Modèle:Lien web, dont une autre version est : Modèle:Article