Abstract:The volume integrals caused by temperature has to be dealt with while the boundary element method is employed to solve thermoelastic mechanics problems,as a result,the advantages of reducing dimensions are extinguished.This paper presents,a numerical approach for solving uncoupled thermoelastic problems with multidomain combinations without calculating volume integral,which is the formation of the numerical formats established for solving heat conduction and corresponding elasticity in accordance with virtual boundary element method (VBEM),respectively,and two ideas are combined,too.The method is available for both 2D and 3D problem,and can also degenerates singledomain problem solving.The method can also be applied to thermoelastic problems of composite materials containing holes or inclusions,arbitrary combinations with different material properties and contact problems.In the end,several numerical examples are given to illustrate the performance of the method,and the results validate the high accuracy and efficiency of the method.