## 1. Introduction

Carbon nanotubes are tubules of carbon, in the dimensional range of *nanometre* (in the order of 10e−9 m), with atoms arranged in a way that gives them exceptional properties. The high aspect ratio, the ratio of length to its diameter, and chirality are the factors that contribute to the same. Carbon nanotubes are classified as single- and multi-walled depending on the number of tubes held concentrically. A considerable amount of research has been dedicated to the study of mechanical properties of both these types using experimental or theoretical methods. Theoretical methods involve the representation of C**─**C bonds as is, conveniently named as atomistic modelling methods, or as continuous structures such as bar, beam and shell. The latter is usually employed to study the vibrational properties as the use of atomistic models require comparatively large computational time and effort, and are often accurate for a particular time period in the loading cycle. The laws of continuum structures have been found to be at the nanometre scale as well [1].

Continuum structures, by definition, consider the model to be continuous and ignore the inter-atomic forces. Eringen and Edelen [2] and Eringen [3, 4] proposed a model that would connect the lattice mechanics with the continuum mechanics, named as nonlocal continuum mechanics. This is achieved by introducing nonlocal factors in the governing equations [5]. As will be seen in the following sections, the introduction of these nonlocal factors have shown improved description of the model. Peddieson et al. [6] extended the nonlocal continuum models to nanotechnology and hence paved the way for these theoretical models to represent the properties of nanomaterials with greater accuracy.

Finite element methods have been employed for analysis of continuum structures since 1950s, as the numerical technique is very adaptable in terms of complex geometries and boundary conditions. The method involves the division of the structure into number of ‘elements’ and using appropriate interpolations functions to express the approximate solutions. The eigenvalue problem hence achieved is solved for system’s natural frequencies and mode shapes.

The study of vibrations of carbon nanotubes, modeled as continuum models has gained considerable attention with beams and shells being used extensively. The elastic model has been presented by Ru [7] to study the column buckling of multi-walled carbon nanotubes, where the adjacent tubes were presented to be interacting through van der Waals interaction forces. Ru [8] also studied the transverse vibrations of single- and multi-walled carbon nanotubes by considering them as Euler-Bernoulli beams. Sufficient attention was given to the effect of van der Waals interaction forces on the vibrational frequency. Yoon et al. [9] studied the non-coaxial vibrational modes of multi-walled carbon nanotubes embedded in an elastic medium and deduced that the concentric geometry of carbon nanotubes will be distorted at ultrahigh frequencies. Nonlinear free vibrations of multi-walled carbon nanotubes embedded in an elastic medium were studied by Fu et al. [10]. The relation between nonlinear free vibration frequency, stiffness of the surrounding medium and van der Waals interaction forces was discussed and amplitude curves presented. Vibrational characteristic analysis of double-walled carbon nanotube (DWCNT) with simply supported boundary condition using multiple Euler-Bernoulli beam model was investigated by Natsuki et al. [11]. It was seen that the amplitudes of inner and outer tubes are non-coaxial and in opposite direction. Natsuki also studied DWCNTs embedded in elastic medium, represented by Winkler spring model, for free supported boundary condition [12] and DWCNTs with different lengths for inner and outer tubes [13]. Xu et al. [14] studied the vibrations of DWCNTs with nonlinear van der Waals interaction forces and concluded that it does not affect the coaxial free vibrations. Elishakoff and Pentaras [15, 16] studied the natural frequencies of DWCNTs and TWCNTs using approximate methods, like Bubnov-Galerkin and Petrov-Galerkin and exact solutions. The coupling between the concentric tubes is considered as van der Waals interaction forces, in both the cases.

Nonlocal continuum mechanics approach has also been used extensively in the recent years. Constitutive relations of nonlocal elasticity to be used in analysis of carbon nanotubes for Euler-Bernoulli, Timoshenko beams and cylindrical shells were presented by Wang and Wang [17]. Extensive attention was given to the shear stress and strain relation for the Timoshenko beams. Wang [18] also studied the wave propagation in carbon nanotubes using nonlocal Euler-Bernoulli and nonlocal Timoshenko beams and concentrated on the effect of nonlocal factor on the mechanical behaviour of carbon nanotubes. Yoon et al. [19] also performed a similar study concluded that the rotary inertia and shear deformation play a significant role when the measured natural frequency is much closer to the critical frequency. One more such study was performed by Wang and Varadan [20], who developed explicit frequency solutions for both single-walled carbon nanotubes (SWCNTs) and double-walled carbon nanotube (DWCNTs) on the basis of nonlocal continuum mechanics and employed both the beam theories. The ratio of local and nonlocal frequencies is used to study the effect of nonlocal factor. Arash and Wang [21] reviewed various studies on the nonlocal continuum theory for carbon nanotubes and graphene sheets. They also compared the nonlocal models to the local equivalents. Reddy [22] reformulated governing equations, based on nonlocal theory, and presented analytical solutions for bending, buckling and vibrations of Euler-Bernoulli, Timoshenko, Reddy and Levinson beams. Further, Reddy and Pang [23] also presented the governing equations of Euler-Bernoulli and Timoshenko beams, based on nonlocal theory, specifically for carbon nanotubes and used these equations to study the vibrations, buckling and static bending. Murmu and Pradhan [24] employed the differential quadrature (DQ) method to study the thermo-mechanical vibration of SWCNT embedded in an elastic medium, represented as Winkler-type elastic foundation, using nonlocal elasticity theory. It was found that the difference between local and nonlocal frequencies is particularly high for low temperature changes. Narendar et al. [25] studied the wave propagation in SWCNTs, under magnetic field, using nonlocal Euler-Bernoulli beam theory. De Rosa and Lippiello [26] analysed the free vibration of DWCNTs, modeled as nonlocal Euler-Bernoulli beams, using cell-discretization method (CDM) and a semi-analytical Rayleigh-Schmidt method. A rare study of forced vibrations of carbon nanotubes was performed by Karaoglu and Aydogyu [27] using both local and nonlocal Euler-Bernoulli beam models, for both SWCNTs and DWCNTs. Khosrozadeh and Hajabasi [28] studied the free vibrations of DWCNTs embedded in an elastic medium with nonlinear interlayer van der Waals interaction forces using nonlocal Euler-Bernoulli beam model. Fang et al. [29] performed a similar study for DWCNTs embedded in elastic medium considering it to be Winkler model.

Timoshenko models have also been used extensively as it considers the rotational inertia and transverse shear. Wang et al. [30] considered it and studied the free vibration of MWCNTs using differential quadrature (DQ) method to solve the governing equations. It is seen that the frequencies are over-predicted by Euler-Bernoulli model when the aspect ratio is small. Yang et al. [31] studied the nonlinear free vibration of SWCNT with the von Karman geometric nonlinearity considering them as nonlocal Timoshenko beams. Nonlinear free vibration of embedded DWCNTs modeled as nonlocal Timoshenko beams was presented by Ke et al. [32] where again differential quadrature (DQ) technique was employed, and the surrounding elastic medium was represented by spring-based Winkler model. Pradhan [33] presented a finite element model using nonlocal Timoshenko beams and Galerkin weighted residual technique. The buckling of carbon nanotubes using Timoshenko model was studied by Zhang et al. [34], and various studies on buckling were extensively reviewed by Wang et al. [35] and Shima [36].

A theory similar to nonlocal theory was presented by Park and Gao [37] using modified couple stress theory where internal material length parameter is considered and a variational model is presented basing on minimum total energy potential for Euler-Bernoulli beam model. It found a relation between deflections and beam thickness. Another nonlocal theory, for nanobeams, was presented by Thai [38] in which the governing equations were derived using Hamilton’s principle and analytical solutions were presented. The results were found to be comparable to Timoshenko and Reddy beam theories even though the formulation is similar to that of Euler-Bernoulli theory. A generalized nonlocal beam theory, to study properties of nanobeams, was presented by Aydogdu [39]. Various beam models were studied as a special case of this beam theory for simply supported nanobeams, and the numerical results were found to be stand valid. Lu et al. [40] compared various nonlocal beam models that are used for carbon nanotubes and made some interesting observations.

Various variational techniques have also been employed with nonlocal theories to study the vibrational behaviour of carbon nanotubes. Adali [41] employed a semi-inverse method and presented one such model. Alshorbagy et al. [42] developed a nonlocal finite element for static analysis of nanobeams. Eltaher et al. [43] analysed the vibrations of nonlocal Euler-Bernoulli nanobeams using finite element method. A similar study was performed by Phadikar and Pradhan [44], where finite element analysis of nanoplates and nanobeams was presented. All these models are found to be efficient for SWCNTs. Ansari and Hemmatnezhad [45] studied the nonlinear vibrations of MWCNTs embedded in an elastic medium using variational iteration method. A finite element analysis of embedded Timoshenko beam was also performed by Hemmatnezhad and Ansari [46], where carbon nanotubes were modeled as nonlocal beam using statically exact shape functions to achieve an eigenvalue problem. Shakouri et al. [47] studied the transverse vibrations of DWCNTs modeled as nonlocal Euler-Bernoulli beams using Galerkin method. Dinçkal [48], more recently, proposed a finite element model that uses exact dynamic shape functions (trigonometric in nature) to form the frequency-dependant dynamic stiffness matrix and hence study the free vibrations of carbon nanotubes, modeled as both nonlocal Euler-Bernoulli and Timoshenko beams. Ehteshami and Hajabasi [49] and Zhang et al. [50] performed an analytical study and presented exact solutions for finding natural frequencies of DWCNTs modeled as nonlocal Euler-Bernoulli beams. The results of both these studies have been considered as benchmark to verify the results obtained in the current study. This chapter exploits the generality of the FEM technique to obtain the natural frequencies of DWCNTs, coupled with van der Waals interaction forces. The model has been proposed and applied to both local (classical) Euler-Bernoulli and nonlocal Euler-Bernoulli beam cases. The results are compared with the values presented in the literature.

## 2. Formulation

This section gives the formulation for the proposed models. Section 2.1 would present the formulation for local Euler-Bernoulli beam and Section 2.2 for nonlocal Euler-Bernoulli beam. As will be seen ahead, the formulation uses Galerkin weighted residual method and Hermite cubic polynomial functions, in both the cases in interest of generality.

### 2.1. DWCNTs based on local Euler-Bernoulli beam

Euler-Bernoulli beam, also known as the classical beam, is known to efficiently represent a carbon nanotube, in the case when the aspect ratio (

where

where

and

The governing equations for a DWCNT are [8]:

where the subscripts 1, 2 for

Following Galerkin weighted residual method, weight functions are introduced as

(7) |

The weak integral form of the governing equations is eliminated by performing integration by parts, twice on Eq. (7), as given below:

The boundary terms in Eq. (8) which vanish with the application of system boundary conditions, represent the shear force and moment of the beams are given by

(9) |

The present work uses Hermite type polynomial approximations to derive the finite element formulation such that the interpolation function is given as:

where

### 2.2. DWCNTs based on nonlocal Euler-Bernoulli beam

The concept of nonlocality, which is accounted for with the introduction of a nonlocal or a small scale factor was extended to the field of nanotechnology by Peddieson et al. [6] where the author applied it to Euler-Bernoulli beam model. The small scale effect or the nonlocal parameter is given as

with *x* axis, are given by:

where

For nonlocal DWCNTs, the governing equations are derived to be [20, 49, 50]:

(15) |

The solution of the equations is in the form

(16) |

To maintain uniformity in mathematical modelling with the formulation proposed for local Euler-Bernoulli model, Galerkin weighted residual model is used for FEM, with

(17) |

Integration by parts is performed on the above set of equations once, and the weak form vanishes with the application of system boundary conditions,

(18) |

The resulting internal shear forces for both the beams are given as:

and the resultant bending moments are given as:

Integration by parts is performed once again on the first terms in Eq. (18) and after rearranging for clarity, the following set of equations is obtained:

(21) |

Using Hermite cubic polynomial functions, given in Eq. (11), the finite element model is developed. The element mass and stiffness matrices for both the local and nonlocal scenarios are explicitly presented in the Appendix A. An FEM code in MATLAB solves the linear eigenvalue problem which is given as:

where **Table 1** gives the element matrices

## 3. Results and discussion

The robustness and generality of proposed FEM model is demonstrated by its application to a double-elastic Euler-Bernoulli beam and then to carbon nanotubes modeled as Euler-Bernoulli beam, both locally and nonlocally. Natural frequencies are calculated for clamped-free, clamped-clamped and simply supported-simply supported boundary conditions. The parameters of the local double elastic Euler-Bernoulli beam are taken as:

These parameters are so selected that the double beam replicates the geometry of carbon nanotubes, i.e., both are concentric and are coupled with an interaction forces between the walls. The natural frequencies of the double-elastic system, in this chapter, are normalized using **Table 2**.

The percentage error has been calculated by **Table 3**. It is to be noted that the values are presented for a single value nonlocal factor, where **Table 1** as that would reduce the nonlocal formulation to local one. The values presented in **Table 1** are compared to the values obtained using the beam parameters in analytical solutions provided by Zhang et al. [50]. It can be seen the error percentages are small hence validating the proposed model.

**Table 4** gives the error percentage between local and nonlocal double-elastic Euler-Bernoulli beam when

It is concluded that a novel double-elastic Euler-Bernoulli beam FEM model has been proposed, and its veracity established for both local and nonlocal cases, the same model will now be applied to carbon nanotubes with the parameters used in Xu et al. [14] and Elishakoff and Pentaras [15], given below:

The natural frequencies of a local DWCNT, modeled as Euler-Bernoulli beams, are presented in **Table 5**. The nonlocal factor **Table 5** confirm that the proposed formulation holds valid even at nanoscale level. Further, negative error percentage denotes that the values obtained using the present FEM model are a little higher than those presented by Ehteshami and Hajabasi [49] (see **Table 6**).

Boundary conditions | [49] | [50] | Present study | Error percentage 1 | Error percentage 2 |
---|---|---|---|---|---|

C─F | 0.166 | 0.16 | 0.16 | 0 | 0.01 |

C─C | 1.06 | 1.06 | 1.06 | −0.04 | 0.04 |

S─S | 0.467 | 0.46 | 0.46 | −0.03 | 0 |

1 | Local | Nonlocal | Error percentage | |
---|---|---|---|---|

C─F | 20 | 0.41 | 0.41 | 0 |

C─C | 2.65 | 2.65 | 0.01 | |

S─S | 1.17 | 1.17 | 0 | |

C─F | 40 | 0.10 | 0.10 | 0 |

C─C | 0.66 | 0.66 | 0 | |

S─S | 0.29 | 0.29 | 0 | |

C─F | 60 | 0.04 | 0.04 | 0 |

C─C | 0.29 | 0.29 | 0 | |

S─S | 0.13 | 0.13 | 0 | |

C─F | 80 | 0.02 | 0.02 | 0 |

C─C | 0.16 | 0.16 | 0 | |

S─S | 0.07 | 0.07 | 0 |

Employing the FEM model proposed in Section 1.2, the natural frequencies of nonlocal DWCNTs modeled as Euler-Bernoulli beam are given in **Table 6**. The natural frequencies thus obtained are compared to the values obtained using the analytical formulations given by Ehteshami and Hajabasi [49] Zhang et al. [50]. The error percentages, between present study and these two studies, are presented as error percentage 1 and 2, respectively. The effect of nonlocal factor on the natural frequencies is best understood when the error percentage between local results and nonlocal results are compared, presented in **Table 7**. It will be noted that the error percentage is very small compared to the error percentage reported in **Table 4**. This difference is attributed to the value of nonlocal factor—the double-elastic beam model considers it to be between 1 and 2 and DWCNTs have it in the order of **Table 8**) that as

It is noted that though many tables in this chapter may have identical values in the columns that are being compared, the indicated percentage error is a nonzero number, or on the contrary, the columns being compared might not have identical values, but the indicated error percentage is zero. A justification for these values stems from the large magnitude of the numbers reported. For the sake of uniformity, we have rounded off the numbers to 2 decimal points.

## 4. Conclusion

This chapter proposed FEM formulations to find the natural frequencies of DWCNTs modeled as, both, local and nonlocal Euler-Bernoulli beams, coupled with van der Waals interaction forces. It uses Galerkin weighted residual approach and Hermite cubic polynomial functions to form the stiffness and mass matrices. The formulation is applied to double-elastic Euler-Bernoulli beam and to DWCNTs, to check its generality. The natural frequencies thus obtained are in agreement with those of the formulations found in the literature. It is seen that as the