饱和砂土液化后大变形机理的离散元细观分析

采用颗粒流软件模拟了饱和砂土在不排水条件下的循环剪切试验,研究了不同因素对液化的影响,并进一步分析了饱和砂土液化后宏观变形的基本规律。在此基础上,从孔隙分布角度解释了砂土液化后的大变形的细观物理机制。通过自编程序对颗粒排列和孔隙分布的演化过程进行定量描述,给出孔隙率标准差作为液化后体积收缩势的度量,并研究了孔隙率标准差与液化后大变形的关系。结果表明,试样初始状态并不会影响液化后的极限状态,仅对初始液化所需循环次数有影响。数值模拟可再现室内试验中的宏观现象,证实了室内试验中饱和砂土液化后的有限剪切大变形是客观真实的材料响应。土体体积收缩势的累积所导致的孔隙均匀化以及土颗粒间自由空隙增大正是饱和砂土液化后循环剪应变逐渐增大的细观机理。孔隙率标准差作为孔隙均匀化的量化指标,与循环剪应变各周次幅值有良好相关性。     

饱和砂土液化后大变形机理的离散元细观分析

采用颗粒流软件模拟了饱和砂土在不排水条件下的循环剪切试验,研究了不同因素对液化的影响,并进一步分析了饱和砂土液化后宏观变形的基本规律。在此基础上,从孔隙分布角度解释了砂土液化后的大变形的细观物理机制。通过自编程序对颗粒排列和孔隙分布的演化过程进行定量描述,给出孔隙率标准差作为液化后体积收缩势的度量,并研究了孔隙率标准差与液化后大变形的关系。结果表明,试样初始状态并不会影响液化后的极限状态,仅对初始液化所需循环次数有影响。数值模拟可再现室内试验中的宏观现象,证实了室内试验中饱和砂土液化后的有限剪切大变形是客观真实的材料响应。土体体积收缩势的累积所导致的孔隙均匀化以及土颗粒间自由空隙增大正是饱和砂土液化后循环剪应变逐渐增大的细观机理。孔隙率标准差作为孔隙均匀化的量化指标,与循环剪应变各周次幅值有良好相关性。     

饱和砂土液化特性的可视化试验研究

针对目前散粒体液化细观机制研究的局限性,对C.K.C型动态三轴仪进行改进,研制了可用于研究散粒体细观组构特征的动三轴可视化试验系统。利用该系统对循环荷载作用下标准砂的颗粒运动规律及细观组构特征进行初步探讨,从颗粒的运动规律、定向性、配位数和孔隙率等细观角度分析了循环荷载作用下饱和砂土发生液化的细观机制,认为液化的发生是土体在循环荷载作用下颗粒不断运动、重新排列的结果,指出颗粒长轴方向、配位数和孔隙率等细观组构参量是反映循环荷载作用下饱和砂土液化宏观响应的重要参量     

循环荷载下砂土液化特性颗粒流数值模拟

利用PFC2D常体积循环双轴试验条件,对砂土在不排水循环荷载作用下的液化特性进行了颗粒流数值模拟,数值模拟按等应力幅加荷方式进行。颗粒流数值模拟的优点在于得到试样液化宏观力学表现的同时,通过不同循环加荷时刻试样内细观组构参量（包括配位数、接触法向分布、粒间法向接触力、粒间切向接触力）的演化规律,分析砂土液化过程中细观组构变化与宏观力学响应之间的内在联系,从而可进一步探讨砂土液化的细观力学机制。数值模拟研究结果表明,砂土液化现象在宏观力学表现上反映为超静孔隙水压力的累积上升和平均有效主应力的不断减小,在细观组构上对应于配位数的累积损失和粒间接触力的不断减小。砂土液化细观机制分析表明,试样配位数的减少与循环加荷过程中组构各向异性滞后于应力各向异性有关。     

饱和层状砂土液化特性的动三轴试验研究

利用GDS动三轴试验系统采用等幅循环应变加载方式对含有不同厚度粉土的饱和层状砂土进行了液化强度试验。分析了均匀砂和含有不同粉粒层厚度的层状砂土在循环荷载作用下的变形和力学特性。试验分析表明：由于含粉粒夹层的层状土特殊的土体结构,其孔隙水压力发展规律与一般的无黏性砂土不同;饱和层状砂土的抗液化强度并不是随着粉粒层厚度的增加而单调增加的,而是存在一个临界点;液化临界剪应变的大小与液化判别标准和循环次数有很大关系。试验结果表明,粉粒夹层对层状砂土的液化特性有很大的影响,且更能模拟自然环境条件下的层状砂土地基液化特性。     

动扭剪试验中砂土液化后流动特性分析

土体在初始液化后仍然可能承受动荷载作用发生大变形。引入流体力学中的剪应变率和表观黏度的概念,对振动扭剪试验中饱和砂土液化后的流动特性进行了分析。分析中将砂土的状态分为0有效应力状态和非0有效应力状态。结果表明,砂土在0有效应力状态下表现出与静扭剪试验类似的“剪切稀化”非牛顿流体的特征,表观黏度随着剪应变率的增大而减小。加载周数对“剪切稀化”状态下的剪应变率幅值有影响,随着加载周数的增加,剪应变率幅值逐渐增大,而对流动曲线的形状没有影响。在非0有效应力状态下,砂土的表观黏度随应变的增大而增大,随孔压比的减小而增大,且所有的试验得到的表观黏度与孔压比具有一致关系。     

纤维加筋砂土抗液化试验与数值模拟

砂土液化是导致重大地震灾害的主要原因之一。本研究探讨了天然纤维加筋砂土在循环荷载作用下的抗液化性能。在不排水条件下,对具有不同纤维含量的加筋砂土试样进行了一系列循环三轴试验,研究了饱和砂土的液化特性以及循环剪应变幅值、纤维含量对饱和砂土抗液化性能的影响。此外,通过模拟已完成的循环三轴试验,建立了二维有限元数值模型,并对具有不同纤维含量的加筋砂土进行了参数标定。研究结果表明:(1)增加循环剪应变幅值将促进超孔压累积,使得滞回曲线斜率和平均有效应力减小速度加快;(2)纤维的存在能够减缓超孔压的累积,随着纤维含量增加,加筋砂土抗液化能力得到明显提高;(3)标定后的本构模型参数能可靠地用于模拟纤维加筋砂土的液化响应。研究结果为饱和砂土抗液化问题与纤维加筋砂土的数值模拟提供了有价值的参考。     

地震历史对砂土抗液化性能影响的试验研究

为探讨地震历史对饱和砂土抗液化性能的影响,通过离心机动力模型试验,观测了遭受过不同强度地震的饱和砂土水平场地在再次地震荷载作用下的响应,包括超静孔隙水压力的发展和土体的变形。研究表明,小地震有利于增强砂土颗粒之间的咬合及结构的稳定性,极大地提高了其抗液化能力,而遭受强地震并且发生液化的砂土,由于在沉积过程中形成的不稳定结构,则可能在将来强度相对较小的地震中再次液化。试验结果说明了地震历史对砂土液化性能有重大影响及进一步开展砂土细观结构研究的重要性。     

砂土力学性质的细观模拟

基于颗粒流理论,对砂土的室内双轴试验和砂土剪切带形成与发展进行了数值模拟,基本再现了砂土试样应力-应变关系。主要研究了颗粒粒径、颗粒摩擦系数等细观参数变化时试样宏观性质的变化情况。对比分析了室内试验和颗粒流数值模拟试验的砂土剪切带特征。同时对比研究了模型试样颗粒粒径、颗粒刚度和摩擦系数等细观参数的变化,对剪切带的形成与发展过程的影响。结果表明：通过颗粒流数值模型试验可以有效模拟砂土剪切带的形成与发展机理。     

平均粒径对砂土剪切特性的影响及细观机理

针对平均粒径对砂土剪切特性的影响作用,结合室内试验和离散元模拟方法对不同平均粒径砂土进行了细观研究。基于3种不同平均粒径砂土的直剪试验结果,通过建立反映砂土剪切试验特征的PFC（particle flow code）颗粒流模型,详细研究了不同粒径砂土在剪切过程中土样体积变化、力链网络、孔隙率和配位数等细观结构参数的变化特征和规律,并从细观角度分析了颗粒粒径对土样宏观剪切特性的影响机制。结果表明：具有不同平均粒径砂土的细观结构参数在剪切过程中存在显著差异,并且其细观参数差异主要集中体现在剪切带处;剪切力学特性的影响主要体现在抗剪强度和剪胀效应方面,砂土平均粒径越大,抗剪强度越高,剪胀效应越明显;具有不同平均粒径的砂土在剪切过程中土颗粒运动规律及剪切带形态变化特征存在一定的差异,平均粒径越大,剪切带内上跨式颗粒占比越大,剪切带厚度越大。     

Large post-liquefaction deformation of sand, part I: physical mechanism, constitutive description and numerical algorithm

This paper presents a theoretical framework for predicting the post-liquefaction deformation of saturated sand under undrained cyclic loading with emphasis on the mechanical laws, physical mechanism, constitutive model and numerical algorithm as well as practical applicability. The revealing mechanism behind the complex behavior in the post-liquefaction regime can be appreciated by decomposing the volumetric strain into three components with distinctive physical background. The interplay among these three components governs the post-liquefaction shear deformation and characterizes three physical states alternating in the liquefaction process. This assumption sheds some light on the intricate transition from small pre-liquefaction deformation to large post-liquefaction deformation and provides a rational explanation to the triggering of unstable flow slide and the post-liquefaction reconsolidation. Based on this assumption, a constitutive model is developed within the framework of bounding surface plasticity. This model is capable of reproducing small to large deformation in the pre- to post-liquefaction regime. The model performance is confirmed by simulating laboratory tests. The constitutive model is implemented in a finite element code together with a robust numerical algorithm to circumvent numerical instability in the vicinity of vanishing effective stress. This numerical model is validated by fully coupled numerical analyses of two well-instrumented dynamic centrifuge model tests. Finally, numerical simulation of liquefaction-related site response is performed for the Daikai subway station damaged during the 1995 Hyogoken-Nambu earthquake in Japan.     

Effects of initial static shear on cyclic resistance and pore pressure generation of saturated sand

In practical engineering, cyclic shear stresses induced by earthquakes, traffic, and waves are superimposed on the initial static shear stress in sand fills or deposits, leading to complex responses of soils such as their deformation characteristics, pore pressure generation, and susceptibility (or cyclic resistance) to liquefaction. To experimentally investigate the undrained cyclic response of saturated sand, a series of triaxial tests were performed, covering a broad range of initial static and cyclic deviatoric stress levels. The results indicate that different stress conditions lead to two types of cyclic behavior: cyclic mobility and residual deformation accumulation. The compressional static stress is beneficial to the cyclic resistance of the dense sand, whereas the extensional static stress is regarded as detrimental as it tended to reduce the cyclic strength. Moreover, by comparing the available test data obtained from the same sand with varying initial densities and confining pressures, the static shear effect on cyclic resistance was found to be dependent on the state of the sand. Compared to the interpretation made using the limiting pore pressure-based criterion, the conventional failure criterion using a cyclic axial strain of 5% may lead to a substantial overestimation of the cyclic resistance, thus resulting in unsafe assessment and design. Hence, by employing the pore pressure criterion, the pore pressure generated in the cyclic tests was investigated and was found to be significantly influenced by the static shear stress. A pore pressure generation model is proposed to obtain the pore pressure characteristics of sand under various static shear stress conditions.     

DEM study of fabric features governing undrained post-liquefaction shear deformation of sand

In an effort to study undrained post-liquefaction shear deformation of sand, the discrete element method (DEM) is adopted to conduct undrained cyclic biaxial compression simulations on granular assemblies consisting of 2D circular particles. The simulations are able to successfully reproduce the generation and eventual saturation of shear strain through the series of liquefaction states that the material experiences during cyclic loading after the initial liquefaction. DEM simulations with different deviatoric stress amplitudes and initial mean effective stresses on samples with different void ratios and loading histories are carried out to investigate the relationship between various mechanics- or fabric-related variables and post-liquefaction shear strain development. It is found that well-known metrics such as deviatoric stress amplitude, initial mean effective stress, void ratio, contact normal fabric anisotropy intensity, and coordination number, are not adequately correlated to the observed shear strain development and, therefore, could not possibly be used for its prediction. A new fabric entity, namely the Mean Neighboring Particle Distance (MNPD), is introduced to reflect the space arrangement of particles. It is found that the MNPD has an extremely strong and definitive relationship with the post-liquefaction shear strain development, showing MNPD’s potential role as a parameter governing post-liquefaction behavior of sand.     

冻融循环作用下砂岩孔隙体积变形模型的建立与分析

为研究冻融作用下砂岩内部孔隙的扩展机制,基于弹性力学和热力学理论,建立了饱水孔隙内水结冰的体积变形计算模型。结合基于核磁共振技术的冻融循环作用下饱和砂岩的孔隙度变化和T₂谱分布的试验结果,通过建立的孔隙体积变形计算模型,分析了冻融作用下砂岩的孔隙体积变形规律。结果表明：冻融作用下饱和砂岩的体积膨胀系数大于8.69%;单次冻融作用对砂岩试样的损伤较小,孔隙体积的变化范围在0.80%~1.03%之间;砂岩的孔隙收缩系数和孔隙弹性变形体积占比的变化规律一致,均为先增大后减小;砂岩试样的孔隙收缩系数在7.92%~8.04%之间,其孔隙弹性变形体积占比平均为97.79%;在冻融过程中大孔隙的弹性变形比小孔隙的弹性变形更小。     

A unified plasticity model for large post-liquefaction shear deformation of sand

Based on previous experimental findings and theoretical developments, this paper presents the formulation and numerical algorithms of a novel constitutive model for sand with special considerations for cyclic behaviour and accumulation of large post-liquefaction shear deformation. Appropriate formulation for three volumetric strain components enables the model to accurately predict loading and load reversal behaviour of sand, fully capturing the features of cyclic mobility. Compliance with the volumetric compatibility condition, along with reversible and irreversible dilatancy, allows for physically based simulation of the generation and accumulation of shear strain at zero effective stress after initial liquefaction. A state parameter was incorporated for compatibility with critical state soil mechanics, enabling the unified simulation of sand at various densities and confining pressures with a same set of parameters. The determination methods for the 14 model parameters are outlined in the paper. The model was implemented into the open source finite-element framework OpenSees using a cutting-plane stress integration scheme with substepping. The potentials of the model and its numerical implementation were explored via simulations of classical drained and undrained triaxial experiments, undrained cyclic torsional experiments, and a dynamic centrifuge experiment on a single pile in liquefiable soil. The results showed the model’s great capabilities in simulating small to large deformation in the pre- to post-liquefaction regime of sand.     

饱和剪胀性砂土液化后流滑现象与机理

基于完全不排水试验的已有研究结果表明,只有在饱和松砂中才可能发生流滑现象,较密实的饱和剪胀性砂土在初始液化后的不排水循环或单调剪切作用条件下,则会呈现出所谓往返活动性或抗剪强度恢复现象。然而,实际土工抗震边值问题中的饱和砂土在许多情况下并不能够被视为处于完全不排水条件,可能由于土性和变形等因素的不同而处于剪切排水甚至剪切吸水条件。通过对饱和砂土液化后常剪切吸水率试验结果的分析指出,剪切吸水条件可能导致剪胀性砂土的抗剪能力显著降低、甚至发生流滑现象。基于可逆性与不可逆性剪切体应变的概念以及3种体应变分量变化的特性规律,推导出剪切吸水条件下的体变约束条件。不同排水条件下的体变约束条件决定了剪胀性砂土应力-应变响应的不同变化,当满足自由吸水剪切条件时饱和砂土会出现具有流动特征的大变形;当满足强制吸水剪切条件时会发生具有非稳定特征的流滑现象。     

Modeling cyclic shearing of sands in the semifluidized state

Liquefaction is associated with the loss of mean effective stress and increase of the pore water pressure in saturated granular materials due to their contractive tendency under cyclic shear loading. The loss of mean effective stress is linked to loss of grain contacts, bringing the granular material to a “semifluidized state” and leading to development and accumulation of large cyclic shear strains. Constitutive modeling of the cyclic stress-strain response in earthquake-induced liquefaction and post-liquefaction is complex and yet very important for stress-deformation and performance-based analysis of sand deposits. A new state internal variable named strain liquefaction factor is introduced that evolves at low mean effective stresses, and its constitutive role is to reduce the plastic shear stiffness and dilatancy while maintaining the same plastic volumetric strain rate in the semifluidized state. This new constitutive ingredient is added to an existing critical state compatible, bounding surface plasticity reference model, that is well established for constitutive modeling of cyclic response of sands in the pre-liquefaction state. The roles of the key components of the proposed formulation are examined in a series of sensitivity analyses. Their combined effects in improving the performance of the reference model are examined by simulating undrained cyclic simple shear tests on Ottawa sand, with focus on reproducing the increasing shear strain amplitude as well as its saturation in the post-liquefaction response.     

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

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

基于数字图像相关法的砂土细观直剪试验及其颗粒流数值模拟

砂土等天然颗粒材料具有的剪胀性、压硬性、各向异性等特殊力学性质受控于其内部的微细结构及其演化,如何利用试验与数值模拟等手段对砂土细观组构进行量化,对深入研究砂土变形机制尤为重要。基于数字图像相关方法的光学测量技术与土工试验相结合的新手段,通过对直剪仪的可视化改造,开发了能初步实现土体宏细观力学性状联合测量的“砂土变形细观瞬时光学测量系统”,并利用该系统对福建标准砂进行了室内细观直剪试验。随后进行了基于颗粒流软件PFC2D直剪数值试验,提取室内试验还难以获得的细观信息,作为对室内试验的补充。通过对宏观力学性质和位移场、应变场、颗粒定向、速度场、颗粒配位数和接触力链等细观的室内与数值试验结果进行剖析,探讨了砂土力学性质的细观机制。     

砂土液化大变形模型在FLAC3D中的开发与应用

通过FLAC3D二次开发平台,在VC++环境下实现了基于砂土液化大变形机制的变动映射中心边界面弹塑性模型的二次开发。基于饱和砂土液化大变形本构理论,该模型符合三维应力空间边界面映射规则,并引入临界状态变量,可实现对不同密度和围压状态下砂土液化大变形分析。针对FLAC3D程序混合离散技术与数据调用模式和模型的体积相容性条件,测试分析了将模型采用不同植入方案的计算稳定性,开发了在FLAC3D混合离散技术下不同子区共享映射中心,进入与离开液化状态保持同步的开发方案,并给出模型开发的关键技术与实施方法。利用开发的模型,对饱和砂土开展了不排水/排水/循环三轴试验与不排水循环扭剪试验模拟,及三维地基的动力反应分析。计算结果表明,模型及所开发程序具有很好的模拟与分析砂土液化后大变形的能力。