曹春建，辛華榮，張德虎，吳忠

(1.河海大學(xué)能源與電氣學(xué)院，江蘇南京 211100;2.江蘇省走馬塘張家港樞紐工程建設(shè)處，江蘇張家港 215600)

導(dǎo)軸承是保證水泵穩(wěn)定運(yùn)行的重要部件，起著承受水泵徑向力、控制主軸擺度的作用。和油潤滑導(dǎo)軸承相比，水潤滑導(dǎo)軸承因密封簡單、安全節(jié)能等優(yōu)點(diǎn)，逐漸應(yīng)用于在大型臥式貫流泵中［1-2］。由于貫流泵轉(zhuǎn)速低、載荷高，軸承潤滑水膜不易形成，使得軸承承載區(qū)磨損嚴(yán)重，直接影響到水泵運(yùn)行的可靠性。為解決上述問題，除了選擇適合的軸承材料外，對(duì)軸承參數(shù)合理的選擇也至關(guān)重要。目前國內(nèi)外學(xué)者對(duì)滑動(dòng)軸承的研究大多采用有限差分法、有限元法等數(shù)值方法對(duì)簡化的雷諾方程進(jìn)行求解［3-5］。這些方法忽略了慣性項(xiàng)、徑向流場變化等因素的影響，當(dāng)軸承形狀較為復(fù)雜時(shí)，其應(yīng)用范圍受到了限制。隨著計(jì)算流體力學(xué)(CFD)的發(fā)展，基于有限體積法的三維CFD技術(shù)以其簡便、精確的優(yōu)點(diǎn)正逐漸應(yīng)用于滑動(dòng)軸承參數(shù)優(yōu)化中［6-8］。本文針對(duì)江蘇省某大型貫流泵水潤滑導(dǎo)軸承，利用Fluent軟件對(duì)不同轉(zhuǎn)速下軸承壓力場進(jìn)行仿真分析，并得出相關(guān)參數(shù)對(duì)軸承潤滑水量及其靜態(tài)性能的影響規(guī)律，為水泵導(dǎo)軸承的參數(shù)優(yōu)化提供一種新的方法和更全面的參考。

1 導(dǎo)軸承模型的建立

1.1 導(dǎo)軸承三維模型及網(wǎng)格劃分

水潤滑導(dǎo)軸承的結(jié)構(gòu)如圖1所示。軸承主要由外部軸承體和樹脂苯酚導(dǎo)軸瓦組成，導(dǎo)軸瓦為徑向整體式結(jié)構(gòu)。供水系統(tǒng)采用恒壓式，清水通過金屬軟管壓入軸承供水端，然后從另一端泄漏，其結(jié)構(gòu)參數(shù)如表1所示。

圖1 水潤滑導(dǎo)軸承結(jié)構(gòu)示意圖

表1 導(dǎo)軸承設(shè)計(jì)參數(shù)

采用Fluent前處理軟件Gambit建立導(dǎo)軸承三維幾何模型并進(jìn)行網(wǎng)格劃分。因模型最小間隙僅為0.05 mm，若直接對(duì)模型進(jìn)行體網(wǎng)格劃分，極易造成網(wǎng)格扭曲率大以及最小體積為負(fù)的情形，影響到計(jì)算收斂速度和結(jié)果的準(zhǔn)確性。為得到高品質(zhì)的網(wǎng)格，本文分別對(duì)模型上下半周、軸向及徑向設(shè)置不同數(shù)量節(jié)點(diǎn)?？紤]到軸承壁面流場變量梯度較大，需對(duì)其附近精細(xì)網(wǎng)格，故對(duì)徑向節(jié)點(diǎn)進(jìn)行加密處理，之后對(duì)面和體依次進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分，最終得到六面體網(wǎng)格(圖2)。

圖2 導(dǎo)軸承網(wǎng)格結(jié)構(gòu)圖

為研究網(wǎng)格數(shù)對(duì)仿真結(jié)果的影響，本文分別以網(wǎng)格數(shù)40.8、81.6和120萬對(duì)飛逸轉(zhuǎn)速下軸承進(jìn)行仿真，所得結(jié)果如表2所示。由表可知，網(wǎng)格數(shù)到達(dá)一定值時(shí)，軸承的承載力、最大壓力及潤滑水量差異較小，綜合考慮計(jì)算時(shí)間和精度，網(wǎng)格數(shù)確定為81.6萬。

表2 不同網(wǎng)格數(shù)下仿真結(jié)果

1.2 模型假設(shè)及邊界條件

a)模型假設(shè)如下:1)潤滑介質(zhì)為牛頓液體，黏度沿膜厚方向無變化，忽略介質(zhì)慣性力;2)潤滑介質(zhì)與壁面無相對(duì)滑移和熱量傳遞，不考慮黏性生成熱;3)大氣壓力為101 kPa，考慮介質(zhì)重力。

b)邊界條件設(shè)置如下:1)供水邊界設(shè)為壓力進(jìn)口，壓力設(shè)置為供水壓力值;2)出水邊界設(shè)為壓力出口，出口壓力等于大氣壓力;3)水膜內(nèi)壁設(shè)為旋轉(zhuǎn)邊界，采用黏性無滑移邊界條件，即水膜內(nèi)壁速度等于軸頸速度;4)水膜外壁設(shè)為固定邊界且無滑移。

1.3 控制方程

軸承處于動(dòng)壓潤滑時(shí)，在軸承發(fā)散楔形區(qū)的水膜因不能承受負(fù)壓而發(fā)生破裂，從而產(chǎn)生空化現(xiàn)象。目前對(duì)動(dòng)壓軸承水膜空化邊界的確定常采用三種條件，即Half-Sommerfeld，Reynolds和JFO［10］?？紤]到 Reynolds和 JFO 邊界條件采用CFD仿真較為困難，故本文采用相對(duì)簡單的Half-Sommerfeld空化邊界條件。采用Half-Sommerfeld邊界條件的不可壓縮定常流動(dòng)控制方程如下［8］:

待水膜壓力求解出之后，利用Fluent的UDF接口將負(fù)壓置零。

2 數(shù)值模擬及分析

文獻(xiàn)［11］認(rèn)為滑動(dòng)軸承內(nèi)部流場雷諾數(shù)Re到1 000后便完全進(jìn)入湍流狀態(tài)，本文引用其觀點(diǎn)，經(jīng)計(jì)算，軸頸的臨界轉(zhuǎn)速為152.8 r/min。當(dāng)軸頸轉(zhuǎn)速低于臨界轉(zhuǎn)速時(shí)采用層流模型(laminar)進(jìn)行計(jì)算，而當(dāng)高于臨界轉(zhuǎn)速則采用RNGk-ε湍流模型;選擇壓力基隱式求解器，近壁面處理采用標(biāo)準(zhǔn)壁面函數(shù)法，壓力-速度耦合用SIMPLE算法，控制方程的離散格式為二階迎風(fēng)格式，模擬精度設(shè)為1×10-4，初始化流場后進(jìn)行迭代計(jì)算。

2.1 水膜壓力場分析

貫流泵轉(zhuǎn)速較低且隨著工況改變而發(fā)生相應(yīng)變化，為研究轉(zhuǎn)速對(duì)軸承潤滑性能的影響，在原設(shè)計(jì)參數(shù)下，分別對(duì)30%低轉(zhuǎn)速(46.4 r/min)、額定轉(zhuǎn)速(154.9 r/min)及飛逸轉(zhuǎn)速(249.8 r/min)進(jìn)行數(shù)值模擬，得到壓力分布情況如圖3:

圖3 不同轉(zhuǎn)速下水膜壓力場分布

由圖3可知，低轉(zhuǎn)速下軸承壓力分布較為均勻，軸向壓力從軸承進(jìn)水端到出水端逐漸降低，隨著轉(zhuǎn)速升高，高低壓集中區(qū)開始出現(xiàn)在軸承最小間隙的兩側(cè)，并分別從供水端和泄水端向另一端逐漸延伸。這是由于軸頸偏心使軸承和軸頸之間形成楔形間隙，當(dāng)軸頸圍繞軸心以較高速度運(yùn)轉(zhuǎn)時(shí)，水膜會(huì)在楔形收斂區(qū)受到軸頸擠壓作用而形成高壓區(qū)，在楔形發(fā)散區(qū)則因水膜速度迅速增大而形成低壓區(qū)，從而形成水膜動(dòng)壓效應(yīng)?？梢钥闯觯~定轉(zhuǎn)速下軸承高壓區(qū)的水膜只覆蓋了承載區(qū)的一半，使得主軸和軸承會(huì)有部分直接接觸，從而加劇軸承承載區(qū)的磨損，可見額定工況下軸承潤滑性能并不理想。當(dāng)水泵處于飛逸轉(zhuǎn)速時(shí)，軸承承載區(qū)的高壓水膜較為完整，從而使得主軸和軸承能完全隔開，此時(shí)軸承潤滑性能較好。綜上可知，提高水泵額定轉(zhuǎn)速有助于增強(qiáng)軸承潤滑性能。

不同轉(zhuǎn)速下軸承中心截面(Z=180 mm)及截面角θ為40°的壓力分布曲線如圖4所示。

圖4 水膜壓力分布曲線圖

從圖4(a)可知，隨著轉(zhuǎn)速升高，以x=0 mm為界，水膜高壓區(qū)的周向壓力逐漸增大而低壓區(qū)則逐漸降低，壓力隨轉(zhuǎn)速升高變化很明顯;從圖4(b)可知，低轉(zhuǎn)速、額定轉(zhuǎn)速下軸承軸向壓力呈線性變化，此時(shí)水在軸承中的流動(dòng)以供水端到泄水端的軸向流動(dòng)為主，滿足設(shè)計(jì)要求。然而飛逸轉(zhuǎn)速下，軸承軸向壓力呈拋物線變化，Z=1 000處的壓力達(dá)到0.82 MPa，遠(yuǎn)高于供水壓力0.5 MPa，這種情況會(huì)使?jié)櫥畯膲毫Ψ逯档奈恢孟蚬┧说沽?，從而降低軸承潤滑水量。實(shí)際上在軸承低壓集中區(qū)，水膜因不能承受低壓而發(fā)生破裂，以致形成空穴，此時(shí)水量降低則很有可能讓水不能及時(shí)充滿空穴而生成大量氣泡，氣泡破裂會(huì)產(chǎn)生巨大破壞力而損傷軸承和軸頸表面。由此可見，額定轉(zhuǎn)速下軸承采用一端供水另一端泄水的方式是合適的，而當(dāng)軸承動(dòng)壓效應(yīng)顯著時(shí)，建議采用中間供水兩端泄水的方式，從而降低空化形成的概率。

2.2 相關(guān)參數(shù)對(duì)軸承靜態(tài)性能的影響

圖5 相關(guān)參數(shù)與承載力的關(guān)系

軸承的剛度及承載力在工程應(yīng)用中是很重要的兩個(gè)參數(shù)［12］，其反映的是軸承靜態(tài)性能。圖5和圖6反映的是其他參數(shù)不變情況下，偏心率為0.6、0.7、0.8、0.9時(shí)，承載力、剛度與轉(zhuǎn)速、供水壓力及半徑間隙之間的關(guān)系曲線。如圖所示，同一偏心率下，軸承承載力及剛度隨轉(zhuǎn)速、供水壓力增大呈非線性增長，隨半徑間隙增大呈非線性減小;其他參數(shù)不變時(shí)，軸承承載力隨偏心率的增大而增大;當(dāng)轉(zhuǎn)速低于額定轉(zhuǎn)速，半徑間隙大于0.3 mm時(shí)，軸承剛度則與偏心率無關(guān)。

2.3 相關(guān)參數(shù)對(duì)軸承潤滑水量的影響

潤滑水量是否充足是保證滑動(dòng)軸承穩(wěn)定運(yùn)行、避免軸承溫度過高的重要前提條件。圖7反映的是其他參數(shù)不變情況下，偏心率為0.6、0.7、0.8、0.9時(shí)，潤滑水量與轉(zhuǎn)速、供水壓力及半徑間隙之間的關(guān)系。由圖7(a)可知，偏心率小于等于0.7時(shí)，潤滑水量與轉(zhuǎn)速無關(guān)，當(dāng)偏心率大于0.7時(shí)，潤滑水量隨轉(zhuǎn)速升高而呈非線性降低;同一轉(zhuǎn)速下，潤滑水量隨著偏心率的增大而增大。由圖7(b)和(c)可知，同一偏心率下，潤滑水量隨供水壓力、半徑間隙的增大而呈非線性增長，可以看出，當(dāng)半徑間隙低于0.2 mm時(shí)，潤滑水量與偏心率無關(guān)。

圖7 相關(guān)參數(shù)與潤滑水量的關(guān)系

3 結(jié)論

1)通過CFD數(shù)值模擬發(fā)現(xiàn)不同轉(zhuǎn)速下軸承水膜壓力分布有較大差異，即隨著轉(zhuǎn)速的升高，高低壓區(qū)的水膜分別從供水端和泄水端向另一端延伸。當(dāng)軸承動(dòng)壓效應(yīng)顯著時(shí)，潤滑水會(huì)倒流至供水端，以致潤滑水量降低;

2)同一偏心率下，軸承承載力與剛度隨轉(zhuǎn)速、供水壓力增大呈非線性增長，隨半徑間隙增大呈非線性降低;潤滑水量隨供水壓力和半徑間隙的增大呈非線性增長;低偏心率下，潤滑水量和轉(zhuǎn)速無關(guān)，高偏心率下，潤滑水量隨轉(zhuǎn)速升高而減小;

3)利用三維CFD技術(shù)可以很好模擬導(dǎo)軸承內(nèi)部流場特性，所得出的結(jié)論為水潤滑導(dǎo)軸承的參數(shù)優(yōu)化提供一定的參考價(jià)值。

［1］彭晉民，朱志宏，王家序.水潤滑軸承的研究現(xiàn)狀及展望［J］.潤滑與密封，2004，(5):124-130.

［2］仇寶云，魏強(qiáng)林，林海江，等.大型水泵導(dǎo)軸承應(yīng)用研究［J］.流體機(jī)械，2006，34(11):12-15.

［3］蘇葒，王小靜，張值明.滑動(dòng)軸承兩種油膜邊界條件的比較［J］.潤滑與密封，2002，(5):2-4.

［4］呂延軍.基于Reynolds邊界的滑動(dòng)軸承動(dòng)力學(xué)系數(shù)的計(jì)算及應(yīng)用［J］.摩擦學(xué)學(xué)報(bào)，2004，24(1):61-65.

［5］Lo Chengying，Wang Chengchi，Lee Yuhan.Performance Analysis of High-speed Spindle Aerostatic Bearings［J］.Tribology International，2005，38:5-14.

［6］Guo Zenglin，Hirano Toshio，Gordon Kirk R.Application of CFD analysis for rotating machinery-Part I:hydrodynamic，hydrostatic bearings and squeeze film damper［J］.Journal of Engineering for Gas Turbines and Power，2005，127(2):445-451.

［7］彭婭玲，張志國，陳汝剛，等.CFD輔助船舶艉部水潤滑軸承設(shè)計(jì)的研究［J］.潤滑與密封，2008，33(5):72-76.

［8］葉曉琰，張軍輝，江偉，等.海水淡化泵水潤滑軸承間隙的優(yōu)化設(shè)計(jì)［J］，排灌機(jī)械工程學(xué)報(bào)，2010，28(2):117-121.

［9］張楚，楊建剛，郭瑞，等.基于兩相流理論的滑動(dòng)軸承流場計(jì)算分析［J］.中國電機(jī)工程學(xué)報(bào)，2010，30(29):80-84.

［10］宋寅，李雪松，顧春偉.一種軸承潤滑油新空化模型［J］.清華大學(xué)學(xué)報(bào)(自然科學(xué)版)，2010，50(7):1047-1052.

［11］段芳莉.橡膠軸承水潤滑機(jī)理的研究［D］.重慶:重慶大學(xué)，2002，3.

［12］張杰，郭宏升，牛D，等.水潤滑高速動(dòng)靜壓滑動(dòng)軸承數(shù)值模擬［J］.農(nóng)業(yè)機(jī)械學(xué)報(bào)，2008，39(6):159-162.