4. 分层断裂力学基础
分层是复合材料层合板最典型的层间失效形式。要判断一条预置裂纹或界面弱区会不会扩展、沿哪个方向扩展、扩展多快，仅靠应力强度已不够，因为层间裂纹尖端处在各向异性、非均质且常伴随纤维桥联的界面环境中。断裂力学用能量释放率 G 描述“裂纹扩展单位面积可释放的弹性能”，用断裂韧性 Gc 描述“界面抵抗扩展的能力”，二者比较即可给出扩展判据。本节围绕 G、Gc、混合模式、界面本构与分层准则展开。
能量释放率：从柔度到 G
对含分层试样做双悬臂梁（DCB）试验时，可把试样视为两片悬臂梁。设加载点位移为 δ、载荷为 P，试样柔度 C=δ/P。线弹性条件下，能量释放率可由柔度法给出：
G = P² / (2b) · dC/da
其中 b 为试样宽度，a 为分层长度。若已知梁理论柔度 C=2a³/(3EI)，则 dC/da=2a²/(EI)，代入得 G= P²a²/(bEI)。这说明：同一载荷下，分层越长，G 越大；同一分层长度下，载荷越大，G 越大。工程中常用修正梁理论（MBT）处理大变形和根部转动，形式为 C=(2(a+Δ)³)/(3EI)，Δ 为根部修正量，可由柔度-裂纹长度曲线拟合得到。
需要强调参照约定：本文约定 G 为裂纹扩展单位面积释放的能量，裂纹扩展方向为分层前沿推进方向，正 G 表示驱动扩展。若采用不同符号约定，必须在报告中明确，避免与闭合积分结果混淆。
断裂韧性：Gc 的测量与含义
断裂韧性 Gc 是材料/界面属性，不是结构响应。DCB 试验中，从载荷-位移曲线偏离线性点或最大载荷点取临界载荷 Pc，代入 G 公式即得 I 型层间断裂韧性 GIc。典型碳纤维/环氧层合板 GIc 约在 0.15～0.35 kJ/m² 量级，具体值随树脂、纤维体积分数、铺层和界面处理变化。示例：某试样 b=20 mm，2h=3 mm，E11=120 GPa，a=40 mm，Pc=120 N，按 MBT 估算 GIc≈0.28 kJ/m²，此值仅作教学示例，不代表实测结果。
Gc 不是常数：它随混合模式、加载速率、温度和湿度变化。因此工程中必须区分 GIc、GIIc 以及混合模式下的 Gc(ψ)，其中 ψ 为混合模式比。
混合模式：I 型、II 型与混合比
层间裂纹可承受张开（I 型）、滑移（II 型）和撕开（III 型）三种模式。实际结构多为 I/II 混合。混合模式比常用两种定义：
- 能量比：GII/GI 或 GII/(GI+GII)；
- 相位角：ψ = arctan(√(GII/GI))，或按应力强度定义的相位角。
混合模式试验常用 MMB（混合模式弯曲）夹具，通过改变加载杠杆比获得不同 ψ。典型规律是：随着 II 型比例增加，界面断裂韧性先略升后趋稳或下降，且裂纹面摩擦和纤维桥联会显著改变 R 曲线。因此，用纯 I 型 GIc 直接判定混合模式分层往往偏保守或偏危险，必须采用与相位角对应的准则。
界面本构：从双线性到牵引-分离
界面本构描述层间界面在损伤前的弹性、损伤起始和损伤演化。最常用的是双线性牵引-分离律（CZM）：
损伤起始： max(⟨t_n⟩/t_n⁰, t_s/t_s⁰, t_t/t_t⁰) = 1
损伤演化： G_c = ∫ t dδ （曲线下面积）
其中 tn、ts、tt 为法向和两个切向牵引，t0 为界面强度，δ 为分离位移，⟨·⟩ 为 Macaulay 括号，表示仅压缩不损伤。双线性律需要输入界面强度 t0、断裂韧性 Gc 和初始刚度 K。K 通常取 105～107 MPa/mm 量级，过大会引起数值振荡，过小会引入虚假柔度。Gc 由试验获得，t0 可由界面剪切或拉伸试验估计。混合模式演化常用 B-K 或幂律准则将 GI、GII 耦合。
分层准则：从判据到工程判断
分层准则把 G 与 Gc 组合成单一判据。常用形式包括：
| 准则 | 表达式 | 适用 |
|---|---|---|
| 幂律准则 | (GI/GIc)α+(GII/GIIc)α=1 | α 常取 1～2，偏保守 |
| B-K 准则 | GIc+(GIIc−GIc)(GII/GT)η=GT | η 由混合模式试验拟合 |
| Reeder 准则 | 分段线性/二次插值 | 需完整 Gc(ψ) 曲线 |
工程判断流程可概括为：先由有限元或解析解得到分层前沿的 GI、GII；再按相位角选取对应 Gc；最后代入准则计算损伤因子。若因子≥1，则判定分层扩展。需要注意，分层前沿的 G 分布并不均匀，边缘处常出现奇异或高值区，网格收敛性和提取方法（VCCT、J 积分、CZM）会显著影响结果。
工程判断中的关键细节
第一，G 的提取方法要与准则匹配。VCCT 适合脆性界面，CZM 适合有损伤演化区的界面，J 积分适合大变形。第二，混合模式相位角定义必须与试验标定一致，否则 Gc(ψ) 曲线会错位。第三，分层前沿的网格密度直接影响 G 峰值，建议在分层前沿做局部加密并做网格收敛性检查。第四，环境因素不可忽略：湿热会降低树脂模量和界面强度，使 GIc 下降，而加载速率提高可能使 Gc 上升。第五，若出现纤维桥联，R 曲线上升，用单一 Gc 会低估韧性，应引入桥联牵引律或等效 R 曲线。
综上，分层断裂力学的核心链条是：由结构响应计算 G → 由试验标定 Gc(ψ) → 用混合模式准则比较 → 判断分层是否扩展。掌握这条链条，才能在有限元建模、试验设计和损伤容限评估中做出一致且可复现的判断。