Integration algorithm is a key to ensure the better application of the dynamic constitutive model. As for the bounding surface model based on swell–shrink rules in this paper, an explicit integration algorithm based on sub-stepping and error control scheme has be described in detail. For being adaptable to the constitutive model, several modifications have been made. The modified explicit algorithm includes two procedures: the first procedure is the second-order forward modified Euler scheme with sub-stepping and error control; the second procedure is the stress correction. However, due to the shrink rule of the bounding surface, the second procedure is not needed at the unloading phase. Finally, several numerical simulations have been conducted for assessing the performance of the explicit algorithm, which demonstrates that the corresponding algorithm is stable and accurate, and can simulate the mechanical characteristics of the clay pretty well.