基于有限元強度折減法重力壩整體穩定及安全評價

池明陽 向文錦

`(1.中國電建集團貴陽勘測設計研究院有限公司，貴州 貴陽 550081； 2.中國水利水電第七工程局有限公司，四川 成都 610081)

?

·水利工程·

基于有限元強度折減法重力壩整體穩定及安全評價

池明陽1向文錦2

`(1.中國電建集團貴陽勘測設計研究院有限公司，貴州 貴陽 550081； 2.中國水利水電第七工程局有限公司，四川 成都 610081)

結合某水庫復雜地質條件下重力壩工程實況，建立了壩體壩基系統三維有限元網格模型，并采用有限元強度折減法，基于位移突變性判據，分析了壩體—壩基系統應力變形規律，求得壩體—壩基整體穩定安全系數達到3.0以上，滿足規范要求。

重力壩，壩基，有限單元法，安全系數

1 工程概況

由于工程區域地質條件的復雜性，為充分考慮壩體—壩基系統協調變形作用下大壩應力變形特點及其整體穩定性，有必要開展壩體壩基系統三維有限元計算。

2 基本原理及方法

2.1 有限元強度折減法

有限元強度折減法的基本原理是在外荷載保持不變的情況下，將材料的抗剪強度指標c，φ按式(1)形式進行折減，用折減后的抗剪強度指標cf，φf取代原c，φ進行計算。

(1)

通過不斷折減材料的強度指標，進行有限元分析，直至邊坡達到臨界破壞，此時得到的折減系數f即為邊坡的安全系數K。

2.2 評價標準

采用有限元強度折減法推求壩基整體安全系數的過程中，需要一個判斷標準：什么時候壩基達到極限狀態，目前常采用突變理論作為失穩判據。

采用有限元強度折減法計算壩基穩定性時，隨著強度折減系數的增大，其安全系數逐漸減小，當壩基失去穩定時，滑體部分將會產生很大的位移，滑體部分將由穩定狀態轉變為運動狀態，其位移和塑性應變不再是定值，而是處于無限塑性流動狀態，滑動巖/土體的位移會有突然變化的現象，基于此提出了突變性判據。

3 有限元計算與分析

3.1 有限元網格模型

根據相關設計資料建立壩體壩基系統三維有限元網格模型，如圖2所示。

3.2 基本計算參數

巖體物理力學參數見表1，表2。

表1 巖體物理力學參數地質建議值

表2 結構面物理力學參數地質建議值

3.3 應力變形分析

1)應力。計算結果表明：a.大壩澆筑后，第一主應力最大值約為1.5 MPa，處于受壓狀態，第三主應力最大值約為5.3 MPa，均出現在模型底部；b.由于巖性較復雜，在巖層屬性變化處，應力存在不連續現象；c.蓄水后，大壩整體仍表現為受壓狀態，僅壩踵區域存在少量的拉應力區，最大拉應力約為0.52 MPa；d.壩體大多處于1 MPa左右的壓應力狀態。

2)變形。計算結果表明：a.澆筑壩體后，整體向下沉降，向上游傾斜，最大沉降量約為4 mm，最大傾向上游量為2.7 mm，均發生在河床中心剖面的壩頂位置；b.蓄水后，壩體表現為向下游變形，最大變形量約為2.82 mm，同樣也發生在河床中心剖面的壩頂位置；c.由于體型基本對稱，蓄水后，大壩兩側均向河床中間變形，變形量約為1 mm。

3.4 安全評價

統計各壩段壩頂特征點順河向位移uy隨折減系數f變化曲線，如圖3所示。基于突變理論得出，折減系數大于3.0后，曲線出現突變，因此判斷，壩體壩基系統整體安全系數大于3.0。

4 結語

根據工程設計資料，建立壩體—壩基系統整體三維有限元網格模型，基于彈塑性理論對其進行了計算分析，得出了壩體—壩基系統應力變形規律；采用有限元強度折減法對壩體—壩基系統整體安全系數進行了分析，結合壩體—壩基屈服區分布情況，分析壩體—壩基整體系統安全系數大于3.0，大壩整體處于穩定狀態。

[1] 陳勝宏.高壩復雜地基及巖石高邊坡穩定分析[M].北京：中國水利水電出版社，2001.

[2] 李 征，張 林，陳 媛，等.復雜地基上重力壩的破壞過程及穩定安全度研究[J].中國農村水利水電，2016(2)：13-15.

[3] 王 河，王志鵬，李永剛.基于有限元超載—折減綜合法的重力壩深層抗滑穩定性分析[J].水電能源科學，2016(2)：8.

[4] 錢聲源，杜 鵑，張 冬，等.基于綜合判據的重力壩深層抗滑穩定性分析[J].水電能源科學，2015(12)：21-23.

[5] 李維樹，張睿琳.建基巖體摩擦系數對重力壩抗滑穩定影響分析[J].地下空間與工程學報，2015(A2)：14-15.

[6] 榮 夢，常曉林，周 偉，等.基于霍克-布朗(Hoek-Brown)準則重力壩抗滑穩定分析[J].武漢大學學報(工學版)，2015(4)：29-30.

[7] 楊志勇，陳建康，吳震宇，等.重力壩深層抗滑穩定計算方法探討[J].四川大學學報(工程科學版)，2015(A1)：19-21.

Gravity overall stability and safety evaluation based on finite element strength reduction

Chi Mingyang1Xiang Wenjin2

(1.GuiyangEngineeringCorporationLimited，PowerChina，Guiyang550081，China； 2.SinohydroBureau7Co.,Ltd，Chengdu610081,China)

It established three-dimensional finite element mesh model combining a reservoir which is a gravity dam under complicated geological conditions, and analyzed the dam body-the dam foundation system’s stress deformation law using the finite element strength reduction basing on the displacement mutant criterion and calculated the safety factors of overall stability in the dam body-the dam foundation system which over 3.0 and them met regulatory requirements.

gravity， dam foundation， FEM， safety factors

1009-6825(2016)26-0209-02

2016-07-03

池明陽(1980- )，男，碩士，高級工程師

TV642.3

A