4. 有限元分析基础：从思想到实战

各位好，我是老张。今天咱们聊聊有限元分析的基础。说实话，这章内容我讲了很多年，但每次备课都会发现新的体会。有限元这东西，说白了就是把一个连续体拆成小块，再拼回去算。听起来简单，但里面的门道可不少。

4.1 有限元法的基本思想

我记得刚入行那会儿，师傅跟我说：“有限元就是‘拆’和‘拼’的艺术。”这么多年下来，我越来越认同这个说法。

具体来说，有限元法分三步走：

1. 离散化——把叶片这个连续体，切成有限个小单元。单元之间通过节点连接。
2. 单元分析——对每个小单元，建立它的刚度矩阵。说白了就是：这个单元在受力时，会怎么变形。
3. 整体求解——把所有单元的刚度矩阵组装起来，加上边界条件，解出每个节点的位移，再反推应力应变。

你想想看，一个几十米长的叶片，内部铺层结构复杂得要命。直接算？不可能。但拆成几万个单元，每个单元内部假设是均匀的，这就好办多了。

4.2 单元类型选择：壳单元 vs 实体单元

做叶片铺层分析，最纠结的就是选什么单元。我个人的习惯是：

单元类型	适用场景	优点	缺点
壳单元（Shell）	薄壁结构，如叶片蒙皮、腹板	计算效率高，适合铺层定义	不能模拟厚度方向的细节应力
实体单元（Solid）	厚壁区域，如根部连接、粘接层	能模拟三维应力状态	计算量大，铺层定义麻烦
连续壳单元（Continuum Shell）	中等厚度，兼顾精度与效率	既有壳的效率，又有实体的精度	对网格质量要求高

嗯，这里要注意：叶片的主结构——蒙皮、腹板、主梁帽——我建议用壳单元。为什么呢？因为叶片是典型的薄壁结构，厚度远小于长度和宽度。壳单元正好擅长处理这种情况。

但根部连接区域，厚度可能达到几十毫米，这时候壳单元就不太准了。我曾经吃过这个亏——用壳单元算根部螺栓连接，结果应力集中系数偏大，搞得我重新建模。后来改用实体单元，结果就靠谱多了。

4.3 网格划分策略

网格划分，说白了就是“怎么切这块蛋糕”。切得好，算得准；切得烂，结果就是垃圾。

我的策略是：

• 全局控制：单元尺寸取叶片长度的1/100~1/200。比如50米的叶片，单元尺寸250~500mm。
• 局部加密：应力集中区域——比如后缘、根部过渡区、螺栓孔周围——网格加密到全局尺寸的1/5~1/10。
• 铺层方向：网格要顺着纤维方向走。你想想看，如果网格歪歪扭扭的，算出来的应力方向都是错的。
• 长宽比：壳单元的长宽比别超过5:1。我见过有人用长条形的单元，结果算出来的刚度偏软。

为什么会这样？因为有限元法的理论基础是“单元内位移插值”。单元形状越规则，插值精度越高。歪瓜裂枣的单元，插值误差就大。

4.4 边界条件与载荷施加

边界条件，就是告诉软件“这个零件是怎么固定的”。载荷，就是“它受了什么力”。这两样搞错了，后面全是白算。

边界条件：

• 叶片根部：通常固定所有自由度（固支）。但实际叶片是通过螺栓连接的，严格来说应该用弹簧约束模拟螺栓刚度。
• 对称边界：如果叶片是左右对称的（实际上不是，因为叶片有预弯），可以取一半模型，在对称面上加对称约束。
• 接触边界：叶片各部件之间（如蒙皮与腹板）用 tie 约束或接触定义。

载荷施加：

• 气动载荷：用压力面施加。我习惯把 CFD 算出来的压力分布，映射到结构网格上。
• 重力载荷：直接加体积力。别忘了叶片内部还有玻纤布、芯材，密度不一样。
• 离心力：旋转产生的。这个在叶片高速旋转时很重要，尤其是靠近根部的区域。
• 温度载荷：铺层固化时的残余应力。这个很多人忽略，但实际影响不小。

4.5 知识体系总览

下面这张图，是我自己总结的有限元分析知识框架。你把它存下来，每次做分析前看一眼，能少走很多弯路。

这张图把本章的核心内容串起来了。你从中心往外看，四个分支对应四个知识点。每个分支下面，又有具体的操作要点。做分析的时候，对着这张图检查一遍，基本不会漏东西。

---

公众号：蓝海资料掘金营，微信 deep3321