Mesh Scalability Concept for Explicit Simulation of Rock Failure

Authors: Dyskin, AV; Caballero, A Download Paper

DOI https://doi.org/10.36487/ACG_repo/808_45

Cite As:
Dyskin, AV & Caballero, A 2008, 'Mesh Scalability Concept for Explicit Simulation of Rock Failure', in Y Potvin, J Carter, A Dyskin & R Jeffrey (eds), SHIRMS 2008: Proceedings of the First Southern Hemisphere International Rock Mechanics Symposium, Australian Centre for Geomechanics, Perth, pp. 185-194, https://doi.org/10.36487/ACG_repo/808_45



Abstract:
Rock failure is a phenomenon which involves the initiation and propagation of (multiple) fractures and the development of surfaces of strain localisation. Direct numerical simulation of these phenomena by introducing the relevant constitutive behaviour in finite elements (or similar computational elements, e.g. particles in the discrete element method) lead to mesh dependence: the computations with successively refined meshes do not stabilise. The mesh dependence is caused by the fact that theoretically, fractures and strain localisation surfaces have zero thickness, while in numerical calculations this thickness is finite controlled by the size of the finite element or particle, or by a characteristic size of the non-local functions in mesh-free methods. Two methods are currently used to overcome this problem. The first is the introduction of a generalised continuum which contributes its own characteristic lengths. The computations stabilise as soon as the size of the computational element becomes considerably smaller than the smallest characteristic length of the continuum. The use of generalised continua as numerical stabilisers leads however to hard questions as to what type of the continuum is to be chosen and how to determine its numerous parameters. In the second method a size of the mesh element is chosen which corresponds to the size of microstructural elements (if this size in known). In this approach the flexibility of choosing computationally efficient element size is lost if the microstructural elements are very small compared with the characteristic size of the geometry to be modelled (e.g. the excavation size) the modelling may become computationally prohibitive, exactly as in the case of particle methods. To deal with mesh dependence we propose the concept of mesh scalability. We use the assumption that, if the element size is sufficiently small, all mesh-dependent quantities should depend upon the size as power functions. The mesh with this property is called scalable and is regarded as sufficiently fine for the simulations. The obtained exponent and the pre-factor can be used to extrapolate the mesh dependent quantities to the actual microstructural sizes. Furthermore, the method permits the identification of stress singularities from finite element computations with different element sizes in similar meshes. The pre-factors are then called the mesh stress intensity factors and can be used to analyse failure instead of the ones referred to in the Linear Fracture Mechanics. We consider examples of application of this philosophy to modelling the stability of rock slopes.

References: