“…The key to the accuracy of the simulation results lies in the establishment of constitutive equations that reflect the evolution of the rock continuum-discrete. At present, scholars widely base their views on the "continuum assumption", using finite element and finite difference methods, through the ideal elastoplastic model [15]; the strain softening plastic model [16]; the plastic model considering the influence of confining pressure on the subsequent yield function [17]; the elastic-plastic damage model [18] based on the shear strength criteria, such as Mohr-Coulomb and Mises shear strength [19,20]; strain energy [21]; and damage variables [22], as criteria, according to the plastic zone, the damage zone, and the energy zone, determined by the boundary of the water-conducting fracture zone. In addition, based on the assumption of discrete blocks, most scholars use the discrete element method to simplify the coal rock matrix into rectangular blocks or spheres and endow it with elastic or rigid mechanical properties, while the joints between matrices often adopt the Coulomb slip model [23,24], according to Newton's second law, Numerical simulation is an important method to study the fractures from mining overlying rocks.…”