4. 疲劳损伤本构模型
复合材料层合结构在循环载荷下，损伤不是突然出现的，而是从基体微裂纹、界面脱粘、层间分层到纤维断裂逐步累积。若只用线弹性本构，仿真无法解释刚度随循环次数下降、迟滞环偏移和最终失稳。疲劳损伤本构模型的任务，就是把“损伤如何演化”与“材料刚度如何退化”写成可计算的数学关系，并嵌入有限元求解流程。
连续损伤力学的基本框架
连续损伤力学（CDM）的核心思想是：用一个标量或张量损伤变量 D 描述材料内部微缺陷引起的承载能力下降。以单轴拉伸为例，有效应力可写为：
σ̃ = σ / (1 - D)其中 σ 为名义应力，σ̃ 为有效应力，D ∈ [0,1]。D=0 表示无损伤，D=1 表示完全失效。在应变等效假设下，损伤材料的本构关系为：
σ = (1 - D) · E₀ · εE₀ 为初始弹性模量。可见损伤变量直接以乘性因子 (1-D) 作用于刚度，这就是刚度退化的数学来源。对于正交各向异性复合材料，更常用的是将 D 扩展为方向性变量，分别描述纤维方向、横向和剪切刚度的退化。
损伤变量定义与刚度退化模型
复合材料中，损伤变量通常按失效模式分别定义。以二维平面应力层合板为例，可定义：
- 纤维拉伸损伤 dft：对应 E₁ 退化，控制纤维方向承载能力；
- 纤维压缩损伤 dfc：对应纤维屈曲或扭结；
- 基体拉伸损伤 dmt：对应横向 E₂ 退化；
- 基体压缩损伤 dmc：对应横向压溃；
- 剪切损伤 ds：对应 G₁₂ 退化。
刚度退化模型一般写成：
E₁ = (1 - d_ft)(1 - d_fc) · E₁₀
E₂ = (1 - d_mt)(1 - d_mc) · E₂₀
G₁₂ = (1 - d_s) · G₁₂₀这种“按模式退化”的方式比单一 D 更贴近复合材料实际，因为横向裂纹对 E₁ 影响很小，却会显著降低 E₂ 和 G₁₂。工程中常用的退化方案有：
| 退化方式 | 刚度折减 | 特点 |
|---|---|---|
| 突然退化 | 直接降至 0.1E₀ 或 0.01E₀ | 简单，但易导致收敛困难 |
| 线性退化 | 随 D 线性下降 | 稳定，但耗散偏大 |
| 指数退化 | E = E₀·exp(-kD) | 平滑，适合疲劳累积 |
疲劳损伤演化方程
疲劳损伤本构的关键在于演化方程。常用形式为：
dD/dN = A · (Δσ / (1 - D))^m其中 Δσ 为循环应力幅，A、m 为材料参数，N 为循环次数。该式说明：损伤速率随有效应力幅增大而加快，且 D 越大，有效应力越高，损伤加速，形成正反馈。对复合材料，可对每种失效模式分别写出演化方程，例如基体拉伸损伤：
dD_mt/dN = A_mt · (Δσ₂ / (1 - D_mt))^m_mt参数 A、m 通常由疲劳试验的 S-N 曲线或刚度退化曲线拟合得到。例如某碳纤维/环氧层合板在 R=0.1、频率 5 Hz 下，横向疲劳试验测得 E₂ 随 N 下降约 15% 时出现饱和，可据此反推 A_mt 与 m_mt。这些数值属于示例性假设，实际应用需以材料试验为准。
在有限元中的实现流程
将上述模型嵌入有限元，通常按以下循环推进：
- 在当前载荷步计算各积分点应变 ε 和名义应力 σ；
- 根据失效准则判断是否进入损伤演化，计算有效应力 σ̃ = σ/(1-D)；
- 用演化方程更新 D，并更新刚度矩阵 E(D)；
- 判断 D 是否达到临界值，若达到则标记失效并进一步退化；
- 进入下一循环块，重复上述过程直至达到目标寿命或结构失稳。
为提升效率，工程中常采用“跳跃循环”策略：在损伤速率较小时，一次计算多个循环块 ΔN，而不是逐循环求解。例如当 dD/dN < 1e-6 时，取 ΔN=1000；当损伤加速后，自动减小 ΔN。这样可在保证精度的同时显著降低计算量。
应用要点与常见问题
在复合材料疲劳仿真中，损伤本构模型能否给出可信结果，取决于三件事：损伤变量是否按失效模式正确解耦、演化参数是否来自对应铺层和载荷比、刚度退化是否与单元类型和网格尺寸匹配。若把单向板的参数直接用于多向层合板，横向约束会显著改变基体裂纹扩展速率，预测寿命可能偏差数倍。