橋梁の風荷重FSI — トラブルシューティングガイド
この記事は統合版に移行しました
より充実した内容を bridge-wind.html でご覧いただけます。
より充実した内容を bridge-wind.html でご覧いただけます。
よくある問題と対策
橋梁断面のCFD解析で困ることは何ですか？
順に見ていこう。
1. Strouhal数が実験と合わない
症状: CFDの渦放出振動数が風洞試験値と10%以上ずれる。
原因: メッシュ分解能不足、乱流モデルの不適切な選択、計算領域が狭い。
対策:
- 桁面付近のメッシュを細分化（最低でも桁厚の1/20以下の要素サイズ）
- LESまたはDESを使用（定常RANS k-εは渦放出を再現できない）
- 下流領域を20B以上に延長
2. フラッタ微係数の散らばりが大きい
症状: 換算風速ごとのフラッタ微係数が振動的で、滑らかな曲線にならない。
原因: 計算時間が短く、空力力の統計的安定が不十分。
対策:
- 各換算風速で最低20振動サイクル以上の計算を実行
- 初期過渡を除去（最初の5サイクルは除外）
- 時間刻みをCFL < 1に設定
3. 自由振動CFDでロックインが再現されない
渦励振のシミュレーションでロックインが出ないのはなぜですか？
原因:
- 構造の質量比（$m^* = m/(\rho_f D^2)$）が大きすぎる設定ミス
- 構造減衰が過大（Scruton数 $Sc = 2m^*\zeta$ が臨界値を超えている）
- 2D計算では3D効果による相関長の影響が無視される
対策:
- 質量比と減衰比を実橋の値と照合
- Scruton数の閾値（一般に $Sc < 5$〜$10$ でVIVが発生）を確認
- 3D計算でスパン方向の位相ずれを考慮
4. バフェティング解析の結果が過大
対策:
- 入口乱流の空間相関長が適切か確認（Davenportの相関係数）
- 空力アドミッタンス関数の取り方を再検討
- モーダル解析のモード数を十分に含める
| 検証項目 | 風洞試験値との目標誤差 |
|---|---|
| $C_D$（抗力係数） | ±10% |
| Strouhal数 | ±5% |
| フラッタ微係数 | ±15% |
| 臨界フラッタ風速 | ±10% |
CFD単独では信頼性が不十分で、風洞試験との照合が不可欠なんですね。
Coffee Break よもやま話
橋梁FSIで「発散しない渦」が出る——LESの落とし穴
橋梁風荷重のFSI解析でLES（ラージ・エディ・シミュレーション）を使うと「渦が出ない、あるいは実験と全く形が違う」という悩みが出ます。多くの場合、原因は入口境界条件の乱流強度不足です。風洞実験の大気境界層は地表面の粗さで自然に乱れていますが、CFDの入口に一様流を与えると下流で整流されてしまい、カルマン渦が弱くなります。合成乱流（SEM法やDFSEM）で現実的な乱流境界層を入口に与えることが改善の鍵です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——橋梁の風荷重FSIの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告