空力弾性解析 — LCO・非線形・縮約モデル

それが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倍以上速くなる。

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

1. 非定常空気力のサロゲートモデル: マッハ数・迎角・振動振幅をパラメータとして、CFD空気力を即座に予測するDNNモデル

2. フラッター境界の予測: 飛行条件（$M$, $q_\infty$, 高度）から直接フラッターマージンを出力する分類器

3. GVTデータからのモデルアップデート: ベイズ推定によるFEMモデルパラメータの同定

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