0
Research Papers

Homogenization and Path Independence of the J-Integral in Heterogeneous Materials OPEN ACCESS

[+] Author and Article Information
Chun-Jen Hsueh

Department of Mechanical and Civil Engineering,
California Institute of Technology,
Pasadena, CA 91125
e-mail: chhsueh@caltech.edu

Kaushik Bhattacharya

Department of Mechanical and Civil Engineering,
California Institute of Technology,
Pasadena, CA 91125
e-mail: bhatta@caltech.edu

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received July 13, 2016; final manuscript received July 20, 2016; published online August 22, 2016. Editor: Yonggang Huang.

J. Appl. Mech 83(10), 101012 (Aug 22, 2016) (3 pages) Paper No: JAM-16-1355; doi: 10.1115/1.4034294 History: Received July 13, 2016; Revised July 20, 2016

The J-integral that determines the driving force on a crack tip is a central concept of fracture mechanics. It is particularly useful since it is path independent in homogeneous materials. However, most materials are heterogeneous at a microscopic scale, and the J-integral is not necessarily path independent in heterogeneous media. In this paper, we prove the existence of an effective J-integral in heterogeneous media, show that it can be computed from the knowledge of the macroscopic or homogenized displacement fields and that it is path independent in macroscopically homogeneous media as long as the contours are large compared to the length scale of the heterogeneities. This result justifies the common engineering use of the J-integral.

FIGURES IN THIS ARTICLE
<>

It has been understood from the work of Griffith [1], Irwin [2], and others that the propagation of cracks is driven by the energy release rate, i.e., the rate of change of elastic energy with respect to crack extension. Rice [3] showed that this energy release rate may be described by the J-integral Display Formula

(1)J=Γti(Wδijukxiσkj)njds

where W is the stored energy density, u is the displacement, σ=W/F the stress, n is the normal to the contour Γ enclosing the crack tip, and t is the tangent to the crack at the crack tip. Conveniently, the J-integral is path independent, i.e., it does not depend on the contour Γ as long as the medium is homogeneous. This makes it extremely useful, since one can choose contours along which it is most convenient to evaluate the integrand.

However, the J-integral is not necessarily path independent in heterogeneous materials. And most materials are heterogeneous at a microscopic scale. Still, the J-integral has proven to be a most useful concept. Typically, one notes that the scale of heterogeneities is small compared to the engineering object of interest, uses homogenization theory to define an effective elastic medium which is homogeneous at the engineering scale, and applies elasticity theory and the J-integral to this effective homogenous medium.

Unfortunately, the relation between a “microscopic” and “macroscopic” J-integral remains open. Specifically, it is not clear whether the J-integral computed with the stress and the strain associated with the heterogeneous medium will converge to the J-integral computed with the stress and the strain associated with the effective homogeneous medium as the contour becomes very large. In short, it is not clear whether one can use the solutions to the homogenized equation to compute an effective J-integral.

Indeed, a casual examination of the expression for the J-integral in Eq. (1) suggests that the microscopic J-integral will in general be different from the macroscopic or homogenized J-integral. While we know from homogenization theory (specifically Hill's Lemma [4]) that the average of the microscopic stored energy density is equal to the macroscopic stored energy density, this does not appear to be true for the second term in the parenthesis in the integrand of Eq. (1). Specifically, the stress and the displacement gradient fluctuate at the microscopic scale in a heterogeneous medium. Therefore, it is generally not true that the product of their averages is equal to the average of their products; in other words, it is unclear if Display Formula

(2)ukxiσkj=ukxiσkj?

This raises an issue of using the J-integral in engineering practice.

Further, the recent decades have seen an attempt to use heterogeneities to enhance fracture toughness. Furthermore, nature exploits microstructure to enhance toughness of nacre and other shells. Finally, the emergence of 3D printing and other methods of additive material synthesis opens the possibility of exploiting carefully controlled heterogeneity for enhancing the toughness of materials. All of this has led to a new interest in understanding the effective toughness of heterogeneous media (see for example, Ref. [5] and references therein).

In this paper, we use homogenization theory in a quasi-periodic setting to show the existence of a macroscopic J-integral, and prove that this is path independent in a macroscopically homogeneous material if the path is large compared to the size of the heterogeneities. The path independence of the J-integral follows from the fact that it is closely related to the configurational stress tensor or the Eshelby energy momentum tensor Display Formula

(3)Cij=Wδijukxiσkj

that satisfies Display Formula

(4)Cijxj=*W*xi

where */*xi denotes the explicit derivative with respect to xi. This equation may be obtained from the equilibrium equation, and is also referred to as the configurational force balance. Integrating this equation over an annular region between two contours and using the divergence theorem leads to the path independence of the J-integral in homogeneous materials. We show that the homogenization of this equation retains the same form leading to a homogenized configurational stress tensor and homogenized J-integral.

We recall homogenization of the variational formulation of elasticity in Sec. 2 and derive our main result in Sec. 3. We conclude in Sec. 4.

Consider a domain Ω (a bounded open set in N) with a heterogeneous elastic medium where the heterogeneities have a length scale ε1=diam(Ω) in dimensionless units. The domain may contain a smooth crack of whose length is O(1). Specifically, we assume that the medium is quasi-periodic so that the stored energy density Display Formula

(5)W̃ε(F,x)=W(F,x,xε)

where F is the displacement gradient, W:N×N×Ω×Y, Y is the unit cube in N and W(F,x,y) is periodic in y for each x, F. This is shown schematically in Fig. 1: the medium appears periodic if we look closely at some point xΩ.

In the case of linear elasticity Display Formula

(6)W(F,x,y)=12Fijijkl(x,y)Fkl

where :Ω×YN4 is the elastic modulus satisfying major and minor symmetries, and (x,y) is periodic in y for each x.

We seek solutions of the equations of elasticity by seeking to minimize the total energy Display Formula

(7)Eε[u]=ΩW̃ε(u,x)dxL[u]=ΩW(u,x,xε)dxL[u]

where L depends on the body force and boundary tractions among all displacements u:Ω3 that satisfy the imposed displacement boundary conditions. This problem is difficult because W and consequently the solution oscillates on the scale of the heterogeneities ε. Homogenization theory [6,7] states that if ε is small enough and if W(F,x,y) is convex in F for each x, y,1 then we can replace the problem above with the following effective problem: minimize Display Formula

(8)E[u]=ΩW¯(u,x)dxL[u]

among all displacements that satisfy the imposed displacement boundary conditions where W¯:N×N×Ω is the effective elastic energy density and may be obtained by solving the following problem for each F, x: minimize Display Formula

(9)W¯(F,x)=YW(F+φ(y),x,y)dy

overall periodic displacement fields φ:YN. Note that the integrand of Eq. (8) is smooth on the scale of ε and thus the solution is also expected to be smooth at that scale. It is also true that the effective stress is given by Display Formula

(10)σ¯ij=W¯Fij

Further, under suitable growth and strict convexity conditions on W, the minimum is attained and unique up to an inessential translation. We call the minimum φF,x. So Display Formula

(11)E[u]=ΩYW(xu+yφxu,x,x,y)dydxL[u]

We use the effective functional (11) to derive an effective configurational force balance. We take L=0 for convenience (i.e., no body force, and traction-free and displacement boundary conditions), but the treatment is easily modified otherwise.

Let u¯ minimize E defined in Eq. (11) for some given boundary conditions. Recall that it is smooth on the scale of ε. We now consider a variation, but by rearranging the domain.2 Consider a family of rearrangments zs:ΩΩ one to one and onto for s[0,1] that satisfy z = x on Ω,detxz>0xΩ, and z0(x)=x. Set u¯s(x)=u¯(zs(x)). Note that u¯s is a family of perturbations of the minimizer u¯ with u¯0=u¯. Therefore, the function Display Formula

(12)f(s):=E[u¯(zs)]

has a minimum at s = 0 and therefore Display Formula

(13)f(0)=0

We now compute f(0).

Set ψs=φxu¯s,x and F¯(x)=xu¯(x), and note that Display Formula

(14)xu¯s=x(u¯(zs(x)))=F¯(zs(x))xzs(x)

So Display Formula

(15)f(s)=ΩYW(F¯(zs(x))xzs+yψs,x,y)dydx
Display Formula
(16)=ΩYW(F¯(z)(zxs)1+yψs,xs,y)Jdydz

where we have changed integration variables from x to z by inverting zs to obtain x=xs(z), and set J=det(zxs). Now, set Display Formula

(17)F¯s=F¯(z)(zxs)1,Gs=F¯s+yψs=F¯s+yφF¯s,xs

We can now calculate Display Formula

(18)f(s)=ΩY(WFijGijs˙+*W*xixi˙+WJ˙J)Jdydz

where we use */*xi to represent the explicit derivative with respect to xi and (q)˙=dq/ds to denote the total derivative of q with respect to s. Recalling the identities Display Formula

(19)A1¯˙=A1A˙A1,detA¯˙=(detA)AijTA˙ij

we obtain Display Formula

(20)G˙ijs=F¯˙ijs+ψ˙isyj=F¯ik(z)(xkszl)1x˙lszm(xmszj)1+ψ˙isyj
Display Formula
(21)J˙J=(xiszj)T(x˙iszj)1

Further, from the unit cell problem of minimizing (9), we can infer that Display Formula

(22)YWFijψ˙isyjdy=0

We substitute Eqs. (20)(22) into Eq. (18), set s = 0 (so that xis/zj=δij) and change integration variables back to xDisplay Formula

(23)f(0)=ΩY((WδikWFijF¯ik)x˙kxj+*W*xkx˙k)dydx

where x˙=x˙s(z)|s=0. We note that F¯ and x˙ are independent of y and therefore we can integrate with respect to y. We obtain by recalling (9)Display Formula

(24)f(0)=Ω((W¯δikW¯FijF¯ik)x˙kxj+*W¯*xkx˙k)dx
Display Formula
(25)=Ω(xj(W¯δikW¯FijF¯ik)+*W¯*xkx˙k)x˙kdx

where we have used the divergence theorem and the fact that x˙=0 on Ω (since x=zs(x) on Ω). Now, since f(0)=0 for all zs and thus arbitrary x˙, we obtain the macroscopic configurational force balanceDisplay Formula

(26)C¯ijxj=*W¯*xi

where Display Formula

(27)C¯ij=W¯δijF¯kiW¯Fkj=W¯δijF¯kiσ¯kj

is the effective configurational stress tensor.

Finally, if the material is macroscopically homogeneous, i.e., W and consequently W¯ is independent of x, then we see from Eq. (26) that Display Formula

(28)0=DC¯ijxjdx=DC¯ijnjdA

using the divergence theorem where D is any domain that is large compared to the size of the heterogeneities (ε). Now, consider a domain with a crack in two dimensions. Given any two contours Γ1 and Γ2 that contain the crack tip, set D to be the annular region between the contours. Now, D=Γ1Γ2 with the outward normal −n and n on the two segments of the boundary. Taking the inner product of the equation above with the tangent to the crack tip, we obtain Display Formula

(29)Γ1ti(W¯δiju¯kxiσ¯kj)njds=Γ2ti(W¯δiju¯kxiσ¯kj)njds

or the path independence of the macroscopic J-integral.

We have shown the existence of a homogenized configuration stress tensor and configurational force balance in a quasi-periodic medium. We have used these to show the existence of an effective J-integral that may be evaluated from the macroscopic displacements and stresses, and that is path independent in macroscopically homogeneous media as long as the contours are large enough.

We gratefully acknowledge the financial support of the U.S. National Science Foundation Award No. DMS-1535083 under the Designing Materials to Revolutionize and Engineer our Future (DMREF) Program.

Griffith, A. , 1921, “ The Phenomena of Rupture and Flow in Solids,” Philos. Trans. R. Soc. A, 221(582–593), pp. 163–198. [CrossRef]
Irwin, G. , 1960, “ Fracture Mechanics,” Proceedings of the First Symposium on Naval Structural Mechanics, J. Goodier and N. Hoff, eds., Pergamon Press, New York, pp. 557–594.
Rice, J. , 1968, “ A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks,” ASME J. Appl. Mech., 35(2), pp. 379–386. [CrossRef]
Hill, R. , 1963, “ Elastic Properties of Reinforced Solids: Some Theoretical Principles,” J. Mech. Phys. Solids, 11(5), pp. 357–372. [CrossRef]
Hossain, M. , Hsueh, C.-J. , Bourdin, B. , and Bhattacharya, K. , 2014, “ Effective Toughness of Heterogeneous Materials,” J. Mech. Phys. Solids, 71, pp. 15–32. [CrossRef]
Braides, A. , 1985, “ Homogenization of Some Almost Periodic Coercive Functional,” Rendiconti. Accademia Nazionale delle Scienze detta dei XL, Vol. 9 (Serie V. Memorie di Matematica e Applicazioni), L'Accademia, Rome, Italy, pp. 313–321.
Müller, S. , 1987, “ Homogenization of Nonconvex Integral Functionals and Cellular Elastic Materials,” Arch. Ration. Mech. Anal., 99(3), pp. 189–212. [CrossRef]
Copyright © 2016 by ASME
View article in PDF format.

References

Griffith, A. , 1921, “ The Phenomena of Rupture and Flow in Solids,” Philos. Trans. R. Soc. A, 221(582–593), pp. 163–198. [CrossRef]
Irwin, G. , 1960, “ Fracture Mechanics,” Proceedings of the First Symposium on Naval Structural Mechanics, J. Goodier and N. Hoff, eds., Pergamon Press, New York, pp. 557–594.
Rice, J. , 1968, “ A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks,” ASME J. Appl. Mech., 35(2), pp. 379–386. [CrossRef]
Hill, R. , 1963, “ Elastic Properties of Reinforced Solids: Some Theoretical Principles,” J. Mech. Phys. Solids, 11(5), pp. 357–372. [CrossRef]
Hossain, M. , Hsueh, C.-J. , Bourdin, B. , and Bhattacharya, K. , 2014, “ Effective Toughness of Heterogeneous Materials,” J. Mech. Phys. Solids, 71, pp. 15–32. [CrossRef]
Braides, A. , 1985, “ Homogenization of Some Almost Periodic Coercive Functional,” Rendiconti. Accademia Nazionale delle Scienze detta dei XL, Vol. 9 (Serie V. Memorie di Matematica e Applicazioni), L'Accademia, Rome, Italy, pp. 313–321.
Müller, S. , 1987, “ Homogenization of Nonconvex Integral Functionals and Cellular Elastic Materials,” Arch. Ration. Mech. Anal., 99(3), pp. 189–212. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

A macroscopic crack in a quasi-periodic heterogeneous medium

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In