圆周运动和离心力 — ω²r加速度从涡轮·旋转机械CAE
1. 圆周运动和旋转机械CAE — 涡轮叶片的离心荷载
喷气发动机、蒸汽涡轮机、涡轮增压器、电机、离心泵——现代机械系统中的许多装置都包含高速旋转的零件。在旋转体的设计中，由旋转产生的离心荷载成为主导的机械荷载，加上热荷载、气体荷载、振动荷载等叠加作用。
圆周运动的物理虽然在高中学过，但要将其与喷气发动机的涡轮盘或风力发电转子的设计相联系，需要理解坐标系的概念、FEM荷载设置、旋转振动分析（Campbell线图）等多个方面。
2. 向心加速度和向心力
做匀速圆周运动的质点虽然速度大小不变，但方向不断改变，因此具有加速度。这个加速度始终指向圆的中心（向心方向），其大小为：
$$a_c = \frac{v^2}{r} = \omega^2 r$$其中 $v$ 是速度 [m/s]，$r$ 是半径 [m]，$\omega$ 是角速度 [rad/s]。对应的向心力为：
$$F_c = \frac{mv^2}{r} = m\omega^2 r$$向心力是"指向中心的牵拉力"——对于用绳子转动的球体，绳子的张力充当向心力；对于绕地球运转的卫星，重力充当向心力；对于在弯道上行驶的汽车，路面的摩擦力充当向心力。
弯道行驶和侧向G
自动车以速度 $v$ 做半径为 $r$ 的弯道行驶时，需要的向心力（摩擦力）为 $F = mv^2/r$。侧向加速度（侧向G）为 $a = v^2/r$。假设干燥路面摩擦系数 $\mu = 0.8$，则极限速度为 $v_{max} = \sqrt{\mu g r}$。半径50 m的弯道，极限速度约为72 km/h。超过这个速度会发生侧滑。对于轮胎–路面的相互作用，需要用Abaqus等进行FEM仿真（包括轮胎模型）来评估行驶稳定性。
3. 角速度·角加速度·线速度的关系
角速度 $\omega$ [rad/s] 是角度的时间变化率：
$$\omega = \frac{d\theta}{dt}$$角加速度 $\alpha$ [rad/s²] 是其时间变化率：
$$\alpha = \frac{d\omega}{dt}$$距离旋转轴 $r$ 处的点的线速度和切向加速度为：
$$v = r\omega, \quad a_t = r\alpha$$RPM和rad/s的换算
$$\omega \,[\text{rad/s}] = \frac{2\pi n}{60} \quad (n\text{ : 转速 [rpm]})$$| 设备 | 典型转速 [rpm] | ω [rad/s] |
|---|---|---|
| 洗衣机（脱水） | 1,200 | ~126 |
| 汽车发动机（常用） | 2,000〜4,000 | 210〜420 |
| 工业电动机 | 1,500〜3,600 | 157〜377 |
| 涡轮增压器 | 100,000〜200,000 | 10,500〜21,000 |
| 喷气发动机高压涡轮 | 12,000〜16,000 | 1,260〜1,680 |
| 牙科钻 | 300,000〜500,000 | 31,400〜52,400 |
4. 离心力（非惯性系的拟态力）和CAE的处理
在惯性系（绝对坐标系）中离心力不存在——没有向心力的话质点只会直进。但是在旋转坐标系（与涡轮机一起旋转的坐标系）中描述运动时，由于坐标系的非惯性性，出现两个拟态力作为补正项："离心力"和"科里奥利力"。
$$\vec{F}_{cf} = -m\vec{\omega}\times(\vec{\omega}\times\vec{r}) = m\omega^2 r\,\hat{r} \quad \text{（离心力：向外）}$$ $$\vec{F}_{Cor} = -2m\vec{\omega}\times\vec{v} \quad \text{（科里奥利力：垂直于速度）}$$*DLOAD, TYPE=CENTRIF 指定角速度和旋转轴，自动计算离心力。
🧑🎓 学生
科里奥利力也需要在旋转体FEM中考虑吗？
🎓 博士
静力分析中不需要，但在求"旋转振动分析"的固有振动数时就必须考虑——如果不考虑科里奥利力，就无法正确评估进行波和后向波的分化（科里奥利效应导致的振动数分裂）。在Abaqus的 *FREQUENCY 中加上 CORIOLIS=YES 参数很重要。
5. 涡轮叶片的离心应力
对于断面均匀、密度恒定的杆，从旋转轴距离 $r_h$〜$r_t$ 的范围，根部（$r = r_h$）处的离心应力为：
$$\sigma_{centrifugal} = \frac{1}{2}\rho\omega^2(r_t^2 - r_h^2)$$其中 $\rho$ 是材料密度 [kg/m³]。比如镍基高温合金（$\rho = 8,200\,\text{kg/m}^3$）的涡轮叶片，$\omega = 1,500\,\text{rad/s}$，$r_t = 0.25\,\text{m}$，$r_h = 0.10\,\text{m}$ 时：
$$\sigma = \frac{1}{2}\times 8200 \times 1500^2 \times (0.25^2 - 0.10^2) \approx 570\,\text{MPa}$$加上热应力和气体荷载，镍基高温合金在1,000℃时的屈服应力降低到约200〜600 MPa，因此材料选择和冷却设计显得极其重要。
锥形叶片的应力优化
实际的涡轮叶片是锥形的，先端处截面积较小。减少先端质量可以降低对根部的离心荷载。通过数值求解Frenet积分来求最优锥形比，再用FEM验证，是实务中的典型流程。
6. 旋转体的轴对称FEM
涡轮盘、飞轮、滚筒等轴对称形状的旋转体，可以用轴对称分析（Axisymmetric Analysis）将3D问题简化为2D（仅需要截面形状），计算成本可以降低到原来的百分之一以下。
轴对称单元中的离心力
向轴对称体输入离心荷载 $\rho\omega^2 r$（径向体积力）时，轴对称2D分析会计算出径向应力 $\sigma_r$、周向应力 $\sigma_\theta$（环向应力）、轴向应力 $\sigma_z$ 三个分量。环向应力 $\sigma_\theta$ 通常是最大的，往往支配着盘的屈服判定。
$$\sigma_\theta \text{（环向应力）} = \rho\omega^2\frac{3+\nu}{8}\left(r_o^2 + r_i^2 + \frac{r_i^2 r_o^2}{r^2} - \frac{1+3\nu}{3+\nu}r^2\right) \quad \text{（均匀盘，解析解）}$$需要3D模型的情况
- 附带叶片等非轴对称结构的情况
- 局部缺陷或损伤的评估
- 键槽、螺栓孔等周向不均匀性很重要的情况
7. 陀螺效应和进动运动
高速旋转的物体拥有角动量 $\vec{L} = I\vec{\omega}$。当施加外部力矩（力的矩）时，角动量向量的方向改变，产生进动运动（Precession）。
$$\vec{\tau} = \frac{d\vec{L}}{dt} = I\vec{\omega} \times \vec{\Omega}_{prec}$$陀螺效应的特点是与直觉相反，比如用手指推陀螺时，陀螺不会倒下，而是改变方向，这就是典型的陀螺效应。
汽车车轮和陀螺力矩
行驶中的汽车车轮高速旋转。当转向盘操作时，车轮轴的方向改变，陀螺力矩作用在悬架上。高转速、大直径车轮（跑车）的陀螺力矩不可忽视，对操纵特性有影响。在FEM多体动力学分析（MBD）中需要考虑陀螺效应来模拟转向。
8. 旋转机械的振动分析和Campbell线图
在旋转机械中，与旋转速度成比例的频率分量与结构的固有振动数相同时，会产生共振（Resonance）。系统地管理这个现象的方法就是Campbell线图（Campbell Diagram）。
Campbell线图的构成
横轴为旋转速度 [rpm或rad/s]，纵轴为振动频率 [Hz或rad/s]。
- 水平线：结构的固有振动数（不依赖于旋转速度，但高速时由于科里奥利效应会分裂）
- 右上倾直线（阶数线）：旋转速度 $n$ 的 $k$ 倍励振频率（$kn/60$ [Hz]）。$k$ 是整数（1, 2, 3...），称为"阶数（Engine Order）"
- 交点：固有振动数与励振频率相等的共振点。如果在运行域内就很危险
临界转速（Critical Speed）
旋转轴的横向振动固有振动数与转速相等时的转速称为临界转速。超过临界转速也能稳定运行的设计（超临界运行）也很常见——蒸汽涡轮轴和高速离心分离机常采用临界转速快速通过、在超临界域运行的方式。
9. 实践：涡轮增压器的破裂分析
涡轮增压器的压缩机轮需要在规定最大转速的120〜150%过速条件下验证破裂强度，这是安全标准（ISO 22728等）要求的试验和分析。
破裂的物理
压缩机轮的最大离心应力（环向应力）达到材料的抗拉强度 $\sigma_B$ 时，轮会急速崩坏。比如铝合金轮（A2024-T6, $\sigma_B \approx 480\,\text{MPa}$），离心应力超过 $\sigma_B$ 的转速就是"破裂转速"。
FEM分析步骤
- 静力分析：输入离心荷载+温度分布（从涡轮侧的热传导）
- 弹塑性分析：考虑大变形和温度依赖的材料特性，进行非线性FEM
- 破裂转速推估：求环向应力达到 $\sigma_B$ 的转速
- 破裂后的碎片运动分析：用显式FEM评估破裂碎片是否会穿透外壳
热应力的复合荷载
涡轮侧受热的轮（最高超过1,000℃）同时承受热膨胀和离心力的复合荷载。高温时屈服应力大幅降低，热膨胀导致尺寸变化，需要进行"热结构耦合分析"同时考虑这两个因素。