not knowing what a quad4 or a quad8 is (possibly NASTRAN plate elements), and guessing that you're calculating stress intensity from a path integral, I'd guess that the difference is the stress shape function (quad8 is parabolic and quad4 linear). and this would also be affected by mesh density (i'd expect that quad4 elements would have a denser mesh than quad8, to offset the element stress assumption).