2、复合材料力学基础（上）：各向异性弹性力学基本方程、单层板的应力-应变关系、正交各向异性材料的工程常数
各位同学，欢迎来到复合材料叶片动态载荷响应模拟课程。
今天咱们聊点硬核的——复合材料力学基础。说实话，这部分内容我当年学的时候也觉得枯燥，一堆张量、矩阵，看着就头疼。但后来做项目才发现，这些基础不牢，后面分析叶片受力、预测失效，根本无从下手。
所以，咱们耐着性子，把这块硬骨头啃下来。
2.1 各向异性弹性力学基本方程
先问大家一个问题：为什么复合材料叶片不能用各向同性材料那套理论？
很简单。你想想看，金属材料，比如铝合金，你往哪个方向拉，它的弹性模量基本一样。但复合材料不一样，纤维方向硬得像钢筋，垂直纤维方向软得像塑料。这种性质，我们叫它各向异性。
描述这种材料的弹性行为，得从广义胡克定律说起。写成矩阵形式，就是：
σ_i = C_ij * ε_j (i, j = 1, 2, ..., 6)这里σ是应力，ε是应变，C是刚度矩阵。注意，下标从1到6，分别对应三个正应力/正应变和三个剪应力/剪应变。
我个人习惯把C矩阵想象成一个“黑箱”。你输入应变，它输出应力。这个黑箱里有多少个独立参数？
理论上，C_ij有6×6=36个元素。但根据应变能密度函数的对称性，C_ij = C_ji，所以独立参数减少到21个。这就是完全各向异性材料，也是最一般的情况。
我在项目中遇到过一种编织复合材料，它的铺层角度很复杂，分析时就得用这21个参数。说实话，测全这21个参数，实验成本高得吓人。所以实际工程中，我们通常做简化。
核心要点：各向异性弹性力学的核心，就是广义胡克定律。21个独立弹性常数，是描述材料最一般情况下的“完全体”。
2.2 单层板的应力-应变关系
好，咱们把问题聚焦一下。复合材料叶片，通常是由一层一层的单层板（也叫铺层）叠起来的。每一层，我们都可以看成是正交各向异性材料。
什么叫正交各向异性？说白了，就是材料有三个互相垂直的弹性对称面。比如单向带，纤维方向是1方向，垂直纤维方向是2方向，厚度方向是3方向。这三个方向，力学性能各不相同，但各自对称。
对于这种材料，刚度矩阵C会简化成什么样子？
嗯，这里要注意。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|你看，非零元素都集中在主对角线附近，像个“带状矩阵”。这就是正交各向异性材料的典型特征。
我曾经犯过一个错，把C13项给忽略了，结果算出来的层间应力完全不对。后来查了半天，才发现是这里出了问题。所以，千万别小看这些耦合项，它们在某些工况下（比如弯曲）影响很大。
避坑指南：在ABAQUS或ANSYS中定义复合材料时，一定要确认材料主方向（1、2、3方向）与你的全局坐标系是否对齐。我曾经因为坐标系没对齐，算出来的叶片变形方向都是反的。
2.3 正交各向异性材料的工程常数
刚才说的C矩阵，是刚度矩阵。但在工程上，我们更习惯用工程常数来描述材料。比如弹性模量E、泊松比ν、剪切模量G。这些常数物理意义明确，实验也好测。
对于正交各向异性材料，工程常数有9个：
- 3个弹性模量：E1, E2, E3（分别对应1、2、3方向）
- 3个泊松比：ν12, ν23, ν31（注意下标顺序，ν12表示1方向受力时，2方向的横向应变与1方向纵向应变的比值）
- 3个剪切模量：G12, G23, G31
它们与刚度矩阵C之间可以互相转换。我给大家一个常用的转换公式（以平面应力状态为例，即σ3=0）：
E1 = C11 - C12^2 / C22
E2 = C22 - C12^2 / C11
ν12 = C12 / C22
G12 = C66反过来，由工程常数求刚度矩阵：
C11 = E1 / (1 - ν12*ν21)
C22 = E2 / (1 - ν12*ν21)
C12 = ν12 * E2 / (1 - ν12*ν21)
C66 = G12注意，这里ν21 = ν12 * E2 / E1，这是由对称性决定的，不是独立常数。
下面这个表格，是我做项目时常用的一组典型碳纤维/环氧树脂单层板的工程常数，供大家参考：
| 工程常数 | 典型值 | 单位 |
|---|---|---|
| E1（纤维方向） | 135 | GPa |
| E2（垂直纤维方向） | 9.5 | GPa |
| G12 | 5.2 | GPa |
| ν12 | 0.30 | - |
你看，E1和E2差了十几倍。这就是复合材料“各向异性”的直观体现。设计叶片时，我们就是利用这种特性，让纤维沿着主受力方向，从而大幅提高结构效率。
重要提醒：工程常数不是随便用的。它们只适用于线弹性、小变形范围。如果材料进入非线性（比如基体开裂），这些常数就失效了。另外，温度、湿度对工程常数也有影响，做高精度分析时一定要考虑。
2.4 本章知识体系总览
为了让大家有个整体印象，我画了一张图，把本章的核心逻辑串起来：
这张图很清楚地展示了我们的逻辑：从最一般的各向异性弹性力学出发，引入正交各向异性假设，得到单层板的简化应力-应变关系，最后用工程常数来落地应用。每一步都有明确的物理意义和工程背景。
好了，这一章的内容就到这里。记住，这些基础概念是后续分析复合材料叶片动态响应的基石。下一章，我们会继续深入，讨论层合板的刚度与强度分析。