Article Preview
TopIntroduction
Since first found in the form of multi-walled carbon nanotubes (MWNTs) by Iijima (1991) carbon nanotubes (CNT) have gained overwhelming attention in many areas of science and engineering for their outstanding and unique properties. CNTs are small size, low density, high stiffness, and high strength materials. CNTs can be used to produce, for example, enhanced composite materials for efficient heat removal, and precise drug delivery. These and other potential applications of CNT have provided strong drivers to develop numerical methods to better predict their mechanical behavior.
Two major types of computational approach can be recognized in the literature. They are the classical molecular dynamics (MD) methods (Iijima et al., 1996) and those that link atomistic constitutive laws to the finite element framework (Tadmor et al., 1996; Miller & Tadmor, 2009; Arroyo & Belytschko, 2002; Pantano et al., 2004; Garg et al., 2007; Tadmor et al., 1999). It is well known that the MD method excels at modeling structural details of a crystal at the lattice level. Its computational cost, however, can be prohibitive for large systems. For example (Tadmor et al., 1999), a MD simulation of a system with the number of atoms in the order of , corresponding to less than a few hundred nanometers in dimension, lends itself to a significant challenge to computing capacity. Given the fact that the mean diameter of an MWNT is in the range of 5 to 100 nm and the length is about several micrometers (Tadmor et al., 1999), the MD method is less than an ideal candidate for routine use to simulate the mechanical and dynamic behavior of CNTs. On the other hand, the results obtained by using the MD method have suggested that the majority of the lattice deforms smoothly and closely obeys continuum elasticity. Yakobson et al. (1996) used the MD method to study CNTs under three basic mechanical loads: axial compression, bending, and torsion. They reported a remarkable synergism between the results of MD and those of macroscopic structural mechanics. Arroyo et al. (2005) used a two-dimensional atomic rope deformation case to study the convergence property of an atomistic based finite element method (or continuum method). It was shown that, with increasing number of atoms simulated, the strain energy of the rope obtained from the continuum method converged to that from a parent atomic model. Consequently, continuum methods with constitutive laws constructed upon atomistic energy functions can play an important role in CNT simulations.
In what follows we present a brief review of two types of reduced-order continuum schemes, i.e. the general continuum method and the elastic model-based continuum method used in the prediction of the mechanical behavior of CNTs. For a more general review readers are referred to Miller and Tadmor (2009) where fourteen different reduced-order schemes were compared for a static Lomer dislocation dipole problem in face centered cubic aluminum.