[論文レビュー] A Flight-Mechanics Solver for Aircraft Inverse Simulations and Application to 3D Mirage-III Maneuver

Osama A. Marzouk|arXiv (Cornell University)|Oct 29, 2024

Aerospace and Aviation Technology被引用数 1

ひとこと要約

本論文は、風圏およびボディ座標系における連立微分代数方程式（DAE）を用いた3次元非線形飛行力学ソルバーを提案し、固定翼機の逆問題シミュレーションを実行する。4次ルンゲ・クッタ法を用いて、所定の飛行パタンに対して制御面の舵角および推力プロファイルを予測する。ミrage-IIIのロールマニューバーを実証例として用い、時間刻みの感度が低く、逆推力を必要としない高精度で安定した結果を得た。

ABSTRACT

The main objective of this paper is to present a general mathematical model and an associated numerical algorithm applicable to an arbitrary fixed-wing fixed-mass aircraft undergoing an arbitrary maneuver, based on the 3D nonlinear coupled differential-algebraic equations of motion, including force, moment, kinematic and constraint equations. The model is formulated to address the inverse simulation problem where a target maneuver is prescribed and the corresponding time dependent patterns of the control variables are solved for to meet this maneuver. The model utilizes two different moving frames of references, namely the body axes and the wind axes. The numerical algorithm features sequential solution of equations in a fully explicit manner. It is straightforward to use the model in a reverse mode, namely the direct simulation problem. The inverse problem may be summarized as follows: Inputs: Time history of desired-trajectory rectangular coordinates relative to the ground-fixed axes. A constraint should be specified, which we arbitrarily chose it to be the bank angle. Also, certain geometric and aerodynamic aircraft data are needed. Outputs: Time history of the control variables (thrust magnitude, elevator angle, rudder angle, ailerons angle), which will satisfy the aimed trajectory. The paper finally applies the presented numerical algorithm to a roll maneuver for the Mirage-III fighter.

研究の動機と目的

• 複雑なマニューバーを遂行する任意の固定翼機に対して汎用的な3次元飛行力学モデルを開発すること。
• 所定の軌道を達成するために必要な時間変化する制御面の舵角および推力を特定することにより、逆問題シミュレーションを解くこと。
• 並進運動を風圏座標系に、回転運動をボディ座標系に定式化することで、計算効率と数値的安定性を確保すること。
• 実戦機（ミrage-III）が継続的ロールマニューバーを実行する状況で、モデルを検証すること。
• スラスト要求、制御面の極端な舵角、迎角範囲を評価するための設計ツールを提供すること。

提案手法

• 6自由度運動における力、モーメント、運動学的および制約方程式を組み合わせた18個の非線形連立DAEを定式化する。
• 線形運動量方程式に風圏座標系を用いることで、数値的剛性を低減し、計算効率を向上させる。
• 角運動量方程式にボディ座標系を用いることで、慣性モーメントを一定に保ち、時間依存項を回避する。
• 4次ルンゲ・クッタ法を用いた逐次的完全明示的数値積分法により、時間に沿ったDAE系の解を求める。
• 実際の力とモーメントのモデル化のため、空気力学的安定性導関数および大気の性質（密度、重力加速度）を組み込む。
• MATLABでアルゴリズムを実装し、ミrage-IIIの6秒間の継続的ロールマニューバーを用いて検証した。

実験結果

リサーチクエスチョン

• RQ1 統一的な3次元飛行力学モデルは、任意の固定翼機のマニューバーに対して逆問題シミュレーションを高精度で解くことができるか？
• RQ2 並進運動に風圏座標系を用いることで、ボディ座標系に比べて数値的安定性がどのように向上するか？
• RQ3 ミrage-III機において、正確な3次元ロールマニューバーを達成するために必要な制御面の舵角および推力プロファイルは何か？
• RQ4 解は時間刻みの大きさにどれほど敏感であり、標準的な計算リソースで達成可能な精度の水準は何か？
• RQ5 モデルはスラスタの領域に進入しないように信頼性高く飛行パラメータを予測できるか？迎角および制御面の舵角の範囲の上限は何か？

主な発見

• ミrage-IIIのロールマニューバーにおける推力プロファイルは対称的であり、全期間にわたり正のまま維持され、逆推力は不要であった。
• ラダーが最大で49.9°の舵角を示し、マニューバー中におけるロール制御における主導的役割を示した。
• 実際の迎角は–6.05°から6.36°の間で変動し、典型的なスラスタ角15°よりも著しく低く、線形空気力学的モデルの妥当性が裏付けられた。
• ピッチ角とヨー角はそれぞれほぼ対称的および反対称的応答を示し、ピッチ・ヘディング平面における安定した位相図が得られた。
• 時間刻みが10⁻⁴ sまで小さくなると解は感度が低く、2×10⁻⁴ sおよび10⁻³ sでもわずかなずれにとどまり、10⁻² sでは顕著な誤差が生じた。
• モデルは、軌道座標およびバンク角の4つの所定入力から、制御面の舵角および推力を含む14個の未知飛行変数を正確に予測できた。