馬 歡 景 卉
(1.西安交通大學(xué)航天航空學(xué)院 陜西西安 710048;2.蘭州大學(xué)管理學(xué)院 甘肅蘭州 730000)
機(jī)械結(jié)合面是機(jī)械裝備中各零部件之間的接觸表面,對于各種高精密設(shè)備,結(jié)合面的接觸剛度對設(shè)備的精度和性能有著重要影響。通常情況下,為了減少磨損,提高設(shè)備使用壽命,結(jié)合面往往都處于潤滑狀態(tài)[1]。機(jī)械結(jié)合面一般可以等效為2個粗糙表面之間的接觸,因此,潤滑狀態(tài)下的粗糙表面間接觸剛度的研究意義重大。
自經(jīng)典赫茲接觸[2]被提出以來,國內(nèi)外學(xué)者對粗糙表面之間的接觸進(jìn)行了廣泛研究[3-10]。其中以GREENWOOD-WILLIAMSON(G-W)模型[3]為代表的統(tǒng)計(jì)學(xué)接觸模型和Persson接觸理論[5]為代表的分形接觸模型被廣泛認(rèn)可。在G-W模型中,兩個粗糙面之間的接觸往往等效為一個剛性平面和一個彈性粗糙表面之間的接觸,粗糙面由曲率相同的球面組成,且其高度服從高斯分布。但由于粗糙表面的多尺度特征,統(tǒng)計(jì)學(xué)模型會受儀器分辨率的影響。為了更好地表征粗糙表面的多尺度特征,PERSSON[5]通過引入放大率概念,發(fā)展了一種基于壓力擴(kuò)散方法的隨機(jī)粗糙表面接觸理論。然而,上述模型僅僅考慮了干接觸條件下的粗糙表面接觸,沒有考慮潤滑對接觸特性的影響。
在實(shí)際工業(yè)應(yīng)用中,各種齒輪、軸承等機(jī)械系統(tǒng)以混合潤滑接觸方式工作,以減少摩擦或磨損。超聲反射法是實(shí)驗(yàn)測量界面接觸剛度的主要手段。GONZALEZ-VALADEZ等[11]利用超聲技術(shù)測量了潤滑界面接觸剛度。然而,基于超聲反射法測得的潤滑粗糙界面的接觸剛度遠(yuǎn)大于實(shí)際值[12]。因此,研究一種理論模型預(yù)測混合潤滑狀態(tài)下的接觸剛度是非常必要的。
本文作者提出了一個混合潤滑狀態(tài)下粗糙表面的接觸剛度等效薄層模型,研究了不同實(shí)際接觸面積下,表面形貌和潤滑劑對法向接觸剛度的影響,為研究混合潤滑狀態(tài)下粗糙界面的接觸特性提供了理論基礎(chǔ)。
潤滑狀態(tài)下的粗糙接觸界面由粗糙峰和周圍的潤滑劑組成,如圖1所示。對于混合潤滑粗糙界面,界面剛度由固體剛度和潤滑劑剛度2部分組成,且其剛度關(guān)系符合彈簧并聯(lián)準(zhǔn)則,即
圖1 混合潤滑狀態(tài)下的粗糙界面
Kn=Ks+Kl
(1)
式中:Kn為潤滑狀態(tài)下粗糙界面總接觸剛度;Ks為固體剛度;Kl為潤滑劑剛度。
對于彈性潤滑界面,假設(shè)等效薄層與接觸固體和平坦剛性表面剛性連接,如圖2所示。因此,單位標(biāo)稱接觸面積的法向接觸剛度[12]可以表示為
圖2 混合潤滑接觸等效薄層示意
(2)
式中:H0為等效薄層厚度,基于粗糙統(tǒng)計(jì)理論,H0一般取值為3倍的均方根粗糙度[13];E0為等效彈性模量,可表示為
E0=φ1Es+φ2Bl
(3)
式中:Es為固體的彈性模量;Bl為潤滑劑的體積模量;φ1、φ2分別為等效彈性模量中固體和潤滑劑的貢獻(xiàn)比例,并滿足
φ1+φ2=1
(4)
因此,固體剛度和潤滑劑剛度可分別表示為
(5)
(6)
文中推導(dǎo)潤滑接觸狀態(tài)下界面接觸剛度的表達(dá)式。首先,基于WANG等[8]提出的干接觸時粗糙表面增量接觸模型推導(dǎo)得到單位標(biāo)稱面積固體接觸剛度表達(dá)式為
(7)
式中:A0為標(biāo)稱接觸面積;A(z)和N(z)分別為當(dāng)粗糙界面分離為z時,粗糙表面的實(shí)際接觸面積和等效圓斑數(shù)量。
將式(7)代入式(5)可得固體部分的貢獻(xiàn)比例φ1為
(8)
將式(8)代入式(4)可得潤滑劑部分的貢獻(xiàn)比例φ2為
(9)
將式(8)和式(9)代入式(3)可得潤滑接觸等效彈性模量E0為
(10)
將式(10)代入式(2)可得潤滑接觸界面法向剛度Kn為
(11)
將式(9)代入式(6)可得潤滑劑剛度Kl為
(12)
盡管固體的彈性模量遠(yuǎn)遠(yuǎn)大于潤滑劑的體積模量,潤滑劑對剛度的貢獻(xiàn)也不可忽略,因?yàn)樵诖植诜宄跏冀佑|時,固體的貢獻(xiàn)比例接近于0,而潤滑劑貢獻(xiàn)比例接近于1。
根據(jù)1.2節(jié)中潤滑接觸狀態(tài)下界面接觸剛度模型,只需粗糙表面的高度節(jié)點(diǎn)信息,便可預(yù)測混合潤滑狀態(tài)下接觸剛度的演化。首先,通過PSD法在MATLAB中生成一個100 μm×100 μm的隨機(jī)粗糙表面,其節(jié)點(diǎn)高度信息如圖3所示。粗糙表面的具體參數(shù)為:σ=3 μm,Es=115.4 GPa,Bl=1.45 GPa,H=0.8,ν=0.3,其中H為赫斯特指數(shù),ν為泊松比?;谝陨闲畔ⅲ玫刃П幽P蛿?shù)值求解混合潤滑狀態(tài)下接觸剛度。如圖4所示,隨著實(shí)際接觸面積的增大,接觸剛度近似線性增加,與G-W模型[3]和Persson接觸理論[5]所得到的線性關(guān)系一致。對于相同的實(shí)際接觸面積,有潤滑劑界面的接觸剛度大于無潤滑劑界面的接觸剛度。固體和潤滑劑分別對總接觸剛度的貢獻(xiàn)比例隨著實(shí)際接觸面積增大的變化規(guī)律如圖5所示。初始接觸時φ2>φ1,潤滑劑作用占主導(dǎo)地位,但隨著實(shí)際接觸面積的增大,固體接觸剛度逐漸超過潤滑劑接觸剛度占據(jù)主導(dǎo)地位。這說明在小載荷下,潤滑劑的接觸剛度對于總接觸剛度具有重要意義。
圖3 粗糙表面形貌
圖4 有無潤滑劑條件下接觸剛度與面積比的關(guān)系
圖5 固體和潤滑劑的貢獻(xiàn)比例與面積比的關(guān)系
為了研究潤滑劑類型對混合潤滑接觸剛度的影響,文中比較了2種不同類型的潤滑劑對接觸剛度的影響。潤滑劑在環(huán)境壓力下的體積模量為
Bl=Bl0e-βkT
(13)
式中:Bl0為環(huán)境壓力下溫度為0 K時的體積模量;βk為泰特-杜利特爾常數(shù);T為絕對溫度。
表1列出了分析中使用的2種潤滑劑的相關(guān)參數(shù)[14]。
表1 環(huán)境壓力下潤滑劑的參數(shù)[14]
圖6(a)和圖6(b)示出了在分別使用T9油和甘油2種潤滑劑下,總接觸剛度、固體接觸剛度和潤滑劑接觸剛度隨實(shí)際接觸面積的變化規(guī)律。可知,總接觸剛度、固體接觸剛度隨著實(shí)際接觸面積的增大而增大,且總接觸剛度總是大于固體接觸剛度,而潤滑劑接觸剛度基本保持不變,并且在小載荷下潤滑劑接觸剛度大于固體接觸剛度。通過對比2種不同的潤滑劑,可以發(fā)現(xiàn)潤滑劑的體積模量越大,混合潤滑接觸剛度增大越明顯,潤滑劑接觸剛度占主導(dǎo)地位的范圍也越大,這與實(shí)際情況是一致的。
圖6 不同潤滑劑下法向接觸剛度與面積比的關(guān)系
為了探究表面形貌對混合潤滑接觸剛度的影響,通過PSD法生成了3種不同粗糙度的隨機(jī)粗糙表面,具體參數(shù)見表2。表面形貌對T9油潤滑界面接觸剛度的影響如圖7所示,不同粗糙度的潤滑界面的法向剛度均隨著實(shí)際接觸面積的增大而增大。在實(shí)際接觸面積相同的情況下,更光滑的表面會產(chǎn)生更大的接觸剛度。
表2 不同表面形貌的具體參數(shù)
圖7 表面形貌對混合潤滑接觸剛度的影響
為了驗(yàn)證文中提出的接觸剛度模型的準(zhǔn)確性,將文中模型與文獻(xiàn)[13]的模型進(jìn)行了比較。選用的參數(shù)為σ=2.04 μm,Es=115.4 GPa,Bl=1.45 GPa,H=0.7,ν=0.3。如圖8所示,當(dāng)實(shí)際接觸面積相同時,文中模型的預(yù)測結(jié)果與文獻(xiàn)[13]的模型的計(jì)算結(jié)果具有良好的一致性,證明了文中提出的混合潤滑狀態(tài)下粗糙界面的法向接觸剛度模型的合理性與有效性。與現(xiàn)有模型相比,文中模型的創(chuàng)新之處在于只需要提供粗糙表面的高度信息,便可以簡單準(zhǔn)確地預(yù)測干接觸和潤滑接觸時接觸剛度隨實(shí)際接觸面積的演化。
圖8 不同模型法向接觸剛度與面積比的關(guān)系
通過提出一種等效薄層模型解決混合潤滑狀態(tài)下粗糙界面接觸剛度的預(yù)測問題,將接觸界面的總剛度等效為固體接觸剛度和潤滑劑接觸剛度之和,從理論上估算了混合潤滑法向接觸剛度與實(shí)際接觸面積的關(guān)系。主要研究結(jié)果如下:
(1)粗糙界面的法向接觸剛度隨法向載荷的增加而增加,且混合潤滑狀態(tài)下的接觸剛度大于干接觸條件下的接觸剛度。
(2)在初始接觸時,粗糙界面的法向剛度主要受潤滑劑性能影響,隨著實(shí)際接觸面積的增大,固體剛度逐漸占據(jù)主導(dǎo)地位,表面形貌對法向接觸剛度的影響變得更加明顯,潤滑劑的影響越來越小。
(3)建立的模型的預(yù)測結(jié)果與經(jīng)典模型計(jì)算結(jié)果具有良好的一致性,驗(yàn)證了建立的模型的正確性。