Next Article in Journal
Nontrivial Topological Properties and Synthesis of Sn2CoS with L21 Structure
Next Article in Special Issue
Semi-Empirical Pseudopotential Method for Graphene and Graphene Nanoribbons
Previous Article in Journal
A Visible-Light-Enhanced Heterogeneous Photo Degradation of Tetracycline by a Nano-LaFeO3 Catalyst with the Assistance of Persulfate
Previous Article in Special Issue
Seed Priming with Single-Walled Carbon Nanotubes Grafted with Pluronic P85 Preserves the Functional and Structural Characteristics of Pea Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparison of Shell Theories for Vibration Analysis of Single-Walled Carbon Nanotubes Based on an Anisotropic Elastic Shell Model

1
Department of Sciences and Methods for Engineering, University of Modena and Reggio Emilia, 42122 Reggio Emilia, Italy
2
Department of Ocean and Mechanical Engineering, Florida Atlantic University, Boca Raton, FL 33431, USA
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(8), 1390; https://doi.org/10.3390/nano13081390
Submission received: 21 March 2023 / Revised: 11 April 2023 / Accepted: 11 April 2023 / Published: 17 April 2023
(This article belongs to the Special Issue Carbon Nanotubes and Nanosheets for Sustainable Solutions)

Abstract

:
In the present paper, a comparison is conducted between three classical shell theories as applied to the linear vibrations of single-walled carbon nanotubes (SWCNTs); specifically, the evaluation of the natural frequencies is conducted via Donnell, Sanders, and Flügge shell theories. The actual discrete SWCNT is modelled by means of a continuous homogeneous cylindrical shell considering equivalent thickness and surface density. In order to take into account the intrinsic chirality of carbon nanotubes (CNTs), a molecular based anisotropic elastic shell model is considered. Simply supported boundary conditions are imposed and a complex method is applied to solve the equations of motion and to obtain the natural frequencies. Comparisons with the results of molecular dynamics simulations available in literature are performed to check the accuracy of the three different shell theories, where the Flügge shell theory is found to be the most accurate. Then, a parametric analysis evaluating the effect of diameter, aspect ratio, and number of waves along the longitudinal and circumferential directions on the natural frequencies of SWCNTs is performed in the framework of the three different shell theories. Assuming the results of the Flügge shell theory as reference, it is obtained that the Donnell shell theory is not accurate for relatively low longitudinal and circumferential wavenumbers, for relatively low diameters, and for relatively high aspect ratios. On the other hand, it is found that the Sanders shell theory is very accurate for all the considered geometries and wavenumbers, and therefore, it can be correctly adopted instead of the more complex Flügge shell theory for the vibration modelling of SWCNTs.

1. Introduction

Since their discovery in 1991 in Japan in the laboratories of the NEC Corporation by Professor Sumio Iijima [1], the study of the vibrations of carbon nanotubes has represented a very demanding challenge for many researchers all over the world.
This interest in carbon nanotubes is due to their extraordinary mechanical properties, in particular the very high elastic modulus and tensile strength, together with their very small diameter, which allows them to reach natural frequencies of the THz order, and therefore leads them to be applied in several high sensitivity electro-mechanical systems, such as resonators, sensors, and oscillators [2,3,4,5,6].
In order to study the vibratory behaviour of carbon nanotubes, three different methods have been proposed: experimental analyses, molecular dynamics simulations, and continuous models.
The experimental analyses, conducted on the basis of resonant Raman spectroscopy, allow to obtain the natural frequencies only of peculiar modes of carbon nanotubes, the so-called “radial breathing modes”, which are characterized by zero waves both longitudinal and circumferential, as undeformed axisymmetric modes [7,8,9]. In addition to this, the main limitation of the experimental analyses is due to their very high technical difficulty and the need to work with scanning or transmission electron microscopes with very high resolution (and therefore very high cost).
Molecular dynamics simulations take into account the discrete nature of carbon nanotubes by modelling bonds and interaction forces between the different carbon atoms based on the fundamental concepts of molecular mechanics [10]. These analyses allow to obtain the natural frequencies of both radial breathing and beam-like modes, where the latter are very important for carbon nanotubes, which have high aspect ratio (i.e., the ratio between length and radius) and, therefore, have a vibratory behaviour very similar to beam structures [11,12]. The main limitation of molecular dynamics simulations results from the high computational effort needed to carry out numerical analyses, especially in the presence of a large number of carbon atoms, a condition that usually occurs in carbon nanotubes, due to the fact that they are long and of reduced thickness (i.e., the ratio between thickness and radius), which leads them to be considered “thin-walled structures”.
As an alternative to the experimental analyses and numerical simulations, researchers have proposed several continuous models, mainly of beam-type or shell-type, to study carbon nanotube vibrations via continuum mechanics [13].
In particular, the continuous beam-type models are able to properly simulate the torsional vibrations of nanotubes [14], whereas the continuous shell-type models are able to also properly simulate the flexural vibrations of nanotubes, which give the highest natural frequencies; they therefore represent more complete models [15].
In general, the main problem in the continuous modelling of carbon nanotubes is the choice of the equivalent parameters that allow to study the discrete nanotubes as continuous structures.
For the shell models, which as previously mentioned are the most accurate ones, Yakobson [16], starting from results of molecular dynamics simulations, proposed equivalent values of tensile and flexural stiffness, and therefore thickness and surface density, thus allowing the study of carbon nanotubes as continuous isotropic cylindrical shells.
Considering an isotropic elastic shell model and using these equivalent parameters, several papers have been published concerning the study of the natural frequencies of single- or multi-walled carbon nanotubes in the framework of thin shell theories; these principally include Donnell-Mushtari, Sanders-Koiter, and Lur’ye-Flügge-Byrne [17,18,19,20].
Readers interested in deepening the peculiarities of these thin shell theories can find the related strain–displacement relationships and equations of motion in the fundamental monographs by Leissa [21], Yamaki [22], Amabili [23], Soedel [24], and Ventsel [25].
An interesting comparison of thin shell theories for vibrations of circular cylindrical shells was carried out by Amabili [26]. It was derived from this work that, among the Donnell, Flügge, and Sanders shell theories, the Donnell shell theory, which has the lowest analytical complexity (i.e., the lowest number of terms in the equations of motion), presents the lowest accuracy in modelling shell vibrations. On the other hand, it was found that the Flügge shell theory, which has the highest analytical complexity, presents the highest accuracy. Finally, the Sanders shell theory, which has an intermediate analytical complexity, presents a relatively high accuracy. Moreover, an examination of cylindrical shell theories, specifically Donnell and Sanders shell theories, for the buckling of carbon nanotubes was performed by Wang et al. [27].
The same continuous models adopted to predict CNT vibrations at the nanoscopic scale can be used at the macroscopic scale to investigate the vibrations of different configurations of circular cylindrical shells, e.g., FGM shells [28,29,30], multi-layer composite shells [31,32], sandwich composite shells [33], and laminated orthotropic shells [34], which are considered in several mechanical components. It is important to stress that, in the passage from the nanoscopic to macroscopic scale, the continuous models lose the anisotropic behaviour due to the inherent chirality of carbon nanotubes.
However, given the intrinsically anisotropic nature of carbon nanotubes, in order to correctly study their vibratory behaviour and also to take into account the dependence of their elastic properties on chirality, it is necessary to adopt an anisotropic model.
To this aim, Chang [35,36], starting from considerations of molecular mechanics, developed a novel and very accurate anisotropic elastic shell model capable of correctly predicting the dependence of the material elastic properties on CNT chirality and dimensions, and therefore, also able to calculate the natural frequencies very similar to those obtained via molecular dynamics simulations.
Adopting this anisotropic elastic model, the linear vibrations of single-walled and multi-walled carbon nanotubes for different geometries and wave numbers were investigated separately via the Donnell [37], Sanders [38], and Flügge [39,40] thin shell theories.
To the knowledge of the authors of this paper, a study on the linear vibrations of SWCNTs based on an anisotropic elastic shell model that compares the values obtained for the natural frequencies by applying the Donnell, Sanders, and Flügge thin shell theories for different geometries and wavenumbers has not yet been published in the literature. In fact, the authors of this paper believe that it could be very useful to investigate the field of applicability and limitation of the three previously indicated thin shell theories in order to identify which of them is able to provide sufficiently accurate results with a relatively low computational effort.
To this aim, in the present paper, the natural frequencies of SWCNTs are obtained in the framework of the Donnell, Sanders, and Flügge shell theories, where the actual discrete SWCNT is modelled via a continuous homogeneous cylindrical shell considering equivalent thickness and surface density. An anisotropic elastic shell model is adopted to take into account the intrinsic chirality effects of CNTs. Simply supported boundary conditions are imposed and the complex method is used to solve the dynamic equations of motion and to obtain the natural frequencies. Vibration modes with different numbers of longitudinal and circumferential waves are studied. SWCNTs with different geometries are analysed.
Taking the more accurate but more complex Flügge theory as a reference, the main objective of this work was to establish whether the simpler Donnell or Sanders shell theories allow to obtain sufficiently accurate natural frequencies and, therefore, can be adopted instead of the Flügge theory to correctly model the linear vibrations of SWCNTS on the basis of an anisotropic elastic shell model.

2. Thin Shell Theories for SWCNTs

In the present paper, the actual discrete SWCNT of Figure 1a is modelled by means of an equivalent continuous elastic thin cylindrical shell, see Figure 1b,c with radius R , length L , and thickness h . A cylindrical coordinate system ( O , x , θ , z ) is adopted, where the origin O of the reference system is located at the centre of one end of the cylindrical 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 as 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 , and t is the time.
In this paper, the natural frequencies obtained by considering three different thin shell theories are compared regarding the linear vibrations of SWCNTs with different geometries and wavenumbers. These theories are based on Love’s first approximation assumptions [21]: (i) the thickness h of the shell is small with respect to the radius of curvature R of the middle surface; (ii) the strains are small; (iii) the transverse normal stress is small; (iv) and the normal to the undeformed middle surface remains straight and normal to the middle surface after the deformation, and undergoes no thickness stretching (Kirchhoff-Love kinematic hypothesis). The considered thin shell theories are: (a) Donnell-Mushtari [23], (b) Sanders-Koiter [24], and (c) Flügge-Lur’ye-Byrne [25]; for all of them, both rotary inertia and shear deformations are neglected.

2.1. Strain–Displacement Relationships

According to the Donnell, Sanders, and Flügge shell theories, the middle surface strains ( ε x , 0 , ε θ , 0 , γ x θ , 0 ) of the shell are related to the displacements ( u , v , w ) in the form [21]:
ε x , 0 = u x     ε θ , 0 = 1 R v θ + w R     γ x θ , 0 = v x + 1 R u θ
From Equation (1), it can be observed that the middle surface strains are expressed in the same form for the three different shell theories considered.
According to Donnell, Sanders, and Flügge, the middle surface changes in the curvature and torsion ( k x , k θ , k x θ ) of the shell are related to the displacements ( u , v , w ) in the form [21]:
k x = 2 w x 2 k θ = 1 R 2 2 w θ 2 + ψ v θ φ w k x θ = 2 R 2 w x θ 1 R 2 1 2 ψ + φ u θ + 1 R 3 2 ψ + φ v x
From Equation (2), it can be noted that the middle surface change in the curvature k θ and torsion k x θ are written in a different form for the three different shell theories, where parameters ( ψ = 0 , φ = 0 ) denote the Donnell shell theory, parameters ( ψ = 1 , φ = 0 ) denote the Sanders shell theory, and parameters ( ψ = 0 , φ = 1 ) denote the Flügge shell theory. Since the Sanders and Flügge shell theories have more terms in the Expansion (2) than the Donnell shell theory, the first two theories can be expected to be more accurate than the third in the modelling of SWCNT linear vibrations.

2.2. Strain Components at an Arbitrary Point of the Shell Surface

According to the Donnell, Sanders and Flügge shell theories, the strain components ( ε x , ε θ , γ x θ ) at an arbitrary point of the surface of the shell are related to the middle surface strains ( ε x , 0 , ε θ , 0 , γ x θ , 0 ) and to the changes in curvature and torsion of the middle surface ( k x , k θ , k x θ ) by the relationships [23]:
ε x = ε x , 0 + z k x     ε θ = ε θ , 0 + z k θ     γ x θ = γ x θ , 0 + z k x θ
where z is the distance of the considered arbitrary point of the shell from the middle surface.
Substituting Equations (1) and (2) into Equation (3), the following is obtained:
ε x = u x z 2 w x 2     ε θ = 1 R v θ + w R + z R 2 2 w θ 2 + ψ v θ φ w γ x θ = v x + 1 R u θ 2 z R 2 w x θ z R 2 1 2 ψ + φ u θ + z R 3 2 ψ + φ v x
where Equation (4) relate the strain components at an arbitrary point of the shell surface ( ε x , ε θ , γ x θ ) to the displacements ( u , v , w ) .

3. Anisotropic Elastic Shell Model

Considering the molecular based anisotropic elastic shell model developed by Chang [35,36], which includes the chirality effects characteristic of SWCNTs, the stress–strain relationships can be written as:
σ x = 1 h Y 11 ε x + Y 12 ε θ + Y 13 γ x θ σ θ = 1 h Y 21 ε x + Y 22 ε θ + Y 23 γ x θ τ x θ = 1 h Y 31 ε x + Y 32 ε θ + Y 33 γ x θ
where ( σ x , σ θ , τ x θ ) are the stress components at an arbitrary point of the shell surface, and Y i j are the anisotropic surface elastic constants of an arbitrary SWCNT, which are defined as [35,36]:
Y i j = 2 3 3 K ρ G l i G l j + 2 K θ a 2 H l i H l j , i , j , l = 1,2 , 3 ( sum over   l )
in which a is the carbon–carbon bond length; ( K ρ , K θ ) are force constants associated with stretching and angular distortion of the carbon-carbon bond, respectively, where these constants can be obtained from quantum (ab initio) mechanics, empirical molecular potential, or fitted to experimental data; and ( G l i , G l j , H l i , H l j ) are elements of matrices G and H , respectively, which are given in detail in Ref. [36].
σ x = 1 h Y 11 u x z 2 w x 2 + Y 12 1 R v θ + w R + z R 2 2 w θ 2 + ψ v θ φ w + Y 13 v x + 1 R u θ 2 z R 2 w x θ z R 2 1 2 ψ + φ u θ + z R 3 2 ψ + φ v x
σ θ = 1 h Y 21 u x z 2 w x 2 + Y 22 1 R v θ + w R + z R 2 2 w θ 2 + ψ v θ φ w + Y 23 v x + 1 R u θ 2 z R 2 w x θ z R 2 1 2 ψ + φ u θ + z R 3 2 ψ + φ v x
τ x θ = 1 h Y 31 u x z 2 w x 2 + Y 32 1 R v θ + w R + z R 2 2 w θ 2 + ψ v θ φ w + Y 33 v x + 1 R u θ 2 z R 2 w x θ z R 2 1 2 ψ + φ u θ + z R 3 2 ψ + φ v x

4. Equations of Motion

The general equations of motion for an arbitrary SWCNT in terms of force ( N x , N θ , N x θ ) and moment ( M x , M θ , M x θ ) resultants are written in the form [22]:
N x x + 1 R N x θ θ 1 2 R 2 M x θ θ ρ h 2 u t 2 = 0 1 R N θ θ + N x θ x + 3 2 R M x θ x + 1 R 2 M θ θ ρ h 2 v t 2 = 0 2 M x x 2 + 2 R 2 M x θ x θ + 1 R 2 2 M θ θ 2 N θ R ρ h 2 w t 2 = 0
where ρ h is the mass density per unit lateral area (i.e., the surface density) of SWCNT.
In the anisotropic elastic shell model, the force and moment resultants are defined based on the stress components in Equation (5), in the form [39]:
N x = h / 2 h / 2 σ x 1 + φ z R d z = Y 11 u x + Y 12 R v θ + w + Y 13 v x + 1 R u θ + φ X 11 R 2 w x 2 X 12 R 1 R 2 2 w θ 2 + w R 2 + X 13 R 2 R 2 w x θ 1 R 2 u θ + 1 R v x
N θ = h / 2 h / 2 σ θ d z = Y 21 u x + Y 22 R v θ + w + Y 23 v x + 1 R u θ
N x θ = h / 2 h / 2 τ x θ 1 + φ z R d z = Y 31 u x + Y 32 R v θ + w + Y 33 v x + 1 R u θ + φ X 31 R 2 w x 2 X 32 R 1 R 2 2 w θ 2 + w R 2 + X 33 R 2 R 2 w x θ 1 R 2 u θ + 1 R v x
M x = h / 2 h / 2 σ x 1 + φ z R z d z = X 11 2 w x 2 + X 12 R 2 2 w θ 2 + ψ v θ φ w + X 13 2 R 2 w x θ 1 R 2 1 2 ψ + φ u θ + 1 R 3 2 ψ + φ v x + φ X 11 R u x + X 12 R 2 v θ + w + X 13 R v x + 1 R u θ
M θ = h / 2 h / 2 σ θ z d z = X 21 2 w x 2 + X 22 R 2 2 w θ 2 + ψ v θ φ w + X 23 2 R 2 w x θ 1 R 2 1 2 ψ + φ u θ + 1 R 3 2 ψ + φ v x
M x θ = h / 2 h / 2 τ x θ 1 + φ z R z d z = X 31 2 w x 2 + X 32 R 2 2 w θ 2 + ψ v θ φ w + X 33 2 R 2 w x θ 1 R 2 1 2 ψ + φ u θ + 1 R 3 2 ψ + φ v x + φ X 31 R u x + X 32 R 2 v θ + w + X 33 R v x + 1 R u θ
where X i j = Y i j h 2 / 12 , with i , j = 1,2 , 3 .
From Equations (11)–(16) it is noted that the integrating functions of the circumferential force N θ and moment M θ resultants are the same for the three different shell theories. Conversely, the integrating functions of the other force and moment resultants are different, since in the Flügge shell theory, they also include the term ( 1 + z / R ) (in this theory the thinness assumption is delayed), whereas this term is ignored in the other two thin shell theories. On the one hand, considering this additional term within the integrating functions of the resultants certainly makes the Flügge shell theory more refined (in fact, it is able to correctly model the vibrations of even relatively thick shells). On the other hand, this considerably increases the number of terms within the expressions of such resultants, and therefore, it strongly increases the computational effort of the numerical analyses. The main goal of this paper is therefore to verify whether a less refined theory, but with lower computational effort, such as the Donnell or Sanders shell theory, can provide sufficiently accurate results in terms of SWCNT natural frequencies compared to those provided by the Flügge shell theory, as investigated in [26] for cylindrical shells.
By substituting the expressions of the force and moment resultants (11)–(16) into the dynamic Equation (10), the equations of motion for the anisotropic elastic shell model are obtained in the form:
Y 11 2 x 2 + 2 Y 13 R φ 3 X 13 2 R 3 2 x θ + Y 33 R 2 + ψ 4 φ X 33 4 R 4 2 θ 2 u + Y 13 + φ X 13 R 2 2 x 2 + Y 12 + Y 33 R ψ 3 X 33 4 R 3 2 x θ + Y 23 R 2 ψ + φ X 23 2 R 4 2 θ 2 v + Y 12 R φ X 12 R 3 x + Y 23 R 2 φ X 23 R 4 θ φ X 11 R 3 x 3 + 1 6 φ X 13 2 R 2 3 x 2 θ + 1 2 φ X 33 R 3 φ X 12 R 3 3 x θ 2 + 1 2 φ X 23 2 R 4 3 θ 3 w = ρ h 2 u t 2
Y 13 + φ 3 X 13 2 R 2 2 x 2 + Y 12 + Y 33 R 3 ψ + 4 φ X 33 4 R 3 2 x θ + Y 23 R 2 X 23 2 R 4 ψ + 2 φ 2 θ 2 u + Y 33 + 9 ψ + 16 φ X 33 4 R 2 2 x 2 + 2 Y 23 R + 6 ψ + 5 φ X 23 2 R 3 2 x θ + Y 22 R 2 + ψ X 22 R 4 2 θ 2 v + Y 23 R φ X 23 R 3 x + Y 22 R 2 φ X 22 R 4 θ ( 3 + 2 φ ) X 13 2 R 3 x 3 X 12 R 2 + 3 + 2 φ X 33 R 2 3 x 2 θ ( 7 + 2 φ ) X 23 2 R 3 3 x θ 2 X 22 R 4 3 θ 3 w = ρ h 2 v t 2
Y 12 R x Y 23 R 2 θ + φ X 11 R 3 x 3 + ψ + 4 φ X 13 2 R 2 3 x 2 θ ψ X 33 R 3 3 x θ 2 ψ + 2 φ X 23 2 R 4 3 θ 3 u + Y 23 R x Y 22 R 2 θ + 3 ψ + 4 φ X 13 2 R 3 x 3 + ψ + φ X 12 R 2 + 3 ψ + 4 φ X 33 R 2 3 x 2 θ + 7 ψ + 6 φ X 23 2 R 3 3 x θ 2 + ψ X 22 R 4 3 θ 3 v + Y 22 R 2 φ X 22 R 4 2 θ 2 X 11 4 x 4 4 X 13 R 4 x 3 θ 2 X 12 + 4 X 33 R 2 4 x 2 θ 2 4 X 23 R 3 4 x θ 3 X 22 R 4 4 θ 4 w = ρ h 2 w t 2

5. Solution Method

In this paper, a complex method is considered to analytically solve the dynamic equations of motions (17)–(19) and to obtain the natural frequencies of SWCNTs. Specifically, a complex variable is used to solve the partial differential Equations (17)–(19) by setting the real and imaginary zero.
In the present work, simply supported boundary conditions are adopted. These boundary conditions, for the complex method, impose the conditions R e v = R e w = 0 at both ends x = ( 0 , L ) of the SWCNT. The displacement field that satisfies these boundary conditions can be written as [39]:
u ( x , θ , t ) = U exp i λ q x cos s θ cos ω t v ( x , θ , t ) = i V exp i λ q x sin s θ cos ω t w ( x , θ , t ) = i W exp i λ q x cos s θ cos ω t
where U , V , W denote the displacement amplitudes along the longitudinal u , circumferential v , and radial w directions, respectively; i is the imaginary unit; λ q is the wavenumber along the longitudinal direction, with λ q = q π / L , where q is the number of longitudinal half-waves and L is the length of the SWCNT; s is the number of circumferential waves; and ω is the circular frequency.
Substituting Equation (20) into Equations (17)–(19), a set of algebraic equations for the displacement amplitudes U , V , W is obtained, which can be written in the form [39]:
E ( λ q , s , ω ) 3 × 3 U V W = 0 0 0
where E is a non-symmetric matrix, whose elements are:
E 11 = Y 11 λ q 2 + 2 Y 13 R φ 3 X 13 2 R 3 λ q s + Y 33 R 2 + ψ 4 φ X 33 4 R 4 s 2 ρ h ω 2
E 12 = Y 13 + φ X 13 R 2 λ q 2 Y 12 + Y 33 R ψ 3 X 33 4 R 3 λ q s Y 23 R 2 ψ + φ X 23 2 R 4 s 2
E 13 = Y 12 R φ X 12 R 3 λ q Y 23 R 2 φ X 23 R 4 s φ X 11 R λ q 3 + ( 1 6 φ ) X 13 2 R 2 λ q 2 s + ( 1 2 φ ) X 33 R 3 φ X 12 R 3 λ q s 2 + 1 2 φ X 23 2 R 4 s 3
E 21 = Y 13 + φ 3 X 13 2 R 2 λ q 2 Y 12 + Y 33 R 3 ψ + 4 φ X 33 4 R 3 λ q s Y 23 R 2 ψ + 2 φ X 23 2 R 4 s 2
E 22 = Y 33 + 9 ψ + 16 φ X 33 4 R 2 λ q 2 + 2 Y 23 R + 6 ψ + 5 φ X 23 2 R 3 λ q s + Y 22 R 2 + ψ X 22 R 4 s 2 ρ h ω 2
E 23 = Y 23 R φ X 23 R 3 λ q + Y 22 R 2 φ X 22 R 4 s + 3 + 2 φ X 13 2 R λ q 3 + X 12 R 2 + ( 3 + 2 φ ) X 33 R 2 λ q 2 s + 7 + 2 φ X 23 2 R 3 λ q s 2 + X 22 R 4 s 3
E 31 = Y 12 R λ q + Y 23 R 2 s φ X 11 R λ q 3 ψ 4 φ X 13 2 R 2 λ q 2 s + ψ X 33 R 3 λ q s 2 ψ + 2 φ X 23 2 R 4 s 3
E 32 = Y 23 R λ q + Y 22 R 2 s 3 ψ + 4 φ X 13 2 R λ q 3 + ( ψ + φ ) X 12 R 2 + ( 3 ψ + 4 φ ) X 33 R 2 λ q 2 s ( 7 ψ + 6 φ ) X 23 2 R 3 λ q s 2 + ψ X 22 R 4 s 3
E 33 = Y 22 R 2 + X 11 λ q 4 φ X 22 R 4 s 2 4 X 13 R λ q 3 s + 2 X 12 + 4 X 33 R 2 λ q 2 s 2 4 X 23 R 3 λ q s 3 + X 22 R 4 s 4 ρ h ω 2
At this point it is useful to remember that by imposing the parameters ( ψ = 0 , φ = 0 ) , we obtain the elements of matrix E for the Donnell shell theory; by imposing the parameters ( ψ = 1 , φ = 0 ) , we obtain the elements of matrix E for the Sanders shell theory; and by imposing the parameters ( ψ = 0 , φ = 1 ) , we obtain the elements of matrix E for the Flügge shell theory. See Expansion (2) for the middle surface change in the curvature k θ and torsion k x θ of the shell.
For the non-trivial solution, the determinant of the set of Equation (21) must be equal to zero [40]:
det E ( λ q , s , ω ) 3 × 3 = 0
Solving Equation (31), we get a third-order algebraic equation in ω 2 ; this last equation provides three different eigenfrequencies for each number of waves ( q , s ) that give three different vibration modes (i.e., longitudinal, torsional and radial modes). Since the highest natural frequency corresponds to the radial vibration mode, only the radial natural frequencies were computed in the numeric results.

6. Numeric Results

In this paper, the natural frequencies of SWCNTs were obtained in the framework of the Donnell, Sanders, and Flügge shell theories. An anisotropic elastic shell model was used to take into account the chirality effects of CNTs. Simply supported boundary conditions were imposed. Vibration modes with different number of waves along the longitudinal and circumferential directions were considered. SWCNTs with different values of radius R and aspect ratio L / R were investigated.
As known from the literature, two relevant open issues related to the continuous modelling of carbon nanotubes are due to their intrinsic anisotropic character and their discrete configuration. To this end, it is very important to adopt parameters and models able to correctly describe the actual molecular structure of carbon nanotubes.
In Table 1, the values of carbon-carbon bond parameters ( a , k ρ , k θ ) and equivalent continuous parameters ( h , ρ ) retrieved from the literature are reported. In particular, parameters k ρ and k θ , which denote force constants correlated to the variance of carbon-carbon bond length a and angle θ , respectively, were adopted to express the anisotropic elastic constants of SWCNTs via the molecular mechanics-based “stick-spiral model” developed by Chang [35].
Moreover, in order to study the dynamics of the actual discrete CNT via a continuous cylindrical shell, an equivalent thickness h , which is derived from MD simulations of CNT energy, and an equivalent mass density ρ , resulting from graphite surface density, were considered; see Ref. [16] for more details.

6.1. Comparison of the Shell Theories with the Results of Molecular Dynamics Simulations

In this section, the natural frequencies of the radial breathing mode ( q = 0 , s = 0 ) of the SWCNT of Table 1 with an aspect ratio L / R = 10 obtained by considering Donnell, Sanders, and Flügge shell theories are compared with the results of molecular dynamics simulations available in literature for different chirality indices ( n , m ) ; see Table 2.
Specifically, the results of the molecular dynamics simulations reported in Table 2 were retrieved from Ref. [10]. In that work, the free vibrations of armchair, zigzag, and chiral SWCNTs with different aspect ratios and diameters were studied via MM3 potential. This potential considers bond stretching, change in angles between adjacent bonds, torsion of the bond, van der Waals forces, and the coupling among stretching, bending, and torsional deformations. In particular, the energy due to the bond stretching has terms that are quadratic, cubic, and quartic in the bond length; thus, the strain energy due to the bond stretching is not an even function of the change in the bond length; see Ref. [10] for more details.
Obviously, the results of molecular dynamics simulations are the most correct since these simulations are able to correctly take into account the actual molecular structure of carbon nanotubes. However, the computational effort of molecular dynamics simulations is very high, in particular when dealing with very long and thin structures (such as carbon nanotubes), which have an extremely large number of atoms, and therefore, these results are available only for a reduced range of CNT geometries. To this end, it seems useful to investigate which shell theory is more accurate in the vibration modelling of carbon nanotubes considered as continuous homogeneous structures.
From Table 2, it can be observed that for all considered chirality indices, by assuming molecular dynamics results as reference, the percentage differences of the Flügge shell theory are the lowest, the percentage differences of the Donnell shell theory are the highest, and the Sanders shell theory gives an intermediate response.
These results could easily be predicted considering the strain–displacement Relationship (2), in which the Donnell shell theory has fewer terms than the other two (and therefore is the least accurate), and observing the force and moment Resultants (11)–(16), in which the Flügge shell theory includes more terms (and therefore is the most accurate).
Even if, among the three considered shell theories, Flügge is obtained to be the most accurate, it is preferable to not adopt this theory, because it has a relatively high computational effort, due to the large number of terms in the expansions of the force and moment resultants. This aspect is relevant not so much in the study of the linear vibrations (natural frequencies) but instead in the analysis of the non-linear vibrations (amplitude-frequency responses), where in the strain–displacement relationships the non-linear terms are also considered, and therefore the computational weight in solving the corresponding dynamic equations of motion becomes much greater.
To this aim, it could be useful to analyse whether the Donnell or Sanders shell theories are sufficiently accurate in the continuous modelling of SWCNT vibrations. Therefore, in the following, the natural frequencies of simply supported SWCNTs with different geometries and wavenumbers will be computed in the framework of the Donnell and Sanders shell theories, where the results of Flügge shell theory will be considered as reference.

6.2. Comparison of the Shell Theories for Different SWCNT Geometries and Wavenumbers

In the present section, the natural frequencies of the simply supported SWCNT of Table 1 obtained by considering the Donnell, Sanders, and Flügge shell theories were compared for different chirality indices ( n , m ) , aspect ratios L / R , numbers of longitudinal half-waves q , and circumferential waves s .
By taking as reference the results of the Flügge shell theory (which are the closest to those obtained from molecular dynamics simulations, see Table 2), we attempted to investigate the fields of applicability and limitations of the Donnell and Sanders shell theories.
Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 show the comparisons between Donnell and Flügge shell theories, and between Sanders and Flügge shell theories, carried out by adopting increasing values of chirality indices (SWCNT radius), specifically ( n = 5 , m = 5 ) (i.e., R = 0.34 nm), ( n = 10 , m = 10 ) (i.e., R = 0.68 nm), and ( n = 25 , m = 25 ) (i.e., R = 1.70 nm), and increasing values of SWCNT aspect ratio L / R = ( 10,15,20 ) are shown.
It must be underlined that this analysis was carried out within 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.5–10.0 nm, though most of the observed single-wall nanotubes have diameters < 2 nm”, see Ref. [9].
As for the Donnell shell theory, when increasing the radius R , the percentage difference with respect to the Flügge shell theory decreased for all longitudinal wavenumbers q , where the maximum value is obtained at q = 1 and decreased with increasing q ; the maximum peak corresponding to q = 1 was always located at s = 2 , whereas the peaks corresponding to q = ( 2 5 ) moved to higher circumferential wavenumbers.
On the other hand, when increasing the aspect ratio L / R , the percentage difference with respect to Flügge shell theory increased for all longitudinal wavenumbers q , where the maximum value was found again at q = 1 and decreased with increasing q , but now the maximum peak corresponding to q = 1 moved from s = 2 to s = 1 (lower circumferential wavenumber), whereas the peaks corresponding to q = ( 2 5 ) were always located at s = 2 .
It should be underlined that the increase in the percentage difference at the longitudinal wavenumber q = 1 obtained by increasing the aspect ratio L / R is extremely higher than the corresponding decrease at the longitudinal wavenumber q = 1 obtained by increasing the radius R (i.e., the effect of aspect ratio is prevalent with respect to radius).
As for Sanders shell theory, increasing both the radius R and the aspect ratio L / R , it was always found that the percentage difference with respect to the Flügge shell theory was relatively low for every number of longitudinal and circumferential waves ( < 1 % ) .
In particular, when increasing the radius R , the maximum percentage difference moved from q = 2 to q = 1 , decreasing its value (from 0.7 % to 0.1 % ); differently, when increasing the aspect ratio L / R , the maximum percentage difference moved from q = 2 to q = 5 , preserving its value ( 0.7 % ) (i.e., no effect of both aspect ratio and radius).
Therefore, by considering the parametric analyses presented in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, it can be observed that the Donnell shell theory cannot be applied for the vibration modelling of SWCNTs with relatively low radius R and relatively high aspect ratio L / R , and for the vibration modelling of modes with relatively low numbers of longitudinal q and circumferential s waves. However, as previously reported, the effect of the aspect ratio is prevalent with respect to the radius (and also to the wavenumbers) in providing the very high percentage difference obtained between the Donnell and Flügge shell theories. This is due to the different expression of the middle surface change in curvature k θ and torsion k x θ of the shell for the two theories (see the strain–displacement Equation (2)), which present more terms in the Flügge than in the Donnell shell theory, and in particular to the middle surface torsion, which is very sensitive to the value of aspect ratio.
Specifically, for the relatively low value of aspect ratio L / R = 10 , the maximum peak of percentage difference was located at the circumferential flexure mode ( q = 1 , s = 2 ) , where the effect of the middle surface change in the curvature k θ was prevalent; see Figure 2. On the other hand, for the relatively high value of aspect ratio L / R = 20 , the maximum peak of percentage difference was located at the beam-like mode ( q = 1 , s = 1 ) , where the effect of the middle surface torsion k x θ prevailed; see Figure 10.
Conversely, by considering the parametric analyses presented in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, it can be observed that Sanders shell theory is able to model with very good accuracy the linear vibrations of SWCNTs for all considered geometries and wavenumbers, and therefore it can be adopted instead of the more complex Flügge shell theory to compute the natural frequencies of SWCNTs. This is due to the fact that, in the Sanders shell theory, the effect of the aspect ratio is not present (as the strain–displacement Equation (2) is very similar to that of the Flügge shell theory), but only the effect of the radius occurs.
In fact, from the numerical simulations carried out in the present work, it was found that for both the Donnell and Sanders shell theories, the difference in the natural frequencies with respect to the Flügge shell theory decreased with CNT radius, see as e.g., Figure 2 compared to Figure 6 (Donnell vs. Flügge), and Figure 3 compared to Figure 7 (Sanders vs. Flügge). This behaviour can be understood by observing once more the expressions of the middle surface change in curvature k θ and torsion k x θ of the shell; see Equation (2). Since in these two expressions, the radius R is located in the denominator, then, by increasing the radius R , their value (i.e., their influence) reduces, and therefore, the natural frequencies of the three theories tend to become closer.
Finally, in Figure 12, six different mode shapes of the simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 10 are shown, where the radial breathing mode ( q = 0 , s = 0 ) corresponds to Rayleigh’s inextensional symmetrical mode (i.e., uniform vibration). Moreover, it also shows the axisymmetric mode ( q = 1 , s = 0 ) (no circumferential waves), beam-like mode ( q = 1 , s = 1 ) (one circumferential wave, characteristic of structures with a very long aspect ratio, such as beams), and shell-like modes ( q = 1 , s = 2 4 ) (two or more circumferential waves, characteristic of three-dimensional thin-walled structures, such as circular cylindrical shells). Such a graphical representation of the vibration modes could be useful for interpreting the previous results and comparisons.

7. Conclusions

In this paper, the natural frequencies of SWCNTs obtained in the framework of the Donnell, Sanders, and Flügge shell theories were compared. An anisotropic elastic shell model was adopted to take into account the intrinsic chirality effects of CNTs. Simply supported boundary conditions were imposed. Vibration modes with different wavenumbers along the longitudinal and circumferential directions were studied. SWCNTs with different values of radius R and aspect ratio L / R were considered. The most important findings of the present paper are reported below.
  • By means of comparisons with the results of molecular dynamic simulations reported in the literature, it was derived that the Flügge shell theory is the most accurate in the computation of the natural frequencies of SWCNTs.
  • Since the Flügge shell theory requires a very high computational effort due to the large number of terms in the equations of motion, it was investigated whether a simpler shell theory is able to model with sufficient accuracy the linear vibrations of SWCNTs.
  • It was found that the Donnell shell theory was not accurate for relatively low longitudinal and circumferential wavenumbers, for relatively low diameters, and for relatively high aspect ratios, and therefore, it is not able to properly model SWCNT vibrations.
  • On the other hand, it was found that the Sanders shell theory was very accurate for all considered geometries and wavenumbers, and therefore, it can be adopted instead of the more complex Flügge shell theory to properly model SWCNT vibrations.
As first application of the results obtained in the present paper, the authors are planning to write a new manuscript on the effects of nonlocal elasticity and strain gradient on the linear vibrations of SWCNTs by considering an anisotropic elastic model in the framework of the Sanders shell theory.
A second relevant application of the findings of the present work regards the nonlinear vibrations of CNTs. It has been proven in the literature that, in the presence of a large number of carbon atoms, i.e., for relatively long or multi-walled carbon nanotubes, molecular dynamics simulations require higher computational effort than all equivalent continuous models, including the Flügge shell theory [3]. Moreover, by carrying out numerical simulations in linear field, it was found that all equivalent continuous models, among which are the Donnell, Sanders, and Flügge shell theories, are equally time costing [23]. The problem with the Flügge shell theory is due to the transition from linear to nonlinear analysis. On the one hand, the additional terms present in the expressions of forces and moments give the Flügge shell theory greater accuracy than that of Sanders and Donnell. On the other hand, these terms lead to very high computational effort in the numerical simulations in the nonlinear field, which are necessary to study the actual dynamic behaviour of carbon nanotubes and to investigate fluid–structure interactions. For this reason, especially for the nonlinear analyses, it is preferable to use the Sanders shell theory, since it is more accurate than the Donnell one (see the results of the present paper) and less computationally expensive than the Flügge one.

Author Contributions

Conceptualization: M.S. and I.E.E.; data curation: M.S. and M.B.; funding acquisition: M.C., R.R. and E.R.; investigation: M.S. and M.B.; methodology: M.S. and I.E.E.; project administration: M.S. and I.E.E.; supervision: M.C., R.R. and E.R.; visualization: M.S. and M.B.; writing—original draft: M.S., M.B. and I.E.E.; writing—review & editing: R.R., M.C. and E.R. All authors have read and agreed to the published version of the manuscript.

Funding

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

Data Availability Statement

All data are available from the authors.

Acknowledgments

Authors M. Strozzi, E. Radi, M. Cocconcelli and R. Rubini are grateful to the Department of Sciences and Methods for Engineering, University of Modena and Reggio Emilia, Reggio Emilia, Italy (Grant 020142_22_FRN_SOSTEGNO_RICERCA_DISMI) for the financial support of this work. The present work is dedicated to the blessed memory of Leonid I. Manevitch, outstanding educator, respected professor, and admired scholar who died, after severe illness, on 20 August 2020, at the age of 82. The great passion and enthusiasm towards scientific research that Manevitch nurtured until the end of his life, in particular towards the much-loved carbon nanotubes, will remain forever in the heart of his collaborators, and will be for them a strong and ever-present motivation to progress in the knowledge of these wonderful structures.

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. Elishakoff, I.E.; Pentaras, D.; Dujat, K.; Versaci, C.; Muscolino, G.; Storch, J.; Bucas, S.; Challamel, N.; Natsuki, T.; Zhang, Y. Carbon Nanotubes and Nanosensors: Vibration, Buckling and Ballistic Impact; John Wiley & Sons: London, UK, 2012. [Google Scholar]
  3. Jorio, A.; Dresselhaus, G.; Dresselhaus, M. Carbon Nanotubes. In Advanced Topics in the Synthesis, Structure, Properties and Applications; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  4. Marulanda, J.M. Carbon Nanotubes: Applications on Electron Devices; InTech Open: Rijeka, Croatia, 2011. [Google Scholar]
  5. Mahar, B.; Laslau, C.; Yip, R.; Sun, Y. Development of Carbon Nanotube-Based Sensors. A Review. IEEE Sens. J. 2007, 7, 266–284. [Google Scholar] [CrossRef]
  6. Hierold, C.; Jungen, A.; Stampfer, C.; Helbling, T. Nano electromechanical sensors based on carbon nanotubes. Sens. Actuators A 2007, 136, 51–61. [Google Scholar] [CrossRef]
  7. Rao, A.; Richter, E.; Bandow, S.; Chase, B.; Eklund, P.C.; Williams, K.; Fang, S.; Subbaswamy, K.R.; Menon, M.; Thess, A. Diameter-Selective Raman Scattering from Vibrational Modes in Carbon Nanotubes. Science 1997, 275, 187–191. [Google Scholar] [CrossRef] [PubMed]
  8. Bandow, S.; Asaka, S.; Saito, Y.; Rao, A.; 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]
  9. 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]
  10. Gupta, S.S.; Bosco, F.G.; Batra, R.C. Wall thickness and elastic moduli of single-walled carbon nanotubes from frequencies of axial, torsional and inextensional modes of vibration. Comput. Mater. Sci. 2010, 47, 1049–1059. [Google Scholar] [CrossRef]
  11. 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]
  12. Duan, W.; 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]
  13. 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]
  14. Arroyo, M.; Belytschko, T. Continuum Mechanics Modeling and Simulation of Carbon Nanotubes. Meccanica 2005, 40, 455–469. [Google Scholar] [CrossRef]
  15. 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]
  16. 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]
  17. 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]
  18. 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]
  19. 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]
  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. Leissa, A.W. Vibration of Shells; Acoustical Society of America: Columbus, OH, USA, 1993. [Google Scholar]
  22. Yamaki, N. Elastic Stability of Circular Cylindrical Shells; Elsevier: Amsterdam, The Netherlands, 1984. [Google Scholar]
  23. Amabili, M. Nonlinear Vibrations and Stability of Shells and Plates; Cambridge University Press: New York, NY, USA, 2008. [Google Scholar]
  24. Soedel, W. Vibrations of Shells and Plates, 3rd ed.; Marcel Dekker: New York, NY, USA, 2004. [Google Scholar]
  25. Ventsel, E.; Krauthammer, T. Thin Plates and Shells: Theory, Analysis, and Applications; The Pennsylvania State University, Marcel Dekker: New York, NY, USA, 2001. [Google Scholar]
  26. Amabili, M. A comparison of shell theories for large-amplitude vibrations of circular cylindrical shells: Lagrangian approach. J. Sound Vib. 2003, 264, 1091–1125. [Google Scholar] [CrossRef]
  27. Wang, C.M.; Tay, Z.Y.; Chowdhury, A.R.; Duan, W.H.; Zhang, Y.Y.; Silvestre, N. Examination of cylindrical shell theories for buckling of carbon nanotubes. Int. J. Struct. Stab. Dyn. 2011, 11, 1035–1058. [Google Scholar] [CrossRef]
  28. Loy, C.T.; Lam, K.Y.; Reddy, J.N. Vibration of functionally graded cylindrical shells. Int. J. Mech. Sci. 1999, 41, 309–324. [Google Scholar] [CrossRef]
  29. Pradhan, S.C.; Loy, C.T.; Lam, K.Y.; Reddy, J.N. Vibration characteristics of functionally graded cylindrical shells under various boundary conditions. Appl. Acoust. 2000, 61, 111–129. [Google Scholar] [CrossRef]
  30. Arshad, S.H.; Naeem, M.N.; Sultana, N. Frequency analysis of functionally graded material cylindrical shells with various volume fractions laws. J. Mech. Eng. Sci. 2007, 221, 1483–1495. [Google Scholar] [CrossRef]
  31. Carrera, E. Historical review of Zig-Zag theories for multi-layered plates and shells. Appl. Mech. Rev. 2003, 56, 287–308. [Google Scholar] [CrossRef]
  32. Dumir, P.C.; Nath, J.K.; Kumari, P.; Kapuria, S. Improved Efficient Zigzag and Third Order Theories for Circular Cylindrical Shells Under Thermal Loading. J. Therm. Stress. 2008, 31, 343–367. [Google Scholar] [CrossRef]
  33. Rahmani, O.; Khalili, S.M.R.; Malekzadeh, K. Free vibration response of composite sandwich cylindrical shell with flexible core. Comp. Struct. 2010, 92, 1269–1281. [Google Scholar] [CrossRef]
  34. Lakis, A.; Selmane, A.; Toledano, A. Non-linear free vibration analysis of laminated orthotropic cylindrical shells. Int. J. Mech. Sci. 1998, 40, 27–49. [Google Scholar] [CrossRef]
  35. 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]
  36. Chang, T. A molecular based anisotropic shell model for single-walled carbon nanotubes. J. Mech. Phys. Solids 2010, 58, 1422–1433. [Google Scholar] [CrossRef]
  37. 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]
  38. 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]
  39. Ghavanloo, E.; Fazelzadeh, S. Vibration characteristics of single-walled carbon nanotubes based on an anisotropic elastic shell model including chirality effect. Appl. Math. Modell. 2012, 36, 4988–5000. [Google Scholar] [CrossRef]
  40. Fazelzadeh, S.; Ghavanloo, E. Nonlocal anisotropic elastic shell model for vibrations of single-walled carbon nanotubes with arbitrary chirality. Comp. Struct. 2012, 94, 1016–1022. [Google Scholar] [CrossRef]
Figure 1. Continuous modelling of SWCNTs: (a) actual discrete SWCNT; (b) geometry of the equivalent continuous circular cylindrical shell; and (c) cross-section of the surface of the equivalent continuous circular cylindrical shell [20].
Figure 1. Continuous modelling of SWCNTs: (a) actual discrete SWCNT; (b) geometry of the equivalent continuous circular cylindrical shell; and (c) cross-section of the surface of the equivalent continuous circular cylindrical shell [20].
Nanomaterials 13 01390 g001
Figure 2. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 2. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g002
Figure 3. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 3. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g003
Figure 4. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 10 , m = 10 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 4. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 10 , m = 10 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g004
Figure 5. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 10 , m = 10 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 5. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 10 , m = 10 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g005
Figure 6. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 25 , m = 25 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 6. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 25 , m = 25 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g006
Figure 7. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 25 , m = 25 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 7. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 25 , m = 25 ) and aspect ratio L / R = 10 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g007
Figure 8. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 15 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 8. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 15 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g008
Figure 9. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 15 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 9. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 15 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g009
Figure 10. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 20 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 10. Percentage differences between the radial natural frequencies obtained via Donnell and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 20 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g010
Figure 11. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 20 . Number of longitudinal half-waves q . Number of circumferential waves s .
Figure 11. Percentage differences between the radial natural frequencies obtained via Sanders and Flügge shell theories (Flügge as the reference). Anisotropic elastic shell model. Simply supported SWCNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 20 . Number of longitudinal half-waves q . Number of circumferential waves s .
Nanomaterials 13 01390 g011
Figure 12. Mode shapes of the simply supported SWNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 10 : (a) Radial breathing mode ( q = 0 , s = 0 ) ; (b) axisymmetric mode ( q = 1 , s = 0 ) ; (c) beam-like mode ( q = 1 , s = 1 ) ; (d) shell-like mode ( q = 1 , s = 2 ) ; (e) shell-like mode ( q = 1 , s = 3 ) ; and (f) shell-like mode ( q = 1 , s = 4 ) .
Figure 12. Mode shapes of the simply supported SWNT of Table 1 with chirality indices ( n = 5 , m = 5 ) and aspect ratio L / R = 10 : (a) Radial breathing mode ( q = 0 , s = 0 ) ; (b) axisymmetric mode ( q = 1 , s = 0 ) ; (c) beam-like mode ( q = 1 , s = 1 ) ; (d) shell-like mode ( q = 1 , s = 2 ) ; (e) shell-like mode ( q = 1 , s = 3 ) ; and (f) shell-like mode ( q = 1 , s = 4 ) .
Nanomaterials 13 01390 g012
Table 1. Mechanical parameters adopted in the anisotropic elastic continuous shell model [16,35].
Table 1. Mechanical parameters adopted in the anisotropic elastic continuous shell model [16,35].
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
Shell thickness h (nm)0.0665
Mass density ρ (kg/m3)11,700
Table 2. Natural frequencies of the radial breathing mode ( q = 0 , s = 0 ) of the SWCNT of Table 1 with aspect ratio. L / R = 10 . Comparisons among Donnell, Sanders, Flügge shell theories and molecular dynamics simulations.
Table 2. Natural frequencies of the radial breathing mode ( q = 0 , s = 0 ) of the SWCNT of Table 1 with aspect ratio. L / R = 10 . Comparisons among Donnell, Sanders, Flügge shell theories and molecular dynamics simulations.
Natural   Frequency   ω R B M (cm−1)Difference (%)
Chirality   Indices   ( n , m ) Molecular Dynamics Simulations [10]Donnell Shell TheorySanders Shell TheoryFlügge Shell Theory
(10, 0)290.4632.381.320.42
(9, 6)221.4963.471.930.43
(8, 8)209.0084.052.250.51
(9, 9)185.8964.032.240.49
(16, 0)181.7473.101.720.04
(10, 10)167.3774.012.230.48
(11, 11)152.2074.002.220.47
(20, 0)145.3633.351.860.11
(22, 7)110.8083.712.060.32
(25, 10)93.1533.561.980.21
(30, 5)88.7503.511.950.16
(36, 5)75.1183.471.930.17
(33, 16)67.2213.561.980.23
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Strozzi, M.; Elishakoff, I.E.; Bochicchio, M.; Cocconcelli, M.; Rubini, R.; Radi, E. A Comparison of Shell Theories for Vibration Analysis of Single-Walled Carbon Nanotubes Based on an Anisotropic Elastic Shell Model. Nanomaterials 2023, 13, 1390. https://doi.org/10.3390/nano13081390

AMA Style

Strozzi M, Elishakoff IE, Bochicchio M, Cocconcelli M, Rubini R, Radi E. A Comparison of Shell Theories for Vibration Analysis of Single-Walled Carbon Nanotubes Based on an Anisotropic Elastic Shell Model. Nanomaterials. 2023; 13(8):1390. https://doi.org/10.3390/nano13081390

Chicago/Turabian Style

Strozzi, Matteo, Isaac E. Elishakoff, Michele Bochicchio, Marco Cocconcelli, Riccardo Rubini, and Enrico Radi. 2023. "A Comparison of Shell Theories for Vibration Analysis of Single-Walled Carbon Nanotubes Based on an Anisotropic Elastic Shell Model" Nanomaterials 13, no. 8: 1390. https://doi.org/10.3390/nano13081390

APA Style

Strozzi, M., Elishakoff, I. E., Bochicchio, M., Cocconcelli, M., Rubini, R., & Radi, E. (2023). A Comparison of Shell Theories for Vibration Analysis of Single-Walled Carbon Nanotubes Based on an Anisotropic Elastic Shell Model. Nanomaterials, 13(8), 1390. https://doi.org/10.3390/nano13081390

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