第二章 复合材料力学基础（一）：各向异性弹性力学基础、应力-应变关系、正交各向异性材料的本构方程
2.1 为什么复合材料力学这么“麻烦”？
各位工程师朋友，咱们做复合材料的，第一道坎就是力学分析。你想想看，以前搞金属，一个弹性模量E、一个泊松比ν，基本就搞定了。但复合材料呢？纤维一个方向，基体另一个方向，铺层角度一变，性能全变了。
我刚开始接触复合材料时，就犯过一个低级错误。当时做一个碳纤维/epoxy的层合板拉伸试验，我拿着金属那套公式去算应力，结果跟实测值差了30%多。后来才发现，我忽略了材料在不同方向上的刚度差异。说白了，复合材料是“看人下菜碟”的——你沿着纤维方向拉，它很硬；你垂直纤维方向拉，它就软了。
这种特性，我们称之为各向异性。今天我们就来啃这块硬骨头。
2.2 各向异性弹性力学基础
先问大家一个问题：一个物体受力后，内部任意一点的应力状态，需要几个分量来描述？
答案是6个：3个正应力（σx, σy, σz）和3个剪应力（τxy, τyz, τzx）。对应的应变也有6个：3个线应变（εx, εy, εz）和3个剪应变（γxy, γyz, γzx）。
对于一般的各向异性材料，这6个应力和6个应变之间，通过一个6×6的刚度矩阵[C]联系起来：
σi = Cij · εj (i, j = 1, 2, ..., 6)这里Cij有36个分量。但别怕，根据能量守恒（应变能密度是正定的），Cij = Cji，所以独立分量最多21个。
21个！这比金属的2个多了整整一个数量级。你想想看，要测出这21个常数，得做多少试验？
核心概念：各向异性材料的弹性常数数量，取决于材料的对称性。对称性越高，独立常数越少。
2.3 应力-应变关系：从广义胡克定律说起
其实，复合材料力学并没有发明新的物理定律。我们用的还是胡克定律，只不过把它推广到了三维、各向异性的情况。
写成矩阵形式，就是：
|σx| |C11 C12 C13 C14 C15 C16| |εx|
|σy| |C12 C22 C23 C24 C25 C26| |εy|
|σz| |C13 C23 C33 C34 C35 C36| |εz|
|τyz| = |C14 C24 C34 C44 C45 C46| |γyz|
|τzx| |C15 C25 C35 C45 C55 C56| |γzx|
|τxy| |C16 C26 C36 C46 C56 C66| |γxy|嗯，看着挺吓人的。但实际工程中，我们很少遇到这种最一般的情况。为什么呢？因为绝大多数复合材料都有一定的对称性。
我个人习惯，拿到一种新材料，先判断它的对称类型。这能省下大量试验时间。
2.4 正交各向异性材料的本构方程
咱们最常遇到的，就是正交各向异性材料。比如单向带（UD层）、机织物层合板，都属于这一类。
什么叫正交各向异性？就是材料有三个互相垂直的弹性对称面。你想想看，一块单向碳纤维板，沿着纤维方向（1方向）、垂直纤维方向（2方向）、厚度方向（3方向），这三个方向上的力学性能都不同，但每个方向内部是均匀的。
这时候，刚度矩阵大大简化了。独立常数从21个降到了9个：
|σ1| |C11 C12 C13 0 0 0 | |ε1|
|σ2| |C12 C22 C23 0 0 0 | |ε2|
|σ3| |C13 C23 C33 0 0 0 | |ε3|
|τ23| = | 0 0 0 C44 0 0 | |γ23|
|τ31| | 0 0 0 0 C55 0 | |γ31|
|τ12| | 0 0 0 0 0 C66| |γ12|看到了吗？正应力与剪应力之间没有耦合了。也就是说，你拉材料的时候，不会产生剪切变形；你剪切材料的时候，也不会产生拉压变形。这在工程上是个好消息。
我的经验：在ABAQUS或ANSYS中定义复合材料时，你只需要输入9个工程常数：E1, E2, E3, ν12, ν13, ν23, G12, G13, G23。软件会自动帮你组装成上面的C矩阵。但要注意，这9个常数不是独立的，它们必须满足一定的热力学约束条件，否则计算会报错。
2.5 工程常数与刚度矩阵的转换
实际工作中，我们更习惯用工程常数（弹性模量、泊松比、剪切模量）来描述材料。对于正交各向异性材料，工程常数与刚度矩阵的转换关系如下：
| 工程常数 | 与刚度矩阵的关系 |
|---|---|
| E1 = 1 / S11 | S = C⁻¹（柔度矩阵） |
| E2 = 1 / S22 | ν12 = -S12 / S11 |
| G12 = 1 / S66 | ν21 = -S12 / S22 |
这里有个坑，我踩过。ν12和ν21不是独立的，它们满足互等关系：ν12 / E1 = ν21 / E2。你输入的时候，如果两个都随便填，算出来的刚度矩阵可能不正定，导致计算发散。
避坑指南：我曾经在做一个风电叶片的结构分析时，因为ν21输错了，导致屈曲载荷算出来偏小20%。后来查了三天才找到原因。所以，工程常数必须满足：E1, E2, E3, G12, G13, G23 > 0，且|νij| < sqrt(Ei / Ej)。
2.6 知识体系总览
为了让大家对本章内容有个整体把握，我画了一张图。这张图把各向异性、正交各向异性、以及工程常数之间的关系串起来了。
2.7 实际应用中的注意事项
好了，理论讲完了。咱们聊聊实际中怎么用。
第一，材料坐标系很重要。你测出来的E1、E2，都是相对于材料主方向的。铺层的时候，如果纤维方向跟你的坐标系对不上，就得做坐标变换。这个我们下一章会详细讲。
第二，厚度方向的性能别忽略。很多做薄板分析的工程师，觉得σz、τyz、τzx都是0，就把3方向忽略了。但如果你做的是厚板或者三维编织复合材料，厚度方向的剪切模量G13、G23对弯曲刚度影响很大。
第三，试验方法要选对。测E1用单向拉伸，测G12用±45°拉伸，测ν12用应变片。我见过有人用三点弯曲去测剪切模量，结果数据完全不能用。
总结一下：正交各向异性材料的本构方程，是复合材料力学分析的基石。你只要掌握了9个工程常数的物理意义和约束条件，就能正确建立材料模型。剩下的，就是坐标变换和层合板理论了。
嗯，今天就先聊到这儿。这些基础概念，你消化一下。后面我们会用这些知识，去分析层合板的刚度、强度，以及如何优化铺层顺序。