Next Article in Journal
Highly Absorbent Antibacterial Hemostatic Dressing for Healing Severe Hemorrhagic Wounds
Previous Article in Journal
A Multi-Level Decision Fusion Strategy for Condition Based Maintenance of Composite Structures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Imperfection Sensitivity of Nonlinear Vibration of Curved Single-Walled Carbon Nanotubes Based on Nonlocal Timoshenko Beam Theory

by
Iman Eshraghi
1,
Seyed K. Jalali
2,* and
Nicola Maria Pugno
3,4,5
1
School of Mechanical Engineering, University of Tehran, Tehran 1193653471, Iran
2
Department of Mechanical Engineering, Kermanshah University of Technology, Kermanshah 6715685438, Iran
3
Laboratory of Bio-Inspired & Graphene Nanomechanics, Department of Civil, Environmental and Mechanical Engineering, Università di Trento, via Mesiano, 77, Trento 38123, Italy
4
Center for Materials and Microsystems, Fondazione Bruno Kessler—via Sommarive 18, Povo, Trento 38123, Italy
5
School of Engineering & Materials Science, Queen Mary University of London—Mile End Road, London E1 4NS, UK
*
Author to whom correspondence should be addressed.
Materials 2016, 9(9), 786; https://doi.org/10.3390/ma9090786
Submission received: 2 August 2016 / Revised: 29 August 2016 / Accepted: 6 September 2016 / Published: 21 September 2016
(This article belongs to the Section Advanced Materials Characterization)

Abstract

:
Imperfection sensitivity of large amplitude vibration of curved single-walled carbon nanotubes (SWCNTs) is considered in this study. The SWCNT is modeled as a Timoshenko nano-beam and its curved shape is included as an initial geometric imperfection term in the displacement field. Geometric nonlinearities of von Kármán type and nonlocal elasticity theory of Eringen are employed to derive governing equations of motion. Spatial discretization of governing equations and associated boundary conditions is performed using differential quadrature (DQ) method and the corresponding nonlinear eigenvalue problem is iteratively solved. Effects of amplitude and location of the geometric imperfection, and the nonlocal small-scale parameter on the nonlinear frequency for various boundary conditions are investigated. The results show that the geometric imperfection and non-locality play a significant role in the nonlinear vibration characteristics of curved SWCNTs.

Graphical Abstract

1. Introduction

Nanostructures have attracted more attention in recent years due to their superior mechanical, electrical, and chemical characteristics [1]. These properties have led to widespread applications of these structures in nano-electro-mechanical-system (NEMS) devices. This growth in development and use of nanostructures necessitates designing, fabricating and characterizing their mechanical behavior. The literature reveals that carbon nanotubes (CNTs) have the potential to be used in these devices. Many of the NEMS devices like oscillators and mass measuring sensors use the vibration characteristics of the CNTs [2,3,4]. This has led to a remarkable number of studies conducted on the characterization of vibration behavior of CNTs [5,6,7,8,9,10]. A review including modeling, simulation and application of vibration behavior of CNTs can be found in [11].
Various methods are utilized to model and analyze the behavior of nano-beams. Among them, nonlocal elasticity of Eringen [12], as a non-classical theory of elasticity which considers the scale effect, has been widely used for mathematical continuum modeling and analysis of bending, buckling and vibration of nanostructures. In this way, known classical theories of beams like Euler-Bernoulli or Timoshenko beam theory can be adopted to model the displacement field of a CNT and nonlocal theory may be employed to modify constitutive stress-strain relationships to take small-scale effect into account. There has been extensive study on modeling and analysis of CNTs using this theory [13,14,15]. Wang and Varadan [16] investigated the free vibration of both single-walled and double-walled nanotubes via nonlocal elasticity. The small-scale effects on the vibration characteristics of carbon nanotubes were explicitly derived in their work, and effects of the length and diameter scale on the vibrational frequency were revealed. Another study by Lu et al. [17] investigated application of nonlocal Euler and Timoshenko beam models for vibration and wave propagation of CNTs and provided the effects of small-scale parameter on the free vibration frequency. Murmu and Pradhan [18] studied thermo-mechanical vibration of a SWCNT embedded in an elastic medium based on nonlocal elasticity theory. They used the DQ method to numerically find thermal vibration response of CNTs. Furthermore, influences of nonlocal small-scale effects, temperature change and elastic medium constant on the frequency values of various modes were investigated. The free vibration response of SWCNTs using the nonlocal elasticity based on Euler–Bernoulli, Timoshenko and Reddy beam theory for various boundary conditions was investigated by Ansari and Sahmani [19]. They obtained fundamental frequencies of nanotubes with different chiralities by implementing molecular dynamics (MD) simulation and comparing the output of their simulations with the results obtained by the nonlocal beam models to identify appropriate values of nonlocal parameter for various types of chirality and boundary conditions.
In spite of significance of nonlinear vibrational response of CNTs, the number of studies which consider large amplitude vibrations are limited. Reddy [20,21] used nonlocal elasticity theory and von Kármán nonlinear strain-displacement relations to reformulate beam and plate theories for analysis of bending, buckling and vibration behavior of these nanostructures. Yang et al. [22] developed nonlinear free vibration equations of SWCNTs using nonlocal Timoshenko beam theory and used DQ procedure to discretize the governing equations. They obtained elastic modulus through molecular mechanics simulation. The influences of nonlocal parameter, length and radius of the nano-beam and type of supports on the nonlinear free vibration characteristics were investigated by using an iterative procedure. Nonlinear vibration analysis of double-walled carbon nanotubes based on nonlocal elasticity theory considering the von Kármán type geometric nonlinearity was studied by Fang et al. [23]. It was shown that nonlocal parameter, aspect ratio and surrounding elastic medium have more effect on the nonlinear non-coaxial vibration than the coaxial vibration of the nanotubes.
All the aforesaid works modeled CNTs as straight beams. However, an examination of the studies concerned with fabrication procedures of CNTs reveals that perfectly straight CNTs without any deviation in geometry and orientation are very difficult to achieve [5]. In a historical view of classical beam theories, these deviations could be modeled as an initial geometric imperfection. Initial geometric imperfection refers to globally or locally unavoidable deviations between the actual shape and desired shape of the structure [24]. It is reported that such imperfections have remarkable effects on nonlinear vibration frequencies of structures [25,26]. There has been quite a few works considering the effects of geometrical imperfection in continuum modeling of CNTs. Farshidianfar and Soltani [27] based on the nonlocal continuum theory analyzed transverse vibration of a curved single-walled carbon nanotube conveying a fluid. They discussed the effects of the flow velocity, nonlocal parameter, the stiffness of the elastic foundation, and the boundary conditions in detail. Ouakad and Younis [28] reported the natural frequencies and mode shapes of initially curved carbon nanotube resonators under electric excitation. The variation of natural frequencies and mode shapes with amplitude of imperfection and the DC electrostatic load using a multimode Galerkin procedure were investigated.
The aim of this paper is to compute the imperfection sensitivity of large amplitude vibration of curved SWCNTs based on a nonlocal Timoshenko beam model. To include geometric nonlinearities, von Kármán type of strain-displacement relation is adopted for small strain and moderate rotations of the beam. An iterative method [22] is employed to solve the discretized eigenvalue equation of the free vibration of the beam. Extensive numerical results are presented to investigate the influence of amplitude and configuration of the geometric imperfection, boundary conditions and the nonlocal small-scale parameter on vibration behavior of a curved SWCNT.

2. Problem Formulation

2.1. Timoshenko Beam Model for Vibration of Curved Single-Walled Carbon Nanotubes (SWCNTs)

The curved SWCNT as a Timoshenko beam model with length L, constant cross section, A , mean wall radius r and thickness h, is illustrated in Figure 1. w ¯ ( x ) denotes the typical initial geometric imperfection of the curved SWCNT.
By using Timoshenko beam theory, displacement components along x- and z- axes are given by
u ( x , z , t ) = u 0 ( x , t ) + z ϕ ( x , t ) w ( x , z , t ) = w 0 ( x , t ) + w ¯ ( x )
where u 0 ( x , t ) and w 0 ( x , t ) are longitudinal and transversal displacements of the straight beam axis, ϕ ( x , t ) , denotes rotation of the cross-section with respect to the straight beam axis and t represents time. For the case of large amplitude vibration, based on nonlinear strain-displacement relations of von Kármán type, the strain components at any point are expressed as [29]:
ε x x = u , x 0 + 1 2 ( w , x 0 ) 2 + w , x 0 w ¯ , x + z ϕ , x γ x z = w , x 0 + ϕ
where ε x x is the normal strain, and γ x z is the shear strain. Note that , x represents differentiation with respect to x. Using the principle of virtual displacements [20] and taking into account geometric initial imperfection, the Euler-Lagrange equations for the beam are given by
N , x + f = I 0 u ¨ 0 Q , x + [ N ( w , x 0 + w ¯ , x ) ] , x + q = I 0 w ¨ 0 M , x Q = I 2 ϕ ¨
And the corresponding immovable boundary conditions at the ends of the beam (x = 0, L) are
u 0 = 0 ,   w 0 = 0 ,   ϕ = 0 Clamped End u 0 = 0 ,   w 0 = 0 ,   M = 0 Hinged End
In Equation (3) a dot above each quantity represents temporal differentiation. f and q are externally applied distributed forces (per length unit) on the beam in the x- and z- directions, respectively. I 0 and I 2 are mass moments of inertia of the beam section given by
I 0 = A ρ ( x , z ) d A I 2 = A z 2 ρ ( x , z )   d A
where ρ ( x , z ) is the beam material density. The resultant moment M, normal force N, and shear force Q can be calculated by integrating the corresponding axial ( σ ) and in-plane shear stress ( τ ) components on the cross-section area as
N = A σ x x   d A M = A z   σ x x   d A Q = A τ x z   d A

2.2. Nonlocal Theory of Nano-Beams

Based on Eringen’s nonlocal elasticity theory, the state of stress at a point in a material is a function of strain field at every point in the continuum [12]. The integral form of the nonlocal theory that defines stress in each point in the material is stated as
σ = V χ ( | Y X | , λ ) S ( Y ) d Y
where X is the location of the reference point and χ ( | Y X | , τ ) is a function that imposes the effect of strain at point Y into the nonlocal constitutive equation. | Y X | denotes the distance between two points and λ is the parameter for small-scale factor. S ( Y ) is the classical stress tensor that is related to local strain field for a Hookean solid using fourth-order elastic moduli tensor C ( X ) according to the following inner product operation of tensors
S ( X ) = C ( X ) : ε ( X )
Equivalent differential form of Equation (7) can be represented by [20]:
( σ ) = S   ,   = 1 μ ˜ 2 2   and   μ ˜ = e 0 a 0
where e 0 is a material constant and a 0 is an internal characteristic length [30] and 2 is the Laplacian operator which is equal to d 2 / d x 2 for a one-dimensional case. Thus the nonlocal constitutive relations of a nano-beam can be written as
σ xx μ ˜ 2 σ , xx xx = E ε xx τ xz μ ˜ 2 τ , xx xz = G γ xz
where E and G are the Young’s modulus and shear modulus of the nano-beam, respectively.

2.3. Derivation of Nonlocal Governing Equations of Motion

Now, integrating Equation (10) on the cross-section area with the assumption of constant material properties and using Equations (2) and (6) for the definition of strain components and resultant forces, yields
N μ ˜ 2 N , x x = E A ( u , x 0 + 1 2 ( w , x 0 ) 2 + w , x 0 w ¯ , x ) M μ ˜ 2 M , x x = E ϕ , x Q μ ˜ 2 Q , x x = K s G A ( w , x 0 + ϕ )
where is the area moment of inertia of the beam cross-section ( = A z 2 d A ) and K s is the shear correction factor. Substituting Equation (3) into Equation (11) and neglecting distributed external forces for the problem of free vibration, stress resultants can be obtained in terms of displacement components as:
N = E A ( u , x 0 + 1 2 ( w , x 0 ) 2 + w , x 0 w ¯ , x ) + μ ˜ 2 I 0 u ¨ , x 0 M = E ϕ , x + μ ˜ 2 { I 0 w ¨ 0 [ N ( w , x 0 + w ¯ , x ) ] , x + I 2 ϕ ¨ , x } Q = K s G A ( w , x 0 + ϕ ) + μ ˜ 2 { I 0 w ¨ , x 0 [ N ( w , x 0 + w ¯ , x ) ] , x x }
Thus, the equations of motion for the nonlocal Timoshenko beam considering nonlinear strain-displacement relations and initial geometric imperfection can be obtained by back substitution of Equation (12) into Equation (3).
E A ( u , x x 0 + w , x 0 w , x x 0 + w , x 0 w ¯ , x x 0 + w , x x 0 w ¯ , x 0 ) = I 0 ( u ¨ 0 μ u ¨ , x x 0 ) K s G A ( w , x x 0 + ϕ , x ) + Λ 1 μ ˜ 2 Λ 2 = I 0 ( w ¨ 0 μ ˜ 2 w ¨ , x x 0 ) I 0 μ ˜ 2 { u ¨ , x x 0 ( w , x 0 + w ¯ , x ) + u ¨ , x 0 ( w , x x 0 + w ¯ , x x ) μ ˜ 2 [ u ¨ , x x x x 0 ( w , x 0 + w ¯ , x ) + u ¨ , x x x 0 ( w , x x 0 + w ¯ , x x ) + u ¨ , x x 0 ( w , x x x 0 + w ¯ , x x x ) + u ¨ , x 0 ( w , x x x x 0 + w ¯ , x x x x ) ] } E ϕ , x x K s G A ( w , x 0 + ϕ ) = I 2 ( ϕ ¨ μ ˜ 2 ϕ ¨ , x x )
where
Λ 1 = E A [ u , xx 0 w ¯ , x + u , x 0 w ¯ , xx + w , xx 0 ( u , x 0 + 3 2 w , x 0 w , x 0 + 3 w , x 0 w ¯ , x + w ¯ , x w ¯ , x ) + w , x 0 ( u , xx 0 + 3 2 w , x 0 w ¯ , xx + 2 w ¯ , xx w ¯ , x ) ] Λ 2 = E A [ u , x x x x 0 w ¯ , x + 3 u , x x x 0 w ¯ , x x + 3 u , x x 0 w ¯ , x x x + u , x 0 w ¯ , x x x x + w , x x x x 0 ( u , x 0 + 3 2 w , x 0 w , x 0 + 3 w , x 0 w ¯ , x + w ¯ , x w ¯ , x ) + w , x x x 0 ( 3 u , x x 0 + 9 w , x x 0 w , x 0 + 9 w , x x 0 w ¯ , x + 9 w , x 0 w ¯ , x x + 6 w ¯ , x x w ¯ , x ) + w , x x 0 ( 3 u , x x x 0 + 3 w , x x 0 w , x x 0 + 9 w , x x 0 w ¯ , x x + 9 w , x 0 w ¯ , x x x + 6 w ¯ , x x w ¯ , x x + 6 w ¯ , x x x w ¯ , x ) + w , x 0 ( u , x x x x 0 + 3 2 w , x 0 w ¯ , x x x x + 2 w ¯ , x x x x w ¯ , x + 6 w ¯ , x x x w ¯ , x x ) ]
The following non-dimensional parameters are introduced to make equations of motion dimensionless:
ξ = x L ,   α = h L ,   ( u * , w * , w ¯ * ) = ( u 0 , w 0 , w ¯ ) h ,   ϕ = ϕ ,   ( I 0 * , I 2 * ) = ( I 0 I 0 , I 2 I 0 h 2 ) ,   ( D 0 * , D 2 * , G * ) = ( E A E A , E E A h 2 , K s G A E A ) ,   μ = μ ˜ L , t * = t L E A I 0
Noting that now a dot superscript denotes differentiation with respect to non-dimensional time ( t * ) and , ξ represents differentiation with respect to non-dimensional parameter, ξ . μ is called the nonlocal parameter. Thus, Equation (13) and the boundary conditions, Equation (4), can be rewritten in the following form.
D 0 * [ u , ξ ξ * + α ( w , ξ * w , ξ ξ * + w , ξ * w ¯ , ξ ξ * + w , ξ ξ * w ¯ , ξ * ) ] = I 0 * ( u ¨ * μ * u ¨ , ξ ξ * ) G * ( w , ξ ξ * + 1 α ϕ , ξ ) + Λ 1 * μ 2 Λ 2 * = I 0 * ( w ¨ * μ 2 w ¨ , ξ ξ * ) I 0 * μ 2 α { u ¨ , ξ ξ * ( w , ξ * + w ¯ , ξ * ) + u ¨ , ξ * ( w , ξ ξ * + w ¯ , ξ ξ * ) μ 2 [ u ¨ , ξ ξ ξ ξ * ( w , ξ * + w ¯ , ξ * ) + u ¨ , ξ ξ ξ * ( w , ξ ξ * + w ¯ , ξ ξ * ) + u ¨ , ξ ξ * ( w , ξ ξ ξ * + w ¯ , ξ ξ ξ * ) + u ¨ , ξ * ( w , ξ ξ ξ ξ * + w ¯ , ξ ξ ξ ξ * ) ] } D 2 * ϕ , ξ ξ 1 α G * ( w , ξ * + 1 α ϕ ) = I 2 * ( Φ ¨ μ 2 ϕ ¨ , ξ ξ )
Λ 1 * = D 0 * α   { u , ξ ξ * w ¯ , ξ * + u , ξ * w ¯ , ξ ξ * + w , ξ ξ * [ u , ξ * + α ( 3 2 w , ξ * w , ξ * + 3 w , ξ * w ¯ , ξ ξ * + w ¯ , ξ * w ¯ , ξ * ) ] + w , ξ * [ u , ξ ξ * + α ( 3 2 w , ξ * w ¯ , ξ ξ * + 2 w ¯ , ξ ξ * w ¯ , ξ * ) ] } Λ 2 * = D 0 * α   { u , ξ ξ ξ ξ * w ¯ , ξ * + 3 u , ξ ξ ξ * w ¯ , ξ ξ * + 3 u , ξ ξ * w ¯ , ξ ξ ξ * + u , ξ * w ¯ , ξ ξ ξ ξ * + w , ξ ξ ξ ξ * [ u , ξ * + α ( 3 2 w , ξ * w , ξ * + 3 w , ξ * w ¯ , ξ * + w ¯ , ξ * w ¯ , ξ * ) ] + w , ξ ξ ξ * [ 3 u , ξ ξ * + α ( 9 w , ξ ξ * w , ξ * + 9 w , ξ ξ * w ¯ , ξ * + 9 w , ξ * w ¯ , ξ ξ * + 6 w ¯ , ξ ξ * w ¯ , ξ * ) ] + w , ξ ξ * [ 3 u , ξ ξ ξ * + α ( 3 w , ξ ξ * w , ξ ξ * + 9 w , ξ ξ * w ¯ , ξ ξ * + 9 w , ξ * w ¯ , ξ ξ ξ * + 6 w ¯ , ξ ξ * w ¯ , ξ ξ * + 6 w ¯ , ξ ξ ξ * w ¯ , ξ * ) ] + w , ξ * [ u , ξ ξ ξ ξ * + α ( 3 2 w , ξ * w ¯ , ξ ξ ξ ξ * + 2 w ¯ , ξ ξ ξ ξ * w ¯ , ξ * + 6 w ¯ , ξ ξ ξ * w ¯ , ξ ξ * ) ] }
Clamped end:
u * = w * = ϕ = 0
Hinged end:
u * = w * = 0 M = D 2 * Φ , ξ D 0 * μ 2 [ w ¯ , ξ * u , ξ ξ * + w ¯ , ξ ξ * u , ξ * + w , ξ * ( u , ξ ξ * + 3 2 α   w , ξ * w ¯ , ξ ξ * + 2 α w ¯ , ξ ξ * w ¯ , ξ * ) + w , ξ ξ * ( u , ξ * + 3 2 α   w , ξ * w , ξ * + 3 α   w , ξ * w ¯ , ξ * + α w ¯ , ξ * w ¯ , ξ * ) ] + I 2 * μ 2 ϕ ¨ , ξ I 0 * μ 4 u ¨ , ξ ξ * ( w , ξ * + w ¯ , ξ * ) I 0 * μ 4 u ¨ , ξ * ( w , ξ ξ * + w ¯ , ξ ξ * ) + 1 α μ 2 I 0 * w ¨ * = 0

2.4. Differential Quadrature Solution Procedure

Application of the DQ method is suggested by many researchers for solving eigenvalue problems [31,32]. Indeed, this approach is based on polynomial approximation of a function at any location in the domain. Then partial derivatives of that function with respect to a domain coordinate are approximated by weighted summation of its values at all selected discrete points in the domain. Thus for the one-dimensional nano-beam considered here, displacement and rotation components and their nth derivatives may be approximated by [33]
{ u * ,   w * ,   ϕ } = q = 1 R L q ( ξ ) { u q * ,   w q * ,   ϕ q }
{ u , ξ n * ,   w , ξ n *   ,   ϕ , ξ n } | ξ = ξ p = q = 1 R C p q n { u q * ,   w q * ,   ϕ q } ,   n = 1 ,   2 ,   ,   R 1 ,   p = 1 ,   2 ,   ,   R
where R is the number of sampling grid points in the domain 0 ξ 1 , L q ( ξ ) is the Lagrange interpolation polynomials and , ξ n represents nth derivative with respect to ξ . The Chebyshev-Gauss-Lobatto grid points are used in this study and are given by [33]:
ξ i = 1 2 ( 1 c o s π i R 1 ) ,   i = 0 ,   1 , , R 1
Recursive formula for weighting coefficients, C p q n , can be found in [33]. Introduction of DQ approximation of spatial derivatives of displacement and rotation components, Equations (20) and (21), into Equation (16) yields
D 0 * [ q = 1 R C p q ( 2 ) u q * + α ( q = 1 R C p q ( 1 ) w q * q = 1 R C p q ( 2 ) w q * + w ¯ , ξ ξ * | ξ = ξ p q = 1 R C p q ( 1 ) w q * + w ¯ , ξ * | ξ = ξ p q = 1 R C p q ( 2 ) w q * ) ] = I 0 * ( u ¨ * | ξ = ξ p μ 2 q = 1 R C p q ( 2 ) u ¨ q * ) G * ( q = 1 R C p q ( 2 ) w q * + 1 α q = 1 R C p q ( 1 ) Φ q * ) + Λ 1 * μ 2 Λ 2 * = I 0 * ( w ¨ * | ξ = ξ p μ 2 q = 1 R C p q ( 2 ) w ¨ q * ) I 0 * μ 2 α { q = 1 R C p q ( 2 ) u ¨ q * ( w ¯ , ξ * | ξ = ξ p + q = 1 R C p q ( 1 ) w q * ) + q = 1 R C p q ( 1 ) u ¨ q * ( w ¯ , ξ ξ * | ξ = ξ p + q = 1 R C p q ( 2 ) w q * ) μ 2 [ q = 1 R C p q ( 4 ) u ¨ q * ( w ¯ , ξ * | ξ = ξ p + q = 1 R C p q ( 1 ) w q * ) + q = 1 R C p q ( 3 ) u ¨ q * ( w ¯ , ξ ξ * | ξ = ξ p + q = 1 R C p q ( 2 ) w q * ) + q = 1 R C p q ( 2 ) u ¨ q * ( w ¯ , ξ ξ ξ * | ξ = ξ p + q = 1 R C p q ( 3 ) w q * ) + q = 1 R C p q ( 1 ) u ¨ q * ( w ¯ , ξ ξ ξ ξ * | ξ = ξ p + q = 1 R C p q ( 4 ) w q * ) ] } D 2 * q = 1 R C p q ( 2 ) ϕ q * 1 α G * ( q = 1 R C p q ( 1 ) w q * + 1 α ϕ | ξ = ξ p ) = I 2 * ( Φ ¨ | ξ = ξ p μ 2 q = 1 R C p q ( 2 ) ϕ ¨ q * )
where
Λ 1 * = D 0 * α   { w ¯ , ξ * | ξ = ξ p q = 1 R C p q ( 2 ) u q * + w ¯ , ξ ξ * | ξ = ξ p q = 1 R C p q ( 1 ) u q * + q = 1 R C p q ( 2 ) w q * [ q = 1 R C p q ( 1 ) u q * + 3 2 α ( q = 1 R C p q ( 1 ) w q * ) 2 + 3 α w ¯ , ξ ξ * | ξ = ξ p q = 1 R C p q ( 1 ) w q * + α ( w ¯ , ξ * | ξ = ξ p ) 2 ] + w ¯ , ξ * | ξ = ξ p [ q = 1 R C p q ( 2 ) u q * + α ( 3 2 w ¯ , ξ ξ * | ξ = ξ p q = 1 R C p q ( 1 ) w q * + 2 w ¯ , ξ * | ξ = ξ p w ¯ , ξ ξ * | ξ = ξ p ) ] } Λ 2 * = D 0 * α   { w ¯ , ξ * | ξ = ξ p q = 1 R C p q ( 4 ) u q * + 3 w ¯ , ξ ξ * | ξ = ξ p q = 1 R C p q ( 3 ) u q * + 3 w ¯ , ξ ξ ξ * | ξ = ξ p q = 1 R C p q ( 2 ) u q * + w ¯ , ξ ξ ξ ξ * | ξ = ξ p q = 1 R C p q ( 1 ) u q * + q = 1 R C p q ( 4 ) w q * [ q = 1 R C p q ( 1 ) u q * + 3 2 α ( q = 1 R C p q ( 1 ) w q * ) 2 + 3 α w ¯ , ξ * | ξ = ξ p q = 1 R C p q ( 1 ) w q * + ( w ¯ , ξ * | ξ = ξ p ) 2 ] + q = 1 R C p q ( 3 ) w q * [ 3 q = 1 R C p q ( 2 ) u q * + 9 α q = 1 R C p q ( 1 ) w q * q = 1 R C p q ( 2 ) w q * + 9 α w ¯ , ξ * | ξ = ξ p q = 1 R C p q ( 2 ) w q * + 9 α w ¯ , ξ ξ * | ξ = ξ p q = 1 R C p q ( 1 ) w q * + 6 α w ¯ , ξ * | ξ = ξ p w ¯ , ξ ξ * | ξ = ξ p ] + q = 1 R C p q ( 2 ) w q * [ 3 q = 1 R C p q ( 3 ) u q * + 3 α ( q = 1 R C p q ( 2 ) w q * ) 2 + 9 α w ¯ , ξ ξ * | ξ = ξ p q = 1 R C p q ( 2 ) w q * + 9 α w ¯ , ξ ξ ξ * | ξ = ξ p q = 1 R C p q ( 1 ) w q * + 6 α ( w ¯ , ξ ξ * | ξ = ξ p ) 2 + 6 α w ¯ , ξ * | ξ = ξ p w ¯ , ξ ξ ξ * | ξ = ξ p ] + q = 1 R C p q ( 1 ) w q * [ q = 1 R C p q ( 4 ) u q * + 3 2 α w ¯ , ξ ξ ξ ξ * | ξ = ξ p q = 1 R C p q ( 1 ) w q * + 2 α w ¯ , ξ * | ξ = ξ p w ¯ , ξ ξ ξ ξ * | ξ = ξ p + 6 α w ¯ , ξ ξ * | ξ = ξ p w ¯ , ξ ξ ξ * | ξ = ξ p ] }
Boundary conditions at the end points of a hinged-hinged beam can be implemented using DQ discretization as
u 𝓀 * = w 𝓀 * = 0 D 2 * q = 1 R C 𝓀 q ( 1 ) Φ q * D 0 * μ 2 [ w ¯ , ξ * q = 1 R C 𝓀 q ( 2 ) u q * + w ¯ , ξ ξ * q = 1 R C 𝓀 q ( 1 ) u q * + q = 1 R C 𝓀 q ( 1 ) w q * ( q = 1 R C 𝓀 q ( 2 ) u q * + 3 2 α w ¯ , ξ ξ * q = 1 R C 𝓀 q ( 1 ) w q * + 2 α w ¯ , ξ ξ * w ¯ , ξ * ) + q = 1 R C 𝓀 q ( 2 ) w q * ( q = 1 R C 𝓀 q ( 1 ) u q * + 3 2 α   ( q = 1 R C 𝓀 q ( 1 ) w q * ) 2 + 3 α w ¯ , ξ * q = 1 R C 𝓀 q ( 1 ) w q * + α w ¯ , ξ * w ¯ , ξ * ) ] + I 2 * μ 2 q = 1 R C 𝓀 q ( 1 ) Φ ¨ q * I 0 * μ 4 q = 1 R C 𝓀 q ( 2 ) u ¨ q * ( q = 1 R C 𝓀 q ( 1 ) w q * + w ¯ , ξ * ) I 0 * μ 4 q = 1 R C 𝓀 q ( 1 ) u ¨ q * ( q = 1 R C 𝓀 q ( 2 ) w q * + w ¯ , ξ ξ * ) + 1 α μ 2 I 0 * w ¨ * = 0
where
𝓀 = 1   at   ξ = 0   and   𝓀 = R   at   ξ = 1
Considering the displacement vector of each grid point as d i = { u i * ,   w i * , Φ i * } T for i = 1 ,   , R , it is possible to write Equation (23) in the following form
[ K ] 3 R × 3 R { X } 3 R × 1 + [ M ] 3 R × 3 R { X ¨ } 3 R × 1 = { 0 } 3 R × 1
where [ K ] is the stiffness matrix that can be decomposed into linear and nonlinear parts ( [ K ] = [ K ] L + [ K ] N L ) . [ M ] is the mass matrix and { X } ( = { d i } T   ,   i = 1 , , R ) is the generalized vector of displacement components. A harmonic motion { X } = { X ˜ } e i ω * t * (where ω * = ω L I 0 / E A is the dimensionless frequency and ω is the vibration frequency of the nano-beam) is considered. By substitution of this expression for { X } in Equation (26) and using the decomposition of stiffness matrix, an eigenvalue equation will be obtained as
( [ K ] L + [ K ] N L ) { X ˜ } ω * 2 [ M ] { X ˜ } = { 0 }
The iterative solution procedure to Equation (27) starts with ignoring the nonlinear part of the stiffness matrix and solving the corresponding linear eigenvalue problem. Thus it is possible to find eigenvalues and eigenvectors of each linear vibration mode. The vector of considered mode shape is appropriately scaled up in a way that the maximum component of lateral displacement is equal to a given dimensionless nonlinear vibration amplitude w m a x * . After that, using the modified eigenvector, the nonlinear part of stiffness matrix is calculated and a new eigenvalue problem is formed to solve for new eigenvalues and eigenvectors. The corresponding eigenvector of considered mode shape is again scaled up and the procedure continues until the relative difference between the eigenvalue of the current iteration and previous one is lower than a predefined value [34].

3. Numerical Results and Discussions

3.1. Results Verification

In this section, the values of fundamental frequency obtained by the present method are compared to those available for various cases in the literature. The number of sampling grid points for all the numerical results presented herein and in the next section is set to R = 11, which is the minimum value obtained by performing convergence studies for nonlinear vibration of a nonlocal Timoshenko nano-beam.
Table 1 presents the comparison between the results of nonlinear to linear fundamental frequency ratio ( ω n l / ω l ) of a Timoshenko beam obtained by the present solution approach and those obtained by the finite element method [35,36]. Various values of dimensionless nonlinear vibration amplitude, w max / I 2 / A , and different boundary conditions are considered. Beam geometrical and mechanical properties are given in Marur [36] and imperfection amplitude η, and nonlocal parameter μ, are set to zero. Close correlation between present iterative DQ approach and finite element method is achieved.
To maintain the confidence in results obtained based on nonlocal constitutive equations, the dimensionless linear frequencies ( ω l 2 ρ A / E I 2 ) of a straight (η = 0) SWCNT are listed in Table 2 for different values of the nonlocal parameter, and various boundary conditions. For possibility of comparison, the SWCNT properties are chosen as follows: radius r = 0.339 nm, tube thickness h = 0.066 nm, Young’s modulus E = 5.5 TPa, Poisson’s ratio υ = 0.19, and shear correction factor Ks = 0.563 [37]. As expected, increasing the nonlocal effect causes a decrease in dimensionless frequency. A noteworthy agreement is observed between the present results and those reported by Yang et al. [22]. Good agreement is observed between the results of current study and those by Wang et al. [37]. However, because they neglected the nonlocal terms in shearing force relation, the discrepancy between the results of two methods increases for higher values of nonlocal parameter.
For the last verification case of the proposed method, Table 3 shows the linear and nonlinear dimensionless frequencies for a straight (η = 0) single-walled nanotube with r = 0.313 nm, L = 5 nm, E = 1.1556 TPa, υ = 0.19, h = 0.34 nm and Ks = 0.563. Nonlocal parameter, μ , and nonlinear vibration amplitude, w m a x * , are set to 0.15 and 0.4, respectively. For both linear and nonlinear frequencies, good agreement exists between the results of current study and those reported in [22].

3.2. Geometric Imperfection Function

The following functional form of initial geometric imperfection is considered for the purpose of numerical computations:
w ¯ * = η   cos π ( 2 ξ γ β ) 2 ( γ β ) ξ 0 ξ ξ 1 w ¯ * = 0 elsewhere
where η is the dimensionless amplitude of geometrical imperfection. Figure 2 shows the schematic of the imperfection function considered in this study.
To be able to study the effects of location and extension of the geometric imperfection on the vibration characteristics of SWCNT, four cases as described in Figure 3 are considered. The “G” type represents a symmetrical full-length imperfection while the local imperfections with variable endpoint locations are denoted by “L1” to “L3” types. It is seen that the imperfection is exactly located at the center of SWCNT for L3 type while for the L1 type the imperfection is the closest case to the end.

3.3. Nonlinear Vibration of Curved SWCNTs

After verification of the formulation and solution procedure, in this section, the results for a SWCNT with Young’s modulus of E = 1.1556   TPa and Poisson’s ratio ν = 0.19 are presented. The geometrical properties are taken as: mean radius, r = 0.5   nm , thickness of the tube h = 0.3   nm and the total length, L = 5   nm unless otherwise specified. The magnitude of shear correction factor for a tubular cross-section is given by K s = 0.563 [22]. Global and local initial geometric imperfections are considered as discussed in Section 3.2 for clamp-clamp (C–C), clamp-hinged (C–H) and hinged-hinged (H–H) boundary conditions.
Figure 4 demonstrates the nonlinear to linear frequency ratio ( ω n l / ω l ) versus maximum nonlinear vibration amplitude w m a x * of a curved SWCNT with “G” type geometric imperfection with and without considering nonlocal parameter μ. Note that for all graphs, although ω n l is the nonlinear frequency of the perfect or imperfect SWCNT, ω l is the linear frequency of the perfect one. It can be observed that for perfect cases (η = 0), symmetry with respect to w m a x * = 0 exists. Thus the beam shows the same vibrational characteristics for positive and negative values of vibration amplitude. This is not the case for an imperfect beam as Figure 4 illustrates. It is mainly because the initial geometric imperfection is not generally symmetric with respect to x-axis which annihilates the symmetry of the problem. The nonlinear to linear frequency ratio increases while the vibration amplitude w m a x * is increased which is expected in large amplitude vibration of structures as a consequence of hard-spring behavior. Increasing nonlocal parameter μ causes a decrease in value of both linear and nonlinear frequencies. As Figure 4 indicates, the nonlinear to linear frequency ratio increases when the nonlocal parameter is increased, which means that this decrease in frequency value is more drastic for the linear frequency. It can also be observed that the rate of increase of frequency ratio is greater when a nonlocal parameter is present for both straight and curved SWCNTs.
The influence of the location of geometric imperfection on the nonlinear frequency ratio of curved SWCNTs with μ = 0.1 is presented in Figure 5 for H–H and C–C boundary conditions. “L1,” “L2” and “L3” are local imperfections with maximum amplitude of η = 0.1. For all local imperfect cases, the hard-spring behavior is observed. It is seen that the local imperfection has the largest effect on the nonlinear frequency ratio when it is located in the center of the single-walled carbon nanotube (L3) and its effect is decreased when the imperfection recedes from the center to the boundary (L1).
Figure 6 focuses on the effect of nonlocal parameter on the nonlinear frequency ratio of SWCNTs. The local case “L3” is selected as the geometric imperfection with maximum amplitude of η = 0.1. Unlike straight SWCNTs (perfect case), curved nanotubes (imperfect cases) do not have an overall increasing trend for nonlinear frequency ratio while the nonlocal parameter μ raises, and there is a certain value of amplitude vibration w m a x * in which the effect of nonlocal parameter gets reversed. Furthermore, it is observed that SWCNTs with C–C boundary condition have the highest nonlinear frequency and the H–H SWCNTs have the lowest, as expected.
Figure 7 depicts imperfection sensitivity of the linear fundamental frequency of a H–H curved SWCNT with “G” type geometric imperfection for different values of nonlocal parameter μ. The imperfection sensitivity indicator is defined as [24]:
S w = ω i m p e r f e c t ω p e r f e c t ω p e r f e c t × 100 %
where ω p e r f e c t and ω i m p e r f e c t are dimensionless fundamental linear frequencies of straight and curved SWCNTs, respectively. Indeed, imperfection sensitivity indicator, S w , presents the percent of changes in the linear frequency because of existence of imperfections. For a constant value of nonlocal parameter μ, it is obvious that the imperfection sensitivity, S w , increases as the imperfection amplitude η is increased. On the other hand, for a constant value of η, increasing the nonlocal parameter causes a decrease in imperfection sensitivity. It means that curved SWCNTs with higher value of nonlocal parameter μ are less sensitive to “G” type geometric imperfections. Besides, imperfection sensitivity indicator can demonstrate the quantity of error due to neglecting imperfections in modeling. For instance, for all values of nonlocal parameter and η < 0.2, imperfection sensitivity and therefore the error of neglecting imperfection of curved SWCNTs is less than 2%.
Investigation on the imperfection sensitivity of the linear fundamental frequency of curved SWCNTs with local (L2) geometry is illustrated in Figure 8. Like the “G” type imperfection, the imperfection sensitivity S w increases when the maximum imperfection amplitude η is raised. For small values of imperfection amplitude η, all values of the nonlocal parameter almost have the same behavior, but after a certain value of η, a sudden large increase in imperfection sensitivity value occurs for higher values of nonlocal parameter μ. It depicts that for the case of local imperfection with large values of nonlocal parameter, initial geometric imperfections may cause considerable change in the fundamental frequency which cannot be neglected.

4. Conclusions

In the present study, based on the nonlocal theory of Eringen, nonlinear vibration of curved SWCNTs with various boundary conditions has been investigated through the DQ solution approach. Curved SWCNTs are considered as a Timoshenko beam model with global and local initial geometric imperfections. Numerical results for a curved SWCNT indicate that the value and the location of geometric imperfection have considerable effect on the linear and nonlinear vibration characteristics of curved nano-beams. Higher imperfection amplitudes clearly increase the vibration frequency of the nano-beam. Furthermore, centrally located geometric imperfections have the largest influence on increasing the frequency ratios of curved SWCNTs. On the other hand, a small-scale parameter may also increase or decrease the nonlinear frequency ratio and sensitivity indicator depending on the amplitude, type and location of the imperfection of the nano-beam.

Acknowledgments

The financial support of Iran Nanotechnology Initiative Council is gratefully acknowledged. N.M.P. is supported by the European Research Council (ERC StG Ideas 2011 BIHSNAM n. 279985, ERC PoC 2015 SILKENE nr. 693670), by the European Commission under the Graphene Flagship (WP14 Polymer Composites, no. 696656).

Author Contributions

Seyed K. Jalali suggested the idea and formulated the problem; Iman Eshraghi solved the equation of motion numerically; Seyed K. Jalali and Iman Eshraghi analyzed the results and wrote the manuscript under the supervision of Nicola Maria Pugno.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CNTsCarbon Nano Tubes
DQDifferential Quadrature
MDMolecular Dynamics
NEMSNano Electro Mechanical System
SWCNTsSingle-Walled Carbon Nano Tubes

References

  1. Bosia, F.; Lepore, E.; Alvarez, N.T.; Miller, P.; Shanov, V.; Pugno, N. Knotted synthetic polymer or carbon nanotube microfibres with enhanced toughness, up to 1400 J/g. Carbon 2016, 102, 116–125. [Google Scholar] [CrossRef]
  2. Ke, C.H.; Pugno, N.; Peng, B.; Espinosa, H.D. Experiments and modeling of carbon nanotube-based NEMS devices. J. Mech. Phys. Solids 2005, 53, 1314–1333. [Google Scholar] [CrossRef]
  3. Ke, C.H.; Espinosa, H.D.; Pugno, N. Numerical analysis of nanotubes based NEMS devices—Part II: Role of finite kinematics, stretching and charge concentrations. J. Appl. Mech. 2005, 72, 519–526. [Google Scholar] [CrossRef]
  4. Hierold, C.H.; Jungen, A.; Stampfer, C.H.; Helbling, T.H. Nano electromechanical sensors based on carbon nanotubes. Sens. Actuators A Phys. 2007, 136, 51–61. [Google Scholar] [CrossRef]
  5. Liew, K.M.; Wang, Q. Analysis of wave propagation in carbon nanotubes via elastic shell theories. Int. J. Eng. Sci. 2007, 45, 227–241. [Google Scholar] [CrossRef]
  6. Strozzi, M.; Smirnov, V.V.; Manevitch, L.I.; Milani, M.; Pellicano, F. Nonlinear vibrations and energy exchange of single-walled carbon nanotubes. Circumferential flexural modes. J. Sound Vib. 2016, 381, 156–178. [Google Scholar] [CrossRef]
  7. Wang, C.Y.; Ru, C.Q.; Mioduchowski, A. Applicability and limitations of simplified elastic shell equations for carbon nanotubes. J. Appl. Mech. 2004, 71, 622–631. [Google Scholar] [CrossRef]
  8. Smirnov, V.V.; Manevitch, L.I.; Strozzi, M.; Pellicano, F. Nonlinear optical vibrations of single-walled carbon nanotubes, Energy exchange and localization of low-frequency oscillations. Physica D 2016, 325, 113–125. [Google Scholar] [CrossRef]
  9. Silvestre, N.; Wang, C.M.; Zhang, Y.Y.; Xiang, Y. Sanders shell model for buckling of single-walled carbon nanotubes with small aspect ratio. Compos. Struct. 2011, 93, 1683–1691. [Google Scholar] [CrossRef]
  10. Amabili, M. A comparison of shell theories for large-amplitude vibration of circular cylindrical shells: Lagrangian approach. J. Sound Vib. 2003, 264, 1091–1125. [Google Scholar] [CrossRef]
  11. Gibson, R.F.; Ayorinde, O.E.; Wen, Y.F. Vibration of carbon nanotubes and there composites: A review. Compos. Sci. Technol. 2007, 67, 1–28. [Google Scholar] [CrossRef]
  12. Eringen, A.C. Nonlocal Continuum Field Theories; Springer: New York, NY, USA, 2002. [Google Scholar]
  13. Barretta, R.; Feo, L.; Luciano, R.; Marotti de Sciarra, F. An Eringen-like model for Timoshenko nanobeams. Compos. Struct. 2016, 139, 104–110. [Google Scholar] [CrossRef]
  14. Barretta, R.; Feo, L.; Luciano, R.; Marotti de Sciarra, F.; Penna, R. Functionally graded Timoshenko nanobeams: A novel nonlocal gradient formulation. Compos. B 2016, 100, 208–219. [Google Scholar] [CrossRef]
  15. Čanadija, M.; Barretta, R.; Marotti de Sciarra, F. On functionally graded Timoshenko nonisothermal nanobeams. Compos. Struct. 2016, 135, 286–296. [Google Scholar] [CrossRef]
  16. Wang, Q.; Varadan, V.K. Wave characteristics of carbon nanotubes. Int. J. Solids Struct. 2006, 43, 254–265. [Google Scholar] [CrossRef]
  17. Lu, P.; Lee, H.P.; Lu, C.; Zhang, P.Q. Application of nonlocal beam models for carbon nanotubes. Int. J. Solids Struct. 2007, 44, 5289–5300. [Google Scholar] [CrossRef]
  18. Murmu, T.; Pradhan, S.C. Thermo-mechanical vibration of a single-walled carbon nanotube embedded in an elastic medium based on nonlocal elasticity theory. Comput. Mater. Sci. 2009, 46, 854–859. [Google Scholar] [CrossRef]
  19. Ansari, R.; Sahmani, S. Small scale effect on vibrational response of single-walled carbon nanotubes with different boundary conditions based on nonlocal beam models. Commun. Nonlinear Sci. 2012, 17, 1965–1979. [Google Scholar] [CrossRef]
  20. Reddy, J.N. Nonlocal nonlinear formulations for bending of classical and shear deformation theories of beams and plates. Int. J. Eng. Sci. 2010, 48, 1507–1518. [Google Scholar] [CrossRef]
  21. Reddy, J.N. Nonlocal theories for bending, buckling, and vibration of beams. Int. J. Eng. Sci. 2007, 45, 288–307. [Google Scholar] [CrossRef]
  22. Yang, J.; Ke, L.L.; Kitipornchai, S. Nonlinear free vibration of single-walled carbon nanotubes using nonlocal Timoshenko beam theory. Physica E 2010, 42, 1727–1735. [Google Scholar] [CrossRef]
  23. Fang, B.; Zhen, Y.; Zhang, C.; Tang, Y. Nonlinear vibration analysis of double-walled carbon nanotubes based on nonlocal elasticity theory. Appl. Math. Model. 2013, 37, 1096–1107. [Google Scholar] [CrossRef]
  24. Kitipornchai, S.; Yang, J.; Liew, K.M. Semi-analytical solution for nonlinear vibration of laminated FGM plates with geometric imperfections. Int. J. Solids Struct. 2004, 41, 2235–2257. [Google Scholar] [CrossRef]
  25. Yamaki, N.; Otomo, K.; Chiba, M. Nonlinear vibrations of a clamped rectangular plate with initial deflection and initial edge displacement-Part 2: Experiment. Thin. Walled Struct. 1983, 1, 101–109. [Google Scholar] [CrossRef]
  26. Calabri, L.; Pugno, N.; Ding, W.; Ruoff, R. Resonance of curved nanowires. J. Phys. Condens. Matter 2006, 18, S2175–S2183. [Google Scholar] [CrossRef]
  27. Farshidianfar, A.; Soltani, P. Nonlinear flow-induced vibration of a SWCNT with a geometrical imperfection. Comput. Mater. Sci. 2012, 53, 105–116. [Google Scholar] [CrossRef]
  28. Ouakad, H.M.; Younis, M.I. Natural frequencies and mode shapes of initially curved carbon nanotube resonators under electric excitation. J. Sound Vib. 2011, 330, 3182–3195. [Google Scholar] [CrossRef]
  29. Amabili, M. Nonlinear vibrations of rectangular plates with different boundary conditions: Theory and experiments. Comput. Struct. 2004, 82, 2587–2605. [Google Scholar] [CrossRef]
  30. Eringen, A.C. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. J. Appl. Phys. 1983, 54, 4703–4710. [Google Scholar] [CrossRef]
  31. Jalali, S.K.; Rastgoo, A.; Eshraghi, I. Large amplitude vibration of imperfect shear deformable nano-plates using non-local theory. J. Solid Mech. 2011, 3, 64–73. [Google Scholar]
  32. Pradhan, S.C.; Murmu, T. Small Scale Effect on the Buckling of Single-Layered Graphene Sheets under Bi-axial Compression via Nonlocal Continuum Mechanics. Comput. Mater. Sci. 2009, 47, 268–274. [Google Scholar] [CrossRef]
  33. Shu, C. Differential Quadrature and Its Application in Engineering; Springer: London, UK, 2000. [Google Scholar]
  34. Malekzadeh, P.; Vosoughi, A.R. DQM large amplitude vibration of composite beams on nonlinear elastic foundations with restrained edges. Commun. Nonlinear Sci. 2009, 14, 906–915. [Google Scholar] [CrossRef]
  35. Bhasyam, G.R.; Prathap, G. Galerkin finite element method for non-linear beam vibrations. J. Sound Vib. 1980, 72, 191–203. [Google Scholar] [CrossRef]
  36. Marur, S.R.; Prathap, G. Non-linear beam vibration problems and simplifications in finite element models. Comput. Mech. 2005, 35, 352–360. [Google Scholar] [CrossRef]
  37. Wang, C.M.; Zhang, Y.Y.; He, X.Q. Vibration of nonlocal Timoshenko beams. Nanotechnology 2007, 18, 105401. [Google Scholar] [CrossRef]
Figure 1. An initially curved beam with annular cross-section representing a SWCNT.
Figure 1. An initially curved beam with annular cross-section representing a SWCNT.
Materials 09 00786 g001
Figure 2. Schematic of initial imperfection of the nano-beam.
Figure 2. Schematic of initial imperfection of the nano-beam.
Materials 09 00786 g002
Figure 3. Various imperfection types considered in numerical results.
Figure 3. Various imperfection types considered in numerical results.
Materials 09 00786 g003
Figure 4. Variation of nonlinear vibration frequency ratio for different combinations of imperfection amplitude η and nonlocal parameter μ for a “G” type imperfection for (a) clamped-clamped (b) clamped-hinged and (c) hinged-hinged endpoint conditions.
Figure 4. Variation of nonlinear vibration frequency ratio for different combinations of imperfection amplitude η and nonlocal parameter μ for a “G” type imperfection for (a) clamped-clamped (b) clamped-hinged and (c) hinged-hinged endpoint conditions.
Materials 09 00786 g004
Figure 5. Variation of nonlinear vibration frequency ratio for (a) clamped-clamped and (b) hinged-hinged endpoint conditions of a nonlocal nano-beam for various “L” types imperfection.
Figure 5. Variation of nonlinear vibration frequency ratio for (a) clamped-clamped and (b) hinged-hinged endpoint conditions of a nonlocal nano-beam for various “L” types imperfection.
Materials 09 00786 g005
Figure 6. Effects of nonlocal parameter μ on the nonlinear vibration frequency ratio of nano-beam with “L3” type of initial imperfection.
Figure 6. Effects of nonlocal parameter μ on the nonlinear vibration frequency ratio of nano-beam with “L3” type of initial imperfection.
Materials 09 00786 g006
Figure 7. Effect of nonlocal parameter μ on sensitivity indicator of a nano-beam with “G” type initial imperfection with hinged-hinged endpoint conditions.
Figure 7. Effect of nonlocal parameter μ on sensitivity indicator of a nano-beam with “G” type initial imperfection with hinged-hinged endpoint conditions.
Materials 09 00786 g007
Figure 8. Effect of nonlocal parameter μ on sensitivity indicator of a nano-beam with “L2” type initial imperfection for (a) clamped-clamped (b) clamped-hinged and (c) hinged-hinged endpoint conditions.
Figure 8. Effect of nonlocal parameter μ on sensitivity indicator of a nano-beam with “L2” type initial imperfection for (a) clamped-clamped (b) clamped-hinged and (c) hinged-hinged endpoint conditions.
Materials 09 00786 g008aMaterials 09 00786 g008b
Table 1. Comparison of nonlinear frequency ratio of a Timoshenko straight beam for different vibration amplitudes.
Table 1. Comparison of nonlinear frequency ratio of a Timoshenko straight beam for different vibration amplitudes.
w max I 2 / A H–HC–CC–H
Ref. [35]Ref. [36]PresentRef. [35]Ref. [36]PresentRef. [35]Ref. [36]Present
1.01.1181.1181.11811.02951.02831.02951.06411.05821.0593
2.01.41411.41351.41431.11271.11051.11281.23181.2151.2182
3.01.80261.80271.80291.23771.23361.23781.46031.43681.4416
4.02.23592.23612.23631.39201.38561.39211.72101.68221.7026
5.02.69232.69252.69281.56591.55741.56601.99951.91801.9862
Table 2. Comparison of dimensionless linear frequency ( ω l 2 ρ A / E I 2 ) of a straight nano-beam with different nonlocal parameter µ and various endpoint conditions.
Table 2. Comparison of dimensionless linear frequency ( ω l 2 ρ A / E I 2 ) of a straight nano-beam with different nonlocal parameter µ and various endpoint conditions.
μH–HC–CC–H
Ref. [37]Ref. [22]PresentRef. [37]Ref. [22]PresentRef. [37]Ref. [22]Present
0.03.09293.09294.44914.44913.78453.7844
0.13.02433.02103.02104.34714.32694.32693.69393.68493.6849
0.32.65382.63852.63853.78953.70323.70323.21153.17243.1724
0.52.28672.26652.26653.24203.13723.13712.74712.69822.6980
0.72.01061.98992.83832.73272.40592.3569
Table 3. Comparison of linear and nonlinear dimensionless frequency of a straight nonlocal nano-beam with μ = 0.15 and w m a x * = 0.4 .
Table 3. Comparison of linear and nonlinear dimensionless frequency of a straight nonlocal nano-beam with μ = 0.15 and w m a x * = 0.4 .
FrequencyH–HC–CC–H
Ref. [22]PresentRef. [22]PresentRef. [22]Present
ω l * 0.42330.42330.80550.80550.60520.6052
ω n l * 0.44050.44350.81880.82190.61970.6236

Share and Cite

MDPI and ACS Style

Eshraghi, I.; Jalali, S.K.; Pugno, N.M. Imperfection Sensitivity of Nonlinear Vibration of Curved Single-Walled Carbon Nanotubes Based on Nonlocal Timoshenko Beam Theory. Materials 2016, 9, 786. https://doi.org/10.3390/ma9090786

AMA Style

Eshraghi I, Jalali SK, Pugno NM. Imperfection Sensitivity of Nonlinear Vibration of Curved Single-Walled Carbon Nanotubes Based on Nonlocal Timoshenko Beam Theory. Materials. 2016; 9(9):786. https://doi.org/10.3390/ma9090786

Chicago/Turabian Style

Eshraghi, Iman, Seyed K. Jalali, and Nicola Maria Pugno. 2016. "Imperfection Sensitivity of Nonlinear Vibration of Curved Single-Walled Carbon Nanotubes Based on Nonlocal Timoshenko Beam Theory" Materials 9, no. 9: 786. https://doi.org/10.3390/ma9090786

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop