一方向流体構造連成

一方向FSIでは、流体解析で得られた圧力・せん断力を構造解析の荷重条件として転送するだけで、構造の変形が流体に影響を与えないと仮定する。つまり情報の流れが「流体→構造」の一方通行だ。

この仮定が妥当なのは、構造の変形が流体場を有意に変えない場合だ。具体的には、

以下のような場面では一方向FSIで十分な精度が得られる。

圧倒的に安い。CFDとFEAを各1回ずつ解くだけだから、双方向FSI（反復的に何十回も解く）と比べて計算時間は1/10-1/100のオーダーだ。設計初期のスクリーニングや、構造強度の概算には非常に効率的だよ。

CFDメッシュ（通常は細かい）からFEAメッシュ（通常は粗い）へのデータマッピングが必要だ。主な手法は以下の通り。

力（圧力・せん断力）の転送ではConservative mappingが重要だ。これは、CFD面上の力の合計がFEA面上の力の合計と厳密に一致することを保証する。Ansys System CouplingやSTAR-CCM+のCo-Simulationではこのオプションが用意されている。

例えば揚力100 NのCFD結果をFEAに転送したとき、FEA上での積分値が95 Nや105 Nになるということ。構造の全体的な力の釣り合いが崩れるから、変形や反力の結果に誤差が入る。特に全体的な荷重が設計の判断基準となる場合（例：ブレードの根元応力）は保存性のあるマッピングが必須だ。

2つのアプローチがある。

1. 瞬時値転送: CFDの各時間ステップの圧力分布をFEAの過渡解析に入力。時刻歴の全てのステップでマッピングを行う。精度は高いが、データ量が膨大になる。

ランダム振動解析（Nastran SOL 111/112）で、圧力PSDを入力として応力のPSDを計算し、Dirlik法やSteinberg法で疲労寿命を推算する。この手法はロケットのフェアリング音響荷重解析などで広く使われている。

一方向FSIならファイルベースでも十分実用的だ。

ファイルベースの場合、FluentからCGNS/EnSight形式で壁面圧力を出力し、Ansys MechanicalやAbaqusにImported Pressureとして読み込む。マッピングの精度はソフトのImport機能に依存するから、転送後の力の総和を確認することが大事だよ。

Workbenchのプロジェクトスキマティックで以下のように構成する。

1. Fluent Analysis（System A）を配置し、CFD解析を完了

Workbenchの強みだね。ただし注意点がある。

STAR-CCM+のExport機能で壁面データをAbaqus INP形式に出力する方法がある。ステップは以下の通り。

非定常の場合はステップごとにINPファイルを出力して、AbaqusのAmplitude tabular dataとして時刻歴を定義する。あるいはCo-Simulation Engineを使って自動的に時刻同期させる方法もあるよ。

簡単なチェック方法がある。一方向FSIで得られた構造の変位 $\delta$ を確認して、

もう一つの指標は無次元変位パラメータ $\Pi$ だ。

その通り。一方向FSIは計算コストが低いから、まずスクリーニングとして使い、双方向FSIの必要性を判断するのが効率的だ。

その通り。CHT（Conjugate Heat Transfer）+ 熱応力解析は一方向FSIの典型例だ。ワークフローは以下の通り。

FluentでCHTを解いて固体領域の温度分布を取得し、それをAnsys Mechanicalに転送して熱応力解析を行う。Workbenchなら温度データの転送が自動化されるよ。

もちろん。Workbenchで同じFluent解析から温度と圧力の両方をMechanicalに転送できる。Imported PressureとImported Body Temperatureの2つのロードがセットアップされる。ターボ機械のブレード解析では、遠心力+空力+熱応力の3つの荷重を同時に考慮することが標準だ。

いくつかの代表例を示そう。

走行風によるドアミラーの振動は運転者の視界に直結するから、自動車メーカーでは重要な評価項目だ。CFDで圧力変動のPSDを取得し、FEAのランダム振動解析で振幅を予測する。目標はミラーの固有振動数を風圧変動の卓越周波数から離すことだよ。

Ansys PyFluentとAnsys PyMechanicalを使えば完全にPythonで自動化できる。

大まかな流れはこうだ。

流速を5条件、迎角を3条件、合計15ケースといったパラメータスタディをforループで回せる。各ケースの最大応力と変位をCSVに集約してPandas/Matplotlibで可視化すれば、設計空間の探索が効率的に行えるよ。

その通り。ロケットの打ち上げ時に発生する高強度音響環境（約140-160 dB）がフェアリング内部の搭載機器に振動荷重を与える。これを音響-構造連成（Vibroacoustics）として解析する。

ワークフローは以下の通りだ。

近年はLESで排気ジェットの近傍場音源を解析し、FWH面積分で遠方場音圧を計算する手法も使われ始めている。ただし周波数範囲が数kHzに及ぶため、メッシュ解像度の要求が非常に厳しい。実務では経験的な音響パワースペクトルを使う方が効率的な場合が多いよ。

建築基準法や各国の設計基準（Eurocode 1、ASCE 7等）では風荷重係数が規定されているけど、特殊形状の建物や周辺建物との干渉効果はCFDで評価することが増えている。

一方向FSIの流れはこうなる。

正確には、CFDの結果を風洞試験の代替として使うには認証が必要だ。多くの国では風洞試験が依然として基準だけど、CTBUH（Council on Tall Buildings and Urban Habitat）のガイドラインでは、CFDと風洞試験の比較検証を条件にCFD結果の設計利用を認めつつある。

Causin et al. (2005) の研究に基づく判断基準を示そう。

その通り。海洋工学のライザー管やスパーブイ、橋梁のケーブルなどは密度比が小さく、VIV（渦励起振動）のLock-in現象が問題になるから双方向FSIが必須だ。一方、航空機の外板パネルや建築物は空気中の構造だから一方向で十分なケースが多い。

よくある原因と対策を整理しよう。

FluentとMechanicalのWorkbench連携なら座標系は自動的に一致するけど、異なるソフト間（STAR-CCM+→Abaqus等）ではCADのImport設定で座標系を揃える必要がある。単位系の違い（m vs mm）にも注意だ。

以下のフィルタリング手法が有効だ。

構造の関心のある固有振動数 $f_n$ に対して、$f_{cutoff} = 5-10 \times f_n$ 程度にする。例えば $f_n = 100$ Hz の構造なら $f_{cutoff} = 500-1000$ Hz。これより高い周波数の圧力変動は構造をほとんど励起しないから、除去しても結果に影響しない。

構造モデルの拘束条件を確認しよう。CFDの圧力は絶対圧（ゲージ圧 + 大気圧）か、ゲージ圧かを確認する。閉じた構造（パイプなど）の内面にゲージ圧を載せる場合、外面に大気圧を載せ忘れると差圧が過大になる。

チェックポイント:

その通り。Operating Pressureが0ならゲージ圧=絶対圧。Operating Pressureが101325 Paなら壁面出力はゲージ圧（大気圧基準）だ。構造解析に渡すときは、この設定を把握しておくことが重要だよ。