何烈永 鄭曉豐 馬國海 王樹勇 周大權(quán) 張剛健

摘要：為研究如何提高內(nèi)燃機的熱效率，深入研究內(nèi)燃機缸內(nèi)燃燒過程，燃燒分析系統(tǒng)的設(shè)計和試驗是關(guān)鍵?；贜I CRIO硬件平臺、LabVIEW虛擬儀器軟件開發(fā)的燃燒分析系統(tǒng)可以實現(xiàn)汽油、柴油發(fā)動機的p-V圖、放熱率、缸內(nèi)溫度、IMEP、COV、爆震因子等參數(shù)的計算。開發(fā)出的系統(tǒng)在汽油發(fā)動機臺架上試驗，對比結(jié)果表明參數(shù)計算算法合理，精度符合要求，界面直觀操作方便。

Abstract： In order to study how to improve the thermal efficiency of internal combustion engine and the combustion process in cylinder， the design and test of combustion analysis system is the key.Based on NI CRIO hardware platform and LabVIEW virtual instrument software development， the combustion analysis system can realize the calculation of p-V diagram， heat release rate， cylinder temperature， IMEP， COV， knock factor and other parameters of gasoline and diesel engines.The system is tested on the gasoline engine bench， and the comparison results show that the parameter calculation algorithm is reasonable， the precision meets the requirements， and the interface is intuitive and easy to operate.

關(guān)鍵詞：汽油、柴油發(fā)動機;燃燒分析;測控系統(tǒng);算法;軟件實現(xiàn)與界面設(shè)計

Key words： gasoline and diesel engines;combustion analysis;measurement and control system;algorithm;software implementation and interface design

中圖分類號：TK435? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標識碼：A? ? ? ? ? ? ? ? ? ? ? ? ?; ? ? ? 文章編號：1674-957X（2021）12-0008-03

0? 引言

隨著國家節(jié)能減排法規(guī)日益加嚴，傳統(tǒng)內(nèi)燃機領(lǐng)域希望提高內(nèi)燃機的熱效率，大幅度實現(xiàn)節(jié)能減排效果?；诖耍瑖鴥?nèi)各大高校、企業(yè)對燃燒分析儀進行自主研發(fā)。由于國外燃燒分析儀價格昂貴，國內(nèi)燃燒分析儀的規(guī)模化應(yīng)用、用戶體驗有待提高，同時計算機、傳感器、信號處理等技術(shù)與時俱進，基于此，文章開發(fā)了一款內(nèi)燃機燃燒分析系統(tǒng)，目前正處于定型生產(chǎn)階段，預(yù)計于2021年年初上市。

1? 項目目標確定及系統(tǒng)組成

本項目的目標是設(shè)計和開發(fā)一款內(nèi)燃機燃燒分析系統(tǒng)，滿足國內(nèi)化生產(chǎn)需求，達到技術(shù)協(xié)議要求的精度、性能的工程目標。該系統(tǒng)主要技術(shù)特征包括：汽油、柴油發(fā)動機單個和每個試驗循環(huán)的p-V圖、p-Φ圖、瞬時燃燒放熱率、累計放熱量、缸內(nèi)溫度、IMEP、COV、最大缸壓、缸壓升高率、爆震信號處理等參數(shù)的計算和圖形化顯示。

該款燃燒分析系統(tǒng)的算法和軟件實現(xiàn)適用于柴油、汽油發(fā)動機臺架試驗。

2? 燃燒分析系統(tǒng)算法

燃燒分析軟件主要由缸壓采集、實時顯示、分析后處理、界面顯示組成。

燃燒分析軟件NI LabVIEW是基于數(shù)據(jù)流驅(qū)動的開發(fā)環(huán)境，可以兼顧其自身強大的數(shù)據(jù)流應(yīng)用能力、虛擬儀器面板軟件開發(fā)方便以及Matlab強大的數(shù)值計算、算法開發(fā)能力等優(yōu)點。且避免了LabVIEW在算法上的復(fù)雜操作和Matlab在工業(yè)采集數(shù)據(jù)流方面的不方便。

2.1 燃燒參數(shù)計算理論

從采集的p-Φ圖計算燃燒放熱率的出發(fā)點是熱力學第一定律。燃料燃燒放出的熱量一部分傳給工質(zhì)，用于增加工質(zhì)的內(nèi)能并對外做功，一部分通過燃燒室壁面散失到冷卻水中。

即氣缸內(nèi)的能量守恒方程進行計算：

式中，QB、Q、Qw分別代表燃料放熱量、？駐U工質(zhì)吸熱量和散熱量，W為工質(zhì)內(nèi)能變化，為工質(zhì)對活塞所作的機械功。

對曲軸轉(zhuǎn)角？漬取微分，得：

在數(shù)值計算中常用差商代替微商，在一個工作循環(huán)的步長曲軸轉(zhuǎn)角？駐？漬內(nèi)：

因此根據(jù)實測的p-Φ圖示功圖及有關(guān)參數(shù)，分別計算出？駐U、？駐W、？駐QW，從而計算？駐？漬步長內(nèi)的燃燒放熱量？駐QB。

為將公式（3）計算數(shù)值化，采用如下方法：

①單換能器方法THE SINGLE TRANSDUCER METHOD[4]——假設(shè)在周期開始關(guān)閉閥門后，不會有質(zhì)量通過系統(tǒng)邊界傳遞。凈放熱速率（包括做功和熱能增加）可表示為如下公式：

公式（4）用于計算單換能器法的熱釋放率Qn，使用測量量P和V及其導數(shù)。其中比熱容公式為[5]：

②單位曲軸轉(zhuǎn)角的工質(zhì)與燃燒室周壁面的換熱量：

其中：n：發(fā)動機轉(zhuǎn)速;A：換熱面積;T：氣缸內(nèi)工質(zhì)瞬時溫度;Tw：壁面的平均溫度;ag：瞬時平均換熱系數(shù)。

2.2 爆震因子

爆震因子是爆震識別的重要方法，Siemens公司設(shè)計的VDO算法因為計算準確可靠而得到廣泛工程應(yīng)用[6]。采用VDO算法計算爆震因子的公式如下：

本文在VDO算法基礎(chǔ)上進行改進應(yīng)用，主要改進措施是采用適配的濾波算法和閾值濾波等方法。當爆震因子大于3（設(shè)定的門檻值）則判定為爆震事件，當爆震事件個數(shù)占總循環(huán)的比值來衡量爆震強度KI。當KI<1%，判定為無爆震;當1%<KI<5%，判定為輕微爆震;當5%<KI<10%，判定為輕爆震;當KI>10%，判定為強爆震。

2.3 IMEP COV

指示平均有效壓力IMEP是對p-Φ圖求積分，反應(yīng)發(fā)動機的燃燒穩(wěn)定性，即Pi的變化系數(shù)，計算公式如下：

其中STD為標準偏差，單位[bar];MEAN為平均值，單位[bar]。

軟件自動計算出每個循環(huán)的IMEP值并顯示曲線，IMEP平均值，樣本標準偏差，得到COV_IMEP。當COV_IMEP小于5%說明發(fā)動機燃燒穩(wěn)定，當COV_IMEP大于5%則燃燒不穩(wěn)定。

3? 數(shù)據(jù)對比

該試驗發(fā)動機為一臺兩缸4沖程排量750cc的汽油發(fā)動機。通過設(shè)置發(fā)動機參數(shù)作為燃燒分析軟件的界面輸入。將自主研發(fā)的燃燒分析系統(tǒng)與奇石樂KiBox和上位機KiBoxCockpit軟件的數(shù)據(jù)進行對比。

3.1 最大爆壓及對應(yīng)角度

每個循環(huán)單缸最大爆壓如圖1所示。最大爆壓及對應(yīng)角度的平均值對比結(jié)果如表1所示。

3.2 缸內(nèi)溫度

發(fā)動機轉(zhuǎn)速為4000rpm下計算得到瞬時缸內(nèi)溫度對比數(shù)據(jù)如圖2所示，每個循環(huán)下最大缸溫如圖3所示。

3.3 爆震因子

通過計算爆震因子以及信號峰值等參數(shù)，每個循環(huán)下爆震因子如圖4所示。

4? 結(jié)論

文章對內(nèi)燃機燃燒分析系統(tǒng)進行仿真計算、軟硬件設(shè)計和相關(guān)試驗開發(fā)驗證。并在汽油發(fā)動機臺架上將自主研發(fā)的燃燒分析系統(tǒng)與奇石樂的Kibox進行對比試驗，結(jié)果表明自研系統(tǒng)新增了缸內(nèi)溫度等參數(shù)，其余參數(shù)計算算法合理，精度符合要求，界面直觀操作方便。

參考文獻：

[1]王兆文，郭凱，武文杰，等.基于實時性優(yōu)化的內(nèi)燃機燃燒分析系統(tǒng)研究[J].農(nóng)業(yè)機械學報，2020，51（01）：372-382.

[2]武文杰.柴油機燃燒放熱率實時分析系統(tǒng)的設(shè)計與開發(fā)[D].武漢：華中科技大學，2018.

[3]朱訪君，吳堅.內(nèi)燃機工作過程數(shù)值計算及其優(yōu)化[M].北京：國防工業(yè)出版社，1997.

[4]Grimm， B. M.， Johnson， R. T.， 1990， “Review of Simple Heat Release Computations”， SAE Technical Paper Series 900445.

[5]Hayes， T. K.， Savage， L. D.， Sorenson， S. C.， 1986， “Cylinder Pressure Data Acquisition and Heat Release Analysis on a Personal Computer”， SAE Technical Paper Series 860029.

[6]康健.基于動態(tài)窗口及統(tǒng)計方式的汽油機爆震識別方法[D].重慶：重慶大學，2010.

[7]陳樹學，劉萱.LabVIEW寶典（第2版）[M].北京：電子工業(yè)出版社，2017.