楚金旺
(中國恩菲工程技術(shù)有限公司,北京 100038)
某尾礦壩的有限元計算分析
楚金旺
(中國恩菲工程技術(shù)有限公司,北京 100038)
基于壩體分區(qū)分層填筑和壩體材料非線性,采取材料分區(qū)和逐級施加荷載的方式建立尾礦壩有限元計算模型,并進(jìn)行計算分析。計算結(jié)果表明該尾礦壩的變形與應(yīng)力與尾礦壩的復(fù)雜地形地質(zhì)條件相適應(yīng),分布較為復(fù)雜;各分區(qū)與壩基之間的變形基本協(xié)調(diào),應(yīng)力從壩頂?shù)綁位饾u增大,應(yīng)力分布符合一般規(guī)律。剪應(yīng)力水平小于1.0,未發(fā)現(xiàn)剪切破壞單元,但壩基含礫粉質(zhì)粘土層的剪應(yīng)力水平相對較高。
尾礦壩; 鄧肯- 張模型; 非線性; 應(yīng)力; 位移; 剪應(yīng)力水平
某尾礦庫是我國一個大型礦山綜合開發(fā)利用項(xiàng)目的配套工程,主要包括尾礦壩、庫內(nèi)排洪設(shè)施、上游攔排洪設(shè)施和回水設(shè)施等。該尾礦壩總壩高135m,一期壩為堆石透水壩,壩頂標(biāo)高2 410m,壩高95m,二期壩為一期壩采用堆石料下游法加高30m,壩頂標(biāo)高2 440m,三期壩為上游法尾礦水力充填筑壩,最終壩頂標(biāo)高2 450m。
該尾礦壩一期壩、二期壩采用堆石料筑壩,三期壩采用尾礦水力充填法筑壩,堆石料和尾礦的變形不僅隨荷載的大小而變化,還與加荷的應(yīng)力路徑相關(guān),應(yīng)力應(yīng)變關(guān)系呈現(xiàn)明顯的非線性。鄧肯- 張非線性模型公式簡單,參數(shù)物理意義明確且容易獲取,為應(yīng)用最廣泛的模型之一。其中E—ν模型以切線彈性模量Et和切線泊松比νt作為計算參數(shù)[1~2],切線彈性模量表達(dá)式為:
(1)
切線泊松比為:
(2)
(3)
對于卸載情況,采用回彈模量Eur進(jìn)行計算:
(4)
式中:pa——單位大氣壓力;
S——剪應(yīng)力水平,反映材料強(qiáng)度發(fā)揮程度;
K、n、Rf、G、F、D、Kur、nur——模型參數(shù),可根據(jù)三軸試驗(yàn)資料整理而得[3~4]。
3.1 計算參數(shù)
根據(jù)尾礦壩壩址工程地質(zhì)、一期堆石壩、二期堆石壩和三期尾礦堆積壩填筑情況對壩體進(jìn)行分區(qū),由室內(nèi)三軸試驗(yàn)得到壩體、壩基的鄧肯E—v模型計算參數(shù),見表1。
表1 壩基壩料鄧肯- 張E—v模型計算參數(shù)
3.2 壩體填筑加載過程與網(wǎng)格剖分
根據(jù)材料分區(qū)和填筑過程,荷載采用逐級施加的方式。按照壩基、一期壩、一期壩內(nèi)尾礦、二期壩、二期壩內(nèi)尾礦、三期尾礦堆壩順序,將尾礦壩及壩基分為39級加載。加載過程模擬和網(wǎng)格剖分見圖1。
圖1 壩體填筑加載過程模擬
根據(jù)上述有限元計算原理、模型參數(shù)和加載過程,對該尾礦壩進(jìn)行了有限元計算分析,尾礦壩的位移、應(yīng)力和應(yīng)力水平的最大值見表2。
表2 尾礦壩有限元計算結(jié)果(最大值)
4.1 變形
正常水位工況下,尾礦壩向下游方向的水平位移最大值為75.0cm,出現(xiàn)在壩體下游2 335m高程與壩基相接處,向上游方向的水平位移最大值為24.0cm,出現(xiàn)在庫區(qū)2 400m高程附近尾粉質(zhì)粘土與尾粉砂相接處;豎直向最大沉降為290.3cm,出現(xiàn)在尾粉質(zhì)粘土庫區(qū)2 390m高程附近處。壩軸線剖面的水平位移最大值為55.0cm,位于河床壩段的2 340m高程處,豎向位移最大值為250.0cm,位于河床壩段的2 350m高程處。尾礦壩正常水位工況下變形等值線圖見圖2a、圖2b。
洪水工況下尾礦壩變形與正常水位工況相比,水位的上升使變形有所增大,但由于水位上升較低,數(shù)值變化甚小,其分布規(guī)律與正常水位時分布規(guī)律一致。向下游方向的水平位移最大值為75.0cm,出現(xiàn)在壩體下游2 335m高程與壩基相接處,向上游方向的水平位移最大值為24.5cm,出現(xiàn)在庫區(qū)2 400m高程附近尾粉質(zhì)粘土與尾粉砂相接處;豎直向最大沉降為291.1cm,出現(xiàn)在尾粉質(zhì)粘土庫區(qū)2 390m高程附近處。壩軸線剖面的位移分布規(guī)律與正常水位時期基本一致,由于上升水位較低,數(shù)值基本不變。尾礦壩洪水工況下變形等值線圖見圖2c、圖2d。
圖2 位移等值線圖(位移單位:cm)
圖3 主應(yīng)力等值線圖(主應(yīng)力單位:MPa)
4.2 應(yīng)力
正常水位工況下,尾礦壩第一主應(yīng)力最大值為2.54MPa,第三主應(yīng)力最大值為0.89MPa,均發(fā)生在二期壩體底部壩軸線位置附近。壩軸線剖面的第一主應(yīng)力極值為2.40MPa,第三主應(yīng)力極值為0.80MPa,均出現(xiàn)在河床壩基。尾礦壩正常水位工況下主應(yīng)力等值線圖見圖3a、圖3b。
洪水工況下,第一、第三主應(yīng)力分布規(guī)律與正常水位時期基本一致。第一主應(yīng)力最大值為2.54MPa,第三主應(yīng)力最大值為0.89MPa,均發(fā)生在二期壩體底部壩軸線位置附近。壩軸線剖面分布規(guī)律與正常水位時期基本相似,第一主應(yīng)力極值為2.40MPa,第三主應(yīng)力極值為0.80MPa,均出現(xiàn)在河床壩基。尾礦壩洪水工況下主應(yīng)力等值線圖見圖3c、圖3d。
4.3 剪應(yīng)力水平
正常水位和洪水位兩種工況下,剪應(yīng)力水平分布規(guī)律基本一致,最大值均約為0.84,位于壩體下游壩基2 310m高程處,尾礦庫區(qū)內(nèi)最大值為0.72,壩體內(nèi)未發(fā)生剪切破壞。兩種工況下,壩軸線剖面剪應(yīng)力水平分布規(guī)律基本一致,從壩頂?shù)綁位饾u增大,最大值約為0.80,未發(fā)生剪切破壞。尾礦壩剪應(yīng)力水平等值線圖見圖4。
圖4 典型剖面剪應(yīng)力水平等值線圖
(1)尾礦壩的變形與應(yīng)力分布較為復(fù)雜,這與尾礦壩的復(fù)雜地形地質(zhì)條件相適應(yīng),如上游壩坡由于尾礦庫的充填減小了堆石的泊松效應(yīng),使得大壩向上游的水平位移較小;由于尾礦的力學(xué)性質(zhì)參數(shù)較低,尾礦的沉降大于堆石壩本身的沉降值,最大變形位于尾礦庫內(nèi)。
(2)尾礦壩計算的最大沉降變形為297.4cm,其中堆石壩的最大沉降變形為255.0cm,分別占壩高的2.2%和1.8%,大壩各分區(qū)與壩基之間的變形基本協(xié)調(diào)。尾礦壩的應(yīng)力從壩頂?shù)綁位饾u增大,應(yīng)力分布符合一般規(guī)律。大主應(yīng)力最大值為2.55MPa,小主應(yīng)力最大值為0.90MPa,位于堆石壩底部的覆蓋層內(nèi)。
(3)尾礦壩的剪應(yīng)力水平均小于1.0,尚未發(fā)現(xiàn)剪切破壞單元。但是壩基特別是含礫粉質(zhì)粘土層的剪應(yīng)力水平相對較高,這與含礫粉質(zhì)粘土層天然孔隙率較高而力學(xué)參數(shù)較低有關(guān)。
[1] 王勖成,邵 敏.有限元法基本原理和數(shù)值方法[M].北京:清華大學(xué)出版社,1997.
[2] 錢家歡,殷宗澤.土工原理與計算(第二版)[M].北京:中國水利水電出版社,1996.
[3] 曹林衛(wèi),彭向和,李 德.變形場和滲流場耦合作用下的尾礦壩靜力穩(wěn)定性分析[J].重慶建筑大學(xué)學(xué)報,2007,29(5).
[4] 殷宗澤.高土石壩的應(yīng)力與變形[J].巖土工程學(xué)報,2009,31(1).
Finite element calculation and analysis of a tailings dam
The finite element calculation model of tailings dam was established and evaluated based on the dam fill in different partitions and layers and nonlinear materials by material partition and loading step by step. The calculation results show that the complex distribution of deformation and stress of the tailings dam adapts to complex topographic and geological conditions, the deformation between the partitions and the dam foundation is basically coordinated, the stress increases gradually from the top of the dam to dam foundation, and the stress distribution accords with the general law. The shear stress level is less than 1.0, and no shear failure unit has been found. However, the shear stress level of the gravel silty clay layer in the dam foundation is relatively high.
tailings dam; Duncan-Zhang model; nonlinear; stress; displacement; shear stress level
TD926.4+1
A
2017-- 06-- 14
楚金旺(1977-),男,山東禹城人,高級工程師,從事尾礦庫設(shè)計工作。
1672-- 609X(2017)04-- 0065-- 04