A bonded-particle model for rock

A numerical model for rock is proposed in which the rock is represented by a dense packing of non-uniform-sized circular or spherical particles that are bonded together at their contact points and whose mechanical behavior is simulated by the distinct-element method using the two- and three-dimensional discontinuum programs PFC2D and PFC3D. The microproperties consist of stiffness and strength parameters for the particles and the bonds. Damage is represented explicitly as broken bonds, which form and coalesce into macroscopic fractures when load is applied. The model reproduces many features of rock behavior, including elasticity, fracturing, acoustic emission, damage accumulation producing material anisotropy, hysteresis, dilation, post-peak softening and strength increase with confinement. These behaviors are emergent properties of the model that arise from a relatively simple set of microproperties. A material-genesis procedure and microproperties to represent Lac du Bonnet granite are presented. The behavior of this model is described for two- and three-dimensional biaxial, triaxial and Brazilian tests and for two-dimensional tunnel simulations in which breakout notches form in the region of maximum compressive stress. The sensitivity of the results to microproperties, including particle size, is investigated. Particle size is not a free parameter that only controls resolution; instead, it affects the fracture toughness and thereby influences damage processes (such as notch formation) in which damage localizes at macrofracture tips experiencing extensile loading.