第二章：复合材料力学基础
各位好，我是老张。干风电叶片结构设计这行十几年了，今天咱们聊聊复合材料力学基础。说实话，这部分内容有点枯燥，但它是后面所有有限元分析的根基。你想想看，连材料的基本力学行为都搞不清楚，怎么去分析叶片铺层？
核心要点：复合材料与金属最大的区别在于——它的性能是方向相关的。金属材料你随便拉，各个方向性能差不多；复合材料呢？顺着纤维方向硬得像钢，垂直方向软得像塑料。这就是各向异性的本质。
2.1 各向异性材料与正交各向异性材料
先说说概念。各向异性材料，说白了就是材料在不同方向上的力学性能不一样。我刚开始接触这行时，总觉得这概念太抽象。直到有一次在车间，看到工人师傅铺层时，我拿边角料试了试——顺着纤维方向掰，费老大劲；垂直方向轻轻一折就断了。嗯，这下印象深了。
正交各向异性是各向异性的一个特例。它有三个互相垂直的弹性对称面。风电叶片用的单向预浸料，就是典型的正交各向异性材料。我习惯用1、2、3三个方向来标记：
- 1方向：纤维方向（纵向），最硬
- 2方向：面内垂直于纤维方向（横向），较软
- 3方向：厚度方向，通常不考虑面内受力
为什么强调这个？因为叶片铺层设计时，我们主要关心1-2平面内的行为。3方向的影响，在薄板假设下通常忽略不计。
我的经验：做叶片有限元分析时，材料属性定义最容易出错的就是方向。我曾经因为搞混了1方向和2方向，算出来的叶片刚度差了30%多。后来我养成了一个习惯——每次定义材料前，先在脑子里过一遍：纤维到底朝哪个方向？
2.2 单层板的应力-应变关系
单层板是复合材料层合板的基本单元。它的应力-应变关系，用广义胡克定律来描述。对于正交各向异性材料，在材料主方向坐标系下，关系式是这样的：
⎧ σ₁ ⎫ ⎡ Q₁₁ Q₁₂ 0 ⎤ ⎧ ε₁ ⎫
⎪ σ₂ ⎪ = ⎢ Q₁₂ Q₂₂ 0 ⎥ ⎪ ε₂ ⎪
⎩ τ₁₂⎭ ⎣ 0 0 Q₆₆⎦ ⎩ γ₁₂⎭其中Q矩阵是刚度矩阵，它的元素由工程常数决定：
| 刚度分量 | 表达式 | 物理意义 |
|---|---|---|
| Q₁₁ | E₁ / (1 - ν₁₂ν₂₁) | 纵向刚度 |
| Q₂₂ | E₂ / (1 - ν₁₂ν₂₁) | 横向刚度 |
| Q₁₂ | ν₁₂E₂ / (1 - ν₁₂ν₂₁) | 泊松耦合项 |
| Q₆₆ | G₁₂ | 剪切刚度 |
这里E₁是纵向模量，E₂是横向模量，G₁₂是面内剪切模量，ν₁₂是主泊松比。这些参数怎么来？要么做实验，要么查材料手册。我个人建议，关键参数一定要实测。手册上的数据只能做参考，实际批次的材料性能可能有10%-20%的波动。
注意：ν₁₂和ν₂₁不是独立的，它们满足关系：ν₁₂/E₁ = ν₂₁/E₂。这个关系在编程实现时很容易搞反，我踩过这个坑。
2.3 经典层合板理论（CLT）简介
CLT是复合材料结构分析的核心工具。它把层合板看作一个整体，用中面的变形来描述整个板的力学行为。说白了，就是假设：
- 直法线假设：变形前垂直于中面的直线，变形后仍然垂直于中面
- 等应变假设：各层之间完美粘接，没有相对滑移
- 平面应力状态：厚度方向的正应力忽略不计
基于这些假设，层合板的本构关系可以写成：
⎧ N ⎫ ⎡ A B ⎤ ⎧ ε⁰ ⎫
⎨ ⎬ = ⎢ ⎥ ⎨ ⎬
⎩ M ⎭ ⎣ B D ⎦ ⎩ κ ⎭这里：
- N：面内力（单位宽度的力）
- M：弯矩（单位宽度的力矩）
- ε⁰：中面应变
- κ：中面曲率
- A：拉伸刚度矩阵
- D：弯曲刚度矩阵
- B：耦合刚度矩阵（拉伸-弯曲耦合）
你可能会问：这个B矩阵有什么用？嗯，它代表了拉伸和弯曲的耦合效应。如果层合板不对称，比如[0/90]和[90/0]两种铺层，B矩阵就不为零。这意味着你拉它一下，它不光伸长，还会弯。叶片设计中，我们通常希望B矩阵尽量小，避免这种耦合。
关键公式：A、B、D矩阵的计算方法：
A_ij = Σ (Q_ij)_k * (z_k - z_{k-1})
B_ij = ½ Σ (Q_ij)_k * (z_k² - z_{k-1}²)
D_ij = ⅓ Σ (Q_ij)_k * (z_k³ - z_{k-1}³)其中k表示第k层，z是各层厚度方向的坐标。
我记得刚学CLT时，总觉得这些公式太繁琐。后来自己写了个小脚本，把铺层参数输进去，自动算A、B、D矩阵。你想想看，手工算一个20层的叶片铺层，那得算到猴年马月去？
2.4 知识体系总览
下面这张图，是我自己总结的本章知识脉络。建议你保存下来，后面学有限元分析时经常回来看看：
这张图把本章的三个核心内容串起来了。你从左边开始看：先理解材料为什么是各向异性的，然后掌握单层板的应力-应变关系，最后用CLT把各层组合起来。说白了，就是从微观到宏观的一个递进过程。
实用建议：刚开始学CLT时，别急着去推导公式。先拿一个简单的[0/90]s对称铺层，手算一遍A、B、D矩阵。算完你就发现，对称铺层的B矩阵确实为零。这个练习我让每个新人都做一遍，效果比看十遍书都好。
好了，这一章的内容就到这里。记住，复合材料力学的核心就三个字：方向性。后面做有限元分析时，你每定义一种材料，每设置一个铺层角度，都要问自己：这个方向对不对？
公众号:蓝海资料掘金营，微信deep3321