📚 课程目录第4章损伤建模的数值框架

4. 损伤建模的数值框架

复合材料层合板的失效不是单一裂纹的瞬时扩展，而是基体开裂、界面脱粘、纤维断裂与层间分层相互耦合的渐进过程。要在有限元中复现这一过程，必须把“材料退化”和“裂纹形成”这两类物理现象分别用合适的数值框架描述。本节围绕连续损伤力学（CDM）、断裂力学（FM）、内聚力模型（CZM）三条主线，说明损伤变量如何演化、本构如何积分，以及隐式与显式求解在损伤问题中的取舍。

一、三类框架的分工与适用边界

连续损伤力学用一组损伤变量 d（0≤d<1）描述材料刚度的渐进退化，典型形式为 σ=(1−d)C:ε。它适合描述基体微裂纹弥散分布、纤维束内损伤累积等“分布型”损伤，网格尺寸对结果影响相对可控，但无法直接给出裂纹面的几何位置。断裂力学以能量释放率 G 为判据，关注裂纹尖端的应力强度因子 K 或 G，适合已存在宏观裂纹的扩展分析，需要预设初始裂纹或依赖网格重划分。内聚力模型则在潜在开裂面上插入零厚度界面单元，用牵引—分离（traction-separation）本构直接刻画界面从承载到完全失效的过程，是层间分层和胶接界面建模的首选。

工程实践中三者常组合使用：层内纤维/基体损伤用 CDM，层间分层用 CZM，宏观裂纹扩展用 VCCT 或 XFEM 辅助判断。选择框架时先问两个问题：损伤是弥散的还是局部的？开裂面是否已知？前者选 CDM，后者选 CZM；若裂纹路径完全未知且需追踪，则考虑断裂力学类方法。

二、损伤变量演化：从判据到软化

CDM 的核心是损伤变量的演化律。以 Hashin 判据为例，纤维拉伸失效的判据为 Ff^t=(σ11/X_T)²+α(τ12/S12)²，当 Ff^t≥1 时损伤启动。损伤变量 d 的演化通常写成等效位移 δeq 的函数：

d = (δeq − δ0) / (δf − δ0)   （线性软化）
δ0 = ε0 · Lc,  δf = 2Gc / (X · Lc)

其中 Lc 是单元特征长度，Gc 是断裂能。这个公式的关键作用是让软化段消耗的能量与网格尺寸无关——即“能量正则化”。若不做这一步，网格越细，损伤局部化越严重，总耗散能随网格加密而下降，结果不收敛于真实值。工程上常取 Lc 为单元体积的立方根或面积的开方，并在 δf 计算中显式引入。

损伤变量的更新必须满足不可逆性：d_n+1=max(d_n, d_new)。若某积分点卸载后重新加载，损伤不能恢复。这一约束在隐式算法中通过历史变量存储实现，在显式算法中则直接比较当前步与上一步的 d 值。

三、本构积分算法：隐式与显式的实现差异

损伤本构的积分本质上是求解 σ_n+1=(1−d_n+1)C:(ε_n+1−ε_n+1^p) 与损伤演化方程的耦合系统。隐式算法在每个增量步内用 Newton 迭代同时求解位移和损伤变量，需要推导一致切线刚度 C^ep=∂σ/∂ε。对于线性软化 CDM，C^ep 可解析给出；对于复杂的 Puck 或 LaRC 判据，则常用数值微分近似。

显式算法（如 Abaqus/Explicit）不组装全局刚度矩阵，而是用中心差分推进：M·a_n=F_ext−F_int，然后更新速度和位移。损伤变量在单元层面直接根据当前应变更新，无需迭代。其代价是稳定时间步受最小单元尺寸和材料波速限制：Δt≤Lmin/cd，cd=√(E/ρ)。对于碳纤维复合材料，cd 约 8000–10000 m/s，若最小单元为 0.5 mm，则 Δt≈5×10⁻⁸ s。模拟 10 ms 的冲击过程需要约 2×10⁵ 步，计算量可观，但显式天然适合处理接触、大变形和单元删除。

选择建议：准静态层合板压缩/拉伸失效优先用隐式（Abaqus/Standard + CDM），配合粘性正则化（viscosity）改善软化段收敛；低速冲击、鸟撞、爆炸等动态问题用显式（Abaqus/Explicit + CDM/CZM），并注意质量缩放不能改变惯性效应主导的物理过程。

四、内聚力模型：界面分层的专用工具

CZM 的核心是牵引—分离本构。以双线性模型为例，法向牵引 tn 随分离位移 δn 先线性上升至界面强度 tn⁰，达到 δn⁰ 后线性软化至零，对应的断裂能 GIc=½tn⁰·δn^f。混合模式（I/II）下常用 BK 准则：

(GI/GIc)^η + (GII/GIIc)^η = 1

其中 η 通常取 1.5–2.5（环氧树脂基复合材料常用 2.0）。在有限元中，CZM 通过零厚度界面单元或接触面粘接行为实现。界面单元的本构积分同样需要处理加载—卸载—再加载路径：卸载时按初始刚度返回原点，再加载时沿软化段继续，损伤变量 dczm=1−t/t⁰ 单调递增。

实际建模时，CZM 参数（tn⁰、GIc、GIIc、η）通常来自 DCB（双悬臂梁）和 ENF（端缺口弯曲）试验。若缺乏试验数据，可参考典型值：碳纤维/环氧层间 tn⁰≈50–60 MPa，GIc≈0.2–0.3 N/mm，GIIc≈0.6–1.0 N/mm。界面单元刚度 Kn 一般取 10⁵–10⁶ MPa/mm，过大导致收敛困难，过小则引入虚假柔度。

五、框架组合与求解策略的工程判断

一个完整的层合板冲击损伤模型通常包含：层内 CDM（纤维/基体损伤）、层间 CZM（分层）、以及可能的断裂力学判据（宏观裂纹）。隐式求解时，CZM 的软化段和 CDM 的损伤局部化都会导致收敛困难，常用粘性正则化（viscosity coefficient 取 10⁻⁵–10⁻³）和线搜索辅助。显式求解则需控制沙漏模式和单元删除后的能量平衡，检查 ALLKE/ALLIE 比值是否合理。

判断结果是否可信，不能只看损伤云图。应检查：损伤耗散能（ALLDMD）是否随加载单调增加且量级合理；分层面积是否与超声 C 扫描或显微观察一致；网格加密后载荷—位移曲线的峰值和软化段是否趋于稳定。若峰值载荷随网格加密持续下降，说明能量正则化未生效或 Lc 取值不当。

常见陷阱：在显式分析中直接使用隐式 CDM 参数而不检查特征长度，会导致损伤带宽度等于单元尺寸，结果严重依赖网格。务必在材料子程序中显式使用 Lc 计算 δf，或在 Abaqus 中启用“Damage Evolution — Energy”类型并指定 Gc。

综上，CDM、FM、CZM 不是互斥选项，而是针对不同损伤尺度和开裂模式的工具组合。理解损伤变量的演化律、本构积分的一致切线、以及隐式/显式求解的稳定性条件，是搭建可靠复合材料损伤模型的基础。