第二章 复合材料力学基础（一）：各向异性材料本构、单层板宏观力学分析、正交各异性材料的工程常数
各位同学，欢迎来到复合材料叶片热力耦合分析实战课程。今天咱们聊聊最基础、也最绕不开的一块内容——复合材料力学基础。
说实话，我刚开始接触复合材料时，也被那一堆张量、矩阵搞得头大。但干这行十几年下来，我慢慢发现，其实核心就几个概念：各向异性、正交各向异性、工程常数。搞懂了这些，后面的热力耦合分析才能站得住脚。
2.1 各向异性材料本构——从张量说起
先问大家一个问题：为什么复合材料比金属难分析？
金属材料，你拉它一下，它就在拉的方向上变形，其他方向基本不动。但复合材料不一样——你沿着纤维方向拉，它很硬；垂直纤维方向拉，它很软。更麻烦的是，你拉它的时候，它还会产生剪切变形。这就是各向异性。
各向异性材料的应力-应变关系，用广义胡克定律来描述：
σ_ij = C_ijkl · ε_kl这里C_ijkl是刚度张量，有81个分量。但别怕，由于对称性，实际上独立的分量只有21个。
关键点：各向异性材料有21个独立弹性常数。这是最一般的情况。我在做三维编织复合材料分析时，就遇到过需要完整21个常数的情况，那叫一个头疼。
写成矩阵形式更直观：
[σ] = [C] · [ε]
其中 [C] 是 6×6 的对称矩阵：
[C] =
[ C11 C12 C13 C14 C15 C16 ]
[ C22 C23 C24 C25 C26 ]
[ C33 C34 C35 C36 ]
[ C44 C45 C46 ]
[ C55 C56 ]
[ C66 ]嗯，这里要注意：矩阵是对称的，所以C_ij = C_ji。这也是为什么81个分量能缩减到21个。
2.2 正交各向异性材料——工程中遇到最多的
说实话，21个常数的各向异性材料在实际工程中很少见。我们遇到最多的，其实是正交各向异性材料。
什么叫正交各向异性？说白了，就是材料有三个互相垂直的对称面。比如单向纤维增强复合材料，沿着纤维方向、垂直纤维方向、厚度方向，这三个方向上的力学性能都不一样，但每个方向内部是均匀的。
正交各向异性材料的刚度矩阵长这样：
[C] =
[ C11 C12 C13 0 0 0 ]
[ C22 C23 0 0 0 ]
[ C33 0 0 0 ]
[ C44 0 0 ]
[ C55 0 ]
[ C66 ]看到了吗？非对角线上的C14、C15、C16这些项全变成0了。独立常数从21个降到了9个。
我的经验：在叶片复合材料分析中，我们通常把单层板当作正交各向异性材料来处理。但要注意，铺层之后整体变成各向异性了——这就是层合板分析的由来。
2.3 工程常数——工程师的语言
刚度矩阵C和柔度矩阵S虽然数学上很漂亮，但工程师们更喜欢用工程常数。为什么？因为直观啊！
正交各向异性材料的工程常数包括：
| 符号 | 名称 | 物理意义 |
|---|---|---|
| E1, E2, E3 | 弹性模量 | 三个主方向上的拉伸刚度 |
| G12, G23, G31 | 剪切模量 | 三个主平面内的剪切刚度 |
| ν12, ν23, ν31 | 泊松比 | 一个方向拉伸时另一个方向的收缩 |
柔度矩阵用工程常数表示为：
[S] =
[ 1/E1 -ν21/E2 -ν31/E3 0 0 0 ]
[ 1/E2 -ν32/E3 0 0 0 ]
[ 1/E3 0 0 0 ]
[ 1/G12 0 0 ]
[ 1/G23 0 ]
[ 1/G31 ]这里有个重要关系：ν12/E1 = ν21/E2，这是由对称性保证的。我曾经在项目里看到有人把这两个泊松比搞混了，结果算出来的变形完全不对。
避坑指南：我曾经犯过一个错误——直接用金属材料的分析方法来处理复合材料。结果呢？应力集中位置完全算错了。记住：复合材料的工程常数必须通过实验测定，不能随便从手册里抄一个值。
2.4 单层板宏观力学分析——平面应力状态
对于复合材料叶片来说，单层板的厚度远小于其他两个方向的尺寸。这时候我们可以做平面应力假设：σ3 = 0, τ23 = 0, τ31 = 0。
在这个假设下，正交各向异性单层板的应力-应变关系简化为：
[σ1] [Q11 Q12 0 ] [ε1]
[σ2] = [Q12 Q22 0 ] [ε2]
[τ12] [ 0 0 Q66] [γ12]这里的[Q]就是缩减刚度矩阵，它和[C]的关系是：
Q11 = C11 - C13²/C33
Q12 = C12 - C13·C23/C33
Q22 = C22 - C23²/C33
Q66 = C66用工程常数表示更直观：
Q11 = E1 / (1 - ν12·ν21)
Q12 = ν12·E2 / (1 - ν12·ν21) = ν21·E1 / (1 - ν12·ν21)
Q22 = E2 / (1 - ν12·ν21)
Q66 = G12核心要点：单层板的宏观力学分析，本质上就是用4个独立常数（E1, E2, ν12, G12）来描述材料的刚度特性。这4个常数，就是复合材料单层板的"身份证"。
我在做叶片铺层设计时，经常需要反复核对这4个常数。有一次，供应商提供的E1和E2数据明显不合理——E1/E2的比值超过了20，这在常规碳纤维/环氧体系中几乎不可能。后来一查，果然是测试方法有问题。
2.5 坐标变换——从材料主方向到全局坐标
单层板的纤维方向（1方向）和叶片的主方向（x方向）通常有个夹角θ。这时候就需要坐标变换。
变换后的刚度矩阵[Q̄]为：
[Q̄] = [T]⁻¹ · [Q] · [T]⁻ᵀ其中[T]是变换矩阵，和θ有关。展开后：
Q̄11 = Q11·cos⁴θ + 2·(Q12+2Q66)·sin²θ·cos²θ + Q22·sin⁴θ
Q̄12 = (Q11+Q22-4Q66)·sin²θ·cos²θ + Q12·(sin⁴θ+cos⁴θ)
Q̄22 = Q11·sin⁴θ + 2·(Q12+2Q66)·sin²θ·cos²θ + Q22·cos⁴θ
Q̄16 = (Q11-Q12-2Q66)·sinθ·cos³θ + (Q12-Q22+2Q66)·sin³θ·cosθ
Q̄26 = (Q11-Q12-2Q66)·sin³θ·cosθ + (Q12-Q22+2Q66)·sinθ·cos³θ
Q̄66 = (Q11+Q22-2Q12-2Q66)·sin²θ·cos²θ + Q66·(sin⁴θ+cos⁴θ)你想想看，本来单层板是正交各向异性的，但一转角度，Q̄16和Q̄26就不为0了——这意味着拉剪耦合效应出现了。这就是为什么复合材料叶片在承受拉伸载荷时，还会产生扭转变形。
我的习惯：在做有限元分析前，我会先用Python写个小脚本，把不同铺层角度下的[Q̄]算一遍，看看耦合项的大小。如果耦合项太大，说明这个铺层方案可能会引起较大的翘曲变形，需要调整。
2.6 本章小结
好了，这一章的内容就到这里。咱们捋一捋：
- 各向异性本构：21个独立常数，最一般的情况
- 正交各向异性：9个独立常数，三个对称面，工程中最常用
- 工程常数：E1, E2, G12, ν12 是单层板的"四大金刚"
- 单层板宏观分析：平面应力假设下，用[Q]矩阵描述刚度
- 坐标变换：角度一变，耦合效应就来了
这些内容看起来有点枯燥，但它们是后续热力耦合分析的基石。下一章我们会继续深入，聊聊层合板的刚度计算和热膨胀效应。到时候你会发现，今天学的这些全都能用上。
公众号:蓝海资料掘金营，微信deep3321