Dynamic response of a non-circular lined tunnel with visco-elastic imperfect interface in the saturated poroelastic medium

A semi-analytical method is proposed to investigate the dynamic response of a non-circular tunnel with visco-elastic imperfect interface in poroelastic medium. Biot’s dynamic theory is used to simulate the saturated poroelastic medium, and the governing equations are solved by reducing them into three Helmholtz equations that the potential functions satisfy. The visco-elastic interface model with elastic and viscosity coefficients is adopted to analyze the interface effect around the non-circular lined tunnel. The analytic solutions of displacements and stresses are expanded in terms of wave functions. Some numerical examples are provided to analyze the effect of visco-elastic interface on the dynamic stress around the tunnel.     

Elastic waves along a cylindrical borehole in a poroelastic medium saturated by two immiscible fluids

The propagation of elastic waves along a cylindrical borehole filled with/without liquid and embedded in an infinite porous medium saturated by two immiscible fluids has been studied. The theory of porous media saturated by two immiscible fluids developed by Tuncay and Corapcioglu (1997) is employed. Frequency equations determining the phase velocity of axial symmetric waves are obtained. It is found that the surface waves along cylindrical borehole are dispersive. The dispersion equation of Rayleigh-type surface waves along the boundary of a poroelastic solid half-space saturated by two immiscible fluids is also obtained. Some special cases have been deduced and the dispersion curves are obtained numerically for a peculiar model. It is found that the density of fluids affects the Rayleigh mode.     

准饱和土中圆柱形衬砌的瞬态动力响应分析

内源瞬态荷载作用下圆柱形孔洞的动力响应解答是土动力学的经典问题之一。已有研究大都假设孔洞周围土体为理想弹性介质或完全饱和多孔介质。然而,实际工程中不存在完全弹性和完全饱和土体。分别视衬砌结构和周围土体为弹性材料和准饱和多孔介质（饱和度 95%）,根据牛顿第二定律、达西定律和Biot波动理论推导出准饱和土体的控制方程。根据边界条件导出衬砌和土体的位移、应力和孔隙压力的Laplace变换空间的解答。利用反Laplace变换数值计算方法,将解答转换为时域解。分析了饱和度对衬砌位移、应力和孔压的影响,结果表明,当95% 99%时,饱和度对径向位移和切向应力的影响较小;99% 100%时,饱和度对径向位移和切向应力的影响较大;但饱和度对孔隙压力的影响远大于对径向位移和切向应力的影响。得出位移、应力和孔压沿径向的衰减规律,当95% 99%时,饱和度对径向位移和切向应力沿径向衰减影响较小,99% 100%时,饱和度对径向位移和切向应力沿径向衰减影响较大,但饱和度对孔压沿径向的衰减影响远大于对径向位移和切向应力沿径向的衰减。     

Scattering of a plane wave by a lined cylindrical cavity in a poroelastic half-plane

The focus of this contribution is to develop a complex variable function method to solve the two-dimensional scattering of plane waves by a lined cylindrical cavity in the poroelastic half-plane. The poroelastic half-plane is based on Biot’s dynamic theory, and the governing equations are solved by reduction to three Helmholtz equations. The lining structure can be treated as an elastic material and decoupled into two Helmholtz equations. Here, the large circle assumption is applied to simulate the half-plane boundary. By using appropriate boundary conditions and continuity conditions, the unknown coefficients in the potentials can be determined. Selected numerical results are presented in this paper.     

Scattering of elastic wave by a partially sealed cylindrical shell embedded in poroelastic medium

Scattering of an elastic wave by a cylindrical shell embedded in poroelastic medium is investigated theoretically with the assumption that the shell material is also a porous elastic medium. The porous medium is modellized via Biot's theory and the scattering by cylindrical shell is expressed by the definition of scattering matrix. The normal mode expansion technique is employed for analysing the scattering field, and the asymptotic solutions of displacements, stresses and pore pressure are derived. Two limiting cases‐scattering by a poroelastic cylinder in Biot medium and a elastic cylindrical shell in elastic medium are obtained from the general solutions. The dispersion curves of displacement amplitude at the interface of shell and medium is compared with the case of elastic shell. The scattering amplitude associated with the fast, slow and transverse waves are identified by numerical simulation. Furthermore, the influence of the poroelastic property of shell material on scattering amplitude is analysed. Copyright © 2005 John Wiley & Sons, Ltd.     

空心圆柱饱和多孔介质热固结问题的解析解

基于考虑热渗效应和等温热流效应的热-水-力耦合的线性热弹性固结控制方程,建立无限长空心圆柱饱和多孔介质热固结问题的一种理论求解方法。该方法先给出Laplace变换域上的解,然后,利用Stehfest法求其数值逆变换。该理论解考虑了空心圆柱体内、外透水界面随时间变化的外力和温度荷载耦合作用过程。最后,通过一算例分析了饱和多孔介质的热固结特征,给出其温度、孔压、位移和应力的演化规律     

衬砌隧道中移动轴向激励作用下两相多孔介质动力响应

为了研究衬砌隧道中地铁振动在饱和地层中的传播情况,采用移动轴向激励模拟技术,建立了隧道-衬砌-两相介质的动力分析模型;利用波函数展开法、傅里叶变换法等,推导了频域内衬砌隧道移动轴向激励作用下两相多孔介质动力响应的解析解,并给出了两相多孔介质临界速度的经验公式;通过离散快速傅里叶逆变换得到了时-空域内两相介质的动力响应。结果表明：无衬砌隧道中移动的轴向常激励作用下,两相介质临界速度只与介质的剪切模量和密度有关,且数值接近其剪切波速的1.1倍;对于衬砌隧道,介质的临界速度随着衬砌剪切模量的增大而增大,随着衬砌密度的增大而减小,衬砌对振动向介质中的传播有一定的削弱作用,衬砌剪切模量与介质剪切模量相差越大,削弱越明显;动力响应频率越接近激振频率,其幅值越大,所对应的临界速度越小。     

饱和土半空间中圆柱形孔洞对平面P波的散射

在 Biot饱和多孔介质动力学理论的基础上,首次建立了求解饱和土半空间中圆柱形孔洞对平面 P 波散射问题的波函数展开法。首先。分析了具有圆柱形孔洞的饱和土半空间场地在平面 P 波入射下产生散射波系,并将入射波和散射波的波函数在圆柱坐标下展开。然后,引入边界条件,求出散射波函数的待定系数,从而,得到饱和土半空间中圆柱形孔洞对平面P波的散射问题的解析解。根据所得的波函数的解,可求解区域内的位移、应力的值,同时,分析了入射波频率、入射角对柱面上的应力集中因子的影响。     

Analytical solutions of a finite two‐dimensional fluid‐saturated poroelastic medium with compressible constituents

An analytical solution is proposed for transient flow and deformation coupling of a fluid‐saturated poroelastic medium within a finite two‐dimensional (2‐D) rectangular domain. In this study, the porous medium is assumed to be isotropic, homogeneous, and compressible. In addition, the point sink can be located at an arbitrary position in the porous medium. The fluid–solid interaction in porous media is governed by the general Biot's consolidation theory. The method of integral transforms is applied in the analytical formulation of closed‐form solutions. The proposed analytical solution is then verified against both exact and numerical results. The analytical solution is first simplified and validated by comparison with an existing exact solution for the uncoupled problem. Then, a case study for pumping from a confined aquifer is performed. The consistency between the numerical solution and the analytical solution confirms the accuracy and reliability of the analytical solution presented in this paper. The proposed analytical solution can help us to obtain in‐depth insights into time‐dependent mechanical behavior due to fluid withdrawal within finite 2‐D porous media. Moreover, it can also be of great significance to calibrate numerical solutions in plane strain poroelasticity and to formulate relevant industry norms and standards. Copyright © 2014 John Wiley & Sons, Ltd.     

Analytical solution of undrained cylindrical cavity expansion in saturated soil under anisotropic initial stress

This paper presents an analytical solution for undrained elasto-plastic cylindrical cavity expansion in saturated soil under anisotropic initial stress. The problem is formulated by assuming small-strain deformation in the elastic zone and large-strain deformation in the plastic zone. Plastic yielding is determined by the Tresca failure criterion and an associated flow rule. Two stress functions are used to describe the stress state in the two zones around the cavity. The elasto-plastic boundary can subsequently be determined by solving the two stress functions under the stress boundary conditions. Additionally, the cavity pressure-expansion relationship, the total stress and the excess pore pressure around the cavity wall under anisotropic initial stress can be easily obtained by application of a numerical integration. The results show that the cavity pressure and excess pore pressure under the isotropic initial total stress (K = 1) are larger than those under the anisotropic initial stress (K > 1 and K < 1), which is true at all states of the expansion. The higher value of K develops lower stress and pore pressure around the cavity wall at the ultimate states. However, the stress and excess pore pressure are not sensitive to the value of K. The present solution may be used for analyzing the uplift capacity of plate anchors in soils and Horizontal Directional Drilling (HDD) problems such as the tunneling, and pipeline installation.     

Analytical study of freezing behavior of a cavity in thermo-poro-elastic medium

The goal of this paper is to present an analytical solution to have a first insight of the impact of ice formation on the surrounding porous rock on underground cavities like reservoir, pipes, tunnels or wellbores. Among the other analytical solutions found in the literature on this topic, the originality of this work resides in the rigorous theoretical framework of poromechanics, which considers the coupling between liquid water and ice crystal under thermodynamic equilibrium. Liquid water transport, thermal conduction, and elastic properties of the phases are also considered. Two analytical solutions are presented, based on a linearization of the system of governing equations. The first one deals with a spherical cavity within an infinite porous medium leading to an exact analytical solution. It allows validating the Stehfest’s algorithm on the numerical inversion of Laplace Transform, used in the second analytical solution, which considers a cylindrical excavation. The validity of this solution is assessed by comparing its results to that issued from a numerical resolution of the nonlinear system of equations. The analytical solution is then ultimately used to identify the influence of key parameters like the thermal/hydraulic conductivities, the amount of ice formed and the thermal dilatation coefficients on the mechanical response of a cylindrical cavity submitted to an internal frost.     

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.     

分数导数模型描述的饱和土桩纵向振动分析

土体具有黏弹性性质,为了更好地考虑饱和土体固相土骨架的黏弹性性质,在分数导数理论和多孔介质理论的基础上,将土体视为液固饱和两相介质,并利用分数导数模型来描述饱和土固相土骨架的应力-应变关系,建立分数导数模型描述的饱和土的控制方程。在三维轴对称情况下,利用势函数和分离变量法研究了分数导数模型描述的饱和土中桩的振动问题。分析了模型参数对饱和土中桩的竖向振动的影响。研究表明,分数导数模型描述的饱和土的控制方程应用范围更广,模型参数对桩的竖向振动有较大的影响。     

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.     

铁路碎石道碴层通风阻力参数试验研究

铁路碎石道碴层能被看作高渗透率的多孔介质,由于其具有较好的通风特性,它的存在改变了寒区铁路路基的顶部热交换条件,对路基温度场产生不可忽视的影响.为了能够实现对寒区铁路路基温度场的更合理预测,基于对铁路碎石道碴层传热和通风特性的研究现状,通过风洞试验对铁路碎石道碴层的通风特性及通风阻力参数进行了试验研究.试验中选用的铁路碎石道碴平均粒径为4.2cm,孔隙率为0.433,对不同空洞风速条件下铁路碎石道碴层中的压力梯度和渗流速度进行了测试分析.结果表明:铁路碎石道碴层中空气渗流速度随着压力梯度的增大而增大,其内部压力梯度与渗流速度之间存在良好的二次非线件关系.通过计算得到铁路碎石道碴层的渗透率和惯性阻力系数,为寒区铁路路基热稳定性的数值研究提供了参数依据.     

低频谐变电偶极源激发的地下电磁场的闭式解

研究偶极子源产生的电磁场在地下或水下无线电通信和地球物理勘探领域有着重要的应用价值。偶极子和介质的作用往往导致Sommerfeld型积分的解。为了深入揭示电磁波在地下的分布与传播规律,本文针对有耗媒质中为获得大探测深度和长距离通信使用的低频发射天线,在准静态近似条件下导出了水平电偶极子产生的下电磁场Sommerfeld型积分解的闭式表达式。这些公式适用于任意场区和除奇异点以外的任何场点,与数值积分相比,可显著提高计算速度。利用闭合表达式,本文还提出了一种新的场区划分方法。