水擊

有壓管道中壓強波動的現象

水擊,是在有壓管道中,液體流速發生急劇變化所引起的壓強大幅度波動的現象。管道系統中閘門急劇啟閉,輸水管水泵突然停機,水輪機啟閉導水葉,室內衛生用具關閉水龍頭,都會產生水擊。水擊壓強水擊的基本問題是最大壓強的計算,最大壓強一般出現在發射波斷面(如閥門處)。水擊壓力波的傳播過程如圖所示。水擊基本方程式一般水擊壓強計算可採用下列方程組,結合起始條件和邊界條件求積分。

機理


封閉管道中的水擊過程是壓力波傳播的過程。假如管壁是剛性的,就水的可壓縮性來說,壓力波同聲波在水中傳播速度是一樣的,可表示為:
,式中K為水的體積彈件模量;ρ為水的密度。在常溫下,,實際上封閉管道中的壓力波傳播速度還受管壁彈性的影響,因而上式變成:計算方法
式中E為管壁材料彈性模量;D和e為管徑和壁厚。
水擊壓力波的傳播過程如圖所示。設從水庫引水的管道長度為L,末端裝有快速啟閉閥門。如果忽略水力損失,則管道末端的初始水頭H等於水庫水頭,管內初始流速。從圖上可以看到閥門瞬間關閉后發生的水擊現象。閥門關閉后,緊靠閥門處管段Δx的流速首先等於零。由於水流的慣性,水體被壓縮,管壁膨脹,水頭增加ΔH並以波速c向上游傳播,使壓頭增至。
當時,壓力波到達進口,此時整個管內壓頭為 ,流速為零。由於管內壓頭比水庫高ΔH,為了保持平衡,管內水體倒流併產生以波速c向下游傳播的降壓反射波,使進口壓頭恢復到初始狀態H。
當時,反射波到達閥門,整個管內壓頭為H,但流速v0朝向上游。由於慣性,在閥門處產生降壓ΔH的反射波,以速度c向上游傳播,使管內壓頭降至 。
當時,整個管內壓頭為 ,流速為零。因管內壓頭比水庫低ΔH,水庫水又流入管內,朝向下游,在進口處產生增壓ΔH的反射波以波速c向下游傳播,進口壓頭回到H。
當時,壓力波到閥門,整個管內流速和水頭恢復到初始狀態,完成了一個壓力振蕩周期,以後水擊現象又重複上述過程。
因為存在水力損失,水擊壓力振幅實際上隨時間衰減,並非保持不變。

正文


水擊可導致管道系統的強烈震動,產生雜訊和氣穴。它是促使管道破裂的最經常的因素。掌握水擊壓強的變化規律對輸水管道的設計,對消減水擊的破壞作用,有很大的實際意義。
理想水擊波動過程設想流速突然變化、波動過程無能量損失的水擊模型,管道長為L,上游接水池,當下游閥門突然關閉時,靠近閥門的液體流速由減低至零,壓強增大。這種減速增壓的增壓波以彈性波速上行傳播,歷時 到達上游末端,以壓強恢復為原始壓強的常壓波向下行傳播,經歷同樣的時間( )回到閥門。此後再以壓強低於原值的減壓波上行傳播,又以常壓波下行回到閥門。完成整個水擊周期,共需時間
。以後將不斷地重複上述波動過程。實際上閥門不是瞬時關閉,因此實際波動可以看作一系列微段時間閥門關閉所導致一系列微小上行波,及微小下行波之和。在波動過程中,粘性作用消耗能量,使波動幅值逐漸趨於消失。
水擊壓強水擊的基本問題是最大壓強的計算,最大壓強一般出現在發射波斷面(如閥門處)。根據上述波動過程可以繪出理想水擊時閥門處的壓強隨時間變化圖。當 時,壓強突然降低,這是由於反射波回到閥門。實際水擊過程仍然存在反射波改變閥門壓強的作用。所以,當閥門關閉尚未到達 的時間時,閥門處流速不斷降低,壓強不斷增高。當閥門關閉時間超過 時,由於反射波已回到閥門,使壓強降低。故有直接水擊和間接水擊的不同計算方法。
當閥門關閉時間 ,產生直接水擊。直接水擊的壓強升高
,可取微段寫能量方程導出下式:
式中
為流速降低值;
為液體密度。
當閥門關閉時間
為間接水擊。間接水擊的計算需要知道流速隨時間變化的關係,可用 時間內的最大流速減低值代入上述求壓強增量
水擊基本方程式一般水擊壓強計算可採用下列方程組,結合起始條件和邊界條件求積分。或在電子計算機上求解。
水擊運動微分方程 和水擊連續性微分方程
在不計摩阻並略去高階微量之後
式中
s為沿管長距離;
a為管軸傾角;
D為管徑;
H為測管水頭高度;λ為沿程阻力係數;ɡ為重力加速度

防止水擊危害的措施

水擊壓強是巨大的,這一巨大的壓強可使管路發生很大的變形,甚至爆炸。為了預防水擊的危害,可在管路上設置空氣室,或安裝具有安全閥性質的水擊消除閥,如安全閥。具體措施如下:
(1)適當加大管徑,限制管道流速。一般在液壓系統中把速度控制在4.5m/s以內,△P不超過5MPa就可以認為是安全的。
(2)正確設計閥口或設置制動裝置,使運動部件制動時速度變化比較均勻。
(3)延長閥門關閉和運動部件制動換向的時間,可採用換向時間可調的換向閥
(4)盡量縮短管長,以減小壓力衝擊波的傳播時間。
(5)在容易發生液壓衝擊的部位採用橡膠軟管或設置蓄能器,以吸收衝擊壓力;也可以在這些部位安裝安全閥,以限制壓力升高。

基本方程


假設管壁不變形,且忽略水力損失項和非線性項,並令(γ為水的比重),則管道非定常壓力流的基本方程(見壓力流)可簡化為:
式中v和H為瞬變流速和水頭;c為波速;x和t為距離和時間;g為重力加速度。上式一般積分的形式為:
式中Φ為順x軸的正向波;F為逆x軸的反向波;為管內的初始流速。
由於計算機的普遍使用,現多採用數值計算方法,如特徵線法,隱式有限差分法和有限體積法。

計算方法


(1)解析法
解析法的基本原理是,將水擊基本方程組化為波動方程並求其通解,結合邊界條件及初始條件可逐步求得任意時刻的任意斷面上的水擊壓強,其物理意義明確,應用簡單,多適用於不計阻力的管路系統。
(2)圖解法
圖解法以解析法導出的關係式為基礎進行圖解。圖解法的複雜程度由管道系統的複雜性來決定,管道越複雜且摩阻損失佔比重越大、水擊波反覆傳播次數越多,圖解法計算過程也越複雜,同時計算精度也較低,限制了對工程水擊問題的研究和解決。
有限差分法(Finite Difference Calculus)用差商來代替控制方程中的導數,解出網格點上的函數值從而確定流體的參數。有限差分法分為顯式和隱式兩種格式。隨著差分方法的發展,產生了眾多不同精度的差分格式,但在收斂性、穩定性等方面,有待解決的問題仍然不斷出現。
(4)特徵線法
該方法將考慮管道摩阻的水擊偏微分方程沿其特徵線變換為全微分方程,再近似的變成差分方程,然後進行數值計算。這種方法便於建立求解微分方程的穩定準則,計算精度高,容易處理邊界條件,可用於複雜系統的計算,編製程序也比較方便。目前特徵線法是水擊計算最主要方法之一。
有限體積法通常用來求解雙曲線方程的情況,比如氣體動力學以及淺水波的問題。採用有限體積法的優勢主要在於(1)能夠保持質量方程和動量方程的完整性,(2)對沒有振蕩的非連續性問題能有很好的解。