Linear complementarity formulation for 3D frictional sliding problems

Frictional sliding on quasi-statically deforming faults and fractures can be modeled efficiently using a linear complementarity formulation. We review the formulation in two dimensions and expand the formulation to three-dimensional problems including problems of orthotropic friction. This formulation accurately reproduces analytical solutions to static Coulomb friction sliding problems. The formulation accounts for opening displacements that can occur near regions of non-planarity even under large confining pressures. Such problems are difficult to solve owing to the coupling of relative displacements and tractions; thus, many geomechanical problems tend to neglect these effects. Simple test cases highlight the importance of including friction and allowing for opening when solving quasi-static fault mechanics models. These results also underscore the importance of considering the effects of non-planarity in modeling processes associated with crustal faulting.

[1]  W. D. Stuart.  Strain softening prior to two‐dimensional strike slip earthquakes , 1979 .

[2]  Jeffrey C. Trinkle,et al.  Complementarity formulations and existence of solutions of dynamic multi-rigid-body contact problems with coulomb friction , 1996, Math. Program..

[3]  Maria Comninou,et al.  The angular dislocation in a half space , 1975 .

[4]  Paul Segall,et al.  Mechanics of discontinuous faults , 1980 .

[5]  M. Gringauz,et al.  On mathematical programming formulation for contact crack problems , 1993 .

[6]  S. Dirkse,et al.  The path solver: a nommonotone stabilization scheme for mixed complementarity problems , 1995 .

[7]  J.A.L. Napier,et al.  The spectral boundary element method: a new window on boundary elements in rock mechanics , 1992 .

[8]  Friedrich Pfeiffer,et al.  Mechanical System Dynamics , 2008 .

[9]  Michele L. Cooke,et al.  Erratum: ``Fracture localization along faults with spatially varying friction'' , 1997 .

[10]  Y. Fung.  A first course in continuum mechanics /2nd edition/ , 1977 .

[11]  Elizabeth Ritz,et al.  Closure of circular arc cracks under general loading: effects on stress intensity factors , 2011 .

[12]  J. Rice,et al.  Elastodynamic analysis for slow tectonic loading with spontaneous rupture episodes on faults with rate‐ and state‐dependent friction , 2000 .

[13]  F. Donath.  EXPERIMENTAL STUDY OF SHEAR FAILURE IN ANISOTROPIC ROCKS , 1961 .

[14]  David D. Pollard,et al.  On the patterns of wing cracks along an outcrop scale flaw: A numerical modeling approach using complementarity , 2008 .

[15]  Y. Kaneko,et al.  Supershear transition due to a free surface in 3-D simulations of spontaneous dynamic rupture on vertical strike-slip faults , 2010 .

[16]  J. Kaven,et al.  Growth of surface-breaching normal faults as a three-dimensional fracturing process , 2007 .

[17]  Georgios E. Stavroulakis,et al.  A complementarity problem formulation of the frictional grasping problem , 2000 .

[18]  Michele L. Cooke,et al.  Solving 3D boundary element problems using constrained iterative approach , 2010 .

[19]  W. F. Brace,et al.  Limits on lithospheric stress imposed by laboratory experiments , 1980 .

[20]  Michael C. Ferris,et al.  Compressional fractures considered as contact problems and mixed complementarity problems , 2000 .

[21]  F. Chester,et al.  Ultracataclasite structure and friction processes of the Punchbowl Fault , 1998 .

[22]  K. Komvopoulos,et al.  A Fractal Analysis of Stiction in Microelectromechanical Systems , 1997 .

[23]  M. Zoback.  State of stress and crustal deformation along weak transform faults , 1991, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[24]  D. J. Andrews,et al.  Bulletin of the Seismological Society of America , 1985 .

[25]  M. Popescu.  Boundary element methods in solid mechanics , 1985 .

[26]  Lei Zhang,et al.  Determining fault slip distribution of the Chi-Chi Taiwan earthquake with GPS and InSAR data using triangular dislocation elements , 2008 .

[27]  David D. Pollard,et al.  Integrating 3-D seismic data, field analogs, and mechanical models in the analysis of segmented normal faults in the Wytch Farm oil field, southern England, United Kingdom , 2001 .

[28]  Amir Sagy,et al.  Evolution of fault-surface roughness with slip , 2007 .

[29]  F. Maerten.  Geomechanics to solve geological structure issues: forward, inverse and restoration modeling , 2010 .

[30]  Leon M Keer,et al.  Modeling slip zones with triangular dislocation elements , 1992 .

[31]  J. Trinkle,et al.  On Dynamic Multi‐Rigid‐Body Contact Problems with Coulomb Friction , 1995 .

[32]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[33]  Gernot Beer,et al.  Programming the Boundary Element Method , 2001 .

[34]  W. Ritz.  Über eine neue Methode zur Lösung gewisser Variationsprobleme der mathematischen Physik. , 1909 .

[35]  Richard W. Cottle,et al.  Linear Complementarity Problem. , 1992 .

[36]  G. Dantzig,et al.  A generalization of the linear complementarity problem , 1970 .

[37]  Kelin Wang,et al.  Mechanics of low‐stress forearcs: Nankai and Cascadia , 1999 .

[38]  T. Wright,et al.  Integrating high‐precision aftershock locations and geodetic observations to model coseismic deformation associated with the 1995 Kozani‐Grevena earthquake, Greece , 2007 .

[39]  W. Griffith,et al.  Rough faults, distributed weakening, and off-fault deformation , 2009 .

[40]  J. Langbein,et al.  Slip on the San Andreas Fault at Parkfield, California, over Two Earthquake Cycles, and the Implications for Seismic Hazard , 2006 .

[41]  M. Zoback,et al.  How faulting keeps the crust strong , 2000 .

[42]  John W. Rudnicki,et al.  Fracture Mechanics Applied to the Earth's Crust , 1980 .

[43]  J. Weeks,et al.  Roughness and wear during brittle faulting , 1988 .

[44]  J. Rudnicki.  Physical models of earthquake instability and precursory processes , 1988 .

[45]  Evert Hoek,et al.  Fracture of Anisotropic Rock , 1964 .

[46]  Peter Wriggers,et al.  A new algorithm for numerical solution of 3D elastoplastic contact problems with orthotropic friction law , 2004 .

[47]  Y. Fung.  A first course in continuum mechanics , 1969 .

[48]  R. Harris.  Numerical Simulations of Large Earthquakes: Dynamic Rupture Propagation on Heterogeneous Faults , 2004 .

[49]  J. Napier,et al.  Symmetric‐Galerkin BEM simulation of fracture with frictional contact , 2003 .