孫留存,胡從川,錢大龍
(中國綠發(fā)投資集團有限公司,北京 100020)
隨著現(xiàn)代化科技與工業(yè)的迅猛發(fā)展,旋轉(zhuǎn)機械設(shè)備成為應(yīng)用最為廣泛的機械設(shè)備之一[1],其工況不但會影響自身機械設(shè)備的運行狀態(tài),還會直接影響后續(xù)生產(chǎn)環(huán)節(jié),因此對旋轉(zhuǎn)機械設(shè)備工況實時監(jiān)測至關(guān)重要。
目前,旋轉(zhuǎn)機械設(shè)備監(jiān)測技術(shù)已受到國內(nèi)外學者的高度重視,國外故障監(jiān)測技術(shù)發(fā)展較早,1965年,快速傅里葉變換的提出與應(yīng)用為故障監(jiān)測奠定技術(shù)基礎(chǔ)[2],而后經(jīng)歷10余年的發(fā)展,頻譜分析和電子測量技術(shù)被引入機械故障監(jiān)測之中,在1985年左右,微型計算機成為機械故障監(jiān)測的重心,商業(yè)化機械故障監(jiān)測系統(tǒng)隨之誕生;國內(nèi)故障監(jiān)測技術(shù)起步較晚,于20世紀70年代末引入國外先進技術(shù),經(jīng)本土化創(chuàng)新后逐漸發(fā)展起來,目前已形成帶有我國自身特色的機械故障監(jiān)測理論與技術(shù)并開發(fā)出眾多機械故障監(jiān)測系統(tǒng)投入使用。陳仁祥等[3]采用連續(xù)小波變換處理信號數(shù)據(jù),獲取時間尺度矩陣并構(gòu)造環(huán)境狀態(tài)空間,利用卷積神經(jīng)網(wǎng)絡(luò)擬合強化學習中Q函數(shù)生成深度Q網(wǎng)絡(luò),將環(huán)境狀態(tài)空間輸入深度Q網(wǎng)絡(luò),輸出狀態(tài)特征表示,通過智能體和環(huán)境狀態(tài)空間交互學習最大化Q函數(shù)值,實現(xiàn)旋轉(zhuǎn)機械設(shè)備故障監(jiān)測;劉頡等[4]采用歸一化后故障信號生成Hankel矩陣,通過奇異值分解獲取信號特征向量構(gòu)建圖數(shù)據(jù),引入GCN提取圖數(shù)據(jù)故障特征,利用softmax分類器實現(xiàn)故障監(jiān)測;侯召國等[5]融合深度殘差網(wǎng)絡(luò)和長短時記憶網(wǎng)絡(luò)提升故障監(jiān)測模型時序信息捕獲能力,在殘差塊中加入Dropout層,增強故障監(jiān)測能力,將預(yù)處理后旋轉(zhuǎn)機械設(shè)備信號輸入網(wǎng)絡(luò)模型,實現(xiàn)旋轉(zhuǎn)機械設(shè)備故障診斷與監(jiān)測。但以上方法沒有降維故障信號,導(dǎo)致訓練時間長、故障監(jiān)測準確率低。
為解決上述方法中存在的問題,引入集合經(jīng)驗?zāi)B(tài),提出一種基于WSN的旋轉(zhuǎn)機械設(shè)備故障時頻監(jiān)測方法。
在傳統(tǒng)旋轉(zhuǎn)機械設(shè)備故障監(jiān)測中通常采用有線方式采集信號,適應(yīng)性和靈活性較差,無法滿足現(xiàn)代智能化工業(yè)的需要,因此,本文方法引入無線傳感器網(wǎng)絡(luò)[6](WSN)用于旋轉(zhuǎn)機械設(shè)備故障信號的采集。
為避免故障監(jiān)測過程中出現(xiàn)單一節(jié)點故障或誤差較大等問題,WSN通常采用多節(jié)點同時監(jiān)測同一設(shè)備關(guān)鍵點的方式,但由此會導(dǎo)致采集的信號數(shù)據(jù)相似度過高的情況,而上位機用戶往往僅需獲取旋轉(zhuǎn)機械設(shè)備工況,并非全部節(jié)點詳細數(shù)據(jù)細節(jié),若傳輸全部無線傳感器節(jié)點數(shù)據(jù)至Sink節(jié)點,不但對故障監(jiān)測沒有實際幫助,還會造成節(jié)點能量的浪費,因此,需要在各節(jié)點向Sink節(jié)點傳輸數(shù)據(jù)前在本地對數(shù)據(jù)預(yù)處理,降低通信代價,提高故障監(jiān)測效率。
引入LEACH路由協(xié)議分簇節(jié)點,構(gòu)建網(wǎng)絡(luò)融合模型,如圖1所示,具體步驟如下所述。
圖1 WSN數(shù)據(jù)融合模型
a.數(shù)據(jù)級融合:各無線傳感器節(jié)點采集故障信號數(shù)據(jù)并處理信號,濾除信號噪聲,提取時頻信息。
b.特征級融合:各成員節(jié)點處理后數(shù)據(jù)傳輸至簇頭節(jié)點,簇頭節(jié)點采用主成分分析法降維數(shù)據(jù)。
c.決策級融合:簇頭節(jié)點將降維后特征數(shù)據(jù)傳輸至Sink節(jié)點,Sink節(jié)點通過免疫聚類RBF神經(jīng)網(wǎng)絡(luò)模型監(jiān)測旋轉(zhuǎn)機械設(shè)備運行狀況。
通過以上流程即可實現(xiàn)信號數(shù)據(jù)的分析和整理,最終經(jīng)由上位機輸出旋轉(zhuǎn)機械設(shè)備運行狀況信息,實現(xiàn)旋轉(zhuǎn)機械設(shè)備監(jiān)測。
目前,旋轉(zhuǎn)機械設(shè)備故障監(jiān)測最常用的方法是分析信號振動信息進而判定故障原因,但由于旋轉(zhuǎn)機械設(shè)備結(jié)構(gòu)和振源較為復(fù)雜,導(dǎo)致采集到的信號中噪聲較多,若想要準確判定旋轉(zhuǎn)機械設(shè)備運行狀況,必須對采集到的信號去噪處理。本文方法引入小波包閾值去噪法處理故障信號,獲取質(zhì)量更高的振動信號用于后續(xù)計算。
小波包閾值去噪法根據(jù)信號和噪聲在不同尺度上小波包系數(shù)特性的區(qū)別對兩者加以判斷,保留原始信號小波包系數(shù),濾除噪聲分量中的小波分量,通過小波包重構(gòu)去噪信號,具體流程如下所述。
(1)
b.依據(jù)事先設(shè)定的熵標準確定最優(yōu)小波包基。
c.確定閾值θ,采用θ處理小波包分解系數(shù),用L表示信號長度,σ表示噪聲標準差,則閾值θ的確定方式為
(2)
傳統(tǒng)小波包閾值法在去噪過程中往往采用統(tǒng)一閾值去噪信號,但該方法會引發(fā)過扼殺的情況,因此,本文方法通過自適應(yīng)閾值選取法確定小波包分解中各層小波包系數(shù)閾值[7],求解j個閾值均值并以該值為最終去噪的全局閾值?。
(3)
由于旋轉(zhuǎn)機械設(shè)備結(jié)構(gòu)復(fù)雜且各部件間傳動具有非線性特征,導(dǎo)致振動信號存在不平穩(wěn)等問題,為故障監(jiān)測帶來一定難度,因此 ,本文引入集合經(jīng)驗?zāi)B(tài)分解法提取旋轉(zhuǎn)機械設(shè)備信號時頻特征,通過時頻特征描述故障信息。采用EEMD提取時頻特征主要流程如下所述。
a.分解時刻t不同狀態(tài)的振動信號,生成主IMF分量ci(t)。
(4)
依據(jù)式(4)可獲取歸一化能量指標,記作tp1=z(1),tp2=z(2),…,tpn=z(N)。
(5)
根據(jù)式(5)可得到IMF矩陣中的奇異熵指標,記作H1,H2,…,Hn。
結(jié)合歸一化能量指標和IMF矩陣奇異譜熵指標,完成時頻特征的提取。
主成分分析法是應(yīng)用最為廣泛的數(shù)據(jù)降維方法之一,將其用于故障監(jiān)測中能夠選取盡可能少的指標表示故障信號特征,避免冗余特征對故障監(jiān)測造成不必要的干擾。
(6)
采用v1,v2,…,vm構(gòu)造特征矩陣V=[v1,v2,…,vm],S=diag(a1,a2,…,am)表示利用a1,a2,…,am建立的對角矩陣,分解協(xié)方差矩陣C,得到C=V·S·VT。
通過時頻特征向量構(gòu)造主成分分析子空間,采用該子空間即可提供P中絕大部分特征信息,實現(xiàn)時頻特征降維。
傳統(tǒng)旋轉(zhuǎn)機械設(shè)備故障監(jiān)測通常根據(jù)振動頻譜分析結(jié)果人工判定故障情況,但由于技術(shù)人員理論和經(jīng)驗的限制,無法避免出現(xiàn)誤判、漏判的情況,為此,本文方法構(gòu)建免疫RBF神經(jīng)網(wǎng)絡(luò)監(jiān)測旋轉(zhuǎn)機械設(shè)備故障,由于RBF神經(jīng)網(wǎng)絡(luò)在訓練中隱含層數(shù)據(jù)中心總數(shù)和位置獲取相對困難,因此引入免疫聚類算法獲取RBF神經(jīng)網(wǎng)絡(luò)相關(guān)參數(shù),提高旋轉(zhuǎn)機械設(shè)備故障監(jiān)測準確率[8-9]。RBF主要結(jié)構(gòu)如圖2所示。
圖2 RBF神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)
圖3 實驗現(xiàn)場
用X={x1,x2,…,xm}表示輸入的降維后故障特征,?xk∈X;Y={y1,y2,…,yn}表示輸出的故障監(jiān)測結(jié)果,?yj∈Y;Z={z1,z2,…,zN}表示隱含層輸出值,?zi∈Z;B={b1,b2,…,bm}表示輸出層偏置,?bj∈B,則zi和yj計算方式為
(7)
構(gòu)建RBF神經(jīng)網(wǎng)絡(luò)模型[10-11]的關(guān)鍵主要有以下3點:
a.確定隱含層節(jié)點中心總數(shù)。
b.確定節(jié)點中心位置。
c.輸出權(quán)值ωij。
因此,本文方法引入免疫聚類算法自適應(yīng)選取隱含層中心總數(shù)和位置,并通過歐氏距離描述向量之間匹配度。免疫聚類算法由克隆選擇和免疫抑制2部分構(gòu)成,主要流程如下所述。
a.隨機生成初始化免疫聚類算法聚類中心C0,設(shè)置最大迭代次數(shù)K。
b.計算xk與C0中全部中心匹配度rij=
c.獲取o個與xk匹配度最高的中心向量,記作Co=[c1,c2,…,co]。
d.克隆Co中全部向量,匹配度越高,則克隆數(shù)量Mc越大,用T表示克隆規(guī)模,則Mc∝T×rij。
e.對克隆后中心向量執(zhí)行變異操作,用η表示學習速率,η與rij呈正相關(guān),得到變異后中心向量集Cp,ck=ck-η(ck-xk)。
h.重復(fù)以上步驟直到迭代次數(shù)達到K。
確定隱含層數(shù)據(jù)中心后對RBF神經(jīng)網(wǎng)絡(luò)加以訓練,在訓練后網(wǎng)絡(luò)中輸入預(yù)處理的旋轉(zhuǎn)機械設(shè)備故障信號數(shù)據(jù),實現(xiàn)故障監(jiān)測。
為驗證本文的旋轉(zhuǎn)機械設(shè)備故障時頻監(jiān)測方法整體有效性,需要進行實驗分析。采用CWRU電氣工程實驗中心的旋轉(zhuǎn)機械設(shè)備振動數(shù)據(jù)作為實驗數(shù)據(jù)集,以MOS管折彎機整型機為實驗對象,對其進行故障監(jiān)測研究,實驗現(xiàn)場傳感器設(shè)置如3圖所示。
分別采用無負載和0.746 kW負載下的8種不同旋轉(zhuǎn)機械設(shè)備狀態(tài)數(shù)據(jù)集作為模型檢測樣本,經(jīng)處理后每種狀態(tài)中包含樣本600個,實驗樣本具體信息如表1所示。
表1 旋轉(zhuǎn)機械設(shè)備狀態(tài)數(shù)據(jù)集信息
為檢驗本文方法、文獻[3]方法和文獻[4]方法的旋轉(zhuǎn)機械設(shè)備故障監(jiān)測能力,測試不同方法的監(jiān)測耗時,結(jié)果如圖4所示。
圖4 監(jiān)測耗時測試結(jié)果
由圖4可以看出,與文獻[3]方法和文獻[4]方法相比,本文方法的監(jiān)測耗時更短,即采用本文方法能夠有效提高信號處理的速度,說明本文方法監(jiān)測模型具有更強的魯棒性。這是因為本文方法引入集合經(jīng)驗?zāi)B(tài)分解法提取旋轉(zhuǎn)機械設(shè)備信號時頻特征,該方法能夠有效抑制振動信號不平穩(wěn)問題,降低故障監(jiān)測復(fù)雜度,提升監(jiān)測效率。
進一步檢測本文方法、文獻[3]方法和文獻[4]方法在相同負載和不同負載下的旋轉(zhuǎn)機械設(shè)備故障監(jiān)測準確率,結(jié)果如圖5所示。
圖5 故障監(jiān)測準確率檢測結(jié)果
由圖5可以看出,本文方法的故障監(jiān)測準確率明顯高于文獻[3]方法和文獻[4]方法,監(jiān)測準確率在95%左右,而文獻[3]方法和文獻[4]方法的監(jiān)測準確率平均值均未達到95%,說明本文方法在旋轉(zhuǎn)機械設(shè)備故障監(jiān)測中具有更強的自適應(yīng)性和魯棒性,因為本文方法引入主成分分析法降維時頻特征,避免冗余特征對故障監(jiān)測造成不必要的干擾,從而提高了本文方法的故障監(jiān)測能力。
近年來,隨著工業(yè)的自動化和智能化發(fā)展,連續(xù)化、大型化以及大載荷、大功率等逐漸成為旋轉(zhuǎn)機械設(shè)備發(fā)展的主流趨勢。旋轉(zhuǎn)機械設(shè)備在使用過程中始終受到多種外界作用,導(dǎo)致其運行狀態(tài)產(chǎn)生一定的改變,進而造成異常、故障等情況的發(fā)生,影響系統(tǒng)的正常運行,因此對旋轉(zhuǎn)機械設(shè)備工況實時監(jiān)測是保障其安全運行的必要條件。為解決目前存在的訓練時間長、故障監(jiān)測準確率低的問題,本文通過小波包閾值去噪法去噪無線傳感器采集的故障信號,利用EEMD提取信號時頻特征,引入主成分分析法降維時頻特征,采用免疫RBF神經(jīng)網(wǎng)絡(luò)完成旋轉(zhuǎn)機械設(shè)備故障監(jiān)測。實驗結(jié)果表明,本文方法能夠有效地縮短訓練時間、提高故障監(jiān)測準確率。