基于級聯(lián)雙穩(wěn)隨機共振和多重分形的機械故障診斷方法研究

郝 研，王太勇，，萬 劍，張 攀

（1.天津大學(xué) 精密儀器與光電子工程學(xué)院，天津 300072；2.天津大學(xué) 機械工程學(xué)院，天津 300072）

分形的本意為不規(guī)則的、分?jǐn)?shù)的、支離破碎的物體，可以理解為局部和整體在某個方面有相似性［1-2］。在自相似分形中，一個分形物體是由單一分形維數(shù)來表征的，但對于大多數(shù)物理現(xiàn)象，單一分形的描述過于簡單化。一些復(fù)雜系統(tǒng)的非線性特征取決于某些物理量的空間分布，多重分形理論由此發(fā)展起來。在分形理論中，分形維數(shù)是一個非常重要的參數(shù)，可以定量地刻畫混沌吸引子的“奇異”程度，在非線性行為的定量描述中得到了廣泛應(yīng)用［3-5］。廣義維數(shù)方法是一種描述多重分形的方法，在多重分形的研究中應(yīng)用的較多［6-8］。眾所周知，噪聲會影響分形維數(shù)的大小，使其無法準(zhǔn)確的反映系統(tǒng)的分形特征。為了得到系統(tǒng)的真實分形維數(shù)，就必須對原始信號進行濾波降噪處理。

隨機共振的概念是Benzi等［9］在研究古氣象冰川問題時提出，后來在信號處理和檢測方面得到了廣泛的發(fā)展［10-12］。在信號降噪方面，隨機共振方法有其獨特之處，利用信號和噪聲在非線性系統(tǒng)中的協(xié)同作用，使其產(chǎn)生共振，將高頻噪聲能量向低頻轉(zhuǎn)移，達(dá)到信號降噪和整形的目的。另外，通過對雙穩(wěn)系統(tǒng)的級聯(lián)，可以進行多次共振降噪，不斷將噪聲能量轉(zhuǎn)移，因此該方法的降噪效果非常明顯。

在機械故障診斷方面，不同的故障狀態(tài)下，非線性因素對機械振動信號的影響是不同的，如果能在一定程度上辨識非線性因素，對于提取機械系統(tǒng)的故障特征和識別系統(tǒng)的故障類型是很有幫助的［13-15］。結(jié)合級聯(lián)雙穩(wěn)隨機共振的降噪性能和廣義維數(shù)的機械故障診斷特點，提出了基于級聯(lián)雙穩(wěn)隨機共振和多重分形的故障診斷方法，即先對含噪信號進行隨機共振降噪，再進行廣義維數(shù)計算分析，從而達(dá)到故障診斷的目的。

1 廣義維數(shù)及相關(guān)性判斷方法

1.1 廣義維數(shù)

多重分形是定義在分形集合上的由多個標(biāo)度指數(shù)的奇異測度所組成的集合，刻畫了分形測度在支集上的分布狀況，可以表示為具有不同維數(shù)的分形子集之并。廣義維數(shù)方法是描述多重分形的一種有效方法，廣義維數(shù)的計算一般采用覆蓋法完成。

用覆蓋法計算廣義維數(shù)的關(guān)鍵是定義覆蓋集合的概率，一般使用頻率來近似概率，可得：

覆蓋法使用尺度為ε相同大小的盒子對整個集合進行覆蓋，需要盒子的總數(shù)為N，設(shè)點落于第i個盒子的概率為pi（ε），給定參數(shù)q，可以計算General信息熵：

從而得到廣義維數(shù)的定義為：

通過改變q的值可得到不同標(biāo)度的子集，由此可知，當(dāng)q=0、q=1和q=2時，分別對應(yīng)盒維數(shù)、信息維數(shù)和關(guān)聯(lián)維數(shù)。廣義維數(shù)的具體計算步驟參見文獻［8］。

1.2 廣義維數(shù)相關(guān)性判斷方法

文獻［8］提出了廣義維數(shù)的相關(guān)性判斷法，該方法針對一維廣義維數(shù)序列，計算待檢測信號廣義維數(shù)序列和各種狀態(tài)下樣本信號廣義維數(shù)序列的相關(guān)系數(shù)，用于判斷故障。

設(shè)共有j種狀態(tài)，每種狀態(tài)對應(yīng)一個廣義維數(shù)序列，q=0，1，…，n-1。待檢測信號為x，對應(yīng)的廣義維數(shù)序列為。則待檢測信號與第j種狀態(tài)的相關(guān)系數(shù)為：

此函數(shù)反映了待檢測信號x和第j種狀態(tài)的關(guān)聯(lián)程度。R（x，j）越大，則待檢測信號x和第j種狀態(tài)的相關(guān)程度就越強，反之，R（x，j）越小，則待檢測信號x和第j種狀態(tài)就越相離。

根據(jù)這種方法，可以先計算得到各種故障信號的廣義維數(shù)序列，將這些廣義維數(shù)序列作為故障樣本。如果要檢驗?zāi)承盘柕墓收希瑒t可以先求得其廣義維數(shù)序列與各故障樣本序列的相關(guān)系數(shù)，然后根據(jù)相關(guān)系數(shù)的大小，判別該信號屬于何種故障類型。

2 級聯(lián)雙穩(wěn)隨機共振

2.1 隨機共振基本原理

在一些非線性系統(tǒng)的作用下，輸入的信號不變，而改變輸入噪聲的強度，輸出信噪比會達(dá)到一個極大值，這樣種現(xiàn)象被稱為隨機共振。

隨機共振最大的特點為：噪聲的引入不僅不影響非線性系統(tǒng)的響應(yīng)質(zhì)量，反而可以加強系統(tǒng)的輸出，提高系統(tǒng)輸出的信噪比。一般的降噪方法都是以抑制或消除噪聲為主，而隨機共振的方法則是以利用噪聲為主。

圖1 雙穩(wěn)隨機共振系統(tǒng)Fig.1 Bistable stochastic resonance system

在噪聲強度D=0時，系統(tǒng)存在臨界值λ。λ可通過雙穩(wěn)勢函數(shù)極點與拐點重合的條件求得。

2.2 級聯(lián)雙穩(wěn)隨機共振

級聯(lián)雙穩(wěn)隨機共振（cascaded bistable stochastic resonance，CBSR）系統(tǒng)就是將若干個雙穩(wěn)系統(tǒng)串聯(lián)相接，圖2是級聯(lián)雙穩(wěn)隨機共振系統(tǒng)的結(jié)構(gòu)圖，其中x1（t）、x2（t）、…等分別是第一級雙穩(wěn)系統(tǒng)U1（x）、第二級雙穩(wěn)系統(tǒng)U2（x）、…等的輸出信號。雙穩(wěn)隨機共振系統(tǒng)可以將高頻能量轉(zhuǎn)移到低頻，使低頻信號能量增強，高頻信號能量減小，而雙穩(wěn)系統(tǒng)的級聯(lián)可以將高頻能量不斷地轉(zhuǎn)移，高頻成分逐漸被慮除，達(dá)到良好的降噪和整形效果。

圖2 級聯(lián)雙穩(wěn)隨機共振系統(tǒng)Fig.2 Cascaded bistable stochastic resonance system

2.3 級聯(lián)雙穩(wěn)隨機共振消噪濾波實驗

建立一個雙穩(wěn)隨機共振系統(tǒng)，對應(yīng)公式（6）中各參數(shù)分別為：a=0.1，b=1，A=1，噪聲強度D=0.6，信號頻率f0=10 Hz，采樣頻率fs=1 000 Hz，采樣點數(shù)為2 048點，二次采樣頻率fsr=8 Hz。含噪正弦信號的波形和頻譜如圖3（a）和圖3（b）所示，圖3（c）和圖3（d）分別為第一級級聯(lián)雙穩(wěn)隨機共振系統(tǒng)的輸出信號和頻譜，圖3（e）和圖3（f）分別為第二級級聯(lián)雙穩(wěn)隨機共振系統(tǒng)的輸出信號和頻譜。從圖中可以看出，高頻噪聲隨著級聯(lián)次數(shù)的增加而減小，時域波形也越來越光滑，這說明級聯(lián)雙穩(wěn)隨機共振系統(tǒng)具有良好的降噪性能。

圖3 含噪正弦信號級聯(lián)雙穩(wěn)隨機共振輸出Fig.3 Output of noisy sinusoidal signals of CBSR

3 實例分析

對6205-2RS軸承進行實驗，在實驗中，軸承有三種狀態(tài)：正常狀態(tài)、內(nèi)圈故障狀態(tài)和外圈故障狀態(tài)。軸承在1 797r/min的額定轉(zhuǎn)速下工作，利用振動加速度傳感器測取軸承的振動信號，在采樣頻率為12 kHz、采樣點數(shù)為20 000點的條件下，得到軸承采樣信號。

建立級聯(lián)雙穩(wěn)隨機共振系統(tǒng)，系統(tǒng)參數(shù)為：a=0.001，b=1，二次采樣頻率fsr=4Hz。對軸承的三種狀態(tài)分別進行級聯(lián)雙穩(wěn)隨機共振降噪處理，圖4、圖5和圖6分別為正常軸承信號、內(nèi)圈故障信號和外圈故障信號經(jīng)過級聯(lián)雙穩(wěn)隨機共振系統(tǒng)的輸出。

軸承在1 797 r/min的額定轉(zhuǎn)速下工作，工作頻率大概在30 Hz左右。對于三種不同的狀態(tài)，信號故障特征主要集中在1 000 Hz以下的低頻段，而超過1 000 Hz的高頻段信號一般作為噪聲進行處理。

采用“sym5”小波基函數(shù)，以極大極小閾值原則分別對三種狀態(tài)的振動信號進行自適應(yīng)濾波處理，所得結(jié)果如圖7所示。

對比濾波后信號的頻譜，在正常狀態(tài)下，雖然小波濾波方法和CBSR濾波方法都可以將高頻噪聲信號慮除，但是CBSR濾波方法加強了低頻信號的能量，使低頻信號的頻譜特點更加鮮明的表現(xiàn)出來。在內(nèi)圈和外圈故障狀態(tài)下，小波濾波方法對高頻信號的濾波效果明顯不如CBSR濾波方法，高頻噪聲信號的能量占有主導(dǎo)地位，無法凸現(xiàn)信號的低頻特性。而CBSR濾波方法不但消除了高頻噪聲，而且加強了低頻信號的能量。對比可知，CBSR濾波方法更適用于消除高頻噪聲，增強低頻信號能量，提高信噪比，有利于凸現(xiàn)故障特點，便于確定故障狀態(tài)。

對于在不同狀態(tài)下得到的一組采樣數(shù)據(jù)，分別對原始數(shù)據(jù)、經(jīng)過小波濾波后的數(shù)據(jù)和經(jīng)過CBSR濾波后的數(shù)據(jù)求取廣義維數(shù)，并將其中的盒維數(shù)、信息維數(shù)和關(guān)聯(lián)維數(shù)提取出來，可以得到表1、表2和表3。這三個經(jīng)過不同濾波方法處理的廣義維數(shù)序列就作為廣義維數(shù)故障樣本，用于相關(guān)系數(shù)計算，以便于判斷故障。

圖4 正常軸承信號級聯(lián)雙穩(wěn)隨機共振輸出Fig.4 Output of normal bearing signals of CBSR

圖5 內(nèi)圈故障信號級聯(lián)雙穩(wěn)隨機共振輸出Fig.5 Output of inner raceway fault signals of CBSR

圖6 外圈故障信號級聯(lián)雙穩(wěn)隨機共振輸出Fig.6 Output of outer raceway fault signals of CBSR

圖7 三種狀態(tài)振動信號的小波濾波輸出Fig.7 Outputs of vibration signals of CBSR in three conditions

表1 各狀態(tài)下原始數(shù)據(jù)的廣義維數(shù)序列Tab.1 The general dimension sequences of signals in different status

表2 小波濾波后各狀態(tài)的廣義維數(shù)序列Tab.2 The general dimension sequences of signals in different status filtered by wavelet

表3 級聯(lián)雙穩(wěn)隨機共振濾波后各狀態(tài)的廣義維數(shù)序列Tab.1 The general dimension sequences of signals in different status filtered by CBSR

任意抽取一個軸承，對其進行測試。采集該軸承的一組振動信號，對其原始數(shù)據(jù)、經(jīng)過小波濾波后的數(shù)據(jù)和經(jīng)過CBSR濾波后的數(shù)據(jù)分別計算廣義維數(shù)，并將其中的盒維數(shù)、信息維數(shù)和關(guān)聯(lián)維數(shù)提取出來，結(jié)果如表4所示。

表4 振動信號的廣義維數(shù)序列Tab.4 The general dimension sequences of signals in unknown status

表5 相關(guān)系數(shù)計算結(jié)果Tab.3 The results of correlation coefficient

使用公式（4），分別計算這三組數(shù)據(jù)的廣義維數(shù)序列與相對應(yīng)的廣義維數(shù)故障樣本之間的相關(guān)系數(shù)，結(jié)果如表5所示。

由表5可知，對原始數(shù)據(jù)進行處理得到的相關(guān)系數(shù)比較接近，無法判斷故障；經(jīng)過小波濾波后，得到的相關(guān)系數(shù)的大小有一定的差距，由于待檢測軸承的廣義維數(shù)序列與外圈故障樣本序列的相關(guān)系數(shù)較大，可以判斷待檢測軸承應(yīng)該處于外圈故障狀態(tài)；經(jīng)過CBSR濾波后，得到的相關(guān)系數(shù)的大小差距更加明顯，由于待檢測軸承的廣義維數(shù)序列與外圈故障樣本序列的相關(guān)系數(shù)最大，且遠(yuǎn)大于其它結(jié)果，可以判斷待檢測軸承的故障狀態(tài)與外圈故障狀態(tài)樣本一致，即該軸承的故障為外圈故障。經(jīng)實際驗證確認(rèn)，該軸承確實處于外圈故障狀態(tài)。

基于上述分析可知，不使用濾波方法，直接對原始數(shù)據(jù)進行分析處理是無法正確識別故障的，這是由于大量噪聲的影響，使得廣義維數(shù)無法正確反映信號的非線性特征，以致無法識別軸承狀態(tài)。雖然使用小波濾波方法和CBSR濾波方法都可以正確的識別軸承狀態(tài)，但是，對比相關(guān)系數(shù)可知，經(jīng)過CBSR濾波后，相關(guān)系數(shù)的大小差距較大，可以更容易的判斷軸承狀態(tài)。雖然小波濾波方法消除了大量的高頻噪聲，但是，原有的低頻信號能量仍然較弱，得到的廣義維數(shù)雖可以反映其非線性特征，但并不十分明顯；CBSR濾波方法不但可以有效地消除高頻噪聲，還能夠利用噪聲加強低頻信號的能量，使其非線性特征更加突出，由此得到的廣義維數(shù)更能反映信號的特點，使故障判斷更加容易。

4 結(jié)論

針對噪聲背景下分形維數(shù)無法準(zhǔn)確反映信號非線性特性和吸引子復(fù)雜情況的問題，提出了先對含噪信號進行級聯(lián)雙穩(wěn)隨機共振濾波，然后再計算其廣義分形維數(shù)的方法。該方法利用了級聯(lián)雙穩(wěn)隨機共振的濾波特性，不但可以逐步消除高頻噪聲，還可以加強低頻信號的能量，提高信噪比。在此基礎(chǔ)上計算信號的廣義維數(shù)，可以得到更準(zhǔn)確的分形維數(shù)，然后提取盒維數(shù)、信息維數(shù)和關(guān)聯(lián)維數(shù)，構(gòu)成廣義維數(shù)序列，采用相關(guān)性判斷法識別機械設(shè)備故障，達(dá)到故障診斷的目的。

［1］ Mandelbrot B B.The fractal geometry of nature［M］.San Francisco：Freeman，1982.

［2］Falconer K J.Fractal geometry：Mathematical foundation and applications［M］. Chichester：John Wiley ＆ Sons Ltd，1990.

［3］胥永剛，何正嘉.分形維數(shù)和近似熵用于度量信號復(fù)雜性的比較研究［J］.振動與沖擊，2003，22（3）：25-27.

［4］吳蘭蘭，劉儉英，文友先.基于分形維數(shù)的玉米和雜草圖像識別［J］.農(nóng)業(yè)機械學(xué)報，2009，40（3）：176-179.

［5］劉小剛，趙慧潔，李 娜.基于多重分形譜的高光譜數(shù)據(jù)特征提取［J］.光學(xué)學(xué)報，2009，29（3）：844-848.

［6］苑 宇，馬孝江.局域波時頻域多重分形在故障診斷中的應(yīng)用［J］.振動與沖擊，2007，26（5）：60-63.

［7］李國賓，段樹林，于洪亮，等.發(fā)動機振動信號特征參數(shù)的多重分形研究［J］.內(nèi)燃機學(xué)報，2008，26（1）：87-91.

［8］徐玉秀，鐘建軍，聞邦椿.旋轉(zhuǎn)機械動態(tài)特性的分形特征及故障診斷［J］.機械工程學(xué)報，2005，41（12）：186-189.

［9］Benzi R，Sutera A，Vulpiana A.The mechanism of stochastic resonance［J］.J Phys A，1981，14（11）：453-457.

［10］何慧龍，王太勇，冷永剛，等.級聯(lián)雙穩(wěn)隨機共振系統(tǒng)非線性濾波特性［J］.吉林大學(xué)學(xué)報（工學(xué)版），2007，37（4）：905-909.

［11］冷永剛，王太勇，郭 焱.二次采樣用于隨機共振從強噪聲中提取弱信號的數(shù)值研究［J］.物理學(xué)報，2003，52（10）：2432-2437.

［12］趙艷菊，王太勇，任成祖，等.強噪聲背景下的經(jīng)驗?zāi)Ｊ椒纸庋芯浚跩］.振動與沖擊，2009，28（3）：149-151.

［13］樊福梅，梁 平，吳庚申.基于分形盒維數(shù)的汽輪機轉(zhuǎn)子振動故障診斷的實驗研究［J］.核動力工程，2006，27（1）：85-89.

［14］李 琳，張永祥，明廷濤.EMD降噪的關(guān)聯(lián)維數(shù)在齒輪故障診斷中的應(yīng)用研究［J］.振動與沖擊，2009，28（4）：145-148.

［15］李 娜，方彥軍.利用關(guān)聯(lián)維數(shù)分析機械系統(tǒng)故障信號［J］.振動與沖擊，2007，26（4）：136-139.