Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA phase-field model of frictional shear fracture in geologic materials
94
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Geologic shear fractures such as faults and slip surfaces involve marked friction along the discontinuities as they are subjected to significant confining pressures. This friction plays a critical role in the growth of these shear fractures, as revealed by the fracture mechanics theory of Palmer and Rice decades ago. In this paper, we develop a novel phase-field model of shear fracture in pressure-sensitive geomaterials, honoring the role of friction in the fracture propagation mechanism. Building on a recently proposed phase-field method for frictional interfaces, we formulate a set of governing equations for different contact conditions (or lack thereof) in which frictional energy dissipation emerges in the crack driving force during slip. We then derive the degradation function and the threshold fracture energy of the phase-field model such that the stress-strain behavior is insensitive to the length parameter for phase-field regularization. This derivation procedure extends a methodology used in recent phase-field models of cohesive tensile fracture to shear fracture in frictional materials in which peak and residual strengths coexist and evolve by confining pressure. The resulting phase-field formulation is demonstrably consistent with the theory of Palmer and Rice. Numerical examples showcase that the proposed phase-field model is a physically sound and numerically efficient method for simulating shear fracture processes in geologic materials, such as faulting and slip surface growth.
rate research
Read More
Cracking of rocks and rock-like materials exhibits a rich variety of patterns where tensile (mode I) and shear (mode II) fractures are often interwoven. These mixed-mode fractures are usually cohesive (quasi-brittle) and frictional. Although phase-field modeling is increasingly used for rock fracture simulation, no phase-field formulation is available for cohesive and frictional mixed-mode fracture. To address this shortfall, here we develop a double-phase-field formulation that employs two different phase fields to describe cohesive tensile fracture and frictional shear fracture individually. The formulation rigorously combines the two phase fields through three approaches: (i) crack-direction-based decomposition of the strain energy into the tensile, shear, and pure compression parts, (ii) contact-dependent calculation of the potential energy, and (iii) energy-based determination of the dominant fracturing mode in each contact condition. We validate the proposed model, both qualitatively and quantitatively, with experimental data on mixed-mode fracture in rocks. The validation results demonstrate that the double-phase-field model -- a combination of two quasi-brittle phase-field models -- allows one to directly use material strengths measured from experiments, unlike brittle phase-field models for mixed-mode fracture in rocks. Another standout feature of the double-phase-field model is that it can simulate, and naturally distinguish between, tensile and shear fractures without complex algorithms.
The analysis of surface wave dispersion curves is a way to infer the vertical distribution of shear-wave velocity. The range of applicability is extremely wide going, for example, from seismological studies to geotechnical characterizations and exploration geophysics. However, the inversion of the dispersion curves is severely ill-posed and only limited efforts have been put into the development of effective regularization strategies. In particular, relatively simple smoothing regularization terms are commonly used, even when this is in contrast with the expected features of the investigated targets. To tackle this problem, stochastic approaches can be utilized, but they are too computationally expensive to be practical, at least, in the case of large surveys. Instead, within a deterministic framework, we evaluate the applicability of a regularizer capable of providing reconstructions characterized by tunable levels of sparsity. This adjustable stabilizer is based on the minimum support regularization, applied before on other kinds of geophysical measurements, but never on surface wave data. We demonstrate the effectiveness of this stabilizer on i) two benchmark - publicly available - datasets at crustal and near-surface scales, ii) an experimental dataset collected on a well-characterized site. In addition, we discuss a possible strategy for the estimation of the depth of investigation. This strategy relies on the integrated sensitivity kernel used for the inversion and calculated for each individual propagation mode. Moreover, we discuss the reliability, and possible caveats, of the direct interpretation of this particular estimation of the depth of investigation, especially in the presence of sharp boundary reconstructions.
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 methods capability for modeling crack growth with frictional contact.
We present a stochastic modeling framework for atomistic propagation of a Mode I surface crack, with atoms interacting according to the Lennard-Jones interatomic potential at zero temperature. Specifically, we invoke the Cauchy-Born rule and the maximum entropy principle to infer probability distributions for the parameters of the interatomic potential. We then study how uncertainties in the parameters propagate to the quantities of interest relevant to crack propagation, namely, the critical stress intensity factor and the lattice trapping range. For our numerical investigation, we rely on an automated version of the so-called numerical-continuation enhanced flexible boundary (NCFlex) algorithm.
We consider a phase-field fracture propagation model, which consists of two (nonlinear) coupled partial differential equations. The first equation describes the displacement evolution, and the second is a smoothed indicator variable, describing the crack position. We propose an iterative scheme, the so-called $L$-scheme, with a dynamic update of the stabilization parameters during the iterations. Our algorithmic improvements are substantiated with two numerical tests. The dynamic adjustments of the stabilization parameters lead to a significant reduction of iteration numbers in comparison to constant stabilization values.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
