張 良，楊 濤

(1.綿陽師范學(xué)院 機電工程學(xué)院，四川 綿陽 621000；2.西南科技大學(xué) 信息工程學(xué)院，四川 綿陽 621010)

基于多重自相關(guān)與包絡(luò)譜分析的滾動軸承故障診斷方法研究*

張 良1，楊 濤2

(1.綿陽師范學(xué)院 機電工程學(xué)院，四川 綿陽 621000；2.西南科技大學(xué) 信息工程學(xué)院，四川 綿陽 621010)

滾動軸承在故障狀態(tài)運行時，傳感器測得的振動信號為非平穩(wěn)、多分量的調(diào)制信號。在故障出現(xiàn)早期，由于調(diào)制信號微弱且含有噪聲，導(dǎo)致故障特征難以識別，采用多重自相關(guān)消除噪聲干擾，提取信號中的周期調(diào)制成分，然后利用Hilbert 變換的包絡(luò)解調(diào)方法獲取故障特征頻率，從而判斷出軸承故障類型。實驗結(jié)果表明，采用多重自相關(guān)與包絡(luò)譜解調(diào)相結(jié)合的方法，能較準(zhǔn)確的提取滾動軸承故障特征頻率，具有一定的工程應(yīng)用價值。

滾動軸承；多重自相關(guān)；Hilbert變換；故障診斷

0 引言

據(jù)統(tǒng)計，由于軸承導(dǎo)致的故障在旋轉(zhuǎn)機械所有故障中占到 1/3 左右[1]，因此對重要軸承進(jìn)行工況監(jiān)視與故障診斷，不但可以減少或杜絕事故發(fā)生，而且還能最大限度地發(fā)揮軸承工作潛力，這在提高生產(chǎn)效益以及保障生產(chǎn)安全方面都具有十分重要的意義[2]。滾動軸承出現(xiàn)局部損傷性故障的特點是微弱的周期信號寬帶沖擊會激起系統(tǒng)結(jié)構(gòu)自身的高頻振動[3]，同時滾動軸承振動信號的頻率范圍較寬，信噪比通常較低，局部損傷性所引起的沖擊比較微弱，直接對振動信號解調(diào)時調(diào)制信息容易淹沒在噪聲中[4]，因此如何從微弱的故障振動信號中提取沖擊所激起的高頻振動信息進(jìn)行解調(diào)分析是軸承故障診斷的關(guān)鍵[5]。目前，用于在低信噪比的條件下提取滾動軸承振動信號包絡(luò)譜的主要方法有形態(tài)分量分析[5-6]、時延相關(guān)解調(diào)[7]、小波降噪[8-9]、局部均值分解[10-11]、經(jīng)驗?zāi)B(tài)分解[12-13]、盲源分離[14-15]等，上述方法在故障振動信號的包絡(luò)譜解調(diào)與提取方面均有一定的效果。本文提出一種基于多重自相關(guān)的包絡(luò)譜滾動軸承故障診斷方法，對原始振動信號首先進(jìn)行多重自相關(guān)運算，最大限度地消除信號中混有的高斯噪聲，然后在對多重自相關(guān)運算結(jié)果做帶通濾波，將濾波截取的高頻共振頻段進(jìn)行希爾伯特變換，得到包絡(luò)譜，從而提取滾動軸承的故障特征頻率，完成了對滾動軸承的故障診斷，通過實例分析驗證了該方法的有效性。

1 理論基礎(chǔ)

1.1 滾動軸承的頻率特征

滾動軸承常見的故障形式有表面損傷、疲勞剝落、過熱燒傷、電化學(xué)腐蝕、膠合及破損等[16]?？蓪⑤S承元件的這些故障歸結(jié)為3種形式:內(nèi)圈故障、外圈故障以及滾動體故障。當(dāng)軸承元件(內(nèi)圈、外圈、滾動體)出現(xiàn)了局部故障時，在轉(zhuǎn)動過程中會與其它相接觸的元件發(fā)生反復(fù)碰撞，從而產(chǎn)生沖擊振動。這種短時周期性的低頻脈沖信號頻率通常稱為滾動軸承的通過頻率，即故障特征頻率[17]。

理想情況下，滾動軸承外圈固定，內(nèi)圈隨轉(zhuǎn)軸旋轉(zhuǎn)，滾動體與滾道之間沒有相對滑動，且軸向、徑向在承受載荷各部件無變形時，由軸承的結(jié)構(gòu)幾何參數(shù)可計算出各軸承元件的理論故障特征頻率[18]，其計算表達(dá)式如下：

內(nèi)圈故障特征頻率fi為：

(1)

外圈故障特征頻率fo為：

(2)

滾動體故障特征頻率foc為：

(3)

保持架故障特征頻率fc：

(4)

上式中各個參數(shù)含義:D為軸承節(jié)徑;d為滾動體直徑;α為軸承接觸角;z為滾動體個數(shù);N為軸承轉(zhuǎn)速。

而在實際的故障診斷中，由于各種干擾因素的影響，雖然理論故障特征頻率值與工程診斷分析得到的檢測值有一定的偏差[19]，但故障特征頻率的提取與分析仍然是判斷軸承故障十分重要和有效的理論依據(jù)。

1.2 多重自相關(guān)

定義rxy(τ)為信號x(t)和y(t)的互相關(guān)函數(shù)，其中τ為延遲時間。

(5)

如果設(shè)y(t)=x(t)，則上面定義的互相關(guān)函數(shù)rxy(τ)變?yōu)樽韵嚓P(guān)函數(shù)rxx(t)：

(6)

rxx(t)反應(yīng)了信號x(t)和自身做了一段延遲時間之后x(t+τ)的相似程度[20]。設(shè)信號x(t)=s(t)+n(t)，其中s(t)為周期信號，n(t)為高斯白噪聲，將x(t)代入式(6)做自相關(guān)運算得：

(7)

將式(7)展開后得到：

rxx(τ)=rss(τ)+rsn(τ)+rnn(τ)

(8)

由相關(guān)函數(shù)性質(zhì)[20]可知，周期信號經(jīng)過自相關(guān)運算后，rss(τ)周期性與原來的周期性函數(shù)s(t)一致，但周期信號與噪聲信號的互相關(guān)函數(shù)rsn(τ)，噪聲信號的自相關(guān)函數(shù)rnn(τ)均趨近于0。因此利用以上結(jié)論，當(dāng)含有高斯噪聲的信號經(jīng)過自相關(guān)運算后，可以對信號進(jìn)行去噪和提取周期成分。

理論上，當(dāng)測量時間T為無窮大時，rsn(τ)、rnn(τ)將趨于0，但在實際工程中，T不可能取到無窮大，所以rsn(τ)、rnn(τ)不等于0，去噪情況不會太理想。采用以下步驟進(jìn)行多次自相關(guān)運算，可以減少高斯噪聲，提高信噪比。

(1)設(shè)

xi(t)=si(t)+ni(t)

(9)

其中，

si(t)=rss(i)(τ)，ni(t) =rsn (i)(τ) +rnn (i)(τ)

(2)把式(9)帶入式(6)做自相關(guān)；

(3)重復(fù)i次步驟(1)和步驟(2)。

這里把周期信號的自相關(guān)函數(shù)rss(τ)i看作新的周期信號，而噪聲和周期信號的互相關(guān)函數(shù)rsn(τ)i、噪聲的自相關(guān)函數(shù)rnn(τ)i看作新的噪聲信號，i為自相關(guān)次數(shù)。對信號xi(t)做多次自相關(guān)可以有效的提高信噪比，從而更好的檢測出淹沒在高斯噪聲中的微弱的周期信號[21-22]。

1.3 包絡(luò)譜分析

采用Hilbert 變換實現(xiàn)對振動信號進(jìn)行包絡(luò)解調(diào)，進(jìn)而進(jìn)行頻譜分析是滾動軸承故障診斷常用的方法[23-24]。

(10)

(11)

z(t)為信號x(t)的解析信號。如果x(t)受到振幅調(diào)制，則：

(12)

(13)

A(f)為振幅調(diào)制信號x(t)的包絡(luò)譜。通過對滾動軸承振動信號進(jìn)行上述包絡(luò)譜的分析，可以提取滾動軸承的故障特征頻率。

2 實例分析

為驗證基于多重自相關(guān)的包絡(luò)譜分析方法在滾動軸承故障診斷中的有效性，實例數(shù)據(jù)來源于美國西儲大學(xué)在圖1所示的實驗平臺上得到故障軸承數(shù)據(jù)庫。從該數(shù)據(jù)庫中選取電機驅(qū)動端 SKF深溝球軸承故障數(shù)據(jù)，其采樣頻率為12kHZ、電機負(fù)載2HP、電機轉(zhuǎn)速1750r/min、故障直徑0.007英寸。驅(qū)動端 SKF軸承尺寸信息如表1所示。

表1 驅(qū)動端軸承參數(shù)

根據(jù)驅(qū)動端軸承的結(jié)構(gòu)幾何參數(shù)可計算出各軸承元件的理論故障特征頻率為表2所示。

表2 驅(qū)動端軸承各部件故障頻率

圖1 美國西儲大學(xué)故障診斷實驗平臺

本文以外圈與內(nèi)圈故障為例，通過表2可得內(nèi)圈故障特征頻率為157.5Hz，外圈故障特征頻率為105Hz。圖2是滾動軸承內(nèi)圈和外圈發(fā)故障時的時頻圖。

圖2 振動信號時頻圖

從圖2可以看出原始信號時域波形復(fù)雜，頻譜成分豐富，故障特征頻率在頻譜圖中沒有突出顯示，已淹沒在噪聲中。利用1重自相關(guān)和3重自相關(guān)分別對原始信號進(jìn)行濾波，結(jié)果如圖 3和圖4所示。

圖3 1重自相關(guān)去噪后的時頻圖

圖4 3重自相關(guān)去噪后的時頻圖

對比圖3和圖4的頻譜圖可明顯看出3重自相關(guān)在去除高斯噪聲的效果優(yōu)于1重自相關(guān)，且濾波后的信號在2.5～4.5kHz頻段有明顯的共振峰，對該頻段進(jìn)行包絡(luò)譜分析，其包絡(luò)譜如圖5所示。

圖5 包絡(luò)譜

從圖5內(nèi)圈包絡(luò)譜中可以明顯的看到156.7Hz、314.9Hz、629.9Hz、786.6Hz、1417Hz 5個譜峰；該頻率值和和內(nèi)圈故障的特征頻率值157.5Hz及其2倍頻315Hz、4倍頻630Hz、5倍頻787.5Hz、9倍頻1417.5Hz值非常接近；同樣在圖5的外圈包絡(luò)譜中也可以明顯的看到105.5Hz、314.9Hz、524.4Hz、629.9Hz 4個譜峰，該頻率值和和外圈故障的特征頻率105HZ及其3倍頻315Hz、5倍頻525Hz、6倍頻630Hz值非常接近，從而可以診斷出滾動軸承出現(xiàn)了內(nèi)圈和外圈故障。

3 結(jié)論

本文將多重自相關(guān)與包絡(luò)譜分析應(yīng)用于滾動軸承故障診斷中。針對滾動軸承的故障特征信號容易被噪聲淹沒這一特點，通過對振動信號取多重自相關(guān)，減小高斯噪聲對有用信號的干擾提高了信噪比；對多重自相關(guān)后得到的降噪信號進(jìn)行Hilbert 變換實現(xiàn)了包絡(luò)解調(diào)，從而較準(zhǔn)確提取了滾動軸承的故障特征頻率，完成對滾動軸承的故障檢測，驗證了該方法在滾動軸承故障診斷中的有效性，具有較好的工程應(yīng)用價值。

[1] 韓捷,張瑞林.旋轉(zhuǎn)機械故障機理及診斷技術(shù)[M].北京:機械工業(yè)出版社,1997.

[2] 陳長征,胡立新,周勃,等.設(shè)備振動分析與故障診斷技術(shù)[M].北京:科學(xué)出版社,2007.

[3] 王晶,陳果,郝騰飛.強噪聲背景下的滾動軸承故障微弱信號檢測新方法[J].軸承,2012(3):42-46.

[4] 郭玉杰,杜新定,石峰,等.滾動軸承故障的自相關(guān)函數(shù)包絡(luò)分析[J].軸承,2013(6):55-58.

[5] 陳向民,于德介,李蓉,基于形態(tài)分量分析和包絡(luò)譜的軸承故障診斷[J].中國機械工程,2014,25 (8):1047-1053.

[6] 李豫川,伍星,遲毅林,等.基于形態(tài)濾波和稀疏分量分析的滾動軸承故障盲分離[J]. 振動與沖擊,2011，30(12):170-174.

[7] 孟濤,廖明夫.利用時延相關(guān)解調(diào)法診斷滾動軸承的故障[J].航空學(xué)報,2004,25(1):41-44.

[8] 唐貴基,王曉龍.可調(diào)品質(zhì)因子小波變換在滾動軸承微弱故障特征提取中的應(yīng)用[J].中國電機工程學(xué)報,2016,36(3):746-754.

[9] 張邦成,陳珉珉,尹曉靜,等.基于諧波小波和 Hilbert 的滾動軸承故障診斷方法[J]. 現(xiàn)代制造工程,2015(6):27-30.

[10] 陳亞農(nóng),郜普剛,何田,等.局部均值分解在滾動軸承故障綜合診斷中的應(yīng)用[J].振動與沖擊,2012,31(3):73-78.

[11] 王志武,孫虎兒,劉維雄.基于局部均值分解和奇異值差分譜的滾動軸承故障診斷研究[J].機械科學(xué)與技術(shù),2014,33(9):1340-1344.

[12] 胡愛軍,馬萬里,唐貴基.基于集成經(jīng)驗?zāi)B(tài)分解和峭度準(zhǔn)則的滾動軸承故障特征提取方法[J] .中國電機工程學(xué)報,2012,32(11):106-111.

[13] 武哲,楊紹普,任彬,等.基于噪聲輔助多元經(jīng)驗?zāi)B(tài)分解和多尺度形態(tài)學(xué)的滾動軸承故障診斷方法[J] . 振動與沖擊,2016,35(4):127-133.

[14] 李曉暉,傅攀.基于EEMD的單通道盲源分離在軸承故障診斷中的應(yīng)用[J].中國機械工程,2014,25(7):924-930.

[15] 明陽.基于循環(huán)平穩(wěn)和盲源分離的滾動軸承故障特征提取方法研究[D].上海:上海交通大學(xué),2013.

[16] 何亮.基于EMD 技術(shù)的滾動軸承故障診斷研究[D]. 大連:大連理工大學(xué),2012.

[17] 趙光昌,謝剛,張清華.滾動軸承早期故障診斷技術(shù)[J]. 軸承,2015(6)：56-61.

[18] 時培明,許帥,李培.基于小波去噪和 EEMD 包絡(luò)解調(diào)分析的滾動軸承故障診斷方法[J]. 現(xiàn)代制造工程, 2015(12):12-17.

[19] 王雷.基于流形學(xué)習(xí)的滾動軸承故障診斷若干方法研究[D].大連:大連理工大學(xué),2013.

[20] 胡廣書.數(shù)學(xué)信號處理理論、算法與實現(xiàn)[M].3版.北京:清華大學(xué)出版社,2012.

[21] 李曉鵬,席繼樓,劉超,等.多重自相關(guān)在微弱低頻電磁信號檢測中的應(yīng)用[J]. 地震地磁觀測與研究，2011,32(4):44-49.

[22] 李一兵,岳欣,楊莘元.多重自相關(guān)函數(shù)在微弱正弦信號檢測中的應(yīng)用[J].哈爾濱工程大學(xué)學(xué)報,2004,25(4):525-528.

[23] McFadden P D，Smith J D . Vibration monitoring of rolling element bearings by the high frequency resonance technique-areview[J].TribologyInternational,1984,17(1):3-10.

[24]RandallR B,Antoni J . Rolling element bearing diagnostics-a tutorial[J]. Mechanical Systems and Signal Processing,2011,25 (2):485-520.

(編輯 李秀敏)

Fault Diagnosis of Rolling Bearings Based on Multi-Layer Autocorrelation and Envelope Analysis

ZHANG Liang1，YANG Tao2

(1.School of Mechanical and Electrical Engineering,Mianyang Normal University,Mianyang Sichuan,621000,China;2.School of Information Engineering,Southwest University of Science and Technology,Mianyang Sichuan 621010,China)

The vibration signal measured by the sensor is non-stationary and multi-component modulation signal when rolling bearings run in a fault condition. It is difficult to identify the characteristic in the presence of early bearing faults because the modulation signal is weak and polluted by noise. The multi-layer autocorrelation is used to eliminate noise while extracting signal cycle modulation component，and the envelope demodulation method based on Hilbert transform is used to obtain fault characteristic frequency and determine the type of bearing failure.The results show that the method of multi-layer autocorrelation and envelope demodulation can extract the characteristic frequency of rolling bearing more accurately, and has certain engineering application value.

rolling bearing ;multi-layer autocorrelation;Hilbert-transformation;fault diagnosis

1001-2265(2017)08-0093-04

10.13462/j.cnki.mmtamt.2017.08.023

2016-11-16

國家自然科學(xué)基金(F011102);四川中煙工業(yè)責(zé)任有限公司公司科技項目(川渝煙工技研[2015]62號)

張良(1982—)，男，四川綿陽人，綿陽師范學(xué)院助教，碩士，研究方向為機電系統(tǒng)故障診斷、工業(yè)自動化，(E-mail)20015601@163.com。

TH133.3；TG506

A