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

(昆明理工大學(xué)云南省內(nèi)燃機(jī)重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650500)

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

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

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

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

1 不均勻度信號(hào)的量化與處理

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

圖1 不均勻度信號(hào)量化流程

圖2 各缸不均勻度信號(hào)

2 各缸均勻性控制策略

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

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

2.1 信號(hào)峰值判斷

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

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

圖4 信號(hào)峰值判斷模塊流程圖

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

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

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

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

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

當(dāng)工況點(diǎn)處于0x12、0x10區(qū)時(shí)，期望噴油量的系數(shù)計(jì)算公式分別為

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

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

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

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

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

在全部不均勻度信號(hào)計(jì)算完成后進(jìn)行修正，修正噴油量QCUC為

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

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

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

圖7 CUC修正油量的計(jì)算

2.4 修正噴油量分配

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

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

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

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

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

圖10 主噴修正噴油量計(jì)算

3 控制策略驗(yàn)證及分析

3.1 修正噴油量分配功能仿真測(cè)試

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

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

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

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

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

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

3.2 臺(tái)架驗(yàn)證及分析

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

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

圖13 發(fā)動(dòng)機(jī)試驗(yàn)臺(tái)架

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

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

圖14 各工況點(diǎn)修正后的時(shí)段信號(hào)

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

4 結(jié)束語

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

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

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