張艷芹+孔祥濱+郭麗麗+程海闊+李銳

摘 要：針對重型裝備制造業(yè)中重型靜壓軸承承載特性的研究，考慮到不同工況下間隙油膜厚度對靜壓軸承承載能力及壓力分布的影響，建立了靜壓軸承間隙油膜三維模型及邊界條件，利用CFD（computational fluid dynamic）原理，應(yīng)用動網(wǎng)格技術(shù)和FLUENT軟件，探討靜壓軸承轉(zhuǎn)速為10r/min以及在空載0t、有載40t、滿載150t不同工況下，油膜厚度變化對壓力場以及油腔壓力值的影響規(guī)律。結(jié)果表明：油腔壓力隨著間隙油膜厚度的減小而增大，當(dāng)油膜減小到一定值時，油腔壓力顯著增加，油膜承載能力顯著增強。

關(guān)鍵詞：靜壓軸承；CFD；動網(wǎng)格；承載特性

DOI：10.15938/j.jhust.2017.06.004

中圖分類號： TH133.3

文獻標(biāo)志碼： A

文章編號： 1007-2683（2017）06-0020-04

Abstract：Research on the bearing characteristics of heavy hydrostatic bearing in heavy equipment manufacturing industry， considering the influence of film thickness on the bearing capacity and pressure distribution of hydrostatic bearing under different working conditions， 3D model and boundary condition of hydrostatic bearing oil film are established. Discussing the influence of oil film thickness variation on pressure field and oil chamber pressure under hydrostatic bearing running in the noload 0t， load 40t， full load 150t and the rotating speed of 10r/min， by using CFD （computational fluid dynamic） principle， dynamic mesh technology and FLUENT software. The results show that the oil cavity pressure increases with the decrease of oil film thickness， when the oil film decreases to a certain value， the pressure and bearing characteristics of oil cavity increases significantly.

Keywords：hydrostatic bearing； CFD； dynamic mesh； loadbearing characteristics

0 引 言

本文研究的重型靜壓軸承用于重型數(shù)控裝備制造業(yè)中，靜壓軸承的原理早在19世紀(jì)就被發(fā)現(xiàn)，但直到20世紀(jì)50年代才在發(fā)達國家興盛起來。我國從50年代后期開始液體靜壓軸承的應(yīng)用研究工作，60年代初開始在金屬切削機床上推廣應(yīng)用[1]。液

體靜壓軸承因其能提供廣泛的油膜厚度以及較高的油膜剛度，具有功率損耗小，低轉(zhuǎn)速下工作平穩(wěn)等特點被廣泛應(yīng)用到重型裝備制造業(yè)中并成為核心部件。

近年來，對于靜壓軸承的研究，國外學(xué)者考慮到靜壓軸承結(jié)構(gòu)參數(shù)變化對軸承承載性能的影響進行了理論與實驗探究。2002年，SATISH對不同油腔形狀的靜壓軸承的靜態(tài)和動態(tài)性能進行了理論分析，探究了其對靜壓軸承承載性能的影響[2]。GRABOVSKII在負(fù)載和轉(zhuǎn)速恒定的條件下運用變積分的方法研究了氣體靜壓推力軸承最大承載能力下的最佳間隙[3]。文[4]通過求解雷諾方程，并用有限元法計算了軸承間隙中速度場和壓力場，研究了節(jié)流器的尺寸和油腔的幾何形狀對靜壓軸承的流態(tài)及壓力分布的影響。文[5]研究了流體在穩(wěn)流和湍流狀態(tài)時，對液體靜壓軸承動靜態(tài)參數(shù)的影響，涉及到剛度和阻尼特性的研究。

國內(nèi)學(xué)者從事靜壓軸承技術(shù)的研究開始于20世紀(jì)50年代后期，文[6-11]對橢圓油腔和扇形油腔的靜壓軸承靜止?fàn)顟B(tài)時的流場進行仿真，指出扇形油腔結(jié)構(gòu)的靜壓軸承要優(yōu)于橢圓油腔。2006年劉賓等采用有限差分法，對徑向空氣軸承的壓力場進行了二維數(shù)值仿真，得到了在不同偏心率下的壓力場分布圖[12]。2009年YU X J等建立了大尺度靜壓軸承液-固態(tài)耦合數(shù)值分析模型，給出了空心同軸和鼓狀內(nèi)襯的最大應(yīng)力位置，得出由于空心同軸和鼓狀內(nèi)襯的變形，導(dǎo)致油膜最薄厚度減少了12.1%[13]。2011年哈爾濱工業(yè)大學(xué)的汪圣飛對徑向推力軸承的靜、動態(tài)性能展開了研究，并運用Fluent軟件進行計算，得到了軸承氣膜中壓強以及氣體流速的分布情況并進行了實驗驗證[14]。王幸福建立流體動力潤滑的數(shù)學(xué)模型對軸承潤滑特性進行研究，得到不同磨損狀態(tài)下的軸承油膜的潤滑特性[15]。上海大學(xué)的李松生、周鵬等利用超高轉(zhuǎn)速電主軸結(jié)構(gòu)，通過改變供油量和轉(zhuǎn)速等狀態(tài)參數(shù)，進行了軸承潤滑性能的試驗研究[16]。唐健研究了靜壓軸承節(jié)流器對軸承承載能力及油膜剛性的影響，采用四種不同的節(jié)流方式進行比較驗證[17]。

從以上國內(nèi)外學(xué)者研究綜述可以看出，對于靜壓軸承的研究主要還是側(cè)重于結(jié)構(gòu)參數(shù)變化對潤滑特性及承載性能的影響，而對于高速重載多油墊靜壓軸承在油膜厚度變化過程中的承載特性的研究較少。本文在團隊前期研究基礎(chǔ)上[18-20]，基于有限體積法和動網(wǎng)格技術(shù)對靜壓軸承間隙油膜進行了動態(tài)模擬，文中研究轉(zhuǎn)速10r/min，分別在空載、有載、滿載工況下，對不同油膜厚度下靜壓軸承的承載性能數(shù)值計算，得到了不同承載條件下油膜厚度變化對壓力場以及油腔壓力值的影響的變化規(guī)律，研究結(jié)果為重型靜壓軸承設(shè)計計算方面的提供理論參考。endprint

1 油膜模型及基本假設(shè)條件

1.1 油膜模型

靜壓軸承間隙油膜整體上采用圓環(huán)形油墊潤滑型式，圓周上共包括24個油腔，并且在圓周上呈周期對稱分布。根據(jù)工廠實際應(yīng)用的尺寸數(shù)據(jù)對重型靜壓軸承間隙油膜三維模型進行建模。靜壓軸承間隙油膜三維模型如圖1所示，其中包括油腔、封油面、回油槽、油膜、入油孔等組成部分。

1.2 基本假設(shè)及計算條件

基本假設(shè)及計算條件：①工作介質(zhì)為不可壓縮流體且流態(tài)為三維非定常流動；②流體流過進油孔，經(jīng)計算雷諾數(shù)Re<1000，故進油孔內(nèi)部流態(tài)為層流；流體進入油腔后，工作臺與導(dǎo)軌之間經(jīng)計算雷諾數(shù)Re<2300，導(dǎo)軌內(nèi)部為層流狀態(tài)；③油膜厚度變化為線性變化；④旋轉(zhuǎn)過程中，不考慮工作臺及底座的熱變形。

2 油膜數(shù)學(xué)模型

文中所研究的重型機床尺寸很大，故將多油墊扇形腔簡化成矩形單腔平面油墊結(jié)構(gòu)來計算，單個油墊結(jié)構(gòu)尺寸如圖2所示，圖中陰影部分面積為有效承載面積，其中R1、R2、R3、R4、φ1、φ2為靜壓工作臺單個油墊的結(jié)構(gòu)尺寸，L、l、B、b為簡化后矩形單腔平面油墊結(jié)構(gòu)尺寸。

通過上式可以看出，運行時靜壓軸承承載特性在采用定量方式供油時，支承的承載能力不僅和靜壓導(dǎo)軌上油腔及封油邊的幾何結(jié)構(gòu)尺寸有關(guān)，還和油膜厚度以及油液的粘度有關(guān)。

3 數(shù)值計算及邊界條件

3.1 邊界條件設(shè)置

在設(shè)置邊界條件時，計算模型為壓力基隱式求解器，流動入口設(shè)置為速度進口邊界，出口設(shè)置選用壓力出口，流動出口靜壓為大氣壓。

根據(jù)工程實際運用UG建模軟件進行三維仿真模型的創(chuàng)建，導(dǎo)軌上周期性沿圓周方向陣列著24個油墊，各個油墊中流體的流動狀態(tài)及流動特性相同，為了便于進行數(shù)值模擬分析，工作臺油膜整體被簡化成周期性油墊模型如圖3所示。

為了研究在高速條件下，靜壓工作臺間隙油膜厚度變化對壓力場的影響，選取比較有特征的截面來進行分析，為此以進油孔軸線為基準(zhǔn)進行周向AA截面，如圖4所示。

3.2 數(shù)值模擬及結(jié)果分析

重型靜壓軸承間隙油膜壓力場反映出支承的承載性能，本文基于有限體積法和動網(wǎng)格技術(shù)對靜壓軸承間隙油膜進行了動態(tài)模擬，得到了不同承載條件下油膜厚度變化對壓力場以及油腔壓力的影響規(guī)律。根據(jù)重型靜壓軸承實際工作中的承載情況，研究轉(zhuǎn)速10r/min，分別在空載、有載、滿載工況下即載荷為0t、40t、150t三種工況時，不同膜厚下靜壓軸承的承載性能。

由于篇幅有限，本文給出空載和滿載兩種載荷工況膜厚為0.12mm和0.04mm的油膜壓力場，油膜壓力場瞬態(tài)計算結(jié)果如圖5～6所示，其余工況見結(jié)果分析折線圖，壓力分布圖中壓力單位為Pa。

由計算結(jié)果可知，空載及滿載條件下，油腔內(nèi)的壓力場呈現(xiàn)均勻分布，壓力值沿著出口方向逐漸減小。

為了更清楚表達出不同載荷下，油膜承載能力與膜厚的關(guān)系，得出油膜厚度變化過程中油腔壓力變化折線圖如圖7所示。

由上述關(guān)系曲線可知，油膜厚度對軸承承載能力影響較大，隨著膜厚的減小油膜壓力逐漸增大，當(dāng)油膜厚度低于0.06mm時，隨著膜厚的減少，油膜承載能力急劇增強。

4 結(jié) 論

1）本文通過采用有限體積法及動網(wǎng)格技術(shù)對靜壓軸承進行數(shù)值模擬計算，得到了靜壓軸承內(nèi)部壓力場的分布，為探究流體內(nèi)部壓力分布提供了可靠的研究方法，提高了計算精度。

2）通過模擬靜壓軸承油膜壓力場可知，空載、有載及滿載條件下，油腔內(nèi)的壓力場呈現(xiàn)均勻分布，壓力值沿著出口方向逐漸減小。

3）通過計算不同油膜厚度下的油腔壓力值可知，油膜厚度對軸承承載能力影響較大，隨著膜厚的減小油膜壓力逐漸增大，當(dāng)油膜厚度低于0.06mm時，隨著膜厚的減少，油膜承載能力急劇增強。

參 考 文 獻：

[1] 鐘洪， 張冠坤. 液體靜壓動靜壓軸承設(shè)計使用手冊[M]. 北京：電子工業(yè)出版社， 2007.

[2] SATISH C S， JAIN S C， BHARUKA D K. Influence of Recess Shape on the Performance of a Capillary Compensated Circular Thrust pad Hydrostatic Bearing [J]. Tribology International， 2002， 35（6）： 347-356.

[3] GRABOVSKII V I. Optimum Clearance of a gas Hydrostatic Thrust Bearing with Maximum Load Capacity [J]. Fluid Dynamics， 2000， 35（4）： 525-533.

[4] CHEN G H， KANG Y， CHANG Y P， et al. Influences of Recess Geometry and Restrictor Dimension on Flow Patterns and Pressure Distribution of Hydrostatic Bearings[C]// American Society of Mechanical Engineers， Proceedings. Ann Arbor， USA， 2007（2）：1045-1053.

[5] DWIVEDI V K， CHAND S， PANDEY K N. Effect of Different Flow Regime on the Static and Dynamic Performance Parameter of Hydrodynamic Bearing[J]. Procedia Engineering， 2013（51）：520-528.endprint

[6] 邵俊鵬，張艷芹，李鵬程.基于FLUENT的靜壓軸承橢圓腔和扇形腔靜止?fàn)顟B(tài)流場仿真[J]. 潤滑與密封， 2007， 32（1）： 93-95.

[7] 邵俊鵬，張艷芹. 重型靜壓軸承扇形腔和圓形腔溫度場數(shù)值模擬與分析[J]. 水動力學(xué)研究與進展A輯， 2009，24（1）：119-123.

[8] 于曉東，周啟慧，王志強. 高速重載靜壓推力軸承溫度場速度特性[J]. 哈爾濱理工大學(xué)學(xué)報， 2014， 19（1）： 1-4.

[9] 于曉東，潘澤，何宇. 重型靜壓推力軸承間隙油膜流態(tài)的數(shù)值模擬[J]. 哈爾濱理工大學(xué)學(xué)報， 2015， 20（6）： 42-45.

[10]于曉東，高春麗，邱志新. 高速重載靜壓推力軸承潤滑性能研究[J]. 中國機械工程，2013，24（23）：3230-3234.

[11]YU Xaiodong， WANG Zhiqiang， MENG Xiuli. Research on Dynamic Pressure of Hydrostatic Thrust Bearing Under the Different Recess Depth and Rotating Velocity [J]. International Journal of Control and Automation，2014，7（2）：439-446.

[12]劉賓， 劉波. 徑向空氣軸承壓力場的數(shù)值分析[J].功能部件， 2006， 23（1）： 93-95.

[13]YU X J， WANG J X， LI Y. Study of Fluidsolid Simulation on the Largescale Hydrostatic Bearing of hollow Coaxial[C]//Materials Science Forum，2009：281-286.

[14]汪圣飛. 徑向推力聯(lián)合氣體靜壓軸承靜動態(tài)性能分析[D]. 哈爾濱：哈爾濱工業(yè)大學(xué)，2011：43-45.

[15]王幸福. 乳化液泵曲柄軸承潤滑特性研究[J]. 現(xiàn)代機械， 2010（5）：33-34.

[16]李松生， 周鵬， 黃曉， 等. 基于油氣潤滑的超高轉(zhuǎn)速電主軸軸承潤滑性能的試驗研究[J]. 潤滑與密封， 2011， 36（10）：25-28.

[17]唐健.不同節(jié)流方式的靜壓軸承承載性能研究[J]. 機床與液壓， 2010， 38（12）： 77-80.

[18]張艷芹， 陳瑤. 四種油腔形狀重型靜壓軸承承載性能理論分析[J]. 哈爾濱理工大學(xué)學(xué)報， 2013，18（2）：68-71.

[19]ZHANG Yanqin， YU Zeyang， CHEN Yao. Simulation and Experimental Study of Lubrication Characteristics of Vertical Hydrostatic Guide Rail[J].高技術(shù)通訊（英文版），2014，20（3）：315-320.

[20]ZHANG Yanqin， FAN Liguo. Simulation and Experimental Analysis on Supporting Characteristics of Multiple Oil Pad Hydrostatic Bearing Disk[J]. Journal of Hydrodynamics，2013，25（2）：236-241.

（編輯：關(guān) 毅）endprint