氷-構造相互作用
氷-構造相互作用の理論基礎
現象の概要
氷と構造物の相互作用って、どういう問題ですか？
北極海の石油プラットフォームや砕氷船が海氷と衝突する際の荷重評価だ。氷は破砕・座屈・クリープなど複雑な破壊挙動を示すため、単純な荷重モデルでは不十分な場合がある。
支配方程式
氷の力学モデルにはどんなものがありますか？
氷は粘弾塑性体としてモデル化されることが多い。Sinha（1978）のモデルでは、
弾性ひずみ、遅延弾性ひずみ（粒界すべり）、粘性クリープひずみの3成分に分解する。高ひずみ速度（衝突時）ではMaxwell-Brittle型モデルも使われる。
構造側は通常の弾塑性FEMだ。氷-構造の接触問題としてペナルティ法やAugmented Lagrangian法で界面を処理する。氷の破砕はelement erosion（要素消去）やCZM（Cohesive Zone Model）で表現する。
氷荷重の経験式もあるんですか？
ISO 19906（Arctic offshore structures）では、接触面積 $A$ に対する氷圧力を、
で与えている。$C_R$はIce Reference Strength、$h$は氷厚だ。数値シミュレーションの検証にこの経験式との比較が有用だ。
氷は「金属より複雑」——温度・速度・塩分で変わる氷の力学特性
構造解析で氷を扱うとき、最初に直面する壁が「氷の材料モデルは何を使えばいいのか」という問題です。氷は温度が-2℃と-20℃では強度が倍以上違い、ひずみ速度が遅いとクリープ（粘性流動）し、速いと脆性破壊します。さらに海氷は塩分を含むため純氷より弱く、同じ「氷」でも特性が大幅に異なります。工学的には「ひずみ速度で破壊モードが遷移する」境界が約10⁻³/s程度とされており、船の衝突（高速）では脆性破壊、橋脚への氷圧（低速）ではクリープが支配的です。この二つの破壊モードを一つの材料モデルで表現することが氷構造相互作用シミュレーションの理論的な核心になっています。
氷-構造相互作用の数値計算手法
離散化手法
氷の破砕をどうやって数値的に扱うんですか？
主に3つのアプローチがある。
| 手法 | 特徴 | 適用 |
|---|---|---|
| FEM + Element Erosion | 破壊基準到達で要素削除 | LS-DYNA, Abaqus/Explicit |
| DEM（離散要素法） | 粒子の集合体として氷を表現 | PFC, YADE |
| SPH | メッシュフリー。破砕の追跡が容易 | LS-DYNA SPH |
| Peridynamics | 非局所モデル。亀裂が自然に発生 | Peridigm |
DEMで氷を扱うのは面白いですね。
DEMでは氷を多数の円盤（2D）や球体（3D）の集合体で表現し、粒子間のボンドが破壊基準を超えると切断される。氷の破砕パターン（radial crack, circumferential crack）を自然に再現できる利点がある。
接触アルゴリズム
氷と構造物の接触はどう処理するんですか？
LS-DYNAでは*CONTACT_AUTOMATIC_SURFACE_TO_SURFACEで摩擦接触を定義する。氷の摩擦係数は温度に強く依存し、$\mu = 0.01$〜$0.3$ の範囲で変動する。
陽解法の時間刻みは最小要素サイズと音速で決まるCourant条件に従う。
氷中の音速は約3,000 m/sなので、要素サイズ0.01 mでは $\Delta t \approx 3 \times 10^{-6}$ s と非常に短くなる。
DEM（個別要素法）で「氷の群れ」を計算する——FEMでは無理な大規模砕氷
砕氷船の実航行では、砕かれた無数の氷片が船体周囲を流れ続けます。この「砕氷片の群れ」をFEM（有限要素法）で一つ一つモデル化すると、要素数が数百万〜数千万になり計算不可能です。そこで使われるのがDEM（Discrete Element Method：個別要素法）です。DEMでは各氷片を「剛体+バネ+ダッシュポット」の単純モデルで表現し、氷片同士・氷片と船体の接触力を効率的に計算します。ノルウェーのSINTEF研究所では100万個以上の氷片を含むシミュレーションを実施し、砕氷船の抵抗力を実測と±15%以内で一致させています。DEMは氷以外にも、砂・岩・粒状材料の解析に広く使われており、氷-構造問題での発展が粒状体力学全体の技術向上にも貢献しています。
氷-構造相互作用の実務適用
モデル構築の手順
氷-構造相互作用のシミュレーションを始める手順を教えてください。
1. 構造物の3D FEモデル作成（鋼構造の場合、シェル要素）
2. 氷板のモデル作成（ソリッド要素。破壊を扱うならerosion設定を追加）
3. 接触定義（面対面接触、摩擦係数の設定）
4. 氷の初期速度・ドリフト速度の設定
5. 材料モデルの定義（氷：Tsai-Wu破壊基準等、構造：弾塑性）
6. 陽解法で実行
氷の材料パラメータ
氷の物性値はどうやって決めるんですか？
海氷の物性は温度、塩分濃度、ひずみ速度に強く依存する。
| パラメータ | 1年氷（-10°C） | 多年氷（-10°C） |
|---|---|---|
| Young率 | 3〜9 GPa | 5〜10 GPa |
| 圧縮強度 | 2〜10 MPa | 5〜15 MPa |
| 引張強度 | 0.5〜2 MPa | 1〜3 MPa |
| Poisson比 | 0.33 | 0.33 |
| 密度 | 900 kg/m³ | 910 kg/m³ |
物性のばらつきが大きいですね。
そのためパラメトリックスタディが必須だ。ISO 19906では特性値として50年再現期間の値を使うことが規定されている。モンテカルロシミュレーションで確率論的評価を行う場合もある。
砕氷船の設計——「氷を割るプロセス」をシミュレーションで最適化する
砕氷船は単純に「ぶつかって割る」のではなく、船首の傾斜角によって氷に「曲げ破壊」を起こすことで効率よく砕氷します。傾斜角が浅すぎると氷が砕けずに船体下に潜り込み、急すぎると氷を圧縮方向に押してしまい非常に高い荷重が生まれます。最適な傾斜角は氷の厚さと強度によって変わるため、実務では「厚さ1.5m、圧縮強度2MPaの海氷を連続砕氷する」という設計条件でシミュレーションを回し、船首形状と推進力の組み合わせを最適化します。ロシアのアトミック号（原子力砕氷船）の設計では、砕氷シミュレーションの計算結果が船首形状の変更に直結し、砕氷能力が初期設計より20%向上したと報告されています。
氷-構造相互作用のソフトウェア比較
ツール比較
氷-構造シミュレーションに使えるソフトウェアは？
LS-DYNAの氷モデルはどう設定するんですか？
*MAT_ISOTROPIC_ELASTIC_FAILURE（MAT_13）が基本だ。最大主応力または最大ひずみで破壊判定し、element erosionで削除する。より高度にはMAT_JOHNSON_HOLMQUIST_CERAMICS（MAT_110）で圧力依存の破壊挙動を表現できる。
OSSの選択肢はありますか？
YADE（DEM）やPeridgm（Peridynamics）が研究用途で使われている。しかし氷の材料モデルの実装は限定的なので、ユーザーサブルーチンの開発が必要になることが多い。
LS-DYNAが氷-構造解析で強い理由——超高速衝突と破壊を一括処理
氷-構造相互作用のシミュレーションで商用ツールとして最も使われているのがLS-DYNAです。その理由は「衝突・破壊・接触の同時処理」が得意なため。氷が船体に衝突する際には、氷の脆性破壊、飛散片の接触、海水中の氷片の浮力——複数の現象が同時に起きます。LS-DYNAは陽解法（explicit）で時間積分するため、衝突のような短時間大変形問題に強く、氷の破壊をSPH（粒子法）で表現し砕氷片の飛散まで追跡できます。一方AbaqusはABAQUS/Explicitで同様の計算ができますが、氷の破壊モデルのライブラリはLS-DYNAの方が充実しています。なお、オープンソースではOpenFOAMの固体力学拡張版（solids4foam）で氷の弾塑性解析もできますが、破壊表現は研究レベルです。
氷-構造相互作用の先端研究
Peridynamicsによる氷の破壊
Peridynamicsって従来のFEMと何が違うんですか？
古典的な連続体力学の偏微分方程式ではなく、積分方程式で定式化する非局所モデルだ。変位の空間微分が不要なため、亀裂が自然に発生・進展する。
$H_\delta$ はhorizon（影響範囲）、$\mathbf{f}$ はbond forceだ。bondのcritical stretchを超えると破断する。
氷海航行シミュレーション
砕氷船が氷を割りながら進む様子もシミュレーションできるんですか？
DEM-FEM連成で砕氷船の抵抗と氷片の飛散パターンを予測する研究が進んでいる。氷をDEM粒子の集合体として表現し、船体をFEMで扱う。preCICEやMpCCIで連成する。計算規模は数百万DEM粒子になることもある。
気候変動と氷荷重
気候変動で北極の氷が減っていますが、シミュレーション研究への影響は？
氷が薄くなる一方で、北極海航路の利用が増えて氷遭遇確率が変化する。確率論的氷荷重評価の重要性が増している。また、多年氷のリッジ（氷の畝）構造が変化しており、従来の経験式の見直しが必要になっている。
北極海航路の開通と氷-構造連成解析の「需要爆発」
地球温暖化による北極海の海氷減少は、北極海航路（NSR：北極海を経由する欧州-アジア最短ルート）の実用化を加速しています。スエズ運河経由より約40%短いこの航路が定期商用化されれば、砕氷能力を持つ商船の需要が急増します。これが氷-構造相互作用解析の研究需要を押し上げており、フィンランド、ノルウェー、カナダ、日本の造船・海洋研究機関が協力して「極海フリート設計基準の近代化」プロジェクトを進めています。特に焦点なのが「マルチイヤーアイス（多年氷）」との衝突——厚さ3m超、圧縮強度が季節氷の2〜3倍というこの氷との衝突シミュレーションは、現行の氷級規則では十分に評価できないとされています。先端技術がそのまま国際規則改定に直結する稀有な分野です。
氷-構造相互作用のトラブル対応
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の実験データベースとの比較も有用だ。接触面積-圧力関係のスケール効果を考慮することが重要だ。
「氷が壊れずに船を突き抜けた」——破壊基準の設定ミスで起きる計算のバグ
氷-構造衝突のFEM解析でよく起きるバグが「氷が破壊されずに船体を貫通してしまう」現象です。原因はほぼ100%「氷の破壊基準の閾値が高すぎる」か「要素削除の設定を忘れている」かのどちらかです。氷は脆性材料なので、引張応力が破壊強度（一般的な海氷で1〜2MPa）を超えたら要素を削除するエロージョン（element erosion）を使うのが標準的な対処法です。ただし要素削除しすぎると「質量が消える」ため、エネルギー収支がおかしくなる副作用があります。経験則では、削除される質量の合計が初期モデルの5%を超えたら結果を信頼してはいけません。また、ひずみ速度依存の破壊強度を使わないと「低速衝突で氷が割れない」という非物理的な結果になるので、速度依存材料モデルの設定は必須です。
なった
詳しく
報告