A contact problem for a poroelastic halfspace containing an embedded inextensible membrane

The paper examines the axisymmetric problem of the indentation of a poroelastic halfspace that is reinforced with an inextensible permeable/impermeable membrane located at a finite depth by a rigid indenter. The constitutive behavior of the poroelastic halfspace is described by the three-dimensional theory of poroelasticity proposed by M.A. Biot. The contact conditions between the indenter and the poroelastic halfspace are varied to accommodate both adhesive/frictionless contact and impermeable/permeable conditions. The formulation of the mixed boundary value problems uses the stress function approaches applicable to semi-infinite domains. Successive applications of Laplace and Hankel integral transforms are used to reduce the mixed boundary value problems to sets of coupled Fredholm integral equations of the second kind. These integral equations are solved using numerical approaches, applicable both for the solution of the systems of coupled equations and for Laplace transform inversion, to examine the time-dependent displacement of the rigid indenter. The analytical-numerical estimates for the time-dependent displacements of the rigid indenter are compared with results obtained using a finite element approach.     

Vertical vibrations of a rigid disk embedded in a poroelastic medium

This paper considers the steady-state vertical vibrations of a rigid circular disk embedded at a finite depth below the free surface of a poroelastic medium. Biot's elastodynamic theory for porous media is used in the analysis. General solutions for axisymmetric poroelastic fields are obtained by using Hankel integral transforms. Analytical solutions for influence functions corresponding to four types of buried axisymmetric excitations are derived. The embedded disk problem is fomulated in terms of a set of coupled integral equations for unknown traction and pore pressure jumps across the disk. The kernel functions of the integral equations are the influence functions corresponding to buried vertical, radial and pore pressure ring loads. The system of integral equations is solved numerically by discretizing the disk into several concentric annular rings. Selected numerical solutions for displacements, vertical stress and pore pressure due to a buried fully flexible disk (uniform pressure) are also presented. The vertical compliances of a rigid disk are examined for different depths of embedment, poroelastic materials and hydraulic boundary conditions. Solutions for traction and pore pressure jumps are also examined. The present results are useful in the study of dynamic response of embedded foundations and anchors in poroelastic soils. Copyright © 1999 John Wiley & Sons, Ltd.     

Transient response of a circular cavity in a poroelastic medium

This paper considers the transient response of a pressurized long cylindrical cavity in an infinite poroelastic medium. To obtain transient solutions, Biot's equations for poroelastodynamics are specialized for this problem. A set of exact general solutions for radial displacement, stresses, pore pressure and discharge are derived in the Laplace transform space by using analytical techniques. Solutions are presented for three different types of prescribed transient radial pressures acting on the surface of a permeable as well as an impermeable cavity surface. Time domain solutions are obtained by inverting Laplace domain solutions using a reliable numerical scheme. A detailed parametric study is presented to illustrate the influence of poroelastic material parameters and hydraulic boundary conditions on the response of the medium. Comparisons are also presented with the corresponding ideal elastic solutions to portray the poroelastic effects. It is noted that the maximum radial displacement and hoop stress at the cavity surface are substantially higher than the classical static solutions and differ considerably from the transient elastic solutions. Time histories and radial variations of displacement, hoop stress, pore pressure and fluid discharge corresponding to a cavity in two representative poroelastic materials are also presented.     

A note on the consolidation settlement of a rigid circular foundation on a poroelastic halfspace

This paper uses Biot's poroelasticity approach to examine the consolidation behaviour of a rigid foundation with a frictionless base in contact with a poroelastic halfspace. The mathematical development of the mixed boundary value problem involves a set of dual integral equations in the Laplace transform domain which cannot be conveniently solved by employing conventional procedures. In this paper, a numerical solution is developed using a scheme where the contact normal stress is approximated by a discretized equivalent. The influence of limiting drainage boundary conditions at the entire surface of the halfspace on the degree of consolidation of the rigid circular foundation is investigated. The results obtained in this study are compared with the corresponding results given in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

层状可压缩岩基三维固结问题的状态空间解

采用状态空间法求解层状可压缩岩基的三维固结问题。首先从直角坐标系下考虑可压缩性的三维Biot固结问题的控制方程出发,通过Laplace-Fourier变换得到状态空间方程,解此方程并通过Cayley-Hamilton定理,得到单层可压缩岩基三维固结问题的传递矩阵;然后利用传递矩阵法,结合层间连续性条件和边界条件,得到了层状可压缩岩基三维固结问题在积分变换域内的解答;最后应用Laplace-Fourier逆变换技术,得到层状可压缩岩基三维固结问题在物理域内的理论解答。编制了相应的计算程序,进行了数值计算与分析,证明了压缩性对岩基固结问题的影响     

Surface load dynamic solution of saturated transversely isotropic multilayer half-space

This paper treats the dynamic response of a multilayered transversely isotropic fluid saturated poroelastic half-space under surface time-harmonic traction. The governing system of partial differential equations is uncoupled with the use of a set of physically meaningful and complete potential functions that decompose different body waves in a saturated poroelastic transversely isotropic medium. After expressing the equations in the Hankel-Fourier domain, a proper algebraic factorization is applied to generate reflection and transmission matrices for decomposed waves. All responses including displacements, stresses, and pore fluid pressure for both general patch load and point load are presented in the form of semi-infinite line integrals. The verification of the method is confirmed with the degeneration of the solutions presented here to the existing solutions for dried both homogeneous and multilayered elastic half-spaces as well as poroelastic half-space. Selected numerical results are depicted to investigate the effects of layering and pore pressure on responses of a transversely isotropic poroelastic medium. The load distribution effects are studied by comparison of the patch and point load responses. Also, resonance notion and effective parameters on this phenomenon such as layering system and anisotropy contrast are discussed. Significant influence of materials and layering configuration on number and amplitude of resonances depicted through the numerical evaluation.     

A high‐frequency open boundary for transient seepage analyses of semi‐infinite layers by extending the scaled boundary finite element method

A high‐frequency open boundary has been developed for the transient seepage analyses of semi‐infinite layers with a constant depth. The scaled boundary finite element equation of pore water pressure is formulated first in the frequency domain. With the eigenvalue problem, the equation can be decoupled into modal equations whose modal dynamic permeability equation can be determined. The continued fraction technique is adopted to formulate the continued fraction solution in the frequency domain. All constants in the solution are determined recursively at the high‐frequency limit. By introducing auxiliary variables and the continued fraction solution to the relationship between the prescribed seepage flow and the pore water pressure in the frequency domain, the open boundary condition is obtained. After transformed to the time domain, the open boundary condition is expressed as a system of fractional differential equations. No convolution integral is required. The accuracy of the analysis results increases with the increasing orders of continued fraction. Copyright © 2015 John Wiley & Sons, Ltd.     

On the indentation of a poroelastic layer

The paper examines the axisymmetric contact problem related to the indentation of a fluid saturated poroelastic layer by a smooth rigid punch. The layer rests in bonded contact with a rigid impermeable base and the surface of the layer is considered to be either permeable or impermeable. The paper develops the integral equations governing the problem for the generalized case where the pore fluid exhibits compressibility. The numerical results presented in the paper illustrate the influence of the relative layer thickness, drainage conditions and the compressibility of the pore fluid on the degree of consolidation settlement of the indenting punch.     

Transient dynamic response of multilayered saturated media subjected to impulsive loadings

This paper presents a stable and efficient method for calculating the transient solution of layered saturated media subjected to impulsive loadings by means of the analytical layer element method. Starting with the field equations based on Biot's linear theory for porous, fluid‐saturated media, and the seepage continuity equation, an analytical layer element for a single layer is established by applying Laplace‐Hankel integral transform. The global stiffness matrix in the transform domain for a layered saturated half‐space subjected to a transient circular patch loading is obtained by assembling the layer elements of each layer. The displacements in the time domain are derived by Laplace‐Hankel inverse transform of the global stiffness matrix. Numerical examples are conducted to verify the accuracy of the method and to demonstrate the influences of type of transient loading, buried depth of loading, permeability, and stratification of materials on the transient response of the multilayered saturated poroelastic media.     

The method of fundamental solution for 3‐D wave scattering in a fluid‐saturated poroelastic infinite domain

By using a complete set of poroelastodynamic spherical wave potentials (SWPs) representing a fast compressional wave P_I, a slow compressional wave P_II, and a shear wave S with 3 vectorial potentials (not all are independent), a solution scheme based on the method of fundamental solution (MFS) is devised to solve 3‐D wave scattering and dynamic stress concentration problems due to inhomogeneous inclusions and cavities embedded in an infinite poroelastic domain. The method is verified by comparing the result with the elastic analytical solution, which is a degenerated case, as well as with poroelastic solution obtained using other numerical methods. The accuracy and stability of the SWP‐MFS are also demonstrated. The displacement, hoop stress, and fluid pore pressure around spherical cavity and poroelastic inclusion with permeable and impermeable boundary are investigated for incident plane P_I and SV waves. The scattering characteristics are examined for a range of material properties, such as porosity and shear modulus contrast, over a range of frequency. Compared with other boundary‐based numerical strategy, such as the boundary element method and the indirect boundary integral equation method, the current SWP‐MFS is a meshless method that does not need elements to approximate the geometry and is free from the treatment of singularities. The SWP‐MFS is a highly accurate and efficient solution methodology for wave scattering problems of arbitrary geometry, particularly when a part of the domain extends to infinity.     

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

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

层状饱和地基上刚性圆板的扭转振动及参数分析

用解析的方法首次研究了层状饱和地基上刚性圆板的扭转振动特性.运用Hankel变换求解饱和介质动力问题的控制方程,按混合边值条件建立了层状饱和地基上刚性圆板扭转振动的第二类Fredholm积分方程.数值算例给出了动力柔度系数和扭转角幅值随无量纲频率的变化曲线,与单相弹性及匀质饱和介质情况进行了对比分析,并进行了参数研究.     

瞬变弹性动力问题的一个边界元法

本文给出了瞬变弹性动力问题的一个边界元法。该法是利用威尔逊——θ法的差分公式,把运动方程化为椭圆型微分方程。根据贝蒂定理和动力点荷载的特解,可获得动力问题的边界积分方程。这个解法是在真实时间域内逐步求解的,不需要使用拉氏变换。一个应力波传播的数值算例证实了该方法使用方便,且解答精度较高。     

A quadratic variation indirect boundary element method for traction boundary-value problems of two-dimensional elastostatics

Boundary integral equations for traction boundary-value problems of two-dimensional elastostatics are derived by the indirect boundary element method. Quadratic variation functions for the representation of geometry, fictitious forces and displacements over each boundary element are described. A system of equations approximating to the boundary integral equations is obtained by a Galerkin formulation in which the integral equation is written at Gauss integration points of elements. The method of computation of the Cauchy principal value is described. Examples of application to the analysis of stress and displacement around underground excavations demonstrate the accuracy and efficiency of the formulation.     

Elasto-dynamic Problem

A boundary element method formulation in the real time domain in the transient elastodynamic problems is presented. The equation motion of the linear isotropic elastiic bodare transformed into the elliptic partial differential eqnations by means of the differential formulations of the Wilson-θ method. By utilizing the particular solution corresponding to a dynamic point force in an infinite medium and the Betti’s theorem the bundary integral equations may be obtain. As this approach formulates and solves the problem in real time domain in conjunction with a time-stepping algorithm, it does not require a nimerical inversion.Which transformes the results from the Laplace or Fourier ransformed domain to the time domain. A numerical example was shown in order to demonstrate the applicability of the present method to the practical problems.     

Deformation of a poroelastic layer overlying an elastic half‐space due to dip‐slip faulting

An analytical solution of the plane strain problem of the deformation of a homogeneous, isotropic, poroelastic layer of uniform thickness overlying a homogeneous, isotropic, elastic half‐space due to two‐dimensional seismic sources buried in the elastic half‐space has been obtained. The integral expressions for the displacements, stresses and pore pressure have been obtained using the stress function approach by applying suitable boundary conditions at the free surface and the interface. The solution obtained is in the Laplace–Fourier transform domain. The case of a vertical dip‐slip line dislocation for the oceanic crust model of Earth is studied in detail. Schapery's formula is used for the Laplace inversion and the extended Simpson's formula for the Fourier inversion. Diffusion of pore pressure in the layer is studied numerically. Contour maps showing the pore pressure in the poroelastic layer have been plotted. The effect of the compressibility of the solid and fluid constituents on pore pressure has also been studied. Copyright © 2015 John Wiley & Sons, Ltd.     

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.     

Vertical dynamic response of a rigid foundation embedded in a poroelastic soil layer

A simplified analytical method is presented for the vertical dynamic analysis of a rigid, massive, cylindrical foundation embedded in a poroelastic soil layer. The foundation is subjected to a time‐harmonic vertical loading and is perfectly bonded to the surrounding soil in the vertical direction. The soil underlying the foundation base is represented by a single‐layered poroelastic soil based on rigid bedrock while the soil at the side of the foundation is modeled as an independent poroelastic layer composed of a series of infinitesimally thin layers. The behavior of the soil is governed by Biot's poroelastodynamic theory and its governing equations are solved by the use of Hankel integral transform. The contact surface between the foundation base and the soil is smooth and fully permeable. The dynamic interaction problem is solved following standard numerical procedures. The accuracy of the present solution is verified by comparisons with the well‐known solutions obtained from other approaches for both the elastodynamic interaction problem and poroelastodynamic interaction problem. Numerical results for the vertical dynamic impedance and response factor of the foundation are presented to demonstrate the influence of nondimensional frequency of excitation, soil layer thickness, poroelastic material parameters, depth ratio and mass ratio on the dynamic response of a rigid foundation embedded in a poroelastic soil layer. Copyright © 2008 John Wiley & Sons, Ltd.     

Seismic wave propagation in laterally inhomogeneous poroelastic media via BIEM

This work addresses in‐plane pressure P and vertically polarized shear SV seismic wave propagation in a finite, laterally inhomogeneous, multilayered poroelastic geological region resting on the homogeneous elastic half‐space. The particular approach followed here is based on a combination of the (i) viscoelastic approximation (isomorphism) to Biot's equations of dynamic poroelasticity and on the (ii) boundary integral equation method (BIEM) using frequency‐dependent fundamental solutions of the governing wave equations. The problem is formulated under plane strain conditions and time‐harmonic motions are assumed. Validation of the viscoelastic isomorphism and verification of the BIEM is done by solution of benchmark examples. These simulation studies reveal that the proposed methodology is able to depict a sensitivity of the seismic signals recovered to the following parameters: (i) poroelastic properties of fluid saturated layers; (ii) lateral geological inhomogeneity; (iii) surface topography and (iv) frequency content and direction of the incident wave. It is concluded that the combination of viscoelastic isomorphism with BIEM software provides an effective numerical tool for evaluating site‐effect phenomena in multilayered, fluid saturated geological regions with complex geometry. The numerical results obtained demonstrate that dynamic poroelasticity interacting with other physical peculiarities of the Earth's surface layers, such as lateral heterogeneity, material properties along the wave path, local geological profile and type of elastic wave, gives rise to complex seismic signals on the free surface at the site of interest. Copyright © 2010 John Wiley & Sons, Ltd.     

Consolidation of a poroelastic half-space with anisotropic permeability and compressible constituents by axisymmetric surface loading

The fully coupled Biot quasi-static theory of linear poroelasticity is used to study the consolidation of a poroelastic half-space caused by axisymmetric surface loads. The fluid and solid constituents of the poroelastic medium are compressible and its permeability in the vertical direction is different from its permeability in the horizontal direction. An analytical solution of the governing equations is obtained by taking the displacements and the pore pressure as the basic state variables and using a combination of the Laplace and Hankel transforms. The problem of an axisymmetric normal load is discussed in detail. An explicit analytical solution is obtained for normal disc loading. Detailed numerical computations reveal that the anisotropy in permeability as well as the compressibilities of the fluid and solid constituents of the poroelastic medium have significant effects on the consolidation of the half-space. The anisotropy in permeability may accelerate the consolidation process and may lead to a dilution in the theoretical prediction of the Mandel-Cryer effect. The compressibility of the solid constituents may also accelerate the consolidation process.