Bending and Buckling of FG-GRNC Laminated Plates via Quasi-3D Nonlocal Strain Gradient Theory

: To improve the structural stiffness, strength and reduce the weight of nanoplate structure, functionally graded (FG) graphene-reinforced nanocomposite (GRNC) laminated plates are exploited in this paper. The bending and buckling behaviors of FG-GRNC laminated nanoplates are investigated by using novel quasi-3D hyperbolic higher order shear deformation plate theory in conjunction with modiﬁed continuum nonlocal strain gradient theory, which considered both length and material scale parameters. The modiﬁed model of Halpin–Tsai is employed to calculate the effective Young’s modulus of the GRNC plate along the thickness direction, and Poisson’s ratio and mass density are computed by using the rule of mixture. An analytical approach of the Galerkin method is developed to solve governing equilibrium equations of the GRNC nanoplate and obtain closed-form solutions for bending deﬂection, stress distributions and critical buckling loads. A detailed parametric analysis is carried out to highlight inﬂuences of length scale parameter (nonlocal), material scale parameter (gradient), distribution pattern, the GPL weight fraction, thickness stretching, geometry and size of GPLs, geometry of the plate and the total number of layers on the stresses, deformation and critical buckling loads. Some details are studied exclusively for the ﬁrst time, such as stresses and nonlocality effect.


Introduction
Functionally graded materials (FGMs) are a novel class of composite materials that possess a gradual variation of constituents through spatial directions.FGM is proposed to reduce the abrupt local stress concentrations between different layers, and is used as a thermal barrier in aircraft, submarines, space station structures and fusion reactors.It has many enhanced properties, such as higher fracture toughness, improved stress spreading, enhanced thermal resistance and inferior stress intensity factors.They are now developed for general use in various fields of engineering.
Due to the remarkable physical, thermal and mechanical properties of graphene, as well as better dispersion, excellent reinforcing nanofillers and relatively low manufacturing cost, lots of studies have begun to turn their attention towards graphene as a functionally graded graphene-reinforced nanocomposite (FG-GRNC).Chen et al. [1] exploited Timoshenko beam theory and mid-plane stretching to predict the post-buckling and nonlinear vibration of multilayer FG-GRNC including porosity effect.Shen et al. [2,3] modeled and analyzed the thermal bending and post-buckling of GRNC laminated plates resting on an elastic foundation and subjected to in-plane temperature variation.Huang et al. [4] studied nonlinear buckling analysis of FG-GRNC shallow arches with elastic rotational constraints under uniform radial load.Song et al. [5] presented bending and buckling analyses of multilayer FG-GRNC polymer composite plates within the framework of the first-order shear deformation theory.The effective Young's modulus of the nanocomposites is estimated through the Halpin-Tsai micromechanics model.Garcia-Macias et al. [6] adopted the two-parameter agglomeration model to estimate the agglomeration effects of CNTs as reinforcement in composite structure and presented effects of restacking of graphene sheets on the bending and free-vibration behaviors of composite plates.Dong et al. [7] investigated the buckling behavior of FG-GRNC porous cylindrical shells with spinning motion and subjected to external axial compressive force and radial pressure in frame of first-order shear deformation theory.Daikh and Megueni [8] presented influences of plate aspect ratio, gradient index and the thermal loading conditions on the critical buckling of FGM sandwich plates modeled by the higher-order shear deformation plate theory.Liu et al. [9,10] investigated bending, buckling and vibration behaviors of an initially stressed FG cylindrical shell and circular plate reinforced with nonuniformly distributed graphene platelets (GPLs) using the state-space formulation.Polit et al. [11] presented impacts of different dispersion patterns for the graphene and porosity, shallowness of the curved beam, thickness ratio and platelet geometry on the bending and elastic stability of a higher-order beam with nonlinear stretching.Yang et al. [12] studied nonlinear in-plane buckling of fixed shallow FG-GRNC arches subjected to uniform radial load and temperature field.Mao and Zhang [13] presented influences of electric potential and axial forces on the buckling and post-buckling of FG-GRNC with piezoelectric properties.Anirudh et al. [14] investigated the mechanical response of porous FG-GRNC curved beams by using a trigonometric shear deformation theory and finite element method.Tam et al. [15,16] investigated nonlinear bending, free vibration and buckling characteristics of FG-GRNC beams containing open edge cracks by using the finite element method.Thai et al. [17,18] studied mechanical behaviors of multilayer FG-GRNC plates based on the four-variable refined plate theory and modified couple stress theory.Zhao et al. [19] presented a comprehensive review for FG-GRNC structures including mechanical properties, existing micromechanics models, technical challenges and future research directions.Eltaher et al. [20] exploited finite element method to present the complex phenomena of the elastic and elastoplastic indentation responses of FGMs in the framework of frictional contact mechanics.Rahimi et al. [21] studied bending and vibration behaviors of cylindrical GRNC shell using the state space technique.Shahgholian et al. [22] explored buckling of a porous cylindrical FG-GRNC shell using first-order shear deformation theory.Hamed et al. [23,24] illustrated the effect of in-plane varying compressive force on critical buckling loads and buckling modes of composite laminated and FG porous sandwich beams rested on elastic foundation using unified higher order beam theory.Wang et al. [25] investigated the nonlinear bending of FG-GRNC plate with dielectric permittivity within the framework of first-order shear deformation plate theory.Daikh et al. [26] studied thermal buckling of power and sigmoidal FG sandwich beams based on the higher-order shear deformation beam.Karami and Shahsavari [27] investigated the forced resonant vibration analysis of FG-GRNC doubly curved nano-shells.Four different geometries of the shells namely spherical, elliptical, hyperbolic and cylindrical have been studied, and the Halpin-Tsai model and a rule of mixture are exploited to estimate the effective material properties.
For a nanomechanics analysis, unfortunately the classical approach of continuum mechanics commonly discarded the effects of micro-level and bypassed it altogether by identifying mechanical properties (i.e., Young's modulus, yield stress, ultimate strength, etc.) directly from macroscopic analysis.When dimensions of the structure are comparable to the material internal length scale, the mechanical behavior of these micro-or nanostructures is highly influenced by the material microstructure.The impact of length scales of microstructures is well proved and described by experiments on nanoscale particular geometries.For that reason, to envisage mechanical responses of structure up to micro and nano size accurately, several continuum mechanics theories have been evolved, such as nonlocal elasticity of Eringen [28,29], strain gradient theory [30,31], couple stress theory [32,33] and surface elasticity theory [34,35] that take into account material length scale parameters.
Nowadays, adopted continuous mechanics theories are recommended to consider the missing contribution of micro/nanoscale by many researchers.Sahmani et al. [36] illustrated the size-dependent nonlinear bending of porous FG-GRNC nanobeam and subjected to the uniform distributed load and an axial compressive load via the nonlocal strain gradient theory of elasticity.Emam et al. [37] studied analytically the post-buckling and vibration response of curved multilayer nanobeams subjected to a pre-stress compressive load by using nonlocal elasticity.Ebrahimi and Barati [38] studied hygro-thermal and sizescale effects on the vibration response of graphene sheets rested on viscoelastic medium by employing nonlocal strain gradient theory.Daikh et al. [39] analyzed analytically the thermal buckling of porous FG sandwich nanoplates resting on a Kerr foundation via nonlocal strain gradient theory.Fattahi et al. [40] developed a nonlocal strain gradient beam model to predict the nonlinear secondary resonance of FG porous micro/nanobeams under periodic excitation by using multiple timescales together with the Galerkin technique.Karami et al. [41] studied the static bending of FG anisotropic nanoplates made of hexagonal beryllium crystals using nonlocal strain gradient theory.Jalaei and Civalek [42] applied nonlocal strain gradient theory to examine thermal effects on the dynamic instability of a graphene sheet under periodic axial load.Daikh et al. [43] developed a comprehensive model based on nonlocal strain gradient constitutive relation to study the bending behavior of cross-ply carbon nanotube reinforced composite (CNTRC) laminated nanobeams under various loading profiles.Ebrahimi and Dabbagh [44] studied the viscoelastic wave propagation of axially motivated double-layered graphene sheets via nonlocal strain gradient theory.Eltaher and Mohamed [45] and Eltaher et al. [46] presented analytically and numerically the size-scale effect of CNT beam structure on the buckling stability and free vibration via doublet mechanics theory.Mohamed et al. [47] exploited the energy-equivalent method to study the size scale effect on the buckling and post-buckling of single-walled carbon nanotube rested on nonlinear elastic foundations.Daikh and Zenkour [48] proposed a refined higher order nonlocal strain gradient theory to predict stresses and deflections of FG sandwich nanoplates resting on Pasternak elastic foundation.Based on nonlocal elasticity, Daikh et al. [49] studied analytically the vibration of FG sandwich nanoplates in the thermal environment using higher shear deformation theory.Torabi et al. [50] studied the dynamic and pull-in instability of FG nanoplates via nonlocal strain gradient theory by using the homotopy as an analytical solution methodology.Xiao and Dai [51] studied the static behavior of a circular nanotube made of FG materials by using nonlocal strain gradient theory and a refined shear model.Chaleshtari et al. [52] presented the mutual influence of geometric parameters and mechanical properties on thermal stresses in composite laminated plates with rectangular holes.Dastjerdi et al. [53] developed a semi analytical solution to study the bending behavior of moderately thick FGM plates in a hygro-thermal environment via a quasi-3D approach with nonlocal constitutive equation.Abo-bakr et al. [54] exploited the Pareto optimality method to obtain the optimal weight under the maximum buckling of FG beam under variable axial load.Melaibari et al. [55,56] studied the free vibration behavior of composite laminated shells reinforced by both randomly (CNTRC) and functionally graded fibers.
Generally, for thin structures, the impact of thickness stretching on the mechanical response can be ignored, but due to the special structures such as multilayered composite plates, which are thick or moderately thick, ignoring the effect of thickness stretching in the analysis of mechanical behaviors can lead to inaccurate calculation results and large errors.Moreover, for the proposed structure, the effect of thickness stretching tends to increase as the total number of layers, distribution pattern and the weight fraction of GPLs increase.For this purpose, several researchers investigated the effect of thickness stretching on the mechanical behavior of advances thick structures.For example, Sobhy and Radwan [57] studied the effect of the thermal environment on the mechanical response of FG nanoplates by employing a new quasi-3D nonlocal hyperbolic plate theory (HSDT).Zenkour [58] proposed a quasi-3D HSDT based on trigonometric functions of in-plane and transverse displacements for an exponentially graded thick rectangular plate.Analysis on the effect of porosity, the biaxially oscillating loading and longitudinal magnetic field on the dynamic instability of viscoelastic nanoplates is performed by Jalaei and Thai [59].Daikh et al. [60] presented a new nonlocal quasi-3D hyperbolic shear theory to study the bending behavior of sigmoid functionally graded sandwich nanoplates posed on variable Winkler elastic foundation.They have also used various quasi-3D theories for static and dynamic stability responses of multilayer FG-CNTRC thick laminated plates.Amir et al. [61] used the quasi-3D HSDT in conjunction with modified couple stress theory to study the vibration behavior of a multi-layered FG porous micro plate.By applying the nonlocal quasi-3D hyperbolic shear deformation theory, Shahraki et al. [62] carried out a buckling and vibration analysis of FG-CNTRC thick nanoplates using the generalized differential quadrature method (GDQM).Buckling, bending and free vibration of FG-CNTRC thick laminated plates is carried out by Khadir et al. [63] and Alazwari [64] using novel quasi-3D HSDT with four-unknowns and five-unknowns, respectively.
Based on the literature survey, it is noted that the available research has not considered a coupling of microstructure and length scale influence on bending, stresses and buckling stability of multilayer FG-GRNC nanoplate.Therefore, the current article will cover this point and present impacts of length scale parameter (nonlocal), material scale parameter (gradient), distribution pattern, the GPL weight fraction, thickness stretching, geometry and size of GPLs, the geometry of the plate and the total number of layers on the stresses, deformation and critical buckling loads.This manuscript is organized as follows: Section 2 illustrates the gradation functions, Halpin-Tsai model, rule of mixtures, different sandwich configurations, kinematic fields, proposed plate theories, nonlocal statin gradient constitutive equation and equilibrium equation of the proposed model.Section 3 focuses on the solution procedure and analytical solutions for different boundary conditions (BCs) by using the Galerkin approach.Section 4 presents verification study and numerical results discussing influences of microstructure and length scale parameters on static deflection, stress distributions and buckling stability of FG-GRNC laminated plates.Most outcomes and conclusions are summarized in Section 5.

Material Distribution and Graduation
Consider a multilayer plate of length a, width b and thickness h, as shown in Figure 1a.Each layer is made from an isotropic polymer matrix and reinforced by graphene nanoplatelets (GPLs), and all layers have the same thickness.The weight fraction of GPLs varies linearly along the thickness from layer to layer.In the present analysis, as shown in Figure 1b, four different patterns of GPL distribution are considered, the uniform distribution of GPLs along the nanoplate thickness (UD) and three FG distributions of GPLs (FG-X GRNC, FG-O GRNC and FG-A GRNC).The FG-X GRNC plates have the maximum weight fraction on the bottom and the top layers and minimum one at the middle of the nanoplate.The FG-O GRNC plates have the opposite GPLs weight fraction toward the FG-X GRNC plates.The GPL weight fraction in the FG-A GRNC plates gradually increases from the top layer to the bottom layer.
FG-X GRNC plates.The GPL weight fraction in the FG-A GRNC plates gradually increases from the top layer to the bottom layer.are expressed as follows [17]: where  is the number of layers (Figure 1c) and g * is the total GPLs weight fraction.Using the Halpin-Tsai model, the effective Young's modulus of kth layer of the nanoplate is stated as [65]: where The weight fraction of GPLs g GPL of kth layer (k = 1, 2, . . ., N L ) for various patterns, are expressed as follows [17]: where N L is the number of layers (Figure 1c) and g * GPL is the total GPLs weight fraction.Using the Halpin-Tsai model, the effective Young's modulus of kth layer of the nanoplate is stated as [65]: where E m and E GPL are Young's modulus of polymer matrix and the GPLs reinforcement, respectively.The effective Poisson's ratio υ e and mass density ρ e of kth layer of the nanocomposite plate are computed using the rule of mixture and are expressed as follows: υ GPL and υ m are Poisson's ratio of the GPLs reinforcement and the epoxy matrix.V (k) GPL is the GPLs volume fraction of kth layer, and GPL is related to the matrix.

Basic Equations
To describe the field of displacement, a novel hyperbolic shear theory is proposed which contains the classical, first or higher-order plate terms as follows (10) in which u, v and w are the displacements along x, y and z directions, respectively.u 0 , v 0 and w 0 are the horizontal and the vertical displacements at the midplane of the plate.ψ x and ψ y are the rotation of the transverse normals around the x and y axes, respectively.Φ(z) is a shape function for the shear distribution.Four functions presented in Table 1 are selected to evaluate the accuracy and reliability of the present analysis.The variation of the used shape functions through the thickness direction is illustrated in Figure 2.  The strain fields of the nanocomposite plate have the following form: The strain fields of the nanocomposite plate have the following form: where the primary strains ε (0) ij and γ (0) ij can be written as functions in terms of derivative of displacement field as

Nonlocal Strain Gradient Elasticity Theory
By including the coupling physical impact of the strain gradient stress and nonlocal elastic stress fields, Lim et al. [70] proposed a function of stresses as Here, σ (0) ij is the classical stress corresponds to strain ε kl and the higher-order stress σ (1) ij corresponds to strain gradient ε kl,x , and they are given by the following form [70]: where C ijkl represent the fourth-order elastic coefficients and l (nm) is defined as the material length scale parameter presented to reflect the significance of the strain gradient stress field.e 0 a and e 1 a (nm) are the nonlocal parameters presented to reflect the significance of the nonlocal elastic stress field.α 0 (x, x , e 0 a) and α 1 (x, x , e 1 a) are the nonlocal attenuation functions developed by Eringen [29].Applying the linear differential operators for e = e 0 = e 1 on both sides of Equation ( 13) and rearranging the expression yield where ∇ 2 represents the Laplacian operator.The total nonlocal strain gradient constitutive relations can be expressed as, Daikh et al. [39] 1 where µ = (ea) 2 and λ = l 2 .
Using the generalized quasi-3D shear deformation theory, the nonlocal strain gradient constitutive stress-strain relations are governed by As the normal strain in z-direction can be discarded (non-stretching), ε (0) zz = 0, and hence However, to include the stretching through thickness, the normal strain ε (0) zz = 0, then Q ij are the 3D elastic constants: where (k) e are the Lamé's coefficients.

Equilibrium Equations
The variational principle is utilized to obtain the equilibrium equations of the GRNC plates which state that δ(U − V) = 0 (22) where U and V are the strain energy and potential energy of the applied loads, respectively.The total strain energy of the FG-GRNC laminated plate can be expressed as The variation of the potential energy of the applied loads can be given as Here, q presents the external transverse applied load.N 0 xx and N 0 yy are the in-plane applied loads in x-direction and y-direction, respectively.
The force and moment resultants may be portrayed by Substituting Equations ( 12) and ( 19) into the variational form in Equation ( 22) yields Then, the equilibrium equation of the GRNC plate is derived as follows: Substituting Equation ( 24) into the constitutive stress-strain equation described by Equation ( 19) give in which the composite stiffnesses can be expressed as Therefore, the governing equilibrium equation of GRNC plates based on the nonlocal strain gradient theory can be stated as

Analytical Solution
Based on the five variable quasi-3D HSDT in conjunction with nonlocal strain gradient theory, an analytical solution is developed by employing the Galerkin approach considering various boundary conditions.Galerkin expressions of displacements can be stated as: where U mn , V mn , W mn , X mn , Z mn and Ψ mn are arbitrary parameters.m and n are mode numbers.The functions X m (x) and Y n (y) that satisfy the simply supported and clamped boundary conditions are expressed in Table 2 as follows: By substituting Equation (34) in Equations ( 28)-( 33), one obtains The elements K ij of the matrix [K] are expressed in detail in Appendix A Table 2.The admissible functions X m (x) and Y n (y) for different boundary conditions.

Boundary Conditions
The Functions X m and Y n The applied load has a sinusoidal distribution, and can be written as [71],

Numerical Results
The present nanocomposite plate is made of a mixture of epoxy as a matrix and GPLs as reinforcement.The plate is exposed to in-plane load in two directions for buckling analysis, and transverse sinusoidal load for the bending analysis.
In this analysis, the thickness, length and width of the nanoplate are h = 20 nm, a = 200 nm and b = 200 nm, and it is made of GNLS with a length a GPL = 3 nm, thickness h GPL = 0.7 nm and width w GPL = 1.8 nm, [72].The normalized critical buckling load, displacements and normalized stresses are evaluated with the following forms.

Verification Analysis
To check the correctness and accuracy of the proposed method of solution and our shape function, a comparative evaluation between the obtained results and some valid predictions from the literature is carried out.The validation of the current model with those obtained by Wu et al. [73], Gholami and Ansari [74] and Thai et al. [18] are presented in Table 3.This analysis is related to the buckling of functionally graded graphene reinforced composite square plates with thickness h = 0.045 m and length a = 0.45 m.Table 3. Dimensionless critical buckling load ×10 −2 of GRC plate subjected to a compressive uniaxial in-plane loading (a GPL = 2.5 µm, h GPL = 1.5 nm, w GPL = 1.5 µm).On the other hand, more validations are made in Table 4 which summarizes the transverse displacement, axial stresses and shear stresses of an FG plate with volume fraction index p = 0 (fully ceramic).Young's modulus of the ceramic constituent (Zirconia ZrO 2 ) is E c = 151 GPa, while Poisson's ratio ν c = 0.3.The stretching effect is considered in these results.Zenkour [76] used the sinusoidal shear deformation theory as shape function and Navier procedure to solve the equilibrium equations.The same procedure is used by Bessaim et al. [77] using the hyperbolic theory.Neves et al. [78] applied Carrera's Unified Formulation and Murakami's Zig-Zag theory for the analysis.From the above results, it can be concluded that the proposed solutions for buckling or bending are in good agreement with the other results.

Parametric Study 4.2.1. Buckling Analysis
In this section, the effect of weight fraction of the constituents on the buckling behavior of the GRNC plates will be analyzed and investigated.Firstly, the effect of weight fraction g * GPL and various distributions of GPLs on the critical buckling load with and without stretching effect will be presented.
Firstly, in Table 5, the effect of g * GPL constituent and the reinforcement distribution on dimensionless buckling load of simply supported plates subjected to uniaxial compressive in-plane loads (χ 1 = −1.χ 2 = 0) is tabulated using various HSDTs.For this the critical buckling load increases.The FG-O GRCL plates have the lowest critical buckling load while the highest values are for the FG-X GRCL plates.The thickness stretching has an important effect on the response of the GRNC plate, especially for plates with 1% of GPLs.Table 6 illustrates the impact of in-plane charges (uniaxial and biaxial) on the critical buckling load of 1% of GPLs plates.It is clear that the obtained values for fully clamped plates (CCCC) are higher than the other plates; however, the lowest loads are for the simply supported plates (SSSS).Additionally, the difference between quasi-3D theories (ε zz = 0) and classical shear deformation theories (ε zz = 0) is significant for the fully clamped plates.
To capture the size-dependent effects, Table 7 summarizes the dimensionless critical buckling load of various patterns of reinforcement under the impact of length scale and nonlocal parameters.The results are reported here for the only weight fraction g * GPL = 1%.As it is reasonable to expect, the buckling response is affected by the nonlocal parameter, with a clear decrease in the stiffness and critical buckling load for an increasing nonlocal parameter µ.Unlike the nonlocality effect, the increase in the length scale parameter leads to an increment in the critical buckling load wherever the GPLs distribution patterns is.The influence of the number of layers N L , length-to-thickness ratio and aspect ratio on buckling response of simply supported GRNC plates subjected to biaxial in-plane compressive loads is presented in Table 8.The results are focused on the only FG-O distribution of GPLs with weight fraction g * GPL = 1%.It can be observed that the increase in the number of layers decreases the buckling load FG-O GRNC plates.In addition, the increase in the length-to-thickness ratio and aspect ratio leads to a decrement in the dimensional buckling load.The size-dependent influence on the critical buckling load of GRNC plates is controlled by including the nonlocal parameter µ and the length-scale parameter λ (see Appendix A).The results of the size effect are plotted in Figure 3. From the figure, the increase in the parameter µ reduces the stiffness of the GRNC plate and therefore leads to a decrement in the critical buckling load.Otherwise, unlike the nonlocality effect, the critical buckling load is increased by increasing the length-scale parameter λ.The size-dependent influence on the critical buckling load of GRNC plates is controlled by including the nonlocal parameter  and the length-scale parameter  (see Appendix A).The results of the size effect are plotted in Figure 3. From the figure, the increase in the parameter  reduces the stiffness of the GRNC plate and therefore leads to a decrement in the critical buckling load.Otherwise, unlike the nonlocality effect, the critical buckling load is increased by increasing the length-scale parameter .Dimensionless critical buckling load in terms of the number of layers for various GPLs reinforcement patterns is presented in Figure 4.The weight fraction is taken as g * GPL = 1%.Similar results are obtained for the patterns FG-A and UD where the number of layers does not change the results.For the FG-O pattern, the critical buckling load decreases by increasing the total number of layers.In contrast, the increase in the total number of layers decreases the buckling load, then tends to be constant for N ≥ 10.Dimensionless critical buckling load in terms of the number of layers for various GPLs reinforcement patterns is presented in Figure 4.The weight fraction is taken as g * = 1%.Similar results are obtained for the patterns FG-A and UD where the number of layers does not change the results.For the FG-O pattern, the critical buckling load decreases by increasing the total number of layers.In contrast, the increase in the total number of layers decreases the buckling load, then tends to be constant for  ≥ 10.Variation of the critical buckling load of GRNC plates for different values of the length-to-width ratio of GPL nanofillers is presented in Figure 5.The four GPLs reinforcement patterns are analyzed.It can be seen that the increase in the thickness-to-length ratio  ℎ ⁄ leads to an increment in critical buckling load.Variation of the critical buckling load of GRNC plates for different values of the lengthto-width ratio of GPL nanofillers is presented in Figure 5.The four GPLs reinforcement patterns are analyzed.It can be seen that the increase in the thickness-to-length ratio a GPL /h GPL leads to an increment in critical buckling load.
of layers for various GPLs reinforcement patterns (g = Variation of the critical buckling load of GRN length-to-width ratio of GPL nanofillers is presented ment patterns are analyzed.It can be seen that the inc  ℎ ⁄ leads to an increment in critical buckling To better comprehend the effect of GPLs geomet ling load versus the width-to-length ratio   ⁄ To better comprehend the effect of GPLs geometry, Figure 6 depicts the critical buckling load versus the width-to-length ratio b GPL /a GPL .The plotted results show that the widthto-length ratio has a significant effect on the critical buckling load for values less than 15, where critical buckling load increase.The results are almost constant for greater values of ratio (≥15) regardless GPLs reinforcement pattern.
Figure 7 plots the critical buckling loads of GRNC plates versus the aspect ratio b/a subjected to compressive biaxial in-plane loads.The coefficient of the in-plane load in x-direction is assumed constant (χ 1 = 1), while the compressive loads in y-direction vary from 0 to 3. Except the case of χ 2 = 0, dimensionless critical buckling load increases as aspect ratio increases.In addition, the intensity of in-plane loads has an important impact on critical buckling load where increasing compressive in-plane loads lead to a decrement in buckling loads.
Figure 8 illustrates the influence of weight faction variation on dimensionless critical buckling load of GRNC plates under various boundary conditions and for various GPLs reinforcement patterns.The dimensionless critical buckling loads increase linearly by increasing a weight fraction g * GPL for all boundary conditions and reinforcement patterns.
width-to-length ratio has a significant effect on the critical buckling than 15, where critical buckling load increase.The results are almos values of ratio (≥ 15) regardless GPLs reinforcement pattern.Figure 7 plots the critical buckling loads of GRNC plates versus subjected to compressive biaxial in-plane loads.The coefficient of th direction is assumed constant ( = 1), while the compressive load from 0 to 3. Except the case of  = 0, dimensionless critical buckli aspect ratio increases.In addition, the intensity of in-plane loads has on critical buckling load where increasing compressive in-plane load in buckling loads.Figure 7 plots the critical buckling loads of GRNC plates versus subjected to compressive biaxial in-plane loads.The coefficient of th direction is assumed constant ( = 1), while the compressive load from 0 to 3. Except the case of  = 0, dimensionless critical buckli aspect ratio increases.In addition, the intensity of in-plane loads has on critical buckling load where increasing compressive in-plane load in buckling loads.

Bending Analysis
In this section, using various HSDTs, comprehensive analysis of displacements and stresses in multilayer FGNC plates are carried out.Table 9 shows the influence of thickness stretching on longitudinal displacement  , transverse displacement  and the stresses  ,  and  of square simply supported FGNC laminated plate with various patterns of reinforcement.Similar results are obtained by comparing our hyperbolic theory with the sinusoidal theory.It can be observed that the thickness stretching have a significant impact on the displacements of the plate.In addition, it can be seen that the

Bending Analysis
In this section, using various HSDTs, comprehensive analysis of displacements and stresses in multilayer FGNC plates are carried out.Table 9 shows the influence of thickness stretching on longitudinal displacement u, transverse displacement w and the stresses σ xx , σ xy and σ xz of square simply supported FGNC laminated plate with various patterns of reinforcement.Similar results are obtained by comparing our hyperbolic theory with the sinusoidal theory.It can be observed that the thickness stretching have a significant impact on the displacements of the plate.In addition, it can be seen that the difference between displacements u with and without the thickness stretching effect is very large for FG-A GRNC laminated plates.In Table 10, the effect of g * GPL constituent and different reinforcements distribution on the dimensionless central deflection of GRNC plates subjected to various boundary conditions is tabulated using the present theory.The increase in the g * GPL reinforced percentage tends to decrease the transverse displacements.The FG-O GRCL plates have the highest central deflections while the lowest values are for the FG-X GRCL plates.The thickness stretching has an important effect on the response of the GRNC plate, especially for the fully clamped plates (CCCC).
Table 11 shows the size-dependent effects on the dimensionless central of GRNC plates with various patterns of reinforcement where the impact of length scale and nonlocal parameters is studied.The results are here reported for the only weight fraction g * GPL = 1%.From this table, the static bending response is affected by the nonlocal parameter, increasing nonlocal parameter µ leads to a decrement in plate stiffness and subsequently an increment in transverse displacement.Unlike the nonlocality effect, the increase in the length scale parameter leads to a decrement in transverse displacement wherever the GPLs distribution patterns remain constant.In Figure 9, dimensionless central deflection versus the total number of layers for various GPL reinforcement patterns is depicted for weight fraction g * GPL = 1%.Unlike the previous results of buckling response obtained in Figure 9, the dimensionless central deflection increases by increasing the total number of layers for FG-O pattern and decreases for the FG-X pattern.In Figure 9, dimensionless central deflection versus the total number of layers for various GPL reinforcement patterns is depicted for weight fraction g * = 1%.Unlike the previous results of buckling response obtained in Figure 9, the dimensionless central deflection increases by increasing the total number of layers for FG-O pattern and decreases for the FG-X pattern.In Figure 10, the impact of nonlocal and small-scale parameters on the central deflection of simply supported GRNC plates for various GPLs reinforcement patterns is investigated.As we demonstrated in the buckling analysis, the inclusion of nonlocal parameter  reduces the plate stiffness, so, transverse displacement increases.In addition, the increase in small-scale parameter  leads to a decrement in transverse displacement.In Figure 10, the impact of nonlocal and small-scale parameters on the central deflection of simply supported GRNC plates for various GPLs reinforcement patterns is investigated.As we demonstrated in the buckling analysis, the inclusion of nonlocal parameter µ reduces the plate stiffness, so, transverse displacement increases.In addition, the increase in small-scale parameter λ leads to a decrement in transverse displacement.In Figure 11, dimensionless central deflection versus the length-to-thickness ratio of the GPLs nanofillers using various GPLs reinforcement patterns is demonstrated.It is clear that the increase in the thickness-to-length ratio  ℎ ⁄ leads to a decrement in dimensionless central deflection.In Figure 11, dimensionless central deflection versus the length-to-thickness ratio of the GPLs nanofillers using various GPLs reinforcement patterns is demonstrated.It is clear that the increase in the thickness-to-length ratio a GPL /h GPL leads to a decrement in dimensionless central deflection.In Figure 11, dimensionless central deflection versus the length-to-thickness ratio of the GPLs nanofillers using various GPLs reinforcement patterns is demonstrated.It is clear that the increase in the thickness-to-length ratio  ℎ ⁄ leads to a decrement in dimensionless central deflection.Figure 12 depicts the dimensionless central deflection influenced by width-to-length ratio b GPL /a GPL .It can be observed that the width-to-length ratio has a significant effect on the dimensionless central deflection for values less than 15, where central deflection decreases.The results are almost constant for greater values of ratio (≥ 15) regardless of the GPL reinforcement pattern.Figure 13 plots the influence of weight fraction g * on the transverse displacem of GRNC plates subjected to various boundary conditions.The highest deflections are SSSS plates and pure Epoxy (g * = 0.0%).Figure 14 shows the variation of the in-plane displacement through-the-thickness of various types of GPLs reinforcements.The plate is assumed simply supported.In the upper-half plate, the in-plane displacements of the FG-O GRNC plates are found to be of the smallest magnitudes and that of the FG-A GRNC plates, of the largest magnitudes.However, in the lower-half plate, the in-plane displacements of the FG-A GRNC plates are found to be of the smallest magnitudes.
In Figure 15, dimensionless transverse displacement w through-the-thickness of the GRNC plate with simply supported boundary conditions for various GPLs reinforcement patterns is investigated.It can be observed that the maximum transverse displacement occurs at a point on the mid-plane of the GRNC plate and their magnitudes for the FG-O pattern are larger than that for the other patterns.Therefore, in order to make the most decrease in the dimensionless transverse displacement w, it is better to put the GNLs nanofiller as far as away from the middle axis of the GRNC plate which creates more flexural rigidity.port Figure 14 shows the variation of the in-plane displacement through-the-thic various types of GPLs reinforcements.The plate is assumed simply supported.In per-half plate, the in-plane displacements of the FG-O GRNC plates are found to b smallest magnitudes and that of the FG-A GRNC plates, of the largest magnitude ever, in the lower-half plate, the in-plane displacements of the FG-A GRNC pl found to be of the smallest magnitudes.In Figure 15, dimensionless transverse displacement  through-the-thickne GRNC plate with simply supported boundary conditions for various GPLs reinfo patterns is investigated.It can be observed that the maximum transverse displ occurs at a point on the mid-plane of the GRNC plate and their magnitudes for th pattern are larger than that for the other patterns.Therefore, in order to make t decrease in the dimensionless transverse displacement  , it is better to put the GN ofiller as far as away from the middle axis of the GRNC plate which creates more rigidity.Dimensionless stresses for different types of the GRNC plate are plotted in Figure 17.It can be found that the stresses σ xx and σ zz are tensile at the top surface and compressive at the bottom surface, unlike the σ xy stresses.For the UD, FG-O and FG-X patterns, the stresses, σ xx , σ xy and σ zz are null at the mid-plane of the plate because of the symmetric distribution of the GPLs nanofillers.For the plates with a symmetric distribution of the GPLs, FG-O and UD patterns, the maximum shear stresses τ xz occur at a point on the mid-plane of the GRNC plate.For the FG-X GRNC plate, considerable reduction of stresses τ xz at the mid-plane of the plate due to the absence of GPLs nanofillers at the central layers.

Conclusions
Novel quasi-3D hyperbolic higher order shear deformation theory in conjunction with modified continuum nonlocal strain gradient theory is used for bending and buck ling analyzes of FG GRNC plates.The governing equations are derived by applying Ham ilton's principle and solved analytically using the proposed Galerkin method.The effec tive Young's modulus of the GRNC plate was estimated using the Halpin-Tsai mode and Poisson's ratio and mass density are computed by using the rule of mixture.Fou different patterns of GPL distribution are considered in this study.The feasibility of th proposed solutions is effectuated by comparing them with the existing analytical solu tions.The numerical examples show that:

Conclusions
Novel quasi-3D hyperbolic higher order shear deformation theory in conjunction with modified continuum nonlocal strain gradient theory is used for bending and buckling analyzes of FG GRNC plates.The governing equations are derived by applying Hamilton's principle and solved analytically using the proposed Galerkin method.The effective Young's modulus of the GRNC plate was estimated using the Halpin-Tsai model, and Poisson's ratio and mass density are computed by using the rule of mixture.Four different patterns of GPL distribution are considered in this study.The feasibility of the proposed solutions is effectuated by comparing them with the existing analytical solutions.The numerical examples show that:

•
Unlike the FG-X GRNC plates, the increase in the number of layers leads to an improvement in the FG-X GRNC plate stiffness; therefore, the critical buckling load increases and the central deflection decreases.

•
Both increasing the weight fraction and the number of layers can improve the stiffness of the GRNC plates.

•
The inclusion of thickness stretching has a significant effect on the response of the GRNC plate, especially for the fully clamped plate.

•
The inclusion of a nonlocal parameter leads to a decrease in buckling loads and an increase in transverse displacements, while the opposite is found when increasing the length-scale parameter.

•
The FG-X reinforcement pattern has the better response under in-plane and transverse loadings, among the other analyzed distributions, due to their higher stiffness.

•
Increasing the length-to-width and length-to-thickness ratios result in increasing the total stiffness of the GRNC plate and consequently lead to higher critical buckling loads and lower transverse displacements.
This work plays a significant role in further studies in composite plates using the 3D elasticity model.

Figure 2 .
Figure 2. Variation of shape functions () through the thickness direction.

Figure 4 .
Figure 4. Dimensionless critical buckling load of simply supported GRNC plate in terms of number of layers for various GPLs reinforcement patterns (g * = 1%).

Figure 4 .
Figure 4. Dimensionless critical buckling load of simply supported GRNC plate in terms of number of layers for various GPLs reinforcement patterns g * GPL = 1% .

Figure 8
Figure8illustrates the influence of weight faction variation on dimensionless critical buckling load of GRNC plates under various boundary conditions and for various GPLs reinforcement patterns.The dimensionless critical buckling loads increase linearly by increasing a weight fraction g * for all boundary conditions and reinforcement patterns.

Figure 9 .
Figure 9. Dimensionless central deflection of simply supported GRNC plate in terms of number of layers for various GPLs reinforcement patterns (g * = 1%).

Figure 9 .
Figure 9. Dimensionless central deflection of simply supported GRNC plate in terms of number of layers for various GPLs reinforcement patterns g * GPL = 1% .

Figure 10 .
Figure 10.Effect of nonlocal and small-scale parameters on the dimensionless central deflection of simply supported GRNC plate for various GPLs reinforcement patterns (g * = 1%).

Figure 10 .
Figure 10.Effect of nonlocal and small-scale parameters on the dimensionless central deflection of simply supported GRNC plate for various GPLs reinforcement patterns g * GPL = 1% .

Figure 10 .
Figure 10.Effect of nonlocal and small-scale parameters on the dimensionless central deflection of simply supported GRNC plate for various GPLs reinforcement patterns (g * = 1%).

Figure 12
Figure 12 depicts the dimensionless central deflection influenced by width-to-length ratio   ⁄ .It can be observed that the width-to-length ratio has a significant effect on the dimensionless central deflection for values less than 15, where central deflection decreases.The results are almost constant for greater values of ratio (≥ 15) regardless of the GPL reinforcement pattern.

Figure 13 plots
Figure 13 plots the influence of weight fraction g *on the transverse displacement of GRNC plates subjected to various boundary conditions.The highest deflections are for SSSS plates and pure Epoxy (g * = 0.0%).

Figure 13 plots
Figure 13  plots the influence of weight fraction g * GPL on the transverse displacement of GRNC plates subjected to various boundary conditions.The highest deflections are for SSSS plates and pure Epoxy g * GPL = 0.0% .

Figure 14 . 28 Figure 15 .
Figure 14.Dimensionless in-plane displacement u through the thickness of simply supported GRNC plate for various GPLs reinforcement patterns g * GPL = 1% .

Figure 16
Figure 16 contains the plots of dimensionless transverse displacement of GRNC with uniformly distributed GPLs by considering the weight fraction impact.The m mum values of the dimensionless transverse displacement are for the pure Epoxy (g * = 1%), whereas the minimum values are for plate with weight fractions g * =

Figure 15 .
Figure 15.Dimensionless transverse displacement w through the thickness of simply supported GRNC plate for various GPLs reinforcement patterns g * GPL = 1% .

Figure 16
Figure 16  contains the plots of dimensionless transverse displacement of GRNC plate with uniformly distributed GPLs by considering the weight fraction impact.The maximum values of the dimensionless transverse displacement are for the pure Epoxy plate g * GPL = 1% , whereas the minimum values are for plate with weight fractions g * GPL = 1%.

Figure 16 containsFigure 16 .Figure 16 .
Figure 16 contains the plots of dimensionless transverse displacement of GRNC p with uniformly distributed GPLs by considering the weight fraction impact.The m mum values of the dimensionless transverse displacement are for the pure Epoxy p (g * = 1%), whereas the minimum values are for plate with weight fractions g * =

n dxdy − 2 (∂ 6 b 0 ∂ 4 X m ∂x 4 ∂ 2 ∂x 4 ∂x 6 ∂ 2 X m ∂x 2 ∂ 4 ∂ 4 2D s 66 a 0 b 0 ∂ 2 X m ∂x 2 ∂ 2 2D s 66 a 0 b 0 ∂ 4 X m ∂x 4 ∂ 2 Y n ∂y 2 b 0 ∂ 2 ∂ 2 0 ∂ 2 X m ∂x 2 ∂ 3 ∂ 2 ∂ 4 b 0 ∂ 2 X m ∂x 2 ∂ 2 Y n ∂y 2 X∂ 2 X m ∂x 2 ∂ 2 ∂ 2 ∂ 4 2 a 0 b 0 ∂ 2 X m ∂x 2 ∂ 2 ∂ 2 ∂ 2 a 0 b 0 ∂ 2 X m ∂x 2 ∂ 2 ∂y 2
D 12 + 2D 66 ) Y n ∂y 6 X m Y n dxdy + a 0 Y n ∂y 2 X m Y n dxdy −2(D 12 + 2D 66 ) Y n X m Y n dxdy + D s 12 + 2D s Y n X m Y n dxdy + D s 11 + D s 12 + 2D s Y n ∂y 4 X m Y n dxdy K 35 = D s Y n ∂y 4 X m Y n dxdy + D s 12 + Y n ∂y 2 X m Y n dxdy −λ D s 12 + X m Y n dxdy + D s n X m Y n dxdy − µ a 0 X m ∂x 2 Y n X m Y n dxdy + Y n ∂y 2 X m Y n dxdy K 41 = B s Y n ∂y 3 X m ∂Y n ∂y dxdy + B s X m ∂x 2 Y n X m Y n dxdy − λ X m ∂x 4 Y n X m Y n dxdy + a 0 Y n ∂y 2 X m Y n dxdy + X m ∂x 2 Y n X m Y n dxdy + X m ∂x 4 Y n X m Y n dxdy + Y n ∂y 2 X m Y n dxdy + X m ∂x 2 Y n X m Y n dxdy − λ Y n ∂y 2 X m Y n dxdy − λ   Y n ∂y 2 X m Y n dxdy + X m Y n dxdy

Table 1 .
Shear shape function distribution.

Table 5 .
Dimensionless critical buckling load of GRC plate subjected to a compressive uniaxial in-plane loading influenced by weight fraction g * GPL .

Table 6 .
Dimensionless critical buckling load of GRNC plate under various boundary conditions subjected to a compressive uniaxial/biaxial in-plane loading g * GPL = 1% .

Table 7 .
Effect of nonlocal and small-scale parameters on dimensionless critical buckling load of GRC plate subjected to a compressive uniaxial in-plane loading g * GPL = 1% .

Table 8 .
Effect of the total number of layers on dimensionless critical buckling load of FG-O GRC plate subjected to a compressive biaxial in-plane loading g * GPL = 1% .

Table 9 .
Dimensionless displacements and stresses of simply supported square GRC plate for various GPLs reinforcement patterns g * GPL = 1% .

Table 10 .
Dimensionless displacements of GRC plate under various boundary conditions influenced by weight fraction g * GPL .

Table 11 .
Effect of nonlocal and small-scale parameters on dimensionless central deflection of GRC plate subjected to a compressive uniaxial in-plane loading g * GPL = 1% .