空力弾性解析

空力弾性はCollarの三角形で説明される。空気力（Aerodynamic Forces）、弾性力（Elastic Forces）、慣性力（Inertial Forces）の三者が結合する問題だ。この三者の相互作用がフラッター、ダイバージェンス、バフェッティングなどの危険な現象を引き起こす。

特にフラッター（flutter）は、空気力と構造弾性の結合による自励振動で、振幅が指数関数的に成長して構造破壊に至る。航空機設計では飛行包絡線内でフラッターが発生しないことの証明が型式証明の必須要件だよ。

典型的な翼型フラッターは、たわみ（plunge, $h$）とねじり（pitch, $\alpha$）の2自由度モデルで説明される。運動方程式は、

非定常空気力が速度の関数なので、速度を上げていくとたわみモードとねじりモードの周波数が近づき（周波数の合流, frequency coalescence）、ある速度でエネルギーの授受が正味プラスになって振動が発散する。この臨界速度がフラッター速度 $V_F$ だ。

古典的にはTheodorsen関数 $C(k)$ を使う。調和振動する翼型の非定常揚力は、

その通り。CFDベースの空力弾性解析では、非定常空気力をCFDから直接計算するから、理論モデルの適用範囲（線形、ポテンシャル流れ）を超えた遷音速フラッターや大振幅振動の予測が可能になる。

大きく3つに分けられる。

実務で最もよく使われるのはCFD + モーダル解析だ。構造を固有モード展開で表現し、各モードの一般化座標 $q_i(t)$ の時間発展をCFDの非定常空気力と連成して解く。

Ansys System Couplingを使えば、Fluent（流体）とAnsys Mechanical（構造）の間でメッシュ変位と面圧力を自動的にマッピングしてくれる。Nastranとの連成ではFluent UDF経由でモーダル座標を受け渡すカスタムワークフローを組むか、専用ツール（例えばMSC FlightLoads, Zona ZAERO）を使う方法がある。

線形フラッター解析のv-g法（velocity-damping method）とp-k法が基本だ。

V-g法: 各速度で調和振動を仮定し、必要な構造減衰 $g$ を求める。$g = 0$ を切る速度がフラッター速度。NastranのSOL 145がこの手法。

CFDの場合は非定常時刻歴をそのまま計算して、応答が減衰するか発散するかを観察する「時間進行法」が最も直接的だ。速度パラメータを段階的に上げて、フラッター境界を特定する。ただし計算コストが高いから、まずは線形理論（DLM: Doublet Lattice Method + Nastran SOL 145/146）でおおよその範囲を絞り、CFDで精緻化するのが効率的な手順だよ。

航空機のフラッター認証では今でもDLMが主力だ。ポテンシャル流れの仮定のもとで非定常空気力を周波数領域で効率的に計算できる。ただし遷音速域では衝撃波の効果を扱えないため、CFDで計算した定常圧力分布でDLMの空気力を補正する「CFD-corrected DLM」が広く使われているよ。

典型的なフローを段階的に示すよ。

Phase 1: 構造モデル構築

その通り。構造モデルの固有振動数とモード形状をGVTで検証するのは型式証明の必須プロセスだ。FEMモデルの固有振動数が実験と5%以上乖離していたらモデルのアップデートが必要になる。

いくつかの標準的なテストケースがある。

45度後退・テーパー比0.66の対称翼型（NACA 65A004）翼で、NASA Langleyの遷音速風洞で様々なマッハ数でフラッター試験が行われた。マッハ数0.499-1.141の範囲でフラッター速度指数 $V_F/b\omega_\alpha$ と周波数比 $\omega_F/\omega_\alpha$ が報告されている。CFDコミュニティでは必ずと言っていいほどこの問題で検証するよ。

用途に応じて使い分けるのが現実的だ。

Ansys Workbenchでのワークフローはこうだ。

極めて重要だ。空力弾性連成は付加質量効果（added mass effect）によって不安定になりやすい。特に軽い構造（翼の質量密度が流体密度に近い場合）では、under-relaxationなしでは連成反復が発散する。0.5程度から始めて、収束を確認しながら上げるのが安全だよ。

preCICEライブラリを使えば、OpenFOAMとCalculiX（またはFEniCS）を連成できる。preCICEはドイツのミュンヘン工科大学/シュトゥットガルト大学で開発された連成ライブラリで、空間マッピング、時間補間、安定化スキーム（Aitken加速、IQN-ILS）を提供する。

Interface Quasi-Newton with Inverse Least Squaresの略で、連成反復の収束を劇的に加速するアルゴリズムだ。過去の連成反復のデータからヤコビアンを近似構成し、準Newton法的に次の推定値を計算する。under-relaxation法と比べて反復回数を半分以下に減らせることが多い。preCICEの大きな強みだよ。

それがLCO（Limit Cycle Oscillation）だ。線形理論ではフラッター速度を超えると振幅が無限に成長するけど、実際には非線形効果（衝撃波の移動、境界層剥離、構造の非線形性）が振幅を制限する。遷音速域では衝撃波の位置が振動に伴って変化し、空気力の非線形性がLCOを引き起こすことが多い。

LCOの振幅予測は線形理論では不可能で、CFDベースの時間進行法が必要だ。CFDで20-50サイクル分の時間進行計算を行い、定常的な振幅に収束することを確認する。

必ずしもそうではない。LCOの振幅が大きければ構造疲労が問題になるし、サブクリティカルLCO（線形フラッター速度以下で発生するLCO）の存在も報告されている。航空機の認証ではあくまでフラッター速度に十分なマージンを持たせることが求められるよ。

ROM（Reduced Order Model）が有力だ。CFDの非定常空気力をコンパクトな伝達関数で近似し、構造方程式と効率的に連成する。

代表的なROM手法を紹介しよう。

CFDの非定常シミュレーションから数百のスナップショット（圧力場など）を取得し、PODで主要な空間モードを抽出する。そしてNavier-Stokes方程式をこのPODモード空間に射影して、数個の常微分方程式に縮約する。元のCFDが100万自由度なのに対して、ROMは10-50自由度のシステムになるから、フラッター速度のパラメータスタディが1000倍以上速くなる。

以下の方面で活発に研究されている。

その通り。認証に使うにはモデルの解釈可能性と信頼性の保証が必要で、現時点ではブラックボックス的なDNNは困難だ。ただし設計初期段階の探索や、風洞試験のオンラインデータ処理には実用化が始まっているよ。

空力弾性連成の発散原因のトップ3を示そう。

翼面が動くと、周囲の流体も連動して動く。この流体の慣性が「付加質量」として構造に作用する。密度比 $\rho_s / \rho_f$ が小さい（構造が軽い）ほど付加質量の影響が大きくなり、明示的な連成スキーム（1ステップで流体→構造→流体）では不安定になる。

対策は以下の通り。

メッシュ変形手法の選択が鍵だ。

大変形が予想される場合はOverset mesh（Chimera法）が最も安全だ。翼まわりのメッシュを独立した体として移動・回転させ、背景メッシュとの重なり領域で補間する。メッシュの変形は不要だから品質劣化がない。

遷音速域で20%の差は珍しくない。DLMはポテンシャル流れの仮定だから、衝撃波の効果が正確に反映されない。一方、CFDの結果も乱流モデルやメッシュ密度に依存する。

判断基準は以下の通りだ。

そう。ただしFAA/EASAもCFDベースのフラッター解析を認証に使うためのガイドライン整備を進めている。将来的にはCFDが主役になる可能性は高いよ。