鄒永湘 倪計(jì)民 張 棟 石秀勇 陳 泓

（同濟(jì)大學(xué)）

1 前言

EGR技術(shù)是通過(guò)將部分廢氣引入進(jìn)氣系統(tǒng)以改變發(fā)動(dòng)機(jī)缸內(nèi)充量成分，從而降低發(fā)動(dòng)機(jī)NOx生成量，是現(xiàn)階段控制柴油機(jī)NOx排放的重要技術(shù)路線[1]。

本文使用仿真軟件BOOST建立某款增壓柴油機(jī)一維流動(dòng)模型，并應(yīng)用試驗(yàn)設(shè)計(jì)（DoE）技術(shù)，對(duì)文丘里管式高壓循環(huán)EGR系統(tǒng)進(jìn)行結(jié)構(gòu)優(yōu)化，提升EGR系統(tǒng)工作能力，使EGR率絕對(duì)值在全工況范圍內(nèi)相對(duì)于優(yōu)化前提高1%～3%。

2 研究樣機(jī)

研究樣機(jī)為某公司一款3L輕型載貨汽車的柴油機(jī)，為2氣門、4缸、四沖程、立式、直列水冷、電控直噴、渦輪增壓柴油機(jī)。原機(jī)主要技術(shù)參數(shù)見(jiàn)表1。

表1 原機(jī)主要技術(shù)參數(shù)

3 柴油機(jī)EGR系統(tǒng)模型的建立

3.1 一維流動(dòng)模型建立

根據(jù)研究樣機(jī)的基本結(jié)構(gòu)參數(shù)，在對(duì)樣機(jī)流動(dòng)過(guò)程進(jìn)行分析的基礎(chǔ)上，應(yīng)用AVL BOOST軟件建立柴油機(jī)原機(jī)一維仿真計(jì)算模型，如圖1所示。在建模過(guò)程中，燃燒模型選用Vibe模型，增壓器模型選用Full Model全參數(shù)模型，傳熱模型選用Woschni 1978模型。

分析認(rèn)為，燃燒模型和增壓器模型的選用對(duì)研究結(jié)果具有重要影響，因此需要對(duì)建模過(guò)程中燃燒模型和增壓器模型的應(yīng)用過(guò)程進(jìn)行詳細(xì)考慮。

3.1.1 Vibe燃燒模型的應(yīng)用

在BOOST軟件中，Vibe燃燒參數(shù)主要包括燃燒起始角、燃燒持續(xù)期和燃燒室形狀影響因子。通過(guò)BOOST軟件內(nèi)置的Burn功能模塊，以柴油機(jī)基本結(jié)構(gòu)參數(shù)和原機(jī)缸內(nèi)壓力試驗(yàn)數(shù)據(jù)為輸入條件，可以精確計(jì)算Vibe燃燒模型的燃燒參數(shù)，能夠滿足研究?jī)?nèi)容的仿真計(jì)算精度和要求。

3.1.2 增壓器Full Model全參數(shù)模型的應(yīng)用

在BOOST軟件中，增壓器全參數(shù)模型的應(yīng)用主要包括壓氣機(jī)Map圖和渦輪機(jī)Map圖的建立。壓氣機(jī)Map圖的參數(shù)包括增壓器轉(zhuǎn)速、壓氣機(jī)效率、壓氣機(jī)空氣流量和增壓比；渦輪機(jī)Map圖的參數(shù)包括膨脹比、渦輪機(jī)流量和渦輪機(jī)效率。在增壓器全參數(shù)模型建立過(guò)程中，根據(jù)增壓器的試驗(yàn)數(shù)據(jù)，分別對(duì)壓氣機(jī)和渦輪機(jī)輸入相應(yīng)的Map圖，以保證增壓器全參數(shù)模型應(yīng)用的準(zhǔn)確性，進(jìn)而滿足研究?jī)?nèi)容的計(jì)算精度要求。

3.2 模型標(biāo)定

原機(jī)模型標(biāo)定結(jié)果如圖2所示。從圖2可以看出，所建模型的計(jì)算結(jié)果曲線與試驗(yàn)結(jié)果曲線的趨勢(shì)吻合，各項(xiàng)誤差均在7%以內(nèi)，可滿足模擬計(jì)算的精度要求。

3.3 高壓循環(huán)EGR系統(tǒng)流動(dòng)模型

為了減少NOx排放以滿足國(guó)Ⅳ排放標(biāo)準(zhǔn)，對(duì)原發(fā)動(dòng)機(jī)進(jìn)行EGR技術(shù)改型。在原機(jī)模型試驗(yàn)驗(yàn)證基礎(chǔ)上，建立帶有高壓循環(huán)EGR系統(tǒng)的柴油機(jī)流動(dòng)模型，如圖3所示。

根據(jù)EGR控制策略，選取13工況法中除外特性和怠速工況以外的9個(gè)點(diǎn)作為計(jì)算工況點(diǎn)，即1 900 r/min（記為 A）、2 300 r/min（記為 B）、2 700 r/min（記為 C）轉(zhuǎn)速下 25%、50%、75%負(fù)荷工況，用 A25、A50、A75、B25、B50、B75、C25、C50、C75 來(lái)表示。

通過(guò)計(jì)算發(fā)現(xiàn)增壓柴油機(jī)采用高壓循環(huán)EGR系統(tǒng)時(shí)，在部分工況下EGR引出點(diǎn)氣體壓力低于引入點(diǎn)壓力，從而產(chǎn)生負(fù)EGR率，最大負(fù)EGR率可達(dá)-5.4%。引入點(diǎn)與引出點(diǎn)的壓差及EGR率如圖4所示。

為了解決該問(wèn)題，在EGR引入位置安裝文丘里管進(jìn)行增壓氣體的減壓，以保證必要的正EGR率[2]。

4 基于DoE的結(jié)構(gòu)布置優(yōu)化

為提高EGR系統(tǒng)的工作能力，并使EGR技術(shù)對(duì)發(fā)動(dòng)機(jī)動(dòng)力性和燃油經(jīng)濟(jì)性造成的不利影響最小化，采用DoE技術(shù)對(duì)文丘里管式高壓循環(huán)EGR系統(tǒng)結(jié)構(gòu)布置方案進(jìn)行多目標(biāo)優(yōu)化設(shè)計(jì)。DoE優(yōu)化設(shè)計(jì)流程如圖5所示。

4.1 因素定義

4.1.1 設(shè)計(jì)參數(shù)定義

文丘里管式高壓循環(huán)EGR系統(tǒng)的主要結(jié)構(gòu)參數(shù)有：

a. 文丘里管喉口直徑D；

b. 文丘里管漸縮段錐角S；

c. 文丘里管漸擴(kuò)段錐角K；

d. 文丘里管混合段長(zhǎng)度L；

e.EGR閥到引入點(diǎn)的距離E。

通過(guò)2水平析因設(shè)計(jì)[3]，利用優(yōu)化軟件SPSS對(duì)設(shè)計(jì)參數(shù)進(jìn)行敏感性分析?；陲@著性t檢驗(yàn)的結(jié)果，綜合分析上述結(jié)構(gòu)參數(shù)對(duì)新鮮空氣進(jìn)氣量和EGR率的影響效應(yīng)，最終選取D、L和E作為主要設(shè)計(jì)因素進(jìn)行優(yōu)化設(shè)計(jì)。設(shè)計(jì)參數(shù)定義如表2所列。

表2 設(shè)計(jì)參數(shù)定義 mm

4.1.2 優(yōu)化目標(biāo)定義

定義優(yōu)化目標(biāo)為發(fā)動(dòng)機(jī)9個(gè)工況下的新鮮空氣進(jìn)氣量和EGR率（即共18個(gè)目標(biāo)），使各工況下的新鮮空氣進(jìn)氣量和EGR率盡可能大，并以提高中高轉(zhuǎn)速及中高負(fù)荷的新鮮空氣進(jìn)氣量和EGR率為主要目的[4]。

4.2 試驗(yàn)方案設(shè)計(jì)

3個(gè)設(shè)計(jì)參數(shù)均選取5個(gè)水平 （最大，3/4，1/2，1/4，最?。┻M(jìn)行方案設(shè)計(jì)，若進(jìn)行全因子試驗(yàn)，則需進(jìn)行5×5×5=125次計(jì)算，為減少仿真次數(shù)，選取田口試驗(yàn)設(shè)計(jì)法[5]進(jìn)行試驗(yàn)設(shè)計(jì)。

4.3 模擬計(jì)算

按照設(shè)定的試驗(yàn)方案，在BOOST軟件中完成模擬計(jì)算，計(jì)算結(jié)果將作為后續(xù)數(shù)學(xué)建模工作的數(shù)據(jù)基礎(chǔ)。

4.4 建立數(shù)學(xué)模型

基于BOOST軟件的仿真試驗(yàn)數(shù)據(jù)，在優(yōu)化軟件OPTIMUS中建立設(shè)計(jì)參數(shù)D、L、E與優(yōu)化目標(biāo)各轉(zhuǎn)速下的新鮮空氣進(jìn)氣量FA、EGR率之間的數(shù)學(xué)模型。

在OPTIMUS中分別建立RBF模型、Kriging模型、泰勒多項(xiàng)式模型和用戶多項(xiàng)式模型，分析對(duì)比各模型的回歸系數(shù)R2[6]，發(fā)現(xiàn)泰勒多項(xiàng)式模型與計(jì)算數(shù)據(jù)擬合度最好，故采用泰勒多項(xiàng)式模型。

4.5 基于泰勒多項(xiàng)式模型的優(yōu)化

4.5.1 全局優(yōu)化

首先采用非支配排序遺傳算法（Non-dominated Sorting Genetic Algorithms，NSGA）進(jìn)行基于數(shù)學(xué)模型的全局尋優(yōu)，根據(jù)目標(biāo)權(quán)重選擇最優(yōu)的Pareto解：D=27.6 mm，L=43 mm，E=519 mm。

4.5.2 局部?jī)?yōu)化

局部尋優(yōu)算法是將全局尋優(yōu)算法的最優(yōu)值作為迭代計(jì)算起點(diǎn)值進(jìn)行尋優(yōu)設(shè)計(jì)，在起點(diǎn)值附近進(jìn)行更加精確的尋優(yōu)，并提升優(yōu)化工作效率，減小工作量[6]。在全局尋優(yōu)的基礎(chǔ)上，基于泰勒多項(xiàng)式模型，利用NSGA算法進(jìn)行局部尋優(yōu)。通過(guò)對(duì)全局尋優(yōu)Pareto解集分析，確定局部尋優(yōu)區(qū)間如表3所列。

表3 局部尋優(yōu)區(qū)間 mm

對(duì)局部尋優(yōu)的Pareto解集進(jìn)行分析，根據(jù)目標(biāo)權(quán)重，選取最優(yōu)的Pareto解：D=28.4 mm，L=46 mm，E=510 mm。

4.5.3 優(yōu)化結(jié)果分析

全局優(yōu)化和局部?jī)?yōu)化結(jié)果如圖6所示。

由優(yōu)化結(jié)果可知，基于泰勒多項(xiàng)式模型的NSGA全局尋優(yōu)對(duì)全工況內(nèi)的新鮮空氣進(jìn)氣量和EGR率都有提升，尤其是對(duì)中高轉(zhuǎn)速、中高負(fù)荷的EGR率提升最為明顯。

通過(guò)兩種尋優(yōu)算法的結(jié)果對(duì)比，可得出如下結(jié)論：

a. 局部尋優(yōu)算法與全局尋優(yōu)算法得出的新鮮空氣進(jìn)氣量與EGR率基本一致。

b. 局部尋優(yōu)算法由于定位更加準(zhǔn)確，因而得出的最優(yōu)值比全局尋優(yōu)算法的稍好。局部尋優(yōu)相對(duì)于全局尋優(yōu)，在大多數(shù)工況內(nèi)EGR率有小幅度的提升，提升率為0.1%～3.0%，低負(fù)荷工況的EGR率有小幅度的降低；在全工況內(nèi)，新鮮空氣流量有小幅度上升，上升幅度在1%以內(nèi)。

c. 相對(duì)于初始方案，局部尋優(yōu)的新鮮空氣進(jìn)氣量和EGR率均有較大程度的改善，新鮮空氣進(jìn)氣量上升1%～4%，EGR率絕對(duì)值上升1%～5%。

因此，以局部尋優(yōu)方案作為該高壓循環(huán)EGR系統(tǒng)最終的結(jié)構(gòu)設(shè)計(jì)方案。

4.6 仿真驗(yàn)證

將局部尋優(yōu)優(yōu)化結(jié)果代入BOOST軟件中和泰勒多項(xiàng)式模型中進(jìn)行計(jì)算，結(jié)果如圖7所示。

由圖7可知，數(shù)學(xué)模型所計(jì)算的優(yōu)化結(jié)果數(shù)據(jù)與發(fā)動(dòng)機(jī)流動(dòng)分析計(jì)算的數(shù)據(jù)誤差很小，各工況的新鮮空氣進(jìn)氣量和EGR率相對(duì)誤差值均在5%以內(nèi)，從而驗(yàn)證了基于數(shù)學(xué)模型優(yōu)化結(jié)果的正確性，進(jìn)而也證明了DoE優(yōu)化過(guò)程的正確性。

5 優(yōu)化方案對(duì)發(fā)動(dòng)機(jī)性能的影響

對(duì)初始方案下、未采用文丘里管EGR系統(tǒng)條件下、優(yōu)化方案（最優(yōu)解）下的發(fā)動(dòng)機(jī)流動(dòng)特性數(shù)據(jù)進(jìn)行對(duì)比，具體數(shù)據(jù)如圖8所示。

由圖8可知：

a. 未采用文丘里管的高壓循環(huán)EGR系統(tǒng)，在低轉(zhuǎn)速區(qū)域由于負(fù)EGR率的存在，使發(fā)動(dòng)機(jī)充氣效率大幅下降，發(fā)動(dòng)機(jī)動(dòng)力性和燃油經(jīng)濟(jì)性相對(duì)原機(jī)惡化嚴(yán)重；中高轉(zhuǎn)速時(shí)壓差較小，EGR率較低。

b. 初始方案喉口直徑過(guò)小，文丘里管內(nèi)流速接近當(dāng)?shù)芈曀俣l(fā)生氣體壅塞現(xiàn)象，進(jìn)氣阻力過(guò)大，影響發(fā)動(dòng)機(jī)新鮮空氣進(jìn)氣量。采用基于數(shù)學(xué)模型的尋優(yōu)計(jì)算后，EGR系統(tǒng)優(yōu)化設(shè)計(jì)方案可以降低發(fā)動(dòng)機(jī)動(dòng)力性和燃油經(jīng)濟(jì)性的損失。相對(duì)于初始設(shè)計(jì)方案，在全工況范圍內(nèi)，優(yōu)化方案動(dòng)力性指標(biāo)提升約在1%～3%，燃油經(jīng)濟(jì)性指標(biāo)改善約在1%～5%。相對(duì)于未采用文丘里管的EGR發(fā)動(dòng)機(jī)，在低轉(zhuǎn)速工況內(nèi)，優(yōu)化方案的發(fā)動(dòng)機(jī)充氣效率、動(dòng)力性和燃油經(jīng)濟(jì)性指標(biāo)提升2%～6%；在高轉(zhuǎn)速區(qū)域，由于文丘里管的節(jié)流損失，發(fā)動(dòng)機(jī)充氣效率、動(dòng)力性和燃油經(jīng)濟(jì)性指標(biāo)有2%以內(nèi)的下降。

6 結(jié)束語(yǔ)

a. 通過(guò)2水平析因設(shè)計(jì)，得出文丘里管喉口直徑對(duì)EGR系統(tǒng)流動(dòng)狀態(tài)影響最為顯著的結(jié)論。

b.采用DoE，技術(shù)通過(guò)科學(xué)合理的設(shè)計(jì)EGR系統(tǒng)仿真計(jì)算方案，正確分析模擬計(jì)算結(jié)果，能夠高效準(zhǔn)確地實(shí)現(xiàn)EGR系統(tǒng)最優(yōu)結(jié)構(gòu)方案的設(shè)計(jì)。

c.文丘里管的采用大大提升了EGR系統(tǒng)的工作能力，優(yōu)化方案的EGR率絕對(duì)值提高了5.5%～11.5%。

d. 結(jié)構(gòu)設(shè)計(jì)優(yōu)化方案相對(duì)于初始方案，其全工況下新鮮空氣進(jìn)氣量上升1%～4%，EGR率絕對(duì)值上升1%～3%，實(shí)現(xiàn)了發(fā)動(dòng)機(jī)動(dòng)力性、燃油經(jīng)濟(jì)性和排放性能的最優(yōu)組合。

1 AddyMajewski W，Mafdi K Khair.Diesel Emissions And Their Control.SAE International，2006.

2 高征.車用增壓中冷柴油機(jī)4氣門及EGR系統(tǒng)流動(dòng)過(guò)程的研究:[學(xué)位論文].上海：同濟(jì)大學(xué)，2005.

3 葉年業(yè),劉潔,倪計(jì)民，等.車用汽油機(jī)流動(dòng)過(guò)程模擬及基于DoE的配氣相位優(yōu)化.內(nèi)燃機(jī)工程.2011（04）.

4 倪計(jì)民.汽車內(nèi)燃機(jī)原理.上海：同濟(jì)大學(xué)出版社，1997.

5 TaguchiG.,YokoyamaY.Taguchimethods.design of experiments，1993.

6 杜倩穎.DoE理論及在發(fā)動(dòng)機(jī)優(yōu)化設(shè)計(jì)中的應(yīng)用研究:[學(xué)位論文].上海：同濟(jì)大學(xué)，2010.