船舶スラミング解析
理論と物理
スラミングの物理
船舶のスラミングってどういう現象ですか？
波浪中を航行する船舶の船底が水面に衝突する際に、極めて短時間に大きな衝撃圧力が発生する現象だ。船首部のbow flare slamming、船底のbottom slamming、船尾のstern slammingに分類される。衝撃圧力はMPaオーダーに達し、局所的な構造損傷やホイッピング振動を引き起こす。
支配方程式
衝撃圧力の理論解はあるんですか？
Wagner（1932）の理論が古典的だ。2Dくさび形状の水面衝突で、
ここで $V$ は衝突速度、$\beta$ はデッドライズ角、$c(t)$ は濡れ幅だ。$\beta \to 0$ で圧力は無限大に発散するため、air cushion効果を考慮した修正モデルが必要になる。
流体側は非圧縮性（衝突速度 $V \ll c_{water}$ のとき）または圧縮性のNavier-Stokes方程式を解く。自由表面はVOF法で追跡する。
構造側は弾塑性FEMだ。衝撃荷重の時間スケール（ミリ秒オーダー）と構造の固有周期の関係で、動的応答か準静的応答かが決まる。
ワグナーの楔入水理論——「水の壁」はなぜ硬いのか
船底が海面を叩くスラミングの瞬間、局所圧力は静水圧の数十〜数百倍にもなる。これを最初に理論化したのが1932年のH.ワグナーで、楔形断面が速度Vで入水するときの最大圧力はp_max ≈ ½ρ(πV/2tanβ)²という関係を導いた（βは半頂角）。重要なのは「入水角が小さいほど圧力が爆発的に増える」という点だ。たとえばβ=5°で入水速度5m/sなら、ピーク圧力は数MPaに達する。FRP製の小型高速艇でスラミングが続くと、船底パネルが「バン！」と音を立てて凹む理由がこれだ。この理論は今でもスラミング設計の出発点として使われ続けている。
各項の物理的意味
- 構造-熱連成項：温度変化による熱膨張が構造変形を誘発し、変形が温度場に影響する。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常の例】夏に線路のレールが伸びて隙間が狭くなる——温度上昇→熱膨張→応力発生の典型例。電子基板がはんだ付け後に反るのも、異なる材料の熱膨張率差による。エンジンのシリンダーブロックは高温部と低温部の温度差で熱応力が発生し、最悪の場合亀裂に至る。
- 流体-構造連成（FSI）項：流体圧力・せん断力が構造を変形させ、構造変形が流体領域を変化させる双方向の相互作用。【日常の例】強風で吊り橋のケーブルが振動する（渦励振）——風の力が構造を揺らし、揺れた構造が風の流れを変え、さらに振動が増幅する。心臓の血流と血管壁の弾性変形、航空機の翼のフラッタ（空力弾性不安定性）も典型的なFSI問題。片方向のみの連成で済む場合もあるが、変形が大きい場合は双方向連成が必須。
- 電磁-熱連成項：ジュール発熱 $Q = J^2/\sigma$ が温度上昇を引き起こし、温度変化が電気抵抗を変化させるフィードバックループ。【日常の例】電気ストーブのニクロム線は電流が流れると発熱（ジュール熱）して赤くなる——温度が上がると抵抗が変わり、電流分布も変化する。IHクッキングヒーターの渦電流発熱、送電線の温度上昇による弛み増加もこの連成の例。
- データ転写項：異なる物理場間のメッシュ不一致を補間で解決。【日常の例】天気予報で「気温のデータ」と「風のデータ」を合わせて体感温度を計算するとき、それぞれの観測地点が異なれば補間が必要——CAEの連成解析でも、構造メッシュとCFDメッシュは一般に一致しないため、界面でのデータ転写（補間）精度が結果の信頼性に直結する。
仮定条件と適用限界
- 弱連成仮定（片方向連成）：一方の物理場が他方に影響するが逆は無視可能な場合に有効
- 強連成が必要なケース：FSIでの大変形、電磁-熱連成での温度依存性が強い場合
- 時間スケールの分離：各物理場の特性時間が大きく異なる場合、サブサイクリングで効率化可能
- 界面条件の整合性：連成界面でのエネルギー・運動量保存が数値的に満たされることを確認
- 適用外ケース：3つ以上の物理場が同時に強く連成する場合、モノリシック手法が必要になることがある
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 熱膨張係数 $\alpha$ | 1/K | 鋼: 約12×10⁻⁶、アルミ: 約23×10⁻⁶ |
| 連成界面力 | N/m²（圧力）またはN（集中力） | 流体側と構造側で力の釣り合いを確認 |
| データ転写誤差 | 無次元（%） | 補間精度はメッシュ密度比に依存。5%以下が目安 |
数値解法と実装
数値手法
スラミングのCFD-FSIではどんな手法が使われますか？
衝撃を伴うため陽解法が基本だ。
| 手法 | 流体 | 構造 | 特徴 |
|---|---|---|---|
| CFD-VOF + FEM陽解法 | OpenFOAM/Fluent | LS-DYNA/Abaqus Explicit | 汎用的。FSIは弱連成が多い |
| SPH + FEM | SPH | FEM | メッシュフリー。飛沫・砕波に強い |
| ALE (LS-DYNA) | ALE fluid | FEM | *CONSTRAINED_LAGRANGE_IN_SOLID |
| BEM + FEM | パネル法 | FEM | 効率的だがスプレー表現不可 |
SPH法はスラミングに向いているんですか？
SPHはメッシュ不要で自由表面の大変形・飛沫を自然に追跡できる。ただし、圧力場に数値振動（ノイズ）が出やすいため、Riemann SPHやδ-SPH等の改良版が使われる。LS-DYNAにSPHソルバーが内蔵されており、FEM構造と直接連成できる。
時空間分解能
衝撃圧力を正確に捕えるにはどのくらいの分解能が必要ですか？
スラミング圧力のピークは0.1〜1 ms程度で消失するから、時間刻みは0.01 ms以下、空間メッシュは衝突面近傍で1〜5 mm以下が必要だ。
| パラメータ | 推奨値 |
|---|---|
| 衝突面メッシュサイズ | 1〜5 mm |
| 時間刻み | 0.01 ms以下 |
| VOF界面解像 | 最低5セル/水膜厚 |
| 圧力サンプリング | 0.001 ms以下 |
水槽試験とCFDの「1000倍の壁」——スラミング解析の現実
スラミングの数値解析は時間刻みの問題と格闘の連続だ。スラミングのピーク圧力は数ミリ秒以内に生じるため、時間刻みΔtを0.01ms以下にしないと解像できない。一方、実際の船が1波を超えるのに数秒かかる。この時間スケールの差が「1000倍の壁」だ。実務では「スラミング単体のCFD」と「船体全体の長時間FEA」を分離し、スラミング圧力をマッピングする片方向連成が主流になっている。水槽試験との比較では、ピーク圧力の±20%以内に収まればまずまず、という実務感覚がある。CFDと実験の差の主因はしばしば「気泡の混入効果」で、これが圧力を10〜30%も下げることがある。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
実践ガイド
解析手順
スラミング解析を実務で進める手順は？
2段階アプローチが一般的だ。
Phase 1: 全船耐航性解析
- ポテンシャル流BEMまたはCFDで船体運動を計算
- スラミング発生条件（相対速度、相対変位）を特定
Phase 2: 局所スラミング解析
- Phase 1の結果から衝突速度・角度を抽出
- 局所的なCFD-FSIまたはSPH-FEMで衝撃圧力・構造応答を計算
なぜ2段階に分けるんですか？
全船スケールでスラミングの時空間分解能を確保すると計算コストが天文学的になるからだ。Phase 1で「いつ、どこで、どのくらいの速度で」衝突するかを特定し、Phase 2でその局所イベントだけを高精度で解析する。
検証データ
検証に使えるベンチマーク問題はありますか？
くさび形状の水面衝突問題が標準的だ。
| ベンチマーク | 形状 | 検証対象 |
|---|---|---|
| Wagner理論 | 2Dくさび | 衝撃圧力分布 |
| Zhao & Faltinsen実験 | 2Dくさび(β=30°) | 圧力時刻歴 |
| Aarsnes弾性くさび | 2D弾性パネル | FSI応答 |
| Luo実験 | 3D船首セクション | 3D圧力分布 |
コンテナ船の「ウィッピング」——スラミング後の船全体が揺れる
スラミングが起きた直後、船体全体がまるでギターの弦を弾いたように縦振動する「ウィッピング」という現象がある。大型コンテナ船では固有周期が2〜4秒程度で、スラミング後に船体中央で数百MPaの曲げ応力スパイクが5〜10回繰り返す。2013年のMOL Comfort沈没事故（コンテナ船が真っ二つに折れた）では、このウィッピング応力が設計想定を上回った可能性が指摘されている。この事故を受けて、国際船級協会（IACS）はコンテナ船の設計基準にウィッピング応力を明示的に組み込むよう規定を改訂した（UR S11A）。スラミング-ウィッピング連成解析が設計に不可欠になった転換点だ。
解析フローのたとえ
風船を膨らませたことがありますか？ あの瞬間、実は高度な流体-構造連成が起きています。内部の空気圧（流体）がゴム壁（構造）を押し広げ→広がった壁が内部の圧力分布を変え→変わった圧力がさらに壁を変形させる…このキャッチボールを計算ステップごとに繰り返すのがFSI解析です。
初心者が陥りやすい落とし穴
「片方向連成で十分でしょ？」——この判断ミスが連成解析で最も危険です。構造の変形が微小なら確かに片方向で足りますが、心臓弁の開閉のように変形が流路を大きく変える場合、片方向では全く話になりません。目安は「変形量が代表長さの1%を超えるか」。超えるなら双方向連成は必須です。片方向で済ませてしまった場合、結果が「もっともらしいけど実は大間違い」になる——これが最も怖いパターンです。
境界条件の考え方
連成界面のデータ交換は「国境の出入国管理」と同じです。各国（物理場）には独自の法律（支配方程式）がありますが、国境（界面）で人や物（力・温度・変位）のやり取りを正確に管理しないと、両国の経済（エネルギーバランス）が崩壊します。メッシュが一致していない場合の補間は「通訳」のようなもの——誤訳（補間誤差）が小さいほど良い結果が得られます。
ソフトウェア比較
ツール比較
スラミング解析に適したソフトウェアは？
LS-DYNAのALEとSPHはどう使い分けるんですか？
ALEはメッシュベースで精度が高いが、飛沫の追跡が困難。SPHは飛沫を自然に捕捉できるが、圧力精度がやや低い。実務ではALEでスラミング圧力を評価し、SPHでグリーンウォーターや飛沫の可視化に使う使い分けが多い。
AbaqusのCEL法は使えますか？
Coupled Euler-Lagrange（CEL）法はEuler領域の流体とLagrange構造を単一ソルバーで連成する。インターフェースの設定が簡単だが、自由表面の捕捉精度はVOF法に劣る場合がある。小規模問題や概念設計段階では有用だ。
船級協会がスラミング解析ツールを「認証」する時代
DNV、Lloyd's Register、NK（日本海事協会）などの船級協会は、近年スラミング解析ツールそのものに対して「検証済みツールリスト」を整備し始めている。LR（ロイズ船級）のCompassシステムでは、スラミング解析の手順書と推奨ツール（LS-DYNA、STAR-CCM+、Abaqusなど）を指定し、それ以外のツールを使う場合は検証報告書の提出を求める。日本では国土交通省・ClassNKが2022年に「高度な数値解析ガイドライン」を改訂し、スラミング-ウィッピング連成解析の受け入れ基準を明文化した。「好きなツールで計算すればよい」時代から「船級が認めたワークフローで解析する」時代へ、業界標準が急速に整備されている。
選定で最も重要な3つの問い
- 「何を解くか」：船舶スラミング解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
エアクッション効果
船底がほぼ水平に衝突する場合、空気の効果が重要だと聞きましたが。
デッドライズ角 $\beta < 3°$ の平底スラミングでは、船底と水面の間に閉じ込められた空気がクッションとして働き、ピーク圧力を大幅に低減する。ただし、空気の圧縮性と逃げ（air escape）の競合で圧力分布が複雑になる。
この場合、二相流（水+空気）の圧縮性を考慮したCFDが必要だ。空気の圧縮によるクッション圧力は、
$\gamma = 1.4$（空気の比熱比）で断熱圧縮を仮定する。
ホイッピング応答への展開
スラミングで船全体が振動するホイッピングはどう解析するんですか？
スラミング衝撃力を全船FEモデルの過渡荷重として入力する。水弾性解析の枠組みで時間領域応答を計算し、VBM（Vertical Bending Moment）のhogおよびsag成分を評価する。ホイッピングによるVBM増分は波浪VBMの30〜50%に達することがある。
確率論的スラミング評価
不規則波中でのスラミング頻度や荷重の確率分布はどう求めるんですか？
耐航性解析で相対運動の統計量を求め、スラミング発生閾値を超える確率を計算する。超過確率は、
$v_c$ は臨界相対速度、$\sigma_v$ は相対速度の標準偏差だ。長期分布にはIACSのルールウェーブ（equivalent design wave）法が使われる。
弾塑性スラミング解析——衝撃で船底は「降伏」する
スラミング荷重が極端に大きいとき、船体構造は弾性域を超えて塑性変形することがある。通常の疲労解析では線形弾性を前提にするが、嵐での一発スラミングは降伏応力を超えることがある。最先端の研究では、CEL（Coupled Eulerian-Lagrangian）法でスラミング流体と弾塑性構造を同時に解き、塑性変形量を直接予測する試みが行われている。DNVとノルウェー科学技術大学（NTNU）の共同プロジェクトでは、波高12mの設計波に対するバルクキャリアの船首スラミングをCEL法で解析し、従来の弾性＋安全率設計より船首外板の最適板厚を約8%削減できることを示した。
トラブルシューティング
圧力ピークのメッシュ依存性
メッシュを細かくするたびにピーク圧力が変わるんですが。
スラミング圧力には本質的な特異性がある（Wagner理論でβ→0でp→∞）。メッシュ収束しない場合はピーク圧力ではなく、力の積分値（衝撃力）や構造応答（変形・応力）で評価する方が安定だ。
実務では圧力×面積の積分値（全力）がメッシュに対してより安定している。また、構造応答を評価対象にすれば、構造の慣性効果で高周波成分がフィルタされる。
SPHの圧力振動
SPHで計算した圧力にノイズが乗ります。
対策は以下の通り。
| 対策 | 効果 |
|---|---|
| δ-SPH法 | 拡散項追加で圧力を平滑化 |
| Riemann SPH | Riemann問題の解で粒子間相互作用を改善 |
| 移動最小二乗法（MLS）による後処理 | 圧力フィールドの平滑化 |
| 粒子数増加 | 分解能向上（ただし計算コスト増大） |
構造応答の過大評価
構造の応力が非現実的に大きい場合はどうしますか？
構造応答がWagner理論の静的荷重換算値の2倍以上になる場合は、動的増幅係数（DAF）の妥当性を確認すべきだ。
スラミング解析の「発散地獄」——圧力スパイクとの闘い
スラミングのCFD解析で最も頭を悩ませるのが「圧力スパイク」だ。気液界面が固体壁に激突する瞬間、計算格子のスケールより小さな現象が生じ、圧力が数値的に発散することがある。よくある症状は「時刻ゼロ付近で圧力が1000倍に跳ね上がって計算が壊れる」というもの。対策として有効なのが「圧縮性流体モデルの導入」だ。液体は非圧縮とみなすことが多いが、衝撃波的な現象では液体のわずかな圧縮性（音速約1500m/s）が圧力を有限に保つ効果がある。また、VOF法での界面の「数値厚み」を1セル以上確保することも収束安定化の基本テクニックだ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——船舶スラミング解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告