基于冗余2 代小波局部梯度譜熵的軸承故障診斷

張 園 李 力

(①三峽大學(xué)科技學(xué)院，湖北 宜昌 443002;②三峽大學(xué)機械與材料學(xué)院，湖北 宜昌 443002)

滾動軸承在工業(yè)生產(chǎn)中應(yīng)用非常廣泛，其狀態(tài)常影響整個設(shè)備的精度、可靠性與壽命［1］。冗余第2 代小波降噪能提供更豐富的特征信息可更好劃分頻帶，但是存在頻帶錯位和誤差積累兩類缺陷，改進的冗余第2 代小波包分解可克服這些不足［2-4］。梯度在信號檢測和圖像識別中得到了很好的應(yīng)用，例如平均梯度反映圖像中微小細(xì)節(jié)特征，而形態(tài)梯度常用來在圖像中進行邊沿檢測，而局部梯度在檢測小波包分解信號突變點方面也得到了一定應(yīng)用［5-7］。本文將改進冗余第2 代小波包分解與局部梯度方法結(jié)合，并提出一種局部梯度譜熵指標(biāo)，對試驗滾動軸承進行故障診斷，能進行較好分類。

1 改進冗余第2 代小波包分解

如圖1 所示其第1 層分解過程，原信號為x 經(jīng)過1 層分解，獲得逼近信號s⌒(1)與細(xì)節(jié)信號d⌒(1)，標(biāo)準(zhǔn)化系數(shù)分別為a1和b1。改進冗余第2 代小波分解得到的逼近信號c(1)與細(xì)節(jié)信號r(1)。用代替原信號得到下一層標(biāo)準(zhǔn)化系數(shù)，依此獲得不同分解層數(shù)下的標(biāo)準(zhǔn)化系數(shù)。

改進冗余第2 代小波包分解結(jié)果總能量與原始信號能量完全相等。因此，通過標(biāo)準(zhǔn)化系數(shù)消除了冗余第2 代小波變換的誤差傳播［6］。

2 基于改進冗余第2 代小波的滾動軸承振動信號的分解

2.1 滾動軸承振動信號的采集

由實驗臺采集的6308 軸承振動信號，其采樣頻率20 kHz。軸承故障類型為:外圈剝落、內(nèi)圈剝落、滾動體剝落和正常［7］。在4 種轉(zhuǎn)速、3 種載荷下運轉(zhuǎn)對軸承進行采樣，如表1 所示不同轉(zhuǎn)速的軸承故障特征頻率。圖2 為轉(zhuǎn)速n1、輕載下的時域圖。

表1 滾動軸承的故障特征頻率(Hz)

2.2 改進冗余第2 代小波的滾動軸承信號分解

如圖3 所示，對圖2a 信號進行冗余第2 代小波3層分解，得到8 頻帶載波，其第2～4 頻帶的與原信號的幅值較接近，其中4 頻帶最接近為200 mV 左右，每隔0.04 s 有一次脈沖，其頻率約等于特征頻率23.7 Hz。

為比較不同頻帶載波和原信號的相關(guān)性，分別計算各載波和原信號的相關(guān)系數(shù)。

對圖2 所示軸承信號的8 個頻帶載波信號與源信號計算其相關(guān)系數(shù)，如表2 所示。對于故障軸承，相關(guān)性最大的頻帶都在第4 頻帶，其值均大于0.75。而正常軸承最大頻帶都為第1 頻帶，其值大約為0.9，說明故障軸承第4 頻帶、正常信號第1 頻帶與源信號高度相似。對多組信號試驗，結(jié)果同上。

表2 滾動軸承8 頻帶互相關(guān)系數(shù)

3 滾動軸承載波信號局部梯度圖

3.1 局部梯度圖

設(shè)備運行狀態(tài)的重要信息信號往往體現(xiàn)為不規(guī)則的突變點，局部梯度可較準(zhǔn)確識別小波包分解后的突變點。

信號設(shè)為x={x(i)，i=1，2，…，n}，n 為數(shù)據(jù)長度，定義第i 點的局部梯度向量為

vi(j)=x(i)-x(i+j)，i=1，…，n，j=-r，…，r

可得到第i 個點相對于其相鄰的第j 個點的變化，整個信號序列的局部梯度向量表示為

信號的局部梯度圖定義為［6］

3.2 滾動軸承頻帶載波的局部梯度圖分析

對圖2 所示信號的最優(yōu)載波求局部梯度，得到圖4a～4c 所示故障軸承都可以看到存在周期性的脈沖，其中圖4a 約隔0.04 s 有1 次脈沖，其頻率正好接近于其特征頻率23.7 Hz，而且最高的兩脈沖間時間間隔約為0.125 s，其頻率接近于轉(zhuǎn)速1 轉(zhuǎn)頻7.7 Hz;圖4b約0.025 s 有1 次脈沖，頻率接近于特征頻率37.9 Hz;圖4c 約0.033 s 有1 次脈沖，其頻率約其特征頻率雙倍頻31.4 Hz;如圖4d 所示，正常軸承的局部梯度圖無周期性的脈沖。

為比較不同轉(zhuǎn)速下故障特征的提取效果，對輕載下n1～n4下不同轉(zhuǎn)速的局部梯度求功率譜，得到圖5～8。如圖5 所示，在n1～n4的局部梯度譜上，分別在圖示24、39、56、80 Hz 及其倍頻處脈沖，接近于外圈剝落的特征頻率及倍頻;同樣，如圖6 所示，在39、61、90、131 Hz 及其倍頻處脈沖，接近于內(nèi)圈剝落的特征頻率及倍頻;同樣，如圖7 所示，在31、51、75、109 Hz及其倍頻處脈沖，接近于滾動體剝落的特征頻率雙倍頻;圖8 中正常軸承局部梯度功率譜高頻段沒有脈沖，其4 種轉(zhuǎn)速下的局部梯度譜脈沖呈隨機分布狀，有些突出的脈沖有可能是轉(zhuǎn)頻的倍頻。綜上所述，采用改進冗余第2 代小波包分解的局部梯度譜可以提取不同轉(zhuǎn)速下的故障特征頻率，正常軸承無缺陷，其局部梯度譜呈隨機分布。

4 基于局部梯度譜熵的滾動軸承故障分類

從以上分析可以看到，最優(yōu)載波局部梯度能夠提取滾動軸承的故障信息的，可以此作為特征指標(biāo)對信號進行分類。

利用局部梯度功率譜的沖擊脈沖可識別故障特征，而信息熵可用來衡量信號不確定性，熵值越小，信息越確定;反之，越大，信息越不確定［8-10］。

如圖9 所示，取不同轉(zhuǎn)速下無載、輕載、重載下的共48 組試驗軸承載波信號，求其局部梯度信息熵，結(jié)果為:正常軸承的局部梯度熵總大于故障軸承，原因是正常軸承信號往往呈隨機分布，不確定性大;但是不同故障類型的軸承之間無法用該指標(biāo)區(qū)分。

如圖10 所示，對其局部梯度譜幅值求熵值，正常軸承遠遠大于故障軸承，而故障軸承的熵值大小順序為:滾動體剝落＞內(nèi)圈剝落＞外圈剝落

由圖5～7 所示，分解后的載波信號局部梯度譜呈周期性波動，且外圈剝落的軸承周期性脈沖最明顯，其信號越確定，熵值越小;內(nèi)圈剝落和滾動體剝落其次，故其熵值稍小。而圖8 所示正常軸承信號局部梯度譜呈隨機波動，信號不確定，故其熵值最大;且從圖9、10中可以看到，n1～n4四種轉(zhuǎn)速下不同載荷下波形，圖9所示的局部梯度熵隨載荷和轉(zhuǎn)速增加，熵值稍有增加;而圖10 所示局部梯度譜熵隨載荷變動不大。比較這兩種指標(biāo)，局部梯度譜熵指標(biāo)比局部梯度熵的試驗，軸承故障分類的效果要好。

5 結(jié)語

對試驗軸承信號用改進冗余第2 代小波分解，得8 頻帶載波。根據(jù)互相關(guān)系數(shù)最大原則獲得最優(yōu)載波。求其局部梯度及功率譜，可從故障軸承局部梯度圖中獲得其特征頻率所對應(yīng)的周期性脈沖，其局部梯度譜可以獲得不同故障特征頻率及其諧波，而正常軸承波形隨機分布。同時，提出基于載波信號的局部梯度譜熵指標(biāo)，可以將試驗工況的軸承信號區(qū)分開。

［1］丁建明，林建輝，楊強，等.基于諧波小波奇異熵的軸承故障實時診斷［J］.中國機械工程，2010，21(1):55 -58.

［2］史斐娜，劉廣璞.基于冗余第二代小波變換的齒輪箱故障診斷［J］.煤礦機械，2014，35(5):248 -250.

［3］高立新，湯文亮，青永剛，等.基于冗余第二代小波的降噪技術(shù)［J］.北京工業(yè)大學(xué)學(xué)報，2008(12):1233 -1237.

［4］羅榮，田福慶，馮昌林，等.改進的冗余第二代小波包及其故障診斷應(yīng)用［J］.華中科技大學(xué)學(xué)報自然科學(xué)版，2014，42(5):40 -46.

［5］李兵，張培林，任國全，等.基于形態(tài)梯度算子的滾動軸承故障特征提取［J］.中北大學(xué)學(xué)報:自然科學(xué)版，2011，32(4):426 -430.

［6］何正佳，袁靜，訾艷陽，等.機械故障診斷的內(nèi)積變換原理與應(yīng)用［M］.北京:科學(xué)出版社，2012:112 -128.

［7］陳超，江濤，徐文學(xué).基于平均梯度的多聚焦圖像融合方法研究［J］.山東建筑大學(xué)學(xué)報，2009，24(2):133 -136.

［8］朱瑜，王殿，王海洋.基于EMD 和信息熵的滾動軸承故障診斷［J］.軸承，2012(6):50 -53.

［9］劉韜，陳進，董廣明.基于頻帶熵的滾動軸承故障診斷研究［J］.振動與沖擊，2014，33(1):132 -135.

［10］劉濤，黃成東.基于信息熵的渦旋壓縮機故障診斷研究［J］.壓縮機技術(shù)，2012，231(1):1 -8.