Abstract In this paper, one of the newest analytical methods, new homotopy perturbation method (NHPM), is considered to solve thermoelasticity equations. Results obtained by NHPM, which does not need small parameters, are compared with the numerical results and a very good agreement is found. This method provides a convenient way to control the convergence of approximation series and adjust convergence regions when necessary. The results reveal that the proposed method is explicit, effective and easy to use.