Next Article in Journal
Analytical Modeling for Mechanical Straightening Process of Case-Hardened Circular Shaft
Previous Article in Journal
Explainable Artificial Intelligence (XAI) and Supervised Machine Learning-based Algorithms for Prediction of Surface Roughness of Additively Manufactured Polylactic Acid (PLA) Specimens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Free Vibration of Single-Walled Carbon Nanotubes Using Nonlocal Truncated Timoshenko-Ehrenfest Beam Theory

by
Maria Anna De Rosa
1,
Maria Lippiello
2,*,
Antonella Onorato
1 and
Isaac Elishakoff
3
1
School of Engineering, University of Basilicata, 85100 Potenza, Italy
2
Department of Structures for Engineering and Architecture, University of Naples “Federico II”, Via Forno Vecchio n°36, 80134 Naples, Italy
3
Department of Ocean and Mechanical Engineering, Florida Atlantic University, Boca Raton, FL 33431-0991, USA
*
Author to whom correspondence should be addressed.
Appl. Mech. 2023, 4(2), 699-714; https://doi.org/10.3390/applmech4020035
Submission received: 27 February 2023 / Revised: 25 April 2023 / Accepted: 26 April 2023 / Published: 12 May 2023

Abstract

:
Carbon nanotubes with their outstanding mechanical, physical and electrical properties have stimulated a significant amount of scientific and technological research due to their uniqueness compared to conventional materials. As a result, an extensive study on their mechanical properties has been conducted, and the static and dynamic behavior of single- walled and multi-walled carbon nanotubes has been examined using Euler-Bernoulli and Timoshenko beam models. The main objective of this paper is to study the free vibration behaviour of single-walled carbon nanotubes (SWCNT) using the nonlocal truncated Timoshenko beam theory. According to the Hamilton principle, the equation of motion of Timoshenko single-walled carbon nanotubes is calculated taking into account the truncated theory; and the general corresponding boundary conditions are derived. Finally, some numerical examples are performed to evaluate the effects of the nonlocal coefficient and the length of the nanotube. The obtained results are validated by comparing them with those found in the literature, and they show the accuracy and efficiency of the developed model. Particularly, the results demonstrate that the present formulation is highly efficient and capable of satisfactorily describing the behavior of nanobeams.

1. Introduction

Carbon nanotubes, with their outstanding mechanical, physical and electrical properties have stimulated extensive research activities in the field of science and technology due to their uniqueness compared to conventional materials. As a result, extensive studies have been conducted to investigate their mechanical properties [1,2,3] and the static and dynamic behavior of single- and multi-walled carbon nanotubes have been carried out.
In the literature, the main existing approaches to study the behavior of nanostructures can be divided into two classes: One at the atomistic level, and the other at the continuous level. Although the molecular dynamics approach can be considered more suitable for the analysis of nano-sized structures, the complexity of the computational process and the time-consuming calculations [4] have attracted little attention compared to continuum models. Among continuum approaches, beam models have proven to be cost-effective. Using the classical Euler-Bernoulli and Timoshenko beam models, extensive studies have been conducted to investigate the mechanical properties of CNTs and elastic beam models have been used to evaluate static and dynamic problems of carbon nanotubes [5,6,7,8,9,10].
An important aspect of nanotube modeling is the scale-effect. Although classical continuum approaches are efficient in the mechanical analysis of CNTs, their applicability in identifying small-scale effects on the mechanical behavior of carbon nanotubes is questionable. Numerous studies have underlined the importance of the scale-effect and have shown that non-local elastic continuum models are more suitable for predicting the structural behavior of nanotubes.
The origins of non-local elasticity theory go back to pioneering works by Eringen [11,12] and Reddy [13]. Using non-local constitutive differential equations, Reddy gives a thorough development of the traditional beam and shear deformation theories in [13] and derives the solutions for, bending, buckling and natural frequency problems for simply-supported beams. Further applications of the nonlocal elasticity theory have been used in the study of the buckling problem and vibration problems in CNTs by applying Euler-Bernoulli beam, plate and shell theories and Timoshenko beam theory [14,15,16,17,18,19,20,21,22,23,24,25,26,27]. The wave dispersion response of elastic nano-sized beams was analytically addressed and the closed-form solution of the phase velocity was determined in [28]. Still based on the nonlocal elasticity beam theory, the effects of the spatial variation of the nonlocal parameter on free frequencies of nanoplates were discussed in [29]. Using nonlocal first-order shear deformation theory with variable nonlocal parameters, the free vibration of the functionally graded doubly curved nanoshells was studied in [30].
The present paper deals with the nonlocal vibrational analysis of the single-walled carbon nanotube (SWCNT). Using the nonlocal elasticity theory, a novel theory for the free vibration analysis of Timoshenko nanotubes is proposed. The main objective of this paper is to investigate the free vibration behaviour of single-walled carbon nanotubes using the nonlocal truncated Timoshenko beam theory. According to Hamilton principle, the equation of motion of Timoshenko single-walled carbon nanotube are calculated and the general corresponding boundary conditions are derived.
The mathematical modeling of CNTs and their solutions play an important role in the field of nanotechnology. As is well-know, it is not always possible to find analytical solutions for all sets of boundary conditions and different geometries. This circumstance has motivated the present research. In addition, since in the literature there are no free vibration solutions for nanotubes with the simplified approach, the authors have chosen to briefly present the key findings and challenges and direct light to possible future research. This paper does not intersect with recent relevant reviews, which reflects its significance to readers.
Recently, a variational model leading to a simplified theory of the dynamic analysis of the Timoshenko beam theory has been proposed, called the truncated Timoshenko theory [31]. Starting from the truncated Timoshenko beam theory, the authors in the present paper have developed a similar theory for Timoshenko nanotubes using geometrical and variational methods. Comparing the classical Timoshenko theory with the truncated theory, a significant simplification of the differential equation can be observed. The novelty of the proposed approach is that it shows a perfect analogy between variational and direct methods for the dynamic analysis of beams. The aim of the proposed formulations is to find the truncated Timoshenko equations and the corresponding boundary conditions and to establish their mathematical similarity with the two different approaches. It is shown that the differential equations and the corresponding boundary conditions, used to solve the dynamic problem of the local and nonlocal truncated Timoshenko equations via the variational formulation, have the same form as those obtained via the direct method. In addition, it is shown that the equation is both simpler and more consistent than the appropriate classical Timoshenko equations extended to include nonlocal stress effects.
Finally, the impacts of the nonlocal coefficient and the length of the nanotube are then assessed using a few numerical examples. By comparing the obtained results with those found in the literature, the findings are validated. Moreover, the outcomes demonstrate how effective the simplified approach is at accurately describing the behavior of nanotubes.

2. Theoretical Formulation: Equations of Motion for Nonlocal Timoshenko Beams Model

Eringen Nonlocal Theory Assumptions and Stress Resultants in Nonlocal Theory

According to Eringen [11,12], the nonlocal elasticity theory assumes that the Cauchy stress state, at a given reference point x, does not depend only on the strain field at this point, but is considered as a function of the strain field at all points of the body. This assumption leads us to the assertion that, in nonlocal elasticity theory, the stress at a point is determined both by the stress at that point and by its spatial derivatives. Eringen attributed this fact to the atomic theory of lattice dynamics and experimental observations on phonon dispersion. For an elastic and isotropic homogeneous body, the nonlocal constitutive behavior is expressed by the following relations:
σ = V α ( | x x | , τ ) t ( x ) d x
t ( x ) = C ( x ) ϵ ( x )
1 τ 2 l e 2 2 σ = t , τ = e 0 l i l e
where σ is the nonlocal stress tensor at point x, t(x) is the conventional stress tensor at point x, α (|x’ − x|, τ ) is the nonlocal coefficient introducing into the constitutive equation the nonlocal effect at the reference point x produced by local strain at the source x’. The Euclidean distance is |x’ − x|, τ is defined as the scale coefficient that incorporates the small-scale factor, where e 0 is a material constant determined experimentally. In addition, l i and l e are the internal and external characteristic lengths, respectively. Finally, t at a point x is related to the strain ϵ at the point by the fourth-order elasticity tensor C and ⊗ denotes the “double-dot product”.
Using Equation (3) we can express the resulting stresses in terms of strains. In contrast to the local theory, the nonlocal constitutive relations lead to differential relationship between resulting stresses and strains. In the following, the nonlocal constitutive relations for isotropic and homogeneous beam are described. In particular, for beams the nonlocal constitutive relations in Equation (3) take the following special form:
σ zz μ 2 2 σ zz z 2 = E e zz , μ = e 0 l i
where σ zz and e zz are the normal stress and normal strain, respectively, E is elasticity modulus, μ is the small scale parameter that incorporates the small scale effect and z is the coordinate measured from the mid-plane of the nanotube.
σ zy μ 2 2 σ zy z 2 = 2 G e zy
where σ zy and e zy are the transverse shear stress and the transverse shear strain, respectively, G is shear modulus and y is the longitudinal coordinate measured from the left-end of the nanotube.
Multiplying Equation (4) by y and integrating the result over the area A yields:
A y σ zz dA μ 2 A 2 y σ zz z 2 dA = E A y 2 ϕ z dA
The integrals with the first member indicate the bending moment, whereas the integral with the second member defines the second moment of area I of the straight section. As a result:
M μ 2 2 M z 2 = EI ϕ z
where M is the bending moment and ϕ is the rotation due to bending. In addition, by integrating Equation (5) over the area, one gets:
Q μ 2 2 Q z 2 = GA κ v z + ϕ
where Q is the shear force, A is the area of the straight section, v the transverse displacement and κ the shear factor.

3. Timoshenko’s Non-Local Truncated Theory: Dynamic Analysis

3.1. Equation of Motion for a Truncated Timoshenko Beam: Euler Method

This section presents the equilibrium of the nanotube element using the Timoshenko-Ehrenfest beam model. For example, if we examine the nanotube element in Figure 1, we should impose the equilibrium of the applied loads at the abscissa z: equilibrium at the vertical translation and at the rotation around the center of the right-hand section of all the forces identified on the elementary ashlar, beginning from its equilibrium.
The contributions for rotational m and translational inertia f I are made separately and considering that
f I = ρ A 2 v t 2
m = ρ I 2 ϕ b t 2
where ρ denotes mass density, t denotes the time variable, the rotational inertial term m is independent of the total rotation ϕ and exclusively depends on the flexural rotation ϕ b = v z :
ϕ = ϕ b + ψ = v z + ψ
The equation for the vertical translation equilibrium is:
Q z = ρ A 2 v t 2
whose first derivative, when substituted in Equation (8), yields:
Q = GA κ v z + ϕ + μ 2 ρ A 3 v t 2 z
hence, Equation (12) can be completed and expressed as follows:
ρ A 2 v t 2 μ 2 ρ A 4 v t 2 z 2 GA κ 2 v z 2 + ϕ z = 0
On the other hand, the rotation’s equilibrium with respect to the center of gravity elementary ashlar’s right face is written:
M z = Q ρ I 3 v t 2 z
whose first derivative, when replaced in Equation (7), gives:
M = EI ϕ z + μ 2 ρ I 4 v t 2 z 2 + ρ A 2 v t 2
The rotation’s equilibrium Equation (15) can be written:
EI 2 ϕ z 2 μ 2 ρ I 5 v t 2 z 3 + ρ I 3 v z t 2 GA κ v z + ϕ = 0
Equations (14) and (17) constitute the set of differential equations of motion for a Timoshenko SWCNT that were obtained from the Timoshenko beam’s truncated theory.

3.2. Equations of Motion for a Truncated Timoshenko Beam: Variational Method

In this section, according to the Hamilton principle, the motion equations for a beam system are derived.
For the local theory, the strain energy Π can be expressed as follows:
Π = 1 2 0 L M ϕ z d z + 0 L Q v z + ϕ d z
It is evident from Equations (13) and (16), which are used to determine the shear stress Q and the bending moment M in the case of non-local elasticity, that the terms, defined in this context as additional to the local theory, are constant quantities that do not change with respect to the force; this means that these additional terms do not contribute to the strain energy but instead become the potential energy of the loads P (see Ref. [23]).
The strain energy Π assumes the following form:
Π = 1 2 0 L EI ϕ z 2 d z + 1 2 0 L GA κ v z + ϕ 2 d z
As can be seen, only the local terms are present in the formulation of the strain energy.
The nonlocal potential energy P d is equal to the work done with changed sign for the corresponding displacement of the inertial forces.
P d = 0 L μ 2 ρ I 4 v t 2 z 2 + ρ A 2 v t 2 ϕ z + μ 2 ρ A 3 v t 2 z v z + ϕ dz
As can be observed, only the moment and shear nonlocal contributions accurately reflect the nonlocal potential energy expression in contrast to the local theory.
Moreover, for the corresponding displacement ϕ , the rotational component P t is expressed in terms of the work of the inertia forces m, with the sign changed:
P t = 0 L m ϕ d z = 0 L ρ I 2 ϕ b t 2 ϕ d z = 0 L ρ I 2 t 2 v z ϕ d z
Finally, the kinetic energy T is given by:
T = 1 2 0 L ρ A v t 2 d z

Equations of Motion for a Timoshenko SWCNT: Truncated Theory

According to the Hamilton principle, the equations of motion for a Timoshenko single-walled carbon nanotube are computed using the truncated theory.
Summing the strain energy Equation (19) and the potential energies Equations (20) and (21) minus the kinetic energy Equation (22), we have:
Π + P d + P t T = 1 2 0 L EI ϕ z ϕ z dz + 0 L μ 2 ρ A 2 v t 2 ρ I 4 v t 2 z 2 ϕ z dz + 0 L μ 2 ρ A 3 v t 2 z ϕ + v z dz + 1 2 0 L GA κ ϕ + v z ϕ + v z dz 1 2 0 L ρ A v t 2 dz 0 L ρ I 3 v z t 2 ϕ dz
Applying the Hamilton principle, Equation (23) becomes:
t 1 t 2 δ Π + P d + P t T dt = t 1 t 2 ( 0 L EI ϕ z δ ϕ z dz + 0 L μ 2 ρ I 4 v t 2 z 2 + μ 2 ρ A 2 v t 2 δ ϕ z dz + 0 L ( μ 2 ρ A 3 v t 2 z δ ϕ + μ 2 ρ A 3 v t 2 z δ v z + GA κ ϕ + v z δ ϕ + GA κ ϕ + v z δ v z ) dz 0 L ρ A v t δ v t dz 0 L ρ I 3 v z t 2 δ ϕ dz ) dt = 0
performing integration by parts and collecting all the terms in the previous equations allows us to write:
ρ A 2 v t 2 μ 2 ρ A 4 v t 2 z 2 GA κ 2 v z 2 + ϕ z = 0
EI 2 ϕ z 2 μ 2 ρ I 5 v t 2 z 3 + ρ I 3 v z t 2 GA κ v z + ϕ = 0
The Equations (25) and (26) represent the set of differential equations of motion for a Timoshenko SWCNT and, as can be seen, have the same form as the Equations (14) and (17) obtained by Euler method.
Finally, the corresponding boundary conditions are given by:
μ 2 ρ A 2 v t 2 μ 2 ρ I 4 v t 2 z 2 + EI ϕ z δ ϕ 0 L = 0
μ 2 ρ A 3 v t 2 z + GA κ v z + ϕ δ v 0 L = 0

3.3. The Solution of Differential Equations System

In order to find a solution to the system of differential equations describing the dynamic behavior of the nanotube, we seek periodic solutions of the form:
v ( z , t ) = v ( z ) e i ω t ϕ ( z , t ) = ϕ ( z ) e i ω t
where ω is the frequency of natural vibration. On substituting Equation (29) into Equations (25) and (26) we have:
ρ A ω 2 v μ 2 ρ A ω 2 2 v z 2 + GA κ 2 v z 2 + ϕ z = 0
EI 2 ϕ z 2 + μ 2 ρ I ω 2 3 v z 3 ρ I ω 2 v z GA κ v z + ϕ = 0
After obtaining ϕ z from Equation (30) and properly substituting it in Equation (31), we arrive at:
EI + μ 2 ρ I ω 2 + EI μ 2 ρ A ω 2 GA κ 4 v z 4 + EI ρ A ω 2 GA κ μ 2 ρ A ω 2 ρ I ω 2 2 v z 2 + ρ A ω 2 v = 0
Equation (32) is the differential equation of motion for a Timoshenko nanotube derived from the fourth order truncated theory in v.
On introducing the following non-dimensional coefficients:
ζ = z L ; Ω 2 = ρ A L 4 ω 2 EI ; η 2 = μ 2 L 2 ; β 2 = EI GA κ L 2 ; α 2 = I A L 2
with ζ [ 0 , 1 ] the dimensionless counterpart of z [ 0 , L ] , Ω 2 the frequency parameter, η 2 scaling effect parameter, β 2 shear deformation parameter and α 2 slenderness ratio, the governing Equation (32) may be rewritten as:
1 Ω 2 α 2 η 2 Ω 2 β 2 η 2 4 v ζ 4 + Ω 2 β 2 + Ω 2 η 2 + Ω 2 α 2 2 v ζ 2 Ω 2 v = 0
The general solution for Equation (34) is given by:
v ( ζ ) = A 1 cos ( a ζ ) + A 2 sin ( a ζ ) + A 3 cosh ( b ζ ) + A 4 sinh ( b ζ )
where:
a = 1 2 p q + q 2 + 4 pr ; b = 1 2 p q + q 2 + 4 pr
p = 1 Ω 2 α 2 η 2 Ω 2 β 2 η 2 q = Ω 2 β 2 + Ω 2 η 2 + Ω 2 α 2 ; r = Ω 2
with A 1 , A 2 , A 3 and A 4 being the integration constants, which are determined using the boundary conditions, and p, q and r are the coefficients of the characteristic polynomial associated with differential Equation (34).
The system of Equations (25) and (26) in dimensionless form is given by:
1 β 2 Ω 2 η 2 2 v ζ 2 + Ω 2 v + 1 β 2 ( ϕ L ) ζ = 0 Ω 2 α 2 η 2 3 v ζ 3 + 1 β 2 Ω 2 α 2 v ζ + 2 ( ϕ L ) ζ 2 1 β 2 ( ϕ L ) = 0
It is possible to derive ϕ ¯ = ϕ L
ϕ ¯ ( ζ ) = Ω 2 α 2 η 2 β 2 + Ω 2 β 4 η 2 β 2 3 v ζ 3 1 + Ω 2 α 2 β 2 + Ω 2 β 4 v ζ

3.4. Boundary Conditions

The quantities defined by Equations (27) and (28) lead to the following boundary conditions in dimensionless form:
1 β 2 + Ω 2 η 2 v ζ 1 β 2 ϕ ¯ v ( ζ ) 0 1 = 0 ϕ ¯ ζ Ω 2 α 2 η 2 2 v ζ 2 + Ω 2 η 2 v ϕ ¯ ( ζ ) 0 1 = 0
By substituting the Equations (35) and (39), that were appropriately calculated, in the boundary conditions (40), we obtain a system of four equations in the four unknowns A i . In order for the system to admit a solution other than the trivial one, the determinant of the coefficient matrix must be set to zero.
The infinite solutions of the transcendental equation obtained by solving the determinant provide the infinite free frequencies of vibration.

3.5. Comparison of the Two Methods

This section compares the boundary conditions and differential equations of motion for the two theories. In particular, we compare the boundary conditions and differential equations of motion for the classical theory (see Equation (18) to Reference [23]) with the truncated theory (see Equation (34) of the current paper).
As demonstrated in Reference [23], the differential equation’s characteristic polynomial for a classical SWCNT has the following form:
p C λ 4 + q C λ 2 + r C = 0
with:
p C = 1 Ω 2 α 2 η 2 Ω 2 β 2 η 2 + Ω 4 α 2 β 2 η 4 ; q C = Ω 2 β 2 + Ω 2 η 2 + Ω 2 α 2 2 Ω 4 α 2 β 2 η 2 ; r C = Ω 2 Ω 4 α 2 β 2 η 2
where p C , q C and r C are the coefficients of the differential equation’s characteristic polynomial. As shown in the present paper, the characteristic polynomial associated with differential Equation (34) based on truncated Timoshenko’s theory is represented by:
p λ 4 + q λ 2 r = 0
with:
p = 1 Ω 2 α 2 η 2 Ω 2 β 2 η 2 ; q = Ω 2 β 2 + Ω 2 η 2 + Ω 2 α 2 ; r = Ω 2
The comparison of the two Equations (42) and (44) shows that the terms p C and p differ by the amount ( α 2 β 2 η 4 Ω 4 ), q C and q differ by the amount ( α 2 2 β 2 η 4 Ω 4 ) and the terms r C and r differ by the amount ( β 2 η 4 Ω 4 ). As you can see, they are all terms multiplicative in 4 v t 4 .
By examining the free frequencies of vibration of a simply-supported beam at both ends, Timoshenko came to the conclusion that the final multiplicative term may be omitted because it would not affect the results, as was previously noted in Reference [32] in which Elishakoff obtained an equation both more consistent and simpler than the Bresse-Timoshenko equation given by the following expression:
EI 4 v x 3 4 + ρ A 2 v t 2 ρ I 1 + E G χ 4 v x 3 2 t 2 + ρ 2 I G χ 4 v t 4 = 0
It can easily be seen from the comparison of the boundary conditions that those relating to displacements and shear remain constant while those relating to rotations and momentum change with respect to a single term. The following equations show how the classical theory and the nonlocal theory for a quantity differ from one another:
1 Ω 2 η 2 α 2 ϕ ¯ ζ + Ω 2 η 2 v ϕ ¯ ( ζ ) 0 1 = 0
Equation (46) is related to the classical theory, while the following equation relates to the nonlocal theory:
ϕ ¯ ζ Ω 2 α 2 η 2 2 v ζ 2 + Ω 2 η 2 v ϕ ¯ ( ζ ) 0 1 = 0
in which it can be noted that Ω 2 η 2 α 2 ϕ ¯ ζ becomes Ω 2 α 2 η 2 2 v ζ 2 .

4. Numerical Examples

In this section, the suggested analytical method is validated in order to assess the impacts of the nonlocal parameter and the length of the nanotube on the free frequencies of vibration. Some numerical examples have been carried out, and the obtained results have been compared with those of papers that have already been published in the literature. All the numerical computations have been performed through in-house software developed in the Mathematica language [33] and the same geometrical features have been used.
For this purpose, Table 1 shows physical and geometrical properties of nanotubes used in Reference [34], which will be used throughout this section. In addition, the following shear factor is introduced:
k = 6 ( 1 + ν ) 1 + c r 2 2 ( 7 + 6 ν ) 1 + c r 2 2 + ( 20 + 12 ν ) c r 2
where c r = (d − 2 h)/d. The external diameter d depends directly on the height h in the ratio d/h, and the span L depends on d in the ratio L/d.

Numerical Comparison between the Two Theories: Conventional and Truncated Theories

A first numerical comparison has been made between the results of the present paper and those obtained by Wu and Lai in Ref. [34]. For different values of the non-dimensional small-scale coefficient η 0 , 2 , 2 , the aspect ratio of L/d [5, 10, 25, 50, 100], with d/h = 3, and for various boundary conditions, the first three nondimensional frequencies Ω i have been calculated. In Table 2 are quoted the corresponding values and a comparison between the present results with those obtained in [34] has been done. As can be easily observed, for L/d = [5, 10], the results of the classical theory and the truncated theory diverge, whereas they coincide for values higher than L/d = 10. Moreover, they demonstrate that with increasing the ratio L/d, the first three nondimensional frequencies increase and it can be seen that if the nonlocal effect η increases the three first nondimensional frequency value decreases. In particular, in order to assess the impact of the small scale parameter on the frequency parameter of the SWCNT, the frequency parameter ratio as R( η ) = Ω ( η )/ Ω ( η = 0) is defined. For various boundary conditions the following remarks may be applicable:
-
For L/d = 5 and d/h = 3, the lowest frequency parameter ratio of the simply-supported SWCNT with η = 2 , 2 are around 0.75 and 0.63, and they are 0.50, 0.36, and 0.36, 0.26 for the second and third modes, respectively. These findings demonstrate that the impact of the small-scale parameter on the frequency parameter of the SWCNT is greater for higher vibration modes than for lower modes.
-
For L/d = 10 and d/h = 3, the nonlocal effect parameter η = 2 and for four boundary conditions (simply-supported (SS), clamped-clamped (CC), clamped-supported (CS) and clamped-free (CF) the frequency parameter ratios for the first vibration mode of the SWCNT are 0.75, 0.72, 0.73 and 1.03, respectively. According to the results, the SWCNT with SS, CC, and CS boundary conditions exhibit a small scale parameter effect on the frequency parameter that is more important than the SWCNT with CF boundary conditions.
In comparison to Table 3 of Ref. [34], the ratio L/d = 5 has been added and the dimensionless free vibration frequencies are calculated using the algorithm developed in Ref. [23] and referred to as additional constraint conditions other than the clamped-free. These values are crucial because, for larger L/d values, we move toward the Euler-Bernoulli theory, for which the classical Timoshenko and truncated Timoshenko theories will coincide; but for lower ratios, we approach Timoshenko’s theory.
Table 2. Three non-dimensional frequency values for η = 0 , 2 , 2, d/h = 3 and L/d [5, 10, 25, 50, 100].
Table 2. Three non-dimensional frequency values for η = 0 , 2 , 2, d/h = 3 and L/d [5, 10, 25, 50, 100].
η [nm]d/hL/d Ω i SSCCCSCF
[34]Present[34]Present[34]Present[34]Present
0351st9.33499.325418.318818.265313.606213.58073.42133.4208
2nd32.756732.445542.397341.937837.766437.370318.687518.6268
3rd63.173961.573071.100769.505867.292565.674144.933544.3155
101st9.72549.724721.101921.095114.893814.89123.49153.4914
2nd37.339437.301554.424854.322845.732245.667021.025821.0189
3rd79.130578.811898.962498.437189.094788.676355.626155.5123
251st9.84609.846022.153522.153315.330315.33023.51203.5120
2nd39.105639.104460.303660.299249.198149.195721.862821.8625
3rd86.979086.9659116.3020116.2720101.2204101.200060.572460.5679
501st9.86379.863722.317722.317715.396115.39613.51503.5150
2nd39.384039.383961.321161.320849.769449.769221.991221.9911
3rd88.350188.3501119.6965119.6940103.4615103.460061.409361.4090
1001st9.86819.868122.359322.359315.412715.41273.51583.5158
2nd39.454739.454761.584361.584349.915749.915722.023622.0236
3rd88.706788.7066120.5979120.5980104.0492104.049061.624861.6248
2 351st7.03397.031413.252113.19389.96029.93813.52323.5225
2nd16.305916.150920.421020.081918.564218.320512.642712.5757
3rd22.577322.004225.423824.561224.052723.341921.259720.8030
101st8.91638.915619.000618.993913.484613.48233.51953.5195
2nd28.154428.125839.840739.752433.981333.929018.675618.6677
3rd48.090147.896558.914358.521453.567153.286640.463640.3540
251st9.69999.699921.746921.746615.066115.06613.51673.5167
2nd36.927936.926756.453556.449346.249946.247721.425721.4255
3rd77.084277.0725101.9404101.908089.192389.172556.691256.6867
501st9.82659.826522.211922.211915.328415.32843.51623.5162
2nd38.799538.799460.264660.264348.968948.968821.879421.8794
3rd85.480085.4792115.3995115.396099.916899.914860.352560.3522
1001st9.85889.858822.333022.333015.395615.39563.51613.5161
2nd38.305939.305961.313661.313649.711249.711221.995521.9955
3rd87.958987.9588119.4686119.4680103.1218103.121061.354661.3546
2351st5.88305.877010.864210.80268.21568.19523.65073.6495
2nd12.318312.201315.321315.043614.014213.81829.76499.7140
3rd16.499316.081218.661417.976117.606317.066216.464516.0494
101st8.28048.279817.409017.402212.406312.40433.54883.5489
2nd23.531823.507932.862832.780928.245628.199316.870916.8630
3rd37.659437.507845.999245.661241.880241.649733.468333.3679
251st9.56019.560121.361121.360914.814914.81483.52133.5213
2nd35.077635.076553.245553.241243.774243.772021.009621.0093
3rd69.944269.933691.843491.810580.634980.615753.467853.4633
501st9.78979.789722.109622.109615.261615.26163.51733.5173
2nd38.240338.240259.259759.259348.205748.205621.769021.7690
3rd82.871982.8710111.5318111.527096.712896.710359.347959.3476
1001st9.84959.849522.306722.306715.378715.37873.51633.5164
2nd39.158739.158761.046361.046349.509249.509221.967521.9675
3rd87.229787.2297118.3702118.3690102.2186102.218061.087961.0879
Figure 2, Figure 3, Figure 4 and Figure 5 show the nondimensional frequency values versus the ratio h/d, for three different values of the nonlocal parameter η = 0, 1, 2 and various boundary conditions. According to Figure 2, Figure 3 and Figure 4 in relation to the cases of the simply-supported, clamped-clamped, clamped-supported SWCNT, for a fixed aspect ratio value L/d = 10 and for 2 < d / h < 10 , the non-dimensional frequency values decrease as the nonlocal effect value increases. In addition, when d/h = 10 and higher values of the non-dimensional frequency are found for the clamped-clamped situation, the frequency curves for η = 0 and η = 1 approach and tend to coincide. Ultimately, the frequency increases for η = 1 and η = 2 and as the ratio d/h varies, whereas it decreases for η = 0. The case of the clamped-free single-walled carbon nanotube is explored in Figure 5. For a fixed aspect ratio value L/d = 10 and for 2 < d / h < 10 , the non-dimensional frequency values versus the ratio d/h, for three alternative values of the nonlocal parameter η = 0, 1, 2, are provided. As can be observed, the nondimensional frequency increases if the nonlocal parameter increases and for the ratio d/h= 10 the frequency curves approach and tend to coincide.
Table 3 compares the results from the truncated Timoshenko theory and those produced using the conventional method, showing the percentage of errors in each case. This inaccuracy grows with the nonlocal effect η and is greater for the clamped-clamped boundary condition. The simply supported nanotube does not change while η changes. From the first to the third frequency, this error increases.

5. Conclusions

In the present paper, the nonlocal truncated Timoshenko single-walled nanotube, with various boundary conditions, has been studied. According to the Eringen’s and nonlocal Timoshenko beam theory, a new theory for the free vibration analysis of the Timoshenko nanotube has been presented and the equations of motion and the general corresponding boundary conditions have been derived.
As is well known, it is not always possible to find analytical solutions for all sets of boundary conditions and varied geometries. This circumstance has motivated the present research. The novelty of the proposed approach has been devoted to finding solutions to the eigenvalue problem of single-walled carbon nanotube by employing the truncated theory in those cases for which, to the authors’ knowledge, no analytical solutions have been provided in the literature.
Starting from the truncated Timoshenko beam theory, the authors have developed the same theory for Timoshenko nanotubes using the geometrical and variational methods. Comparing the classical Timoshenko theory with the truncated theory, we have shown a considerable simplification of the differential equation. Some numerical examples have shown the effectiveness of the proposed approach, even through a comparison against results in the literature. In particular, the results have demonstrated that the present formulation is very efficient and able to describe the behavior of nanobeams in a satisfactory way.
Finally, it has been explained how the length of the nanotube and nonlocal parameters affect the three initial frequencies of SWCNT. The main points are:
-
the first three nondimensional frequencies decrease with increasing of nonlocal effect;
-
the first three nondimensional frequencies increase with increasing the ratio L/d;
-
the effect of the small scale parameter on the frequency parameter of the SWCNT with SS, CC and CS boundary conditions is more significant than that of the SWCNT with CF boundary conditions.
The present approach can be applied to analyze the stability and dynamics behaviour of single-walled and multi-walled carbon nanotubes (MWCNTs) with different boundary conditions. The first results obtained, and reported in the degree thesis [35], are significant and in perfect agreement with the results of the literature. In the upcoming publication, which will be released, the effects of surface tension on the natural frequency of nanotubes will also be investigated.
Moreover, the exact nonlocal truncated Timoshenko beam solutions presented herein should be useful to engineers who are designing nano-beams and nanotubes. In addition, the exact solutions serve as reference results for confirming numerical vibration solutions derived from other mathematical models and methods.

Author Contributions

Conceptualization, M.A.D.R. and M.L.; methodology, M.A.D.R. and M.L.; software, A.O.; validation, M.A.D.R., M.L. and A.O.; formal analysis, M.A.D.R. and A.O.; investigation, A.O.; resources, M.A.D.R. and M.L.; data curation, M.A.D.R. and A.O.; writing—original draft preparation, M.L.; writing—review and editing, M.A.D.R. and M.L.; visualization, I.E.; supervision, M.A.D.R. and I.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Krishnan, A.; Dujardin, E.; Ebbesen, T.W.; Yianilos, P.N.; Treacy, M.M.J. Young’s modulus of single-walled nanotubes. Phys. Rev. B 1998, 58, 14013–14019. [Google Scholar] [CrossRef]
  2. Demczyk, B.G.; Wang, Y.M.; Cumings, J.; Hetman, M.; Han, W.; Zettl, A.; Ritchie, R.O. Direct mechanical measurement of the tensile strength and elastic modulus of multiwalled carbon nanotubes. Mat. Sci. and Eng. A 2002, 334, 173–178. [Google Scholar] [CrossRef]
  3. Ruoff, R.S.; Qian, D.; Liu, W.K. Mechanical properties of carbon nanotubes: Theoretical predictions and experimental measurements. C. R. Phys. 2003, 4, 993–1008. [Google Scholar] [CrossRef]
  4. Ansari, R.; Rouhi, H.; Rad, A.N. Vibrational analysis of carbon nanocones under different boundary conditions: An analytical approach. Mech. Res. Commun. 2014, 56, 130–135. [Google Scholar] [CrossRef]
  5. Elishakoff, I.; Versaci, C.; Muscolino, G. Clamped-free double-walled car-bon nanotube-based mass sensor. Acta Mech. 2011, 219, 29–43. [Google Scholar] [CrossRef]
  6. Elishakoff, I.; Pentaras, D. Fundamental natural frequencies of double-walled nanotubes. Int. J. Sound Vib. 2009, 322, 652–664. [Google Scholar] [CrossRef]
  7. Ru, C.Q. Column buckling of multiwalled carbon nanotubes with interlayer radial displacements. Phys. Rev. B 2000, 62, 16962–16967. [Google Scholar] [CrossRef]
  8. Wang, Q.; Hu, T.; Chen, Q.J. Bending instability characteristics of double walled nanotubes. Phys. Rev. B 2005, 71, 045403. [Google Scholar] [CrossRef]
  9. Zhang, Y.Y.; Wang, C.M.; Tan, V.B.C. Buckling of multiwalled carbon nanotubes using Timoshenko beam theory. J. Eng. Mech. 2006, 132, 952–958. [Google Scholar] [CrossRef]
  10. Rafiee, R.; Moghadam, R.M. On the modeling of carbon nanotubes: A critical review. Compos. Part B Eng. 2014, 56, 435–449. [Google Scholar] [CrossRef]
  11. Eringen, A.C. Nonlocal polar elastic continua. Int. J. Eng. Sci. 1972, 10, 1–16. [Google Scholar] [CrossRef]
  12. Eringen, A.C. On differential equations of non local elasticity and solutions of screw dislocation and surface-waves. J. Appl. Phys. 1983, 54, 4703–4710. [Google Scholar] [CrossRef]
  13. Reddy, J.N. Nonlocal theories for bending, buckling and vibration of beams. Int. J. Eng. Sci. 2007, 45, 288–307. [Google Scholar] [CrossRef]
  14. Wang, L.; Hu, H. Flexural wave propagation in single-walled carbon nanotubes. Phys. Rev. B 2005, 71, 195412. [Google Scholar] [CrossRef]
  15. Lu, P.; Lee, H.P.; Lu, C.; Zhang, P.Q. Dynamic properties of exural beams using a nonlocal elasticity model. Appl. Phys. 2006, 99, 073510. [Google Scholar] [CrossRef]
  16. De Rosa, M.A.; Lippiello, M. Free vibration analysis of DWCNTs using CDM and Rayleigh-Schimdt based on nonlocal Euler-Bernoulli beam theory. Sci. World J. 2014, 2014, 194529. [Google Scholar] [CrossRef]
  17. Reddy, J.N.; Pang, S.D. Nonlocal continuum theories of beams for the analysis of carbon nanotubes. J. Appl. Phys. 1998, 103, 023511. [Google Scholar] [CrossRef]
  18. De Rosa, M.A.; Lippiello, M. Hamilton principle for SWCN and a modified approach for nonlocal frequency analysis of nanoscale biosensor. Int. J. Recent Sci. Res. 2015, 6, 2355–2365. [Google Scholar]
  19. Kucuk, I.; Sadek, I.S.; Adali, S. Variational principles for multiwalled carbon nanotubes undergoing vibrations on nonlocal Timoshenko beam theory. J. Nanomat. 2010, 2010, 461252. [Google Scholar] [CrossRef]
  20. De Rosa, M.A.; Lippiello, M.; Martin, H.D. Free vibrations of a cantilevered SWCNT with distributed mass in the presence of non local effect. Sci. World J. 2015, 2015, 825342. [Google Scholar] [CrossRef]
  21. De Rosa, M.A.; Lippiello, M. Nonlocal frequency analysis of embedded singlewalled carbon nanotube using the Differential Quadrature Method. Compos. Part B Eng. 2016, 84, 41–51. [Google Scholar] [CrossRef]
  22. De Rosa, M.A.; Lippiello, M.; Martin, H.D.; Piovan, M.T. Nonlocal frequency analysis of nanosensors with different boundary conditions and attached distributed biomolecules: An approximate method. Acta Mech. 2016, 227, 2323–2342. [Google Scholar] [CrossRef]
  23. De Rosa, M.A.; Lippiello, M. Nonlocal Timoshenko frequency analysis of single-walled carbon nanotube with attached mass: An alternative Hamiltonian approach. Compos. Part B Eng. 2017, 111, 409–418. [Google Scholar] [CrossRef]
  24. De Rosa, M.A.; Lippiello, M.; Tomasiello, S. Differential quadrature solutions for the nonconservative instability of a class of single-walled carbon nanotubes. Eng. Comput. 2018, 35, 251–267. [Google Scholar] [CrossRef]
  25. Auciello, N.M.; De Rosa, M.A.; Lippiello, M.; Tomasiello, S. Non-conservative instability of cantilevered nanotube via cell discretization method. Springer Proc. Math. Stat. 2018, 248, 13–24. [Google Scholar]
  26. De Rosa, M.A.; Lippiello, M.; Auciello, N.M.; Martin, H.D.; Piovan, M.T. Variational method for non-conservative instability of a cantilever SWCNT in the presence of variable mass or crack. Arch. Appl. Mech. 2021, 91, 301–316. [Google Scholar] [CrossRef]
  27. De Rosa, M.A.; Lippiello, M.; Babilio, E.; Ceraldi, C. Nonlocal Vibration Analysis of a Nonuniform Carbon Nanotube with Elastic Constraints and an Attached Mass. Materials 2021, 14, 3445. [Google Scholar] [CrossRef]
  28. Faghidian, S.A.; Tounsi, A. Dynamic characteristics of mixture unified gradient elastic nanobeams. Facta Univ. Ser. Mech. Eng. 2022, 20, 539–552. [Google Scholar] [CrossRef]
  29. Vinh, P.V.; Tounsi, A. The role of spatial variation of the nonlocal parameter on the free vibration of functionally graded sandwich nanoplates. Eng. Comp. 2022, 38, 4301–4319. [Google Scholar] [CrossRef]
  30. Vinh, P.V.; Tounsi, A. Free vibration analysis of functionally graded doubly curved nanoshells using nonlocal first-order shear deformation theory with variable nonlocal parameters. Thin-Walled Struct. 2022, 174, 109084. [Google Scholar] [CrossRef]
  31. De Rosa, M.A.; Lippiello, M.; Elishakoff, I. Variational Derivation of Truncated Timoshenko-Ehrenfest Beam Theory. J. Appl. Comput. Mech. 2022, 8, 996–1004. [Google Scholar]
  32. Elishakoff, I. An Equation Both More Consistent and Simpler Than Bresse-Timoshenko Equation. In Advances in Mathematical Modeling and Experimental Methods for Materials and Structures; Gilat, R., Sills-Banks, L., Eds.; Springer: Berlin/Heidelberg, Germany, 2009; pp. 249–254. [Google Scholar]
  33. Wolfram, S. The Mathematica 8; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  34. Wu, C.P.; Lai, W.W. Free vibration of an embedded single-walled carbon nanotube with various boundary conditions using the RMVT-based nonlocal Timoshenko beam theory and DQ method. Phys. E 2015, 68, 8–21. [Google Scholar] [CrossRef]
  35. Onorato, A. A New Theory for Dynamic Analysis of Timoshenko SWCNT. Degree Thesis, University of Basilicata “School of Engineering”, Potenza, Italy, 2023. [Google Scholar]
Figure 1. Translational and rotational equilibrium of the elementary ashlar.
Figure 1. Translational and rotational equilibrium of the elementary ashlar.
Applmech 04 00035 g001
Figure 2. The nondimensional frequency values versus the ratio d/h, for three different values of the nonlocal parameter η = 0, 1, 2: simply-supported single-walled nanotube case.
Figure 2. The nondimensional frequency values versus the ratio d/h, for three different values of the nonlocal parameter η = 0, 1, 2: simply-supported single-walled nanotube case.
Applmech 04 00035 g002
Figure 3. The nondimensional frequency values versus the ratio d/h, for three different values of the nonlocal parameter η = 0, 1, 2: clamped-clamped single-walled nanotube case.
Figure 3. The nondimensional frequency values versus the ratio d/h, for three different values of the nonlocal parameter η = 0, 1, 2: clamped-clamped single-walled nanotube case.
Applmech 04 00035 g003
Figure 4. The nondimensional frequency values versus the ratio d/h, for three different values of the nonlocal parameter η = 0, 1, 2: clamped-supported single-walled nanotube case.
Figure 4. The nondimensional frequency values versus the ratio d/h, for three different values of the nonlocal parameter η = 0, 1, 2: clamped-supported single-walled nanotube case.
Applmech 04 00035 g004
Figure 5. The nondimensional frequency values versus the ratio d/h, for three different values of the nonlocal parameter η = 0, 1, 2: clamped-free single-walled nanotube case.
Figure 5. The nondimensional frequency values versus the ratio d/h, for three different values of the nonlocal parameter η = 0, 1, 2: clamped-free single-walled nanotube case.
Applmech 04 00035 g005
Table 1. Geometrical and material properties adopted in numerical experiments.
Table 1. Geometrical and material properties adopted in numerical experiments.
SWCNT PropertiesSymbolValueUnit
Cross-sectional area A π (d − h) hm 2
Thicknessh0.34 × 10 9 m
Moment of inertia I 1/8 π   ( h d 3 3 h 2 d 2 + 4 h 3 d 2 h 4 ) m 4
Density ρ 2300kg/m 3
Young’s modulusE1000GPa
Poisson’s ratio ν 0.19-
Shear modulusG E ( 2 + 2 ν ) GPa
Table 3. Percentage error between the traditional SWCNT and the result of the truncated theory, for L/d = 5, 10 per ( η = 0 , 2 ) , 2 .
Table 3. Percentage error between the traditional SWCNT and the result of the truncated theory, for L/d = 5, 10 per ( η = 0 , 2 ) , 2 .
η [nm]L/d Ω i SSCCCSCF
% Error *
051st0.1020.2920.1870.016
2nd0.9501.0841.0490.325
3rd2.5342.2432.4051.375
101st0.0080.0320.0170.002
2nd0.1020.1870.1430.033
3rd0.4030.5310.4700.205
2 51st0.0350.4400.2220.021
2nd0.9511.6611.3130.530
3rd2.5383.3932.9552.148
101st0.0080.0350.0170.000
2nd0.1020.2220.1540.042
3rd0.4030.6670.5240.271
251st0.1020.5670.2480.033
2nd0.9501.8131.3990.521
3rd2.5343.6723.0682.521
101st0.0070.0390.0160.003
2nd0.1020.2490.1640.047
3rd0.4030.7350.5500.300
* % error = approx - exact exact · 100 .
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

De Rosa, M.A.; Lippiello, M.; Onorato, A.; Elishakoff, I. Free Vibration of Single-Walled Carbon Nanotubes Using Nonlocal Truncated Timoshenko-Ehrenfest Beam Theory. Appl. Mech. 2023, 4, 699-714. https://doi.org/10.3390/applmech4020035

AMA Style

De Rosa MA, Lippiello M, Onorato A, Elishakoff I. Free Vibration of Single-Walled Carbon Nanotubes Using Nonlocal Truncated Timoshenko-Ehrenfest Beam Theory. Applied Mechanics. 2023; 4(2):699-714. https://doi.org/10.3390/applmech4020035

Chicago/Turabian Style

De Rosa, Maria Anna, Maria Lippiello, Antonella Onorato, and Isaac Elishakoff. 2023. "Free Vibration of Single-Walled Carbon Nanotubes Using Nonlocal Truncated Timoshenko-Ehrenfest Beam Theory" Applied Mechanics 4, no. 2: 699-714. https://doi.org/10.3390/applmech4020035

APA Style

De Rosa, M. A., Lippiello, M., Onorato, A., & Elishakoff, I. (2023). Free Vibration of Single-Walled Carbon Nanotubes Using Nonlocal Truncated Timoshenko-Ehrenfest Beam Theory. Applied Mechanics, 4(2), 699-714. https://doi.org/10.3390/applmech4020035

Article Metrics

Back to TopTop