Nonlinear Vibration of a Functionally Graded Nanobeam Based on the Nonlocal Strain Gradient Theory considering Thickness Effect

In this work, a nonlocal strain gradient beam model considering the thickness effect is developed to study the nonlinear vibration response of a functionally graded nanobeam. e governing equation of the functionally graded nanobeam is derived by using the Euler–Bernoulli beam theory with von Kármán’s nonlinear strain-gradient relationship and the Hamilton principle. e expression of the nonlinear frequency for the functionally graded nanobeam with pinned-pinned boundary conditions is obtained with the help of Galerkin technique and the Hamiltonian approach. e obtained results show that the effect of thickness is very important for the size-dependent vibration response of the functionally graded nanobeam; the nonlinear vibration response of the nanobeam depends not only on the material length scale parameter and nonlocal parameter but also on the slenderness ratio. Effects of the slenderness ratio and the power-law index on the vibration response of the functionally graded nanobeam are also investigated and discussed. e numerical results show that the nonlocal parameter reduces the nonlinear frequency of the functionally graded nanobeam, while the material length scale parameter increases the nonlinear frequency of the functionally graded nanobeam. e slenderness ratio leads to an increase in the nonlinear frequency of the functionally graded nanobeam, while the power-law index leads to a decrease in the nonlinear frequency of the functionally graded nanobeam.


Introduction
Because of wide applications in areas such as micro-/nanoelectromechanical systems (MEMs/NEMs), biosensors, micro-/nanosensors, and nanoactuators, the behaviors of micro-/nanostructures have attracted attention of many scientists. Experiment observations showed that the sizedependent influences of the micro-/nanostructures are very important and cannot be ignored [1,2]. Some higher-order elasticity theories that can capture the size-dependent effects on the behaviors of the micro-/nanostructures have been developed.
e nonlocal elasticity theory introduced by Eringen [3,4] states that the stress at a point is considered as a function of strains at all points in the continuum body. is theory is widely used to analyze responses of the structures at the nanoscale. Based on the nonlocal elasticity theory, a number of works on the static and dynamic behaviors of beams have been published [5][6][7][8][9][10][11]. e vibration and stability problems of the carbon nanotubes based on the nonlocal elasticity theory have also attracted attention of many authors [12][13][14][15][16][17][18][19][20]. It can be seen that the nonlocal elasticity theory has been employed very effectively for modeling the nanostructures; however, research results show that only the stiffness-softening effect is observed by using this theory. e strain gradient theory (SGT) [21][22][23][24][25] states that the total stress should consider some additional strain gradient terms by modeling small-scaled materials as atoms with the higher-order deformation mechanism instead of collection of points. By modifying the SGT, Yang et al. [26] proposed the modified couple stress theory (MCST); this theory considers the strain energy density as a function of both the strain tensor conjugated with the stress tensor and the curvature tensor conjugated with the couple stress tensor. In 2003, with the similar approach as Yang et al. [26], Lam et al. [27] introduced the modified strain gradient theory (MSGT) with three material length scale parameters. Other than the nonlocal elasticity theory, the SGT, MCST, and MSGT play stiffness-hardening effect on responses of the micro-/ nanostructures. Based on the SGT and MCST, many works related to the static and dynamic responses of beams were carried out [28][29][30][31][32][33][34][35][36][37].
It can be seen that the nonlocal elasticity theory and the SGT describe size-dependent nano-/micromechanical characteristics of materials in two different manners. In 2015, Lim et al. [38] combined the nonlocal elasticity theory and the strain gradient theory into a generalized higher-order elasticity theory called the nonlocal strain gradient theory (NSGT). Effects of the two length scales including the material length scale parameter and the nonlocal parameter on the mechanical and physical behaviors of size-dependent structures are described in the NSGT. e NSGTcan be considered the most generalized theory of the higher-order elasticity theories that can describe the behaviors of size-dependent structures. Many articles related to the static and dynamic behaviors of micro-/nanostructures have been published based on the NSGT. Based on the NSGT, Li and Hu [39] investigated the buckling problem of size-dependent nonlinear beams. A size-dependent sinusoidal shear deformation beam model was developed by Lu et al. [40] for analyzing the free vibration response of nanobeams. e NSGT was used by Li et al. [41] to examine size-dependent effects on critical flow velocity of fluid-conveying microtubes. Bahaadini et al. [42] investigated free vibration and instability behaviors of a nanotube conveying nanoflow based on the NSGT and the Timoshenko beam theory. e nonlinear dynamical behavior of a fluid-conveying viscoelastic microtube was studied by Ghayesh and Farokhi using the NSGT and the Euler-Bernoulli beam theory [43]. e NSGT was employed by Arefi et al. to investigate the size-dependent bending response of a sandwich porous nanoplate integrated with piezomagnetic face-sheets [44] and the magneto-electroelastic vibration response of a porous functionally graded core sandwich nanoplate with piezomagnetic face-sheets and resting on an elastic foundation [45].
First discovered by scientists in Japan in 1984, the functionally graded materials (FGMs) have some advantages over the traditional composite materials. FGMs are increasingly being used extensively in the construction of structures that work in environments requiring high temperature resistance, corrosion resistance, improved stress spreading, and inferior stress intensity factors [46]. ere is a large number of works analyzing mechanical behaviors of functionally graded (FG) structures. Ansari et al. [47] investigated free vibration characteristics of FG microbeams based on the strain gradient Timoshenko beam theory. Based on the MCST and the Euler-Bernoulli beam theory, Arbind and Reddy [48] studied nonlinear vibration responses of FG microbeams. Nonlinear free vibration and postbuckling behaviors of FG beams resting on the nonlinear elastic layer were examined by Fallah and Aghdam [49]. Bending and vibration responses of the nonlocal FG beams were reported by ai and Vo [50] using various higher-order shear deformation beam theories. e NSGT [38] was also used very effectively to analyze the vibration response of FG beams. Simsek [51] investigated the nonlinear free vibration of a FG nanobeam based on the NSGT and the Euler-Bernoulli beam theory. Based on the NSGT, size-dependent nonlinear Euler-Bernoulli and Timoshenko beam models were developed by Li et al. [52] to study bending and free vibration characteristics of FG beams. e NSGT was employed by Allam and Radwan [53] to investigate bending, buckling, and vibration responses of viscoelastic FG curved nanobeam resting on an elastic foundation. Surface effects on nonlinear vibration behavior of an electrostatic FG nanoresonator were investigated by Esfahani et al. [54] using the NSGT and the Euler-Bernoulli beam theory. e effect of the neutral surface on the electroelastic analysis of FG piezoelectric plate resting on the Winkler-Pasternak foundation was investigated by Arefi et al. [55] in the framework of the two-variable sinusoidal shear deformation theory. Arefi et al. [56] studied the free vibration of a sandwich nanoplate including the FG core and piezoelectric face-sheets by using the two-variable sinusoidal shear deformation theory. e NSGT [38] is being used more and more extensively in examining the size-dependent effects on the mechanical behaviors of micro-/nanostructures. However, in previous publications, the size-dependent effects are often assumed to be neglected in the thickness direction. Recently in 2018, Li et al. [57] showed that the size-dependent effect of thickness on the mechanical behaviors of structures is very important. In [57], the authors developed a nonlocal strain gradient beam model incorporating the thickness to investigate the buckling problem of nanobeams. When considering the effect of thickness, the stiffness-softening and stiffnesshardening effects depend not only on the ratio of the stressgradient parameter to the strain gradient parameter but also on the geometric feature (slenderness ratio) [57]. Nonlinear free vibration response of homogeneous nanobeams was investigated by Chen et al. [58] using the NSGT incorporating the thickness effect.
According to authors' knowledge, nonlinear vibration behavior of the FG nanobeam based on the NSGT and considering the influence of the thickness is not yet announced. us, the nonlinear vibration response of the FG nanobeam based on the NSGT incorporating the thickness effect will be investigated in this paper. e governing equation of the FG nanobeam is derived by using the Euler-Bernoulli beam theory with von Kármán's nonlinear strain-gradient relationship and the Hamilton principle. With the help of Galerkin technique and the Hamiltonian approach, the expression of nonlinear frequency for the FG nanobeam with pinned-pinned boundary conditions is obtained in a closed form. e results show that the effect of thickness is very important to the size-dependent vibration response of the FG nanobeam; the nonlinear vibration response of the FG nanobeam depends not only on the material length scale parameter and nonlocal parameter but also on the slenderness ratio. Effects of the slenderness ratio 2 Advances in Civil Engineering and the power-law index on the vibration response of the FG nanobeam are also investigated and discussed in this paper.

Modeling of the Functionally Graded Nanobeam.
A FG nanobeam of length L, width d, and height h is considered as in Figure 1. e material properties of the nanobeam vary as a function of the thickness coordinate. Based on the rule of mixture, the material properties P of the nanobeam including Young's modulus E and the material density ρ are expressed as where P U and P L are the material properties at the upper and lower surfaces of the nanobeam, respectively, and V U and V L are the volume fractions at the upper and lower surfaces of the nanobeam, respectively. e volume fractions V U and V L are supposed as Here, k is the power-law index which prescribes the material variation and z ms is the thickness coordinate from the geometry neutral surface. From equations (1)-(3), the material properties can be obtained as Because the material properties of the FG nanobeam are not symmetric about its geometry neutral surface, the geometry and the physical neutral surface of the FG nanobeam are not the same. e position of the physical neutral surface can be determined as [51] c � h/2 − h/2 E z ms z ms dz ms For a rectangular FG nanobeam, the position of the physical neutral surface can be expressed as [52] c � As in Figure 1, the thickness coordinate from the physical neutral surface, z, is us, in the new coordinate system, by substituting equation (7) into equation (4), the material properties of the nanobeam can be expressed as Based on equation (8), Young's modulus E and the material density ρ of the FG nanobeam are determined as It can be seen from equations (9) and (10)

2.2.
e Nonlocal Strain Gradient eory. Based on the NSGTproposed by Lim et al. [38], the total stress of the beam is defined as where ∇ � z/zx is the differential operator and σ xx and σ (1) xx are the nonlocal stress and the higher-order nonlocal stress, respectively, which are expressed as In equations (12) and (13), E is Young's modulus, l m is the material length scale parameter or the strain gradient length scale parameter, ε xx is the classical strain, ε xx,x is the strain gradient, α 0 and α 1 are the two nonlocal kernel functions satisfying the conditions given by Eringen [3,4], and e 0 a and e 1 a are the two nonlocal parameters.
For the one-dimensional problem, by applying the linear operators L i � 1 − (e i a) 2 ∇ 2 for i � 0, 1 into equation (11), the general nonlocal strain gradient constitutive equation can be obtained as [38] 1 in which ∇ 2 � z 2 /zx 2 is the Laplace operator. When considering e 0 � e 1 � e, equation (14) can be reduced as By letting l m � 0, equation (15) becomes which is the nonlocal constitutive equation for the nonlocal elasticity theory [3,4].
Geometry neutral surface

Advances in Civil Engineering
Also, on the contrary, the strain gradient constitutive equation for the strain gradient theory [25] can be obtained by letting ea � 0:

e Governing Equation for a FG Nanobeam considering ickness
Effect. e displacements of the FG nanobeam based on the Euler-Bernoulli beam theory take the form According to von Kármán's geometrical relationship, the nonzero strain of the nanobeam is given as where u and w are the axial and transverse displacements of any point on the physical neutral surface, respectively. We can see from equation (19) that the strain of the nanobeam depends not only on x but also on z; thus, Li et al. [57] proposed that Laplace operator ∇ 2 � z 2 /zx 2 in Section 2.2 should be replaced by ∇ 2 � z 2 /zx 2 + z 2 /zz 2 . According to the suggestion by Li et al. [57], the size-dependent effect of thickness will be considered for analyzing the mechanical behavior of the FG nanobeam. By considering the thickness effect, the total stress of the FG nanobeam can be expressed as [57] where σ xx is the nonlocal stress and σ (1) xxx and σ (1) xxz are the higher-order nonlocal stresses. e virtual strain energy can be given as [57] δU Substituting the expression of the strain in equation (19) into equation (21), we get e following resultants are introduced: where N and M are the resultant force and moment for axial and transverse directions, respectively, N x and N z are the higher-order resultant forces for axial and transverse directions, respectively, and M x is the higher-order resultant moment. Employing equations (22) and (23), the virtual strain energy can be obtained as For the FG nanobeam, the resultants can be extended from the results in [57] as follows: where e virtual kinetic energy of the FG nanobeam taking into account both axial and transverse motions is given as where m 0 and m 2 are the mass moments of inertia which are defined as e virtual work caused by the distributed transverse Hamilton's principle states that Substituting equations (24), (31), and (33) into equation (34) and after some mathematical manipulations, we can obtain the following equations: e boundary conditions at x � 0 and x � L can be expressed as From equations (36) and (37), we get e expression for N * can be extended from the result for the homogeneous beam obtained by Chen et al. [58] as follows: From equations (27)- (29), we obtain Advances in Civil Engineering Now, substituting equations (44) and (45) into equation (42), the governing equation of motion for the FG nanobeam can be obtained as follows: In the case of pinned-pinned boundary conditions, the classical and nonclassical boundary conditions for the nanobeam at x � 0 and x � L are [57] w Equation (46) is the governing equation of motion for the FG nonlocal strain gradient nanobeam considering the thickness effect. For the homogeneous nanobeam, i.e., D xx � EI, A xx � EA, m 0 � ρA, and the second mass moment of inertia m 2 � ρI is ignored, equation (46) becomes the governing equation of motion for the homogeneous nonlocal strain gradient nanobeam considering the thickness effect developed by Chen et al. [58]: When the thickness effect is not taken into account, equation (46) is reduced to the model of the FG Euler-Bernoulli beam developed by Li and Hu [52]: For convenience, the following dimensionless variables are introduced: For the purpose of studying the free vibration, the transverse distributed load is ignored, and using equation (50), the equation of motion for the FG nanobeam can be rewritten in a dimensionless form as follows: where 6

Advances in Civil Engineering
e classical and nonclassical boundary conditions at x � 0 and x � 1 become Now, we will apply the Galerkin method to convert partial differential equation (51) into the ordinary differential one. By this method, the displacement function is assumed to have the following form: in which V(t) is an unknown time-dependent function and needs to be determined, and ϕ(x) is a shape function that satisfies the boundary conditions of the nanobeam. For the pinned-pinned nanobeam, the shape function can be chosen as Applying the Galerkin method, equation (51) is reduced to the following ordinary differential equation: where coefficients c 1 and c 2 are given as where ϕ (6) � d 6 ϕ/dx 6 , ϕ (4) � d 4 ϕ/dx 4 , ϕ ″ � d 2 ϕ/dx 2 , and ϕ ′ � dϕ/dx.

Solution Procedure
In this section, we will find the approximate solution of equation (56). As can be seen that equation (56) is a Duffing oscillator, there are many techniques to find the approximate solution of this equation. e Hamiltonian approach proposed by He [59] can be applied to find the approximate solution of equation (56). First, we assume that equation (56) satisfies the following initial conditions: where V 0 � w(0.5) is the dimensionless maximum vibration amplitude of the nanobeam. According to the Hamiltonian approach [59], the Hamiltonian of equation (31) can be constructed as Integrating equation (60) with respect to t from 0 to T/4, we get: where T is the period of oscillation.
Assume that the solution of equation (56) can be given as Here, ω is the frequency of oscillation. Substituting equation (62) into equation (61), we get e approximate frequency of oscillation can be obtained by setting [59] z zV 0 From equations (63) and (64), the amplitude-frequency relationship of equation (56) can be obtained as follows:

Advances in Civil Engineering
Substituting the shape function in equation (55) into the expressions of coefficients c 1 and c 2 in equations (57) and (58) and calculating the integrals, the amplitude-frequency relationship for the FG nanobeam with pinned-pinned boundary conditions is obtained as Equation (66) is the nondimensional nonlinear frequency of the FG nanobeam considering the thickness effect. e nondimensional linear frequency of the FG nanobeam can be obtained from equation (66) by letting V 0 � 0 as follows: (67)

Results and Discussion
In this section, various numerical examples are presented and discussed to verify the accuracy of the present results in predicting the free nonlinear vibration responses of the FG nanobeams. A nanobeam composed of aluminum (metal) and alumina (ceramic) is considered. e material properties of aluminum and alumina are shown in Table 1 [50]. e variation of the physical neutral surface to the powerlaw index k is presented in Table 2. Table 3 shows the first nondimensional natural frequencies (linear frequencies) of the FG classical beams for different values of the power-law index k and the slenderness ratio (L/h). It can be observed that the results obtained for the Euler-Bernoulli beam theory are in good agreement with the results given by ai and Vo [50].
From equation (67), we can see that the nondimensional natural frequency of the FG nanobeam depends not only on the material length scale parameter (β � l m /L) and the nonlocal parameter (α � ea/L) but also on the slenderness ratio (c � L/h). e material length scale parameter increases the natural frequency of the FG nanobeam, while the nonlocal parameter decreases the natural frequency of the FG nanobeam. For two cases of α � 0.01, β � 0.01 and α � 0.02, β � 0.03, Tables 4 and 5 show the first nondimensional natural frequencies of the FG nanobeams for different values of the power-law index k and the slenderness ratio (L/h), respectively. We can see that the slenderness ratio L/h has a great influence on the natural frequency of the FG nanobeam, and the slenderness ratio increases the natural frequency of the FG nanobeam. On the contrary, the Table 1: e material properties of aluminum and alumina [50].

Materials
Young's modulus E (GPa)     power-law index k decreases the natural frequency of the FG nanobeam. With each specific value of the power-law index and the slenderness ratio, the natural frequency of the FG nanobeam considering the thickness effect is always greater than the natural frequency of the FG nanobeam without considering the thickness effect.
To investigate effects of the material length scale parameter β � l m /L and the nonlocal parameter α � ea/L on the vibration responses of the FG nanobeam, we introduce the scale ratio c and the frequency ratio ω ratio , respectively, which are defined as  where ω L0 is the linear frequency of the homogeneous beam using the classical elasticity theory.
Effects of the material length scale parameter (β � l m /L) and the nonlocal parameter (α � ea/L) on the nonlinear vibration response of the FG nanobeam are presented in Figures 2-7. Figures 2-5 show the variations of the nonlinear frequencies and the frequency ratios of the FG nanobeams to the nonlocal parameter and the material length scale parameter. It can be concluded that both nonlinear frequency and frequency ratio of the nanobeam decrease as the nonlocal parameter increases and increase as the material length scale parameter increases. e effects of both nonlocal and material length scale parameters are integrated in the NSGT; thus, to have a better observation of the effects of these two parameters on the vibration response of the nanobeam, variations of the  nonlinear frequencies and frequency ratios of the nanobeam to the scale ratio c are investigated. Figures 6 and 7 show the variations of the nonlinear frequencies and the frequency ratios of the FG nanobeams, respectively, to the scale ratio c for some different values of the material length scale parameter. As expected, when the thickness effect is not considered (Figures 6(a) and 7(a)), when c � 1 (i.e., ea � l m ), the nonlinear frequencies and frequency ratios of the FG classical beams are equal to the ones of the FG nonlocal strain gradient beams; the FG nanobeam shows a stiffnesssoftening effect if c > 1 (i.e., ea > l m ) and shows a stiffnesshardening effect if c < 1 (i.e., ea < l m ). ese results completely agree with the results obtained by Li and Hu [52]. As in Figures 6(a) and 7(a), without considering the thickness effect, the nonlinear frequencies and frequency ratios of the FG nanobeams increase as the material length scale parameter (β � l m /L) increases if c < 1, and the nonlinear frequencies and frequency ratios of the FG nanobeams  decrease as the material length scale parameter (β � l m /L) increases if c > 1. However, when the thickness effect is considered as shown in Figures 6(b) and 7(b), we see that the slenderness ratio L/h has a great influence on the size-dependent vibration response of the FG nanobeam. From  Figures 6(b) and 7(b), it can be observed that, for a selected specific value of the power-law index k � 2, the slenderness ratio L/h � 15 and the initial amplitude V 0 � 0.01. e FG nanobeam shows a stiffness-softening effect if c > 18 (i.e., ea > 18l m ) and a stiffness-hardening effect if c < 18 (i.e., ea < 18l m ); the nonlinear frequencies and frequency ratios of the FG nonlocal strain gradient beams are equal to the ones of the FG classical beams if ea � 18l m . From Figure 7, we see that the frequency ratio of the FG nanobeam considering the thickness effect is always larger than the frequency ratio of the FG nanobeam without considering the thickness effect. e frequency ratio of the FG nanobeam decreases when the scale ratio increases.
Effect of the slenderness ratio L/h on the nonlinear vibration response of the FG nanobeams is also investigated Nonlinear frequency, ω  ese figures show the variations of the nonlinear frequencies and frequency ratios of the FG nanobeams to the scale ratio c for some differential values of the slenderness ratio L/h. It can be concluded that the slenderness ratio L/h leads to an increase in the nonlinear frequencies and the frequency ratios of the FG nanobeams.
However, with the specific value of the slenderness ratio, both nonlinear frequencies and frequency ratios of the FG nanobeams decrease when the scale ratio increases.
Finally, the effect of the power-law index k on the nonlinear vibration response of the FG nanobeams is examined. Figures 10 and 11 show the variations of the  nonlinear frequencies and frequency ratios of the FG nanobeams to the scale ratio c for some various values of the power-law index k, respectively, while Figures 12 and 13 show the variations of the nonlinear frequencies and frequency ratios of the FG nanobeams to the slenderness ratio L/h for some various values of the power-law index k, respectively. It can be observed that both nonlinear frequencies and the frequency ratios of the FG nanobeams decrease when the power-law index k increases. However, for each specific value of the power-law index k, the nonlinear frequencies and the frequency ratios of the FG nanobeam decrease as the scale ratio c increases and increase as the slenderness ratio L/h increases. e nonlinear frequencies of the FG nanobeam considering the thickness effect are always larger than the ones without considering the thickness effect.

Conclusion
A nonlocal strain gradient Euler-Bernoulli beam model considering the effect of the thickness is developed to analyze the size-dependent vibration response of the functionally graded nanobeam. With the help of Galerkin's technique and the Hamiltonian approach, the amplitude-frequency relationship of the pinned-pinned FG nanobeam is obtained in a closed form.
When considering the thickness effect, the size-dependent vibration response of the nanobeam depends not only on the scale ratio c � ea/l m (the ratio of nonlocal parameter ea to material length scale parameter l m ) but also on the slenderness ratio L/h. e thickness affects greatly the sizedependent vibration response of the FG nanobeam.
Effects of the slenderness ratio L/h and the power-law index k on the nonlinear vibration response of the FG nanobeam are also investigated. It is found that the slenderness ratio leads to an increase in the linear and nonlinear frequencies of the FG nanobeams, while the power-law index k leads to a decrease in the linear and nonlinear frequencies of the FG nanobeams.
When the thickness effect is not considered, the nonlinear frequencies and frequency ratios of the FG nanobeams increase as the material length scale parameter (β � l m /L) increases if ea < l m , and the nonlinear frequencies and frequency ratios of the FG nanobeams decrease as the material length scale parameter (β � l m /L) increases if ea > l m . However, when the thickness effect is considered, the slenderness ratio L/h has a great influence on the size-dependent vibration response of the FG nanobeam. For a selected specific value of the power-law index k � 2, the slenderness ratio L/h � 15 and the initial amplitude V 0 � 0.01. e FG nanobeam shows a stiffness-softening effect if c > 18 (i.e., ea > 18l m ) and a stiffness-hardening effect if c < 18 (i.e., ea < 18l m ); the nonlinear frequencies and frequency ratios of the FG nonlocal strain gradient beams are equal to the ones of the FG classical beams if ea � 18l m .

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.