Highlights• A framework for concurrent atomistic/continuum simulation of graphene is developed.• The TBDM is re-formulated based on the Nosé-Hoover thermostat.• The AIREBO potential is used, which is carefully weighted in the bridging zone.• The Q1STs solid-shell element is modified to include the energy scaling function.• A thermodynamically rigorous high-order continuum constitutive model is parameterized.
AbstractRecently developed three-layer-mesh bridging domain method (TBDM) enhanced the conventional bridging domain method (BDM) by (1) mitigating the temperature cooling effect on the atoms in the bridging domain, and (2) employing a mesh-independent physics-based discrimination between thermal and mechanical atomic motions. In this paper, we present the new enhancements for the TBDM to achieve an appropriate framework for concurrently coupled atomistic-continuum simulation of graphene. To capture the mechanical behavior of graphene accurately, we employed the adaptive intermolecular reactive bond order (AIREBO) potential in the atomistic model, which is carefully weighted by the atomic energy scaling function in the bridging zone. In the continuum model, a thermodynamically rigorous high-order continuum description, considering the symmetries of graphene, is used which is parameterized using full molecular dynamics (MD) simulations. To accurately capture the bending behavior of graphene, a recently developed explicit finite-deformation solid-shell element is used to discretize the continuum domain, and its formulation is modified to include the continuum energy scaling function. To achieve realistic constant-temperature condition (canonical ensemble), the Nosé-Hoover thermostat is used in the full MD domain and also as local thermostats in the bridging domain. 5-value Gear predictor-corrector time integrator is implemented, which is well-suited to be used with the Nosé-Hoover thermostat. Accordingly, the TBDM formulation is modified to work with this time integrator. Some modifications are also made in the TBDM formulation to increase the robustness of the multiscale simulations. Finally, the effectiveness of the proposed A. Sadeghirad et al. / Comput. Methods Appl. Mech. Engrg. 294 (2015) 278-298 279 multiscale method for graphene is demonstrated by running in-plane shear, out-of-plane bending, and nanoindentation simulations and comparing the results with those obtained from full MD and full finite-element simulations.