▎ 摘　　要

This paper deals with a numerical investigation on the bi-nonlinear response of functionally graded porous (FGP) plates reinforced by graphene platelets (GPLs) with surface bonded piezoelectric layers. Both the porosity distribution type and GPL dispersion pattern of the GPLs reinforced FGP plates are assumed to have three forms: non-uniformly symmetric, non-uniformly asymmetric and uniform. The effective material properties of the GPLs reinforced FGP plates are estimated by using the closed-cell Gaussian random field (GRF) scheme, the Halpin-Tsai micromechanics model and the rule of mixture. The electroelastic nonlinear constitutive equations and von Karman type nonlinear strain-displacement relations are applied to consider both material nonlinearity and geometric nonlinearity. Subsequently, the isogeometric analysis (IGA) method in combination with Reddy's higher-order shear deformation theory (HSDT) is adopted to derive the bi-nonlinear isogeometric finite element motion equations of the structures based on the Hamilton's principle and the Total Lagrangian (TL) incremental formulation. Through the comparison studies, the accuracy and effectiveness of present method are verified, and the comprehensive assessment for bi-nonlinear natural frequency and static bending response of the plates are conducted to prove the significance of this work.