Dynamic response of a pile embedded in a porous medium subjected to plane SH waves

In this paper, the frequency domain dynamic response of a pile embedded in a porous medium subjected to SH seismic waves is investigated. The surrounding porous medium of the pile is described by Biot’s poro-elastic theory, while the pile embedded in the porous medium is treated as a beam and described by a beam vibration theory. Using the Hankel transformation method, the fundamental solution for a half-space porous medium subjected to a horizontal circular patch load is established. According to the fictitious pile methodology, the second kind of Fredholm integral equation for the pile is established in terms of the obtained fundamental solution and free wave field. The solution of the integral equation yields the dynamic response of the pile to plane SH waves. Numerical results indicate that the parameters of the porous medium, the pile and incident waves have considerable influences on the dynamic response of the pile and the porous medium.     

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.     

Analytical solutions for dynamic pressures of coupling fluid–porous medium–solid due to SV wave incidence

This paper presents the results of theoretical investigation on the dynamic coupling of an ideal fluid‐porous medium‐elastic half‐space system subjected to SV waves to study the effect of sediment on the seismic response of dams for reservoirs that are deposited with a significant amount of sediment after a long period of operation. The effects of the porous medium and the incident wave angle on dynamic pressures in the overlying ideal fluid are analyzed, and the reflection and transmission coefficients of the wave at the material interfaces are derived using an analytical solution in terms of displacement potentials. The numerical test of modeling shows that the dynamic pressures significantly depend on the properties of porous medium. The fully saturated porous medium reduces the response peaks slightly, while the partially saturated porous medium causes a considerable increase in the resonance peaks. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysis of pile groups in a poroelastic medium subjected to horizontal vibration

This work investigates the dynamic response of pile groups embedded in a poroelastic medium subjected to horizontal loading. The dynamic response is analyzed using the Muki and Sternberg Method. The load transfer problem is formulated in terms of a second-kind Fredholm integral. The dynamic impedance of the pile groups is computed using the pile–soil–pile dynamic interaction factors. The shear force, bending moment and pore pressure is obtained using the superposition method. Numerical results indicate that the pile flexibility ratio and the pile distance have considerable influence on the dynamic response of the piles and the poroelastic medium.     

移动荷载作用下饱和土中单桩的动力响应

根据Biot动力理论,采用Fourier和Hankel变换方法得到了半空间饱和土受移动载荷及土体内受垂直简谐载荷作用下频域内基本解。根据虚拟桩法,得到了移动载荷作用下桩基的第2类Fredholm积分方程,并应用IFFT方法得到时间、空间域内单桩的动力响应。数值结果表明,移动荷载会引起桩身的负摩擦力;桩身最大轴力、孔压随移动荷载速度增加而增大;此外,在桩上端部会出现孔压集中现象。     

饱和土体中排桩对移动荷载的被动隔振效果分析

利用Muki和Sternberg的虚拟桩法,研究了饱和土体地基中排桩对移动荷载引起振动的被动隔振效果。隔振桩作为一维杆,饱和土体满足Biot理论。利用已有的移动荷载作用下的饱和土体的自由波场解及饱和土体内部受竖向圆形分布荷载作用下的基本解,建立了频域内土-桩的第2类Fredholm积分方程。通过Fourier逆变换得到时间域内评价隔振效果的振幅比。与已知文献结果相比较,验证了方法的正确性。数值结果表明,荷载速度对排桩的隔振效果有一定影响,即在相同隔振系统情况下,单排桩对低速荷载引起振动的隔振效果比高速移动荷载效果好。同时,较高速时的最佳桩长比低速时要短。     

Vertical amplitude reduction of Rayleigh waves by a row of piles in a poroelastic half‐space

A row of rigid piles is addressed as the countermeasures for isolating Rayleigh waves in a poroelastic half‐space. The complex characteristic equations for Rayleigh waves are derived via Biot's theory and their existence conditions are given. The piles are modeled as Euler–Bernoulli beams with longitudinal displacements and the diffracted field by each pile is constructed only with Rayleigh waves. Six infinite linear systems of algebraic equations are obtained in terms of the equilibrium of forces and continuity of displacements at the pile–soil interfaces. The systems are subsequently solved in the complex least‐squares sense. The influence of certain pile and soil characteristics such as the permeability of poroelastic soil, spacing between the piles and length of the piles on the isolating performance of a pile barrier is investigated. Computed results show that the permeability of poroelastic soil displays a significant effect on the vertical amplitude reduction of Rayleigh waves. Copyright © 2009 John Wiley & Sons, Ltd.     

层状饱和土体中排桩对简谐荷载隔振效果分析

根据Biot理论,利用已有的传递、透射矩阵法得到层状饱和土体内部受竖向圆形分布荷载作用下的基本解,再由Muki虚拟桩法,建立了频域内层状土-桩的第二类Fredholm积分方程。通过离散方法求解积分方程,得到了评价隔振效果的振幅比。与已知文献结果比较,验证了方法的正确性。数值结果表明：对于同种类型的振源,采用相同的隔振系统,在上软下硬的多层土体中的隔振效果比在上硬下软的土层中要好。采用较长的桩或刚度较大的桩,或桩之间的距离加密都可得到好的隔振效果。     

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.     

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.     

Stresses due to vertical subsurface loading for an inhomogeneous cross‐anisotropic half‐space

In this article, we present the solutions for the stresses induced by four different loads associated with an axially loaded pile in a continuously inhomogeneous cross‐anisotropic half‐space. The planes of cross‐anisotropy are parallel to the horizontal surface of the half‐space, and the Young's and shear moduli are assumed to vary exponentially with depth. The four loading types are: an embedded point load for an end‐bearing pile, uniform skin friction, linear variation of skin friction, and non‐linear parabolic variation of skin friction for a friction pile. The solutions for the stresses due to the pile load are expressed in terms of the Hankel integral and are obtained from the point load solutions of the same inhomogeneous cross‐anisotropic half‐space which were derived recently by the authors (Int. J. Rock Mech. Min. Sci. 2003; 40 (5):667–685). A numerical procedure is proposed to carry out the integral. For the special case of homogeneous isotropic and cross‐anisotropic half‐space, the stresses predicted by the numerical procedure agree well with the solutions of Geddes and Wang (Geotechnique 1966; 16 (3):231–255; Soils Found. 2003; 43 (5):41–52). An illustrative example is also given to investigate the effect of soil inhomogeneity, the type and degree of soil anisotropy, and the four different loading types on the vertical normal stress. The presented solutions are more realistic in simulating the actual stratum of loading problem in many areas of engineering practice. Copyright © 2004 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.     

Interaction of a cylindrical thin‐walled pile with an inhomogeneous isotropic elastic medium

This paper presents a rigorous analysis for the static interaction of a cylindrical thin‐walled pile with an inhomogeneous isotropic elastic half‐space under vertical, horizontal, and torsional forces individually applied at the top of pile. The inhomogeneity is specified with the exponential variation of shear modulus along depth of the embedding medium, and the Poisson's ratio is assumed to be constant. By means of a set of Green's functions for pile and soil medium and satisfying the compatibility conditions between the 2 interacting media, the formulation is reduced to coupled Fredholm integral equations. Using the adaptive‐gradient elements, capable of capturing the singular stress transfer at both ends of the pile, a numerical procedure is developed and utilized for evaluating the relevant integral equations and studying the inhomogeneity effect on the soil‐pile interaction responses. The analysis results have been validated for different soil‐pile modulus ratios under axial load and for a Poisson's ratio of 0.3 under lateral load. The procedure does not consider the nonlinear behavior of the soil medium or plastic yielding in the pile section, and the impact of the unreliable results for the case of high Poisson's ratio is not examined.     

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

BEM analysis of the interaction factor for vertically loaded dissimilar piles in saturated poroelastic soil

This paper focuses on an analysis by the boundary element method (BEM) of the pile-to-pile interaction for pile groups with dissimilar piles of different pile lengths embedded in saturated poroelastic soil. The behaviour of the poroelastic homogeneous soil is governed by Biot’s consolidation equations. The pile–soil system is decomposed into extended soil and fictitious piles. Considering the compatibility of vertical strain between fictitious piles and soil, the second kind of Fredholm integral equations were obtained to predict the axial force and settlement along pile shafts numerically. For the analysis of the interaction factor, two loading conditions for a two-dissimilar-pile system were proposed: (a) only one pile is loaded and (b) each pile is subjected to a load proportional to the pile length. Furthermore, the two-pile system was extended to pile groups with a rigid cap to capture the optimum design where each pile shares the same loading at the pile heads. The optimum results require shortening the peripheral piles and elongating internal piles, and the consolidation effect needs to be considered due to the adjustment of loading distribution among piles.     

层状地基中单桩负摩擦问题积分方程解法

利用Biot固结理论和积分方程方法研究了表面有堆载的层状地基中单桩负摩擦问题。根据层状饱和土的圆形载荷基本解得出了单桩在圆形均布载荷作用下在时间域内的第二类Fredholm积分方程组。运用Laplace变换对上述积分方程组进行简化。再结合传递和刚度矩阵传递到各个层中去,对变换域内的积分方程采用Schapery 逆变换方法得到时域内单桩的近似积分方程。求解积分方程组并进行相应的数值逆变换,就可得出层状地基中的单桩在表面圆形均布载荷作用下的位移、轴力、孔压和桩侧摩阻力随时间的变化情况。计算结果表明,桩侧剪力和孔压分层明显。     

Analysis of pile‐raft foundations under complex loads in layered soils

Considering there is hardly any concerted effort to analyze the pile‐raft foundations under complex loads (combined with vertical loads, horizontal loads and moments), an analysis method is proposed in this paper to estimate the responses of pile‐raft foundations which are subjected to vertical loads, horizontal loads and moments in layered soils based on solutions for stresses and displacements in layered elastic half space. Pile to pile, pile to soil surface, soil surface to pile and soil surface to soil surface interactions are key ingredients for calculating the responses of pile‐raft foundations accurately. Those interactions are fully taken into account to estimate the responses of pile‐raft foundations subject to vertical loads, horizontal loads and moments in layered soils. The constraints of the raft on vertical movements, horizontal movements and rotations of the piles as well as the constraints of the raft on vertical movements and horizontal movements of the soils are considered to reflect the coupled effect on the raft. The method is verified through comparisons with the published methods and FEM. Then, the method is adopted to investigate the influence of soil stratigraphy on pile responses. The study shows that it is necessary to consider the soil non‐homogeneity when estimating the responses of pile‐raft foundations in layered soils, especially when estimating the horizontal responses of pile‐raft foundations. The horizontal loads and the moments have a significant impact on vertical responses of piles in pile‐raft foundations, while vertical loads have little influence on horizontal responses of piles in pile‐raft foundations in the cases of small deformations. The proposed method can provide a simple and useful tool for engineering design. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

饱和土地下源u-P形式解答动力响应计算

一直以来,由Biot孔隙弹性动力方程得到的饱和土地下源Green函数都是u-w形式（u为固相介质位移,w为流相相对于固相的平均位移）。应用两相介质纵波解耦理论,得到了饱和土半空间地下点源荷载的u-P形式（P为孔压）Green函数频域解答;克服了u-w形式Green函数在边界元（BEM）积分时的增根影响。再由Hankel反演,结合Somigliana表象积分,完成BEM计算。并以计算结果分析了地下集中力作用时,饱和土位移、孔压、排水量等动力特性,这对地铁等交通工程、地震工程、土-结构动力相互作用（SSI）的响应计算都具有较重要应用价值。     

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.