In order to simulate the movement and dynamic response of adjacent blocks in an underground cavern group under seismic loads, a dynamic anti-sliding model, considering the bond strength between the blocks, was established. The effectiveness of the anti-sliding model in comparison to other models was analyzed using a numerical model. In the practical engineering examples, the calculated root mean square values of horizontal displacement in the X direction of two adjacent blocks of type I were 2. 16 cm and 2. 02 cm when the classical anti-sliding model was used, and 2. 00 cm and 1. 94 cm when the modified anti-sliding model was used throughout the earthquake process. The corresponding simulated results of two adjacent blocks of type II were 3. 68 cm and 3. 65 cm when the classical model was used and 3. 55 cm and 3. 54 cm when the modified model was used. Hence, it is demonstrated that adjacent blocks considering bonding characteristics show much better stability, and the proposed model has high application value.