SUMMARYWe present an incremental quasi-static contact algorithm for path-dependent frictional crack propagation in the framework of the extended finite element (FE) method. The discrete formulation allows for the modeling of frictional contact independent of the FE mesh. Standard Coulomb plasticity model is introduced to model the frictional contact on the surface of discontinuity. The contact constraint is borrowed from nonlinear contact mechanics and embedded within a localized element by penalty method. Newton-Raphson iteration with consistent linearization is used to advance the solution. We show the superior convergence performance of the proposed iterative method compared with a previously published algorithm called 'LATIN' for frictional crack propagation. Numerical examples include simulation of crack initiation and propagation in 2D plane strain with and without bulk plasticity. In the presence of bulk plasticity, the problem is also solved using an augmented Lagrangian procedure to demonstrate the efficacy and adequacy of the standard penalty solution.
We present a stabilized extended finite element formulation to simulate the hydraulic fracturing process in an elasto-plastic medium. The fracture propagation process is governed by a cohesive fracture model, where a trilinear traction-separation law is used to describe normal contact, cohesion and strength softening on the fracture face. Fluid flow inside the fracture channel is governed by the lubrication equation, and the flow rate is related to the fluid pressure gradient by the 'cubic' law. Fluid leak off happens only in the normal direction and is assumed to be governed by the Carter's leak-off model. We propose a 'local' U-P (displacement-pressure) formulation to discretize the fluid-solid coupled system, where volume shape functions are used to interpolate the fluid pressure field on the fracture face. The 'local' U-P approach is compatible with the extended finite element framework, and a separate mesh is not required to describe the fluid flow. The coupled system of equations is solved iteratively by the standard Newton-Raphson method. We identify instability issues associated with the fluid flow inside the fracture channel, and use the polynomial pressure projection method to reduce the pressure oscillations resulting from the instability. Numerical examples demonstrate that the proposed framework is effective in modeling 3D hydraulic fracture propagation.where the fluid pressure is uniformly equal to zero. Therefore, the traction on dry fracture face S d is purely governed by the traction-separation (cohesive) lawwhere t c represents the cohesive traction. We note that until now, we have not assumed any particular constitutive law for the bulk material. In other words, the previous equations are valid for both linear and nonlinear bulk material laws.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.