This paper applies the variational iteration method to solve the Cauchy problem arising in one dimensional non-linear thermo-elasticity. The advantage of this method is to overcome the difficulty of calculation of Adomian’s polynomials in the Adomian’s decomposition method. The numerical results of this method are compared with the exact solution of an artificial model to show the efficiency of the method. The approximate solutions show that the variational iteration method is a powerful mathematical tool for solving non-linear problems.