This article presents the nonlinear investigation of the thermal and mechanical buckling of orthotropic annular/circular single-layer/bilayer nanoplate with the Pasternak and Winkler elastic foundations based on the nonlocal strain gradient theory. The stability equations of the graphene plate are derived using higher-order shear deformation theory (HSDT) and first-order shear deformation theory (FSDT) considering nonlinear von Karman strains. Furthermore, this paper analyses the nonlinear thermal and mechanical buckling of the orthotropic bilayer annular/circular nanoplate. HSDT provides an appropriate distribution for shear stress in the thickness direction, removes the limitation of the FSDT, and provides proper precision without using a shear correction coefficient. To solve the stability equations, the differential quadratic method (DQM) is employed. Additionally, for validation, the results are checked with available papers. The effects of strain gradient coefficient, nonlocal parameter, boundary conditions, elastic foundations, and geometric dimensions are studied on the results of the nondimensional buckling loads. Finally, an equation is proposed in which the thermal buckling results can be obtained from mechanical results (or vice versa).