TY - JOUR
T1 - Accuracy of quasicontinuum approximations near instabilities
AU - Dobson, M.
AU - Luskin, M.
AU - Ortner, C.
PY - 2010/10
Y1 - 2010/10
N2 - The formation and motion of lattice defects such as cracks, dislocations, or grain boundaries, occurs when the lattice configuration loses stability, that is, when an eigenvalue of the Hessian of the lattice energy functional becomes negative. When the atomistic energy is approximated by a hybrid energy that couples atomistic and continuum models, the accuracy of the approximation can only be guaranteed near deformations where both the atomistic energy as well as the hybrid energy are stable. We propose, therefore, that it is essential for the evaluation of the predictive capability of atomistic-to-continuum coupling methods near instabilities that a theoretical analysis be performed, at least for some representative model problems, that determines whether the hybrid energies remain stable up to the onset of instability of the atomistic energy. We formulate a one-dimensional model problem with nearest and next-nearest neighbour interactions and use rigorous analysis, asymptotic methods, and numerical experiments to obtain such sharp stability estimates for the basic conservative quasicontinuum (QC) approximations. Our results show that the consistent quasi-nonlocal QC approximation correctly reproduces the stability of the atomistic system, whereas the inconsistent energy-based QC approximation incorrectly predicts instability at a significantly reduced applied load that we describe by an analytic criterion in terms of the derivatives of the atomistic potential.
AB - The formation and motion of lattice defects such as cracks, dislocations, or grain boundaries, occurs when the lattice configuration loses stability, that is, when an eigenvalue of the Hessian of the lattice energy functional becomes negative. When the atomistic energy is approximated by a hybrid energy that couples atomistic and continuum models, the accuracy of the approximation can only be guaranteed near deformations where both the atomistic energy as well as the hybrid energy are stable. We propose, therefore, that it is essential for the evaluation of the predictive capability of atomistic-to-continuum coupling methods near instabilities that a theoretical analysis be performed, at least for some representative model problems, that determines whether the hybrid energies remain stable up to the onset of instability of the atomistic energy. We formulate a one-dimensional model problem with nearest and next-nearest neighbour interactions and use rigorous analysis, asymptotic methods, and numerical experiments to obtain such sharp stability estimates for the basic conservative quasicontinuum (QC) approximations. Our results show that the consistent quasi-nonlocal QC approximation correctly reproduces the stability of the atomistic system, whereas the inconsistent energy-based QC approximation incorrectly predicts instability at a significantly reduced applied load that we describe by an analytic criterion in terms of the derivatives of the atomistic potential.
KW - Atomistic-to-continuum coupling
KW - Defects
KW - Fracture
KW - Quasicontinuum method
KW - Sharp stability estimates
UR - http://www.scopus.com/inward/record.url?scp=79955654949&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=79955654949&partnerID=8YFLogxK
U2 - 10.1016/j.jmps.2010.06.011
DO - 10.1016/j.jmps.2010.06.011
M3 - Article
AN - SCOPUS:79955654949
SN - 0022-5096
VL - 58
SP - 1741
EP - 1757
JO - Journal of the Mechanics and Physics of Solids
JF - Journal of the Mechanics and Physics of Solids
IS - 10
ER -