临界转速的计算Word格式文档下载.docx
- 文档编号:4524519
- 上传时间:2023-05-03
- 格式:DOCX
- 页数:20
- 大小:229.98KB
临界转速的计算Word格式文档下载.docx
《临界转速的计算Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《临界转速的计算Word格式文档下载.docx(20页珍藏版)》请在冰点文库上搜索。
建立微分方程,据动能是转速的函数计算转速
原理简单,易于
理解
如果假设的振型
不准确会带来误
差
特征方程法
将通用的指数解带入微分方程,得到以临界转速为解的多项式
方程难以求解,
应用不多
数值积分法
(前进法)
以数值积分的方法求解支撑系统的运动微分方程,从初始条件开始,以步长很小
的时间增量时域积分,逐步推算出轴系的运动
唯一能模拟非线性系统的计算方法,在校核其他方法及研究非线性对临界转速的影响方面很有价值
计算量较大,必须有足够的积分
步数
注:
I.Stodola斯托多拉法
2.Prohl-Myklestad莫克来斯塔德法
传递矩阵法
基本原理:
传递矩阵法的基本原理是,去不同的转速值,从转子支撑系统的一端开始,循环
进行各轴段截面状态参数的逐段推算,直到满足另一端的边界条件。
披第四堂的方怯卜把转子简化为集息质量復堂后,传递矩阵注是把爲统分为圜盘*轴段和支般等若干个典型的单无成部幷*用力学方怯建立这些部件两規截面伏态向战间的传递关系,禅利用连壤条件就可求得转子在任-截面的状卷向疑与雄始栽面的狀态向趙阿的关索,通过对能满足岂界条件的涡动頻曝的搜索'
就疔得出转子系统的各阶临界转速。
1+算在琳定自转角建度时转子在不平撕质最激席F的握动,就可求得转子的不平衡陥]应。
为此先讨港备料典型部件的传递矩阵“这歧典型部件可以呈画盘、咕段、弹性支承或其姐合等◎'
优点:
对于多支撑多元盘的转子系统,通过其特征值问题或通过建立运动微分方程的方法求
解系统的临界转速和不平衡响应,矩阵的维数随着系统的自由度的增加而增加,计算量往往
较大:
采用传递矩阵法的优点是矩阵的维数不随系统的自由度的增加而增大,且各阶临界转
速计算方法相同,便于程序实现,所需存储单元少,这就使得传递矩阵法成为解决转子动力学问题的一个快速而有效的方法。
缺点:
求解高速大型转子的动力学问题时,传递矩阵法,保留传递矩阵的所有优点,理想的方法,但目前还没有普遍推广。
轴段划分:
首先根据支撑系统中刚性支撑(轴承)的个数划分跨度。
在整个轴段内,凡是轴承、集中质量、轮盘、联轴器等所在位置,以及截面尺寸、材料有变化的地方都要划分为轴段截面。
若存在变截面轴,应简化为等截面轴段,这是因为除了个别具有特殊规律的变截面轴段外,其他的变截面轴段的传递矩阵特别复杂。
传递矩阵:
4.轴段传递矩阵
每段起始状态参数和终端状态参数的转换方程,根据是否考虑转轴的分布质量,可以建
立两种轴段传递矩阵
1当考虑轴段的分布质量时:
起始和终端的转换方程是均质等截面杆的振动弹性方
程:
2不考虑转轴的分布质量时建立的传递矩阵
X、
■1
L
«
12
a12L^111
Q
1
a22
a22L^21
e
M
2j
ki
卫
1一
IQ”
其中,a11,a12,a21,a22为该轴段的影响系数,根据材料力学:
L3
3EJ
L2
12=「21=•
2EJ
ail和
a12是终端的剪力和弯矩在终端引起的挠度,
a21和
EJ
a22是终端的剪力和弯矩在终端引起的转角
4.各轴段间的传递矩阵
从前一轴段的终端到下一轴段的始端,如果中间没有独立的结构单元,则状态参数不发生变化,传递矩阵是单位矩阵;
两者之间有独立的结构单元时,用前一轴段的终端矩阵乘以此单元的矩阵,即的下一单元的始端矩阵。
独立的结构单元大概可以分为以下四种:
b.
mip
通过转动盘时为
1000
01丄0
6,其中,ch为铰链的力矩刚性系数
0010
]0001一
Xo(i书=Xki
=0(i1)=ki=i
0(i+)kii其中,球头联轴器未知的相对转角,在以后此跨度的计算
M0你书=M总=0
Qo(iqi)=Qki
中,用B代替上个跨度中消除的参数,从而使未知变量的个数不变。
4.初始条件:
第一跨度0截面的初始条件根据约束条件和轴的载荷分析来确定,在所有四
个状态参数中,或有两个为零,两个是未知的,或只有两个是独立的,其他的参数可以用这两个独立参数表示。
这就意味着在计算过程中所有各段的起始端和末端的状态参数都是两个未知数的线性函数。
最主要的是末端的状态参数也总是或两个为0,或可以用
两个参数来表示,因此末端的四个参数方程可以简化为两个具有两个未知数的齐次方
程。
5.临界转速的确定:
转子临界转速的确定可以用"
瞎子爬山”、对分法等来确定。
选取某
个P值,写出所有轴段的传递矩阵,然后根据初始端的边界条件选取合适的初始参数矩阵。
从转子的起始端逐段推算其状态参数,在每个跨度的终端,按照条件进行参数的消除和变换,
最终递推到末端时,可以得到两个含有两个未知数的齐次方程。
假设齐次方程的系数行列式
为0,着计算转速就是临界转速;
若行列式不为零,则重新选取临界转速计算。
将各阶临界转速带入重新计算可得各段始、末端的参数,从而作出振型图。
计算过程中,可以将第一跨
度的初始截面的某个状态参数设为1,以后各截面的参数值是相对于1的比例值。
(上摄第和则方袪汁算转子临界转速时.转于念出现正进动和反进动q由十陀螺效应的柞用,Mi疳转子门转角速贋I旳捉高,氏进动固青频率降低,而正进动M]有频率将提帕抿据临界转速询定文,阪只対正进动IS冇顼率进行分忻在后魅理申官先剔除血固有饮率*惦后分析各諦模态模犁论疋冋一阶抵熨的
幼固有频執改变转f自转角邃厦血计和出新的"
・嚴启曲出曲线Z即转子坎贝尔曲线'
如禺示月根据临界转速辽SL当曲求临界睥速缶嗖
临界转速计算:
△L・L-r!
单圆盘转子的临界转速和不平衡响应
早期的旋转机械比较简单,可以把转子看做是圆盘装在无重的弹性转轴上,而转轴的两
端则由完全刚性即不变形的轴承及轴承座支撑,这种模型成为刚性支撑。
1.1涡动的定义
通常转轴的两支点在同一水平线上,转轴未变形时,转子的轴线处于水平位置,(实际
上由于盘的重力作用,即使在静止时,转轴也会变形,而不是处于水平位置),由于转
子的静变形交小,对转子的运动的影响可以忽略不计。
有时为了避开静变形,可以考虑让转轴的两支点在同一垂直线上。
图J.i图ia
假设转子以角速度Q做等速转动,当处于正常运转时,轴线是直的,如果在他的一
侧添加一横向冲击,则因转轴有弹性而发生弯曲振动,涡动就是研究这种性质的运动。
假设圆盘的质量为m,他所受到的力是转轴的弹性恢复力F=-kr,其中k为转轴的刚度
系数,R=oo'
圆盘的运动微分方程:
帼£
二Fz=亠
r
m¥
=F*二-
F
y
令・
mJ-k/m
a,2)
则有
世*硏夷=i)1§
*心卜-Jj
(1.3)
它的解呵写作
x=Xcos<
djay-Ys\n(0).
2叭)〔
(L4)
由式1.4可知,圆盘或转轴的中心o'
在相互垂直的两个方向作频率同为Wn的简谐振动。
在一般情况下,振幅X和Y是不相同的,式1.4确定点轨迹为一椭圆,o'
的这种运动成为“涡动”,自然频率Wn称为进动角速度。
U(K3)式变为
其解为
N=琴+炒(|=y/
£
+wjs=0
ZMB岸诃£
谊7扌,
6
(1,5>
(1.6>
其中B1和B2都是复数,
由起始的横向冲击决定。
第
「项是半径为
B1的反时针运动,
运动方向和转动角速度相同,成为正进动。
第二项是半径为B2的顺时针运动,运动方向和
转轴的转动方向相反,成为反进动。
圆盘中心o'
的涡动就是这两种进动的合成。
由于起始
条件的不同,转子中心的涡动可能出现以下情况:
1.B1!
=0,B2=0涡动为正进动,轨迹为圆,半径为B1;
2.B仁0,B2!
=0涡动为正进动,轨迹为圆,半径为B2;
3.B仁B2轨迹为直线;
4.B1!
=B2轨迹为椭圆,B1>
B2时为正涡动;
B1<
B2时为反涡动。
由以上讨论可知,圆盘或转轴的中心的进动或涡动属于自然振动,他的频率就是圆盘没
有转动时,转轴弯曲振动的自然频率。
1.2圆盘的偏心质量引起的振动、临界转速
假设转盘的质量为m,偏心距为£
,角速度为w,设离心力的初始相位a为0,则在某一时刻t,离心力矢量和x轴的夹角为wt,此时离心力在X和丫向的投影为:
2
Fx=m:
.cost
Fy二m;
2sint
Fx和Fy分别是各自方向上的周期性变化的力,频率和转盘的频率相同,在这种交变力的作用下,转子在X和丫方向也将做周期性运动,假设两个方向上阻尼
和刚度相同,则转子的运动微分方程:
mx'
'
ex'
kx二Fconwt
my'
cy'
ky二Fsinwt
其解为:
x(t)二Zcos(wt-)y(t)=Zsin(wt-)ZF
代人(1.13〉式后,可玻得振幅
圆盘或转轴中心对于不平的响应为
(AUS)
钟一谔
结论:
1.只考虑强迫振动时,轴心的响应频率和偏心质量的激振频率相同,在转速小于临界转速时且不考虑阻尼时,相位也相同,轴心和质心在一条直线上;
当转速大于临界转速时且不考虑阻尼时,相位相差180°
。
2.当考虑转子的涡动时,运动比较复杂;
3.不平衡矢量所在的位置成为重点,振动矢量所在的位置成为高点,高点比重点滞后的角
度成为滞后角,当令阻尼比为0时,0为0,说明滞后角是由阻尼引起的;
4.转子存在偏心,运行的过程中又出现动挠度,当转速小于临界转速时,挠度和F即偏心
方向相同,使终偏心增大;
当转速等于临界转速时,出现共振;
当转速大于临界转速时,挠度方向和偏心方向相反,使终偏心减小,转子振动趋于平稳,这种现象成为自动对心;
1.3等截面转子的振动
并不是所有的转子系统都可以简化为具有刚性支撑的单轮盘转子系统模型,对于均质、等截面转子,如果按照集中质量处理,将不能反映真实振动特性。
均质、等截面转子系统的运动规律可以用一个偏微分方程表示,该偏微分方程含有时间和轴向位置两个自变量,因此
可以确定任意轴线位置在任意时刻的位置,利用均质、等截面转子模型研究得出的结论对一
般转子也是适用的。
运动方程:
如图上图所示的两端简支的等截面转子,设其密度为P,截面面积为A,弯
曲刚度为EI,分布干扰力在xoz和yoz平面分别为Fx(z,t)Fy(z,t),则转子的振动可以用以下一组微分方程组成:
:
A"
xz,tEI几(严=fxz,t
t:
z
Q2:
EI5严7即
.t:
令分布干扰力为
o,即可得到转子的自由振动微分方程:
nz
xz,t-'
Dnsinconwnt-;
nn=1l
其中,固有频率w
2_.2
n■-
n:
】z
Dn和-:
n分别为n阶自由振动的振型和初相位
sin节为振型函数,
由上式可知转子的自由振动是一系列简谐振动的合成,这些简谐振动有以下特点:
1.固有频率和振型函数是一一对应的;
2•振型函数sin反映了转子轴线上各点位移的相对比例关系,无论振幅Dn如何
l
变化,这种比例关系不会变化;
3.振型是由转子-支撑系统自身的特点决定的,所以又称为固有振型,不同类型的转子
系统的振型函数不同,上述的是均质等截面转子的振型函数。
有关振型的基本概念:
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.赖柴定理:
质点系对固定点的动量距矢量断点的速度等于外力系对同一点的主距。
当圆盘不在两支撑的中点而偏于一边时,转轴变形后,圆盘的轴线和两指点的连线AB
有一夹角0。
设圆盘的自转角速度为w,极转动惯量为Jp,则圆盘对质心的动量距为:
H=Jpw,根据右手定则,它与AB连线的夹角也为0。
设转轴涡动的频率为Wn,则圆盘中心o'
与轴线AB所构成的平面绕AB轴有进动角速度。
由于进动,圆盘的动量距H将不断变化,因此动量距矢量的终点将具有速度U,根据赖柴定理(质点系对固定点的动量距矢量断点的速度等于外力系对同一点的主距),而圆盘重力距等于0,显然和动量距矢量终点的速度相等的外力距只可能是轴承的动反力F1、F2产生的力矩;
力矩-M(根据作用和反作用)称为陀螺力矩,它是圆盘施加与转轴的力矩,相当于弹性力矩。
在正进动(0<
0<
n/2)时,它是转轴的变形减小,从而提高了转轴的弹性刚度,即提高了转子的临界转速;
在反进动(n/2<
n)时,它是转轴的变形增大,从而降低了
转轴的弹性刚度,即降低了转子的临界转速。
当机械中的高速转动部件的对称轴被迫在空间中改变方位时,即对称轴被迫进动时,转
动部件必须对约束作用一个附加力偶,这种现象称为陀螺效应。
当陀螺效应严重时,可能使机械产生故障,尤其是轴承。
1.4弹性支撑对转子临界转速的影响
Jeffcott转子:
这种转子模型是对真实转子的简化,刚性支承的单盘转子,单盘位于支承的中间,分析临界转速和陀螺力矩等,是转子动力学的基础。
假设盘在平面内运动,不考虑轮盘的偏转,轴是无重轴。
2000000000000.000000000982
1.基本参数:
截面惯性矩:
J=4.91E-10,弹性模量E=2E11,右端质量m3=0.096g=0.1kg,
2.
两个盘的质量m1=m2=0.8kg,.=叶『r22/2=5.73E-4,Id=Ip/2=2.86E-4。
各轴段的传递矩阵:
第一段:
l=0.045m,J=9.82E-10,a11=1.55E-7,a12=a2仁5.16E-6,a22=2.29E-4,
10.0455.16x10$7.73x10*
012.29x10*5.16X10-6
0010.045
0001一
第二段:
l=0.11m,J=9.82E-10,a11=2.26E-6,a12=a21=3.08E-5,a22=5.6E-4,
0.11
3.08x10,
1.128x10》]
5.6x10鼻
3.08x10
第三段:
_10.155.73x10」2.865x10』「
45
017.64x105.73x10
0010.15
第四段:
X0J
901
初始参数列阵为:
01,令X01=1,选取P=2050r/s,用第一段的矩阵乘此矩阵,即
Xk1=1■0.045%1=0
可得此段的终端参数:
纵1=日01=-22.22
Mk1=0
Qk1=°
X°
2=0
日02=-22.22
M02=0
Q02=-RA
[Xk2=-2.4442-1.128x10」RaI
IBk2=-22.22—3.08x10」Ra
Mk2=—0.11RA
一Qk2=—RA
通过转盘的传递矩阵:
第二段的始端参数列阵为:
,用第二段的传递矩阵乘此列阵,得终端参
数:
3
1.2x10’
01
3362x10
*03=—2.4442—1.128x10》Ra〕|e°
3=—22.22—3.08x10*Ra
4
M03=-2.67x10-0.147Ra
=-8.22X106-4.79Ra
第三段的终端参数:
,用此矩阵乘以第二段的终端矩阵得第三段的起始参数:
PXk3=-30.854-2.789x10'
RA〕久3=-513.62-4.18x103raMk3=-1.26x10°
-0.8655Ra
-Qk3
=-8.22x10°
-4.79RA
通过转盘的传递矩阵:
]3.362x106
,用此矩阵乘以第三段的终端矩阵得第四段的起始参数:
「X04=-30.854-2.789x10-5Ra〕I%=-513.62-4.18x10」RaM04=-1.876X106-1.367RA
]Qo4=-1.12X108-98.556Ra一
Xk4=-219.47-2.27x10-4Ra
杠4=—5013.78—4.219x10-3Ra
Mk4=-1.4196x107—7.72Ra
-Qk4=—1.12X108—12.21Ra_
零,从而可得RA=-9.668E5,第五段的起始矩阵为:
X05=0
入5--200
M05二-6.44x106
.Qo5=—8..27X107
=6.524E5.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 临界 转速 计算