This paper presents an investigation on thermo-electro-mechanical nonlinear low-velocity impact behaviors of the geometrically imperfect functionally graded (FG) graphene-reinforced composite (GRC) beam with surface-bonded piezoelectric layers. Both uniformly distributed and FG patterns of graphene nanoplatelets (GPLs) are considered along the thickness of the GRC host beam. The effective Young’s modulus is calculated by the Halpin–Tsai model. The Poisson ratio and mass density are calculated by the rule of mixture. The modified nonlinear Hertz contact law is employed to predict the impact force between the spherical impactor and the geometrically imperfect GRC piezoelectric beam during impacting. By considering the first-order shear deformation theory and von Kármán nonlinear displacement–strain relationship, the nonlinear governing equations are obtained by Hamilton principle and dispersed by differential quadrature (DQ) method. The Newmark-β method associated with Newton–Raphson iterative process is adopted to parametrically identify the impact force and the dynamic response of the system. The effects of geometric imperfection, weight fraction and distribution pattern of GPLs, temperature variation, thickness of piezoelectric layer and impactor’s initial velocity on nonlinear low-velocity impact behaviors of geometrically imperfect GRC beams are discussed in detail. Our results illustrate that the coupling effect of geometric imperfection and thermo-electro-mechanical load has a significant effect on the nonlinear low-velocity impact behavior of GRC beam, and GPLs distributing into the piezoelectric layers is better for reducing the impact response of geometrically imperfect GRC piezoelectric beam.