用薄层法分析层状地基中各种基础的阻抗函数

薄层法是分析和模拟弹性波在层状介质中传播的一种半解析半数值方法.采用薄层单元和傍轴边界分别模拟层状地基和弹性半空间.利用薄层法位移基本解和容积法推导了层状地基中基础-地基动力相互作用方程及块式基础、桩基础和承台群桩基础阻抗函数的统一计算公式.通过计算半无限弹性地基中桩基础、块式基础和二层地基上基础的阻抗函数验证了方法的适用性.进而计算了某实际层状地基中承台群桩基础的阻抗函数,并与试验结果进行对比,两者吻合较好.本文方法可用于分析弹性层状地基中各种基础的阻抗函数.     

精度可控地基阻抗力的一种时域差分计算方法

将集总参数模型与时域差分递归的滤波器模型相结合,提出了一种可完整考虑基础阻抗函数中奇异项和正则项且精度可控的时域地基阻抗力的计算方法,数值算例表明该方法在感兴趣的频域范围内对频率相关的阻抗函数拟合良好.      

用薄层法分析层状地基中条形基础的阻抗函数

薄层法是分析和模拟弹性波在层状介质中传播的一种半解析半数值方法,以往薄层法在柱坐标体系中建立求解方程,并通过直角坐标系和柱坐标系的转化关系而得到直角坐标系中的解答,本文从直角坐标直接推导了层状地基在无限长线荷载作用下的计算公式,求解了层状地基在垂直方向和水平方向上无限长线激振荷载作用下薄层法位移基本解.结合容积法得出了层状地基中基础-地基动力相互作用方程及条形基础阻抗函数的计算公式.本文计算了半无限弹性地基以及基岩上覆盖层在无限长简谐线荷载作用下的位移反应,计算了半无限弹性地基以及基岩上覆盖层地基中明置与埋置条形基础的地基阻抗函数.计算结果与已有的研究结果的比较表明两者吻合较好,验证了本文方法的适用性.     

基于弹性半空间Green函数的明置条形基础阻抗解答

基于弹性半空间理论,在Lamb解的基础上推导了非轴对称的简谐线激振荷载及条形均布激振荷载作用下的地基Green函数。然后将明置条形基础与地基的接触面分割成若干个条形单元,根据刚体位移所决定的各单元位移,运用上述Green函数求解各单元的接触力,将所求得的力叠加即为地基阻抗函数。通过分段积分法及Cauchy主值积分对涉及的多值广义积分进行数值处理,计算了明置条形基础的垂直阻抗,观察了其收敛性,亦计算了其水平和回转动柔度,并与半解析半数值的薄层法的结果进行了比较,表现出很好的一致性。最后通过算例讨论了土体泊松比对明置条形基础阻抗的影响。     

层状地基中埋管地基阻抗函数的分析方法

本文首先利用波函数展开法推导了弹性半空间中埋管的地基轴向阻抗函数计算公式.然后采用层状地基中无限长线激振荷载薄层法基本位移解,结合容积法求解了层状地基中埋管的地基轴向、垂直及水平阻抗函数.分别用两种方法计算了弹性全空间、弹性半空间地基内埋管的地基轴向阻抗函数,两者的计算结果符合良好,验证了用薄层法求解的可行性.本文还利用薄层法分析了弹性半空间内埋管的埋深对地基水平、垂直及轴向阻抗函数的影响;计算了弹性半空间内埋管的轴向刚度系数并与我国及日本的相关规范对比分析;计算分析了上海典型层状地基内地铁隧道的地基阻抗函数.     

一种三维饱和土-基础-结构动力相互作用分析方法

地震波入射时土与结构动力相互作用分析是地震工程领域的重要问题.由于问题的复杂性,以往的研究大多考虑地基土为干土情形.而实际工程中,土体中经常充满孔隙水,土层往往是含水层或部分含水层.孔隙水对土层的地震反应影响较大,进而影响支撑于其上的基础和上部结构的响应.因此,有必要考虑饱和土-基础-结构动力相互作用问题.基于此,土体采用Biot模型,利用集中质量显式有限元方法并结合多次透射人工边界进行模拟;结构经有限元离散后,采用Newmark隐式时步积分方法进行分析,可通过ANSYS等有限元软件实现;基础假定为刚性,采用显式时步积分进行求解;土体和结构(及基础)可分别采用不同的时间步距;通过FORTRAN程序实现了三维饱和土-基础-结构系统在地震作用下的整体分析.从饱和多孔介质动力微分方程出发可知,干土是饱和土的特殊情形,通过将流体体积模量及孔隙率置为零,可将饱和土退化到干土,从而将干土统一到饱和土的计算框架中,通过土-结构相互作用算例对此进行了验证,进一步实现了干土与饱和土混合情形时的土-结构动力相互作用分析,使得问题的模拟更贴近实际工程(如考虑地下水位情形).通过算例对比了饱和土地基、干土地基、干土覆盖饱和土地基(考虑地下水位)情形时,土-结构相互作用对基础和结构响应的影响,结果表明地下水位对基础和结构响应的影响较大.     

层状地基任意形状刚性基础动力响应求解

提出了基于积分变换、对偶方程与精细积分算法求解多层地基任意形状刚性基础的动力刚度问题. 首先在频率波数域内圆柱坐标体系中利用圆形微元的对称与反对称特性建立多层地基中格林影响函数的波动方程,然后将应力和位移关系表示成对偶形式进行精细积分求解以提高计算精度和稳定性. 再将任意形状刚性基础与地基的交界面离散化为一系列圆形微元,利用格林影响函数建立其平动与转动动力刚度的矩阵方程. 该求解方法高效、准确并且计算稳定,适于任意复杂多层地基任意形状基础动力刚度的计算.      

自然单元法在Winkler地基薄板计算中的应用

自然单元法是一种基于Voronoi图及Delaunay三角形剖分图,以自然邻接点插值为试函数的一种无网格数值方法.本文以目前该方法中自然邻接点的Laplace插值形函数为基础,求出了其一阶及二阶导函数,建立了Winkler地基上薄板弯曲挠度的自然单元法求解控制方程,并编制了相应的计算程序,通过算例分析表明了本文方法的可行性和有效性.     

用粘性边界有限元法分析弹性半无限地基中的动力问题

本文使用设置粘性边界单元的有限元方法，分析了简谐集中激振力产生的地表位移反应，刚性基础及桩基础的阻抗函数等半无限地基中的动力问题。计算结果与其他数值分析方法结果的比较表明，粘性边界单元的有限元方法适用于分析弹性半无限地基中的动力问题。本文还讨论了有限元网格尺寸及模型大小对计算结果的影响。     

广义有限差分法在含阻抗边界空腔声学分析中的应用

声学分析在噪声控制、室内隔音等工程计算中有着重要的作用. 由于现实生活中的声学模型往往伴随着吸声材料, 因此分析含阻抗边界条件的声学问题显得十分必要. 广义有限差分法是一种新型区域型无网格数值离散方法, 该方法基于多元函数泰勒级数展开式和加权最小二乘拟合, 将控制方程中未知参量的各阶偏导数表示为相邻节点函数值的线性组合. 本文首次将广义有限差分法应用于含阻抗边界条件空腔声学问题的分析中, 建立了空腔声场问题的广义有限差分法数值离散格式. 与传统算法相比, 所建立的数值模型具有无需网格剖分和数值积分、计算精度高、适用于大规模声学分析等优点. 通过具有解析解的经典算例, 研究了总节点数目和局部支撑点数目对数值结果的影响, 得到了最大计算频率与节点间距之间关系的经验公式. 此外, 将广义有限差分法应用于无解析解的二维和三维复杂声学模型, 并与COMSOL Multiphysics软件所得的有限元结果进行了比较分析. 数值实验表明, 该算法是一种高效、精确、稳定、收敛的数值模拟方法, 在含阻抗边界空腔声学分析中具有广阔的应用前景.      

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

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

TORSIONAL VIBRATIONS OF RIGID CIRCULAR PLATE ON TRANSVERSELY ISOTROPIC SATURATED SOIL

An analytical method was presented for the torsional vibrations of a rigid disk resting on transversely isotropic saturated soil. By Hankel transform, the dynamic governing differential equations for transversely isotropic saturated poroelastic medium were solved. Considering the mixed boundary-value conditions, the dual integral equations of torsional vibrations of a rigid circular plate resting on transversely isotropic saturated soil were established. By appropriate transform, the dual integral equations were converted into a Fredholm integral equation of the second kind. Subsequently, the dynamic compliance coefficient, the torsional angular amplitude of the foundation and the contact shear stress were expressed explicitly. Selected examples were presented to analyse the influence of saturated soil's anisotropy on the foundation's vibrations.     

3-D DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED SOILS

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

饱和土中端承桩纵向振动特性研究

基于饱和多孔介质理论研究了三维轴对称条件下端承桩在饱和土中的纵向耦合振动. 首先通 过引入势函数对Biot动力固结方程解耦，采用算子分解方法及分离变量法求得饱和土层振动 解，进而利用桩土完全耦合条件得到桩土系统定解. 然后对桩顶的频率和时域响应进行参数 研究，结果表明，桩的长径比和桩土模量比对桩顶动力响应有较大影响，而渗透力对其影响 较小. 最后将幅频及时域反射波理论拟合结果与桩基实测结果加以对比，结果表明理论曲线 与实测曲线规律一致.     

Lamb's integral formulas of two-phase saturated medium for soil dynamic with drainage

When dynamic force is applied to a saturated porous soil,drainage is common.In this paper,the saturated porous soil with a two-phase saturated medium is simulated,and Lamb's integral formulas with drainage and stress formulas for a two-phase saturated medium are given based on Biot's equation and Betti's theorem(the reciprocal theorem).According to the basic solution to Biot's equation,Green's function Gij and three terms of Green's function G4i,Gi4,and G44 of a two-phase saturated medium subject to a concentrated force on a spherical coordinate are presented.The displacement field with drainage,the magnitude of drainage,and the pore pressure of the center explosion source are obtained in computation.The results of the classical Sharpe's solutions and the solutions of the two-phase saturated medium that decays to a single-phase medium are compared.Good agreement is observed.     

The Non-Axisymmetrical Dynamic Response of Transversely Isotropic Saturated Poroelastic Media

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｄｙｎａｍｉｃｒｅｓｐｏｎｓｅｏｆｓａｔｕｒａｔｅｄｐｏｒｏｕｓｍｅｄｉｕｍｉｓｔｈｅｂａｓｉｃｉｓｓｕｅ,ｗｈｉｃｈｉｓａｐｐｌｉｃａｂｌｅｉｎｇｅｏｐｈｙｓｉｃａｌｅｘｐｌｏｒａｔｉｏｎ ,ｓｅｉｓｍｏｌｏｇｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ ,ｏｉｌｅｘｐｌｏｒａｔｉｏｎ ,ｄｉｓｐｏｓａｌｏｆｎｕｃｌｅａｒｗａｓｔｅｓ,ｈｙｄｒａｕｌｉｃｅｎｇｉｎｅｅｒｉｎｇａｎｄｃｉｖｉｌｅｎｇｉｎｅｅｒｉｎｇ .Ｓｉｎｃｅｔｈｅｃｏｎｓｔｉｔｕｔｉｖｅｒｅｌａｔｉｏｎｓａｎｄｗａｖｅ…     

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.     

Acoustical characterization of fluid-saturated porous media with local heterogeneities: Theory and application

A linear dynamic model of fully saturated porous media with local (either microscopic or mesoscopic) heterogeneities is developed within the context of Biot’s theory of poroelasticity. Viscoporoelastic behavior associated with local fluid flow is characterized by the notion of the dynamic compatibility condition on the interface between the solid and the fluid. Complex, frequency-dependent material parameters characterizing the viscoporoelasticity are derived. The complex properties can be obtained through determining the quasi-static poroelastic parameters, the properties of individual constituents, and the relaxation time of the dynamic compatibility condition on the interface. Relationships among various quasi-static poroelastic parameters are developed. It is shown that local fluid flow mechanism is significant only in the porous media with local heterogeneities. The relaxation time of the compatibility condition on the interface depends upon the details of local structure of porous media that control local fluid pressure diffusion. The new model is used to describe the velocity dispersion and attenuation in fully saturated porous media. The proposed model provides a theoretical framework to simulate the acoustical behavior of fully saturated porous media over a wide range of frequencies without making any explicit assumption about the structure of local heterogeneities.     

Application of results on Eshelby tensor to the determination of effective poroelastic properties of anisotropic rocks-like composites

The present work is devoted to the determination of the macroscopic poroelastic properties of anisotropic elastic porous materials saturated by a fluid under pressure. It makes use of the theoretical results provided by Withers [Withers, P.J., 1989. The determination of the elastic field of an ellipsoidal inclusion in a transversely isotropic medium, and its relevance to composite materials. Philosophical Magazine A 59 (4), 759–781.] for the problem of an ellipsoidal inclusion embedded in a transversely isotropic elastic medium. The particular case of a spherical inclusion is very important for rock-like composites such as argillite and shales. The implementation of these results in a micromechanical theory of poroelasticity allows to quantify the effects of the solid matrix anisotropy and of pore space on the effective poromechanical properties. Closed form expressions of Biot tensor and of Biot modulus are presented as well as numerical applications for anisotropic shales.