第二章 复合材料力学基础：各向异性材料、正交各向异性、单层板宏观力学分析、经典层合板理论
各位工程师朋友，大家好。我是你们这门课的老朋友。今天咱们要啃的这块骨头，是复合材料结构设计的“地基”——力学基础。说实话，我当年刚接触复合材料时，也被这些“各向异性”、“正交各向异性”搞得头大。但后来我发现，只要把几个核心概念理清楚，后面做轻量化设计就会顺风顺水。
这一章，我会带着大家从最基础的“各向异性”讲起，一路走到“经典层合板理论”。嗯，别被这些名词吓到，咱们一步步来。
2.1 各向异性材料 vs 正交各向异性材料
先问大家一个问题：为什么复合材料比金属难分析？
答案很简单——因为金属是各向同性的，你从哪个方向拉它，它的反应都一样。但复合材料不一样。你想想看，碳纤维沿着纤维方向强度极高，但垂直方向就弱很多。这就是典型的各向异性。
各向异性材料：材料在不同方向上的力学性能完全不同。比如单晶、木材、以及我们最关心的复合材料单层板。
正交各向异性材料：这是各向异性的一个特例。它有三个互相垂直的弹性对称面。说白了，就是材料在三个正交方向上的性能各自独立，但互不影响。
我在项目中遇到过一位同事，他拿各向同性的公式去算复合材料板的变形，结果差了30%以上。嗯，这就是没搞清楚材料属性的后果。
核心区别总结：
- 各向同性：1个弹性常数（E, ν）就能描述
- 正交各向异性：需要9个独立弹性常数
- 各向异性：需要21个独立弹性常数
2.2 单层板的宏观力学分析
单层板，说白了就是复合材料的最基本单元——一层纤维加一层基体。我们分析它时，通常把它看作正交各向异性材料。
我个人习惯用应力-应变关系来切入。对于正交各向异性的单层板，在材料主方向（1-纤维方向，2-垂直纤维方向，3-厚度方向）上，本构关系可以写成：
| ε1 | | 1/E1 -ν21/E2 -ν31/E3 0 0 0 | | σ1 |
| ε2 | | -ν12/E1 1/E2 -ν32/E3 0 0 0 | | σ2 |
| ε3 | | -ν13/E1 -ν23/E2 1/E3 0 0 0 | | σ3 |
| γ23| = | 0 0 0 1/G23 0 0 | | τ23|
| γ31| | 0 0 0 0 1/G31 0 | | τ31|
| γ12| | 0 0 0 0 0 1/G12| | τ12|这里E1、E2、E3是三个方向的弹性模量，G12、G23、G31是剪切模量，ν是泊松比。注意，ν12和ν21不是独立的，它们满足：ν12/E1 = ν21/E2。
我记得刚入行时，总记不住这个关系。后来我给自己编了个口诀：“大模量配小泊松，小模量配大泊松”。虽然不严谨，但挺好用。
2.3 经典层合板理论（CLT）
好了，单层板搞清楚了，那多层板怎么分析？这就是经典层合板理论（CLT）的用武之地。
CLT的核心假设其实很简单：
- 层间粘结完美，没有滑移
- 变形前垂直于中面的直线，变形后仍然垂直于中面（直法线假设）
- 厚度方向的应变忽略不计
基于这些假设，我们可以把层合板的合力和合力矩与中面应变和曲率联系起来。这就是著名的ABD矩阵：
| N | | A B | | ε⁰ |
| M | = | B D | | κ |其中：
- A矩阵：拉伸刚度矩阵，决定层合板的抗拉能力
- D矩阵：弯曲刚度矩阵，决定层合板的抗弯能力
- B矩阵：耦合刚度矩阵，描述拉伸和弯曲之间的耦合效应
这里我要特别提醒大家注意B矩阵。我曾经在设计一个对称层合板时，不小心把铺层顺序写反了，结果B矩阵非零，导致板一受拉就翘曲。嗯，那次教训挺深刻的。
避坑指南：
我曾经在计算非对称层合板时，忽略了B矩阵的影响，结果有限元分析结果和实验对不上。后来发现，B矩阵引起的耦合效应在薄板中尤其明显。所以，如果你的层合板不对称，一定要考虑B矩阵。
2.4 知识体系框架图
为了让大家更直观地理解这一章的知识结构，我画了一张图。嗯，这张图是我自己总结的，希望能帮你把零散的知识点串起来。
2.5 实战中的小技巧
最后，分享几个我在实战中总结的小技巧：
- 铺层顺序：尽量保持对称铺层，这样可以消除B矩阵的耦合效应。我一般会在设计初期就检查铺层对称性。
- 材料方向：在有限元软件中定义材料方向时，一定要和实际铺层方向一致。我曾经见过有人把0°和90°搞反了，结果分析结果完全不对。
- ABD矩阵计算：我建议用手算验证一下简单层合板的ABD矩阵，再交给软件。这样能帮你建立直觉。
个人经验：
我在做风电叶片轻量化设计时，经常用CLT快速估算层合板的刚度。虽然现在有限元软件很强大，但手算ABD矩阵能让你对设计有更深刻的理解。说白了，工具是死的，脑子是活的。
好了，这一章的内容就到这里。记住，力学基础是复合材料设计的“内功”，内功练好了，后面学轻量化设计才能事半功倍。