非饱和渗流问题的自适应欠松弛变量变换方法

›› 2012, Vol. 33 ›› Issue (9): 2857-2862.

非饱和渗流问题的自适应欠松弛变量变换方法

程勇刚¹，常晓林¹，李典庆¹，陈曦²

1.武汉大学水资源与水电工程科学国家重点实验室，武汉 430072；2.北京交通大学土木建筑工程学院，北京 100044

收稿日期:2012-05-10 出版日期:2012-09-11 发布日期:2012-09-12
作者简介:程勇刚，男，1977年生，博士，讲师，主要从事水工结构、岩土工程数值方法等方面的研究工作
基金资助:
国家自然科学基金项目（No. 51028901）；高等学校博士学科点专项科研基金资助课题（No. 20110141120019）；中国博士后科学基金资助项目（No. 20110491214）。

Unsaturated seepage analysis using temporal adaptive transformation method with under-relaxation

CHENG Yong-gang¹, CHANG Xiao-lin¹, LI Dian-qing¹, CHEN Xi²

1. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China; 2. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China

Received:2012-05-10 Online:2012-09-11 Published:2012-09-12

摘要/Abstract

摘要： 在使用有限元方法求解非饱和土渗流问题时，土-水特征曲线和渗透率函数的强烈非线性经常会造成计算中出现迭代不收敛、计算误差大等问题。基于变量变换的思想，结合时间步长自适应技术提出了一种求解非饱和渗流问题的新方法——欠松弛RFT变换方法（ATUR1）。ATUR1方法通过变量变换，大大降低了Richards方程中未知数在空间和时间上的非线性程度，从而改善这种非线性所带来的计算收敛困难和精度差等问题。欠松弛技术的引入减少了迭代过程中的振荡现象，进一步提高了非线性迭代计算的效率。时间步长自适应技术则有效地控制整个计算过程的误差。数值算例结果说明，ATUR1可以有效地提高计算效率和精度，是一种准确有效的计算方法。

关键词: 非饱和渗流, 有限元, 变量变换, 欠松弛, 时间步长自适应

Abstract: Finite element method is often used for solving unsaturated soil seepage problem. However, numerical difficulties such as slow convergence rate and inaccurate results are commonly existed in numerical analysis because of the highly nonlinear soil-water characteristic curve and hydraulic conductivity function. Based on a rational function transformation (RFT) method, a new temporal adaptive transformation method with under-relaxation (ATUR1) is proposed to solve the unsaturated soil seepage problem. It is shown that the transformation in ATUR1 can reduce the nonlinearity of solution profiles of the Richards' equation both spatially and temporally; thus the inefficiencies in the numerical solution process caused by these nonlinearities can be overcome. The under-relaxation technique can further increase the efficiency by curbing the oscillation tendency during iterations; and the temporal errors during the whole simulation can be controlled effectively by the automatic temporal adaptive scheme. Numerical examples show that the proposed ATUR1 method is a robust and efficient numerical method for practical unsaturated soil seepage simulation.

Key words: unsaturated seepage, finite elements, variable transformation, under-relaxation, temporal adaptive

中图分类号:

P 641.2

程勇刚，常晓林，李典庆，陈曦 . 非饱和渗流问题的自适应欠松弛变量变换方法[J]. , 2012, 33(9): 2857-2862.

CHENG Yong-gang , CHANG Xiao-lin , LI Dian-qing , CHEN Xi . Unsaturated seepage analysis using temporal adaptive transformation method with under-relaxation[J]. , 2012, 33(9): 2857-2862.

参考文献

相关文章 15

[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.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed