The linear and nonlinear isogeometric finite element models of an axially functionally graded graphene platelet-reinforced composite (AFG-GPLRC) curved beam are established within the framework of the third-order shear deformation beam theory (TSDT) and von-Kármán’s nonlinear geometric relation. The AFG-GPLRC curved beams can be seen as composite structures in which the graphene platelets (GPLs) are continuously distributed in the matrix along the length direction of the curved beam according to different patterns. The modified Halpin-Tsai parallel model and the rule of mixture are implemented to predict the effective Young’s modulus and mass density as well as Poisson’s ratio, respectively. Hamilton's principle, TSDT, and von-Kármán’s strain-displacement relation are combined to derive the governing partial differential equation of motion and corresponding boundary conditions. Furthermore, the Non-Uniform Rational B-splines (NURBS)-based isogeometric analysis (IGA) approach together with a direct iterative technique are utilized to solve the nonlinear governing equation. The accuracy and efficiency of the proposed IGA framework are confirmed by comparing corresponding numerical solutions with other available results. The parametric investigations, such as the curved beam’s geometric parameters, boundary conditions, and GPL’s distribution patterns, on the nonlinear bending and vibration responses of the AFG-GPLRC curved beams are carried out by some illustrative examples.