In the present research, finite element solutions of thermal post-buckling load-bearing strength of functionally graded (FG) sandwich shell structures are reported by adopting a higher-order shear deformation type kinematics. For the numerical calculation, nine nodes are considered for each element. A specialized MATLAB code is developed incorporating the present mathematical model to evaluate the numerical buckling temperature. The Green-Lagrange nonlinear strain is adopted for the formulation of the sandwich structure. The eigenvalue equation of the FG sandwich structure is solved to predict the postbuckling temperature values of the structure. Moreover, three kinds of temperature distributions across the panel thickness are assumed, viz., uniform, linear and nonlinear. In addition, the properties are described using the power law distributions. The numerical solutions are first validated and, subsequently, the impact of alterations of structural parameters, viz., the curvature ratios, core-face thickness ratios, support conditions and power law index (n Z ) including the amplitude ratio on the thermal post-buckling response of FG sandwich curved panels have been studied in details. The investigation reveals different interesting outcomes, which may help for future references for the analysis and design of the graded sandwich structure.