In this paper we analyze, from the numerical point of view, a dynamic thermoelastic problem. Here, the so-called exact heat conduction model with a delay term is used to obtain the heat evolution. Thus, the thermomechanical problem is written as a coupled system of partial differential equations, and its variational formulation leads to a system written in terms of the velocity and the temperature fields. An existence and uniqueness result is recalled. Then, fully discrete approximations are introduced by using the classical finite element method to approximate the spatial variable and the implicit Euler scheme to discretize the time derivatives. A priori error estimates are proved, from which the linear convergence of the algorithm could be derived under suitable additional regularity conditions. Finally, a two-dimensional numerical example is solved to show the accuracy of the approximation and the decay of the discrete energy.