A half-space saturated poro-elastic medium subjected to a moving point load

In this paper, an analytical solution for the dynamic response of a half-space porous medium subjected to a moving point load is derived. In the model, the displacements of the solid skeleton and the pore pressure are expressed in terms of two scalar potentials and one vectorial potential. Based on Biot’s theory, the frequency domain Holmholtz equations for the potentials are derived through the Fourier transformation with respect to time. The general solutions for the potentials are derived through the Fourier transformation with respect to the horizontal coordinates. Numerical results suggest that moving loads have very complicated effects on the dynamic response of the porous medium. Generally speaking, a moving load with a high speed will generate a larger response in the porous medium than a static or a lower speed load.     

Effect of imperfect bonding on axisymmetric elastodynamic response of a lined circular tunnel in poroelastic soil due to a moving ring load

A two-dimensional elasticity analysis for steady-state axisymmetric dynamic response of an arbitrarily thick elastic homogeneous hollow cylinder of infinite length, which is imperfectly bonded to the surrounding fluid-saturated permeable formation, subject to an axially moving ring load, is presented. The problem solution is derived by using Biot’s dynamic theory of poroelasticity in conjunction with double Fourier transformation with respect to time (frequency) and axial coordinate (axial wave number). The analytical results are illustrated with numerical examples in which a concrete tunnel lining of uniform wall thickness is imperfectly bonded to a surrounding water-saturated poroelastic formation of soft/stiff frame characteristic. Numerical solutions for the radial shell mid-plane and formation displacements are calculated by analytical (numerical) inversion of the Fourier transformation with respect to the frequency (axial wave number). Primary attention is focused on the influence of bonding condition at the liner/soil interface, formation material type, and load velocity on the system’s dynamic response. Limiting cases are considered and good agreements with the solutions available in the literature are obtained.     

移动荷载作用下饱和土地基中的波动特性分析

基于Biot波动方程,经过Fourier变换和逆Fourier变换后可获得波数-频率域以及时间-空间域的解析解。通过数值分析的手段研究了移动荷载作用下饱和多孔弹性地基中波的传播特性。重点就弥散曲线、多谱勒效应、波的成分和动力响应频率等几个特性进行了分析,发现饱和土地基由于比弹性地基多了一项流体介质,波动特性明显差异于弹性介质。     

Dynamic response of an axially loaded rigid sphere embedded in a saturated poroelastic medium

In this paper, an analytical solution for the response of a rigid sphere embedded in a full space poroelastic medium subjected to a dynamic lateral load is derived. The solution is obtained using Biots theory for acoustic waves. In this solution, the displacements of the solid skeleton and the pore pressure are expressed in terms of three scalar potentials. These potentials correspond to the wave velocities of the slow and fast compressional waves and to the shear wave. The governing equation for the dynamic motion is expressed in the frequency domain using Fourier transformation. Different boundary and load conditions were investigated. Curves showing variation in the fluid pressure and solid displacements with the loads frequency were plotted in non-dimensional forms.     

Response of railway track system on poroelastic half-space soil medium subjected to a moving train load

Based on the dynamic poroelastic theory of Biot, dynamic responses of a track system and poroelastic half-space soil medium subjected to moving train passages are investigated by the substructure method. The whole system is divided into two separately formulated substructures, the track and the ground, and the rail is described by introducing the Green function for an infinitely long Euler beam subjected to the action of moving axle loads of the train and the reactions of the sleeper. Sleepers are represented by a continuous mass and the effect of the ballast is considered by introducing the Cosserat model for granular medium. Using the double Fourier transform, the governing equations of motion are then solved analytically in the frequency-wave-number domain. The time domain responses are evaluated by the inverse Fourier transform computation for a certain train speed. Computed results show that the shape of the rail displacements of the elastic and poroelastic soil medium are in good agreement with each other of the low train velocity, but the result of the poroelastic soil medium is significantly different to that of the elastic soil medium for the high train velocity which is higher than Rayleigh-wave speed in the soil. The influence of the soil intrinsic permeability on soil responses is discussed with great care in both time domain and frequency domain. The dynamic responses of the soil medium are considerably affected by the fluid phase as well as the load velocity.     

条形荷载下梯度非均匀非饱和土的动力响应分析

基于多孔介质混合物理论, 建立了梯度非均匀非饱和土地基模型, 研究了条形荷载作用下梯度非均匀非饱和土地基的动力响应问题. 通过傅里叶积分变换和Helmholtz矢量分解原理, 获得频域内非饱和土地基动力响应问题的通解, 结合回传射线矩阵法和边界条件, 求解获得了非均匀非饱和土层中位移、应力以及孔隙压力的计算列式. 假设沿深度方向梯度非均匀非饱和土的物理力学性质按幂函数连续变化, 通过数值傅里叶逆变换得到了非均匀非饱和土地基中的应力、位移以及孔隙压力等物理量的数值解, 分析讨论了土体非均匀性对非饱和土介质动力响应的影响规律. 结果表明: 土体非均匀性显著改变了非饱和土中竖向位移、正应力和孔隙压力在其深度方向上的振动模态, 其中孔隙气压在其深度方向的振动频率随着梯度因子的增加而不断增大, 波峰值不断靠近地表处附近; 竖向位移随着梯度因子的增大不断减小; 正应力和孔隙水压随着梯度因子的增大先增大后减小, 并且土体非均匀程度越高, 正应力与孔隙水压的幅值越大.      

层状横观各向同性饱和土的非轴对称动力响应

通过方位角的Fourier变换,将圆柱坐标系下横观各向同性饱和土的Biot非轴对称波动方 程转化为一组一阶常微分方程组. 然后基于径向Hankel变换,建立问题的状态方程;求解状态方程后,得到传递矩阵. 进而利用传递矩阵,结合饱和层状地基的边界条件、排水条件及层间接触和连续条件,求解 了任意震源力作用下层状横观各向同性饱和地基频域动力响应问题. 时域解可通过频率的Fourier积分得到.     

Dynamic stress intensity factors around a cylindrical crack in an infinite elastic medium subject to impact load

This paper investigates transient stresses around a cylindrical crack in an infinite elastic medium subject to impact loads. Incoming stress waves resulting from the impact load impinge on the crack in a direction perpendicular to the crack axis. In the Laplace transform domain, by means of the Fourier transform technique, the mixed boundary value equations with respect to stresses and displacements were reduced to two sets of dual integral equations. To solve the equations, the differences in the crack surface displacements were expanded in a series of functions that are zero outside the crack. The boundary conditions for the crack were satisfied by means of the Schmidt method. Stress intensity factors were defined in the Laplace transform domain and were numerically inverted to physical space. Numerical calculations were carried out for the dynamic stress intensity factors corresponding to some typical shapes assumed for the cylindrical crack.     

Complex variable function method for the scattering of plane waves by an arbitrary hole in a porous medium

Based on the complex variable function method, a new approach for solving the scattering of plane elastic waves by a hole with an arbitrary configuration embedded in an infinite poroelastic medium is developed in the paper. The poroelastic medium is described by Biot's theory. By introducing three potentials, the governing equations for Biot's theory are reduced to three Helmholtz equations for the three potentials. The series solutions of the Helmholtz equations are obtained by the wave function expansion method. Through the conformal mapping method, the arbitrary hole in the physical plane is mapped into a unit circle in the image plane. Integration of the boundary conditions along the unit circle in the image plane yields the algebraic equations for the coefficients of the series solutions. Numerical solution of the resulting algebraic equations yields the displacements, the stresses and the pore pressure for the porous medium. In order to demonstrate the proposed approach, some numerical results are given in the paper.     

A 2.5-D dynamic model for a saturated porous medium. Part II: Boundary element method

The 3-D boundary integral equation is derived in terms of the reciprocal work theorem and used along with the 2.5-D Green’s function developed in Part I [Lu, J.F., Jeng, D.S., Williams, S., submitted for publication. A 2.5-D dynamic model for a saturated porous medium: Part I. Green’s function. Int. J. Solids Struct.] to develop the 2.5-D boundary integral equation for a saturated porous medium. The 2.5-D boundary integral equations for the wave scattering problem and the moving load problem are established. The Cauchy type singularity of the 2.5-D boundary integral equation is eliminated through introduction of an auxiliary problem and the treatment of the weakly singular kernel is also addressed. Discretisation of the 2.5-D boundary integral equation is achieved using boundary iso-parametric elements. The discrete wavenumber domain solution is obtained via the 2.5-D boundary element method, and the space domain solution is recovered using the inverse Fourier transform. To validate the new methodology, numerical results of this paper are compared with those obtained using an analytical approach; also, some numerical results and corresponding analysis are presented.     

粘性阻尼土中变截面桩的纵向振动特性与应用研究

考虑土体轴对称波动效应,对变截面桩在任意激振力作用下的纵向振动特性进行了研究。假定桩为竖直、弹性、变截面体,土为线性粘弹性体,其材料阻尼为粘性阻尼。利用拉普拉斯变换,将定解问题转化到拉普拉斯域内求解,通过引入势函数并结合阻抗函数的传递性,得到了拉普拉斯域内的桩顶阻抗函数解析解,进而可得到频域内的桩顶阻抗函数和速度导纳的解析解,利用卷积定理和傅里叶逆变换,求得了半正弦脉冲激振力作用下桩顶速度时域响应半解析解。基于所得解对桩的纵向振动特性进行了分析,重点讨论了桩身截面变化情况对速度导纳曲线和反射波曲线的影响,得到了许多重要结论。     

Evaluation of Green’s function for vertical point-load excitation applied to the surface of a layered semi-infinite elastic medium

The topic of this paper is to show that the integrals of infinite extent representing the surface displacements of a layered half-space loaded by a harmonic, vertical point load can be reduced to integrals with finite integration range. The displacements are first expressed through wave potentials and the Hankel integral transform in the radial coordinate is applied to the governing equations and boundary conditions, leading to the solutions in the transformed domain. After the application of the inverse Hankel transform it is shown that the inversion integrands are symmetric/antimetric in the transformation parameter and that this characteristic is preserved for any number of layers. Based on this fact the infinite inversion integrals are reduced to integrals with finite range by choosing the suitable representation of the Bessel function and use of the fundamental rules of contour integration, permitting simpler analytical or numerical evaluation. A numerical example is presented and the results are compared to those obtained by the CLASSI program.     

Three Dimensional Vibration Transmission of a Plate on a Layer of Porous Medium Due to a Moving Load

Three-dimensional (3D) transmission of vibration in an infinite elastic thin plate on a layer of poroviscoelastic medium, due to a harmonic, rectangular moving load, is investigated theoretically based on Biot’s theory. The material of the medium is idealized as a uniform, fully saturated poroviscoelastic layer on bedrock. By introducing four scalar potential functions and Helmholtz decomposition theorem, analytical solutions of stress, displacement, and pore pressure with and without thin plate are derived using Fourier transform technique. Numerical results are obtained with the help of inverse Fourier transform and are used to analyze the influence of load velocity, porosity, permeability, relative stiffness of plate versus ground, and the thickness of plate on the vibration. Furthermore, the results are compared with the available dynamic response results of a non-moving load on a layer of viscoelastic material.     

3D scattering by multiple cylindrical cavities buried in an elastic formation

This paper presents the three-dimensional scattering field obtained when multiple cylindrical circular cavities of infinite length buried in a homogeneous elastic medium, are subjected to dilatational point loads placed at some point in the medium. The solution is formulated using boundary elements for a wide range of frequencies and spatially harmonic line loads, which are then used to obtain time series by means of (fast) inverse Fourier transforms into space-time. The method and the expressions presented are implemented and validated by applying them to a cylindrical circular cavity buried in an infinite homogeneous elastic medium subjected to a dilatational point load, for which the solution is calculated in closed form. The boundary elements method is then used to evaluate the wave-field elicited by a dilatational point load source in the presence of a different number of cylindrical oval cavities. Simulation analyses using this idealized model are then applied to the study of wave propagation patterns in the vicinity of these inclusions. The amplitude of the wavefield in the frequency vs axial-wavenumber domain is presented, allowing the recognition, identification, and physical interpretation of the variation of the wavefield.     

3-D DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED SOILS

Introduction Thedynamicanalysisofsaturatedporoelasticmediaisofbroadmeaninginmanyaspectssuchasseismology,earthquakeengineering,soilmechanicsandgeophysics,etc.Thetransverselyisotropicsaturatedsoilsissolid fluidcoupledtwophasemediainwhichthe soilskeletondisp…     

Flow of Maxwell fluids in porous media

In this paper we analyze the flow of a Maxwell fluid in a rigid porous medium using the method of volume averaging. We first present the local volume averaged momentum equation which contains Darcy-scale elastic effects and undetermined integrals of the spatial deviations of the pressure and velocity. A closure problem is developed in order to determine the spatial deviations and thus obtain a closed form of the momentum equation that contains a time-dependent permeability tensor. To gain some insight into the effects of elasticity on the dynamics of flow in porous media, the entire problem is transformed to the frequency domain through a temporal Fourier transform. This leads to a dynamic generalization of Darcy's law. Analytical results are provided for the case in which the porous medium is modeled as a bundle of capillary tubes, and a scheme is presented to solve the transformed closure problem for a general microstructure.     

3D dynamic stress intensity factors around two parallel square cracks in an infinite elastic medium under impact load

Summary  Transient stresses around two parallel cracks in an infinite elastic medium are investigated in the present paper. The shape of the cracks is assumed to be square. Incoming shock stress waves impinge upon the two cracks normal to tzheir surfaces. The mixed boundary value equations with respect to stresses and displacements are reduced to two sets of dual integral equations in the Laplace transform domain using the Fourier transform technique. These equations are solved by expanding the differences in the crack surface displacements in a double series of a function that is equal to zero outside the cracks. Unknown coefficients in the series are calculated using the Schmidt method. Stress intensity factors defined in the Laplace transform domain are inverted numerically to the physical space. Numerical calculations are carried out for transient dynamic stress intensity factors under the assumption that the shape of the upper crack is identical to that of the lower crack. Received 2 February 2000; accepted for publication 10 May 2000     

Dynamic Green’s Functions for an Anisotropic Multilayered Poroelastic Half-Space

The dynamic responses of an anisotropic multilayered poroelastic half-space to a point load or a fluid source are studied based on Stroh formalism and Fourier transforms. Taking the boundary conditions and the continuity of the materials into consideration, the three-dimensional Green’s functions of generalized concentrated forces (force and fluid source) applied at the free surface, interface and in the interior of a layer are derived in the Fourier transformed domain, respectively. The actual solutions in the frequency domain can further be acquired by inverting the Fourier transform. Finally, numerical examples are carried out to verify the presented theory and discuss the Green’s fields due to three cases of a concentrated force or a fluid source applied at three different locations for an anisotropic multilayered poroelastic half-space.

     

Ground vibration due to a high-speed moving harmonic rectangular load on a poroviscoelastic half-space

The transmission of vibrations in the ground, due to a high-speed moving vertical harmonic rectangular load, is investigated theoretically. The problem is three-dimensional and the interior of the ground is modelled as a totally or partially saturated porous viscoelastic half-space, using the complete Biot theory. The solutions in the transformed domain are obtained using a double Fourier transform on the surface spatial variables. A modified hysteretic damping model defined in the wavenumber domain is used, first presented by Lefeuve-Mesgouez et al. [Lefeuve-Mesgouez, G., Le Houédec, D., Peplow, A.T., 2000. Vibration in the vicinity of a high-speed moving harmonic strip load. Journal of Sound and Vibration 231(5) 1289–1309]. Numerical results for the displacements of the solid and fluid phases, over the surface of the ground and in depth, are presented for loads moving with speeds up to and beyond the Rayleigh wave speed of the medium.     

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.