The main aim of this paper was to develop an advanced processing method for analyzing of anisotropic thermoelastic metal and alloy discs with holes. In the boundary element method (BEM), the heat impact is expressed as an additional volume integral in the corresponding boundary integral equation. Any attempt to integrate it directly will necessitate domain discretization, which will eliminate the BEM’s most distinguishing feature of boundary discretization. This additional volume integral can be transformed into the boundary by using branch-cut redefinitions to avoid the use of additional line integrals. The numerical results obtained are presented graphically to show the effects of the transient and steady-state heat conduction on the quasi-static thermal stresses of isotropic, orthotropic, and anisotropic metal and alloy discs with holes. The validity of the proposed technique is examined for one-dimensional sensitivity, and excellent agreement with finite element method and experimental results is obtained.