Modeling hydraulic cracks and inclusion interaction using XFEM

Published on Sep 1, 2018in Underground Space
· DOI :10.1016/J.UNDSP.2018.04.005
Bo He1
Estimated H-index: 1
(Leibniz University of Hanover),
Bo He6
Estimated H-index: 6
(Leibniz University of Hanover),
Xiaoying Zhuang55
Estimated H-index: 55
(Leibniz University of Hanover)
Sources
Abstract
Abstract In this study, we employ a coupled hydromechanical model to study the hydraulic fracture propagation path in porous media under the influence of existing pressurized voids. The hydraulic fracturing field study reveals that the existing natural voids and cracks alter the local properties of the porous media and influence the fracture propagation pattern. We incorporate these phenomena into the presented hydromechanical model, which is constructed from the mass and momentum balance equations for saturated porous media. The extended finite element method (XFEM) is applied for modeling the fluid flow through discrete cracks. The nonlinear hydromechanical equations are solved by the Newton–Raphson scheme with an implicit time integration procedure. Finally, numerical examples are presented and compared with experimental results. It is found that the fracture propagation path is significantly influenced by the existing pressurized voids and essential properties of the porous media; that is, the crack trends to propagate towards the pressurized voids.
📖 Papers frequently viewed together
1 Citations
2 Citations
206 Citations
References48
Newest
#1Huilong Ren (Bauhaus University, Weimar)H-Index: 11
#2Xiaoying Zhuang (Tongji University)H-Index: 55
Last. Timon Rabczuk (Ton Duc Thang University)H-Index: 101
view all 3 authors...
Abstract In this paper, we present a dual-horizon peridynamics formulation which allows for simulations with dual-horizon with minimal spurious wave reflection. We prove the general dual property for dual-horizon peridynamics, based on which the balance of momentum and angular momentum in PD are naturally satisfied. We also analyze the crack pattern of random point distribution and the multiple materials issue in peridynamics. For selected benchmark problems, we show that DH-PD is less sensitive...
297 CitationsSource
#1Huilong Ren (Bauhaus University, Weimar)H-Index: 11
#2Xiaoying Zhuang (Tongji University)H-Index: 55
Last. Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
view all 4 authors...
Summary In this paper, we develop a dual-horizon peridynamics (DH-PD) formulation that naturally includes varying horizon sizes and completely solves the ‘ghost force’ issue. Therefore, the concept of dual horizon is introduced to consider the unbalanced interactions between the particles with different horizon sizes. The present formulation fulfills both the balances of linear momentum and angular momentum exactly. Neither the ‘partial stress tensor’ nor the ‘slice’ technique is needed to ameli...
333 CitationsSource
#1Huilong Ren (Bauhaus University, Weimar)H-Index: 11
#2Xiaoying Zhuang (Leibniz University of Hanover)H-Index: 55
Last. Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
view all 3 authors...
We propose a new peridynamic formulation with shear deformation for linear elastic solid. The key idea lies in subtracting the rigid body rotation part from the total deformation. Based on the strain energy equivalence between classic local model and non-local model, the bond force vector is derived. A new damage rule of maximal deviatoric bond strain for elastic brittle fracture is proposed in order to account for both the tensile damage and shear damage. 2D and 3D numerical examples are tested...
35 CitationsSource
#1Pedro M. A. Areias (University of Évora)H-Index: 37
#2Mohammed A. Msekh (Bauhaus University, Weimar)H-Index: 12
Last. Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
view all 3 authors...
Abstract We propose a crack propagation algorithm which is independent of particular constitutive laws and specific element technology. It consists of a localization limiter in the form of the screened Poisson equation with local mesh refinement. This combination allows the capturing of strain localization with good resolution, even in the absence of a sufficiently fine initial mesh. In addition, crack paths are implicitly defined from the localized region, circumventing the need for a specific ...
198 CitationsSource
#1Pedro M. A. Areias (University of Évora)H-Index: 37
#2Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
Last. Pedro P. Camanho (University of Porto)H-Index: 70
view all 3 authors...
Abstract In the context of plane fracture problems, we introduce an algorithm based on our previously proposed rotation of edges but now including the injection of continuum softening elements directly in the process region. This is an extension of the classical smeared (or regularized) approach to fracture and can be seen as an intermediate proposition between purely cohesive formulations and the smeared modeling. Characteristic lengths in softening are explicitly included as width of injected ...
142 CitationsSource
#1F. Amiri (Bauhaus University, Weimar)H-Index: 4
#2Cosmin Anitescu (Bauhaus University, Weimar)H-Index: 15
Last. Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
view all 5 authors...
In this paper, we develop a method based on local maximum entropy shape functions together with enrichment functions used in partition of unity methods to discretize problems in linear elastic fracture mechanics. We obtain improved accuracy relative to the standard extended finite element method at a comparable computational cost. In addition, we keep the advantages of the LME shape functions, such as smoothness and non-negativity. We show numerically that optimal convergence (same as in FEM) fo...
146 CitationsSource
#1Pedro M. A. Areias (University of Évora)H-Index: 37
#2Daniel Dias-da-Costa (UC: University of Coimbra)H-Index: 21
Last. Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
view all 4 authors...
Within the theme of ductile fracture in metals, we propose an algorithm for FEM-based computational fracture based on edge rotations and smoothing of complementarity conditions. Rotation axes are the crack front nodes in surface discretizations and each rotated edge affects the position of only one or two nodes. Modified edge positions correspond to the predicted crack path. To represent softening, porous plasticity in the form of the Rousselier yield function is used. The finite strain integrat...
37 CitationsSource
#1Pedro M. A. Areias (University of Évora)H-Index: 37
#2Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
Last. Pedro P. Camanho (University of Porto)H-Index: 70
view all 3 authors...
We propose alternative methods for performing FE-based computational fracture: a mixed mode extrinsic cohesive law and crack evolution by edge rotations and nodal reposition. Extrinsic plastic cohesive laws combined with the discrete version of equilibrium form a nonlinear complementarity problem. The complementarity conditions are smoothed with the Chen-Mangasarian replacement functions which naturally turn the cohesive forces into Lagrange multipliers. Results can be made as close as desired t...
78 CitationsSource
#1Pedro M. A. Areias (University of Évora)H-Index: 37
#2Timon Rabczuk (KU: Korea University)H-Index: 101
Last. Daniel Dias-da-Costa (UC: University of Coimbra)H-Index: 21
view all 3 authors...
Abstract We propose an alternative, simpler algorithm for FEM-based computational fracture in brittle, quasi-brittle and ductile materials based on edge rotations. Rotation axes are the crack front edges (respectively nodes in surface discretizations) and each rotated edge affects the position of only one or two nodes. Modified positions of the entities minimize the difference between the predicted crack path (which depends on the specific propagation theory in use) and the edge or face orientat...
188 CitationsSource
#1Pedro M. A. Areias (University of Évora)H-Index: 37
#2Timon Rabczuk (Bauhaus University, Weimar)H-Index: 101
SUMMARY We propose a simple and efficient algorithm for FEM-based computational fracture of plates and shells with both brittle and ductile materials on the basis of edge rotation and load control. Rotation axes are the crack front nodes, and each crack front edge in surface discretizations affects the position of only one or two nodes. Modified positions of the entities maximize the modified mesh quality complying with the predicted crack path (which depends on the specific propagation theory i...
209 CitationsSource
Cited By3
Newest
#1Leilei Chen (Xinyang Normal University)
#1Leilei Chen (Xinyang Normal University)H-Index: 11
Last. Haojie Lian (TUT: Taiyuan University of Technology)H-Index: 13
view all 6 authors...
This paper simulates the propagation of hydraulic fractures under constant fluid pressure using the isogeometric analysis based on boundary element methods (BEM) i.e. the so called isogeomtric BEM. Due to the seamless integration between computer aided design (CAD) and numerical analysis, the isogeomtric BEM eliminates the meshing procedure and geometric errors, and explicitly parametrizes the crack propagation path. Graded knot insertion near the crack tip is used to effectively capture the str...
Source
#1Heng Zhang (CUMT: China University of Mining and Technology)H-Index: 5
#2Cai Ping Lu (CUMT: China University of Mining and Technology)H-Index: 3
Last. Hongyu Wang (UWA: University of Western Australia)H-Index: 6
view all 6 authors...
Abstract To uncover the failure and instability mechanism of coal-rock combination under loading, it is very crucial to investigate crack development and energy evolution. In this work, the uniaxial compressive and Brazilian split tests of coal-rock combinations were first performed. Secondly, the corresponding numerical models were built using the UDEC-Trigon method. Based on numerical results, the number, length and macroscopic area of cracks during failure and instability were determined, and...
1 CitationsSource
#1Yanjun Chang (Xida: Guangxi University)H-Index: 1
#2Qiao Song (Xida: Guangxi University)H-Index: 2
Last. Zhanguang Zheng (Xida: Guangxi University)H-Index: 1
view all 5 authors...
Based on the ABAQUS platform for finite element analysis, the extended finite element method (XFEM) considering elastoplastic constitutive relationship is developed, by which the displacement discontinuity across the crack surface and the strong nonlinearity near the crack tip can be described more accurately. The strip specimens with unilateral cracks and central cracks under uniaxial tension are simulated using the XFEM and the FEM, respectively. The J-integral across the crack on each increme...
2 CitationsSource