4. 有限元建模基础：单元类型选择、网格划分策略、边界条件设置
做风电支撑结构疲劳分析，说白了就是跟三个东西打交道：单元类型、网格质量、边界条件。这三样搞不定，后面算得再花哨也是白搭。我这些年踩过的坑，十有八九都跟它们有关。
4.1 单元类型选择：梁、壳、实体，怎么选？
先说说单元类型。你想想看，一个风电塔筒，几十米高，壁厚才几十毫米。用实体单元去算？网格数量直接爆炸，算个疲劳得等到下个月。所以，选对单元类型，是效率与精度的平衡艺术。
4.1.1 梁单元（Beam Element）
梁单元，说白了就是一根线。它只关心轴线位置，截面属性用实常数定义。我个人习惯在初步设计阶段用梁单元，或者做整机模态分析时用。速度快，内存占用小。
- 适用场景：塔筒整体弯曲、整机模态、载荷传递路径分析。
- 典型单元：ANSYS中的BEAM188/BEAM189，Abaqus中的B31/B32。
- 注意：梁单元算不出局部应力集中。焊缝细节？别指望它。
4.1.2 壳单元（Shell Element）
壳单元是我做塔筒疲劳分析的主力。它兼顾了精度和效率。塔筒壁厚远小于直径，用壳单元模拟非常合适。我在项目中遇到过，用壳单元算出的热点应力，跟实测数据误差能控制在5%以内。
- 适用场景：塔筒筒体、法兰连接段、门洞加强区。
- 典型单元：ANSYS中的SHELL181/SHELL281，Abaqus中的S4R/S8R。
- 关键参数：注意单元积分点数量。减缩积分单元（如S4R）要注意沙漏控制。
4.1.3 实体单元（Solid Element）
实体单元，精度最高，代价也最大。我一般只在三种情况下用：螺栓连接细节、法兰根部过渡区、以及焊缝趾部应力集中分析。嗯，这里要注意，实体单元对网格质量要求极高，稍不注意就出现应力奇异。
- 适用场景：螺栓预紧力分析、焊缝细节、法兰连接非线性接触。
- 典型单元：ANSYS中的SOLID185/SOLID186，Abaqus中的C3D8R/C3D10。
- 避坑指南：我曾经用一阶六面体单元算螺栓，结果应力值一直不收敛。后来换成二阶单元，问题才解决。所以，实体单元尽量用二阶。
4.2 网格划分策略：不是越密越好
网格划分，很多人觉得越密越准。其实不然。网格太密，计算量上去了，但疲劳寿命结果可能只差1%。得不偿失。我的策略是：关键区域加密，非关键区域粗放。
4.2.1 全局网格控制
对于塔筒，我通常先设定一个全局网格尺寸。比如塔筒壁厚30mm，我设全局尺寸为100mm。这样整体网格不会太密。
- 长宽比：壳单元长宽比控制在5:1以内，实体单元控制在3:1以内。
- 翘曲度：壳单元翘曲角不超过5度。否则计算精度会下降。
- 雅可比：实体单元雅可比不低于0.7。低于这个值，单元畸变严重。
4.2.2 局部网格细化
关键区域必须细化。比如门洞周围、焊缝趾部、法兰连接处。我习惯用过渡网格，从细网格平滑过渡到粗网格。
4.2.3 网格质量检查清单
| 检查项 | 合格标准 | 我的经验值 |
|---|---|---|
| 长宽比 | ≤5:1 | 尽量≤3:1 |
| 翘曲角 | ≤5° | ≤2°更安全 |
| 雅可比 | ≥0.7 | ≥0.8 |
| 最小角度 | ≥20° | ≥30° |
| 最大角度 | ≤160° | ≤150° |
4.3 边界条件设置：别让模型飞起来
边界条件设置，是有限元分析中最容易出错的地方。我见过太多人，模型建得漂亮，网格画得完美，结果边界条件一设，算出来的结果完全不对。
4.3.1 基础约束
塔筒底部，通常是固定约束。但要注意：固定约束不是万能的。如果塔筒底部是法兰连接，你应该用螺栓预紧力+接触来模拟，而不是简单地把所有节点自由度全锁死。
- 固定约束：约束所有自由度（UX, UY, UZ, ROTX, ROTY, ROTZ）。适用于基础环或灌浆段。
- 简支约束：只约束平动自由度，释放转动。适用于某些简化模型。
- 弹簧约束：模拟基础刚度。我建议用实测或地勘数据来标定弹簧刚度。
4.3.2 载荷施加
疲劳分析中，载荷通常是时间序列。我习惯用多点约束（MPC）或刚性连接来施加风载荷和波浪载荷。直接加载在节点上？容易产生局部应力集中，影响疲劳结果。
4.3.3 接触与非线性
法兰连接、螺栓预紧、垫片接触，这些都需要设置接触对。我建议用罚函数法，摩擦系数取0.15-0.2。注意：接触分析收敛性差，需要耐心调整。
4.4 知识体系总览
下面这张图，是我总结的有限元建模核心逻辑。你把它记在心里，建模时就不会跑偏。
好了，以上就是有限元建模的基础三要素。单元类型选对了，网格策略跟上了，边界条件设准了，你的疲劳分析就成功了一大半。剩下的，就是耐心调试和验证了。