空力弾性 — CAE用語解説

良い質問だ。航空宇宙工学の分野では、空力弾性は「空気力、弾性力、慣性力の相互作用によって生じる現象」と厳密に定義される。代表的な現象はフラッター（発散振動）やダイバージェンス（静的発散）だ。例えば、飛行機の主翼が特定の速度で振動を始め、振幅が発散して破壊に至るフラッターは、1940年のタコマナローズ橋の崩壊でも知られる。

その「別々に解く」というのがポイントで、完全に連成した単一の方程式は存在しない。構造の運動方程式と流体のNavier-Stokes方程式を、界面での変位と圧力/せん断応力を条件として結合して解く。構造の運動は、例えば2自由度の翼断面モデルでは次のように表される。

まさにそこが核心だ。流体の支配方程式を解く領域の境界、つまり翼表面の形状と速度が、構造の変位と速度によって時間とともに変化する。これが「移動境界」問題を生み出す。非定常の空気力を計算するには、例えばシンプルな線形理論では、シアーズ関数やテオドーセン関数といった既知の関数を使うが、実際のCAEでは、メッシュを変形させながらNavier-Stokes方程式を逐次解く。

その通り。線形系として扱える範囲では、系の固有値の実部が正になる条件が不安定性の閾値だ。例えば、ダイバージェンス速度 V_D は、

ほとんどの実務では「パーティションド（分割）連成解法」が使われる。つまり、構造ソルバー（Abaqus/Standardなど）と流体ソルバー（Fluent, Star-CCM+など）を別々に走らせ、界面でデータ（変位⇔圧力）を交換する。交換のタイミングには、明示的スキームと暗黙的スキームがある。明示的スキームは計算が簡単だが、安定性に問題があり、小さな時間ステップが必要だ。

そのリスクが常にある。対策は主に二つだ。一つは「動的メッシュ」技術。ANSYS Fluentでは、拡張性の高いスプリング近似法や、再メッシュを伴う動的レイヤー法を使う。もう一つは「任意ラグランジュ・オイラー（ALE）定式化」で、メッシュを構造の変位に合わせて変形させつつ、流体方程式を解く。変形が大きすぎる場合は、メッシュを一旦作り直す「リメッシュ」が必要になる。

鋭い指摘だ。構造の振動（例えば10-100Hz）に比べ、流体の渦の剥離などの現象はより高い周波数成分を含む。一般的に、流体側がボトルネックとなり、より小さな時間ステップを要求する。実務では、流体解析のCFL条件（例えばCFL数<1）から決まるステップ幅、例えば1e-4秒程度を、連成解析全体のステップとして採用することが多い。構造解析はそのステップで十分解ける。

その感覚は正しい。最初にやるべきは「非連成検証」だ。手順はこうだ：1) 構造のみの固有値解析（モード解析）で、対象とする振動モード（例えば第1ねじりモード、第2曲げモード）とその固有振動数を確認する。2) 剛体モデルで非定常CFD解析を行い、空力荷重の周波数応答を確認する。それぞれのソルバーが単体で正しく動くことを確かめてから、連成に進む。

それでは非現実的な過渡応答が生じる。実用的な手順は「段階的起動」だ。まず、構造を固定した状態で定常CFD解析を収束させ、初期の流れ場（圧力分布）を構築する。次に、その流れ場を初期条件として、構造の自由度を解放して非定常連成解析を開始する。これにより、いきなり大きな振動が発生するのを防げる。ANSYSのWorkbenchでFluentとTransient Structuralを連成させる場合、この設定は「System Coupling」コンポーネントで行う。

実験が最良ではあるが、必ずしも必須ではない。いくつかのベンチマークがある。例えば、NASAの「Standard Dynamics Model」や、AGARD（Advisory Group for Aerospace Research & Development）が定めた一連のフラッター試験ケース（AGARD 445.6翼など）だ。これらの幾何形状、材料特性、流れ条件は公開されており、あなたの解析結果（フラッター速度や周波数）を既存の文献値や他ソルバーの結果と比較できる。また、メッシュ依存性調査（メッシュを2倍細かくして結果が変わらないか確認）は必須だ。

連成解析の発散は、以下の順で切り分ける。1) **データ転送**: インターポレーション設定がおかしく、力や変位が異常値になっていないか。2) **時間ステップ**: 明示的スキームを使っている場合、ステップ幅が大きすぎないか。CFL数は1以下か。3) **メッシュ変形**: 動的メッシュの設定が不適切で、メッシュの歪みが大きくなりすぎていないか。4) **物理モデル**: ダンピング（構造減衰）を設定し忘れていないか。現実の構造には必ず減衰があるが、モデルに入れないと発散しやすくなる。

いや、選択肢はいくつかある。大きく分けて、**統合型ソルバー**と**連成プラットフォーム型**だ。統合型は、COMSOL Multiphysicsが代表的で、一つの環境内で構造と流体（および他の物理場）を直接連成して解く。もう一つは、ANSYS（Fluent + Mechanical）、Siemens（Star-CCM+ + Simcenter 3D）、Dassault（XFlow + Abaqus）のように、専用の連成環境（ANSYS System Coupling, SIMULIA Co-Simulation Engine）で異なるソルバーを結びつける方式だ。

数値解法の根幹が違う。COMSOLは「モノリシック解法」に近く、構造と流体の未知数を一つの大規模なマトリックスに集約し、直接連立方程式として解く。利点は強連成がデフォルトで安定性が高いこと。欠点はメモリ消費が大きく、専用の流体/構造ソルバーに比べて個別の機能（乱流モデルなど）が劣ることがある。一方、ANSYSのパーティションド連成は、各ソルバーの高度な機能を活かせるが、連成の安定性設定（under-relaxationなど）にユーザーの経験が求められる。

歴史的に、専用の社内開発コードも多いが、商用ソフトでは**MSC Nastran**が構造側で非常に強い。そのAeroelasticソリューション（SOL 144, SOL 146）は、線形空力理論（ダブルト・ラティス法など）と組み合わせたフラッター解析のデファクトスタンダードだ。高精度な非線形解析が必要な場合は、Nastran（構造）とCFDソルバー（例えばCFD++、OVERFLOW）を連成させる。SiemensのStar-CCM+も、その多物理機能（共形メッシュによるFSI）で、航空エンジンや自動車のミラー振動などの分野で採用が増えている。

あるが、ハードルは高い。典型的な組み合わせは、構造ソルバーに**CalculiX**（Abaqus-like）や**Code_Aster**、流体ソルバーに**OpenFOAM**を使う。連成には**preCICE**というライブラリが注目されている。preCICEは、異なるソルバー間のデータ交換と時間ステップ管理を担う。ただし、商用ソフトのようなGUIベースの統合環境はなく、全ての設定（メッシュマッピング、スキーム選択）を自分で行う必要があり、研究者や高度なユーザー向けだ。

これは「剛体運動モードの発生」または「初期過渡応答の暴走」が疑われる。まず、構造モデルに適切な拘束がかけられているか確認せよ。わずかなローディングでも、剛体モードがあれば無限大に動く。次に、先述の「段階的起動」ができているか。定常流れ場から始めていないと、翼にいきなり大きな非対称な力が加わり、暴走する。ANSYS System Couplingでは、「Coupling Initialization」で定常結果を読み込む設定がある。

動的メッシュの設定不足だ。Fluentの場合、1) 「Dynamic Mesh」設定で、変形する領域に適切な方法（Smoothing, Remeshing, Layering）を割り当てる。2) 「Smoothing」のスプリング定数係数を下げる（デフォルト1.0を0.5などに）。3) 「Remeshing」を有効にし、メッシュの歪みやサイズの閾値を設定する。また、構造の変位が大きすぎないか確認し、必要なら物理的なストッパー（接触条件）を構造モデルに追加する。

それは「安定」な状態だ。考えられる原因は二つ。第一に、**解析している流速がフラッター速度より十分に低い**。速度を段階的に上げていき、振動の減衰率（ログ減衰）がゼロになる点を探す必要がある。第二に、**数値的減衰が大きすぎる**。流体ソルバーの時間離散化スキーム（例えば1次風上差分）や構造ソルバーの数値減衰（アルファ法のパラメータ）が、物理的な減衰以上に効いている可能性がある。まずは、既知のベンチマークケースでソルバー設定を検証せよ。

ボトルネックを特定する必要がある。まず、**流体解析**が遅いことがほとんどだ。対策：1) メッシュ数を必要最小限に減らす（境界層は重要だが、遠方場は粗く）。2) 時間ステップ幅を可能な限り大きくする（ただしCFL条件は守る）。3) 並列計算コア数を増やす（FluentやStar-CCM+は並列効率が高い）。次に、**データ転送の頻度**を見直す。連成ステップごとにデータ交換するのではなく、流体側で数ステップ進めてから平均化した荷重を構造側に渡す「サブサイクリング」が有効な場合がある。ただし、精度とのトレ