3. 层间应力解析方法

上一章我们建立了层间损伤的物理图像：分层往往不是从结构中部开始，而是从自由边、孔边、缺口或铺层突变处萌生。要预测它何时发生，必须回答一个更基础的问题——层间应力到底有多大、分布在哪里。本章围绕层间应力的解析方法展开，重点讲清自由边效应的成因、层间正应力与剪应力的分布规律，以及经典解析解与工程简化公式的适用边界。

一、为什么经典层合板理论"看不见"层间应力

经典层合板理论（CLT）建立在克希霍夫假设之上：垂直于中面的直线在变形后仍保持直线且长度不变。由此得到的位移场只含面内分量，应变沿厚度方向是分片线性的，应力则通过每层的偏轴刚度矩阵换算。问题在于，CLT 假设层与层之间是理想连续体，没有考虑层间应力的平衡。当相邻铺层刚度差异很大时（例如 0° 与 90° 相邻），两层在面内的泊松收缩和剪切变形不一致，界面处必然产生相互约束，这种约束只能由层间应力来传递。CLT 无法给出这些应力，因此必须引入更高阶的解析模型。

层间应力有三个分量：层间正应力 σ_z（剥离应力）、层间剪应力 τ_xz 与 τ_yz。其中 σ_z 对分层最危险，因为复合材料厚度方向强度极低；τ_xz 与 τ_yz 则驱动面内剪切型分层。理解它们的分布，是判断分层起点的关键。

二、自由边效应：应力集中的物理机制

考虑一块有限宽层合板，在自由边处（如拉伸试样的侧边），面内应力 σ_x 必须降为零，因为边界上没有外力。但层合板内部各层的 σ_x 分布并不均匀：0° 层承担大部分载荷，90° 层承担很少。在自由边附近，这种面内应力的重新分配只能通过层间应力来完成。于是，在自由边附近一个很小的区域内（通常为几个铺层厚度），层间应力迅速升高，形成边界层效应。

这个边界层的特征长度与铺层厚度、铺层顺序和材料刚度有关。对于典型的碳纤维/环氧层合板，边界层宽度约为 1～2 倍单层厚度。这意味着：如果单层厚度为 0.125 mm，自由边附近的应力集中区只有 0.125～0.25 mm 宽，但峰值应力可能达到 CLT 面内应力的 10%～30%，甚至更高。这就是为什么分层总是从自由边开始。

三、层间正应力与剪应力的分布规律

在自由边附近，层间应力的分布具有几个鲜明特征。第一，峰值位置：σ_z 的峰值通常出现在自由边处或距自由边极近的位置，而 τ_xz 的峰值可能略向内偏移。第二，衰减长度：应力从自由边向内迅速衰减，衰减长度约为 1～2 倍铺层厚度。第三，铺层顺序影响：将 90° 层分散到多个界面（如 [0/90/0/90]s 而非 [0/90₂]s）可以显著降低单个界面上的层间应力峰值，因为每个界面承担的约束减小了。第四，厚度方向位置：σ_z 在层合板中面附近通常较小，在靠近外表面的界面处较大；τ_xz 则在中面附近可能较大。

这些规律直接指导铺层设计：避免将大角度差铺层直接相邻，尽量采用分散铺层，在自由边附近增加铺层厚度或设置软化层，都是工程上常用的降低层间应力的手段。

四、经典解析解：Pipes-Pagano 与边界层理论

层间应力的经典解析解中，最具代表性的是 Pipes 和 Pagano 在 1970 年代提出的边界层解。他们针对 [0/90]s 层合板在面内拉伸下的自由边问题，通过引入层间应力函数，得到了层间应力的封闭形式解。其核心结论是：层间应力在自由边附近呈指数衰减，峰值与铺层刚度差成正比，与铺层厚度成反比。

Pipes-Pagano 解的形式可以简化为：

σz(x) ≈ σ0 · f(ξ) · (E1 - E2) / E1
τxz(x) ≈ σ0 · g(ξ) · (G12 - G23) / G12
其中 ξ = x / h，h 为单层厚度，f(ξ)、g(ξ) 为衰减函数

这个解虽然基于线弹性、小变形和各向异性层合板的假设，但它首次定量揭示了自由边效应的本质，至今仍是理解层间应力的基础。后续的 Wang-Pillai、Whitney-Pagano 等解进一步考虑了层间正应力与剪应力的耦合、多层界面相互作用以及非对称铺层的影响。

五、工程简化公式与适用边界

在工程实践中，直接求解上述解析解仍然繁琐。因此发展了一系列简化公式，用于快速估算层间应力峰值。常用的简化思路包括：

• 基于面内应力的比例估算：对于 [0/90]s 层合板，层间剪应力峰值可近似取为面内剪应力的 10%～20%，层间正应力峰值取为面内正应力的 5%～15%。具体系数取决于铺层顺序和材料体系。
• 基于刚度差的估算：σ_z,max ≈ k · σ₀ · (E₁ - E₂) / E₁，其中 k 为经验系数，通常取 0.1～0.3。
• 基于边界层宽度的估算：边界层宽度 δ ≈ 2h，峰值应力在 δ 范围内按线性或指数衰减。

需要特别强调的是，这些简化公式的适用边界非常明确：它们主要适用于线弹性、无损伤、自由边为直边的层合板，且铺层为对称均衡铺层。对于含孔、缺口、冲击损伤或非线性材料行为的情况，简化公式只能作为初步估算，必须借助有限元或实验验证。

六、从解析到工程判断：一个示例

假设有一块 [0/90]s 碳纤维/环氧层合板，单层厚度 h = 0.125 mm，承受面内拉伸应力 σ₀ = 500 MPa。材料参数为 E₁ = 135 GPa，E₂ = 10 GPa，G₁₂ = 5 GPa，层间剪切强度 S₁₃ = 60 MPa，层间拉伸强度 Z_t = 40 MPa。

用简化公式估算：σ_z,max ≈ 0.2 × 500 × (135 - 10) / 135 ≈ 92.6 MPa。这个值已经超过了层间拉伸强度 40 MPa，说明在 500 MPa 面内应力下，自由边处可能已经发生分层萌生。而 τ_xz,max ≈ 0.15 × 500 × (5 - 3.5) / 5 ≈ 22.5 MPa，低于剪切强度 60 MPa，剪切分层风险较低。因此，该层合板的主要失效模式是自由边处的剥离型分层。

这个示例说明：解析估算的价值不在于精确预测，而在于快速判断哪种层间应力分量是主导的、是否接近强度极限，从而指导后续的有限元建模和试验设计。实际工程中，通常将解析估算作为第一道筛选，再用有限元细化自由边附近的网格，最后用试验验证。

总结来说，层间应力的解析方法提供了从面内应力到层间应力的桥梁。自由边效应是层间应力集中的根源，σ_z 和 τ_xz 的分布规律决定了分层的起始位置和类型，经典解析解给出了定量框架，而工程简化公式则提供了快速估算的工具。理解这些方法的适用边界，是正确使用它们的前提。