On poro-hyperelastic shear

The paper examines the problem of the shear of a porous hyperelastic material, the pore space of which is saturated with an incompressible fluid. Poro-hyperelasticity provides a suitable approach for modelling the mechanical behaviour of highly deformable materials in engineering applications and particularly soft tissues encountered in biomechanical applications. Unlike with the infinitesimal theory of poroelasticity, the application of pure shear generates pore fluid pressures that dissipate with time as fluid migrates either from or into the pore space due to the generated fluid pressure gradients. The analytical results provide benchmark problems that can be used to examine the accuracy of computational approaches.     

Material instabilities of anisotropic saturated multiphase porous media

This paper analyses the material instability of fully saturated multiphase porous media. On account of the fact that anisotropic mechanical behaviours are widely observed in saturated and partially saturated geomaterials, the anisotropic constitutive model developed by Rudnicki for geomaterials is used to model the anisotropic mechanical behaviour of the solid skeleton of saturated porous geomaterials in axisymmetric compression test. The inertial coupling effect between solid skeleton and pore fluid is also taken into account in dynamic cases. Conditions for static instability (strain localisation) and dynamic instability (stationary discontinuity and flutter instability) of fully saturated porous media are derived. The critical modulus, shear band angle for strain localisation, and the bound within which flutter instability may occur are given in explicit forms. The effects of material parameters on material instability are investigated in detail by numerical computations.     

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.     

The elastic moduli estimation of the solid-water mixture

Conceptually, the undrained elastic constants estimated by the poroelasticity theory should be identical to the effective moduli of the two-phase composite of a porous material saturated with pore water. Here we show numerically that the undrained elastic constants determined by an effective moduli estimate are almost identical with those calculated by poroelasticity theory, and if pore shapes are not exactly known and the porosity is around 50%, estimating the elastic constant as the average value of its Voigt and Reuss bounds is reasonably accurate. This is the situation in bone and dentin, the materials that are our primary intended application. This result will hold for situations in which the totally enclosed water phase is constrained to small deformations by virtue of its confinement. Importantly, in this work we assume that water is an isotropic elastic solid with a shear modulus that is 10^?4 times the bulk modulus of the water. Note that it is compressible, but almost incompressible with a Poisson’s ratio of 0.4999.     

Numerical simulation of a fluid-saturated soil under a dynamically loaded pile

The behaviour of the soil under a dynamically loaded pile toe is studied. The soil is modelled as a fluid saturated porous continuum. The constitutive behaviour of the solid skeleton is described by the elasto-plastic model of Drücker-Prager. The wave propagation is simulated with a dynamical finite-element program.A two-phase model of soil gives useful information about effective stress and pore pressure in the soil. In saturated soil the main wave under the pile toe propagates more downards than in dry soil, due to the higher compressional stiffness in saturated soil. The plastic zone under the pile toe propagates with the velocity of the fast compressional wave. The pore fluid influences the plasticity strongly and can be expected to affect pile driving too.The distribution of effective stress and pore pressure under the pile toe depends on the permeability of the soil and cannot be calculated uniquely from a single-phase calculation. Therefore, a nonlinear soil cannot be modelled correctly as a conventional single-phase material.     

不可压饱和多孔弹性简支梁的动力响应

在杆件弯曲小变形的假定下,考虑杆件的侧向变形因素,根据多孔介质理论,本文首先建立了不可压饱和多孔弹性梁弯曲变形时动力响应的控制方程。其次,基于所建立的控制微分方程,利用变量分离法,研究了两端可渗透的饱和多孔弹性简支梁在梁中间集中载荷作用下的动力响应,得到了不同物性参数下简支梁动态弯曲时挠度和孔隙流体压力等效力偶等随时间的响应曲线。研究发现由于孔隙流体和固相骨架的相互作用,不可压饱和多孔弹性梁挠度的动力响应具有粘性特征,同时,随着时间的增加,饱和多孔弹性梁的挠度、弯矩等最终趋于经典弹性梁的静挠度、弯矩,此时,孔隙流体压力为零,梁的固相骨架承担所有的外载荷。     

一维流体饱和粘弹性多孔介质层的动力响应

本文研究了不可压流体饱和粘弹性多孔介质层的一维动力响应问题。基于粘弹性理论和多孔介质理论,在流相和固相微观不可压、固相骨架服从粘弹性积分型本构关系和小变形的假定下,建立了不可压流体饱和粘弹性多孔介质层一维动力响应的数学模型,利用Laplace变换,求得了原初边值问题在变换空间中的解析解,并利用Laplace逆变换的Crump数值反演方法,得到原动力响应问题的数值解。数值研究了饱和标准线性粘弹性多孔介质层的动力响应,分析了固相位移、渗流速度、孔隙压力及固相有效应力等的响应特征。结果表明,与不可压流体饱和弹性多孔介质相同,不可压流体饱和粘弹性多孔介质中亦只存在一个纵波,并且固相骨架的粘性对动力行为有显著的影响。     

Finite element simulation of deformation behaviour of cellular rubber components

This paper presents a new prospect of investigating the mechanical behaviour of cellular rubber using a porous hyperelastic material model within the framework of homogenisation method to consider pore volume fraction. There are number of hyperelastic material models to describe the behaviour of homogeneous elastomer, but very few to characterise the complex properties of cellular rubber. The analysis of dependence of material behaviour on pore density using the new material model is supported with experiments to characterise the actual material behaviour. The finite element simulations are then followed by compression load tests to validate the material model.     

Scattered fracture of porous materials with brittle skeleton

A model of damage accumulation in a porous medium with a brittle skeleton saturated with a compressible fluid is formulated in the isothermal approximation. The model takes account of the skeleton elastic energy transformation into the surface energy of microcracks. In the case of arbitrary deformations of an anisotropic material, constitutive equations are obtained in a general form that is necessary and sufficient for the objectivity and thermodynamic consistency principles to be satisfied. We also formulate the kinetics equation ensuring that the scattered fracture dissipation is nonnegative for any loading history. For small deviations from the initial state, we propose an elastic potential which permits describing the principal characteristics of the behavior of a saturated porous medium with a brittle skeleton. We study the acoustic properties of the material under study and find their relationship with the strength criterion depending on the accumulated damage and the material current deformation. We consider the problem of scattered fracture of a saturated porous material in a neighborhood of a spherical cavity. We show that the cavity failure occurs if the Hadamard condition is violated.     

Fluid flow and fluid shear stress in canaliculi induced by external mechanical loading and blood pressure oscillation

The paper studies the problem of fluid flow and fluid shear stress in canaliculi when the osteon is subject to external mechanical loading and blood pressure oscillation. The single osteon is modeled as a saturated poroelastic cylinder. Solid skeleton is regarded as a poroelastic transversely isotropic material. To get near-realistic results, both the interstitial fluid and the solid matrix are regarded as compressible. Blood pressure oscillation in the Haverian canal is considered. Using the poroelasticity theory, an analytical solution of the pore fluid pressure is obtained. Assuming the fluid in canaliculi is incompressible, analytical solutions of fluid flow velocity and fluid shear stress with the Navier-Stokes equations of incompressible fluid are obtained. The effect of various parameters on the fluid flow velocity and fluid shear stress is studied.     

Effects of large deformation and material nonlinearity on spherical indentation of hyperelastic soft materials

Up to now, the indentation of hyperelastic soft materials has not been completely understood. In this paper, the spherical indentation on hyperelastic soft solids was systematically investigated through theoretical analysis and finite element method (FEM). The validation and application of the Hertzian load-displacement relation for indentation of hyperelastic soft materials were clarified, the effects of large deformation and material nonlinearity on spherical indentation of hyperelastic soft materials were analyzed and discussed. It was found that the complicated indentation behaviors of hyperelastic soft solids mainly depended on the coupling interactions of large deformation and material nonlinearity. Besides, we proposed two new nonlinear elastic contact models to separate the effects of large deformation and material nonlinearity on spherical indentation responses of hyperelastic soft solids. Our efforts might help to enhance the understanding of hyperelastic indentation problems and provided necessary instructions for the mechanical characterization of hyperelastic soft materials.     

Convective cooling/heating induced thermal stresses in a fluid saturated porous medium undergoing local thermal non-equilibrium

Local thermal non-equilibrium (LTNE) may have profound effects on the pore pressure and thermal stresses in fluid saturated porous media under transient thermal loads. This work investigates the temperature, pore pressure, and thermal stress distributions in a porous medium subjected to convective cooling/heating on its boundary. The LTNE thermo-poroelasticity equations are solved by means of Laplace transform for two fundamental problems in petroleum engineering and nuclear waste storage applications, i.e., an infinite porous medium containing a cylindrical hole or a spherical cavity subjected to symmetrical thermo-mechanical loads on the cavity boundary. Numerical examples are presented to examine the effects of LTNE under convective cooling/heating conditions on the temperature, pore pressure and thermal stresses around the cavities. The results show that the LTNE effects become more pronounced when the convective heat transfer boundary conditions are employed. For the cylindrical hole problem of a sandstone formation, the thermally induced pore pressure and the magnitude of thermal stresses are significantly higher than the corresponding values in the classical poroelasticity, which is particularly true under convective cooling with moderate Biot numbers. For the spherical cavity problem of a clay medium, the LTNE effect may become significant depending on the boundary conditions employed in the classical theory.     

Revisiting poromechanics in the context of microporous materials

Poromechanics offers a consistent theoretical framework for describing the mechanical response of porous solids, fully or partially saturated with a fluid phase. When dealing with fully saturated microporous materials, which exhibit pores of the nanometre size, aside from the fluid pressure acting on the pore walls additional effects due to adsorption and confinement of the fluid molecules in the smallest pores must be accounted for. From the mechanical point of view, these phenomena result into volumetric deformations of the porous solid: the so-called “swelling” phenomenon. The present work investigates how the poromechanical theory should be refined in order to describe adsorption and confinement induced swelling in microporous solids. Firstly, we report molecular simulation results that show that the pressure and density of the fluid in the smallest pores are responsible for the volumetric deformation of the material. Secondly, poromechanics is revisited in the context of a microporous material with a continuous pore size distribution. Accounting for the thermodynamic equilibrium of the fluid phase in the overall pore space, the new formulation introduces an apparent porosity and an interaction free energy. We use a prototype constitutive relation relating these two quantities to the Gibbs adsorption isotherm, and then calculate the induced deformation of the solid matrix. Agreement with experimental data found in the literature is observed. As an illustrating example, we show the predicted strains in the case of adsorption of methane on activated carbon.     

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.     

A general approach for defining the macroscopic free energy density of saturated porous media at finite strains under non-isothermal conditions

A general approach is proposed for defining the macroscopic free energy density function (and its complement, the free enthalpy) of a saturated porous medium submitted to finite deformations under non-isothermal conditions, in the case of compressible fluid and solid constituents. Reference is made to an elementary volume treated as an ‘open system’, moving with the solid skeleton. The proposed free energy depends on the generalised strains (namely an appropriate measure of the strain of the solid skeleton and the variation in fluid mass content) and the absolute temperatures of the solid and fluid phases (which are assumed to differ from each other for the sake of generality). This macroscopic energy proves to be a potential for the generalised stresses (namely the associated measure of the total stress and the free enthalpy of the pore fluid per unit mass) and the entropies of the solid and fluid phases. In contrast with mixture theories, the resulting free energy is not the simple sum of the free energies of the single constituents. Two simplified cases are examined in detail, i.e. the semilinear theory (originally proposed for isothermal conditions and extended here to non-isothermal problems) and the linear theory. The proposed approach paves the way to the consistent non-isothermal-hyperelastic-plastic modelling of saturated porous media with a compressible fluid and solid constituents.     

Dynamic Model of a Periodic Medium with Double Porosity

The paper presents a unified mathematical approach for describing the dynamic stressstrain state of mechanical structures from heterogeneous materials possessing a double coupled system of pore channels filled with fluid. New dynamic equations describing the oscillations of poroelastic systems based on the developed model of a continuous medium with additional degrees of freedom in the form of various pressures of the components constituting the liquid phase of the material are obtained. The equations and the method of obtaining them have a greater degree of generalization than those encountered in the literature. Theoretical results can be used to study the propagation of vibrations in fractured geological rocks saturated with liquid, to develop technical systems of new structural materials with a porous structure, for the analysis of micro streams of fluid in the hierarchical system of microporous bone tissue.     

Acoustical characterization of fluid-saturated porous media with local heterogeneities: Theory and application

A linear dynamic model of fully saturated porous media with local (either microscopic or mesoscopic) heterogeneities is developed within the context of Biot’s theory of poroelasticity. Viscoporoelastic behavior associated with local fluid flow is characterized by the notion of the dynamic compatibility condition on the interface between the solid and the fluid. Complex, frequency-dependent material parameters characterizing the viscoporoelasticity are derived. The complex properties can be obtained through determining the quasi-static poroelastic parameters, the properties of individual constituents, and the relaxation time of the dynamic compatibility condition on the interface. Relationships among various quasi-static poroelastic parameters are developed. It is shown that local fluid flow mechanism is significant only in the porous media with local heterogeneities. The relaxation time of the compatibility condition on the interface depends upon the details of local structure of porous media that control local fluid pressure diffusion. The new model is used to describe the velocity dispersion and attenuation in fully saturated porous media. The proposed model provides a theoretical framework to simulate the acoustical behavior of fully saturated porous media over a wide range of frequencies without making any explicit assumption about the structure of local heterogeneities.     

基于孔隙弹性耦合的水力压裂数值模拟研究

系统地阐述基于多孔介质渗流-损伤耦合原理,进行水力压裂的FEM数值实现方法.基本架构为:(1)引入孔隙流体压力膨胀系数将孔隙流体压力与应力场进行耦合;(2)基于损伤局部化模型,提出裂缝张开度表达式;(3)提出水力压裂引起的多孔介质水-力学属性的各向异性表达式;(4)提出全流量加载的耦合分析方案.作为实例,模拟三维地层水压裂缝扩展形态,通过比较模型的数值解和经典理论解,验证该方法的正确性.