Jeffcott 转子基本方程
$$N_c = \frac{60}{2\pi}\sqrt{\frac{k}{m}}\text{ rpm}$$
不平衡响应振幅:
$$A = \frac{e \cdot r^2}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}$$
$r = \omega/\omega_c$(转速比)
输入转子质量、轴刚度、偏心量和阻尼比,计算临界转速Nc。不平衡响应曲线与进动轨迹动画帮助您直观理解共振行为,适用于透平机械设计。
$$N_c = \frac{60}{2\pi}\sqrt{\frac{k}{m}}\text{ rpm}$$
不平衡响应振幅:
$$A = \frac{e \cdot r^2}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}$$
$r = \omega/\omega_c$(转速比)
这个模拟器基于最经典的Jeffcott转子模型。它把复杂的转子简化成一个装在无质量弹性轴中央的集中质量点。这个模型虽然简单,却能抓住临界转速现象的核心。
$$N_c = \frac{60}{2\pi}\sqrt{\frac{k}{m}}\quad \text{[rpm]}$$$N_c$:临界转速(转/分)
$k$:轴的刚度(N/m),轴越硬,k值越大
$m$:转子的质量(kg)
知道了临界转速在哪,我们还需要知道经过它时振动有多大。这由“不平衡响应振幅”公式描述,它告诉你振幅如何随转速比和阻尼变化。
$$A = \frac{e \cdot r^2}{\sqrt{(1-r^2)^2 + (2\zeta r)^2}}$$$A$:转子的振动振幅(m)
$e$:转子的质量偏心距(m),制造和安装中无法避免的“瑕疵”
$r$:转速比,$r = n / N_c$(工作转速/临界转速)
$\zeta$:阻尼比,衡量系统消耗振动能量能力的无量纲数
航空发动机与燃气轮机:这些高速旋转机械的转子必须快速通过多个临界转速才能达到工作转速。设计时需精确计算各阶临界转速,并利用挤压油膜阻尼器等部件提供足够的阻尼,确保安全过临界。
大型离心压缩机与泵:在石油化工领域,这类设备通常工作在亚临界区(工作转速低于一阶临界转速)。工程师使用此工具进行初步估算,确保工作转速有足够的隔离裕度(如API标准要求的±20%),防止因基础沉降或结垢导致质量变化而意外落入共振区。
电机与发电机:对于柔性转子(工作转速高于一阶临界转速),启动和停机过程是危险期。通过计算,可以制定安全的升速/降速方案,避免在临界转速附近长时间停留,防止振动过大损坏轴承或绕组。
故障诊断与动平衡:现场机组振动异常时,工程师会测量其振动随转速变化的曲线(坎贝尔图)。若在特定转速出现峰值,很可能就是碰到了临界转速。这有助于判断是设计问题,还是因部件脱落、结焦导致转子质量分布改变,从而引发了共振。
首先,要纠正“临界转速是绝对不可逾越的危险速度”这一误解。虽然在该转速下进行稳态运行确实危险,但许多高速旋转机械是以“超临界运行”为前提设计的。关键在于通过合理的阻尼设计,在启停过程中快速穿越该转速区间。例如,某涡轮分子泵的一阶临界转速约在10,000 rpm附近,但其工作转速高达80,000 rpm。停机时虽会在10,000 rpm区间出现短暂振动增大,但借助阻尼设计可确保安全穿越。
其次,关于仿真中“等效质量”与“等效刚度”的设置。这是现场最令人困扰的问题之一。例如,考虑轴心处集中质量(m)时,并非直接代入转子的实际质量,还需叠加轴自身质量的一部分(通常为1/3至1/2)。刚度(k)也往往不仅取决于轴刚度,轴承与支撑结构的刚度常起主导作用。需特别注意轴承类型(滚珠轴承 vs 滑动轴承)及油膜刚度的显著影响。
最后,关于“偏心量(e)归零即可消除振动”的理想化观点。现实中,加工装配误差、材料不均匀性、运行中的热变形及磨损必然导致不平衡。通过工具调整e值观察振幅变化,旨在预测动平衡校正的效果。例如定性理解:将偏心量从10μm减小至5μm(进行平衡修正),共振点振幅可降低约一半。相比追求完全归零,将振幅控制在允许范围内的成本意识更为重要。
临界转速的计算思路,本质上与诸多工程领域的“共振问题”异曲同工。首当其冲的是结构动力学。建筑与桥梁在风载或地震作用下的振动分析,同样基于固有频率(对应临界转速)与强迫振动力(对应不平衡力)的关系。模态分析方法也用于求解转子复杂的振动模态。
另一领域是控制工程,尤其是“振动控制”方向。仿真中增大阻尼比(ζ)可抑制共振峰值的现象,类似于主动控制(主动振动控制)中“提高反馈增益以增强稳定性”的操作。实际上,采用磁轴承的转子系统,其控制系统设计直接决定阻尼特性,并影响临界转速穿越的难易程度。
此外,该问题与材料力学·疲劳分析也密切相关。共振时的大幅振动会在轴内引发交变应力,大幅缩短材料疲劳寿命。例如,泵轴在临界转速附近长期运行后早于设计寿命断裂的案例并不罕见。振动分析获得的应力幅值与循环次数,可直接作为疲劳寿命预测的输入数据。
熟悉Jeffcott转子模型后,下一步可拓展至多自由度系统。实际转子更接近带多个叶盘的“连续体”,可通过多质量-弹簧模型建模,并用矩阵建立运动方程($[M]\{\ddot{x}\} + [C]\{\dot{x}\} + [K]\{x\} = \{F\}$)。求解该方程可获得第一阶、第二阶、第三阶…等多阶临界转速及对应振型(挠曲形态),这是绘制“坎贝尔图”的基础。
数学层面,建议深入理解频率响应函数的推导过程。这实质是二阶线性常微分方程的强迫振动解。建议从无阻尼解入手,逐步加入阻尼项,再引入不平衡强迫力项…通过教材梳理此流程。理解这一原理后,便能透彻理解工具输出曲线的形态(为何r=1处出现峰值,以及阻尼如何使峰值圆钝化)。
与工程实践紧密相关的下一主题是“不稳定振动”机理。临界转速类共振不称为“不稳定”。旋转机械特有的“油膜涡动”与“高速失稳”等现象,会在特定转速区间引发发散性振动。这些属于阻尼将振动能量从“吸收”转为“供给”的自激现象,是超临界运行设计的核心难点之一。建议从判别稳定性的“瑞利稳定性判据”概念入手学习。