In this paper, a set of internal stress integral equations is derived for solving thermoelastic problems.A jump term and a strongly singular domain integral associated with the temperature of the material are produced in these equations. The strongly singular domain integral is then regularizedusing a semi-analytical technique. To avoid the requirement of discretizing the domain into internal cells,domain integrals included in both displacement and internal stress integral equations are transformedinto equivalent boundary integrals using the radial integration method (RIM). Two numerical examples for 2D and 3D, respectively, are presented to verify the derived formulations.