Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
To purchase hard copies of this book, please contact the representative in India:
CBS Publishers & Distributors Pvt. Ltd.
www.cbspd.com
|
customercare@cbspd.com
Department of Mechanical Engineering, American University of Beirut, Lebanon
*Address all correspondence to: najib01@idm.net.lb
1. Introduction
The intended applications of carbon nanotubes have steadily increased since their discovery by Ijima [4]. They range from the nanoscale, as in the tip of an atomic electron microscope, to the macroscale, as in the preliminary design of the space elevator cable.
In modeling CNTs, molecular as well as quantum mechanics have been the primary tools for analysis. Also, closed form expressions were developed to study the response of CNTs in different environments [3,8,12].
On the other hand, some attempted to use structural mechanics, and built corresponding finite element models, to study the behavior of CNTs, as evident in several publications [5-8]. Some of these publications simplified the property relations between molecular mechanics and structural mechanics. They assumed the structural bending stiffness EI/a to be a constant and set it equal to the molecular bond bending stiffness C.
However, in [1], with a simple proof, we showed that the main assumption used by various authors to equate the element bending stiffness to the bond bending stiffness (C= EI/a) does not hold. In addition, in our previous publications (Kasti [1,2]), we related some of the mechanical properties used in molecular and structural mechanics.
In [1], we derived an expression for the axial deformation of zigzag CNTs that accounts for the axial and bending structural stiffnesses under simple tension. While molecular mechanics uses the bond angle between two bonds to describe bond bending deformations, structural mechanics uses the bending within one 3D frame element for this definition. Comparing the deformation equation in structural mechanics to the equivalent equation derived for molecular mechanics, leads us to a “consistent” frame bending stiffness for an infinitely long zigzag CNT. For large diameter tubes, the frame bending stiffness tends to half the bond bending stiffness. This later case is representative of a graphene sheet. For small diameters, EI/a changes with the bond bending stiffness C, the torsional angle φ and the lattice translational index n. The expression for the axial deformation was then expanded to include the torsional stiffness of 3D frame elements. This provided us with an equation for the “elastic modulus” of an infinitely long zigzag CNT. It included the axial, bending and torsional deformations. Then, in [2], we extended the work to the case of simple torsion of zigzag carbon nanotubes. The expression relating the molecular bond bending stiffness C and the structural bending stiffness EI/a was derived. It was found to be different from the case of simple tension. The structural bond bending stiffness was both load and chirality dependent. However, for the particular configuration of a graphene sheet, the relation of simple tension was recovered, namely EI/a=C/2. We concluded the work by presenting the expression for the deformation of the tube when axial, bending and torsional structural stiffnesses are accounted for. We noticed in this case of simple torsion that the axial stiffness couples with the bending and torsion stiffnesses, unlike simple tension.
In a recent paper (Chen et al., [9]), the radial elastic modulus of the original Molecular Structure Mechanics model (MSM) was compared to the one from the Molecular Dynamics (MD) simulation. In that paper, it was pointed to the fact that a modification to the original MSM model was suggested in our previous paper (Kasti, [1]).
In this chapter, we extend our previous work to armchair carbon nanotubes under simple tension. In addition, we summarize the equivalent results for zigzag CNTs under simple tension and torsion.
We start with a brief review of molecular and structural mechanics and we refer to the work of Chang and Gao [3]. Then, the relation between the structural bending stiffness EI/a and the molecular bond bending stiffness C is derived for the case of simple tension. This shows that EI/a depends on the bond bending stiffness C and the torsional angle φ. In the limit of an infinite tube radius, which represents a graphene sheet, we recover the previous relation, i.e., EI/a tends to C/2. Finally, an expression for the Young’s modulus is presented that accounts for the axial, bending and torsional deformations. We conclude the chapter with numerical simulations that validate the results.
This section deals with the molecular and structural mechanics formulations of bond energies. Also, a short review of the finite element method is presented as it applies to the modeling of carbon nanotubes subjected to mechanical loading.
2.1. Characterization of the Atomic Structure and Molecular/Structural Mechanics of Carbon Nanotubes (CNTs)
The geometry of a CNT could be described with the pair (n,m), the lattice translational indices, and the bond length a. In general, the diameter d of a CNT is defined using the expression
d=3aπn2+ m2+ nm E1
For zigzag CNTs (Fig. 1a), the value of m is zero. In this case, the diameter is simply given by the formulad=3anπ. However, for armchair CNTs (Fig. 1b), the value of m is equal to n. In this case, the diameter is given by the formulad=3anπ.
The bond energies between carbon atoms include the stretching Ua, bending Ub and torsional Ut energies. For small distortions from equilibrium, these energies could take the forms:
Ua=12K(r–ro)2Ub=12C(Θ–Θo)2Ut=12Ct(Ф–Фo)2E2
where K, C, and Ct are the bond stretching, bending and torsional stiffnesses.
r, Θ, and Ф are the stretched position, bending and torsional angles, respectively.
The subscript “o” refers to the initial equilibrium configuration.
As far as structural mechanics, the linear elastic deformation is assumed to be the combination of axial, bending and torsional deformations. Their corresponding strain energies are expressed as:
where EA, EI, and GJ are the axial, bending and torsional stiffnesses; and u, θ and φ are the axial, bending and torsional deformations, respectively.
Figure 1.
Carbon Nanotubes: (a) Zigzag, (b)Armchair.
2.2. Review of the Finite Element Method for 3D Space Frames
For linear elastic behavior of 3D space frames, the axial, bending and torsional strain energies can be expressed as in equation (3).
The bonding between two carbon atoms is modeled by placing a 3D space frame element between them, Fig. 2. When this procedure is repeated throughout the tube, a finite element mesh is obtained with the carbon atoms becoming the nodes in the mesh.
Each node is assumed to have six degrees of freedom, three translational and three rotational.
Figure 2.
Space frame element.
The stiffness matrix K relating the degrees of freedom to their corresponding forces and moments at both ends of a 3D frame element is a 12x12 matrix.
Chang and Gao derived closed form expressions for carbon nanotubes subjected to simple tensile loading using molecular mechanics.
Representing units of zigzag and armchair carbon nanotubes are shown in Fig. 3 with α and β being the internal angles. Equivalent equations to the ones of Chang and Gao will be derived in the next section for the armchair CNT.
Figure 3.
a) Two units of a zigzag CNT. (b) One unit of an armchair CNT.
4. Axial and Bending Stiffnesses of Armchair CNTs: Molecular versus Structural Mechanics
We start by expressing the results of Chang and Gao [3] in a more suitable form using the principle of minimum total potential energy. We will split the approach into bond stretching and bond bending deformations.
4.1. Bond Stretching – Molecular/Structural Mechanics
Due to bond stretching, it is easy to verify thatK=EA/a.
We start with the molecular energy expression
Π=∑112K(u−uo)2−∑2FΔE5
where F is the axial load applied to a single carbon atom, Δ is the deflection of the end of the tube, Σ1 is a summation over the number of bonds and Σ2 is the summation over the number of atoms with applied loads.
Expression (5) takes the following form in structural mechanics:
Π=∑112(EA/a)(u−uo)2−∑2FΔE6
Since both K and EA/a are conjugate to the axial deformation between carbon atoms in the energy equation, they represent the same axial stiffness. Thus,
K=EA/aE7
To determine the tube deformation, we let nu be the number of vertical units of Fig. 3 in a carbon nanotube and nv equal to (2nu–1). The deflection at the end of the tube due to axial bond deformations can be expressed as:
Δ=2nvFsin2(α/2)(EA/a)E8
4.2. Bond Bending – Molecular Mechanics
Due to bond bending, the total potential energy is written as:
∏=∑112C(α−α0)2+∑212C(β−β0)2−∑3FΔE9
Let nx be the number of units along the circumference and nv the number of vertical units. Then, for an infinite cylinder with no end effects,i.e., all units have the same deformation, Π will be equal to:
Since cosβ = -cos(π/2n)cos(α/2) and cosφ = tan(α/2)/tan(β), we get dβ/dα = cosφ/2, where φ is the torsion angle between the planes of adjacent units of an armchair nanotube (Fig. 3).
The vertical deformation of the tube can be expressed as:
Δ=H−H0=(2nv)a(sin(α/2)−sin(α0/2))E10
where Ho,H are the initial and current heights of the tube, respectively. Differentiating Eq. (10) with respect to α
dΔ/dα=12(2nv)acos(α/2)
Solving for Δα by minimizing the total potential energy, we get
Δα=α−α0=Facos(α/2)/(2C+Ccos2ϕ)E11
Substituting Δα in Eq. (10) above, we get the following expression for the vertical deflection of an infinitely long CNT (due to bond bending only)
Δ=(2nv)Fa2cos2(α2)4C+2Ccos2ϕE12
In the next Section, we will derive an equivalent expression in terms of the material properties of structural mechanics. This will allow us to deduce a relation between EI/a and C.
4.3. Bond Bending - Structural Mechanics: Armchair Carbon Nanotube under Simple Tension
Let O and A be two atoms on the CNT with coordinates (Rcos(-θ/2), Rsin(-θ/2), 0 ) and ( Rcos(-π/n+θ/2), Rsin(-π/n+θ/2), a sin(α/2) ), respectively. The angle θ is given by 2R2(1-cosθ)=a2.
Also, let u be the unit vector from O to A expressed as:
u=u1i+u2j+u3kE13
where i, j and k are the unit vectors in the Cartesian coordinate system.
For an infinite cylinder and due to symmetry, the radial and tangential rotations at O and A are zero.
In addition, due to multiple symmetries (for n=4,8,12,16,..), we assume the displacement and rotation fields at O and A take the following forms:
δro=δiδrA=δsin(−π/n)i+δcos(−π/n)j+ΔkE14
and
θo=θkθA=−θkE15
where ( δro, δrA ) and ( θo, θA ) are the displacements and rotations vectors at O and A respectively. In this work, small displacements and rotations are assumed. Without loss of generality, the tensile load at each carbon atom is assigned a value of one.
Following a similar procedure to Kasti [1], one can show that θ is equal to zero.
And, the only force in the inclined member OA in Fig. 3 is the vertical force F and the moment is (1/2)F.a.sin(α/2).
Due to multiple symmetries (for n=4,8,12,16,..), and after some simplifications, the axial deflection at the end of the whole CNT is given by:
Δ=(2nv)Fa2cos2(α2)12(EI/a)E16
and the corresponding elastic modulus Ys
Ys=8EIsin(α2)a3cos2(α2)E17
Zigzag CNTs under simple tension
The corresponding axial deformation at the end of a zigzag CNT under simple tension is given by [1]:
Δ=(nv+1)Fa2sin2α(EI/a)[24−36γ/(1+γ)]E18
withγ=2(1–cos(π/n))/(3–2(R/a)2sin2(π/n))E19
Zigzag CNTs under simple torsion
The tangential deformation at the end of a zigzag CNT unit under simple torsion is given by following expression for ΔB [2]:
5. Discussion and Validation of the Results obtained in Section 4 – Armchair CNTs
Comparing the molecular mechanics expression Eq. (12) with the structural mechanics Eq. (16), we obtain
EIa=C6[2+cos2ϕ]E21
Thus, in general, the bending stiffness to be used in the structural mechanics varies with the bond bending stiffness C and torsional angle φ.
For long CNT tubes with large diameters, cos2(φ) → 1 and EI/a → C/2, which is the same result derived by Kasti [1] for the particular case of a graphene sheet.
To validate the closed form solution Eq. (16), we compared the axial deformation of member OA (Fig. 3) and the change in radius to the results from a finite element model in ABAQUS [14]. The results are shown in Tables 1 and 2 below for C=1.42 nN.nm.rad-2. The accuracy obtained is excellent.
Lattice Translational Index, n
Molecular/Structural Mechanics
ABAQUS
4
5.9844
5.9849
8
5.9477
5.9470
12
5.9358
5.9358
16
5.9310
5.9306
Table 1.
Vertical deformation of member OA (x10-4).
Lattice Translational Index, n
Molecular/Structural Mechanics
ABAQUS
4
-1.3538
-1.3539
8
-2.6402
-2.6401
12
-3.9383
-3.9383
16
-5.2403
-5.2402
Table 2.
Change in radius (x10-3).
Zigzag CNTs under simple tension
For zigzag CNT under simple tension, the equivalent stiffness is given by [1]:
EIa=4C+8Ccos2ϕ24−36γ/(1+γ)E22
Thus, in general, the bending stiffness to be used in the structural mechanics varies with the bond bending stiffness C, torsional angle φ and γ.
For long CNT tubes with large diameters, cosφ → 1, γ→0 and EI/a → C/2, which is the same result for the particular case of a graphene sheet.
Zigzag CNTs under simple torsion
For a zigzag CNTs under simple torsion, the corresponding stiffness is given by [2]:
EIa=C6[1+2cos2(π2n)].E23
Thus, in general, the bending stiffness to be used in structural mechanics varies with the bond bending stiffness C and lateral translational index n.
For long CNT tubes with large diameters, n→ ∞, cos(π2n)→1 and EI/a → C/2.
6. Deformation of Armchair CNTs due to Axial, Bending and Torsional Structural Stiffnesses
In Sections 4 and 5, a closed form expression was developed for the deformation of infinitely long armchair CNT under simple tension. It included the axial and bending stiffnesses of 3D frame elements. In this Section, we study the effect of the torsional stiffness of 3D space frames.
6.1. Bond Bending and Torsion - Structural Mechanics: Armchair Nanotube under Simple Tension
Similar work to Kasti [1] will show that the torsional stiffness does not enter the expression for the deformation of an infinitely long armchair carbon nanotube under simple tension.
When the axial, bending and torsional deformations are combined, we obtain the following formula:
Δ=2nvFsin2(α/2)(EA/a)+(2nv)Fa2cos2(α2)12(EI/a)E24
To validate the closed form expression of Eq. (24), we compared the vertical deformation of member OA (Fig. 3) and the change in radius to the results from a finite element model in ABAQUS. The results are shown in Tables 3 and 4 below for K=652nN.nm-1 and C=1.42 nN.nm.rad-2. The accuracy obtained is excellent.
Lattice Translational Index, n
Molecular/Structural Mechanics
ABAQUS
4
1.7686
1.7686
8
1.7500
1.7500
12
1.7461
1.7461
16
1.7446
1.7446
Table 3.
Vertical deformation of member OA (x10-3).
Lattice Translational Index, n
Molecular/Structural Mechanics
ABAQUS
4
-0.53131
-0.53136
8
-0.96001
-0.96001
12
-1.4088
-1.4088
16
-1.8634
-1.8634
Table 4.
Change in radius (x10-3).
The contribution of each of the bond stiffnesses (axial, bending and torsion) to the total vertical deformation of an armchair carbon nanotube is shown in the following example.
Two long carbon nanotubes (40 armchair carbon units) with lattice translational indices “n” equal to 4 and 16, respectively, are modeled using MSC/Nastran [15]. The tubes are supported at the bottom and subjected to tensile loading at the top.
The resulting vertical deformations are compared to the closed form solution of Eq. (24), as shown in Fig. 4. In spite of the difference in boundary conditions between the closed form solution and the finite element modeling, the errors in the results are less than 3%.
Zigzag CNTs under simple tension
Going through the same manipulations as for an armchair CNT, the vertical deformation at the end of a zigzag CNT under simple tension that accounts for bending and torsional deformations can be expressed as [1]:
Δ=(nv+1)Fa2sin2α(EI/a)[24−36γ/(1+γ)]E25
where F is the load applied at a carbon atom, nu is the number of vertical units of Fig. 3 and nv = (2nu –1). However, in this case, γ takes on the following expression:
Vertical deformations of armchair carbon nanotubes with lattice translational indices of 4 and 16, respectively. Three cases are considered: 1) bending stiffness only, 2) bending + torsion and 3) axial+bending+torsional.
Zigzag CNTs under simple torsion
When the torsional stiffness of 3D frame elements is included, in addition to the axial and bending deformations, the tangential deformation ΔB of a zigzag CNT under simple torsion takes the following form [2]:
ΔB=NTKNDTKDE27
where N and D are 6x1 vectors function of n, the lattice translational index, and a, the bond length. KN and KD are 6x1 vectors of structural stiffnesses. These vectors could be expressed as:
Thus, in this case of negligible torsional stiffness, the axial and bending deformations are decoupled.
6.2. Elastic Modulus of an Armchair CNT
Similar to the previous work by Kasti [1], an elastic modulus Ys (Tpa.nm) could be defined that doesn’t include the thickness of CNTs, and is expressed as:
Ys=Ft2πR(Δ/L)=2L3aΔE28
where Ft (equal to 2nF with F=1) is the total load applied, L and R are the length and radius of the tube, respectively.
For a finite length cylinder, the elastic modulus obtained from the closed form expressions of Eq. (28) and ABAQUS are compared in Table 5 for a lattice translational index of 4 and variable tube length. The following values of stiffnesses were used:
The closed form results compare very well with the values from ABAQUS.
LatticeTranslationalIndex, n
nv=40,60(Eq. 28)
nv=40(ABAQUS)
nv=60(ABAQUS)
4
358.12
358.284
358.307
Table 5.
Elastic Modulus (nN.nm-1).
Zigzag CNTs under simple torsion
Similar to the definition of an elastic modulus Ys (Tpa.nm) that doesn’t include the thickness of CNTs [1], an elastic shear modulus Gs (Tpa.nm) could be defined as [2]:
Gs=T(φL)*2πR3=nFR(φL)*2πR3=n(φL)*2πR2E29
where T is the torsional moment applied to the tube, F is the tangential load applied to a single carbon atom which can be taken as unity, L and R are the length and radius of the tube, respectively.
Relations between the structural bending stiffness EI/a and the molecular bond bending stiffness C for carbon nanotubes were derived for the cases of simple tension and torsion. In addition, expressions for the deformations and “Young’s moduli” of these nanotubes were presented that account for the axial, bending and torsional effects.
References
1.KastiN.2007Zigzag Carbon Nanotubes: Molecular/Structural Mechanics and the Finite Element MethodInt. J. of Solids and Structures4469146929
3.ChangT.GaoH.2003Size-dependent elastic properties of a single-walled carbon nanotube via a molecular mechanics modelJ. of the Mechanics and Physics of Solids5110591074
4.IjimaS.1991Helical microtubules of graphite CarbonNature3545658
5.LiC.ChouT.2003A structural mechanics approach for the analysis of carbon nanotubesInt. J. of Solids and Structures4024872499
6.ToC.2006Bending and shear moduli of single-walled carbon nanotubesFinite Elements in Analysis and Design42404413
7.TserpesK. I.PapanikosP.2005Finite element modeling of single-walled carbon nanotubesComposites B36468477
8.XiaJ. R.GamaB. A.GillespieJ. W. Jr.2005An analytical molecular structural mechanics model for the mechanical properties of carbon nanotubesInt. J. of Solids and Structures4230753092
9.ChenW-H.ChengH-C.LiuY-L.2010Radial mechanical properties of single-walled carbon nanotubes using modified molecular structure mechanicsComputational Materials Science47985993
10.LiC.ChouT-W.2004Elastic Properties of single-walled carbon nanotubes in transverse directionsPhysical Reviews B697073401
11.WangX-F.XuZ-J.ZhuZ-Y.2007Reversible mechanical bistability of carbon nanotubes under radial compressionChemical Physics334144147
12.ChangT.GengJ.GuoX.2006Proceedings of the Royal Society A4622523
13.GallagherR.1975Finite Element Analysis: Prentice Hall
14.ABAQUS- Finite Element Analysis Program.
15.MSC/Nastran- Finite Element Analysis Program.
Written By
Najib A. Kasti
Submitted: 09 March 2012Published: 27 February 2013