Three‐dimensional time‐harmonic fundamental solutions for a fluid‐saturated poroelastic half‐space with partially permeable free surface

By virtue of a pair of scalar potentials for the displacement of the solid skeleton and the pore fluid pressure field of a saturated poroelastic medium, an alternative solution method to the Helmholtz decomposition is developed for the wave propagation problems in the framework of Biot's theory. As an application, a comprehensive solution for three‐dimensional response of an isotropic poroelastic half‐space with a partially permeable hydraulic free surface under an arbitrarily distributed time‐harmonic internal force field and fluid sources is developed. The Green's functions for the poroelastic fields, corresponding to point, ring, and disk loads, are reduced to semi‐infinite complex‐valued integrals that can be evaluated numerically by an appropriate quadrature scheme. Analytical and numerical comparisons are made with existing elastic and poroelastic solutions to illustrate the quality and features of the solution. Copyright © 2016 John Wiley & Sons, Ltd.     

Displacement ring load Green's functions for saturated porous transversely isotropic tri‐material full‐space

Based on the Biot's poroelastic theory and using scalar potential functions both the ring load and point load displacement Green's functions for a transversely isotropic saturated porous full‐space composed of an upper half‐space, a finite thickness middle layer and a lower half‐space is analytically presented for the first time. It is assumed that each region consists of a different transversely isotropic material. The equations of poroelastodymanics in terms of the solid displacements and the pore fluid pressure are uncoupled with the help of two scalar potential functions, so that the governing equations for the potential functions are either a second order wave equation or a repeated wave‐heat transfer equation of sixth order. With the aid of Fourier expansion with respect to circumferential direction and Hankel integral transforms with respect to the radial direction in cylindrical coordinate system, the response is determined in the form of line integrals in the real space, followed by theorem of inverse Hankel integral transforms. The solutions degenerate to a single phase elastic material, and the results are compared with previous studies, where an excellent agreement may be observed with the results provided in the literature. Some examples of displacement Green's functions are finally given to illustrate the solution. Copyright © 2016 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Dynamic response of road pavement resting on a layered poroelastic half‐space to a moving traffic load

This paper is dedicated to study the dynamic response of a thin‐plate resting on a layered poroelastic half‐space under a moving traffic load. Based on the dynamic poroelastic theory of Biot, the general solutions of the homogeneous poroelastic foundation are obtained by Fourier translation. By using the transmission and reflection matrices method in the frequency domain, the equivalent stiffness of the layered poroelastic half‐space is presented. Kirchhoff's hypotheses are applied to obtain the vertical displacement of the thin plate. By using the inverse Fourier transform, the time domain solution is obtained. As an example of three layers, the influences of the load velocity, the material properties of poroelastic layers, and the flexural rigidity of the plate on the response of the pavement system are examined. Analyses show that a soft intermediate layer results in the significant increase of vertical displacement of road pavement. Comparison with the existing work validates the present model. Copyright © 2013 John Wiley & Sons, Ltd.     

Fundamental solutions for a fluid‐saturated,transversely isotropic,poroelastic solid

The fundamental solutions were obtained for step‐like point forces acting in three orthogonal directions and an instantaneous fluid point source in a fluid‐saturated, porous, infinite solid of transversely isotropic elasticity and permeability. After expressing the governing equations in the form of matrix in the Laplace space, we employed Kupradze's method together with the triple Fourier transforms. This method reduces the simultaneous partial differential equations with respect to three displacement components and a pore fluid pressure to a differential equation in terms of only one potential scalar function, which can be operationally solved in the transformed space. After the Laplace inversion of the potential, the residue theorem was applied to its Fourier inverse transform with respect to one of the transformation variables. The Fourier transforms with respect to two other variables were rewritten into the Hankel transforms. Copyright © 2002 John Wiley & Sons, Ltd.     

Dynamic response to fluid extraction from a poroelastic half‐space

In this study, the dynamic response of a poroelastic half‐space to a point fluid sink is investigated using Biot's dynamic theory of poroelasticity. Based on Biot's theory, the governing field equations are re‐formulated in frequency domain with solid displacement and pore pressure. In a cylindrical coordinate system, a method of displacement potentials for axisymmetric displacement field is proposed to decouple the Biot's field equations to three scalar Helmholtz equations, and then the general solution to axisymmetric problems are obtained. The full‐space fundamental singular solution for a point sink is also derived using potential methods. The mirror‐image method is finally applied to construct the fundamental solution for a point sink buried in a poroelastic half‐space. Furthermore, a numerical study is conducted for a rock, that is, Berea sandstone, as a representative example. Copyright © 2013 John Wiley & Sons, Ltd.     

Linear coupled analysis of desiccation shrinkage in a double‐layer partially saturated medium: semi‐explicit solutions

Solutions are presented for the problem of isothermal dessiccation shrinkage in a double‐layer porous partially saturated medium. The rheological model taken into account is linear poroelastic. Hence the analysis is mainly focused on hydromechanical coupling effects and contrasts of mechanical and hydraulic properties between two materials: a low thickness skin comprised between the outer boundary and the reference porous material. Three one‐dimensional ideal structures are taken into account: a wall of finite thickness (cartesian geometry), a thick cylinder and a thick sphere. The solution of the time‐dependent problem is arrived at by applying Laplace transforms to the field variables. Exact solutions are obtained in Laplace transform space using Mathematica^© to solve the field equations whilst taking into account the continuity equations at the interface and the boundary conditions. The Talbot's modified algorithm has been performed to invert the Laplace transform solutions. A bibliographical and numerical study shows that this method is remarkably precise, stable and close to the analytical inversion. Results are presented using poroelastic data representative of a concrete material and involve a strong coupling effect between hydraulical and mechanical behaviours. A first approach elastic modelling of degradation process have been presented using a thin outer layer. Apart from emphasising the semi‐explicit solution utility due to accurate speed calculation, this paper deals with more complex problems than those which can be solved using purely analytical solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

二维饱和多孔介质因点汇诱发比奥固结的解析解

给出了有限二维饱和多孔介质因点汇诱发的Biot固结的一个解析解。其中假设多孔介质为均匀各向同性和线弹性,假设孔隙压力场符合第1类边界条件,数学模型采用可压缩多孔介质模型。利用傅里叶和拉普拉斯变换及相应反演获得了双重无穷项级数和形式的精确解。然后特别探讨了定流量点汇诱发的稳态解析解,并用文献现有解析解进行了验证。所提出的解析解适合于验证数值解,并可用于深入分析有限二维多孔介质的流-固耦合行为。     

A solution around a backfilled cavity in a low‐permeability poroelastic medium with application in in situ heating tests

Thermo‐hydro‐mechanical responses around a cylindrical cavity drilled or excavated in a low‐permeability formation are studied when the cavity is subjected to a time‐dependent thermal loading. The cavity is considered backfilled after it is supported by casing or lining. Solutions of temperature, pore water pressure, stress, and displacement responses are analytically formulated based on Biot's consolidation theory with the assumption that the backfilling material, supporting material, and surrounding low‐permeability formation are poroelastic media. The solution is expressed in Laplace space, and numerical inversion techniques are used to find field variables in the real‐time domain. After the solution is verified with the numerical results, it is applied in a large‐scale in situ heating test – PRACLAY heating test – for a predictive reference calculation and an extensive parametric study. Another medium‐scale in situ heating test – ATLAS III heating test – is also analyzed using the solution, which provides reasonable agreement with measurements. The new analytical solution proves to be a convenient tool for a good understanding of the resulting coupled thermo‐hydro‐mechanical behavior and is therefore valuable for the interpretation of measured data in engineering practices and for a rational design of potential radioactive waste repositories. Copyright © 2016 John Wiley & Sons, Ltd.     

Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.     

Dynamic responses of a pile embedded in a layered poroelastic half‐space to harmonic lateral loads

A single pile embedded in a layered poroelastic half‐space subjected to a harmonic lateral load is investigated in this study. Based on Biot's theory, the frequency domain fundamental solution for a horizontal circular patch load applied in the layered poroelastic half‐space is derived via the transmission and reflection matrices method. Utilizing Muki and Sternberg's method, the second kind of Fredholm integral equation describing the dynamic interaction between the layered half‐space and the pile subjected to a top harmonic lateral load is constructed. The proposed methodology is validated by comparing results of this paper with some existing results. Numerical results show that for a two‐layered half‐space, the thickness of the upper softer layer has pronounced influences on the dynamic response of the pile and the half‐space. For a three‐layered half‐space, the presence of a softer middle layer in the layered half‐space will enhance the compliance for the pile significantly, while a stiffer middle layer will diminish the dynamic compliance of the pile considerably. Copyright © 2009 John Wiley & Sons, Ltd.     

Torsional wave propagation in non‐homogeneous layer between non‐homogeneous half‐spaces

The study of surface wave in a layered medium has their possible application in geophysical prospecting. In the present work, dispersion equation for torsional wave in an inhomogeneous isotropic layer between inhomogeneous isotropic half‐spaces has been derived. Two cases are discussed separately for torsional wave propagation in inhomogeneous layer between homogeneous and non‐homogeneous half‐spaces, respectively. Further, two possible modes for torsional wave propagation are obtained in case of inhomogeneous layer sandwiched between non‐homogeneous half‐spaces. Closed form solutions for displacement in the layer and half‐spaces are obtained in each case. The study reveals that the layer width, layer inhomogeneity, frequency of inhomogeneity, as well as inhomogeneity in the half‐space has significant effect on the propagation of torsional surface waves. Displacement and implicit dispersion equation for torsional wave velocities are expressed in terms of Heun functions and their derivatives. Effects of inhomogeneity on torsional wave velocity are also discussed graphically by plotting the dimensionless phase velocity against dimensionless and scaled wave number for different values of inhomogeneity parameter. Copyright © 2012 John Wiley & Sons, Ltd.     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Scaled boundary finite‐element analysis of a non‐homogeneous axisymmetric domain subjected to general loading

The scaled boundary finite‐element method is derived for elastostatic problems involving an axisymmetric domain subjected to a general load, using a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system. The method is particularly well suited to modelling unbounded problems, and the formulation allows a power‐law variation of Young's modulus with depth. The efficiency and accuracy of the method is demonstrated through a study showing the convergence of the computed solutions to analytical solutions for the vertical, horizontal, moment and torsion loading of a rigid circular footing on the surface of a homogeneous elastic half‐space. Computed solutions for the vertical and moment loading of a smooth rigid circular footing on a non‐homogeneous half‐space are compared to analytical ones, demonstrating the method's ability to accurately model a variation of Young's modulus with depth. Copyright © 2003 John Wiley & Sons, Ltd.     

列车荷载下上覆弹性层横观各向同性饱和地基的动力响应

基于Biot多孔弹性介质的波动理论,构造了轨道、道碴、枕木及弹性层的横观各向同性饱和地基在列车荷载下的动力计算模型。利用Fourier变换,得到了列车荷载作用下横观各向同性软土地基上弹性层动力响应的解析结果。利用离散Fourier逆变换得到数值计算结果,分析了荷载速度、地基的各向异性参数、弹性层刚度系数及厚度对位移和孔压响应的影响。分析结果表明：弹性层对控制地基振动作用显著,地表振动幅值随荷载速度的增加而增大,软土的横向弹性模量对地表振动及土中孔压有较大影响。     

Dynamic responses of unsaturated half‐space soil to a moving harmonic rectangular load

The problem of the dynamic responses of a semi‐infinite unsaturated poroelastic medium subjected to a moving rectangular load is investigated analytical/numerically. The dynamic governing equations are obtained with consideration of the compressibility of solid grain and pore fluid, inertial coupling, and viscous drag as well as capillary pressure in the unsaturated soil, and they can be easily degraded to the complete Biot's theory. Using the Fourier transform, the general solution for the equations is derived in the transformed domain, and then a corresponding boundary value problem is formulated. By introducing fast Fourier transform algorithm, the unsaturated soil vertical displacements, effective stresses, and pore pressures induced by moving load are computed, and some of the calculated results are compared with those for the degenerated solution of saturated soils and confirmed. The influences of the saturation, the load speed, and excitation frequency on the response of the unsaturated half‐space soil are investigated. The numerical results reveal that the effects of these parameters on the dynamic response of the unsaturated soil are significant.     

Dynamic responses of a poroelastic half‐space from moving trains caused by vertical track irregularities

The vibrations of railway tracks on a poroelastic half‐space generated by moving trains are investigated through a vehicle–track–ground coupling model. The theoretical model incorporates a vehicle, a track, and a fully saturated poroelastic half‐space soil medium. The source of vibration excitation is divided into two components: the quasi‐static loads and the dynamic loads. The quasi‐static loads are related to the static component of the axle loads, whereas the dynamic loads are due to the dynamic wheel–rail interaction. A linear Hertizian contact spring is introduced between each wheelset and the rail to consider the dynamic loads. Biot's dynamic theory is used to characterize the poroelastic half‐space soil medium. Using the Fourier transform, the governing equations for the track–ground system are solved and the numerical results are presented for a single axle vehicle model. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium are investigated. In addition, the different roles of the moving axle loads and the roughness‐induced dynamic loads are identified. It is concluded that the vibration level of the free field off the track predicted by the poroelastic soil medium is smaller than that predicted by the elastic soil medium for vehicle speed below the Rayleigh wave speed of the poroelastic half‐space, whereas it is larger for vehicle speed above the Rayleigh wave speed. The dynamic loads play an important role in the dynamic responses of the track–ground system. Copyright © 2010 John Wiley & Sons, Ltd.     

横观各向同性饱和地基三维瞬态Lamb问题

研究了横观各向同性饱和土地基在地表动力荷载作用下的三维瞬态响应。基于饱和多孔介质的三维Biot波动理论,利用Laplace变换,建立圆柱坐标系下横观各向同性饱和土的波动方程;解耦波动方程后,根据算子理论,并借助Fourier展开和Hankel变换技术,得到瞬态荷载作用下,饱和土介质的土骨架位移和应力、孔隙水相对位移和孔隙水压力的一般解;利用一般解,给出横观各向同性饱和地基在地表集中荷载激励下的瞬态Lamb问题的解答。数值算例结果表明,采用各向同性饱和介质的动力学模型,不能准确描述具有明显各向异性特性的饱和土地基的瞬态动力特性。     

Analytical study of mine closure behaviour in a poro‐visco‐elastic medium

This paper is interested in the hydro‐mechanical behaviour of an underground cavity abandoned at the end of its service life. It is an extension of a previous study that accounted for a poro‐elastic behaviour of the rock mass (Int. J. Comput. Geomech. 2007; DOI: 10.1016/j.compgeo.2007.11.003 ). Deterioration of the lining support with time leads to the transfer of the loading from the exterior massif to the interior backfill. The in situ material has a poro‐visco‐elastic constitutive behaviour while the backfill is poro‐elastic, both saturated with water. This loading transfer is accompanied by an inward cavity convergence, thereby compressing the backfill, and induces an outward water flow. This leads to a complex space–time evolution of pore pressures, displacements and stresses, which is not always intuitive. In its general setting, a semi‐explicit solution to this problem is developed, using Laplace transform, the inversion being performed numerically. Analytical inversion leading to a quasi‐explicit solution in the time domain is possible by identifying the characteristic creep and relaxation times of volumetric strains with those of the deviatoric strains, on the basis of a parametric study. A few numerical examples are given to illustrate the hydro‐mechanical behaviour of the cavity and highlight the influence of key parameters (e.g. stiffness of backfill, lining deterioration rate, etc.). Further studies accounting for more general material behaviours for the backfill and external ground are ongoing. Copyright © 2008 John Wiley & Sons, Ltd.     

Quasi-static deformation due to two-dimensional seismic sources embedded in an elastic half-space in welded contact with a poroelastic half-space

The Biot linearized theory of fluid saturated porous materials is used to study the plane strain deformation of a two-phase medium consisting of a homogeneous, isotropic, poroelastic half-space in welded contact with a homogeneous, isotropic, perfectly elastic half-space caused by a two-dimensional source in the elastic half-space. The integral expressions for the displacements and stresses in the two half-spaces in welded contact are obtained from the corresponding expressions for an unbounded elastic medium by applying suitable boundary conditions at the interface. The case of a long dip-slip fault is discussed in detail. The integrals for this source are solved analytically for two limiting cases: (i) undrained conditions in the high frequency limit, and (ii) steady state drained conditions as the frequency approaches zero. It has been verified that the solution for the drained case (ω → 0) coincides with the known elastic solution. The drained and undrained displacements and stresses are compared graphically. Diffusion of the pore pressure with time is also studied.