Vibration Analysis of a Unimorph Nanobeam with a Dielectric Layer of Both Flexoelectricity and Piezoelectricity

In this study, for the first time, free and forced vibrational responses of a unimorph nanobeam consisting of a functionally graded base, along with a dielectric layer of both piezoelectricity and flexoelectricity, is investigated based on paradox-free local/nonlocal elasticity. The formulation and boundary conditions are attained by utilizing the energy method Hamilton’s principle. In order to set a comparison, the formulation of a model in the framework of differential nonlocal is first presented. An effective implementation of the generalized differential quadrature method (GDQM) is then utilized to solve higher-order partial differential equations. This method can be utilized to solve the complex equations whose analytic results are quite difficult to obtain. Lastly, the impact of various parameters is studied to characterize the vibrational behavior of the system. Additionally, the major impact of flexoelectricity compared to piezoelectricity on a small scale is exhibited. The results show that small-scale flexoelectricity, rather than piezoelectricity, is dominant in electromechanical coupling. One of the results that can be mentioned is that the beams with higher nonlocality have the higher voltage and displacement under the same excitation amplitude. The findings can be helpful for further theoretical as well as experimental studies in which dielectric material is used in smart structures.


Introduction
The big data era demands developments in micro/nano-electromechanical systems (M/NEMS) including nanosensors [1,2], nanotransistors [3], nanoharvesters [4,5], and nanoelectronics. The electromechanical coupling in dielectric materials plays a crucial role in ascertaining the performance of many M/NEMS devices. The electromechanical coupling was mainly related to piezoelectricity resulting from a uniform strain of the dielectric materials. However, recently, the coupling effect due to the nonuniform strain, i.e., strain gradient, has become increasingly attractive to researchers in exploring the potential of small-scale smart structures, since the gradient of the strain becomes more predominant with reduced sizes, especially at a sub-micron and nanoscale [6,7]. Thus, to capture the size effect, as the classic continuum theories are insufficient, a large amount of effort has been devoted to enriching the continuum theories by introducing an additional material length scale. The integral nonlocal formulation accounts for the size effect [8,9], which assumes that the strain in a particular point of a structure results not only from local forces applied at that point but also from other forces applied to other regions of the domain. In the effort of enhancing the conventional elastic theory, the differential form of nonlocal theory [10], largely due to less complexity together with a reduction in computational costs compared to the integral form, has been proposed. Based on the differential form of nonlocal elasticity, in sensors and actuators. Additionally, using the strain gradient sensing of flexoelectric materials, a sensor for detecting crack growth was presented [7]. Further, the enhancement impact of considering the flexoelectricity-resulting from nonuniform strain-on the energy harvesting of piezo as well as non-piezo materials was investigated [40]. Yan and Jiang [41] presented a study to show the effect of flexoelectricity on the bending of nanobeams under electrical as well as mechanical loading. In another paper, by using Timoshenko beam theory, they explored the flexoelectricity impacts on the dynamic and static responses of simply supported nanobeams [42]. In addition, the nonlocal theory has been used to capture the size effect. For instance, Sidhardh and Ray [43], by employing the finite element method, examined the static bending of a two-layered nanobeam, including a layer with a flexoelectric effect as an actuator. In their report, elasticity was utilized as the size-dependent theory. Additionally, based on strain gradient theory and isotropic flexoelectric theory, the vibration analysis of microplates considering the microscopic electrical field, polarization gradient, and strain gradient effects were examined [44]. In these studies, the extensive flexoelectric characteristics on the system's dynamic and static responses was exhibited.
In this work, the free and forced vibration, as well as the output voltage of a unimorph nanobeam including a functionally graded (FG) base along with a piezo-flexoelectric dielectric layer, is modeled for the first time, according to paradox-free elasticity or two-phase local/nonlocal theory. The formulation and boundary conditions are extracted by utilizing the energy method, i.e., Hamilton's principle. The nonlocal differential elastic beam formulations are presented and the GDQM is employed. The performance of the proposed model is presented by studying the influence of the input material parameters, boundary conditions, and structural dimensions. The obtained numerical results indicate the possibility of our proposed approach in characterizing the vibrational behavior of the piezo-flexoelectric bimorph nanobeam.

Problem Formulation
The vibration of a unimorph nanobeam energy harvesting including a dielectric layer with the flexoelectricity and piezoelectricity as well as a functionally graded (FG) base is examined. The Euler beam theory as well as two-phase local/nonlocal elasticity are employed. Figure 1 exhibits the schematic view for the nanobeam with L, b, h, and h d as its length, width, the thickness of FG base, and thickness of the dielectric layer, respectively. Open electric circuit condition is applied, and the voltage difference between the upper and lower dielectric layers is measured. strain gradient effect of flexoelectric material for sensing and energy harvesting [37,38]. Jiang et al. [39] reviewed the flexoelectricity in various materials as well as the application of flexoelectricity in sensors and actuators. Additionally, using the strain gradient sensing of flexoelectric materials, a sensor for detecting crack growth was presented [7]. Further, the enhancement impact of considering the flexoelectricity-resulting from nonuniform strain-on the energy harvesting of piezo as well as non-piezo materials was investigated [40]. Yan and Jiang [41] presented a study to show the effect of flexoelectricity on the bending of nanobeams under electrical as well as mechanical loading. In another paper, by using Timoshenko beam theory, they explored the flexoelectricity impacts on the dynamic and static responses of simply supported nanobeams [42]. In addition, the nonlocal theory has been used to capture the size effect. For instance, Sidhardh and Ray [43], by employing the finite element method, examined the static bending of a two-layered nanobeam, including a layer with a flexoelectric effect as an actuator. In their report, elasticity was utilized as the size-dependent theory. Additionally, based on strain gradient theory and isotropic flexoelectric theory, the vibration analysis of microplates considering the microscopic electrical field, polarization gradient, and strain gradient effects were examined [44]. In these studies, the extensive flexoelectric characteristics on the system's dynamic and static responses was exhibited.
In this work, the free and forced vibration, as well as the output voltage of a unimorph nanobeam including a functionally graded (FG) base along with a piezo-flexoelectric dielectric layer, is modeled for the first time, according to paradox-free elasticity or two-phase local/nonlocal theory. The formulation and boundary conditions are extracted by utilizing the energy method, i.e., Hamilton's principle. The nonlocal differential elastic beam formulations are presented and the GDQM is employed. The performance of the proposed model is presented by studying the influence of the input material parameters, boundary conditions, and structural dimensions. The obtained numerical results indicate the possibility of our proposed approach in characterizing the vibrational behavior of the piezo-flexoelectric bimorph nanobeam.

Problem Formulation
The vibration of a unimorph nanobeam energy harvesting including a dielectric layer with the flexoelectricity and piezoelectricity as well as a functionally graded (FG) base is examined. The Euler beam theory as well as two-phase local/nonlocal elasticity are employed. Figure 1 exhibits the schematic view for the nanobeam with L, b, h, and hd as its length, width, the thickness of FG base, and thickness of the dielectric layer, respectively. Open electric circuit condition is applied, and the voltage difference between the upper and lower dielectric layers is measured.  Here, the electrical enthalpy for the dielectric layer with both the piezoelectricity and flexoelectricity is as follows [45] in which c ijkl denotes the elastic constant and the dielectric constant tensor is κ kl . Additionally, µ ijkl , b ijkl as well as e kij are flexoelectric constants, nonlocal electrical coupling, and piezoelectric constants, respectively. In addition, E l and E k are electrical fields, and ε ij represent the strain tensor of the beam. Now, according to Euler-Bernoulli theory, the displacement field for the neutral axis of the nanobeam is as follows [26] in which U x and U z are displacement fields in x and z directions. Additionally, w(x, t) and ∂w(x,t) ∂x are transverse deflection and rotation of the neutral plane of the beam. Further, the corresponding strains of the nanobeam are [26] where ε xx and ε xx,z show the elastic strain and the gradient of the strain, respectively.
Here, the stress-strain relations associated with the functionally graded base is as follows [46] σ where E(z) is elastic modulus associated with the FG layer and σ b xx is the stress tensor related to the base layer. Next, the following equation represents the normal and shear stress associated with piezoelectric nanobeams with the flexoelectric effect considered [47].
where E z = − ∂φ ∂z and φ is electrical potential, and σ d xx as well as τ xxz represent the normal and second-order stress in the dielectric layer. Additionally, superscript d shows the dielectric layer. The dielectric nanobeam's electrical displacement and quadrupolar contribution are as follows [47].
in which D z and Q zz are electrical displacement and quadrupolar contribution. Now, with free electrical charges being zero, the following can be written based on Gauss's law [48].
Here, by using Equation (6) into Equation (7), the equation below can be yielded.
In addition, by considering the following boundary conditions for electrical potential Equation (8) can be solved for electrical potential as follows [48,49] where η = κ 33 b 33 based on Equation (10), E z and E z,z can be obtained as follows.
Additionally, the functionally graded material properties used in the base of the beam are in which ρ cm and ρ m denote the mass density related to metal and ceramic in the FG layer, respectively, and E cm and E m represent the ceramic and metal elastic modulus, respectively. Additionally, ν m and ν cm are Poisson's ratio of metal and ceramic, and n represents the FG power index. Now, the strain Π s and kinetic Π k energy of the nanobeam made of two layers, including a dielectric layer, considering the flexoelectricity and piezoelectricity and an elastic base are shown in the following equation. where Here, the energy method is utilized to attain the equations in addition to geometrical end conditions.
The following equation, representing the governing equation associated with the vibration of a two-layered nanobeam that contains an FG base and a dielectric layer can be achieved using Equation (14).
Moreover, the end conditions are

Two-Phase Local/Nonlocal Theory
The stress-strain relation defined in the two-phase local/nonlocal theory is [26] t in which C, t(x), ε(x), x, α(x, x, κ), ξ, V, and k denote the fourth-order elasticity tensor, Cauchy stress tensor in the two-phase state, the strain tensor, reference point, kernel function, local phase fraction factor, domain volume, and nonlocal parameter, respectively. Thus, two-phase elasticity can be written as follows in which T(x, t) is a function of the local quadrupolar contribution. It should be mentioned that by setting ξ = 0 in Equation (18), the pure nonlocal theory can be attained. Now, by utilizing the transformation presented by Polyanin et al. [50], which was firstly used for two-phase elasticity by Fernandez et al. [26], the integral formulation is written in the differential format as follows Here, it is vital to mention that this transformation can be used only when two constitutive boundary conditions-CBCs-are satisfied. Therefore, the CBCs at the tips of the nanobeam are given.
Now, the stress-strain relations and the electrical displacements can be rewritten in two-phase form using Equation (18) [43,51]. Now-by using Equations (21)-(22)-the two-phase bending moment, two-phase higherorder bending moment in differential form are derived and presented in Appendices A.1 and A.2.

Solution Procedure
Generalized Differential Quadrature Method In this section, the solution procedure that is utilized to extract the forced and free vibrational response of a nanobeam which is made of two layers, including an FG base and a layer made of dielectric material considering flexoelectric and piezoelectric effects, is introduced. The GDQM from [33] is used. Based on this method, the momentum and higher-order momentum in addition to the deflection are the independent variables which can be rewritten using the separation of variables as follows.
Now, by employing the fourth-order GDQM, the n-th derivative of any function such as (x) can be written as follows.
In which Λ (r) jl (x i ) is Hermit interpolation, and ns denotes the number of sampling points. Additionally, the location corresponded to the grid points is x i . Here, the discrete form of the momentum, higher-order momentum, and displacement field can be given as the following equation together with Figure 2.
In which M j , P j , W j are the discretized form of bending moment, higher-order bending moment, and lateral displacement. In which j M , j P , j W are the discretized form of bending moment, higher-order bending moment, and lateral displacement Additionally, the distribution of grid points is selected by the following method.
( ) ( ) Now, the Hermit interpolation can be introduced as Additionally, the distribution of grid points is selected by the following method.
Now, the Hermit interpolation can be introduced as where where E j (x i ) are the Lagrange interpolation and the constant, a ri , b ri , and c ri , can be found in Appendix A.3. Additionally, the Lagrange interpolation can be written as where Additionally, the higher-order derivative associated with the Lagrange interpolation can be obtained in the following equation.
Now, the discretized formulation and boundary conditions presented in Appendix A.4 can be written in a matrix format as follows.
in which b and d show boundary and domain grid points, respectively. Additionally, the discretized form of formulation associated with the differential form of the purely nonlocal model is presented in Appendix A.5. Next, the forced vibration analysis results are obtained using Newmark-beta and GIQM procedure, as explained in Ref. [1].
Additionally, the damping matrix is where where ξ i is the damping ratio and ω i as well as ω j are the natural frequencies.
In addition, the harmonic force applied to the system for the forced vibrational results can be defined as follows.
F ext = f 0 sin(Ωt) (34) in which f 0 and Ω denote the amplitude and the frequency of the force.

Numerical Results and Discussion
This section presents the numerical results related to the vibration of a unimorph nanobeam, which is made of a dielectric layer with consideration of flexoelectric and piezoelectric effects as well as an FG layer based on two-phase elasticity. It should be mentioned that, in this study, it is supposed that the nanobeam is made of BaTiO 3 with material properties such as c 11 = 167.55 (GPa), ρ = 6020 kg/m 3 , e 31 = 4.43 C/m 2 , b 33 = 1.265 × 10 −17 Jm 2 /VC , κ 33 = 1.265 × 10 −8 (C/Vm), and µ 1133 = 5 × 10 −8 (C/m) [52]. Additionally, the mechanical properties related to the FG layer are ν m = 0.3, ν cm = 0.2, E m = 70(GPa), E cm = 200(GPa), ρ m = 2702 kg/m 3 , and ρ cm = 5700 kg/m 3 . Additionally, the damping ratio is ξ i = 0.05 and b = h + h d . It should be mentioned that from [52], the value of flexoelectric coefficient is 1-10 V, in this work we choose f 31 = 3.9526 V, and the flexoelectric coefficient is computed as Firstly, in Table 1, to verify the present formulation in the two-phase framework, the first vibrational frequency ratio of a nanobeam with flexoelectric effect is obtained by eliminating the FG base, and compared the results to those from [42]. In this table, the other constants are ξ = 1 and L = 40 h d . Additionally, the frequency ratio in the reference is defined as follows.  Next, to study the validity of the two-phase theory formulation, in Table 2, the dimensionless frequencies of a nanobeam with omitting the effect of piezoelectricity and flexoelectricity are tabulated and compared with those of Ref. [26]. The other constant which plays a role in determining the values is k = 0.05 L. The results are presented for simply supported (SS) and clamped-free (CF) end conditions. It should be noted that the exact results are obtained using the equations and solution procedure presented in the reference. However, the result which is entitled from Ref. [26] is extracted approximately from the figure presented in this article. As can be seen from these two tables, the close agreement between the results acquired through the presented formulation and solution method confirms the credibility of these two in exploring the vibration behavior of a two-layered nanobeam, including an FG base and a dielectric layer with the flexoelectricity and piezoelectricity. In addition, it should be mentioned that two-phase theory can have good agreement with the results of molecular dynamic models [25], making it a good tool to investigate small-sized structures. Additionally, by setting ξ = 1 the results associated with the classic continuum are obtained, showing the validity of the current continuum theory.
Here, the frequency ratio, which is defined below, for different nanobeam lengths and boundary conditions are tabulated in Table 3. The other constants used in this table are k = 0.05 L, ξ = 0.1, n = 0, and h d = 0.005 L .
Frequency of the nanobeam with flexoelectricity Frequency of the nanobeam without flexoelectricity (36)  It can be understood from this table that the impact of flexoelectricity on the vibrational frequency of a unimorph system can be intensified by reducing the length of the nanobeam, so much so that the highest frequency ratio, regardless of the boundary condition, occurs in the cases with the lowest length. Now, the impact of slenderness related to the FG base of the two-layered nanobeam with an FG layer and a dielectric layer on its frequency ratio is examined for various FG index numbers in Figure 3. Additionally, other constants are k = 0.05 L, ξ = 0.1, and h d = 0.005 L . As can be seen from these two tables, the close agreement between the results acquired through the presented formulation and solution method confirms the credibility of these two in exploring the vibration behavior of a two-layered nanobeam, including an FG base and a dielectric layer with the flexoelectricity and piezoelectricity. In addition, it should be mentioned that two-phase theory can have good agreement with the results of molecular dynamic models [25], making it a good tool to investigate small-sized structures. Additionally, by setting 1 ξ = the results associated with the classic continuum are obtained, showing the validity of the current continuum theory.
Here, the frequency ratio, which is defined below, for different nanobeam lengths and boundary conditions are tabulated in Table 3  It can be understood from this table that the impact of flexoelectricity on the vibrational frequency of a unimorph system can be intensified by reducing the length of the nanobeam, so much so that the highest frequency ratio, regardless of the boundary condition, occurs in the cases with the lowest length. Now, the impact of slenderness related to the FG base of the two-layered nanobeam with an FG layer and a dielectric layer on its frequency ratio is examined for various FG index numbers in Figure 3 Figure 3 exhibits that intensifying the base layer's thickness ratio diminishes the value of the frequency ratio of the nanobeam. In other words, the thicker base causes the effect of flexoelectricity to reduce. Additionally, it is worth noting that the cases with a higher FG index can possess a higher frequency ratio, which means that the flexoelectricity is higher if the FG power index is higher. Further, the higher the FG power index is, the higher the metal phase's contribution in the base is. It can be seen that the frequency ratio is dominated by the base plate and thickness ratio of the beam rather than end conditions. Next, the impact of considering flexoelectricity in the current model, specifically smallscale, on the vibration response of the system, is shown in Figure 4. In this figure, the nondimensional vibration frequency of the beam is plotted against the value of the thickness ratio of the base for three models which consider flexoelectricity and piezoelectricity, only flexoelectricity, and only piezoelectricity. Additionally, the non-dimensional vibration frequency is obtained through the following equation.   Similar to the previous figure, it can be seen that the effect of flexoelectricity rises with reducing the thickness of the nanobeam. Additionally, it can be concluded that, as the model with flexoelectricity and piezoelectricity yields a similar frequency to the one that only considers flexoelectricity, and hence, flexoelectricity is considered the dominant electromechanical coupling on small scales.
Here, Figure 5 deals with the impact of the thickness ratio of the dielectric layer with respect to the base layer's thickness on the vibrational frequency of the system. Additionally, other constants in this figure are k = 0.05 L, ξ = 0.1, and h = 0.01 L . The quite interesting results indicate that increasing the dielectric layer's thickness causes the system's vibrational frequency to increase with a high slope when the thickness ratio is far from 0.5. By increasing the value, the increment rate diminishes to a point in which the vibration frequency slightly decreases. Now, the effect of two-phase elasticity parameters is investigated on the non-dimensional vibration of a two-layered nanobeam, including an FG base and a dielectric layer.
To this aim, in Figures 6 and 7, respectively, the influence of ξ and / k L on the fundamental vibration frequency of the two-layered system was analyzed for different end conditions and values for the FG index. Additionally, the other parameters that can affect the results are ( ) Additionally, in Figure 6, the arrowheads show the dimensionless vibration frequency obtained by the differential form of purely nonlocal. Further, in Figure 8, the non-dimensional frequencies are obtained according to the differential form of purely nonlocal and are plotted against various nonlocality.
ω ω Now, the effect of two-phase elasticity parameters is investigated on the non-dimensional vibration of a two-layered nanobeam, including an FG base and a dielectric layer. To this aim, in Figures 6 and 7, respectively, the influence of ξ and k/L on the fundamental vibration frequency of the two-layered system was analyzed for different end conditions and values for the FG index. Additionally, the other parameters that can affect the results are L = 50 (nm), h d = 0.005 L , and h = 0.01 L . Additionally, in Figure 6, the arrowheads show the dimensionless vibration frequency obtained by the differential form of purely nonlocal. Further, in Figure 8, the non-dimensional frequencies are obtained according to the differential form of purely nonlocal and are plotted against various nonlocality.
In addition, it can be understood that escalating ξ can cause the vibrational frequency of the nanobeam to increase, which is more observable in C-C and C-F end conditions. This means that the frequency of the unimorph is lower in the cases that nonlocal forces have bigger impacts and lower values of ξ. Additionally, as previously stated, the nanobeam in which the FG base has a higher index value can have a higher frequency, despite the boundary conditions and ξ values. FG index shows the function based on which the FG between two surface layers of ceramic and metal forms. Furthermore, it should be mentioned that, according to the purely nonlocal model, by decreasing the value of ξ toward 0, the vibration frequency of the system in the two-phase framework should lead to the values obtained by the differential form of purely nonlocal. Thus, as it can be seen in this figure, except for the cases with SS end conditions, the results of the two-phase model are significantly different from those attained via differential form of purely nonlocal elasticity, showing the inconsistency of this model due to the lack of additional end condition in the transformation from integral to differential model. In addition, it should be mentioned that studies have shown that not all of the forces between the particle of small-scale structures are purely nonlocal, and it is more a combination of local and nonlocal forces [53]. Thus, the two-phase theory is a more reliable mathematical model in order to investigate small-scale structures. Further, the value of ξ depends on the various geometrical parameters of the structure, which should be determined in every case specifically.
(a) (b) Figure 5. Variation of fundamental vibration frequency of two-layered nanobeams against hd/h for various n and (a) C-C and (b) C-F boundary conditions. Now, the effect of two-phase elasticity parameters is investigated on the non-dimensional vibration of a two-layered nanobeam, including an FG base and a dielectric layer.
To this aim, in Figures 6 and 7, respectively, the influence of ξ and / k L on the fundamental vibration frequency of the two-layered system was analyzed for different end conditions and values for the FG index. Additionally, the other parameters that can affect the results are ( ) 50 Additionally, in Figure 6, the arrowheads show the dimensionless vibration frequency obtained by the differential form of purely nonlocal. Further, in Figure 8, the non-dimensional frequencies are obtained according to the differential form of purely nonlocal and are plotted against various nonlocality. In addition, it can be understood that escalating ξ can cause the vibrational frequency of the nanobeam to increase, which is more observable in C-C and C-F end conditions. This means that the frequency of the unimorph is lower in the cases that nonlocal Interesting results of Figure 7 show the importance of using two-phase elasticity in order to study the small-scale structures with consideration to the piezoelectric and flexoelectric coupling, as by increasing the nonlocality of nanobeam, the value of frequency diminishes. In other words, by comparing the results of Figures 7 and 8, decreasing frequency shows the softening effect which the nonlocal integral elasticity proposed and differential form of purely nonlocal theory is insufficient to capture, specifically in the C-F end condition. Additionally, similar to the previous figure, this softening effect can be more observable in C-F, C-C, and C-S end conditions compared to S-S one. The only end condition for which differential form of purely nonlocal can produce a comparatively accurate result is S-S, as the effect of additional boundary conditions is the least in this condition. Interesting results of Figure 7 show the importance of using two-phase elasticity in order to study the small-scale structures with consideration to the piezoelectric and flexoelectric coupling, as by increasing the nonlocality of nanobeam, the value of frequency diminishes. In other words, by comparing the results of Figures 7 and 8, decreasing frequency shows the softening effect which the nonlocal integral elasticity proposed and differential form of purely nonlocal theory is insufficient to capture, specifically in the C-F end condition. Additionally, similar to the previous figure, this softening effect can be more observable in C-F, C-C, and C-S end conditions compared to S-S one. The only end condition for which differential form of purely nonlocal can produce a comparatively accurate result is S-S, as the effect of additional boundary conditions is the least in this condition. Here, in Figure 9, the impact of the flexoelectric coefficient on the vibrational frequencies of a unimorph system containing an FG base in addition to a dielectric layer is examined. In these figures, the variation of vibration frequencies Here, in Figure 9, the impact of the flexoelectric coefficient on the vibrational frequencies of a unimorph system containing an FG base in addition to a dielectric layer is Here, in Figure 9, the impact of the flexoelectric coefficient on the vibrational frequencies of a unimorph system containing an FG base in addition to a dielectric layer is examined. In these figures, the variation of vibration frequencies The results from this figure indicate that the vibrational frequency of the unimorph system is higher, provided that the flexoelectric coefficient associated with the dielectric layer possesses a higher positive or negative value. Additionally, the softening effect due to nonlocality can be observed in this figure, as the vibrational frequency is lower, despite the value of 31 μ , for beams with higher nonlocality. Now, by putting a harmonic load with an excitation frequency of 0.85f1 to 1.15 f1-in which f1 represents the fundamental vibration frequency of the system-the displacement response and output voltage of the system are obtained. Additionally,  The results from this figure indicate that the vibrational frequency of the unimorph system is higher, provided that the flexoelectric coefficient associated with the dielectric layer possesses a higher positive or negative value. Additionally, the softening effect due to nonlocality can be observed in this figure, as the vibrational frequency is lower, despite the value of µ 31 , for beams with higher nonlocality. Now, by putting a harmonic load with an excitation frequency of 0.85 f 1 to 1.15 f 1 -in which f 1 represents the fundamental vibration frequency of the system-the displacement response and output voltage of the system are obtained. Additionally, f 0 = 2 × 10 −9 N m in a condition that L = 50 (nm), h d = 0.005 L , and h = 0.01 L . Additionally, it should be noted that the deflection FRFs are presented for the endpoint in CF and the middle point for the CC boundary conditions. Next, the impact of different parameters on the frequency response function (FRF) associated with nanobeam's deflection and output voltage are investigated.
Firstly, the deflection FRF for two types of boundary conditions for three different models are presented in Figure 10. The three models used are the beams considering only piezoelectric coupling, flexoelectric coupling, and a combination of both electromechanical couplings. It can be seen that flexoelectricity is the dominant form of electromechanical coupling, as the result related to the case that considers both of the couplings has a similar outcome to the one only considering the flexoelectricity. Also, as previously shown, the vibrational frequency of the system with flexoelectricity has a higher resonant frequency than the one in which piezoelectricity is the only considered form of electromechanical coupling. Therefore, these results indicate that in order to reach a better model for dielectric materials on a small scale, flexoelectric coupling must be considered the dominant form of electromechanical coupling. Additionally, in Figures 11 and 12, respectively, the impact of the base's slenderness on the displacement and voltage FRF of the unimorph system with two boundary conditions is examined. The interesting results in these figures indicate that intensifying h/L of the nanobeam can cause the deflection of the beam in addition to the output voltage to diminish, regardless of the boundary condition. However, the deflection and output voltage are higher provided that the softer end condition is utilized. The other notable result is that the output voltage and displacement peak occur in higher frequencies by decreasing h/L. Additionally, in Figures 11 and 12, respectively, the impact of the base's slenderness on the displacement and voltage FRF of the unimorph system with two boundary conditions is examined. The interesting results in these figures indicate that intensifying h/L of the nanobeam can cause the deflection of the beam in addition to the output voltage to diminish, regardless of the boundary condition. However, the deflection and output voltage are higher provided that the softer end condition is utilized. The other notable result is that the output voltage and displacement peak occur in higher frequencies by decreasing h/L.
Next, in Figures 13 and 14, the displacement and voltage FRFs are plotted for different nonlocality and boundary conditions. In addition, other constants are equal to ξ = 0.1, n = 1, L = 50 (nm), h d = 0.005 L , and h = 0.01 L . It can be seen that the peaks related to deflection FRF and voltage FRF are moving towards lower frequencies due to softening effect of nonlocality. Therefore, as the cases with a higher nonlocality are softer, their deflection corresponding to these cases is higher, meaning they possess higher energy and can produce higher voltages. Similar to the previous figure, the output voltage and deflection are higher in softer boundary conditions. Lastly, the displacement and voltage FRFs for different ξ are plotted in Figures 15 and 16. In addition, the other constants are as follows k/L = 0.05, n = 1, L = 50 (nm), h d = 0.005 L , and h = 0.01 L . As expected, increasing the value leads to increasing the system's natural frequency, moving the deflection and voltage peak in FRFs towards higher frequencies.
Consequently, as the system is less nonlocal by intensifying ξ, the system is stiffer and produces less voltage, regardless of the boundary conditions. Additionally, in Figures 11 and 12, respectively, the impact of the base's slenderness on the displacement and voltage FRF of the unimorph system with two boundary conditions is examined. The interesting results in these figures indicate that intensifying h/L of the nanobeam can cause the deflection of the beam in addition to the output voltage to diminish, regardless of the boundary condition. However, the deflection and output voltage are higher provided that the softer end condition is utilized. The other notable result is that the output voltage and displacement peak occur in higher frequencies by decreasing h/L. It can be seen that the peaks related to deflection FRF and voltage FRF are moving towards lower frequencies due to softening effect of nonlocality. Therefore, as the cases with a higher nonlocality are softer, their deflection corresponding to these cases is higher, meaning they possess higher energy and can produce higher voltages. Similar to the previous figure, the output voltage and deflection are higher in softer boundary conditions.
(a) (b) It can be seen that the peaks related to deflection FRF and voltage FRF are moving towards lower frequencies due to softening effect of nonlocality. Therefore, as the cases with a higher nonlocality are softer, their deflection corresponding to these cases is higher, meaning they possess higher energy and can produce higher voltages. Similar to the previous figure, the output voltage and deflection are higher in softer boundary conditions. As expected, increasing the value leads to increasing the system's natural frequency, moving the deflection and voltage peak in FRFs towards higher frequencies. Consequently, as the system is less nonlocal by intensifying ξ , the system is stiffer and produces less voltage, regardless of the boundary conditions.  As expected, increasing the value leads to increasing the system's natural frequency, moving the deflection and voltage peak in FRFs towards higher frequencies. Consequently, as the system is less nonlocal by intensifying ξ , the system is stiffer and produces less voltage, regardless of the boundary conditions.

Conclusions
This research presents a study on the frequency response of a nanobeam made of a FG base along with a dielectric layer considering both flexoelectric and piezoelectric couplings. The effect of the small size is considered by means of a paradox-free elasticity named two-phase local/nonlocal theory. Then, the equations of motions and end condi-

Conclusions
This research presents a study on the frequency response of a nanobeam made of a FG base along with a dielectric layer considering both flexoelectric and piezoelectric couplings. The effect of the small size is considered by means of a paradox-free elasticity named two-phase local/nonlocal theory. Then, the equations of motions and end conditions are extracted using the Euler-Bernoulli beam theory, variational energy method, and GDQM. The results of this study are validated by employing other published articles. Finally, a parametric investigation is presented, the most important results of which are: The smaller the nanobeams are, the dominant effect of flexoelectricity over piezoelctricty can be observed more clearly.
The lower values of ξ the system are closer to being purely nonlocal, causing the frequency to diminish. On the other hand, increasing ξ and h/L leads to a reduction in the peak values of voltage and displacement FRFs.
By intensifying the FG index, the vibrational frequency as well as the flexoelectric effect of the dielectric layer can be increased.
The higher the nonlocality, the higher the voltage and displacement FRFs peaks. The two phase instead of differential nonlocal possess no paradoxes, making a reliable model by which the behavior of MEMS and NEMS can be studied. These findings suggest promising applications in nanoenergy transduction, nanogenerator, nanosensing and nanoactuation.  The ceramic and metal elastic modulus ν m and ν cm Poisson's ratio of metal and ceramic N FG power index Π s and Π k Strain and kinetic energy of the nanobeam C, t(x), ε(x), x, α(x, x, κ), The fourth-order elasticity tensor, Cauchy stress ξ, V, and k tensor in the two-phase state, the strain tensor, reference point, kernel function, local phase fraction factor, domain volume, as well as nonlocal parameter, respectively Λ (r) Hermit interpolation x i The location corresponded to the grid points M j , P j , W j Discretized form of bending moment, higher-order bending moment, and lateral displacement E j (x i ) Lagrange interpolation f 0 and Ω The amplitude and the frequency of the force Equation (20). should be employed in order to obtain the differential form of the twophase bending moment, the first part of Equation (22). Therefore, the following is proposed.
Now, the second part of Equation (A1) is rewritten as follows.
T(x, t) = EI e f f  (20), the differential form of the two-phase bending moment can be attained.