›› 2010, Vol. 31 ›› Issue (S2): 383-389.

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

复杂水力路径下非饱和流动的数值分析

陈 盼,韦昌富,李 幻,陈 辉,魏厚振   

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

    国家自然科学基金(No. 19872211);岩土力学与工程国家重点实验室重点方向性项目(Q110601);中国科学院“百人计划”择优项目资助

Numerical analysis of unsaturated flow in porous media under complicated hydraulic path

CHEN Pan, WEI Chang-fu, LI Huan, CHEN Hui, WEI Hou-zhen   

  1. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
  • Received:2010-07-29 Online:2010-12-10 Published:2011-01-12

摘要:

多孔介质中的流动问题,与孔隙介质的特性,含水量状态以及含水量的变化历史密切相关。基于毛细循环滞回理论模型,考虑含水量变化历史对土水特征关系的影响,在开发的U-DYSAC2有限元程序中进行了相应的数值实施。在试验给定的初边值条件下进行了非饱和渗流模拟分析,并将模拟结果与实测数据比较,表明在压力边界条件反复变化下,考虑滞回效应能获得更接近实测的结果,证实该模型在模拟各种循环变化条件下非饱和土渗流初边值问题的适用性与必要性。对入渗重分布反复变化条件下非饱和土柱流动的数值模拟表明,考虑滞回与不考虑滞回条件下,含水量、孔隙水压力和湿峰的迁移的预测在入渗后的重分布过程差异较大。考虑滞回效应时,土柱上部的脱湿速率、下部的吸湿速率比不考虑滞回时要低。从而证实了非饱和多孔介质中的土水状态依赖于含水量变化,而且强烈依赖于土体的水力路径变化。因此,循环边界条件变化下,毛细滞回效应在非饱和渗流模拟中的影响显著,必须加以考虑。

关键词: 渗流分析, 毛细滞回, 非饱和流动, 有限元方法, 数值模拟

Abstract:

Unsaturated flow problem in a porous medium is closely related to the characteristics of porous medium, soil moisture state and its drying/wetting history. Based on capillary cyclic hysteresis model which considers the effect of the history of water content on the relation of soil water characteristics, a modified U-DYSAC2 finite element program is developed to simulate the unsaturated seepage problems. Unsaturated flow analysis is done under experimental conditions. Compared the simulated results with experimental data from the literature, the results with hysteresis effect are much nearer to the data from the tests undergoing repeating changes of pressure boundary conditions. Hence the applicability and necessity of the numerical model is verified being capable for modeling the intial boundary problems of unsaturated seepage. Subsequently, an example is shown for flow analysis in a soil column subjected to complicated infiltration redistribution history. The results show that the discrepancy of water content, pore water pressure and movement rate of wetting front in the column is large during redistribution process considering hysteresis effect and non-hysteresis effect respectively. It is shown that rate of drying in uppper layer and rate of wetting in lower layer considering hysteresis effect than that of nonhysteresis analysis. The simulated results confirm that soil moisture state in unsaturated porous media is not only dependent on water content, but also related to hydraulic path in porous medium. Hence, under cyclic change of boundary conditions, it is necessary to consider hysteresis effect in modeling unsaturated seepage problems.

Key words: seepage analysis, capillary hysteresis, unsaturated flow, finite element method, numerical simulation

中图分类号: 

  • TD 324
[1] 李翻翻, 陈卫忠, 雷江, 于洪丹, 马永尚, . 基于塑性损伤的黏土岩力学特性研究[J]. 岩土力学, 2020, 41(1): 132-140.
[2] 夏 坤, 董林, 蒲小武, 李璐, . 黄土塬地震动响应特征分析[J]. 岩土力学, 2020, 41(1): 295-304.
[3] 郭院成, 李明宇, 张艳伟, . 预应力锚杆复合土钉墙支护体系增量解析方法[J]. 岩土力学, 2019, 40(S1): 253-258.
[4] 闫国强, 殷跃平, 黄波林, 张枝华, 代贞伟, . 三峡库区巫山金鸡岭滑坡成因机制与变形特征[J]. 岩土力学, 2019, 40(S1): 329-340.
[5] 刘红岩. 宏细观缺陷对岩体力学特性及边坡稳定影响研究[J]. 岩土力学, 2019, 40(S1): 431-439.
[6] 金爱兵, 刘佳伟, 赵怡晴, 王本鑫, 孙浩, 魏余栋, . 卸荷条件下花岗岩力学特性分析[J]. 岩土力学, 2019, 40(S1): 459-467.
[7] 韩征, 粟滨, 李艳鸽, 王伟, 王卫东, 黄健陵, 陈光齐, . 基于HBP本构模型的泥石流动力过程SPH数值模拟[J]. 岩土力学, 2019, 40(S1): 477-485.
[8] 吴锦亮, 何吉, . 岩质边坡动态开挖模拟的复合单元模型[J]. 岩土力学, 2019, 40(S1): 535-540.
[9] 吴凤元, 樊赟赟, 陈剑平, 李军, . 基于不同侵蚀模型的高速崩滑碎屑 流动力过程模拟分析[J]. 岩土力学, 2019, 40(8): 3236-3246.
[10] 孙峰, 薛世峰, 逄铭玉, 唐梅荣, 张翔, 李川, . 基于连续损伤的水平井射孔-近井筒三维破裂模拟[J]. 岩土力学, 2019, 40(8): 3255-3261.
[11] 穆锐, 浦少云, 黄质宏, 李永辉, 郑培鑫, 刘 旸, 刘 泽, 郑红超, . 土岩组合岩体中抗拔桩极限承载力的确定[J]. 岩土力学, 2019, 40(7): 2825-2837.
[12] 金俊超, 佘成学, 尚朋阳. 基于应变软化指标的岩石非线性蠕变模型[J]. 岩土力学, 2019, 40(6): 2239-2246.
[13] 张 聪, 梁经纬, 阳军生, 曹 磊, 谢亦朋, 张贵金, . 堤坝脉动注浆浆液扩散机制及应用研究[J]. 岩土力学, 2019, 40(4): 1507-1514.
[14] 严 健, 何 川, 汪 波, 蒙 伟, . 高地温对隧道岩爆发生的影响性研究[J]. 岩土力学, 2019, 40(4): 1543-1550.
[15] 李世俊, 马昌慧, 刘应明, 韩玉珍, 张 彬, 张 嘎, . 离心模型试验与数值模拟相结合研究 采空边坡渐进破坏特性[J]. 岩土力学, 2019, 40(4): 1577-1583.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 于小军,施建勇,徐杨斌. 考虑各向异性的软黏土扰动状态本构模型[J]. , 2009, 30(11): 3307 -3312 .
[2] 高 玮. 基于蚁群聚类算法的岩石边坡稳定性分析[J]. , 2009, 30(11): 3476 -3480 .
[3] 林 杭,曹 平,李江腾,江学良,何忠明. 层状岩质边坡破坏模式及稳定性的数值分析[J]. , 2010, 31(10): 3300 -3304 .
[4] 徐新生,燕 彬. 群桩与筏板的竖向动力相互作用研究[J]. , 2009, 30(1): 281 -284 .
[5] 徐 彬,殷宗泽,刘述丽. 膨胀土强度影响因素与规律的试验研究[J]. , 2011, 32(1): 44 -50 .
[6] 杨 洋, 姚海林,卢 正. 蒸发条件下路基对气候变化的响应模型及影响因素分析[J]. , 2009, 30(5): 1209 -1214 .
[7] 沈 扬,闫 俊,张朋举,费仲秋. 主应力方向变化路径下等压固结粉土强度特性差异和能量评价方法研究[J]. , 2011, 32(S1): 118 -123 .
[8] 楚锡华. 颗粒材料数值样本的坐标排序生成技术[J]. , 2011, 32(9): 2852 -2855 .
[9] 王小平,夏 雄. 岩土类材料率相关性及硬化-软化特性模型研究[J]. , 2011, 32(11): 3283 -3287 .
[10] 林鲁生 ,蒋 刚 ,白世伟 ,刘祖德 . 土体抗剪强度参数取值的统计分析方法[J]. , 2003, 24(2): 277 -280 .