›› 2012, Vol. 33 ›› Issue (1): 295-300.

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

近饱和条件下非饱和多孔介质渗流过程的数值分析

陈 盼,韦昌富,王吉利,伊盼盼,曹华峰   

  1. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,武汉 430071
  • 收稿日期:2010-05-26 出版日期:2012-01-10 发布日期:2012-01-17
  • 通讯作者: 韦昌富,男,1966年生,博士,研究员,博士生导师,主要从事非饱和土力学与多相孔隙介质多场耦合问题方面的研究工作。 E-mail: cfwei@whrsm.ac.cn E-mail:pchen@whrsm.ac.cn
  • 作者简介:陈盼,男,1984年生,博士,主要从事非饱和土渗流与边坡稳定性方面的研究工作。
  • 基金资助:

    国家自然科学基金资助项目(No. 10872211,No. 11072255),岩土力学与工程国家重点实验室重点方向性项目资助(No. kzcxz-yw-150)。

Numerical analysis of seepage processes in unsaturated porous media under nearly saturated conditions

CHEN Pan,WEI Chang-fu,WANG Ji-li,YI Pan-pan,CAO Hua-feng   

  1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
  • Received:2010-05-26 Online:2012-01-10 Published:2012-01-17

摘要: 在非饱和多孔介质渗流分析中,近饱和条件下物理模型与数值模型之间的差异会导致数值不稳定问题。为解决这一问题,并保证模拟结果的可靠性,提出了3种方法,并在有限元分析程序U-DYSAC2中分别进行了程序代码的实施。通过数值试验与试验数据的比较,证实了在近饱和条件下土-水特征曲线和水力传导函数的高非线性可引起数值收敛性、稳定性和精度问题,而且在不同条件下含水率和基质吸力的预测结果差异明显。在3种方法中,修正的Van Genuchten模型(MVGM)方法对含水率的预测较为准确,而Line方法对基质吸力的预测较为合理。因此,解决在分析近饱和条件下非饱和多孔介质渗流问题时,为获得接近真实的模拟结果,采用合适的数值方法进行预测是非常关键的。

关键词: 近饱和状态, 土-水特征, 渗流分析, 非饱和多孔介质, 数值稳定性

Abstract: In analyzing the seepage processes in unsaturated porous media, numerical instability may occur due to the difference between the adopted physical and numerical models under nearly saturated conditions. To resolve the issue and improve the numerical results, three methods are proposed and then implemented into the finite element code U-DYSAC2. Numerical simulations are compared with the experimental data available in the literature, showing that the problems related to numerical stability, convergence and accuracy can be induced by the high nonlinearity of soil water characteristics and hydraulic conductivity functions under nearly saturated conditions; and significant numerical discrepancies can appear in predicting the water content and matric suction of the porous media. Among the three proposed methods, modified Van Genuchten model (MVGM) method yields the best results in calculating the moisture content; while the Line method achieves best prediction for matric suction. Consequently, in analyzing the seepage processes in unsaturated porous media, it is crucial to introduce a proper analysis procedure such that reliable numerical results can be achieved

Key words: near saturated state, soil moisture characteristics, seepage flow analysis, unsaturated porous media, numerical stability

中图分类号: 

  • TU 411
[1] 陈仁朋, 王朋飞, 刘鹏, 程威, 康馨, 杨微, . 路基煤矸石填料土-水特征曲线试验研究[J]. 岩土力学, 2020, 41(2): 372-378.
[2] 陶高梁, 吴小康, 甘世朝, 肖衡林, 马 强, 罗晨晨, . 不同初始孔隙比下非饱和黏土渗透性 试验研究及模型预测[J]. 岩土力学, 2019, 40(5): 1761-1770.
[3] 郑国锋, 郭晓霞, 邵龙潭, . 基于状态曲面的非饱和土强度准则及其验证[J]. 岩土力学, 2019, 40(4): 1441-1448.
[4] 汪华斌, 李建梅, 金怡轩, 周 博, 周 宇, . 降雨诱发边坡破坏数值模拟两个关键问题 的解决方法[J]. 岩土力学, 2019, 40(2): 777-784.
[5] 李明玉, 孙文静. 黏土掺入生物炭后的持水特性及其影响机制[J]. 岩土力学, 2019, 40(12): 4722-4730.
[6] 陶高梁,李 进,庄心善,肖衡林,崔惜琳,徐维生. 利用土中水分蒸发特性和微观孔隙分布规律确定SWCC残余含水率[J]. , 2018, 39(4): 1256-1262.
[7] 毕 骏,谌文武,戴鹏飞,林高潮, . 校正系数对不同形式的Van Genuchten方程各拟合参数的影响[J]. , 2018, 39(4): 1302-1310.
[8] 陶高梁,柏 亮,袁 波,甘世朝. 土-水特征曲线与核磁共振曲线的关系[J]. , 2018, 39(3): 943-948.
[9] 黎澄生,孔令伟,柏 巍,安 然,李甜果, . 土-水特征曲线滞后阻塞模型[J]. , 2018, 39(2): 598-604.
[10] 刘星志,刘小文,陈 铭,谷明晗. 基于3个不等粒径颗粒接触模型的土-水特征曲线[J]. , 2018, 39(2): 651-656.
[11] 方瑾瑾,冯以鑫,邵生俊,. 真三轴条件下原状Q3黄土的土-水特征[J]. , 2017, 38(9): 2597-2604.
[12] 毕 骏,谌文武, . 参数m的形式和最大基质吸力点对Van Genuchten方程各参数的影响[J]. , 2017, 38(12): 3604-3612.
[13] 蒋泽锋 ,朱大勇,. 强降雨条件下具有张裂缝边坡临界滑动场[J]. , 2016, 37(S2): 25-34.
[14] 黄启迪,赵成刚,蔡国庆,. 基于热力学并考虑体变及滞回效应的土-水特征曲线模型[J]. , 2016, 37(7): 1857-1867.
[15] 曹春山 ,吴树仁 ,潘 懋,梁昌玉,. 工程切坡诱发黄土滑坡成因机制研究[J]. , 2016, 37(4): 1049-1060.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘清秉,项 伟,张伟锋,崔德山. 离子土壤固化剂改性膨胀土的试验研究[J]. , 2009, 30(8): 2286 -2290 .
[2] 况雨春,伍开松,杨迎新,马德坤. 三牙轮钻头破岩过程计算机仿真模型[J]. , 2009, 30(S1): 235 -238 .
[3] 杜文琪,王 刚. 土工结构地震滑动位移统计分析[J]. , 2011, 32(S1): 520 -0525 .
[4] 鄢治华,刘志伟,刘厚健. 黄河阶地上某电厂高边坡参数选取及其工程治理[J]. , 2009, 30(S2): 465 -468 .
[5] 许振浩 ,李术才 ,李利平 ,侯建刚 ,隋 斌 ,石少帅. 基于层次分析法的岩溶隧道突水突泥风险评估[J]. , 2011, 32(6): 1757 -1766 .
[6] 江 权 ,冯夏庭 ,周 辉 ,赵 阳 ,徐鼎平 ,黄 可 ,江亚丽. 层间错动带的强度参数取值探讨[J]. , 2011, 32(11): 3379 -3386 .
[7] 温世清 ,刘汉龙 ,陈育民. 浆固碎石桩单桩荷载传递特性研究[J]. , 2011, 32(12): 3637 -3641 .
[8] 李顺群 ,高凌霞 ,柴寿喜. 冻土力学性质影响因素的显著性和交互作用研究[J]. , 2012, 33(4): 1173 -1177 .
[9] 钟 声 ,王川婴 ,吴立新 ,唐新建 ,王清远. 点状不良地质体钻孔雷达响应特征 ——围岩及充填效应正演分析[J]. , 2012, 33(4): 1191 -1195 .
[10] 孟 振,陈锦剑,王建华,尹振宇. 砂土中螺纹桩承载特性的模型试验研究[J]. , 2012, 33(S1): 141 -145 .