劉笑凡，余 刃，王 石，王天舒，謝旭陽

基于HHT的電動泵便攜式故障檢測裝置

劉笑凡，余 刃，王 石，王天舒，謝旭陽

（海軍工程大學(xué)核科學(xué)與技術(shù)學(xué)院，武漢 430033）

為了有效地檢測電動泵運行異常狀況，本文在振動信號時域、頻域分析的基礎(chǔ)上，引入時頻分析方法Hilbert-Huang變換（Hilbert-Huang Transform, HHT），適當選取對各種故障模式敏感的特征參數(shù)，構(gòu)成特征參數(shù)集。通過對電動泵振動信號各特征參數(shù)在不同故障模式下的表現(xiàn)進行分析，實現(xiàn)對電動泵運行狀態(tài)和故障模式的判斷。采用某數(shù)據(jù)中心網(wǎng)站的電動泵軸承故障數(shù)據(jù)作為樣本數(shù)據(jù)驗證了方法的有效性。結(jié)果表明，通過引入HHT后的特征參數(shù)集，可以更有效地對數(shù)據(jù)中的四種運行狀態(tài)進行判斷?；谒岢龅奶卣鲄?shù)集和分析方法，開發(fā)了電動泵便攜式故障檢測裝置。

振動信號 時頻分析 HHT 故障診斷 電動泵

0 引言

電動泵作為船舶動力裝置的重要輔助機械，對其安全運行起著重要作用。在船舶動力裝置中，由于運行環(huán)境和安裝空間的限制，對一些電動泵的監(jiān)測手段有限，難以及時檢測和發(fā)現(xiàn)其異常運行狀態(tài)，從而錯過最佳處理時間。

電動泵運轉(zhuǎn)時會產(chǎn)生振動，當出現(xiàn)異常的時候，其振動信號必然也會出現(xiàn)異常，因而通過監(jiān)測電動泵運行的振動信號判斷電動泵運行是否出現(xiàn)異常，是一種較為常用的監(jiān)測手段[1]。電動泵振動信號的特征提取和故障識別是故障診斷中的兩個重要部分。由于其振動信號具有非線性和非平穩(wěn)性的特點，傳統(tǒng)的以傅里葉變換為基礎(chǔ)的方法難以取得較好的分析效果[2]。鑒于傳統(tǒng)方法的局限性，以及船舶動力裝置的特殊性，本文提出基于HHT的核動力裝置電動泵故障診斷方法。希爾伯特-黃變換（Hilbert-Huang Transform, HHT），是由Norden E. Huang（黃鍔）等人與1998年提出的，并引入了Hilbert譜的概念和Hilbert譜分析的方法，在非線性、非平穩(wěn)信號的故障診斷中得到了廣泛引用[3,4]。

在進行電動泵運行狀態(tài)判斷時，需要適當選取若干能夠表征電動泵運行狀態(tài)的特征參數(shù)，構(gòu)成電動泵狀態(tài)檢測和故障診斷的特征參數(shù)集。選取的特征參數(shù)集應(yīng)對各種故障運行狀態(tài)具有足夠的敏感性，并應(yīng)保證通過分析其中各特征參數(shù)在不同運行狀態(tài)（含正常和故障狀態(tài)）下的表現(xiàn)，能夠有效區(qū)分各種運行狀態(tài)，從而達到異常運行狀態(tài)監(jiān)測和故障診斷的目的。本文針對船舶動力裝置電動泵缺乏方便有效的故障檢測手段問題，研究了電動泵振動信號的HHT分析方法，提出了綜合時域、頻域和時頻域分析的特征參數(shù)集，并利用美國凱斯西儲大學(xué)軸承數(shù)據(jù)中心網(wǎng)站的電動泵軸承故障數(shù)據(jù)作為樣本數(shù)據(jù)，進行運行狀態(tài)分析和故障判定的實驗驗證。開發(fā)了便攜式電動泵故障檢測裝置，實現(xiàn)對船舶動力裝置中電動泵運行狀態(tài)的監(jiān)測。

1 基于HHT的振動信號處理方法

不論是在完全時域或是完全頻域上處理非平穩(wěn)信號，都會造成信號局部信息的缺失。而對于電動泵運行狀態(tài)監(jiān)測，獲取其振動信號的局部特性又是非常重要的，因而引入信號的時頻分析方法對振動信號進行處理。HHT就是一種信號的時頻分析方法。

HHT主要由兩部分組成：第一部分為經(jīng)驗?zāi)B(tài)分解（Empirical Mode Decomposition, EMD）；第二部分為對任意信號經(jīng)過EMD后得到的各本征模函數(shù)（Intrinsic Mode Function, IMF）進行Hilbert變換，從而推導(dǎo)得出Hilbert時頻譜和Hilbert邊際譜。該方法與以傅里葉變換理論為基礎(chǔ)的時頻分析方法有較大不同，能很好地解釋以傅里葉變換為基礎(chǔ)的方法所不能解釋的現(xiàn)象。

1.1 經(jīng)驗?zāi)B(tài)分解

EMD方法的基本思想是：利用采集到的時域信號上下包絡(luò)線的平均值確定瞬時平衡位置，進而提取本征模函數(shù)（IMF）。IMF反映了信號中存在的不同特征的成分，由于電動泵發(fā)生局部損傷故障時振動信號表現(xiàn)出高頻調(diào)幅信號與低頻諧波相疊加的特征，可以推斷，EMD能夠分離出信號中反映損傷特征的高頻調(diào)幅信號，從而為故障診斷提供依據(jù)[3]。一個IMF必須滿足兩個條件：

1）函數(shù)在整個時間范圍內(nèi)，局部極值點和過零點的數(shù)目相等，或者最多相差一個；

2）任意時刻的局部最大值的包絡(luò)線（上包絡(luò)線）和局部最小值的包絡(luò)線（下包絡(luò)線）的平均值必須為0。

EMD方法必須基于以下假設(shè)條件：

1）數(shù)據(jù)至少有2兩個極值，一個極大值和一個極小值；

2）數(shù)據(jù)的局部時域特性是由極值點間的時間尺度唯一確定。

EMD的具體步驟如下[9]：

1）求出原始數(shù)據(jù)()的局部極大值和極小值，利用三次樣條插值連接，分別得到極大值包絡(luò)線x()和極小值包絡(luò)線x()；

2）對每個時刻局部極大值x()和極小值x()取平均，得到瞬時平均值()；

3）用原始數(shù)列()減去瞬時平均值得到一個去掉低頻的新數(shù)列()，即

4）檢查()是否滿足IMF的兩個條件，若滿足，則()為一個IMF；若不滿足，將()作為原始數(shù)列重復(fù)上述3個步驟，直到滿足兩個條件為止，這樣就得到了第一個IMF，記為c()。一般來說，c()代表了原始數(shù)列中的高頻部分，也稱c()為原始數(shù)列的一個IMF；

5）將原始數(shù)列減去c()得到一個去掉高頻組分的差值數(shù)列r()，對r()進行上述平穩(wěn)化處理過程可以得到第2個IMF分量c()，如此重復(fù)下去直到最后一個序列r()不可再分解（為單調(diào)函數(shù)或常數(shù)）為止，此時r()代表原始數(shù)列的均值或趨勢。

最后，原始的數(shù)據(jù)序列即由這些IMF分量以及一個趨勢或均值表示：

1.2 對本征模函數(shù)做Hilbert變換

將IMF做快速傅里葉變換（FFT）可以得到振動信號的Hilbert譜。

對IMF做Hilbert變換:

式中，

上式中定義的瞬時振幅和瞬時相位反映了信號的瞬時特性。在此基礎(chǔ)上定義的瞬時頻率為：

在概念上，瞬時頻率可以理解為一個標準正弦波在局部最佳逼近被分析信號的頻率值。在統(tǒng)計上它與傅里葉頻率是相容的，傅里葉頻率是全局量，而瞬時頻率是局部表述方式；傅里葉頻率是獨立的，而瞬時頻率是時間的函數(shù)。

由以上可以看出，對IMF所做的Hilbert變換得出的振幅和頻率都是時間的函數(shù)，那么將每個IMF的解析信號求和，即得到原始信號()的解析信號：

式中a表示瞬時頻率下的振幅，ω表示瞬時頻率。

稱上式為()的Hilbert譜，記為()。將其對時間進行積分得到邊際譜：

將其平方對積分即得到瞬時能量譜，對時間積分即得到能量譜：

2 特征參數(shù)集的選取

將振動信號轉(zhuǎn)換為頻域信號是判斷電動泵運行狀態(tài)的一種常用方法，這是因為電動泵高速的周期性轉(zhuǎn)動。當出現(xiàn)故障時，其振動信號的變化在頻域上的反映尤為明顯。除頻域分析以外，振動信號的時域分析也能夠一定程度地反映電動泵的運行狀態(tài)。為了增加判斷的準確性，本文綜合應(yīng)用了信號的時域、頻域和時頻域的分析結(jié)果，構(gòu)成判斷運行狀態(tài)的特征參數(shù)。針對時域和頻域信號，主要選取了以下特征參數(shù)：

1）振動烈度

振動烈度通常用于表征振動水平，一般使用機械設(shè)備測點振動信號的均方根值表示。其計算公式如下：

其中表示所測信號的長度，()表示振動速度。

2）峭度

峭度是無量綱參數(shù)，是描述某變量所有取值分布形態(tài)陡峭程度的統(tǒng)計量。它與電動泵泵體的參數(shù)（如轉(zhuǎn)速、尺寸、載荷等）無關(guān)，但是對沖擊信號非常敏感。因而這種特點特別適用于表面損傷類故障、尤其是早期故障的監(jiān)測。

峭度的計算方法如下：

在實際應(yīng)用中，通常將峭度值做減3處理，使得正態(tài)分布的峭度為0。

3）頻譜

頻譜是判斷振動信號是否運行穩(wěn)定的重要參數(shù)，一旦電動泵出現(xiàn)非常明顯的振動異常，可以通過頻譜信號直接進行判斷，同時也可通過頻譜信號中出現(xiàn)的倍頻信號對運行狀態(tài)進行判斷。頻譜信號是通過對采集到的振動信號做快速傅里葉變換（FFT）得到的，其基本運算方法如下：

上式中，()為時域信號，F()為時域信號經(jīng)過變換后得到的頻域信號。

4）功率譜

最終，本文選取了振動信號的IMF分量、Hilbert邊際譜和能量譜、振動烈度、峭度、頻譜、功率譜等參數(shù)構(gòu)成電動泵狀態(tài)監(jiān)測的特征參數(shù)集。此特征參數(shù)集中的各參數(shù)在電動泵的不同運行狀態(tài)下會有不同的表現(xiàn)形式。將其作為判斷電動泵運行的依據(jù)，與已知的正?；蚬收线\行狀態(tài)所對應(yīng)的特征參數(shù)集進行對比，即可判斷電動泵當前的運行狀態(tài)。

3 特征參數(shù)集對不同運行狀態(tài)的敏感性分析

3.1 特征參數(shù)的計算

為驗證所選取的特征參數(shù)集對電動泵各類故障的敏感性，以及用其進行故障判定的有效性，本文采用來自美國凱斯西儲大學(xué)軸承數(shù)據(jù)中心網(wǎng)站的軸承正常和故障數(shù)據(jù)[10]進行分析。該數(shù)據(jù)為1.47 kW 電機，轉(zhuǎn)速在1730-1797 rpm 范圍波動下運行時的振動數(shù)據(jù)。其中，所選取的驅(qū)動端軸承振動數(shù)據(jù)的采樣頻率是48000 Hz，數(shù)據(jù)包括正常運行、內(nèi)圈故障、滾珠故障、外圈故障等四種狀態(tài)滾動軸承運行的振動信號，每組信號約記錄有48萬個數(shù)據(jù)。表1給出了振動信號的時域和頻域特征參數(shù)計算結(jié)果，表2給出了振動信號的時頻域特征信號計算結(jié)果，其中的Hilbert譜和對IMF做Hilbert變換的結(jié)果，是對經(jīng)驗?zāi)B(tài)分解中得出的能夠有效表征電動泵運行狀態(tài)的IMF做譜分析和Hilbert變換得到的。

表1 振動信號的時域、頻域分析結(jié)果

表2 振動信號的時頻分析結(jié)果

3.2 利用特征參數(shù)集對電動泵運行狀態(tài)的分析

分析表1中的振動信號時域、頻域參數(shù)計算結(jié)果，可以發(fā)現(xiàn)，正常狀態(tài)下的振動烈度和功率譜中的最大功率較其他運行狀態(tài)要小一個數(shù)量級左右。因而，通過振動烈度和功率譜分析可以甄別電動泵的正常運行狀態(tài)和故障運行狀態(tài)。但是僅僅通過這兩種特征參數(shù)不能有效區(qū)分三種故障狀態(tài)，所以需要綜合分析特征參數(shù)集的各特征參數(shù)，才能夠判斷電動泵所處的運行狀態(tài)。

分析表1中的各運行狀態(tài)的峭度值，可以將滾珠故障和外圈故障與其他兩種運行狀態(tài)區(qū)別開，但即使結(jié)合振動烈度、頻譜、功率譜等參數(shù)也不能有效區(qū)分滾珠故障和外圈故障兩種狀態(tài)。分析表2 中IMF做Hilbert變換的結(jié)果，可以看出，引入HHT后，變換后的波形可以明顯地區(qū)分滾珠故障和外圈故障。而此時通過分析表2 中IMF的Hilbert頻譜，可看出內(nèi)圈故障狀態(tài)下頻域內(nèi)的最大振幅比正常狀態(tài)下的最大振幅高出約一個數(shù)量級，因而可以通過IMF的Hilbert頻譜區(qū)分正常運行狀態(tài)與內(nèi)圈故障狀態(tài)。而經(jīng)過EMD后得到的IMF作為一種重要的中間量，也可以反映電動泵的運行狀態(tài)。通過上述分析可以得出，當時域、頻域分析無法明確區(qū)分電動泵的運行狀態(tài)時，通過引入HHT可以更加精確地對電動泵的各種運行狀態(tài)進行甄別。

綜上，振動烈度、峭度、頻譜、功率譜、經(jīng)EMD得到的IMF、IMF的Hilbert頻譜以及IMF的Hilbert變換等是對上述四種狀態(tài)敏感的特征參數(shù)，可以有效地對電動泵不同的運行狀態(tài)進行甄別。

4 便攜式故障診斷裝置的設(shè)計

在船舶航行中，為實現(xiàn)便捷地對動力系統(tǒng)中主要電動泵的運行狀態(tài)進行監(jiān)測，本文設(shè)計開發(fā)了便攜式電動泵運行狀態(tài)監(jiān)測裝置。該裝置的硬件結(jié)構(gòu)包括傳感器、數(shù)據(jù)采集裝置和平板計算機三個部分，其基本結(jié)構(gòu)如圖1所示，設(shè)備如圖2所示。在監(jiān)測過程中，將傳感器用磁鐵吸附在泵體上，采集電動泵的振動信號，故本裝置的應(yīng)用無需對原設(shè)備作任何改動。本裝置采用加速度傳感器拾取振動信號，通過NI-9234振動信號采集模塊采集振動信號數(shù)據(jù)，并采用NI-9181數(shù)據(jù)采集卡，通過以太網(wǎng)將振動信號數(shù)據(jù)傳輸給平板計`算機。采用LABVIEW開發(fā)了振動信號的采集與分析軟件，實現(xiàn)對振動信號采集的啟?？刂?，并對采集到的振動信號進行時域、頻域和時頻特征參數(shù)的分析計算，從而得到前文所述的特征參數(shù)集，提供給使用人員。使用人員根據(jù)計算結(jié)果，結(jié)合其經(jīng)驗，實現(xiàn)對電動泵運行狀態(tài)的分析判斷。

圖1 硬件總體結(jié)構(gòu)圖

5 結(jié)論

通過電動泵運行時產(chǎn)生的振動信號來判斷電動泵的運行狀態(tài)，是電動泵運行狀態(tài)監(jiān)測的一種常用的手段。在狀態(tài)監(jiān)測中，僅通過振動信號的時域和頻域分析法并不能對所有的運行狀態(tài)進行甄別。為此，本文引入了振動信號的時頻分析法HHT，結(jié)合時域和頻域分析法，構(gòu)成了能夠有效判別電動泵各種正常和故障運行狀態(tài)的特征參數(shù)集。采用來自美國凱斯西儲大學(xué)軸承數(shù)據(jù)中心網(wǎng)站的電動泵軸承故障數(shù)據(jù)作為樣本數(shù)據(jù)，對利用所選特征參數(shù)集進行各種運行狀態(tài)甄別的有效性進行了驗證。在此基礎(chǔ)上設(shè)計了便攜式故障診斷裝置的設(shè)計，通過振動信號的時域、頻域和時頻分析三種手段，實現(xiàn)了對船舶動力裝置電動泵運行狀態(tài)的監(jiān)測。

圖2 加固型便攜式計算機

[1] 徐卓飛, 劉凱, 張海燕等. 基于經(jīng)驗?zāi)Ｊ椒纸夂椭髟治龅臐L動軸承故障診斷方法研究[J]. 振動與沖擊, 2014, 33(23): 133-139.

[2] 許凡, 方彥軍, 張榮. 基于EEMD模糊熵的PCA-GG滾動軸承聚類故障診斷[J]. 計算機集成制造系統(tǒng), 2016, 22(11): 2631-2642.

[3] 高強, 杜小山, 范虹等. 滾動軸承故障的EMD診斷方法研究[J]. 振動工程學(xué)報, 2007, 20(1): 15-18.

[4] 熊衛(wèi)華, 趙光宙. 基于希爾伯特-黃變換的變壓器鐵心振動特性分析[J]. 電工技術(shù)學(xué)報, 2006, 21(8):9-13.

[5] 劉世金. Hilbert-Huang變換及其應(yīng)用研究[J]. 高師理科學(xué)刊, 2012, 32(4): 40-42.

[6] 鄒巖崑. 局域波分析的理論方法研究及應(yīng)用[D]. 大連理工大學(xué), 2004.

[7] 李博. 基于LabVIEW的HHT方法實現(xiàn)及在軸承故障診斷中的應(yīng)用研究[D]. 太原理工大學(xué), 2008.

[8] 高東璇. 基于BP神經(jīng)網(wǎng)絡(luò)的齒輪故障模式識別研究[D]. 長安大學(xué), 2010.

[9] 陳一賢. HHT方法分析[D]. 浙江大學(xué), 2007.

[10] csegroups.case.edu /bearingdatacenter /pages /48k-drive -end-bearing-fault-data

Portable Fault Detection Device for Electric Pumps Based on HHT

Liu Xiaofan, Yu Ren, Wang Shi, Wang Tianshu, Xie Xuyang

(Collage of Nuclear Energy Science and Engineering, Naval University of Engineering, Wuhan 430033, China)

TP206.3

A

1003-4862（2018）12-0012-06

2018-07-13

劉笑凡（1989-），男，研究生。研究方向：核動力控制與運行。E-mail: 1508969066@qq.com