強連成と弱連成
考え方：情報をどう、いつ受け渡すか
概要
流体と構造の連成解析で、「弱連成で十分か、強連成にすべきか」と聞かれたんですが、判断の基準がわかりません。
まず2つの軸を分けて考えよう。1つ目は情報の向き。流体→構造だけの一方向か、構造の変形が流体にも返る双方向か。2つ目は受け渡しのタイミング。1ステップに1回ずつ交互に解いて相手の結果を1ステップ遅れで使うのが弱連成、同じステップ内で両者が一致するまで反復するか、まとめて解くのが強連成だよ。どちらが必要かは、相互作用がどれだけ強いかで決まる。流体構造連成では、それを表す分かりやすい指標がある。
連成の分類
| 分類 | やり方 | 向いている問題 |
|---|---|---|
| 一方向 | 片方を解いて結果をもう片方に渡すだけ | 変形が小さく流れに影響しない（熱応力、風荷重の強度評価） |
| 双方向・弱連成 | 1ステップに1回ずつ交互に解く（陽的） | 相互作用が弱い（空気中の翼のフラッタなど） |
| 双方向・強連成（分離型） | 各ステップ内で一致するまで反復 | 相互作用が強い（液体と軽い構造） |
| 双方向・強連成（一体型） | 全体を1つの連立方程式で解く | 相互作用が非常に強い、または反復が収束しにくい |
付加質量という目安
構造が加速すると、まわりの流体も一緒に動かす必要があるため、構造はあたかも質量が増えたように振る舞う。これが付加質量 $m_a$ で、構造の質量 $m_s$ との比 $r$ が相互作用の強さのよい目安になる。薄い殻の場合はおおよそ次の比になる（$\rho_f$・$\rho_s$ は流体・構造の密度、$R$ は代表寸法、$t$ は板厚）。
| 例 | 付加質量比の目安 | 連成の必要性 |
|---|---|---|
| 空気中の航空機の翼（アルミ、厚さ3 mm、代表寸法0.5 m） | 0.07 | 弱連成で可 |
| 血液の中の動脈壁（厚さ0.5 mm、半径4 mm） | 7.7 | 強連成が必要 |
血流シミュレーションが教えた教訓
航空機のフラッタ解析では、流体と構造を交互に解く弱連成が長く使われ、うまくいっていた。ところが2000年代に入って血管の中の血流と血管壁の変形を同じ方法で解こうとすると、時間刻みを小さくするほど早く発散するという奇妙な現象に直面した。原因は、血液と血管壁の密度がほぼ同じで壁が薄いため、付加質量が構造の質量を上回ることにあった。この「付加質量不安定」の解明をきっかけに、分離型の強連成の手法や緩和法の研究が大きく進んだ。
計算例：1自由度モデル
弱連成（陽的）が発散する条件
弱連成では、構造は1ステップ前の加速度から計算した流体力を受け取る。すると加速度の誤差が1ステップごとに $-r$ 倍になり、$r > 1$ では時間刻みによらず発散する。質量1・ばね定数1の構造で200ステップ計算し、加速度の振幅が何倍になったかを調べた。
| 付加質量比 r | Δt = 0.2 | Δt = 0.05 | Δt = 0.01 |
|---|---|---|---|
| 0.9 | 0.28倍（安定） | 0.54倍（安定） | 0.12倍（安定） |
| 1.1 | 約670万倍 | 約6.4万倍 | 約2,300倍 |
| 2.0 | 10⁵⁸倍 | 10⁵⁶倍 | 10⁵⁵倍 |
時間刻みを小さくしても、r が1を超えると結局発散するんですね。
そこが付加質量不安定のやっかいなところ。普通の数値不安定なら時間刻みを小さくすれば止まるけど、これは連成のやり方そのものが原因だから、刻みをいくら小さくしても治らない。表は同じ200ステップでの比較だけど、誤差は1ステップごとに増えるので、刻みを小さくすると同じ物理時間の中のステップ数が増えて、時間で見るとかえって早く発散する。現実の問題では、さっきの表のように密度比と板厚からおおよその r を見積もって、1に近いか超えるなら最初から強連成にするのが安全だよ。
強連成の反復と緩和
強連成では同じステップ内で加速度を反復して更新する。そのまま置き換えると r > 1 ではやはり発散するので、更新量に係数 $\omega$ をかける（緩和）。1自由度モデルで、相対誤差 $10^{-8}$ まで収束する反復回数を比べた。
| 付加質量比 r | 緩和なし | 固定 ω = 0.5 | Aitken |
|---|---|---|---|
| 0.2 | 12 | 21 | 2 |
| 1.0 | 収束せず | 1 | 1 |
| 2.0 | 収束せず | 29 | 2 |
| 10 | 収束せず | 収束せず | 2 |
固定の係数は、ちょうどよい値（この例では $1/(1+r)$）から外れると遅くなったり発散したりする。Aitken 法は直前2回の残差から係数を自動で決めるので、r が大きくても速く収束する。1自由度の線形問題なので2回で収束しているが、実際の問題では数回〜十数回の反復が必要になることが多い。
実務での選び方
判断の手順
- 構造の変形が流れ（または温度・電磁場）に効くかを考える。効かなければ一方向連成で十分。
- 双方向なら、付加質量比などの指標で相互作用の強さを見積もる。
- 弱い（r が十分小さい）なら弱連成で始め、時間刻みを変えて結果が変わらないことを確かめる。
- 強い、または弱連成で振動・発散するなら、サブ反復付きの強連成（Aitken 緩和など）に切り替える。収束しにくければ一体型を検討する。
流体構造連成以外の例
- 熱と構造：変形で接触状態や隙間が変わり、熱の伝わり方が変わる場合は双方向。熱膨張が小さく接触が変わらなければ一方向で十分。
- 電磁と熱：温度で導電率や透磁率が大きく変わる誘導加熱では、双方向の連成が必要。
- 時間スケールの違い：熱は遅く電磁は速いなど、時間スケールが大きく違う場合は、速い側を周期平均して遅い側に渡す方法がよく使われる。
「連成したのに結果が変わらない」は良い知らせ
双方向の強連成で解き直したら、一方向で解いた結果とほとんど変わらなかった――これは時間の無駄ではなく、一方向で十分だったという証拠になる。連成解析では、まず簡単な方法で解き、より厳密な方法で結果が変わるかを確かめる、という順番が効率的だ。結果が変わらなければ、次からは簡単な方法を使ってよい根拠ができる。逆に大きく変われば、どの効果が効いているかを理解する手がかりになる。
よくある間違い
間違いと対策
| 間違い | 起きること | 対策 |
|---|---|---|
| 発散したら時間刻みを小さくする | 付加質量不安定は治らない | 付加質量比を見積もり強連成へ |
| 緩和係数を固定のまま | 反復が遅い・発散する | Aitken など適応的な緩和 |
| サブ反復の収束判定がゆるい | エネルギーが増えて不安定 | 界面の変位・力の両方で判定 |
| メッシュ間の受け渡しを確かめない | 力の合計が保存されない | 受け渡し前後で合力を比較 |
| 最初から最も厳密な方法で解く | 計算時間がかかりすぎる | 簡単な方法から段階的に |
反復や緩和の考え方を試してみたいです。
ガウス・ザイデル法シミュレーターとSOR緩和法シミュレーターで緩和係数による収束の速さの違い、ニュートン・ラフソン法シミュレーターで非線形反復の収束を確かめられるよ。連成の方式の詳しい比較は分離型と一体型と交互（スタッガード）連成、メッシュ間のデータの受け渡しは異なるメッシュ間のデータ受け渡しの記事にまとめてある。
関連トピック
なった
詳しく
報告