李 源，吳亞鋒，樸雅慶

（西北工業(yè)大學(xué) 陜西 西安 710072）

二沖程發(fā)動機具有功率大、旋轉(zhuǎn)平穩(wěn)、振動小、噪音低、維護簡單等優(yōu)點，特別適用于低速、中低空的小型無人機、動力翼傘等飛行器的應(yīng)用。

缸內(nèi)直噴（GDI）發(fā)動機將噴油嘴置于氣缸之內(nèi)，通過高壓燃油泵將燃料精確地噴注進燃燒室內(nèi)，使其霧化的更加細致，并根據(jù)不同的點火時刻將燃油混合氣點燃，具有節(jié)省燃油、減少廢氣排放、提升動力性能。

目前國內(nèi)的二沖程直噴系統(tǒng)的研究剛剛起步，主要圍繞二沖程摩托車發(fā)動機進行的。對于無人機二沖程發(fā)動機直噴系統(tǒng)，為了提高缸內(nèi)過量充氣系數(shù)，改善進氣效率，提高動力性能，通常需要對進氣系統(tǒng)進行優(yōu)化設(shè)計。 GT-Power是發(fā)動機工作過程模擬計算軟件，現(xiàn)被世界上大多數(shù)發(fā)動機和汽車制造廠商使用。該軟件采用有限體積法進行流體計算，計算步長自動可調(diào)，并且有強大的輔助建模前處理工具，自帶有豐富的燃燒模型和優(yōu)化設(shè)計功能，可較好地滿足二沖程發(fā)動機直噴系統(tǒng)建模和分析的需要。

文中針對一臺用于某型無人機直噴發(fā)動機，基于其結(jié)構(gòu)和尺寸，使用GT-Power軟件建立了仿真模型。通過運用亥姆霍茲理論設(shè)計諧振進氣系統(tǒng)結(jié)構(gòu)，并利用軟件自帶的優(yōu)化工具，分析了諧振進氣系統(tǒng)對發(fā)動機動力性能的影響。

1 直噴系統(tǒng)結(jié)構(gòu)及建模

1.1 直噴系統(tǒng)結(jié)構(gòu)

直噴系統(tǒng)工作循環(huán)有：空氣經(jīng)由空氣濾清器、節(jié)氣門、進氣管、簧片閥、曲軸箱，在掃氣口掃氣的作用下，進入氣缸的空氣與燃油混合，按照設(shè)定的點火提前角將油氣混合氣點燃，廢氣經(jīng)掃氣口掃氣排出。直噴系統(tǒng)主要參數(shù)如表1所示。

表1 直噴系統(tǒng)主要參數(shù)Tab.1 Main parameters of the direct injection system

1.2 直噴系統(tǒng)建模

基于直噴發(fā)動機結(jié)構(gòu)和尺寸，使用GT-Power軟件建立了曲軸箱模塊、氣缸模塊、噴射系統(tǒng)模塊、進排氣系統(tǒng)模塊。

燃燒計算仿真選擇了準(zhǔn)維燃燒模型[1－2]，該模型可以考慮氣缸形狀、點火時間及燃料屬性等因素，較SIWiebe燃燒模型能更準(zhǔn)確地預(yù)測噴油、點火等參數(shù)對燃燒放熱率的影響。

進排氣系統(tǒng)流動計算從N－S方程出發(fā)，并簡化為一維控制體積進行計算[3]。傳熱計算選用Woschi模型[4]，在摩擦損失壓力計算中模型采用Chen－Flynn的關(guān)系式，摩擦損失壓力為活塞速度和峰值壓力的函數(shù)[5]。

水平對置發(fā)動機屬于V型發(fā)動機，在曲軸箱中氣缸布置角度設(shè)置為180度。在建模時，還要明確如下的參數(shù)設(shè)置[6]：

1）網(wǎng)格的粗細（離散化長度）應(yīng)趨于均勻，對進/排氣系統(tǒng)離散長度設(shè)置分別為0.4×缸徑/0.5×缸徑。

2）節(jié)氣門、簧片閥的參考流通面積設(shè)置為定值，即零升程時，流量系數(shù)為零。

3）由于進排氣流動、傳熱等情況很復(fù)雜，本文使用推薦進氣端口的溫度為450 K；排氣端口的溫度為550 K。

4）循環(huán)初始條件相同，環(huán)境溫度和壓力在海拔高度0 m和1 000 m設(shè)定為300 K/0.1 MPa和294 K/0.09 MPa。

1.3 直噴發(fā)動機模型的驗證

針對已建立的模型系統(tǒng)，通過對其氣缸壓力和發(fā)動機扭矩進行仿真，并與試驗結(jié)果進行了對比，圖1為發(fā)動機在轉(zhuǎn)速5 000 r/min時的缸壓對比曲線。圖2為不同轉(zhuǎn)速工況點的發(fā)動機扭矩對比曲線。

圖1 缸壓仿真結(jié)果與試驗結(jié)果對比Fig.1 Comparison of the cylinder pressure simulation results and test results

圖2 扭矩仿真結(jié)果與試驗結(jié)果對比Fig.2 Comparison of simulation results and experimental results of torque

從圖1中可以看出，在燃燒階段仿真計算與試驗的缸壓曲線基本重合，表明了所建模型能較準(zhǔn)確反映發(fā)動機內(nèi)部燃燒過程的基本特征。由圖2可以看出，最大扭矩的仿真值和試驗值分別為62 N·m和63 N·m。兩者在3 600 r/min至4 600 r/min時誤差較大，在3 800 r/min時達到最大誤差4.6%，均小于5%，表明了所建模型能夠較準(zhǔn)確模擬直噴發(fā)動機外部動力特征，能夠用于對發(fā)動機進氣系統(tǒng)的設(shè)計和優(yōu)化。

2 諧振進氣系統(tǒng)結(jié)構(gòu)設(shè)計

利用亥姆霍茲原理，將諧振器并聯(lián)安裝于進氣口/節(jié)氣門與空氣濾之間，諧振器與氣道的連接方式主要有共軛式和插入式，增加諧振器的進氣系統(tǒng)在特定的轉(zhuǎn)速內(nèi)，進氣壓力波與其門的進氣周期協(xié)調(diào)，增高進氣歧管壓力，從而增加進氣量。

2.1 諧振器管徑計算

僅考慮管道摩擦的影響，Kastner[7]給出了最小管徑的估算公式：

式中：S為行程；d為缸徑（mm）；n為轉(zhuǎn)速。

由上式，將S=52 mm，d=65 mm和n=3 000 r/min代入，可得到dmin=25.8mm。但綜合考慮管道摩擦、波動效應(yīng)和安裝要求，最終取dmin=26.8 mm。

2.2 諧振器容積計算

1）空氣中音速 c（m/s）與溫度 t（℃）的關(guān)系：

2）發(fā)動機轉(zhuǎn)速n與諧振頻率fn的關(guān)系：

i為發(fā)動機缸數(shù)；n為發(fā)動機轉(zhuǎn)速（r/min），τ為沖程數(shù)。

3）諧振頻率計算：

式中，S 為插接管面積，S=π（d/2）2，mm2；LP為插接管當(dāng)量長度，mm；V為諧振器容積，m3。則可得諧振器容積：

2.3 進氣管管長優(yōu)化方案

在確立諧振器的結(jié)構(gòu)參數(shù)后，還需要對進氣管部分進行優(yōu)化。本文為了降低改進設(shè)計的難度及模具等的開發(fā)成本，只對進氣總管和歧管長度進行選定，分別設(shè)計縮短10%、20%和加長10%、20%4個方案。

3 設(shè)計優(yōu)化方案仿真分析

3.1 諧振器對動力性能的影響

增加諧振器后，在轉(zhuǎn)速區(qū)（2 000～4 000 r/min）和轉(zhuǎn)速區(qū)（5 000～6 000 r/min）充氣效率有所提高，如圖3所示。這是由于進氣流在進氣關(guān)閉時，進氣瞬時發(fā)生逆反回流，緊接著的下個進氣行程吸入新的氣流，二者在進氣道沖撞，增加諧振器緩沖了該頻率區(qū)氣流的諧振干擾，提高了缸內(nèi)充氣量。

圖3 諧振器對充氣效率的影響Fig.3 Effect of resonator on charging efficiency

3.2 進氣總管長度對扭矩和充氣效率的影響

增加諧振器，扭矩和充氣效率隨著進氣總管長度變化，如圖4、圖5所示。由圖可知，在低轉(zhuǎn)速區(qū)和轉(zhuǎn)速在（4 000～5 000 r/min）影響明顯。結(jié)合無人機設(shè)計要求，最終選定縮短20%進氣總管長度的方案。

圖4 扭矩隨進氣總管長度變化曲線Fig.4 Torque curve with the intake manifold length change

3.3 進氣歧管長度對扭矩和充氣效率的影響

通過增加諧振器，改變進氣總管長度，扭矩和充氣效率得到提高。進氣歧管的長度變化也會影響發(fā)動機的動力性能，其扭矩和充氣效率變化分別如圖6、7所示。

在低轉(zhuǎn)速區(qū)（2 000～3 000 r/min）和中轉(zhuǎn)速區(qū)（4 000～5 000 r/min）對發(fā)動機動力性能及充氣效率影響明顯。對于縮短10%和20%歧管長度的方案，在中低轉(zhuǎn)速通過進氣諧振作用，進氣壓力波峰值增大，使得進氣量增多，回流空氣量減少。

4 結(jié)束語

文中通過對一無人機用四缸兩沖程汽油直噴發(fā)動機的仿真，其研究結(jié)果表明可以利用亥姆霍茲原理可以對諧振進氣系統(tǒng)進行優(yōu)化設(shè)計。增加諧振器后充氣效率有所提高，進氣管長度優(yōu)化后扭矩有所提高，燃油消耗率亦有所降低，解決了無人機進氣系統(tǒng)設(shè)計優(yōu)化問題，并使得無人機動力性能進一步提高。

圖5 充氣效率隨進氣總管長度變化曲線Fig.5 The volumetric efficiency curve with the intake manifold length change

圖6 扭矩隨進氣歧管長度變化曲線Fig.6 Torque curve changes with the length of the intake manifold

圖7 充氣效率隨進氣歧管長度變化曲線Fig.7 The charging efficiency curve changes with the length of the intake manifold

[1]伍小軍.缸內(nèi)直噴汽油機的準(zhǔn)維燃燒模型研究 [D].武漢:華中科技大學(xué)，2002.

[2]李向榮.內(nèi)燃機燃燒科學(xué)與技術(shù)[M].北京:北京航空航天大學(xué)出版社，2012.

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

[4]周松，王銀燕.內(nèi)燃機工作過程仿真技術(shù)[M].北京:北京航空航天大學(xué)出版社，2012.

[5]戴正興，宋義忠，崔毅，等.汽油機可變氣門正時技術(shù)仿真及方案分析[J].汽車技術(shù)，2008（3）:22-26.DAI Zhengxing，SONG Yizhong，CUI Yi，et al.Gasoline engine with variable valve timing technology simulation and program analysis[J].Automotive Technology，2008（3）:22-26.

[6]黎幸榮.發(fā)動機排氣噪聲模擬及測試方法的研究[D].合肥:合肥工業(yè)大學(xué)，2010.

[7]侯獻軍，鞏學(xué)軍，方丹，等.LJ276M電噴汽油機進氣系統(tǒng)設(shè)計及優(yōu)化[J].武漢理工大學(xué)學(xué)報：交通科學(xué)與工程版，2009，33（4）:635-638.HOU Xian-jun，GONG Xue-jun，F(xiàn)ANG Dan，et al.LJ276M EFIgasoline engine intake system design and optimization[J].Wuhan University of Technology：Transportation Science&Engineering Edition，2009，33（4）:635-638.