Finite Elements & Fracture Mechanics

Stress Intensity Factor

Note: In a Zencrack analysis the j-integral is always calculated (if available in the interfaced f.e. code) and is the default method for providing fracture mechanics parameters to drive crack growth calculations. However, the displacement based stress intensity factors described here can be used to drive crack growth calculations if desired. For Abaqus and Ansys interfaces the interaction integral calculations for stress intensity factors may also be requested and used to drive crack growth calculations.

The stress intensity factor distribution along each crack front is always calculated when the results of a finite element analysis are processed by Zencrack by converting nodal displacements to stress intensity factors and energy release rates. This is a process which is valid for linear elastic isotropic materials. If there are multiple loads steps in the f.e. analysis, for example due to the application of a thermal transient, then results for each step are processed giving the variation of stress intensity factors through the load history.

Zencrack triad

Representation of a local triad system at a crack front location

The calculations are based on the relative opening, sliding and tearing displacements derived from an orthogonal set of axes at each crack front node, as shown at one crack position in the image for a symmetry model. These relative displacements are used to calculate the stress intensity factors using equations derived from the Westergaard solution for the stress field around a crack tip. The equations that are used are valid for linear elastic isotropic materials.

The relative opening displacement is also an important quantity in determining a local "open" or "closed" status for a crack front node. This status has implications for crack growth integration.

If the analysis also includes a j-integral evaluation, a conversion of the j-integral to mode I stress intensity factor is also reported. This is strictly only valid if there are no mode II or mode III effects present.

For Abaqus and Ansys interfaces the interaction integral calculations for stress intensity factors may also be requested.

The following may be calculated:

  • Ki, Kii and Kiii from nodal displacements using crack tip opening displacements (CTOD) - only valid for linear elastic isotropic materials
  • Ki, Kii and Kiii from integration integral calculations (Abaqus and Ansys)
  • equivalent energy release rate from displacement-based or interaction integral Ki, Kii and Kiii
  • Ki from j-integral energy release rate (mode I conversion)
  • relative opening displacements to define a proper "open" or "closed" status at crack front and crack face nodes
  • the local crack growth direction from the local maximum energy release rate
  • the local crack growth direction from the stress intensity factors.

Technical

Finite Elements & Fracture Mechanics

More in this category