基于流變特性的面板砂礫石壩的有限元分析

付勇智

(黑河市水務(wù)局，黑龍江 黑河 164300)

?

基于流變特性的面板砂礫石壩的有限元分析

付勇智

(黑河市水務(wù)局，黑龍江 黑河 164300)

面板壩具有完整型好，適應(yīng)性強，抗震能力突出等特點，目前在國內(nèi)外大壩建設(shè)中占有很大比例。隨著越來越多的工程中應(yīng)用砂礫石和軟巖等材料，適宜的應(yīng)力變形問題對大壩的安全施工、運行具有及其重要的意義。文章研究面板砂礫石壩在滲透變形的條件下，對大壩進(jìn)行三維有限元分析，得出壩體在竣工期和蓄水期的不同工況下的位移和應(yīng)力變化。通過計算結(jié)果為高面板砂礫石壩的流變效應(yīng)和抗震設(shè)計提供必要的理論依據(jù)和工程參考數(shù)值，所獲取的計算成果對類似工程的設(shè)計具有一定的參考應(yīng)用價值。

流變特性；面板；砂礫石壩；有限元；分析

0 前 言

混凝土面板堆石壩是指以堆石料或者砂礫石料作為壩體主要填筑材料，防滲體系由面板、防浪墻、趾板和防滲帷幕構(gòu)成的一完整防滲體系。和傳統(tǒng)的土石壩相比，面板壩的施工方法簡便、安全系數(shù)高和適應(yīng)各種地質(zhì)等優(yōu)點[1]。在面板砂礫石壩的應(yīng)力變形有限元分析中，采用非線性彈性方法對于壩體高度較低時可以得到較為合理的結(jié)論。文章在考慮砂礫石流變作用影響的基礎(chǔ)上，準(zhǔn)確的進(jìn)行有限元靜力模擬，對于提高面板砂礫石壩設(shè)計的準(zhǔn)確性具有一定的現(xiàn)實意義[2-3]。

1 面板砂礫石壩的靜力有限元分析方法

1.1 砂礫石料靜力計算本構(gòu)模型

文章采用較廣泛的鄧肯張E-B模型模擬材料進(jìn)行計算，該模型的基本原理如下[4-5]：

在初始加載中，切線彈性模量Et為：

Et=Ei(1-RfSl)2

(1)

式中：Ei為初始彈性模量：

(2)

式中：σ3為小主應(yīng)力。

(3)

式中：Rf破壞比，(σ1-σ3)f根據(jù)摩爾庫倫準(zhǔn)則確定，定義為砂礫石料強度。

(4)

St是主應(yīng)力差在實際和破壞時的比值，反應(yīng)出應(yīng)力水平：

Sl=(σ1-σ3)/(σ1-σ3)f

(5)

卸荷模量Eur為：

Eur=KurPa(σ3-Pa)nur

(6)

式中：nur為卸荷后加荷時的彈模指數(shù)，Kur為彈?；鶖?shù)。

(7)

式中：Bt為切線體積模量；m為初始體積模量指數(shù)；Kb為初始體積模量基數(shù)。

加卸荷函數(shù)的確定表達(dá)式：

(8)

(9)

由公式可以推導(dǎo)在鄧肯張E-B模型中，通過三軸試驗確定K、Kb、Kur、c、φ、m、n和Rf，將各參數(shù)代入(6)和(7)確定Bt和Eur，進(jìn)而得出剛度矩陣[D],最終計算應(yīng)力應(yīng)變的增量。

1.2 混凝土面板和趾板的本構(gòu)模型

國內(nèi)外學(xué)者已通過大量工程實例得出：當(dāng)混凝土材料所受荷載導(dǎo)致的變形未達(dá)到破壞狀態(tài)時，有限元分析采用彈塑性模型或是線彈性模型相差不大[6]。因為面板混凝土材料和趾板一般是線彈性的，但在工程實際中常常是三向應(yīng)力加載，故認(rèn)定變形也是非線性的。大量工程試驗表明，當(dāng)所用材料抗壓強度是混凝土材料的30%-50%時，彈性變形居主導(dǎo)地位，應(yīng)力變形曲線呈現(xiàn)線性變化。

1.3 面板與墊層接觸面以及面板接縫的有限元模型

目前接觸面常用的接觸單元分為兩類：無厚度單元和薄層單元[7]。面板砂礫石壩是填筑要求分為多個分區(qū)，各分區(qū)材料區(qū)分，性能區(qū)別。當(dāng)大壩蓄水后，水壓力現(xiàn)施加面板上，然后傳遞給壩體，進(jìn)而形成一定的擠壓或開裂。在接觸面兩種不同特性的材料會出現(xiàn)變形不一致。為了更貼切的表明不同性質(zhì)的材料間相互作用，在接觸面需要設(shè)置接觸單元，文章采用無厚度單元。

2 工程實例

2.1 工程概況

黑泉水電樞紐站位于青海省湟水支流的寶庫河上，距里青海省會西寧僅75km，流域面積1044km2，水電站裝機3×4MW，總庫容1.82億m3，多年平均氣溫2.8℃，是兼具發(fā)電和調(diào)洪功效的大(2) 型水利樞紐工程[8]。壩高123.5m，正常蓄水位2887.5m，壩頂長438m?？偺钪?44萬m3，砂礫石料為421萬m3。防浪墻高4.7m，壩頂高程2894.5m，，上游壩坡比1∶1.55，下游設(shè)三級馬道，黑泉水利樞紐由大壩、右岸溢洪道、左岸發(fā)電洞、導(dǎo)流洞、電站構(gòu)成[9]。黑泉面板砂礫石壩標(biāo)準(zhǔn)剖面圖見圖1。

2.2 有限元模型

壩基設(shè)置在壩體與基巖接觸面，上游設(shè)置到趾板，下游為1.5倍壩體高度，左右岸設(shè)置在與基巖接觸面。壩基約束采用固定鉸形式，兩岸與基巖的接觸面采用法向約束形式。整體模型主要為六面體八節(jié)點單元，局部設(shè)置為四面體四節(jié)點單元。得到剖分后壩體單元5792個，接觸面單元312個，節(jié)點總數(shù)6708個[10]。根據(jù)材料不同性質(zhì)共設(shè)置6個區(qū)進(jìn)行分析，面—面接觸單元設(shè)置在混凝土面板與墊層料處。三維有限元模型如圖2。

圖1 黑泉面板砂礫石壩標(biāo)準(zhǔn)剖面圖

圖2 黑泉面板砂礫石壩三維有限元模型

2.3 計算工況

文章采用兩種工況來分析面板砂礫石壩的變形與應(yīng)力變化規(guī)律。在竣工期，面板砂礫石壩按設(shè)計高度填筑到2894.5高程時壩體應(yīng)力變形規(guī)律，所受荷載僅為大壩自重。在蓄水期模擬中，所受靜荷載(不考慮動荷載浪壓力等)為上游水壓力和壩體自重，考慮流透變形的情況下，模擬出面板砂礫石壩在水荷載作用下變形和應(yīng)力變化規(guī)律。

2.4 計算結(jié)果

本工程實例靜力和流變部分的應(yīng)力應(yīng)變等值線圖均為標(biāo)準(zhǔn)斷面的結(jié)果，其符號規(guī)定：變形結(jié)果中，垂向位移以向下為負(fù)；堆石體橫斷面上的位移以指向下游為正，向上游為負(fù)；壩體和面板的應(yīng)力結(jié)果中以受壓為負(fù)，以受拉為正。

圖3 竣工期壩體橫斷面豎向位移等值線圖

圖4 竣工期壩體橫斷面水平向位移等值線圖

圖5 蓄水期壩體橫斷面豎向位移等值線圖

圖6 蓄水期壩體的橫斷面水平向位移等值線圖

從以上結(jié)果可以得出：大壩竣工時的最大豎向位移0.37m，最大位移處于中部偏下處，蓄水期在水壓力作用下，最大豎向值為0.39m，增大了約6.5%，在竣工期壩體的上游位移為0.126m，下游位移為0.11m，幾乎都位于壩體偏下游位置。蓄水期的上游位移0.063m，下游位移0.126m，圖6可以得出上游位移與竣工期相比變化較大，且等值線圖變化分布整體偏向下游。

圖7 竣工期壩體橫斷面大主應(yīng)力等值線圖

圖8 竣工期壩體橫斷面小主應(yīng)力等值線圖

圖9 蓄水期壩體橫斷面大主應(yīng)力等值線圖

圖10 蓄水期壩體橫斷面小主應(yīng)力等值線圖

以上應(yīng)力變化曲線圖的有限元計算結(jié)果可以得出：竣工期的大壩大主應(yīng)力為2.62MPa，蓄水后為2.77 MPa，增大了約5.7%；竣工期的小主應(yīng)力為1.08 MPa，蓄水期為1.22 MPa，增大約13.4%；不同工況下的分析結(jié)果顯示蓄水后壩體的主應(yīng)力均有不同程度增加，但分布規(guī)律未發(fā)生較大變化，峰值都位于壩體下部壩軸線附近，小主應(yīng)力的分布在壩體中部比較少；比較明顯的是由于水壓力，上游局部壩坡的應(yīng)力值均有所增加。

3 結(jié) 論

經(jīng)過對工程實例在不同工況的三維有限元分析，得出以下結(jié)論：砂礫石料是比較好的筑壩材料，特點是壓縮性低，變形模量高，開采難度小等優(yōu)點，在工程的實際應(yīng)用中，可以著重考慮，但對于防滲體系的建立問題不容忽視。本工程在考慮砂礫石的滲流作用下，通過對非線性模型鄧肯張E-B進(jìn)行二次開發(fā)，通過大型有限元數(shù)值仿真軟件ADINA來模擬大壩在不同工況下的應(yīng)力變形規(guī)律研究，計算結(jié)果較為合理的模擬出了大壩實際工程施工中和蓄水后的變形規(guī)律。為壩體在應(yīng)力變形方面的控制設(shè)計提供更為合理的參數(shù)，同時也可為類似工程興建和安全運行提供一定的參考。

[1]蔣國澄,傅志安,鳳家驥.混凝土面板壩工程[M].武漢:湖北科學(xué)技術(shù)出版社,1997：35-43.

[2]劉宏,韓文喜,張倬元.砂礫石土料的壓實特性[J].三峽大學(xué)學(xué)報：自然科學(xué)版,2002,24(04):297-299.

[3]顧淦臣,張振國.鋼筋混凝土面板堆石壩三維非線性有限元動力分析[J].水力發(fā)電學(xué)報,1988(01):26-45

[4]岑威鈞.堆石料亞塑性本構(gòu)模型及面板堆石壩數(shù)值分析[D].南京:河海大學(xué),2005.

[5]中華人民共和國水利部.GB50218—94工程巖體分級標(biāo)準(zhǔn)[S].北京：中國計劃出版社，1995.

[6]羅先啟，葛修潤.混凝土面板堆石壩應(yīng)力應(yīng)變分析方法研究[M].北京：中國水利水電出版社，2007：32-67.

[7]沈鳳生,陳慧遠(yuǎn),潘家錚.混凝土面板堆石壩的蓄水變形分析[J].巖土工程學(xué)報,1990,12(01) :74-81.

[8]何輝.面板堆石壩應(yīng)力變形特性分析與程序設(shè)計[D].天津:天津大學(xué),2005.

[9]青海省水利水電勘測設(shè)計研究院.黑泉混凝土面板砂礫石壩的設(shè)計[J].水利水電技術(shù),2000,31(01):71-73.

[10]蔣國澄.中國混凝土面板堆石壩20年[M].北京:中國水利水電出版社,2005:542.

Finite Element Analysis of Facing Sandy-gravel Dam based on Rheologic Characteristics

FU Yong-zhi

(Heihe Rural Water Affairs Bureau, Heihe 164300, China)

Deck dam owns the characteristics of complete dam type, high suitability and outstanding seismic capacity and accounts for large proportion in domestic and foreign dam construction. With the material application of sandy gravel and soft rock for more and more projects, suitable stress deformation will be meaningful for safe construction and operation of dam. In this paper, three-dimensional finite element was analyzed for deck sandy-gravel dam under the condition of seepage and deformation, drawing the conclusion of displacement and stress change in various project situations of completion period and water storage period. To supply necessary theory accordance and project reference value for rheologic effect and seismic design of high deck sandy-gravel dam shows that the calculation results obtained have a certain applying value for similar project designs.

rheological characteristics;face slab;sandy-rock dam; finite element; analysis

1007-7596(2017)05-0019-04

2017-04-12

付勇智(1982-)，男，黑龍江黑河人，工程師，從事水資源管理工作。

TV641.4

B