Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 517695, 9 pages

http://dx.doi.org/10.1155/2013/517695

Research Article

Nonlinear Analysis of Electrically Actuated Carbon Nanotube Resonator Using a Novel Discretization Technique

Hassen M. Ouakad,¹ Fehmi Najar,² and Oussama Hattab²

¹Mechanical Engineering Department, King Fahd University of Petroleum and Minerals, P.O. Box 31261, Dhahran, Saudi Arabia

²Applied Mechanics and Systems Research Laboratory, Tunisia Polytechnic School, University of Carthage, BP 743, 2078 La Marsa, Tunisia

Received 21 April 2013; Revised 22 August 2013; Accepted 22 August 2013

Academic Editor: Sarp Adali

Copyright © 2013 Hassen M. Ouakad et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Static as well as dynamic analyses have been performed on clamped-clamped carbon nanotube (CNT) resonator. The nonlinear CNT model is investigated with a novel discretization technique: a differential quadrature method (DQM) to discretize the spatial variables and a finite difference method (FDM) for limit-cycle solutions. Parametric study is performed by varying the electric load, as well as the initial curvature (due to fabrication). It is found that the pull-in voltage decreases nonlinearly with initial curvature and linearly with residual stresses. The eigenvalue problem is also solved to obtain the bending natural frequencies of the system as function of the DC voltage as well as the initial curvature of the CNT. Frequency-response curves near some selected resonant frequencies are plotted to better understand the nanotubes' dynamic behavior. Different linear and nonlinear phenomena are depicted such as dynamic pull-in, hardening, and softening behavior and veering of the odd modes. We have found that even when exciting the CNT near its first natural frequency, the vibration mode located at the veering process significantly alters the CNT's motion and hence may decrease its overall quality factor.

1. Introduction

Since their discovery [1], carbon nanotubes (CNTs) are in the center of interest of researchers thanks to their revolutionary characteristics such as high Young’s modulus and low density leading to light structures with high mechanical performance. In addition to that, promising applications in several fields like optics, electronics, and resonant sensors are investigated. Their natural frequencies are found to be very high compared to other structures, which permit to design high frequency resonators [2].

Static and dynamic analyses of CNTs rise difficult issues in solving the highly nonlinear governing equation. In this case the nonlinearity is both of the geometric type (mid-plane stretching) and of the loading type (electric load) [3]. Several methods were attempted to solve the problem using the Galerkin procedure along with shooting techniques [4] or molecular dynamics simulations [5, 6].

Previous studies have focused on mechanical as well as electric response of CNTs exploring their interesting characteristics such as Young’s modulus and pull-in properties. Experimental investigations of the vibration of electrostatically actuated CNTs have been conducted to characterize their effective Young modulus (Poncharal et al. [7]). On the other hand, others have found that Van der Waals effect on static response of CNTs is negligible for gap widths greater than 3 nm [5, 6]. Using the reduced order model (ROM) with the Galerkin method, Ouakad and Younis [4] have investigated the effect of geometric parameters and initial curvature on pull-in and the dynamic response of CNTs. They are found to hold a hardening behavior in almost all of the cases until applied voltage is high enough (near pull-in) to change the dynamic response into softening type. On the other hand, cantilever beams are found to hold softening behavior [4]. Natural frequencies variation with initial curvature has revealed the existence of phenomenon called modes veering. Ouakad and Younis [4] investigated the effect of van der Waals forces on the fundamental natural frequency of CNTs as well and found them significant for CNTs of very small gap widths.

Using the nonlocal Timoshenko model, the effect of temperature change on mode shapes is investigated [8]. Studying the dynamic response of an array of nanotubes has revealed that interaction between tubes effects considerably the mode shapes and their stability [9]. Nonlinear bifurcations have been shown to result in a transition from quasi-periodic to chaotic behavior in dynamic response [10]. Double-walled carbon nanotubes (DWCNTs) have been analyzed by Fang et al. [11] determining the effect of different nonlinearities on their dynamic response and by Hajnayeb and Khaddem [12] who found that both layers vibrate with the same frequency. Huang [13] has investigated the effects of adhesion temperature, tube length, and peeling velocity on the dynamic behavior of a DWCNT on an Au substrate.

Many researchers have fabricated and tested systems based on CNTs such as torsional electromechanical systems [14], strain sensors [15], biosensors [16, 17], and pH sensing system [18]. For example, Cho et al. [19] have designed a cantilever nanotube which has been tested successfully using multiscale analysis with comparison to experimental results. The cantilever nanotube is found to be sensitive to geometric parameters changes altering the dynamic response from hardening behavior to softening behavior. Graham et al. [20] have applied the 2D Fourier transform spectroscopy to semiconducting SWCNTs determining spectral and time-domain map of exciton-phonon assisted excitations. Several studies have focused on the energetic behavior of CNTs. As examples, Peng et al. [21] have studied experimentally using TEM the dissipation of energy due to air damping and contact loss. Greaney and Grossman [22] have used molecular dynamics to simulate energy transfer between nanotubes and have found that efficient energy transfer is allowed by sharp resonance effects.

The objective of the paper is to solve both static and dynamic nonlinear problems using novel discretization technique: differential quadrature method and finite difference method. These methods have been applied successfully to MEMS electrostatic actuators in order to investigate their nonlinear dynamics [23–25]. The considered model accounts for geometrical through the midplane stretching and also takes into account the fully nonlinear form of the electrostatic force. In this paper, we present the proposed ROM by applying the DQM. Then, free-vibration problem under the actuation of DC load is solved in order to predict natural frequencies and mode shapes of the system. Finally, we attempt to give a response to dynamic transformations occurring during veering phenomenon by plotting the frequency-response curves under AC and DC loading around different excitation frequencies.

2. Problem Formulation and Reduced-Order Model

2.1. Problem Formulation

A hollow cylindrical Euler-Bernoulli beam of length  nm, shell thickness  nm, radius nm,  nm is considered in this paper to model an electrostatic CNT resonator initially curved in the direction of the lower electrode [26] (Figure 1). It is assumed to have Young’s modulus TPa [7] and a mass density [27]. The mathematical formulation of the static and dynamic behaviors of the clamped-clamped CNT in its nondimensional form is given by [4]

where the longitudinal coordinates and time are normalized with respect to the length and the time constant , respectively. is the nondimensional transverse deflection toward the substrate, the initial shape function is [28, 29] and

where is the cross-section area, is the moment of inertia, is the natural frequency, is the axial force modeling the residual stress, and is the quality factor.

Figure 1: Schematic of the CNT resonator.

2.2. Reduced-Order Model Using the Differential Quadrature Method

The differential quadrature method (DQM) is used to solve the space dependent partial differential equation by transforming it into ordinary differential equations describing the motion of the CNT with respect to time at preselected grid points . The derivatives of the deflection with respect to space variable are expressed as a weighted linear sum of the deflection at all grid points. Following [23], we discretize the integral terms using the Newton-Cotes formula at the same grid points. Using the boundary conditions given in (2), we end up with the following ODEs,

and the coefficients and are given in [23].

3. Static Analysis

3.1. Convergence and Validation

In this section, we solve the algebraic system obtained by dropping the time dependent terms in (4). First, we perform a convergence study in which the static response of the CNT is calculated by varying the number of grid points varied from 5 to 15. The results are compared to those obtained by Ouakad and Younis [29] and calculated using the Galerkin method. In Figure 2, we show that the number of grid points required to have converged results is 11. Hence, in all of the following results, we will be using at least 11 points in the DQM to assure convergence.

Figure 2: Variation of the normalized static deflection with DC voltage for different numbers of grid points.

3.2. Effect of Residual Stress

Residual stresses have thermal origins due to fabrication techniques of the CNT [30]. We model this effect here by adding an axial force constant along the length of the CNT. We study here the influence of these residual stresses on the static response of the CNT for both compression and tension cases. It is always assumed to be lower than the first buckling load .

In Figure 3 we show that residual stresses have a relatively small effect on the static response of the CNT. Indeed an increase of the residual stress induces an increase of the pull-in voltage for different values of the initial curvature.

Figure 3: Variation of the pull-in voltage with the residual stresses at different initial curvatures.

3.3. Effect of Initial Curvature

We investigate next the effect of the initial curvature into the static response of the CNT. The initial curvature is supposed to vary from 100 nm to 200 nm.

The results in Figure 4 are in a good agreement with results reported in [29]. In Figure 5, we show that pull-in voltage is lowered when initial curvature is higher. This is due to the fact that in the initially curved configuration, the beam is closer to the electrode which means that, at the same level of voltage, electrostatic force is higher compared to elastic force.

Figure 4: Variation of the normalized static deflection with DC voltage for different initial curvatures .

Figure 5: Variation of the pull-in voltage with the initial curvatures .

4. Natural Frequencies and Mode Shapes of the Curved CNT under DC Voltage

Assuming that superposition principle applies to the CNT deflection, we split it into static part due to DC voltage and dynamic part due to AC voltage assumed to be harmonic, that is, . For the DQM form, we replace the latter expression by its value at node ; that is, . Then we use (4) and drop the damping and forcing AC terms. The electrostatic force is linearized using the Taylor series to the first order around its static solution [25]. Finally, we obtain the following eigenvalue problem of order :

which can be reduced to the following matrix form:

where . The natural frequencies are the square roots of the eigenvalues of the matrix , and the mode shapes are the corresponding eigenvectors.

4.1. Effect of Initial Curvature on the Natural Frequencies

At a given applied DC voltage, we investigate the effect of variation of initial curvature on natural frequencies described by the parameter . The results are summarized in Figure 6.

Figure 6: Variation of the 9 first natural frequencies with at V.

We can see clearly that even frequencies are insensitive to the variation of curvature, whereas the odd frequencies vary with slack. One can see that odd frequencies do not intersect, but they diverge in a manner called frequency or mode veering [29]. A frequency veering occurs when the loci of two frequencies approach each other when a parameter is varied and then veer away when being too close (similar to the repulsive charges process). Adding to that, we can also clearly see that the odd frequencies intersect the even ones for certain levels of slack offering many possibilities of internal resonances and exchange of energy among higher- and lower-order modes [31].

The same phenomenon of mode veering is observed in Figure 7. This behavior will be further investigated in the next section by exploring the dynamic behavior of the CNT before and after veering.

Figure 7: Variation of the 9 first natural frequencies with DC voltage at nm.

4.2. Effect of DC Voltage on the Natural Frequencies

For initial curvature nm, we vary the applied DC voltage and investigate its effect on the natural frequencies. The results are summarized in Figure 7. It shows that natural frequencies of the CNT increase as function of DC voltage exhibiting a hardening-type behavior. This behavior perseveres until pull-in where the natural frequencies drop down. Also from the same figure, two consecutive odd natural frequencies approach each other as the DC voltage is varied and then diverge; this is known as mode veering.

4.3. CNT Mode Shapes Interchangeability

To further clarify the veering process with the DC voltage depicted in Figure 7, we plot the eigenfunctions (mode shapes) of the CNT natural frequencies for  V, Figure 8. In a frequency veering process, the eigenfunctions associated with the eigenvalues on each locus before veering are interchanged during the veering. For example, in Figure 8, we can clearly see that mode number 5 is the one that has the shape of the fundamental clamped-clamped beam mode along with its fifth mode. Now, and as seen before in Figure 6, the shape of the first mode is transferred into the different odd modes, from the higher to the lower modes as the DC voltage is increased till it is transmitted to the lowest fundamental frequency, which eventually drops to zero at pull-in, a thing very similar to what was reported in [29] when varying the CNT slack level.

Figure 8: Mode shaped of the CNT without the initial curvature for  nm,  V, and 41 grid points.

5. Dynamic Response and Limit-Cycle Solution of the CNT

In this section, we present the dynamic behavior of the CNT by generating their frequency-response curves for different sets of parameters. For this, we have to solve the dynamic equation (4) obtained by the DQM. For grid points and using the boundary conditions and the symmetry of the problem, we end up with ODEs describing the motion of the system. To seek limit-cycle solutions, we assume that the periodic orbits have the same period than the excitation (). The time is discretized along the assumed period, and the ODEs obtained by the DQM are discretized in time using the FDM. In addition, we enforce the condition that the first and last solutions along the orbit are equal; this will satisfy the periodicity condition of the limit-cycle solution.

To construct the limit-cycle solutions and the frequency-response curves of the CNT, we seek periodic orbits of the system with period . In implementing the FDM, we rescale the time so that the period is unity and discretize the orbit by points; thus, . To guarantee periodicity of the orbit, we impose the condition that the first and last points of the orbit (points 1 and ) are identical. We end up with the following system of nonlinear algebraic equation, which is solved using the Newton-Raphson method [25]:

where , , and .

5.1. Frequency-Response Curves of the CNT

The frequency-response curves are obtained by plotting the variation of the maximum deflection of the CNT at its midpoint as function of the excitation frequency. This excitation is obtained by a DC load superimposed to an AC harmonic load around one of the natural frequencies of the CNT.

In Figures 9 and 10, we plot the frequency responses of CNT near the 1st and 9th natural frequencies for initial curvature nm and frequency responses of CNT near 1st and 11th natural frequencies for initial curvature nm. We observe linear-like shape frequency-response curves in all cases due to low applied voltages (V, V). The results are in agreement with the Galerkin method results reported in [29].

Figure 9: Frequency-response curves near the 1st and 9th natural frequencies when nm. Results are shown for V, V, and .

Figure 10: Frequency-response curves of nm slack showing peaks at the 1st and 11th natural frequencies. Results are shown for V, V, and .

5.2. Modes’ Veering

We study modes veering phenomenon by investigating the dynamic behavior of the CNT before and after veering. According to Figure 6, the first veering between the first and third modes occurs at approximately nm, and the second veering occurs between the third and fifth modes at approximately nm.

Concerning the first veering between and modes, we excite the CNT near its fundamental frequency and look at the response simultaneously around the first and third modes. Two peaks in the Frequency-response curve are obtained and showed in Figures 9 and 10. What we wanted to show in Figures 9 and 11 is that the linear dynamic response is significant in the neighborhood of the first natural frequency () and the frequencies that are located on the veering straight line, as predicted from Figure 6. Those frequencies are in the case of 100 nm curvature and in the case of 200 nm curvature. This is actually one consequence of the mode veering where the system energy is transferred from one mode to another through the veering process. To further clarify this point, we will adopt the same previous procedure for the second veering point (between the third and the fifth natural frequencies in Figures 11 and 12) while exciting near the third natural frequency and getting the dynamic responses around the third and fifth vibration modes. In both cases after veering, the amplitude peak at the higher mode is increased compared to that before veering. The results prove that the dynamic amplitude of the CNT is a sort of exchange between two consecutive odd modes laying on the veering line. This will definitely has an effect on the system quality factor, and we may have a sort of energy exchange among the respective modes.

(a)

(b)

Figure 11: Frequency response near the first and third natural frequencies when first mode is excited before veering at (a) nm and after veering (b) nm.

(a)

(b)

Figure 12: Frequency response near the third and fifth natural frequencies when third mode is excited before veering at (a) nm and after veering (b) nm.

6. Conclusions

In this paper, the nonlinear dynamics of a CNT have been investigated using a novel discretization technique: FDM and DQM. In the static part, CNT responses to DC electric load were calculated solving the static equation using DQM. We have studied the static CNT behavior and how it is affected by residual stress as well as initial curvature. The natural frequencies of the CNT have been calculated for different DC voltages and initial curvatures. This study revealed the existence of the phenomenon of modes veering as both DC voltage and initial curvature are varied.

Dynamic analysis was also investigated for the CNT response to harmonic electric load near the natural frequencies. Frequency-response curves have been plotted and validated under different excitation conditions.

Veering of the modes has been analyzed by comparing the dynamic behavior of the CNT before and after veering. Veering of the modes with initial curvature as parameter is found to exchange dynamic properties between modes. In fact, the maximum deflection at resonant frequency of lower mode before veering is higher than that after veering, and maximum deflection at resonant frequency of higher mode increases considerably during veering of the modes. This proves that energy of the system is transferred from one mode to the other during veering.

References

S. Iijima, “Helical microtubules of graphitic carbon,” Nature, vol. 354, pp. 56–58, 1991.

J. H. Bak, Y. D. Kim, S. S. Hong et al., “High-frequency micromechanical resonators from aluminium-carbon nanotube nanolaminates,” Nature Materials, vol. 7, no. 6, pp. 459–463, 2008.

H. W. C. Postma, I. Kozinsky, A. Husain, and M. L. Roukes, “Dynamic range of nanotube- and nanowire-based electromechanical systems,” Applied Physics Letters, vol. 86, no. 22, Article ID 223105, pp. 1–3, 2005.

H. M. Ouakad and M. I. Younis, “Nonlinear dynamics of electrically actuated carbon nanotube resonators,” AMSE Journal of Computational and Nonlinear Dynamics, vol. 5, no. 1, pp. 1–13, 2010.

M. Dequesnes, S. V. Rotkin, and N. R. Aluru, “Parameterization of continuum theories for single wall carbon nanotube switches by molecular dynamics simulations,” Journal of Computational Electronics, vol. 1, no. 3, pp. 313–316, 2002.

M. Dequesnes, Z. Tang, and N. R. Aluru, “Static and dynamic analysis of carbon nanotube-based switches,” AMSE Journal of Engineering Materials and Technology, vol. 126, no. 3, pp. 230–237, 2004.

P. Poncharal, Z. L. Wang, D. Ugarte, and W. A. de Heer, “Electrostatic deflections and electromechanical resonances of carbon nanotubes,” Science, vol. 283, no. 5407, pp. 1513–1516, 1999.

R. Ansari and H. Ramezannezhad, “Nonlocal Timoshenko beam model for the large-amplitude vibrations of embedded multiwalled carbon nanotubes including thermal effects,” Physica E, vol. 43, no. 6, pp. 1171–1178, 2011.

A. Isacsson and J. M. Kinaret, “Parametric resonances in electrostatically interacting carbon nanotube arrays,” Physical Review B, vol. 79, no. 16, Article ID 165418, 2009.

M. A. Hawwa and H. M. Al-Qahtani, “Nonlinear oscillations of a double-walled carbon nanotube,” Computational Materials Science, vol. 48, no. 1, pp. 140–143, 2010.

B. Fang, Y.-X. Zhen, C.-P. Zhang, and Y. Tang, “Nonlinear vibration analysis of double-walled carbon nanotubes based on nonlocal elasticity theory,” Applied Mathematical Modelling, vol. 37, no. 3, pp. 1096–1107, 2012.

A. Hajnayeb and S. E. Khadem, “Nonlinear vibration and stability analysis of a double-walled carbon nanotube under electrostatic actuation,” Journal of Sound and Vibration, vol. 331, no. 10, pp. 2443–2456, 2012.

P.-H. Huang, “Atomistic simulations of shearing friction and dynamic adhesion of double-walled carbon nanotubes on Au substrates,” Composites Science and Technology, vol. 72, no. 5, pp. 599–607, 2012.

A. R. Hall, S. Paulson, T. Cui, J. P. Lu, L. C. Qin, and S. Washburn, “Torsional electromechanical systems based on carbon nanotubes,” Reports on Progress in Physics, vol. 75, no. 11, Article ID 116501, 2012.

P. Boissy, J. Genest, J. Patenaude et al., “Carbon nanotubes (CNTs) based strain sensors for a wearable monitoring and biofeedback system for pressure ulcer prevention and rehabilitation,” in Proceedings of the 33rd Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '11), pp. 5824–5827, IEEE, Boston, Mass, USA, September 2011.

C. Boero, J. Olivo, S. Carrara, and G. de Micheli, “A self-contained system with CNTs-based biosensors for cell culture monitoring,” IEEE Journal on Emerging and Selected Topics in Circuits and Systems, vol. 2, no. 4, pp. 658–671, 2012.

H. Lorenz, J. Fritzsche, A. Das, et al., “Advanced elastomer nano-composites based on CNT-hybrid filler systems,” Composites Science and Technology, vol. 69, no. 13, pp. 2135–2143, 2009.

Z. Dong, U. C. Wejinya, and S. N. S. Chalamalasetty, “Development of CNT-ISFET based pH sensing system using atomic force microscopy,” Sensors and Actuators A, vol. 173, no. 1, pp. 293–301, 2012.

H. Cho, B. Jeong, M.-F. Yu, A. F. Vakakis, D. M. McFarland, and L. A. Bergman, “Nonlinear hardening and softening resonances in micromechanical cantilever-nanotube systems originated from nanoscale geometric nonlinearities,” International Journal of Solids and Structures, vol. 49, no. 15-16, pp. 2059–2065, 2012.

M. W. Graham, T. R. Calhoun, A. A. Green, M. C. Hersam, and G. R. Fleming, “Two-dimensional electronic spectroscopy reveals the dynamics of phonon-mediated excitation pathways in semiconducting single-walled carbon nanotubes,” Nano Letters, vol. 12, no. 2, pp. 813–819, 2012.

H. B. Peng, C. W. Chang, S. Aloni, T. D. Yuzvinsky, and A. Zettl, “Microwave electromechanical resonator consisting of clamped carbon nanotubes in an abacus arrangement,” Physical Review B, vol. 76, no. 3, Article ID 035405, 2007.

P. A. Greaney and J. C. Grossman, “Nanomechanical energy transfer and resonance effects in single-walled carbon nanotubes,” Physical Review Letters, vol. 98, no. 12, Article ID 125503, 2007.

F. Najar, S. Choura, S. El-Borgi, E. M. Abdel-Rahman, and A. H. Nayfeh, “Modeling and design of variable-geometry electrostatic microactuators,” Journal of Micromechanics and Microengineering, vol. 15, no. 3, pp. 419–429, 2005.

F. Najar, S. Choura, E. M. Abdel-Rahman, S. El-Borgi, and A. Nayfeh, “Dynamic analysis of variable-geometry electrostatic microactuators,” Journal of Micromechanics and Microengineering, vol. 16, no. 11, pp. 2449–2457, 2006.

F. Najar, A. H. Nayfeh, E. M. Abdel-Rahman, S. Choura, and S. El-Borgi, “Nonlinear analysis of MEMS electrostatic microactuators: primary and secondary resonances of the first mode,” Journal of Vibration and Control, vol. 16, no. 9, pp. 1321–1349, 2010.

V. Sazonova, Y. Yalsh, I. Üstünel, D. Roundy, T. A. Arlas, and P. L. McEuen, “A tunable carbon nanotube electrochemical oscillator,” Nature, vol. 431, no. 7006, pp. 284–287, 2004.

M.-F. Yu, “Fundamental mechanical properties of carbon nanotubes: current understanding and the related experimental studies,” ASME Journal of Engineering Materials and Technology, vol. 126, no. 3, pp. 271–278, 2004.

H. M. Ouakad and M. I. Younis, “Dynamic response of slacked single-walled carbon nanotube resonators,” Nonlinear Dynamics, vol. 67, no. 2, pp. 1419–1436, 2012.

H. M. Ouakad and M. I. Younis, “Natural frequencies and mode shapes of initially curved carbon nanotube resonators under electric excitation,” Journal of Sound and Vibration, vol. 330, no. 13, pp. 3182–3195, 2011.

A. Bhushan, M. M. Inamdar, and D. N. Pawaskar, “Investigation of the internal stress effects on static and dynamic characteristics of an electrostatically actuated beam for MEMS and NEMS application,” Microsystem Technologies, vol. 17, no. 12, pp. 1779–1789, 2011.

H. Ouakad and M. I. Younis, “Internal resonance in electrically actuated MEMS arches,” in Proceeding of the 7th European Nonlinear Dynamics Conference (ENOC '11), Rome, Italy, 2011.