A method of solving three-dimensional coupled quasi-static problems of linear thermoelasticity is presented. Any initial-boundary value problem contains a coupled initial state. The concept is based on extending a region occupied by the considered body over the whole space where the space value thermoelastic problem is formulated. Both the Helmholtz displacement potentials and temperature from the space value Helmholtz problem are convolutionally determined. All the outside component parts of the convolution products are here defined as fictitious ones. The solution of space value thermoelastic problem includes fictitious displacement-temperature components. Capacities of approximate fictitious components are calculated from the boundary condition contracted to the finite time interval. Approximate solution to the primary thermoelasticity problem is obtained by spatial contracting of the solution to the space value thermoelastic problem.