Equivalent Mechanical Models for Sloshing

Résumé Technique : Modèles Mécaniques Équivalents pour le Ballonnement (Sloshing)

Énoncé du Problème
Le ballonnement des ergols demeure un défi critique pour la stabilité et le contrôle des engins spatiaux, particulièrement dans deux régimes distincts : les environnements de microgravité nécessitant une haute stabilité de pointage (par exemple, les missions interplanétaires) et les conditions de forte accélération ($g$) rencontrées lors du lancement, de la rentrée et de l'atterrissage. Bien que les modèles de pendule et de masse-ressort-amortisseur soient des outils standards pour les analystes de commande afin d'évaluer les perturbations induites par le ballonnement, une formulation mathématique rigoureuse reliant ces représentations, tenant compte spécifiquement des effets des forces longitudinales nominales dans les régimes à forte accéléation (high-$g$), est nécessaire pour une conception robuste de la guidance, de la navigation et du contrôle (GNC).

Méthodologie
L'article présente une dérivation rigoureuse de modèles mécaniques équivalents pour le ballonnement, progressant des dynamiques non linéaires vers des représentations linéarisées adaptées à l'analyse de commande.

1. Formulation Non Linéaire :

   • L'étude commence par un pendule sphérique unique de masse $m_P$ et de longueur $\ell$ attaché à un corps rigide $B$ de masse $m_B$.
   • En utilisant une approche Newtonienne, les auteurs dérivent les équations du mouvement couplées pour le système à 8 degrés de liberté (translation/rotation du corps rigide + angles du pendule).
   • La dérivation traite explicitement la force de tension $N$ et l'élimine pour réduire les variables du système, résultant en une forme matricielle compacte $\mathcal{M}\ddot{\mathbf{x}} = \mathbf{u}$.
   • La formulation est étendue pour accommoder des pendules multiples, représentant plusieurs réservoirs ou modes de ballonnement, en utilisant une matrice de masse à structure par blocs pour exploiter la parcimonie en vue d'une solution numérique efficace.
   • Un modèle d'amortissement visqueux est incorporé en ajoutant des termes proportionnels aux vitesses angulaires, qui affectent la dynamique du pendule sans modifier directement les équations du corps rigide.

2. Linéarisation et Analyse Modale :

   • Reconnaissant que les systèmes de commande embarqués estiment généralement le centre de masse en supposant une distribution d'ergols rigide, les auteurs reformulent la dynamique par rapport à un centre de masse global nominal ($\bar{G}$) plutôt qu'au centre de masse du corps rigide.
   • Le système est linéarisé autour d'un équilibre caractérisé par une attitude constante alignée avec le référentiel inertiel et une force longitudinale nominale prescrite $\bar{F}_z$.
   • Deux cas sont analysés : des forces définies dans le référentiel inertiel et des forces définies dans le référentiel lié au corps.
   • À travers des transformations de coordonnées et des manipulations algébriques, la dynamique linéarisée est transposée dans une « forme modale ». Cela implique de définir de nouvelles coordonnées modales ($\eta_\theta, \eta_\phi$) qui découplent la dynamique du pendule en deux oscillateurs plans indépendants.
   • Le système résultant est exprimé sous la forme d'un système MIMO de second ordre ($M\delta\ddot{\mathbf{x}} + K\delta\mathbf{x} = \delta\mathbf{u}_{ext}$), identifiant explicitement la matrice de participation modale et les termes de raideur induits par la force nominale.

3. Équivalence avec les Modèles Masse-Ressort-Amortisseur :

   • L'article dérive les équations linéarisées pour un système masse-ressort-amortisseur où la masse de ballonnement est contrainte de se déplacer latéralement via des ressorts de raideur $k$.
   • En fixant la raideur du ressort $k = m_P \omega_0^2$ (où $\omega_0$ est la fréquence naturelle dérivée du modèle de pendule), les auteurs démontrent que la dynamique linéarisée masse-ressort-amortisseur est mathématiquement identique à la dynamique linéarisée du pendule.
   • L'analyse souligne que la présence d'une force longitudinale nominale $\bar{F}_z$ est essentielle pour que le pendule présente un comportement oscillatoire analogue à un ressort ; en l'absence de cette force ($\bar{F}_z=0$), le modèle de pendule nécessite une modification (par exemple, l'ajout d'un troisième degré de liberté modal) pour représenter le ballonnement en 0-g.

Résultats Clés

• Équivalence Mathématique : L'étude établit formellement l'équivalence exacte entre le modèle de pendule linéarisé et le modèle classique masse-ressort-amortisseur sous des conditions de forte accélération, clarifiant comment la force nominale agit comme le mécanisme de rappel dans la formulation du pendule.
• Représentation Modale : La dérivation produit une représentation modale claire où les modes de ballonnement sont découplés de la dynamique du corps rigide dans la matrice de raideur, le couplage n'apparaissant que dans la matrice de masse via la matrice de participation modale.
• Validation :
  • Validation Non Linéaire : Le modèle non linéaire a été validé par rapport à une implémentation MATLAB Simscape Multibody (moyeu rigide avec quatre pendules). Les simulations temporelles ont montré un excellent accord, avec des erreurs maximales inférieures à 0,002 % de l'amplitude de la réponse, les écarts étant attribués uniquement aux différences entre les schémas d'intégration numérique.
  • Validation Linéaire : Les diagrammes de Bode comparant le modèle linéarisé analytique avec le modèle Simscape linéarisé ont démontré un excellent accord sur toute la plage de fréquences pour les entrées de force dans les référentiels inertiel et lié au corps. De légers écarts à basse fréquence dans le cas du référentiel inertiel ont été attribués aux limites de précision numérique du processus de linéarisation de Simscape.

Signification et Revendications
L'article prétend fournir une base mathématique rigoureuse pour la modélisation du ballonnement des ergols dans les environnements à forte accélération. Sa contribution primaire est la dérivation explicite des équations du mouvement non linéaires pour un système multi-pendule et la linéarisation subséquente qui prend en compte les forces longitudinales nominales. En prouvant l'équivalence mathématique entre les modèles de pendule et masse-ressort-amortisseur, ce travail valide l'utilisation de ce dernier pour la conception de commande tout en l'ancrant dans la dynamique plus fondamentale du pendule. Les auteurs soulignent que les modèles proposés sont bien posés pour des applications spatiales réalistes où les angles de ballonnement restent dans des limites raisonnables, et la validation confirme la correction des formulations dérivées tant pour la simulation temporelle que pour l'analyse fréquentielle.