Engineers require accurate determination of the configurational force at the crack tip for fracture fatigue analysis and accurate crack propagation. However, obtain- ing highly accurate crack tip configuration force values is challenging with numer- ical methods requiring knowledge of the stress field around the crack tip a priori. In this thesis, the symmetric interior penalty discontinuous Galerkin finite element method is combined with a residual based a posteriori error estimator which drives a hp-adaptive mesh refinement scheme to determine accurate solutions of the stress field about about the crack. This facilitates the development of a novel method to calculate the crack tip configurational force that is accurate, requires no a priori knowledge of the stress field about the crack tip with, its error bound by an error estimator which is calculated a posteriori. Benchmark values of the crack tip con- figurational force are presented for problems containing multiple mixed mode cracks in both isotropic and anisotropic materials. Additionally, the hp-adaptivity is com- bined with a mathematical analysis of the stress field at the crack tip to critique the convergence and limitations of other methods in the literature to calculate the crack tip configurational force. Two methods for staggered quasi-static crack prop- agation are also presented. An rp-adaptive method which is simple to implement and computationally inexpensive, element edges aligned with the crack propagation path with the exploitation of the discontinuous Galerkin edge sti↵ness terms exist- ing along element interfaces to propagate a crack. The second method is denoted the hpr-adaptive method which combines the accurate computation of the crack tip configuration force with r-adaptivity to produce a computationally expensive but accurate method to propagate multiple cracks simultaneously. Further, for indeter- minate systems, an average boundary condition that restrains rigid body motion and rotation is introduced to make the system determinate
