第四章 经典层合板理论入门：单层板的本构关系、层合板的ABD矩阵、刚度与柔度矩阵、耦合效应

各位同行，欢迎来到第四章。

说实话，经典层合板理论（CLT）是咱们做铺层设计的“内功心法”。你铺层铺得再花哨，如果不懂ABD矩阵，那就像开车不看仪表盘——迟早要出事。我当年刚入行时，觉得这东西太理论，结果第一次做非对称铺层设计，产品一上夹具就扭成了麻花。嗯，从那以后，我再也不敢小看CLT了。

4.1 单层板的本构关系：从胡克定律说起

先聊单层板。说白了，单层板就是一层预浸料固化后的薄片。它很薄，但受力情况并不简单。

我们假设单层板处于平面应力状态。什么意思？就是厚度方向的应力σ₃、τ₂₃、τ₃₁都为零。这个假设在薄板分析里基本成立。

那么，单层板的正轴应力-应变关系可以写成：

| σ₁ |   | Q₁₁  Q₁₂  0   | | ε₁ |
| σ₂ | = | Q₁₂  Q₂₂  0   | | ε₂ |
| τ₁₂|   | 0    0    Q₆₆ | | γ₁₂|

这里的Q矩阵，就是缩减刚度矩阵。Q₁₁、Q₁₂这些系数，由材料的工程常数决定：

Q₁₁ = E₁ / (1 - ν₁₂ν₂₁)
Q₂₂ = E₂ / (1 - ν₁₂ν₂₁)
Q₁₂ = ν₁₂E₂ / (1 - ν₁₂ν₂₁) = ν₂₁E₁ / (1 - ν₁₂ν₂₁)
Q₆₆ = G₁₂

我个人习惯，每次拿到新材料数据，第一件事就是把Q矩阵算出来。为什么？因为后面所有计算都基于它。你想想看，如果E₁和E₂搞反了，那整个设计就全错了。

4.2 从单层到层合板：ABD矩阵的诞生

单层板搞清楚了，那多层板怎么办？

经典层合板理论的核心假设是：变形后，层合板的横截面仍然保持平面。说白了，就是各层之间没有滑移，变形协调。

基于这个假设，层合板内任意一点的应变可以写成：

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

其中ε⁰是中面应变，κ是中面曲率。z是距离中面的坐标。

然后，我们把每一层的应力沿厚度方向积分，就得到了层合板的合力与合力矩。结果就是著名的ABD矩阵：

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

这里：

• A矩阵（拉伸刚度）：描述面内拉力与中面应变的关系
• D矩阵（弯曲刚度）：描述弯矩与曲率的关系
• B矩阵（耦合刚度）：描述拉弯耦合效应——这是重点！

A、B、D矩阵的计算公式如下：

Aᵢⱼ = Σ (Q̄ᵢⱼ)ₖ · (zₖ - zₖ₋₁)
Bᵢⱼ = ½ Σ (Q̄ᵢⱼ)ₖ · (zₖ² - zₖ₋₁²)
Dᵢⱼ = ⅓ Σ (Q̄ᵢⱼ)ₖ · (zₖ³ - zₖ₋₁³)

注意，这里的Q̄是偏轴刚度矩阵。因为每一层的纤维方向可能不同，需要把正轴刚度通过坐标变换转到层合板的整体坐标系下。

4.3 刚度矩阵与柔度矩阵：一对“冤家”

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

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

其中：

• a = A⁻¹ + A⁻¹B(D - BA⁻¹B)⁻¹BA⁻¹
• b = -A⁻¹B(D - BA⁻¹B)⁻¹
• d = (D - BA⁻¹B)⁻¹

看着复杂是吧？其实用起来很简单。你给一个载荷N、M，就能算出变形ε⁰、κ。反过来，你给一个变形，也能算出需要的载荷。

我个人习惯，在做铺层优化时，经常看柔度矩阵的d项。d₁₁大，说明这个层合板在x方向容易弯曲；d₆₆大，说明容易扭转。这些信息对判断结构刚度特性非常直观。

4.4 耦合效应：B矩阵的“爱恨情仇”

耦合效应是层合板设计中最有意思，也最容易出问题的地方。

B矩阵非零，意味着：

• 拉弯耦合：面内拉力会产生弯曲变形
• 拉扭耦合：面内拉力会产生扭转变形
• 弯扭耦合：弯矩会产生扭转变形

为什么会这样？你想想看，如果层合板不对称，比如上面是0°层，下面是90°层。你拉它时，0°层很硬，90°层很软，上下变形不一致，自然就弯了。

消除耦合的最简单方法：对称铺层。只要铺层关于中面对称，B矩阵就为零。这是工程中最常用的方法。

但有时候，我们故意保留耦合。比如：

• 设计自适应变形结构，利用拉弯耦合实现“一拉就弯”
• 设计螺旋桨叶片，利用弯扭耦合实现被动变桨距

4.5 知识体系总览

下面这张图，是我自己总结的CLT知识框架。每次做设计前，我都会在脑子里过一遍这个流程：

4.6 实战中的几点提醒

最后，分享几个我在项目中积累的小经验：

1. 计算ABD矩阵时，注意单位统一。 我见过有人用mm算厚度，用Pa算模量，结果A矩阵的单位是N/mm，D矩阵的单位是N·mm，最后对不上。建议全部用国际单位制。
2. B矩阵的数值大小，可以反映耦合的强弱。 如果B矩阵元素比A矩阵小两个数量级，那耦合效应基本可以忽略。但如果只差一个数量级，就要小心了。
3. 柔度矩阵的d项，是判断弯曲稳定性的关键。 如果d₁₁和d₂₂相差太大，说明层合板在两个方向上的弯曲刚度严重不匹配，容易发生屈曲。
4. 不要迷信对称铺层。 对称铺层虽然消除了B矩阵，但有时候会限制设计自由度。比如你想实现特定的热变形，非对称铺层反而更合适。

好了，这一章的内容就到这里。记住，理论是死的，应用是活的。多算几个例子，多看看不同铺层下的ABD矩阵变化，慢慢你就会有感觉了。

---