3. 经典层合板理论回顾
在复合材料层合板的结构分析中，经典层合板理论（Classical Lamination Theory, CLT）是连接单层材料性能与整体结构响应的核心桥梁。它建立在克希霍夫（Kirchhoff）薄板假设之上：垂直于中面的直线在变形后仍保持直线且垂直于中面，忽略横向剪切变形。基于这一假设，层合板内任意一点的位移可由中面位移和转角唯一确定，从而将三维问题简化为二维问题。
一、从单层到层合板：ABD矩阵的物理含义
经典层合板理论的核心成果，是将层合板的内力（面内力 N 与弯矩 M）与中面应变 ε⁰ 和曲率 κ 联系起来，写成如下矩阵形式：
[ N ] [ A B ] [ ε⁰ ]
[ M ] = [ B D ] [ κ ]
其中 A、B、D 三个子矩阵分别称为拉伸刚度矩阵、耦合刚度矩阵和弯曲刚度矩阵。它们由各单层的偏轴刚度 Q̄ 和铺层几何位置共同决定：
A_ij = Σ Q̄_ij^(k) · (z_k − z_{k−1})
B_ij = ½ Σ Q̄_ij^(k) · (z_k² − z_{k−1}²)
D_ij = ⅓ Σ Q̄_ij^(k) · (z_k³ − z_{k−1}³)
式中 z_k 为第 k 层相对于中面的坐标，Q̄_ij^(k) 为该层在整体坐标系下的偏轴刚度。A 矩阵反映层合板抵抗面内拉伸/剪切的能力；D 矩阵反映其抵抗弯曲/扭转的能力；而 B 矩阵则描述面内变形与弯曲变形之间的耦合——这是层合板区别于各向同性板的最重要特征之一。
二、应力—应变关系的层级传递
层合板分析遵循“单层应力→单层应变→中面应变→层合板内力”的传递链。对第 k 层，其面内应力为：
σ_x^(k) Q̄_11 Q̄_12 Q̄_16 ε_x⁰ + z·κ_x
σ_y^(k) = Q̄_12 Q̄_22 Q̄_26 · ε_y⁰ + z·κ_y
τ_xy^(k) Q̄_16 Q̄_26 Q̄_66 γ_xy⁰ + z·κ_xy
可以看到，每层的应力不仅取决于中面应变，还与曲率 κ 和该层距中面的距离 z 有关。这正是弯曲问题中表层应力通常大于内层的原因。偏轴刚度 Q̄ 由单层正轴刚度 Q 经角度 θ 旋转得到，旋转公式中包含 cos⁴θ、sin⁴θ 及 cos²θsin²θ 等项，因此铺层角度对刚度的影响是非线性的。
三、铺层角度对刚度与耦合效应的影响
铺层角度是层合板设计中最灵活、最有效的调控手段。以典型的 T300/环氧单层为例，其正轴工程常数约为 E₁ = 135 GPa、E₂ = 9 GPa、G₁₂ = 5 GPa、ν₁₂ = 0.3。当铺层角度从 0° 变化到 90° 时，偏轴模量 E_x 会从约 135 GPa 急剧下降到约 9 GPa，变化超过一个数量级。这意味着仅通过改变角度，就能在很大范围内“定制”层合板的刚度。
更值得关注的是耦合效应。B 矩阵是否为零，取决于铺层是否关于中面对称：
- 对称层合板（如 [0/90]ₛ、[±45]ₛ）：B = 0，面内变形与弯曲变形互不耦合，拉伸时不产生弯曲，弯曲时不产生面内伸缩。
- 反对称层合板（如 [0/90]₂）：B ≠ 0，存在拉—弯耦合，拉伸会伴随弯曲或扭转。
- 非对称层合板：B ≠ 0，且可能出现拉—剪耦合（A₁₆、A₂₆ ≠ 0）和弯—扭耦合（D₁₆、D₂₆ ≠ 0）。
此外，A₁₆、A₂₆ 是否为零取决于铺层中 +θ 与 −θ 是否成对出现。例如 [±45]ₛ 中 +45° 与 −45° 层成对，A₁₆ = A₂₆ = 0；而 [0/45/90] 这类非平衡铺层则会产生拉—剪耦合，拉伸时伴随剪切变形。
四、工程实践中的判断与操作
在实际铺层设计中，工程师通常按以下逻辑使用 ABD 矩阵：
- 确定设计目标：需要面内刚度为主（如蒙皮受拉），还是弯曲刚度为主（如翼面抗弯），或两者兼顾。
- 选择铺层角度集合：常用 0°、±45°、90° 四种角度。0° 提供轴向刚度，90° 提供横向刚度，±45° 提供剪切刚度并改善损伤容限。
- 检查对称性：优先采用对称铺层使 B = 0，避免固化后因拉—弯耦合产生翘曲变形。若必须非对称，需在固化工艺中增加补偿或使用防翘曲模具。
- 检查平衡性：确保 +θ 与 −θ 层成对出现，使 A₁₆ = A₂₆ = 0，避免拉伸—剪切耦合导致的面内变形不协调。
- 计算 ABD 并验证：用上述积分公式计算 A、B、D，代入层合板本构方程求解中面应变和曲率，再回代到各层应力公式，检查最大应力是否低于许用值。
以 [0/±45/90]ₛ 为例，该铺层既对称又平衡，B = 0 且 A₁₆ = A₂₆ = 0，面内行为接近正交各向异性板，是工程中最常用的“准各向同性”或“均衡”铺层形式之一。而 [0/45/90] 非对称非平衡，虽然面内刚度可能相近，但固化后翘曲风险显著增大，且拉伸时会出现面内剪切变形，对连接孔和边缘应力分布不利。
理解 ABD 矩阵的构成逻辑，本质上就是理解“每一层的角度和位置如何通过积分累加，最终决定层合板的整体刚度与耦合行为”。这是后续铺层优化、失效分析和工艺补偿的共同基础。