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
, length
, and thickness
. A cylindrical coordinate system
is adopted, where the origin
of the reference system is located at the centre of one end of the cylindrical shell. Three displacements are present: longitudinal
, circumferential
, and radial
, where the radial displacement
is assumed as positive outward,
are the longitudinal and angular coordinates of an arbitrary point on the middle surface of the shell,
is the radial coordinate along the thickness
, and
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
of the shell is small with respect to the radius of curvature
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
of the shell are related to the displacements
in the form [
21]:
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
of the shell are related to the displacements
in the form [
21]:
From Equation (2), it can be noted that the middle surface change in the curvature and torsion are written in a different form for the three different shell theories, where parameters denote the Donnell shell theory, parameters denote the Sanders shell theory, and parameters 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
at an arbitrary point of the surface of the shell are related to the middle surface strains
and to the changes in curvature and torsion of the middle surface
by the relationships [
23]:
where
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:
where Equation (4) relate the strain components at an arbitrary point of the shell surface
to the displacements
.
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
at both ends
of the SWCNT. The displacement field that satisfies these boundary conditions can be written as [
39]:
where
denote the displacement amplitudes along the longitudinal
, circumferential
, and radial
directions, respectively;
is the imaginary unit;
is the wavenumber along the longitudinal direction, with
, where
is the number of longitudinal half-waves and
is the length of the SWCNT;
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
is obtained, which can be written in the form [
39]:
where
is a non-symmetric matrix, whose elements are:
At this point it is useful to remember that by imposing the parameters , we obtain the elements of matrix for the Donnell shell theory; by imposing the parameters , we obtain the elements of matrix for the Sanders shell theory; and by imposing the parameters , we obtain the elements of matrix for the Flügge shell theory. See Expansion (2) for the middle surface change in the curvature and torsion of the shell.
For the non-trivial solution, the determinant of the set of Equation (21) must be equal to zero [
40]:
Solving Equation (31), we get a third-order algebraic equation in ; this last equation provides three different eigenfrequencies for each number of waves 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 and aspect ratio 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
and equivalent continuous parameters
retrieved from the literature are reported. In particular, parameters
and
, which denote force constants correlated to the variance of carbon-carbon bond length
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
, 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
of the SWCNT of
Table 1 with an aspect ratio
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
; 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
, aspect ratios
, numbers of longitudinal half-waves
, and circumferential waves
.
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
(i.e.,
nm),
(i.e.,
nm), and
(i.e.,
nm), and increasing values of SWCNT aspect ratio
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 , the percentage difference with respect to the Flügge shell theory decreased for all longitudinal wavenumbers , where the maximum value is obtained at and decreased with increasing ; the maximum peak corresponding to was always located at , whereas the peaks corresponding to moved to higher circumferential wavenumbers.
On the other hand, when increasing the aspect ratio , the percentage difference with respect to Flügge shell theory increased for all longitudinal wavenumbers , where the maximum value was found again at and decreased with increasing , but now the maximum peak corresponding to moved from to (lower circumferential wavenumber), whereas the peaks corresponding to were always located at .
It should be underlined that the increase in the percentage difference at the longitudinal wavenumber obtained by increasing the aspect ratio is extremely higher than the corresponding decrease at the longitudinal wavenumber obtained by increasing the radius (i.e., the effect of aspect ratio is prevalent with respect to radius).
As for Sanders shell theory, increasing both the radius and the aspect ratio , 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 .
In particular, when increasing the radius , the maximum percentage difference moved from to , decreasing its value (from to ); differently, when increasing the aspect ratio , the maximum percentage difference moved from to , preserving its value (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
and relatively high aspect ratio
, and for the vibration modelling of modes with relatively low numbers of longitudinal
and circumferential
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
and torsion
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
, the maximum peak of percentage difference was located at the circumferential flexure mode
, where the effect of the middle surface change in the curvature
was prevalent; see
Figure 2. On the other hand, for the relatively high value of aspect ratio
, the maximum peak of percentage difference was located at the beam-like mode
, where the effect of the middle surface torsion
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
and torsion
of the shell; see Equation (2). Since in these two expressions, the radius
is located in the denominator, then, by increasing the radius
, 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
and aspect ratio
are shown, where the radial breathing mode
corresponds to Rayleigh’s inextensional symmetrical mode (i.e., uniform vibration). Moreover, it also shows the axisymmetric mode
(no circumferential waves), beam-like mode
(one circumferential wave, characteristic of structures with a very long aspect ratio, such as beams), and shell-like modes
(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 and aspect ratio 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.