Numerical Simulation of Crack Propagation and Branching Behaviors in Heterogeneous Rock-like Materials

Abstract

The characterization and understanding of crack evolution in non-uniform geological structures are crucial for predicting the mechanical response of rock-like materials or structures under varying loading conditions. In this study, an improved Peridynamic model with a degree of heterogeneity characterized by random pre-breaking “bonds” coefficients is introduced to capture the intricacies of crack initiation, propagation, and branching behaviors in heterogeneous rock-like materials. MATLAB discrete programs for heterogeneous material models and PD simulation programs based on the FORTRAN language were developed. The effectiveness of the heterogeneous PD model in simulating crack propagation and branching patterns in heterogeneous materials has been verified through dynamic and static (quasi-static) loading cases with pre-notch. The different levels of heterogeneity not only affect the direction of crack propagation but also determine the crack deflection direction and branching patterns. The crack propagation path appears to possess obvious asymmetry in the crack propagation direction. As the load applied continues to increase, the asymmetric multi-crack branching phenomenon will occur. The higher the level of heterogeneity, the more complex the behaviors of crack propagation and branching become. This research provides valuable insights into the interplay of material heterogeneity and crack evolution, offering a foundation for improved numerical simulations and contributing to the broader field of geomechanics.