駱東松 薛 鑫
(蘭州理工大學電氣工程與信息工程學院 蘭州 730050)
運行可靠始終是旋轉(zhuǎn)設(shè)備的第一安全底線,也是企業(yè)發(fā)展的根本保障。齒輪箱作為旋轉(zhuǎn)設(shè)備中較復雜的結(jié)構(gòu)部件,工作運行環(huán)境復雜,故障種類繁多,在服役期內(nèi),受徑向和軸向負荷突變,引起嚴重損傷,直接影響設(shè)備的安全運行,甚者引發(fā)安全事故。因此,設(shè)備的在線檢測和故障判斷對齒輪箱的維護和或更換具有重大的變被動為主動的戰(zhàn)略意義。軸承箱作為復雜機械部件,故障數(shù)據(jù)獲取存在困難,且不同設(shè)備之間故障數(shù)據(jù)存在較大差異,導致常用的故障監(jiān)測方法因故障數(shù)據(jù)的缺乏以及系統(tǒng)開發(fā)的困難程度,成為制約齒輪箱設(shè)備狀態(tài)監(jiān)測的一個重要原因。
支持向量數(shù)據(jù)描述(SVDD)是對支持向量機(SVM)的一種延伸。與傳統(tǒng)的SVM 算法不同,SVDD 是一種單值分類方法,就是可以對某一類數(shù)據(jù)進行分類,SVDD 的支持向量決定的最小超球體能夠包含被分類數(shù)據(jù)中的絕大部分,對樣本依賴度小,不會出現(xiàn)過擬合,具有計算速度快、魯棒性強、適合小樣本數(shù)據(jù)的處理,更適合在嵌入式設(shè)備中的應用,同時在文本分類、垃圾郵件識別、圖像分類、蛋白質(zhì)分類中被廣泛應用。
本文提出了一種經(jīng)PSO 優(yōu)化后的SVDD 算法,該算法在本文用來對齒輪箱故障進行診斷。利用經(jīng)典的振動信號處理方法,對不同故障狀態(tài)下齒輪箱振動數(shù)據(jù)進行降噪與特征參數(shù)的提取,用SVDD算法將處理后的振動數(shù)據(jù)先求取支持向量,然后利用支持向量求取最小超球體,進而求取特征參數(shù)所存在的特征域,構(gòu)造能夠表征羅茨風機齒輪箱在不同故障下運行狀態(tài)的特征域,對相關(guān)結(jié)果進一步分析總結(jié)出本文提出的方法是可行的[1]。
SVDD 早期首先是由TAX 等人結(jié)合了SVM 思想提出來的,是低維空間向高維空間的延伸。在機械故障診斷、圖像檢測等方面得到了廣泛應用,因為有了核函數(shù)的引入,導致SVDD 的性能將會直接受到核參數(shù)選擇的影響,要想達到理想的效果只有選擇較為合適的核參數(shù)才能實現(xiàn)[14]。SVDD是一種單值分類算法,可以對單組數(shù)據(jù)進行空間向量描述。
最小超球體示意見圖1,其中:a 是最小超球體的球心,R 是最小超球體的半徑。假設(shè)樣本集空間中包含N 個樣本xi(i=1,2,…,N),在實際工程當中,由于外界噪音的干擾,在構(gòu)造超求體時,會將一些奇異值包含在內(nèi),使得分類準確率下降或構(gòu)造的超球體并非為最小超球體。為減少奇異值的影響,引入松弛變量ξi[2],對應優(yōu)化問題的數(shù)學描述為
圖1 二維空間超球體示意
式(1)中:C 為懲罰因子,用來平衡目標樣本和非目標樣本的比例。
為方便求解,引入拉格朗日乘子。
式(2)中:αi≥0,βi≥0 為引入的拉格朗日乘子。將LP(a,R,α,β)由L代替,則式(2)的對偶形式為
對該優(yōu)化問題進行求解可以得出αi的值,當αi≠0 時,其對應的樣本即為支持向量[2],超球體的球心與半徑的求解公式為
對于任何一個測試樣本點z,如果滿足:
通常情況下,在低維空間中數(shù)據(jù)的分布不會呈球狀,故引入徑向基核函數(shù)K(xi·xi)代替式(3)中的內(nèi)積運算,將低維空間中的非線性問題轉(zhuǎn)化為高維空間中的線性問題[15],求解過程變得簡單。即
測試樣本z為目標樣本只需滿足:
粒子群優(yōu)化算法(Particle Swarm optimization,PSO)又翻譯為粒子群算法。PSO 初始化為一群隨機粒子(隨機解),然后通過迭代找到最優(yōu)解。假設(shè)n 維空間中第i 個粒子的位置和速度分別為Si=(Si1,Si2,…Sid) 和Vi=(Vi1,Vi2,…Vid) ,為了找到全局最優(yōu)解,算法進行迭代計算,每次迭代計算中,粒子通過跟蹤個體最優(yōu)解Pbest,Pi=(Pi1,Pi2,…Pid)與全局最優(yōu)解(gbest)來更新自己[14]。
粒子通過式(8)更新自己的速度和位置:
其中:ω為慣性權(quán)重因子,c1和c2為大于0 的學習因子,r1和r2為隨機數(shù),其分布區(qū)間為0~1。
在PSO 算法中ω是最重要的可調(diào)參數(shù),t和t+1 分別代表迭代次數(shù)當前值和下次值。由于慣性權(quán)重越大越有利于PSO算法跳出局部極小點[14],防止算法進入局部循環(huán),提高算法算法收斂速率,采用可變化的權(quán)重因子,權(quán)重因子ω在定區(qū)間內(nèi)進行變化,以適應不同時期權(quán)重因子的最優(yōu)值。ω隨算法迭代次數(shù)的變化公式為
其中:ωmax-ωmin表示權(quán)重因子的取值空間,t代表當前迭代次數(shù),tmax代表最大迭代次數(shù)。
優(yōu)化PSO 算法的參數(shù)組[K,α],將其初始種群數(shù)量初始化為10,最大迭代次數(shù)初始化為30,學習因子c1和c2都初始化為1.5,根據(jù)圖2 可以看出PSO 在尋求最優(yōu)解的過程中,局部極小值EP隨種群進化,當?shù)螖?shù)為4時,算法收斂。
圖2 PSO種群進化
針對某熱源廠大型羅茨風機減速齒輪箱故障診斷中存在的響應時間長,響應精度低等問題。本文提出將粒子群算法PSO 與支持向量數(shù)據(jù)描述SVDD 相結(jié)合,通過兩種算法的結(jié)合,發(fā)揮各算法的優(yōu)勢,而后對齒輪箱進行故障診斷。故障診斷過程如圖3所示。
圖3 故障診斷過程
步驟1:對廠區(qū)羅茨風機運行振動數(shù)據(jù)進行采集并進行特征提取,樣本數(shù)據(jù)集包括:正常樣本數(shù)據(jù)集D1、點蝕樣本數(shù)據(jù)集D2、點蝕磨損樣本數(shù)據(jù)集D3、斷齒樣本數(shù)據(jù)集D4、磨損樣本數(shù)據(jù)集D5、斷齒磨損樣本數(shù)據(jù)集D6。
步驟2:利用優(yōu)化后的PSO-SVDD 算法對樣本數(shù)據(jù)所包含的正常樣本和已知故障類型樣本數(shù)據(jù)構(gòu)建超球面,用做后面兩步故障類型判斷的依據(jù)。
步驟3:利用數(shù)據(jù)集D1構(gòu)建的故障診斷模型對測試樣本進行第一輪診斷,判別式為函數(shù)式(7)。如果滿足則說明羅茨風機運行正常;反之,則說明樣羅茨風機存在故障。
步驟4:第一輪診斷后部分測試數(shù)據(jù)顯示該齒輪箱存在故障,對存在故障信號的測試數(shù)據(jù)進行第二輪診斷。分別將測試數(shù)據(jù)用數(shù)據(jù)集D2構(gòu)建的故障診斷模型進行故障診斷,判別式依然為函數(shù)式(7),若輸出有故障,則實際數(shù)據(jù)并無該故障。同理,依次用不同故障樣本數(shù)據(jù)構(gòu)建的超球面對測試數(shù)據(jù)進行故障診斷,當輸出顯示無故障,則測試數(shù)據(jù)存在此類故障,并進行故障顯示。
本次設(shè)計應用于某熱源廠2#羅茨風機的變速齒輪箱中,模型訓練數(shù)據(jù)集通過設(shè)計的下位數(shù)據(jù)采集器采集得來,設(shè)備額定功率為75KW,額定轉(zhuǎn)速為800r/min,驅(qū)動端軸承型號為SKF6205 數(shù)據(jù)集采集轉(zhuǎn)速為額定轉(zhuǎn)速的70%,即羅茨風機的變頻器頻率為35Hz,振動信號采樣頻率為12KHz,16 位AD采樣,采樣時長8s,將采集到的數(shù)據(jù)集劃分為60組,其中40 組作為訓練數(shù)據(jù),取SVDD 算法的拒絕率為百分之一,徑向基核函數(shù)的σ參數(shù)為5,進而求出支持向量、球心和半徑。其余20 組作為測試數(shù)據(jù)。齒輪箱狀態(tài)包括正常、點蝕、點蝕磨損、斷齒、磨損以及斷齒磨損。經(jīng)PSO-SVDD診斷后得到的故障診斷結(jié)果如圖4所示。
圖4 PSO-SVDD算法診斷結(jié)果
傳統(tǒng)的SVDD 算法,作為典型的單值分類器,其缺點也是顯而易見,通過該模型只能簡單判斷出設(shè)備在運行過程中是否出現(xiàn)故障,但對于出現(xiàn)故障的具體類型卻不能夠做出判斷,即使如此簡單也已經(jīng)對事故的預防起到了重大的作用。本文提出的通過將單一故障樣本假設(shè)為正常樣本進行的第二次故障診斷,通過邏輯運算,可以很好地解決SVDD 算法因其單值分類這一特性所造成的無法對機器的故障類型做出更準確的判斷這一缺憾,完成了齒輪箱相應的故障診斷,此法亦可用于其他部件故障診斷中。
齒輪箱的故障診斷在大型風機系統(tǒng)中是不可或缺的一環(huán),但僅憑單一的傳感器及傳統(tǒng)的齒輪箱故障診斷方法很難準確診斷出具體故障類型,影響設(shè)備安全運行以及企業(yè)的經(jīng)濟效益。
本文提出的利用PSO-SVDD 單值分類器對采集到的振動數(shù)據(jù)進行兩輪診斷,通過對結(jié)果的邏輯運算,分析出齒輪箱故障。在用SVDD 模型進行故障診斷時,核函數(shù)的選擇是重中之重。本文選取的徑向基核函數(shù)相比于其他核函來說較為簡單,編程難度較低。試驗結(jié)果表明,與傳統(tǒng)的旋轉(zhuǎn)設(shè)備故障診斷技術(shù)相比,該診斷算法更加簡單,對硬件要求更加簡單,適合便攜式故障監(jiān)測設(shè)備開發(fā)所需要的智能診斷算法。
本課題由于傳感器安裝位置、傳感器誤差和特征提取誤差導致的數(shù)據(jù)間的沖突忽略不計。