張 穎，張學(xué)良，姜 來，溫淑花，陳永會，蘭國生，劉麗琴(太原科技大學(xué)機械工程學(xué)院，太原 030024)

據(jù)美國材料試驗學(xué)會ASTM的分析報告指出，全世界約有40%以上的各類能源被摩擦消耗掉，摩擦?xí)鹨幌盗械奈锢砗突瘜W(xué)變化，產(chǎn)生熱、噪音、振動、污染，甚至?xí)l(fā)危險。研究和預(yù)測機械系統(tǒng)的摩擦同樣是非常重要的。而靜摩擦因數(shù)是研究中摩擦的一個十分重要的參數(shù)。因此，關(guān)于結(jié)合面靜摩擦因數(shù)的研究具有十分重要的意義。Chang等[1]人基于經(jīng)典的G-W接觸模型[2]，通過考慮彈性接觸點承受的切向載荷的能力，提出了金屬表面間的靜摩擦因數(shù)計算模型。Kogut和Etsion通過研究發(fā)現(xiàn)，結(jié)合面彈-塑性接觸點同樣具有承受切向載荷的能力，如果僅僅考慮彈性接觸點承受切向載荷的靜摩擦模型，顯然對靜摩擦因數(shù)的預(yù)測是不精確的。同時以上的研究都是建立在二維分形模型基礎(chǔ)上的。應(yīng)用二維的分形模型計算三維粗糙表面的接觸特性，存在著理論上的不合理性[3]。鑒于此，文獻[4]提出了一種基于三維接觸分形理論的結(jié)合面靜摩擦因數(shù)模型，但該模型沒有考慮彈-塑性變形階段對靜摩擦因數(shù)的影響。

本文基于三維接觸分形理論及文獻[5]的彈塑性變形機制，綜合考慮結(jié)合面微凸體的三種變形機制，從而建立了一種考慮彈-塑性變形的結(jié)合面靜摩擦因數(shù)三維分形模型，并進行了仿真計算分析。

1 結(jié)合面三維分形接觸模型

1.1 結(jié)合面三維建模

結(jié)合面實質(zhì)就是兩個粗糙表面的接觸，在研究的過程中，為了簡化通常將結(jié)合面的接觸表示為一個剛性平面與一個粗糙平面的接觸。當結(jié)合面受到擠壓，微凸體將發(fā)生變形時，單個微凸體和剛性表面的彈性接觸如圖1所示。

結(jié)合面的復(fù)合彈性模量為:

式中，E*為兩接觸材料的復(fù)合彈性模量；E1、E2為兩接觸材料的彈性模量；v1、v2為兩接觸材料的泊松比。

在該模型中假設(shè)了微凸體的分布是各向同性的，

圖1 單個微凸體與剛性平面彈性接觸示意圖

微凸體之間的相互作用忽略不計，并且沒有宏觀的變形，變形過程只發(fā)生在微凸體上。

根據(jù)文獻[3]的研究結(jié)果可知，當法向變形量δ處于不同范圍時，微凸體將發(fā)生完全彈性變形、彈-塑性變形及完全塑性變形。

單個微凸體的圓形接觸面積a、截斷面積a′和法向接觸載荷Pe，當微凸體的變形在完全彈性變形階段時可以表示為:

(1)

ae=πRδ

(2)

其中下標e表示發(fā)生彈性變形。

當微凸體的變形在彈-塑性變形階段時，彈-塑性變形區(qū)域的范圍為δc≤δ≤110δc，并將彈-塑性區(qū)域進一步劃分為兩個區(qū)域：彈-塑性變形第一區(qū)域(δc≤δ≤6δc)和彈-塑性變形第二區(qū)(6δc≤δ≤110δc).微凸體彈-塑性變形階段法向接觸載荷、接觸面積和變形量之間的關(guān)系為[5]:

(3)

(4)

(5)

(6)

式中下標ep1、ep2分別表示發(fā)生彈-塑性第一區(qū)變形、彈-塑性第二區(qū)變形。

彈性臨界變形量δc為[5]:

(7)

當微凸體的變形在完全塑性變形時，微凸體的變形量在δp>110δc范圍內(nèi)，微凸體在塑性變形階段的接觸面積和法向接觸載荷Pp可以表示為:

Pp=Hap

(8)

ap=a′

(9)

(10)

粗糙表面上的橫截微凸體大小分布函數(shù)為[3]:

(11)

結(jié)合面真實接觸面積Ar通過積分得到:

(12)

1.2 結(jié)合面法向載荷計算模型

微凸體的變形量δ是由W-M函數(shù)所決定的，其值取決于余弦函數(shù)的波峰和波谷的幅值，即[3]:

(13)

微凸體的變形量δ、微接觸點的截半徑r′、微凸體的曲率半徑R的幾何關(guān)系為:

(R-δ)2+(r′)2=R2

(14)

又因為δ<<2R，所以:

(15)

將式(13)代入式(15)得[3]：

(16)

結(jié)合面法向載荷與真實接觸面積之間的關(guān)系式為:

(17)

當2

(18)

當D=2.5時，

(19)

根據(jù)接觸點的實際接觸面積a與其截面積a′之間的關(guān)系a′=2a，有:

(20)

(21)

將式(20)、式(21)帶入式(18)、式(19)得：當2

(22)

當D=2.5時，

(23)

將式(22)、式(23)無量綱化，得：當2

(24)

當D=2.5時，

(25)

2 考慮彈塑性變形的靜摩擦因數(shù)三維分形模型

根據(jù)文獻[6]可知，當微凸體的接觸界面達到完全屈服時，結(jié)合面單個微凸體承受的切向載荷:

完全彈性區(qū)和彈-塑性一區(qū)的微凸體的接觸面積和接觸載荷分別為式(1)～式(4).因此完全彈性區(qū)和彈-塑性一區(qū)的微凸體的切向載荷分別為：

完全彈性區(qū)和彈-塑性一區(qū)的微凸體，在結(jié)合面處于臨界滑動狀態(tài)時，接觸界面都將達到完全屈服，所以此時微凸體所受到的切向的總載荷Q就是最大靜摩擦力[6]。

(26)

(27)

(28)

將式(20)、式(21)帶入式(26)得：

(29)

將式(27)、式(28)帶入式(29)得：當2

(30)

當D=2.5時，

(31)

將式(30)、式(31)無量綱化，得：

當2

(32)

當D=2.5時，

(33)

圖2 D對f的影響

靜摩擦因數(shù)的計算公式為：

(34)

將式(24)、式(25)、式(32)、式(33)帶入式(34)得結(jié)合面靜摩擦因數(shù)。

3 模型仿真

圖3 P*對f的影響(K=2.8,φ=1.0,G*=1.0×10-5)

圖4 G*對f的影響(K=2.8,φ=1.0)

圖5 忽略與包含彈塑性過渡變形機制的兩種模型仿真結(jié)果對比

由仿真結(jié)果可見

1)由圖2可知，結(jié)合面靜摩擦因數(shù)隨著分形維數(shù)的增大，先增大再減小。

2)由圖3可知，結(jié)合面靜摩擦因數(shù)隨著無量綱法向接觸總載荷P*的增大而增大，并且呈凸弧非線性關(guān)系。

3)由圖4可知，結(jié)合面靜摩擦因數(shù)隨著無量綱分形特征長度尺度參數(shù)G*的增大而減小。

文獻[4]中的靜摩擦因數(shù)三維分形模型與本文考慮彈-塑性變形機制的摩擦因數(shù)三維分形模型仿真結(jié)果對比如圖5所示。由圖5可以看出，忽略彈-塑性與包含彈-塑性模型仿真結(jié)果表明，靜摩擦因數(shù)都是隨著分形維數(shù)的增大，先增大然后減??；隨著無量綱法向總載荷的增大而增大。兩種模型的曲線走勢接近，但由于考慮了彈-塑性變形機制，摩擦因數(shù)的值均小于忽略彈-塑性的摩擦因數(shù)的值。

4 結(jié)論

(1)本文提出了一種基于三維分形理論，考慮彈-塑性過渡變形機制的結(jié)合面靜摩擦因數(shù)三維分形模型，該模型可以用于進行結(jié)合面靜摩擦因數(shù)的預(yù)測。

(2)結(jié)合面靜摩擦因數(shù)隨著分形維數(shù)的增大，先增大后減??；隨著分形特征長度尺度參數(shù)的增大而減?。浑S著法向總接觸載荷的增大而增大。

(3)考慮彈-塑性過渡變形機制的結(jié)合面靜摩擦因數(shù)要小于不考慮彈-塑性變形機制的結(jié)合面靜摩擦因數(shù)。

參考文獻：

[1] CHANG W R，ETSION I，BOGY D B.Static Friction Coefficient model for metallic rough surfaces[J].J Tribol ASME，1988，110(1)：57-63.

[2] GREENWOOD J A，WILLIAMSON J B P.Contact of nominally flat surfaces[J].Proceedings of The Royal Society of London，1966，A295：300-317.

[3] YAN W，KOMVOPOULOS K.Contact analysis of elastic-plastic fractal surfaces[J].J Appl Phys，1998，84：3617-3624.

[4] 殷東華，張學(xué)良，溫淑花.結(jié)合面靜摩擦因數(shù)三維分形模型[J].太原科技大學(xué)學(xué)報，2013，34(6)：451-455.

[5] LIOU J L.The Theoretical study for microcontact model with variable topography parameters[D].Tiwan：National Cheng Kung University，2006.

[6] 尤晉閩，陳天寧.結(jié)合面靜摩擦系數(shù)的統(tǒng)計模型[J].振動與沖擊，2010，29(12)：26-29.

[7] 張學(xué)良.機械結(jié)合面動態(tài)特性及應(yīng)用[M].北京：中國科技出版社，2002.

[8] KOGUT L，ETSION I.A Static friction model for elastic-plastic contacting rough surfaces[J].ASME Journal of Tribology，2004，126(1)：34-40.

[9] KOGUT L，ETSION I.A semi-analytical solution for the sliding inception of a spherical contact [J].ASME Journal of Tribology，2003，125(3)：499-506.

[10] 蘭國生，張學(xué)良，丁紅欽，等.基于分形理論的結(jié)合面靜摩擦因數(shù)改進模型[J].農(nóng)業(yè)機械學(xué)報，2012，43(1)：213-218.