氷-構造相互作用 — トラブルシューティングガイド
この記事は統合版に移行しました
より充実した内容を ice-structure-interaction.html でご覧いただけます。
より充実した内容を ice-structure-interaction.html でご覧いただけます。
氷-構造相互作用 — トラブルシューティングガイド
Element Erosionによるエネルギー散逸
要素消去で氷を削除するとエネルギーバランスがおかしくなるんですが。
要素消去時に残留内部エネルギーが系から失われる。この影響を確認するために、
- エネルギーバランス（全エネルギー保存のモニタリング）
- 消去された要素のエネルギー総和を記録
- eroded internal energyが全エネルギーの5%以下であることを確認
LS-DYNAではglstatファイルでエネルギーバランスを確認できる。
接触不安定
氷と構造の接触面でノードが貫通してしまいます。
ペナルティ法のスケーリングファクターが不足している。対策は、
| 問題 | 対策 |
|---|---|
| ノード貫通 | SLSFAC（ペナルティスケーリング）を増加。0.1→1.0 |
| 接触振動 | 粘性減衰を追加。VDCパラメータ |
| 不整合なメッシュ密度 | 接触面のメッシュサイズを揃える |
| 非物理的な接触力 | *CONTACT_FORCE_TRANSDUCER_PENALTYで力を確認 |
結果の妥当性はどう検証するんですか？
ISO 19906の経験式による氷荷重推定値と比較する。また、Sodhi & Hatten（1980）やTimcoの実験データベースとの比較も有用だ。接触面積-圧力関係のスケール効果を考慮することが重要だ。
Coffee Break よもやま話
「氷が壊れずに船を突き抜けた」——破壊基準の設定ミスで起きる計算のバグ
氷-構造衝突のFEM解析でよく起きるバグが「氷が破壊されずに船体を貫通してしまう」現象です。原因はほぼ100%「氷の破壊基準の閾値が高すぎる」か「要素削除の設定を忘れている」かのどちらかです。氷は脆性材料なので、引張応力が破壊強度（一般的な海氷で1〜2MPa）を超えたら要素を削除するエロージョン（element erosion）を使うのが標準的な対処法です。ただし要素削除しすぎると「質量が消える」ため、エネルギー収支がおかしくなる副作用があります。経験則では、削除される質量の合計が初期モデルの5%を超えたら結果を信頼してはいけません。また、ひずみ速度依存の破壊強度を使わないと「低速衝突で氷が割れない」という非物理的な結果になるので、速度依存材料モデルの設定は必須です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——氷-構造相互作用の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
この記事の評価
ご回答ありがとうございます！
参考に
なった
なった
もっと
詳しく
詳しく
誤りを
報告
報告