Discrete element modeling of static liquefaction of shield tunnel face in saturated sand

The failure mechanism and limit support pressure are essential factors for the stability assessment of tunnel face during shield tunneling. Previous studies were mostly based on the assumption that soil fails when the plastic limit is reached. The static liquefaction of saturated sand at pre-failure state often largely reduces shear strength; a stable tunnel face could be destabilized and even collapses abruptly. This paper aims to explore the static liquefaction failure of tunnel face in saturated sands by 3D discrete element simulation. The particle shape is considered by using clump composed of two identical spheres, and the micro-material parameters are calibrated by fitting against triaxial tests. The tunnel face is initially supported by the earth pressure at rest, and then the pressure is reduced by applying a displacement on the free surface until soil body completely collapses. The second-order work during the unloading process is calculated, and a negative value denotes the initiation of static liquefaction. The distribution of negative second-order work point shows the initiation and propagation of static liquefaction. The ultimate failure mode is shown to be composed of a sliding wedge and an overlying chimney, and the chimney is confined to a local area for soil arching. The experienced stress path of tunnel face verifies the existence of static liquefaction instability, indicating the stability analysis should consider static liquefaction rather than simply using the conventional plastic limit state analysis.

     

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

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

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

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

煤造石油与煤成油理论关系研究进展综述

综述了煤造石油的研究进展。煤加氢液化研究结果表明挥发分、Ｈ与Ｃ比值、煤级、催化剂是影响煤炭液化的重要因素;对煤造石油与石油玩机成因说、煤造石油与煤一理论之间的联系进行了评述,现代含煤岩系煤成油研究得出的结论与煤造石油实践得出的结果近似一致,一些大油田分布带与大构造活动区带近于一致,由此使人们得出了类似于Ｆ－Ｔ法合成的石油无机成因理论。     

Characteristics of the unloading process of rocks under high initial stress

The unloading process of rocks under high initial stress is complex, and verifying the mechanism of the unloading process in the field or in a laboratory is not straightforward. In this study, the unloading process of rocks under high initial stress was characterised by a mathematical physics model, which was then implemented in the finite element program LS-DYNA for analysis. In particular, the implicit and explicit methods were performed in sequence in the finite element simulation of rocks with initial stresses. In the numerical simulation, the characteristics of the dynamic unloading process of rocks were investigated for various peak initial stresses, initial stress release paths and initial stress release rates (ISRRs). The numerical results indicated that the rock failure could be induced by the release of the initial stress; furthermore, there is a relationship between the magnitude of the unloading failure and the peak initial stresses, the initial stress release paths and the ISRRs. When the initial stresses were at the same level, the equivalent initial stress release rate (EISRR) was introduced to quantitatively describe the characteristics of the unloading process. Using the numerical results, the unloading failure process was characterised, and a method for the static stress initialisation-dynamic unloading of rock was developed.     

蒙古南贝尔凹陷白垩纪地震作用的证据

在区域构造背景研究和岩心观察的基础上,在蒙古南贝尔凹陷白垩纪地层中识别出震积岩。震积岩的主要标志是发育各种类型的软沉积变形构造,包括微同沉积断裂、卷曲变形构造、液化砂岩脉、震积枕状及伴生构造等,这是首次在蒙古南贝尔凹陷发现与地震作用有关的沉积特征。在塔21 9井中发现完整的震积岩垂向序列,序列自下而上分为下伏未震层、砂岩脉（砂火山）、震动卷曲变形层、内碎屑角砾岩层、阶梯状微断层和震裂缝层、震浊积岩层和上覆未震层。如此完整地震序列在同一口井岩芯中同时出现实属罕见,为进一步深入研究古地震作用提供了良好的物质组合。     

液化型路堤边坡动力数值模拟分析

液化型路堤边坡动力稳定性问题涉及岩土工程与工程地震两个学科领域,是边坡工程与砂土液化的交叉课题。采用天然地震记录为输入条件,应用Finn本构关系模型,运用有限差分法,对填土+砂土+卵砾土地层组合的路堤边坡进行了全时程动力分析,探讨了地震作用下路堤边坡的液化初步规律和稳定性。数值模拟结果表明:地震作用引起了路基饱和砂土有效应力急剧减小,并导致路基砂土液化,引起路堤变形破坏。孔隙水压力的积累与消散不仅与地震记录序列存在对应关系,也与砂土所处的位置和深度有密切关系。地表变形破坏主要表现为路堤顶面发生震陷和拉裂破坏,坡底面产生挤压隆起变形。地面以下的变形破坏主要包括土体剪切破坏和深部砂土液化引起的侧向流动破坏。     

液化场地桩基桥梁震害响应大型振动台模型试验研究

采用大型振动台进行液化场地桩基桥梁震害响应模型试验,很好再现了自然地震触发场地液化及结构破坏的各种宏观现象。0.15gEl Centro波输入下,上部砂层局部液化,桩-柱墩加速度主要表现为低频反应,桩动应变幅值自下而上很快增大、到达地表则大幅度减小。0.5gEl Centro波输入下,整个砂层全部液化,桩被折断且加速度也主要表现为低频反应,桩动应变幅值自下而上很快增大、到达地表则大幅度减小。砂层液化与否,对桩-柱墩动力反应影响很大。     

高地应力区砂岩在卸荷条件下的变形参数劣化试验研究

以锦屏某高边坡砂岩所赋存的高地应力环境为基础,开展峰前卸围压力学特性试验,研究了应力-应变全过程曲线、变形特征及参数劣化效应。试验结果表明,相对于加载试验,岩样卸荷破坏表现出明显的脆性特征,相同初始围压下卸荷破坏所需偏应力较加载试验少;卸荷过程横向应变和体积应变急剧增大,卸荷方向表现出明显的扩容特征;当初始围压小于某一值（30 MPa）,岩样破坏所需卸荷量随初始围压增大而增大,当初始围压达到一定值（40 MPa时）,脆性特征更加明显,很少的卸荷量即可引起岩样破坏。以表征岩石卸荷程度的卸荷量H为参变量,分别对高、低围压下卸荷过程中变形参数（变形模量和泊松比）随卸荷量变化关系进行拟合,得到高、低地应力区卸荷过程中变形劣化参数的关系式。试验结果和计算方法对高应力区边坡开挖的稳定性评价具有一定指导意义。     

地下岩体洞室群地震响应的加、卸载响应比分析

将加、卸载响应比理论引入地下岩体洞室群地震响应分析中。以输入的地震加速度作为加、卸载,以洞室围岩的加速度作为响应,合理地确定了加、卸载响应区段。通过有限差分程序FLAC3D建立了白鹤滩水电站13号机组剖面数值分析模型,选用汶川地震波进行动力时程计算,探讨地下岩体洞室群的地震动力稳定性,研究结果表明：在100年超越概率2%的地震作用下,白鹤滩水电站地下洞室群围岩响应比峰值范围在4.05～11.52之间;结合应力及位移分析,主厂房拱顶及层间错动带C4在尾调室上游边墙出露部位的围岩均进入了非线性变形状态,但响应比在整个时程中并没有趋于无穷大,围岩没有发生失稳破坏;错动带两侧的围岩产生了一定的相对位移,会对尾调室的边墙稳定性产生不利影响。该研究方法可应用于一般地下岩体结构的地震稳定性分析中。     

复合桩基重复加卸载过程中侧向约束桩变位规律试验研究

复合桩基上堆载可能引起侧向约束桩侧移,导致侧向约束桩和复合桩基失稳破坏。采用室内模型试验研究复合桩基重复加、卸载过程中侧向约束桩水平变位规律。结果表明：（1）复合桩基上加载时,侧向约束桩的桩身侧移沿深度先增大、后减小、存在峰值,峰值随加载增大而增大,出现在距离地面0.4倍地面以下桩长处,峰值位移与桩顶（或地面）处位移的比例系数和发生侧移的桩身长度随荷载增大而增大;（2）复合桩基上重复加、卸载时,侧向约束桩的桩顶侧移随桩间距增大而增大,弹性变形则变小。中桩侧移比边桩大。桩顶侧移增长率随重复加载次数增加而减小。若荷载超过首次加载的最大荷载,侧向复合地基上加载-桩水平位移曲线将回到首次加载曲线的延长线,具有记忆效应。第 次卸载曲线与第i次卸载曲线线型相似且相互平行。每次卸载的初、中期,荷载的减小不影响桩顶侧移,只有卸载到最后1～2级时,侧移才开始减小,直到永久塑性变形。侧向约束桩的顶部水平位移回弹曲线特征与岩石或土在垂直加、卸载作用下的回弹曲线相似。     

Liquefaction analysis and soil improvement in Beydag dam

Beydag dam is under construction on Kucukmenderes River for irrigation purposes. Due to the scarcity of core material and liquefaction of alluvium at the dam site, the original design was changed to Roller Compacted Concete (RCC) from rockfill dam with claycore. Although the new design was safer, it nearly doubled the cost of the dam, so the owner, State Hydraulic Works of Turkey, (DSI) set out to find more economical but equally safe alternative. Since jet-grouting is a cheap ground improvement tool in Turkey, such an alternative was developed for the ground improvement against liquefaction together with concrete face rockfill dam sitting on top of improved ground. This paper presents a detailed discussion of how the new alternative was developed and evaluated: it discusses the determination of jet grouting pattern, the placement of jet grouted blocks, and the assesment of liquefaction. On one hand the soil cement strength of jetgrout columns, internal friction angle of alluvium and rockfill were important in determining the dimensions of the blocks, on the other hand the location of the blocks were highly affected by the areas where liquefaction occurred. One of the most important parameter that has a considerable influence in delineating the boundary betweeen liquefaction and non-liquefaction was the value of stress reduction coefficient (r _d), being primarily sensitive to the weight of overburden, which is calculated by the height from the face of dam to the depth where the calculation was made. This approach is justified by two-dimensional ground response analysis. Most importantly, this paper shows that there exists an alternative solution for building dams on the liquefaction prone sites without removing alluvium by using a well known jet grouting technique for improving ground at only selected places.     

鄂尔多斯盆地三叠系延长组震积岩特征研究

震积岩是具有古地震记录的岩层。根据大量岩心观察,鄂尔多斯盆地三叠系延长组广泛发育震积岩沉积标志,如液化砂岩脉、阶梯状小断层、液化卷曲变形、砂球枕构造、泥岩撕裂屑、滑揉构造等。详细分析了各种震积岩特征及其分布,初步建立了鄂尔多斯盆地延长组震积岩垂向沉积序列,包括:A段为下伏未震层;B段为微同沉积断裂层;C段为振动液化卷曲变形层;D段为液化砂岩脉和砂岩墙段;E段为砂球枕及碎块层;F段为液化均一层和G段上覆未震层。震积岩的首次发现为研究盆地周边构造演化强度和期次提供了佐证,同时震积岩也是一种潜在的有利储层。     

表层液化地基中结构与桩惯性相互作用简化分析方法

传统桩基础结构抗震设计中一般认为忽略结构与桩基础相互作用影响使得结构所受惯性作用计算结果偏于保守,但对几次强震中地震动反应谱和桩基础结构震害特征的总结与分析表明,液化作用和结构与桩基础相互作用可能导致桩基础结构惯性作用增强。通过可液化地基中桩基础结构的动力离心模型试验手段,分析了地基液化对结构与桩惯性相互作用的影响,发展了适用于表层液化地基中的结构与桩惯性相互作用的简化计算方法,并利用该方法对文中及Wilson等的离心模型试验进行分析,得到的结构加速度与试验结果吻合较好,从而验证了方法的有效性。     

水平液化场地土表位移简化理论解答

水平液化场地地表往返位移的求解是基础和地下工程抗震设计的迫切需求。利用双层模型模拟实际水平场地,提出了可考虑液化层存在下的土表位移简化计算方法并给出频域理论解答,同时采用逐循环累计方法给出水平场地任意荷载下土表位移时域解答。采用适于水平场地孔压增量模型,并逐循环修正土层模量,模拟液化引起的土层的非线性过程。振动台试验结果和提出棋型的计算结果吻合,验证了所提方法能够反映土体液化对土表位移影响的基本过程。提出的方法和解答物理意义明确,可代表土体液化对土表位移影响的基本形态,并可用于频域及时域的无量纲分析。     

山西云岗石窟侏罗系地震液化砂岩柱的发现及其大地构造意义

本文主要报道在华北克拉通北缘云岗石窟景区内侏罗系地层中发现的20余个由地震液化形成的砂岩柱。地震导致的砂土液化是非常普遍的自然现象,地震波的震动使埋藏在地下未固结的饱和砂质沉积物迅速变为流体并在巨大压力下喷涌至地表,形成砂火山或泥火山。2008年汶川8. 0级大地震和2012年新西兰6. 2级地震都产生了严重的液化现象,在对地表造成严重破坏的同时,分别在地表形成了一系列的溢出丘和形态完好的砂火山。古地震之后,液化砂质沉积物往往在上涌的通道内固结成岩,形成穿层的岩柱或岩管。近年来,美国科罗拉多大峡谷相继发现很多地震液化形成的中生代巨型碎屑岩柱。2018年,本文作者在云岗石窟景区内的侏罗系云岗组中发现了20多条直径15~20cm左右、高可达2m以上穿层“侵入”的砂岩柱,有的在纵向上呈串珠状排列,有的呈不规则的树枝状自下向上伸展,膨大狭缩、分支复合现象非常明显,砂岩柱内部结构均一,不见任何层理。砂岩柱顶端的泥质围岩石中有明显向上牵引、拖拽的痕迹。这些砂岩柱没有任何硅化木的特征,也不是普通成岩作用形成的结核,与暴雨、泥石流、重力滑塌等作用均无任何关系,是古地震液化作用在地层内保存的遗迹,故称为液化砂岩柱。地震液化砂岩柱发育的层位位于李振宏等人(2014)在宁武—静乐盆地侏罗系云岗组顶部的凝灰质泥晶碳酸盐层位(160 Ma)之下约40m处,初步判定形成这些液化砂岩柱的古地震发生于160. 75 Ma~160 Ma之间的晚侏罗世。古地震的形成过程与燕山运动引起的地壳运动有直接的关联,是华北克拉通内与燕山运动伴生的古地震在地层中留下的最直接证据。     

Paleoseismic event recorded in the Upper Cretaceous Nenjiang Formation in southeastern area of the Songliao Basin (NE China)

The Yaojiachezhan (YJCZ) section and NGN1, Zk3389 and SK1-S wells, which are located in the southeastern part of the Songliao Basin, NE China, were used to analyse paleoseismic events recorded in the Upper Cretaceous first member of the Nenjiang Formation. Two representative seismic sequences are recognised in these lacustrine deposits. The seismic sequence observed in the YJCZ section, from the base to top, includes step micro-faulting, seismic fissures, liquefaction sandstone veins, brecelated structures, liquefied fold-deformed laminations and sand pillow structures. Whereas in NGN1 core, from the base to top, includes micro-faulting, seismic fissures, sand pillow structures, liquefied fold-deformed laminations and liquefaction sandstone veins. These two seismic sequences together reflect the processes resulting from a strong seismic event from initiation to culmination, then decline and recession. The occurrence of brecciated structures reflects the earthquake magnitude in the YJCZ section area probably reached M7.5, while liquefied fold-deformed laminations in the NGN1 core area reflects an earthquake magnitude of M7.0. Empirical relations between earthquake magnitude and the maximum distance of liquefaction deformation structure sites from the epicentres show that the longest distance between the earthquake source area for the YJCZ section and NGN1 well are <100?km, whereas the distance for the Zk3389 and SK1-S wells is >100?km. The trigger source area for this paleoseismic event is most likely the Yilan-Yitong Fault, meaning that the Yilan-Yitong Fault was active during the deposition of the first member of the Nenjiang Formation.     

近远场地震动作用下地铁车站结构地基液化效应的振动台试验

开展了近场和远场地震动作用下3跨3层地铁车站结构地基液化效应的振动台模型试验,测试了地铁车站结构的加速度、应变、水平位移反应和地基土孔隙水压力、加速度、震陷及其作用于模型结构侧墙的动土压力反应。分析和总结了地铁车站结构地基液化效应特征,结果表明：模型结构对其周围地基土孔隙水压力场的分布有明显影响,结构两侧和底部地基土中的孔压峰值小于相同深度离结构较远地基土中的孔压峰值;地基土中孔压的消散速度自下而上呈逐渐减慢的趋势;地震动作用过程中,模型结构产生向上的相对运动,强地震动作用时模型结构上浮现象明显;模型结构侧墙受到的动土压力随深度增大而减小,输入地震动特性对动土压力的大小有显著影响。     

软岩三轴加-卸载试验的破坏特征及抗压强度取值方法研究

软岩三轴加、卸载试验中的抗压强度参数取值是一个非常关键的问题,目前还没有比较确定的方法,通常根据经验或工程要求确定。基于此,结合一组砂质泥岩的三轴加、卸载试验进行了详细研究,试验结果表明:(1)当围压大于10 MPa时,不论是三轴加载试验还是卸载试验,砂质泥岩均表现出明显的应变硬化特性;(2)砂质泥岩在三轴压缩和三轴卸荷两种应力路径的力学响应差别明显,在卸荷条件下通常只存在一条完整的剪切破坏面,而压缩条件下岩样没有明显的控制破裂面,并且剪胀效应明显;(3)将半对数法引入到软岩三轴加、卸载试验数据分析中,可以比较方便地确定软岩三轴加、卸载试验的抗压强度;(4)基于三轴加、卸载试验过程中的应力、应变增量分析,将软岩变形破坏过程分为3个阶段:平滑阶段、规则跳跃阶段和急剧跳跃阶段,各阶段划分的意义明确,可以作为抗压强度确定时的重要参考,因此,在软岩三轴试验数据分析时,可以综合采用这两种方法。研究提出的抗压强度确定方法可操作性强,在其他软岩三轴试验数据分析中值得参考借鉴。     

Liquefaction potential of coastal slopes induced by solitary waves

Tsunami runup and drawdown can cause liquefaction failure of coastal fine sand slopes due to the generation of high excess pore pressure and the reduction of the effective over burden pressure during the drawdown. The region immediately seaward of the initial shoreline is the most susceptible to tsunami-induced liquefaction failure because the water level drops significantly below the still water level during the set down phase of the drawdown. The objective of this work is to develop and validate a numerical model to assess the potential for tsunami-induced liquefaction failure of coastal sandy slopes. The transient pressure distribution acting on the slope due to wave runup and drawdown is computed by solving for the hybrid Boussinesq—nonlinear shallow water equations using a finite volume method. The subsurface pore water pressure and deformation fields are solved simultaneously using a finite element method. Two different soil constitutive models have been examined: a linear elastic model and a non-associative Mohr–Coulomb model. The numerical methods are validated by comparing the results with analytical models, and with experimental measurements from a large-scale laboratory study of breaking solitary waves over a planar fine sand beach. Good comparisons were observed from both the analytical and experimental validation studies. Numerical case studies are shown for a full-scale simulation of a 10-m solitary wave over a 1:15 and 1:5 sloped fine sand beach. The results show that the soil near the bed surface, particularly along the seepage face, is at risk to liquefaction failure. The depth of the seepage face increases and the width of the seepage face decreases with increasing bed slope. The rate of bed surface loading and unloading due to wave runup and drawdown, respectively, also increases with increasing bed slope. Consequently, the case with the steeper slope is more susceptible to liquefaction failure due to the higher hydraulic gradient. The analysis also suggests that the results are strongly influenced by the soil permeability and relative compressibility between the pore fluid and solid skeleton, and that a coupled solid/fluid formulation is needed for the soil solver. Finally, the results show the drawdown pore pressure response is strongly influenced by nonlinear material behavior for the full-scale simulation.