田紅亮 趙春華 朱大林 秦紅玲

（三峽大學(xué) 機械與材料學(xué)院，湖北 宜昌 443002）

1 文獻［1]的1個缺陷

傳統(tǒng)的結(jié)合部方法采用彈簧－阻尼器［2－6]（springdamper）模型，該模型能模擬結(jié)合部兩側(cè)零件的相對運動，但該模型的一個缺點是理論預(yù)測結(jié)果與實驗數(shù)據(jù)相差較大.

為克服彈簧－阻尼器模型的缺陷，進一步提高理論預(yù)測精度，田紅亮和李斌等［1]采用如圖1所示的固體1維表面形貌首次提出機械結(jié)構(gòu)固定結(jié)合部虛擬材料的動力學(xué)建模方法，將2種配對材料的微觀接觸凹凸部分看成一種虛擬的各向同性材料（認為無論沿任何方向，固體的力學(xué)性能都是相同的），首先，通過一些理論的方法得到該虛擬材料的彈性模量、Poisson比、厚度和密度；然后，將這些參數(shù)嵌入到有限元軟件中，可將整個機械結(jié)構(gòu)的動力學(xué)特性通過分形維數(shù)D（1＜D＜2）、分形粗糙度G、近似理論模態(tài)等參數(shù)表現(xiàn)出來；最后，為驗證虛擬材料參數(shù)是否具有工程實際可操作性及是否有效，可將近似理論模態(tài)與實驗?zāi)B(tài)進行比較，例如將兩者的相似振型進行定性比較，將兩者相應(yīng)的固有頻率、阻尼比進行定量比較，給出相對誤差.因此該方法存在以下1個缺陷：虛擬的各向同性材料假設(shè)不能完全反映實際工程材料的特征，事實上，實際工程材料大多為各向異性（沿不同方向力學(xué)性能不同的材料）.

圖1 固體表面形貌

本文分析三維各向異性分形表面特征，將分形維數(shù)推廣到三維的一般情況（2＜D＜3），區(qū)分彈、塑性接觸時單峰的實際接觸面積（都用a表示）.

2 三維各向異性分形表面特征

更起作用通用的單變量標量WM函數(shù)模擬加權(quán)、隨機重疊的隆起部狀表面.單變量標量 WM函數(shù)［7]為

一個三維分形表面在各平面方向顯示隨機性，其高度函數(shù)可用式（1）的實部表示

式中，ρ和θ為平面極坐標；γ＞1為決定輪廓空間頻率密度的參數(shù)；M為構(gòu)造表面重疊隆起部的個數(shù)；Am為控制表面幾何各向異性的量值為與采樣長度L有關(guān)的波數(shù)；n為頻率指數(shù)；2＜D＜3為表面的分形維數(shù)；αm為在方位角的方向偏置隆起部的任意角度；Φm，n為均勻分布的隨機相位，且

將式（3）～式（7）代入式（2），得

式中，G為分形粗糙度；nmax為與6晶格距離截止長度Ls有關(guān)的頻率指數(shù)上限，且

式（8）的表面經(jīng)驗功率譜對應(yīng)的等高線示意圖如圖2所示，其中kx、ky分別為沿x、y軸方向的空間頻率，單位是m－1.可以看出，功率譜的等高線明顯表明具有8重對稱性，盡管下文的圖3（a）、圖3（b）都無特定的方位角結(jié)構(gòu).

圖2 表面經(jīng)驗功率譜的等高線

當γ＝1.5、M＝10、Φm，n＝0、G＝1.36×10－5μm時，式（8）模擬各向異性表面如圖3所示.

圖3 WM函數(shù)模擬隆起部狀表面

可以看出，分形維數(shù)D越大，表面形貌越光滑.圖3（b）中表面高、低頻率組分的幅值比大于圖3（a）中表面的情況.

當M＝1，則m＝1時，式（8）退化為

橫截微接觸點的半徑為r′，單峰波形式的最長波長為2r′［8].類比式（9），可得基波長單峰對應(yīng)的頻率指數(shù)為

頻率指數(shù)式（11）在式（10）中對應(yīng)的一個余弦項為

當Φ1，n0＝0、x＝r′時，由式（12）可得單峰干涉量［9]為

等效球體單峰與剛性平面的彈性接觸如圖4所示.

圖4 等效球體單峰與剛性平面的彈性接觸

對圖4使用勾股定理，得

式中，R為單峰的曲率半徑.

當δ?2R［10]時，可得下列近似式

將式（13）代入式（15），得

進一步可改寫為

式中，a′＝πr′2為微接觸點的橫截面積.

Hertz接觸理論預(yù)測微接觸點的平均壓應(yīng)力［11]為

施加在一個微接觸點上的彈性力為

式中，a＝πr2為微接觸點的實際接觸面積；r為微接觸點實際接觸面積的半徑.

將以下關(guān)系式

代入式（19），得

將式（17）代入式（21），得

完全塑性流時的準則［12]為

式中，b為待定無量綱常數(shù).

將式（17）代入式（23），得

式（24）除以式（25），得

式中，a′c為臨界微接觸點面積，且

從式（26）可以看出，當a′＞a′c即δ＜δR時，微凸體彈性變形；當a′≤a′c即δ≥δR時，微凸體完全塑性變形，此時施加在一個完全塑性微接觸點上的力為

Pullen和 Williamson［13]用一剛性平面壓縮一粗糙金屬表面，該實驗證明：微接觸點上的塑性流使非接觸區(qū)域表面高度均勻增加.因為對于小和中等接觸壓應(yīng)力，總的實際接觸面積比表觀接觸面積小得多，a略微大于a′.CEB模型［14]預(yù)測一個塑性變形微接觸點的接觸面積為a＝πRδ（2－δR／δ），當δ?δR并利用式（15）時，a≤2πRδ＝πr′2＝a′，于是式（28）可近似為

由式（22）可得臨界微接觸點面積a′c對應(yīng)的彈性力為

由式（29）可得臨界微接觸點面積a′c對應(yīng)的塑性力為

兩表面接近時，假設(shè)施加在一個微接觸點上的法向力連續(xù)變化，即令 ΔFe（a′c）＝ΔFp（a′c），得

將式（27）代入式（32），得

將式（33）代入式（27），得

由式（19）得

由式（29）得

圖5 無量綱應(yīng)力－應(yīng)變關(guān)系

Mandelbrot提出［15]，大于特定面積s′的面積的島嶼個數(shù)N服從冪定律關(guān)系

式中，Ds為海岸線的分形維數(shù)，且

海洋水平線上的島嶼與粗糙表面上的橫截微凸體具有緊密的分形相似，類比式（37）可得粗糙表面上的橫截微凸體個數(shù)為

由式（39）可得粗糙表面上的橫截微凸體尺度分布函數(shù)為

將式（38）代入式（40），得

3 彈塑性接觸分析

當a′L＞a′c時，彈性、完全塑性微接觸點都存在，總的彈性變形力、總的塑性變形力、總的變形力、總的實際彈性接觸面積、總的實際塑性接觸面積、總的實際接觸面積分別為

將式（22）、式（41）代入式（42），得

進一步定積分，得

將式（29）、式（41）代入式（43），得

將式（49）～式（51）代入式（44），得

4 塑性接觸分析

等效粗糙表面的總橫截面積為

當a′L≤a′c時，所有微凸體都處于完全塑性變形，總的彈性變形力、總的塑性變形力、總的變形力、總的實際彈性接觸面積、總的實際塑性接觸面積、總的實際接觸面積分別為

5 結(jié) 論

分析三維各向異性分形表面特征，將分形維數(shù)推廣到三維的一般情況（2＜D＜3），給出彈、塑性接觸時單峰的實際接觸面積（都用a表示）與單峰的橫截面積a′之間的關(guān)系.

［1]Tian Hongliang，Li Bin，Liu Hongqi，et al.A New Method of Virtual Material Hypothesis－Based Dynamic Modeling on Fixed Joint Interface in Machine Tools［J].ELSEVIER International Journal of Machine Tools ＆Manufacture，2011，51（3）：239－249.

［2]Ibrahim R A，Pettit C L.Uncertainties and Dynamic Problems of Bolted Joints and Other Fasteners［J].Journal of Sound and Vibration，2005，21（3－5）：857－936.

［3]Jiang Shuyun，Zheng Yunjian，Zhu Hua.A Contact Stiffness Model Machined Plane Joint Based on Fractal Theory［J].ASME Journal of Tribology，2010，132（1）：011401－1－011401－7.

［4]劉 恒，劉 意，王為民.接觸界面法向剛度等效的新方法［J].機械工程學(xué)報，2011，47（17）：37－43.

［5]蘭國生，張學(xué)良，丁紅欽，等.基于分形理論的結(jié)合面改進接觸模型［J].農(nóng)業(yè)機械學(xué)報，2011，42（10）：217－223，229.

［6]尤晉閩，陳天寧.基于分形接觸理論的結(jié)合面法向接觸參數(shù)預(yù)估［J].上海交通大學(xué)學(xué)報，2011，45（9）：1275－1280.

［7]Ausloos M，Berman D H.A Multivariate Weierstrass－Mandelbrot Function［J].Proceedings of the Royal Society of London，1985，Series A Mathematical and Physical Sciences 400（1819）：331－350.

［8]Majumdar A，Bhushan Bharat.Role of Fractal Geometry in Roughness Characterization and Contact Mechanics of Surfaces［J].ASME Journal of Tribology，1990，112（2）：205－216.

［9]Yan W，Komvopoulos Kyriakos.Contact Analysis of E－lastic－Plastic Fractal Surfaces［J].Journal of Applied Physics，1998，84（7）：3617－3624.

［10]田紅亮，朱大林，秦紅玲.固定接觸界面法向靜彈性剛度［J].應(yīng)用力學(xué)學(xué)報，2011，28（3）：318－322.

［11]田紅亮，朱大林，秦紅玲.兩彈塑性非赫茲接觸粗糙表面溫升的分形模型［J].三峽大學(xué)學(xué)報：自然科學(xué)版，2010，32（4）：77－83.

［12]Greenwood J A，Williamson J B P.Contact of Nominally Flat Surfaces［J].Proceedings of the Royal Society of London，1966，Series A Mathematical and Physical Sciences 295（1442）：300－319.

［13]Pullen J，Williamson J B P.On the Plastic Contact of Rough Surfaces［J].Proceedings of the Royal Society of London，1972，Series A Mathematical and Physical Sciences 327（1569）：159－173.

［14]Chang W R，Etsion I，Bogy David B.An Elastic－Plastic Model for the Contact of Rough Surfaces［J].ASME Journal of Tribology，1987，109（2）：257－263.

［15]Peitgen Heinz－Otto.Beno?t B Mandelbrot（1924－2010）［J].Science，2010，330（6006）：926.