采用环剪仪对超固结黏土抗剪强度特性的研究

利用大型高速环剪仪,对不同超固结比、法向应力和剪切速率下的超固结饱和黏土的峰值强度和残余强度特性进行了研究,并对循环荷载作用前后残余强度的变化进行了初步探讨。试验结果表明,（1）超固结比对超固结黏土的峰值强度和残余强度有着明显的影响;（2）在剪切速率相同的条件下,土体达到残余强度时的位移取决于现存的应力状态,而与应力历史无关;（3）剪切速率越大,峰值强度随之增大,达到稳定残余强度时的剪切位移也随之增加,但剪切速率的变化对残余强度值几乎没有影响;（4）在循环荷载作用下残余强度不同程度的降低,最大降低幅度达12.2 %;当土体剪切面为不规则剪切带时,施加循环荷载后出现残余强度上升的现象。     

桩柱结构在地震作用下非线性响应研究

基于试验基础上建立的经典弹塑性模型--剑桥模型能够准确描述正常固结土的应力-应变关系。当土体的应力历史上经历过卸载或受到循环交变荷载作用即进入超固结状态,它作为土的应力历史的反映,相比正常固结土受力特性有着显著的差异。为研究超固结因素对土体加载特性的影响,在引入能考虑超固结状态影响的下负荷面剑桥模型后,通过三轴压缩和剪切试验对处于超固结状态下土体的受力特性进行了对比分析,并对循环剪切加载下的应力-应变关系以及超固结比的演化规律进行了研究。结果表明,下负荷面剑桥模型能准确反映超固结因素对土体力学特性的影响,相比原状土有着更高的屈服强度。而通过数值模拟自由场地基在地震作用下的动力响应可以看出,超固结因素对地基的动力响应起到了不可忽略的影响,尤其在强震下更需要考虑其影响。在自由场地基地震动力响应基础上,通过对桩柱结构桩-土耦合系统在地震作用下非线性动力响应的模拟对土体非线性以及超固结因素的影响进行了对比研究,研究表明：土体的非线性因素能显著降低结构振动响应中的高频成分,由于土体在交变加载下很快进入超固结状态,相对于剑桥模型,下负荷面剑桥模型在考虑超固结因素后土体的承载性能显著提高,尤其在强震作用下超固结因素带来的影响更加明显,因此,建议对桩基结构物地震响应研究考虑超固结因素影响,以提高桩基结构物地震响应模拟的精确度和可靠性。     

土体各向异性边界面模型隐式算法在ABAQUS软件中的实现

基于Wheeler土体各向异性旋转硬化法则,结合边界面理论,构造一个能够反映土体初始各向异性及加载后应力诱发各向异性的边界面本构模型,并借助ABAQUS软件提供的UMAT子程序接口,采用隐式积分算法--图形返回算法实现。通过对正常固结状态下（OCR=1）高岭土试样三轴不排水剪切试验进行模拟,并将模拟结果与ABAQUS自带的修正剑桥模型模拟结果进行了比较分析,表明本模型的模拟结果能够反映土体在偏压加载过程中产生的各向异性现象。在此基础上,采用本模型对中等超固结（OCR=4）高岭土试样三轴不排水剪切试验进行模拟,并再次与ABAQUS自带的修正剑桥模型模拟结果进行比较,表明本模型能够较好地反映中等超固结土在小应变情况下的非线性特性。相比于经典弹塑性模型,如修正剑桥模型,本模型的模拟结果更符合中等超固结土的变形特性。     

不同固结状态下黏土抗剪强度与剪切速率的关系

利用直剪仪针对不同固结状态下的饱和黏性土,在不同剪切速率条件下进行了系统的试验研究,探讨了剪切速率与剪应力-剪切位移曲线变化规律,分析了剪切速率对强度参数的影响及原因。试验结果表明:对于正常固结土,剪切速率越大,峰值剪应力越小,内摩擦角越小。不同的固结状态下,剪切速率对抗剪强度参数影响不同。超固结比为2时,内摩擦角随着剪切速率的增大而减小,黏聚力随着剪切速率的增大而增大;超固结比为3时,内摩擦角随剪切速率的影响较小,黏聚力随着剪切速率的增大而增大。另外,从黏聚力和内摩擦力的角度,分析了不同剪切速率条件下土体抗剪强度变化的主要控制因素。最后,从孔隙水压力的角度分析了不同剪切速率对抗剪强度的影响。在相同的法向应力下,对于正常固结土,不同剪切速率引起的剪切带周围孔隙水压力变化量与破坏剪应力变化量成正比关系;对于超固结土,黏聚力变化量减去破坏剪应力变化量的差值与孔隙水压力的增量成正比。     

上海粘性土剪切带倾角的试验研究

采用可量测局部侧向应变的平面应变仪，对上海地区第（2）层褐黄色粉质粘土和第（6）层暗绿色粉质粘土的原状土样，进行固结不排水平面应变压缩试验，观察土样的局部化变形过程，量测剪切带倾角，研究剪切带倾角随固结压力的变化情况，计算剪切带倾角的传统土力学理论解并与实测值进行比较，根据剪切带内局部剪胀角的近似估算，初步研究剪切带内的变形机理，认为剪切带的形成与孔隙水的运动密切相关。     

基于黏弹性理论的沉桩再固结问题分析

假设桩周土体为饱和黏弹性介质,采用Burgers流变模型进行描述,同时考虑竖向和径向固结,建立了固结控制方程。根据不排水和自由排水情况,将边界条件分为3类并分别得到超孔隙水压力消散的级数解答,该解答能够为孔压静力触探反求固结系数提供一定的理论依据。在此基础上编制了应用程序,对Burgers流变模型中主要参数进行了分析。结果表明,地基表面自由排水、桩端地基不排水条件下,在一定深度以内的桩周土体的固结速度随深度降低,但超过某一范围后固结速度趋于稳定;上、下边界均自由排水条件下,固结速度随深度增加呈现下降、稳定、升高;上、下边界均不排水条件时,孔压消散速度不随深度变化,可简化为本解答仅考虑径向固结的特例。同时土体的流变特性对超孔隙水压力消散的影响比较显著,流变参数G1/1的变化使超孔隙水压力趋于某不为0的定值,且该值随G1/1比值的增大而增加;其他参数不变时,土体剪切刚度比G1/G2的增大会引起孔压消散速度的下降。     

用快剪试验推求先期固结压力

土体先期固结压力是判断土体应力历史的关键性指标。依据Wroth在一维固结条件下提出的不排水强度与超固结比（OCR）的关系式,利用快速直剪的剪切强度计算临界状态孔隙水压力系数Λ₀,从而计算出OCR以及先期固结压力。对南京某黏土的重塑饱和黏土的试验分析结果表明,快剪试验确定的超固结比与设计超固结比基本一致,因此可作为求解土体先期固结压力的一个途径。      

超固结粘土的二元介质模型

为了分析超固结粘土边坡的变形和稳定的需要,在岩土破损力学的框架内建议了一个适用于这类土的二元介质模型。该模型考虑了土体内在的不均一性,把它看作由结构块和结构带组成的复合体,两者共同分担外荷载,模型包含10参数,并拟定了这些参数的测定方法。通过模拟三轴试验的计算表明,该模型可以反映London粘土的应力-应变特征。     

不同应力路径下超固结黏土试样变形局部化分析

基于改进伏斯列夫面超固结黏土三维本构模型,利用有限元软件ABAQUS材料子程序接口,采用回映应力更新算法,实现了该模型在有限元分析中的应用。通过该模型与比奥固结理论的耦合,对超固结比为8的超固结黏土在三轴压缩、三轴伸长及平面应变应力条件下的变形局部化问题,进行了水-土耦合弹塑性有限元分析。分析结果表明：剪切带带内、带外点经历不同应力路径;剪切带带外单元经历了体缩、剪胀及被吸水体缩过程,而剪切带带内单元一直保持剪胀趋势;剪切带的形成伴随着剪胀,剪切带内、外出现了负的孔压,且孔压的分布也具有局部化特性。关于剪切带带内、外的孔隙水压及体变变化趋势与剪切速率有关,而平面应变介于三轴压缩与三轴伸长之间,但平面应变较早出现剪切带。孔隙水的迁移速度影响剪切带带内单元的剪胀,进而影响剪切带的形成及发展;而围压和弱单元位置也对剪切带的形成也有影响。     

非等温分布条件下考虑半透水边界时饱和黏土的一维固结解析解

温度的变化会导致土体的物理力学性质改变,且在一些实际工程中,饱和黏土会处于非等温分布状态。为此,针对非等温分布条件下饱和黏土的一维固结问题,考虑了更具普遍性的半透水边界,通过某些假定推导了单级线性加荷形式下饱和黏土一维固结控制方程,并利用分离变量法求解得到了控制方程的解析解。通过将所提解析解分别与已有解析解和有限差分解展开对比分析,验证了所提解答的正确性。基于所提解析解,利用某一算例分析了温度梯度、半透水边界参数及加荷时间对固结性状的影响。结果表明：温度梯度 M 越大,土体的渗透性越大,土体的固结速率越快;半透水边界参数 R1和 R2越大,相同时间内土体的超孔隙水压力越小,土体的平均固结度越大;土体的平均固结度随加荷时间 tc 的增大而减小,这主要是由于加荷阶段所施加的外荷载小于最终荷载,但加荷时间tc的延长可一定程度减小土体中产生的最大超孔隙水压力。     

不同应力路径下剪切带的数值模拟

采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明：超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。     

Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model

A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd.     

重塑超固结上海软土力学特性及弹塑性模拟

对典型上海软土重塑样进行了围压不变和平均主应力不变的三轴排水剪切试验,得到重塑上海软土在不同初始超固结比和围压条件下的应力-应变关系,弄清了超固结比、围压以及应力路径对重塑上海软土的变形和强度特性的影响;根据土体的应力-应变曲线得到重塑上海软土的临界状态应力比及内摩擦角。采用姚仰平等建议的基于伏斯列夫面的超固结土本构模型,并根据等向压缩及三轴排水剪切试验确定其模型参数,对保持围压和平均主应力不变的三轴压缩试验进行了模型预测。预测结果表明,此超固结土本构模型能较好地反映重塑超固结上海软土的变形和强度特性。     

Prediction of stresses and strains around model tunnels with adjacent embedded walls in overconsolidated clay

This paper presents the results of finite element analyses carried out using different constitutive models for overconsolidated clay: the Modified Cam clay model and the Three-Surface Kinematic Hardening (3-SKH) model. These analyses are evaluated against data from an extensive series of physical model tests examining the influence of an embedded wall placed near a tunnel on ground movements and tunnel stability. It is shown that for heavily overconsolidated soils reasonable predictions of both deformations and failure can be obtained from kinematic hardening models such as the 3-SKH model, which allow plastic deformation inside a Modified Cam clay state boundary surface.     

Experimental study of strain localization in sensitive clays

For evaluation of slope stability in materials displaying strain-softening behavior, knowledge concerning the failed state material response is of importance. Here, soft sensitive clay is studied. Such clays behave contractant at failure, which for undrained conditions yields a strain-softening behavior governed by the generation of excess pore water pressure. Strain softening is further linked with material instability and the phenomenon of strain localization. In the case of shear band formation, internal pore pressure gradients are then expected to be present for globally undrained conditions in the sensitive clay due to its low permeability. In the present study, this hypothesis and its implications on the global response and shear band properties are investigated. Utilizing an experimental setup with a modified triaxial cell allowing for shear band formation, the effect of varying the displacement rate is studied. Onset of strain localization is interpreted to occur just before or at the peak shear strength. A strong rate dependency of the softening response is observed. Increasing displacement rates give raised brittleness in terms of the slope of the global softening curve due to accumulating pore pressure. Also, reduced shear band thickness and a shear band inclination approaching 45° are obtained for increasing rates. In the context of slope failure in such materials, the rate dependency in the post-peak state opens up for a large variation in behavior, all depending on time as an important factor.     

The Speeton Clay of North Yorkshire, England: An investigation of its geotechnical properties

The Lower Cretaceous in North Yorkshire is represented by a marine succession termed the Speeton Clay. The principal minerals in these clays consist of quartz, mica and clay minerals, although pyrite, glauconite and siderite are notable in some beds. The proportions of clay minerals change within the Speeton Clay and generally illite becomes more important in the younger horizons.

The Speeton Clay is a silty clay which generally has a high plasticity, the latter being influenced by the proportion of clay fraction present. All the clays tested were inactive although they would appear to have a potential for expansion. Neither the plasticity nor the moisture content were affected by depth.

In terms of its undrained shear strength, most of the Speeton Clay could be regarded as a stiff clay. It was also insensitive with generally a low slake-durability index. Its coefficients of volume compressibility were characteristic of overconsolidated and heavily overconsolidated days.     

Cut Slopes in Clayey Soils: Consolidation and Overall Stability by Finite Element Method

When a cut slope in a saturated clay is undertaken, a transient water flow occurs and stress transferences from the water to the soil skeleton take place in time (consolidation). Mainly in strongly overconsolidated clays, these stress transferences may determine swelling of soil and therefore reduction of its shear strength in time. However, the lowering of the water level associated to the cut increases effective mean stress, which may therefore counterbalance the above-mentioned effect. In the paper, the behaviour of a cut slope in an overconsolidated clay is analysed by a finite element program that incorporates the Biot consolidation theory (coupled analysis), with constitutive relations simulated by the p–q–θ critical state model. In addition, the variation in time of the overall stability is assessed with a computer program that uses the finite element results and formulations of the critical state soil mechanics. In order to achieve a more complex geotechnical interpretation of the problem, the analysis in time of the excess pore pressures, effective stresses, displacements and stress levels is also presented. Finally, comparisons of stability results are analysed by changing some parameters, namely the problem geometry (weight of excavated soil) and the over-consolidation ratio of the clay.     

粘土的双硬化模型

本文中建议了一个包含塑性体应变和偏应变两个硬化参数的封闭屈服面模型，可以同时应用于正常固结土和超固结土。该模型沿用了剑桥模型的参数测定方法，而大大拓展其适用范围。通过正常固结试样、超固结试样和应力路线转折条件下计算与试验曲线的对比，说明该模型有较广泛的适应性。