旋转体动力学分析
理论与物理
旋转体振动的特殊性
老师，旋转体的振动和静止结构有什么不同？
旋转体有3种特殊效应：
1. 离心力引起的刚度变化（旋转硬化） — 已在预应力模态中讨论过
2. 科里奥利力 — 在旋转坐标系中观察时，作用于振动物体的表观力
3. 陀螺效应 — 旋转物体方向改变时产生的效应
所以是增加了科里奥利力和陀螺效应。
运动方程：
- $[G]$ — 陀螺矩阵（反对称。与速度成比例的力）
- $[K_\sigma]$ — 几何刚度（离心力引起的旋转硬化）
- $[K_c]$ — 离心软化（spin softening）
$[G]$ 是反对称！它和阻尼矩阵 $[C]$ 在方程中的位置相同，但是非对称的，对吧？
$[G]$ 是反对称的（$[G]^T = -[G]$），它既不向系统输入能量，也不从系统提取能量。但会引起模态分裂（前行波和后退波频率不同的现象）。
坎贝尔图
绘制旋转速度 $\Omega$ 与 固有频率 $f$ 关系的坎贝尔图是旋转体振动分析的基本工具。
特征：
- 前行波（Forward Whirl） — 与旋转方向相同的振动。$\Omega$ 增加时频率上升
- 后退波（Backward Whirl） — 与旋转方向相反的振动。$\Omega$ 增加时频率下降
- 激励线（$f = n\Omega$） — 与不平衡激励（$n=1$）或高阶激励的交点即为共振点
前行波和后退波发生分裂…这就是陀螺效应吗？
是的。静止状态下频率相同的两个模态，随着旋转速度的增加会发生分裂。涡轮叶片或转子设计中需要准确预测这种分裂。
临界转速
临界转速（critical speed）对应激励线与固有频率线的交点处的旋转速度。尤其是不平衡激励（$n=1$）的交点最为重要。
基本原则是避开临界转速运行吗？
基本如此。API 617（压缩机）和API 612（蒸汽轮机）标准中，要求在临界转速上下设置±15%的裕度。
总结
我来梳理一下旋转体的振动。
要点：
- 3种特殊效应 — 离心力刚度变化、科里奥利力、陀螺效应
- 陀螺矩阵 $[G]$ 是反对称的 — 引起模态前行波/后退波分裂
- 坎贝尔图 — 旋转速度 vs. 固有频率。用于识别共振条件
- 临界转速 — 与激励线的交点。通过±15%的裕度来避开
- 需要进行复特征值分析 — 陀螺效应导致特征值为复数
陀螺效应的发现与应用历史
陀螺效应由法国物理学家莱昂·傅科于1852年在研究陀螺稳定性时发现。陀螺仪（Gyroscope）这个名字也是他命名的。旋转体角动量抵抗力矩的性质在20世纪初被应用于无线电罗盘，至今仍是INS（惯性导航系统）的核心技术。FEM中的陀螺效应表现为反对称陀螺矩阵[G]。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您是否有过急刹车时身体向前冲的经验？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是基于“缓慢施力，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 和 $\nabla \cdot \sigma$。拉弹簧时能感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长”的特性就是杨氏模量 $E$，它决定了刚度。常见的误解是：“刚度高=强度高”。不对，刚度是“不易变形”，强度是“不易破坏”，这是两个不同的概念。
- 外力项（载荷项）：体积力 $f_b$（如重力）和表面力 $f_s$（压力、接触力等）。可以这样理解——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓预紧力…这些都是外力。这里容易犯的错误是：搞错载荷方向。本想施加“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时，确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他弦。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力或弦的内部摩擦转化成了热能。汽车的减震器也是同样原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷·弹性模量也需统一为MPa/N系单位 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系单位时为tonne/mm³（钢的密度约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N，保持统一 |
数值解法与实现
FEM中的旋转体分析
如何用FEM进行旋转体的振动分析？
Nastran
```
SOL 107 $ 复特征值（含陀螺效应）
CEND
CMETHOD = 10
BEGIN BULK
RFORCE, 100, 1, , 100., 0., 0., 1. $ 旋转速度100 rad/s，绕z轴旋转
```
SOL 107用于包含陀螺矩阵的复特征值分析。RFORCE指定旋转速度和轴。
Abaqus
```
*STEP
*COMPLEX FREQUENCY, CORIOLIS=ON
20, ,
*DLOAD
element_set, CENTRIF, omega_squared, x, y, z
*END STEP
```
CORIOLIS=ON 启用陀螺效应。
Ansys
```
/SOLU
ANTYPE, MODAL
MODOPT, QRDAMP, 20
CORIOLIS, ON, , ON ! 启用科里奥利/陀螺效应
OMEGA, , , 100. ! 绕z轴旋转100 rad/s
SOLVE
```
在Ansys中是用CORIOLIS命令来启用陀螺效应啊。
Ansys的QRDAMP法是将实模态投影后，再求解包含阻尼/陀螺效应的复特征值。对于大规模模型很高效。
坎贝尔图的创建
步骤：
1. 设置10~20个旋转速度点（从0到最大运行速度）
2. 在每个速度下进行预应力模态+复特征值分析
3. 绘制旋转速度 vs. 固有频率图
4. 叠加绘制激励线（$f = \Omega, 2\Omega, 3\Omega, ...$）
5. 将交点识别为临界转速
Ansys Workbench里是自动化的吧。
Ansys Workbench的“Rotordynamics Analysis”可以自动生成坎贝尔图。它将旋转速度的参数化扫描和结果绘图集成在一起。
总结
我来梳理一下旋转体的数值方法。
要点：
- 需要进行复特征值分析 — 陀螺矩阵 $[G]$ 导致复特征值
- SOL 107（Nastran）, COMPLEX FREQUENCY CORIOLIS=ON（Abaqus）, QRDAMP+CORIOLIS（Ansys）
- 自动生成坎贝尔图 — Ansys最为便捷
- 分析多个不同旋转速度的工况 — 参数化研究
Campbell图的创建与临界转速的识别
Campbell图是以固有频率为纵轴、转速为横轴，绘制各转速下的固有频率曲线以及转速的整数倍（发动机阶次）直线的图。交点即为“临界转速（Critical Speed）”，设计上需要避开该转速运行或快速通过共振区。ISO 10816标准要求将坎贝尔图用于旋转机械设计标准中的临界转速评估。