# Nonlocal Theory to Analyse Nanotube Structures under Tension

Authors: Hamidreza Yazdani Sarvestani; Ali Naghashpour
Journal: Journal of Nanomaterials
 (2015)

---

## Abstract

Because nanocomposites have found augmented use in many industries, the analytical solutions are required to be developed. This paper presents the development of a new analytical method for studying nanotube structures under tension using layer-wise and Eringen theories. Two opposite ends of tubes are subjected to normal forces. Nonlocal governing differential equations are derived and presented. The theoretical developments determine the effect of the geometric and nonlocal constitutive relations for single-walled nanotubes (SWNTs), double-walled nanotubes (DWNTs), and multiwalled nanotubes (MWNTs) under tension loading. It is observed that all displacement components increase with the increase in the nonlocal parameter.


---

## Body

## 1. Introduction

Carbon nanotubes (CNTs) have been used as multifunctional materials while analytical and numerical analyses are required to be performed for CNT structures. Naghashpour and Hoa [1, 2] dispersed CNTs into polymer composites to produce multifunctional polymer composites. Significant amount of experimental work on the use of CNTs in composites is being conducted while analytical and numerical analyses are required to be performed to obtain the properties and behaviors of nanocomposites.The vast majority of structural theories are derived using the constitutive assumptions that the stress at a point depends only on the strain at that point. Many investigations were made on the development of analytical solutions to analyse the composite structures [3, 4] while it remains to be observed whether these solutions can be applied to make analyses for nanostructures. Continuum based analyses have been generally utilized for the formulation of different nanostructures. Two main reasons are mentioned. One is experimental investigations in nanoscale that are difficult. The other is molecular dynamic simulations being highly computationally expensive for nanostructures. Over the past decade, some researchers have applied classical continuum mechanics such as Euler-Bernoulli theory, Timoshenko beam theory, and Kirchhoff’s plate theory to predict the behavior of nanostructures [5, 6].The nonlocal constitutive behaviour developed by Eringen [7–10] is based on the hypothesis that the stress at a point is a function of strains at all points in the continuum. The nonlocal elasticity theory was used for linear isotropic nonlocal solids by Eringen and Kim [11], in such a way that the nonlocal theory differs from the local theory in the stress-strain constitutive relations. Since the classical continuum elasticity is a scale-free theory, the use of classical continuum models may be not able to be relied on in the analysis of structural elements in nanoscale such as CNTs and graphene sheets. There are different modified classical continuum theories which capture size effects such as couple stress theory [12], strain gradient elasticity theory [13], modified couple stress theory [14, 15], and nonlocal elasticity theory [7–11]. Among all size-dependent theories, the nonlocal elasticity theory has been frequently used in the theoretical investigations of structures at small scale [16–21]. Heireche et al. [18] analysed the wave propagation in CNTs by developing a single-elastic beam model using nonlocal elasticity. By considering small-scale effect, the buckling behavior of nanoscale circular plates under uniform radial compression was investigated [20]. Mohammadi et al. [21] investigated the free vibration behavior of circular and annular graphene sheet by utilizing elasticity theory. They also derived the governing equations for single-layered graphene sheets (SLGS). Reddy [22–24] used different beam theories including those of Euler-Bernoulli, Timoshenko, and Levinson to analyze bending, buckling, and vibration of nonlocal beams. In his study, different displacement functions were selected in the first step and then all steps were redone when deriving beam equilibrium and equations of motion. Recently, a finite element method was developed for a nonlocal Timoshenko beam model [25]. The model was based on the key idea that nonlocal effects include long-range volume forces and moments exchanged by nonadjacent beam segments.It is found from the presented works that there is a lack of research for developing a modeling method that analyzes the mechanical response of nanotubes under extension. Even though a lot of researchers have considered nanotube as a sheet for their analysing, a tube structure is considered here to develop a displacement based theory to analysis of nanotubes. In the current work, we are studying the deformation of nanotube structures under axial forces based on the nonlocal elasticity theory in which the small-scale effects are taken into account. The displacement based governing equations of layer-wise theory are developed. Then, the nonlocal constitutive relations are used to express the stress resultants in terms of the generalized displacements. Finally, the equations are solved for specified axial force. The formulations show that the nonlocal parameter has prominent effect on the behavior of nanotube structures. It is anticipated that the results of the present work would be helpful for designing composite structures using SWNT or/and DWNT or/and MWNT. In addition, since direct experimental measurements are impractical due to the very small size of CNTs, the developed theoretical approach provides good alternative.


## 2. Theoretical Formulation

### 2.1. Displacement Field

SWNT, DWNT, and MWNT with mean radiusR and thickness h are subjected to a torque as shown in Figures 1, 2, and 3, respectively. The cylindrical coordinates (x,θ,r) are placed at the middle wall of the nanotubes so that x and r are the axial and radial coordinate, respectively. The integration of the appropriate linear strain-displacement relations of elasticity within cylindrical coordinate system will yield the following displacement field for the kth wall:(1a)u1kx,θ,r=xrC5kcos⁡θ+C4ksin⁡θ+C6kx+ukθ,r,(1b)u2kx,θ,r=xC1kcos⁡θ-C2ksin⁡θ-C3kr-12x2C4kcos⁡θ-C5ksin⁡θ+vkθ,r,(1c)u3kx,θ,r=xC1ksin⁡θ+C2kcos⁡θ-12x2C5kcos⁡θ+C4ksin⁡θ+wkθ,r,where u1kx,θ,r, u2kx,θ,r, and u3kx,θ,r represent the displacement components in the x, θ, and r directions, respectively, of a material point located at (x,θ,r) in the kth wall of the nanotube in Figures 1, 2, and 3. In order to satisfy the interfacial continuities of the displacement components, it is necessary for the integration constants appearing in (1a), (1b), and (1c) to be the same for all walls. Thus, relations (1a), (1b), and (1c) are presented as(2a)u1kx,θ,r=xrC5cos⁡θ+C4sin⁡θ+C6x+ukθ,r,(2b)u2kx,θ,r=xC1cos⁡θ-C2sin⁡θ-C3r-12x2C4cos⁡θ-C5sin⁡θ+vkθ,r,(2c)u3kx,θ,r=xC1sin⁡θ+C2cos⁡θ-12x2C5cos⁡θ+C4sin⁡θ+wkθ,r.Moreover, in (2a) ukθ,r can be replaced by -C1rsinθ+ukθ,r; it can be verified that the terms involving C1 in (2a), (2b), and (2c) correspond to an infinitesimal rigid-body rotation. These terms will, therefore, be ignored in the following developments since they will generate no strain. Similarly, it can be readily shown that the terms involving C2 must also be eliminated since they represent another rigid-body rotation of the tube. Furthermore, as long as the loading conditions at the two ends of the nanotube are identical, the constant C4 must vanish in order to satisfy the symmetry condition in deformation u3kx,θ,r=u3k-x,-θ,r. It is thus concluded that the most general form of the displacement field for the kth wall of a nanotube is given as(3a)u1kx,θ,r=C5xrcos⁡θ+C6x+ukθ,r,(3b)u2kx,θ,r=-C3rx+12C5x2sin⁡θ+vkθ,r,(3c)u3kx,θ,r=-12C5x2cos⁡θ+wkθ,r.Figure 1
Schematic sketch of SWNT under axial force and the location of coordinate axes.Figure 2
Schematic sketch of DWNT and the location of coordinate axes.Figure 3
Schematic sketch of MWNT and the location of coordinate axes.


### 2.2. Layer-Wise Theory (LWT)

Various theories such as the equivalent single-layer theories are unable to precisely represent the local phenomena in laminated composites. But then, the LWTs, which allow each layer of the laminate to act like a real three-dimensional layer, are able to present excellent results in definition of the localized phenomena. In LWT, the displacement components of a generic point in the laminate are conveniently given as(4a)u1x,θ,z=ukx,θΦkz,(4b)u2x,θ,z=vkx,θΦkz,(4c)u3x,θ,z=wkx,θΦkzk=1,2,…,N+1,with k, here and in what follows, being a dummy index implying summation of terms from k=1 to k=N+1. In (4a), (4b), and (4c), u1, u2, and u3 denote the displacement components in the x, θ, and r directions, respectively. Also, ukx,θ, vkx,θ, and wkx,θ represent the displacements of the points initially located on the kth wall of the nanotube in the x, θ, and r directions, respectively. Furthermore, Φkz is the global Lagrangian interpolation function that is used for the discretization of the displacement through thickness. The linear global interpolation function is defined as (5)ϕkz=0z≤zk-1ψk-12zzk-1≤z≤zkψk1zzk≤z≤zk+10z≥zk+1k=1,2,…,N+1,where ψkjj=1,2 are the local Lagrangian linear interpolation functions within the kth wall which are defined as(6)ψk1z=1hkzk+1-z,ψk2z=1hkz-zk,where hk is the thickness of the kth wall.Based on the reduced elasticity displacement field in (3a), (3b), and (3c), the LWT displacement field in (4a), (4b), and (4c) is rewritten as(7a)u1kx,θ,z=C5xR+zcos⁡θ+C6x+UkθΦkz,(7b)u2kx,θ,z=-C3R+zx+12C5x2sin⁡θ+VkθΦkz,(7c)u3kx,θ,z=-12C5x2cos⁡θ+WkθΦkz.By introducing r=R+z and considering thin-shell assumption (i.e., 1+z/R≈1), the strain-displacement relations are given as(8)εx=∂u1∂x,εθ=1R∂u2∂θ+u3R,γxθ=∂u2∂x+1R∂u1∂θ,γθz=1R∂u3∂θ+∂u2∂z-u2R,γxz=∂u3∂x+∂u1∂z,εz=∂u3∂z.Substitution of (7a), (7b), and (7c) into the strain-displacement relations (8) yields the following results:(9)εx=R+zC5cos⁡θ+C6,εθ=Vk′+WkΦkR,εz=WkΦk′,γθz=Wk′-VkΦkR+VkΦk′,γxz=UkΦk′,γxθ=-C3R+z+Uk′ΦkR.In (9) and what follows, a prime indicates an ordinary differentiation with respect to an appropriate variable (i.e., either θ or z). The equilibrium equations of a straight tube with N numerical walls are obtained by employing (9) in the principle of minimum total potential energy. The results are, in general, 3(N+1) local equilibrium equations corresponding to 3(N+1) unknown functions Uk, Vk, and Wk and three global equilibrium equations associated with the three parameters C3, C5, and C6.According to the principle of minimum total potential energy at the equilibrium configuration of a body the variation of the total potential energyΠ of the body must vanish. That is,(10)δΠ≡δU+δV=0,where δU is the variation of total strain energy of the body as follows:(11)δU=∫-aa∫-h/2h/2∫-ππσxδεx+σθδεθ+σzδεz+σθzδγθz+σxzδγxz+σxθδγxθdθdzdx.In (9), V is negative of the work done on the body by the specified external forces. Here, V=-2F0·u1x=±a,θ,z=-2F0C6a and, therefore, ∂V=-2F0∂C6a. Also, the variations of strains in (11) are found as(12)δεx=δC5R+zcos⁡θ+δC6,δεθ=δVk′+δWkΦkR,δεz=δWkΦk′,δγθz=δWk′-δVkΦkR+δVkΦk′,δγxz=δUkΦk′,δγxθ=-δC3R+z+δUk′ΦkR.Employing the fundamental lemma of calculus of variations the equilibrium equations and the associated boundary conditions of a laminate under tension are obtained as(13a)δUk:Qxk-1RdMxθkdθ=0,(13b)δVk:Qθk-1RdMθkdθ+Rθk=0,(13c)δWk:1RMθk-dRθkdθ+Nzk=0,(14a)δC3:∫-ππ∫-h/2h/2σxθR+zdzdθ=0,(14b)δC5:∫-ππ∫-h/2h/2σxR+zcos⁡θdzdθ=0,(14c)δC6:∫-ππ∫-h/2h/2Rσxdzdθ=F0,where F0 indicates the prescribed value of axial force applied at the two ends of the nanotube. Also, the following traction-free boundary conditions must be satisfied: (15)Rθk=Qxk=Nzk=0atz=±h2,where the generalized stress and moment resultants are defined as(16a)Nzk,Qxk,Qθk=∫-h/2h/2σz,σxz,σθzΦk′dz,(16b)Mθk,Mxθk,Rθk=∫-h/2h/2σθ,σxθ,σθzΦkdz.It is to be noted that, in (15) and (16a) and (16b), the superscript k refers to the kth wall in the nanotube.


## 2.1. Displacement Field

SWNT, DWNT, and MWNT with mean radiusR and thickness h are subjected to a torque as shown in Figures 1, 2, and 3, respectively. The cylindrical coordinates (x,θ,r) are placed at the middle wall of the nanotubes so that x and r are the axial and radial coordinate, respectively. The integration of the appropriate linear strain-displacement relations of elasticity within cylindrical coordinate system will yield the following displacement field for the kth wall:(1a)u1kx,θ,r=xrC5kcos⁡θ+C4ksin⁡θ+C6kx+ukθ,r,(1b)u2kx,θ,r=xC1kcos⁡θ-C2ksin⁡θ-C3kr-12x2C4kcos⁡θ-C5ksin⁡θ+vkθ,r,(1c)u3kx,θ,r=xC1ksin⁡θ+C2kcos⁡θ-12x2C5kcos⁡θ+C4ksin⁡θ+wkθ,r,where u1kx,θ,r, u2kx,θ,r, and u3kx,θ,r represent the displacement components in the x, θ, and r directions, respectively, of a material point located at (x,θ,r) in the kth wall of the nanotube in Figures 1, 2, and 3. In order to satisfy the interfacial continuities of the displacement components, it is necessary for the integration constants appearing in (1a), (1b), and (1c) to be the same for all walls. Thus, relations (1a), (1b), and (1c) are presented as(2a)u1kx,θ,r=xrC5cos⁡θ+C4sin⁡θ+C6x+ukθ,r,(2b)u2kx,θ,r=xC1cos⁡θ-C2sin⁡θ-C3r-12x2C4cos⁡θ-C5sin⁡θ+vkθ,r,(2c)u3kx,θ,r=xC1sin⁡θ+C2cos⁡θ-12x2C5cos⁡θ+C4sin⁡θ+wkθ,r.Moreover, in (2a) ukθ,r can be replaced by -C1rsinθ+ukθ,r; it can be verified that the terms involving C1 in (2a), (2b), and (2c) correspond to an infinitesimal rigid-body rotation. These terms will, therefore, be ignored in the following developments since they will generate no strain. Similarly, it can be readily shown that the terms involving C2 must also be eliminated since they represent another rigid-body rotation of the tube. Furthermore, as long as the loading conditions at the two ends of the nanotube are identical, the constant C4 must vanish in order to satisfy the symmetry condition in deformation u3kx,θ,r=u3k-x,-θ,r. It is thus concluded that the most general form of the displacement field for the kth wall of a nanotube is given as(3a)u1kx,θ,r=C5xrcos⁡θ+C6x+ukθ,r,(3b)u2kx,θ,r=-C3rx+12C5x2sin⁡θ+vkθ,r,(3c)u3kx,θ,r=-12C5x2cos⁡θ+wkθ,r.Figure 1
Schematic sketch of SWNT under axial force and the location of coordinate axes.Figure 2
Schematic sketch of DWNT and the location of coordinate axes.Figure 3
Schematic sketch of MWNT and the location of coordinate axes.


## 2.2. Layer-Wise Theory (LWT)

Various theories such as the equivalent single-layer theories are unable to precisely represent the local phenomena in laminated composites. But then, the LWTs, which allow each layer of the laminate to act like a real three-dimensional layer, are able to present excellent results in definition of the localized phenomena. In LWT, the displacement components of a generic point in the laminate are conveniently given as(4a)u1x,θ,z=ukx,θΦkz,(4b)u2x,θ,z=vkx,θΦkz,(4c)u3x,θ,z=wkx,θΦkzk=1,2,…,N+1,with k, here and in what follows, being a dummy index implying summation of terms from k=1 to k=N+1. In (4a), (4b), and (4c), u1, u2, and u3 denote the displacement components in the x, θ, and r directions, respectively. Also, ukx,θ, vkx,θ, and wkx,θ represent the displacements of the points initially located on the kth wall of the nanotube in the x, θ, and r directions, respectively. Furthermore, Φkz is the global Lagrangian interpolation function that is used for the discretization of the displacement through thickness. The linear global interpolation function is defined as (5)ϕkz=0z≤zk-1ψk-12zzk-1≤z≤zkψk1zzk≤z≤zk+10z≥zk+1k=1,2,…,N+1,where ψkjj=1,2 are the local Lagrangian linear interpolation functions within the kth wall which are defined as(6)ψk1z=1hkzk+1-z,ψk2z=1hkz-zk,where hk is the thickness of the kth wall.Based on the reduced elasticity displacement field in (3a), (3b), and (3c), the LWT displacement field in (4a), (4b), and (4c) is rewritten as(7a)u1kx,θ,z=C5xR+zcos⁡θ+C6x+UkθΦkz,(7b)u2kx,θ,z=-C3R+zx+12C5x2sin⁡θ+VkθΦkz,(7c)u3kx,θ,z=-12C5x2cos⁡θ+WkθΦkz.By introducing r=R+z and considering thin-shell assumption (i.e., 1+z/R≈1), the strain-displacement relations are given as(8)εx=∂u1∂x,εθ=1R∂u2∂θ+u3R,γxθ=∂u2∂x+1R∂u1∂θ,γθz=1R∂u3∂θ+∂u2∂z-u2R,γxz=∂u3∂x+∂u1∂z,εz=∂u3∂z.Substitution of (7a), (7b), and (7c) into the strain-displacement relations (8) yields the following results:(9)εx=R+zC5cos⁡θ+C6,εθ=Vk′+WkΦkR,εz=WkΦk′,γθz=Wk′-VkΦkR+VkΦk′,γxz=UkΦk′,γxθ=-C3R+z+Uk′ΦkR.In (9) and what follows, a prime indicates an ordinary differentiation with respect to an appropriate variable (i.e., either θ or z). The equilibrium equations of a straight tube with N numerical walls are obtained by employing (9) in the principle of minimum total potential energy. The results are, in general, 3(N+1) local equilibrium equations corresponding to 3(N+1) unknown functions Uk, Vk, and Wk and three global equilibrium equations associated with the three parameters C3, C5, and C6.According to the principle of minimum total potential energy at the equilibrium configuration of a body the variation of the total potential energyΠ of the body must vanish. That is,(10)δΠ≡δU+δV=0,where δU is the variation of total strain energy of the body as follows:(11)δU=∫-aa∫-h/2h/2∫-ππσxδεx+σθδεθ+σzδεz+σθzδγθz+σxzδγxz+σxθδγxθdθdzdx.In (9), V is negative of the work done on the body by the specified external forces. Here, V=-2F0·u1x=±a,θ,z=-2F0C6a and, therefore, ∂V=-2F0∂C6a. Also, the variations of strains in (11) are found as(12)δεx=δC5R+zcos⁡θ+δC6,δεθ=δVk′+δWkΦkR,δεz=δWkΦk′,δγθz=δWk′-δVkΦkR+δVkΦk′,δγxz=δUkΦk′,δγxθ=-δC3R+z+δUk′ΦkR.Employing the fundamental lemma of calculus of variations the equilibrium equations and the associated boundary conditions of a laminate under tension are obtained as(13a)δUk:Qxk-1RdMxθkdθ=0,(13b)δVk:Qθk-1RdMθkdθ+Rθk=0,(13c)δWk:1RMθk-dRθkdθ+Nzk=0,(14a)δC3:∫-ππ∫-h/2h/2σxθR+zdzdθ=0,(14b)δC5:∫-ππ∫-h/2h/2σxR+zcos⁡θdzdθ=0,(14c)δC6:∫-ππ∫-h/2h/2Rσxdzdθ=F0,where F0 indicates the prescribed value of axial force applied at the two ends of the nanotube. Also, the following traction-free boundary conditions must be satisfied: (15)Rθk=Qxk=Nzk=0atz=±h2,where the generalized stress and moment resultants are defined as(16a)Nzk,Qxk,Qθk=∫-h/2h/2σz,σxz,σθzΦk′dz,(16b)Mθk,Mxθk,Rθk=∫-h/2h/2σθ,σxθ,σθzΦkdz.It is to be noted that, in (15) and (16a) and (16b), the superscript k refers to the kth wall in the nanotube.


## 3. Nonlocal Theory (Constitutive Relations)

According to Eringen [7–9], the stress field at a point x in an elastic continuum depends not only on the strain field at the point but also on strains at all other points of the body. Eringen applied this fact to the atomic theory of lattice dynamics and experimental observations on phonon dispersion. The integral constitutive relations in an equivalent differential form are presented based on Eringen [7–9] as(17)1-μ∇2σij=Cijmnεmn,where ∇2 and μ=e0a2 are the Laplacian operator and nonlocal parameter, respectively. In addition, e0 and a are a material constant and internal characteristic length, respectively. From both consideration of (17) and Hook’s law, it is written that(18)1-μ∇2σxnlσθnlσznlσθznlσxznlσxθnlk=C-11C-12C-1300C-16C-12C-22C-2300C-26C-13C-23C-3300C-36000C-44C-450000C-45C-550C-16C-26C-3600C-66kεxεθεzγθzγxzγxθk,where C-ijkrepresent the off-axis stiffnesses. Multiplying the relation (18) by Φk′ and Φk and integrating them yield(19a)∫-h/2h/2σθσθzσxθkΦkdz-μ∇2∫-h/2h/2σθσθzσxθkΦkdz=C-12C-22C-2300C-26000C-44C-450C-16C-26C-3600C-66k·∫-h/2h/2εxεθεzγθzγxzγxθkΦkdz,(19b)∫-h/2h/2σzσθzσxzkΦk′dz-μ∇2∫-h/2h/2σzσθzσxzkΦk′dz=C-13C-23C-3300C-36000C-44C-450000C-44C-450k·∫-h/2h/2εxεθεzγθzγxzγxθkΦk′dz.By considering (16a) and (16b) and strain displacement relations (Equation (9)), the stress resultants are given as(20)1-μ∇2Nzk,Mθk,Mxθk=B36jk,D26kj,D66kjUj′R+B23jk,D22kj,D26kjVj′R+A33kj,B23kj,B36kjWj+B23jk,D22kj,D26kjWjR-A36k,B26k,B66kRC3-B-36k,D26k,D66kC3+A13k,B12k,B16kRC5cos⁡θ+B-13k,D12k,D16kC5cos⁡θ-A13k,B12k,B16kC6,1-μ∇2Qxk,Qθk,Rθk=A55kj,A45kj,B45kjUj+A45kj,A44kj,B44kjVj-B45jk,B44jk,D44kjVjR+B45jk,B44jk,D44kjWj′R,where the wall rigidities in (20) are defined as(21)Apqkj,B-pqkj,Dpqkj=∑i=1N∫zizi+1C-pqiϕk′ϕj′,ϕkϕj′,ϕkϕjdz,Apqk,Bpqk,B-pqk,Dpqk=∑i=1N∫zizi+1C-pqiϕk',ϕk,ϕk'z,ϕkzdz,k,j=1,2,…,N+1.The nonlocal displacement equilibrium equations within LWT are obtained by substituting (20) into (13a), (13b), and (13c):(22a)δUk:D66kjR21-λ2/z2μUj′′-A55kj1-λ2/z2μUj+D26kjR21-λ2/z2μVj′′-A45kj-B45jk/R1-λ2/z2μVj+D26kj/R2+B36kj-B45jk/R1-λ2/z2μWj′=B16k+D16k/R1+μ/z2C5sin⁡θ+zkRμ0.5e-zk/μ0.5θ-ezk/μ0.5θ,(22b)δVk:D26kjR21-λ2/z2μUj′′-A45kj-B45kj/R1-λ2/z2μUj+D22kjR21-λ2/z2μVj′′-A44kj-B44kj+B44jk/R+D44kj/R21-λ2/z2μVj+B23kj-B44jk/R+D44kj+D22kj/R21-λ2/z2μWj′=B12k+D12k/R1+μ/z2C5sin⁡θ,(22c)δWk:B45kj-B36jk/R-D26kj/R21-λ2/z2μUj′+B44kj-B23jk/R-D44kj+D22kj/R21-λ2/z2μVj′+D44kjR21-λ2/z2μWj′′+D22kj/R2+B23jk+B23kj/R+A33kj1-λ2/z2μWj=-B-36k+B26k+RA36k+D26kRC3+A13k+B12kRC6+B-13k+B12k+RA13k+D12k/R1+μ/z2C5cos⁡θ+1Re-zk/μ0.5θ+ezk/μ0.5θk=1,2,…,N+1.Also, the global equilibrium equations of the nanotube are expressed in terms of displacement functions by substituting (9) into (19a), (19b) and the subsequent results into (14a), (14b), and (14c).


## 4. Analytical Solution for Nanotube Structure

The system in (22a), (22b), and (22c) shows 3(N+1) coupled ordinary differential equations with constant coefficients which may be displayed in a matrix form as(23)Mη′′+Kη=F1C+F2,where(24a)η=UT,VT,W-TT,U=U1,U2,…,UN+1T,C=C3,C5,C6T,V=V1,V2,…,VN+1T,W-=W-1,W-2,…,W-N+1T,(24b)W-j=∫θWjdθ.The coefficient matrices [M], [K], [F1], and [F2] in (23) are defined in the Appendix. It can readily be confirmed that the general solution of (23) may be presented as(25)η=ψsinh⁡λθk+K-1F1C+K-1F2and [sinh(λθ)] is a 3(N+1)×3(N+1) diagonal matrix. That is,(26)
sinh
⁡
λ
θ
=
diag
⁡
sinh
⁡
λ
1
θ
,
sinh
⁡
λ
2
θ
,
…
,
sinh
⁡
λ
3
N
+
1
θ
.Also, [ψ] and λ12,λ22,…,λ3(N+1)2 are the model matrix and eigenvalues of (-M-1[K]), respectively. Matrix k is an unknown vector representing 3(N+1) integration constants. The constants Cj(j=3,5,6) must be calculated within LWT analysis. Therefore, the boundary conditions in (15) are first imposed to find the vector k in terms of the unknown parameters Cj(j=3,5,6). These constants are then calculated in terms of the specified axial force F0 by satisfaction of the global equilibrium conditions in (14a), (14b), and (14c).


## 5. Results and Discussion

The numerical results are discussed for single-walled carbon nanotube (SWCNT), double-walled carbon nanotube (DWCNT), and multiwalled carbon nanotube (MWCNT) under extension. The mechanical properties of SWCNT, DWCNT, and MWCNT are given in Table1 [26–28]. Furthermore, the displacement components are normalized as u-i=ui/u0, where u0=L∗σ0/E and σ0=(F0)/(π/4∗(OD2-ID2). An innermost radius of 8.5 nm is assumed for SWCNT, DWCNT, and MWCNT while assuming thickness of SWCNT (heff) to be 0.34 nm [29]. Also, L is the length of nanotubes. A 10-layer MWCNT is considered here in which R/(N∗heff)=2. L/h is defined as an aspect ratio.Table 1
Mechanical properties of SWCNT, DWCNT, and MWCNT [26–28].
Nanotubes
Mechanical properties
E (TPa)
v
SWCNT
1.80
0.19
DWCNT
1.33
0.19
MWCNT
0.95
0.19The nanotube lengths, thickness, and nonlocal parameters play significant roles on determining the nondimensional displacement components. These parameters are studied numerically.The numerical results for SWCNT, DWCNT, and MWCNT under axial force (F0) at the point in which θ=270° are presented in Tables 2, 3, and 4. Table 2 shows the effects of length, thickness, and nonlocal parameter (μ) on the nondimensional longitudinal displacement component (u-x). It is found from Table 2 that the longitudinal displacement component (u-x) increases with increasing the thickness while CNT length is kept constant. For SWCNT, the longitudinal displacement component (u-x) increases as the length of SWCNT increases. Observing Table 2, the longitudinal displacement component increases while the length and thickness of CNT are kept constant.Table 2
Comparison of nondimensional displacement component in the longitudinal direction (u-x).
L
/
h
µ
SWCNT
DWCNT
MWCNT (×101)
100
0
1.082
1.276
2.170
1.0
2.014
2.376
4.040
2.0
2.052
2.421
4.116
3.0
2.201
2.597
4.415
4.0
2.657
3.135
5.329
5.0
3.021
3.564
6.060
200
0
1.136
1.340
2.279
1.0
2.114
2.495
4.242
2.0
2.154
2.542
4.322
3.0
2.311
2.727
4.635
4.0
2.789
3.292
5.596
5.0
3.172
3.743
6.363
500
0
1.150
1.358
2.308
1.0
2.142
2.527
4.297
2.0
2.182
2.575
4.378
3.0
2.341
2.762
4.696
4.0
2.826
3.334
5.669
5.0
3.213
3.791
6.445Table 3
Comparison of nondimensional displacement component in the circumferential direction (u-θ).
L
/
h
µ
SWCNT (×10−2)
DWCNT (×10−2)
MWCNT (×10−1)
100
0
1.110
1.309
2.226
1.0
2.131
2.514
4.274
2.0
2.265
2.672
4.543
3.0
2.579
3.043
5.173
4.0
3.001
3.540
6.018
5.0
3.260
3.846
6.539
200
0
1.165
1.375
2.337
1.0
2.237
2.640
4.488
2.0
2.378
2.806
4.770
3.0
2.707
3.195
5.432
4.0
3.150
3.717
6.318
5.0
3.423
4.039
6.866
500
0
1.180
1.393
2.368
1.0
2.266
2.674
4.546
2.0
2.409
2.842
4.832
3.0
2.743
3.236
5.502
4.0
3.190
3.765
6.401
5.0
3.467
4.091
6.955Table 4
Comparison of nondimensional displacement component in the radial direction (u-z).
L
/
h
µ
SWCNT (×10−4)
DWCNT (×10−4)
MWCNT (×10−3)
100
0
1.354
1.597
2.716
1.0
1.887
2.226
3.785
2.0
2.265
2.672
4.543
3.0
2.402
2.834
4.818
4.0
3.012
3.554
6.042
5.0
3.223
3.803
6.465
200
0
1.421
1.677
2.851
1.0
1.981
2.337
3.974
2.0
2.378
2.806
4.770
3.0
2.522
2.976
5.059
4.0
3.162
3.731
6.344
5.0
3.384
3.993
6.788
500
0
1.440
1.699
2.889
1.0
2.007
2.368
4.026
2.0
2.409
2.842
4.832
3.0
2.554
3.014
5.125
4.0
3.203
3.780
6.426
5.0
3.428
4.045
6.876The effects of length, thickness, and nonlocal parameter (μ) on nondimensional circumferential displacement components (u-θ) are studied in Table 3. It is observed from Table 3 that the circumferential displacement component (u-θ) increases with increasing the thickness while CNT length is kept constant. The circumferential displacement component (u-θ) increases with increasing the length. Observing Table 3, the circumferential displacement component (u-θ) increases as the nonlocal parameter (μ) increases where length and thickness of CNT are fixed.Table4 reveals investigation for the effects of length, thickness, and nonlocal parameter (μ) on nondimensional radial displacement components (u-z). Table 4 shows the radial displacement component (u-z) increases as the thickness of nanotube increases. For the case of constant thickness, the radial displacement component (u-z) increases as CNT length increases. It is clear from Table 4 that the increase of nonlocal parameter (μ) causes increase of the radial displacement component (u-z).


## 6. Conclusions

We have developed a new nonlocal method for nanotube structures subjected to tension. Equations of motion using layer-wise theory were derived based on Eringen’s differential constitutive equations of nonlocal elasticity. Variational statements of the theory were also presented to facilitate direct development of the nonlocal displacement models of layer-wise theory. The equations of motion were then analytically solved for nanotubes under axial force to bring out the effect of nonlocal parameter. All displacement components increase with the increase in the nonlocal parametersμ. Remark that the displacement components increase with increasing the length and thickness of CNT. Also, it is observed that the displacement components for MWCNT and DWCNT are greater than those for SWCNT. The values for displacement component in the longitudinal direction (u-x) are greater than those for the displacement components in the circumferential (u-θ) and radial (u-z) directions for SWCNT, DWCNT, and MWCNT.



---

*Source: 635287-2015-06-22.xml*