The J -integral depends on the contour you are using to calculate it. You experience convergence if the plastic zone around crack tip is fully included in the domain you are considering. The mesh around the tip should be refined in order to accurately model the singularity and the stress/strain gradients. I suggest this report: