基于CFD的地震液化研究新进展

综述了近年来关于液化土体流体动力学特征的试验发展状况,以及基于计算流体动力学（简称CFD）的地震液化数值模拟现状,重点介绍了目前比较活跃的可以较高精度模拟液化土体流动状态的三次伪质点数值方法（简称CIP法）。通过对CFD和传统固体力学在地震液化研究中的应用比较,指出了应用CFD的三大优势,即土体大变形问题、液化土体参数分析以及液化土体中结构物的变形应用CFD分析,均可获得较好的结果。进一步提出,在地震液化应用中,未来CFD的发展应该考虑整合液化前的土体性状研究和地震液化中桩-土-结构物的综合分析。     

液化地基中桩的破坏机理研究进展

地震诱发的地基液化对桩基础的破坏极大,液化地基中桩的破坏机理是岩土地震工程中的一个重要研究课题。目前地基液化时桩土结构系统的地震性态尚没有认识充分,已有的研究内容较多局限于桩身材料的强度破坏方面,难以考虑液化土体侧向流动、基桩屈曲失稳、以及土与结构动力相互作用等复杂因素的影响。本文重点加强以下3个方面的深入探讨和研究:(1)液化地基中桩的屈曲失稳;(2)液化地基中桩基破坏的数值模拟新方法;(3)液化地基中桩-土-结构的动力相互作用分析。     

液化流动条件下隧道结构动力分布场研究

将液化流动的土体视为流体,运用流体力学原理,基于矢量符号运算法进行液化场地的动力场分析,求解出动力场解析解。采用ABAQUS/CFD进行液化场地流体动力学有限元分析,根据流动基本控制方程,计算出动力场的数值解。分析结果表明,液化土体横向流动时,隧道结构周围的应力场既包括由表面压力引起的压力阻力,也包括由剪应力引起的摩擦阻力;隧道结构周围的液化土体流动速度很小,但在隧道结构下方存在流动速度加强区;理论计算的动力场解析解大于有限元计算的数值解,但总体的分布变化趋势基本一致;隧道结构附近的应力场变化较为密集,所受到的应力主要分布在液化土体流动的迎面位置。     

液化砂土流动效应的振动台试验研究

采用振动台激励使饱和砂土发生液化,并侧向拖拽埋入砂土中的铝管,模拟液化土体与管体相对运动以分析液化砂土流动的力学效应。引入流体力学理论与方法,推导出以拉力反算表观黏度的表达式以及液化土体作为流体对管壁作用的黏滞剪切力。分析和比较了振时拖动、振后拖动下土体的流体性质及其流动效应的率相关性和孔压相关性,探讨了砂土密实度对土体流动效应的影响。结果表明,土体初始密实度与液化后土体的表观黏度正相关;液化土体的表观黏度以及因流动产生的黏滞剪切力与孔压反相关;液化砂土流动产生的黏滞剪切力具有强烈的率相关性。针对可液化场地中的结构抗震分析,应考虑土-结构率相关相互作用。     

矩形位移模式下被动桩砂土绕流的模型试验研究

被动桩桩周土体易发生绕桩流动,但绕桩流动土体的位移特征及其影响因素迄今还缺乏深入研究。采用自主研制的模型试验装置,对矩形位移模式下砂土绕圆形被动桩流动过程进行详细测试与分析,研究了桩径、加载速度及土体密度对砂土的特征位移和桩身受力的影响。试验结果表明,砂土绕桩流动过程可分为迎土侧土体隆起、桩周土体绕流和背土侧土体沉陷3个阶段。桩径或土体密度越大,桩周砂土越不易发生绕桩流动;桩周土体移动速度越大,则越易发生绕桩流动。分析结果进一步表明,土体滑动面的位置可根据桩身绕流阻力曲线的拐点来判断,在土体滑动面以上,模型桩绕流阻力沿桩长沿深度大致呈线性增长,其值始终小于沈珠江绕流阻力理论值。     

考虑黏性的液化土中水平振动桩基桩顶阻抗研究

建立了三维黏性流体-桩-土体相互作用分析模型,对简谐激振水平动荷载作用下的液化土中桩基振动响应问题进行解析研究。将桩周液化土体视为黏性不可压缩流体,建立流体运动方程,利用亥姆霍兹分解和分离变量法并结合流体边界条件和桩-流体位移、速度连续条件及桩身边界条件,求得了黏性流体动压力及流体速度势解析表达式,从而得到桩身阻力表达式。用饱和多孔介质模型模拟饱和未液化土层,在已有饱和未液化土层振动响应解析解的基础上,推导得出上覆黏性液化流体,下层土体为饱和未液化土中水平振动桩基桩顶阻抗解析解。与已有的水中悬臂梁自由振动解析解对比,验证提出的模型解的正确性,最后分析了流体黏滞系数、桩长、桩土模量比对桩顶阻抗的影响。结果表明,忽略液化土体的黏性特征会高估桩基础桩顶的刚度阻抗,低估其阻尼阻抗。     

抗液化排水刚性桩沉桩过程的土压力响应

抗液化排水刚性桩是一种将刚性桩与竖向排水体相结合的新桩型。基于某建筑桩基工程,开展了抗液化排水刚性桩和不含排水体的普通刚性桩的沉桩对比现场试验,采用了动态土压力传感器实时监测沉桩过程中桩周土体内产生的土压力响应,对比了排水桩与普通桩沉桩对桩周土体水平方向应力及有效应力影响的差异。试验结果表明：抗液化排水刚性桩能够有效减小沉桩过程对桩周深部可液化土体的扰动,在桩身近侧（距桩心0.6 m）深部埋深（-15 m）位置,排水桩的水平土压力响应峰值仅为普通桩的1/4;排水桩能够有效降低沉桩对可液化土层有效应力的影响,使桩周土体更加稳定;在单次沉桩过程中,对于浅部埋深（-5 m）,排水桩对桩周土压力峰值的影响作用较小,对于存在可液化土层的深部埋深（-10、-15 m）,排水桩对土压力峰值的有效影响半径可达4倍桩径。现场试验数据为抗液化排水刚性桩的桩间距选择提供了有力的设计参考依据。     

强度弱化条件下饱和砂土地基中桩-土相互作用p-y曲线研究

在地震荷载作用下,可液化土层中的桩基础往往会由于地基土体液化而发生破坏。在此过程中即使土体最终没有达到完全液化,但由于超孔隙水压力的存在,饱和砂土会发生强度弱化,也会导致土体对桩身水平抗力的降低。此时若不考虑超孔隙水压力对土抗力的影响,仍然采用API规范中的p-y曲线对桩基础进行设计,结果将偏于危险。针对这种情况,首先利用竖向-扭转双向耦合剪切仪对饱和砂土进行了循环扭剪动强度试验,研究了不同弱化状态下饱和砂土的动力特性和弱化参数;然后基于浅层处改进的土楔体理论模型推导极限土抗力公式,并结合深层处的绕桩流动破坏理论模型,得到了任意深度处不同孔压比下的极限土抗力,进而构造了不同弱化状态下饱和砂土地基中桩-土相互作用的p-y曲线。通过研究发现:表征土体强度弱化状态的孔压比对桩-土相互作用的极限土抗力的影响非常显著,孔压比越大,土体强度弱化程度越严重,饱和砂土的极限土抗力值就越小,即横向受荷桩对周围土体的作用随着土体强度弱化程度的增加而降低,反之则增大。     

排水刚性桩群桩抗液化性能的振动台试验研究

排水刚性桩是一种将竖向排水体与刚性桩相结合的新型抗液化处理措施。为研究排水刚性桩群桩的抗液化作用效果,开展了桩顶竖向荷载作用下排水刚性桩处理可液化地基的振动台试验研究,分析了地基土体的超孔压响应、加速度响应及桩顶结构的水平位移响应,并与未设置排水体的普通桩群桩工况进行对比。结果表明：加载开始后,排水桩桩身排水通道有大量超孔隙水排出,普通桩桩身没有排水现象。采用排水桩时超孔压比峰值比普通桩中减小12%,孔压消散稳定后超孔压比减小13%左右,排水桩桩身的排水通道对超孔压的消散作用集中在振动作用的前期。排水桩桩顶的侧向永久位移较普通桩桩顶侧向永久位移减小约27%。试样土体液化前,剪应力-应变滞回圈包络面积较大,土体呈现一定的剪胀特性。液化后,排水桩的剪应力-应变滞回圈的割线模量更大。上述试验结果均表明了排水刚性桩在抗液化方面的有效作用。     

可液化场地碎石桩复合地基地震动力响应分析

采用砂土液化大变形弹塑性本构模型分析可液化砂土,采用模量随应力与应变变化的等效非线性模型增量形式分析碎石桩,应用FLAC3D有限差分软件对地震动力作用下可液化场地碎石桩复合地基进行三维动力响应分析。模拟分析了在地震作用下碎石桩刚度效应和排水效应对加固处理可液化场地的抗液化效果,从初始小变形到液化后大变形的变形发展,超静孔压累积与消散,及桩与土的变形与应力分配变化等。结果表明,所用模型与方法可合理描述可液化场地碎石桩复合地基在地震作用下场地的动力响应特性和抗液化效果;在地震作用下可液化场地中桩周土体与碎石桩体的竖向应力与水平向剪切应力向碎石桩体集中,竖向有效应力比可降至约1/6～1/3;桩周土体与桩体为非协调变形,剪应变比可达7～10;碎石桩抗液化影响范围约为2.5～3倍桩径,对超过3.5倍桩径范围影响较小;碎石桩与砂土渗透系数比大于100时对降低砂土中超静孔隙水压影响明显;碎石桩对场地的加密效应可显著降低超静孔隙水压力,而碎石桩刚度则对超静孔隙水压力变动影响较小,但有助于减低地面加速度响应峰值。     

横向扩展液化土中单桩-土系统的非线性分析

将地震液化场地分为地表的上覆未液化土层、底部的未液化基层以及夹在两者之间的液化土层,基于桩-土相互作用的非线性Winkler模型,考虑桩弯曲的非线性弯矩-曲率本构关系和桩的几何非线性变形,建立了液化土层横向扩展下桩非线性大挠度变形的基本控制方程,并利用打靶法进行了数值求解。同时,给出了桩线弹性小变形情形下的解析解。通过与非线性有限元解和线弹性小变形解析解的比较,验证了文中打靶法的有效性和可靠性。用数值方法分析了液化土层横向扩展对桩力学性能的影响,结果表明：非线性桩-土相互作用和桩材料非线性效应强于桩的几何非线性效应,随着液化土层横向扩展位移的增加,几何非线性效应逐渐增大,此时,应采用完全非线性模型进行桩力学行为的分析。     

液化地基侧向流动引起的桩基础破坏分析

20世纪60年代以来,流动地基中的桩基础的震害现象和抗震设计受到了工程师和研究者的广泛关注。对侧向流动地基中桩基础的一些典型震害现象和其可能震害原因的归纳和分析,显示目前研究仍不能完全解释侧向流动地基中桩基础的震害现象。选取新泻地震中昭和大桥桩基础破坏案例采用p-y方法进行计算,分析地基侧向流动引起的桩基础破坏的影响因素。计算结果显示,合理地描述液化砂土的p-y曲线模型在侧向流动地基桩基础分析中起到关键作用。对于侧向流动地基中桩基础的震害机制的进一步理解和抗震设计,有赖于更为合理和有效地液化砂土中的p-y模型的发展。     

Minimum principle and related numerical scheme for simulating initial flow and subsequent propagation of liquefied ground

The problem of predicting the evolution of liquefied ground, modelled as a viscoplastic material, is addressed by combining a minimum principle for the velocity field, which characterizes such an evolution, and a time step integration procedure. Two different numerical schemes are then presented for the finite element implementation of this minimum principle, namely, the regularization technique and the decomposition‐co‐ordination method by augmented Lagrangian. The second method, which proves more accurate and efficient than the first, is finally applied to simulate the incipient flow failure and subsequent spreading of a liquefied soil embankment subject to gravity. The strong influence of liquefied soil residual shear strength on reducing the maximum amplitude of the ground displacement is particularly emphasized in such an analysis. Copyright © 2005 John Wiley & Sons, Ltd.     

文克尔地基模型下液化土桩基水平振动响应分析

在桩基顶部承受竖向荷载作用的条件下,将完全液化后的上层土体视为流体,将桩基等效为欧拉-伯努利梁模型,探讨了桩底嵌固时桩基顶部的水平振动阻抗。运用流体动力方程模拟顶部液化土层的运动,运用文克尔地基模拟下部非液化分层土的运动。结合传递矩阵法,利用液化土与非液化分层土交界面处的位移、转角和内力连续条件,得到桩基顶部和底部的相关位移?内力表达关系式。根据桩基底部的嵌固条件,求得桩顶阻抗的表达式。与已有文献解进行对比,验证了分析过程的正确性。对阻抗影响条件进行参数分析,表明液化深度、轴力和流体密度大小对桩顶阻抗有不同的影响。     

Test Investigation on Liquefied Deformation Structure in Saturated Lime-Mud Composites Triggered by Strong Earthquakes

Three identical model boxes were made from transparent plexiglass and angle iron. Using the method of sinking water and according to the sedimentary rhythm of saturated calcium carbonate(lime-mud) intercalated with cohesive soil,calcites with particle sizes diameters of ≤ 5 μm,10–15 μm and 23–30 μm as well as cohesive soil were sunk alternatively in water of three boxes to build three test models,each of which has a specific size of calcite. Pore water pressure gauges were buried in lime-mud layers at different depths in each model,and connected with a computer system to collect pore water pressures. By means of soil tests,physical property parameters and plasticity indices(Ip) were obtained for various grain-sized saturated lime-muds. The lime-muds with Ip ranging from 6.3 to 8.5(lower than 10) are similar to liquid saturated silt in the physical nature,indicating that saturated silt can be liquefied once induced by a strong earthquake. One model cart was pushed quickly along the length direction of the model so that its rigid wheels collided violently with the stone stair,thus generating an artificial earthquake with seismic wave magnitude greater than VI degree. When unidirectional cyclic seismic load of horizontal compression-tension-shear was imposed on the soil layers in the model,enough great pore water pressure has been accumulated within pores of lime-mud,resulting in liquefaction of lime-mud layers. Meanwhile,micro-fractures formed in each soil layer provided channels for liquefaction dewatering,resulting in formation of macroscopic liquefaction deformation,such as liquefied lime-mud volcanoes,liquefied diapir structures,vein-like liquefied structures and liquefied curls,etc. Splendid liquefied lime-mud eruption lasted for two to three hours,which is similar to the sand volcano eruption induced by strong earthquake. However,under the same artificial seismic conditions,development of macroscopic liquefied structures in three experimental models varied in shape,depth and quantity,indicating that excess pore water pressure ratios at initial liquefaction stage and complete liquefaction varied with depth. With size increasing of calcite particle in lime-mud,liquefied depth and deformation extent increase accordingly. The simulation test verifies for the first time that strong earthquakes may cause violent liquefaction of saturated lime-mud composed of micron-size calcite particles,uncovering the puzzled issue whether seafloor lime-mud can be liquefied under strong earthquake. This study not only provides the latest simulation data for explaining the earthquake-induced liquefied deformations of saturated lime-mud and seismic sedimentary events,but also is of great significance for analysis of foundation stability in marine engineering built on the soft calcium carbonate layers in neritic environment.     

Study of shaking table test of seismic subsidence loess landslides induced by the coupling effect of earthquakes and rainfall

Light rain or moderate rain is the most common meteorological event in the rainy season in the loess area of China, so the probability of landslide hazards induced by the coupling effect of earthquakes and rainfall under the condition of light rain or moderate rain is relatively higher than that under heavy rain. To study the dynamic response characteristics and instability mechanism of loess slopes by the coupling effect of earthquakes and rainfall under the conditions of moderate rain and light rain, a low-angle slope model test of a large-scale shaking table after 10 mm of rainfall was carried out. By gradually increasing the dynamic loading, the evolution of the macroscopic deformation and the instability failure mode of the slope model are observed; the temporal and spatial trends of the amplification effect, acceleration spectrum, pore pressure and soil pressure are analyzed; and the failure mechanism of the slope is determined. The results showed that the amplification effect increased along the slope surface upward, and a strong amplification effect appeared at the front of the top of the slope. Because of the stronger dynamic stress action on the upper part of the slope, the immersed soil in the upper part of the slope experienced seismic subsidence deformation, the saturation in the seismic subsidence soil increased, and the water content temporarily increased locally. With the further increase in the loading intensity, a large number of tension cracks were generated in the seismic subsidence area, and water infiltrated down along the cracks and the wetting range expanded under dynamic action. The range of seismic subsidence and cracks further extended to the deep part of the slope. Under the reciprocating action of the subsequent ground motion, the swing amplitude of the soil mass in the seismic subsidence area, which is divided by a large number of cracks in the upper part of the slope, increased further, resulting in the further reduction in the residual strength of the seismic subsidence soil mass located at the crack tip due to the pull and shear action. Finally, under the combined action of gravity and dynamic force, the upper soil mass in the seismic subsidence area dragged the lower soil mass in the seismic subsidence area downward because the sliding force is greater than the residual strength of the soil mass, which induced a seismic subsidence-type loess landslide. Under the coupling effect of earthquakes and rainfall, the instability mode and mechanism of this landslide are significantly different from those of liquefaction-type landslides.

     

Determining the origin of soft‐sediment deformation structures: a case study from Upper Carboniferous delta deposits in south‐west Wales,UK

Determining the cause of sediment mobilization is a major problem; possible triggers include earthquakes, sediment loading and wave action. A detailed sedimentological and palaeoenvironmental analysis of soft‐sediment deformation in Upper Carboniferous deltaic deposits in SW Wales, UK, shows that two styles of deformation occur. Type A (syndepositional convolute stratification) affects most sandstone beds and was generated by rapid sedimentation. Type B (localized sand‐in‐sand pseudonodules) incorporates beds that already contained Type A deformation, and developed when the substrate was liquefied by disturbance due to movement on a near‐surface gravity slide. Neither type of deformation was triggered by seismic events.     

轴力作用下液化土中端承桩的水平振动特性

将地震液化场地土层分为非液化表层土、中部的液化土层和底部的基层,基于饱和多孔介质理论和Novak薄层法,研究轴向压力作用下液化黏弹性土层中端承桩的水平动力特性。利用Helmholtz分解和变量分离法,得到液化土层对桩水平振动的阻抗。利用矩阵传递法,在频率域得到轴力作用下液化土层中端承桩简谐振动的解析解和桩头复刚度的表达式,并进行参数研究,分析轴力、桩-土模量比、桩长径比、液-固耦合系数等对桩头动力刚度和阻尼的影响。结果表明,在轴力作用下,不同长径比、桩-土模量比、液-固耦合系数时的动力刚度绝对值均比无轴力作用时减小,但随频率的变化趋势相同;轴力对桩水平振动的动力阻抗影响显著,随着轴力的增加,桩的水平振动动力刚度因子的绝对值减小,若轴力继续增大,其绝对值趋近于0,桩发生失稳破坏;桩长径比和桩土模量比对桩的水平振动动力阻抗有显著的影响,而液-固耦合系数的影响较小。     

微倾斜场地中地铁地下结构周围地基液化与变形特性振动台模型试验研究

强震中场地砂土液化产生的土层侧移对地面建筑结构和地下生命线工程造成了严重的破坏。可以预见,微倾斜液化场地的土层侧移也将对地铁地下结构的地震安全造成严重的威胁。鉴于此,开展了微倾斜（倾角为6o）可液化场地中两层三跨地铁地下车站结构与区间隧道连接部位地震反应的大型振动台模型试验研究。结果表明：微倾斜可液化场地中地铁车站结构两侧地基出现了明显的非对称液化分布特征,坡体下方水平土层比上方水平土层更易液化;因坡体内土体液化沿坡向下流滑引起了下方水平土层发生了明显的地面抬升,总体上坡体段内的地面侧移量最大,下方水平土层地面侧移量次之,坡体上方水平土层地面侧移量最小。同时,在试验过程中也发现,隧道和车站结构之间发生了明显的差异上浮,可能会造成连接部位附近结构的应力集中或加重该部位的地震破坏。     

苏北盆地高邮凹陷戴南组震积岩特征及其石油地质意义

高邮凹陷戴南组沉积时期,在中-新生代构造运动的影响下,古地震活动强烈,岩层中发育了一系列与古地震有关的准同生变形构造。诸如液化岩脉、液化角砾岩、枕状层、环状层理及震褶层等液化挤压变形构造;又可见阶梯状断层、震裂岩、自碎角砾岩、内碎屑副角砾等脆性变形构造。另外,在岩芯观察的过程中发现,这些地震作用变形构造层往往与地震作用下的重力流沉积岩层相伴生,在此基础上通过总结区域联井和单井序列,刻画出了研究区震积岩的垂向序列。结合研究区油气成藏特征,从油气运移、储集和封堵三个方面对地震事件的油气地质意义进行了讨论,在此基础之上,提供了一个新的勘探方向,打破以往"环凹找油"的勘探思路,认为深凹带砂体也具有重要的勘探价值。