A phase‐field method for modeling cracks with frictional contact

We introduce a phase-field method for continuous modeling of cracks with frictional contacts. Compared with standard discrete methods for frictional contacts, the phase-field method has two attractive features: (1) it can represent arbitrary crack geometry without an explicit function or basis enrichment, and (2) it does not require an algorithm for imposing contact constraints. The first feature, which is common in phase-field models of fracture, is attained by regularizing a sharp interface geometry using a surface density functional. The second feature, which is a unique advantage for contact problems, is achieved by a new approach that calculates the stress tensor in the regularized interface region depending on the contact condition of the interface. Particularly, under a slip condition, this approach updates stress components in the slip direction using a standard contact constitutive law, while making other stress components compatible with stress in the bulk region to ensure non-penetrating deformation in other directions. We verify the proposed phase-field method using stationary interface problems simulated by discrete methods in the literature. Subsequently, by allowing the phase field to evolve according to brittle fracture theory, we demonstrate the proposed method's capability for modeling crack growth with frictional contact.

[1]  John E. Dolbow,et al.  A Nitsche stabilized finite element method for frictional sliding on embedded interfaces. Part II: Intersecting interfaces , 2013 .

[2]  Ronaldo I. Borja,et al.  Mechanical aspects of thrust faulting driven by far-field compression and their implications for fold geometry , 2007 .

[3]  Ronaldo I. Borja,et al.  Extended finite element framework for fault rupture dynamics including bulk plasticity , 2013 .

[4]  A. Combescure,et al.  A mixed augmented Lagrangian‐extended finite element method for modelling elastic–plastic fatigue crack growth with unilateral contact , 2007 .

[5]  L. Wong,et al.  Crack Coalescence in Molded Gypsum and Carrara Marble: Part 1. Macroscopic Observations and Interpretation , 2009 .

[6]  S. Nemat-Nasser,et al.  Compression‐induced nonplanar crack extension with application to splitting, exfoliation, and rockburst , 1982 .

[7]  Ronaldo I. Borja,et al.  Finite deformation formulation for embedded frictional crack with the extended finite element method , 2010 .

[8]  R. Borja,et al.  Pore‐scale modeling of deformation and shear band bifurcation in porous crystalline rocks , 2016 .

[9]  John E. Dolbow,et al.  A phase-field formulation for dynamic cohesive fracture , 2018, Computer Methods in Applied Mechanics and Engineering.

[10]  Ethan T. Coon,et al.  A Nitsche-extended finite element method for earthquake rupture on complex fault systems , 2011 .

[11]  R. Juanes,et al.  Phase field model of fluid‐driven fracture in elastic media: Immersed‐fracture formulation and validation with analytical solutions , 2017 .

[12]  Richard A. Regueiro,et al.  Embedded strong discontinuity finite elements for fractured geomaterials with variable friction , 2007 .

[13]  Danna Zhou,et al.  d. , 1934, Microbial pathogenesis.

[14]  Ted Belytschko,et al.  An extended finite element method for modeling crack growth with frictional contact , 2001 .

[15]  Amir R. Khoei,et al.  An enriched finite element algorithm for numerical computation of contact friction problems , 2007 .

[16]  David Wells,et al.  The deal.II library, Version 9.0 , 2018, J. Num. Math..

[17]  J. Oden,et al.  Algorithms and numerical results for finite element approximations of contact problems with non-classical friction laws , 1984 .

[18]  Mervyn Paterson,et al.  Experimental rock Deformation: the Brittle Field (2nd ed) , 2005 .

[19]  J. C. Simo,et al.  An augmented lagrangian treatment of contact problems involving friction , 1992 .

[20]  P. Wriggers,et al.  Crack face contact for a hexahedral-based XFEM formulation , 2012 .

[21]  Christian Miehe,et al.  Thermodynamically consistent phase‐field models of fracture: Variational principles and multi‐field FE implementations , 2010 .

[22]  Christian Miehe,et al.  A phase field model for rate-independent crack propagation: Robust algorithmic implementation based on operator splits , 2010 .

[23]  J. Choo,et al.  Liquid CO2 Fracturing: Effect of Fluid Permeation on the Breakdown Pressure and Cracking Behavior , 2018, Rock Mechanics and Rock Engineering.

[24]  R. Borja,et al.  On the pore‐scale mechanisms leading to brittle and ductile deformation behavior of crystalline rocks , 2015 .

[25]  Michael Kaliske,et al.  A phase-field crack model based on directional stress decomposition , 2018, Computational Mechanics.

[26]  J. C. Simo,et al.  An analysis of strong discontinuities induced by strain-softening in rate-independent inelastic solids , 1993 .

[27]  Jean-Jacques Marigo,et al.  Morphogenesis and propagation of complex cracks induced by thermal shocks , 2013 .

[28]  John E. Dolbow,et al.  A robust Nitsche’s formulation for interface problems , 2012 .

[29]  Christian Miehe,et al.  Phase field modeling of fracture in multi-physics problems. Part I. Balance of crack surface and failure criteria for brittle crack propagation in thermo-elastic solids , 2015 .

[30]  Ronaldo I. Borja,et al.  Stabilized low-order finite elements for frictional contact with the extended finite element method , 2010 .

[31]  Daniel Marceau,et al.  The adapted augmented Lagrangian method: a new method for the resolution of the mechanical frictional contact problem , 2012 .

[32]  J. C. Simo,et al.  A perturbed Lagrangian formulation for the finite element solution of contact problems , 1985 .

[33]  WaiChing Sun,et al.  Computational thermomechanics of crystalline rock, Part I: A combined multi-phase-field/crystal plasticity approach for single crystal simulations , 2018, Computer Methods in Applied Mechanics and Engineering.

[34]  Mary F. Wheeler,et al.  Pressure and fluid-driven fracture propagation in porous media using an adaptive finite element phase field model , 2016 .

[35]  Luca Heltai,et al.  The deal.II library, Version 9.2 , 2020, J. Num. Math..

[36]  Ronaldo I. Borja,et al.  A finite element model for strain localization analysis of strongly discontinuous fields based on standard Galerkin approximation , 2000 .

[37]  John E. Dolbow,et al.  On strategies for enforcing interfacial constraints and evaluating jump conditions with the extended finite element method , 2004 .

[38]  Christian Miehe,et al.  A phase‐field model for chemo‐mechanical induced fracture in lithium‐ion battery electrode particles , 2016 .

[39]  H. H. Einstein,et al.  Crack Coalescence in Molded Gypsum and Carrara Marble: Part 2—Microscopic Observations and Interpretation , 2008 .

[40]  A. Ingraffea,et al.  FINITE ELEMENT MODELS FOR ROCK FRACTURE MECHANICS , 1980 .

[41]  Christian Linder,et al.  A variational framework to model diffusion induced large plastic deformation and phase field fracture during initial two-phase lithiation of silicon electrodes , 2016 .

[42]  Sia Nemat-Nasser,et al.  Compression‐induced microcrack growth in brittle solids: Axial splitting and shear failure , 1985 .

[43]  Ronaldo I. Borja,et al.  An extended finite element framework for slow‐rate frictional faulting with bulk plasticity and variable friction , 2009 .

[44]  R. Borja Assumed enhanced strain and the extended finite element methods: A unification of concepts , 2008 .

[45]  T. Laursen Computational Contact and Impact Mechanics , 2003 .

[46]  R. Taylor,et al.  Lagrange constraints for transient finite element surface contact , 1991 .

[47]  WaiChing Sun,et al.  Coupled phase-field and plasticity modeling of geological materials: From brittle fracture to ductile flow , 2018 .

[48]  Phillipp Bergmann,et al.  Fundamentals Of Structural Geology , 2016 .

[49]  Z. T. Bieniawski,et al.  Brittle fracture propagation in rock under compression , 1965 .

[50]  Cv Clemens Verhoosel,et al.  A phase-field description of dynamic brittle fracture , 2012 .

[51]  P. Wriggers,et al.  FINITE ELEMENT FORMULATION OF LARGE DEFORMATION IMPACT-CONTACT PROBLEMS WITH FRICTION , 1990 .

[52]  D. Owen,et al.  Computational model for 3‐D contact problems with friction based on the penalty method , 1992 .

[53]  Peter Wriggers,et al.  Computational Contact Mechanics , 2002 .

[54]  Li Shen,et al.  Identifying the crack path for the phase field approach to fracture with non-maximum suppression , 2016 .

[55]  Jean-Jacques Marigo,et al.  Regularized formulation of the variational brittle fracture with unilateral contact: Numerical experiments , 2009 .

[56]  Jean H. Prevost,et al.  Frictional slip plane growth by localization detection and the extended finite element method (XFEM) , 2011 .

[57]  Ronaldo I. Borja,et al.  A contact algorithm for frictional crack propagation with the extended finite element method , 2008 .

[58]  WaiChing Sun,et al.  Cracking and damage from crystallization in pores: Coupled chemo-hydro-mechanics and phase-field modeling , 2018, Computer Methods in Applied Mechanics and Engineering.

[59]  Vinh Phu Nguyen,et al.  Length scale and mesh bias sensitivity of phase-field models for brittle and cohesive fracture , 2019, Engineering Fracture Mechanics.

[60]  T. Wong,et al.  Brittle faulting and cataclastic flow in Tavel limestone , 2010 .

[61]  Isaac Harari,et al.  An efficient finite element method for embedded interface problems , 2009 .

[62]  W. Bangerth,et al.  deal.II—A general-purpose object-oriented finite element library , 2007, TOMS.

[63]  Joshua A. White Anisotropic damage of rock joints during cyclic loading: constitutive framework and numerical integration , 2014 .