基于峭度與IMF能量融合特征和LS-SVM的齒輪故障診斷研究

王建國,楊云中,秦　波,劉永亮（內(nèi)蒙古科技大學(xué)機械工程學(xué)院，內(nèi)蒙古　包頭014010）

?

基于峭度與IMF能量融合特征和LS-SVM的齒輪故障診斷研究

王建國,楊云中,秦波,劉永亮
（內(nèi)蒙古科技大學(xué)機械工程學(xué)院，內(nèi)蒙古包頭014010）

摘要：針對齒輪振動信號非線性非平穩(wěn)特性，為避免傳統(tǒng)時頻方法在表征設(shè)備狀態(tài)時的不足，提出一種基于融合峭度與IMF能量特征和LS-SVM的齒輪故障診斷方法。首先，對齒輪振動信號在EMD分解；然后，提取包含主要故障信息的IMF分量的峭度特征和能量特征，組成融合特征向量；最后，將齒輪正常、齒根裂紋、斷齒3種狀態(tài)下的融合特征向量輸入到LS-SVM，通過訓(xùn)練好的LS-SVM對齒輪狀態(tài)進行分類識別。仿真實驗結(jié)果表明：該方法能準(zhǔn)確識別齒輪的工作狀態(tài)，且與BP神經(jīng)網(wǎng)絡(luò)、SVM相比，有著更高的故障識別效率，可用于齒輪信號的故障診斷。

關(guān)鍵詞：IMF分量；峭度和能量特征；最小二乘支持向量機；故障診斷

0 引　言

齒輪是用來改變轉(zhuǎn)速和傳遞動力的常用零件，在現(xiàn)代機械設(shè)備中有著廣泛應(yīng)用。由于齒輪本身工作環(huán)境比較惡劣，故極易受到損害，出現(xiàn)故障。齒輪故障信號由機械系統(tǒng)正常信號、周期沖擊性故障信號、干擾信號和噪聲信號等疊加而成，是一種典型的非平穩(wěn)信號[1]。因此，更好的提取出表征故障信息的特征成為齒輪故障診斷的關(guān)鍵。

文獻[1]引入EMD能量熵，提取能量特征輸入到支持向量機進行齒輪故障診斷；文獻[2]將峭度與小波包能量相結(jié)合，并輸入到概率神經(jīng)網(wǎng)絡(luò)模型，對齒輪早期故障進行降噪和診斷；文獻[3]提出了基于BP神經(jīng)網(wǎng)絡(luò)的齒輪故障診斷方法，其提取信號的偏斜度、峭度和陡度作為特征向量。然而上述方法在表征設(shè)備狀態(tài)時都有一定的不足。由于EMD是一種處理非平穩(wěn)信號的有效方法[4]，通過對齒輪信號的EMD分解后提取IMF分量的峭度特征和能量特征，可有效表征齒輪故障信息。相比于BP神經(jīng)網(wǎng)絡(luò)和SVM，最小二乘支持向量機（LS-SVM）在小樣本下具有結(jié)構(gòu)簡單、泛化能力強、學(xué)習(xí)速度快和分類準(zhǔn)確度高的特點[5]；因此，本文提出了基于峭度與IMF能量融合特征和LS-SVM的齒輪故障診斷方法，將提取的融合特征作為特征向量輸入LS-SVM，對齒輪工作狀態(tài)和故障類型進行分類。

1 特征向量構(gòu)建

1.1 IMF能量特征

經(jīng)驗?zāi)B(tài)分解（EMD）是Huang于1998提出的一種用于非線性、非平穩(wěn)信號分析的新方法,將信號自適應(yīng)地分解為若干內(nèi)蘊模式函數(shù)之和。EMD方法基于信號的局部時間特征進行自適應(yīng)的分解，克服了用無意義的諧波分量來表示非線性、非穩(wěn)態(tài)信號的缺點[6]，非常適合于非線性、非穩(wěn)態(tài)信號的分析。EMD分解后獲得的IMF分量必須滿足以下兩個條件：1）在整個數(shù)據(jù)序列中，極大值與極小值點個數(shù)之和與過零點個數(shù)相等，且相差不超過1個；2）在任意時間點上，由局部極大值與局部極小值點構(gòu)成的包絡(luò)均值為零或近似為零。

經(jīng)驗?zāi)B(tài)分解的具體步驟如下：

1）采用3次樣條插值擬合原始信號為x（t）的上下包絡(luò)線。得到上下包絡(luò)線的均值為m（t），包絡(luò)線均值減去原始信號得到數(shù)據(jù)序列h（t），即：

2）數(shù)據(jù)序列h（t）如果滿足上述IMF分量條件，則得到信號的第1個IMF分量，即h（t）；如不滿足上述條件，則重新對h（t）進行以上篩選過程，直到得到符合要求的IMF分量。

3）得到信號的第1個IMF分量和該分量分離后的余項為

對余項r（t）重復(fù)上述篩選過程，由此得到一系列的IMF分量ci（t）（i=1，2，…，n），最終的余項rn（t）。原始信號可以重構(gòu)為

式中：n——IMF分量個數(shù)；

rn——殘差函數(shù)，是一個單調(diào)函數(shù)。

EMD能夠?qū)⑿盘柗纸鉃槿舾蓚€具有不同頻段范圍的IMF分量，而當(dāng)故障發(fā)生時，不同故障類型會導(dǎo)致不同頻率范圍的信號發(fā)生變化，從而引起相應(yīng)頻段內(nèi)的能量發(fā)生變化。對齒輪振動信號x（t）進行EMD分解，可以計算得到各IMF分量的能量E1，E2，…，En。根據(jù)EMD分解的正交性，忽略殘余分量,可知原始信號的總能量應(yīng)該等于計算的n個IMF分量的能量之和。由于各個IMF分量包含不同的頻率成分,且具有不同的能量，因此構(gòu)成了信號的能量特征向量E=[E1，E2，…，En]。由此可以得到能量熵的定義：

其中pi=Ei/E，表示第i個本征模式分量能量在總能量中的比重。

計算齒輪正常、齒根裂紋和斷齒3種狀態(tài)下的EMD能量熵值，結(jié)果如圖1所示。由圖可知，齒輪的工作狀態(tài)和故障類型不同，得到的能量熵值也就不同；因此，可以通過能量特征作為特征向量來判斷齒輪的工作狀態(tài)和故障類型。

圖1 齒輪不同狀態(tài)下的能量熵值

1.2峭度特征

信號的峭度指標(biāo)是無量綱參數(shù)，對于振動中的沖擊信號非常敏感，適用于表面損傷類故障[7]。齒輪發(fā)生故障時，會產(chǎn)生明顯的沖擊，導(dǎo)致不同齒輪不同狀態(tài)的峭度值不同。計算齒輪正常、齒根裂紋和斷齒3種狀態(tài)下的峭度值，結(jié)果如圖2所示。由圖可知，齒輪的工作狀態(tài)和故障類型不同，其得到的峭度值也就不同；因此，可以通過峭度特征作為特征向量來判斷齒輪的工作狀態(tài)和故障類型。

圖2 齒輪不同狀態(tài)下的峭度值

1.3 構(gòu)建特征向量

首先將采集到的振動信號通過EMD分解，選出包含主要故障信息的前m個IMF分量作為研究對象。然后計算選出的IMF分量能量（i=1，2，…，m），得到能量特征向量T=[E1/E，E2/E，…，Em/E]，其中；接著計算選出的IMF分量的峭度，得到峭度特征向量K=[K1，K2，…，Km]。最后組成融合特征向量F=[E1/E，E2/E，…，Em/E，K1，K2，…，Km]。

2 基于LS-SVM的齒輪故障診斷模型的建立

2.1 最小二乘支持向量機

最小二乘支持向量機（least squares support vector machine，LS-SVM），采用最小二乘線性系統(tǒng)作為損失函數(shù)，代替?zhèn)鹘y(tǒng)SVM采用的二次規(guī)劃方法，簡化了計算過程，運算速度也明顯提高[8-9]。LS-SVM分類的具體實現(xiàn)過程如下：

假設(shè)訓(xùn)練樣本為（xi，yi），xi∈Rn，yi∈R，其中xi是輸入向量，yi是輸出向量。將輸入向量ψ（x）通過非線性映射從原始樣本空間映射到一個高維特征空間，在高維特征空間中，求取最優(yōu)分類超平面，使分類間隔最大[10]。為描述這個分類超平面，構(gòu)造線性決策函數(shù)為

式中：ω——權(quán)向量；

b——常數(shù)。

根據(jù)結(jié)構(gòu)風(fēng)險最小化的原則，將分類問題轉(zhuǎn)化為如下優(yōu)化問題

式中：γ——懲罰因子；

ε——松弛因子。

定義拉格朗日函數(shù)：

其中，拉格朗日乘子αi∈R，根據(jù)KKT條件，對上式進行優(yōu)化，即對ω，b，ε，α的偏導(dǎo)數(shù)為0，得：

通過求解式（9），用最小二乘法求出α和b，可得到分類決策函數(shù)：

其中K（x，xi）為核函數(shù)，滿足Mercer條件，對應(yīng)于特征空間的點積。本文選用性能較好的徑向基核函數(shù)，其表達式為

2.2 基于LS-SVM的診斷模型構(gòu)建

將齒輪振動信號進行EMD分解得到IMF分量，提取其融合特征作為特征向量輸入到最小二乘支持向量機進行模式識別，從而判斷齒輪的狀態(tài)[11]。具體流程如圖3所示。

其具體步驟如下：

1）采集齒輪正常、齒根裂紋、斷齒故障的振動信號若干組。

2）通過EMD分解，提取其融合特征向量。

3）選用“一對一”算法的多分類支持向量機，建立3個最小支持向量機（LS-SVM1，LS-SVM2，LSSVM3）組成多故障分類器，將其特征向量作為輸入，對LS-SVM進行訓(xùn)練和分類測試[12]。

圖3 基于LS-SVM的齒輪故障診斷模型流程圖

圖4 試驗裝置結(jié)構(gòu)圖

圖5 齒輪箱結(jié)構(gòu)簡圖

3 仿真實驗

本文的試驗裝置結(jié)構(gòu)如圖4所示，其主要由電動機、聯(lián)軸器、行星齒輪箱、平行軸齒輪箱、可編程磁力控制器等組成。采用ZonicBook/618E振動測量系統(tǒng)采集數(shù)據(jù)，3個加速度傳感器分別安裝在平行軸齒輪箱中間軸左端軸承蓋的垂直徑向、水平徑向和軸向上。平行軸齒輪箱的結(jié)構(gòu)如圖5所示，齒輪箱中有兩對齒輪副，4個齒輪均為直齒輪，齒數(shù)分別為Z1= 100，Z2=29，Z3=90，Z4=36。

表1 數(shù)據(jù)集

試驗中，電機轉(zhuǎn)速設(shè)置為2100r/min，采樣頻率為5120 Hz，采樣時間為5s。齒輪正常、齒根裂紋和斷齒故障3種狀態(tài)下，每種狀態(tài)各測試20組。各組中任意的10組數(shù)據(jù)作為訓(xùn)練樣本，剩余10組數(shù)據(jù)作為測試樣本。數(shù)據(jù)集的詳細描述見表1。

表2 齒輪不同狀態(tài)下的部分特征向量

表3 LS-SVM測試結(jié)果

表4 LS-SVM與SVM、BP神經(jīng)網(wǎng)絡(luò)性能比較

首先對訓(xùn)練和測試樣本數(shù)據(jù)進行EMD分解得到其IMF分量，選用了前6個包含主要故障信息的IMF分量，求取3種狀態(tài)信號的IMF分量的峭度特征和能量特征，組成一個10×12的特征向量矩陣作為LS-SVM的輸入。表2為訓(xùn)練樣本中齒輪不同狀態(tài)下的部分特征向量。將訓(xùn)練樣本特征向量輸入到由3個LS-SVM組成的多分類器中進行訓(xùn)練，然后將10組測試樣本輸入到分類器中進行識別，為區(qū)分有無故障，定義輸出“+1”表示正常狀態(tài)，“-1”表示故障狀態(tài)。區(qū)分完有無故障后，為識別具體的故障類型，定義輸出“+1”表示齒根裂紋故障，“-1”表示斷齒故障，測試結(jié)果見表3。從表中可看出，LS-SVM能夠?qū)y試樣本進行正確率很高的故障診斷。說明本文提出的齒輪故障診斷方法有效。

試驗選取了10組原始信號進行處理，屬于小樣本情況。在小樣本情況下，LS-SVM具有分類的訓(xùn)練時間短，收斂速度快以及測試準(zhǔn)確度高的特點，以其作為分類器，性能比BP神經(jīng)網(wǎng)絡(luò)、SVM分類器要好。表4通過實驗對3種分類器進行了性能比較，從表中可看出，LS-SVM在小樣本情況下具有良好的分類準(zhǔn)確度，能準(zhǔn)確進行齒輪的故障診斷。

4 結(jié)束語

本文通過對齒輪振動信號進行EMD分解，提取包含主要故障信息的IMF分量峭度特征和能量特征作為融合特征向量輸入LS-SVM，通過訓(xùn)練好的LSSVM對齒輪狀態(tài)進行分類識別。實驗結(jié)果表明，該方法能準(zhǔn)確識別齒輪的工作狀態(tài)，且與BP神經(jīng)網(wǎng)絡(luò)、SVM相比，有著更高的故障識別效率，可用于齒輪信號的故障診斷。

參考文獻

[1]張超，陳建軍，郭迅.基于EMD能量熵和支持向量機的齒輪故障診斷方法[J].振動與沖擊，2010，29（10）：216-220.

[2]陳水宣，鄒俊，易際明，等.基于峭度和小波包能量特征的齒輪箱早期故障診斷及抗噪研究[J].機械傳動，2012，36（11）：9-14.

[3]吳偉，鄭娟.基于BP神經(jīng)網(wǎng)絡(luò)的齒輪故障診斷[J].機械研究與應(yīng)用，2015，28（1）：88-93.

[4]于德介，程軍圣，楊宇. Hilbert-Huang變換在齒輪故障診斷中的應(yīng)用[J].機械工程學(xué)報，2005，41（6）：102-106.

[5]肖志勇，楊小玲，劉愛倫.基于改進EMD和LSSVM的機械故障診斷[J].自動化儀表，2008，29（6）：24-26.

[6]唐貴基，龐爾軍，王曉龍.基于EMD的齒輪箱齒輪故障診斷的研究[J].機床與液壓，2013，41（13）：188-190.

[7]韓清凱，于曉光.基于振動分析的現(xiàn)代機械故障診斷原理及應(yīng)用[M].北京：科學(xué)出版社，2010：49-51.

[8]黃俊，潘宏俠，都衡.基于EMD近似熵和LSSVM的齒輪箱故障診斷研究[J].組合機床與自動化加工技術(shù)，2014（3）：111-113.

[9]韓永杰，郝偉，董辛?xí)F.基于EMD和沖擊脈沖的齒輪故障診斷[J].機械設(shè)計與制造，2009（10）：222-224.

[10] ABBASIONAS，RAFSANJANIA A，F(xiàn)ARSHIDIANFARB A，et al. Rolling element bearings multi-fault classification based on the wavelet denoising and support vector machine[J]. Mechanical Systems and Signal Processing，2007（21）：2933-2945.

[11] LI Y B，XU M Q，WEI Y，et al. An improvement EMD method based on the optimized rational Hermite interpolation approach and its application to gear fault diagnosis[J]. Measurement，2015（63）：330-345.

[12] YANG D L，LIU Y L，LI S B，et al. Gear fault diagnosis based on support vector machine optimized by artificial bee colony algorithm [J]. Mechanism and Machine Theory，2015，90（8）：219-229.

（編輯：劉楊）

Gear fault diagnosis research based on kurtosis and IMF energy feature fusion and least squares support vector machine

WANG Jianguo，YANG Yunzhong，QIN Bo，LIU Yongliang
（Mechanical Engineering School，Inner Mongolia University of Science & Technology，Baotou 014010，China）

Abstract:Gear vibration signals have nonlinear and non-stationary characteristics. To avoid the disadvantages of existing time and frequency domain methods in the characterization of equipment state，this paper has been proposed a gearbox fault diagnosis method based on kurtosis and IMF energy feature fusion and least squares support vector machine. First，the gear vibration signals were decomposed by the EMD method. Second，the IMF components which contain major fault information were extracted and their energy and kurtosis feature calculated as fusion vectors. Third，the fusion feature vectors of three teeth conditions，viz.，normal，root crack and broken，were input to the least squares support vector machine（LS-SVM）to classify and identify gearbox faults. The simulation results show that this method can accurately identify the gear working state and more efficient in fault identification compared with BP neural network and SVM. It can be used for diagnosing gear signal faults.

Keywords:IMF component；kurtosis and energy feature；LS-SVM；fault diagnosis

文獻標(biāo)志碼：A

文章編號：1674-5124（2016）04-0093-05

doi：10.11857/j.issn.1674-5124.2016.04.020

收稿日期：2015-06-01；收到修改稿日期：2015-07-20

基金項目：國家自然科學(xué)基金項目（21366017）；內(nèi)蒙古科技廳高新技術(shù)領(lǐng)域科技計劃重大項目（20130302）

作者簡介：王建國（1958-），男，內(nèi)蒙古呼和浩特市人，教授，碩士生導(dǎo)師，博士，研究方向為機電系統(tǒng)智能診斷與復(fù)雜工業(yè)過程建模的優(yōu)化。