Forced Transverse Vibration of a Closed Double Single-Walled Carbon Nanotube System Containing a Fluid with Effect of Compressive Axial Load

Based on the Rayleigh beam theory, the forced transverse vibrations of a closed double single-walled carbon nanotube (SWCNT) system containing a fluid with a Pasternak layer in-between are investigated. It is assumed that the two single-walled carbon nanotubes of the system are continuously joined by a Pasternak layer and both sides of SWCNTs containing a fluid are closed. The dynamic responses of the system caused by arbitrarily distributed continuous loads are obtained. The effect of compressive axial load on the forced vibrations of the double single-walled carbon nanotube system is discussed for one case of particular excitation loading. The properties of the forced transverse vibrations of the system are found to be significantly dependent on the compressive axial load. The steady-state vibration amplitudes of the SWCNT decrease with increasing of length of SWCNT. Vibrations caused by the harmonic exciting forces are discussed, and conditions of resonance and dynamic vibration absorption are formulated.The SWCNT-type dynamic absorber is a new concept of a dynamic vibration absorber (DVA), which can be applied to suppress excessive vibrations of corresponding SWCNT systems.


Introduction
Carbon nanotubes (CNTs) are becoming increasingly promising materials for modern technology applications (nanocomposites, nanodevices, nanoelectronics, and nanomedicine) due to their exceptional thermomechanical, chemical, and optoelectronic properties.The superior strength and unusual deformability of CNTs have been a subject of numerous recent studies.For example, the classical Euler and Timoshenko elastic-beam models have been used effectively to study the deformation of CNTs, including static deflections, column buckling, resonant frequencies, and modes, as well as wave propagation such as Wong et al. [1], Poncharal et al. [2], Yoon et al. [3], Li and Chou [4], Ru [5], Wang and Hu [6], Zhang et al. [7], and Wang et al. [8].
Recently, carbon nanotube (CNT) filled with fluids has been of a great interest.This is because CNTs promise many new applications in nanobiological devices and nanomechanical systems such as fluid conveyance and drug delivery.
Chang and Lee [9] studied the free vibration of a singlewalled carbon nanotube containing a fluid flow by using the Timoshenko beam model.They analyzed the effects of the flow velocity and the aspect ratio of length to diameter on the vibration frequency and mode shape of the SWCNT.Their results showed the effects of rotary inertia and transverse shear deformation result in a reduction of the vibration frequencies, especially for higher modes of vibration and short nanotubes.
The systems connected by an elastic layer constitute one group of such mechanical structures which are commonly encountered in mechanical, construction, and aeronautical industry.Such systems are important as they give the initial approximation of the solution and dynamic behavior of the system at slight motion.The problem of two elastically connected beams joined by the Winkler elastic layer emerged in order to determine the conditions for the behavior of the system acting as a dynamic absorber in technical practice.
A mathematical model was developed by Seelig and Hoppmann [10].They investigated the problem of an impulse load effect on a beam and produced a system of partial differential equations describing its vibration.The obtained theoretical and experimental results confirmed a sound approximation of an analytical solution obtained for slender beams at small transverse motions using the Euler-Bernoulli theory.
Oniszczuk [11,12] analyzed the problem of free and forced vibration of two elastically connected Euler-Bernoulli beams.He determined analytical solutions for eigenfrequencies, amplitudinous functions, and vibration modes.He discussed the effect of stiffness which the elastic interlayer had on the frequencies and amplitudes of the system.Also, the conditions for the occurrence of resonance and the behavior of the system as a dynamic absorber are determined.
The analysis of the system composed of two connected beams was carried on by Zhang et al. [13,14].In their work, free and forced vibrations by two elastically connected Euler-Bernoulli beams affected by axial compression forces are investigated.They presented analytical solutions for natural frequencies of the system in the function of axial compression force impact and their effect on the vibration amplitude.They determined the codependency between the system's critical force and the Euler critical load in the function of an axial force of the other beam.
Stojanović et al. [15] analyzed free vibration and static stability of two elastically connected beams with Winkler elastic layer in-between with the influence of rotary inertia and transverse shear.They described the motion of the system by a homogeneous set of two partial differential equations.Also, they determined the natural frequencies and associated amplitude ratios of an elastically connected double-beam complex system and the analytical solution of the critical buckling load.
Stojanović and Kozić [16] discussed the case of forced vibration of two elastically connected beams with Winkler elastic layer in-between and the effect of axial compression force on the amplitude ratio of system vibration for three types of external forcing (arbitrarily continuous harmonic excitation, uniformly continuous harmonic excitation, and concentrated harmonic excitation).They determined general conditions of resonance and dynamic vibration absorption.In the other paper, Stojanović et al. [17] discussed the analytic analysis of static stability of a system consisting of three elastically connected Timoshenko beams on an elastic foundation.They provided expressions for critical force of the system under the influence of elastic Winkler layers.
In an attempt to find a physically close and mathematically simple representation of an elastic foundation for these materials, Pasternak proposed a foundation model consisting of a Winkler foundation with shear interactions.This may be accomplished by connecting the ends of the vertical springs to a beam consisting of incompressible vertical elements, which deforms only by transverse shear.

The Beam Model for SWCNTs
The Rayleigh beam theory, which has taken rotary inertia into account, is applied to analyze forced transverse vibration of  a closed double SWCNT system containing a fluid with effect of compressive axial load with Pasternak interlayer.Figure 1 shows a double SWCNT system containing a fluid with a Pasternak layer in-between with length of , which was subjected to axial compressions  1 and  2 that are positive in compression, which do not change with time and arbitrarily distributed transverse continuous loads  1 and  2 that are positive when they act downward.
An element of deflected differential Rayleigh beam of length  with Pasternak layer between two cross sections is shown in Figure 2.
Let the functions of longitudinal and transverse motion of the Rayleigh beam system be as follows:  1 (, , ),  1 (, , ),  2 (, , ), and  2 (, , ).We have Deformation in the  direction in the function of motion and the stress-deformation relation according to Hooke's law are as follows: Virtual work of inertial forces is expressed as Virtual work of internal forces is Virtual work of external forces is given by Based on the principle of virtual work   +   +   = 0,  = 1, 2, and (3) to ( 5), we have By successive application of Green's theorem on expressions (6), we obtain the set of coupled differential equations of forced transverse vibration of a double SWCNT system containing a fluid with effect of compressive axial load with Pasternak layer in-between: where If we introduce  0 = 0 in the equations, the effect of shear foundation will be omitted and we can obtain the vibration equations of a double SWCNT system containing a fluid with Winkler elastic layer in-between.Equations ( 7) can be reduced to fourth-order partial differential equations for forced vibration of a double SWCNT system model: where The initial conditions in general form and boundary conditions for simply supported SWCNT of the same length  are assumed as follows: (0, ) =    (0, ) =   (, ) =    (, ) = 0,  = 1, 2. (13)

Solution of Equations
Equations ( 9) and ( 10) represent forced vibrations of a double SWCNT system.The natural frequencies and the corresponding mode shapes of the system should be obtained by solving the undamped free vibration with appropriate boundary conditions.Assuming time harmonic motion and using separation of variables, the solutions to ( 9) and ( 10) with the governing boundary conditions ( 13) can be written in the form where Substitution of ( 14) into ( 9) and ( 10) yields ordinary differential equations for the double SWCNT system where

Forced Vibration of a Double SWCNT System
First, we should solve the undamped free vibration of the double SWCNT system The solutions of (18) can be assumed to have the following forms: Substituting ( 19) into (18) results in the following system of homogeneous algebraic equations for the unknown constants: When the determinant of the coefficients in (20) vanishes, nontrivial solutions for the constants   and   can be obtained, which yield the following frequency (characteristic) equation: The roots of (21) are where For each of the natural frequencies, the associated amplitude ratio of vibration modes of the two beams is given by Applying the classical modal expansion method, the forced vibrations of a double SWCNT system subjected to arbitrarily distributed continuous loads can be determined.Following the above analysis for the undamped free transverse vibration, particular solutions of nonhomogeneous differential equations ( 9) and (10) representing forced vibrations of a double SWCNT system model can be assumed in the following: Introduction of ( 25) and ( 26) into ( 9) and ( 10) results in By multiplying relations ( 27) and (28) by the eigenfunction   and then integrating them with respect to  from 0 to  and using orthogonality condition where now we have From ( 24) and (31) after some algebra we obtain From (32), we obtain where Also,  1 =  1 / 1  and  2 =  2 / 2 .By combining (25), (26), and (33), the forced vibrations of a double SWCNT system with Pasternak layer in-between which contains a fluid can be described by Now these general solutions (35) and (36) are used to find the vibrations of a double SWCNT system containing a fluid.
In the following, we conduct an analysis of forced vibrations for the case of uniformly distributed continuous harmonic load.For simplicity of further analysis, it is assumed that only one of the two carbon nanotubes is subjected to the exciting load.Without loss of generality, we suppose (see Figure 3) Substituting (37) into (34), we obtain where Substituting ( 38) and (39) into ( 35) and (36) gives where ] . (42) Ignoring the free response, the forced vibrations of the double SWCNT system can be obtained by The following fundamental conditions of resonance and dynamic vibration absorption have practical significance: (1) Resonance (2) Dynamic vibration absorption

Numerical Results and Discussion
For simplicity, it is assumed that both SWCNTs are geometrically and physically identical and the elastic interlayer is Winkler.We assumed that the fluid in the SWCNT is water.The values of the parameters characterizing properties of the system are shown in Table 1: If the axial compressions vanish, for the double SWCNT system we have where Also, To determine the effect of compressive axial load on the steady-state vibration amplitudes  1 and  2 of the system, the results under compressive axial load and those without axial load are compared.We introduce the relation We use nondimensional ratio where  =  2 / 2 is the smallest load at which the SWCNT ceases to be in stable equilibrium under axial compression.With the vibration mode number  = 3 and the exciting frequency Ω = 0.6  , the effects of compressive axial load on the steady-state vibration amplitudes  1 and  2 of the SWCNT are represented by the ratios  1 and  2 , shown in Figures 4 and 5, respectively.As can be seen, the ratio  1 decreases with the increase of the axial compression, which implies that the magnitude of the steady-state vibration amplitude  1 becomes smaller when the axial compression increases, and the ratio  2 increases with the increase of the axial compression, which implies that the magnitude of the steady-state vibration amplitude  2 becomes larger when the axial compression increases.Numerical values of the ratios  1 and  2 for different axial compression ratio  of the two SWCNTs are shown in Table 2.It can be observed that the effect of compressive axial load on the magnitude of  1 is almost independent of the axial compression ratio  of the two SWCNTs whereas it is significantly dependent on the magnitude of  2 .On the other hand, with the axial compression ratio  = 0.5 and the exciting frequency Ω = 0.6  , steady-state vibration amplitudes  1 and  2 are shown in Figures 6  and 7 for different mode shape number , respectively.It can be seen that, with the same axial compression, the ratios  1 and  2 diminish with the increasing vibration mode number , which implies that the magnitudes of the steady-state vibration amplitudes  1 and  2 get smaller when the vibration mode number  becomes larger.
Table 3 shows the effects of compressive axial load on the difference between the steady-state vibration amplitude ratios  1 and  2 of the SWCNT for different mode number .It can be observed that the differences between ratios  1 and  2 of the SWCNT increase with increasing the dimensionless parameter .
Table 4 shows the effects of length of SWCNT on the difference between the steady-state vibration amplitude ratios  1 and  2 of the SWCNT for mode number  = 5 and axial compression ratio  = 0.5.It can be observed that the differences between ratios  1 and  2 of the SWCNT decrease with increasing of length of SWCNT.
We accounted the results for three cases of shear foundation modulus of Pasternak layer  0 = 0, 100, 1000.With the vibration mode number  = 3 and the exciting frequency   5. Numerical results of the ratios  1 and  2 show that the difference with ratios  1 and  2 decreases with the increase of shear modulus of the Pasternak layer for different axial compression ratio , but it can be observed that the effect of shear foundation modulus of Pasternak layer on the magnitude of  1 is almost independent of the axial compression ratio  of the two beams whereas it is significantly dependent on the magnitude of  2 .

Conclusions
Based on the Rayleigh beam theory, the forced transverse vibrations of a closed double single-walled carbon nanotube (SWCNT) system containing a fluid with a Pasternak layer in-between, under compressive axial loading for one case of particular excitation loading, are studied.The dynamic responses of the system caused by arbitrarily distributed continuous loads are obtained.The effect of compressive axial load and length of SWCNT on the forced vibrations of the double single-walled carbon nanotube system is discussed for one case of particular excitation loading.The magnitudes of the steady-state vibration amplitudes of the beam are dependent on the axial compression and length of SWCNT.The effect of compressive axial load on the magnitude of  1 is almost independent of the axial compression ratio  of the two SWCNTs whereas it is significantly dependent on the magnitude of  2 .Also, it is shown with the same axial compression that the ratios  1 and  2 diminish with the increasing vibration mode number , which implies that the magnitudes of the steady-state vibration amplitudes  1 and  2 get smaller when the vibration mode number  becomes larger.The ratio  1 decreases with the increase of the axial compression, which implies that the magnitude of the steady-state vibration amplitude  1 becomes smaller when the axial compression increases, and the ratio  2 increases with the increase of the axial compression, which implies that the magnitude of the steady-state vibration amplitude  2 becomes larger when the axial compression increases.The differences between ratios  1 and  2 of the SWCNT decrease with increasing of length of SWCNT.The SWCNT-type dynamic absorber is a new concept of a dynamic vibration absorber (DVA), which can be applied to suppress excessive vibrations of corresponding SWCNT systems.At a given vibration amplitude, increase of the nanotube length leads to decrease of the frequencies.Also we concluded that the ratios  1 and  2 decrease with increasing of the shear foundation modulus of Pasternak layer ( 0 ).

Figure 2 :
Figure 2: Deflected differential Rayleigh beam element with Pasternak layer in-between.

Figure 3 :
Figure 3: Double single-walled carbon nanotube system subjected to harmonic uniform distributed continuous load.

2 Figure 7 :
Figure 7: Relationship between ratio  2 and dimensionless parameter  for different mode number .
= (, ):Transverse displacements of the carbon nanotubes /: G l o b a lr o t a t i o no ft h ec r o s ss e c t i o n   (): Unknown time function   (): Known mode shape function for simply supported carbon nanotube   () ( = , ): Unknown time function corresponding to the natural frequencies     : Kronecker delta function   : E i g e n f u n c t i o n (): A m p l i t u d eo ft h el o a d Ω: E x c i t i n g f r e q u e n c y o f t h e l o a d  0 1 ,  0 2 : Steady-state vibration amplitudes of the two carbon nanotubes without axial compression   : Natural frequencies of the system : Cross section area of the hollow cylindrical nanotube : Y o u n g ' s m o d u l u s : Second moment of inertia of the nanotube   : Mass density of the carbon nanotube   : Mass density of the fluid : W i n k l e r f o u n d a t i o n m o d u l u s  0 : S h e a r f o u n d a t i o n m o d u l u s : Inner diameter   : Th i c k n e s s .

Table 1 :
Values of the parameters characterizing properties of the system.

Table 2 :
Steady-state vibration amplitudes ratios  1 and  2 of SWCNT for different compressive axial ratio .

Table 3 :
Steady-state vibration amplitudes ratios  1 and  2 for different mode shape number .

Table 4 :
Steady-state vibration amplitudes ratios  1 and  2 for different length of SWCNT, ., the effects of Pasternak layer on the steady-state vibration amplitudes  1 and  2 of the Rayleigh beam are represented by the ratios  1 and  2 , shown in Table5.Numerical values of the ratios  1 and  2 for different shear foundation modulus of Pasternak layer with different axial compression ratio  of the two beams are shown in Table

Table 5 :
Effects of shear foundation modulus of Pasternak layer on the steady-state vibration amplitudes ratios  1 and  2 for different compressive axial ratio .