任學(xué)平，王金磊，馬文生

（1.內(nèi)蒙古科技大學(xué) 機(jī)械工程學(xué)院，包頭 014010；2.哈爾濱工業(yè)大學(xué) 機(jī)電工程學(xué)院，哈爾濱 150001）

設(shè)備故障診斷的一個(gè)重要內(nèi)容就是對(duì)故障進(jìn)行定位，判定出現(xiàn)故障的部件。然而設(shè)備的各個(gè)部件之間是相互聯(lián)系的，因此設(shè)備具有故障傳播特性，即當(dāng)傳感器測(cè)得某一部件的振動(dòng)信號(hào)出現(xiàn)異常時(shí)，并不能馬上確定其為故障源，因?yàn)榇斯收犀F(xiàn)象有可能是其他部件的故障傳播過(guò)來(lái)導(dǎo)致的，這就需要對(duì)故障進(jìn)行尋因定位，避免發(fā)生誤判，因此故障定位的方法就顯得尤為重要。文獻(xiàn)[1]提出的利用相關(guān)分析進(jìn)行故障定位的方法只適用于傳播路徑單一的簡(jiǎn)單結(jié)構(gòu)系統(tǒng)，并且需要知道故障源的大概位置，而且沒(méi)有提及到信號(hào)的預(yù)處理方法；文獻(xiàn)[2]運(yùn)用的時(shí)域平均技術(shù)可以提高信噪比作為預(yù)處理方法，但需要在測(cè)量振動(dòng)加速度的同時(shí)測(cè)量一個(gè)轉(zhuǎn)速同步脈沖信號(hào)以確定定整周期；文獻(xiàn)[3]提出的應(yīng)用希爾伯特變換提取相位信號(hào)的方法沒(méi)有形成一個(gè)系統(tǒng)的方法，即沒(méi)有前處理和后續(xù)處理。本文根據(jù)故障的傳播特性，提出了一種利用相位特征傳播特性來(lái)進(jìn)行故障定位的方法。本方法包括信號(hào)的預(yù)處理，相位的提取和后處理和分析實(shí)現(xiàn)故障定位，形成了一套系統(tǒng)的方法。最后對(duì)軋機(jī)主傳動(dòng)系統(tǒng)故障診斷的實(shí)例分析驗(yàn)證了本方法的有效性。

1 基于相位傳播特性的分析方法

1.1 自適應(yīng)消噪技術(shù)

自適應(yīng)消噪技術(shù)(ANC)[4]是一種以維納濾波為基礎(chǔ)對(duì)混有噪聲的信號(hào)進(jìn)行提取處理的方法。其原理是最小平方濾波器加上反饋的相關(guān)回路,構(gòu)成自適應(yīng)消噪器。

LMS（Least Mean Square）算法是自適應(yīng)橫向?yàn)V波器中求最佳權(quán)系數(shù)的簡(jiǎn)單有效的遞歸算法,它不需要求相關(guān)矩陣，也不涉及矩陣求逆，而是運(yùn)用了最優(yōu)化的數(shù)學(xué)算法最陡下降法。

1.2 希爾伯特變換提取相位特征信息

對(duì)一調(diào)制信號(hào)進(jìn)行希爾伯特變換[5]，然后用原信號(hào)和希爾伯特變換得到的信號(hào)構(gòu)成一個(gè)解析信號(hào)，再對(duì)解析信號(hào)移頻、反正切變換就得到信號(hào)的相位估計(jì)。

通過(guò)上述方法計(jì)算得到的相位是信號(hào)的折疊相位，即在求值的時(shí)候超過(guò)2π的相位是進(jìn)行了模的2π運(yùn)算的，這樣在經(jīng)過(guò)反余弦求出來(lái)的相位在(0，2π)之間，與真實(shí)相位不等，即產(chǎn)生了相位折疊。此時(shí)需要加上2π的整數(shù)倍才能得到無(wú)折疊的瞬時(shí)相位。

1.3 振動(dòng)信號(hào)的平滑處理方法

利用傳感器測(cè)得的振動(dòng)信號(hào)中往往含有噪聲信號(hào)。噪聲信號(hào)除了有周期性的噪聲信號(hào)外，還有不規(guī)則的隨機(jī)干擾信號(hào)；由于隨即干擾信號(hào)的頻帶較寬，有時(shí)高頻成分所占的比例還比較大，使得采集離散數(shù)據(jù)得到的振動(dòng)曲線上呈現(xiàn)有許多毛刺，很不平滑。為了消弱干擾信號(hào)的影響，提高振動(dòng)曲線光滑度，常常需要對(duì)采樣數(shù)據(jù)進(jìn)行平滑處理。本文主要介紹五點(diǎn)三次平滑法[6]。

五點(diǎn)三次平滑法是利用最小二乘法原理對(duì)數(shù)據(jù)進(jìn)行三次最小二乘多項(xiàng)式平滑的方法。

1.4 振動(dòng)信號(hào)的相位傳播特性

齒輪箱的振動(dòng)系統(tǒng)是一個(gè)相當(dāng)復(fù)雜的非線性系統(tǒng)[7]。在研究齒輪箱故障時(shí)，通常將齒輪傳動(dòng)副進(jìn)行簡(jiǎn)化。齒輪傳動(dòng)副作為一個(gè)振動(dòng)系統(tǒng)，其動(dòng)力學(xué)方程為：

式中x——沿嚙合線上齒輪相對(duì)位移；

C——齒輪嚙合阻尼；

k(t)——嚙合剛度；

M——當(dāng)量質(zhì)量；

F(t)——外界激勵(lì)。

由齒輪箱的振動(dòng)機(jī)理可知，故障具有傳播特性。根據(jù)故障的傳播特性，相位特征的傳播與故障的位置和測(cè)點(diǎn)位置之間的關(guān)系有關(guān)[8]。例如一振動(dòng)測(cè)試的測(cè)點(diǎn)布置如圖1所示。若部件(1)出現(xiàn)故障，部件(2)無(wú)故障，由于測(cè)點(diǎn)1位于故障位置的前面，則由測(cè)點(diǎn)1測(cè)得的振動(dòng)信號(hào)的相位特征為超前，如圖2a所示；反之，由于測(cè)點(diǎn)2和測(cè)點(diǎn)3位于故障位置的后面，則由測(cè)點(diǎn)2或測(cè)點(diǎn)3測(cè)得的振動(dòng)信號(hào)的相位特征為滯后，如圖2b所示。

圖1 測(cè)點(diǎn)布置示意圖Fig.1 Measurement point layout schematic diagram

圖2 相位躍變波形圖Fig.2 Phase jump waveform diagram

由此可見(jiàn)，只要通過(guò)傳感器測(cè)得設(shè)備的振動(dòng)信號(hào)，然后對(duì)信號(hào)處理分析得到相位波形，就能根據(jù)相位波形的躍變特征來(lái)進(jìn)行故障定位?；诖嗽恚覀兛梢酝卣沟烬X輪箱之外部件出現(xiàn)故障的定位。

1.5 基于相位特征傳播特性的故障定位方法

本文提出的基于相位特征傳播特性的故障定位方法其具體步驟為：1）先通過(guò)自適應(yīng)消噪技術(shù)除去振動(dòng)信號(hào)中的強(qiáng)相關(guān)噪聲；2）再應(yīng)用希爾伯特變換提取相位信息；3）然后通過(guò)五點(diǎn)三次平滑方法去除非相關(guān)噪聲得到信號(hào)的相位波形；4）最后根據(jù)相位波形的躍變特征實(shí)現(xiàn)故障定位。

本文中的信號(hào)處理方法均可通過(guò)MATLAB 7軟件編程實(shí)現(xiàn)[9]。

2 診斷實(shí)例

針對(duì)某鋼廠軋機(jī)設(shè)備主傳動(dòng)系統(tǒng)工況異常，對(duì)軋機(jī)進(jìn)行了現(xiàn)場(chǎng)振動(dòng)測(cè)試，軋機(jī)主傳動(dòng)系統(tǒng)結(jié)構(gòu)如圖3所示?；邶X輪箱的箱體軸承座端是振動(dòng)易測(cè)試部位，就此選用4個(gè)加速度傳感器分別對(duì)齒輪座的輸入軸和輸出軸同時(shí)進(jìn)行振動(dòng)測(cè)試，并對(duì)4個(gè)測(cè)點(diǎn)進(jìn)行編號(hào)，其中1號(hào)、3號(hào)加速度傳感器放置方向?yàn)閺较颍?號(hào)、4號(hào)加速度傳感器放置方向?yàn)檩S向，測(cè)點(diǎn)布置如圖3所示。

齒輪座主要性能參數(shù)：齒輪及軸的轉(zhuǎn)動(dòng)頻率為1.67～4.58 Hz，齒輪嚙合頻率為46.76～128.24 Hz。分析頻率2 000 Hz，采樣點(diǎn)數(shù)8129，采樣長(zhǎng)度0.8 s，N因子5.12，由此可得采樣頻率10 240 Hz。軋機(jī)齒輪座振動(dòng)信號(hào)時(shí)域圖如圖4所示：

圖3 軋機(jī)主傳動(dòng)系統(tǒng)結(jié)構(gòu)簡(jiǎn)圖Fig.3 Mill main transmission system structure diagram

圖4 振動(dòng)信號(hào)的時(shí)域圖Fig.4 Vibration signal time-domain diagram

由圖4所示的時(shí)域信號(hào)可以看出，振動(dòng)信號(hào)有明顯的周期沖擊現(xiàn)象，這說(shuō)明系統(tǒng)已經(jīng)出現(xiàn)了故障，但是否是齒輪座出現(xiàn)了故障，還不能確定，因?yàn)榇斯收闲盘?hào)有可能是其他部件傳播過(guò)來(lái)的，所以需要對(duì)時(shí)域信號(hào)近一步處理。

首先用自適應(yīng)消噪技術(shù)除去強(qiáng)相關(guān)噪聲。根據(jù)前文所介紹的自適應(yīng)消噪處理方法，選擇2號(hào)傳感器測(cè)得的振動(dòng)信號(hào)作為主輸入D(n)，4號(hào)傳感器測(cè)得的振動(dòng)信號(hào)作為參考輸入X(n)，運(yùn)用LMS算法進(jìn)行自適應(yīng)消噪，得到最終輸出Z(n)的時(shí)域圖。如圖5所示：

再對(duì)自適應(yīng)消噪處理后的信號(hào)進(jìn)行希爾伯特變換，得到振動(dòng)信號(hào)的折疊相位，如圖6(a)所示，再進(jìn)行變換得到無(wú)折疊的瞬時(shí)相位，如圖6(b)所示：

因?yàn)樽赃m應(yīng)消噪無(wú)法去除非相關(guān)噪聲，所以圖中有許多毛刺出現(xiàn)。再對(duì)相位信號(hào)進(jìn)行五點(diǎn)三次平滑處理，如圖7 a所示。

用同樣的方法對(duì)4號(hào)傳感器測(cè)得的振動(dòng)信號(hào)進(jìn)行處理，最終得到相位波形如圖7 b所示。

圖5 自適應(yīng)消噪處理后的振動(dòng)信號(hào)時(shí)域圖Fig.5 Vibration signal time-domain diagram after ANC

由圖7可以看出，相位均出現(xiàn)了明顯的超前，相位超前的趨勢(shì)相同，并且兩測(cè)點(diǎn)的最大相位躍變都出現(xiàn)在0.506 6 s處。由此可判定故障應(yīng)發(fā)生在兩傳感器測(cè)試點(diǎn)之外的下一級(jí)上，即齒輪箱輸出軸側(cè)萬(wàn)向接軸或軋輥出現(xiàn)了故障。

在后續(xù)的設(shè)備檢修中重點(diǎn)檢查齒輪箱輸出軸側(cè)部件的異常，最后發(fā)現(xiàn)是萬(wàn)向接軸出現(xiàn)了故障——萬(wàn)向接軸存在間隙導(dǎo)致軋機(jī)齒輪箱的異常振動(dòng)。這與上述分析結(jié)果相符，從而說(shuō)明通過(guò)相位特征傳播特性進(jìn)行故障定位方法的有效性。

圖6 希爾伯特變換得到的信號(hào)相位圖Fig.6 Signal phase diagram after Hilbert transform

圖7 平滑處理后振動(dòng)信號(hào)的相位圖Fig.7 Signal phase diagram after smooth processing

3 結(jié)語(yǔ)

本文提出了基于相位特征傳播特性的故障定位方法，通過(guò)自適應(yīng)消噪技術(shù)和希爾伯特變換提出相位信息，再經(jīng)過(guò)平滑處理得到信號(hào)的相位特征，實(shí)現(xiàn)故障定位；通過(guò)對(duì)軋機(jī)主傳動(dòng)系統(tǒng)故障診斷實(shí)例分析判斷本方法的準(zhǔn)確性，得到如下結(jié)論：

（1）通過(guò)自適應(yīng)消噪技術(shù)對(duì)軋機(jī)信號(hào)的精化，消除了時(shí)變的相關(guān)干擾噪聲的影響，有效地提高了信噪比。五點(diǎn)三次平滑法對(duì)信號(hào)的后處理，剔除了隨機(jī)非相關(guān)噪聲，使得信號(hào)更加接近真實(shí)值，從而可以更容易地發(fā)現(xiàn)故障特征。希爾伯特變換對(duì)軋機(jī)信號(hào)的處理提取出了信號(hào)的相位信息，計(jì)算方法簡(jiǎn)單，易于通過(guò)計(jì)算機(jī)軟件實(shí)現(xiàn)；

（2）針對(duì)單一的信號(hào)處理方法很難準(zhǔn)確地實(shí)現(xiàn)復(fù)雜設(shè)備系統(tǒng)的故障定位這一問(wèn)題，本文提出的利用相位特征傳播特性進(jìn)行故障定位的方法理論較為完善，更加成系統(tǒng)，操作性強(qiáng)，為故障方向定位提供了一條途徑。

[1]佟德純.工程機(jī)械診斷技術(shù)講座[J].工程機(jī)械與維修，2001，1-12.

[2]余紅英，閆宏偉，潘宏俠.齒輪振動(dòng)信號(hào)分解及其在故障診斷中的應(yīng)用[J].振動(dòng)、測(cè)試與診斷，2005 25(2)：109-113.

[3]王秉仁.希爾伯特變換在機(jī)械故障診斷中的應(yīng)用[J].水利電力機(jī)械，2000(1)：43-45.

[4]王 平，廖明夫.齒輪箱中滾動(dòng)軸承的故障診斷[J].航空動(dòng)力學(xué)報(bào)，2005 20(6)：944-950.

[5]Frederick W.King.Hilbert transforms.Cambridge University Press.2008.

[6]孫苗鐘.基于MATLAB的振動(dòng)信號(hào)平滑處理方法[J].電子測(cè)量技術(shù)，2007 30(6)：55-57.

[7]丁 康，李巍華，朱小勇.齒輪及齒輪箱故障診斷實(shí)用技術(shù)[M].北京：機(jī)械工業(yè)出版社，2005.

[8]張新華.齒輪故障識(shí)別[J].機(jī)械強(qiáng)度，1999 21(4)：307-309.

[9]飛思科技產(chǎn)品研究發(fā)中心.MATLAB7輔助信號(hào)處理技術(shù)與應(yīng)用[M].北京：電子工業(yè)出版社，2005.