第三章 结构力学基础：材料力学、疲劳分析、模态分析、有限元方法简介
各位工程师朋友，大家好。今天我们来聊聊结构力学。
说实话，搞风电仿真，结构力学是绕不开的坎。叶片转起来，塔筒立在那，机舱里几十吨的设备——这些都得靠结构来扛。我刚开始做风电那会儿，总觉得流体仿真才是核心，结构嘛，差不多就行。结果有一次算出来的叶片变形量，跟实测差了快一倍。嗯，从那以后我再也不敢小看结构力学了。
核心观点：结构力学是风电多物理场仿真的骨架。没有它，你的流体、电磁、控制算得再漂亮，也是空中楼阁。
3.1 材料力学：风电结构的“基本功”
材料力学，说白了就是研究材料在外力作用下怎么变形、怎么破坏的。风电里最常用的几个概念，我给大家捋一捋。
3.1.1 应力与应变
应力是单位面积上的力，应变是变形量除以原长。公式很简单：
σ = F / A （应力 = 力 / 截面积）
ε = ΔL / L （应变 = 变形量 / 原长）我在项目中遇到过一个问题：塔筒连接法兰的螺栓，算出来的应力明明在许用范围内，可用了两年就断了。后来一查，是应力集中——螺栓孔边缘的应力比平均值高了好几倍。所以啊，光看平均应力是不够的，得关注局部。
3.1.2 弹性模量与泊松比
弹性模量E，就是应力-应变曲线的斜率。钢的E大约是210 GPa，玻璃钢的E就低得多，大概20-40 GPa。泊松比ν，描述的是材料在拉伸时横向收缩的程度，钢材大约0.3，复合材料可以到0.4以上。
个人经验：风电叶片用的复合材料，弹性模量是各向异性的——顺着纤维方向和垂直方向差很多。我建议做仿真时，一定要用正交各向异性材料模型，别偷懒用各向同性。
3.1.3 强度理论
风电结构承受的是多轴应力状态，不是简单的拉压。常用的强度理论有：
- 最大拉应力理论（第一强度理论）：适合脆性材料，比如铸铁
- 最大切应力理论（第三强度理论）：适合塑性材料，比如钢材
- 畸变能密度理论（第四强度理论）：工程上最常用，也叫von Mises屈服准则
你想想看，塔筒在风载和自重作用下，既有弯曲又有扭转，应力状态很复杂。我一般用von Mises等效应力来评估，简单可靠。
3.2 疲劳分析：风电的“隐形杀手”
风电结构最怕什么？不是一次性的强风，而是日复一日的循环载荷。叶片每转一圈，应力就循环一次。20年设计寿命，要转几千万次。这就是疲劳问题。
3.2.1 S-N曲线与疲劳寿命
S-N曲线描述的是应力幅值S与循环次数N的关系。钢的S-N曲线有个疲劳极限——低于某个应力幅，理论上可以无限循环。但复合材料没有明显的疲劳极限，这就麻烦了。
| 材料 | 疲劳极限（MPa） | 循环基数 |
|---|---|---|
| 结构钢Q345 | 约170 | 10^7 |
| 铝合金6061 | 约70 | 10^7 |
| 玻璃钢（叶片用） | 无明显极限 | — |
避坑指南：我曾经在叶片疲劳分析中，直接用材料的S-N曲线算寿命，结果偏危险。后来发现，实际结构有焊接、螺栓连接、几何突变，这些地方的疲劳强度比母材低很多。一定要考虑细节系数，比如焊接系数、尺寸系数、表面粗糙度系数。
3.2.2 疲劳累积损伤理论
风电承受的载荷是变幅的，不是恒幅。今天风大，明天风小，怎么办？用Miner线性累积损伤法则：
D = Σ (ni / Ni)
其中：
ni — 第i级应力幅的实际循环次数
Ni — 第i级应力幅对应的疲劳寿命
D — 累积损伤，当D ≥ 1时，结构失效这个公式看着简单，但实际用起来有坑。Miner法则假设损伤是线性累积的，但实际材料有载荷顺序效应——先大后小和先小后大，结果不一样。我一般会留点安全裕度，把D的限值取0.8-0.9。
3.3 模态分析：别让结构“共振”
模态分析，说白了就是找结构的“固有频率”和“振型”。风电最怕什么？共振。叶片转速产生的激振频率，如果跟塔筒的固有频率重合，那振幅会越来越大，很快就把结构搞坏了。
3.3.1 单自由度系统
最简单的模型：质量-弹簧-阻尼系统。固有频率：
ωn = √(k/m) （无阻尼固有频率）
fn = ωn / (2π) （固有频率，单位Hz）我记得刚入行时，有个老工程师跟我说：“塔筒的固有频率，一定要避开叶片通过频率的1P和3P。”1P是叶片转速频率，3P是三个叶片轮流通过塔筒的频率。这个经验我一直用到现在。
3.3.2 多自由度系统与有限元模态
实际风电结构是连续体，有无数个自由度。有限元法把结构离散成有限个单元，求解特征值问题：
[K]{φ} = ω² [M]{φ}
其中：
[K] — 刚度矩阵
[M] — 质量矩阵
{φ} — 振型向量
ω — 固有频率实用技巧：做模态分析时，我建议至少提取前10阶模态。但要注意，高阶模态的精度受网格密度影响很大。我曾经用粗网格算第5阶模态，频率差了15%。后来加密网格，才跟实测对上了。
3.4 有限元方法简介：仿真的“发动机”
有限元方法，简称FEM，是现代工程仿真的核心。它的思想很简单：把复杂结构拆成一个个小单元，每个单元用简单的方程描述，然后组装起来求解。
3.4.1 基本步骤
- 前处理：几何建模、网格划分、材料属性、边界条件
- 求解：组装刚度矩阵、施加载荷、求解方程组
- 后处理：查看应力、变形、模态等结果
你想想看，一个风电叶片可能有几十万个单元，每个单元有6个自由度（3个平移+3个转动），那总自由度就是上百万。求解这么大的方程组，没有计算机根本不可能。
3.4.2 单元类型选择
| 单元类型 | 适用场景 | 自由度 |
|---|---|---|
| 梁单元 | 塔筒、主轴 | 6/节点 |
| 壳单元 | 叶片蒙皮、机舱罩 | 6/节点 |
| 实体单元 | 轮毂、连接件 | 3/节点 |
注意：单元类型选错了，结果可能完全不对。我见过有人用实体单元算薄壁塔筒，结果弯曲刚度偏大，因为实体单元有“剪切闭锁”问题。薄壁结构，用壳单元更合适。
3.4.3 网格质量
网格质量直接影响计算精度和收敛性。几个关键指标：
- 长宽比：最好小于5，别超过10
- 扭曲角：小于60度
- 雅可比：大于0.7
我个人的习惯是，关键区域（比如螺栓孔、焊缝）网格加密，其他地方可以粗一些。但过渡要平滑，别突然从细网格跳到粗网格，那样会产生数值误差。
知识体系总览
下面这张图，是我自己总结的结构力学在风电仿真中的知识框架，大家可以参考一下：
这张图把四个核心模块串起来了。材料力学是基础，疲劳分析是寿命保障，模态分析是动态安全，有限元方法是实现工具。缺一不可。
好了，这一章的内容就到这里。结构力学的东西很多，但抓住这几个核心，你就能在风电仿真中站稳脚跟。下一章我们聊聊流体力学基础，那是另一个有意思的世界。
公众号:蓝海资料掘金营，微信deep3321