A numerical contact algorithm in saturated porous media with the extended finite element method

In this paper, a coupled hydro-mechanical formulation is developed for deformable porous media subjected to crack interfaces in the framework of extended finite element method. Governing equations of the porous medium consist of the momentum balance of the bulk together with the momentum balance and continuity equations of the fluid phase, known as [inline-graphic not available: see fulltext] formulation. The discontinuity in fractured porous medium is modeled for both opening and closing modes that results in the fluid flow within the fracture, and/or contact behavior at the crack edges. The fluid flow through the fracture is assumed to be viscous and is modeled by employing the Darcy law in which the permeability of fracture is obtained using the cubic law. The contact condition in fractured porous medium is handled by taking the advantage from two different algorithms of LATIN method and penalty algorithm. The effect of contact on fluid phase is employed by considering no leak-off from/into the porous medium. The nonlinearity of coupled equations produced due to opening and closing modes is carried out using an iterative algorithm in the Newton–Raphson procedure. Finally, several numerical examples are solved to illustrate the performance of proposed X-FEM method for hydro-mechanical behavior of fractured porous media with opening and closing modes.