Skip to main content

The J integral

Submitted by Zhigang Suo on

For a crack in an elastic body subject to a load, the elastic energy stored in the body is a function of two independent variables: the displacement of the load, and the area of the crack. The energy release rate is defined by the partial derivative of the elastic energy of the body with respect to the area of the crack.

This definition of the energy release rate assumes that the body is elastic, but invokes no field theory. Indeed, the energy release rate can be determined experimentally by measuring the load-displacement curves of identically loaded bodies with different areas of the cracks. No field need be measured.

Many materials, however, can be modeled with a field theory of elasticity. When a material is modeled by such a field theory, the energy release rate can be represented in terms of field variables by an integral, the J integral.

This lecture describes the J integral, along with examples of calculation. Uses of the J integral are often better appreciated in the context of individual applications, which we will describe in later lectures.

The J integral can be developed for both linear and nonlinear elastic theories. The nonlinear elastic theory will be used in class, and the linear elastic theory will be used in a homework problem.

These notes belong to a course on fracture mechanics

Attachment Size
J integral 2016 04 08.pdf 1.34 MB

Hi Zhigang,

Thanks a lot for this handy and useful lecture notes on J-integral.

In engineering applications, we often need to calculate stress intensity factors
for a corner crack. A convenient approach is to deduce them from J- or similar
integrals. Numerically, such integral does not appear to be path-independent.
And I read somewhere (unfortunately cannot recall where) that theoretically the
path-independence can only be guaranteed for a 2-D crack.

Do you have any comments on the 3-D J-integral?

 

Xiaobo

Sun, 03/21/2010 - 03:35 Permalink

Several points come to mind:

  1. If the material is inhomogenous only in the direction normal to the faces of the crack, the J integral is path-independent.  For example, the J intergal can be used to calculate the energy release rate of a crack on the interface between two materials.
  2. If the material is piece-wise homogeneous, when the tip of a crack is inside a homogeneous material, then the J integral along a contour inside the homogeneous material equals the energy release rate.
  3. If the material is piece-wise homogeneous, when the tip of the crack is on the interface, but the faces of the are off the interface, then the energy release rate is not well defined, and the J integral is path dependent.
  4. Nothing comes to mind of any use of the J integral when it is path-dependent. 
Fri, 04/02/2010 - 23:57 Permalink

Dear Reza:  You will get any useful result by using the J-integral around the point where two cracks intersect.  For example, here is a paper on a crack impinging upon an interface.  The J-integral is not useful here.  There are many examples like this.

R. Huang, J.H. Prévost, Z.Y. Huang, and Z. Suo, "Channel-cracking of thin films with the extended finite element method". Engineering Fracture Mechanics,70, 2513-2526 (2003).

Wed, 01/12/2011 - 23:33 Permalink

Dear zhigang

I have read the article (J.R. Rice, Mathematical analysis in the mechanics of fracture, Chapter 3)  that proves G=J for elastic materails. and I faced a problem,  why in two-dimensional crack we take A* any finite region in which the crack tip is imbedded?

can you please help me.

Sh.Eskandari

Wed, 05/12/2010 - 12:59 Permalink

Dear Kejie

Yes it is obvious, but when we are going to prove that energy releae rate is equivalent to J integral , there is an area integral that it should be taken over all the body but in Prof. Rice article for some reason that i didn't get it we take the area as any finite region containing the crack tip.

if you can just read this part of the article and help with this.

thank you

Shahin Eskandari 

 

Thu, 05/13/2010 - 12:05 Permalink

Hi Shahin Eskandari,

Sorry for the late reply, somehow I missed this thread. I have not read Prof. Rice's paper but I did see some proof in the textbooks.  My understanding so far is, since the integral is path-indepedent (i.e., it's a constant for any finite region containing the singular crack tip), we can take any area to prove G=J. 

This might be deviated from your question? 

Kejie

Sun, 05/16/2010 - 18:22 Permalink

Dear Kejie

Thank you for your reply.

I think we can not use path-independency, becouse in the proovement  we start from the rate of potential energy and then we show that it is equivalence to J-integra.

Anyway, can you tell me about the textbooks that you have seen the proovement.

Thanks.

Shahin Eskandari

Tue, 05/18/2010 - 06:46 Permalink

Dear Everyone

I am working on Fracture Mechanics, and facing problem to evaluate  J-Integral in 2D and 3D crack model and as well as in composite model by ANSYS.

If anyone have idea about it please send me tutorials or steps for the same problem.

Your help will be heartly appreciated.

Thanks.........

email ID: himanshu [at] iitp.ac.in

               hpiitp [at] gmail.com

Sun, 05/16/2010 - 08:50 Permalink

1.
Relationships between J-integral and the strain energy evaluated in a
finite volume surrounding the tip of sharp and blunt V-notches


International
Journal of Solids and Structures
, Volume 44, Issues 14-15, July
2007
, Pages 4621-4645

F. Berto, P. Lazzarin
2.
Elastic stress distributions for hyperbolic
and parabolic notches in round shafts under torsion and uniform
antiplane shear loadings


International Journal of
Solids and Structures
, Volume 45, Issues 18-19, September
2008
, Pages 4879-4901

M. Zappalorto, P. Lazzarin, J.R. Yates
3.
Some advantages derived from the use of the
strain energy density over a control volume in fatigue strength
assessments of welded joints


International Journal of
Fatigue
, Volume 30, Issue 8, August 2008, Pages
1345-1357


P. Lazzarin, F. Berto, F.J. Gomez, M. Zappalorto
4.
Fracture of U-notched specimens under mixed
mode: Experimental results and numerical predictions


Engineering
Fracture Mechanics
, Volume 76, Issue 2, January 2009,
Pages 236-249

F.J. Gómez, M. Elices, F. Berto, P. Lazzarin
5.
A review of the volume-based strain energy
density approach applied to V-notches and welded structures


Theoretical
and Applied Fracture Mechanics
, Volume 52, Issue 3, December
2009
, Pages 183-194

F. Berto, P. Lazzarin
6.
Local fatigue strength parameters for welded
joints based on strain energy density with inclusion of small-size
notches


Engineering Fracture Mechanics, Volume
76, Issue 8
, May 2009, Pages 1109-1130

D. Radaj, F.
Berto, P.
Lazzarin
7.
Effect of the thickness on elastic
deformation and quasi-brittle fracture of plate components


Engineering
Fracture Mechanics
, In Press, Corrected
Proof
, Available online 14 April 2010

A.
Kotousov, P.
Lazzarin, F. Berto, S. Harding
8.
Fatigue design of complex welded structures

International
Journal of Fatigue
, Volume 31, Issue 1, January 2009,
Pages 59-69

B. Atzori, P. Lazzarin, G. Meneghetti, M. Ricotta

 

Michele Ciavarella, Politecnico di BARI - Italy, Rector's delegate.
http://poliba.academia.edu/micheleciavarella
Editor, Italian Science Debate, www.sciencedebate.it
Associate Editor, Ferrari Millechili Journal, http://imechanica.org/node/7878

Sun, 05/16/2010 - 21:55 Permalink

soton.ac.uk
[PDF]
S Filippi, M Ciavarella, P Lazzarin -
International Journal of Fracture, 2002 - Springer

S.
FILIPPI1, M. CIAVARELLA2 and P. LAZZARIN3,∗ 1Department of Mechanical
Engineering

- University of Padova, Via Venezia 1, 35100 Padova (Italy); 2CNR-IRIS,
Computational Mechanics

of Solids, Str. Crocefisso 2B, 70125 Bari (Italy); 3Department of
Management and ...

 

Michele Ciavarella, Politecnico di BARI - Italy, Rector's delegate.
http://poliba.academia.edu/micheleciavarella
Editor, Italian Science Debate, www.sciencedebate.it
Associate Editor, Ferrari Millechili Journal, http://imechanica.org/node/7878

Sun, 05/16/2010 - 21:56 Permalink

8.

Duality and symmetry lost in solid mechanics

Comptes
Rendus Mécanique
, Volume 336, Issues 1-2, January-February
2008
, Pages 12-23

Huy Duong Bui

 Preview   Purchase PDF (193 K)
| Related
Articles

 

Michele Ciavarella, Politecnico di BARI - Italy, Rector's delegate.
http://poliba.academia.edu/micheleciavarella
Editor, Italian Science Debate, www.sciencedebate.it
Associate Editor, Ferrari Millechili Journal, http://imechanica.org/node/7878

Sun, 05/16/2010 - 21:58 Permalink

Dear all,

I am a new Abaqus user. My work involve with fracture modelling by using Abaqus.

For 2D Edge Crack problem, I have requested 5 values of J contour integral in abaqus, but i can not plot the J-integral path (J - integral versus Radius). I know abaqus will choose the path of each contour integral automatically. But How to show the path of each contour integral?

Thank you so much,

Sutham A.

Sat, 11/20/2010 - 03:13 Permalink

Hi All,

I have a question in mind for a while and unfortunately

have not yet found a clear answer to it. I would highly appreciate if you

provide me with your idea about this question:

 Is J-integral valid inside the process zone ahead of a

sharp crack? and if is valid, Is it equal to the value obtained from a contour

outside the process zone?

 Thanks

Ahmad Khayer D.

Wed, 01/05/2011 - 23:31 Permalink

To be precise the J-integral is only path independent for isothermal, non-linear elastic problems where there are no body forces. For monotonic loading where there is plasticity then it is path independent only under deformation plasticity. For in-elastic problems where there you hasve cyclic loading the J-integral is  not path dependent. In such situations the crack tip can be described by a related path integral called T*, that was first formulated by Satya Atluri. When J is path independent then J=T*. The importance of the J integral is related to its ability to characterise the near tip stresses and strains for a non-linear elastic material. The so called HRR solution. The experimental work of Professor Albert Kobayashi has shown that this is not true for cyclic loading.

 As an aside many people mistakenly think that the J-Integral was developed by Jim Rice. This is not true. It was first proposed in the early 50's by Professor Eshelby at the University of Sheffield. Prof Rice popularised it many years later. The 3D variant was first developed in a joint paper by Professors George Sih and George Irwin, both at Lehigh, shortly thereafter. You can chase down the relevant references with a simple Google search. The reference to the original Sih-Irwin paper is in one of Professor Sih's papers in his Journal, viz: Theoretical and Applied Fracture Mechanics.

 Hope this helps

 

Rhys  Jones

 

Thu, 01/06/2011 - 12:55 Permalink

Dear Colleagues,

My problem is that if I define a zero-area path around the crack tip of the mode-II ENF specimen and
calculate the traction vectors, there is a negative sign (2nd page)
that I do not understand and this is what leads to the wrong result. Can
You take a look at my note and tell me what is wrong?

page 1

page 2

Thanks

Andras

Sat, 10/15/2011 - 11:59 Permalink

Hi

i'm modeling crack with EFG method and using M-integral to calculate stress intensity and i have some quastions:

 is domain integral form of j-integral still path-independence?

what is the best 'q' function?

what domain size is good for the M-integral?

i found diffrence between answers with several domain size is it mean that m-integral is path-dependence?

and my bigest problem is sometimes i got worse answer with finer mesh!!! (also with uniform distribution)

and in some domain size  my enriched basis got worse answer than the ordinary basis!!!

is my model unstable or something like that???

what should i do???

i checked everything!!! my code works well what do you suggest????

 

 

 

Mon, 05/27/2013 - 06:52 Permalink

Dear Prof. Suo,

I have read your posted outlines for fracture meachanis and they are really helpful. Thank you for sharing them! But one problem is that I couldn't find all the outlines for all the listed topics. I could not find any pdf file for this one, for example. Are you still posting those outline?

Best,

Thu, 11/06/2014 - 05:12 Permalink

Hi Zhigang,

Thank you for the useful explanation. I would like to ask you if you can model a crack on a 3D curved surface like a pipe. How could this be modeled ? can i still use the contour integral method?

 

Thank you,

Diana.

Mon, 12/22/2014 - 21:46 Permalink

Dear Prof. Suo,

Rice introduced the J integral concept which is equal to the energy release rate for a uniform, linear or nonlinear elastic material free of body forces and subjected to a 2d deformation field (plane strain, plane stress etc.).

 

In ANSYS and Abaqus, the implementation is always referenced to the work of Shih (1986). There has been numerous publications on the use of J-integral to estimate Gc (critical energy release rate) for cracks which are located at the interface of dissimilar materials.

I could not find the relevant information if this approach is applicable, if yes, under what  circumstances? Could you please give me some information or references about the implementation of Jintegral at the interface of two materials?

I am looking forward to hearing from you,

Yalcin

 

Wed, 01/21/2015 - 07:17 Permalink