劉竹麗 陳赟 伊元忠

摘 要：為研究閘門流激振動(dòng)問題，采用雙向流固耦合方法，以ANSYS-Workbench為工具，建立了閘門和流場(chǎng)的有限元模型，將閘門和流場(chǎng)的接觸面設(shè)置為流固耦合接觸面，進(jìn)行流場(chǎng)和結(jié)構(gòu)場(chǎng)間數(shù)據(jù)的雙向交換，對(duì)不同開度下閘門流激振動(dòng)情況進(jìn)行了數(shù)值模擬，得到了閘門不同開度下的位移、應(yīng)力情況。結(jié)果表明：閘門振動(dòng)位移及應(yīng)力隨閘門開度的增大而逐漸減小。

關(guān)鍵詞：平面閘門;流激振動(dòng);流固耦合;數(shù)值模擬

中圖分類號(hào)：TV663+.4 ? 文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2020.11.016

Abstract：In order to study the flow-induced vibration of the gate， the finite element model of the gate and flow field were established by ANSYS-Workbench. The contact surfaces between the gate and flow field were set as fluid structure interaction interface to exchange data between flow field and structure field. Based on two-way fluid-structure interaction， the research conducted numerical simulation of the flow-induced vibration of the gate with different openings and obtained the distribution of displacement and stress at each opening of the gate. The results show that the displacement and stress of the gate gradually decrease with the increase of the opening degree.

Key words： plane gate; flow-induced vibration; fluid structure interaction; numerical simulation

閘門是水利樞紐工程的重要組成部分，保證閘門正常運(yùn)行對(duì)水利工程至關(guān)重要。閘門在動(dòng)水作用下易發(fā)生振動(dòng)，劇烈振動(dòng)會(huì)影響閘門的安全運(yùn)行，甚至導(dǎo)致閘門破壞，給人民的生命財(cái)產(chǎn)安全帶來巨大威脅，閘門的流激振動(dòng)問題日漸受到重視。目前研究閘門振動(dòng)問題的主要方法有原型觀測(cè)法[1]、模型試驗(yàn)法[2-3]、數(shù)值分析法[4-6]。古華等[7]得到了流體對(duì)閘門振動(dòng)特性的影響規(guī)律及流固耦合計(jì)算時(shí)液相長(zhǎng)度的合理取值范圍。胡劍杰等[8]研究了影響水流作用下弧面三角閘門自振特性的因素，得到了閘門干模態(tài)以及在不同門前水體寬度和高度下的濕模態(tài)。薛惠芳[9]研究了流固耦合對(duì)閘門結(jié)構(gòu)動(dòng)態(tài)特性的影響，分析了平面閘門在無水和有水狀態(tài)下的振動(dòng)特性。目前，閘門振動(dòng)特性方面的研究較多，但關(guān)于流激振動(dòng)下閘門的應(yīng)力及位移情況的研究較少。

本文基于ANSYS Workbench，通過流固耦合方法對(duì)閘門不同開度下的流激振動(dòng)情況及開度對(duì)閘門振動(dòng)位移和應(yīng)力的影響進(jìn)行仿真分析，以期為閘門設(shè)計(jì)及閘門的流激振動(dòng)仿真分析提供參考。

1 流固耦合理論基礎(chǔ)

1.1 結(jié)構(gòu)動(dòng)力學(xué)方程

1.3 流固耦合理論

流固耦合的求解過程較為復(fù)雜，需要通過流固耦合交界面進(jìn)行數(shù)據(jù)交換，實(shí)現(xiàn)流場(chǎng)與結(jié)構(gòu)場(chǎng)的耦合模擬[10]。計(jì)算時(shí)，單獨(dú)建立流場(chǎng)和結(jié)構(gòu)場(chǎng)的控制方程，在一個(gè)時(shí)間步中分別求解，流場(chǎng)計(jì)算結(jié)果通過耦合面?zhèn)鬟f給結(jié)構(gòu)場(chǎng)，結(jié)構(gòu)場(chǎng)計(jì)算結(jié)果通過耦合面?zhèn)鬟f給流場(chǎng)。雙向流固耦合分別對(duì)流場(chǎng)和結(jié)構(gòu)場(chǎng)進(jìn)行求解，既可利用相互獨(dú)立的兩個(gè)領(lǐng)域內(nèi)各自算法的優(yōu)勢(shì)，計(jì)算簡(jiǎn)便，精度提高，又通過耦合面實(shí)現(xiàn)了不同計(jì)算域間的數(shù)據(jù)傳遞，考慮了流場(chǎng)和結(jié)構(gòu)場(chǎng)間的相互影響，計(jì)算結(jié)果準(zhǔn)確，符合工程實(shí)際。流場(chǎng)與結(jié)構(gòu)場(chǎng)間數(shù)據(jù)交換類型見表1。

2 模型描述

2.1 閘門模型

某水利樞紐泄洪閘事故檢修閘門為潛孔式平面鋼閘門，孔口尺寸為8.0 m×10.0 m（寬×高），閘底高程為1 750 m，正常蓄水位為1 820 m。平面閘門為焊接結(jié)構(gòu)，閘門門葉布置10道主梁、8道縱梁及4道邊梁，主輪布置在內(nèi)外邊梁之間。閘門材料為Q345C，材料屬性見表2。

為避免外部軟件所建模型導(dǎo)入ANSYS時(shí)因接口問題而造成扭曲、多面、丟面的現(xiàn)象，本文直接在ANSYS Workbench的DM模塊中建立閘門實(shí)體模型，而后抽取模型的中面，并根據(jù)實(shí)際情況進(jìn)行適當(dāng)修補(bǔ)，生成閘門面模型。根據(jù)閘門結(jié)構(gòu)特點(diǎn)，選取SHELL 181單元[11]對(duì)閘門中面模型進(jìn)行網(wǎng)格劃分，得到閘門有限元模型。SHELL 181單元為殼體單元，適合分析中等厚度的板殼結(jié)構(gòu)，具有大變形、大扭轉(zhuǎn)的功能。該單元有4個(gè)節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)有6個(gè)自由度。平面閘門有限元模型共有78 107個(gè)節(jié)點(diǎn)、80 422個(gè)單元。閘門跨度方向?yàn)閤向，閘門高度方向?yàn)閥向，上游至下游的流向方向?yàn)閦向，閘門兩側(cè)受到跨度方向和水流方向約束，閘門頂部受到豎直方向約束。

2.2 流場(chǎng)模型

在建立平面閘門模型的基礎(chǔ)上，在ANSYS Workbench中建立流場(chǎng)模型。上游庫(kù)區(qū)長(zhǎng)度為5倍的流道進(jìn)口高度，為節(jié)省計(jì)算資源，庫(kù)區(qū)只建立下半部分，流道長(zhǎng)160 m，見圖1。流場(chǎng)模型建立完成后，通過ICEM對(duì)模型進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分。庫(kù)區(qū)左側(cè)和庫(kù)區(qū)頂部為壓力入口：庫(kù)區(qū)左側(cè)為靜壓分布，具體壓力值通過UDF命令設(shè)置;庫(kù)區(qū)頂部設(shè)置為受30 m水深壓力。流道出口為大氣壓力出口。

2.3 流固耦合接觸設(shè)置

將閘門與水流接觸面以及水流與閘門接觸面分別在結(jié)構(gòu)場(chǎng)計(jì)算和流場(chǎng)計(jì)算中設(shè)置為流固耦合面，并在System Coupling中對(duì)這對(duì)流固耦合面設(shè)置數(shù)據(jù)交換。流場(chǎng)向結(jié)構(gòu)場(chǎng)傳遞力，結(jié)構(gòu)場(chǎng)向流場(chǎng)傳遞位移。計(jì)算時(shí)長(zhǎng)為2 s。

3 閘門計(jì)算結(jié)果分析

對(duì)相對(duì)開度為0.1～0.8的閘門流激振動(dòng)進(jìn)行模擬，以閘門相對(duì)開度為0.2的工況為例，分析流激振動(dòng)作用下閘門的位移及應(yīng)力分布情況。

3.1 流激振動(dòng)下閘門振動(dòng)位移

圖2為閘門相對(duì)開度為0.2時(shí)的位移云圖，閘門在流激振動(dòng)作用下最大位移出現(xiàn)在閘門面板跨度方向中軸線處，圖中Max位置。位移由閘門面板中部向左右兩側(cè)逐漸減小，且沿閘門面板中軸線對(duì)稱，位移由閘門中部向上部也逐漸減小，位移分布層次明顯。閘門變形發(fā)生在跨度方向中軸線附近，形式為彎曲。圖3為閘門最大位移剖面，由圖3可知，面板壁薄且受力面較大，沿順?biāo)鞣较虬枷菝黠@。閘門位移在初始階段于3.407～4.768 mm范圍內(nèi)波動(dòng)，計(jì)算0.4 s后穩(wěn)定在4.07 mm左右。

3.2 流激振動(dòng)下閘門振動(dòng)應(yīng)力

相對(duì)開度為0.2時(shí)，閘門在流激振動(dòng)作用下的應(yīng)力分布見圖4。由圖4可知，最大應(yīng)力出現(xiàn)在閘門主梁、縱梁與面板接觸位置，其余部位應(yīng)力相對(duì)較小，應(yīng)力分布基本沿閘門面板跨度方向的中軸線對(duì)稱。閘門應(yīng)力在初始階段于158.8～185.8 MPa范圍內(nèi)波動(dòng)，并在計(jì)算0.35 s后穩(wěn)定在172 MPa左右。

3.3 開度對(duì)閘門位移和應(yīng)力的影響

（1）位移情況。隨著閘門開度的增大，最大位移點(diǎn)逐漸向閘門面板下方移動(dòng)。本文采用隨機(jī)統(tǒng)計(jì)學(xué)方法分析動(dòng)位移，用位移均值+3倍均方差來估算各開度下的最大位移，閘門最大位移隨開度變化曲線見圖5。由圖5可知，閘門最大位移出現(xiàn)在相對(duì)開度為0.1時(shí)，為4.608 mm，閘門最大位移隨開度的增大逐漸減小，但相對(duì)開度為0.4時(shí)的最大位移僅次于相對(duì)開度為0.1時(shí)的，比相對(duì)開度為0.2、0.3時(shí)都要大。圖5 閘門最大位移隨相對(duì)開度變化曲線

相對(duì)開度為0.1～0.8時(shí)閘門位移隨時(shí)間變化曲線見圖6，閘門位移波動(dòng)范圍隨閘門開度的增大而減小，且達(dá)到穩(wěn)定后的位移也隨開度的增大而減小。在相對(duì)開度為0.4時(shí)，位移初始波動(dòng)幅度比相對(duì)開度為0.2、0.3時(shí)大，且閘門位移一直處于波動(dòng)的狀態(tài)。因此，延長(zhǎng)相對(duì)開度0.4時(shí)的計(jì)算時(shí)間至5 s。由圖7可看出，閘門開度為0.4時(shí)，位移后期在3.90～4.25 mm范圍波動(dòng)，直至計(jì)算5 s仍處于波動(dòng)狀態(tài)，估算得到此工況下最大位移為4.429 mm。

圖8為相對(duì)開度為0.2時(shí)主梁8在流激振動(dòng)下的位移云圖，最大位移出現(xiàn)在主梁跨中部位，兩側(cè)位移較小，位移由中部向兩側(cè)逐漸減小，且位移基本沿閘門面板中軸線對(duì)稱分布，位移朝向下游方向。

各個(gè)主梁在各開度下的最大位移均出現(xiàn)在跨中處，圖9為主梁1～10在不同開度時(shí)的最大位移曲線。由圖9可知，主梁1～5的最大位移都隨閘門開度的增大而逐漸減小，而主梁6～10在相對(duì)開度為0.4時(shí)的最大位移相比相對(duì)開度為0.3時(shí)明顯增大，開度大于0.4后，主梁最大位移隨開度的增大逐漸減小?？芍鄬?duì)開度為0.4時(shí)閘門的位移情況相對(duì)較復(fù)雜。

圖10為各開度下主梁最大位移變化曲線。由圖10可知，相對(duì)開度為0.1時(shí)，主梁位移最大。在相對(duì)開度為0.5～0.8時(shí)，主梁最大位移隨著主梁編號(hào)的增大而逐漸增大，開度越小，主梁位移越大。在相對(duì)開度為0.1～0.3時(shí)，主梁位移隨主梁編號(hào)的增大而增大，在主梁編號(hào)為8時(shí)出現(xiàn)小幅下降。相對(duì)開度為0.4時(shí)，主梁位移增大較快，主梁7～10的最大位移接近于相對(duì)開度為0.1時(shí)的。

（2）應(yīng)力情況。同樣采用隨機(jī)統(tǒng)計(jì)學(xué)方法分析動(dòng)應(yīng)力，用應(yīng)力均值+3倍均方差來估算各開度下的最大應(yīng)力，閘門最大應(yīng)力隨開度變化曲線見圖11。閘門最大應(yīng)力出現(xiàn)在相對(duì)開度為0.3時(shí)，為185 MPa。閘門最大應(yīng)力隨開度的增大而逐漸減小，這是因?yàn)殡S著閘門開度的增大，閘門與水體直接接觸面積減小，閘門垂直受力面積變小，閘門所受的水體作用也在減小，所以受力也逐漸減小。

閘門相對(duì)開度為0.1～0.8時(shí)應(yīng)力隨時(shí)間變化曲線見圖12。由圖12可知，閘門應(yīng)力的波動(dòng)幅度隨開度的增大而減小，且波動(dòng)一般只出現(xiàn)在初始階段，后期逐漸達(dá)到穩(wěn)定。但在相對(duì)開度為0.4時(shí)，閘門應(yīng)力一直處于波動(dòng)的狀態(tài)，其最大值也比相對(duì)開度為0.5時(shí)大。因此，延長(zhǎng)相對(duì)開度為0.4時(shí)的計(jì)算時(shí)間至5 s，由圖13可看出，閘門開度為0.4時(shí)，應(yīng)力在153～171 MPa間波動(dòng)，直至計(jì)算5 s仍處于波動(dòng)狀態(tài)，估算得到此工況下的最大應(yīng)力為174.32 MPa。結(jié)合圖7分析，相對(duì)開度為0.4時(shí)位移、應(yīng)力均處于波動(dòng)狀態(tài)，此工況是較危險(xiǎn)工況，在閉門時(shí)應(yīng)盡量快速通過此開度，避免較長(zhǎng)時(shí)間停留。

圖14為主梁1～10在不同開度時(shí)后翼緣跨中處應(yīng)力變化曲線。由圖14可知，主梁1～5的應(yīng)力皆隨閘門開度的增大而減小，而主梁6～10在相對(duì)開度為0.4時(shí)的最大應(yīng)力比相對(duì)開度為0.3時(shí)明顯增大，開度大于0.4后，主梁6～8應(yīng)力隨開度增大逐漸減小，主梁9、10在相對(duì)開度為0.6時(shí)最大應(yīng)力與相對(duì)開度為0.5時(shí)基本持平，之后應(yīng)力隨開度增大逐漸減小。

圖15為各開度下主梁跨中處應(yīng)力變化曲線。由圖15可知，在相對(duì)開度為0.6～0.8時(shí)，主梁應(yīng)力隨著主梁編號(hào)的增大而逐漸增大，且開度越大，主梁應(yīng)力越小。在相對(duì)開度為0.4～0.5時(shí)，主梁應(yīng)力隨編號(hào)增大而逐漸增大，到主梁9時(shí)出現(xiàn)小幅下降。在相對(duì)開度為0.1～0.4時(shí)主梁5、主梁8應(yīng)力比周圍主梁大，原因是主梁5和主梁8無前翼緣，容易產(chǎn)生較大應(yīng)力。相對(duì)開度為0.4時(shí)，主梁應(yīng)力增大較快，主梁6～10的最大應(yīng)力接近甚至超過了相對(duì)開度為0.1時(shí)的應(yīng)力。

4 結(jié) 論

通過流固耦合方法對(duì)閘門流激振動(dòng)進(jìn)行模擬，分別對(duì)不同開度下的流激振動(dòng)情況進(jìn)行研究，發(fā)現(xiàn)閘門在流激振動(dòng)作用下的位移和應(yīng)力與閘門開度存在一定關(guān)系，得到如下結(jié)論：

（1）各開度下最大位移都出現(xiàn)在閘門面板跨度方向的中軸線位置，且位移由閘門面板中部向左右兩側(cè)逐漸減小，位移沿閘門面板中軸線基本呈對(duì)稱狀態(tài)。面板壁薄且受力面大，沿水流方向會(huì)出現(xiàn)明顯凹陷。閘門最大位移隨著閘門開度的增大而逐漸減小，且出現(xiàn)最大位移的位置隨著閘門開度的增大逐漸向閘門下部移動(dòng)。

（2）各開度下閘門最大應(yīng)力均出現(xiàn)在閘門主梁、縱梁與面板接觸的位置，應(yīng)力基本沿閘門面板中軸線對(duì)稱。隨著閘門開度的增大，閘門垂直受力面積逐漸減小，最大應(yīng)力也逐漸減小。

（3）閘門最大位移為4.608 mm，與計(jì)算跨度的比值為1/1 953，遠(yuǎn)小于規(guī)范[12]規(guī)定的允許值（1/750）;最大應(yīng)力為185 MPa，在允許范圍內(nèi)，閘門強(qiáng)度、剛度滿足要求，能夠安全運(yùn)行。

參考文獻(xiàn)：

[1] 王桂青.陡河水庫(kù)輸水洞工作閘門振動(dòng)觀測(cè)分析[J].海河水利，2016（1）：48-50.

[2] 孔繁臣，王永潭，張喜峰，等.豐滿水電站重建工程潛孔弧形閘門流激振動(dòng)模型試驗(yàn)研究[J].水利水電技術(shù)，2016，47（6）：76-79.

[3] 吳杰芳，張林讓，陳敏中，等.三峽大壩導(dǎo)流底孔閘門流激振動(dòng)水彈性模型試驗(yàn)研究[J].長(zhǎng)江科學(xué)院院報(bào)，2001，18（5）：76-79.

[4] 劉權(quán)，張燎軍，關(guān)超年.高水頭平面閘門流激振動(dòng)的數(shù)值模擬[J].水電能源科學(xué)，2014，32（1）：176-179，190.

[5] 潘樹軍，王新.大型平面鋼閘門流激振動(dòng)模型試驗(yàn)與數(shù)值模擬[J].水電能源科學(xué)，2011，29（8）：148-151.

[6] 鮑秋慧，胡兆榮，孫蓬勃.升臥式臥倒閘門的自振特性[J].人民黃河，2013，35（2）：114-115，119.

[7] 古華，嚴(yán)根華.水工閘門流固耦合自振特性數(shù)值分析[J].振動(dòng)、測(cè)試與診斷，2008，28（3）：242-246.

[8] 胡劍杰，胡友安，王煦.基于ANSYS的弧面三角閘門自振特性研究[J].三峽大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，36（6）：24-27.

[9] 薛惠芳.基于ANSYS的平面閘門流固耦合振動(dòng)特性的有限元分析[J].安徽農(nóng)業(yè)科學(xué)，2012，40（8）：5007-5009.

[10] HSU M C，BAZILEVS Y. Fluid-Structure Interaction Modeling of Wind Turbines： Simulating the Full Machine[J]. Computational Mechanics， 2012， 50（6）： 821-833.

[11] 王新敏.ANSYS工程結(jié)構(gòu)數(shù)值分析[M].北京：人民交通出版社，2007：11.

[12] 中華人民共和國(guó)水利部.水利水電工程鋼閘門設(shè)計(jì)規(guī)范：SL74—2013[S].北京：中國(guó)水利水電出版社，2013：28.

【責(zé)任編輯 張華巖】