基于聲信號的滾動軸承故障診斷*

李宏亮，黃　民，高　宏，馬　超

(北京信息科技大學(xué) 機電工程學(xué)院，北京　100192)

?

基于聲信號的滾動軸承故障診斷*

李宏亮，黃民，高宏，馬超

(北京信息科技大學(xué) 機電工程學(xué)院，北京100192)

對現(xiàn)有振動傳感器存在安裝困難等問題，提出了基于聲信號的滾動軸承故障診斷方法；利用聲學(xué)傳感器，采用非接觸方式采集故障軸承聲信號；通過Matlab設(shè)計線性相位FIR濾波器，得到所需頻帶的沖擊響應(yīng)序列；通過Hilbert變換，得到被故障源調(diào)制信號的包絡(luò)信號；通過對包絡(luò)信號進(jìn)行頻譜分析，確定滾動軸承的故障類型。利用該方法對實驗臺故障軸承進(jìn)行診斷，得到了準(zhǔn)確的診斷結(jié)果，驗證了方法的可行性與有效性。

滾動軸承；故障診斷；FIR濾波器；Hilbert變換

0　引言

對于滾動軸承的故障診斷，最關(guān)鍵、最困難的問題之一是如何對滾動軸承的故障特征頻率進(jìn)行正確的提取，它直接關(guān)系到滾動軸承故障診斷的準(zhǔn)確性和預(yù)報的可靠性。為了解決這一問題，國內(nèi)外學(xué)者主要通過采集滾動軸承運行過程中的振動信號、聲發(fā)射信號，運用信號分析與處理的理論和技術(shù)手段來獲取更多信息。目前，應(yīng)用在滾針軸承故障診斷中信號分析與處理的數(shù)學(xué)方法有基于傅里葉變換頻譜分析方法、小波變換和EMD等分析方法[1-3]。其中，Hilbert解調(diào)技術(shù)作為重要的包絡(luò)分析方法也被廣泛應(yīng)用[4-8]。然而，振動信號、聲發(fā)射信號的采集都為接觸式。在實際應(yīng)用中，有些設(shè)備并不能為傳感器安裝提供合適的安裝場合，如處于高溫、高腐蝕、高濕度等較惡劣環(huán)境下設(shè)備；滾動軸承運轉(zhuǎn)過程中位置需要移動的設(shè)備等。所以，使得基于振動信號、聲發(fā)射信號的采集與分析方法受到限制，需通過非接觸式方式進(jìn)行數(shù)據(jù)采集。而聲信號同振動信號一樣，同樣包含著滾動軸承運行狀態(tài)的重要信息。

因此，本文提出了一種基于聲信號、采用Hilbert包絡(luò)解調(diào)技術(shù)的滾動軸承故障診斷方法。通過非接觸式進(jìn)行數(shù)據(jù)采集，不影響設(shè)備正常操作，是一種簡易、快捷的診斷方法。通過對故障軸承聲信號進(jìn)行Hilbert包絡(luò)解調(diào)分析，找出滾動軸承故障特征頻率，從而確定滾動軸承故障類型。

1　滾動軸承故障聲信號診斷原理

當(dāng)滾動軸承發(fā)生故障時，由故障引起的激振力會激發(fā)軸承振動，在軸承振動的同時，往往還會激發(fā)出聲音信號，如果滾動軸承停止運轉(zhuǎn)，振動也會隨之停止，聲音也會消失，這就說明滾動軸承的振動產(chǎn)生了聲音。其實，聲音就是一種機械波，是物體機械振動通過彈性媒介向外傳播的結(jié)果，通過彈性媒介向外輻射，形成聲場。通過傳聲器，采用非接觸式方式采集滾動軸承運轉(zhuǎn)過程中的聲信號，其同振動信號一樣，同樣能夠反映出滾動軸承運行的狀態(tài)信息。其聲學(xué)檢測系統(tǒng)如圖1所示。

圖1　滾動軸承故障聲學(xué)檢測系統(tǒng)

圖1所示為滾動軸承故障聲學(xué)檢測系統(tǒng)，其主要通過非接觸式測量，將傳聲器正對滾動軸承，放置在承輻射出聲場范圍內(nèi)的某點，采集該點處的聲壓信號。該方法的主要特點是信號易取，采集方法和設(shè)備簡單；無需粘連傳感器，可對運動設(shè)備(如火車滾動軸承)進(jìn)行監(jiān)測。彌補了傳統(tǒng)基于振動信號和聲發(fā)射信號檢測軸承運行狀態(tài)中傳感器無法安裝等問題。

2　滾動軸承故障聲信號診斷方法

基于聲信號滾動軸承故障診斷方法同其他診斷方法基本一致，只是獲取信號的方式不同。當(dāng)滾動軸承某一部件出現(xiàn)故障，采集到的聲信號往往都是被調(diào)制了的信號，在診斷過程中，只需要對被故障源調(diào)制了的信號進(jìn)行分析，其它的信號可作為噪聲信號進(jìn)行濾除。因此，可通過線性相位FIR濾波器對原始聲信號進(jìn)行處理，然后對濾波處理后的信號進(jìn)行Hilbert變換，獲得Hilbert包絡(luò)信號，通過對Hilbert包絡(luò)信號進(jìn)行譜分析，判斷滾動軸承是否出現(xiàn)故障。其中，基于Hilbert包絡(luò)解調(diào)分析方法提取滾動軸承故障特征頻率原理圖如圖2所示。

圖2　Hilbert包絡(luò)解調(diào)分析方法原理圖

2.1線性相位FIR濾波器設(shè)計

線性相位FIR濾波器是指相位響應(yīng)為頻率的線性函數(shù)的FIR濾波器[9]。如果用振幅響應(yīng)和相位響應(yīng)表示濾波器的頻率響應(yīng)，可表示為：

H(ejw)=A(w)eφ(w)

(1)

線性相位響應(yīng)表示為：

φ(w)=α-wτ

(2)

其中：A(w)為振幅響應(yīng)，α為常數(shù)，表示初始相位；

在設(shè)計線性相位FIR濾波器時，需要充分利用其振幅響應(yīng)，因為不同的振幅響應(yīng)對濾波器的類型是有要求的。不同的類型，對應(yīng)著不同的沖擊響應(yīng)。根據(jù)不同的沖擊響應(yīng)序列可分為4類，如表1所示。

表1　4類線性相位FIR濾波器的特點

通過Matlab軟件、采用窗函數(shù)設(shè)計方法進(jìn)行Ⅰ型線性相位FIR濾波器設(shè)計。其中: 帶通濾波器上、下截止頻率為13kHz和14kHz，過渡帶寬為50Hz，阻帶衰減AS為60dB。設(shè)計出的濾波器沖擊響應(yīng)和幅度響應(yīng)如圖3所示。

圖3　濾波器的沖擊響應(yīng)和幅度響應(yīng)

圖3a為濾波器的沖擊響應(yīng)，圖3b為濾波器的幅度響應(yīng)。沖擊響應(yīng)為偶對稱，類型為帶通濾波器，通帶范圍內(nèi)無衰減，阻帶范圍內(nèi)衰減較大，通過原始信號與沖擊響應(yīng)進(jìn)行卷積運算，便可得濾波后的信號。

2.2Hilbert變換原理

Hilbert變換[10]是通過將實信號表示為復(fù)信號，可以由此研究信號的包絡(luò)、瞬時相位和瞬時頻率。設(shè)信號為x(t)，則其包絡(luò)信號a(t)為：

3　實驗信號分析

3.1模擬實驗平臺的搭建

實驗所用設(shè)備為模擬試驗臺，其中包括：驅(qū)動部分、加載部分和故障軸承部分。加載部分通過螺栓進(jìn)行加載。如圖4所示。

圖4　模擬實驗平臺

所采用的軸承為圓柱滾子軸承(N1004)，基本參數(shù)和理論故障特征頻率見表2、表3。其中，表3中fr為回轉(zhuǎn)軸(或內(nèi)環(huán))的旋轉(zhuǎn)頻率。

表2　滾動軸承基本參數(shù)

表3　滾動軸承理論故障特征頻率

滾動軸承故障聲學(xué)檢測系統(tǒng)為非接觸式、單通道采集，對2個軸承分別做一處內(nèi)圈故障和外圈故障。當(dāng)回轉(zhuǎn)軸轉(zhuǎn)速為493r/min、加載工況下、采樣頻率為51.2kHz的情況下，采樣時間2s。采集到滾動軸承內(nèi)、外圈故障的原始聲信號如圖5所示。

圖5　內(nèi)、外圈故障的原始聲信號

圖5a為滾動軸承內(nèi)圈故障的原始聲信號，圖5b為滾動軸承外圈故障的原始聲信號。根據(jù)表3可計算出內(nèi)、外圈理論故障特征頻率分別為58.05Hz和40.55Hz。

3.2內(nèi)、外圈故障分析

對滾動軸承內(nèi)、外圈故障原始聲信號做帶通濾波處理，并對濾波處理后的時域信號進(jìn)行Hilbert變換，根據(jù)變換結(jié)果求包絡(luò)信號。濾波處理后的時域信號和Hilbert包絡(luò)信號分別如圖6、圖7所示。

圖6　內(nèi)圈濾波處理后的時域信號及其Hilbert包絡(luò)信號

圖7　外圈濾波處理后的時域信號及其Hilbert包絡(luò)信號

圖6a、圖7a為帶通濾波處理后的滾動軸承內(nèi)、外圈故障聲信號；圖6b、圖7b為滾動軸承內(nèi)、外圈故障聲信號經(jīng)濾波處理后的Hilbert包絡(luò)信號。對滾動軸承內(nèi)、外圈Hilbert包絡(luò)信號進(jìn)行頻譜分析，分析得到的結(jié)果分別如圖8所示。

圖8a為內(nèi)圈故障Hilbert包絡(luò)譜分析圖，圖8b為外圈故障Hilbert包絡(luò)譜分析圖。由于采樣頻率較高，使得Hilbert包絡(luò)譜分析頻帶較寬。而滾動軸承故障特征頻率主要集中在低頻部分，故圖8所示頻帶范圍為30Hz到400Hz，并分別在此頻帶范圍內(nèi)找出內(nèi)、外圈故障軸承Hilbert包絡(luò)譜分析結(jié)果的前五階極大值，在圖8中用*號標(biāo)出，其頻率和幅值分別見表4。

表4　Hilbert包絡(luò)譜前五階極大值

由表4可知：內(nèi)、外圈的Hilbert包絡(luò)譜在頻率58.5Hz和41Hz處能量最大，內(nèi)、外圈理論故障特征頻率為58.05Hz和40.55Hz，二者相差0.45Hz，而Hilbert包絡(luò)譜的分辨率為0.5Hz。從分析結(jié)果來看，可判定軸承內(nèi)、外圈出現(xiàn)了故障。由此可見，利用滾動軸承運行過程中的聲信號，采用包絡(luò)解調(diào)方法也能實現(xiàn)對滾動軸承早期故障的準(zhǔn)確診斷。

4　結(jié)束語

本文應(yīng)用FIR濾波、Hilbert變換以及頻譜分析對軸承聲音信號進(jìn)行分析，通過實驗數(shù)據(jù)分析能夠判斷軸承的狀態(tài)，該方法彌補了傳統(tǒng)的基于振動和聲發(fā)射信號檢測軸承狀態(tài)中傳感器安裝的問題。采用非接觸方式采集滾動軸承運行的聲信號，操作簡單、快捷，不影響設(shè)備正常運行；利用Hilbert包絡(luò)解調(diào)同樣能夠有效的將故障特征從復(fù)雜的、被調(diào)制的聲信號中提取出來；通過以實驗室故障軸承為例，驗證了該方法的可行性與有效性；但同振動信號相比，滾動軸承聲信號受環(huán)境噪聲干擾較大、信噪比相對較低，在某些測試場合，單一聲學(xué)傳感器不能滿足測試需求，因此，可通過聲陣列等其他方式進(jìn)行聲信號采集與故障診斷，提高信噪比。

[1] 付元華,羅仁澤. 小波包樣本熵灰色關(guān)聯(lián)度軸承故障診斷[J]. 組合機床與自動化加工技術(shù),2015(7):128-130,134.

[2] 王磊,張清華,馬春燕,等. 基于GA和EMD包絡(luò)譜分析的滾動軸承故障診斷[J]. 組合機床與自動化加工技術(shù),2014(4):53-56,59.

[3] 姜濤,袁勝發(fā).基于改進(jìn)小波神經(jīng)網(wǎng)絡(luò)的滾動軸承故障診斷[J].華中農(nóng)業(yè)大學(xué)學(xué)報,2014,33(1):131-136.

[4] 郭福平,段志宏,孫志偉. 基于包絡(luò)譜分析的滾動軸承滾動體故障聲發(fā)射診斷研究[J]. 組合機床與自動化加工技術(shù),2015(2):104-106,109.

[5] 代拓,花云浩,王林.基于Hilbert 變換的滾動軸承故障診斷[J].煤礦機械,2012,33(2):255-256.

[6] 黃中華,謝雅.基于Hilbert 變換的滾動軸承內(nèi)環(huán)和外環(huán)故障診斷[J].中南大學(xué)學(xué)報(自然科學(xué)版),2011,42(7):1992-1996.

[7] 陳永會,姜旭,郭山國,等.基于小波分析和Hilbert 變換的滾動軸承故障診斷[J].機械設(shè)計,2010,27(8): 91-94.

[8] 楊超,李亦滔.基于信號預(yù)處理和Hilbert變換的滾動軸承故障診斷[J].華東交通大學(xué)學(xué)報,2012,29(4): 1-4.

[9] 姚天任.數(shù)字信號處理[M].北京：清華大學(xué)出版社,2011.

[10] 何正嘉,陳進(jìn),王太勇,等.機械故障診斷理論及應(yīng)用[M].北京:高等教育出版社,2012.

(編輯李秀敏)

A Fault Diagnosis of Rolling Bearings Based on Acoustic Signals

LI Hong-liang,HUANG Min,GAO Hong,MA Chao

(School of Eleclromechanical Engineering,Beijing Information Science & Technology University,Beijing 1000192,China)

To solve difficult of installing the vibration sensor, a fault diagnosis method for rolling bearings was proposed in this paper, which is based on acoustic signals. Acoustic sensor was used to acquisition acoustic signals of the fault bearing by the non-contact way. The impulse response sequence of Linear-phase FIR filter was calculated using Matlab. The envelope signal was obtained by Hilbert transform for the fault source modulated signal. The fault type of rolling bearing was determined by analyzing the frequency spectrum of the signal envelope. The method was applied to fault diagnosis for rolling bearing，and the feasibility and effectiveness of method were verified through examples of laboratory.

rolling bearing;fault diagnosis;linear-phase FIR filter;Hilbert transform

1001-2265(2016)07-086-03DOI:10.13462/j.cnki.mmtamt.2016.07.024

2015-08-31；

2015-09-19

國家科技重大專項資助項目(2013ZX04011012)

李宏亮(1989—)，男，北京人，北京信息科技大學(xué)碩士研究生，研究方向為機電系統(tǒng)故障診斷，(E-mail)shjshlhl@126.com。

TH17；TG506

A