田奇勇 黃文俊 饒 杰 朱曉農(nóng) 王賢震

（1.合肥通用機(jī)械研究院；2.安徽安風(fēng)風(fēng)機(jī)有限公司）

0 引言

流體動(dòng)壓滑動(dòng)軸承承載能力強(qiáng)、運(yùn)行平穩(wěn)、抗振性好、壽命長(zhǎng)，廣泛應(yīng)用于壓縮機(jī)、汽輪機(jī)、離心機(jī)等設(shè)備。根據(jù)流體動(dòng)壓潤滑理論，影響滑動(dòng)軸承油膜壓力分布的參數(shù)主要有軸承寬徑比、相對(duì)間隙、徑向載荷、潤滑油工作溫度和主軸轉(zhuǎn)速等。傳統(tǒng)的軸承參數(shù)計(jì)算方法是利用已有的寬徑比的軸承靜動(dòng)態(tài)參數(shù)，首先假定平均工作溫度，按所選潤滑油確定粘度，并根據(jù)已知的軸承尺寸和運(yùn)轉(zhuǎn)參數(shù)，計(jì)算出軸承承載量系數(shù)，再結(jié)合軸承寬徑比按已知的圖表插值求得軸承工作狀態(tài)下的各參數(shù)值，這種方法不太方便且有較大誤差。

本文采用有限差分法對(duì)液體動(dòng)壓滑動(dòng)軸承的靜態(tài)參數(shù)進(jìn)行計(jì)算，編制出通用性較強(qiáng)的實(shí)用計(jì)算程序，以代替?zhèn)鹘y(tǒng)的圖表計(jì)算，縮短了設(shè)計(jì)周期，提高了設(shè)計(jì)效率，對(duì)流體動(dòng)壓滑動(dòng)軸承的設(shè)計(jì)及故障診斷均具有較大幫助。

1 Reynolds基本方程

流體動(dòng)壓滑動(dòng)軸承結(jié)構(gòu)示意圖見圖1，二維穩(wěn)態(tài)等溫層流不可壓縮流體Reynolds方程[6]：

式中，h為徑向軸承油膜厚度；μ為潤滑油工作溫度下動(dòng)力粘度；p為油膜壓力；U為軸承處軸徑線速度。

圖1 動(dòng)壓軸承結(jié)構(gòu)示意圖Fig.1 Dynamic pressure bearing structure diagram

研究動(dòng)壓徑向滑動(dòng)軸承的關(guān)鍵問題就是對(duì)Reynolds方程進(jìn)行求解，由于二維Reynolds方程的復(fù)雜性，一般都進(jìn)行各種假設(shè)以將其簡(jiǎn)化。比較典型的是將軸承假設(shè)為無限寬或無限窄，此時(shí)方程轉(zhuǎn)化成一維穩(wěn)態(tài)流動(dòng)雷諾方程，目前對(duì)這兩種情況下的雷諾方程可以得到較為精確的解析解，對(duì)研究徑向滑動(dòng)軸承的各種性能有一定的指導(dǎo)意義，如這些解析解指出了軸承性能的趨向，而且確定了軸承性能的上限和下限，然而實(shí)際軸承是有限寬的，需求解完整的二維流動(dòng)雷諾方程[7]。

為了便于計(jì)算機(jī)運(yùn)算，將Reynolds方程無量綱化，取軸承處軸半徑r作為“相對(duì)單位”來度量x，得φ=x r，φ為從豎直位置開始計(jì)算的軸承偏角，其取值范圍根據(jù)具體軸瓦結(jié)構(gòu)來確定；對(duì)軸向坐標(biāo)z，選取軸承寬度L的一半作為相對(duì)單位，則z方向的無量綱坐標(biāo)λ=z(L/2)，坐標(biāo)原點(diǎn)放在寬度中央，則λ的區(qū)間是-1≤λ≤1。無量綱油膜厚度H=h c=1+εcosφ，c為軸承半徑間隙，ε為偏心率，則得到無量綱化形式為：

式中,D為軸徑直徑；p=p p0為無量綱化油膜壓力，p0=2Ω μ ψ2;ψ為半徑間隙比。

2 數(shù)值求解方法

2.1 差分法求解Reynolds方程原理

采用有限差分法將式（2）離散，將軸瓦的油膜劃分為m×n個(gè)方格（m為圓周φ方向，n為寬度λ方向），將網(wǎng)格節(jié)點(diǎn)按所在列數(shù)和行數(shù)順序編號(hào)，沿φ方向的列數(shù)用i編號(hào)，即從1到m+1；沿λ方向的行數(shù)用j編號(hào)，即從1到n+1；每個(gè)節(jié)點(diǎn)的位置用(i,j)二維編號(hào)表示，用各個(gè)節(jié)點(diǎn)上的壓力值構(gòu)成各階差商，近似取代Reynolds方程中的導(dǎo)數(shù)，則近似方程為：

將式（3）差分表達(dá)式代入式（2）中得：

從式（3）和式（4）將Reynolds方程化為一組代數(shù)方程，根據(jù)(i,j)節(jié)點(diǎn)周圍四節(jié)點(diǎn)上的壓力值來計(jì)算中間節(jié)點(diǎn)的壓力值，由此解出各節(jié)點(diǎn)上的壓力值，得到一組離散的壓力數(shù)值，近似表達(dá)出油膜中的壓力分布，然后根據(jù)這組壓力值，用相應(yīng)的數(shù)值積分，求得軸承的承載力、阻力、流量、溫升等靜態(tài)性能值。

2.2 求解邊界條件

求解Reynolds方程時(shí)，邊界條件的處理對(duì)壓力分布的影響很大，因此，邊界條件的確定是一個(gè)十分重要的問題。常用的邊界條件有Sommerfeld邊界條件，半Sommerfeld邊界條件和Reynolds邊界條件，對(duì)于有限長(zhǎng)動(dòng)壓軸承，公認(rèn)Reynolds邊界條件比較符合實(shí)際[5-7]。Reynolds邊界條件認(rèn)為油膜不連續(xù)，壓力油膜的終點(diǎn)在最小油膜厚度后油楔發(fā)散區(qū)的某個(gè)位置φp處，該處油膜出現(xiàn)自然破裂，本文計(jì)算采用Reynolds邊界條件。

2.3 程序化計(jì)算流程與分析

傳統(tǒng)方法是利用已有的軸承寬徑比下的靜動(dòng)態(tài)參數(shù)圖表，首先假定軸承的工作平均溫度，按所選潤滑油確定油的粘度，并根據(jù)已知的軸承尺寸和運(yùn)轉(zhuǎn)參數(shù)，計(jì)算出軸承承載量系數(shù)，再由和軸承寬徑比L D按已知的圖表插值或擬合求得軸承工作狀態(tài)下的各參數(shù)值，其目的在于核算工作狀態(tài)下最小油膜厚度、功耗、供油量和溫升等是否在容許范圍內(nèi)。這種插值或擬合會(huì)引入誤差，影響結(jié)果的判斷，給設(shè)計(jì)帶來諸多不便。

本文編程計(jì)算流程見圖2，直接利用已知的運(yùn)轉(zhuǎn)參數(shù)對(duì)外載荷進(jìn)行迭代計(jì)算，避免了插值運(yùn)算帶來的誤差，求解精度較高。

圖2 計(jì)算流程圖Fig.2 Calculation flow chart

3 算例及結(jié)果分析

以文獻(xiàn)[6]第五章表5-4中寬徑比L D=1.0的圓柱軸承和表5-5中L D=0.5及橢圓比為0.5的橢圓軸承為例。圓柱軸承各參數(shù)計(jì)算結(jié)果見表1，偏位角θ的最大偏差不超過0.30%，承載量系數(shù)的偏差最大不超過0.1%；橢圓軸承各參數(shù)計(jì)算結(jié)果見表2，偏位角θ的最大偏差不超過0.30%，承載量系數(shù)的偏差最大不超過0.6%。本文的計(jì)算結(jié)果與文獻(xiàn)[6]的結(jié)果吻合良好，從而驗(yàn)證了本文所述方法的計(jì)算程序的精度較高。

文獻(xiàn)[6]第五章例1（部分計(jì)算結(jié)果見表1），某機(jī)械用流體動(dòng)壓徑向滑動(dòng)軸承，軸承處載荷F=31 882.5N，轉(zhuǎn)速n=3 600r/min，軸承公稱直徑D=152mm，寬徑比L D=1.0，相對(duì)間隙ψ=0.002，潤滑油工作粘度為18.006×10-3Pa·s。

表1 圓柱軸承各參數(shù)計(jì)算結(jié)果對(duì)比[6]Tab.1 The comparison of θandmcalculation results in cylindrical bearing

表1 圓柱軸承各參數(shù)計(jì)算結(jié)果對(duì)比[6]Tab.1 The comparison of θandmcalculation results in cylindrical bearing

偏心率ε θ>m本文文獻(xiàn)[6]本文文獻(xiàn)[6]0.3 55.791 6 55.857 7 0.676 3 0.675 7 0.4 50.057 1 50.116 6 1.025 3 1.025 4 0.5 45.559 5 45.442 3 1.520 7 1.520 6 0.6 40.791 4 40.774 5 2.293 9 2.296 0 0.7 35.894 5 35.878 6 3.666 7 3.663 8

表2 橢圓軸承各參數(shù)計(jì)算結(jié)果對(duì)比[6]Tab.2 The comparison of θandF>mcalculation results in elliptical bearing

表2 橢圓軸承各參數(shù)計(jì)算結(jié)果對(duì)比[6]Tab.2 The comparison of θandF>mcalculation results in elliptical bearing

偏心率ε θm本文文獻(xiàn)[6]本文文獻(xiàn)[6]0.3 89.630 5 89.850 4 0.072 7 0.073 1 0.4 88.585 7 88.580 9 0.105 0 0.105 6 0.5 87.070 1 87.023 9 0.146 3 0.145 5 0.6 84.601 3 84.683 3 0.201 1 0.202 3 0.7 81.551 9 81.232 1 0.285 6 0.286 6

按傳統(tǒng)計(jì)算方法，根據(jù)已知條件計(jì)算出承載量系數(shù)=1.626，按表1插值計(jì)算得到軸承偏心率ε=0.513 7，偏位角θ=44.806 9°，然后再插值計(jì)算出其它性能參數(shù)，這樣會(huì)帶來插值計(jì)算誤差。本文方法程序計(jì)算得到軸承偏心率ε=0.516 8，偏位角θ=44.778 6°，同時(shí)可獲得最小油膜厚度、油膜破裂位置和摩擦功耗等參數(shù)。軸承下瓦無量綱油膜厚度和無量綱油膜壓力分布分別見圖3和圖4所示。

圖3 圓柱軸承無量綱油膜厚度分布Fig.3 Non-dimensional of cylindrical bearing oil film thickness distribution

圖4 圓柱軸承無量綱油膜壓力分布Fig.4 Non-dimensionless oil film pressure distribution of cylindrical bearing

通過圖4可知，沿周向無量綱油膜壓力先增大后減小，在某一點(diǎn)達(dá)到最大值，并且呈非線性；沿寬度方向以中央最大呈拋物線狀。當(dāng)?shù)玫接湍毫?，可以通過對(duì)壓力分布積分求得流量系數(shù)、阻力系數(shù)、摩擦功耗、油膜剛度和油膜阻尼系數(shù)等[8-10]。

4 結(jié)論

采用有限差分法求解二維Reynolds方程可以計(jì)算出軸承油膜厚度和油膜壓力分布，同時(shí)可直接得到軸承偏心率、偏位角、最小油膜厚度、油膜破裂位置和摩擦阻力等參數(shù)，即采用數(shù)值計(jì)算可以方便、快捷地得到不同參數(shù)下軸承的各項(xiàng)性能值，且軸承參數(shù)可以任選，避免了計(jì)算傳統(tǒng)方法中根據(jù)圖表的插值或擬合所產(chǎn)生的誤差。本方法的使用還可以較好地掌握軸承的實(shí)際工作狀態(tài)，對(duì)實(shí)際設(shè)備的運(yùn)行維護(hù)、軸承的選型設(shè)計(jì)有較強(qiáng)的指導(dǎo)意義。

[1]黃民毅.液體動(dòng)壓滑動(dòng)軸承設(shè)計(jì)的數(shù)值計(jì)算[J].四川工業(yè)學(xué)院學(xué)報(bào)，1998，17（3）:37-41.

[2]張祖立.液體動(dòng)壓徑向滑動(dòng)軸承靜態(tài)特性參數(shù)的計(jì)算機(jī)程序設(shè)計(jì)與計(jì)算[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào)，1999，30（4）:440-443.

[3]孟繁娟，杜永平.徑向滑動(dòng)軸承油膜壓力分析[J].軸承，2008（1）：23-25.

[4]王增勝，穆國華，牛月蘭，等.基于Matlab的動(dòng)壓軸承油膜特性仿真研究[J].機(jī)械，2013，40（2）：9-10，23

[5]謝帆，荊建平，萬召，等.基于有限差分法的徑向滑動(dòng)軸承油膜壓力分布計(jì)算[J].潤滑與密封，2012，37（2）：12-15.

[6]西安交通大學(xué)透平壓縮機(jī)教研室編著.離心式壓縮機(jī)強(qiáng)度[M].北京：機(jī)械工業(yè)出版社，1980.

[7]敏政，王樂，魏志國，等.基于MATLAB技術(shù)的滑動(dòng)軸承油膜壓力分布的模擬[J].潤滑與密封，2008，33（8）：51-53.

[8]康寧，遲啟明.基于N-S方程的徑向滑動(dòng)軸承油膜承載力的計(jì)算[J].潤滑與密封，2010，35（10）：10-11.

[9]李元生，敖良波，李磊，等.滑動(dòng)軸承動(dòng)力特性系數(shù)動(dòng)態(tài)分析方法[J].機(jī)械工程學(xué)報(bào)，2010，46（21）：48-53.

[10]李強(qiáng)，許偉偉，王振波，等.滑動(dòng)軸承動(dòng)力特性的數(shù)值計(jì)算方法[J].中國石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，38（5）：165-172.