In this work, a new mathematical model of thermoelasticity theory has been considered in the context of a new consideration of heat conduction with fractional order theory. A functionally graded isotropic unbounded medium is considered subjected to a periodically varying heat source in the context of space-time non-local generalization of three-phase-lag thermoelastic model and Green-Naghdi models, in which the thermophysical properties are temperature dependent. The governing equations are expressed in Laplace-Fourier double transform domain and solved in that domain. Then the inversion of the Fourier transform is carried out by using residual calculus, where poles of the integrand are obtained numerically in complex domain by using Laguerre’s method and the inversion of Laplace transform is done numerically using a method based on Fourier series expansion technique. The numerical estimates of the thermal displacement, temperature and thermal stress are obtained for a hypothetical material. Finally, the obtained results are presented graphically to show the effect of non-local fractional parameter on thermal displacement, temperature and thermal stress. A comparison of the results for different theories (three-phase-lag model, GN model II, GN model III) is presented and the effect of non-homogeneity is also shown. The results, corresponding to the cases, when the material properties are temperature independent, agree with the results of the existing literature.