求解渗流自由面的复合单元全域迭代法ξ.pdfVIP

( ) 　应用数学和力学 ,第 20 卷 第 10 期 1999 年 10 月 应用数学和力学编委会编 　　　　　　　　　　 　 　　Applied Mathematics and Mechanics 重 庆 出 版 社 出 版 文章编号 :1000-0887 ( 1999) 10-1045-06 求解渗流自由面的复合单元全域迭代法 陈洪凯 , 　唐红梅 , 　肖盛燮 (重庆交通学院 河海建筑工程系 ,重庆 400074) (戴世强推荐) 摘要 : 　渗流自由面是岩体水力学研究的重点和难点 , 也是边坡稳定性综合研究的重要内容 ·　基 于对渗流自由面数值求解方法的分析 , 提出了复合单元全域迭代法 , 并编制了有限元计算程序 ·　 计算果表明 , 该法具有程序处理比较简单、计算结果比较精确等特性 , 一般迭代 5 次左右便可收 敛 ·　 关 　键 　词 : 　岩体边坡 ; 　渗流自由面 ; 　复合单元全域迭代法 ; 　验证 中图分类号: 　P6422 　　　文献标识码 : 　A 引　　言 在边坡、土坝、地下洞室及地下水运动等渗流分析中 ,均涉及渗流自由面问题 ·　渗流自由 面是岩体水力学研究的重点和难点之一[1～3 ] ,也是边坡稳定性综合研究的重要内容 ·　岩体边 坡中渗流自由面的确定 ,可以明确地下水在边坡中的赋存特性及其运动过程 ·　有限元法、边界 元法及离散元法等是求解无压渗流自由面主要的数值计算方法 ·　 在求解渗流自由面的数值计算方法中 ,使用效果最好的方法当属有限元法[4 ] ,传统的有限 元方法属于网格变动法 ,即在每次迭代以后计算网格作一次变动 ,节点坐标发生改变 ·　由于渗 流自由面作为渗流域的自然边界面是待定的 ,是一个非线性问题 , 因而需要迭代求解 ·　但是, 网格变动法存在着比较显著的缺陷: 1) 当初始渗流自由面和最终渗流自由面位置相差较大时 ,会使计算单元发生畸变 ,乃至 与相邻单元发生交替、重叠 , 以至于在计算过程中常需要对渗流域进行重新剖分计算; 2) 当自由面附近渗流介质不均一 ,尤其有水平介质层时网格变动会破坏介质分区 ,程序 处理十分困难 ; 3) 当渗流域内有结构物时 ,网格变动常会改变结构的边界条件 ,影响计算精度 ; 4) 网格变动过程中 ,每一次迭代计算网格均要随自由面的变动而变动 ,总体渗透矩阵需 要重新生成 ,需要大量机时; 5) 在研究渗流与应力耦合作用中 , 由于应力分析经常包括渗流虚区, 因而不能用同一网 格进行耦合作用分析 ·　 自从1973 年Neumann 提出用网格不变法进行渗流分析以来 ,Desai[5 ] 提出了残余流量法、 收稿日期 : 　1998-11-23 ; 修订日期 : 　1999-01-08 ( ) ( ) 基金项目 : 　国家自然科学基金资助项目 ; 国家“八五”重点科技攻关项目 ( ) 作者简介 : 　陈洪凯 1964～ ,副教授 ,博士. 1045 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 1046 陈 　洪 　凯 　　　唐 　红 　梅 　　　肖 　盛 　燮 张有天提出了初流量法[6 ] 、速宝玉提出了节