基于微凸體連續(xù)變形理論的結(jié)合面切向剛度模型

王世軍，李志濤，韓子銳，趙金娟，李鵬陽

(1.西安理工大學(xué) 機(jī)械與精密儀器工程學(xué)院，陜西西安710048；2.西安理工大學(xué) 印刷包裝與數(shù)字媒體學(xué)院，陜西西安710048)

結(jié)合面的接觸剛度包含法向接觸剛度和切向接觸剛度，接觸研究中通常單獨(dú)分析法向剛度和切向剛度，應(yīng)用中則需要同時(shí)考慮法向剛度和切向剛度的影響。

傳統(tǒng)的法向接觸模型，如基于統(tǒng)計(jì)分析的G-W模型[1]和基于分形理論的M-B模型[2]，以及張學(xué)良[3]和尤晉閩[4]的模型，都沒有考慮彈塑性過渡階段的接觸變形，只把彈性與塑性接觸變形計(jì)入模型。Zhao[5]首次給出了一個(gè)法向接觸模型(ZMC模型)，該模型含有彈性、彈塑性以及塑性。隨后繆小梅[6]得到考慮微凸體彈塑性過渡階段的法向接觸剛度模型，王世軍[7]則基于微凸體彈塑性有限元接觸模型，建立了微凸體在彈性、彈塑性、塑性變形階段統(tǒng)一的接觸變量變化規(guī)律的表達(dá)式。何聯(lián)格[8]進(jìn)一步將彈塑性過渡階段分為三個(gè)不同的小階段，得到了更加準(zhǔn)確的法向接觸剛度模型。

在切向剛度的研究中，張學(xué)良[9]首次建立了具有尺度獨(dú)立性的結(jié)合部切向剛度分形模型，但是在該模型中僅僅將接觸變形過程劃分為完全彈性和完全塑性階段。CHEN[10]基于分形理論建立了曲面切向接觸剛度的分形模型，但是對(duì)微凸體的變形過程的劃分與張學(xué)良[9]相同,并未考慮微凸體變形過程中的彈塑性過渡階段。李小彭[11]在分形理論和摩擦學(xué)原理的基礎(chǔ)上，給出了考慮摩擦因素的切向剛度模型，此切向剛度模型中也僅僅考慮了微凸體在完全彈性過程中產(chǎn)生的切向接觸剛度。田紅亮等[12]改進(jìn)分形幾何理論并且嚴(yán)格應(yīng)用Hertz法向接觸力學(xué)得到了一個(gè)改進(jìn)的切向接觸剛度模型，微凸體的剛度分析也只考慮完全彈性階段。劉鵬等[13]對(duì)考慮摩擦因素以及球面接觸時(shí)切向接觸剛度進(jìn)行了研究，在得到的模型中仍然將接觸變形劃分成不連續(xù)的完全彈性變形過程以及完全塑性變形過程，但仍然未考慮接觸變形過程中的彈塑性過渡階段。上述切向接觸模型的建立都將微凸體的變形過程劃分為完全彈性和完全塑性兩個(gè)不連續(xù)的變形階段，僅僅考慮了微凸體在彈性階段的所產(chǎn)生的切向剛度。

近年來微凸體法向接觸特性的研究表明[14-16]，在逐漸增大法向載荷過程中，微凸體中的塑性區(qū)從無到有，由小變大，彈性區(qū)則由大變小，接觸面上微凸體的材料性質(zhì)由彈性逐漸變?yōu)樗苄裕⑼贵w的變形過程從完全彈性階段到完全塑性階段需要經(jīng)過一個(gè)彈塑性過渡階段。在這個(gè)過渡階段，微凸體在接觸面的切線方向上仍然具有彈性。微凸體切向接觸剛度分析中應(yīng)當(dāng)與法向接觸剛度分析一樣[5-8]計(jì)入彈塑性過渡階段的切向接觸剛度，不應(yīng)該只考慮彈性階段的切向接觸剛度[9-13]。

本文依據(jù)Hertz理論，分形幾何理論和切向接觸理論，首次考慮彈塑性過渡階段中屈服區(qū)域沒有擴(kuò)大至整個(gè)接觸表面前的彈塑性過渡階段所產(chǎn)生的切向剛度，推導(dǎo)出一個(gè)新的接觸表面的切向剛度模型。

1 結(jié)合面接觸變形基本理論

1.1 結(jié)合面微凸體接觸變形

結(jié)合面在相互接觸時(shí)，是非常復(fù)雜的，因此為了方便建模，可以將接觸進(jìn)行等效，即互相接觸表面為粗糙表面和光滑的剛性平面。在整個(gè)粗糙表面接觸時(shí)，將整個(gè)表面上的單個(gè)微凸體等效為球體，當(dāng)承受法向載荷P的作用時(shí)，其接觸狀態(tài)見圖1。其中R為微凸體等效的曲率半徑，ω為微凸體的法向變形量，P為施加在微凸體上的法向載荷，r為接觸區(qū)域的半徑。

圖1 單個(gè)微凸體與剛性平面的接觸狀態(tài)Fig.1 Contact state of a single asperity with a rigid plane

1.1.1完全彈性接觸狀態(tài)

在法向變形量ω<ωec的情況下，微凸體此時(shí)的接觸狀態(tài)為完全彈性變形。法向接觸載荷Fe和微凸體的臨界變形量ωec為[15]：

(1)

(2)

(3)

式中：D表示分形維數(shù)；G表示分形粗糙度；H表示兩材料中較軟材料的硬度，H=2.8σy，這里σy表示較軟材料的屈服強(qiáng)度。K表示硬度系數(shù)，K=0.454+0.41v，這里v表示兩材料中較軟材料的Possion比。E表示兩接觸材料的當(dāng)量彈性模量，E1和E2分別表示材料各自的彈性模量，v1和v2分別表示材料各自的Possion比。

1.1.2完全塑性接觸狀態(tài)

根據(jù)文獻(xiàn)[3]的研究，當(dāng)法向變形量ω≥ωpc=110ωec時(shí)(ωpc為發(fā)生塑性變形的臨界法向變形量)，微凸體進(jìn)入完全塑性階段，此時(shí)平均接觸壓力等于較軟材料的硬度，單個(gè)微凸體法向接觸載荷FP為[15]：

FP=2πHRω

(4)

1.1.3彈塑性過渡接觸狀態(tài)

當(dāng)法向變形量ωec≤ω<ωpc時(shí)，Kogut[16]認(rèn)為微凸體處于彈塑性過度階段，發(fā)生彈塑性變形，并且可將該階段分為三個(gè)小階段，即ωec≤ω<6ωec、6ωec≤ω<68ωec和68ωec≤ω<110ωec。當(dāng)法向變形量處于ωec≤ω<6ωec時(shí)，微凸體的屈服區(qū)域均發(fā)生在接觸面下方，塑性變形發(fā)生在接觸表面之下的微小體積內(nèi)，塑性區(qū)域被大量的彈性區(qū)包圍，且實(shí)際接觸面積上每一點(diǎn)均發(fā)生彈性變形。由于接觸面積上都為彈性變形，因此在該階段，實(shí)際接觸表面還具有抵抗切向變形的能力，即具有切向剛度。隨著法向變形量的不斷增大，當(dāng)ω=6ωec時(shí)，屈服區(qū)域由接觸表面下方擴(kuò)大至接觸表面。當(dāng)法向變形量處于6ωec≤ω<68ωec時(shí)，法向接觸載荷不斷增大，屈服區(qū)域在接觸表面隨之增大。由于在該階段實(shí)際接觸表面還是存在彈性接觸，因此也具有一定的抵抗切向變形的能力，只不過是隨著法向接觸載荷的增大以及屈服區(qū)域的增大，抵抗切向變形的能力不斷減小，即切向剛度在該階段不斷減小。當(dāng)法向變形量增大至ω=68ωec時(shí)，屈服區(qū)域擴(kuò)大至整個(gè)接觸表面，此時(shí)以及之后的階段，由于整個(gè)接觸表面都是塑性變形，不再具有抵抗切向變形的能力，也就不再具有切向剛度。但是此時(shí)結(jié)合面的接觸壓力仍舊小于兩材料中較軟材料的硬度。當(dāng)微凸體頂端的變形量增大至ω=110ωec時(shí)，結(jié)合面接觸壓力和兩材料中較軟材料的硬度相等。此時(shí)結(jié)合面的微凸體進(jìn)入完全塑性階段。微凸體在彈塑性過渡變形階段中三個(gè)不同階段的法向接觸載荷分別為如下[15]。

1) 當(dāng)ωec≤ω<6ωec時(shí)：

(5)

2) 當(dāng)6ωec≤ω<68ωec與68ωec≤ω<110ωec時(shí)：

(6)

1.2 結(jié)合面法向接觸載荷模型

根據(jù)文獻(xiàn)[15]可知，結(jié)合面的法向接觸載荷模型如下。

1) 當(dāng)1

(7)

2) 當(dāng)1

(8)

(9)

2 結(jié)合面切向接觸載荷模型

由于接觸表面已經(jīng)發(fā)生塑性變形的微凸體，不能承受切向載荷，由此在建立切向接觸載荷模型時(shí)，不能包含整個(gè)接觸表面已經(jīng)發(fā)生完全塑性變形的彈塑性過渡中第三小階段以及完全塑性階段。

根據(jù)文獻(xiàn)[17]可得，結(jié)合面的單個(gè)微凸體所能承受的切向接觸載荷為：

(10)

式中：v為兩材料中較軟材料的Possion比，F(xiàn)為單個(gè)微凸體所承受的法向接觸載荷。

根據(jù)文獻(xiàn)[2]可知，結(jié)合面微接觸點(diǎn)的面積分布密度函數(shù)為：

(11)

所有處于完全彈性階段的微凸體所能承受的切向接觸載荷為：

(12)

將式(1)、(10)和(11)代入到式(12)得如下式子。

1)當(dāng)D≠1.5時(shí)：

(13)

2) 當(dāng)D=1.5時(shí)：

(14)

所有處于彈塑性過渡第一階段的微凸體所能承受的切向接觸載荷為：

(15)

將式(5)、(10)和(11)代入到式(15)，得：

(16)

所有處于彈塑性過渡第二階段的微凸體所能承受的切向接觸載荷為：

(17)

將式(6)、(10)和(11)代入到式(17),得：

(18)

3 結(jié)合面切向剛度模型

結(jié)合面的真實(shí)接觸面積為：

(19)

將式(11)代入式(19)，得：

(20)

單個(gè)微凸體與平面接觸時(shí)的切向接觸剛度可表示為[18]：

(21)

式中：μ為摩擦系數(shù)。

處于彈塑性過渡第三階段以及完全塑性階段的微凸體接觸表面已經(jīng)發(fā)生完全塑性變形，不再具有接觸剛度，由此結(jié)合面切向接觸剛度Kt可用下式來表示:

(22)

將式(11)和式(21)代入式(22)可得：

(23)

對(duì)式(23)進(jìn)行量綱一化：

(24)

4 分析結(jié)果的討論

選取當(dāng)量Possion比ν=0.2,摩擦系數(shù)μ=0.2，φ分別取1.0,0.1,0.01，量綱一分形粗糙度參數(shù)G*取1×10-8,1×10-9,1×10-10,1×10-11,1×10-12，而D=1.1～1.9。數(shù)字仿真計(jì)算結(jié)果見圖2～7。

對(duì)圖2～7進(jìn)行分析可得如下結(jié)果。

切向接觸剛度隨著D的變化出現(xiàn)先增大后減小的趨勢(shì)與之前的研究[13,19]結(jié)果相同，但是之前的研究結(jié)果并未給出達(dá)到最大切向接觸剛度對(duì)應(yīng)的D是隨著φ的增大呈現(xiàn)出減小的趨勢(shì)。

圖2 D對(duì)的影響 (φ=1,G*=1×10-9)

圖3 D與φ共同對(duì)的影響(G*=1×10-9)Fig.3 Influence of D and φ on (G*=1×10-9)

圖4 G*對(duì)的影響(φ=0.01)Fig.4 Influence of G* on (φ=0.01 )

圖5 φ對(duì)的影響 (G*=1×10-9)

圖6 h對(duì)的影響 (D=1.6,φ=0.01,G*=1×10-9)Fig.6 Influence of h on (D=1.6,φ=0.01,G*=1×10-9)

圖7 實(shí)際彈塑性接觸面積占實(shí)際接觸面積的比例圖Fig.7 Ratio of actual elastoplastic contact area to actual contact area

7) 微凸體處于彈塑性過渡階段的實(shí)際接觸面積占整個(gè)接觸過程的實(shí)際接觸面積的比例與D間的關(guān)系見圖7。該結(jié)果與文獻(xiàn)[6]的趨勢(shì)基本一致。由圖7可以看出，當(dāng)D較大或較小時(shí)，微凸體處于彈塑性過渡階段的實(shí)際接觸面積占整個(gè)接觸過程的實(shí)際接觸面積的比例較大，因此不能忽略彈塑性過渡階段對(duì)接觸特性的影響。

5 理論計(jì)算與試驗(yàn)驗(yàn)證

為驗(yàn)證本文建立的考慮彈塑性過渡階段的切向剛度模型的優(yōu)越性，設(shè)計(jì)如圖8所示的接觸試驗(yàn)?zāi)Ｐ汀?/p>

圖8 實(shí)驗(yàn)裝置圖Fig.8 Experimental device

為將結(jié)合部孤立，螺栓兩側(cè)采用橡膠墊將兩塊板與螺栓隔離開，并且將整個(gè)結(jié)構(gòu)與地面也通過橡膠墊隔離開。螺栓上粘有應(yīng)變片用來測(cè)量螺栓拉力，在平板上裝有IC壓電式加速度傳感器，其靈敏度為100 mv/g，量程為50 g，分辨率為0.000 2 g，采樣頻率范圍為0.7～10 000 Hz，采用阿爾泰USB2085動(dòng)態(tài)數(shù)據(jù)采集儀進(jìn)行實(shí)驗(yàn)數(shù)據(jù)的采集與記錄，該采集儀具有32通道并且可以進(jìn)行數(shù)模轉(zhuǎn)化，能夠連續(xù)地采集所需要的離散數(shù)據(jù)點(diǎn)的信息。

實(shí)驗(yàn)裝置見圖9，在邊界位置布置8個(gè)位移傳感器，其分布位置為圖9中1～8，并在內(nèi)部布置9個(gè)激振點(diǎn)，其分布位置為圖9中1～9。激振力信號(hào)通過信號(hào)放大器與采集卡反饋給計(jì)算機(jī)，并且利用位移傳感器測(cè)量輸出信號(hào)，通過采集信號(hào)的頻域分析，得到試樣的固有頻率。

圖9 結(jié)合面固有頻率測(cè)定實(shí)驗(yàn)裝置圖Fig.9 Experimental device for the measurement of the natural frequency of the joint

有限元模型見圖10，上下試樣之間的接觸層被等效為一種均質(zhì)的各向同性的虛擬材料，虛擬材料的等效厚度為309.6 μm。

圖10 包含虛擬接觸層的接觸試樣有限元模型Fig.10 Finite element model for contact sample containing virtual layer

有限元模型采用8節(jié)點(diǎn)的等參單元，單元總數(shù)為11 140個(gè)，其中虛擬接觸層的單元總數(shù)為852個(gè)。有限元模型的邊界條件為自由狀態(tài)，虛擬材料的法向參數(shù)通過文獻(xiàn)[8]的法向接觸剛度模型推導(dǎo)得出，而切向參數(shù)通過本文模型以及文獻(xiàn)[12]建立的切向剛度模型分別推導(dǎo)得出。將所獲得的虛擬接觸層材料參數(shù)通過APDL程序引入所建立的有限元模型，通過本文模型、文獻(xiàn)[12]模型得到的有限元仿真結(jié)果與試驗(yàn)所得到在不同擰緊力矩下的前三階的固有頻率的對(duì)比曲線，見圖11。

圖11 有限元模型固有頻率與實(shí)驗(yàn)固有頻率對(duì)比Fig.11 Comparison of the natural frequencies of FE model and the experimental natural frequencies

由圖11可以看出，不論虛擬材料結(jié)合部還是試驗(yàn)測(cè)量得到的固有頻率都是隨著擰緊力矩的增大而增大，這是因?yàn)榻Y(jié)合部接觸壓力增大，導(dǎo)致結(jié)合部接觸剛度增大，使得系統(tǒng)的固有頻率升高。

從圖11還可以看出，在不同階的固有頻率曲線中，本文所建立的模型得到的固有頻率總是比文獻(xiàn)[12]模型所建立的固有頻率更加地靠近試驗(yàn)測(cè)得的固有頻率曲線，這是因?yàn)楸疚乃⒌哪Ｐ突谖⑼贵w的連續(xù)變形理論，考慮彈塑性過渡階段產(chǎn)生的切向剛度，比文獻(xiàn)[12]所建立的模型更加地符合實(shí)際接觸情況，所以得到的數(shù)據(jù)更加準(zhǔn)確，精度也更高。因此本文提出的由切向接觸剛度模型來預(yù)測(cè)固定結(jié)合部的動(dòng)態(tài)特性會(huì)更加準(zhǔn)確。

6 結(jié) 論

1) 基于微凸體連續(xù)變形理論，首次考慮了將彈塑性過渡階段中屈服區(qū)域沒有擴(kuò)大至整個(gè)接觸表面前的彈塑性過渡階段所產(chǎn)生的切向剛度計(jì)入粗糙表面切向接觸剛度模型。在文獻(xiàn)[17]的基礎(chǔ)上，建立了包含微凸體在彈塑性過渡階段中前兩個(gè)小階段的切向剛度分形模型。

2) 分析結(jié)果顯示，隨著分形維數(shù)的增大，切向接觸剛度出現(xiàn)先增大后減小的趨勢(shì)。但是在不同的較軟材料的屈服強(qiáng)度與兩材料的復(fù)合彈性模量的比值下，出現(xiàn)最大切向接觸剛度對(duì)應(yīng)的分形維數(shù)是不同的。φ增大時(shí)使得切向接觸剛度發(fā)生改變的分形維數(shù)出現(xiàn)減小。

3) 仿真結(jié)果顯示：當(dāng)分形維數(shù)較大或較小時(shí)，彈塑性過渡階段的實(shí)際接觸面積占整個(gè)接觸過程的實(shí)際接觸面積的比例較大，不能忽略彈塑性過渡階段對(duì)接觸特性的影響。此外還發(fā)現(xiàn)彈塑性過渡階段產(chǎn)生的切向接觸剛度的大小取決于表面分形維數(shù)的大小。因?yàn)榇植诒砻鎻椝苄赃^渡階段的實(shí)際接觸面積占整個(gè)接觸過程的實(shí)際接觸面積的比例隨著分形維數(shù)的增加出現(xiàn)先減小后增大的趨勢(shì)，當(dāng)分形維數(shù)較大或較小時(shí)，忽略彈塑性過渡階段會(huì)使得切向接觸剛度過小。

4) 本文所建立的粗糙表面微凸體切向剛度分形模型與文獻(xiàn)[12]的模型進(jìn)行比較，與試驗(yàn)結(jié)果更為接近。說明了由于考慮了微凸體彈塑性過渡階段的彈性，使得建立的模型得到的切向剛度更加合理。