基于AutoBANK的丁塢水庫壩坡穩定性分析

劉永,馬飛,張坤強,張風友

基于AutoBANK的丁塢水庫壩坡穩定性分析

劉永1,馬飛1,張坤強2*,張風友3

1. 德州市水利事業發展中心, 山東 德州 253014 2. 山東農業大學水利土木工程學院, 山東 泰安 271018 3. 青州市規劃編制研究中心, 山東 青州 262500

水庫壩坡的穩定直接關系到水庫樞紐工程的安全性，本文以丁塢水庫為工程實例，采用先進的有限元分析技術和全自動網格劃分手段，利用AutoBANK對該水庫壩坡穩定進行分析。這一技術手段為今后工程的模擬提供了一定思路。

水庫;壩坡; 穩定性

1 工程概況

樂陵市位于魯北平原，德州市東北部，全市總面1172 km2，轄9鎮3鄉4個辦事處，現狀人口69萬人，耕地面積6.78萬hm2。樂陵市當地水資源嚴重缺乏，根據可供水量和需水量分析計算，保證率50%時現狀年、2020年、2030年分別缺水1764.21萬m3、1289.98萬m3、242.44萬m3，缺水率分別6.77%、4.56%、0.84%，已經成為經濟社會發展的“瓶頸”。因而，興建平原水庫，利用非農田灌溉期蓄水備用，解決水資源緊缺問題已是當務之急。

丁塢水庫不僅是樂陵市城鄉供水的骨干構架，也是樂北引黃調度的樞紐。水庫位于樂陵市丁塢鎮西部，馬家溝以南，趙橋溝與蛤蟆劉溝之間，水庫近似呈瓶子狀，總庫容901.5萬m3，水庫總占地面積130.38 hm2。建設后可有效發揮引黃效益，保證城鄉居民用水安全。

2 壩坡穩定分析計算

2.1 計算程序

AutoBANK軟件是針對我國的水利水電行業要求而設計的有限元計算法計算程序，利用圖形化的界面，對水工建筑物的滲流、變形、應力、穩定等進行一體化計算，而且具有嚴密的計算過程。可進行位移、應力應變分析，并根據輸出的位移、應力、變形等數據，得到主應力隨空間和時間的變化規律，并自助生成任意點、任意斷面的內力數值和力學參數分布圖。

2.2 水庫圍壩設計斷面

操作步驟：啟動AutoBANK軟件——工程里面的材料屬性定義——繪制單元格，同時選擇各單元的材料數據——添加模型邊界條件，荷載——保存AutoCAD圖形文件，形成數據文件——求解——最后查看各項計算結果。

圍壩全長4.69 km，為碾壓式均質壩。壩頂寬7.0 m，上、下游壩坡均為1:3，壩頂高程17.45 m，防浪墻頂高程18.05 m，壩高7.0 m，在上游壩腳處設混凝土防滑齒墻，深1.0 m，寬0.6 m。水庫近似呈矩形，水庫占地130.38 hm2。設計蓄水位為16.30 m，相應庫容為901.5萬m3，死水位7.95 m，死庫容134萬m3，庫底高程6.15 m。調節庫容767.5萬m3。

庫區原地面平均高程10.45 m，下挖4.3 m，向上筑7.0 m，挖出的土方部分用于筑壩，多余土方運出另做它用。壩址地形稍有起伏，地面高程變化不大，根據壩址地質剖面圖，選用標準斷面作計算斷面。壩基為砂壤土和壤土，下部以砂壤土為主，計算中，對壩基各土層計算指標進行綜合選取。

2.3 穩定計算

按施工期，穩定滲流期，水庫水位降落期三種工況，有防滲塑膜計算。穩定滲流期，上游水位16.04 m進行上游壩坡穩定計算；施工期，圍壩剛建成，水庫未蓄水時進行上、下游壩坡穩定計算；水位驟降期，庫水位由正常蓄水位16.04驟降到死水位7.95 m，計算上游壩坡穩定。計算條件見表1。

表 1 壩坡穩定分析表

壩坡穩定計算采用瑞典法、畢肖普法，利用AutoBANK邊坡穩定分析程序進行計算，其中各種計算工況抗剪強度指標的設計采用值及飽和容重、粘聚力、內摩擦角等參數參照《工程地質、水文地質勘察報告》。

式中：—滑弧半徑；∑—某一土條的總高度；—土條底面中點處切線與水平線的夾角；—壩坡外水位高出土條底面中點的距離；o—水容重；U—在條塊底面中點處的孔隙水壓力值；、—土的強度指標(根據所采用的計算方法確定)；∑W(∑W)—計算抗滑力時單位寬度土條的重量。

邊坡穩定采用的條件按表1選定。土的抗剪強度指標的測定和應用見表2，物理力學指標按地址勘探提供的土工試驗成果建議值選取，壩體及壩基指標見表3。

表 2 土的抗剪強度指標的測定和應用

表 3 壩體及壩基指標

2.4 結果分析

分別取丁塢水庫東、西、南、北四個方向的斷面作為計算標準斷面。計算成果如下：

東壩段標準斷面計算成果如圖1所示：施工期，有效應力法，無滲流，瑞典法，0 g，2.27538(左)，2.25542(右)。

正常運行期，有效應力法，(水位線)設計水位，瑞典法，0 g，3.89826(左)，2.05903(右)。

水位降落期，有效應力法，(水位線)設計水位～(水位線)施工期水位線，瑞典法，0 g，2.1853(左)，2.25063(右)。

圖 1 東壩段標準斷面計算成果圖

西壩段標準斷面計算成果如圖2所示：

圖 2 西壩段標準斷面計算成果圖

施工期，有效應力法，無滲流，瑞典法，0 g，2.07116(左)，2.07047(右)。

正常運行期，有效應力法，(水位線)設計水位，瑞典法，0 g，3.6387(左)，2.47564(右)。

水位降落期，有效應力法，(水位線)設計水位～(水位線)施工期水位線，瑞典法，0 g，2.25485(左)，2.40549(右)。

北壩段標準斷面計算成果圖如圖3所示：

圖 3 北壩段標準斷面計算成果圖

施工期，有效應力法，無滲流，瑞典法，0 g，2.13311(左)，2.11629(右)。

正常運行期，有效應力法，(水位線)設計水位，瑞典法，0 g，3.5203(左)，2.26911(右)。

水位降落期，有效應力法，(水位線)設計水位～(水位線)施工期水位線，瑞典法，0 g，2.09153(左)，2.20267(右)。

南壩段標準斷面計算成果如圖4所示：

圖 4 南壩段標準斷面計算成果圖

施工期，有效應力法，無滲流，瑞典法，0 g，1.97283(左)，1.95418(右)。

正常運行期，有效應力法，(水位線)設計水位，瑞典法，0 g，3.62228(左)，2.27189(右)。

水位降落期，有效應力法，(水位線)設計水位～(水位線)施工期水位線，瑞典法，0 g，2.13805(左)，2.23177(右)。

庫區壩坡穩定計算成果見表4。

表 4 庫區壩坡穩定計算成果表

滲流穩定期所得安全穩定系數為無考慮邊坡外側靜水壓力的結果。完建期和滲流穩定期采用的計算方法均為有效應力法。由以上計算結果可以看出，迎水坡、背水坡壩坡在正常運行和非常運行工況下計算的抗滑穩定系數均大于《碾壓式土石壩設計規范》正常運行和非常運行允許抗滑穩定系數，滿足規范要求。

3 結 語

通過以上分析可以看出：相比以往的傳統計算方法，AutoBANK有限元分析軟件能較好地模擬壩坡斷面，且結果更具可靠性，能較準確反映實際工程情況。計算結果詳細精確，可為大壩的設計提供可靠的依據，并可輸出精確直觀的圖形文件成果。軟件與AutoCAD相配合使用，還能使計算操作更加便捷。采用有限元進行分析不僅提高了水利水電工程技術人員關于大壩穩定計算的效率和質量，而且能夠為設計提供詳實的基礎數據，為今后的工程模擬提供了一定思路。

[1] 徐潞珩,陳瀚,黃冠霖.基于Autobank的邊坡坡角對邊坡穩定性影響研究[J].科技創新與應用,2013(30):236-237

[2] 杜守來.AutoBank軟件在土壩滲流穩定計算中的應用[J].現代農業科技,2012(6):252-253

[3] 馬飛,馬兆會,李永春.基于ANSYS的孟家閘閘室穩定有限元分析[J].水利規劃與設計,2013(11):54-55

[4] 朱峰,湯洪潔.提防工程設計若干問題分析[J].水利規劃與設計,2013(6):55-57

[5] 陳璽,孫平,李守義,等.考慮地震與潰壩洪水共同作用的土石壩壩坡穩定分析方法[J].水利學報,2017,48(12):1499-1506

[6] 林淏圣.Autobank軟件邊坡穩定計算存在的問題與解決方法[J].水利科技,2018(3):42-44

[7] 占平.基于AUTOBANK的土石壩穩定性研究[J].水利技術監督,2019(1):36-39

[8] 田國民,尹望軍,趙佳作.桃源水電站庫區岸坡穩定分析及運行現狀調查[J].大壩與安全,2020(5):40-43

[9] 肖秀麗.Autobank軟件在軟基邊坡穩定計算中深滑弧問題的探討[J].福建水利發電,2020(1):12-13

[10] 楊清青.水位降落對小型土石壩上游壩坡穩定的影響[J].工程技術研究,2020,5(12):275-276

Analysis on the Stability of Dingwu Reservoir Dam Slope Based on AutoBANK

LIU Yong1, MA Fei1, ZHANG Kun-qiang2*, ZHANG Feng-you3

1.253014,2.271018,3.262500,

The stability of the dam slope is directly related to the safety of the reservoir pivot project. This paper takes Dingwu Reservoir as an engineering example, adopts advanced finite element analysis technology and automatic grid division method, and uses AutoBANK to analyze the dam slope stability of the reservoir. This technical means provides some ideas for future engineering simulation.

Reservoir; dam slope; stability

TV62

B

1000-2324(2022)01-0145-05

10.3969/j.issn.1000-2324.2022.01.022

2021-08-02

2021-09-23

山東省自然基金(ZR2018PEE027)

劉永(1986-),女,高級工程師,研究方向:水利工程管理. E-mail:261179971@qq.com

Author for corrspondence. E-mail:zkq0129@163.com