王樹圣,郭紅,李瑞珍,張紹林
(鄭州大學(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)行分析。
高速動(dòng)壓滑動(dòng)軸承的結(jié)構(gòu)示意圖如圖1所示。圖中,F(xiàn)為軸承載荷;φ為周向角度;O為軸承中心;Oj為軸頸中心;e為軸承與軸頸間的偏心距;ω為軸頸角速度。
圖1 軸承結(jié)構(gòu)示意圖
由于計(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)
軸承采用常用潤(rùn)滑油,其黏度隨溫度升高而降低,壓力對(duì)黏度的影響很小,可以忽略不計(jì)[11]。引入Reynolds黏溫方程[12]
μ=μ0e-α(T-T0),
(2)
式中:μ為潤(rùn)滑油黏度,Pa·s;α為黏溫指數(shù),℃-1;T為油膜溫度,℃;T0為潤(rùn)滑油初始溫度,取40 ℃。
量綱一的油膜厚度為
(3)
式中:h為油膜厚度,mm;ε為偏心率。
變密度、變黏度的量綱一的Reynolds方程為
(4)
壓力邊界條件選取Reynolds邊界條件,表述為
(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ù)。
高速動(dòng)壓滑動(dòng)軸承的結(jié)構(gòu)參數(shù)及運(yùn)行參數(shù)見表1。
表1 高速動(dòng)壓滑動(dòng)軸承的結(jié)構(gòu)參數(shù)及運(yùn)行參數(shù)
轉(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)
量綱一的承載力隨轉(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 量綱一的端泄量
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ì)量綱一的端泄量影響較大。