The Scientific World Journal

Volume 2014 (2014), Article ID 194529, 13 pages

http://dx.doi.org/10.1155/2014/194529

Research Article

Free Vibration Analysis of DWCNTs Using CDM and Rayleigh-Schmidt Based on Nonlocal Euler-Bernoulli Beam Theory

Maria Anna De Rosa¹ and Maria Lippiello²

¹School of Engineering, Viale dell’Ateneo Lucano 10, 85100 Potenza, Italy

²Department of Structures for Engineering and Architecture, Via Forno Vecchio 36, 80134 Naples, Italy

Received 14 August 2013; Accepted 20 October 2013; Published 25 February 2014

Academic Editors: G. J. Gibbons, P. Lonetti, G. Nikas, and B. F. Yousif

Copyright © 2014 Maria Anna De Rosa and Maria Lippiello. 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.

Abstract.

The free vibration response of double-walled carbon nanotubes (DWCNTs) is investigated. The DWCNTs are modelled as two beams, interacting between them through the van der Waals forces, and the nonlocal Euler-Bernoulli beam theory is used. The governing equations of motion are derived using a variational approach and the free frequencies of vibrations are obtained employing two different approaches. In the first method, the two double-walled carbon nanotubes are discretized by means of the so-called “cell discretization method” (CDM) in which each nanotube is reduced to a set of rigid bars linked together by elastic cells. The resulting discrete system takes into account nonlocal effects, constraint elasticities, and the van der Waals forces. The second proposed approach, belonging to the semianalytical methods, is an optimized version of the classical Rayleigh quotient, as proposed originally by Schmidt. The resulting conditions are solved numerically. Numerical examples end the paper, in which the two approaches give lower-upper bounds to the true values, and some comparisons with existing results are offered. Comparisons of the present numerical results with those from the open literature show an excellent agreement.

1. Introduction

Carbon nanotubes (CNTs) constitute a prominent example of nanomaterials and nanostructures which have stimulated extensive research activities in science and engineering field. It is well known that CNTs are hollow cylindrical tubules composed of concentric graphitic shells with diameters on the scale of nanometers and, based on the number of walls, are designed as single-walled, double-walled and multiwalled nanotubes. In particular, the double-walled nanotubes can be viewed as two concentrically nested seamless grapheme cylinders bonded together by van der Waals forces.

Their discovery, since the publication of Iijima’s paper [1] in 1991, has attracted much attention by many researches and several studies have shown that the carbon nanotubes possess extraordinary mechanical, physical, and electrical properties. As has been pointed out in the literature, extensive studies have been conducted on studying the mechanical properties of single-walled carbon nanotubes (SWNTs) or multiwalled carbon nanotubes (MWNTs), and several investigations have been performed by employing computational and experimental methods [2, 3]. In this context, an excellent review article on the mechanical properties of nanotubes was published by Ruoff et al. [4], in which experimental measurements as well as theoretical predictions can be found.

In the earlier studies, the investigations on carbon nanotubes have mainly focused on numerous experiments [5] although these texts, at nanoscale, are very cumbersome. In addition, two different theoretical approaches have been developed: the atomistic and the continuum mechanics. Among the methods of atomistic simulations, the classical molecular dynamic (MD) simulations is the most common method in investigating the behaviour of CNTs [6, 7]. Since this approach is restricted to small-scale systems, the continuum modelling is considered to be a more appropriate method of investigating the structural behaviour of nanotubes. In the literature, there exist a lot of studies on analysing the bending, buckling, and postbuckling problems of nanotubes using Euler-Bernoulli [8, 9] and Timoshenko [10] beam models.

In recent years, due to the remarkable importance of nanostructures for many engineering and medical devices, research interest has grown on evaluating the vibrational properties of carbon nanotubes. The literature concerning the vibrational properties in CNTs is very rich and it is devoted to the dynamic problem of single and multiwalled carbon nanotubes. In this topic, the state of the art can be found in a review work by Gibson et al. [11]. In this paper, the authors report a coherent yet concise review of as many of these publications as possible and the main themes treated are the modelling and simulation of vibrating nanotubes. Extensive studies have been conducted to investigate the vibration behaviour by means of molecular dynamic simulations: for example, Ansari et al. in [12] have analyzed the vibrations characteristics of SWCNTs and DWCNTs under various layerwise boundary conditions at different lengths. The analysis performed and the results obtained show that the natural frequency of carbon nanotubes is strongly dependent on their boundary conditions especially when tubes are shorter in length. Moreover, several researchers implemented the elastic models of beams to study the dynamic problems, such as vibration and wave propagation, of carbon nanotubes [13, 14]. Xu et al. in [15, 16], for example, studied the free vibration of DWCNTs which consist of two coaxial single-walled CNT with interacting each other by the interlayer van der Waals forces. Therefore, the inner and outer tubes are modeled as two individual elastic beams and by using the Euler-Bernoulli beam model the exact solutions for natural frequencies, at different boundary conditions, have been derived. Also Elishakoff and Pentaras [17] deal with the evaluation of fundamental natural frequencies of DWCNTs under various boundary conditions and the expressions for the natural frequencies have been derived by applying the Bubnov-Galerkin and Petrov-Galerkin methods. In the paper by Natsuki et al. [18], instead, a theoretical approach to vibration characteristic analysis of CNTs with simply supported boundary condition has been presented. Applying the Euler-Bernoulli beam theory, the authors obtained the resonant frequencies of DWCNTs with different vibrational modes and they showed that the resonant frequencies decrease with increasing length of nanotubes.

Although the classical continuum methods are efficient in performing mechanical analysis of CNTs, their applicability to identify the small-scale effects on carbon nanotubes mechanical behaviours is questionable. The importance of size effect has been pointed out in a number of studies where the size dependence of the properties of nanotubes has been investigated. For example, Sun and Zhang, in [19], discussed the scarce applicability of continuous models to nanotechnology and proposed a semicontinuum model in studying nano-materials. The authors demonstrated that the values of the Young’s modulus and Poisson’s ratios depend on the number of atom layers in the thickness direction. These results show that the nanostructures and nanomaterials cannot be homogenized into a continuum. At this point, the nonlocal elastic continuum models are more pertinent in predicting the structural behaviour of nanotubes because of being capable of taking in the small-scale effects. It is well known that the nonlocal elasticity theory assumes that the stress state, at a given reference point, is considered to be a function of the strain field at all points of the body. The origins of the nonlocal theory of elasticity go to pioneering works, published in early 80s, by Eringen [20]. In [21] Reddy reports a complete development of the classical and shear deformation beam theories using the nonlocal constitutive differential equations and derived the solutions for bending, buckling, and natural frequencies problems of simply supported beams.

In recent years, many researchers have applied the nonlocal elasticity concept for the bending, buckling, and vibration analysis of nanostructures. Peddieson et al. [22] have used nonlocal Euler-Bernoulli model for static analysis of nanobeams and particular attention is paid to cantilever beams which are often used as actuators in small-scale systems. Further applications of the nonlocal elasticity theory have been employed in studying the buckling problem [23, 24] and vibration problems, by applying Euler-Bernoulli beam and shell theories and Timoshenko beam theory, in CNTs [23–32].

It is worth mentioning that, in the literature, most of the attention has been focused on deriving the variational formulation of equations and boundary conditions for single- and double-walled nanotubes undergoing vibrations with nonlocal elastic continuum methods. Reddy and Pang [33] reformulated the equation of motion of the Euler-Bernoulli and Timoshenko beam theories, using the nonlocal differential constitutive relations of Eringen. Following this approach, the equations of motion are used to evaluate the static bending, vibration, and buckling responses of beams with various boundary conditions. Adali [34] proposed a continuum model for studying the mechanical behaviour of multiwalled carbon nanotubes under compressive loads; the nonlocal theory of Euler-Bernoulli beams has been employed and the results are extended to multiwalled nanotubes subjected to transverse vibrations.

In this paper the free vibration frequencies of coaxial DWCNTs are detected, using two different approaches. The first method has already been used by the authors [35] and by Raithel and Franciosi [36] for different structural problems, and it has been properly modified for the title problem. The nanotube is reduced to a set of rigid bars, linked together by elastic cells, where masses and stiffnesses are supposed to be concentrated. The resulting discrete system is simple enough to allow to take into account nonlocal effects, constraint elasticities, and van der Walls forces, and a classical eigenvalue problem is reached, which can be easily handled by Mathematica [37].

The second proposed approach belongs to the semianalytical methods, and more precisely it can be considered an optimized version of the classical Rayleigh quotient, as proposed originally by Schmidt and then employed for various eigenvalue problems [38–40]. Basically, the trial function in the Rayleigh quotient is allowed to depend on parameters, and the resulting quotient is properly optimized. The resulting conditions can be solved numerically.

Numerical examples end the paper, in which the two approaches give lower-upper bounds to the true values, and some comparisons with existing values are offered.

2. Theoretical Approach

The beam structure, under consideration in Figure 1, is a concentric system of two nanotubes of cylindrical shape of length , with Young’s modulus and mass density . For each nanotube the cross-sectional area and the moment of inertia are defined, where the index refers to the order of the nanotubes with the inner tube indicated by and the outer tube by . For the DWCNTs, the main point in the analysis is the consideration of van der Waals (vdW) forces between the inner and outer tubes: the interaction pressure at any point between any two adjacent tubes depends on the difference of their deflections at that point. To take in to account the vdW forces, one defines the interaction coefficient between the inner and outer nanotubes, which can be estimated approximately as [14]

where is the inner radius of the wall and is the carbon-carbon bond length.

Figure 1: Geometry of double-walled carbon nanotubes (DWCNTs).

As already said, the small-scale effect is taken into account by using the nonlocal theory for Euler-Bernoulli beams, so that the parameter is introduced, where is a constant which has to be experimentally determined for each material. In turn, is an internal characteristic length, as already defined.

In order to analyze the dynamic behaviour of the system under consideration, the governing equations of motion, by considering the vdW forces and the small-scale effect, have been derived using a variational approach:

where and , with , as already defined above, are rotational and translational stiffness, respectively, at , while, analogously, and are rotational and translational stiffness at , respectively. In the above equations, the abscissa represents the spatial coordinate while is the time; in (2), (3), and (4) denotes the kinetic energy, is the strain energy and is the potential energy due to vdW forces between the two nanotubes.

3. Discretization of DWCNTs by means of CDM Method

In this section the so-called “cell discretization method” (CDM), employed to analyze the dynamic behaviour of structure under consideration, is discussed. As already said, the two nanotubes are reduced to a set of rigid bars with the same length , linked together by elastic cells (see Figure 2). The moment of inertia and the cross-sectional area with will be evaluated at the cells abscissae, obtaining the concentrated stiffness , and the concentrated masses , for the inner and outer tubes, respectively. Both these quantities can be organized into the so-called unassembled stiffness diagonal matrix and the unassembled mass diagonal matrix , with dimension (), , for each of two nanotubes.

Figure 2: Structural system CDM.

In this way, the structure is reduced to a classical holonomic system, with degrees of freedom, in particular, vertical displacements , for inner tube, and vertical displacements , for outer tube, at the cells abscissae will be conveniently assumed as Lagrangian coordinates and will be organized into the -dimensional vector . Moreover, for the inner and outer nanotubes the rotations of the rigid bars can be calculated as a function of the Lagrangian coordinates as follows:

or, in matrix form: and where is a rectangular transfer matrix with rows and columns.

The relative rotations between the two faces of the elastic cells are given by

or in matrix form: for inner rigid bar and for outer rigid bar, where is another rectangular transfer matrix with rows and columns.

The strain energies, with , (the first terms of (3)), are given by

and they are concentrated at the cells of the inner and outer tubes, respectively.

The strain energies should be expressed as functions of the Lagrangian coordinates as follows:

so that, the total strain energy can be expressed as

with and . The global assembled stiffness matrix assumes the following form:

The last term of (3), as function of the Lagrangian coordinates, assumes the following form:

or:

By assembling the terms of the (12), one gets

where and . The assembled mass matrix assumes the following form

The kinetic energy, (2), can be expressed as the sum of the following form terms:

or, in assembled form

The global assembled mass matrix is given by

The strain energy due to the vdW forces, (4), can be expressed as

where , so that (18) becomes

with equal to

The terms of matrix are given by

and the terms of the matrix assume, instead, the following form:

and are two matrices with rows and columns and have half-bandwidths equal to 2 and they can be organized into a matrix with rows and columns; so that the matrix takes the following form:

Finally, the strain energy terms of the flexible constraints at the ends are given by

with , so that the assembled stiffness matrix must be modified as follows:

The rotational stiffness of the constraints of each nanotube can be taken into account by summing up the corresponding flexibilities with the flexibilities of the rigid bars; for example, for the end constraints one gets

These terms will be organized in two matrices and furnished in (9).

Finally, the equation of motion can be written as

where is the global assembled stiffness matrix

and the global assembled mass matrix

4. Numerical Comparisons

In order to show the potentialities of the proposed approach (CDM), several numerical examples have been performed, using a general code developed in Mathematica [37], and the obtained results are compared with those of available works in the literature and listed in bibliography. In the present study, the vibration analysis was carried out for DWCNTs with the same and different boundary conditions between inner and outer nanotubes, respectively. Some numerical comparisons have been performed with reference to Adali’s paper [34] in which the fundamental frequencies of clamped-free double-walled nanotubes have been computed by Rayleigh-Ritz (R-R) method.

As first numerical example, the free frequencies of vibration of simply supported-simply supported DWCNTs have been calculated using as approximation function and with the vertical displacement of the inner and outer nanotubes given by

which satisfy the boundary conditions at the simply supported ends. According to the to the Rayleigh-Ritz method, the following equation is obtained, corresponding to (45) of the paper [34]:

where

with . Substituting (32) into (31), the two fundamental frequencies values are deduced. As first case, setting the interaction coefficient equal to zero, one obtains the first two frequencies values, , for inner and outer nanotubes, respectively. In the second one, putting  TPa (31) gives the first fundamental frequency for DWCNTs.

The simply supported-simply supported nanotubes, under consideration, have the following mechanical and geometric properties:  TPa,  g/cm³,  nm,  nm,  nm, and  nm, where and are the average radii of the inner and outer nanotubes and is the thickness of the two nanotubes. In Table 1 the numerical comparisons with the results of the proposed method (CDM) and those obtained applying Rayleigh-Ritz method are reported. As it can be seen, the fundamental frequencies show that there is an excellent agreement among the results obtained by the two different numerical procedures.

Table 1: Numerical comparison among R-R and CDM of simply supported DWCNTs: in columns 3 and 5 the first two dimensionless fundamental frequencies and , with , are reported, while in column 4, the first frequency , with  TPa, is listed.

Method

0
R-R
9.870
15.759
18.025

CDM
9.870
15.759
18.025

0.1
R-R
9.416
15.035
17.197

CDM
9.416
15.035
17.197

0.2
R-R
8.357
13.344
15.262

CDM
8.357
13.344
15.262

0.3
R-R
7.182
11.468
13.117

CDM
7.182
11.468
13.117

0.4
R-R
6.146
9.813
11.224

CDM
6.146
9.813
11.224

0.5
R-R
5.300
8.463
9.680

CDM
5.300
8.463
9.680

This system has been already solved by Reddy and Pang [33], in the absence of van der Waals interactions, using an exact approach. On the other hand, the van der Waals forces have been taken into account in [34], where an approximate Rayleigh-Ritz method has been adopted, using a single-term trial function.

In order to check the correctness of the numerical calculations of CDM, a numerical comparison with the results given by [34] is proposed applying an optimized version of the classical Rayleigh quotient, as proposed originally by Schmidt and then employed for various eigenvalue problems in [40]. In the Rayleigh-Schmidt (R-S) method, the components of displacement and rotation as suitable analytical approximation functions are assumed, in which one or more unknown parameters are presented. For example, the two-beam displacement approximation functions can be expressed as

The kinetic and strain energies and the elastic energy associated to the vdW forces can be expressed as

After same algebra, one can write

Substituting in appropriate way (33) into (34), one gets

Therefore, depends on the unknown parameters , , , which, in turn, can be obtained by minimizing (35). In fact, the properties of the Rayleigh quotient allow us to obtain the minimizing parameters by putting equal to zero the first derivatives of .

Let us consider a clamped-clamped double-walled nanotubes having the same geometric and mechanical properties of the Example 1. The analysis is carried out applying the Rayleigh-Ritz and Rayleigh-Schmidt methods and assuming the following approximation functions:

In Table 3, a numerical comparison with the results given by CDM and those obtained by R-R and R-S methods is considered. As shown, the CDM results are nearer to the Rayleigh-Schmidt values than to the Rayleigh-Ritz results.

In Table 4, the fundamental frequencies of a clamped-supported double-walled nanotube are considered, for the cases and  TPa. In the first case the nondimensional frequencies and have been calculated and in the second one the value of the first frequency has been determined.

Applying Rayleigh-Ritz and Rayleigh-Schmidt methods, the approximation functions assume the following form:

As one can see, also this numerical example confirms that the CDM results are nearer to the Rayleigh-Schmidt values than to the Rayleigh-Ritz results.

A further numerical example, of the first free frequency of vibration and for various scaling effect parameter (), is illustrated in Table 5. The different values of have been chosen so that possible comparisons with other references can be deduced. The results are presented for a single-walled nanotube with various boundary conditions at two ends which are of a variety of combinations, namely, simply supported-simply supported (SS-SS), clamped-simply supported (CL-SS), clamped-clamped (CL-CL), and clamped-free (CL-FR). The results given by CDM have been obtained neglecting the vdW forces. The numerical comparison has been done between the values of Tables 1–4, of the papers [24, 30] and the results given by [23]. As one can note, there is an excellent agreement between the obtained results for SS-SS, CL-SS, and CL-CL cases. In the CL-Fr single-walled nanotube case and for scaling effect parameter , the calculations provide values higher than those obtained for and this is impossible. The exact values are given by the proposed method (CDM), for scaling effect parameter , and reported in the last column of Table 5.

In Tables 6 and 7, the free frequencies values for clamped-sliding end (CL-SL) and sliding end-simply supported (SL-SS) double-walled nanotubes, having the same mechanical and geometric properties of the previous examples, are reported and obtained by the CDM.

In the Table 6, the clamped-sliding end (CL-SL) double-walled nanotube case is treated. The two first nondimensional frequencies have been obtained for and with . The case of sliding end-simply supported (SL-SS) for and is reported in Table 7.

All previous numerical examples show that the nondimensional frequencies decrease as the small-scale parameter increases as observed in similar studies on the free vibrations of SWCNTs [23] and DWCNTs (Tables 1, 2, 3, and 4) using the nonlocal theory.

Table 2: Numerical comparison among [33, 34] and CDM of clamped-free DWCNTs: in columns 3 and 5 the first two dimensionless fundamental frequencies and , with , are reported, while in column 4 the first frequency , with  TPa, is listed.

Method

0
[34]
3.664
5.850
6.692

[33]
3.516
—
6.422

CDM
3.516
5.614
6.421

0.1
[34]
3.568
5.697
6.517

[33]
3.531
—
6.449

CDM
3.531
5.638
6.448

0.2
[34]
3.320
5.302
6.064

[33]
3.579
—
6.537

CDM
3.570
5.702
6.520

0.3
[34]
3.002
4.793
5.483

[33]
3.669
—
6.700

CDM
3.615
5.772
6.602

Table 3: Numerical comparison among R-R, R-S, and CDM of clamped-clamped DWCNTs: in columns 3 and 5, the first two dimensionless fundamental frequencies and , with , are reported, while in column 4 the first frequency , with  TPa, is listed.

Method

0
R-R
22.793
36.394
41.628

R-S
22.410
35.783
—

CDM
22.373
35.722
40.859

0.1
R-R
21.427
34.212
39.132

R-S
21.137
33.750
—

CDM
21.109
33.704
38.551

0.2
R-R
18.449
29.458
33.694

R-S
18.303
29.225
—

CDM
18.289
29.202
33.402

0.3
R-R
15.422
24.625
28.166

R-S
15.359
24.525
—

CDM
15.353
24.515
27.001

0.4
R-R
12.934
20.652
23.622

R-S
12.907
20.609
—

CDM
12.905
20.605
21.139

0.5
R-R
11.005
17.571
20.098

R-S
10.992
17.552
—

CDM
10.991
17.550
17.273

Table 4: Numerical comparison among R-R, R-S, and CDM of clamped-simply supported DWCNTs: in columns 3 and 5, the first two dimensionless fundamental frequencies and , with , are reported, while in column 4, the first frequency , with  TPa, is listed.

Method

0
R-R
15.799
25.227
28.855

R-S
15.723
25.105
—

CDM
15.418
24.618
28.158

0.1
R-R
14.901
23.801
27.224

R-S
14.846
23.705
—

CDM
14.600
23.311
26.662

0.2
R-R
12.928
20.642
23.611

R-S
12.893
20.587
—

CDM
12.746
20.351
23.278

0.3
R-R
10.876
17.366
19.863

R-S
10.857
17.336
—

CDM
10.777
17.208
19.682

0.4
R-R
9.162
14.628
16.732

R-S
9.151
14.612
—

CDM
9.106
14.540
16.630

0.5
R-R
7.818
12.483
13.228

R-S
7.881
12.473
—

CDM
7.784
12.428
14.216

Table 5: Numerical comparison between the results obtained with CDM and [23, 24, 30] of a SWCNT and for different values of parameter .

Method
SS-SS
Cl-SS
Cl-Cl
Cl-FR

0
[23]
9.8696
15.4182
22.3733
3.5160

[30]
9.8697
15.4182
22.3733
3.5160

[24]
9.8696
15.4182
22.3733
3.5160

CDM
9.8696
15.4180
22.3728
3.5160

0.5
[23]
5.3003
7.7837
10.9914
4.0882

[30]
5.3001
7.7835
10.9912
4.0881

CDM
5.3002
7.7837
10.9913
3.5874

0.7
[23]
4.0854
5.9362
8.3483
—

[30]
4.0852
5.9362
8.3483
—

CDM
4.0854
5.9362
8.3482
—

Table 6: Numerical results for clamped-sliding end DWCNTs: in columns 1 and 3, the first two dimensionless fundamental frequencies and , with , are reported, while in column 2, the first frequency , with  TPa, is listed.

0
5.593
8.931
10.215

0.1
5.509
8.797
10.062

0.2
5.278
8.428
9.640

0.3
4.949
7.902
9.038

0.4
4.575
7.304
8.355

0.5
4.197
6.702
7.666

Table 7: Numerical results for sliding end-simply supported DWCNTs: in columns 1 and 3, the first two dimensionless fundamental frequencies and , with , are reported, while in column 2, the first frequency , with  TPa, is listed.

0
2.467
3.940
4.506

0.1
2.438
3.892
4.452

0.2
2.354
3.759
4.299

0.3
2.232
3.564
4.077

0.4
2.090
3.336
3.816

0.5
1.941
3.099
3.545

In Table 8, a numerical comparison is illustrated between the results given by CDM and Rayleigh-Ritz methods and the results given by [27, 28, 32]. The structural system under consideration is a double-walled nanotube having the following mechanical and geometric properties:  TPa,  g/cm³,  nm,  nm,  nm,  nm, and the vdW interaction coefficient is equal to  TPa.

Table 8: Numerical comparison between the results obtained with CDM and [27, 28, 32] of a DWCNT and for different values of parameter .

Boundary condition

Simply supported

     [32]
3.099
—
—

     [27]
3.14
6.27
9.35

     [28]
3.141
6.265
9.276

    C.D.M.
3.141
6.265
8.275

    R-R
3.141
—
—

     [32]
3.026
—
—

     [27]
3.07
5.78
8.01

     [28]
3.068
5.770
7.976

    C.D.M.
3.068
5.780
8.036

    R-R
3.068
—
—

Clamped

     [32]
4.482
—
—

     [27]
4.73
7.82
10.82

     [28]
4.726
7.796
10.654

    C.D.M.
4.726
7.796
10.653

    R-S
4.732
—
—

     [32]
4.359
—
—

     [27]
4.59
7.12
9.19

     [28]
4.590
7.105
9.123

    C.D.M.
4.593
7.137
9.251

    R-S
4.596
—
—

Propped

     [32]
3.802
—
—

     [27]
3.93
7.05
10.09

     [28]
3.925
7.035
9.981

    C.D.M.
3.925
7.035
9.981

    R-S
3.965
—
—

     [32]
3.701
—
—

     [27]
3.82
6.45
8.60

     [28]
3.819
6.444
8.557

    C.D.M.
3.820
6.463
8.643

    R-S
3.853
—
—

Cantilever

     [32]
1.88
4.69
7.82

     [27]
1.875
4.690
7.797

    C.D.M.
1.875
4.690
7.797

     [27]
1.88
4.55
7.13

     [28]
1.879
4.544
7.111

    C.D.M.
1.879
4.547
7.143

The numerical calculations have been performed for simply supported-simply supported, clamped-clamped, clamped-simply supported, and clamped-free double-walled nanotubes and setting the small-scale parameter equal to .

With reference to the paper [27], the free frequencies of double-walled carbon nanotubes (DWCNTs), , with and , are calculated. The structural system is modeled following the nonlocal Euler beams theory and using the Galerkin approach. In [32], the effects of small-scale parameters on the vibrations of DWCNTs, embedded in elastic medium and based on nonlocal Timoshenko theory, are examined in detail. Finally, in [28] the fundamental free frequencies of DWCNTs have been found by means of two analytical approaches in which solving the coupled governing equations of the motion are solved.

The numerical comparisons listed in Table 8 show that there is an excellent agreement between the CDM values and those obtained in [28].

In the following numerical examples the cases of double-walled carbon nanotubes with different boundary conditions between the inner and outer tubes is considered and the obtained results are compared with the values reported in [15]. Moreover, a numerical comparison is given for the case of clamped-clamped nanotubes between the CDM and the Rayleigh-Schmidt values. The nanotubes under consideration have the following mechanical and geometric properties:  TPa,  g/cm³,  nm,  nm,  nm, and  nm, where and are the average radii of the inner and outer nanotubes and is the thickness of the two nanotubes. The vdW interaction coefficient is equal to  GPa and the nonlocal effects are neglected. In Table 9, the first six free frequencies for different boundary conditions are reported. As one can see, the obtained results with CDM method are lower than those by [15].

Table 9: First free frequencies of vibrations of DWCNTs with different boundary conditions between inner and outer tubes, (for ).

Boundary condition
Method

Fr-Fr inner
[15]
1.040
2.84
5.140
7.890
8.130
8.380
9.350

Cl-Cl outer
CDM
1.021
2.786
5.760
6.270
6.669
7.915
8.083

Fr-Fr inner
[15]
0.170
1.040
2.890
5.290
6.550
7.890
8.170

Cl-Fr outer
CDM
0.163
1.024
2.833
5.210
6.491
7.860
8.041

SS-SS inner
[15]
1.050
2.840
5.180
7.290
7.890
8.240
9.080

Cl-Cl outer
CDM
1.0219
2.789
5.104
7.212
7.914
8.243
9.010

Cl-Cl inner
[15]
1.080
2.940
5.490
7.900
8.130
8.240
—

Cl-Cl outer
CDM
1.058
2.881
5.394
7.925
8.00
8.255
9.400

R-S
1.061
—
—
—
—
—
—

Recently, Hemmatnezhad and Ansari [31] have furnished a finite element formulation for the free vibration analysis of embedded double-walled carbon nanotubes based on nonlocal Timoshenko beam theory. In their analysis a numerical comparison is offered with the results given by [15] and reported in Table 4 in [31] for the case of DWCNTs so constrained: the simply supported-simply supported inner tube and clamped-clamped outer tube. In Table 10 the previous numerical comparison is proposed introducing the results given by CDM and differential quadrature method (DQM), developed in [36, 41].

Table 10: Numerical comparison among [15, 31] CDM and DQM and the first three free frequencies of vibrations of DWCNTs are reported (for ).

Mode
[31]
[15]
CDM
DQM

14

0.9097
1.05
1.0298
1.0299

2.3601
2.84
2.7895
2.7896

4.1481
5.18
5.1043
5.1045

18

0.5751
0.64
0.6276
0.6276

1.4648
1.75
1.7202
1.7203

2.6153
3.36
3.3036
3.3039

24

0.3340
0.36
0.3550
0.3549

0.8390
1.00
0.9770
0.9767

1.5162
1.94
1.9046
1.9046

28

0.2482
0.27
0.2614
0.2612

0.6203
0.73
0.7200
7.1905

1.1261
1.43
1.4074
1.4066

As can be observed the first three free frequencies values given by [15] overestimate the frequencies and the obtained results employing the CDM and DQM methods are lower than those by [15] while they greater than those by [31].

In Table 11, the nonlocal effects influence on the first three free frequencies of vibration are investigated, putting the small-scale parameter equal to .

Table 11: Numerical comparison among [31] and CDM and the first three frequencies of vibrations of DWCNTs are reported (for and ).

Mode

[31]
CDM
[31]
CDM

14

0.8678
0.997
0.7103
0.866

1.7807
2.381
1.1764
1.726

2.7046
4.046
1.7232
2.583

18

0.5471
0.604
0.4328
0.524

1.0966
1.459
0.7206
1.045

1.6772
2.453
1.0936
1.563

24

0.3171
0.340
0.2446
0.295

0.6247
0.822
0.4089
0.588

0.9606
1.382
0.6379
0.879

28

0.2355
0.250
0.1800
0.2168

0.4610
0.604
0.3013
0.432

0.7103
1.016
0.4749
0.647

The numerical calculations relative to a DWCNTs, having various lengths (14 nm, 18 nm, 24 nm, and 28 nm), are performed and the results compared with the CDM method and FEM approach given by [31]. As shown, the results given by CDM method are greater than those by [31].

In the last Tables 12 and 13 the nondimensional free frequencies values are reported for different boundary conditions among inner and outer nanotubes, respectively, and for different values of the small-scale parameter . The numerical calculations have been performed using the geometrical and physical data given by [34]. In particular the first nondimensional frequency value, , is reported. As shown, the fundamental frequencies decrease with increasing values of the small-scale parameter .

Table 12: Fundamental frequency of vibration of DWCNTs with different boundary conditions between inner and outer tubes and for .

SS-SS inner

Cl-SS outer
SS-SS inner

Cl-Cl outer
SS-SS inner

Cl-Fr outer
SS-SS inner

Cl-Sl outer

0
24.567
35.558
24.564
35.408

0.1
23.304
33.683
23.303
33.684

0.2
20.346
29.187
20.346
29.188

0.3
17.204
24.504
17.204
24.504

0.4
14.537
20.597
14.537
20.597

0.5
12.426
17.543
12.426
17.543

Table 13: Fundamental frequency of vibration of DWCNTs with different boundary conditions between inner and outer tubes and for .

Cl-Fr inner

Cl-SS outer
Cl-Fr inner

Cl-Cl outer
Cl-Fr inner

SS-Sl outer
Cl-Fr inner

Cl-Sl outer

0
24.618
35.577
8.231
8.916

0.1
23.311
33.694
8.676
8.796

0.2
20.351
29.195
8.322
8.427

0.3
17.208
24.510
7.809
7.901

0.4
14.540
20.601
7.225
7.304

0.5
12.428
17.547
6.634
6.702

5. Conclusions

Coaxial DWCNTs are modeled as two beams interacting between then through van der Waals forces, and nonlocal Euler-Bernoulli beam theory is employed in order to calculate the free vibration frequencies of the system. Two different numerical approaches are used in order to perform numerical comparisons. In the first method, the system, under consideration, has been modeled as a set of rigid bars linked together by elastic cells, where masses and stiffnesses are supposed to be concentrated. The resulting finite degree of freedom has allowed taking into account nonlocal effects, constraint elasticities, and van der Walls forces. The second proposed approach belongs to the semianalytical methods, and more precisely it can be considered an optimized version of the classical Rayleigh quotient, as proposed originally by Schmidt and then employed for various eigenvalue problems.

Several numerical examples have been treated in detail, comparing numerical and approximate results from the literature, and the proposed approaches have furnished excellent results.

More particularly, emphasis has been given to the influence of the small-scale parameter, of the length of the nanotubes and of the various boundary conditions on the free vibration frequency behaviour.

In the author’s opinion, the first method will be particularly useful—because of its intrinsic simplicity—in the future analysis of mass nanosensor and of nanotube in the presence of soil and follower forces.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

S. Iijima, “Helical microtubules of graphitic carbon,” Nature, vol. 354, no. 6348, pp. 56–58, 1991.

A. Krishnan, E. Dujardin, T. W. Ebbesen, P. N. Yianilos, and M. M. J. Treacy, “Young's modulus of single-walled nanotubes,” Physical Review B, vol. 58, no. 20, pp. 14013–14019, 1998.

B. G. Demczyk, Y. M. Wang, J. Cumings et al., “Direct mechanical measurement of the tensile strength and elastic modulus of multiwalled carbon nanotubes,” Materials Science and Engineering A, vol. 334, no. 1-2, pp. 173–178, 2002.

R. S. Ruoff, D. Qian, and W. K. Liu, “Mechanical properties of carbon nanotubes: theoretical predictions and experimental measurements,” Comptes Rendus Physique, vol. 4, no. 9, pp. 993–1008, 2003.

M. M. J. Treacy, T. W. Ebbesen, and J. M. Gibson, “Exceptionally high Young's modulus observed for individual carbon nanotubes,” Nature, vol. 381, no. 6584, pp. 678–680, 1996.

J. Tersoff and R. S. Ruoff, “Structural properties of a carbon-nanotube crystal,” Physical Review Letters, vol. 73, no. 5, pp. 676–679, 1994.

M. B. Nardelli, B. I. Yakobson, and J. Bernholc, “Brittle and ductile behavior in carbon nanotubes,” Physical Review Letters, vol. 81, no. 21, pp. 4656–4659, 1998.

C. Q. Ru, “Column buckling of multiwalled carbon nanotubes with interlayer radial displacements,” Physical Review B, vol. 62, no. 24, pp. 16962–16967, 2000.

Q. Wang, T. Hu, and Q. J. Chen, “Bending instability characteristics of double walled nanotubes,” Physical Review B, vol. 71, no. 4, pp. 045403–045412, 2005.

Y. Y. Zhang, C. M. Wang, and V. B. C. Tan, “Buckling of multiwalled carbon nanotubes using Timoshenko beam theory,” Journal of Engineering Mechanics, vol. 132, no. 9, pp. 952–958, 2006.

R. F. Gibson, E. O. Ayorinde, and Y.-F. Wen, “Vibrations of carbon nanotubes and their composites: a review,” Composites Science and Technology, vol. 67, no. 1, pp. 1–28, 2007.

R. Ansari, S. Ajori, and B. Arash, “Vibrations of single- and double-walled carbon nanotubes with layerwise boundary conditions: a molecular dynamics study,” Current Applied Physics, vol. 12, no. 3, pp. 707–711, 2012.

J. Yoon, C. J. Ru, and A. Mioduchowski, “Non-coaxial resonance of an isolated multiwall carbon nanotubes,” Physical Review B, vol. 66, pp. 233402–233409, 2002.

J. Yoon, C. Q. Ru, and A. Mioduchowski, “Vibration of an embedded multiwall carbon nanotube,” Composites Science and Technology, vol. 63, no. 11, pp. 1533–1542, 2003.

K.-Y. Xu, E. C. Aifantis, and Y.-H. Yan, “Vibrations of double-walled carbon nanotubes with different boundary conditions between inner and outer tubes,” Journal of Applied Mechanics, vol. 75, no. 2, pp. 0210131–0210139, 2008.

K. Y. Xu, X. N. Guo, and C. Q. Ru, “Vibration of a double-walled carbon nanotube aroused by nonlinear intertube van der Waals forces,” Journal of Applied Physics, vol. 99, no. 6, pp. 064303–064307, 2006.

I. Elishakoff and D. Pentaras, “Fundamental natural frequencies of double-walled carbon nanotubes,” Journal of Sound and Vibration, vol. 322, no. 4-5, pp. 652–664, 2009.

T. Natsuki, Q.-Q. Ni, and M. Endo, “Analysis of the vibration characteristics of double-walled carbon nanotubes,” Carbon, vol. 46, no. 12, pp. 1570–1573, 2008.

C. T. Sun and H. Zhang, “Size-dependent elastic moduli of platelike nanomaterials,” Journal of Applied Physics, vol. 93, no. 2, pp. 1212–1218, 2003.

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.

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.

J. Peddieson, G. R. Buchanan, and R. P. McNitt, “Application of nonlocal continuum models to nanotechnology,” International Journal of Engineering Science, vol. 41, no. 3-5, pp. 305–312, 2003.

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.

S. C. Pradha and J. K. Phadikar, “Bending, buckling and vibration analyses of nonhomogeneous nanotubes using GDQ and nonlocal elasticity theory,” Structural Engineering and Mechanics, vol. 33, no. 2, pp. 193–213, 2009.

Q. Wang and V. K. Varadan, “Vibration of carbon nanotubes studied using nonlocal continuum mechanics,” Smart Materials and Structures, vol. 15, no. 2, pp. 659–666, 2006.

R. Ansari and S. Sahmani, “Small scale effect on vibrational response of single-walled carbon nanotubes with different boundary conditions based on nonlocal beam models,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 4, pp. 1965–1979, 2012.

A. Shakouri, R. M. Lin, and T. Y. Ng, “Free flexural vibration studies of double-walled carbon nanotubes with different boundary conditions and modeled as nonlocal Euler beams via the Galerkin method,” Journal of Applied Physics, vol. 106, no. 9, Article ID 094307, pp. 094307–094316, 2009.

H. Ehteshami and M. A. Hajabasi, “Analytical approaches for vibration analysis of multi-walled carbon nanotubes modeled as multiple nonlocal Euler beams,” Physica E, vol. 44, no. 1, pp. 270–285, 2011.

R. Ansari and H. Rouhi, “Analytical treatment of the free vibration of single-walled carbon nanotubes based on the nonlocal Flugge shell theory,” Journal of Engineering Materials and Technology, vol. 134, no. 1, Article ID 011008, pp. 011008–011016, 2012.

C. M. Wang, Y. Y. Zhang, and X. Q. He, “Vibration of nonlocal Timoshenko beams,” Nanotechnology, vol. 18, no. 10, Article ID 105401, pp. 105401–105409, 2007.

M. Hemmatnezhad and R. Ansari, “Finite element formulation for the free vibration analysis of embedded double-walled carbon nanotubes based on nonlocal Timoshenko beam theory,” Journal of Theoretical and Applied Physics, vol. 44, no. 5, pp. 270–285, 2013.

L. L. Ke, Y. Xiang, J. Yang, and S. Kitipornchai, “Nonlinear free vibration of embedded double-walled carbon nanotubes based on nonlocal Timoshenko beam theory,” Computational Materials Science, vol. 47, no. 2, pp. 409–417, 2009.

J. N. Reddy and S. D. Pang, “Nonlocal continuum theories of beams for the analysis of carbon nanotubes,” Journal of Applied Physics, vol. 103, no. 2, Article ID 023511, pp. 023511–023526, 2008.

S. Adali, “Variational principles for transversely vibrating multiwalled carbon nanotubes based on nonlocal euler-bernoulli beam model,” Nano Letters, vol. 9, no. 5, pp. 1737–1741, 2009.

M. A. De Rosa and M. Lippiello, “Natural vibration frequencies of tapered beams,” Engineering Transactions, vol. 57, no. 1, pp. 44–66, 2009.

A. Raithel and C. Franciosi, “Dynamic analysis of arches using Lagrangian approach,” Journal of Structural Engineering, vol. 110, no. 4, pp. 847–858, 1984.

S. Wolfram, The Mathematica 8, Wolfram Media/Cambridge University Press, 2010.

C. W. Bert, “Application of a version of the Rayleigh technique to problems of bars, beams, columns, membranes, and plates,” Journal of Sound and Vibration, vol. 119, no. 2, pp. 317–326, 1987.

P. A. A. Laura, “Optimization of variational methods,” Ocean Engineering, vol. 22, no. 3, pp. 235–250, 1995.

M. A. De Rosa and C. Franciosi, “The optimized Rayleigh method and mathematica in vibrations and buckling problems,” Journal of Sound and Vibration, vol. 191, no. 5, pp. 795–808, 1996.

M. A. De Rosa and M. Lippiello, “Non-classical boundary conditions and DQM for double-beams,” Mechanics Research Communications, vol. 34, no. 7-8, pp. 538–544, 2007.