馮燕,俞小莉,劉震濤

?

應用熱彈塑性理論的結(jié)合部法向載荷分形模型

馮燕1,2,俞小莉1,劉震濤1

粗糙表面接觸特性參數(shù)是影響機械產(chǎn)品靜、動態(tài)特性的關(guān)鍵因素,而法向接觸剛度與法向接觸載荷是機械結(jié)構(gòu)靜態(tài)特性研究的重要內(nèi)容,且法向載荷對接觸剛度具有明顯影響[1]。早先基于實驗統(tǒng)計學參數(shù)建立的接觸模型,由于基礎(chǔ)數(shù)據(jù)受測試儀器和樣本限制,其計算結(jié)果具有不確定性[2]。工程表面大多具有分形特征,應用分形理論對結(jié)合部的接觸行為進行研究,是接觸理論發(fā)展的必然趨勢[3]。針對結(jié)合部法向接觸載荷,Goerke[4]、溫淑花[5]、Jiang[6]、蘭國生[7]、宿月文[8]、田紅亮[9-10]等學者先后進行了一些研究。但是,這些研究都局限于純機械領(lǐng)域,沒有考慮結(jié)合部溫度及熱應力變化帶來的影響,而實際工作中,以內(nèi)燃機為代表的機械產(chǎn)品往往處于變化著的熱環(huán)境中,溫度變化對接觸情況的影響是不容忽視的。

目前,溫度因素對彈塑性接觸力學行為的影響仍缺乏研究。劉天祥等建立了考慮屈服點溫度相關(guān)效應的粗糙表面熱彈塑性接觸無網(wǎng)格法數(shù)值計算模型,但該研究受限于儀器分辨率和采樣長度[11]。

綜合考慮以上情況,本文基于各向異性分形幾何理論的法向接觸載荷模型[9](以下簡稱TZQ模型),將其與熱彈塑性接觸理論建立的粗糙表面熱力學特性相結(jié)合,引入機械載荷比例系數(shù),建立了一種新的熱彈塑性接觸法向載荷分形模型,并分析了結(jié)合部溫差對熱彈塑性接觸法向載荷的影響規(guī)律。

1 建立模型

機械結(jié)構(gòu)的固定結(jié)合部問題從本質(zhì)上來說是兩個粗糙表面的接觸問題,可簡化為粗糙表面與剛性平面的接觸[6],如圖1所示。從微觀上說,粗糙表面存在大量微凸體,接觸僅發(fā)生在微凸體上。因此,研究結(jié)合部整體的法向載荷問題,需要以單個微凸體的熱彈塑性接觸研究為基礎(chǔ)。

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

1.1 單個微凸體熱彈塑性接觸的法向載荷

機械零部件之間一般通過螺栓、螺釘、銷釘、鍵等方式連接,當所處環(huán)境溫度發(fā)生變化或接觸的兩個零部件之間受熱狀況改變時,由于受原本連接的約束限制,在兩個接觸面間就會產(chǎn)生熱應力。熱應力問題與無熱載荷的應力分析問題相比,除了增加一項以初應變形式出現(xiàn)的溫度載荷外,其他完全相同。所以,可以將單個微凸體承受的法向載荷p分為機械載荷pm和由熱應力引起的載荷pt,即p=pm+pt。

在彈性變形階段,單個微凸體承受的法向彈性載荷

(1)

式中:由熱應力引起的法向載荷pte=σae,其中σ為由熱應力引起的單個微凸體微接觸截面上所承受的壓應力,ae為彈性變形階段單個微凸體的實際微接觸截面積。在接觸過程中,根據(jù)熱應力概念,σ=αE′ΔT[12],所以

(2)

(3)

根據(jù)文獻[9],單個微凸體承受的法向彈性機械載荷

(4)

將式(2)～式(4)代入式(1),得到考慮熱應力時單個微凸體承受的法向彈性載荷

(5)

在塑性變形階段,單個微凸體承受的法向塑性載荷

(6)

同理,由熱應力引起的法向載荷

(7)

式中:ap為塑性變形階段單個微凸體的實際微接觸截面積,此時ap與單個微凸體微接觸截面積a′的關(guān)系為[9]

(8)

由文獻[9]可知,單個微凸體承受的法向塑性機械載荷

(9)

式中:K為硬度與屈服強度之比,稱為相關(guān)因子,K=H/σy,其中H為較軟材料的硬度,σy為較軟材料的屈服強度。

將式(7)～式(9)代入式(6),得到考慮熱應力時單個微凸體承受的法向塑性載荷

(10)

為給出結(jié)合部熱彈塑性接觸問題法向載荷分形模型的另一種更簡潔的表達形式,引入機械載荷比例系數(shù)。若pte≠0,ptp≠0,設(shè)pme=f1αE′ΔTa′,pmp=f2αE′ΔTa′,其中f1為彈性變形階段機械載荷相比于熱載荷的比例系數(shù)(f1≠0),f2為塑性變形階段機械載荷相比于熱載荷的比例系數(shù)(f2≠0),比例系數(shù)越大則機械載荷占總載荷的比值越大。需要說明的是,f1和f2分別為彈、塑性變形階段某一特定情況下的比例系數(shù),其值可能會隨著工況變化而改變。包含機械載荷比例系數(shù)的單個微凸體承受的法向彈、塑性載荷分別為

(11)

(12)

1.2 整體熱彈塑性接觸的法向總載荷

微凸體微接觸截面積的分布函數(shù)為[9]

(13)

1

(14)

其定義域不同于文獻[9]中式(46)的定義域,因為D=1時,圖形不具備分形特征;D=2時,Ψ可取任意值,沒有意義。圖2為給定D的取值、根據(jù)信賴域反射算法繪制出的D-Ψ曲線。

圖2 D-Ψ曲線

(15)

將式(5)、(10)、(13)代入式(15),可求得整體熱彈塑性接觸的法向總載荷。

當D≠1.5時

(16)

式中

當D=1.5時

(17)

式中

將式(11)～式(13)重新代入式(15),得到考慮機械載荷比例系數(shù)的整體熱彈塑性接觸法向總載荷的表達式

(18)

需要強調(diào)的是:式(16)、(17)和式(18)是整體熱彈塑性接觸法向總載荷的兩種不同的表達形式,與文獻[5,7,9]中模型的最大不同之處在于,式(16)、(17)和式(18)中包含了結(jié)合部的溫度變化量ΔT,能反映溫度變化對接觸部位整體法向總載荷的影響。

當α=0或ΔT=0時,pte=0,ptp=0,此時僅存在機械載荷,式(18)沒有意義。除此之外,式(18)都適用。作為結(jié)合部熱彈塑性接觸問題法向載荷分形模型的另一種表達形式,式(18)的意義和好處是提供了另一種計算和分析的方法,相比于不包含機械載荷比例系數(shù)的模型,它更簡潔,計算更方便。然而,式(18)僅適用于考慮熱應力的情況,對于純機械的情況并不適用,在應用范圍上,與常見的純機械載荷情況下接觸問題的經(jīng)典分形模型(如TZQ模型)互不包含,且沒有交集。從建模過程可以看出,不論是否存在熱載荷,描述結(jié)合部整體熱彈塑性接觸法向總載荷的式(16)、(17)都適用。

2 特例驗證

以下通過特例來驗證本文所建模型的有效性。

當α=0時,熱載荷為0,代入式(16)、(17)可得到僅存在機械載荷時結(jié)合部的整體法向總載荷。

當D≠1.5時

(19)

當D=1.5時

(20)

當ΔT=0時,也不存在熱載荷,將其代入式(16)、(17),可以得到同樣的結(jié)果。

式(19)、(20)即文獻[9]中不考慮熱應力的結(jié)合部整體法向總載荷的經(jīng)典模型。

3 計算結(jié)果對比與分析

下面將本文模型的計算結(jié)果與不考慮熱應力的TZQ模型的計算結(jié)果進行對比和分析。

3.1 兩鑄鐵接觸

在工程機械的靜摩擦接觸范疇內(nèi),廣泛存在鑄鐵-鑄鐵接觸副,本節(jié)以此為例進行對比計算。

表1 兩鑄鐵接觸的整體法向總載荷(TZQ模型)

由于材料的彈性模量和線膨脹系數(shù)會隨溫度升高而改變,故參照鑄鐵的變物性參數(shù)取值,選用在不同工作溫度下的彈性模量與線膨脹系數(shù),其余參數(shù)取值同上。以實際工作中兩鑄鐵材料結(jié)合部最常見的0—500 K溫差為變量范圍,得到整體熱彈塑性接觸法向總載荷隨ΔT變化的關(guān)系曲線,如圖3所示,圖中各條線均為曲率很小的曲線。各條線的起點,即ΔT=0時的p值,就是表1中TZQ模型的計算結(jié)果。為增強對比效果,以G=10-11m時TZQ模型的計算結(jié)果為代表,同時考慮材料彈性模量的變物性特點,通過計算,在圖3中添加了反映該結(jié)果的曲線。

圖3 兩鑄鐵接觸的計算結(jié)果對比

3.2 鑄鐵-不銹鋼接觸

在內(nèi)燃機中,缸蓋與機體常采用鑄鐵材料,與之接觸的缸墊材料常為不銹鋼,故本節(jié)以鑄鐵與不銹鋼接觸為例進行對比計算。

當ΔT=0時,取E′=6.5×1010Pa,α=1.13×10-5K-1,K=0.4,σy=2.5×108Pa。對應不同的G值,根據(jù)TZQ模型計算得到的鑄鐵-不銹鋼接觸整體法向總載荷見表2。

表2 鑄鐵-不銹鋼接觸的法向總載荷(TZQ模型)

選用不同工作溫度下鑄鐵和不銹鋼的彈性模量與線膨脹系數(shù),其余參數(shù)同上,采用本文模型和TZQ模型分別計算,得到了鑄鐵-不銹鋼接觸的法向載荷,見圖4。圖中各條線的起點,即是表2中TZQ模型的計算結(jié)果。

圖4 鑄鐵-不銹鋼接觸的計算結(jié)果對比

通過以上兩個實例可見,即便考慮變物性參數(shù)的影響,TZQ模型的計算結(jié)果受ΔT變化的影響仍很小,且隨著溫差增大,與本文模型間的計算結(jié)果差距增大。這是因為TZQ模型沒有考慮結(jié)合部溫度變化產(chǎn)生的熱載荷對計算結(jié)果帶來的影響,在溫差顯著增大的情況下,其計算結(jié)果不符合客觀預期。

4 溫差對法向總載荷的影響

采用考慮機械載荷比例系數(shù)的熱彈塑性接觸法向載荷分形模型,假設(shè)機械載荷比例系數(shù)不隨ΔT改變,研究溫差對熱彈塑性接觸法向總載荷p的影響規(guī)律。

圖5 溫差ΔT對法向總載荷p的影響曲線

由數(shù)字仿真計算結(jié)果可見,熱彈塑性接觸法向總載荷隨著結(jié)合部溫差的增大而增大,且比例系數(shù)f1、f2越大,法向總載荷增大的幅度也越大。在圖5所示的4組數(shù)據(jù)中,在其他參數(shù)相同的情況下,當f1=1、f2=0.2時,得到的法向總載荷最小。這是因為結(jié)合部溫差增大導致微凸體的變形也增大,繼而使熱彈塑性接觸法向總載荷增大。

5 結(jié) 論

(1)本文建立的不包含機械載荷比例系數(shù)的結(jié)合部熱彈塑性接觸問題法向載荷分形模型,既適用于考慮熱應力的情況,也適用于純機械情況,是傳統(tǒng)純機械模型在基礎(chǔ)理論和應用范圍上的擴展。

(2)本文建立的考慮機械載荷比例系數(shù)的熱彈塑性接觸問題法向載荷分形模型,僅適用于考慮熱應力的情況,作為模型的另一種表達形式,提供了另一種計算和分析的方法,其形式更簡潔,計算更方便。

(3)熱彈塑性接觸法向載荷隨著結(jié)合部溫差的增大而增大。

(4)本文基于熱彈塑性接觸的結(jié)合部法向載荷分形模型綜合考慮了溫差、接觸材料線膨脹系數(shù)、機械載荷比例系數(shù)等參數(shù)的影響,從而可用于計算和分析工程實際中大量存在的結(jié)合部溫度發(fā)生改變的接觸情況,其結(jié)果更符合客觀預期。

參考文獻:

[1] 尤晉閩, 陳天寧. 結(jié)合面法向動態(tài)參數(shù)的分形模型 [J]. 西安交通大學學報, 2009, 43(9): 91-94. YOU Jinmin, CHEN Tianning. Fractal model for normal dynamic parameters of joint surfaces [J]. Journal of Xi’an Jiaotong University, 2009, 43(9): 91-94.

[2] GREENWOOD J A, WILLIMSON J B P. Contact of nominally flat surfaces [J]. Proc R Soc Lond: A, 1996, 295(1442): 300-319.

[3] 姬翠翠, 朱華. 粗糙表面分形接觸模型的研究進展 [J]. 潤滑與密封, 2011, 36(9): 114-119. JI Cuicui, ZHU Hua. Research progress on M-B fractal contact model [J]. Lubrication Engineering, 2011, 36(9): 114-119.

[4] GOERKE D, WILLNER K. Normal contact of fractal surfaces [J]. Wear, 2008, 264(7): 589-598.

[5] 溫淑花, 張學良, 武美先, 等. 結(jié)合面法向接觸剛度分形模型建立與仿真 [J]. 農(nóng)業(yè)機械學報, 2009, 40(11): 197-202. WEN Shuhua, ZHANG Xueliang, WU Meixian, et al. Fractal model and simulation of normal contact stiffness of joint interfaces and its simulation [J]. Transactions of the Chinese Society of Agricultural Machinery, 2009, 40(11): 197-202.

[6] JIANG Shuyun, ZHENG Yunjian, ZHU Hua. A contact stiffness model of machined plane joint based on fractal theory [J]. ASME Journal of Tribology, 2010, 132: 011401.

[7] 蘭國生, 張學良, 丁紅欽, 等. 基于分形理論的結(jié)合面改進接觸模型 [J]. 農(nóng)業(yè)機械學報, 2011, 42(10): 217-223. LAN Guosheng, ZHANG Xueliang, DING Hongqin, et al. Modified contact model of joint interfaces based on fractal theory [J]. Transactions of the Chinese Society of Agricultural Machinery, 2011, 42(10): 217-223.

[8] 宿月文, 陳渭, 朱愛斌, 等. 分形特征表面接觸磨損模擬 [J]. 西安交通大學學報, 2013, 47(7): 52-56. SU Yuewen, CHEN Wei, ZHU Aibin, et al. Contact and wear simulation between fractal surfaces [J]. Journal of Xi’an Jiaotong University, 2013, 47(7): 52-56.

[9] 田紅亮, 鐘先友, 秦紅玲, 等. 依據(jù)各向異性分形幾何理論的固定結(jié)合部法向接觸力學模型 [J]. 機械工程學報, 2013, 49(21): 108-122. TIAN Hongliang, ZHONG Xianyou, QIN Hongling, et al. Normal contact mechanics model of fixed joint interface adopting anisotropic fractal geometrical theory [J]. Chinese Journal of Mechanical Engineering, 2013, 49(21): 108-122.

[10]田紅亮, 陳從平, 方子帆, 等. 應用改進分形幾何理論的結(jié)合部切向剛度模型 [J]. 西安交通大學學報, 2014, 48(7): 46-52. TIAN Hongliang, CHEN Congping, FANG Zifan, et al. Tangential stiffness model for joint interface adopting the revised fractal geometric theory [J]. Journal of Xi’an Jiaotong University, 2014, 48(7): 46-52.

[11]劉天祥, 劉更, 謝琴, 等. 粗糙表面熱彈塑性接觸問題的無網(wǎng)格法 [J]. 機械工程學報, 2009, 45(1): 20-26. LIU Tianxiang, LIU Geng, XIE Qin, et al. Mesh-free method for thermal elasto-plastic contact of rough surfaces [J]. Chinese Journal of Mechanical Engineering, 2009, 45(1): 20-26.

[12]肖永寧, 潘克煜, 韓國埏. 內(nèi)燃機熱負荷和熱強度 [M]. 北京: 機械工業(yè)出版社, 1988: 162.

[13]成大先. 機械設(shè)計手冊: 第1卷 [M]. 北京: 化學工業(yè)出版社, 2008: 38.

[本刊相關(guān)文獻鏈接]

褚坤明,王家序,蒲偉,等.混合潤滑狀態(tài)下濾波減速器的嚙合沖擊分析與修形方法.2015,49(7):140-148.[doi:10.7652/xjtuxb201507023]

周安安,陳天寧,王小鵬,等.自仿射接觸點及其在分形接觸理論中的應用.2015,49(6):8-15.[doi:10.7652/xjtuxb2015 06002]

陳衛(wèi)華,陳天寧,王小鵬,等.纖維多孔金屬的流阻率分形模型研究.2015,49(6):132-137.[doi:10.7652/xjtuxb201506 021]

陳兵奎,梁棟,彭帥,等.共軛曲線齒輪嚙合管齒面的幾何及接觸特性分析.2015,49(3):85-94.[doi:10.7652/xjtuxb 201503015]

田紅亮,趙美云,鄭金華,等.新的柔性結(jié)合部法向接觸剛度和接觸阻尼方程.2015,49(1):118-126.[doi:10.7652/xjtuxb201501020]

田紅亮,陳從平,方子帆,等.應用改進分形幾何理論的結(jié)合部切向剛度模型.2014,48(7):46-52.[doi:10.7652/xjtuxb 201407009]

宿月文,陳渭,朱愛斌,等.分形特征表面接觸磨損模擬.2013,47(7):52-56.[doi:10.7652/xjtuxb201307010]

田紅亮,趙春華,朱大林,等.整個螺栓結(jié)合部的法向連接動剛度及試驗驗證.2012,46(9):31-36.[doi:10.7652/xjtuxb201209007]

方兵,張雷,趙繼,等.軸承結(jié)合部動態(tài)參數(shù)識別與等效分析模型的研究.2012,46(11):69-74.[doi:10.7652/xjtuxb2012 11014]

張學良,溫淑花,蘭國生,等.平面結(jié)合面切向接觸阻尼分形模型及其仿真.2011,45(5):74-77.[doi:10.7652/xjtuxb 201105014]

尤晉閩,陳天寧.結(jié)合面法向動態(tài)參數(shù)的分形模型.2009,43(9):91-94.[doi:10.7652/xjtuxb200909020]

(編輯 葛趙青)

(1.浙江大學動力機械及車輛工程研究所,310027,杭州;2.浙江同濟科技職業(yè)學院機電系,311231,杭州)

針對現(xiàn)有法向接觸載荷分形模型沒有考慮熱應力的影響、不適用于結(jié)合部溫度變化的問題,基于熱彈塑性理論和接觸分形理論,建立了熱彈塑性接觸法向載荷分形模型。該模型是對傳統(tǒng)模型在基礎(chǔ)理論和應用范圍上的拓展,綜合考慮了溫差、線膨脹系數(shù)、機械載荷比例系數(shù)等參數(shù)的影響,從而可用于計算和分析工程實際中大量存在的結(jié)合部溫度發(fā)生改變的接觸情況,其結(jié)果更符合客觀規(guī)律。此外,針對考慮熱應力的接觸問題,引入機械載荷比例系數(shù),給出了模型的另一種表達形式。通過特例驗證了模型的有效性,并與前人的模型進行了對比分析,結(jié)果顯示:模型間的計算結(jié)果差距隨結(jié)合部溫差的增大而增大。最后,通過數(shù)字仿真揭示了結(jié)合部溫差對熱彈塑性接觸法向載荷的影響規(guī)律——熱彈塑性接觸法向載荷隨結(jié)合部溫差的增大而增大。

熱彈塑性;法向接觸載荷;結(jié)合部;比例系數(shù);分形模型

A Normal Load Fractal Model of Joint Interface Based on Thermal Elasto-Plastic Theory

FENG Yan1,2,YU Xiaoli1,LIU Zhentao1

(1. Power Machinery & Vehicular Engineering Institute, Zhejiang University, Hangzhou 310027, China; 2. Department of Mechanical and Electrical Engineering, Zhejiang Tongji Vocational College of Science and Technology, Hangzhou 311231, China)

Without considering the influence of heat, the existing fractal contact models are not applicable to analyze the contact issues when the temperature changes. Based on thermal elasto-plastic theory and contact fractal theory, a new fractal model of normal contact load is established in this study. This model expands the basic theory and applications of traditional models, and comprehensively considers the effects of such parameters as temperature difference, linear expansion coefficient, mechanical load scale factor, etc. Therefore this model can help calculate and analyze the real contact conditions between fixed contact surfaces when the temperature changes, and it shows better conformity with objective laws. Besides, another expression of the model is proposed by introducing mechanical load scale factor. Then the model is validated through a special case, and comparative study is conducted between this model and TZQ model. The results show that the difference between the model calculations increases with the junction temperature difference. In addition, the effect of temperature difference on the normal contact load is also analyzed. The results reveal that the normal contact load of thermal elasto-plastic contact increases with the temperature difference at joint interface.

thermal elasto-plastic; normal contact load; joint interface; scale factor; fractal model

2015-01-15。 作者簡介:馮燕(1985—),女,博士生,浙江同濟科技職業(yè)學院講師;劉震濤(通信作者),男,副教授。 基金項目:國家高技術(shù)研究發(fā)展計劃資助項目(2012AA111709)。

時間:2015-06-19

http:∥www.cnki.net/kcms/detail/61.1069.T.20150619.1649.002.html

10.7652/xjtuxb201509004

TH123

A

0253-987X(2015)09-0018-06