王樹(shù)圣，郭紅，李瑞珍，張紹林

(鄭州大學(xué) 機(jī)械工程學(xué)院，鄭州 450001)

隨著旋轉(zhuǎn)機(jī)械向高速重載方向發(fā)展，動(dòng)壓滑動(dòng)軸承的油膜溫度場(chǎng)和黏度場(chǎng)分布越來(lái)越不均勻，熱效應(yīng)將顯著影響其油膜特性[1-2]，同時(shí)又因沖擊和壓力下降等因素導(dǎo)致原本潤(rùn)滑油中溶解的氣體釋放出來(lái)或是液態(tài)潤(rùn)滑油變?yōu)橛蜌夂笮纬蓺庋ǎ治鰸?rùn)滑油膜時(shí)考慮氣穴的影響更接近實(shí)際油膜特征[3-4]，所以在研究軸承特性時(shí)綜合考慮氣穴和熱效應(yīng)影響更加符合實(shí)際工況。文獻(xiàn)[5]在理論和試驗(yàn)的基礎(chǔ)上，采用熱流連續(xù)邊界條件提出了滑動(dòng)軸承一維溫度場(chǎng)穩(wěn)態(tài)計(jì)算模型。文獻(xiàn)[6]采用耦合求解的方法分析了滑動(dòng)軸承腔內(nèi)潤(rùn)滑油的溫度分布，并進(jìn)行了流場(chǎng)等效溫度的合理確定。文獻(xiàn)[7]研究了氣穴對(duì)徑推浮環(huán)軸承壓力場(chǎng)和靜特性的影響，得出氣穴使油膜壓力峰值、承載力和摩擦功耗下降的結(jié)論。文獻(xiàn)[8]應(yīng)用ANSYS中的CFX模塊在不同偏心率、寬徑比和轉(zhuǎn)速下分析了氣穴現(xiàn)象對(duì)滑動(dòng)軸承摩擦性能的影響。文獻(xiàn)[9]研究發(fā)現(xiàn)，表面形貌和熱效應(yīng)對(duì)燃油泵滑動(dòng)軸承的潤(rùn)滑性能存在直接影響，且偏心率越大，影響越明顯。文獻(xiàn)[10]提出了氣液兩相流潤(rùn)滑時(shí)軸承性能的數(shù)值求解方法，并討論了流體性質(zhì)對(duì)軸承主要運(yùn)行參數(shù)的影響，發(fā)現(xiàn)氣穴對(duì)軸承性能有顯著影響?，F(xiàn)基于上述文獻(xiàn)建立相關(guān)數(shù)學(xué)模型，并采用有限差分法進(jìn)行求解，進(jìn)而對(duì)不同工況下計(jì)入氣穴和熱效應(yīng)影響時(shí)高速動(dòng)壓滑動(dòng)軸承的壓力分布和靜特性參數(shù)進(jìn)行分析。

1 數(shù)學(xué)模型

高速動(dòng)壓滑動(dòng)軸承的結(jié)構(gòu)示意圖如圖1所示。圖中，F(xiàn)為軸承載荷；φ為周向角度；O為軸承中心；Oj為軸頸中心；e為軸承與軸頸間的偏心距；ω為軸頸角速度。

圖1 軸承結(jié)構(gòu)示意圖

1.1 氣油兩相流當(dāng)量密度

由于計(jì)入氣穴的影響，動(dòng)壓滑動(dòng)軸承的潤(rùn)滑已變?yōu)闅庥蛢上嗔鳚?rùn)滑。為方便研究，假設(shè)油液中的氣泡很小且均勻分布，氣泡隨潤(rùn)滑油以相同速度流動(dòng)，則兩相流可視為單相擬Newton流體，其密度與油膜壓力有關(guān)。

不考慮溫度變化對(duì)密度的影響，潤(rùn)滑流體量綱一的當(dāng)量密度[7]為

(1)

1.2 黏溫關(guān)系

軸承采用常用潤(rùn)滑油，其黏度隨溫度升高而降低，壓力對(duì)黏度的影響很小，可以忽略不計(jì)[11]。引入Reynolds黏溫方程[12]

μ=μ0e-α(T-T0)，

(2)

式中：μ為潤(rùn)滑油黏度，Pa·s；α為黏溫指數(shù)，℃-1；T為油膜溫度，℃；T0為潤(rùn)滑油初始溫度，取40 ℃。

1.3 油膜厚度

量綱一的油膜厚度為

(3)

式中：h為油膜厚度，mm；ε為偏心率。

1.4 Reynolds方程及壓力邊界條件

變密度、變黏度的量綱一的Reynolds方程為

(4)

壓力邊界條件選取Reynolds邊界條件，表述為

(5)

1.5 能量方程及邊界條件

絕熱流動(dòng)假設(shè)下變密度潤(rùn)滑油膜的量綱一的能量方程為

(6)

潤(rùn)滑油膜的溫度邊界條件為φ=0,T=T0,μ=μ0。

采用有限差分法聯(lián)立求解密度-壓力方程、黏溫方程、Reynolds方程和能量方程，將壓力和溫度進(jìn)行多次迭代修正，得到油膜密度場(chǎng)、黏度場(chǎng)、溫度場(chǎng)和壓力場(chǎng)分布，進(jìn)而積分求得油膜各靜特性參數(shù)。

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

高速動(dòng)壓滑動(dòng)軸承的結(jié)構(gòu)參數(shù)及運(yùn)行參數(shù)見(jiàn)表1。

表1 高速動(dòng)壓滑動(dòng)軸承的結(jié)構(gòu)參數(shù)及運(yùn)行參數(shù)

2.1 油膜的溫度分布和壓力分布

轉(zhuǎn)速2×104r/min、含氣率0.1、偏心率0.3時(shí)的油膜溫度場(chǎng)分布如圖2所示。沿著軸頸旋轉(zhuǎn)方向，潤(rùn)滑油黏度使油質(zhì)點(diǎn)不斷消耗由軸頸供給的機(jī)械功，摩擦功耗轉(zhuǎn)化為熱量，并且在建立能量方程時(shí)假設(shè)潤(rùn)滑油絕熱流動(dòng)，故油膜溫度從初始邊界逐漸升高，在破裂邊位置達(dá)到最大值。油膜軸向溫度分布如圖3所示。由于絕熱流動(dòng)的假設(shè)和封油邊的影響，軸承兩端面處(λ=±1)的油膜溫度高于軸承中部(λ=0)的油膜溫度，且隨著周向角度的增大，二者溫差越大。由于兩端面處采用的差商格式不同導(dǎo)致軸承兩端溫度略有差異，這在圖3中φ=0.6π和φ=0.8π處顯現(xiàn)出來(lái)。

圖2 油膜溫度場(chǎng)

圖3 軸向溫度分布

轉(zhuǎn)速為2×104r/min、偏心率0.3時(shí)油膜中部(λ=0)的壓力分布及周向密度分別如圖4、圖5所示。氣穴和熱效應(yīng)均會(huì)使油膜壓力峰值下降，且隨著轉(zhuǎn)速的提高以及含氣率和偏心率的增大，而下降得越明顯；計(jì)入氣穴影響時(shí)油膜提前破裂，油膜動(dòng)壓區(qū)變小。隨著含氣率增大，量綱一的密度減小，且含氣率較高時(shí)，密度的降幅有所收窄。油膜最大壓力的變化情況如圖6所示，計(jì)入氣穴和熱效應(yīng)影響使油膜量綱一的最大壓力減小。

圖4 周向壓力分布

圖5 量綱一的周向密度

圖6 量綱一的最大壓力(n=2×104 r/min)

2.2 油膜靜特性參數(shù)

量綱一的承載力隨轉(zhuǎn)速、含氣率、偏心率的變化情況如圖7所示，在氣穴和熱效應(yīng)影響下量綱一的承載力減小，且轉(zhuǎn)速越高、偏心率和含氣率越大，減小程度越明顯。對(duì)比可知，小偏心率(ε≤0.4)下轉(zhuǎn)速和含氣率對(duì)量綱一的承載力影響不大，而在大偏心率下轉(zhuǎn)速對(duì)量綱一的承載力影響較大。

圖7 量綱一的承載力

量綱一的摩擦力變化情況如圖8所示。計(jì)入氣穴和熱效應(yīng)影響會(huì)使量綱一的摩擦力減小，但通過(guò)對(duì)比可知，計(jì)入氣穴和溫度影響后，含氣率變化對(duì)量綱一的摩擦力影響幾乎可以忽略，而轉(zhuǎn)速對(duì)量綱一的摩擦力影響相對(duì)較明顯。

圖8 量綱一的摩擦力

計(jì)入氣穴和熱效應(yīng)影響時(shí)油膜量綱一的端泄量如圖9所示。在二者影響下，量綱一的端泄量增大。轉(zhuǎn)速對(duì)量綱一的端泄量影響很小，而量綱一的端泄量增大主要是由于計(jì)入了氣穴影響。

圖9 量綱一的端泄量

3 結(jié)論

1)計(jì)入氣穴和熱效應(yīng)影響后，高速動(dòng)壓滑動(dòng)軸承油膜壓力峰值隨轉(zhuǎn)速、含氣率和偏心率的增大而下降得越明顯，且在氣穴的影響下油膜提前破裂，動(dòng)壓區(qū)變小。

2)在氣穴和熱效應(yīng)影響下，高速動(dòng)壓滑動(dòng)軸承量綱一的承載力和摩擦力減小，量綱一的端泄量增大。其中，量綱一的承載力隨轉(zhuǎn)速、含氣率和偏心率的增大而減幅越大；轉(zhuǎn)速對(duì)量綱一的承載力和摩擦力影響較大，而氣穴對(duì)量綱一的端泄量影響較大。