垂直受荷群桩在半空间饱和土中的稳态反应

研究了均质半空间饱和土中群桩在垂直稳态荷载作用下的动力反应问题。半空间饱和土采用Biot提出的三维波动原理。将桩看作是一维的弹性杆单元。利用圆形简谐载荷作用下的Biot固结方程的基本解和桩土之间的变形协调条件得到饱和土中群桩的第二类Fredholm积分方程,群桩的动力阻抗采用动力相互作用因子的方法。分析桩距等参数对群桩阻抗、轴力和孔压等的影响。本文的方法可以用于计算层状饱和土中群桩的动力反应问题。     

非饱和土-隧道系统动力响应计算的波函数法

针对非饱和地基土中埋置隧道的三维动力响应计算问题, 提出了波函数法.采用无限长的Flügge薄壁圆柱壳模拟圆形隧道衬砌,采用流、固、气组成的三相介质模拟非饱和地基土体.分别采用分离变量法以及Helmholtz矢量分解定理求解薄壁圆柱壳的振动控制方程与非饱和土的波动方程.根据隧-土交界面与地表面处的应力、位移以及孔隙流体压力等边界条件,利用平面波与柱面波的转换性质,实现了隧道内作用单位简谐载荷时隧道衬砌与土体系统动力响应的耦合求解.通过与既有单相弹性介质2.5维有限元-边界元法、两相饱和多孔介质2.5维有限元-边界元法以及三相非饱和介质Pip in Pip半解析法的计算结果进行对比, 验证了本文计算方法的可靠性. 最后,基于该方法, 通过算例分析了不同饱和度下非饱和土-隧道系统的动力响应特征.结果表明, 饱和度对土体动位移与超孔隙水压力的幅值响应有较大影响.该方法的非饱和地基土参数退化后,也可用来计算和分析饱和地基土或单相弹性地基土与隧道系统的动力响应.      

A discrete element model for simulating saturated granular soil

A numerical model is developed to simulate saturated granular soil, based on the discrete element method. Soil particles are represented by Lagrangian discrete elements, and pore fluid, by appropriate discrete elements which represent alternately Lagrangian mass of water and Eulerian volume of space. Macro-scale behavior of the model is verified by simulating undrained biaxial compression tests. Micro-scale behavior is compared to previous literature through pore pressure pattern visualization during shear tests. It is demonstrated that dynamic pore pressure patterns are generated by superposed stress waves. These pore-pressure patterns travel much faster than average drainage rate of the pore fluid and may initiate soil fabric change, ultimately leading to liquefaction in loose sands. Thus, this work demonstrates a tool to roughly link dynamic stress wave patterns to initiation of liquefaction phenomena.     

准饱和土体中圆形衬砌对弹性波的散射

采用Vardoulakis和Beskos提出的准饱和土体的波动控制方程,根据Helmholtz矢量分解定理,得到了准饱和土中P1波(快压缩波)、P2波(慢压缩波)和S波(剪切波)的波数的势函数表达式.将准饱和土体和圆形衬砌视为各向同性的均质体,运用波函数展开法将入射波、散射波和折射波的势函数展开成Fourier-Bessel函数的级数形式,根据准饱和土体与衬砌边界处应力和位移连续及衬砌内完全自由的边界条件,得到了平面P1波入射下,准饱和土体内深埋圆形衬砌的散射系数和折射系数的理论解,通过数值计算分析了饱和度对准饱和土体和衬砌的DSCF(动应力集中因子)及准饱和土体的PPCF(孔压集中因子)的影响规律,结果表明:准饱和土体的DSCF随着饱和度的增大而减小,衬砌的DSCF基本不受饱和度的影响,而准饱和土体的PPCF则随着饱和度的增大而增大.     

弹性波在准饱和土和弹性土界面的反射与透射

考虑土颗粒和孔隙流体的压缩性,根据准饱和土体的变形连续条件和平衡条件,推导了P1波(快压缩波)、P2波(慢压缩波)和S波的计算公式．根据Snell定理,研究了平面P1波从准饱和土体入射到弹性土层时在界面上的反射和透射,并绘制了反射系数、透射系数和界面应力、位移随入射角的变化曲线,结果表明:饱和度是影响界面效应的一个重要因素,当饱和度低于99％时,反射和透射系数以及应力和位移基本相等,变化不大,但当饱和度为100％时,虽然饱和度比99％提高了1％,但反射和透射系数以及应力和位移值比饱和度为99％时明显降低．     

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

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

一维流体饱和粘弹性多孔介质层的动力响应

本文研究了不可压流体饱和粘弹性多孔介质层的一维动力响应问题。基于粘弹性理论和多孔介质理论,在流相和固相微观不可压、固相骨架服从粘弹性积分型本构关系和小变形的假定下,建立了不可压流体饱和粘弹性多孔介质层一维动力响应的数学模型,利用Laplace变换,求得了原初边值问题在变换空间中的解析解,并利用Laplace逆变换的Crump数值反演方法,得到原动力响应问题的数值解。数值研究了饱和标准线性粘弹性多孔介质层的动力响应,分析了固相位移、渗流速度、孔隙压力及固相有效应力等的响应特征。结果表明,与不可压流体饱和弹性多孔介质相同,不可压流体饱和粘弹性多孔介质中亦只存在一个纵波,并且固相骨架的粘性对动力行为有显著的影响。     

表面堆载作用下群桩负摩擦研究

利用Biot固结理论和积分方程方法研究了表面有堆载的群桩负摩擦问题。根据基本解得出了群桩在圆形均布载荷作用下在时间域内的第二类Fredholm积分方程组。运用Laplace变换对上述积分方程组进行简化，求解上述积分方程组并进行相应的数值逆变换就可得出群桩在表面圆形均布载荷作用下的变形、轴力、孔压和桩侧摩阻力随时间的变化情况。     

Blast-induced waves in soil around a tunnel

The paper is concerned with the numerical modelling of waves in saturated granular soil around a cylindrical tunnel with a tubbing-type lining. The waves are induced by a short-term large-amplitude load on the tunnel wall. The soil is assumed to be fully saturated or contain a small amount (a few volume percent) of free gas in the pore water. The behaviour of the solid skeleton is described by a constitutive relation of the hypoplasticity theory. The presence of gas makes the soil stiffness stronger pressure dependent and leads to the formation of a shock front in the soil in the immediate vicinity of the tunnel. Free gas in the pore water is shown to be an adverse factor from the viewpoint of both the residual state of the soil and the deformation of the lining. A permanent reduction of the effective stresses in the soil caused by the wave is stronger for the soil with gas. For a sufficiently large loading amplitude, the wave may lead to the vanishing of the effective stresses and liquefaction of the soil. It is also shown that the presence of gas increases the displacement and the velocity of the tubbings.     

一维渗透力与浮力

为简化分析,针对一维渗流问题,研究提出了土力学中渗透力和浮力的两种推导方法,作为传统的宏观尺度孔隙水隔离体法的有益补充.弹性力学平衡微分方程、土力学Terzaghi有效应力方程和流体力学简化Bernoulli方程构成本文分析渗透力和浮力的3个基本方程.在基本方程基础上,容易推导出相应的骨架和孔隙水两种隔离体的平衡微分方程,从而在静力平衡范畴内揭示渗透力和浮力的内涵.单位体积饱和土体的渗透力,源于总水头压力的梯度,而浮力则源于位置水头压力在竖向的梯度,这两者统一于骨架或孔隙水的平衡微分方程.实际工程关注的有效应力计算问题,一般可以直接应用3个基本方程来确定;只有在简化条件下可按渗透力和浮力计算土体中有效应力分布规律.还讨论了若干研究热点问题,重点探讨了当前一种渗透力新定义j=nγ_wi在形式上的合理性以及在实际应用中可能存在的风险,并验证了一维渗透力的一种经典精细化表述结果中考虑渗流速度时间导数的严谨性.在土力学渗透力和浮力问题研究中应重视和正确应用Terzaghi有效应力方程.     

饱和土中球空腔的瞬态动力响应

采用工程上通用的饱和土力学模型,考虑流体与固全之间的耦合作用,利用Ｌａｐｌａｃｅ变换求解了饱和土中球空腔的瞬态动力响应问题,得到了变换域内的解析解,借助数值Ｌａｐｌａｃｅ换求解了饱和土中球空腔瞬态动力响应的位移、应力及孔压的变化规律,为分析地下结构瞬态动力响应提供了一种有效的方法,模型符合工程实际,有一定的工程应用价值。     

桩侧扰动分区饱和土中桩竖向振动特性研究

考虑沉桩过程中桩周土受到扰动产生软化或硬化的现象,建立了饱和土中径向分区桩竖向振动简化模型,通过分离变量方法,得到了桩竖向振动频域解析解和时域半解析解。分析了分区土对桩竖向振动特性的影响,数值计算结果表明,土层模量比、分区半径和桩底支承系数对桩竖向振动有较明显的影响。研究结果还表明,在较强的振动作用下,由于分区软化或强...     

标准线性固体土模型条件下桩的动力响应分析与应用

在桩侧土为标准线性固体模型条件下，建立了成层土中考虑桩端和桩侧土作用的有限长桩，桩顶受纵向激励情况下的定解问题，用拉氏变化(Laplace)求得桩顶阻抗函数和速度频率响应的函数，然后利用卷积定理和傅立叶(Fourier)逆交换求得任意激振力作用下桩顶速度响应的半解析解，并研究了模型中各个主要参数对桩顶响应及速度响应曲线的影响，并分析了缺陷桩和不同土参数对曲线的影响，得到了许多新的有意义的结论。     

饱和土地基在简谐荷载作用下的解析分析

基于多孔介质混合物理论,用解析的方法研究了不可压饱和土地基受到简谐荷载作用下的动力响应问题。利用Fourier积分变换求解耦合方程组,得到了二维饱和土介质在简谐荷载作用下的通解。针对表面透水的具有下卧基岩的饱和土层以及半无限饱和土地基的边界条件,获得了固体骨架位移、孔隙流体位移、固体骨架有效应力以及孔隙流体压力的积分形式解答,并通过数值算例分析了饱和土地基在简谐荷载作用下的响应。     

黏弹性准饱和土中球空腔动力特性

在频率域内研究了内水压力作用下分数导数型黏弹性准饱和土中球空腔的稳态动力响应. 通 过引入与孔隙率有关的应力系数合理地确定了介质和孔隙水共同承担的内水压力值. 将土骨 架和衬砌分别视为具有分数导数本构模型的黏弹性体和多孔柔性材料, 基于Biot两相介质模 型, 通过引入位移势函数解耦得到了内水压力作用下分数导数型黏弹性准饱和土中半封闭球 形空腔的位移、应力和孔隙水压力解析表达式. 考察了物性和几何各参数对球形空腔动力响 应的影响, 结果表明: 分数导数本构模型更合理地描述了土体的动力学行为; 饱和度对应力 和孔隙水压力影响较大, 而对位移影响较小.     

Compressional waves in fluid-saturated porous solid containing a penny-shaped crack

Diffraction of normal compression waves by a penny-shaped crack in a fluid-saturated porous medium is investigated. Two wave types are considered, namely, compressional wave of the first kind, and the second kind. The former, also known as fast wave, propagates primarily through the solid, whereas the latter or slow wave, propagates mainly in the fluid. Each wave propagates in the medium along with induced wave of the same type in the companion constituent of the material. Application of Biot’s theory in conjunction with integral transform technique reduces the problem to a mixed boundary-value problem whose solution is in turn governed by a Fredholm integral equation of the second kind. Near-field and far-field solutions are obtained in terms of the dynamic stress-intensity factor and the scattering cross section, respectively. They are of particular importance to the linear elastic fracture mechanics (LEFM) and in the scattering theory of elastic waves. The mode I stress-intensity factors are computed numerically for a set of selected material property values, and shown graphically for various mass density and viscosity-to-permeability ratios. The obtained results reveal significant impact of the presence of pore fluid upon the stress-intensity factors, both magnitudes and frequencies at their peak values. The influence of the fluid is also observed from the calculated scattering cross sections of the scattered far-field. Accuracy of the present solution procedure is verified by comparing the numerical results with existing results in the limiting case of dry elastic materials.     

Body waves in fractured porous media saturated by two immiscible newtonian fluids

A study of body waves in fractured porous media saturated by two fluids is presented. We show the existence of four compressional and one rotational waves. The first and third compressional waves are analogous to the fast and slow compressional waves in Biot's theory. The second compressional wave arises because of fractures, whereas the fourth compressional wave is associated with the pressure difference between the fluid phases in the porous blocks. The effects of fractures on the phase velocity and attenuation coefficient of body waves are numerically investigated for a fractured sandstone saturated by air and water phases. All compressional waves except the first compressional wave are diffusive-type waves, i.e., highly attenuated and do not exist at low frequencies.Now at Izmir Institute of Technology, Faculty of Engineering, Gaziosmanpasa Bulvari, No.16, Cankaya, Izmir, Turkey.     

粘弹性饱和土体中半封闭圆形隧洞的动力响应分析

基于Biot波动方程,研究分析了粘弹性饱和土体中半封闭圆形隧洞的动力响应问题．假定衬砌材料为多孔介质,引入了更符合工程实际的半封闭边界条件．通过引入势函数,在Laplace变换域中得到隧洞边界上作用轴对称荷载和流体压力条件下应力、位移和超孔隙水压力的解答．利用Laplace数值逆变换得到时域中的解,分析了隧洞边界的半透水特性对隧洞动力响应问题的影响,结果表明：隧洞边界的半透水特性对应力、位移场的变化和超孔隙水压力的消散有很大的影响,透水和不透水下条件的解仅是本文的两个特例。     

Waves in partially saturated porous media

The propagation of compressional waves in a porous medium is investigated in case the pore liquid contains a small volume fraction of gas. The effect of oscillating gas bubbles is taken into account by introducing a frequency-dependent fluid bulk modulus, which is incorporated in the Biot theory. Using a shock tube technique, new experimental data are obtained for a porous column subjected to a pressure step wave. An oscillatory behaviour is observed, consisting of two distinct frequency bands, which is predicted by the theoretical analysis.     

Coupled hydro-mechanical effects in a poro-hyperelastic material

Fluid-saturated materials are encountered in several areas of engineering and biological applications. Geologic media saturated with water, oil and gas and biological materials such as bone saturated with synovial fluid, soft tissues containing blood and plasma and synthetic materials impregnated with energy absorbing fluids are some examples. In many instances such materials can be examined quite successfully by appeal to classical theories of poroelasticity where the skeletal deformations can be modelled as linear elastic. In the case of soft biological tissues and even highly compressible organic geological materials, the porous skeleton can experience large strains and, unlike rubberlike materials, the fluid plays an important role in maintaining the large strain capability of the material. In some instances, the removal of the fluid can render the geological or biological material void of any hyperelastic effects. While the fluid component can be present at various scales and forms, a useful first approximation would be to treat the material as hyperelastic where the fabric can experience large strains consistent with a hyperelastic material and an independent scalar pressure describes the pore fluid response. The flow of fluid within the porous skeleton is defined by Darcy's law for an isotropic material, which is formulated in terms of the relative velocity between the pore fluid and the porous skeleton. It is assumed that the form of Darcy's law remains unchanged during the large strain behaviour. This approach basically extends Biot's theory of classical poroelasticity to include finite deformations. The developments are used to examine the poro-hyperelastic behaviour of certain one-dimensional problems.