4. 经典层合板理论

经典层合板理论（Classical Lamination Theory, CLT）是层合板刚度分析的基础工具。它把单层板的应力—应变关系沿厚度方向逐层积分，最终把整块层合板的面向内力与中面应变、曲率联系起来。理解 CLT 的关键，不是记住三个矩阵的公式，而是弄清它们各自“管什么”：A 管面内拉压剪切，D 管弯曲扭转，B 管面内与弯曲之间的耦合。本节从基尔霍夫假设出发，推导 A、B、D 的积分形式，说明其物理意义，并给出适用边界。

4.1 基尔霍夫假设与位移场

CLT 建立在经典薄板理论上，其核心假设为：

• 直法线假设：变形前垂直于中面的直线，变形后仍为直线且保持垂直，即忽略横向剪切变形；
• 厚度不变：板厚方向的应变 ε_z≈0，各层之间无厚度方向拉压；
• 小变形：位移与应变满足几何线性关系；
• 层间完全粘结：相邻层界面位移连续，无滑移与脱层。

取板中面为 z=0，z 轴向上为正。设中面内位移为 u₀(x,y)、v₀(x,y)，中面挠度为 w₀(x,y)。由直法线假设，距中面 z 处的位移为：

u(x,y,z) = u0(x,y) - z · ∂w0/∂x
v(x,y,z) = v0(x,y) - z · ∂w0/∂y
w(x,y,z) = w0(x,y)

由此得到面内应变的线性分布：

εx = εx0 + z·κx
εy = εy0 + z·κy
γxy = γxy0 + z·κxy

其中 ε⁰ 为中面应变，κ 为中面曲率与扭率。这一步是 CLT 的“分水岭”：应变沿厚度是 z 的线性函数，因此整块板的变形可由中面量完全描述。

4.2 单层应力与层合板内力

第 k 层材料主方向与层合板坐标可能不一致，需用偏轴刚度矩阵 Q̄_k 描述其面内应力—应变关系：

σx       εx
σy  = Q̄k · εy
τxy      γxy

把 σ = Q̄_k(ε⁰ + zκ) 沿厚度积分，定义面内力 N 与弯矩 M：

N = ∫ σ dz ，M = ∫ σ·z dz

由于 Q̄_k 在每层内为常数，积分可逐层求和，得到 CLT 的核心方程：

N   A  B   ε0
   =       ·
M   B  D   κ

其中：

A = Σ Q̄k · (zk - zk-1)
B = ½ Σ Q̄k · (zk² - zk-1²)
D = ⅓ Σ Q̄k · (zk³ - zk-1³)

z_k 为第 k 层上表面坐标，z_k-1 为下表面坐标，求和覆盖全部 n 层。A、B、D 均为 3×3 对称矩阵，共 18 个独立分量。

4.3 A、B、D 的物理意义

A 矩阵描述层合板抵抗面内拉压与剪切的能力。它由各层偏轴刚度乘以层厚累加得到，因此只与各层的厚度和角度有关，与铺层顺序无关。A₁₁ 是 x 向拉压刚度，A₆₆ 是面内剪切刚度，A₁₆、A₂₆ 反映拉剪耦合，通常通过均衡铺层（如 ±θ 成对出现）使其为零。

B 矩阵是面内—弯曲耦合刚度。它来自积分中 z 的一次项，因此对铺层关于中面的对称性极为敏感。若层合板关于中面对称，正负 z 的贡献相互抵消，B=0，此时拉伸不引起弯曲、弯曲不引起面内伸缩。若铺层非对称，B≠0，会出现“拉弯耦合”：受拉时板会翘曲，受弯时中面会伸缩。工程上除非有特殊需求，一般优先采用对称铺层以消除 B 矩阵。

D 矩阵描述层合板抵抗弯曲与扭转的能力。它来自 z 的三次项，因此对离中面较远的层赋予更高权重——把高模量层放在外层能显著提高 D。D₁₁、D₂₂ 是主弯曲刚度，D₆₆ 是扭转刚度，D₁₆、D₂₆ 反映弯扭耦合。D 矩阵对铺层顺序非常敏感，是铺层优化设计的主要对象。

4.4 一个可手算的示例

取一块四层对称层合板 [0/90]_s，单层厚度 t=0.125 mm，总厚 0.5 mm。设单层正轴模量 E₁=140 GPa，E₂=10 GPa，G₁₂=5 GPa，ν₁₂=0.3。由于铺层对称，B=0。计算 A 与 D 时，只需分别对 0° 层和 90° 层的 Q̄ 矩阵按厚度加权。

0° 层的 Q̄ 即正轴 Q，90° 层的 Q̄ 相当于把 1、2 方向互换。以 A₁₁ 为例：

A11 = Q11(0°)·(2t) + Q11(90°)·(2t)
    = 141.3 GPa × 0.25 mm + 10.1 GPa × 0.25 mm
    ≈ 37.85 GPa·mm

D₁₁ 则需按 z³ 加权。0° 层位于最外与最内，对 D 贡献大；90° 层靠近中面，贡献小。这正是“外层材料决定弯曲刚度”的量化体现。实际设计中，若希望提高 D₁₁，可把 0° 层移到外层，形成 [0₂/90]_s 等铺层。

4.5 适用边界与常见误用

CLT 的适用边界必须明确：

• 薄板假设：板厚与最小跨度之比通常应小于 1/10。厚板中横向剪切变形不可忽略，需用一阶剪切变形理论（FSDT）或高阶理论。
• 小变形：大挠度问题需引入几何非线性，CLT 的线性应变—位移关系不再成立。
• 层间完全粘结：CLT 假设层间无滑移，因此无法预测自由边分层、层间应力集中等失效模式。这些需借助层间应力分析或三维有限元。
• 忽略横向正应力：σ_z≈0 在远离载荷与边界处成立，在集中载荷或厚截面处误差较大。
• 材料线弹性：CLT 本身不限制材料模型，但工程中常配合线弹性本构使用；若考虑塑性或损伤，需在单层本构中引入相应模型。

4.6 从 CLT 到铺层设计

CLT 给出的 A、B、D 是铺层设计的“刚度指纹”。设计流程通常是：先根据载荷确定所需的面内刚度 A 与弯曲刚度 D，再选择单层材料与厚度，最后通过调整铺层角度、顺序和对称性来逼近目标刚度。由于 A 与顺序无关、D 与顺序强相关，可以先定各层角度与层数满足 A，再通过交换层位微调 D，同时用对称与均衡条件把 B 和拉剪耦合压到零。

需要强调的是，CLT 是“刚度”层面的理论，它不直接给出强度。层合板强度需在得到各层应变后，逐层代入单层强度准则（如 Tsai-Wu、Hashin）判断。因此，CLT 是层合板分析的起点而非终点：它把三维问题降为二维中面问题，使后续逐层应力计算与失效判断成为可能。