Systems and methods are disclosed for determining stress intensity factors. In one or more embodiments, the method can include the steps of defining a crack tip zone about one or more crack tips of one or more arbitrarily shaped cracks in an arbitrarily shaped solid. The one or more crack tip zones can be constrained within a finite element model representation mesh of the arbitrarily shaped solid to provide one or more constrained crack tip zones. The combination of the finite element model representation mesh and the one or more constrained crack tip zones can be processed to determine the stress intensity factor for each of the one or more arbitrarily shaped cracks.