3. 复合材料力学基础（二）：层合板的刚度与柔度矩阵、铺层顺序对力学性能的影响、热残余应力简介
各位好，欢迎来到第三讲。上一讲我们聊了单层板的力学行为，那是复合材料的基础砖块。今天我们要把这些砖块砌成墙——也就是层合板。说白了，层合板就是把不同角度铺层的单层板叠在一起，用树脂固化成一个整体。这其中的力学行为，比单层板复杂得多。
我个人习惯把层合板理论看作是「搭积木」：先搞清楚每一块积木（单层板）的脾气，再研究它们叠在一起后整体怎么受力、怎么变形。今天我们就来拆解这个过程。
3.1 层合板的刚度矩阵与柔度矩阵
先问一个问题：一块由0°、90°、±45°铺层组成的层合板，你从外面拉它一下，它内部各层的应力分布是怎样的？
嗯，这个问题其实很关键。单层板有自己的刚度矩阵[Q]，但层合板作为一个整体，我们需要一个「等效」的刚度矩阵。这个矩阵通常用[A]、[B]、[D]三个子矩阵来描述。
- [A]矩阵——面内刚度矩阵：描述层合板在面内拉伸、压缩、剪切时的整体刚度。说白了，就是你从板子侧面拉它，它抵抗变形的能力。
- [B]矩阵——耦合刚度矩阵：这个最让人头疼。它描述了弯曲和拉伸之间的耦合效应。你拉它一下，它不光伸长，还会弯。我在项目中遇到过一块用[0/90]对称铺层的板子，拉的时候居然翘起来了，当时我还以为是夹具没装好，后来一查，原来是[B]矩阵不为零导致的。
- [D]矩阵——弯曲刚度矩阵：描述层合板抵抗弯曲变形的能力。
这三个矩阵的计算公式如下：
A_ij = Σ (Q_ij)_k * (z_k - z_{k-1})
B_ij = (1/2) * Σ (Q_ij)_k * (z_k² - z_{k-1}²)
D_ij = (1/3) * Σ (Q_ij)_k * (z_k³ - z_{k-1}³)其中，k代表第k层，z_k是第k层底面的坐标。你想想看，这其实就是把每一层的刚度按厚度方向积分。
有了[A]、[B]、[D]矩阵，层合板的本构关系就可以写成：
[N] = [A] [ε⁰] + [B] [κ]
[M] = [B] [ε⁰] + [D] [κ]这里[N]是面内力合力，[M]是弯矩合力，[ε⁰]是中面应变，[κ]是中面曲率。
至于柔度矩阵，其实就是刚度矩阵的逆矩阵。但在实际工程中，我更习惯用刚度矩阵直接算，因为柔度矩阵的物理意义不如刚度矩阵直观。当然，如果你要做变形分析，柔度矩阵会更方便。
重要概念：对称铺层（即铺层关于中面对称）可以消除[B]矩阵，即消除拉弯耦合。这是工程中最常用的设计手段之一。我建议你在设计层合板时，优先考虑对称铺层，除非有特殊需求。
3.2 铺层顺序对力学性能的影响
同样的铺层比例，比如都是50%的0°层和50%的90°层，但铺层顺序不同，力学性能会差很多。为什么会这样？
我举个例子。假设有两块板：
- 板A：[0/90/0/90]s（交替铺层）
- 板B：[0/0/90/90]s（集中铺层）
这两块板的[A]矩阵是一样的，因为总厚度和铺层比例相同。但[D]矩阵不同，因为铺层顺序影响了各层对弯曲刚度的贡献。
具体来说：
- 弯曲刚度：板B的0°层集中在表面，所以它的弯曲刚度更大。因为[D]矩阵中，远离中面的层贡献更大（z³项）。
- 层间应力：板A的交替铺层会导致更大的层间剪切应力，尤其是在自由边附近。我在做疲劳分析时，发现交替铺层的板子在自由边更容易出现分层，这就是原因。
- 热残余应力：铺层顺序也会影响固化后的热残余应力分布。这个我们后面会讲。
避坑指南：我曾经设计过一块风电叶片的主梁帽，为了追求弯曲刚度，把0°层全部堆在表面。结果固化后，板子翘曲得厉害，根本没法用。后来我改成了对称铺层，并把部分0°层分散到中间，问题才解决。所以，铺层顺序不是随便排的，要综合考虑刚度、热应力、工艺性。
铺层顺序的设计原则，我总结了几条：
- 对称铺层优先：避免拉弯耦合，减少翘曲。
- 避免相同角度层集中：比如不要连续铺6层0°层，否则层间应力会很大。
- ±45°层成对出现：这样可以消除剪切耦合效应。
- 表面层用±45°或织物：提高抗冲击和抗磨损能力。
3.3 热残余应力简介
复合材料是在高温下固化，然后冷却到室温的。这个过程中，由于纤维和基体的热膨胀系数不同，以及各铺层方向不同，会产生热残余应力。
你想想看，碳纤维的热膨胀系数几乎是零，而环氧树脂的热膨胀系数大约是50×10⁻⁶/°C。固化温度通常是180°C，冷却到室温20°C，温差160°C。这个温差下，树脂想收缩，但纤维不让它收缩，于是内部就产生了应力。
热残余应力的计算其实不复杂：
{σ_th} = [Q] * ({α} * ΔT - {ε⁰})其中{α}是热膨胀系数向量，ΔT是温差，{ε⁰}是中面应变（由整体平衡条件确定）。
对于层合板，热残余应力会导致：
- 翘曲变形：如果铺层不对称，冷却后板子会弯。
- 微裂纹：尤其是90°层，热残余拉应力可能超过基体的强度，导致横向裂纹。
- 影响疲劳寿命：热残余应力会叠加到外载荷引起的应力上，加速疲劳损伤。
注意：热残余应力不是一成不变的。在疲劳载荷作用下，它会逐渐释放（应力松弛），尤其是当温度接近玻璃化转变温度时。所以，在做疲劳分析时，不能简单地把热残余应力当作恒定值处理。
我记得有一次做叶片疲劳试验，发现某块试件的疲劳寿命比预期短了很多。排查了很久，最后发现是固化工艺出了问题，导致热残余应力比设计值大了30%。从那以后，我每次做设计都会要求工艺部门提供实际的固化曲线，而不是用理论值。
3.4 知识体系总览
下面这张图是我自己画的，把今天讲的内容串了起来。你可以看到，从单层板到层合板，核心就是[A]、[B]、[D]三个矩阵。铺层顺序和热残余应力，都是通过影响这三个矩阵来影响整体性能的。
这张图把今天的内容串起来了。从单层板的[Q]矩阵出发，通过层合板理论得到[A][B][D]矩阵，然后铺层顺序和热残余应力作为两个关键影响因素，最终决定了层合板的力学性能。
好了，今天就到这里。记住，层合板的设计不是简单的堆叠，而是一门平衡的艺术——刚度、强度、热应力、工艺性，都要兼顾。下一讲我们会深入疲劳裂纹扩展的机理，到时候这些基础知识都会用上。