In order to overcome the inefficiency of existing numerical calculation programs in calculating the contact slip law of complex fracture surfaces

a new numerical calculation method was proposed by absorbing the idea of boundary element method. In this algorithm

all grids were located at the boundary of rock blocks

and the interaction force between rock blocks and the displacement of rock blocks were calculated by using the explicit difference method. The algorithm was validated by two simulation experiments

i. e. the ball sliding on the parabolic surface and the dumbbell sliding on the inclined surface with different inclinations. The results show that: the algorithm can accurately describe the contact between objects and accurately calculate the normal displacement between blocks. The algorithm can accurately judge the "sliding state" and "stable state" of the object. The error between the numerical and analytical results of friction force is less than 10?10. When calculating the variation of bearing capacity during slip of fracture surface

the calculation results of the new method are consistent with those of the finite element method

but the calculation efficiency of the new method is significantly improved.