基于ANSYS 的某溢流雙曲拱壩應(yīng)力分析

黃德才 丁澤霖

（1.貴州省黔西南州望謨縣水利局，貴州望謨 552300;2.華北水利水電學(xué)院，河南鄭州 450011）

基于ANSYS 的某溢流雙曲拱壩應(yīng)力分析

黃德才1丁澤霖2

（1.貴州省黔西南州望謨縣水利局，貴州望謨 552300;2.華北水利水電學(xué)院，河南鄭州 450011）

基于三維有限元計算軟件ANSYS對某混凝土雙曲拱壩壩體應(yīng)力分析，通過分析得出，拱壩最大位移發(fā)生在壩體表孔位置，但是壩頂開設(shè)表孔對壩體應(yīng)力的影響是局部的，主要集中在孔口周邊，而對壩體下部幾乎沒有影響，并總結(jié)規(guī)律，為工程設(shè)計、施工和加固處理提供了重要科學(xué)依據(jù)。

水工結(jié)構(gòu)，應(yīng)力分析，有限元，拱壩

1 概述

拱壩由于其經(jīng)濟性、安全性、抗超載能力強等優(yōu)越性，在現(xiàn)代水利水電工程中越來越多的被采用，而且隨著現(xiàn)代工程經(jīng)驗的不斷積累和科技水平的進步，拱壩向著壩型更薄、壩體更高的方向發(fā)展［1］。由于拱壩主要是依靠左右壩肩山體的支撐來維持穩(wěn)定，拱壩壩身強度直接決定工程的正常運行和安全性。由于拱壩通常采用表、深孔聯(lián)合泄流的形式，破壞了拱圈的連續(xù)性，影響拱壩整體的應(yīng)力變形規(guī)律［2，3］，因此，為了保證高拱壩工程的穩(wěn)定和安全性，有必要針對高拱壩表孔對壩體應(yīng)力變形的影響問題進行研究，有限元法是解決上述問題的重要途徑之一［4］。

ANSYS是目前應(yīng)用最為廣泛的大型通用有限元軟件，其在大壩、橋梁、隧道及房屋建筑等土木工程中得到了最為廣泛的運用，是目前世界上唯一通過ISO 9001認證的有限元軟件［5］。因此，分析選用該軟件來對拱壩進行有限元靜力分析。本文以某高拱壩壩體應(yīng)力分析為研究對象，采用非線性有限元的ANSYS軟件分析表孔影響下拱壩壩體應(yīng)力分布情況，為工程設(shè)計、施工和加固處理提供了重要科學(xué)依據(jù)。

2 工程概述

某水電站攔河大壩為拋物線雙曲拱壩，壩頂表孔溢洪道，壩頂高程507.0 m，壩底高程417.0 m，最大壩高90 m，壩頂寬5 m，壩底寬17.5 m，采用壩頂溢流，溢流段布置在壩中部，堰頂高程500 m，溢流凈寬24 m，自由跌落水墊消能。拱壩體型如圖1，圖2所示。壩體材料為R200混凝土，容重2.4 t/m3，彈模2.8E6，泊松比0.167，線膨脹系數(shù)9×10－6/℃，導(dǎo)溫系數(shù)3m2/月，多年平均氣溫19.1℃。壩基巖體容重2.7 t/m3，彈模1.8E6，泊松比0.27，線膨脹系數(shù)10×10－6/℃，導(dǎo)溫系數(shù)3 m2/月。

圖1 整體三維有限元網(wǎng)格全視圖

3 有限元模型的建立

在有限元計算中，建模時作了一定的簡化，考慮了溢流堰的開口形式，模型的單元剖分在ANSYS軟件中進行，壩體采用Solid65單元，巖體采用Solid45單元。型兩岸邊界垂直河向約束，底面為三向固定約束，上、下游邊界順河向約束。計算中采用的直角坐標(biāo)為:X軸方向與水流方向垂直，以向左岸為正;Y軸方向平行水流方向，以向下游為正;高度方向為Z軸方向，向上為正，壩前、壩后、左右岸拱端及底部地基取1倍壩高左右。整體三維有限元網(wǎng)格圖、壩體有限元網(wǎng)格圖分別如圖1，圖2所示。

圖2 壩體三維有限元網(wǎng)格下游立視圖

4 計算結(jié)果分析

4.1 變形分析

壩體上部位移大于下部位移，表孔處位移大于拱端位移，拱冠上部位移大于下部位移，徑向位移大于切向位移。壩體位移總體上對稱性較好。壩體徑向變位整體趨向下游，其中，最大位移為1 cm，發(fā)生在溢流堰頂區(qū)域，如圖3所示。

圖3 下游壩面Y方向位移等值線圖

4.2 應(yīng)力分析

在孔口局部應(yīng)力集中現(xiàn)象十分明顯，上、下游的孔口周邊均出現(xiàn)小范圍的拉應(yīng)力，但對壩體影響不大。壩體上、下游壩面主應(yīng)力分布均近于對稱。壩體下游面各種正應(yīng)力基本均為壓應(yīng)力，最大主壓應(yīng)力出現(xiàn)在下游壩面419 m高程壩趾處，為7.36 MPa;壩體上游面部分為較小壓應(yīng)力，拉應(yīng)力主要分布在上游壩面的拱端與基巖接觸區(qū)，最大主拉應(yīng)力為1.17 MPa，在463.875 m高程上游面左拱端處，上、下游應(yīng)力第一主應(yīng)力分布圖如圖4，圖5所示。

圖4 上游壩面第一主應(yīng)力等值線圖

圖5 下游壩面第一主應(yīng)力等值線圖

5 結(jié)語

文章采用ANSYS的三維有限元建模，針對高拱壩壩頂開設(shè)大表孔泄水對壩體應(yīng)力變形的影響問題進行研究，分析了拱壩在荷載作用下各個部位的應(yīng)力分布狀況，主要分析結(jié)論可歸納如下:

1）壩體整體的變形趨勢是向下游方向，最大位移為1 cm，發(fā)生在溢流堰頂區(qū)域;

2）拱壩表孔處局部應(yīng)力集中現(xiàn)象十分明顯，上、下游的孔口周邊均出現(xiàn)小范圍的拉應(yīng)力，但對壩體影響不大。

［1］潘家錚，何璟.中國大壩50年［M］.北京：中國水利水電出版社，2000.

［2］朱伯芳.拱壩設(shè)計與研究［M］.北京：中國水利水電出版社，2002.

［3］鄒爽，李小可.大表孔拱壩應(yīng)力分析［J］.貴州水力發(fā)電，2007，21（3）：44-46.

［4］劉先珊，周創(chuàng)兵，王軍.復(fù)雜條件下高拱壩應(yīng)力及壩肩穩(wěn)定分析［J］.巖土力學(xué)，2008，29（1）：225-229.

［5］張朝暉.ANSYS 11.0有限元分析理論與工程應(yīng)用［M］.北京：電子工業(yè)出版社，2008.

Stress analysis on an overflow hyperbolic arch dam based on ANSYS

HUANG De-cai1DING Ze-lin2

（1.Water Conservancy Bureau ofWangmo County，Southwest Qian of Guizhou Province，Wangmo 552300，China; 2.North China University ofWater Resources and Electric Power，Zhengzhou 450011，China）

This papermade stress analysis on a concrete hyperbolic arch dam based on three-dimensional finite element calculation software ANSYS，through the analysis gained，the largestdisplacementofarch dam occurred in the surface hole position of dam，but the influence ofopening surface holes in dam head to dam stresswere local，mainly concentrated in the outlets circum，had almostno effect to the dam.And summarized the laws，provided important scientific basis to engineering design，construction and reinforcement treatment.

hydraulic structure，stress analysis，finite element，arch dam

TV541

A

10.13719/j.cnki.cn14-1279/tu.2012.19.063

1009-6825（2012）19-0249-02

2012-05-03

黃德才（1971-），男，工程師