馬富強，張笮娜，李海波

(貴州省水利水電勘測設(shè)計研究院有限公司，貴陽 550002)

1 概 述

某水庫是貴州省“三位一體”綜合規(guī)劃的重點水源工程，工程主要功能是工業(yè)和城市供水，兼顧灌溉，年可供水量4 275×104m3/a。擋水建筑物為C25碾壓混凝土拋物線型雙曲拱壩，布置在“一線天”的非典型峽谷型河谷地帶，河谷異常狹窄，底寬6m，正常水位時的河谷寬85m，寬高比0.61。壩軸線長137.9m，壩頂寬6m，底寬24m，壩頂高程1 290.50m，河床建基面開挖高程1 137.5m，最大壩高153m，拱冠厚高比0.157[1]，是貴州省內(nèi)水利行業(yè)第一高薄拱壩，也是國內(nèi)最窄深的高薄雙曲拱壩。該水庫拱壩應(yīng)力采用多拱梁法進行分析計算，進而初步確定壩體體形，設(shè)計采用《拱壩體形優(yōu)化設(shè)計——模型、方法與程序》[2]ADASO電算程序計算。將拱冠梁分為9拱19梁，作用于壩體的主要荷載有水荷載、泥沙壓力、自重和溫度荷載等。

多拱梁法通過拱、梁劃分及多層分載求解，可以初步擬定壩體體型及核算壩肩穩(wěn)定及拱座應(yīng)力，其計算成果作為衡量拱壩強度安全的主要標(biāo)準(zhǔn)，為壩體結(jié)構(gòu)初步設(shè)計打下基礎(chǔ)。但多拱梁法由于采用結(jié)構(gòu)力學(xué)的基本假定，計算成果不免帶有一定的近似性，且由于受計算假定的限制，多拱梁分載法難以取得滿意的成果，而有限元法可以反映各種復(fù)雜的因素，取得比較接近于實際情況的成果。隨著計算機技術(shù)的發(fā)展和有限元計算軟件的推廣和普及，工程設(shè)計人員利用有限元進行拱壩應(yīng)力分析已不存在制約性因素，采用有限元進行拱壩應(yīng)力計算是拱壩設(shè)計的一個大趨勢，從而為拱壩結(jié)構(gòu)體型設(shè)計、細(xì)部結(jié)構(gòu)設(shè)計、溫控設(shè)計，以及壩體分區(qū)、基礎(chǔ)處理、溫控措施的選用等提供技術(shù)支持，對類似建設(shè)條件的同類型擋水建筑物設(shè)計及應(yīng)力復(fù)核分析，具有一定的參考價值。

2 應(yīng)力復(fù)核工況選取

按照文獻【3】第7.2.1條規(guī)定，對于該項目150m級超高拱壩應(yīng)力除采用多拱梁法進行分析計算外，同時還應(yīng)采用三維有限元法對拱壩進行應(yīng)力復(fù)核計算。按照文獻【4】第1.0.5條規(guī)定，該水庫大壩抗震設(shè)防類別為乙類；根據(jù)文獻【4】第7.1.2條的規(guī)定，該項目大壩強度抗震分析方法以靜動力拱梁分載法為基本方法，無需進行有限元動力分析。因此，有限元應(yīng)力復(fù)核針對大壩靜力工況開展。

3 有限元建模及荷載組合

3.1 有限元模型

3.1.1 坐標(biāo)系的選取

整個有限元模型建立在笛卡爾坐標(biāo)系下，X軸為垂直于河流方向，Y軸為順河流方向，Z軸為豎直方向。

3.1.2 計算域的選取

壩前、壩后、左右岸拱端及底部地基取1倍壩高左右。

3.1.3 邊界條件

在壩體和地基構(gòu)成的整個計算域內(nèi)，對基巖的上下游面施加Y方向的約束，左右岸邊界均施加X方向的約束，基巖底部施加X、Y及Z方向的約束。

3.1.4 有限元網(wǎng)格模型

按壩體的實際體型進行建模，由于溢洪道開孔較小，忽略其影響。在有限元模型建立時，充分考慮壩體的形狀、基巖材料的分區(qū)以及多種荷載的施加區(qū)域等。同時，壩基及壩體近似假定為各同向同性、均質(zhì)、連續(xù)的線彈性體。壩體體形共分為76 964個單元、80 908個節(jié)點，整體三維有限元網(wǎng)格圖見圖1。

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

3.2 荷載組合及計算參數(shù)

3.2.1 荷載組合

根據(jù)規(guī)范，混凝土拱壩計算荷載分為基本運用和特殊運用兩類組合，根據(jù)本項目實際情況，擬定4種工況進行計算。分別如下：

1)基本組合。

工況一(基本組合①)：正常水位1 287m+壩體自重+設(shè)計正常溫降+揚壓力+泥沙壓力

工況二(基本組合②)：死水位1 213m+壩體自重+設(shè)計正常溫升+揚壓力+泥沙壓力

工況三(基本組合③)：設(shè)計洪水位1 287m+壩體自重+設(shè)計正常溫升+揚壓力+泥沙壓力

2)特殊組合。

工況四(特殊組合)：校核洪水位1 288.89m+壩體自重+設(shè)計正常溫升+揚壓力+泥沙壓力

3.2.2 應(yīng)力計算參數(shù)

應(yīng)力計算參數(shù)見表1。

表1 應(yīng)力計算參數(shù)表

3.3 應(yīng)力控制標(biāo)準(zhǔn)

按照文獻【3】第7.3.2條規(guī)定，采用有限元法計算時，應(yīng)計算有限元等效應(yīng)力[5]，且求得的壩最大壓應(yīng)力和最大拉應(yīng)力要滿足規(guī)范要求的應(yīng)力控制指標(biāo)。

壩體允許最大壓應(yīng)力：基本組合為6.25MPa、非抗震特殊組合為7.14MPa。

壩體允許最大拉應(yīng)力：基本組合為1.5MPa、非抗震特殊組合為2.0MPa。

4 有限元計算結(jié)果及分析

根據(jù)計算結(jié)果整理如下：應(yīng)力圖的單位為MPa，壓應(yīng)力為負(fù)，拉應(yīng)力為正，均為下游立視圖。

4.1 工況一(基本組合①)

正常水位1 287m+壩體自重+設(shè)計正常溫降+揚壓力+泥沙壓力。見圖2、圖3。

圖2 組合①上游面第一主應(yīng)力圖

圖3 組合①下游面第三主應(yīng)力圖

4.2 工況二(基本組合②)

死水位1 213m+壩體自重+設(shè)計正常溫升+揚壓力+泥沙壓力。見圖4、圖5。

圖4 組合②上游面第一主應(yīng)力圖

圖5 組合②下游面第三主應(yīng)力圖

4.3 工況三(基本組合③)

正常蓄水水位1 287m+壩體自重+設(shè)計正常溫升+揚壓力+泥沙壓力。見圖6、圖7。

圖6 組合③上游面第一主應(yīng)力圖

圖7 組合③下游面第三主應(yīng)力圖

4.4 工況四(特殊組合)

校核洪水位1 288.89m+壩體自重+設(shè)計正常溫升+揚壓力+泥沙壓力。見圖8、圖9。

圖8 特殊組合上游面第一主應(yīng)力圖

圖9 特殊組合下游面第三主應(yīng)力圖

根據(jù)上述有限元計算成果圖，分析如下：

工況1：大壩上游面主拉應(yīng)力表現(xiàn)較為突出，基本上控制在1.5MPa以內(nèi)，但在靠近基礎(chǔ)區(qū)域存在應(yīng)力集中現(xiàn)象，極值為4.51MPa。大壩下游面主壓應(yīng)力表現(xiàn)較為突出，極值為-5.25MPa。

工況2：大壩上游面主拉應(yīng)力表現(xiàn)較為突出，基本上控制在1.5MPa以內(nèi)，但在靠近基礎(chǔ)區(qū)域存在應(yīng)力集中現(xiàn)象，極值為2.64MPa。大壩下游面主壓應(yīng)力表現(xiàn)較為突出，極值為-3.06MPa。

工況3：大壩上游面主拉應(yīng)力表現(xiàn)較為突出，基本上控制在1.5MPa以內(nèi)，但在靠近基礎(chǔ)區(qū)域存在應(yīng)力集中現(xiàn)象，極值為3.56MPa。大壩下游面主壓應(yīng)力表現(xiàn)較為突出，極值為-6.30MPa。

工況4：大壩上游面主拉應(yīng)力表現(xiàn)較為突出，基本上控制在1.5MPa以內(nèi)，但在靠近基礎(chǔ)區(qū)域存在應(yīng)力集中現(xiàn)象，極值為3.59MPa。大壩下游面主壓應(yīng)力表現(xiàn)較為突出，極值為-6.38MPa。

5 有限元等效應(yīng)力結(jié)果及分析

5.1 有限元應(yīng)力極值成果分析

根據(jù)上述應(yīng)力圖可知，通過ANSYS計算出的主拉應(yīng)力與主壓應(yīng)力的極值均相對較高，可以看出較高的拉應(yīng)力及壓應(yīng)力基本上出現(xiàn)在靠近基礎(chǔ)部位。采用三維彈性有限單元法計算拱壩應(yīng)力時，靠近基礎(chǔ)部位存在明顯的應(yīng)力集中問題，本壩上游面有節(jié)點1 415個，有限元主拉應(yīng)力超過控制指標(biāo)約3%，有限元主壓應(yīng)力均未超過控制指標(biāo)；本壩下游面有節(jié)點1 415個，有限元主拉應(yīng)力僅工況2有出現(xiàn)超過控制指標(biāo)的情況(占0.5%)，有限元主壓應(yīng)力超過控制指標(biāo)的最多約2%。在本次有限元應(yīng)力極值成果分析中，同時計算了有限元等效應(yīng)力。

5.2 有限元等效應(yīng)力成果分析

根據(jù)有限元應(yīng)力的分布情況，采用有限元等效應(yīng)力計算分析方式，對計算數(shù)據(jù)進行整理、修正及研究分析后，結(jié)果如下：基本組合下拉應(yīng)力最大值為1.24MPa，壓應(yīng)力最大值為4.65MPa；特殊組合下拉應(yīng)力最大值為1.3MPa，壓應(yīng)力最大值為-4.62MPa，均滿足規(guī)范應(yīng)力標(biāo)準(zhǔn)控制要求。

6 結(jié) 語

對于極窄峽谷超高薄拱壩，通過三維有限元結(jié)構(gòu)應(yīng)力仿真計算，可以有效指導(dǎo)壩體體型設(shè)計[6]、分縫設(shè)計及溫控設(shè)計[7]，以達到優(yōu)化壩體結(jié)構(gòu)、節(jié)約工程量、節(jié)省工程投資的目的。同時，還解決了常規(guī)多拱梁法未能解決的其他諸多問題。另外，有限元應(yīng)力仿真計算，也是對拱壩應(yīng)力分析最為有效的一種復(fù)核驗證方式，對同類型擋水建筑物設(shè)計，具有一定的借鑒價值。