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　Karman　nonlinear　displacement-strain　relationship,　the　nonlinear　governing　equations　are　obtained　by　Hamilton　principle　and　dispersed　by　differential　quadrature　(DQ)　method.　The　Newmark-beta　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.