RBF神經(jīng)網(wǎng)絡(luò)在旋轉(zhuǎn)機(jī)械故障診斷中的應(yīng)用

陳 松

（安徽建筑大學(xué) 機(jī)械與電氣工程學(xué)院，合肥 230601）

旋轉(zhuǎn)機(jī)械故障的診斷可以通過測量振動信號、振聲、溫度、油液或光譜來實現(xiàn).對于振動信號有時域特征提取、時序特征提取、頻譜特征提取、時頻特征提取、高階譜特征提取和小波頻帶特征提取等方法.對于各種特征量的識別，應(yīng)用較多的算法有：模糊集合、專家系統(tǒng)、神經(jīng)網(wǎng)絡(luò)等［1-3］.故障特征與故障模式并不是簡單的一一對應(yīng)關(guān)系，其構(gòu)成的故障特征空間比較復(fù)雜，常常不是線性可分的，而神經(jīng)網(wǎng)絡(luò)能夠映射任意復(fù)雜的非線性關(guān)系，具有自學(xué)習(xí)、自組織、自適應(yīng)等特性，并且有極強(qiáng)的容錯和聯(lián)想能力、較快的計算速度，所以神經(jīng)網(wǎng)絡(luò)被廣泛用于機(jī)械故障診斷識別中.

目前廣泛應(yīng)用的是BP神經(jīng)網(wǎng)絡(luò)，它能反映診斷過程的本質(zhì)，靜態(tài)模式識別能力強(qiáng)，可實現(xiàn)任意復(fù)雜的判決表面，具有自學(xué)習(xí)及自適應(yīng)能力［4］.但BP網(wǎng)絡(luò)采用的是梯度下降的搜索算法，這就不可避免地出現(xiàn)了網(wǎng)絡(luò)收斂速度慢、容易陷入局部極小、學(xué)習(xí)結(jié)果受初始權(quán)值分布影響較大、結(jié)果不穩(wěn)定等問題.Elman神經(jīng)網(wǎng)絡(luò)作為一種重要的動態(tài)反饋性學(xué)習(xí)模型，通過增加結(jié)構(gòu)單元對各個樣本的相互關(guān)聯(lián)進(jìn)行記憶，實現(xiàn)動態(tài)建模，該模型具有極強(qiáng)的復(fù)雜模式動態(tài)映射能力，更是一種極具潛力的復(fù)雜故障模式辨識工具.但由于Elman神經(jīng)網(wǎng)路的學(xué)習(xí)過程與前饋神經(jīng)網(wǎng)絡(luò)類似，難免會出現(xiàn)收斂速度慢和易收斂到局部極小的缺陷，導(dǎo)致故障模式辨識結(jié)果不穩(wěn)定［5］.

徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)（Radical Basis Function神經(jīng)網(wǎng)絡(luò)）是一種高效的前饋式神經(jīng)網(wǎng)絡(luò)，它具有其他前向網(wǎng)絡(luò)所不具有的最佳逼近性能和全局最優(yōu)特性，可以避免陷入局部極小的可能，并且結(jié)構(gòu)簡單，訓(xùn)練速度快.它也是一種可以廣泛應(yīng)用于模式識別、非線性函數(shù)逼近等領(lǐng)域的神經(jīng)網(wǎng)絡(luò)模型.所以本文將徑向基神經(jīng)網(wǎng)絡(luò)用于旋轉(zhuǎn)機(jī)械實現(xiàn)故障診斷，使用Matlab軟件編程并得出識別結(jié)果.

1 徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)算法

RBF網(wǎng)絡(luò)結(jié)構(gòu)如圖1所示.輸入層到隱層為權(quán)值為1的固定連接，隱層到輸出層為權(quán)值為W 的線性鏈接.隱層神經(jīng)元基函數(shù)常采用高斯型徑向基函數(shù).

圖1 RBF網(wǎng)絡(luò)結(jié)構(gòu)

隱層函數(shù)：

輸出層函數(shù)：

式中，cr為第r個隱層節(jié)點的數(shù)據(jù)中心，σr是第r個隱層節(jié)點的數(shù)據(jù)方差，X是神經(jīng)網(wǎng)絡(luò)的輸入向量，Wmr是隱層到輸出層的連接權(quán).

徑向基神經(jīng)網(wǎng)絡(luò)算法步驟如下：

1）從輸入向量中選一組初始中心值cr，初始化連接權(quán)值Wmr，計算方差值σ＝dmax／R，其中dmax是最大的距離；R是cr的數(shù)量；

2）用聚類分析中的K－NN算法求取cr和σr；

3）計算輸出層權(quán)值：

計算網(wǎng)絡(luò)輸出

式中，em（n）＝y(tǒng)m（n）-dm（n），dm（n）為樣本目標(biāo)輸出；μW是學(xué)習(xí)參數(shù)的學(xué)習(xí)步長.

4）計算網(wǎng)絡(luò)誤差并判斷是否收斂.

5）若收斂，訓(xùn)練結(jié)束.

6）否則，更新網(wǎng)絡(luò)學(xué)習(xí)參數(shù)cr，σr和Wmr，并轉(zhuǎn)到步驟（2）［6-7］.

2 機(jī)械故障的輸入和輸出樣本設(shè)置

本設(shè)計中可以識別的故障包括：轉(zhuǎn)子不平衡，角度不對中，平行不對中，油膜渦動，油膜振蕩，喘振，軸彎曲，轉(zhuǎn)子松動.有的文獻(xiàn)中，對上述故障進(jìn)行診斷使用頻譜能量作為特征量，難以區(qū)分角度不對中與平行不對中、油膜渦動與油膜振蕩［8-9］；有的文獻(xiàn)中，對上述故障進(jìn)行診斷是識別軸心軌跡的特征，難以區(qū)分轉(zhuǎn)子不平衡與軸彎曲、喘振與轉(zhuǎn)子松動［10-11］.本文綜合使用頻譜能量，軸心軌跡特征作為機(jī)械故障的特征量，可以更有效地區(qū)分上述故障類型.

本設(shè)計使用MOBIUS公司的轉(zhuǎn)子機(jī)械故障仿真儀器，轉(zhuǎn)子轉(zhuǎn)速為1800r／min，設(shè)置各種機(jī)械故障，觀察振動信號頻譜以及軸心軌跡.其中，角度不對中對應(yīng)振動信號頻譜如圖2所示，角度不對中對應(yīng)軸心軌跡如圖3所示.

圖2 角度不對中對應(yīng)的振動信號頻譜

圖3 角度不對中對應(yīng)的軸心軌跡

轉(zhuǎn)子不平衡軸心軌跡為橢圓形，角度不對中軸心軌跡為香蕉形，平行不對中軸心軌跡為外八字形，油膜渦動軸心軌跡為內(nèi)八字形，油膜振蕩軸心軌跡為花形，軸彎曲軸心軌跡為橢圓形，喘振和轉(zhuǎn)子松動軸心軌跡非常紊亂，變化不定.在本設(shè)計中，提取的軸心軌跡的特征包括：孔洞數(shù)，質(zhì)心偏移程度，凹入個數(shù)，交叉點個數(shù)，各種故障對應(yīng)的這些特征量見表1.

表1 轉(zhuǎn)子故障對應(yīng)的軸心軌跡特征量數(shù)值

孔洞數(shù)、交叉點個數(shù)和凹入個數(shù)大于10的都記為10，以防止數(shù)值過大對其他參數(shù)造成影響.質(zhì)心偏移距離指圖形的質(zhì)心與軸心之間的距離，若質(zhì)心與軸心重合，質(zhì)心偏移程度記為0，若距離很大，質(zhì)心偏移程度記為1.

將（0.01～0.39）f、（0.40～0.49）f、（0.51～0.99）f、f、2f、（3～5）f等6個頻率段的頻譜能量作為特征頻率.根據(jù)頻譜信號能量計算公式（5），計算各個頻段頻譜能量［12］.其中W 表示能量，ω1～ω2為頻率范圍.

樣本特征在輸入神經(jīng)網(wǎng)絡(luò)之前必須進(jìn)行歸一化，原始的數(shù)據(jù)幅值大小不一，有時相差懸殊，如果直接使用，測量值大的波動就會壟斷了神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過程，使其不能反映小的測量值的變化.對頻譜能量進(jìn)行歸一化處理，得到表格2前6列的數(shù)值.對軸心軌跡特征量進(jìn)行歸一化處理，得到表格2后4列的數(shù)值.表格2為神經(jīng)網(wǎng)絡(luò)輸入樣本，表格3為神經(jīng)網(wǎng)絡(luò)輸出樣本.

表2 神經(jīng)網(wǎng)絡(luò)輸入樣本

表3 神經(jīng)網(wǎng)絡(luò)輸出樣本

3 試驗與分析

分別訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)、Elman神經(jīng)網(wǎng)絡(luò)和RBF神經(jīng)網(wǎng)絡(luò).BP網(wǎng)絡(luò)分別使用trainlm、traingdx、traingd 3種訓(xùn)練函數(shù)進(jìn)行試驗，trainlm為Levenberg-Marquardt算法，對于中等規(guī)模的BP神經(jīng)網(wǎng)絡(luò)有最快的收斂速度，訓(xùn)練中的計算量相對較少，但需要較大內(nèi)存量；traingdx（自適應(yīng)動量梯度下降法）為帶適應(yīng)學(xué)習(xí)率和動量因子的梯度遞減法；traingd為梯度下降訓(xùn)練函數(shù)，沿網(wǎng)絡(luò)性能參數(shù)的負(fù)梯度方向調(diào)整網(wǎng)絡(luò)的權(quán)值和閾值.訓(xùn)練次數(shù)為3000，訓(xùn)練目標(biāo)為0.001，學(xué)習(xí)速率為0.1，其中trainlm算法對應(yīng)的Matlab程序如下：

Elman神經(jīng)網(wǎng)絡(luò)訓(xùn)練次數(shù)為1000，訓(xùn)練目標(biāo)為0.001，對應(yīng)的主要程序如下：

RBF神經(jīng)網(wǎng)絡(luò)對應(yīng)的主要程序如下，其中GOAL為均方誤差，SPREAD為徑向基函數(shù)的分布密度，神經(jīng)元的最大數(shù)目為10，兩次顯示之間所添加的神經(jīng)元數(shù)目為1.

對每個網(wǎng)絡(luò)的識別能力進(jìn)行測試，輸入轉(zhuǎn)子不平衡、角度不對中、平行不對中對應(yīng)的特征量數(shù)據(jù)：

使用RBF神經(jīng)網(wǎng)絡(luò)識別結(jié)果為：

將網(wǎng)絡(luò)輸出與各模式閾值進(jìn)行比較，如果輸出結(jié)果大于預(yù)設(shè)閾值，則此故障發(fā)生，否則不發(fā)生該故障.本例中，閾值設(shè)為0.9，由網(wǎng)絡(luò)輸出結(jié)果，可以看出，3個故障模式對應(yīng)節(jié)點的輸出與待識別樣本的故障類型一致.從圖4可以看出訓(xùn)練誤差為7.21995×10-30，訓(xùn)練步數(shù)為7步.使用norm函數(shù)計算出識別誤差為0.0243.

圖4 訓(xùn)練誤差與步長

各種神經(jīng)網(wǎng)絡(luò)對應(yīng)的訓(xùn)練次數(shù)，訓(xùn)練誤差，識別誤差見表4.

表4 各種神經(jīng)網(wǎng)絡(luò)的性能

4 結(jié) 論

1）使用Matlab編程，得出識別結(jié)果，比較了幾種神經(jīng)網(wǎng)絡(luò)在機(jī)械故障識別中的性能.使用BP神經(jīng)網(wǎng)絡(luò)，trainlm學(xué)習(xí)算法可以識別，訓(xùn)練次數(shù)少，僅9次；traingdx學(xué)習(xí)算法訓(xùn)練次數(shù)多；traingd訓(xùn)練3000次，沒有識別.Elman神經(jīng)網(wǎng)絡(luò)計算次數(shù)多，時間長.RBF神經(jīng)網(wǎng)絡(luò)，訓(xùn)練次數(shù)少，速度快，準(zhǔn)確度高，識別誤差小.

2）使用頻譜能量和新的軸心軌跡特征作為機(jī)械故障的特征量，提供的信息更全面，可以有效地區(qū)別各種故障類型.

［1］鄧 堰.轉(zhuǎn)子故障智能診斷中的特征提取與選擇技術(shù)研究［D］.南京：南京航空航天大學(xué)，2008.

［2］張 鍵.機(jī)械故障診斷技術(shù)［M］.北京：機(jī)械工業(yè)出版社，2012.

［3］沈慶根，鄭水英.設(shè)備故障診斷［M］.北京：化學(xué)工業(yè)出版社，2009.

［4］鮑久圣.BP神經(jīng)網(wǎng)絡(luò)在機(jī)械故障診斷中的應(yīng)用［J］.現(xiàn)代制造工程，2005（1）：113－114.

［5］陳法法，湯寶平，黃慶卿.免疫遺傳優(yōu)化Elman神經(jīng)網(wǎng)絡(luò)的旋轉(zhuǎn)機(jī)械故障診斷［J］.重慶大學(xué)學(xué)報，2012（5）：7－13.

［6］汪慶華，王敬濤，鄧東花.基于徑向基神經(jīng)網(wǎng)絡(luò)的旋轉(zhuǎn)機(jī)械故障診斷［J］.現(xiàn)代電子技術(shù)，2010（18）：141－142.

［7］田景文，高美娟.人工神經(jīng)網(wǎng)絡(luò)算法研究及應(yīng)用［M］.北京：北京理工大學(xué)出版社，2006：41－43.

［8］肖 坤，原思聰，王 丹.SOM神經(jīng)網(wǎng)絡(luò)在旋轉(zhuǎn)機(jī)械故障診斷中的應(yīng)用［J］.機(jī)械設(shè)計與制造，2010（11）：44－46.

［9］欒美潔，許飛云，賈民平.旋轉(zhuǎn)機(jī)械故障診斷的神經(jīng)網(wǎng)絡(luò)方法研究［J］.噪聲與振動控制，2008（2）：85－87.

［10］肖圣光.基于虛擬儀器的軸心軌跡分析儀的研制［D］.重慶：重慶大學(xué)，2009.

［11］周訓(xùn)強(qiáng).旋轉(zhuǎn)機(jī)械軸心軌跡的提純、特征提取與自動識別研究［D］.重慶：重慶大學(xué)，2010.

［12］曾禹村.信號與系統(tǒng)［M］.北京：北京理工大學(xué)出版社，2010.