(1.西安科技大學工程訓練中心 陜西西安 710054;2.西安科技大學理學院 陜西西安 710054)
船舶艉軸與橡膠軸承摩擦副在運行時,軸承啟閉時出現的過載、重載及螺旋槳的重力引起的艉軸彎曲,均會引起軸與橡膠軸承的局部接觸,使得軸與軸承處于干摩擦狀態下,導致潤滑不足,磨損加劇,引起系統的振動及噪聲,造成系統的損傷事故。
國內外學者對金屬軸頸-橡膠軸承摩擦副在潤滑條件下的摩擦磨損性能進行了大量的研究,取得了一些成果,如國外學者BRODSKII和CHARRIER等[1-2]發現了高分子材料磨損橡膠的特殊現象。近年來,國內外學者應用試驗和數值模擬的手段分析了表面形貌對彈流潤滑區域、混合潤滑區域和接觸區域的潤滑特性及摩擦磨損特性的影響。何仁洋、柳瓊俊等[3-4]通過實驗方法探討了干摩擦和邊界潤滑條件下橡膠對金屬的摩擦磨損機制;REN、KRUPKA等[5-6]討論了表面粗糙對潤滑膜形成的影響,認為通過合理設計粗糙度的表面紋理可提高表面潤滑性能;WANG等[7]通過反復試驗點接觸下表面粗糙度對滑動摩擦的影響,發現橫向紋理比縱向紋理的摩擦因數小;LITWIN等[8-10]采用試驗方法研究了不同表面形貌和不同材料下的水潤滑軸承流體動力學潤滑特性;CHOO等[11]開發了一種新的試驗方法來研究表面粗糙度對潤滑膜分布的影響;ZHANG等[12]分析了表面粗糙度對彈塑性接觸力學性能的影響;范曉夢等[13]分析了瞬變載荷作用下水潤滑軸承彈流潤滑性能分析;張闊[14]利用數值模擬方法研究了材料的表面形貌對混合潤滑下摩擦性能的影響。最近,楊勇等人[15]基于快速傅立葉變換的數值仿真方法生成了給定參數分布的非高斯曲面,求得了不同表面粗糙度特征下的齒輪接觸區域潤滑特性;劉明勇等[16]探討了二維規則下表面粗糙度的波長、幅值及工況變化下的潤滑特性;GHERCA、WANG等[17-18]通過數值模擬和實驗相結合的方法,通過優化推力軸承的表面紋理提高了軸承的流體動力學及摩擦潤滑特性;ZHU和HUANG[19]基于軸承曲率的形貌參數對摩擦磨損性能做了探討。
上述研究主要集中在金屬軸頸表面粗糙度對軸承潤滑膜的形成及潤滑性能的影響,但研究金屬軸頸表面粗糙度對橡膠軸承-軸頸摩擦副干摩擦特性影響的文章較少。為研究軸頸表面粗糙度對水潤滑橡膠軸承在干摩擦狀態下摩擦特性的影響,本文作者采用試驗與理論計算相結合的方法,探討金屬軸頸-橡膠材料摩擦副干摩擦狀態下摩擦因數的計算方法,為水潤滑橡膠軸承摩擦性能的研究提供依據。
船舶艉軸與橡膠軸承在相對運動時,當軸承供水系統出現故障或軸承過載時,軸頸與橡膠軸承之間會發生局部接觸。從微觀來看,任何真實的工程表面都不是絕對光滑的,而是存在一定的粗糙度[20],這個粗糙度使得兩物體之間的接觸變成微凸體之間的接觸。
圖1所示為粗糙剛性表面與橡膠材料的接觸模型,橡膠在一個粗糙表面上滑動,在此情況下,摩擦力主要取決于金屬表面的粗糙度,橡膠表面粗糙度所起作用很小,因此問題可以簡化為光滑橡膠與粗糙金屬之間的摩擦。由于在局部接觸時橡膠軸承和軸頸之間存在納米尺度的薄膜[21-22],可假設橡膠材料和軸頸之間無黏著力,此時滯后摩擦起主要作用,即滑動橡膠材料在周期變形過程中產生的能量損耗,此時問題就轉化為計算橡膠材料的變形和能量損耗[20]。

圖1 粗糙剛性表面與橡膠材料的接觸模型
表面形貌通常以表面粗糙度的分布形式表征,因此金屬的表面粗糙度分布是研究的重點。

圖2 微凸體的表面分布模型

(1)
一個微凸體的特征直徑[20]:
(2)
微凸體接觸的單位體積能量損耗[20]:
(3)
摩擦因數[20]:
μ=ζ
(4)
μ≈z
(5)
以某水潤滑橡膠軸承主軸為測量對象,利用TIME3230表面粗糙度測量儀對金屬軸徑表面的微凸體位置分布參數進行測量。軸頸的材料為40Cr,表面粗糙度為Ra0.8 μm,調質處理,表面磨削。
采用圖3所示TIME3230表面粗糙度測量儀對沿加工紋理方向、垂直于加工紋理方向及與加工紋理方向呈45°夾角方向的金屬軸頸表面的微凸體高度分布進行測量。

圖3 TIME3230表面粗糙度測量儀
微凸體位置分布參數的測量原理如圖4所示,粗糙度形狀測量儀采用大量程的差動電感傳感器拾取被測信號,測量軸表面時,將傳感器放在軸表面上,由驅動器帶動傳感器沿軸表面做等速滑行,驅動器中的導軌保證軌跡的精度。傳感器通過銳利的金剛石觸針感受被測表面的輪廓信號,此時軸表面輪廓起伏引起觸針產生垂直位移,并輸出運行軌跡垂直位移的坐標。
測量方法:軸頸低速運動,在恒定載荷下,采用表面粗糙度形狀測量儀沿軸頸的加工紋理方向、垂直加工紋理方向及與加工紋理呈45°夾角方向共3個方向進行數據采集。每個方向采集4組數據,每隔0.001 μm設置一個數據采集點采集數據一次,取樣長度為0.8 mm,評定長度取5倍的取樣長度。最終得到12組金屬軸徑表面微凸體位置分布數據。
取沿軸頸加工紋理方向的原始測量數據2組,如圖5所示。它包含一定的隨機噪聲,同時,機械加工零件表面存在較大的各向異性,這些都使得測量的表面數據應用較為困難。因此,采用Savitzky-Golay算法進行平滑去噪,有效去除測量信號中的脈沖噪聲和白噪聲,去噪后的結果如圖6所示,它是一系列隨機的表面波,且可以看出去噪后的數據保有效保留了原始數據的特征。

圖5 沿加工紋理方向的粗糙度原始數據分布曲線

圖6 去噪后的沿加工紋理方向的粗糙度數據分布曲線
任何連續的時序或信號,都可以表示為不同頻率的正弦波信號的無限疊加,因此傅立葉變換可以將表面粗糙度儀測得的隨機表面波分解為有一定規律的正弦波分量的組合。利用Matlab軟件中的傅立葉擬合函數將去噪后的原始數據進行擬合變換,經過傅立葉擬合變換后重新構造的軸頸表面粗糙度分布表達式f(x)為
(6)
式中:A為波幅;k為特征波波長,令λ=1/k為波長系數。
采用7階傅立葉變換,將去噪后的數據進行變換處理,得到粗糙度分布結果如圖7所示。

圖7 采用傅立葉變換后的沿加工紋理方向的粗糙度數據分布曲線
觀察采用傅立葉變換后的曲線分布可發現,新擬合的曲線與去噪后的原始數據非常接近,可以真實地反映粗糙軸表面的粗糙度輪廓特點。因此軸頸的表面粗糙度分布可以看作是7個正弦函數和7個余弦函數的疊加,因此它們對摩擦因數的貢獻相加等于[20]:

(7)
不同潤滑狀態下的摩擦因數如表1[24]所示。將采用傅立葉變換后的各項系數代入摩擦因數計算表達式(7),得μ=0.167 2。由表1可以判斷此時的潤滑狀態處于干摩擦狀態,這與實際情況吻合。

表1 不同潤滑狀態下的摩擦因數
3.4 粗糙度幅值和特征波波長系數對摩擦因數的影響
在軸頸轉速一定的情況下,利用式(7)計算了在不同粗糙度幅值0.5A、A、1.5A、2A、2.5A和不同特征波波長系數0.5λ、λ、1.5λ、2.0λ、2.5λ下的摩擦因數,結果如圖8所示。

圖8 不同幅值和特征波波長系數下的摩擦因數
由圖8可以看出:隨粗糙度函數幅值從0.5A逐漸增加到2.5A,摩擦因數從0.130 9增加到0.543 4,這是由于隨著微凸體粗糙度函數幅值的增大,加工表面越粗糙,微凸體特征高度的正負高差越大,微凸體表面均方根l變大,使得表面輪廓梯度的平方根z也增大,摩擦因數增大;隨粗糙度波長系數從0.5λ逐漸增加到2.5λ,摩擦因數從0.157 2增加到0.496 8,這是由于隨著特征波波長系數的增大,表面微凸體的密度變大,表面間的接觸面積增大,表面輪廓梯度的平方根z也增大,使得摩擦因數增大。
在實際工程中,表面紋理方向多與潤滑流動方向呈一定夾角。采用表面儀測量了與加工紋理方向呈45°夾角和垂直加工紋理方向的粗糙度,并進行傅立葉變換構造處理,利用公式(7)計算出這2個不同方向下的摩擦因數,并與順加工紋理方向的摩擦因數進行比較,如表2所示。

表2 不同紋理下的摩擦因數
由表2可以看出:潤滑流動方向順著加工紋理方向時摩擦因數最小,垂直加工紋理方向時摩擦因數最大。這是由于潤滑流動方向順加工紋理方向時,微凸體特征高度正負高差較小,同時微凸體密度也較小,波動的幅值和特征波長均較小,使得z變小,摩擦因數變小;而垂直加工紋理時,微凸體特征高度正負高差較大,同時微凸體密度也較大,波動的幅值和特征波長均較大,使得z變大,摩擦因數變大;當與加工紋理方向呈45°時,微凸體特征高度和密度均居于二者之間,因此其摩擦因數大于順加工紋理方向的摩擦因數,小于垂直于加工紋理方向的摩擦因數。
因此,選擇合理的軸承表面粗糙度函數的幅值和波長可以提高金屬軸頸-橡膠軸承摩擦副的摩擦潤滑性能;同時,沿加工紋理加工、裝配金屬軸頸-橡膠軸承摩擦副,可降低摩擦因數,從而有效避免軸承的打滑擦傷失效。
(1)在干摩擦狀態下,水潤滑橡膠軸承-軸頸副的摩擦因數隨表面粗糙度函數幅值的增大而增大,且隨著粗糙度函數特征波長系數的增大非線性增加。潤滑流動方向順著加工紋理方向時摩擦因數最小,潤滑流動方向與加工紋理方向呈45°時摩擦因數次之,潤滑流動方向垂直加工紋理方向時摩擦因數最大。
(2)利用傅立葉變換重新構造表面粗糙度分布模型,并依據理論計算橡膠材料-軸頸副的摩擦因數,計算結果與實際情況吻合,驗證了模型的正確性。
(3)研究表明:選擇合理的軸承表面粗糙度函數的幅值和波長可以提高金屬軸頸-橡膠軸承摩擦副的摩擦潤滑性能;同時,沿加工紋理加工、裝配金屬軸頸-橡膠軸承摩擦副,可降低摩擦因數,從而有效避免軸承的打滑擦傷失效。上述試驗研究可為后續研究分析水潤滑橡膠軸承的動力學特性提供理論依據。