A multiscale FEM-MD coupling method for investigation into atomistic-scale deformation mechanisms of nanocrystalline metals under continuum-scale deformation

This study aims to develop a multiscale bridging method for investigating nanocrystalline metals based on macro-scale deformation. For this purpose, we propose a hierarchical multiscale computational method that can focus on some of the elements in a finite element model for scale bridging to atomistic-scale models. This method assumes that atomistic-scale nanocrystalline models are related to the integration points in a finite element and deform based on the macro-scale deformation. Nanocrystalline aluminum was chosen for the validation of the multiscale method. The finite element method (FEM) and the molecular dynamics (MD) method were used for continuum-scale and atomistic-scale simulations, respectively. We utilized the notion of the CauchyBorn rule (CBR) for communicating deformation information from the continuum scale to the atomistic scale. We studied three different cases with two nanocrystalline models and two loading cases to compare differences resulting from crystal structures and loading. Based on the crystal structure change during relaxation, nonequilibrium grain boundaries (NEGBs) were shown to play a role as deformation mechanisms in the plastic regime and induce the onset and migration of crystal defects, including deformation twins, as reported in the experiment. Furthermore, the crystal orientation dependence of the onset of crystal defects was confirmed by the comparison of the results from the two different nanocrystalline models. The qualitative agreement of the results with experimental observations is also confirmed. The proposed ‘FEM-MD’ method can bridge a large-scale gap, for example, from a nano-scale to a continuum-scale such that an MD model can be coupled to a millimeter or centimeter scale compared to other embedding methods. The present method is ideal for investigating the dislocation behavior of nanocrystalline materials, which contain multi-grained nanostructure at finite temperature, undergoing various loading scenarios at the macro-scale.