The paper deals with a finite element solution of transient thermoelasticity
problems. For each time step the system of linear algebraic equations is solved using the conjugate gradient method
preconditioned by incomplete factorization of the matrix derived from the original matrix. The time step is chosen adaptively. The results of numerical tests are presented. A procedure for the solution of large practical problems is proposed.