论《DQB》程序和《STAB》程序在土石坝安全鉴定中的应用 来源：邵阳市水利水电勘测设计院 作者：王飚 日期：[2006-2-14] 土石坝二向稳定及非稳定渗流计算程序《DQB》，系由南京水利科学研究院水工所李祖贻，陈平等到同志编制，用FORTRAN语言在TQ-16机及IBM-PC/XT机实现。该程序既可用于稳定渗流分析，又可用于非稳定渗流分析，并能适用于均质、心墙、斜墙土坝不同排水型式的变化。程序采用自动剖分单元，数据准备工作量小，算题速度快，是土石坝分析的有效工具之一，经1995年水电总局考核通过，列为在水电系统推广应用的土石坝计算程序包十个程序之一。  土石坝边坡稳定分析程序《STAB》是根据水利水电科学研究院陈祖煜同志所编M-16机土石坝边坡稳定分析程序《STAB》中的简化法。  编者应用以上二个程序于数座土石安全鉴定，均取得了满意效果。  1. 程序的使用范围及功能 《DQB》程序可用来计算土石坝上游坝壳水位降落期的非稳定渗流和具有不同排水型式的均质、心墙、斜墙土坝的稳定渗流，以及任意过流断面的渗流量。该程序具有自动剖分功能，只要给出剖分信息，单元及结点住处即可由程序自动形成，并计算给出自由表面线（浸润线）位置，全部结点水头值，不同百分数的等势线等计算成果。  《STAB》程序可以同时用瑞典法，毕肖普法和改良瑞典圆弧法（罗厄法或工和师兵团法）算出贺圆弧滑面的安全系数，并找出相应于毕肖普法的最小安全系数及相应的滑弧位置。该程序能用总应力法或有效应力法计算施工期，稳定渗流其，库水位骤降期以及以上各期遇有地震时的边坡稳定分析，强度有音强度及组合强度，使用者对程序功能的选用可通过功能控制变量的赋值来实现。  2. 计算原理  2．1《DQB》程序计算原理  2．1．1数学方法与边界条件  土坝二向渗流总是在一定条件下，稳定渗流是求解拉普拉斯议程，非稳定渗流在土体可压缩时求解拉普拉斯方程，同时自由面作渗流量补给边界。对上述方程用有限单元作渗流场离散，引用三结点的三角形单元和线性插值函数，线性代数方程组改进平方根肢求解。离散后所得的线性代数方程组为：  （[K]+1/△t[S]）{h}t+△t+1/△t  {h*}t+△t-1/△t[S]{h}t-1/△t  {h*}t={F}-[D]{q}(1)  (1)式为可压缩土体的非稳定渗流有限计算公式。当式中矩阵[S]=0时，得不可压缩土体的非稳定渗流公式： [K]{h}t+△t+1/△t  {h*}t+△t-1/△t  {h*}t={F}-[D]{q}(2)  在不计时间项，且[S]，  矩阵等于零，得稳定渗流有限