We give an analysis of a continuation algorithm for the numerical solution of the force-based quasicontinuum equations. The approximate solution of the force-based quasicontinuum equations is computed by an iterative method using an energy-based quasicontinuum approximation as the preconditioner. The analysis presented in this paper is used to determine an efficient strategy for the parameter step size and number of iterations at each parameter value to achieve a solution to a required tolerance. We present computational results for the deformation of a Lennard-Jones chain under tension to demonstrate the necessity of carefully applying continuation to ensure that the computed solution remains in the domain of convergence of the iterative method as the parameter is increased. These results exhibit fracture before the actual load limit if the parameter step size is too large.
Bibliographical noteFunding Information:
This work was supported in part by DMS-0304326, by the University of Minnesota Supercomputing Institute and by the Department of Energy under Award Number DE-FG02-05ER25706.
- Atomistic to continuum