旗・膜のFSI解析 — トラブルシューティングガイド
より充実した内容を flag-membrane-fsi.html でご覧いただけます。
メッシュ変形による計算破綻
旗の変位が大きくなるとメッシュが潰れて止まります。
ALE法の限界だ。対策は以下の通り。
1. リメッシング頻度の増加: Ansys Fluentではremeshing intervalを小さく設定
2. overset meshの利用: STAR-CCM+やAnsys Fluentのoverlapping grid
3. IB法への切り替え: メッシュ固定で大変形に対応
4. スプリングアナロジーの強化: メッシュスムージングのばね定数を調整
膜の初期形状問題
初期状態で膜が自重で垂れ下がる状態をどう設定しますか？
prestress解析が必要だ。まず構造単体で自重+初期張力による静的解析を実施し、変形後の形状を初期形状としてFSI解析に引き継ぐ。AbaqusではIMPORT機能、AnsysではICSTATEで実現できる。
数値的な膜のめくれ
膜が非物理的に折れ曲がったり、自己交差したりすることがあります。
膜要素の自己接触判定が不足している場合に起きる。対策は、
- 自己接触条件（self-contact）を追加する
- 膜の曲げ剛性を微小でも設定する（純粋な膜要素ではなくシェル要素を使う）
- 安定化項（Baumgarte stabilization等）を追加する
LS-DYNAの場合、*CONTACT_AUTOMATIC_SINGLE_SURFACEで自己接触を設定できる。
「旗が反対方向にしかなびかない」——初期条件の非対称が招く計算バグ
フラッグFSIの計算で「旗が一方向にしか振れず、実際の翻翻りが再現できない」という症状で詰まる人が多いです。原因はほとんどの場合、初期条件の対称性が完璧すぎることです。旗を完全に平らな状態（変位ゼロ）から計算を始めると、流れも完全に対称になり、旗が上下どちらに振れるかの「決め手」がなくなります。対称系は微小な数値ノイズで崩れるはずですが、計算機の丸め誤差では十分な非対称性が得られないことがある。対策は「旗に0.01L（旗長の1%）程度の初期たわみを与える」こと。たったこれだけで計算が一気に安定します。Wind tunnel実験でも旗は「僅かな非対称」から翻翻りが始まるので、物理的にも正しい扱いです。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——旗・膜のFSI解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告