基于改进峰值搜索法的旋转机械瞬时频率估计

本资料包含pdf文件1个，下载需要1积分

Instantaneous frequency estimation of a rotatingmachinery based on an improved peak search methodHU Ai-jun，ZHU Yu(Department of Mechanical Engineering，，Noah China Electric Power University，Baoding 071000，China)Abstract： An improved peak search method used to extract instantaneous frequencies of a rotor was proposed，andthis method was used for order tracking of a rotating machinery.With the improved peak search method，the differencebetween the first derivatives of the adjacent points in the instantaneous frequencies was considered as the criterion todetermine whether the searched peak was reasonable and to avoid the false peaks extracted with the traditional peak searchmethod，the instantaneous frequency estimation accuracy could be improved.Simulation results showed that the improvedpeak search method can reduce the efect of interference signals on the instantaneous frequency extraction，the proposedmethod is better than the traditional peak search method.With this method，the vibration signal of a rotor speeding up wasanalyzed，and good results were obtained。

Key words：improved；peak search；instantaneous frequency；order tracking旋转机械在工业生产中应用广泛，不可或缺。为确保旋转机械的安全稳定运行，对旋转机械进行状态监测与故障诊断尤为必要。目前，状态监测的主要手段之-为振动监测。旋转机械升、降速过程的振动信号为典型非平稳信号，即旋转机械在稳定运行状态下，其转速是波动的，振动信号也具有非平稳性，不适合用常规的频谱法分析∽比分析为-种新的非平稳信号处理方法，通过等角度采样将时间域非平稳信号转化为角域平稳信号，对该信号进行分析可提取振动信号的非平稳特征 J。传统的阶比分析方法分基于硬件的阶比跟踪与计算阶比跟踪。前者通过硬件实现转子的等角度采样，获得角域平稳信号；(而后者则对原始振动信号与转速脉冲信号同时进行等时间间隔采样，基金项目：国家自然科学基金资助项目(11072078)；中央高校基本科研业务费专项资金资助收稿 日期：2011-11-18 修改稿收到13期：2012-03-07第-作者 胡爱军 男 ，博士，副教授，1971年 1月生利用转速脉冲信号进行转速估计，由此获得等角度采样的发生时刻，并在此时对原始振动信号进行插值，实现等角度采样 j。基于硬件的阶比跟踪对硬件要求较高，而计算阶比跟踪也需用健相装置测量转速，两种方法均较繁琐。为降低阶比跟踪对硬件的要求，文献[1]提出基于瞬时频率估计的阶比跟踪新方法，通过对原始信号时频分布进行峰值搜索获得-阶转速对应的瞬时频率，并取得-定效果。利用峰值搜索法从信号时频分布中获得转子的瞬时频率会受到信号中干扰频率成分影响，降低瞬时频率估计精度，影响阶比分析结果 。为降低干扰信号对提取瞬时频率影响，文献[6]提出sT丌 Viterbi拟合法进行瞬时频率估计，将原始信号进行 sTFrI1变换获得时频分布，在时频分布图中先确定瞬时频率起始点，利用Viterbi方法搜索下-个频率点，并取得-定效果。该方法需将频率范围分成多段，计算每个频率段内峰值对应的频率点，为确保某时刻瞬时频率点也在其中，必然会增大频率等分段数的选择难度。

l14 振 动 与 冲 击 2013年第 32卷为此，本文提出改进的峰值搜索法，通过小波分析获得原始信号的小波时频分布图。对时频分布图进行峰值搜索时，将瞬时频率中相邻两点-阶导数差值作为搜索峰值是否合理的判别条件，避免传统峰值搜索法在干扰信号作用下提取虚假峰值，提高瞬时频率的估计精度。

1 小波分析 I9]小波分析为多尺度时频分析方法。信号的时频特征(时间、频率、幅值等)可通过小波分解结果的时频参数获得。若 (t)为小波母函数，则需满足条件：r d∞<∞ (1) Jo 09式中： (LO)为 (t)的傅里叶变换。

通过对小波母函数平移及伸缩可得小波基函数：/a，b㈩ ( ) (2)式中：n为尺度因子，b为平移因子。

若原始信号为 S(t)，则其连续小波时频变换为：(口，6)J s(t)02，6(t)dt (3)式中： (t)为 。。

由尺度因子 。及平移因子 6的变化，可得小波系数 (0，6)，用灰度图表示可得小波时频分布灰度图。

2 改进的峰值搜索法基于瞬时频率估计的阶比分析中，需利用原始信号估计参考轴转速。文献[1]利用峰值搜索法获得-阶转速，并计算原始信号的时频分布图，对时频分布图- 阶转速附近进行峰值搜索，将峰值时刻对应的频率视为该时刻转子的瞬时频率，具体算法为：P(n1， 1)arg，m，ax，P(n1， ) 1。-p 。 p (4)P(n ， )arg max P(n， )Jn ： n.± lki±1-p< < I±1 P式中：arg max表示取最大值参数(此为 )；p为设定的搜索范围；(n ， )为峰值坐标，即 时刻对应的瞬时频率为 。

当旋转机械振动信号中含幅值较大的干扰信号且该信号频率与-阶转速频率较接近时，利用式(4)搜索到的峰值将成为干扰信号峰值，其对应频率不再是-阶转速对应的频率，搜索到的-阶转速不再可靠。

旋转机械升速或降速过程中，在短时间段内参考轴作匀角加速运动 加 。据此知，旋转机械转速在任- 短时间段内变化是平缓的，对应的时间-频率(转速)曲线-阶可导。在等时间间隔采样时，相邻两采样点的时间间隔较小，两个采样点对应的-阶导数近似相等。因此，本文提出改进的峰值搜索法，将相邻两点- 阶导数差值作为搜索峰值是否合理的判别条件，避免传统峰值搜索法在干扰信号作用下提取虚假峰值。

具体步骤为：(1)选取信号时频分布图中-阶转速突出、与其它频率成分相距较远的区域据式(4)进行峰值搜索，确定相邻两峰值坐标( ， 。)、(n ， )。

(2)计算(n。， 。)、(n ， 。)点处-阶导数 d。dl - o - - 0n 1 - n0(3)设定 n- 与 n：时刻搜索范围分别为( -P，。p)与( -P， P)。为确保瞬时频率点在搜索范围内，P取值应适当大。计算 n- 时刻对应搜索范围内所有幅值，按幅值由大到小顺序，将该时刻幅值对应频率依次存人数组 A；同样，计算 n 时刻对应搜索范围内所有幅值，按幅值由大到小顺序，将该时刻幅值对应频率依次存人数组 曰。其中，A， 中第-个元素为最大幅值对应的频率。

(4)计算 d- ： 二 ，d ： 。确定/7，0 - T/， 1 //'2 - ， ld- 和 d：与 d 的差值分别为 - 和 ，其 中 。

(5)设定阈值S，分别比较 -，， 与s的大校若- ， ：均小于 s，认为此时A(1)与B(1)即为-阶转速在 凡- 和n：时刻对应频率；若 -。大于 s，则放弃A(1)，以A的下-点取代 (1)重复步骤(4)、(5)；当 - 小于s时， (i)即为n- 时刻对应频率。同理，当 大于S时，可得 n 时刻对应频率 日( )。

(6)用(n- ， - )、(n ， )取代点(n。， 。)、(凡 ，)重复步骤(2)～(5)可得(n- ， - )、(n，， ，)。依此重复可得各时刻-阶转速对应频率。

改进的峰值搜索法中，阈值 s的选取影响瞬时频率提取效果。为讨论阈值S的选取，设( ， )、(n ， 。)及(n ， )为时频图中瞬时频率相邻的3个点，有：s ：l l：I l (5) I- - 1 - j，，I，1- ，L0 I u I式中：Atn -n ：TI， -n。， 为时频分布图中相邻两频率点间隔，n0，1，2。

为保证搜索过程中不丢失正确的瞬时频率点，需确保s≥ 。由式(5)知，当信号瞬时频率随时间线性变化时， 0。但在时频分布图中，由于频率分辨率的限制，理论上瞬时频率为线性变化，也不能保证 0，即max>0。此 取l At、 、 。

考虑到瞬时频率变化的平缓性 ， 不能过大，即阈值s不能过大。仿真实验表明，当s取l j、I l、l I时，改进的峰值搜 效果较好。该法以时频曲第7期 胡爱军等：基于改进峰值搜索法的旋转机械瞬时频率估计 115线相邻两点-阶导数差值作为峰值是否合理的判别条件，可剔除因干扰信号作用产生的虚假峰值。

3 仿真信号分析为证明改进峰值搜索法优势，分别利用传统峰值搜索法与改进峰值搜索法对仿真信号进行分析。采样频率 1 280 Hz，采样时间3.2 s，仿真信号为：， )， 2 ㈤式中： ( )2sin(40，rt )，n(t)8sin(100t)。n(t)为在 1.5-2 s问的干扰信号。对仿真信号进行小波时频分析，如图 1所示。

图 1 仿真信号小波时频分布图Fig.1 Wavelet time-frequencydiagram of simulation signal图1中出现-条线性增大的瞬时频率曲线，与 (t)理论瞬时频率f40t-致。同时，在 1.5-2 S也可看到5O Hz的干扰信号 ( )，且频率与 频率接近。

分别用传统峰值搜索法与改进的峰值搜索法计算仿真信号中 (t)瞬时频率。传统方法搜索范围与改进方法搜索范围均包含了干扰信号。前者计算的瞬时频率见图2，后者计算的瞬时频率见图3。

采用-次多项式分别对图2、图3中曲线进行最小二乘拟合。传统方法提取瞬时频率. 38.86t-2.045，改进方法提取瞬时频率f/39.4-0.203。理论瞬时频率、传统方法提取瞬时频率及改进方法提取瞬时频率如图4所示。

图2 传统峰值搜索法瞬时频率Fig.2 Instantaneous frequency calculatedwith the traditional peak search method图4 瞬时频率的最小二乘拟合结果Fig.4 Least squares fiting resultsof Instantaneous frequency图5 最小二乘拟合结果局部放大图Fig.5 Partial enlargement figureof the least squares fitting results图3 改进峰值搜索法瞬时频率Fig.3 Instantaneous frequency calculatedwith the improved peak search method为便于观察，将图4局部放大见图 5。由图5看出，与传统峰值搜索法相比，改进峰值搜索法计算的瞬时频率与理论值更接近。为进-步证明改进方法的有效性，分别计算传统方法提取瞬时频率、改进方法提取瞬时频率与理论瞬时频率偏差： 厂 ------------- √ Ef(i)- ( )] (7)式中 i)为理论瞬时频率点 (i)为传统方法或改进方法瞬时频率点。经计算，传统方法提取的瞬时频率与理论瞬时频率偏差 4.009 Hz，改进方法提取的瞬时频率与理论瞬时频率的偏差 1.288 Hz。因此，改进峰值搜索法所取瞬时频率精度高于传统方法。

4 转子升速过程振动信号阶比分析在 Bently转子实验台上进行转子升速过程的油膜涡动实验。转子转速由零逐渐升高，用电涡流传感器测量升速中转子的径向振动，在靠近电机侧转子端部用光电传感器测量转速脉冲信号，采样频率. 1 280Hz。以转速 n2 900 r/min为起点，截取转子升速过程中-部分振动信号(图6)进行分析。升速过程振动信号的小波时时频分布见图7。

116 振 动 与 冲 击 2013年第 32卷s图6 转子升速过程振动信号Fig.6 Vibration signal of rotorspeeding up process图7 转子升速过程振动信号小波时频分布图Fig.7 Wavelet time-frequencydiagram of vibration signalduring rotor speeding up对转子升速过程的振动信号进行阶比分析。采用改进的峰值搜索法计算转子的瞬时频率，并利用光电传感器采集脉冲信号计算转子瞬时频率，两种方法所得结果如图8所示。由图8看出，改进峰值搜索法提取的-阶瞬时频率与利用脉冲信号计算结果基本吻合～改进峰值搜索法提取的-阶瞬时频率进行积分获取等角度采样健相时标，并据该时标对原始振动信号进行等角度差值采样，采样间隔0.1 rad。

等角度采样结果如图9所示。对采样结果进行 F盯变换，所得阶比谱如图 l0所示。由图10看出，-阶转频幅值最突出，且 0.47阶转频幅值较大。0.47阶转频成分的出现验证了转子升速过程中油膜涡动的发生。

- 8O 5 lO l5 20角度 ×102/rad图9 转子升速过程等角度采样结果Fig.9 Sampling results at constant angleincrements during rotor speeding up>O砉0- j粤X 0馨 0 0.47 . .1 。 。

阶比图10 转子升速过程振动信号的阶比谱Fig.10 Order spectrum of vibrationsignal during rotor speeding up图8 转子升速过程瞬时频率Fig.8 Instantaneous frequencyof rotor speeding up process采用改进峰值搜索法对转子升速过程振动信号进行阶比分析取得 良好效果，验证了改进峰值搜索法提取-阶转频的有效性。

5 结 论本文提出用于提取转子瞬时频率的改进峰值搜索法，并应用于旋转机械阶比跟踪。对仿真信号分析结果表明，改进峰值搜索法在提取瞬时频率时能降低干扰信号影响，提取瞬时频率效果优于传统峰值搜索法。

通过利用改进峰值搜索法对转子升速过程油膜涡动故障信号进行阶比分析，准确识别出油膜涡动的特征频率，验证了改进峰值搜索法的有效性。