李 玲，王晶晶，史小輝，楊亞蘭，阮曉光

(1.西安建筑科技大學(xué)機(jī)電工程學(xué)院，西安 710055；2.西安陜鼓動(dòng)力股份有限公司，西安 710075)

結(jié)合面在機(jī)床運(yùn)行過程中起著傳遞運(yùn)動(dòng)、載荷和能量的重要作用，同時(shí)也是引起機(jī)床非線性和不確定性的根源[1?3]。研究表明，結(jié)合面的剛度約占機(jī)床中總剛度的60%～80%，其剛度值的大小和變化規(guī)律，直接影響著機(jī)床的靜、動(dòng)態(tài)力學(xué)性能和精度保持性水平[4?5]。因此，開展粗糙表面接觸剛度特性的研究，對(duì)整機(jī)性能的評(píng)估與動(dòng)態(tài)優(yōu)化具有重要的理論指導(dǎo)意義。

Greenwood和Williamson[6]是較早開展接觸理論研究的先驅(qū)者。他們將粗糙表面上微凸體高度分布看成隨機(jī)函數(shù)，并基于Hertz理論提出統(tǒng)計(jì)學(xué)微觀彈性接觸模型(GW模型)，奠定了接觸力學(xué)基礎(chǔ)。為進(jìn)一步完善粗糙表面的接觸理論，Chang等[7]考慮微凸體變形過程中體積守恒原則，建立了粗糙表面的彈-塑性接觸模型(CEB模型)。但該模型卻忽略了微凸體由彈性到塑性變形的中間過程，使得在臨界屈服點(diǎn)處接觸載荷出現(xiàn)不連續(xù)現(xiàn)象。Zhao等[8]為彌補(bǔ)CEB模型的不足，將微凸體的變形分為彈性變形、彈塑性變形和完全塑性變形階段，并利用多項(xiàng)式數(shù)學(xué)插值法建立了粗糙表面的彈塑性接觸模型(ZMC模型)。該模型增強(qiáng)了接觸載荷和平均接觸壓力在彈性屈服點(diǎn)與塑性臨界接觸點(diǎn)的連續(xù)性，但卻忽略了光滑性，導(dǎo)致建立的接觸剛度模型在2個(gè)臨界點(diǎn)處均是不連續(xù)的。Kogut和Estion[9?10]利用有限元法，將微凸體的彈塑性接觸變形區(qū)細(xì)分為2個(gè)子區(qū)域，建立出更高精度的彈塑性接觸模型(KE模型)。但在這兩個(gè)子區(qū)域的轉(zhuǎn)換點(diǎn)處和彈塑性變形的開始與結(jié)束點(diǎn)處，接觸載荷也是不連續(xù)的，從而引起其剛度模型的不連續(xù)性。Brake[11]采用三次Hermite多項(xiàng)式插值函數(shù)，提出了一種在不同接觸狀態(tài)之間，接觸載荷光滑且連續(xù)變化的粗糙表面彈塑性接觸模型，但由于其多項(xiàng)式階次較高，獲得的接觸載荷與變形量之間的關(guān)系曲線出現(xiàn)振蕩，從而導(dǎo)致使用其所建立的接觸剛度模型也會(huì)出現(xiàn)振蕩，甚至出現(xiàn)負(fù)剛度的情況。

上述模型除了各自所存在的缺陷外，還忽略了微凸體所附著的基體的變形情況。對(duì)于一些含有硬涂層的粗糙表面，其接觸特性主要由軟基體上的硬涂層材料決定，此時(shí)較軟基體的變形對(duì)粗糙表面的接觸特性有著顯著的影響[1,12 ?13]。在之前的研究中[1]，僅對(duì)純彈性接觸下的基體變形進(jìn)行了分析，但未包含微凸體的彈塑性和完全塑性變形階段。為此，本文基于微凸體在彈塑性接觸變形過程中，其接觸載荷和接觸剛度具有連續(xù)且光滑的特性，首先對(duì)Hermite多項(xiàng)式直接插值函數(shù)進(jìn)行了改進(jìn)，提出了一種在微凸體的整個(gè)變形過程中，其接觸剛度滿足連續(xù)且光滑變化的模型。其次，引入了基體的彈性變形，提高了傳統(tǒng)接觸剛度建模的精度。然后，基于統(tǒng)計(jì)學(xué)方法，將單個(gè)微凸體-基體系統(tǒng)模型拓展到整個(gè)粗糙表面上，建立了考慮基體變形的連續(xù)光滑結(jié)合面剛度模型。最后，與GW、ZMC、KE、Brake模型進(jìn)行對(duì)比分析，驗(yàn)證了本文接觸剛度模型的正確性。

1 單個(gè)微凸體-基體系統(tǒng)的接觸剛度

基于以下4點(diǎn)假設(shè)開展結(jié)合面剛度建模研究。

1)將圖1(a)兩個(gè)粗糙表面的接觸，等效為圖1(b)所示的光滑剛性平面與等效粗糙表面的接觸。

圖1 粗糙表面接觸模型Fig.1 Rough surface contact model

2)微凸體發(fā)生接觸變形前后的體積保持不變，其變形剛度為ka；基體的變形為彈性變形，其變形剛度為kb，如圖1(c)所示。

3)微凸體的分布是各向同性的，且各個(gè)微凸體間的變形相互獨(dú)立。

4)微凸體高度服從高斯分布，所有微凸體的峰頂都等效成半球體，且具有相同的峰頂曲率半徑。

圖1(c)為單個(gè)微凸體-基體系統(tǒng)的受力分析。當(dāng)剛性平面被施加一法向載荷f時(shí)，微凸體產(chǎn)生δa的變形，同時(shí)基體也發(fā)生δb宏觀變形，z為微凸體變形前的高度，d為基體平均平面與剛性平面間的距離。考慮基體變形影響下系統(tǒng)總變形量δ可表示為：

基于Chang等[7]的研究，微凸體由彈性變形轉(zhuǎn)變?yōu)閺椝苄宰冃螘r(shí)的臨界接觸變形量可表示為：

微凸體開始發(fā)生完全塑性變形時(shí)的臨界變形量δp與初始屈服點(diǎn)臨界變形量δc的關(guān)系為[15]：

基于式(1)～式(3)，在考慮基體變形的同時(shí)，對(duì)微凸體的彈性、彈塑性以及完全塑性接觸變形進(jìn)行逐一分析。為便于描述，下文依次按彈性變形、塑性變形和彈塑性變形3個(gè)階段建模分析。

1.1 彈性接觸變形階段

當(dāng)微凸體變形量0＜δa＜δc時(shí)，根據(jù)Hertz接觸理論，微凸體的彈性接觸載荷fae可表示為：

式中，δae為微凸體彈性階段的接觸變形。由剛度的定義，可得單個(gè)微凸體發(fā)生彈性變形時(shí)的接觸剛度：

當(dāng)微凸體受到載荷時(shí)，微凸體和基體的接觸區(qū)是一個(gè)半徑為rb的圓形區(qū)域(如圖1(c)所示)，且在該區(qū)域內(nèi)各點(diǎn)受到的壓力各不相同，最大接觸壓力P0發(fā)生在接觸中心，可表示為[1]：

距離接觸區(qū)域中心的距離r不同，基體表面發(fā)生的法向位移不同，當(dāng)基體受到的最大壓力P0不斷變化時(shí)，基體表面的法向位移可表示為：

聯(lián)立方程式(6)和式(8)，可獲得基體的彈性剛度：

根據(jù)力的相互作用原理和疊加原理，可得：

將式(4)、式(9)和式(11)代入式(10)，可得微凸體變形δae和系統(tǒng)總變形δ之間的關(guān)系：

由式(11)可知，微凸體變形δae是系統(tǒng)總變形δ的函數(shù)，令δae=g(δ)，則式(13)可重新表示為：

為獲得δae與δ的直接關(guān)系，以便下文進(jìn)行基體變形影響下粗糙表面的接觸建模與分析，在之前的工作中已采用不動(dòng)點(diǎn)迭代法對(duì)式(14)進(jìn)行了數(shù)值求解，獲得彈性接觸時(shí)，考慮基體變形影響下微凸體接觸變形δae的計(jì)算表達(dá)式[16]：

1.2 完全塑性接觸變形階段

當(dāng)微凸體變形量δa＞δp時(shí)，微凸體將發(fā)生完全塑性接觸變形，此時(shí)微凸體的接觸載荷fap和接觸剛度kap可分別表示為[17]：

式中，δap表示完全塑性階段的微凸體接觸變形。同理，根據(jù)力的相互作用原理，可得：

聯(lián)立式(9)、式(11)、式(17)和式(18)，可得微凸體發(fā)生完全塑性變形時(shí)的變形量δap與系統(tǒng)總變形δ之間的關(guān)系為：

1.3 混合彈塑性接觸變形階段

當(dāng)微凸體變形量δc≤δa＜δp時(shí)，微凸體變形處于混合彈塑性階段。通過研究發(fā)現(xiàn)，改進(jìn)Brake模型[11]的Hermite多項(xiàng)式直接插值函數(shù)，來彌補(bǔ)插值函數(shù)階次較高而引起微凸體接觸參數(shù)振蕩行為，可實(shí)現(xiàn)微凸體的接觸載荷連續(xù)光滑變化。利用對(duì)數(shù)變換可將混合彈塑性區(qū)域插值間區(qū)間(1，110)減小為(0，4.7≈ln110)，進(jìn)而很大程度降低了接觸行為的振蕩性，最終可得彈塑性接觸變形階段的法向接觸載荷和法向接觸剛度分別為：

式中：δaep為混合彈塑性階段的微凸體接觸變形；C1、C2和C3為Hermite多項(xiàng)式插值函數(shù)的系數(shù)項(xiàng)，分別表示為：其中，q=δp/δc=110，表示塑性臨界接觸變形與彈性臨界變形之比。

同理，在微凸體的彈塑性接觸變形階段中引入基體變形的影響。根據(jù)力的相互作用原理，有：

聯(lián)立式(9)、式(11)、式(20)和式(22)，可得：

為簡(jiǎn)化式(23)的數(shù)值求解過程，忽略參數(shù)值較小的C3項(xiàng)，式(23)可重新表示為：

進(jìn)一步化簡(jiǎn)為：

式中，λ=R/4rb為確定微凸體幾何形狀的基本參數(shù)。

類似地，令δaep=g(δ)，則式(25)可表示為：

采用不動(dòng)點(diǎn)迭代法可得式(26)的近似解，迭代公式為：

式中，n=0,1,···,∞，迭代的次數(shù)越多，得到的結(jié)果越精確。經(jīng)過數(shù)值計(jì)算發(fā)現(xiàn)，只需迭代2次得到的近似解就非常接近于數(shù)值解，且兩者之間的誤差不超過5%。因此微凸體的變形量δa可以近似地表示為一個(gè)關(guān)于總變形δ的函數(shù)，即：

圖2顯示了微凸體變形量δa數(shù)值解與近似解的對(duì)比。由圖可知，對(duì)于不同的微凸體幾何參數(shù)λ，數(shù)值解與近似解之間的誤差很小，從而證明用近似解代替數(shù)值解的有效性。從圖中可以發(fā)現(xiàn)，對(duì)于相同的系統(tǒng)變形量δ，微凸體變形量δa隨著λ的增大而減小，而基體的變形逐漸增大。此時(shí)，基體的接觸變形在微凸體-基體系統(tǒng)中起主導(dǎo)作用。

圖2 微凸體變形量數(shù)值解與近似解的對(duì)比Fig.2 Comparison of numerical solution and approximate solution of asperity deformation

2 考慮基體變形的接觸剛度模型

為研究結(jié)合面的微觀接觸特性，需將單個(gè)微凸體-基體系統(tǒng)模型擴(kuò)展到整個(gè)粗糙表面上，從而建立一種具有連續(xù)光滑特性并考慮基體變形的粗糙表面接觸新模型。根據(jù)統(tǒng)計(jì)學(xué)方法和微凸體高度服從高斯分布，假設(shè)粗糙表面名義接觸面積上微凸體的總數(shù)為N，則不同接觸間隙d下發(fā)生接觸的微凸體數(shù)量可表示為[14]：

式中：N=ηAn，η為微凸體密度，An表示名義接觸面積；φ(z)為微凸體高度分布的概率密度函數(shù)。對(duì)于某一給定的表面接觸間隙d，整個(gè)變形區(qū)內(nèi)粗糙表面的接觸載荷Fnew可表示為：

式中，F(xiàn)e、Fp、Fep分別表示彈性接觸載荷、塑性接觸載荷及彈塑性接觸載荷，分別表示如下：

將式(4)和式(15)代入式(32)，式(16)和式(19)代入式(33)，式(20)和式(28)代入式(34)，即可分別獲得基體變形影響下粗糙表面的彈性接觸載荷、塑性接觸載荷及彈塑性接觸載荷，從而將式(31)進(jìn)一步表示為：

在微凸體的三個(gè)變形階段中，考慮基體變形δb影響的無量綱微凸體接觸變形可根據(jù)式(15)、式(19)和式(28)分別表示為：

無量綱概率密度函數(shù)φ(z)，可表示為：

式中，z*=z/σ，表示無量綱的微凸體高度參數(shù)，σs為微凸體高度的標(biāo)準(zhǔn)差，與粗糙表面高度的標(biāo)準(zhǔn)差σ的關(guān)系為[14]：

式中，γ=ηRσ為粗糙表面的綜合粗糙度參數(shù)。

同理，整個(gè)粗糙表面的法向接觸剛度Knew可表示為：

類似地，將式(5)、式(15)、式(17)、式(19)、式(21)和式(28)分別代入式(40)即可獲得考慮基體變形影響下的粗糙表面接觸剛度，則有：

3 分析與討論

3.1 接觸剛度模型對(duì)比分析

表1 四種不同粗糙表面形貌參數(shù)[18]Table1 Four different rough surface topography parameters

利用式(34)和式(38)，可得不同粗糙表面下接觸剛度與接觸載荷之間的關(guān)系，如圖3所示。由圖3可知，本文模型與GW、ZMC、KE、Brake模型計(jì)算的接觸剛度曲線均一致，且均隨著接觸壓力的增大而增大，從而表明本文剛度模型的合理性，可適用于不同粗糙形貌表面的接觸剛度特性分析。

通過對(duì)比圖3(a)～圖3(d)可知：隨著表面粗糙度的增大，不同接觸模型在相同載荷下的接觸剛度增大，會(huì)引起粗糙表面間真實(shí)接觸面積減少，從而使得粗糙表面的接觸剛度增大。因此，可以通過降低表面粗糙度來改善粗糙表面的接觸剛度特性。

此外，從圖3還發(fā)現(xiàn)：①對(duì)于光滑接觸表面，Brake模型預(yù)測(cè)的接觸剛度最大，且不同模型之間的差距較小，尤其是當(dāng)接觸壓力較小時(shí)，它們之間的差距很小；②對(duì)于中等粗糙接觸表面，依然是Brake模型的結(jié)果最大，但此時(shí)本文未考慮基體變形的模型最接近于GW模型結(jié)果；③對(duì)于粗糙和非常粗糙的表面，GW模型預(yù)測(cè)的接觸剛度最大，且其與其他模型的計(jì)算結(jié)果存在較大的差異。

圖3 不同粗糙表面模型下接觸剛度與法向接觸壓力的關(guān)系Fig.3 Relationship between contact stiffness and normal contact pressure for different rough surface models

3.2 基體變形對(duì)接觸剛度的影響

圖4顯示了不同粗糙表面形貌下基體變形對(duì)表面接觸剛度的影響。由圖可知，對(duì)于相同的無量綱表面接觸間隙d*，與未考慮基體變形下的接觸剛度值相比，考慮基體變形的接觸剛度較小。這是因?yàn)?，在未考慮基體變形的接觸模型中，認(rèn)為微凸體的變形等于系統(tǒng)的總變形，此時(shí)獲得的微凸體變形量不是其真實(shí)的變形量，包含基體的變形量。換言之，未考慮基體的宏觀變形時(shí)，計(jì)算得到的微凸體變形量大于其實(shí)際的變形，從而導(dǎo)致計(jì)算的接觸剛度偏大，而考慮基體變形時(shí)預(yù)測(cè)的接觸剛度更準(zhǔn)確。

圖4 不同粗糙表面形貌下基體變形對(duì)接觸剛度的影響Fig.4 The influence of matrix deformation with different rough surface morphologies on contact stiffness

對(duì)比圖4(a)～圖4(d)可知，隨著表面粗糙度的增大，兩者之間的差距逐漸減小，即基體變形的影響減弱了，直到表面非常粗糙時(shí)，兩者之間的差距很小，這與實(shí)際情況是一致的。從圖中還可以發(fā)現(xiàn)，隨著表面無量綱接觸間隙的減小，兩者之間的差距遞增，這是因?yàn)楸砻娼佑|間隙的減小意味著接觸載荷的增大，此時(shí)基體變形越大。

此外，從圖4中還可以獲得，表面接觸剛度隨表面間隙的增大而減小，且其減小的速率呈先快后慢的變化趨勢(shì)，這是因?yàn)楫?dāng)表面粗糙度與外載荷一定時(shí)，隨著表面接觸間隙的增大，表面間發(fā)生接觸微凸體的數(shù)目減小，引起粗糙表面真實(shí)接觸面積的減小，從而導(dǎo)致接觸剛度的減小。此外，當(dāng)無量綱表面接觸間隙d*=d/σ增大到3時(shí)，此時(shí)粗糙表面的接觸剛度非常小，趨近于0，這是因?yàn)閷?duì)于微凸體高度服從高斯分布的粗糙表面，微凸體高度分布在[?3σ，3σ]的概率為99.73%，即幾乎所有的微凸體高度z均小于3σ，所以當(dāng)d*=3時(shí)，粗糙表面間接觸的微凸體數(shù)目極少，從而導(dǎo)致表面接觸剛度很小，幾乎趨近于0。

3.3 基體變形對(duì)接觸載荷的影響

圖5表明，不同粗糙表面形貌下，接觸載荷與接觸間隙之間的關(guān)系，即隨著表明粗糙度的增加，接觸載荷在增加，而考慮基體變形和未考慮基體變形差距在減小，且隨著表面粗糙度的增大，基體變形的影響逐漸減小，其產(chǎn)生機(jī)理與3.2中接觸剛度與接觸間隙之間的規(guī)律相同。

綜上所述，基體變形直接影響結(jié)合面的接觸特性，尤其在粗糙度較小時(shí)，對(duì)結(jié)合面的接觸特性影響很大，需要重點(diǎn)考慮。因此，對(duì)于一些含有硬涂層的粗糙表面，其接觸特性主要由軟基體上的硬涂層材料決定，此時(shí)較軟基體的變形對(duì)粗糙表面的接觸特性有著顯著的影響。

4 結(jié)論

本文為提高了結(jié)合面?zhèn)鹘y(tǒng)接觸剛度建模的精度，引入了基體的彈性變形，并將單個(gè)微凸體-基體系統(tǒng)模型拓展到整個(gè)粗糙表面上，建立了考慮基體變形的連續(xù)光滑結(jié)合面剛度模型，通過與GW、ZMC、KE、Brake模型進(jìn)行對(duì)比分析，驗(yàn)證了本文剛度模型的正確性。主要結(jié)論如下：

(1)通過對(duì)比分析發(fā)現(xiàn)，針對(duì)相對(duì)光滑表面，GW、ZMC、KE、Brake模型與本文模型所得的接觸剛度之間的差異較小。而當(dāng)表面粗糙度增大時(shí)，GW模型與其他幾種模型的差異逐漸增大，此時(shí)本文模型結(jié)果介于Brake模型和ZMC模型之間，表明本文模型的正確性。

(2)考慮基體變形時(shí)所獲得的接觸剛度和接觸載荷比忽略基體變形時(shí)的要小，且隨著表面粗糙度的減小，基體變形的影響快速增加。因此，當(dāng)表面粗糙度較小時(shí)，基體變形直接影響結(jié)合面的接觸特性不容忽略。

(3)接觸載荷與表面粗糙度是影響結(jié)合面接觸剛度的兩個(gè)主要因素，隨著載荷的增大或表面粗糙度的減小，接觸剛度隨之遞增。因此，實(shí)際工程中可通過優(yōu)化接觸載荷與表面粗糙度參數(shù)，來匹配設(shè)計(jì)結(jié)合面的接觸剛度和接觸阻尼值，獲得性能更加優(yōu)越動(dòng)態(tài)力學(xué)性能。