Material Point Method Simulation of Soil Landslide Considering Interface Friction Effect

Authors

• ZUO Zhipeng, Nanchang Urban Planning & Design Institute Group Co ., Ltd., Nanchang , Jiangxi 330038 , ChinaFollow
• HU Shengliang, Nanchang Urban Planning & Design Institute Group Co ., Ltd., Nanchang , Jiangxi 330038 , China
• LI Maowen, Nanchang Urban Planning & Design Institute Group Co ., Ltd., Nanchang , Jiangxi 330038 , China
• TONG Chenxi, School of Civil Engineering , Central South University , Changsha , Hunan 410075 , China
• LI Haichao, College of Transportation Science and Engineering , Civil Aviation University of China , Tianjin 300300 , ChinaFollow

Corresponding Author

李海潮，男，博士，讲师.E-mail：haichao821@outlook.com

Abstract

During the high-speed and long-runout sliding failure process of large landslides, severe friction and heat generation in the sliding zone induce an interface friction effect between the debris flow and the bedrock, influencing the final accumulation morphology of the landslide. In this paper, a two-dimensional landslide model of the Yangbaodi landslide in Shenzhen was established using the material point method. A contact algorithm was introduced through a multi-background grid method to conduct numerical simulations on the initiation, sliding, and accumulation processes of the landslide. The results show that after initiation, the landslide moves over a long distance in a flow-like state along the bedrock. Influenced by the topographic and geomorphic characteristics of the bedrock, local accumulation is formed in some gently inclined areas. The final accumulation morphology is relatively close to the calculation results of methods such as the particle finite element method, which preliminarily verifies the effectiveness of the method in this paper. Calculations changing the friction coefficient between the debris flow and the bedrock reveal that the interface friction effect significantly influences the movement behavior of the debris flow. The maximum sliding velocity and farthest sliding distance decrease as the friction coefficient increases. However, the soil deformation of the debris flow near the bedrock area is more obviously affected by friction, and the equivalent plastic strain conversely increases as the friction coefficient increases. At the same time, a large number of discontinuous plastic deformation protrusions are formed, influencing the final accumulation morphology of the debris flow. The research results can provide an effective technical means for analyzing the failure mechanism of natural landslides under the interface friction effect.

Publication Date

6-27-2026

DOI

10.14048/j.issn.1671-2579.2026.03.002

First Page

9

Last Page

18

Submission Date

June 2026

Recommended Citation

Zhipeng, ZUO; Shengliang, HU; Maowen, LI; Chenxi, TONG; and Haichao, LI (2026) "Material Point Method Simulation of Soil Landslide Considering Interface Friction Effect," Journal of China & Foreign Highway: Vol. 46: Iss. 3, Article 2.
DOI: 10.14048/j.issn.1671-2579.2026.03.002
Available at: https://zwgl1980.csust.edu.cn/journal/vol46/iss3/2

Reference

[1]中国自然资源报.去年我国共成功预报地质灾害 534起[EB/OL ] 北京：中国自然资源，2021 -01-18.China Natural Resources News.China successfully predicted 534 geological disasters last year [EB/OL ],2021 -01-18.
[2]HE X,LIANG D,BOLTON M D.Run-out of cut-slope landslides:mesh-free simulations [J].Géotechnique,2018,68(1):50-63.
[3]LIN C,PASTOR M,YAGUE A,et al.A depth-integrated SPH model for debris floods:Application to Lo Wai (Hong Kong )debris flood of August 2005 [J].Géotechnique,2019,69(12):1035 -1055.
[4]吴兵,梁瑶,李杰,等.注浆微型钢管桩组合结构在公路滑坡处治中的应用 [J].中外公路,2025,45(4):17-24.WU Bing,LIANG Yao,LI Jie,et al.Application of grouted steel pipe micropile composite structure in highway landslide treatment [J].Journal of China & Foreign Highway,2025,45(4):17-24.
[5]芮国荣,许建聪,孙钧.软岩公路隧道下穿古滑坡稳定性与风险分析 [J].公路,2024,69(3):399-404.RUI Guorong,XU Jiancong,SUN Jun.Stability and risk analysis of soft rock highway tunnel passing through ancient landslide [J].Highway,2024,69(3):399-404.
[6]ZHANG W J,ZHENG H,JIANG F Y,et al.Stability analysis of soil slope based on a water-soil-coupled and parallelized smoothed particle hydrodynamics model [J].Computers and Geotechnics,2019,108:212-225.
[7]黄洪超.某高速公路巨型滑坡稳定性分析及处置方案研究[J].中外公路,2023,43(4):24-29.HUANG Hongchao.Study on stability analysis and treatment of a expressway landslide [J].Journal of China & Foreign Highway,2023,43(4):24-29.
[8]许强,黄润秋.5⋅12汶川大地震诱发大型崩滑灾害动力特征初探 [J].工程地质学报,2008,16(6):721-729.XU Qiang,HUANG Runqiu.Kinetics charateristics of large landlides triggered by may 12th Wenchuan earthquake [J].Journal of Engineering Geology,2008,16(6):721-729.
[9]张卫杰,余瑞华,陈宇,等.强度指标影响下滑坡运动特征及参数反分析 [J].岩土工程学报,2022,44(12):2304 -2311,F0004.ZHANG Weijie,YU Ruihua,CHEN Yu,et al.Post-failure movement characteristics and parameter back-analysis for landslides considering effect of strength parameters [J].Chinese Journal of Geotechnical Engineering,2022,44(12):2304 -2311,F0004.
[10]沈佳轶,陈前,库猛,等.基于 CEL的灵敏性黏土斜坡渐进破坏数值模拟研究 [J].岩土工程学报,2022,44(12):2297 -2303.SHEN Jiayi,CHEN Qian,KU Meng,et al.Numerical simulation of progressive failure of sensitive clay slopes using CEL method [J].Chinese Journal of Geotechnical Engineering,2022,44(12):2297 -2303.
[11]LI W C,LI H J,DAI F C,et al.Discrete element modeling of a rainfall-induced flowslide [J].Engineering Geology,2012,149:22-34.
[12]ZHANG X,KRA BBENHOFT K,SHENG D C,et al.Numerical simulation of a flow-like landslide using the particle finite element method [J].Computational Mechanics,2015,55(1):167-177.
[13]SULSKY D,ZHOU S J,SCHREYER H L.Application of a particle-in-cell method to solid mechanics [J].Computer Physics Communications,1995,87(1/2):236-252.
[14]SULSKY D,KAUL A.Implicit dynamics in the material-point method [J].Computer Methods in Applied Mechanics and Engineering,2004,193(12/13/14):1137 -1170.
[15]SULSKY D,GONG M.Improving the material-point method [M/OL ]//WEINBERG K,WALLIN M.Innovative numerical approaches for multi-field and multi-scale problems.Cham:Springer International Publishing,2016:217-240.https://doi.org/10.1007 /978-3-319-39022 -2_10.
[16]LI X P,TANG X,ZHAO S X,et al.MPM evaluation of the dynamic runout process of the giant Daguangbao landslide [J].Landslides,2021,18(4):1509 -1518.
[17]FERNÁNDEZ F,VARGAS E,MULLER A L,et al.Material point method modeling in 3D of the failure and run-out processes of the Daguangbao landslide [J].Acta Geotechnica,2024,19(7):4277 -4296.
[18]ZHANG Y F,ZHANG W G,WANG L Q,et al.Mechanism of the high-speed and long-run-out landslide considering the evolution of the frictional heat in the sliding zone [J].Natural Hazards,2024,120(4):3299 -3317.
[19]朱晨光,刘春,许强,等.滑坡滑带摩擦热离散元数值模拟研究 [J].工程地质学报,2019,27(3):651-658.ZHU Chenguang,LIU Chun,XU Qiang,et al.Discrete element numerical simulation research on friction heat in sliding zone of the landslide [J].Journal of Engineering Geology,2019,27(3):651-658.
[20]周翠英,刘镇,薛翊国,等.关于红层灾变基础研究的若干思考 [J].工程地质学报,2023,31(3):689-705.ZHOU Cuiying,LIU Zhen,XUE Yiguo,et al.Some thoughts on basic research of red beds disaster [J].Journal of Engineering Geology,2023,31(3):689-705.
[21]BARDENHAGEN S G,KOBER E M.The generalized interpolation material point method [J].Computer Modeling in Engineering and Sciences,2004,5(6):477-496.
[22]BARDENHAGEN S G,BRACKBILL J U,SULSKY D.The material-point method for granular materials [J].Computer Methods in Applied Mechanics and Engineering,2000,187(3/4):529-541.
[23]MA J,WANG D,RANDOLPH M F.A new contact algorithm in the material point method for geotechnical simulations [J].International Journal for Numerical and Analytical Methods in Geomechanics,2014,38(11):1197 -1210.
[24]赵春宏,戴福初.深圳某填土滑坡破坏机理研究 [J].中国地质灾害与防治学报,2007,18(2):1-8.ZHAO Chunhong,DAI Fuchu.Study on failure mechanism of a fill slope in Shenzhen [J].The Chinese Journal of Geological Hazard and Control,2007,18(2):1-8.