柴 凱，張梅軍，黃 杰，陳 灝

（解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院，江蘇 南京 210007）

0 引言

隨著科技的飛速發(fā)展，現(xiàn)代機(jī)械對旋轉(zhuǎn)機(jī)構(gòu)的效率、性能、轉(zhuǎn)速和安全可靠性等多個方面提出越來越高的要求。當(dāng)轉(zhuǎn)子發(fā)生不對中、碰摩和油膜渦動、質(zhì)量不平衡等故障時，往往會導(dǎo)致整個機(jī)器的故障，輕則機(jī)器停止運(yùn)轉(zhuǎn)，造成經(jīng)濟(jì)損失，重則對人民生命安全造成嚴(yán)重威脅［1］。目前針對轉(zhuǎn)子故障診斷的研究方法分三類：基于知識的方法、基于解析模型的方法和基于信號處理的方法。人工智能、專家系統(tǒng)、人工神經(jīng)網(wǎng)絡(luò)和支持向量機(jī)等智能診斷也推動了轉(zhuǎn)子故障診斷技術(shù)向自動化、智能化方向發(fā)展［2］。奇異值分解（SVD）是線性代數(shù)中一種重要的矩陣分解，是矩陣分析中正規(guī)矩陣酉對角化的推廣，目前在濾波器設(shè)計、信號消噪、特征提取和弱信號分離、小波和小波包變換后續(xù)處理中得到了有效地應(yīng)用［3］。而利用支持向量機(jī)（SVM）對轉(zhuǎn)子故障進(jìn)行診斷，不僅能解決傳統(tǒng)診斷對操作人員專業(yè)知識要求高的問題，而且相比其他智能診斷方法，它在解決小樣本、非線性和高維模式識別中具有獨(dú)到的優(yōu)勢［4］。本文通過小波分解獲得不同故障下的分支信號，利用奇異值分解得到轉(zhuǎn)子的故障特征值，研究支持向量機(jī)在不同核函數(shù)和參數(shù)時對轉(zhuǎn)子不平衡、不對中、碰摩和油膜渦動等故障進(jìn)行有效識別的方法。

1 基于SVD和SVM的轉(zhuǎn)子故障診斷實驗平臺

研究采用的轉(zhuǎn)子實驗平臺示意圖及測點(diǎn)布置如圖1 所示。

圖1 轉(zhuǎn)子實驗平臺示意圖與測點(diǎn)布置

2 振動信號的小波分解

故障特征的數(shù)目是由小波分解的層數(shù)來確定的，本次模擬正常以及不平衡、不對中、碰摩和油膜渦動4種故障。圖2 為奇異值分解濾波后的轉(zhuǎn)子不對中信號時域圖，圖3 為轉(zhuǎn)子不對中信號的4層小波分解圖，a4，d4，d3，d2，d1分別代表原始信號中從低頻到高頻各個階段的分支信號，作為奇異值分解矩陣的各個行向量。經(jīng)過奇異值分解后得到5個奇異值，分別賦給A1，A2，A3，A4，A5，成為原始信號中每個頻率范圍故障信息的代表值。

3 振動信號的奇異值分解

本次實驗中，支持向量機(jī)訓(xùn)練樣本數(shù)為100個，每個信號經(jīng)過小波分解后得到5個奇異值S1，S2，S3，S4，S5，故共有100×5個需提取的奇異值。每種故障的4個模擬實驗信號奇異值見表1。由表1可知每個故障類別的奇異值有獨(dú)特的特征，這正是奇異值分解技術(shù)能用于提取故障診斷特征值的關(guān)鍵。由于電渦流傳感器所測位置不同，每種故障的4個特征值中，第1和第3行的特征值相似，第2和第4行相似。其中單數(shù)行是軸向數(shù)據(jù)，偶數(shù)行是徑向數(shù)據(jù)。各種故障類型中，軸向振動沒有徑向振動強(qiáng)烈。

圖2 奇異值分解濾波后的轉(zhuǎn)子不對中信號時域圖

圖3 轉(zhuǎn)子不對中信號4層小波分解

4 基于支持向量機(jī)的故障分類和診斷

4.1 支持向量機(jī)不同核函數(shù)對診斷結(jié)果的影響

支持向量機(jī)常用的核函數(shù)有線性核函數(shù)、多項式核函數(shù)、RBF核函數(shù)。本次用于訓(xùn)練的樣本數(shù)有5×20共100個，測試樣本數(shù)為5×10共50個。正常、不平衡、不對中、碰摩和油膜渦動五種狀態(tài)標(biāo)簽分別為1，2，3，4，5，通過采用各種核函數(shù)來比較預(yù)測進(jìn)度以及完成整個過程的速度，從而找出何種核函數(shù)最佳，測試結(jié)果如圖4 、圖5 、圖6 所示。

圖4 為線性核函數(shù)的診斷分類圖，有兩個正常樣本被誤判成碰摩，一個不平衡樣本卻判成了正常，兩個碰摩樣本誤判為正常，測試精度為90%（45／50），訓(xùn)練和預(yù)測所用時間為2.510 580s。

表1 各種故障類別的奇異值

圖4 線性核函數(shù)的診斷分類圖

圖5 多項式核函數(shù)的診斷分類圖

圖5 為多項式核函數(shù)的診斷分析圖，相比線性核函數(shù)，多項式核函數(shù)對轉(zhuǎn)子故障診斷精度達(dá)到100%，診斷精度明顯提升，但訓(xùn)練和預(yù)測所用時間為27.354 179s，診斷速率明顯下降。

圖6 為RBF核函數(shù)的診斷分析圖，測試精度為100%，訓(xùn)練和預(yù)測所用時間為0.007 853s，相比前兩者的診斷精度和診斷速度，RBF函數(shù)有更佳的綜合性能。

通過上面比較可知：線性核函數(shù)的分類準(zhǔn)確率并未達(dá)到100%，說明轉(zhuǎn)子故障屬于非線性問題。采用多項式核函數(shù)和徑向基核函數(shù)所得分類準(zhǔn)確率相同，但是考慮到多項式核函數(shù)有多項式次數(shù)、常數(shù)項、系數(shù)3個參數(shù)，而徑向基核函數(shù)只有1個與高斯函數(shù)寬度相關(guān)的參數(shù)，顯然較少的參數(shù)數(shù)量使得在參數(shù)選擇或參數(shù)尋優(yōu)上需要的計算量要小得多，所以在同樣的分類性能上，本次設(shè)計選擇RBF核函數(shù)。

圖6 RBF核函數(shù)的診斷分類圖

4.2 基于支持向量機(jī)測試樣本診斷研究

選定支持向量機(jī)各項參數(shù)后，就能對所設(shè)計的診斷系統(tǒng)進(jìn)行檢測和驗證。在本次實驗中，運(yùn)用100個樣本對構(gòu)建的支持向量機(jī)模型進(jìn)行訓(xùn)練后，隨機(jī)選取7個采集的未知故障實驗?zāi)M信號來進(jìn)行測試，其實驗結(jié)果見圖7 。

驗證故障信號組成為：1個正常信號、2個不平衡信號、1個不對中、1個碰摩和2個油膜渦動，均診斷正確。由實驗結(jié)果可知，基于SVD和SVM的轉(zhuǎn)子故障診斷系統(tǒng)對驗證信號進(jìn)行了準(zhǔn)確的分類。

5 結(jié)論

本文提出了基于SVD和SVM的轉(zhuǎn)子故障診斷方法，得到結(jié)論如下：

（1）對幾種故障類型的轉(zhuǎn)子信號分別做4層小波分解，獲得從低頻到高頻各個階段的分支信號，作為奇異值分解矩陣的行向量。經(jīng)奇異值分解得到奇異值，結(jié)合支持向量機(jī)實現(xiàn)分類診斷。

（2）通過對支持向量機(jī)不同核函數(shù)對診斷正確率和訓(xùn)練和預(yù)測所用時間的影響研究，獲得了診斷轉(zhuǎn)子故障的最優(yōu)核函數(shù)，實現(xiàn)了對未知故障實驗?zāi)M信號的測試和準(zhǔn)備分類，為轉(zhuǎn)子故障分析診斷提供了新的有效途徑。

圖7 測試樣本診斷分類圖

［1］殷蘇民，董曉明.電機(jī)轉(zhuǎn)子動平衡機(jī)及其去重策略的研究［J］.機(jī)械設(shè)計與制造，2010（12）：251-252.

［2］王凱，張永祥，李軍等.基于支持向量機(jī)的齒輪故障診斷方法研究［J］.振動與沖擊，2006，25（6）：97-99.

［3］王計生，喻俊馨，黃惟公.小波包分析和支持向量機(jī)在刀具故障診斷中的應(yīng)用［J］.振動、測試與診斷，2008，28（3）：273-276.

［4］于德介，楊宇，程軍圣，等.一種基于SVM和EMD的齒輪故障診斷方法［J］.機(jī)械工程學(xué)報，2005，41（1）：140-144.

［5］袁勝發(fā)，褚福磊.支持向量機(jī)及其在機(jī)械設(shè)備故障診斷中的應(yīng)用［J］.振動與沖擊，2007，26（11）：29-34.