楊 勝, 彭 顯 國
(1.中國電建集團成都勘測設(shè)計研究院有限公司,四川 成都 610072;2.國電大渡河流域水電開發(fā)有限公司,四川 成都 610041)
目前我國西南地區(qū)涌現(xiàn)出一大批已建和待建的高拱壩,泄洪量不斷攀升,壩頂開設(shè)表孔尺寸日益增大。壩頂開設(shè)大孔口,造成壩頂不能形成連續(xù)性拱圈,且壩頂一般較薄,設(shè)置弧形閘門又需要一定的空間,從而要求在表孔的上、下游增設(shè)支承水平梁及閘墩,使得表孔結(jié)構(gòu)變得十分復(fù)雜。此外,孔口會造成局部應(yīng)力集中,可能導(dǎo)致孔口局部混凝土開裂和裂縫的發(fā)展,進而影響壩體的正常運行,對大壩的安全造成影響,是大壩設(shè)計中較為關(guān)注的問題之一。拱壩應(yīng)力狀態(tài)呈典型三維分布特征,孔口簡化為二維計算必然損失精度,無法真實地描述孔口周圍的應(yīng)力分布規(guī)律,分析其影響因素。筆者以某工程的表孔為研究主導(dǎo),建立了不同表孔結(jié)構(gòu)形式的拱壩模型,基于ANSYS 中的子模型技術(shù)進行了三維有限元計算分析,以了解不同孔口尺寸、不同形式的水平梁及閘墩情況下表孔局部應(yīng)力分布情況及其對壩體應(yīng)力變形的影響,從而為高拱壩表孔的布置及優(yōu)化設(shè)計提供了參考依據(jù),對節(jié)省工程投資具有重要意義。
某水電站樞紐位于高山峽谷區(qū),攔河大壩為對數(shù)螺旋線雙曲薄拱壩,最大壩高250 m。壩身泄洪采用表孔、深孔(底孔)聯(lián)合泄流方式。三個泄洪表孔布置于壩頂中部,由進口段、平直段、WES曲線段和反弧段(左、右表孔包括鼻坎段)或斜直段(中表孔包括跌坎)組成。三個進口共用一扇平板檢修門,每孔布置一扇弧形工作門,單孔凈寬13 m;消能工為橢圓及蝙蝠型舌型挑(跌)鼻坎。
拱壩整體模型包括整個壩體、地基及表孔閘墩結(jié)構(gòu)等,沒有考慮壩體廊道及深孔(底孔)的影響。子模型選取表孔所在位置的四個壩段,并通過比較切割邊界上的應(yīng)力結(jié)果與整體模型相應(yīng)位置的結(jié)果是否一致驗證了子模型邊界的合理性。在保證計算精度的前提下,對表孔結(jié)構(gòu)的主要研究部位盡量精確模擬,對次要部位合理簡化。整體模型、子模型網(wǎng)格剖分基本采用八節(jié)點六面體單元,部分通過四面體實體單元進行過渡。拱壩整體、子模型有限元模型見圖1、2。
圖1 拱壩整體有限元模型示意圖
圖2 子模型有限元模型示意圖
計算工況為正常蓄水位、設(shè)計洪水位兩種,主要考慮的荷載有水壓力、揚壓力、泥沙壓力、壩體自重、弧門推力以及溫度荷載等,荷載組合情況見表1。其中水荷載簡化計算,即在上、下游壩面施加靜水壓力,壩底及壩肩施加揚壓力(壩底揚壓力的計算考慮了防滲帷幕及排水管的作用);表孔泄流時根據(jù)水工模型試驗所得水面線簡化施加流道水壓力;弧門擋水時,弧門上游均為相應(yīng)水頭的靜水壓力;檢修情況考慮進口事故檢修門擋水,檢修門上游為相應(yīng)水頭的靜水壓力,檢修門后無水壓力。根據(jù)設(shè)計施工資料計算得到壩體溫度場(泄水時孔內(nèi)為孔口庫水溫邊界,檢修時孔內(nèi)是氣溫邊界),進而計算出溫度引起的初應(yīng)變,而后求得相應(yīng)的初應(yīng)變引起的等效結(jié)點溫度荷載,最后按通常的求解應(yīng)力方法求得溫度應(yīng)力。
表1 荷載組合表
通過對壩體開設(shè)表孔及不開孔模型有限元計算成果進行對比分析可知:壩頂開設(shè)表孔后,由于拱向作用被削弱,孔口周邊的X、Y、Z向應(yīng)力受到很大的影響。上、下游的孔口周邊均出現(xiàn)小范圍的X、Z向拉應(yīng)力;孔口附近的Y向應(yīng)力總的趨勢是壓應(yīng)力較無孔時有所減小,且越接近孔口降低越多,局部出現(xiàn)拉應(yīng)力??卓趥?cè)壁三向均出現(xiàn)大范圍的拉應(yīng)力,但其值不大;孔口角沿應(yīng)力集中現(xiàn)象較明顯,Y、Z向壓應(yīng)力均有較大幅度的增加;孔口各向位移均略有增大?,F(xiàn)從開孔尺寸(不同寬度、高度)、孔口形式(水平梁、不連續(xù)牛腿)等方面分別進行計算分析,進一步研究了表孔應(yīng)力變形分布及其對壩體應(yīng)力變形的影響。由于兩種工況下壩體應(yīng)力變形分布規(guī)律基本一致,故筆者僅對工況1下的計算結(jié)果進行分析。
為了探討開孔寬度(孔口總寬度)對壩體應(yīng)力變形的敏感性分析,分別對壩頂開設(shè)3孔、5孔、7孔、9孔幾種情況進行計算分析。每孔尺寸為13 m×17 m(寬×高),孔之間設(shè)有閘墩,孔口總的寬度占頂拱弧長比例分別為8.36%,13.93%,19.5%,25.1%。各方案孔口關(guān)鍵部位的應(yīng)力值見表2。從計算成果分析可知,壩頂設(shè)置不同數(shù)目的孔口,對壩體下部的影響很小,隨著孔口寬度(總的寬度)的增大,壩體上部X、Z向位移逐步增大,而Y向位移基本不變。
開設(shè)不同數(shù)目孔口時,孔口周邊及附近應(yīng)力均有變化。5孔方案壩體上、下游面孔口周邊及孔口側(cè)壁的拉應(yīng)力值最大;表孔數(shù)目由3孔增加到5孔時,孔口周邊及附近的拉應(yīng)力略有增大;由5孔增到7、9孔時,拉應(yīng)力值逐步降低,但Y向應(yīng)力降幅很小??梢?,總孔口寬度占頂拱弧長比例較小時,加大孔口寬度使得孔口局部的拉應(yīng)力值增大;但是,當(dāng)比例增到20%左右時,再增加總寬度,使得頂部拱圈進一步被削弱,傳遞的力減小,孔口局部應(yīng)力反而開始降低。
為了進一步研究壩體表孔尺寸對壩體應(yīng)力變形分布規(guī)律的影響,以壩頂開設(shè)3孔的情況為例,分別計算孔口高度為7 m、12 m、17 m、22 m四種方案的壩體及表孔局部應(yīng)力變形分布,各方案孔口關(guān)鍵部位的應(yīng)力值見表3。從計算成果可知,改變表孔的孔口高度對壩體下部的影響很小。隨著孔口高度的增大,壩體上部X、Y、Z三向位移均逐步增大,但最大值出現(xiàn)的部位一致??卓诟叨葹?2 m時壩體位移最大,X、Y向拉應(yīng)力最大,Z向壓應(yīng)力最大,而孔口高度為7 m時,Z向拉應(yīng)力最大(1.5 MPa)。
孔口高度由7 m漸變?yōu)?2 m時,孔口局部應(yīng)力變化趨勢為X、Y向應(yīng)力略有增大,Z向應(yīng)力減小,Y向拉應(yīng)力區(qū)域分布位置變化很大。當(dāng)孔口高度與壩高之比達到一定程度(22/250=8.8%)時,進一步增加孔口高度,Y向拉應(yīng)力仍有增加,但X、Z向拉應(yīng)力基本不變。
表2 不同孔口數(shù)目時孔口關(guān)鍵部位應(yīng)力值表 /MPa
注:壓應(yīng)力為負,拉應(yīng)力為正。
表3 不同孔口高度時孔口關(guān)鍵部位應(yīng)力值表 /MPa
注:壓應(yīng)力為負,拉應(yīng)力為正。
表孔流道、上下游懸臂、閘墩、水平大梁、不連續(xù)牛腿等區(qū)域的應(yīng)力變形是設(shè)計中極為關(guān)注的問題,關(guān)系到結(jié)構(gòu)設(shè)計、孔口配筋、體型優(yōu)化等,不同孔口形式勢必影響到孔口及周邊應(yīng)力變形的分布。筆者分別以水平梁式孔口方案和不連續(xù)牛腿式方案為例對此進行探討。
設(shè)置水平梁及閘墩后,表孔邊墩側(cè)壁的X、Y向拉應(yīng)力大幅度增加,主要分布于邊墩與水平梁交接部位的上游及下部,X向最大拉應(yīng)力為3.5 MPa,Y向最大拉應(yīng)力為2.5 MPa;下游懸臂與邊墩交接處外側(cè)Y向拉應(yīng)力增大,最大拉 應(yīng) 力 為
1.9 MPa;Z向應(yīng)力變化較小。表孔其余部位應(yīng)力變化不大。與水平梁方案相比,不連續(xù)牛腿方案的邊墩過流側(cè)壁X向應(yīng)力降低,由3.5 MPa減小到2.7 MPa,且拉應(yīng)力區(qū)下移;Y向拉應(yīng)力明顯降低,最大拉應(yīng)力由2.5 MPa減小到1.7 MPa,且拉應(yīng)力區(qū)域明顯減小,主要分布于牛腿與邊墩交接處及牛腿下部;邊墩進、出口外側(cè)壁的Y向拉應(yīng)力明顯增大,最大拉應(yīng)力為2.9 MPa,比水平梁方案增大1 MPa??梢?,兩種方案各有優(yōu)越性及不足,應(yīng)根據(jù)泄水、交通、金結(jié)安裝等要求綜合考慮采用。
研究結(jié)果表明:拱壩壩頂開設(shè)表孔及增設(shè)水平梁和閘墩對壩體應(yīng)力變形影響是局部的,主要集中在孔口周邊及附近,而對壩體的下部幾乎沒有影響。這是由于孔高與壩高之比(17 m/250 m)較小之故,僅為0.07。增設(shè)弧門支撐結(jié)構(gòu)及閘墩增加了孔口附近壩體的厚度,能緩解孔口周邊的應(yīng)力狀態(tài),明顯改善了孔口對壩體的削弱作用,但表孔邊墩側(cè)壁的X、Y向拉應(yīng)力大幅度增加,水平梁上產(chǎn)生了一定的應(yīng)力集中,孔口周邊強度需進一步加強,水平梁及不連續(xù)牛腿方案各有優(yōu)勢,在滿足泄水、交通、金屬結(jié)構(gòu)安裝要求的前提下,有待進一步優(yōu)化表孔結(jié)構(gòu)形式,改善應(yīng)力狀態(tài)。
參考文獻:
[1] 陳 進,黃 薇,丁 茜.大壩矩形孔口應(yīng)力狀態(tài)分析[J].水力發(fā)電,2000.46(12):21-24.
[2] 趙小蓮,張仲卿.增設(shè)孔口牛腿對雙曲拱壩應(yīng)力的影響[J].廣西大學(xué)學(xué)報(自然科學(xué)版),2003.28(4):333-336.
[3] 車 軼,宋玉普.混凝土高拱壩孔口三維非線性分析[J].大連理工大學(xué)學(xué)報,2003.23(2):218-221.
[4] 王進廷,楊 劍,王吉煥,等.二灘拱壩應(yīng)力仿真及參數(shù)敏感性分析[J].水利學(xué)報,2007.38(7):832-837.