非稳定渗流作用下土石坝边坡稳定性研究.pdfVIP

10 内蒙 古水 利 2012年第6期(总第 142期) 【分析研究 】 非稳定渗流作用下土石坝边坡稳定性研究 刘春林 ，刘淑英 ，韩 凌。 (1．内蒙古自治区水利水电勘测设计院，内蒙古 呼和浩特010020；2．内蒙古自治区水土保持工作站，内蒙古 呼和浩特o10020； 3．内蒙古河套灌区义长灌域管理局，内蒙古 五原015100) [摘 要] 文章以考虑压缩性的非稳定性渗流模型为基础，引入有限元法与渗流力作用原理，实现 对土石坝边坡稳定性的数值模拟，并引入工程实例验证 了该法的合理性。 [关键词] 非稳定渗流；有限元；边坡稳定性 中图分类号 ：TV223 文章标识码：B 文章编号：1009—0088(2012)06—0010—02 1 非稳定渗流场数学模型 IC『iCi ciej CiCm] c c． c c c 。cm l l 。 。 c c,n％ J (k )+ (k )=s () ， 水头边界h jr=h(x、z、t) ： ： 三02：] 流量边界k，_Ohlr2=一q(x L0 1 2J 、 z、t) 式中 △一三角形单元的面积 ； 自由面边界 h=z，q= c。s。 一 给水度 ； k 、k一 初始边界h I =h(x，z，0) z 渗透系数； 式中 k、k、k一渗透系数； 其余参数可由三角形单元三节点的坐标值求 S一单位贮水量 。 得 。 2 非稳定渗流场有限元求解 3 土石坝边坡稳定计算 非稳定渗流场数学模型的总体有限元方程为： 首先应用有限元法求解单元渗透力，然后根据极 ([K] )lht +1 P『] 限平衡理论建立起土石坝边坡稳定的分析表达式。直 接应用渗流场有限元法分析中所划分单元和计算结点 {hft 一1 … t=… (2) 水头值，持续计算单元渗透力，统一编制程序，进行土 坡稳定性分析计算，既提高了计算精度，又保持了计算 总系数矩阵[K]、[S]、[P]和常数列向量 [F]中 的连贯性 。 的典型元素都是对各单元求和，即Ki善kij，Sij 一 R{Zcl+∑[(r1+rJ)COSOtrJsine~]△tg中 ，。、 ∑[(r1+rJ)sino~+rJ(2rcoso~一Rcoso~)]△r k ，P kpiJ，F k． 式中 e一滑动土体有效粘聚力； 式中 m一单元数；