李永康 王 軍 廉自生

1.太原理工大學(xué)機(jī)械與運(yùn)載工程學(xué)院，太原，030024

2.煤礦綜采裝備山西省重點(diǎn)實(shí)驗(yàn)室，太原，030024

0 引言

斯特封在高壓、高速工況下表現(xiàn)出良好的密封特性，將其引入柱塞泵高壓高速柱塞副，研究其在不同工況下的摩擦與泄漏特性，對提高柱塞副的密封性能具有重要意義。為揭示往復(fù)密封的潤滑與密封機(jī)理，盡可能降低泄漏量、減小摩擦力，許多學(xué)者采用不同方法對多種密封結(jié)構(gòu)進(jìn)行了研究。蔡智媛等[1]基于ANSYS建立了O形往復(fù)密封圈的有限元模型，分析了結(jié)構(gòu)參數(shù)對密封圈最大von Mises應(yīng)力和最大平面剪應(yīng)力的影響規(guī)律。JOHANNESSON[2]、歐陽小平等[3]利用IHL流體動力逆解法研究了往復(fù)密封的摩擦力、泄漏量等密封性能參數(shù)。SALANT等[4]基于混合彈流潤滑理論研究了表面粗糙度對U形圈往復(fù)密封性能的影響。王冰清等[5]基于軟彈流潤滑模型研究了格萊圈的往復(fù)密封性能，結(jié)果表明：較小的密封粗糙度具有較好的密封性能，增大密封壓力會導(dǎo)致摩擦力和泄漏量增大。XIANG等[6]結(jié)合IHL理論和混合彈流潤滑理論提出了一種收斂速度更快的混合彈流潤滑計(jì)算模型，然而，IHL流體動力逆解法未考慮往復(fù)運(yùn)動中油膜流體壓力與密封圈彈性變形的流-固耦合作用，以及摩擦副表面粗糙度和空化等因素，使計(jì)算摩擦力和泄漏量與實(shí)際存在偏差。此外，采用混合彈流潤滑理論的研究大多集中于O形密封、VL型密封等密封結(jié)構(gòu)，對應(yīng)用于高壓高速柱塞副的斯特封往復(fù)密封性能的研究鮮有涉及。

本文綜合考慮密封圈表面粗糙度和流體空化效應(yīng)，建立斯特封的混合彈流潤滑模型，利用ANSYS仿真得到密封圈在安裝和密封壓力作用下的接觸壓力分布規(guī)律。在此基礎(chǔ)上，基于雷諾方程和接觸方程，采用MATLAB程序計(jì)算得到密封壓力和柱塞桿往復(fù)速度對斯特封的摩擦力、泄漏量及油膜厚度分布的影響。最后，搭建往復(fù)運(yùn)動試驗(yàn)臺，驗(yàn)正所構(gòu)建模型與計(jì)算結(jié)果的準(zhǔn)確性。

1 柱塞副斯特封的密封模型

1.1 柱塞副斯特封的結(jié)構(gòu)

斯特封一般由超彈性O(shè)形圈和改良的聚四氟乙烯階梯形滑環(huán)組成，將斯特封引入高壓柱塞泵高速柱塞副的結(jié)構(gòu)如圖1所示。斯特封位于與柱塞配合的缸套溝槽內(nèi)，斯特封兩端有高壓擋圈，用于減小密封區(qū)域的密封壓力，提高密封性能。根據(jù)柱塞與缸套的相對運(yùn)動方向，將柱塞的運(yùn)動分為內(nèi)行程運(yùn)動與外行程運(yùn)動，柱塞向低壓側(cè)運(yùn)動定義為外行程運(yùn)動，此時的柱塞泵吸液，反之則定義為內(nèi)行程，此時的柱塞泵排液。

圖1 斯特封應(yīng)用于高壓高速柱塞副示意圖

1.2 斯特封的接觸力學(xué)分析

物理模型中，缸套與柱塞桿材料為不銹鋼，O形圈材料為橡膠，擋圈材料為聚四氟乙烯(PTFE)。有限元模型中，缸套、柱塞桿與擋圈采用Linear Isotropic材料模型，其中，缸套、柱塞桿的彈性模量為210 GPa，泊松比為0.3；PTFE滑環(huán)的彈性模量為600 MPa，泊松比為0.4。描述O形圈橡膠高度非線性的多種模型中，應(yīng)用最為廣泛的是雙參數(shù)的Mooney-Rivlin超彈性模型[1,5-7]：

W=C10(I1-3)+C01(I2-3)+1/[d(J-1)2]

(1)

式中，W為應(yīng)變能密度函數(shù)，I1、I2為應(yīng)變張量不變量；材料常數(shù)C10=0.432 95，C01=1.515 08；d、J為反映材料變形程度的參數(shù)，d=0.001 027。

鑒于超彈性體橡膠材料的高度非線性，以及求解接觸壓力分布的復(fù)雜性，本文采用ANSYS研究密封接觸壓力，具體步驟如下：

(1)有限元建模。本文建立的斯特封有限元模型的幾何參數(shù)為Hallite公司的R16ML00250型密封幾何參數(shù)。密封結(jié)構(gòu)、載荷和邊界條件關(guān)于柱塞均呈軸對稱，利用ANSYS建立軸對稱的斯特封有限元模型，如圖2a所示。

(2)網(wǎng)格劃分。采用ANSYS中的Quad/Free自動網(wǎng)格劃分，尺寸設(shè)置0.1 mm。為提高計(jì)算精度與迭代收斂性，在滑環(huán)邊緣與柱塞桿可能的接觸部位進(jìn)行局部網(wǎng)格加密，并進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證。

(3)加載載荷及邊界條件。利用ANSYS進(jìn)行有限元仿真時，載荷及邊界條件分兩個載荷步加載：①柱塞桿向上移動一定位移(模擬安裝預(yù)壓縮量)；②在密封圈流體側(cè)的未接觸區(qū)施加均布載荷(模擬密封壓力的軸向壓縮)。通過對比施加均布載荷邊界節(jié)點(diǎn)上的密封壓力與接觸壓力，確定施加均布載荷的邊界節(jié)點(diǎn)。當(dāng)密封壓力大于接觸壓力時，施加均布載荷的邊界節(jié)點(diǎn)往空氣一側(cè)逐點(diǎn)迭代滲透。均布載荷邊界節(jié)點(diǎn)上的密封壓力與接觸壓力相差最小時，該節(jié)點(diǎn)為密封壓力施加的邊界點(diǎn)。缸套約束條件為所有節(jié)點(diǎn)位移都為零。

(a)網(wǎng)格劃分圖

斯特封在過盈安裝和密封壓力的組合作用下，O形圈產(chǎn)生大變形而擠壓滑環(huán)，滑環(huán)與柱塞的接觸界面產(chǎn)生接觸壓力，其分布如圖2b所示，圖中，L為接觸長度。從圖2b中可看到，斯特封在密封接觸區(qū)的接觸壓力急劇升高，沿外行程方向的接觸線壓力逐漸降低。

2 混合潤滑計(jì)算模型

斯特封的混合潤滑模型建立的基本假設(shè)如下：①柱塞往復(fù)運(yùn)動且具有光滑表面，斯特封表面粗糙且保持靜止；②密封接觸面粗糙峰的剪切徑向變形不影響其周向變形；③為能計(jì)算內(nèi)行程的內(nèi)泄漏量，模型假設(shè)斯特封低壓一側(cè)始終充滿液體[5-6]。

2.1 流體力學(xué)控制方程

基于上述假設(shè)，接觸密封面的流體力學(xué)作用規(guī)律采用考慮界面粗糙度和空化效應(yīng)的雷諾方程[4-6]：

(2)

(3)

截?cái)嘤湍ず穸菻T可表示為

(4)

其中，erf(·)為高斯誤差函數(shù)。

2.2 接觸力學(xué)模型

密封圈表面粗糙度與柱塞的相互作用可采用Greenwood-Williamson接觸模型來表示[7]：

(5)

(6)

采用非接觸式光學(xué)顯微鏡測量得到滑環(huán)表面均方根粗糙度值σ=3.11 μm，如圖3所示。

(a)三維輪廓

密封區(qū)域內(nèi)密封圈的變形量非常小，根據(jù)小變形理論，密封區(qū)任一點(diǎn)的油膜厚度的變化等于接觸區(qū)上其他載荷在該點(diǎn)引起的位移之和[4,6，8-10]，則密封區(qū)任一點(diǎn)油膜厚度為

(7)

p=pf+pc-psc

其中，Hs為靜態(tài)油膜厚度，由靜態(tài)接觸壓力psc代入式(4)線性回歸得到；Iik為由有限元分析軟件ANSYS計(jì)算得到[4,8]的影響系數(shù)矩陣I元素，表示第k個節(jié)點(diǎn)上施加單位載荷時，第i個節(jié)點(diǎn)上產(chǎn)生的變形。

2.3 計(jì)算流程

本文采用有限體積法把雷諾方程離散成代數(shù)方程組，采用基于Gauss消元法的Thomas算法(TDMA算法[11])求解上述代數(shù)方程組，得到油膜流體壓力與油膜厚度。該潤滑模型耦合了彈性變形、油膜流體壓力和粗糙度接觸壓力等因素，能更好地反映斯特封與柱塞桿接觸面的潤滑性能。雷諾方程經(jīng)離散后，轉(zhuǎn)化為如下形式的代數(shù)方程[7]：

AiΦi=BiΦi+1+CiΦi-1+Di

(8)

其中，Ai、Bi、Ci、Di(i=1,2,…,n)為離散方程的計(jì)算系數(shù)，定義C1=Bn=0。TDMA算法的求解過程分為消元和回代兩步。消元時，從系數(shù)矩陣的第二行起，逐一將每行中的非零元素消去一個，使原來的三元方程轉(zhuǎn)化為二元方程。消元進(jìn)行到最后一行時，該二元方程就轉(zhuǎn)化為一元方程，可立即得出未知量的值。然后逐一往前回代，由各二元方程解出其他未知值。消元的目的是要得到遞推方程：

Φi-1=Pi-1Φi+Qi-1

(9)

將式(9)兩邊同時乘以Ci并加上式(8)，得到

AiΦi+CiΦi-1=BiΦi+1+CiΦi-1+Di+CiPi-1Φi+CiQi-1

(10)

整理得

(11)

要計(jì)算Pi和Qi，就要計(jì)算Pi-1和Qi-1，最終要求P1和Q1。P1和Q1可由左端點(diǎn)的離散方程確定：

A1Φ1=B1Φ2+C1Φ0+D1

(12)

而C1Φ0=0，所以

(13)

當(dāng)消元進(jìn)入最后一行時，有

Φn=PnΦn+1+Qn

(14)

而PnΦn+1=0，所以

Φn=Qn

(15)

由式(15)出發(fā)，利用式(8)、式(11)及式(13)逐個回代，得出Φi(i=n-1,n-2,…, 1)，得到油膜壓力，進(jìn)而得到油膜厚度和粗糙度接觸壓力，最后計(jì)算得到反映摩擦特性的摩擦力和密封特性的泄漏量。上述過程在MATLAB軟件中采用迭代解法求解，具體計(jì)算流程如圖4所示。

圖4 數(shù)值計(jì)算流程圖

流體平均剪應(yīng)力

(16)

式中，φf、φfs、φfp為剪切流動因子[12]；μ為潤滑油黏度；v為柱塞速度。

微凸峰接觸剪應(yīng)力

τc=-fpcu/|u|

(17)

式中，f為經(jīng)驗(yàn)?zāi)Σ烈驍?shù)。

總摩擦力包含流體的黏性摩擦力和粗糙度微凸峰之間的接觸摩擦力：

(18)

式中，D為柱塞桿直徑。

泄漏率可表示為

(19)

2.4 模型驗(yàn)證

柱塞副結(jié)構(gòu)與計(jì)算參數(shù)如下：柱塞桿直徑D=25 mm，常壓黏度μ0=0.0396 Pa·s，壓黏系數(shù)α=2108Pa-1[4-6,10,13]，密封壓力ps=40 MPa，經(jīng)驗(yàn)?zāi)Σ烈驍?shù)f=0.036。

摩擦力易于測量，本文通過試驗(yàn)往復(fù)密封摩擦力與數(shù)值計(jì)算的摩擦力進(jìn)行對比，驗(yàn)證數(shù)值計(jì)算模型的可靠性。圖5a所示的高壓缸套組件由工作腔、透蓋、驅(qū)動單元、力傳感器和位移傳感器等部件組成。高壓缸套內(nèi)的2個斯特封對稱安裝，柱塞運(yùn)動過程中，始終有一個密封圈處于內(nèi)行程，一個密封圈處于外行程。圖5b所示的往復(fù)密封試驗(yàn)臺主要由變頻電機(jī)、手動泵、臥式柱塞泵、高壓往復(fù)密封裝置與蓄能器等組成。實(shí)驗(yàn)測試柱塞摩擦力的傳感器型號為JSG2000，激光位移傳感器型號為CD22-100VM12。

(a)缸套組件

變頻器調(diào)節(jié)電機(jī)轉(zhuǎn)速，通過三柱塞臥式柱塞泵將旋轉(zhuǎn)運(yùn)動轉(zhuǎn)化為柱塞的往復(fù)運(yùn)動；采用手動泵對高壓缸套系統(tǒng)加壓，模擬不同的密封壓力工況，蓄能器可在一定時間內(nèi)使得密封腔體的壓力保持恒定；激光位移傳感器測量柱塞位移，計(jì)算其加速度；將力傳感器采集到的載荷減去柱塞桿慣性力即可得到密封圈摩擦力。

圖6所示為柱塞往復(fù)速度v=28.6 mm/s時，不同密封壓力下的斯特封與柱塞桿間的摩擦力計(jì)算結(jié)果與試驗(yàn)結(jié)果。由圖6可以看到，內(nèi)外行程中，密封所受摩擦力隨密封壓力的增大而明顯增大，且內(nèi)行程摩擦力略微大于外行程摩擦力。密封接觸區(qū)總摩擦力的試驗(yàn)結(jié)果與計(jì)算結(jié)果的平均誤差小于3.8%，表明采用混合潤滑模型研究斯特封的摩擦與密封特性具有較高的可靠性。

圖6 試驗(yàn)與仿真摩擦力

3 摩擦與密封性能研究

研究接觸型往復(fù)動密封的密封特性時，必須考慮摩擦副的摩擦特性。接觸界面油膜太薄，泄漏量較小，但摩擦力過大不僅影響柱塞桿響應(yīng)速度，還會加速密封件的摩擦磨損。摩擦產(chǎn)生的熱量若不能及時擴(kuò)散，會嚴(yán)重縮短密封圈使用壽命。若接觸界面的油膜厚度較大，則摩擦磨損減小，但泄漏量會明顯增大，導(dǎo)致柱塞泵的容積效率顯著下降。因此，不同工藝參數(shù)下往復(fù)密封的摩擦力與泄漏量是其摩擦與密封性能的重要指標(biāo)。外行程中，受柱塞桿速度影響，泄漏量只包含向腔外的Couette流，構(gòu)成外行程的外泄漏，始終為正；內(nèi)行程中，除了速度引起向腔內(nèi)的剪切流，還存在由密封梯度引起的向腔外的Poiseuille流，泄漏量為剪切流與壓差流的流量差值(泄漏值為正表示內(nèi)泄漏，為負(fù)表示外泄漏)。整個柱塞副的凈泄漏量為外行程泄漏量減去內(nèi)行程泄漏量。

圖7所示為柱塞副在內(nèi)外行程中時的靜態(tài)接觸壓力psc、粗糙度接觸壓力pc、油膜流體壓力pf和量綱一油膜厚度H沿接觸線的分布。由圖7可知，內(nèi)外行程中的量綱一油膜厚度均小于4，表明該條件下的斯特封處于混合潤滑狀態(tài)，因此，有必要開展基于混合潤滑條件的斯特封摩擦與密封特性的研究。

(a)內(nèi)行程(ps=30 MPa,σ=1.2 μm,v=0.6 m/s)

內(nèi)行程中，油膜流體壓力從流體一側(cè)到空氣一側(cè)逐漸減小，靜態(tài)接觸壓力與粗糙度接觸壓力呈現(xiàn)出相似的分布，即先迅速增大，后逐漸減小，這表明油膜流體壓力與粗糙度接觸壓力共同影響密封性能，但粗糙度接觸壓力占主導(dǎo)。外行程中，密封壓力基本為零，靜態(tài)接觸壓力與粗糙度接觸壓力完全一致，這表明外行程的密封性能主要由表面粗糙度決定，而流體動力效應(yīng)基本可以忽略。

3.1 密封壓力的影響

圖8 不同密封壓力油膜厚度分布(σ=2.4 μm,v=0.6 m/s)

由圖8可知，在密封接觸區(qū)，油膜厚度先迅速減小，后沿接觸區(qū)先緩慢增大，直至空氣側(cè)時迅速增大。此外，隨著密封壓力的升高，最小油膜厚度無明顯差異，但潤滑油膜接觸區(qū)長度明顯增大，該結(jié)果與PENG等[13]得到的結(jié)果相類似。

圖9給出了不同密封壓力下，往復(fù)密封內(nèi)外行程中密封摩擦力與泄漏量的變化規(guī)律。由圖9可知，密封壓力的升高導(dǎo)致內(nèi)行程摩擦力近似線性增大，但外行程則無明顯變化。這是由于內(nèi)行程中，靜態(tài)接觸壓力隨密封壓力的升高而升高；外行程中，密封壓力為常壓，摩擦力保持不變。

(a)摩擦力

此外，隨著密封壓力的升高，內(nèi)行程泄漏量增大，并表現(xiàn)為外泄漏，外行程泄漏量無變化。這是由于密封壓力的增大使得內(nèi)行程向外的Poiseuille流增加，但向內(nèi)的Couette流無明顯變化，使得內(nèi)行程泄漏量隨密封壓力的升高而增大；而外行程時只有柱塞運(yùn)動引起向外的Couette流，基本不受密封壓力的影響，所以外行程泄漏量保持不變。

3.2 往復(fù)速度的影響

柱塞泵的高速化直接體現(xiàn)在柱塞往復(fù)速度(影響密封性能的重要因素)的升高。由圖10可知，隨著往復(fù)速度的升高，油膜厚度增大且接觸線分布基本一致。這是由于往復(fù)速度的增大增強(qiáng)了流體動壓效應(yīng)，使油膜厚度增大。

圖10 不同往復(fù)速度下的油膜厚度分布(ps=30 MPa,σ=1.2 μm)

(a)摩擦力

由圖11可知，隨著往復(fù)速度增大，內(nèi)行程中的摩擦力線性減小，外行程中的摩擦力基本無變化。這是由于內(nèi)行程中，往復(fù)速度的增大增強(qiáng)了流體的動壓效應(yīng)，使油膜厚度增大、摩擦力減小。外行程中的摩擦力主要源于密封摩擦副粗糙峰接觸，速度變化引起的摩擦力基本保持不變。

柱塞速度升高時，外行程泄漏量基本無變化。內(nèi)行程的泄漏量存在某臨界速度vc，速度v

4 結(jié)論

(1)斯特封與柱塞的接觸界面處于混合潤滑狀態(tài)，油膜壓力與粗糙度接觸壓力共同影響斯特封的密封性能，但粗糙度接觸摩擦占主導(dǎo)。計(jì)算摩擦力與試驗(yàn)?zāi)Σ亮φ`差小于3.8%，采用混合潤滑模型能更好地模擬高壓流體作用下斯特封的密封特性。

(2)隨著密封壓力的升高，油膜接觸長度與內(nèi)行程的摩擦力明顯增大，但外行程摩擦力基本無變化。隨著密封壓力的升高，內(nèi)行程泄漏量和凈泄漏量增大，而外行程時泄漏量無明顯變化。

(3)在剪切流與壓差流綜合作用下，柱塞副的臨界往復(fù)速度使內(nèi)行程中的泄漏方向發(fā)生改變。當(dāng)往復(fù)速度小于臨界值時，隨著往復(fù)速度的增大，外泄漏明顯減?。划?dāng)往復(fù)速度大于臨界值時，向內(nèi)的剪切流量大于向外的壓差流量，內(nèi)行程泄漏量為正值，且隨著速度的增加明顯減小，表現(xiàn)出斯特封較好的高速密封性。