Issue 61

H. Mazighi et alii, Frattura ed Integrità Strutturale, 61 (2022) 154-175; DOI: 10.3221/IGF-ESIS.61.11

(15)

  p n.u d = .(p u) d Ω -

p n.u d

f

f

f

N F Ω

N E Ω

E

The pressure load on the surface of the fracture is then formulated as a body force applied to the entire domain. To consider the crack propagation, we introduce the phase-field in Eqn. (15):

(16)

 p n.u d = g( ) .(p u) d Ω -  

 g( )p n.u d

.(p u) d Ω -

f

f

f

f

N E Ω

E

N

Thus, the external energy ψ S becomes:

  

s

 Ω ψ = f.u d Ω - g( ) .(p u) d Ω + g p n.u d + t.u d     N N f f Ω Ω Ω

(17)

By substituting Eqns. (5), (12), and (17) into Eqn. (1), the Fréchet derivative of the total potential energy ψ yields the following Euler equations [65]:                  2 2 + c 0 ε f f 0 G -l =2 1-k 1- ψ ε +p .u+u. p l (18)

and:

     f σ +p g +f=0

(19)

where the Cauchy stress tensor is given by:

 

 

 

  

  +

σ =g

λ tr ε I+2 με + λ tr ε I+2 μ

(20)

+

-

The regime flow inside the porous media material is based on Darcy’s law, which is derived from the Navier-Stokes equation. It describes a linear relationship between the velocity v (m/s) and the gradient of pressure p f (Pa). Lomiz [69] and Louis [70] carried out that the penetration of the fluid through a rock fracture follows the Cubic law developed using the parallel plates approach [71].

f p v 

 

(21)

k

with  the dynamic viscosity (Pa.s) and k (m 2 ) the permeability of the porous medium of the fluid respectively. Darcy’s law is only valid for low velocities when the regime is laminar [69,70]. We account for the irreversibility of the crack propagation process by adopting the following strain-history field H + (u, p f , t) [72]:           + + f ε f f x 0,t H u,p ,t =max ψ ε +p .u+u. p  (22)

Eqn. (18) can then be rewritten as:

         2 2 + c 0 0 G -l =2 1-k 1- H l

(23)

Thereby, the Euler-Lagrangian equation becomes:

158

Made with FlippingBook - Online Brochure Maker