2. 复合材料力学基础与损伤表征

复合材料损伤有限元建模的可靠性，首先取决于对材料本构、层合板响应和损伤判据的准确描述。与各向同性金属不同，单层复合材料在纤维方向与垂直纤维方向上的刚度差异可达一个数量级，且面内剪切与横向拉压行为彼此耦合。若在有限元中直接套用金属的弹性模型，即使网格与边界条件完全正确，位移和应力结果也会出现系统性偏差。因此，本节从单层正交各向异性本构出发，逐步建立工程常数、坐标变换、层合板理论、刚度退化与失效判据之间的完整链条。

正交各向异性本构与工程常数

连续纤维增强单层板在宏观上可视为正交各向异性材料，具有三个相互垂直的弹性对称面。在材料主坐标系（1为纤维方向，2为面内垂直纤维方向，3为厚度方向）下，三维应力—应变关系由广义胡克定律给出：

σ = C · ε
其中 C 为 9 个独立常数构成的刚度矩阵，
工程常数形式常用 5 个面内参数描述：
E1  纵向弹性模量
E2  横向弹性模量
ν12 主泊松比
G12 面内剪切模量
G23 面外剪切模量（三维分析时需补充）

以典型碳纤维/环氧单层板为例，可取 E1=135 GPa、E2=9.5 GPa、ν12=0.30、G12=5.0 GPa、G23=3.4 GPa。这些数值仅作为教学示例，实际建模应以材料数据表或试验结果为准。由工程常数可组装柔度矩阵 S，再求逆得到刚度矩阵 C。需要注意，正交各向异性材料的泊松效应满足 Maxwell 互等关系 ν12/E1 = ν21/E2，因此 ν21 并非独立参数，不能随意赋值。

坐标变换与铺层方向

层合板中每一层的纤维方向与整体坐标系（x、y、z）通常不一致，需要将材料主坐标系下的刚度矩阵通过转轴公式变换到整体坐标系。设铺层角 θ 为纤维方向与整体 x 轴的夹角，则变换后的面内刚度矩阵 Q̄ 由 Q 矩阵经方向余弦矩阵 T 变换得到：

Q̄ = T⁻¹ · Q · T⁻ᵀ
θ = 0°   纤维沿 x 轴
θ = 90°  纤维沿 y 轴
θ = ±45° 纤维斜向，产生拉剪耦合

变换后会出现 Q̄16、Q̄26 等耦合项，它们使正应力与剪应变相互耦合。这正是 ±45° 铺层在拉伸时产生剪切变形、在剪切时产生正应力的根源。建模时若忽略铺层角或把 θ 与 −θ 混用，层合板的耦合刚度会完全错误，因此必须在有限元铺层定义中逐层核对角度与材料方向。

层合板理论与刚度矩阵

层合板由多层单层板按不同角度叠合而成，经典层合板理论（CLT）假设层间完全粘结、变形沿厚度线性分布。其核心是三个刚度矩阵：

• A 矩阵：面内刚度，由各层 Q̄ 沿厚度积分得到，决定面内载荷下的应变。
• B 矩阵：拉弯耦合刚度，当铺层关于中面不对称时不为零，拉伸会引起弯曲。
• D 矩阵：弯曲刚度，决定弯矩与曲率的关系。

以 [0/90]₂ₛ 对称层合板为例，由于铺层关于中面对称，B 矩阵为零，面内拉伸不会产生弯曲耦合，这是工程中优先采用对称铺层的原因。在有限元中，壳单元或连续壳单元通常直接输入铺层序列、每层厚度和材料方向，由求解器自动计算 A、B、D 矩阵。建模时应检查铺层顺序是否与设计一致，尤其是 ±45° 层的正负号顺序，它直接影响耦合项。

损伤变量与刚度退化

复合材料损伤不是单一裂纹，而是纤维断裂、基体开裂、界面脱粘等多种模式的累积。连续损伤力学用损伤变量 d 描述材料刚度的退化程度，d=0 表示完好，d=1 表示完全失效。常见的退化方式包括：

• 纤维损伤：d_f 控制 E1 的退化，纤维断裂后纵向承载能力急剧下降。
• 基体损伤：d_m 控制 E2 和 G12 的退化，基体开裂后横向和剪切刚度降低。

在有限元实现中，损伤变量通常由等效应变或断裂能演化方程驱动。例如，基体拉伸损伤可写为 d_m = (ε − ε₀)/(ε_f − ε₀)，其中 ε₀ 为损伤起始应变，ε_f 为完全失效应变。退化后的刚度回代到层合板刚度矩阵中，形成非线性迭代。需要特别注意的是，刚度退化过猛会导致收敛困难，工程中常采用粘性正则化或延迟损伤起始来改善数值稳定性。

失效准则与损伤起始判据

判断损伤是否起始，需要选择合适的失效准则。常用准则包括：

以 Hashin 准则为例，纤维拉伸失效判据为 (σ11/X_T)² + (τ12/S12)² ≥ 1，基体拉伸失效判据为 (σ22/Y_T)² + (τ12/S12)² ≥ 1。其中 X_T、Y_T 分别为纵向和横向拉伸强度，S12 为面内剪切强度。当某一层满足判据时，该层对应损伤变量开始演化，刚度按退化模型折减。有限元中通常在每个积分点独立判断，因此同一层不同位置可能先后进入损伤状态，形成损伤带扩展。

从本构到损伤的建模闭环

综合来看，复合材料损伤有限元建模的流程可以概括为：定义单层工程常数 → 组装刚度矩阵 → 按铺层角变换 → 计算层合板 A/B/D 矩阵 → 求解整体位移与应变 → 反算各层应力 → 代入失效准则判断损伤起始 → 更新损伤变量与刚度 → 迭代至收敛。这一链条中任何一环出错，都会导致损伤起始位置和扩展路径偏离实际。因此，在建模初期用简单算例（如单层板拉伸、[0/90]ₛ 层合板弯曲）验证本构与变换的正确性，是保证后续复杂结构分析可信的必要步骤。