A novel single-step method is proposed for the analysis of dynamic response of visco-elastic structures containing non-smooth contactable interfaces. In the method, a two-level algorithm is employed for dealing with a nonlinear boundary condition caused by the dynamic contact of interfaces. At the first level, an explicit method is adopted to calculate nodal displacements of global viscoelastic system without considering the effect of dynamic contact of interfaces and at the second level, by introducing contact conditions of interfaces, a group of equations of lower order is derived to calculate dynamic contact normal and shear forces on the interfaces. The method is convenient and efficient for the analysis of problems of dynamic contact. The accuracy of the method is of the second order and the numerical stability condition is wider than that of other explicit methods.
