Slope stability analysis using the virtual element method and shear strength reduction technique
LIN Shan1，GUO Yukui2，SUN Guanhua1，YANG Yongtao1，ZHANG Guohua1
(1. State Key Laboratory of Geomechanics and Geotechnical Engineering，Institute of Rock and Soil Mechanics，Chinese Academy
of Sciences，Wuhan，Hubei 430071，China；2. CCCC Second Highway Consultants Co.，Ltd.，
Wuhan，Hubei 430056，China)
Abstract：A new procedure for elasto-plastic mechanical problems with VEM is proposed，by considering the specialties of incremental elasto-plastic calculation and bilinear projection operator. The shear strength reduction technique together with the inequality is implemented into the VEM to investigate the stability of slopes. Based on the proposed procedure，the deformation of homogeneous and heterogeneous slopes under different shear strength reduction coefficients is solved. In addition，the mesh dependency of the model is tested. A series of numerical examples are employed to test the correctness and effectiveness of the method. The results show that：(1) the arbitrary polygonal elements can be adopted and the calculation accuracy is independent of the shape of the element. The VEM is simple to implement and has a reliable mathematical foundation. (2) The proposed procedure is able to accurately predict the factor of safety of slopes，and capture the failure patterns of slopes. In addition，it can provide evaluation basis for slopes(landslides) design.
[1] 尉迟衍春. 复杂边坡多失效模式的确定方法及其工程应用[硕士学位论文][D]. 青岛：青岛理工大学，2018.(YUCHI Yanchun. Method for determining multiple failure modes of complex slope and its engineering application[M. S. Thesis][D]. Qingdao：Qingdao University of Technology，2018.(in Chinese))
[2] 郑颖人. 岩土数值极限分析方法的发展与应用[J]. 岩石力学与工程学报，2012，31(7)：1 297–1 316.(ZHENG Yingren. Development and application of numerical limit analysis for geological materials[J]. Chinese Journal of Rock Mechanics and Engineering，2012，31(7)：1 297–1 316.(in Chinese))
[3] ZHENG H，LIU D F，LI C G. Slope stability analysis based on elasto-plastic finite element method[J]. International Journal for Numerical Methods in Engineering，2010，64(14)：1 871–1 888.
[4] 陈国庆，黄润秋，石豫川，等. 基于动态和整体强度折减法的边坡稳定性分析[J]. 岩石力学与工程学报，2014，33(2)：243–256. (CHEN Guoqing，HUANG Runqiu，SHI Yuchuan，et al. Stability analysis of slope based on dynamic and whole strength reduction methods[J]. Chinese Journal of Rock Mechanics and Engineering，2014，33(2)：243–256.(in Chinese))
[5] 王勖成，邵 敏. 有限单元法基本原理和数值方法[M]. 北京：清华大学出版社有限公司，2003：1–776.(WANG Xucheng，SHAO Min. Fundamental and numerical methods of finite element method[M]. Beijing：Tsinghua University Press，2003：1–776.(in Chinese))
[6] ZIENKIEWICZ O C H C L R W. Associated and non-associated visco-plasticity and plasticity in soil mechanics[J]. Geostechnique，1975，25(4)：671–689.
[7] GRIFFITHS D，LANE P. Slope stability analysis by finite elements[J]. Geotechnique，1999，49(3)：387–403.
[8] MATSUI T，SAN K-C. Finite element slope stability analysis by shear strength reduction technique[J]. Soils and Foundations，1992，32(1)：59–70.
[9] MANZARI M T，NOUR M A. Significance of soil dilatancy in slope stability analysis[J]. Journal of Geotechnical and Geoenvironmental Engineering，2000，126(1)：75–80.
[10] MUKHLISIN M，BAIDILLAH M R，IBRAHIM A，et al. Effect of soil hydraulic properties model on slope stability analysis based on strength reduction method[J]. Journal of the Geological Society of India，2014，83(5)：586–594.
[11] 袁 维，李小春，王 伟，等. 一种双折减系数的强度折减法研究[J]. 岩土力学，2016，37(8)：2 222–2 230.(YUAN Wei，LI Xiaochun，WANG Wei，et al. Study on strength reduction method based on double reduction parameters[J]. Rock and Soil Mechanics，2016，37(8)：2 222–2 230.(in Chinese))
[12] YUAN W，BAI B，LI X C，et al. A strength reduction method based on double reduction parameters and its application[J]. Journal of Central South University，2013，20(9)：2 555–2 562.
[13] CHENG Y M，LANSIVAARA T，WEI W B. Two-dimensional slope stability analysis by limit equilibrium and strength reduction methods[J]. Computers and Geotechnics，2007，34(3)：137–150.
[14] 宋二祥. 土工结构安全系数的有限元计算[J]. 岩土工程学报，1997，19(2)：1–7.(SONG Erxiang. Finite element analysis of safety factor for soil structures[J]. Chinese Journal of Geotechnical Engineering，1997，19(2)：1–7.(in Chinese))
[15] 连镇营，韩国城，孔宪京. 强度折减有限元法研究开挖边坡的稳定性[J]. 岩土工程学报，2001，23(4)：407–411.(LIAN Zhenying，HAN Guocheng，KONG Xianjing，et al. Stability analysis of excavation by strength reduction FEM[J]. Chinese Journal of Geotechnical Engineering，2001，23(4)：407–411.(in Chinese))
[16] 周翠英，刘祚秋，董立国，等. 边坡变形破坏过程的大变形有限元分析[J]. 岩土力学，2003，24(4)：644–652.(ZHOU Cuiying，LIU Zuoqiu，DONG Liguo，et al. Large deformation FEM analysis of slopes failure[J]. Rock and Soil Mechanics，2003，24(4)：644–652.(in Chinese))
[17] 栾茂田，武亚军，年廷凯. 强度折减有限元法中边坡失稳的塑性区判据及其应用[J]. 防灾减灾工程学报，2003，23(3)：1–8.(LUAN Maotian，WU Yajun，NIAN Tingkai，et al. A criterion for evaluating slope stability based on development of plastic zone by shear strength reduction FEM[J]. Journal of Disaster Prevention and Mitigation Engineering，2003，23(3)：1–8.(in Chinese))
[18] 冯 慧. 时空间稳定节点积分算法及其在磁–力耦合问题中的应用研究[博士学位论文][D]；长沙：湖南大学，2016.(FENG Hui. A temporal stable nodal integration method and the application in magnetic-mechanical coupling problem[Ph. D. Thesis][D]. Changsha：Changsha University，2016.(in Chinese))
[19] 杨 洪. 膨胀岩边坡稳定性数值分析的若干问题研究[硕士学位论文][D]. 武汉：长江科学院，2007.(YANG Hong. Study on several problems of numerical analysis of swelling rocks slope stability[M. S. Thesis][D]. Wuhan：Changjiang River Scientific Research Institute，2007.(in Chinese))
[20] 庄晓莹. 基于无网格方法的节理岩体边坡稳定性分析研究[硕士学位论文][D]. 上海：同济大学，2007.(ZHUANG Xiaoying. A study on jointed rock slope stability analysis using meshless method[M. S. Thesis][D]. Shanghai：Tongji University，2007.(in Chinese))
[21] BENEDETTO M F，BERRONE S，SCIALò S. A globally conforming method for solving flow in discrete fracture networks using the Virtual Element Method[J]. Finite Elements in Analysis and Design，2016，109：23–36.
[22] BREZZI F，MARINI L D. Virtual element methods for plate bending problems[J]. Computer Methods in Applied Mechanics and Engineering，2013，253(1)：455–462.
[23] VACCA G，LOURENCO B D V. Virtual element methods for parabolic problems on polygonal meshes[J]. Computers and Mathematics with Applications，2015，31(6)：2 110–2 134.
[24] VEIGA L B D，MARINI L D. Virtual Elements for Linear Elasticity Problems[J]. Siam Journal on Numerical Analysis，2013，51(2)：794–812.
[25] CHEN L. Equivalence of weak Galerkin methods and virtual element methods for elliptic equations[J/OL]. https：//www.researchgate.net/ publication/273704128_Equivalence_of_Weak_Galerkin_Methods_and_Virtual_Element_Methods_for_Elliptic_Equations.
[26] 郑 宏，李春光，李焯芬，等. 求解安全系数的有限元法 [J]. 岩土工程学报，2002，24(5)：626–628(ZHENG Hong，LI Chunguang，LEE C F，et al. Finite element method for solving the factor of safety[J]. Chinese Journal of Geotechnical Engineering，2002，24(5)：626–628.(in Chinese)).
[27] ORTIZ-BERNARDIN A，RUSSO A，SUKUMAR N. Consistent and stable meshfree Galerkin methods using the virtual element decomposition[J]. International Journal for Numerical Methods in Engineering，2017，112(7)：655–684.
[28] WRIGGERS P，RUST W T，REDDY B D. A virtual element method for contact[J]. Computational Mechanics，2016，58(6)：1 039–1 050.
[29] GAIN A L，TALISCHI C，PAULINO G H. On the virtual element method for three-dimensional linear elasticity problems on arbitrary polyhedral meshes[J]. Computer Methods in Applied Mechanics and Engineering，2014，282：132–160.
[30] 李 宁，郭双枫，姚显春. 再论岩质高边坡稳定性分析方法[J]. 岩土力学，2018，39(8)：397–406.(LI Ning，GUO Shuangfeng，YAO Xianchun. Further study of stability analysis methods of high rock slopes[J]. Rock and Soil Mechanics，2018，39(8)：397–406.(in Chinese))
[31] 王 栋，年廷凯，陈煜淼. 边坡稳定有限元分析中的三个问题[J]. 岩土力学，2007，(11)：2 309–2 313.(WANG Dong，NIAN Tingkai，CHEN Yumiao. Three problems in slope stability analysis with finite element method[J]. Rock and Soil Mechanics，2007，28(11)：2 309–2 313.(in Chinese))
[32] 刘祚秋，周翠英，董立国，等. 边坡稳定及加固分析的有限元强度折减法[J]. 岩土力学，2005，26(4)：558–561.(LIU Zuoqiu，ZHOU Cuiying，DONG Liguo，et al. Slope stability and strengthening analysis by strength reduction FEM[J]. Rock and Soil Mechanics，2005，26(4)：558–561.(in Chinese))
[33] 赵尚毅，郑颖人，时卫民，等. 用有限元强度折减法求边坡稳定安全系数[J]. 岩土工程学报，2002，24(3)：343–346.(ZHAO Shangyi，ZHENG Yingren，SHI Weimin，et al. Analysis on safety factor of slop by strength reduction FEM[J]. Chinese Journal of Geotechnical Engineering，2002，24(3)：343–346.(in Chinese))
[34] 夏志皋. 塑性力学[M]. 上海：同济大学出版社，1991：1–268.(XIA Zhigao. Plastic mechanics[M]. Shanghai：Tongji University Press，1991：1–268.(in Chinese))
[35] ZHENG H，THAM L G，LIU D. On two definitions of the factor of safety commonly used in the finite element slope stability analysis[J]. Computers and Geotechnics，2006，33(3)：188–195.
[36] CHANDLER R J. The inclination of talus，arctic talus terraces，and other slopes composed of granular materials[J]. Journal of Geology，1973，81(1)：1–14.
[37] 李世海，汪远年. 三维离散元土石混合体随机计算模型及单向加载试验数值模拟[J]. 岩土工程学报，2004，26(2)：172–177.(LI Shihai，WANG Yuannian. Stochastic model and numerical simulation of uniaxial loading test for rock and soil blending by 3D-DEM[J]. Chinese Journal of Geotechnical Engineering，2004，26(2)：172–177.(in Chinese))
[38] CHEN W F. Limit analysis and soil plasticity[M]. Amsterdam：Elsevier，1975：1–606.
[39] 徐文杰，胡瑞林，岳中崎. 土–石混合体随机细观结构生成系统的研发及其细观结构力学数值试验研究[J]. 岩石力学与工程学报，2009，28(8)：1 652–1 665.(XU Wenjie，HU Ruilin，YUE Zhongqi. Development of random mesostructure generating system of soil-rock mixture and study of its mesostructural mechanics based on numerical test[J]. Chinese Journal of Rock Mechanics and Engineering，2009，28(8)：1 652–1 665.(in Chinese))
2019, Vol. 38 
