曹 飛,，王 瑾,，郭前崗，周西峰
?。暇┼]電大學 自動化學院,，江蘇 南京 210023）
摘 要： 針對現(xiàn)有機車車輪超聲檢測系統(tǒng)無法準確區(qū)分其故障類型,，探討了一種基于小波包變換與BP神經(jīng)網(wǎng)絡相結合的方法來識別基于超聲檢測機車車輪的故障類型。該方法對機車車輪的超聲檢測的回波信號進行小波包分解,，并提取了頻譜能量特征值,，由此構造不同類型機車車輪內(nèi)部傷損對應的特征向量，然后利用改進的BP神經(jīng)網(wǎng)絡算法訓練隨機抽取的樣本空間,，并對測試樣本數(shù)據(jù)進行了識別驗證,。通過實驗驗證了該方法能夠快速準確地識別出機車車輪不同類型的故障對應的超聲檢測信號。
關鍵詞： 機車車輪,；超聲檢測,；小波包；BP神經(jīng)網(wǎng)絡,；故障診斷
0 引言
隨著我國鐵路建設的高速發(fā)展,，越來越多的機車車輛進入了現(xiàn)有高速鐵路系統(tǒng)。機車車輪是列車車輛行車的重要部件,。目前針對列車車輪的無損檢測的方法和設備種類很多,，超聲檢測是國內(nèi)外使用最為廣泛、使用頻率最高的一項無損檢測技術[1],。目前機車車輪的超聲檢測中依然多以A型探傷儀人工掃查為主,，檢測結果多憑借操作人員的個人經(jīng)驗，摻雜了較多的主觀因素[2],。本文主要是利用小波包分解對機車車輪的超聲檢測結果進行分析,，然后結合神經(jīng)網(wǎng)絡原理對其故障類型歸類，進而實現(xiàn)對故障類型自動識別,。
根據(jù)超聲探傷原理,，判別材料內(nèi)部缺陷的主要依據(jù)來源于超聲檢測的回波信號，其回波信號中包含了時間（聲程）和超聲波聲壓信息,，一般根據(jù)經(jīng)驗并結合被測材料的性質(zhì)來分析和估計材料的內(nèi)部情況,。而利用小波包多分辨分析結合頻譜能量的方法可以提取到超聲檢測過程中回波信號的特征值，神經(jīng)網(wǎng)絡又具有優(yōu)秀的分類識別能力,，使用一個隱含層就能正確完成超聲信號的分類識別,，兩者的結合能夠很好地解決缺陷識別的問題[3-4]。
1 超聲信號特征提取
1.1 信號小波包分解
正交小波變換主要用來對信號的低頻部分進行分解,，所以小波變換只可以表征以低頻信息為主要成分的信號,，而不能用來表示含有很多細節(jié)（高頻）的信息。而神經(jīng)網(wǎng)絡需要對缺陷信號的多個頻譜特征學習,，所以對超聲信號提取頻譜特征的前提是要保留高頻和低頻的頻率特性,，小波包分解具有同時提取信號的高頻部分和低頻部分的能力。根據(jù)定義，超聲信號可使用以下雙尺度方程進行小波包分解：
上式中hk和gk分別為低通實濾波器和高通濾波器,。對超聲檢測的回波信號進行小波包分解實際上是使用高低通濾波器組對回波信號進行分解,，每次分解的過程中都會將信號分解為高低兩個頻率，即將頻率[0,，2jπ]的信號分成低頻[0,，2（j-1）π]和高頻[2（j-1）π，2jπ]兩部分,，并對每個部分進行二次抽樣,。接著對高低頻部分做同樣的分解過程，直到滿足要求,。經(jīng)過每次分解后各頻段的數(shù)值都會減少1/2,，選擇合適的分解層數(shù)便能得到合適的頻帶[5]。小波包的分解構成了一個分解二叉樹,，本文對車輪的超聲檢測信號進行3層小波包分解的小波包樹結構如圖1所示,，信號經(jīng)過每一層都分解為高低頻兩個部分，到第三層時共產(chǎn)生8個節(jié)點,。
1.2 超聲信號的特征提取
列車車輪的超聲檢測中,，如果車輪內(nèi)部存在傷損，則其超聲檢測的回波信號的幅頻特性會出現(xiàn)不同程度的變化,，具體為某一頻率段中的超聲檢測回波信號的譜能量會增強,。各頻率成分信號的能量中均包含著超聲檢測到的內(nèi)部傷損信息,。因此結合小波包多分辨率分析與各頻段能量提取的方法,，通過“能量—缺陷”的信號特征提取可以提取列車車輪超聲檢測信號中的車輪內(nèi)部傷損特征。提取超聲信號特征向量的步驟如下[4-7]：
?。?）小波包分解：對得到的超聲檢測數(shù)據(jù)使用式（1）進行3層小波包分解,。
（2）信號重構：對經(jīng)過步驟（1）的小波包分解得到的8個頻率段內(nèi)的序列進行重構,，得到8個小波包重構信號,，重構信號中包含了超聲信號中沿超聲探頭中心頻率的從低到高的各頻率段的信息。
?。?）各節(jié)點頻譜能量計算：設S3,，i（i=0，1,，…,，7）為第三層小波包節(jié)點的重構信號，先求出各重構信號的幅度頻譜值SF（3,，i）（i=0,，1，…，7）,，然后分別求取各頻譜的向量二范數(shù),，即各頻帶對應的能量值：
其中|SF（3，i）（j）|表示各節(jié)點重構信號的單邊頻幅離散點的幅值,。
?。?）構造特征向量：將各經(jīng)過小波包分解的節(jié)點的頻譜能量特征值構成特征向量E=[E3，0,，E3,，1，E3,，2,，E3，3,，E3,，4，E3,，5,，E3，6,，E3,，7]，若頻譜能量特征值過大,，分析和使用過程將變得繁雜,，所以需要對提取的特征向量作歸一化處理，即：T3,，i=E3,，i/‖E‖2。
歸一化之后的特征向量為：T=[T3,，0,，T3，1,，T3,，2，T3,，3,，T3，4,，T3,，5，T3，6,，T3,，7]。利用上述方法對檢測的多組列車輪的超聲檢測信號分別提取了譜能量特征值,，并構造其特征向量,。經(jīng)過實驗比較后選取db6作為小波函數(shù)進行分解，其分解出不同頻率段的頻譜差異性較大,，不同頻段的頻譜能量特征比較明顯,。
2 超聲信號類型分類識別
2.1 BP神經(jīng)網(wǎng)絡
BP神經(jīng)網(wǎng)絡是一種按誤差逆?zhèn)鞑ニ惴ㄓ柧毜那梆伨W(wǎng)絡，其具有優(yōu)秀的自學習,、自適應能力[7],。BP神經(jīng)網(wǎng)絡能夠用來學習大量的輸入/輸出的模式映射關系，通過反向傳播調(diào)整網(wǎng)絡參數(shù)中的權值和閾值,，最終實現(xiàn)誤差平方和最小,。BP神經(jīng)網(wǎng)絡的拓撲結構一般包括輸入層、隱層和輸出層,，如圖2所示,。其學習過程分為正向傳播過程和反向傳播過程，其中正向傳播時,，神經(jīng)元的狀態(tài)只能影響下一層神經(jīng)元網(wǎng)絡,，神經(jīng)網(wǎng)絡在運算過程中，要先將輸入信息向前傳播到隱節(jié)點上,，再把隱節(jié)點的輸出傳播到輸出節(jié)點,，層與層之間的聯(lián)系都要經(jīng)過激活函數(shù)運算；當沒有輸出期望結果時,，說明實際輸出與期望輸出之間存在一定的誤差,，需要將誤差繼續(xù)返回反向傳播的過程,，同時修改各層的權值,，重新逐層計算。正向傳播和反向傳播不斷執(zhí)行,，直到滿足誤差要求,，訓練學習的過程才結束。
2.2 BP網(wǎng)絡參數(shù)選擇
若設某神經(jīng)網(wǎng)絡輸入為X,，輸入層含有n個神經(jīng)元,；其中隱含層的輸出為Y，并且含有q個神經(jīng)元,；輸出向量為O,，輸出層含有m個神經(jīng)元，目標矢量為T；隱層和輸出層對應的激活函數(shù)為f1和f2,，對應的權值矩陣分別為W1和W2,。則對于隱層有：
根據(jù)缺陷頻譜能量特征向量的維數(shù)，將對應的輸入節(jié)點數(shù)目設定為8個,，根據(jù)檢測輸出的結果,，共定義了3種缺陷，用兩位二進制數(shù)即可表示缺陷類型,，定義（0,，0）為車輪踏面表面裂紋；（1,，0）為車輪橫孔缺陷,；（0，1）為車輪周向裂紋缺陷[3,，8],。隱層節(jié)點個數(shù)的確定主要根據(jù)經(jīng)驗公式，結合多次實驗的結果,，本文選取如下經(jīng)驗公式確定隱層節(jié)點的個數(shù)[9]：
其中,，m為隱層節(jié)點個數(shù)，n為輸入層節(jié)點個數(shù),，l為輸出層節(jié)點數(shù)[9],。
2.3 改進的BP網(wǎng)絡算法
標準BP算法在應用過程中存在以下缺陷：（1）容易形成局部極小卻達不到全局最優(yōu)；（2）訓練次數(shù)過多,，使得學習效率低且收斂速度慢,；（3）隱節(jié)點的選取缺乏理論性的指導。
為了更準確快速地實現(xiàn)超聲信號特征值的識別,，本文在標準BP算法的基礎上引入一種改進的BP算法——彈性算法（RPROP）,。RPROP算法改進的意義是消除梯度幅度的影響，在權值修正過程中,，減少幅值對權值的影響,，權值的大小取決于與幅值無關的修正值。與其他改進的BP算法相比,，RPROP算法在模式識別中使用時收斂速度最快[10-11],。
3 實驗分析
3.1 機車車輪超聲檢測系統(tǒng)
國內(nèi)某公司自研的超聲檢測系統(tǒng)主要由超聲檢測探頭、數(shù)據(jù)采集單元和顯示控制主機部分組成,。其硬件組成和測試過程如圖3所示,。
3.2 機車車輪缺陷信號的提取
系統(tǒng)中的探頭主要使用中心頻率為2.5 MHz的雙晶縱波探頭，縱波探頭對周向缺陷較敏感,，主要用來探測機車車輪的踏面剝離,、內(nèi)部孔狀裂紋和輪對周方向裂紋等,。數(shù)據(jù)采集單元激勵和采集超聲波信號，超聲回波信號經(jīng)過高速A/D芯片采集,，并經(jīng)過濾波和壓縮處理后不失真地恢復為原信號,，從啟振點起一組回波信號由248個數(shù)據(jù)點組成。顯示控制主機部分將數(shù)據(jù)采集單元處理過的數(shù)據(jù)進行顯示和管理,。本文實驗中使用的數(shù)據(jù)是通過系統(tǒng)在機車樣板輪上選擇不同類型的缺陷點采集的多組數(shù)據(jù),。
3.3 車輪超聲檢測信號特征提取
實驗中使用機車車輪超聲探傷系統(tǒng)對樣板輪進行了掃查，樣板輪在表面和內(nèi)部畫有刻槽以模擬車輪的踏面表面剝離和車輪內(nèi)部孔狀裂紋及內(nèi)部周向裂紋等故障,，探傷系統(tǒng)上傳的采樣超聲回波信號經(jīng)過濾波和數(shù)據(jù)壓縮,，其對應的時域圖如圖4所示。根據(jù)超聲波傳播的特點,，當超聲波信號遇到內(nèi)部缺陷（不同密度的物質(zhì)）后會有回波信號反射,，從時域圖上只能分辨出缺陷故障的聲程，不太容易從直觀上區(qū)分缺陷故障的類型,。不同類型的缺陷故障對應的小波包分解頻譜圖和譜能量柱狀圖如圖5,、6所示，分解后的頻譜圖中不同缺陷的頻譜已有較大的區(qū)別,，再提取出不同的頻譜能量特征可以用來區(qū)分不同的故障類型,。將頻譜能量作為特征量用于神經(jīng)網(wǎng)絡的訓練學習。
3.4 BP神經(jīng)網(wǎng)絡自診斷
實驗中針對機車車輪的不同故障類型,，從每種類型的故障隨機選取5組實驗數(shù)據(jù)作為訓練樣本和2組數(shù)據(jù)作為測試樣本,。設定網(wǎng)絡訓練次數(shù)為1 000，目標誤差為0.000 01,。其中測試樣本的故障診斷測試結果和訓練誤差曲線分別如表1和圖7所示,。從圖中可以看出基于RPROP算法的改進的神經(jīng)網(wǎng)絡具有較快的收斂速度，只需經(jīng)過訓練29步就能夠很好地識別列車車輪超聲檢測損傷,，達到了預期結果,。
4 結束語
本文介紹了一種基于改進BP算法的神經(jīng)網(wǎng)絡在機車車輪超聲檢測中識別車輪故障類型的方法，該方法通過小波包對實測的超聲檢測信號進行分解,，提取其特征值向量,，利用改進BP算法的神經(jīng)網(wǎng)絡進行適當?shù)挠柧毢螅梢杂脕碜R別車輪內(nèi)部的傷損類型,。
超聲檢測信號的歸類和識別方法中,，目前提及較多的是提升小波分析提取特征值的方法和基于LM（Levenberg-Marquadt）算法的BP神經(jīng)網(wǎng)絡的歸類方法[8]。提升小波分析對專業(yè)要求較高,，分解過程復雜，不適合普及,；LM學習算法在處理歸類識別的有關問題時,，處理能力很弱,。本文采用小波包分解并結合基于RPROP算法的神經(jīng)網(wǎng)絡能夠有效地區(qū)分不同基于超聲檢測的機車車輪的故障類型，后期在實踐中可以將此方法嵌入超聲檢測系統(tǒng)的上位機軟件,，使系統(tǒng)具有自動識別和歸類機車車輪故障類型的能力,，此方法具有重要的實際意義。
參考文獻
[1] 宋天民,，張國福,，劉峰，等.超聲檢測[M].北京：中國石化出版社,，2012.
[2] 馬晨.基于PC機的機車輪軸超聲探傷儀實現(xiàn)[J].現(xiàn)代電子技術,，2010，33（5）：138-141.
[3] 周西峰,，索會迎,，郭前崗，等.基于小波包-BP網(wǎng)絡的超聲檢測缺陷類型識別[J].解放軍理工大學學報（自然科學版）,，2012,，13（5）：521-526.
[3] 張冬雨，劉小方,，楊劍,，等.基于小波包變換的復合材料分層缺陷信號特征分析[J].兵工自動化，2009,，28（11）：56-58.
[4] 孫延奎.小波分析及其應用[M].北京：機械工業(yè)出版社,，2005.
[5] 師小紅，敦怡,，徐章遂,，等.基于“能量-缺陷”的金屬基復合材料缺陷信號特征提取[J].計算機測量與控制，2006,，14（1）：109-110.
[6] 桂中華,，韓鳳琴.小波包特征熵神經(jīng)網(wǎng)絡在尾水管故障診斷中的應用[J].中國電機工程學報，2005,，25（4）：100-102.
[7] 候媛彬,，杜京義，汪梅.神經(jīng)網(wǎng)絡[M].西安：西安電子科技大學出版社,，2007.
[8] Li Jian,， Zhan Xianglin， Jin Shijiu. An automatic flaw classification method of ultrasonic nondestructive testing for pipeline girth welds[C]. Proceeding of the 2009 IEEE International Conference on Information and Automation,，2009：22-25.
[9] 韓立群.人工神經(jīng)網(wǎng)絡教程[M].北京：北京郵電大學出版社,，2006.
[10] 周品.MATLAB神經(jīng)網(wǎng)絡設計與應用[M].北京：清華大學出版社，2013.
[11] RIEDMILLER M,， BRAUN H. A direct adaptive method for faster backpropagation learning： the RPROP algorithm[C]. Neural Networks,， IEEE International Conference on. IEEE,， 1993： 586-591.