渦励振（VIV）解析 — トラブルシューティングガイド
より充実した内容を vortex-induced-vibration.html でご覧いただけます。
連成計算が発散する場合
VIVの連成計算を回したら発散してしまいました。何が原因でしょうか？
VIV-FSI解析の発散は主に以下の原因で起きる。
| 症状 | 原因 | 対策 |
|---|---|---|
| 1ステップ目で発散 | 界面マッピングの不整合 | 流体・構造の界面メッシュを確認。ノードの対応関係を検証 |
| 数ステップ後に発散 | added mass不安定性 | 強連成に切り替え。緩和係数を0.1〜0.3に下げる |
| 振幅成長後に発散 | メッシュ変形過大 | リメッシング頻度を上げる。overset meshを検討 |
| 特定流速でのみ発散 | ロックイン領域での共振 | 時間刻みを細かくし、サブイテレーション回数を増やす |
added mass不安定性って具体的にはどういうことですか？
流体の付加質量が構造質量と同程度以上（$m^* \lesssim 1$）になると、弱連成では人工的なエネルギーが界面に蓄積して発散する。Causin et al.（2005）の解析で理論的に示されている。対策は強連成＋緩和、またはRobin-Neumann条件の適用だ。
Strouhal数が文献値と合わない
計算で得られたStrouhal数が0.2から大きくずれるんですが。
チェックポイントを整理しよう。
- ブロッケージ比: 計算領域幅に対する円柱直径の比が5%以上だと壁の影響でStが変わる。$D/H < 0.03$ を確保する
- 時間分解能: FFTのサンプリングが不足すると周波数の同定精度が落ちる。Welch法で窓関数を適用し、十分な時間データを使う
- 境界条件: 出口境界が近すぎると後流が拘束される。$30D$ 以上離す
- 2D vs 3D: Re > 200 では3D効果が顕著。2D計算はStを過大評価する傾向がある
疲労評価への接続
VIV解析の結果を疲労寿命評価に使うにはどうすればいいですか？
時刻歴応力をレインフロー法でカウントし、S-N曲線とMinerの線形累積損傷則で疲労損傷度を算出する。
DNV-RP-C203のS-N曲線を使い、$D_{fatigue} < 1/DFF$（DFF: Design Fatigue Factor、通常3〜10）を満足させる。応力集中係数（SCF）の設定が結果を大きく左右するから注意が必要だ。
CFD-FSIの結果をそのまま疲労評価に流せるパイプラインがあると便利ですね。
AnsysではFluent→Mechanical→nCode DesignLifeの連携が可能だ。STAR-CCM+もFE-Safeとの連携インターフェースを持っている。ただし、CFDの時刻歴は非常に長いデータになるから、サイクルカウントの効率化が実務上の課題だよ。
「設計範囲外のはずの風速で振動する」——VIVの見落としパターン
VIVのトラブルでよくあるのが「設計上はロック・インが起きないはずの流速で振動している」というケース。原因を調査すると、流体が構造周囲で加速しており「局所流速」が設計の「平均流速」より40〜50%高かったという事例がある。構造物が地面や壁面に近いとベルヌーイ効果で流速が高まり、また群集した構造物（複数の煙突が並ぶ場合）では上流の構造が下流側の流れパターンを激変させる。もう一つの見落としは「モード形状の変化」で、実際の境界条件が設計想定（ピン-ピン支持など）と異なり、固有振動数が大幅にずれているケースだ。VIVトラブルの際は「本当の局所流速」と「本当の固有振動数」の確認から始めるのが鉄則だ。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——渦励振（VIV）解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告