An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model
-
摘要:
降雨条件下,边坡土体吸水软化造成抗剪强度降低是引发边坡失稳的重要因素之一。雨水在边坡土体中的入渗过程是动态变化的,土体含水率(饱和度)也会随之发生持续变化,因此,研究土体抗剪强度在动态环境下的变化规律,对实时评价降雨工况下土体边坡的稳定性具有重要意义。文章以贵州省尖山营地区的膨胀土为研究对象,基于非饱和土强度理论,结合室内剪切试验数据分析了膨胀土的抗剪强度指标与土体含水率的关系,提出了以含水率为变量的抗剪强度公式;基于ABAQUS有限元软件平台的二次开发技术编写了USDFLD子程序,结合室内试验得出的土体抗剪强度变化规律,分析了降雨条件下土体含水率变化引起的强度衰减与土体重度变化对膨胀土边坡稳定性的影响。结果表明:随膨胀土含水率增加,土中颗粒吸水膨胀,颗粒间距缩小,矿物颗粒间的黏结作用逐渐增大,直至集合体空隙充水,吸力减小,黏结作用迅速减弱;膨胀土的抗剪强度总体呈递减趋势,其中,黏聚力先增后减,内摩擦角递减。应用ABAQUS有限元软件的二次开发技术可实现土体参数随含水率的动态变化而变化,为降雨条件下非饱和土边坡的有限元分析提供了一种新的技术手段。
Abstract:Under the condition of rainfall, the decrease in shear strength of a slope is one of the important factors that lead to the slope instability. Infiltration of rainwater into the slope is a dynamic process, and the water content or saturation of soil will change continuously. Therefore, it is of great significance to study the shear strength of soil under dynamic environment for real-time evaluation of the stability of soil slope under rainfall conditions. Based on the strength theory of unsaturated soil, this paper takes the expansive soil in the Jianshanying area of Guizhou Province as the object, combines the shear test data to establish the fitting relationship between the shear strength index and water content of the expansive soil, and puts forward the shear strength formula with water content as the variable. Based on the secondary development technology of ABAQUS platform, USDFLD subroutine is compiled to import the shear strength formula of the Jianshanying expansive soil, and the influence of strength attenuation and gravity change caused by the change of soil moisture content on the stability of the expansive soil slope under the condition of rainfall are also explored. The results show that with the increasing water content of the expansive soil, the particles in the soil expand by the absorbing water, the particle spacing decreases, and the adhesion between the mineral particles increases gradually. Until the water filling suction of the aggregate space decreases, the adhesion decreases rapidly, and the shear strength of the expansive soil decreases generally, in which the cohesion increases first and then decreases, and the internal friction angle decreases. The secondary development technology of ABAQUS software can realize the relationship of soil parameters with the change of water content, which provides a new technical means for the finite element analysis of unsaturated soil slope under the condition of rainfall.
-
表 1 土体材料参数
Table 1. Soil material parameters
重度/
弹性模量/
MPa泊松比
有效黏聚力/
kPa有效内摩擦角/
(°)饱和渗透系数/
18 100 0.27 14.02 13 0.000369 表 2 Van-Genuchten模型参数
Table 2. Van-Genuchten model parameters
α n 0.4413 0.0648 0.002368 1.911 表 3 边坡稳定系数
Table 3. Slope stability coefficient
降雨历时/h 稳定系数 说明 导入子程序 不导入子程序 0 1.45 1.45 初始工况 12 − 0.994 塑性面贯通 21 0.997 − 24 0.927 0.845 失稳 -
[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