›› 2012, Vol. 33 ›› Issue (S1): 237-243.

• 数值分析 • 上一篇    下一篇

非饱和渗流Richards方程数值求解的欠松弛方法

陈 曦1,于玉贞2,程勇刚3   

  1. 1. 北京交通大学 土木建筑工程学院,北京 100044;2. 清华大学 水沙科学与水利水电工程国家重点实验室,北京 100084; 3. 武汉大学 水资源与水电工程科学国家重点实验室,武汉 430072
  • 收稿日期:2002-05-28 出版日期:2012-09-11 发布日期:2012-09-21
  • 作者简介:陈曦,男,1977年生,博士,副教授,主要从事岩土工程数值计算方法等方面的研究
  • 基金资助:

    国家自然科学基金资助项目(No. 50879039);973计划课题资助(No. 2010CB732103);中央高校基本科研业务费专项资金资助(No. 2011JBM070)

Under-relaxation methods for numerical solution of Richards equation of variably saturated flow

CHEN Xi1, YU Yu-zhen2, CHENG Yong-gang3   

  1. 1. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China; 2. State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China; 3. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
  • Received:2002-05-28 Online:2012-09-11 Published:2012-09-21

摘要: 非饱和土渗流理论是岩土工程问题的基础理论,在土石坝渗流、污染物传输、冻土渗流相变和边坡稳定分析等领域有着广泛的应用。非饱和土渗流Richards方程的数值求解过程中,某些参数如水力传导系数计算不当可能引起非线性方法,如Picard方法或Newton方法的迭代收敛震荡,从而导致非线性迭代方法收敛缓慢和精度降低。为了消除或降低迭代收敛震荡对求解精度和计算性能的影响,目前主要采用欠松弛方法。通过一维入渗算例和二维非均质土坝渗流算例演示已有欠松弛方法的局限性,进而提出新的短项混合欠松弛法,并对其实用性和可靠性进行验证。

关键词: 非饱和渗流, Richards方程, 有限元, 迭代收敛震荡, 欠松弛法

Abstract: Variably saturated flow theory is the fundamental theory of geotechnical engineering; and it has wide applications to seepage of earth rockfill dam, contaminant transport, seepage and phase transformation of frozen soil, and slope stability analysis. In numerical solution of Richards equation for variably saturated flow, inappropriate computation of some parameters, particularly for the hydraulic conductivity, may lead to convergence oscillation of nonlinear iterative methods such as Picard method or Newton method, and consequently result in slow convergence and inaccurate solution. To eliminate or reduce the influence of iterative convergence oscillation, the under-relaxation method is often employed. By using 1D problem and 2D heterogeneous earth dam, the limitations of available under-relaxation methods are demonstrated. Furthermore, a novel short-term mixed under-relaxation method is proposed with verified practicability and reliability.

Key words: variably saturated flow, Richards equation, finite elements, iterative convergence oscillation, under-relaxation

中图分类号: 

  • TU 443
[1] 孙锐, 杨峰, 阳军生, 赵乙丁, 郑响凑, 罗静静, 姚捷, . 基于二阶锥规划与高阶单元的 自适应上限有限元研究[J]. 岩土力学, 2020, 41(2): 687-694.
[2] 叶观宝, 郑文强, 张 振, . 大面积填土场地中摩擦型桩负摩阻力分布特性研究[J]. 岩土力学, 2019, 40(S1): 440-448.
[3] 赵密, 欧阳文龙, 黄景琦, 杜修力, 赵旭, . P波作用下跨断层隧道轴线地震响应分析[J]. 岩土力学, 2019, 40(9): 3645-3655.
[4] 张海廷, 杨林青, 郭芳, . 基于SBFEM的层状地基埋置管道动力 响应求解与分析[J]. 岩土力学, 2019, 40(7): 2713-2722.
[5] 曹洪, 胡瑶, 骆冠勇. 滤管两端均不在含水层层面的承压不 完整井近似计算方法研究[J]. 岩土力学, 2019, 40(7): 2774-2780.
[6] 王翔南, 李全明, 于玉贞, 喻葭临, 吕禾, . 基于扩展有限元法对土体滑坡破坏过程的模拟[J]. 岩土力学, 2019, 40(6): 2435-2442.
[7] 周小文, 程 力, 周 密, 王 齐, . 离心机中球形贯入仪贯入黏土特性[J]. 岩土力学, 2019, 40(5): 1713-1720.
[8] 梅慧浩, 冷伍明, 聂如松, 刘文劼, 伍晓伟, . 重载铁路路基面动应力峰值随机分布特征研究[J]. 岩土力学, 2019, 40(4): 1603-1613.
[9] 吴顺川, 马 骏, 程 业, 成子桥, 李建宇, . 平台巴西圆盘研究综述及三维启裂点研究[J]. 岩土力学, 2019, 40(4): 1239-1247.
[10] 谢 强, 田大浪, 刘金辉, 张建华, 张志斌, . 土质边坡的饱和−非饱和渗流分析及特殊应力修正[J]. 岩土力学, 2019, 40(3): 879-892.
[11] 邱 敏, 袁 青, 李长俊, 肖超超, . 基于孔穴扩张理论的黏土不排水抗剪强度 计算方法对比研究[J]. 岩土力学, 2019, 40(3): 1059-1066.
[12] 郑安兴, 罗先启, 陈振华, . 基于扩展有限元法的岩体水力劈裂耦合模型[J]. 岩土力学, 2019, 40(2): 799-808.
[13] 王冬勇, 陈曦, 于玉贞, 吕彦楠, . 基于二阶锥规划有限元增量加载法的条形浅基础极限承载力分析[J]. 岩土力学, 2019, 40(12): 4890-4896.
[14] 刘庆彬, 潘懋, 刘洁, 郭艳军, 张小双, 姚健鹏, 李芳玉, . 基于ParaView的Abaqus有限元输出结果的 可视化与虚拟现实[J]. 岩土力学, 2019, 40(12): 4916-4924.
[15] 刘锋涛, 张绍发, 戴北冰, 张澄博, 林凯荣, . 边坡稳定分析刚体有限元上限法的锥规划模型[J]. 岩土力学, 2019, 40(10): 4084-4091.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 谭贤君,陈卫忠,杨建平,杨春和. 盐岩储气库温度-渗流-应力-损伤耦合模型研究[J]. , 2009, 30(12): 3633 -3641 .
[2] 魏 星,王 刚,余志灵. 交通荷载下软土地基长期沉降的有限元法[J]. , 2010, 31(6): 2011 -2015 .
[3] 温世亿,李静,苏霞,姚雄. 复杂应力条件下围岩破坏的细观特征研究[J]. , 2010, 31(8): 2399 -2406 .
[4] 毛 宁,张尧亮. 经验公式简便求法典型实例[J]. , 2010, 31(9): 2978 -2982 .
[5] 刘 杰,李建林,屈建军,陈 星,李剑武,骆世威. 基于卸荷岩体力学的大岗山坝肩边坡水平位移发育的多因素影响分析[J]. , 2010, 31(11): 3619 -3626 .
[6] 李伟华,赵成刚,杜楠馨. 软弱饱和土夹层对地铁车站地震响应的影响分析[J]. , 2010, 31(12): 3958 -3963 .
[7] 韩现民. 西格二线关角隧道浅埋砂层段施工技术及力学效应研究[J]. , 2010, 31(S2): 297 -302 .
[8] 刘用海,朱向荣,常林越. 基于Casagrande法数学分析确定先期固结压力[J]. , 2009, 30(1): 211 -214 .
[9] 祝 磊,洪宝宁. 粉状煤系土的物理力学特性[J]. , 2009, 30(5): 1317 -1322 .
[10] 周春梅,章泽军,徐大杰,王生维,李先福. 古构造应力场数值模拟及危险性预测研究[J]. , 2009, 30(7): 2141 -2146 .