李雙路 訚耀保 張鑫彬 王曉露 傅俊勇
1.同濟大學機械與能源工程學院,上海,201804 2.上海航天控制技術研究所,上海,201109
高端液壓閥按結構主要分為滑閥、錐閥、球閥和剪切閥等形式,其中,液壓滑閥的閥口具有薄壁孔口特征,閥口流量受油液黏度、溫度等因素影響較小,因此在比例閥、伺服閥等高端液壓控制元件中大量使用,其中全周邊液壓滑閥因面積梯度大、閥芯質量小、控制特性好,應用最為廣泛。液壓滑閥的形貌形性對伺服控制系統的精確控制有決定性作用,出廠時對其閥口銳邊有非常高的要求,但在服役過程中,滑閥不可避免地受到油液中顆粒物的沖蝕,造成閥口處閥芯閥套的材料流失并產生圓角,引起滑閥性能出現不可逆的演化過程。
高速流體攜帶固體粒子(顆粒物)對靶材(對應本文的閥芯、閥套)沖擊而造成材料表面流失的現象即沖蝕磨損,沖蝕磨損的理論研究始于20世紀60年代,最初研究塑性材料和脆性材料的沖蝕破壞形式。美國加州大學伯克利分校Finnie首先提出了塑性材料的微切削理論,認為當磨粒劃過靶材表面時,如同一把微型刀具將材料切除而產生沖蝕磨損,該理論適用于低攻角下塑性材料受剛性磨粒沖蝕的分析,但計算高攻角下的沖蝕磨損誤差較大;1963年殼牌公司Bitter提出了變形磨損理論,認為當粒子垂直撞擊壁面的沖擊力超過靶材的屈服強度時會造成材料發生塑性變形,產生裂紋并引起靶材的體積流失。變形磨損理論完善了在高攻角下塑性材料的沖蝕,塑性材料總的沖蝕磨損率為變形磨損和切削磨損的代數和。后來Grant、Forder、Edwards等基于Finnie和Bitter的模型提出了沖蝕磨損率的不同計算方法,拓展了沖蝕理論在不同環境下的適用范圍[1-2],其中Edwards的模型由于對沖蝕預測的精確度較高且形式簡單,被廣泛應用于氣固、液固及氣液固流動中。
沖蝕會造成液壓閥口形狀的變化。文獻[3]研究了油氣運輸中節流閥的沖蝕現象,發現選擇合適的抗蝕材料可以使閥芯壽命明顯延長;文獻[4]通過加速磨損試驗測量并研究沖蝕磨損前后伺服閥節流邊的輪廓變化,分析了節流邊壓差等因素對磨損的影響規律。在理論計算方面,文獻[5-6]通過計算流體力學方法研究了噴嘴擋板伺服閥中滑閥副節流邊、噴嘴和擋板的沖蝕磨損以及前置級沖蝕對伺服閥零偏的影響;文獻[7]研究了閥芯閥套沖蝕區域隨閥口開度的變化情況,發現來流面比回流面的沖蝕更加嚴重;文獻[8]提出了偏轉板伺服閥前置級射流盤劈尖沖蝕磨損的預測方法,實現不同油液污染度下的劈尖沖蝕形貌預測;文獻[9]建立了射流管不對中時的前置級沖蝕模型,并分析劈尖角度和沖蝕率之間的關系。上述研究主要集中在磨損位置的確定以及初始沖蝕率,但未考慮沖蝕過程的動態變化。文獻[10]通過仿真和實驗對比,發現仿真過程中忽略幾何模型變化會導致沖蝕磨損率計算不準確。文獻[11]考慮沖蝕形貌演化對沖蝕過程的影響,建立了滑閥節流銳邊沖蝕磨損深度和磨損輪廓的定量預測模型,并得到了顆粒尺寸、節流壓差、滑閥開度、液流方向等關鍵因素對滑閥沖蝕的影響。上述研究主要集中在高端液壓閥沖蝕預測的實驗探究和數值仿真技術方面,高端液壓閥沖蝕預測的工程應用尚不多見。
本文考慮顆粒物尺寸以及顆粒物撞擊閥口的概率,提出一種適用于全周邊滑閥閥口的沖蝕圓角計算模型,建立閥口沖蝕圓角同顆粒物尺寸、質量流量、撞擊速度、沖擊角度、閥口開度間的關系,進一步得到閥口流量、壓差對沖蝕過程的影響,結合閥控缸的負載和運動速度,分析滑閥沖蝕圓角的定量計算方法,并研究滑閥沖蝕過程中閥口形貌特征及性能演化規律。
圖1所示為某伺服閥滑閥結構,采用全周邊四邊滑閥形式,圖中閥芯位于左位,高壓油液從P口經過節流口1進入閥腔S1,負載輸出的油液從閥腔S2經節流口3從T口流出,此時負載處于伸出狀態;當閥芯位于右位時,1、3口關閉,2、4口打開,此時負載處于縮回狀態。由于伺服閥壓降大,故4個節流口受到高速射流的沖刷,其中的顆粒物不斷撞擊節流邊造成閥口材料的流失,導致各個閥口出現沖蝕磨損,嚴重影響伺服閥的控制精度。

圖1 全周邊四邊滑閥Fig.1 Full-circumference four-sided slide valve
為便于閥口沖蝕的分析,根據閥腔油液流動規律,作如下假設:①閥口沖蝕后的輪廓為四分之一圓弧,且油液流出閥腔與油液流入閥腔造成的沖蝕磨損相同;②假設液壓閥來流油液中顆粒物運動的角度是固定的且顆粒物撞擊壁面的角度僅與顆粒物直徑有關;③由于閥口尺寸小,故忽略顆粒物在閥口處與壁面碰撞后反彈造成的沖蝕;④流體中的顆粒物分布均勻。
在傳統的液壓閥沖蝕數值仿真中,顆粒物被視為不占據空間的質點,顆粒物撞擊閥口的概率不受顆粒物尺寸的影響,但在高端電液伺服閥中,滑閥前的過濾精度一般為10~20 μm,閥口開度常在數十微米甚至更小,兩者大小接近,顆粒物的尺寸對其撞擊閥口概率的影響不可忽略。按照節流口寬度和顆粒物尺寸的相對大小可將顆粒物通過閥口的狀態分成圖2所示的兩種情況:第一種情況如圖2a所示,節流口寬度xk較大,顆粒物直徑dp較小,部分顆粒物沒有與閥口壁面發生碰撞便直接流向下游;第二種情況如圖2b所示,節流口寬度xk較小,顆粒物直徑dp較大,顆粒物必然會撞擊到閥芯或者閥套。

(a)節流口寬度相對 顆粒物尺寸較大 (b)節流口寬度相對 顆粒物尺寸較小圖2 顆粒物通過閥口的兩種狀態Fig.2 Two states in which particles pass through thevalve port
假設顆粒物的固定運動角度為θ,閥芯和閥套的沖蝕圓角半徑分別是rs和rb,閥芯位移為xv,閥芯閥套徑向間隙為s。對于圖2a所示的第一種情況,圖中1號、5號顆粒物分別與閥芯和閥套沖蝕邊界外緣發生碰撞,2號、4號顆粒物處于碰撞的臨界點,3號顆粒物不會與閥芯和閥套發生碰撞,油液中的全部顆粒物都處于1號顆粒物至5號顆粒物范圍之內(假設還有顆粒物在1號顆粒物和5號顆粒物之外,那么閥芯和閥套的沖蝕邊界將向外擴張)。定義顆粒物撞擊到閥芯上的區域寬度是bs,撞擊到閥套上的區域寬度是bb,而在bn寬度范圍內顆粒物會直接流向下游而不與閥芯或者閥套碰撞。
顆粒物通過閥口時,主要幾何尺寸關系如下:
bb=rb(1-sinθ)+dp/2
(1)
bs=rs(1-cosθ)+dp/2
(2)
(3)
假設顆粒物分布均勻且運動角度一致,顆粒物撞擊到閥芯或閥套上的概率等于相應區域寬度占總寬度的比值。故直徑為dp的顆粒物撞擊到閥芯和閥套上的概率Ps(dp)、Pb(dp)分別為
(4)
(5)
顆粒物撞擊到閥芯或閥套上的概率反映了撞擊到閥芯或閥套上的顆粒物占流經閥口總的顆粒物的比值。
閥芯和閥套的計算方法一致,以閥芯為例,介紹其沖蝕圓角的計算模型。為計算顆粒物沖蝕引起的沖蝕圓角,需首先計算閥芯表面在單位時間、單位面積上的質量損失,即沖蝕率。根據Edwards的研究,沖蝕率與撞擊到閥口的顆粒物數量成正比,與閥芯的受沖蝕面積成反比,并與撞擊速度、沖擊角度、顆粒物直徑等因素有關[7],可表示為
(6)

閥芯和閥套的材料主要為440C不銹鋼,結合砂粒沖擊碳鋼表面的研究數據[7],顆粒直徑函數C(dp)取經驗值1.8×10-9;沖擊角函數f(α)采用分段函數描述,當沖擊角α為0°、20°、30°、45°和90°時,f(α)分別為0、0.8、1、0.5和0.4;相對速度函數b(v)取2.41[12]。
撞擊到閥芯上的顆粒物質量流率可以用油液中顆粒物的質量流率Rmass乘以碰撞概率表示,故式(6)可寫為
(7)
在t至t+Δt時間段內,m種不同直徑的顆粒物沖蝕造成閥芯的體積損失Vs(Δt)可表示為
(8)
式中,ρs為閥芯材料密度。
Δt時間段內,閥芯的沖蝕圓角半徑從rs(t)變化至rs(t+Δt),閥芯的體積損失Vs(Δt)與圓角的關系有
(9)
式中,Dv為閥芯直徑。
根據式(1)~式(9),可由顆粒物數量、撞擊速度、沖擊角度、顆粒物直徑以及閥口開度等因素得到任意時刻閥芯的沖蝕圓角,同理可得任意時刻閥套的沖蝕圓角。
文獻[11]研究表明,在液壓滑閥中,顆粒的撞擊速度、質量流量、沖擊角度等因素與閥口的壓差、閥口開度、顆粒物直徑等因素有關。其中,顆粒物的撞擊速度主要受顆粒物直徑和壓差的影響,直徑增大,撞擊速度減小。閥口壓差每增大2倍,顆粒物撞擊速度增大1.3倍左右,即顆粒物的撞擊速度為
vdp(Δp)=vp0×1.3lb(Δp/7 000 000)
(10)
式中,Δp為閥口壓差,Pa;vp0為7 MPa壓差下顆粒物的撞擊速度,m/s。
直徑為10 μm、20 μm、60 μm的顆粒物在7 MPa壓差下撞擊速度分別約為41.6 m/s、32.2 m/s、20 m/s,可近似擬合為
(11)
文獻[11]研究發現沖擊角度受顆粒物尺寸影響最為明顯,與閥口開度及閥口壓差的關系不大,因此沖擊角度可看成顆粒物直徑的唯一函數。10 μm、20 μm、60 μm的顆粒物沖擊角度分別約為10°、15°、22°,可近似擬合為
α(dp)=6.712 9ln(dp)+87.453
(12)
由于油液顆粒物的尺寸較小,沖擊角一般不超過20°,在此范圍內,可視沖擊角函數與沖擊角之間為線性關系,即
(13)
而顆粒物質量流率Rmass與通過閥口的流量Q存在如下關系:
(14)
式中,n100 mL為每100 mL油液中直徑為dp的顆粒物數量;ρp為顆粒物密度。
由此建立了單個閥口沖蝕圓角與閥口壓差、流量和閥口開度之間的關系。
四邊滑閥工作時,各個閥口壓差和流量以及閥口開度等與液壓缸負載、活塞運動速度、液壓缸幾何尺寸參數等相關,本節結合圖3所示的伺服閥控對稱缸分析四邊滑閥的沖蝕輪廓變化規律。

圖3 四邊滑閥控對稱缸Fig.3 Four-sided slide valve controlledsymmetrical cylinder
液壓缸的負載FL和兩腔壓力之間的關系為
FL=(p1-p2)A
(15)
式中,p1、p2分別是左右兩腔油液壓力;A為液壓缸油液作用面積。
當活塞桿向右運動時,忽略閥口2和4的流動,閥口1和3節流口寬度、前后壓力和活塞運動速度間的關系分別為
(16)
(17)
式中,Cd為流量系數,取0.63;xk1、xk3分別為閥口1的節流口寬度和閥口3的節流口寬度;ps為供油壓力;vri為活塞向右運動速度;ρ為油液密度。
閥口1的節流口寬度xk1和閥口3的節流口寬度xk3之間存在如下關系:
(18)
(19)
式中,r1、r3分別為閥口1和閥口3的沖蝕圓角半徑;xvr為活塞向右運動時的閥芯位移。
同理,當活塞桿向左運動時,閥口2和4節流口寬度、前后壓力和活塞運動速度間的關系為
(20)
(21)
式中,vle為活塞向左運動速度。
閥口2的節流口寬度xk2和閥口4的節流口寬度xk4的關系如下:
(22)
(23)
式中,r2、r4分別為閥口2和閥口3的沖蝕圓角半徑;xvl為活塞向左運動時的閥芯位移。
在時間T內,活塞桿處于向左運動的時間tle和向右運動的時間tri分別為
(24)
(25)
根據式(15)~式(25)可分別得到指定負載大小、負載流量下兩腔的壓力以及閥芯位移。
閥口輪廓的演化規律計算流程如圖4所示,首先根據負載狀態和閥口初始形貌特征,得到各個閥口的壓差和流速,并得到顆粒物的撞擊速度、角度、碰撞概率等;然后計算各個閥口的沖蝕率,進而得到時間Δt內的沖蝕面積,并得到Δt后的閥口形貌;閥口形貌改變后,為了保證負載需求,閥口開度發生變化,會進一步影響到閥口流動狀態,造成沖蝕演化過程發生變化,以此進行迭代計算,最終獲得任意時刻的四邊滑閥在指定負載條件下沖蝕形貌。

圖4 閥口輪廓的演化規律計算流程Fig.4 The calculation process of the evolution law ofthe valve port profile
為方便對閥控對稱缸動力機構4個閥口的沖蝕過程進行分析,選取基本參數如表1所示的某閥控缸動力機構。假設顆粒物撞擊到閥芯閥套上的概率相同,即閥芯閥套的沖蝕圓角大小相同,統一記為ri(i為閥口編號)。油液中的顆粒物直徑為10 μm,污染物濃度為NAS6級(16 000/100 mL)。假設新閥初始加工圓角半徑為1 μm,零遮蓋。負載大小為pL,活塞桿向左和向右的運動速度分別是vle和vri。

表1 閥控缸動力機構基本參數Tab.1 Basic parameters of valve-controlled cylinder
取液壓缸活塞桿向左向右運動速度一致,負載大小為0.5psA,結合沖蝕輪廓計算方法可得不同時間下閥口各沖蝕圓角,計算時取Δt為1 h。圖5所示為液壓滑閥4個閥口的沖蝕圓角半徑隨服役時間變化的理論結果。液壓滑閥工作2000 h后,閥口1和閥口3的沖蝕圓角半徑約為30 μm,閥口2和閥口4的沖蝕圓角半徑約為84 μm。可以看出,閥口2和4的沖蝕圓角始終比閥口1和3的沖蝕圓角更大。這是由于雖然對稱缸往復運動的速度相等,流經4個閥口的流量相同,但是控制活塞桿縮回時的閥口2和閥口4承受的壓降小于控制活塞桿伸出時的閥口1和閥口3承受的壓降,因此閥口2和閥口4處顆粒物的撞擊速度更大,沖蝕更加嚴重。

圖5 各閥口沖蝕圓角半徑變化曲線Fig.5 Erosion fillet of valve ports
液壓滑閥閥口沖蝕圓角的變化導致其靜態特性發生改變。圖6~圖8分別反映了控制對稱缸的四邊滑閥在服役不同時間后的壓力特性、壓力增益和零偏位移以及空載流量特性。取閥口1打開的方向為正方向(即閥芯向左為正),由于閥芯的位置變化范圍大,故需要考慮閥口形式的變化,取閥口形式轉變的臨界閥芯位置xv0=-(rs+rb+s)。當閥口大于臨界閥芯位置時,閥口為薄壁孔口,流動狀態為湍流;當閥口小于臨界閥芯位置時,閥口為環形縫隙,流動狀態為層流。

圖6 服役不同時間后的壓力特性Fig.6 Pressure characteristic after different service time

圖7 壓力增益和零偏位移隨服役時間變化曲線Fig.7 Pressure gain and null bias after differentservice time

圖8 服役不同時間后的空載流量特性Fig.8 No-load flow characteristic after differentservice time
壓力特性曲線±40%額定供油壓力處兩點連線的斜率為壓力增益kp[13],可以發現,隨著閥口的磨損,滑閥的壓力增益顯著下降,服役2000 h后,壓力增益僅有560 MPa/mm,不足新閥的16%。磨損同時使得滑閥的零位發生了變化,壓力特性曲線右移,零偏位移Δxv不斷加大,液壓滑閥工作2000 h后,零偏約為0.031 mm。零偏同時使流量特性曲線右移,由于閥口沖蝕圓角的影響,使零位附件閥口的面積梯度變大,流量增益略微增大。
圖9所示為服役不同時間后的泄漏情況,可以看出,閥口的沖蝕對滑閥的泄漏也造成了嚴重的影響,對于本例零開口四邊滑閥,雖然其初始泄漏量僅為1.01 L/min,但是閥口的輕微磨損就對泄漏產生嚴重影響,液壓滑閥服役500 h、1000 h、2000 h后的零位泄漏量qe分別達到了6.54 L/min、9.90 L/min、15.44 L/min,從全壽命周期服役的角度來看,滑閥應進行正重疊的設計以減小泄漏量。

圖9 服役不同時間后的泄漏曲線Fig.9 Leakage curves after different service times
四邊滑閥的零偏位移可以通過橋路平衡原理進行計算,滑閥的4個節流邊可等效為電橋中的電阻,構成圖10所示的等效液壓橋路。根據惠斯通電橋平衡原理,液壓橋路平衡時的條件為相對橋臂的液阻值乘積相等,即

圖10 四邊滑閥的等效液壓橋路Fig.10 Equivalent hydraulic bridge of four-sidedslide valve
(26)
(27)
式中,Ri為4個閥口形成的液阻;Δpi、Δqi、Aki分別為4個閥口的壓差、流量和節流口寬度。
由此可得

(28)
式(28)反映了閥芯處于任意位置下四邊滑閥的2個控制壓力與4個閥口面積之間的關系。無論是空載流量零位還是斷載壓力零位下,都有qL=0,p1=p2,液壓橋路的壓力狀態相等,此時等式(28)的等號左側等于1。為使液壓缸滿足平衡條件,需等號右側等于1,令MS1=Ak1/Ak4,表示與閥腔S1的兩個閥口面積比,MS2=Ak2/Ak3,表示與閥腔S2的兩個閥口面積比。當MS1/MS2等于1時,四邊滑閥的等效液壓橋路處于平衡狀態,此時的閥芯位移即為四邊滑閥的零偏位移。因此,四邊滑閥的零偏位移可以通過兩個閥腔的閥口面積比曲線確定,圖11所示為上述控制對稱缸的四邊滑閥工作2000 h后兩個閥口面積比與閥芯位置之間的關系。隨著閥芯位移增大,閥口1和4的面積比MS1不斷增大,閥口2和3的面積比MS2不斷減小,在某一位置下,兩條曲線相交于一點,該點對應的橫坐標即為四邊滑閥的零偏位移。

圖11 閥芯處于不同位置時兩對閥口面積比Fig.11 The area ratio of two pairs of valve ports whenthe spool is in different positions
對某型電液伺服閥進行了200 h耐久性試驗并通過氣動配磨曲線分析其閥芯閥套的沖蝕圓角半徑,圖12所示為耐久性試驗后閥芯上出現的沖蝕磨損痕跡。按照GJB 3370—1998[13]要求,試驗時供油壓力為21 MPa,回油壓力為0.6 MPa,工作介質為航空10號液壓油,油液清潔度等級約為NAS5級,閥芯和閥套材料為440C不銹鋼,伺服閥的輸入信號如表2所示。

圖12 伺服閥閥芯沖蝕磨損痕跡Fig.12 Erosion wear of servo valves pool

表2 耐久性試驗輸入信號Tab.2 Durability test input signal
通過測量試驗前后滑閥副的氣動配磨曲線來定量評估滑閥的沖蝕磨損程度。圖13為滑閥氣動配磨曲線測繪裝置示意圖。伺服閥供油口P和回油口T均接通恒壓氣源;閥腔S1和S2分別通過浮子流量計接通大氣;旋動調節螺釘改變閥芯位置;以千分表讀數為橫坐標,相應的浮子流量計讀數為縱坐標,即可分別得到4個節流口氣體流量-閥芯位移曲線。

圖13 滑閥副氣動配磨曲線測繪裝置示意圖Fig.13 Schematic diagram of spool pneumaticmatching curve measuring device
通過閥口圓角分析方法,得到閥口1至4的正重疊量分別是15.8 μm、15.0 μm、13.7 μm、13.6 μm,200 h耐久性試驗后4個閥口的圓角半徑分別是7.6 μm、5.1 μm、9.6 μm、8.8 μm。試驗前后氣體流量-閥芯位移理論曲線和試驗結果的對比如圖14所示,擬合效果很好,表明沖蝕圓角的分析結果準確可靠。

圖14 試驗前后滑閥副的氣動配磨曲線Fig.14 Match grinding curve before and after thelife test
根據耐久性試驗條件進行滑閥沖蝕圓角半徑的理論計算,取顆粒物平均直徑為10 μm,每百毫升數量8000個,閥口壓降為10.2 MPa。由于耐久性試驗中伺服閥的輸入信號不固定,閥芯位置始終處于動態變化中,為了便于計算,首先分析閥芯位移對平均沖蝕率和撞擊概率的影響,結果如圖15所示。可以看出,由于壓差固定,顆粒物的撞擊速度不變,當閥芯位移較小時,顆粒物撞擊概率不變,平均沖蝕率主要受顆粒物質量流率影響,隨著閥芯位移的增大,顆粒物質量流量增大,因而平均沖蝕率增大;當閥芯位移較大時,隨著閥芯位移的增加,顆粒物質量流率依舊增大,但顆粒物的撞擊概率不斷減小,因此平均沖蝕率基本不變。由于耐久性試驗中伺服閥的開度較大,因此取閥芯位移為100 μm進行簡化計算,得到閥口的沖蝕圓角半徑及其與試驗結果的對比如表3所示。

圖15 閥芯位移對平均沖蝕率和撞擊概率的影響Fig.15 Influence of spool displacement on averageerosion rate and impact probability

表3 沖蝕圓角半徑的試驗結果與理論結果對比Tab.3 Comparison of test and theoreticalresults of erosion fillet
通過表3的對比結果可以發現,理論計算得到4個閥口的沖蝕圓角大小相同,同閥口3和閥口4的沖蝕圓角試驗結果接近,相對誤差不超過10%,但閥口1和2的計算結果比試驗結果略大。進一步分析可知,閥口1和2中油液從供油口流入閥腔,閥口3和4油液從閥腔流出至回油口,流入閥腔的閥口沖蝕圓角比從流出閥腔的閥口沖蝕圓角略小,表明閥口的沖蝕與閥口油液流動方向有關。本文尚未考慮油液流動方向對沖蝕過程的影響,造成了計算結果存在一定偏差,有待進一步研究。此外,由于目前耐久性試驗僅抽樣一臺伺服閥,結果存在一定的隨機性,后續有待通過概率統計與數學抽樣問題理論研究,以及耐久性測試探討統計學規律以及隨機抽樣的影響。總體看來,理論計算結果同試驗結果基本一致。
(1)考慮顆粒物尺寸以及顆粒物撞擊閥芯閥套概率,建立了全周邊液壓滑閥沖蝕圓角計算模型。在閥控缸動力機構中,可通過閥控缸結構尺寸、負載大小、活塞運動速度等參數得到各閥口壓降、流量和閥口開度,并進一步取得撞擊到閥口的顆粒物數量、撞擊速度、沖擊角度,進而定量計算各閥口的沖蝕圓角大小。
(2)滑閥閥口磨損后,將會導致四邊滑閥產生零偏,零偏位移可以通過惠斯通電橋平衡原理求出。滑閥的沖蝕會導致壓力增益顯著降低,泄漏量顯著增大,零位附近的流量增益升高。