基于有限元的堤坝渗流稳定分析论文.pdf

第47卷第11期山西建筑Voi.47No.11

・25・2021年6月SHANXIARCHITECTUREJu/.2021

DOI：1.13719//cn/b1009-6525.2021.11.057

基于有限元的堤坝渗流稳定分析

刘忠在

（内蒙古大学公路工程试验检测中心，内蒙古呼和浩特010070）

摘要:为了研究干旱区重力式土坝的稳定性，基于有限元软件,从堤坝稳态渗流角度，并考虑重力和渗流的影响，进行了水位线、

渗透系数和土的种类等因素对堤坝稳定性的影响分析,为干旱区相类似的工程提供一定的参考。

关键词:堤坝，渗流，稳定性，流固耦合

中图分类号:TV139.11文献标识码：A文章编号：1029-6525（2221）11-2155-23

堤坝的稳定性关系到整个水利工程的安全，早在20世海拔较高,南边是地形开阔的平原区。整体地势北高南低。

纪，许多学者就对坝身内部渗流进行了大量研究,法国工程黄河流经于此，有了黄河流域的锦上添花，为该处带来了丰

师达西⑴利用砂土进行透水试验,得到了达西定律,成为渗富的水资源。

流发展的一个标志性阶段。通过对当地土类土质等进行调查,结合实际情况得到

在渗流分析的过程中，首先要考虑的是浸润面的位置下列数据,坝体材料渗透系数为化二々二1x12一m/d,y二

问题,只有确定了水位自由面才可以进行渗流计算,现在对19.1kN/m3，在浸润面之上土体抗剪强度为e二21.5。,二

于浸润面的确定方法主要有:数值解析法、数值模拟法和饱57.3kPa,在浸润面之下土体抗剪强度为©二15.2°,=

和一非饱和渗流模型方法等一些其他方法2。王媛⑶为解36.54kPao

决初流量法解不稳定,率先提出了区域识别函数,对其解进2数值建模

行调整，提出了改进的初流量法。彭华等⑷利用全新的虚通过现场调查,根据工程概况,建立三维有限元模型分

单元法进行计算，增加了计算的速度,减小了计算的误差,析计算渗流情况下的堤坝边坡稳定性分析。首先对堤坝进

使得解的精度提升了一个层次。同时给出了一种新的弃单行建模,采用饱和非饱和理论进行计算，在浸润面之上选用

元网格法,通过逐渐迭代来推算出真实的浸润线以及渗流土的非饱和抗剪强度，浸润面之下用饱和土的抗剪强度。

区域的大小。陈洪凯⑸通过对三峡工程边坡岩体发生渗流在除了对浸润面之下采用非饱和的抗剪强度，同时将渗流

过程的研究，提出了全新的整合单元全区域逼近法。张培力平均分配到浸润面之下的各单元上，堤坝渗流得到的浸

文⑹通过有限元模拟了边坡孔隙水压力的变化，改进了边润面如图1所示。

坡边界条件转化的问题，可以定量分析降雨导致的边坡失

稳。王均星等日］提出了流动形单元法，利用区域覆盖的方

式区分浸润面，用于有自由面的渗流更加简便有效。黄燕

宏〔5］利用热传导控制方程和渗流微分方程的相似性，通过

APDL语言设计了程序进行渗流计算。

SzaPoan/Mccraing/：对瞬态非稳定流利用有限差分法

进行了数值模拟。R1A.freeze：12：对地下含水层的饱和一非图1堤坝渗流浸润面图

饱和渗流进行了分析，在原有问题的基础上,通过数值模拟模型设置的主要参数如表1所示。

的方法进行了改进。之后随着计算机的高速发展，数值模表1模型参数设置

参数名称大小含义

拟的方法开始在工程建设中得到大量应用。目前，国内对