[論文レビュー] Attitude Maneuvers of a Rigid Spacecraft in a Circular Orbit

Taeyoung Lee, Melvin Leok+1 more|ArXiv.org|Sep 14, 2005

Spacecraft Dynamics and Control参考文献 2被引用数 8

ひとこと要約

本稿は、回転行列をSO(3)に用い、特異点を回避するためのLie群変分積分法を用いて、円軌道上を飛行する剛体宇宙機の大きな角度の姿勢遷移を、幾何学的かつグローバルに有効なモデルとして提示する。非線形境界値問題および最小インパルス最適制御問題を解くために、根を求めることと逐次二次計画法の計算手法を構築し、so(3)のリー代数において感度微分を導出する。主たる貢献は、宇宙機の姿勢力学の内在的構造を保つ、数値的に安定で幾何学的に正確な手法である。

ABSTRACT

A global model is presented that can be used to study attitude maneuvers of a rigid spacecraft in a circular orbit about a large central body. The model includes gravity gradient effects that arise from the non-uniform gravity field and characterizes the spacecraft attitude with respect to the uniformly rotating local vertical local horizontal coordinate frame. An accurate computational approach for solving a nonlinear boundary value problem is proposed, assuming that control torque impulses can be applied at initiation and at termination of the maneuver. If the terminal attitude condition is relaxed, then an accurate computational approach for solving the minimal impulse optimal control problem is presented. Since the attitude is represented by a rotation matrix, this approach avoids any singularity or ambiguity arising from other attitude representations such as Euler angles or quaternions.

研究の動機と目的

• 円軌道上を飛行する剛体宇宙機の大きな角度の姿勢遷移を、重力勾配効果を考慮して、グローバルに有効で幾何学的なモデルとして開発すること。
• Euler角やクォータニオンに内在する特異点や曖昧さを回避するため、SO(3)に属する回転行列を用いて、従来の姿勢表現の限界を克服すること。
• 2種類の姿勢遷移問題を定式化・解法する：(1) 境界値問題による固定終端姿勢、(2) 終端条件を緩和した最小インパルス制御。
• 状態遷移行列を用いて解析的感度微分を計算することで、幾何的構造を保ちながら計算の正確性を確保すること。
• 初期姿勢から軌道の1/4を通過する姿勢遷移の数値解を提示し、手法の有効性を実証すること。

提案手法

• 宇宙機の姿勢を、本体固定座標系と局所垂直局所水平（LVLH）座標系の間の変換を表す回転行列 $ R^{bl} \in \mathrm{SO(3)} $ で表現する。
• SO(3)の幾何的構造を保存し、不変量を保存するLie群変分積分法を用いて、連続的な姿勢力学をモデル化する。
• 重力ポテンシャルから導出される重力勾配モーメントを組み込み、LVLH座標系内での非一様重力場効果を捉える。
• 状態遷移行列を用いて解析的感度微分を計算する根を求めることによる非線形境界値問題の解法。
• 性能指標および制約の勾配をso(3)で解析的に導出し、逐次二次計画法（SQP）を用いて最小インパルス最適制御問題を解法する。
• 感度微分をリー代数 $ \mathfrak{so}(3) $ の言語で表現することで、幾何的整合性を確保し、Euler角やクォータニオンに内在する特異点を回避する。

実験結果

リサーチクエスチョン

• RQ1 Euler角やクォータニオンに内在する特異点や曖昧さを回避しつつ、円軌道上を飛行する剛体宇宙機の大きな角度の姿勢遷移を、どのように幾何学的に正確にモデル化できるか？
• RQ2 宇宙機の姿勢制御における非線形境界値問題を解くために、最も正確で幾何学的に整合性のある計算的手法は何か？
• RQ3 与えられた終端制約を満たしつつ、姿勢遷移に要する最小総インパルスをどのように計算できるか？
• RQ4 性能指標および制約の感度微分が最適化プロセスの収束性と正確性にどのように寄与するか？
• RQ5 Lie群変分積分法を用いることで、長時間シミュレーションにおいてSO(3)の内在的幾何的構造はどの程度正確に保存されるか？

主な発見

• 最小インパルス最適制御問題は、制約違反が $ 4.80 \times 10^{-14} $ にまで低減され、高い数値的正確性を示した。
• 最適化された初期角運動量は $ \Pi_0 = [-2.915, -2.347, -2.734]^T $ であり、終端角運動量は $ -\Pi_N = [-0.343, 2.686, 2.734]^T $ となり、所望の姿勢が達成された。
• 最終的な回転行列 $ R^{bl}_N $ により得られた終端方向ベクトル $ \Lambda_N = [0, 0, -1]^T $ は、座標系の向きの補正後、所望の $ \Lambda_{N_d} = [0, -1, 0]^T $ と一致した。
• 性能指標は $ \mathcal{J} = 6.771 $ となり、この値は遷移に要する総インパルスの大きさを表した。
• この手法は幾何学的に正確で、数値的にも効率的であり、$ \mathfrak{so}(3) $ で解析的に導出した感度微分により特異点を回避した。
• 重力勾配がゼロ（$ \omega_0 = 0 $）の場合、このフレームワークは自由剛体のケースに還元され、古典的力学と整合性があることが確認された。