陳 帥,黃晉英
(中北大學(xué) 機(jī)械工程學(xué)院,太原 030051)
滾動(dòng)軸承作為機(jī)械設(shè)備的重要部件之一,因其復(fù)雜的結(jié)構(gòu)和惡劣的工作環(huán)境,在長(zhǎng)時(shí)間運(yùn)行后可能會(huì)出現(xiàn)各種各樣的故障,結(jié)果可能導(dǎo)致企業(yè)生產(chǎn)效率降低嚴(yán)重時(shí)甚至出現(xiàn)生命安全問(wèn)題[1]。所以,對(duì)其故障診斷方法進(jìn)行系統(tǒng)全面的探究具有十分重要的實(shí)際意義。如果發(fā)生了故障,相關(guān)的振動(dòng)信號(hào)往往會(huì)具有非平穩(wěn)以及非線性等方一系列特征?,F(xiàn)有的非線性、非平穩(wěn)性以熵值為特征參數(shù)的信號(hào)分析方法主要有:近似熵、排列熵、模糊熵、樣本熵等,但上述方法均為基于時(shí)間序列的單一尺度分析,分類效果較差[2]。之后又提出了多尺度熵算法,其優(yōu)點(diǎn)是通過(guò)多尺度粗?;^(guò)程使得熵?cái)U(kuò)展到多個(gè)時(shí)間尺度上,從多尺度表征故障信號(hào)的特征信息,進(jìn)而解決了單一尺度熵的局限性,減小了偏差[3]。但是當(dāng)數(shù)據(jù)序列較長(zhǎng)時(shí),會(huì)出現(xiàn)運(yùn)算速度慢、故障信息丟失嚴(yán)重等問(wèn)題,基于此,學(xué)者們提出一種精細(xì)復(fù)合多尺度散布熵(RCMDE)算法[4],該方法能夠明顯改善在粗?;^(guò)程中信息丟失問(wèn)題。近年來(lái),神經(jīng)網(wǎng)絡(luò)已經(jīng)獲得了廣泛的應(yīng)用,它在非線性關(guān)系數(shù)據(jù)的處理領(lǐng)域中尤其顯著的優(yōu)勢(shì)。概率神經(jīng)網(wǎng)絡(luò)(PNN)的主要特點(diǎn)是有良好的并行處理、自學(xué)習(xí)以及自組織等方面的能力,可以在不同的精度水平之下向各種連續(xù)非線性函數(shù)進(jìn)行逼近[5]。在此基礎(chǔ)上,本文以概率神經(jīng)網(wǎng)絡(luò)(也就是PNN)以及精細(xì)復(fù)合多尺度散布熵(也就是RCMDE)為基礎(chǔ)首次提出了一種新的滾動(dòng)軸承故障診斷分析方法,在研究的過(guò)程中,還和MDE-PNN法展開(kāi)全面的分析與對(duì)比。所得到的結(jié)果表明,所提方法既能夠?qū)崿F(xiàn)軸承故障的有效識(shí)別分類,且診斷效果優(yōu)于MDE-PNN方法,平均分類精度達(dá)到了97.65%。
散布熵是度量時(shí)間序列復(fù)雜性的非線性動(dòng)力分析方法,對(duì)于長(zhǎng)度為N的時(shí)間序列x={x1,x2,…,xN},其計(jì)算步驟如下[6]:
1)利用正態(tài)分布函數(shù)將時(shí)間序列x映射到y(tǒng)={yj,j=1,2,…,N},yi∈(0,1),如式(1)所示。
式(1)中,σ和μ分別表示標(biāo)準(zhǔn)差和均值。
2)線性變換將y映射到[1,2,…,c],如式(2)所示。
式(2)中,c是類別個(gè)數(shù),int()是取整函數(shù)。
3)計(jì)算嵌入向量zm,ci:
式(3)中,m是嵌入維數(shù),d是時(shí)間延遲。
7)計(jì)算標(biāo)準(zhǔn)散布熵,公式如式(7)所示。
RCMDE計(jì)算過(guò)程如下[7]:
1)多尺度粗粒化。給定原始數(shù)據(jù)u,其第k個(gè)粗粒化序列的計(jì)算如式(8)所示。
式(8)中,τ是尺度因子。
2)對(duì)每個(gè)尺度τ,RCMDE值定義如式(9)所示。
它為徑向基網(wǎng)絡(luò)當(dāng)中的重要分支,主要的優(yōu)點(diǎn)就是分類精度比較高,可以有效地利用簡(jiǎn)單的線性模型訓(xùn)練非線性模型,有效避免了BP神經(jīng)網(wǎng)絡(luò)的局部最小值問(wèn)題[8]。
PNN模型分為4層結(jié)構(gòu),第1層是輸入層,第2層是具有N個(gè)神經(jīng)元的隱含層(N是訓(xùn)練數(shù)據(jù)的數(shù)量),第i類樣本的第j個(gè)神經(jīng)元的輸入與輸出關(guān)系由式(10)確定[9]。
式(10)中是平滑因子,d為數(shù)據(jù)維度。
第3層是求和層,其把上一層中為同類的隱含層神經(jīng)元的輸出按式(11)平均:
式(11)中,fi為第i類的概率。
第4層是神經(jīng)元的輸出層,如式(12)所示。
式(12)中,z為輸出層的輸出。
在本研究中將會(huì)通過(guò)西儲(chǔ)大學(xué)數(shù)據(jù)集來(lái)進(jìn)行全面的分析與驗(yàn)證。所采用的實(shí)驗(yàn)臺(tái)詳見(jiàn)本文的圖1,其中有電子控制器、功率測(cè)試分析設(shè)備、扭矩傳感器以即電動(dòng)機(jī)等多個(gè)重要的組成部分。軸承選用SKF6205-2RS,為了模擬軸承常見(jiàn)的4種狀態(tài),采用電火花加工方式對(duì)軸承進(jìn)行不同程度的單點(diǎn)損傷。選擇采樣頻率12kHz、電機(jī)轉(zhuǎn)速為1797r/min的驅(qū)動(dòng)端軸承數(shù)據(jù)進(jìn)行分析,每種狀態(tài)選擇80個(gè)樣本,共計(jì)320個(gè)樣本,每個(gè)樣本長(zhǎng)度為4096,軸承狀態(tài)說(shuō)明見(jiàn)表1。圖2為軸承時(shí)域波形圖。
圖1 實(shí)驗(yàn)臺(tái)
圖2 軸承4種狀態(tài)時(shí)域波形圖
表1 軸承的4種狀態(tài)
計(jì)算320個(gè)樣本的RCMDE值,設(shè)置算法相關(guān)參數(shù):最大尺度因子τ=20,嵌入維數(shù)m=2,延遲時(shí)間d=1。每種狀態(tài)不同樣本計(jì)算的熵值均值如圖3所示,可知,正常軸承振動(dòng)信號(hào)RCMDE值遠(yuǎn)大于3種故障狀態(tài)RCMDE值,并且變化趨勢(shì)明顯不同。對(duì)于排在前面的三個(gè)尺度,正常情況下,RCMDE值會(huì)隨著尺度逐漸提高而相依地提升,如果是出于故障的狀態(tài)中,那么它將會(huì)是始終下降的,尤其是由第三個(gè)尺度至第十個(gè)的過(guò)程中,將出現(xiàn)一定程度的波動(dòng)現(xiàn)象,在此之后將比較的穩(wěn)定。從整體上來(lái)看,RCMDE值與尺度因子成反比,在尺度因子達(dá)到11之后逐漸趨于平穩(wěn)。
圖3 振動(dòng)信號(hào)RCMDE均值
此外,還與MDE法展開(kāi)全面的對(duì)比分析,以闡述RCMDE法所具有的主要優(yōu)勢(shì)。計(jì)算每種狀態(tài)信號(hào)不同樣本MDE值,所得結(jié)果如圖4所示。對(duì)比圖3、圖4,可以看出:隨著尺度因子τ的增大,RCMDE值變化較MDE值變化相對(duì)平緩,從圖4可知,外圈和內(nèi)圈故障從第5個(gè)尺度因子之后基本重合,不利于下一步分類器的識(shí)別分類,所以,選擇RCMDE值作為特征參數(shù)進(jìn)行軸承4種狀態(tài)的特征提取較為合理。
圖4 前20個(gè)尺度MDE均值
實(shí)驗(yàn)隨機(jī)選取每種狀態(tài)的50個(gè)樣本作為訓(xùn)練樣本,30個(gè)樣本作為測(cè)試樣本,為了減小實(shí)驗(yàn)偶然性,對(duì)RCMDEPNN和MDE-PNN兩種模型分別進(jìn)行10次實(shí)驗(yàn),結(jié)果如圖5所示。根據(jù)圖5可知,基于RCMDE-PNN的故障診斷方法分類效果明顯優(yōu)于比MDE-PNN方法,分類準(zhǔn)確率達(dá)到97.65%,而MDE-PNN的分類準(zhǔn)確率為92.85%,分類準(zhǔn)確率提高了4.8%。
圖5 實(shí)驗(yàn)結(jié)果
1)引入一種檢驗(yàn)振動(dòng)信號(hào)復(fù)雜性的算法-精細(xì)復(fù)合多尺度散布熵(RCMDE),與MDE方法相比,RCMDE解決了單一尺度熵的局限性,實(shí)現(xiàn)了多尺度表征故障信號(hào)的特征信息,有利于下一步的故障識(shí)別分類;
2)RCMDE-PNN模型的軸承故障診斷準(zhǔn)確率為97.65%,而MDE-PNN模型診斷準(zhǔn)確率僅為92.85%,本文所提RCMDE-PNN方法,提升其分類準(zhǔn)確率達(dá)4.8%,并且對(duì)它的優(yōu)勢(shì)與有效性進(jìn)行了全面的驗(yàn)證。