变转速下f细化阶次全息谱分析

变转速下f细化阶次全息谱分析

为了监控旋转机械故障，实现多通道多信息融合的方法主要有三种：全谱理论、全矢谱理论和全谱理论。而全息谱通过对相互垂直两方向振动信号各倍频幅值、相位信息融合以椭圆形式直观表达转子的运动轨迹,能较好对转子横向裂纹、不平衡、弯曲、动静部件径向摩擦、油膜涡动、轴承座松动等故障进行识别。转速稳定情况,对互相垂直的两方向振动信号进行FFT分析可得转速频率的倍频分量幅值及相位信息。但在极端环境及特殊工况下,旋转机械转速极不稳定,所得振动信号非平稳,如风电机组在正常运行时转速不稳定,不能通过FFT分析直接提取各倍频分量。为解决变转速下全息谱分析问题,魏玉果等将阶次跟踪引入全息谱分析,通过阶次重采样获得各阶次分量幅值及相位信息,再进行全息谱分析,称阶次全息谱。与传统全息谱分析因各倍频幅值、相位存在误差需校正相同,阶次全息谱分析在阶次谱计算中也存在误差,同样需校正。转速不稳定时由于进行阶次重采样,阶次分量会很明确,阶次全息谱只需计算部分阶次;传统的全息谱校正方法不适合阶次全息谱校正,而FT细化校正法适合对局部进行细化校正,且在密集阶次谱情况下也能获得较高精度,较适合阶次全息谱校正,由此提出FT细化校正阶次全息谱分析方法。1基于三对插值的阶次全谱1.1采用插值滤波器重采样,获得所采样点阶次跟踪阶次跟踪技术通过振动信号、转速信号结合实现振动信号等角度采样。对阶次跟踪技术的研究主要分两大类:①硬件式阶次跟踪,直接通过模拟设备实现对振动信号等角度采样;②计算式阶次跟踪(ComputedOrderTrackingCOT)。硬件式阶次跟踪设备较昂贵且安装不便,故本文采用计算阶次跟踪技术中较精确的三次样条插值计算阶次跟踪。即用三次样条插值获得等角度采样点序列号,再据采集的X,Y方向振动信号序列用插值滤波重采样完成阶次跟踪,具体步骤如下:(1)用同一采样频率fs同时采样获得互相垂直的振动信号序列X(n),Y(n)及转速脉冲信号序列S(n)。(2)对转速脉冲信号进行三次样条插值,据所选角域采样频率Os(即等角度采样时每转采样点数)及每转所采脉冲数求出所需等角度重采样序列号T(n)。(3)据所得等角度重采样序列号T(n),用插值滤波器对原始采样数据插值以获得信号等角度采样式中:i=0,1,2,3…n为重采样后序列号i,k的取值在T(i)附近点;hs(t)为内插函数。由此获得X,Y方向角域平稳信号。(4)求出X,Y两方向阶次谱,获得各阶次幅值及相位fx,Ax,ax,fy,Ay,ay。将两方向振动信号各阶次写成以t为参数的方程为(5)消去式(4)、(5)参数t,可得各倍频全息谱方程为进而可计算出全息谱中重要参数,如长短轴、离心率、正余弦系数等。1.2旋转机械转速变化用一组在二倍频含密集频谱的非平稳模拟信号X(t),Y(t)进行阶次全息谱分析,表达式为方便模拟,取w=at2/2,旋转机械工频即转速频率为at,本文用加速度a=100r/s2模拟旋转机械转速变化。取采样频率1024Hz、采集1024点见图1。由于全息谱为对各倍频分析,而对图1中转速变化信号不能直接提取各倍频。因此对其以60的角域重采样频率进行基于三次样条插值的阶次重采样,截取前1024点获得重采样角域稳定信号见图2。对重采样后X、Y方向信号进行阶次全息谱分析见图3。2nt详细校正阶的全集2.1阶次全息谱图的计算步骤FFT+FT法为在不增加采样点数情况下,先对信号进行FFT求得全景谱,再针对所需细化局部用连续FT细化校正计算获得局部细化、精度极高的频谱。本文提出FT细化校正阶次全息谱分析方法,先通过三次样条插值阶次跟踪获得等角度重采样的角域信号,再通过FT细化校正获得各阶次的精确幅值及相位信息,据幅值、相位值得出阶次全息谱图,具体计算流程为(1)基于三次样条插值的阶次跟踪重采样后获得角域采样频率Os阶,角域工频为1阶,角域平稳信号为x(n),y(n)。(2)对x(n),y(n)进行N点FFT分析,找出0~1阶间幅值最大的对应阶次作为分倍频分析阶次。计算出各倍频对应的阶次谱序列号k,按式(9)、(10)分别计算各倍频FT细化校正的起止频率f1,f2:(3)针对实际信号确定细化倍数D,一般取20~50之间。据f1,f2,D计算出阶次间隔Δf:(4)将{f1,f1+Δf,f1+2Δf,…,f1+lΔf=F2}分别代入式(12)、(13),即可获得细化校正后阶次谱实、虚部值,取该阶次成分中幅值最大的作为该倍频幅值,即可分别求出X,Y方向各倍频校正后的幅值及相位Ax,ax,Ay,ay:获得精确幅值、相位后,即可获得精确阶次全息谱图。FT细化校正阶次全息谱分析计算流程见图4。2.2全息谱图校正对式(7)、(8)模拟信号进行FT连续细化校正阶次全息谱分析,得阶次全息谱图见图5。对式(7)、(8)模拟信号进行阶次跟踪后基于比值校正法进行阶次全息谱分析,所得阶次全息谱见图6。由图3与图5、图6对比看出,不校正直接计算所得各倍频全息谱图区别较大。由表1看出,直接计算误差较大,因此计算全息阶次谱校正较重要。由图5、图6对比看出,在1X,3X,4X处基本相同,但在含密集谱的2X处,图5椭圆更扁平些且初始点位置不同,即两幅图离心率、初始相角不同。因密集频谱影响,直接比值校正法所得幅值、相位均有一定误差,反映到全息谱参数上即为正余弦系数与离心率不同。由表1结果分析知,在2X处直接比值校正所得X,Y向正余弦参数结果与理论值绝对差分别为0.0974,0.1001,0.0835,0.0467,离心率绝对差为0.0047,而FT细化校正法所得X,Y向正余弦参数结果与理论值绝对差分别为0.0080,0.0006,0.0080,0.0015,离心率结果基本相同。故FT细化校正法较比值校正法精度高,所得全息谱亦更精确。3振动信号采集为验证本文方法的可行性与正确性,对转子实验台进行变速实验测试见图7,共采集4路信号,第1路为转速脉冲信号,第2路为X向振动信号,第3路为Y向振动信号,第4路为实验台振动加速度信号。用512Hz采样率进行同步采样,采集两组数据分别为正常、加偏心质量。加偏心质量测得的转速信号转换为标准脉冲信号见图8,加偏心质量测得的X,Y方向振动信号见图9。用采到的脉冲信号及X,Y向振动信号,据本文方法进行重采样,所得加偏心质量重采样后X,Y向信号见图10。4测定转子故障的光谱分析本文针对极端环境、特殊工况只能采集到转速不稳定振动信号、阶次全息谱分析精度不高问题,提出FT细化校正阶次全息谱分析,并对模拟信号及转子试验台采集信号分析,结论如下:(1)对非平稳信号进行阶次全息谱分析时,可用三次样条插值方法进行阶次重采样,FT细化校正法进行频谱计算,能在存在密集普情况下获得较精确的全息谱。(2)对转子试验台信号进行阶次全息谱分析与直接阶次谱分析对比充分说明该方法的有效性。通过与正常信号分析对比验证FT细化校正阶次全息谱分析方法能正确判断转子试验台故障。对原X,Y方向信号及重采样后X,Y方向信号进行频谱与阶次谱分析,其局部见图11、图12。由图11看出,转速频率在30~40Hz波动,此时则不能选择工频,但重采样后即可较好提取到1阶频率。对基于三次样条插值跟踪重采样信号直接进行阶次全息谱分析所得阶次全息谱见图13。用FT细化校正方法对重采样数据进行阶次全息谱分析,见图14。由图14看出,工频振幅较大,且为椭圆,而转子失衡时全息谱特征除工频椭圆外基本无其它倍频分量,工频椭圆形状较扁,可