王 宇 張永祥 姚曉山

（1.海軍工程大學(xué)動(dòng)力工程學(xué)院 武漢 430033）（2.空軍預(yù)警學(xué)院軍械技術(shù)系 武漢 430021）

1 引言

拉缸故障是船用柴油機(jī)極易發(fā)生的故障之一，發(fā)生嚴(yán)重拉缸故障時(shí)，不僅會(huì)使柴油機(jī)損壞，還會(huì)引起人員安全事故。因此柴油機(jī)拉缸故障一直引起人們的高度關(guān)注，并進(jìn)行了廣泛研究。早在二十世紀(jì)八十年代，國內(nèi)外就對(duì)拉缸故障進(jìn)行了理論分析和試驗(yàn)研究。近年來，許多專家學(xué)者從拉缸產(chǎn)生機(jī)理、拉缸的信號(hào)特征、降低拉缸的活塞環(huán)材料選配等方面進(jìn)行了基礎(chǔ)性研究。文獻(xiàn)［1］為降低柴油機(jī)拉缸的風(fēng)險(xiǎn)，研究了不同活塞環(huán)材料對(duì)拉缸的影響，從材料的微觀組織、摩擦系數(shù)上進(jìn)行分析，發(fā)現(xiàn)了拉缸不同階段的機(jī)理。文獻(xiàn)［2］應(yīng)用聲發(fā)射技術(shù)監(jiān)測拉缸故障，發(fā)現(xiàn)在磨損試驗(yàn)機(jī)上可以通過聲發(fā)射信號(hào)和小波包分析監(jiān)測拉缸，但在實(shí)際應(yīng)用上還存在困難。文獻(xiàn)［3］通過磨損試驗(yàn)機(jī)試驗(yàn)研究，揭示了摩擦副材料的金相組織、摩擦力、及摩擦副表面形貌對(duì)拉缸的影響。文獻(xiàn)［4］應(yīng)用三體磨粒磨損機(jī)理模擬發(fā)動(dòng)機(jī)拉缸故障，對(duì)測量的發(fā)動(dòng)機(jī)機(jī)體振動(dòng)信號(hào)進(jìn)行時(shí)域、頻域以及時(shí)頻域分析，發(fā)現(xiàn)拉缸時(shí)將激起2kHz～6kHz頻率上的振動(dòng)；通過分析振動(dòng)信號(hào)的最大值、能量、峭度等參數(shù)，發(fā)現(xiàn)拉缸時(shí)會(huì)有明顯變化。宋力喆［5］等對(duì)柴油機(jī)拉缸進(jìn)行了動(dòng)力學(xué)仿真，研究表明當(dāng)發(fā)生拉缸故障時(shí)，活塞經(jīng)過上、下止點(diǎn)其加速度信號(hào)會(huì)發(fā)生劇烈震動(dòng)，產(chǎn)生沖擊信號(hào)。Hsiao Yeh Chu［6］等通過試驗(yàn)研究了以金剛石納米粒子作為潤滑油添加劑的潤滑油當(dāng)添加劑濃度不同時(shí)其抗拉缸的能力大小。

盡管國內(nèi)外已有眾多學(xué)者對(duì)拉缸故障進(jìn)行分析研究，但直至今日尚未見到用于現(xiàn)場監(jiān)測柴油機(jī)拉缸故障的商業(yè)產(chǎn)品，其根本原因是：一柴油機(jī)是一個(gè)復(fù)雜的往復(fù)式機(jī)械設(shè)備，結(jié)構(gòu)復(fù)雜、振源眾多，且在時(shí)域、頻域與拉缸產(chǎn)生的振動(dòng)、聲信號(hào)重疊，拉缸信號(hào)難以從強(qiáng)背景噪聲中分離出來；二拉缸故障實(shí)機(jī)模擬困難、風(fēng)險(xiǎn)較大。本文為實(shí)現(xiàn)船舶柴油機(jī)拉缸故障診斷，采用已發(fā)生拉缸的6-135柴油機(jī)缸套活塞模擬拉缸故障，通過測量柴油機(jī)機(jī)體的振動(dòng)信號(hào)，采用加窗［7］、EEMD［8］信號(hào)分析方法，識(shí)別診斷柴油機(jī)拉缸故障。

2 加窗和集合經(jīng)驗(yàn)?zāi)B(tài)分解的基本原理

2.1 加窗

柴油機(jī)振動(dòng)激勵(lì)源眾多，對(duì)拉缸故障信號(hào)影響較大的振源主要包括活塞撞擊缸套、氣閥開啟和關(guān)閉、氣缸壓力等，活塞撞擊主要發(fā)生在上下止點(diǎn)附近，由于側(cè)壓力換向而產(chǎn)生。對(duì)于6135柴油機(jī)，某一缸的信號(hào)不僅與自身振動(dòng)激勵(lì)有關(guān)，還受到其他缸振動(dòng)激勵(lì)的影響。圖1為6135柴油機(jī)主要振源激勵(lì)相位圖。

圖1 6135柴油機(jī)主要振源激勵(lì)相位圖

圖中0°為壓縮上止點(diǎn)，進(jìn)氣閥開啟角和排氣閥關(guān)閉角大小為20°，進(jìn)氣閥關(guān)閉角和排氣閥開啟角大小為50°。由圖1可以看出，信號(hào)由于發(fā)生重疊而難以分離。如某缸做功沖程不僅包含自身振動(dòng)源信號(hào)，還受到其他缸的氣閥起落沖擊和活塞撞擊激勵(lì)的影響。因此需要運(yùn)用加窗的方法提取有效信號(hào)段進(jìn)行分析。研究表明：柴油機(jī)拉缸最容易發(fā)生在做功沖程，從圖1可以看出，在柴油機(jī)的壓縮上止點(diǎn)后60°～80°這一段信號(hào)受其他缸影響較小，基本只受到自身激勵(lì)信號(hào)的作用，再考慮到信號(hào)重疊干擾和拉缸信號(hào)獲取的難易程度，所以選擇這一段信號(hào)進(jìn)行柴油機(jī)拉缸故障診斷。為了獲取壓縮上止點(diǎn)后60°～80°這一段信號(hào)，設(shè)計(jì)一矩形窗函數(shù)，保留該段信號(hào)，其他曲軸轉(zhuǎn)角段信號(hào)為零。

2.2 集合經(jīng)驗(yàn)?zāi)B(tài)分解

EEMD是在EMD［9］的基礎(chǔ)上針對(duì)EMD的模態(tài)混疊問題做的改進(jìn)，其基本原理是在EMD的分解過程中加入不同幅值的高斯白噪聲，由于白噪聲是隨機(jī)產(chǎn)生，互不相關(guān)，在多次平均之后噪聲就相互抵消，可以在一定程度上抑制了EMD的模態(tài)混疊問題，進(jìn)而提高了信號(hào)分解的準(zhǔn)確性。對(duì)于EEMD來說最重要的是添加的白噪聲的幅值以及集合平均次數(shù)這兩個(gè)參數(shù)：若是白噪聲幅值過大，會(huì)大大增加計(jì)算量或者使分解信號(hào)的殘余噪聲過多；若是過小則不能起到抑制模態(tài)混疊的作用；集合平均次數(shù)越大也會(huì)增加計(jì)算量，文獻(xiàn)［10］通過利用相關(guān)系數(shù)、相對(duì)均方根誤差準(zhǔn)則和信噪比分析，自適應(yīng)的確定了最優(yōu)的白噪聲幅值系數(shù)K為0.2倍的信號(hào)標(biāo)準(zhǔn)差和集合平均次數(shù)m為50，不但能抑制模態(tài)混疊現(xiàn)象，而且計(jì)算量也較小。EEMD分解的算法流程如下：

1）給原始信號(hào)x（t）加上一組隨機(jī)的白噪聲N1（t），得到總體X1（t）：

2）對(duì)得到的新信號(hào)X1（t）進(jìn)行EMD分解，得到一組基本模態(tài)分量（Intrinsic mode function，IMF）C1j（t），j=1-n。

式中C1j（t）為第一組EMD分解得到的第j個(gè)IMF，r1n為分解后的剩余分量。

3）給原始信號(hào)x（t）加上不同的白噪聲Ni（t），i=2-m，組合后形成新的信號(hào)Xi（t），重復(fù)式（1）、（2）的步驟。

4）對(duì)所得到的m組IMF對(duì)應(yīng)取均值，并作為最終得到的IMF。

3 試驗(yàn)與分析

3.1 試驗(yàn)過程

為了研究柴油機(jī)拉缸故障診斷方法，我們選擇以6-135柴油機(jī)為研究對(duì)象，該柴油機(jī)與一臺(tái)100kW的電力測功器相連，柴油機(jī)啟動(dòng)時(shí)由電力測功器倒拖，在柴油機(jī)工作時(shí)作為發(fā)電機(jī)改變柴油機(jī)負(fù)載。柴油機(jī)拉缸故障模擬非常困難，常用的方法有增大活塞從而減小活塞與缸套間的間隙、調(diào)節(jié)水溫或者在潤滑油內(nèi)添加磨粒，但這些方法不易控制，且拉缸故障風(fēng)險(xiǎn)很大，難以把握，為了真實(shí)模擬柴油機(jī)拉缸故障，選用一套已經(jīng)發(fā)生拉缸的活塞、缸套作為柴油機(jī)拉缸試驗(yàn)件，以確保柴油機(jī)工作時(shí)產(chǎn)生拉缸。為了避免拉缸時(shí)引起柴油機(jī)整機(jī)損壞，試驗(yàn)時(shí)分兩個(gè)階段并逐步加載。第一階段為6-135柴油機(jī)各缸正常試驗(yàn)，測量獲取柴油機(jī)在不同轉(zhuǎn)速、不同負(fù)荷工況時(shí)的振動(dòng)、轉(zhuǎn)速、上止點(diǎn)等信號(hào)。第二階段，首先是由電力測功器倒拖柴油機(jī)，測量振動(dòng)等信號(hào)，然后分別在900r/min、1200r/min和1500r/min時(shí)測量空載和部分負(fù)荷工況時(shí)柴油機(jī)振動(dòng)等信號(hào)。為防止柴油機(jī)產(chǎn)生嚴(yán)重拉缸，采用逐漸增加柴油機(jī)負(fù)荷，且只加載到額定負(fù)荷的20%。為了驗(yàn)證柴油機(jī)各試驗(yàn)工況是否產(chǎn)生拉缸，每次試驗(yàn)后用內(nèi)窺鏡檢驗(yàn)。為了便于柴油機(jī)機(jī)體振動(dòng)信號(hào)的測量，將拉缸的缸套和活塞安裝于柴油機(jī)第一缸。根據(jù)6-135柴油機(jī)的結(jié)構(gòu)，柴油機(jī)拉缸信號(hào)主要由缸套與機(jī)體的上凸緣傳遞至機(jī)體，所以在第一缸活塞上止點(diǎn)對(duì)應(yīng)的機(jī)體位置上粘貼了加速度傳感器轉(zhuǎn)接塊，并在垂直和平行于缸套軸線方向安裝了2個(gè)加速度傳感器，傳感器的型號(hào)為BK4534B，頻寬為0.2kHz～12.8kHz，信號(hào)采集采用的是B & K公司的3053-B-120型數(shù)據(jù)采集模塊，采樣頻率為65536Hz。

在試驗(yàn)過程中我們測取了多個(gè)工況的加速度信號(hào)，由于篇幅所限，這里僅給出部分典型工況的數(shù)據(jù)。圖2為6-135柴油機(jī)900r/min、空載工況正常和拉缸狀態(tài)下的機(jī)體加速度信號(hào)，為了便于觀察將橫坐標(biāo)由時(shí)間轉(zhuǎn)變?yōu)榻嵌取?/p>

從圖2中可以看出多缸機(jī)加速度信號(hào)比較復(fù)雜，主要包含了活塞撞擊激勵(lì)信號(hào)、氣閥起落激勵(lì)信號(hào)以及缸內(nèi)壓力信號(hào)等，信號(hào)重疊嚴(yán)重，很難從圖中直接判斷出柴油機(jī)拉缸故障。根據(jù)文獻(xiàn)［11］和文獻(xiàn)［12］，活塞撞擊信號(hào)和氣閥落座信號(hào)都是高頻信號(hào)，信號(hào)能量主要集中在2kHz～7kHz，最高達(dá)10kHz。

圖2 正常與拉缸狀態(tài)下的加速度信號(hào)

3.2 試驗(yàn)數(shù)據(jù)分析

利用加窗方法對(duì)各工況的加速度信號(hào)進(jìn)行處理，并對(duì)加窗后的加速度信號(hào)進(jìn)行能量計(jì)算和頻譜分析，加窗位置在第一缸壓縮上止點(diǎn)后60°-80°（數(shù)據(jù)點(diǎn)數(shù)為188），數(shù)據(jù)分析長度為50個(gè)工作周期。

圖3為6-135柴油機(jī)在900r/min、空載工況時(shí)正常與拉缸狀態(tài)下的加窗后加速度信號(hào)及其頻譜圖，表1為加窗后對(duì)應(yīng)工況的能量。

圖3 正常與拉缸狀態(tài)下加窗后加速度信號(hào)及頻譜圖

表1 典型工況下拉缸與正常時(shí)加窗信號(hào)能量值

拉缸信號(hào)是高頻信號(hào)，頻帶較寬，從圖3和表1中可以發(fā)現(xiàn)，在同一工況下，柴油機(jī)發(fā)生拉缸故障后，其故障缸的信號(hào)在2kHz-10kHz頻率段的幅值明顯大于該缸在正常狀態(tài)下的，能量也較正常情況時(shí)信號(hào)的能量大，大約為4倍。試驗(yàn)數(shù)據(jù)分析表明：對(duì)于拉缸狀態(tài)下加窗后的能量比正常狀態(tài)下高，高頻信號(hào)能量增加突出，可以作為柴油機(jī)拉缸故障診斷的依據(jù)。

同理，我們對(duì)上述工況的加速度信號(hào)按2至6缸的壓縮上止點(diǎn)后60°-80°分別進(jìn)行加窗處理，并計(jì)算能量和頻譜分析，此時(shí)曲柄轉(zhuǎn)角的0°位置為各缸的壓縮上止點(diǎn)。結(jié)果見圖4和表2。

圖4 900r/min、空載工況時(shí)拉缸狀態(tài)下部分缸加窗后加速度信號(hào)圖及頻譜圖

表2 典型工況下拉缸狀態(tài)時(shí)各缸的加窗信號(hào)能量

試驗(yàn)數(shù)據(jù)分析表明：柴油機(jī)單缸發(fā)生拉缸故障時(shí)該缸加速度信號(hào)在2kHz～10kHz頻率段的幅值以及信號(hào)能量都大于同一工況下其他未拉缸的信號(hào)，能量大約是未拉缸信號(hào)的10倍。所以，可以通過對(duì)加窗信號(hào)的頻譜分析和能量計(jì)算初步實(shí)現(xiàn)故障缸的定位。

為了進(jìn)一步對(duì)拉缸故障進(jìn)行精確診斷，對(duì)1200r/min、100N·m工況下拉缸狀態(tài)加窗后的加速度信號(hào)分別進(jìn)行EMD和EEMD分解，圖5（a）、（b）分別為經(jīng)EMD和EEMD分解后的前三階IMF。

圖5 加窗后加速度信號(hào)EMD和EEMD的分解結(jié)果

表3 典型工況下正常和拉缸時(shí)加窗信號(hào)前3階IMF能量

從圖5可以看出，信號(hào)經(jīng)EMD分解的IMF與實(shí)際的拉缸信號(hào)相比相差較大，而EEMD分解的IMF與實(shí)際拉缸信號(hào)較為一致。因?yàn)椴裼蜋C(jī)振源眾多，拉缸信號(hào)受影響因素較多，導(dǎo)致EMD在分解時(shí)容易產(chǎn)生模態(tài)混疊現(xiàn)象，故EMD不適合處理拉缸信號(hào)。而EEMD可以克服EMD模態(tài)混疊，能夠用于柴油機(jī)拉缸故障診斷。依次將加窗后的加速度信號(hào)進(jìn)行EEMD分解，分解得到16個(gè)IMF分量和一個(gè)殘差，分析發(fā)現(xiàn)對(duì)于拉缸信號(hào)影響最大的是前三階IMF。計(jì)算前三階IMF的能量值并對(duì)由前3階IMF重構(gòu)的信號(hào)進(jìn)行頻譜分析，具體結(jié)果如表3和圖6。

圖6 1200r/min、100N·m工況下拉缸與正常時(shí)重構(gòu)信號(hào)的頻譜圖

從表3和圖6可以發(fā)現(xiàn)：一旦柴油機(jī)發(fā)生拉缸故障，則故障缸加速度信號(hào)加窗后經(jīng)EEMD分解的前三階IMF能量都比同一工況下未發(fā)生拉缸時(shí)的能量大，大約為6～9倍，且由前3階IMF重構(gòu)信號(hào)的頻譜圖可以發(fā)現(xiàn)，柴油機(jī)拉缸時(shí)故障缸重構(gòu)信號(hào)的高頻信號(hào)明顯增加，大約增加2-3倍，主要增加段為8.5kHz～11kHz和13kHz～16kHz段。所以通過對(duì)加窗后的信號(hào)進(jìn)行EEMD分解，并計(jì)算前三階IMF的能量值以及對(duì)前3階IMF重構(gòu)信號(hào)進(jìn)行頻譜分析，能夠?qū)崿F(xiàn)柴油機(jī)拉缸故障的診斷。

同理，我們對(duì)1缸發(fā)生拉缸時(shí)的其余2至6缸加速度信號(hào)加窗后進(jìn)行EEMD分解，得到其前3階IMF能量分布以及前3階IMF重構(gòu)信號(hào)頻譜圖，如表4和圖7。

表4 1200r/min、100N·m工況下拉缸時(shí)加窗信號(hào)前3階IMF能量值

試驗(yàn)數(shù)據(jù)分析表明：當(dāng)柴油機(jī)發(fā)生拉缸故障時(shí)，故障缸的加速度信號(hào)經(jīng)加窗和EEMD分解后的前三階IMF的能量都比其他缸信號(hào)的大，大約為10-20 倍，重構(gòu)信號(hào)在 8.5kHz～11kHz和 13kHz～16kHz段的幅值明顯增加，大約增加2-3倍。所以在柴油機(jī)發(fā)生拉缸故障時(shí)，可以通過計(jì)算加窗和EEMD分解的前三階IMF能量以及對(duì)重構(gòu)信號(hào)頻譜分析實(shí)現(xiàn)故障缸的定位。

圖7 1200r/min、100N·m工況時(shí)拉缸狀態(tài)下部分缸重構(gòu)信號(hào)頻譜圖

綜上所述，當(dāng)船用柴油機(jī)發(fā)生拉缸故障時(shí)，對(duì)加窗后加速度信號(hào)進(jìn)行EEMD分解，在計(jì)算前3階IMF的能量及對(duì)重構(gòu)信號(hào)進(jìn)行頻譜分析，能夠更加準(zhǔn)確有效地實(shí)現(xiàn)拉缸故障診斷。

4 結(jié)語

1）采用已產(chǎn)生拉缸的活塞缸套裝入柴油機(jī)某缸進(jìn)行拉缸試驗(yàn)研究，是一種行之有效的方法。

2）合理選取加窗參數(shù)，可有效提取柴油機(jī)的拉缸故障信號(hào)，減小其他振源對(duì)拉缸故障信號(hào)的影響。

3）船用柴油機(jī)發(fā)生拉缸時(shí)，其加窗段的信號(hào)能量要比不拉缸時(shí)的能量大；拉缸缸對(duì)應(yīng)的加窗信號(hào)能量要比其他未發(fā)生拉缸缸對(duì)應(yīng)的加窗信號(hào)能量大。拉缸缸經(jīng)EEMD分解后的重構(gòu)信號(hào)在8.5kHz～11kHz和13kHz～16kHz段的幅值明顯增大，據(jù)此可實(shí)現(xiàn)船用多缸柴油機(jī)某缸拉缸故障的判別和定位。

4）與EMD相比，EEMD可以克服其模態(tài)混疊現(xiàn)象，而且利用前3階模態(tài)（IMF）的能量分布以及重構(gòu)信號(hào)的頻譜分析可以更加有效地進(jìn)行船用柴油機(jī)拉缸故障的診斷。