“山”形自閉式水封水密特性仿真計算與試驗(yàn)研究

楊斌 周雨微 王蒂

摘要： 為改善傳統(tǒng)充壓伸縮式水封機(jī)械加壓方式所帶來的閘門結(jié)構(gòu)復(fù)雜和不耐久的問題，設(shè)計了一種利用水封在庫水作用下自身的大變形主動完成封水的自閉式新型水封，從而取消機(jī)械加壓，這種水封設(shè)計的關(guān)鍵是對各種工況下水封的變形進(jìn)行準(zhǔn)確預(yù)估。基于水封橡膠材料的非線性張量本構(gòu)關(guān)系、幾何大變形非線性特征、接觸非線性特征、以某尖頭“山”形水封為例，利用連續(xù)介質(zhì)有限元計算方法對該種水封進(jìn)行了自閉的仿真計算，并以室內(nèi)試驗(yàn)加以驗(yàn)證。實(shí)踐表明：“山”形橡膠水封能實(shí)現(xiàn)較好的自閉，且大大削減了閘門充壓水封裝置的建造成本、便利了工程維護(hù)、提高了耐久性，具有較大的應(yīng)用價值。

關(guān) 鍵 詞： 自閉式水封; 水封變形; 封頭接觸應(yīng)力; 仿真計算; 連續(xù)介質(zhì)力學(xué)

中圖法分類號： ?TV131.6

文獻(xiàn)標(biāo)志碼： ?A

DOI： 10.16232/j.cnki.1001-4179.2022.03.025

?? 0 引 言

充壓伸縮式水封是近20 a才發(fā)展起來的一種高壓閘門水封，目前采用的型式是在水封座板下聯(lián)通一根充氣管，通過機(jī)械充氣加壓使封頭外伸與止水面板接觸緊密，從而達(dá)到止水的目的? [1] 。這種機(jī)械加壓方式使得閘門結(jié)構(gòu)復(fù)雜，加壓設(shè)備和管道需要進(jìn)行維護(hù)與更換?；诖耍O(shè)計出了一種水封裝置取消了加壓設(shè)備和管道，可實(shí)現(xiàn)直接利用庫水壓力作為水封的常用背壓進(jìn)行“自閉式”封水（見圖1），該水封命名為自閉式水封。

這種封水方式利用水封橡膠材料的大變形特點(diǎn)，而對于大變形，事前要作準(zhǔn)確估計和精確計算，以避免設(shè)計失誤。

1 自閉式水封仿真計算基本理論

水封橡膠材料的大變形涉及材料? [2] 、幾何? [3] 、接觸? [3] 三重非線性，是計算的關(guān)鍵所在? [4-7] 。水封橡膠材料的本構(gòu)關(guān)系見式（1），材料非線性指應(yīng)力張量 S 與變形張量 B 不是線性關(guān)系，因此通過擬合應(yīng)變能函數(shù) W e 進(jìn)行計算? [8-11] 。幾何非線性見式（2），指變形的幾何分解中存在位移對質(zhì)點(diǎn)隨體坐標(biāo)偏導(dǎo)的二次項(xiàng)? [12-14] 。

S =2[（ ?W e ?I 1 +I 1 ?W e ?I 2 ） B - ?W e ?I 2? B ?2]-p I? （1）

式中： p為靜水壓力; I 為單位張量; S 為應(yīng)力張量; B 為變形張量;W e 為應(yīng)變能函數(shù)。

ε? ij =ε L? ij +ε N? ij = 1 2 （ ?u j ?ξ i + ?u i ?ξ j ）+ 1 2 （ ?u k ?ξ i? ?u k ?ξ j ） （2）

式中： ε? ij 為9個應(yīng)變張量分量，i、j、k為指標(biāo)序號，符合愛因斯坦求和約定;ε L? ij 、ε N? ij 為9個線性應(yīng)變張量分量和9個非線性應(yīng)變張量分量;u i、u j為位移矢量分量;ξ i、ξ j 為物質(zhì)坐標(biāo)系的質(zhì)點(diǎn)隨體坐標(biāo)。

幾何非線性與材料非線性耦合在一起，處理方法為將應(yīng)變線性部分與非線性部分分離，如式（3）～（5），再代入應(yīng)變能函數(shù)三項(xiàng)Ogden模型（8） （因計算要求有相當(dāng)高的精度，故采用了多項(xiàng)級數(shù)模型）和邊界條件，構(gòu)成有限元模型（9） 。

ε L=SymbolQC@ Lu =SymbolQC@ LNa ?e= B ?L a ?e （3）

ε N= 1 2? Aθ = 1 2? AGa ?e= B ?N a ?e （4）

ε=ε L+ε N=（ B ?L+ B ?N） a ?e=? B ???～? a ?e （5）

其中算子：

SymbolQC@ L =?? ???ξ 1? 0 0 0? ???ξ 3?? ???ξ 2 ?0? ???ξ 2? 0? ???ξ 3? 0? ???ξ 1 ?0 0? ???ξ 3?? ???ξ 2?? ???ξ 1? 0? ??T? （6）

矩陣：

A =?? ???u ???T ??ξ 1? 0 0 0? ???u ???T? ?ξ 3?? ???u ???T? ?ξ 2 ?0? ???u ???T? ?ξ 2? 0? ???u ???T? ?ξ 3? 0? ???u ???T? ?ξ 1 ?0 0? ???u ???T? ?ξ 3?? ???u ???T? ?ξ 2?? ???u ???T? ?ξ 1? 0? ??T

θ =?? ??u? ?ξ 1 ????u? ?ξ 2?? ??u? ?ξ 3?? ??T

G =?? ?N 1 ?ξ 1 I 3? ?N 2 ?ξ 1 I 3 …? ?N m ?ξ 1 I 3? ?N 1 ?ξ 2 I 3? ?N 2 ?ξ 2 I 3 …? ?N m ?ξ 2 I 3? ?N 1 ?ξ 3 I 3? ?N 2 ?ξ 3 I 3 …? ?N m ?ξ 3 I 3? ??T? ?（7）

式中：? ε 為應(yīng)變張量; ε ?L、 ε ?N為線性應(yīng)變張量部分和非線性應(yīng)變張量部分;N i為形函數(shù);? a ??e為單元節(jié)點(diǎn)位移矢量;? B ??L、? B ??N為線性應(yīng)變張量部分和非線性應(yīng)變張量部分分別對? a ??e的轉(zhuǎn)換矩陣; θ 為位移梯度矢量; G 為 θ 對? a ??e 的轉(zhuǎn)換矩陣。

W e=? n i? μ n α n （λ? a n? 1+λ? a n? 2+λ? a n? 3-3） （8）

式中： λ 1、λ 2、λ 3為伸長比;μ n、α n 為參數(shù)。

a∫? V 0? B? T D ?B? ?～ ?d V=∫? V 0?? N ???T p 0 d V+∫? A? 0t??? N ???T q 0 d A （9）

式中：? D 為柔度矩陣;p 0為體力;q 0為面力;V 0為積分體;A? 0t? 為積分面。

變形接觸非線性法向不嵌入，切向采用罰函數(shù)。

2 水密性仿真計算

如前所述，自閉式水封無法調(diào)節(jié)背部充壓（因背部充壓就是庫水壓力）從而調(diào)整水封變形形態(tài)，因此在水封設(shè)計時就必須對其在使用過程中的變形形態(tài)有著相當(dāng)準(zhǔn)確的預(yù)計。這種預(yù)計由仿真計算來實(shí)現(xiàn)，關(guān)注的計算指標(biāo)是封頭的接觸寬度和接觸應(yīng)力（肢體與背板和壓板的接觸應(yīng)力也要保證不漏水）? [15-17] 。封頭在庫水壓力背推作用下外伸，與止水面板的接觸應(yīng)力大于庫水壓力作為自閉封水條件（否則封頭會被庫壓水沖開從而脫離止水面板而漏水），自閉完全依賴于橡膠水封的變形。要實(shí)現(xiàn)期望的水封變形，需要通過大量的試算找出水封最合適的斷面，這一過程便是自閉式水封設(shè)計。

以某一尖頭“山”形水封為例（見圖2）來說明自閉式水封的仿真計算，圖2斷面是通過多次試算后選擇的。通過材料實(shí)驗(yàn)測定的橡膠材料 Ogden 應(yīng)變能函數(shù)模型參數(shù)見表1? [13] ，橡膠和鋼板的接觸，法向取為非嵌入，切向采用罰函數(shù)，罰數(shù)取0.5? [18-19] 。通過大變形有限元仿真計算，得到水密性仿真計算結(jié)果見表2。

由表2（右側(cè)為計算封水區(qū)）可知各種工況下的止水水頭。當(dāng)止水面板與封頭的初始間隙量為2 mm時（即預(yù)壓-2mm），庫水壓力為3～70 m水頭時，封頭與止水面板的接觸應(yīng)力均大于庫水壓力，水封能止住水;而止水面板的預(yù)壓量為0時，止水水頭可以提高到120 m水頭;其他預(yù)壓量下的該型水封能封住工程所有工況下的水頭。分析封水規(guī)律可知該水封當(dāng)初始間隙量為2 mm、預(yù)壓量為0時均出現(xiàn)高壓時封頭與止水面板分離的現(xiàn)象（接觸應(yīng)力為0），說明這2種情況只能應(yīng)用于較低水頭的工程;而預(yù)壓量不小于3 mm時，封頭持續(xù)封水保壓，接觸應(yīng)力與背壓（庫壓等于背壓）基本成線性關(guān)系。在保證接觸應(yīng)力滿足封水的要求下，封頭和止水面板有一定的接觸寬度是封水的更佳保障。經(jīng)計算可知：在封水區(qū)段壓力越大接觸寬度也越寬，高壓時自閉式水封封水具有一定的保障，封頭接觸寬度大，封頭與止水面板不易脫開。在最大壓力工況下（庫壓和背壓為300 m水頭，止水面板預(yù)壓量6 mm），水封橡膠材料第一主應(yīng)力最大值為4.154 MPa（見圖3），小于水封橡膠材料強(qiáng)度（7 MPa），水封材料不會因開裂而漏水。

3 自閉式水封水密性試驗(yàn)

為了驗(yàn)證自閉式水封的封水效果，進(jìn)行了自閉式水封的封水試驗(yàn)。水封斷面變形特性在專用大型水封試驗(yàn)臺上進(jìn)行（見圖4）。試驗(yàn)裝置主要由高壓水箱、承壓面板（活動臺板）、加壓機(jī)、穩(wěn)壓器和各種控制、測量儀器儀表組成。承壓面板

（1 600 mm×1 600 mm）用以模擬閘門止水座板，可保持水封與止水面板（止水對偶面）之間所需的不同間隙;高壓水箱結(jié)構(gòu)形狀為長圓形，用來模擬庫水壓力，它由一臺高壓水泵和一臺穩(wěn)壓器及相應(yīng)儀表控制，一連通管向水封底座用于充壓使水封封頭外伸，另一連通管向封頭側(cè)部用于模擬庫水壓力使封頭與止水面板脫離。水封試件采取框形試件，以保障庫水封閉。

給水封施加壓力，底部充壓和側(cè)部施壓同時進(jìn)行并保持一致，觀察止水面板的不漏水區(qū)間。如果狀態(tài)從射流過渡到僅有水珠漏出到?jīng)]有水珠漏出，則進(jìn)入到封水區(qū)間;如果狀態(tài)從干燥過渡到有水珠漏出到射流，則為離開封水區(qū)間。記錄該型水封在不同初始間隙量下的自封閉壓力區(qū)間（見圖5），并與仿真計算結(jié)果進(jìn)行比較。

將試驗(yàn)結(jié)果與仿真計算結(jié)果進(jìn)行比較可知，試驗(yàn)結(jié)果的自封閉壓力區(qū)間略微小于仿真計算結(jié)果，封水規(guī)律體現(xiàn)出一致性。當(dāng)存在初始間隙量時封水高壓上不去（超過43 m水頭就封不住）、封水低壓下不來（低于13 m水頭也封不住）;當(dāng)封頭改為預(yù)壓時，則該型水封封水性能良好，試驗(yàn)中能封住的最高水頭為2.2 MPa。

4 結(jié) 語

自閉式水封或?qū)⒁云洫?dú)特的應(yīng)用優(yōu)勢取代傳統(tǒng)機(jī)械充壓伸縮式水封。通過準(zhǔn)確的仿真計算，闡述了“山”形自閉式水封的封水原理、封水規(guī)律和封水有效水頭范圍（水封肢體的水密性不再贅述），并通過試驗(yàn)驗(yàn)證了該水封水密性性能的可靠性。

該水封使整體水封裝置更加簡單，理論上提高了整體裝置的耐久性，但水封自閉后的保壓性、耐久性還需要進(jìn)一步研究和確認(rèn)。由于該水封封水對水頭有一定的要求，故該水封不適宜于太低水頭的庫水封水。

參考文獻(xiàn)：

[1]? ?楊斌，干量，劉禮華，等.三峽升船機(jī)閘門氣壓伸縮水封水密特性仿真設(shè)計[J].河海大學(xué)學(xué)報（自然科學(xué)版），2018，46（6）：558-564.

[2] ALAN N G.Engineering with rubber-how to design rubber component[M].Munich：Carl Hanser Verlag，2001.

[3] MAYAW P.Micromechanical analysis of fatigue cracks in cord-rubber composites[J].Composite Structures，2001，54（2）：459-465.

[4] 劉禮華，王蒂，歐珠光，等.橡膠類止水材料超彈性性能的研究與應(yīng)用[J].四川大學(xué)學(xué)報（工程科學(xué)版），2011，43（4）：199-204.

[5] 雷艷，劉禮華，陳五一，等.大截面承載傳力水封的非線性計算研究[J].水力發(fā)電學(xué)報，2010，29（2）：228-233.

[6] 付明偉，楊斌，劉禮華，等.拍式浮體閘門充壓水封仿真計算研究[J].人民長江，2018，49（15）：80-84.

[7] 劉禮華，王蒂，歐珠光，等.基于Ogden公式的充壓伸縮式水封非線性計算研究[J].固體力學(xué)學(xué)報，2011，32（1）：102-107.

[8] NAGENDRAV， AKSHANTALA R T.A mechanistic based model for predicting fatigue life of composite laminates[J].Materials Sciences and Engineering，2000，285（1）：303-313.

[9] BOYCE M C，ARRUD E M. Constitutive models of rubber elasticity：a review[J].Rubber Chemistry and Technology，2000，73（3）：511-516.

[10]? 匡震邦.非線性連續(xù)介質(zhì)力學(xué)[M].西安：西安交通大學(xué)出版社，1989.

[11] YANG B，LIU L H，ZHANG S K.Fitting theory and method of hyper-elastic material’s nonlinear constitutive relations base on BP neural network[J].Applied Mechanics and Materials，2012，137（1）：36-41.

[12] 楊斌.橡膠水封的變形計算理論與應(yīng)用研究[D].武漢：武漢大學(xué)，2014.

[13] 王葉青.典型橡膠彈性元件的性能分析及其實(shí)驗(yàn)工裝設(shè)計[D].株洲：湖南工業(yè)大學(xué)，2012.

[14] 王者超，李術(shù)才，喬麗蘋，等.大型地下石油洞庫自然水封性應(yīng)力–滲流耦合分析[J].巖土工程學(xué)報，2013，35（8）：1535-1543.

[15] 薛小香，吳一紅，李自沖，等.高水頭平面閘門P型水封變形特性及止水性能研究[J].水力發(fā)電學(xué)報，2012，31（1）：56-61，92.

[16] 趙海峰，梁小飛，楊維平.充壓伸縮式水封系統(tǒng)在錦屏一級大壩工程的應(yīng)用[J].水電與抽水蓄能，2016，2（2）：65-67，82.

[17] 鄭宏偉，張燎軍，關(guān)超年，等.高水頭平面鋼閘門水封接觸性態(tài)研究[J].水電能源科學(xué)，2014，32（1）：172-175.

[18] 李自沖，薛小香，吳一紅，等.高水頭平面閘門P型水封非線性有限元模擬[J].水利水電技術(shù)，2012，43（4）：33-36.

[19] 張廣川.阿爾塔什水利樞紐工程2號深孔泄洪洞鏈輪閘門水封型式設(shè)計研究[J].水電站機(jī)電技術(shù)，2016，39（10）：41-43.

（編輯：胡旭東）

Simulation calculation and experimental study on water tightness ?characteristics of mountain-like self-closing water seal

YANG Bin? 1，2 ，ZHOU Yuwei 2，WANG Di 3

（ 1.School of Civil & Architectural Engineering，Wuhan University，Wuhan 430072，China; 2.Hubei Water Resources Technical College，Wuhan430070，China; 3.Changjiang Survey，Planning，Design and Research，Co.，Ltd.，Wuhan 430010，，China ）

Abstract：

In order to overcome the defect of the gate structures complexity and the non-durability caused by the traditional pressurization method of the pressure-filled telescopic water seal，a self-closing water seal that utilizes the large deformation of the water seal under reservoir water was designed，thereby eliminating mechanical pressure.The key of this kind of water seal design is to accurately predict the deformation of water seal under various working conditions.Based on the non-linear tensor constitutive relation of the water-sealed rubber material，the non-linear characters of the geometric large deformation and the non-linear characteristic of contact，taking a sharp mountain-like water seal as an example，the self-closing simulation calculation of this kind of water seal was carried out by using continuous medium finite element method，which was also verified by indoor experiments.Simulation and test results showed that the mountain-like rubber water seal can often achieve better self-closure，greatly reduce the construction cost of gate pressure water seal device，facilitate engineering maintenance and improve durability，which has great application value.

Key words：

self-closing water seal;water seal deformation;contact stress at seal head;simulation calculation;continuum mechanics