A stochastic continuum damage model for dynamic fracture analysis of quasi-brittle materials using asynchronous Spacetime Discontinuous Galerkin (aSDG) method

Abstract

The microstructural design has an essential effect on the fracture response of brittle materials. We present a stochastic bulk damage formulation to model dynamic brittle fracture. This model is compared with a similar interfacial model for homogeneous and heterogeneous materials. The damage models are rate-dependent, and the corresponding damage evolution includes delay effects. The evolution equation specifies the rate at which damage tends to its quasi-static limit. The relaxation time of the model introduces an intrinsic length scale for dynamic fracture and addresses the mesh sensitivity problem of earlier damage models with much less computational efforts. The ordinary differential form of the damage equation makes this remedy quite simple and enables capturing the loading rate sensitivity of strain-stress response. A stochastic field is defined for material cohesion and fracture strength to involve microstructure effects in the proposed formulations. The statistical fields are constructed through the Karhunen-Loeve (KL) method.An advanced asynchronous Spacetime Discontinuous Galerkin (aSDG) method is used to discretize the final system of coupled equations. Local and asynchronous solution process, linear complexity of the solution versus the number of elements, local recovery of balance properties, and high spatial and temporal orders of accuracy are some of the main advantages of the aSDG method.Several numerical examples are presented to demonstrate mesh insensitivity of the method and the effect of boundary conditions on dynamic fracture patterns. It is shown that inhomogeneity greatly differentiates fracture patterns from those of a homogeneous rock, including the location of zones with maximum damage. Moreover, as the correlation length of the random field decreases, fracture patterns resemble angled-cracks observed in compressive rock fracture. The final results show that a stochastic bulk damage model produces more realistic results in comparison with a homogenizes model