4. 有限元建模基础：单元类型选择、网格划分、边界条件与连接定义

做叶片复合材料结构分析，有限元建模是绕不开的一关。说实话，我见过太多人在这上面栽跟头——模型建得花里胡哨，一算就飞了，或者算出来结果自己都不敢信。今天咱们就聊聊建模的几个核心问题。

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

这是第一个要拍板的事。我个人习惯是：能选壳，就不选实体。为什么？

叶片是典型的薄壁结构，壁厚和长度比通常在1:100以上。壳单元用中面来模拟，计算量小，结果也够准。我做过一个项目，用壳单元算的模态和实测只差了3%。

但实体单元也不是没用。比如你要分析叶片根部的螺栓连接，或者铺层之间的层间应力，壳单元就抓瞎了。这时候得老老实实用实体。

对比项	壳单元	实体单元
适用场景	薄壁结构、整体变形、模态	局部细节、层间应力、接触
计算成本	低（一个节点6自由度）	高（一个节点3自由度，但层数多）
建模难度	低（中面提取即可）	高（需要处理厚度方向网格）
典型单元	S4R, S8R, S9R5	C3D8R, C3D20R

4.2 网格划分策略

网格这玩意儿，说白了就是「粗了不准，细了算不动」。怎么平衡？我有个原则：关键区域加密，非关键区域放水。

叶片上哪些是关键区域？

• 前缘和后缘：曲率大，应力集中，网格要密
• 根部过渡区：几何变化剧烈，至少加密2倍
• 铺层变化处：比如从单向布换到双轴布，那里应力会突变

非关键区域呢？比如叶片中部的腹板区域，应力梯度小，网格可以粗一些。我一般控制在20-30mm一个单元。

嗯，这里要注意：不要用四面体网格做壳单元分析。四面体太硬了，算出来的位移偏小，应力偏大。我吃过这个亏，后来老老实实切六面体。

4.3 边界条件设定

边界条件设错了，后面全白干。叶片分析常见的边界条件就两种：

1. 根部固支：所有自由度全约束。这是最常用的，模拟叶片装在轮毂上。
2. 离心力加载：给整个模型加转速，模拟旋转状态。

我遇到过一个问题：有人把根部约束设成了「铰支」，结果算出来的频率低了一大截。你想想看，叶片根部明明是螺栓拧死的，怎么会是铰支呢？

4.4 接触与连接定义

叶片不是一整块材料，它是铺层叠起来的。铺层之间怎么连接？这是个好问题。

我的做法是：

• 铺层之间：用共节点。说白了就是网格节点对齐，直接传力。这样最准，也最简单。
• 叶片和金属件（比如螺栓）：用接触对。主面选金属，从面选复合材料。摩擦系数设0.3左右。
• 胶接区域：用 cohesive 单元。这个比较高级，能模拟胶层脱粘。

为什么铺层之间不用接触？因为接触是非线性的，算起来慢。而且铺层之间本来就是固接，用共节点更合理。

好了，单元类型、网格、边界、连接，这四个东西搞明白了，有限元模型就立住了。剩下的就是调参数、跑计算、看结果。下一节咱们聊聊材料属性和铺层定义，那又是另一门学问。