This paper considers the numerical implementation of two-dimensional thermoviscoelastic waves. The elastic collision of an aluminum cylinder with a two-layer plate of aluminum and iron is considered. In work [1] the difference schemes and algorithm of their realization are given. The most complete reviews of the main methods of calculation of transients in deformable solids can be found in [2, 3, 4], which also indicates the need and importance of generalized studies on the comparative evaluation of different methods and identification of the areas of their most rational application. In the analysis and physical interpretation of numerical results in this work it is also useful to use a priori information about the qualitative behavior of the solution and all kinds of information about the physics of the phenomena under study. Here is the stage of evolution of contact resistance of collision – plate, stress profile.