張 敏 ，蔡振宇 ，包珊珊

（1.西南交通大學(xué)機械工程學(xué)院，四川 成都 610031；2.西南交通大學(xué)軌道交通運維技術(shù)與裝備四川省重點實驗室，四川 成都 610031）

滾動軸承廣泛存在于機械結(jié)構(gòu)中，作為重要的旋轉(zhuǎn)結(jié)構(gòu)，是最容易出現(xiàn)故障的部件之一.滾動軸承振動信號具有非線性與非穩(wěn)定的特性，單從時域或頻域進行故障診斷較困難[1].于德介等[2]首次引入固有模態(tài)函數(shù)（intrinsic mode function，IMF）分量應(yīng)用于機械故障診斷，本文簡稱模態(tài)函數(shù)，將M 距離函數(shù)和支持向量機（support vector machine，SVM）進行結(jié)合實現(xiàn)故障模式識別，但經(jīng)驗?zāi)B(tài)分解（empirical mode decomposition，EMD）[3]分解出具有原始信號特征的IMFS 容易出現(xiàn)模態(tài)混疊，掩蓋原始信號的真實特征[4-5].Wu 和Huang 等[6]在EMD 的基礎(chǔ)上改進提出集合經(jīng)驗?zāi)B(tài)分解（ensemble empirical mode decomposition，EEMD），引入高斯白噪聲能夠有效地避免模態(tài)混疊[7]，分解出的IMFS 具有完整信息，但IMF 所含信息有限，還需進一步處理.Hilbert 包絡(luò)解調(diào)是一種運用廣泛的信息解調(diào)技術(shù)，但在確定帶通濾波參數(shù)時需要豐富的經(jīng)驗，限制了該技術(shù)的發(fā)展[8]，若直接在模態(tài)函數(shù)上進行Hilbert 變換則可避免帶通濾波參數(shù)的確定.

多分類SVM 通過核參數(shù)將特征向量映射到高維空間實現(xiàn)分類，參數(shù)的確定影響整個分類效果.何青等[9]用果蠅優(yōu)化算法優(yōu)化SVM 參數(shù)，但測試集量少，大數(shù)據(jù)下還有待驗證.煙花算法（firework algorithm，F(xiàn)WA）[10]是Tan 和Zhu 在2010年提出的一種新型進化算法，具有很強的優(yōu)化求解能力，應(yīng)用在優(yōu)化SVM參數(shù)時比粒子群算法（particle swarm optimization，PSO）迭代時間快且更準確[11]，近年來逐漸受到研究者的關(guān)注.

本文提出利用EEMD 將原始數(shù)據(jù)分解成IMFS，再對IMF 進行Hilbert 變換，避免帶通濾波參數(shù)的確定，同時可對模態(tài)函數(shù)進行信息解調(diào)分析；將模態(tài)函數(shù)及其變換后的解調(diào)信息進行統(tǒng)計特征提取并降維處理；最后采用新型煙花算法優(yōu)化SVM 分類參數(shù)，代入數(shù)據(jù)實現(xiàn)快速、準確的故障診斷.

1 數(shù)據(jù)特征提取

1.1 集合經(jīng)驗?zāi)B(tài)分解

EEMD 是通過給目標信號加入一定幅值的高斯白噪聲，利用白噪聲頻譜均衡分布的特點來均衡噪聲，有效地解決了EMD 出現(xiàn)模態(tài)混疊的現(xiàn)象.適用于各種非線性和非穩(wěn)定的信號處理，步驟如下：

步驟1隨機生成均值 μ為0，標準差 σ為en的高斯白噪聲nm(t)，t為時間，設(shè)定原始參數(shù)，m表示第m次分解，1 ≤m≤M，M表示EMD 分解次數(shù).

步驟2將高斯白噪聲nm(t)加入待處理的信號y(t)中 ，得到信號為ym(t)，t為時間，即

步驟3對處理后的信號ym(t)進行EMD 分解，得到S個IMF 分量cs,m(t)，

式中：rm(t)為 第m次分解得到的余項；cs,m(t)為進行第m次EMD 分解后的第s個IMF 分量，由頻率從高到低排列.

步驟4若m＜M，則返回執(zhí)行步驟（2），M=m+1，到M次終止；

步驟5計算進行M次EMD 分解后得到每個IMF 分量的均值，根據(jù)不相關(guān)隨機序列統(tǒng)計均值為0，消除加入高斯白噪聲對真實IMF 分量的影響，即：

步驟6將(t)作為EEMD 最終得到的第s個IMF 分量.

1.2 Hilbert 包絡(luò)解調(diào)原理

信號進行EEMD 分解得到排序降次的IMF，本文選取前幾個作為研究對象，后面殘存的以低頻噪聲為主，不予考慮.對選取的(t)進行Hilbert 變換，解調(diào)原理如下：

濾波器公式h(t)為

式中：δ(t)為脈沖函數(shù).

解析信號公式q(t)為

1.3 統(tǒng)計特征提取與降維

將主要的IMF 及其對應(yīng)的解調(diào)信息作為提取特征對象，為有效提取特征值，挑選最合適的統(tǒng)計特征值來對處理的信息進行特征提取，本文選取的統(tǒng)計特征值依次是均值、峰值、極差、標準差、偏度、峰度、變異系數(shù)和平方和.綜合上述的所有統(tǒng)計特征量，特征維度較高，本文采用核主元分析（kernel principal component analysis，KPCA）對特征維度進行壓縮提取，實現(xiàn)特征降維，便于后面快速故障診斷.由于篇幅問題，本文不對KPCA 降維具體內(nèi)容進行論述，詳細參考文獻[12].

2 模式識別

2.1 支持向量機

支持向量機通過核函數(shù)來實現(xiàn)線性不可分向線性可分的轉(zhuǎn)化，研究表明徑向基核函數(shù)K(xi,yi)在SVM 中表現(xiàn)出良好的泛化能力[13]，將輸入向量從原來的空間映射到高維特征空間P，并在該特征空間P內(nèi)建立優(yōu)化超平面.分類線方程為[14-15]

式中：(xi,yi)為 訓(xùn)練樣本；ω為權(quán)值；x為輸入向量值；b為閾值；l為向量的個數(shù).

根據(jù)Karush-Kuhn-Tucker，優(yōu)化各個系數(shù)得到最優(yōu)決策函數(shù)為

式中：ai為Lagrange 系數(shù).

最后得到適應(yīng)度函數(shù)為

式中：C為懲罰參數(shù)；σ為核參數(shù).

由式（10）可知進行SVM 分類模型構(gòu)建時，性能的關(guān)鍵因素在于參數(shù)C和 σ的選取.

2.2 FWA-SVM 模型

FWA 是將每個煙花都當作一個解空間中的可行解，通過爆炸產(chǎn)生煙花點作為全局搜索的可行解.通過每個煙花的適應(yīng)度值來變化爆炸半徑和爆炸數(shù)，適應(yīng)度值越小的點爆炸范圍越小，爆炸數(shù)越多，適應(yīng)度值大的則相反.煙花算法核心包括爆炸算子、變異操作、映射規(guī)則、選擇策略4 個部分[16].

假定待求解的優(yōu)化問題形式為minf(u)∈R,u∈Ω，Ω為可行解.煙花算法對SVM 參數(shù)優(yōu)化步驟如下：

步驟1在解空間內(nèi)隨機生成N個初始un，有N個煙花；

步驟2計算每個煙花的適應(yīng)度值與它們的爆炸半徑Bn和爆炸花火個數(shù)Qn；

式中：ymin=min{f(un)},為這次迭代中的最優(yōu)值，也為最小值，n=1,2,...,N；ymax=max{f(un)}為當前迭代中的最劣值，也為最大值；H為爆炸火花數(shù)的大??；ε為機器最小量，避免出現(xiàn)零操作.

為了避免適應(yīng)度值優(yōu)的或者劣的產(chǎn)生過多或者過少的爆炸火花，文獻[10]對火花個數(shù)做出如下的限制：

式中：Sn為火花個數(shù)；a、b為常數(shù).

步驟3產(chǎn)生爆炸火花，集合DC具有z個維度，z=round(D×rand(0,1))，D為煙花un的維數(shù).在DC中的每個維度k下進行爆炸操作，再經(jīng)過越界處理將Tnk保存到火花種群中.

式中：h為偏移量；unk為第n個煙花在第k維上的位置；Tnk為unk爆炸操作后的火花位置.

步驟4進行高斯變異操作，每個維度通過式（16）進行高斯變異，再經(jīng)過越界處理保存到高斯種群當中.

式中：e～N(1,1)表示均值為1，方差為1 的高斯分布.

步驟5選擇操作，在所有得到的種群中挑選最好的一個，另外N-1 個則通過輪盤賭法進行選擇如式（17）、（18）.

式中：P(un)為R(un)的 概率；R(un)為個體un與其它個體的距離之和；d(un-uj)是 指任意兩個個體un和uj之間的歐式距離.

步驟6判斷是否滿足終止迭代條件，滿足則結(jié)束輸出最優(yōu)SVM 參數(shù)，不滿足就繼續(xù)迭代.

2.3 基于EEMD-Hilbert 和FWA-SVM 的滾動軸承故障診斷

基于EEMD-Hilbert 特征提取和FWA-SVM 故障診斷具體步驟如下：

步驟1原始信號樣本提取.獲得設(shè)備運行所采集到的數(shù)據(jù)樣本，對樣本進行分組.

步驟2EEMD 分解處理.將樣本數(shù)據(jù)進行EEMD 分解得到若干個IMF 函數(shù).

步驟3Hilbert 變換處理.對選取的IMF 進行Hilbert 變換，獲得瞬時頻率.

步驟4統(tǒng)計特征提取.對IMF 分量和解調(diào)信息分別計算其對應(yīng)的統(tǒng)計特征值.

步驟5KPCA 特征降維.提取出有用的信息特征實現(xiàn)特征降維.

步驟6FWA 尋SVM 最佳參數(shù).將訓(xùn)練樣本帶入FWA-SVM 進行訓(xùn)練，得到最佳的SVM 分類參數(shù).

步驟7模式識別.將訓(xùn)練好的FWA-SVM 模型進行測試集故障分類，輸出結(jié)果.

基于EEMD-Hilbert 特征提取和FWA-SVM 故障診斷流程如圖1所示.

圖1 故障診斷流程Fig.1 Flowchart for troubleshooting

3 案例分析

為驗證算法模型的可行性和有效性，本文采用美國Case Western Reserve University 電氣工程實驗室的滾動軸承實驗平臺數(shù)據(jù)[17]，實驗中測試的軸承是由SFK 公司生產(chǎn)的6205-2RS 深溝球軸承.

3.1 數(shù)據(jù)特征樣本

選取滾動軸承在轉(zhuǎn)速為1797 r/min，采樣頻率為12 kHz 情況下的正常工作狀態(tài)，內(nèi)圈故障、滾動體故障和外圈故障3 種故障狀態(tài)，其中外圈故障選擇發(fā)生在6 點鐘方向上.同時考慮3 種故障狀態(tài)下的3 種損傷尺寸，分別為0.178、0.356、0.533 mm，具體數(shù)據(jù)如表1.

表1 軸承故障樣本Tab.1 Bearing failure samples

將信號進行EEMD 分解，得到正常、內(nèi)圈故障、滾動體故障和外圈故障，故障點直徑為0.178 mm的4 個IMF 分量與其對應(yīng)的原始信號如圖2所示.

3.2 各種方法結(jié)果對比

3.2.1 不同信號處理方法的比較研究

方法a 將信號進行EEMD 和Hilbert 變換，通過統(tǒng)計特征量提取，再進行KPCA 降維，將特征數(shù)據(jù)代入PSO 優(yōu)化SVM 參數(shù)模型內(nèi)進行分類，記為EEMD_H模型.方法b 與方法a 類似，但不進行Hilbert 變換，記為EEMD 模型.方法c 對原始信號進行EMD 分解[18]，后面和方法a 一樣，記為EMD_H 模型.PSO 初始參數(shù)為20 種群數(shù)，400 的最大迭代數(shù).結(jié)果如表2（數(shù)據(jù)結(jié)果為5 次平均值）所示，圖3是同一信號分別進行EEMD 和EMD 分解之后選擇首個IMF 進行圖像化的波形圖.

進行EEMD 分解的迭代時間比進行EMD 分解時間更短，且正確率要低4%左右.圖3說明EEMD可以較好避免EMD 的模態(tài)混疊現(xiàn)象.不加Hilbert變換處理的數(shù)據(jù)比加Hilbert 變換數(shù)據(jù)要差5%，說明加Hilbert 變換特征更全面，達到更好的數(shù)據(jù)特征提取.以上證明了EEMD_H 方法能夠更精確地提取原始信號的特征信號，達到高效、準確的特征提取，證明了該特征提取方法的有效性與可行性.

3.2.2 算法參數(shù)尋優(yōu)收斂性對比研究

圖2 原始信號和 IMF1～IMF4 分量波形圖Fig.2 Original signal and IMF1- IMF4 component waveforms

表2 3 種分類結(jié)果Tab.2 Classification results for the three methods

利用目前成熟的粒子群算法和遺傳算法（genetic algorithm，GA）優(yōu)化SVM 參數(shù)構(gòu)造的模型與煙花算法優(yōu)化參數(shù)進行對比，所選的粒子群算法和遺傳初始參數(shù)種群數(shù)為20，最大迭代數(shù)為200；煙花算法初始參數(shù)煙花種群數(shù)量為20，最大迭代數(shù)為200.結(jié)果如表3和圖4所示，其中適應(yīng)度值是原始算法中適應(yīng)度實值的絕對值.

由表3與圖4可知，F(xiàn)WA 優(yōu)化SVM 參數(shù)比PSO更早達到最大值，且FWA 比PSO、GA 得到的訓(xùn)練集正確率更高，說明煙花算法比PSO、GA 在搜索域內(nèi)能快速、準確得到最優(yōu)的核函數(shù)參數(shù).GA 算法源于搜索域較廣，容易陷入局部最優(yōu)，間接說明FWA的迭代收斂能力更強.另外FWA 與PSO 出現(xiàn)此類差別原因需從算法構(gòu)造思路不同處進行分析：（1）煙花算法通過分布式信息共享，根據(jù)分布在不同區(qū)域煙花的適應(yīng)度值決定爆炸強度大小和輻射范圍，但PSO 是單項流動，搜索迭代過程是跟隨當代最優(yōu)解；（2）煙花算法中的高斯變異在變異中選出的不同維，而維度上的位移是相同的，保證了某些維度之間一些聯(lián)系，而PSO 中各維變異是不相同的；（3）煙花算法中的高斯變異每代都要進行，而粒子群中每隔一定的迭代次數(shù)才運行1 次.

3.2.3 不同模型的故障診斷效果對比研究

利用遺傳算法GA 和粒子群算法PSO 分別優(yōu)化SVM 參數(shù)構(gòu)造的模型與煙花算法FWA 優(yōu)化SVM 參數(shù)進行對比，分別記為GA 模型、PSO 模型和FWA 模型.各算法初始種群數(shù)都為20，最大迭代數(shù)為400，其他都選各自最合適的參數(shù).將上述處理好的訓(xùn)練集和測試集代入構(gòu)造好的模型，進行結(jié)果分析.結(jié)果如表4（數(shù)據(jù)結(jié)果為5 次平均值）和圖5所示.

從表4和圖5可以看出，在時間上FWA 優(yōu)化SVM 參數(shù)的平均時間只要14 s，相對較好迭代時間效果GA 縮短了100 s 多，而PSO 迭代時間較長，效果較差；FWA 優(yōu)化的模型也在正確率上體現(xiàn)出優(yōu)勢，比GA 高0.4%，比PSO 高0.2%.以上說明FWA能夠?qū)崿F(xiàn)快速、準確地對SVM 參數(shù)優(yōu)化，因此驗證了煙花算法能夠很好地優(yōu)化SVM 參數(shù)，構(gòu)建的模式識別模型能準確地進行故障診斷.

圖3 兩種方法首個IMFFig.3 Two methods for the first IMF chart

表3 FWA、PSO、GA 對SVM 參數(shù)尋優(yōu)Tab.3 FWA，PSO，and GA for SVM parameter optimisation

圖4 3 種算法SVM 參數(shù)迭代對比Fig.4 Comparison of SVM parameters for the three algorithms

表4 3 種分類結(jié)果Tab.4 Classification results for the three algorithms

圖5 3 種分類結(jié)果對比Fig.5 Comparison of the three classification results

4 結(jié)束語

提出一種采用集合經(jīng)驗?zāi)B(tài)分解、Hilbert 變換的特征提取方法，并利用煙花算法優(yōu)化支持向量機分類參數(shù)的滾動軸承故障診斷方法.將振動信號進行EEMD 分解，避免了EMD 分解的模態(tài)混疊現(xiàn)象，準確提取信號特征；通過對IMF 分量進行Hilbert 變換獲得頻域統(tǒng)計特征，綜合了Hilbert 處理數(shù)據(jù)的優(yōu)點，避免帶通濾波參數(shù)的選取，實現(xiàn)原始信號特征更準確提取；最后利用煙花算法優(yōu)化SVM 參數(shù)，比傳統(tǒng)的遺傳算法、粒子群算法識別率更高，迭代時間更少，實現(xiàn)了快速、準確的故障診斷.案例分析的結(jié)果證明，運用該算法模型在故障診斷方面是一條可行的途徑.