吳超琦 王 軍 沈潔陽 張英杰 劉海鷗 尉 玉
(中國航發(fā)商用航空發(fā)動機(jī)有限責(zé)任公司 上海 200241)
附件齒輪箱作為航空發(fā)動機(jī)傳動系統(tǒng)的重要單元體之一,能夠在發(fā)動機(jī)啟動階段提供扭矩[1]。該齒輪箱采用多級輪系,其傳動鏈復(fù)雜、傳遞功率大,齒輪工作時,齒面之間的高副接觸和滑動摩擦?xí)a(chǎn)生大量的熱量。熱量若不能及時散發(fā)出去會使齒輪表面溫度急劇增加,齒面磨損,進(jìn)而導(dǎo)致齒輪失效,齒輪的失效將直接影響到整個航空發(fā)動機(jī)的運轉(zhuǎn)[2]。
為確保齒輪箱正常工作,噴嘴需向齒輪嚙合區(qū)的齒面提供足夠的潤滑油量,航空發(fā)動機(jī)傳動齒輪屬于高速齒輪,多采用噴油潤滑的方式對齒輪進(jìn)行潤滑和冷卻[3]。王延忠等[4]通過數(shù)值模擬的方法分析了齒輪旋轉(zhuǎn)流場對射流的影響,并通過實驗進(jìn)行了驗證。錢雪凌等[5]利用Fluent軟件對復(fù)雜齒輪箱噴射潤滑進(jìn)行了研究,得到了噴嘴布置位置不同對于潤滑效果的影響。MASSINI等[6]利用VOF模型模擬了潤滑油射流撞擊到單個齒輪面的過程,評估了射流潤滑的阻力矩。MEHDIZADEN等[7]利用CFD方法計算了齒輪的風(fēng)阻并得到了齒輪周圍高速旋轉(zhuǎn)下的流場。PAYRI等[8]研究指出噴嘴內(nèi)部的幾何形狀不同,會影響其對高速齒輪的冷卻性能。蘇媛媛等[9]利用一維仿真軟件Flowmaster對滑油系統(tǒng)各處噴嘴和節(jié)流嘴尺寸給出了建議,以匹配流量要求。葛玉柱等[10]利用Flowmaster軟件對齒輪箱潤滑系統(tǒng)進(jìn)行了建模,并結(jié)合試驗驗證了系統(tǒng)設(shè)計的準(zhǔn)確性。王超等人[11]基于Flowmaster軟件對發(fā)動機(jī)整個潤滑系統(tǒng)進(jìn)行了分析,確定了各處噴嘴的尺寸,匹配了各油道的壓力分布。
綜上,目前對齒輪箱噴油潤滑的研究主要集中在齒輪嚙合區(qū)的潤滑和一維系統(tǒng)仿真兩方面。然而航空發(fā)動機(jī)傳動齒輪箱內(nèi)部結(jié)構(gòu)復(fù)雜,除噴嘴和齒輪兩大部件外,還設(shè)有滑油管路,管路上安裝噴嘴;并且為了空間避讓,滑油管路還會設(shè)有折彎,整體結(jié)構(gòu)極為復(fù)雜,齒輪箱內(nèi)部滑油管路的布置方式對齒輪箱的供油和潤滑起著重要作用。本文作者以某航空發(fā)動機(jī)的附件齒輪箱“連通折彎式”復(fù)雜油路結(jié)構(gòu)作為研究對象,通過三維流體分析,預(yù)測發(fā)動機(jī)典型供油壓力下,齒輪箱噴嘴的流量、流量系數(shù)和流向等關(guān)鍵參數(shù),并建立了該油路的壓力-流量數(shù)學(xué)模型。
圖1給出了附件齒輪箱的供油模型圖。齒輪箱總共設(shè)置有6個噴嘴為齒輪嚙合區(qū)和齒輪軸承提供潤滑和冷卻,確保附件齒輪箱工作運轉(zhuǎn)。圖2所示是原方案噴嘴2內(nèi)部結(jié)構(gòu)局部放大圖。
圖1 齒輪箱油路模型
圖2 噴嘴2內(nèi)部結(jié)構(gòu)局部放大
圖3給出了附件齒輪箱油路結(jié)構(gòu)原方案的流體域模型和體網(wǎng)格模型,油路結(jié)構(gòu)較為復(fù)雜,呈“連通折彎式多噴嘴結(jié)構(gòu)”,采用全六面體網(wǎng)格對油路進(jìn)行網(wǎng)格劃分。
圖3 原方案的計算模型與體網(wǎng)格
潤滑油物性參數(shù)取自與附件齒輪箱流量試驗溫度一致的80 ℃下所對應(yīng)的參數(shù),密度為948.7 kg/m3,黏度為0.007 8 Pa·s。計算域入口采用壓力入口(Pressure Inlet)邊界類型,給定表壓0.3 MPa,計算域出口采用壓力出口(Pressure Outlet)邊界類型,給定背壓0。
采用壓力基求解器進(jìn)行穩(wěn)態(tài)求解,用Realizablek-ε湍流模型[12-13]封閉方程組,壁面函數(shù)采用增強(qiáng)型壁面函數(shù)[14](Enhanced Wall Treatment),壓力-速度耦合采用SIMPLE算法[15]。
為了驗證網(wǎng)格無關(guān)性,以附件齒輪箱油路原方案的流體域模型為研究對象,劃分了不同疏密性的體網(wǎng)格模型,體網(wǎng)格的總數(shù)分別是425萬、650萬和792萬,入口和出口指定相同的邊界條件。計算結(jié)果發(fā)現(xiàn),3個不同的體網(wǎng)格方案對應(yīng)的入口體積流量分別是9.42、9.74、9.77 L/min??梢?,體網(wǎng)格總數(shù)從650萬增加至792萬時計算結(jié)果已變化不大。結(jié)合計算資源和時間考慮,文中選取體網(wǎng)格總數(shù)為650萬的方案進(jìn)行仿真計算和結(jié)果分析。
對于附件齒輪箱油路性能評估方式是在特定的工況點下對其開展流量試驗,考察進(jìn)油口的潤滑油體積流量是否滿足要求?;谖闹薪⒌哪P停瑢⒃囼灥玫降牧髁拷Y(jié)果和計算結(jié)果進(jìn)行比較,如表1所示,可以看到二者吻合較好,說明了文中建立模型的準(zhǔn)確性。
表1 計算結(jié)果和試驗結(jié)果比較
對于附件齒輪箱油路模型,在結(jié)構(gòu)設(shè)計完成初期,并未開展流量試驗,因此將原方案計算結(jié)果與設(shè)計要求做了對標(biāo)校核。根據(jù)設(shè)計要求,進(jìn)油口的潤滑油流量需達(dá)到11.5 L/min,將仿真計算得到的體積流量和設(shè)計要求進(jìn)行了對比,如表2所示。可見原方案的計算結(jié)果和設(shè)計要求相差較大,計算結(jié)果比設(shè)計要求的體積流量明顯偏低,說明該附件齒輪箱的油路結(jié)構(gòu)設(shè)計無法滿足要求,因此需進(jìn)行三維流場分析,找出體積流量偏低的原因。
表2 原方案計算結(jié)果和設(shè)計要求對比
圖4給出了該附件齒輪箱油路的壓力分布云圖,可看出整個油路壓力變化顯著的地方有兩處:一處是所有噴嘴的噴孔段,壓力出現(xiàn)了突然下降的情況,這是由于噴點區(qū)域流通面積突然減小,局部阻力損失突然變得很大,導(dǎo)致壓力損失較大導(dǎo)致的;另一處是圖4中黑框所示的區(qū)域,在噴嘴2兩側(cè)的油路出現(xiàn)了壓力損失的情況。
圖4 原方案的壓力云圖
提取不同噴嘴進(jìn)口側(cè)和出口側(cè)的壓力損失,如表3所示??梢钥闯?,噴嘴2進(jìn)、出口側(cè)管路的壓力損失較大,說明潤滑油在噴嘴2處出現(xiàn)了非常大局部阻力損失,該處噴嘴結(jié)構(gòu)需要改進(jìn)優(yōu)化。
表3 各個噴嘴進(jìn)、出口側(cè)的壓力損失
為了找出噴嘴2局部區(qū)域壓力損失異常大的原因,分析了原方案整個附件齒輪箱油路的內(nèi)部流動軌跡,結(jié)果如圖5所示??芍?,當(dāng)潤滑油從入口流入油路后,在管路段的流線較為平直,因為管路段流通面積不變,受到的流動阻力很小,故流線的狀態(tài)接近于“平推流”,不會產(chǎn)生任何渦旋;但是當(dāng)潤滑油流入各個噴嘴時,由于受到噴嘴壁面的影響,流動狀態(tài)發(fā)生了改變,在噴嘴內(nèi)部出現(xiàn)了渦旋的流動狀態(tài)。
圖5 原方案的內(nèi)部流線
為了更加清晰地查看不同噴嘴流入和流出側(cè)的流動情況,在圖5中將紅框中4個區(qū)域的局部流動情況進(jìn)行了放大。對比4個區(qū)域的流場變化情況可以看到,A區(qū)域中的流動劇烈,流出側(cè)的管路內(nèi)形成的渦流清晰可見,速度相比其他3個區(qū)域明顯更大;B、C和D區(qū)域中的流線速度變化相對較小,流出側(cè)的管路內(nèi)幾乎沒有旋流。
圖6給出了噴嘴2流入和流出側(cè)中心截面的流線圖??梢钥吹?,當(dāng)潤滑油流入、流出時,管路流入和流出側(cè)的有效流通面積較小,速度突然增加,局部阻力損失驟然增大,形成了較為強(qiáng)烈的渦旋,流動耗散較大,壓力損失較大。因此噴嘴2的流出側(cè)和流入側(cè)的壓力損失很大,影響了整個附件齒輪箱油路的入口流量。因此,需對噴嘴2的局部結(jié)構(gòu)進(jìn)行優(yōu)化。
圖6 噴嘴2流入和流出側(cè)的流線圖
原方案的三維仿真計算分析結(jié)果表明,由于噴嘴2區(qū)域的流通面積不足導(dǎo)致局部壓力損失較大,影響了整個附件齒輪箱油路的入口流量,故對噴嘴2的局部結(jié)構(gòu)進(jìn)行調(diào)整,將噴嘴2端面A在原方案的基礎(chǔ)上向下移動4.5 mm,增加油路有效流通面積,如圖7所示。
圖7 優(yōu)化方案中噴嘴2內(nèi)部結(jié)構(gòu)局部放大
將優(yōu)化方案的模型采用相同的邊界條件進(jìn)行計算,結(jié)果如表4所示。附件齒輪箱油路的進(jìn)口潤滑油流量為10.8 L/min,相比起原方案的9.74 L/min有了較大增幅,且計算結(jié)果和設(shè)計要求更加接近。
表4 優(yōu)化方案計算結(jié)果和設(shè)計要求對比
圖8給出了附件齒輪箱油路優(yōu)化方案的壓力云圖,可以看到,在各個噴嘴左右側(cè)的進(jìn)出口管路段上,并沒有出現(xiàn)很明顯的壓力變化,壓力損失較大的區(qū)域仍在每個噴嘴的噴孔段處,這是由于流通面積突然減小導(dǎo)致的。但是在優(yōu)化方案中并沒有出現(xiàn)噴嘴2兩側(cè)的管路段上壓力顯著變化的情況。
圖8 優(yōu)化方案的壓力云圖
表5給出了優(yōu)化方案中各個噴嘴左右兩側(cè)的進(jìn)出口管路壓力損失??梢钥闯?,各個噴嘴兩側(cè)管路的壓力損失逐步減小,整個附件齒輪箱油路的壓力分布更加合理。
表5 優(yōu)化方案中各個噴嘴進(jìn)、出口側(cè)的壓力損失
圖9給出了附件齒輪箱油路優(yōu)化方案的內(nèi)部流線圖,可以看到,A區(qū)域流入側(cè)和流出側(cè)的內(nèi)部流線速度變化不大,在進(jìn)入和流出側(cè)的管路內(nèi),原方案中的渦旋流動已消失,流動狀態(tài)更加接近于理想的“平推流”狀態(tài);B、C 和D 3個區(qū)域的流動狀態(tài)和A區(qū)域內(nèi)的流動狀態(tài)相似。
圖9 優(yōu)化方案的內(nèi)部流線
圖10給出了優(yōu)化方案噴嘴2流入和流出側(cè)中心截面的速度矢量圖??梢娫黾恿魍娣e后,局部阻力損失降低,速度的最大值隨之降低,流入和流出側(cè)管道內(nèi)的流動比較平直;而在噴嘴2內(nèi)部還存在著局部區(qū)域的渦旋流動,這是由于流動受到噴嘴結(jié)構(gòu)和流出側(cè)管路的位置限制所致。
圖10 噴嘴2流入和流出側(cè)的流線圖
從以上對附件齒輪箱油路優(yōu)化方案的內(nèi)部流動分析可以看出,原方案油路入口供油量距離設(shè)計要求相差較大是由于噴嘴2結(jié)構(gòu)設(shè)計不合理,流道內(nèi)的流通面積過小,局部阻力損失較大。基于三維仿真分析,提出了相應(yīng)有效的優(yōu)化方案。優(yōu)化方案的油路入口供油量和設(shè)計要求較為相近,驗證優(yōu)化方案的有效性。
附件齒輪箱的流量試驗僅針對典型的進(jìn)口壓力,即0.3 MPa下的潤滑油體積流量試驗,對其他不同進(jìn)口壓力下的潤滑油體積流量未進(jìn)行試驗。前文已經(jīng)通過試驗驗證了文中模型的準(zhǔn)確性,基于此,通過改變?nèi)肟谶吔鐥l件,采用模型計算得到不同進(jìn)口壓力下對應(yīng)的入口潤滑油體積流量,如表6所示。
表6 不同進(jìn)口壓力對應(yīng)的潤滑油體積流量
計算的入口壓力范圍是0.05~0.4 MPa,總共8個壓力點,對8個壓力點下對應(yīng)的滑油體積流量曲線進(jìn)行擬合,如圖11所示,得到擬合方程如下:
圖11 不同進(jìn)口壓力下潤滑油體積流量曲線
Y=-32.238X2+38.193X+2.324 3
式中:Y對應(yīng)附件齒輪箱油路的入口潤滑油流量;X對應(yīng)油路的進(jìn)口壓力。
從圖11可知,該曲線的擬合度為0.998 8,擬合度較好。附件齒輪箱油路在不同進(jìn)口壓力下對應(yīng)潤滑油體積流量呈現(xiàn)二次函數(shù)關(guān)系,在0.05~0.4 MPa的壓力范圍內(nèi),對應(yīng)的潤滑油最小和最大體積流量分別是4.02和12.54 L/min。
基于三維仿真分析,對該附件齒輪箱“連通折彎式”復(fù)雜油路結(jié)構(gòu)建立了壓力-流量的數(shù)學(xué)模型,該模型為不同進(jìn)口壓力下的潤滑油體積流量選擇,以及為附件齒輪箱性能評估提供了數(shù)據(jù)支撐。
對某航空發(fā)動機(jī)附件齒輪箱呈“連通折彎式”的復(fù)雜油路結(jié)構(gòu)模型進(jìn)行三維仿真分析,基于所建立的模型分析了齒輪箱進(jìn)口潤滑油流量不足的原因并提出優(yōu)化方案。主要結(jié)論如下:
(1)該齒輪箱因噴嘴結(jié)構(gòu)的設(shè)計不當(dāng),導(dǎo)致齒輪箱油路結(jié)構(gòu)中的局部流通面積較小,使得該區(qū)域內(nèi)產(chǎn)生渦流,局部阻力損失變大,造成油路進(jìn)口處的滑油流量偏小。
(2)通過對齒輪箱油路中噴嘴結(jié)構(gòu)的優(yōu)化,增大附件齒輪箱油路流通面積,減小了油路的壓力損失,提高了油路進(jìn)口的潤滑油流量,使其滿足設(shè)計要求。
(3)優(yōu)化后的附件齒輪箱油路結(jié)構(gòu)中,壓力損失較大的區(qū)域出現(xiàn)每個噴嘴的噴孔段處,但各個噴嘴兩側(cè)管路段的壓力變化并不大,整個油路的壓力分布更合理,各個管路段內(nèi)的流動狀態(tài)接近于理想的“平推流”形式。
(4)通過三維仿真計算,建立了附件齒輪箱油路結(jié)構(gòu)的壓力-流量數(shù)學(xué)模型,發(fā)現(xiàn)入口壓力為0.05~0.4 MPa時,入口壓力和體積流量呈二次函數(shù)關(guān)系。