考虑气相硬化的准饱和黏土弹塑性模型

假定准饱和黏土由饱和基质（土骨架和水）和封闭气体组成,采用饱和基质塑性体应变和气体塑性体应变作为硬化参数以反映准饱和黏土的塑性硬化行为,在综合考虑气体溶解系数随温度和水中含盐量变化、预固结应力随封闭气体的变化等关键因素的基础上,基于临界状态土力学理论提出了可以反映封闭气体气压变化对准饱和黏土力学影响的弹塑性本构模型。模型共计10个材料参数,均可通过等向压缩试验和三轴剪切试验获得。通过与已有试验数据对比分析,结果显示本模型可以较好地模拟准饱和黏土在不排水应力路径下的应力-应变关系、孔隙水压力发展以及饱和度的演变规律;所建立模型形式简单,参数易于确定,为准饱和黏土地基中岩土工程问题的设计、计算及灾变控制提供了重要的理论基础。     

基于修正剑桥模型的三轴试验仿真模拟

基于虚拟现实思想,研制了三轴试验的计算机虚拟仿真课件。介绍了使用Authorware平台开发研制的三轴试验的计算机仿真课件以及太沙基三轴固结过程中孔隙水压力计算公式,详细论述了基于修正剑桥模型的三轴试验土样在不排水剪切过程中超孔隙水压力、应力-应变关系、有效应力和总应力的推导过程,并在课件中编程实现。通过利用课件对三轴试验进行模拟操作,不仅能使学生进行系统的试验前的操作训练,而且能够使学生更深刻地掌握土力学的基本原理。     

含气软黏土的不排水抗剪强度计算模型

含气软黏土广泛分布在世界五大洲,大量离散气泡的存在对软黏土的不排水抗剪强度具有强化或损伤作用,然而目前仍缺乏能综合考虑气相强化与气相损伤竞争机制的含气土不排水剪切强度解析表达式。结合临界状态土力学理论,基于笔者提出的含气土屈服函数,建立了v-lnp¢ 空间下临界状态线截距与土体含气特性(u_{w 0},y₀)的数学关系,进而推导了含气软黏土在三轴应力状态下不排水抗剪强度s_u的理论显式表达式。同时,采用所提出的含气软黏土不排水抗剪强度理论计算公式,分别对比了3种典型的含气软黏土在不同含气特性下的不排水抗剪强度理论预测与试验结果,包括马来西亚高岭土、Combwich 黏土和高岭黏土。另外,对比分析了4种不同的含气软黏土不排水抗剪强度理论计算模型的预测误差百分比,证明了提出的基于改进的含气软黏土本构模型的不排水抗剪强度理论计算式的合理性,能够更加准确地同时考虑含气软黏土的气相强化与气相损伤双重效应。     

考虑循环载荷下饱和黏土软化的损伤边界面模型研究

研究表明,循环载荷作用下饱和黏土将发生软化,其机制主要有两个：一是孔压的积累;二是土体原有结构的不断损伤和新结构的不断重塑。针对上述机制,基于广义各向同性硬化准则建立了考虑饱和黏土循环软化的损伤单面模型。该模型在有效应力空间中引入损伤变量,表征土体结构的损伤和重塑程度,在连续的循环加载下,损伤不断累积,边界面则随着损伤的累积不断收缩,以模拟饱和黏土刚度和强度的软化;以应力反向点作为边界面的广义各向同性硬化中心和映射法则的映射中心,灵活地选择塑性模量的插值公式以模拟塑性变形和孔压的累积以及应力-应变的滞回特性。应用该模型对不排水循环三轴试验进行模拟,并且考查了循环周次、循环应力水平和固结历史对饱和黏土循环软化特性的影响,并与相关试验比较,验证了模型的有效性。     

DEM study on the undrained mechanical behavior of gassy sand

Exploitation of gas hydrate requires an improved understanding of the behavior of gassy soil. In this study, a numerical simulation technique is proposed that uses the discrete element method (DEM) to evaluate the undrained mechanical behavior of gassy sand, which consists of three phases—sand particles, water, and gas. The governing equation of gassy sand under undrained condition is derived and incorporated into the DEM simulation, in which gas solution and exsolution are considered. Undrained isotropic unloading and undrained triaxial compression tests on gassy sand samples with dissolved methane and carbon dioxide are simulated. Reasonable consistency is found when the simulation results are compared with experimental observation and theoretical calculation. The simulation highlights the significant influence of the gas solubility on the undrained behavior of gassy sand. The micro-responses of the observed phenomenon are also investigated.

     

饱和粉土静态液化性能试验研究

通过固结不排水三轴压缩试验,分析了围压、固结比和干密度等因素对饱和粉土静态液化特性的影响。试验结果表明,在干密度较小时,饱和粉土的偏应力-应变曲线呈现明显的硬变软化型,随轴向应变增大超静孔隙水压力增加、有效应力减小而发生静态液化,当干密度达到1.58 g· cm-3时,饱和粉土的偏应力-应变曲线表现出硬变硬化现象,超静孔隙水压力为负值或接近0,饱和粉土不再发生静态液化,即饱和粉土存在静态液化的干密度临界值;其他条件不变,随着围压、固结比或干密度的增大,偏应力峰值和残余强度均增大,静态液化势降低;根据有效应力路径建立了流滑面以作为饱和粉土稳定区与非稳定区的分界面。     

Cone penetration-induced pore pressure distribution and dissipation

The excess pore water pressure distribution (u) induced by the penetration of a piezocone into clay and its dissipation behaviour have been investigated by laboratory model tests, theoretical analysis and numerical simulation. Based on the results of the tests and the analysis, a semi-theoretical method has been proposed to predict the piezocone penetration-induced pore pressure distribution in the radial direction from the shoulder of the cone. The method can consider the effect of the undrained shear strength (s_u), over-consolidation ratio (OCR) and rigidity index (I_r) of the soil. With a reliably predicted initial distribution of u and the measured curve of dissipation of pore water pressure at the shoulder of the cone (u₂), the coefficient of consolidation of the soil in the horizontal direction (c_h) can be back-fitted by analysis of the pore pressure dissipation. Comparing the back-fitted values of c_h with the values directly estimated by a previously proposed method indicates that the previously proposed method can be used reliably to estimate c_h values from non-standard dissipation curves (where u₂ increases initially and then dissipates with time).     

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

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

SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

不排水条件下砾石土的应变率效应

为研究应变率对砾石土应力、应变、强度特性的影响和应变率效应产生机制,对砾石和黏土混合备制的饱和砾石土心墙料,在应变率0.028～0.690 %/min范围内进行了不同围压下的固结不排水三轴压缩试验,定量地分析了三轴不排水抗剪强度的应变率效应。试验结果表明,不同的应变率下砾石土应力-应变曲线具有相似性;应变率增大10倍,不排水抗剪强度增大7.4%～29.9%;砾石土的应变率效应是剪切产生的孔隙水压力和土料本身黏滞性共同作用的结果,在不同的围压下二者发挥的程度不同     

Analysis of undrained cavity expansion in elasto‐plastic soils with non‐linear elasticity

A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd.     

冲击作用下饱和土性状的试验研究

在三轴条件下,对饱和土（砂土和黏土）进行排水与不排水条件下的冲击试验及冲击后再固结试验,对比研究了不同渗透性土在不同排水条件下的冲击动力响应和冲击后再固结性状。结果表明：饱和黏土不排水冲击时的孔隙水压力随冲击击数增加而升高并逐渐稳定,排水冲击时的孔隙水压力则是先达到峰值然后有所下降;砂土不排水冲击时的冲击能量对孔隙水压力影响最明显;饱和砂土不排水冲击时的轴向应变与冲击击数呈近似线性关系,饱和黏土冲击及饱和砂土排水冲击则呈近二次曲线关系;饱和砂土不排水冲击后再固结阶段的孔隙水压力立即消散为0,同时体变迅速增大到一定值;饱和黏土在冲击后再固结阶段的孔隙水压力在一定时间内逐渐消散完毕,同时体变逐渐增大;饱和黏土排水冲击时,冲击阶段产生的体变占冲击引起总体变的39%～49%,冲击后再固结阶段产生的体变占51%～61%;砂土和黏土的总体变均表现为排水冲击明显大于不排水冲击,改善冲击时的排水条件有利于提高加固效果。     

基于以气驱水技术的含水地层灌浆改性试验研究

当对含水地层进行灌浆改性时,通常会受到土体孔隙中自由水的不利影响。为此,考虑将以气驱水技术应用于含水地层化学灌浆中,以提高改性效果。对饱和黏土与粉质砂土开展以气驱水可行性试验,结果显示在饱和粉质砂土中充气驱替孔隙中自由水具有较好的可行性,在饱和黏土中则难以顺利进行。这是因为土体内部孔隙大小会直接影响初始进气值Pc,同时充气压力亦不能超过极限压力Pmax,否则土体内部会发生劈裂破坏,形成大范围空气通道;研制含水地层以气驱水化学灌浆模拟装置,以粉质砂土为研究对象进行对比试验,试验结果显示在以气驱水技术的辅助下进行化学灌浆可更加显著地改善原有土体的渗透性、无侧限抗压强度与微观结构,电子显微镜扫描也显示出相似结果,这是因为通过充气能够在含水地层中产生空气隔幕,为浆液提供良好的灌注、凝硬环境。     

Sensitive bounding surface constitutive model for structured clay

The main purpose of the paper is to present a relatively simple, yet realistic, constitutive model for simulations of structured sensitive clays. The proposed constitutive model can simulate 1‐D and isotropic consolidation, and drained and undrained shear response of sensitive structured clay. The proposed sensitive bounding surface model is based on concepts from the modified Cam clay model 8 and bounding surface plasticity 27 , with the addition of a simple degradation law. The key material parameters are M, λ, κ, and ν from the modified Cam clay framework, h from the bounding surface framework to model a smoothed elasto‐plastic transition, and ω_v, ω_q, and S_sr to model softening associated with destructuration. The model has separate parameters to model destructuration caused by volumetric strain and deviatoric strain. The model is capable of modeling unusual behavior of strain softening during 1‐D compression (i.e., a reduction of effective stress as void ratio decreases). A good match between test results and the model simulation is demonstrated. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling thermo-elastic–viscoplastic behaviour of marine clay

Marine clay supporting high-temperature offshore structure susceptible to random movements (such as a buckling high-temperature pipeline) experiences variable shearing rates at an elevated temperature that is higher than the marine environment (typically 4 °C). This practically implies that the undrained shear strength (s_u) of the marine clay being routinely characterized in situ by penetrometers at a constant rate under an isothermal condition (4 °C) should be carefully corrected, by taking into account the temperature and rate dependency. To date, the combined effects of rate and temperature on the undrained shear behaviour of marine clay are merely investigated experimentally and theoretically. This study presents the development of an anisotropic thermo-elastic–viscoplastic model and a series of temperature- and rate-controlled triaxial tests for validation purpose. Compared to the modified Cam-Clay model, the proposed model only introduces three new parameters to characterize the temperature dependency, rate dependency and the inherent anisotropy of K₀-consolidated marine clay. The predictive capability of the model has been validated by the triaxial test results. Based on the new model, an explicit equation is formulated for quantifying the temperature- and rate-dependent s_u of marine clay. Calculation charts are also developed to quantify s_u of marine clay with different plasticity indexes under various strain rates and temperatures.

     

描述饱和黏性土应变软化的弹塑性双面模型

常规三轴压缩试验中具有较强结构性的黏性土在围压较低时其应力−应变关系会呈现应变软化现象,一般还伴有塑性变形,通常土体内部结构损伤是应变软化产生的主要原因。考虑到采用经典塑性理论描述材料的应变软化不仅会违背 Drucker 的稳定性假设,而且也不能描述卸载塑性。因此,基于修正剑桥模型及 Li 和 Meissner 提出的塑性硬化准则,建立了一个描述饱和黏性土不排水应变软化的弹塑性双面模型。该模型以应力−应变曲线的峰值点分界,将应变硬化和应变软化分别作为独立的加载事件进行分析,同时引入新的结构性参数表征剪切过程中土体结构损伤导致的塑性刚度衰退。对不同固结状态饱和结构性黏土的三轴固结不排水压缩试验结果的模拟表明,所建模型能够较好地描述饱和黏性土的不排水应变软化特性。     

长期循环荷载作用下排水条件对饱和软黏土动力特性影响

由于软黏土渗透系数很小,目前针对饱和软黏土的循环荷载试验绝大多数都是在不排水条件下进行的。但对交通荷载而言,软黏土承受经年累月的长期循环荷载作用,实际的排水条件应该是部分排水的。基于此,通过对温州原状饱和软黏土进行不同循环应力比下的不排水和部分排水大周数（50 000次）循环三轴试验,分析了长期循环荷载作用下排水条件对饱和软黏土动力特性的影响。研究结果表明与不排水条件相比,部分排水条件下软黏土试样的动力特性表现出很大的不同。随着循环次数的增加,部分排水条件下的孔压随先增加后降低,存在一个峰值;经过较大的循环次数后,回弹应变逐渐减小,应力-应变滞回圈逐渐缩短,面积也明显减小。为了准确预测交通荷载作用下的长期沉降,采用部分排水条件进行试验是十分必要的。     

Monotonic and cyclic characteristics of K0-Consolidated saturated soft clay under a stress path involving a variable confining pressure

In certain field conditions such as offshore projects under wave loads or embankments under traffic loads, both the vertical and horizontal stresses are variable. However, previous investigations rarely considered the variation in horizontal stress. To better understand the characteristics of natural saturated soft clay, a series of monotonic and cyclic triaxial tests with a K₀-consolidation state were carried out under a variable confining pressure (VCP) stress path. The development of axial strain, pore water pressure and effective stress path is analysed. The results show that with the increase in η (the ratio of the variation in the mean effective principal stress to that of the deviatoric stress), the undrained shear strength (q_f) decreases continuously. The pore water pressure generation is slightly improved under a stress path with increasing confining pressure. Based on the test results, a unified formula was established to predict the pore water pressure under VCP stress paths. The unique p–q–e relationship of normally consolidated clay in monotonic VCP triaxial tests was also demonstrated. Under VCP stress paths, the amplitude of the pore pressure increases, and the effective stress path tilts more sharply to the right. Moreover, a unified formula was established that can provide a good reference for predicting effective stress paths under cyclic VCP triaxial tests.

     

A multisurface kinematic hardening model for the behavior of clays under combined static and undrained cyclic loading

In dynamic geotechnical problems, soils are often subjected to a combination of sustained static and fast cyclic loading. Under such loading conditions, saturated and normally consolidated clays generally experience a build-up of excess pore water pressure along with a degradation of stiffness and strength. If the strength of the soil falls below the static stress demand, a self-driven failure is triggered. In this paper, a constitutive model is presented for the analysis of such problems, based on a general multisurface plasticity framework. The hardening behavior, the initial arrangement of the surfaces, and the nonassociated volumetric flow rule are defined to capture important aspects of cyclic clay behavior. This includes nonlinear hysteretic stress-strain behavior, the effect of anisotropic consolidation, and the generation of excess pore water pressure during undrained cyclic loading along with a degradation of stiffness and strength. The model requires nine independent parameters, which can be derived from standard laboratory tests. A customized experimental program has been performed to validate the model performance. The model predictions show a good agreement with test results from monotonic and cyclic undrained triaxial tests, in particular with respect to the strain-softening response and the number of loading cycles to failure. A procedure for a general stress-space implicit numerical implementation for undrained, total stress-based finite element analyses is presented, including the derivation of the consistent tangent operator. Finally, a simulation of the seismic response of a submarine slope is shown to illustrate a possible application of the presented model.     

Anisotropic viscoplastic modeling of rate‐dependent behavior of clay

The objective of this study is to derive an effective stress‐based constitutive law capable of predicting rate‐dependent stress–strain, stress path and undrained shear strength and creep behavior. The flow rule used in the MIT‐E3 model and viscoplasticity theory is employed in the derivation. The model adopts the yield surface capable of representing the yield behavior of the Taipei silty clay and assumes that it is initially symmetric about the K₀‐line. A method is then developed to compute the gyration and expansion of the loading surface to simulate the anisotropic behavior due to the principal stress rotation after shear. There are 11 parameters required for the model to describe the soil behavior and six of them are exactly the same as those used in the Modified Cam‐clay model. The five additional parameters can be obtained by parametric studies or conventional soil tests, such as consolidation tests, triaxial compression and extension tests. Finally, verification of the model for the anisotropic behavior, creep behavior and the rate‐dependent undrained stress–strain and shear strength of the Taipei silty clay is conducted. Copyright © 2010 John Wiley & Sons, Ltd.