一種基于BPNN和SVM-PDE的旋轉(zhuǎn)機(jī)械變工況預(yù)警方法

崔錦淼， 胡明輝,2， 馮 坤,2*， 賀 雅,2， 石?；?/p>

(1.北京化工大學(xué) 發(fā)動(dòng)機(jī)健康監(jiān)控及網(wǎng)絡(luò)化教育部重點(diǎn)實(shí)驗(yàn)室，北京 100029； 2.北京化工大學(xué) 高端機(jī)械裝備健康監(jiān)控與自愈化北京市重點(diǎn)實(shí)驗(yàn)室，北京 100029； 3.中國(guó)石化銷(xiāo)售股份有限公司華南分公司，廣東 廣州 510180)

旋轉(zhuǎn)機(jī)械是工業(yè)應(yīng)用最為廣泛的機(jī)械設(shè)備之一，如電機(jī)、離心泵、離心壓縮機(jī)、燃?xì)廨啓C(jī)和航空發(fā)動(dòng)機(jī)等，涉及到國(guó)民經(jīng)濟(jì)各個(gè)發(fā)展領(lǐng)域[1]，由于在高溫、高壓、變負(fù)載環(huán)境下長(zhǎng)時(shí)間連續(xù)運(yùn)行，故障多發(fā)。為及時(shí)發(fā)現(xiàn)設(shè)備故障以避免安全事故、減小設(shè)備維修成本，需實(shí)時(shí)監(jiān)測(cè)和評(píng)估設(shè)備運(yùn)行狀態(tài)，及時(shí)對(duì)設(shè)備故障預(yù)警。

振動(dòng)監(jiān)測(cè)是目前常用的旋轉(zhuǎn)機(jī)械狀態(tài)監(jiān)測(cè)手段。傳統(tǒng)的監(jiān)測(cè)方法是專(zhuān)業(yè)技術(shù)人員結(jié)合相關(guān)振動(dòng)標(biāo)準(zhǔn)和現(xiàn)場(chǎng)經(jīng)驗(yàn)設(shè)定報(bào)警閾值實(shí)現(xiàn)設(shè)備振動(dòng)故障的預(yù)警，一般取所有工況中的最大測(cè)量值作為設(shè)備振動(dòng)烈度，以設(shè)備可長(zhǎng)期運(yùn)行的基線(xiàn)值上限25%作為報(bào)警限值，即固定報(bào)警閾值[2]。該方法雖能在一定程度上避免某些惡性故障的發(fā)生，但由于判別特征單一且未考慮設(shè)備在時(shí)變工況下運(yùn)行導(dǎo)致的振動(dòng)水平變化，對(duì)于燃?xì)廨啓C(jī)、航空發(fā)動(dòng)機(jī)等工況變化頻繁復(fù)雜的旋轉(zhuǎn)機(jī)械，容易造成高工況虛警和低工況漏警的問(wèn)題。

國(guó)內(nèi)外研究學(xué)者在設(shè)備報(bào)警閾值設(shè)定方面開(kāi)展了大量的研究。魏中青等[3]根據(jù)旋轉(zhuǎn)機(jī)械發(fā)生突變故障時(shí)故障部件特征頻率幅值突變的故障機(jī)理，基于高斯分布3σ法則對(duì)故障特征頻率幅值設(shè)定報(bào)警閾值，實(shí)現(xiàn)了旋轉(zhuǎn)機(jī)械突變故障的快速預(yù)警。程晶等[4]基于高斯分布置信度為0.95的單側(cè)置信區(qū)間分別設(shè)定設(shè)備振動(dòng)總能量和1/3倍頻程帶級(jí)能量特征的報(bào)警閾值，該閾值可對(duì)艦船設(shè)備故障進(jìn)行有效預(yù)警。李會(huì)鵬等[5]利用Bayes隨機(jī)加權(quán)法估計(jì)小樣本數(shù)據(jù)振動(dòng)有效值的均值和方差，然后基于高斯分布的3σ法則設(shè)定有效值閾值，解決了樣本數(shù)據(jù)不足時(shí)閾值難設(shè)定的問(wèn)題。然而，上述預(yù)警方法中考慮的特征較為單一，難以全面涵蓋設(shè)備各狀態(tài)振動(dòng)信息和實(shí)現(xiàn)早期故障預(yù)警。因此，Jiang等[6]基于設(shè)備正常狀態(tài)下振動(dòng)信號(hào)峰峰值、有效值等多個(gè)參數(shù)的波動(dòng)特性，設(shè)定了多個(gè)特征報(bào)警閾值，結(jié)果表明，多特征報(bào)警閾值能夠有效提高風(fēng)電齒輪箱故障預(yù)警準(zhǔn)確率。明廷鋒等[7]通過(guò)相關(guān)分析選取對(duì)故障敏感的特征參數(shù)并基于高斯分布3σ法則設(shè)定報(bào)警閾值，應(yīng)用時(shí)依據(jù)報(bào)警效果自適應(yīng)調(diào)整閾值，相較于固定閾值應(yīng)用價(jià)值更高。上述報(bào)警閾值的設(shè)定均基于高斯分布3σ法則，但對(duì)于大型復(fù)雜設(shè)備，其故障形式多樣，振動(dòng)特征往往不屬于高斯分布且特征分布復(fù)雜未知；并且分別對(duì)每個(gè)特征設(shè)定閾值，忽略了各個(gè)特征間的交互作用。

核概率密度估計(jì)法是常用的非參數(shù)估計(jì)法，可解決特征分布未知問(wèn)題，且其中的核函數(shù)可對(duì)多元參數(shù)進(jìn)行融合，因此該方法可廣泛應(yīng)用于設(shè)備預(yù)警閾值估計(jì)中。劉晗等[8]基于Parzen窗估計(jì)多元監(jiān)測(cè)參數(shù)的聯(lián)合概率密度函數(shù)以設(shè)定報(bào)警限，并隨數(shù)據(jù)積累更新調(diào)整閾值，實(shí)現(xiàn)了風(fēng)機(jī)機(jī)組故障預(yù)警。Desforges等[9]基于核概率密度估計(jì)求解數(shù)據(jù)概率密度，依據(jù)測(cè)量數(shù)據(jù)概率密度與正常數(shù)據(jù)概率密度偏離度對(duì)電機(jī)故障有效預(yù)警。Zang等[10]提出一種基于多核概率密度估計(jì)的多變量報(bào)警閾值設(shè)定方法，相較于依據(jù)3σ法則設(shè)定的閾值，有效降低了虛警率和漏警率。但上述概率密度方法中窗函數(shù)寬度難以確定，選取過(guò)大或過(guò)小都易造成較大誤差，很難應(yīng)用于大樣本數(shù)據(jù)。為了解決該問(wèn)題，張炤等[11]提出了基于支持向量機(jī)(Support Vector Machine,SVM)的概率密度估計(jì)(Probability Density Estimation,PDE)方法，該方法利用SVM求解樣本概率密度，在相同求解精度下解決了窗函數(shù)寬度選取的問(wèn)題，相較于核概率密度估計(jì)方法，該方法工程應(yīng)用價(jià)值更高，更適合于估計(jì)振動(dòng)特征分布特性設(shè)定振動(dòng)閾值。

考慮到變工況對(duì)設(shè)備振動(dòng)水平的影響，胡雷等[12]利用相關(guān)向量機(jī)擬合振動(dòng)特征與工況參數(shù)間的函數(shù)，然后設(shè)定隨工況參數(shù)變化的自適應(yīng)閾值，實(shí)現(xiàn)了設(shè)備不同工況下的故障預(yù)警。顧煜炯等[13]手動(dòng)劃分設(shè)備運(yùn)行工況區(qū)間區(qū)分設(shè)備運(yùn)行工況，然后分別設(shè)定各個(gè)運(yùn)行工況的報(bào)警閾值，實(shí)現(xiàn)了風(fēng)電齒輪箱在變工況下的故障預(yù)警。劉長(zhǎng)良等[14]則利用模糊C均值聚類(lèi)方法識(shí)別設(shè)備運(yùn)行工況，然后設(shè)定各工況的預(yù)警閾值，相較于固定報(bào)警閾值降低了虛警率。但上述變工況閾值設(shè)定過(guò)程中，當(dāng)設(shè)備工況參數(shù)個(gè)數(shù)過(guò)多、運(yùn)行工況區(qū)間過(guò)多時(shí)，上述工況識(shí)別方法難以有效擬合工況參數(shù)與工況標(biāo)記間的函數(shù)關(guān)系。BP神經(jīng)網(wǎng)絡(luò)(BP Neural Network,BPNN)具有強(qiáng)大的擬合能力，可擬合任何非線(xiàn)性函數(shù)，且網(wǎng)絡(luò)節(jié)點(diǎn)數(shù)可依據(jù)工況參數(shù)個(gè)數(shù)自適應(yīng)調(diào)整，不受限于工況參數(shù)維度與運(yùn)行工況區(qū)間個(gè)數(shù)，更具應(yīng)用價(jià)值。

因此，針對(duì)目前相關(guān)預(yù)警方法和設(shè)備運(yùn)行工況識(shí)別方法存在的問(wèn)題，提出基于BPNN和SVM-PDE的旋轉(zhuǎn)機(jī)械變工況預(yù)警方法。首先利用BPNN準(zhǔn)確識(shí)別設(shè)備運(yùn)行工況；然后提取各工況振動(dòng)信號(hào)多維特征，并基于SVM-PDE設(shè)定各工況振動(dòng)閾值進(jìn)行預(yù)警；最后利用雙轉(zhuǎn)子試驗(yàn)臺(tái)實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證所提出的方法的有效性。

1 方法理論

1.1 BPNN原理

圖1 典型三層BPNN結(jié)構(gòu)

BPNN基于BP算法反向更新權(quán)重訓(xùn)練網(wǎng)絡(luò)，原理如下[15]。

① 網(wǎng)絡(luò)前向計(jì)算。

y=f(w·x+b)

(1)

式中,“·”為矩陣乘積;x為網(wǎng)絡(luò)輸入;w為連接權(quán)重;b為偏置，可視為輸入為1的神經(jīng)元所對(duì)應(yīng)的權(quán)重；f(·)為激勵(lì)函數(shù)；y為網(wǎng)絡(luò)輸出。

② 誤差反向計(jì)算。

δi=wi+1·δi+1

(2)

式中，δi為第i層誤差項(xiàng)；wi+1為第i+1層的權(quán)重；δi+1為第i+1層誤差項(xiàng)。網(wǎng)絡(luò)最后一層的誤差計(jì)算與選取的目標(biāo)函數(shù)相關(guān)，如常用的目標(biāo)函數(shù)為最小化均方誤差：

(3)

式中，m為訓(xùn)練樣本總數(shù)；ti為目標(biāo)值；yi為網(wǎng)絡(luò)輸出；L為網(wǎng)絡(luò)誤差。則網(wǎng)絡(luò)最后一層誤差項(xiàng)δlast為

(4)

③ 權(quán)重梯度計(jì)算。

(5)

④ 權(quán)重更新。

(6)

利用BPNN識(shí)別設(shè)備運(yùn)行工況時(shí)，輸入層為表征設(shè)備運(yùn)行工況的參數(shù)，輸出層為人為給定的工況標(biāo)簽，通過(guò)BPNN擬合輸入與輸出即工況參數(shù)和工況標(biāo)簽間的函數(shù)關(guān)系，實(shí)現(xiàn)工況參數(shù)到工況標(biāo)簽的非線(xiàn)性映射，繼而實(shí)現(xiàn)多工況參數(shù)、多工況區(qū)間下的設(shè)備運(yùn)行工況自動(dòng)識(shí)別。

1.2 SVM-PDE原理

張炤、溫珂?zhèn)?、徐玉兵等人利用SVM求解線(xiàn)性算子方程以求解樣本概率密度，樣本的概率密度函數(shù)僅與樣本的邊緣即支持向量有關(guān)，相較于核概率密度估計(jì)方法計(jì)算量更小，且解決了窗函數(shù)寬度選取的問(wèn)題，該方法原理如下[11,16-17]。

利用SVM求解線(xiàn)性算子方程可轉(zhuǎn)換為對(duì)樣本概率密度進(jìn)行回歸估計(jì)，即

(7)

式中，k(·)為核函數(shù);xi為第i個(gè)樣本;xj為第j個(gè)樣本;k(xi,xj)為樣本xi、xj經(jīng)核函數(shù)投影值;βj為系數(shù)向量β對(duì)應(yīng)的第j個(gè)系數(shù)；m為樣本總數(shù)；p(xi,β)為樣本xi的概率密度。

其中常用的多項(xiàng)式核函數(shù)為

(8)

結(jié)合式(7)和式(8)可得樣本的分布函數(shù)為

(9)

(10)

為防止樣本中噪聲干擾引起過(guò)擬合，引入松弛因子ξ(ξ≥0)和不敏感損失常數(shù)ε，保證SVM擬合的概率分布F(xi,β)與經(jīng)驗(yàn)概率分布F(xi)保持一致性，則有：

(11)

因此，利用SVM求解線(xiàn)性算子方程時(shí)，SVM相應(yīng)的二次規(guī)劃問(wèn)題可轉(zhuǎn)換為

(12)

式中不敏感損失常數(shù)ε可依據(jù)經(jīng)驗(yàn)公式計(jì)算，即

(13)

通過(guò)求解式(12)中的二次規(guī)劃問(wèn)題，可解得系數(shù)向量β，結(jié)合核函數(shù)矩陣K(xi,xj)，代入式(7)即解得各樣本的概率密度值。

故障數(shù)據(jù)特征分布與正常數(shù)據(jù)特征分布存在一定的差異，即故障數(shù)據(jù)各維特征與正常數(shù)據(jù)各維特征間的距離大于各正常數(shù)據(jù)各維特征間的距離，因此結(jié)合式(7)和式(8)，基于正常數(shù)據(jù)計(jì)算的故障數(shù)據(jù)概率密度值會(huì)小于正常數(shù)據(jù)概率密度值。以正常數(shù)據(jù)概率密度下邊界作為報(bào)警限，當(dāng)待預(yù)測(cè)樣本概率密度值小于概率密度閾值時(shí)，即可視為異常，反之則視為正常，基于此可對(duì)故障預(yù)警。

1.3 特征提取

從設(shè)備振動(dòng)信號(hào)時(shí)域波形中提取出有效表征設(shè)備運(yùn)行狀態(tài)的特征參數(shù)，是實(shí)現(xiàn)設(shè)備狀態(tài)監(jiān)測(cè)、故障預(yù)警和診斷的關(guān)鍵[18]。本文共選取了21個(gè)時(shí)域、頻域統(tǒng)計(jì)指標(biāo)，計(jì)算公式如表1所示。

表1 振動(dòng)信號(hào)特征計(jì)算公式

上述振動(dòng)信號(hào)特征之間難免會(huì)存在信息重疊問(wèn)題，并且特征維度過(guò)高，不利于后續(xù)處理分析。主成分分析(Principal Component Analysis,PCA)是常用的無(wú)監(jiān)督線(xiàn)性降維方法[19]，其將高維數(shù)據(jù)向低維空間投影以實(shí)現(xiàn)降維，并保留數(shù)據(jù)的主要信息。所以可利用PCA對(duì)提取到的特征參數(shù)進(jìn)行數(shù)據(jù)降維，提取二次特征，減少數(shù)據(jù)冗余。

2 提出的方法

提出的基于BPNN和SVM-PDE的旋轉(zhuǎn)機(jī)械變工況預(yù)警方法原理如圖2所示，主要包含以下步驟。

圖2 方法原理圖

(1) 設(shè)備運(yùn)行工況識(shí)別。

確定表征設(shè)備運(yùn)行工況的參數(shù)X，并根據(jù)設(shè)備振動(dòng)水平隨運(yùn)行工況表征參數(shù)的變化，劃分設(shè)備運(yùn)行工況區(qū)間，標(biāo)定設(shè)備運(yùn)行各工況Y。將參數(shù)X和標(biāo)定工況Y劃分為訓(xùn)練集和測(cè)試集。構(gòu)建BPNN模型，利用訓(xùn)練集迭代訓(xùn)練網(wǎng)絡(luò)，測(cè)試集測(cè)試模型工況識(shí)別效果，進(jìn)而實(shí)現(xiàn)多工況參數(shù)、多工況區(qū)間下設(shè)備運(yùn)行工況的快速準(zhǔn)確識(shí)別。

(2) 振動(dòng)數(shù)據(jù)特征提取。

按照表1所示的特征計(jì)算公式，分別從振動(dòng)信號(hào)中提取各個(gè)特征指標(biāo)，并構(gòu)建出高維特征集。利用PCA對(duì)特征集降維，得到低維特征矩陣?；诓襟E①中標(biāo)定的工況，劃分低維特征矩陣，得到各工況下的低維特征矩陣，實(shí)現(xiàn)各工況振動(dòng)數(shù)據(jù)的特征提取。

(3) 變工況報(bào)警閾值設(shè)定。

基于各工況正常振動(dòng)狀態(tài)的低維特征矩陣，求解式(12)的二次規(guī)劃問(wèn)題，得到正常數(shù)據(jù)的系數(shù)向量β，進(jìn)而按式(7)和式(8)計(jì)算各工況下的正常數(shù)據(jù)概率密度值，并以概率密度值邊界即下限作為報(bào)警限。考慮到數(shù)據(jù)的波動(dòng)特性，將報(bào)警限下調(diào)3倍標(biāo)準(zhǔn)差。

(4) 報(bào)警狀態(tài)判別。

基于步驟①的工況識(shí)別模型識(shí)別待判定數(shù)據(jù)所處工況，依據(jù)式(8)計(jì)算待判定數(shù)據(jù)與該工況正常數(shù)據(jù)的核函數(shù)矩陣，依據(jù)步驟③中解得的該工況正常數(shù)據(jù)系數(shù)向量β和式(7)，計(jì)算出待判定數(shù)據(jù)的概率密度值。比較其與該工況振動(dòng)概率閾值的大小，當(dāng)該概率密度值大于閾值時(shí)，判定設(shè)備處于正常狀態(tài)；反之，判定設(shè)備處于報(bào)警狀態(tài)，基于此實(shí)現(xiàn)故障預(yù)警。

3 實(shí)驗(yàn)驗(yàn)證

3.1 雙轉(zhuǎn)子試驗(yàn)臺(tái)

本研究使用的正常數(shù)據(jù)和故障數(shù)據(jù)取自北京化工大學(xué)診斷與自愈工程研究中心的雙轉(zhuǎn)子試驗(yàn)臺(tái)。試驗(yàn)采集系統(tǒng)如圖3所示，試驗(yàn)臺(tái)主要由電主軸(外)、聯(lián)軸器、軸承座、1#軸承、2#軸承和軸承座、聯(lián)軸器、電主軸(內(nèi))組成。本實(shí)驗(yàn)通過(guò)改變與中介軸承內(nèi)外圈相連主軸的轉(zhuǎn)速模擬旋轉(zhuǎn)機(jī)械變工況運(yùn)行狀態(tài)，分別采集內(nèi)外圈轉(zhuǎn)速在0/300 r/min、300/600 r/min、600/1200 r/min、900/1500 r/min、1500/1800 r/min下中介軸承在正常狀態(tài)、內(nèi)圈故障、外圈故障、滾動(dòng)體故障狀態(tài)下的內(nèi)、外主軸轉(zhuǎn)速信號(hào)和2#軸承座處的加速度信號(hào)，用于驗(yàn)證所提出方法的有效性，其中加速度信號(hào)的采樣頻率為20480 Hz。由于本實(shí)驗(yàn)未模擬雙轉(zhuǎn)子設(shè)備隨負(fù)載變化情況，故僅以?xún)?nèi)外轉(zhuǎn)軸轉(zhuǎn)速作為工況標(biāo)定參數(shù)劃分5個(gè)運(yùn)行工況(如表2所示)，實(shí)際應(yīng)用中還應(yīng)考慮負(fù)載變化情況。

圖3 雙轉(zhuǎn)子實(shí)驗(yàn)采集系統(tǒng)

表2 工況標(biāo)定表

3.2 固定閾值預(yù)警分析

由于本實(shí)驗(yàn)主要模擬滾動(dòng)軸承故障，加速度信號(hào)對(duì)滾動(dòng)軸承故障敏感，故采用加速度信號(hào)設(shè)定報(bào)警閾值，實(shí)際應(yīng)用中應(yīng)結(jié)合設(shè)備情況選取位移、速度、加速度確定報(bào)警閾值[20]。若采用固定閾值預(yù)警法，即不考慮設(shè)備時(shí)變工況影響，各個(gè)工況下使用同一閾值進(jìn)行報(bào)警，選取試驗(yàn)臺(tái)內(nèi)外主軸轉(zhuǎn)速300/600 r/min無(wú)故障運(yùn)行狀態(tài)下振動(dòng)加速度數(shù)據(jù)有效值平均值的1.25倍作為報(bào)警閾值[21]，即報(bào)警閾值為2.4654 m/s2。報(bào)警效果如圖4所示，在低工況(工況1、工況2)下，內(nèi)圈故障、外圈故障和滾動(dòng)體故障的振動(dòng)有效值均低于閾值，會(huì)出現(xiàn)漏警；類(lèi)似的，在高工況(工況3～工況5)下則有嚴(yán)重的虛警現(xiàn)象。顯然，固定報(bào)警閾值難以滿(mǎn)足設(shè)備時(shí)變工況下故障準(zhǔn)確預(yù)警的需求。

此外，進(jìn)一步對(duì)各工況分別設(shè)定振動(dòng)加速度有效值報(bào)警閾值，即以各工況正常數(shù)據(jù)振動(dòng)加速度有效值平均值的1.25倍作為該工況振動(dòng)閾值。報(bào)警效果如圖4所示，仍存在故障值小于閾值的漏警現(xiàn)象，說(shuō)明振動(dòng)加速度有效值不能全面反映出設(shè)備運(yùn)行狀態(tài)，僅以振動(dòng)加速度有效值單一特征設(shè)定報(bào)警閾值亦無(wú)法滿(mǎn)足設(shè)備預(yù)警需求。

圖4 振動(dòng)有效值固定閾值預(yù)警效果圖

3.3 所提方法預(yù)警分析

3.3.1 工況識(shí)別

從已人為標(biāo)定的工況1～工況5的轉(zhuǎn)速信號(hào)中各選取2000組數(shù)據(jù)，并以1∶1的比例進(jìn)行劃分，分別作為BPNN的訓(xùn)練數(shù)據(jù)和測(cè)試數(shù)據(jù)。

BPNN參數(shù)設(shè)置如表3所示，為2-10-1的3層簡(jiǎn)單網(wǎng)絡(luò)。利用同工況的訓(xùn)練集和測(cè)試集訓(xùn)練并測(cè)試BPNN模型，工況識(shí)別效果如圖5所示，BPNN能準(zhǔn)確自動(dòng)識(shí)別設(shè)備運(yùn)行工況。由于實(shí)驗(yàn)限制，本實(shí)驗(yàn)只考慮內(nèi)外轉(zhuǎn)軸轉(zhuǎn)速，未考慮負(fù)載變化，只模擬5個(gè)運(yùn)行工況，但基于BPNN的工況識(shí)別方法可通過(guò)調(diào)節(jié)網(wǎng)絡(luò)節(jié)點(diǎn)參數(shù)在多工況參數(shù)、多工況區(qū)間下自動(dòng)準(zhǔn)確識(shí)別設(shè)備運(yùn)行工況，故推廣性強(qiáng)。

表3 BPNN參數(shù)

圖5 工況識(shí)別效果圖

3.3.2 特征提取

選取5個(gè)工況下各120組正常狀態(tài)的振動(dòng)數(shù)據(jù)設(shè)定振動(dòng)閾值，根據(jù)表1中的計(jì)算公式提取特征。其中，本實(shí)驗(yàn)信號(hào)最大分析頻率為10240 Hz，將頻譜以2000 Hz為間隔劃分為0～2000 Hz、2001～4000 Hz、4001～6000 Hz、6001～8000 Hz、80001～10000 Hz共5個(gè)頻帶，可求得5個(gè)頻帶能量占比特征，共為25維特征。為了降低數(shù)據(jù)冗余度，利用PCA對(duì)提取的25維特征進(jìn)行降維。按照?qǐng)D5中識(shí)別出的工況區(qū)間，對(duì)各工況特征降維，得到各工況下低維特征矩陣，工況2的低維特征如圖6所示(篇幅有限，其他工況低維特征不予展示)，在軸承各狀態(tài)下，振動(dòng)特征間差異明顯，即提取的特征對(duì)軸承狀態(tài)敏感，可有效反映出雙轉(zhuǎn)子試驗(yàn)臺(tái)運(yùn)行狀態(tài)。

圖6 工況2降維特征效果圖

3.3.3 報(bào)警閾值設(shè)定

利用上述提取的各工況低維特征，基于SVM-PDE分別計(jì)算各工況下正常數(shù)據(jù)的概率密度值，并設(shè)定各工況下的振動(dòng)概率閾值。首先基于正常數(shù)據(jù)求解式(12)的二次規(guī)劃問(wèn)題，得到正常數(shù)據(jù)的系數(shù)向量β；然后結(jié)合式(7)和式(8)計(jì)算正常數(shù)據(jù)概率密度值，最后以正常數(shù)據(jù)概率密度值邊界即下限作為振動(dòng)報(bào)警閾值?？紤]到數(shù)據(jù)的波動(dòng)特性，將報(bào)警限下調(diào)3倍標(biāo)準(zhǔn)差，求得的各工況振動(dòng)閾值如表4所示。

表4 各工況振動(dòng)閾值表

3.3.4 報(bào)警狀態(tài)判別

基于測(cè)試數(shù)據(jù)，結(jié)合正常數(shù)據(jù)系數(shù)向量β及式(7)和式(8)計(jì)算測(cè)試數(shù)據(jù)的概率密度值，基于表4中的振動(dòng)概率閾值和測(cè)試數(shù)據(jù)的概率密度值判別測(cè)試數(shù)據(jù)狀態(tài)。其中測(cè)試數(shù)據(jù)組成如表5所示，包括各工況下的60組正常數(shù)據(jù)、180組內(nèi)圈故障數(shù)據(jù)、180組外圈故障數(shù)據(jù)和180組滾動(dòng)體故障數(shù)據(jù)。

表5 測(cè)試數(shù)據(jù)列表

各工況下的振動(dòng)概率閾值與測(cè)試數(shù)據(jù)的概率密度值間的大小關(guān)系如圖7所示，正常數(shù)據(jù)概率值均大于閾值，故障數(shù)據(jù)概率值均小于閾值，按照?qǐng)?bào)警狀態(tài)的判別規(guī)則，待判別數(shù)據(jù)的概率密度值小于振動(dòng)概率閾值會(huì)被判定為報(bào)警狀態(tài)，反之則為正常狀態(tài)。從圖7可以看出，該方法在各工況下虛警率、漏警率均為0，有效解決了固定閾值造成的虛警、漏警問(wèn)題，體現(xiàn)出所提出的方法的有效性。

圖7 基于BPNN和SVM-PDE的變工況振動(dòng)閾值預(yù)警效果圖

此外，當(dāng)設(shè)定振動(dòng)概率閾值的正常數(shù)據(jù)不足時(shí)，因數(shù)據(jù)無(wú)法涵蓋設(shè)備正常運(yùn)行時(shí)全部的振動(dòng)狀態(tài)，可能會(huì)產(chǎn)生部分虛警。因此在實(shí)際應(yīng)用中，應(yīng)隨設(shè)備正常數(shù)據(jù)庫(kù)的不斷完善，調(diào)整、更新振動(dòng)概閾值以降低虛警率，實(shí)現(xiàn)設(shè)備的準(zhǔn)確預(yù)警。

4 結(jié)束語(yǔ)

針對(duì)固定閾值預(yù)警方法易造成振動(dòng)監(jiān)測(cè)嚴(yán)重虛警、漏警問(wèn)題，開(kāi)展基于BPNN和SVM-PDE的旋轉(zhuǎn)機(jī)械變工況預(yù)警方法研究，將傳統(tǒng)SVM核函數(shù)改造為概率密度函數(shù)，提出SVM-PDE方法設(shè)定閾值進(jìn)行故障預(yù)警，該方法具有以下優(yōu)點(diǎn)。

① 基于BPNN的設(shè)備運(yùn)行工況識(shí)別方法可在多工況參數(shù)、多工況區(qū)間下自動(dòng)準(zhǔn)確識(shí)別設(shè)備所處運(yùn)行工況，為設(shè)備變工況預(yù)警奠定基礎(chǔ)。

② 結(jié)合信號(hào)處理方法從時(shí)域、頻域充分挖掘信號(hào)中有用信息，并利用SVM核函數(shù)融合多維特征信息，相較于單一特征更能反映設(shè)備運(yùn)行狀況。

③ 基于SVM-PDE設(shè)定閾值的故障預(yù)警方法相較于依據(jù)經(jīng)驗(yàn)對(duì)單一特征設(shè)定固定閾值的預(yù)警方法能在有效控制虛警率前提下，大幅降低故障漏警率；實(shí)驗(yàn)數(shù)據(jù)下實(shí)現(xiàn)了零虛警、零漏警，對(duì)于變工況旋轉(zhuǎn)機(jī)械振動(dòng)監(jiān)測(cè)具有較好的工程應(yīng)用價(jià)值。