油罐自力式液位控制閥流場分析

王 銘,張世富,張起欣,張冬梅,陳 暢

(中國人民解放軍后勤工程學院 a.研究生大隊； b.國家救災應急裝備工程技術研究中心；c.軍事供油工程系， 重慶 401331

?

油罐自力式液位控制閥流場分析

王 銘a,張世富b,張起欣c,張冬梅b,陳 暢a

(中國人民解放軍后勤工程學院 a.研究生大隊； b.國家救災應急裝備工程技術研究中心；c.軍事供油工程系， 重慶 401331

為分析油罐自力式液位控制閥主閥關閥過程，對閥芯運動過程中的受力進行分析并建立微分方程，得出關閥過程中閥芯的速度及加速度公式，分析了關閥速度與閥門工作壓力、節流孔大小、閥芯內腔平均直徑的關系。借助Fluent軟件對關閥過程中閥門的速度場和壓力場進行研究，并進一步分析了閥芯受力。研究結果對于主閥結構的優化及錐閥閥芯結構的設計具有參考意義。

油罐；閥門；液位控制；受力分析；仿真

油料的液位控制一般指對儲液容器液位進行控制調節，使其保持在某一高度下，避免發生溢油及其他安全事故[1-3]。目前，各類液位控制裝置廣泛應用于油罐車、計量罐、地面及地下儲油罐等各類容器，在減小工人工作強度、保障生產安全方面發揮了巨大的作用，是現代油料生產及儲存計量中必不可少的裝置[4]。為解決液位控制問題，研制了一種全機械自力式高液位控制先導閥，并對其進行動力學分析，利用Fluent軟件對受控主閥關閥過程進行數值模擬。

1 高液位控制閥工作原理

高液位控制閥的工作原理如圖1所示。閥門正常工作時，手柄被拉至最右端，真空閥閥芯打開，啟閉導閥開啟。閥門入口處的高壓流體經流道①進入真空閥，經過文氏管時產生負壓，外界空氣經氣道③進入真空閥，和流體一起經流道⑥排出。由于氣道②與③聯通，真空執行器的膜片右側氣室壓力下降較小，鋼珠受到真空執行器殼體與錐芯軸的卡阻，與鎖桿位置被固定。主閥入口處流體經節流孔進入閥芯內腔，再由流道④、⑤流入主閥流道。在節流孔的節流減壓作用下，閥芯內腔的壓力小于閥芯前液體的壓力，閥門一直處于打開的狀態。當受油油箱的液位達到氣道③的下沿時，氣道被阻塞，膜片右側氣室壓力下降，膜片帶動錐芯軸右移，鋼珠的卡阻消失，彈簧推動真空閥閥芯左移直至關閉。在手柄的作用下，鎖桿左移，同時啟閉導閥關閉，閥芯內腔液體不再流動，其壓力與閥芯前端壓力相同，在主閥閥芯彈簧的作用下，閥芯上移直至閥門關閉。

2 關閥過程動力學分析

關閥過程從啟閉導閥完全關閉控制流道開始，閥芯前腔液體通過節流孔進入閥芯內腔，液體在閥芯內腔堆積，推動閥芯向右運動。主閥的結構如圖2所示。

1.主閥； 2.節流孔； 3.啟閉導閥 ；4.真空閥 ；5.手柄；

6.真空執行器 ；7.錐芯軸 ；8.鋼珠； 9.鎖桿

圖1 高液位控制閥原理

圖2 主閥結構

閥芯關閥運動方程[5-6]：

P2S2+Fk-P3S3-P1S1-Ft-

(1)

式中：P1、S1為閥芯前腔壓強(Pa)及受壓面積(m2)；P2、S2為閥芯內腔壓強(Pa)及受壓面積(m2)；P3、S3為閥芯側腔壓強(Pa)及受壓面積(m2);Fk為彈簧彈力(N)；Ft為穩態液動力(N)；Fw為瞬態液動力(N)；f為閥芯與閥體間的摩擦力(N)；M為閥芯質量(kg);x為閥芯位移(m);β為黏性阻尼系數；t為關閥時間(s)。

(2)

瞬態液動力Ft是因油料流速和流動方向的變化作用在閥芯上的力，可根據動量方程求出：

(3)

式中：Q為主閥口流量(m3/s);α′為動量修正能系數，對于湍流一般取1。

閥門流量方程為

(4)

式中：C1為閥門的流量系數；A1為閥門出口的垂直過流面積(m2)；ΔP1為閥門前后壓差(Pa)，由閥門工作條件決定，假設閥門出口壓力為P0，則ΔP1=P1-P0； ρ為油料密度(kg/m3)。

根據閥芯頂部結構的幾何關系，垂直過流面積為

(5)

式中：x0為閥門工作位置距關閥位置距離，即關閥總位移;θ為閥芯錐角。

根據質量守恒原理，Q進=Q出，得出v進與v出的關系為

(6)

將式(4)、(6)代入式(3)，得出瞬態液動力的表達式：

(7)

根據閥芯前部錐形結構，油料通過節流孔進入閥芯內腔的流量為

(8)

式中：C2為節流孔的流量系數；A2為節流孔垂直過流面積(m2)；ΔP2為節流孔前后壓差(Pa)，等于P1-P2。

(9)

由于閥芯的運動速度很小，可以認為閥芯側腔的壓力P3與導壓孔處流道內流體的壓力相同。取從閥門入口到導壓孔處的流道為控制體，主閥流量也可表示為

(10)

式中：C3為控制體的流量系數；A3為導壓孔處主閥流道的垂直過流面積(m2)；ΔP3為控制體前后壓差(Pa)，等于P1-P3。

將式(4)與式(10)聯立，得出：

P3=(1-k)P1+kP0

(11)

彈簧彈力

(12)

式中：x1為閥門處于關閉狀態時彈簧的壓縮量；K為彈簧的彈性系數。

將式(7)、(8)、(9)、(11)、(12)代入式(2)可得

(13)

為簡化方程，令：

則式(2)可表示為

(14)

式(14)是二階非線性常微分方程，求解可得：

(15)

將式(15)對時間進行求導，得到加速度的表達式：

(16)

式(15)、 (16)說明：此高液位閥閥芯的在關閥過程中速度隨關閥行程不斷增加；關閥運動是一個受到多方面因素綜合影響的運動：閥芯內腔平均直徑和閥芯的質量越大，速度和加速度越小；閥門前后壓力和節流孔的面積越大，速度和加速度越大。

3 數值模擬

3.1 建立模型及網格劃分

此閥門結構較為復雜，為簡化計算，在仿真時利用二維模型代替三維結構，對流動情況較為簡單的部位的結構進行簡化，如閥芯側腔、閥芯內腔等。利用Solidworks軟件建立二維模型，利用ICED軟件進行網格的劃分。閥芯內腔采用非結構網格，閥其余部分采用結構網格，并在節流孔入口處設置interface，在流動情況較為復雜的流道入口處加密網格。模型采用動網格模擬，網格重構采用Smoothing、Layering、Remeshing 3種方法。在Fluent軟件中對初始網格(圖3(a))進行預覽，得到關閥過程中各個時刻的網格劃分情況。關閥過程中幾個典型位置的網格情況如圖3所示。

3.2 基本控制方程

連續性方程：

(17)

式中：ρ為流體的密度(kg/m3) ;t為時間(s)；u、υ、ω為x、y、z方向的速度分量(m/s)。

動量方程：

(18)

(19)

(20)

式中:p為流體的靜壓強(Pa)；τij為應力張量的9個分量(Pa)，i、j分別可取x、y、z;f為體積力的3個坐標分量(N)。

能量方程：

(21)

圖3 關閥行程特殊位置網格

3.3 Fluent軟件自定義函數(UDF)的編寫與求解器參數設置

利用UDF宏DEFINE_CG_MOTION定義閥芯在每一個時間步的位移。閥芯單個時間步的速度增量可按式(17)進行計算：

dv=F·dt/M

(22)

式中：dv為dt時間內閥芯速度的變化量；F為閥芯所受合力；dt為時間步長度；M為閥芯質量。

設置靜態變量對速度變化量進行累加，利用速度與時間步長度的乘積定義單個時間步的位移。

利用UDF宏DEFINE_EXECUTE_AT_END定義在每個時間步的迭代計算完畢后輸出流量、壓力等相關的參數。

選擇基于壓力的分離求解器(SIMPLEC算法)作為模型求解器。選擇標準k-ε湍流模型，流體材料設置為水。導入自定義函數(UDF)，并以此定義相關區域為剛體(Rigid Body)運動。根據邊界運動設置deforming類型運動邊界。設置閥門入口為壓力入口邊界條件，根據閥門實際應用的壓力，設置入口壓力為106Pa。設置閥門出口為壓力出口邊界條件，壓力設置為105Pa，并將所有的interface進行配對。時間步設置為10-5s。其他設置采用默認值，初始化后進行計算。

3.4 仿真結果及分析

仿真得到了閥門流量隨關閥行程變化的曲線，如圖4所示，可見流量隨閥門開度的減小而減小。通過積分計算，得出關閥過程中的流量總計為0.022 5 m3。

圖4 閥門流量隨關閥行程變化關系

關閥行程中特殊位置的速度云圖如圖5所示。從圖中可以看出：在閥芯前腔關閥行程前期，最大速度出現在流道最為狹窄的區域；節流孔處速度較小，這是由于閥芯運動速度較慢，也說明了此閥關閥平穩。從壓力云圖還可以看出：閥芯前腔壓力沒有出現大的波動。

圖5 關閥行程中特殊位置的速度云圖

仿真得到了閥芯所受合力隨關閥行程變化的曲線，如圖 6 所示。閥芯所受合力決定了其運動規律，從圖中可以發現：閥芯所受合力隨關閥行程先上升后下降，并一直處于波動的狀態。這說明閥芯的運動為速度恒大于0的加速運動，加速度先增大后減小。在關閥過程中出現了水擊，但由于關閥速度較慢，所以壓力波動并不大。

壓力關閥行程中特殊位置的壓力云圖如圖7所示。從圖中可以看出：閥芯前腔和閥芯內腔的壓力較高，這是由于閥門入口處的速度較小，所以壓降較小；壓力在流道入口前端開始下降，并隨著閥門的關閉下降越來越快，這是因為流道前端的局部阻力越來越大；閥芯內腔的壓力與導壓孔外流體壓力基本保持一致，且在關閥過程中一直小于0，這是由于流道內流體速度增大造成的。閥芯錐頂部平面壓力較大，明顯要高于錐角部分，從速度云圖(圖5)中也可以看出：由于受到與水流方向垂直的錐頂平面的阻礙，閥芯錐頂平面及前部區域速度較小，所以壓力相對于錐角部分小。閥門流道部位出現了極低的壓力，說明在流道處流體會發生氣化。

圖6 閥芯受力隨關閥行程變化

圖7 關閥行程中特殊位置的壓力云圖

4 結束語

通過理論分析及仿真，分析了主閥關閉過程中的流場及壓力情況，說明此自力式閥門能利用液體壓力自動關閥，且關閥平穩，沒有出現較大的壓力波動，不會因水擊損壞閥門結構。閥芯的關閥運動方程表明：閥芯內腔平均直徑和閥芯的質量越大，閥芯速度和加速度越小；閥門前后壓力和節流孔的面積越大，閥芯速度和加速度越大。Fluent軟件仿真結果表明：關閥過程出現了水擊，但壓力波動不大，說明可以適當增大節流孔直徑或減小閥芯質量。

[1] 王繼承.液位報警器在密閉計量罐上的應用[J].中國石油和化工標準與質量，2013，33(3)：97.

[2] 斯建華.重油罐液位控制方案的改進[J].銅業工程,1999(2):17-19.

[3] 郭玉瓊.高線B潤滑站氣壓貯油罐液位控制的改進[J].浙江冶金,2003(4):49-51.

[4] 黃增進,梁朝旭.關于儲罐液位報警值設定的探討[J].廣東化工,2015,42(10):157-158.

[5] 張超.減速頂壓力閥的特性分析[J].減速頂與調速技術，2004，81(4)：3-6.

[6] 唐曉寅.工程流體力學[M].重慶：重慶大學出版社，2007.

(責任編輯 劉 舸)

Flow Field Analysis of Self-Operated High Level Control Valve of Oil Tank

WANG Minga, ZHANG Shi-fub, ZHANG Qi-xinc, ZHANG Dong-meib, CHEN Changa

(a.Graduate Student Brigade; b.Research Center of National Disaster Emergency Equipment Engineering Technology; c.Department of Military Oil Supply Engineering,Logistic Engineering University of PLA, Chongqing 401331, China)

In order to study the shutoff process of high level valve for oil tank, the formula of the velocity and acceleration of the valve element, which is inferred from the analysis of forces executing on valve and differential equation, based on which the relationship between working pressure, orifice size, valve element lumen mean diameter and shutoff velocity of valve element were analyzed. Software Fluent was used to simulate the velocity and pressure field during the shutoff process. Further analysis of the force executing on valve element was done, which is of good value for the majorization of valve structure and the design of valve element of cone valve.

oil tank; valve; level control; force analysis; simulation

2016-11-28

王銘(1992—)，男，碩士研究生，主要從事工程流體力學研究，E-mail:1832778276@qq.com。

王銘,張世富,張起欣，等.油罐自力式液位控制閥流場分析[J].重慶理工大學學報(自然科學)，2017(4):82-87.

format：WANG Ming, ZHANG Shi-fu, ZHANG Qi-xin,et al.Flow Field Analysis of Self-Operated High Level Control Valve of Oil Tank[J].Journal of Chongqing University of Technology(Natural Science)，2017(4):82-87.

10.3969/j.issn.1674-8425(z).2017.04.013

TH816

A

1674-8425(2017)04-0082-06