In this paper, the finite element formulation for the thermal analysis of quasi-static, uncoupled, homogeneous, isotropic and linear viscoelastic problems is presented based on the principle of virtual work. The viscoelastic material is assumed to be thermorheologically simple, which is well known material property in a large class of high polymers. The variational formulation and the finite element equation in matrix from are derived. Effective generation and storage of the hereditary stiffness matrices are given in detail especially for the case of the steady state temperature distribution T=T(x). Some numerical examples are given and compared with published results to show the versatility of the derived finite element formulations.