Research Article  Open Access
A Nonlocal Model for Carbon Nanotubes under Axial Loads
Abstract
Various beam theories are formulated in literature using the nonlocal differential constitutive relation proposed by Eringen. A new variational framework is derived in the present paper by following a consistent thermodynamic approach based on a nonlocal constitutive law of gradienttype. Contrary to the results obtained by Eringen, the new model exhibits the nonlocality effect also for constant axial load distributions. The treatment can be adopted to get new benchmarks for numerical analyses.
1. Introduction
Carbon nanotubes (CNTs) are a topic of major interest both from theoretical and applicative points of view. This subject is widely investigated in literature to describe smallscale effects [1–4], vibration and buckling [5–13], and nonlocal finite element analysis [14–18]. A comprehensive review on applications of nonlocal elastic models for CNTs is reported in [19] and therein references. Buckling of triplewalled CNTs under temperature fields is dealt with in [20]. An alternative methodology is based on an atomisticbased approach [21] which predicts the positions of atoms in terms of interactive forces and boundary conditions. The standard approach to analyze CNTs under axial loads consists in solving an inhomogeneous secondorder ordinary differential equation providing the axial displacement field, see, for example, [22]. The known term of the differential equation is the sum of two contributions. The former describes the local effects linearly depending on the axial load. The latter characterizes the smallscale effects depending linearly on the second derivative along the rod axis of the axial load. This model is thus not able to evaluate smallscale effects due to, constant axial loads per unit length. This approach, commonly adopted in literature, is based on the following nonlocal linearly elastic constitutive law proposed by Eringen [23]: where is a material constant, is the internal length, is the Young modulus, is the normal stress, the apex is second derivative along the rod axis, and is the axial elongation. Indeed, integrating on the rod cross section domain and imposing that the axial force is equal to the resultant of normal stress field we get the differential equation where , with being first derivative along the rod axis of the axial displacement field , where is the rod length and denotes the cross section area. Since the equilibrium prescribes that the first derivative of is opposite to the axial load , we infer the wellknown differential equation (see, e.g., [7]) as follows: Note that the nonlocal contribution vanishes for constant loads . In the present paper, an alternative nonlocal constitutive behavior is adopted to assess smallscale effects in nanotubes also for constant axial loads. The corresponding axial displacement field is shown to be governed by a fourthorder inhomogeneous differential equation. Boundary conditions are naturally inferred by performing a standard localization procedure of a variational problem formulated by making recourse to thermodynamic restrictions see, for example, [24–26], according to the geometric approach illustrated in [27–30]. As an example, the displacement field of nanotubes under constant axial loads per unit length is evaluated in the appendix. Vibration and buckling effects are not the subject of this paper and will be addressed in a forthcoming paper.
2. Nonlocal Variational Formulation
Let be the threedimensional spatial domain of a straight rod subjected to axial loads. An apex stands for th derivative along the rod centroidal axis. Kinematic compatibility between axial elongations and axial displacements is expressed by the differential equation . Denoting by a dot the timerate, the following noteworthy relations hold true The differential equation of equilibrium turns out to be . Boundary equilibrium prescribes that at the end cross sections act axial loads equal to for and to for . Let us now consider a nonlocal constitutive model of gradienttype defined by assigning the following elastic energy functional per unit volume: with being nonlocal parameter. Relation (5) is similar to the elastic energy density proposed in [31] where a homogeneous quadratic functional including also mixed terms is assumed. The elastic energy time rate is hence expressed by the formula where are the static variables conjugating with the kinematic variables and . The static variable is the scalar counterpart of the socalled double stress tensor [31]. By imposing the thermodynamic condition (see, e.g., [32–34]) where is the normal stress, we infer the relation The relevant differential and boundary equations are thus obtained as shown hereafter. Substituting the expression of the rates and in terms of the axial displacement of the cross section at abscissa , we get the formulae with axial force (static equivalence condition on the cross sections) and for . Thermodynamic condition (9) provides the axial contribution
3. Differential and Boundary Equations of Elastic Equilibrium
Resorting to Green’s formula, a standard localization procedure provides the differential and boundary equations corresponding to the variational conditions inferred in Section 2, as follows. A direct computation gives
Substituting into the variational condition (11), a suitable localization provides the relevant differential equation and boundary conditions
These conditions can be conveniently expressed in terms of the axial displacement field as follows. A direct evaluation of the scalar functions for and of their derivatives gives with . Accordingly, the boundary and differential conditions of elastic equilibrium (13) and (14) take the form
4. Example
Let us consider a straight rod subject to a constant axial load as depicted in Figure 1. End cross sections and are assumed to be hinged and simply supported, respectively. As illustrated in Section 3, the computation of the rod axial displacement field involves the following cross section geometric and elastic properties: area , Young modulus , and nonlocal parameter . By setting and , the differential equation of elastic equilibrium is as follows: The general integral takes thus the form (see the appendix) with The evaluation of the constants is carried out by imposing the following boundary conditions (see also Section 3): Resorting to the expressions of the derivatives and for , and having , a direct computation provides the algebraic system A further condition can be obtained by imposing that the scalar field is a particular solution of the differential equation (17), whence it follows that . The remaining constants are given by the formulae having . A plot of the rod axial displacement field for different values of the nonlocal parameter is provided in Figure 2. It is apparent that the rod becomes stiffer if the nonlocal parameter increases. The evaluated axial displacement at the free end of the rod provides the same value independently of the nonlocal parameter. Such a value coincides with the displacement of the point if a local model is considered. Moreover, the limit of the axial displacement field for tending to plus infinity can be evaluated to get the lower bound Hence, large values of the nonlocal parameter provide a displacement field which tends to a linear one, see Figure 2, for . Further, the limit value of the axial displacement for and , obtained by (25), yields nm which coincides with the axial displacement at for any value of the nonlocal parameter , see Figure 3 and Table 1.

The upper bound of the axial displacement is provided by the local solution (i.e., ) The axial displacement evaluated for by (26) yields the value nm which coincides with the axial displacement at for any value of , see Figure 3 and Table 1. For the considered model, the upper and lower bounds of the axial displacement field are given by (25) and (26). The axial displacement at the middle point of the rod and the maximum axial displacement as functions of the nonlocal parameter are depicted in Figure 3. The corresponding numerical values of and are listed in Table 1.
It is worth noting that equilibrium prescribes that axial force must be a linear function, confirmed by the blue diagram in Figure 4 obtained as difference between the local contribution (dashed line) and the nonlocal one (continuous thin line), according to (14)_{1} for any value of .
5. Conclusions
The outcomes of the present paper may be summarized as follows. (i) Linearly elastic carbon nanotubes under axial loads have been investigated by a nonlocal variational approach based on thermodynamic restrictions. The treatment provides an effective tool to evaluate smallscale effects in nanotubes subject also to constant axial loads, a goal not achievable by the Eringen model commonly adopted in literature as motivated in Section 1. (ii)Relevant boundary and differential conditions of elastic equilibrium have been inferred by a standard localization procedure. Such a procedure provides, in a consistent way, the relevant class of boundary conditions for the nonlocal model. (iii)The present approach yields a firm thermodynamic procedure to derive different nonlocal models for CNTs by suitable specializations of the elastic energy. (iv)Exact solutions of carbon nanotubes subject to a constant axial load have been obtained. An advantage of the proposed procedure consists in providing an effective tool to be used as a benchmark for numerical analyses. Finally, a range to which any nonlocal solution must belong is analytically evaluated.
Appendix
The procedure to solve the ordinary differential equation with , being constant coefficients and being a continuous function, is summarized as follows. Let us consider the homogeneous differential equation and the relevant characteristic (algebraic) equation . The roots of the polynomial are with multiplicity 2, with multiplicity 1 and with multiplicity 1. The general integral of (A.2) is thus expressed by the formula with denoting exponential function and for . The general integral of (A.1) is writen therefore as where is a particular solution of (A.1). It is worth noting that, for defined by a polynomial of degree , the solution can be looked for by setting with for .
Acknowledgments
The authors were supported by the “Polo delle Scienze e delle Tecnologie,” University of Naples Federico II, through the research project FARO. Useful hints and precious comments by anonymous reviewers are also gratefully acknowledged.
References
 Q. Wang and K. M. Liew, “Application of nonlocal continuum mechanics to static analysis of micro and nanostructures,” Physics Letters A, vol. 363, no. 3, pp. 236–242, 2007. View at: Publisher Site  Google Scholar
 M. Aydogdu, “A general nonlocal beam theory: its application to nanobeam bending, buckling and vibration,” Physica E, vol. 41, no. 9, pp. 1651–1655, 2009. View at: Publisher Site  Google Scholar
 Ö. Civalek and Ç. Demir, “Bending analysis of microtubules using nonlocal EulerBernoulli beam theory,” Applied Mathematical Modelling, vol. 35, no. 5, pp. 2053–2067, 2011. View at: Publisher Site  Google Scholar
 M. A. De Rosa and C. Franciosi, “A simple approach to detect the nonlocal effects in the static analysis of EulerBernoulli and Timoshenko beams,” Mechanics Research Communications, vol. 48, pp. 66–69, 2013. View at: Publisher Site  Google Scholar
 M. Aydogdu, “Axial vibration analysis of nanorods (carbon nanotubes) embedded in an elastic medium using nonlocal elasticity,” Mechanics Research Communications, vol. 41, pp. 34–40, 2012. View at: Publisher Site  Google Scholar
 M. A. KazemiLari, S. A. Fazelzadeh, and E. Ghavanloo, “Nonconservative instability of cantilever carbon nanotubes resting on viscoelastic foundation,” Physica E, vol. 44, no. 78, pp. 1623–1630, 2012. View at: Publisher Site  Google Scholar
 H.T. Thai and T. P. Vo, “A nonlocal sinusoidal shear deformation beam theory with application to bending, buckling, and vibration of nanobeams,” International Journal of Engineering Science, vol. 54, pp. 58–66, 2012. View at: Publisher Site  Google Scholar
 M. A. Eltaher, S. A. Emam, and F. F. Mahmoud, “Static and stability analysis of nonlocal functionally graded nanobeams,” Composite Structures, vol. 96, pp. 82–88, 2013. View at: Publisher Site  Google Scholar
 S. A. Emam, “A general nonlocal nonlinear model for buckling of nanobeams,” Applied Mathematical Modelling, vol. 37, no. 1011, pp. 6929–6939, 2013. View at: Publisher Site  Google Scholar
 B. Fang, Y.X. Zhen, C.P. Zhang, and Y. Tang, “Nonlinear vibration analysis of doublewalled carbon nanotubes based on nonlocal elasticity theory,” Applied Mathematical Modelling, vol. 37, no. 3, pp. 1096–1107, 2013. View at: Publisher Site  Google Scholar
 S. A. M. Ghannadpour, B. Mohammadi, and J. Fazilati, “Bending, buckling and vibration problems of nonlocal Euler beams using Ritz method,” Composite Structures, vol. 96, pp. 584–589, 2013. View at: Publisher Site  Google Scholar
 M. Şimşek and H. H. Yurtcu, “Analytical solutions for bending and buckling of functionally graded nanobeams based on the nonlocal Timoshenko beam theory,” Composite Structures, vol. 97, pp. 378–386, 2013. View at: Publisher Site  Google Scholar
 B. L. Wang and K. F. Wang, “Vibration analysis of embedded nanotubes using nonlocal continuum theory,” Composites B, vol. 47, pp. 96–101, 2013. View at: Publisher Site  Google Scholar
 J. K. Phadikar and S. C. Pradhan, “Variational formulation and finite element analysis for nonlocal elastic nanobeams and nanoplates,” Computational Materials Science, vol. 49, no. 3, pp. 492–499, 2010. View at: Publisher Site  Google Scholar
 C. M. C. Roque, A. J. M. Ferreira, and J. N. Reddy, “Analysis of Timoshenko nanobeams with a nonlocal formulation and meshless method,” International Journal of Engineering Science, vol. 49, no. 9, pp. 976–984, 2011. View at: Publisher Site  Google Scholar
 F. F. Mahmoud, M. A. Eltaher, A. E. Alshorbagy, and E. I. Meletis, “Static analysis of nanobeams including surface effects by nonlocal finite elements,” Journal of Mechanical Science and Technology, vol. 26, no. 11, pp. 3555–3563, 2012. View at: Publisher Site  Google Scholar
 S. C. Pradhan, “Nonlocal finite element analysis and small scale effects of CNTs with Timoshenko beam theory,” Finite Elements in Analysis and Design, vol. 50, pp. 8–20, 2012. View at: Publisher Site  Google Scholar
 M. A. Eltaher, A. E. Alshorbagy, and F. F. Mahmoud, “Vibration analysis of EulerBernoulli nanobeams by using finite element method,” Applied Mathematical Modelling, vol. 37, no. 7, pp. 4787–4797, 2013. View at: Publisher Site  Google Scholar
 B. Arash and Q. Wang, “A review on the application of nonlocal elastic models in modeling of carbon nanotubes and graphenes,” Computational Materials Science, vol. 51, no. 1, pp. 303–313, 2012. View at: Publisher Site  Google Scholar
 Y. Yan, W. Q. Wang, and L. X. Zhang, “Nonlocal effect on axially compressed buckling of triplewalled carbon nanotubes under temperature field,” Applied Mathematical Modelling, vol. 34, no. 11, pp. 3422–3429, 2010. View at: Publisher Site  Google Scholar
 R. Rafiee and R. M. Moghadam, “On the modeling of carbon nanotubes: a critical review,” Composites B, vol. 56, pp. 435–4490, 2014. View at: Publisher Site  Google Scholar
 J. N. Reddy, “Nonlocal theories for bending, buckling and vibration of beams,” International Journal of Engineering Science, vol. 45, no. 28, pp. 288–307, 2007. View at: Publisher Site  Google Scholar
 A. C. Eringen, “On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves,” Journal of Applied Physics, vol. 54, no. 9, pp. 4703–4710, 1983. View at: Publisher Site  Google Scholar
 F. Marotti de Sciarra, “Variational formulations, convergence and stability properties in nonlocal elastoplasticity,” International Journal of Solids and Structures, vol. 45, no. 78, pp. 2322–2354, 2008. View at: Publisher Site  Google Scholar
 G. Romano, M. Diaco, and R. Barretta, “Variational formulation of the first principle of continuum thermodynamics,” Continuum Mechanics and Thermodynamics, vol. 22, no. 3, pp. 177–187, 2010. View at: Publisher Site  Google Scholar
 F. Marotti De Sciarra, “Hardening plasticity with nonlocal strain damage,” International Journal of Plasticity, vol. 34, pp. 114–138, 2012. View at: Publisher Site  Google Scholar
 G. Romano and R. Barretta, “Covariant hypoelasticity,” European Journal of Mechanics A, vol. 30, no. 6, pp. 1012–1023, 2011. View at: Publisher Site  Google Scholar
 G. Romano and R. Barretta, “On Euler's stretching formula in continuum mechanics,” Acta Mechanica, vol. 224, no. 1, pp. 211–230, 2013. View at: Publisher Site  Google Scholar
 G. Romano and R. Barretta, “Geometric constitutive theory and frame invariance,” International Journal of NonLinear Mechanics, vol. 51, pp. 75–86, 2013. View at: Publisher Site  Google Scholar
 G. Romano, R. Barretta, and M. Diaco, “Geometric continuum mechanics,” Meccanica, 2013. View at: Publisher Site  Google Scholar
 R. D. Mindlin, “Microstructure in linear elasticity,” Archive for Rational Mechanics and Analysis, vol. 16, no. 1, pp. 51–78, 1964. View at: Publisher Site  Google Scholar
 F. Marotti de Sciarra, “Novel variational formulations for nonlocal plasticity,” International Journal of Plasticity, vol. 25, no. 2, pp. 302–331, 2009. View at: Publisher Site  Google Scholar
 F. Marotti de Sciarra, “On nonlocal and nonhomogeneous elastic continua,” International Journal of Solids and Structures, vol. 46, no. 34, pp. 651–676, 2009. View at: Publisher Site  Google Scholar
 F. Marotti de Sciarra, “A nonlocal model with strainbased damage,” International Journal of Solids and Structures, vol. 46, no. 2223, pp. 4107–4122, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Raffaele Barretta and Francesco Marotti de Sciarra. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.