论《dqb》程序和《stab》程序在土石坝安全鉴

《论《dqb》程序和《stab》程序在土石坝安全鉴》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、论《DQB》程序和《STAB》程序在土石坝安全鉴摘要：土石坝二向稳定及非稳定渗流计算程序《DQB》，系由南京水利科学研究院水工所李祖贻、陈平等同志编制，用FORTRAN语言在TQ-16机及IBM-PC/XT机实现。该程序既可用于稳定渗流分析，又可用于非稳定渗流分析，并能适用于均质、心墙、斜墙土坝不同排水型式的变化。程序采用自动剖分单元，数据准备工作量小，算题速度快，是土石坝分析的有效工具之一，经1995年水电总局考核通过，列为在水电系统推广应用的土石坝计算程序包十个程序之一。关键词：DQB程序STAB程序土石坝安全鉴定1前言土石坝二

2、向稳定及非稳定渗流计算程序《DQB》，系由南京水利科学研究院水工所李祖贻、陈平等同志编制，用FORTRAN语言在TQ-16机及IBM-PC/XT机实现。该程序既可用于稳定渗流分析，又可用于非稳定渗流分析，并能适用于均质、心墙、斜墙土坝不同排水型式的变化。程序采用自动剖分单元，数据准备工作量小，算题速度快，是土石坝分析的有效工具之一，经1995年水电总局考核通过，列为在水电系统推广应用的土石坝计算程序包十个程序之一。土石坝边坡稳定分析程序《STAB》是根据水利水电科学研究院陈祖煜同志所编M-16机土石坝边坡稳定分析程序《STAB》中的

3、简化法。编者应用以上二个程序于数座土石坝安全鉴定，均取得了满意效果。2程序的使用范围及功能《DQB》程序可用来计算土石坝上游坝壳水位降落期的非稳定渗流和具有不同排水型式的均质、心墙、斜墙土坝的稳定渗流，以及任意过流断面的渗流量。该程序具有自动部分功能，只要给出剖分信息、单元及结点信息即可由程序自动形成，并计算给出自由表面线（浸润线）位置，全部结点水头值、不同百分数的等势线等计算成果。《STAB》程序可以同时用瑞典法、毕肖普法和改良瑞典圆弧法（罗厄法或工程师兵团法）算出圆弧滑裂面的安全系数，并找出相应于毕肖普法的最小安全系数及相应的滑

4、弧位置。该程序能用总庆力法或有效应力法计算施工期、稳定渗流期、库水位骤降期以及以上各期遇有地震时的边坡稳定分析，强度有单强度及组合强度，使用者对程序功能的通过功能控制变量的赋值来实现。3计算原理3.1《DQB》程序计算原理3.1.1数学方法与边界条件土坝二向渗流问题在一定条件下，稳定渗流是求解拉普拉斯议程，非稳定渗流在土体可压缩时求固结议程，在土体不可压缩时求解拉普拉斯议程，同时自由面作渗流量补给边界。对上述议程用有限单元作渗流场离散，引用三结点的三角形单元和线性插值函数，线性代数方程组改进平方根法求解。离散后所得的线性代数议程组为

5、：（1）（1）式为可压缩土体的非稳定渗流有限单元计算公式。当式中矩阵[S]=0时，得不可压缩土体的非稳定渗流公式：（2）在不计时间项，且[S]、[P]矩阵等于零，得稳定渗流有限单元计算公式：（3）在坝基无潜流时，其中[D]、项为零。土坝渗流边界条件是：上、下游水位以下的入流和出流面及自由渗出段，其水头是已知的，属第一类边界；渗流自由面和不透水层面属第二类边界。稳定渗流的自由面和不透水层面一样，没有流量从该面流入和流出；而非稳定渗流则有流量从自由面流进坝体，此流量是取二连续自由面之间的一块水体表示补给流量。此外自由面上尚应满足其水头等

6、于位置高程的条件。3.1.2渗流量的计算渗流量是采用中线法计算的，可计算任意过流断面的流量，渗流量计算公式为：（4）上式是具有方向性的，计算流量断面是取划分单元的初始坐标的方向，故计算时要规定其正向，然后相加得总渗流量。3.1.3渗流自由面的确定首先假定一渗流自由面位置。原则上依照渗流概念按“简化法”计算浸润线，同时对于不同排水设施的渗出段长度a，以及排水起点处的渗流水深h0作出粗估，尽量接近实际浸润线位置，以减少浸润线的修改迭代次数，节省机时，同时也可避免出错。然后计算机程序将自由面结点计算水头值h*与其Z坐标比较，直到满足

7、A-

8、h*

9、＜ε（ε为给定的计算精度），此时则获正确的理论自由面位置。自由选代时，沿结点线上下调正移动，为此要求结点线上结点按顺序从上到下由小到大编号。为避免假定的自由面位置过高，或自由面穿过非均质区及非稳定渗流自由面变化范围时，使计算不能持续进行下去，程序采用丢结点（单元）的方法来处理，即自由面结点的计算水头h*小于其下结点Z坐标时，将该结点丢弃，取其下结点为自由面结点继续进行计算。3.1.4渗出点的确定自由面与坝坡的交点即渗出点，采用沿坡面滑动或二次曲线相交法求得，根据不同情况选用。渗出点一般是作为未知水头结点处理，也可按已知水头结点

10、处理。调正渗出点的过程中可以由程序自动改变信息，反复试求得出渗出点的正确位置。程序的关键在于合理地确定渗流场的计算范围、边界性质及单元信息。3.2《STAB》程序计算原理3.2.1程序采用的强度指标（1）不固结不排水（三轴仪）或不固结