張惠軍

(上?？睖y設(shè)計研究院，上海 200434)

在高面板堆石壩工程實際施工中，由于堆石料自身特性以及施工期進(jìn)度與度汛的需要，一般采用分期填筑的方式完成壩體和面板的施工。面板堆石壩彈塑性有限元分析計算體現(xiàn)了堆石料復(fù)雜變形的特點，逐級加載方法反映了材料剛度在各加載階段的變化，較好地模擬了壩體施工過程。

1 工程概況

水布埡樞紐工程位于湖北省巴東縣境內(nèi)的清江干流上，屬一等大(1)型工程。水庫正常蓄水位高程為400m，相應(yīng)庫容43.12億m3，校核洪水位高程為404m，水庫總庫容為45.8億m3，具有突出的多年調(diào)蓄能力。水布埡面板堆石壩［1］作為樞紐工程的主體建筑物之一，其最大壩高為233m，壩頂長度為674.66m，壩頂寬為12m。大壩壩前坡度1∶1.4，壩后平均坡度1∶1.46。壩體斷面從上游至下游依次分為蓋重區(qū)、鋪蓋區(qū)、墊層區(qū)、過渡區(qū)、主堆石區(qū)、次堆石區(qū)和下游堆石區(qū)等7個主要填筑區(qū)，面板沿壩軸線方向分縫分塊為寬厚不等的58塊，其中受壓區(qū)面板寬為16.0m，而受拉區(qū)面板寬為8.0m，板厚為0.3m～1.1m不等。趾板與基巖間設(shè)有錨筋連接，結(jié)構(gòu)形式為壩前設(shè)標(biāo)準(zhǔn)板，下接防滲板。

2 面板堆石壩三維彈塑性有限元計算

2.1 有限元模型

以水布埡面板堆石壩為實例建立的三維有限元網(wǎng)格［2］共包括8090個單元，8292個節(jié)點。網(wǎng)格單元主要采用八節(jié)點六面體等參單元計算，為適應(yīng)基巖開挖面的不規(guī)則變化，壩體與基巖結(jié)合部由退化的五面體和四面體單元填充。有限元網(wǎng)格如圖1所示。

面板為鋼筋混凝土，密度ρ=2.4t/m3，彈性模量E=25.5GPa，泊松比υ=0.167，面板單元采用線彈性本構(gòu)模型計算。壩體材料計算參數(shù)［3］見表 1。

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

2.2 堆石料本構(gòu)關(guān)系

堆石料作為面板堆石壩主要的填筑材料，其力學(xué)性質(zhì)決定了壩體的基本工作狀態(tài)，亦對防滲面板和周邊縫止水產(chǎn)生重要影響。由于堆石料等粗粒土具有非線性變形、剪脹等特性，選擇簡單可行，貼近實際的本構(gòu)模型，因此具有重要的意義。土體本構(gòu)彈塑性模型國內(nèi)工程中較常用的有沈珠江院士的南水雙屈服面模型［4-6］，采用雙屈服面假定，既反映了堆石料的剪切變形及壓縮變形機理，同時又能較好適用于復(fù)雜的應(yīng)力路徑和應(yīng)力狀態(tài)。

2.3 分級加載

面板堆石壩三維有限元計算荷載共分為21級逐級加載，為了模擬施工填筑及蓄水過程，其中壩體施工填筑分為18級，面板鋪裝分為1級，蓄水分為2級加載。分級加載過程如表2所示。

表2 分級加載

面板堆石壩壩體單元按施工分期填筑厚度來逐層加載，待壩體填筑完畢后，再加載面板單元，開始蓄水分兩級加載。堆石料本構(gòu)關(guān)系采用南水模型，設(shè)置薄層單元模擬堆石體和面板間的相互作用，采用中點增量法的三維彈塑性有限元分析來模擬施工和蓄水過程中的荷載逐級遞增。

3 計算結(jié)果分析

通過對計算結(jié)果分析，為了反映壩體竣工期和蓄水期的變形和應(yīng)力分布情況，截取壩體軸線中心附近處典型剖面來做比較。

3.1 壩體變形分析

竣工期和蓄水期的壩體水平變形見圖2，豎直變形見圖3。可從對比圖中看到，位移等值線大體分布規(guī)律相同，蓄水期的水平位移和豎直位移都較竣工期有幅度較小的增大。蓄水后壩體豎直變形最大為1.843m，發(fā)生在壩體中軸線最大典型剖面附近處，壩體向上游的最大位移為0.232m，向下游位移為0.194m。

3.2 壩體應(yīng)力分析

竣工期和蓄水期的壩體第一主應(yīng)力等值線見圖4，第三主應(yīng)力等值線見圖5。從對比圖中可見兩者分布規(guī)律一致，計算結(jié)果相近，竣工期壩體的第一和第三主應(yīng)力值分別為3.11MPa和1.41MPa，蓄水期壩體的第一和第三主應(yīng)力值分別為3.56MPa和1.60MPa。主應(yīng)力最大值均出現(xiàn)在壩體中部偏上游處。

4 結(jié)語

堆石料作為面板堆石壩壩體填筑的主要材料，具有明顯的非線性和剪脹等性質(zhì)，雙屈服面模型能較好地體現(xiàn)這種受力特性。逐級加載方法能夠很好地模擬壩體結(jié)構(gòu)本身隨施工過程的變化，反映出施工過程中各階段的應(yīng)力和應(yīng)變情況。但計算結(jié)果與實際觀測結(jié)果相比有不小的差別［1］。1)蓄水期水平位移和豎直變形相較竣工期有增大。由于蓄水后靜水壓垂直作用在面板上并且傳遞到堆石體，壩體部分區(qū)域水平位移等值線趨向大壩下游，豎直變形等值線有向大壩上游移動趨勢。2)在竣工期和蓄水期兩種工況下，同一剖面壩體主應(yīng)力分布和變化規(guī)律基本相同，結(jié)果趨近一致，蓄水期同一點的主應(yīng)力較竣工期增大明顯，且應(yīng)力隨著壩體升高逐漸減小，其最大主應(yīng)力均位于壩底中部。

［1］楊啟貴，劉 寧，孫 役，等.水布埡面板堆石壩筑壩技術(shù)［M］.北京:中國水利水電出版社，2010.

［2］張惠軍.機群環(huán)境下面板堆石壩施工加載有限元并行計算方法研究［D］.大連:大連理工大學(xué)建設(shè)工程學(xué)部，2011.

［3］王柏樂.中國當(dāng)代土石壩工程［M］.北京:中國水利水電出版社，2004.

［4］朱百里，沈珠江.計算土力學(xué)［M］.上海:上?？茖W(xué)技術(shù)出版社，1990:289-306.

［5］章為民，沈珠江.混凝土面板堆石壩三維彈塑性有限元分析［J］.水利學(xué)報，1992(4):75-78.

［6］錢家歡，殷宗澤.土工原理與計算［M］.第2版.北京:中國水利水電出版社，1996.