The nonlinear post-buckling response of functionally graded (FG) copper matrix plates enforced by graphene origami auxetic metamaterials (GOAMs) is investigated in the current work.The auxetic material properties of the plate are controlled by graphene content and the degree of origami folding, which are graded across the thickness of the plate. The material properties of the GOAM plate are evaluated using genetic micro-mechanical models. Governing nonlinear eigenvalue problems for the post-buckling response of the GOAM composite plate are derived using the virtual work principle and a four-variable nonlinear shear deformation theory. A novel differential quadrature method (DQM) algorithm is developed to solve the nonlinear eigenvalue problem. Detailed parametric studies are presented to explore the effects of graphene content, folding degree, and GO distribution patterns on the post-buckling responses of GOAM plates. Results show that high tunability in post-buckling characteristics can be achieved by using GOAM. Functionally Graded Graphene Origami Auxetic Metamaterials (FG-GOAM) plates can be used in aerospace structures to improve their structural performance and response