Finite element methods for the non-linear mechanics of crystalline sheets and nanotubes

The formulation and nite element implementation of a nite 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 twodimensional 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 nite elements. A smooth discrete representation of the surface is required, and subdivision nite elements, proposed for thin-shell analysis, are used. A detailed set of numerical experiments, in which the continuum/ nite 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