#### Abstract

A fully gradient elasticity model for bending of nanobeams is proposed by using a nonlocal thermodynamic approach. As a basic theoretical novelty, the proposed constitutive law is assumed to depend on the axial strain gradient, while existing gradient elasticity formulations for nanobeams contemplate only the derivative of the axial strain with respect to the axis of the structure. Variational equations governing the elastic equilibrium problem of bending of a fully gradient nanobeam and the corresponding differential and boundary conditions are thus provided. Analytical solutions for a nanocantilever are given and the results are compared with those predicted by other theories. As a relevant implication of applicative interest in the research field of nanobeams used in nanoelectromechanical systems (NEMS), it is shown that displacements obtained by the present model are quite different from those predicted by the known gradient elasticity treatments.

#### 1. Introduction

Nanostructures are nowadays widely used as main components of many micro- and nanoelectromechanical systems (MEMS and NEMS) [1–5]. Typical MEMS and NEMS consist of thin beams made of metals, polymers, traditional silicon-based materials, or functionally graded materials [6–11].

The mechanical behaviour of nanostructures can be studied by both theoretical and experimental methods; see, for example, [12] for a review. Generally, conducting controlled experiments at the nanoscale is difficult and expensive. Accordingly, theoretical modelling is the main tool for studying the properties of nanostructures.

Nanobeams can be modelled by using an atomistic approach [13, 14] or continuum mechanics-based models which are extensively used due to computational simplicity and effectiveness. One of the basic assumptions in continuum models is that the lattice structure is neglected and it is replaced by a continuum medium.

In this framework, several theories have been introduced to consider small-scale effects at the micro- and nanoscale. Unlike classical (local) elastic models, nonlocal elastic models introduce length-scale parameters in the constitutive relations. On the basis of an analogy presented in [15, 16], nonlocal effects on nanorods and nanobeams, formulated according to the Eringen model, can be simulated by prescribing suitable axial and curvature distortions on corresponding local rods and beams. Accordingly a general procedure is provided to establish if nonlocal nanorods and nanobeams are free of small-scale effects.

Several models based on the Euler-Bernoulli or Timoshenko beam theories have been proposed in the literature in order to develop effective nonlocal nanobeam models for bending, buckling, and vibrations; see, for example, nonlocal elasticity model [17–24] couple stress theory [25, 26], modified couple stress theory [27–29], gradient plasticity for strain softening materials [30–33], gradient elasticity model [34–40], and FE analysis of finite random composite bodies [41, 42].

In particular, the existing gradient elasticity model (GM) for Euler-Bernoulli nanobeams encompasses the first derivative of the axial strain in the beam axis direction (see, e.g., the recent contribution by Aifantis and coworkers [40, 43]) and disregards the other nonvanishing component of the axial strain gradient.

The aim of this paper is to propose a fully gradient elasticity model (FUGM) for Euler-Bernoulli nanobeams, starting from a thermodynamic approach, in which all the components of the gradient of the axial strain are taken into account to model the nonlocal behaviour of the nanobeam. The solutions of the FUGM in terms of transverse displacements can thus be compared with the corresponding ones evaluated by the GM. As a result, the influence of the gradient component, which is usually disregarded in the GM, on the bending behaviour of nanobeams can be clearly enlightened.

The starting point of the present work is the definition of the Helmholtz free energy in terms of the axial strain and of its gradient. Accordingly a fully gradient elasticity model (FUGM) for bending of Euler-Bernoulli nanobeams is provided based on a nonlocal thermodynamic approach; see, for example, [44, 45]. Two length-scale parameters are introduced in the model and are associated with the nonvanishing components of the strain gradient of the axial strain.

Nonlocal thermodynamics allows us to build up a reliable methodology to derive the nonlocal variational formulation of the elastic equilibrium for the FUGM which yields the differential equation with the relevant boundary conditions.

As an example, a nanocantilever with a concentrated load at its tip is investigated using the FUGM. Results given by the FUGM are compared with those evaluated by the GM and it is shown that the effects of the gradient component are usually disregarded in the GM; that is, the derivative of the axial strain in the transverse direction is significant. The difference between the results of the FUGM and the results of the GM will be quantitatively shown and analysed. In particular, the nanobeam solved by the FUGM becomes stiffer than the one solved by the GM.

#### 2. Nonlocal Thermodynamics for Euler-Bernoulli Nanobeams

The internal energy density in a nonlocal elastic material can be assumed in the formwhere the kinematic internal variable is the gradient of the strain tensor and is the entropy.

Denoting by the Helmholtz free energy defined by means of the Legendre transform, the first law of thermodynamics for isothermal processes and for a nonlocal behaviour (see, e.g., [44, 46]) can be formulated as follows:where and is the nonlocal stress tensor field of the body . The superscript dot indicates the differentiation with respect to time and the symbol denotes the single (double) index saturation. The energy balance in (2) can be written pointwise in in the formwhere the thermodynamic scalar variable is the nonlocal residual function which accounts for the energy exchanges between neighbour particles (see, e.g., [47]). Since nonlocal effects due to elastic deformations are confined into the body, the residual fulfils the insulation condition:The second principle of nonlocal thermodynamics for isothermal processes is written in the local form everywhere in , where is the internal entropy production rate per unit volume; see, for example, [48]. Accordingly the nonlocal Clausius-Duhem inequality for isothermal processes follows from relation (3) in the formwhere the presence of the nonlocal residual function guarantees the nonnegativeness of the dissipation and accounts for material nonlocality.

The body energy dissipation is provided by integrating relation (5) to getExpanding relation (6), dissipation (5) is pointwise vanishing according to the reversible nature of the model. As a consequence, for any admissible deformation mechanism, it results inIn the next section, a fully gradient elasticity theory for Euler-Bernoulli nanobeams is presented starting from (7) and taking into account all the components of the strain gradient.

#### 3. Fully Gradient Elasticity Model for Nanobeams

Let us consider a homogeneous isotropic nanobeam of length . The -coordinate is taken along the length of the beam and the -coordinate along the thickness and the -coordinate is taken along the width of the beam. The geometry and the applied loads of the nanobeam are such that the displacements along the axes are functions of the - and -coordinates. It is further assumed that the displacement is identically zero. The cross-sectional area and the second moment of area about the -axis areThe proposed fully gradient elasticity model (FUGM) for Euler-Bernoulli nanobeams is based on the following classical displacement field:where is the transverse displacement and the apex denotes the -derivative of the function along the nanobeam axis . The rotation of the cross-section of the nanobeam is .

Accordingly the nonvanishing kinematically compatible deformation is given by the axial strain and the related nonvanishing components of the strain gradient arewhere denotes the bending curvature of the nanobeam. In the sequel, for simplicity, the subscript of strains will be dropped.

##### 3.1. Variational Formulation of Elastic Equilibrium

The proposed FUGM for Euler-Bernoulli nanobeams is governed by the following expression of the free energy:where is Young’s modulus. The coefficients and incorporate small-scale effects, where is a material constant and are the material length-scales associated with the variation of the normal strain in the longitudinal direction and in the transversal direction . Accordingly, the nonlocal model depends on two parameters and .

Note that the nonlocal model tends to the local model in the limit of vanishing nonlocal parameters as shown in the sequel.

Using expression (12) of the free energy, the thermodynamic requirement (7) provides the variational formulation:where the static variables and denote the axial stresses which are, respectively, duals of the strain and of the strain derivative . Thus, the axial stresses and are explicitly given bywhere and .

#### 4. Bending Solution

The explicit expression of the variational formulation for the FUGM can be recovered by (13). In fact substituting the kinematically compatible relation (10) and the expression of the strain derivatives (11) in (13) we get the nonlocal variational formulation associated with the FUGM: where the stress resultant moments areThe differential equilibrium equation and the boundary conditions corresponding to the proposed nonlocal FUGM can be obtained by applying the integration by parts to the nonlocal variational formulation (15) to getand the boundary conditions areThe classical differential equilibrium relation can be recovered by integrating by parts the l.h.s. of (15) to get , where is the distributed transversal load. The boundary conditions at provide the relations and , where is the shear force and are the transverse force and couple, respectively.

Differential equation (17) and boundary conditions (18) of the FUGM can be reformulated in terms of the transverse displacement by expressing the bending moments in the following forms:so that the differential equilibrium equation for the nonlocal FUGM can be obtained by substituting (19) into (17) and (18) to getand the related boundary conditions areThe bending moment can then be obtained in terms of the transverse displacement by performing an integration by parts of the second term at the r.h.s. of (15). Hence, we haveand the boundary condition at is fulfilled due to (18)_{3}.

The boundary conditions to be imposed at the end of a nanobeam to determine the six unknown coefficients associated with the FUGM are reported hereafter for external constraints of engineering interest:

(i) simply supported end at the point

(ii) clamped end at the point

(iii) free end at the point loaded by a couple and a force :

*Remark 1. *If the parameter is vanishing, the governing equation (20) and the boundary conditions (21) pertaining to the FUGM degenerate to the corresponding relations of the nonlocal Euler-Bernoulli nanobeam theory based on the gradient elasticity model (GM). Moreover, if both the material parameters and are equal to zero, the governing equations (20) and boundary conditions (21) reduce to those of the classical (local) Euler-Bernoulli beam model.

#### 5. Example: Cantilever Nanobeam

The solutions of the FUGM obtained in the previous sections are specialized for a nanocantilever subjected to a concentrated load at its tip in terms of the six unknown coefficients with introduced below. The geometric constants are the length , the width , and the height .

To this end it is convenient to introduce the following dimensionless parameter related to the length-scale parameter by the expression where is the radius of gyration.

The solution of the sixth-order differential equation (20) isand bending moment (22) is given byThe six unknown coefficients are evaluated by means of the boundary conditions provided by relations (24) at and by relations (25) at , where and .

Accordingly the six unknown coefficients of the FUGM appearing in the expression of transverse displacement (27) and of bending moment (28) are given byThe lower and upper bounds and of the nanocantilever transverse displacement can then be evaluated by taking the limit of for and , respectively, and are given byHence, the displacement of the FUGM belongs to the strip bounded by the functions and . Note that the lower bound coincides with the classical (local) displacement of the considered nanocantilever and the upper bound depends on , that is, the length-scale parameter . Moreover, the nonlocal displacement tends to vanish for .

The bending moment (28) of the considered statically determined nanobeam reduces to its classical (local) counterpart and the shear force is . A general theoretical motivation of this result is reported in [20].

*Remark 2. *The expression of the transverse displacement pertaining to the GM is obtained from (27) of the FUGM by setting (i.e., ) into the expressions of the six coefficients (29).

Let us now introduce the following dimensionless quantities:Thus, the parameter (see (26)) can be rewritten as and its plot is reported in Figure 1 in terms of the dimensionless ratio .

For computations, we assume the following values of the dimensionless parameters: and which yields .

The transverse dimensionless deflection of the nanocantilever is plotted in Figure 2 for with and for with .

Note that the transverse dimensionless deflection obtained by the FUGM(i)with and coincides with the classical (local) solution ;(ii)with and coincides with the one obtained by the GM with the corresponding ;(iii)for a given and for increasing values of the nanocantilever becomes stiffer than the solution obtained by the GM.

The dimensionless upper bound of the transverse deflection is independent of the length-scale parameters and is reported in Figure 2 with the black dot-dashed line. The dimensionless lower bound depends on the parameter and is reported in Figure 2 with a coloured dot-dashed line for the considered values of . Hence, for a given , the transverse dimensionless deflection of the FUGM belongs to the strip bounded by and by , corresponding to the chosen , for any value of the dimensionless length-scale parameter .

The transverse dimensionless deflection of the nanocantilever is reported in Figure 3 for and . The vanishing of the length-scale parameter ensures that the GM provides the same solution of the classical (local) model. On the contrary, the FUGM induces a nonlocal behaviour due to the presence of the parameter as shown in Figure 3. As a consequence, the nanocantilever becomes stiffer than the local one for increasing values of the nonlocal parameter .

The maximum dimensionless deflection of the nanocantilever obtained by the FUGM is plotted in Figure 4 in terms of the dimensionless parameter for different values of . The magnitude of ranges in the interval , , and all the plots have a horizontal asymptote for ; that is, .

Denoting by the maximum dimensionless deflection of the nanocantilever obtained by the GM, the comparison between and is reported in Table 1 where the variation between the two maximum dimensionless deflections is evaluated according to formula .

The maximum dimensionless deflection in the classical (local) Euler-Bernoulli model is and is recovered from the FUGM by setting and ; that is, .

Table 1 clearly shows that the maximum dimensionless deflection evaluated by the FUGM quickly decreases for increasing with respect to the corresponding one of the GM. Analogously, for a given , the FUGM becomes stiffer than the corresponding GM for increasing .

#### 6. Conclusions

A fully gradient nanobeam model (FUGM) is presented for the static bending of nanobeams based on a consistent thermodynamic approach. The governing equations and the related high-order boundary conditions are derived by using a variational formulation of the nonlocal model. The proposed theory can be specialized to recover the gradient elasticity beam model (GM).

Analytical solutions of the proposed nonlocal fully gradient nanobeam model in terms of transverse displacements are provided for nanocantilevers.

Present results are compared with the corresponding ones obtained using the gradient elasticity model to establish the validity of the present formulation.

The novelty of the present paper arises from the several points. Firstly, as stated above, a more general nonlocal nanobeam model that contains existing beam theories (classical and gradient elasticity theories) as special cases is proposed. Two nonlocal parameters are introduced in the FUGM in order to consider the nonlocal contributions in the beam axis direction and in the transverse direction.

Secondly, the following implications of technical interest are reported below.(1)The proposed FUGM shows that an increasing of the nonlocal parameter associated with the strain gradient component in the transverse direction provides an increasing of the nanostructural stiffness so that the deflection of a nanobeam reduces with respect to the corresponding one of the local model; see Figures 2, 3, and 4.(2)The existing gradient elasticity model (GM) of Euler-Bernoulli nanobeams overestimates the displacement solution when compared with the new gradient elasticity model (FUGM) as reported in Figures 2-3.(3)The proposed FUGM provides upper and lower bounds for displacements in terms of the parameter . The upper bound is independent of the small-scale parameters while the lower bound depends on the nonlocal parameter .(4)The results given in Table 1 can be useful benchmarks for other researchers to compare their outcomes in the future.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

The present study has been developed within the activities of the research program FARO 2012, Compagnia San Paolo, Polo delle Scienze e delle Tecnologie, University of Naples Federico II.