空力弾性解析 — トラブルシューティング
より充実した内容を aeroelasticity.html でご覧いただけます。
連成計算が発散する
CFD-CSD連成計算が発散してしまうんですが、何が原因ですか？
空力弾性連成の発散原因のトップ3を示そう。
| 原因 | 症状 | 対策 |
|---|---|---|
| 付加質量効果 | メッシュ変位が発散 | Under-relaxationを0.3-0.5に下げる |
| メッシュ変形破綻 | 負体積セル発生 | Laplacian smoothingにdiffusivity=inverseDist使用 |
| 時間刻み不整合 | 解が振動的に発散 | 流体・構造の時間刻みを一致させる |
付加質量効果って具体的にどういうことですか？
翼面が動くと、周囲の流体も連動して動く。この流体の慣性が「付加質量」として構造に作用する。密度比 $\rho_s / \rho_f$ が小さい（構造が軽い）ほど付加質量の影響が大きくなり、明示的な連成スキーム（1ステップで流体→構造→流体）では不安定になる。
対策は以下の通り。
1. 陰的連成: 各時間ステップ内でFSIサブイテレーションを回す（System Couplingのcoupling iterations）
2. Aitken加速: 動的にunder-relaxation factorを最適化する
3. Robin-Neumann法: 界面条件にRobin条件を導入して安定化
メッシュ変形の品質劣化
翼の変形量が大きいとき、メッシュが潰れてしまうんですが...
メッシュ変形手法の選択が鍵だ。
| 手法 | 大変形耐性 | 計算コスト | 実装ソフト |
|---|---|---|---|
| Spring analogy | 低〜中 | 低 | Fluent |
| Diffusion-based smoothing | 中〜高 | 中 | Fluent, STAR-CCM+ |
| RBF（Radial Basis Function） | 高 | 高 | OpenFOAM（RBFMotionSolver） |
| Overset mesh | 非常に高 | 高 | Fluent, STAR-CCM+ |
大変形が予想される場合はOverset mesh（Chimera法）が最も安全だ。翼まわりのメッシュを独立した体として移動・回転させ、背景メッシュとの重なり領域で補間する。メッシュの変形は不要だから品質劣化がない。
フラッター速度がモデル間で大きく異なる
DLMとCFDでフラッター速度が20%以上違うんですが、どちらを信じるべきですか？
遷音速域で20%の差は珍しくない。DLMはポテンシャル流れの仮定だから、衝撃波の効果が正確に反映されない。一方、CFDの結果も乱流モデルやメッシュ密度に依存する。
判断基準は以下の通りだ。
1. まず、CFDの定常圧力分布が実験データ（もしあれば）と一致するか確認
2. CFDのフラッター結果が格子収束していることを確認（2水準以上のメッシュ）
3. 遷音速域ではCFD結果をベースに判断するが、DLMとの差分の原因を理解しておく
4. 認証ではDLMを基本とし、CFDを補正データとして活用するのが業界標準
最終的な認証はDLMベースなんですね。
そう。ただしFAA/EASAもCFDベースのフラッター解析を認証に使うためのガイドライン整備を進めている。将来的にはCFDが主役になる可能性は高いよ。
フラッター解析が「収束したのに答えが正反対」になる落とし穴
空力弾性解析のトラブルで意外と多いのが、固有モードの符号や位相を間違えたまま計算を進めてしまうケースだ。フラッター速度の計算では、流体と構造の連成行列の固有値虚部の符号が重要で、符号を誤ると「安定」と「不安定」が逆転してしまう。「解析では余裕があると出たのに実機試験でフラッターが発生した」という事態は、担当者にとって悪夢そのもの。連成解析は符号・位相の確認が命です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——空力弾性解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告