Thermal strength is one of the main thermophysical characteristics of structural materials. For homogeneous systems it is determined by the strength characteristics of the material. While for inhomogeneous systems, in particular, multiphase ones, it is necessary to consider the nature of the microstructure. Heat resistant real materials such as steels are known to be multi-phase systems. One of the mechanisms of their destruction is associated with the presence of propagating heat fluxes that generate thermal stresses. The aim of this paper is to evaluate the patterns of the formation of spatial distributions of thermal stresses in matrix systems of round inclusions characterized by different mutual disposition. The spatial distributions of thermal stresses in a two-phase material characterized by a matrix structure with round inclusions are investigated. For the numerical solution of the problem of stationary thermal conductivity the finite element method with discretization of the medium by triangular elements is used. It was found that at certain points in the medium the values of thermal stresses are ten times higher than the average for the material. It is shown that the spatial distribution and the local magnitude of the temperature gradient depend on the shape of the particles of the phase components and the values of their thermal conductivities. It is considered that the elastic moduli of inclusion and matrix differ little from each other.