Equivalent Mechanical Models for Sloshing

Technische Samenvatting: Equivalent Mechanische Modellen voor Sloshing

Probleemstelling
Propellant sloshing (het klotsen van brandstof) blijft een kritieke uitdaging voor de stabiliteit en besturing van ruimtevaartuigen, met name in twee verschillende regimes: microzwaartekrachtomgevingen die een hoge richtingsstabiliteit vereisen (bijv. interplanetaire missies) en situaties met hoge versnelling ($g$) tijdens lancering, terugkeer in de atmosfeer en landing. Hoewel pendelgebaseerde en massa-veersysteem-dempermodellen standaardinstrumenten zijn voor controle-analisten om sloshing-geïnduceerde verstoringen te beoordelen, is een rigoureuze wiskundige formulering nodig die deze representaties verbindt, specifiek rekening houdend met de effecten van nominale longitudinale krachten in hoge-$g$-regimes, voor robuuste guidance, navigation, and control (GNC) ontwerp.

Methodologie
Het artikel presenteert een rigoureuze afleiding van equivalente mechanische modellen voor sloshing, waarbij wordt geprogresseerd van nietlineaire dynamica naar geliniariseerde representaties die geschikt zijn voor controle-analyse.

1. Nietlineaire Formulering:

   • De studie begint met een enkele sferische pendel met massa $m_P$ en lengte $\ell$, bevestigd aan een star lichaam $B$ met massa $m_B$.
   • Met behulp van een Newtoniaanse benadering leiden de auteurs de gekoppelde bewegingsvergelijkingen af voor het 8-vrijheidsgraden systeem (translatie/rotatie van het starre lichaam + pendelhoeken).
   • De afleiding behandelt expliciet de trekkracht $N$ en elimineert deze om de systeemvariabelen te reduceren, wat resulteert in een compacte matrixvorm $\mathcal{M}\ddot{\mathbf{x}} = \mathbf{u}$.
   • De formulering wordt uitgebreid om meerdere pendels te accommoderen, die meerdere tanks of sloshing-modi vertegenwoordigen, waarbij gebruik wordt gemaakt van een blokgestructureerde massamatrix om de ijdelheid te benutten voor een efficiënte numerieke oplossing.
   • Een viskeus dempingsmodel wordt geïntegreerd door termen toe te voegen die proportioneel zijn aan de hoeksnelheden, welke de pendeldynamica beïnvloeden zonder de starre-lichaamvergelijkingen direct te wijzigen.

2. Linearisatie en Modale Analyse:

   • In het besef dat onboard controlesystemen doorgaans het zwaartepunt schatten uitgaande van een starre brandstofverdeling, herformuleren de auteurs de dynamica ten opzichte van een nominaal globaal zwaartepunt ($\bar{G}$) in plaats van het zwaartepunt van het starre lichaam.
   • Het systeem wordt gelinieariseerd rond een evenwicht dat wordt gekenmerkt door een constante attitude uitgelijnd met het inertiële kader en een voorgeschreven nominale longitudinale kracht $\bar{F}_z$.
   • Twee gevallen worden geanalyseerd: krachten gedefinieerd in het inertiële kader en krachten gedefinieerd in het lichaamsvast kader.
   • Door middel van coördinatentransformaties en algebraïsche manipulatie wordt de gelinieariseerde dynamica gegoten in een "modale vorm". Dit houdt in dat nieuwe modale coördinaten ($\eta_\theta, \eta_\phi$) worden gedefinieerd die de pendeldynamica ontkoppelen in twee onafhankelijke planaire oscillatoren.
   • Het resulterende systeem wordt uitgedrukt als een tweede-orde lineair MIMO-systeem ($M\delta\ddot{\mathbf{x}} + K\delta\mathbf{x} = \delta\mathbf{u}_{ext}$), waarbij de modale participatiematrix en de stijfheidstermen geïnduceerd door de nominale kracht expliciet worden geïdentificeerd.

3. Equivalentie met Massa-Veer-Demper Modellen:

   • Het artikel leidt de gelinieariseerde vergelijkingen af voor een massa-veer-dempersysteem waarbij de sloshing-massa zijdelings wordt beperkt door veren met stijfheid $k$.
   • Door de veerstijfheid $k = m_P \omega_0^2$ te stellen (waarbij $\omega_0$ de natuurlijke frequentie is afgeleid van het pendelmodel), demonstreren de auteurs dat de gelinieariseerde massa-veer-demperdynamica wiskundig identiek is aan de gelinieariseerde pendeldynamica.
   • De analyse benadrukt dat de aanwezigheid van een nominale longitudinale kracht $\bar{F}_z$ essentieel is voor de pendel om oscillerend gedrag te vertonen analoog aan een veer; bij afwezigheid van deze kracht ($\bar{F}_z=0$) vereist het pendelmodel een modificatie (bijv. het toevoegen van een derde modale vrijheidsgraad) om 0-g sloshing te representeren.

Belangrijkste Resultaten

• Wiskundige Equivalentie: De studie vestigt formeel de exacte equivalentie tussen het gelinieariseerde pendelmodel en het klassieke massa-veer-dempermodel onder hoge-$g$-condities, waarbij wordt verduidelijkt hoe de nominale kracht fungeert als het herstellend mechanisme in de pendelformulering.
• Modale Representatie: De afleiding levert een heldere modale representatie waarbij de sloshing-modi ontkoppeld zijn van de starre-lichaam-dynamica in de stijfheidsmatrix, waarbij koppeling alleen optreedt in de massamatrix via de modale participatiematrix.
• Validatie:
  • Nietlineaire Validatie: Het nietlineaire model werd gevalideerd tegen een MATLAB Simscape Multibody implementatie (starre hub met vier pendels). Tijddomein-simulaties vertoonden uitstekende overeenstemming, met maximale fouten onder de 0,002% van de responsamplitude, waarbij discrepanties uitsluitend werden toegeschreven aan verschillen in numerieke integratieschema's.
  • Lineaire Validatie: Bode-plots die de analytische gelinieariseerde modellen vergelijken met de gelinieariseerde Simscape-plant toonden uitstekende overeenstemming over het frequentiebereik voor zowel inertiële als lichaamsvasse krachtinputs. Kleine laagfrequente discrepanties in het inertiële kader-geval werden toegeschreven aan de numerieke nauwkeurigheidslimieten van het Simscape-linearisatieproces.

Betekenis en Claims
Het artikel beweert een rigoureus wiskundig fundament te bieden voor het modelleren van propellant sloshing in omgevingen met hoge versnelling. De primaire bijdrage is de expliciete afleiding van de nietlineaire bewegingsvergelijkingen voor een multi-pendelsysteem en de daaropvolgende linearisatie die rekening houdt met nominale longitudinale krachten. Door de wiskundige equivalentie tussen het pendel- en het massa-veer-dempermodel te bewijzen, valideert het werk het gebruik van dat laatste voor controledesign, terwijl het geworteld blijft in de meer fundamentele pendeldynamica. De auteurs benadrukken dat de voorgestelde modellen goed gedefinieerd zijn voor realistische ruimtevaarttoepassingen waarbij sloshing-hoeken binnen redelijke grenzen blijven, en de validatie bevestigt de correctheid van de afgeleide formuleringen voor zowel tijddomein-simulatie als frequentiedomein-analyse.