The energy-based growth formulation and accompanying simulation technique introduced in Part I of this series is generalized in this work to predict arbitrary, mixed-mode, non-planar crack evolution. The implementation uses a novel basis-function approach to generate a crack extension expression, rather than relying on the local, point-by-point approach described in Part I. The basis-function expression dampens the effect of numerical noise on crack growth predictions that could produce numerically unstable simulation results. Two simulations are presented to demonstrate the technique’s ability to capture both general non-planar behavior, as well as local mixed-mode phenomena, e.g. “factory-roof” formation, along the crack front.