3. 疲劳损伤仿真力学基础

复合材料疲劳损伤仿真的可靠性，首先取决于力学基础是否自洽。与各向同性金属不同，复合材料层合板的疲劳问题同时涉及单层各向异性、铺层角度、层间应力与多模式失效，因此仿真链条必须从应力应变分析出发，经层合板理论传递到单层应力状态，再由失效准则判定损伤起始与模式，最后接入疲劳寿命预测框架。本节按这一逻辑展开，并给出可操作的实现路径。

一、应力应变分析：从三维退化到平面应力

单层板在材料主方向（1为纤维方向，2为面内垂直纤维方向，3为厚度方向）下是正交各向异性材料。工程仿真中通常采用平面应力假设，即 σ₃=τ₂₃=τ₁₃=0，此时单层本构关系为：

σ₁   = Q₁₁·ε₁ + Q₁₂·ε₂
σ₂   = Q₁₂·ε₁ + Q₂₂·ε₂
τ₁₂  = Q₆₆·γ₁₂

其中：
Q₁₁ = E₁/(1-ν₁₂ν₂₁)
Q₂₂ = E₂/(1-ν₁₂ν₂₁)
Q₁₂ = ν₁₂E₂/(1-ν₁₂ν₂₁)
Q₆₆ = G₁₂
ν₂₁ = ν₁₂·E₂/E₁

以典型碳纤维/环氧单层为例（示例参数，非实测）：E₁=135 GPa，E₂=9.5 GPa，G₁₂=5.0 GPa，ν₁₂=0.30。代入可得 Q₁₁≈136.6 GPa，Q₂₂≈9.6 GPa，Q₁₂≈2.9 GPa，Q₆₆=5.0 GPa。这些参数是后续所有层合板计算的输入，必须与材料批次、纤维体积分数和固化工艺对应，不能随意套用。

当铺层角度为 θ 时，需将材料主方向刚度矩阵通过转轴公式变换到层合板整体坐标系（x-y）。变换后得到 Q̄ 矩阵，其中 Q̄₁₆、Q̄₂₆ 表征拉剪耦合。例如 θ=45° 时，Q̄₁₆ 与 Q̄₂₆ 显著非零，这正是 ±45° 铺层在剪切载荷下产生复杂应力状态的原因。仿真中若忽略转轴变换或角度符号约定，会直接导致应力分布错误。

二、层合板理论：经典层合板理论（CLT）

经典层合板理论基于 Kirchhoff 假设：直法线保持直且垂直于中面，忽略厚度方向剪切与正应力。层合板中面应变 ε⁰ 与曲率 κ 通过 ABD 矩阵与合力 N、合力矩 M 关联：

[N]   [A B] [ε⁰]
[M] = [B D] [κ ]

A_ij = Σ Q̄_ij·(z_k - z_{k-1})
B_ij = ½·Σ Q̄_ij·(z_k² - z_{k-1}²)
D_ij = ⅓·Σ Q̄_ij·(z_k³ - z_{k-1}³)

A 为拉伸刚度，D 为弯曲刚度，B 为拉弯耦合刚度。对称铺层（如 [0/90]ₛ）时 B=0，可避免固化后翘曲。以 [0/90/±45]ₛ 层合板为例，共8层，单层厚度0.125 mm，总厚1 mm。通过 ABD 矩阵求解给定载荷下的 ε⁰ 和 κ，再逐层计算应变：

ε_k = ε⁰ + z_k·κ
σ_k = Q̄_k · ε_k

这一步是疲劳仿真的核心：每个单层的应力状态不同，疲劳损伤必然从应力最危险的单层开始。仿真中必须输出每层的 σ₁、σ₂、τ₁₂ 随载荷步的变化，而不是只看层合板整体响应。

三、失效准则：从静强度到疲劳判据

复合材料失效模式包括纤维拉伸/压缩、基体拉伸/压缩、纤维-基体剪切等。常用准则有最大应力/最大应变准则、Tsai-Hill、Tsai-Wu 和 Hashin 准则。Hashin 准则因能区分失效模式，在疲劳仿真中应用最广。其二维形式为：

纤维拉伸 (σ₁≥0)：
(σ₁/X_T)² + (τ₁₂/S₁₂)² = 1

纤维压缩 (σ₁<0)：
(σ₁/X_C)² = 1

基体拉伸 (σ₂≥0)：
(σ₂/Y_T)² + (τ₁₂/S₁₂)² = 1

基体压缩 (σ₂<0)：
(σ₂/2S₂₃)² + [(Y_C/2S₂₃)²-1]·(σ₂/Y_C) + (τ₁₂/S₁₂)² = 1

示例强度参数：X_T=1500 MPa，X_C=1200 MPa，Y_T=50 MPa，Y_C=200 MPa，S₁₂=70 MPa。若某层应力为 σ₁=600 MPa，σ₂=20 MPa，τ₁₂=30 MPa，则纤维拉伸判据值=(600/1500)²+(30/70)²=0.16+0.18=0.34<1，基体拉伸判据值=(20/50)²+(30/70)²=0.16+0.18=0.34<1，均未失效。但疲劳载荷下，即使初始判据值远小于1，循环累积仍会导致损伤演化。

四、疲劳寿命预测框架

复合材料疲劳寿命预测主要有三类框架：S-N 曲线法、剩余强度/剩余刚度退化法和损伤力学法。工程仿真中常采用“疲劳损伤变量 + 材料性能退化 + 失效准则”的耦合框架。其基本流程为：

1. 通过 CLT 计算每个单层在每个循环载荷下的应力幅值 σ_max、σ_min 和应力比 R。
2. 将单层应力代入疲劳判据（如 Hashin 疲劳准则或 S-N 曲线），计算该层的疲劳损伤增量 ΔD。
3. 累积损伤 D=ΣΔD，当 D≥1 时判定该层失效。
4. 失效层的刚度按退化准则折减（如纤维失效折减 E₁ 至 10%，基体失效折减 E₂、G₁₂ 至 20%），重新计算层合板应力分布。
5. 重复上述过程，直到层合板整体失去承载能力，记录总循环数 N_f。

以 [0/90/±45]ₛ 层合板在 R=0.1、σ_max=300 MPa 下的示例计算为例：初始状态 90° 层横向应力最高，基体损伤最先累积；当 90° 层 D 达到 1 后，其 E₂ 和 G₁₂ 折减，载荷向 0° 层转移，0° 层纤维损伤开始加速。最终失效由 0° 层纤维断裂控制。这一“基体先损伤、纤维后断裂”的顺序，是复合材料疲劳的典型特征，也是仿真中必须捕捉的关键过程。

五、仿真实现中的关键约定与检查点

在编写或使用疲劳仿真脚本时，以下约定必须统一并在文档中明确：

• 坐标系与角度符号：1轴为纤维方向，2轴为面内垂直纤维方向，3轴为厚度方向；铺层角 θ 从 x 轴逆时针转到 1 轴为正。全文及所有图示保持一致。
• 应力符号：拉伸为正，压缩为负。Hashin 准则中纤维压缩和基体压缩公式不同，符号判断错误会导致失效模式误判。
• 单层厚度与总厚：示例中单层 0.125 mm，8层总厚 1 mm。若使用不同厚度，ABD 矩阵和应力计算全部需要重算。
• 刚度退化准则：纤维失效后 E₁ 折减至 10%，E₂ 和 G₁₂ 折减至 20%；基体失效后 E₂ 和 G₁₂ 折减至 20%，E₁ 不变。退化系数需与试验或文献依据对应，不能随意设定。
• 收敛判断：每次迭代后检查最大位移或载荷是否收敛，若层合板已失去承载能力则终止并输出 N_f。

综上，复合材料疲劳损伤仿真的力学基础是一条从单层本构到层合板响应、再到单层失效判定和疲劳累积的完整链条。每一步的输入参数、坐标系约定和判据选择都会影响最终寿命预测结果。实际建模时，建议先用手算或小规模脚本验证单层应力和失效判据值，再接入疲劳循环和刚度退化迭代，避免在复杂模型中直接调试导致问题定位困难。