坝基渗透系数对均质土坝下游坝坡稳定影响.docVIP

坝基渗透系数对均质土坝下游坝坡稳定影响摘要：本文针对坝基渗透系数对均质土坝下游坝坡稳定性影响的问题,基于二维有限元法建立了相应的模型,运用GeoStudio岩土工程软件，模拟分析了在不同的面板渗透系数下浸润线的位置以及下游坝坡的稳定性,论证了防渗系统的有效性。 关键词：坝基渗透系数；均质土坝；有限元；下游坝坡稳定；GeoStudio软件 中图分类号：TV223 文献标识码： A文章编号：2095-2104（2012）01-0020-02 一 前言 我国已修建了8.6万余座水库，其中绝大部分水库都是采用土石坝。由于修建时期经济水平低、筑坝技术不成熟等原因，许多土石坝都存在工程质量差、病险隐患多和工程管理不善等严重问题，其中渗透变形是土石坝的主要病害。根据土石坝工程安全事故的统计，其中有37.1%的安全问题是由渗透变形引起的。因此在土石坝的研究领域，关于土石坝渗流性态的研究具有重要意义。目前研究土石坝渗流安全性态的方法有数值计算和实测资料分析两种，本文结合湖北某水库主坝工程，采用以实测资料分析为主、数值计算为辅的综合分析方法，对土石坝渗流安全性态进行了研究。本文采用有限元分析法，对土坝地基的透水层厚度对其抗倾覆稳定能力进行了计算研究。 二 工程概况 本文研究均质土坝的稳定渗流。拟定的相关参数如下：坝底黄河高程为500m，坝基厚度为20m，坝基宽为145m，坝体高为25m（即坝顶高程 525m）坝顶宽为10m，坝底宽为125m，上游水位20m（即黄河高程520m），下游水位10m。坝体渗透系数k=2.5m/d，坝基渗透系数有五种工况k1=0.50m/d、k2=0.55m/d、k3=0.60m/d、k4=0.65m/d、k5=0.70m/d。算得坝体稳定渗流场分布，从而分析坝基渗透系数对均质土坝下游坝坡稳定的影响。 三 有限单元理论计算 有限元法即是利用有限个单元的集合体代替连续的渗流场。单元间的结合点称为节点。选择简单的函数关系近似地表示单元上的水头分布，最后解得渗流场节点处满足一定精度的水头值。 首先对渗流场作剖分，即用一些假想的线将整个渗流场划分为有限个小区域。在划分单元时选择最常用的三角形、正方形单元，在渗流坡降变化较大的部位，即在坝体坝脚段加密网格，其他网格按从两端向中间渐变的方式划分。 渗流计算是在已知定解条件下解渗流基本方程，以求得渗流场水头分布和渗流量等渗流要素，它是工程设计的重要内容。有限单元法是解数学物理问题的一种数值解法。由于它所依据的理论具有普遍性.使其能成功地用于研究许多领域中的问题。就流体力学来说，它是最早又最成功地用于地下水渗流问题。有限单元法能系统地编成计算程序，很方便地处理复杂的边界条件，非均质土层，以及第二类流量边界不需专门处理而能自动满足，故目前在渗流计算中得到广泛地应用。 考虑土和水的压缩性，符合达西定律的二向非均质各向异性土体渗流. 其基本方程为 本文计算中不考虑土和水不可压缩，上式变为 即土坝稳定渗流的基本方程。式中：h是水头函数；x、z是空间坐标；kx、kz是以x、z轴为主轴方向的渗透系数；Ss是单位贮存量。 稳定渗流的有限单元法计算公式 [K]{h}={F} [K]总渗透矩阵；{F}是已知常数项，由已知节点水头得出。 四 软件在工程中应用 4.1 GeoStudio解决问题的基本流程 环境设置（工作空间的大小、坐标纸的定义）→创建有限元模型（创建或读入有限元模型）→属性定义（定义材料属性、渗透系数等）→网格划分（采用1m网格）→设置边界条件→分析求解→查看结果（查看分析结果、检查结果是否正确）。 4.2 创建有限元模型 根据已建立的模型，导入GeoStudio，并划分网格，如图1所示。 图 1计算网格图 4.3 分析解算。得到各种工况下浸润线的直观图。 1）坝基渗透系数为0.50m/d的渗流分布图，如图2所示。 图 2 渗流分布图 2）坝基渗透系数为0.55m/d的渗流分布图，如图3所示。 图 3 渗流分布图 3）坝基渗透系数为0.60m/d的渗流分布图，如图4所示。 图 4 渗流分布图 4）坝基渗透系数为0.65m/d的渗流分布图，如图5所示。 图 5 渗流分布图 5）坝基渗透系数为0.70m/d的渗流分布图，如图6所示。 图 6 渗流分布图 五 稳定计算分析 出口坡降是校核下游坝坡稳定安全的重要指标。 对于下部浸没段。 五种工况计算点为，下游水位水面点，因此。于是为一常数K。浸没段坡降为，。计算结果如表一所示。 表1各工况下渗出点坡降