A phase-field crack model is developed for numerical analysis of thermal buckling and postbuckling behavior of a functionally graded (FG) graphene platelet-reinforced composite (FG-GPLRC) plate with a central crack. The inclined central crack is represented according to a stable, effective phase-field formulation (PFF) by introducing a virtual crack rotation. The problem is formulated using first-order shear deformation theory (SDT) incorporated with von Kármán geometric nonlinearity. And it is approximated by combining regular Laplace interpolation functions and crack-tip singular functions in the framework of the 2D extended natural element method (XNEM). Troublesome shear locking is suppressed by applying the concept of the MITC (mixed-interpolated tensorial components)3+ shell element to the present numerical method. The results demonstrate the effectiveness of this method in accurately predicting the critical buckling temperature rise (CBTR) and the thermal postbuckling path. In addition, the parametric results reveal that the CBTR and postbuckling path of the FG-GPLRC plate are distinct from those of the FG carbon nanotube-reinforced composite (FG-CNTRC) plate and remarkably affected by the parameters associated with the crack and graphene platelet (GPL).