紅黏土塑性混凝土防滲墻與覆蓋層變形協(xié)調及應力變形分析

李續(xù)楠　楊具瑞　張正安

摘要利用ABAQUS有限元模擬軟件，對云南某黏土斜墻土石壩進行有限元分析，探究壩基中紅黏土塑性混凝土防滲墻對覆蓋層影響的變化規(guī)律。結果表明，在校核洪水位下，0.3、0.4和0.5 m防滲墻的鉛垂向位移與覆蓋層協(xié)調性較好；0.3 m厚防滲墻與覆蓋層協(xié)調性稍差于0.4、0.5 m厚防滲墻；墻頂垂直應力值為墻底垂直應力值的58%左右，分布較為均勻，有利于墻體的安全；大小主應力值均以壓應力為主，其中大主應力值在墻體與基巖連接處出現(xiàn)極小范圍的拉力區(qū)，但拉應力值遠小于抗拉強度0.19 MPa，小主應力最大值為-0.62 MPa，遠小于抗壓強度值1.69 MPa。

關鍵詞紅黏土塑性混凝土防滲墻；斜墻壩；有限元法；變形協(xié)調；滲流

中圖分類號TV640.31文獻標識碼

A文章編號0517-6611（2017）22-0162-05

AbstractABAQUS finite element simulation software was used to analyze the finite element analysis of a sloping core earth rock dam in Yunnan，and the influence of the plastic clay in the dam foundation was studied.The results showed that the vertical displacement of 0.3 ，0.4 and 0.5 m cutoff walls was better than that of the cover layer under the check flood level.The 0.3 m thick cutoff wall and the cover were slightly less than 0.4，0.5 m thick wall； the vertical stress value of the wall top was about 58% of the vertical stress value of the wall，the distribution was more uniform，which was favorable to the safety of the wall.The principal stress values were mainly compressive stress，The tensile stress was much smaller than the tensile strength of 0.19 MPa，and the maximum principal stress was -0.62 MPa，which was much smaller than the compressive strength value of 1.69 MPa.

Key wordsRed clay plastic concrete seepage wall；Inclinedwall dam；Finite element method；Deformation coordination；Seepage

水利工程是農(nóng)業(yè)的命脈，水庫的滲漏對農(nóng)田灌溉有著不可忽視的影響。1989年起，以塑性混凝土防滲墻為防滲體的工程已達到了60余座，均解決了水利工程中的各類問題，且運行情況良好。沈珠江等[1]通過對三峽二期圍堰中塑性混凝土防滲墻的有限元模擬，表明塑性混凝土防滲墻在復雜斷面覆蓋層中的位移和應力情況均在安全范圍內；于玉貞等[2]通過對比分析小浪底工程圍堰下擬建的塑性混凝土防滲墻和剛性混凝土防滲墻，表明塑性混凝土防滲墻具有較好的變形能力和較低的應力水平；宋博等[3]通過對不同荷載下的塑性混凝土滲透試驗，表明塑性混凝土應力水平在30%～80%時，滲透系數(shù)不會顯著增大，防滲效果良好；李創(chuàng)團[4]通過廣西鳳亭河水庫塑性混凝土防滲墻試驗的研究，表明塑性混凝土在只摻入黏土的情況下，可以滿足工程中抗壓強度及抗?jié)B性的要求，能夠適應土層周圍的變形，大大降低工程造價。基于國內對單獨摻入云南紅黏土的塑性混凝土應用較少、數(shù)值模擬較為匱乏，筆者結合云南省某斜墻土石壩，采用三維有限元法對蓄水期不同水位下[5]的紅黏土塑性混凝土防滲墻進行有限元模擬分析，對比分析紅黏土塑性混凝土防滲墻在校核洪水位、正常蓄水位及死水位下與兩側覆蓋層之間的變形協(xié)調性以及防滲效果。

1物理模型、模型參數(shù)及ABAQUS模型

1.1物理模型

云南省石林縣路美邑新壩水庫，其工程規(guī)模為小（2）型蓄水工程，樞紐工程等別為Ⅴ等。經(jīng)大壩安全評價鑒定路美邑新壩水庫大壩為三類壩。壩頂高程為1 734.66～1 735.41 m，此次復核高程為1 735.00 m，最大壩高12.50 m， 現(xiàn)狀壩頂長約220.00 m，壩頂寬5.00～8.00 m，此次復核寬為6.00 m，上下游壩坡坡均比為1∶2.0。紅黏土塑性混凝土防滲墻位于斜墻下端中部，厚度為0.3 m，墻頂插入斜墻0.5 m，墻底深入基巖0.5 m。該研究主要模擬防滲墻在校核洪水位、正常蓄水位、死水位下與覆蓋層之間的變形協(xié)調性及抗?jié)B能力。壩體及紅黏土塑性混凝土防滲墻，坐落于黏土礫石沖洪積層和二疊系上統(tǒng)峨眉山玄武巖上。壩體的典型橫、縱剖面圖如圖1和圖2所示。

1.2計算模型參數(shù)該研究中ABAQUS計算模型的參數(shù)均在云南淩屹工程設計有限公司針對該工程的協(xié)助下試驗得出，對模型合理分區(qū)后并考慮濕化效應的材料參數(shù)如表1所示。

1.3三維有限元計算模型

在數(shù)值模擬中，壩體及防滲墻的應力變形受到了覆蓋層、基巖以及兩岸山體的影響，為了使計算結果更加接近實際情況，得到的結果更加準確，在建立計算模型過程中，除了壩體之外還應包括兩岸山體、覆蓋層以及基巖。此次計算所應用的軟件為非線性功能強大的ABAQUS，采用Hypermesh軟件建立物理模型以及劃分網(wǎng)格，導入ABAQUS/standard模塊中完成計算。計算網(wǎng)格單元總個數(shù)140 312個，單元警告數(shù)203個，網(wǎng)格警告率為0.14%，網(wǎng)格整體質量較好，模型整體三維有限元網(wǎng)格如圖3所示。

1.4關鍵部位網(wǎng)格、分級填筑、地應力平衡及流固耦合

1.4.1關鍵部位網(wǎng)格質量。

網(wǎng)格的質量對計算結果有著很大的影響，該研究中的計算網(wǎng)格單元數(shù)量140 312，其中在關鍵部位，即斜墻下的紅黏土塑性混凝土防滲墻及其周圍覆蓋層進行網(wǎng)格的加密，在加密過程中，通過初次計算以及再次加密網(wǎng)格密度來比較關鍵部位計算結果誤差，初步判斷網(wǎng)格密度是否合適[6]，前后兩次垂直應力計算結果誤差為3%，小于規(guī)范誤差最小值5%，可以初步認定網(wǎng)格的精度可達到基本要求。

1.4.2地應力平衡。

在壩體填筑前，由于覆蓋層、基巖的自重產(chǎn)生的應力[7]，往往要考慮覆蓋層、基巖對紅黏土塑性混凝土防滲墻以及壩體的影響，同時又要求施加地應力后，填筑壩體前的覆蓋層以及基巖不產(chǎn)生變形，因此需要對覆蓋層以及基巖進行地應力的平衡。覆蓋層以及基巖地應力平衡后，在保證其各應力基本不變的情況下，位移量由自重情況下的厘米量級減小至10-6 m量級，接近于零，可見地應力平衡已滿足要求。地應力平衡后的壩體位移如圖4所示。

1.4.3施工蓄水的模擬。

在三維有限元計算分析中，將整個模型分為11個分析步，控制各層單元的生死，通過分級加載的計算方式，實現(xiàn)對工程分級填筑的模擬[8]，施工過程中壩基覆蓋層、防滲墻先于壩體填筑，斜墻與壩體填料同時分級填筑上升，在模擬施工完成后，分別模擬蓄水至校核洪水位、正常蓄水位及死水位。具體模擬過程如表2所示。

1.4.4ABAQUS中流固耦合的實現(xiàn)。

目前土石壩的有限元模擬分析，往往單獨進行應力場和滲流場的模擬，雖然簡單易行，但和實際工程不符，經(jīng)前人研究得出，單獨考慮應力場和滲流場的情況下，所得計算結果與應力場、滲流場耦合作用下的計算結果有明顯差異。該研究在計算模擬中考慮應力場和滲流場的相互作用下，通過選擇單元類型、改變孔隙率進而改變滲透系數(shù)的方法，實現(xiàn)應力場和滲流場的耦合分析。

2本構模型

紅黏土塑性混凝土是相對剛性混凝土摻入了大量的云南本地紅黏土，使得紅黏土塑性混凝土與壩體填土和覆蓋層等有著相似的本構關系和變形能力。目前在水利工程中主要采用鄧肯張E-B模型實現(xiàn)對土石壩的模擬。相比摩爾庫倫等其他本構模型，鄧肯張E-B模型能夠更加準確地反映出應

力變形情況。在參數(shù)的選取上，鄧肯張E-B模型有著成熟經(jīng)

驗且確定較為方便，故此次模擬中選取鄧肯張E-B模型作為土石壩體及防滲墻的本構模型，其本構模型公式如下[9]：

Et=Kpaσ3pan1-Rf（σ1-σ3）（1-sinφ）2c cosφ+2σ3 sin φ2 （1）

Eur=Kurpaσ3panur （2）

Bt=Kbpaσ3pam（3）

式中，Et為切線變形模量；K為切線模量基數(shù)；n為切模量指數(shù)；c為凝聚力；σ1為土體的大主應力；σ3為土體的小主應力；φ為內摩擦角；Rf為破壞比；Eur為切線彈性模量；Kur、nur為由試驗確定的系數(shù)，確定方法與K、n相似；Bt為體積彈性模量系數(shù)；Kb為K的基數(shù)；m為體積模量指數(shù)；pa為大氣壓。

3防滲墻與覆蓋層變形協(xié)調性分析

為了便于對比分析紅黏土塑性混凝土防滲墻在不同工況下的變形情況，該研究選取壩體縱剖面進行應力變形的對比分析。該研究主要研究上游水庫分別蓄水至校核洪水位、正常蓄水位及死水位時，紅黏土塑性混凝土防滲墻與兩側覆蓋層的變形協(xié)調性和抗?jié)B能力。此次計算模擬中，規(guī)定順河向位移向下游為正，鉛錘向位移以向上為正，各個主應力以拉為正，反之為負。

3.1防滲墻與覆蓋層鉛垂向協(xié)調

由圖5可見，

在校核洪水位下，防滲墻厚度分別為0.3、0.4和0.5 m時，墻頂鉛錘向位移值分別為-2.22、-2.10和-2.05 cm，墻底鉛錘向位移值分別為-0.68、-0.69和0.69 cm，位移值隨著高程的增加而增大，防滲墻的最大位移值均發(fā)生在墻頂，最小值發(fā)生在墻底。上、下游覆蓋層整體位移值與防滲墻較為接近[10]，協(xié)調性較好，由于上游壩址的拖拽，在墻頂附近協(xié)調性稍差，但在合理范圍內。隨著防滲墻厚度的增加，同一高程點的防滲墻位移值逐漸減小，且隨著高程的增加，減小的趨勢越明顯。

3.2防滲墻與覆蓋層順河向協(xié)調

由圖6可知，在校核洪水位下，防滲墻及覆蓋層的順河向位移值最大值發(fā)生在1 718.5 m附近，其中0.3、0.4和0.5 m防滲墻順河向位移最大值分別為-0.54、-0.54和-0.53 cm，最小值均發(fā)生在墻頂附近，分別為-0.10、-0.13和-0.14 cm。其中防滲墻中部位移均大于墻頂和墻底的位移值。隨著墻體厚度的增加，此時防滲墻變形能力變差，同一高程點的順河向位移值隨之減小，但在墻頂附近，隨著墻體厚度的增加，在上游壩趾的作用下，墻體頂端受到拖拽作用減小，故順河向位移值有所增大。

3.3防滲墻變形及應力

圖7為校核洪水位下，防滲墻的鉛垂向沉降值與順河向位移計算結果。由圖8可見，在校核洪水位下作用，防滲墻厚度分別為0.3、0.4和0.5 m時，墻體內部的垂直應力最大值分別為-0.61、-0.57和-0.56 MPa，均發(fā)生在墻底，最小值分別為-0.37、-0.34和-0.32 MPa，均發(fā)生在墻頂，墻頂垂直應力值分別為墻底的60%、58%和57%。墻體頂部垂直應力為墻底的50%～60%，垂向應力值分布較為均勻[11]，與傳統(tǒng)混凝土小于10%范圍相比，有著明顯區(qū)別。其原因是傳統(tǒng)的剛性混凝土防滲墻變形能力較差，與兩側覆蓋層產(chǎn)生巨大沉降差，如小浪底工程中的防滲墻與覆蓋層沉降差高達28.2 cm，在防滲墻和墻體之間產(chǎn)生巨大的拖拽力，導致防滲墻底部垂向應力值遠大于頂部應力值。當防滲墻在不同水位作用下，墻體在校核洪水位、正常蓄水位及死水位下的垂直應力值隨著水位的下降而減小。在相同水位作用下，隨著防滲墻墻體厚度的增加，垂直應力值也隨之減小，表明厚度的增加有利于減小斜墻下部的紅黏土塑性混凝土防滲墻垂直應力。校核洪水位下，防滲墻厚度分別為0.3、0.4和0.5 m時，大主應力的最大值分別為-0.06、-0.62和-0.77 MPa，大主應力只在防滲墻與兩端的極小區(qū)域內出現(xiàn)拉應力，但僅有0.04 MPa，遠遠低于該研究工程中的紅黏土塑性混凝土抗拉破壞值0.19 MPa。在防滲墻縱剖面的上的小主應力，同樣以壓應力為主，最大值均發(fā)生在墻底。在校核洪水位下，防滲墻厚度分別為0.3、0.4和0.5 m時，小主應力的最大值分別為-0.58、-0.54和-0.52 MPa，最大值均為達到該研究工程中的抗壓強度破壞值-1.69 MPa，故防滲墻是安全的[12]。

4結論

在校核洪水位作用下，鉛垂向位移值隨著墻體厚度的增加，協(xié)調性變差；鉛垂向位移值隨著墻體厚度的增加，協(xié)調性變好；墻頂與墻底垂直應力比值隨著墻體厚度的增加，有減小的趨勢，表明墻體厚度的增加不利于墻體的安全；0.3、0.4和0.5 m的大小主應力值均為超過破壞值，故防滲墻是安全的。

參考文獻

[1]

沈珠江，劉松濤.三峽二期高土石圍堰應力應變分析研究[J].人民長江，1996，27（10）：1-4.

[2] 于玉貞，濮家騮，劉鳳德.土石壩基礎塑性混凝土防滲墻材料力學特性研究[J].水利學報，1995（8）：21-27.

[3] 宋博，何江濤，劉長禮.不同荷載作用下塑性混凝土滲透性試驗研究[J].南水北調與水利科技，2015，13（2）：309-313.

[4] 李創(chuàng)團.鳳亭河水庫主壩塑性混凝土防滲墻試驗研究[J].防滲技術，1999，5（4）：13-15.

[5] 滕彥磊.塑性混凝土防滲墻應力變形有限元分析[D].鄭州：鄭州大學，2010.

[6] 費康，張建偉.ABAQUS在巖土工程中的應用[M].北京： 中國水利水電出版社，2013.

[7] 代汝林，李忠芳，王姣.基于ABAQUS的初始地應力平衡方法研究[J].重慶工商大學學報（自然科學版），2012，29（9）：76-81.

[8] 何芳嬋.土石壩施工過程應力變形仿真分析[D].鄭州：鄭州大學，2007.

[9] DUNCAN J M，CHANG C Y.Nonlinear analysis of stress and strain in soils[J].Asce soil mechanics & foundation division journal，1970，96（5）：1629-1653.

[10] 鄧明基.東平湖圍壩塑性混凝土防滲墻性能研究[D].北京：清華大學，2005.

[11] 王清友，王綦正.塑性混凝土防滲墻結構非線性分析及其設計[J].水力發(fā)電，1992（8）：18-23.

[12] 司政，陳堯隆，李守義.土石壩壩基塑性混凝土防滲墻應力變形分析[J].水力發(fā)電，2008，34（2）：32-35.