3. 力学基础（二）：经典层合板理论（CLT）、层合板的刚度与柔度矩阵、耦合效应分析

各位，欢迎来到第三讲。上一讲我们聊了单层板的应力应变关系，那是基础中的基础。今天我们要进入一个更实用、也更“复合材料特色”的领域——经典层合板理论，也就是大家常说的CLT。

说实话，我刚入行那会儿，觉得CLT就是一堆矩阵运算，枯燥得很。直到有一次，我设计的一个无人机机翼蒙皮，在测试时出现了明显的“弯扭耦合”，整个翼尖像麻花一样扭起来，我才意识到——嗯，这东西不学透，是要出大问题的。

3.1 经典层合板理论的基本假设

CLT不是凭空想出来的，它建立在几个关键假设上。你想想看，如果假设不成立，后面的计算全是白搭。

• 直法线假设：变形前垂直于中面的直线，变形后仍为直线，且仍垂直于中面。说白了，就是忽略横向剪切变形。
• 等应变假设：各层之间完美粘接，没有滑移。层间应变是连续的。
• 薄板假设：板的厚度远小于其他两个方向的尺寸。一般厚度比跨度小于1/10时适用。
• 小变形假设：位移远小于厚度，应变远小于1。

3.2 层合板的应变-位移关系

有了假设，我们就可以建立几何方程。设中面位移为 \( u_0, v_0, w \)，那么板内任意一点的位移可以写成：

u(x,y,z) = u_0(x,y) - z * ∂w/∂x
v(x,y,z) = v_0(x,y) - z * ∂w/∂y
w(x,y,z) = w(x,y)

对应的应变分量就是：

ε_x = ε_x^0 + z * κ_x
ε_y = ε_y^0 + z * κ_y
γ_xy = γ_xy^0 + z * κ_xy

其中，上标0表示中面应变，κ表示曲率或扭曲率。我个人习惯把κ称为“弯曲变形量”，这样好记。

3.3 层合板的合力与合力矩

单层板我们关心应力，但到了层合板层面，我们更关心合力和合力矩。为什么？因为层合板是由很多层叠起来的，每一层的应力都不一样，我们需要把它们“积分”成一个宏观的力。

定义如下：

N_x = ∫ σ_x dz      (单位宽度上的力，N/m)
N_y = ∫ σ_y dz
N_xy = ∫ τ_xy dz

M_x = ∫ σ_x * z dz   (单位宽度上的弯矩，N·m/m)
M_y = ∫ σ_y * z dz
M_xy = ∫ τ_xy * z dz

积分区间是从板的下表面到上表面。这里要注意，N_xy和M_xy虽然带下标xy，但它们不是扭矩，而是面内剪切力和扭转力矩。

3.4 经典层合板的本构方程——ABD矩阵

好了，重头戏来了。把应变-位移关系代入应力-应变关系，再积分，就得到了层合板的本构方程：

[ N ]   = [ A ] [ ε^0 ] + [ B ] [ κ ]
[ M ]     [ B ] [ ε^0 ] + [ D ] [ κ ]

写成矩阵形式更清楚：

| N_x |   | A11 A12 A16 | | ε_x^0 |   | B11 B12 B16 | | κ_x |
| N_y | = | A12 A22 A26 | | ε_y^0 | + | B12 B22 B26 | | κ_y |
| N_xy|   | A16 A26 A66 | | γ_xy^0|   | B16 B26 B66 | | κ_xy|

| M_x |   | B11 B12 B16 | | ε_x^0 |   | D11 D12 D16 | | κ_x |
| M_y | = | B12 B22 B26 | | ε_y^0 | + | D12 D22 D26 | | κ_y |
| M_xy|   | B16 B26 B66 | | γ_xy^0|   | D16 D26 D66 | | κ_xy|

这三个矩阵的含义：

• A矩阵（拉伸刚度）：描述面内拉伸/剪切与面内应变的关系。
• D矩阵（弯曲刚度）：描述弯曲/扭转与曲率的关系。
• B矩阵（耦合刚度）：描述拉伸与弯曲之间的耦合效应。这是复合材料特有的！

3.5 耦合效应分析——为什么B矩阵这么重要？

为什么会存在耦合？说白了，是因为层合板在厚度方向上不对称。比如，你把两层0°的碳纤维铺在一边，两层90°的铺在另一边，那这个板就是不对称的。一拉，它就会弯。

常见的耦合类型有三种：

1. 拉弯耦合（B11, B12, B22）：面内正应力引起弯曲变形。
2. 拉扭耦合（B16, B26）：面内剪切引起扭转变形。
3. 弯扭耦合（D16, D26）：弯曲引起扭转，或者扭转引起弯曲。

我给大家画个图，直观理解一下：

3.6 刚度矩阵与柔度矩阵的转换

在实际工程中，我们有时需要从力/力矩反求应变/曲率，这时候就要用到柔度矩阵。把ABD矩阵求逆即可：

| ε^0 |   = | a  b | | N |
| κ  |     | bT d | | M |

其中：

• a = A^-1 + A^-1 B (D - B A^-1 B)^-1 B A^-1 （面内柔度）
• d = (D - B A^-1 B)^-1 （弯曲柔度）
• b = - A^-1 B d （耦合柔度）

看着复杂，但别怕。现在有限元软件和计算工具都能自动算。我建议你手算一次简单的对称层合板（比如[0/90]s），感受一下矩阵运算的过程，这对理解物理意义很有帮助。

3.7 工程应用中的注意事项

讲到这里，我想分享几个实际项目中的经验：

• 对称铺层是首选：除非有特殊需求（比如需要预变形），否则尽量设计对称铺层。B矩阵为零，分析简单，制造变形也小。
• 注意D16和D26：即使是对称层合板，D16和D26也可能不为零。这意味着你施加纯弯矩时，板会产生扭转。我在做直升机旋翼时，就利用这个特性设计了特定的弯扭耦合，实现了被动减振。
• 温度效应：CLT也可以扩展到热-力耦合分析。记得加上热膨胀项，特别是碳纤维和环氧树脂的热膨胀系数差异很大，固化冷却后会产生残余应力。

3.8 小结

经典层合板理论是复合材料结构设计的基石。你掌握了ABD矩阵，就掌握了层合板的“力学身份证”。

记住三个关键点：

1. CLT基于直法线假设，适用于薄板
2. ABD矩阵完整描述了层合板的拉伸、弯曲和耦合行为
3. 耦合效应是复合材料特有的，设计时要善用或避免

下一讲我们会深入讨论强度理论，包括Tsai-Wu准则和Hashin准则。但今天的内容，我建议你找几个简单的铺层，手算一下ABD矩阵，再和有限元结果对比。相信我，这一关过了，后面的路就好走了。

---