毛璐明，田威，楊玉群

（1.陜西省電力公司電力經(jīng)濟技術(shù)研究院，陜西西安 710075；2.長安大學建筑工程學院，西安 710061）

深厚覆蓋層上土石圍堰的應力變形研究

毛璐明1，田威2，楊玉群1

（1.陜西省電力公司電力經(jīng)濟技術(shù)研究院，陜西西安 710075；2.長安大學建筑工程學院，西安 710061）

以修建于深厚覆蓋層上的土石圍堰為研究對象，利用有限元軟件ANSYS分步模擬堰體填筑和大壩基坑開挖過程。重點研究圍堰竣工及基坑開挖的過程中堰體與混凝土防滲墻的應力變形特征。計算結(jié)果表明：整個施工過程中堰體與防滲墻沿三個坐標軸方向的位移分布規(guī)律良好，數(shù)值均在變形合理范圍內(nèi)，且各應力特征值均未超過設計標準，該圍堰與防滲墻從施工期到運行期都是安全穩(wěn)定的。

土石圍堰；混凝土防滲墻；深厚覆蓋層；有限元；應力變形

圍堰是臨時擋水建筑物，修建的目的是為了使水利樞紐建筑物能在干地上進行施工和檢修。它的受力及變形狀態(tài)對整個工程的施工安全起決定作用[1-2]。

20世紀90年代，我國興建的長江三峽深水高土石圍堰被國務院三峽建設委員會初步設計審查核心專家組確定為唯一的必須進行單項技術(shù)設計的臨時工程[3]。該圍堰從研究、設計到建設運行，歷經(jīng)了18年[4-5]，其中試驗研究與設計分析階段就長達11年。可見圍堰在大型水利工程中的重要性，它的設計與施工正越來越受到人們的關(guān)注。

我國西部地區(qū)水能資源豐富，隨著“西電東送”等重大項目的落實實施，越來越多的大型、特大型水電工程坐落于此，于是大量的土石圍堰工程也破土動工。但該區(qū)域普遍存在壓縮性較高的深厚覆蓋層，大量的工程經(jīng)驗證明[6-9]，修建在深厚覆蓋層上的土石圍堰堰體及其防滲墻的應力變形特性是影響整個圍堰工程安全的關(guān)鍵，因此在施工的全過程中它們的穩(wěn)定性都需要特別關(guān)注。

本文以西部某二等大（2）型水電站的上游圍堰為例，通過模擬堰體分層填筑的過程，研究堰體和混凝土防滲墻的應力變形特性，為圍堰能否安全運行作出評價。

1 工程概況

某水電站位于西部一江的中游段，裝機容量510 MW。大壩為混凝土重力壩，正常蓄水位3 310.0 m，

壩頂高程3 314.0 m，大壩基礎(chǔ)高程3 198.0 m，最大壩高116.0 m，壩頂長377.60 m，采用壩后式地面廠房。本工程為二等大（2）型，主要建筑物為2級，相應導流建筑物為4級。壩址區(qū)為深切的高山峽谷，兩岸基巖出露，巖性以花崗巖為主。壩址河床覆蓋層最大厚度約45.0 m，枯期河床水面寬約120.0 m，水深約4.0 m。

該工程二期大壩基坑上游圍堰采用土工膜斜墻土石圍堰，堰頂高程3 280.0 m，最大堰高39.0 m，堰基采用封閉式砼防滲墻，防滲墻最大深度52.0 m。下游采用心墻土石圍堰，堰頂高程3 260.0 m，最大堰高16.0 m，堰基采用封閉式混凝土防滲墻，防滲墻最大深度48 m。其典型剖面如圖1所示。

圖1 圍堰典型剖面Fig.1 Typical sections of the cofferdam

2 堰體三維有限元模型與計算參數(shù)

2.1 基本假定

1）覆蓋層土體、堰體堆石料為連續(xù)的彈塑性體，符合Mohr-Coulomb條件[10-11]，采用八節(jié)點四面體實體單元。

2）混凝土防滲墻及基巖視為線彈性體，均選用八節(jié)點四面體單元。

3）封閉式混凝土防滲墻與覆蓋層土體，以及堰體堆石料界面采用面-面的摩擦接觸單元。

4）上游的水壓力以梯度荷載的形式施加于堰坡及覆蓋層上。

2.2 有限元模型尺寸及邊界條件

有限元模型計算區(qū)域：上游邊界為上游圍堰防滲帷幕向上游延伸200 m，左岸邊界為左岸堰肩以左100.0 m，右岸邊界為右岸堰肩以右100.0 m，下游邊界自上游圍堰坡腳向下游延165 m，達到壩軸線處。垂直向自弱風化層下限向下延50 m。坐標采用笛卡兒直角坐標系，坐標原點位于圍堰軸線與堰體典型剖面交點的0標高處。以橫河向為x軸，指向左岸為正向；以順河向為z軸，指向下游為正向；以垂直向為y軸，垂直向上為正。計算模型的約束條件為：覆蓋層表面為自由邊界，不受任何約束；計算模型底部邊界為3個方向的全約束；其余邊界面均施加法向約束。模型如圖2所示。其中單元數(shù)218 717個，結(jié)點數(shù)39 609個。

圖2 有限元網(wǎng)格示意圖Fig.2 FEM mesh of the cofferdam

2.3 計算參數(shù)的選取

上游圍堰三維應力應變計算區(qū)域內(nèi)共涉及到13種材料，它們分別為：截流戧堤、堰體堆石料、噴混凝土、堰坡干砌石、帷幕灌漿、導流明渠混凝土、混凝土防滲墻、覆蓋層2（alQ4）、覆蓋層1（fglQ3）、基巖Ⅳ、基巖Ⅲ1、基巖Ⅲ2和基巖Ⅱ。具體參數(shù)見表1。

表1 材料參數(shù)Tab.1 Material parameters

防滲墻與覆蓋層、基巖，以及堰體堆石料界面的接觸面單元的具體計算參數(shù)如表2所示。

3 填筑與計算過程

利用ANSYS軟件中“殺死”和“激活”單元的功

能，對大壩基坑開挖的彈塑性應力變形進行研究，數(shù)值模擬的具體過程如下：

表2 計算采用的接觸面模型參數(shù)Tab.2 Interface model parameters in the calculation

第一步：施工截流戧堤和防滲墻施工平臺；

第二步：施工混凝土防滲墻；

第三步：施工高程3 260 m以下上游圍堰，堰前及導流明渠中按水位高程3 248.00 m施加水壓力；

第四步：施工高程3 264 m以下上游圍堰；

第五步：施工高程3 268 m以下上游圍堰；

第六步：施工高程3 272 m以下上游圍堰；

第七步：施工高程3 276 m以下上游圍堰；

第八步：施工高程3 280 m以下上游圍堰、土工膜及堰坡干砌石（圍堰竣工）；

第九步：堰后基坑開挖至高程3 234 m，堰前水位上升至高程3 277.90 m；

第十步：堰后基坑開挖至高程3 219 m；

第十一步：堰后基坑開挖至高程3 210 m；

第十二步：堰后基坑開挖至高程3 198 m。

4 計算結(jié)果分析

計算結(jié)果按照圍堰竣工與大壩基坑開挖完成兩個階段進行整理。ANSYS軟件中規(guī)定拉應力為正，壓應力為負。主要分析堰體及防滲墻的變形規(guī)律，以及安全穩(wěn)定性。

4.1 堰體變形及應力

從圖3中可以看出，圍堰碾壓到設計高程3 280 m時，堰體的最大沉降量約為34.46 cm。最大沉降量發(fā)生在約1/3堰高處，堰體的沉降曲線呈以堰體中心線為中心的同心圓圈分布。垂向變形規(guī)律與一次加載，壩頂處沉降量最大明顯不同。這是由于圍堰逐級填筑，假定變形在施工中瞬時完成，則下部結(jié)構(gòu)的自重不影響上部結(jié)構(gòu)的變形。當堰體施工到某一高程，這個高程以下土重引起的沉降已經(jīng)發(fā)生，此高程上各點如果發(fā)生位移，僅僅是其上土的自重所引起的。壩頂以上不再有荷載，故豎向位移為零。

圖3 圍堰竣工時沉降等值線圖Fig.3 The subsidence isogram upon the cofferdam completion

擋水后開始進行大壩基坑開挖。如圖4所示，開挖完成時圍堰的最終沉降量為38.83 cm，較圍堰竣工時有所增加，不過仍小于堰高的1%，在土石圍堰變形的合理范圍內(nèi)。

圖4 大壩基坑開挖完成時沉降等值線圖Fig.4 The subsidence isogram upon the excavation completion

堰體上游的側(cè)向水平位移指向上游，下游水平位移指向下游，基本呈對稱分布，見圖5（a）、（b）。

圖5 圍堰順河向位移等值線圖Fig.5 The displacement contour of the cofferdam

堰體竣工時，上游水平位移最大值10.46 cm，下游水平位移最大值14.76 cm。擋水后由于水壓力的作用，整個堰體均發(fā)生向下游的水平變形，指向上游的位移明顯減小，下游水平位移有所增大。大壩基坑開挖完成時，上游水平位移最大值減小到5.35 cm，下游水平位移最大值增加了3.34 cm。

如圖6（a）、（b）所示，堰體從竣工到大壩開挖完成，最大應力水平從0.49逐步增大到0.76。根據(jù)摩爾庫倫準則可以判定從填筑碾壓到大壩基坑開挖的整個過程中，堰體邊坡都是安全的。

4.2 混凝土防滲墻的應力變形

圖7反映出從圍堰填筑至高程3 272 m到大壩基坑開挖完成，混凝土防滲墻沿x軸方向的變形基本規(guī)律。施工的整個過程中，防滲墻x、y和z 3個方向的位

移分布規(guī)律良好。橫河x軸方向的位移左右岸對稱分布，最大值出現(xiàn)在防滲墻頂部，最大位移量約為-0.52 cm/0.50 cm。

圖6 圍堰應力水平等值線圖Fig.6 The stress level contour of the cofferdam

圖7 防滲墻x方向的最大水平位移Fig.7 The max displacement in x direction of the cutoff wall

垂直向最大沉降量出現(xiàn)在防滲墻頂部，呈現(xiàn)以防滲墻頂部中心為圓心的半圓形同心圓分布如圖8所示?？⒐r垂直向最大沉降量為2.85 cm。在模擬大壩基坑開挖的過程中覆蓋層卸荷，使得防滲墻豎向變形產(chǎn)生少量的回彈，最大沉降量降至2.02 cm。

圖8 防滲墻垂向位移等值線圖Fig.8 The vertical displacement contour of the cutoff wall

觀察圖9可以發(fā)現(xiàn)，防滲墻順河水平方向受力相當于一塊文克萊彈性地基板，因此順河水平方向位移呈現(xiàn)以防滲墻中心為圓心的同心圓，最大位移位置在防滲墻2/3高度處。圍堰竣工時，由于堰體自重壓縮覆蓋層產(chǎn)生側(cè)向位移，導致順河水平方向位移最大值為11.03 cm，指向上游。圖10反映了計算過程六至十二步，防滲墻沿z軸方向的變形規(guī)律。從圖中可以看出第九步圍堰擋水后，水壓力有把防滲墻推向下游的作用，因此防滲墻指向上游的位移顯著減小，最終為5.1 cm。之后基本保持不變，說明基坑開挖對防滲墻z軸方向的變形影響很小。

圖9 防滲墻z向位移等值線圖Fig.9 The displacement contour in z direction of the cutoff wall

圖10 防滲墻z方向的最大水平位移Fig.10 The max displacement in z direction of the cutoff wall

由于防滲墻兩側(cè)地層不均勻，因此，防滲墻應力分布比較復雜，規(guī)律性不明顯。各應力特征值見表3。

表3 防滲墻應力極值Tab.3 Maximum values of stress

根據(jù)表3可以看出：從圍堰碾壓至設計高程到大壩基坑開挖完成的整個過程中，防滲墻的壓應力均未超過C15混凝土的軸心抗壓強度設計值7.2 MPa，且未出現(xiàn)拉應力，產(chǎn)生拉裂縫的幾率很小。防滲墻的最大應力水平也遠遠小于1，根據(jù)摩爾庫倫準則可以判定在整個施工過程中防滲墻的受力是安全的[6]。

5 結(jié)論

本文通過對某深厚覆蓋層上土石圍堰及防滲

墻的受力變形分析，可以得出以下結(jié)論：

1）在數(shù)值模擬整個施工工程中，圍堰的最大沉降量小于堰高的1%，在堰體豎向變形的合理范圍內(nèi)。由于深厚覆蓋層具有一定壓縮性，因此堰體最大沉降出現(xiàn)的位置下移到約1/3堰高處。

2）從堰體填筑碾壓到大壩基坑開挖的整個過程中，堰體的應力水平逐步增大，但始終沒有超過0.8，根據(jù)摩爾庫倫準則可以判定，堰體邊坡都是安全的。

3）施工的整個過程中，防滲墻x、y、z 3個方向的位移分布規(guī)律良好，墻體所受應力均在合理范圍內(nèi)，最大應力水平僅為0.15，可以判定防滲墻的施工在理論上是安全的。

[1]鄭守仁，王世華，夏仲平，等.導流截流及圍堰工程[M].北京：中國水利電力出版社，2004.

[2]石蹈波.大型水電站高土石圍堰安全研究[D].南京：河海大學，2007.

[3]沈珠江，劉松濤.三峽二期高土石圍堰應力應變分析研究[J].人民長江，1996，27（10）：1-4.SHEN Zhujiang，LIU Songtao.Study on stress and strain of TGP StageⅡcofferdam[J].Yangtze River，1996，27（10）：1-4（in Chinese）.

[4]包承綱.三峽工程二期圍堰的理論與實踐[R].武漢:長江科學院，2000.

[5]劉松濤，曾玲，龔壁衛(wèi).三峽工程二期圍堰防滲墻結(jié)構(gòu)非線性應力應變研究[J].長江科學院院報，1992，9（2）: 1-10.LIU Songtao，ZENG Ling，GONG Biwei.A research on non-linear stress-strain for the 2nd-stage cofferdam with concrete diaphragm wall of the three gorges project[J].Journal of Yangtze River Scientific Research Institute，1992，9（2）:1-10（in Chinese）.

[6]徐晗，汪明元.深厚覆蓋層300 m級超高土質(zhì)心墻壩應力變形特征[J].巖土力學，2008，29（增）:64-68.XU Han，WANG Mingyuan.3D stress and deformation characteristics of 300 m level high earth core dams with thick overburden layer[J].Rock and Soil Mechanics，2008，29（S）:64-68（in Chinese）.

[7]盧小春.深厚覆蓋層上深水高土石圍堰結(jié)構(gòu)性態(tài)研究[D].武漢：武漢大學，2011.

[8]熊歡，王清友，高希章，等.沙灣水電站一期圍堰塑性混凝土防滲墻應力變形分析[J].水力發(fā)電學報，2010，29（2）：197-203.XIONG Huan，WANG Qingyou，GAO Xizhang，et al.First stage cofferdam of Shawan hydropower station of plastic concrete anti-seepage wall stress deformation analysis[J].Journal of Hydroelectric Power Generation，2010，29（2）：197-203（in Chinese）.

[9]張祥，車明杰，樊艷欣，等.深厚覆蓋層地基高土石圍堰應力變形分析[J].電網(wǎng)與清潔能源，2015，31（12）: 121-125.ZHANG Xiang，CHE Mingjie，F(xiàn)AN Yanxin，et al.Stress and deformation analysis of high earth-rock cofferdam on deep overburden foundation[J].Power System and Clean Energy，2015，31（12）:121-125（in Chinese）.

[10]陶紀南.莫爾-庫侖準則的應用范圍及臨界應力狀態(tài)探討[J].巖土工程學報，1990，12（5）:76-83.TAO Jinan.The scope of application of Mohr-Coulomb criteria and research of critical stress state[J].Chinese Journal of Geotechnical Engineering，1990，12（5）:76-83（in Chinese）.

[11]劉祖典，黨發(fā)寧.土的彈塑性理論基礎(chǔ)[M].西安:世界圖書出版公司，2002.

Study on Stress and Deformation of Earth-Rock Cofferdam on Deep Overburden

MAO Luming1，TIAN Wei2，YANG Yuqun1
（1.Shaanxi Electric Power Corporation Economic Research Institute，Xi’an 710075，Shannxi，China；2.School of Civil Engineering，Chang’an University，Xi’an 710061，Shaanxi，China）

This paper takes the earth-rock cofferdam on the deep overburden as its main study object and numerically simulates the cofferdam filling and dam excavation process using FEM software ANSYS.The research is mainly focused on the stress and deformation characteristic of the cofferdam and concrete cutoff wall during completion and excavation processes.The results show that displacements of the cofferdam and cutoff wall are well distributed in three coordinate axes during the entire construction with all the values within the acceptable dis-formation limit and stress characteristic values not exceeding the design standard，and the cofferdam and cutoff wall are safe and stable both in construction period and operation period.

earth-rock cofferdam；concrete cutoff wall；deep overburden；finite element method；stress and deformation

2016-03-15。

毛璐明（1985—），女，工程師，主要從事電網(wǎng)工程建設等相關(guān)工作。

（編輯 黃晶）

陜西省自然科學基金資助項目（2015JM5160）。

Project Supported by the Natural Science Foundation of Shaanxi Province（2015JM5160）.

1674-3814（2016）09-0141-05

O319.56

A