Finite element implementation of a geometrically and physically nonlinear consolidation model

• Published in: Continuum mechanics and thermodynamics 2023-07, Vol.35 (4), p.1291-1308
• Main Authors: Artamonova, Nina B., Sheshenin, Sergey V.
• Format: Article
• Language: English
• Subjects: Analysis; Boundary value problems; Classical and Continuum Physics; Conservation laws; Continuum mechanics; Coordinates; Deformation; Differential equations; Engineering Thermodynamics; Environmental law; Equilibrium; Fluid pressure; Geometric nonlinearity; Heat and Mass Transfer; Iterative methods; Laws, regulations and rules; Mathematical analysis; Mathematical models; Mechanics; Original Article; Permeability; Physics; Physics and Astronomy; Porosity; Porous materials; Porous media; Solid phases; Structural Materials; Theoretical and Applied Mechanics; Time integration; Viscoelasticity
• ISSN: 0935-1175; 1432-0959
• DOI: 10.1007/s00161-022-01124-5

The paper presents a rather general formulation of a porous medium deformation coupled with a fluid flowing through the pores within the framework of physical and geometric nonlinearity. The boundary-value problem is formulated in terms of the solid phase displacement increment, fluid pressure and porosity increments in the form of differential and variational equations. The equations were derived from the general conservation laws of Continuum Mechanics using spatial averaging over a representative volume element (RVE). The model takes into account the porosity and permeability evolutions during deformation process. The equations of filtration and porosity evolution are formulated in the material coordinate system related to the solid phase, according to the idea of Arbitrary Lagrangian–Eulerian (ALE) approach. The linearization of variational equations was done using Gateaux differentiation technique. The proper finite elements were used for the spatial discretization of the saddle system of equations to satisfy well-known Ladyzhenskaya–Babuška–Brezzi (LBB) correctness condition. A generalization of the implicit time integration scheme with internal iterations at each time step according to the Uzawa method is employed. The convergence of the iterative process is partly theoretically studied. The formulation is numerically implemented in the form of a self-made computer code. Examples of calculations are given.