论《DQB》程序和《STAB》程序在土石坝安全鉴定中的应用

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

2、院水工所李祖贻、陈平等同志编制，用FRTRAN语言在TQ-16机及IB-P/XT机实现。该程序既可用于稳定渗流分析，又可用于非稳定渗流分析，并能适用于均质、心墙、斜墙土坝不同排水型式的变化。程序采用自动剖分单元，数据准备工作量小，算题速度快，是土石坝分析的有效工具之一，经1995年水电总局考核通过，列为在水电系统推广应用的土石坝计算程序包十个程序之一。土石坝边坡稳定分析程序?STAB?是根据水利水电科学研究院陈祖煜同志所编-16机土石坝边坡稳定分析程序?STAB?中的简化法。编者应用以上二个程序于数座土石坝平安鉴定，均获得了满意效果。2程序的使用范围及功能?DQB?程序可用来计算土石坝上游坝壳

3、水位降落期的非稳定渗流和具有不同排水型式的均质、心墙、斜墙土坝的稳定渗流，以及任意过流断面的渗流量。该程序具有自动局部功能，只要给出剖分信息、单元及结点信息即可由程序自动形成，并计算给出自由外表线浸润线位置，全部结点水头值、不同百分数的等势线等计算成果。?STAB?程序可以同时用瑞典法、毕肖普法和改进瑞典圆弧法罗厄法或工程师兵团法算出圆弧滑裂面的平安系数，并找出相应于毕肖普法的最小平安系数及相应的滑弧位置。该程序能用总庆力法或有效应力法计算施工期、稳定渗流期、库水位骤降期以及以上各期遇有地震时的边坡稳定分析，强度有单强度及组合强度，使用者对程序功能的通过功能控制变量的赋值来实现。3计算原理3.

4、1?DQB?程序计算原理3.1.1数学方法与边界条件土坝二向渗流问题在一定条件下，稳定渗流是求解拉普拉斯议程，非稳定渗流在土体可压缩时求固结议程，在土体不可压缩时求解拉普拉斯议程，同时自由面作渗流量补给边界。对上述议程用有限单元作渗流场离散，引用三结点的三角形单元和线性插值函数，线性代数方程组改进平方根法求解。离散后所得的线性代数议程组为：11式为可压缩土体的非稳定渗流有限单元计算公式。当式中矩阵S=0时，得不可压缩土体的非稳定渗流公式：2在不计时间项，且S、P矩阵等于零，得稳定渗流有限单元计算公式：3在坝基无潜流时，其中D、项为零。土坝渗流边界条件是：上、下游水位以下的入流和出流面及自由渗出

5、段，其水头是的，属第一类边界；渗流自由面和不透水层面属第二类边界。稳定渗流的自由面和不透水层面一样，没有流量从该面流入和流出；而非稳定渗流那么有流量从自由面流进坝体，此流量是取二连续自由面之间的一块水体表示补给流量。此外自由面上尚应满足其水头等于位置高程的条件。3.1.2渗流量的计算渗流量是采用中线法计算的，可计算任意过流断面的流量，渗流量计算公式为：4上式是具有方向性的，计算流量断面是取划分单元的初始坐标的方向，故计算时要规定其正向，然后相加得总渗流量。3.1.3渗流自由面确实定首先假定一渗流自由面位置。原那么上按照渗流概念按“简化法计算浸润线，同时对于不同排水设施的渗出段长度a，以及排水起

6、点处的渗流水深h0作出粗估，尽量接近实际浸润线位置，以减少浸润线的修改迭代次数，节省机时，同时也可防止出错。然后计算机程序将自由面结点计算水头值h*与其Z坐标比拟，直到满足|A-h*|为给定的计算精度，此时那么获正确的理论自由面位置。自由选代时，沿结点线上下调正挪动，为此要求结点线上结点按顺序从上到下由小到大编号。为防止假定的自由面位置过高，或自由面穿过非均质区及非稳定渗流自由面变化范围时，使计算不能持续进展下去，程序采用丢结点单元的方法来处理，即自由面结点的计算水头h*小于其下结点Z坐标时，将该结点丢弃，取其下结点为自由面结点继续进展计算。3.1.4渗出点确实定自由面与坝坡的交点即渗出点，采

7、用沿坡面滑动或二次曲线相交法求得，根据不同情况选用。渗出点一般是作为未知水头结点处理，也可按水头结点处理。调正渗出点的过程中可以由程序自动改变信息，反复试求得出渗出点的正确位置。程序的关键在于合理地确定渗流场的计算范围、边界性质及单元信息。3.2?STAB?程序计算原理3.2.1程序采用的强度指标1不固结不排水三轴仪或不固结快剪直剪仪简称Q剪，相应指标uu，uu；2固结排水三轴仪或固结慢剪直剪仪简称S剪，相应指标d、d；3固结不排水三轴仪或固结快剪直剪仪简称R剪，相应指标u、u；4不固结不排水或固结不排水，测孔压三轴仪，简称；5常用的现场强度试验所获得的强度指标为=0的值，用qu代表，成果与R

8、剪指标等效，qu值为现场被试验地点坐标x，y的函数，程序要求将qux，y按坐标网格输入，然后自动进展内插，确定滑弧面上各点的qu值即值。3.2.2有效应力法有效应力法的强度表达式为：=+-utan5强度指标和是根据、S、剪切试验确定的，、S、试验结果分别用于施工期、稳定渗流期、库水位降落期。用有效应力法计算滑弧稳定时，孔隙水压力u值事先已求出，可按网格逐点输入，程序自动内插，确定滑弧面上各点的孔压值。假如孔隙水压在某一区域内按静压分布，或浸润线较平缓，可近似假定等势线垂直，即孔压按静压分布，那么只需输入浸润线位置即可，这种情况一般用稳定渗流期。3.2.3总应力法总应力法是将土样在试验室模拟坝体

9、实际工况条件下测定其强度指标，把超孔隙水压包涵在抗剪强度指标内的一种方法。本程序的决庆力法能计算施工期及库水位骤降期的坝坡稳定分析，施工期除用Q剪指标外，计算时与有效应力法一样。库水位骤降时的总应力法有以下两种：1采用现场快剪强度指标qu；2采用RS组合强度包线。3.2.4程序中地震惯性力的计算按?水工建筑物抗震设计标准?SDJ10-78第16条和第23条规定计算。3.2.5程序图形处理约定土石坝断面被看成是由一系列边界限构成，土层不同界限及外边坡均看作边界限。程序规定任一铅直线均不得同时碰到两根浸润线，但在滑弧可能到达的范围内必须碰到一根浸润线；浸润线边界限都不允许是垂直线。当所计算的滑动边坡外有水时，程序要求将水位线以下的外边坡线作为浸润线。3.2.6滑动圆弧控制参数及布置方式圆弧可计算一个或多个。计算多个圆弧时有两种布置方式：一是圆弧始终通过上、下交点，半径那么按一定步长变化；二是以圆心及滑弧深度为中心，给予它们变化的步长后，程序将自动算出所有组合情况下的平安系数。假如某圆弧不与坝坡相交，程序将自动识别，并打印出信息。4应用在1998年邵阳市病险水库大坝平安鉴定工作中，编者成功地将?DQB?程序和?STAB?程序应用到邵东流光岭水库和大圳灌区东风水库两座中型大坝的平安鉴定中。由于以上二程序数据准备工作量小，计算速度快，节约了大量的时间