The complexity of geo-materials’ behaviors surpass the capabilities of conventional numerical methods alone, particularly when realistically modeling rock structures. Coupling methods can enhance the realism of numerical modeling. Discontinuous deformation analysis (DDA) and displacement discontinuity method (DDM) are combined to model block displacement and crack propagation mechanisms in a blocky rock mass. DDA calculates stresses, strains, rotations, and displacements of the blocks, while DDM predicts crack propagation paths based on specified boundary conditions. The displacements result from DDA are transformed into stresses, followed by investigating crack propagation mechanisms within each block using Kelvin's solution. Boundary stresses are adjusted as cracks propagate, updating stress boundary conditions in DDA. This iterative process continues until crack propagation halts or a new block emerges. In this paper, the complete failure process of a rock slab with one pre-existing non-persistent joint located on the crest of a rock slope and compression test on the rock-like specimen with a single random crack with a 45° slope angle is simulated using the proposed numerical approach. The proposed numerical solution is compared with existing numerical results. This comparison confirms the accuracy and effectiveness of the proposed procedure, as the predicted crack propagation paths from the coupled numerical method align well with the corresponding numerical results from the rock-like samples.