Equivalent Mechanical Models for Sloshing

技術要約：スロッシングに対する等価な機械モデル

問題提起
推進剤のスロッシング（液面揺動）は、宇宙機の安定性と制御における極めて重要な課題であり、特に、高いポインティング安定性が求められる微小重力環境（例：惑星間航行ミッション）と、打ち上げ、再突入、着陸時に遭遇する高加速度（$g$）条件という、2つの異なる領域において重要となる。振り子モデルや質量・ばね・ダンパモデルは、制御解析者がスロッシングによる摂動を評価するための標準的なツールであるが、これらの表現を、特に高-$g$領域における公称縦方向力の効果を考慮して厳密に結びつける数学的定式化、すなわち、ロバストな誘導・航法・制御（GNC）設計のための数学的定式化が必要である。

手法
本論文は、非線形動力学から制御解析に適した線形表現へと進展する、スロッシングの等価な機械モデルの厳密な導出を提示している。

1. 非線形定式化:

   • 本研究は、質量 $m_P$、長さ $\ell$ の単一の球形振り子が、質量 $m_B$ の剛体 $B$ に取り付けられている状態から始まる。
   • 著者は、ニュートン力学の手法を用い、8自由度システム（剛体の並進・回転 ＋ 振り子の角度）の結合運動方程式を導出している。
   • 導出過程では、張力 $N$ を明示的に扱い、変数削減のためにこれを排除することで、コンパクトな行列形式 $\mathcal{M}\ddot{\mathbf{x}} = \mathbf{u}$ を得ている。
   • この定式化は、複数のタンクまたは複数のスロッシングモードを表す複数の振り子に対応できるよう拡張されており、数値解法を効率化するためにブロック構造化された質量行列を利用している。
   • 粘性減衰モデルは、角速度に比例する項を加えることで組み込まれており、これらは剛体の方程式を直接変更することなく、振り子の動力学に影響を与える。

2. 線形化およびモード解析:

   • オンボード制御システムは通常、推進剤の分布が剛体であると仮定して重心を推定するため、著者は、剛体の重心ではなく、公称のグローバル重心（$\bar{G}$）を基準とした動力学の再定式化を行っている。
     kapsamında、慣性系に整列した一定の姿勢と、規定の公称縦方向力 $\bar{F}_z$ によって特徴付けられる平衡状態の周りでシステムを線形化している。
   • 力の定義として、慣性系で定義される場合と、機体固定座標系で定義される場合の2つのケースを分析している。
   • 座標変換と代数的な操作を通じて、線形化された動力学を「モード形式」へとキャストしている。これには、振り子の動力学を2つの独立した平面振動子へとデカップリングする、新しいモード座標（$\eta_\theta, \eta_\phi$）の定義が含まれる。
   • 結果として得られるシステムは、2次線形MIMOシステム（$M\delta\ddot{\mathbf{x}} + K\delta\mathbf{x} = \delta\mathbf{u}_{ext}$）として表され、モード参加行列と、公称力によって誘発される剛性項を明示的に特定している。

3. 質量・ばね・ダンパモデルとの等価性:

   • 本論文では、スロッシング質量がばね剛性 $k$ を持つばねを介して横方向に移動するように制約された、質量・ばね・ダンパシステムの線形化方程式を導出している。
   • ばね剛性を $k = m_P \omega_0^2$（ここで $\omega_0$ は振り子モデルから導かれる自然周波数）と設定することで、著者らは、線形化された質量・ばね・ダンパの動力学が、線形化された振り子の動力学と数学的に同一であることを示している。
   • この解析は、振り子定式化において、公称縦方向力 $\bar{F}_z$ が存在することが、ばねのような振動挙動を示すために不可欠であることを強調している。この力（$\bar{F}_z=0$）が存在しない場合、振り子モデルは（例：0-gスロッシングを表現するために第3のモード自由度を追加するなど）修正を必要とする。

主な結果

• 数学的等価性: 本研究は、高加速度条件下における線形化された振り子モデルと古典的な質量・ばね・ダンパモデルの間の正確な等価性を正式に確立しており、公称力が振り子定式化においてどのように復元メカニズムとして機能するかを明らかにしている。
• モード表現: 導出されたモード表現は、剛性行列においてスロッシングモードが剛体動力学からデカップリングされ、結合がモード参加行列を介して質量行列にのみ現れる、簡潔な形式となっている。
• 検証:
  • 非線形検証: 非線形モデルは、MATLAB Simscape Multibody実装（4つの振り子を持つ剛体ハブ）に対して検証された。時間領域シミュレーションは優れた一致を示し、最大誤差は応答振幅の0.002%未満であった。不一致は、数値積分スキームの違いのみに起因している。
  • 線形検証: 解析的な線形化モデルと、線形化されたSimscapeプラントを比較したボード線図は、慣性系および機体固定系の両方の力入力に対して、全周波数範囲にわたって優れた一致を示した。慣性系の場合における軽微な低周波の不一致は、Simscapeの線形化プロセスにおける数値精度の限界に起因している。

意義および主張
本論文は、高加速度環境における推進剤スロッシングをモデリングするための厳密な数学的基礎を提供することを主張している。その主要な貢献は、複数の振り子システムに対する非線形運動方程式の明示的な導出と、それに続く、公称縦方向力を考慮した線形化である。振り子モデルと質量・ばね・ダンパモデルの間の数学的等価性を証明することにより、本研究は、より基礎的な振り子動力学に基づきつつ、制御設計のための後者の使用を正当化している。著者らは、提案されたモデルがスロッシング角が妥当な範囲内に留まる現実的な宇宙船アプリケーションに対して適切であることを強調しており、検証によって、時間領域シミュレーションおよび周波数領域解析の両方における導出された定式化の正当性が確認されている。