傅煜 王哲 臧鵬飛 孫晨樂
(同濟大學新能源汽車工程中心,上海 201804)
自由活塞直線發電機換氣品質影響因素研究*
傅煜王哲臧鵬飛孫晨樂
(同濟大學新能源汽車工程中心,上海201804)
基于自由活塞直線發電機系統搭建零維Matlab/Simulink模型和三維AVL Fire模型相耦合的仿真體系,并驗證其精確性。基于該仿真體系分析了自由活塞發動機的運動特性和換氣過程,并研究氣口高度和氣道角度對換氣品質的影響規律。結果表明,大掃排氣口高度、65°主掃氣道傾角和50°副掃氣道仰角有利于提高掃氣效率;小掃排氣口高度、小主掃氣道傾角和大副掃氣道仰角有利于提高缸內湍動能。
主題詞:自由活塞直線發電機換氣品質氣口高度氣道角度
自由活塞直線發電機(FPLG)作為一種新型的發電裝置,由起動發電一體化的直線電機和取消了曲柄連桿機構的自由活塞發動機(FPE)耦合而成。該系統具有壓縮比可變、燃料適應性強、摩擦損耗低等優勢,受到國內外研究機構的廣泛關注[1~4]。研究人員在研究過程中發現,良好、高效的換氣系統是保證FPLG穩定著火的基礎,更是保證系統穩定運行的前提。然而該系統的FPE由傳統二沖程汽油機(TSE)改造而來,其回流掃氣式和氣口式的換氣系統存在換氣過程時間短、新鮮充量氣口倒灌和短路損失等缺點,使換氣品質降低,影響FPE動力性和經濟性,進一步影響系統穩定性。
本文基于FPLG的運行特點,建立零維Matlab/Simu-link模型和三維AVL FIRE模型相耦合的仿真體系,研究氣口高度和氣道角度對FPLG換氣品質的影響。
2.1FPLG系統簡介
FPLG采用雙活塞對置的形式,直線電機位于兩個水平對置的FPE中間,電機動子和發動機活塞連桿機構固定連接,系統結構如圖1所示。

圖1 FPLG系統結構
系統起動時,直線電機處于電動模式,其動子帶動活塞連桿組件以一定頻率直線往復運動,當活塞運動的頻率及行程滿足著火條件時,發動機噴油與點火;當發動機穩定燃燒后,直線電機切換為發電模式,此時活塞連桿組件帶動電機動子直線往復運動,使直線電機持續發電,實現化學能-機械能-電能的依次轉化。
為深入研究FPLG,設計開發了試驗樣機[5],其中FPE由傳統TSE去除曲柄連桿機構并重新設計掃氣箱改造而來,同時改進了供油方式和點火系統。根據FPE的結構與性能及FPLG的起動與發電要求,設計研發永磁三相直線ISG電機,將其與兩邊的FPE相耦合[6]。臺架試驗的主要設備及傳感器參數如表1所列。

表1 試驗用主要設備和傳感器參數
2.2零維Matlab/Simulink模型
基于Matlab/Simulink平臺搭建FPLG系統的仿真模型,以獲得系統的運動特性,為后續的AVL FIRE三維仿真模型提供初始條件和邊界條件。
為驗證Simulink仿真模型的準確性,在直線電機起動力為300 N、起動力換向位置為極限位置前4.5 mm、循環噴油量為1.68 mg、直線電機電磁負載331N、點火位置為極限位置前3 mm的工況下,將仿真結果與樣機試驗結果進行對比,起動階段的活塞位移曲線對比如圖2所示。可見,在系統起動過程中,仿真曲線和實測曲線總體趨勢基本一致,誤差小于6.0%。因此,Simulink仿真模型具有一定的精確性。

圖2 起動階段活塞位移曲線對比
2.3三維AVL FIRE模型
基于發動機實機在Catia中建立換氣系統幾何模型,包括燃燒室、排氣道、掃氣箱。考慮到掃氣箱內氣體流動不作為研究重點,因此將不規則的掃氣箱簡化為等體積圓柱體[7]。鑒于整個幾何模型為對稱模型,為了減少計算量,取模型的一半進行分析。
將幾何模型導入Hypermesh軟件進行二維網格化以獲得高質量面網格。在仿真軟件AVL FIRE中導入面網格以生成靜網格和動網格。動網格包括燃燒室和部分掃氣箱,分別模擬系統運行時燃燒室和掃氣箱的容積變化。動網格和靜網格之間的接觸部分是實際流體聯通區域,由FIRE軟件中的Arbitrary Connect功能實現[8]。最終生成的三維模型由26萬個六面體網格組成。
三維仿真邊界條件如表2所列。部分初始條件如湍動能和耗散率等通過經驗公式計算或取經驗值來初步確定[9]。壓力、溫度和密度等初始條件則根據零維仿真和試驗初步確定,并在三維模型中進行循環迭代計算,直到初始條件的誤差不超過5%。所確定的關鍵初始條件如表3所列。

表2 關鍵邊界參數

表3 關鍵初始條件
缸內流動的仿真結果在目前的試驗條件下難以直接驗證,因此采用整個循環的缸壓結果來間接評價CFD模型的可靠性。將三維CFD計算所得的缸壓曲線與零維仿真和試驗結果進行對比,如圖3所示。可見,3條缸壓曲線總體變化趨勢相近,其峰值誤差小于5.4%,因此三維CFD仿真具有一定可靠性。
該系統的發動機采用回流掃氣方式,掃氣口布置形式為兩個對稱的主掃氣口和兩個對稱的副掃氣口。主副掃氣口各自在縱橫方向上有一定傾斜角(圖4),新鮮充量能沿氣缸壁和活塞頂流動,形成回流并將廢氣擠出。

圖3 三維CFD仿真缸壓驗證結果

圖4 主副掃氣道角度示意
在驗證三維CFD仿真可靠性的基礎上,獲得系統在換氣過程中掃排氣口的質量流量以及缸內廢氣系數和廢氣質量,分別如圖5和圖6所示。在圖5中正值代表氣流流入氣缸,負值代表相反方向(下文同)。

圖5 掃排氣口質量流量示意

圖6 缸內廢氣質量和殘余廢氣系數示意
從排氣口率先開啟到掃氣口開啟為換氣過程的自由排氣階段(94°-126°),此時氣缸內充斥著高溫高壓廢氣,與排氣道形成較高的壓力差,促使缸內廢氣迅速流出。由于缸壓迅速下降至排氣道壓力,在自由排氣階段后期缸內廢氣不再排出,甚至出現排氣口倒灌的現象。此階段共排出47.2%的廢氣。
從掃氣口開啟到掃氣口關閉為強制掃氣階段(126°-219°),掃氣箱內的壓縮氣體與缸內廢氣再次形成壓力差,促使新鮮充量流入氣缸。由于主掃氣道的氣口面積較大,提供的橫向掃氣氣流較為強勁,并與副掃氣道涌出的縱向掃氣氣流協同作用于氣缸,在掃氣口開啟后的20°內迅速將廢氣擠壓至排氣口排出。在強制掃氣階段的中后期再次出現掃氣口和排氣口倒灌的現象。此階段共排出41.6%的廢氣。
從掃氣口關閉到排氣口關閉為過后排氣階段(219°-249°),缸內壓力和溫度與排氣道趨于一致,氣體流動較前兩個階段放緩,共排出3.9%的廢氣。
由于取消了曲柄連桿機構,FPE與TSE的活塞位移曲線呈現不同的特征,如圖7所示,其中Y軸零點處表示氣缸壁上沿位置。可見,FPE的上止點和壓縮比比TSE的高,且在上止點附近的活塞瞬時加速度大,造成活塞在上止點附近停留時間較短。另外,FPE的排氣口開啟相位角較晚,而關閉相位角較早,使得在相同轉速下換氣時間相對較少。

圖7 TSE和FPE活塞位移曲線
為對換氣品質進行量化評價設立3個性能指標:掃氣效率、捕獲率和換氣終了時的缸內湍動能。掃氣效率表明了缸內廢氣被新鮮充量清掃干凈的程度,與峰值缸壓和發動機功率成正比;捕獲率表明新鮮充量從排氣口直接逃逸的情況,用于輔助評價換氣品質;湍動能則是表征氣流運動劇烈程度的重要參數,較大的湍動能有利于提高燃燒時的火焰傳播速度,縮短燃燒持續期,提高放熱率,提升發動機性能[10]。定義掃氣效率和捕獲率公式為:


式中,ηs為掃氣效率;φt為捕獲率;m1為換氣結束后缸內新鮮充量質量;m0為換氣結束后缸內氣體總質量;m2為換氣過程中通過掃氣口的新鮮充量質量。
將傳統發動機在上述FPE的工況(轉速1 800 r/min、壓縮比11)下進行仿真計算并對比換氣品質,結果如表4所列。可見,TSE的掃氣效率略高于FPE,而兩者的捕獲率基本持平,這可能是由兩者活塞位移曲線的不同造成的,TSE在換氣階段的曲軸轉角范圍較大,在同轉速下掃氣更充分;FPE在換氣結束后的缸內湍動能較大,這可能是由其活塞平均速度較高,使缸內空氣運動較為劇烈所致。

表4 兩種發動機換氣品質對比
4.1排氣口高度的影響
排氣口的高度對FPE的換氣過程和燃燒過程具有重要影響,不僅影響掃氣相位角和換氣過程時長,還影響氣口式發動機的有效行程和壓縮比。為探究排氣口高度對換氣過程的影響,設定16 mm、18 mm、20 mm(原排氣口高度)、22 mm、24mm 5種高度。經迭代計算,5種三維模型在換氣開始時的部分初始參數如表5所列。可見,換氣過程開始時的缸壓和缸溫隨著排氣口高度的減小而降低。這是因為排氣口高度的降低使排氣口延遲開啟,使做功行程中缸壓和缸溫下降更多。

表5 不同排氣口高度下的初始條件
圖8為排氣口高度對換氣品質的影響。可見,隨著排氣口高度的增大,掃氣效率逐漸上升,捕獲率逐漸下降。這是因為延遲開啟的排氣口延長了換氣時間,使掃氣更加充分,但也使更多的新鮮充量流出氣缸,造成捕獲率降低。而排氣口高度越小,其換氣結束后的缸內湍動能越大,有利于燃燒過程的快速進行。圖9為3種排氣口高度下通過排氣口的質量流量趨勢圖。可見,在強制排氣階段(126°-219°)排氣道質量流量變化趨勢基本相同,而在自由排氣階段和過后排氣階段,排氣口高度越大,排出氣體越多。另外,在自由排氣階段后期出現了排氣口倒灌現象,其倒灌量隨著排氣口高度的增大即自由排氣時間的延長而增多。

圖8 排氣口高度對換氣品質的影響

圖9 排氣口高度對排氣口質量流量的影響
4.2掃氣口高度的影響
掃氣口高度決定了掃氣口面積,直接影響掃氣氣流在掃氣道內的流通能力。由于副掃氣道的流量較小,故本次仿真對副掃氣口不作改動。對主掃氣口設定9 mm、10 mm、11 mm(原氣口高度)、12 mm、13mm 5種高度,用以探究掃氣口高度對換氣過程的影響。計算結果如圖10和圖11所示。可見,隨著主掃氣口高度增大,掃氣開啟時刻提前,質量流量和掃氣量增大,使得掃氣充分,進而掃氣效率逐步提高。換氣結束后的缸內湍動能則隨著氣口高度的增大而遞減;而捕獲率則呈現波動的趨勢,影響規律不明顯。

圖10 主掃氣口高度對換氣品質的影響
掃氣道布置角度是換氣系統結構的重要參數,決定換氣過程中缸內空氣的運動方式。主掃氣道傾角主要形成缸內渦流,副掃氣道仰角則可形成缸內滾流,兩者共同作用以組織缸內氣流運動并影響換氣品質。

圖11 主掃氣口高度對主掃氣口質量流量的影響
5.1主掃氣道傾角的影響
為探究不同主掃氣道傾角的影響規律,分別建立主掃氣道傾角為45°、50°、55°、60°、65°、70°和75°時的三維模型,在FIRE中進行仿真計算。圖12為主掃氣道傾角對換氣品質的影響。可見,掃氣效率隨著主掃氣道傾角的增大呈先上升后下降的趨勢,并在65°傾角時掃氣效率達到最大值;而捕獲率和湍動能則隨著傾角的增大而逐步下降。傾角越小,缸內湍動能越大,越有利于提高燃燒速度和加快火焰傳播。圖13為不同主掃氣道傾角下通過主掃氣口的氣體質量流量曲線。隨著傾角增大,質量流量就越大,通過掃氣口的新鮮充量也越多。這是因為在氣口面積不變的情況下,氣體的流通能力與傾角成正比。

圖12 主掃氣道傾角對換氣過程的影響

圖13 主掃氣道傾角對主掃氣口質量流量的影響
圖14為不同主掃氣道傾角的三維模型截面1在不同曲軸轉角時的廢氣分布云圖和流場跡線圖。由圖14可知,活塞下行剛開啟掃氣口時大量新鮮充量進入氣缸,在缸內形成大尺度滾流。隨著掃氣口的逐漸開啟,滾流分裂為兩個小滾流,右側小滾流隨著活塞面下沉并破碎消失,左側小滾流則上升并在下止點時重新形成占據整個氣缸的大尺度滾流。從該過程可以看出,傾角越小,左側滾流形成得越快,右側滾流消失得也越快。自下止點起,由于掃氣箱和缸內壓力已持平,掃氣氣流消失甚至出現倒灌現象,由左側滾流發展而來的大尺度滾流將始終占據整個氣缸,并隨著活塞上行直到上止點。其中可見50°傾角的滾流中心上升的比較快。縱觀整個換氣過程可知,主掃氣道傾角越小,缸內氣體流動越劇烈,說明湍流速度越大,與圖12中湍動能的變化趨勢相符。

圖14 不同主掃氣道的廢氣分布云圖和流場跡線圖
缸內廢氣分布情況基本與氣流流向相對應。由于換氣過程前期掃氣氣流沖入氣缸造成氣流紊亂,缸內廢氣分布很不均勻,然而在換氣過程中后期,大尺度滾流的存在使得缸內廢氣分布趨于一致。
5.2副掃氣道仰角的影響
為探究副掃氣道仰角對換氣過程的影響規律,分別建立仰角為35°、40°、45°、50°和55°的三維模型,仿真結果如圖15所示。可見,隨著副掃氣道仰角的增大,掃氣效率先增大再減小,并在50°時具有最大值;捕獲率呈先增大再減小的趨勢,并在40°處出現最大值;而換氣結束后的缸內湍動能則隨著仰角的增大而逐漸增大。圖16為不同副掃氣道仰角下,新鮮充量通過副掃氣口的質量流量變化,可見仰角越大,副掃氣道的氣體流通能力越好,掃氣氣體和倒灌氣體的質量流量都隨之越大。

圖15 副掃氣道仰角對換氣品質的影響

圖16 副掃氣道仰角對副掃氣口質量流量的影響
圖17為不同傾角的三維模型的截面2在不同曲軸轉角時的廢氣分布云圖和流場跡線圖。掃氣口剛打開時,缸內未形成有效大尺度滾流。此時從較大角度的副掃氣道出來的氣流流向更靠近氣缸中心,能更快速的排擠缸內廢氣。這股掃氣氣流流至缸頂后在右側形成了小型滾流。而在缸頂左側出現的小型滾流則主要是由主掃氣道涌出的氣流經左側缸壁流至頂端形成。隨著活塞下行,氣缸右上角形成了較大滾流,其中心逐漸向左轉移,并在下止點時形成占據整個氣缸的大尺度滾流。在此期間,較大仰角模型中的滾流受掃氣氣流的影響,形成稍慢。然而在下止點之后這些大滾流趨于一致,并無太大差別。

圖17 不同副掃氣道的廢氣分布云圖和流場跡線圖
a.由于活塞運動特性的不同,在相同轉速和壓縮比下,相比于傳統二沖程發動機,自由活塞發動機的掃氣效率略低,而缸內湍動能較強。
b.隨著排氣口高度的增加,掃氣效率上升,捕獲率和湍動能則下降;隨著主掃氣口高度的增加,掃氣效率上升,捕獲率波動,而湍動能下降。
c.隨著主掃氣道傾角的增大,掃氣效率先上升后下降,捕獲率和湍動能則下降;隨著副掃氣道仰角的增大,掃氣效率和捕獲率先上升后下降,而湍動能逐漸上升。
d.大掃排氣口高度、65°主掃氣道傾角和50°副掃氣道仰角有利于提升峰值缸壓和功率;小掃排氣口高度、小主掃氣道傾角和大副掃氣道仰角有利于缸內空氣運動程度的加劇和燃燒過程的加速。
1Rinderknecht F,Kock F.A High Efficient Energy Converter for a Hybrid Vehicle Concept.26th Electric Vehicle Symposium,Los Angeles,2012.
2Kosaka H,Akita T,Moriya K,et al.Development of free piston engine linear generator system part 1-investigation of fundamental characteristics.SAE Technical Papers of SAE 2014 World Congress and Exhibition,Detroit,2014-01-1203.
3李慶峰.自由活塞式內燃發電機的研究:[學位論文].上海:上海交通大學,2011.
4欒延龍.直線發動機設計開發及起動優化試驗研究:[學位論文].上海:同濟大學,2013.
5尹兆雷.自由活塞直線發電機系統起動性能及控制策略研究:[學位論文].上海:同濟大學,2015.
6孫晨樂,王哲,鄧俊,等.直線式增程器用直線ISG電機結構設計與優化.汽車工程,2014(9):1132~1139.
7姜一通.自由活塞內燃發電系統換氣過程研究:[學位論文].北京:北京理工大學,2015.
8毛金龍.自由活塞直線發電機工作過程數值仿真試驗研究:[學位論文].北京:北京理工大學,2011.
9黃霞.自由活塞內燃發電機仿真:[學位論文].上海:上海交通大學,2013.
10孫晶晶.汽油機進氣及燃燒過程的多維數值模擬研究:[學位論文].北京:北京交通大學,2008.
(責任編輯晨曦)
修改稿收到日期為2016年6月1日。
Research on Factors Affecting the Scavenging Quality of Free Piston Linear Generator
Fu Yu,Wang Zhe,Zang Pengfei,Sun Chenle
(Clean Energy Vehicle Engineering Center,Tongji University,Shanghai 201804)
Simulation architecture coupled with zero-dimensional model Matlab/Simulink and three-dimensional model AVL Fire was built based on free piston linear generator(FPLG)and its accuracy passed experimental test.The dynamic characteristics and the scavenging process of FPLG were analyzed based on this simulation architecture and the effect of key structure parameters of the scavenging system including port height and port angle on scavenging quality was studied as well.The results show that,large height of scavenging and exhaust ports,tangential inclination of main scavenging port at 65°and axial inclination of sub scavenging port at 50°can increase scavenging efficiency;small height of scavenging and exhaust ports,small tangential inclination of main scavenging port and large axial inclination of sub scavenging port can enhance the turbulent kinetic energy in cylinder.
Free piston linear generator,Scavenging quality,Port height,Port angle
U464.2
A
1000-3703(2016)09-0047-06
上海市科委項目(14DZ1202802)。