Geotechnical stability analysis considering strain softening using micro-polar continuum finite element method

Abstract

Geotechnical stability analyses based on classical continuum may lead to remarkable underestimations on geotechnical safety. To attain better estimations on geotechnical stability, the micro-polar continuum is employed so that its internal characteristic length (4) can be utilized to model the shear band width. Based on two soil slope examples, the role of internal characteristic length in modeling the shear band width of geomaterial is investigated by the second-order cone programming optimized micro-polar continuum finite element method. It is recognized that the underestimation on factor of safety (FOS) calculated from the classical continuum tends to be more pronounced with the increase of 4. When the micro-polar continuum is applied, the shear band dominated by 4 is almost kept unaffected as long as the adopted meshes are fine enough, but it does not generally present a slip surface like in the cases from the classical continuum, indicating that the micro-polar continuum is capable of capturing the non-local geotechnical failure characteristic. Due to the coupling effects of 4 and strain softening, softening behavior of geomaterial tends to be postponed. Additionally, the bearing capacity of a geotechnical system may be significantly underestimated, if the effects of 4 are not modeled or considered in numerical analyses.

Keywords

slope stability / strain localization / non-local geotechnical failure / micro-polar continuum / internal characteristic length

References

[1] GriffithsD V, laneP A. Slope stability analysis by finite elements [J]. Géotechnique, 1999, 49(3): 387-403

[2] ZhaoL-h, LiL, YangF, LuoQ, LiuX. Upper bound analysis of slope stability with nonlinear failure criterion based on strength reduction technique [J]. Journal of Central South University of Technology, 2010, 17(4): 836-844

[3] KrabbenhoftK, LyaminA V. Strength reduction finite-element limit analysis [J]. Géotechnique Letters, 2015, 5(4): 250-253

[4] DengD-p, LiL, ZhaoL-H. Limit equilibrium analysis for stability of soil nailed slope and optimum design of soil nailing parameters [J]. Journal of Central South University, 2017, 24(11): 2496-2503

[5] LiuY, ZhangW-g, ZhangL, ZhuZ-r, HuJ, WeiH. Probabilistic stability analyses of undrained slopes by 3D random fields and finite element methods [J]. Geoscience Frontiers, 2018, 9(6): 1657-1664

[6] ZhangR, LongM-x, LanT, ZhengJ-l, GeoffC. Stability analysis method of geogrid reinforced expansive soil slopes and its engineering application [J]. Journal of Central South University, 2020, 27(7): 1965-1980

[7] ChenX, WuY-k, YuY-z, LiuJ-k, XuX F, renJ. A two-grid search scheme for large-scale 3-D finite element analyses of slope stability [J]. Computers and Geotechnics, 2014, 62: 203-215

[8] TschuchniggF, SchweigerH F, sloanS W. Slope stability analysis by means of finite element limit analysis and finite element strength reduction techniques. Part II: Back analyses of a case history [J]. Computers and Geotechnics, 2015, 70: 178-189

[9] WangD-y, ChenX, YuY-z, JieY-x, LyuY-nan. Stability and deformation analysis for geotechnical problems with nonassociated plasticity based on second-order cone programming [J]. International Journal of Geomechanics, 2019, 19(2): 04018190

[10] MühlhausH B, vardoulakisI. The thickness of shear bands in granular materials [J]. Géotechnique, 1987, 373271-283

[11] TejchmanJ, WuW. Numerical study on patterning of shear bands in a Cosserat continuum [J]. Acta Mechanica, 1993, 99(1–4): 61-74

[12] TejchmanJ, WuW. Numerical simulation of shear band formation with a hypoplastic constitutive model [J]. Computers and Geotechnics, 1996, 18(1): 71-84

[13] GracioJ J. The double effect of grain size on the work hardening behaviour of polycrystalline copper [J]. Scripta Metallurgica et Materialia, 1994, 31(4): 487-489

[14] SharbatiE, NaghdabadiR. Computational aspects of the Cosserat finite element analysis of localization phenomena [J]. Computational Materials Science, 2006, 38(2): 303-315

[15] VoyiadjisG Z, alsalehM I, alshibliK A. Evolving internal length scales in plastic strain localization for granular materials [J]. International Journal of Plasticity, 2005, 21(10): 2000-2024

[16] TangH-x, SunF-b, ZhangY-p, DongY. Elastoplastic axisymmetric Cosserat continua and modelling of strain localization [J]. Computers and Geotechnics, 2018, 101: 159-167

[17] WangD-y, ChenX, LvY-n, TangC. Geotechnical localization analysis based on Cosserat continuum theory and second-order cone programming optimized finite element method [J]. Computers and Geotechnics, 2019, 114103118

[18] ChangJ-f, ChuX-h, XuY-jie. Finite-element analysis of failure in transversely isotropic geomaterials [J]. International Journal of Geomechanics, 2015, 15(6): 04014096

[19] ChuJ, LoS C R, LeeI K. Strain-softening behavior of granular soil in strain-path testing [J]. Journal of Geotechnical Engineering, 1992, 1182191-208

[20] ReadH E, hegemierG A. Strain softening of rock, soil and concrete: A review article [J]. Mechanics of Materials, 1984, 3(4): 271-294

[21] ConteE, SilvestriF, TronconeA. Stability analysis of slopes in soils with strain-softening behaviour [J]. Computers and Geotechnics, 2010, 37(5): 710-722

[22] ZhangK, CaoP, BaoR. Progressive failure analysis of slope with strain-softening behaviour based on strength reduction method [J]. Journal of Zhejiang University SCIENCE A, 2013, 14(2): 101-109

[23] RafieiR, MartinC D. Factor of safety of strain-softening slopes [J]. Journal of Rock Mechanics and Geotechnical Engineering, 2020, 12(3): 473-483

[24] LiX-k, TangH-xiang. A consistent return mapping algorithm for pressure-dependent elastoplastic Cosserat continua and modelling of strain localisation [J]. Computers & Structures, 2005, 83(1): 1-10

[25] ReissnerE. On a variational theorem in elasticity [J]. Journal of Mathematics and Physics, 1950, 29(1–4): 90-95

[26] KrabbenhoftK, KarimM R, lyaminA V, sloanS W. Associated computational plasticity schemes for nonassociated frictional materials [J]. International Journal for Numerical Methods in Engineering, 2012, 90(9): 1089-1117

[27] ZhangX, ShengD-c, SloanS W, bleyerJ. Lagrangian modelling of large deformation induced by progressive failure of sensitive clays with elastoviscoplasticity [J]. International Journal for Numerical Methods in Engineering, 2017, 112(8): 963-989

[28] MOSEK, A. The MOSEK C optimizer for API C8.0.0.94 [OB/OL]. [2018-12-25] http://docs.mosek.com/8.0/capi/index.html.

[29] DAVIS, E H. Theories of plasticity and failure of soil masses [M]// Soil Mechanics: Selected Topics, 1968: 341–354.

[30] ChenX, WangD-y, YuY-z, LvY-nan. A modified Davis approach for geotechnical stability analysis involving non-associated soil plasticity [J]. Géotechnique, 2020, 70(12): 1109-1119