Modeling large-deforming fluid-saturated porous media using an Eulerian incremental formulation. (November 2017)

• Record Type: Journal Article
• Authors: Rohan, Eduard; Lukeš, Vladimír
• Is Part Of: Advances in engineering software. Volume 113(2017)
• Page Start: 84
• Page End: 95
• Publication Date: 2017-11
• Subjects: Large deformation -- Fluid saturated hyperelastic porous media -- Updated Lagrangian formulation -- Biot model -- Finite element method; Computer-aided engineering -- Periodicals; Engineering -- Computer programs -- Periodicals; Engineering -- Software -- Periodicals; Periodicals; 620.0028553
• DOI: 10.1016/j.advengsoft.2016.11.003
• Languages: English
• ISSNs: 0965-9978

Abstract

Highlights: Darcy flow saturating the large deforming hyperelastic porous material is described in the deformed configuration. A consistent Eulerian incremental formulation is based on the differentiation of the residual function. The proposed linearization scheme leads to a Biot-type model suitable for the two-scale homogenization. Time discretization using backward finite differences, predictor-corrector algorithm is proposed. Numerical model is implemented in the SfePy code, validated using consolidation tests with transient dynamic effects. Abstract: The paper deals with modeling fluid saturated porous media subject to large deformation. An Eulerian incremental formulation is derived using the problem imposed in the spatial configuration in terms of the equilibrium equation and the mass conservation. Perturbation of the hyperelastic porous medium is described by the Biot model which involves poroelastic coefficients and the permeability governing the Darcy flow. Using the material derivative with respect to a convection velocity field we obtain the rate formulation which allows for linearization of the residuum function. For a given time discretization with backward finite difference approximation of the time derivatives, two incremental problems are obtained which constitute the predictor and corrector steps of the implicit time-integration scheme. Conforming mixed finite element approximation in space is used. Validation of the numerical model implemented in the SfePy code is reported for an isotropic medium with a hyperelastic solid phase. The proposed linearization scheme is motivated by the two-scale homogenization which will provide the local material poroelastic coefficients involved in the incremental formulation.