Stress intensity factor calculation based on the work of external forces

The external forces method is a numerical method for K calculation based on the finite element method. It uses the work of the external forces W for the calculation of the energy release rate and is particularly advantageous when that forces are applied far from the crack front. The method was applied to a corner crack geometry with the objective of studying its accuracy. Good results were obtained for a wide range of virtual crack displacements (0.03% < Δa/a < 6%) considering 4 values of W along with a polynomial regression of order 3. For that choice of parameters the inaccuracy of K is mainly due to FEM errors. A great sensitivity of K to FEM errors was observed, however accurate values of K were obtained, with errors lower than 2 percent. So, the use of the external forces method for the calculation of K is recommended, considering its simplicity and accuracy.

[1]  Z. Bažant,et al.  Surface singularity and crack propagation , 1979 .

[2]  S. Chan,et al.  On the Finite Element Method in Linear Fracture Mechanics , 1970 .

[3]  D. M. Parks A stiffness derivative finite element technique for determination of crack tip stress intensity factors , 1974 .

[4]  H. G. Delorenzi Energy release rate calculations by the finite element method , 1985 .

[5]  H. deLorenzi,et al.  On the energy release rate and the J-integral for 3-D crack configurations , 1982 .

[6]  T. K. Hellen On the method of virtual crack extensions , 1975 .

[7]  M. Kanninen,et al.  A finite element calculation of stress intensity factors by a modified crack closure integral , 1977 .

[8]  Ray Kai Leung Su,et al.  Analytical solution for mode I crack orthogonal to free surface , 1996 .

[9]  A. Ingraffea,et al.  Stress‐intensity factor computation in three dimensions with quarter‐point elements , 1980 .

[10]  J. Benthem State of stress at the vertex of a quarter-infinite crack in a half-space , 1977 .

[11]  D. Sherman,et al.  On quarter-point three-dimensional finite elements in linear elastic fracture mechanics , 1989 .

[12]  J. Rice A path-independent integral and the approximate analysis of strain , 1968 .

[13]  K. N. Shivakumar,et al.  A virtual crack-closure technique for calculating stress intensity factors for cracked three dimensional bodies , 1988, International Journal of Fracture.

[14]  G. De Roeck,et al.  Strain energy release rate formulae for 3D finite element , 1995 .

[15]  Brian Moran,et al.  Energy release rate along a three-dimensional crack front in a thermally stressed body , 1986, International Journal of Fracture.

[16]  A. Needleman,et al.  A COMPARISON OF METHODS FOR CALCULATING ENERGY RELEASE RATES , 1985 .

[17]  A. Bakker,et al.  THREE‐DIMENSIONAL CONSTRAINT EFFECTS ON STRESS INTENSITY DISTRIBUTIONS IN PLATE GEOMETRIES WITH THROUGH‐THICKNESS CRACKS , 1992 .

[18]  D. M. Tracey,et al.  The natural isoparametric triangle versus collapsed quadrilateral for elastic crack analysis , 1976, International Journal of Fracture.

[19]  T. Sato,et al.  THREE-DIMENSIONAL J-INTEGRAL CALCULATIONS OF PART-THROUGH SURFACE CRACK PROBLEMS , 1983 .