Numerical modeling of the failure process of the heterogeneous karst rock mass using the DDA-SPH method

The sanding process caused by karstification in dolomite creates a special sandy dolomite stratum,where the frequent catastrophic instability of the surrounding rock occurred during tunnel construction.In this study,the micro-origin and macro-performance of the sandy dolomite stratum are first discussed.Then,a numerical model based on the coupling method between the discontinuous deformation analysis and smoothed particle hydrodynamics is proposed to depict the heterogeneous dolomite formation with different sanding degrees.Following,the mechanical behaviors of the heterogeneous dolomite samples under uniaxial compression are studied after calibrating the numerical parameters with the two single materials sampled from the tunnel site respectively.Further,the instability disasters of the dolomite surrounding rock with different sanding degrees are reproduced,and the failure behaviors of tunnels are explained with respect to the stress distribution and plastic zone.The obtained results show that the rotation and dislocation of the remained dolomite block contribute to the unsmooth stress–strain curve and deterioration in uniaxial compressive strength.However,the block serves as the skeleton in the transmission of field stress in underground space,which improves the stability of the formation.