王莉娜， 陳國定， 孫恒超

西北工業(yè)大學(xué) 機電學(xué)院， 西安 710072

軸承腔油滴沉積特性及油膜流動特征分析

王莉娜， 陳國定*， 孫恒超

西北工業(yè)大學(xué) 機電學(xué)院， 西安 710072

軸承腔的潤滑和換熱設(shè)計依賴于對腔內(nèi)油氣兩相流動和換熱狀態(tài)的準確理解。針對先前研究工作的不足，在包含油滴碰撞腔壁熱量交換、沉積熱量以及油膜溫度即考慮溫度效應(yīng)的條件下，開展了油滴沉積特性及油膜流動特征分析。首先分析了單個油滴碰撞腔壁沉積特性，確定了油滴的沉積質(zhì)量、動量和熱量。其次在考慮油滴尺寸分布的條件下，通過離散油滴尺寸范圍的方式確定了腔內(nèi)所有油滴碰撞腔壁的沉積油膜質(zhì)量、動量和熱量。最后借助力學(xué)平衡和質(zhì)量守恒條件，以上述參數(shù)為基礎(chǔ)計算了腔壁油膜溫度以及速度和厚度的分布情況。計算結(jié)果表明，隨著轉(zhuǎn)子轉(zhuǎn)速的增加，油滴的質(zhì)量、動量和熱量沉積率有所降低，一次沉積油膜質(zhì)量、沉積油膜動量、沉積油膜熱量以及油膜速度均有所增加，而油膜厚度降低。通過與試驗數(shù)據(jù)的對比表明，提出的軸承腔油滴沉積特性及油膜流動特征分析方法是較為可靠的，并且考慮溫度效應(yīng)使油膜厚度的計算更為準確。分析工作為軸承腔的潤滑和換熱設(shè)計提供了一定的參考依據(jù)。

航空發(fā)動機； 軸承腔； 油滴； 沉積； 油膜； 對流換熱； 溫度

在航空發(fā)動機軸承腔中，潤滑和冷卻滾動軸承的潤滑油被高速旋轉(zhuǎn)的軸承離散成大量尺寸各異的油滴，油滴在腔內(nèi)空氣中高速運動，最終碰撞到腔壁上形成油膜，油膜沿著腔壁流動并從排油口流出軸承腔。為了防止?jié)櫥蛷妮S承腔靜止/旋轉(zhuǎn)界面的間隙泄漏，通過將壓縮空氣從間隙引入軸承腔的方式實現(xiàn)滑油密封。軸承腔內(nèi)空氣、潤滑油以及機械構(gòu)件的溫度差異也比較大。這樣，軸承腔內(nèi)形成了“油滴-空氣-油膜”共存的油氣兩相流動和換熱狀態(tài)。航空發(fā)動機高轉(zhuǎn)速、大推重比及長壽命等發(fā)展趨勢對軸承腔的設(shè)計也提出了新的要求。軸承腔的潤滑和換熱設(shè)計很大程度上依賴于對腔內(nèi)復(fù)雜油氣兩相流動和換熱狀態(tài)的準確理解。因此，近年來很多學(xué)者開展了這方面的研究。

就有關(guān)軸承腔油滴沉積及油膜流動分析的研究而言，Wittig等[1]最先設(shè)計了軸承腔試驗臺，利用超聲波技術(shù)測量了軸承腔壁面油膜厚度，并使用熱電偶測量了油氣介質(zhì)的溫度。Farrall等確定了軸承腔中油滴碰撞腔壁的飛濺與沉積過渡準則，給出了油滴碰撞腔壁后不同變化狀態(tài)的定量描述[2]，并進一步探討了油滴初始條件對油滴碰撞腔壁沉積作用的影響[3-4]。劉登等[5]開展了軸承腔中運動油滴與腔壁碰撞狀態(tài)辨識與特征分析，提出了運動油滴與腔壁斜碰撞條件下的狀態(tài)判斷準則。陳國定等[6]則分析了軸承腔油滴尺寸分布與軸承腔結(jié)構(gòu)和工況參數(shù)的關(guān)系，結(jié)合油滴尺寸分布并且通過油滴運動及沉積特征分析，確定了軸承腔中沉積油膜質(zhì)量和動量。劉亞軍等[7]分析了軸承腔腔壁油膜受到的重力以及空氣剪切力等作用，通過力學(xué)平衡條件分別建立了低轉(zhuǎn)速和高轉(zhuǎn)速條件下腔壁油膜流動分析模型。陳國定等[8]將軸承腔油滴運動和沉積特征以及軸承腔腔壁油膜流動等進行了串行分析，提出了油氣兩相流動狀態(tài)下軸承腔腔壁油膜物理特征研究的一種新途徑。Glahn和Wittig[9]使用激光多普勒測速儀測量了軸承腔腔壁油膜流動速度，建立了腔壁油膜流動速度分析模型，接著Glahn和Wittig[10]從油膜流動理論分析出發(fā)建立了油膜與腔壁對流換熱系數(shù)理論分析方法。趙靜宇等[11]對軸承腔油膜流動進行了非穩(wěn)態(tài)數(shù)值計算，討論了不同轉(zhuǎn)速條件下軸承腔油膜的運動特征。胡劍平等[12]提出了軸承腔腔壁油膜的二維計算模型，對腔壁油膜厚度進行了數(shù)值模擬。然而迄今有關(guān)軸承腔油滴沉積分析工作，很少涉及油滴碰撞腔壁熱量交換以及沉積熱量的分析，這使油滴沉積質(zhì)量和動量以及后續(xù)油膜流動特征分析的準確性有所降低；加之缺少油滴碰撞腔壁沉積熱量的輸入條件，也限制了軸承腔腔壁油膜流動溫度分析工作的實施。

借助之前油滴運動和溫度分析工作確定的油滴碰撞腔壁速度和溫度等參數(shù)，在包含油滴碰撞腔壁沉積熱量與油膜流動溫度的條件下開展了軸承腔油滴沉積特性及油膜流動特征分析。首先確定了區(qū)分單個油滴碰撞腔壁沉積或飛濺的過渡準則，并根據(jù)熱量守恒條件確定了單個油滴的沉積熱量。其次，確定了描述軸承腔中油滴尺寸的質(zhì)量分布密度函數(shù)，并確定了軸承腔中所有油滴碰撞腔壁的沉積油膜質(zhì)量、動量和熱量。繼而確定了油膜流動的平均溫度，完善了油膜流動分析模型并計算了腔壁油膜流動厚度和速度。借助相關(guān)試驗結(jié)果驗證了油膜流動分析的合理性以及考慮油滴碰撞腔壁沉積熱量以及油膜流動溫度影響的必要性。分析工作完善了軸承腔油滴沉積特性及油膜流動特征理論分析體系，有助于提升軸承腔的潤滑和換熱設(shè)計水平。

1 軸承腔單個油滴的沉積特性分析

軸承腔中空氣受轉(zhuǎn)子剪切作用后切向速度很高。油滴以較高的速度自軸承甩出在空氣中運動，受到空氣拖曳力、離心力及與空氣熱量交換作用后，油滴的運動軌跡、速度和溫度會發(fā)生明顯的改變。孫恒超等[13-14]給出了油滴在腔中運動軌跡、速度和溫度的分析方法。本文的工作是在這些工作的基礎(chǔ)上展開的。

軸承腔中油滴在運動至腔壁位置處，會以一定的速度碰撞軸承腔腔壁，并因油滴的直徑、速度以及潤滑油的物性參數(shù)等的不同而產(chǎn)生不同的碰撞狀態(tài)。軸承腔壁面油膜的厚度是很薄的，通常約為0.5 mm[1]，厚度很薄的油膜對油滴沉積特性影響較小，因此分析油膜碰撞腔壁時忽略了油膜的存在。Mundo等[15]研究了液滴碰撞固體壁面產(chǎn)生的沉積現(xiàn)象，根據(jù)他的研究結(jié)果可推斷軸承腔中油滴與腔壁碰撞的狀態(tài)。

Mundo等[16]通過觀測液滴與固體壁面碰撞沉積現(xiàn)象，建立了液滴碰撞固體壁面的沉積準則。根據(jù)其研究結(jié)論可以確定軸承腔中油滴與腔壁碰撞的沉積準則，并通過無量綱飛濺參數(shù)B表示，即

(1)

式中：Wed和Red分別為油滴入射腔壁時的韋伯?dāng)?shù)和雷諾數(shù)，表達式為

(2)

(3)

式中：ρd、μd和σ分別為油滴入射腔壁時的密度、動力黏度和表面張力，將油滴入射腔壁時的溫度Td代入潤滑油物理性質(zhì)隨溫度的變化關(guān)系式即可確定這些參數(shù)；dd和udr分別為油滴入射腔壁時的直徑和徑向速度。Td和udr可以通過文獻[13-14]給出的數(shù)值方法計算獲得。先前研究工作中未考慮溫度效應(yīng)，確定潤滑油物性參數(shù)是以油滴從軸承甩出時的溫度Td0為參考溫度的。本文確定油滴與腔壁碰撞沉積準則時，考慮了油滴自軸承向腔壁運動過程中溫度變化的實際情況，包含了溫度變化對潤滑油物性、油滴直徑以及油滴軌跡和速度的影響，這種包含溫度效應(yīng)的油滴沉積特性分析方法提升了確定油滴入射腔壁狀態(tài)參數(shù)的計算精度，也可以在一定程度上提升油滴沉積特性分析的準確性。

圖1 軸承腔油滴與腔壁碰撞沉積準則Fig. 1 Impact deposition criteria for droplet/housing interaction in a bearing chamber

油滴與腔壁碰撞后不破裂，油滴的全部質(zhì)量轉(zhuǎn)移至腔壁形成油膜，油滴的徑向動量在碰撞時被消耗，切向動量轉(zhuǎn)移至腔壁油膜，油滴的熱量及碰撞腔壁時從腔壁吸收的熱量轉(zhuǎn)移至腔壁油膜。

油滴與腔壁碰撞產(chǎn)生破裂，破裂油滴直接沉積部分的質(zhì)量轉(zhuǎn)移至腔壁形成油膜，其徑向動量被消耗，但切向動量轉(zhuǎn)移至腔壁油膜。同樣地，沉積部分的熱量及碰撞腔壁時從腔壁吸收的熱量轉(zhuǎn)移至腔壁油膜。碰撞腔壁飛濺形成的二次油滴中，部分再次沉積為油膜，另一部分從通風(fēng)口排出軸承腔。一般而言二次油滴的質(zhì)量比較少，油滴碰撞腔壁飛濺的二次油滴的速度相比于碰撞腔壁前油滴的速度也會小很多，并且二次油滴的速度是雜亂無章的，故可認為二次油滴的動量相互抵消，二次油滴再次沉積時沒有動量轉(zhuǎn)移至腔壁油膜。二次油滴再次沉積時，沉積部分的熱量及碰撞腔壁時從腔壁吸收的熱量也會轉(zhuǎn)移至腔壁油膜。

為了更好地評估油滴與腔壁碰撞沉積質(zhì)量、動量和熱量的變化情況，在此引入油滴碰撞腔壁的質(zhì)量沉積率、動量沉積率和熱量沉積率的概念。

油滴與腔壁碰撞后，不包含飛濺的二次油滴再次沉積的質(zhì)量，沉積部分的質(zhì)量與碰撞前油滴質(zhì)量的比值為油滴一次質(zhì)量沉積率η1。計入二次油滴再次沉積的質(zhì)量后，沉積的總質(zhì)量與碰撞前油滴質(zhì)量的比值為油滴質(zhì)量沉積率η。Chen等[8]根據(jù)Farrall等[4]的分析確定了油滴質(zhì)量沉積率η1和η的表達式為

(4)

(5)

式中：εL為二次油滴再次沉積質(zhì)量與二次油滴總質(zhì)量的比值[8]。

油滴動量沉積率是油滴與腔壁碰撞后，沉積部分轉(zhuǎn)移到油膜的動量與油滴初始動量的比值。分析中不考慮二次油滴的動量轉(zhuǎn)移，油滴動量沉積率為

(6)

式中：udtz為油滴碰撞腔壁時的切向速度；ud0為油滴自軸承甩出時的速度。

Pasandideh等[17]分析了液滴碰撞熱壁面時對壁面的冷卻現(xiàn)象，并給出了液滴碰撞熱壁面時對壁面冷卻效率與液滴參數(shù)的關(guān)系。根據(jù)Pasandideh等的研究方法可以分析軸承腔中油滴碰撞軸承腔腔壁時的熱量沉積率。定義油滴碰撞腔壁時對腔壁的冷卻效率為

(7)

式中：qc為油滴碰撞腔壁時與腔壁的熱交換量。

油滴碰撞腔壁時對腔壁的冷卻效率εc與油滴碰撞腔壁時的工況參數(shù)之間的關(guān)系為[17]

(8)

式中：Prd為油滴的普朗特數(shù)，與潤滑油的運動黏度和導(dǎo)溫系數(shù)有關(guān)；θa為潤滑油與腔壁的接觸角。

根據(jù)油滴碰撞腔壁時對腔壁的冷卻效率可以推導(dǎo)油滴碰撞腔壁時的熱量沉積率。

1) 油滴碰撞腔壁不破裂(B≤57.7)

(9)

油滴在撞擊壁面時，垂直于壁面的動能會轉(zhuǎn)化為熱能，忽略油滴碰撞腔壁時潤滑油與空氣間的熱交換量以及黏性耗散熱后，根據(jù)熱量守恒條件，即油滴碰撞腔壁時熱量的變化等于潤滑油/腔壁的熱交換量與油滴垂直于壁面的動能之和，其表達式為

(10)

通過式(10)可以獲得油滴碰撞腔壁不破裂時，油滴的熱量沉積率為

(11)

2) 油滴碰撞腔壁破裂(B>57.7)

(12)

同樣根據(jù)熱量守恒條件，可以確定

(1-η)mdcvl(Tsd)Tsd-mdcvl(Td)Td=

(13)

通過式(13)可以獲得油滴碰撞腔壁破裂時，油滴的熱量沉積率為

(14)

目前還缺少關(guān)于油滴碰撞腔壁后形成的二次油滴的溫度Tsd的理論和試驗研究，其他研究領(lǐng)域中也未見關(guān)于液滴撞壁形成的二次液滴溫度研究的報道?？梢酝茢嗟氖嵌斡偷蔚臏囟萒sd介于撞壁前油滴溫度Td和壁面溫度Tw之間，因此假定Tsd=(Td+Tw)/2。

2 考慮油滴尺寸分布的油滴沉積及油膜流動分析

第1節(jié)分析了單個油滴碰撞腔壁時的沉積質(zhì)量、動量和熱量，然而軸承腔腔壁油膜的形成是軸承甩射的所有油滴沉積累積作用的結(jié)果。軸承甩射的油滴數(shù)量極多、尺寸不均勻，首先需要確定軸承腔中油滴尺寸的數(shù)學(xué)描述方式。

2.1 軸承腔油滴尺寸分布

根據(jù)文獻[6]的分析，軸承腔潤滑油在潤滑軸承后被分解成的大量油滴的質(zhì)量分布密度函數(shù)可以用R-R函數(shù)描述為

(15)

(16)

(17)

2.2 考慮油滴尺寸分布的油滴沉積分析

進行沉積油膜質(zhì)量、動量和熱量分析時需要考慮油滴尺寸分布后確定所有油滴與腔壁碰撞沉積的質(zhì)量、動量和熱量，將油滴分布直徑離散后，逐一確定每個直徑區(qū)間的沉積質(zhì)量、動量和熱量，最后累加獲得全部油滴的沉積油膜特征。

1) 沉積油膜質(zhì)量

將軸承甩射的油滴的直徑范圍[dmin,dmax]劃分為J個等距區(qū)間，第j個子區(qū)間直徑范圍為[dj,dj+1]?？梢愿鶕?jù)油滴質(zhì)量分布密度函數(shù)fm確定第j個子區(qū)間中油滴的質(zhì)量為

(18)

式中：Mb為時間段tb內(nèi)軸承甩射的油滴的總質(zhì)量，Mb=mltb，ml為潤滑油的質(zhì)量流量。

第j個子區(qū)間一次沉積油膜質(zhì)量Mc1j和沉積油膜質(zhì)量Mcj分別為

(19)

(20)

式中：η1j和ηj分別為區(qū)間[dj,dj+1]加權(quán)平均直徑dj+0.5的油滴的一次質(zhì)量沉積率和質(zhì)量沉積率。

時間tb內(nèi)一次沉積油膜質(zhì)量Mc1和沉積油膜質(zhì)量Mc分別是所有油滴直徑區(qū)間的一次沉積油膜質(zhì)量Mc1j和沉積油膜質(zhì)量Mcj之和，其表達式為

(21)

(22)

進一步可以確定單位時間的一次沉積油膜質(zhì)量mc1和沉積油膜質(zhì)量mc:

(23)

(24)

2) 沉積油膜動量

與確定沉積油膜質(zhì)量類似，可以通過離散油滴直徑區(qū)間的方式確定油滴碰撞腔壁時的沉積油膜動量。根據(jù)第j個子區(qū)間中油滴的質(zhì)量Mj以及直徑dj+0.5的油滴碰撞腔壁時的動量沉積率ξj可以確定第j個子區(qū)間中油滴的沉積油膜動量為

(25)

同樣地，疊加所有直徑區(qū)間油滴的沉積油膜動量即可確定所有油滴的沉積油膜動量，其表達式為

(26)

并且，單位時間內(nèi)沉積油膜動量為

(27)

3) 沉積油膜熱量

與確定沉積油膜質(zhì)量和動量的方式相同，可以確定考慮油滴尺寸分布后所有油滴的最終的單位時間沉積油膜熱量q。在此，直接給出其計算表達式為

(28)

2.3 軸承腔油膜流動狀態(tài)分析

在沉積油膜質(zhì)量、動量和熱量分析基礎(chǔ)上，可以進行軸承腔腔壁油膜速度、厚度和溫度的分析。根據(jù)油滴碰撞腔壁后單位時間的沉積油膜熱量q可確定油膜的平均溫度為

(29)

王軍[18]根據(jù)軸承腔腔壁油膜流動的力學(xué)平衡和質(zhì)量守恒條件確定了腔壁油膜流動厚度和平均速度。根據(jù)其分析方法確定軸承腔腔壁油膜厚度和速度的求解方程組為

(30)

(31)

本文提出的軸承腔腔壁油膜厚度、速度以及溫度分析方法考慮了油滴在腔中運動以及沉積過程中的溫度效應(yīng)，并且本小節(jié)分析油膜流動厚度和速度時也考慮了油膜溫度對潤滑油密度和動力黏度的影響。這種處理方式有助于提升對油膜流動狀態(tài)理論預(yù)測的準確性。

3 計算結(jié)果與討論

進行油滴碰撞腔壁沉積特性以及油膜流動分析時，除討論工況參數(shù)的影響外，各個參數(shù)的數(shù)值如下：轉(zhuǎn)子轉(zhuǎn)速ns=10 000 r/min，軸承腔高度hb=28 mm，轉(zhuǎn)子半徑rs=62 mm，軸承腔寬度wb=20 mm，進氣溫度Tgi=105 ℃，潤滑油溫度Tli=60 ℃，使用的潤滑油為4109號航空潤滑油，其物性參數(shù)隨溫度的變化關(guān)系可以參閱文獻[19]，腔壁溫度Tw=210 ℃，轉(zhuǎn)子壁面溫度Ts=95 ℃。

3.1 單個油滴碰撞腔壁沉積特性

圖2是轉(zhuǎn)子轉(zhuǎn)速對直徑dd=100 μm的單個油滴沉積特性的影響。從圖2(a)中可以看出，油滴的一次質(zhì)量沉積率雖小于1，但是較為接近于1。這說明在所涉及的工況條件下，油滴碰撞腔壁時雖然產(chǎn)生了破裂，但破裂形成的二次油滴的總質(zhì)量很少。而質(zhì)量沉積率幾乎等于1，這說明油滴碰撞腔壁時，幾乎所有質(zhì)量都沉積至腔壁上。油滴沉積至腔壁的動量約為油滴從軸承甩出時初始動量的一半。從考慮溫度效應(yīng)和不考慮溫度效應(yīng)的對比可以看出，考慮溫度效應(yīng)后，油滴的質(zhì)量沉積率和動量沉積率均略有降低。這或許是因為考慮溫度效應(yīng)后，油滴被空氣加熱，使其直徑有所增大，大直徑的油滴碰撞腔壁破裂程度加劇，致使油滴沉積質(zhì)量減少，加之考慮溫度效應(yīng)后油滴運動速度更低，所以油滴沉積動量會更少。此外，考慮溫度效應(yīng)會影響潤滑油的黏度，并進一步對油滴的沉積特性產(chǎn)生一定影響。

從圖2(b)中可以看出，由于碰撞腔壁時與腔壁存在強烈的換熱作用，油滴的熱量明顯增加，在碰撞腔壁后，油滴沉積的熱量約為其從軸承甩出時初始熱量的2.8倍。

從圖中油滴沉積特性隨轉(zhuǎn)子轉(zhuǎn)速的變化可以看出，隨著轉(zhuǎn)子轉(zhuǎn)速的增加，油滴的質(zhì)量、動量和熱量沉積率均有所降低。推斷其原因或許是，轉(zhuǎn)子轉(zhuǎn)速增加后，油滴速度增加、動量增大，與腔壁碰撞更為劇烈，油滴破裂程度增加，所以一次質(zhì)量沉積率降低。一次質(zhì)量沉積率的降低也會使動量沉積率降低。轉(zhuǎn)子轉(zhuǎn)速增加后，油滴在空氣中停留時間變短，從空氣中吸收的熱量減少，碰撞腔壁時吸收的熱量變化不大，所以熱量沉積率會因運動過程從空氣吸收熱量的降低而減少。

圖2 轉(zhuǎn)子轉(zhuǎn)速對單個油滴沉積特性的影響Fig. 2 Effect of shaft speed on single oil droplet deposition

3.2 考慮油滴尺寸分布的油滴沉積特性

圖3所示為轉(zhuǎn)子轉(zhuǎn)速對沉積油膜特征的影響。從圖3(a)可以看出，考慮油滴尺寸分布的油滴碰撞腔壁一次沉積油膜質(zhì)量與沉積油膜質(zhì)量較為接近。這說明軸承腔中油滴首次碰撞腔壁時幾乎全部沉積，破裂飛濺形成的二次油滴的質(zhì)量很小。從具體數(shù)據(jù)來看，就所涉及的工況條件，每秒約有1.2 g的潤滑油從通風(fēng)口流出軸承腔。而沉積油膜動量約為潤滑油從軸承甩出時動量的一半。對比考慮溫度效應(yīng)和不考慮溫度效應(yīng)的結(jié)果，考慮溫度效應(yīng)后除一次沉積油膜質(zhì)量略有降低外，沉積油膜質(zhì)量和沉積油膜動量變化均不是很大。但只有考慮溫度效應(yīng)才可獲得沉積油膜熱量，從圖3(b)可以看出，油滴碰撞腔壁時從腔壁吸收了較多的熱量。

從圖3(a)和圖3(b)中可以看出轉(zhuǎn)子轉(zhuǎn)速對沉積油膜特征的影響。隨著轉(zhuǎn)子轉(zhuǎn)速的增加，一次沉積油膜質(zhì)量和沉積油膜動量均有所增加。這是因為轉(zhuǎn)子轉(zhuǎn)速增加后，油滴的平均直徑更小，小直徑油滴碰撞腔壁后破裂形成二次油滴的質(zhì)量更少，所以一次沉積油膜質(zhì)量增加。轉(zhuǎn)速增加后，油滴速度是增加的，所以沉積時的動量增加，因此沉積油膜動量增加。沉積油膜質(zhì)量隨轉(zhuǎn)速幾乎沒有變化，這是因為絕大部分油滴的質(zhì)量沉積率均接近1，因此沉積油膜質(zhì)量也就接近于潤滑油流量。雖然圖2(b)中表現(xiàn)出由于轉(zhuǎn)子轉(zhuǎn)速增加、油滴停留時間變短使單個油滴的熱量沉積率降低，但是轉(zhuǎn)速增加后，油滴平均直徑變小，潤滑油離散成的油滴與氣相介質(zhì)的接觸面積增大，并且后者占主導(dǎo)因素，所以整體上油滴從氣相介質(zhì)中吸收了更多的熱量。相應(yīng)地，油滴數(shù)量增加、油滴直徑降低，均會使油滴碰撞腔壁時，從腔壁吸收的總熱量增加，所以沉積油膜熱量也明顯增加。

圖3 轉(zhuǎn)子轉(zhuǎn)速對沉積油膜特征的影響Fig. 3 Effect of shaft speed on deposition oil film characteristic

3.3 軸承腔壁面油膜流動厚度、速度和溫度

圖4 轉(zhuǎn)子轉(zhuǎn)速對軸承腔腔壁油膜流動的影響Fig. 4 Effect of shaft speed on oil film flow of bearing chamber housing

Gorse等[20]利用軸承腔試驗臺并借助電容傳感器測量了軸承腔圓周坐標θ為60°、135°、240°、270°和315°五處位置的油膜厚度。試驗中軸承腔的轉(zhuǎn)子半徑為64 mm、腔高10 mm、腔寬15 mm。在潤滑油量為100 L/h、進氣量為10 g/s、進油和進氣溫度均為100 ℃、轉(zhuǎn)速12 000 r/min的條件下測量的油膜厚度如圖5所示。作者根據(jù)其試驗條件，開展了軸承腔壁面油膜厚度的理論計算，并且同時給出了是否考慮油滴運動以及油滴沉積過程中溫度效應(yīng)的兩種計算結(jié)果，理論計算值與試驗值的比較如圖5所示。從圖中可以看出，提出的軸承腔壁面油膜厚度分析方法預(yù)測的油膜厚度與試驗值較為接近，考慮溫度效應(yīng)的計算值與不考慮溫度效應(yīng)大約有7%的差異，并且考慮溫度效應(yīng)使計算值整體來看更為接近試驗值。這也說明，在開展軸承腔油滴運動和油膜沉積分析時考慮潤滑油的溫度效應(yīng)是有必要的。在θ=315°處，考慮溫度效應(yīng)后理論結(jié)果偏離試驗結(jié)果的原因或許為θ=315°接近排油口，排油口潤滑油和空氣的流出會對油膜流動存在一定的干擾。

圖5 軸承腔壁面油膜厚度計算值和試驗值[18]的比較Fig. 5 Comparison of oil film thickness of bearing chamber housing between calculation data and experimental data[18]

4 結(jié) 論

1) 隨著轉(zhuǎn)子轉(zhuǎn)速的增加，油滴碰撞腔壁的質(zhì)量、動量和熱量沉積率均有所降低；一次沉積油膜質(zhì)量、沉積油膜動量以及沉積油膜熱量均有所增加；油膜的平均流動速度增加、油膜厚度降低。

2) 較之不考慮溫度效應(yīng)，考慮溫度效應(yīng)后油滴的質(zhì)量沉積率和動量沉積率均略有降低、一次沉積油膜質(zhì)量略有降低。

3) 通過與試驗數(shù)據(jù)的比較，證明提出的分析模型能較為準確地預(yù)測軸承腔壁面油膜的流動狀態(tài)；并且考慮溫度效應(yīng)使計算值更為接近試驗值，這也說明在開展軸承腔油滴運動和油膜沉積分析時考慮潤滑油的溫度效應(yīng)是有必要的。

然而，本文建立軸承腔油滴沉積特性分析方法時考慮的情況仍較為理想，對軸承腔十分復(fù)雜的油滴沉積現(xiàn)象做了一些必要的簡化。雖然模型具有一定的有效性，但是該模型仍然存在進一步深入改進之處并需在今后工作中予以關(guān)注，例如，進行油滴碰撞腔壁沉積分析時需要討論油膜的影響、利用試驗手段準確地獲得油滴沉積部分的溫度以及二次油滴的溫度等。

[1] WITTIG S, GLAHN A, HIMMELSBACH J. Influence of high rotational speeds on heat transfer and oil film thickness in aero-engine bearing chambers[J]. Journal of Engineering for Gas Turbines and Power, 1994, 116(2): 395-401.

[2] FARRALL M, HIBBERD S, SIMMONS K. Modelling oil droplet/film interaction in an aero-engine bearing chamber[C]//9th International Conference on Liquid Atomization and Spray Systems. Sorrento: Institute for Liquid Atomization and Spray System, 2003: 1-8.

[3] FARRALL M, SIMMONS K, HIBBERD S, et al. A numerical model for oil film flow in an aeroengine bearing chamber and comparison to experimental data[J]. Journal of Engineering for Gas Turbines and Power, 2006, 128(1): 111-117.

[4] FARRALL M, HIBBERD S, SIMMONS K. The effect of initial injection conditions on the oil droplet motion in a simplified bearing chamber[J]. Journal of Engineering for Gas Turbines and Power, 2008, 130(1): 1-7.

[5] 劉登， 陳國定， 方龍， 等. 運動油滴/固體壁面斜碰撞的狀態(tài)辨識及特征分析[J]. 航空學(xué)報, 2015, 36(4): 1359-1366.

LIU D, CHEN G D, FANG L, et al. State identification and characteristics analysis of oil droplet’s oblique impact onto solid wall[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(4): 1359-1366 (in Chinese).

[6] 陳國定， 陳薄， 王軍， 等. 軸承腔中油滴尺寸分布和沉積特性[J]. 機械工程學(xué)報, 2011, 47(10): 139-144.

CHEN G D, CHEN B, WANG J, et al. Oil droplet size distribution and deposition properties in bearing chamber[J]. Journal of Mechanical Engineering, 2011, 47(10): 139-144 (in Chinese).

[7] 劉亞軍， 陳國定， 王軍. 高轉(zhuǎn)速下軸承腔內(nèi)壁油膜流動建模[J]. 航空動力學(xué)報, 2010, 25(8): 1900-1905.

LIU Y J, CHEN G D, WANG J. Flow modeling of oil film on bearing chamber wall at high-speed[J]. Journal of Aerospace Power, 2010, 25(8): 1900-1905 (in Chinese).

[8] CHEN G D, SUN H C, WANG J. Research into configuration and flow of wall oil film in bearing chamber based on droplet size distribution[J]. Chinese Journal of Aeronautics, 2011, 24(3): 355-362.

[9] GLAHN A, WITTIG S. Two-phase air/oil flow in aero engine bearing chambers: Characterization of oil film flows[J]. Journal of Engineering for Gas Turbines and Power, 1996, 118(3): 578-583.

[10] GLAHN A, WITTIG S. Two-phase air/oil flow in aero-engine bearing chambers-assessment of an analytical prediction method for the internal wall heat transfer[J]. International Journal of Rotating Machinery, 1999, 5(3): 155-165.

[11] 趙靜宇， 劉振俠， 胡劍平， 等. 軸承腔內(nèi)壁油膜運動特性的數(shù)值研究[J]. 推進技術(shù), 2013, 35(1): 25-32.

ZHAO J Y, LIU Z X, HU J P, et al. Numerical study on motion characteristics of lubricating oil film of bearing chamber wall[J]. Journal of Propulsion Technology, 2013, 35(1): 25-32 (in Chinese).

[12] 胡劍平， 劉振俠， 呂亞國， 等. 軸承腔壁面液膜厚度數(shù)值模擬與實驗測量[J]. 航空計算技術(shù), 2013, 43(2): 68-71.

HU J P, LIU Z X, LU Y G, et al. Comparison of simulation and measurement for film thickness in bearing chamber[J]. Aeronautical Computing Technique, 2013, 43(2): 68-71 (in Chinese).

[13] 孫恒超， 陳國定， 王莉娜， 等. 軸承腔油滴含率及油滴相與空氣能量傳遞分析[J]. 航空學(xué)報, 2016, 37(3): 1074-1087.

SUN H C, CHEN G D, WANG L N, et al. Oil droplets fractions and oil droplets/air energy transfer in bearing chamber[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(3): 1074-1087 (in Chinese).

[14] SUN H C, CHEN G D, ZHANG Y H, et al. Theoretical and experimental study on the motion and thermal states of oil droplet in a bearing chamber[J]. Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering (in press).

[15] MUNDO C, SOMMERFELD M, TROPEA C. Droplet-wall collisions: Experimental studies of the deformation and breakup process[J]. International Journal of Multiphase Flow, 1995, 21(2): 151-173.

[16] MUNDO C, TROPEA C, SOMMERFELD M. Numerical and experimental investigation of spray characteristics in the vicinity of a rigid wall[J]. Experimental Thermal and Fluid Science, 1997, 15(3): 228-237.

[17] PASANDIDEH F M, AZIZ S D, CHANDRA S, et al. Cooling effectiveness of a water drop impinging on a hot surface[J]. International Journal of Heat and Fluid Flow, 2001, 22(2): 201-210.

[18] 王軍. 航空發(fā)動機軸承腔油滴運動及油膜形成研究[D]. 西安： 西北工業(yè)大學(xué), 2010: 58-59.

WANG J. Study of oil droplet motion and the form of oil film in an aeroengine bearing chamber[D]. Xi’an: Northwestern Polytechnical University, 2010: 58-59 (in Chinese).

[19] 林基恕. 航空發(fā)動機設(shè)計手冊 第12冊 傳動及潤滑系統(tǒng)[M]. 北京： 航空工業(yè)出版社， 2002: 456-457.

LIN J S. Design handbook of aero engine, volume 12, transmission and lubrication systems[M]. Beijing: Aviation Industrial Press, 2002: 456-457 (in Chinese).

[20] GORSE P, BUSAM S, DULLENKOPF K. Influence of operating condition and geometry on the oil film thickness in aeroengine bearing chambers[J]. Journal of Engineering for Gas Turbines and Power, 2006, 128(1): 103-110.

王莉娜女， 博士研究生。主要研究方向： 航空發(fā)動機密封與潤滑。

E-mail: wangxiweigood@163.com

陳國定男， 博士， 教授， 博士生導(dǎo)師。主要研究方向： 潤滑與密封、 機電系統(tǒng)熱分析和現(xiàn)代機械設(shè)計理論與方法。

Tel.: 029-88493929

E-mail: gdchen@nwpu.edu.cn

孫恒超男， 博士研究生。主要研究方向： 航空發(fā)動機潤滑系統(tǒng)設(shè)計。

E-mail: shc361@163.com

URL：www.cnki.net/kcms/detail/11.1929.V.20151126.1429.018.html

Characteristicsanalysisofoildropletdepositionandoilfilmflowinabearingchamber

WANGLi’na,CHENGuoding*,SUNHengchao

SchoolofMechanicalEngineering,NorthwesternPolytechnicalUniversity,Xi’an710072,China

Thelubricationandheattransferdesignsofanaeroenginebearingchamberdependontheunderstandingaboutthestatesofoil/airtwo-phaseflowandheattransferinchamber.Toovercometheweaknessofpreviousstudy,weinvestigatesintothecharacteristicsofoildropletdepositionandoilfilmflowinabearingchamberconsideringtheheattransferofoildroplet/housinginteraction,depositionheatandoilfilmtemperature,namely,theeffectoftemperature.Firstly,thedepositioncharacteristicsofasingleoildropletareanalyzed,andthedepositionmass,momentumandheatareobtained.Secondly,withtheoildropletsizedistribution,thedepositionoilfilmmass,momentumandheatofallthedropletsimpactingonhousingareobtainedbydispersingthediameterrangeofthedroplets.Finally,withtheaboveparameters,theoilfilmtemperature,averagevelocityandthicknessarecalculatedbymeansofforcebalanceandmassconservationconditions.Theresultsshowthatwithanincreasingshaftspeed,theoildropletdepositionratesofmass,momentumandheatdecrease,andthefirstdepositionoilfilmmass,depositionoilfilmmomentum,depositionoilfilmheatandoilfilmvelocityincrease,whileoilfilmthicknessdecreases.Thecomparisonofoilfilmthicknessbetweentheoreticaldataandexperimentaldatashowsthattheproposedtheoreticalmodelisreliableandtheconsiderationoftemperatureeffectimprovesthecalculationaccuracyaboutoilfilmthickness.Thepresentworkprovidessome

forthedesignofbearingchamber.

aeroengine;bearingchamber;droplets;deposition;oilfilm;heattransfer;temperature

2015-09-30;Revised2015-10-20;Accepted2015-11-10;Publishedonline2015-11-261429

NationalNaturalScienceFoundationofChina(51275411)

.Tel.:029-88493929E-mailgdchen@nwpu.edu.cn

2015-09-30;退修日期2015-10-20;錄用日期2015-11-10; < class="emphasis_bold">網(wǎng)絡(luò)出版時間

時間:2015-11-261429

www.cnki.net/kcms/detail/11.1929.V.20151126.1429.018.html

國家自然科學(xué)基金 (51275411)

.Tel.:029-88493929E-mailgdchen@nwpu.edu.cn

王莉娜， 陳國定， 孫恒超． 軸承腔油滴沉積特性及油膜流動特征分析J． 航空學(xué)報,2016,37(10):3159-3169.WANGLN,CHENGD,SUNHC.CharacteristicsanalysisofoildropletdepositionandoilfilmflowinabearingchamberJ.ActaAeronauticaetAstronauticaSinica,2016,37(10):3159-3169.

http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

10.7527/S1000-6893.2015.0302

V233.4

A

1000-6893(2016)10-3159-11