# Rayleigh-Ritz Vibrational Analysis of Multiwalled Carbon Nanotubes Based on the Nonlocal Flügge Shell Theory

Authors: H. Rouhi; M. Bazdid-Vahdati; R. Ansari
Journal: Journal of Composites
 (2015)

---

## Abstract

A nonlocal elastic shell model considering the small scale effects is developed to study the free vibrations of multiwalled carbon nanotubes subject to different types of boundary conditions. Based on the nonlocal elasticity and the Flügge shell theory, the governing equations are derived which include the interaction of van der Waals forces between adjacent and nonadjacent layers. To analytically solve the problem, the Rayleigh-Ritz method is employed. In the present analysis, different combinations of layerwise boundary conditions are taken into account. Some new intertube resonant frequencies and the associated noncoaxial vibrational modes are identified owing to incorporating circumferential modes into the shell model.


---

## Body

## 1. Introduction

Nanotechnology is quickly becoming one of the fastest developing fields of research. This is largely due to numerous applications of carbon-based nanostructures in different fields. Among such nanostructures, carbon nanotubes (CNTs) discovered by Iijima [1] in 1991 have attracted the attention of researchers worldwide [2–4].The theoretical models based on continuum mechanics can be efficiently employed for the analysis of nanostructures. Applications of the classical (local) continuum modeling to the study of CNTs have been suggested in several papers [5–11]. Because the local continuum models are incapable of capturing the influence of small scale of nanostructures, they become controversial to implement for the analysis of such nanoscale systems. So as to consider the small size effects, some other researchers have proposed the application of nonlocal continuum mechanics to study nanostructured materials.The use of nonlocal elasticity to nanotechnology was first proposed by Peddieson et al. [12]. After that, nonlocal continuum models have been utilized by many researchers owing to their computational efficiency as well as their accuracy [13–22]. The vibrational analysis of CNTs based on the nonlocal continuum mechanics has been the focus of considerable research. Wang and Varadan [13] investigated the vibration of single-walled and double-walled CNTs (SWCNTs) based on a nonlocal beam model. It was shown that the nonlocal results are in good agreement with the experimental results. Li and Kardomateas [14] investigated the vibrations of multiwalled carbon nanotubes (MWCNTs) subject to simply supported boundary conditions via a nonlocal shell model. It was concluded that the beam models may be insufficient for analyzing the dynamics of nanotubes due to not considering the circumferential mode. Narendar and Gopalakrishnan [16] presented the influence of nonlocal scaling parameter on the wave propagation in MWCNTs using the nonlocal Timoshenko beam model. Arash and Ansari [20] utilized a nonlocal shell model to study the free vibration behavior of SWCNTs. The results indicated that although uncertainty exists in defining nanotube wall thickness, their developed nonlocal shell model is able to predict the results obtained from the molecular dynamics (MD) simulations.Using shell models in the analysis of CNTs has attracted the attention of several researchers. There are different shell theories. The reader is referred to [23] for a comprehensive summary and discussion on shell theories. A pioneering work in this field is related to Love’s paper [24]. Using the first approximation theory of Love, Flügge shell theory [25] was established. The Flügge shell theory is according to the Kirchhoff-Love assumption for thin elastic shells. By means of this theory, the kinematical relations of the middle surface of a cylindrical shell can be derived. Donnell’s simplified theory can be also obtained by neglecting few terms in the equations of Flügge shell theory.The Flügge shell theory is very accurate for thin shells. In addition, CNT is in fact a crystalline membrane with one atom thickness. From this viewpoint, the Flügge shell theory assumptions made in representing a CNT as a homogeneous thin shell can be satisfactory. Also, the Flügge shell theory is frequently used for the analysis of CNTs due to the relatively accurate results in spite of its theoretical simplicity. A marked increase in the volume of research papers dealing with the application of Flügge shell theory to different problems of CNTs stands as testimony to this claim [26–29]. In [27], the axisymmetric and beamlike vibrations of MWCNTs under simply supported boundary conditions were studied using a Flügge-type shell model.The majority of the previous studies are limited to CNTs with all edges simply supported for which exact analytical solution exists. The boundary conditions can significantly affect the mechanical behavior of nanostructures. In addition, for MWCNTs, different combinations of layerwise boundary conditions are possible to be considered. In this regard, Xu et al. [30] commented that “the relevance of the existing model in which both tubes have the same boundary conditions for the vibration of double-walled CNTs is questionable.”In this work, a nonlocal Flügge shell model is developed to study the vibrations of MWCNTs with different boundary conditions in a layerwise manner. The analytical Rayleigh-Ritz method is applied to the variational statement derived from the nonlocal Flügge-type vibration equations. Some new noncoaxial intertube resonances for MWCNTs are also predicted in this paper which is a matter of great technical importance in the study of vibrations of nested nanotubes.


## 2. Elastic Shell Model in Nonlocal Elasticity

Based on Eringen’s nonlocal theory [31, 32], the nonlocal stress tensor σ at point x is given by(1)σx=∫Vλx-x′,μtx′dVx′,∀x∈V,where t denotes the macroscopic stress tensor and λx-x′,μ is the nonlocal modulus or attenuation function whose arguments are the Euclidean distance x-x′ and a material constant μ=e0a/l with a as an internal characteristic length, l as an external characteristic length (e.g., crack length, wavelength), and e0 as material constant. The differential form of (1) is(2)1-e0a2∇2σ=t,where e0a is the nonlocal parameter. Hooke’s generalized law states that(3)t=S:ϵin which S is the fourth-order elasticity tensor and “:” denotes the double dot product. Using (3) leads to(4a)σxx-e0a2∇2σxx=E1-ν2εxx+νεθθ,(4b)σθθ-e0a2∇2σθθ=E1-ν2εθθ+νεxx,(4c)σxθ-e0a2∇2σxθ=E21+νγxθ,where E is Young’s modulus of the material and ν is the Poisson ratio. x and θ are longitudinal and angular circumferential coordinates. Each tube of MWCNTs is described as an individual cylindrical shell of radius R, length L, and thickness h, as shown in Figure 1. According to the classic shell theory, the three-dimensional displacement components ux, uθ, and uz in x, θ, and z directions, respectively, are of the form(5a)uxx,θ,z,t=ux,θ,t-z∂w∂xx,θ,t,(5b)uθx,θ,z,t=vx,θ,t-z∂w∂θx,θ,t,(5c)uzx,θ,z,t=wx,θ,z,t,where u, v, and w are the reference surface displacements. The kinematical relations are given by(6)
ε
x
x
ε
θ
θ
γ
x
θ
=
∂
u
∂
x
1
R
∂
v
∂
θ
+
w
R
∂
v
∂
x
+
1
R
∂
u
∂
θ
+
z
-
∂
2
w
∂
x
2
-
1
R
2
∂
2
w
∂
θ
2
+
w
-
2
R
∂
2
w
∂
x
∂
θ
+
1
R
∂
v
∂
x
-
1
R
2
∂
u
∂
θ
.The stress and moment resultants can be obtained by(7a)N=Nxx,Nθθ,NxθTdz=∫-h/2h/2σxx,σθθ,σxθTdz,(7b)M=Mxx,Mθθ,MxθTdz=∫-h/2h/2σxx,σθθ,σxθTzdz,where h is the shell thickness.Figure 1
Schematic of an MWCNT.Moreover, the stress and moment resultants are obtained as [25](8)
N
x
x
N
θ
θ
N
x
θ
M
x
x
M
θ
θ
M
x
θ
-
e
0
a
2
∇
2
N
x
x
N
θ
θ
N
x
θ
M
x
x
M
θ
θ
M
x
θ
=
E
h
1
-
ν
2
∂
∂
x
ν
E
h
1
-
ν
2
1
R
∂
∂
θ
ν
E
h
1
-
ν
2
1
R
ν
E
h
1
-
ν
2
∂
∂
x
E
h
1
-
ν
2
1
R
∂
∂
θ
E
h
1
-
ν
2
1
R
E
h
2
1
+
ν
1
R
∂
∂
θ
E
h
2
1
+
ν
∂
∂
x
0
0
D
ν
R
2
∂
∂
θ
-
D
∂
2
∂
x
2
+
ν
R
2
∂
2
∂
θ
2
0
D
1
R
2
∂
∂
θ
-
D
ν
∂
2
∂
x
2
+
1
R
2
∂
2
∂
θ
2
0
D
R
1
-
ν
∂
∂
x
-
D
R
1
-
ν
∂
2
∂
x
∂
θ
u
v
win which D is the bending rigidity.The governing equations on the basis of the Flügge shell theory are given as [25](9a)∂Nxx∂x+1R∂Nθx∂θ=ρh∂2u∂t2,(9b)1R∂Nθθ∂θ+∂Nxθ∂x+1R2∂Mθθ∂θ+1R∂Mxθ∂x=ρh∂2v∂t2,(9c)∂2Mxx∂x2+1R2∂2Mθθ∂θ2+2R∂2Mxθ∂x∂θ-NθθR+p=ρh∂2w∂t2,where p is the pressure exerted on tube i through the van der Waals (vdW) interaction forces, as indicated in Figure 2. The vdW model is formulated as [33](10)pvdWi=wi∑j=1Ncij-∑j=1Ncijwj,i=1,2,…,N,in which vdW coefficients cij representing the pressure increment contributing to layer i from layer j are given by(11)cij=1001πεσ123a4Eij13-1120πεσ69a4Eij7Rj,where a=1.42Å is the C-C bond length, ε is the depth of the potential, σ is a parameter that is determined by the equilibrium distance and Rj is the radius of jth layer, and Eijm denotes the elliptic integral defined as(12)Eijm=Rj+Ri-m∫0π/2dθ1-Kijcos2⁡θm/2;here m is an integer and the coefficient Kij is given by(13)Kij=4RjRiRj+Ri2.Figure 2
Cross-sectional view of an MWCNT under the vdW interactions.


## 3. Field Equations

Forith layer of an MWCNT, by the use of (8), (9a), (9b), and (9c) can be stated in terms of the three field variables u(i),v(i),w(i),(i=1,2,…,N):(14a)l11iui+l12ivi+l13iwi=ρhu¨i-e0a2ρhu¨xxi+1Ri2u¨θθi,(14b)l21iui+l22ivi+l23iwi=ρhv¨i-e0a2ρhv¨xxi+1Ri2v¨θθi,(14c)l31iui+l32ivi+l33iwi=ρhw¨i+wi∑j=1j≠iNcij-∑j=1j≠iNcijwj-e0a2ρhw¨xxi+1Ri2w¨θθi+wxxi+1Ri2wθθi∑j=1j≠iNcij-∑j=1j≠iNcijwxxj+1Ri2wθθj,where lpq(i)(p,q=1,2,3) are the partial operators, which are given in the Appendix.


## 4. Solution Using the Rayleigh-Ritz Method

### 4.1. Variational Form of the Field Equations

According to the semi-inverse method [34–36], variational trial-functional Πu1,u2,…,uN,v1,v2,…,vN,w1,w2,…,wN is constructed as(15)Π=∑i=1NΠiui,vi,w1,w2,…,wN,where(16)Πiui,vi,w1,w2,…,wN=ΠKiui,vi,wi+ΠMiui,vi,wi+ΠvdWiw1,w2,…,wNin which(17)ΠKiui,vi,wi=12∫0T∫ΩEh1-ν2u,xi+1RiνEh1-ν2v,θi+wiu,xi+Eh21+νu,θiRi+v,xi2+νEh1-ν2u,xi+1RiEh1-ν2v,θi+wi1Riv,θi+wi+Dw,xxi+νRi2w,θθiw,xxi+DRi21Ri2w,θθi+νw,xxiw,θθi+21-νDRi2w,xθi-v,xi2-2DRi2νw,xxi+1Ri2w,θθiv,θidΩdt,ΠMiui,vi,wi=12∫0T∫Ωρhu˙i2-e0a2u˙,xxi+1Ri2u˙,θθiu˙i+v˙i2-e0a2v˙,xxi+1Ri2v˙,θθiv˙i+w˙i2-e0a2w˙,xxi+1Ri2w˙,θθiw˙idΩdt,ΠvdWiw1,w2,…,wN=12∑j=1N∫0T∫Ωcijwi-wj2-e0a2wiw,xxi+1Ri2w,θθi+wjw,xxj+1Ri2w,θθjdΩdt.


### 4.2. Implementation of the Rayleigh-Ritz Method with Beam Functions

The field variables corresponding to theith tube, that is, ui, vi, and wi, can be expressed in the form(18a)uix,θ,t=Ai∂ϕx∂xcos⁡nθsin⁡ωt,(18b)vix,θ,t=Biϕxsin⁡nθsin⁡ωt,(18c)wix,θ,t=Ciϕxcos⁡nθsin⁡ωt,where Ai, Bi, and Ci are the constant parameters denoting the amplitudes of the vibrations, n is the circumferential wave number, ω is the angular frequency of vibrating CNT, and ϕx is the axial function which satisfies the geometric boundary conditions. The axial function ϕ(x) is chosen as the characteristics beam function as(19)ϕx=α1cosh⁡λmxL+α2cos⁡λmxL-ςmα3sinh⁡λmxL+α4sinλmxLin which αii=1,…,4 are constants with value 0, 1, or -1 depending on the tube ends, λm is the roots of the transcendental equations obtained from the CNT boundary conditions, and ςm is the parameters corresponding to λm. The geometric boundary conditions for clamped (C), free (F), and simply supported (S) boundary conditions can be expressed mathematically in terms of ϕx as follows.Clamped boundary condition is(20)ϕx=ϕ′x=0.Free boundary condition is(21)ϕ′′x=ϕ′′′x=0.Simply supported boundary condition is(22)ϕx=ϕ′′x=0.Parametersαi, λm, and ςm that are chosen according to the CNT boundary conditions are given in Table 1. Substituting (18a), (18b), and (18c) into (15) and then minimizing the energy functional Π with respect to unknown coefficients Ai, Bi, and Ci result in the following algebraic equations:(23)∂Π∂Ai=0,∂Π∂Bi=0,∂Π∂Ci=0,i=1,2,…,N.The above equations can be recast in the form of a generalized eigenvalue problem as(24)K-ω2MX=0,where K is the stiffness matrix, M is the mass matrix of the CNT, and X is given by(25)XT=A1,B1,C1,A2,B2,C2,…,AN,BN,CN.Table 1
Values ofαi, λm, and ζm for SS, CC, FF, CS, CF, and FS boundary conditions.
Boundary condition
α
i  i=1,…,4
Values forλm
ζ
m
SS
α
1
=
0
,
α
2
=
0 α3=0,α4=-1
m
π
1
CC
α
1
=
1
,
α
2
=
-
1 α3=1,α4=-1
cos
⁡
λ
m
⁡
cosh
⁡
λ
m
=
1
cosh
⁡
λ
m
⁡
-
cos
⁡
λ
m
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
m
FF
α
1
=
1
,
α
2
=
1 α3=1,α4=1
cos
⁡
λ
m
⁡
cosh
⁡
λ
m
=
1
cosh
⁡
λ
m
⁡
-
cos
⁡
λ
m
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
m
CS
α
1
=
1
,
α
2
=
-
1 α3=1,α4=-1
tan
⁡
λ
m
⁡
=
tanh
⁡
λ
m
cosh
⁡
λ
m
⁡
-
cos
⁡
λ
m
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
m
CF
α
1
=
1
,
α
2
=
-
1 α3=1,α4=-1
cos
⁡
λ
m
⁡
cosh
⁡
λ
m
=
-
1
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
m
cosh
⁡
λ
m
⁡
+
cos
⁡
λ
m
FS
α
1
=
1
,
α
2
=
1 α3=1,α4=1
tan
⁡
λ
m
⁡
=
tanh
⁡
λ
m
cosh
⁡
λ
m
⁡
-
cos
⁡
λ
m
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
mBy solving this eigenvalue problem, the natural frequencies of MWCNTs can be obtained and the associated eigenvectors yield the corresponding mode shapes.


## 4.1. Variational Form of the Field Equations

According to the semi-inverse method [34–36], variational trial-functional Πu1,u2,…,uN,v1,v2,…,vN,w1,w2,…,wN is constructed as(15)Π=∑i=1NΠiui,vi,w1,w2,…,wN,where(16)Πiui,vi,w1,w2,…,wN=ΠKiui,vi,wi+ΠMiui,vi,wi+ΠvdWiw1,w2,…,wNin which(17)ΠKiui,vi,wi=12∫0T∫ΩEh1-ν2u,xi+1RiνEh1-ν2v,θi+wiu,xi+Eh21+νu,θiRi+v,xi2+νEh1-ν2u,xi+1RiEh1-ν2v,θi+wi1Riv,θi+wi+Dw,xxi+νRi2w,θθiw,xxi+DRi21Ri2w,θθi+νw,xxiw,θθi+21-νDRi2w,xθi-v,xi2-2DRi2νw,xxi+1Ri2w,θθiv,θidΩdt,ΠMiui,vi,wi=12∫0T∫Ωρhu˙i2-e0a2u˙,xxi+1Ri2u˙,θθiu˙i+v˙i2-e0a2v˙,xxi+1Ri2v˙,θθiv˙i+w˙i2-e0a2w˙,xxi+1Ri2w˙,θθiw˙idΩdt,ΠvdWiw1,w2,…,wN=12∑j=1N∫0T∫Ωcijwi-wj2-e0a2wiw,xxi+1Ri2w,θθi+wjw,xxj+1Ri2w,θθjdΩdt.


## 4.2. Implementation of the Rayleigh-Ritz Method with Beam Functions

The field variables corresponding to theith tube, that is, ui, vi, and wi, can be expressed in the form(18a)uix,θ,t=Ai∂ϕx∂xcos⁡nθsin⁡ωt,(18b)vix,θ,t=Biϕxsin⁡nθsin⁡ωt,(18c)wix,θ,t=Ciϕxcos⁡nθsin⁡ωt,where Ai, Bi, and Ci are the constant parameters denoting the amplitudes of the vibrations, n is the circumferential wave number, ω is the angular frequency of vibrating CNT, and ϕx is the axial function which satisfies the geometric boundary conditions. The axial function ϕ(x) is chosen as the characteristics beam function as(19)ϕx=α1cosh⁡λmxL+α2cos⁡λmxL-ςmα3sinh⁡λmxL+α4sinλmxLin which αii=1,…,4 are constants with value 0, 1, or -1 depending on the tube ends, λm is the roots of the transcendental equations obtained from the CNT boundary conditions, and ςm is the parameters corresponding to λm. The geometric boundary conditions for clamped (C), free (F), and simply supported (S) boundary conditions can be expressed mathematically in terms of ϕx as follows.Clamped boundary condition is(20)ϕx=ϕ′x=0.Free boundary condition is(21)ϕ′′x=ϕ′′′x=0.Simply supported boundary condition is(22)ϕx=ϕ′′x=0.Parametersαi, λm, and ςm that are chosen according to the CNT boundary conditions are given in Table 1. Substituting (18a), (18b), and (18c) into (15) and then minimizing the energy functional Π with respect to unknown coefficients Ai, Bi, and Ci result in the following algebraic equations:(23)∂Π∂Ai=0,∂Π∂Bi=0,∂Π∂Ci=0,i=1,2,…,N.The above equations can be recast in the form of a generalized eigenvalue problem as(24)K-ω2MX=0,where K is the stiffness matrix, M is the mass matrix of the CNT, and X is given by(25)XT=A1,B1,C1,A2,B2,C2,…,AN,BN,CN.Table 1
Values ofαi, λm, and ζm for SS, CC, FF, CS, CF, and FS boundary conditions.
Boundary condition
α
i  i=1,…,4
Values forλm
ζ
m
SS
α
1
=
0
,
α
2
=
0 α3=0,α4=-1
m
π
1
CC
α
1
=
1
,
α
2
=
-
1 α3=1,α4=-1
cos
⁡
λ
m
⁡
cosh
⁡
λ
m
=
1
cosh
⁡
λ
m
⁡
-
cos
⁡
λ
m
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
m
FF
α
1
=
1
,
α
2
=
1 α3=1,α4=1
cos
⁡
λ
m
⁡
cosh
⁡
λ
m
=
1
cosh
⁡
λ
m
⁡
-
cos
⁡
λ
m
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
m
CS
α
1
=
1
,
α
2
=
-
1 α3=1,α4=-1
tan
⁡
λ
m
⁡
=
tanh
⁡
λ
m
cosh
⁡
λ
m
⁡
-
cos
⁡
λ
m
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
m
CF
α
1
=
1
,
α
2
=
-
1 α3=1,α4=-1
cos
⁡
λ
m
⁡
cosh
⁡
λ
m
=
-
1
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
m
cosh
⁡
λ
m
⁡
+
cos
⁡
λ
m
FS
α
1
=
1
,
α
2
=
1 α3=1,α4=1
tan
⁡
λ
m
⁡
=
tanh
⁡
λ
m
cosh
⁡
λ
m
⁡
-
cos
⁡
λ
m
sinh
⁡
λ
m
⁡
-
sin
⁡
λ
mBy solving this eigenvalue problem, the natural frequencies of MWCNTs can be obtained and the associated eigenvectors yield the corresponding mode shapes.


## 5. Numerical Results and Discussion

The mechanical properties and thickness of each tube of MWCNTs are assumed to beE=1 TPa, D=0.85 eV, h=0.34 nm, ν=0.27, and ρ=2300 Kg/m3. The configuration of layerwise boundary conditions, for example, will be denoted by (SS/CF/CC), where SS is related to the outermost tube and CC is related to the innermost tube. Also, for given intertube mode number p, for convenience, the frequency associated with the nth axial mode will be denoted by ωn(p).


### 5.1. Validation of the Present Approach

For validation, the present results are compared with the MD results given in [37]. The first resonant frequencies of clamped and cantilever SWCNTs against nanotube aspect ratio (L/d) are depicted in Figures 3 and 4, respectively. The nonlocal parameter e0a should be calibrated. The calibrated values for e0a associated with clamped-clamped and clamped-free boundary conditions are 1.98 and 2 nm, respectively. It means that the calibrated value of e0a depends on end conditions. Moreover, Figure 3 provides a comparison between the nonlocal shell model and its local counterpart for clamped end conditions. It is observed that the local shell model (e0a=0) overestimates the frequencies, particularly for small aspect ratios. As the aspect ratio increases, resonant frequencies decrease and the small scale effect diminishes so that the frequency envelopes tend to converge.Figure 3
First resonant frequencies from the present shell model and MD simulation [37] for clamped single-walled CNT (R=0.34 nm).Figure 4
First resonant frequencies from the present shell model and MD simulation [37] for cantilever single-walled CNT (R=0.34 nm).


### 5.2. Explanatory Examples

Example 1.
The fundamental resonant frequencies of a triple-walled CNT obtained by the present nonlocal shell model versus nanotube aspect ratio are graphed in Figure5. Two different layerwise boundary conditions have been considered. The values of nonlocal parameter are varied from 0 to 1 nm. From this figure, one can see that the effects of the boundary conditions and small scale become more significant for shorter nanotubes. Moreover, it is observed that the frequencies decrease by increasing e0a. As the ratio of length-to-innermost radius increases, resonant frequencies tend to decrease and the effects of small length scale and CNT end conditions diminish.Figure 5
Fundamental resonant frequency curves for a triple-walled CNT corresponding to various end conditions and nonlocal parameters (R1=8.5 nm).Example 2.
Since the effective thicknesses of nanotubes are scattered in the range of 0.066–0.34 nm, Figure6 is presented to demonstrate the influence of thickness variation on the fundamental resonant frequencies of a triple-walled CNT with clamped end conditions (CC/CC/CC). From this figure, the frequency difference can be observed due to this effect in both local and nonlocal models. Furthermore, resonant frequencies of CNTs with h=0.066 nm are higher than those of CNTs with h=0.34 nm when aspect ratio decreases.Figure 6
Effect of the nanotube thickness variation on the fundamental resonant frequencies of a fully clamped triple-walled CNT (R1=0.35 nm).Example 3.
Figure7 depicts the natural frequency of a simply supported triple-walled CNT versus circumferential mode number. Different values of nonlocal parameter have been considered ranging from 0 to 1.5 nm. Once again, it is found that the natural frequency diminishes with increasing nonlocal parameter. Physical interpretation is that the small scale effects in the nonlocal model make nanotubes more flexible. It is further observed that the magnitude of decrease in natural frequencies corresponding to higher circumferential modes is considerably higher than those corresponding to lower ones. This reveals that the influence of the small scale becomes more important for shorter wavelength at higher modes.Figure 7
Variation of natural frequency with circumferential mode number for a triple-walled CNT with simply supported end conditions (R1=8.5nm,L/R1=5).Example 4.
Figure8 shows the frequency ratio of a triple-walled (SS/SS/SS) CNT for several length-to-innermost radius ratios. It can be found that the effects of the small length scale are more significant for shorter length CNTs. For example, in the case of frequencies of nanotube with L/R1=5 a relative error equal to 26.5% for nonlocal parameter e0a=2 nm is obtained. This relative error reduces to about 9% as L/R1 increases by 20.Figure 8
Effect of the small length scale on the frequency ratio for a simply supported triple-walled CNT corresponding to various length-to-innermost radius ratios (R1=8.5 nm).Example 5.
To show the ability of the present shell model in predicting new intertube frequencies and corresponding noncoaxial vibrational modes, mode shapes of a simply supported triple-walled CNT are shown in Figure9. As revealed in this figure, noncoaxial vibrational modes are predictable. Furthermore, noncoaxial vibrational modes may shift to the ones corresponding to higher circumferential mode numbers as the radius of MWCNT increases. The three-dimensional vibrational mode shape of a simply supported triple-walled CNT associated with the first intertube and the fifth axial and circumferential modes is also plotted in Figure 10. This figure is accompanied by a cross-sectional view in the middle of the nanotube.Figure 9
Mode shapes of a simply supported triple-walled CNT predicted by the shell model (R1=2 nm, L/R1=10).Figure 10
Three-dimensional mode shape associated with the first intertube and the fifth axial and circumferential modes for a triple-walled CNT with simply supported boundary conditions (R1=8.5 nm, L/R1=5).Example 6.
This example provides a comparison between the present shell model and the beam model given by Xu et al. [30]. The first three natural frequencies of a double-walled CNT with different layerwise boundary conditions obtained by the present shell model and by the beam model [30] versus nanotube length are tabulated in Table 2. It is observed that, for long double-walled CNTs for which the beam-like vibrations are dominant, the two models agree reasonably well. However, the results obtained by the beam model for double-walled CNTs of finite length are overestimated due to not taking circumferential modes into consideration. Moreover, this reduction in the size of nanotubes brings the effects of small scale and boundary conditions into focus and accordingly makes the present nonlocal shell model more preferable than the local beam counterpart.Table 2
Comparison between the present shell model and the beam model given in [30] for the first three natural frequencies of a double-walled CNT with different boundary conditions (R1=0.35nm,e0a=0).
Frequencies (THz)
Length (nm)
14
16
18
20
24
26
28
30
CC/SS
ω
1
a
0.8760
0.6876
0.5531
0.4539
0.3208
0.2750
0.2383
0.2084
ω
1
b
1.05
0.81
0.64
0.52
0.36
0.31
0.27
0.23
ω
2
a
2.2109
1.7654
1.4385
1.1925
0.8548
0.7365
0.6407
0.5623
ω
2
b
2.84
2.20
1.75
1.43
1.00
0.85
0.73
0.64
ω
3
a
3.8679
3.1476
2.6038
2.1848
1.5944
1.3828
1.2098
1.0666
ω
3
b
5.18
4.17
3.36
2.76
1.94
1.66
1.43
1.25
CC/FF
ω
1
a
0.9290
0.7276
0.5843
0.4790
0.3381
0.2897
0.2509
0.2194
ω
1
b
1.04
0.8
0.53
0.51
0.36
0.31
0.26
0.23
ω
2
a
2.3230
1.8484
1.5027
1.2437
0.8898
0.7661
0.6662
0.5844
ω
2
b
2.84
2.2
1.75
1.42
0.99
0.85
0.73
0.64
ω
3
a
4.0609
3.2857
2.7080
2.2666
1.6491
1.4290
1.2493
1.1008
ω
3
b
5.14
4.16
3.36
2.76
1.94
1.66
1.43
1.25
CC/CC
ω
1
a
0.9277
0.7269
0.5839
0.4788
0.3380
0.2896
0.2509
0.2193
ω
1
b
1.08
—
—
—
—
—
—
—
ω
2
a
2.3176
1.8455
1.5010
1.2427
0.8893
0.7658
0.6660
0.5842
ω
2
b
2.94
—
—
—
—
—
—
—
ω
3
a
4.0420
3.2756
2.7021
2.2630
1.6475
1.4279
1.2484
1.1002
ω
3
b
5.49
—
—
—
—
—
—
—
aBased on the present shell model.
bBased on the beam model given by Xu et al. [30].Example 7.
Presented graphically in Figure11 are the frequency ratios related to a simply supported CNT with different number of walls. It is seen that the effects of the small length scale are dependent on the number of tubes so that the relative error in resonant frequencies decreases by increasing the number of walls. For e0a=2 nm as an example, the relative errors corresponding to the double-walled and five-walled CNTs are approximately 19% and 12.5%, respectively.Figure 11
Effect of the small length scale on the frequency ratio for a simply supported CNT with different number of walls (RN=8.5 nm, L/RN=10).


## 5.1. Validation of the Present Approach

For validation, the present results are compared with the MD results given in [37]. The first resonant frequencies of clamped and cantilever SWCNTs against nanotube aspect ratio (L/d) are depicted in Figures 3 and 4, respectively. The nonlocal parameter e0a should be calibrated. The calibrated values for e0a associated with clamped-clamped and clamped-free boundary conditions are 1.98 and 2 nm, respectively. It means that the calibrated value of e0a depends on end conditions. Moreover, Figure 3 provides a comparison between the nonlocal shell model and its local counterpart for clamped end conditions. It is observed that the local shell model (e0a=0) overestimates the frequencies, particularly for small aspect ratios. As the aspect ratio increases, resonant frequencies decrease and the small scale effect diminishes so that the frequency envelopes tend to converge.Figure 3
First resonant frequencies from the present shell model and MD simulation [37] for clamped single-walled CNT (R=0.34 nm).Figure 4
First resonant frequencies from the present shell model and MD simulation [37] for cantilever single-walled CNT (R=0.34 nm).


## 5.2. Explanatory Examples

Example 1.
The fundamental resonant frequencies of a triple-walled CNT obtained by the present nonlocal shell model versus nanotube aspect ratio are graphed in Figure5. Two different layerwise boundary conditions have been considered. The values of nonlocal parameter are varied from 0 to 1 nm. From this figure, one can see that the effects of the boundary conditions and small scale become more significant for shorter nanotubes. Moreover, it is observed that the frequencies decrease by increasing e0a. As the ratio of length-to-innermost radius increases, resonant frequencies tend to decrease and the effects of small length scale and CNT end conditions diminish.Figure 5
Fundamental resonant frequency curves for a triple-walled CNT corresponding to various end conditions and nonlocal parameters (R1=8.5 nm).Example 2.
Since the effective thicknesses of nanotubes are scattered in the range of 0.066–0.34 nm, Figure6 is presented to demonstrate the influence of thickness variation on the fundamental resonant frequencies of a triple-walled CNT with clamped end conditions (CC/CC/CC). From this figure, the frequency difference can be observed due to this effect in both local and nonlocal models. Furthermore, resonant frequencies of CNTs with h=0.066 nm are higher than those of CNTs with h=0.34 nm when aspect ratio decreases.Figure 6
Effect of the nanotube thickness variation on the fundamental resonant frequencies of a fully clamped triple-walled CNT (R1=0.35 nm).Example 3.
Figure7 depicts the natural frequency of a simply supported triple-walled CNT versus circumferential mode number. Different values of nonlocal parameter have been considered ranging from 0 to 1.5 nm. Once again, it is found that the natural frequency diminishes with increasing nonlocal parameter. Physical interpretation is that the small scale effects in the nonlocal model make nanotubes more flexible. It is further observed that the magnitude of decrease in natural frequencies corresponding to higher circumferential modes is considerably higher than those corresponding to lower ones. This reveals that the influence of the small scale becomes more important for shorter wavelength at higher modes.Figure 7
Variation of natural frequency with circumferential mode number for a triple-walled CNT with simply supported end conditions (R1=8.5nm,L/R1=5).Example 4.
Figure8 shows the frequency ratio of a triple-walled (SS/SS/SS) CNT for several length-to-innermost radius ratios. It can be found that the effects of the small length scale are more significant for shorter length CNTs. For example, in the case of frequencies of nanotube with L/R1=5 a relative error equal to 26.5% for nonlocal parameter e0a=2 nm is obtained. This relative error reduces to about 9% as L/R1 increases by 20.Figure 8
Effect of the small length scale on the frequency ratio for a simply supported triple-walled CNT corresponding to various length-to-innermost radius ratios (R1=8.5 nm).Example 5.
To show the ability of the present shell model in predicting new intertube frequencies and corresponding noncoaxial vibrational modes, mode shapes of a simply supported triple-walled CNT are shown in Figure9. As revealed in this figure, noncoaxial vibrational modes are predictable. Furthermore, noncoaxial vibrational modes may shift to the ones corresponding to higher circumferential mode numbers as the radius of MWCNT increases. The three-dimensional vibrational mode shape of a simply supported triple-walled CNT associated with the first intertube and the fifth axial and circumferential modes is also plotted in Figure 10. This figure is accompanied by a cross-sectional view in the middle of the nanotube.Figure 9
Mode shapes of a simply supported triple-walled CNT predicted by the shell model (R1=2 nm, L/R1=10).Figure 10
Three-dimensional mode shape associated with the first intertube and the fifth axial and circumferential modes for a triple-walled CNT with simply supported boundary conditions (R1=8.5 nm, L/R1=5).Example 6.
This example provides a comparison between the present shell model and the beam model given by Xu et al. [30]. The first three natural frequencies of a double-walled CNT with different layerwise boundary conditions obtained by the present shell model and by the beam model [30] versus nanotube length are tabulated in Table 2. It is observed that, for long double-walled CNTs for which the beam-like vibrations are dominant, the two models agree reasonably well. However, the results obtained by the beam model for double-walled CNTs of finite length are overestimated due to not taking circumferential modes into consideration. Moreover, this reduction in the size of nanotubes brings the effects of small scale and boundary conditions into focus and accordingly makes the present nonlocal shell model more preferable than the local beam counterpart.Table 2
Comparison between the present shell model and the beam model given in [30] for the first three natural frequencies of a double-walled CNT with different boundary conditions (R1=0.35nm,e0a=0).
Frequencies (THz)
Length (nm)
14
16
18
20
24
26
28
30
CC/SS
ω
1
a
0.8760
0.6876
0.5531
0.4539
0.3208
0.2750
0.2383
0.2084
ω
1
b
1.05
0.81
0.64
0.52
0.36
0.31
0.27
0.23
ω
2
a
2.2109
1.7654
1.4385
1.1925
0.8548
0.7365
0.6407
0.5623
ω
2
b
2.84
2.20
1.75
1.43
1.00
0.85
0.73
0.64
ω
3
a
3.8679
3.1476
2.6038
2.1848
1.5944
1.3828
1.2098
1.0666
ω
3
b
5.18
4.17
3.36
2.76
1.94
1.66
1.43
1.25
CC/FF
ω
1
a
0.9290
0.7276
0.5843
0.4790
0.3381
0.2897
0.2509
0.2194
ω
1
b
1.04
0.8
0.53
0.51
0.36
0.31
0.26
0.23
ω
2
a
2.3230
1.8484
1.5027
1.2437
0.8898
0.7661
0.6662
0.5844
ω
2
b
2.84
2.2
1.75
1.42
0.99
0.85
0.73
0.64
ω
3
a
4.0609
3.2857
2.7080
2.2666
1.6491
1.4290
1.2493
1.1008
ω
3
b
5.14
4.16
3.36
2.76
1.94
1.66
1.43
1.25
CC/CC
ω
1
a
0.9277
0.7269
0.5839
0.4788
0.3380
0.2896
0.2509
0.2193
ω
1
b
1.08
—
—
—
—
—
—
—
ω
2
a
2.3176
1.8455
1.5010
1.2427
0.8893
0.7658
0.6660
0.5842
ω
2
b
2.94
—
—
—
—
—
—
—
ω
3
a
4.0420
3.2756
2.7021
2.2630
1.6475
1.4279
1.2484
1.1002
ω
3
b
5.49
—
—
—
—
—
—
—
aBased on the present shell model.
bBased on the beam model given by Xu et al. [30].Example 7.
Presented graphically in Figure11 are the frequency ratios related to a simply supported CNT with different number of walls. It is seen that the effects of the small length scale are dependent on the number of tubes so that the relative error in resonant frequencies decreases by increasing the number of walls. For e0a=2 nm as an example, the relative errors corresponding to the double-walled and five-walled CNTs are approximately 19% and 12.5%, respectively.Figure 11
Effect of the small length scale on the frequency ratio for a simply supported CNT with different number of walls (RN=8.5 nm, L/RN=10).


## 6. Concluding Remarks

Considering various layerwise boundary conditions, this paper probed the free vibrations of MWCNTs based on a nonlocal elastic shell model. Using the Flügge shell theory, the displacement field equations coupled by vdW forces were derived. The variational form of the Flügge type equations was constructed to which the analytical Rayleigh-Ritz method was applied. Among the more significant conclusions to be obtained, the following findings may be summarized from the present study:(i)
The efficiency of the present shell model was checked by the MD simulation and nonlocal parameters were calibrated for clamped and clamped-free SWCNTs.(ii)
The small scale effects in the nonlocal continuum model reduce the frequencies of CNT as competed to the predictions of classical model.(iii)
It was shown that the significance of the small size effects on the natural frequencies of MWCNTs is dependent on the geometric parameters of CNTs, vibrational modes, boundary conditions, and number of walls.(iv)
The elastic beam model tends to overestimate the resonant frequencies of CNTs as compared to its shell counterpart, due to not incorporating circumferential mode number into the model, especially for shorter CNTs.(v)
The results generated provide valuable information concerning new noncoaxial modes affecting the properties of MWCNTs.



---

*Source: 750392-2015-12-27.xml*