汽轮发电机组振动特性及处理对策——云南电力试验研究院(集团)公司穆钢前言汽轮发电机组振动状况对于机组的安全运行有至关重大的影响,轻则产生噪音(旋转机械的谐波干扰),影响机组的安全使用寿命;重则破坏机组的正常运行,损坏机组部件危及安全生产;更为严重者,可能会造成整个机组的灾难性破坏。
汽轮发电机组振动状况是国家综合工业实力的表征;机械制造技术水平的重要标志;发电企业生产管理水平的重要体现。
汽轮发电机组振动处理经历过三个阶段:原始型-经验型-故障诊断型。
汽轮发电机组振动是一项复杂的工程技术问题,其难点在于:同一振动现象,所产生的原因却各不相同,就当代的技术水平而言,还不能确定机组振动在现象-产生原因-处理措施三个方面完全意义上的一一对应关系;同一振动原因对于不同的机组,处理对策却各不相同。
汽轮发电机组各种各样的振动原因,可归纳为三个方面的问题(1)干扰力和干扰力矩的作用;(2)轴承油膜的不稳定和汽隙中的压力不稳定;(3)转子-轴承-基础系统中,支承系统的弹性变形和不稳定性。
在分析、处理机组振动故障的过程中,能否正确判断影响机组振动诸因素中的关键因素是属于哪一个方面的问题,是否有一个正确的原则性判论,是成功的处理机组振动故障的关键之所在。
转子振动的基本特性一、临界转速与转子固有频率相对应的转速,称为转子的临界转速。
数学表达式:ω= √K / MK:转子的刚度;M:转子的质量。
表达式表明:转子的结构和质量决定之后,其自由振动的固有频率就确定了,转子有一系列的固有频率,同理,转子就有一系列的临界转速。
二、共振转子以转速n旋转时,由于不平衡质量所产生的惯性离心力也以转速n的频率施加于转子本身,使转子受到周期性的干扰力。
当干扰力对转子的作用频率与转子的固有振动频率相同时,即ω =ωn ,称为共振。
机组发生共振时,振幅激剧增加振动强烈,当改变转速时,共振引起的振动将随之减小或消失。
工作转速避开共振转速的安全裕量:柔性转子 1.4 n1< n < 0.7n2 ;刚性转子 n1=(1.25~1.8)n 。
三、正进动、反进动考虑静弯曲(偏心矩e)时的单轮盘转子在x、y方向上的振动微分方程式:M d2x+kxx = Meω2 cosωt dt2M d2y+kyx = Meω2 sinωt dt2上式是二阶非齐次方程,其通解为:x=c1cosω1t+c2sinω1t+Axcosωty=c1sinω1t+c2cosω1t+Aysinωt式中的前二项是对应齐次方程(无阻尼情况下)的通解,在有阻尼情况下,这二项将迅速消失,在外干扰力作用下的振动,实际上只有第三项,即:x = A x cosωty = A y sin ωt将x 和y 代入微分方程,解得:如果x 、y 方向的刚度相同,则A X = A y ,于是有:可见,轮盘中心点(S )的运动轨迹是一个圆,圆心在坐标原点0,圆的半径为:A x =Me ω2K x 一M ω2A y =Me ω2K y 一M ω2x =Me ω2cos ωtK 一M ω2y =Me ω2sin ωtK 一M ω2r = √x 2+y 2 = AA =Me ω2K 一M ω2转子的运动实际上可以看成两种运动的合成:一是轮盘绕其轴心(S )等速转动,角速度ω;另一是轮盘中心(S )点又绕其静平衡时的位置0点作圆周运动,矢量OS 的转动角速度也为ω。
上述结论建立了二个重要概念:一是由于偏心离心力的作用,轴被弯成固定的弯弓形状,就象弓绕其弦转动一样称为弓状迴转运动,这样的运动也称为涡动。
二是在转子旋转时,轴及轮盘以角速度ω绕其轴线旋转,这一轴线是转子弹性曲线的切线,而旋转着的弹性曲线在横截面处的位移轨迹,可以是圆,也可以是椭圆。
正进动:弹性曲线以与轮盘相同的角速度和相同的方向旋转,称为正进动。
反进动:弹性曲线以与轮盘相同的角速度旋转,但旋转方向相反,称为反进动。
四、振动位移、速度、加速度振动位移:物体振动时的瞬态幅值大小(通常用x 表示);振动速度:单位时间内物体振动时的瞬态幅值大小(通常用v 表示); 振动加速度:单位时间内物体振动速度的变化率(通常用a 表示)。
数学表达式:x = A 0 cos ωt五、幅频特性从微分方程的推导可知:V =dx = A 0ω cos(ωt+π/2)dta =dV = A 0ω2cos(ωt+π)dtA =Meω2=e (ω2 / KM)=e ( ω2ω12)K一Mω21-(ω2/ KM) 1-( ω2ω12)上式表达的几个极重要的概念,即转子的幅频特性。
幅频特性:(1)e ↘→ A ↘机组做动平衡的意义;(2)ω→ω1时,A ↑;(3)ω=ω1时,A→∞共振的危险;(4)ω>ω1之后,A ↘(逐渐减小),最后→ e , 即振幅放大系数β=A / e= -1,转子自动定心,运转趋于平稳,这是挠性转子可靠运转的原理,也是高频干扰力对于固有频率很低的振动系统,实际上不产生强迫振动的道理。
(5)ω<ω1时,干扰力与振动位移是同方向;(6)ω>ω1时,干扰力与振动位移是反方向,也就是说,在转速变化过程中,偏心激振力和振动位移的相差有一突变。
幅频特性曲线:βω/ω1六、共振方程ω=ω1时,方程式M d2x+kxx = Meω2 cosωt dt2变为:d2x+ω12 x = eω12 cosωt (1)dt2特解:x = At cos( ω1t + a )将其代入方程式(1),求得满足方程式(1)的A和a:A = e ω1/ 2 ; a = -π/2 .于是,得到方程式(1)的特解:x = e ω1 t cos( ω1t +π/2) / 2上式称为共振方程,其意义:(1)共振时,振幅是随时间t的延续而增大的,不论e多么小(干扰力多么小),机组在共振转速停留是很危险的。
(2)在启动机组时,当升速到接近临界转速时,应加快增速,在临界转速附近应很快冲过。
这样,转子的振幅还来不及增大就已超过临界转速,使转子进入自动定心的平稳运行过程。
(3)解释了实测临界转速高于计算临界转速的原因。
(4)解释了现场对同一台机组多次实测临界转速,在数值上存在着差别的七、阻尼力、机械滞后角、相频特性阻尼对转子振动的影响:转子振动时,通常存在着各种阻碍运动的阻力,这些阻力称为阻尼,它对转子振动起着衰减和抑制作用。
阻尼有不同的来源,如:(1)干摩擦;(2)材料内摩擦;(3)粘滞阻尼。
只有在粘滞阻尼下的強迫振动才是简谐运动,才可用简单的数学方程作详细分析。
单自由度有阻尼強迫振动的微分方程式:M d2x=-kx- cdx+Meω2cosωt dt2dt(惯性力、阻尼力和弹性恢复力与偏心不平衡离心力相平衡构成稳定系统。
)d2x= - kX -c*dx+eω2cosωtdt2M M dt令:2D=C/M ; ω12 = K/M式中:C: 阻尼比例系数,它表示速度为1个单位时所受到阻力的大小,用公斤.秒/厘米表示。
D:阻尼系数,表示单位时间内能量消耗率之一半。
微分方程由下式所代替:d2x+2D dx+ω12 x = eω2cosωt (1)dt2dt上式微分方程的通解,应是本微分方程的一个特解与下列齐次方程通解的化数和。
d2x+2D dx+ω12 x = 0 (2)dt2dt在有阻尼的情况下,齐次方程(2)式的通解将会消失,只剰下方程(1)式的特解这一部分,特解:x = A COS(ωt - φ) 将x 带入微分方程解得:A =eω2ω12[ ( 1-ω2)2 + 4 D2ω2]1/2ω12ω14φ = tg-1 (2Dω2)ω121-ω2ω12tgφ =2Dω2ω12 -ω2A:振动过程中的实际最大振幅;φ:机械滞后角,它表示机械振动中,由于惯性效应的存在,振幅(位移)始终滞后于引起振动的扰动力F一个角度,该角度和振动系的自振频率、ω/ω1的比值及系统的阻尼有关。
分析A的数学表达式,可得出如下结论:共振时,ω=ω1,共振振幅是 A = e / (2D/ω1)= e /2 rc, rc=D/ω1称为相对阻尼系数,在有阻尼的情况下,转子在共振转速时,其振幅不是无穷大,而是一有限值。
相频特性:(1)ω<ω时,φ的变化范围在0—900之间(tgφ为正);1(2)ω=ω1时,φ = 90O(tgφ为∞);(3)ω>ω1时,φ的变化范围在900—1800之间(tgφ为负)。
相频特性曲线:ππ/ω/ωn(φ-ω/ωn 关系曲线)D:阻尼系数,表示单位时间内能量消耗率之一半。
φ-ω/ωn 关系曲线表明:ω<ωn D大—φ大;ω=ωn φ与D无关;ω>ωn D大—φ小或D小—φ大。
做平衡时:重型机组D大(一般设计为16),φ取小点;轻型机组D小,φ取小点。
注:D是状态参数,是一个变量,将系统阻尼D作为非线性影响因素之一分析考虑,就不难解释影响系数对同一台机组都不能通用,更不用说对同型机组通用的原因。
八、力平衡、旋转矢量平面汇交力系:各力作用线位于同一平面,且汇交于一点。
平面汇交力系平衡的几何条件:力之多边形成闭合ΣF=0,或闭合积分为0。
力平衡:按平面汇交力系力平衡条件建立的微分方程式——惯性力(Mω2x)、弹性恢复力(kx)、阻尼力(cωx)与由于偏心矩的存在而产生的不平衡离心力(F= Me ω2)相平衡,即、稳定振动系统的力平衡条件,否则就是一个发散系统。
旋转矢量:位移、速度、加速度的图形表示 x ω2xωx(1)ω < ω12= x * k/M∵ xω2 < x ω1∴ Mω2x < xk2x Kx FCωMω2x(2) ω = ω12x∵ Mω2x = Mω12x = Mx*k/M = Mω2x∴ Kx = Mω1F Kx FCωxω2x Mωx(3) ω > ω1显然 M ω2x > kx = M ω12xFKxM ω2M ω2八、阻力特性(阻力对振动的影响)A =eω2ω12 [ ( 1-ω2)2 + 4 D 2ω2 ]1/2ω12ω14β=1 [ ( 1-ω2 )2 + 4 r c 2ω2 ]1/2ω12ω12r c :相对阻尼系数。
式中: D 2 = r c 2ω12微分方程特解:x = A COS (ωt - φ) x =βA 0 COS (ωt - φ)β: 振动过程中的振幅放大系数(β=A/A 0)。
A 0 = eω2=Me ω2ω12kA:转轴在偏心离心力F= Meω2作用下的静挠度。
分析β得知:ω/ω1很小时,β→ 1;ω/ω1很大时,β→ 0;ω/ω1→ 1时,β↑↑,此时阻尼对β的影响很大,随着rc→0,β→∞。
βω/ωn阻尼特性:(1)ω << ω1 和ω>>ω1时,阻尼对振动的影响都很小。