岩土力学 ›› 2019, Vol. 40 ›› Issue (2): 777-784.doi: 10.16285/j.rsm.2017.1457

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

降雨诱发边坡破坏数值模拟两个关键问题 的解决方法

汪华斌1,李建梅1,金怡轩2,周 博1,周 宇1   

  1. 1. 华中科技大学 土木工程与力学学院,湖北 武汉 430074;2. 武汉市城市建设投资开发集团有限公司,湖北 武汉 430063
  • 收稿日期:2017-07-11 出版日期:2019-02-11 发布日期:2019-02-19
  • 作者简介:汪华斌,男,1968年生,博士,教授,博士生导师,主要从事边坡工程、岩土体与地下结构相互作用等方面的研究
  • 基金资助:
    国家自然科学基金(No. 41372296,No. 41672267)。

The numerical methods for two key problems in rainfall-induced slope failure

WANG Hua-bin1, LI Jian-mei1, JIN Yi-xuan2, ZHOU Bo1, ZHOU Yu1   

  1. 1. School of Civil Engineering and Mechanics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China; 2. Wuhan Urban Construction Investment Development Group Co., Ltd., Wuhan, Hubei 430063, China
  • Received:2017-07-11 Online:2019-02-11 Published:2019-02-19
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (41372296, 41672267)

摘要: 研究降雨因素对边坡稳定性的影响,最有效的方法是室内外试验和数值模拟方法,最为重要的是设置与实际情况接近的降雨入渗边界与动态流量边界下初始孔隙水压力分布问题。为有效地解决这两方面的问题,首先将现有的降雨数值边界施加方法所得模拟结果与室内试验结果进行了对比分析,指出了传统方法的不足之处,并在此基础上提出了新的降雨数值边界条件,即在土体表面施加一层很薄的“空气单元”,表述的是一个动态的边界条件。降雨入渗边界施加在“空气单元”表面,而不直接施加在土体单元表面,应用该方法完成的入渗土柱数值模拟结果更加接近现场和试验测得的降雨入渗率,为有效模拟降雨诱发边坡破坏的边界条件提供准确、真实可行的方法。同时,由于相对渗透系数和负孔隙水压力均为饱和度的函数,通过设置流体的抗拉强度来赋予节点负孔隙水压力模拟基质吸力,在数值模拟中每一个计算步根据节点饱和度自动更新负孔隙水压力和相对渗透系数,可实现FLAC软件中非饱和土的有效渗流分析。此外,为解决非饱和区土体施加负孔隙水压力所带来的问题,提出数值模拟过程中施加动态的流量边界,通过模拟自然界中地下水的补给,很好地解决使用FLAC进行数值模拟时动态流量边界条件下初始孔隙水压力的分布问题。

关键词: 降雨入渗边界, 初始孔隙水压力, 空气单元, 动态流量边界, 渗流分析, 数值模拟

Abstract: The effect of rainfall on the failure of slope can be analysed through experimental testing, field investigation, and numerical simulations. In numerical simulation for rainfall-triggered slope failure, how accurate boundary conditions of infiltration are initialized and initial pore water pressure is obtained in a dynamic flowing boundary condition are the most important aspects. To solve these problems, firstly, disadvantages of traditional numerical boundary in infiltration are found through numerical simulation and experimental results. Thus, new boundary conditions of infiltration are presented, in which the rainfall infiltration is processed through air units, as a dynamic boundary condition, and then into unsaturated soils. Considering these new infiltration conditions, a variety of numerical simulation is carried out in a soil column test for infiltration. It is shown that the results from simulation are consistent with those from experimental tests and field monitoring by others. The proposed method is proved to be effective and accurate for the rainfall-triggered slope failure. Meanwhile, the tensile strength for water at nodes is set to express suction in unsaturated soils. Because the relative permeability coefficient and suction are the function of saturation, the changes of saturation in each node can update the relative permeability coefficient and suction in each step of simulation. Then the seepage is developed using FISH language in the FLAC software. Moreover, a dynamic flowing boundary condition is proposed to consider the presence of suction in unsaturated soils. In this condition, the natural rainfall infiltration is simulated through the change of groundwater, and the initial pore water is precisely obtained after the development of seepage modules in the FLAC software.

Key words: rainfall infiltration boundary, initial pore water pressure, air units, dynamic flowing boundary, seepage analysis, numerical simulation

中图分类号: 

  • TU 432
[1] 毛浩宇, 徐奴文, 李彪, 樊义林, 吴家耀, 孟国涛, . 基于离散元模拟和微震监测的白鹤滩水电站左岸地下厂房稳定性分析[J]. 岩土力学, 2020, 41(7): 2470-2484.
[2] 史林肯, 周辉, 宋明, 卢景景, 张传庆, 路新景, . 深部复合地层TBM开挖扰动模型试验研究[J]. 岩土力学, 2020, 41(6): 1933-1943.
[3] 张振, 张朝, 叶观宝, 王萌, 肖彦, 程义, . 劲芯水泥土桩承载路堤渐进式失稳破坏机制[J]. 岩土力学, 2020, 41(6): 2122-2131.
[4] 苏杰, 周正华, 李小军, 董青, 李玉萍, 陈柳. 基于偏振特性的下孔法剪切波到时判别问题探讨[J]. 岩土力学, 2020, 41(4): 1420-1428.
[5] 杨高升, 白冰, 姚晓亮, . 高含冰量冻土路基融化固结规律研究[J]. 岩土力学, 2020, 41(3): 1010-1018.
[6] 马秋峰, 秦跃平, 周天白, 杨小彬. 岩石剪切断裂面接触算法的开发与应用[J]. 岩土力学, 2020, 41(3): 1074-1085.
[7] 李康, 王威, 杨典森, 陈卫忠, 亓宪寅, 谭彩. 周期振荡法在低渗透测量中的应用研究[J]. 岩土力学, 2020, 41(3): 1086-1094.
[8] 李翻翻, 陈卫忠, 雷江, 于洪丹, 马永尚, . 基于塑性损伤的黏土岩力学特性研究[J]. 岩土力学, 2020, 41(1): 132-140.
[9] 夏 坤, 董林, 蒲小武, 李璐, . 黄土塬地震动响应特征分析[J]. 岩土力学, 2020, 41(1): 295-304.
[10] 郭院成, 李明宇, 张艳伟, . 预应力锚杆复合土钉墙支护体系增量解析方法[J]. 岩土力学, 2019, 40(S1): 253-258.
[11] 闫国强, 殷跃平, 黄波林, 张枝华, 代贞伟, . 三峡库区巫山金鸡岭滑坡成因机制与变形特征[J]. 岩土力学, 2019, 40(S1): 329-340.
[12] 刘红岩. 宏细观缺陷对岩体力学特性及边坡稳定影响研究[J]. 岩土力学, 2019, 40(S1): 431-439.
[13] 金爱兵, 刘佳伟, 赵怡晴, 王本鑫, 孙浩, 魏余栋, . 卸荷条件下花岗岩力学特性分析[J]. 岩土力学, 2019, 40(S1): 459-467.
[14] 韩征, 粟滨, 李艳鸽, 王伟, 王卫东, 黄健陵, 陈光齐, . 基于HBP本构模型的泥石流动力过程SPH数值模拟[J]. 岩土力学, 2019, 40(S1): 477-485.
[15] 吴锦亮, 何吉, . 岩质边坡动态开挖模拟的复合单元模型[J]. 岩土力学, 2019, 40(S1): 535-540.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!