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

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

单调荷载下砂土变形过程数值模拟及细观机制研究

采用离散单元法的颗粒流理论,模拟了松砂和密砂在单调荷载作用下的变形过程,研究了砂土渐进破坏过程中的宏观力学行为和细观组构参量的演化规律。采用PFC的FISH语言开发了细观组构统计程序,通过记录加载不同时刻试样的细观参量,如配位数、接触法向分布、粒间法向接触力、切向接触力等的演化,分析了砂土变形过程中细观组构变化与宏观力学响应之间的内在联系。应用表征上述量的组构参数研究了砂土的诱发各向异性,探讨了松砂剪缩、密砂剪胀的细观机制。研究成果对于揭示砂土变形的细观机制以及建立砂土的细观力学模型具有重要意义。     

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

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

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

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

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

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

基于细观统计的砂土摩擦特性与破坏机制研究

砂土是一种典型的颗粒材料,其力学特性和变形机制是细观颗粒在荷载作用下的宏观表现,鉴于砂土具有显著的颗粒性和碎散性,使得对其进行细观颗粒层面的研究非常必要。将土颗粒简化为空间椭球体,假设方向角满足正态分布规律,应用细观统计的方法,从颗粒层面对砂土的摩擦特性和破坏机制进行研究。砂土剪切面上的摩擦特性主要是颗粒在剪切荷载作用下,克服垂直于剪切面应力抬升所遇到的阻力。颗粒运动方向与剪切面滑移方向的夹角可视为单个颗粒的摩擦角,颗粒摩擦角满足与方向角同样的正态分布规律;砂土的剪切破坏是颗粒重排列和应力重分布的过程,该过程中一部分颗粒摩擦角不断弱化,而承载颗粒数量逐渐增多,土体内部应力分布渐趋均匀;随着垂直于剪切面荷载的增大,砂土内摩擦角随之弱化减小,内摩擦角与荷载之间具有对应关系,可进行量化计算。     

胶新铁路砂土液化区路基沉降规律研究

地震液化常给人们带来巨大损失,而剪切振动和循环荷载作用下的动力学效应常被认为是地震液化的主要原因,人们对剪切荷载作用下饱和砂土的液化问题进行了较多的研究,而对循环荷载作用下砂土液化的动力学效应研究较少。胶新铁路在DK39+000开始为高地震烈度区,DK283+550~DK283+770分布有地震可液化层,工程修建后列车动荷载的影响将会有诱发砂土液化的可能性。为了研究通车前自然沉降特征和通车后循环荷载作用下的路基沉降变形规律,本文在具体分析了砂土液化的概念和准则判别的基础上,重点分析了砂土液化区路基沉降特征,包括测试断面竖向分层沉降变形特征分析和路基水平位移特征分析。最后在试验的基础上,从理论上给出了循环荷载下砂土的本构关系。     

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

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

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

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

典型巴东组泥质粉砂岩双轴压缩试验细观组构特征

基于泥质粉砂岩室内双轴压缩试验，建立PFC_2D颗粒流数值模型，以此来探究泥质粉砂岩破坏的细观机理。考虑组成泥质粉砂岩试样的颗粒形状，根据电镜扫描图勾选出5种典型颗粒形状，与圆形颗粒一起生成给定孔隙率的稳定数值试样。颗粒间选用平行粘结接触模型，选取弹性模量、泊松比、峰值应力分别对饱和、天然状态的泥质粉砂岩试样进行不同围压下的细观参数标定，然后进行双轴压缩试验模拟，分析试样在双轴压缩试验过程中颗粒法向接触力、切向接触力、配位数、孔隙率等细观组构参数的分布特征和演化规律。试验结果表明:双轴压缩试验的细观参数标定可不考虑抗剪强度指标黏聚力c和内摩擦角ϕ值的影响。破坏前后，试样各方向统计范围内都存在法向接触力和切向接触力。围压的存在影响试样的起始配位数、孔隙率和试样破坏后的稳定配位数和孔隙率，对配位数和孔隙率的改变速率影响很小。试样空间孔隙率的演化在一定程度上反映了试样破坏时内部结构的变迁，能更加直观地反映试样的破坏模式。      

饱和砂土的剪切波速与其抗液化强度关系研究

根据饱和砂土剪切波速与其抗液化强度的相关性原理,利用剪切波速与振动三轴联合实验装置,进行了控制饱和砂土初始剪切波速的振动液化实验,依据实验结果建立了剪切波速与抗液化强度的定量关系。最后用现场勘查数据对此定量关系进行验证,结果表明：该关系式对实际 66 个未液化地点的判别准确率达到 81.2 %;对 108 个实际液化地点的判别准确率达到 62.8 %;平均判别准确率达到 69.5 %。     

考虑组构变化效应的饱和密砂循环本构模型

为了真实地描述饱和密砂在循环加载过程中的变形行为,需要引入考虑剪胀阶段组构变化的宏观参量。在已有的基于状态参量的本构模型基础上,引入反映组构变化的剪胀内变量,简称组构-剪胀内变量z。以相变线PTL作为参考线,采用基于相变的状态参量判断砂土在初始时刻和任意时刻体积变形的变化趋势,并通过z对剪胀比d的影响,考虑反向加载过程中塑性变形的累积,建立了一个针对饱和密砂的循环加载的弹塑性本构模型。该模型根据试验现象将已有模型中的塑性剪切模量区分为首次加载模量与再加载模量,能较好地模拟排水情况下砂土循环加载的胀-缩变化过程。最后,针对密砂的三轴排水情况,利用文中模型进行预测,并把预测结果与试验结果进行比较,结果表明该模型能够总体反映砂土循环加载的变形行为。     

适用于砂土循环加载分析的边界面塑性模型

基于临界状态土力学框架,建立了一个适用于砂土排水循环加载的边界面塑性模型。采用了考虑虚拟峰值应力比的偏应变硬化准则,初始加载阶段应力点位于边界面上,反向加载阶段以历史最大屈服面作为边界面,同时实现了对密砂软化现象的模拟和对历史所受最大应力的记忆。边界面采用修正的椭圆形,引入考虑密度与应力水平的状态相关剪胀函数,采用非相关联流动法则和以应力反向点作为映射中心的径向映射准则。模型仅有10个参数,通过常规三轴试验即可确定,并且使用一套参数可以模拟不同围压、密度的单调和循环加载情况。分别对饱和砂土的单调、循环排水三轴试验进行模拟,结果表明,该模型能够合理地反映饱和砂土排水条件下的应力-应变特性。     

A modification to dense sand dynamic simulation capability of Pastor–Zienkiewicz–Chan model

A modification to the nonlinear Pastor–Zienkiewicz–Chan (PZC) constitutive model without any change in the number of model parameters is introduced in order to simulate stiffness degradation of dense sands at dynamic loading. The PZC model is based on generalized plasticity and was verified by good prediction of liquefaction and undrained behavior of saturated sand. The PZC is a robust model that can predict drained dynamic behavior of sands, especially stiffness increase in loose sand at reloading of dynamic loading. Yet, this model does not show stiffness degradation of dense sand at reloading. The modification is made through modifying the stress memory factor, H _DM, which is multiplied by the plastic modulus, H _L. This modification does not influence reloading behavior of loose sand. The modified PZC model is verified via results of drained cyclic tests. Two cyclic triaxial tests on loose and dense specimens, along with two cyclic plane strain tests on dense sand are utilized for validation. The model simulation shows that the modified PZC model is able to predict the stiffness degradation of dense sand at reloading well.     

K0固结条件下砂土的循环剪切特性试验研究

利用土工静力-动力液压三轴-扭转多功能剪切仪,针对福建标准砂Dr = 30 %,在K0固结条件下,进行不同围压下的循环三轴试验,探讨了K0固结过程和K0固结条件下砂土的孔隙水压力特性与强度特性,并与原有均等固结条件和非均等固结条件下的循环三轴试验进行对比分析。通过分析和比较表明,砂土的K0固结试验过程可以分为两个阶段,第1个阶段是初期的缓慢变形阶段,第2个阶段是后期的稳定发展阶段,围压对静止土压力系数K0的测定有一定影响。不同固结过程对砂土动孔隙水压力的发展有很大的影响,K0固结条件下试样出现明显的剪胀现象。另外,K0固结条件下砂土动强度高于非均等固结条件（Kc = 1.67）下砂土动强度,低于非均等固结条件（Kc =2.94）下砂土动强度。     

Analysis of liquefaction susceptibility of nearly saturated sands

A new constitutive formulation for simulating the behaviour of nearly saturated sands under seismic loads is presented. The formulation is based on combining the Henry's law for dissolution of gas in water, the ideal or perfect gas law and the law of conservation of mass. The effects of transient air dissolution in water on the compressibility of partially saturated soils are also taken into account. The model was calibrated based on numerical simulations of isotropically consolidated cyclic triaxial tests conducted on partially saturated samples of Toyoura sand. A multi‐yield plasticity soil constitutive model implemented in the finite element code DYNAFLOW was used for these numerical simulations. It is shown that the formulation proposed here is able to reasonably predict the soil cyclic undrained behaviour at various degrees of saturation (95% and higher). Copyright © 2006 John Wiley & Sons, Ltd.     

Experimental study of the grading characteristic effect on the liquefaction resistance of various graded sands and gravelly sands

The liquefaction vulnerability of fine- to coarse-graded saturated sand and gravelly sand is evaluated by stress control cyclic triaxial laboratory tests. Cyclic triaxial test on reconstituted specimens is carried out with relative density of 45 %. Cyclic triaxial tests performed by sinusoidal wave form with 1-Hz frequency and with the effective consolidation pressure equal to 100 kPa. Cyclic stress ratio which caused initial liquefaction in 15 cycles recognized as the liquefaction resistance in cyclic triaxial test of present study. The result had been used to find a relationship between liquefaction resistance and the grading characteristics (e.g., coefficient of curvature and coefficient of uniformity) of various graded sands and gravelly sands. Then, it is realized that a relationship between liquefaction resistance and any of the sizes (i.e., D ₁₀, D ₃₀, D ₅₀, D ₆₀) appears to be more practical.     

Constitutive model for cyclic behaviour of cohesionless sands

This paper presents a constitutive model for describing the stress-strain response of sands under cyclic loading. The model, formulated using the critical state theory within the bounding surface plasticity framework, is an upgraded version of an existing model developed for monotonic behaviour of cohesionless sands. With modification of the hardening law, plastic volumetric strain increment and unloading plastic modulus, the original model was modified to simulate cyclic loading. The proposed model was validated against triaxial cyclic loading tests for Fuji River sand, Toyoura sand and Nigata sand. Comparison between the measured and predicted results suggests that the proposed modified model can capture the main features of cohesionless sands under drained and undrained cyclic loading.     

Liquefaction characteristics of silts

Summary The paper presents a summary of the results of cyclical triaxial load testing of samples of silt and silty sands. The paper emphasizes differences in behaviour observed between reconstituted triaxial samples of clean sand, of sand containing 10, 20, 30 and 60% silt, pure silt, and undisturbed samples of silt and silty sand. An important observation is that the mechanisms of deformation for sit are different for reconstituted and undisturbed samples, the undisturbed sample having a specific geological structure which seems to slow down the excess pore water pressure accumulation, but which still results in cyclic deformations regularly increasing from the very beginning of the test and rapidly reaching high levels. The other important observation is that fine-grained noncohesive soils such as silts and silty sands can be as, or even more, susceptible to liquefaction as clean sands. Test results on samples of sand containing 10, 20 or 30% of silt indicate lesser resistance to liquefaction than pure sand samples. The paper shows the difficulty in identifying a representative parameter to compare the behaviour of silts and silty sands with pure sand, and it seems that more research will be needed in this area.