Influence of specimen size, material ductility, and porosity on process zone in mode I fracture of rock

Numerical three-point bending tests were conducted on a softening material model to study the effect of specimen size, material ductility and porosity on the size and shape of the fracture process zone. A bonded particle model (BPM) was used for simulation of the rock. The particles at the contact points were allowed to follow a softening behavior to be able to capture the initiation and development of the process zone. Five different beam sizes of 20 (height) x 60 (span), 40 x 120, 80 x 240, 160 x 480, and 320 x 960 mm2 were used. All beams had a notch at their center to study the mode I fracturing of rock. For each specimen size, six different realizations were introduced to study the effect of particle arrangement on the induced damage zone. The material ductility was controlled and modified by introducing different slopes for the post peak behavior of the contact points between the particles. The shape of the process zone was obtained by calculating both the width and the length of the process zone at the peak load. The results suggest that as the specimen size increases, the process zone expands in its size. In addition, the results indicate that for quasi-brittle materials, the length to width ratio of the process zone is greater compared to that of ductile materials. Furthermore, it is shown that the sizes of the representative elemental volumes corresponding to the width and the length of the process zone may not be identical in the quasi-brittle materials. To change the porosity, some particles were deleted from the domain of analysis. It is shown that both porosity and ductility affect the size of the process zone in mode I fracture. While for more brittle materials, the increase in porosity leads to the greater fracture process zone, this is not the case for more ductile materials. In this case, the process zone can shrink in size with increase in porosity.