層間剥離解析 — トラブルシューティングガイド
より充実した内容を delamination.html でご覧いただけます。
層間剥離解析のトラブル
層間剥離解析でよくあるトラブルを教えてください。
CZMは非線形解析の中でも特に収束が難しい。
CZMで収束しない
剥離が進展し始めると収束しなくなります。
対策（段階的に）：
1. 粘性正則化 — *DAMAGE STABILIZATION で $\eta = 10^{-5}$ 程度
2. 増分を小さく — 最小増分 $10^{-10}$
3. 安定化法 — *STATIC, STABILIZE
4. 陽解法に切り替え — Abaqus/Explicit
陽解法に切り替えることが多いですか？
衝撃解析は元々陽解法だから問題ないが、準静的な問題（DCB試験等）でも陽解法が使われることがある。質量スケーリングで安定時間増分を大きくし、運動エネルギー/全エネルギー < 5%を確認する。
剥離面積が実験と合わない
FEMの剥離面積が実験のCスキャンと一致しません。
確認項目：
- 界面強度が正しいか — 高すぎると剥離が遅れ、低すぎると早く広がる
- $G_c$ が正しいか — DCB/ENF試験のシミュレーションで検証済みか
- CZMメッシュが十分細かいか — プロセスゾーンに3要素以上
- 層間の位置が正しいか — 実験で剥離が起きる層間にCZMを配置しているか
「どの層間に剥離が起きるか」の予測も難しそうですね。
全層間にCZMを入れれば自動的に最弱の界面が剥離するが、計算コストが膨大。実務では繊維角が急変する層間（例: 0°/90°界面）に優先的にCZMを配置する。
初期剛性の問題
CZMの初期剛性が適切でないとどうなりますか？
初期剛性の目安：$K = \alpha E_3 / t_{ply}$、$\alpha = 10 \sim 100$。AbaqusのデフォルトはPENALTY STIFFNESSの自動計算。手動設定の場合はこの範囲で試す。
まとめ
層間剥離のトラブル対処、整理します。
「DCB/ENFが合わなければ先に進めない」。層間剥離解析の鉄則ですね。
クーポンレベルの検証なしに構造レベルの予測を信じてはいけない。これは複合材解析だけでなく、全てのFEM解析に共通する原則だ。
CZMデラミネーション解析で界面剛性ペナルティが大きすぎる場合
CZMのペナルティ剛性（cohesive stiffness K=E/t₀）が大きすぎるとFEMの条件数が悪化し収束しなくなる。K=10⁶〜10⁷ N/mm³（CFRP積層板で界面厚t₀=0.01mmの場合）が実用的な範囲で、これを超えると数値的不安定が生じる。収束しない場合はKを1桁下げ、要素サイズをコヒーレントゾーン長さの1/3以下に細分化してから再試行するとよい。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——層間剥離解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告