[Objective] Among numerical simulation methods for blasting fragmentation, the continuous medium simulation method has high efficiency, but its mechanical mechanisms are not rigorous and errors are significant when dealing with discontinuous problems; the discontinuous deformation analysis (DDA) method performs well for discontinuous problems, but when the fragment size becomes too small, excessively long computation time and non-convergence are likely to occur. This study aims to propose a numerical simulation method for blasting fragmentation that considers both computational efficiency and mechanical rationality. [Methods] Field tests were conducted to reveal the formation characteristics of blasting fragmentation, and the necessity of selecting appropriate numerical simulation methods for different fragment size ranges of blasting fragmentation was clarified. A continuous-discontinuous numerical simulation method for blasting fragmentation based on LS-DYNA+DDA coupling was proposed. In the near-field region of the blast hole, a continuous medium numerical simulation based on LS-DYNA was used to improve the computational efficiency of the crushing zone. In the middle- and far-field regions, a discontinuous method based on DDA was used to achieve discontinuous characterization of blasting fragmentation. The accuracy of using stress and velocity components as coupling parameters was compared. Finally, the LS-DYNA+DDA coupling method was validated based on the mining and blasting practice of Zhoushan Green Petrochemical Mine. [Results] Through field experiments and numerical simulation, it was determined that small-sized fragments were mainly concentrated within a very small range near the blast hole. The continuous medium method could efficiently simulate the distribution of small-sized fragments while ensuring accuracy. It was more reasonable to use DDA method to simulate the fragmentation of medium- and large-sized fragments. Using peak velocity as the coupling parameter between different methods could reduce the pressure loss during computation transmission. [Conclusion] Based on the measured results, comparison and validation between existing numerical simulation methods and the proposed LS-DYNA+DDA coupling method show that the proposed method improves the accuracy of blasting fragmentation prediction and has advantages in balancing the mechanical rationality of fragmentation mechanisms and computational efficiency. However, this method is currently applied in limited engineering scenarios, and its prediction efficiency needs further summary and optimization for different lithologies and blasting parameters.