周鵬 李子競 佀慶濤 張建銳 張正興

（一汽解放汽車有限公司商用車開發(fā)院，長春 130011）

主題詞：發(fā)動機起停系統(tǒng) 直噴輔助起動 發(fā)動機動力學模型 能量損失 試驗標定

1 前言

在城市工況中，汽車經(jīng)常處于怠速狀態(tài)，增加油耗的同時帶來了排放問題。近幾年，發(fā)動機起停功能（Start-Stop）在乘用車上的應用正逐漸增多，其能夠有效降低整車污染物排放量和燃油消耗量[1]。汽油機直噴起動是Start-Stop停機后重新起動的方式之一，其是指在發(fā)動機短暫停機后，依靠膨脹缸火花塞點燃缸內(nèi)混合氣，推動活塞下行，以減輕起動機的負擔。

通過試驗可獲得成功起動所需的合理噴油點火時刻，并以此制定噴油策略。然而整個起動過程缸內(nèi)狀態(tài)復雜多變，絕大部分參數(shù)均為隨時間變化的瞬變量，因此，發(fā)動機試驗過程中難以記錄從靜態(tài)到動態(tài)基于曲軸轉(zhuǎn)角的缸壓、放熱率等參數(shù)。此外，點火過程較難控制，不易獲得期望的初始著火狀態(tài)，起動過程存在發(fā)動機轉(zhuǎn)速慢、氣缸漏氣及傳熱過程作用明顯等問題，因而根據(jù)試驗數(shù)據(jù)進行起動過程分析較為困難。

本文通過物理方法建立發(fā)動機動力學模型，計算等熵過程缸內(nèi)狀態(tài)，考慮溫度和時間影響下的傳熱量對起動過程的影響，計算不同起動位置下成功起動所需的最低燃燒爆發(fā)壓力或膨脹缸所需的最低機械能，以及影響起動過程的能量損失分配等參數(shù)，并對影響起動過程的關(guān)鍵因素進行了分析。

2 汽油機直噴起動系統(tǒng)建模

本文以某4缸汽油機為研究對象，基于表1所示的基本參數(shù)建立動力學模型[2]，為獲得試驗無法獲取的數(shù)據(jù)，對起動過程中的能量進行分解，將基于時間的數(shù)據(jù)轉(zhuǎn)化為基于曲軸轉(zhuǎn)角的數(shù)據(jù)，再根據(jù)物理原理進行建模，包括等熵膨脹壓縮模型、活塞及氣門系統(tǒng)摩擦阻力模型、漏氣模型、傳熱模型等。

表1 發(fā)動機基本參數(shù)

發(fā)動機起動過程中，分別存在1個壓縮缸和1個膨脹缸，另外2個氣缸處于氣門開啟狀態(tài)。建模過程集中在點火后的1個沖程內(nèi)，此時發(fā)動機阻力來自壓縮缸內(nèi)壓縮氣體和摩擦阻力，動力來自膨脹缸的燃料燃燒。模型方程組[3]為：

式中，Ti為指示扭矩；Pi為缸內(nèi)壓力；S為氣缸截面積；r為曲柄長；αi、βi分別為曲軸和連桿與活塞中線的夾角；i為氣缸編號；Tf為阻力扭矩；I為發(fā)動機轉(zhuǎn)動慣量；α為膨脹缸曲軸轉(zhuǎn)角。

2.1 發(fā)動機起動過程阻力模型

發(fā)動機在起動過程中需要克服氣缸壓縮阻力、泵氣損失和摩擦阻力矩。為建立可靠的阻力模型，需要對起動時阻力矩進行分解，通過發(fā)動機反拖摩擦功試驗，可以分解各部件以及泵氣帶來的阻力，低速平均摩擦有效壓力部分因難以測量，通過線性回歸獲得。

2.1.1 水溫和轉(zhuǎn)速因素

水溫和轉(zhuǎn)速均影響摩擦阻力。實際發(fā)動機起動后第1個沖程內(nèi)，轉(zhuǎn)速處于0～300 r/min，且水溫不確定，本文通過反拖試驗測定低轉(zhuǎn)速時發(fā)動機阻力矩與轉(zhuǎn)速和水溫的關(guān)系，如圖1所示，其可修正阻力模型中水溫和轉(zhuǎn)速對阻力計算的影響，也可據(jù)此將水溫作為直噴起動限制條件。

圖1 低轉(zhuǎn)速發(fā)動機阻力矩

2.1.2 泵氣損失

起動時2個缸處于換氣過程，為較易獲得期望的活塞位置，試驗在節(jié)氣門開度為5°～12°下進行。與該開度下節(jié)氣門泵氣損失相比，進、排氣門泵氣損失可以忽略，同時可以忽略起動時低排氣流量下渦輪和后處理的阻力。借鑒節(jié)氣門-進氣歧管-氣缸充氣系統(tǒng)平均值模型，進氣沖程結(jié)束時進氣歧管壓力近似等于缸內(nèi)壓力，發(fā)動機停機時經(jīng)歷了斷油后若干沖程的換氣，缸內(nèi)殘余廢氣可近似認為全部排出，因此有：

式中，Pim為進氣歧管壓力；Ra為空氣氣體狀態(tài)常數(shù)；Tim為進氣歧管溫度；Vim為進氣歧管容積；mth為通過節(jié)氣門的流量；map為進入氣缸的流量；Vd為發(fā)動機排量；ev為充氣效率；n為發(fā)動機轉(zhuǎn)速；Pamb為環(huán)境壓力；θ為節(jié)氣門開度。

根據(jù)Hendricks理論[4]計算缸內(nèi)充量，通過臺架標定獲得節(jié)氣門開度與節(jié)氣門前、后壓比及進氣流量的關(guān)系。曲軸箱壓力默認為0.1 MPa，結(jié)合缸內(nèi)壓力可計算任意曲軸轉(zhuǎn)角的泵氣阻力矩。

2.1.3 附件及曲軸阻力矩

勻速轉(zhuǎn)動扭矩扳手測量發(fā)動機轉(zhuǎn)速接近零時發(fā)動機曲軸扭矩，其隨曲軸轉(zhuǎn)角變化，然后通過拆下皮帶輪重復該過程，兩者扭矩相減獲得附件的阻力矩。如圖2所示，附件阻力矩最小值為8.00 N·m，最大值為11.60 N·m，隨曲軸變化規(guī)律不明顯，可取平均值9.60 N·m代替。

圖2 附件摩擦阻力測量

附件摩擦阻力矩由機油泵、水泵、燃油泵、發(fā)電機等運行時產(chǎn)生的阻力矩構(gòu)成。通過反拖試驗發(fā)現(xiàn)附件及曲軸摩擦阻力矩與發(fā)動機轉(zhuǎn)速呈線性關(guān)系。

2.1.4 活塞摩擦阻力矩

通過發(fā)動機反拖試驗，獲得活塞循環(huán)平均摩擦有效壓力與發(fā)動機轉(zhuǎn)速的關(guān)系，計算單缸活塞-氣缸壁摩擦阻力矩，再對4個缸的合力矩關(guān)于曲軸轉(zhuǎn)角積分求得活塞-氣缸壁平均有效摩擦功，計算結(jié)果與摩擦功試驗吻合[5-6]。

2.1.5 氣門機構(gòu)阻力矩

氣門機構(gòu)摩擦阻力矩存在于凸輪、挺桿、搖臂和氣門支點之間，通過反拖試驗和進排、氣門凸輪軸扭矩測量，氣門機構(gòu)的摩擦阻力矩不與發(fā)動機轉(zhuǎn)速成純單調(diào)關(guān)系[7]。

2.2 漏氣量計算

活塞環(huán)與氣缸壁、氣門與氣門座等部件之間并不完全密封，在發(fā)動機正常運轉(zhuǎn)過程中存在一定漏氣量，其在允許范圍內(nèi)不會對性能產(chǎn)生影響。在起動過程中，由于低轉(zhuǎn)速時間長，漏氣量從膨脹缸缸壓建立開始對輸出功有影響，故其不能被忽略。氣缸漏氣速度與缸壓呈函數(shù)關(guān)系[8-9]，可通過試驗獲得活塞靜止狀態(tài)下缸壓并擬合缸壓下降速率與時間關(guān)系獲得。由于測量過程壓縮缸沒有著火，近似認為其溫度變化不大，忽略傳熱對缸壓的影響。因此，漏氣壓降速率與壓力的關(guān)系為：

式中，p為缸壓；t為時間；f(t)為擬合關(guān)系；Δp為漏氣導致的壓降差分值。

2.3 最高爆發(fā)壓力

膨脹缸燃燒產(chǎn)生的壓力是起動過程唯一的動力來源。通過試驗發(fā)現(xiàn)，燃燒的最高壓力只維持在一定范圍，很難保持穩(wěn)定并進行控制，所以著火是一個概率事件。影響燃燒的因素較多，理想狀態(tài)下的最高爆發(fā)壓力可通過當量空燃比下的化學反應計算[10]：

式中，U為熱力學能；H為焓；V為燃燒時缸內(nèi)容積；R為混合氣氣體狀態(tài)常數(shù)，與組分有關(guān)；m為缸內(nèi)組分總質(zhì)量；T為缸內(nèi)燃燒溫度；R0為摩爾氣體狀態(tài)常數(shù)；c(CO2)、c(H2O)、c(N2)為各組分的體積分數(shù)；M(CO2)、M(H2O)、M(N2)為各組分摩爾質(zhì)量；pmax為最高爆發(fā)壓力；p1、p2分別為燃燒前、后缸內(nèi)壓力；n1、n2分別為燃燒前、后缸內(nèi)成分摩爾數(shù)；T1、T2分別為燃燒前、后缸內(nèi)溫度。

圖3所示為混合氣初始溫度與理論燃燒最高爆發(fā)壓力的關(guān)系，最高爆發(fā)壓力隨混合氣初始溫度升高而降低。

實際燃燒最高爆發(fā)壓力遠低于理論值，燃油非當量比燃燒，在靜止氣缸中噴射燃油，其蒸發(fā)與空氣混合質(zhì)量相差較大，因此空燃比一般為0.6～0.8，同時活塞環(huán)處漏氣也會對最高爆發(fā)壓力的建立產(chǎn)生影響。

圖3 最高爆發(fā)壓力

2.4 傳熱量計算

起動過程缸壓變化范圍較小，按定值近似標定傳熱系數(shù)[11]，傳熱面積根據(jù)曲軸轉(zhuǎn)角計算。發(fā)動機缸內(nèi)的溫度受傳熱影響按等容計算，忽略壓力對溫度的影響。圖4所示為計算缸壓與實際缸壓的對比結(jié)果，定值傳熱系數(shù)會對活塞運行前半段造成影響，較高的壓力溫度使傳熱增加[11]，計算壓力高于實際缸壓，后半段相反，但在本文工況下影響較小，可近似計算。

圖4 計算缸壓與實際缸壓對比

3 模型的輸入與求解

模型輸入為發(fā)動機基本參數(shù)，如缸徑、曲柄長、發(fā)動機起動初始位置，基于時間的壓縮缸與膨脹缸缸壓數(shù)據(jù)和發(fā)動機摩擦阻力矩通過試驗采集。

求解高階微分方程的解析解較為復雜，通過編寫程序，以0.2°～1°曲軸轉(zhuǎn)角為步長，用歐拉法將連續(xù)性方程離散化求數(shù)值解[12]。其結(jié)果可將扭矩隨機時間變化轉(zhuǎn)為隨曲軸轉(zhuǎn)角變化，如圖5、圖6所示，使分析基于曲軸轉(zhuǎn)角起動過程的影響因素成為可能。之后建立一種基于曲軸轉(zhuǎn)角且具備預測能力的物理模型，進而研究最優(yōu)控制方法。

4 試驗驗證

試驗臺架如圖7所示，選取29個工況點進行試驗，停機位置72°～132°、空燃比0.5～0.9、水溫80 ℃，而水溫作為摩擦阻力矩的重要影響因素單獨進行了摩擦測量試驗和不同水溫下的起動計算。

圖5 基于時間的扭矩

圖6 基于曲軸轉(zhuǎn)角的扭矩

圖7 試驗臺架

選取其中10個起始位置作為試驗計算分析點，缸壓和起始位置由試驗獲得，通過對比活塞最終停止位置是否一致可以驗證模型的準確性。觀察缸壓數(shù)據(jù)，挑選沒有發(fā)生反轉(zhuǎn)的工況點進行試驗，并與計算停機位置進行對比，如表2所示，壓縮上止點0表示活塞初始位置和停止位置。

表2 試驗與計算活塞停止時曲軸轉(zhuǎn)角對比 °CA

為確保模型精度，對每個工況下的計算缸壓與燃燒分析儀采集的缸壓進行對比，如圖8所示，結(jié)果表明，模型能較好地反映活塞的真實運動情況，可用于計算直噴起動過程。

圖8 缸內(nèi)壓力隨時間變化對比

5 計算結(jié)果分析

5.1 摩擦損失對快速起動的影響分析

在最高爆發(fā)壓力和起始曲軸轉(zhuǎn)角（上止點后90°）相同的情況下，分別對水溫為54℃和80℃的2種工況進行對比，測得54℃水溫下平均摩擦阻力近似為80℃水溫下的0.85倍，前者通過下止點成功起動，后者到達165°后反彈，試驗過程如圖9所示，在建立最高爆發(fā)壓力后活塞運動的前半程存在摩擦阻力矩與燃氣燃燒產(chǎn)生的動力矩，而隨著活塞下行，缸壓迅速降低，使得摩擦阻力的作用突顯。

圖9 不同摩擦阻力矩的影響

5.2 最高爆發(fā)壓力的影響分析

起動過程中能量耗散途徑較多，但最高爆發(fā)壓力是唯一動力源，對成功起動至關(guān)重要。為求得使發(fā)動機順利起動的最高爆發(fā)壓力，從試驗數(shù)據(jù)中挑選初始位置相同（上止點后90°）且最高爆發(fā)壓力為0.4 MPa、0.5 MPa、0.62 MPa和0.67 MPa的工況點，過程如圖10所示，結(jié)果表明，只有最高爆發(fā)壓力為0.67 MPa的工況中活塞通過了下止點，通過計算得到在上止點后90°成功起動所需要的最高爆發(fā)壓力極限值為0.67 MPa。

5.3 膨脹缸排氣門開啟時刻影響分析

膨脹缸排氣門開啟使缸壓迅速降至環(huán)境壓力，為分析膨脹缸排氣門開啟時刻對活塞運動距離的影響，分別設(shè)置5種不同膨脹缸排氣門正時，其初始爆發(fā)壓力和摩擦阻力矩相同，計算曲軸角速度，結(jié)果如圖11所示，排氣門開啟時刻越早，膨脹缸能量損失越大，活塞運動距離越短，但通常排氣門開啟時刻不會過早，所以相較于摩擦損失和最高爆發(fā)壓力，排氣正時對起動過程的影響較小。

圖10 最高爆發(fā)壓力的影響

圖11 膨脹缸排氣門開啟時刻的影響

5.4 轉(zhuǎn)動慣量影響分析

轉(zhuǎn)動慣量對發(fā)動機起動時間影響較大，為考察其對活塞運動距離的影響，分別設(shè)置5種轉(zhuǎn)動慣量，其都具有相同初始爆發(fā)壓力和摩擦阻力矩，不考慮漏氣和傳熱，計算曲軸角速度，結(jié)果如圖12所示，5種工況下的活塞均先加速后減速，在相近位置到達最低值。根據(jù)能量守恒定律，燃料的化學能被各種損失消耗，導致運動距離相近，若考慮傳熱和漏氣，活塞運動速度和距離會隨著其增加而減小。

圖12 轉(zhuǎn)動慣量的影響

5.5 活塞起動初始位置影響分析

本文的直噴起動策略是由初始活塞位置計算缸內(nèi)空氣量，按照固定空燃比計算噴油量。在忽略著火瞬間油氣混合質(zhì)量、傳熱、漏氣等因素的情況下，假設(shè)相同空燃比具有相同爆發(fā)壓力，活塞越接近下止點起動，有效膨脹比越低，在到達下止點前更容易保持較高缸壓。然而曲柄連桿的結(jié)構(gòu)決定了缸壓轉(zhuǎn)化扭矩的關(guān)系，扭矩轉(zhuǎn)化最易點可以通過計算獲得。而到達下止點的損失和與下止點的距離成正比。但是由于排氣門的提前開啟，膨脹缸過于接近下止點會使有效功較低。

為尋找最佳活塞初始位置，分別設(shè)置9種活塞起始位置，且具有相同最高爆發(fā)壓力，取消排氣門開啟過程，計算曲軸角速度，其結(jié)果如圖13所示，初始位置在100°～130°范圍內(nèi)壓縮缸較易越過上止點。

圖13 曲軸起動初始位置影響

5.6 起動過程能量分配

為考察活塞不同起始位置時能量的分配情況，選取起動位置72°～132°共10個工況，分別計算起動過程中的能量分配。計算結(jié)果如圖14所示，膨脹缸輸出功是扭矩與曲軸轉(zhuǎn)角的積分結(jié)果，缸壓扭矩轉(zhuǎn)換關(guān)系使膨脹缸輸出在曲軸轉(zhuǎn)角72°～98°范圍內(nèi)較為顯著。72°～132°范圍內(nèi)摩擦損失比例逐漸增加，壓縮缸消耗比例逐漸減少。

圖14 不同曲軸轉(zhuǎn)角起動過程能量分配

選取單一典型工況點作為模型輸入，活塞初始位置為壓縮上止點后90°，最高爆發(fā)壓力0.62 MPa，計算結(jié)果如圖15所示，對直噴輔助起停過程中扭矩損失進行分解，對起動過程的阻礙作用從高到低依次為摩擦損失、壓縮缸壓縮氣體損失、傳熱和漏氣損失，泵氣損失忽略不計，基于此可進行降摩擦或控制方法的改進。

圖15 扭矩損失詳細分解

6 結(jié)論

a.直噴輔助起動過程中動力損耗部分來自摩擦，其中氣門凸輪機構(gòu)和活塞環(huán)、活塞裙部與氣缸壁之間的摩擦占比最大，其次是附件、曲軸和泵氣損失。

b.最高爆發(fā)壓力是起動過程唯一的動力來源，對成功起動起關(guān)鍵作用。

c.排氣門開啟時刻越提前，膨脹缸能量損失越大，但排氣門開啟時刻提前量不會過大，氣門正時對起動過程的影響程度與摩擦和最高爆發(fā)壓力相比較小。

d.膨脹缸起動最佳初始位置的計算需綜合膨脹比、缸壓扭矩轉(zhuǎn)換關(guān)系、摩擦阻力以及氣門正時等因素。