Macroscopic constitutive equations of thermo-poroviscoelasticity derived using eigenstrains

Macroscopic constitutive equations for thermo-viscoelastic processes in a fully saturated porous medium are re-derived from basic principles of micromechanics applicable to solid multi-phase materials such as composites. Simple derivations of the constitutive relations and the void occupancy relationship are presented. The derivations use the notion of eigenstrain or, equivalently, eigenstress applied to the separate phases of a porous medium. Governing coupled equations for the displacement components and the fluid pressure are also obtained.     

Theory of thin thermoelastic rods made of porous materials

In this paper, we consider thin rods modeled by the direct approach, in which the rod-like body is regarded as a one-dimensional continuum (i.e., a deformable curve) with a triad of rigidly rotating orthonormal vectors attached to each material point. In this context, we present a model for porous thermoelastic curved rods, having natural twisting and arbitrary shape of cross-section. To describe the porosity, we employ the theory of elastic materials with voids. The basic laws of thermodynamics are applied directly to the one-dimensional continuum, and the nonlinear governing equations are established. We formulate the constitutive equations and determine the structure of constitutive tensors. We prove the uniqueness of solution to the boundary-initial-value problem associated with the deformation of porous thermoelastic rods in the framework of linear theory. Then, we show the decoupling of the bending-shear and extension-torsion problems for straight porous rods. Using a comparison with three-dimensional equations, we identify and give interpretations to the relevant fields introduced in the direct approach. Finally, we consider the case of orthotropic materials and determine the constitutive coefficients for deformable curves in terms of three-dimensional constitutive constants by means of comparison between simple solutions obtained in the two approaches for porous thermoelastic rods.     

An elastoplasticity model in porous media theories

In this article, porous media theories are referred to as mixture theories extended by the well-known concept of volume fractions. This approach implies the diverse field functions of both the porous solid matrix and the pore fluid to be represented by average functions of the macroscale.The present investigations are based on a binary model of incompressible constituents, solid skeleton, and pore liquid, where, in the constitutive range, use is made of the second-grade character of general heterogeneous media. Within the framework of geometrically finite theories, the paper offers a set of constitutive equations for the solid matrix, the viscous pore liquid and the different interactions between the constituents. The constitutive model applies to saturated as well as to empty solid materials, taking into account the physical nonlinearities based on elasto-plastic solid deformations. In particular, the constitutive model concentrates on granular materials like soil or concrete, where the elastic deformations are usually small and the plastic range is governed by kinematically hardening properties.     

Wave propagation in fractured porous media

A theory of wave propagation in fractured porous media is presented based on the double-porosity concept. The macroscopic constitutive relations and mass and momentum balance equations are obtained by volume averaging the microscale balance and constitutive equations and assuming small deformations. In microscale, the grains are assumed to be linearly elastic and the fluids are Newtonian. Momentum transfer terms are expressed in terms of intrinsic and relative permeabilities assuming the validity of Darcy's law in fractured porous media. The macroscopic constitutive relations of elastic porous media saturated by one or two fluids and saturated fractured porous media can be obtained from the constitutive relations developed in the paper. In the simplest case, the final set of governing equations reduce to Biot's equations containing the same parameters as of Biot and Willis.Now at Izmir Institute of Technology, Anafartalar Cad. 904, Basmane 35230, Izmir, Turkey.     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY (Ⅱ)—LINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

IntroductionInpart(Ⅰ )ofthework[1],byuseofmixturetheory ,thenonlinearconstitutiveequationsandthefieldequationsofunsaturatedsoilwereconstructed ,andthecompleteequationsforthethermodynamicsystemofunsaturatedsoilwasformed .Inthispart,thelinearconstitutiveequationsandfieldequationsofunsaturatedsoilareobtainedthroughlinearizingnonlinearequations,andthelinearequationsarewrittenintheformssimilartoBiot’sequationsforsaturatedporousmedia .ItisprovedthatDarcy’slawissuitabletodescribethemotionofliquid…     

On mechanics of saturated granular materials

A continuum theory of saturated granular materials is formulated. The basic balance laws for the solid phase as well as for the fluid phase are presented. The constitutive equations are derived and the basic equations of motion of the solid and fluid continua are obtained. Several cases of interest, such as incompressible granules saturated with liquids are discussed. It is shown that the theory contains, as its special cases, the Mohr-Coulomb criterion for a granular material as well as Darcy's law of flow through porous media.     

多孔功能梯度材料Timoshenko梁的自由振动分析

基于Timoshenko梁理论研究多孔功能梯度材料梁(FGMs)的自由振动问题.首先,考虑多孔功能梯度材料梁的孔隙率模型,建立了两种类型的孔隙分布.其次,基于Timoshenko梁变形理论,给出位移场方程、几何方程和本构方程,利用Hamilton原理推导多孔功能梯度材料梁的自由振动控制微分方程,并进行无量纲化,然后应用微分变换法(DTM)对无量纲控制微分方程及其边界条件进行变换,得到含有固有频率的等价代数特征方程.最后,计算了固定-固定(C-C)、固定-简支(C-S)和简支-简支(S-S)三种不同边界下多孔功能梯度材料梁自由振动的无量纲固有频率.将其退化为均匀材料与已有文献数据结果对照,验证了正确性.讨论了孔隙率、细长比和梯度指数对多孔功能梯度材料梁无量纲固有频率的影响.     

Study of vibrational characteristics of poroelastic mechanical systems

The paper deals with theoretical problems of analysis of forced harmonic vibrations in liquid-saturated porous structures. The differential equations of motion written for the vector of the solid phase displacements and the liquid phase pressure are derived from the equations of phase component dynamics and the constitutive equations of anisotropic continuum. An example of transverse vibrations of a porous framing is used to study the influence of material constants on the dynamic characteristics of a poroelastic system. It is shown that an increase in the excitation frequency significantly increases the effect of inertial interaction between the phases of the poroelastic material, especially for the amplitudes of the liquid pressure in the pores. Thus, to obtain exact solutions of problems of poroelastic material dynamics, it is necessary to take into account all types of interaction between the solid and liquid phases of heterogenous materials.     

大变形下初始斜交异性本构方程

采用材料主轴法,建立了初始斜交异性材料在变形构形（Euler描述）下的斜交异性本构方程,以及在初始构形（Lagrange描述）下的形式。具体给出了斜交异性线弹性材料方程的显式,它在Lagrange描述下形式简洁,可方便地用于有限元计算。文中指出,在变形构形下是线弹性的材料,在Lagrange描述下其本构方程一般已成为非线性,我们称之为本构转换非线性。这种非线性在实际的有限元计算中还未引起重视。为理论简明,本构方程是对二维给出的。     

Discontinuous bifurcation analysis of thermodynamically consistent gradient poroplastic materials

In this work, analytical and numerical solutions of the condition for discontinuous bifurcation of thermodynamically consistent gradient-based poroplastic materials are obtained and evaluated. The main aim is the analysis of the potentials for localized failure modes in the form of discontinuous bifurcation in partially saturated gradient-based poroplastic materials as well as the dependence of these potentials on the current hydraulic and stress conditions. Also the main differences with the localization conditions of the related local theory for poroplastic materials are evaluated to perfectly understand the regularization capabilities of the non-local gradient-based one. Firstly, the condition for discontinuous bifurcation is formulated from wave propagation analyses in poroplastic media. The material formulation employed in this work for the spectral properties evaluation of the discontinuous bifurcation condition is the thermodynamically consistent, gradient-based modified Cam Clay model for partially saturated porous media previously proposed by the authors. The main and novel feature of this constitutive theory is the inclusion of a gradient internal length of the porous phase which, together with the characteristic length of the solid skeleton, comprehensively defined the non-local characteristics of the represented porous material. After presenting the fundamental equations of the thermodynamically consistent gradient based poroplastic constitutive model, the analytical expressions of the critical hardening/softening modulus for discontinuous bifurcation under both drained and undrained conditions are obtained. As a particular case, the related local constitutive model is also evaluated from the discontinuous bifurcation condition stand point. Then, the localization analysis of the thermodynamically consistent non-local and local poroplastic Cam Clay theories is performed. The results demonstrate, on the one hand and related to the local poroplastic material, the decisive role of the pore pressure and of the volumetric non-associativity degree on the location of the transition point between ductile and brittle failure regimes in the stress space. On the other hand, the results demonstrate as well the regularization capabilities of the non-local gradient-based poroplastic theory, with exception of a particular stress condition which is also evaluated in this work. Finally, it is also shown that, due to dependence of the characteristic lengths for the pore and skeleton phases on the hydraulic and stress conditions, the non-local theory is able to reproduce the strong reduction of failure diffusion that takes place under both, low confinement and low pore pressure of partially saturated porous materials, without loosing, however, the ellipticity of the related differential equations.     

Thermodynamic theory for porous piezoelectric materials

In the context of a temperature rate-dependent formulation of thermodynamics for porous piezoelectric materials, we derive the basic equations of the linear theory by discussing restrictions on constitutive equations with the help of an antropy production inequality proposed by Green and Laws. Some basic theorems concerning reciprocity and uniqueness are also established.

---

Sommario Nell'ambito di una teoria della termoelasticità generalizzata si stabiliscono le equazioni di base per mezzi piezoelecttrici porosi. Si studiano, inoitre, alcuni teoremi di reciprocità e unicità.

     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY (Ⅱ)-LINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

The linear constitutive equations and field equations of unsaturated soils were obtained through linearizing the nonlinear equations given in the first part of this work. The linear equations were expressed in the forms similar to Biot’s equations for saturated porous media. The Darcy’s laws of unsaturated soil were proved. It is shown that Biot’s equations of saturated porous media are the simplification of the theory. All these illustrate that constructing constitutive relation of unsaturated soil on the base of mixture theory is rational.     

A Micropolar Theory of Porous Media: Constitutive Modelling

The extension of the classical mixture theory by the concept of volume fractions leads to the theory of porous media. In this article, the theory of porous media is generalised to micropolar constituents. The kinematic relations and the balance equations for a porous medium are developed without restricting the number of constituents. Based on the entropy inequality, the general form of the constitutive equations are derived for a binary medium consisting of a porous elastic skeleton saturated by a viscous pore-fluid. Both constituents are assumed to be compressible. Handling the saturation constraint by a Lagrangian multiplier leads to a compatibility of the proposed model to so-called hybrid and incompressible models.     

不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

On localization in saturated porous continua

Localization in elastic-plastic saturated porous media is investigated here using a linear perturbation approach. The adopted localization criterion corresponds to unbounded growth of perturbations. The critical conditions are compared with those obtained by a classical band analysis. While for one phase materials these conditions coincide, in the present context the linear perturbation approach leads in the limit of unbounded growth to the singularity of the undrained acoustic tensor, while the band analysis leads to the singularity of the drained acoustic tensor. Some general results clarifying the hierarchy of these two conditions are established for a quite general set of constitutive equations.     

Infiltration of initially dry, deformable porous media

The present paper studies the infiltration of an incompressible liquid in an initially dry (or partially dry), deformable sponge-like material made of an incompressible constituent in the slug-flow approximation having in mind the application to some industrial processes involving flow through sponge-like materials and, in particular, some composite materials manufacturing processes. The resulting initial-boundary value problem is of Stefan type, with suitable interface conditions and evolution equations describing the position of the interfaces delimiting the saturated region within the porous material. Different models are then suggested in the saturated region, depending on the importance of the inertial terms and on the constitutive equation for the stress. Comparison of the simulation with known experimental results is satisfactory.     

微孔洞生长对板材拉伸过程中变形局部化的影响

本文在作者提出的含孔洞材料下限本构方程的基础上，采用了初始缺陷带模型对微孔洞生长及分布对板材拉伸过程中变形局部影响进行了，分析着重研究了细观损演化规律对变形局部化模式及临界应变的影响，并成功预测了ＡＩＳＩ４３４０钢板材拉伸试件变形局部化换稳为及失稳方向。     

The constitutive equations of finite strain poroelasticity in the light of a micro-macro approach

After recalling the constitutive equations of finite strain poroelasticity formulated at the macroscopic level, we adopt a microscopic point of view which consists of describing the fluid-saturated porous medium at a space scale on which the fluid and solid phases are geometrically distinct. The constitutive equations of poroelasticity are recovered from the analysis conducted on a representative elementary volume of porous material open to fluid mass exchange. The procedure relies upon the solution of a boundary value problem defined on the solid domain of the representative volume undergoing large elastic strains. The macroscopic potential, computed as the integral of the free energy density over the solid domain, is shown to depend on the macroscopic deformation gradient and the porous space volume as relevant variables. The corresponding stress-type variables obtained through the differentiation of this potential turn out to be the macroscopic Boussinesq stress tensor and the pore pressure. Furthermore, such a procedure makes it possible to establish the necessary and sufficient conditions to ensure the validity of an ‘effective stress’ formulation of the constitutive equations of finite strain poroelasticity. Such conditions are notably satisfied in the important case of an incompressible solid matrix.     

On the theory of porous elastic rods

We consider the direct approach to the theory of rods, in which the thin body is modelled as a deformable curve with a triad of rigidly rotating orthonormal vectors attached to every material point. In this context, we employ the theory of elastic materials with voids to describe the mechanical behavior of porous rods. First, we derive the dynamical nonlinear field equations of the model. Then, in the framework of linear theory, we prove the uniqueness of the solution to the associated boundary-initial-value problem. We identify the relevant field quantities from the theory of directed curves by comparison with the three-dimensional equations of straight porous rods. Finally, for orthotropic and homogeneous rods, we determine the constitutive coefficients in terms of the three-dimensional elasticity constants by solving several problems in the two different approaches.