Research Article | Open Access
A Nonlocal Model for Carbon Nanotubes under Axial Loads
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 gradient-type. 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.
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 small-scale 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  and therein references. Buckling of triple-walled CNTs under temperature fields is dealt with in . An alternative methodology is based on an atomistic-based approach  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 second-order ordinary differential equation providing the axial displacement field, see, for example, . 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 small-scale effects depending linearly on the second derivative along the rod axis of the axial load. This model is thus not able to evaluate small-scale 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 : 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 well-known differential equation (see, e.g., ) as follows: Note that the nonlocal contribution vanishes for constant loads . In the present paper, an alternative nonlocal constitutive behavior is adopted to assess small-scale effects in nanotubes also for constant axial loads. The corresponding axial displacement field is shown to be governed by a fourth-order 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 three-dimensional 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 time-rate, 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 gradient-type 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  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 so-called double stress tensor . 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
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 .
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 small-scale 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.
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 .
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.
- Q. Wang and K. M. Liew, “Application of nonlocal continuum mechanics to static analysis of micro- and nano-structures,” Physics Letters A, vol. 363, no. 3, pp. 236–242, 2007.
- 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.
- Ö. Civalek and Ç. Demir, “Bending analysis of microtubules using nonlocal Euler-Bernoulli beam theory,” Applied Mathematical Modelling, vol. 35, no. 5, pp. 2053–2067, 2011.
- M. A. De Rosa and C. Franciosi, “A simple approach to detect the nonlocal effects in the static analysis of Euler-Bernoulli and Timoshenko beams,” Mechanics Research Communications, vol. 48, pp. 66–69, 2013.
- 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.
- M. A. Kazemi-Lari, S. A. Fazelzadeh, and E. Ghavanloo, “Non-conservative instability of cantilever carbon nanotubes resting on viscoelastic foundation,” Physica E, vol. 44, no. 7-8, pp. 1623–1630, 2012.
- 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.
- 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.
- S. A. Emam, “A general nonlocal nonlinear model for buckling of nanobeams,” Applied Mathematical Modelling, vol. 37, no. 10-11, pp. 6929–6939, 2013.
- B. Fang, Y.-X. Zhen, C.-P. Zhang, and Y. Tang, “Nonlinear vibration analysis of double-walled carbon nanotubes based on nonlocal elasticity theory,” Applied Mathematical Modelling, vol. 37, no. 3, pp. 1096–1107, 2013.
- 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.
- 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.
- B. L. Wang and K. F. Wang, “Vibration analysis of embedded nanotubes using nonlocal continuum theory,” Composites B, vol. 47, pp. 96–101, 2013.
- 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.
- 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.
- 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.
- 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.
- M. A. Eltaher, A. E. Alshorbagy, and F. F. Mahmoud, “Vibration analysis of Euler-Bernoulli nanobeams by using finite element method,” Applied Mathematical Modelling, vol. 37, no. 7, pp. 4787–4797, 2013.
- 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.
- Y. Yan, W. Q. Wang, and L. X. Zhang, “Nonlocal effect on axially compressed buckling of triple-walled carbon nanotubes under temperature field,” Applied Mathematical Modelling, vol. 34, no. 11, pp. 3422–3429, 2010.
- R. Rafiee and R. M. Moghadam, “On the modeling of carbon nanotubes: a critical review,” Composites B, vol. 56, pp. 435–4490, 2014.
- J. N. Reddy, “Nonlocal theories for bending, buckling and vibration of beams,” International Journal of Engineering Science, vol. 45, no. 2-8, pp. 288–307, 2007.
- 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.
- F. Marotti de Sciarra, “Variational formulations, convergence and stability properties in nonlocal elastoplasticity,” International Journal of Solids and Structures, vol. 45, no. 7-8, pp. 2322–2354, 2008.
- 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.
- F. Marotti De Sciarra, “Hardening plasticity with nonlocal strain damage,” International Journal of Plasticity, vol. 34, pp. 114–138, 2012.
- G. Romano and R. Barretta, “Covariant hypo-elasticity,” European Journal of Mechanics A, vol. 30, no. 6, pp. 1012–1023, 2011.
- G. Romano and R. Barretta, “On Euler's stretching formula in continuum mechanics,” Acta Mechanica, vol. 224, no. 1, pp. 211–230, 2013.
- G. Romano and R. Barretta, “Geometric constitutive theory and frame invariance,” International Journal of Non-Linear Mechanics, vol. 51, pp. 75–86, 2013.
- G. Romano, R. Barretta, and M. Diaco, “Geometric continuum mechanics,” Meccanica, 2013.
- R. D. Mindlin, “Micro-structure in linear elasticity,” Archive for Rational Mechanics and Analysis, vol. 16, no. 1, pp. 51–78, 1964.
- F. Marotti de Sciarra, “Novel variational formulations for nonlocal plasticity,” International Journal of Plasticity, vol. 25, no. 2, pp. 302–331, 2009.
- F. Marotti de Sciarra, “On non-local and non-homogeneous elastic continua,” International Journal of Solids and Structures, vol. 46, no. 3-4, pp. 651–676, 2009.
- F. Marotti de Sciarra, “A nonlocal model with strain-based damage,” International Journal of Solids and Structures, vol. 46, no. 22-23, pp. 4107–4122, 2009.
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.