基于多物理場耦合的航空發(fā)動機燃油泵液膜密封仿真分析與試驗*

任寶杰 郝木明 胡廣陽 王銘章 蘇志善 李勇凡

(1.東營海森密封技術(shù)有限責任公司 山東東營 257000；2.中國石油大學(華東) 新能源學院 山東青島 266580；3.中國航空發(fā)動機集團航空發(fā)動機動力傳輸重點實驗室 遼寧沈陽 110015；4.中國航發(fā)西安動力控制科技有限公司 陜西西安 710077)

機械密封作為旋轉(zhuǎn)機械的重要部件，在工業(yè)上得到了廣泛應(yīng)用[1-2]。隨著現(xiàn)代工業(yè)的發(fā)展，機械密封的使用環(huán)境越來越苛刻，對機械密封的性能提出了更高的要求。對于接觸式機械密封，高參數(shù)工況下機械密封存在潤滑失效問題，導致整機故障[3]。

非接觸式液膜密封因其兼具密封性和良好的潤滑特性，能夠?qū)崿F(xiàn)在高參數(shù)工況下的長壽命安全運行。在實際運轉(zhuǎn)中，非接觸式液膜密封端面的潤滑流體膜處于微米級，與密封環(huán)受到熱力影響產(chǎn)生的變形處于同一量級，因此在對高參數(shù)液膜密封進行分析時，應(yīng)考慮熱-流-固多物理場的耦合效應(yīng)。QIU和KHONSARI[4]基于CFD-ACE+得到了螺旋槽機械密封副的溫度和壓力分布。WU等[5]采用FLUENT建立密封環(huán)、周圍沖洗液和其他密封組件的三維流固耦合模型，分析了機械密封的流場和溫度分布規(guī)律。楊丹丹等[6]建立了高溫油泵機械密封的三維力熱耦合分析模型，指出力熱耦合變形中熱變形占主導地位。朱維兵等[7]采用ANSYS Workbench對斜線槽上游泵送機械密封進行單向流固耦合分析，分析了工況參數(shù)對密封環(huán)變形的影響。

為提高非接觸式機械密封的性能，相關(guān)學者針對端面槽型優(yōu)化開展了大量研究工作。在單列螺旋槽的基礎(chǔ)上，LAI[8]提出了人字槽和Y形槽等雙列螺旋槽結(jié)構(gòu)；丁雪興等[9]采用多目標優(yōu)化理論對人字槽密封和螺旋槽密封進行結(jié)構(gòu)優(yōu)化，結(jié)果表明人字槽密封的綜合性能高于螺旋槽密封；MENG等[10]采用三維熱流模型，分析了人字槽結(jié)構(gòu)參數(shù)對密封性能的影響。總體而言，與單列螺旋槽液膜密封相比，人字槽液膜密封同時具有上游泵送和下游泵送的能力，具有更大的流體膜厚度及流體膜剛度，穩(wěn)定性優(yōu)于單列螺旋槽[11-12]。黃偉峰等[13-14]發(fā)現(xiàn)了機械密封在高轉(zhuǎn)速下的溫升和變形對密封性能的削弱效應(yīng)，總結(jié)了機械密封特性參數(shù)變化規(guī)律并應(yīng)用于上游泵送機械密封的設(shè)計優(yōu)化。YANG等[15]基于熱彈流分析方法，對三角形織構(gòu)進行數(shù)值計算后發(fā)現(xiàn)該織構(gòu)形槽型對于密封性能具有明顯提升作用。孫鑫暉等[16]在考慮熱流體動力學特性基礎(chǔ)上，對槽數(shù)、槽深等參數(shù)對于密封性能的影響進行計算分析。

本文作者以某型航空發(fā)動機燃油增壓泵燃油側(cè)機械密封為研究對象，針對燃油側(cè)機械密封溫升過高導致密封失效的問題進行密封端面的結(jié)構(gòu)設(shè)計?；诹髯儗W、摩擦學與潤滑理論、材料學、傳熱學等學科的基本原理，建立了適用于航空發(fā)動機燃油增壓泵用機械密封多物理場仿真分析模型，并通過試驗驗證了人字槽機械密封的性能滿足燃油泵運行要求。

1 理論模型

1.1 物理模型

燃油側(cè)機械密封動靜環(huán)物理模型如圖1所示。

圖1 機械密封結(jié)構(gòu)

圖2所示為人字槽液膜密封端面結(jié)構(gòu)示意圖，人字槽型端面結(jié)構(gòu)特征參數(shù)為槽數(shù)N與槽深hg，槽數(shù)N預設(shè)計范圍為12～24，槽深hg預設(shè)計范圍為16.0～40.0 μm。

圖2 密封端面結(jié)構(gòu)

1.2 考慮粗糙度的雷諾方程

考慮表面粗糙度的、描述機械密封端面液膜壓力分布的控制方程如式(1)所示[13]。

(1)

1.3 粗糙表面微凸體彈塑性接觸力模型

隨微凸體法向變形量逐漸增大，將依次發(fā)生彈性變形、彈塑性變形和塑性變形，同時也產(chǎn)生相應(yīng)的彈性接觸力、彈塑性接觸力和塑性接觸力。描述包括所述3種接觸力階段的理論模型如式(2)所示。

(2)

式中：pt(h)、pet(h)、pept(h)、ppt(h)為總接觸壓力、彈性接觸壓力、彈塑性接觸壓力、塑性接觸壓力，Pa；w為微凸體法向變形量；w1、w2為初始屈服點臨界法向變形量、完全塑性變形的臨界法向變形量；ys為表面凸峰平均高度線和表面平均高度線之間的距離；z為微凸體的高度(相對于表面凸峰平均高度線)；φs(z)為微凸體高度分布的概率密度函數(shù)(以表面凸峰平均高度線為基準)；h為以表面平均高度線為基準的兩表面平均距離；R為微凸體的曲率半徑；H為軟材料的硬度，Pa；η為微凸體面積密度。

1.4 能量方程

描述機械密封端面液膜溫度分布的控制方程(平均能量方程)如式(3)所示。

(3)

1.5 熱傳導方程

非接觸式機械密封動環(huán)、靜環(huán)和液膜的溫度計算方法與接觸式機械密封動靜環(huán)溫度計算方法相一致，影響密封環(huán)溫度場的主要因素有2個：一是密封環(huán)旋轉(zhuǎn)生成熱量，二是密封靜動環(huán)與液膜及外界流體進行熱量交換。熱量主要由兩部分組成：一是攪拌熱，由于主軸旋轉(zhuǎn)使得流體流動攪拌而產(chǎn)生；二是摩擦熱，因動環(huán)旋轉(zhuǎn)并且液體具有黏性而在端面間隙產(chǎn)生。由于對流傳熱是密封環(huán)同流體熱量交換的主要方式，故而影響環(huán)溫度分布的主要參數(shù)之一便是對流傳熱系數(shù)。

在整個密封系統(tǒng)中主要涉及到的熱傳遞方式有以下幾種：

(1)密封動靜環(huán)端面與間隙流體之間是流-固耦合傳熱；

(2)密封動靜環(huán)內(nèi)部同種介質(zhì)之間通過熱傳導方式傳遞熱量；

(3)密封環(huán)內(nèi)外側(cè)流體與環(huán)的內(nèi)外徑之間是對流傳熱。

溫度場和變形的研究對象是動環(huán)和靜環(huán)，極坐標下彈塑性體的傳熱方程如下：

(4)

式中：T為溫度，℃；Q/υ為內(nèi)部生成熱；cp為比熱容，J/(kg·℃)；k為傳熱系數(shù)，W/(m·℃)。

傳熱方程變?yōu)?/p>

(5)

密封環(huán)的周圍是密封介質(zhì)，它們之間的傳熱方式是對流傳熱。對流邊界條件要求垂直于邊界的傳導熱流等于通過對流傳熱帶走的熱量。

因此對流傳熱邊界滿足如下方程：

(6)

對于耦合端面，要求傳遞的總熱量與生成熱相等。則耦合邊界條件數(shù)學表達式為

(7)

式中：α為對流傳熱系數(shù)，W/(m2·℃)；T∞為環(huán)境溫度，℃；k1、k2為結(jié)構(gòu)1和2的傳熱系數(shù)，W/(m·℃)；T1、T2為結(jié)構(gòu)1和2的溫度，℃。

1.6 對流傳熱系數(shù)

機械密封的傳熱計算主要是指密封動、靜環(huán)與周圍介質(zhì)間傳熱效應(yīng)的計算，其對流傳熱系數(shù)α目前基本均采用經(jīng)驗公式或半經(jīng)驗公式進行計算。

(1)靜環(huán)外圓面與介質(zhì)之間：

Nu=0.023εlRe0.8Pr0.4

Nu=2αl/λl

Re=2vl/υl

(8)

式中:εl為修正系數(shù)，一般可取εl=2；l為靜環(huán)外圓面與密封箱內(nèi)周壁間的徑向間隙；λl為液體介質(zhì)的導熱系數(shù)；υl為液體介質(zhì)的運動黏度；v為靜環(huán)外周周圍介質(zhì)的軸向流速。

(2)動環(huán)外圓面與介質(zhì)之間：

(9)

2 數(shù)值求解方法

采用的仿真計算流程如圖3所示。該計算流程以開啟力與閉合力相平衡為最終的收斂判別依據(jù)，同時包括液膜壓力場收斂、液膜溫度場收斂、密封環(huán)端面溫度分布收斂等3個子判別依據(jù)；滿足上述各判別依據(jù)之后，端面間隙形狀將自行滿足收斂。

圖3 仿真計算流程

表1給出了文中所采用的密封結(jié)構(gòu)參數(shù)和工況參數(shù)。

表1 密封面幾何結(jié)構(gòu)與工況參數(shù)

采用數(shù)字式黏度計測量所購油品在不同溫度下的動力黏度，利用所測數(shù)據(jù)擬合介質(zhì)的黏溫關(guān)系方程(黏溫特性模型)，以用于仿真計算。油品的實測數(shù)據(jù)及擬合公式如圖4所示。

圖4 黏溫特性曲線

表2給出了密封環(huán)溫度場及力熱變形計算相關(guān)物性參數(shù)。

表2 密封動靜環(huán)物性參數(shù)

機械密封密封環(huán)的力、熱邊界條件序號如圖5所示，密封環(huán)計算力、熱邊界條件序號所對應(yīng)的邊界類型如表3所示。

圖5 密封環(huán)邊界序號

表3 密封環(huán)力、熱邊界條件類型

3 計算結(jié)果與分析

設(shè)定轉(zhuǎn)速為9 210 r/min，密封外徑壓力設(shè)定0.60 MPa，環(huán)境溫度設(shè)定25 ℃，密封內(nèi)徑壓力設(shè)定0.10 MPa,計算所得燃油側(cè)平端面接觸式機械密封與人字槽液膜密封性能特征參數(shù)見圖6，人字槽液膜密封計算結(jié)果如圖7所示。

由圖6可知，在計算工況下該人字槽型密封結(jié)構(gòu)泄漏量約0.138 mL/min，空化率約0.97%，端面平均溫度約117.6 ℃，計算泄漏量滿足相關(guān)設(shè)計要求，這與該密封端面結(jié)構(gòu)所具有的上游泵送效果密切相關(guān)。

圖6 平端面密封與人字槽型密封性能特征參數(shù)

與平端面接觸式機械密封進行對比，雖然人字槽型液膜密封泄漏量大于平端面密封，但是人字槽型密封摩擦扭矩與摩擦功率均小于平端面接觸式機械密封，人字槽型密封端面平均溫度低于平端面密封端面平均溫度約50 ℃，端面溫度的降低可相當程度減小密封介質(zhì)高溫變質(zhì)、積碳發(fā)生和密封環(huán)高溫熱裂、皰疤的可能性。另外，與平端面接觸式機械密封相比，由于人字槽型液膜密封摩擦副之間存在一層薄膜流體，在運行工況發(fā)生波動時，薄膜流體具有的剛度特性可實現(xiàn)端面膜厚的自動調(diào)整，以適應(yīng)外界工況的變化；在密封閉合力突然增大時，密封端面薄膜流體的存在可起到減沖緩阻作用，避免密封摩擦副可能出現(xiàn)的接觸撞擊，有利于密封運行穩(wěn)定性的提高。

分析圖7(a)可以得出，雙環(huán)帶人字槽型密封端面間流體膜的壓力分布有2個低壓區(qū)域，分別出現(xiàn)在人字槽的外徑側(cè)槽的槽根處及內(nèi)徑側(cè)槽的槽根處，最低壓力為流體在計算溫度下的飽和蒸汽壓；在人字槽型的岔口尖端處出現(xiàn)高壓區(qū)，其最大壓力約為1.80 MPa。人字槽型液膜密封相當于在同一密封端面開設(shè)不同的兩組螺旋槽，在動環(huán)旋轉(zhuǎn)以及內(nèi)外徑壓差的作用下，端面流體具有較高的動能和壓能，當流體在隨著動環(huán)端面旋轉(zhuǎn)的情況下遇到人字槽岔口尖端壁面，降速升壓，故此處出現(xiàn)高壓區(qū)。由于端面高壓的存在使密封動、靜環(huán)端面間充滿一層厚度為微米級的薄膜流體而維持密封端面摩擦副的非接觸狀態(tài)，從而避免了端面固體直接接觸而引起的密封副摩擦磨損現(xiàn)象，降低了密封環(huán)因摩擦磨損而出現(xiàn)熱裂、皰疤、崩邊等密封失效行為的可能性；且由于該結(jié)構(gòu)通過端面開設(shè)的動壓槽在旋轉(zhuǎn)條件下的黏性剪切作用,把高壓側(cè)泄漏到低壓側(cè)的被密封介質(zhì)再反輸至高壓側(cè)，降低了被密封介質(zhì)由高壓側(cè)向低壓側(cè)的泄漏，大大提高了密封運行可靠性、穩(wěn)定性。

圖7(b)所示為預設(shè)計結(jié)構(gòu)參數(shù)下人字槽型密封端面流體溫度分布?？擅黠@地看出，密封端面流體溫度基本在120～140 ℃左右，靠近外徑側(cè)槽區(qū)溫度略低于靠近內(nèi)側(cè)處槽區(qū)溫度，這是由于靠近內(nèi)徑槽區(qū)面積較大，流體動壓效應(yīng)強，流體與槽壩的摩擦熱高于靠近外徑處溫度。另外雖然該結(jié)構(gòu)密封端面存在一層薄膜流體避免了密封副的直接固體接觸，相對減少了摩擦生熱功率并使密封端面溫度降低，但由于運行工況的高轉(zhuǎn)速與密封燃油介質(zhì)具有較大的黏度，故密封端面溫度仍較大于密封內(nèi)外徑處介質(zhì)溫度。

基于能量守恒，對密封動靜環(huán)及周圍流體進行傳熱計算，得到人字槽型密封動環(huán)、靜環(huán)溫度分布，如圖7(c)、(d)所示。可以發(fā)現(xiàn)，密封動靜環(huán)溫度最大值均位于密封摩擦副接觸端面處，但動環(huán)溫度最值低于靜環(huán)溫度相應(yīng)最值，這是由于動環(huán)材料9Cr18的導熱系數(shù)高于靜環(huán)材料M298的導熱系數(shù)，增大了熱量在動環(huán)內(nèi)部的傳遞，且由于動環(huán)隨轉(zhuǎn)軸的運動增大了動環(huán)與周圍流體的對流換熱，故雖然密封端面流體黏性剪切熱量流入動環(huán)量大于密封靜環(huán)，但動環(huán)內(nèi)大部分區(qū)域溫度仍略低于靜環(huán)溫度。同時可看出，動環(huán)內(nèi)溫度梯度分布明顯大于靜環(huán)內(nèi)溫度梯度分布，這也是由動環(huán)材料9Cr18的導熱系數(shù)高于靜環(huán)材料M298的導熱系數(shù)所決定的。另外，雖然密封外徑側(cè)對流傳熱系數(shù)高于密封內(nèi)徑側(cè)，但由于密封外徑側(cè)介質(zhì)溫度大于內(nèi)徑處溫度，故密封動靜環(huán)外徑側(cè)溫度仍均相對高于內(nèi)徑側(cè)溫度。

圖7(e)、(f)所示為人字槽型密封動環(huán)、靜環(huán)軸向力熱變形?？梢园l(fā)現(xiàn)，密封動環(huán)最大變形發(fā)生在密封摩擦副接觸端面外徑處，該處邊界條件為液膜壓力與溫度分布；密封靜環(huán)最大變形發(fā)生在靜環(huán)遠離接觸端面處，該處為軸向約束邊界條件；密封靜環(huán)最小變形發(fā)生在密封摩擦副接觸端面內(nèi)徑處，該處邊界條件為液膜壓力與溫度分布。分析可知，由于密封動環(huán)材料9Cr18的彈性模量遠大于密封靜環(huán)材料M298的彈性模量，因此密封動環(huán)發(fā)生力熱變形能力弱于密封靜環(huán)，故雖然在密封摩擦副接觸端面處具有相同的力熱邊界條件，但靜環(huán)端面最小變形量6.26 μm仍大于動環(huán)端面最大變形量3.72 μm。根據(jù)密封動、靜環(huán)端面變形計算可知，密封環(huán)變形量與密封端面流體膜厚、端面動壓槽深皆為微米量級，密封環(huán)變形的發(fā)生將通過對端面流體膜厚的影響而進一步對密封性能產(chǎn)生影響。因此在模擬分析時采用了熱流固雙向耦合的方法，即通過計算軟件MatLab與ANSYS的相互調(diào)用，實現(xiàn)密封端面流場、溫度場、密封環(huán)變形的雙向計算，提高了該仿真計算的精確性。

圖7 人字槽型密封熱流固耦合計算結(jié)果

4 試驗驗證分析

試驗總時長2 h，為穩(wěn)態(tài)試驗過程，轉(zhuǎn)速9 210 r/min，介質(zhì)壓力為0.5 MPa。對泄漏量、工裝處溫度、工裝處介質(zhì)溫度及軸套溫度進行測量和記錄，每30 min記錄一次數(shù)據(jù)。試驗前后密封環(huán)對比如圖8所示。

圖8 試驗前后密封環(huán)端面情況

從圖8(b)中可以發(fā)現(xiàn)：動環(huán)端面磨痕呈周向布置，均勻、無斷續(xù)痕跡；端面有明晰的3個區(qū)域，對應(yīng)外側(cè)壩區(qū)、中間槽區(qū)、內(nèi)徑壩區(qū)；對應(yīng)中間槽區(qū)部分磨痕較兩側(cè)為重，符合啟停過程中槽邊刮削所致，但無任何崩邊、掉塊等接觸式密封常見問題。如圖8(b)所示，動環(huán)內(nèi)徑側(cè)劃痕明顯，這與圖7(f)的靜環(huán)變形方向及規(guī)律相對應(yīng)，靜環(huán)軸向變形形成沿內(nèi)徑到外徑的收斂形錐度，導致動環(huán)內(nèi)徑處接觸比壓更大，磨痕相對明顯。這一定程度驗證理論計算結(jié)果的準確性。

如圖9所示是溫度和泄漏量隨時間的變化。在試驗過程中，各處溫度穩(wěn)定，無較大波動。未觀察到有泄漏產(chǎn)生，泄漏量始終為0，可滿足密封使用需求。

圖9 溫度和泄漏量隨時間的變化

5 結(jié)論

(1)以描述粗糙表面的流體動壓潤滑方程、熱效應(yīng)能量方程、粗糙表面微凸體彈塑性接觸力方程、固體傳熱學方程及固體彈性力學方程等基本理論為指導，建立了適用于航空發(fā)動機燃油增壓泵用機械密封的熱流固多物理場耦合理論模型，基于計算軟件MATLAB與ANSYS的耦合求解原理確定了仿真計算流程，給出了機械密封動靜環(huán)的計算物理模型、計算工況參數(shù)、密封環(huán)計算邊界條件，并使用數(shù)字式黏度計測量、擬合分析出燃油黏溫關(guān)系方程。

(2)依據(jù)所建立的燃油增壓泵機械密封熱流固耦合理論模型，開展了人字槽型密封的仿真分析，與平端面接觸式機械密封進行對比，人字槽型密封端面平均溫度降低約50 ℃，有利于密封運行穩(wěn)定性的提高。

(3)對人字槽型機械密封進行了驗證試驗。密封測試過程中，各處溫度始終保持在穩(wěn)定范圍，且未觀察到有泄漏產(chǎn)生，泄漏量始終為0，滿足航空發(fā)動機燃油增壓泵對密封性能的要求。