用經(jīng)驗(yàn)?zāi)Ｊ椒纸夂椭貥?biāo)極差分析法進(jìn)行摩擦振動(dòng)信號(hào)分析

李精明+魏海軍+梅立強(qiáng)+許德志

摘要：

為提取和表征摩擦振動(dòng)的特征信號(hào)，在摩擦磨損試驗(yàn)機(jī)上進(jìn)行以船用柴油機(jī)缸套和活塞環(huán)為材質(zhì)的摩擦副摩擦磨損試驗(yàn).應(yīng)用經(jīng)驗(yàn)?zāi)Ｊ椒纸猓‥mpiricalModeDecomposition，EMD）對(duì)非線性、非平穩(wěn)的摩擦振動(dòng)信號(hào)進(jìn)行分解，獲得若干個(gè)本征模式函數(shù)（IntrinsicModeFunction，IMF）.從中選擇反映摩擦振動(dòng)特征的IMF重新合成摩擦振動(dòng)特征信號(hào)，應(yīng)用重標(biāo)極差分析法對(duì)摩擦振動(dòng)特征信號(hào)進(jìn)行分析，得到線性回歸譜和Hurst指數(shù).結(jié)果表明，EMD能夠?qū)崿F(xiàn)摩擦振動(dòng)特征信號(hào)的提取，重標(biāo)極差分析法可以分析摩擦振動(dòng)信號(hào)的漸變過程，提取摩擦振動(dòng)信號(hào)的特征.該方法可為基于摩擦振動(dòng)信號(hào)的機(jī)械摩擦副摩擦磨損行為的研究提供新的途徑.

關(guān)鍵詞：

經(jīng)驗(yàn)?zāi)Ｊ椒纸猓‥MD）；重標(biāo)極差分析；Hurst指數(shù)；摩擦振動(dòng)

中圖分類號(hào)：U664.121；TH117.1

文獻(xiàn)標(biāo)志碼：A 收稿日期：20150814 修回日期：20151113

0引言

摩擦振動(dòng)是機(jī)械運(yùn)動(dòng)摩擦副在摩擦磨損過程中產(chǎn)生的現(xiàn)象，摩擦振動(dòng)信號(hào)中包含了許多反映系統(tǒng)摩擦學(xué)特征和摩擦狀態(tài)的信息[1]，如摩擦力矩、摩擦因數(shù)、磨損表面形貌、磨粒形貌[2]等.許多學(xué)者從這些摩擦學(xué)參數(shù)入手來研究摩擦磨損狀態(tài)，但通過摩擦力矩、摩擦因數(shù)、磨損表面形貌來提取摩擦學(xué)特征存在一定的困難[3]，而通過磨粒形貌進(jìn)行的分析程序繁瑣、低效，且與研究人員的經(jīng)驗(yàn)關(guān)系很大[4].機(jī)械設(shè)備摩擦振動(dòng)信號(hào)的獲取容易，可在不影響機(jī)械設(shè)備運(yùn)行的情況下實(shí)時(shí)在線采集，因此摩擦振動(dòng)的研究是摩擦副的在線監(jiān)測(cè)和狀態(tài)識(shí)別的基礎(chǔ).

經(jīng)驗(yàn)?zāi)Ｊ椒纸猓‥mpiricalModeDecomposition，EMD）是針對(duì)非線性、非平穩(wěn)信號(hào)提出的一種方法[5]，能自適應(yīng)地把信號(hào)按照自身的內(nèi)在特性正交地分解成一系列本征模函數(shù)（IntrinsicModeFunction，IMF）.該方法已成功地應(yīng)用于信號(hào)處理[6]、機(jī)械故障診斷[79]等的研究.重標(biāo)極差分析法是英國(guó)水文學(xué)家HURST在1951年提出的，它是一種無參數(shù)統(tǒng)計(jì)方法，可以判別數(shù)據(jù)是否具有長(zhǎng)期非函數(shù)周期時(shí)間序列與隨機(jī)序列特性[10].該方法被應(yīng)用到巖土力學(xué)[11]、股市分析[12]等領(lǐng)域，取得了很好的效果.許多學(xué)者從定量的角度，應(yīng)用時(shí)頻譜圖[13]、諧波小波[14]、混沌吸引子[15]等提取摩擦振動(dòng)信號(hào)特征，取得了一定的成果.本文對(duì)以船用柴油機(jī)缸套和活塞環(huán)為材質(zhì)的摩擦副在摩擦磨損試驗(yàn)中獲得的摩擦振動(dòng)信號(hào)進(jìn)行EMD處理，獲得若干個(gè)IMF分量，再從中選擇反映摩擦振動(dòng)特征的IMF分量重新合成摩擦振動(dòng)特征信號(hào).應(yīng)用重標(biāo)極差分析法對(duì)摩擦振動(dòng)特征信號(hào)進(jìn)行分析，得到線性回歸譜和Hurst指數(shù)，實(shí)現(xiàn)特征參數(shù)對(duì)摩擦副摩擦磨損狀態(tài)的表征，為基于摩擦振動(dòng)信號(hào)的機(jī)械摩擦副摩擦磨損行為的研究提供了新的途徑.

1試驗(yàn)部分

1.1試驗(yàn)方法與材料

試驗(yàn)設(shè)備采用CFTI型摩擦磨損試驗(yàn)機(jī)，圖1為該試驗(yàn)機(jī)的原理.銷試樣和盤試樣是采用線切割

方法從船用柴油機(jī)的活塞環(huán)和缸套上截取的.銷

試樣固定在主軸下端，盤試樣固定在臺(tái)架上，電機(jī)驅(qū)動(dòng)的偏心輪連桿機(jī)構(gòu)使臺(tái)架和盤試樣做往復(fù)運(yùn)動(dòng).載荷由加載機(jī)構(gòu)經(jīng)銷試樣施加于盤試樣，振動(dòng)信號(hào)采用美國(guó)壓電公司生產(chǎn)的356A16型三軸加速度傳感器測(cè)量，每個(gè)加速度單位（g）輸出100mV，量程為±50g，檢測(cè)位置為盤試樣下方.銷試樣材質(zhì)為合金鑄鐵，主要成分為Fe，C，Si，Mo，P等，硬度HV600～680GPa，原始表面粗糙度Ra=0.693μm，截面尺寸3mm×4mm.盤試樣接觸表面為弧面，材質(zhì)為合金鑄鐵，主要成分為Fe，C，Si，Mn，P等，硬度HV320～420GPa，原始表面粗糙度Ra=1.652μm，截面圓直徑為30mm，高為10mm.施加的力為50N，驅(qū)動(dòng)電機(jī)轉(zhuǎn)速為600r/min，銷試樣盤試樣相對(duì)運(yùn)動(dòng)的平均速度為0.1m/s，試驗(yàn)時(shí)間為600min.潤(rùn)滑方式為滴油潤(rùn)滑，潤(rùn)滑介質(zhì)為MobilgardM440船用潤(rùn)滑油，黏度142cSt（40℃）.

1.2摩擦振動(dòng)信號(hào)的采集

摩擦振動(dòng)信號(hào)的采集裝置為比利時(shí)LMS公司生產(chǎn)的SCADAS型前端數(shù)據(jù)采集系統(tǒng).該系統(tǒng)采集的是法向的摩擦振動(dòng)信號(hào)，采樣頻率為25600Hz.每0.16s所采集的數(shù)據(jù)會(huì)自動(dòng)生成一個(gè)文本文件，存入計(jì)算機(jī).

圖2為摩擦磨損試驗(yàn)在初期、中期、末期采集的摩擦振動(dòng)信號(hào)時(shí)域波形，

3個(gè)子圖中所用數(shù)據(jù)分別為在試驗(yàn)的

第1，第300，第600個(gè)1min內(nèi)選取的一個(gè)0.16s的數(shù)據(jù).從圖2可以看出，所獲得的摩擦振動(dòng)信號(hào)是非線性、非平穩(wěn)信號(hào)，無規(guī)律可循，微弱的摩擦振動(dòng)信號(hào)湮沒于背景噪聲之中，需要去噪才能提取正確的摩擦振動(dòng)特征.

2摩擦振動(dòng)信號(hào)的EMD處理

EMD被認(rèn)為是以傅里葉變換為基礎(chǔ)的線性和穩(wěn)態(tài)頻譜分析的一個(gè)重大突破.該方法依據(jù)數(shù)據(jù)自身的時(shí)間尺度特征進(jìn)行信號(hào)分解，可將任意信號(hào)分解為若干個(gè)本征模式分量與一個(gè)余項(xiàng)之和

，在處理非線性、非平穩(wěn)數(shù)據(jù)方面具有非常明顯的優(yōu)勢(shì)

.一個(gè)IMF必須滿足以下兩個(gè)條件：①在整個(gè)時(shí)間范圍內(nèi)，函數(shù)局部極值點(diǎn)的數(shù)目與過零點(diǎn)的數(shù)目必須相等，或最多相差1個(gè)；②任意時(shí)刻點(diǎn)，函數(shù)局部最大值的包絡(luò)（上包絡(luò)線）與局部最小值的包絡(luò)（下包絡(luò)線）的平均值必須為零.算法實(shí)現(xiàn)步驟如下：

（1）設(shè)

待分析信號(hào)為x（t），首先找出信號(hào)的所有極大值和極小值，然后采用三次樣條插值算法分別獲得信號(hào)的上包絡(luò)線xH（t）和下包絡(luò)線xL（t），信號(hào)所有數(shù)據(jù)點(diǎn)位于上、下包絡(luò)線之間.計(jì)算出兩條包絡(luò)線的平均值

若x1（r）滿足IMF的兩個(gè)標(biāo)準(zhǔn)，則可以作為第一個(gè)IMF并記為c1（t），否則作為原始信號(hào)，重復(fù)步驟（1）和（2），直到得到第一個(gè)IMF分量.

（3）計(jì)算剩余信號(hào)r1（t）=x（t）-c1（t），將其作為新的信號(hào)，重復(fù)步驟（1）和（2），直到提取出信號(hào)的所有IMF分量.當(dāng)剩余的信號(hào)為非振蕩的單調(diào)函數(shù)或小于預(yù)定值的常數(shù)時(shí)，認(rèn)為分解完畢.原始的數(shù)據(jù)最終分解為n個(gè)IMF和一個(gè)剩余的rn（t），即

應(yīng)用EMD對(duì)摩擦振動(dòng)信號(hào)進(jìn)行分析.圖3是對(duì)摩擦磨損試驗(yàn)初期的摩擦振動(dòng)信號(hào)運(yùn)用EMD得到的8個(gè)IMF分量c1～c8和1個(gè)殘差r8（限于篇幅，未列出摩擦振動(dòng)信號(hào)中期、末期的EMD結(jié)果）.從圖3

可以看出，得到的各個(gè)IMF分量包含了不同的尺

度，可使摩擦振動(dòng)信號(hào)的特征在不同的分辨率下顯現(xiàn)出來，高頻率和小振幅集中出現(xiàn)在前幾個(gè)分量中.EMD是一種主成分提取方法，分解結(jié)果中前幾個(gè)IMF分量包含原始信號(hào)的主要信息.摩擦振動(dòng)信號(hào)具有振幅小、頻率高等特征[14]，分解得到的IMF分量c1，c2和c3滿足這些特征，故從摩擦振動(dòng)信號(hào)EMD結(jié)果中選取c1，c2和c3重新合成試驗(yàn)初期的摩擦振動(dòng)特征信號(hào).試驗(yàn)中期、末期的摩擦振動(dòng)信號(hào)也按此方法合成，得到去噪的摩擦振動(dòng)特征信號(hào).圖4為摩擦磨損試驗(yàn)初期、中期、末期信號(hào)經(jīng)EMD處理后適當(dāng)選取IMF分量合成的摩擦振動(dòng)特征信號(hào).從圖4可以看出，IMF分量合成的信號(hào)幅值呈現(xiàn)逐漸減小的變化趨勢(shì)，呈現(xiàn)出清晰的摩擦振動(dòng)的沖擊信息.

3摩擦振動(dòng)信號(hào)重標(biāo)極差分析法

3.1重標(biāo)極差分析法理論及算法

3.2摩擦振動(dòng)信號(hào)Hurst指數(shù)分析

應(yīng)用重標(biāo)極差分析法對(duì)經(jīng)EMD去噪后獲得的摩擦振動(dòng)特征信號(hào)進(jìn)行分析，摩擦振動(dòng)數(shù)據(jù)采集系統(tǒng)在0.16s內(nèi)采集的信號(hào)有4096個(gè)，均分序列的長(zhǎng)度A取16～1024均分的19個(gè)長(zhǎng)度.圖5為應(yīng)用R/S分析法分析摩擦振動(dòng)初期、中期、末期信號(hào)得到的雙對(duì)數(shù)線性回歸譜，圖中星點(diǎn)、圓點(diǎn)、方塊點(diǎn)分別為計(jì)算摩擦振動(dòng)初期、中期、末期信號(hào)得到的R/S的對(duì)數(shù)值，直線1，2，3為對(duì)應(yīng)信號(hào)的回歸線.從圖中可以看出，初期、中期、末期回歸線的斜率（即Hurst指數(shù)）分別為0.1589，0.1886，0.2425，均在0到0.5之間.這表明摩擦磨損試驗(yàn)獲得的摩擦振動(dòng)信號(hào)是反持久性序列，摩擦振動(dòng)未來的趨勢(shì)表現(xiàn)為與之前的相反.從初期到末期，信號(hào)回歸線的斜率呈現(xiàn)逐漸增大的變化，表明反持續(xù)性逐漸減小，序列的突變跳躍逆轉(zhuǎn)性趨于緩和.

摩擦振動(dòng)的Hurst指數(shù)變化與磨合過程摩擦副磨損表面的變化密切相關(guān).磨合初期，銷試樣和盤試樣的表面粗糙度較大，銷試樣Ra=0.693μm，盤試樣Ra=1.652μm，相對(duì)運(yùn)動(dòng)過程中兩表面磨損嚴(yán)重，摩擦振動(dòng)劇烈.從圖4中摩擦振動(dòng)初期IMF分量合成的摩擦振動(dòng)特征信號(hào)可以看出，初期的摩擦振動(dòng)主峰明顯，平均振幅較大.這表明摩擦副相對(duì)運(yùn)動(dòng)產(chǎn)生的能量較大，摩擦磨損處于不穩(wěn)定狀態(tài)，Hurst指數(shù)最小，反持續(xù)性最強(qiáng).隨著磨合的進(jìn)行，摩擦副表面的粗糙度逐漸下降，激發(fā)的摩擦振動(dòng)也減弱，磨合過程中摩擦表面消耗的能量減弱，摩擦磨損逐漸趨于平衡，Hurst指數(shù)增大.磨合后期，摩擦副進(jìn)入穩(wěn)定磨損狀態(tài)，摩擦振動(dòng)主峰不明顯，平均振幅減小，Hurst指數(shù)最大，反持續(xù)性減弱.

3.3改變?cè)囼?yàn)條件的摩擦振動(dòng)信號(hào)分析

為驗(yàn)證重標(biāo)極差分析法提取摩擦振動(dòng)信號(hào)特征的有效性，改變?cè)囼?yàn)條件，另制備一對(duì)銷試樣和盤試樣，試驗(yàn)方法與材料與前述試驗(yàn)相同，測(cè)得銷試樣和盤試樣的原始表面粗糙度Ra分別為0.627μm和1.586μm，施加力30N，試驗(yàn)時(shí)間480min.按相同的數(shù)據(jù)分析方法得到的線性回歸譜見圖6.從圖6可以看出，初期、中期、末期信號(hào)的Hurst指數(shù)分別為0.2065，0.2090，0.2641.這表明Hurst指數(shù)隨試驗(yàn)條件的變化而變化，能夠敏感地反映試驗(yàn)條件的變化.因此，可以將本方法用于機(jī)械故障診斷，比如柴油機(jī)活塞環(huán)斷裂、活塞環(huán)與缸套間潤(rùn)滑不良等故障的檢測(cè).

以上分析表明，摩擦學(xué)系統(tǒng)為非線性耗散的動(dòng)力系統(tǒng)，可以用重標(biāo)極差分析法線性回歸譜和Hurst

指數(shù)來表征摩擦振動(dòng)信號(hào)的特征，識(shí)別摩擦振動(dòng)的狀態(tài).

4結(jié)論

摩擦振動(dòng)信號(hào)包含摩擦副磨合磨損的狀態(tài)信息，本文利用經(jīng)驗(yàn)?zāi)Ｊ椒纸猓‥MD）和重標(biāo)極差分析法研究摩擦副摩擦磨損過程中的摩擦振動(dòng)信號(hào)，結(jié)論為：

（1）應(yīng)用EMD對(duì)試驗(yàn)獲得的摩擦振動(dòng)原始信號(hào)進(jìn)行分解，根據(jù)摩擦振動(dòng)的特征適當(dāng)選擇本征模式函數(shù)（IMF）分量進(jìn)行重構(gòu)，可以獲得反映摩擦副摩擦振動(dòng)特征的特征信號(hào).

（2）應(yīng)用重標(biāo)極差分析法能有效地分析摩擦副摩擦振動(dòng)的漸變過程，可以用重標(biāo)極差分析得到的線性回歸譜及Hurst指數(shù)表征摩擦振動(dòng)信號(hào)的特征.

參考文獻(xiàn)：

[1]

葛世榮，朱華.摩擦學(xué)復(fù)雜系統(tǒng)及其問題的量化研究方法[J].摩擦學(xué)學(xué)報(bào)，2002，22（5）：405408.

[2]劉燁，王永彬，王紅麗，等.基于統(tǒng)計(jì)和分形理論的氣缸套磨損表面形貌表征[J].車用發(fā)動(dòng)機(jī)，2013，205（2）：8992.

[3]ZHUHua，GEShirong，CAOXucheng.Thechangesoffractaldimensionsoffrictionalsignalsintherunninginwearprocess[J].Wear，2007，263（7）：15021507.

[4]徐金龍，粟斌，程琛.光譜和鐵譜及自動(dòng)磨粒激光技術(shù)在油液監(jiān)測(cè)中的應(yīng)用[J].潤(rùn)滑與密封，2009，34（9）：105109.

[5]HUANGNE，SHENZ，LONGSR，etal.TheempiricalmodedecompositionandtheHilbertspectrumfornonlinearandnonstationarytimeseriesanalysis[J].ProceedingoftheRoyalSociety：AMathematicalPhysical&EngineeringSciences，1998，454（1）：903995.

[6]GUOKejian，ZHANGXingang，LIHongguang，etal.ApplicationofEMDmethodtofrictionsignalprocessing[J].MechanicalSystemsandSignalProcessing，2008，22（1）：284259.

[7]LIYujun，PWTse，YANGXin，etal.EMDbasedfaultdiagnosisforabnormalclearancebetweencontactingcomponentsinadieselengine[J].MechanicalSystemsandSignalProcessing，2010，24（1）：193210.

[8]余德介，程軍圣.EMD方法在齒輪故障診斷中的應(yīng)用[J].湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2002，29（6）：4851.

[9]江國(guó)和，趙開琦，王志剛，等.基于EMD的Hilbert變換的柴油機(jī)氣缸套磨損故障診斷[J].上海海事大學(xué)學(xué)報(bào)，2014，35（3）：8084.DOI：10.13340/j.jsmu.2014.03.015.

[10]HURSTHE.Longtermstorageofreservoirs[J].TransactionsofAmericanSocietyofCivilEngineers，1951，116（3）：770808.

[11]賀可強(qiáng)，孫林娜，王思敬.滑坡位移分形參數(shù)Hurst指數(shù)及其在堆積層滑坡預(yù)報(bào)中的應(yīng)用[J].巖土力學(xué)與工程學(xué)報(bào)，2009，28（6）：11071115.

[12]莊新田，莊新路，田瑩.Hurst指數(shù)及股市的分形結(jié)構(gòu)[J].東北大學(xué)學(xué)報(bào)（自然科學(xué)版），2003，24（9）：861965.

[13]黃朝明，于洪亮，關(guān)德林，等.摩擦振動(dòng)時(shí)頻圖像特征提取[J].振動(dòng)與沖擊，2012，31（7）：4549.

[14]李國(guó)賓，任宗英，王宏志，等.摩擦振動(dòng)信號(hào)諧波小波包特征提取[J].摩擦學(xué)學(xué)報(bào)，2011，31（5）：452456.

[15]SUNDi，LIGuobin，WEIHaijun，etal.Experimentalstudyonthechaoticattractorevolvementofthefrictionvibrationinarunninginprocess[J].TribologyInternational，2015（88）：290297.