Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.