王文歡，王細(xì)洋，萬在紅

(南昌航空大學(xué)，南昌 330063)

0 引言

齒輪箱作為重要的傳動部件，被廣泛應(yīng)用于各類機(jī)械設(shè)備中，因此對其進(jìn)行故障診斷與分類的研究有著重要意義，是提高齒輪箱可靠性、減少停機(jī)檢修時間的重要一環(huán)。振動信號常常作為介質(zhì)來間接判斷齒輪的當(dāng)前運(yùn)行狀態(tài)，通過對比正常狀態(tài)和故障狀態(tài)振動信號的某些特征量，作為是否維修的決策依據(jù)，還可以判斷不同類型的故障，如不平衡、轉(zhuǎn)子碰磨、齒面磨損等。當(dāng)齒輪出現(xiàn)故障時，會在頻域上有所表現(xiàn)，出現(xiàn)各式各樣的調(diào)制邊頻帶，包含有豐富的齒輪故障信息，如果想對這部分頻帶進(jìn)行細(xì)致分析，提取出影響結(jié)果的細(xì)微差別，則需要足夠高的頻率分辨率。細(xì)化譜分析是一種提高某特定頻率段的分辨率的分析方法，能更精確的描述局部頻率，Hoyer［1］和謝明［2］等在算法的改進(jìn)和細(xì)化倍數(shù)的提高上作了很多的工作，使細(xì)化譜的性能更加優(yōu)良。在細(xì)化譜分析法的應(yīng)用上Hongzhong Ma［3］、段禮祥等［4］研究了細(xì)化譜分析法來診斷電動馬達(dá)和泵閥的故障，孫偉等［5］把細(xì)化譜分析法用在齒輪故障診斷中，他們的共同缺點(diǎn)是需要有豐富經(jīng)驗(yàn)的人員在較長時間做出決策，且缺乏對故障種類的判斷。

隱馬爾科夫模型具有很強(qiáng)的特征分類能力，已在語音識別等多個領(lǐng)域有成功的應(yīng)用［6］，它所具有的兩層隨機(jī)結(jié)構(gòu)與基于振動的故障診斷機(jī)理較為接近，通過可觀測的振動信號來估計(jì)設(shè)備的隱藏狀態(tài)。Baruah 等［7］詳細(xì)地討論了HMM 應(yīng)用在金屬切削刀具故障診斷和預(yù)診斷中一些的具體問題，并對算法作了修改，使之效率更高。Boutros等［8］研究了離散隱馬爾科夫模型對銑刀的磨損狀態(tài)的分類，同時也研究對軸承不同故障點(diǎn)的識別，兩者成功率均較高。Zhu Kunpeng 等［9］把HMM 用在高速精細(xì)機(jī)床刀具磨損的監(jiān)控中。鑒于軸承，各種刀具等信號與齒輪振動信號的類似性，將HMM 應(yīng)用到齒輪的故障診斷中。

本試驗(yàn)將細(xì)化譜分析方法與HMM 的特征分類方法相結(jié)合，引入到齒輪箱故障分類中。該方法首先利用時域同步平均法對信號預(yù)處理，然后變換到頻域，再對基頻及其倍頻處做細(xì)化譜處理，提取基頻、倍頻及其邊頻帶幅值作為特征向量，輸入到HMM 中，建立故障分類器，然后自動對故障類型作出識別。

1 HMM 建模

HMM 是一個統(tǒng)計(jì)學(xué)模型，是一個雙重隨機(jī)的過程:一重隨機(jī)為馬爾科夫鏈描述的狀態(tài)之間的轉(zhuǎn)移，是不可觀測的;二重隨機(jī)為狀態(tài)中觀測量的隨機(jī)分布，可觀測。HMM 包含有限個隱狀態(tài)，狀態(tài)根據(jù)其觀察量的概率分布產(chǎn)生一系列觀察量序列。一個HMM 可由如下參數(shù)描述［6］:

1)N:模型中馬爾科夫鏈狀態(tài)的個數(shù)。記N個狀態(tài)為θ1，θ2，…，θN，t 時刻馬爾科夫鏈的當(dāng)前狀態(tài)為q1，有qt∈(θ1，θ2，…，θN)。在實(shí)際應(yīng)用中常常有著具體的物理意義。

2)M:每個狀態(tài)中觀察值的數(shù)目。記M 個觀察值為v1，v2，…，vM，記t 時刻的觀察值為ot，ot∈(v1，v2，…，vM)。觀察變量對應(yīng)這模型的輸出，觀察序列是在第二層隨機(jī)的基礎(chǔ)上給出，由觀察變量組成。

3)π:初始狀態(tài)概率分布向量。π=(π1，π2，…，πN)，式中πi=P(qt=θi)。

4)A:狀態(tài)轉(zhuǎn)移概率矩陣。A=(aij)N×N，式中aij=P(qt+1=θj｜q=θi)，1≤i，j≤N。

5)B:觀察值概率矩陣。B=(bjk)N×M，式中bjk=P(ot=vk｜qt=θj)，1≤j≤N，1≤k≤M。

一個HMM 可以用以上5 個參數(shù)表示，記為λ=(N，M，π，A，B)，N 和M 在HMM 中對模型訓(xùn)練影響不大，一般選定后可將其作為一個常量，所以簡寫為λ=(π，A，B)。

實(shí)驗(yàn)中采集的數(shù)據(jù)是一系列離散的點(diǎn)，通常不會把這些點(diǎn)直接作為HMM 狀態(tài)的特征量，數(shù)據(jù)太多會導(dǎo)致主要特征淹沒其中，計(jì)算量會明顯增大;因此，需要選擇一個對故障敏感的特征量，以前的研究表明頻率成分的變化是故障敏感的，可以用細(xì)化譜分析特定頻帶提取特征量，進(jìn)行故障分類。

2 時域同步平均和細(xì)化譜算法

2.1 時域同步平均(TSA)

時域同步平均法是一種很有效的信號預(yù)處理方法，在齒輪的振動信號的分析中應(yīng)用廣泛，它可以從復(fù)雜的齒輪箱外部采集的振動信號中提取目標(biāo)齒輪振動信號。TSA 的原理是:通過與目標(biāo)齒輪的轉(zhuǎn)頻同步，多段信號累加，再平均而得到TSA信號，可以消減非目標(biāo)齒輪以外的信號和背景噪聲，提高信噪比，達(dá)到提取目標(biāo)齒輪振動信號的目的。

假設(shè)有一段多圈的采樣信號S(1，2，…)，它的TSA 信號表示為

其中，L 為平均的圈數(shù)，Ti為一圈采樣的點(diǎn)數(shù)［8］。因?yàn)辇X輪箱里不同齒輪轉(zhuǎn)一圈采樣的點(diǎn)數(shù)不同，也就是Tt不同，就可以對不同的齒輪作TSA 處理，提取出目標(biāo)齒輪的信號。

2.2 復(fù)解析帶通濾波器的細(xì)化譜理論

1)構(gòu)造FIR 非遞歸復(fù)解析帶通濾波器。

需要構(gòu)造的復(fù)解析帶通濾波器的幅頻特性如圖1 所示，其中f 為采樣頻率的一半，通帶寬為ω2-ω1，ωe為通帶的中心頻率，由于缺少負(fù)頻率部分，令正頻率通帶處的頻率幅值為2。為了構(gòu)造解析帶通濾波器，需先作一個實(shí)低通濾器，其沖擊響應(yīng)函數(shù)為

再對低通濾波器進(jìn)行復(fù)移頻，即可得到該解析濾波器，其沖擊響應(yīng)為

圖1 幅頻特性Fig.1 Amplitude-frequency characteristic

2)復(fù)解析帶通濾波器復(fù)調(diào)制細(xì)化譜分析方法［2］。

記原始信號的采樣頻率為fs，N 為傅里葉變換分析的點(diǎn)數(shù)，D 為細(xì)化倍數(shù)，M 為濾波器的半階數(shù)，采樣序列為x(n)(n=0，1，…，DN +2M)。則實(shí)現(xiàn)復(fù)解析帶通濾波的細(xì)化譜分析法的步驟是:首先構(gòu)造復(fù)解析帶通濾波器寬度為fs/D，再隔D點(diǎn)選抽一點(diǎn)、移頻，然后作N 點(diǎn)譜分析，最后把頻帶f1～f2用N 條獨(dú)立譜線來表示。

具體分解為以下5 個步驟:

①給定細(xì)化倍數(shù)與中心頻率。按實(shí)際需要給定細(xì)化頻段的細(xì)化倍數(shù)D 與中心頻率fe。

②構(gòu)造復(fù)解析帶通濾波器，濾波器的沖擊響應(yīng)為h0(n)，寬度為fs/D。

③選抽濾波。采用前一步構(gòu)造好的濾波器h0(n)對序列信號x(n)做選抽濾波，選抽出N點(diǎn)，選抽比為D。選抽的點(diǎn)數(shù)與FFT 運(yùn)算的點(diǎn)數(shù)都是N 點(diǎn)。用復(fù)解析帶通濾波器可以將濾波和選抽兩個過程結(jié)合起來計(jì)算，以提高效率。實(shí)信號經(jīng)過濾波器后，變成頻率在通帶以內(nèi)的復(fù)解析信號。記選抽后的復(fù)信號為g(m)=y(Dm)。

④復(fù)調(diào)制與移頻。對選抽后的復(fù)信號g(m)需要做復(fù)調(diào)制移頻，即將細(xì)化后的起始頻率移到零頻率點(diǎn)，得到g'(m)=g(m)e-jω1m，(m=0，1，…，N-1)。

⑤對g'(n)作N 點(diǎn)的FFT 與譜分析，不需要進(jìn)行頻率范圍的調(diào)整，就可以得到具有N 條獨(dú)立譜線的細(xì)化后的頻譜。

3 HMM 故障診斷試驗(yàn)

試驗(yàn)采用單級減速齒輪箱為實(shí)驗(yàn)對象，對4種狀態(tài)進(jìn)行了實(shí)驗(yàn)，即3 種齒輪常見故障狀態(tài)和1 種正常狀態(tài)。整個裝置由同步轉(zhuǎn)速為1 750 r/min、額定動率為4.5 kW 的交流電動機(jī)驅(qū)動，并在終端連接一電動機(jī)作為負(fù)載，其額定功率較電動機(jī)大(圖2)。電機(jī)的轉(zhuǎn)速由數(shù)字矢量驅(qū)動單元控制，其負(fù)載也是通過類似單元控制發(fā)電機(jī)電流來控制負(fù)載大小。采用基于虛擬儀器平臺搭建的數(shù)據(jù)采集系統(tǒng)。

對各種故障進(jìn)行模擬實(shí)驗(yàn)，主動輪齒數(shù)為21，從動論齒數(shù)為70。預(yù)加載的扭力為73 N·m。人為制造故障特征后，再放到實(shí)驗(yàn)臺上運(yùn)行，每隔10 min 作一次采集，采用編碼器控制采集，使每次開始采集的相位都一致，數(shù)據(jù)的采集寬度為10 s，采樣頻率為20 kHz。信號經(jīng)10 倍增益、A/D 轉(zhuǎn)換，再存儲。每個數(shù)據(jù)文件含20 萬個采樣點(diǎn)。實(shí)驗(yàn)臺集成加速度式振動傳感器、編碼器器、扭矩傳感器及速度傳感器等。根據(jù)參數(shù)計(jì)算得到嚙合頻率為612.5 Hz。

圖2 實(shí)驗(yàn)系統(tǒng)Fig.2 Test setup

3.1 數(shù)據(jù)的預(yù)處理和特征的提取

實(shí)驗(yàn)中的傳感器安裝在齒輪箱表面，齒輪徑向的加速傳感器采集的信號對齒輪故障最為敏感，故以此傳感器采集的數(shù)據(jù)為研究對象。采集得到的原始信號因含有噪聲且衰減嚴(yán)重，需要用時域同步平均技術(shù)來提高信噪比。圖3 為TSA預(yù)處理信號，從原始信號的頻域圖中看到頻譜比較雜亂，經(jīng)過時域同步平均后單個齒輪的嚙合頻率非常明顯，經(jīng)FFT 變換后求得嚙合頻率為612.7 Hz，與理論的頻率一致，并且出現(xiàn)二倍頻和三倍頻，對基頻和各級倍頻處做細(xì)化譜(Zoom-FFT)處理(圖4)，可以看出在主要頻率附近出現(xiàn)了不同大小的對稱或者不對稱的調(diào)制邊頻帶，根據(jù)之前的研究［4-5］，齒輪箱的狀態(tài)與基頻及各級倍頻及其邊頻帶有一定的關(guān)系。齒輪箱狀態(tài)由正常到故障，邊頻帶也會從無到有，最后逐漸增大，不同種類的故障會使頻譜出現(xiàn)不同的種類和大小調(diào)制邊頻帶。

圖3 TSA 預(yù)處理Fig.3 TSA preprocessing

圖4 信號的細(xì)化譜分析Fig.4 The ZOOMFFT spectrum of signal

3.2 HMM 建模

以試驗(yàn)采集的數(shù)據(jù)作為基礎(chǔ)，采用3.1 節(jié)介紹的數(shù)據(jù)預(yù)處理方法，從中提取特征向量來訓(xùn)練和驗(yàn)證HMM。試驗(yàn)采用一階齊次馬爾科夫模型描述隱狀態(tài)，假設(shè)特征量彼此相互獨(dú)立，每一個模型代表一種齒輪狀態(tài)。采用EM(Expectation Maximization)算法為每種故障狀態(tài)和正常狀態(tài)建立模型，通過Viterbi 算法計(jì)算特征向量在4 個HMM 下的輸出概率，以此判斷該特征向量所屬的模型，圖5 給出了處理過程模型。

3.3 試驗(yàn)分析

試驗(yàn)中共有4 組數(shù)據(jù)對應(yīng)正常、不平衡、裂紋和點(diǎn)蝕4 種狀態(tài)。每組采集的文件個數(shù)不一樣，分別采集122、112、172、180 個。每組數(shù)據(jù)的一半用于訓(xùn)練，另一半用于驗(yàn)證。對每個文件的數(shù)據(jù)先進(jìn)行TSA 去噪處理，后作FFT 變換，最后進(jìn)行細(xì)化譜分析提取基頻、倍頻及其邊頻帶幅值，歸一化后，統(tǒng)一量化處理，量化精度為40。

用交叉驗(yàn)證尋找最優(yōu)狀態(tài)數(shù)，當(dāng)HMM 的隱狀態(tài)數(shù)選為5 時誤差最小。設(shè)置收斂誤差為，最大迭代次數(shù)為50，4 種模型的迭代次數(shù)與每次迭代的對數(shù)似然概率見圖6，隨著迭代的進(jìn)行，4 種模型的對數(shù)似然概率也在不斷增加，在25 步內(nèi)達(dá)到收斂誤差，收斂快速。

圖5 特征提取和HMM 診斷模型Fig.5 Feature Extraction and HMM based on diagnostic model

圖6 HMM 訓(xùn)練圖Fig.6 Training graph of HMM

各種狀態(tài)的模型建立以后，用各組數(shù)據(jù)的另一半來診斷故障，驗(yàn)證模型的準(zhǔn)確性。4 組數(shù)據(jù)的另一半提取特征向量后分別輸入到建立好的4種模型中，模擬實(shí)驗(yàn)測試的統(tǒng)計(jì)結(jié)果見表1。從表中可知，在有限次的實(shí)驗(yàn)中，出現(xiàn)了極少數(shù)的不成功診斷，例如在正常的61 次實(shí)驗(yàn)中，有58 次作了正確的診斷，2 次誤診斷為裂紋，1 次誤診斷為點(diǎn)蝕，成功率為95.08%，分類結(jié)果較為理想。在其他故障的診斷中，也只出現(xiàn)了極少次的診斷錯誤，平均成功率達(dá)到了95.56%，在實(shí)際的診斷中通過多次采集數(shù)據(jù)和診斷，基本可以確定故障類型，診斷精度較高。

表1 故障診斷實(shí)驗(yàn)結(jié)果Table 1 Result of experiment

4 結(jié)論

1)以齒輪的振動信號作為研究對象，利用時域同步平均法準(zhǔn)確地抽取目標(biāo)齒輪的特征信號。

2)引入了基于細(xì)化譜的特征提取方法，對感興趣的頻帶進(jìn)行細(xì)化處理，能更準(zhǔn)確地反應(yīng)頻譜的微小變化。

3)應(yīng)用HMM 進(jìn)行了齒輪的故障種類識別，采用交叉驗(yàn)證獲取最優(yōu)狀態(tài)數(shù)，并通過齒輪箱試驗(yàn)系統(tǒng)驗(yàn)證了該方法的有效性。

［1］Hoyer E A，Stork R F.The zoom FFT using complex modulation［C］.Acoustics，Speech，and Signal Processing，IEEE International Conference on ICASSP '77，1977，2:78-81.

［2］謝明，丁康.基于復(fù)解析帶通濾波器的復(fù)調(diào)制細(xì)化譜分析的算法研究［J］.振動工程學(xué)報(bào)，2002，15(4):479-483.

［3］Ma H Z，Xu Q，Song J G，et al.The application of zoom FFT technique to the extraction of fault character of induction motor［C］.Condition Monitoring and Diagnosis，2008.CMD 2008.International Conference on，2008:221-225.

［4］段禮祥，張來斌，王朝暉，等.往復(fù)泵泵閥故障的細(xì)化譜診斷法［J］.儀器儀表學(xué)報(bào)，2004，25(4 增刊):568-570.

［5］孫偉，王細(xì)洋，徐英帥.基于復(fù)解析帶通濾波器的ZOOMFFT法應(yīng)用于齒輪故障診斷［J］.組合機(jī)床與自動化加工技術(shù)，2011(10):57-59.

［6］Rabiner L.A tutorial on hidden Markov models and selected application in speech recognition［J］.Proceedings of the IEEE，1989，77(2);257-286.

［7］Baruah P，Chinnam R B.HMMs for diagnostics and prognostics in machining processes［J］.International Journal of Production Research，2005，43(6):1275-1293.

［8］Boutros T，Liang M.Detection and diagnosis of bearing and cutting tool faults using hidden Markov models［J］.Mechanical Systems and Signal Processing，25(6):2102-2124.

［9］Zhu K P，Wong Y S，Hong G S.Multi-category micro-milling tool wear monitoring with continuous hidden Markov models［J］.Mechanical Systems and Signal Processing，2009，23(2):547-560.