Slope stability analysis based on a multigrid method using a nonlinear 3D finite element model

Abstract

The rigid-body limit equilibrium method cannot reflect the actual stress distribution in a rock mass, and the finite-element-based strength reduction method also has some problems with respect to convergence. To address these problems, a multi-grid method was adopted in this study to establish a structural grid for finite element computation and a slip surface grid for computing slope stability safety factors. This method can be used to determine the stability safety factor for any slip surface or slide block through a combination of nonlinear finite element analysis and limit equilibrium analysis. An ideal elastic–plastic incremental analysis method based on the Drucker–Prager yield criterion was adopted in the nonlinear finite element computation. Elasto-plastic computation achieves good convergence for both small load steps and large load steps and can increase computation precision to a certain extent. To increase the scale and accuracy of the computation, TFINE, a finite element parallel computation program, was used to analyze the influence of grid density on the accuracy of the computation results and was then applied to analysis of the stability of the Jinping high slope. A comparison of the results with results obtained using the rigid-body limit equilibrium method showed that the slope stability safety factors determined using finite element analysis were greater than those obtained using the rigid-body limit equilibrium method and were in better agreement with actual values because nonlinear stress adjustment was considered in the calculation.

Keywords

slope / stability / multi-grid method / nonlinear / finite element method

References

[1] Bishop A W. The use of the slip circle in the stability analysis of slopes. Geotechnique, 1955, 5(1): 7-17

[2] Duncan J M, 0. State of the art: limit equilibrium and finite-element analysis of slopes. Journal of Geotechnical Engineering, 1996, 122(7): 577-596

[3] Zuyu C, Xiaogang W, Jian Y, Zhixin J, Yujie W. Rock slope stability analysis: Theory methods and programs. Beijing: China Water Power Press, 2005 (in Chinese)

[4] Griffiths D V, Lane P A. Slope stability analysis by finite elements. Geotechnique, 1999, 49(3): 387-403

[5] Lenchman J B, Griffiths D V. Analysis of the progression of failure of the earth slopes by finite elements. In: Proceedings of Sessions of Geo-Denver 2000-Slope Stability 2000. GSP 101, 289: 250-265

[6] Liu Y R. Yang Q, Zhu L. Abutment stability analysis of arch dam based on 3D nonlinear finite element method. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(1): 3222-3228 (in Chinese)

[7] Jiang G L, Magnan J P. Stability analysis of embankments: Comparison of limit analysis with methods of slices. Geotechnique, 1997, 47(4): 857-872

[8] Dawson E M, Roth W H, Drescher A. Slope stability analysis by strength reduction. Geotechnique, 1999, 49(6): 835-840

[9] Zhao S Y, Zheng Y R, Shi W M, Wang J L. Analysis on safety factor of slope by strength reduction FEM. Chinese Journal of Geotechnical Engineering, 2002, 24(3): 343-346 (in Chinese)

[10] Yang Q, Zhu L, Xue L J. Application of limit equilibrium method to stability analysis of Jinping high slope based on 3D multi-grid method. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(Supp.2): 5313-5318 (in Chinese)

[11] Liu Y R, Yang Q, Xue L J, Zhou W Y. Rock Slop Stability Analysis with Nonlinear Finite Element Method. In: Cai M F, Wang J, eds. Boundaries of Rock Mechanics. London: Taylor & Francis Group, 2008, 503-507

[12] Hjiaj M, Fortin J, de Saxce G. A complete stress update algorithm for the non-associated Drucker-Prager model including treatment of the apex. International Journal of Engineering Science, 2003, 41(10): 1109-1143

[13] Yang Q, Yang X J, Chen X. On integration algorithms for perfect plasticity based on Drucker-Prager criterion. Engineering Mechanics, 2005, 22(4): 15-19, 47 (in Chinese)

[14] Yang Q, Chen X, Zhou W Y. A practical 3D elastic-plastic incremental method in FEM based on D-P yield criteria. Chinese Journal of Geotechnical Engineering, 2002, 24(1): 16-20 (in Chinese)

[15] Chen X. Yang Q, Huang Y S. Sub-incremental method for perfect elasto-plastic material based on D-P yield criteria. Chinese Journal of Rock Mechanics and Engineering, 2002, Suppl 2: 2465-2469 (in Chinese)

[16] Yang Q. Chen X, Zhou W Y. Elastic-plastic basis of geotechnical engineering reinforcement analysis. Rock and Soil Mechanics, 2005, 26(4): 553-557 (in Chinese)

[17] Schreyer H L, Kulak R F, Kramer J M. Accurate numerical solutions for elastoplastic models. Journals of Pressure Vessel Technology, 1979, 101(2): 226-234

[18] Ortiz M, Popov E P. Accuracy and stability of integration algorithms for elastoplastic constitutive relations. International Journal for Numerical Methods in Engineering, 1985, 21(9): 1561-1576

[19] Liu Y R, Zhou W Y, Yang Q. A distributed memory parallel element-by-element scheme based on Jacobi-conditioned conjugate gradient for 3-D finite element analysis. Finite Elements in Analysis and Design, 2007, 43(6-7): 494-503

Acknowledgements

The research reported in this paper was supported by the National Science Foundation of China (Grant Nos. 50823005 and 50709014) and subsidized by China National Funds for Distinguished Young Scientists (Grant No. 50925931).

RIGHTS & PERMISSIONS

2014 Higher Education Press and Springer-Verlag Berlin Heidelberg