完整版小型土石坝运用有限元渗流计算进行边坡稳定分析.docVIP

(完整版)小型土石坝运用有限元渗流计算进行边坡稳定分析

(完整版)小型土石坝运用有限元渗流计算进行边坡稳定分析

(完整版)小型土石坝运用有限元渗流计算进行边坡稳定分析

PAGEI

小型土石坝运用有限元渗流计算进行边坡稳定分析

摘要：安阳电厂丰安桥贮灰场基本充满,厂方拟新填子坝增加库容.另由于输送粉煤灰使用的是水力排放法，场区滞存大量积水，使坝体出现不同程度渗漏，威胁坝体安全,并造成坝后土地沼泽化。欲解决问题,先须进行渗流计算——是在已知定解条件下解渗流基本方程，以求得渗流场水头分布和渗流量等渗流要素，是工程设计的重要内容。电子计算机海量计算能力和数值计算方法的发展，特别是有限单元法提出后，推进了渗流数学模型的发展，为渗流计算提供了有效的方法。本次渗流计算的目的就是通过二维有限单元法进行渗流计算，分析渗流场的分布情况，通过对不同工况的计算，提出截渗方案.圆弧滑动的有限元法,以三角形单元的渗透力取代垂直条块周边的水压力，正确考虑了渗流方向及其力矩的影响，提高了计算精度。连续计算渗流与滑坡的统一程序使计算过程简便迅速。

关键词：渗流计算有限单元法浸润线圆弧滑动法

1有限单元法运用于渗流计算的目的和方法

土石坝水库失事事故中约1/4由渗流问题引起。安阳电厂贮灰场是在老河道一端筑坝而成,由于输送粉煤灰采用水力排放法,灰场滞存大量积水与原地下水形成较大的水头差,使坝体出现不同程度的渗漏。在这样的状况下使用场区粉煤灰为材料填筑子坝，就必须进行精确的渗流计算分析，因粉煤灰的物理力学特性决定了灰坝的安全与贮灰场的渗流分布，渗流问题解决不好，将导致不可估量的损失。

渗流计算的目的就是通过有限元渗流计算分析贮灰场渗流场的分布情况，得到的浸润线位置为坝体的边坡稳定及静、动力分析提供渗透压力，通过对不同工况的计算，提出截渗方案。

2渗流计算初始条件

2.1岩土工程勘察数据物理力学性质分析数据

2.2设定计算条件

水位条件

名称

子坝材料渗透系数(cm/s）

地基砂砾石层

4。05?10—2

初期坝

地基粘质粉土

3。70?10-5

上游水位:86m

坝体砂质粉土

2。29?10-5

下游地下水位：68m

原状粉煤灰

4.28?10—4

粉煤灰子坝

击实粉煤灰

1.21?10—4

上游水位：91m

碎石桩

5.00?10—1

下游地下水位：68m

排渗体

5。00?10—1

2.3设定计算工况

（SEQ编号一\r11）初期坝渗流计算，作为与实际情况校核之用。

（SEQ编号一2）加高工况渗流计算.①粉煤灰子坝基础加碎石桩对渗流场分布的影响。桩间距2m，桩长5m，桩径0。6m。②粉煤灰子坝基础不加碎石桩对渗流场分布的影响。③粉煤灰子坝基础加碎石桩、子坝下加排渗体时对渗流场分布的影响.间距2m，桩长5m,桩径0。6m。

3渗流计算的步骤与方法

3.1基本假设

（1)渗流服从达西定律；

（2)不考虑土体和水的压缩性，渗透时土体的孔隙大小和孔隙率不变；

(3）土体内水的饱和度不变。

渗流计算是在已知定解条件下解渗流基本方程，以求得渗流场水头分布和渗流量等渗流要素，但是无压渗流有渗流自由面(浸润线）,且非稳定渗流自由面随库水位升降而变动，加之一般渗流场有不同程度的非均质和各向异性,几何形状和边界条件较复杂，数学求解过程困难重重，仅能对一些简单流动情况获得解析解.然而，计算机海量计算能力配合有限元法处理复杂工程问题，尤其有关工程问题中连续体的力学问题，能解决过去解析方法难于解决的问题,使得渗流计算问题一下子迎刃而解,以渗流分析计算为例，有限单元法能系统地编成计算程序，很方便地处理复杂的边界条件，非均质土层，以及第二类流量边界不需专门处理而能自动满足，优势明显.

3.2一般步骤

将坝体有限边界区域按照材料分层概化模型，进行工况模拟计算，就可以形成渗流数学模型，即可进行有限元法一般5步骤计算：(1）物体的离散化；（2）挑选形函数或插值函数;（3）确定单元的性质；(4)组成物体的总性质方程组;(5）解方程组。本次算例采用的《PCSSL渗流有限元计算程序》就是较好的将有限元计算5步骤有效整合的程序，使用人员仅需要依据工程实际，填入参数即可得到贴近实际的计算结果。

渗流演算支配方程与定解条件

一般以水头h的分布为研究对象，对于垂直剖面上的二向渗流问题，按

（3—1）

渗流的基本微分方程式及其定解条件如下：

(3—2)

初始条件:h（x，z,0)=h0（x,z)

边界条件:①水头边界：h|=h（x,z,t)

②流量边界:k｜=-q（x，z,t)

3.3试算原则

渗流计算的基本及重要环节即