In this paper, the authors develop Reissner's mixed variational theorem (RMVT)-based finite layer methods for the three-dimensional (3D) coupled thermoelastic analysis of simply supported, functionally graded, doubly curved (DC) shells with temperature-independent material properties. A two-phase composite material is considered to form the shell, and its material properties are assumed to obey a power-law distribution of the volume fractions of the constituents through the thickness direction of the shell. The effective material properties are estimated using the Mori-Tanaka scheme. The accuracy and convergence rate of these RMVT-based finite layer methods are validated by comparing their solutions with the quasi 3D and accurate two-dimensional solutions available in the literature.