This work is devoted to the formulation and numerical solution of coupled boundary value problems of thermoelasticity with respect to deformations. Two equivalent coupled boundary value problems of thermoelasticity with respect to strain and temperature are proposed. The first consists of six differential equations of thermoelasticity found within the framework of the compatibility conditions of Saint-Venant deformations and the heat influx equation with the corresponding initial and boundary conditions. In the second case, the first three of the six differential equations of thermoelasticity are replaced with three differentiated equations of motion. The validity of the formulated two boundary value problems of thermoelasticity is justified by comparing their numerical, obtained by the variable direction method and recurrence relations, as well as solving a similar related problem regarding displacements.