Modified Nonlocal Strain Gradient Elasticity for NanoRods and Application to Carbon Nanotubes

Nowadays, the modified nonlocal strain gradient theory provides a mathematically well-posed and technically reliable methodology to assess scale effects in inflected nano-structures. Such an approach is extended in this paper to investigate the extensional behavior of nano-rods. The considered integral elasticity model, involving axial force and strain fields, is conveniently shown to be equivalent to a nonlocal differential problem equipped with constitutive boundary conditions. Unlike treatments in the literature, no higher-order boundary conditions are required to close the nonlocal problem. Closed-form solutions of elastic nano-rods under selected loadings and kinematic boundary conditions are provided. As an innovative implication, Young’s moduli of Single-Walled Carbon Nanotubes (SWCNT) weare assessed and compared with predictions of Molecular Dynamics (MD). New benchmarks for numerical analyses were also detected.


Introduction
Components of Nano-Electro-Mechanical-Systems (NEMS), such as sensors and actuators, are usually modeled as nano-rods and nano-beams.It is well known that methods of local continuum mechanics cannot be adopted for such elements.Molecular Dynamics (MD) simulations are time-consuming and micro/nano-scaled experiments are usually difficult to implement.Hence, nonlocal continuum models have been developed for predicting the size-dependent mechanical behavior of nano-structures.
In the framework of nonlocal elasticity, Eringen's strain-driven differential model [1] has been widely adopted in the literature (see, e.g., [2][3][4][5][6]).It is worth noting that recent papers on the strain-driven nonlocal model [7][8][9][10] prove that, if the bending field is expressed as convolution of elastic curvature with an averaging kernel assuming an exponential expression, a solution of this problem exists only if the bending field satisfies constitutive boundary conditions.Accordingly, it is shown in [11] that the nonlocal elastostatic problem is ill-posed in all cases of applicative interest, as acknowledged in the literature (see, e.g., [12][13][14][15][16][17][18]).A modified version of the Eringen integral model is proposed in [19] and has been recently applied to inflected nanobeams in [20].
The gradient elasticity theory [21] assumes that a material at the nano-scale is modeled via gradient terms.Many works investigate the small-scale effects on the static and dynamic behaviors of rods, beams and plates and the effect of stiffness enhancement has been often reported in these strain gradient models (see, e.g., [22][23][24][25][26]).
Recently, the Eringen's integral law [1] has been combined with the strain gradient elasticity in [27] to formulate a higher-order nonlocal theory, thus collecting nonlocal theory and strain gradient theory into a single model.
The procedure consists in considering the integral nonlocal gradient method for structural problems defined on bounded domains equivalent to a differential law of higher-order than the one of the classical local problem.Therefore, additional non-classical suitable boundary conditions must be added to solve the nonlocal strain gradient elastostatic problem.
To solve the problem, different choices have been followed in the literature.Two usual choices consist in imposing higher-order boundary conditions pertaining to the strain gradient theory, of kinematic type [38] or static type [39].It is worth noting that the structural behavior is greatly influenced by such choices.
The nonlocal strain gradient theory with higher-order boundary conditions has been recently adopted in [40] to study nano-rods in tension.The closed-form solutions for predicting the axial displacement and the variation of the Young's modulus have been derived for four different nano-rods differing by the choice of the higher-order boundary conditions.In particular, for each nano-rod, the nonlocal parameters have been set to match the variation of the Young's modulus obtained by the MD simulations.
The choice of the higher-order boundary conditions is disputed in the literature and is considered an open question (see, e.g., [41]).In the context of nano-beams subjected to flexure, a recent contribution [42] provides a definite solution to this issue.In fact, the non-classical boundary conditions to be imposed to solve the elastostatic problem of nonlocal strain gradient inflected nano-beams are given by constitutive boundary conditions (CBC) that naturally follow from the nonlocal strain gradient integral model.The consistent nonlocal strain gradient strategy has been successfully applied to free vibrations of nano-beams in [43].
In the present paper, the structural behavior of nano-rods in tension is formulated in the framework of the modified nonlocal strain gradient (NSG) theory.The expressions of the CBC for nano-rods are explicitly provided and it is shown that no unmotivated higher-order boundary conditions have to be prescribed to solve the nonlocal structural problem.
In addition, the variation of the Young's modulus provided by MD simulations is recovered based on the NSG model for nano-rods developed in the present paper.In particular, carbon nanotubes are effectively described by NSG nano-rods with the usual boundary conditions, that is clamped at the one end and with a tensile force at the other end.As illustrated in Appendix A, the same result is obtained if a doubly-clamped nano-rod with an imposed axial displacement at one end is considered.
Finally, numerical analyses are presented as benchmark examples for applications and experimental tests on nonlocal nano-rods.

Modified Nonlocal Strain Gradient Law for Rods
Let us consider a functionally graded (FG) straight nano-rod of length L, the x-coordinate is taken along the length of the nano-rod with the y-coordinate along the thickness and the z-coordinate along the width of the nano-rod.The local Young's modulus E of the FG nano-rod continuously changes in the thickness direction y, so that the Young's elastic modulus at the point y is E(y) and the elastic area is A E = Ω E (y) dA, being Ω the nano-rod cross-section.
In the modified nonlocal strain gradient (NSG) model for FG nano-rods, the axial force N is defined in terms of elastic axial strain ε el and of its derivative (1) The smoothing kernels φ 0 and φ 1 depend on two non-dimensional nonlocal parameters λ 0 > 0 and λ 1 > 0. The scale parameter l > 0, characteristic of the strain gradient elasticity [21], was subsequently introduced in [27] to make dimensionally homogeneous the convolutions in Equation (1) and to describe the importance of higher-order strain gradient fields.
Introducing the following fields the modified nonlocal strain gradient elastic law (Equation ( 1)) can be rewritten as As proven in the next proposition, the modified nonlocal strain gradient integral relation (Equation ( 4)) for FG nano-rods is equivalent to a suitable differential law with constitutive boundary conditions.Proposition 1 (Constitutive equivalence for FG nano-rods).The modified nonlocal strain gradient constitutive law (Equation (4)) equipped with the bi-exponential kernel (Equation (2)) with x ∈ [0, L], is equivalent to the differential relation subject to the following two constitutive boundary conditions (CBC) Proof.Since the bi-exponential averaging function is given by a direct evaluation provides the first derivative of the convolutions (Equation (3)) The second derivative of the convolutions (Equation ( 3)) follows from Equation (9) using Equation (3) to get Subtracting the third derivative of Equation ( 9) 2 from Equation (9) 1 , it turns out to be so that, recalling Equation ( 5) and rearranging the terms, we have and Equation ( 6) is recovered.
The CBC in Equation ( 7) can be recovered as follows.
Using Equation (9) 2 , Equation ( 5) can be rewritten in the form and using Equations (9) 1 and (7) 2 , the first derivative of Equation (5) becomes The CBC (Equation ( 7)) of modified nonlocal strain gradient nano-rods follows by evaluating Equations ( 13) and ( 14) at nano-rod boundary points x = 0 and x = L.In fact, we have at x = 0 so that Equation (15) provides the relations and the CBC in Equation ( 7) 1 is recovered.Analogously, setting x = L in Equation ( 14), we get so that Equation ( 17) provides the relations and the CBC in Equation ( 7) 2 is recovered.Conversely, sufficient condition can be inferred from the uniqueness of the solution of Equation ( 6) consequent to the fact that the associated homogeneous equations admit only the trivial solution under the homogeneous boundary conditions for (19) 1 and for (19) 2 respectively.

Elastic Equilibrium Problem
Let us consider a FG nano-rod subject to a distributed axial load q(x) per unit length in the interval [0, L] and to concentrated axial forces F at the end cross-sections x = 0 and x = L.
Differential condition of equilibrium can be written as with the boundary conditions N (x, L c , l) = ∓F at x = 0 and x = L.
The axial displacement at the abscissa x along the nano-rod axis is denoted by u(x) and the kinematically compatible axial strain has the form In the sequel, elastic ε el and kinematically compatible ε strains are assumed to be coincident.Exact solutions according to the proposed modified nonlocal strain gradient (NSG) model for FG nano-rods can be performed by the following steps.

•
Step 1: Solve the equilibrium Equation ( 22) to get the expression of the axial force • Step 2: Solve the second-order differential Equation (6) in the form obtaining the expression of the elastic axial strain ε el of the nano-rod in terms of three integration constants (a 1 , a 2 , and a 3 ) to be determined.

•
Step 3: Solve the first-order differential in Equation ( 23) in terms of the axial displacement u of the nano-rod to get the expression of u in terms of four integration constants (a 1 , a 2 , a 3 , and a 4 ) to be determined.

•
Step 4: Determine the four integration constants (a 1 , a 2 , a 3 , and a 4 ) by imposing the two CBC given by Equation ( 7) in terms of the axial displacement u and the two classical boundary conditions at the nano-rod end points x = 0 and x = L by specifying u or N.
It is worth noting that, in statically determinate rods, the axial force N can be obtained by Equation (24) by imposing the classical static boundary conditions.

Closed-Form Solutions for FG Nano-Rods
Closed-form elastic solutions for FG nano-rods with a clamped end at x = 0 and a free end at x = L and with both clamped ends are presented hereafter.The applied loads are a uniform load p, a concentrated force F at x = L (or an imposed axial displacement δ at x = L for doubly-clamped rods).Kinematic and static boundary conditions are enforced to the FG nano-rod ends according to classical rod theory and, in addition, the constitutive boundary conditions Equation ( 26) are imposed according to the proposed NSG model.Hence, the axial displacement u can be recovered following Steps 1-4 in Section 3.
In the sequel, the abbreviations CF and CC stand for clamped-free andclamped-clamped, respectively.Moreover, let us assume that the elastic area A E is constant along the nano-rod axis x.
For completeness sake, the FG nano-rod constraints, the considered applied load and the related boundary conditions are reported in Table 1.
Table 1.Fundamental schemes with corresponding boundary conditions.

FG Nano-Rod Constraints-Applied Load Boundary Conditions
Classical Constitutive To provide a non-dimensional analysis of FG nano-rods, the following non-dimensional variable ξ and the non-dimensional characteristic parameters λ and µ are adopted in the examples The non-dimensional axial displacement ū depends on the kind of load applied to the nano-rod according to Table 2.
Table 2. Non-dimensional axial displacement ū in terms of the applied load.

Non-Dimensional Axial Displacement
Applied Load Case I: CF FG Nano-Rod with a Concentrated Load at the Free End Let us consider a FG nano-rod of length L with a clamped end at x = 0 and a free end at x = L subject to a concentrated load F at the free end.
Following the steps in Section 3, the axial force N can be evaluated by means of the equilibrium equation so that Equation (24), with the boundary condition N(L) = F , yields N(x) = F .For simplicity, we drop the dependence on the nonlocal characteristic parameters (L c , l).Then, the analytical solution of the NSG model for the FG nano-rod is obtained from the following nonlocal differential equation under the CBC (Equation ( 26)) and the kinematic boundary condition in Equation ( 27) Hence, the solution of the differential equation (Equation ( 29)) with the boundary conditions in Equations ( 30) and (31) provides the axial displacement where u e is the rod axial displacement of the local model The maximum displacement takes place at x = L and can be obtained from Equation ( 32) The classical (local) displacement u e of the FG nano-rod is provided by L c → 0.Moreover, the limit nano-rod displacement u ∞ for l → +∞ is given by It is important to note that the NSG model for FG nano-rods exerts a softening effect, with respect to the local behavior, in terms of the nonlocal parameter L c .It is of interest that the maximum displacement u(L) of the NSG rod model does not depend on the gradient parameter l and the maximum axial displacement tends to the one of the classical (local) rod if L c → 0.
The effects of the non-dimensional characteristic parameters λ and µ on the elastic response of nano-rods are examined in Figures 1 and 2.
Figure 1a,b show the non-dimensional axial displacement ū in terms of the gradient non-dimensional parameter µ for λ = 0.4 and λ = 0.8, respectively.The local non-dimensional axial displacement ūe is recovered by λ → 0 + and is reported with the dot dashed line.The limit non-dimensional axial displacement ū∞ follows from Equation (35) and is reported with the dotted line in terms of the non-dimensional nonlocal parameter λ.
The comparison between nonlocal FG nano-rods and classical (local) FG rods in Figure 1 highlights the increment of the axial displacement ū due to the behavior of the nonlocal model.The parameter λ has the effect of increasing the axial displacement, i.e. a larger λ involves greater axial displacements ū for a given value of the non-dimensional gradient parameter µ.
Figure 2a shows the non-dimensional axial displacement ū of the FG nano-rod in terms of the nonlocal non-dimensional parameter λ for µ = 0.15.The local non-dimensional axial displacement ūe is recovered by λ → 0 + and is reported with the dot dashed line.The limit nano-rod non-dimensional axial displacement ū∞ in terms of λ is plotted with the dotted line.
Figure 2b shows the non-dimensional axial displacement ū of the FG nano-rod in terms of the nonlocal non-dimensional parameter λ for µ = 0.15 (thick lines) and µ = 0.30 (dotted lines).The non-dimensional maximum axial displacement ū(1) increases for increasing values of the nonlocal parameter λ and is independent of the non-dimensional gradient parameter µ.The limit nano-rod non-dimensional axial displacement ū∞ is plotted with the black dotted line.The 3D plot of the non-dimensional maximum axial displacement ū(1) for the proposed NSG method versus the non-dimensional characteristic parameters λ and µ is reported in Figure 3.The horizontal plane is the non-dimensional local maximum axial displacement ūe (1) = 1.
The innovative nonlocal model exhibits a hardening behavior in terms of the non-dimensional characteristic parameter λ and the maximum axial displacement does not depend on the non-dimensional gradient parameter µ.

Reduced Young's Modulus
In applications, the evaluation of the Young's modulus of micro-and FG nano-rods is of great interest.Hence, we evaluate the reduced rigidity K r from Equation (34) to get Note that the classical rigidity K r = A E can be recovered if L c → 0.
Using the reduced rigidity in Equation ( 36), the axial displacement of the modified nonlocal strain gradient model can be calculated by adopting the classical (local) analysis for rods.Let us consider a constant Young's elastic modulus E(y) = E so that the elastic area is A E = Ω E (y) dA = EA.Hence, the reduced Young's modulus E r can be defined from Equation ( 36) as depending on the nonlocal parameter L c .The upper bound of the nano-rod reduced Young's modulus E r for L → ∞ is provided by the local Young's modulus E. Analogously, for L c → 0, we recover the local Young's modulus E and for L c → +∞, the nano-rod reduced Young's modulus tends to vanish.
To make a comparison, we considered the data presented in [40] for a SWCNT of armchair (10,10).The diameter d of the SWCNT (n, m) can be calculated by where the carbon-carbon bond length is a = 0.142 nm.The effective thickness of the considered SWCNT is t = 0.34 nm and the classical (local) Young's modulus assumed in [40,44] is E = 909.5 GPa.The Young's modulus predicted by the MD simulations is reported in [40,44].Note that the strain gradient elastic theory with high-order boundary conditions was adopted in [40] and the necessity to impose higher-order (non-classical) boundary conditions has the effect that the classical boundary conditions of the classical theory (such as free and clamped boundary conditions) may no longer be meaningful for the modified nonlocal strain gradient rod.As a result, following [28,35], one needs to take into account further boundary conditions involving higher-order stress and strain distributions.Accordingly, four nonlocal nano-rods are considered in [40], depending on the considered higher-order boundary conditions, and for each of them the small-scale parameters are set to match the results of the MD simulations.
On the contrary, the proposed NSG model has no higher-order boundary condition to add so that the classical definitions of external constraints must not be modified.Hence, a unique model of rod can be considered and a unique value of the nonlocal parameter has to be set to match the results of the MD simulations.
In Figure 4 the results provided by Equation ( 37) are plotted together with the MD data versus the SWCNT length.Upper and lower bounds of the nonlocal parameter L c = 0.04272 nm and L c = 0.06942 nm are reported in Figure 4a to include the values provided by the MD simulations.A good agreement between the Young's modulus obtained by the NSG model and the MD results could be obtained by setting the nonlocal parameter L c = 0.0534 nm, as shown in Figure 4b.
The NSG model provides values of the reduced Young's modulus E r tending to the classical (local) value E for increasing values of the SWCNT length L. On the contrary, the values of the MD simulations appear to be constant for values of the SWCNT's length greater than 27 nm.
The small-scale effect on displacements can be clearly observed in Figure 5 where the non-dimensional maximum displacement u(L)/L, pertaining to the NSG model, is plotted versus the SWCNT's length L for an applied force F = 1 nN and the nonlocal parameter L c ranging in the set {0.04272, 0.0534, 0.06942} nm.As can be seen, the small-scale effect on the displacement can be observed when the length of SWCNT is small and the small-scale effect increases for increasing values of the nonlocal parameter L c .If L → +∞, the non-dimensional maximum displacement u(L)/L tends to the corresponding local one 1/A E = 7.59117 × 10 −4 .(10,10) having diameter d = 1.356 nm, effective thickness t = 0.34 nm and classical (local) Young's modulus E = 909.5 GPa simulated by a CF nano-rod with a concentrated load F = 1 nN at the free end.Plot of the non-dimensional maximum displacement u(L)/L of the NSG model vs. the SWCNT's length L for the nonlocal parameter L c ranging in the set {0.04272, 0.0534, 0.06942} nm.The non-dimensional local maximum displacement (black dotted line) is 1/A E = 7.59117 × 10 −4 .
In [45], a continuum mechanics model has been proposed to predict the effective wall thickness of a SWCNT and to calculate its Young's modulus.The deformation of a central long SWCNT in a bundle of SWCNTs, subjected to an external pressure, has been considered in plane-strain and has been modeled as a thin ring with a mean radius R , thickness t in the radial direction and a unit width in the axial direction.Hence, it has been obtained that the radius of the nanotube is R = 0.7066 nm and the predicted thickness is t = 0.0617 nm so that the related Young's modulus is E = 4880 GPa.Considering the nonlocal parameter (dependent on the longitudinal atom spacing in armchair CNTs [46,47]) L c = 0.0534 nm as previously calibrated by means of the MD data, the variation law of Young's modulus obtained by the proposed NSG model for the SWCNT investigated in [45] is reported in Figure 6, with radius R = 0.7066 nm, thickness t = 0.0617 nm and Young's modulus E = 4880 GPa.The NSG model provides values of the reduced Young's modulus E r tending to the value E = 4880 GPa for increasing values of the SWCNT length L. It is apparent that, for a nanotube with L = 4 nm, the reduced Young's modulus is E r = 4753 GPa, achieving 97.4% of the value of E .Hence, we can conclude that, for the SWCNT considered in [45], the variation of the Young's modulus E r is really small if the length of the SWCNT is greater that 4 nm.The small-scale effect on displacements is observed in Figure 7, in which the non-dimensional maximum displacement u(L)/L of the NSG model is plotted versus the SWCNT's length L for an applied force F = 1 nN and the nonlocal parameter L c ranging in the set (0.04272, 0.0534, 0.06942) nm.The small-scale effect on the displacement is apparent and the small-scale effect increases for increasing values of the nonlocal parameter L c .The non-dimensional maximum displacement u(L)/L tends to the corresponding local displacement 1/AE = 7.48069 × 10 −4 if L tends to +∞ .

Case II: CF FG Nano-Rod Subject to a Uniformly Distributed Axial Load
Let us now consider a FG nano-rod with a clamped end at x = 0 and a free end at x = L subject to a uniformly distributed axial load q(x) = p.
Following the steps reported in Section 3, the axial force N can be evaluated by means of the equilibrium equation so that Equation (24), with the boundary condition N(L) = 0, yields N(x) = p(L − x).Hence, the closed-form solution of the NSG model for the considered FG nano-rod is obtained by the following nonlocal differential equation under the CBC in Equation ( 26) and the classical boundary condition in Equation ( 27) Hence, the axial displacement is where u e is the axial displacement of the local model The maximum axial displacement of the CF FG nano-rod is attained at the free end x = L and is given by It is apparent from Equation (42) that the axial displacement u depends on the nonlocal parameters l and L c .On the contrary, the maximum axial displacement u(L) (see Equation (44)) is independent of the nonlocal gradient parameter l.
The limit maximum axial FG nano-rod displacement for the nonlocal parameter L c → 0 is given by the classical (local) displacement u e (L) = pL 2 2A E of the FG nano-rod.The non-dimensional axial displacement ū of the FG nano-rod versus the FG nano-rod non-dimensional length ξ is reported in Figure 8 in terms of the nonlocal non-dimensional parameter λ for µ = 0.15 (thick lines) and µ = 0.30 (dotted lines).The non-dimensional maximum axial displacement ū(1) increases for increasing values of the nonlocal parameter λ and is independent of the non-dimensional gradient parameter µ.The limit FG nano-rod non-dimensional axial displacement for the non-dimensional nonlocal parameter λ tending to vanish is given by ūlim and is reported with magenta thick line for µ = 0.15 and dotted line for µ = 0.30.Note that ūlim (1) = ūe (1) = 0.5.
The limit non-dimensional axial displacement ū∞ for µ → 0 is obtained from Equation (42) and is given by ū∞ The corresponding plot is reported with the black dotted line in Figure 8 for the considered values of the non-dimensional nonlocal parameter λ.Let us consider a fully clamped FG nano-rod subject to a uniformly distributed axial load q(x) = p.Following the steps reported in Section 3, the axial force is N (x) = −px + a 1 so that the analytical solution of the FG nano-rod is provided as under the CBC in Equation ( 26) and the classical boundary condition in Equation ( 27) Hence, the axial displacement is where u e is the axial displacement of the local model Note that we have a 1 = pL/2 so that the axial force of the NSG method is N(x) = p(L − 2x)/2 and coincides to the local one.
The non-dimensional axial displacement ū of the FG nano-rod versus the FG nano-rod non-dimensional length ξ is reported in Figure 9 in terms of the nonlocal non-dimensional parameter λ for µ = 0.15 (thick lines) and µ = 0.30 (dotted lines).The non-dimensional midspan axial displacement ū(1/2) increases for increasing values of the nonlocal parameter λ for a given µ.The limit FG nano-rod non-dimensional axial displacement for the non-dimensional nonlocal parameter λ tending to vanish follows from Equation ( 50) and is reported with magenta thick line for µ = 0.15 and dotted line for µ = 0.30.The limit non-dimensional axial displacement ū∞ is obtained from Equation (50) for µ → ∞ and is the vanishing one, i.e., ū∞ (ξ) = 0.The corresponding plot is reported with the black dotted line in Figure 9.
The 3D plot of the non-dimensional maximum axial displacement ū(1/2) for the proposed NSG method versus the non-dimensional characteristic parameters λ and µ is reported in Figure 10.It is apparent that the NSG method stiffness or soften the nano-rod depending on the values of the non-dimensional nonlocal and gradient parameters (λ, µ).
The cuts of the 3D plot for given values of the non-dimensional parameters λ and µ are provided in Figure 11.In particular, the plots of the non-dimensional maximum axial displacement ū(1/2) for the NSG method versus the non-dimensional characteristic parameter λ for µ = 0.2 and µ = 0.6 are provided in Figure 11a.It is apparent that for µ = 0.2 and λ < 0.070156 or µ = 0.6 and λ < 0.4 the NSG model is stiffer than the local model and for µ = 0.2 and λ > 0.070156 or µ = 0.6 and λ > 0.4 the NSG model softens the nano-rod.The limit values of the non-dimensional maximum axial displacement ū(1/2) for λ → 0 + are given by 0.0915228 for µ = 0.2 and by 0.0281989 for µ = 0.6.The plots of the non-dimensional maximum axial displacement ū(1/2) for the NSG method versus the non-dimensional characteristic parameter µ for λ = 0.2 and λ = 0.6 are provided in Figure 11b.It is immediate to note that for λ = 0.2 and µ < 0.374166 or λ = 0.6 and µ < 0.812404 the NSG model is stiffer than the local model and for λ = 0.2 and µ > 0.374166 or λ = 0.6 and µ > 0.812404 the NSG model softens the nano-rod.The limit values of the non-dimensional maximum axial ū(1/2) for µ → 0 + are given by 0.265 for λ = 0.2 and by 0.785 for λ = 0.6.for the NSG method vs. the non-dimensional characteristic parameters µ for λ = 0.2 and λ = 0.6.Finally, a standard numerical analysis of Equation (50) shows that the non-dimensional maximum axial displacement ū(1/2) for the proposed NSG method is attained at the midspan and is reported in Table 3 in terms of the non-dimensional nonlocal parameter λ and gradient parameter µ.It worth noting that the maximum displacement is attained at the midspan ξ = 1/2 of the FG nano-rod independent of the values of λ and µ.The non-dimensional maximum axial displacement of the classical (local) model is ūe (1/2) = 1/8 = 0.125.Accordingly, the non-dimensional maximum axial displacements ū(1/2) in terms of the pairs (λ, µ), which are less than the non-dimensional maximum classical axial displacements ūe (1/2) = 0.125 of the FG nano-rod, are reported in italic in Table 3.
Hence, the italic values of the non-dimensional maximum axial displacement ū(1/2) of the NSG method are smaller than the one of the classical (local) model, thus Table 3 allows one to identify the corresponding pairs (λ, µ) having the effect of stiffen or soften the FG nano-rod with respect the classical (local) behavior.

Conclusions
FG elastic nano-rods under tension have been investigated by the modified nonlocal strain gradient (NSG) theory.The new formulation contains a nonlocal parameter and a material length scale parameter to incorporate the scaling effects of nonlocal stress and microstructure-dependent strain gradient.In comparison to other strain-driven methodologies, the new proposal has been shown to be well-posed and does not require higher-order boundary conditions.In fact, in addition to the classical static and kinematic boundary conditions, closure of the NSG model has to be carried out by prescribing suitable constitutive boundary conditions.Closed-form nonlocal solutions of FG clamped-free and clamped-clamped nano-rods have been provided, exhibiting stiffening or softening effects depending on the values of nonlocal and gradient parameters.Single-Walled Carbon Nanotubes (SWCNT) of armchair (10,10) were modeled as NSG nano-rods, showing that the new approach could capture the small-scale behavior of Young's modulus as predicted by the MD simulations.The nonlocal parameter was thus tuned to characterize Young's modulus vs. the SWCNT length.

Figure 1 .
Figure1.CF nano-rod with a concentrated load at the free end.Plot of the non-dimensional axial displacement ū, limit non-dimensional axial displacement ū∞ (black dotted line) and local non-dimensional axial displacement ūe (black dot dashed line) vs. the non-dimensional nano-rod axis for the gradient non-dimensional parameter µ in the set {0.10, 0.15, 0.20, 0.25} and: (a) nonlocal non-dimensional parameter λ = 0.4; and (b) nonlocal non-dimensional parameter λ = 0.8.

Figure 2 .
Figure 2. CF nano-rod with a concentrated load at the free end.Plot of the non-dimensional axial displacement ū, limit non-dimensional axial displacement ū∞ (black dotted line) and local non-dimensional axial displacement ūe (black dot dashed line) in terms of the nonlocal non-dimensional parameter λ in the set {0.2, 0.4, 0.8, 1.0} and: (a) gradient non-dimensional parameter µ = 0.15; and (b) gradient non-dimensional parameter µ in the set {0.15, 0.30}.The nonlocal model coincides with the classical (local) model of rods for non-dimensional characteristic parameters λ tending to vanishing, i.e., and λ → 0 + .

Figure 3 .
Figure 3. CF nano-rod with a concentrated load at the free end.3D plot of the non-dimensional maximum axial displacement ū(1) vs. the non-dimensional characteristic parameters λ and µ.The non-dimensional local maximum axial displacement ūe (1) = 1 is the horizontal plane.

Figure 4 .Figure 5 .
Figure 4. SWCNT of armchair (10, 10) having diameter d = 1.356 nm, effective thickness t = 0.34 nm and classical (local) Young's modulus E = 909.5 GPa simulated by a CF nano-rod with a concentrated load at the free end.(a) Plot the variation of the Young's modulus obtained by the NSG method together with the MD data versus the SWCNT length.The upper bound is given for the nonlocal parameter L c = 0.04272 nm and the lower bound is given for the nonlocal parameter L c = 0.06942 nm.(b) Variation of the Young's modulus obtained by the NSG model and the MD data for the nonlocal parameter L c = 0.0534 nm.

Figure 6 .
Figure6.Plot the variation of the Young's modulus obtained by the NSG method versus the SWCNT length for the nonlocal parameter L c = 0.0534 nm simulated by a CF nano-rod with a concentrated load at the free end.The data of the SWCNT are: diameter d = 1.4132 nm, effective thickness t = 0.0617 nm and Young's modulus E = 4880 GPa.

Figure 7 .
Figure 7. Plot of the non-dimensional maximum displacement u(L)/L of the NSG model vs. the SWCNT's length L for the nonlocal parameter L c ranging in the set (0.04272, 0.0534, 0.06942) nm.The non-dimensional local maximum displacement (black dotted line) is 1/AE = 7.48069 × 10 −4 .The data of the SWCNT are: diameter d = 1.4132 nm, effective thickness t = 0.0617 nm and Young's modulus E = 4880 GPa.

Figure 8 .
Figure8.CF nano-rod with a uniformly distributed axial load.Plot of the non-dimensional axial displacement ū, limit non-dimensional axial displacement ū∞ (black dotted line) and local non-dimensional axial displacement ūe (black dot dashed line) vs. the nano-rod non-dimensional length ξ in terms of the nonlocal non-dimensional parameter λ in the set {0.2, 0.4, 0.8, 1.0} and gradient non-dimensional parameter µ in the set {0.15, 0.30}.

4. 4 .
Case III: CC FG Nano-Rod Subject to a Uniformly Distributed Axial Load

Figure 9 .
Figure9.CC nano-rod subject to a uniformly distributed axial load.Plot of the non-dimensional axial displacement ū, limit non-dimensional axial displacement ū∞ (black dotted line) and local non-dimensional axial displacement ūe (black dot dashed line) vs. the nano-rod non-dimensional length ξ in terms of the nonlocal non-dimensional parameter λ in the set {0.2, 0.4, 0.8, 1.0} and gradient non-dimensional parameter µ in the set {0.15, 0.30}.

Table 3 .
Non-dimensional maximum axial displacement ū(1/2) for the proposed NSG method in terms of the non-dimensional parameters λ and µ.