劉 琦 林瑞霖
(海軍工程大學船舶與動力學院 武漢 430033)
柴油機的缸內(nèi)壓力直接反映柴油機的工作性能,包括進氣過程的好壞、噴油定時是否準確、氣閥間隙是否合適、氣缸的氣密性的好壞、燃料燃燒的完整程度等[1]。因此,測取壓力信號可以對柴油機的工作狀態(tài)進行檢測。通過缸內(nèi)壓力信號對柴油機進行故障診斷,主要是通過從缸內(nèi)壓力信號中提取出缸內(nèi)溫度、壓力升高率、燃燒放熱率等隨曲軸轉(zhuǎn)角的變化規(guī)律[2]。但是,通過壓力傳感器測取的原始缸內(nèi)壓力信號是混有很多噪聲的電壓信號,能否有效地從中提取出隨曲軸轉(zhuǎn)角變化的比較光滑的壓力信號,對判斷柴油機工作是否正常和診斷柴油機的故障類型至關(guān)重要。
1)等曲軸轉(zhuǎn)角化
對壓力信號的采樣是通過壓力傳感器所采集的等時間間隔的信號,而對柴油機工作狀態(tài)的分析是通過等曲軸轉(zhuǎn)角的壓力信號實現(xiàn)的[3],所以,必須對等時間間隔采樣的信號進行等曲軸轉(zhuǎn)角處理。壓力傳感器在柴油機曲軸的各個輪齒轉(zhuǎn)動過程中所采集的信號個數(shù)為
其中,Z為柴油機曲軸一圈輪齒的個數(shù);fs為壓力傳感器的采樣頻率;ni為曲軸在經(jīng)過不同輪齒時的瞬時轉(zhuǎn)速;x(Zi)為不同輪齒時壓力傳感器所采集的信號個數(shù)。
根據(jù)各個輪齒間采樣點的信號,通過Newten插值公式求出等曲軸轉(zhuǎn)角的壓力傳感器信號,Newten插值算法如下[4]:
采集信號的精度取決于壓力傳感器的頻率和柴油機轉(zhuǎn)速,轉(zhuǎn)速越低,采樣頻率越高,信號的準度越高。
2)多循環(huán)平均
柴油機的做功過程是通過燃料的燃燒產(chǎn)生高溫高壓氣體,推動活塞做功。對于四沖程柴油機,每一個做功循環(huán)都經(jīng)歷吸氣、壓縮、燃燒做功和排氣沖程,根據(jù)這個順序周而復始地進行工作,每個工作循環(huán)都經(jīng)歷度曲軸轉(zhuǎn)角[5]。即使在穩(wěn)定工況下,柴油機工作中的各個不同循環(huán),缸內(nèi)壓力信號也會出現(xiàn)波動[6]。因此,為了減小隨機誤差和噪聲干擾,必須對測取的缸內(nèi)壓力信號進行多循環(huán)平均。多循環(huán)平均的算法如下:其中,θ為曲軸轉(zhuǎn)角;Vi(θ)為隨曲軸轉(zhuǎn)角變化的壓力傳感器信號;N為循環(huán)數(shù);ˉV(θ)為平均化處理后的壓力傳感器信號。
3)標定
壓力傳感器輸出的是電壓信號,而要得到缸內(nèi)壓力信號,需要對壓力傳感器進行標定,根據(jù)標定結(jié)果,將采集的電壓信號轉(zhuǎn)化為壓力信號。壓力傳感器的輸出信號隨壓力的大小呈線性變化,利用最小二乘法對壓力傳感器進行標定,標定算法為,設
求滿足
的α、b,其中α、b分別為壓力傳感器的放大倍數(shù)和零點漂移,V、P(V)分別為壓力傳感器輸出的電壓信號和所要得出的缸內(nèi)壓力信號。
4)光順處理
壓力傳感器采集的缸內(nèi)壓力信號之中混有很多高頻的噪聲干擾,即使進行了多個循環(huán)的平均,壓力信號的毛刺依然很大;計算放熱率要求P-θ曲線是可導的,但是由于壓力值采集誤差,以及其他誤差等原因,常常使壓力數(shù)據(jù)并不滿足可導的條件。由于以上原因,用數(shù)學方法計算的dQ-dθ曲線一般會產(chǎn)生較大的震蕩[7]。所以,必須將平均化以后的壓力信號進行光順處理。
光順處理的方法很多,可采用樣條函數(shù)、最小二乘法、正交多項式法、五點三次平滑法等。五點三次平滑法的算法如下[8]:
以上即為對壓力傳感器所采集的原始信號的處理流程。
本實驗是在6135非增壓四沖程柴油機上進行的,轉(zhuǎn)速置為1000r/min,負載400Nm,6135柴油機 的 發(fā) 火 順 序 為 1-5-3-6-2-4,壓 力 傳 感 器 采 用HEM375耐高溫壓力傳感器,采樣頻率為65536Hz,實驗采集了2號缸40個循環(huán)的壓力信號進行分析。
1)在轉(zhuǎn)速1000r/min,負載400Nm的工況下,采集的40個循環(huán)的缸內(nèi)壓力原始數(shù)據(jù)和單個循環(huán)的缸內(nèi)壓力原始數(shù)據(jù)分別如圖1和圖2。從圖1可以看出,即使在穩(wěn)定工況下,各個不同循環(huán)的缸內(nèi)壓力信號波動依然很大,其中最高爆發(fā)壓力的波動達到31.6%。單缸的缸內(nèi)壓力信號在不同做功階段受噪聲干擾均比較大,其中以低壓段的表現(xiàn)最為明顯,如圖2,這主要由于高壓段壓力信號占主導,而低壓段噪聲信號的比例相對高壓段高。
圖1 40個循環(huán)的缸內(nèi)壓力信號
圖2 單個循環(huán)的缸內(nèi)壓力信號
2)按照等曲軸轉(zhuǎn)角化的算法,將等采樣頻率的壓力傳感器信號轉(zhuǎn)化為等曲軸轉(zhuǎn)角的信號,結(jié)果如圖3,可以看出,經(jīng)過等曲軸轉(zhuǎn)角化后,更加能清楚地反映柴油機在不同做功階段的壓力變化規(guī)律。
圖3 等曲軸轉(zhuǎn)角化的壓力傳感器信號
圖4 40個循環(huán)平均的壓力傳感器信號
3)按照多循環(huán)平均算法,對40個循環(huán)的壓力傳感器信號作平均化處理,結(jié)果如圖4,可以看出,經(jīng)過多循環(huán)平均處理后,信號的隨機誤差明顯減小,信號相對平滑。
4)按照光順算法,分別對平均化以后的壓力傳感器信號進行一次、二次、三次光順處理,如圖5??梢钥吹剑S著光順次數(shù)的增加,信號越來越平滑[9],但是三次光順以后,信號在低壓段(包括壓縮段)的噪聲干擾仍然很大。為有效減少信號在低壓段的噪聲干擾,對信號進行六次光順,信號在低壓和高壓段均比較平滑,但卻使信號的一些特征點不明顯,如壓縮拐點消失,信號的峰值誤差也加大,信號的峰值壓力誤差達到5%,如圖5。針對壓力信號在高壓段有用信號占主導,而在低壓段噪聲干擾大的特點,本文選擇對高壓段和低壓段分別采用三次光順和六次光順的方法進行處理,結(jié)果如圖6,可以看到,處理后的結(jié)果既保持了原始信號的特征點,又使峰值的誤差比較小,峰值壓力誤差為0.88%。
圖5 不同光順次數(shù)的缸內(nèi)壓力信號對比
5)按照最小二乘算法,對壓力傳感器進行標定。所采用的標定點如表1。
最小二乘法標定曲線如圖7,由圖可見,壓力傳感器的線性度比較好,標定的結(jié)果為,
P(V)=5.2193V+1.0238 (7)
利用式(7)的標定結(jié)果,對壓力傳感器所測量的壓力信號進行標定,得到最終的缸內(nèi)壓力隨曲軸轉(zhuǎn)角變化的結(jié)果,如圖8所示。
表1 標定點的選取
圖6 分段光順的缸內(nèi)壓力信號
圖7 最小二乘法標定曲線
圖8 壓力曲線標定結(jié)果
6)本文對比了不同光順方法對壓力升高率計算的影響,如圖9,可以看出,壓力信號未經(jīng)過光順處理,得出的壓力升高率曲線毛刺比較大[10],經(jīng)過五點三次光順后,壓力升高率曲線高壓段毛刺較小,而低壓段毛刺依然很大,分段光順后,以壓縮終點為零度轉(zhuǎn)角,計算的缸內(nèi)壓力升高率曲線在低壓和高壓段均比較平滑,并且根據(jù)分段光順的壓力曲線,由理想氣體狀態(tài)方程計算的壓縮段缸內(nèi)溫度曲線也比較平滑。
圖9 不同光順方法對壓力升高率計算的影響以及根據(jù)分段光順壓力計算的壓縮階段溫度
1)對壓力傳感器測量的原始缸內(nèi)壓力信號經(jīng)過等曲軸轉(zhuǎn)角化、多循環(huán)平均、光順處理和標定以后,可以得到滿足故障診斷需要的缸內(nèi)壓力數(shù)據(jù)。
2)采用分段光順方法對壓力數(shù)據(jù)進行處理可以使壓力曲線在高壓和低壓段均比較平滑,而且保持了壓力曲線的特征點,減小了峰值誤差。
[1]郭江華,侯馨文,陳國鈞,等.船舶柴油機故障診斷技術(shù)研究[J].中國航海,2005(4):75
[2]黃碧華,裘崇偉,謝友柏.柴油機磨損狀態(tài)監(jiān)測及故障診斷專家系統(tǒng)知識庫建立的研究[J].摩擦學學報,1994,45:36
[3]潘峰.船舶主推進柴油機示功圖識別專家系統(tǒng)的研究[D].武漢理工大學碩士學位論文,2007(5):15
[4]何漢林,梅家斌.數(shù)值分析[M].北京:科學出版社,2007(2):130
[5]周龍寶,高宗英.內(nèi)燃機學[M].北京:機械工業(yè)出版社,2009(2):16
[6]林瑞霖,孫云嶺.基于瞬時轉(zhuǎn)速相關(guān)分析的柴油機工作均勻性檢驗[J].內(nèi)燃機,2006,8(4):19
[7]曹國豪.船舶柴油機示功圖測量系統(tǒng)研制[D].武漢理工大學碩士學位論文,2009(5):6
[8]孫怡,潛偉建,夏志瀾,等.在柴油機放熱曲線計算中兩個關(guān)鍵問題的解決方法[J].海軍工程大學學報,2001,6(13):86
[9]葉劍.135柴油機燃燒過程數(shù)值模擬[D].大連理工大學碩士學位論文,2008(12):32
[10]Ryuichi Kimura,Noboru Nakai,Hiroshi Mizutani.Diagnostic System of Diesel Engine by Statistical Vibration Analysis[R].Bulletin of the M.E.S.J,1966,24(1):14~20