TY - JOUR

T1 - A large deformation formulation for fluid flow in a progressively fracturing porous material

AU - Irzal, F.

AU - Remmers, J.J.C.

AU - Huyghe, J.M.R.J.

AU - Borst, de, R.

PY - 2013

Y1 - 2013

N2 - A general numerical model has been developed for fluid flow in a progressively fracturing porous medium subject to large deformations. The fluid flow away from the crack is modelled in a standard manner using Darcy’s relation. In the discontinuity a similar relation is assumed for the fluid flow, but with a different permeability to take into account the higher porosity within the crack due to progressive damage evolution. The crack is described in a discrete manner by exploiting the partition-of-unity property of finite element shape functions. The nucleation and the opening of micro-cracks are modelled by a traction-separation relation. A heuristic approach is adopted to model the orientation of the cracks at the interfaces in the deformed configuration. A two-field formulation is derived, with the solid and the fluid velocities as unknowns. The weak formulation is obtained, assuming a Total Lagrangian formulation. This naturally leads to a set of coupled equations for the continuous and for the discontinuous parts of the mixture. The resulting discrete equations are nonlinear due to the cohesive-crack model, the large-deformation kinematic relations, and the coupling terms between the fine scale and the coarse scale. The capabilities of the model are shown at the hand of some example problems.

AB - A general numerical model has been developed for fluid flow in a progressively fracturing porous medium subject to large deformations. The fluid flow away from the crack is modelled in a standard manner using Darcy’s relation. In the discontinuity a similar relation is assumed for the fluid flow, but with a different permeability to take into account the higher porosity within the crack due to progressive damage evolution. The crack is described in a discrete manner by exploiting the partition-of-unity property of finite element shape functions. The nucleation and the opening of micro-cracks are modelled by a traction-separation relation. A heuristic approach is adopted to model the orientation of the cracks at the interfaces in the deformed configuration. A two-field formulation is derived, with the solid and the fluid velocities as unknowns. The weak formulation is obtained, assuming a Total Lagrangian formulation. This naturally leads to a set of coupled equations for the continuous and for the discontinuous parts of the mixture. The resulting discrete equations are nonlinear due to the cohesive-crack model, the large-deformation kinematic relations, and the coupling terms between the fine scale and the coarse scale. The capabilities of the model are shown at the hand of some example problems.

U2 - 10.1016/j.cma.2012.12.011

DO - 10.1016/j.cma.2012.12.011

M3 - Article

VL - 256

SP - 29

EP - 37

JO - Computer Methods in Applied Mechanics and Engineering

JF - Computer Methods in Applied Mechanics and Engineering

SN - 0045-7825

ER -