基于W-M分形函數(shù)的三維粗糙表面摩擦生熱研究

,,,,,

(大連交通大學(xué) 機(jī)械工程學(xué)院,遼寧 大連 116028)

大功率風(fēng)電制動(dòng)器具有制動(dòng)轉(zhuǎn)速高、制動(dòng)力矩大的特點(diǎn),其制動(dòng)過(guò)程中大部分動(dòng)能通過(guò)摩擦作用轉(zhuǎn)化為熱能,制動(dòng)器摩擦副表面將產(chǎn)生大量摩擦熱.制動(dòng)閘片由于局部高溫和應(yīng)力集中的原因，材料屬性發(fā)生改變,造成閘片不均勻摩擦損耗的加劇,影響了風(fēng)電制動(dòng)器的制動(dòng)性能,并降低了制動(dòng)閘片材料的利用率,因此,如能對(duì)制動(dòng)過(guò)程中摩擦接觸的微觀過(guò)程進(jìn)行深入研究,即可有效預(yù)測(cè)摩擦副接觸表面的力-熱分布狀態(tài),進(jìn)而為改善制動(dòng)摩擦副的工作性能,提高制動(dòng)閘片的使用壽命提供理論依據(jù).

制動(dòng)過(guò)程中制動(dòng)盤與閘片接觸表面的摩擦接觸,其微觀實(shí)質(zhì)是兩個(gè)粗糙表面上一系列不規(guī)則微凸體的相互接觸過(guò)程.實(shí)際滑動(dòng)摩擦過(guò)程中,由于真實(shí)接觸面積遠(yuǎn)小于名義接觸面積且不連續(xù),造成摩擦副間很小的真實(shí)接觸面積上承擔(dān)很大的實(shí)際載荷,接觸微凸體將發(fā)生彈塑性變形,并形成“熱點(diǎn)”[1].這些“熱點(diǎn)”就是制動(dòng)摩擦副的局部閃溫.局部溫度過(guò)高,及由此導(dǎo)致的摩擦副材料屬性變化會(huì)引起兩接觸表面狀態(tài)的變化.粗糙表面之間的滑動(dòng)摩擦實(shí)質(zhì)上是一個(gè)復(fù)雜的非線性接觸問(wèn)題.粗糙表面的形貌和接觸特性對(duì)制動(dòng)過(guò)程摩擦、磨損和傳熱都有著重要的影響.因此,研究粗糙表面輪廓形貌以及其對(duì)滑動(dòng)摩擦中的溫度與應(yīng)力的影響有重要的應(yīng)用價(jià)值.

一直以來(lái),國(guó)內(nèi)外學(xué)者就粗糙表面輪廓形貌進(jìn)行了表征，并對(duì)其接觸摩擦的過(guò)程進(jìn)行了大量的研究.Archard[2]在其研究中首先體現(xiàn)出了分形的思想,將多尺度粗糙表面看作較大尺度的球形微凸體上承載了一簇較小尺度的球形微凸體.Mandelbort[3]首先提出了分形這一概念,并創(chuàng)立了分形幾何學(xué),用于解釋那些不規(guī)則的、破碎的、參差不齊的和斷裂的形狀;其次在Weierstrass函數(shù)的基礎(chǔ)上,提出了一種分形曲線函數(shù)的表達(dá)式,稱為Weierstrass-Mandelbrot分形函數(shù)(簡(jiǎn)稱W-M分形函數(shù)).Mandelbort提出的分形理論被逐漸應(yīng)用到摩擦學(xué)的研究中.Bhushan[4]基于W-M分形函數(shù)提出了M-B分形模型,用W-M分形函數(shù)來(lái)模擬粗糙表面的輪廓線,把模擬兩粗糙表面的接觸簡(jiǎn)化為一等效粗糙表面與一理想剛性光滑平面的接觸,提出了分形接觸模型.葛世榮等[5]研究了通過(guò)磨削、車削等加工方法得到的表面輪廓曲線,發(fā)現(xiàn)粗糙表面具有明顯的分形特征,提出分形維數(shù)與表面粗糙度之間呈負(fù)指數(shù)關(guān)系,并定義了特征粗糙度這一概念.魏龍等[6]考慮微凸體的變形特征和摩擦作用的影響，建立了滑動(dòng)摩擦表面的分形接觸模型,采用一個(gè)三次多項(xiàng)式來(lái)表達(dá)彈塑性變形微凸體的接觸壓力與接觸面積的關(guān)系,并推導(dǎo)出不同臨界條件下微凸體的真實(shí)接觸面積.鄧可月等[7]利用分形理論對(duì)表面形貌的分形特點(diǎn)進(jìn)行研究,在建立W-M分形函數(shù)模型的前提下，對(duì)表面輪廓形貌進(jìn)行二維及三維的模擬仿真分析,得出分形維數(shù)是影響表面輪廓曲面形貌的主要參數(shù).韓傳軍等[8]建立了一個(gè)含球形微凸體粗糙表面與理想平面的滑動(dòng)接觸模型,探究了特定形狀微凸體在摩擦過(guò)程中的應(yīng)力與溫度變化規(guī)律.本文基于W-M分形理論建立粗糙面微接觸的滑動(dòng)計(jì)算模型,將兩個(gè)粗糙表面簡(jiǎn)化為一個(gè)分形粗糙表面與一個(gè)光滑平面的組合,建立粗糙面的滑動(dòng)接觸有限元計(jì)算模型,結(jié)合大功率風(fēng)電制動(dòng)器高速重載的制動(dòng)工況,在不同分形維數(shù)、相對(duì)運(yùn)動(dòng)速度及施加載荷的邊界條件下,模擬并分析了制動(dòng)過(guò)程中粗糙表面的摩擦生熱及熱應(yīng)力變化規(guī)律.分析結(jié)果為進(jìn)一步研究制動(dòng)過(guò)程中微觀摩擦機(jī)理研究、粗糙表面摩擦副的閃點(diǎn)溫度、摩擦副接觸表面的力-熱分布狀態(tài)提供了參考依據(jù).

1 基于W-M分形函數(shù)的粗糙表面模型建立

W-M分形函數(shù)可以準(zhǔn)確地模擬和重構(gòu)具有分形特性的分形表面,且處處連續(xù)但不可微,其本身具有自相似性,是用于表示隨機(jī)輪廓的一種典型函數(shù)[9].其適用于工程表面的數(shù)學(xué)模型表達(dá)式為

(1)

式中:Z(x)為隨機(jī)表面形貌的高度;x為輪廓的位置坐標(biāo);D為分形維數(shù),它描述函數(shù)Z(x)在所有尺度上的不規(guī)則性;G為特征尺度系數(shù),它反映Z(x)幅值大小,并決定Z(x)的具體尺寸；γ為輪廓的空間頻率，對(duì)于服從正態(tài)分布的隨機(jī)輪廓,為適用于高頻譜密度及相位的隨機(jī)性,一般γ的取值為1.5.基于W-M分形函數(shù)的三維分形表面的函數(shù)模型為

ysinBn)+An] (2

(2)

式中:Cn為尺度系數(shù)，是服從均值為0、方差為1的正態(tài)分布的隨機(jī)數(shù);An與Bn為相互獨(dú)立且服從[0,2π]均勻分布的隨機(jī)數(shù);Ds為理論分形維數(shù);n為自然序列數(shù).

以W-M分形函數(shù)公式為理論基礎(chǔ),借助Matlab軟件進(jìn)行編程,在固定某些特定參數(shù)的條件下,對(duì)分形維數(shù)進(jìn)行變換,模擬出不同分形維數(shù)下的粗糙表面形貌,如圖1所示.固定參數(shù)中尺度系數(shù)取定值C=0.01,自然序列數(shù)為n=1,2,…,100.在選取橫縱坐標(biāo)x,y時(shí),由于分形函數(shù)公式的特性,要求x≠0,y≠0.因此，橫縱坐標(biāo)的選取區(qū)間定在區(qū)間[1, 2]內(nèi).在區(qū)間內(nèi),x方向與y方向上每間距0.04上取一個(gè)點(diǎn),取樣點(diǎn)的數(shù)量為x方向和y方向各26個(gè),總采樣點(diǎn)數(shù)為676個(gè)；再根據(jù)An和Bn的兩個(gè)隨機(jī)數(shù),求出每個(gè)采樣點(diǎn)的z坐標(biāo)值,進(jìn)行樣條擬合并最終繪制出三維狀態(tài)下不同分形維數(shù)的微觀表面輪廓形貌.

圖1 不同分形維數(shù)三維粗糙表面輪廓形貌Fig.1 Three-dimensional rough surface profile topography at different fractal dimension

圖1(a)～圖1(e)分別為分形維數(shù)2.1,2.3,2.5,2.7和2.9時(shí)獲得的三維粗糙表面輪廓形貌.從圖1可以看出,粗糙表面輪廓曲面具有不規(guī)則性,呈現(xiàn)出高低不平的形狀,不同分形維數(shù)下的曲面凹凸程度差別較明顯,不僅顯露出了輪廓形狀的隨機(jī)性,還體現(xiàn)出了表面形貌的復(fù)雜性.隨著分形維數(shù)值的逐漸增大,粗糙表面輪廓形貌趨于復(fù)雜化,表面凹凸程度趨于密集化,且高度幅值變化頻率也越來(lái)越快,這說(shuō)明了分形維數(shù)是影響表面輪廓曲面形貌的主要參數(shù),這為通過(guò)改變分形維數(shù)建立不同粗糙程度的摩擦表面模型提供了理論依據(jù).

2 摩擦生熱有限元模型的建立

2.1 三維粗糙表面模型的建立

由于粗糙表面輪廓形貌的不規(guī)則性,若將此三維模型直接用于有限元模擬,不僅需要增加網(wǎng)格數(shù)量導(dǎo)致運(yùn)算量顯著上升,而且很容易造成接觸分析不收斂的情況,因此，需要對(duì)粗糙表面輪廓進(jìn)行光滑處理.運(yùn)用Matlab中的樣條擬合功能,將分形粗糙表面輪廓擬合成光滑表面輪廓形貌,擬合后表面形貌的凹凸程度趨勢(shì)并沒有發(fā)生太大改變,但表面光滑程度較高,適合三維模型的建立和有限元仿真.選取其中一組分形維數(shù)Ds=2.7的三維粗糙表面輪廓,將其擬合后的表面各點(diǎn)數(shù)據(jù)導(dǎo)入到Creo中生成三維模型,粗糙體表面的輪廓形貌如圖2所示.

圖2 擬合后分形維數(shù)Ds=2.7的粗糙表面三維模型Fig.2 Three-dimensional model of rough surface with Ds=2.7 fractal dimension after fitting

2.2 邊界條件與有限元計(jì)算模型

2.2.1熱邊界條件

在有限元接觸分析的設(shè)定中,接觸面單元不能穿透目標(biāo)面,而目標(biāo)面單元可以穿透接觸面,且目標(biāo)面一般選擇標(biāo)準(zhǔn)為平面、網(wǎng)格較粗和面積較大的面,而接觸面一般都選擇硬度偏軟的面.根據(jù)這一性質(zhì),選擇制動(dòng)盤為光滑表面作為目標(biāo)面,而制動(dòng)閘片為粗糙表面作為接觸面.在制動(dòng)過(guò)程中,摩擦生熱產(chǎn)生的能量為

(3)

式中:FH是摩擦耗散能轉(zhuǎn)換成能量的比例系數(shù),當(dāng)認(rèn)為摩擦全部轉(zhuǎn)化為熱量的情況下，F(xiàn)H的值為1;τ為等效摩擦應(yīng)力;v為滑動(dòng)率.分配在接觸面和目標(biāo)面的摩擦生熱能量分別為qc和qt,即

(4)

式中:FW為接觸面與目標(biāo)面之間的熱量分配權(quán)重因子,qc和qt的關(guān)系為

(5)

式中:ρc,ρt分別為接觸面和目標(biāo)面的密度;Cc,Ct分別為接觸面和目標(biāo)面的比熱容;kc,kt分別為接觸面和目標(biāo)面的熱導(dǎo)率[10].在代入制動(dòng)盤(目標(biāo)面)和制動(dòng)閘片(接觸面)材料屬性計(jì)算得熱分配系數(shù)分別為31%和69%.

在制動(dòng)過(guò)程接觸的熱邊界條件中，熱傳導(dǎo)占主導(dǎo)地位,摩擦接觸區(qū)域溫度Tc=Tt,其中，接觸面和目標(biāo)面中，接觸部分的熱傳導(dǎo)方程為

(6)

接觸面與目標(biāo)面未接觸部分主要的熱邊界條件為熱對(duì)流,其中，接觸面和目標(biāo)面的非接觸部分的熱對(duì)流方程為

(7)

式中:Tα為流體介質(zhì)的溫度,即剛接觸時(shí)空氣的溫度;h12為接觸面與流體介質(zhì)的對(duì)流換熱系數(shù);h13為目標(biāo)面與流體介質(zhì)的對(duì)流換熱系數(shù).經(jīng)計(jì)算,將對(duì)流換熱系數(shù)設(shè)定為13.35 N·m/(h·km2).

2.2.2有限元計(jì)算模型

為簡(jiǎn)化模型,提高分析效率,對(duì)摩擦生熱模型作如下基本假設(shè):① 兩摩擦件的材料組成均勻,且各向同性;② 制動(dòng)閘片材料(銅基粉末冶金)和制動(dòng)盤材料(Q345-B)的密度、比熱、熱導(dǎo)率、熱膨脹系數(shù)、彈性模量、泊松比等各項(xiàng)物理參數(shù)均為常數(shù),且不受溫度變化影響;③ 由于滑動(dòng)過(guò)程時(shí)間極短,因此，僅考慮熱傳導(dǎo)和熱對(duì)流的影響,忽略熱輻射的作用;④ 忽略粗糙表面的磨損及其影響,認(rèn)為動(dòng)能完全轉(zhuǎn)化為摩擦熱而被摩擦副吸收;⑤ 摩擦過(guò)程中,符合庫(kù)侖摩擦定律,摩擦系數(shù)保持不變.兩摩擦件材料的熱物理參數(shù)如表1所示.

表1 摩擦件材料物理參數(shù)Tab.1 Friction material physical parameters

為簡(jiǎn)化運(yùn)算,制動(dòng)盤和制動(dòng)閘片均采用局部模型,對(duì)兩者的三維模型進(jìn)行網(wǎng)格劃分,劃分后的有限元模型如圖3所示.圖3中:上方模型為分形粗糙表面采用四面體單元?jiǎng)澐志W(wǎng)格,共計(jì)有19 392個(gè)四面體單元(C3D4T);下方模型為光滑表面,理想光滑平面體外形規(guī)則,為減少運(yùn)算量,選用六面體單元(C3D8T)劃分網(wǎng)格,共計(jì)有9 600個(gè)六面體單元.由于接觸計(jì)算及彈塑性計(jì)算的復(fù)雜性,考慮到計(jì)算成本,滑動(dòng)距離為L(zhǎng)=20 mm.運(yùn)動(dòng)過(guò)程分為2個(gè)階段:第1階段在分形表面上施加載荷,使分形粗糙表面往下運(yùn)動(dòng),擠壓光滑表面,在光滑表面的底面限制Z方向自由度,并對(duì)粗糙體的側(cè)表面分別施加法向約束;第2階段對(duì)光滑平面分別施加不同的速度載荷.

圖3 制動(dòng)過(guò)程分形粗糙體與光滑表面有限元模型Fig.3 Finite element model of fractal aperture and smooth surface in braking process

3 摩擦生熱仿真結(jié)果分析

制動(dòng)過(guò)程中,各種因素共同作用影響摩擦生熱.為探究不同因素各自對(duì)于摩擦生熱的影響程度,在保證其他參數(shù)不變的情況下,分別改變單一變量以研究其對(duì)制動(dòng)過(guò)程中摩擦生熱的影響.因此，選取不同分形維數(shù)、滑動(dòng)速度以及施加載荷的條件，對(duì)制動(dòng)過(guò)程分形粗糙表面的摩擦生熱進(jìn)行探究,結(jié)合大功率風(fēng)電制動(dòng)器制動(dòng)的實(shí)際工況,選取分形維數(shù)Ds分別為2.1,2.3,2.5和2.7,運(yùn)動(dòng)速度v分別為10,20,30和40 m/s,施加載荷P分別為10,20,30和40 MPa,以探究不同單一變量變化對(duì)分形粗糙表面摩擦生熱的影響.

3.1 分形維數(shù)對(duì)摩擦生熱的影響分析

圖4(a)～圖4(d)分別是施加載荷為30 MPa時(shí),分形維數(shù)Ds為2.1,2.3,2.5和2.7不同情況下，摩擦初期的表面熱區(qū)分布形貌圖.圖4(a)為分形維數(shù)Ds=2.1時(shí)的溫度分布云圖,從圖中可以看出其熱區(qū)分布形式主要以塊狀熱區(qū)為主,點(diǎn)狀熱區(qū)的數(shù)量?jī)H占極少數(shù),而在制動(dòng)初期其最大閃點(diǎn)溫度為146.7 ℃.圖4(b)為分形維數(shù)Ds=2.3時(shí)的溫度分布云圖,相較于圖4(a)其塊狀熱區(qū)的數(shù)量有所減少,而點(diǎn)狀熱區(qū)的數(shù)量呈增長(zhǎng)的趨勢(shì),其制動(dòng)初期的最大閃點(diǎn)溫度為200.7 ℃.圖4(c)為分形維數(shù)Ds=2.5時(shí)的溫度分布云圖,塊狀熱區(qū)的接觸面積與數(shù)量都有大幅度減少,而點(diǎn)狀熱區(qū)的數(shù)量也超過(guò)了塊狀熱區(qū)的數(shù)量.其制動(dòng)初期的最大閃點(diǎn)溫度為276.2 ℃.圖4(d)為分形維數(shù)Ds=2.7時(shí)的溫度分布云圖,其熱區(qū)基本上都為點(diǎn)狀熱區(qū),而制動(dòng)初期的最大閃點(diǎn)溫度為317.0 ℃.從圖4的溫度分布云圖可以分析出,當(dāng)分形維數(shù)較小時(shí),摩擦初期形成的熱區(qū)以塊狀熱區(qū)居多,而點(diǎn)狀熱區(qū)的數(shù)量偏少,隨著分形維數(shù)的增大,點(diǎn)狀熱區(qū)的數(shù)量增多,而塊狀熱區(qū)的數(shù)量減少,形成這種趨勢(shì)的原因是當(dāng)分形維數(shù)較小時(shí),表面形貌復(fù)雜程度低,分形粗糙表面單位面積內(nèi)的波峰波谷數(shù)量也較少,表面更為平坦,因此，更容易形成塊狀熱區(qū).而隨著分形維數(shù)增大,表面形貌的復(fù)雜程度也隨之上升,表面凹凸程度趨于密集化,使得單位面積內(nèi)的波峰波谷數(shù)量也隨之增加,點(diǎn)狀熱區(qū)形成的概率隨之增大.此外,隨著分形維數(shù)增大,粗糙體接觸界面的最高溫度也呈現(xiàn)上升趨勢(shì).這是因?yàn)榉中尉S數(shù)較小時(shí),接觸部分大多是面接觸,接觸部分最大接觸壓力較小,因此，摩擦初期摩擦生熱的最大溫度值較小;當(dāng)分形維數(shù)較大時(shí),面接觸部分減少,點(diǎn)接觸部分增加,導(dǎo)致最大接觸壓力值也隨之增大,其接觸的微凸體的最大閃點(diǎn)溫度也隨之升高.由圖4還可看出,雖然分形粗糙曲面在構(gòu)造時(shí)就具有隨機(jī)性的特點(diǎn),導(dǎo)致摩擦生熱形成的熱區(qū)分布趨勢(shì)也存在著一定的隨機(jī)性,但分形維數(shù)數(shù)值的變化還是會(huì)對(duì)摩擦生熱產(chǎn)生較大影響.

圖4 不同分形維數(shù)情況下粗糙表面熱區(qū)分布圖Fig.4 Rough surface hot zone distribution at different fractal dimension

3.2 運(yùn)動(dòng)速度對(duì)摩擦生熱的影響分析

圖5為不同相對(duì)運(yùn)動(dòng)速度下,滑動(dòng)時(shí)間為500 μs內(nèi)的分形粗糙表面最頂端微凸體節(jié)點(diǎn)(節(jié)點(diǎn)814)的溫度隨時(shí)間變化歷程圖.由圖5可以發(fā)現(xiàn),在制動(dòng)過(guò)程中,溫度隨著時(shí)間的變化呈上升趨勢(shì).根據(jù)溫度曲線變化的趨勢(shì),可以將其分為2個(gè)階段:第1個(gè)階段(0～20 μs)為急劇增長(zhǎng)階段；第2階段(20～500 μs)為緩慢增長(zhǎng)階段.從圖5中可以看出,相對(duì)運(yùn)動(dòng)速度越高,急劇增長(zhǎng)階段微凸體的溫升就越大,緩慢增長(zhǎng)階段微凸體的最高溫度值也越大.當(dāng)速度分別為10,20,30和40 m/s時(shí),粗糙峰最大溫度分別為215.295 ℃,388.176 ℃,539.303 ℃和691.053 ℃.

圖5 分形粗糙體最高微凸體(節(jié)點(diǎn)814)溫度-時(shí)間曲線圖Fig.5 Temperature-time curve of highest asperity (Node 814) in fractal aperture

不同相對(duì)運(yùn)動(dòng)速度下,粗糙體非接觸面區(qū)域節(jié)點(diǎn)(節(jié)點(diǎn)504)溫度-時(shí)間變化曲線如圖6所示.由圖6可見,非接觸節(jié)點(diǎn)的溫度隨著制動(dòng)時(shí)間的推進(jìn)，呈近似線性增長(zhǎng)的趨勢(shì).當(dāng)相對(duì)運(yùn)動(dòng)速度越大時(shí),曲線的斜率也越大,這是因?yàn)樗俣仍娇煳⑼贵w的最高溫度越大,微凸體與非接觸部分溫差越大,導(dǎo)致非接觸面區(qū)域的熱傳遞速率越快,其結(jié)論也完全符合傅里葉定律.

3.3 施加載荷對(duì)摩擦生熱的影響分析

圖7(a)～圖7(d)分別是運(yùn)動(dòng)速度為20 m/s,分形維數(shù)Ds=2.7條件下,施加載荷為10,20,30和40 MPa的分形粗糙表面接觸壓力云圖.從圖7(a)中可以看出,當(dāng)施加載荷為10 MPa時(shí),加載后的制動(dòng)初期最大接觸壓力為1 588 MPa,實(shí)際接觸點(diǎn)的數(shù)量為7個(gè).在圖7(b)中,施加載荷為20 MPa,加載后的制動(dòng)初期最大接觸壓力為1 714 MPa,實(shí)際接觸點(diǎn)較圖7(a)有所增加,其數(shù)量大約為14個(gè).在圖7(c)中,加載后的制動(dòng)初期最大接觸壓力為1 653 MPa,實(shí)際接觸點(diǎn)較圖7(b)有小幅增加,其數(shù)量大約為18個(gè)左右.圖7(d)為施加載荷為40 MPa的接觸壓力云圖,其在加載后的制動(dòng)初期最大接觸壓力為1 564 MPa,實(shí)際接觸點(diǎn)數(shù)量在26個(gè)左右.根據(jù)圖7可以看出:當(dāng)施加的載荷越大,真實(shí)的接觸面積越大,接觸點(diǎn)數(shù)量越多,但接觸點(diǎn)數(shù)量呈非線性增長(zhǎng)趨勢(shì);且實(shí)際接觸點(diǎn)的接觸壓力遠(yuǎn)大于施加載荷,最大接觸壓力數(shù)值處于1 600 MPa左右;隨著施加載荷的增加,表面最大接觸壓力出現(xiàn)了一定程度起伏的情況,造成此情況的原因是真實(shí)的粗糙表面形狀具有隨機(jī)性與不規(guī)則性,使得接觸面積并不是隨著載荷增大而等比例的增大,這就導(dǎo)致接觸壓力并不隨載荷增大而增大,而是在表面最大接觸壓力附近小幅波動(dòng).

圖6 分形粗糙體非接觸面(節(jié)點(diǎn)504)溫度-時(shí)間曲線圖Fig.6 Temperature-time curve of non-contact surface (Node 504)in fractal aperture

圖8(a)～圖8(d)分別是施加載荷為10,20,30和40 MPa的分形粗糙表面的溫度云圖.從圖8(a)中可以看出,當(dāng)施加載荷P=10 MPa在制動(dòng)過(guò)程結(jié)束時(shí),其最大閃溫為443.7 ℃,并形成以最高閃溫點(diǎn)為中心的一整塊熱區(qū),且溫度梯度的趨勢(shì)明顯.在圖8(b)中，當(dāng)施加載荷P=20 MPa時(shí),相較于圖8(a)其熱區(qū)的范圍有所擴(kuò)張,熱區(qū)數(shù)量也有一定增幅,其制動(dòng)末期的最大閃點(diǎn)溫度為456.0 ℃.在圖8(c)中,施加載荷后熱區(qū)的接觸面積與數(shù)量都有一定量的增加,其制動(dòng)末期的最大閃點(diǎn)溫度為453 ℃.圖8(d)為施加載荷P=40 MPa時(shí)的溫度分布云圖,可以看出,在制動(dòng)末期熱區(qū)數(shù)量的分布范圍廣且密集,且不少獨(dú)立熱區(qū)在最后合并成了整塊熱區(qū),其最大閃溫為454.1 ℃.從圖8中可以看出,當(dāng)施加載荷越大時(shí),真實(shí)接觸面積隨之增大,在摩擦過(guò)程中,形成的熱區(qū)數(shù)量越多,隨著時(shí)間的推移,熱區(qū)的面積也逐漸擴(kuò)張.但在不同載荷作用下,其微凸體的最高溫度基本都處于450.0 ℃左右,這是由于最大閃溫點(diǎn)一般都是分形粗糙曲面的最高點(diǎn),也是分形粗糙曲面最大接觸壓力點(diǎn).由前文分析可知,不同接觸表面最大接觸壓力值雖有較小波動(dòng),但是其值相差不大,因此，最高閃溫點(diǎn)由于摩擦生熱所產(chǎn)生的熱量也相近,最高閃溫點(diǎn)的溫度也近似相同.但由于載荷增大、熱區(qū)增多的原因,次閃溫點(diǎn)的數(shù)量也會(huì)增多,且受到熱區(qū)增多的影響,非接觸區(qū)域受到更多熱區(qū)熱傳導(dǎo)的作用,因此，載荷越大非接觸區(qū)域單位時(shí)間內(nèi)溫升越顯著.綜上所述,施加載荷的大小對(duì)于微凸體的最高閃溫點(diǎn)的溫度影響并不大,但會(huì)對(duì)熱區(qū)的數(shù)量以及一些次閃溫點(diǎn)和非接觸點(diǎn)的溫度造成一定的影響.

圖7 不同施加載荷下分形粗糙表面接觸壓力分布云圖Fig.7 Fractal rough surface contact pressure distribution at different applied loads

圖8 不同施加載荷下分形粗糙表面溫度分布云圖Fig.8 Fractal rough surface temperature distribution at different applied loads

4 結(jié)論

本文基于W-M分形函數(shù),建立了制動(dòng)接觸摩擦副三維粗糙表面模型,將兩個(gè)粗糙表面簡(jiǎn)化為一個(gè)分形粗糙表面與一個(gè)光滑平面的組合.結(jié)合熱傳遞理論及風(fēng)電制動(dòng)過(guò)程的實(shí)際工況,建立了制動(dòng)盤與閘片粗糙表面的滑動(dòng)摩擦有限元模型,模擬并分析了制動(dòng)過(guò)程中粗糙表面的摩擦生熱及熱應(yīng)力變化規(guī)律,得到如下結(jié)論:

(1) 粗糙表面的分形維數(shù)對(duì)摩擦生熱有較大影響,當(dāng)分形維數(shù)較小時(shí),摩擦初期形成的熱區(qū)以塊狀熱區(qū)數(shù)量較多,而點(diǎn)狀熱區(qū)的數(shù)量偏少.隨著分形維數(shù)的增大,點(diǎn)狀熱區(qū)的數(shù)量增多,而塊狀熱區(qū)的數(shù)量減少,且粗糙體接觸界面的最高溫度也呈現(xiàn)上升的趨勢(shì).

(2) 分形粗糙表面最高溫度隨運(yùn)動(dòng)時(shí)間的變化過(guò)程可分為急劇增長(zhǎng)和緩慢增長(zhǎng)兩個(gè)階段.相對(duì)運(yùn)動(dòng)速度越高,接觸區(qū)域最頂端的微凸體節(jié)點(diǎn)在急劇增長(zhǎng)階段的溫升越大,在緩慢增長(zhǎng)階段的最高溫度值也越大,并且非接觸面區(qū)域的溫升速率也越快.

(3) 施加載荷越大,接觸點(diǎn)數(shù)量越多,真實(shí)的接觸面積越大,且實(shí)際接觸點(diǎn)的接觸壓力遠(yuǎn)大于施加載荷,但接觸壓力不隨著載荷增大而成比例增大.隨著施加載荷的增加,微凸體的最高閃溫點(diǎn)的溫度變化幅度不大,但會(huì)影響熱區(qū)的數(shù)量大小與次閃溫點(diǎn)和非接觸點(diǎn)的溫度.

參考文獻(xiàn):

[1] 黃健萌,高誠(chéng)輝,李友遐,等.工程粗糙表面接觸摩擦熱力學(xué)研究進(jìn)展[J].中國(guó)工程機(jī)械學(xué)報(bào),2007(4):490-496.

HUANG J M,GAO C H,LI Y X,et al.Advances of frictional contact thermodynamics for fractal rough surfaces[J].Chinese Journal of Construction Machinery,2007(4):490-496.

[2] ARCHARD J F.Elastic deformation and the laws of friction[J].Proceedings of the Royal Society A:Mathematical Physical and Engineering Sciences,1957,243(1233):190-205.

[3] MANDELBROT B B.Fractals:form,chance,and dimension[M].San Francisco:W.H.Freeman and Company,1977.

[4] BHUSHAN M B.Fractal model of elastic-plastic contact between rough surfaces[J].ASME Journal of Tribology,1991,113:1-11.

[5] 葛世榮,TONDER K.粗糙表面的分形特征與分形表達(dá)研究[J].摩擦學(xué)學(xué)報(bào),1997,17(1):73-80.

GE S R,TONDER K.The fractal behavior and fractal characterization of rough surfaces[J].Tribology,1997,17(1):73-80

[6] 魏龍,劉其和,張鵬高.基于分形理論的滑動(dòng)摩擦表面接觸力學(xué)模型[J].機(jī)械工程學(xué)報(bào),2012,48(17):106-113.

WEI L,LIU Q H,ZHANG P G.Sliding friction surface contact mechanics model based on fractal theory[J].Chinese Journal of Mechanical Engineering,2012,48(17):106-113.

[7] 鄧可月,劉政,鄧居軍,等.W-M函數(shù)模型下表面輪廓形貌的變化規(guī)律[J].機(jī)械設(shè)計(jì)與制造,2017(1):47-50.

DENG K Y,LIU Z,DENG J J,et al.Variation of surface profile topography based on W-M function model[J].Machinery Design and Manufacture,2017(1):47-50.

[8] 韓傳軍,張杰,梁政.粗糙表面在滑動(dòng)過(guò)程中的摩擦生熱研究[J].西南石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,37(2):159-164.

HAN C J,ZHANG J,LIANG Z.Study on frictional heating of rough surface in the sliding process[J].Journal of Southwest Petroleum University (Science and Technology Edition),2015,37(2):159-164.

[9] 李成貴,董申.三維表面微觀形貌的表征趨勢(shì)[J].中國(guó)機(jī)械工程,2000,11(5):488-492.

LI C G,DONG S.The trends of characterizing 3D surface microtopography[J].China Mechanical Engineering,2000,11(5):488-492.

[10] 張靖周.高等傳熱學(xué)[M].北京:科學(xué)出版社,2009.

ZHANG J Z.Advanced heat transfer[M].Beijing:Science Press,2009.