Discrete Element Modeling of the Influence of Void Size and Distribution on the Mechanical Behavior of Rock

The effect of porosity on mechanical behavior of rock has been studied in the past using both experimental and numerical techniques. One important aspect of this study that has not received much attention is the role of void size and distribution on the mechanical behavior of rock. In this study, the CA2 computer program that is a hybrid discrete-finite element program for two-dimensional simulation of geomaterials is used. Numerical specimens with different macro-void sizes are loaded uniaxially. The Young?s modulus, Poisson?s ratio, uniaxial strength and crack initiation stress are studied. The numerical results are compared with the database of experimental tests on Plaster of Paris specimens containing Styrofoam spheres with different void sizes and porosities. It is shown that the macro void spatial distribution is partly responsible for the scatter of the experimental data. It is also shown that not-only the porosity but also the pore size can affect the deformability, crack initiation stress and rock strength.