李燕波

(新疆水利水電勘測設(shè)計研究院，烏魯木齊 830000)

綜合性水利樞紐工程是國民經(jīng)濟(jì)和社會發(fā)展的重要基礎(chǔ)設(shè)施，土石壩作為當(dāng)?shù)夭牧蠅?，具有可就地取材、結(jié)構(gòu)簡單、對地基要求低、施工技術(shù)簡單等特點，在我國得到廣泛的推廣和建設(shè)[1]。2019年初，全國水利工作會議上提出“水利工程補(bǔ)短板、水利行業(yè)強(qiáng)監(jiān)管”的基礎(chǔ)論調(diào)，自此水利行業(yè)設(shè)計進(jìn)入了“精細(xì)設(shè)計、理論研究”的新時代[2]。在以往的工程設(shè)計中，關(guān)于土石壩沉降變形，往往是簡單的根據(jù)公式或經(jīng)驗估算壩體沉降量。鑒于此，本文以新疆某水利樞紐工程為例，分別應(yīng)用公式法、二維有限元法和三維有限元法分析土石壩的沉降變形特點，結(jié)合已建工程實例，驗證計算的可靠性，以期為類似工程提供借鑒。

1 工程概況

該工程位于新疆某河山區(qū)中游河段上，是該河上的控制性工程，起著龍頭水庫的重要作用，主要承擔(dān)防洪、灌溉和發(fā)電的工程任務(wù)。水庫總庫容為0.927×108m3，正常蓄水位為1 894.0 m，死水位1 873.0 m，電站總裝機(jī)容量為26 MW，多年平均年發(fā)電量為0.777×108kW·h。工程由碾壓式瀝青混凝土心墻壩、溢洪道、泄洪沖沙洞、發(fā)電引水系統(tǒng)及電站廠房等組成，工程為Ⅲ等中型工程，最大壩高75 m，為2級建筑物，各泄水建筑物及發(fā)電引水建筑物、電站廠房為3級建筑物，次要建筑物為4級建筑物，臨時建筑物為5級建筑物。地震設(shè)防烈度Ⅷ度。壩體典型橫剖面見圖1。

圖1 壩體典型橫剖面

2 公式法計算土石壩沉降變形

根據(jù)《碾壓式土石壩設(shè)計規(guī)范》(SL274-2001)中的規(guī)定，壩體和壩基沉降計算應(yīng)采用分層總和法。本工程壩基可壓縮層厚度約為15 m，最大斷面底寬約260 m，壩基可壓縮層厚度與壩底寬度的比值為0.06，小于0.1，故可不考慮壩體荷載引起的附加應(yīng)力在壩基的應(yīng)力擴(kuò)散，取壩頂以下的最大壩體自重應(yīng)力作為壩基的附加應(yīng)力[3]。采用非黏性土的公式計算沉降量：

式中：S∞為壩體或壩基的最終沉降量；pi為第i計算土層由壩體荷載產(chǎn)生的豎向應(yīng)力(本工程每層土體計算厚度取5.0 m)；Ei為第i計算土層的變形模量。

變形模量E只能通過原位試驗獲得，但部分工程由于工期和經(jīng)費(fèi)的限制，只能獲得土體的壓縮模量Es，可根據(jù)下式進(jìn)行轉(zhuǎn)換[4]：

式中：μ為土體的泊松比。

根據(jù)地質(zhì)勘察成果，該工程壩殼料、壩基的物理力學(xué)參數(shù)見表1。

表1 壩殼料、壩基的物理力學(xué)參數(shù)表

公式法計算壩體沉降計算結(jié)果見表2。

表2 公式法計算壩體沉降

根據(jù)計算結(jié)果可知，最大壩高斷面壩體總沉降量為17.38 cm，壩基的總沉降量為6.58 cm，合計總沉降量為23.96 cm，沉降量與壩高的比值為0.32%。

3 有限元計算參數(shù)選取

靜三軸試驗固結(jié)圍壓分別為500、1 000和1 500 kPa，研究材料的應(yīng)力-應(yīng)變關(guān)系、體變-應(yīng)變關(guān)系，以確定有限元靜力計算鄧肯張E-B模型參數(shù)。試驗設(shè)備選用500和1 500 kN的高精度大型液壓靜動兩用三軸儀，見圖2。

圖2 高精度大型液壓靜動三軸儀

根據(jù)本工程材料靜力試驗結(jié)果，確定壩殼料、過渡料和壩基料的鄧肯張E-B模型參數(shù)見表3。

表3 靜力試驗鄧肯張E-B模型參數(shù)

4 二維有限元計算法

計算工況分為完建期工況和蓄水期工況，完建期只受重力作用，蓄水期受重力、水壓力等作用。選取壩體最大橫剖面作為計算斷面，計算模型網(wǎng)格劃分見圖3。

圖3 計算模型網(wǎng)格劃分圖

模型的邊界條件為：基巖底部施加固定位移約束，模型水平向上下游覆蓋層及基巖長度約等于1倍壩高，在上下游截斷面上，靜力計算只施加水平法向約束，允許其沉降。計算模擬施工的過程，荷載分級按壩體填筑次序進(jìn)行，第1步計算地應(yīng)力平衡；第2至第11步模擬壩體填筑，每層填筑土體厚度5.0 m；第12步模擬蓄水。二維有限元計算成果匯總見表4。

表4 二維有限元計算成果匯總

由計算結(jié)果可知，蓄水期的壩體最大水平位移向下游15.8 cm，發(fā)生在瀝青混凝土心墻墻高約1/2附近；最大垂直位移(沉降)為36.8 cm，約占壩高的0.49%，小于壩體高度的1%，發(fā)生在壩軸線2/3壩高附近。二維有限元法計算壩體沉降云圖見圖4、圖5。

圖4 完建期壩體沉降云圖

圖5 蓄水期壩體沉降云圖

5 三維有限元計算法

5.1 邊界條件及計算模型

三維有限元計算載荷步共分為48級，其中填筑分34級，蓄水分14級。大壩填筑完成之后開始蓄水，蓄至正常蓄水位1 894 m。瀝青混凝土心墻壩填筑采用分層激活單元法模擬，蓄水時水下部分施加浮托力，同時水壓力以面力的形式作用在瀝青混凝土心墻和基座上。其中，竣工期為第34步，滿蓄期為第48步。

壩體三維有限元網(wǎng)格見圖6，瀝青混凝土心墻有限元網(wǎng)格見圖7。三維有限元模型單元數(shù)為158 364個，節(jié)點數(shù)為174 954個。筑壩料和瀝青混凝土心墻料采用鄧肯張E-B模型，混凝土基座采用C30混凝土材料，ρ取值2.40 g/cm3，E取值30 GPa，υ取值0.20。心墻與過渡料、基座與心墻之間設(shè)置接觸面單元，計算參數(shù)參考同類工程[5-6]，見表5。

圖6 壩體三維有限元網(wǎng)格(單元：158 364個；節(jié)點：174 954個)

5.2 計算結(jié)果

通過大壩靜力有限元計算，典型斷面順河向和豎向位移分布見圖8-圖11。

圖7 心墻三維有限元網(wǎng)格(單元：5 949個；節(jié)點：11 944個)

表5 接觸面模型參數(shù)

圖8 大壩典型斷面竣工期順河向位移(cm，順河向為正)

圖9 大壩典型斷面滿蓄期順河向位移(cm，順河向為正)

圖11 大壩典型斷面滿蓄期豎向位移(cm，沉降為負(fù))

分析可知，竣工期大壩沉降最大值為22.10 cm；滿蓄期大壩沉降最大值為19.90 m?？⒐て?，順河向位移最大值為1.85 cm(向上游)和6.80 cm(向下游)；滿蓄期，由于水壓力作用，向上游變形區(qū)域和數(shù)值減小，向下游變形區(qū)域和數(shù)值增大，最大值為0.65 cm(向上游)和10.70 cm(向下游)。

6 工程類比

新疆水利水電勘測設(shè)計研究院近年來設(shè)計完成了數(shù)座碾壓式瀝青混凝土心墻壩，表6為幾座地質(zhì)條件、壩體結(jié)構(gòu)設(shè)計與本工程相近大壩的實際監(jiān)測資料。

表6 參考工程碾壓式瀝青混凝土心墻壩的豎向位移(磁環(huán)式沉降儀)

以上5座碾壓式瀝青混凝土心墻壩的沉降量與壩高的比值在0.18%～0.47%的范圍內(nèi)，本工程采用公式法、二維有限元法和三維有限元法計算得到的沉降量與壩高的比值分別為0.32%、0.49%和0.29%，計算得到的大壩沉降值分布規(guī)律與已建工程相近，符合工程經(jīng)驗，故認(rèn)為以上計算結(jié)果是相對準(zhǔn)確、合理的。在工程設(shè)計過程中，可根據(jù)3種計算方法的計算結(jié)果及已建工程實際，結(jié)合施工期現(xiàn)場監(jiān)測的結(jié)果，進(jìn)一步優(yōu)化設(shè)計。

7 結(jié) 語

本文以新疆某砂礫石土石壩為例，采用500、1 000 kN和1 500 kN高精度大型液壓靜動三軸儀測定了筑壩材料的靜力試驗鄧肯張E-B模型參數(shù)，分別采用公式法、二維有限元法和三維有限元法分析了土石壩的沉降變形特點，根據(jù)近年新疆已建的5座類似工程的實測數(shù)據(jù)驗證了計算結(jié)果的可靠性，為類似工程壩體沉降變形計算提供了借鑒。