榮明星

（黑龍江科技學院電氣與信息工程學院，黑龍江哈爾濱 150027）

0 引言

隨著現(xiàn)代工業(yè)化進程的發(fā)展，電動機已成為了現(xiàn)代工業(yè)的主要動力設(shè)備，如果電動機因故障而停運將會對生產(chǎn)以至人身安全造成嚴重的影響。因此有必要建立一套成熟的故障診斷系統(tǒng)。電動機的振動信號最能全面的反應(yīng)電動機的運行狀態(tài)［1-2］。由于電動機振動信號是非平穩(wěn)隨機信號［3-4］，往往含有大量時變短時突發(fā)性質(zhì)的成分，傳統(tǒng)的信號分析方法已不能有效地提取電動機的故障特征，還可能將含有故障信息的微弱信號作為噪聲去除。所以采用小波包變換將非平穩(wěn)信號的各頻率成分區(qū)分開，小波包變換具有優(yōu)良的時頻局部特性和多分辨率處理能力，能準確的捕捉到非平穩(wěn)局部突變信號中的故障特征。神經(jīng)網(wǎng)絡(luò)所具備的自適應(yīng)，非線性映射能力等優(yōu)點［5］，使它在故障診斷領(lǐng)域有很大的應(yīng)用空間，因此可以用神經(jīng)網(wǎng)絡(luò)模型來進行故障分類識別。

1 電動機故障診斷方法

電動機故障診斷一般分為信號檢測、特征提取、狀態(tài)識別、診斷決策4個步驟。而使用傳感器從被診斷的設(shè)備中獲得原始信息是診斷的第一步，根據(jù)實驗條件設(shè)計適當?shù)碾妱訖C故障診斷方案。本文研究電動機轉(zhuǎn)子在正常，轉(zhuǎn)子不對中，軸承碰磨狀態(tài)下對應(yīng)的電動機振動信號，提取有價值的實驗數(shù)據(jù)。以這些數(shù)據(jù)作為神經(jīng)網(wǎng)絡(luò)的輸入，建立適合的BP神經(jīng)網(wǎng)絡(luò)并在MATLAB上進行仿真，從而研究小波變換和神經(jīng)網(wǎng)絡(luò)在電動機故障診斷領(lǐng)域中的有效結(jié)合及其應(yīng)用前景。

2 小波包分析與故障特征提取

2.1 小波變換基本分析［6］

小波是一種持續(xù)時間很短的波，其必須滿足一定容許條件。同傅里葉分析不同，小波分析的基不是唯一存在的，所有滿足小波基條件的函數(shù)都可以作為小波函數(shù)，所以要依據(jù)小波容許性條件確定基本小波。

設(shè)ψ（t）∈L2（R），其傅里葉變換為ψ（ω），如果滿足:

則稱ψ（t）為母小波或基本小波。將ψ（t）進行某種伸縮和平移，就得到:

式中:ψaτ（t）為連續(xù)小波基函數(shù);a為尺度因子;τ為位移因子。

然后，給定函數(shù)f（t）小波變換為:

實際運用中，信號要經(jīng)過計算機分析、處理，必須將信號f（t）變?yōu)殡x散時間序列。對連續(xù)小波變換而言，含有的尺度因子和平移因子是連續(xù)的，也必須離散，就有了離散小波變換DWT。將尺度因子進行冪級數(shù)離散化，令a=，a0＞0，j∈Z，位移因子離散化，即，其中，ao和τ0為大于零的實常數(shù)，j和k為整數(shù)，則得離散小波基函數(shù)為:

通常取常數(shù)a0=2，τ0=1，所以小波基函數(shù)簡化為

給定函數(shù)f（t）離散小波變換可以表示為:

2.2 小波包提取故障特征信號

小波包分析不僅能對尺度空間即低頻部分進行分解，還能對小波空間即高頻部分進行分解，有很大的優(yōu)越性，所以得到了廣泛運用。小波包分解可以將原信號分解為在不同頻帶上的投影，對這些頻帶內(nèi)的信號進行分析稱之為頻帶分析技術(shù)［7］。小波包分析能夠?qū)㈩l帶進行多層次劃分并能根據(jù)被分析信號的特征，自適應(yīng)的選擇相應(yīng)的頻帶，使之與信號頻譜相匹配，從而提高時頻分辨率，使故障特征提取能在更加細化的頻帶內(nèi)進行。圖1以3層小波包分解說明小波包分解過程。

圖1 三層小波包分解原理圖

圖1中（0，0）代表原始信號S，（1，0）代表小波包分解的第一層低頻系數(shù)，（1，1）代表小波包分解的第一層高頻系數(shù)，（3，0）代表第0個節(jié)點系數(shù)，其他均類推。

利用小波包對信號進行重構(gòu)可以根據(jù)需要把干擾和噪聲頻段清零進行信號重構(gòu)。因而只要分解時將信號中的特征信號和噪聲分解到不同頻段上，就可以得到去除干擾噪聲的重構(gòu)特征信號［8］。下面給出小波包特征向量的構(gòu)造過程［9］:

1）任給信號f（t）∈L2（R）的小波包系數(shù)可通過式（3）得到

4）由于小波變換后的ejn得總能量較大，不利于模式分類，對ejn進行如下運算:

構(gòu)造特征向量，當 j=3 時，T=［E30，E31，E32，E33，E34，E35，E36，E37］在仿真分析中，給出了故障特征提取的具體數(shù)據(jù)，作為神經(jīng)網(wǎng)絡(luò)的訓練數(shù)據(jù)和測試數(shù)據(jù)。

3 神經(jīng)網(wǎng)絡(luò)與故障識別

3.1 BP 神經(jīng)網(wǎng)絡(luò)

神經(jīng)網(wǎng)絡(luò)是在現(xiàn)代神經(jīng)學研究成果的基礎(chǔ)上提出的模擬人腦某些特性由大量人工神經(jīng)元組成的并行分布式網(wǎng)絡(luò)，具有存儲經(jīng)驗知識和使之可用的特性?；贐P算法的神經(jīng)網(wǎng)絡(luò)是目前應(yīng)用最廣泛的模型屬于前饋型網(wǎng)絡(luò)，又稱誤差反傳算法網(wǎng)絡(luò)，包括輸入層、隱含層和輸出層。它的學習規(guī)則是通過反向傳播來調(diào)整網(wǎng)絡(luò)的權(quán)值和閥值使網(wǎng)絡(luò)誤差的平方和最小。包括正向傳播與反向傳播，在正向傳播中學習樣本送入輸入層，經(jīng)隱含層運算后傳至輸出層，如果在輸出層沒有得到期望結(jié)果，則計算輸出層的誤差變化值，然后進入誤差反向傳播階段，誤差信號沿著原來的連接從輸出層返回輸入層并逐層調(diào)整連接權(quán)值，使誤差達到最小。

以電動機為診斷對象，而故障原因和狀態(tài)之間沒有明確的線性關(guān)系，不利于用精確的數(shù)學模型來描述，可以將電動機的故障狀態(tài)數(shù)據(jù)作為樣本特征向量輸入到神經(jīng)網(wǎng)絡(luò)而且給定期望輸出。神經(jīng)網(wǎng)絡(luò)通過修改神經(jīng)元連接權(quán)值進行訓練學習，達到允許的誤差范圍后停止訓練，就可以將其他故障特征向量輸入到訓練好的神經(jīng)網(wǎng)絡(luò)，從而實現(xiàn)故障的識別。

3.2 BP神經(jīng)網(wǎng)絡(luò)設(shè)計

利用BP神經(jīng)網(wǎng)絡(luò)設(shè)計一個故障識別器，識別電動機的運行狀態(tài)。第一步，以小波包分解提取的電動機狀態(tài)信號的特征向量作為輸入樣本，如下:

給定相應(yīng)的期望輸出為Xid，用兩個輸出神經(jīng)元結(jié)點，神經(jīng)元輸出（0，1）表示轉(zhuǎn)子正常，（1，0）表示轉(zhuǎn)子不對中，（1，1）表示軸承碰磨;第二步，根據(jù)Kolmogorov定理確定用三層網(wǎng)絡(luò)結(jié)構(gòu)，輸入層，隱含層和輸出層的神經(jīng)元結(jié)點數(shù)分別為N，2N+1和M，其中N=8，M=2。中間層神經(jīng)元選擇S型正切函數(shù)為激活函數(shù)，輸出層神經(jīng)元選擇S型對數(shù)函數(shù)為激活函數(shù);第三步，編寫B(tài)P神經(jīng)網(wǎng)絡(luò)程序進行訓練，直至訓練結(jié)果能夠滿足要求，選擇新的測試樣本作為訓練完畢的神經(jīng)網(wǎng)絡(luò)的輸入，測試網(wǎng)絡(luò)性能。過程流程如圖2 所示［10］。

圖2 BP神經(jīng)網(wǎng)絡(luò)算法框圖

4 Matlab軟件仿真實現(xiàn)

用加速度傳感器采集電動機原始振動信號，得到6組樣本信號數(shù)據(jù)，如圖3所示。

圖3 樣本信號

將上述信號經(jīng)過A/D轉(zhuǎn)換和離散化處理后，按特征信號提取步驟利用Matlab編寫小波包信號分析程序，對上述樣本信號進行特征提取，得到6組樣本信號特征向量，如表1。

表1 樣本信號特征向量

按照BP神經(jīng)網(wǎng)絡(luò)設(shè)計步驟，利用Matlab編寫B(tài)P網(wǎng)絡(luò)程序，輸入樣本信號特征向量和期望輸出訓練神經(jīng)網(wǎng)絡(luò)，訓練完畢后期望輸出和實際輸出如表2。

表2 期望輸出和實際輸出

在誤差允許范圍內(nèi)期望輸出和實際輸出基本一致，表明該網(wǎng)絡(luò)組成的故障狀態(tài)識別器訓練成功。再同樣利用加速度傳感器采集另外3組測試信號，經(jīng)過A/D轉(zhuǎn)換和離散化后的測試信號如圖4所示，用于測試訓練完畢的神經(jīng)網(wǎng)絡(luò)性能。

圖4 測試信號

利用小波包對測試信號進行分析，提取的特征向量如表3所示。

表3 測試信號特征向量

將測試的3組向量輸入到神經(jīng)網(wǎng)絡(luò)進行故障狀態(tài)識別，測試結(jié)果如表4所示。

表4 測試輸出

測試結(jié)果符合信號對應(yīng)的狀態(tài)，證明了利用小波包分析能夠有效的提取故障特征向量，并用經(jīng)過訓練的BP神經(jīng)網(wǎng)絡(luò)能夠?qū)﹄姍C故障狀態(tài)做出準確的診斷。

6 結(jié)語

本文將小波包分析和BP神經(jīng)網(wǎng)絡(luò)結(jié)合起來用于電動機的故障診斷。該方法可以有效的用于電動機的故障診斷，利用小波包分析故障信號，克服了傅里葉變換分析信號對故障信號中的微弱和奇異信號成分難以提取的缺點，進一步用神經(jīng)網(wǎng)絡(luò)的訓練方法進行故障狀態(tài)識別，實驗結(jié)果證明了該方法的準確性和有效性，可以實際應(yīng)用到電動機的故障診斷系統(tǒng)研究應(yīng)用中，從而對電動機的維護和維修工作提供及時的信息，節(jié)省大量的人力和資金。

［1］彭文季，羅興鑄，郭鵬程.基于第2代小波的水電機組振動信號預(yù)處理［J］.中國電機工程學報，2007，27（30）:25-29.

［2］顏秋容，劉欣，尹建國.基于小波理論的電力變壓器振動信號特征研究［J］.高電壓技術(shù)，2007，33（l）:4-8.

［3］馮志鵬，李學軍，褚福磊.基于平穩(wěn)小波包分解的水輪機非平穩(wěn)振動信號希爾伯特譜分析［J］.中國電機工程學報，2006，26（12）:32-37.

［4］鐘秉林，黃仁.機械故障診斷學［M］.2版.北京:機械工業(yè)出版社，2002.

［5］王仲生.智能故障診斷與容錯控制［M］.北京:國防工業(yè)出版社，2005.

［6］劉明才.小波分析及其應(yīng)用［M］.北京:清華大學出版社，2005.

［7］莫琦，陳立定，馮太合.小波包在旋轉(zhuǎn)機械故障振動信號處理中的應(yīng)用［J］.計算技術(shù)與自動化，2005，24（2）:14-17.

［8］廖傳軍，羅曉莉，李學軍.小波包在聲發(fā)射信號特征提取中的應(yīng)用［J］.電子測量與儀器學報，2008，22（4）.

［9］錢華明，馬吉臣，施麗娟.基于小波神經(jīng)網(wǎng)絡(luò)的導航傳感器故障診斷［J］.電子測量與儀器學報，2007，21（增刊）.

［10］王平.基于小波神經(jīng)網(wǎng)絡(luò)的電機軸承故障診斷［D］.山西:太原理工大學，2005.