1.3 仿真信號分析
設(shè)置仿真信號x(t)由4個分量疊加,如下式所示
式中 x1(t)為余弦信號,x2(t)和x3(t)為調(diào)頻信號,x4(t)為高斯白噪聲。
分別采用MEWT和原始EWT對x(t)進(jìn)行分解,并取與x(t)相關(guān)性較強(qiáng)的前3層進(jìn)行重構(gòu)分析,其時頻譜分別如圖1和2所示。由圖可見,與原始EWT相比,MEWT能較為準(zhǔn)確地分解仿真信號,對噪聲魯棒性較強(qiáng)。
1.4 綜合評價指標(biāo)
文獻(xiàn)[12]提出一種加權(quán)評價指標(biāo),有效改進(jìn)了VMD的虛假模態(tài)問題,本文將其引入MEWT中,其表達(dá)式如下:
式中 φ,η,ζ為比例系數(shù),0<φ,η,ζ<1。Kr表示峭度,其值越大,表明信號周期性沖擊成分比重越大;rxy為相關(guān)系數(shù),其值越大,表明各模態(tài)與原始信號相關(guān)程度越大;exy為能量比,可剔除MEWT因能量泄漏引入的附加成分??紤]到軸承振動信號的特性,經(jīng)反復(fù)試驗并結(jié)合文獻(xiàn)[12]和粒子群優(yōu)化算法,最終取各指標(biāo)為相同權(quán)值,并選擇P值最大的前4個模態(tài)分量進(jìn)行重構(gòu)。
2 IFractalNet
2.1 FractalNet簡介
CNN作為深度學(xué)習(xí)的重要模型,主要由卷積層、池化層、全連接層和分類層等組成,一般架構(gòu)如圖3所示。
CNN隨層數(shù)加深,梯度消失現(xiàn)象明顯,進(jìn)而造成識別率降低,且其最優(yōu)結(jié)構(gòu)難以確定。FractalNet[13]不依賴殘差框架,通過選擇合適的子路徑集合提升模型表現(xiàn),較好地解決了CNN存在的問題,結(jié)構(gòu)如圖4所示。
其中,紅色卷積層Convolution為基礎(chǔ)層;綠色Join層為相加操作。fC(z)中,C為列數(shù),z為輸入,C=1表示一個基礎(chǔ)層。fC+1(z)如圖4(b)所示,在右邊疊加兩個fC(z),左邊連接一個基礎(chǔ)層。以此類推,當(dāng)C取4時,f4(z)可作為一個block塊。圖4(d)所示為網(wǎng)絡(luò)連接5個block,block之間以池化層連接,最后是預(yù)測層。令block個數(shù)為B,每個block列數(shù)為C,則FractalNet深度為B×2C-1。FractalNet采取drop?path正則化,如圖5所示,可有效防止網(wǎng)絡(luò)過擬合。網(wǎng)絡(luò)訓(xùn)練時,mini?batch之間交叉使用局部采樣和全局采樣,圖5(a)和(c)為局部采樣方法,對join層的輸入進(jìn)行dropout操作,要求至少保證要有一個輸入,圖5(b)和(d)為全局采樣方法,對于整個網(wǎng)絡(luò)來說,只選擇一條路徑,且限制為某個單獨(dú)列,F(xiàn)ractalNet詳細(xì)訓(xùn)練算法見文獻(xiàn)[13]。
2.2 改進(jìn)損失函數(shù)
為有效緩解因軸承故障樣本不平衡造成故障識別率低的問題,鑒于文獻(xiàn)[14],本文改進(jìn)原始FractalNet的損失函數(shù),對樣本較多的故障類別賦以較小權(quán)重,對樣本較少的類型賦以較大的權(quán)值,如下式所示:
式(17)?(22)中 nc為第c類故障樣本個數(shù),vc為第c類故障樣本相對于其他類樣本的不平衡程度,C和Q分別為總故障類別數(shù)和樣本總數(shù),Pc(q)為樣本輸出值,wcK,l為連接權(quán)值,yq為樣本y在q類下的類標(biāo),L為網(wǎng)絡(luò)層數(shù),λ為正則項系數(shù),R(W)為偏差權(quán)重矩陣,為非卷積層之間的權(quán)值矩陣。
文獻(xiàn)[15]利用Fisher懲罰項,使深層自編碼網(wǎng)絡(luò)學(xué)習(xí)到的隱層特征對振動信號的結(jié)構(gòu)變化具有可區(qū)分性。本文將Fisher懲罰項引入FractalNet。Fisher準(zhǔn)則表達(dá)式如下
式中 JB為類間離散度,JW為類內(nèi)離散度,SB為類間離散度矩陣,SW為類內(nèi)離散度矩陣,w為投影方向。本文將FractalNet看作是從振動信號到相似矩陣的一種非線性映射關(guān)系,對于軸承振動信號樣本集X={x1,x2,,xN},模型輸出如下
2.3 改進(jìn)激活函數(shù)
CNN中最常用的ReLU激活函數(shù)雖然計算效率較高,但很容易產(chǎn)生“神經(jīng)元死亡”現(xiàn)象,造成網(wǎng)絡(luò)性能下降;PReLU,LReLU等激活函數(shù)一定程度改進(jìn)了ReLU的缺陷,但在零點(diǎn)處不光滑且效率不高。高斯誤差線性單元(Gaussian Error Linear Unit,GELU)通過對網(wǎng)絡(luò)的輸入進(jìn)行隨機(jī)正則化。與ReLU不同,GELU依據(jù)輸入的幅值大小進(jìn)行加權(quán),在語音識別任務(wù)上效果明顯優(yōu)于ReLU,GELU近似公式如下
其在正值區(qū)間內(nèi)是非線性的,并且在所有點(diǎn)處都有曲率,而ReLU在正值區(qū)間是線性的,因此缺乏曲率;GELU在正值區(qū)間的變化更為顯著,在反向傳播時可以更有效地更新梯度,在負(fù)值區(qū)域,GELU函數(shù)的值普遍更接近x軸,因此具有更佳的單邊抑制效果,可更好地緩解深層網(wǎng)絡(luò)的梯度消失或梯度爆炸現(xiàn)象。因此本文將其引入FractalNet。綜上,本文整體識別算法步驟如下:
(1)采集不同故障的滾動軸承振動信號;
(2)進(jìn)行樣本劃分,并對每個樣本進(jìn)行MEWT分解,利用加權(quán)評價指標(biāo)選擇故障特征信息明顯的imfs分量并重構(gòu);
(3)將訓(xùn)練樣本輸入IFractalNet進(jìn)行訓(xùn)練;
(4)通過測試樣本驗證基于MEWT和IFractalNet模型的識別性能。圖6為識別流程。
3 實驗驗證
3.1 實驗數(shù)據(jù)描述
為驗證提出算法的可行性和有效性,進(jìn)行實際滾動軸承故障識別實驗。本次實驗所用的測試平臺由Spectra Quest公司設(shè)計生產(chǎn),該平臺可以開展各類滾動軸承或滑動軸承的故障診斷及加速退化實驗,如圖7所示。通過杭州億恒科技有限公司的MI6008型數(shù)據(jù)采集儀、美國PCB公司的352C33型ICP加速度傳感器和筆記本電腦采集滾動軸承的振動信號。測試平臺由交流感應(yīng)電動機(jī)、電機(jī)速度控制器、測試軸承、加載系統(tǒng)等組成,該平臺是為在不同工況(即不同負(fù)載和轉(zhuǎn)速)下對軸承進(jìn)行故障識別試驗而設(shè)計的。負(fù)載由加載系統(tǒng)產(chǎn)生,應(yīng)用于被測軸承殼體,轉(zhuǎn)速由交流感應(yīng)電動機(jī)的速度控制器設(shè)定和保持。測試軸承為LDK UER204滾動軸承。
為了采集被測軸承的振動信號,兩個ICP 352C33型加速度計位于被測軸承外殼上的90°位置,即一個安裝在水平軸上,另一個安裝在垂直軸上。采樣頻率設(shè)置為25.6 kHz,每次采樣總共記錄40000個數(shù)據(jù)點(diǎn)。使用電火花技術(shù)在滾動軸承內(nèi)圈、外圈和滾動體上加工故障直徑為0.16和0.3 mm的環(huán)槽,共設(shè)計了7種不同的軸承故障工況,如表1所示。
訓(xùn)練IFractalNet需要大量樣本數(shù)據(jù),因此本文采用文獻(xiàn)[16]提出的方法對振動信號進(jìn)行有重疊樣本分割,從而實現(xiàn)樣本擴(kuò)充。最后得到每種工況下8000個樣本,隨機(jī)選取70%作為訓(xùn)練樣本,剩余作為測試樣本。圖8為7種工況的時域圖和頻域圖。
由時域圖可知,復(fù)合故障沖擊成分不明顯,部分沖擊淹沒在噪聲中,振動情況較為復(fù)雜;由頻域圖可知對于同種故障類型但不同工況的故障,從頻譜上難以區(qū)分。
3.2 軸承實際振動信號分解
以轉(zhuǎn)速30 Hz、負(fù)載12 kN工況下軸承外圈故障振動信號為例,分別采用原始EWT和MEWT對其進(jìn)行分解,結(jié)果如圖9和10所示,原始EWT模態(tài)混疊嚴(yán)重,干擾較多;MEWT分解模態(tài)數(shù)少于原始EWT,根據(jù)加權(quán)評價指標(biāo),對重構(gòu)信號進(jìn)行時頻分析,如圖11和12所示。由下式計算求得軸承外圈故障特征頻率約為120 Hz。
式中 d為滾子直徑,D為節(jié)圓直徑,為接觸角,Z為滾子數(shù),fr為轉(zhuǎn)頻30 Hz,從MEWT時頻譜中可以比較清晰地看出故障外圈故障頻率以及倍頻,驗證了MEWT的有效性。
3.3 故障識別結(jié)果與分析
采用深度置信網(wǎng)絡(luò)(deep belief network,DBN)、深層稀疏自動編碼器(deep sparse auto?encoder,DSAE)、原始EWT+IFractalNet,MEWT+CNN和MEWT+原始FractalNet等方法進(jìn)行對比分析。其中,DBN和DSAE有兩種輸入:一種為原始振動信號,一種為經(jīng)MEWT分解并重構(gòu)后的振動信號(2000維),各方法的參數(shù)如下:
方法1(本文方法):MEWT+IFractalNet,IFractalNet結(jié)構(gòu)如圖4(d)所示,卷積核寬度為3,高度為1;輸入特征圖數(shù)量為32;池化層寬度為2;卷積核移動步長為1;為避免尺寸的變化,在輸入矩陣的邊界上加入全0填充,卷積層與池化層之間都存在批歸一化層與激活層。
方法2(DBN):DBN結(jié)構(gòu)為2000?1000?500?200?50?7,每個限制玻爾茲曼機(jī)的學(xué)習(xí)率、動量因子分別為0.04,0.1。
方法3(MEWT+DBN):DBN結(jié)構(gòu)同方法2。
方法4(DSAE):DSAE的結(jié)構(gòu)為2000?1000?500?200?50?7,每個稀疏自編碼器的學(xué)習(xí)率為0.14,動量因子為0.1,迭代次數(shù)為200。
方法5(MEWT+DSAE):DSAE結(jié)構(gòu)同方法4。
方法6(原始EWT+IFractalNet):IFractalNet結(jié)構(gòu)同方法1。
方法7(MEWT+CNN):CNN結(jié)構(gòu)同文獻(xiàn)[6]。
方法8(MEWT+原始FractalNet):FractalNet同方法1,但其損失函數(shù)和激活函數(shù)同文獻(xiàn)[13]。
表2列出了測試階段的平均識別準(zhǔn)確率與標(biāo)準(zhǔn)差,圖13列出了在每次試驗中采用不同方法測試的詳細(xì)識別結(jié)果。
由表2可知,與其他方法相比,本文方法具有更高的故障識別準(zhǔn)確率和更小的標(biāo)準(zhǔn)差,平均識別準(zhǔn)確率達(dá)99.78%,標(biāo)準(zhǔn)差僅為0.08。方法2與方法3、方法4與方法5、方法1與方法6的對比結(jié)果表明MEWT的分解降噪效果優(yōu)于原始EWT;方法1與方法7、方法8的對比結(jié)果表明,IFractalNet的特征提取和識別結(jié)果優(yōu)于CNN和原始FractalNet。本文提出方法利用MEWT自適應(yīng)劃分滾動軸承振動信號頻譜的邊界,從而確定分解模態(tài)然后利用加權(quán)評價指標(biāo)進(jìn)行信號重建,并改進(jìn)FractalNet的損失函數(shù)和激活函數(shù),相較于其他方法具有更高的故障識別能力。圖14給出了本文提出方法的第1次測試結(jié)果的混淆矩陣,可以看出,復(fù)合故障工況f和g的識別準(zhǔn)確率較低。
實際情況中,正常工況的樣本比例要大于故障工況樣本比例,進(jìn)一步研究本文提出方法在不平衡樣本比例下的故障識別性能。共設(shè)置5組不同比例的不平衡樣本,正常工況和各故障工況的訓(xùn)練樣本數(shù)量比例分別設(shè)置為8000∶8000,8000∶7000,8000∶6000,8000∶5000和8000∶4000,共進(jìn)行10次實驗,取平均結(jié)果,同時與方法6、方法7和方法8進(jìn)行對比分析,識別準(zhǔn)確率如表3所示。
由表3可知,在不平衡訓(xùn)練樣本情況下,本文方法(方法1)的故障識別效果仍然優(yōu)于其他幾種方法;在第5組中,即當(dāng)不平衡比例達(dá)8000∶4000時,本文方法的平均識別準(zhǔn)確率達(dá)到95.14%。綜上,4種方法的識別準(zhǔn)確率均隨不平衡比例的增大整體呈現(xiàn)下降趨勢,但本文方法在面對不平衡訓(xùn)練樣本時表現(xiàn)更加優(yōu)異。
以組5為例,計算了本文方法和方法8基于不平衡樣本的準(zhǔn)確率(P)、召回率(R)和F1值,如下式所示
式中 TP為正樣本中被模型判斷為正的個數(shù),F(xiàn)P為負(fù)樣本中被模型判斷為正的個數(shù),F(xiàn)N為正樣本中被模型判斷為負(fù)的個數(shù),F(xiàn)1值范圍在[0,1]之間,0代表最差,1代表最好。表4列出了組5中兩種方法的計算結(jié)果,可知,組5中本文方法的P,R和F1指標(biāo)值均較高,類似的結(jié)果在其他組中也有明顯體現(xiàn)。對比結(jié)果進(jìn)一步驗證了本文方法在面對不平衡樣本時故障識別的優(yōu)越性和有效性。
作為深層網(wǎng)絡(luò)傳遞信息的“門”,激活函數(shù)對于深層網(wǎng)絡(luò)至關(guān)重要。本文研究ReLU,Swish和GELU函數(shù)對IFractalNet識別性能的影響,圖15為不同激活函數(shù)下IFractalNet的訓(xùn)練精度隨迭代次數(shù)的變化曲線,可見8000次迭代后,GELU函數(shù)的識別準(zhǔn)確率比ReLU的訓(xùn)練精度提高了7.9%,比Swish提高了5.8%,表明GELU函數(shù)具有較好的收斂精度。
4 結(jié) 論
本文提出一種基于MEWT和IFractalNet的滾動軸承故障識別方法,主要結(jié)論如下:
(1)提出一種形態(tài)學(xué)EWT分解方法,利用形態(tài)濾波提取振動信號頻譜中具有顯著特征的局部極大值作為頻譜劃分邊界,自適應(yīng)地分割滾動軸承振動信號的頻譜,并利用輔助白噪聲優(yōu)化EWT,實現(xiàn)更精確的EWT自適應(yīng)模態(tài)分解,仿真信號和實際信號實驗表明MEWT的分解降噪效果優(yōu)于原始EWT,對噪聲的魯棒性更強(qiáng)。
(2)針對CNN層數(shù)越深,梯度消失現(xiàn)象越明顯,進(jìn)而造成識別率越低的缺陷,將FractalNet引入滾動軸承故障識別領(lǐng)域,并針對不平衡數(shù)據(jù)集改進(jìn)FractalNet的損失函數(shù),對樣本較多的故障類別賦以較小權(quán)重,對樣本較少的故障類型賦以較大的權(quán)值;將Fisher懲罰項引入FractalNet的損失函數(shù),使網(wǎng)絡(luò)學(xué)習(xí)到的隱層特征對振動信號的結(jié)構(gòu)變化具有可區(qū)分性;為緩解ReLU函數(shù)造成的“神經(jīng)元死亡”現(xiàn)象,將GELU作為IFractalNet的激活函數(shù),實驗表明IFractalNet的自動特征提取與故障識別能力均優(yōu)于CNN和原始FractalNet,在不平衡訓(xùn)練樣本集上也取得較好的效果。
參考文獻(xiàn):
[1] 王奉濤, 鄧? 剛, 王洪濤. 基于EMD和SSAE的滾動軸承故障診斷方法[J].振動工程學(xué)報, 2019, 32(2): 368-376.
Wang Fengtao, Deng Gang, Wang Hongtao, et al. A? ?rolling bearing fault diagnosis method based on EMD? ?and SSAE[J]. Journal of Vibration Engineering,? ?2019, 32(2): 368-376.
[2] Pan H Y, Yang Y, Li X, et al. Symplectic geometry mode decomposition and its application to rotating machinery compound fault diagnosis[J]. Mechanical Systems and Signal Processing, 2019, 114: 189-211.
[3] Gilles J. Empirical wavelet transform[J]. IEEE Transactions on Signal Processing, 2013, 61(16): 3999-4010.
[4] Jurgen S. Deep learning in neural networks: An overview[J]. Neural Networks, 2015, 61(1): 85-117.
[5] Chen Z Q, Li C, Sanchez R V. Gearbox fault identification and classification with convolutional neural networks[J]. Shock and Vibration, 2015, (2):1-10.
[6] Wang L H, Zhao X P, Wu J X, et al. Motor fault diagnosis based on short-time Fourier transform and convolutional neural network[J]. Chinese Journal of Mechanical Engineering, 2017, 30(6): 1357-1368.
[7] He K, Zhang X, Ren S, et al. Deep residual learning for image recognition[C]. 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 2016: 770-787.
[8] Zhao M, Kang M, Tang B, et al. Deep residual networks with dynamically weighted wavelet coefficients for fault diagnosis of planetary gearboxes[J]. IEEE Transactions on Industrial Electronics, 2018, 65(5): 4290-4300.
[9] Yeh Jia-Rong, Shieh Jiann-Shing, Huang Norden E. Complementary ensemble empirical mode decomposition: A novel noise enhanced data analysis method[J]. Advances in Adaptive Data Analysis, 2011, 2(2):135-156.
[10] 姜海旭, 張 科, 王靖宇, 等. 基于形態(tài)變分模態(tài)分解和JRD的航天器異常狀態(tài)識別[J]. 西北工業(yè)大學(xué)學(xué)報, 2018, 36(1): 20-27.
Jiang Haixu, Zhang Ke, Wang Jingyu, et al. Spacecraft anomaly recognition based on morphological variational mode decomposition and JRD[J]. Journal of Northwestern Polytechnical University, 2018, 36(1): 20-27.
[11] 劉 姝. 數(shù)學(xué)形態(tài)學(xué)在信號處理方面的應(yīng)用研究[D].大連: 大連理工大學(xué),2006.
Liu Shu.Applying the mathematical morphology in signal processing[D]. Dalian: Dalian University of Technology, 2006.
[12] 何洋洋, 呂躍剛, 劉俊承. 基于VMD與粒子濾波的滾動軸承故障診斷[J].可再生能源, 2019, 37(1): 126-131.
He Yangyang, Lü Yuegang, Liu Juncheng. Fault diagnosis of rolling bearing based on VMD and particle filter[J]. Renewable Energy Resources,2019, 37(1): 126-131.
[13] Larsson G, Maire M, Shakhnarovich G. FractalNet: Ultra-deep neural networks without residuals[J]. 2016, arXiv:1605.07648.
[14] 馮新?lián)P, 張巧榮, 李慶勇. 基于改進(jìn)型深度網(wǎng)絡(luò)數(shù)據(jù)融合的滾動軸承故障識別[J]. 重慶大學(xué)學(xué)報, 2019, 42(2): 52-62.
Feng Xinyang, Zhang Qiaorong, Li Qingyong. Fault recognition of rolling bearing based on improved deep networks with data fusion in unbalanced data sets[J]. Journal of Chongqing University, 2019, 42(2): 52-62.
[15] 陳志剛, 杜小磊, 王衍學(xué), 等.改進(jìn)集成深層自編碼器在軸承故障診斷中的應(yīng)用[J].控制與決策, 2021, 36(1): 135-142.
Chen Zhigang, Du Xiaolei, Wang Yanxue, et al. Application of improved ensemble deep auto-encoder in bearing fault diagnosis[J]. Control and Decision, 2021, 36(1): 135-142.
[16] Qu J X, Zhang Z S, Gong T. A novel intelligent method for mechanical fault diagnosis based on dual-tree complex wavelet packet transform and multiple classifier fusion[J]. Neurocomputing, 2016, 171: 837-853.
作者簡介: 杜小磊(1993-),男,碩士研究生。電話: 13260417530;E-mail:13260417530@163.com
通訊作者: 陳志剛(1979-),男,副教授。電話:13521909643;E-mail:gangzi22@163.com