In the process of blasting rock breaking, the existence of joints and fissures in rock mass affects the blasting crushing effect. Some finite element numerical simulation software can realize the blasting rock breaking simulation of jointed rock mass, while it is difficult to realize the quantitative description of cracks. A numerical model of jointed rock mass was established using LS-DYNA software to analyze the obstruction and reflection effects of joint surfaces on explosion stress waves. The images of blast-induced cracks were binarized, and crack extraction from the binary images was achieved using the Hough transform method. Statistical results of the crack images were subsequently analyzed through the maximum likelihood estimation method. The results show that: The existence of joints seriously hinders the propagation of blasting stress wave, and the peak stress is reduced by about 60%. The rock mass between blast hole and joint is subjected to stronger stress superposition due to reflection, and the crushing effect is more remarkable; In the initial stage of explosion, the occurrence probability of small cracks is large, and at the end, it is dominated by long cracks. In this process, many short cracks expand into long cracks, and the distribution of crack angles is relatively average. The distribution of long cracks with joint angles of 45° and 60° is more uniform, and the extreme value is not obvious. When the joints are 15° and 30°, the crack development in the vertical direction is obvious. The Hough transform method realizes the quantitative result statistics of crack images in the numerical simulation results.

To study the response of mode Ⅰ cracks under impact loading, the dynamic fracture test on the single cleavage triangle (SCT) granite specimens with the lateral opening was carried out using a split Hopkinson pressure bar test system. The Hough transform method was used to quantitatively describe the distributions of the surface crack lengths and angles, and the relationship between them and the absorbed energy of the specimen was analyzed. The three-dimensional (3D) point cloud data of the fracture surface was obtained by using a 3D surface topography detection device, and a fracture surface reconstruction method based on the threshold detection of the fitting surface was proposed, which effectively solved the error problem caused by the appearance of error points in the 3D point cloud data of the fracture surface. The relationship between fracture surface roughness and absorbed energy was also discussed. The results show that the exponential distribution of crack length λ tends to increase, and the distribution of crack angles is more even with the increase in absorbed energy. There are relatively few cracks in the horizontal direction, and the cracks have obvious directionality when the dissipated energy of the specimen is lower, while the bending degree of the cracks is higher and the connectivity is better when the dissipated energy is higher. The fracture surface reconstruction method can complete the fracture surface processing better at 0.03-0.08 of the points in the x and y directions and the threshold of 0.25. The roughness statistical parameters of fracture surfaces A and B showed a decreasing trend with the increase in dissipated energy.