2. 复合材料疲劳损伤基础

复合材料疲劳与金属疲劳最大的区别在于：金属通常可视为各向同性材料，疲劳裂纹萌生与扩展主要由单一应力幅控制；而纤维增强复合材料是各向异性的，其疲劳行为强烈依赖铺层方向、载荷方向、应力比以及纤维/基体/界面的协同响应。理解这些差异，是后续建立疲劳损伤仿真模型的前提。

一、各向异性疲劳特征

以典型的碳纤维/环氧层合板为例，沿纤维方向（0°）的疲劳强度远高于垂直纤维方向（90°）。在单向层合板中，0°方向拉伸疲劳往往由纤维断裂主导，寿命较长；90°方向拉伸疲劳则由基体开裂主导，寿命显著缩短。对于多向层合板，如[0/±45/90]s，不同铺层角承担不同应力分量，损伤往往先在偏轴层或横向层中萌生，再逐步向其他层扩展。

各向异性还带来两个重要现象：一是拉-拉、拉-压、压-压疲劳的损伤机制不同，压-压载荷下可能出现纤维微屈曲或分层；二是剪切耦合效应，面内剪切与正应力之间存在耦合，使疲劳损伤路径不再是单一裂纹面，而是多模式损伤的组合。因此，在仿真中不能简单套用金属的等效应力准则，而应分别跟踪纤维拉伸、纤维压缩、基体拉伸、基体压缩和面内剪切等损伤变量。

二、主要损伤模式

复合材料疲劳损伤通常包含四种基本模式：基体开裂、界面脱粘、分层和纤维断裂。基体开裂是最早出现的损伤，表现为横向层中沿纤维方向的微裂纹；界面脱粘发生在纤维与基体结合处，削弱载荷传递能力；分层是层间裂纹，常由自由边效应或冲击损伤引发；纤维断裂是最终失效的主导模式，一旦大量纤维断裂，层合板承载能力迅速下降。

这些损伤模式并非独立发生，而是相互耦合。例如，基体开裂会诱发层间应力集中，促进分层扩展；分层又会使相邻层失去支撑，加速纤维屈曲和断裂。在疲劳仿真中，常用内变量或损伤变量描述各模式的演化，并通过刚度退化模型将损伤与材料性能下降联系起来。

三、S-N曲线及其各向异性表现

S-N曲线描述应力幅与疲劳寿命之间的关系，通常以应力幅为纵轴、失效循环数为横轴（对数坐标）。对于复合材料，S-N曲线不是一条单一曲线，而是随铺层方向、应力比和载荷类型变化的一组曲线。以T300/环氧单向板为例，在应力比R=0.1的拉-拉疲劳下，0°方向在约10⁶次循环时仍能保持约60%的静强度；而90°方向在相同寿命下仅能保持约25%～30%的静强度。±45°层合板的疲劳寿命通常介于两者之间，且对剪切应力更为敏感。

在仿真中，S-N曲线常用于两种方式：一是作为寿命预测的经验输入，通过插值或拟合得到不同应力水平下的失效循环数；二是作为损伤演化模型的标定依据，将S-N曲线转化为损伤变量随循环数的演化规律。需要注意的是，复合材料的S-N曲线在双对数坐标下往往呈现分段线性特征，低周与高周区域的斜率不同，拟合时应分段处理。

四、剩余刚度与剩余强度概念

剩余刚度E(n)定义为材料在经历n次循环后仍具有的刚度，通常用初始刚度E₀归一化表示。在疲劳过程中，基体开裂和分层会使刚度逐渐下降，但下降速率通常较慢；当损伤进入纤维断裂阶段时，刚度下降加速。剩余刚度模型常写为E(n)/E₀ = 1 − f(D)，其中D为损伤变量，f为退化函数。常见的退化函数包括线性、幂律和指数形式，具体形式需根据材料体系标定。

剩余强度σ(n)定义为材料在经历n次循环后仍能承受的极限应力。与刚度不同，剩余强度在疲劳前期下降较慢，接近寿命终点时急剧下降。剩余强度模型常写为σ(n)/σ₀ = 1 − (n/N)^α，其中α为材料常数，N为失效循环数。该模型可用于判断当前损伤状态是否接近最终失效，也可作为仿真中失效判据的输入。

五、仿真中的关键处理

在疲劳损伤仿真中，通常需要将上述概念转化为可计算的模型。第一步是定义损伤变量，例如用d_f、d_m、d_s分别表示纤维、基体和剪切损伤；第二步是建立损伤演化方程，将循环数、应力状态与损伤变量联系起来；第三步是耦合刚度退化，将损伤变量映射到材料刚度矩阵；第四步是设置失效判据，当剩余强度低于当前应力或损伤变量达到临界值时判定失效。

需要特别注意的是，各向异性使损伤演化具有方向性，不同铺层的损伤速率不同，仿真中应逐层计算，而不是对整个层合板取平均。此外，分层损伤通常需要单独处理，可采用界面单元或内聚力模型描述层间裂纹的萌生与扩展。只有将各向异性特征、损伤模式、S-N曲线和剩余刚度/强度模型统一在一个框架下，才能较准确地模拟复合材料疲劳损伤的演化过程。