Two-Dimensional Mode I Crack Propagation in Saturated Ionized Porous Media Using Partition of Unity Finite Elements

Shales, clays, hydrogels, and tissues swell and shrink under changing osmotic conditions, which may lead to failure. The relationship between the presence of cracks and fluid flow has had little attention. The relationship between failure and osmotic conditions has had even less attention. The aim of this research is to study the effect of osmotic conditions on propagating discontinuities under different types of loads for saturated ionized porous media using the finite element method (FEM). Discontinuous functions are introduced in the shape functions of the FEM by partition of unity method, independently of the underlying mesh. Damage ahead of the crack-tip is introduced by a cohesive zone model. Tensile loading of a crack in an osmoelastic medium results in opening of the crack and high pressure gradients between the crack and the formation. The fluid flow in the crack is approximated by Couette flow. Results show that failure behavior depends highly on the load, permeability, (osmotic) prestress and the stiffness of the material. In some cases it is seen that when the crack propagation initiates, fluid is attracted to the crack-tip from the crack rather than from the surrounding medium causing the crack to close. The results show reasonable mesh-independent crack propagation for materials with a high stiffness. Stepwise crack propagation through the medium is seen due to consolidation, i.e., crack propagation alternates with pauses in which the fluid redistributes. This physical phenomenon challenges the numerical scheme. Furthermore, propagation is shown to depend on the osmotic prestressing of the medium. This mechanism may explain the tears observed in intervertebral disks as degeneration progresses.

[1]  D. S. Dugdale.  Yielding of steel sheets containing slits , 1960 .

[2]  G. I. Barenblatt.  THE MATHEMATICAL THEORY OF EQUILIBRIUM CRACKS IN BRITTLE FRACTURE , 1962 .

[3]  J. Rice,et al.  Some basic stress diffusion solutions for fluid‐saturated elastic porous media with compressible constituents , 1976 .

[4]  Pieter A. Vermeer,et al.  An accuracy condition for consolidation by finite elements , 1981 .

[5]  N. C. Huang,et al.  Hydraulic fracturing of a saturated porous medium-I: General theory , 1985 .

[6]  Robert L. Taylor,et al.  The patch test for mixed formulations , 1986 .

[7]  D. A. Spence,et al.  Transport of magma and hydrothermal solutions by laminar and turbulent fluid fracture , 1986 .

[8]  Y Lanir,et al.  Biorheology and fluid flux in swelling tissues. I. Bicomponent theory for small deformations, including concentration effects. , 1987, Biorheology.

[9]  O. Tervonen,et al.  Lumbar disc degeneration. Correlation between CT and CT/discography. , 1990, Acta radiologica.

[10]  Anthony R. Ingraffea,et al.  Simulation of hydraulic fracture propagation in poroelastic rock with application to stress measurement techniques , 1991 .

[11]  Emmanuel M Detournay,et al.  Plane strain analysis of a stationary hydraulic fracture in a poroelastic medium , 1991 .

[12]  Kenneth Runesson,et al.  Discontinuous Displacement Approximation for Capturing Plastic Localization , 1993 .

[13]  M. Ortiz,et al.  Computational modelling of impact damage in brittle materials , 1996 .

[14]  I. Babuska,et al.  The Partition of Unity Method , 1997 .

[15]  M. Battié,et al.  The influence of occupation on lumbar degeneration. , 1999, Spine.

[16]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .

[17]  Francisco Armero,et al.  An analysis of strong discontinuities in a saturated poro-plastic solid , 1999 .

[18]  Jonas Larsson,et al.  Localization Analysis of a Fluid Saturated Elastoplastic Porous Medium Using Regularized Discontinuities , 2000 .

[19]  W. Hutton,et al.  Does Long-Term Compressive Loading on the Intervertebral Disc Cause Degeneration? , 2000, Spine.

[20]  N. Broom,et al.  Biomechanical factors influencing nuclear disruption of the intervertebral disc. , 2001, Spine.

[21]  S. McGill,et al.  Intervertebral disc herniation: studies on a porcine model exposed to highly repetitive flexion/extension motion with compressive force. , 2001, Clinical biomechanics.

[22]  Garth N. Wells,et al.  Discontinuous modelling of strain localisation and failure , 2001 .

[23]  T. Belytschko,et al.  Extended finite element method for cohesive crack growth , 2002 .

[24]  Luciano Simoni,et al.  Cohesive fracture mechanics for a multi‐phase porous medium , 2003 .

[25]  A. Needleman,et al.  A cohesive segments method for the simulation of crack growth , 2003 .

[26]  Sally Roberts,et al.  Degeneration of the intervertebral disc , 2003, Arthritis research & therapy.

[27]  R. Lin,et al.  Geometric and morphological changes of the intervertebral disc under fatigue testing. , 2003, Clinical biomechanics.

[28]  Emmanuel M Detournay,et al.  The near-tip region of a fluid-driven fracture propagating in a permeable elastic solid , 2003, Journal of Fluid Mechanics.

[29]  M. Battié,et al.  Lumbar Disc Degeneration: Epidemiology and Genetic Influences , 2004, Spine.

[30]  N. Broom,et al.  Mechanically Induced Disruption of the Healthy Bovine Intervertebral Disc , 2004, Spine.

[31]  Michel Quintard,et al.  Upscaling Fractured Heterogeneous Media: Permeability and Mass Exchange Coefficient , 2006 .

[32]  Jjc Joris Remmers,et al.  Discontinuities in materials and structures: A unifying computational approach , 2006 .

[33]  Frank P T Baaijens,et al.  Influence of Osmotic Pressure Changes on the Opening of Existing Cracks in 2 Intervertebral Disc Models , 2006, Spine.

[34]  Franz-Josef Ulm,et al.  A micromechanical analysis of damage propagation in fluid-saturated cracked media , 2006 .

[35]  Julien Réthoré,et al.  A two‐scale approach for fluid flow in fractured porous media , 2006 .

[36]  M. Adams,et al.  What is Intervertebral Disc Degeneration, and What Causes It? , 2006, Spine.

[37]  Jan Carmeliet,et al.  A coupled discrete-continuum approach to simulate moisture effects on damage processes in porous materials , 2006 .

[38]  Bernhard A. Schrefler,et al.  On adaptive refinement techniques in multi-field problems including cohesive fracture , 2006 .

[39]  R. Al-Khoury,et al.  A computational model for fracturing porous media , 2007 .

[40]  Bernhard A. Schrefler,et al.  Mesh adaptation and transfer schemes for discrete fracture propagation in porous materials , 2007 .

[41]  A. Needleman,et al.  The simulation of dynamic crack propagation using the cohesive segments method , 2008 .

[42]  René de Borst,et al.  Challenges in computational materials science: Multiple scales, multi-physics and evolving discontinuities , 2008 .