The nonlocal nonlinear buckling of a double layer graphene sheet (DLGS) covered by zinc oxide (ZnO) piezoelectric layers is investigated in this study. The surrounding circumstances of the system are considered as a Pasternak foundation including spring constants and a shear layer. Graphene sheets are subjected to longitudinal magnetic field and biaxial forces. On the other hand, the ZnO piezoelectric layer is subjected to an electric field. Eringen’s nonlocal theory is used for considering small-scale effects. Classical plate theory (CPT) is employed to model the plates. Nonlinear Von-Karman theory, the energy method and Hamilton’s principle are utilized to derive the size dependent governing equations. The known numerical differential quadrature method (DQM) is applied to obtain a nonlocal nonlinear buckling load. The detailed parametric study is conducted focusing on the effects of magnetic field strength, the dimensions of plates, small-scale effects and the intensity of the stiffness matrix on the nonlocal nonlinear buckling load of system. Results indicate that intensifying magnetic field makes the system more stable. Furthermore, increase in thickness of both piezoelectric and graphene layers makes the system stiffer, and consequently the buckling load becomes larger. The results of this study might be useful for the designing and manufacturing of graphene-based structures in micro or nanoelectromechanical systems.