基于多重概率機器算法的齒輪箱耦合故障診斷

鐘建華，鐘 敏，施 雯

（1.福州大學(xué) 機械工程及自動化學(xué)院，福州 350108；2.福安職業(yè)技術(shù)學(xué)校，寧德 355000）

0 引言

旋轉(zhuǎn)設(shè)備被廣泛用于風(fēng)力發(fā)電、直升機、工程機械等大型復(fù)雜機械裝備中，工作環(huán)境惡劣經(jīng)常導(dǎo)致其齒輪和軸承部件出現(xiàn)嚴(yán)重磨損或疲勞裂紋等故障同時發(fā)生情況[1]。由于不同的故障激勵源相互作用，導(dǎo)致復(fù)雜的故障耦合形態(tài)，增加了故障診斷的難度。例如：齒輪斷齒的同時還存在軸承疲勞剝落故障，技術(shù)人員在現(xiàn)場排查中更換了斷齒的齒輪，由于個人經(jīng)驗不足未發(fā)現(xiàn)疲勞剝落的軸承故障，將造成機組在復(fù)機運行不久后的二次停機待修，進(jìn)一步增加運行與維護(hù)成本。

目前少部分研究人員在旋轉(zhuǎn)機械領(lǐng)域?qū)︸詈瞎收险归_研究并取得初步成果。Patel等人[2,3]研究轉(zhuǎn)子不平衡、裂紋與動-靜碰摩耦合故障的動力學(xué)特性。焦衛(wèi)東等人[4]利用時-頻分析方法研究齒輪表面破壞、機座松動、齒輪表面破壞+機座松動耦合故障情況下故障征提取，然后進(jìn)一步利用灰色關(guān)聯(lián)分析對耦合故障特性進(jìn)行定量評價，從而實現(xiàn)復(fù)雜耦合故障模式的識別。劉洋等人[5]和許琦等人[6]采用頻譜諧波分析方法研究軸承不對中、松動、裂紋已經(jīng)碰摩擦等耦合故障診斷問題。韓東穎等人[7]提出基于Empirical Mode Decomposition（EMD）經(jīng)驗?zāi)B(tài)分解方法對旋轉(zhuǎn)機械耦合故障信號進(jìn)行分析。

本文提出多重概率機器算法聯(lián)合診斷方法，識別齒輪箱耦合故障中單故障的數(shù)量與類型。該方法同時提取振動信號與聲音信號在集合經(jīng)驗?zāi)B(tài)分解（EEMD）下本征模態(tài)函數(shù)（IMF）的能量特征，并結(jié)合信號的時域統(tǒng)計特征，在此故障特征基礎(chǔ)上訓(xùn)練多重概率分類器（SBELM），為平衡各概率分類器的性能差異，利用粒子群算法（PSO）優(yōu)化分類器參數(shù)與閾值，最終構(gòu)建旋轉(zhuǎn)設(shè)備耦合故障診斷的新方法。

1 相關(guān)理論

1.1 基于EEMD的能量特征提取

EEMD通過在原始信號中多次添加頻率均勻的高斯白噪聲，經(jīng)過多次對IMF分量求平均以消除加入的噪聲影響，可有效降低由EMD極值分布不均造成的模態(tài)混疊問題[8,9]，具體步驟如下：

1）向信號添加一組均勻的高斯白噪聲，生產(chǎn)新的信號x(t)；

2）對新信號進(jìn)行EMD，得到n個本征模態(tài)函數(shù)（intinsic mode function，IMF）和一個余項：

其中xj對應(yīng)的第j個本征模態(tài)函數(shù)；rn為分解到最后剩余的余項；

3）迭代并重復(fù)步驟1）和2），迭代次數(shù)通常設(shè)置高于100次，最后得到的結(jié)果：

式中xi(t)為第i次迭代產(chǎn)生的新信號；xij對應(yīng)的第i次迭代時EMD產(chǎn)生第j個IMF；rin為第i次迭代EMD分解的余項；

4）對步驟3所得到的各個IMF和余項去均值得到最終的IMF和余項：

其中N為迭代次數(shù)；

5）計算各IMF與信號之間的相關(guān)系數(shù)，選取適合的IMF：

式中，Ii(t)為第i個IMF，x(t)為原始信號；分別是x(t)和Ii(t)的平均值；n為IMF的個數(shù)；值越大表示x(t)和Ii(t)之間相關(guān)性越高，從而選取出適合的IMF個數(shù)，剔除剩余IMF和余項；

6）計算剩余各IMF的能量特征：

式中tΔ為時間間隔，m為信號點數(shù)，Ei為第i個IMF的能量特征值；

經(jīng)過EEMD處理后，原始數(shù)據(jù)在不同的時間尺度上局部特征信號被抽取形成IMF分量，其能量特征進(jìn)一步反應(yīng)了不同特征尺度隨機序列的能量分布情況。

1.2 基于稀疏性貝葉斯極限學(xué)習(xí)機（SBELM）的模式識別

稀疏貝葉斯極限學(xué)習(xí)機（Sparse Bayesian Extreme Learning Machine）通過輸出層權(quán)值參數(shù)W上定義0均值高斯分布來獲得稀疏的模型解，使ELM具有良好的稀疏性能[10]。隱含層輸出H是SBELM的輸入：

其中θ是輸入與隱含層之間的權(quán)值；b為隱含層節(jié)點閾值。在處理分類問題時，通過多項概率分布方程給出輸出概率分布：

本文在處理Multi-classification（多類別標(biāo)簽/多重故障數(shù)據(jù)）將采用一對一（one-versus-one）策略。一對一策略消除了多對一策略產(chǎn)出的灰色不確定區(qū)間影響，提高分類精度如圖1所示。具體方法描述請參考[11]。

圖1 一對一與一對多策略比較

2 診斷流程

第一步是信號預(yù)處理。首先通過EEMD對采集到的振動和聲音信號進(jìn)行分解，并利用相關(guān)性分析（CC）如表1所示，其中I8和I9明顯小于其他，本文選取前7個IMF，并提取IMF的能量特征；同時，為豐富故障特征的數(shù)量提高分類精度，提取原信號的6個常用的時域統(tǒng)計特征值，如表2所示。由此，IMF能量特征與時域統(tǒng)計特征合成13個提取的故障特征值，并將提取的故障特征拆分成訓(xùn)練數(shù)據(jù)集、驗證數(shù)據(jù)集和測試數(shù)據(jù)集。

表1 振動/聲音號相關(guān)性分析

表2 信號時域統(tǒng)計指標(biāo)

圖2 多重概率聯(lián)合診斷機制

第二步是訓(xùn)練并搭建多重概率分類器聯(lián)合診斷機制（MPC）。利用訓(xùn)練數(shù)據(jù)集訓(xùn)練振動信號分類器和聲音信號分類器；

第三步是參數(shù)優(yōu)化利用PSO優(yōu)化分類器的權(quán)重值(w)，平衡分類器之間的性能差異，利用式(9)構(gòu)建多重概率分類器診斷機制，從而輸出公平公正的診斷結(jié)果；并優(yōu)化閾值，選取正確的單故障數(shù)量與類型。

式中Pi表示第i個單故障發(fā)生的概率；wj-opt為第j個概率分類器權(quán)重值；ρj-i表示第j個概率分類器下的第i個單故障發(fā)生的概率值。

第四步性能測試。利用測試數(shù)據(jù)集，測試多重概率聯(lián)合診斷機制模型。

3 實驗平臺與數(shù)據(jù)采集

實驗設(shè)備為美國SpectraQuest設(shè)計的動力傳動系統(tǒng)實驗平臺（DDS）。實驗平臺的布置如圖3所示，該實驗臺可模擬直齒和斜齒的齒面磨損、輪齒裂紋和缺齒，也可模擬滾動軸承內(nèi)圈、外圈、滾動體等故障。可引入單一故障，或同時引入多個單故障耦合，研究其相互間的耦合效應(yīng)。并可通過3馬力交流變頻驅(qū)動電機編程自定義速度來加載，來模擬不同負(fù)載情況下對設(shè)備運行狀態(tài)的影響。本文主要分析四種常見的齒輪箱故障，以及兩種耦合性故障如圖5和表3所描述。

單/耦合故障的振動和聲音信號的采樣頻率8092Hz，負(fù)載分三擋分別是1、2和3馬力加載。單/耦合故障各采集100組，拆分成訓(xùn)練、驗證和測試數(shù)據(jù)集如表4所示。

圖3 動力傳動系統(tǒng)實驗平臺（DDS）

圖4 齒輪與軸承故障類型

圖5 故障振動和聲音信號圖

表3 故障類型描述

表4 信號數(shù)據(jù)拆分

4 實驗結(jié)果分析

4.1 不同特征提取方法的比較

為了驗證本文中故障特征提取方法的有效性，與常用的特征提取算法比較，結(jié)果如表5所示。

表5 信號特征提取方法比較

上述表中，PNN、RVM和SBELM算法中參數(shù)采用默認(rèn)值；MPC是由兩個SBELM組成，它們分別是振動信號和聲音信號訓(xùn)練而成，權(quán)重值與閾值分別采用默認(rèn)值（0.5和1）。結(jié)果表明，相比較于其他特征提取方法，本文基于EEMD+TDSF的故障特征提取方法更佳有效的提高分類器的識別準(zhǔn)確性。

4.2 參數(shù)優(yōu)化結(jié)果討論

為了驗證粒子群算法PSO對診斷模型中的參數(shù)優(yōu)化的有效性，將優(yōu)化結(jié)果與默認(rèn)參數(shù)結(jié)果箱比較，結(jié)果如表6所示。

表6 參與優(yōu)化結(jié)果比較（單/耦混合測試）

對比結(jié)果表明，PSO對參數(shù)的優(yōu)化提高1.35%診斷精度，因此本文的多重概率聯(lián)合診斷機制是在PSO參數(shù)優(yōu)化結(jié)果上建立的。

4.3 診斷方法比較結(jié)果討論

為了進(jìn)一步驗證本文提出的多重概率分類器聯(lián)合診斷模型（MPC）的有效性，MPC分別與PNN、RVM和SBELM算法在單故障、耦合故障和混合故障（包含單/耦合故障）三類測試數(shù)據(jù)的分類結(jié)果進(jìn)行比較，比較結(jié)果如表7所示。

表7 MPC與PNN、RVM和SBELM比較

圖6 SBELM1（振動信號）對35個C5數(shù)據(jù)估計結(jié)果

圖7 SBELM2（聲音信號）對35個C5數(shù)據(jù)估計結(jié)果

由上述表格比較結(jié)果表明，本文提出的MPC診斷模型比較PNN，RVM和SBELM算法在耦合故障分類準(zhǔn)確性上分別有2.79%，2%和1.47%提高，并且在單故障和混合故障也達(dá)到96.32%和92.56%。由此可見，MPC相比其他分類器更能有效的識別出單/耦合故障。MPC中的SBELM1和SBELM2概率分類器對混合故障中的35個耦合故障C5（齒面破損+軸承內(nèi)圈破損）分類結(jié)果如圖6和圖7所示。圖6和圖7對可能發(fā)生的故障進(jìn)行大小排序，大于閾值0.628（紅色線）被認(rèn)定為故障存在，MPC正確判斷出耦合故障C5是由齒面破損和軸承內(nèi)圈磨損組成。工程師也可根據(jù)排序結(jié)果進(jìn)行故障先后排查，降低設(shè)備在開機后停機待修事件發(fā)生。

5 結(jié)論

為了解決工程中二次停機待修問題，本文基于概率機器算法SBELM，提出一種新的多重概率分類器聯(lián)合診斷模型（MPC）。首先，通過EEMD對振動與聲音信號預(yù)處理，提取IMF能量和原信號的時域統(tǒng)計特征；再利用提取振動和聲音信號特征分別訓(xùn)練SBELM，并優(yōu)化分類器參數(shù)與閾值，從而構(gòu)建多重概率分類器聯(lián)合診斷機制。實驗結(jié)果表明：本文提出的多重概率分類器聯(lián)合診斷框架，不僅能正確識別出單故障類型，同時還能有效識別耦合故障中的單故障數(shù)量以及類型。診斷結(jié)果將為工程師提供依據(jù)，逐一排查可能發(fā)生的故障，從而避免復(fù)機后二次停機問題。