Abstract
This article presents a virtual boundary meshless with Trefftz method for calculation of the two-dimensional steady-state heat conduction crack problem, which is based on the idea of multidomain combination. The virtual source function is constructed by the radial basis functions approximation in the subdomain not containing the crack. The temperature or heat flux is calculated by crack special solution of Trefftz function in the crack subdomain. Thus, the proposed method has the advantages of both the boundary-type meshless method and the Trefftz method. No nodes or elements are distributed on the crack, because the boundary condition of the crack is automatically satisfied by employing the crack special solution of the Trefftz function. The heat stress intensity factor is simply and directly calculated.