3. 层合板宏观力学基础
前面我们把单层板当作正交各向异性材料处理，得到了它的工程弹性常数与偏轴刚度。但真实结构很少只用一层，而是把若干不同角度的单层按一定顺序叠合、固化成一个整体，这就是层合板。层合板宏观力学要回答的核心问题是：给定各单层的材料、厚度和铺层角，如何用一个整体的刚度关系描述它在面内载荷与弯曲载荷下的响应。经典层合板理论（CLT）给出的答案，就是著名的 A、B、D 三个刚度矩阵。
一、经典层合板理论的三条基本假设
CLT 建立在克希霍夫（Kirchhoff）薄板假设之上，可归纳为三条：
- 直法线假设：变形前垂直于中面的直线，变形后仍为直线且仍垂直于变形后的中面，即忽略横向剪切变形，因此 γxz = γyz = 0。
- 厚度不变假设：板厚方向的应变可忽略，εz ≈ 0，各层之间完全粘接、无滑移、无脱层。
- 面内位移沿厚度线性分布：由前两条可推出，板内任意一点的位移由中面位移与中面转角唯一确定。
由此得到层合板内任意 z 高度处的应变场：
εx(z) = εx⁰ + z·κx
εy(z) = εy⁰ + z·κy
γxy(z) = γxy⁰ + z·κxy其中 ε⁰ 是中面应变，κ 是中面曲率。这个线性关系是理解 A、B、D 的钥匙：面内应变与弯曲曲率通过 z 耦合在一起，而 z 的幂次决定了积分后出现 0 次、1 次、2 次三类刚度。
二、A、B、D 矩阵的定义与物理含义
把每层的偏轴刚度 Q̄ 沿厚度积分，得到层合板的内力—应变关系：
N = A·ε⁰ + B·κ
M = B·ε⁰ + D·κ
A_ij = Σ Q̄_ij·(z_k − z_{k−1})
B_ij = ½ Σ Q̄_ij·(z_k² − z_{k−1}²)
D_ij = ⅓ Σ Q̄_ij·(z_k³ − z_{k−1}³)三个矩阵的物理意义截然不同：
- A 矩阵（面内刚度）：描述中面拉伸/剪切与面内合力 N 的关系，单位是 N/m，只与各层 Q̄ 和厚度有关，与铺层顺序无关。
- D 矩阵（弯曲刚度）：描述曲率与弯矩 M 的关系，单位是 N·m，由于含 z³，对铺层位置极其敏感，把高模量层放得离中面越远，D 越大。
- B 矩阵（耦合刚度）：描述面内与弯曲之间的耦合，单位是 N。它含 z²，因此只有当铺层关于中面不对称时才不为零。B ≠ 0 意味着拉板会同时产生弯曲，或压弯会同时产生面内伸缩，这是层合板最反直觉、也最容易出问题的效应。
三、铺层顺序如何改变刚度
用一个具体示例说明。取 T300/环氧单层，E₁ = 135 GPa，E₂ = 10 GPa，G₁₂ = 5 GPa，ν₁₂ = 0.3，单层厚 0.125 mm，共 8 层，总厚 1 mm。比较三种铺层：
| 铺层方案 | B 矩阵 | A₁₁ (MN/m) | D₁₁ (N·m) | 典型行为 |
|---|---|---|---|---|
| [0₂/90₂]s（对称） | ≈ 0 | 约 78 | 约 6.5 | 拉伸不引起弯曲，工程首选 |
| [0/90/0/90]s（对称） | ≈ 0 | 约 78 | 约 4.2 | A 相同，D 因 0° 层靠近中面而下降 |
| [0₂/90₂/0₂/90₂]（不对称） | 非零 | 约 78 | — | 拉伸即翘曲，固化后易变形 |
可以看到：A 只取决于各角度的层数与厚度，与排列顺序无关；而 D 强烈依赖 0° 层离中面的距离。把 0° 层放到外侧（如 [0₂/90₂]s），D₁₁ 明显高于把 0° 层放在内侧的方案，抗弯能力更强。这正是工程上常把主要承载方向铺层放在外层的原因。
四、工程实现中的几个关键点
在有限元中建立层合板时，通常不直接输入 A、B、D，而是把每层的材料、厚度、角度按顺序赋给壳单元或实体单元，由求解器沿厚度积分得到等效刚度。但理解 A、B、D 能帮助我们判断结果是否合理：
- 若模型是对称铺层，拉伸载荷下不应出现明显面外位移；若出现，检查铺层定义是否真的对称。
- 若 B 非零，纯拉伸会产生弯曲，纯弯曲会产生面内伸缩，此时用壳单元必须开启耦合项，否则结果偏差很大。
- 比较不同铺层方案时，A 相同不代表性能相同，D 的差异可能使抗弯能力相差 50% 以上。
总结本节：CLT 用三条假设把三维问题降为二维，A、B、D 分别对应面内、耦合、弯曲刚度。A 由层数与厚度决定，D 由铺层位置决定，B 由对称性决定。铺层顺序不是细节，而是直接决定层合板是“拉弯独立”还是“拉弯耦合”的关键设计变量。掌握这三者的来源与量级，才能在有限元建模时正确设置铺层、合理解读位移与应力结果。