Next Article in Journal
Conducting Polymer Metallic Emerald: Magnetic Measurements of Nanocarbons/Polyaniline and Preparation of Plastic Composites
Next Article in Special Issue
Physical Processes Occurring in Dispersed Media with Carbon Nanomaterials under the Influence of Ultrasonification
Previous Article in Journal
Lauric Acid Treatments to Oxidized and Control Biochars and Their Effects on Rubber Composite Tensile Properties
Previous Article in Special Issue
Stability of Carboxyl-Functionalized Carbon Nanotubes in Simulated Cement Pore Solution and Its Effect on the Compressive Strength and Porosity of Cement-Based Nanocomposites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applicability and Limitations of Ru’s Formulation for Vibration Modelling of Double-Walled Carbon Nanotubes

Department of Sciences and Methods for Engineering, University of Modena and Reggio Emilia, Via Giovanni Amendola 2, 42122 Reggio Emilia, Italy
Submission received: 13 September 2022 / Revised: 27 October 2022 / Accepted: 31 October 2022 / Published: 2 November 2022
(This article belongs to the Collection Novel Applications of Carbon Nanotube-Based Materials)

Abstract

:
In this paper, a comparison is conducted between two different formulations of the van der Waals interaction coefficient between layers, as applied to the vibrations of double-walled carbon nanotubes (DWCNTs); specifically, the evaluation of the natural frequencies is achieved through Ru’s and He’s formulations. The actual discrete DWCNT is modelled by means of a couple of concentric equivalent continuous thin cylindrical shells, where Donnell shell theory is adopted to obtain strain-displacement relationships. In order to take into account the chirality effect of DWCNT, an anisotropic elastic shell model is considered. Simply supported boundary conditions are imposed and the Rayleigh–Ritz method is used to obtain approximate natural frequencies and mode shapes. A parametric analysis considering different values of diameters and numbers of waves along longitudinal and circumferential directions is performed by adopting Ru’s and He’s formulations. From the comparisons, it is evident that Ru’s formulation provides unsatisfactory results for relatively low values of diameters and relatively high numbers of circumferential waves with respect to the more accurate He’s formulation. This behaviour is observed for every number of longitudinal half-waves. Therefore, Ru’s formulation cannot be used for the vibration modelling of DWCNTs in a large range of diameters and wavenumbers.

1. Introduction

Carbon nanotubes (CNTs) were discovered in 1991 inside the NEC Corporation laboratories by S. Iijima, who studied the synthesis of fullerenes and first prepared a new type of carbon structures, referred to as CNTs, which were described as “helical microtubules of graphitic carbon” [1]. Since CNTs have very high Young’s modulus and tensile strength, together with very small diameter, then they can reach the natural frequencies of the THz order and therefore can be used as ultra-high sensitivity resonators in many electro-mechanical devices, such as sensors, charge detectors and oscillators [2]. In addition, the extremely high electrical and thermal conductivity of CNTs, together with the very high transparency in a visible range, makes them promising candidates for innovative applications including heat exchangers and energy conversion systems, with the role of advanced catalyst supports or electrodes in solar and fuel cells [3,4].
Due to their relevant applications in several mechanical components, vibrations of CNTs have attracted much attention from Researchers. Specifically, CNT vibrations have been studied by means of experiments, atomistic mechanic simulations and continuum mechanic models. Resonant Raman spectroscopy (RRS) starts from the experimental measurement of CNT diameter by means of atomic force microscopy and investigates the atomic structure, chirality and natural frequencies of CNTs [5,6,7]. However, due to high technological complexity and costs, experimental methods cannot be considered efficient approaches to studying the mechanical behaviour of CNTs. Molecular dynamics (MD) simulations consider CNT atoms as interacting point-like masses, where the vibrations of the free atoms are recorded for a certain time duration at a fixed temperature [8,9,10]. However, due to the high computational effort, atomistic mechanics methods cannot be easily applied to the structural simulation of CNTs, which incorporate a large number of carbon atoms. Since theoretical models based on continuum mechanics are computationally more efficient than MD simulations and do not present the technological complexity and high costs of RRS, then continuous elastic models are the most commonly adopted methods for the study of CNT vibrations. In these models, it is assumed that the effective discrete structure of CNTs can be replaced by an equivalent homogeneous elastic structure with a continuous distribution of mass and stiffness, not considering their intrinsic atomic nature and so reducing the number of degrees of freedom [11,12,13].
In order to carefully investigate the validity of continuum mechanic approaches for the study of CNT vibrations, in the cases when the actual discrete CNT is modelled via an equivalent continuous thin cylindrical shell, some relevant open issues should be taken into account and properly discussed.
The first relevant issue to be taken into account in the modelling of CNTs as continuous elastic shells is given by anisotropy. CNTs are usually modelled as isotropic elastic shells [14,15,16,17,18,19,20,21]; however, they have a chirality-induced anisotropic behaviour due to their discrete nature that cannot be properly captured by using an isotropic shell model. To this aim, Chang et al. [22,23] developed a molecular mechanics model, called “stick-spiral model”, able to correctly predict the chirality and size-dependent elastic properties of CNTs. They found explicit expressions for longitudinal Young’s modulus and Poisson’s ratio, circumferential Young’s modulus and Poisson’s ratio, and longitudinal shear modulus in the case of chiral CNTs. This molecular-based anisotropic elastic shell model including chirality effects was validated via comparisons with RRS and MD data, proving that the classical relationship of isotropic elastic continuum mechanics between Young’s and shear modulus is not valid for CNTs. Ghavanloo and Fazelzadeh [24,25], starting from the results of [23], developed an anisotropic elastic shell model including chirality effects to study the vibration characteristics of CNTs. By applying the Flügge thin shell theory and complex method, they obtained the natural frequencies of radial breathing, torsional, longitudinal and radial modes under different CNT diameters and chiralities. Finally, they validated their model by means of comparisons with experimental RRS and numerical MD simulation data from the literature. In the present paper, the molecular mechanics-based anisotropic elastic shell model developed by Chang [22,23], and then refined by Ghavanloo and Fazelzadeh [24,25], for SWCNTs is extended to DWCNTs. Readers interested in deepening vibrations and stability of circular cylindrical shells can refer to the fundamental books of Leissa [26], Yamaki [27], Soedel [28] and Ventsel [29], where the strain-displacement relationships of the most important thin shell theories together with the related equations of motion and natural frequencies under different boundary conditions are reported. CNTs are divided into two main classes: single-walled carbon nanotubes (SWCNTs), given by graphene sheets rolled into cylinders, and multi-walled carbon nanotubes (MWCNTs), composed by concentric SWNTs, where the different cylinders are connected by van der Waals interaction forces.
In addition to anisotropy, in the specific case of MWCNTs, another relevant issue to be taken into consideration in the modelling of CNTs as continuous elastic shells is given by the simulation of van der Waals interaction forces between the concentric SWCNTs. The two most relevant formulations of van der Waals interaction forces in MWCNTs were developed by Ru and He. Both of them are based on the Lennard–Jones pair potential, which takes into account carbon–carbon potential depth, equilibrium separation distance and atom distance, where van der Waals interaction forces are obtained by deriving Lennard–Jones pair potential with respect to the carbon–carbon atom distance [30]. The pressure exerted on the ith SWCNT, as a function of the radial displacement, is obtained by integrating van der Waals interaction forces over the entire MWCNT.
Ru [31,32,33] proposed a linear relationship between pressure due to van der Waals interaction forces and radial displacement for the buckling and vibration analysis of MWCNTs in which the interaction coefficient is constant, i.e., it is not dependent on the radius of the individual SWCNT. In this model, the value of the interaction coefficient can be estimated using data given in [34], which are derived from the known van der Waals interaction between a carbon atom and a flat graphite sheet, and the van der Waals pressure between the inner and outer tube at any point is assumed to depend on the interlayer spacing at that point. By adopting Ru’s formulation, Elishakoff and Pentaras [35,36] obtained fundamental natural frequencies of DWCNTs considering different boundary conditions and shell theories by adopting an isotropic elastic shell model. In particular, the results of Ref. [35], obtained by considering simply supported boundary conditions and Donnell shell theory, are adopted to validate the model proposed in the present paper in the isotropic case for Ru’s formulation. Other relevant analytical studies on the linear vibrations of DWCNTs under different boundary conditions and geometries considering Ru’s formulation are reported in Refs. [37,38,39,40]. However, it is clear that this simplified formulation is not accurate in modelling van der Waals interaction forces.
In order to accurately describe van der Waals interactions in MWCNTs, He [41,42,43] proposed a linear relationship between pressure due to van der Waals interaction forces and radial displacement where the interaction coefficient is radius dependent. By considering this refined formulation, natural frequencies of MWCNTs were analysed for different geometries and chiralities obtaining an excellent agreement with MD simulations. In particular, the results of Ref. [42], obtained by considering simply supported boundary conditions and the Donnell shell theory, are adopted to validate the model proposed in the present paper in the isotropic case for He’s formulation. Strozzi and Pellicano [44,45,46,47] studied linear and nonlinear vibrations of MWCNTs considering He’s formulation. The effect of geometry, boundary conditions and wavenumbers along longitudinal and circumferential directions on the natural frequencies were analysed. Applicability and limitations of the Donnell shell theory, with respect to the more accurate Sanders shell theory, were considered. The effect of nonlinear van der Waals interaction forces on the nonlinear response was investigated. Other interesting results on the linear vibrations of DWCNTs considering He’s formulation are shown in Refs. [48,49,50], where the influence of boundary conditions and surrounding elastic medium on the resonant natural frequencies in case of noncoaxial vibration is studied. Moreover, in Refs. [51,52,53,54,55], some recent studies on vibrations and stability of DWCNTs, based on molecular dynamics simulations, non-local effect and anisotropic elastic model, are reported.
In the present paper, the comparison is conducted between Ru’s and He’s formulations of the van der Waals interaction coefficient for the evaluation of the natural frequencies of DWCNTs. The actual discrete DWCNT is modelled by means of a couple of concentric equivalent continuous cylindrical shells, where the Donnell shell theory is considered to obtain strain-displacement relationships. In order to take into account the chirality effects of DWCNTs, an anisotropic elastic shell model is used. Simply supported boundary conditions are imposed and the Rayleigh–Ritz method is adopted to obtain approximate natural frequencies and mode shapes. In the beginning, the results of the proposed model in the isotropic form are compared with data retrieved from the literature to check the validity of the shape functions used in the modal expansions of the displacement field. Then, the values of van der Waals interaction coefficient from Ru’s and He’s formulations are compared for different DWCNT diameters to verify the higher accuracy of He’s formulation in simulating the actual anisotropic behaviour of DWCNTs. Starting from these results, considering an anisotropic elastic shell model and adopting Donnell shell theory, natural frequencies of DWCNTs with different diameters and wavenumber along longitudinal and circumferential directions obtained via Ru’s and He’s formulations are compared. Assuming the more accurate He’s formulation as a reference, the applicability and limitations of Ru’s formulation for the vibration modelling of DWCNTs are investigated.

2. Donnell Shell Theory for DWCNTs

In this paper, the actual discrete DWCNT shown in Figure 1a is modelled by means of a couple of concentric equivalent continuous cylindrical shells with van der Waals interactions between the layers. In Figure 1b,c, a single continuous cylindrical shell with radius R , length L and thickness h is shown; in order to describe the displacement field, a cylindrical coordinate system ( O , x , θ , z ) is adopted, due to the axisymmetric symmetry of the cylindrical shell, where the origin O of the reference system is located at the centre of one end of the shell. Three displacements are present: longitudinal u ( x , θ , t ) , circumferential v ( x , θ , t ) and radial w ( x , θ , t ) , where the radial displacement w is assumed to be positive outward; ( x , θ ) are the longitudinal and angular coordinates of an arbitrary point on the middle surface of the shell; z is the radial coordinate along the thickness h of the shell; t is the time variable. Starting from this continuum modelling, in the following, displacement field and strain-displacement relations for DWCNTs will be derived.

2.1. Displacement Field

In order to avoid very fast time scales, which are due to the infinitesimal dimensions of CNTs, and therefore very hard convergence problems in the numerical integration, in this work, all dimensionless parameters are considered.
The dimensionless displacement field ( u ˜ i , v ˜ i , w ˜ i ) of the i-th cylindrical shell is written as [44]:
u ˜ i = u i R i v ˜ i = v i R i i = 1 , 2 w ˜ i = w i R i
where ( u i , v i , w i ) is the corresponding dimensional displacement field and R i is the radius of the i -th shell.

2.2. Strain-Displacement Relationships

In this paper, the Donnell shell theory is adopted to obtain linear strain-displacement relationships.
The dimensionless middle surface strains of the i -th cylindrical shell ( ε ˜ x , 0 , i , ε ˜ θ , 0 , i , γ ˜ x θ , 0 , i ) are written in the form [44]:
ε ˜ x , 0 , i = α i u ˜ i η ε ˜ θ , 0 , i = v ˜ i θ + w ˜ i i = 1 , 2 γ ˜ x θ , 0 , i = u ˜ i θ + α i v ˜ i η
where η = x / L is the dimensionless longitudinal coordinate of the shell and α i = R i / L .
The dimensionless middle surface changes in curvature and torsion of the i -th cylindrical shell ( k ˜ x , i , k ˜ θ , i , k ˜ x θ , i ) are written as [44]:
k ˜ x , i = α i 2 2 w ˜ i η 2 k ˜ θ , i = 2 w ˜ i θ 2 i = 1 , 2 k ˜ x θ , i = 2 α i 2 w ˜ i η θ
The dimensionless strains in a generic point of the external surface of the i -th cylindrical shell ( ε ˜ x , i , ε ˜ θ , i , γ ˜ x θ , i ) are related to the corresponding dimensionless middle surface strains and changes in curvature and torsion by means of the following relationships [44]:
ε ˜ x , i = ε ˜ x , 0 , i + ζ i k ˜ x , i ε ˜ θ , i = ε ˜ θ , 0 , i + ζ i k ˜ θ , i i = 1 , 2 γ ˜ x θ , i = γ ˜ x θ , 0 , i + ζ i k ˜ x θ , i
where ζ i = z i / R i is the dimensionless radial coordinate of the i -th cylindrical shell.

3. Anisotropic Elastic Shell Model for DWCNTs

Actually, CNTs are discrete systems composed of carbon atoms with covalent bonds, and therefore, they are intrinsically non-isotropic. Accordingly, it was proven that the use of an isotropic elastic shell model can lead to relevant errors in CNT vibration analysis, in particular to natural frequencies very different from those obtained by means of resonant Raman experimental spectroscopy or molecular dynamics numerical simulations, see Ref. [24]. On the basis of these findings, in the present paper, an anisotropic elastic shell model is adopted for DWCNT vibration analysis.

3.1. Stress–Strain Relationships

First of all, it should be stressed that the present study is conducted by assuming the plane stress hypothesis ( σ z = 0 ), which is peculiar to thin-walled structures.
In the case of anisotropic elastic material, the dimensionless stresses ( σ ˜ x , i , σ ˜ θ , i , τ ˜ x θ , i ) of the i -th cylindrical shell are related to the corresponding dimensionless strains by means of the following constitutive Equation [46]:
σ ˜ x , i = Y ˜ 11 , i ε ˜ x , i + Y ˜ 12 , i ε ˜ θ , i + Y ˜ 13 , i γ ˜ x θ , i σ ˜ θ , i = Y ˜ 21 , i ε ˜ x , i + Y ˜ 22 , i ε ˜ θ , i + Y ˜ 23 , i γ ˜ x θ , i i = 1 , 2 τ ˜ x θ , i = Y ˜ 31 , i ε ˜ x , i + Y ˜ 32 , i ε ˜ θ , i + Y ˜ 33 , i γ ˜ x θ , i
where the dimensionless material parameters Y ˜ j k , i , which substitute Young’s modulus and Poisson’s ratio of the corresponding isotropic model, represent anisotropic surface elastic constants, in the form [46]:
Y ˜ j k , i = ( G ˜ l j , i G ˜ l k , i + 2 μ H ˜ l j , i H ˜ l k , i ) j , k , l = 1 , 2 , 3   ( sum   over   l ) i = 1 , 2
which allow the chirality effects to be considered in the anisotropic elastic constitutive Equation (5), where ( G ˜ l j , i , G ˜ l k , i , H ˜ l j , i , H ˜ l k , i ) are dimensionless constants and μ is a dimensionless constant parameter [46]:
μ = K θ K ρ a 2
where K ρ and K θ are force constants associated with the stretching and angular distortion of the carbon–carbon bond, respectively, which can be derived from molecular mechanics analyses or experimental data, see Ref. [23] for more details, and a is the carbon–carbon bond length.
The corresponding matrices G ˜ i and H ˜ i can be calculated as follows [46]:
G ˜ i = B ˜ i 1 ( I ˜ D ˜ i F ˜ i ) ,           H ˜ i = Q ˜ i F ˜ i i = 1 , 2
where I ˜ is the identity matrix, matrix F ˜ i is given by [46]:
F ˜ i = [ U ˜ i B ˜ i 1 D ˜ i ( 2 μ V ˜ i A ˜ i + W ˜ i ) ] 1 U ˜ i B ˜ i 1 i = 1 , 2
and matrices ( A ˜ i , B ˜ i , D ˜ i , U ˜ i , V ˜ i , W ˜ i , Q ˜ i ) are given by [46]:
A ˜ i = { A ˜ j k , i } = { cos ω j p , i cos ω k p , i } j , k , p = 1 , 2 , 3   ( sum   over   p ) i = 1 , 2
B ˜ i = 1 3 r i 2 + r i s i + s i 2 ( ( 2 r i + s i ) cos ϕ 1 , i ( r i s i ) cos ϕ 2 , i ( r i + 2 s i ) cos ϕ 3 , i 3   s i sin ϕ 1 , i 3 ( r i s i ) sin ϕ 2 , i 3   r i sin ϕ 3 , i ( 2 r i + s i ) sin ϕ 1 , i ( r i s i ) sin ϕ 2 , i ( r i + 2 s i ) sin ϕ 3 , i )
D ˜ i = 1 3 r i 2 + r i s i + s i 2 ( ( 2 r i + s i ) sin ϕ 1 , i ( r i s i ) sin ϕ 2 , i ( r i + 2 s i ) sin ϕ 3 , i 3   s i cos ϕ 1 , i 3 ( r i + s i ) cos ϕ 2 , i 3   r i cos ϕ 3 , i ( 2 r i + s i ) cos ϕ 1 , i ( r i s i ) cos ϕ 2 , i ( r i + 2 s i ) cos ϕ 3 , i )
U ˜ i = ( sin ϕ 1 , i sin ϕ 2 , i sin ϕ 3 , i cos ϕ 1 , i cos ϕ 2 , i cos ϕ 3 , i s i cos ϕ 1 , i ( r i + s i ) cos ϕ 2 , i r i cos ϕ 3 , i )
V ˜ i = ( cos ϕ 1 , i cos ϕ 2 , i cos ϕ 3 , i sin ϕ 1 , i sin ϕ 2 , i sin ϕ 3 , i 0 0 0 )
W ˜ i = ( 0 0 0 0 0 0 s i sin ϕ 1 , i ( r i + s i ) sin ϕ 2 , i r i sin ϕ 3 , i )
Q ˜ i = { Q ˜ j k , i } = { cos ω k j , i }                   j , k = 1 , 2 , 3 i = 1 , 2
where [46]:
cos ω j k , i = { ( cos ϕ j , i sin ϕ p , i cos φ k , i sin ϕ j , i cos ϕ p , i ) / sin θ k , i 0                       j k p j = k
and ( r i , s i ) are the chirality indices of the i -th carbon nanotube, which describe its rolling angle (i.e., direction), and therefore, define its radius by means of the formulation [18]:
R i = 3   a 2 π   r i 2 + r i s i + s i 2
The structural parameters of the i -th SWCNT, i.e., chiral angles ( ϕ 1 , i , ϕ 2 , i , ϕ 3 , i ) and torsion angles ( φ 1 , i , φ 2 , i , φ 3 , i ) , which can be calculated by means of Equation [47]:
ϕ 1 , i = arccos 2 r i + s i 2 r i 2 + r i s i + s i 2 ϕ 2 , i = 4 π 3 + ϕ 1 , i ϕ 3 , i = 2 π 3 + ϕ 1 , i
φ 1 , i = π r i 2 + r i s i + s i 2 cos ϕ 1 , i φ 2 , i = π r i 2 + r i s i + s i 2 cos ( π 3 + ϕ 1 , i ) φ 3 , i = π r i 2 + r i s i + s i 2 cos ( π 3 ϕ 1 , i )
and bond angles ( θ 1 , i , θ 2 , i , θ 3 , i ) , which are written in the following form [47]:
cos θ j , i = sin ϕ k , i sin ϕ p , i cos φ j , i + cos ϕ k , i cos ϕ p , i j , k , p = 1 , 2 , 3   ( sum   over   p )                                     i = 1 , 2
are illustrated in Figure 2, see Ref. [22] for more details.
It should be stressed that the chirality of a generic SWCNT is indicated by the chiral angle ϕ 1 , see Equation (19) and Figure 2. In particular, the two limiting cases of carbon nanotubes are ( r , 0 ) , which gives chiral angle ϕ 1 = 0 and is defined “zigzag” CNT, and ( r , r ) , which gives chiral angle ϕ 1 = π / 6 and is defined “armchair” CNT. These two definitions are based on the different geometry of the carbon atoms along the circumferential direction of the nanotube. Moreover, zigzag and armchair CNTs are considered as achiral nanotubes, due to their high geometric symmetry, whereas SWCNTs with a chiral angle within the range 0 < ϕ 1 < π / 6 are referred to as “chiral” CNTs, see Ref. [22] for details.

3.2. Force and Moment Resultants

Starting from the anisotropic elastic constitutive Equation (5), the corresponding dimensionless in-plane force resultants per unit length of the i -th shell ( N ˜ x , i , N ˜ θ , i , N ˜ x θ , i ) can be written as a function of the corresponding dimensionless middle surface strains and of the dimensionless material parameters Y ˜ j k , i in the form:
N ˜ x , i = β i / 2 + β i / 2 σ ˜ x , i d ζ i = Y ˜ 11 , i ε ˜ x , 0 , i + Y ˜ 12 , i ε ˜ θ , 0 , i + Y ˜ 13 , i γ ˜ x θ , 0 , i N ˜ θ , i = β i / 2 + β i / 2 σ ˜ θ , i d ζ i = Y ˜ 21 , i ε ˜ x , 0 , i + Y ˜ 22 , i ε ˜ θ , 0 , i + Y ˜ 23 , i γ ˜ x θ , 0 , i N ˜ x θ , i = β i / 2 + β i / 2 τ ˜ x θ , i d ζ i = Y ˜ 31 , i ε ˜ x , 0 , i + Y ˜ 32 , i ε ˜ θ , 0 , i + Y ˜ 33 , i γ ˜ x θ , 0 , i
where β i = h / R i is the thickness ratio of the i -th cylindrical shell.
Similarly, the dimensionless moment resultants per unit length of the i -th cylindrical shell ( M ˜ x , i , M ˜ θ , i , M ˜ x θ , i ) can be written as a function of the corresponding dimensionless middle surface changes in curvature and torsion and of the dimensionless material parameters Y ˜ j k , i in the form:
M ˜ x , i = β i / 2 + β i / 2 σ ˜ x , i ζ i d ζ i = β i 12 ( Y ˜ 11 , i k ˜ x , i + Y ˜ 12 , i k ˜ θ , i + Y ˜ 13 , i k ˜ x θ , i ) M ˜ θ , i = β i / 2 + β i / 2 σ ˜ θ , i ζ i d ζ i = β i 12 ( Y ˜ 21 , i k ˜ x , i + Y ˜ 22 , i k ˜ θ , i + Y ˜ 23 , i k ˜ x θ , i ) M ˜ x θ , i = β i / 2 + β i / 2 τ ˜ x θ , i ζ i d ζ i = β i 12 ( Y ˜ 31 , i k ˜ x , i + Y ˜ 32 , i k ˜ θ , i + Y ˜ 33 , i k ˜ x θ , i )
Finally, the dimensionless transverse shear forces per unit length of the i -th cylindrical shell ( Q ˜ x , i , Q ˜ θ , i ) , which are obtained directly from the expressions of the dimensionless moment resultants per unit length, can be written as:
Q ˜ x , i = β i 2 12 ( α i ( Y ˜ 11 , i k ˜ x , i η + Y ˜ 12 , i k ˜ θ , i η + Y ˜ 13 , i k ˜ x θ , i η ) + Y ˜ 31 , i k ˜ x , i θ + Y ˜ 32 , i k ˜ θ , i θ + Y ˜ 33 , i k ˜ x θ , i θ ) Q ˜ θ , i = β i 2 12 ( α i ( Y ˜ 31 , i k ˜ x , i η + Y ˜ 32 , i k ˜ θ , i η + Y ˜ 33 , i k ˜ x θ , i η ) + Y ˜ 21 , i k ˜ x , i θ + Y ˜ 22 , i k ˜ θ , i θ + Y ˜ 23 , i k ˜ x θ , i θ )
The previous equations will be used in the following in the expressions of the natural boundary conditions imposed on the i -th cylindrical shell, which involve forces and moments resultants.

3.3. Elastic Strain Energy

In the case of homogeneous, anisotropic and elastic material, assuming plane stress hypothesis ( σ z = 0 ) , the dimensionless elastic strain energy U ˜ i of the i -th cylindrical shell is expressed as a function of the dimensionless stresses ( σ ˜ x , i , σ ˜ θ , i , τ ˜ x θ , i ) and strains ( ε ˜ x , i , ε ˜ θ , i , γ ˜ x θ , i ) as follows [44]:
U ˜ i = 1 2 β i 0 1 0 2 π β i / 2 + β i / 2 ( σ ˜ x , i ε ˜ x , i + σ ˜ θ , i ε ˜ θ , i + τ ˜ x θ , i γ ˜ x θ , i ) d η d θ d ζ
By inserting Equations (4) and (5) into Equation (25), the dimensionless elastic strain energy U ˜ i can be rewritten as a function of the dimensionless middle surface strains ( ε ˜ x , 0 , i , ε ˜ θ , 0 , i , γ ˜ x θ , 0 , i ) and changes in curvature and torsion ( k ˜ x , i , k ˜ θ , i , k ˜ x θ , i ) , and of the dimensionless material parameters Y ˜ j k , i , and it becomes:
                U ˜ i = 1 2 0 1 0 2 π ( Y ˜ 11 , i ε ˜ x , 0 , i 2 + Y ˜ 12 , i ε ˜ x , 0 , i ε ˜ θ , 0 , i + Y ˜ 13 , i ε ˜ x , 0 , i γ ˜ x θ , 0 , i + Y ˜ 21 , i ε ˜ x , 0 , i ε ˜ θ , 0 , i +                 + Y ˜ 22 , i ε ˜ θ , 0 , i 2 + Y ˜ 23 , i ε ˜ θ , 0 , i γ ˜ x θ , 0 , i + Y ˜ 31 , i ε ˜ x , 0 , i γ ˜ x θ , 0 , i + Y ˜ 32 , i ε ˜ θ , 0 , i γ ˜ x θ , 0 , i + Y ˜ 33 , i γ ˜ x θ , 0 , i 2 ) d η d θ + β i 2 24 0 1 0 2 π ( Y ˜ 11 , i k ˜ x , i 2 + Y ˜ 12 , i k ˜ x , i k ˜ θ , i + Y ˜ 13 , i k ˜ x , i k ˜ x θ , i + Y ˜ 21 , i k ˜ x , i k ˜ θ , i                 + Y ˜ 22 , i k ˜ θ , i 2 + Y ˜ 23 , i k ˜ θ , i k ˜ x θ , i + Y ˜ 31 , i k ˜ x , i k ˜ x θ , i + Y ˜ 32 , i k ˜ θ , i k ˜ x θ , i + Y ˜ 33 , i k ˜ x θ , i 2 ) d η d θ
It should be observed that, in Equation (26), the first term on the right-hand side is associated to the middle surface strains of the shell, and therefore, it represents the stretching (tensile/compression) energy, while the second term is associated to the middle surface changes in curvature and torsion of the shell, and therefore it represents the bending (flexural) energy [18].
The dimensionless elastic strain energy of a DWCNT is given by [44]:
U ˜ = i = 1 2 δ i U ˜ i
where δ i = R i / R 1 and R 1 is the radius of the inner SWCNT.

3.4. Kinetic Energy

The most important issue in the determination of the dimensionless kinetic energy is the definition of the dimensionless time variable. In the present work, the dimensionless time variable τ is obtained by introducing the reference dimensional circular frequency ω 0 , in the form:
τ = ω 0 t ω 0 = Y ρ m h R 1 2
where t is the corresponding dimensional time variable, ρ m is the mass density of the material and
Y = 2 K ρ 3 3
is a reference dimensional surface elastic constant.
Starting from the definition of dimensionless time variable τ (28), the dimensionless velocity field of the i -th cylindrical shell ( u ˜ i , v ˜ i , w ˜ i ) is written as [44]:
u ˜ i = d u ˜ i d τ = u ˙ i R i ω 0 v ˜ i = d v ˜ i d τ = v ˙ i R i ω 0 i = 1 , 2 w ˜ i = d w ˜ i d τ = w ˙ i R i ω 0
where ( u ˙ i , v ˙ i , w ˙ i ) is the corresponding dimensional velocity field.
Neglecting the rotary inertia effect, the dimensionless kinetic energy of the i -th cylindrical shell is given by [44]:
T ˜ i = 1 2 δ i 2 0 1 0 2 π ( u ˜ i 2 + v ˜ i 2 + w ˜ i 2 ) d η d θ i = 1 , 2
It can be noted that, differently from the elastic strain energy (26), the kinetic energy (31) does not depend on the material properties, as the dimensionless material parameters Y ˜ j k , i are not present and the mass density ρ m is included in the definition of the reference circular frequency ω 0 ), and therefore, it assumes the same expression for both isotropic and anisotropic elastic material.
The dimensionless kinetic energy of a DWCNT is given by [44]:
T ˜ = i = 1 2 δ i T ˜ i

3.5. Van der Waals Interaction Forces

In the present paper, the actual discrete DWCNT is modelled by means of a couple of concentric equivalents continuous thin cylindrical shells, and van der Waals interaction energy potential between the two layers can be estimated by the Lennard–Jones model [30]:
V L J ( a ) = 4 ε ( ( σ a ) 12 ( σ a ) 6 )
where ε is the carbon–carbon potential depth, σ is the carbon–carbon equilibrium separation distance and a is the carbon–carbon bond length (i.e., the distance between the interacting atoms).
The van der Waals interaction force is obtained by deriving Lennard–Jones pair potential (33), with respect to the carbon–carbon bond length [41]:
F ( a ) = d V L J ( a ) d a = 24 ε σ ( 2 ( σ a ) 13 ( σ a ) 7 )
where the negative value represents the attractive force between a pair of atoms, whereas the positive value represents the corresponding repulsive force.
Since only the infinitesimal buckling of CNTs is of interest, then the van der Waals force (34) can be estimated via the Taylor expansion to the first order around the equilibrium position prior to buckling, in the form [41]:
F ( a ) = F ( a 0 ) + d F ( a 0 ) d a 0 ( a a 0 )
where a 0 is the initial carbon–carbon atom distance of the two layers prior to bucking.
By substituting Equation (34) into expansion (35), it is derived [41]:
F ( a ) = 24 ε σ ( 2 ( σ a 0 ) 13 ( σ a 0 ) 7 ) 24 ε σ 2 ( 26 ( σ a 0 ) 14 7 ( σ a 0 ) 8 ) ( a a 0 )
At the equilibrium position prior to buckling, the initial van der Waals interaction force F ( a 0 ) is very small, and therefore, can be neglected.
By integrating Equation (36) over the entire DWCNT, it is obtained the pressure p i exerted on the ith SWCNT due to the van der Waals interactions between the layers ( i , j ) , expressed as a function of the radial displacements ( w i , w j ) [41]:
p i ( η , θ ) = c i j ( w i w j ) i , j = 1 , 2 i j
where c i j is the van der Waals interaction coefficient between the layers ( i , j ) .
In the literature, several formulations for the van der Waals interaction coefficient are reported; in particular, it can be expressed by adopting Ru’s and He’s formulations.

3.5.1. Van der Waals Interaction Coefficient: Ru’s Formulation

In Ru’s formulation, the van der Waals interaction coefficient between the two layers is estimated starting from the energy potential per unit area between a carbon atom and a flat graphite sheet, as reported by Girifalco and Lad [34], and it is expressed in the form [31]:
c R u = 3.2 × 10 2   erg / cm 2 0.16 a 2
where 1   erg / cm 2 = 10 3   N / m .
Note that the equilibrium distance between a carbon atom and a flat graphite sheet is around 0.34 nm [34], so the initial pressure between the inner and outer tubes is zero if the initial interlayer spacing is 0.34 nm. In this case, any increase (decrease) in the interlayer spacing causes an attractive (repulsive) van der Waals interaction, and then coefficient c should be negative (positive) [32].
From expression (38), it can be observed that, in Ru’s formulation, the van der Waals interaction coefficient represents a constant parameter, which is assumed to depend on the distance between the interacting atoms a , while it is not dependent on CNT diameter. Specifically, the curvature effect has been ignored because it is usually small for the MWNTs of the larger innermost radius [5,6,7].

3.5.2. Van der Waals Interaction Coefficient: He’s Formulation

In He’s formulation, the van der Waals interaction coefficient is derived via the direct integration of Equation (36), and it is expressed in the form [41]:
c i j , H e = ( 1001 π ε σ 12 3 a 4 E i j 13 1120 π ε σ 6 9 a 4 E i j 7 ) R j i , j = 1 , 2 i j
From expression (39), it can be observed that, in He’s formulation, the van der Waals interaction coefficient depends on the radius R j of the jth layer. In addition, for a specific DWCNT geometry, it is obtained that c i j c j i , i.e., the two coefficients are different, with the relationship c i j R i = c j i R j .
The elliptical integral E i j m of expression (39) is written as [41]:
E i j m = ( R j + R i ) m 0 π / 2 d θ ( 1 k i j cos 2 θ ) m / 2 i , j = 1 , 2 i j m = 7 , 13
and the geometric coefficient k i j of Equation (40) is written as [41]:
k i j = 4 R j R i ( R j + R i ) 2 i , j = 1 , 2 i j

3.5.3. Van der Waals Interaction Energy

The van der Waals interaction energy of the i-th cylindrical shell, which models the SWCNT, is expressed in dimensionless form as follows [44]:
V ˜ i = 1 2 δ i 0 1 0 2 π p ˜ i ( η , θ ) w ˜ i d η d θ i = 1 , 2
The dimensionless van der Waals interaction energy of a DWCNT is given by [44]:
V ˜ = i = 1 2 δ i V ˜ i

4. Vibration Modelling of DWCNTs

In the case of DWCNTs, a modal vibration can be formally written in terms of the three dimensionless displacements of the i -th cylindrical shell ( u ˜ i , v ˜ i , w ˜ i ) as [44]:
u ˜ i ( η , θ , τ ) = U ˜ i ( η , θ ) f ˜ i ( τ ) v ˜ i ( η , θ , τ ) = V ˜ i ( η , θ ) f ˜ i ( τ ) i = 1 , 2 w ˜ i ( η , θ , τ ) = W ˜ i ( η , θ ) f ˜ i ( τ )
where ( U ˜ i , V ˜ i , W ˜ i ) are the corresponding three dimensionless components of the modal shape and f ˜ i ( τ ) is the corresponding dimensionless time law, which is supposed to be the same for the three different displacements (i.e., synchronous motion hypothesis).
The three modal shape components ( U ˜ i , V ˜ i , W ˜ i ) can be expressed in different ways. Usually, they are expanded by means of a double series in terms of harmonic functions along the longitudinal and circumferential directions, see for example Ref. [35]. In the present work, they are expanded by means of a double mixed series in terms of m -th degree Chebyshev orthogonal polynomials T m * ( η ) = T m ( 2 η 1 ) along the longitudinal direction η and harmonic functions ( cos n θ , sin n θ ) along the circumferential direction θ [44]:
U ˜ i ( η , θ ) = m = 0 M u n = 0 N U ˜ i , m , n T m * ( η ) cos n θ V ˜ i ( η , θ ) = m = 0 M v n = 0 N V ˜ i , m , n T m * ( η ) sin n θ i = 1 , 2 W ˜ i ( η , θ ) = m = 0 M w n = 0 N W ˜ i , m , n T m * ( η ) cos n θ
where m is the number of longitudinal half-waves (maximum values ( M u , M v , M w ) ), n is the number of circumferential waves (maximum value N ) and ( U ˜ i , m , n , V ˜ i , m , n , W ˜ i , m , n ) are unknown coefficients, which can be obtained by imposing the boundary conditions.
This last is the most relevant difference due to the choice of Chebyshev polynomials instead of harmonic functions to model DWCNT vibrations along the longitudinal direction; actually, the boundary conditions are automatically satisfied by harmonic functions, while they are not automatically satisfied by Chebyshev polynomials, and therefore, in the second approach, they must be imposed on the unknown coefficients ( U ˜ i , m , n , V ˜ i , m , n , W ˜ i , m , n ) , as in the present work.

4.1. Boundary Conditions

In this paper, simply supported DWCNTs are considered. The corresponding boundary conditions on the two ends of the two SWCNTs are expressed in the form [45]:
v ˜ i ( η , θ ) = 0 w ˜ i ( η , θ ) = 0 N ˜ x , i ( η , θ ) = 0 M ˜ x , i ( η , θ ) = 0 η = 0 , 1       i = 1 , 2
which represent both geometric (on the displacements ( v ˜ i , w ˜ i ) ) and natural (on the force N ˜ x , i and moment M ˜ x , i resultants) conditions, where the dependency on the dimensionless time variable τ is not taken into account (modal vibration analysis).
By taking into account Equations (22) and (23) for the dimensionless force and moment resultants, the two natural boundary conditions (46) can be rewritten as a function of the dimensionless middle surface strains and changes in curvature and torsion of the i -th cylindrical shell in the following form:
N ˜ x , i ( η , θ ) = Y ˜ 11 , i ε ˜ x , 0 , i + Y ˜ 12 , i ε ˜ θ , 0 , i + Y ˜ 13 , i γ ˜ x θ , 0 , i = 0 M ˜ x , i ( η , θ ) = β i 12 ( Y ˜ 11 , i k ˜ x , i + Y ˜ 12 , i k ˜ θ , i + Y ˜ 13 , i k ˜ x θ , i ) = 0 η = 0 , 1       i = 1 , 2
Since the boundary conditions are related only to the dimensionless longitudinal coordinate η and do not depend on the circumferential coordinate θ (i.e., they are axisymmetric), then in the expressions (47) of the dimensionless force and moment resultants only the dimensionless middle surface strain and change in curvature along the longitudinal direction must be considered, and therefore, expressions (47) become:
N ˜ x , i ( η ) = Y ˜ 11 , i ε ˜ x , 0 , i = Y ˜ 11 , i α i u ˜ i η = 0 M ˜ x , i ( η ) = β i 12 Y ˜ 11 , i k ˜ x , i = β i 12 Y ˜ 11 , i α i 2 2 w ˜ i η 2 = 0 η = 0 , 1       i = 1 , 2
Starting from the expressions for the geometric boundary conditions (46) and the natural boundary conditions (48), the following equations in terms of the dimensionless displacements ( u ˜ i , v ˜ i , w ˜ i ) related to the longitudinal coordinate η of the i -th cylindrical shell are derived:
v ˜ i ( η ) = 0 w ˜ i ( η ) = 0 u ˜ i η ( η ) = 0 2 w ˜ i η 2 ( η ) = 0 η = 0 , 1       i = 1 , 2
By substituting expressions (49) into Equations (44) for the dimensionless displacements, the following conditions in terms of expansions (45) for the corresponding modal shape components are derived [45]:
V ˜ i ( η ) = m = 0 M v V ˜ i , m , n T m * ( η ) = 0 W ˜ i ( η ) = m = 0 M w W ˜ i , m , n T m * ( η ) = 0 U ˜ i , η ( η ) = m = 0 M u U ˜ i , m , n T m , η * ( η ) = 0 W ˜ i , η η ( η ) = m = 0 M w W ˜ i , m , n T m , η η * ( η ) = 0 η = 0 , 1 i = 1 , 2
where ( · ) , η = ( · ) / η and ( · ) , η η = 2 ( · ) / η 2 .
The linear algebraic system given by Equation (50) can then be solved analytically in terms of the coefficients ( U ˜ i , 1 , n , U ˜ i , 2 , n , V ˜ i , 0 , n , V ˜ i , 1 , n , W ˜ i , 0 , n , W ˜ i , 1 , n , W ˜ i , 2 , n , W ˜ i , 3 , n ) , for n [ 0 , N ] . Therefore, in the specific case of simply supported DWCNTs, eight different dimensionless coefficients for each SWCNT are obtained.

4.2. Rayleigh–Ritz Method

In the case of DWCNTs, the maximum number of variables needed to properly describe a generic vibration mode ( m , n ) , with m longitudinal half-waves and n circumferential waves, is given by N p = 2 × ( M u + M v + M w + 3 p ) , where 2 is the number of concentric SWCNTs, M u = M v = M w is the maximum number of longitudinal half-waves considered and p is the number of equations needed to satisfy the boundary conditions: as previously demonstrated, in the case of DWCNTs with simply supported boundary conditions it is imposed p = 8 .
A convergence analysis was performed in order to obtain the correct order of Chebyshev orthogonal polynomials (i.e., the maximum number of longitudinal half-waves) to be adopted: since it was obtained that M u = M v = M w = 11 allows the modal shapes with a number of longitudinal half-waves until m = 5 to be correctly modelled, and m = 5 is the maximum number of longitudinal half-waves considered in the numerical analyses of the present work, then M u = M v = M w = 11 provides accurate results with relatively reduced computational effort.
For a multi-mode analysis including different values of circumferential waves n , the number of degrees of freedom of the mechanical system is given by N max = N p × ( N + 1 ) , where N is the maximum number of circumferential waves considered.
Equations (44) are inserted in the expressions of elastic strain energy (26), kinetic energy (31) and van der Waals interaction energy (42) in order to calculate the value of Rayleigh quotient R ( q ˜ ) , where q ˜ is a vector containing all the unknown coefficients of expansions (45), expressed as [46]:
q ˜ = [ U ˜ i , m , n V ˜ i , m , n W ˜ i , m , n ] i = 1 , 2
After imposing the stationarity on the Rayleigh quotient, the following classic eigenvalue problem in dimensionless form is obtained [46]:
( ω ˜ 2 M ˜ + K ˜ ) q ˜ = 0
which provides approximate dimensionless circular frequencies (eigenvalues ω ˜ j ) and modal shapes (eigenvectors q ˜ j ), with j = ( 1 , 2 , , N max ) , where M ˜ and K ˜ are the dimensionless mass and stiffness matrices, respectively.
The approximate modal shape of the j -th mode of the i -th shell is provided by the expansions (45), where coefficients ( U ˜ i , m , n , V ˜ i , m , n , W ˜ i , m , n ) are replaced with coefficients ( U ˜ i , m , n ( j ) , V ˜ i , m , n ( j ) , W ˜ i , m , n ( j ) ) , which are the components of the j -th eigenvector q ˜ of Equation (52), and the vector function [47]:
W ˜ ( j ) ( η , θ ) = [ U ˜ i ( j ) ( η , θ ) V ˜ i ( j ) ( η , θ ) W ˜ i ( j ) ( η , θ ) ] i = 1 , 2
is the approximation of the j -th eigenfunction vector of the original problem.
Finally, it must be underlined that, in order to obtain the corresponding dimensional natural frequencies f j of DWCNTs, the dimensionless circular frequencies (eigenvalues) ω ˜ j have to be multiplied by the dimensional constant term ω 0 / 2 π (with dimensions s−1).

5. Numerical Results

In the present work, natural frequencies of DWCNTs obtained via Ru’s (38) and He’s (39) formulations of the van der Waals interaction coefficient between the two concentric SWCNTs are compared. The Donnell shell theory is applied to obtain strain-displacement relationships. An anisotropic elastic shell model is adopted to take into account the chirality effects of CNTs. Simply supported boundary conditions are considered. Vibration modes with different numbers of waves along the longitudinal and circumferential directions are analysed. DWCNTs with different values of diameters are investigated.
In Table 1 the values of carbon–carbon (C-C) bond parameters ( a , k ρ , k θ ) , C-C distance parameters ( ε , σ ) and CNT equivalent parameters ( h , ρ m ) , retrieved from the pertinent literature, are reported.
In particular, parameters k ρ and k θ , which denote force constants related to the variance of C-C bond length a and angle θ , respectively, and are used in the present work to compute the anisotropic elastic properties of CNTs, were derived in the molecular mechanics “stick-spiral model” developed by Chang et al. [22].
Moreover, parameters ε and σ , which denote C-C potential depth and equilibrium distance, respectively, and are used in this work to simulate the van der Waals interaction forces, were determined by He et al. [41] via buckling analysis of MWCNTs.
Finally, parameter h , which represents CNT equivalent thickness, and is adopted in the present analysis to study the actual discrete DWCNT as a couple of concentric equivalent continuous cylindrical shells, was derived from MD simulations of CNT energy by Yakobson et al. [14].

5.1. Validation of the Proposed Model in the Isotropic Case

In order to justify the validity of the proposed model, let us consider the corresponding isotropic case, in which the SWCNT is assumed to be an isotropic elastic thin shell, with tensile rigidity C (in-plane stiffness), bending rigidity D (flexural stiffness), surface density ρ s (mass density ρ m per unit lateral area h ) and Poisson’s ratio ν reported in Table 2, where chirality effects are neglected, see Ref. [14] for more details.
To this purpose, in this Section, the results of the present elastic shell model in isotropic form, which considers Chebyshev orthogonal polynomials along the longitudinal direction and harmonic functions along the circumferential direction, are compared with the ones of the isotropic elastic shell models of Refs. [35,42], which adopt harmonic functions along both longitudinal and circumferential directions, and use Ru’s and He’s formulations for van der Waals interaction coefficient, respectively.
This check is important because there can be relevant differences between the results due to the different shape functions used in the modal expansions of the displacement field, and the accuracy of the natural frequencies depends strongly on the choice of modal shape functions. Moreover, it must be stressed that, in order to perform a correct analysis, the same thin shell theory, i.e., the Donnell shell theory, is adopted.
In Table 3 and Table 4, the natural frequencies of a simply supported DWCNT with inner radius R 1 = 5   nm and aspect ratio L / R 2 = 10 are reported. The Donnell shell theory and isotropic elastic shell model are adopted. The van der Waals interaction coefficient is expressed via Ru’s (Table 3) and He’s (Table 4) formulations. Circumferential flexure modes ( n = 2 ) are analysed. Comparisons between present and Ref. [35] models (see Table 3), and present and Ref. [42] models (see Table 4), are performed.
From Table 3 and Table 4, it is observed that the present elastic shell model in isotropic form is accurate both for Ru’s and for He’s formulations, as it gives natural frequencies with less than 1% difference, and therefore, the validity of the shape functions adopted in the modal expansions of the displacement field is confirmed.

5.2. Ru’s vs. He’s Formulations of the van der Waals Interaction Coefficient for DWCNTs

Before comparing the natural frequencies of DWCNTs obtained by considering an anisotropic elastic shell model via Ru’s and He’s formulations, it is useful to investigate the values of the corresponding interaction coefficients under different values of CNT diameter, in order to justify the assumption of He’s formulation as a reference in the comparisons.
In Figure 3, the effect of the DWCNT inner radius R 1 on the value of the van der Waals interaction coefficient c obtained by assuming the C-C bond parameters ( a , ε , σ ) of Table 1 and considering Ru’s (38) and He’s (39) formulations is shown.
From this Figure, it can be observed that, for a very large inner radius ( R 1 > 10   nm ) , the two different interaction coefficients c 12 and c 21 deriving from He’s formulation tend asymptotically to the same value ( c 1.09 ) . This is due to the fact that as the inner radius increases, the outer radius also increases, while the distance between the two layers remains constant; therefore, the percentage difference between the inner and outer radius decreases, and the values of the two coefficients converge.
Moreover, it can be observed that this value is quite close to the constant one assumed by the interaction coefficient from Ru’s formulation ( c 1 ) . This is because, for a sufficiently large inner radius, the interaction coefficient no longer depends on the radius, but it assumes a constant value, which is similar to Ru’s and He’s formulations.
This behaviour, however, is not important for the modelling of CNTs, since it lies outside the range commonly assumed for the radius of CNTs. In fact, it is reported that “a single-wall nanotube is defined by a cylindrical graphene sheet with a diameter of about 0.7–10.0 nm, though most of the observed single-wall nanotubes have diameters < 2 nm”, see Ref. [24]. For this reason, it is important to focus the attention on the comparison between the interaction coefficients, and the corresponding natural frequencies of DWCNTs, from Ru’s and He’s formulations in this range of inner radius.
From Figure 3 it can be observed that, for a sufficiently small inner radius ( R 1 < 1   nm ) , such as commonly assumed by CNTs, the two interaction coefficients c 12 and c 21 from He’s formulation have very different values. This is because, in this range, the distance between the layers is significant with respect to the relative radii; therefore, the percentage difference between the inner and outer radius is very high, and the values of the two interaction coefficients diverge.
Moreover, it can be observed that these two values are very different compared to the constant one assumed by the interaction coefficient derived from Ru’s formulation. This is due to the fact that, for a sufficiently small inner radius, the interaction coefficient strongly depends on the radius itself, and therefore it assumes very distant values for Ru’s and He’s formulations.
As stated above, He’s formulation, with respect to Ru’s one, taking into account the effect of CNT diameter on the value of the van der Waals interaction coefficient, is able to simulate more correctly the actual CNT behaviour, which is anisotropic, that is the elastic properties of the material depend on CNT chirality and size, and therefore is more accurate.
This is especially true for a relatively small CNT diameter ( < 2   nm ) , where it is known from the literature that the surface elastic constants assume values that strongly depend on chirality and size, and therefore it is necessary to adopt an anisotropic elastic model. On the other hand, the surface elastic constants tend asymptotically to constant values coinciding with those of the relative constants of graphite as CNT diameter increases, and therefore it is possible to adopt an isotropic elastic model, see Ref. [25] for more details.
In order to take into account chirality and size effects of CNTs, in the following, an anisotropic elastic shell model will be used to obtain the natural frequencies of DWCNTs. Moreover, since He’s formulation is able to simulate more properly the actual anisotropic behaviour of CNTs, then in the following, it will be used as a reference to investigate the applicability and limitations of Ru’s formulation for vibration modelling of DWCNTs.

5.3. Applicability and Limitations of Ru’s Formulation for DWCNTs

In the present Section, starting from the previously reported results, natural frequencies of DWCNTs obtained via Ru’s and He’s formulations of the van der Waals interaction coefficient between the layers are compared. The Donnell shell theory and anisotropic elastic shell model are applied. Different values of diameters and wavenumbers along the longitudinal and circumferential directions are considered.
In Table 5, comparisons between the natural frequencies of the simply supported DWCNT of Table 1 with inner radius R 1 = 5   nm and aspect ratio L / R 2 = 10 are shown. Circumferential flexure modes ( n = 2 ) are analysed. The number of longitudinal half-waves 1 m 3 is studied. Vibration modes with a prevalent longitudinal u , circumferential v or radial w displacement component, which corresponds to a primarily axial, torsional or flexural vibration, respectively, are investigated.
From these comparisons, it can be observed that the percentage difference between the results of the two formulations is concentrated in the highest natural frequencies corresponding to the modal shape with prevalent radial displacement w , while elsewhere the percentage difference is negligible. The same result is found by comparing the natural frequencies of the present model in the isotropic form obtained via Ru’s and He’s formulations of the van der Waals interaction coefficient in Table 3 and Table 4.
Starting from these results, in the following, only the highest radial natural frequencies will be considered to investigate the differences between the two formulations (Ru’s vs. He’s) adopted for the simulation of the van der Waals interaction coefficient.
It should be underlined that, in a previous work of the same author, it was obtained that the differences in the results due to the different shell theory (Sanders vs. Donnell) considered to model the linear vibrations of DWCNTs are concentrated in the lowest radial natural frequencies, see Ref. [46] for more details.
Let us extend the previous analysis of Table 1 by considering the range of longitudinal half-waves 0 m 5 , the range of circumferential waves 0 n 8 and the range of inner radius values 0.5   nm R 1 10   nm . In addition, on the basis of the previous results, let us focus our attention on the highest radial natural frequencies.
In Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, the curves of the percentage differences between the natural frequencies of the simply supported DWCNT of Table 1 obtained by considering Ru’s and He’s formulations are shown, where the more accurate He’s formulation is assumed as the reference. Vibration modes with a different number of waves along the longitudinal and circumferential directions ( m , n ) are analysed. Different values of inner radius R 1 are adopted. The aim of the present study is to investigate both applicability and limitations of Ru’s formulation for vibration modelling of DWCNTs.
In Figure 4, a simply supported DWCNT with inner radius R 1 = 0.5   nm is considered. From this Figure, it is observed that the percentage difference assumes the initial value of 1.5 % for n = 0 . Then, it decreases to the local minimum of 0.4 % for n = 1 and increases to the local maximum of 14 % for n = 2 . Finally, for n 3 the percentage difference decreases very quickly with increasing n .
In Figure 5, a simply supported DWCNT with inner radius R 1 = 1   nm is analysed. From this Figure it is observed that the percentage difference assumes the initial value of 4 % for n = 0 . Then, it decreases to the local minimum of 1.5 % for n = 1 and increases to the local maximum of 10.5 % for n = 2 . Finally, for n 3 the percentage difference decreases quickly with increasing n .
In Figure 6, a simply supported DWCNT with inner radius R 1 = 2   nm is investigated. From this Figure, it is observed that the percentage difference assumes the initial value of 6 % for n = 0 . Then, it decreases to the local minimum of 1 % for n = 2 and increases to the local maximum of 8 % for n = 4 . Finally, for n 5 the percentage difference decreases very slowly with increasing n .
In Figure 7, a simply supported DWCNT with inner radius R 1 = 5   nm is examined. From this Figure, it is observed that the percentage difference assumes the initial value of 5.5 % for n = 0 , which is maintained quasi-constant for n = ( 1 , 2 ) . Then, it decreases to the local minimum of 2 % for n = 4 and increases to the local maximum of 6 % for n = 6 , which is maintained quasi-constant with increasing n .
In Figure 8, a simply supported DWCNT with inner radius R 1 = 10   nm is studied. From this Figure, it is observed that the percentage difference assumes the initial value of 5 % for n = 0 , which is maintained quasi-constant for 1 n 4 . Then, it decreases slowly for n = ( 5 , 6 ) and finally very quickly with increasing n .
From Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 it is noted that, for all considered values of inner radius R 1 , the behaviour is almost the same for every number of longitudinal half-waves m ; therefore, the percentage differences of Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 can be rightly synthesized in Figure 9, which shows the percentage differences between the highest radial natural frequencies obtained via Ru’s (38) and He’s (39) formulations for the inner radius range 0.5   nm R 1 10   nm and the number of longitudinal half-waves m = 1 .
From Figure 9 it can be observed that, as the inner radius R 1 increases, the value of the initial percentage difference for n = 0 increases, converging to the final value of 5 % for R 1 = 10   nm .
Again, as the inner radius R 1 increases, the oscillation amplitude of the portion of the percentage difference curve between the local minimum and maximum decreases, where the local minimum and maximum tend to approach, moving towards increasing numbers of circumferential waves n .
Finally, as the inner radius R 1 increases, the percentage difference curves tend to converge to the quasi-constant curve corresponding to R 1 = 10   nm , without maximum and minimum values and with the percentage difference close to 5%.
Since, in Section 5.2, it was derived that, for a sufficiently large inner radius, the interaction coefficient no longer depends on the radius, but it assumes a constant value, which is similar to Ru’s and He’s formulations, then, in the following, the limit value of applicability of Ru’s formulation for the vibration modelling of DWCNTs will be imposed at 5%, which is the value of the quasi-constant percentage difference curve corresponding to R 1 = 10   nm .
From Figure 9, the following comments for increasing values of the inner radius R 1 can be made:
  • for n = ( 0 , 1 ) the percentage difference is initially very low ( < 2 % ) , then it increases, finally stabilizing at 5 % ;
  • for n = 2 , the percentage difference is initially extremely high ( 15 % ) , then it decreases by assuming a minimum value of 2 % and finally increases, stabilizing at 5 % ;
  • for n = 3 , the percentage difference is initially very high ( 10 % ) , then it reduces, stabilizing at 5 % ;
  • for n = 4 , the percentage difference is initially relatively low ( 4 % ) , then it increases until a maximum value of 10 % , it decreases until a minimum value of 2 % and finally it increases, stabilizing at 5 % ;
  • for 5 n 8 , the percentage difference is initially very low ( < 2 % ) , then it increases until a maximum value of 8 % and finally decreases at a value lower than 5 % .
Therefore, by assuming the more accurate He’s formulation as a reference, considering the Donnell shell theory and adopting an anisotropic elastic shell model, from Figure 9, with respect to the number of circumferential waves n and as a function of the value of the inner radius R 1 , the following ranges of applicability of Ru’s formulation for vibration modelling of DWCNTs are identified (percentage difference 5 % ):
  • n = 0 (axisymmetric modes), for every value of R 1 ;
  • n = 1 (beam-like modes), for every value of R 1 ;
  • n = 2 (circumferential flexure modes), for R 1 2   nm ;
  • n = 3 (shell-like modes), for R 1 5   nm ;
  • n = ( 4 , 5 ) (shell-like modes), for R 1 0.5   nm and for R 1 5   nm ;
  • n = 6 (shell-like modes), for R 1 0.5   nm and for R 1 10   nm ;
  • n = ( 7 , 8 ) (shell-like modes), for R 1 1   nm and for R 1 10   nm ;
while the corresponding ranges of limitation of Ru’s formulation for vibration modelling of DWCNTs are (percentage difference > 5 % ):
  • n = 2 (circumferential flexure modes), for R 1 1   nm ;
  • n = 3 (shell-like modes), for R 1 2   nm ;
  • n = ( 4 , 5 ) (shell-like modes), for 1   nm R 1 2   nm ;
  • n = 6 (shell-like modes), for 1   nm R 1 5   nm ;
  • n = ( 7 , 8 ) (shell-like modes), for 2   nm R 1 5   nm .
Similarly, from Figure 9, with respect to the value of inner radius R 1 and as a function of the number of circumferential waves n , the following ranges of applicability of Ru’s formulation for the vibration modelling of DWCNTs are identified (percentage difference 5 % ):
  • R 1 = 0.5   nm (very low inner radius), for n 1 and for n 4 ;
  • R 1 = 1   nm (low inner radius), for n 1 and for n 7 ;
  • R 1 = 2   nm (intermediate inner radius), for n 2 ;
  • R 1 = 5   nm (high inner radius), for n 5 ;
  • R 1 = 10   nm (very high inner radius), for every number of n ;
while the corresponding ranges of limitation of Ru’s formulation for vibration modelling of DWCNTs are (percentage difference > 5 % ):
  • R 1 = 0.5   nm (very low inner radius), for n = ( 2 , 3 ) ;
  • R 1 = 1   nm (low inner radius), for 2 n 6 ;
  • R 1 = 2   nm (intermediate inner radius), for n 3 ;
  • R 1 = 5   nm (high inner radius), for n 6 .
As previously verified, the vibrational mode shapes, and the corresponding natural frequencies, of DWCNTs are sensitive to the van der Waals interactions. Actually, from Figure 3 it can be noted that the maximum difference between the values of Ru’s and He’s van der Waals interaction coefficient is present at the inner radius R1 = 0.5 nm. So, the investigation of the critical order number (m, n) of the mode shape with which the difference between Ru’s and He’s formulations cannot be ignored should be carried out at that value of the inner radius, see Figure 4. First of all, from Figure 4 it can be observed that the behaviour is almost the same for every number of longitudinal half-waves m, and therefore, the analysis will be focused on the critical number of circumferential waves n. By imposing the limit value of percentage difference of 5%, below which Ru’s formulation can be considered as sufficiently accurate, it is found that Ru’s formulation cannot be applied in the range of the number of circumferential waves n = (2,3), which denotes the critical order number of the mode shapes.
In the end, we want to validate the previous assumption of Section 5.2, i.e., that it is necessary to assume an anisotropic model for R 1 < 1   nm (relatively small inner radius), while it can be adopted an isotropic model for R 1 > 5   nm (large inner radius).
In Figure 10, the percentage difference in curves between the highest radial natural frequencies obtained via He’s formulation by adopting isotropic and anisotropic elastic shell models (anisotropic model as reference) are reported. He’s formulation is applied since it takes into account the effect of the radius, and therefore, is able to simulate the actual anisotropic behaviour of CNTs more correctly than Ru’s one. A simply supported DWCNT with an aspect ratio L / R 2 = 10 is analysed. The Donnell shell theory is adopted. Circumferential flexure modes ( n = 2 ) are studied. Different values of inner radius R 1 and numbers of longitudinal half-waves m are investigated.
From this Figure, it is noted that, for a relatively small inner radius R 1 < 1   nm , the percentage difference is significant ( 3.5 % ) , and therefore, the anisotropic model must be adopted; conversely, for a large inner radius R 1 > 5   nm , the percentage difference is negligible ( < 0.2 % ) , and therefore, within this range, the isotropic model can also be correctly adopted. Again, this behaviour is almost the same for every number of longitudinal half-waves m .

6. Conclusions

In the present paper, the comparison is conducted between Ru’s and He’s formulations of the van der Waals interaction coefficient for the evaluation of the natural frequencies of DWCNTs. The actual discrete DWCNT is modelled by means of a couple of concentric equivalent continuous cylindrical shells, where the Donnell shell theory is considered to obtain strain-displacement relationships and an anisotropic elastic shell model is adopted to take into account the chirality effects. Simply supported boundary conditions are imposed and Rayleigh–Ritz method is used to obtain approximate natural frequencies and mode shapes. Increasing values of the inner ratio R 1 are considered. Different numbers of longitudinal half-waves m and circumferential waves n are investigated.
The main findings are reported below:
  • for a small inner radius, the interaction coefficient strongly depends on the radius, and therefore it assumes very different values for Ru’s and He’s formulations;
  • for a large inner radius, the interaction coefficient no longer depends on the radius, and it assumes a constant value, which is similar for Ru’s and He’s formulations;
  • He’s formulation is able to simulate the actual anisotropic behaviour of CNTs more accurately than Ru’s formulation, and is assumed as reference for the comparisons;
  • the difference in the results between Ru’s and He’s formulations is localised in the highest radial natural frequencies;
  • for axisymmetric and beam-like modes, Ru’s formulation can be applied for every value of the inner radius;
  • for shell-like modes with a relative low number of circumferential waves, Ru’s formulation cannot be applied for a small inner radius;
  • for shell-like modes with a high number of circumferential waves, Ru’s formulation can be applied for every value of the inner radius;
  • for a large inner radius, Ru’s formulation can be applied for every number of circumferential waves;
  • the behaviour is almost the same for every number of longitudinal half-waves.
In all papers available in the literature that are focused on the analysis of MWCNT vibrations, the van der Waals interaction coefficient is always modelled by considering either Ru’s or He’s formulation. The main novelty of the present work is that, for the first time, the results of Ru’s and He’s formulations are compared with each other. From this comparison, which is carried out by varying DWCNT geometry and wavenumber, it is derived when it is possible to apply the less accurate but simpler Ru’s formulation and when it is necessary to use the more accurate but also more complex He’s formulation.

Funding

This research was funded by Department of Sciences and Methods for Engineering (University of Modena and Reggio Emilia) grant number 020142_22_FRN_SOSTEGNO_RICERCA_DISMI.

Acknowledgments

The Author is grateful to the Department of Sciences and Methods for Engineering of the University of Modena and Reggio Emilia (Grant 020142_22_FRN_SOSTEGNO_RICERCA_DISMI) for the financial support of the present work, which is dedicated to the blessed memory of Professor Leonid I. Manevitch, outstanding educator, respected teacher, admired scholar.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Iijima, S. Helical microtubules of graphitic carbon. Nature 1991, 354, 56–58. [Google Scholar] [CrossRef]
  2. Ibrahim, K.S. Carbon nanotubes, properties and applications: A review. Carbon Lett. 2013, 14, 131–144. [Google Scholar] [CrossRef] [Green Version]
  3. Khan, D.; Ali, Z.; Asif, D.; Panjwani, M.K.; Khan, I. Incorporation of carbon nanotubes in photoactive layer of organic solar cells. Ain Shams Eng. J. 2021, 12, 897–900. [Google Scholar] [CrossRef]
  4. Luo, C.; Xie, H.; Wang, Q.; Luo, G.; Liu, C. A Review of the Application and Performance of Carbon Nanotubes in Fuel Cells. J. Nanomater. 2015, 2015, 4. [Google Scholar] [CrossRef] [Green Version]
  5. Rao, A.M.; Richter, E.; Bandow, S.; Chase, B.; Eklund, P.C.; Williams, K.A.; Fang, S.; Subbaswamy, K.R.; Menon, M.; Thess, A.; et al. Diameter-Selective Raman Scattering from Vibrational Modes in Carbon Nanotubes. Science 1997, 275, 187–191. [Google Scholar] [CrossRef] [Green Version]
  6. Bandow, S.; Asaka, S.; Saito, Y.; Rao, A.M.; Grigorian, L.; Richter, E.; Eklund, P.C. Effect of the Growth Temperature on the Diameter Distribution and Chirality of Single-Wall Carbon Nanotubes. Phys. Rev. Lett. 1998, 80, 3779–3782. [Google Scholar] [CrossRef]
  7. Dresselhaus, G.; Dresselhaus, M.S.; Hafner, J.H.; Hunter, M.; Jorio, A.; Lieber, C.M.; McClure, T.; Saito, R. Structural (n, m) Determination of Isolated Single-Wall Carbon Nanotubes by Resonant Raman Scattering. Phys. Rev. Lett. 2001, 86, 1118. [Google Scholar] [CrossRef]
  8. Gupta, S.S.; Bosco, F.G.; Batra, R.C. Breakdown of structural models for vibrations of single-wall zigzag carbon nanotubes. J. Appl. Phys. 2009, 106, 063527. [Google Scholar] [CrossRef] [Green Version]
  9. Cheng, H.C.; Liu, Y.L.; Wu, C.; Chen, W.H. On radial breathing vibration of carbon nanotubes. Comput. Methods Appl. Mech. Eng. 2010, 199, 2820–2827. [Google Scholar] [CrossRef]
  10. Duan, W.H.; Wang, C.; Zhang, Y. Calibration of nonlocal scaling effect parameter for free vibration of carbon nanotubes by molecular dynamics. J. Appl. Phys. 2007, 101, 024305. [Google Scholar] [CrossRef]
  11. Odegard, G.M.; Gates, T.S.; Nicholson, L.M.; Wise, K.E. Equivalent-Continuum Modeling of Nano-Structured Materials. Compos. Sci. Technol. 2002, 62, 1869–1880. [Google Scholar] [CrossRef] [Green Version]
  12. Arroyo, M.; Belytschko, T. Continuum Mechanics Modeling and Simulation of Carbon Nanotubes. Meccanica 2005, 40, 455–469. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, P.; Huang, Y.; Geubelle, P.H.; Hwang, K. On the continuum modeling of carbon nanotubes. Acta Mech. Sin. 2002, 18, 528–536. [Google Scholar] [CrossRef]
  14. Yakobson, B.I.; Brabec, C.J.; Bernholc, J. Nanomechanics of Carbon Tubes: Instabilities beyond Linear Response. Phys. Rev. Lett. 1996, 76, 2511–2514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wang, C.; Ru, C.Q.; Mioduchowski, A. Applicability and Limitations of Simplified Elastic Shell Equations for Carbon Nanotubes. J. Appl. Mech. 2004, 71, 622–631. [Google Scholar] [CrossRef]
  16. Silvestre, N.; Wang, C.; Zhang, Y.; Xiang, Y. Sanders shell model for buckling of single-walled carbon nanotubes with small aspect ratio. Compos. Struct. 2011, 93, 1683–1691. [Google Scholar] [CrossRef]
  17. Silvestre, N. On the accuracy of shell models for torsional buckling of carbon nanotubes. Eur. J. Mech.-A/Solids 2012, 32, 103–108. [Google Scholar] [CrossRef]
  18. Strozzi, M.; Manevitch, L.I.; Pellicano, F.; Smirnov, V.V.; Shepelev, D.S. Low-frequency linear vibrations of single-walled carbon nanotubes: Analytical and numerical models. J. Sound Vib. 2014, 333, 2936–2957. [Google Scholar] [CrossRef] [Green Version]
  19. Strozzi, M.; Smirnov, V.V.; Manevitch, L.I.; Pellicano, F. Nonlinear vibrations and energy exchange of single-walled carbon nanotubes. Radial breathing modes. Compos. Struct. 2018, 184, 613–632. [Google Scholar] [CrossRef] [Green Version]
  20. Strozzi, M.; Pellicano, F. Nonlinear Resonance Interaction between Conjugate Circumferential Flexural Modes in Single-Walled Carbon Nanotubes. Shock Vib. 2019, 2019, 3241698. [Google Scholar] [CrossRef]
  21. Strozzi, M.; Smirnov, V.V.; Manevitch, L.I.; Pellicano, F. Nonlinear normal modes, resonances and energy exchange in single-walled carbon nanotubes. Int. J. Non-Linear Mech. 2020, 120, 103398. [Google Scholar] [CrossRef]
  22. Chang, T.; Geng, J.; Guo, X. Prediction of chirality- and size-dependent elastic properties of single-walled carbon nanotubes via a molecular mechanics model. Proc. R. Soc. A Math. Phys. Eng. Sci. 2006, 462, 2523–2540. [Google Scholar] [CrossRef]
  23. Chang, T. A molecular based anisotropic shell model for single-walled carbon nanotubes. J. Mech. Phys. Solids 2010, 58, 1422–1433. [Google Scholar] [CrossRef]
  24. Ghavanloo, E.; Fazelzadeh, S. Vibration characteristics of single-walled carbon nanotubes based on an anisotropic elastic shell model including chirality effect. Appl. Math. Model. 2012, 36, 4988–5000. [Google Scholar] [CrossRef]
  25. Fazelzadeh, S.A.; Ghavanloo, E. Nonlocal anisotropic elastic shell model for vibrations of single-walled carbon nanotubes with arbitrary chirality. Compos. Struct. 2012, 94, 1016–1022. [Google Scholar] [CrossRef]
  26. Leissa, A.W. Vibration of Shells; Acoustical Society of America: Columbus, OH, USA, 1993. [Google Scholar]
  27. Yamaki, N.; Simitses, G.J. Elastic Stability of Circular Cylindrical Shells. J. Appl. Mech. 1985, 52, 501–502. [Google Scholar] [CrossRef]
  28. Soedel, W.; Qatu, M.S. Vibrations of Shells and Plates, Third Edition. J. Acoust. Soc. Am. 2005, 117, 1683. [Google Scholar] [CrossRef]
  29. Ventsel, E.; Krauthammer, T.; Carrera, E. Thin Plates and Shells: Theory, Analysis, and Applications. Appl. Mech. Rev. 2002, 55, B72–B73. [Google Scholar] [CrossRef]
  30. Jones, J.E. On the determination of molecular fields.—I. From the variation of the viscosity of a gas with temperature. Proc. R. Soc. London. Ser. A Math. Phys. Sci. 1924, 106, 441–462. [Google Scholar] [CrossRef] [Green Version]
  31. Ru, C.Q. Elastic Models for Carbon Nanotubes. Encycl. Nanosci. Nanotechnol. 2004, 2, 731–744. [Google Scholar]
  32. Ru, C.Q. Axially compressed buckling of a double-walled carbon nanotube embedded in an elastic medium. J. Mech. Phys. Solids 2001, 49, 1265–1279. [Google Scholar] [CrossRef]
  33. Ru, C.Q. Column buckling of multiwalled carbon nanotubes with interlayer radial displacements. Phys. Rev. B 2000, 62, 16962–16967. [Google Scholar] [CrossRef]
  34. Girifalco, L.A.; Lad, R.A. Energy of Cohesion, Compressibility, and the Potential Energy Functions of the Graphite System. J. Chem. Phys. 1956, 25, 693–697. [Google Scholar] [CrossRef]
  35. Pentaras, D.; Elishakoff, I. Effective Approximations for Natural Frequencies of Double-Walled Carbon Nanotubes Based on Donnell Shell Theory. J. Nanotechnol. Eng. Med. 2011, 2, 021013. [Google Scholar] [CrossRef]
  36. Elishakoff, I.; Pentaras, D. Fundamental natural frequencies of double-walled carbon nanotubes. J. Sound Vib. 2009, 322, 652–664. [Google Scholar] [CrossRef]
  37. Pentaras, D.; Elishakoff, I. Free vibration of triple-walled carbon nanotubes. Acta Mech. 2011, 221, 239–249. [Google Scholar] [CrossRef]
  38. Ehteshami, H.; Hajabasi, M.A. Analytical approaches for vibration analysis of multi-walled carbon nanotubes modeled as multiple nonlocal Euler beams. Phys. E Low-Dimens. Syst. Nanostructures 2011, 44, 270–285. [Google Scholar] [CrossRef]
  39. Xu, K.Y.; Guo, X.N.; Ru, C.Q. Vibration of a double-walled carbon nanotube aroused by nonlinear intertube van der Waals forces. J. Appl. Phys. 2006, 99, 064303. [Google Scholar] [CrossRef]
  40. Xu, K.-Y.; Aifantis, E.C.; Yan, Y.-H. Vibrations of Double-Walled Carbon Nanotubes With Different Boundary Conditions Between Inner and Outer Tubes. J. Appl. Mech. 2008, 75, 021013. [Google Scholar] [CrossRef]
  41. He, X.; Kitipornchai, S.; Liew, K. Buckling analysis of multi-walled carbon nanotubes: A continuum model accounting for van der Waals interaction. J. Mech. Phys. Solids 2005, 53, 303–326. [Google Scholar] [CrossRef]
  42. He, X.Q.; Eisenberger, M.; Liew, K.M. The effect of van der Waals interaction modeling on the vibration characteristics of multiwalled carbon nanotubes. J. Appl. Phys. 2006, 100, 124317. [Google Scholar] [CrossRef]
  43. He, X.Q.; Kitipornchai, S.; Wang, C.; Liew, K. Modeling of van der Waals force for infinitesimal deformation of multi-walled carbon nanotubes treated as cylindrical shells. Int. J. Solids Struct. 2005, 42, 6032–6047. [Google Scholar] [CrossRef] [Green Version]
  44. Strozzi, M.; Pellicano, F. Linear vibrations of triple-walled carbon nanotubes. Math. Mech. Solids 2018, 23, 1456–1481. [Google Scholar] [CrossRef]
  45. Strozzi, M.; Gendelman, O.; Elishakoff, I.; Pellicano, F. Applicability and Limitations of Simplified Elastic Shell Theories for Vibration Modelling of Double-Walled Carbon Nanotubes. C 2021, 7, 61. [Google Scholar] [CrossRef]
  46. Strozzi, M.; Elishakoff, I.E.; Manevitch, L.I.; Gendelman, O.V. Applicability and limitations of Donnell shell theory for vibration modelling of double-walled carbon nanotubes. Thin-Walled Struct. 2022, 178, 109532. [Google Scholar] [CrossRef]
  47. Strozzi, M.; Smirnov, V.V.; Pellicano, F.; Kovaleva, M. Nonlocal anisotropic elastic shell model for vibrations of double-walled carbon nanotubes under nonlinear van der Waals interaction forces. Int. J. Non-Linear Mech. 2022, 146, 104172. [Google Scholar] [CrossRef]
  48. Natsuki, T.; Ni, Q.-Q.; Endo, M. Analysis of the vibration characteristics of double-walled carbon nanotubes. Carbon 2008, 46, 1570–1573. [Google Scholar] [CrossRef] [Green Version]
  49. Natsuki, T.; Lei, X.-W.; Ni, Q.-Q.; Endo, M. Free vibration characteristics of double-walled carbon nanotubes embedded in an elastic medium. Phys. Lett. A 2010, 374, 2670–2674. [Google Scholar] [CrossRef]
  50. Kumar, B.R. Investigation on mechanical vibration of double-walled carbon nanotubes with inter-tube van der Waals forces. Adv. Nano Res. 2018, 6, 135–145. [Google Scholar]
  51. Afsharirad, F.; Mousanezhad, S.; Biglari, H.; Rahmani, O. Molecular dynamics of axial interwall van der Waals force and mechanical vibration of double-walled carbon nanotubes. Mater. Today Commun. 2021, 28, 102708. [Google Scholar] [CrossRef]
  52. Asghar, S.; Naeem, M.N.; Hussain, M. Non-local effect on the vibration analysis of double walled carbon nanotubes based on Donnell shell theory. Phys. E Low-Dimens. Syst. Nanostructures 2020, 116, 113726. [Google Scholar] [CrossRef]
  53. Gul, U.; Aydogdu, M. Noncoaxial vibration and buckling analysis of embedded double-walled carbon nanotubes by using doublet mechanics. Compos. Part B: Eng. 2018, 137, 60–73. [Google Scholar] [CrossRef]
  54. Senthilkumar, V. Axial vibration of double-walled carbon nanotubes using double-nanorod model with van der Waals force under Pasternak medium and magnetic effects. Vietnam J. Mech. 2022, 44, 29–43. [Google Scholar] [CrossRef]
  55. Natsuki, J.; Wu, P.; Jiang, H.; Natsuki, T. Dynamic analysis of double-walled carbon nanotubes embedded in elastic medium under a nanoparticle delivery. Diam. Relat. Mater. 2022, 128, 109194. [Google Scholar] [CrossRef]
Figure 1. Continuum modelling of a DWCNT. (a) Actual discrete DWCNT; (b) geometry of the equivalent continuous thin circular cylindrical shell; (c) cross-section of the surface of the equivalent continuous shell.
Figure 1. Continuum modelling of a DWCNT. (a) Actual discrete DWCNT; (b) geometry of the equivalent continuous thin circular cylindrical shell; (c) cross-section of the surface of the equivalent continuous shell.
Carbon 08 00059 g001
Figure 2. Schematic illustration of a ( r , s ) single-walled carbon nanotube. (a) Global structure with a zoom of a single representative atom; (b) side view of the local structure; (c) top view of the local structure. Reprinted with permission from Ref. [22], 2006, The Royal Society.
Figure 2. Schematic illustration of a ( r , s ) single-walled carbon nanotube. (a) Global structure with a zoom of a single representative atom; (b) side view of the local structure; (c) top view of the local structure. Reprinted with permission from Ref. [22], 2006, The Royal Society.
Carbon 08 00059 g002
Figure 3. Effect of the inner radius R1 (10−9 m) on the value of van der Waals interaction coefficient c (1020 N/m3) of the DWCNT of Table 1. Comparison between Ru’s (38) and He’s (39) formulations.
Figure 3. Effect of the inner radius R1 (10−9 m) on the value of van der Waals interaction coefficient c (1020 N/m3) of the DWCNT of Table 1. Comparison between Ru’s (38) and He’s (39) formulations.
Carbon 08 00059 g003
Figure 4. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 0.5 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Figure 4. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 0.5 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Carbon 08 00059 g004
Figure 5. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 1 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Figure 5. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 1 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Carbon 08 00059 g005
Figure 6. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 2 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Figure 6. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 2 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Carbon 08 00059 g006
Figure 7. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 5 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Figure 7. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 5 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Carbon 08 00059 g007
Figure 8. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 10 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Figure 8. Percentage differences between the highest radial natural frequencies of the simply supported DWCNT of Table 1 with inner radius R1 = 10 nm and aspect ratio L/R2 = 10 from Ru’s (38) and He’s (39) formulations. Donnell shell theory. Anisotropic elastic shell model. Number of longitudinal half-waves m. Number of circumferential waves n.
Carbon 08 00059 g008
Figure 9. Percentage differences between the highest radial natural frequencies obtained via Ru’s (38) and He’s (39) formulations (He as reference). Donnell shell theory. Anisotropic elastic model. Simply supported DWCNT of Table 1 with m = 1 longitudinal half-wave and aspect ratio L/R2 = 10. Inner radius R1. Number of circumferential waves n.
Figure 9. Percentage differences between the highest radial natural frequencies obtained via Ru’s (38) and He’s (39) formulations (He as reference). Donnell shell theory. Anisotropic elastic model. Simply supported DWCNT of Table 1 with m = 1 longitudinal half-wave and aspect ratio L/R2 = 10. Inner radius R1. Number of circumferential waves n.
Carbon 08 00059 g009
Figure 10. Percentage differences between the highest radial natural frequencies obtained via isotropic and anisotropic elastic shell models (anisotropic model as reference). Donnell shell theory. He’s formulation (39). Simply supported DWCNT of Table 1 and Table 2 with aspect ratio L/R2 = 10 and n = 2. Inner radius R1. Number of longitudinal half-waves m.
Figure 10. Percentage differences between the highest radial natural frequencies obtained via isotropic and anisotropic elastic shell models (anisotropic model as reference). Donnell shell theory. He’s formulation (39). Simply supported DWCNT of Table 1 and Table 2 with aspect ratio L/R2 = 10 and n = 2. Inner radius R1. Number of longitudinal half-waves m.
Carbon 08 00059 g010
Table 2. Mechanical parameters adopted in the isotropic elastic shell model.
Table 2. Mechanical parameters adopted in the isotropic elastic shell model.
Tensile rigidity C   ( J / m 2 ) 360
Bending rigidity D   ( J ) 1.362 × 10 19
Surface density ρ s   ( kg / m 2 ) 7.718 × 10 7
Poisson’s ratio ν 0.19
Table 3. Natural frequencies (THz) of a simply supported DWCNT with innermost radius R1 = 5 nm and aspect ratio L/R2 = 10. Donnell shell theory. Isotropic elastic shell model. Van der Waals interaction coefficient expressed on the basis of Ru’s formulation. Circumferential flexure modes (n = 2). Comparisons between present and Ref. [35] models.
Table 3. Natural frequencies (THz) of a simply supported DWCNT with innermost radius R1 = 5 nm and aspect ratio L/R2 = 10. Donnell shell theory. Isotropic elastic shell model. Van der Waals interaction coefficient expressed on the basis of Ru’s formulation. Circumferential flexure modes (n = 2). Comparisons between present and Ref. [35] models.
Natural Frequency (THz)Diff. (%)
Mode
(m,n)
Prevalent DisplacementPrevalent
Mode Shape
Present
Model (Ru)
Ref. [35]
Model (Ru)
(1,2)wR10.01630.01620.62
(2,2)wR10.05110.05090.39
(3,2)wR10.10250.10200.49
(1,2)uR20.85080.84380.83
(2,2)uR20.88500.87750.85
(1,2)uR10.90700.89990.79
(3,2)uR20.93340.92540.86
(2,2)uR10.94460.93750.76
(3,2)uR11.00500.99770.73
(1,2)vR21.31141.30800.26
(2,2)vR21.36011.35670.25
(3,2)vR21.43711.43340.26
(1,2)vR11.53841.53710.08
(2,2)vR11.56921.56790.08
(3,2)vR11.62141.62010.08
(1,2)wR12.70062.70000.02
(2,2)wR12.70112.70050.02
(3,2)wR12.70212.70150.02
Table 4. Natural frequencies (THz) of a simply supported DWCNT with innermost radius R1 = 5 nm and aspect ratio L/R2 = 10. Donnell shell theory. Isotropic elastic shell model. Van der Waals interaction coefficient expressed on the basis of He’s formulation. Circumferential flexure modes (n = 2). Comparisons between present and Ref. [42] models.
Table 4. Natural frequencies (THz) of a simply supported DWCNT with innermost radius R1 = 5 nm and aspect ratio L/R2 = 10. Donnell shell theory. Isotropic elastic shell model. Van der Waals interaction coefficient expressed on the basis of He’s formulation. Circumferential flexure modes (n = 2). Comparisons between present and Ref. [42] models.
Natural Frequency (THz)Diff. (%)
Mode
(m,n)
Prevalent DisplacementPrevalent
Mode Shape
Present
Model (He)
Ref. [42]
Model (He)
(1,2)wR10.01630.01630.00
(2,2)wR10.05110.05090.39
(3,2)wR10.10250.10210.39
(1,2)uR20.85080.84400.81
(2,2)uR20.88520.87820.80
(1,2)uR10.90710.89980.81
(3,2)uR20.93390.92660.79
(2,2)uR10.94470.93720.80
(3,2)uR11.00520.99720.80
(1,2)vR21.31571.31300.21
(2,2)vR21.36421.36200.16
(3,2)vR21.44091.43800.20
(1,2)vR11.53871.53600.18
(2,2)vR11.56961.56700.17
(3,2)vR11.62181.61900.17
(1,2)wR12.78312.78300.00
(2,2)wR12.78352.78400.02
(3,2)wR12.78432.78400.01
Table 5. Natural frequencies (THz) of the simply supported DWCNT of Table 1 with inner radius R1 = 5 nm and aspect ratio L/R2 = 10. Donnell shell theory. Anisotropic elastic shell model. Circumferential flexure modes (n = 2). Van der Waals interaction forces. Comparisons between Ru’s (38) and He’s (39) formulations.
Table 5. Natural frequencies (THz) of the simply supported DWCNT of Table 1 with inner radius R1 = 5 nm and aspect ratio L/R2 = 10. Donnell shell theory. Anisotropic elastic shell model. Circumferential flexure modes (n = 2). Van der Waals interaction forces. Comparisons between Ru’s (38) and He’s (39) formulations.
Natural Frequency (THz)Diff. (%)
Mode
(m,n)
Prevalent DisplacementPrevalent
Mode Shape
Anisotropic
Model (Ru)
Anisotropic
Model (He)
(1,2)wR10.01620.01620.00
(2,2)wR10.05100.05100.00
(3,2)wR10.10230.10230.00
(1,2)uR20.85930.85940.01
(2,2)uR20.89360.89430.08
(1,2)uR10.91600.91610.01
(3,2)uR20.94210.94370.17
(2,2)uR10.95380.95410.03
(3,2)uR11.01491.01520.03
(1,2)vR21.29751.30860.85
(2,2)vR21.34621.35670.77
(3,2)vR21.42301.43270.68
(1,2)vR11.52481.52570.06
(2,2)vR11.55501.55600.06
(3,2)vR11.60631.60740.07
(1,2)wR12.63132.85077.70
(2,2)wR12.63192.85107.69
(3,2)wR12.63292.85177.67
Table 1. Mechanical parameters adopted in the anisotropic elastic shell model.
Table 1. Mechanical parameters adopted in the anisotropic elastic shell model.
C-C Bond Parameters
C-C bond length a   ( nm ) 0.142
C-C bond elongation K ρ   ( nN / nm ) 742
C-C bond angle variance K θ   ( nN · nm ) 1.42
C-C distance parameters
C-C potential depth ε   ( 10 22   J ) 4.755
C-C equilibrium separation distance σ   ( nm ) 0.3407
CNT equivalent parameters
Thickness h   ( nm ) 0.066
Mass density ρ m   ( kg / m 3 ) 11,700
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Strozzi, M. Applicability and Limitations of Ru’s Formulation for Vibration Modelling of Double-Walled Carbon Nanotubes. C 2022, 8, 59. https://doi.org/10.3390/c8040059

AMA Style

Strozzi M. Applicability and Limitations of Ru’s Formulation for Vibration Modelling of Double-Walled Carbon Nanotubes. C. 2022; 8(4):59. https://doi.org/10.3390/c8040059

Chicago/Turabian Style

Strozzi, Matteo. 2022. "Applicability and Limitations of Ru’s Formulation for Vibration Modelling of Double-Walled Carbon Nanotubes" C 8, no. 4: 59. https://doi.org/10.3390/c8040059

APA Style

Strozzi, M. (2022). Applicability and Limitations of Ru’s Formulation for Vibration Modelling of Double-Walled Carbon Nanotubes. C, 8(4), 59. https://doi.org/10.3390/c8040059

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop