I need to calculate the SIF in ABAQUS explicit. It seems contour integral and XFEM stationary crack are not available in explicit module to calculate SIF. Is there any other way to get this parameter in ABAQUS/explicit?
Thanks for your reply. But, my question is I don't know if there is any tool or procedure to obtain the stress intensity in Abaqus/Explicit module! However, for Abaqus/Standard one can use J-integral to obtain the stress intensities.
I do not think there is a direct tool, but there are methods to extract SIFs from the stress field - this usually requires a very good resolution of the stresses, probably you should use quadratic quarter-node elements.
If you calculate the stresses on a line or plane remote from the hole, you can use these and compare with tabulated analytical solutions for holes.
But real cracks are not holes, and needs to be repaired more often. Real cracks have infinite stresses close, and are analysed with LEFM, ie J-integrals.
Many codes say that inclusions can be regarded as holes relying on that the steel is ductile.
To obtain J-integrals similar to tabulated K-values, I think it must be a fine grid if you do not use elements that gives the 1/r dependency. From Dr Baekers answer, it would be possible.
What is your purpose with the J-integrals? Is it to compare values at different meshes?
I am not sure if ABAQUS can do that. However, we have developed a simple mathematical framework to extract the SIFs for pure mode I, II, and mixed mode problems. It works for both the isotropic and orthotropic mediums and also holds for homogeneous and non-homogeneous materials. It can be implemented in a MATLAB, C++, etc environment.
Link to the paper: http://www.sciencedirect.com/science/article/pii/S0013794414002239