Simulation of instabilities in thin nanostructures by a perturbation approach

Y. Cong (1) , Julien Yvonnet (2) , H. Zahrouni (1)

1 LEM3 - Laboratoire d'Etude des Microstructures et de Mécanique des Matériaux
2 MSME - Laboratoire de Modélisation et Simulation Multi Echelle

A new numerical method is proposed to simulate instabilities in thin atomistic structures in quasi-static regime. In contrast with previous approaches based on energy minimization or Newton-Raphson methods, the present technique uses a series expansion of atomistic displacements with respect to a loading path parameter, truncated at high orders. The nonlinear set of equations de ned by minimizing the potential energy of the discrete system with respect to nuclei positions is then transformed into a sequence of linear sets of equations, which can be solved e ciently. The solution can be described along very large loading steps without correction, resulting in a signi cant reduction of matrices to be inverted. Finally, the treatment of limit points and snap-back/snap-through arising when instabilities occur is simpli ed due to a continuous description with respect tothe loading path parameter. The method is applied to the analysis of single carbon atom layers nanostructures like graphene sheets or nanotubes in traction or compression regimes. Accuracy and e ciency of the technique is demonstrated by comparisons with iterative Newton procedures.

Domaines: Mécanique [physics.med-ph]

Y. Cong, Julien Yvonnet, H. Zahrouni. Simulation of instabilities in thin nanostructures by a perturbation approach. Computational Mechanics, 2014, 53 (4), pp.739-750. ⟨10.1007/s00466-013-0927-7⟩. ⟨hal-00958886⟩