基于分形理論的機(jī)械密封干摩擦?xí)r端面接觸特性的研究

陳 志，蔡 垚，顧燦鴻

(四川大學(xué) 化學(xué)工程學(xué)院，四川 成都 610065)

機(jī)械密封為端面密封，其依賴兩個(gè)固體表面間的接觸和相對(duì)運(yùn)轉(zhuǎn)，以防止流體泄漏。常規(guī)設(shè)計(jì)中通常按照接觸端面的名義面積來計(jì)算其端面比壓。然而，機(jī)械加工表面在本質(zhì)上都是粗糙的，動(dòng)、靜環(huán)端面的實(shí)際接觸僅發(fā)生在部分微凸體上，導(dǎo)致實(shí)際接觸面積僅占名義接觸面積很小的一部分[1]。干摩擦是機(jī)械密封最惡劣的潤滑狀態(tài)，其接觸狀態(tài)直接決定了密封的性能與壽命[2]。張永振[3]指出在干摩擦狀態(tài)下，摩擦熱量與速度、接觸壓力存在冪函數(shù)的關(guān)系，外載荷和速度對(duì)端面的接觸、摩擦磨損行為的影響尤為突出。所以，研究外載荷和速度對(duì)機(jī)械密封干摩擦?xí)r端面接觸特性的影響對(duì)了解其密封機(jī)理和端面磨損具有重要意義。

目前，對(duì)相對(duì)運(yùn)動(dòng)摩擦端面的動(dòng)態(tài)過程的研究仍然缺乏直接觀察和檢測評(píng)估的手段，因此，為了準(zhǔn)確描述摩擦端面的接觸特性，需要建立合理的接觸模型，前人對(duì)此進(jìn)行了大量研究。Majumdar等[4–5]首先基于2維W–M分形函數(shù)建立了彈塑性接觸分形模型，探究了總載荷與真實(shí)接觸面積的關(guān)系式；Wang等[6–7]在M–B模型的基礎(chǔ)上進(jìn)行改進(jìn)，得到了M–B分形接觸修正模型；孫見君[8]、魏龍[9]、丁雪興[10]等參照M–B模型建立了機(jī)械密封摩擦端面的分形模型，通過2維理論計(jì)算對(duì)密封端面泄漏、真實(shí)接觸面積、溫度分布進(jìn)行了研究?；诜中卫碚摻⒌慕佑|模型在各個(gè)領(lǐng)域都得到了廣泛運(yùn)用[11–16]。可見，利用分形理論可以實(shí)現(xiàn)對(duì)粗糙面的唯一表征，不受尺度的影響。

干摩擦的宏觀表現(xiàn)行為是眾多接觸摩擦學(xué)行為耦合作用的結(jié)果[3]。粗糙表面上的接觸壓力、熱應(yīng)力、溫度場是相互影響的，因此，考慮微凸體間的相互作用及熱–力耦合對(duì)分析摩擦端面的接觸特性尤為重要。孫星星等[17]考慮密封動(dòng)靜環(huán)的熱力變形，數(shù)值求解了密封環(huán)的溫度分布，但忽略了粗糙面微凸體的影響。房桂芳等[18]依據(jù)分形理論提出了機(jī)械密封端面摩擦熱的耦合計(jì)算方法。Li等[19]提出了考慮磨損和熱變形的接觸機(jī)械密封分形泄漏模型。黃健萌[20–21]、丁雪興[22]等建立了粗糙面/平面的1維滑移接觸模型，通過有限元模擬求解了滑動(dòng)摩擦過程中微凸體變形及熱–力耦合對(duì)粗糙面應(yīng)力、溫度變化的影響。總的來說，對(duì)于機(jī)械密封摩擦端面接觸的研究主要在理論計(jì)算方面，模擬計(jì)算也僅采用了直線滑動(dòng)摩擦的方式；對(duì)于相互接觸的表面而言，直線運(yùn)動(dòng)與回轉(zhuǎn)運(yùn)動(dòng)的摩擦學(xué)特性是不同的[3]，因此采用3維回轉(zhuǎn)接觸模型進(jìn)行仿真將更符合實(shí)際。

針對(duì)碳石墨與碳化鎢配對(duì)的機(jī)械密封，首先，作者基于W–M分形函數(shù)建立了模擬其摩擦副端面的3維粗糙實(shí)體與理想光滑剛體的轉(zhuǎn)動(dòng)摩擦接觸模型；然后，在考慮熱–力耦合及微凸體間相互作用的基礎(chǔ)上，通過ABAQUS有限元分析軟件求解了其瞬時(shí)干摩擦?xí)r密封端面間的接觸特性，并討論不同外載荷、滑動(dòng)速度對(duì)其特性的影響，并針對(duì)API682中串聯(lián)式機(jī)械密封的第2級(jí)接觸式干摩擦機(jī)械密封（抑制機(jī)械密封）進(jìn)行了分析，以期為探索機(jī)械密封在干摩擦?xí)r的發(fā)熱、摩擦行為和密封設(shè)計(jì)提供借鑒。

1 計(jì)算模型的建立

1.1 粗糙表面及轉(zhuǎn)動(dòng)接觸模型的建立

表征粗糙接觸界面的形貌是摩擦學(xué)的熱點(diǎn)，Yan等[23]的研究給出了能夠描述粗糙表面微觀特征的3維W–M函數(shù)：式中：Z(x，y)為3維粗糙表面形貌的高度；D和G分別為分形維數(shù)和分形粗糙度；L為取樣長度；γ為表面輪廓頻率密度的特征參數(shù)，通常取1.5；M為曲面褶皺的重疊數(shù)，值為10；φm,n為隨機(jī)相位；nmax為最大頻率，其值為：

式中，Ls為截止長度（近似取材料的原子間距）。

對(duì)于同一粗糙表面，Z(x，y)可由與尺度無關(guān)的分形參數(shù)D與G唯一確定，因此式（1）具有在不同長度尺度上描述粗糙表面的能力[23]，并且隨著D增大，G減小，表面粗糙度降低，趨于光滑[24–25]。本文中，D取2.35，G取2.45×10–9m，此時(shí)表面粗糙度Ra的大小為0.8 μm。

由于接觸問題屬于典型的邊界非線性問題，且粗糙表面的特征尺度在微米級(jí)，考慮到計(jì)算成本，研究中將微元體的模型尺寸也控制在微米級(jí)，其中，靜環(huán)模型：外徑DAi=260 μm，內(nèi)徑DAo=380 μm，厚度hA=60 μm；動(dòng)環(huán)模型：外徑DBi=250 μm，內(nèi)徑DBo=390 μm，厚度hB=60 μm。通過增大轉(zhuǎn)速的數(shù)量級(jí)可以實(shí)現(xiàn)模型微元體的真實(shí)線速度V，根據(jù)文獻(xiàn)[26]粗糙表面輪廓曲線的局部與整體的自相似性可知，微元體的接觸行為可以反映實(shí)際端面之間的接觸行為。

機(jī)械密封中，動(dòng)、靜環(huán)材料常采用硬軟配對(duì)，摩擦副為硬質(zhì)合金（WC）對(duì)碳石墨（M106K），兩者的彈性模量相差1～2個(gè)數(shù)量級(jí)[22]，見表1。因此，研究中將靜環(huán)（實(shí)體A）與動(dòng)環(huán)（實(shí)體B）的接觸簡化為分形粗糙表面與剛性理想光滑平面的接觸[9]；同時(shí)，由于密封環(huán)模型軸對(duì)稱，為了減少計(jì)算時(shí)間，靜環(huán)采用1/8周期模型，如圖1所示。其中：實(shí)體A中，A2為頂面，A1為底面（分形粗糙面），A3、A4為側(cè)面，A5為內(nèi)圓柱面，A6為外圓柱面；實(shí)體B中，B1為頂面（接觸面），B2為底面，B3為內(nèi)圓柱面，B4為外圓柱面；RP為參考點(diǎn)；徑向X、周向Y、軸向Z（實(shí)體A、B厚度方向）為坐標(biāo)系。

表1 摩擦副材料參數(shù)[22]Tab. 1 Material parameters of friction pair[22]

圖1 分形接觸模型Fig. 1 Fractalcontact model

1.2 位移邊界條件

計(jì)算時(shí)設(shè)置兩個(gè)時(shí)間段：在第1個(gè)時(shí)間段內(nèi)，將均布載荷由0線性增加到P，加載到實(shí)體A的A2面上，同時(shí)對(duì)其4個(gè)側(cè)面施加沿X、Y方向的位移約束。對(duì)于實(shí)體B，約束其Z方向的位移，同時(shí)將實(shí)體B用參考點(diǎn)RP（0，0，0）進(jìn)行剛性約束（設(shè)置實(shí)體B為剛性體，同時(shí)可以通過RP的運(yùn)動(dòng)控制實(shí)體B的運(yùn)動(dòng)），通過約束RP的6個(gè)自由度使B保持靜止，讓兩實(shí)體建立穩(wěn)定的接觸關(guān)系。在第2個(gè)時(shí)間段內(nèi)，保持A2面上載荷P不變，同時(shí)給定實(shí)體B的參考點(diǎn)RP繞Z軸的轉(zhuǎn)動(dòng)速度ω，使兩實(shí)體間發(fā)生轉(zhuǎn)動(dòng)摩擦。

1.3 模型假設(shè)

由于計(jì)算條件的限制，本文在計(jì)算中做出如下假設(shè)：

1）滑動(dòng)摩擦過程遵循庫侖定律，在計(jì)算時(shí)間范圍內(nèi)認(rèn)為摩擦系數(shù)保持不變，不受外載荷、滑動(dòng)速度變化的影響。

2）構(gòu)成實(shí)體A、B的材料各向同性，且在瞬態(tài)計(jì)算中，由于溫升不大，假設(shè)材料的熱物性參數(shù)不會(huì)改變。

3）忽略材料磨損的影響，認(rèn)為摩擦消耗的功全部轉(zhuǎn)化為熱能，并被摩擦副吸收，則局部熱流密度可表示為：

式中，q為熱流密度，f為摩擦系數(shù)，Pb為粗糙表面微元體的接觸壓力，V為相對(duì)滑動(dòng)速度，t為時(shí)間，r為半徑。

4）考慮摩擦熱在兩實(shí)體間的分配，認(rèn)為接觸區(qū)域?yàn)槔硐霟醾鲗?dǎo)，即接觸點(diǎn)處兩實(shí)體的瞬時(shí)溫度相等，這樣摩擦熱流密度就能根據(jù)材料的熱物性參數(shù)及散熱環(huán)境在兩實(shí)體間自由分配[20]，同時(shí)忽略熱輻射的影響。

1.4 熱邊界條件

實(shí)體A、B的熱傳導(dǎo)方程可表示為：

式中，ρ為實(shí)體A、B材料的密度，c為比熱容，k為熱導(dǎo)率，T為溫度，t為時(shí)間，r為半徑。

初始條件：TA=TB=298 K，t=0。

耦合條件（假設(shè)4））：TA=TB，qA+qB=q（實(shí)體A、B接觸點(diǎn)）。

接觸面（A1、B1）的熱邊界條件為：忽略熱輻射，非接觸區(qū)域的傳熱為對(duì)流換熱；接觸區(qū)域?yàn)槟Σ翢崃鞯妮斎隱21]：

式中，接觸點(diǎn)處g(m)=1，非接觸點(diǎn)處g(m)=0。

由于本文研究的是瞬態(tài)問題，在很短的時(shí)間內(nèi)，熱量只在兩實(shí)體接觸面的表層及亞表層傳遞，因此可以認(rèn)為密封環(huán)非接觸面（A2、B2）是絕熱的；同樣地，干摩擦條件下可以認(rèn)為實(shí)體A及實(shí)體B的側(cè)面A3～A6、B3～B4是絕熱的。

1.5 求解方法

采用ABAQUS/Explicit求解計(jì)算模型，對(duì)于復(fù)雜的接觸問題，顯示方法求解效率較高且結(jié)果可靠。其采用中心差分的方法求解運(yùn)動(dòng)方程，用一個(gè)時(shí)間段的運(yùn)動(dòng)學(xué)條件計(jì)算下一個(gè)時(shí)間段的運(yùn)動(dòng)學(xué)條件，通過自動(dòng)劃分足夠小的時(shí)間增量段來確保模型狀態(tài)穩(wěn)定，從而獲得高精度的解[27]。求解前需要設(shè)置接觸對(duì)，將剛性光滑面B1設(shè)置為主面，分形粗糙面A1設(shè)為從面，采用罰函數(shù)法作為求解接觸問題的算法，兩表面間的相對(duì)滑動(dòng)通過有限滑移公式描述。

考慮熱–力耦合計(jì)算，使用C3D8T六面體熱–力耦合單元對(duì)整個(gè)計(jì)算模型進(jìn)行網(wǎng)格劃分，并對(duì)實(shí)體A的粗糙面（A1）的表層及亞表層進(jìn)行網(wǎng)格細(xì)化，對(duì)光滑實(shí)體B使用較粗的網(wǎng)格劃分方式。為了證明求解結(jié)果與網(wǎng)格大小無關(guān)，對(duì)模型進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證分析，由于主要研究對(duì)象是粗糙表面，所以采用調(diào)整粗糙表面網(wǎng)格大小的方式驗(yàn)證接觸面積的變化。定義無量綱接觸面積為Ar/Aa，其中，Ar為真實(shí)接觸面積，Aa為粗糙面A1的名義接觸面積。ABAQUS軟件中可設(shè)置接觸面積的監(jiān)控，計(jì)算完成后輸出Ar隨時(shí)間的變化，則無量綱接觸面積隨網(wǎng)格尺寸的變化如圖2（a）所示。可以發(fā)現(xiàn)，隨著網(wǎng)格尺寸減小，無量綱接觸面積趨于穩(wěn)定，實(shí)際粗糙面網(wǎng)格尺寸為1.736 μm，可以較好地保證計(jì)算結(jié)果的準(zhǔn)確性。網(wǎng)格劃分完成后，光滑實(shí)體單元數(shù)目為10 400，粗糙實(shí)體單元數(shù)為20 475，如圖2（b）、（c）所示。

圖2 網(wǎng)格劃分Fig. 2 Mesh generation

2 算例分析

2.1 計(jì)算參數(shù)選擇

動(dòng)環(huán)、靜環(huán)的材料如表1所示，兩實(shí)體的干摩擦系數(shù)選擇為f=0.25[28]，靜環(huán)材料的屈服強(qiáng)度為σy=200 MPa[9]。在初期計(jì)算中，外載荷P過小時(shí)，數(shù)據(jù)變化不明顯，故P選取變化范圍為0.2～40.0 MPa；滑動(dòng)速度V為光滑實(shí)體B外徑處的線速度，變化范圍為20～50 m/s。

在模擬計(jì)算過程中不考慮材料非線性問題，只考慮幾何非線性與接觸邊界非線性，由于計(jì)算過程難以收斂，同時(shí)粗糙表面實(shí)際處于不斷磨損的狀態(tài)中，表面形貌也不斷更新[24]，因此，只求解轉(zhuǎn)動(dòng)發(fā)生后60 μs以內(nèi)的瞬態(tài)過程。為了使兩實(shí)體發(fā)生平穩(wěn)接觸，將轉(zhuǎn)動(dòng)前的計(jì)算時(shí)間設(shè)置為40 μs，外載荷以線性增加的方式施加在A2面上直到載荷等于P為止，之后保持載荷不變使兩實(shí)體發(fā)生相對(duì)轉(zhuǎn)動(dòng)，速度為V；密封環(huán)初始溫度與環(huán)境溫度設(shè)置為298 K。每次計(jì)算改變外載荷P、滑動(dòng)速度V兩個(gè)參數(shù)中的一個(gè)值，計(jì)算得到密封環(huán)不同工況下的7組摩擦數(shù)據(jù)。

2.2 結(jié)果分析與討論

2.2.1 真實(shí)接觸面積

定義無量綱載荷為P/E′，其中：

得到本文分形模型與魏龍分形模型[9]、M–B分形模型[5]理論值的對(duì)比如圖3所示。

圖3 本文分形模型與其他模型理論值的對(duì)比Fig. 3 Comparison of the theoretical values of fractal model in this paper with other models

從圖3中可以發(fā)現(xiàn)，隨著外載荷增大，真實(shí)接觸面積近似呈線性增大，外載荷增加3倍時(shí)，接觸面積增大了2.33倍，這是因?yàn)橥廨d荷的增加導(dǎo)致粗糙表面受力增加，接觸面主峰微凸體發(fā)生彈塑性變形，次峰微凸體逐漸參與接觸，增加了接觸面積。同時(shí)，可以看出：隨著載荷的增大，模擬分析的結(jié)果與理論值之間具有相同的變化趨勢，在載荷較低時(shí)，它們間差異較?。ù?lián)式機(jī)械密封第二級(jí)接觸式干運(yùn)轉(zhuǎn)機(jī)械密封恰恰是在低載荷下工作）；隨著載荷的進(jìn)一步增大，它們間的誤差也逐漸增大（無量綱載荷為0.001 5時(shí)，M–B模型計(jì)算得到的無量綱接觸面積是本文模擬計(jì)算的1.364倍，是魏龍模型的1.186倍）。產(chǎn)生這一差異的原因在于：魏龍分形模型、M–B分形模型的建立是基于2維粗糙輪廓的，而不是真實(shí)的3維粗糙表面，它們未考慮粗糙表面相鄰微凸體之間的相互作用。文獻(xiàn)[29]的研究表明，魏龍分形模型、M–B分形模型計(jì)算得到的真實(shí)接觸面積與實(shí)驗(yàn)值相比較大。本文所采用的3維有限元模擬考慮了粗糙表面相鄰微凸體之間的相互作用，所以其結(jié)果更接近實(shí)際，從而驗(yàn)證了模擬分析的正確性。

圖4為P=30 MPa時(shí)，不同滑動(dòng)速度下，滑動(dòng)0.3 μs及60.0 μs時(shí)粗糙面的無量綱接觸面積?？梢园l(fā)現(xiàn)，不同滑動(dòng)時(shí)間下，接觸面積隨速度增加的斜率相同，速度越大，接觸面積越大。但在滑動(dòng)速度增加1.5倍的情況下，接觸面積僅增大為原來的1.002倍，因此，粗糙面的真實(shí)接觸面積的大小主要取決于外載荷，而滑動(dòng)速度影響較小。

圖4 滑動(dòng)階段無量綱接觸面積隨滑動(dòng)速度的變化Fig. 4 Non-dimensional contact area vs sliding velocity when ring sliding

2.2.2 接觸壓力

摩擦磨損是表面行為，因此，研究粗糙表面的真實(shí)接觸面積、接觸壓力、溫度的分布極其重要。圖5（a）～（d）所示為V=30 m/s時(shí)，P=10、20、30、40 MPa下，滑動(dòng)60 μs時(shí)刻分形粗糙面的接觸壓力分布云圖。由圖5可以看出：隨著外載荷增大，真實(shí)接觸面積Ar上較高的接觸壓力節(jié)點(diǎn)在增加，粗糙表面上最大接觸壓力值位于節(jié)點(diǎn)13538；當(dāng)外載荷為10、20、30、40 MPa時(shí)，節(jié)點(diǎn)13538的接觸壓力分別為415.88、416.64、423.25和431.59 MPa。對(duì)比不同載荷下各個(gè)接觸區(qū)域上接觸壓力的變化，可以發(fā)現(xiàn)外載荷增大3倍后，節(jié)點(diǎn)13538上接觸壓力值增加的幅度并不是很大，這說明相對(duì)運(yùn)動(dòng)的摩擦副的接觸壓力峰值達(dá)到一定值后變化緩慢，該值對(duì)外載荷的變化不敏感，應(yīng)與摩擦副的材料力學(xué)特性有關(guān)。

粗糙表面節(jié)點(diǎn)13538的接觸壓力隨載荷增加變化緩慢，其原因是接觸壓力在達(dá)到一定值后，材料發(fā)生彈塑性變形，次峰微凸體進(jìn)入接觸并承擔(dān)了額外的載荷，使得節(jié)點(diǎn)13538所在微凸體受力得到緩解。表現(xiàn)為粗糙表面的接觸壓力峰值維持在一個(gè)定值范圍內(nèi)，而其他接觸節(jié)點(diǎn)的接觸壓力繼續(xù)增加。換句話說，在本文的研究參數(shù)范圍內(nèi)，接觸壓力峰值是“自限性”的，這種特征保證端面在較大的外載荷下的接觸是由點(diǎn)接觸到局部面接觸、一個(gè)節(jié)點(diǎn)接觸到眾多節(jié)點(diǎn)接觸。更多的節(jié)點(diǎn)參與接觸，端面受力條件得到改善，這種“自限性”對(duì)于密封環(huán)端面微凸體受力狀況的改善是有利的。

圖5 不同外載荷下滑動(dòng)60 μs時(shí)粗糙表面接觸壓力PB云圖Fig. 5 Contourmap of contact pressure on the rough surface when sliding at 60 μs under different external loads

圖6（a）～（d）所示為P=30 MPa時(shí)，V=20、30、40、50 m/s下，滑動(dòng)60 μs時(shí)刻分形粗糙面的接觸壓力分布云圖。

由圖6可以發(fā)現(xiàn)，在外載荷相同時(shí)，接觸壓力的分布基本相同，最大接觸壓力分別為430.34、423.25、428.86、420.26 MPa，均位于粗糙表面節(jié)點(diǎn)13538，并隨著速度增加具有微小的下降趨勢，這是因?yàn)榛瑒?dòng)期間速度增大時(shí)引起的微凸體局部摩擦熱流密度變大，微凸體熱變形更大使接觸節(jié)點(diǎn)局部接觸面積略微增大，導(dǎo)致節(jié)點(diǎn)最大接觸壓力略有減小。

圖6 不同速度下滑動(dòng)60 μs時(shí)粗糙表面接觸壓力Pb云圖Fig. 6 Contourmap of contact pressure on rough surface when sliding at 60 μs under different sliding velocity

2.2.3 溫度分布

圖7（a）～（d）所示為V=30 m/s，P=10、20、30、40 MPa下，滑動(dòng)60 μs時(shí)刻分形粗糙面的溫度分布云圖。由于粗糙表面微凸體的形狀大小、高低不一，滑動(dòng)接觸正是發(fā)生在這些不連續(xù)的接觸微凸體上。摩擦所引起的熱量產(chǎn)生、釋放也發(fā)生這些區(qū)域，使得接觸節(jié)點(diǎn)處溫度較周圍較高，形成“熱點(diǎn)”。

由圖7可以看出，P=10 MPa時(shí)高溫區(qū)主要分布在粗糙表面微凸體接觸點(diǎn)位置，且端面上存在明顯的溫度梯度分布，會(huì)有熱應(yīng)力產(chǎn)生。隨著載荷增大，粗糙表面接觸微凸體及接觸點(diǎn)越多，在滑動(dòng)瞬間形成的“熱點(diǎn)”也越多，接觸中心區(qū)域的溫度也越高，粗糙表面整體溫度也增高。另外，可以發(fā)現(xiàn)，不同載荷下最高溫度值所在的節(jié)點(diǎn)不同，不是單純地位于接觸壓力最大的位置。這表明粗糙表面的最高溫度還受到局部熱流密度的影響；即最高溫度受到局部接觸壓力和局部滑動(dòng)速度和局部熱傳遞的共同影響。

圖7 不同外載荷下滑動(dòng)60 μs時(shí)粗糙表面溫度云圖Fig. 7 Contourmap of temperature on rough surface when sliding at 60 μs under different external loads

隨著外載荷的增加，雖然粗糙表面最大接觸壓力改變較小，但載荷變化使得粗糙表面接觸點(diǎn)的數(shù)量明顯增多，滑動(dòng)過程中粗糙表面上的摩擦熱源增多，導(dǎo)致了端面整體溫度升高，故材料的耐摩擦磨損性能也受到影響。

圖8所示為外載荷為10 MPa、滑動(dòng)速度為30 m/s，摩擦滑動(dòng)300 μs條件下的粗糙面各節(jié)點(diǎn)的溫度變化曲線（節(jié)點(diǎn)分布如圖7（a）所示）。由圖8可以看出：0～40 μs為外載荷加載階段，粗糙面各節(jié)點(diǎn)溫度保持不變；當(dāng)兩表面開始滑動(dòng)，40～42 μs內(nèi)由于速度突增及熱傳導(dǎo)的滯后性出現(xiàn)了溫度急升；隨后42～340 μs，在自身熱傳導(dǎo)的作用下，各節(jié)點(diǎn)溫度開始緩慢上升，溫度上升斜率也是一致的，且滑動(dòng)60 μs的溫升斜率與滑動(dòng)300 μs的溫升斜率也相同。這也驗(yàn)證了本文采用摩擦滑動(dòng)60 μs的計(jì)算時(shí)間進(jìn)行研究是可行的。

圖8 摩擦滑動(dòng)300 μs內(nèi)粗糙面節(jié)點(diǎn)溫度隨時(shí)間的變化Fig. 8 Temperature change of rough surface node with time within 300 μs of friction sliding

API682中規(guī)定串聯(lián)式機(jī)械密封第2級(jí)接觸式干運(yùn)轉(zhuǎn)機(jī)械密封的介質(zhì)壓力應(yīng)控制在0.07 MPa以內(nèi)。為了闡釋該密封在有氣體冷卻的配合下實(shí)現(xiàn)25 000 h運(yùn)轉(zhuǎn)的可能性，本文模擬計(jì)算了外載荷0.2 MPa、滑動(dòng)速度30 m/s條件下滑動(dòng)60 μs時(shí)粗糙面溫度分布云圖，如圖9所示。由圖9可以看到，摩擦滑動(dòng)60 μs后，粗糙面的最高溫度僅從環(huán)境溫度25 ℃升高到25.42 ℃，溫升僅有0.42 ℃，粗糙面溫度變化不明顯。圖10為滑動(dòng)速度30 m/s時(shí)，不同外載荷下粗糙面最高溫度的變化，隨著外載荷的增加，最高溫度是線性增大的，與式（3）描述的規(guī)律相符?？梢娊档蜋C(jī)械密封的外載荷可以有效防止其干摩擦?xí)r石墨環(huán)溫升過高，API682標(biāo)準(zhǔn)中建議介質(zhì)壓力≤0.07 MPa是有必要的。

圖9 0.2 MPa下摩擦滑動(dòng)60 μs時(shí)粗糙表面溫度云圖Fig. 9 Contourmap of temperature on rough surface when sliding at 60 μs under 0.2 MPa

圖10 不同外載荷下滑動(dòng)60 μs時(shí)粗糙面最高溫度的變化Fig. 10 Change of maximum temperature on rough surface under different load when sliding at 60 μs

圖11（a）～（d）所示為P=30 MPa，V=20、30、40、50 m/s下，滑動(dòng)60 μs時(shí)刻分形粗糙面的溫度分布云圖。由圖11可以發(fā)現(xiàn)，在相同的外載荷下，對(duì)應(yīng)的粗糙表面最高溫度值分別為44.84、49.12、52.99和57.52 ℃，隨著滑動(dòng)速度增加溫度呈上升趨勢，粗糙表面的溫度不均勻趨勢明顯加劇，各區(qū)域溫度值均有所提高。根據(jù)式（3），速度增加，摩擦發(fā)熱增加，粗糙表面的溫度梯度也隨之變大，熱應(yīng)力也會(huì)加劇。

圖11 不同速度下滑動(dòng)60 μs時(shí)粗糙表面溫度云圖Fig. 11 Contourmap of temperature on rough surface when sliding at 60 μs under different sliding velocity

3 結(jié) 論

在P=0.2～40 MPa、V=20～50 m/s，f=0.25和粗糙表面Ra=0.8 μm的工況下，針對(duì)碳石墨–硬質(zhì)合金配對(duì)摩擦副在干摩擦狀態(tài)下的接觸特性進(jìn)行了分析，研究發(fā)現(xiàn)：

1）粗糙端面真實(shí)接觸面積的大小主要取決于外載荷。外載荷增加真實(shí)接觸面積近似線性增大，但滑動(dòng)速度影響較小。

2）隨著外載荷與滑動(dòng)速度的增加，石墨粗糙表面微凸體接觸壓力峰值在415～432 MPa范圍內(nèi)，表現(xiàn)出一種“自限性”。

3）粗糙端面溫度分布不均勻，局部溫度較高；端面最高溫度并非位于接觸壓力最大的位置，而是受到局部接觸壓力和局部滑動(dòng)速度及局部熱傳遞的共同影響。

4）外載荷、速度增大均會(huì)使端面整體溫度增加。

5）API682標(biāo)準(zhǔn)中要求抑制機(jī)械密封的介質(zhì)壓力≤0.07 MPa是有必要的。降低介質(zhì)壓力可以降低外載荷，從而可以有效防止干摩擦?xí)r密封環(huán)溫升過高。