楊通強,鄭海起,龔烈航,唐力偉

(1.解放軍理工大學 工程兵工程學院,江蘇 南京 210007;2.軍械工程學院一系,河北石 家莊 050003)

當前,階次分析技術(shù)是旋轉(zhuǎn)機械狀態(tài)監(jiān)測和故障診斷領(lǐng)域的1種基礎(chǔ)性分析方法,應用非常廣泛.由于可以避免轉(zhuǎn)速不穩(wěn)所帶來的“頻率模糊”現(xiàn)象[1],所以即使最終采用其他信號處理算法,也往往要基于其核心技術(shù)——角域采樣獲得數(shù)據(jù),而后進行分析.常用的階次分析方法有2種:硬件階次分析法和計算階次分析法.前者需通過硬件實現(xiàn)同步采樣,而后者則通過數(shù)值計算方法實現(xiàn)信號的角域重采樣,具有采樣率設置靈活、成本低、傳感器安裝方便等優(yōu)點,應用更為普遍,被近期絕大多數(shù)文獻[2～4]所采用.

與時域采樣類似,角域采樣中也不可避免地會遇到采樣率不足、采樣率過高及抗混濾波器設置等問題,本文依據(jù)采樣定理,分析了角域采樣率的設置原則,給出了計算公式;計算階次分析實驗中既涉及到時域采樣又有角域采樣.轉(zhuǎn)速變化的影響使得1種采樣相對另1種采樣必然為變采樣率的采樣,而跟蹤轉(zhuǎn)速變化的變采樣率抗混濾波又難以實現(xiàn),為此本文分析總結(jié)了只用1次定帶寬抗混濾波——時域抗混濾波,避開角域濾波的采樣率設置方法,并通過仿真實驗進行了驗證說明.

1 計算階次分析

計算階次分析法是近年來發(fā)展起來的1種先進技術(shù),其主旨在于通過信號處理算法將等時間間隔采樣信號轉(zhuǎn)化為等角度采樣信號,即同步采樣信號.與硬件階次分析法相比,此方法既可降低成本、放寬傳感器安裝限制又能提高角域采樣率,優(yōu)點顯著.具體實現(xiàn)過程如下:

首先要設定參考軸的角加速模式來確定重采樣的時間點.一般認為在1個小的時間段內(nèi)參考軸作勻角加速運動(特殊情況下加速模式已知的除外),前人[5]的研究證明了該假設的可行性.

在此前提下,參考軸的轉(zhuǎn)角θ可以表達為如下形式:

式中:b0,b1,b2為待定系數(shù);t為時間點.

對等時間間隔采樣信號,將3個依次到達的脈沖時間點t1,t2,t3和轉(zhuǎn)角增量Δφ帶入式(1),得

由式(2)求出bi(i=0,1,2)后,將其代入(1)式可得

根據(jù)式(3)便可求出等角度采樣時恒定角增量Δθ所對應的需要插值的時間點.再根據(jù)所求出的時間點,利用插值算法對振動信號進行插值,可以求出振動信號角域里對應于采樣時間點的幅值,再對其進行快速傅里葉變換(FFT),結(jié)果就是重采樣后的階次譜圖.

2 角域采樣率準則

由上節(jié)可見,角域重采樣技術(shù)是計算階次分析的關(guān)鍵步驟,若該環(huán)節(jié)出現(xiàn)誤差必將對后續(xù)分析產(chǎn)生較大影響,因此為了避免階次域混迭現(xiàn)象,本節(jié)重點研究階次采樣率問題.

由于階次分析與傳統(tǒng)頻率分析相比,基本數(shù)學變換理論依據(jù)均為傅里葉變換,只是分析對象由時域換為角域,所以下面類比時域采樣定理推導過程分析角域采樣準則:

設1個上限階次為Om的有限帶寬的連續(xù)角域信號 x(θ),其階次域表達式為X(O),且當｜O｜＞Om時,X(O)=0,如圖1a,b所示.若用1個脈沖序列δT(θ)去乘 x(θ)時,則其乘積 x(θ)δT(θ)是1個間距為 Tθ的脈沖序列,其在相應的瞬時具有與x(θ)值相等的強度,如圖1c,e所示,即取樣后的信號為

均勻脈沖序列δT(θ)的傅里葉變換δT(O)為

式中:δOs(O)為頻域內(nèi)以Os為周期的均勻脈沖序列.δT(O)也是一個均勻脈沖序列,如圖1d所示,每個脈沖的間隔為

根據(jù)卷積定理,采樣后信號 xs(θ)的傅立葉變換為 Xs(O),如圖1f所示,即

它是由X(O)與脈沖序列OsδOs(O)卷積而得,即原信號 x(θ)的頻譜每隔Os=2π/Tθ周期性地重復1次.顯然,只有滿足

時Xs(O)才能包含x(θ)的全部信息,周期出現(xiàn)的X(O)才不會產(chǎn)生首尾重疊現(xiàn)象,即不混迭.

綜上可見:1個在階次Om以上無階次分量的有限帶信號,可以由它在小于或等于π/Om均勻角度間隔Δθ上的取值唯一地加以確定.即

以上說明,如果在某一階次Om以上x(θ)傅里葉變換等于零時,則關(guān)于 x(θ)的全部信息均包含在它的采樣間隔小于π/Om的均勻采樣信號里.信號 x(θ)每隔Δθ角度被采樣1次,或者說以大于或等2Om采樣率進行采樣,這些采樣值 x(n)包含了 x(θ)在每一個θ角度值的信息.

3 角域重采樣實驗采樣率設置方法分析

圖1 角域采樣原理Fig.1 Effect in the order domain of sampling in the angular domain

設某轉(zhuǎn)速為ω的轉(zhuǎn)速機械的時域信號為X(ω,t),對于變轉(zhuǎn)速機械,ω也是時間t的函數(shù).對其進行時域采樣,采樣頻率為fs,則采樣脈沖間隔為

時域采樣信號為

式中:δT(t)為以T為周期的單位沖擊信號序列.根據(jù)采樣定理,抗混濾波器的上限頻率為

那么,在角域內(nèi)可將此信號看作變采樣率采樣信號,且角域采樣率隨ω的變化而變化:

相應地,角域內(nèi)原來的時域抗混濾波器的上限階次也成為ω的函數(shù),為

因此,欲對時域采樣信號進行采樣階次為定值O—s的重采樣時,首先要對該時域信號進行時域內(nèi)的跟蹤濾波.根據(jù)式(8)得O—m=O—s/2,代入式(14)可得上限頻率為

以上分析可見,計算階次分析實驗中既涉及到時域采樣又有角域采樣.轉(zhuǎn)速變化的影響使得1種采樣相對另1種采樣必然為變采樣率采樣,而跟蹤轉(zhuǎn)速變化的變采樣率抗混濾波又難以實現(xiàn),下面分析總結(jié)只用1次定帶寬抗混濾波——時域抗混濾波,避免跟蹤抗混濾波的采樣率設置方法.

(1)由式(15)及時域采樣定理,欲保留所需信息,則 fs需滿足

表明時域采樣率應根據(jù)最大轉(zhuǎn)速和預期采樣階次來設置,否則會采樣不足.

(2)由式(13)知,在角域重采樣時,若不進行抗混濾波,為了滿足采樣準則,則需Os≥max(2πfs/ω)即

4 算例分析

設某變轉(zhuǎn)速信號y為

式中:ω=20t,t=[1,3].

從信號表達式可見其中僅含有軸頻的15階信息.假設欲采用先時域采樣,再作角域重采樣和階次分析來分析該階信息,即Om=15,O—s≥2Om,取其最小值O—s=2Om代入式(17),確定時域采樣頻率為 fs≈287 Hz.現(xiàn)分別采用2種時域采樣率對其進行采樣,結(jié)果如圖2和圖3.圖2為采樣率滿足該條件時信號采樣結(jié)果,圖中信號表現(xiàn)為規(guī)則的掃頻信號,準確地表達了原信號信息;圖3為根據(jù)平均轉(zhuǎn)速計算所得采樣率245 Hz時的采樣結(jié)果,顯然信號由于采樣率不足而產(chǎn)生了嚴重失真.說明了式(17)的正確性.

圖2 采樣率大于287 Hz時的采樣信號Fig.2 Sampled signal when the sampling frequency is above 287 Hz

圖3 采樣率為245 Hz時的采樣信號Fig.3 Sampled signal when the sampling frequency is 245 Hz

下面對式(10)進行驗證:

設上算例信號中含有40階高階信號,則信號表達式變?yōu)?/p>

仍以分析15階信息為目的,由式(18)得:min(O—s)≈90階.

圖4 正常重采樣時的階次譜Fig.4 Order spectrum of the normal sampled signal

圖5 采樣率過低時的階次譜Fig.5 Order spectrum of the undersampled signal

5 結(jié)論

計算階次分析技術(shù)作為旋轉(zhuǎn)機械狀態(tài)監(jiān)測和故障診斷領(lǐng)域的1種基礎(chǔ)性分析方法,應用日益廣泛.角域重采樣是計算階次分析的關(guān)鍵步驟,其采樣精度將直接影響分析結(jié)果.角域采樣率的設置需滿足Os≥2Om;計算階次分析實驗中,轉(zhuǎn)速變化的影響使得1種采樣相對另1種采樣必然為變采樣率采樣.本文分析總結(jié)了只用1次定帶寬抗混濾波——時域抗混濾波,避開變帶寬跟蹤濾波的采樣率設置方法:時域采樣率應根據(jù)最大轉(zhuǎn)速和預期采樣階次來設置,fs≥max(ω/2π,否則采樣不足;角域重采樣率則應根據(jù)最小轉(zhuǎn)速和時域采樣率來確定,Os≥2πfs/min(ω),否則會產(chǎn)生階次混迭.

[1]朱繼梅.非穩(wěn)態(tài)振動信號分析[J].振動與沖擊,2000,19(1):86-87.ZHU Jimei.Nonstationary vibration signal analy sis[J].Journal of Vibration and Shock.2000,19(1):86-87.

[2]郭瑜,秦樹人.旋轉(zhuǎn)機械非穩(wěn)定信號的偽轉(zhuǎn)速跟蹤階比分析[J].振動與沖擊,2004,23(1):61-64.G UO Yu,QING Shuren.Pseudo-speed tracking order analysis for non-stationary vibration signal of rotating machinery[J].Journal of Vibration and Shock.2004,23(1):61-64

[3]楊通強,鄭海起,唐力偉.基于角域平均和倒階次譜分析的齒輪箱故障診斷[J].機械科學與技術(shù),2005,25(4):452-454.YANG T ong qiang,ZHENG Haiqi,TA NG Liwei.Gearbox fault diagnosis based on angle domain averaging and order cepstrum analysis[J].Mechanical Science and Technology.2006,25(4):452-454

[4]李輝 鄭海起唐力偉.瞬時頻率估計的齒輪箱升降速信號階次跟蹤[J].振動、測試與診斷,2007,27(2):125-128.LI Hui,ZHENG Haiqi,TANG Liwei.Order tracking of speed-up signal for gearbox based on instantaneous frequency estimation[J].Journal of Vibration,Measurement&Diagnosis.2007,27(2):125-128.

[5]F YFE K R,MUNCK E D S.A naly sis of computed order tracking[J].Mechanical Sy stems and Signal Processing,1997,11(2):187-205.