The paper proposes a new iterative approach for the numerical solution of partially coupled problems of thermo elasticity for isotropic bodies. The tasks are considered in two settings, dynamic and static. A discrete analogue of the boundary value problem is compiled on the basis of the finite difference method and an iterative process is performed, which allows you to find the values of the desired functions. It is assumed that, at the zero approximation, the values of the desired functions in the internal nodes are trivial. The essence of the proposed algorithm is demonstrated by numerically solving the one-dimensional dynamic and two-dimensional static thermo elasticity problems. The proposed algorithm can be applied for the numerical solution of thermoplastic coupled and unbound problems for isotropic and anisotropic bodies.