岩土力学 ›› 2021, Vol. 42 ›› Issue (12): 3345-3355.doi: 10.16285/j.rsm.2021.0608

• 基础理论与实验研究 • 上一篇    下一篇

基于多相物质点法的土−水耦合动力模型

唐雄1, 2,李新坡1, 2,姚军1, 2,孙玉莲1, 2   

  1. 1. 中国科学院水利部成都山地灾害与环境研究所,四川 成都 610041;2. 中国科学院大学,北京 101499
  • 收稿日期:2021-04-21 修回日期:2021-07-05 出版日期:2021-12-13 发布日期:2021-12-14
  • 作者简介:唐雄,男,1997年生,硕士研究生,主要从事岩土工程、滑坡数值模拟方面的研究。
  • 基金资助:
    国家重点研发计划(No.2020YFC1512005);国家自然科学基金(No.41877291)。

Soil-water coupling dynamic model based on multiphase material point method

TANG Xiong1, 2, LI Xin-po1, 2, YAO Jun1, 2, SUN Yu-lian1, 2   

  1. 1. Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu, Sichuan 610041, China; 2. University of Chinese Academy of Sciences, Beijing 101499, China
  • Received:2021-04-21 Revised:2021-07-05 Online:2021-12-13 Published:2021-12-14
  • Supported by:
    This work was supported by the National Key R&D Program of China(2020YFC1512005) and the National Natural Science Foundation of China(41877291).

摘要: 作为无网格法的一种,物质点法(MPM)在计算中结合了有网格法和无网格法的优势,避免了有网格法在求解大变形问题时出现的网格畸变现象,近年来逐渐应用于大变形岩土工程问题和多相流问题中。基于标准物质点法的原理与算法,研究了多相物质点法的两种形式:适用于饱和土的两相单点格式;适用于非饱和土,考虑基质吸力的两相单点格式。为了验证多相物质点法模拟水力耦合多相流问题的准确性和适用性,模拟了Terzaghi一维固结问题、一维压缩波的传播问题和一维渗流问题3个土力学多相算例,并将所得到的数值解与理论解进行了对比分析。模拟结果表明,计算所得到的数值解与理论解基本保持一致,多相物质点法能够准确地描述饱和孔隙介质的水力特性和非饱和孔隙介质的渗流特性。最后,采用多相物质点法模拟了非饱和边坡降雨入渗问题,分析了边坡在降雨入渗作用下基质吸力的变化和大变形破坏全过程。

关键词: 多相物质点法, 两相单点格式, 水力耦合问题

Abstract: As a kind of meshless method, the material point method (MPM) combines the best features of the meshed method and meshless method in the calculation, and avoids the mesh distortion phenomenon of the meshed method in solving large deformation problems. In the last decades, MPM has been successfully applied in a number of large deformation geotechnical problems and multiphase flow problems. Based on the theory and algorithm of the standard material point method, two multiphase MPM formulations were studied: the two-phase single-point formulation for saturated soil, and the two-phase single-point formulation with suction for unsaturated soil. To examine the accuracy and applicability of the multiphase MPM when simulating coupled hydro-mechanical multiphase flow problems, Terzaghi’s one-dimensional consolidation problem, one-dimensional compression wave propagation problem and one-dimensional seepage problem were simulated and the results were compared to the analytical solution. The simulation results show that the numerical solution is generally consistent with the theoretical solution, and the multiphase MPM can accurately describe the hydraulic characteristics of saturated porous media and the seepage characteristics of unsaturated porous media. Finally, the multiphase MPM was used to simulate the problem of rainfall infiltration in unsaturated slope, and the change of matrix suction and the whole process of large deformation failure of the slope under rainfall infiltration were analyzed.

Key words: multiphase material point method, two-phase single-point formulation, hydro-mechanical coupling problem

中图分类号: 

  • TU452
[1] 王磊, 刘怀谦, 谢广祥, 袁秋鹏, 陈礼鹏, . 含瓦斯煤孔裂隙结构精细表征及强度劣化机制[J]. 岩土力学, 2021, 42(12): 3203-3216.
[2] 李春元, 左建平, 张勇, . 深部开采底板破坏与基本顶岩梁初次垮断的 联动效应[J]. 岩土力学, 2021, 42(12): 3301-3314.
[3] 周超彪, 刘东, 姜清辉, . 拉剪作用下类岩石试样的力学 特性与损伤破坏机制[J]. 岩土力学, 2021, 42(12): 3335-3344.
[4] 杨建华, 章伟鹏, 姚池, 张小波, 周创兵. 基于变分模态分解的深埋洞室爆破开挖围岩 振动特性分析[J]. 岩土力学, 2021, 42(12): 3366-3375.
[5] 陈念, 王成虎, 高桂云, 王璞, . 基于应力多边形与钻孔崩落的白鹤滩 右岸厂房区地应力场特征研究[J]. 岩土力学, 2021, 42(12): 3376-3384.
[6] 杨亮, 杨永涛, 郑宏, . 相场数值流形法模拟岩石裂纹扩展[J]. 岩土力学, 2021, 42(12): 3419-3427.
[7] 宋磊博, 亢倩倩, 杜时贵, 钟振, 王刚, 王兴开, 韩观胜, 赵金帅, . 基于表面微凸体磨损特征演化规律的节 理剪切强度各向异性机制研究[J]. 岩土力学, 2021, 42(9): 2331-2343.
[8] 张超, 杨楚卿, 白允. 岩石类脆性材料损伤演化分析及其模型方法研究[J]. 岩土力学, 2021, 42(9): 2344-2354.
[9] 徐鼎平, 郭广涛, 夏跃林, 柳秀洋, 江权, 李邵军, 李治国, . 高应力强卸荷下双江口花岗岩岩爆中间 主应力效应宏细观试验研究[J]. 岩土力学, 2021, 42(9): 2375-2386.
[10] 蔡灿, 张沛, 孙明光, 杨迎新, 谢松, 蒲治成, 杨显鹏, 高超, 谭政博, . 油气钻井中的分离式冲击−切削复合破岩机制研究[J]. 岩土力学, 2021, 42(9): 2535-2544.
[11] 张建聪, 江权, 郝宪杰, 丰光亮, 李邵军, 汪志林, 樊启祥, . 高应力下柱状节理玄武岩应力−结 构型塌方机制分析[J]. 岩土力学, 2021, 42(9): 2556-2568.
[12] 杜崧, 肖明, 陈俊涛, . 洞室块体危险性分析的突变级数评价法研究[J]. 岩土力学, 2021, 42(9): 2578-2588.
[13] 蒋水华, 欧阳苏, 冯泽文, 康青, 黄劲松, 杨志刚, . 基于结构面参数概率分布更新的 节理岩质边坡可靠性分析[J]. 岩土力学, 2021, 42(9): 2589-2599.
[14] 陈世杰, 肖明, 王小威, 陈俊涛, . 陡倾角断层下地下洞室地震破坏特性数值分析[J]. 岩土力学, 2021, 42(9): 2600-2610.
[15] 贾蓬, 杨其要, 刘冬桥, 王述红, 赵永, . 高温花岗岩水冷却后物理力学特性及微观破裂特征[J]. 岩土力学, 2021, 42(6): 1568-1578.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘 晓,唐辉明,刘 瑜. 基于集对分析和模糊马尔可夫链的滑坡变形预测新方法研究[J]. , 2009, 30(11): 3399 -3405 .
[2] 胡大伟,周 辉,谢守益,张 凯,邵建富,冯夏庭. 大理岩破坏阶段Biot系数研究[J]. , 2009, 30(12): 3727 -3732 .
[3] 师旭超,韩 阳. 卸荷作用下软黏土回弹吸水试验研究[J]. , 2010, 31(3): 732 -736 .
[4] 原喜忠,李 宁,赵秀云,杨银涛. 东北多年冻土地区地基承载力对气候变化敏感性分析[J]. , 2010, 31(10): 3265 -3272 .
[5] 白 冰,李小春,石 露,唐礼忠. 弹塑性应力-应变曲线的斜率恒等式及其验证和应用[J]. , 2010, 31(12): 3789 -3792 .
[6] 唐利民. 地基沉降预测模型的正则化算法[J]. , 2010, 31(12): 3945 -3948 .
[7] 蔡辉腾,危福泉,蔡宗文. 重庆主城区粉质黏土动力特性研究[J]. , 2009, 30(S2): 224 -228 .
[8] 宋 玲 ,刘奉银 ,李 宁 . 旋压入土式静力触探机制研究[J]. , 2011, 32(S1): 787 -0792 .
[9] 金解放 ,李夕兵 ,殷志强 ,邹 洋. 循环冲击下波阻抗定义岩石损伤变量的研究[J]. , 2011, 32(5): 1385 -1393 .
[10] 周延军 ,耿应春 ,王贵宾 ,唐洪林 ,李祖奎. 深部地层岩石力学性质测试与分析研究[J]. , 2011, 32(6): 1625 -1630 .