SUMMARYIn this paper, the boundary element method (BEM) for solving quasi-static uncoupled thermoelasticity problems in materials with temperature dependent properties is presented. The domain integral term, in the integral representation of the governing equation, is transformed to an equivalent boundary integral by means of the dual reciprocity method (DRM). The required particular solutions are derived and outlined. The method ensures numerically efficient analysis of thermoelastic deformations in an arbitrary geometry and loading conditions. The validity and the high accuracy of the formulation is demonstrated considering a series of examples. In all numerical tests, calculation results are compared with analytical and/or finite element method (FEM) solutions.