Vertical deformation of member OA (x10^{-4}).

## 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 (* = 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

*, the torsional angle*C

*and the lattice translational index*φ

*. 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*n

*and the structural bending stiffness*C

*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

*=*EI/a

*. 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.*C/2

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

*is derived for the case of simple tension. This shows that*C

*depends on the bond bending stiffness*EI/a

*and the torsional angle*C

*In the limit of an infinite tube radius, which represents a graphene sheet, we recover the previous relation, i.e.,*φ.

*tends to*EI/a

*. 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.*C/2

## 2. Bond Energies and the Finite Element Method

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*,

*), the lattice translational indices, and the bond length*m

*In general, the diameter*a.

*of a CNT is defined using the expression*d

For zigzag CNTs (Fig. 1a), the value of * m*is zero. In this case, the diameter is simply given by the formula

*is equal to*m

*. In this case, the diameter is given by the formula*n

The bond energies between carbon atoms include the stretching _{a}, bending _{b}and torsional _{t}energies. For small distortions from equilibrium, these energies could take the forms:

where * K*,

*, and*C

C

_{t}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*,

*, and*EI

*are the axial, bending and torsional stiffnesses; and*GJ

*,*u

*and*θ

*are the axial, bending and torsional deformations, respectively.*φ

### 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.

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.

## 3. Work of Chang and Gao [3]

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.

## 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 that

We start with the molecular energy expression

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:

Since both * K*and

*are conjugate to the axial deformation between carbon atoms in the energy equation, they represent the same axial stiffness. Thus,*EA/a

To determine the tube deformation, we let _{u}be the number of vertical units of Fig. 3 in a carbon nanotube and _{v}equal to _{u}_{.} The deflection at the end of the tube due to axial bond deformations can be expressed as:

### 4.2. Bond Bending – Molecular Mechanics

Due to bond bending, the total potential energy is written as:

Let _{x}be the number of units along the circumference and _{v}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:

Minimizing the total potential energy with respect to α gives

Since cosβ = -cos(π/2n)cos(α/2) and cos* φ*= tan(

*/2)/tan(*α

*), we get d*β

*/d*β

*= cos*α

*, where*φ/2

*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:

where _{o,}* H*are the initial and current heights of the tube, respectively. Differentiating Eq. (10) with respect to

α

Solving for * Δα*by minimizing the total potential energy, we get

Substituting * Δα*in Eq. (10) above, we get the following expression for the vertical deflection of an infinitely long CNT (due to bond bending only)

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

*be two atoms on the CNT with coordinates (*A

*cos(-*R

*sin(-*θ/2), R

*) and (*θ/2), 0

*cos(-*R

*)*π/n+θ/2

*sin(-*, R

*)*π/n+θ/2

*(*, a sin

*) ), respectively. The angle*α/2

*is given by*θ

2R

^{2}

(1-cosθ)=a

^{2}.

Also, let * u*be the unit vector from

*to*O

*expressed as:*A

where * i*,

*and*j

*are the unit vectors in the Cartesian coordinate system.*k

For an infinite cylinder and due to symmetry, the radial and tangential rotations at * O*and

*are zero.*A

In addition, due to multiple symmetries (for * n*=4,8,12,16,..), we assume the displacement and rotation fields at

*and*O

*take the following forms:*A

and

where ( _{o}_{A}) and ( _{o}_{A}) are the displacements and rotations vectors at * O*and

*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.*A

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:

and the corresponding elastic modulus Y_{s}

The corresponding axial deformation at the end of a zigzag CNT under simple tension is given by [1]:

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

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, ^{2}* (φ)*→ 1 and

*→*EI/a

*, which is the same result derived by Kasti [1] for the particular case of a graphene sheet.*C/2

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.

4 | 5.9844 | 5.9849 |

8 | 5.9477 | 5.9470 |

12 | 5.9358 | 5.9358 |

16 | 5.9310 | 5.9306 |

4 | -1.3538 | -1.3539 |

8 | -2.6402 | -2.6401 |

12 | -3.9383 | -3.9383 |

16 | -5.2403 | -5.2402 |

For zigzag CNT under simple tension, the equivalent stiffness is given by [1]:

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

*, which is the same result for the particular case of a graphene sheet.*C/2

For a zigzag CNTs under simple torsion, the corresponding stiffness is given by [2]:

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*→ ∞,

*→*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:

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

*=1.42 nN.nm.rad*C

^{-2}. The accuracy obtained is excellent.

4 | 1.7686 | 1.7686 |

8 | 1.7500 | 1.7500 |

12 | 1.7461 | 1.7461 |

16 | 1.7446 | 1.7446 |

4 | -0.53131 | -0.53136 |

8 | -0.96001 | -0.96001 |

12 | -1.4088 | -1.4088 |

16 | -1.8634 | -1.8634 |

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%.

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]:

where * F*is the load applied at a carbon atom,

n

_{u}is the number of vertical units of Fig. 3 and

n

_{v}= (2

n

_{u}–1). However, in this case, γ takes on the following expression:

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]:

where * N*and

*are 6x1 vectors function of*D

*, the lattice translational index, and*n

*, the bond length.*a

*and*KN

*are 6x1 vectors of structural stiffnesses. These vectors could be expressed as:*KD

* KN*= [ (EI/L)

^{3}(EI/L)

^{2}(GJ/L) (EA/L)(GJ/L)

^{2}(EA/L)(EI/L)(GJ/L) (EI/L)(GJ/L)

^{2}(EA/L)(EI/L)

^{2}]

^{T}

* KD*= [ (EI/L)

^{4}(GJ/L)(EI/L)

^{3}(EA/L)(GJ/L)

^{2}(EI/L) (EA/L)(EI/L)

^{2}(GJ/L) (EI/L)

^{2}(GJ/L)

^{2}(EA/L)(EI/L)

^{3}]

^{T}

* N*= [ N1 N2 N3 N4 N5 N6 ]

^{T,}

*[ D1 D2 D3 D4 D5 D6]*D =

^{T}

For example, for n=4,

* N*= [ -21.1721 -24.53 -0.0015 -0.0969 -0.4074 -0.0832 ]

^{T}

* D =*[ 0. –2.3062x10

^{3}-0.4408 -76.1719 497.0279 -62.4643 ]

^{T}

One point worth mentioning is that when the torsional stiffness is neglected, i.e._{B}from Eq. (27) takes the form:

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 _{s}(Tpa.nm) could be defined that doesn’t include the thickness of CNTs, and is expressed as:

where _{t}(equal to 2nF with F=1) is the total load applied, * L*and

*are the length and radius of the tube, respectively.*R

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:

* EA/a*= 652 nN.nm

^{-1},

*= 0.875 nN.nm.rad*EI/a

^{-2}and

*= 0.278 nN.nm.rad*GJ/a

^{-2.}

The closed form results compare very well with the values from ABAQUS.

_{v}=40,60 | _{v}=40 | _{v}=60 | |

4 | 358.12 | 358.284 | 358.307 |

Similar to the definition of an elastic modulus _{s}(Tpa.nm) that doesn’t include the thickness of CNTs [1], an elastic shear modulus Gs (Tpa.nm) could be defined as [2]:

where * T*is the torsional moment applied to the tube,

*is the tangential load applied to a single carbon atom which can be taken as unity,*F

*and*L

*are the length and radius of the tube, respectively.*R

## 7. Conclusions

Relations between the structural bending stiffness * EI/a*and the molecular bond bending stiffness

*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.*C