Energetic Formulation of Large-Deformation Poroelasticity

The modeling of coupled fluid transport and deformation in a porous medium is essential to predict the various geomechanical process such as CO2 sequestration, hydraulic fracturing, and so on. Current applications of interest, for instance, that include fracturing or damage of the solid phase, require a nonlinear description of the large deformations that can occur. This paper presents a variational energy-based continuum mechanics framework to model large-deformation poroelasticity. The approach begins from the total free energy density that is additively composed of the free energy of the components. A variational procedure then provides the balance of momentum, fluid transport balance, and pressure relations. A numerical approach based on finite elements is applied to analyze the behavior of saturated and unsaturated porous media using a nonlinear constitutive model for the solid skeleton. Examples studied include the Terzaghi and Mandel problems; a gas-liquid phase-changing fluid; multiple immiscible gases; and unsaturated systems where we model injection of fluid into soil. The proposed variational approach can potentially have advantages for numerical methods as well as for combining with data-driven models in a Bayesian framework.

A computational model of interaction between material and interface cracks

A quasi-static model for numerical solution of initiation and propagation of cracks along interfaces or inside materials is developed. The two types of cracks are modelled by the material damage theory with two independent damage parameters introduced. For cracks at the interface, in fact represented by contact of construction components, cohesive or adhesive contact is considered, for which several computational relationships based on energetic formulation exist. Accordingly, the appropriate modelling of bulk damage also includes energy consideration. In terms of cracks, it leads to so called diffuse cracks. The computational approach is referred to as phase field models. These will cause damage in a very narrow band representing the actual crack. The computational analysis provides stress-strain quantities and the damage variables to simulate both interface and material cracks. The proposed mathematical approach has a variational form based on an energetic formulation looking for a kind of weak solution. The solution is approximated by a time stepping procedure, a finite element code, and it utilizes quadratic programming algorithms.

Lateral-torsional buckling of compressed and highly variable cross section beams

In the critical state of a beam under central compression a flexural-torsional equilibrium shape becomes possible in addition to the fundamental straight equilibrium shape and the Euler bending. Particularly, torsional configuration takes place in all cases where the line of shear centres does not correspond with the line of centres of mass. This condition is obtained here about a z-axis highly variable section beam; with the assumptions that shear centres are aligned and line of centres is bound to not deform. For the purpose, let us evaluate an open thin wall C-cross section with flanges width and web height linearly variables along z-axis in order to have shear centres axis approximately aligned with gravity centres axis. Thus, differential equations that govern the problem are obtained. Because of the section variability, the numerical integration of differential equations that gives the true critical load is complex and lengthy. For this reason, it is given an energetic formulation of the problem by the theorem of minimum total potential energy (Ritz-Rayleigh method). It is expected an experimental validation that proposes the model studied.

An Energetic Approach to Rate Independent Delamination with Cohesive Contact - An SGBEM Implementation

A mathematical model for analysis of contact delamination problems has been developed and implemented into the program Matlab by means of the Symmetric Galerkin Boundary Element Method (SGBEM). The SGBEM numerical algorithm enables to exploit an energetic formulation which governs interface rupture. A rate-independent model of the delamination process is obtained, considering an interface damage variable. A numerical algorithm has been proposed to provide maximallydissipative local solution which yields numerically stable time-discretization. The developed 2-dimensional sample example of mathematical model demonstrates the model behaviour and its suitability in many aspects of engineering practise.