Complexity of geomaterial’s behavior is beyond the capabilities of conventional numerical methods alone for realistically model rock structures. Coupling of numerical methods can make the numerical modeling more realistic. Discontinuous Deformation Analysis (DDA) and Displacement Discontinuous Method (DDM) are hybridized for modeling block displacement and crack propagation mechanism in a blocky rock mass. DDA is used to compute the displacements of the blocks, and DDM is used to predict the crack propagation paths due to the specified boundary conditions. The displacements obtained from DDA are converted into stress and considering Kelvin's solution of the problem the crack propagation mechanism within each block is investigated. Boundary stresses are updated due to crack propagation and new stress boundary conditions in DDA. This cycle continued until crack propagation stopped or a new block formed. Numerical solutions of the experimental rock samples including two random cracks with crack 1 fixed and crack 2 created with different angles and one crack with a slope angle of 30 degrees are compared with the existing experimental and numerical results. This comparison validates the accuracy and effectiveness of the proposed procedure because crack propagation paths predicted are in good agreement with the corresponding experimental results of rock samples.