Aperçu Hero Physique
Article

Chargement de l'article...

... • Par Étude de Physique
Image de couverture
Problème de Physique Avancé : Mouvement d'une Toupie Symétrique
NIVEAU : Master / École d'Ingénieurs DOMAINE : Mécanique Classique THÈME : Dynamique du Solide (Effets Gyroscopiques)

Mouvement d'une Toupie Symétrique

Sonde ASTRA-V en Salle Blanche
1. Présentation du Système PhysiqueÉNONCÉ DU PROBLÈME
📝 Contexte de l'Étude : Anomalie en Salle Blanche

📋 Mémorandum Interne - Département d'Ingénierie Spatiale (CNES)

  • Destinataire : Pôle de Dynamique des Solides
  • Projet : Sonde Spatiale ASTRA-V
  • Statut : 🔴 URGENT - Phase de test au sol

Nous effectuons actuellement une série de tests décisifs en salle blanche sur le module de pointage principal de la sonde spatiale. Le système repose sur un rotor de stabilisation monté sur un pivot parfait.

Physiquement, ce système se comporte exactement comme une toupie symétrique pesante en rotation.

Le problème identifié :
Soumis au champ de pesanteur terrestre de notre laboratoire, ce dispositif subit un couple gravitationnel continu. Ce phénomène perturbe son alignement de référence, générant une oscillation dangereuse.

L'enjeu de votre mission (Exercice de niveau Ingénieur) :

Avant d'autoriser l'intégration définitive du module pour le vol spatial, vous devez impérativement modéliser son mouvement analytique. Il est vital d'anticiper l'évolution temporelle de deux phénomènes clés de la dynamique des solides :

  • 🔄 La précession (rotation de l'axe).
  • ↕️ La nutation (oscillation de l'axe).

Une modélisation précise nous évitera une perte de contrôle de l'orientation du satellite, tout en validant la résistance thermique du composant.

🔬
Objectif Analytique (Corrigé Pas-à-Pas)

Votre étude s'appuiera rigoureusement sur le formalisme de Lagrange. Votre mission se divise en trois axes majeurs :

  1. Modéliser la cinématique et la dynamique de rotation du rotor.
  2. Identifier les constantes du mouvement du système gyroscopique.
  3. Déterminer la condition mathématique stricte qui garantit une précession régulière (un mouvement fluide, sans aucune variation de nutation).
📐 SCHÉMA CINÉMATIQUE DU DISPOSITIF
X Y Z u z G x O ψ θ φ
Légende 1 : Repère fixe galiléen (salle d'essai) et repère mobile lié au rotor.
Légende 2 : Représentation des angles d'Euler \( (\psi, \theta, \phi) \) caractérisant la toupie.
📌
Hypothèse Simplificatrice : Le Modèle Idéal

Pour cet exercice d'application, nous adopterons un cadre de modélisation pur :

  • Liaison sphérique parfaite : Le point de pivot \( O \) à la base du rotor est un point fixe absolu, exempt de tout couple de frottement mécanique.
  • Symétrie idéale : La répartition des masses du rotor garantit une symétrie de révolution absolue autour de son axe principal, simplifiant le calcul de la matrice d'inertie.
2. Données Physiques, Constantes & Documentation

Pour garantir la plus stricte rigueur analytique attendue à ce niveau, vous devez extraire vos données d'entrée directement depuis la documentation officielle du constructeur.

📚 Cadre Théorique de Référence
  • Mécanique Analytique : Utilisation du Lagrangien et des coordonnées généralisées.
  • Cinématique du Solide : Repérage tridimensionnel via les angles d'Euler.
🗄️
RESSOURCE INGÉNIEUR : Extrait de la fiche technique (Série ASTRA)

Lisez attentivement ce tableau pour extraire les paramètres inertiels spécifiques au modèle étudié dans ce corrigé : la Sonde ASTRA-V (Modèle Type C).

Modèle de Rotor Masse totale \( m \) \(\text{kg}\) Bras de levier \( l \) \(\text{m}\) Inertie Axiale \( I_{\Delta} \) \( (\text{kg} \cdot \text{m}^2) \) Inertie Transv. \( I_{\text{T}} \) \( (\text{kg} \cdot \text{m}^2) \) Limite thermique \(\text{tr/min}\)
Sonde ASTRA-III (Type A) 1,8 0,10 \( 1,0 \cdot 10^{-2} \) \( 3,2 \cdot 10^{-2} \) 12 000
Sonde ASTRA-IV (Type B) 2,1 0,11 \( 1,2 \cdot 10^{-2} \) \( 3,6 \cdot 10^{-2} \) 10 500
Sonde ASTRA-V (Type C) 2,5 0,12 \( 1,5 \cdot 10^{-2} \) \( 5,0 \cdot 10^{-2} \) 8 500
⏱️ CONDITIONS INITIALES & CONSTANTES
  • 🚀
    Variables Cinématiques (Coordonnées Généralisées)
    • Angle de précession \( \psi(t) \) : Rotation globale autour de l'axe vertical absolu.
    • Angle de nutation \( \theta(t) \) : Inclinaison de l'axe de révolution du rotor.
    • Angle de rotation propre \( \phi(t) \) : Rotation ultra-rapide du rotor sur lui-même.
  • 🌍
    Environnement de Test
    • Champ de pesanteur : \( g = \mathbf{9,81} \text{ m} \cdot \text{s}^{-2} \) (Salle blanche)

3. Stratégie de Résolution & Travail Demandé

Afin de livrer un diagnostic infaillible aux ingénieurs, vous allez procéder à une modélisation analytique pas-à-pas. L'étude reposera sur l'expression de l'énergie du système pour en déduire les équations fondamentales du mouvement.

1

Question 1 : Construction du Modèle Énergétique

  • Exprimez l'énergie cinétique \( E_{\text{cinétique}} \) et l'énergie potentielle \( E_{\text{potentielle}} \) du rotor en utilisant exclusivement les angles d'Euler \( (\psi, \theta, \phi) \).
  • Déduisez-en l'expression complète et rigoureuse du Lagrangien \( \mathcal{L} \) du système gyroscopique.
2

Question 2 : Lois de Conservation et Symétries

  • Identifiez les coordonnées cycliques du système en observant attentivement la structure de votre Lagrangien.
  • Déduisez-en les expressions des moments conjugués associés et expliquez leur signification physique profonde.
  • Démontrez formellement pourquoi l'énergie mécanique totale \( E_{\text{mécanique}} \) de la sonde se conserve au cours du temps.
3

Question 3 : Condition de Précession Régulière

Les ingénieurs souhaitent stabiliser le rotor sur une trajectoire parfaite, avec un angle de nutation constant \( \theta = \theta_{0} \). Dans ce régime stationnaire, la vitesse de nutation \( \dot{\theta} \) et l'accélération \( \ddot{\theta} \) sont strictement nulles. L'axe du rotor décrit alors un cône parfait.

🎯 Visualisation de la trajectoire cible (Cône de Précession) dψ/dt Z z θ0 = 30° O
  • Appliquez l'équation de Lagrange relative à la coordonnée \( \theta \).
  • Démontrez que la vitesse de précession constante \( \dot{\psi} \) obéit obligatoirement à une équation polynomiale du second degré.
4

Question 4 : Limite de Stabilité Thermique et Gyroscopique

🚨 ALERTE DU MODÈLE : Le modèle mathématique de la précession régulière n'admet de solutions physiques que sous une condition critique. Si le rotor ralentit excessivement, la gravité l'emporte et le système s'effondre.

  • Déterminez la condition mathématique stricte pesant sur la vitesse de rotation axiale \( \omega_{z} \) (souvent assimilée à la rotation propre \( \dot{\phi} \) à haute vitesse) pour qu'une précession régulière soit physiquement possible.
  • Calculez numériquement cette vitesse critique minimale (en \( \text{tr/min} \)) pour la sonde ASTRA-V, en ciblant un angle d'inclinaison \( \theta_{0} = 30^{\circ} \).
  • Concluez en ingénieur : Comparez votre résultat avec la limite thermique tolérée par le constructeur (voir fiche technique) et validez ou rejetez la faisabilité du test au sol en tenant compte de la précision des chiffres significatifs.
CORRECTION

Mouvement d'une Toupie Symétrique

1
Construction du Modèle Énergétique
🎯 Objectif Mathématique & Physique

L'objectif de cette première étape est d'établir le bilan énergétique complet du rotor ASTRA-V dans le référentiel du laboratoire. La définition rigoureuse de ce modèle d'énergie (le Lagrangien) est le socle absolu sans lequel aucune équation de dynamique spatiale ne peut être calculée.

🧠 Analyse Préliminaire : La Méthode

Face à un solide en rotation 3D contraint par un point fixe (la liaison sphérique en \( O \)), l'approche newtonienne classique nécessiterait le calcul extrêmement complexe des forces de réaction de l'axe. En ingénierie spatiale, nous contournons cet obstacle par le raisonnement pas-à-pas suivant :

  • 📐 1. Identifier les degrés de liberté : Le système est libre de pivoter dans les 3 directions de l'espace.
  • 🧭 2. Définir le paramétrage : Nous utilisons les angles d'Euler pour repérer l'orientation.
  • 🔋 3. Calculer les énergies : Évaluer l'énergie liée au mouvement (cinétique) et à la gravité (potentielle).
  • ⚖️ 4. Fusionner : Soustraire ces deux énergies pour obtenir la fonction d'état de Lagrange.
Loi Fondamentale / Postulat : Formalisme Lagrangien et Cinématique

Pour un système mécanique conservatif (sans dissipation d'énergie par frottement), la dynamique est pilotée par la fonction de Lagrange \( \mathcal{L} \) :

\[ \mathcal{L} = E_{\text{cinétique}} - E_{\text{potentielle}} \]

Détail des grandeurs physiques :

  • 🔸 \( \mathcal{L} \) : Lagrangien du système (en \( \text{J} \)).
  • 🔸 \( E_{\text{cinétique}} \) : Énergie cinétique totale de rotation (en \( \text{J} \)).
  • 🔸 \( E_{\text{potentielle}} \) : Énergie potentielle de pesanteur (en \( \text{J} \)).

L'énergie cinétique d'un solide en pure rotation autour d'un point fixe s'obtient par le produit matriciel entre sa vitesse de rotation et son tenseur d'inertie :

\[ E_{\text{cinétique}} = \frac{1}{2} \vec{\Omega}^{T} \cdot \bar{\bar{I}} \cdot \vec{\Omega} \]

Détail des grandeurs physiques :

  • 🔸 \( \vec{\Omega} \) : Vecteur vitesse de rotation instantanée (en \( \text{rad} \cdot \text{s}^{-1} \)).
  • 🔸 \( \bar{\bar{I}} \) : Matrice d'inertie du solide évaluée au point de pivot (en \( \text{kg} \cdot \text{m}^{2} \)).
📐 Équations de Base

Pour rendre le calcul soluble, nous projetons le vecteur rotation \( \vec{\Omega} \) dans le repère mobile lié aux axes du rotor. Les composantes issues des angles d'Euler sont :

\[ \begin{aligned} \omega_{x} &= \dot{\psi} \sin\theta \sin\phi + \dot{\theta} \cos\phi \\ \omega_{y} &= \dot{\psi} \sin\theta \cos\phi - \dot{\theta} \sin\phi \\ \omega_{z} &= \dot{\psi} \cos\theta + \dot{\phi} \end{aligned} \]

L'énergie potentielle dépend, quant à elle, uniquement de l'altitude verticale \( z_{G} \) du centre de masse :

\[ E_{\text{potentielle}} = m g z_{G} \]
Traçabilité des Données / Conditions Initiales

Pour modéliser numériquement notre système, nous sourçons nos valeurs depuis la documentation technique officielle :

  • ⚖️ Masse du rotor : \( m = \mathbf{2,5} \text{ kg} \) (Source : Extrait ASTRA, Type C, Col. 2).
  • 📏 Bras de levier : \( l = \mathbf{0,12} \text{ m} \) (Source : Extrait ASTRA, Type C, Col. 3).
  • 🔄 Inertie Axiale : \( I_{\Delta} = \mathbf{0,015} \text{ kg} \cdot \text{m}^{2} \) (Source : Extrait ASTRA, Type C, Col. 4).
  • 🔀 Inertie Transversale : \( I_{\text{T}} = \mathbf{0,050} \text{ kg} \cdot \text{m}^{2} \) (Source : Extrait ASTRA, Type C, Col. 5).
  • 🌍 Champ de pesanteur : \( g = \mathbf{9,81} \text{ m} \cdot \text{s}^{-2} \) (Source : Constantes Locales de la salle blanche).

Astuce Mathématique / Méthodologique

Avant de manipuler l'énergie cinétique, développons le produit matriciel. Puisque le repère mobile est parfaitement aligné avec les axes de symétrie du rotor, la matrice d'inertie est strictement diagonale. Le produit complexe se simplifie immédiatement en une somme scalaire classique et facile à manipuler : 👉 \( \frac{1}{2} (I_{x}\omega_{x}^{2} + I_{y}\omega_{y}^{2} + I_{z}\omega_{z}^{2}) \).

Développement Analytique
1. Étape A : Simplification de l'énergie cinétique

💡 Raisonnement Pédagogique : La sonde ASTRA-V étant un rotor axisymétrique, ses inerties transversales sont rigoureusement égales (\( I_{x} = I_{y} = I_{\text{T}} \)) et son inertie axiale est isolée (\( I_{z} = I_{\Delta} \)).

L'énergie cinétique se réécrit donc en regroupant les termes transversaux :

\[ E_{\text{cinétique}} = \frac{1}{2} I_{\text{T}} \left( \omega_{x}^{2} + \omega_{y}^{2} \right) + \frac{1}{2} I_{\Delta} \omega_{z}^{2} \]

💡 Démonstration du terme transversal \( (\omega_{x}^{2} + \omega_{y}^{2}) \) : Élevons les composantes d'Euler au carré et additionnons-les pour exploiter l'identité trigonométrique fondamentale \( \cos^{2}\phi + \sin^{2}\phi = 1 \).

\[ \begin{aligned} \omega_{x}^{2} + \omega_{y}^{2} &= \left( \dot{\psi} \sin\theta \sin\phi + \dot{\theta} \cos\phi \right)^{2} + \left( \dot{\psi} \sin\theta \cos\phi - \dot{\theta} \sin\phi \right)^{2} \\ &= \left( \dot{\psi}^{2} \sin^{2}\theta \sin^{2}\phi + 2\dot{\psi}\dot{\theta}\sin\theta\sin\phi\cos\phi + \dot{\theta}^{2} \cos^{2}\phi \right) \\ &\quad + \left( \dot{\psi}^{2} \sin^{2}\theta \cos^{2}\phi - 2\dot{\psi}\dot{\theta}\sin\theta\cos\phi\sin\phi + \dot{\theta}^{2} \sin^{2}\phi \right) \end{aligned} \]

Les doubles produits s'annulent strictement. En factorisant, l'angle propre \( \phi \) disparaît :

\[ \begin{aligned} \omega_{x}^{2} + \omega_{y}^{2} &= \dot{\psi}^{2} \sin^{2}\theta \left( \sin^{2}\phi + \cos^{2}\phi \right) + \dot{\theta}^{2} \left( \cos^{2}\phi + \sin^{2}\phi \right) \\ &= \dot{\psi}^{2} \sin^{2}\theta \times (1) + \dot{\theta}^{2} \times (1) \\ &= \dot{\psi}^{2} \sin^{2}\theta + \dot{\theta}^{2} \end{aligned} \]

L'énergie cinétique totale du rotor prend ainsi sa forme définitive :

\[ E_{\text{cinétique}} = \frac{1}{2} I_{\text{T}} \left( \dot{\psi}^{2} \sin^{2}\theta + \dot{\theta}^{2} \right) + \frac{1}{2} I_{\Delta} \left( \dot{\psi} \cos\theta + \dot{\phi} \right)^{2} \]
2. Étape B : Démonstration de l'énergie potentielle

💡 Raisonnement Pédagogique : Dans le triangle rectangle formé par le pivot \( O \), le centre de masse \( G \) et l'axe vertical absolu, la fonction trigonométrique cosinus permet d'isoler l'altitude.

\[ \begin{aligned} z_{G} &= l \cos\theta \\ &\Rightarrow E_{\text{potentielle}} = m g l \cos\theta \end{aligned} \]
\[ \begin{aligned} \textbf{Bilan Étape 1 : } \mathcal{L} &= \frac{1}{2} I_{\text{T}} \left( \dot{\psi}^{2} \sin^{2}\theta + \dot{\theta}^{2} \right) + \frac{1}{2} I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right)^{2} - m g l \cos\theta \end{aligned} \]
Analyse Dimensionnelle / Sens Physique

Chaque bloc de ce Lagrangien s'exprime obligatoirement en Joules (\( \text{J} \)). Un regard averti révèle un fait de la plus haute importance : les variables de position \( \psi \) (précession) et \( \phi \) (rotation propre) n'apparaissent nulle part de manière explicite. Ce sont des coordonnées dites "cycliques", ce qui annonce de très puissantes propriétés de stabilité pour la suite.

Piège Classique

L'erreur de modélisation fatale consiste à confondre la vitesse de rotation propre (\( \dot{\phi} \)) avec la vitesse axiale totale (\( \omega_{z} \)). N'oubliez jamais le couplage cinématique \( \dot{\psi} \cos\theta \)! Il démontre un phénomène physique fascinant : faire pivoter le repère (précession) force mécaniquement l'axe de la toupie à tourner sur lui-même.

❓ Question Fréquente

Pourquoi s'acharne-t-on à projeter le vecteur rotation dans le repère mobile du rotor ?
Si nous travaillions dans le repère de la salle, la toupie s'inclinerait : la répartition globale de ses masses changerait à chaque milliseconde. La matrice d'inertie dépendrait du temps, rendant tout calcul insoluble. En "voyageant" sur le rotor, la géométrie est gelée, et la matrice d'inertie devient une constante fiable.

2
Lois de Conservation et Symétries
🎯 Objectif Mathématique & Physique

L'objectif de cette deuxième étape est d'exploiter les propriétés de symétrie repérées dans le Lagrangien pour extraire les constantes du mouvement. Ces invariants sont les "garde-fous" garantissant la stabilité orbitale inébranlable du gyroscope.

🧠 Analyse Préliminaire : La Puissance de Noether

La marche à suivre pour déverrouiller les lois de conservation est dictée par le célèbre Théorème de Noether :

  • 🔎 1. Observer le Lagrangien pour lister les variables spatiales "invisibles" (cycliques).
  • 🧮 2. Dériver \( \mathcal{L} \) par rapport à la vitesse de ces variables pour obtenir le moment conjugué.
  • 🔒 3. Conclure formellement sur la conservation éternelle de ce moment.
Loi Fondamentale / Postulat : Lois de Conservation

Le moment conjugué \( p_{q} \) associé à une coordonnée généralisée \( q \) se calcule par la dérivée partielle directe :

\[ p_{q} = \frac{\partial \mathcal{L}}{\partial \dot{q}} \]

L'équation d'Euler-Lagrange relie la variation temporelle de ce moment à la présence spatiale de la coordonnée dans l'énergie :

\[ \frac{d}{dt} \left( p_{q} \right) = \frac{\partial \mathcal{L}}{\partial q} \]

La règle d'or : Si \( \frac{\partial \mathcal{L}}{\partial q} = 0 \) (variable cyclique), le moment conjugué ne subit aucune variation. Il devient une constante absolue.

📐 Équations de Base

Pour la sonde ASTRA-V, nous avons identifié deux coordonnées cycliques fondamentales (\( \psi \) et \( \phi \)). Leurs équations de conservation s'écrivent :

\[ \frac{\partial \mathcal{L}}{\partial \psi} = 0 \Rightarrow p_{\psi} = \frac{\partial \mathcal{L}}{\partial \dot{\psi}} = \text{constante} \]
\[ \frac{\partial \mathcal{L}}{\partial \phi} = 0 \Rightarrow p_{\phi} = \frac{\partial \mathcal{L}}{\partial \dot{\phi}} = \text{constante} \]
Traçabilité des Données / Constantes
  • Variables Cycliques : Angles \( \psi \) et \( \phi \) (Source : Étape 1, observation visuelle de \( \mathcal{L} \)).
  • Variable Non-Cyclique : Angle de nutation \( \theta \) (Source : Étape 1, visible dans la gravité \( -mgl\cos\theta \)).

Astuce Mathématique / Méthodologique

La dérivation par rapport à la vitesse de précession \( \dot{\psi} \) demande une attention extrême. Cette vitesse contamine deux blocs du Lagrangien (transversal et axial). Il faut appliquer la règle d'or de la dérivation en chaîne : \( \frac{\partial}{\partial x} [u(x)^{2}] = 2 u(x) \cdot \frac{\partial u}{\partial x} \).

Développement Analytique
1. Étape A : Extraction du moment conjugué propre

💡 Raisonnement : On dérive partiellement le Lagrangien par rapport à \( \dot{\phi} \). Seul le bloc d'inertie axiale réagit à cette opération.

\[ \begin{aligned} p_{\phi} &= \frac{\partial}{\partial \dot{\phi}} \left[ \frac{1}{2} I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right)^{2} \right] \\ &= \frac{1}{2} I_{\Delta} \times \left[ 2 \left( \dot{\phi} + \dot{\psi} \cos\theta \right)^{1} \right] \times \frac{\partial}{\partial \dot{\phi}} \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \\ &= I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \times (1) \\ &\Rightarrow p_{\phi} = I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \end{aligned} \]
2. Étape B : Extraction du moment conjugué de précession

💡 Raisonnement : On distribue la dérivation sur les deux blocs d'énergie cinétique contenant \( \dot{\psi} \).

\[ \begin{aligned} p_{\psi} &= \frac{\partial}{\partial \dot{\psi}} \left[ \frac{1}{2} I_{\text{T}} \dot{\psi}^{2} \sin^{2}\theta \right] + \frac{\partial}{\partial \dot{\psi}} \left[ \frac{1}{2} I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right)^{2} \right] \\ &= \left[ \frac{1}{2} I_{\text{T}} \left( 2 \dot{\psi} \right) \sin^{2}\theta \right] + \left[ \frac{1}{2} I_{\Delta} \times 2 \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \times \frac{\partial}{\partial \dot{\psi}} \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \right] \\ &= I_{\text{T}} \dot{\psi} \sin^{2}\theta + I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \times (\cos\theta) \\ &\Rightarrow p_{\psi} = I_{\text{T}} \dot{\psi} \sin^{2}\theta + I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \cos\theta \end{aligned} \]
3. Étape C : Démonstration de l'invariance de l'énergie

💡 Raisonnement : Le pivot étant idéal (sans friction) et la salle blanche immobile, le Lagrangien ne dérive pas de manière explicite avec le temps.

\[ \frac{\partial \mathcal{L}}{\partial t} = 0 \iff \frac{d}{dt} \left( E_{\text{mécanique}} \right) = 0 \]
\[ \begin{aligned} \textbf{Bilan Étape 2 : } p_{\phi} &= \text{constante} \\ p_{\psi} &= \text{constante} \\ E_{\text{mécanique}} &= \text{constante} \end{aligned} \]
Analyse Dimensionnelle / Sens Physique

Les grandeurs \( p_{\phi} \) et \( p_{\psi} \) ont la dimension physique d'un moment cinétique (\( \text{kg} \cdot \text{m}^{2} \cdot \text{s}^{-1} \)).
🎯 Le paramètre \( p_{\phi} \) est la rigidité gyroscopique du rotor. C'est le bouclier qui empêche la sonde de basculer.
🎯 Le paramètre \( p_{\psi} \) garantit la conservation du moment cinétique projeté sur l'axe vertical du laboratoire (axe exempt de tout couple de contrainte externe).

Piège Classique

Ne commettez pas l'hérésie de déclarer que la nutation \( \theta \) est cyclique. L'angle \( \theta \) brille par sa présence dans l'énergie potentielle gravitationnelle. Le moment transversal \( p_{\theta} \) est donc perpétuellement agressé par la gravité, autorisant la toupie à plonger vers le sol.

❓ Question Fréquente

Pourquoi s'encombrer de \( p_{\phi} \) au lieu de dire que la vitesse \( \dot{\phi} \) est constante ?
C'est le génie de la dynamique des solides ! Ce n'est pas la vitesse \( \dot{\phi} \) seule qui est protégée par la nature, mais le bloc matriciel complet \( I_{\Delta}(\dot{\phi} + \dot{\psi} \cos\theta) \). Si la toupie penche (variation de \( \theta \)), sa vitesse de rotation interne subit une accélération compensatoire fulgurante pour maintenir \( p_{\phi} \) intouchable.

3
Condition de Précession Régulière
🎯 Objectif Mathématique & Physique

L'objectif est d'extraire la loi de pilotage de la précession régulière. Nous allons mathématiquement forcer l'axe du rotor à balayer un cône parfait de hauteur immuable, et déterminer l'équation d'équilibre des forces qui régit cette orbite.

🧠 Analyse Préliminaire : La Stabilité Forcée

Le régime stationnaire de précession régulière exige une approche systématique :

  • ⚓ 1. Figer l'espace : Imposer mathématiquement que l'inclinaison \( \theta \) est une constante temporelle absolue.
  • ⚖️ 2. Dresser le bilan : Écrire l'équation d'Euler-Lagrange spécifique à l'axe de nutation.
  • 🧹 3. Nettoyer l'équation : Injecter l'annulation des dérivées pour réduire le chaos différentiel en un élégant polynôme algébrique.
Loi Fondamentale / Postulat : Équation d'Euler-Lagrange

Le destin temporel de l'angle d'inclinaison \( \theta \) est gouverné par l'affrontement des couples (gravitationnels et inertiels) :

\[ \frac{d}{dt} \left( \frac{\partial \mathcal{L}}{\partial \dot{\theta}} \right) - \frac{\partial \mathcal{L}}{\partial \theta} = 0 \]

Détail des grandeurs physiques :

  • 🔸 \( \frac{\partial \mathcal{L}}{\partial \dot{\theta}} \) : Moment cinétique transversal \( p_{\theta} \) tendant à faire chuter la toupie.
  • 🔸 \( \frac{\partial \mathcal{L}}{\partial \theta} \) : Somme de toutes les forces généralisées agissant sur l'axe de basculement.
📐 Équations de Base

Pour modéliser l'exigence de la salle de contrôle, la position est verrouillée à sa valeur initiale :

\[ \theta = \theta_{0} \quad \text{(Constante absolue au cours du temps)} \]

Automatiquement, la vitesse de plongée est annihilée :

\[ \dot{\theta} = \frac{d\theta}{dt} = 0 \text{ rad} \cdot \text{s}^{-1} \]

Et l'accélération transversale disparaît du spectre :

\[ \ddot{\theta} = \frac{d^{2}\theta}{dt^{2}} = 0 \text{ rad} \cdot \text{s}^{-2} \]
Traçabilité des Données / Paramètres
  • 🎯 Angle cible de précession : \( \theta_{0} = \mathbf{30^{\circ}} \) (Source : Énoncé de la mission).
  • ⚙️ Vitesse axiale constante : L'inertie axiale dynamique étant constante, nous utiliserons le raccourci classique : \( \omega_{z} = \dot{\phi} + \dot{\psi}\cos\theta_{0} \).

🛑 Aération Pédagogique : L'écrasement de la dérivée

L'annulation radicale de \( \dot{\theta} \) est une bénédiction calculatoire. Le terme dynamique \( \frac{d}{dt} \left( \frac{\partial \mathcal{L}}{\partial \dot{\theta}} \right) \) correspond à la dérivée de \( I_{\text{T}} \dot{\theta} \), soit \( I_{\text{T}} \ddot{\theta} \). Puisque \( \ddot{\theta} = 0 \), toute cette complexité mathématique s'évapore littéralement de l'équation ! Il ne reste qu'à extraire la dérivée spatiale simple.

Développement Analytique
1. Étape A : Dérivation de la force généralisée sur \( \theta \)

💡 Raisonnement : Nous dérivons le Lagrangien complet par rapport à la position spatiale \( \theta \), avant de geler le système avec la condition stationnaire.

\[ \begin{aligned} \frac{\partial \mathcal{L}}{\partial \theta} &= \frac{\partial}{\partial \theta} \left[ \frac{1}{2} I_{\text{T}} \dot{\psi}^{2} \sin^{2}\theta \right] + \frac{\partial}{\partial \theta} \left[ \frac{1}{2} I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right)^{2} \right] - \frac{\partial}{\partial \theta} \left[ m g l \cos\theta \right] \end{aligned} \]

Exécution stricte des dérivations :

  • 🔹 Terme 1 : \( \frac{\partial}{\partial \theta} (\sin^{2}\theta) = 2\sin\theta\cos\theta \)
  • 🔹 Terme 2 : Dérivation en chaîne avec \( \frac{\partial}{\partial \theta} (\dot{\phi} + \dot{\psi} \cos\theta) = -\dot{\psi}\sin\theta \)
  • 🔹 Terme 3 : \( \frac{\partial}{\partial \theta} (\cos\theta) = -\sin\theta \)
\[ \begin{aligned} \frac{\partial \mathcal{L}}{\partial \theta} &= \frac{1}{2} I_{\text{T}} \dot{\psi}^{2} (2 \sin\theta \cos\theta) + \frac{1}{2} I_{\Delta} \times 2 \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \times (-\dot{\psi}\sin\theta) - m g l (-\sin\theta) \\ &= I_{\text{T}} \dot{\psi}^{2} \sin\theta \cos\theta - I_{\Delta} \left( \dot{\phi} + \dot{\psi} \cos\theta \right) \dot{\psi} \sin\theta + m g l \sin\theta \end{aligned} \]
2. Étape B : Formulation du Polynôme d'Équilibre Gyroscopique

💡 Raisonnement : L'équation de mouvement s'est écroulée en \( 0 - \frac{\partial \mathcal{L}}{\partial \theta} = 0 \). Figeons le système avec \( \theta = \theta_{0} \) et compactons la vitesse sous la forme \( \omega_{z} \).

\[ \begin{aligned} -\left( I_{\text{T}} \dot{\psi}^{2} \sin\theta_{0} \cos\theta_{0} - I_{\Delta} \omega_{z} \dot{\psi} \sin\theta_{0} + m g l \sin\theta_{0} \right) &= 0 \end{aligned} \]

Écartons l'éventualité triviale d'un rotor rigoureusement vertical (\( \sin\theta_{0} \neq 0 \)). En divisant l'égalité intégrale par \( -\sin\theta_{0} \), le polynôme directeur émerge :

\[ \begin{aligned} I_{\text{T}} \dot{\psi}^{2} \left( \frac{\sin\theta_{0} \cos\theta_{0}}{-\sin\theta_{0}} \right) - I_{\Delta} \omega_{z} \dot{\psi} \left( \frac{\sin\theta_{0}}{-\sin\theta_{0}} \right) + m g l \left( \frac{\sin\theta_{0}}{-\sin\theta_{0}} \right) &= 0 \\ - I_{\text{T}} \cos\theta_{0} \dot{\psi}^{2} + I_{\Delta} \omega_{z} \dot{\psi} - m g l &= 0 \end{aligned} \]

Multiplions par \( -1 \) pour obtenir la forme canonique noble du second degré, centrée sur l'inconnue spatiale \( \dot{\psi} \) :

\[ \begin{aligned} \left( I_{\text{T}} \cos\theta_{0} \right) \dot{\psi}^{2} - \left( I_{\Delta} \omega_{z} \right) \dot{\psi} + m g l &= 0 \end{aligned} \]
\[ \textbf{Bilan Étape 3 : } A \dot{\psi}^{2} + B \dot{\psi} + C = 0 \quad \text{avec} \quad A = I_{\text{T}} \cos\theta_{0} \]
Analyse Dimensionnelle / Sens Physique

Cette parabole dresse le théâtre tragique de la lutte des couples mécaniques.
🔴 Le terme quadratique \( A\dot{\psi}^{2} \) symbolise la force centrifuge aveugle qui cherche à fracasser l'axe à l'horizontale.
🟢 Le terme linéaire \( B\dot{\psi} \) (généré par le couplage de Coriolis) incarne le couple gyroscopique de la toupie.
⚫ Le terme absolu \( C = mgl \) est le châtiment gravitationnel qui tire la sonde vers le crash.
Si ce triumvirat de forces s'annule strictement, l'astre mécanique est stabilisé.

Piège Classique

Résistez au chant des sirènes vous incitant à redévelopper la vitesse axiale \( \omega_{z} \) (en réécrivant \( \dot{\phi} + \dot{\psi}\cos\theta_{0} \)) en plein milieu de l'arène algébrique ! La grandeur \( I_{\Delta}\omega_{z} \) a été couronnée constante inaliénable à la Question 2. La sceller sous le label \( B \) préserve la pureté de l'équation du second degré et vous sauve d'un désastre calculatoire.

❓ Question Fréquente

Pourquoi une même toupie offrirait-elle deux vitesses de précession possibles (les deux racines du polynôme) ?
C'est la signature de la physique non-linéaire. Sous une même impulsion moteur, le cosmos autorise deux chorégraphies stables :
1. Une précession lente (le mode commun, sous le joug de la gravité).
2. Une précession rapide (le mode furtif de l'ingénierie avancée, où l'inertie centrifuge dicte sa loi).

4
Limite de Stabilité Thermique et Gyroscopique
🎯 Objectif Mathématique & Physique

L'heure de vérité. Nous devons connecter la pureté de la théorie à la brutalité de l'ingénierie. Nous allons chiffrer le seuil d'effondrement gravitationnel de la sonde ASTRA-V, et statuer définitivement sur la survie du composant face aux contraintes thermiques du tir.

🧠 Analyse Préliminaire : La Survie du Modèle

La précession régulière a exigé une équation du type \( A \dot{\psi}^{2} + B \dot{\psi} + C = 0 \). Mais la réalité matérielle impose une sanction cruelle :

  • 🔍 1. Condition d'existence : La vitesse de précession \( \dot{\psi} \) doit obligatoirement être un nombre réel (\( \mathbb{R} \)).
  • 📉 2. Couperet du discriminant : Une parabole ne crache des racines réelles que si son discriminant \( \Delta_{\text{math}} \) est positif. S'il plonge dans le négatif, l'axe s'effondre dans le chaos de la nutation.
  • ⚙️ 3. Isoler le paramètre moteur : Poser l'inéquation du discriminant nous livrera sur un plateau d'argent le seuil de basculement de la vitesse moteur \( \omega_{z} \).
Loi Fondamentale / Postulat : La Barrière de Stabilité

Le couperet mathématique de survie s'écrit formellement :

\[ \Delta_{\text{math}} = B^{2} - 4 A C \ge 0 \]

Détail du code génétique de la parabole :

  • 🔸 \( A = I_{\text{T}} \cos\theta_{0} \) (Directeur inertiel).
  • 🔸 \( B = -I_{\Delta} \omega_{z} \) (Bouclier gyroscopique).
  • 🔸 \( C = m g l \) (Faille gravitationnelle).
📐 Évaluation de la Perturbation

Le point de non-retour est franchi lorsque le couple gyroscopique dépose les armes face à l'alliance de la gravité et de la centrifugation. L'inéquation de survie s'énonce ainsi :

\[ \left( - I_{\Delta} \omega_{z} \right)^{2} - 4 \left( I_{\text{T}} \cos\theta_{0} \right) \left( m g l \right) \ge 0 \]
Traçabilité des Données / Certification Numérique

Invocation intégrale des paramètres du banc d'essai pour le scellement du vol ASTRA-V :

  • ⚖️ Masse : \( m = \mathbf{2,5} \text{ kg} \) (Source : Extrait Constructeur, Ligne Type C, Col. 2).
  • 📏 Levier : \( l = \mathbf{0,12} \text{ m} \) (Source : Extrait Constructeur, Ligne Type C, Col. 3).
  • 🔄 Inertie Axiale : \( I_{\Delta} = \mathbf{0,015} \text{ kg} \cdot \text{m}^{2} \) (Source : Extrait Constructeur, Ligne Type C, Col. 4).
  • 🔀 Inertie Transversale : \( I_{\text{T}} = \mathbf{0,050} \text{ kg} \cdot \text{m}^{2} \) (Source : Extrait Constructeur, Ligne Type C, Col. 5).
  • 🎯 Angle cible : \( \theta_{0} = \mathbf{30^{\circ}} \) (Source : Cahier des charges de la mission).
  • 🔥 Plafond thermique : 8 500 tr/min (Source : Alerte d'incandescence, Tableau Constructeur, Col. 6).
  • 🌍 Pesanteur : \( g = \mathbf{9,81} \text{ m} \cdot \text{s}^{-2} \) (Source : Constantes du Labo).

Astuce Mathématique / La Méthode de Transparence

L'opération se fera en deux temps chirurgicaux, en gardant à l'esprit la règle d'or des chiffres significatifs (les données d'entrée en comptent deux, notre résultat final devra donc comporter 2 chiffres significatifs) :
1. Extraction de la variable vitale : Isoler le régime moteur \( \omega_{z} \).
2. La loi de la conversion pure : Transformer le résultat du Système International (\( \text{rad} \cdot \text{s}^{-1} \)) en un langage industriel lisible par les ingénieurs motoristes (\( \text{tr/min} \)), en déroulant étape par étape l'origine des facteurs multiplicatifs (interdiction absolue de toute conversion "magique").

Analyse de l'Effet (Calcul de la Limite Critique)
1. Étape A : Isolation algébrique de la vitesse critique

💡 Raisonnement : Évacuons le bloc de soustraction vers la droite de l'inégalité. Les variables physiques massiques et d'inertie garantissent la stricte positivité du bloc de droite (le cosinus d'un angle aigu tel que 30° étant lui-même positif). L'extraction de la racine carrée est donc mathématiquement sanctifiée.

\[ \begin{aligned} \left( I_{\Delta} \omega_{z} \right)^{2} &\ge 4 I_{\text{T}} m g l \cos\theta_{0} \\ \omega_{z}^{2} &\ge \frac{4 I_{\text{T}} m g l \cos\theta_{0}}{I_{\Delta}^{2}} \\ &\Rightarrow \omega_{z} \ge \frac{2 \sqrt{I_{\text{T}} m g l \cos\theta_{0}}}{I_{\Delta}} \end{aligned} \]
2. Étape B : Application Numérique (A.N.) en Système International

💡 Raisonnement : Blindez le calcul du bloc interne à la racine avant de diviser, pour annuler tout risque d'erreur de frappe en chaîne sur la calculatrice.

\[ \begin{aligned} \omega_{z,\text{minimale}} &= \frac{2 \sqrt{0,050 \times 2,5 \times 9,81 \times 0,12 \times \cos(30^{\circ})}}{0,015} \end{aligned} \]

Résolution du bloc radiculaire (avec \( \cos(30^{\circ}) \approx 0,8660 \)) :

\[ \begin{aligned} \omega_{z,\text{minimale}} &\approx \frac{2 \sqrt{0,050 \times 2,943 \times 0,8660}}{0,015} \\ &\approx \frac{2 \sqrt{0,1274}}{0,015} \\ &\approx \frac{2 \times 0,3570}{0,015} \\ &\approx 47,59... \text{ rad} \cdot \text{s}^{-1} \\ &\Rightarrow \omega_{z,\text{minimale}} \approx \mathbf{48} \text{ rad} \cdot \text{s}^{-1} \quad \text{*(arrondi à 2 c.s.)*} \end{aligned} \]
3. Étape C : Transparence absolue de la conversion industrielle

💡 Raisonnement : Repartons de la valeur non arrondie (\( 47,59 \)) pour éviter de propager l'erreur. La norme constructeur est en tours par minute (\text{tr/min}). Posons les fondations temporelles et spatiales de la conversion :
🌍 1 révolution circulaire (tour) vaut un angle cinématique de \( 2\pi \) radians. Équivalence spatiale : \( 1 \text{ rad} = \frac{1}{2\pi} \text{ tour} \).
⏱️ 1 minute rassemble 60 secondes. Équivalence temporelle : \( 1 \text{ s} = \frac{1}{60} \text{ minute} \).

Assimilons ces deux contraintes de transfert d'échelle sur notre valeur radiale brute :

\[ \begin{aligned} V_{\text{critique}} &= 47,59... \text{ rad} \cdot \text{s}^{-1} \\ &= 47,59... \times \left( \frac{1 \text{ tour}}{2\pi} \right) \times \left( \frac{60}{1 \text{ minute}} \right) \\ &= 47,59... \times \left( \frac{60}{2\pi} \right) \text{ tr} \cdot \text{min}^{-1} \\ &\approx 47,59... \times 9,549 \\ &\approx 454,4... \text{ tr/min} \end{aligned} \]

En respectant rigoureusement la contrainte des 2 chiffres significatifs, nous fixons la limite industrielle.

\[ \begin{aligned} V_{\text{critique}} &\Rightarrow \mathbf{4,5 \cdot 10^{2}} \text{ tr/min} \quad \text{*(soit } 450 \text{ tr/min)} \end{aligned} \]
\[ \textbf{Bilan Étape 4 : } \omega_{z,\text{critique}} \ge \mathbf{4,5 \cdot 10^{2}} \text{ tr/min} \]
Analyse Dimensionnelle / Sens Physique

L'inéquation \( \omega_{z} \ge 450 \text{ tr/min} \) démontre brillamment que la stabilité inertielle est un barrage d'énergie pure. Plus l'engin est lourd et perché haut sur son mât (puissance de destruction du terme \( mgl \)), plus le "bouclier gyroscopique" requis pour lutter contre la gravité est colossal. La valeur extraite (450 tr/min) rassure : grâce à l'architecture optimisée de ses matrices d'inertie (\( I_{\text{T}} \) et \( I_{\Delta} \)), l'ASTRA-V parvient à s'auto-stabiliser à un régime moteur d'une extrême douceur pour l'industrie spatiale.

Piège Classique d'Application Numérique

Le fléau de l'ingénieur reste l'intégration des constantes. Injecter aveuglément l'inertie axiale de la table \( 1,5 \cdot 10^{-2} \) sous la forme \( 1,5 \) ou \( 0,15 \) en plein cœur de la division finale pulvérisera votre diagnostic. Veillez à la frapper sous la forme chirurgicale \( 0,015 \).
La faute suprême résiderait dans le traitement de l'angle \( \theta_{0} = 30 \) : le laisser s'évaluer en radians dans la mémoire de votre puce de calcul ferait imploser le bloc de la racine carrée.

❓ Question Fréquente

Diagnostic Final : Conformément à la requête du CNES, le test au sol de la sonde ASTRA-V peut-il se poursuivre en toute sécurité thermique ?
Absolument. La modélisation mécanique exige que le rotor tourne à une vitesse supérieure à environ 450 tr/min pour maintenir son inclinaison sans tomber. La documentation limite la température du moteur à un seuil correspondant à 8 500 tr/min. La vitesse minimale requise représente à peine 5 % de la capacité maximale autorisée. L'équipe d'ingénierie dispose donc d'une marge colossale pour réaliser ce vol de test dans un confort matériel total. La théorie confirme l'expérience spatiale : GO FOR TEST.

Bilan Ingénieur : Analyse de Stabilité Orbitale
SEUIL CRITIQUE (Δ = 0) 4,5 × 10² tr/min LIMITE THERMIQUE 8 500 tr/min dψ/dt rad/s ωz tr/min PRÉCESSION RAPIDE Dominance Centrifuge PRÉCESSION LENTE Dominance Gravitationnelle ZONE D'EFFONDREMENT
✓ VALIDATION : La plage nominale de test se déploie confortablement dans le vide inertiel stable, à des milliers de rotations/minute de la rupture thermique.
Mouvement d'une Toupie Symétrique