Based on the first-order shear deformation (FOSD) hypothesis, a geometrically nonlinear finite element formulation is developed for static and dynamic analysis of carbon nanotube reinforced magneto-electro-elastic (CNT-MEE) plates. The nonlinear finite element model takes the fully geometrically nonlinear strain–displacement relations with large rotations. The nonlinear dynamic governing equations are derived using Hamilton’s principle. The electric and magnetic potentials are assumed to be varied only along the thickness direction of CNT-MEE plates. Four functionally graded (FG) pattens, FG-U, FG-X, FG-O and FG-V are computationally studied for understanding the reinforcement efficiency of the magneto-electro-elastic matrix. The proposed geometrically nonlinear model is verified by comparing with the computational results from the literature. Furthermore, the volume fractions of CNT reinforcement material, distribution forms, reinforcement angles and volume fractions of MEE matrix material are studied by systematically in the computational simulation. These results can be used as a benchmark for designing and analyzing the magneto-electro-elastic properties of CNT-MEE structures with functionally graded reinforcements in extreme environments.