基于有限元软件自定义本构模型的膨胀土边坡降雨入渗分析

饶鸿, 王金淑, 赵志明, 吴光, 冯涛. 基于有限元软件自定义本构模型的膨胀土边坡降雨入渗分析[J]. 水文地质工程地质, 2021, 48(1): 154-162. doi: 10.16030/j.cnki.issn.1000-3665.202002020
引用本文: 饶鸿, 王金淑, 赵志明, 吴光, 冯涛. 基于有限元软件自定义本构模型的膨胀土边坡降雨入渗分析[J]. 水文地质工程地质, 2021, 48(1): 154-162. doi: 10.16030/j.cnki.issn.1000-3665.202002020
RAO Hong, WANG Jinshu, ZHAO Zhiming, WU Guang, FENG Tao. An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model[J]. Hydrogeology & Engineering Geology, 2021, 48(1): 154-162. doi: 10.16030/j.cnki.issn.1000-3665.202002020
Citation: RAO Hong, WANG Jinshu, ZHAO Zhiming, WU Guang, FENG Tao. An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model[J]. Hydrogeology & Engineering Geology, 2021, 48(1): 154-162. doi: 10.16030/j.cnki.issn.1000-3665.202002020

基于有限元软件自定义本构模型的膨胀土边坡降雨入渗分析

  • 基金项目: 四川省科技计划项目资助(2019YFG0460);中铁二院工程集团有限责任公司资助(勘-16-贵南-31(061))
详细信息
    作者简介: 饶鸿(1996-),男,硕士研究生,主要从事工程地质和特殊岩土体工程研究。E-mail: 578341707@qq.com
    通讯作者: 赵志明(1977-),男,讲师,主要从事铁路工程地质及工程岩土体稳定性研究。E-mail: zmzhao@swjtu.edu.cn
  • 中图分类号: P642.13;P642.2

An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model

More Information
  • 降雨条件下,边坡土体吸水软化造成抗剪强度降低是引发边坡失稳的重要因素之一。雨水在边坡土体中的入渗过程是动态变化的,土体含水率(饱和度)也会随之发生持续变化,因此,研究土体抗剪强度在动态环境下的变化规律,对实时评价降雨工况下土体边坡的稳定性具有重要意义。文章以贵州省尖山营地区的膨胀土为研究对象,基于非饱和土强度理论,结合室内剪切试验数据分析了膨胀土的抗剪强度指标与土体含水率的关系,提出了以含水率为变量的抗剪强度公式;基于ABAQUS有限元软件平台的二次开发技术编写了USDFLD子程序,结合室内试验得出的土体抗剪强度变化规律,分析了降雨条件下土体含水率变化引起的强度衰减与土体重度变化对膨胀土边坡稳定性的影响。结果表明:随膨胀土含水率增加,土中颗粒吸水膨胀,颗粒间距缩小,矿物颗粒间的黏结作用逐渐增大,直至集合体空隙充水,吸力减小,黏结作用迅速减弱;膨胀土的抗剪强度总体呈递减趋势,其中,黏聚力先增后减,内摩擦角递减。应用ABAQUS有限元软件的二次开发技术可实现土体参数随含水率的动态变化而变化,为降雨条件下非饱和土边坡的有限元分析提供了一种新的技术手段。

  • 加载中
  • 图 1  尖山营膨胀土标准环刀试样

    Figure 1. 

    图 2  抗剪强度指标拟合曲线

    Figure 2. 

    图 3  计算流程图

    Figure 3. 

    图 4  膨胀土边坡计算模型

    Figure 4. 

    图 5  土水特征曲线

    Figure 5. 

    图 6  降雨过程含水率曲线

    Figure 6. 

    图 7  降雨前后含水率云图

    Figure 7. 

    图 8  土体位移曲线及总位移矢量图

    Figure 8. 

    图 9  边坡塑性贯通面示意图

    Figure 9. 

    表 1  土体材料参数

    Table 1.  Soil material parameters

    重度/
    弹性模量/
    MPa
    泊松比
    有效黏聚力/
    kPa
    有效内摩擦角/
    (°)
    饱和渗透系数/
    18 100 0.27 14.02 13 0.000369
    下载: 导出CSV

    表 2  Van-Genuchten模型参数

    Table 2.  Van-Genuchten model parameters

    α n
    0.4413 0.0648 0.002368 1.911
    下载: 导出CSV

    表 3  边坡稳定系数

    Table 3.  Slope stability coefficient

    降雨历时/h 稳定系数 说明
    导入子程序 不导入子程序
    0 1.45 1.45 初始工况
    12 0.994 塑性面贯通
    21 0.997
    24 0.927 0.845 失稳
    下载: 导出CSV
  • [1]

    李国荣, 陈文婷, 朱海丽, 等. 青藏高原东北部黄土地区降雨入渗对土质边坡稳定性的影响研究[J]. 水文地质工程地质,2015,42(2):105 − 111. [LI Guorong, CHEN Wenting, ZHU Haili, et al. Influence of rain infiltration to stability of loess slope in the northeast Qinghai-Tibet Plateau[J]. Hydrogeology & Engineering Geology,2015,42(2):105 − 111. (in Chinese with English abstract)

    [2]

    NG C W W, SHI Q. A numerical investigation of the stability of unsaturated soil slopes subjected to transient seepage[J]. Computers and Geotechnics,1998,22(1):1 − 28. doi: 10.1016/S0266-352X(97)00036-0

    [3]

    CAI F, UGAI K, WAKAI A, et al. Effects of horizontal drains on slope stability under rainfall by three-dimensional finite element analysis[J]. Computers and Geotechnics,1998,23(4):255 − 275. doi: 10.1016/S0266-352X(98)00021-4

    [4]

    韦立德, 陈从新, 徐健, 等. 考虑渗流和锚固作用的强度折减有限元法研究[J]. 岩石力学与工程学报,2008,27(增刊2):3471 − 3476. [WEI Lide, CHEN Congxin, XU Jia, et al. Study on strength reduction FEM considering seepage and bolt[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(Sup2):3471 − 3476. (in Chinese with English abstract)

    [5]

    陈善雄, 陈守义. 考虑降雨的非饱和土边坡稳定性分析方法[J]. 岩土力学,2001,22(4):447 − 450. [CHEN Shanxiong, CHEN Shouyi. Analysis of stability of unsaturated soil slope due to permeation of rainwater[J]. Rock and Soil Mechanics,2001,22(4):447 − 450. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-7598.2001.04.019

    [6]

    麻官亮, 邵玉刚. 基于Spencer法的三维极限平衡边坡稳定性研究[J]. 水文地质工程地质,2016,43(3):118 − 123. [MA Guanliang, SHAO Yugang. Study on three-dimensional limit equilibrium for slope stability based on Spencer method[J]. Hydrogeology & Engineering Geology,2016,43(3):118 − 123. (in Chinese with English abstract)

    [7]

    李宁, 刘冠麟, 许建聪, 等. 降雨条件下边坡有限元强度折减法计算平台开发及其应用[J]. 水文地质工程地质,2018,45(3):63 − 70. [LI Ning, LIU Guanlin, XU Jiancong, et al. Development of shear strength reduction method for slope stability analysis under rainfall conditions and its application[J]. Hydrogeology & Engineering Geology,2018,45(3):63 − 70. (in Chinese with English abstract)

    [8]

    王协群, 张有祥, 邹维列, 等. 降雨入渗条件下非饱和路堤变形与边坡的稳定数值模拟[J]. 岩土力学,2010,31(11):3640 − 3644. [WANG Xiequn, ZHANG Youxiang, ZOU Weilie, et al. Numerical simulation for unsaturated road-embankment deformation and slope stability under rainfall infiltration[J]. Rock and Soil Mechanics,2010,31(11):3640 − 3644. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-7598.2010.11.046

    [9]

    李绍红, 朱建东, 王少阳, 等. 考虑降雨类型的基岩型浅层边坡稳定性分析方法[J]. 水文地质工程地质,2018,45(2):131 − 135. [LI Shaohong, ZHU Jiandong, WANG Shaoyang, et al. Stability analysis methods for the bedrock shallow slope considering rainfall types[J]. Hydrogeology & Engineering Geology,2018,45(2):131 − 135. (in Chinese with English abstract)

    [10]

    张晨阳, 张泰丽, 张明, 等. 东南沿海地区玄武岩残积土雨水运移特征及滑坡失稳数值模拟[J]. 水文地质工程地质,2019,46(4):42 − 50. [ZHANG Chenyang, ZHANG Taili, ZHANG Ming, et al. Rainfall infiltration characteristics and numerical simulation of slope instability in the basalt residual soil in the coastal area of Southeast China[J]. Hydrogeology & Engineering Geology,2019,46(4):42 − 50. (in Chinese with English abstract)

    [11]

    BISHOP A W, BLIGHT G E. Some aspects of effective stress in saturated and partly saturated soils[J]. Géotechnique,1963,13(3):177 − 197. doi: 10.1680/geot.1963.13.3.177

    [12]

    FREDLUND D G, MORGENSTERN N R, WIDGER R A. The shear strength of unsaturated soils[J]. Canadian Geotechnical Journal,1978,15(3):313 − 321. doi: 10.1139/t78-029

    [13]

    GAN J K M, FREDLUND D G, RAHARDJO H. Determination of the shear strength parameters of an unsaturated soil using the direct shear test[J]. Canadian Geotechnical Journal,1988,25(3):500 − 510.

    [14]

    卢肇钧. 非饱和土抗剪强度的探索研究[J]. 中国铁道科学,1999,20(2):10 − 16. [LU Zhaojun. Explorations on the suctional shear strength of unsaturated soils[J]. China Railway Science,1999,20(2):10 − 16. (in Chinese with English abstract)

    [15]

    缪林昌, 仲晓晨, 殷宗泽. 膨胀土的强度与含水量的关系[J]. 岩土力学,1999,20(2):71 − 75. [MIAO Linchang, ZHONG Xiaochen, YIN Zongze. The relationship between strength and water content of expensive soil[J]. Rock and Soil Mechanics,1999,20(2):71 − 75. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-7598.1999.02.015

    [16]

    郭倩怡, 谷天峰, 吴熠哲. 永靖非饱和黄土抗剪强度试验研究[J]. 水文地质工程地质,2015,42(6):103 − 107. [GUO Qianyi, GU Tianfeng, WU Yizhe. A test study of shear strength of unsaturated loess in Yongjing County[J]. Hydrogeology & Engineering Geology,2015,42(6):103 − 107. (in Chinese with English abstract)

    [17]

    卢肇钧, 吴肖茗, 孙玉珍, 等. 膨胀力在非饱和土强度理论中的作用[J]. 岩土工程学报,1997,19(5):22 − 29. [LU Zhaojun, WU Xiaoming, SUN Yuzhen, et al. The role of swelling pressure in the shear strength theory of unsaturated soils[J]. Chinese Journal of Geotechnical Engineering,1997,19(5):22 − 29. (in Chinese with English abstract)

    [18]

    童富果, 田斌, 刘德富. 改进的斜坡降雨入渗与坡面径流耦合算法研究[J]. 岩土力学,2008,29(4):1035 − 1040. [TONG Fuguo, TIAN Bin, LIU Defu. A coupling analysis of slope runoff and infiltration under rainfall[J]. Rock and Soil Mechanics,2008,29(4):1035 − 1040. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-7598.2008.04.034

    [19]

    吴梦喜. 饱和-非饱和土中渗流Richards方程有限元算法[J]. 水利学报,2009,40(10):1274 − 1279. [WU Mengxi. Finite-element algorithm for Richards’ equation for saturated-unsaturated seepage flow[J]. Journal of Hydraulic Engineering,2009,40(10):1274 − 1279. (in Chinese with English abstract) doi: 10.3321/j.issn:0559-9350.2009.10.018

    [20]

    费康, 彭劼. ABAQUS岩土工程实例详解[M]. 北京: 人民邮电出版社, 2017: 173-174. [

    FEI Kang, PENG Jie. Detailed explanation of ABAQUS geotechnical engineering example[M]. Beijing: Posts & Telecom Press, 2017: 173-174. (in Chinese)]

    [21]

    VAN GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal,1980,44(5):892 − 898. doi: 10.2136/sssaj1980.03615995004400050002x

  • 加载中

(9)

(3)

计量
  • 文章访问数:  2351
  • PDF下载数:  96
  • 施引文献:  0
出版历程
收稿日期:  2020-02-16
修回日期:  2020-03-31
刊出日期:  2021-01-15

目录