Abstract: Solid thermal stress caused by temperature has always been one of the important factors affecting the life of solid materials.The boundary element method, using boundary integrals instead of domain integrals, reduces the dimensionality of complex issue and the computational complexity. The non-uniform rational B-splines (NURBS) curve is used to construct the boundary of the calculation model, and the isogeometric idea is used to divide the unit, so as to overcome the geometric error caused by the unit discretization. Because of the existence of thermal stress, the constructed isogeometric boundary integral equation contains the domain integration term. The radial basis function method is used to convert the domain integration to the boundary integration. The advantage of the dimensionality reduction can be fully presented by using isogeometric boundary element method. The model calculation shows that the proposed calculation method is more reliable than the traditional boundary element method.