Numerical Study of Dynamic Mechanical Behaviour of Rocks with a Flaw Under Different Infilling Conditions Using a Strain Rate Enhanced Continuum Damage Model
摘要
In this study, the dynamic mechanical behaviour of flawed rock specimens featuring various flaw inclinations and infill materials (cement paste and epoxy resin) is numerically investigated under different strain rates. An existing coupled damage-plastic constitutive model is improved to capture the rate-dependent response of rocks. This enhanced model is incorporated into the finite element software (ABAQUS) via a user-defined material subroutine (UMAT) and subsequently calibrated against experimental data. The experimental study involves conducting dynamic compression tests on synthetic rock specimens using in-house testing facilities. Finite element simulations are then performed on synthetic rock specimens with various flaw infill conditions to evaluate the stress–strain response and track the progression of localised damage. The model captures the response of rock under impact loading across different strain rates, showing good agreement with experimental results. For flawed specimens (unfilled, cement-filled, and epoxy-filled), it accurately predicts changes in stress–strain behaviour with flaw orientation, including peak stress, stiffness, and post-peak response. The simulation results show the spatial evolution of stress concentration with varying flaw orientation. As the flaw angle increases from 0° to 90°, the stress concentration shifts from the flaw tips to regions away from the tips, leading to a steady rise in peak stress, consistent with experimental observations.