文章快速检索     高级检索
  北京化工大学学报(自然科学版)  2017, Vol. 44 Issue (6): 81-86   DOI: 10.13543/j.bhxbzr.2017.06.013
0

引用本文  

马波, 杨灵, 沈伟. 透平机械转子系统不对中故障定量诊断监测方法研究[J]. 北京化工大学学报(自然科学版), 2017, 44(6): 81-86. DOI: 10.13543/j.bhxbzr.2017.06.013.
MA Bo, YANG Ling, SHEN Wei. A method for quantitative diagnosis and monitoring of rotor misalignment faults in turbine machinery[J]. Journal of Beijing University of Chemical Technology (Natural Science), 2017, 44(6): 81-86. DOI: 10.13543/j.bhxbzr.2017.06.013.

基金项目

压缩机技术国家重点实验室开放基金(gjzdsys2014)

第一作者

马波, 男, 1977年生, 副研究员.

通信联系人

沈伟, E-mail:shenweibuct@163.com

文章历史

收稿日期:2017-03-16
透平机械转子系统不对中故障定量诊断监测方法研究
马波 1, 杨灵 1, 沈伟 2     
1. 北京化工大学 机电工程学院 诊断与自愈工程研究中心, 北京 100029;
2. 环境保护部核与辐射安全中心, 北京 100082
摘要:针对透平机械转子系统不对中故障程度难以定量诊断的问题,提出一种基于转子的状态监测数据建立转子平行不对中量与节点振动通频幅值的关系曲线来定量诊断转子平行不对中故障的方法。对平行不对中状态下的转子进行受力分析得到转子不对中激励力,通过建立的转子试验台有限元数值分析模型计算得到不同不对中激励力下的转子节点振动通频幅值,进而通过分析转子节点振动通频幅值与不对中量的对应关系曲线及转子运行正常数据对关系曲线进行修正;根据监测得到的转子节点振动通频幅值反推出转子平行不对中量,达到对平行不对中故障定量诊断的目的。最后试验验证了定量诊断方法的可行性。
关键词转子平行不对中    定量诊断    故障机理    状态监测    
A method for quantitative diagnosis and monitoring of rotor misalignment faults in turbine machinery
MA Bo1 , YANG Ling1 , SHEN Wei2     
1. Diagnosis and Self-Recovery Engineering Research Center, College of Mechanical and Electrical Engineering, Beijing University of Chemical Technology, Beijing 100029;
2. Nuclear and Radiation Safety Center, MEP, Beijing 100082, China
Abstract: It is difficult to diagnose quantitatively the degree of misalignment in a turbine rotor system. A quantitative method for diagnosis of parallel misalignment based on monitoring data, utilizing the relationship between the misalignment of the rotor and the vibration amplitude of the node, is proposed in this work. The steps involved in the method are as follows. The rotor misalignment force is obtained by analyzing the state of the parallel misalignment rotor. A rotor model is then established to analyze the variation of rotor node vibration for different misalignments based on dynamics theory. The relationship between the vibration amplitude of the node with the misalignment is then analyzed, and modified according to the data obtained with the normal operation of the rotor. By monitoring the vibration amplitude of the node, the parallel misalignment of the rotor can be derived, and the aim of diagnosing the faults in the parallel misalignment is thus achieved. Using a test rig, the vibration of the rotor with different misalignments was analyzed by experiment. The experimental results show that the proposed method of monitoring and assessment is valid.
Key words: rotor parallel misalignment    quantitative diagnosis    fault mechanism    monitoring    
引言

转子不对中是透平机械的频发故障之一,尤其是大型多轴系的汽轮发电机组和透平压缩机组,准确评估不对中故障的严重程度即转子不对中量的大小,对于制定检维修措施和优化生产计划排程具有重要意义,其中诊断与定量评估不对中故障的严重程度尤为关键。

对旋转机械进行定量诊断的方法主要有3种。一是基于运行状态的故障诊断,如姚红良等[1]提出了基于谐波分量的转子系统碰摩故障定量诊断方法,鞠华等[2]针对轴承的故障定量诊断使用支持向量回归的方法进行研究,Shi等[3]提出了基于电机定子电流信号的小波系数频带分析方法精确诊断电动机转子断条故障。这些方法通过建立振动信号与故障特征的关系进行故障诊断,但振动分析的信号是针对特定测点在瞬间采集的一段波形,如果发生故障,在某一状态下,故障信息会被淹没在其他信息中,因此基于某时刻状态波形的方法就不能很好地区分这些故障。二是基于运行状态变化过程的诊断方法,如陈非等[4]提出针对旋转机械振动故障基于过程的定量诊断方法,通过采用多通道多转速下的信息熵矩阵来描述振动过程的变化规律以达到定量诊断的目的;费成巍等[5]提出了利用转子振动故障的过程功率谱熵特征进行定量诊断的技术。这些方法通过分析整个振动过程的数据可以弥补状态分析方法的不足,但是难以准确判断故障发生的程度。三是基于模型的诊断方法,如于涛等[6]针对转子碰磨故障提出了基于模态扩展与谐波分解的诊断方法,采用模态扩展的方法对监测得到的机组测点的振动量进行处理,得到机组各个部分的振动量;Markert等[7]利用模态扩展法得到的旋转机械各部分振动值来推测得到故障力大小,利用最小二乘法进行计算确定等效力的位置;Zhao等[8]提出一种基于模态分解和近似熵值法的滚动轴承剥落故障的定量诊断方法,实现了滚动轴承剥落故障的定量诊断。这些方法通过计算转子的模态参数,利用模态扩展的方法估算故障力来确定故障位置,但通过力进行故障程度判断对大型设备及机组的诊断准确性不高。

综上所述,目前对旋转机械故障尤其是轴承故障和转子不平衡、不对中、碰摩等故障的定量诊断方法研究已经取得了较大的进展,但仍存在诊断准确度不高的问题。因此本文提出基于转子动力学模型和透平机械转子状态监测数据的转子平行不对中故障定量诊断方法,先进行转子故障类型判断,确定为平行不对中故障后再根据本文方法进行定量诊断,最后通过试验对该方法的可行性进行了验证。该方法避免了直接理论推导得到的数据与实际运行数据相差大、难以匹配且诊断结果不准确的问题。

1 转子平行不对中动力学模型理论及受力分析

轴承座及联轴器结合处将对转子系统产生激振力,同时考虑不平衡和重力的作用,可以得到转子系统具有不对中故障时的微分方程[9]

$ \mathit{\boldsymbol{M}}\ddot a + \mathit{\boldsymbol{C}}\dot a + \mathit{\boldsymbol{K}}a = \mathit{\boldsymbol{Q}}\left( t \right) + {\mathit{\boldsymbol{F}}_{\rm{r}}} + \mathit{\boldsymbol{G}} $ (1)

其中,M为系统质量矩阵,C为系统刚度矩阵,K为系统阻尼矩阵,Q为不平衡力,Fr为广义不对中力,G为重力,a为位移。简化的动力学模型如图 1所示。

图 1 转子平行不对中故障系统动力学模型 Fig.1 Dynamic model of the rotor misalignment model

根据牛顿第二定律,转子的运动方程为:

$ \left\{ \begin{array}{l} {m_1}{{\ddot x}_1} - {k_1}\left( {{x_1} - {x_2}} \right) + {k_1}\left( {{x_1} - {x_3}} \right) + {c_1}{{\dot x}_1} = \\ \;\;\;\;\;\;\;{F_x} + {m_1}{\omega ^2}e\cos \omega t\\ {m_1}{{\ddot y}_1} - {k_1}\left( {{y_1} - {y_2}} \right) + {k_1}\left( {{y_1} - {y_3}} \right) + {c_1}{{\dot x}_1} = \\ \;\;\;\;\;\;\;{F_y} + {m_1}{\omega ^2}e\sin \omega t\\ {k_1}\left( {{x_1} - {x_2}} \right) = 2{k_2}{x_3} + 2{c_2}{{\dot x}_3}\\ {k_1}\left( {{y_1} - {y_3}} \right) = 2{k_2}{y_3} + 2{c_2}{{\dot y}_3} \end{array} \right. $ (2)

其中,ml为转盘的质量,k1c1为转盘的刚度和阻尼;k2c2为右侧轴承的刚度和阻尼;k3c3为左侧轴承的刚度和阻尼;联轴器的激振力在xy方向上的分力为FxFy

2 平行不对中故障诊断方法 2.1 诊断步骤

本文提出的平行不对中故障定量诊断方法是基于转子动力学模型分析计算和转子平行不对中试验,得到转子节点振动通频幅值与平行不对中量的对应关系,并结合透平机械转子状态监测数据来进行平行不对中故障定量诊断。转子平行不对中故障定量诊断的主要步骤如下:

(1) 对有平行不对中故障的转子进行受力分析得到转子不对中力;

(2) 建立转子试验台的有限元数值分析模型,计算模型节点在不对中激励力下的振动响应;

(3) 分析建立平行不对中量不同时,转子节点运行的各幅值与平行不对中量的关系曲线;

(4) 利用不对中试验数据来修正模拟得到的关系曲线,比较修正后模拟曲线与试验曲线的相关性,找出相关性最好的曲线;

(5) 根据得到的曲线及监测得到的转子运行数据反推得到转子的平行不对中量。

2.2 转子平行不对中激励力计算

不对中分为联轴器不对中和轴承不对中两种,而联轴器不对中又分为平行位移不对中、角度位移不对中以及混合不对中3种[10]。平行不对中理论上包含两种情况:一是安装时法兰中心错开不在同一轴线上,文献[11]中针对这种情况推导出不对中公式模型;二是安装时对中良好,后期由于支承松动导致的不对中使得法兰错开,造成螺旋或轴弯曲。当联轴器的两转子轴线之间发生不对中故障时,两个半联轴器均绕自己的中心转动。

研究表明,出现平行不对中时联轴器外壳的轴心线运动轨迹为一圆柱体[12]。进行联轴器受力分析时将连接位置看作刚性,即无任何变形,则在不对中时转子运动过程中发生变形的是联轴器法兰以及转轴。

图 2为平行不对中故障模型。图中O1O2为两个半联轴器的中心,当联轴器平行不对中时,连接螺栓为P,从动轴与法兰连接点为A,初始状态时AO2重合,不对中量ey

图 2 平行不对中故障模型 Fig.2 The parallel misalignment model

旋转角度α=ωt后的偏心量为e′,同时P旋转到P′,A旋转到A′,O2旋转到O2′。由于已假设连接螺栓P为刚性,则联轴器的力由从动轴与法兰变形产生,联轴器受力所对应的变形量为AO2之间的距离在半径方向的投影,即

$ e = P{O_1} - P{O_2} $ (3)
$ e' = P'{O_1} - P'{O_2} $ (4)

因为PO2′=PO2,且e一般都很小,故可近似认为O2O2′垂直于PO1,即

$ \left\{ \begin{array}{l} e' = P'{O_1} - P'{O_2} = P'{O_1} - P'{{O'}_2} = {{O'}_2}{O_1}\\ {{O'}_2}{O_1} = {O_2}{O_1}\cos \alpha = e\cos \alpha \end{array} \right. $ (5)
$ {O_1}{O_2} = {O_1}A' $ (6)

则转动角度α后偏心距的变化量为

$ \Delta e = e - e' = e - e\cos \alpha $ (7)

将公式(5)、(6)带入公式(7),得到

$ \Delta e = e - e\cos \alpha = {O_1}A' - {{O'}_2}{O_1} = A'{{O'}_2} $ (8)

即变形量为AO2之间距离在半径方向的投影。

设联轴节和轴的刚度为k,则FxFy由公式(9)计算:

$ \left\{ \begin{array}{l} F = k\Delta e = ke - ke\cos \alpha \\ {F_x} = F\cos \alpha = ke\cos \alpha - ke{\cos ^2}\alpha = \\ \;\;\;\;\;\;\;\;ke\cos \alpha - \frac{1}{2}ke\cos 2\alpha + \frac{1}{2}ke\\ {F_y} = F\sin \alpha = ke\sin \alpha - ke\cos \alpha \sin \alpha = \\ \;\;\;\;\;\;\;\;ke\sin \alpha - \frac{1}{2}ke{\sin ^2}2\alpha \end{array} \right. $ (9)

式(9)中k取3.6×105N/m,e取0、0.05、0.15、0.20cm。

在计算动力学响应时,将联轴器的激振力施加于联轴器所在节点位置来代替不对中故障下的受力[13],计算得到转子节点在不同e下的振动响应。

2.3 有限元数值分析模型建立

在单跨转子系统中,电机是系统的驱动部件,转轴是被驱动部件,两者之间通过联轴器相连。对转子试验台建立有限元数值分析的转轴由beam188单元代替,圆盘由Mass21单元代替,轴承由Com214单元代替。转轴和轴承的单元参数设置如表 1所示,转子模型在网格划分后的有限元数值分析模型如图 3所示。

下载CSV 表 1 转轴和轴承单元参数 Table 1 Setting of the rotor parameters
图 3 转轴有限元模型 Fig.3 Finite element model of the rotating shaft
2.4 响应结果及分析

动力学响应分析的计算参数设置为:转子转速6000r/min,求解时间0.2s,求解步长1000步。计算出节点在y方向位移随时间的变化关系,选取传感器所在位置节点振动响应的通频幅值分析,同时对该节点的数据进行傅立叶变换得到节点的频率变化。e为0.5mm时节点振动值随时间及节点频率变化情况如图 4所示。可以看出在转子出现平行不对中故障时,节点振动幅值随着时间的推移呈现周期性变化,且出现了二倍频幅值较高的情况。

图 4 e=0.5mm时节点振动示意图和节点频谱图 Fig.4 The vibrational waveform and the spectrum of the node with a misalignment of 0.5mm

分别模拟e为0.150mm、0.200mm、0.250mm、0.325mm、0.425mm和0.500mm时转轴的节点振动情况,得到节点在不同e下的各幅值变化如表 2所示。可以看出随着e的增大,各幅值的大小也相应的增加,其中通频幅值的增幅最大,二倍频次之,一倍最小。出现平行不对中故障后各幅值呈相似变化规律。

下载CSV 表 2 不对中量与各幅值关系(模拟) Table 2 Corresponding relationship between fault quantity and amplitude value (simulation)

由于不对中激励力为周期力,可以通过多项式对其进行拟合。6种不对中量下的模拟节点各幅值采用二次多项式拟合得到的公式如表 3所示,拟合曲线如图 5所示,决定系数R2由公式(10)计算。由图 5中点线关系及表 3R2的大小可以看出,拟合公式准确度很高,拟合效果好。

下载CSV 表 3 各关系拟合公式 Table 3 Relationship formulas
图 5 各幅值与e的关系曲线(模拟) Fig.5 The change of amplitude value with e (simulation)
$ {R^2} = \frac{{{{\left( {\sum\limits_{i = 1}^n {\left( {{x_i} - {{x'}_i}} \right)\left( {{y_i} - {{y'}_i}} \right)} } \right)}^2}}}{{\sum\limits_{i = 1}^n {{{\left( {{x_i} - {{x'}_i}} \right)}^2}} \sum\limits_{i = 1}^n {{{\left( {{y_i} - {{y'}_i}} \right)}^2}} }} $ (10)
3 试验验证 3.1 试验台及相关参数

本文所用转子试验台为电机带动单转子运行,电机与转子由刚性联轴器连接。采用电涡流传感器测量转子的振动位移,传感器灵敏度为4V/mm,采样频率5120Hz,采样点数2048,放置在距离轴承支座10mm的位置,此位置为从该试验台有限元模型中选取的分析节点。

3.2 试验方案

通过调节转盘所在转轴的高度模拟平行不对中故障,在转盘所在轴的2个支撑轴承的4个底部同时添加相同厚度的垫片,添加垫片的过程中要确保轴承支座仅有上下的偏移而不出现左右偏移,同时通过螺栓连接的联轴器两半也仅出现上下的错动,不出现左右的错动。

共设计4组试验,分别是1组正常试验和3组平行不对中试验,其中e分别为0.05mm、0.15mm和0.20mm。正常试验中连接转子所在轴和电机所在轴的两片刚性联轴器完全重合,转子转速设置为3000r/min,与动力学响应分析时的转速一致。不对中试验中由传感器采集到不同e下转子节点的振动值,通过分析转子节点在不同平行不对中故障下的振动变化情况及其与e的变化关系,验证定量诊断平行不对中故障方法的可行性和准确性。

3.3 结果与讨论

4次试验转子上节点的振动通频、一倍频和二倍频的幅值如表 4所示。随着e逐渐增大,通频幅值、一倍频幅值和二倍频幅值也随着增大,并且通频幅值变化最大,二倍频次之,一倍频最小,与模拟结果的变化规律基本相同。虽然一倍频幅值出现了先上升后下降的情况,相对于无故障时幅值仍是上升的,原因可能是一倍频在不对中量的一定范围内更加敏感。

下载CSV 表 4 不对中量与各幅值关系(试验) Table 4 Relationship between fault quantity and amplitude value (experiment)

对试验得到的各幅值也采用二次多项式进行拟合,拟合得到的各幅值与e的关系曲线如图 6所示,可以看出拟合度很高。

图 6 各幅值与e的关系曲线(试验) Fig.6 Change of amplitude value with e (experiment)

在模拟计算中,转子在正常情况下没有不对中激励力,故节点无初始振动值;而试验中转子在正常不对中可忽略的情况下,转轴由于转动会出现振动,故在e=0时的模拟与试验得到的通频幅值出现偏差,需要进行修正;而一倍频和二倍频的幅值没有出现偏差,无需修正。通过添加试验过程中正常状态下的初始通频幅值得到修正后的通频幅值随e变化的拟合公式:

$ y = 1.19{x^2} + 427.33x + 126.34 $ (11)

将试验中的e带入模拟拟合得到的公式中进行计算,得到不同e在模拟时的各幅值,再对模拟数据和试验数据得到的曲线进行比较。将各幅值的试验得到的关系曲线与模拟得到的3条关系曲线进行对比可知,通频幅值的模拟曲线与试验曲线相关度最高,其幅值大小相差最小,故利用e与通频幅值关系曲线进行不对中量定量诊断更为准确。修正后的通频幅值模拟数据和试验数据曲线如图 7所示。

图 7 试验曲线与修正曲线对比图 Fig.7 Relationship between the fault and the two frequency multiplication

图 7看出,建模结果与试验结果的振动值相差小于10%(表 5),误差在可接受范围内,由此验证了转子平行不对中故障下节点振动通频幅值与不对中量关系曲线的准确性,以及利用关系曲线和节点振动通频幅值进行平行不对中故障定量诊断的可行性。

下载CSV 表 5 不同e下试验与模拟的通频幅值误差 Table 5 Experimental and simulated vibration values of nodes for different values of e
4 结束语

本文提出了一种基于转子动力学模型和旋转机械监测运行状态的转子平行不对中故障的定量诊断方法,并试验验证了该方法的可行性。通过实时监测得到的转子节点振动数据以及由建模分析并修正的节点振动通频幅值随不对中量的关系曲线,本文方法可对转子平行不对中故障程度进行定量诊断,实现在监测转子运行过程中对平行不对中故障的定量监测,同时能提高平行不对中故障定量诊断的准确性。

参考文献
[1]
姚红良, 韩清凯, 李凌轩, 等. 基于谐波分量的转子系统碰摩故障定量诊断方法[J]. 机械工程学报, 2012, 48(5): 43-48.
Yao H L, Han Q K, Li L X, et al. Method for detecting rubbing fault in rotor system based on harmonic components[J]. Journal of Mechanical Engineering, 2012, 48(5): 43-48. (in Chinese)
[2]
鞠华, 沈长青, 黄伟国, 等. 基于支持向量回归的轴承故障定量诊断应用[J]. 振动、测试与诊断, 2014, 34(4): 767-771.
Ju H, Shen C Q, Huang W G, et al. Application of bearing fault diagnosis based on support vector regression[J]. Journal of Vibration, Measurement and Diagnosis, 2014, 34(4): 767-771. (in Chinese)
[3]
Shi P, Chen Z, Vagapov Y, et al. A new diagnosis of broken rotor bar fault extent in three phase squirrel cage induction motor[J]. Mechanical Systems and Signal Processing, 2014, 42(1): 388-403.
[4]
陈非, 黄树红, 张燕平, 等. 基于过程的旋转机械振动故障定量诊断方法[J]. 动力工程, 2008, 28(4): 543-547.
Chen F, Huang S H, Zhang Y P, et al. Aquantitative diagnosis method of vibration faults of rotating machinery based on process[J]. Journal of Power Engineering, 2008, 28(4): 543-547. (in Chinese)
[5]
费成巍, 柏树生, 白广忱, 等. 转子振动故障的过程功率谱熵特征分析与定量诊断[J]. 航空发动机, 2011, 37(1): 48-51.
Fei C W, Bai S S, Bai G C, et al. Feature analysis and quantitative diagnosis of process power spectrum entropy for rotor vibration failure[J]. Aeroengine, 2011, 37(1): 48-51. (in Chinese)
[6]
于涛, 韩清凯, 孙伟, 等. 基于模态扩展与谐波分解的转子碰摩故障精确诊断[J]. 东北大学学报:自然科学版, 2006, 27(5): 524-527.
Yu T, Han Q K, Sun W, et al. Rotor rub fault diagnosis based on modal expansion and harmonic decomposition[J]. Journal of Northeastern University:Natural Science, 2006, 27(5): 524-527. (in Chinese)
[7]
Bach H, Hiller R, Markert R. Representation of rotor-stator-rub in terms of equivalent forces for model based diagnostics[J]. Acoustical and Vibratory Surveillance Methods and Diagnostic Techniques, 1998, 3: 723-732.
[8]
Zhao S F, Liang L, Xu G H, et al. Quantitative diagnosis of a spall-like fault of a rolling element bearing by empirical mode decomposition and the approximate entropy method[J]. Mechanical Systems and Signal Processing, 2013, 40(1): 154-177. DOI:10.1016/j.ymssp.2013.04.006
[9]
韩捷. 齿式联接不对中转子的故障物理特性研究[J]. 振动工程学报, 1996, 9(3): 297-301.
Han J. Study on the physical characteristics of the fault of the rotor with misalignment[J]. Journal of Vibration Engineering, 1996, 9(3): 297-301. (in Chinese)
[10]
刘占生, 赵广, 龙鑫. 转子系统联轴器不对中研究综述[J]. 汽轮机技术, 2007, 49(5): 321-325.
Liu Z S, Zhao G, Long X. Survey of the research on coupling with misalignment of rotary machinery[J]. Turbine Technology, 2007, 49(5): 321-325. (in Chinese)
[11]
Mehrjou M R, Mariun N, Marhaban M H, et al. Rotor fault condition monitoring techniques for squirrel-cage induction machine-a review[J]. Mechanical Systems and Signal Processing, 2011, 25(8): 2827-2848. DOI:10.1016/j.ymssp.2011.05.007
[12]
陈果, 李兴阳. 航空发动机整机振动中的不平衡-不对中-碰摩耦合故障研究[J]. 航空动力学报, 2009, 24(10): 125-132.
Chen G, Li X Y. Study on imbalance-misalignment-rubbing coupling faults in aero engine vibration[J]. Journal of Aerospace Power, 2009, 24(10): 125-132. (in Chinese)
[13]
高兴军, 赵恒华. 大型通用有限元分析软件ANSYS简介[J]. 辽宁石油化工大学学报, 2004, 24(3): 94-98.
Gao X J, Zhao H H. Large general finite element analysis software ANSYS[J]. Journal of Liaoning University of Petroleum & Chemical Technology, 2004, 24(3): 94-98. (in Chinese)