基于阶比分析的转子系统油膜涡动故障诊断研究

第54卷第3期2012年6月汽轮机技术TURBINETECHNOLOGYVO154NO3JUN2012基于阶比分析的转子系统油膜涡动故障诊断研究胡爱军，朱瑜华北电力大学机械工程系，保定071003摘要针对转子升速过程中油膜涡动信号的非平稳特点，将阶比分析应用于油膜涡动故障诊断中。首先对转子升速时的振动信号进行时域采样，再对时域信号进行等角度采样，将其转化为角域平稳信号，然后对角域信号进行频谱分析识别油膜涡动故障频率。对实测油膜涡动信号分析表明，阶比分析能有效地避免传统频谱分析方法中的“频率模糊”现象，能够准确识别转子升速过程中油膜涡动故障频率。关键词油膜涡动；阶比分析；频率模糊分类号TK288文献标识码A文章编号1001588420120343223433ANALYSISOFOILWHIRLDIAGNOSISBASEDONORDERANALYSISFORROTORSYSTEMHUAIJUNZHUYUMECHANICALENGINEERINGDEPARTMENT，NOAHCHINAELECTRICPOWERUNIVERSITY，BAODING071003，CHINAABSTRACTINORDERTOPROCESSTHENONSTATIONARYOILWHIRLSIGNALTHATISPRODUCEDDURINGROTORSSPEEDINGUP，THEORDERANALYSISISUSEDFORTHEOILWHIRLDIAGNOSISFIRST，THEVIBRATIONSIGNALSATSPEEDINGUPISSAMPLEDATCONSTANTTIMEINCREMENTSINTIMEDOMAIN，ANDTHENTHEDATAISAMPLEDATCONSTANTANGLEINCREMENTS，THETIMEDOMAINSIGNALISCHANGEDINTOANGLEDOMAINSIGNAL，THEN，THEANGLEDOMAINSIGNALWILLBEANALYZEDWITHSPECTRUMTODETECTTHEOI1WHIRLDIAGNOSISEQUENCYTHEPRACTICALOILWHIRLSIGNALISANALYZEDANDRESULTSSHOWTHATORDERANALYSISCANAVOIDFREQUENCYAMBIGUITYOFTHETRADITIONALSPECTRUMSANDORDERANALYSISCANIDENTIFYOILWHIRLDIAGNOSISEQUENEYACCURATELYKEYWORDSOILWHIRL；ORDERANALYSIS；FREQUENCYAMBIGUITY0前言汽轮机等大型旋转机械通常采用滑动轴承来支承转子，在转子升速过程中，通常会发生油膜涡动的故障，使系统不能正常工作。及时准确地识别转子升速过程中的油膜涡动故障，并采取相应的处理措施，对确保转子系统安全稳定运行具有重要意义。长期以来，对旋转机械的故障诊断，主要采用振动信号频谱分析的方法，通过分析振动的频谱图来判断转子系统故障。转子升速过程中油膜涡动发生时，时域采样信号为典型的非平稳信号，直接对时域采样信号进行频谱分析，就会产生“频率模糊”现象J。文献3指出，对于旋转机械因转速变化而引起的非平稳振动信号，采用阶比分析要比小波分析、短时傅立叶变换等分析方法更具优势。而阶比分析的关键在于获得等角度增量的采样信号。传统的等角度采样依靠硬件，不仅成本较高，而且较为繁琐。文献4提出了一种基于瞬时频率估计的阶比跟踪方法，通过对原始信号的时频分布进行峰值搜索估计转子的瞬时频率，依靠转子的瞬时频率对时域信号进行等角度采样，取得了较好的效果。为此，本文将采用基于瞬时频率估计的方法实现对首先通过小波时频分析获得时域信号的小波时频分布图，并通过峰值搜索法提取转子的瞬时频率，依据瞬时频率进行等角度采样，将角域信号进行频谱分析，提取油膜涡动故障频率。1小波时频分析小波分析是多尺度的时频分析方法。信号的时频特征时间、频率、幅值等可通过小波分解结果的时频参数来获得。若T为小波母函数，则T要满足条件RD。1J0式中，是T的傅立叶变换。通过对小波母函数进行平移和伸缩可得到小波基函数1F124A、A式中，U为尺度因子；B为平移因子。若原始信号为ST，则信号ST的连续小波时频变换为U，BI5T6TDT3式中，T是LF，T的复共轭。随着尺度因子U和平移因子B的变化，可得到小波系数A，B。利用灰度图表示U，6可得到小波时频分布灰收稿日期20110908作者简介胡爱军1971一男，副教授，主要从事状态监测与故障诊断研究。224汽轮机技术第54卷度图。2峰值搜索法阶比分析的关键是获得恒定角增量的等角度采样信号。基于瞬时频率的阶比分析依靠转子瞬时频率对时域信号进行等角度差值采样。文献4给出了一种转子瞬时频率估计的峰值搜索法，并指出当信号的一阶阶比成分可以和其它频率成分分开时，利用峰值搜索法可以获得较好的结果。利用峰值搜索法提取转子瞬时频率，首先对时域信号进行时频分析获得时频分布图，在时频分布图的一阶转速附近进行峰值搜索从而获得转子瞬时频率。峰值搜索法的具体算法为RPNL，K1A唱，MAXPN1，蜘“P4LP，KARGMAXPN，NIKI1一PKKILP式中，ARGMAX表示取最大值的参数在此为K；P为设定的搜索范围；、K为峰值坐标，即N时刻对应的瞬时频率为K。依据式4进行峰值搜索时，搜索的起始位置一般选择瞬时频率较为明显，且与其它频率成分相距较远地方进行峰值搜索。油膜涡动故障阶比分析步骤对转子升速过程油膜涡动信号进行阶比分析，具体步骤为1对转子升速过程油膜涡动时域信号进行小波时频分析，获得时域信号的时频分布图。2采用峰值搜索法对时频分布图进行峰值搜索，提取转子瞬时频率。3利用多项式对峰值搜索法提取的瞬时频率进行最小二乘拟合，得到拟合瞬时频率L厂T。4确定等角度采样间隔，并对最小二乘拟合获得的瞬时频率，T积分R2丌JTDT，以5R0式中，为等角度采样的键相时标；N为等角度采样时刻序列号；TO为等角度采样间隔。由式5即可获得等角度采样的键相时标。5依据等角度采样键相时标对时域信号进行插值采样，获得角域信号。其中，差值采样选用拉格朗日线性差值。6对角域信号进行频谱分析，识别油膜涡动故障频率。4实例分析在旋转机械故障模拟的BENTLY转子实验台上进行转子升速过程的油膜涡动实验，转子实验台的一阶临界转速在1900RMIN左右。实验时，通过转子实验台的转速控制装置调节转子转速，使转子转速由零逐渐升高，最高转速不超过转子一阶临界转速的2倍。在转子中部用电涡流传感器测量转子升速过程的转子的径向振动，其中采样频率1280HZ，截取转子升速过程中部分振动信号进行分析。转子升速过程的振动信号如图1所示，升速过程振动信号的幅值谱如图2所示。_6000趔7000名、一80000510时间，S图1转子升速过程振动信号图2升速过程振动信号幅值谱由图2可看出，升速过程振动信号的幅值谱出现了明显的频率模糊现象，不利于油膜涡动故障频率的识别。对转子升速过程振动信号进行阶比分析，转子升速过程振动信号的小波时频分布图如图3所示。苣槲图3转子升速过程振动信号的时频分布图由图3可看出，转子升速过程中，一阶转频最为突出，和其它阶频率成分没有交叉，区分明显，并且近似呈线性变化。在时频分布图中，以一阶转频的最右侧为起始位置进行峰值搜索，提取转子的瞬时频率。考虑到时频分布图中转子瞬时频率近似线性变化，采用一次多项式对瞬时频率进行最小二乘拟合。瞬时频率的峰值搜索结果及最小二乘拟合结果如图4所示。设定等角度采样间隔1RAD，对图4中的转子瞬时频率进行积分获取等角度采样的键相时标，并根据键相时标对第3期胡爱军等基于阶比分析的转子系统油膜涡动故障诊断研究225料骚时间，图4瞬时频率峰值搜索及最小二乘拟合结果图1所示的时域信号进行等角度差值采样，等角度采样结果如图5所示。对图5所示的等角度采样信号进行幅值谱分析即可得到阶比谱，其中等角度采样频率2R_2R。转子升速过程振动信号的阶比谱如图6所示。与转子升速过程振动信号的幅值谱相比，振动信号的阶比谱不存在“频率模糊”现象，各阶频率成分十分清晰。在阶比谱中，047倍转频处出现了较大的幅值，047倍转频恰为油膜涡动故障特征频率。由此可见，通过对转子升速过程振动信号的阶比分析，准确识别出了油膜涡动故障频率。5结论本文介绍了一种转子瞬时频率提取的方法峰值搜索法，并依据该方法对转子升速过程油膜涡动信号进行了阶比分析，结果表明，阶比分析可有效地避免传统频谱分析所产生“频率模糊”现象，能够准确识别出油膜涡动故障特征频率。阶比分析在转子升速过程油膜涡动故障诊断中具有一定的应用价值。2参考文献唐贵基，向玲，朱永利基于HHT的旋转机械油膜涡动和油膜振荡特征分析FJ中国电机工程学报，2O08，2827781杨志坚，丁康，杨茜基于频谱校正理论的阶比跟踪分析500055006000_6500坚一7000点趔罂750080000500100015002000角度TAD图5转子升速过程振动信号等角度采样结果阶比图6转子升速过程振动信号的阶比谱J机械工程学报，2009，451241453郭瑜，秦树人，梁玉前时频分析阶比跟踪技术J重庆大学学报，2002，25517204郭瑜，秦树人，汤宝平，等基于瞬时频率估计的旋转机械阶比跟踪J机械工程学报，2003，39332365冯志鹏，刘立，等基于小波时频框架分解方法的滚动轴承故障诊断J振动与冲击，2008，2721LO一1146邹剑，陈进，等转子裂纹仿真研究中的小波时频分析方法J应用力学学报，2002，19410一L37JZOU，JCHEN，YPPUWAVELETTTIMEFREQUENCYANALYSISOFTORSIONALVIBRATIONSINROTORSYSTEMWITHATRANSVERSECRACKJCOMPUTERSANDSTRUCTURES，2004，82711811187上接第222页1迎面风速在16MS到22MS之间变化时，压力变化不大，且当迎面风速在22MS时，排汽压力波动范围接近额定值，机组运行经济性最好。2冬季环境温度较低时排汽流量增大对压力变化影响很小，且压力低于最低运行压力。此时机组应采取防冻措施。夏季环境温度较高时，排汽压力较高，机组运行负荷增加，此时应采取喷水降温措施。3计算结果与实际偏差不大，说明参数优化结果较设计值好。可为同类1000MW空冷机组在变工况下选择合适的运行值提供参考依据。参考文献1丁尔谋发电厂空冷技术M北京水利电力出版社，1992严俊杰，张春雨直接空冷系统变工况特性的理论研究J热能动力工程，2000，6601603马义伟空冷器设计与应用M哈尔滨哈尔滨工业大学出版社，1998刘志昌工程流体力学M天津天津科学技术出版社，1995杨世铭传热学M北京高等教育出版社，1998周兰欣，杨靖600MW直接空冷