龐牧華，李瑞忠

(1.安徽省·水利部淮河水利委員會水利科學研究院，安徽 合肥 230088； 2.山東電力建設(shè)第三工程有限公司，山東 青島 266100)

1 概述

拱壩是一種先進的壩型，在外荷載作用下，以其三面受到巖體約束的高次超靜定殼體結(jié)構(gòu)，借助拱的作用，將大部分的壓力以水平推力方式傳至壩肩兩側(cè)的巖體，少部分荷載靠自重承擔[1]，它能適當?shù)乩貌牧系膹姸葋慝@得自身的穩(wěn)定，是混凝土壩中最為經(jīng)濟和安全的一種壩型，因而得到了廣泛的應(yīng)用。

拱壩的壩體應(yīng)力分析基本方法為拱梁分載法，1級、2級拱壩和高拱壩或情況比較復(fù)雜的拱壩，除應(yīng)采用拱梁分載法計算外，還應(yīng)進行線彈性有限元法分析，必要時可采用非線性有限元法進行分析[2]。近年來隨著仿真技術(shù)的日益發(fā)展，使得計算精度不斷提高，利用ANSYS三維有限元分析軟件可將分析方法精細化、計算機化。軟件的積分功能能夠滿足高精度積分的要求，可以充分地考慮有限元計算帶來的應(yīng)力集中問題，實現(xiàn)拱壩的等效應(yīng)力計算[3]，最終使得應(yīng)力分析結(jié)果在精度、可靠度和運算速度上均取得很大的提升。

文章依據(jù)拱壩有限元等效應(yīng)力計算的原理，基于ANSYS有限元分析軟件，建立三維有限元模型，對某砌石拱壩壩體應(yīng)力進行分析，壩體以及壩基材料均假定為線彈性材料，單元類型采用Solid45單元，可以得到正常蓄水工況、設(shè)計工況和校核工況下壩體主拉應(yīng)力和主壓應(yīng)力及應(yīng)力分布情況。

2 分析方法

使用有限元法計算過程中，得到的拱壩壩體應(yīng)力是整體坐標系中的應(yīng)力，可以在水平拱圈的節(jié)點建立局部坐標系，原點設(shè)在中心線上，局部坐標系與整體坐標系的關(guān)系詳見文獻[4]。在拱的中心線上對梁的水平界面取單位寬度，沿厚度方向?qū)α旱膽?yīng)力和彎矩進行積分，即可得到梁的水平截面上的內(nèi)力。通過材料力學的計算方法，可以得到壩體內(nèi)部拱和梁的等效應(yīng)力及剪應(yīng)力，并計算出壩體的主應(yīng)力，詳見文獻[4]，在ANSYS中計算等效應(yīng)力方法詳見文獻[5]。

3 工程算例

3.1 工程概況

某水庫位于皖南山區(qū)，工程級別為Ⅲ級，大壩為3級建筑物，采用30年一遇設(shè)計洪水標準，200年一遇校核洪水標準。水庫正常蓄水位為553 m，死水位為534 m，設(shè)計洪水位為555.40 m，校核洪水位為556.50 m。攔河壩壩體為漿砌石雙曲拱壩，壩頂高程557.0 m，壩頂弦長65.5 m，厚高比0.19，壩底高程527.0 m，最大壩高30 m，壩頂厚2.5 m，壩底厚6 m。拱壩投產(chǎn)運行多年，運行良好，目前拱壩下游中間處存在滲水現(xiàn)象。

壩址巖層均為堅硬的條痕狀混合巖，巖體內(nèi)沒有風化軟弱帶和斷裂破碎帶，整個壩址巖層不具明顯狀特征，巖石可視為成分和結(jié)構(gòu)均一的塊狀體。在右壩肩上方偏下游的NEE向陡傾節(jié)理密集帶向山體內(nèi)斷續(xù)延伸，此節(jié)理不影響壩肩穩(wěn)定。壩區(qū)未見斷裂構(gòu)造，細小節(jié)理、裂隙走向多垂直河床流向，且傾角較陡，節(jié)理面后期膠結(jié)呈閉合狀。壩基巖層風化較淺，左壩肩及左河床壩基上部巖層以強風化為主，往下為弱風化，強～弱風化帶厚度一般3 m～5 m；右壩肩及右河床壩基上部巖層以弱風化為主，弱風化厚度一般為1 m～2 m。壩區(qū)河谷右岸有少量風化卸荷裂隙分布，卸荷裂隙深度一般小于2 m。

工程區(qū)地震烈度為Ⅵ級，地震動峰值加速度為0.05g。

3.2 計算模型

本次壩體應(yīng)力分析的有限元計算模型包括漿砌石壩體及地基，建模范圍為壩體上下游及兩岸各取150 m，壩基取100 m，壩體線性有限元法單元類型采用Solid45單元進行模擬(見圖1)。計算模型基礎(chǔ)約束采用地基底部三向約束，上下游面及左右側(cè)采用法向約束；x坐標軸為水流方向，定義順水流向為正方向，y坐標軸為壩體的高度方向，以高度增加為正方向，水平面內(nèi)垂直水流方向為z軸，以指向右岸為正，坐標系原點位于壩體上游面右側(cè)位置。

3.3 材料參數(shù)

壩體及壩基材料均假定為線彈性材料，見表1。大壩有限元模型劃分材料見圖2。

表1 壩體材料及壩基巖體計算參數(shù)

3.4 應(yīng)力控制標準

根據(jù)《砌石壩設(shè)計規(guī)范》和《混凝土拱壩設(shè)計規(guī)范》規(guī)定，用有限元法計算時，應(yīng)補充計算。按“有限元等效應(yīng)力”求得的壩體主拉應(yīng)力和主壓應(yīng)力，應(yīng)符合下列應(yīng)力控制指標的規(guī)定[6]：

1)容許壓應(yīng)力：對于基本荷載組合，安全系數(shù)采用3.5；對于非地震情況特殊荷載組合，安全系數(shù)采用3.0。2)容許拉應(yīng)力：對于基本荷載組合，拉應(yīng)力不得大于1.5 MPa；對于非地震情況特殊荷載組合，拉應(yīng)力不得大于2.0 MPa。

3.5 計算工況及荷載組合

砌石拱壩荷載組合及各工況見表2。

表2 計算工況及荷載組合

工況一為正常蓄水位工況，上游水位為553.00 m，工況二為設(shè)計水位工況，上游水位為555.40 m，工況三為校核水位工況，上游水位為556.50 m。三種工況壩體下游均按照無水考慮。

4 計算結(jié)果與分析

本次砌石拱壩模型的有限元計算分析過程中，規(guī)定正值代表拉應(yīng)力，負值代表壓應(yīng)力，應(yīng)力單位為Pa。

4.1 工況一應(yīng)力分布

工況一主要研究砌石拱壩在正常運行期壩體應(yīng)力的分布情況，計算成果提供了砌石拱壩壩體上下游壩面第一、第三主應(yīng)力云圖，見圖3～圖6。

拱壩正常運行期，壩體出現(xiàn)第一主應(yīng)力較大值的范圍為：上游面的左右岸壩肩及壩踵處，下游面的壩體中心處；出現(xiàn)第三主應(yīng)力較大值的范圍為：上游面的壩肩及壩基處，下游面的壩肩、壩體中心及壩趾處。

4.2 工況二應(yīng)力分布

工況二主要研究砌石拱壩在遭遇設(shè)計洪水工況應(yīng)力的分布情況，計算成果提供了砌石拱壩壩體上下游壩面第一、第三主應(yīng)力云圖，見圖7～圖10。

拱壩遭遇設(shè)計洪水位期，壩體出現(xiàn)第一主應(yīng)力較大值的范圍為：上游面的左右岸壩肩及壩踵處，下游面的壩面處；出現(xiàn)第三主應(yīng)力較大值的范圍為：上游面的壩肩及壩基處，下游面的壩肩、壩體中部及壩趾處。

4.3 工況三計算成果

工況三主要研究砌石拱壩在遭遇校核洪水工況應(yīng)力的分布情況，計算成果提供了砌石拱壩上下游壩面第一、第三主應(yīng)力云圖，見圖11～圖14。

拱壩遭遇校核洪水位期，壩體出現(xiàn)第一主應(yīng)力較大值的范圍為：上游面出現(xiàn)在左右岸壩肩及壩踵處，下游面出現(xiàn)在壩面處；出現(xiàn)第三主應(yīng)力較大值的范圍為：上游面壩肩及壩基處，下游面壩肩、壩體中部及壩趾處。

4.4 計算成果分析

由以上計算結(jié)果分析可得，砌石拱壩在各工況下，壩體出現(xiàn)第一主應(yīng)力時，上游面出現(xiàn)在左右岸壩肩及壩踵處，下游面出現(xiàn)在壩面處，其中上游壩踵處出現(xiàn)最大拉應(yīng)力；壩體出現(xiàn)第三主應(yīng)力時，上游面出現(xiàn)在壩肩及壩基處，下游面出現(xiàn)在壩肩、壩體中部及壩趾處，其中下游壩趾處出現(xiàn)最大壓應(yīng)力；隨水位增高最大主拉應(yīng)力和主壓應(yīng)力均呈增大趨勢。經(jīng)等效后的上游壩踵處拉應(yīng)力正常蓄水位工況為1.35 MPa，設(shè)計洪水位工況為1.71 MPa，校核洪水位工況為1.82 MPa；下游壩趾處壓應(yīng)力正常蓄水位工況為2.52 MPa，設(shè)計洪水位工況為2.69 MPa，校核洪水位工況為3.22 MPa。最大拉應(yīng)力、壓應(yīng)力的位置及變化符合拱壩受力規(guī)律，因而該種方法的計算結(jié)果是可靠的。

5 結(jié)語

文章根據(jù)拱壩有限元等效應(yīng)力計算的原理，通過ANSYS有限元分析軟件，對某砌石拱壩正常運行工況、設(shè)計水位工況及校核水位工況下壩體應(yīng)力進行分析，得到結(jié)論如下：三種工況下壩體主拉應(yīng)力和主壓應(yīng)力分布區(qū)域大致相同，且隨擋水水位增加壩體主拉應(yīng)力和主壓應(yīng)力呈增加趨勢，壩體最大拉應(yīng)力出現(xiàn)在上游壩踵處，最大壓應(yīng)力出現(xiàn)在下游壩趾處。根據(jù)計算成果可知該砌石拱壩壩體應(yīng)力滿足應(yīng)力控制標準的規(guī)定，計算結(jié)果反映了砌石拱壩的實際工作狀態(tài)，為今后砌石拱壩的設(shè)計及除險加固措施選擇提供了幫助。