基于動網格6DOF方法的液壓閥流場特性分析

郭曉霞,黃家海,,權　龍,王勝國

(1.太原理工大學 新型傳感器與智能控制教育部與山西省重點實驗室,太原 030024;

2.太重榆次液壓工業有限公司 技術中心,山西 榆次 030600)

?

基于動網格6DOF方法的液壓閥流場特性分析

郭曉霞1,黃家海1,2,權龍1,王勝國2

(1.太原理工大學 新型傳感器與智能控制教育部與山西省重點實驗室,太原 030024;

2.太重榆次液壓工業有限公司 技術中心,山西 榆次 030600)

摘要:利用計算流體動力學軟件FLUENT中的6DOF模塊模擬插裝閥閥芯開啟過程中的非定常流場,閥芯運動由其所受流體力確定,當閥芯受力達到平衡時,閥芯停止運動。將閥芯處于穩定時閥芯位移、控制腔壓力、閥芯所受穩態液動力的仿真計算結果與理論公式計算值進行了比較,并對閥芯開啟過程中的瞬態壓力場、速度場進行了分析。結果表明,該仿真計算與理論計算值誤差較小,具有較高的可靠性;與以往研究方法相比,該方法與閥芯實際運動更加接近。本研究能夠較真實地模擬液壓閥開啟過程中的瞬態流場特性,為進一步研究液壓閥性能提供重要理論依據。

關鍵詞:插裝閥;有限元分析;流體動力學軟件-6DOF;液壓流場；理論計算

液壓閥是液壓系統中的重要元件,用于控制系統的液流方向、流量和壓力,其性能優劣在很大程度上決定了液壓系統的使用性能。對于液壓閥而言,傳統的研究方法主要借助于理論推導與實驗研究。近年來,隨著計算機技術和計算流體動力學(CFD,Computational Fluid Dynamics)理論的發展,針對液壓閥內部流場的數值分析和可視化研究逐漸成為液壓閥領域的重要方向之一。

以往的數值模擬多局限于定常研究或直接給定閥芯運動速度[1-3]。如，太原理工大學鄭淑娟通過動網格方法自定義閥芯運動速度,以此來獲得閥芯開啟、關閉過程受到的瞬態液動力、流量系數與通過閥流量、閥開口度之間的定量關系;哈爾濱工業大學劉華坪采用UDF函數給定閥芯速度隨時間的變化規律,來模擬4種閥門動態過程中的流動狀態和閥體受力情況等等。而實際工作中,閥芯的運動規律無法用既定的代數式給出,直接給定閥芯運動速度所獲得的流場存在失真問題,無法獲得閥芯真實的位移曲線與實際的平衡位置。近年來,已有學者開始研究閥芯在流場受力下的運動狀態,多數采用UDF(用戶自定義函數)定義閥芯的受力狀況,以此來實現閥芯的運動。如，解放軍理工大學的何曉暉教授[4]利用動網格技術及UDF函數對液壓球閥二維模型啟閉過程中的液動力進行了研究,得到了閥芯啟閉過程中的流場壓力分布；蘭州理工大學的蘇華山博士[5]基于動網格技術分析了加油泵溢流閥內流體振動及產生噪聲的原因；Dong-A大學的SONG[6]使用CFX軟件的動網格模型分析了溢流閥在不同超壓下的流場特性。以上研究均需要編寫較為復雜且容易出錯的UDF函數,需花費較多時間。

為了探索液壓閥瞬態流場特性,使閥芯在受力不平衡的狀態下自行開啟,受力平衡后處于穩定狀態,筆者采用FLUENT中的動網格6DOF模型來模擬閥芯運動過程,避免了編寫較為復雜的UDF函數。為了驗證此方法的可行性,選用Valvistor閥作為研究對象,將閥芯作為運動實體,閥芯運動由其所受的流體力決定。

1Valvistor閥工作原理及數學建模

Valvistor閥工作原理如圖1所示。當先導閥口關閉時,主閥進油腔油液經反饋節流槽流入并充滿控制腔,上下兩腔壓力相同,由于彈簧力及上下腔面積差作用,主閥處于關閉狀態。當先導閥打開后,控制腔油液經先導閥流至主閥出口,控制腔壓力下降,主閥芯向上移動；當流經反饋節流槽的流量與先導閥的流量相同時,達到穩態,主閥芯停止移動。

閥芯動力平衡方程為:

(1)

1-先導閥；2-閥套；3-閥芯圖1　Valvistor閥工作原理Fig.1　Working principle of valvistor valve

符號含義符號含義Ff液動力Cdx主閥口流量系數Aa閥芯小端面積ρ流體密度pi入口壓力wx主閥口面積梯度Ac閥芯大端面積θ射流角po出口壓力qm通過主閥口的流量pc控制腔壓力¨x閥芯加速度m閥芯質x閥芯速度Bm阻尼系數x閥芯位移k復位彈簧剛度

忽略閥芯與閥套之間的摩擦力及閥芯重力,因其相對于流體力較小。其中,

(2)

(3)

則:

(4)

當流場處于穩態時,通過主閥口的流量為:

(5)

由式(5)可得:

(6)

其中,面積梯度wx為:

(7)

2數值計算方法

2.1幾何與網格劃分

根據實際流域建立液壓閥的流道模型,如圖2所示。因不知先導閥的具體尺寸,無法對先導閥進行建模,因此,根據文中所使用實驗數據,依據

將先導閥閥口簡化為小圓柱閥口。將小圓柱口的尺寸不斷修改,并將多組穩態模擬結果與實驗數據進行對比之后,確定先導閥閥口為直徑1.1 mm,高0.1 mm的圓柱閥口,如圖2中OUTLET-2.

圖2　1/2流體域Fig.2　Half of the fluid domain

對于復雜幾何體而言,生成六面體網格較為困難。為減少劃分網格的時間,采用四面體網格進行網格劃分。為了更好地捕捉液壓閥的內部流場,對壓力、速度梯度較大的局部區域進行網格細化。

Fluent中的動網格技術需要確保拓撲結構不發生變化,運動物體不接觸到邊界,因動網格技術的此項限制,在閥芯與閥套之間有一個0.05 mm的預開口量來確保流體域是一個整體。

網格劃分如圖3所示,其中99.99%的網格質量高于0.35,可認為網格質量符合精度要求。

2.2假設

對Valvistor閥進行解析計算時,對上述模型作如下假設:

1) 假設Valvistor閥為理想錐閥,即閥芯與閥套配合精確無徑向間隙,且閥芯只有平移運動,沒有旋轉運動;

2) 流體為不可壓縮、恒定牛頓流體;

3) 不考慮熱傳遞,流體保持恒溫;

4) 忽略流體重力。

圖3　流域網格Fig.3　Mesh of the flow domain

2.3計算方法

瞬態模擬采用標準k-ε湍流模型;進出口條件為壓力邊界,進口壓力為4.25 MPa,出口壓力為0.7 MPa;激活動網格6DOF模型。因本次模擬采用四面體網格,選定動網格更新方法為彈性光順法與局部網格重構法;壓力-速度耦合方式選用PISO格式;設置時間步長為0.001 ms,仿真時間0.04 s。

2.46DOF模型

因閥芯的運動由作用在其上的流體力、閥芯重力以及其他力共同決定,其運動規律難以用簡單的代數式給出。Fluent所提供的6DOF模型可計算物體在流體作用下的運動學特征(速度、加速度、力等),求出其運動軌跡。該模型將運動部件作為剛體,忽略其變形。

為了計算閥芯在液壓力作用下的運動姿態,在6DOF模型中采用UDF宏DEFINE-SDOF-PROPERTIES定義剛體的質量、動量、轉動慣量等,利用該函數對模型進行簡化,只考慮閥芯沿Z軸的平移運動,約束其在X,Y方向上移動,約束其繞X,Y,Z軸的轉動。

3仿真結果驗證

圖4　閥芯位移曲線Fig.4　Displacement curve of the valve core

圖4為閥芯位移曲線,因閥口處有0.05 mm的預開口量,且沒有考慮先導閥突然開啟的情況,所以位移曲線沒有超調現象。由位移曲線可知,在閥芯開啟的初期階段,閥芯運動速度較大；隨著開啟過程的進行,閥芯運動速度逐漸減小；在0.02 s之前閥芯已經趨于穩定,閥芯位移0.660 5 mm。

圖5為閥芯穩定后,主閥口處的流線圖,由圖可知,射流角約為51°.

圖5　主閥口處流線圖Fig.5　Streamline of the main valve

取Cdx為0.77[7],θ為51°,當閥芯穩定時,通過主閥口的質量流量為1.454 1 kg/s,則由公式(6)計算得到閥芯位移x=0.696 8 mm,理論計算與仿真計算的相對誤差

根據閥芯壁面壓力分布，可以計算閥芯受到的穩態液動力。仿真計算得到的穩態液動力為75.9 N。為驗證本文計算方法,將仿真計算結果與理論計算得到的結果進行比較。由理論計算式(2)推導出的結果為82.3 N，理論計算與仿真計算結果的相對誤差為7.8 %。

由控制腔壓力公式(4)可得:

仿真計算得到的控制腔壓力為2.3 MPa左右,如圖6所示。

由以上計算可知,理論計算與仿真計算結果的誤差較小,可以認為仿真計算結果與理論計算結果相互吻合。因此，該模型具有較高精度和可靠性,能真實模擬液壓閥開啟過程中的流動狀態。

4動態流場分析

受篇幅所限,以下分析只給出t=0.001 s (S=0.12 mm),t=0.01 s(S=0.52 mm),t=0.02 s(S=0.66 mm)時的仿真計算結果。

4.1壓力分析

圖7為閥芯運動過程中不同瞬時的流場壓力云圖。t=0.001 s時刻,閥口前端形成高壓區,過流面積減小處形成低壓區,此時閥芯下端面所受液壓力大于閥芯上端面液壓力,閥芯所受合力方向向上,向趨于打開方向運動;t=0.01 s時刻,油液經主閥口流入閥腔,主閥口處壓力迅速降低,在拐角處出現低壓區;t=0.02 s時刻,閥芯上下端面受力平衡,運動速度近似于0 m/s,位移達到最大值。

從圖中可以看出,隨著閥芯向上運動,閥芯位移逐漸增大,閥芯上下端面壓力逐漸趨于平衡;不同瞬時,閥腔內的壓力分布不斷變化,但閥內最高壓力、最低壓力位置基本相同,壓力分布規律相似;閥口處壓力迅速降低,存在較大的壓力損失;在閥芯錐部上下端與閥座拐角處的位置出現低壓區,在高壓環境下很容易產生氣穴、氣蝕現象。

圖7　不同瞬時壓力分布云圖Fig.7　Pressure distribution contours under different transient times

4.2速度分析

圖8　不同瞬時流場速度云圖Fig.8　Speed contours of flow field under different transient times

圖8給出了Valvistor閥閥芯運動過程中不同瞬時的流場速度云圖。從圖中可以看出,不同瞬時,閥腔內的速度分布不斷變化,但閥內最大速度位置基本相同,即閥口處、上節流口處與控制腔出口處;隨著閥芯向上運動,開口量逐漸增大,相對流速高的區域有所擴大,流速高的區域貼著閥芯側面流動,在閥芯下端面中部流速接近0的區域也在縮小。由不同瞬時的流場速度云圖可看出,Valvistor閥內部流場前后對稱,可簡化為1/2模型,以減少計算資源;結合壓力分布云圖可以看出,最大流速所在區域也是壓降最大的區域,同時流場內的最小負壓值也位于附近。

5結論

1) 利用計算流體動力學(CFD)軟件Fluent 14.5進行了Valvistor閥三維瞬態流場仿真。結合瞬態模擬、動網格技術、動態方程,將閥芯位移、控制腔壓力、閥芯受到穩態液動力的數值計算結果與理論計算值進行比較,結果表明,該液壓閥模型建模正確,該方法具有較高精度和可靠性,能真實模擬液壓閥開啟過程中的流動狀態。因液壓閥的動態特性難以用穩態模擬獲得,此方法對具有復雜結構的液壓閥,在設計階段進行性能評估具有重要意義,之后將利用此模型研究不同壓差下的動態特性。

2) 液體在流動的過程中,過流面積突變處引起流速和壓力的相應變化,在閥芯錐部上下端與閥座拐角處的位置出現低壓區,在高壓環境下很容易產生氣穴、氣蝕現象。

3) 由速度云圖可知,閥內流場前后對稱,可簡化為1/2模型,減少計算資源。參考文獻：

[1]鄭淑娟,權龍,陳青.閥芯運動過程液壓錐閥流場的CFD計算與分析[J].農業機械學報,2007,38(1):168-172.

[2]劉華坪,陳浮,馬波.基于動網格與UDF技術的閥門流場數值模擬[J].汽輪機技術,2008,50(2):106-108.

[3]宋治偉,崔寶玲,尚照輝，等.截止閥啟閉過程內部瞬態流動特性[J].工程熱物理學報,2012,33(6):957-960.

[4]何曉暉,孫宏才,程健生，等.基于動網格的液壓閥閥芯啟閉中的液動力分析[J].解放軍理工大學學報(自然科學),2011,12(5):491-495.

[5]蘇華山,楊國來,張立強，等.加油機溢流閥流體振動噪聲分析與優化[J].振動與沖擊,2013,32(23):130-134.

[6]SONG X G,LIN W,PARK Y C.Transient analysis of a spring-loaded pressure safety valve using computational fluid dynamics(CFD)[J].Journal of Pressure Vessel Technology,2010,132(5):054501.

[7]VAUGHAN N D,JOHNSTON D N,EDGE K A.Numerical simulation of fluid flow in poppet valves[J].Proceedings of the Institution of Mechanical Engineers,Part C:Journal of Mechanical Engineering Science,1992,206(2):119-127.

(編輯：龐富祥)

Flow Field Characteristic Analysis of Valve Based on Dynamic Mesh 6DOF Technique

GUO Xiaoxia1,HUANG Jiahai1,2,QUAN Long1,WANG Shengguo2

(1.KeyLabofAdvancedTransducersandIntelligentControlSystem,MinistryofEducationandShanxiProvince,TaiyuanUniversityofTechnology,Taiyuan,Shanxi030024,China;2.TechnologyCenterofTaiyuanHeaveyMachineryGroupYuciHydraulicIndustry,Yuci,Shanxi030600，China)

Abstract:The numerical simulation and visualization research of internal flow field of hydraulic valve have become one of the important directions in the field of fluid power transmission. But the present numerical analyses are limited to the steady model or constant spool movement velocity,thus impossible to simulate real transient flow field. In this paper, the transient flow field of the cartridge valve in the opening process is simulated by combining computational fluid dynamics software FLUENT with 6DOF model of dynamic mesh technique. The movement of the spool is decided by flow force which exerts on it. The spool stops moving when the forces are balanced.This method makes it easier to realize transient simulation and prevents writing more complex UDF functions. By comparing numerical simulation results about the valve displacement, pressure of the control chamber and the steady flow force with the theoretical formula calculation values,the results show that the numerical results are the same with theoretical value.So the proposed method is feasible to simulate the flow state in the opening process of hydraulic valve. And on this basis, the transient pressure field and velocity field of valve camber in the opening process of the valve are researched,which provides important theoretical basis for further research on hydraulic valve performance.

Key words:dynamic mesh technique;6DOF;fluent;valve;flow field

中圖分類號:TH137

文獻標識碼:A

DOI：10.16355/j.cnki.issn1007-9432tyut.2016.01.003

基金項目：國家自然科學基金資助項目:有源、流量閉環先導級為驅動的電液流量控制理論與方法(511753620);新型先導脈寬調制的電液比例節流與調速復合控制方法及系統集成(51205271);高等學校博士學科點專項科研基金(20121402120002)

作者簡介：郭曉霞(1990-),女,山西陽泉人,碩士生,主要從事液壓閥的CFD仿真研究,(E-mail)372472953@qq.com通訊作者：黃家海,男,副教授,主要從事液壓傳動及理論研究，(E-mail)huangjiahai@tyut.edu.cn

收稿日期：2015-08-29

文章編號:1007-9432(2016)01-0011-05