于 浩

(江門市科禹水利規(guī)劃設(shè)計咨詢有限公司，廣東 江門 529000)

1 大型水閘結(jié)構(gòu)靜動力分析

1.1 水閘結(jié)構(gòu)分析概述

土基上開敞式水閘閘室底板的應(yīng)力分析可采用反力直線分布法或彈性地基梁法。巖基上開敞式水閘閘室底板的應(yīng)力分析可按基床系數(shù)法計算。

開敞式水閘閘室底板的應(yīng)力可按閘門門檻的上下游分別進行計算，并計入閘門門檻切口處分配于閘墩和底板的不平衡剪力。開敞式或胸墻與閘墩簡支連接的胸墻式水閘，其閘墩應(yīng)力分析方法應(yīng)根據(jù)閘門形式確定。平面閘門閘墩應(yīng)力分析可采用材料力學方法，弧形閘門閘墩的應(yīng)力分析宜采用彈性力學方法。涵洞式、雙層式或胸墻與閘墩固支連接的胸墻式水閘，其閘室結(jié)構(gòu)應(yīng)力可按彈性地基上的整體框架結(jié)構(gòu)進行計算[1]。

受力條件復雜的大型水閘閘室結(jié)構(gòu)宜視為整體結(jié)構(gòu)采用空間有限元法進行應(yīng)力分析，必要時應(yīng)經(jīng)結(jié)構(gòu)模型試驗驗證。

1.2 水閘結(jié)構(gòu)應(yīng)力分析方法比較

水閘底板應(yīng)力分析中的反力直線分布法和彈性地基梁法，把水閘的空間結(jié)構(gòu)簡化為平面問題，采用“截板成梁”的方法進行計算，這種計算方法忽略了各閘墩處變位不等的重要因素，誤差較大，因而不宜在大中型水閘設(shè)計中采用[2]。

平面閘門閘墩的受力條件主要是偏心受壓，可假定閘墩固定于底板上的懸臂梁，其應(yīng)力狀況可采用材料力學的方法進行分析，弧形閘門閘墩的受力條件比較復雜，不只是偏心受拉，而且還受扭，是一塊一邊固定、三邊自由的彈性矩形板，其應(yīng)力狀況宜采用彈性力學的方法進行分析。因此，對于大型水閘弧形閘門閘墩，宜按閘門支座作用力的實際位置，加大閘門支座處的閘墩厚度并考慮閘底板的約束作用等邊界條件采用有限單元法進行應(yīng)力分析。

1.3 靜力學分析平衡方程

靜力分析適用于緩慢加載的情況，如果完成加載時間遠大于結(jié)構(gòu)的被激活振動周期，則認為是緩慢加載，此種情況結(jié)構(gòu)的響應(yīng)慢，不會產(chǎn)生震蕩，慣性力可以忽略不計。且遵循如下平衡方程式(1)：

[K]{u}={F}

(1)

式中：[K]為剛度矩陣；{u}為位移矩陣；{F}為荷載列矩陣。

1.4 動力學分析平衡方程

動力學分析適用于快速加載、沖擊碰撞的情況，這種情況下慣性力和阻尼的影響不能忽略。動力學分析需要考慮結(jié)構(gòu)的慣性，所以材料參數(shù)必須定義密度及材料的彈性模量、泊松比。動力學問題遵循如下平衡方程式(2)：

[M]{ü}+[C]{ü}+[K]{u}={F(t)}

(2)

式中：[M]為質(zhì)量矩陣；[M]{ü}反映了慣性效應(yīng)；[C]為阻尼矩陣，[C]{ü}反映了阻尼效應(yīng)，t為時間列矩陣；F(t)為行載時間函數(shù)矩陣。

2 某大型水閘閘室體型結(jié)構(gòu)

某水利樞紐由土石壩、泄洪閘、河床式水電站等建筑物組成。泄洪閘采用平底寬頂堰，堰體建基面高程1062.3 m(最低的齒槽處高程為1060.3 m)，閘頂前部高程1078.7 m、后部高程1079.8 m，順水流方向長40.0 m。

泄洪閘閘室布置采用兩孔一聯(lián)，單孔凈寬14.0 m，中墩厚3.0 m，縫墩厚2.5 m，每聯(lián)長度36.0 m，閘底板厚度2.7 m?；⌒伍l門前、后設(shè)檢修閘門。閘墩頂部設(shè)交通橋與壩頂公路相連?；A(chǔ)為深厚砂礫石層[3-7]。

2.1 閘室荷載及組合

設(shè)計荷載如下：結(jié)構(gòu)自重和永久設(shè)備重、水重；上、下游及縫間靜水壓力(水庫正常蓄水位：1076 m，正常尾水為：1065.02)；揚壓力；浪壓力；地震荷載。

基本荷載組合包括：完建工況、正常擋水工況。

特殊荷載組合包括：檢修工況、正常擋水+地震工況。

由于完建工況、正常擋水工況和檢修工況的結(jié)構(gòu)計算都屬于靜力分析且過程相似。本文只對正常擋水工況下水閘的靜力特性進行有限元分析。

2.2 荷載及模型的簡化處理

由于計算關(guān)心的是閘底板和閘墩的結(jié)構(gòu)應(yīng)力，可忽略交通橋、檢修橋支座處的局部應(yīng)力，為方便地建立有限元模型及施加荷載，有效地控制單元劃分數(shù)目，需要對模型及荷載進行簡化處理。

模型簡化處理：為方便單元劃分把閘墩上的臺階簡化為斜面，把弧形墩頭按矩形處理，不考慮交通橋、檢修橋的局部支座。縫面處的兩道止水合并按一道考慮，弧形門槽按直線近似處理，不考慮下游檢修門槽。

荷載的簡化處理：弧形閘門自重折算為線荷載作用在閘底板頂面；弧形閘門水推力根據(jù)劃分的網(wǎng)格節(jié)點折算為表面壓力和節(jié)點力施加在牛腿上；橋梁自重及相關(guān)荷載按作用在閘墩上的面荷載進行折算；油泵房的荷載折算為面荷載施加在中墩頂中部；對揚壓力做適當簡化，把底板底面沿水流方向分成14個面，每個表面施加表面壓力來模擬揚壓力。

基礎(chǔ)簡化：基礎(chǔ)為深厚砂礫石基礎(chǔ)，根據(jù)工程經(jīng)驗，上、下游方向以及深度方向取15 m范圍，垂直水流方向與閘室同寬。

坐標系設(shè)置：采用平面直角坐標系，閘底板頂面Y坐標值為0，其他高程相應(yīng)折減，尺寸單位按米計。

3 結(jié)果與分析

運行Autodesk Simulation Mechanical程序，執(zhí)行CAD實體模型數(shù)據(jù)交換，分析類型設(shè)置為線性靜力分析。采用8節(jié)點三維塊體單元進行網(wǎng)格剖分?；A(chǔ)部分視為彈性基礎(chǔ)，其自重引起的變形忽略不計，質(zhì)量密度不賦予數(shù)值；閘室材料屬性為C25混凝土，密度為2400 kg/m3，泊松比為0.167，彈性模量為28 GPa。基礎(chǔ)底面為固定約束，基礎(chǔ)上游面、下游面X向約束，基礎(chǔ)沿壩軸線方向的兩個面Z向約束，閘室底面和基礎(chǔ)的接觸采用黏合處理。

交通橋、檢修橋、油泵房相關(guān)荷載折算為面荷載施加在閘墩頂部，面荷載數(shù)值見表1。

表1 墩頂面荷載 kN·m-2

對底板底面的14個面域施加揚壓力，每個面域?qū)?yīng)的揚壓力值見表2，底板底面面域的劃分見圖1。

表2 底板底面每個區(qū)域揚壓力值 kN·m-2

圖1 閘室底板面域的劃分

弧形閘門的重力折算80 kN施加到底板對應(yīng)的節(jié)點上?；⌒伍l門的水推力傳遞到一個牛腿的推力為6309 kN，牛腿受力面積為3.64 m2，折算為表面壓力1733 kN/m2施加到牛腿表面。

由線性靜力分析結(jié)果可知，由于閘墩固結(jié)于底板，牛腿固結(jié)于閘墩，閘墩與底板的銜接處、閘墩與牛腿的連接處呈現(xiàn)應(yīng)力集中現(xiàn)象。但最大拉應(yīng)力、壓應(yīng)力均小于混凝土的強度運行范圍內(nèi)[8]。

將靜力分析的水閘數(shù)值模型分析類型設(shè)置為荷載剛化固有頻率分析，模態(tài)數(shù)量為10，固有頻率分析結(jié)果見表3。

表3 水閘固有頻率分析結(jié)果

將靜力分析的水閘數(shù)值模型分析類型設(shè)置為響應(yīng)譜分析，單元劃分、材料、約束條件、荷載等與靜力分析相同，受到Z向和Y向的地震荷載作用，地震加速度頻譜數(shù)據(jù)見表4。

表4 地震加速度響應(yīng)頻譜

由線性動力分析結(jié)果可知，增加水閘的階數(shù)，固有頻率隨之增加，水閘的自振特性于一般規(guī)律相符。水閘在Z和Y兩向荷載作用下，峰值位移及峰值應(yīng)力均為超過規(guī)范的允許值。

4 結(jié) 論

本文以某大型水閘實例為研究對象，通過Autodesk Simulation Mechanical三維有限元分析軟件，構(gòu)建三維數(shù)值模型，對水閘結(jié)構(gòu)的靜動力特性進行數(shù)值分析。以期了解結(jié)構(gòu)內(nèi)部的應(yīng)力分布狀態(tài)，為水閘的結(jié)構(gòu)設(shè)計提供合理、可靠的理論依據(jù)，對結(jié)構(gòu)應(yīng)力集中部位加大配筋或優(yōu)化結(jié)構(gòu)減少應(yīng)力集中。