转子轴承系统失稳转速及稳定裕度的研究.docxVIP

转子轴承系统失稳转速及稳定裕度的研究 近年来，旋转动力学主要研究领域之一是旋转轴承系统的稳定性。评价稳定最直接的方法是通过计算和实验考察系统的失稳转速。我们使用了两种模型进行计算并对计算结果进行了比较。第一种是用多质量模型计算系统阻尼过零点的转速(即失稳转速);第二种是用单质量模型计算不同转速下的系统阻尼。计算结果表明,两种模型计算出的失稳转速值十分吻合。本文主要说明后者。 1 转子两端轴颈的几何模型 汽轮发电机组的发电机转子及其实验模型的刚度质量分布基本上是对称的。支持转子两端轴颈的轴承在结构形式、几何参数及运行条件上也完全一样,在这种情况下,图1所示的单质量数学模型会比一般的质量、刚度分布不对称转子或两端支撑轴承不相同的情况有更好的精度。 2.1 油膜刚度t结构 对图1所示的转子轴承系统简化模型,集中质量及轴颈运动的动力学方程可以写成下列矩阵的形式: [m0000m0000m0000m]?{??x??y??xB??yB}+[0000000000BxxBxy00ByxByy]?{˙x˙y˙xB˙yB}+[ks0-ks00ks0-ks-ks0ks+kxxkxy0-kskyxks+kyy]?{xyxByB}=0(1) 式中: {2m转子质量ks转子的等价刚度k轴承油膜刚度系数B轴承油膜阻系数 ;设(1)的解为 {xB=XBeυωtyB=ˉYBeυωt{x=Xeυωty=Yeυωt(2) 其中: {υ=ξ+jηω=转子的角速度j=√-1 ;则: {˙xB=XBυωeυωt??xB=XBυ2ω2eυωt{˙yB=YBυωeυωt??yB=YBυ2ω2eυωt{˙x=Xυωeυωt??x=Xυ2ω2eυωt{˙y=Yυωeυωt??y=Yυ2ω2eυωt 代入(1)并整理,我们有∶{(mω2υ2+k˙s)X-ksXB=0(mω2υ2+ks)ˉY-ksˉYB=0(3){-ksX+(ks+kxx+υωBxx)XB+(kxy+υωBxy)ˉYB=0-ksX+(ks+kyy+υωByy)ˉYB+(kyx+υωByx)XB=0(4) 将(3)代入(4),并整理: {(-k2smυ2ω2+ks+ks+kxx+υωBxx)XB+(kxy+υωBxy)ˉYB=0(kyx+υωByx)XB+(-k2smυ2ω2+ks+ks+kyy+υωByy)ˉYB=0(5) 欲使(5)式有非零解,则需: |-ksmυ2ω2+ks+ks+kxx+υωBxxkxy+υωBxykyx+υωByx-ksmυ2ω2+ks+ks+kyy+υωByy|=0(6) 将(6)式展开并整理:D′6υ6+D′5υ5+D′4υ4+D′3υ3+D′2υ2+D′1υ+D′0=0 (7) 其中: {D′0=ks2c′0D′1=ks2ωc′1D′2=ks2ω2c′2+2mksω2c′0+mks2ω2c′4D′3=2mksω2c′1+mks2ω3c′3D′4=m2ω4c′0+2mksω4c′4+m2ksω2c′4+m2ks2ω4D′5=m2ω5c′1+m2ksω5c′3D′6=m2ω6c′2{c′0=kxxkyy-kxykyxc′1=kxxByy+kyyBxx-kxyByx-kyxBxyc′2=BxxByy-BxyByxc′3=Bxx+Byyc′4=kxx+kyy 将(7)无量纲化:D′0=ks2c′0=ks(kxxkyy-kxykyx)=ks2(μωl/ψ3)2(kˉxxkˉyy-kˉxykˉyx)=ks2(μωl/ψ3)2c0 其中:ψ=cR;c0=(kˉxxkˉyy-kˉxykˉyx);C—轴承半径间隙;R—轴承半径;μ—润滑油粘度;L—轴承长度; kˉ—无量纲油膜刚度系数。 (7)式两边同除以ks2(μωl/ψ3)2得:D6υ6+D5υ5+D4υ4+D3υ3+D2υ2+D1υ+D0=0 (8) (8)式即为单质量转子轴承系统模型的特征方程。 其中: {D0=c0D1=c1D2=c2+(ωωk1)2(2c0+cωk12gsfc4)D3=(ωωk1)2(2c1+cωk12gsfc3)D4=2(ωωk1)2c2+(ωωk1)4?(c0+cωk12gsfc4+(cωk12gsf)2)D5=(ωωk1)4(c1+cωk12gsfc3)D6=(ωωk1)4c2 c0=kxxkyy-kxykyxc1=kxxByy+kyyBxx-kxyByx-kyxBxyc2=BxxByy-BxyByxc3=Bxx+Byyc4=kxx+kyyωk1=2πks2m= 刚性支承一临界转速(RAD/s) sf=mgμlRω g=重力加速度 1.2 构造约束方程 由于kˉ、Bˉ、ω、ωk1、sf等均为实数,式(8)中的各系数Di(i=0,1,2,…6)亦均为实数。由实系高次代数方程根的理论可知,特征方程(8)