As stated, the main challanges with modelling fracture mechanics problems, via Finite-Element [FE] and Boundary-Element [BE] analysis occur from the presence of an infinite-stress at the crack-tip. Sufficient accuracy can typically not be achieved, even by utilizing a very fine-mesh or applying higher-order elements. Due to the fact that the nature of the singularity of the crack-tip is known, unique crack-tip elements can and should be modelled that have built-in singularity. These elements have previously been utilized and modelled by others, via Finite Element Analysis and other methods to determine the Stress Intensity Factor. Other researchers, in the recent past, have clearly shown that the crack-tip elements are actually unnecessary, due to the fact that the discussed singularity can be easily incorporated by moving critical nodes on the isoparametric quadratic elements, so that the midpoint of the subject element is typically moved to the quarter point closer to the crack-tip, in order to create the desired stress singularity. Through the useage of these elements, the amount of mesh-refinement that is needed and necessary at the crack-tip, in order to obtain accurate values of the stress intensity factor, is significantly reduced.
How much mesh-refinement for three dimensional (3-D) structural-components is typically suggested at or near the crack-tip to generate verry accurate results. Are you utilizing any other specific methods to refine the mesh even further that work for you?