(1. State Key Laboratory of Intelligent Construction and Healthy Operation and Maintenance of Deep Underground Engineering，China University of Mining and Technology，Xuzhou，Jiangsu 221116，China；2. State Key Laboratory of Coal Mine Safety Technology，China Coal Technology and Engineering Group Shenyang Research Institute，Fushun，Liaoning 113122，China；
3. Research Institute of Exploration and Development，Tarim Oilfield Company，PetroChina，Korla，Xinjiang 841000，China)
Abstract：In order to reflect the compressive shear failure mode of rock in the process of compression，and to solve the problem that the original multi-scale rock damage mechanics constitutive model cannot describe the shear weakening，a new nonlinear mesoscopic damage-plastic coupling constitutive model is established in the framework of multi-scale rock damage mechanics by introducing local Mohr-Coulomb plastic yield function. The analytical solution of the constitutive model under conventional triaxial compression is derived，and the calibration method of all parameters of the constitutive model is given. The optimized semi-implicit decoupling algorithm is used to write the constitutive model into LS-DYNA material library，and the convergence and calculation speed of the algorithm are verified by unit test. The triaxial compressive mechanical properties of limestone are simulated by using the proposed meso-mechanical damage model. The numerical simulation results are in good agreement with the experimental data. The influences of the constitutive model parameters on the peak stress，peak strain and brittle-ductile transition of rock are obtained by the sensitivity analysis of constitutive model parameters. The established constitutive model does not need to introduce additional hardening/softening functions，significantly reduces the number of model parameters，and the parameters have a strict calibration method. The correlation flow law can be used to accurately describe the axial，transverse and shear nonlinear mechanical behavior. The convergence of the optimized semi-implicit decoupling algorithm is good，and the calculation process of the compiled LS DYNA solver is stable，and the calculation speed is significantly faster than that of the traditional phenomenological model of rock，which is of great significance for the practical application of the constitutive model.
[1] ZHU H H，YAN J X，LIANG W H. Challenges and development prospects of ultra-long and ultra-deep mountain tunnels[J]. Engineering，2019，5(3)：384–392.
[2] WALTON G，HEDAYAT A，KIM E，et al. Post-yield strength and dilatancy evolution across the brittle-ductile transition in indiana limestone[J]. Rock Mechanics and Rock Engineering，2017，50(7)：1 691–1 710.
[3] XIAO Y M，QIAO Y F，HE M C，et al. A unified strain-hardening and strain-softening elastoplastic constitutive model for intact rocks[J]. Computers and Geotechnics，2022，148：104772.
[4] HORII H，NEMAT-NASSER S. Compression-induced microcrack growth in brittle solids：Axial splitting and shear failure[J]. Journal of Geophysical Research：Solid Earth，1985，90(B4)：3 105–3 125.
[5] ZHANG J C. Experimental and modelling investigations of the coupled elastoplastic damage of a quasi-brittle rock[J]. Rock Mechanics and Rock Engineering，2018，51(2)：465–478.
[6] ZHU Q Z，KONDO D，SHAO J F. Micromechanical analysis of coupling between anisotropic damage and friction in quasi brittle materials：Role of the homogenization scheme[J]. International Journal of Solids and Structures，2008，45(5)：1 385–1 405.
[7] CHEN L，WANG C P，LIU J F，et al. Damage and plastic deformation modeling of Beishan granite under compressive stress conditions[J]. Rock Mechanics and Rock Engineering，2015，48(4)：1 623–1 633.
[8] QIU S L，FENG X T，XIAO J Q，et al. An experimental study on the pre-peak unloading damage evolution of marble[J]. Rock Mechanics and Rock Engineering，2014，47(2)：401–419.
[9] ZHOU H，ZHANG K，FENG X. A coupled elasto-plastic-damage mechanical model for marble[J]. Science China-Technological Sciences，2011，54：228–234.
[10] ZHOU H，BIAN H B，JIA Y，et al. Elastoplastic damage modeling the mechanical behavior of rock-like materials considering confining pressure dependency[J]. Mechanics Research Communications，2013，53：1–8.
[11] HU D W，ZHU Q Z，ZHOU H，et al. A discrete approach for anisotropic plasticity and damage in semi-brittle rocks[J]. Computers and Geotechnics，2010，37(5)：658–666.
[12] PARISIO F，SAMAT S，LALOUI L. Constitutive analysis of shale：A coupled damage plasticity approach[J]. International Journal of Solids and Structures，2015，75–76：88–98.
[13] SHAO J F，JIA Y，KONDO D，et al. A coupled elastoplastic damage model for semi-brittle materials and extension to unsaturated conditions[J]. Mechanics of Materials，2006，38(3)：218–232.
[14] ZHU Q，SHAO J. Micromechanics of rock damage：Advances in the quasi-brittle field[J]. Journal of Rock Mechanics and Geotechnical Engineering，2017，9(1)：29–40.
[15] ZHU Q Z，SHAO J F，KONDO D. A micromechanics-based thermodynamic formulation of isotropic damage with unilateral and friction effects[J]. European Journal of Mechanics A-Solids，2011，30(3)：316–325.
[16] BASISTA M，GROSS D. The sliding crack model of brittle deformation：An internal variable approach[J]. International Journal of Solids and Structures，1998，35(5–6)：487–509.
[17] 王体富，周宗红. I–II型复合裂纹对脆性岩石受压破坏的影响[J]. 有色金属：矿山部分，2022，74(4)：88–94.(WANG Tifu，ZHOU Zonghong. Influence of mode I–II compound crack on the failure of brittle rock under compression[J]. Nonferrous Metals：Mine Section，2022，74(4)：88–94.(in Chinese))
[18] ZHAO L，ZHU Q，XU W，et al. A unified micromechanics-based damage model for instantaneous and time-dependent behaviors of brittle rocks[J]. International Journal of Rock Mechanics and Mining Sciences，2016，84：187–196.
[19] SHEN W Q，SHAO J F，KONDO D，et al. A micro–macro model for clayey rocks with a plastic compressible porous matrix[J]. International Journal of Plasticity，2012，36：64–85.
[20] 赵伦洋，赖远明，牛富俊，等. 基于增量变分理论的多孔岩石多尺度弹塑性本构模型[J]. 岩石力学与工程学报，2022，41(增2)：3 163–3 173.(ZHAO Lunyang，LAI Yuanming，NIU Fujun，et al. A multi-scale elastoplastic constitutive model of porous rock based on incremental variational theory[J]. Chinese Journal of Rock Mechanics and Engineering，2022，41(Supp.2)：3 163–3 173.(in Chinese))
[21] 朱其志，刘海旭，王 伟，等. 北山花岗岩细观损伤力学本构模型研究[J]. 岩石力学与工程学报，2015，34(3)：433–439.(ZHU Qizhi，LIU Haixu，WANG Wei，et al. A micromechanical constitutive damage model for Beishan granite[J]. Chinese Journal of Rock Mechanics and Engineering，2015，34(3)：433–439.(in Chinese))
[22] ZHU Q Z，SHAO J F. A refined micromechanical damage-friction model with strength prediction for rock-like materials under compression[J]. International Journal of Solids and Structures，2015，60–61：75–83.
[23] BUDIANSKY B，O?CONNELL R J. Elastic moduli of a cracked solid[J]. International Journal of Solids and Structures，1976，12(2)：81–97.
[24] ZHU Q Z，SHAO J F. Multiscale modelling of true triaxial behaviors of brittle rocks[J]. Secondary，2015：1 961–1 966.
[25] 谢 发. 基于摩尔库伦准则岩石材料UMAT子程序二次开发[硕士学位论文][D]. 哈尔滨：哈尔滨工业大学，2014.(XIE Fa. Development of umat subroutine for rock baserd on mohr-coulomb criterion[M. S. Thesis][D]. Harbin：Harbin Institute of Technology，2014.(in Chinese))
[26] SLOAN S，BOOKER J. Removal of singularities in tresca and mohr-coulomb yield functions[J]. Communications in Applied Numerical Methods，1986，2(2)：173–179.
[27] 朱其志. 多尺度岩石损伤力学[M]. 北京：科学出版社，2019：150–160.(ZHU Qizhi. Multi-scale rock damage mechanics[M]. Beijing：Science Press，2019：150–160.(in Chinese))
[28] SIMO J C，TAYLOR R L. Consistent tangent operators for rate-independent elastoplasticity[J]. Computer Methods in Applied Mechanics and Engineering，1985，48(1)：101–118.
[29] ZHANG J C. Experimental and modelling investigations of the coupled elastoplastic damage of a quasi-brittle rock[J]. Rock Mechanics and Rock Engineering，2018，51(2)：465–478.
[30] GENG D，DAI N，GUO P，et al. Implicit numerical integration of highly nonlinear plasticity models[J]. Computers and Geotechnics，2021，132：1–10.
[31] FU B，FENG D C，JIANG H. A new family of explicit model-based integration algorithms for structural dynamic analysis[J]. International Journal of Structural Stability and Dynamics，2019，19(6)：1950053.
[32] LI J，YU K，LI X. An identical second-order single step explicit integration algorithm with dissipation control for structural dynamics[J]. International Journal for Numerical Methods in Engineering，2021，122(4)：1 089–1 132.
[33] TU X，ANDRADE J E，CHEN Q. Return mapping for nonsmooth and multiscale elastoplasticity[J]. Computer Methods in Applied Mechanics and Engineering，2009，198(30)：2 286–2 296.
[34] ZHU Q Z，ZHAO L Y，SHAO J F. Analytical and numerical analysis of frictional damage in quasi brittle materials[J]. Journal of the Mechanics and Physics of Solids，2016，92：137–163.
[35] ZHANG J C. Experimental and modelling investigations of the coupled elastoplastic damage of a quasi-brittle rock[J]. Rock Mechanics and Rock Engineering，2018，51(2)：465–478.
[36] ZHU Q，SHAO J. A micro-mechanics based plastic damage model for quasi-brittle materials under a large range of compressive stress[J]. International Journal of Plasticity，2018，100：156–176.
[37] 谢和平，鞠 杨，黎立云. 基于能量耗散与释放原理的岩石强度与整体破坏准则[J]. 岩石力学与工程学报，2005，24(17)：3 003–3 010. (XIE Heping，JU Yang，LIN Liyun. Criteria for strength and structural failure of rocks based on energy dissipation and energy release principles[J]. Chinese Journal of Rock Mechanics and Engineering. 2005，24(17)：3 003–3 010.(in Chinese))
2024, Vol. 43 
