Equivalent Mechanical Models for Sloshing

Technische Zusammenfassung: Äquivalente mechanische Modelle für Sloshing (Schwappen)

Problemstellung
Das Schwappen von Treibstoffen bleibt eine kritische Herausforderung für die Stabilität und Regelung von Raumfahrzeugen, insbesondere in zwei unterschiedlichen Regimen: Mikrogravitationsumgebungen, die eine hohe Punktgenauigkeit erfordern (z. B. interplanetare Missionen), und Bedingungen mit hoher Beschleunigung ($g$), wie sie während des Starts, des Wiedereintritts und der Landung auftreten. Während pendelbasierte oder Masse-Feder-Dämpfer-Modelle Standardwerkzeuge für Steuerungsanalysten zur Bewertung von durch Sloshing induzierten Störungen sind, ist eine rigorose mathematische Formulierung, die diese Darstellungen verbindet – insbesondere unter Berücksichtigung der Auswirkungen nominaler longitudinaler Kräfte in High-$g$-Regimen – für ein robustes Guidance, Navigation, and Control (GNC)-Design notwendig.

Methodik
Die Arbeit präsentiert eine rigorose Herleitung äquivalenter mechanischer Modelle für Sloshing, ausgehend von nichtlinearen Dynamiken hin zu linearisierten Repräsentationen, die für die Steuerungsanalyse geeignet sind.

1. Nichtlineare Formulierung:

   • Die Studie beginnt mit einem einfachen kugelförmigen Pendel der Masse $m_P$ und der Länge $\ell$, das an einen starren Körper $B$ der Masse $m_B$ gekoppelt ist.
   • Unter Verwendung eines Newtonschen Ansatzes leiten die Autoren die gekoppelten Bewegungsgleichungen für das 8-Freiheitsgrade umfassende System (Translation/Rotation des starren Körpers + Pendelwinkel) ab.
   • Die Herleitung behandelt explizit die Seilspannung $N$ und eliminiert diese, um die Systemvariablen zu reduzieren, was zu einer kompakten Matrixform $\mathcal{M}\ddot{\mathbf{x}} = \mathbf{u}$ führt.
   • Die Formulierung wird erweitert, um mehrere Pendel aufzunehmen, welche mehrere Tanks oder Sloshing-Moden repräsentieren, wobei eine blockstrukturierte Massenmatrix genutzt wird, um die Sparsity (Dünnbesetztheit) für eine effiziente numerische Lösung auszunutzen.
   • Ein viskoses Dämpfungsmodell wird durch Hinzufügen von Termen implementiert, die proportional zu den Winkelgeschwindigkeiten sind, welche die Pendeldynamik beeinflussen, ohne die Starrkörpergleichungen direkt zu verändern.

2. Linearisierung und Modalanalyse:

   • In Anerkennung der Tatsache, dass bordgestützte Steuerungssysteme den Schwerpunkt unter der Annahme einer starren Treibstoffverteilung schätzen, formulieren die Autoren die Dynamik relativ zu einem nominalen globalen Schwerpunkt ($\bar{G}$) statt zum Schwerpunkt des starren Körpers um.
   • Das System wird um ein Gleichgewicht linearisiert, das durch eine konstante Orientierung parallel zum Inertialrahmen und eine vorgegebene nominale longitudinale Kraft $\bar{F}_z$ charakterisiert ist.
   • Zwei Fälle werden analysiert: Kräfte, die im Inertialrahmen definiert sind, und Kräfte, die im körperfesten Rahmen definiert sind.
   • Durch Koordinatentransformationen und algebraische Manipulationen werden die linearisierten Dynamiken in eine „Modalform“ überführt. Dies beinhaltet die Definition neuer modaler Koordinaten ($\eta_\theta, \eta_\phi$), die die Pendeldynamik in zwei unabhängige planare Oszillatoren entkoppeln.
   • Das resultierende System wird als ein zweiter Ordnung MIMO-System ($M\delta\ddot{\mathbf{x}} + K\delta\mathbf{x} = \delta\mathbf{u}_{ext}$) ausgedrückt, wobei die modale Partizipationsmatrix und die durch die nominale Kraft induzierten Steifigkeitsterme explizit identifiziert werden.

3. Äquivalenz mit Masse-Feder-Dämpfer-Modellen:

   • Die Arbeit leitet die linearisierten Gleichungen für ein Masse-Feder-Dämpfer-System her, bei dem die Slosh-Masse über Federn mit der Steifigkeit $k$ lateral begrenzt ist.
   • Durch Setzen der Federsteifigkeit $k = m_P \omega_0^2$ (wobei $\omega_0$ die aus dem Pendelmodell abgeleitete Eigenfrequenz ist) zeigen die Autoren, dass die linearisierten Masse-Feder-Dämpfer-Dynamiken mathematisch identisch mit den linearisierten Pendeldynamiken sind.
   • Die Analyse hebt hervor, dass die Anwesenheit einer nominalen longitudinalen Kraft $\bar{F}_z$ essenziell ist, damit das Pendel eine oszillatorische Bewegung analog zu einer Feder ausführt; in Abwesenheit dieser Kraft ($\bar{F}_z=0$) erfordert das Pendelmodell eine Modifikation (z. B. das Hinzufügen eines dritten modalen Freiheitsgrades), um 0-g Sloshing darzustellen.

Hauptergebnisse

• Mathematische Äquivalenz: Die Studie stellt formal die exakte Äquivalenz zwischen dem linearisierten Pendelmodell und dem klassischen Masse-Feder-Dämpfer-Modell unter High-$g$-Bedingungen her und klärt dabei, wie die nominale Kraft als Rückstellmechanismus in der Pendelformulierung wirkt.
• Modale Repräsentation: Die Herleitung liefert eine klare modale Repräsentation, in der die Slosh-Moden in der Steifigkeitsmatrix von der Starrkörperdynamik entkoppelt sind, wobei die Kopplung nur in der Massenmatrix über die modale Partizipationsmatrix auftritt.
• Validierung:
  • Nichtlineare Validierung: Das nichtlineare Modell wurde gegen eine MATLAB Simscape Multibody Implementierung (starrer Hub mit vier Pendeln) validiert. Zeitbereichssimulationen zeigten eine exzellente Übereinstimmung mit maximalen Fehlern von unter 0,002 % der Antwortamplitude, wobei Diskrepanzen allein auf Unterschiede in den numerischen Integrationsverfahren zurückgeführt wurden.
  • Lineare Validierung: Bode-Diagramme, die das analytische linearisierte Modell mit der linearisierten Simscape-Anlage vergleichen, zeigten eine exzellente Übereinstimmung über den gesamten Frequenzbereich sowohl für Inertialrahmen- als auch für körperfeste Krafteingänge. Geringfügige Abweichungen im niedrigen Frequenzbereich im Fall des Inertialrahmens wurden auf die numerischen Genauigkeitsgrenzen des Simscape-Linearisierungsprozesses zurückgeführt.

Bedeutung und Ansprüche
Die Arbeit beansprucht, eine rigorose mathematische Grundlage für die Modellierung von Treibstoff-Sloshing in Umgebungen mit hoher Beschleunigung zu liefern. Ihr primärer Beitrag ist die explizite Herleitung der nichtlinearen Bewegungsgleichungen für ein Mehrpendelsystem und die anschließende Linearisierung unter Berücksichtigung nominaler longitudinaler Kräfte. Durch den Nachweis der mathematischen Äquivalenz zwischen dem Pendel- und dem Masse-Feder-Dämpfer-Modell validiert die Arbeit die Verwendung des letzteren für das Steuerungsdesign und bettet es gleichzeitig in die grundlegendere Pendeldynamik ein. Die Autoren betonen, dass die vorgeschlagenen Modelle für realistische Raumfahrtanwendungen, bei denen die Slosh-Winkel innerhalb vernünftiger Grenzen bleiben, wohldefiniert sind, und die Validierung bestätigt die Korrektheit der abgeleiteten Formulierungen sowohl für die Zeitbereichssimulation als auch für die Frequenzbereichsanalyse.