Rock initiation and propagation simulation under compression-shear loading using continuous-discontinuous cellular automaton method

A continuous-discontinuous cellular automaton method is developed for rock initiation and propagation simulations, in which the level set method, discontinuous enrichment shape functions and discontinuous cellular automaton are combined. No remeshing is needed for crack growth analysis, and all calculations are restricted to cells without an assembled global stiffness matrix. The frictional contact theory is employed to construct the contact model of normal pressure and tangential shear on crack surfaces. A discontinuous cellular automaton updating rule suitable for frictional contact of rock is proposed simultaneously with Newton's iteration method for nonlinear iteration. Besides, a comprehensive fracturing criterion for brittle rock under compression-shear loading is developed. The accuracy and effectivenesss of the proposed method is proved by numerical simulation.