A discrete element analysis for general failure behavior of basalt

Discrete element method is a commonly used numerical method in rock mechanics and rock engineering. In this research, the flat-joint contact bond model is used to reproduce the macro-properties of basalt on an assembly of particles. A series of laboratory tests, including uniaxial compression test, Brazilian tension test, and fracture toughness test are used to simulate progressive degradation of rock mass strength in a size-independent model with a unique set of micro-properties. The work involves a detailed analysis for estimating compressive elastic modulus, Poisson's ratio, bimodularity, mode I fracture toughness, crack closure, crack initiation, and crack damage thresholds, as well as the tensile and compressive strengths. The constructed flat-jointed material yields a realistic uniaxial compression to tensile strength ratio and nonlinear stress-strain response through sustaining adequate particle interlocking and embedded pre-existing microcracks. The role of sample size on the compressive elastic modulus, crack initiation and crack damage stress limits is demonstrated in detail. A model size-independency analysis was also performed to examine the effect of particle size on the results. The numerical results presented in this study are within 10% of the laboratory results.