第二章 复合材料力学基础（一）：各向异性弹性力学基础、单层板的应力-应变关系、经典层合板理论入门
各位好，欢迎来到《叶片复合材料疲劳裂纹扩展分析》课程。今天咱们聊聊复合材料力学最基础的部分——各向异性弹性力学、单层板的应力-应变关系，以及经典层合板理论入门。
说实话，我刚接触复合材料那会儿，也被这些概念搞得头大。但后来发现，只要把基础打牢，后面分析疲劳裂纹扩展就顺理成章了。好，咱们开始。
2.1 各向异性弹性力学基础
先问大家一个问题：为什么复合材料力学比金属力学复杂那么多？
答案很简单——各向异性。金属材料，你从哪个方向拉，弹性模量基本一样。但复合材料不一样，纤维方向和非纤维方向，力学性能天差地别。
我当年做风电叶片分析时，第一次拿到材料数据表，看到E1、E2、G12这些符号，愣是懵了半天。后来才明白，这就是各向异性材料的“身份证”。
2.1.1 广义胡克定律
在复合材料力学中，应力-应变关系用广义胡克定律描述：
σ_i = C_ij * ε_j (i, j = 1, 2, ..., 6)其中C_ij是刚度矩阵，有36个分量。但别怕，实际没那么多。根据能量守恒和对称性，C_ij = C_ji，所以最多21个独立常数。
嗯，这里要注意：对于正交各向异性材料（比如单层板），独立常数只有9个。对于横观各向同性材料（比如单向纤维增强复合材料），更少，只有5个。
关键点：复合材料力学分析的第一步，就是确定材料属于哪种对称类型，然后才能正确使用本构关系。
2.1.2 工程常数与刚度矩阵的转换
实际工程中，我们更习惯用工程常数（E1, E2, ν12, G12等）来描述材料。这些常数和刚度矩阵之间可以互相转换。
举个例子，对于正交各向异性单层板，面内刚度矩阵Q的表达式为：
Q11 = E1 / (1 - ν12 * ν21)
Q22 = E2 / (1 - ν12 * ν21)
Q12 = ν12 * E2 / (1 - ν12 * ν21)
Q66 = G12其中ν21 = ν12 * E2 / E1，这个关系式我建议你记牢，我在做参数反演时经常用到。
个人经验：我建议你在做实验数据拟合时，先检查工程常数是否满足热力学稳定性条件：E1, E2, G12 > 0，且|ν12| < sqrt(E1/E2)。否则，你的材料数据可能有问题。
2.2 单层板的应力-应变关系
单层板是层合板的基本组成单元。说白了，就是一层纤维加一层基体。分析清楚了单层板，层合板的分析就水到渠成。
2.2.1 主方向坐标系下的应力-应变关系
在材料主方向坐标系（1-2坐标系，1方向沿纤维方向）下，单层板的应力-应变关系为：
|σ1| |Q11 Q12 0 | |ε1|
|σ2| = |Q12 Q22 0 | |ε2|
|τ12| | 0 0 Q66| |γ12|这个公式看起来简单，但背后有讲究。Q矩阵中的元素都是工程常数的函数。我在做叶片铺层设计时，经常需要根据纤维体积分数调整这些参数。
2.2.2 偏轴方向下的应力-应变关系
实际结构中，纤维方向往往和全局坐标系不一致。这时候就需要坐标变换。
假设纤维方向和x轴夹角为θ，那么偏轴刚度矩阵Q̄可以通过下式得到：
Q̄ = T⁻¹ * Q * T其中T是应变变换矩阵。展开后，Q̄的每个元素都是θ的函数。我记得有一次做叶片根部连接区分析，就因为角度搞错了，结果应力分布完全不对。后来花了整整两天排查，才发现是坐标变换矩阵写反了。
避坑指南：我曾经在坐标变换时搞混了应力变换和应变变换矩阵。记住：应力变换用Tσ，应变变换用Tε，两者不同。Tε = (Tσ⁻¹)ᵀ。这个关系一定要搞清楚。
2.3 经典层合板理论入门
经典层合板理论（CLT）是复合材料结构分析的核心工具。它把层合板看作一个整体，用中面应变和曲率来描述板的变形。
2.3.1 基本假设
CLT基于以下几个假设：
- 各层之间完美粘接，没有相对滑移
- 变形前垂直于中面的直线，变形后仍为直线且垂直于中面
- 层合板厚度方向的应变忽略不计
说白了，就是认为层合板在厚度方向不变形，且各层变形协调。这些假设在薄板分析中基本成立，但厚板或局部应力集中区域就要小心了。
2.3.2 本构关系：ABD矩阵
层合板的本构关系可以写成：
|N| |A B| |ε⁰|
|M| = |B D| |κ |其中：
- N：面内力合力
- M：弯矩合力
- ε⁰：中面应变
- κ：中面曲率
- A：拉伸刚度矩阵
- D：弯曲刚度矩阵
- B：耦合刚度矩阵
这里我要特别强调B矩阵。如果层合板不对称，B矩阵非零，就会出现拉弯耦合效应。你想想看，拉一下板子，它居然弯了。这在叶片设计中是要尽量避免的。
关键点：对称铺层可以消除拉弯耦合（B=0）。我在做叶片主梁设计时，一定会检查铺层顺序是否对称。不对称的铺层会导致叶片在挥舞方向产生附加弯矩，影响疲劳寿命。
2.3.3 层合板刚度计算步骤
实际计算层合板刚度时，我一般按以下步骤来：
- 确定每层的材料常数和厚度
- 计算每层在主方向下的Q矩阵
- 根据铺层角度，计算每层在全局坐标系下的Q̄矩阵
- 逐层积分，得到A、B、D矩阵
这个流程看起来简单，但手动计算很繁琐。我建议你写个小程序或者用Excel VBA来实现。我自己就写了一个，每次做铺层优化时直接调用，省了不少时间。
2.4 本章知识体系
为了让大家更直观地理解本章的知识结构，我画了一张图：
这张图展示了本章的核心逻辑：从各向异性弹性力学出发，建立单层板的应力-应变关系，再通过经典层合板理论，最终得到层合板的整体力学行为。每一步都是下一步的基础，环环相扣。
2.5 小结
好，咱们来捋一捋今天的内容：
- 各向异性弹性力学：复合材料力学的理论基础，核心是广义胡克定律和刚度矩阵
- 单层板应力-应变关系：主方向和偏轴方向的本构关系，坐标变换是关键
- 经典层合板理论：ABD矩阵是核心，对称铺层可以消除拉弯耦合
这些内容看起来有点枯燥，但确实是后续疲劳裂纹扩展分析的基础。我当年做叶片疲劳寿命评估时，就是因为层合板刚度算错了，导致裂纹扩展速率预测偏差很大。后来重新梳理了ABD矩阵的计算，才把问题解决。
嗯，今天就到这里。这些基础概念，建议大家多动手算几个例子，光看是记不住的。
实用建议：我建议你找一组简单的铺层数据（比如[0/90/0]对称铺层），手动计算一下ABD矩阵。算完之后，再用有限元软件验证一下。这样印象会深刻很多。