黃粉蓮，王鵬，申立中，王旭

(昆明理工大學(xué)云南省內(nèi)燃機重點實驗室，云南 昆明 650500)

發(fā)動機工作過程中，各氣缸內(nèi)的循環(huán)噴油量、燃燒特性和器件磨損狀況等存在差異，導(dǎo)致各氣缸輸出的扭矩不一致。各缸輸出扭矩不平衡會導(dǎo)致發(fā)動機曲軸扭振加劇、NVH特性劣化，嚴重影響柴油機的工作穩(wěn)定性和可靠耐久性能[1]。柴油機各缸均勻性控制(Cylinder Uniformity Control，簡稱“CUC”)基于發(fā)動機瞬時轉(zhuǎn)速波動情況計算各缸修正噴油量，從而減小各缸扭矩輸出的差異性。隨著柴油機高壓共軌噴射技術(shù)的發(fā)展，柴油機噴油量的精確計算和準確控制得以實現(xiàn)。各缸均勻性控制是確保柴油機穩(wěn)定運行的基本要求，是改善循環(huán)內(nèi)各缸燃燒的均勻性、提高柴油機綜合性能的關(guān)鍵技術(shù)。

針對各缸均勻性的研究，國內(nèi)外專家學(xué)者從多角度對該問題進行了研究分析。Van Nieuwstadt等[2]開發(fā)了一種針對因制造偏差和老化情況不同而導(dǎo)致各氣缸不均勻性的方案，提出了幾種用于平衡各缸扭矩輸出的傳感器應(yīng)用方法和算法。Cavina等[3]提出一種基于發(fā)動機瞬時轉(zhuǎn)速來分析和計算發(fā)動機失火及工作不均勻的方法，根據(jù)扭矩與速度的關(guān)系，在考慮扭振因素的情況下，構(gòu)造一個失火檢測參數(shù)。研究表明：根據(jù)發(fā)動機速度波動來診斷扭矩輸出不均勻性和失火具有較高的可行性。Schiefer等[4]通過缸內(nèi)壓力傳感器測量氣缸內(nèi)的燃燒壓力來表征各缸扭矩的輸出情況。研究表明：在所有量化各缸扭矩輸出不均勻的方法中，直接測量各氣缸的燃燒壓力是最準確的，其缺點是在氣缸蓋上打孔安裝缸內(nèi)壓力傳感器會干擾氣缸內(nèi)的氣體運動，降低氣缸蓋的強度，增加制造成本。Macián[5]提出了一種基于線性模型的不均勻性控制算法，該算法是將傅里葉變換的低頻區(qū)和噴油量不均勻聯(lián)系起來，通過模型反算得到各缸噴油量差異。Solmaz[6]采用一種數(shù)學(xué)模型來研究失火和循環(huán)變化對內(nèi)燃機曲軸轉(zhuǎn)速波動的影響，模型中考慮了缸內(nèi)氣體壓力、發(fā)動機摩擦、產(chǎn)生線性和圓周運動部件的慣性效應(yīng)對模型的影響，結(jié)果表明，單缸或多缸失火都會導(dǎo)致曲軸轉(zhuǎn)速的波動。Junichi[7]將相同曲軸轉(zhuǎn)角內(nèi)產(chǎn)生的平均扭矩作為修正依據(jù)，設(shè)計了一套發(fā)動機扭矩平衡算法。

國內(nèi)的研究工作主要以瞬時轉(zhuǎn)速相鄰波峰(或波谷)之間的差值，或者對各齒的瞬時轉(zhuǎn)速進行統(tǒng)計分析來量化不均勻度。胡重慶等[8]提出了一種基于瞬時轉(zhuǎn)速統(tǒng)計分析的失火檢測方法，通過調(diào)整、濾波、切分與重采樣后，將瞬時轉(zhuǎn)速序列的處理問題轉(zhuǎn)化為多變量的分析問題，采用主分量分析方法獲得各缸輸出扭矩的不均勻度。周泉[9]提出了一種基于帶通濾波的各缸工作不均勻性檢測與控制方法，當各缸輸出扭矩出現(xiàn)不均勻時，頻率為最低Z/2(Z為發(fā)動機的氣缸數(shù))的波的振幅將加大。李建秋等[10]利用瞬時轉(zhuǎn)速波動分量的波峰差值作為反映發(fā)動機工作不均勻的指標，建立了發(fā)動機各缸均勻性反饋控制算法。王永庭等[11]通過比較各氣缸做功沖程波峰和波谷之間的轉(zhuǎn)速差來判斷各氣缸工作是否均勻，并以此為基礎(chǔ)進行了柴油機各缸噴油量不均勻調(diào)節(jié)的硬件在環(huán)仿真研究。結(jié)果表明，可通過改變單缸噴油量改善柴油機的工作不均勻性。劉世元等[12-13]提出了采用轉(zhuǎn)速波動信號在線監(jiān)測發(fā)動機各缸工作不均勻性和失火故障的方法，提出了多種量化各缸扭矩輸出不均勻的方法。

綜上所述，各缸均勻控制的關(guān)鍵是通過檢測不均勻度信號對發(fā)動機各缸進行均勻度控制，根據(jù)不均勻度信息計算各缸各次噴射的修正噴油量，達到消除各缸輸出扭矩不均勻的目的。本研究基于柴油機燃油噴射控制機制，設(shè)計了多缸柴油機各缸均勻性控制算法，采用Matlab/Simulink軟件建立了圖形化的各缸均勻性控制模型，通過模型在環(huán)測試和臺架試驗對控制策略功能進行驗證。

1 不均勻度信號的量化與處理

為了能夠得到不均勻度信號，需要對時段信號進行預(yù)處理和不均勻度信號提取。不均勻度信號的量化與處理流程如圖1所示。時段信號的預(yù)處理包括缺齒處理、重采樣和放大處理，不均勻信號的提取包括一倍基頻和二倍基頻下的帶通濾波。缺齒處理通過判斷當前信號與上一信號的比例，來確定當前信號是否為缺齒信號，用于去除缺齒信號對不均勻度量化的干擾。重采樣處理的作用是在保證信號準確的前提下減小數(shù)據(jù)量，重采樣因子越大，數(shù)據(jù)量越小，信號越不準確。設(shè)置重采樣因子為1～5進行仿真，當重采樣因子為3時，需要處理的數(shù)據(jù)量最少，計算結(jié)果的準確性較好，因此選取最佳重采樣因子為3。放大處理的功能是將時段信號放大，保證帶通濾波后的結(jié)果不會失真。通過帶通濾波濾出所需頻率的信號，對兩個穩(wěn)定后的濾波結(jié)果進行疊加，得到不均勻度信號(見圖2)。

圖1 不均勻度信號量化流程

圖2 各缸不均勻度信號

2 各缸均勻性控制策略

各缸均勻性控制策略是根據(jù)量化與處理后的不均勻度信息，計算各缸各次噴射的修正噴油量，達到消除各缸輸出扭矩不均勻的目的。各缸均勻性控制策略主要包括信號峰值判斷、控制區(qū)域判斷及系數(shù)計算、噴油修正與限制、修正噴油量分配等功能模塊，軟件架構(gòu)見圖3。

圖3 各缸均勻性控制軟件架構(gòu)

2.1 信號峰值判斷

從不均勻度量化模塊輸出的不均勻度信號可能出現(xiàn)過大或過小的情況。當不均勻度信號過小時，可能是由于采樣率過低或當前各缸扭矩輸出較為均勻；而當不均勻度信號過大時，可能是發(fā)動機出現(xiàn)其他故障導(dǎo)致的，例如噴油器損壞。修正噴油量計算及限制模塊不能采用這兩種信號作為計算修正噴油量的依據(jù)。因此，需要對不均勻度信號進行評估，并控制不均勻度信號是否進入修正噴油量計算及限制模塊，流程圖見圖4。

若不均勻度過小，則不啟用控制；若不均勻度過大，則不啟用控制，并上報故障信息；若不均勻度處于合適范圍內(nèi)，啟用各缸均勻性控制。上閾值和下閾值是通過兩個標定量確定的，可以通過臺架標定改變控制策略的啟用范圍。

圖4 信號峰值判斷模塊流程圖

2.2 控制區(qū)域判斷及系數(shù)計算

為了保證控制效果，對于不同的發(fā)動機工況，各缸均勻性控制策略采用不同的控制方法。將各缸均勻性控制的運行區(qū)間分為閉環(huán)控制區(qū)域、開環(huán)控制區(qū)域和非控制區(qū)域(見圖5)。 0x11為閉環(huán)控制區(qū)域；0x31和0x21是僅受發(fā)動機平均轉(zhuǎn)速影響的開環(huán)控制區(qū)；0x10和0x12是僅受期望油量影響的開環(huán)控制區(qū)；0x30、0x32、0x20、0x23是受發(fā)動機平均轉(zhuǎn)速和期望噴油量同時影響的開環(huán)控制區(qū)；其余部分為非控制區(qū)。

圖5 各缸均勻性控制策略工作控制劃分

當工況處于開環(huán)控制區(qū)時，計算可能出現(xiàn)兩種情況：一是平均轉(zhuǎn)速和期望噴油量都處在開環(huán)控制區(qū)，此時系數(shù)為轉(zhuǎn)速系數(shù)cn乘以油量系數(shù)cq；二是兩者只有一個處在開環(huán)控制區(qū)，另一個處在閉環(huán)控制區(qū)內(nèi)，此時系數(shù)為轉(zhuǎn)速系數(shù)cn或油量系數(shù)cq。如果工況點處在閉環(huán)控制區(qū)，兩者系數(shù)都為1。

當工況點處于0x31、0x21區(qū)時，平均轉(zhuǎn)速的系數(shù)計算公式分別為

當工況點處于0x12、0x10區(qū)時，期望噴油量的系數(shù)計算公式分別為

式中：cn為轉(zhuǎn)速系數(shù)；cq為油量系數(shù)；nGMin為CUC管理轉(zhuǎn)速最小值；nGMax為CUC管理轉(zhuǎn)速最大值；nCMin為CUC控制轉(zhuǎn)速最小值；nCMax為CUC控制轉(zhuǎn)速最大值；qGMin為CUC管理油量最小值；qGMax為CUC管理油量最大值；qCMin為CUC控制油量最小值；qCMax為CUC控制油量最大值。

2.3 噴油量修正與限制

修正噴油量計算與限制模塊的功能是根據(jù)不均勻度信號計算修正噴油量，并對計算結(jié)果進行限制。采用PI算法對修正噴油量計算進行閉環(huán)控制，使得各缸不均勻度信號的最低Z/2個階次的合成波幅值為0。假設(shè)在每個不均勻信號進入對應(yīng)氣缸的PI積分器后就進行修正，每次修正值都是基于第一個進入PI積分器的不均勻信號。在第k+1個信號進入PI控制器時，當前偏差為

e(k)=[T(k)-T(0)]。

式中：T(k)為進入PI積分器的信號；T(0)為第一個進入PI積分器的信號?；诘谝粋€不均勻度信號的修正噴油量u(k)為

由于假設(shè)每個不均勻度信號進入PI積分器后就進行修正，所以基于k-1個不均勻度信號的修正噴油量Δu(k)為

Δu(k)=u(k)-u(k-1)。

在全部不均勻度信號計算完成后進行修正，修正噴油量QCUC為

為了防止過多或過少的修正噴油量添加到實際噴射油量中，控制策略要對PI積分器的計算結(jié)果進行限制。修正油量監(jiān)控與限制曲線見圖6，該曲線通過發(fā)動機臺架試驗標定得到。在不同控制區(qū)域，修正噴油量限制曲線的斜率是不同的。針對不同的期望噴油量，修正噴油量的限制值也不同。為了對應(yīng)不同控制區(qū)域的監(jiān)控需求，將修正噴油量的范圍分為A、B、C、D 4個區(qū)域。D區(qū)為監(jiān)控不啟用區(qū)域，期望噴油量在該區(qū)域內(nèi)時，監(jiān)控功能不啟用。C區(qū)為修正噴油量正常區(qū)域，在該區(qū)域內(nèi)，修正噴油量為合理值，監(jiān)控功能不向診斷系統(tǒng)管理(Diagnostic System Management，DSM)報告錯誤。當修正噴油量出現(xiàn)在A區(qū)或B區(qū)時，則認為發(fā)動機出現(xiàn)故障，向DSM報告故障。

圖6 修正油量監(jiān)控與限制曲線

修正噴油量的最大值是由修正噴油量限制曲線決定的，輸入變量為期望噴油量?？刂七壿嬕妶D7。首先，在I環(huán)節(jié)計算出修正噴油量后，與修正噴油量限制曲線的輸出值進行比較，取二者最小值。之后，修正噴油量限制曲線的輸出值與PI積分器的計算結(jié)果進行比較，取二者最小值作為當前計算的修正噴油量。

圖7 CUC修正油量的計算

2.4 修正噴油量分配

完成修正噴油量的計算后，需要對修正噴油量進行分配，保證每一缸修正后的各次噴射符合修正前的比例。修正噴油量分配計算需要保證修正噴油量不會干擾正常的噴射。如果修正噴油量大于總期望噴油量，則認為當前修正噴油量不可信，修正噴油量不能被添加到各次射中。修正噴油量分配模塊邏輯圖見圖8。

模型中設(shè)計了冗余程序，防止修正噴油量計算錯誤。若修正噴油量大于等于期望噴油量，則認為當前修正噴油計算存在錯誤，不進行各缸均勻性控制；若修正噴油量小于期望噴油量，則認為當前修正噴油量正確，可以將修正噴油量分配到各次噴射中。本研究中試驗機型的噴射結(jié)構(gòu)包含6次噴射(預(yù)噴3、預(yù)噴2、預(yù)噴1、主噴、后噴2和后噴1)。預(yù)噴3、2、1的修正油量與后噴2、1的修正油量計算相同，修正油量的計算邏輯見圖9(以預(yù)噴3為例)。

圖8 修正噴油量分配模塊邏輯圖

圖9 預(yù)噴、后噴修正噴油量計算

主噴修正噴油量計算邏輯見圖10，將預(yù)噴和后噴的修正噴油量從總的修正噴油量中減去即可。在完成全部噴射的修正噴油量計算后，將計算結(jié)果保存在數(shù)組內(nèi)，以供其他功能模塊調(diào)用。

圖10 主噴修正噴油量計算

3 控制策略驗證及分析

3.1 修正噴油量分配功能仿真測試

針對閉環(huán)控制區(qū)和開環(huán)控制區(qū)對修正噴油量分配功能進行模型在環(huán)測試。期望噴油量由扭矩結(jié)構(gòu)與扭矩-油量轉(zhuǎn)換模塊計算得到。各次噴射的期望噴油量由噴油控制模塊計算得到。當發(fā)動機平均轉(zhuǎn)速為900 r/min、期望噴油量為10 mg/hub(毫克每沖程)時，控制策略處于閉環(huán)控制區(qū)，各次噴射的噴油量見表1設(shè)置。仿真得到各次噴射的修正噴油量(見圖11)，各缸的修正噴油量按比例分配到各次噴射中。

表1 閉環(huán)控制區(qū)各次噴射的期望噴油量

圖11 閉環(huán)控制區(qū)各次噴射的修正噴油量

當發(fā)動機平均轉(zhuǎn)速為750 r/min、期望噴油量為1.6 mg/hub時，控制策略處于開環(huán)控制區(qū)，各次噴射的噴油量設(shè)置見表2。計算得到各次噴射的修正噴油量(見圖12)。各缸的修正噴油量按比例分配到各次噴射中。

表2 開環(huán)控制區(qū)各次噴射的期望噴油量

圖12 開環(huán)控制區(qū)各次噴射的修正噴油量

3.2 臺架驗證及分析

臺架試驗采用D25電控高壓共軌柴油機，基本參數(shù)見表3，試驗臺架見圖13。將各缸均勻性控制策略模型與發(fā)動機其他控制系統(tǒng)模型進行集成，通過自動代碼功能與編譯軟件，與底層軟件進行集成與編譯，生成ECU可執(zhí)行文件(.elf、hex、s19)及標定測量文件(.a2l)，將可執(zhí)行文件通過專業(yè)的調(diào)試工具刷寫到ECU中即可實現(xiàn)對柴油機的有效控制，采用INCA標定軟件對CUC策略進行調(diào)節(jié)。

表3 D25柴油機基本參數(shù)

圖13 發(fā)動機試驗臺架

為了確?？刂撇呗栽陂]環(huán)控制區(qū)內(nèi)的任意工況下都可以進行修正，進行4個工況點的測試，分別為低轉(zhuǎn)速低負荷(900 r/min，10 mg/hub)、低轉(zhuǎn)速高負荷(900 r/min，50 mg/hub)、高轉(zhuǎn)速低負荷(1 200 r/min，15 mg/hub)和高轉(zhuǎn)速高負荷(1 200 r/min，50 mg/hub)。并設(shè)置第4缸的實際噴油量為其他氣缸的95%。修正噴油量計算結(jié)果見表4，各缸工況點修正后的時段信號見圖14。

表4 各工況下的修正噴油量 mg/hub

圖14 各工況點修正后的時段信號

通過分析圖14和表4可知，控制策略可以根據(jù)發(fā)動機運行工況實時計算各缸的修正噴油量，并按照各次噴射期望噴油量的比例將修正噴油量分配到各次噴射中。表4中各修正噴油量變化與圖13的變化趨勢吻合，經(jīng)過修正后的時段信號變得更加均勻，柴油機各缸扭矩輸出不均勻度減小，說明各缸均勻性控制策略設(shè)計合理、有效。

4 結(jié)束語

通過對時段信號預(yù)處理和不均勻度信號提取，獲得量化的不均勻度信息。最佳重采樣因子為3時，可以有效減小計算資源的占據(jù)，并保證計算結(jié)果的準確性。設(shè)計了判斷結(jié)果延遲功能，保證了時段信號完全被處理成不均勻度信號，并減小計算方法反復(fù)變化帶來的影響。

選取單缸扭矩輸出不均勻的情況分別對閉環(huán)控制區(qū)和開環(huán)控制區(qū)進行修正噴油量分配的模型在環(huán)測試，修正噴油量能夠按照各次噴射期望噴油量的比例被分配到各次噴射中。修正噴油量分配功能能夠保證修正噴油量不會改變各次噴射油量的比例。

針對發(fā)動機轉(zhuǎn)速為800～1 300 r/min、噴油量為2～50 mg/hub的工況范圍進行臺架試驗，結(jié)果表明，控制策略可以根據(jù)發(fā)動機運行工況實時計算各缸的修正噴油量，各缸扭矩輸出不均勻度減小，說明各缸均勻性控制策略設(shè)計合理、有效。