(河北工業(yè)大學(xué)機(jī)械工程學(xué)院 天津 300130)

浮環(huán)軸承以其可靠性高、制造容易、成本低廉等優(yōu)點(diǎn)，被廣泛地應(yīng)用于渦輪增壓器、燃?xì)廨啓C(jī)等高速轉(zhuǎn)子系統(tǒng)中[1]。然而由于轉(zhuǎn)軸與浮動(dòng)環(huán)的旋轉(zhuǎn)使內(nèi)外油膜承受剪切力產(chǎn)生摩擦功耗，導(dǎo)致潤(rùn)滑油溫度升高，而溫度的上升一方面直接降低潤(rùn)滑油的黏度，另一方面浮環(huán)受熱膨脹間接導(dǎo)致內(nèi)外間隙的變化，從而影響浮環(huán)軸承潤(rùn)滑靜特性。

近年來(lái)，許多學(xué)者對(duì)浮環(huán)軸承潤(rùn)滑特性進(jìn)行了深入研究。秦超等人[2]基于雷諾方程并結(jié)合隨機(jī)粗糙模型建立粗糙形狀的浮環(huán)軸承模型，采用有限差分法對(duì)模型進(jìn)行求解，得到浮環(huán)軸承潤(rùn)滑過(guò)程中的油膜厚度和油膜壓力分布。王軍事等[3]采用有限差分法聯(lián)立求解內(nèi)、外油膜的Reynolds方程、膜厚方程和浮環(huán)彈性變形方程，得到在不同轉(zhuǎn)速和偏心率下浮環(huán)的彈性變形量，研究浮環(huán)彈性變形對(duì)浮環(huán)軸承潤(rùn)滑特性的影響。

DELIGANT等[4-5]在不同的轉(zhuǎn)速下對(duì)渦輪增壓器浮環(huán)軸承進(jìn)行了磨損能耗與扭矩測(cè)試試驗(yàn)，并通過(guò)建立CFD模型，計(jì)算得到流體摩擦能耗，最后通過(guò)實(shí)驗(yàn)驗(yàn)證了所建模型的正確性。

當(dāng)轉(zhuǎn)子系統(tǒng)高速旋轉(zhuǎn)時(shí)，浮環(huán)軸承內(nèi)外油膜溫度都會(huì)升高，在黏溫效應(yīng)的影響下，會(huì)使得潤(rùn)滑油黏度降低，進(jìn)而影響浮環(huán)軸承的潤(rùn)滑靜特性。ANDRéS和KERTH[6]建立了熱效應(yīng)對(duì)浮環(huán)軸承動(dòng)靜特性參數(shù)影響的理論模型，并進(jìn)行了實(shí)驗(yàn)驗(yàn)證，指出軸承間隙將會(huì)受到熱效應(yīng)的影響，進(jìn)而會(huì)影響浮動(dòng)環(huán)轉(zhuǎn)速以及其內(nèi)外油膜壓力分布。裴世源等[7]建立了浮環(huán)軸承的穩(wěn)態(tài)熱流體動(dòng)力潤(rùn)滑模型，計(jì)算了典型工況下軸承的動(dòng)靜特性參數(shù)，研究了等溫、導(dǎo)熱和絕熱情況下環(huán)速比、溫升等關(guān)鍵參數(shù)隨轉(zhuǎn)速的變化規(guī)律，指出基于導(dǎo)熱模型的計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合良好，驗(yàn)證了熱效應(yīng)對(duì)于浮環(huán)軸承的靜動(dòng)特性具有顯著影響。張文靜等[8]建立了更為精確的渦輪增壓器浮環(huán)軸承系統(tǒng)溫度預(yù)測(cè)模型，并與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，證明了溫度預(yù)測(cè)模型的準(zhǔn)確性。張立偉[9]考慮了熱效應(yīng)對(duì)潤(rùn)滑油黏度和油膜間隙的影響，得到浮環(huán)徑向變形隨離心力、油膜溫度變化而變化的規(guī)律。

由于浮環(huán)是一個(gè)很輕的薄壁零件，當(dāng)內(nèi)外油膜溫度升高時(shí)，浮環(huán)受熱膨脹，會(huì)影響內(nèi)外油膜間隙，進(jìn)而影響浮環(huán)軸承潤(rùn)滑靜特性。郭紅等人[10]建立浮環(huán)的有限元模型，通過(guò)ANSYS求解得到浮環(huán)所受的內(nèi)外油膜應(yīng)力以及浮環(huán)在應(yīng)力作用下產(chǎn)生的應(yīng)變。結(jié)果表明:當(dāng)主軸轉(zhuǎn)速不太高、偏心率不大時(shí)，浮環(huán)的變形量極其微小，對(duì)軸承油膜厚度影響可以忽略不計(jì)，將浮環(huán)看作剛體可以滿足要求；而當(dāng)轉(zhuǎn)速較高、偏心率大于0.6時(shí)，必須考慮浮環(huán)變形量對(duì)油膜厚度的影響?？嫡佥x等[11]分析了浮環(huán)彈性變形對(duì)浮環(huán)軸承的穩(wěn)定性和內(nèi)外油膜壓力分布的影響，指出隨著偏心率的增大，擠壓油膜系統(tǒng)的最大正油膜壓力出現(xiàn)的位置向氣穴油膜系統(tǒng)的正負(fù)壓區(qū)交界處移動(dòng)。

綜上所述，學(xué)者們對(duì)浮環(huán)軸承的熱效應(yīng)和浮環(huán)彈性變形等方面進(jìn)行了深入研究，但并沒有考慮進(jìn)油溫度以及浮環(huán)受熱膨脹變形對(duì)浮環(huán)軸承潤(rùn)滑靜特性的影響。本文作者建立了計(jì)入熱效應(yīng)的浮環(huán)軸承流體動(dòng)壓潤(rùn)滑模型，利用數(shù)值差分法聯(lián)立求解雷諾方程，同時(shí)求解了浮環(huán)軸承的端泄流量，將計(jì)算結(jié)果與文獻(xiàn)[6]中所公布的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比，驗(yàn)證了理論模型的正確性。并在此基礎(chǔ)上，考慮浮環(huán)受熱膨脹變形的影響，研究了不同供油溫度對(duì)浮環(huán)軸承的潤(rùn)滑靜特性的影響。

1 浮環(huán)軸承潤(rùn)滑理論

圖1描述了浮環(huán)軸承中元件的坐標(biāo)關(guān)系，可以看出軸頸、浮環(huán)和軸瓦三者的中心不重合。軸頸與浮環(huán)、軸瓦與浮環(huán)之間存在楔形間隙，且三者存在相對(duì)運(yùn)動(dòng)，順著軸頸轉(zhuǎn)動(dòng)方向，楔形間隙由大變小的區(qū)域稱為收斂區(qū)，由小變大的區(qū)域稱為發(fā)散區(qū)。表征浮環(huán)軸承內(nèi)外油膜壓力pi、po的雷諾方程分別為

(1)

(2)

式中：0≤x≤2πRi,o，0≤z≤Li,o;ωJ和ωR分別表示軸頸和浮環(huán)的角速度；ηi和ηo分別表示內(nèi)、外油膜潤(rùn)滑油黏度；RJ和RR分別表示軸頸和浮環(huán)半徑；hi和ho分別表示內(nèi)層和外層油膜厚度：

hi=ci+eicosθi

(3)

ho=co+eocosθo

(4)

其中：ei和eo分別表示軸頸中心相對(duì)于浮環(huán)中心、浮環(huán)中心相對(duì)于軸瓦中心偏心距；ci和co分別表示內(nèi)、外油膜間隙。

圖1 浮環(huán)軸承的徑向截面結(jié)構(gòu)示意圖Fig 1 Schematic diagram of radial cross section of floating ring bearing

1.1 浮環(huán)軸承油膜熱平衡方程

內(nèi)外油膜的摩擦功耗可由式(5)計(jì)算得出。摩擦功耗所產(chǎn)生的熱量將以對(duì)流、傳導(dǎo)和輻射3種方式與外界進(jìn)行熱交換。對(duì)于浮環(huán)軸承來(lái)說(shuō)，輻射散熱量很小，幾乎可以忽略不計(jì)，熱量主要以對(duì)流和傳導(dǎo)的方式傳播[11]。

(5)

式中：Ψt、Ψο為內(nèi)、外油膜摩擦功耗;L為軸承長(zhǎng)度；μ為潤(rùn)滑油動(dòng)力黏度；Ro表示浮動(dòng)環(huán)外徑。

1.2 浮環(huán)熱變形方程

式(6)所示為能量平衡方程。內(nèi)油膜摩擦功耗的一部分通過(guò)內(nèi)膜端泄以對(duì)流方式帶走，另一部分以傳導(dǎo)方式與軸頸和浮動(dòng)環(huán)進(jìn)行熱交換。外油膜摩擦功耗的一部分通過(guò)外膜端泄以對(duì)流方式帶走，另一部分以傳導(dǎo)方式與浮動(dòng)環(huán)和軸瓦進(jìn)行熱交換。對(duì)于浮動(dòng)環(huán)本身，其吸收的熱量等于擴(kuò)散的熱量。

(6)

式中：QSIDEi、QSIDEo為內(nèi)、外油膜端泄流量；RB、Ri為軸瓦半徑、浮動(dòng)環(huán)內(nèi)徑；λJ、λR、λB為潤(rùn)滑油與軸頸、浮動(dòng)環(huán)、軸瓦之間的熱傳導(dǎo)率；ΔTi、ΔTo、ΔTJ、ΔTR、ΔTB為內(nèi)、外油膜、軸頸、浮動(dòng)環(huán)、軸瓦的溫升；cp為潤(rùn)滑油比熱容；ρ為潤(rùn)滑油密度。

式(6)中的溫升值均為未知量，無(wú)法直接通過(guò)能量平衡方程求出，須補(bǔ)充其他約束條件。根據(jù)熱力學(xué)原理，軸承傳導(dǎo)散熱與對(duì)流散熱之比取決于材料的熱擴(kuò)散率與摩擦副的接觸情況，由式(7)表示。

(7)

式中：ki、ko為內(nèi)、外油膜傳導(dǎo)對(duì)流散熱比；α為熱擴(kuò)散率。

將式(7)與能量平衡方程(6)聯(lián)立，即可求出ΔTi、ΔTo、ΔTJ、ΔTR、ΔTB等5個(gè)部分的溫升。根據(jù)其中固體部分溫升值與材料的熱膨脹系數(shù)，即可計(jì)算浮環(huán)軸承新的內(nèi)外間隙。

(8)

式中：γJ、γR、γB為軸頸、浮動(dòng)環(huán)、軸瓦的熱膨脹數(shù)；cio、coo為浮環(huán)軸承原內(nèi)、外間隙。

浮環(huán)軸承的端泄溫升取決于內(nèi)外間隙的溫升和端泄流量，由式(9)表示。

ΔTSIDE=(ΔTiQSIDEi+ΔToQSIDEo)/(QSIDEi+QSIDEo)

(9)

式中：ΔTSIDE為端泄溫升。

受端泄溫升的影響，浮環(huán)軸承間隙內(nèi)潤(rùn)滑油的溫黏關(guān)系不能被忽略。由于精確計(jì)算油膜中的溫度分布十分復(fù)雜，文中采用有效黏度法推導(dǎo)潤(rùn)滑油黏度[12]。該方法認(rèn)為軸承中在入口溫度和出口溫度之間存在一個(gè)有效溫度，如式(10)所示。通過(guò)Reynolds溫黏關(guān)系模型，式(11)，計(jì)算該溫度下的有效黏度，作為計(jì)算軸承潤(rùn)滑性能的依據(jù)。

Tef=Tin+0.8(Tout-Tin)

(10)

式中：Tef、Tin、Tout分別為有效溫度、入口溫度和出口溫度。

μ=μ′e-β(T-T′)

(11)

式中：T′為參考溫度；μ′為潤(rùn)滑油在參考溫度下的動(dòng)力黏度；β為溫黏指數(shù)。

2 數(shù)值計(jì)算方法

采用數(shù)值差分求解雷諾方程，得到內(nèi)外油膜壓力，進(jìn)而求得內(nèi)外油膜摩擦功耗，聯(lián)立求解方程(6—11)，得到浮環(huán)軸承內(nèi)、外間隙和端泄溫升。其中壓力收斂判據(jù)為

(12)

求解流程圖如圖2所示。

2 計(jì)算程序流程圖Fig 2 Flow chart of the calculation program

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

3.1 理論模型驗(yàn)證

為了驗(yàn)證理論模型的正確性，文中根據(jù)文獻(xiàn)[6]中所公布的浮環(huán)軸承結(jié)構(gòu)和潤(rùn)滑參數(shù)，對(duì)不同軸頸轉(zhuǎn)速下的浮動(dòng)環(huán)轉(zhuǎn)速以及端泄溫升進(jìn)行了仿真計(jì)算，并將所得結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比分析。表1列出了所用的主要仿真參數(shù)。

表1 主要仿真參數(shù)Table 1 Main simulation parameters

圖3所示為文中擬合及文獻(xiàn)實(shí)驗(yàn)[6]得到的浮環(huán)軸承環(huán)速比在1 000～12 000 rad/s軸頸角速度范圍內(nèi)的變化情況。可見，文獻(xiàn)[6]實(shí)驗(yàn)結(jié)果與擬合結(jié)果之間的相對(duì)誤差僅為1.8%～5.1%。

圖3 不同轉(zhuǎn)速下浮環(huán)軸承環(huán)速比Fig 3 The speed ratio of floating ring and journal under different speeds

采用圖3給出的環(huán)速比，通過(guò)計(jì)算得到浮環(huán)軸承端泄流量溫升在1 000～12 000 rad/s軸頸角速度范圍內(nèi)的變化情況，并與文獻(xiàn)[6]實(shí)驗(yàn)結(jié)果進(jìn)行比較，如圖4所示。文獻(xiàn)[6]實(shí)驗(yàn)結(jié)果與仿真計(jì)算結(jié)果之間的相對(duì)誤差僅為1.3%～4.2%，表明文中理論模型計(jì)算出的流量溫升趨勢(shì)基本反應(yīng)了實(shí)際端泄流量溫升的變化。

圖4 不同轉(zhuǎn)速下的端泄流量溫升Fig 4 Temperature rise of side leakage under different speeds

經(jīng)過(guò)與文獻(xiàn)[6]中所公布實(shí)驗(yàn)數(shù)據(jù)的對(duì)比，驗(yàn)證了所建理論模型的正確性。在此基礎(chǔ)上,進(jìn)一步研究了供油溫度對(duì)浮環(huán)軸承潤(rùn)滑靜特性參數(shù)的影響。通過(guò)曲線擬合得到浮環(huán)軸承環(huán)速比與軸徑角速度的關(guān)系如同圖3所示，其關(guān)系式為

Y=3.166X-0.2682

(13)

3.2 供油溫度對(duì)浮環(huán)軸承潤(rùn)滑靜特性的影響

本節(jié)討論供油溫度的改變對(duì)浮環(huán)軸承的摩擦功耗、端泄溫升、內(nèi)外間隙、浮動(dòng)環(huán)轉(zhuǎn)速與端泄流量的影響。

圖5、6所示為浮環(huán)軸承摩擦功耗與供油溫度之間的關(guān)系。

圖5 供油溫度對(duì)外油膜摩擦功耗的影響Fig 5 Influences of oil supply temperature on friction loss of outer oil film

圖6 供油溫度對(duì)內(nèi)油膜摩擦功耗的影響Fig 6 Influences of oil supply temperature on friction loss of inner oil film

由圖5、6可知：內(nèi)外油膜摩擦功耗隨著供油溫度的提高而降低；受溫黏關(guān)系的影響，在較低的溫度范圍內(nèi)，潤(rùn)滑油的黏度隨溫度升高而急劇下降，而在較高的溫度范圍內(nèi)，潤(rùn)滑油黏度隨溫度的變化變得平緩；另外，供油溫度的升高，使內(nèi)外油膜摩擦功耗之比呈現(xiàn)下降趨勢(shì)。這是因?yàn)殡S著供油溫度的升高，潤(rùn)滑油黏度降低，當(dāng)軸頸旋轉(zhuǎn)時(shí)，潤(rùn)滑油受到的剪切力減小，因此產(chǎn)生的摩擦功耗也逐漸降低。

圖7所示為供油溫度對(duì)浮環(huán)軸承端泄溫升的影響?？梢姡憾诵箿厣S軸頸轉(zhuǎn)速的升高而增大，而溫升隨軸頸轉(zhuǎn)速變化的幅度在采用較低溫度的潤(rùn)滑油時(shí)表現(xiàn)得更為劇烈；當(dāng)軸頸角速度達(dá)到12 000 rad/s的情況下，供油溫度為10和30 ℃時(shí)，端泄溫升分別為25和17 ℃；繼續(xù)提高供油溫度至50和70 ℃時(shí)，端泄溫升隨軸頸轉(zhuǎn)速的變化幅度相對(duì)較小，尤其在5 000 rad/s以下的軸頸角速度范圍，溫升的區(qū)別不再明顯。同樣的規(guī)律也體現(xiàn)在了軸承間隙的變化上。

圖7 供油溫度對(duì)端泄溫升的影響Fig 7 Influence of oil supply temperature on the temperature rise of side leakage

圖8、9所示為供油溫度導(dǎo)致的內(nèi)外間隙的變化情況。可見，采用較低溫度的油液潤(rùn)滑時(shí)，其內(nèi)外間隙隨軸頸轉(zhuǎn)速的變化率相對(duì)較高。由于與內(nèi)間隙相比，浮環(huán)軸承的外間隙一般較大，因此外間隙的變化受熱效應(yīng)的影響更為明顯。

圖8 供油溫度對(duì)內(nèi)油膜間隙變化率的影響Fig 8 Influences of oil supply temperature on change rate of inner clearance

圖9 供油溫度對(duì)外間隙變化率的影響Fig 9 Influences of oil supply temperature on change rate of outer clearance

圖10所示為供油溫度對(duì)浮動(dòng)環(huán)轉(zhuǎn)速的影響?？梢钥闯觯焊…h(huán)與軸頸角速度之比隨著供油溫度的升高而增大，當(dāng)隨著供油溫度繼續(xù)升高，環(huán)速比的增長(zhǎng)趨勢(shì)變緩。這是因?yàn)楣┯蜏囟忍岣?，?rùn)滑油黏度降低，摩擦阻力減小，因而浮動(dòng)環(huán)轉(zhuǎn)速越快，但隨著供油溫度繼續(xù)提高，潤(rùn)滑油黏度降低變緩，且浮環(huán)受到的剪切力也降低，則環(huán)速比變化也變緩。

圖11、12所示為浮環(huán)軸承在不同供油溫度下的端泄流量的變化規(guī)律。在采用較高的溫度供油時(shí)，潤(rùn)滑油的黏度降低，內(nèi)外油膜的端泄流量將得到提升。由圖11可見，當(dāng)采用較低溫度的油液潤(rùn)滑，在較高的軸頸轉(zhuǎn)速下，內(nèi)油膜端泄流量下降明顯，有可能造成在內(nèi)間隙內(nèi)難以形成完整的油膜，而導(dǎo)致貧油潤(rùn)滑現(xiàn)象的發(fā)生。這里需要指出的是，文中只考慮了軸承本身的摩擦生熱，沒有計(jì)入其他熱源以及潤(rùn)滑油溫度對(duì)軸承冷卻性能的影響，因此不能得出供油溫度越高，潤(rùn)滑狀況越好的結(jié)論。

圖10 供油溫度對(duì)浮動(dòng)環(huán)轉(zhuǎn)速的影響Fig 10 Influences of oil supply temperature on the speed of floating ring

圖11 供油溫度對(duì)內(nèi)油膜端泄流量的影響Fig 11 Influences of oil supply temperature on side leakage of inner oil film

圖12 供油溫度對(duì)外油膜端泄流量的影響Fig 12 Influences of oil supply temperature on side leakage of outer oil film

4 結(jié)論

(1)隨著軸頸轉(zhuǎn)速的提升，浮環(huán)軸承的摩擦功耗、端泄溫升會(huì)增加；浮環(huán)軸承內(nèi)間隙增加，外間隙減小，浮動(dòng)環(huán)與軸頸角速度之比呈下降趨勢(shì)；內(nèi)油膜端泄流量在整個(gè)軸頸轉(zhuǎn)速的范圍呈現(xiàn)先上升后下降的趨勢(shì)，而外油膜端泄流量受軸頸轉(zhuǎn)速的影響不大。

(2)隨著供油溫度的提高，浮環(huán)軸承的摩擦功耗、端泄溫升都減??；浮環(huán)軸承內(nèi)間隙減小，外間隙增加，浮動(dòng)環(huán)與軸頸角速度之比呈上升趨勢(shì)；內(nèi)外油膜端泄流量都增加。

(3)采用過(guò)低溫度的油液潤(rùn)滑時(shí)，隨著轉(zhuǎn)速的增加，外油膜端泄流量變化不大，而內(nèi)油膜端泄流量明顯下降，有可能導(dǎo)致貧油潤(rùn)滑現(xiàn)象的發(fā)生。