We analyze a method for the approximation of exact controls of a second order infinite dimensional system with bounded input operator. The algorithm combines Russell’s “stabilizability implies controllability” principle and a finite elements method of order θ with vanishing numerical viscosity. We show that the algorithm is convergent for any initial data in the energy space and that the error is of order θ for sufficiently smooth initial data. Both results are consequences of the uniform exponential decay of the discrete solutions guaranteed by the added viscosity and improve previous estimates obtained in the literature. Several numerical examples for the wave and the beam equations are presented to illustrate the method analyzed in this article.