The formulation and finite element implementation of a finite deformation continuum theory for the mechanics of crystalline sheets is described. This theory generalizes standard crystal elasticity to curved monolayer lattices by means of the exponential Cauchy-Born rule. The constitutive model for a two-dimensional continuum deforming in three dimensions (a surface) is written explicitly in terms of the underlying atomistic model. The resulting hyper-elastic potential depends on the stretch and the curvature of the surface, as well as on internal elastic variables describing the rearrangements of the crystal within the unit cell. Coarse grained calculations of carbon nanotubes (CNTs) are performed by discretizing this continuum mechanics theory by finite elements. A smooth discrete representation of the surface is required, and subdivision finite elements, proposed for thin-shell analysis, are used. A detailed set of numerical experiments, in which the continuum/finite element solutions are compared to the corresponding full atomistic calculations of CNTs, involving very large deformations and geometric instabilities, demonstrates the accuracy of the proposed approach. Simulations for large multi-million systems illustrate the computational savings which can be achieved.
CitationArroyo, M.; Belytschko, T. Finite element methods for the non-linear mechanics of crystalline sheets and nanotubes. "International journal for numerical methods in engineering", Gener 2004, vol. 59, núm. 3, The definitive version is available at http://www3.interscience.wiley.com/journal/106569279/abstract, p. 419-456.
All rights reserved. This work is protected by the corresponding intellectual and industrial property rights. Without prejudice to any existing legal exemptions, reproduction, distribution, public communication or transformation of this work are prohibited without permission of the copyright holder. If you wish to make any use of the work not provided for in the law, please contact: email@example.com