馮昌林,羅 榮,劉 輝
(1.海軍裝備研究院,北京 100161;2.海軍工程大學(xué)兵器工程系,湖北武漢 430033;3.海南省軍區(qū)軍械修理所,海南定安 571231)
巴特沃斯小波在軸承早期故障診斷中的應(yīng)用
馮昌林1,羅 榮2,劉 輝3
(1.海軍裝備研究院,北京 100161;2.海軍工程大學(xué)兵器工程系,湖北武漢 430033;3.海南省軍區(qū)軍械修理所,海南定安 571231)
設(shè)計(jì)了基于巴特沃斯數(shù)字濾波器的巴特沃斯小波及其濾波器組,并從小波的卷積型定義出發(fā),利用零相位濾波原理實(shí)現(xiàn)了沒(méi)有下抽樣環(huán)節(jié)的具有零相位特性的巴特沃斯小波變換快速算法,成功消除濾波器的非線性相位對(duì)分解結(jié)果的影響。相對(duì)于常用正交小波濾波器組,巴特沃斯小波濾波器組不僅具有更好的頻帶分離特性,還可通過(guò)改變尺度參數(shù)來(lái)改變對(duì)信號(hào)的頻帶劃分的方式,從而將原信號(hào)高頻帶分解在更精細(xì)的頻帶范圍內(nèi),克服了常用正交小波濾波器組頻帶分離特性不夠好,并總按二進(jìn)方式劃分頻帶、頻率分析精度顯得不足的固有缺陷。實(shí)例分析證明:巴特沃斯小波濾波器組相比于其他常用的正交小波濾波器組具有更好的頻帶分離特性和更強(qiáng)的微弱故障特征提取能力,非常適合于軸承早期故障診斷。
故障診斷;巴特沃斯小波;包絡(luò)解調(diào);軸承
小波分析是近些年才廣泛應(yīng)用于振動(dòng)信號(hào)處理的時(shí)頻分析方法。它能同時(shí)提供振動(dòng)信號(hào)的時(shí)域和頻域的局部化信息,具有多尺度和“數(shù)學(xué)顯微”特性,這使得小波分析非常適合于機(jī)械故障特征的提?。?]。在機(jī)械故障診斷領(lǐng)域所用的小波變換通常采用經(jīng)典的離散正交小波變換,并利用著名的Mallat算法實(shí)現(xiàn)其快速運(yùn)算。然而,Mallat算法存在下采樣環(huán)節(jié),分解的長(zhǎng)度隨分解層數(shù)的增加而依次減半,時(shí)間分辨率越來(lái)越低,導(dǎo)致信號(hào)失真,難以提供準(zhǔn)確的頻域局部化信息,不利于機(jī)械故障信號(hào)的特征提取。若用Mallat算法分解之后再對(duì)感興趣的頻帶進(jìn)行單支重構(gòu),不僅增加了計(jì)算量[2],還會(huì)使信號(hào)產(chǎn)生畸變;同時(shí),Mallat算法分解后的各個(gè)子帶的重構(gòu)信號(hào)是存在頻率混疊的,并不是原信號(hào)的真實(shí)分量[3],這極大地影響了分析結(jié)果的可讀性,頻率混疊嚴(yán)重時(shí)還會(huì)造成誤判;最后,經(jīng)典離散正交小波對(duì)應(yīng)的濾波器組是共軛正交濾波器組,其對(duì)信號(hào)的頻帶總是按二進(jìn)方式劃分,缺乏靈活性,且僅相當(dāng)于振動(dòng)和噪聲信號(hào)分析中的1/1頻程分析,對(duì)高頻段的頻帶劃分非常粗糙,頻率分辨率較低,對(duì)提取高頻段的故障特征很不利。
為了更好地將小波變換應(yīng)用到機(jī)械故障診斷中,首先應(yīng)該設(shè)計(jì)幅頻特性好(即幅頻響應(yīng)具有陡峭的頻率截止特性)且截止頻率可以隨尺度變化而改變的功率互補(bǔ)小波濾波器組。從理論上講,與常用的db小波、雙正交小波等FIR濾波器相比,IIR濾波器具有更好的阻帶衰減,用較低的階數(shù)就可以獲得比高階FIR濾波器更好的幅頻特性,可以構(gòu)造出幅頻特性更好的小波濾波器組[4];另外,使用IIR濾波器可以方便實(shí)現(xiàn)對(duì)頻帶劃分方式更加靈活的小波變換,克服經(jīng)典離散小波濾波器組總按1/1頻程分析信號(hào)的不足。其次應(yīng)該摒棄經(jīng)典的Mallat算法,采取沒(méi)有下抽樣環(huán)節(jié)的離散小波變換算法。文獻(xiàn)[5]研究的卷積小波變換快速算法,文獻(xiàn)[6-7]研究的非抽樣小波變換等均是沒(méi)有下抽樣環(huán)節(jié)小波變換算法。筆者直接從小波的卷積型定義出發(fā),利用功率互補(bǔ)巴特沃斯濾波器組和零相位濾波原理實(shí)現(xiàn)了另外一種沒(méi)有下抽樣環(huán)節(jié)、具有零相位特性、且頻帶劃分方式更加靈活的巴特沃斯小波變換,并驗(yàn)證了零相位巴特沃斯小波變換在實(shí)際故障診斷中的有效性。
從小波的卷積定義出發(fā),對(duì)尺度參數(shù)s和時(shí)間參數(shù)t離散化,得到的離散小波變換的計(jì)算公式為
取不同尺度sj的展縮小波ψsj(n)與信號(hào)x(n)進(jìn)行卷積,就可得到信號(hào)的各階小波系數(shù)。由于每個(gè)展縮小波是一個(gè)隨尺度增大,中心頻率向低頻方向移動(dòng),頻域窗寬減小的帶通濾波器,所以小波變換的實(shí)質(zhì)是由展縮小波系{ψsj(n)}構(gòu)成的濾波器組對(duì)信號(hào)進(jìn)行多頻帶分解。又由于{ψsj(n)}構(gòu)成的濾波器組的頻率通帶是由高至低連續(xù)分布的,在實(shí)際應(yīng)用中為減少計(jì)算量,小波分解可采用基于濾波器組的遞歸分解形式實(shí)現(xiàn),如圖1所示。實(shí)現(xiàn)小波變換的關(guān)鍵是尋找一組截止頻率隨尺度變化的且功率互補(bǔ)的濾波器組{Hj(ω),Gj(ω)},即滿足下列條件的{Hj(ω),Gj(ω)}:
|Hj(ω)|2+|Gj(ω)|2=1,且ω<ωj時(shí),|Gj(ω)|2=0;ω≥ωj時(shí),|Gj(ω)|2=1。
功率互補(bǔ)的巴特沃斯濾波器組設(shè)計(jì)已有成熟的理論[8],其與db30小波對(duì)頻帶的劃分比較如圖2所示。
由圖2可知,20階巴特沃斯小波濾波器組的幅頻響應(yīng)具有更陡峭的頻率截止特性,其可以通過(guò)改變?yōu)V波器組的階數(shù)調(diào)整頻率截止特性,抑制濾波器組中的頻帶重疊現(xiàn)象,因此頻帶分離能力比db30小波好,更適于微弱故障特征的提?。煌瑫r(shí)巴特沃斯小波還可以根據(jù)信號(hào)特點(diǎn)靈活選擇尺度參數(shù)來(lái)改變頻帶劃分方式,當(dāng)尺度參數(shù)選為2j/3時(shí),巴特沃斯小波濾波器組對(duì)高頻段頻帶的劃分比db30小波濾波器組的二進(jìn)劃分更為精細(xì),克服了二進(jìn)方式劃分頻帶的不足,特別適合于高頻段故障特征的提取。
圖1是小波分解的遞歸分解形式,若圖1中的濾波器組為巴特沃斯小波濾波器組,則可得巴特沃斯小波變換的快速算法:
式中:Sj(n)為尺度j上的近似信號(hào);Wj(n)為尺度j上的細(xì)節(jié)信號(hào);hj(n),gj(n)為巴特沃斯小波濾波器組;J為最大分解尺度。遞推計(jì)算的初值為S0(n)=x(n)。
由于巴特沃斯濾波器的非線性相位,經(jīng)巴特沃斯小波變換得到的各級(jí)小波分量產(chǎn)生的相移隨所處頻帶不同而不同,給判斷各級(jí)小波分量中的模極大值點(diǎn)或過(guò)零點(diǎn)帶來(lái)了麻煩;而且即使是同一尺度上的小波分量,由于存在一個(gè)帶寬,該帶寬范圍內(nèi)的不同頻率成分也會(huì)有不同的相移,從而使小波分量產(chǎn)生一定程度的畸變。對(duì)一般的小波變換而言,如果小波函數(shù)的形狀不是對(duì)稱的,也會(huì)存在這樣的問(wèn)題[9]。解決這一問(wèn)題的方法是采用零相位濾波。零相位濾波的原理很簡(jiǎn)單,將一個(gè)濾波器H(z)和它的逆濾波器H(z-1)串聯(lián)就可以形成零相位濾波[10]。為了消除相位延遲對(duì)小波分量的影響,利用零相位濾波原理得到的具有零相位巴特沃斯小波變換快速算法為
式中rh(n)、rg(n)分別為hj(n)、gj(n)的自相關(guān)函數(shù)。
用一個(gè)故障軸承內(nèi)圈故障的振動(dòng)信號(hào)(如圖3所示)來(lái)比較20階巴特沃斯小波和db30小波對(duì)信號(hào)的頻帶分離特性。分別采用20階巴特沃斯小波(采用零相位濾波快速算法)和db30小波對(duì)信號(hào)分解至第3層。為便于比較兩者的頻帶分離特性,巴特沃斯小波也采用二進(jìn)劃分的方式進(jìn)行分解,即尺度參數(shù)取為2j。這樣巴特沃斯小波和db30小波均得到以下4個(gè)頻帶:W1(理想頻率范圍:3~6 kHz)、W2(理想頻率范圍:1.5~3kHz)、W3(理想頻率范圍:0.75~1.5kHz)和S3(理想頻率范圍:0~0.75kHz),結(jié)果分別如圖4和圖5所示。
由圖可知,db30小波分解得到的相鄰頻帶頻率范圍均有較大重合,巴特沃斯小波分解得到的各子帶頻率范圍的重疊比db30小波分解得到的各子帶要小得多,基本上達(dá)到了對(duì)頻帶的理想二進(jìn)劃分,這對(duì)微弱故障特征提取是有利的??梢姡碗A的巴特沃斯小波頻帶分離特性明顯優(yōu)于高階的db系列小波。
圖6為采用未進(jìn)行零相位濾波處理的巴特沃斯小波變換快速算法(見式2)對(duì)圖3信號(hào)進(jìn)行3層分解的結(jié)果。
由圖6可知,各頻帶開始部分有一段值恒為零,且其長(zhǎng)度隨層數(shù)的增加而增加,這表明未進(jìn)行零相位濾波處理的巴特沃斯小波變換快速算法得到的各頻帶,均存在較大的相位延遲,且層數(shù)越大,相位延遲也越大。由圖5可知,采用零相位巴特沃斯小波變換快速算法得到各子帶不存在任何相位延遲,因此零相位巴特沃斯小波變換快速算法確實(shí)實(shí)現(xiàn)了零相位濾波,避免了巴特沃斯小波非線性相位導(dǎo)致的小波分量畸變,可以克服相位延遲給故障特征提取帶來(lái)的不利影響,較好地解決小波的正交性與對(duì)稱性不能同時(shí)滿足的矛盾。
滾動(dòng)軸承的滾動(dòng)體既有自轉(zhuǎn)又有公轉(zhuǎn),并且它的早期故障沖擊信號(hào)傳至傳感器的過(guò)程中所受的干擾比其他元件多得多,因此滾動(dòng)體早期故障特征比其他元件更加微弱,提取也更加困難[11]。
筆者就利用滾動(dòng)軸承滾動(dòng)體早期故障振動(dòng)信號(hào)來(lái)分析比較20階巴特沃斯小波與db40小波對(duì)軸承早期故障特征提取的效果。
實(shí)驗(yàn)數(shù)據(jù)來(lái)源于文獻(xiàn)[12]中軸承故障模擬實(shí)驗(yàn)臺(tái),圖7為滾動(dòng)體早期故障的時(shí)域振動(dòng)信號(hào)及其頻譜,分別用兩種小波對(duì)該信號(hào)進(jìn)行分解的結(jié)果如圖8和圖9所示。從圖中可以看出,巴特沃斯小波分解得到的子帶W3中出現(xiàn)了比較清晰的等間隔的沖擊特征,其沖擊時(shí)間發(fā)生間隔為8.8ms,即對(duì)應(yīng)于故障特征頻率114Hz。
這表明巴特沃斯小波已提取到隱藏在強(qiáng)噪聲和其他干擾背景下的滾動(dòng)體微弱故障特征。但db40小波分解得到的各個(gè)子帶均沒(méi)有出現(xiàn)這種現(xiàn)象,這就說(shuō)明巴特沃斯小波對(duì)微弱故障特征的提取能力要強(qiáng)于db40小波。
圖10為db40小波分解與巴特沃斯小波分解的各子帶的包絡(luò)譜比較。從圖中可以看出,db40小波分解的各子帶信號(hào)的包絡(luò)譜中看不見滾動(dòng)體故障特征頻率114Hz的存在,而巴特沃斯小波分解的W3、W4子帶信號(hào)的包絡(luò)譜中,可以很清晰的看見滾動(dòng)體故障特征頻率114.3Hz,這進(jìn)一步說(shuō)明了巴特沃斯小波對(duì)微弱故障特征的提取能力要強(qiáng)于db40小波。
通過(guò)以上的滾動(dòng)體早期故障特征的提取效果比較,可知巴特沃斯小波有非常強(qiáng)的微弱故障特征提取能力,非常適合于機(jī)械早期故障診斷。巴特沃斯小波對(duì)滾動(dòng)體早期故障特征的提取效果優(yōu)于db40小波的原因,除了巴特沃斯小波濾波器擁有優(yōu)良的頻帶分離特性外,更主要的是巴特沃斯小波可以根據(jù)信號(hào)的特點(diǎn)靈活地選擇尺度參數(shù)來(lái)將原信號(hào)高頻帶分解在更精細(xì)的頻帶范圍內(nèi)。在本次故障特征提取效果比較中,db40小波用尺度參數(shù)2j進(jìn)行分解,其分析精度相當(dāng)于1/1倍頻程,而巴特沃斯小波選用尺度參數(shù)2j/6進(jìn)行分解,其分析精度相當(dāng)于1/6倍頻程。對(duì)比圖8與圖9可知,db40對(duì)高頻頻帶的分解精度太低,幾乎將其全部分解在W1、W2子帶中,這樣W1、W2子帶的包絡(luò)中必定混入較多的噪聲和干擾以致它們包絡(luò)譜中的滾動(dòng)體故障特征頻率淹沒(méi)在噪聲和干擾中而無(wú)法辨識(shí);而巴特沃斯小波將軸承振動(dòng)信號(hào)的高頻頻帶進(jìn)行了相當(dāng)高精度的分解,將其分解為W1、W2、W3和W44個(gè)頻帶,這樣把含有噪聲和干擾較少而含故障沖擊信號(hào)所調(diào)制的高頻固有振動(dòng)信號(hào)較多的高頻子帶W3、W4分離出來(lái),因此W3、W4子帶的包絡(luò)譜中就明顯出現(xiàn)滾動(dòng)體的故障特征頻率。
1)設(shè)計(jì)了基于巴特沃斯數(shù)字濾波器的巴特沃斯小波及其濾波器組。該濾波器組可以通過(guò)改變其階數(shù)方便地調(diào)整其頻率截止特性,從而抑制頻帶重疊現(xiàn)象,還可通過(guò)改變尺度參數(shù)來(lái)改變對(duì)信號(hào)的頻帶劃分的方式從而將原信號(hào)高頻帶分解在更精細(xì)的頻帶范圍,因此巴特沃斯小波克服了通常正交小波幅頻特性不好,并且總按二進(jìn)方式劃分頻帶以致高頻帶頻率分析精度不足的缺陷。
2)提出了零相位巴特沃斯小波快速分解算法,成功消除分解過(guò)程中存在相位延遲,克服了非線性相位給故障特征提取帶來(lái)的不利影響。
3)通過(guò)頻帶分離特性比較實(shí)例證實(shí)了巴特沃斯小波比其他常用的一些正交小波具有更好的頻帶分離特性。通過(guò)軸承滾動(dòng)體早期故障特征提取效果的對(duì)比分析證明巴特沃斯小波比其他一些常用的正交小波有更強(qiáng)的微弱故障特征提取能力。除了擁有優(yōu)良的頻帶分離特性外,更主要的是巴特沃斯小波可以根據(jù)信號(hào)的特點(diǎn)靈活地選擇尺度參數(shù)來(lái)將原信號(hào)分解在更精細(xì)的頻帶范圍內(nèi),因此,巴特沃斯小波具有更強(qiáng)的微弱特征提取能力,可以廣泛應(yīng)用于機(jī)械早期故障診斷。
(References)
[1]張中民,盧文祥,楊叔子.基于小波系數(shù)包絡(luò)譜的滾動(dòng)軸承故障診斷[J].振動(dòng)工程學(xué)報(bào),1998,11(1):65-69.
ZHANG Zhongmin,LU Wenxiang,YANG Shuzi.Application of the envelope spectrum of wavelet coefficients to fault diagnosis in roller bearing[J].Journal of Vibration Engineering,1998,11(1):65-69.(in Chinese)
[2]ZHAO Xuezhi,YE Bangyan.Convolution wavelet packet transform and its applications to signal processing[J].Digital Signal Processing,2010,20:1352-1364.
[3]楊建國(guó).小波分析及其工程應(yīng)用[M].北京:機(jī)械工業(yè)出版社,2005:47-55.
YANG Jianguo.Wavelet analysis and its engineering applications[M].Beijing:China Machine Press,2005:47-55.(in Chinese)
[4]黃云志,徐科軍.基于IIR小波濾波器的渦街流量計(jì)數(shù)字信號(hào)處理系統(tǒng)[J].儀器儀表學(xué)報(bào),2007,28(12):2142-2146.
HUANG Yunzhi,XU Kejun.Digital signal processing system based on IIR wavelet filter for vortex flow meter[J].Chinese Journal of Scientific Instrument,2007,28(12):2142-2146.(in Chinese)
[5]趙學(xué)智,陳統(tǒng)堅(jiān),葉邦彥,等.內(nèi)積型和卷積型小波變換對(duì)信號(hào)處理效果的研究[J].機(jī)械工程學(xué)報(bào),2004,40(3):55-60.
ZHAO Xuezhi,CHEN Tongjian,YE Bangyan,et al.Research on signal processing effect of both inner product type of wavelet transformation and convolution type of wavelet transformation[J].Chinese Journal of Mechanical Engineering,2004,40(3):55-60.(in Chinese)
[6]萬(wàn)書亭,呂路勇,何玉靈.基于提升模式非抽樣小波變換的滾動(dòng)軸承故障診斷方法研究[J].振動(dòng)與沖擊,2009,28(1):170-173.
WAN Shuting,LYU Luyong,HE Yuling.Fault diagnosis method of rolling bearing based on undecimated wavelet transformation of lifting scheme[J].Journal of Vibration and Shock,2009,28(1):170-173.(in Chinese)
[7]DUAN Chendong,HE Zhengjia,JIANG Hongkai.New method for weak fault feature extraction based on second generation wavelet transform and its application[J].Chinese Journal of Mechanical Engineering,2004,17(4):543-547.
[8]何嶺松,王峻峰.巴特沃斯小波變換算法在故障診斷中的應(yīng)用[J].華中理工大學(xué)學(xué)報(bào):自然科學(xué)版,2000,28(10):31-33.
HE Lingsong,WANG Junfeng.Butterworth based on discrete wavelet transform and its application in machine fault diagnosis[J].Journal of Huazhong University of Science and Technology:Natural Science Edition,2000,28(10):31-33.(in Chinese)
[9]ZHANG Xi,Toshinori Yoshikawa.Design of two-channel IIR linear phase PR filter banks[J].Signal Processing.1999,72(3):167-175.
[10]郭瑜,秦樹人.無(wú)轉(zhuǎn)速計(jì)旋轉(zhuǎn)機(jī)械升降速振動(dòng)信號(hào)零相位階比跟蹤濾波[J].機(jī)械工程學(xué)報(bào),2004,40(3):50-54.
GUO Yu,QIN Shuren.Tacholess order tracking filtering for run-up or coast down viberation signal of rotating machinnery based on zero-phase distortion digital filtering[J].Chinese Journal of Mechanical Engineering,2004,40(3):50-54.(in Chinese)
[11]張進(jìn),馮志鵬,褚福磊.滾動(dòng)軸承故障特征的時(shí)間-小波能量譜提取方法[J].機(jī)械工程學(xué)報(bào),2011,47(17):44-49.
ZHANG Jin,F(xiàn)ENG Zhipeng,CHU Fulei.Extraction of rolling bearing fault feature based on time-wavelet energy spectrum[J].Chinese Journal of Mechanical Engineering,2011,47(17):44-49.(in Chinese)
[12]HUANG Yixiang,LIU Chengliang,Xuan F Zha,et al.A lean model for performance assessment of machinery using second generation wavelet packet transform and fisher criterion[J].Expert Systems with Applications,2010,37(5):3815-3822.
Application of Butterworth Wavelet in Bearing Incipient Fault Diagnosis
FENG Changlin1,LUO Rong2,LIU Hui3
(1.Naval Academy of Armament,Beijing 100161,China;2.Department of Weaponry Engineering,Naval University of Engineering,Wuhan 430033,Hubei,China;3.Ordnance Repairing Department,Hainan Province Military District,Ding’an 571231,Hainan,China)
The Butterworth wavelet and its filter banks based on Butterworth digital filter were designed.The fast algorithm of the Butterworth wavelet transform was established from the viewpoint of the convolution definition of the wavelet by use of the Butterworth wavelet filter banks and zero-phase filter.The effect of nonlinear phase on the decomposition result was overcome in the fast algorithm.Compared with the other orthogonal wavelets filter banks,Butterworth wavelet filter banks not only have better band separation capability,but also can change band separation mode by means of changing the scale to decompose the signal into the smaller frequency range.As a result,the filter banks overcome the inherent defect that orthogonal wavelets filter banks separate the band only by use of the dyadic mode.The instance analysis to the incipient fault data showed that,compared with the other orthogonal wavelet,the Butterworth wavelet possesses the better band separation characteristics and can extract the weaker fault feature more effectively.Thus the Butterworth wavelet can widely apply to the bearing incipient fault diagnosis.
fault diagnosis;Butterworth wavelet;envelope demodulation;bearing
TH165
A
1673-6524(2014)02-0076-07
2013-08-12;
2013-12-30
馮昌林(1983-),男,博士,工程師,主要從事艦炮彈藥武器論證與研究。E-mail:fcl_325@126.com