第三章：复合材料力学基础（二）——经典层合板理论（CLT）

各位好，欢迎来到复合材料叶片分析的第三讲。

上一章我们聊了单层板的应力-应变关系，那是基础。但实际工程中，谁会用单层板做叶片？没人。我们用的都是层合板——把好几层不同方向的单层板叠在一起，用树脂固化成一个整体。

这一章，我们就来啃这块硬骨头：经典层合板理论（CLT）。说实话，我刚入行时觉得这玩意儿就是一堆矩阵乘法，枯燥得很。直到有一次做叶片扭转分析，算出来的变形跟实测差了30%，我才老老实实回去翻CLT的推导。嗯，有些坑，不踩不知道。

3.1 基本假设——CLT的“三条军规”

任何理论都有假设。CLT的假设，说白了就是三条：

1. 直法线假设：变形前垂直于中面的直线，变形后仍然垂直于中面，且长度不变。这跟梁理论里的平截面假设是一个道理。
2. 层间完美粘接：各层之间没有滑移，也没有脱粘。变形是连续的。
3. 平面应力状态：每层都处于平面应力状态，即σ_z = τ_xz = τ_yz = 0。

你可能会问：第三条靠谱吗？叶片那么厚，怎么可能没有厚度方向的应力？

问得好。实际上，在叶片远离自由边的区域，这个假设是成立的。但在自由边附近，层间应力会变得很复杂，那是另一门学问了。我当年做的一个风电叶片根部连接区分析，就吃了这个亏——自由边效应导致的分层，差点让叶片在台架试验时就废了。所以记住：CLT适用于板的中部区域，边缘区域要小心。

3.2 运动学关系——从应变到位移

有了假设，我们就可以建立位移场了。

设中面的位移为u₀、v₀、w₀，那么板内任意一点（距离中面z处）的位移为：

u(x,y,z) = u₀(x,y) - z · ∂w₀/∂x
v(x,y,z) = v₀(x,y) - z · ∂w₀/∂y
w(x,y,z) = w₀(x,y)

这里的关键是那个“-z·∂w/∂x”项。它代表了弯曲引起的面内位移。你想想看，板弯曲时，上表面被压缩，下表面被拉伸，就是这个项在起作用。

对应的应变场为：

εₓ = εₓ⁰ + z · κₓ
εᵧ = εᵧ⁰ + z · κᵧ
γₓᵧ = γₓᵧ⁰ + z · κₓᵧ

其中：

• ε⁰ —— 中面应变（拉伸/剪切）
• κ —— 中面曲率（弯曲/扭转）

这个形式非常优雅：应变 = 中面应变 + 曲率 × 厚度坐标。线性分布，简单明了。

3.3 本构关系——从应变到应力

有了每层的应变，我们就可以用上一章讲的偏轴刚度矩阵Q̅来计算每层的应力了。

对于第k层：

σₖ = Q̅ₖ · εₖ

注意，这里的εₖ是总应变，包含了中面应变和弯曲应变两部分。而且，由于各层的铺层角度不同，Q̅ₖ也不同，所以应力沿厚度方向是分段线性、不连续的。

这一点很重要。我见过不少新手直接把各层的应力平均，结果算出来的强度完全不对。记住：应力在层间是跳跃的，因为材料方向变了。

3.4 合力与合力矩——从应力到力

应力是微观量，我们工程师关心的是宏观的力和力矩。

把各层的应力沿厚度积分，就得到：

• 合力 Nₓ, Nᵧ, Nₓᵧ （单位：N/m）
• 合力矩 Mₓ, Mᵧ, Mₓᵧ （单位：N·m/m）

积分过程我就不手写了，直接给结果——这就是CLT最核心的ABD矩阵：

| N |   | A  B |   | ε⁰ |
|   | = |      | · |    |
| M |   | B  D |   | κ  |

展开写就是：

Nₓ     A₁₁ A₁₂ A₁₆   B₁₁ B₁₂ B₁₆     εₓ⁰
Nᵧ  =  A₁₂ A₂₂ A₂₆   B₁₂ B₂₂ B₂₆  ·  εᵧ⁰
Nₓᵧ    A₁₆ A₂₆ A₆₆   B₁₆ B₂₆ B₆₆     γₓᵧ⁰
─────────────────────────────────────
Mₓ     B₁₁ B₁₂ B₁₆   D₁₁ D₁₂ D₁₆     κₓ
Mᵧ  =  B₁₂ B₂₂ B₂₆   D₁₂ D₂₂ D₂₆  ·  κᵧ
Mₓᵧ    B₁₆ B₂₆ B₆₆   B₁₆ B₂₆ D₆₆     κₓᵧ

这三个子矩阵的物理意义：

• A矩阵（面内刚度）：抵抗拉伸和剪切
• D矩阵（弯曲刚度）：抵抗弯曲和扭转
• B矩阵（耦合刚度）：连接拉伸和弯曲——这是层合板独有的特性

3.5 铺层顺序的影响——别小看顺序

同样的铺层角度和厚度比例，只是改变叠放顺序，结果可能天差地别。

举个例子：

铺层方案	铺层顺序（从下到上）	B矩阵	弯曲性能
方案A	[0/90/0/90]	非零	弯扭耦合
方案B	[0₂/90₂]	非零（更大）	耦合更强
方案C	[0/90]ₛ（对称）	零	无耦合

看到了吗？对称铺层（方案C）的B矩阵为零，拉伸和弯曲互不干扰。而非对称铺层（方案A、B）则存在耦合。

我曾经遇到一个案例：某型无人机机翼，设计时用了非对称铺层来产生弯扭耦合，结果加工时工人把顺序搞反了，导致机翼在飞行中反向扭转，差点失控。所以，铺层顺序不是小事，图纸上怎么写，工艺上就得怎么做。

3.6 刚度矩阵与柔度矩阵的互逆

有了ABD矩阵，我们可以求逆得到柔度矩阵：

| ε⁰ |   | a  b |   | N |
|    | = |      | · |   |
| κ  |   | bᵀ d |   | M |

其中a、b、d分别是A、B、D的逆矩阵的对应子块。注意，不是简单的A⁻¹、D⁻¹，因为B矩阵的存在，整个6×6矩阵的逆需要整体计算。

柔度矩阵在工程中很有用。比如，你给叶片施加一个已知的弯矩M，想算曲率κ，直接用d矩阵就行，不用去解方程组。

3.7 知识体系框架

下面这张图，是我自己总结的CLT知识脉络，帮你理清思路：

3.8 小结

这一章我们走完了CLT的完整推导路径：

• 从假设出发，建立运动学关系
• 通过本构关系得到各层应力
• 积分得到合力/合力矩
• 最终浓缩为ABD矩阵

CLT是复合材料结构分析的基石。后面的有限元实现、失效分析、优化设计，全都建立在这个框架之上。我个人建议，如果你刚开始学，不妨亲手算一个简单的[0/90]ₛ层合板的ABD矩阵，用Excel或者Python都行。算一遍，比看十遍都管用。