宋子龍,張曉俊,梁義維,張 杰

(太原理工大學(xué) 機(jī)械與運(yùn)載工程學(xué)院,山西 太原 030024)

0 引 言

滑閥是液壓控制系統(tǒng)中的重要組成部分,而其中的閥芯位移改變流體的流向及閥口的開度,其性能對(duì)整個(gè)液壓系統(tǒng)關(guān)系重大。

液動(dòng)力是影響滑閥性能的重要因素之一[1]。當(dāng)流體流經(jīng)液壓閥的閥口時(shí),其流速和方向會(huì)發(fā)生變化,從而導(dǎo)致流體的動(dòng)量發(fā)生變化,進(jìn)而對(duì)閥芯產(chǎn)生一個(gè)軸向力(即液動(dòng)力),過大的液動(dòng)力會(huì)影響液壓閥的操縱力及穩(wěn)定性,從而對(duì)整個(gè)液壓系統(tǒng)的性能產(chǎn)生影響[2]。

許多學(xué)者采用計(jì)算流體力學(xué)(CFD)對(duì)滑閥內(nèi)部流場(chǎng)的特性及其液動(dòng)力進(jìn)行了研究[3-5]。張曉俊[6]采用對(duì)滑閥閥芯壁面壓力分布的表面積進(jìn)行積分的方式,研究了滑閥的穩(wěn)態(tài)液動(dòng)力。AMIRANTE R等人[9]提出了一種減小比例閥的液動(dòng)力,并保證其所需閥口開啟力不變的方法。AMIRANTE R等人[10]分析了流量對(duì)滑閥液動(dòng)力的影響,并對(duì)此進(jìn)行了試驗(yàn),以驗(yàn)證仿真分析的結(jié)果。廖瑤瑤等人[11]建立了液壓支架用換向閥流道模型,并提出了解決換向閥復(fù)雜流道液動(dòng)力的方法。湯志勇等人[12]提出了一種滑閥的閥芯固定,而其閥套運(yùn)動(dòng)的方法,并將其用于對(duì)滑閥的穩(wěn)態(tài)液動(dòng)力進(jìn)行補(bǔ)償。白瓊[13]對(duì)型槽、V型槽、K型槽和全周開口的水液壓滑閥穩(wěn)態(tài)液動(dòng)力進(jìn)行了計(jì)算,得出了全周開口液動(dòng)力最小,V型槽液動(dòng)力最大的結(jié)論。袁俊等人[14]利用動(dòng)網(wǎng)格分析了電液伺服閥二維模型在不同時(shí)刻的受力情況,為滑閥的設(shè)計(jì)提供了參考。BURECEK A等人[15]通過試驗(yàn)測(cè)量和數(shù)學(xué)模擬的方式,對(duì)比例閥的液動(dòng)力、壓力和閥芯位置的關(guān)系進(jìn)行了研究。

上述研究人員對(duì)滑閥[16]所進(jìn)行的研究多數(shù)著眼于穩(wěn)態(tài)液動(dòng)力。然而在大流量、小開口、閥芯高速運(yùn)動(dòng)時(shí),滑閥的瞬態(tài)液動(dòng)力也會(huì)較大。因此,在上述研究中,研究人員普遍都忽略了閥口氣穴現(xiàn)象對(duì)滑閥液動(dòng)力的影響,導(dǎo)致在滑閥開口度較小時(shí),其壓力值與實(shí)際壓力不符。

以非全周開口滑閥為研究對(duì)象,筆者采用兩相流模型對(duì)其內(nèi)部流場(chǎng)進(jìn)行分析,通過仿真計(jì)算得出滑閥穩(wěn)態(tài)液動(dòng)力和瞬態(tài)液動(dòng)力,進(jìn)行具體分析,并通過液動(dòng)力對(duì)閥套結(jié)構(gòu)進(jìn)行優(yōu)化,為液壓閥的設(shè)計(jì)提供理論依據(jù)。

1 網(wǎng)格模型及其計(jì)算

1.1 幾何模型

非全周開口內(nèi)流式滑閥結(jié)構(gòu)示意圖如圖1所示。

圖1 非全周開口內(nèi)流式滑閥結(jié)構(gòu)示意圖

圖1中,滑閥由閥體、閥套和閥芯組成,閥套入口處加工4個(gè)方形孔進(jìn)行節(jié)流,出口處為4個(gè)圓形孔,不產(chǎn)生節(jié)流作用。

x為閥芯運(yùn)動(dòng)過程中所研究的閥的開口值,其主要尺寸為:a=18 mm,b=10 mm,c=d=6.5 mm,e=3 mm,f=6.5 mm,D=12 mm,m=18 mm。

1.2 網(wǎng)格劃分

由于滑閥內(nèi)部流道比較復(fù)雜,筆者采用SolidWorks軟件和Geometry軟件建立滑閥流體域三維模型,如圖2所示。

圖2 滑閥流體域三維模型

由于滑閥三維模型的結(jié)構(gòu)具有對(duì)稱性,為了減小計(jì)算量和提高計(jì)算精度,可取其半剖模型進(jìn)行研究。

筆者使用前處理軟件ICEM CFD進(jìn)行網(wǎng)格劃分。為了便于計(jì)算收斂,大部分區(qū)域采用六面體網(wǎng)格,部分區(qū)域采用四面體網(wǎng)格;考慮到閥口附近壓力和速度梯度較大,故進(jìn)行細(xì)化處理。

同時(shí),筆者對(duì)不同數(shù)量的網(wǎng)格進(jìn)行無關(guān)化檢驗(yàn),即分別選取38 128、62 452和106 575的網(wǎng)格數(shù)量進(jìn)行檢驗(yàn),結(jié)果發(fā)現(xiàn)其對(duì)計(jì)算的結(jié)果沒有影響。

流體域半剖網(wǎng)格模型如圖3所示。

圖3 流體域半剖網(wǎng)格模型

1.3 計(jì)算條件

筆者采用速度入口和壓力出口進(jìn)行仿真計(jì)算,并且通過給定進(jìn)口速度來控制流量的大小。半剖截面設(shè)為對(duì)稱面,其余與流體接觸壁面為靜止壁面;流體設(shè)為牛頓流體,不可壓縮,流動(dòng)狀態(tài)為紊流,采用標(biāo)準(zhǔn)k-ε模型;流體介質(zhì)選為液壓油,密度為850 kg/m3,運(yùn)動(dòng)黏度為0.02 kg/m·s,出口壓力為0.12 MPa。計(jì)算模型和計(jì)算條件與實(shí)際工況相對(duì)應(yīng)。仿真過程中選擇適當(dāng)?shù)那匪沙谝蜃?保證計(jì)算收斂,從而得到準(zhǔn)確合理的解。

考慮到實(shí)際情況中的氣穴現(xiàn)象,筆者在Fluent中的多相流模型里選擇氣穴模型,采用氣穴模型進(jìn)行計(jì)算可以模擬流場(chǎng)里的氣泡在壓降下形成氣穴。研究中運(yùn)用Fluent中的動(dòng)網(wǎng)格技術(shù)來模擬閥芯移動(dòng)的過程,采用動(dòng)態(tài)層技術(shù)計(jì)算,運(yùn)動(dòng)區(qū)域設(shè)置為剛體運(yùn)動(dòng)。

當(dāng)閥芯的運(yùn)動(dòng)速度為0.5 m/s時(shí),閥芯運(yùn)動(dòng)不同位置的網(wǎng)格示意圖如圖4所示。

圖4 閥芯運(yùn)動(dòng)不同位置網(wǎng)格示意圖

2 仿真計(jì)算及分析

2.1 壓力場(chǎng)和速度場(chǎng)

閥芯運(yùn)動(dòng)使閥口打開時(shí),不同流量和閥口開度的滑閥穩(wěn)態(tài)流場(chǎng)壓力云圖和速度云圖,如圖5所示。

當(dāng)閥芯的運(yùn)動(dòng)速度為0.1 m/s時(shí),滑閥的瞬態(tài)流場(chǎng)壓力云圖和速度云圖如圖6所示。

圖5 滑閥穩(wěn)態(tài)流場(chǎng)壓力云圖和速度云圖

圖6 滑閥瞬態(tài)流場(chǎng)壓力云圖和速度云圖

由圖5和圖6可知:

油液流經(jīng)閥口時(shí)會(huì)產(chǎn)生壓降,這種現(xiàn)象在低壓區(qū)域也比較明顯;當(dāng)油液剛進(jìn)入閥內(nèi)時(shí),其速度基本穩(wěn)定;當(dāng)油液經(jīng)過閥口時(shí),由于過流面積突然減小,其流速增大,由伯努利方程可知,閥口處的壓力減小;相同開口和流量情況下,瞬態(tài)流場(chǎng)的壓降和射流速度明顯小于穩(wěn)態(tài)流場(chǎng);

隨著閥口開度的增大,閥口處的低壓區(qū)域明顯減小,閥口處前后壓差變小,壓力損失減少,最大壓力減小,最小壓力升高,閥內(nèi)整體壓力分布改變,因此閥口開度可以改變閥內(nèi)壓力分布情況。隨著開口度的增大,油液流經(jīng)閥口的通流面積增大,閥口的節(jié)流作用減弱,壓力升高,閥口處的射流速度減小,入口射流角減小;

當(dāng)閥口的開度一定時(shí),流量越大,滑閥流場(chǎng)的最大壓力和閥口處的壓力差都會(huì)變大,壓力損失也會(huì)增加。因?yàn)閴毫ψ兇?所以液流經(jīng)過閥口時(shí)的射流速度明顯增大,同時(shí)射流角變小。

2.2 兩相流仿真分析

由于油液流經(jīng)閥口時(shí)的射流速度較大,會(huì)產(chǎn)生油液和壁面分離的現(xiàn)象,主流的油液會(huì)帶動(dòng)壁面附近的油液以一定的能量流動(dòng),然后回流。主流的油液動(dòng)能很大,回流的油液會(huì)因回流不充分使流場(chǎng)產(chǎn)生低壓區(qū)域,區(qū)域內(nèi)的絕對(duì)壓力比油液和空氣的分離壓力小,因此,油液中有氣泡產(chǎn)生,也即出現(xiàn)氣穴現(xiàn)象。

流量為60 L/min,閥口開度為0.5 mm時(shí)的氣相體積分?jǐn)?shù)云圖如圖7所示。

圖7 氣相體積分?jǐn)?shù)云圖

2.3 液動(dòng)力

液動(dòng)力是因?yàn)橛鸵毫鹘?jīng)閥口時(shí)流速和方向發(fā)生變化,引起油液的動(dòng)量發(fā)生變化,從而使閥芯受到的一種附加作用力。通常,液動(dòng)力可以采用動(dòng)量定理來計(jì)算。

根據(jù)動(dòng)量定理,可得滑閥的穩(wěn)態(tài)液動(dòng)力為:

FS=ρqv1cosθ1-ρqv2cosθ2

(1)

通常認(rèn)為流體出口射流角θ2為90°,因此,此處的液動(dòng)力公式可寫為:

FS=ρqv1cosθ1=2CqCvSΔpcosθ1

(2)

式中:Cq—流量系數(shù);Cv—流速系數(shù);Δp—進(jìn)出口壓差;S—節(jié)流口過流面積;θ1—入口射流角,最理想為69°。

滑閥的瞬態(tài)液動(dòng)力動(dòng)量公式為:

(3)

采用CFD仿真計(jì)算得到的穩(wěn)態(tài)液動(dòng)力,實(shí)際就是在閥芯靜止時(shí),對(duì)閥芯兩壁面壓力的面積進(jìn)行積分得到的差值。而瞬態(tài)液動(dòng)力則是在閥芯運(yùn)動(dòng)狀態(tài)下,對(duì)閥芯兩壁面壓力的面積進(jìn)行積分得到的差值。

不同流量下,隨著閥口開度變化的穩(wěn)態(tài)液動(dòng)力的仿真曲線如圖8所示。

圖8 穩(wěn)態(tài)液動(dòng)力仿真曲線正值—液動(dòng)力指向閥口關(guān)閉方向;負(fù)值—液動(dòng)力指向閥口打開方向

由圖8可知:在閥口開度一定時(shí),穩(wěn)態(tài)液動(dòng)力隨著流量的增大而變大。因?yàn)榱髁吭龃髸?huì)使進(jìn)口處液流流速和進(jìn)出口壓差增大,所以由式(2)可知,其穩(wěn)態(tài)液動(dòng)力也隨之增大;

在流量一定時(shí),閥口開度增大會(huì)導(dǎo)致進(jìn)出口壓差減小,進(jìn)而穩(wěn)態(tài)液動(dòng)力也減小。閥口開度在0～1 mm時(shí),穩(wěn)態(tài)液動(dòng)力較大,減小的幅度也較為明顯。因此,在閥口開度微小的時(shí)候,穩(wěn)態(tài)液動(dòng)力的變化比較大;在閥口開度為1 mm～2.5 mm時(shí),穩(wěn)態(tài)液動(dòng)力較小,減小幅度比較緩慢;

在閥口開度不斷增大的過程中:當(dāng)開口在0～2 mm內(nèi)時(shí),穩(wěn)態(tài)液動(dòng)力雖然一直在減小,但始終為正值,即液動(dòng)力方向?yàn)橼呌陂y口關(guān)閉方向;當(dāng)閥口開度在2 mm～2.5 mm時(shí),穩(wěn)態(tài)液動(dòng)力繼續(xù)減小,直到其值變?yōu)樨?fù)值,即液動(dòng)力方向?yàn)橼呌陂y口打開方向。這是因?yàn)殡S著閥口開度的增大,閥口處的前后壓差減小,入口處的射流速度和射流角減小,當(dāng)其開口達(dá)到一定程度時(shí),射流角接近90°,因此,由式(1)可知,液動(dòng)力會(huì)出現(xiàn)負(fù)值。

隨著閥口開度的變化,不同的流量和閥芯運(yùn)動(dòng)速度的瞬態(tài)液動(dòng)力仿真曲線,如圖9所示。

圖9 瞬態(tài)液動(dòng)力仿真曲線

從圖9中可知:

影響滑閥瞬態(tài)液動(dòng)力的主要因素有閥口開度、閥芯速度以及流量;和穩(wěn)態(tài)液動(dòng)力一樣,瞬態(tài)液動(dòng)力與滑閥的閥口開度成反比(即隨著開口度的增大,瞬態(tài)液動(dòng)力減小);在閥口開度和閥芯移動(dòng)速度不變時(shí),瞬態(tài)液動(dòng)力隨著流量的增大而變大;當(dāng)流量一定時(shí),閥芯運(yùn)動(dòng)速度越快,瞬態(tài)液動(dòng)力就越大。

對(duì)比瞬態(tài)液動(dòng)力和穩(wěn)態(tài)液動(dòng)力的仿真曲線圖可知:

在相同條件下,滑閥穩(wěn)態(tài)液動(dòng)力比瞬態(tài)液動(dòng)力大很多;由式(1,2)可知,穩(wěn)態(tài)液動(dòng)力和瞬態(tài)液動(dòng)力都與閥口處的壓差有關(guān);

由圖(8,9)可看出:瞬態(tài)流場(chǎng)閥口處的壓差明顯小于穩(wěn)態(tài)流場(chǎng),因此瞬態(tài)液動(dòng)力比穩(wěn)態(tài)液動(dòng)力小。但是,在閥芯運(yùn)動(dòng)速度很快、流量較大,并且閥口開度很小時(shí),其瞬態(tài)液動(dòng)力也會(huì)比較大,不可以忽略不計(jì)。因此,在對(duì)滑閥進(jìn)行研究和設(shè)計(jì)時(shí),瞬態(tài)液動(dòng)力也要加以重視。

3 閥套結(jié)構(gòu)優(yōu)化及分析

液動(dòng)力是影響液壓閥性能的因素之一,它不僅影響閥的驅(qū)動(dòng)力,也會(huì)使閥產(chǎn)生自激振動(dòng),降低閥的操縱穩(wěn)定性,對(duì)整個(gè)液壓系統(tǒng)的穩(wěn)定性造成影響,甚至?xí)档烷y的使用壽命。因此,對(duì)滑閥的液動(dòng)力進(jìn)行優(yōu)化補(bǔ)償研究具有重要的意義。

筆者提出一種針對(duì)滑閥閥套結(jié)構(gòu)的液動(dòng)力優(yōu)化方案,即將閥套進(jìn)油孔處的4個(gè)方孔設(shè)計(jì)為斜孔,油液進(jìn)入閥套的方向?yàn)閮A向于閥口打開的方向,傾斜角為θ。

θ=20°時(shí)的閥套半剖模型,如圖10所示。

圖10 θ=20°閥套半剖模型

當(dāng)傾斜角為θ=10°、θ=20°和θ=30°時(shí),筆者分別建立滑閥的三維模型,并在同等條件(優(yōu)化前)下,對(duì)其進(jìn)行仿真計(jì)算。

θ=20°和θ=0°(優(yōu)化前),閥芯開啟時(shí)的液動(dòng)力仿真曲線,如圖11所示。

圖11 優(yōu)化前后液動(dòng)力仿真曲線

由圖11可知:

(1)將閥套進(jìn)油孔設(shè)計(jì)為斜孔后,可以有效降低滑閥的液動(dòng)力,且流量越大,優(yōu)化效果越明顯;

(2)優(yōu)化前后液動(dòng)力仿真曲線的斜率很相似,都是在閥口開度達(dá)到1 mm之前,液動(dòng)力減小較快,在1 mm之后變得緩慢。

不同之處在于,優(yōu)化前的液動(dòng)力隨著開口增大而減小,方向?yàn)橼呌陂y口關(guān)閉方向,在開口2.2 mm時(shí)接近零點(diǎn),液動(dòng)力變?yōu)樨?fù)值,方向趨于閥口打開,且隨著開度的增大而變大;在開度達(dá)到1 mm時(shí),優(yōu)化后的穩(wěn)態(tài)液動(dòng)力接近于0;在閥口開度0.8 mm時(shí),瞬態(tài)液動(dòng)力接近于0。這是因?yàn)閮?yōu)化后閥套的傾斜方向?yàn)殚y口打開方向,因此進(jìn)油方向也傾向于閥口打開方向,相較于優(yōu)化前,壁面1受到的力更大,優(yōu)化后的瞬態(tài)液動(dòng)力和穩(wěn)態(tài)液動(dòng)力絕對(duì)值明顯減小。

當(dāng)流量為60 L/min,閥芯以0.1 m/s的速度打開閥口時(shí),不同傾斜角的滑閥液動(dòng)力仿真曲線如圖12所示。

圖12 不同傾斜角的液動(dòng)力仿真曲線

由圖12可知,不同傾斜角對(duì)滑閥液動(dòng)力的優(yōu)化效果各有差異[17，18]:

當(dāng)傾斜角為θ=10°時(shí)和θ=0°時(shí),在閥口開度小于0.5 mm時(shí),優(yōu)化前后滑閥液動(dòng)力非常接近。這是因?yàn)樵谛￠_度時(shí),開度變化對(duì)壁面1壓力分布的影響不是很明顯;

θ=20°時(shí),滑閥瞬態(tài)液動(dòng)力和穩(wěn)態(tài)液動(dòng)力的絕對(duì)值最小;

θ=30°時(shí),在閥口開度較小時(shí),滑閥液動(dòng)力很小;且在開口為0.5 mm左右時(shí),滑閥穩(wěn)態(tài)液動(dòng)力就改變了方向;開口在0.3 mm左右時(shí),隨著閥口開度的增大,滑閥瞬態(tài)液動(dòng)力幾乎呈線性增大,在開口為2.5 mm時(shí),其值接近10 N。

當(dāng)流量為60 L/min,閥芯以0.1 m/s的速度打開閥口時(shí),隨著傾斜角的變化,不同閥口開度的的液動(dòng)力仿真曲線,如圖13所示。

圖13 不同閥口開度的液動(dòng)力仿真曲線

由圖13可知:

不同閥口開度下,隨著傾斜角的增大,滑閥穩(wěn)態(tài)液動(dòng)力和瞬態(tài)液動(dòng)力的絕對(duì)值都先減小而后再增大;當(dāng)傾斜角在15°～20°時(shí),在各閥口開度下,穩(wěn)態(tài)液動(dòng)力和瞬態(tài)液動(dòng)力的絕對(duì)值最小,表明其優(yōu)化效果最好(其中,當(dāng)穩(wěn)態(tài)液動(dòng)力最大值為7 N時(shí),滑閥的瞬態(tài)液動(dòng)力最大值為4 N);

優(yōu)化后的滑閥液動(dòng)力降低了很多,但隨著閥口開度的增大,液動(dòng)力方向變?yōu)橼呌陂y口打開方向的時(shí)間比優(yōu)化前早,且改變方向的液動(dòng)力也比優(yōu)化前大。

產(chǎn)生以上情況的原因在于:(1)優(yōu)化后的滑閥液動(dòng)力最大值小于優(yōu)化前;(2)當(dāng)閥芯運(yùn)動(dòng)打開閥口時(shí),趨于閥口打開方向的液動(dòng)力會(huì)有助于閥芯的開啟,減小閥芯運(yùn)動(dòng)的阻力(但是液動(dòng)力也不可過大,否則會(huì)產(chǎn)生自激振動(dòng),影響閥的操縱穩(wěn)定性)。

4 結(jié)束語(yǔ)

筆者利用Fluent軟件對(duì)非全周開口內(nèi)流式滑閥內(nèi)部的瞬態(tài)流場(chǎng)和穩(wěn)態(tài)流場(chǎng)進(jìn)行了研究,計(jì)算了閥芯運(yùn)動(dòng)狀態(tài)下的瞬態(tài)流場(chǎng),以及閥芯靜止時(shí)的穩(wěn)態(tài)流場(chǎng),并對(duì)其液動(dòng)力進(jìn)行了分析,提出了一種閥套結(jié)構(gòu)的優(yōu)化方法。

研究結(jié)果表明:

(1)液動(dòng)力與閥口開度和流量大小有關(guān),瞬態(tài)液動(dòng)力還受到閥芯速度的影響;相比于穩(wěn)態(tài)液動(dòng)力,瞬態(tài)液動(dòng)力較小;閥芯高速運(yùn)動(dòng)、大流量、小開口時(shí),滑閥瞬態(tài)液動(dòng)力也較大,所以其影響不可忽略;

(2)當(dāng)閥口開度增大到一定程度時(shí),液動(dòng)力會(huì)出現(xiàn)負(fù)值,即傾向于閥口開啟的方向;

(3)在對(duì)滑閥閥套進(jìn)行優(yōu)化時(shí),將閥套4個(gè)進(jìn)油孔的方孔設(shè)計(jì)為斜孔,對(duì)不同傾斜角的液動(dòng)力進(jìn)行對(duì)比可知,當(dāng)傾斜角在15°～20°時(shí),其液動(dòng)力的優(yōu)化效果最好。

因?yàn)榛y的閥腔不對(duì)稱,所以閥芯在運(yùn)動(dòng)過程中,閥桿會(huì)因受力不均而產(chǎn)生徑向不平衡的卡緊力。因此,在后續(xù)的工作中,筆者將對(duì)因閥桿受力不均產(chǎn)生的卡緊力進(jìn)行深入研究。