Time-Domain Two-Way Fluid–Structure Interaction Analysis of the HIRENASD Transonic Wing: Static Aeroelastic Response and Flutter Prediction

技術要約：HIRENASD 遷音速翼における時間領域二方向流体–構造相互作用解析

問題提起
輸送機翼の構造的完全性は、遷音速領域における複雑なエアロエラスティック（空気弾性）結合、特に破滅的な破壊につながる可能性のあるフラッターによって支配される。Doublet Lattice Method (DLM) のような線形周波数領域法は、移動衝撃波、衝撃波誘起境界層剥離、および強い流体–構造相互作用を含む、遷音速流に固有の重要な非線形現象を捉えるための忠実度が不足していることが多い。High Reynolds Number Aero-Structural Dynamics (HIRENASD) プロジェクトは、超臨界翼（BAC3-11 エアフォイル）に関する包括的な実験データベースを提供しているが、市販の汎用 CFD-CSD ソルバーを用いたこれらの構成の検証に関するオープンな文献には、依然として大きな空白が存在する。さらに、複数のマッハ数にわたるフラッター境界の系統的な時間領域特定、および自由流密度とマッハ数がフラッター動圧にどのように共同で影響するかという定量的な分析は、このベンチマークに対して完全には取り組まれていない。

手法
本研究では、ANSYS Workbench 内に実装された、時間領域、疎結合（loosely-coupled）の二方向流体–構造相互作用（FSI）フレームワークを採用している。アプローチは、以下の要素間で解を分割する：

• CFD: 密度ベースのソルバーを用いた ANSYS Fluent を使用し、SST $k-\omega$ 乱流モデルを用いたレイノルズ平均 ナビエ–ストークス（RANS）方程式を利用する。流体ドメインは、$y^+ \approx 1$ を達成するためにプリズム層を備えた非構造格子テトラメッシュ（約60万セル）を特徴とする。
• CSD: HIRENASD 翼構造（ニッケル・マレージング鋼）の六面体有限要素モデル（41,923 要素）を用いた ANSYS Mechanical を使用する。モード解析には Block-Lancos アルゴリズムを用い、過渡積分には Newmark-$\beta$ 法を用いる。
• 結合: ANSYS System Coupling が、力の保存のための General Grid Interface (GGI) および変位マッピングのための Smart Bucket インターポレーションを介したデータ交換を管理する。メッシュ変形は、局所的な再メッシュを伴うスプリング・アナロジー・スムージングによって処理される。

本研究は、NASA FUN3D の空気力学的ベンチマークおよび HIRENASD 実験による構造周波数に対してフレームワークを検証する。静的エアロエラスティック解析は、迎角 $-1.5^\circ$ から $+4.5^\circ$ の範囲において、マッハ数 ($M$) 0.7 および 0.8 で実施される。動的フラッター解析は、2 つの異なる飛行条件（$M=0.8$ ($\rho = 0.2525$ kg/m³) および $M=0.7$ ($\rho = 0.50$ kg/m³)）において、自由流速度を系統的に変化させることで、亜臨界（減衰）、臨界（中立）、および超臨界（発散）の応答レジームを特徴付け、フラッター境界を特定する。

主な貢献

1. フレームワークの検証: 本研究は、排他的に商用ツール（ANSYS Fluent, Mechanical, System Coupling）のみを使用して、遷音速エアロエラスティシティを実現する、再現可能で産業的にアクセス可能な FSI ワークフローを確立することに成功した。
2. 静的エアロエラスティック・ベンチマーキング: 本フレームワークは、定常的な空気力学的係数 ($C_L, C_D, C_M$) を NASA FUN3D 参照データに対して 1.87% から 3.52% 以内の誤差で再現し、基本曲げ周波数 (26.504 Hz) を実験による地上振動試験結果の 1.94% 以内で予測する。
3. 動的フラッター特性の決定: 本研究は、減衰振動から発散振動への遷移を明示的に特徴付けることにより、時間領域による系統的なフラッター境界の決定を提供する。また、テストされた両方のマッハ数レジームにおいて、フラッターが単一モードの曲げメカニズムによって支配されていることを特定した。
4. スケーリングに関する物理的洞察: 本研究は、マッハ数、自由流密度、およびフラッター動圧の関係を定量化し、フラッター速度はより低い密度（高高度）において高くなる可能性がある一方で、フラッター動圧は $q_f \propto \rho V_f^2$ というスケーリング則により、より高い密度において著しく高くなることを示した。

結果

• 静的応答: $M=0.8$ において、揚力傾斜は約 0.107/degree である。最大翼端たわみは、迎角 $4.5^\circ$ において 12 mm（無次元化 $\delta/b = 0.0093$）に達する。$M=0.8$ の抗力極曲線は鋭い衝撃波誘起抗力上昇を示し、一方 $M=0.7$ の極曲線はより放物線状であり、誘導抗力が支配的である。
• モード解析: 第1次面外曲げモードは 26.504 Hz と予測された（実験値に対し 1.94% の誤差）。高次モードは、有限要素モデルにおける構造的簡略化（例：ファスナーのコンプライアンスの欠如）に起因する、系統的な過大評価（9–13%）を示している。
• フラッター境界:
  • $M=0.8$ ($\rho = 0.2525$ kg/m³) において: フラッター速度 ($V_f$) は 266.56 m/s、フラッター動圧 ($q_f$) は 8.97 kPa である。フラッター周波数は約 22 Hz であり、これは遷音速的な空気力学的剛性の修正による、真空中の構造周波数からの 16.9% の減少を表している。
  • $M=0.7$ ($\rho = 0.50$ kg/m³) において: フラッター速度は 230.46 m/s、フラッター動圧は 13.27 kPa である。フラッター周波数は約 25.5 Hz であり、構造周波数からの減少は 3.8% と小さく、遷音速の非線形性が弱いことを示している。
• スケーリング則: $M=0.7$ におけるフラッター速度は $M=0.8$ よりも低いにもかかわらず、自由流密度が 98% 高いため、フラッター動圧は 47.9% 高くなっている。これは、次元的なスケーリング $q_f \propto \rho V_f^2$ を裏付けている。

意義および主張
本論文は、検証された ANSYS Workbench フレームワークが、独自のインハウスコードやプロプライエタリなソルバーを必要とせずに、輸送機カテゴリーの翼の予備的なフラッター境界評価のための、実用的かつ産業に関連する経路を提供すると主張している。研究は、商用ツールのみを使用することが、産業デザイン環境への即時の転用を可能にすることを強調している。

著者らは、フラッター動圧の密度依存的なスケーリングの定量的確認が、風洞モデルの結果をフルスケールの飛行条件へ外挿するための、科学的に厳密な基礎を提供すると断言している。特に、高度依存のフラッターリスクに関するものである。本研究は、一定のマッハ数であっても、低高度（高密度）ではフラッターが発生しないが、巡航高度（低密度）では有効な空気力学的剛性の減少により、フラッター境界に近づく可能性があることを浮き彫りにしている。研究は、特定されたフラッター速度は、AGARD 445.6 翼の検証済みフレームワークとの比較に基づき、5–8% の手法的な不確実性を伴うと結論付けており、観察された単一モード曲げメカニズムは、この特定の構成における曲げモードとねじれモードの間の大きな周波数分離と一致している。