求解堤坝渗流场的罚函数无单元法 - 河海期刊.pdfVIP

第 卷第 期 河 海 大 学 学 报 （自然 科 学 版 ） )+ $ 15 . )+ E . $ 年 月 （ ） !- $ ;=7:45 ?49 @:9A67B9CD E4C=745 (F96:F6B ;4: .!- 求解堤坝渗流场的罚函数无单元法 ， $ $ ! ) 沈振中，陈小虎 ，吴越建 （ 河海大学水利水电工程学院，江苏 南京 ； 长江勘测规划设计研究院枢纽设计处，湖北 武汉 ； $. !$%- ! . /)$ 中国水利水电第 工程局 分局，四川 成都 ） ) . # ) +$$$) 摘要：从变分原理出发，采用罚函数法处理渗流边界条件，推导了用无单元法求解有自由面的堤坝 渗流问题的基本方程和迭代格式，给出了罚因子选取的具体表达式，并计算分析了某不透水地基上 土坝的稳定渗流场. 与电模拟试验成果的比较表明，文中所提出的计算方法和程序是正确的，且应 用罚函数无单元法求解堤坝渗流问题，可以得到较高的计算精度和计算效率. 关键词：堤坝；渗流；边界条件；罚函数法；无单元法；变分原理 中图分类号： 文献标识码： 文章编号： （ ） 01!!) ./ 2 $ $%- !- $ // , ! ! ! 有自由面的渗流问题普遍存在于堤坝、岩土边坡、地下洞室等实际工程中，目前常采用有限单元法求解. 由于单元网格的限制，复杂渗流问题的有限元前处理工作量十分庞大 无单元法是一种新兴的工程数值计算 . 方法，它采用滑动最小二乘法近似场函数，摆脱了单元的限制，计算过程中可以任意布置节点，因此在求解渗 流自由面等可动边界问题时优势明显. 迄今为止，无单元法的研究主要集中在裂纹扩展追踪、变形和应力分 ［ ］ ［ ］ $ , + % 析等领域，且研究成果较多 ! ，并已发展了三维无单元法 ! . 但是，在渗流问题求解方面的研究尚不多见. ［ ］ ［ ］ $ $$ 曾清红等 将无单元法用于求解稳定渗流问题，李广信等 提出了有自由面渗流的无单元法，但他们都没 ［ ］ 有考虑如何准确施加渗流边界条件的问题 $! . 李晓春 将无