Mathematical Modelling of CNT to Derive its Different Parameter to Be used in Finite Element Analysis.

DOI : 10.17577/IJERTCONV3IS17008

Download Full-Text PDF Cite this Publication

Text Only Version

Mathematical Modelling of CNT to Derive its Different Parameter to Be used in Finite Element Analysis.

P K Mehta

Dept. of Mechanical Engineering,

City Engineering College, Bangalore, Karnataka, India

Abstract: Carbon Nano Tubes have significant mechanical and electrical properties which have been subjected to many theoretical and experimental studies. Experimental studies are quite cumbersome and expensive, hence researchers are moving towards theoretical study of CNT. In this connection, modeling and analysis of CNT has given a tremendous help because of its result approximate to exact result. For modeling of CNT in FEM software, we need different dimensional parameter. In this paper, effort has been made to extract different dimensional parameter for CNT which can be used to model the CNT in FEM Software.

Keywords- CNT, Chirality, FEM.


    Nanotechnology refer to the understanding and control of matter at the atomic levels, at the length scale of approximately 1 to 100 nanometers, where unique phenomena enable novel applications. Nanotechnologies are the design, characterization, production and application of structure, devices and systems by controlling shape and size at nanometer scale.

    One of the fundamental components in nanotechnology is Carbon Nanotubes made of Carbon atoms arranged in a hexagonal lattice, whose diameter ranges from one to tens of nanometer and that are from hundreds up to thousands nanometers in lengths.

    CNTs have incredible physical properties, both mechanical and electrical, and they are the basis of most of the new nano technological innovations. FEM and Molecular mechanics share a common ground of energy minimization, with respect to the nodes and to the discrete atoms, respectively.


    A SWCNT can be thought of as a sheet of graphite (a hexagonal lattice of carbon) rolled into an endless cylinder. It can be constructed by wrapping up a single sheet of graphite such that two equivalent sites of the hexagonal lattice coincide. The wrapping vector, C, which defines the relative location of the two sites, is specified by a pair of integer (n,m) that relate C to the two unit vectors, a1 and a2 and is given by the following equation.

    Ch = na1 +ma2 (1)

    Where a1 and a2 are unit vector,

    Where acc =distance between nearest neighbour carbon atom.

    Let a be the unit length of grapheme honeycomb lattice, then


    A tube is called armchair if n equals m and zigzag in the case m=0. All other tubes are of the chiral type and have a finite wrapping angle with 0< < 30

    Fig 1

    If we equate the length of chirality vector to the circumference of the tube, we will get the diameter of the tube.


    The magnitude of the chirality vector can also be determined in terms of the two lattice parameters, n and m, and the nearest-neighbor carbon distance:


    For an arm chair CNT (n,n) , Ch = 3nacc; For a zigzag CNT Ch = 3nacc

    Thus the diameter of the tube associated with an arbitrary chirality vector can be defined as

    Fig 2

    The translational vector T perpendicular to the chiral vector Ch is expressed as




    The chiral angle ,, determines the degree of helicity to the lattice, with a zig-zag (n; 0) lattice defined as =0.This value can also be determined in terms of the lattice parameters.

    Tangent of the chiral angle


    Where, d is is the highest common divisor of n and m.

    The translational vector is


    The length of translational vector is the NT unit cell length

    L along the nanotube axis direction:


    For arm chair CNT

    For ZigZag CNT

    Number of hexagons N in a nanotube unit cell is


    For arm chair CNT , m=n

    For ZigZag CNT, m=0


CNTs carbon atoms are bonded together with covalent bonds forming an hexagonal lattice. These bonds have a characteristic bond length a CC and bond angle in the 3D space. The displacement of individual atoms under an external force is constrained by the bonds. Therefore, the total deformation of the nanotube is the result of the interactions between the bonds. By considering the bonds as connecting load-carrying elements, and the atoms as joints of the connecting elements, CNTs may be simulated as space-frame structures.

Fig 3

By treating CNTs as space-frame structures, their mechanical behavior can be analyzed using classical structural mechanics methods.

3.1 Finite element modeling with Representative Volume Element (RVE)

In RVE approach, a single nanotube with surrounding matrix material is modeled, with properly applied boundary and interface conditions to account for the effects of the surrounding materials it is possible to study the

nanocomposite based on the single Representative volume element. This RVE model can be employed to study the interactions of the nanotube with the matrix. It is possible to model and evaluate the effective material properties of the nanocomposite.

According to the shape of the cross section, three types of representative volume elements are possible : cylindrical, square, and hexagonal RVEs shown in

Fig. 4 . Specifically, the cylindrical RVE is used to model carbon nanotubes with different diameters. The square RVE is applied when carbon nanotubes are arranged evenly in a square pattern, while the hexagonal RVE is adopted when carbon nanotubes are in a hexagonal pattern. Compared to the square RVE, the cylindrical RVE tends to overestimate the effective Youngs module of carbon nanotube/polymer composites. This can be explained by the fact that a cylindrical RVE overestimates the volume fraction of carbon nanotube due to the negligence of the small amount of polymer matrix.

Fig.4 Three types of Representative Volume Elements (RVEs) (a) Cylindrical RVE (b) Square (c) Hexagonal


[1]. Rappe AK, Casemit CJ, Colwell KS, Goddard WA, Skiff WM. UFF, a full periodic-table force-field for molecular mechanics and molecular dynamics simulations. J Am Chem Soc 1992;114:10024 35.

[2]. K.I. Tserpes, P. Papanikos, Composites B 36, 468 (2005)

  1. Chang T. Gao H Size dependent elastic properties of a single walled carbon nanotube via a molecular mechanics model. Journal of the Mechanics and Physics of Solids 2003. 51:1059 1074

  2. Herron N, Thorn DL. Nanoparticles: uses and relationships to molecular clusters. Adv Mater 1998;10:117384.

  3. Song SY, Youn JR. Modelling of effective elastic properties of polymer based carbon nanotube composites. Polymer 2006; 47, 1741- 1748.

  4. M. Endo, T. Hayashi, Y.A. Kim, M. Terrones, M.S. Dresselhaus, Applications of carbon nanotubes in the twenty-first century, Phil. Trans. R. Soc. Lond. A, 362, 2223-2238, 2004.

  5. Edelstein, A.S. and Cammarata, R.C.: Nanomaterials: Synthesis, Properties andApplications. Bristol: Institute of Physics Publishing, 1996.

  6. Salvetat JD., Rubio A., Mechanical properties of carbon nanotubes: a fiber digest for beginners, Carbon, 2002; 40:172934.

  7. Noor, A.K.; Anderson, M.S., and Greene, W.H.: Continuum Models for Beam- and Platelike Lattice Structures. AIAA Journal. Vol. 16, No. 12, 1978, pp. 1219-1228.

Leave a Reply