减饱和松砂静态液化的水-气两相流耦合分析

减饱和法是一种通过减小饱和砂土地基中的饱和度,提高地基抗液化强度的新方法。基于水-气两相流反应与土体骨架变形的耦合模拟方法,建立了单调加载条件下减饱和砂土静态液化的数值分析模型。开展了减饱和松砂的三轴不排水试验数值模拟研究,与室内试验结果对比,发现两相流模型能够准确描述减饱和砂土加载过程中的应力-应变关系、应力路径及孔隙水压力增长规律,验证了两相流模拟方法的正确性。数值分析结果还表明,加载过程中减饱和松砂中的饱和度会增加,直至达到一个稳定值,当围压一定时,减饱和松砂加载结束时的饱和度与初始饱和度呈线性关系;且砂土中气体会在荷载作用下被压缩,使得减饱和砂土在不排水条件下发生剪缩。计算发现,当砂土饱和度从100%减小到94.5%时,孔隙水压力系数B值会减小约80%,最大孔隙压力值会降低40%~50%,不排水剪切强度提高2.0~2.5倍,残余强度会提升10倍以上,由此可知,此为减饱和法抗液化的主要机制,而基质吸力不是减饱和法提高砂土抗剪强度的主要原因。     

减饱和砂土缓倾场地的液化性状分析

减饱和法是一种通过减小砂土地基的饱和度,从而提高地基抗液化强度的新方法。基于减饱和砂土中流体模量同步更新的改进算法对减饱和砂土离心机振动台试验进行了数值分析,并与单一流体模量的简化算法进行了对比分析。结果表明：由于改进算法中考虑了因孔压变化引起的等效流体模量的变化,其计算结果更接近试验结果,而简化算法低估了减饱和砂土的孔压积累。基于改进算法开展了不同饱和度、倾斜角度的缓倾场地上液化变形的数值模拟研究,分析发现超孔隙水压力增长的速度及其峰值随着饱和度的增加而增大,饱和度从100%降低至96.4%,同一深度处的超孔压峰值降低约20%～65%,加速度响应的峰值也有明显的降低;沿地基深度0.75 m到9.00 m,侧向位移减少约20%～50%,表明饱和度的降低对抑制倾斜场地上可液化砂土层的侧向变形有显著效果,随着地基深度的减小,饱和度对于侧向位移的影响越来越明显。     

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

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

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

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

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

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

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

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

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

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

夯扩挤密碎石桩加固液化砂土地基的动力数值分析

砂土液化问题一直是土动力学与岩土地震工程研究领域的重要课题之一。基于南水北调中线某工程,通过现场和室内试验获取土体的物理力学参数,利用岩土数值分析软件FLAC3D对夯扩挤密碎石桩加固干渠液化砂土地基进行了动力数值分析。结果表明,由于夯扩挤密碎石桩的排水作用,干渠底部饱和砂土地基中的超静孔隙水压力和孔压比与加固前相比明显减小;干渠渠道底部饱和砂土中的监测曲线表明,随着地震荷载持续时间的增加,饱和砂土地基中超静孔隙水压力和孔压比峰值较加固前大幅值降低,且时程曲线达到峰值之后也由加固前的基本保持不变改为迅速消减降低;由于夯扩挤密碎石桩的排水和挤密作用,有效消除了干渠渠道底部以及渠堤坡面外侧平台至坡脚底部砂土层的液化现象,加固后干渠底部饱和砂土地基中没有液化现象产生。     

剪胀性砂土边界面模型的研究

剪胀性对于砂土,尤其是中密以及密实砂土,是一个非常显著的特性。相变线是剪胀性砂土的特征曲线,能够反映砂土的围压以及初时孔隙比对变形特性的影响。本文在边界面塑性理论的框架内,把相变状态参量引入到剪胀方程以及塑性硬化模量中,建立了一个能够描述砂土剪胀性以及循环特性的本构模型。本模型采用一套参量可以模拟不同初时孔隙比、不同围压、排水(或不排水)条件下单调(或循环)加载的应力-应变特性。验证表明本模型数值计算与试验结果相吻合。     

安哥拉Quelo砂场地非饱和渗流试验与计算

在安哥拉某Quelo砂土场地开展了现场试坑浸水试验,实测了水在地基土中的渗透和消散过程;通过室内测试得到砂土的主要物理参数,采用改进的Kovács粒径方法估计土-水特征曲线,Van Genuchen模型估算渗透系数曲线,数值模拟了与现场试验情况相同的水渗透和消散过程;对比了实测结果和数值计算结果的异同。研究结果表明浸水条件下自由水面之下的砂土远未达到完全饱和,数值计算结果虽能较好反映浸水渗透的宏观规律以及停水后水的消散过程,但要较准确预测水的浸润范围,需要建立经验的“饱和体积含水率”取值方法,本砂土场地“饱和体积含水率”取对应饱和度65%的体积含水率较为适宜。     

常剪应力路径下含气砂土的三轴试验

天然气水合物完全分解时,产生的气体使得能源土孔隙压力急速增加,有效应力减小,进而引起土体液化破坏。此时深海能源土斜坡的应力状态与静力液化失稳过程可简化为含气土在常剪应力排水(或不排水)应力路径下的破坏问题。以此为背景,提出了制备含气砂土试样的改进充气管法,并开展了含气砂土的常剪应力路径三轴试验。22组试验结果表明:同一孔隙比的含气密砂在不同围压与常剪应力下具有相同的失稳线;含气砂土试样失稳时的应力比和体变均随初始相对密实度的增大而增大;含气密砂在常剪应力路径下饱和度对失稳特征影响的规律性在排水与不排水条件下均不明显,但在不排水条件下含气砂土的孔压(或体变)对变形的敏感性降低;含气密砂在常剪应力路径到达失稳点之后,排水条件下是瞬变的液化鼓胀破坏,不排水条件下是渐变的剪切破坏。     

A hypoplastic constitutive model for debris materials

Debris flow is a very common and destructive natural hazard in mountainous regions. Pore water pressure is the major triggering factor in the initiation of debris flow. Excessive pore water pressure is also observed during the runout and deposition of debris flow. Debris materials are normally treated as solid particle–viscous fluid mixture in the constitutive modeling. A suitable constitutive model which can capture the solid-like and fluid-like behavior of solid–fluid mixture should have the capability to describe the developing of pore water pressure (or effective stresses) in the initiation stage and determine the residual effective stresses exactly. In this paper, a constitutive model of debris materials is developed based on a framework where a static portion for the frictional behavior and a dynamic portion for the viscous behavior are combined. The frictional behavior is described by a hypoplastic model with critical state for granular materials. The model performance is demonstrated by simulating undrained simple shear tests of saturated sand, which are particularly relevant for the initiation of debris flows. The partial and full liquefaction of saturated granular material under undrained condition is reproduced by the hypoplastic model. The viscous behavior is described by the tensor form of a modified Bagnold’s theory for solid–fluid suspension, in which the drag force of the interstitial fluid and the particle collisions are considered. The complete model by combining the static and dynamic parts is used to simulate two annular shear tests. The predicted residual strength in the quasi-static stage combined with the stresses in the flowing stage agrees well with the experimental data. The non-quadratic dependence between the stresses and the shear rate in the slow shear stage for the relatively dense specimens is captured.     

饱和土体静态液化失稳理论预测

通过在屈服准则和剪胀性方程中引入材料状态参数,建立了一个与材料状态相关的本构模型,用于模拟不同初始孔隙比和围压下土体的应力-应变关系。基于二阶功理论,建立了饱和土体静态液化失稳理论判别准则。通过预测一系列饱和松砂三轴不排水试验结果,验证了所建立的本构模型及判别预测准则的合理性。分析结果表明,静态液化发生于土体硬化阶段,静态液化触发伴随着剪应力达到峰值,其后剪应力降低且孔隙水压力持续增长。进一步地,模拟了充分降解的城市固体废弃物在不排水条件下的应力-应变特性,并预测到了潜在失稳线及静态液化触发点。     

Pore pressure tension cut off in two-phase finite element analyses

In the coupled finite element analysis of the two-phase, saturated soils the pore pressure could decrease below the atmospheric pressure in particular when the plastic dilation of the skeleton is accounted for. The consequent reduction of the bulk modulus of the liquid phase is here modelled in a simplified manner through a tension cut off on the pore pressure. Two alternative approaches are outlined in the framework of the finite element method. The first one introduces a unilateral constraint on the nodal pore pressures and involves manipulations on matrices of the assembled mesh, whilst the second approach requires calculations solely at the element level. The effects of the pore pressure tension cut off are illustrated through some undrained test examples in elastic–plastic regime.     

Effect of fines content and void ratio on the saturated hydraulic conductivity and undrained shear strength of sand–silt mixtures

The hydraulic conductivity represents an important indicator parameter in the generation and redistribution of excess pore pressure of sand–silt mixture soil deposits during earthquakes. This paper aims to determine the relationship between the undrained shear strength (liquefaction resistance) and the saturated hydraulic conductivity of the sand–silt mixtures and how much they are affected by the percentage of low plastic fines (finer than 0.074 mm) and void ratio of the soil. The results of flexible wall permeameter and undrained monotonic triaxial tests carried out on samples reconstituted from Chlef river sand with 0, 10, 20, 30, 40, and 50 % non-plastic silt at an effective confining pressure of 100 kPa and two initial relative densities (D _r = 20, 91 %) are presented and discussed. It was found that the undrained shear strength (liquefaction resistance) can be correlated to the fines content, intergranular void ratio and saturated hydraulic conductivity. The results obtained from this study reveal that the saturated hydraulic conductivity (k _sat) of the sand mixed with 50 % low plastic fines can be, in average, four orders of magnitude smaller than that of the clean sand. The results show also that the global void ratio could not be used as a pertinent parameter to explain the undrained shear strength and saturated hydraulic conductivity response of the sand–silt mixtures.     

循环荷载下软黏土不排水强度弱化分析的动力有限元ABAQUS实现

在波浪循环荷载作用下饱和软黏土地基出现孔隙水压力升高,并导致不排水强度弱化,严重影响防波堤的承载性能。考虑静偏应力影响,基于最大孔隙水压力发展模型和正常固结软黏土不排水强度公式,推导出软黏土不排水强度随循环荷载作用次数和应力水平变化的动态折减规律。结合软黏土不排水强度动态折减规律和M-C屈服准则,在有限元软件ABAQUS上实现软黏土不排水强度循环弱化分析的数值开发和动力运算过程。运用该动力有限元方法对天津港防波堤地基软黏土的动、静三轴试验进行数值模拟运算。结果表明,最大孔隙水压力发展曲线以及循环荷载作用后不排水强度的数值预测结果与动三轴试验结果吻合良好。另外,动力有限元方法（DFEM）能够表示土体强度在循环荷载作用下的具体弱化过程。     

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

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

A continuum‐discrete model using Darcy's law: formulation and verification

This paper presents a numerical scheme for fluid‐particle coupled discrete element method (DEM), which is based on poro‐elasticity. The motion of the particles is resolved by means of DEM. While within the proposition of Darcian regime, the fluid is assumed as a continuum phase on a Eulerian mesh, and the continuity equation on the fluid mesh for a compressible fluid is solved using the FEM. Analytical solutions of traditional soil mechanics examples, such as the isotropic compression and one‐dimensional upward seepage flow, were used to validate the proposed algorithm quantitatively. The numerical results showed very good agreement with the analytical solutions, which show the correctness of this algorithm. Sensitivity studies on the effect of some influential factors of the coupling scheme such as pore fluid bulk modulus, volumetric strain calculation, and fluid mesh size were performed to display the accuracy, efficiency, and robustness of the numerical algorithm. It is revealed that the pore fluid bulk modulus is a critical parameter that can affect the accuracy of the results. Because of the iterative coupling scheme of these algorithms, high value of fluid bulk modulus can result in instability and consequently reduction in the maximum possible time‐step. Furthermore, the increase of the fluid mesh size reduces the accuracy of the calculated pore pressure. This study enhances our current understanding of the capacity of fluid‐particle coupled DEM to simulate the mechanical behavior of saturated granular materials. Copyright © 2014 John Wiley & Sons, Ltd.