表面織構(gòu)排布方式對徑向滑動軸承摩擦性能的影響*

王保民 趙瑞平 房文博 閆瑞祥 劉洪芹

(蘭州理工大學(xué)機(jī)電工程學(xué)院 甘肅蘭州 730050)

由于其結(jié)構(gòu)簡單、接觸面積大、摩擦學(xué)性能好，滑動軸承在旋轉(zhuǎn)機(jī)械中得到了廣泛的應(yīng)用[1]。然而，隨著載荷和轉(zhuǎn)速的增加，摩擦和磨損是導(dǎo)致滑動軸承失效的主要原因[2]。因此，提高承載力和減小摩擦一直是滑動軸承設(shè)計和研究的主要內(nèi)容。表面織構(gòu)技術(shù)通過在摩擦副表面制備出具有一定尺寸和排列的凹坑、凹痕或凸包等圖案的點陣，能有效改變潤滑表面之間的摩擦學(xué)性能，已成為改善滑動軸承摩擦性能的一種有效手段[3]。

HAMILTON等[4]于1966年最先提出表面織構(gòu)的概念，引起國內(nèi)外學(xué)者的極大關(guān)注。SUH等[5]通過實驗發(fā)現(xiàn)，在真空及無潤滑條件下，溝槽表面和光滑表面相比，摩擦因數(shù)低并且穩(wěn)定；同時還揭示了表面溝槽結(jié)構(gòu)能夠捕捉和容納磨屑，從而防止產(chǎn)生三體磨損的減磨機(jī)制。尹明虎等[6]對比分析了圓形、方形、三角形微織構(gòu)對徑向滑動軸承性能的影響，指出設(shè)計合理的微織構(gòu)能夠有效提高滑動軸承的潤滑性能。LIN等[7-9]研究了考慮空化效應(yīng)的大面積織構(gòu)界面滑移現(xiàn)象，結(jié)果表明在織構(gòu)結(jié)構(gòu)設(shè)計中，新型織構(gòu)及滑移邊界條件可以作為傳統(tǒng)邊界條件的補(bǔ)充。鄭文斌等[10]研究發(fā)現(xiàn)，排布合理的織構(gòu)能夠顯著提高軸承剛度以及轉(zhuǎn)子系統(tǒng)穩(wěn)定性。彭龍龍等[11]數(shù)值分析了矩形、菱形、球形和Koch雪花形等4種不同形狀的表面織構(gòu)對軸承湍流動能的影響，提出了一種遞增多尺度織構(gòu)，其承載能力優(yōu)于遞減多尺度織構(gòu)。MARIAN等[12]基于熱流體動力學(xué)模型分析了帶有方形凹坑的止推軸承潤滑油膜厚度與承載能力的關(guān)系，并對織構(gòu)幾何參數(shù)進(jìn)行了優(yōu)化。蘇華和余志雄[13]研究發(fā)現(xiàn)，在織構(gòu)型滑動軸承的基礎(chǔ)上增加溝槽可以改變油膜分布，降低端泄量。馮慧慧等[14]通過實驗表明，階梯槽織構(gòu)能有效降低水潤滑止推軸承摩擦力。紀(jì)敬虎等[15]分析了在不同收斂比參數(shù)下，局部微凹坑織構(gòu)對軸瓦潤滑油膜壓力分布的影響，并指出收斂比與織構(gòu)參數(shù)共同影響軸承承載能力，在進(jìn)行織構(gòu)參數(shù)設(shè)計時需要考慮收斂比的影響。王麗麗等[16]通過數(shù)值仿真和實驗，分析了微織構(gòu)尺寸對軸承摩擦磨損性能的影響，結(jié)果表明合適的微織構(gòu)尺寸可以有效提高滑動軸承的潤滑性能、降低摩擦因數(shù)及減小磨損。于如飛和陳渭[17]基于Winkler彈性基礎(chǔ)模型，計算了軸承表面的變形量，發(fā)現(xiàn)與軸承表面的彈性變形相比，表面織構(gòu)對徑向軸承靜特性的影響更加顯著，且織構(gòu)可以明顯改變油膜壓力的分布情況。趙遠(yuǎn)方等[18]利用激光打標(biāo)機(jī)加工不同排列方式凹坑織構(gòu)，用自主搭建的摩擦磨損實驗臺測量摩擦副表面力矩，結(jié)果表明無織構(gòu)表面磨損遠(yuǎn)大于微織構(gòu)表面，且軸向間隔分布織構(gòu)減磨效果最好。羅敏敏等[19]通過銷盤實驗?zāi)M滑動軸承織構(gòu)化摩擦副表面工況，分析了軸徑和軸套表面橢圓凹坑織構(gòu)對鉆頭滑動軸承摩擦學(xué)性能的影響，結(jié)果顯示橢圓凹坑同時交錯分布于銷盤試樣表面時潤滑減磨效果最優(yōu)，橢圓凹坑同時平行分布于銷盤試樣表面時潤滑減磨效果最差。

以上研究為滑動軸承表面織構(gòu)分析和設(shè)計奠定了重要的理論依據(jù)和試驗基礎(chǔ)，然而，現(xiàn)有研究主要探究了織構(gòu)形狀、結(jié)構(gòu)參數(shù)等特征對滑動軸承摩擦性能的影響及機(jī)制，關(guān)于織構(gòu)排布位置及不同位置織構(gòu)排布方式對軸承摩擦性能影響的研究較少。鑒于此，本文作者以HZ040滑動軸承為研究對象，采用MRF方法和Mixture模型建立了表面織構(gòu)化滑動軸承數(shù)值分析模型，運(yùn)用壓力耦合方程的改進(jìn)半隱式SIMPLEC(Semi-Implicit Method for Pressure Linked Equations Consistent)算法進(jìn)行求解，對比分析了表面織構(gòu)的軸向和周向排布方式以及表面形狀等在不同位置對軸承承載力與摩擦因數(shù)的影響，設(shè)計滑動軸承復(fù)合型織構(gòu)排布方式，對比分析了復(fù)合型織構(gòu)對滑動軸承摩擦性能的提升效果。

1 控制方程

采用基于Navier-Stokes(N-S)方程的CFD方法對織構(gòu)化徑向滑動軸承摩擦性能進(jìn)行分析，N-S方程表達(dá)式為

(1)

(2)

式中：t是時間；u是潤滑油速度；F是每單位質(zhì)量流體的質(zhì)量力；p是潤滑油壓力；ρ是潤滑油密度。

滑動軸承油膜厚度方程為

h=h1+h2

(3)

h1=c+ecosθ

(4)

(5)

式中：e為軸頸偏心距；c為軸承間隙；h3為織構(gòu)處油膜厚度。

采用壓力耦合方程的改進(jìn)半隱式SIMPLEC算法進(jìn)行求解，得到軸承油膜的壓力和切應(yīng)力。由壓力對軸承內(nèi)壁面面積的積分求得軸承承載力

Fs=?pdxdy

(6)

由切應(yīng)力對面積的積分求得軸承油膜的摩擦力

Ff=?τdxdy

(7)

摩擦因數(shù)即為摩擦力與承載力之比

f=Ff/Fs

(8)

2 數(shù)值計算模型

2.1 幾何模型

文中以HZ040徑向滑動軸承為研究對象，軸承具體結(jié)構(gòu)參數(shù)見表1。由于軸承在實際工況中施加載荷后有偏心的存在，流體域由寬變窄、再由窄變寬，具有一定黏度的潤滑油流入楔形間隙后由于楔形效應(yīng)在最小油膜厚度兩側(cè)形成對稱分布的高低壓區(qū)域，如圖1所示。因此，定義流體域由寬變窄的一側(cè)高壓流體域為“收斂楔”區(qū)域，流體域由窄變寬的另一側(cè)負(fù)壓流體域為“發(fā)散楔”區(qū)域[20]。

圖1 油膜壓力分布Fig.1 Pressure distribution of the oil film

表1 滑動軸承結(jié)構(gòu)參數(shù)Table 1 Structural parameters of the journal bearing

將油入口定義為起始0°，沿內(nèi)壁面順時針方向，以最大與最小油膜厚度為分界線，從50° (最大油膜厚度位置)起，每隔60°劃分一個織構(gòu)區(qū)域，將收斂楔區(qū)域均分為三部分，即收斂楔入口(50°～110°)、收斂楔中部(110°～170°)、收斂楔出口(170°～230°)，230°(最小油膜厚度位置)至290°為發(fā)散楔區(qū)域，在軸承腔內(nèi)共劃分出4個區(qū)域，如圖2、3所示。

圖2 滑動軸承腔結(jié)構(gòu)示意Fig.2 Structure of journal bearing cavity

圖3 滑動軸承流體域展開示意Fig.3 Expansion of fluid domain inside journal bearing

2.2 網(wǎng)格劃分

由于軸承腔內(nèi)潤滑油最小油膜厚度僅為16 μm，而軸承寬度為40 mm，寬徑比過大，為保證網(wǎng)格不發(fā)生畸變，保持運(yùn)算穩(wěn)定性，全局采用六面體結(jié)構(gòu)網(wǎng)格，如圖4所示。

圖4 流體域網(wǎng)格模型Fig.4 The mesh model of bearing fluid domain

2.3 邊界條件及求解器設(shè)置

位于軸承頂部的潤滑油入口采用壓力入口，進(jìn)油壓力設(shè)置為0.2 MPa。軸承腔體兩端為潤滑油出口，與外界大氣相連，潤滑油出口采用壓力出口，設(shè)置壓力為大氣壓。軸承腔內(nèi)表面設(shè)置為旋轉(zhuǎn)壁面，以模擬軸徑的旋轉(zhuǎn)運(yùn)動，外壁面設(shè)置為固定壁面，軸承轉(zhuǎn)速設(shè)置為5 000 r/min。油氣材料參數(shù)見表2。

表2 油氣材料參數(shù)Table 2 Material parameters of oil and air

采用Mixture模型模擬油氣二相流，因網(wǎng)格長寬比較大，為利于計算結(jié)果的收斂性，采用壓力耦合方程的改進(jìn)半隱式SIMPLEC算法進(jìn)行求解。為保證計算精確性，動量方程采用二階迎風(fēng)格式進(jìn)行離散，體積分?jǐn)?shù)采用QUICK (Quadratic Upwind Interpolation of Convective Kinematics)格式離散。殘差值設(shè)置為1×10-4，同時監(jiān)測出入口油相體積變化率，變化率低于1×10-5，兩判據(jù)均滿足則判定結(jié)果收斂。在求得油膜壓力、切應(yīng)力、速度分布后，通過積分計算可以獲得軸承的承載力和摩擦力以及摩擦因數(shù)等特性參數(shù)。

2.4 網(wǎng)格無關(guān)性驗證及模型驗證

運(yùn)用ICEM CFD對流體域模型進(jìn)行網(wǎng)格劃分，為了在保證精確性前提下盡可能提高計算速度，需進(jìn)行網(wǎng)格無關(guān)性驗證，分別于徑向、周向與軸向進(jìn)行網(wǎng)格層數(shù)加密，設(shè)置7組對照。如表3所示，當(dāng)徑向網(wǎng)格層數(shù)大于7時計算誤差小于1%，因此，最終選定后續(xù)計算中網(wǎng)格數(shù)均大于901 440。

表3 不同網(wǎng)格數(shù)目下最大油膜壓力Table 3 Maximum oil film pressure of different number of grids

通過與文獻(xiàn)[16]中的仿真和實驗數(shù)據(jù)進(jìn)行比較，來驗證文中計算方法的有效性。表4給出文獻(xiàn)[16]中的軸承參數(shù)，軸承模型如圖5所示。表5、圖6所示為2種方法計算結(jié)果。由表5可見，文中方法計算的最大壓力和承載能力文獻(xiàn)[16]計算結(jié)果相差很小。由圖6也可以看出，仿真結(jié)果與文獻(xiàn)[16]中的實驗數(shù)據(jù)基本一致。這表明文中提出的滑動軸承計算方法有較高的可靠性。

圖5 文獻(xiàn)[16]中的軸承結(jié)構(gòu)Fig.5 Bearing structure in literature 16

圖6 文中計算結(jié)果與文獻(xiàn)[16]中實驗結(jié)果比較Fig.6 Comparison of the calculation results in this paper and the experimental results in literature 16

表4 文獻(xiàn)[16]中軸承的結(jié)構(gòu)參數(shù)和邊界條件Table 4 Structural parameters and boundary conditions of bearings in literature 16

表5 文中計算結(jié)果與文獻(xiàn)[16]對比Table 5 Comparison of calculation results in this paper and in the literature 16

3 結(jié)果與分析

3.1 織構(gòu)軸向和周向排布對軸承性能的影響

3.1.1 織構(gòu)軸向排布對軸承性能的影響

為分析織構(gòu)軸向排列對滑動軸承承載能力和摩擦因數(shù)的影響，在保證織構(gòu)軸向總寬度相等的情況下，設(shè)計了4種不同的矩形織構(gòu)，如圖7所示，分別排布在軸向位置Ⅰ(收斂楔入口)、Ⅱ(收斂楔中部)、Ⅲ(收斂楔出口)、Ⅳ(發(fā)散楔)?？棙?gòu)分布周向間距為9.5°，織構(gòu)結(jié)構(gòu)如圖8所示，具體參數(shù)及排布見表6。軸承油膜承載力與摩擦因數(shù)如圖9所示，可以看出：

圖7 矩形織構(gòu)軸向排列示意Fig.7 Schematic of the axial arrangement of rectangular textures：(a)groove texture；(b)short groove texture；(c)rectangular texture；(d)square texture

圖8 矩形織構(gòu)剖面圖Fig.8 Sectional drawing of rectangular texture

圖9 織構(gòu)軸向排布對軸承摩擦特性的影響Fig.9 Effect of texture axial arrangement on friction properties of journal bearing：(a)bearing capacity；(b)friction coefficient

表6 不同微織構(gòu)參數(shù)及軸向排布Table 6 The parameters and axial arrangement of rectangular textures

(1)在位置 Ⅰ (收斂楔入口)軸向排布矩形織構(gòu)時，油膜承載力明顯提高，摩擦因數(shù)略低于無織構(gòu)軸承，方形織構(gòu)和短溝槽織構(gòu)效果更好。這是因為收斂楔入口處的壓力相對較低，當(dāng)潤滑油從織構(gòu)出口流出時，軸承間隙突然減小，壓力增大。潤滑油在表面織構(gòu)中的擴(kuò)散和收斂引起的動壓效應(yīng)起主導(dǎo)作用。LI等[21]指出，表面織構(gòu)可以提高承載力的一個主要原因是潤滑劑流入表面凹槽內(nèi)部形成擴(kuò)散和收斂從而引起動壓效果。

(2)當(dāng)矩形織構(gòu)軸向排布在位置 Ⅱ (收斂楔中部)時，油膜壓力隨著油膜厚度的減小而增大。在宏觀上滑動軸承的動壓效應(yīng)以及微觀上織構(gòu)處的動壓效應(yīng)的共同作用下，排布4種織構(gòu)表現(xiàn)出不同的摩擦性能，短溝槽織構(gòu)降低承載能力，增加摩擦因數(shù)，矩形織構(gòu)和方形織構(gòu)同時增加承載能力和摩擦因數(shù)，溝槽織構(gòu)增加了承載能力并略微降低了摩擦因數(shù)。

(3)當(dāng)矩形織構(gòu)軸向排布在位置Ⅲ(收斂楔出口)時，承載能力下降，摩擦因數(shù)大幅增加。這是因為收斂楔的出口處于油膜厚度最薄的高壓區(qū)，導(dǎo)致油膜承載能力下降，摩擦因數(shù)增加。

(4)當(dāng)矩形織構(gòu)軸向排列在位置Ⅳ(發(fā)散楔)時，溝槽和短溝槽織構(gòu)降低了承載能力和摩擦因數(shù)。矩形織構(gòu)降低承載能力，增加了摩擦因數(shù)，而方形紋理增加了承載能力，降低了摩擦因數(shù)。

3.1.2 織構(gòu)周向排布對軸承性能的影響

如圖10所示，在保證織構(gòu)周向總寬度相等的情況下，分別在周向位置Ⅰ、Ⅱ、Ⅲ、Ⅳ設(shè)計4種不同的矩形織構(gòu)。具體織構(gòu)參數(shù)見表7，織構(gòu)化滑動軸承的承載能力和摩擦因數(shù)如圖11所示。

圖10 矩形織構(gòu)周向排列示意Fig.10 Schematic of the circumferential arrangement of rectangular textures：(a)groove texture；(b)short groove texture；(c)rectangular texture；(d)square texture

圖11 織構(gòu)周向排布對軸承摩擦特性的影響Fig.11 Effect of texture circumferential arrangement on friction properties of journal bearing：(a)bearing capacity：(b)friction coefficient

表7 不同微織構(gòu)參數(shù)及周向排布Table 7 The parameters and circumferential arrangement of rectangular textures

由圖11可見，在位置 Ⅰ (收斂楔入口)排布周向織構(gòu)，潤滑油動壓效應(yīng)使軸承承載力明顯提高，而對摩擦因數(shù)影響很??；在位置 Ⅱ (收斂楔中部)排布周向織構(gòu)使軸承承載力明顯提高，溝槽和短溝槽織構(gòu)使摩擦因數(shù)降低，而矩形和方形織構(gòu)使摩擦因數(shù)增大；在位置 Ⅲ(收斂楔出口)處排布周向織構(gòu)，對軸承承載力影響較小，而使摩擦因數(shù)明顯增加，這說明在高壓區(qū)不適合布置周向織構(gòu)，尤其是方形等“密布型”織構(gòu)；在位置 Ⅳ (發(fā)散楔)周向排布短溝槽和方形織構(gòu)使軸承承載力小幅增加、摩擦因數(shù)略有降低，而溝槽和矩形織構(gòu)使軸承承載力小幅降低、摩擦因數(shù)略有增加。

3.2 織構(gòu)形狀對軸承性能的影響

為分析織構(gòu)形狀對軸承承載力與摩擦因數(shù)的影響，對表8所示7種形狀織構(gòu)進(jìn)行對比分析。7種形狀織構(gòu)排布的面積率一致，每個微織構(gòu)單元面積一致，軸向分布20排，周向分布7列，織構(gòu)深度均為20 μm。

表8 不同形狀的表面織構(gòu)參數(shù)Table 8 Parameters of surface texture with different shape

在位置Ⅰ、Ⅱ、Ⅲ、Ⅳ分別排布這7種不同形狀織構(gòu)時，軸承油膜承載力與摩擦因數(shù)如圖12所示?？梢钥闯觯?/p>

圖12 織構(gòu)形狀和周向排布對軸承摩擦特性的影響Fig.12 Effect of texture shape and circumferential arrangement on friction properties of journal bearing：(a)bearing capacity；(b)friction coefficient

(1)當(dāng)7種織構(gòu)排布在位置 Ⅰ (收斂楔入口)時，油膜的承載能力明顯提高。7種織構(gòu)中，橢圓、倒橢圓、三角形和倒三角形織構(gòu)略微增加摩擦因數(shù)，菱形和圓形織構(gòu)略微降低摩擦因數(shù)，而十字形織構(gòu)可以顯著提高承載能力，降低摩擦因數(shù)。

(2)當(dāng)織構(gòu)排布在位置 Ⅱ (收斂楔中部)時，7種織構(gòu)中，菱形、倒橢圓、十字形織構(gòu)可以提高承載能力，降低摩擦因數(shù)，而其他形狀織構(gòu)則增加了承載能力和摩擦因數(shù)。

(3)當(dāng)7種織構(gòu)排布在位置 Ⅲ (收斂楔出口)時，承載力下降，摩擦因數(shù)顯著增加，其中倒三角形織構(gòu)最為明顯。

(4)當(dāng)織構(gòu)排布在位置 Ⅳ (發(fā)散楔)時，7種織構(gòu)中，圓形、倒三角形和十字織構(gòu)增加了承載能力和摩擦因數(shù)，橢圓和倒橢圓織構(gòu)降低了承載能力和摩擦因數(shù)，菱形和三角形織構(gòu)增加了承載能力并降低了摩擦因數(shù)。相比較而言，在位置 Ⅳ 處布置三角形織構(gòu)更有利于改善滑動軸承的摩擦性能。

3.3 新型復(fù)合型織構(gòu)對軸承性能的影響

綜合以上分析可以得出，位置 Ⅰ (收斂楔入口)布置十字形織構(gòu)，位置 Ⅱ (收斂楔中部)布置菱形、倒橢圓形、周向溝槽和短溝槽、十字形織構(gòu)，位置 Ⅳ (發(fā)散楔)布置三角形織構(gòu)，可以獲得較好的軸承摩擦性能，而位置 Ⅲ (收斂楔出口)不宜布置織構(gòu)。根據(jù)以上分析，設(shè)計了5種復(fù)合織構(gòu)，如表9所示，其徑向深度為20 μm。復(fù)合織構(gòu)化滑動軸承的承載能力和摩擦因數(shù)如圖13所示。

圖13 不同復(fù)合型織構(gòu)軸承的摩擦特性Fig.13 Friction properties of journal bearing with different compound textures

表9 復(fù)合型織構(gòu)參數(shù)Table 9 Parameters of the composite textures

由圖13可見，與無織構(gòu)滑動軸承相比，A、C、D、E型復(fù)合織構(gòu)分別提高了承載能力8.33%、15.69%、15.66%、17.16%，分別降低了摩擦因數(shù)1.25%、3.75%、4.17%、2.92%；而B型復(fù)合織構(gòu)承載力增加了10.78%，同時摩擦因數(shù)增加了1.46%?？梢缘贸鼋Y(jié)論，C、D、E型復(fù)合織構(gòu)可以顯著提高滑動軸承的摩擦性能。

4 結(jié)論

以HZ040徑向滑動軸承為研究對象，基于N-S方程，采用CFD方法分析表面織構(gòu)排布對軸承摩擦性能的影響。主要結(jié)論如下：

(1)受流體動壓效應(yīng)的影響，在位置 Ⅰ (收斂楔入口)軸向、周向排布不同形狀的織構(gòu)，油膜承載能力顯著提高，但減摩效果不明顯。其中排布十字形織構(gòu)時可以顯著提高承載能力。

(2)在位置 Ⅱ (收斂楔中部)周向排布矩形織構(gòu)，增強(qiáng)了流體域內(nèi)潤滑油的擠壓強(qiáng)度，潤滑油的動壓效應(yīng)更加明顯，因此更有利于改善滑動軸承的摩擦性能，尤其是溝槽和短溝槽織構(gòu)。與其他形狀織構(gòu)相比，菱形、倒橢圓形和十字形織構(gòu)可以更明顯地提高滑動軸承的摩擦性能。

(3)在位置 Ⅲ (收斂楔出口)排布織構(gòu)不利于軸承摩擦性能的提高，因此不適宜排布織構(gòu)；在位置 Ⅳ (發(fā)散楔)布置三角形織構(gòu)有利于改善軸承的摩擦學(xué)性能。

(4)與單一織構(gòu)相比，合理的復(fù)合織構(gòu)設(shè)計能更有效地改善軸承的摩擦性能，在五種復(fù)合織構(gòu)中，當(dāng)在位置 Ⅰ (收斂楔入口)布置十字形織構(gòu)，位置 Ⅱ (收斂楔中部)布置周向溝槽、短溝槽和十字形織構(gòu)，位置 Ⅳ (發(fā)散楔)布置三角形織構(gòu)能顯著改善摩擦性能。