Discrete element model for the analysis of fluid-saturated fractured poro-plastic medium based on sharp crack representation with embedded strong discontinuities

• Published in: Computer methods in applied mechanics and engineering Vol. 298; pp. 407 - 427
• Main Authors: Nikolic, Mijo, Ibrahimbegovic, Adnan, Miscevic, Predrag
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01-01-2016
• Subjects: Discrete element method; Embedded discontinuity; Heterogeneous materials; Localized failure; Poroplasticity; Poroplasticity; Heterogeneous materials; Embedded discontinuity; Discrete element method; Localized failure
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2015.10.009

In this work, we present a discrete beam lattice model capable of simulating localized failure in a heterogeneous fluid-saturated poro-plastic solid. Coupling conditions between the solid and fluid phase are governed by the Biot’s porous media theory, enhanced with the fluid flow through cracks. The basis for development of discrete 2D plane strain model representation of heterogeneous material consisting of material grains, is an assembly of Voronoi cells that are kept together by cohesive links in terms of Timoshenko beams. Localized failure of saturated medium is enabled through embedded discontinuities positioned in cohesive links where Timoshenko beam’s longitudinal and transversal directions possess enhanced kinematics representing failure modes I and II. The model can also take into account the fracture process zone with pre-existing microcracks coalescence prior to the localized failure, which is described by the poro-plasticity formulation. Fluid flow is spread across the lattice network by Darcy’s law in terms of continuous pore pressures, with a special care taken in computing the lattice permeability parameters. This is brought by the proposed discrete model which can capture the initial material heterogeneities, the fracture process zone initiation and the localized failure mechanisms. The macro-cracks can be deflected throughout the medium by the heterogeneities and exhibit irregular geometries. The model theoretical formulation and computational procedure are presented in sufficient detail. Several numerical simulations are given to illustrate a very satisfying performance of the presented discrete model.