， ，

(蘭州理工大學(xué) 能源與動力工程學(xué)院， 甘肅 蘭州 730050)

引言

液壓閥是液壓系統(tǒng)中非常重要的元件，主要通過控制流體的壓力、流量和流動方向來滿足工作要求，使各類執(zhí)行元件實(shí)現(xiàn)不同的動作[1]。液壓控制閥的內(nèi)部結(jié)構(gòu)比較復(fù)雜，主要由閥體、閥芯、操縱控制機(jī)構(gòu)等主要零部件組成?；y類的閥芯是圓柱形，通過閥芯在閥體孔內(nèi)的滑動來改變液流通路開口的大小，以實(shí)現(xiàn)對液流壓力、流量及方向的控制。非全周開口滑閥具有水力半徑大，抗阻塞的特點(diǎn)，其流量范圍大，易得到較小的穩(wěn)定流量，在液壓比例閥和伺服閥中得到了廣泛應(yīng)用[2]，節(jié)流槽滑閥的特性分析對液壓閥的性能提升起著很重要的作用。

近幾年隨著計(jì)算機(jī)科學(xué)技術(shù)的不斷發(fā)展以及計(jì)算流體力學(xué)理論的豐富。人們借助CFD技術(shù)對液壓閥復(fù)雜內(nèi)部流場進(jìn)行數(shù)值模擬和可視化分析，成為液壓領(lǐng)域的熱點(diǎn)。在液壓技術(shù)中，流體仿真主要研究方向有：氣穴現(xiàn)象分析，流場速度壓力流量分析和液動力分析等。文獻(xiàn)[3-6]利用Fluent軟件分別對伺服閥、電磁閥、滑閥式換向閥以及液壓比例閥的三維扇形區(qū)域內(nèi)部流場進(jìn)行模擬仿真，得到壓力與流量，穩(wěn)態(tài)液動力與閥口開度之間的關(guān)系。文獻(xiàn)[7]采用CFX軟件對三維溢流閥模型的內(nèi)部流場進(jìn)行模擬仿真，并結(jié)合流體動力學(xué)方程對溢流閥的瞬態(tài)流場特性和動態(tài)特性進(jìn)行了研究。目前對液壓閥內(nèi)部流場的仿真分析主要是針對全周開口滑閥即傳統(tǒng)滑閥,關(guān)于節(jié)流槽滑閥又稱為非全周開口滑閥的研究比較少。本研究以非全周開口液壓閥滑閥為研究對象，針對滑閥的內(nèi)部流場進(jìn)行三維穩(wěn)態(tài)仿真模擬。采用有限元分析軟件ANSYS的前處理模塊GAMBIT，分析計(jì)算模塊Fluent以及后處理模塊TECPLOT，得到了滑閥內(nèi)部流場的速度壓力、流量特性以及流量系數(shù)隨閥口開度的變化規(guī)律。

1 計(jì)算流體力學(xué)基礎(chǔ)

任何物體的運(yùn)動都要遵循一些運(yùn)動規(guī)律，流體運(yùn)動也不例外，因此滑閥內(nèi)部流體的流動也要遵循質(zhì)量守恒定律、動量守恒定律、能量守恒定律。

1.1 流體動力學(xué)的連續(xù)性方程

連續(xù)性方程就是我們所說的質(zhì)量守恒方程，質(zhì)量守恒定律是流體力學(xué)中最基本的定律。在一定的時(shí)間內(nèi)，某一特定平面上流入的流體質(zhì)量與流出的相等，得到流體流動時(shí)的質(zhì)量守恒微分方程[8]：

(1)

式中：ux、uy、uz分別為x、y、z，3個(gè)方向上的速度分量，m/s；t為時(shí)間，s；ρ為流體密度，kg/m3。

1.2 流體動力學(xué)的動量方程

(2)

式中：Su,Sv和Sw是動量守恒方程的廣義源項(xiàng)，Su=Fu+sx，Sv=Fv+sy，Sw=Fw+sz，其中，sx，sy和sz表達(dá)式：

1.3 湍流模型

流體的流動狀態(tài)可以分為層流和湍流，而湍流在自然界中很常見，大部分工程中的流動問題都是湍流模型[9]。

湍流模型主要有隨機(jī)性、有渦性、擴(kuò)散性和耗散性等特征。當(dāng)流體通過液壓閥時(shí)，在液壓閥內(nèi)部由于節(jié)流口的作用，使流體流動的時(shí)候雷諾數(shù)很大，因此液壓閥內(nèi)部流場的流動狀態(tài)是湍流。

2 液壓閥閥口流場的分析

2.1 幾何模型的建立

液壓閥模型的正確建立對后續(xù)的仿真計(jì)算有著非常重要的意義。為了使模型很好的反應(yīng)出流場的流動情況，在建模的時(shí)候?qū)δＰ鸵M(jìn)行適當(dāng)?shù)暮喕?。采用三維設(shè)計(jì)軟件對滑閥閥口開度分別為：10 mm、16 mm、22 mm時(shí)進(jìn)行建模，主閥芯直徑為120 mm,流道油液進(jìn)口直徑為140 mm，為了能夠清楚的看出在不同閥口開度下流場的運(yùn)動情況，取每個(gè)閥口開度下的1/4來進(jìn)行流場仿真，其中開度為16 mm時(shí)的三維模型如圖1所示。

圖1 開口為16 mm時(shí)三維模型

2.2 采用GAMBIT進(jìn)行網(wǎng)格劃分如圖2所示

圖2 開口為16 mm時(shí)網(wǎng)格劃分

2.3 邊界條件的設(shè)置及相關(guān)的流體參數(shù)

(1) 本研究中所使用的液壓滑閥是比較理想的，裝配精確并且加工良好，沒有什么縫隙和瑕疵；

(2) 所使用的流體不可壓縮，是理想流體；

(3) 流體與壁面的接觸邊界是無滑移壁面；

(4) 采用的邊界條件為壓力入口和壓力出口。操作環(huán)境在標(biāo)準(zhǔn)大氣壓下，其中取進(jìn)口壓力為2.1 MPa， 液壓油密度為889 kg/m3，油液動力黏度為0.03556 kg/(m·s)。

3 液壓閥內(nèi)部流場的計(jì)算結(jié)果

3.1 壓力分布云圖

在保證其它條件相同的情況下分別對閥口開度為10 mm、16 mm、22 mm時(shí)進(jìn)行流場仿真，壓力分布云圖如圖3～圖5所示。

圖3 閥口開度10 mm,壓差2.1 MPa時(shí)壓力分布圖(Pa)

圖4 閥口開度16 mm,壓差2.1 MPa時(shí)壓力分布圖(Pa)

圖5 閥口開度22 mm,壓差2.1 MPa時(shí)壓力分布圖(Pa)

從壓力分布圖以及壓力值可以看出，當(dāng)流體流經(jīng)節(jié)流槽閥口時(shí)，壓力下降比較快。隨著閥口開度的增大，流場的最低壓力減小，閥口開度為16 mm和 22 mm 時(shí)流場的最低壓力相差不大，但是與10 mm相比，最低壓力還是有所降低。壓力的變化范圍也隨著閥口開度的增大而增大。由于流體流經(jīng)節(jié)流槽時(shí)流速很高，根據(jù)動量守恒定律，該處的壓力會很低，甚至有負(fù)壓會產(chǎn)生，可能會產(chǎn)生氣穴現(xiàn)象，當(dāng)液壓系統(tǒng)中有氣穴現(xiàn)象產(chǎn)生時(shí)，會造成液壓系統(tǒng)不穩(wěn)定。

3.2 速度分布云圖

以閥口開度為10 mm、16 mm時(shí)的速度分布圖為例，如圖6，圖7所示。

圖6 閥口開度10 mm,壓差2.1 MPa時(shí)速度分布圖(m/s)

圖7 閥口開度16mm,壓力差2.1MPa時(shí)速度分布圖(m/s)

從上圖可以看出，在相同的壓差條件下，閥口開度為10 mm和16 mm時(shí)，滑閥流場處的最大速度降低,并且隨著閥口的增大，閥口處的最大速度也隨之降低。流體在進(jìn)口處的速度基本不變，但在經(jīng)過閥口時(shí)，速度發(fā)生急劇變化，經(jīng)過閥腔時(shí)流速產(chǎn)生擴(kuò)散。在不同壓差條件下，對開度為16 mm時(shí)的該滑閥進(jìn)行了內(nèi)部流場仿真，得到：在相同的開口條件下，隨著壓差的增大，閥口的流速會隨之增大。

4 液壓閥壓力流量特性分析

4.1 不同閥口開度下的流量特性分析

滑閥的壓力流量特性是指流經(jīng)滑閥的流量與閥口開度以及閥口壓差三者之間的關(guān)系。由流量公式得到流量系數(shù)Cd的公式：

(6)

式中：qν為流量，m3/s；A為閥口過流面積，m2；ρ為液壓油密度，kg/m3；Δp為閥口進(jìn)出口壓力差，MPa。

根據(jù)仿真結(jié)果得到該滑閥在不同開度下的壓力-流量特性曲線圖，如圖8所示。

圖8 不同閥口開度下的壓力-流量特性曲線

從圖中可以看出通過閥口的流量與開度呈線性分布的關(guān)系，表明在滑閥開啟的過程中流量與閥口開度成正比。

通過式(6)我們得到不同閥口開度下的壓力-流量系數(shù)曲線，如圖9所示。

從圖9可以看出：節(jié)流槽閥口的流量系數(shù)與閥口的開度大小密切相關(guān)的。隨著閥口開度的減小，流量系數(shù)逐漸增大。

圖9 不同閥口開度下的壓力-流量系數(shù)曲線

4.2 與半圓形槽滑閥的比較

半圓形槽示意圖如圖10所示。

圖10 不同節(jié)流槽形狀下的壓力-流量系數(shù)曲線圖

在相同條件下，對閥口開度為16 mm的半圓形槽滑閥進(jìn)行模擬仿真，將得到的數(shù)據(jù)代入式(6)中得到流量系數(shù)值，并且與之前的矩形節(jié)流槽比較，結(jié)果如圖11所示。從圖中看出：滑閥流量系數(shù)與節(jié)流槽的形狀有關(guān)；在閥口開度相同的條件下，半圓形槽滑閥與矩形節(jié)流槽滑閥相比較，流量系數(shù)變化較平緩，并且隨著進(jìn)出口壓差的增大都有緩慢下降趨于穩(wěn)定的趨勢。

圖11 半圓形槽示意圖

5 結(jié)論

本研究對液壓閥滑閥進(jìn)行了內(nèi)部流場的仿真分析，得出如下結(jié)論：

(1) 在進(jìn)出口壓差相同的條件下，流體流經(jīng)節(jié)流槽閥口時(shí)，壓力下降較快；隨著閥口開度從10 mm增加到6 mm時(shí)，流場的最低壓力有所增加，閥口開度為16 mm 和22 mm時(shí)流場的最低壓力相差不大，但是與10 mm相比，最低壓力有所下降.壓力的變化范圍也隨著閥口開度的增大而增大。

(2) 進(jìn)出口壓差相同的情況下，隨著閥口開度的增大，閥口的最大速度隨之降低；在相同的開口條件下，隨著進(jìn)出口壓差的增大，閥口的流速會隨之增大；流體在進(jìn)口處的速度基本相同，但在經(jīng)過閥口時(shí)，速度發(fā)生急劇變化，經(jīng)過閥腔時(shí)流速產(chǎn)生擴(kuò)散。

(3) 滑閥節(jié)流槽閥口的流量系數(shù)與閥口的開度大小，進(jìn)出口壓力差以及節(jié)流槽形狀緊密相關(guān)。在相同壓差條件下，隨著閥口開度的減小，流量系數(shù)逐漸增大；在相同的閥口開度下，隨著進(jìn)出口壓差的增大，流量系數(shù)在逐漸減小，半圓形槽滑閥與矩形節(jié)流槽滑閥相比較，流量系數(shù)變化較平緩，并且，隨著進(jìn)出口壓差的增大都有緩慢下降趨于穩(wěn)定的趨勢。

參考文獻(xiàn)：

[1] 李壯云.液壓元件與系統(tǒng)(第三版)[M]. 北京：機(jī)械工業(yè)出版社，2011.

[2] 冀宏，張繼環(huán)，王東升，等.滑閥矩形節(jié)流槽閥口的流量系數(shù)[J].蘭州理工大學(xué)學(xué)報(bào)，2010，36(3)：1-2.

[3] 呂庭英，黃效國，何康寧.基于Fluent的液壓伺服閥液動力研究[J].機(jī)床與液壓，2011,39(13)：131-132.

[4] GEE Soo Lee， HYUNG Jin Sung， HYUN Chul Kim, et al. Flow Force Analysis of a Variable Force Solenoid Valve for Automatic Transmissions[J]. Journal of Fluids Engineering，2010，132(3)：782-793.

[5] 周元春，楊曙東，王金林，等. 基于Fluent的大通徑滑閥壓力流量特性研究[J].機(jī)床與液壓，2011，39(19)：103-105.

[6] Jose R.Valdes，Mario J. Miana，Jose L. Nunez, et al. Reduced Order Model for Estimation of Fluid and Flow Forces in Hydraulic Proportional Valves [J]. Energy Conversion and Management，2008，49(8)：1517-1529.

[7] XUE Guan Song，LIN Wang， YONG Chul Park. Transient Analysis of a Spring-Loaded Pressure Safety Valve Using Computational Fluid Dynamics(CFD)[J]. Journal of Pressure Vessel Technology，2010，132(5)：156-164.

[8] 金曉宏，李遠(yuǎn)慧.流體力學(xué)[M].北京：中國電力出版社，2011.

[9] 潘永成，王勇，等.基于CFD的調(diào)節(jié)閥內(nèi)流場流動特性研究[J].機(jī)床與液壓，2011，(1)：5-7.