A multiscale method is presented which couples a molecular dynamics approach for describing fracture at the crack tip with an extended finite element method for discretizing the remainder of the domain. After recalling the basic equations of molecular dynamics and continuum mechanics the discretization is discussed for the continuum subdomain where the partition-of-unity property of finite element shape functions is used, since in this fashion the crack in the wake of its tip is naturally modelled as a traction-free discontinuity. Next, the zonal coupling method between the atomistic and continuum models is described, including an assessment of the energy transfer between both domains for a one-dimensional problem. It is discussed how the stress has been computed in the atomistic subdomain, and a two-dimensional computation is presented of dynamic fracture using the coupled model. The result shows multiple branching, which is reminiscent of recent results from simulations on dynamic fracture using cohesive-zone models.
|Title of host publication||Computer Methods in Mechanics Lectures of the CMM 2009|
|Editors||M. Kuczma, K. Wilmanski|
|Place of Publication||Berlin|
|Publication status||Published - 2010|
|Name||Advanced structured materials|
Aubertin, P., Rethore, J., & Borst, de, R. (2010). A multiscale molecular dynamics / extended finite element method for dynamic fracture. In M. Kuczma, & K. Wilmanski (Eds.), Computer Methods in Mechanics Lectures of the CMM 2009 (pp. 211-237). (Advanced structured materials). Berlin: Springer. https://doi.org/10.1007/978-3-642-05241-5_12