差動式快開閥開閥過程動力學分析與仿真

2013-09-05 孫振華 后勤工程學院軍事供油工程系

  以電磁閥為先導(dǎo)閥,自力式差壓控制的快開閥用于模擬產(chǎn)生管線發(fā)生破裂時的瞬間泄露工況,在設(shè)計時間內(nèi)快速、穩(wěn)定的開啟是該閥研制的關(guān)鍵技術(shù)難點。分析了差動式快開閥的機構(gòu)及開啟原理,以流體力學基本理論為基礎(chǔ),建立了壓降隨開度變化的關(guān)系,完成了閥門開啟過程動力學分析與仿真計算,得到了不同閥門開度下的開啟時間、速度及加速度,閥芯運動過程中的驅(qū)動力逐漸減小,開啟時間、開啟速度及加速度三者隨開度的變化呈非線性變化關(guān)系。計算結(jié)果表明,閥體開閥時間與實驗實測數(shù)據(jù)吻合良好,開閥過程的速度及加速度與實驗結(jié)果一致,為工程設(shè)計中主要部件的優(yōu)化研究提供了理論依據(jù)。

  泄漏是管線常見的事故工況,管線泄漏不僅造成油料的損失,影響管線的正常輸送,甚至會導(dǎo)致管線輸送中斷,直接影響戰(zhàn)爭后勤油料保障;大量的油料泄漏還會造成嚴重的環(huán)境污染,危及人民的生命財產(chǎn)安全。因此,必須對管線使用和管理人員進行泄漏事故工況的針對性訓(xùn)練,提高處置泄漏事故的能力。

  為此,可設(shè)計研制一套訓(xùn)練輔助裝置,模擬產(chǎn)生瞬間泄漏等事故工況。該裝置突發(fā)泄漏工況的模擬關(guān)鍵在于快開閥的研制。目前國內(nèi)外高速開關(guān)電磁閥控制流量較小,單一依托電磁閥無法滿足該訓(xùn)練輔助裝置控制需求。本文研究的快開閥采用電磁控制先導(dǎo)閥,主閥利用差壓控制原理,既可有效提高控制流量,其開啟響應(yīng)時間也能滿足控制需求。由于該閥在輸油現(xiàn)場使用,不宜采用電氣原理,而應(yīng)采用科學、穩(wěn)定、可靠的液體與機械控制技術(shù)來實現(xiàn)高速開啟,其中閥體的開閥過程是研制該閥的關(guān)鍵性技術(shù)。本文重點介紹了差動式快開閥開閥過程的動力學原理、閥體結(jié)構(gòu)和仿真計算,并對計算結(jié)果進行了分析和討論。

1、閥體結(jié)構(gòu)與開閥過程動力學分析

  閥體結(jié)構(gòu)如圖1所示?扉_閥未動作時,電磁先導(dǎo)閥關(guān)閉,流體經(jīng)孔1流入并充滿閥芯內(nèi)部空腔,流體對閥芯底部壓力等于閥芯內(nèi)部流體壓力,彈簧預(yù)壓力等于底座對閥芯的支持力,閥芯受力平衡。當電磁先導(dǎo)閥通電開啟,流體經(jīng)孔2流出。由于孔2直徑遠大于孔1直徑,流體在閥芯內(nèi)部無法聚壓,閥芯底部受力大于腔內(nèi)流體壓力、彈簧力與流體粘性摩擦力的合力,使閥芯向上運動,閥門開啟,大部分流體從閥芯右側(cè)流出。閥門的開啟程度由限開度機構(gòu)控制。

快開閥結(jié)構(gòu)簡圖

圖1 快開閥結(jié)構(gòu)簡圖

閥內(nèi)空腔劃分圖

圖2 閥內(nèi)空腔劃分圖

2、閥體開閥時間計算

  如圖2所示,閥內(nèi)空腔直徑不同,為了便于分析,將閥內(nèi)空腔(有流體流經(jīng)的部分)分為a,1,2,3,4,5幾個部分,流體在每個部分流動都相當于在圓管內(nèi)流動,只是長度和直徑有所區(qū)別。a部分流量、上下壓差、長度、直徑分別表示為Qa,ΔPa,La,da;1部分可表示為Q1,ΔP1,L1,d1,其他依次類推。

  閥門在開啟過程中,空腔內(nèi)的流體經(jīng)孔2排出,腔3內(nèi)流體排出時間即為閥門開啟時間。

  2.1、閥芯內(nèi)腔流量計算

  主閥閥體是傾斜的,閥芯空腔內(nèi)劃分的各個部分可看成為非水平管。非水平管道層流時,流體重力只影響位能,不影響流動特性。因此,在只討論速度分布的情況下,閥芯空腔內(nèi)各個部分流體速度分布可按水平管處理,根據(jù)水平管的N-S方程,經(jīng)進一步計算,水平管內(nèi)過該截面的流量可表示為:

差動式快開閥開閥過程動力學分析與仿真(1)

  式中,d為水平管的直徑。

  2.2、閥芯空腔各部分壓差計算

  根據(jù)圓管流量式(1)得到閥空腔內(nèi)6個部分的流量,分別為:

差動式快開閥開閥過程動力學分析與仿真

  在閥芯向上運動的過程中,a,1,2部分的流體液面不會發(fā)生變化,認為流量沒有變化,即Qa=Q1=Q2。根據(jù)流量關(guān)系得:

差動式快開閥開閥過程動力學分析與仿真 (2)

  閥芯向上運動過程中,腔3內(nèi)流體經(jīng)孔2中排出,腔3內(nèi)液位會有明顯變化;腔4,5內(nèi)液位不會產(chǎn)生變化,直到閥門達到最大開度。若dt時間內(nèi),閥芯沿軸線方向上升dx,即腔3內(nèi)有dx長度的流體流出,考慮到從腔2內(nèi)進入的流量,可表示為:

差動式快開閥開閥過程動力學分析與仿真 (3)

  由于腔4,5內(nèi)的凈剩流量沒有變化,則,Q3=Q4=Q5,由6個部分流量表達式及式(3),可以得到下面關(guān)系式:

差動式快開閥開閥過程動力學分析與仿真

  2.3、閥芯腔內(nèi)壁面粘性摩擦力計算

  閥門在開啟過程中,流體受擠壓經(jīng)孔2中排出,閥芯內(nèi)腔壁面會對流體產(chǎn)生一個阻礙其運動的力,這個力就是流體對閥芯內(nèi)腔內(nèi)壁面產(chǎn)生的粘性摩擦力。

  根據(jù)牛頓內(nèi)摩擦定律,可得管中流體產(chǎn)生的切應(yīng)力為:

差動式快開閥開閥過程動力學分析與仿真

  在管壁處τ=τ0,r=R=d/2,則有

差動式快開閥開閥過程動力學分析與仿真

  于是流體對管壁的摩擦力可表示為:

差動式快開閥開閥過程動力學分析與仿真 (4)

  根據(jù)式(4)可以得到,流體對閥空腔內(nèi)6個部分的摩擦力為:

差動式快開閥開閥過程動力學分析與仿真 (5)

  2.4、壓差力計算

  為便于計算壓差產(chǎn)生的力,將閥芯內(nèi)腔部分由軸線方向向兩側(cè)劃

差動式快開閥開閥過程動力學分析與仿真

  幾個區(qū)域,區(qū)域劃分以此前劃分的6個部分的直徑大小(da  經(jīng)計算,各部分壓差力分別為:

差動式快開閥開閥過程動力學分析與仿真 (6)
差動式快開閥開閥過程動力學分析與仿真

圖3 計算壓差力時閥內(nèi)空腔劃分圖

  2.5、閥門開啟時間分析

  閥門在開啟過程中,彈簧受力壓縮,彈簧壓力f的大小可近似認為等于流體壓差力與流體對閥芯內(nèi)壁面粘性摩擦力的合力,表示為:

差動式快開閥開閥過程動力學分析與仿真 (7)

  將式(5)、(6)代入式(7),得:

差動式快開閥開閥過程動力學分析與仿真 (8)

  將式(2)代入式(8),進一步化簡,得:

差動式快開閥開閥過程動力學分析與仿真

  (9) 式中:

差動式快開閥開閥過程動力學分析與仿真

  彈簧產(chǎn)生的壓力f可表示為:

差動式快開閥開閥過程動力學分析與仿真

  式中:k表示彈簧的彈性系數(shù);x表示主閥開度,即閥芯沿軸線上升的距離;

差動式快開閥開閥過程動力學分析與仿真

  表示彈簧的預(yù)壓量。將上式代入式(9),得:

差動式快開閥開閥過程動力學分析與仿真 (10)

  閥門在開啟過程中,a部分流速很大,壓差ΔPa比其他幾個部分的壓差大得多,可近似認為a部分的壓差ΔPa同閥芯上下壓差ΔP(閥內(nèi)6個部分壓差之和)相等,故式(10)又可表示為:

 (11)

  閥芯向上運動過程中,閥芯上、下的壓力差ΔP與主閥開度相關(guān),可表示為ΔP=f(x),其中,x為主閥開度,函數(shù)關(guān)系f需通過實驗測得。則式(11)可變化為:

差動式快開閥開閥過程動力學分析與仿真 (12)
差動式快開閥開閥過程動力學分析與仿真 (13)

  對式(13)積分,得:

差動式快開閥開閥過程動力學分析與仿真""

  當x=0時,t=0,則

差動式快開閥開閥過程動力學分析與仿真

  上式可進一步簡化,得閥門開啟時間為:

差動式快開閥開閥過程動力學分析與仿真(14)

  本文所研究的閥門跟截止閥原理相同,根據(jù)研究,閥門開啟時壓降與主閥開度之間的關(guān)系ΔP=f(x)滿足雙曲線變化關(guān)系,可設(shè):

差動式快開閥開閥過程動力學分析與仿真(15)

  式中:P0為閥門開啟前主閥底部壓力,近似為管線運行壓力;A,B,C是閥門壓降與開度變化關(guān)系的系數(shù),可通過實驗測得。

  閥芯上、下的壓降隨開度變化關(guān)系可經(jīng)實驗測得。在式(14)中代入不同的開度x值,可得不同開度情況下的閥門開啟時間。

  2.6、分析結(jié)果

  簡化閥門結(jié)構(gòu),并進行一定的假設(shè),采用數(shù)值模擬的方法得到式(15)的系數(shù)值;式(14)中,除閥門開度x值外,其余都是已知數(shù),代入不同的主閥開度x值,就可得到不同開度情況下的閥門開啟時間理論值。樣機實驗中,利用瞬變數(shù)據(jù)采集系統(tǒng),通過位移傳感器,閥芯向上每運動2mm記錄一次時間數(shù)據(jù),得到不同開度下閥門開啟時間實驗值。結(jié)果見表1。

表1 閥門開啟時間理論值、實驗值對比

差動式快開閥開閥過程動力學分析與仿真

  為了直觀地反映理論值與實驗值之間的關(guān)系,將部分數(shù)據(jù)擬合成曲線,見圖4?梢钥闯,閥門開啟時間隨開度的變化近似呈線性變化,理論值與實驗值吻合良好。

3、仿真計算分析

  3.1、開啟速度分析

  閥門的開啟速度即為腔3內(nèi)液柱長度沿軸線方向的減小速度。根據(jù)式(12)得,閥門開啟速度為:

差動式快開閥開閥過程動力學分析與仿真(16)

  式(16)中代入不同的閥門開度x值,即為不同開度時的閥門瞬時速度。

  3.2、開啟加速度分析

  閥門開啟加速度是閥門開啟速度對時間的倒數(shù),即加速度

差動式快開閥開閥過程動力學分析與仿真

  將閥門開啟速度關(guān)系式(16)對時間求導(dǎo),得:

差動式快開閥開閥過程動力學分析與仿真 (17)
差動式快開閥開閥過程動力學分析與仿真

圖4 閥門開啟時間理論值與實驗值對比曲線

  式(17)中代入不同的閥門開度x值,即為不同開度時的閥門瞬時加速度。

  3.3、分析結(jié)果

  將上述理論值與實驗值進行對比,結(jié)果見表2,3。為了直觀地反映理論值與實驗值之間的關(guān)系,將部分數(shù)據(jù)擬合成曲線,見圖5,6。由圖5,6可以看出,閥門開啟過程的速度隨開度的變化逐漸增大,開度越大,速度變化越平穩(wěn);加速度隨開度的變化逐漸減小。理論值與實驗值吻合較好。

表2 閥門開啟瞬時速度理論值、實驗值對比

閥門開啟瞬時速度理論值、實驗值對比

表3 閥門開啟瞬時加速度理論值、實驗值對比

閥門開啟瞬時加速度理論值、實驗值對比

圖5 閥門開啟瞬時速度理論值與實驗值對比曲線

閥門開啟瞬時速度理論值與實驗值對比曲線

圖6 閥門開啟瞬時加速度理論值與實驗值對比曲線

  由速度和加速度變化規(guī)律分析可得,閥門開啟過程閥芯運動穩(wěn)定,未受到較大沖擊,達到設(shè)計要求,滿足工程學應(yīng)用需求。

4、結(jié)論

  1)利用N-S方程解決了閥門開啟過程中液體流速的精確解問題,從而求出了閥芯在不同開度下的開啟時間;

  2)給出了開閥過程的仿真計算原理、模型、過程和計算結(jié)果;

  3)計算結(jié)果表明,閥門開啟到不同開度時,閥芯開啟時間與實驗實測數(shù)據(jù)吻合良好,開啟過程的速度及加速度與實驗結(jié)果一致。