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.     

New Conservation Laws of Energy and C-D Inequalities in Continua Without Microstructure

ＩｎｔｒｏｄｕｃｔｉｏｎＣｏｎｔｉｎｕｕｍｍｅｃｈａｎｉｃｓｉｓｎｏｔｏｎｌｙａｎｏｌｄａｎｄｂｕｔａｌｓｏａｙｏｕｎｇｓｃｉｅｎｔｉｆｉｃｄｉｓｃｉｐｌｉｎｅ.Ｉｔｃｏｎｓｉｓｔｓｏｆｓｏｍｅｆｕｎｄａｍｅｎｔａｌｌａｗｓ,ｗｈｉｃｈａｒｅｖａｌｉｄｆｏｒａｌｌｂｏｄｉｅｓｉｒｒｅｓｐｅｃｔｉｖｅｏｆｔｈｅｉｒｓｈａｐｅｓ,ｃｏｎｓｔｉｔｕｔｉｏｎｓａｎｄｃｏｎｓｔｉｔｕｔｉｖｅｒｅｌａｔｉｏｎｓ,ｗｈｉｃｈｍｕｓｔｒｅｆｌｅｃｔｔｈｅｎａｔｕｒｅｏｆｔｈｅｍａｔｅｒｉａｌａｎｄｔｈｅｃｏｎｓｔｉ…     

Waves in Fractal Media

The term fractal was coined by Benoît Mandelbrot to denote an object that is broken or fractured in space or time. Fractals provide appropriate models for many media for some finite range of length scales with lower and upper cutoffs. Fractal geometric structures with cutoffs are called pre-fractals. By fractal media, we mean media with pre-fractal geometric structures. The basis of this study is the recently formulated extension of continuum thermomechanics to such media. The continuum theory is based on dimensional regularization, in which we employ fractional integrals to state global balance laws. The global forms of governing equations are cast in forms involving conventional (integer-order) integrals, while the local forms are expressed through partial differential equations with derivatives of integer order. Using Hamilton??s principle, we derive the equations of motion of a fractal elastic solid under finite strains. Next, we consider one-dimensional models and obtain equations governing nonlinear waves in such a solid. Finally, we study shock fronts in linear viscoelastic solids under small strains. In all the cases, the derived equations for fractal media depend explicitly on fractal dimensions and reduce to conventional forms for continuous media with Euclidean geometries upon setting the dimensions to integers.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CON-TINUUM THEORIES (Ⅱ)—MICROMORPHIC CONTINUUM THEORY AND COUPLE STRESS THEORY

IntroductionThispaperisadirectcontinuationofRef.[1 ] .InitthecoupledconservationlawofenergypresentedinRef.[2 ]wasextendedandtherathercompletesystemsofbasicbalancelawsandequationsformicropolarcontinuumtheoryhavebeenconstitutedbycombiningtherenewedresultsandthetraditionalconservationlawsofmassandmicroinertiaandtheentropyinequality .Thepurposeofthispaperistorestablishthesystemsofbasicbalancelawsandequationsformicromorphiccontinuumtheoryandcouplestresstheoryviadirecttransitionsandreductionsfromth…     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CON-TINUUM THEORIES (Ⅱ)—MICROMORPHIC CONTINUUM THEORY AND COUPLE STRESS THEORY

The purpose is to reestablish the balance laws of momentum, angular momentum and energy and to derive the corresponding local and nonlocal balance equations for micromorphic continuum mechanics and couple stress theory. The desired results for micromorphic continuum mechanics and couple stress theory are naturally obtained via direct transitions and reductions from the coupled conservation law of energy for micropolar continuum theory, respectively. The basic balance laws and equations for micromorphic continuum mechanics and couple stress theory are constituted by combining these results derived here and the traditional conservation laws and equations of mass and microinertia and the entropy inequality. The incomplete degrees of the former related continuum theories are clarified. Finally, some special cases are conveniently derived. Foundation items: the National Natural Science Foundation of China (10072024); the Research Foundation of Liaoning Education Committee (990111001) Biography: DAI Tian-min (1931≈)     

多孔材料中声波的传播与演化

采用两相多孔介质的拉格朗日模型来描述一种理论流体充填的多孔弹性固体材料,其中孔隙度的变化满足一个附加的平衡方程。     

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

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

A theory of finite elastic consolidation

Presented in this paper is a general theory describing the consolidation of a porous elastic soil. The formulation allows for the occurrence of finite geometry changes and finite elastic strains during the consolidation process. The governing equations have been cast in a rate form and the laws which determine deformation and pore fluid flow, i.e. Hooke's law and Darcy's law, are presented in a frame indifferent manner. A numerical technique is described that provides an approximate solution to the governing equations. The theory and the solution technique are illustrated by several examples of practical interest.     

成层饱和介质平面波斜入射问题的一维化时域方法

地震波斜入射下自由场的输入是大型结构抗震分析中亟待解决的问题之一,尤其是成层饱和多孔介质自由场问题,由于问题的复杂性,目前研究甚少. 本文基于Biot提出的饱和多孔介质动力方程,建立了一种新的求解平面波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法. 该方法首先根据Snell定律将饱和多孔介质二维空间问题转化为一维时域问题,通过对深度方向的有限元离散,得到饱和多孔介质波动问题的一维化有限元方程,然后采用单相弹性介质精确人工边界条件模拟基岩半空间的波动辐射和输入特征,通过考虑基岩与饱和多孔介质间透水或不透水边界条件以及不同饱和多孔介质交界面边界条件,形成基岩上覆成层饱和介质系统的整体有限元方程,最后采用中心差分法与Newmark平均加速度近似格式相结合的方法对时间进行离散,得到节点的动力时程的显式表达. 典型场地的地震反应分析表明,本文方法的计算结果与传递矩阵法结合傅里叶变换的计算结果完全吻合,证明了其有效性.      

An analytical homogenization method for heterogeneous porous materials

The objective of this work is to develop an analytical homogenization method to estimate the effective mechanical properties of fluid-filled porous media with periodic microstructure. The method is based on the equivalent inclusion concept of homogenization applied earlier for solid–solid mixture. It is assumed that porous media are described by the poroelastic constitutive law developed by Biot where porosity is a material parameter. By solving the governing equations of poroelasticity in Fourier transformed domain, the relation between periodic strain and eigenstrain in porous media is established. This relation is subsequently used in an average consistency condition involving both solid and fluid phase stresses and strains. The geometry of the porous microstructure is captured in the g-integral. This homogenization method can also be applied to estimate the equivalent properties of solid–fluid mixture where a pure solid and fluid can be modeled by assuming very low and high porosity, respectively. Several examples are considered to establish this new method by comparing with other existing analytical and numerical methods of homogenization. As an application, poroelastic properties of cortical bone fibril are estimated and compared with previously computed values.     

轴向扩散下简支饱和多孔弹性梁的大挠度分析

基于多孔介质理论和弹性梁的大挠度理论,并考虑轴向变形,在孔隙流体仅沿轴向扩散的假设下,建立了微观不可压饱和多孔弹性梁大挠度弯曲变形的一维非线性数学模型.在此基础上,忽略饱和多孔弹性梁的轴向应变,并利用Galerkin截断法,研究了两端可渗透的简支饱和多孔弹性梁在突加横向均布载荷作用下的拟静态弯曲,给出了饱和多孔梁弯曲时挠度、弯矩和轴力以及孔隙流体压力等效力偶等沿轴线的分布曲线.揭示了大挠度非线性和小挠度线性模型的结果差异,指出大挠度非线性模型的结果小于相应小挠度线性模型的结果,并且这种差异随着载荷的增大而增大.计算表明:当无量纲载荷参数q>5时,应该采用大挠度非线性数学模型进行研究.     

Double porosity in fluid-saturated elastic media: deriving effective parameters by hierarchical homogenization of static problem

We propose a model of complex poroelastic media with periodic or locally periodic structures observed at microscopic and mesoscopic scales. Using a two-level homogenization procedure, we derive a model coherent with the Biot continuum, describing effective properties of such a hierarchically structured poroelastic medium. The effective material coefficients can be computed using characteristic responses of the micro- and mesostructures which are solutions of local problems imposed in representative volume elements describing the poroelastic medium at the two levels of heterogeneity. In the paper, we discus various combinations of the interface between the micro- and mesoscopic porosities, influence of the fluid compressibility, or solid incompressibility. Gradient of porosity is accounted for when dealing with locally periodic structures. Derived formulae for computing the poroelastic material coefficients characterize not only the steady-state responses with static fluid, but are relevant also for quasistatic problems. The model is applicable in geology, or in tissue biomechanics, in particular for modeling canalicular-lacunar porosity of bone which can be characterized at several levels.     

Nonlinear fluid flow laws for orthotropic porous media

Nonlinear fluid flow laws for orthotropic porous media are written in invariant tensor form. As usual in the theory of fluid flow through porous media [1, 2], the equations contain the flow velocity up to the second power. Expressions that determine the nonlinear resistances to fluid flow are presented and it is shown that, on going over from linear to nonlinear flow laws, the asymmetry effect may manifest itself, that is, the fluid flow characteristics may differ along the same straight line in the positive and negative directions. It is shown that, as compared with the linear fluid flow law for orthotropic media when for three symmetry groups a single flow law is sufficient, in nonlinear laws the anisotropy manifestations are much more variable and each symmetry group must be described by specific equations. A system of laboratory measurements for finding the nonlinear flow characteristics for orthotropic porous media is considered.     

Frequency domain fundamental solutions for a poroelastic half-space

In frequency domain, the fundamental solutions for a poroelastic half-space are re-derived in the context of Biot＇s theory. Based on Biot＇s theory, the governing field equations for the dynamic poroelasicity are established in terms of solid displacement and pore pressure. A method of potentials in cylindrical coordinate system is proposed to decouple the homogeneous Biot＇s wave equations into four scalar Helmholtz equations, and the general solutions to these scalar wave equations are obtained. After that, spectral Green＇s functions for a poroelastic full-space are found through a decomposition of solid displacement, pore pressure, and body force fields. Mirror-image technique is then applied to construct the half-space fundamental solutions.Finally, transient responses of the half-space to buried point forces are examined.     

Wave propagation in 3-D poroelastic media including gradient effects

In the framework of the theory of mixtures, the governing equations of motion of a fluid-saturated poroelastic medium including microstructural (for both the solid and the fluid) and micro-inertia (for the solid) effects are derived. This is accomplished by appropriately combining the conservation of mass and linear momentum equations with the constitutive equations for both the solid and the fluid constituents. The solid is assumed to be gradient elastic, that is, its stress tensor depends on the strain and the second gradient of strain tensor. The fluid is assumed to have an analogous behavior, that is, its stress tensor depends on the pressure and the second gradient of pressure. A micro-inertia term in the form of the second gradient of the acceleration of the solid is also included in the equations of motion. The equations of motion in three dimensions are seven equations with seven unknowns, the six displacement components for the solid and the fluid and the pore-fluid pressure. Because of the microstructural effects, the order of these equations is two degrees higher than in the classical case. Application of the divergence and the rot operations on these equations enable one to study the propagation of plane harmonic waves in the infinitely extended medium separately in the form of dilatational and rotational dispersive waves. The effects of the microstructure and the micro-inertia on the dispersion curves are determined and discussed.     

Low-Frequency Asymptotic Analysis of Seismic Reflection from a Fluid-Saturated Medium

Reflection of a seismic wave from a plane interface between two elastic media does not depend on the frequency. If one of the media is poroelastic and fluid-saturated, then the reflection becomes frequency-dependent. This paper presents a low-frequency asymptotic formula for the reflection of seismic plane p-wave from a fluid-saturated porous medium. The obtained asymptotic scaling of the frequency-dependent component of the reflection coefficient shows that it is asymptotically proportional to the square root of the product of the reservoir fluid mobility and the frequency of the signal. The dependence of this scaling on the dynamic Darcy’s law relaxation time is investigated as well. Derivation of the main equations of the theory of poroelasticity from the dynamic filtration theory reveals that this relaxation time is proportional to Biot’s tortuosity parameter.     

A computational framework for fluid–porous structure interaction with large structural deformation

We study the effect of poroelasticity on fluid–structure interaction. More precisely, we analyze the role of fluid flow through a deformable porous matrix in the energy dissipation behavior of a poroelastic structure. For this purpose, we develop and use a nonlinear poroelastic computational model and apply it to the fluid–structure interaction simulations. We discretize the problem by means of the finite element method for the spatial approximation and using finite differences in time. The numerical discretization leads to a system of non-linear equations that are solved by Newton’s method. We adopt a moving mesh algorithm, based on the Arbitrary Lagrangian–Eulerian method to handle large deformations of the structure. To reduce the computational cost, the coupled problem of free fluid, porous media flow and solid mechanics is split among its components and solved using a partitioned approach. Numerical results show that the flow through the porous matrix is responsible for generating a hysteresis loop in the stress versus displacement diagrams of the poroelastic structure. The sensitivity of this effect with respect to the parameters of the problem is also analyzed.

     

Hyperelastic Multiphase Porous Media with Strain-Dependent Retention Laws

This article presents poroelastic laws accounting for a retention behavior dependent also on porosity, as suggested by experimental evidence. Motivated by the numerical formulation of the corresponding boundary-value problem presented in a companion article, these constitutive equations employ displacements and fluid pressures as primary variables. The thermodynamic admissibility of the proposed rate laws for stress and fluid contents is assessed by means of symmetry and Maxwell conditions obtained from the Biot theory. In the case of strain-dependent saturation, the two elasticity tensors describing the drained response in saturated and unsaturated conditions, respectively, are proven to be in general not coincident, with their difference depending on capillary pressure and porosity. Furthermore, it is shown that besides the stress decomposition proposed by Coussy, also the stress split proposed by Lewis and Schrefler is consistent with the Biot framework. The former decomposition is obtained for retention laws depending only on capillary pressure, as expected. The Lewis–Schrefler split is proven to be consistent with retention models depending also on porosity. In these developments, the compressibility of all the phases is taken into account, in order to assess the thermodynamic consistency of an extension of the Biot’s coefficient to partially saturated anisotropic porous media.     

ACOUSTOELASTIC THEORY FOR FLUID-SATURATED POROUS MEDIA

Based on the finite deformation theory of the continuum and poroelastic theory, the aeoustoelastic theory for fluid-saturated porous media （FSPM） in natural and initial coordi- nates is developed to investigate the influence of effective stresses and fluid pore pressure on wave velocities. Firstly, the assumption of a small dynamic motion superimposed on a largely static pre- deformation of the FSPM yields natural, initial, and final configurations, whose displacements, strains, and stresses of the solid-skeleton and the fluid in an FSPM particle could be described in natural and initial coordinates, respectively. Secondly, the subtraction of initial-state equations of equilibrium from the final-state equations of motion and the introduction of non-linear constitu- rive relations of the FSPM lead to equations of motion for the small dynamic motion. Thirdly, the consideration of homogeneous pre-deformation and the plane harmonic form of the small dynamic motion gives an acoustoelastic equation, which provides analytical formulations for the relation of the fast longitudinal wave, the fast shear wave, the slow shear wave, and the slow longitudinal wave with solid-skeleton stresses and fluid pore-pressure. Lastly, an isotropic FSPM under the close-pore jacketed condition, open-pore jacketed condition, traditional unjacketed condition, and triaxial condition is taken as an example to discuss the velocities of the fast and slow shear waves propagating along the direction of one of the initial principal solid-skeleton strains. The detailed discussion shows that the wave velocities of the FSPM are usually influenced by the effective stresses and the fluid pore pressure. The fluid pore-pressure has little effect on the wave velocities of the FSPM only when the components of the applied initial principal solid-skeleton stresses or strains are equal, which is consistent with the previous experimental results.