4. 复合材料层合板经典理论回顾
断裂与损伤分析不能脱离层合板的整体力学响应。经典层合板理论（Classical Lamination Theory, CLT）是连接单层材料性能与层合板宏观刚度、应变的桥梁，也是理解层间应力与自由边效应的起点。本节回顾 CLT 的基本假设、ABD 矩阵的构成与物理意义，并说明为什么 CLT 无法直接给出层间应力，以及自由边效应如何成为分层损伤的诱因。
一、CLT 的基本假设与位移场
CLT 建立在 Kirchhoff 薄板假设之上：
- 板厚远小于面内尺寸，忽略横向剪切变形，即 γyz = γxz = 0；
- 厚度方向正应变忽略，即 εz ≈ 0；
- 直法线假设：变形前垂直于中面的直线，变形后仍保持直线且垂直于变形后的中面。
由此，板内任意一点的位移可写成中面位移与转角的关系：
u(x,y,z) = u₀(x,y) − z · ∂w₀/∂x
v(x,y,z) = v₀(x,y) − z · ∂w₀/∂y
w(x,y,z) = w₀(x,y)其中 u₀、v₀、w₀ 为中面位移。对应面内应变为中面应变与曲率的线性组合：
ε = ε⁰ + z · κ这里 ε⁰ 为中面应变向量，κ 为中面曲率向量。该式是 CLT 的核心运动学关系，也是后续 ABD 矩阵推导的基础。
二、ABD 矩阵：层合板刚度矩阵
单层在材料主方向上的应力—应变关系由折减刚度矩阵 Q 描述。将单层刚度转到层合板整体坐标系后得到 Q̄。对第 k 层，其面内力与弯矩贡献为：
N = Σ Q̄ₖ (zₖ − zₖ₋₁) ε⁰ + ½ Σ Q̄ₖ (zₖ² − zₖ₋₁²) κ
M = ½ Σ Q̄ₖ (zₖ² − zₖ₋₁²) ε⁰ + ⅓ Σ Q̄ₖ (zₖ³ − zₖ₋₁³) κ整理成标准形式：
[ N ] [ A B ] [ ε⁰ ]
[ M ] = [ B D ] [ κ ]其中：
- A 矩阵：面内刚度，反映层合板抵抗面内拉伸、压缩与剪切的能力；
- B 矩阵：耦合刚度，反映面内变形与弯曲变形的耦合，对称铺层时 B = 0；
- D 矩阵：弯曲刚度，反映层合板抵抗弯曲与扭转的能力。
ABD 矩阵的物理意义在于：它把铺层顺序、铺层角度、单层厚度与材料性能整合为一个整体刚度描述。对于断裂损伤分析，ABD 矩阵决定了载荷在层间的分配，也决定了各层应力水平，从而影响首层失效位置与损伤演化路径。
三、层间应力的来源与 CLT 的局限
CLT 假设直法线保持，因此横向剪应变 γxz、γyz 为零，横向正应力 σz 也被忽略。但在实际层合板中，相邻铺层刚度不匹配，面内应力在层间必须通过横向剪应力传递，才能满足层间应力连续与平衡条件。因此，层间应力 σz、τxz、τyz 是真实存在的，且往往是分层的直接驱动力。
CLT 给出的面内应力在层间是不连续的，无法直接满足三维平衡。要获得层间应力，通常采用两种后处理方式：
- 逐层积分平衡方程：从面内应力出发，沿厚度方向积分平衡方程，得到横向剪应力与正应力；
- 高阶剪切变形理论或三维有限元：放弃直法线假设，直接求解三维应力场。
对于工程断裂分析，常用第一种方式快速估算层间应力，再用断裂力学参数（如能量释放率 G）判断分层是否扩展。
四、自由边效应：分层的起点
自由边效应是指层合板在自由边界附近，层间应力急剧升高并形成边界层的现象。以 [0/90]ₛ 层合板受面内拉伸为例：0°层与 90°层横向刚度差异大，在自由边处，90°层试图收缩而 0°层约束它，导致层间剪应力 τxz 与层间正应力 σz 在边界附近迅速增大。该应力集中区域通常只有 1~2 倍板厚宽，但峰值可达到面内应力的数倍。
自由边效应的关键特征：
- 应力集中局限于边界层，远离自由边后迅速衰减；
- 铺层角度差越大、层间刚度不匹配越严重，效应越强；
- 层间正应力 σz 是 I 型分层的主要驱动力，层间剪应力则驱动 II/III 型分层；
- 自由边处的初始分层一旦形成，会沿层间界面扩展，导致层合板刚度下降与最终失效。
五、从 CLT 到断裂损伤分析的衔接
CLT 与 ABD 矩阵提供了层合板的整体刚度与逐层面内应力，是断裂损伤分析的输入。层间应力与自由边效应则揭示了 CLT 的边界，并指出分层损伤的优先位置。实际分析流程通常为：
- 用 CLT 计算层合板中面应变与曲率；
- 逐层计算面内应力 σx、σy、τxy；
- 通过平衡方程后处理得到层间应力 σz、τxz、τyz；
- 识别自由边附近的高应力区，判断分层起始位置；
- 结合断裂力学准则（如 G ≥ Gc）评估分层扩展。
这一链条把经典层合板理论与断裂损伤力学连接起来：ABD 矩阵决定应力分配，层间应力决定分层驱动力，自由边效应决定损伤起始位置。理解这三者的关系，是后续学习分层断裂准则、界面损伤模型与层合板失效分析的必要基础。