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.

3D simulation of the micro indentation process and relative problems research

In this paper the finite element method was used to simulate micro-scale indentation process. The several standard indenters were simulated with 3D finite element model. The emphasis of this paper was the differences between 2D axisymmetric cone model and

Sedimentation Process of Saturated Sand under Impact Loading

The initial small inhomogeneity of saturated sand could be amplified during the sedimentation process after liquefaction, and cracks could be observed in the sand column. Layers of fine sand could also be found at the exact place where cracks developed and disappeared. The phenomena and the whole process were experimentally shown by X-rays images. To account for the phenomena, a linearized stability analysis of the sedimentation of saturated sand was conducted; however, it did not produce a satisfactory result. A three-phase flow model describing the transportation of fine sand is presented in this paper. It is shown that such a kind of erosion/deposition model was qualitatively in good agreement with the experimental observation.

Molecular dynamics simulation on dislocation emission process

A correlative reference model for a computer simulation of molecular dynamics is proposed in this paper. Based on this model, a flexible displacement boundary scheme is naturally introduced and the dislocations emitted from a crack tip are presumed to continuously pass through the border of an inner discrete atomic region to pile up at an outer continuum region. The simulations for a Mo crystal show that the interaction between a crack and emitted dislocations results in the decrease in local stress intensity factor gradually.