Dynaflow Research Group

  • Engineering
    • Projects
    • Mechanical
    • CFD & Surge Analysis
    • Pulsations & Vibrations
    • Fiberglass Engineering
    • Geothermal Wells
  • Software
    • BOSfluids
    • BOSpulse
    • FEPipe
    • NozzlePRO
    • FEATools
    • Jive
    • ISOtracer
    • BOSview
    • BOS B31
  • Training
  • Career
  • News
  • Contact
Skip to content

Dynaflow Research Group

Advanced Solutions in Engineering
  • Engineering
    • Projects
    • Mechanical
    • CFD & Surge Analysis
    • Pulsations & Vibrations
    • Fiberglass Engineering
    • Geothermal Wells
  • Software
    • Flow Analysis
      • BOSfluids
      • BOSpulse
    • Structural Analysis
      • FEPipe
      • NozzlePRO
      • FEATools
    • Academic
      • Jive
    • Tools
      • ISOtracer
      • BOSview
      • BOS B31
  • Training
  • Career
  • News
  • contact

The J-integral Approach in Elastic-Plastic Fatigue Analysis

7 Jan 2022

The J-integral Approach in Elastic-Plastic Fatigue Analysis by Tammo Dukker

In a fatigue analysis one would like to investigate the susceptibility of a structural component to crack propagation under cyclic loading to determine the number of cycles until a failure occurs. Common techniques applied in a fatigue analysis follow from the results of a linear elastic analysis with limited plasticity assumed. For cracks in critical components under high stress however, high plasticity could be present and an elastic plastic analysis will provide a more accurate estimate. Rice, J (1968) presented a novel technique within this regime by introducing a path independent integral around a crack tip: the J-integral [1]. This J-integral is representative of a potential energy release with respect to crack extension and was coined by Rice as [1]:

\[J = \int_\Gamma (Wdy-\vec{T}\cdot\frac{\partial\vec{u}}{dx}ds) \]

Where $\Gamma$ is a curve surrounding the notch tip starting on the bottom crack surface and ending on the top crack surface as in figure 1, $\vec{T}$ is the traction vector defined with the outward normal from $\Gamma$: $T_i = σ_{ij}n_j$, $\vec{u}$ is the displacement vector and W is the strain-energy density defined by

\[W = \int_0^\epsilon \sigma_{ij}d\epsilon_{ij}\]

where $\epsilon_{ij}$ is the infinitesimal strain tensor. In essence, as Bowles (1980) notes, the J-integral can be viewed as the plastic version of the linear stress intensification factor K used in linear fatigue analysis [2]. In the linear region they are related by

\[\Delta J = \frac{\Delta K ^2(1-\nu^2)}{E}\]

Where E and $\nu$ are the young’s modulus and Poisson’s ratio respectively.

Figure 1: Crack with a path Γ of the J-integral
Figure 1: Crack with a path Γ of the J-integral
The power of the J-integral lies in its path independence. Any path around the crack tip as defined in figure 1 yields the same value for J. Thus the stress field around a crack tip delivers one unique value for the J-integral at that crack tip. With finite element analysis the evaluation of the J-integral through the first equation above is relatively simple seeing as the numerics provides us with a stress and a strain field in our geometry. With the J-integral in hand we can then link our values of J to crack tip propagation rate with empirical data. We draw the analogy with linear theory where one uses Paris’s law to relate stress intensification factor to crack growth rate:

\[\frac{da}{dN} = C\cdot(\Delta K)^m\]

as we saw earlier J is in essence the nonlinear analog of K so we get:

\[ \frac{da}{dN} = C\cdot(\Delta J)^m\]

Here the fitting parameters C and m are material and geometry dependent.  Dowling (1976) for example gives the formula for fatigue crack growth rate in centre cracked A533B steel as [3]:

\[\frac{da}{dN} = 2.5\cdot 10^{-8}(\Delta J)^{1.56}\]

With this, one can determine the number of cycles dependent on the tolerable crack size in your system; thereby completing a nonlinear, plastic fatigue analysis of your system.

Figure 2: Equivalent von Mises Stress Distribution
Figure 3: Crack bottom fillet weld radially through shell

A recent application of this technique performed by Dynaflow for a client was in the fatigue analysis of a half pipe heat exchanger jacket with the finite element software Ansys. A 2D axis-symmetric model was built and under the design operating conditions high stresses and some plastic straining was observed between the coils as is showcased in figures 2 & 4. A number of crack locations were trialed two of which are shown in figures 3 & 5. Within Ansys one can define the crack tips and crack propagation directions and directly calculated the J-integral around these crack tips and thus, with an appropriate formula, calculate the number of cycles until failure if a crack were present at that location in the half pipe heat exchanger.

Figure 4: Plastic strain distribution
Figure 5: Crack top fillet weld through coil

For more information, please contact:

Edwin Schimmel
Project Engineer

Phone: +31 85 058 0046

References:

  • [1] James R Rice. A path independent integral and the approximate analysis of strain concentration by notches and cracks. Journal of Applied Mechanics, 1968.
  • [2] CQ Bowles. The j integral approach in elastic-plastic fatigue crack propagation. In Advances in Fatigue Science and Technology, pages 467–488. Springer, 1989.
  • [3] Norman E Dowling. Geometry effects and the j-integral approach to elastic-plastic fatigue crack growth. In Cracks and fracture. ASTM International, 1976.

    Note on author: Tammo Dukker is an engineering intern within our consultancy branch at Dynaflow Research Group. He is pursuing an MSc in Applied Physics at TU Delft with a focus on fluid mechanics. He has a keen interest in the application of fluid mechanics in the energy industry and offshore & coastal domains.

  • Related news

    3 June 2022
    DRG

    DRG welcomes David Gonzalez Pomares

    READ INFO
    2 June 2022

    FRP Pipe Stress Analysis Webinar – Use of the pipe stress envelope according to the design standard ISO 14692

    READ INFO
    13 May 2022

    Exhange Program Dynaflow and Stressman Engineering

    READ INFO
    • Dynaflow Dynaflow Research Group is ISO:9001 certified.
      ISO 9001:2015
      certified
    • HEAD OFFICE

      Laan van Oversteen 20
      6th floor
      2289 CX Rijswijk
      The Netherlands

    • MIDDLE EAST OFFICE

      Bay X2 Tower, Unit 3109
      Jumeirah Lakes Towers, cluster X
      P.O. Box 73001, Dubai
      United Arab Emirates

    • CONTACT

      Phone: +31 (0)85 058 0046
      E-mail: infoaanvraag@dynaflow.com

    • © Dynaflow Research Group BV

      • Terms & Conditions
      • Privacy Policy
      • Cookie policy (EU)
      • About us
      • Contact