Static and free vibration analysis of carbon nano wires based on Timoshenko beam theory using differential quadrature method

ARTICLES

Static and free vibration analysis of carbon nano wires based on Timoshenko beam theory using differential quadrature method

Maziar Janghorban^*

Department of Mechanical Engineering, Shiraz Branch, Islamic Azad University, Shiraz, Iran

ABSTRACT

Static and free vibration analysis of carbon nano wires with rectangular cross section based on Timoshenko beam theory is studied in this research. Differential quadrature method (DQM) is employed to solve the governing equations. From the knowledge of author, it is the first time that free vibration of nano wires is investigated. It is also the first time that differential quadrature method is used for bending analysis of nano wires.

Keywords: carbon nano wires, Timoshenko beam theory, differential quadrature method, free vibration, static analysis.

1 INTRODUCTION

A nano wire (NW) is an extremely thin wire with a diameter on the order of a few nanometers (nm) or less, where 1 nm = 10 ^-9 meters. Nano wires have many interesting properties that are not seen in bulk or 3-D materials. Many different types of nano wires exist, including metallic (e.g., Ni, Pt, Au), semiconducting (e.g., InP, Si, GaN, etc.), and insulating (e.g., SiO₂, TiO₂). Nano structures such as single and double walled carbon nanotubes (CNTs) [3], nano plates [10, 12, 14] and nano rods [1] are widely studied in past decade but static and dynamic analysis of nano wires are rarely investigated. Surface effects on the elastic behavior of static bending nano wires were studied by using a comprehensive Timoshenko beam model by Jiang and Yan [9]. Fu and Zhang [7] were established a continuum elastic model for core-shell nano wires with weak interfacial bonding. Critical buckling loads and resonant frequencies of simply supported nano wires are obtained by using the Ritz method. Song and Huang [19] presented a model of surface stress effects on bending behavior of nano wires based on the incremental deformation theory. Differential quadrature (DQ) was introduced by Bellman et al [2] as a simple and highly efficient numerical technique in 1971. This method generalized and simplified further by Quan and Chang [15] and Shu and Richards [18] through explicit expressions of the weighting coefficients associated with their derivatives. This method was used several times for investigation of micro and nano structures. Civalek et al [5] applied the differential quadrature method to the equations of motion and bending of Euler-Bernoulli beam using the nonlocal elasticity theory for cantilever microtubules. Wang et al [22] concerned with the use of the Timoshenko beam model for free vibration analysis of multi-walled carbon nanotubes by using differential quadrature method. Pradhan and Murmu [13] developed a single nonlocal beam model and applied to investigate the flapwise bending-vibration characteristics of a rotating nano cantilever using DQ method. Murmu and Pradhan [11] also presented the buckling analysis of a single-walled carbon nanotube embedded in elastic medium based on nonlocal elasticity and Timoshenko beam theory and using DQM. In this paper, static and free vibration analysis of carbon nano wires based on Timoshenko beam theory using differential quadrature method is investigated. From the knowledge of author, two kinds of nano wire are exits, one with rectangular cross section and the other with circular cross section. In this study, nano wires with rectangular cross sections are studied.

2 DIFFERENTIAL QUADRATURE METHOD

The numerical solution of partial differential equations plays a considerable role in the areas of engineering. The differential quadrature method is one which can satisfy the above purpose. Differential quadrature is characterized by approximating the derivatives of a function using a weighted sum of functions values on a set of selected discrete points. In order to show the mathematical representations of DQM, consider a function f(x, y) having its field on a rectangular domain 0 < x < a and 0 < y < b. let, in the given domain, the function values be known or desired on a grid of sampling points. According to DQ method, the r^th derivative of a function f(x, y) can be approximated as [2],

where i = 1, 2, ..., N_x , j = 1, 2, ..., N_y and r = 1, 2, ..., N_x - 1.

From this equation one can deduce that the important components of DQ approximations are weighting coefficients and the choice of sampling points. In order to determine the weighting coefficients a set of test functions should be used in Eq. (1). For polynomial basis functions DQ, a set of Lagrange polynomials are employed as the test functions. The weighting coefficients for the first-order derivatives in ξ-direction are thus determined as

where L_x is the length of domain along the x-direction and M(x_i) = (x_i - x_k) .

The weighting coefficients of second order derivative can be obtained as,

In a similar way, the weighting coefficients for y-direction can be obtained. The weighting coefficient of the third and fourth order derivatives (C_ij, D_ij) can be computed easily from (B_ij) by

3 GOVERNING EQUATIONS

The importance of shear deformation and rotary inertia in the description of the dynamic response of beams is well documented and an improved theory was given by Timoshenko already in 1921 [20]. Consider a constant thickness nano wire with rectangular cross section with length L, thickness h and width b. In order to calculate the deflections and natural frequencies of carbon nano wires, Timoshenko beam theory is used. The Timoshenko model is an extension of the Euler-Bernoulli model by taking into account two additional effects: shearing force effect and rotary motion effect. For Timoshenko beam theory, the displacement field is assumed to be as follow,

The nonzero strains of the Timoshenko beam theory are [16],

The equations for a Timoshenko beam theory are given by,

Where W is the deflection of beam, θ is the rotation of the normal line. We also have,

Assuming the sinusoidal motion in time, Eq. (4) yields

Contrary to most other publications, here, these two differential equations of second order are used instead of rearranging these to one differential equation of fourth order. Using the DQ-rules for the spatial derivatives, the DQ-analogs of the governing Eqs. (9) become

Based on Eq. (7) a computer code is developed in MATLAB for static and free vibration analysis. For free vibration analysis, the DQ discretized form of the equations of motion and the related boundary conditions can be expressed in the matrix form as,

In the above equations, K_dd, K_db, K_bd and K_bb refers to domain-domain, domain-boundary, boundary-domain and boundary-boundary stiffness matrix, respectively. The mass matrix [M] is obtained from equations of motion and those of the stiffness matrices [K_ii] (i = b, d) are obtained from the domain and boundary conditions. Solving the eigenvalue system of below equation, the natural frequencies will be obtained.

Where [KK] = [K_dd] - [K_db][K_bb]^-1[K_bd]. For static analysis, deflections can easily calculated from the equation below,

where {F} is the force vector. It is mentioned that in this case, the matrix [KK] can be written directly and we do not need to use [K_ii] (i = b, d).

The boundary condition used in this study is defined as,

The discretized form of boundary condition can be obtained by,

4 NUMERICAL RESULTS

Let us consider a carbon nano wire having a length L = 400 nm, thickness h=40 nm, width b= 2h, Young modulus E = 76 GPa [9], density ρ = 1220 Kg/m³ and shear correction factor for rectangular cross section k = 5(1 + v)/(6 + 5v) [9]. The moments of inertia and the areas are calculated using I = (1/12)(bh³) and A = bh. The numerical results can be shown in the dimensionless quantities defined as follow:

In order to show the accuracy of present work, our results are compared with other works for single walled carbon nano tubes and then the results for carbon nano wires are presented. For static analysis, In Table 1, the results are compared with those of other numerical solutions for a simply supported single walled carbon nanotube. The convergency of the method is also investigated, too. One can easily find that even 12 grid points is enough to have accurate results. For free vibration analysis, to validate the present formulation and the computer program developed by the author, using the DQ method, our natural frequencies are compared with Reddy [16], Heireche et al [8] and Demir et al [6] in Table 2. Very good agreement is achieved between the results of present work and other methodologies. For carbon nano wires, the effects of width to thickness, length to width and length to thickness ratios on the deflections of carbon nano wires subjected to uniformly distributed loading are shown in Figs. 2, 3 and 4. One can see that increasing the length to thickness, length to width and thickness to width ratios will increase the deflections of carbon nano wires. The shapes of carbon nano wires under mechanical loading are as we expected and the maximum deflections appear at the middle of nanowire. The convergency of free vibration analysis is studied in Fig. 5. It seems that only nine grid points are enough for investigation of carbon nano wires. It is also shown that for higher modes, more grid points are needed.

From the Tables 1 and 2 and Fig. 5, it may be observed that the differential quadrature method is one of the fastest numerical tools for solving governing equations. The effects of width to thickness and length to thickness ratios on natural frequencies are investigated in Figs. 6 and 7. It is shown that for increasing the natural frequencies of carbon nano wires, one may increase the length to thickness and thickness to width ratios. The increase in the values of the deflections and natural frequencies as the width to thickness and length to thickness ratios are varied is evident of its importance contribution to the mathematical model.

5 CONCLUSION

Deflections and natural frequencies of carbon nano wires with constant thickness have been obtained by using a Timoshenko beam theory. The differential quadrature method was employed to convert the governing differential equations into a linear system. The accuracy and convergency of present work were investigated. The effects of different parameters are studied, too. Some of the results of present study is,

1. Increasing the length to thickness, length to width and thickness to width ratios will increase the deflections of carbon nano wires.

2. Increase the length to thickness and thickness to width ratios will increase the deflections of carbon nano wires

Received 30 May 2011;

In revised form 4 Sep 2011

• [1] M. Aydogdu. Axial vibration of the nano rods with nonlocal continuum rod model. Physica E, 41:8614, 2009.
• [2] R.E. Bellman and J. Casti. Differential quadrature and long term integration. J. Math. Anal. Appl., 34:235238, 1971.
• [3] Ö. Civalek and B. Akgööz. Static analysis of single walled carbon nanotubes (SWCNT) based on Eringen's nonlocal elasticity theory. Int. J. Eng. Appl. Sciences, 1:4756, 2009.
• [4] Ö. Civalek and C. Demir. Bending analysis of microtubules using nonlocal Euler-Bernoulli beam theory. Applied Mathematical Modelling, 2010. In Press.
• [5] Ö. Civalek, C. Demir, and B. Akgööz. Free vibration and bending analyses of cantilever microtubules based on nonlocal continuum model. Math. Comput. Appl., 15:289298, 2010.
• [6] Ç. Demir, Ö. Civalek, and B. Akgööz. Free vibration analysis of carbon nanotubes based on shear deformable beam theory by discrete singular convolution technique. Math. Comput. Appl., 15:5765, 2010.
• [7] Y. Fu and P. Zhang. Buckling and vibration of core-shell nano wires with weak interfaces. Mechanics Research Communications, 37:622626, 2010.
• [8] H. Heireche, A. Tounsi, A. Benzair, M. Maachou, and E.A.D. Bedia. Sound wave propagation in single-walled carbon nanotubes using nonlocal elasticity. Physica E, 40:27912799, 2008.
• [9] L.Y. Jiang and Z. Yan. Timoshenko beam model for static bending of nano wires with surface effects. Physica E, 42:22742279, 2010.
• [10] E. Jomehzadeh and A.R. Saidi. Decoupling the nonlocal elasticity equations for three dimensional vibration analysis of nano-plates. Compos Struct, 93(2):101520, 2010.
• [11] T. Murmu and S.C. Pradhan. Buckling analysis of a single-walled carbon nanotube embedded in an elastic medium based on nonlocal elasticity and Timoshenko beam theory and using DQM. Physica E, 41:12321239, 2009.
• [12] T. Murmu and S.C. Pradhan. Vibration analysis of nano plates under uniaxial pretressed conditions via nonlocal elasticity. J Appl Phys, 106:104301, 2009.
• [13] S.C. Pradhan and T. Murmu. Application of nonlocal elasticity and DQM in the flapwise bending vibration of a rotating nano cantilever. Physica E, 42:19441949, 2010.
• [14] S.C. Pradhan and J.K. Phadikar. Nonlocal elasticity theory for vibration of nano plates. J Sound Vib., 325:20623, 2009.
• [15] J.R. Quan and C.T. Chang. New insights in solving distributed system equations by the quadrature method I: Analysis. Comput. Chem. Eng., 13:779788, 1989.
• [16] J.N. Reddy. Nonlocal theories for bending, buckling and vibration of beams. Int. J. Eng. Sciences, 45:288307, 2007.
• [17] J.N. Reddy and S.D. Pang. Nonlocal continuum theories of beams for the analysis of carbon nanotubes. J. Appl. Phys., 103:023511, 2008.
• [18] C. Shu and B.E. Richards. Application of generalized differential quadrature to solve two-dimensional incompressible Navier-Stokes equations. Int. J. Numer. Meth. Fluids, 15:791798, 1992.
• [19] F. Song and G.L. Huang. Modeling of surface stress effects on bending behavior of nano wires: Incremental deformation theory. Physics Letters A, 373:39693973, 2009.
• [20] S.P. Timoshenko. On the correction for shear of the differential equation for transverse vibrations of prismatic bars. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 41(6):744746, 1921.
• [21] ucsdnews.ucsd.edu
  » link
• [22] C.M. Wang, V.B.C. Tan, and Y.Y. Zhang. Timoshenko beam model for vibration analysis of multi-walled carbon nanotubes. Journal of Sound and Vibration, 294:10601072, 2006.

Publication Dates

• Publication in this collection
  17 Apr 2012
• Date of issue
  2011

History

• Received
  30 May 2011
• Accepted
  04 Sept 2011