第三章 经典层合板理论：基本假设、层合板刚度矩阵、ABD矩阵推导
各位工程师朋友，这一章咱们来啃一块硬骨头——经典层合板理论。说实话，我刚入行那会儿，看到ABD矩阵也是一头雾水。但干复合材料结构设计这行，这玩意儿就像吃饭的筷子，绕不开。
经典层合板理论，简称CLPT。它其实是从各向同性板的Kirchhoff理论延伸过来的。你想想看，单层板我们搞明白了，但实际结构都是多层叠在一起的。怎么把单层板的性能叠加成层合板的性能？这就是CLPT要解决的问题。
核心思想：层合板的整体行为，由各单层的材料属性、铺层角度和铺层顺序共同决定。说白了，就是“1+1≠2”的典型例子。
3.1 基本假设
任何理论都有假设前提。CLPT的假设，我个人总结为“四条铁律”：
- 直法线假设：变形前垂直于中面的直线，变形后仍为直线，且仍垂直于中面。这意味着层合板厚度方向的剪切变形被忽略了。
- 等应变假设：各层之间完美粘接，没有相对滑移。层间应变是连续的。
- 平面应力状态：厚度方向的正应力σz忽略不计。因为层合板通常很薄，面外应力很小。
- 小变形假设：位移远小于板厚，几何线性。
注意：我曾经在一个风电叶片项目中，忽略了第一条假设，结果算出来的变形和实测差了30%。后来才发现，厚叶片必须考虑剪切变形，CLPT不适用。所以，薄板用CLPT，厚板得用一阶剪切变形理论。
3.2 层合板的应变-位移关系
基于上述假设，层合板内任意一点的应变可以写成：
εx = εx⁰ + z·κx
εy = εy⁰ + z·κy
γxy = γxy⁰ + z·κxy其中：
- εx⁰, εy⁰, γxy⁰ —— 中面应变（膜应变）
- κx, κy, κxy —— 中面曲率（弯曲应变）
- z —— 该点到中面的距离
嗯，这里要注意：z的符号。中面以上为正，以下为负。我刚开始做铺层设计时，经常搞混正负号，结果算出来的ABD矩阵符号全反了。
3.3 单层板的应力-应变关系
对于第k层，在材料主方向（1-2坐标系）下，应力-应变关系为：
|σ1| |Q11 Q12 0 | |ε1|
|σ2| = |Q12 Q22 0 | |ε2|
|τ12| | 0 0 Q66| |γ12|其中Qij是缩减刚度系数，由工程常数决定：
Q11 = E1 / (1 - ν12·ν21)
Q22 = E2 / (1 - ν12·ν21)
Q12 = ν12·E2 / (1 - ν12·ν21)
Q66 = G12但实际铺层时，各层的材料主方向往往和整体坐标系（x-y）有个夹角θ。这时候就需要坐标变换。我记得第一次手动推导变换矩阵时，推了整整两页纸，最后发现一个符号错了，全部重来。
变换后的刚度矩阵为：
|σx| |Q̄11 Q̄12 Q̄16| |εx|
|σy| = |Q̄12 Q̄22 Q̄26| |εy|
|τxy| |Q̄16 Q̄26 Q̄66| |γxy|Q̄ij是变换后的刚度系数，它是铺层角度θ的函数。具体公式我就不列了，大家查教材都有。但我想强调一点：Q̄16和Q̄26的存在，意味着拉剪耦合。这是复合材料特有的现象，各向同性材料没有。
3.4 ABD矩阵的推导
好，重头戏来了。ABD矩阵是怎么来的？
层合板的内力（N）和内力矩（M）是各层应力沿厚度积分的结果：
Nx = ∫σx dz, Ny = ∫σy dz, Nxy = ∫τxy dz
Mx = ∫σx·z dz, My = ∫σy·z dz, Mxy = ∫τxy·z dz把应力用应变表示，应变再用中面应变和曲率表示，然后逐层积分。因为各层材料不同，积分要分段进行。最终得到：
|N| |A B| |ε⁰|
|M| = |B D| |κ |这就是大名鼎鼎的ABD矩阵方程。其中：
- A矩阵（拉伸刚度）：Aij = Σ(Q̄ij)k · (zk - zk-1)
- B矩阵（耦合刚度）：Bij = ½ · Σ(Q̄ij)k · (zk² - zk-1²)
- D矩阵（弯曲刚度）：Dij = ⅓ · Σ(Q̄ij)k · (zk³ - zk-1³)
这里zk是第k层上表面的z坐标。我建议你写代码时，把层序号从下往上编号，这样不容易乱。
个人经验：B矩阵如果非零，说明层合板存在拉弯耦合。也就是说，你拉它，它会弯。这在很多结构中是不希望出现的。所以对称铺层（B=0）是工程中最常用的。我设计飞机尾翼蒙皮时，第一件事就是检查B矩阵是否为零。
3.5 知识体系结构图
下面这张图，是我梳理的CLPT核心逻辑。你顺着箭头走一遍，就能把整个理论串起来：
3.6 工程应用中的注意事项
理论讲完了，说点实际的。我在做复合材料结构设计时，ABD矩阵有几个关键点要盯住：
| 矩阵项 | 物理意义 | 工程关注点 |
|---|---|---|
| A11, A22 | x、y方向的拉伸刚度 | 决定层合板的整体刚度，铺层越多越厚，A越大 |
| A66 | 剪切刚度 | ±45°铺层贡献最大，0°和90°贡献很小 |
| B矩阵 | 拉弯耦合 | 对称铺层时B=0，非对称铺层会导致翘曲 |
| D11, D22 | x、y方向的弯曲刚度 | 外层铺层对D贡献大，因为z³权重 |
| D66 | 扭转刚度 | ±45°铺层在外层时扭转刚度最大 |
避坑指南：我曾经设计一个无人机机翼，为了减重用了非对称铺层。结果固化后，机翼自己弯了5mm。这就是B矩阵不为零导致的固化变形。从那以后，但凡没有特殊要求，我都用对称铺层。
3.7 小结
经典层合板理论，说白了就是三件事：
- 搞清楚各层的应力-应变关系（Q̄矩阵）
- 沿厚度积分得到内力和内力矩
- 组装成ABD矩阵，建立载荷与变形的关系
ABD矩阵是复合材料结构设计的基石。你把它搞透了，后面的强度分析、失效判定、铺层优化，都是在这个基础上搭积木。我建议你找个简单的算例，手算一遍ABD矩阵，再用有限元软件验证一下。这样印象最深。
小技巧：手算时，先把各层的z坐标列清楚，再逐层计算A、B、D的贡献。用Excel表格做，比手算快得多，也不容易出错。我至今还保留着当年做的那个Excel模板。
公众号:蓝海资料掘金营，微信deep3321