# Nonlinear Resonance Interaction between Conjugate Circumferential Flexural Modes in Single-Walled Carbon Nanotubes

Authors: Matteo Strozzi; Francesco Pellicano
Journal: Shock and Vibration
 (2019)

---

## Abstract

This paper presents an investigation on the dynamical properties of single-walled carbon nanotubes (SWCNTs), and nonlinear modal interaction and energy exchange are analysed in detail. Resonance interactions between two conjugate circumferential flexural modes (CFMs) are investigated. The nanotubes are analysed through a continuous shell model, and a thin shell theory is used to model the dynamics of the system; free-free boundary conditions are considered. The Rayleigh–Ritz method is applied to approximate linear eigenfunctions of the partial differential equations that govern the shell dynamics. An energy approach, based on Lagrange equations and series expansion of the displacements, is considered to reduce the initial partial differential equations to a set of nonlinear ordinary differential equations of motion. The model is validated in linear field (natural frequencies) by means of comparisons with literature. A convergence analysis is carried out in order to obtain the smallest modal expansion able to simulate the nonlinear regimes. The time evolution of the nonlinear energy distribution over the SWCNT surface is studied. The nonlinear dynamics of the system is analysed by means of phase portraits. The resonance interaction and energy transfer between the conjugate CFMs are investigated. A travelling wave moving along the circumferential direction of the SWCNT is observed.


---

## Body

## 1. Introduction

Carbon nanotubes (CNTs) were first discovered in 1991 by Iijima in the laboratories of the NEC Corporation in Japan. These “helical microtubules of graphitic carbon” were observed in the form of needles comprising several coaxial tubes of graphitic sheets by means of a transmission electron microscopy (TEM). They were grown on the negative end of the carbon electrode used in the arc-discharge synthesis of fullerenes and were detected by high-resolution electron micrographs [1].Starting from their discovery, CNTs have attracted increasing interest from researchers all over the world because of their outstanding mechanical properties, in particular, very high Young’s modulus (of the order of 1-2 TPa) and tensile strength (up to 100 GPa) [2–4].The combination of the two previous mechanical properties allows CNTs to reach natural vibration frequencies of the order of THz: these ultrahigh frequencies have led to experiment the application of CNTs and their composites as ultrahigh sensitivity resonators within nano-electro-mechanical-systems (NEMS), such as oscillators, sensors, and charge detectors [5–10].Several methodologies have been adopted in order to better investigate the mechanical properties of CNTs and their application in many industrial fields. Experimental tests, numerical simulations, and continuum mechanics studies have been performed. Single-walled and multi‐walled carbon nanotubes (MWCNTs) have been considered.Experiments based on the resonant Raman spectroscopy (RRS) technology have been carried out to obtain atomic structure and chirality of SWCNTs considering the frequency spectrum of the CNT and applying the theory of the resonant transitions. This method starts from the measurement of the SWCNT diameter obtained by means of atomic force microscopy (AFM) and investigates natural frequency and energy of the SWCNT radial breathing mode (RBM) [11–15].However, owing to their high technological complexity and costs, experimental methods cannot be considered as efficient approaches to investigate the mechanical behaviour of CNTs.Molecular dynamics simulations (MDS) have been performed to get numerically the time evolution of SWCNT as a system of atoms treated as point-like masses interacting with one another according to an assumed potential energy which describes the atomic bond interactions along the CNT. The vibrations of the free atoms of the SWCNT are recorded for a certain duration at fixed temperature and the corresponding natural frequencies are computed by the Discrete Fourier Transform (DFT) [4, 16–19]. Modelling the SWCNT as an equivalent space frame-like discrete structure, on the basis of atomistic structural mechanics (MDS), is computationally very expensive and this method cannot be easily applied to the structural simulation of CNTs, which incorporate a large number of atoms.Theoretical models based on continuum mechanics are computationally more efficient than MDS. In such kind of modelling, CNTs are treated as continuous homogeneous structures, ignoring their intrinsic atomic nature; this allows the achievement of a strong reduction of the number of degrees of freedom. In order to investigate the vibration properties of SWCNTs and MWCNTs, both beam-like [20–32] and shell-like [33–37] equivalent continuous models have been proposed in the past.In particular, the analogy between circular cylindrical shell and CNT nanostructures led to extensive application of continuous shell models for the CNT vibration analysis. To deepen the knowledge on shells, the fundamental books of Leissa [38] and Yamaki [39] are suggested to the reader. In these works, nonlinear dynamics and stability of shells with different shapes, material properties, and boundary conditions are treated; the most important thin shell theories (Donnell, Sanders-Koiter, Flugge-Lur’e-Byrne) are reported, numerical and experimental results are described.In order to study the CNT dynamics through continuous models, equivalent mechanical parameters (i.e., Young’s modulus, Poisson’s ratio, and mass density) must be applied. Indeed, we are going to simulate the intrinsically discrete CNT with its continuous twin: a circular cylindrical shell; in order to be dynamically equivalent all shell parameters must be suitably set, this was achieved comparing the equations of the classical shell theory for tensile and bending rigidity with the MDS results [40].One of the most controversial topics in modelling the discrete CNTs as continuous cylindrical shells is denoted by the inclusion of the size effects, i.e., surface stresses, strain gradients, and nonlocalities, into the elastic shell theory considered. It was found in the literature [41] that the size effects influence the vibration characteristics of the CNTs in the higher region of their frequency spectrum (vibration modes with relatively high number of longitudinal half-waves). Therefore, to obtain relatively high resonant frequencies correctly, nonlocal elastic shell models were developed, where the size effects are taken into account in the stress-strain relationships [42–44].Both theoretical and experimental aspects of nonlinear dynamics and stability of shells are widely explored in the literature [45–49]. In these works, the effect of the geometry, boundary conditions, fluid-structure interaction, material distribution, geometric imperfections, and external loads on the nonlinear vibrations of circular cylindrical shells is investigated. Particular attention was paid to the choice of modes used in the expansion of the displacements; the general goal is to simulate the nonlinear dynamic behaviour of the shell with a minimum number of degrees of freedom: these models were obtained by means of convergence analyses of the nonlinear response of the system, see [50–52] for more details.It must be stressed that the application of the CNTs and their composites as ultra-high-sensitivity resonators within NEMS requires deep investigation of the nonstationary dynamics phenomena, i.e., energy transfers, internal resonances, and propagating waves, in the CNTs, considered as strongly nonlinear systems. To this aim, the nonlinear resonance interaction and energy exchange between vibration modes in SWCNTs were extensively studied in the literature [53–57].An interesting nonstationary dynamic phenomenon, a travelling wave response, was observed both numerically and experimentally in circular cylindrical shells, with and without fluid-structure interaction [58, 59]. This phenomenon is due to the coexistence of three specific conditions: (i) axisymmetric structure, (ii) nonlinear vibrations, and (iii) conjugate modes. Shells with circular cylindrical geometry are symmetric with respect to the longitudinal axis; in presence of large amplitude displacements, they are subjected to nonlinear vibrations; conjugate modes are couples of modes having the same mode shape and frequency, but angularly shifted of a quarter of period in circumferential direction. Therefore, in circular cylindrical shell, the nonlinear resonance interaction between conjugate modes can give rise to a pure travelling wave moving along the circumferential direction of the shell. To the authors’ best knowledge, the present paper is the first work in which the travelling wave response in SWCNT, considered as “graphene sheet rolled into a seamless tube” [1], is observed in nonlinear field. Other relevant linear and nonlinear dynamic phenomena of CNTs can be found in [60–62].In this paper, the nonlinear oscillations and energy distribution of SWCNTs are investigated. The resonance interaction between two CFMs having the same natural frequency and with mode shapes shifted ofπ/4 (conjugate modes) is analysed. The Sanders-Koiter shell theory is applied to model the nonlinear dynamics of the system in case of finite amplitude vibration. The CNT deformation is given in terms of longitudinal, circumferential, and radial displacement field. Free-free boundary conditions are considered. In the linear analysis, the three displacements are expanded by using a double mixed series and the Rayleigh–Ritz method is applied to get approximated eigenfunctions. In the nonlinear analysis, the three displacements are re‐expanded by considering the approximated eigenfunctions of the linear analysis, and an energy method based on the Lagrange equations is used to reduce the partial differential equations of motion to a system of ordinary differential equations, which is then solved by means of the implicit Runge–Kutta numerical scheme. The model is validated by comparing the natural frequencies with the results of experiments and molecular dynamics simulations found in the literature. A convergence analysis is carried out to select the vibration modes to be added to the two directly excited CFMs providing a proper description of the nonlinear response of the system. The evolution in time of the nonlinear energy distribution of the conjugate CFMs over the SWCNT surface is studied. Phase portraits describing the nonlinear dynamics of the system are found. The effect of the conjugate mode resonance interaction on the nonlinear behaviour of the SWCNT is investigated. The participation of both conjugate modes gives a pure travelling wave moving circumferentially around the CNT.


## 2. Sanders-Koiter Shell Theory

In this work, the Sanders-Koiter elastic thin shell theory is used to model the SWCNT dynamics.It must be stressed that in the present paper, the contribution of the size effects (surface stresses, strain gradients, and nonlocalities) is not taken into account in the stress-strain relationships.Since size effects have a relevant influence in the higher region of the SWCNT frequency spectrum (vibration modes with relatively high number of longitudinal half-waves and resonant frequencies), then they can be neglected in this manuscript, which is focused on the lower region of the SWCNT frequency spectrum (conjugate CFMs with relatively low number of longitudinal half-waves and resonant frequencies). Therefore, in the present work, the local continuum mechanics neglecting the size effects can be applied correctly without loss of accuracy (see ref. [41] for more details).In Figure1(a), a picture of a discrete SWCNT is added to better describe the transition from the actual atomistic structure of the CNT to the equivalent continuous shell model adopted in this work.Figure 1
Equivalent continuum structure of a SWCNT. (a) Discrete SWCNT [19]. (b) Geometry of the equivalent circular cylindrical shell. (c) Cross section of the circular cylindrical shell surface.In Figures1(b) and 1(c), a circular cylindrical shell having radius R, length L, and thickness h is represented; a cylindrical coordinate system (O; x, θ, z) is considered where the origin O of the reference system is located at the centre of one end of the shell. In Figures 1(b) and 1(c), three displacements are represented: longitudinal u (x, θ, t), circumferential v (x, θ, t), and radial w (x, θ, t); the radial displacement w is considered positive outward; x and θ 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; t is the time.


### 2.1. Strain-Displacement Relationships

Since one of the most important issues of CNT continuum modelling is given by the infinitesimal dimensions, which cause very fast time scales and therefore induce hard computational troubles during numerical integration, then, in this work, the governing parameters, i.e., shell displacements and time, are transformed into dimensionless form.The dimensionless displacement fieldu˜,v˜,w˜ of the shell can be written in the form [37](1)u˜=uR,v˜=vR,w˜=wR,where u,v,w is the dimensional displacement field and R is the radius of the shell.The Sanders-Koiter shell theory is based on Love’s first approximation, which considers the following five hypotheses [38]: (i) the thickness of the shell is very small in comparison with the other dimensions, (ii) the deformations are small, (iii) the transverse normal stress is small with respect to the other stresses and can be neglected (σz = 0), (iv) the normal to the reference undeformed shell surface remains straight and normal to the deformed shell surface, and (v) the normal to the reference undeformed surface suffers no extension during deformation, i.e., no thickness stretching is present (Kirchhoff–Love kinematic hypothesis).The consequences of the last two assumptions are that the transverse shear the deformations of the shell are neglected (γxz = γθz = 0); in addition, also the rotary inertia of the shell is neglected.The dimensionless middle surface strains of the shellε˜x,0,ε˜θ,0,γ˜xθ,0 are expressed as [54](2)ε˜x,0=α∂u˜∂η+12α2∂w˜∂η2+18α∂v˜∂η−∂u˜∂θ2,(3)ε˜θ,0=∂v˜∂θ+w˜+12∂w˜∂θ−v˜2+18∂u˜∂θ−α∂v˜∂η2,(4)γ˜xθ,0=∂u˜∂θ+α∂v˜∂η+α∂w˜∂η∂w˜∂θ−v˜,where the linear terms are inserted into square brackets, the nonlinear terms are inserted into curly brackets, η = x/L is the dimensionless longitudinal coordinate of the shell and α = R/L.The separation of the linear and nonlinear terms in equations (2)–(4) can be useful in order to properly insert the strain contributions into the equation of the dimensionless elastic strain energy of the shell in the following linear and nonlinear vibration analyses.The dimensionless middle surface changes in curvature and torsion of the shellk˜x,k˜θ,k˜xθ are linear terms given by [54](5)k˜x=−α2∂2w˜∂η2,k˜θ=∂v˜∂θ−∂2w˜∂θ2,k˜xθ=−2α∂2w˜∂η∂θ+32α∂v˜∂η−12∂u˜∂θ.


### 2.2. Force and Moment Resultants

The dimensionless force (N˜x,N˜θ,N˜xθ,Q˜x,Q˜θ) and moment (M˜x,M˜θ,M˜xθ) resultants per unit length of the shell are expressed as [37](6)N˜x=11−ν2ε˜x,0+νε˜θ,0,N˜θ=11−ν2ε˜θ,0+νε˜x,0,N˜xθ=121+νγ˜xθ,0,(7)M˜x=β121−ν2k˜x+νk˜θ,M˜θ=β121−ν2k˜θ+νk˜x,M˜xθ=β241+νk˜xθ,(8)Q˜x=β2121−ν2αk˜x,x+νk˜θ,x+1−ν2k˜xθ,θ,Q˜θ=β2121−ν21−ν2αk˜xθ,x+k˜θ,θ+νk˜x,θ,where ν is Poisson’s ratio of the shell and β = h/R. These expressions of the dimensionless force and moment resultants will be used into the equations of the natural-type boundary conditions.


### 2.3. Elastic Strain Energy

In the case of homogeneous and isotropic material, the dimensionless elastic strain energyU˜ of a cylindrical shell, by neglecting the transverse normal stress σz (Love’s first approximation), can be written as a function of the dimensionless strains ε˜x,0,ε˜θ,0,γ˜xθ,0 and k˜x,k˜θ,k˜xθ in the form [37](9)U˜η,θ,τ=1211−ν2∫01∫02πε˜x,02+ε˜θ,02+2νε˜x,0ε˜θ,0+1−ν2γ˜xθ,02dηdθ+β212∫01∫02πk˜x2+k˜θ2+2νk˜xk˜θ+1−ν2k˜xθ2dηdθ,where the first term of the right-hand side of equation (9) is the membrane energy (also referred to stretching energy) and the second one is the bending energy.Another relevant issue of the thin-shell modelling of SWCNTs is given by the choice of an isotropic or anisotropic model; indeed, CNTs are discrete systems, i.e., they are intrinsically nonisotropic. However, it was proven in the past that the use of isotropic models for SWCNTs does not lead to significant errors. Therefore, in the present work, an isotropic elastic shell model is considered.


### 2.4. Kinetic Energy

The dimensionless time variableτ is defined introducing a reference dimensional circular frequency ω0 (i.e., the lowest extensional circular frequency of a ring in plane strain) in the form [56](10)τ=ω0t,ω0=E1−ν2ρR2,where t is the dimensional time and E and ρ are Young’s modulus and mass density of the shell, respectively.The dimensionless velocity fieldu˜′,v˜′,w˜′ of the shell is given by [56](11)u˜′=du˜dτ=u˙Rω0,v˜′=dv˜dτ=v˙Rω0,w˜′=dw˜dτ=w˙Rω0,where u˙,v˙,w˙ is the dimensional velocity field.The dimensionless kinetic energyT˜ of a cylindrical shell (rotary inertia is neglected) is [56](12)T˜η,θ,τ=12γ∫01∫02πu˜′2+v˜′2+w˜′2dηdθ,where γ=ρR2ω02/E.


### 2.5. Total Energy

The dimensionless total energyẼ of the shell can be expressed as the sum of the dimensionless elastic strain energy Ũ and kinetic energy T˜ in the following form:(13)E˜η,θ,τ=U˜η,θ,τ+T˜η,θ,τ.


## 2.1. Strain-Displacement Relationships

Since one of the most important issues of CNT continuum modelling is given by the infinitesimal dimensions, which cause very fast time scales and therefore induce hard computational troubles during numerical integration, then, in this work, the governing parameters, i.e., shell displacements and time, are transformed into dimensionless form.The dimensionless displacement fieldu˜,v˜,w˜ of the shell can be written in the form [37](1)u˜=uR,v˜=vR,w˜=wR,where u,v,w is the dimensional displacement field and R is the radius of the shell.The Sanders-Koiter shell theory is based on Love’s first approximation, which considers the following five hypotheses [38]: (i) the thickness of the shell is very small in comparison with the other dimensions, (ii) the deformations are small, (iii) the transverse normal stress is small with respect to the other stresses and can be neglected (σz = 0), (iv) the normal to the reference undeformed shell surface remains straight and normal to the deformed shell surface, and (v) the normal to the reference undeformed surface suffers no extension during deformation, i.e., no thickness stretching is present (Kirchhoff–Love kinematic hypothesis).The consequences of the last two assumptions are that the transverse shear the deformations of the shell are neglected (γxz = γθz = 0); in addition, also the rotary inertia of the shell is neglected.The dimensionless middle surface strains of the shellε˜x,0,ε˜θ,0,γ˜xθ,0 are expressed as [54](2)ε˜x,0=α∂u˜∂η+12α2∂w˜∂η2+18α∂v˜∂η−∂u˜∂θ2,(3)ε˜θ,0=∂v˜∂θ+w˜+12∂w˜∂θ−v˜2+18∂u˜∂θ−α∂v˜∂η2,(4)γ˜xθ,0=∂u˜∂θ+α∂v˜∂η+α∂w˜∂η∂w˜∂θ−v˜,where the linear terms are inserted into square brackets, the nonlinear terms are inserted into curly brackets, η = x/L is the dimensionless longitudinal coordinate of the shell and α = R/L.The separation of the linear and nonlinear terms in equations (2)–(4) can be useful in order to properly insert the strain contributions into the equation of the dimensionless elastic strain energy of the shell in the following linear and nonlinear vibration analyses.The dimensionless middle surface changes in curvature and torsion of the shellk˜x,k˜θ,k˜xθ are linear terms given by [54](5)k˜x=−α2∂2w˜∂η2,k˜θ=∂v˜∂θ−∂2w˜∂θ2,k˜xθ=−2α∂2w˜∂η∂θ+32α∂v˜∂η−12∂u˜∂θ.


## 2.2. Force and Moment Resultants

The dimensionless force (N˜x,N˜θ,N˜xθ,Q˜x,Q˜θ) and moment (M˜x,M˜θ,M˜xθ) resultants per unit length of the shell are expressed as [37](6)N˜x=11−ν2ε˜x,0+νε˜θ,0,N˜θ=11−ν2ε˜θ,0+νε˜x,0,N˜xθ=121+νγ˜xθ,0,(7)M˜x=β121−ν2k˜x+νk˜θ,M˜θ=β121−ν2k˜θ+νk˜x,M˜xθ=β241+νk˜xθ,(8)Q˜x=β2121−ν2αk˜x,x+νk˜θ,x+1−ν2k˜xθ,θ,Q˜θ=β2121−ν21−ν2αk˜xθ,x+k˜θ,θ+νk˜x,θ,where ν is Poisson’s ratio of the shell and β = h/R. These expressions of the dimensionless force and moment resultants will be used into the equations of the natural-type boundary conditions.


## 2.3. Elastic Strain Energy

In the case of homogeneous and isotropic material, the dimensionless elastic strain energyU˜ of a cylindrical shell, by neglecting the transverse normal stress σz (Love’s first approximation), can be written as a function of the dimensionless strains ε˜x,0,ε˜θ,0,γ˜xθ,0 and k˜x,k˜θ,k˜xθ in the form [37](9)U˜η,θ,τ=1211−ν2∫01∫02πε˜x,02+ε˜θ,02+2νε˜x,0ε˜θ,0+1−ν2γ˜xθ,02dηdθ+β212∫01∫02πk˜x2+k˜θ2+2νk˜xk˜θ+1−ν2k˜xθ2dηdθ,where the first term of the right-hand side of equation (9) is the membrane energy (also referred to stretching energy) and the second one is the bending energy.Another relevant issue of the thin-shell modelling of SWCNTs is given by the choice of an isotropic or anisotropic model; indeed, CNTs are discrete systems, i.e., they are intrinsically nonisotropic. However, it was proven in the past that the use of isotropic models for SWCNTs does not lead to significant errors. Therefore, in the present work, an isotropic elastic shell model is considered.


## 2.4. Kinetic Energy

The dimensionless time variableτ is defined introducing a reference dimensional circular frequency ω0 (i.e., the lowest extensional circular frequency of a ring in plane strain) in the form [56](10)τ=ω0t,ω0=E1−ν2ρR2,where t is the dimensional time and E and ρ are Young’s modulus and mass density of the shell, respectively.The dimensionless velocity fieldu˜′,v˜′,w˜′ of the shell is given by [56](11)u˜′=du˜dτ=u˙Rω0,v˜′=dv˜dτ=v˙Rω0,w˜′=dw˜dτ=w˙Rω0,where u˙,v˙,w˙ is the dimensional velocity field.The dimensionless kinetic energyT˜ of a cylindrical shell (rotary inertia is neglected) is [56](12)T˜η,θ,τ=12γ∫01∫02πu˜′2+v˜′2+w˜′2dηdθ,where γ=ρR2ω02/E.


## 2.5. Total Energy

The dimensionless total energyẼ of the shell can be expressed as the sum of the dimensionless elastic strain energy Ũ and kinetic energy T˜ in the following form:(13)E˜η,θ,τ=U˜η,θ,τ+T˜η,θ,τ.


## 3. Vibration Analysis

In order to study the linear and nonlinear vibrations of the shell, a two-step energy based procedure is adopted: (i) in the linear analysis, the three displacements are expanded by using a double mixed series, the potential and kinetic energies are developed in terms of the series free parameters, and the Rayleigh–Ritz method is considered to obtain approximated eigenfunctions; (ii) in the nonlinear analysis, the three displacements are re‐expanded by using the approximated eigenfunctions of the linear analysis, the potential and kinetic energies are developed in terms of modal coordinates and the Lagrange equations are applied to obtain a system of nonlinear ordinary differential equations of motion, which is solved numerically.


### 3.1. Linear Vibration Analysis

The linear vibration analysis is carried out by considering only the quadratic terms in the equation of the dimensionless elastic strain energy (9) and adopting the same procedure of [37, 45].A modal vibration, i.e., a synchronous motion, can be formally written in the form [37](14)u˜η,θ,τ=U˜η,θf˜τ,v˜η,θ,τ=V˜η,θf˜τ,w˜η,θ,τ=W˜η,θf˜τ,where u˜η,θ,τ, v˜η,θ,τ, w˜η,θ,τ are the three dimensionless displacements of the cylindrical shell; U˜η,θ, V˜η,θ, W˜η,θ describe the corresponding dimensionless mode shape; and f˜τ is the dimensionless time law, which is supposed to be the same for each displacement in the modal vibration analysis (synchronous motion hypothesis).The dimensionless mode shapeU˜,V˜,W˜ is expanded by considering a double mixed series in terms of m-th order dimensionless Chebyshev polynomials T˜m∗η along the longitudinal direction and harmonic functions cosnθ,sinnθ along the circumferential direction; following [37], the expansion reads(15)U˜η,θ=∑m=0Mu∑n=0NU˜m,nT˜m∗ηcosnθ,(16)V˜η,θ=∑m=0Mv∑n=0NV˜m,nT˜m∗ηsinnθ,(17)W˜η,θ=∑m=0Mw∑n=0NW˜m,nT˜m∗ηcosnθ,where T˜m∗η=T˜m2η−1, m denotes the degree of the Chebyshev polynomial and n is the number of nodal diameters.It must be noted thatU˜m,n,V˜m,n,W˜m,n are unknown dimensionless coefficients which are found by imposing the boundary conditions, i.e., by applying constraints to the expansions (15)–(17).Due to the axial symmetry and the isotropy of the system, in the absence of imperfections and in the case of axisymmetric boundary conditions (e.g., simply, clamped, and free), the harmonic functions are orthogonal with respect to the linear operator of the shell in the circumferential directionθ and the mode shapes are characterized by a specific number of nodal diameters n; therefore, the double series of equations (15)–(17) can be reduced to a single series and the linear mode shape can assume the following simplified expression [50]:(18)U˜η,θ=∑m=0MuU˜m,nT˜m∗ηcosnθ,Vη,θ=∑m=0MvV˜m,nTm∗ηsinnθ,Wη,θ=∑m=0MwW˜m,nTm∗ηcosnθ,n=0,1,2,…,N.The dimensionless mode shapeU˜c,V˜c,W˜c of the corresponding conjugate mode is given by(19)U˜cη,θ=∑m=0Mu∑n=1NU˜m,nT˜m∗ηsinnθ,(20)V˜cη,θ=−∑m=0Mv∑n=1NV˜m,nT˜m∗ηcosnθ,(21)W˜cη,θ=∑m=0Mw∑n=1NW˜m,nT˜m∗ηsinnθ.From expansions (19)–(21), it can be noted that the mode shape of the conjugate mode is rotated by an angle equal to π/(2n) with respect to the mode shape described in expansions (15)–(17). Moreover, expansions (19)–(21) do not include the axisymmetric modes (n = 0), which are taken into account in the previous ones (15)–(17).


### 3.2. Boundary Conditions

Free-free boundary conditions are given by the following equations [37](22)N˜x=0,N˜xθ+M˜xθR=0,Q˜x+1R∂M˜xθ∂θ=0,M˜x=0,forη=0,1,where the expressions of the dimensionless force N˜x,N˜xθ,Q˜x and moment M˜x,M˜xθ resultants per unit length are given in equations (6)–(8).It can be observed from equations (22) that the free-free boundary conditions are of natural type, i.e., they involve forces and moments only. Therefore, it is not strictly necessary that the expansions (15)–(17) respect such constraints; indeed, it is well known that the Rayleigh–Ritz method can be properly applied to obtain approximated eigenfunctions even if the natural boundary conditions are not respected.In the present paper, focused on free-free SWCNTs, no boundary equations are imposed and no constraints are applied to the unknown coefficients of expansions (15)–(17).


### 3.3. Rayleigh–Ritz Method

The maximum number of variables needed for describing a vibration mode withn nodal diameters is Np = Mu + Mv + Mw + 3 − p, where Mu = Mv = Mw denotes the degree of the Chebyshev polynomials and p describes the number of equations needed to satisfy the boundary conditions.In this paper, since no boundary equations are imposed for the free-free SWCNTs,p is taken equal to zero. Moreover, a specific convergence analysis is carried out in order to select the degree of the Chebyshev polynomials: degree 11 is found suitably accurate (Mu = Mv = Mw = 11). Therefore, in the present paper, it is imposed Np = 36.For a multi‐mode analysis including different values of nodal diametersn, the number of degrees of freedom of the system is computed by the relation Nmax = Np × (N + 1), where N represents the maximum value of the nodal diameters n considered.Equations (14) are inserted into the expressions of elastic strain energy (9) and kinetic energy (12) to compute the Rayleigh quotient Rq˜, where q˜ is a vector containing the unknown coefficients of expansions (15)–(17) [45](23)q˜=⋮U˜m,mV˜m,nW˜m,n⋮.After imposing the stationarity to the Rayleigh quotient, the eigenvalue problem is obtained [45]:(24)−ω2M˜+K˜q˜=0,which furnishes approximate natural frequencies (eigenvalues ωj) and mode shapes (eigenvectors q˜j), with j = (1, 2, …, Nmax).The approximated mode shape of thej-th mode is given by the equations (15)–(17) where coefficients U˜m,n,V˜m.n,W˜m,n are substituted with U˜m,nj,V˜m,nj,W˜m,nj, which denote the components of the j-th eigenvector q˜j of the equation (24).The approximation of thej-th eigenfunction vector of the original problem is given by [45](25)W˜jη,θ=U˜jη,θV˜jη,θW˜jη,θ.Thej-th eigenfunction vector (25) can be normalized, as shown, for example, in [45], imposing the following:(26)maxmaxUjη,θ,maxVjη,θ,maxWjη,θ=1.The previous normalization has the following simple physical meaning. Suppose that the dominant direction of vibration of thej-th mode shape is along the thickness h of the shell (radial direction w); then, the third component of the j-th eigenfunction vector (25) has a maximum amplitude higher than the other two components. Therefore, if the j-th eigenfunction vector is normalised in order to obtain a maximum amplitude equal to unity, then the dimensionless time law (modal coordinate) associated to the third component of the j-th mode shape gives the maximum amplitude of vibration along the radial direction of the cylindrical shell. This normalization results to be particularly useful in the nonlinear vibration analysis, where the nonlinearity appears when the amplitude of the radial vibration is of the order of the thickness, i.e., the dimensionless amplitude is of the order of 1; in such situation, all Lagrangian variables are of the order of 1 or less, and this strongly improves the numerical accuracy and efficiency.


### 3.4. Nonlinear Vibration Analysis

In the nonlinear vibration analysis, the full expression of the elastic strain energy (9) (quartic terms leading to cubic nonlinearity) is applied and the same procedure of Refs. [54, 56] is followed.The three dimensionless displacementsu˜η,θ,τ, v˜η,θ,τ, w˜η,θ,τ are re‐expanded by using the linear mode shapes U˜j,nη,θ, V˜j,nη,θ, W˜j,nη,θ, and their conjugates U˜cj,nη,θ, V˜cj,nη,θ, W˜cj,nη,θ, which are now available from the modal analysis of the previous section(27)u˜η,θ,τ=∑j=1Nu∑n=1NU˜j,nη,θf˜u,j,nt+U˜cj,nη,θf˜u,j,n,ct,(28)v˜η,θ,τ=∑j=1Nv∑n=1NV˜j,nη,θf˜v,j,nτ+V˜cj,nη,θf˜v,j,n,cτ,(29)w˜η,θ,τ=∑j=1Nw∑n=1NW˜j,nη,θf˜w,j,nτ+W˜cj,nη,θf˜w,j,n,cτ,where the dimensionless time laws f˜u,j,nτ,f˜u,j,n,cτ,f˜v,j,nτ,f˜v,j,n,cτ,f˜w,j,nτ,f˜w,j,n,cτ are unknown functions, the index j is used for ordering the modes with increasing natural frequency and the index n indicates the number of nodal diameters of the mode shape of the shell.From expansions (27)–(29), it can be observed that, in the nonlinear analysis, the time synchronicity is relaxed: for each mode j and for each component u˜,v˜,w˜, different time laws f˜τ are allowed.


### 3.5. Lagrange Equations

The expansions (27)–(29) are inserted into the equations of dimensionless elastic strain energy (9) and kinetic energy (12); the dimensionless Lagrange equations of motion, in the case of free vibrations, can be expressed in the form(30)ddτ∂T˜∂q˜i′+∂U˜∂q˜i=0,i∈1,Nmax,where the maximum number of degrees of freedom Nmax depends on the number of vibration modes considered in the expansions (27)–(29).The dimensionless Lagrangian coordinatesq˜i,q˜i′,q˜i″ can be written as(31)q˜i=qiR,q˜i′=dq˜idτ=q˙iω0R,q˜i″=dq˜i′dτ=q¨iω02R,i∈1,Nmax,where q˜i correspond to the previous dimensionless modal coordinates (f˜u,j,nτ, f˜u,j,n,cτ, f˜v,j,nτ, f˜v,j,n,cτ, f˜w,j,nτ,f˜w,j,n,cτ), and qi,q˙i,q¨i are the dimensional Lagrangian coordinates.By substituting the dimensionless vector functionsF˜q˜i=∂U˜/∂q˜i and M˜q˜i″=d/∂T˜/∂q˜i′/dτ in equation (30), where M˜ is the dimensionless mass matrix, it can be obtained that (see [56])(32)M˜q˜i″+F˜q˜i=0,i∈1,Nmax.Introducing the vector functionF˜x,i=M˜−1F˜q˜i in the equation (32), the equations of motion for free vibrations can be expressed in the following form:(33)q˜i″+F˜x,i=0,i∈1,Nmax.The set of nonlinear ordinary differential equations (33), completed by the modal initial conditions on displacements and velocities, is then solved numerically by means of the implicit Runge–Kutta iterative method with proper accuracy, precision, and number of steps.


## 3.1. Linear Vibration Analysis

The linear vibration analysis is carried out by considering only the quadratic terms in the equation of the dimensionless elastic strain energy (9) and adopting the same procedure of [37, 45].A modal vibration, i.e., a synchronous motion, can be formally written in the form [37](14)u˜η,θ,τ=U˜η,θf˜τ,v˜η,θ,τ=V˜η,θf˜τ,w˜η,θ,τ=W˜η,θf˜τ,where u˜η,θ,τ, v˜η,θ,τ, w˜η,θ,τ are the three dimensionless displacements of the cylindrical shell; U˜η,θ, V˜η,θ, W˜η,θ describe the corresponding dimensionless mode shape; and f˜τ is the dimensionless time law, which is supposed to be the same for each displacement in the modal vibration analysis (synchronous motion hypothesis).The dimensionless mode shapeU˜,V˜,W˜ is expanded by considering a double mixed series in terms of m-th order dimensionless Chebyshev polynomials T˜m∗η along the longitudinal direction and harmonic functions cosnθ,sinnθ along the circumferential direction; following [37], the expansion reads(15)U˜η,θ=∑m=0Mu∑n=0NU˜m,nT˜m∗ηcosnθ,(16)V˜η,θ=∑m=0Mv∑n=0NV˜m,nT˜m∗ηsinnθ,(17)W˜η,θ=∑m=0Mw∑n=0NW˜m,nT˜m∗ηcosnθ,where T˜m∗η=T˜m2η−1, m denotes the degree of the Chebyshev polynomial and n is the number of nodal diameters.It must be noted thatU˜m,n,V˜m,n,W˜m,n are unknown dimensionless coefficients which are found by imposing the boundary conditions, i.e., by applying constraints to the expansions (15)–(17).Due to the axial symmetry and the isotropy of the system, in the absence of imperfections and in the case of axisymmetric boundary conditions (e.g., simply, clamped, and free), the harmonic functions are orthogonal with respect to the linear operator of the shell in the circumferential directionθ and the mode shapes are characterized by a specific number of nodal diameters n; therefore, the double series of equations (15)–(17) can be reduced to a single series and the linear mode shape can assume the following simplified expression [50]:(18)U˜η,θ=∑m=0MuU˜m,nT˜m∗ηcosnθ,Vη,θ=∑m=0MvV˜m,nTm∗ηsinnθ,Wη,θ=∑m=0MwW˜m,nTm∗ηcosnθ,n=0,1,2,…,N.The dimensionless mode shapeU˜c,V˜c,W˜c of the corresponding conjugate mode is given by(19)U˜cη,θ=∑m=0Mu∑n=1NU˜m,nT˜m∗ηsinnθ,(20)V˜cη,θ=−∑m=0Mv∑n=1NV˜m,nT˜m∗ηcosnθ,(21)W˜cη,θ=∑m=0Mw∑n=1NW˜m,nT˜m∗ηsinnθ.From expansions (19)–(21), it can be noted that the mode shape of the conjugate mode is rotated by an angle equal to π/(2n) with respect to the mode shape described in expansions (15)–(17). Moreover, expansions (19)–(21) do not include the axisymmetric modes (n = 0), which are taken into account in the previous ones (15)–(17).


## 3.2. Boundary Conditions

Free-free boundary conditions are given by the following equations [37](22)N˜x=0,N˜xθ+M˜xθR=0,Q˜x+1R∂M˜xθ∂θ=0,M˜x=0,forη=0,1,where the expressions of the dimensionless force N˜x,N˜xθ,Q˜x and moment M˜x,M˜xθ resultants per unit length are given in equations (6)–(8).It can be observed from equations (22) that the free-free boundary conditions are of natural type, i.e., they involve forces and moments only. Therefore, it is not strictly necessary that the expansions (15)–(17) respect such constraints; indeed, it is well known that the Rayleigh–Ritz method can be properly applied to obtain approximated eigenfunctions even if the natural boundary conditions are not respected.In the present paper, focused on free-free SWCNTs, no boundary equations are imposed and no constraints are applied to the unknown coefficients of expansions (15)–(17).


## 3.3. Rayleigh–Ritz Method

The maximum number of variables needed for describing a vibration mode withn nodal diameters is Np = Mu + Mv + Mw + 3 − p, where Mu = Mv = Mw denotes the degree of the Chebyshev polynomials and p describes the number of equations needed to satisfy the boundary conditions.In this paper, since no boundary equations are imposed for the free-free SWCNTs,p is taken equal to zero. Moreover, a specific convergence analysis is carried out in order to select the degree of the Chebyshev polynomials: degree 11 is found suitably accurate (Mu = Mv = Mw = 11). Therefore, in the present paper, it is imposed Np = 36.For a multi‐mode analysis including different values of nodal diametersn, the number of degrees of freedom of the system is computed by the relation Nmax = Np × (N + 1), where N represents the maximum value of the nodal diameters n considered.Equations (14) are inserted into the expressions of elastic strain energy (9) and kinetic energy (12) to compute the Rayleigh quotient Rq˜, where q˜ is a vector containing the unknown coefficients of expansions (15)–(17) [45](23)q˜=⋮U˜m,mV˜m,nW˜m,n⋮.After imposing the stationarity to the Rayleigh quotient, the eigenvalue problem is obtained [45]:(24)−ω2M˜+K˜q˜=0,which furnishes approximate natural frequencies (eigenvalues ωj) and mode shapes (eigenvectors q˜j), with j = (1, 2, …, Nmax).The approximated mode shape of thej-th mode is given by the equations (15)–(17) where coefficients U˜m,n,V˜m.n,W˜m,n are substituted with U˜m,nj,V˜m,nj,W˜m,nj, which denote the components of the j-th eigenvector q˜j of the equation (24).The approximation of thej-th eigenfunction vector of the original problem is given by [45](25)W˜jη,θ=U˜jη,θV˜jη,θW˜jη,θ.Thej-th eigenfunction vector (25) can be normalized, as shown, for example, in [45], imposing the following:(26)maxmaxUjη,θ,maxVjη,θ,maxWjη,θ=1.The previous normalization has the following simple physical meaning. Suppose that the dominant direction of vibration of thej-th mode shape is along the thickness h of the shell (radial direction w); then, the third component of the j-th eigenfunction vector (25) has a maximum amplitude higher than the other two components. Therefore, if the j-th eigenfunction vector is normalised in order to obtain a maximum amplitude equal to unity, then the dimensionless time law (modal coordinate) associated to the third component of the j-th mode shape gives the maximum amplitude of vibration along the radial direction of the cylindrical shell. This normalization results to be particularly useful in the nonlinear vibration analysis, where the nonlinearity appears when the amplitude of the radial vibration is of the order of the thickness, i.e., the dimensionless amplitude is of the order of 1; in such situation, all Lagrangian variables are of the order of 1 or less, and this strongly improves the numerical accuracy and efficiency.


## 3.4. Nonlinear Vibration Analysis

In the nonlinear vibration analysis, the full expression of the elastic strain energy (9) (quartic terms leading to cubic nonlinearity) is applied and the same procedure of Refs. [54, 56] is followed.The three dimensionless displacementsu˜η,θ,τ, v˜η,θ,τ, w˜η,θ,τ are re‐expanded by using the linear mode shapes U˜j,nη,θ, V˜j,nη,θ, W˜j,nη,θ, and their conjugates U˜cj,nη,θ, V˜cj,nη,θ, W˜cj,nη,θ, which are now available from the modal analysis of the previous section(27)u˜η,θ,τ=∑j=1Nu∑n=1NU˜j,nη,θf˜u,j,nt+U˜cj,nη,θf˜u,j,n,ct,(28)v˜η,θ,τ=∑j=1Nv∑n=1NV˜j,nη,θf˜v,j,nτ+V˜cj,nη,θf˜v,j,n,cτ,(29)w˜η,θ,τ=∑j=1Nw∑n=1NW˜j,nη,θf˜w,j,nτ+W˜cj,nη,θf˜w,j,n,cτ,where the dimensionless time laws f˜u,j,nτ,f˜u,j,n,cτ,f˜v,j,nτ,f˜v,j,n,cτ,f˜w,j,nτ,f˜w,j,n,cτ are unknown functions, the index j is used for ordering the modes with increasing natural frequency and the index n indicates the number of nodal diameters of the mode shape of the shell.From expansions (27)–(29), it can be observed that, in the nonlinear analysis, the time synchronicity is relaxed: for each mode j and for each component u˜,v˜,w˜, different time laws f˜τ are allowed.


## 3.5. Lagrange Equations

The expansions (27)–(29) are inserted into the equations of dimensionless elastic strain energy (9) and kinetic energy (12); the dimensionless Lagrange equations of motion, in the case of free vibrations, can be expressed in the form(30)ddτ∂T˜∂q˜i′+∂U˜∂q˜i=0,i∈1,Nmax,where the maximum number of degrees of freedom Nmax depends on the number of vibration modes considered in the expansions (27)–(29).The dimensionless Lagrangian coordinatesq˜i,q˜i′,q˜i″ can be written as(31)q˜i=qiR,q˜i′=dq˜idτ=q˙iω0R,q˜i″=dq˜i′dτ=q¨iω02R,i∈1,Nmax,where q˜i correspond to the previous dimensionless modal coordinates (f˜u,j,nτ, f˜u,j,n,cτ, f˜v,j,nτ, f˜v,j,n,cτ, f˜w,j,nτ,f˜w,j,n,cτ), and qi,q˙i,q¨i are the dimensional Lagrangian coordinates.By substituting the dimensionless vector functionsF˜q˜i=∂U˜/∂q˜i and M˜q˜i″=d/∂T˜/∂q˜i′/dτ in equation (30), where M˜ is the dimensionless mass matrix, it can be obtained that (see [56])(32)M˜q˜i″+F˜q˜i=0,i∈1,Nmax.Introducing the vector functionF˜x,i=M˜−1F˜q˜i in the equation (32), the equations of motion for free vibrations can be expressed in the following form:(33)q˜i″+F˜x,i=0,i∈1,Nmax.The set of nonlinear ordinary differential equations (33), completed by the modal initial conditions on displacements and velocities, is then solved numerically by means of the implicit Runge–Kutta iterative method with proper accuracy, precision, and number of steps.


## 4. Numerical Results

The numerical analyses carried out in this section are obtained considering the equivalent mechanical parameters listed in Table1 (see [40]).Table 1
Equivalent mechanical parameters of the free-free SWCNT [40].
Young’s modulusE (TPa)
5.5
Poisson’s ratioν
0.19
Mass densityρ (kg/m3)
11700
Thicknessh (nm)
0.066In Tables2 and 3, the vibrational approach adopted in this paper, based on the Sanders-Koiter thin shell theory and the equivalent mechanical parameters reported in Table 1, is validated in the linear field (natural frequencies) by means of comparisons with experimental studies and numerical simulations from the pertinent literature. From these comparisons, it can be observed that the present continuous model is in very good agreement with the results of RRS [11] and MDS [17].Table 2
Natural frequencies of the radial breathing mode (m = 0, n = 0) of the free-free SWCNT of Table 1 for different values of radius R and same value of aspect ratio L/D. Comparisons between Sanders-Koiter elastic shell theory (SKT) (present model) and experimental resonant Raman spectroscopy (RRS) [11].
Radius
Aspect ratio
Natural frequency (THz)
Difference (%)
R (nm)
L/D
SKT (present model)
RRS [11]
0.509
10
6.905
7.165
3.63
0.518
10
6.785
7.105
4.50
0.527
10
6.669
6.865
2.85
0.569
10
6.177
6.295
1.87
0.705
10
5.025
5.276
4.76
0.708
10
4.964
5.216
4.83
0.718
10
4.895
5.066
3.37
0.733
10
4.788
4.947
3.21
0.746
10
4.711
4.917
4.19
0.764
10
4.594
4.797
4.23
0.770
10
4.559
4.737
3.76
0.782
10
4.494
4.677
3.91
0.799
10
4.393
4.617
4.85
0.814
10
4.318
4.527
4.62
0.822
10
4.271
4.437
3.74
0.847
10
4.150
4.317
4.04Table 3
Natural frequencies of the beam-like mode (m = 0, n = 1) of the clamped‐free SWCNT of Table 1 for same value of radius R and different values of aspect ratio L/D. Comparisons between Sanders‐Koiter elastic shell theory (SKT) (present model) and numerical molecular dynamics simulations (MDS) [17].
Radius
Aspect ratio
Natural frequency (THz)
Difference (%)
R (nm)
L/D
SKT (present model)
MDS [17]
0.339
5.26
0.217
0.212
2.36
0.339
5.62
0.191
0.188
1.60
0.339
5.99
0.169
0.167
1.20
0.339
6.35
0.151
0.150
0.67
0.339
6.71
0.136
0.136
≤0.01
0.339
7.07
0.123
0.123
≤0.01
0.339
7.44
0.112
0.112
≤0.01
0.339
7.80
0.102
0.102
≤0.01
0.339
8.16
0.094
0.094
≤0.01
0.339
8.52
0.086
0.086
≤0.01
0.339
8.89
0.080
0.080
≤0.01
0.339
9.25
0.074
0.074
≤0.01
0.339
9.61
0.069
0.069
≤0.01
0.339
9.98
0.064
0.064
≤0.01
0.339
10.34
0.060
0.060
≤0.01
0.339
10.70
0.056
0.056
≤0.01
0.339
11.06
0.053
0.053
≤0.01It must be underlined that this validation in linear field confirms the correctness of the continuous shell model proposed in the present paper with regard to the choice of the values of the equivalent mechanical parameters and to the previously assumed hypotheses (no size effects, isotropic shell).In Table4 and Figures 2–4, natural frequencies and mode shapes of the free-free SWCNT of Table 1 with radius R = 0.39 nm and length L = 3.0 nm obtained by applying the Sanders-Koiter thin shell theory are reported.Table 4
Natural frequencies of the free-free SWCNT of Table1 with radius R = 0.39 nm and length L = 3.0 nm.
Vibration mode (m,n)
Natural frequency (THz)
(0,2)
1.17609
(1,2)
1.21558
(2,2)
1.52195
(3,2)
2.32386
(0,4)
6.37264
(1,4)
6.42757
(2,4)
6.61777
(3,4)
6.95070
(0,0)
9.04466
(1,0)
9.28746
(2,0)
8.60968
(3,0)
8.79534Figure 2
Mode shapes of the free-free SWCNT of Table1 with radius R = 0.39 nm and length L = 3.0 nm. (a) Mode (m,n) = (0,2). (b) Mode (m,n) = (1,2). (c) Mode (m,n) = (2,2). (d) Mode (m,n) = (3,2).
(a)
(b)
(c)
(d)Figure 3
Mode shapes of the free-free SWCNT of Table1 with radius R = 0.39 nm and length L = 3.0 nm. (a) Mode (m,n) = (0,4). (b) Mode (m,n) = (1,4). (c) Mode (m,n) = (2,4). (d) Mode (m,n) = (3,4).
(a)
(b)
(c)
(d)Figure 4
Mode shapes of the free-free SWCNT of Table1 with radius R = 0.39 nm and length L = 3.0 nm. (a) Rayleigh inextensional mode (m,n) = (0,0). (b) Love inextensional mode (m,n) = (1,0) [38]. (c) Mode (m,n) = (2,0). (d) Mode (m,n) = (3,0).
(a)
(b)
(c)
(d)In the following, a further index is used for defining a mode shape and its longitudinal feature, i.e., the indexm, which indicates the number of nodal circumferences. It should be reminded that the mode shapes with m = 0 (no longitudinal half-waves) and m = 1 (one-half longitudinal half-wave) of Figure 4 correspond to the Rayleigh and the Love inextensional modes, respectively [38].In Figures5 and 6, mode shapes of four different pairs of conjugate modes of the free-free SWCNT of Table 1 having radius R = 0.39 nm and length L = 3.0 nm are reported; the conjugate modes have the same natural frequency and mode shape angularly shifted by π/2n. In particular, the conjugate modes (m,n) = (1,2), (m,n) = (1,2,c) with natural frequency f = 1.21558 THz and the conjugate modes (m,n) = (3,2), (m,n) = (3,2,c) with natural frequency f = 2.32386 THz have the respective mode shapes rotated by π/4, while the conjugate modes (m,n) = (1,4), (m,n) = (1,4,c) with natural frequency f = 6.42757 THz and the conjugate modes (m,n) = (3,4), (m,n) = (3,4,c) with natural frequency f = 6.9507 THz have the respective mode shapes rotated by π/8.Figure 5
Mode shape comparisons of the free-free SWCNT of Table1 with radius R = 0.39 nm and length L = 3.0 nm. (a) Conjugate modes (m,n) = (1,2), (b) (m,n) = (1,2,c): natural frequency f = 1.21558 THz, mode shapes rotated by π/4, (c) conjugate modes (m,n) = (3,2), and (d) (m,n) = (3,2,c): natural frequency f = 2.32386 THz, mode shapes rotated by π/4.
(a)
(b)
(c)
(d)Figure 6
Mode shape comparisons of the free-free SWCNT of Table1 with radius R = 0.39 nm and length L = 3.0 nm. (a) Conjugate modes (m,n) = (1,4), (b) (m,n) = (1,4,c): natural frequency f = 6.42757 THz, mode shapes rotated by π/8, (c) conjugate modes (m,n) = (3,4), and (d) (m,n) = (3,4,c): natural frequency f = 6.95070 THz, mode shapes rotated by π/8.
(a)
(b)
(c)
(d)


### 4.1. Linear Analysis

In order to better understand the nonlinear effect of the resonance interaction between the conjugate CFMs considered, an initial linear study is carried out by taking into account only the linear terms of the dimensionless middle surface strains (2)–(4) into the expression of the dimensionless elastic strain energy (9).The three displacementsu˜η,θ,τ, v˜η,θ,τ, w˜η,θ,τ are expanded considering the approximated linear mode shapes U˜j,nη,θ, U˜j,n,cη,θ, V˜j,nη,θ, V˜j,n,cη,θ, W˜j,nη,θ, W˜j,n,cη,θ of the two conjugate modes (1,2), (1,2,c) and the following linear modal expansion is obtained:(34)u˜η,θ,τ=U˜1,2η,θf˜u,1,2τ+U˜1,2,cη,θf˜u,1,2,cτ,v˜η,θ,τ=V˜1,2η,θf˜v,1,2τ+V˜1,2,cη,θf˜v,1,2,cτ,w˜η,θ,τ=W˜1,2η,θf˜w,1,2τ+W˜1,2,cη,θf˜w,1,2,cτ,where f˜u,1,2τ, f˜u,1,2,cτ, f˜v,1,2τ, f˜v,1,2,cτ, f˜w,1,2τ, f˜w,1,2,cτ are unknown linear modal coordinates and the maximum number of degrees of freedom is Nmax = 6.The modal initial conditionsf˜u,1,20, f˜u,1,2,c0, f˜v,1,20, f˜v,1,2,c0, f˜w,1,20, f˜w,1,2,c0 imposed on the linear modal coordinates are given by(35)f˜u,1,20=0.1,f˜v,1,20=0.1,f˜w,1,20=0.1,f˜u,1,2,c0=0.05,f˜v,1,2,c0=0.05,f˜w,1,2,c0=0.05.It should be highlighted that, in the present analysis, the initial velocities are taken equal to zero (the modal initial conditions are imposed only on the displacements).The system of linear ordinary differential equations of motion (33), completed with the initial conditions (35), is then solved through an implicit Runge–Kutta numerical method in order to obtain time histories, frequency spectra, and phase portraits describing the linear steady-state dynamic response of the SWCNT.


#### 4.1.1. Time Histories, Frequency Spectra, and Phase Portraits

The natural frequency and dimensionless time period of the two conjugate modes (1,2), (1,2,c) are f(1,2) = f(1,2,c) = 1.21558 THz and τ(1,2) = τ(1,2,c) = 2 × π × f(0,2)/f(1,2) = 6.07905, respectively, where f(0,2) = 1.17609 THz is the natural frequency of the reference (i.e., lowest frequency) mode (0,2).The time histories of the two conjugate modes (1,2), (1,2,c) are reported in Figure 7(a); this figure clearly proves that the vibration is a periodic motion, i.e., the CNT is vibrating in linear conditions.Figure 7
Time histories (a) and frequency spectra (b) of the two conjugate modes (1,2), (1,2,c). Free-free SWCNT of Table 1. Modal expansion (34). Initial conditions (35). Time period τ(1,2) = 6.07905.
(a)
(b)The corresponding frequency spectra are shown in Figure7(b), where the horizontal axis represents the normalised circular frequency ω/ω(0,2). The spectra are perfectly clean, with a single spike at the frequency corresponding to the eigenfrequency of the selected mode, as it is expected in linear field.It must be stressed that these analyses were performed in order to compare the results of the linear vibrations with the following results of the nonlinear vibrations.The two-dimensional phase portrait of the mode (1,2) in the plane displacement-velocity is shown in Figure8(a). By calculating the first (lowest) root of the equation f˜w,1,2τ=f˜w,1,20=0.1, the time τ = 6.07905 is found, which is equal to the time period τ(1,2) (the motion is periodic).Figure 8
Phase portraits of the conjugate modes (1,2) (a) and (1,2,c) (b). Free-free SWCNT of Table 1. Modal expansion (34). Initial conditions (35). Time period τ(1,2) = 6.07905.
(a)
(b)Similarly, the two-dimensional phase portrait of the conjugate vibration mode (1,2,c) in the plane displacement-velocity is reported in Figure 8(b). Calculating the first (lowest) root of the equation f˜w,1,2,cτ=f˜w,1,2,c0=0.05, the same time τ = 6.07905 of the mode (1,2) is obtained (the conjugate modes have the same time period).


#### 4.1.2. Total Energy Distribution

The evolution in time of the total energy distributionẼ (η, τ) of the conjugate modes (1,2), (1,2,c) in the linear field along the SWCNT axis, by assuming the circumferential coordinate θ = 0, is shown in Figure 9. An initial energy is imposed on the edge sections along the SWCNT axis (η = (0,1)): it can be noted that this energy is preserved in amplitude throughout the time.Figure 9
Contour plot (a) and 3D plot (b) of the total energy distributionẼ (η, τ) of the conjugate CFMs (1,2), (1,2,c) along the SWCNT axis in time. Circumferential coordinate θ = 0. Free-free SWCNT of Table 1.
(a)
(b)The total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the CNT surface (η, θ) for different time values throughout the time period τ(1,2) is shown in Figures 10 and 11. From these figures, it can be observed that the total energy distribution over the carbon nanotube surface is periodic along the circumferential direction θ and symmetric along the longitudinal direction η.Figure 10
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ) for different time values throughout the reference period τ(1,2). (a) τ = 0, (b) τ = τ(1,2)/11, (c) τ = 2 × τ(1,2)/11, (d) τ = 3 × τ(1,2)/11, (e) τ = 4 × τ(1,2)/11, (f) τ = 5 × τ(1,2)/11, (g) τ = 6 × τ(1,2)/11, (h) τ = 7 × τ(1,2)/11, (i) τ = 8 × τ(1,2)/11, (j) τ = 9 × τ(1,2)/11, (k) τ = 10 × τ(1,2)/11, (l) τ = τ(1,2).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Figure 11
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ) for different time values throughout the reference period τ(1,2). (a) τ = 0, (b) τ = τ(1,2)/11, (c) τ = 2 × τ(1,2)/11, (d) τ = 3 × τ(1,2)/11, (e) τ = 4 × τ(1,2)/11, (f) τ = 5 × τ(1,2)/11, (g) τ = 6 × τ(1,2)/11, (h) τ = 7 × τ(1,2)/11, (i) τ = 8 × τ(1,2)/11, (j) τ = 9 × τ(1,2)/11, (k) τ = 10 × τ(1,2)/11, and (l) τ = τ(1,2).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)


#### 4.1.3. Total Energy Conservation

After integration along the longitudinalη and circumferential θ coordinates, the evolution in time of the elastic strain and kinetic energy of the conjugate CFMs (1,2), (1,2,c) is found, see Figure 12(a). The evolution in time of the total energy is shown in Figure 12(b): since no damping is included in this model, then the dynamical system is conservative over the time.Figure 12
Evolution in time of the energy of the conjugate modes (1,2), (1,2,c). (a) Comparison of the two energies. (b) Total energy conservation.
(a)
(b)In the following, the previous analyses carried out for the conjugate CFMs (1,2), (1,2,c) in the linear field will be extended to the nonlinear field: nonlinear time histories, frequency spectra, and phase portraits will be found. Moreover, the nonlinear total energy distribution over the CNT surface in the time will be obtained. The aim of these analyses is the investigation of the phenomena arising in presence of nonlinear resonance interaction and energy exchange between the two conjugate CFMs.


### 4.2. Nonlinear Analysis

In the nonlinear vibration analysis, the full expression of the dimensionless elastic strain energy (9), containing terms up to the fourth order (cubic nonlinearity), is considered.


#### 4.2.1. Convergence Analysis

The first step of the nonlinear study is the convergence analysis in terms of the modal expansion (27)–(29): it is carried out by selecting the vibration modes to be inserted into the modal expansion providing an accurate description of the actual nonlinear behaviour of the SWCNT.In particular, the aim of this convergence analysis is to find a modal expansion with the minimum number of degrees of freedom that allows the nonlinear behaviour of the sum of the time histories of the two conjugate modes (1,2), (1,2,c) to be correctly described.An initial two-mode approximation involving only the conjugate modes (1,2) and (1,2,c) is applied (6 dof model); the convergence is then checked by adding suitable modes to the conjugate ones, i.e., asymmetric and axisymmetric modes, in the modal expansion (27)–(29), as shown in Table 5.Table 5
Nonlinear convergence analysis of the sum of the time histories of the two conjugate modes (1,2), (1,2,c). Nonlinear modal expansion (27)–(29). Modal initial conditions (36). 22 dof reference model. ERRORRMS% in the time domain for the 6 dof model, 8 dof model, 14 dof model, and 20 dof model. Free-free SWCNT of Table 1.
(j,n)
(1,2)
(1,2,c)
(1,4)
(1,4,c)
(3,2)
(3,2,c)
(1,0)
(3,0)
ERRORRMS%
6 dof model
u
,
v
,
w
u
,
v
,
w
—
—
—
—
—
—
12.8
8 dof model
u
,
v
,
w
u
,
v
,
w
—
—
—
—
u
,
w
—
8.63
14 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
—
—
u
,
w
—
4.54
20 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
w
—
2.23
22 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
w
u
,
w
—The modal initial conditionsf˜u,j,n0,f˜v,j,n0,f˜w,j,n0 applied in this convergence analysis are(36)f˜u,1,20=0.1,f˜v,1,20=0.1,f˜w,1,20=0.1,f˜u,1,2,c0=0.05,f˜v,1,2,c0=0.05,f˜w,1,2,c0=0.05,f˜u,j,n0=0.001,f˜v,j,n0=0.001,f˜w,j,n0=0.001.The convergence is reached by means of a 22 dof model, described in Table5, which is assumed as a reference model; the percentage root mean square error in the time domain (ERRORRMS%) with respect to the reference model is used as a comparison parameter for the other nonlinear models.In Figure13, the nonlinear convergence analysis of the sum of the time histories of the conjugate modes (1,2), (1,2,c) obtained applying the nonlinear modal expansion (27)–(29) with the modal initial conditions (36) for the different nonlinear models of Table 5 is shown. The 22 dof nonlinear model is assumed as a reference model.Figure 13
Nonlinear convergence analysis of the sum of the time histories of the two conjugate modes (1,2), (1,2,c). Nonlinear modal expansion (27)–(29). Modal initial conditions (36). Nonlinear models of Table 5. 22 dof reference model. (a) 6 dof model vs. 22 dof model, (b) 8 dof model vs. 22 dof model, (c) 14 dof model vs. 22 dof model, (d) 20 dof model vs 22 dof model, (e) Comparisons, and (f) Comparisons.
(a)
(b)
(c)
(d)
(e)
(f)From this convergence analysis, it can be found that the smallest model able to predict the nonlinear dynamics of the system with an acceptable accuracy is the 14 dof model (ERRORRMS% = 4 ÷ 5). The main weakness of the 6–8 dof models is the insufficient number of asymmetric and axisymmetric modes, which are very important for properly modelling bending deformation (asymmetric modes) and circumferential stretching (axisymmetric modes) during the modal vibration, see [51, 52].Since the 14 dof model provides accurate results with the minimal computational effort, then, in the following, this model will be used for studying the nonlinear oscillations of the two conjugate modes (1,2), (1,2,c). By inserting the following into modal expansion (27)–(29), the resulting nonlinear modal expansion is written in the form (37):(i)
modes (1,2), (1,2,c), (1,4), (1,4,c), (1,0) along the longitudinal displacement u(ii)
modes (1,2), (1,2,c), (1,4), (1,4,c) along the circumferential displacement v(iii)
modes (1,2), (1,2,c), (1,4), (1,4,c), (1,0) along the radial displacement w(37)u˜η,θ,τ=U˜1,2η,θf˜u,1,2τ+U˜1,2,cη,θf˜u,1,2,cτ+U˜1,4η,θf˜u,1,4τ+U˜1,4,cη,θf˜u,1,4,cτ+U˜1,0η,θf˜u,1,0τ,v˜η,θ,τ=V˜1,2η,θf˜v,1,2τ+V˜1,2,cη,θf˜v,1,2,cτ+V˜1,4η,θf˜v,1,4τ+V˜1,4,cη,θf˜v,1,4,cτ,w˜η,θ,τ=W˜1,2η,θf˜w,1,2τ+W˜1,2,cη,θf˜w,1,2,cτ+W˜1,4η,θf˜w,1,4τ+W˜1,4,cη,θf˜w,1,4,cτ+W˜1,0η,θf˜w,1,0τ,where the longitudinal and radial displacements have five terms, the circumferential displacement has four terms, and the total number of degrees of freedom of the model is equal to (Nmax = 14).


#### 4.2.2. Modal Initial Conditions

In this section, the minimum value of the modal initial conditions (36) which activates the nonlinear response of the conjugate mode (1,2,c) is investigated. The nonlinear modal expansion (37) is used (14 dof model). The initial energy on the modes (1,2), (1,4), (1,4,c), (1,0) is fixed. The value of the initial energy imposed on the mode (1,2,c) is increased until its nonlinear activation, see Table 6.Table 6
Modal initial conditions imposed on the conjugate modes (1,2), (1,2,c).
Case
Modal initial conditions on the modes (1,2), (1,2,c)
Modal initial condition ratio
A
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
∞
B
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.00005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
2000
C
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.0001
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
1000
D
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.0005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
200
E
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.001
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
100
F
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
20In Figure14, the time histories of the conjugate modes (1,2), (1,2,c) are shown. In Figures 14(a)–14(c) (cases A, B, C of Table 6), the mode (1,2,c) is not activated. In Figures 14(d)–14(e) (cases D, E of Table 6), a weak activation of the mode (1,2,c) arises after suitable long time. In Figure 14(f) (case F of Table 6), the nonlinear response of the mode (1,2,c) is strongly activated and it vibrates with a large amplitude: this is caused by an energy transfer probably due to a nonlinear resonance interaction between the two conjugate CFMs.Figure 14
Nonlinear time histories of the conjugate modes (1,2), (1,2,c) corresponding to the different modal initial conditions listed in Table 6. Nonlinear modal expansion (37). 14 dof model. (a) Case A: f˜w,1,20=0.1, f˜w,1,2,c0=0, (b) Case B: f˜w,1,20=0.1, f˜w,1,2,c0=0.00005, (c) Case C: f˜w,1,20=0.1, f˜w,1,2,c0=0.0001, (d) Case D: f˜w,1,20=0.1, f˜w,1,2,c0=0.0005, (e) Case E: f˜w,1,20=0.1, f˜w,1,2,c0=0.001, (f) Case F: f˜w,1,20=0.1, f˜w,1,2,c0=0.005.
(a)
(b)
(c)
(d)
(e)
(f)From Table6 and Figure 14, it is obtained that the modal initial conditions activating the nonlinear response of the conjugate mode (1,2,c) are given by (Case F)(38)f˜u,1,20=0.1,f˜v,1,20=0.1,f˜w,1,20=0.1,f˜u,1,2,c0=0.005,f˜v,1,2,c0=0.005,f˜w,1,2,c0=0.005,f˜u,j,n0=0.001,f˜v,j,n0=0.001,f˜w,j,n0=0.001.


#### 4.2.3. Time Histories, Frequency Spectra, and Phase Portraits

The time histories of the conjugate modes (1,2), (1,2,c) of the CNT are shown in Figure 15(a), the corresponding frequency spectra are presented in Figure 15(b). These plots clearly demonstrate that the vibration is a chaotic motion (no time period), i.e., the CNT is vibrating in nonlinear conditions, where the resonance interaction gives rise to strong energy exchange between the conjugate CFMs.Figure 15
Time histories (a) and frequency spectra (b) of the two conjugate modes (1,2), (1,2,c). Free-free SWCNT of Table 1. Nonlinear modal expansion (37). Modal initial conditions (38).
(a)
(b)Phase portraits of the mode (1,2) for different times are shown in Figure16. By calculating the first (lowest) root of the equation f˜w,1,2τ=f˜w,1,20=0.1, the time τ = 0.77834 is obtained, which is not a time period.Figure 16
2D and 3D phase portraits of the mode (1,2). Free-free SWCNT of Table1. Nonlinear modal expansion (37). Modal initial conditions (38). (a) Time τ = 5. (b) Time τ = 20. (c) Time τ = 50.
(a)
(b)
(c)Similarly, phase portraits of the conjugate mode (1,2,c) for different times are reported in Figure 17. Calculating the first (lowest) root of the equation f˜w,1,2,cτ=f˜w,1,2,c0=0.005, the time τ = 0.78436 is found, which is different from the time value previously obtained for the mode (1,2).Figure 17
2D and 3D phase portraits of the conjugate mode (1,2,c). Free-free SWCNT of Table 1. Nonlinear modal expansion (37). Modal initial conditions (38). (a) Time τ = 5. (b) Time τ = 20. (c) Time τ = 50.
(a)
(b)
(c)From Figures16 and 17 (phase portraits), it is confirmed the non‐periodicity of the motion of the CNT in nonlinear field (i.e., chaotic motion) observed in Figure 15 (time histories and frequency spectra).


#### 4.2.4. Total Energy Distribution

The evolution in time of the total energy distributionẼ (η, τ) of the conjugate modes (1,2), (1,2,c) along the SWCNT axis, considering the circumferential coordinate θ = 0, is shown in Figure 18. An initial energy is imposed on the edge sections of the SWCNT (η = (0,1)): differently from the linear case (Figure 9), this energy is not preserved in amplitude throughout the time.Figure 18
Contour plot (a) and 3D plot (b) of the total energy distributionẼ (η, τ) of the conjugate CFMs (1,2), (1,2,c) along the SWCNT axis in time. Circumferential coordinate θ = 0. Free-free SWCNT of Table 1.
(a)
(b)


#### 4.2.5. Travelling Wave Response

In the present section, the energy exchange between the conjugate CFMs (1,2), (1,2,c) is studied. By considering the time histories of the two modes in the time range τ = (80, 90) (Figure 19(a)), it can be observed the presence of the phase difference Δϕ = π between the two nonlinear responses in the specific time range τ = (84.27, 85.17), see the enlargement of Figure 19(b).Figure 19
Time histories of the two conjugate modes (1,2), (1,2,c). (a) Time range τ = (80, 90). (b) Time range τ = (84.27, 85.17).
(a)
(b)The evolution of the total energy distributionẼ (η, θ, τ) of the two conjugate CFMs (1,2), (1,2,c) over the CNT surface (η, θ) in the time range τ = (84.27, 85.17) is reported in Figures 20–25: differently from the linear case (Figures 10 and 11), the total energy distribution over the SWCNT surface evolves in the time in a complex pattern, where the total energy symmetry along the longitudinal direction and periodicity along the circumferential direction are still preserved.Figure 20
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57). (a) τ = 84.27, (b) τ = 84.33, (c) τ = 84.39, (d) τ = 84.45, (e) τ = 84.51, (f) τ = 84.57.
(a)
(b)
(c)
(d)
(e)
(f)Figure 21
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87). (a) τ = 84.57, (b) τ = 84.63, (c) τ = 84.69, (d) τ = 84.75, (e) τ = 84.81, and (f) τ = 84.87.
(a)
(b)
(c)
(d)
(e)
(f)Figure 22
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17). (a) τ = 84.87, (b) τ = 84.93, (c) τ = 84.99, (d) τ = 85.05, (e) τ = 85.11, and (f) τ = 85.17.
(a)
(b)
(c)
(d)
(e)
(f)Figure 23
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57). (a) τ = 84.27, (b) τ = 84.33, (c) τ = 84.39, (d) τ = 84.45, (e) τ = 84.51, and (f) τ = 84.57.
(a)
(b)
(c)
(d)
(e)
(f)Figure 24
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87). (a) τ = 84.57, (b) τ = 84.63, (c) τ = 84.69, (d) τ = 84.75, (e) τ = 84.81, and (f) τ = 84.87.
(a)
(b)
(c)
(d)
(e)
(f)Figure 25
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17). (a) τ = 84.87, (b) τ = 84.93, (c) τ = 84.99, (d) τ = 85.05, (e) τ = 85.11, and (f) τ = 85.17.
(a)
(b)
(c)
(d)
(e)
(f)From Figures20–25, it can be observed that the nonlinear resonance interaction between the two conjugate modes, in presence of the phase difference Δϕ = π, gives rise to an energy transfer between the two modes which generates a pure travelling wave moving circumferentially around the CNT. In Figures 20 and 23, the forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57), corresponding to the initial ascending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is presented. In Figures 21 and 24, the backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87), corresponding to the intermediate descending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is described. In Figures 22 and 25, the forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17), corresponding to the final ascending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is shown.It is reported in the literature for the circular cylindrical shells that the presence of one couple of modes having the same mode shape but different angular orientation (such as the conjugate modes (1,2), (1,2,c) of the present paper), the first one described by cos nθ (see radial displacement (17), directly excited mode), the second one described by sin nθ (see radial displacement (21), slightly perturbed mode), in case of nonlinear resonance interaction between the modes, can lead to the appearance of a pure travelling wave moving along the circumferential direction of the shell. This phenomenon, which was observed also experimentally, is due to the axial symmetry of the structure and it occurs only in nonlinear field (that is a relevant difference vs. linear vibrations), see [58, 59] for more details. To the authors’ best knowledge, the travelling wave phenomenon was not still observed for the SWCNTs, which are modelled in the present paper as continuous circular cylindrical shells in the framework of a thin shell theory.


## 4.1. Linear Analysis

In order to better understand the nonlinear effect of the resonance interaction between the conjugate CFMs considered, an initial linear study is carried out by taking into account only the linear terms of the dimensionless middle surface strains (2)–(4) into the expression of the dimensionless elastic strain energy (9).The three displacementsu˜η,θ,τ, v˜η,θ,τ, w˜η,θ,τ are expanded considering the approximated linear mode shapes U˜j,nη,θ, U˜j,n,cη,θ, V˜j,nη,θ, V˜j,n,cη,θ, W˜j,nη,θ, W˜j,n,cη,θ of the two conjugate modes (1,2), (1,2,c) and the following linear modal expansion is obtained:(34)u˜η,θ,τ=U˜1,2η,θf˜u,1,2τ+U˜1,2,cη,θf˜u,1,2,cτ,v˜η,θ,τ=V˜1,2η,θf˜v,1,2τ+V˜1,2,cη,θf˜v,1,2,cτ,w˜η,θ,τ=W˜1,2η,θf˜w,1,2τ+W˜1,2,cη,θf˜w,1,2,cτ,where f˜u,1,2τ, f˜u,1,2,cτ, f˜v,1,2τ, f˜v,1,2,cτ, f˜w,1,2τ, f˜w,1,2,cτ are unknown linear modal coordinates and the maximum number of degrees of freedom is Nmax = 6.The modal initial conditionsf˜u,1,20, f˜u,1,2,c0, f˜v,1,20, f˜v,1,2,c0, f˜w,1,20, f˜w,1,2,c0 imposed on the linear modal coordinates are given by(35)f˜u,1,20=0.1,f˜v,1,20=0.1,f˜w,1,20=0.1,f˜u,1,2,c0=0.05,f˜v,1,2,c0=0.05,f˜w,1,2,c0=0.05.It should be highlighted that, in the present analysis, the initial velocities are taken equal to zero (the modal initial conditions are imposed only on the displacements).The system of linear ordinary differential equations of motion (33), completed with the initial conditions (35), is then solved through an implicit Runge–Kutta numerical method in order to obtain time histories, frequency spectra, and phase portraits describing the linear steady-state dynamic response of the SWCNT.


### 4.1.1. Time Histories, Frequency Spectra, and Phase Portraits

The natural frequency and dimensionless time period of the two conjugate modes (1,2), (1,2,c) are f(1,2) = f(1,2,c) = 1.21558 THz and τ(1,2) = τ(1,2,c) = 2 × π × f(0,2)/f(1,2) = 6.07905, respectively, where f(0,2) = 1.17609 THz is the natural frequency of the reference (i.e., lowest frequency) mode (0,2).The time histories of the two conjugate modes (1,2), (1,2,c) are reported in Figure 7(a); this figure clearly proves that the vibration is a periodic motion, i.e., the CNT is vibrating in linear conditions.Figure 7
Time histories (a) and frequency spectra (b) of the two conjugate modes (1,2), (1,2,c). Free-free SWCNT of Table 1. Modal expansion (34). Initial conditions (35). Time period τ(1,2) = 6.07905.
(a)
(b)The corresponding frequency spectra are shown in Figure7(b), where the horizontal axis represents the normalised circular frequency ω/ω(0,2). The spectra are perfectly clean, with a single spike at the frequency corresponding to the eigenfrequency of the selected mode, as it is expected in linear field.It must be stressed that these analyses were performed in order to compare the results of the linear vibrations with the following results of the nonlinear vibrations.The two-dimensional phase portrait of the mode (1,2) in the plane displacement-velocity is shown in Figure8(a). By calculating the first (lowest) root of the equation f˜w,1,2τ=f˜w,1,20=0.1, the time τ = 6.07905 is found, which is equal to the time period τ(1,2) (the motion is periodic).Figure 8
Phase portraits of the conjugate modes (1,2) (a) and (1,2,c) (b). Free-free SWCNT of Table 1. Modal expansion (34). Initial conditions (35). Time period τ(1,2) = 6.07905.
(a)
(b)Similarly, the two-dimensional phase portrait of the conjugate vibration mode (1,2,c) in the plane displacement-velocity is reported in Figure 8(b). Calculating the first (lowest) root of the equation f˜w,1,2,cτ=f˜w,1,2,c0=0.05, the same time τ = 6.07905 of the mode (1,2) is obtained (the conjugate modes have the same time period).


### 4.1.2. Total Energy Distribution

The evolution in time of the total energy distributionẼ (η, τ) of the conjugate modes (1,2), (1,2,c) in the linear field along the SWCNT axis, by assuming the circumferential coordinate θ = 0, is shown in Figure 9. An initial energy is imposed on the edge sections along the SWCNT axis (η = (0,1)): it can be noted that this energy is preserved in amplitude throughout the time.Figure 9
Contour plot (a) and 3D plot (b) of the total energy distributionẼ (η, τ) of the conjugate CFMs (1,2), (1,2,c) along the SWCNT axis in time. Circumferential coordinate θ = 0. Free-free SWCNT of Table 1.
(a)
(b)The total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the CNT surface (η, θ) for different time values throughout the time period τ(1,2) is shown in Figures 10 and 11. From these figures, it can be observed that the total energy distribution over the carbon nanotube surface is periodic along the circumferential direction θ and symmetric along the longitudinal direction η.Figure 10
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ) for different time values throughout the reference period τ(1,2). (a) τ = 0, (b) τ = τ(1,2)/11, (c) τ = 2 × τ(1,2)/11, (d) τ = 3 × τ(1,2)/11, (e) τ = 4 × τ(1,2)/11, (f) τ = 5 × τ(1,2)/11, (g) τ = 6 × τ(1,2)/11, (h) τ = 7 × τ(1,2)/11, (i) τ = 8 × τ(1,2)/11, (j) τ = 9 × τ(1,2)/11, (k) τ = 10 × τ(1,2)/11, (l) τ = τ(1,2).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Figure 11
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ) for different time values throughout the reference period τ(1,2). (a) τ = 0, (b) τ = τ(1,2)/11, (c) τ = 2 × τ(1,2)/11, (d) τ = 3 × τ(1,2)/11, (e) τ = 4 × τ(1,2)/11, (f) τ = 5 × τ(1,2)/11, (g) τ = 6 × τ(1,2)/11, (h) τ = 7 × τ(1,2)/11, (i) τ = 8 × τ(1,2)/11, (j) τ = 9 × τ(1,2)/11, (k) τ = 10 × τ(1,2)/11, and (l) τ = τ(1,2).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)


### 4.1.3. Total Energy Conservation

After integration along the longitudinalη and circumferential θ coordinates, the evolution in time of the elastic strain and kinetic energy of the conjugate CFMs (1,2), (1,2,c) is found, see Figure 12(a). The evolution in time of the total energy is shown in Figure 12(b): since no damping is included in this model, then the dynamical system is conservative over the time.Figure 12
Evolution in time of the energy of the conjugate modes (1,2), (1,2,c). (a) Comparison of the two energies. (b) Total energy conservation.
(a)
(b)In the following, the previous analyses carried out for the conjugate CFMs (1,2), (1,2,c) in the linear field will be extended to the nonlinear field: nonlinear time histories, frequency spectra, and phase portraits will be found. Moreover, the nonlinear total energy distribution over the CNT surface in the time will be obtained. The aim of these analyses is the investigation of the phenomena arising in presence of nonlinear resonance interaction and energy exchange between the two conjugate CFMs.


## 4.1.1. Time Histories, Frequency Spectra, and Phase Portraits

The natural frequency and dimensionless time period of the two conjugate modes (1,2), (1,2,c) are f(1,2) = f(1,2,c) = 1.21558 THz and τ(1,2) = τ(1,2,c) = 2 × π × f(0,2)/f(1,2) = 6.07905, respectively, where f(0,2) = 1.17609 THz is the natural frequency of the reference (i.e., lowest frequency) mode (0,2).The time histories of the two conjugate modes (1,2), (1,2,c) are reported in Figure 7(a); this figure clearly proves that the vibration is a periodic motion, i.e., the CNT is vibrating in linear conditions.Figure 7
Time histories (a) and frequency spectra (b) of the two conjugate modes (1,2), (1,2,c). Free-free SWCNT of Table 1. Modal expansion (34). Initial conditions (35). Time period τ(1,2) = 6.07905.
(a)
(b)The corresponding frequency spectra are shown in Figure7(b), where the horizontal axis represents the normalised circular frequency ω/ω(0,2). The spectra are perfectly clean, with a single spike at the frequency corresponding to the eigenfrequency of the selected mode, as it is expected in linear field.It must be stressed that these analyses were performed in order to compare the results of the linear vibrations with the following results of the nonlinear vibrations.The two-dimensional phase portrait of the mode (1,2) in the plane displacement-velocity is shown in Figure8(a). By calculating the first (lowest) root of the equation f˜w,1,2τ=f˜w,1,20=0.1, the time τ = 6.07905 is found, which is equal to the time period τ(1,2) (the motion is periodic).Figure 8
Phase portraits of the conjugate modes (1,2) (a) and (1,2,c) (b). Free-free SWCNT of Table 1. Modal expansion (34). Initial conditions (35). Time period τ(1,2) = 6.07905.
(a)
(b)Similarly, the two-dimensional phase portrait of the conjugate vibration mode (1,2,c) in the plane displacement-velocity is reported in Figure 8(b). Calculating the first (lowest) root of the equation f˜w,1,2,cτ=f˜w,1,2,c0=0.05, the same time τ = 6.07905 of the mode (1,2) is obtained (the conjugate modes have the same time period).


## 4.1.2. Total Energy Distribution

The evolution in time of the total energy distributionẼ (η, τ) of the conjugate modes (1,2), (1,2,c) in the linear field along the SWCNT axis, by assuming the circumferential coordinate θ = 0, is shown in Figure 9. An initial energy is imposed on the edge sections along the SWCNT axis (η = (0,1)): it can be noted that this energy is preserved in amplitude throughout the time.Figure 9
Contour plot (a) and 3D plot (b) of the total energy distributionẼ (η, τ) of the conjugate CFMs (1,2), (1,2,c) along the SWCNT axis in time. Circumferential coordinate θ = 0. Free-free SWCNT of Table 1.
(a)
(b)The total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the CNT surface (η, θ) for different time values throughout the time period τ(1,2) is shown in Figures 10 and 11. From these figures, it can be observed that the total energy distribution over the carbon nanotube surface is periodic along the circumferential direction θ and symmetric along the longitudinal direction η.Figure 10
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ) for different time values throughout the reference period τ(1,2). (a) τ = 0, (b) τ = τ(1,2)/11, (c) τ = 2 × τ(1,2)/11, (d) τ = 3 × τ(1,2)/11, (e) τ = 4 × τ(1,2)/11, (f) τ = 5 × τ(1,2)/11, (g) τ = 6 × τ(1,2)/11, (h) τ = 7 × τ(1,2)/11, (i) τ = 8 × τ(1,2)/11, (j) τ = 9 × τ(1,2)/11, (k) τ = 10 × τ(1,2)/11, (l) τ = τ(1,2).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Figure 11
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ) for different time values throughout the reference period τ(1,2). (a) τ = 0, (b) τ = τ(1,2)/11, (c) τ = 2 × τ(1,2)/11, (d) τ = 3 × τ(1,2)/11, (e) τ = 4 × τ(1,2)/11, (f) τ = 5 × τ(1,2)/11, (g) τ = 6 × τ(1,2)/11, (h) τ = 7 × τ(1,2)/11, (i) τ = 8 × τ(1,2)/11, (j) τ = 9 × τ(1,2)/11, (k) τ = 10 × τ(1,2)/11, and (l) τ = τ(1,2).
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)


## 4.1.3. Total Energy Conservation

After integration along the longitudinalη and circumferential θ coordinates, the evolution in time of the elastic strain and kinetic energy of the conjugate CFMs (1,2), (1,2,c) is found, see Figure 12(a). The evolution in time of the total energy is shown in Figure 12(b): since no damping is included in this model, then the dynamical system is conservative over the time.Figure 12
Evolution in time of the energy of the conjugate modes (1,2), (1,2,c). (a) Comparison of the two energies. (b) Total energy conservation.
(a)
(b)In the following, the previous analyses carried out for the conjugate CFMs (1,2), (1,2,c) in the linear field will be extended to the nonlinear field: nonlinear time histories, frequency spectra, and phase portraits will be found. Moreover, the nonlinear total energy distribution over the CNT surface in the time will be obtained. The aim of these analyses is the investigation of the phenomena arising in presence of nonlinear resonance interaction and energy exchange between the two conjugate CFMs.


## 4.2. Nonlinear Analysis

In the nonlinear vibration analysis, the full expression of the dimensionless elastic strain energy (9), containing terms up to the fourth order (cubic nonlinearity), is considered.


### 4.2.1. Convergence Analysis

The first step of the nonlinear study is the convergence analysis in terms of the modal expansion (27)–(29): it is carried out by selecting the vibration modes to be inserted into the modal expansion providing an accurate description of the actual nonlinear behaviour of the SWCNT.In particular, the aim of this convergence analysis is to find a modal expansion with the minimum number of degrees of freedom that allows the nonlinear behaviour of the sum of the time histories of the two conjugate modes (1,2), (1,2,c) to be correctly described.An initial two-mode approximation involving only the conjugate modes (1,2) and (1,2,c) is applied (6 dof model); the convergence is then checked by adding suitable modes to the conjugate ones, i.e., asymmetric and axisymmetric modes, in the modal expansion (27)–(29), as shown in Table 5.Table 5
Nonlinear convergence analysis of the sum of the time histories of the two conjugate modes (1,2), (1,2,c). Nonlinear modal expansion (27)–(29). Modal initial conditions (36). 22 dof reference model. ERRORRMS% in the time domain for the 6 dof model, 8 dof model, 14 dof model, and 20 dof model. Free-free SWCNT of Table 1.
(j,n)
(1,2)
(1,2,c)
(1,4)
(1,4,c)
(3,2)
(3,2,c)
(1,0)
(3,0)
ERRORRMS%
6 dof model
u
,
v
,
w
u
,
v
,
w
—
—
—
—
—
—
12.8
8 dof model
u
,
v
,
w
u
,
v
,
w
—
—
—
—
u
,
w
—
8.63
14 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
—
—
u
,
w
—
4.54
20 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
w
—
2.23
22 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
w
u
,
w
—The modal initial conditionsf˜u,j,n0,f˜v,j,n0,f˜w,j,n0 applied in this convergence analysis are(36)f˜u,1,20=0.1,f˜v,1,20=0.1,f˜w,1,20=0.1,f˜u,1,2,c0=0.05,f˜v,1,2,c0=0.05,f˜w,1,2,c0=0.05,f˜u,j,n0=0.001,f˜v,j,n0=0.001,f˜w,j,n0=0.001.The convergence is reached by means of a 22 dof model, described in Table5, which is assumed as a reference model; the percentage root mean square error in the time domain (ERRORRMS%) with respect to the reference model is used as a comparison parameter for the other nonlinear models.In Figure13, the nonlinear convergence analysis of the sum of the time histories of the conjugate modes (1,2), (1,2,c) obtained applying the nonlinear modal expansion (27)–(29) with the modal initial conditions (36) for the different nonlinear models of Table 5 is shown. The 22 dof nonlinear model is assumed as a reference model.Figure 13
Nonlinear convergence analysis of the sum of the time histories of the two conjugate modes (1,2), (1,2,c). Nonlinear modal expansion (27)–(29). Modal initial conditions (36). Nonlinear models of Table 5. 22 dof reference model. (a) 6 dof model vs. 22 dof model, (b) 8 dof model vs. 22 dof model, (c) 14 dof model vs. 22 dof model, (d) 20 dof model vs 22 dof model, (e) Comparisons, and (f) Comparisons.
(a)
(b)
(c)
(d)
(e)
(f)From this convergence analysis, it can be found that the smallest model able to predict the nonlinear dynamics of the system with an acceptable accuracy is the 14 dof model (ERRORRMS% = 4 ÷ 5). The main weakness of the 6–8 dof models is the insufficient number of asymmetric and axisymmetric modes, which are very important for properly modelling bending deformation (asymmetric modes) and circumferential stretching (axisymmetric modes) during the modal vibration, see [51, 52].Since the 14 dof model provides accurate results with the minimal computational effort, then, in the following, this model will be used for studying the nonlinear oscillations of the two conjugate modes (1,2), (1,2,c). By inserting the following into modal expansion (27)–(29), the resulting nonlinear modal expansion is written in the form (37):(i)
modes (1,2), (1,2,c), (1,4), (1,4,c), (1,0) along the longitudinal displacement u(ii)
modes (1,2), (1,2,c), (1,4), (1,4,c) along the circumferential displacement v(iii)
modes (1,2), (1,2,c), (1,4), (1,4,c), (1,0) along the radial displacement w(37)u˜η,θ,τ=U˜1,2η,θf˜u,1,2τ+U˜1,2,cη,θf˜u,1,2,cτ+U˜1,4η,θf˜u,1,4τ+U˜1,4,cη,θf˜u,1,4,cτ+U˜1,0η,θf˜u,1,0τ,v˜η,θ,τ=V˜1,2η,θf˜v,1,2τ+V˜1,2,cη,θf˜v,1,2,cτ+V˜1,4η,θf˜v,1,4τ+V˜1,4,cη,θf˜v,1,4,cτ,w˜η,θ,τ=W˜1,2η,θf˜w,1,2τ+W˜1,2,cη,θf˜w,1,2,cτ+W˜1,4η,θf˜w,1,4τ+W˜1,4,cη,θf˜w,1,4,cτ+W˜1,0η,θf˜w,1,0τ,where the longitudinal and radial displacements have five terms, the circumferential displacement has four terms, and the total number of degrees of freedom of the model is equal to (Nmax = 14).


### 4.2.2. Modal Initial Conditions

In this section, the minimum value of the modal initial conditions (36) which activates the nonlinear response of the conjugate mode (1,2,c) is investigated. The nonlinear modal expansion (37) is used (14 dof model). The initial energy on the modes (1,2), (1,4), (1,4,c), (1,0) is fixed. The value of the initial energy imposed on the mode (1,2,c) is increased until its nonlinear activation, see Table 6.Table 6
Modal initial conditions imposed on the conjugate modes (1,2), (1,2,c).
Case
Modal initial conditions on the modes (1,2), (1,2,c)
Modal initial condition ratio
A
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
∞
B
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.00005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
2000
C
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.0001
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
1000
D
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.0005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
200
E
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.001
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
100
F
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
20In Figure14, the time histories of the conjugate modes (1,2), (1,2,c) are shown. In Figures 14(a)–14(c) (cases A, B, C of Table 6), the mode (1,2,c) is not activated. In Figures 14(d)–14(e) (cases D, E of Table 6), a weak activation of the mode (1,2,c) arises after suitable long time. In Figure 14(f) (case F of Table 6), the nonlinear response of the mode (1,2,c) is strongly activated and it vibrates with a large amplitude: this is caused by an energy transfer probably due to a nonlinear resonance interaction between the two conjugate CFMs.Figure 14
Nonlinear time histories of the conjugate modes (1,2), (1,2,c) corresponding to the different modal initial conditions listed in Table 6. Nonlinear modal expansion (37). 14 dof model. (a) Case A: f˜w,1,20=0.1, f˜w,1,2,c0=0, (b) Case B: f˜w,1,20=0.1, f˜w,1,2,c0=0.00005, (c) Case C: f˜w,1,20=0.1, f˜w,1,2,c0=0.0001, (d) Case D: f˜w,1,20=0.1, f˜w,1,2,c0=0.0005, (e) Case E: f˜w,1,20=0.1, f˜w,1,2,c0=0.001, (f) Case F: f˜w,1,20=0.1, f˜w,1,2,c0=0.005.
(a)
(b)
(c)
(d)
(e)
(f)From Table6 and Figure 14, it is obtained that the modal initial conditions activating the nonlinear response of the conjugate mode (1,2,c) are given by (Case F)(38)f˜u,1,20=0.1,f˜v,1,20=0.1,f˜w,1,20=0.1,f˜u,1,2,c0=0.005,f˜v,1,2,c0=0.005,f˜w,1,2,c0=0.005,f˜u,j,n0=0.001,f˜v,j,n0=0.001,f˜w,j,n0=0.001.


### 4.2.3. Time Histories, Frequency Spectra, and Phase Portraits

The time histories of the conjugate modes (1,2), (1,2,c) of the CNT are shown in Figure 15(a), the corresponding frequency spectra are presented in Figure 15(b). These plots clearly demonstrate that the vibration is a chaotic motion (no time period), i.e., the CNT is vibrating in nonlinear conditions, where the resonance interaction gives rise to strong energy exchange between the conjugate CFMs.Figure 15
Time histories (a) and frequency spectra (b) of the two conjugate modes (1,2), (1,2,c). Free-free SWCNT of Table 1. Nonlinear modal expansion (37). Modal initial conditions (38).
(a)
(b)Phase portraits of the mode (1,2) for different times are shown in Figure16. By calculating the first (lowest) root of the equation f˜w,1,2τ=f˜w,1,20=0.1, the time τ = 0.77834 is obtained, which is not a time period.Figure 16
2D and 3D phase portraits of the mode (1,2). Free-free SWCNT of Table1. Nonlinear modal expansion (37). Modal initial conditions (38). (a) Time τ = 5. (b) Time τ = 20. (c) Time τ = 50.
(a)
(b)
(c)Similarly, phase portraits of the conjugate mode (1,2,c) for different times are reported in Figure 17. Calculating the first (lowest) root of the equation f˜w,1,2,cτ=f˜w,1,2,c0=0.005, the time τ = 0.78436 is found, which is different from the time value previously obtained for the mode (1,2).Figure 17
2D and 3D phase portraits of the conjugate mode (1,2,c). Free-free SWCNT of Table 1. Nonlinear modal expansion (37). Modal initial conditions (38). (a) Time τ = 5. (b) Time τ = 20. (c) Time τ = 50.
(a)
(b)
(c)From Figures16 and 17 (phase portraits), it is confirmed the non‐periodicity of the motion of the CNT in nonlinear field (i.e., chaotic motion) observed in Figure 15 (time histories and frequency spectra).


### 4.2.4. Total Energy Distribution

The evolution in time of the total energy distributionẼ (η, τ) of the conjugate modes (1,2), (1,2,c) along the SWCNT axis, considering the circumferential coordinate θ = 0, is shown in Figure 18. An initial energy is imposed on the edge sections of the SWCNT (η = (0,1)): differently from the linear case (Figure 9), this energy is not preserved in amplitude throughout the time.Figure 18
Contour plot (a) and 3D plot (b) of the total energy distributionẼ (η, τ) of the conjugate CFMs (1,2), (1,2,c) along the SWCNT axis in time. Circumferential coordinate θ = 0. Free-free SWCNT of Table 1.
(a)
(b)


### 4.2.5. Travelling Wave Response

In the present section, the energy exchange between the conjugate CFMs (1,2), (1,2,c) is studied. By considering the time histories of the two modes in the time range τ = (80, 90) (Figure 19(a)), it can be observed the presence of the phase difference Δϕ = π between the two nonlinear responses in the specific time range τ = (84.27, 85.17), see the enlargement of Figure 19(b).Figure 19
Time histories of the two conjugate modes (1,2), (1,2,c). (a) Time range τ = (80, 90). (b) Time range τ = (84.27, 85.17).
(a)
(b)The evolution of the total energy distributionẼ (η, θ, τ) of the two conjugate CFMs (1,2), (1,2,c) over the CNT surface (η, θ) in the time range τ = (84.27, 85.17) is reported in Figures 20–25: differently from the linear case (Figures 10 and 11), the total energy distribution over the SWCNT surface evolves in the time in a complex pattern, where the total energy symmetry along the longitudinal direction and periodicity along the circumferential direction are still preserved.Figure 20
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57). (a) τ = 84.27, (b) τ = 84.33, (c) τ = 84.39, (d) τ = 84.45, (e) τ = 84.51, (f) τ = 84.57.
(a)
(b)
(c)
(d)
(e)
(f)Figure 21
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87). (a) τ = 84.57, (b) τ = 84.63, (c) τ = 84.69, (d) τ = 84.75, (e) τ = 84.81, and (f) τ = 84.87.
(a)
(b)
(c)
(d)
(e)
(f)Figure 22
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17). (a) τ = 84.87, (b) τ = 84.93, (c) τ = 84.99, (d) τ = 85.05, (e) τ = 85.11, and (f) τ = 85.17.
(a)
(b)
(c)
(d)
(e)
(f)Figure 23
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57). (a) τ = 84.27, (b) τ = 84.33, (c) τ = 84.39, (d) τ = 84.45, (e) τ = 84.51, and (f) τ = 84.57.
(a)
(b)
(c)
(d)
(e)
(f)Figure 24
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87). (a) τ = 84.57, (b) τ = 84.63, (c) τ = 84.69, (d) τ = 84.75, (e) τ = 84.81, and (f) τ = 84.87.
(a)
(b)
(c)
(d)
(e)
(f)Figure 25
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17). (a) τ = 84.87, (b) τ = 84.93, (c) τ = 84.99, (d) τ = 85.05, (e) τ = 85.11, and (f) τ = 85.17.
(a)
(b)
(c)
(d)
(e)
(f)From Figures20–25, it can be observed that the nonlinear resonance interaction between the two conjugate modes, in presence of the phase difference Δϕ = π, gives rise to an energy transfer between the two modes which generates a pure travelling wave moving circumferentially around the CNT. In Figures 20 and 23, the forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57), corresponding to the initial ascending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is presented. In Figures 21 and 24, the backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87), corresponding to the intermediate descending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is described. In Figures 22 and 25, the forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17), corresponding to the final ascending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is shown.It is reported in the literature for the circular cylindrical shells that the presence of one couple of modes having the same mode shape but different angular orientation (such as the conjugate modes (1,2), (1,2,c) of the present paper), the first one described by cos nθ (see radial displacement (17), directly excited mode), the second one described by sin nθ (see radial displacement (21), slightly perturbed mode), in case of nonlinear resonance interaction between the modes, can lead to the appearance of a pure travelling wave moving along the circumferential direction of the shell. This phenomenon, which was observed also experimentally, is due to the axial symmetry of the structure and it occurs only in nonlinear field (that is a relevant difference vs. linear vibrations), see [58, 59] for more details. To the authors’ best knowledge, the travelling wave phenomenon was not still observed for the SWCNTs, which are modelled in the present paper as continuous circular cylindrical shells in the framework of a thin shell theory.


## 4.2.1. Convergence Analysis

The first step of the nonlinear study is the convergence analysis in terms of the modal expansion (27)–(29): it is carried out by selecting the vibration modes to be inserted into the modal expansion providing an accurate description of the actual nonlinear behaviour of the SWCNT.In particular, the aim of this convergence analysis is to find a modal expansion with the minimum number of degrees of freedom that allows the nonlinear behaviour of the sum of the time histories of the two conjugate modes (1,2), (1,2,c) to be correctly described.An initial two-mode approximation involving only the conjugate modes (1,2) and (1,2,c) is applied (6 dof model); the convergence is then checked by adding suitable modes to the conjugate ones, i.e., asymmetric and axisymmetric modes, in the modal expansion (27)–(29), as shown in Table 5.Table 5
Nonlinear convergence analysis of the sum of the time histories of the two conjugate modes (1,2), (1,2,c). Nonlinear modal expansion (27)–(29). Modal initial conditions (36). 22 dof reference model. ERRORRMS% in the time domain for the 6 dof model, 8 dof model, 14 dof model, and 20 dof model. Free-free SWCNT of Table 1.
(j,n)
(1,2)
(1,2,c)
(1,4)
(1,4,c)
(3,2)
(3,2,c)
(1,0)
(3,0)
ERRORRMS%
6 dof model
u
,
v
,
w
u
,
v
,
w
—
—
—
—
—
—
12.8
8 dof model
u
,
v
,
w
u
,
v
,
w
—
—
—
—
u
,
w
—
8.63
14 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
—
—
u
,
w
—
4.54
20 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
w
—
2.23
22 dof model
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
v
,
w
u
,
w
u
,
w
—The modal initial conditionsf˜u,j,n0,f˜v,j,n0,f˜w,j,n0 applied in this convergence analysis are(36)f˜u,1,20=0.1,f˜v,1,20=0.1,f˜w,1,20=0.1,f˜u,1,2,c0=0.05,f˜v,1,2,c0=0.05,f˜w,1,2,c0=0.05,f˜u,j,n0=0.001,f˜v,j,n0=0.001,f˜w,j,n0=0.001.The convergence is reached by means of a 22 dof model, described in Table5, which is assumed as a reference model; the percentage root mean square error in the time domain (ERRORRMS%) with respect to the reference model is used as a comparison parameter for the other nonlinear models.In Figure13, the nonlinear convergence analysis of the sum of the time histories of the conjugate modes (1,2), (1,2,c) obtained applying the nonlinear modal expansion (27)–(29) with the modal initial conditions (36) for the different nonlinear models of Table 5 is shown. The 22 dof nonlinear model is assumed as a reference model.Figure 13
Nonlinear convergence analysis of the sum of the time histories of the two conjugate modes (1,2), (1,2,c). Nonlinear modal expansion (27)–(29). Modal initial conditions (36). Nonlinear models of Table 5. 22 dof reference model. (a) 6 dof model vs. 22 dof model, (b) 8 dof model vs. 22 dof model, (c) 14 dof model vs. 22 dof model, (d) 20 dof model vs 22 dof model, (e) Comparisons, and (f) Comparisons.
(a)
(b)
(c)
(d)
(e)
(f)From this convergence analysis, it can be found that the smallest model able to predict the nonlinear dynamics of the system with an acceptable accuracy is the 14 dof model (ERRORRMS% = 4 ÷ 5). The main weakness of the 6–8 dof models is the insufficient number of asymmetric and axisymmetric modes, which are very important for properly modelling bending deformation (asymmetric modes) and circumferential stretching (axisymmetric modes) during the modal vibration, see [51, 52].Since the 14 dof model provides accurate results with the minimal computational effort, then, in the following, this model will be used for studying the nonlinear oscillations of the two conjugate modes (1,2), (1,2,c). By inserting the following into modal expansion (27)–(29), the resulting nonlinear modal expansion is written in the form (37):(i)
modes (1,2), (1,2,c), (1,4), (1,4,c), (1,0) along the longitudinal displacement u(ii)
modes (1,2), (1,2,c), (1,4), (1,4,c) along the circumferential displacement v(iii)
modes (1,2), (1,2,c), (1,4), (1,4,c), (1,0) along the radial displacement w(37)u˜η,θ,τ=U˜1,2η,θf˜u,1,2τ+U˜1,2,cη,θf˜u,1,2,cτ+U˜1,4η,θf˜u,1,4τ+U˜1,4,cη,θf˜u,1,4,cτ+U˜1,0η,θf˜u,1,0τ,v˜η,θ,τ=V˜1,2η,θf˜v,1,2τ+V˜1,2,cη,θf˜v,1,2,cτ+V˜1,4η,θf˜v,1,4τ+V˜1,4,cη,θf˜v,1,4,cτ,w˜η,θ,τ=W˜1,2η,θf˜w,1,2τ+W˜1,2,cη,θf˜w,1,2,cτ+W˜1,4η,θf˜w,1,4τ+W˜1,4,cη,θf˜w,1,4,cτ+W˜1,0η,θf˜w,1,0τ,where the longitudinal and radial displacements have five terms, the circumferential displacement has four terms, and the total number of degrees of freedom of the model is equal to (Nmax = 14).


## 4.2.2. Modal Initial Conditions

In this section, the minimum value of the modal initial conditions (36) which activates the nonlinear response of the conjugate mode (1,2,c) is investigated. The nonlinear modal expansion (37) is used (14 dof model). The initial energy on the modes (1,2), (1,4), (1,4,c), (1,0) is fixed. The value of the initial energy imposed on the mode (1,2,c) is increased until its nonlinear activation, see Table 6.Table 6
Modal initial conditions imposed on the conjugate modes (1,2), (1,2,c).
Case
Modal initial conditions on the modes (1,2), (1,2,c)
Modal initial condition ratio
A
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
∞
B
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.00005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
2000
C
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.0001
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
1000
D
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.0005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
200
E
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.001
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
100
F
f
˜
w
,
1,2
0
=
0.1
f
˜
w
,
1,2
,
c
0
=
0.005
f
˜
w
,
1,2
0
/
f
˜
w
,
1,2
,
c
0
=
20In Figure14, the time histories of the conjugate modes (1,2), (1,2,c) are shown. In Figures 14(a)–14(c) (cases A, B, C of Table 6), the mode (1,2,c) is not activated. In Figures 14(d)–14(e) (cases D, E of Table 6), a weak activation of the mode (1,2,c) arises after suitable long time. In Figure 14(f) (case F of Table 6), the nonlinear response of the mode (1,2,c) is strongly activated and it vibrates with a large amplitude: this is caused by an energy transfer probably due to a nonlinear resonance interaction between the two conjugate CFMs.Figure 14
Nonlinear time histories of the conjugate modes (1,2), (1,2,c) corresponding to the different modal initial conditions listed in Table 6. Nonlinear modal expansion (37). 14 dof model. (a) Case A: f˜w,1,20=0.1, f˜w,1,2,c0=0, (b) Case B: f˜w,1,20=0.1, f˜w,1,2,c0=0.00005, (c) Case C: f˜w,1,20=0.1, f˜w,1,2,c0=0.0001, (d) Case D: f˜w,1,20=0.1, f˜w,1,2,c0=0.0005, (e) Case E: f˜w,1,20=0.1, f˜w,1,2,c0=0.001, (f) Case F: f˜w,1,20=0.1, f˜w,1,2,c0=0.005.
(a)
(b)
(c)
(d)
(e)
(f)From Table6 and Figure 14, it is obtained that the modal initial conditions activating the nonlinear response of the conjugate mode (1,2,c) are given by (Case F)(38)f˜u,1,20=0.1,f˜v,1,20=0.1,f˜w,1,20=0.1,f˜u,1,2,c0=0.005,f˜v,1,2,c0=0.005,f˜w,1,2,c0=0.005,f˜u,j,n0=0.001,f˜v,j,n0=0.001,f˜w,j,n0=0.001.


## 4.2.3. Time Histories, Frequency Spectra, and Phase Portraits

The time histories of the conjugate modes (1,2), (1,2,c) of the CNT are shown in Figure 15(a), the corresponding frequency spectra are presented in Figure 15(b). These plots clearly demonstrate that the vibration is a chaotic motion (no time period), i.e., the CNT is vibrating in nonlinear conditions, where the resonance interaction gives rise to strong energy exchange between the conjugate CFMs.Figure 15
Time histories (a) and frequency spectra (b) of the two conjugate modes (1,2), (1,2,c). Free-free SWCNT of Table 1. Nonlinear modal expansion (37). Modal initial conditions (38).
(a)
(b)Phase portraits of the mode (1,2) for different times are shown in Figure16. By calculating the first (lowest) root of the equation f˜w,1,2τ=f˜w,1,20=0.1, the time τ = 0.77834 is obtained, which is not a time period.Figure 16
2D and 3D phase portraits of the mode (1,2). Free-free SWCNT of Table1. Nonlinear modal expansion (37). Modal initial conditions (38). (a) Time τ = 5. (b) Time τ = 20. (c) Time τ = 50.
(a)
(b)
(c)Similarly, phase portraits of the conjugate mode (1,2,c) for different times are reported in Figure 17. Calculating the first (lowest) root of the equation f˜w,1,2,cτ=f˜w,1,2,c0=0.005, the time τ = 0.78436 is found, which is different from the time value previously obtained for the mode (1,2).Figure 17
2D and 3D phase portraits of the conjugate mode (1,2,c). Free-free SWCNT of Table 1. Nonlinear modal expansion (37). Modal initial conditions (38). (a) Time τ = 5. (b) Time τ = 20. (c) Time τ = 50.
(a)
(b)
(c)From Figures16 and 17 (phase portraits), it is confirmed the non‐periodicity of the motion of the CNT in nonlinear field (i.e., chaotic motion) observed in Figure 15 (time histories and frequency spectra).


## 4.2.4. Total Energy Distribution

The evolution in time of the total energy distributionẼ (η, τ) of the conjugate modes (1,2), (1,2,c) along the SWCNT axis, considering the circumferential coordinate θ = 0, is shown in Figure 18. An initial energy is imposed on the edge sections of the SWCNT (η = (0,1)): differently from the linear case (Figure 9), this energy is not preserved in amplitude throughout the time.Figure 18
Contour plot (a) and 3D plot (b) of the total energy distributionẼ (η, τ) of the conjugate CFMs (1,2), (1,2,c) along the SWCNT axis in time. Circumferential coordinate θ = 0. Free-free SWCNT of Table 1.
(a)
(b)


## 4.2.5. Travelling Wave Response

In the present section, the energy exchange between the conjugate CFMs (1,2), (1,2,c) is studied. By considering the time histories of the two modes in the time range τ = (80, 90) (Figure 19(a)), it can be observed the presence of the phase difference Δϕ = π between the two nonlinear responses in the specific time range τ = (84.27, 85.17), see the enlargement of Figure 19(b).Figure 19
Time histories of the two conjugate modes (1,2), (1,2,c). (a) Time range τ = (80, 90). (b) Time range τ = (84.27, 85.17).
(a)
(b)The evolution of the total energy distributionẼ (η, θ, τ) of the two conjugate CFMs (1,2), (1,2,c) over the CNT surface (η, θ) in the time range τ = (84.27, 85.17) is reported in Figures 20–25: differently from the linear case (Figures 10 and 11), the total energy distribution over the SWCNT surface evolves in the time in a complex pattern, where the total energy symmetry along the longitudinal direction and periodicity along the circumferential direction are still preserved.Figure 20
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57). (a) τ = 84.27, (b) τ = 84.33, (c) τ = 84.39, (d) τ = 84.45, (e) τ = 84.51, (f) τ = 84.57.
(a)
(b)
(c)
(d)
(e)
(f)Figure 21
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87). (a) τ = 84.57, (b) τ = 84.63, (c) τ = 84.69, (d) τ = 84.75, (e) τ = 84.81, and (f) τ = 84.87.
(a)
(b)
(c)
(d)
(e)
(f)Figure 22
Contour plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17). (a) τ = 84.87, (b) τ = 84.93, (c) τ = 84.99, (d) τ = 85.05, (e) τ = 85.11, and (f) τ = 85.17.
(a)
(b)
(c)
(d)
(e)
(f)Figure 23
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57). (a) τ = 84.27, (b) τ = 84.33, (c) τ = 84.39, (d) τ = 84.45, (e) τ = 84.51, and (f) τ = 84.57.
(a)
(b)
(c)
(d)
(e)
(f)Figure 24
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87). (a) τ = 84.57, (b) τ = 84.63, (c) τ = 84.69, (d) τ = 84.75, (e) τ = 84.81, and (f) τ = 84.87.
(a)
(b)
(c)
(d)
(e)
(f)Figure 25
3D plot of the total energy distributionẼ (η, θ, τ) of the conjugate modes (1,2), (1,2,c) over the SWCNT surface (η, θ). Forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17). (a) τ = 84.87, (b) τ = 84.93, (c) τ = 84.99, (d) τ = 85.05, (e) τ = 85.11, and (f) τ = 85.17.
(a)
(b)
(c)
(d)
(e)
(f)From Figures20–25, it can be observed that the nonlinear resonance interaction between the two conjugate modes, in presence of the phase difference Δϕ = π, gives rise to an energy transfer between the two modes which generates a pure travelling wave moving circumferentially around the CNT. In Figures 20 and 23, the forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.27, 84.57), corresponding to the initial ascending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is presented. In Figures 21 and 24, the backward travelling wave from θ = 2π to θ = 0 in the time range τ = (84.57, 84.87), corresponding to the intermediate descending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is described. In Figures 22 and 25, the forward travelling wave from θ = 0 to θ = 2π in the time range τ = (84.87, 85.17), corresponding to the final ascending part of the nonlinear time history of the mode (1,2) of Figure 19(b), is shown.It is reported in the literature for the circular cylindrical shells that the presence of one couple of modes having the same mode shape but different angular orientation (such as the conjugate modes (1,2), (1,2,c) of the present paper), the first one described by cos nθ (see radial displacement (17), directly excited mode), the second one described by sin nθ (see radial displacement (21), slightly perturbed mode), in case of nonlinear resonance interaction between the modes, can lead to the appearance of a pure travelling wave moving along the circumferential direction of the shell. This phenomenon, which was observed also experimentally, is due to the axial symmetry of the structure and it occurs only in nonlinear field (that is a relevant difference vs. linear vibrations), see [58, 59] for more details. To the authors’ best knowledge, the travelling wave phenomenon was not still observed for the SWCNTs, which are modelled in the present paper as continuous circular cylindrical shells in the framework of a thin shell theory.


## 5. Concluding Remarks

In the present paper, the nonlinear resonance interaction between two conjugate CFMs of a SWCNT with free-free boundary conditions is analysed. The Sanders-Koiter shell theory is used to model the dynamics of the shell. The Rayleigh–Ritz method is applied to obtain approximated eigenfunctions. The Lagrange equations are used to get the system of nonlinear ordinary differential equations of motion, which is numerically solved by means of the implicit Runge–Kutta iterative method.(1)
Numerical results from the linear analysisThe model proposed in this paper is validated in linear field by comparing the natural frequencies obtained with data retrieved from the pertinent literature: from these comparisons, it can be observed that the present model is in very good agreement with the results of experimental analyses (RRS) and structural simulations (MDS).The time histories, frequency spectra, and phase portraits obtained demonstrate a periodic vibration of the SWCNT in linear field.The total energy distribution over the carbon nanotube surface is periodic along the circumferential direction and symmetric along the longitudinal direction, where the initial energy imposed on the edge sections of the SWCNT is preserved in amplitude throughout the time; the total energy conservation of the considered undamped system is verified.(2)
Numerical results from the nonlinear analysisA convergence analysis is carried out to select the vibration modes to be added to the two directly excited CFMs providing an accurate description of the actual nonlinear behaviour of the SWCNT.The modal expansion obtained involves, in addition to the conjugate modes (1,2) and (1,2,c), also the conjugate asymmetric modes (1,4), (1,4,c) and the axisymmetric mode (1,0) (14 dof model): it is confirmed the fundamental role of the additional asymmetric and axisymmetric modes to properly predict the nonlinear dynamics of the system.The modal initial conditions activating the nonlinear response of the mode (1,2,c) are investigated: this activation gives rise to an energy exchange between the two conjugate CFMs.The time histories, frequency spectra, and phase portraits obtained demonstrate a chaotic vibration of the SWCNT in nonlinear field.The total energy distribution over the nanotube surface evolves in a complex pattern, maintaining symmetry along the longitudinal direction and periodicity along the circumferential direction, where the initial energy on the edge sections of the CNT is no longer preserved in amplitude in the time.The energy exchange between the two conjugate CFMs generates a pure travelling wave moving circumferentially around the SWCNT.



---

*Source: 3241698-2019-04-24.xml*