A hybrid finite-discrete element method for modelling cracking processes in sandy mudstone containing a single edge-flaw under cyclic dynamic loading

Rock mass deformation and failure are macroscopic manifestations of crack initiation, propagation, and coalescence. However, simulating the transition of rocks from continuous to discontinuous media under cyclic dynamic loading remains challenging. This study proposes a hybrid finite-discrete element method (HFDEM) to model crack propagation, incorporating a frequency-dependent cohesive-zone model. The mechanical properties of standard sandy mudstone under quasi-static and cyclic dynamic loading were simulated using HFDEM, and the method's reliability was verified through experimental comparison. The comparative analysis demonstrates that HFDEM successfully captures crack interaction mechanisms and accurately simulates the overall failure behavior of specimens. Additionally, the effects of pre-existing flaw inclination angle and dynamic loading frequency on rock failure mechanisms were investigated. The numerical results reveal that rock samples exhibit significantly higher compressive strength under dynamic loading compared to quasi-static loading, with compressive strength increasing with higher cyclic dynamic load frequencies. Furthermore, by analyzing the strength characteristics, crack propagation, and failure modes of the samples, insights into the failure mechanisms of rocks under different frequency loads were obtained. This study provides valuable insights into crack development and failure of rocks under seismic loads, offering guidance for engineering practices.