Abstract
We establish an analytical model to investigate the surface effects on the vibration and instability of fluidconveying nanotubes embedded in viscoelastic medium. Based on nonlocal elastic theory and EulerBernoulli beam theory, the vibration equation of fluidconveying nanotubes is established with considering three typical boundary conditions. The effects of both inner and outer surface layers on the nanotubes are taken into consideration and the KelvinVoigt model is introduced. The results show that the boundary conditions of system, the damping and elastic coefficient of the surrounding medium, thickness of nanotubes and aspect ratios have significant effects on the dynamic behaviors of the nanotubes. The damping parameter of the viscoelastic foundation causes an obvious reduction of the critical flow velocity. For smaller tube thickness, larger aspect ratio or higher elastic parameter of surrounding foundation, the stability of the nanotubes may be greatly enhanced. This article might be helpful for the design and improvement of nanotubes for fluidconveying applications embedded in elastic medium in nanoelectromechanical systems and microelectromechanical systems.
1. Introduction
For the wide applications of nanotubes/nanopipes in nanoelectromechanical systems (NEMS) and microelectromechanical systems (MEMS), many attempts have been made to study their properties [1, 2]. Due to the hollow cylindrical geometry structure of nanotubes, they are expected to be used as nanofluidic devices and systems, with potential application in fluid filtration devices, biomimetic selective transport of ions, targeted drug delivery devices and so on [35].
As a small scale fluidstructure interaction system, the fluidconveying nanotubes show higher sensitivity to the vibration characteristics. The research of flowinduced vibration and instability is of fundamental significance. In the past two decades, the literature on this topic has been constantly expanding. Many researchers have done much work in this area [612]. The systems they studied range from straight to curved configurations of nanotubes with either single wall or multi walled, linear and nonlinear vibration [1316].
In classical continuum mechanics, the surface energy is usually small and negligible. But when materials and structures shrink to nanometers, it plays a critical role in their static or dynamical behavior due to the high specific surface area of nanomaterials [17]. He and Lilley [18] studied the effects of surface tension on the static bending of nanowires. Farshi et al. [19] researched the vibration of nanotubes with consideration of surface effects utilizing Timoshenko beam model. Gheshlaghi and Hasheminejad [20] investigated the surface effects on the nonlinear vibration of nanobeams based on the EulerBernoulli beam theory including the von Kármán geometric nonlinearity.
The vibration and buckling instability of fluidconveying nanotubes by accounting for surface layers were analyzed by Wang [21], who found that the natural frequencies and critical flow velocity are strongly related to the surface effect. With the increasing ratio of surface layer to volume (nanotube wall) at the nanoscale, the vibration characteristics of nanotubes conveying fluid were shown to exhibit distinct size dependences. Furthermore, the surface effect on the buckling and postbuckling behavior of a nanobeam with supported ends containing flowing fluid is investigated using a nonliner model [22].
In this paper, an analytical model is established to invesitigate the surface effect on the vibration and instability of fluidconveying nanotubes embedded in visoelastic medium by nonlocal Eulerbernoulli beam model. It will be shown that the current model will be helpful for designing the nanotubes for fluidconveying applications embedded in viscoelastic medium in NEMS and MEMS.
2. Fundamental equations
The circular nanotube embedded in viscoelastic medium is modeled as a straight and slender nanobeam with length $L$, as shown in Fig. 1. From reference [12] we get the vibration equation of fluidconveying nanotube embedded in viscoelastic medium based on nonlocal EulerBernoulli beam theory as follow:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}({e}_{0}a{)}^{2}\frac{{\partial}^{2}}{\partial {x}^{2}}\left[m\frac{{\partial}^{2}w}{\partial {t}^{2}}+{m}_{f}\left({u}^{2}\frac{{\partial}^{2}w}{\partial {x}^{2}}+2u\frac{{\partial}^{2}w}{\partial x\partial t}+\frac{{\partial}^{2}w}{\partial {t}^{2}}\right)+kw+c\frac{\partial w}{\partial t}\right]=0,$
where $w(x,t)$ is the transverse displacements of the nanotube along the $x$axis at time $t$, $EI$ is the flexural rigidity of the tube without surface layers, $m$ and ${m}_{f}$ are the mass of nanotube and the mass of the fluid per unit length respectively, $u$ is the flow velocity of the internal fluid. $k$ and $c$ are the elastic stiffness and damping coefficient of the surrounding viscoelastic medium based on the KelvinVoigt model, respectively. ${e}_{0}a$ represents the effect of nonlocal elasticity, and when it equals to zero, the equation reduces to classical condition.
Fig. 1Schematic diagram of a) a fluidconveying nanotube embedded in viscoelastic medium; b) the circular cross section with surface layers
a)
b)
Surface effects on the mechanical behaviour of nanomaterials can be examined by considering surface energy and surface residual stresses. The surface energy mainly acts as an increased effect on bending rigidity, the effect with surface residual stress acts as distributed transverse loads and can be calculated by LaplaceYoung equations [19]. The additional flexural rigidity of nanotube due to the inner and outer surface layers may be written as [19]:
where ${E}_{s}$ and ${t}_{0}$ are Young’s modulus and thickness of both surface layers, ${d}_{i}$ and ${d}_{o}$ are the inner and outer diameters of the nanotube, respectively, as shown in Fig. 1(b). The distributed transverse loads induced by the surface residual stress is $q\left(x\right)={\mathrm{\Pi}}_{0}\left({\partial}^{2}w/\partial {x}^{2}\right)$ and the constant ${\mathrm{\Pi}}_{0}$ is given by:
where ${\tau}_{0}$ is the residual surface tension.
Considering the influence of two surface layers, the equation of motion of the fluidconveying nanotube embedded in viscoelastic medium can be written as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}({e}_{0}a{)}^{2}\frac{{\partial}^{2}}{\partial {x}^{2}}\left[\left({m}_{f}{u}^{2}{\mathrm{\Pi}}_{0}\right)\frac{{\partial}^{2}w}{\partial {x}^{2}}+2{m}_{f}u\frac{{\partial}^{2}w}{\partial x\partial t}+\left(m+{m}_{f}\right)\frac{{\partial}^{2}w}{\partial {t}^{2}}+kw+c\frac{\partial w}{\partial t}\right]=0.$
Three typically boundary conditions are considered:
(i) A beam simply supported at both ends, i.e. pinnedpinned condition (PP):
(ii) A beam clamped at both ends, i.e. clampedclamped condition (CC):
(iii) A beam clamped at one end and simply supported at the other end, i.e. clampedpinned condition (CP):
For simplicity, introducing dimensionless variables as follows:
$K=\frac{k{L}^{4}}{EI},\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}C=\frac{c{L}^{2}}{\sqrt{EI({m}_{c}+{m}_{f})}},\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{e}_{n}=\frac{{e}_{0}a}{L},\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\alpha =\frac{\mathrm{\hslash}}{EI},\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\gamma =\frac{{\mathrm{\Pi}}_{0}{L}^{2}}{{\pi}^{2}EI}.$
Eq. (4) can be written in the dimensionless form:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{e}_{n}^{2}\frac{{\partial}^{4}\eta}{\partial {\xi}^{2}\partial {\tau}^{2}}+2\sqrt{\beta}U\frac{{\partial}^{2}\eta}{\partial \xi \partial \tau}C{e}_{n}^{2}\frac{{\partial}^{3}\eta}{\partial {\xi}^{2}\partial \tau}+C\frac{\partial \eta}{\partial \tau}+K\eta +\frac{{\partial}^{2}\eta}{\partial {\tau}^{2}}=0.$
The boundary conditions are converted to the following:
PP condition:
CC condition:
CP condition:
3. Numerical results
We use the Differential Quadrature Method (DQM) to obtain the solutions of Eqs. (9)(12). For calculations convenience, the nanotube with circular crosssection area is assumed to be made of anodic alumina with crystallographic of [1 1 1] direction. The parameters may be chosen as follows:$E=$ 70 GPa, ${\rho}_{t}=$ 2700 kg/m^{3} (mass density of the tube), ${E}_{s}{t}_{0}=$ 5.1882 N/m, ${\tau}_{0}=$ 0.9108 N/m [19], ${\rho}_{w}=$ 1000 kg/m^{3} (mass density of the internal water), ${e}_{0}a=$ 1 nm, $c=\mathrm{}$1.02×10^{4} Pa·S, ${d}_{i}=$ 20 nm. According to the work by Farshi et al. [19], for nanostructures made of anodic alumina material, surface effect maybe the most important factor associated with the sizedependent properties.
Fig. 2 shows the variation of the first natural frequency with respect to the dimensionless flow velocity, for different boundary conditions and thicknesses ($t=({d}_{o}{d}_{i})/2$) of the nanotube. It can be seen that the natural frequency decreases as the flow velocity increases and the nanotube buckles when the flow velocity increase to a critical value. The results demonstrate that as the bending stiffness of the nanotube rises from PP to CC, the corresponding natural frequencies and the dimensionless critical flow velocities increase. It can be observed that surface effects generally increase the natural frequency and critical flow velocity comparing to the counterpart without surface effects especially for smaller tube thickness. The thickness has a negative influence on the natural frequency, i.e., the smaller the thickness, the larger the natural frequency.
Fig. 2Natural frequency as functions of dimensionless flow velocity
Fig. 3Critical flow velocity as functions of tube thickness
Fig. 3 shows the surface effects on the critical flow velocity, for various values of tube thickness and different aspect ratios. It is found that the influence of surface layers on the critical flow velocity is significant, especially for nanotubes with smaller tube thickness or larger aspect ratio. For example, in the case of $L/{d}_{o}=$ 20 and $t=$ 0.5 nm, the critical flow velocity predicted with surface layers is about 3.51 times than that predicted without surface layers, which implies that the surface layers enhance the stability of fluidconveying nanotubes. The influence of surface layers on critical flow velocity decreases as the tube thickness increases.
Fig. 4Critical flow velocity as functions of dimensionless damping coefficient
Fig. 5Critical flow velocity as functions of dimensionless Winkler parameter
The effect of the surrounding viscoelastic medium on the critical flow velocity of nanotube is shown in Fig. 4 and Fig. 5 for different aspect ratios. The dimensionless parameters $C$ and $K$ are used to demonstrate the effects of the viscoelastic behavior of the foundation. It can be concluded from Figs. 4 and 5 that the critical velocity is significantly influenced by the damping and elastic parameter especially for smaller aspect ratio. By increasing the dimensionless damping coefficient $C$ or decreasing the dimensionless elastic parameter $K$, the dimensionless critical flow velocity decreases. This implies that a higher damping property dissipates the vibrational energy quickly and the resonant frequency vanishes at a lower flow velocity, and the elastic parameter enhance the stability of the fluidconveying nanotube.
4. Conclusions
In summary, an analytical model is developed to investigate the surface effects on the vibration and stability analysis of fluidconveying nanotubes embedded in viscoelastic medium. The results show that the damping and elastic coefficient of medium, thickness of nanotubes and aspect ratios have significant effects on the dynamic behaviors of the nanotubes. This article might be helpful for the design and improvement of nanotubes for fluidconveying applications embedded in elastic medium in NEMS and MEMS.
References

Gibson R. F., Ayorinde E. O., Wen Y. F. Vibrations of carbon nanotubes and their composites: a review. Composites Science and Technology, Vol. 67, 2007, p. 128.

Li C., Thostenson E. T., Chou, T.W. Sensors and actuators based on carbon nanotubes and their composites: a review. Composites Science and Technology, Vol. 68, Issue 6, 2008, p. 12271249.

Gadd G. E., Blackford M., Moricca S., Webb N., Evans P. J., et al. The world’s smallest gas cylinder. Science, Vol. 277, 1997, p. 933936.

Yang Z., Zhang Y., Yang Y., Sun L., Han D., et al. Pharmacological and toxicological target organelles and safe use of singlewalled carbon nanotubes as drug carriers in treating Alzheimer disease. NanomedicineNanotechnology Biology and Medicne, Vol. 6, 2010, p. 427441.

Gao Y. H., Bando Y. Nanotechnology: Carbon nanothermometer containing gallium. Nature, Vol. 415, 2002, p. 599599.

Yoon J., Ru C. Q., Mioduchowski A. Vibration and instability of carbon nanotubes conveying fluid. Composites Science and Technology, Vol. 65, 2005, p. 13261336.

Yoon J., Ru C. Q., Mioduchowski, A. Flowinduced flutter instability of cantilever carbon nanotubes. International Journal of Solids and Structures, Vol. 43, 2006, p. 33373349.

Chang W. J., Lee H. L. Free vibration of a singlewalled carbon nanotube containing a fluid flow using the Timoshenko beam model. Physical Letters A, Vol. 373, 2009, p. 982985.

Khosravian N., RafiiTabar H. Computational modeling of the flow of viscous fluids in carbon nanotubes. Journal of Physics D: Applied physics, Vol. 40, 2007, p. 70467052.

Yan Y., He X. Q., Zhang L. X., Wang Q. Flowinduced instability of doublewalled carbon nanotubes based on an elastic shell model. Journal of Applied Physics, Vol. 102, 2007, p. 044307.

Wang L. Dynamical behaviors of doublewalled carbon nanotubes conveying fluid accounting for the role of small length scale. Computational Materials Science, Vol. 45, 2009, p. 584588.

Zhen Y. X., Fang, B. Thermalmechanical and nonlocal elastic vibration of singlewalled carbon nanotubes conveying fluid. Computational Materials Science, Vol. 49, 2010, p. 276282.

Xia W., Wang L. Vibration characteristics of fluidconveying carbon nanotubes with curved longitudinal shape. Computational Materials Science, Vol. 49, 2010, p. 99103.

Ghavanloo E., Rafiei M., Daneshmand F. Inplane vibration analysis of curved carbon nanotubes conveying fluid embedded in viscoelastic medium. Physics Letters A, Vol. 375, 2011, p. 19941999.

Wang Y. Z., Li F. M. Nonlinear free vibration of nanotube with small scale effects embedded in viscous matrix. Mechanics Research Communications, Vol. 60, 2014, p. 4551.

Zhen Y. X., Fang B. Nonlinear vibration of fluidconveying singlewalled carbon nanotubes under harmonic excitation. International Journal of NonLinear Mechanics, Vol. 76, 2015, p. 4855.

She H., Wang B. A geometrically nonlinear finite element model of nanomaterials with consideration of surface effects. Finite Elements in Analysis and Design, Vol. 45, 2009, p. 463467.

He J., Lilley C. M. Surface effect on the elastic behavior of static bending nanowires. Nano Letters, Vol. 8, Issue 7, 2008, p. 17981802.

Farshi B., Assadi A., Aliniaziazi A. Frequency analysis of nanotubes with consideration of surface effects. Applied Physics Letters, Vol. 96, 2010, p. 093105.

Gheshlaghi B., Hasheminejad S. M. Surface effects on nonlinear free vibration of nanobeams. Composites: Part B, Vol. 42, 2011, p. 934937.

Wang L. Vibration analysis of fluidconveying nanotubes with consideration of surface effects. Physica E, Vol. 43, 2010, p. 437439.

Wang L. Surface effect on buckling configuration of nanobeams containing internal flowing fluid: A nonlinear analysis. Physica E, Vol. 44, 2012, p. 808812.
About this article
This work is supported by National Natural Science Foundation of China under Grant No. 11402087 and The Fundamental Research Funds for the Central Universities (No. 2016MS64).