層間剥離解析 — 理論と支配方程式
層間剥離とは
先生、「層間剥離(デラミネーション)」は複合材の最も危険な破壊モードですか?
その通り。層間剥離は積層板の層と層の間が剥がれる破壊で、複合材の最も一般的かつ危険な損傷モードだ。表面からは見えない「内部損傷」であり、検出が困難なことが危険性を高める。
表面から見えない! それは怖いですね。
低速衝撃(落下物、工具の落下等)で表面にほとんど痕跡がないのに、内部で広範囲に層間剥離が広がっている場合がある。これがBVID(Barely Visible Impact Damage)と呼ばれ、航空機の複合材設計で最も重要な設計条件だ。
層間剥離のメカニズム
層間剥離は層間のせん断応力($\tau_{xz}, \tau_{yz}$)と剥離応力($\sigma_z$)が界面の強度を超えたときに発生する。
破壊力学的には3つのモードがある:
| モード | 応力 | 変形 | 代表試験 |
|---|---|---|---|
| Mode I(開口) | $\sigma_z$(引張) | 層間が開く | DCB |
| Mode II(面内せん断) | $\tau_{xz}$ | 層がずれる | ENF / 4ENF |
| Mode III(面外せん断) | $\tau_{yz}$ | 層がねじれる | ECT |
Mode Iが最も危険ですか?
Mode Iの臨界エネルギー解放率 $G_{Ic}$ が最も低い($\approx 0.1 \sim 0.3$ kJ/m²)。Mode IIの $G_{IIc}$ はその2〜4倍。だからMode Iの開口が先行することが多い。
エネルギー解放率
層間剥離の進展条件はエネルギー解放率(Energy Release Rate, ERR)で記述される:
$G$ は現在のエネルギー解放率、$G_c$ は臨界値(材料特性)。
混合モード(Mode I + Mode II が同時に作用)の場合は混合モード基準を使う:
$\alpha = \beta = 1$ がBenzeggagh-Kenane(BK)基準の特殊ケース。
FEMでの層間剥離モデル化
FEMで層間剥離をどうモデル化しますか?
2つの主要手法:
1. VCCT(Virtual Crack Closure Technique)
亀裂先端の節点力と開口変位からエネルギー解放率を計算する。既存の亀裂の進展を追跡する手法。
2. CZM(Cohesive Zone Model)
界面にコヒーシブ要素を配置し、応力-開口変位の構成則(トラクション-セパレーション則)で剥離を表現。亀裂の核生成と進展の両方を扱える。
CZMのほうが汎用的ですか?
CZMは亀裂の位置を事前に仮定する必要がない(コヒーシブ要素を全界面に配置すればよい)。VCCTは既存亀裂の進展のみ。CZMが現在の主流だ。
まとめ
層間剥離の理論を整理します。
要点:
- 層間剥離は複合材の最も危険な損傷 — 表面から見えない(BVID)
- 3つの破壊モード — Mode I(開口)、II(せん断)、III(ねじり)
- エネルギー解放率 $G \geq G_c$ で進展 — 混合モード基準
- CZM(コヒーシブゾーンモデル)が主流 — 核生成+進展を両方扱える
- VCCTは既存亀裂の進展のみ — CZMより制限的だが計算は軽い
BVID(目に見えない衝撃損傷)が航空機の設計を支配するとは…。「見えない敵」と戦う設計ですね。
まさにそう。航空機の複合材設計は「最悪のBVIDが存在する前提」で行われる。だからCAI強度が設計許容値になる。
タイタニック号と安全率の教訓
「不沈」と謳われたタイタニック号は、低温でのリベット材の脆性破壊が沈没の一因とされています。現代の破壊力学CAEでは、温度依存の材料特性と応力拡大係数を計算して「その温度で本当に大丈夫か?」を事前に検証できます。技術の進歩は、過去の悲劇から学んだ結果です。
各項の物理的意味
- 慣性項(質量項):$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか? あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項(弾性復元力):$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね? あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか? 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い:「剛性が高い=強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項(荷重項):体積力 $f_b$(重力など)と表面力 $f_s$(圧力、接触力など)。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」(体積力)、タイヤが路面を押す力は「表面だけにかかる力」(表面力)。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗:荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項:レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか? いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら? 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m(メートル) | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa(パスカル)= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元(m/m) | 工学歪みと対数歪みの区別に注意(大変形時) |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³(= 10⁻⁹ tonne/mm³ for 鋼) |
| 力 $F$ | N(ニュートン) | mm系ではN、m系ではNで統一 |
数値例:片持ち梁の先端荷重(L=1m, 断面50×100mm, 鋼材E=210GPa, P=1kN)
最大たわみ δ = PL³/(3EI) = 1000×1000³/(3×210000×4,166,667) ≈ 0.381 mm 最大応力 σ = PL×(h/2)/I ≈ 12.0 MPa(降伏応力235MPaに対して安全率19.6)
メッシュ密度を変えた収束性の確認:
ポイント:要素数を4倍にしても結果は0.5%しか変わらない→8,000要素で十分収束。これが「メッシュ収束性」の確認です。
構造解析の収束問題や計算コストに課題を感じていませんか? — Project NovaSolverは、実務者が日々直面するこうした課題の解決を目指す研究開発プロジェクトです。
層間剥離解析の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →