临界转速的计算

一、临界转速分析的目的

临界转速分析的主要目的在于确定转子支撑系统的临界转速，并按照经验或有关的技术规定，将这些临界转速调整，使其适当的远离机械的工作转速，以得到可靠的设计。 例如设计地面旋转机械时，如果工作转速低于其一阶临界转速Nc1，应使N<0.75Nc1, 如果工作转速高于一阶临界转速，应使1.4Nck要较为准确的确定出转子支撑系统的临界转速，必须注意以下两点

1. 所选择的计算方法的数学模型和边界条件要尽可能的符合系统的实际情况。

2. 原始数据的（系统支撑的刚度系数和阻尼系数）准确度，也是影响计算结果准确度

的重要因素。

3. 适当的考虑计算速度，随着转子支撑系统的日益复杂，临界转速的计算工作量越来

越大，因此选择计算方法的效率也是需要考虑的重要因素。

三、常用的计算方法 名称 1. 矩阵迭代法 （Stodola 斯托多拉） 2. 3. 4. 5. 原理 假定一阶振型挠曲弹性线并选择试算速度 计算转子涡动惯性载荷，并用此载荷计算挠性曲线 以计算得到的挠性曲线和适当调整的转速重新循环计算 当计算曲线和初始曲线吻合的时的转速即为一临转速 高阶临界转速方法同，但需利用正交条件消除低阶弹性线成分，否则计算错误 将四个状态参1. 逐段推算法 （传递矩阵法） （Prohl-Myklestad） 3. 4. 根据与其相邻轴段在该截面处的约束条件，得到下个轴段的状态参数 换个转速重复计算，直到计算得到的状态参数满足边界条件，此时的转速即为临界转速 机械制造业中是最为通用、发展最为完善的方法 如果假设的振型能量法 （Rayleigh-Ritz） 1. 以能量守恒原理为理论基础，根据轴系中的最大应变能等于最大的动能，建立微分方程，据动能是转速的函数计算转速 原理简单，易于不准确会带来误理解 差 方程难以求解，特征方程法 将通用的指数解带入微分方程，得到以临界转速为解的多项式 应用不多 唯一能模拟非线性系统的计算方法，在校核数值积分法 （前进法） 以数值积分的方法求解支撑系统的运动微分方程，从初始条件开始，以步长很小其他方法及研的时间增量时域积分，逐步推算出轴系的运动 究非线性对临界转速的影响方面很有价值 步数 须有足够的积分计算量较大，必计算转速，比较盲目 2. 划分转轴为若干等截面段，选择试算转速 从转轴的一端算起，计算另一端的四个状态参数（挠度、转角、弯矩、剪力） 数写成矩阵的形式，计算方便，在各类旋转根据经验或有关的计术资料选择收敛较快，一阶高阶临界转速精临界转速结果度差，计算复杂 较为准确 优点 缺点 注：1.Stodola 斯托多拉法 2.Prohl-Myklestad莫克来斯塔德法

不得用于商业用途

仅供个人参考

传递矩阵法

基本原理：传递矩阵法的基本原理是，去不同的转速值，从转子支撑系统的一端开始，循环进行各轴段截面状态参数的逐段推算，直到满足另一端的边界条件。

优点：对于多支撑多元盘的转子系统，通过其特征值问题或通过建立运动微分方程的方法求解系统的临界转速和不平衡响应，矩阵的维数随着系统的自由度的增加而增加，计算量往往较大：采用传递矩阵法的优点是矩阵的维数不随系统的自由度的增加而增大，且各阶临界转速计算方法相同，便于程序实现，所需存储单元少，这就使得传递矩阵法成为解决转子动力学问题的一个快速而有效的方法。

缺点：求解高速大型转子的动力学问题时，有可能出现数值不稳定现象。今年来提出的Riccati传递矩阵法，保留传递矩阵的所有优点，而且在数值上比较稳定，计算精度高，是一种比较理想的方法，但目前还没有普遍推广。

轴段划分：首先根据支撑系统中刚性支撑（轴承）的个数划分跨度。在整个轴段内，凡是轴承、集中质量、轮盘、联轴器等所在位置，以及截面尺寸、材料有变化的地方都要划分为轴段截面。若存在变截面轴，应简化为等截面轴段，这是因为除了个别具有特殊规律的变截面轴段外，其他的变截面轴段的传递矩阵特别复杂。

传递矩阵：

4. 轴段传递矩阵

每段起始状态参数和终端状态参数的转换方程，根据是否考虑转轴的分布质量，可以建立两种轴段传递矩阵

① 当考虑轴段的分布质量时：起始和终端的转换方程是均质等截面杆的振动弹性方

程：

② 不考虑转轴的分布质量时建立的传递矩阵

X1L120122M001Qki00012L-11X22L-21L1

MQ0i其中，a11,a12,a21,a22为该轴段的影响系数，根据材料力学：

L3113EJL2，a11和a12是终端的剪力和弯矩在终端引起的挠度，a21和12212EJL22EJ不得用于商业用途

仅供个人参考

a22是终端的剪力和弯矩在终端引起的转角

4. 各轴段间的传递矩阵

从前一轴段的终端到下一轴段的始端，如果中间没有的结构单元，则状态参数不发生变化，传递矩阵是单位矩阵；两者之间有的结构单元时，用前一轴段的终端矩阵乘以此单元的矩阵，即的下一单元的始端矩阵。的结构单元大概可以分为以下四种：

a. 通过点质量时为：

1002mpi000100，其中，mi为点质量，p为系统的固有频率 010001b. 通过转动盘时为

1002miP0wIpIdP21*pId010000，其中，mi盘的质量，Ip盘的极转动1001惯量，Id盘的直径转动惯量，w盘的转动角速度 c. 通过弹性铰链时为

1000010001ch1000，其中，ch为铰链的力矩刚性系数 01000，其中，co弹性支座的刚性系数，如果没有弹性约束则ch=0. 01d. 通过具有弹性约束的弹性支座时为

100c001ch01004. 各跨度间的传递矩阵

a. 通过刚性支座的传递

刚性支座是一个跨度的结束，在支座处的横向位移为0，所以：

X0（i1）Xki00(i1)ki其中，Ri为支座的反作用力，在以后整个跨度的计算中，MM0(k1)kiQ0(i1)QkiRi此反作用力代替前一跨度中被消除的参数（挠度），而未知参数的个数不变。 b. 通过球头联轴器的传递

球头联轴器也是一个跨度的结束，在此处的弯矩为0.所以：

不得用于商业用途

仅供个人参考

X0（i1）Xki0(i1)kii其中，球头联轴器未知的相对转角，在以后此跨度的计算MM0ki0(k1)Q0(i1)Qki中，用θ代替上个跨度中消除的参数，从而使未知变量的个数不变。

4. 初始条件：第一跨度0截面的初始条件根据约束条件和轴的载荷分析来确定，在所有四

个状态参数中，或有两个为零，两个是未知的，或只有两个是的，其他的参数可以用这两个参数表示。这就意味着在计算过程中所有各段的起始端和末端的状态参数都是两个未知数的线性函数。最主要的是末端的状态参数也总是或两个为0，或可以用两个参数来表示，因此末端的四个参数方程可以简化为两个具有两个未知数的齐次方程。

5. 临界转速的确定：转子临界转速的确定可以用“瞎子爬山”、对分法等来确定。选取某个P值，写出所有轴段的传递矩阵，然后根据初始端的边界条件选取合适的初始参数矩阵。从转子的起始端逐段推算其状态参数，在每个跨度的终端，按照条件进行参数的消除和变换，最终递推到末端时，可以得到两个含有两个未知数的齐次方程。假设齐次方程的系数行列式为0，着计算转速就是临界转速；若行列式不为零，则重新选取临界转速计算。将各阶临界转速带入重新计算可得各段始、末端的参数，从而作出振型图。计算过程中，可以将第一跨度的初始截面的某个状态参数设为1，以后各截面的参数值是相对于1的比例值。

临界转速计算：

不得用于商业用途

仅供个人参考

不得用于商业用途

仅供个人参考

单圆盘转子的临界转速和不平衡响应

早期的旋转机械比较简单，可以把转子看做是圆盘装在无重的弹性转轴上，而转轴的两端则由完全刚性即不变形的轴承及轴承座支撑，这种模型成为刚性支撑。 1.1 涡动的定义

通常转轴的两支点在同一水平线上，转轴未变形时，转子的轴线处于水平位置，（实际上由于盘的重力作用，即使在静止时，转轴也会变形，而不是处于水平位置），由于转子的静变形交小，对转子的运动的影响可以忽略不计。有时为了避开静变形，可以考虑让转轴的两支点在同一垂直线上。

假设转子以角速度Ω做等速转动，当处于正常运转时，轴线是直的，如果在他的一侧添加一横向冲击，则因转轴有弹性而发生弯曲振动，涡动就是研究这种性质的运动。假设圆盘的质量为m，他所受到的力是转轴的弹性恢复力F= -kr，其中k为转轴的刚度系数，R=oo’，圆盘的运动微分方程：

由式1.4可知，圆盘或转轴的中心o’在相互垂直的两个方向作频率同为Wn的简谐振动。在一般情况下，振幅X和Y是不相同的，式1.4确定点轨迹为一椭圆，o’的这种运动成为“涡动”，自然频率Wn称为进动角速度。

不得用于商业用途

仅供个人参考

其中B1和B2都是复数，由起始的横向冲击决定。第一项是半径为B1的反时针运动，运动方向和转动角速度相同，成为正进动。第二项是半径为B2的顺时针运动，运动方向和转轴的转动方向相反，成为反进动。圆盘中心o’的涡动就是这两种进动的合成。由于起始条件的不同，转子中心的涡动可能出现以下情况：

①. B1!=0，B2=0 涡动为正进动，轨迹为圆，半径为B1； ②. B1=0，B2!=0 涡动为正进动，轨迹为圆，半径为B2； ③. B1=B2轨迹为直线；

④. B1!= B2轨迹为椭圆， B1>B2时为正涡动；B1由以上讨论可知，圆盘或转轴的中心的进动或涡动属于自然振动，他的频率就是圆盘没有转动时，转轴弯曲振动的自然频率。

1.2 圆盘的偏心质量引起的振动、临界转速

Ywεwt观察方向Xm

假设转盘的质量为m，偏心距为ε，角速度为w，设离心力的初始相位α为0，则在某一时刻t，离心力矢量和x轴的夹角为wt，此时离心力在X和Y向的投影为：

Fxm2costFymsint2

Fx和Fy分别是各自方向上的周期性变化的力，频率和转盘的频率相同，在这种交变力的作用下，转子在X和Y方向也将做周期性运动， 假设两个方向上阻尼和刚度相同，则转子的运动微分方程：

mx''cx'kxFconwt

my''cy'kyFsinwt不得用于商业用途

仅供个人参考

其解为：

Zx(t)Zcos(wt)y(t)Zsin(wt)Fw1-wn2arctan2w/wn2寇胜利 2w/wn221w/wn钟一谔

结论：

1. 只考虑强迫振动时，轴心的响应频率和偏心质量的激振频率相同，在转速小于临界转速

时且不考虑阻尼时，相位也相同，轴心和质心在一条直线上；当转速大于临界转速时且不考虑阻尼时，相位相差180°。

2. 当考虑转子的涡动时，运动比较复杂；

3. 不平衡矢量所在的位置成为重点，振动矢量所在的位置成为高点，高点比重点滞后的角

度成为滞后角，当令阻尼比为0时，φ为0，说明滞后角是由阻尼引起的；

4. 转子存在偏心，运行的过程中又出现动挠度，当转速小于临界转速时，挠度和F即偏心

方向相同，使终偏心增大；当转速等于临界转速时，出现共振；当转速大于临界转速时，挠度方向和偏心方向相反，使终偏心减小，转子振动趋于平稳，这种现象成为自动对心；

1.3等截面转子的振动

并不是所有的转子系统都可以简化为具有刚性支撑的单轮盘转子系统模型，对于均质、等截面转子，如果按照集中质量处理，将不能反映真实振动特性。均质、等截面转子系统的运动规律可以用一个偏微分方程表示，该偏微分方程含有时间和轴向位置两个自变量，因此可以确定任意轴线位置在任意时刻的位置，利用均质、等截面转子模型研究得出的结论对一般转子也是适用的。

不得用于商业用途

仅供个人参考

Y 轴Z轴

运动方程：如图上图所示的两端简支的等截面转子，设其密度为P，截面面积为A，弯曲刚度为EI，分布干扰力在xoz和yoz平面分别为Fx(z,t) Fy(z,t)，则转子的振动可以用以下一组微分方程组成：

2xz,t4x(z,t)AEIfz,tx令分布干扰力为0，t2z4即可得到转子的自由2yz,t4y(z,t)AEIfz,ty振动微分方程： t2z4X 轴2xz,t4x(z,t)AEI024tz

其解为：

xz,tDnsinn1nzconwntnlEIA

n22其中，固有频率wn2lsinnz为振型函数，Dn和n分别为n阶自由振动的振型和初相位l由上式可知转子的自由振动是一系列简谐振动的合成，这些简谐振动有以下特点：

①. 固有频率和振型函数是一一对应的； 不得用于商业用途

仅供个人参考

②. 振型函数sinnz反映了转子轴线上各点位移的相对比例关系，无论振幅Dn如何l变化 ，这种比例关系不会变化；

③. 振型是由转子-支撑系统自身的特点决定的，所以又称为固有振型，不同类型的转子

系统的振型函数不同，上述的是均质等截面转子的振型函数。 有关振型的基本概念：

a) 节点：轴线上某一点的振型函数值称为该点的振型值，振型值为0的点成为节

点，阶数越高节点越多，N阶振型的节点数为N-1；

b) 对称性：对于两端简支的等截面转子，奇数阶振型是对称的，而偶数阶振型是

反对称的。因此。在两支座间，奇数阶振型相位相同而偶数阶振型相位相反； c) 正交性：转子的不同阶振型间具有正交性，即第m阶振型和第n阶振型的乘

积在轴长上的积分为0。

d) 理论上，转子的1、2、3阶固有频率的比值是1：4：9，实际1、2阶固有频

率间的比值为1：3左右；

e) 理论上，转子-支撑系统经过临界点时，相位变化180°，实际上由于阻尼的

存在，在临界转速处相位一般变化90°，即振动矢量和不平衡矢量间的滞后角为90°。

f) 如下图所示，由于阻尼的存在，转子中心对不平衡质量的响应在w=Wn处不

仅不是无限大值，而且不是最大值，最大值发生在w略微大于Wn时。对于实际的转子系统，有时通过在升速或降速的过程中测量响应的办法来确定转子的临界转速，常常把这个过程中的最大值即峰值的转速作为临界转速。有图可知，通过测量所获得的临界转速在升速上略大于实际的临界转速，而在降速时这略小于实际的临界转速。

1.3 陀螺力矩

基本概念：

1. 对质点的动量距：质点Q的动量对于点O的距定义为质点对于点O的动

量距，其值为点O到质点Q的矢量差乘以动量：Mo(mv)=Rxmv，方向按照右手定则判定。

2. 对轴的动量距：质点Q的动量在xoy面内的投影mv(xy)对与O点的距定

义为质点Q对Z轴的动量距；

3. 刚体对轴的转动惯量：刚体的转动惯量是刚体转动时惯性的度量，它等于

刚体内各质点的质量与该质点到轴的垂直距离平方的乘机的和。

4. 赖柴定理：质点系对固定点的动量距矢量断点的速度等于外力系对同一点

的主距。 不得用于商业用途

仅供个人参考

Mf2UF1 F2f1

当圆盘不在两支撑的中点而偏于一边时，转轴变形后，圆盘的轴线和两指点的连线AB有一夹角φ。设圆盘的自转角速度为ｗ，极转动惯量为Jp，则圆盘对质心的动量距为：H=Jpw，根据右手定则，它与AB连线的夹角也为φ。设转轴涡动的频率为Wn，则圆盘中心o’与轴线AB所构成的平面绕AB轴有进动角速度。由于进动，圆盘的动量距H将不断变化，因此动量距矢量的终点将具有速度U，根据赖柴定理（质点系对固定点的动量距矢量断点的速度等于外力系对同一点的主距），而圆盘重力距等于0，显然和动量距矢量终点的速度相等的外力距只可能是轴承的动反力F1、F2产生的力矩；

力矩-M（根据作用和反作用）称为陀螺力矩，它是圆盘施加与转轴的力矩，相当于弹性力矩。在正进动（0<φ<π/2）时，它是转轴的变形减小，从而提高了转轴的弹性刚度，即提高了转子的临界转速；在反进动（π/2<φ<π）时，它是转轴的变形增大，从而降低了转轴的弹性刚度，即降低了转子的临界转速。

当机械中的高速转动部件的对称轴被迫在空间中改变方位时，即对称轴被迫进动时，转动部件必须对约束作用一个附加力偶，这种现象称为陀螺效应。当陀螺效应严重时，可能使机械产生故障，尤其是轴承。

1.4 弹性支撑对转子临界转速的影响

Jeffcott转子：这种转子模型是对真实转子的简化，刚性支承的单盘转子，单盘位于支承的中间，分析临界转速和陀螺力矩等，是转子动力学的基础。 假设盘在平面内运动，不考虑轮盘的偏转，轴是无重轴。

不得用于商业用途

仅供个人参考

临界转速计算：200000000000 0.000000000982

1. 基本参数：截面惯性矩：J=4.91E-10，弹性模量E=2E11，右端质量m3=0.096g=0.1kg，

两个盘的质量m1=m2=0.8kg，Ipm1r1r222/2=5.73E-4，Id=Ip/2=2.86E-4。

2. 各轴段的传递矩阵：

第一段：l=0.045m，J=9.82E-10，a11=1.55E-7，a12=a21=5.16E-6，a22=2.29E-4，

10.0455.16x1064012.29x100010007.73x1085.16x106 0.0451第二段：l=0.11m，J=9.82E-10，a11=2.26E-6，a12=a21=3.08E-5，a22=5.6E-4，

10.113.08x1051.128x10645015.6x103.08x10 0010.110001第三段：l=0.15m，J=9.82E-10，a11=5.73E-6，a12=a21=5.73E-5，a22=7.E-4，

10.155.73x1054017.x100010002.865x1065.73x105

0.151第四段：l=0.11m，J=9.82E-10，a11=2.26E-6，a12=a21=3.08E-5，a22=5.6E-4，

10.113.08x1051.128x10645015.6x103.08x10 0010.110001第五段：l=0.01，J=9.82E-10，a11=1.70E-9，a12=a21=2.55E-7，

10.012.55x107015.09x105a22=5.09E-50010008.5x10102.54x107

0.011X0101初始参数列阵为：，令X01=1，选取P=2050r/s，用第一段的矩阵乘此矩阵，即

00不得用于商业用途

仅供个人参考

010Xk110.04522.22k101 可得此段的终端参数：Mk10Qk10X02022.2202，用第二段的传递矩阵乘此列阵，得终端参 第二段的始端参数列阵为：

M020Q-RA02Xk22.44421.128x106RA522.223.08x10Rk2A数：， Mk20.11RAQRk2A通过转盘的传递矩阵：

100101.2x103603.362x10001000，用此矩阵乘以第二段的终端矩阵得第三段的起始参数： 01X032.44421.128x106RA522.223.08x10RA03M032.67x1040.147RA 6Q8.22x104.79R03AXk3-30.854-2.7x105RA4513.624.18x10Rk3A 第三段的终端参数：6Mk31.26x100.8655RA6Qk38.22x104.79RA通过转盘的传递矩阵：

100101.2x103603.362x10001000，用此矩阵乘以第三段的终端矩阵得第四段的起始参数： 01X04-30.854-2.7x10-5RA4513.624.18x10RA04 6M041.876x101.367RA8Q1.12x1098.556RA04

不得用于商业用途

仅供个人参考

-4Xk4-219.47-2.27x10RA-35013.784.219x10Rk4A第四段的终端参数：，该段末端为刚性支撑，估位移为Mk41.4196x1077.72RA8Qk41.12x1012.21RA零，从而可得RA=-9.668E5，第五段的起始矩阵为：

X05020005

6M056.44x107Q8..27x10RB05Xk53.718.5x1010RB-7-548.81-2.54x10Rk5B第五段的终端矩阵：，点质量的传递矩阵为： 6Mk57.367x100.01RB7Q8..27x10Rk5B10054.2x1001000010X3.718.5x1010RB0-7-548.81-2.54x10R0B，转子终端的传递矩阵为：，最后一6M7.367x100.01RB07Q8.426x10R1B段是一个悬臂，其右端没有外力和力矩，因此上式的后两式为零。系数行列式展开：Λ

=6.524E5.

不得用于商业用途

仅供个人参考

不得用于商业用途

仅供个人参考

仅供个人用于学习、研究；不得用于商业用途。

For personal use only in study and research; not for commercial use.

Nur für den persönlichen für Studien, Forschung, zu kommerziellen Zwecken verwendet werden. Pour l 'étude et la recherche uniquement à des fins personnelles; pas à des fins commerciales.

только для людей, которые используются для обучения, исследований и не должны использоваться в коммерческих целях.

以下无正文

不得用于商业用途