▎ 摘　　要

The present study investigated the nonlinear harmonic vibration of functionally graded multilayer graphene nanoplatelet (GPL)-reinforced nanocomposite beams on the basis of the third-order shear deformation theory. The GPL volume fraction shows a layer-wise change, while in each individual layer GPLs are uniformly dispersed in the matrix. The effective Young's moduli of the GPL-reinforced nanocomposite (GPLRC) beams were estimated through the Halpin-Tsai micromechanics model. The mass densities as well as the effective Poisson's ratios of the GPLRC beams were predicted by the rule of mixture. The nonlinear partial differential equations of motion were discretized by means of the Galerkin procedure. A parametric study was carried out by using the multiple scales method to examine the effects of GPL distribution pattern, weight fraction, geometry, and size on the nonlinear response of the primary, secondary, and combination resonances. Results show that an addition of a very low weight fraction of GPL nanofillers significantly reduces the primary, superharmonic, subharmonic, and combinational resonant responses of the beams. The square-shaped GPLs with fewer graphene layers are the most favorable reinforcements.