第3章：铺层设计理论基础：经典层合板理论（CLT）简介、应力-应变关系、ABD矩阵的物理意义

各位工程师朋友，欢迎来到铺层设计的理论核心。说实话，经典层合板理论（CLT）是复合材料结构设计的基石。我当年刚入行时，觉得这东西太抽象，直到亲手算错一个铺层顺序，导致叶片在试验中提前失效……嗯，从那以后，我老老实实把CLT啃透了。

3.1 经典层合板理论（CLT）—— 到底在讲什么？

CLT说白了，就是把多层不同方向的单层板“叠”在一起，然后分析它在受力时的行为。你想想看，每一层纤维方向不同，有的0°、有的±45°、有的90°，它们共同工作时，变形必须协调——这就是CLT要解决的核心问题。

我个人习惯把CLT看作“复合材料的梁理论”。它基于三个基本假设：

• 直法线假设：变形前垂直于中面的直线，变形后仍为直线且垂直于中面。
• 等应变假设：各层之间完美粘接，没有滑移，应变沿厚度方向线性变化。
• 平面应力状态：每层厚度方向的应力忽略不计（σz = 0）。

3.2 应力-应变关系：从单层到层合板

我们先从单层板说起。每一层都有自己的材料主方向（1方向沿纤维，2方向垂直纤维）。在材料主坐标系下，应力-应变关系很简单：

σ1 = Q11·ε1 + Q12·ε2
σ2 = Q12·ε1 + Q22·ε2
τ12 = Q66·γ12

这里的Q矩阵是“缩减刚度矩阵”。我在项目中遇到过一个问题：有人直接用材料主方向的Q矩阵去算全局坐标下的应力，结果完全不对。记住，铺层方向不同，必须通过坐标变换把Q矩阵转到全局坐标系下。

坐标变换后的刚度矩阵记为Q̄（带横杠的Q）。公式我就不列了，你记住一点：Q̄是铺层角度θ的函数。0°层和90°层的Q̄完全不同。

3.3 ABD矩阵——层合板的“身份证”

好了，现在我们把所有单层叠起来。层合板整体受力时，中面会产生应变ε⁰和曲率κ。每一层的应变由这两部分叠加而成：

ε(z) = ε⁰ + z·κ

其中z是该层到中面的距离。然后对每一层积分，得到层合板的总内力和总弯矩。最终结果就是大名鼎鼎的ABD矩阵：

{N}   [A  B] {ε⁰}
{M} = [B  D] {κ }

这里：

• A矩阵（面内刚度矩阵）：描述拉伸和剪切变形。说白了，就是拉一块板子，它有多“硬”。
• D矩阵（弯曲刚度矩阵）：描述弯曲和扭转变形。板子抗弯能力怎么样，就看它。
• B矩阵（耦合刚度矩阵）：这是复合材料特有的！它把拉伸和弯曲耦合在一起。你拉一下板子，它不光伸长，还会弯——这就是B矩阵在作怪。

3.4 ABD矩阵的物理意义——一张图说清楚

为了让你更直观地理解，我画了一张图。这张图展示了层合板在三种基本载荷下的变形模式：

3.5 如何计算ABD矩阵？—— 一个简单例子

假设我们有一个两层板：一层0°（厚度0.5mm），一层90°（厚度0.5mm）。材料是T300/环氧，工程常数如下：

参数	E1 (GPa)	E2 (GPa)	ν12	G12 (GPa)
数值	135	9.5	0.30	5.2

计算步骤：

1. 先算每层的Q矩阵（材料主方向）。
2. 通过坐标变换得到Q̄矩阵（全局方向）。0°层Q̄ = Q，90°层需要旋转90°。
3. 对每层积分，累加得到A、B、D矩阵。

我习惯用Python写个小脚本，几行代码就搞定：

import numpy as np

# 材料常数
E1, E2, nu12, G12 = 135e9, 9.5e9, 0.30, 5.2e9
nu21 = nu12 * E2 / E1

# Q矩阵（材料主方向）
Q11 = E1 / (1 - nu12*nu21)
Q22 = E2 / (1 - nu12*nu21)
Q12 = nu12 * E2 / (1 - nu12*nu21)
Q66 = G12

Q = np.array([[Q11, Q12, 0],
              [Q12, Q22, 0],
              [0,   0,   Q66]])

# 坐标变换函数（角度转弧度）
def Q_bar(theta_deg):
    theta = np.radians(theta_deg)
    c, s = np.cos(theta), np.sin(theta)
    T = np.array([[c**2, s**2, 2*c*s],
                  [s**2, c**2, -2*c*s],
                  [-c*s, c*s, c**2 - s**2]])
    # 注意：这里用的是应变变换矩阵的逆
    return T @ Q @ T.T

# 两层：0°和90°，各0.5mm
layers = [(0, 0.5e-3), (90, 0.5e-3)]
z = -0.5e-3  # 中面在z=0，底层z=-0.5mm

A = np.zeros((3,3))
B = np.zeros((3,3))
D = np.zeros((3,3))

for theta, t in layers:
    Qb = Q_bar(theta)
    z0 = z
    z1 = z + t
    A += Qb * (z1 - z0)
    B += Qb * (z1**2 - z0**2) / 2
    D += Qb * (z1**3 - z0**3) / 3
    z = z1

print("A矩阵 (N/m):\n", A)
print("B矩阵 (N):\n", B)
print("D矩阵 (N·m):\n", D)

运行结果你会发现，B矩阵不为零——因为铺层不对称。这就是为什么我前面强调对称铺层的重要性。

3.6 小结

这一章我们聊了CLT的三个核心：基本假设、应力-应变关系、ABD矩阵。我个人觉得，理解ABD矩阵的物理意义比记住公式更重要。你想想看，当你拿到一个铺层方案时，能一眼看出它的A、B、D矩阵大概是什么样子，这才是真功夫。

我曾经带过一个新人，他算出来的B矩阵很大，我问他：“你觉得这个铺层拉一下会怎样？”他想了半天说：“会弯。”对了！这就是理解了物理意义。理论是为工程服务的，别把它当成数学题来做。

---

公众号:蓝海资料掘金营，微信deep3321