Evaluation of an anisotropic elastoplastic–viscoplastic bounding surface model for clays

An anisotropic time-dependent bounding surface model for clays is developed by generalizing a previous time-independent model that adopts a flexible bounding surface. It is based on the framework for coupled elastoplasticity–viscoplasticity for clays and Perzyna’s overstress theory. Three viscoplastic parameters were introduced and explained in detail. The model was validated against undrained creep tests for both isotropically and anisotropically consolidated clays, undrained and drained stress relaxation tests on some undisturbed clays, and undrained triaxial tests with varying strain rates on natural Hong Kong marine deposit clay. The general agreement between the model simulations and test results was satisfactory. The varying effects of lower-level parameters were discussed on the undrained multistage stress relaxation response for normally consolidated soils which had been ignored in literature. The flexibility of the model in capturing the shear strengths, which is the unique feature of the current model, was shown in the simulations of time-dependent triaxial tests on Taipei silty clay. All the simulations show that the proposed model is a relatively practical model considering both anisotropy and time dependency of clays.     

各向异性弹黏塑性模型在ABAQUS中的研发

通过将Perzyna过应力理论与临界状态理论相结合,并引入Wheeler旋转硬化法则,提出一个能描述土体初始各向异性及应力诱发各向异性的三维弹黏塑性本构模型。模型考虑流变发生的下限,在三维应力空间,模型存在形状相似的静屈服面及动态加载面。采用缩放形式的幂函数。本构模型数值算法采用回映算法,借助ABAQUS软件UMAT子程序接口实现。并通过对三轴不排水蠕变试验的模拟,确定合适的积分步长。此后,分别对三轴不排水蠕变试验及常应变率三轴不排水剪切试验进行了模拟。模拟中通过设置不同参数值,可将模型退化为各向同性模型,并对这两种模拟结果进行了比较。模拟结果表明：(1) 对于三轴不排水蠕变,在低剪应力水平下,各向同性模型和各向异性模型模拟的结果相差不大,而在高剪应力水平下,各向异性模型模拟结果更接近试验结果;(2) 对于常应变率加载试验的模拟,模型合理反映了土体不排水强度随着加载速率的增大而增大现象。     

Delayed plastic model for time‐dependent behaviour of materials

A delayed plastic model, based on the theory of plasticity, is proposed to represent the time‐dependent behaviour of materials. It is assumed in this model that the stress can lie outside the yield surface and the conjugate stress called static stress is defined on the yield surface. The stress–strain relation is calculated based on the plastic theory embedding the static stress. Thus, the stress–strain relation of the model practically corresponds to that of the inviscid elastoplastic model under fairly low rate deformation. The delayed plastic model is coupled with the Cam‐clay model for normally consolidated clays. The performance of the model is then examined by comparing the model predictions with reported time‐dependent behaviour of clays under undrained triaxial conditions. It is shown that the model is capable of predicting the effect of strain rate during undrained shear and the undrained creep behaviour including creep rupture. Copyright © 2001 John Wiley & Sons, Ltd.     

原状和重塑海洋粘土经历动载后的静强度衰减

根据海洋粉质粘土原状和重塑土样的动三轴试验结果, 比较和分析了两种土样应力、应变、孔隙水压力和不排水抗剪强度行为, 得到了土样的静不排水抗剪强度衰减与波浪荷载作用下土样产生的动应变以及平均累积孔压之间的关系。 并将波浪荷载作用使土样内孔压升高, 有效应力降低, 形成拟超固结现象的理论, 应用到土样不排水抗剪强度衰减与平均累积孔压之间关系的分析;同时与重塑土样的超固结静态剪切试验结果进行比较, 得到了土样在波浪荷载作用后的归一化不排水抗剪强度与拟超固结比之间的关系式。 建议以少量原状土样, 配合大量重塑土样的动三轴试验结果, 实现对实际海洋粘土地基在波浪荷载作用后的静不排水抗剪强度衰化规律的评估。     

Undrained Shear Strength of Normally Consolidated and Overconsolidated Clays from Pressuremeter Tests: A Case Study

Undrained shear strength (s_u) of foundation soil of Marquette interchange near Milwaukee, Wisconsin was evaluated from the results of a number of pressuremeter tests conducted on normally consolidated (NC) organic silts and overconsolidated (OC) silty clay. The s_u-values were interpreted from traditional closed-form methods. The pressuremeter geometry and test sequence as well as response of the soil profiles were also simulated using axisymmetric finite element (FE) method with Cam-Clay soil model. The Cam-Clay model parameters were estimated from laboratory tests on undisturbed soil samples. Results show that the s_u estimated from the rate of cavity pressure change with volumetric strain (referred to as direct traditional method) is almost twice the s_u estimated from an indirect traditional method that estimates s_u from shear modulus, in situ horizontal stress, and ultimate cavity pressure obtained from the cavity pressure curves. The s_u-values predicted from the FE models are lower than those estimated from the traditional methods and shows that the assumption of infinite pressuremeter length in traditional methods results in overprediction of undrained shear strength by a factor of 1.5 for NC clay and 2.2 for OC clay. The results of finite element analysis considering Cam-Clay soil model and finite length for pressuremeters suggest the undrained shear strength of 63 ± 7 kPa for NC organic silt and 259 ± 68 kPa for OC silty clay.     

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

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

粉质黏土的真三轴试验及强度特性研究

为研究复杂应力路径下粉质黏土强度特性,利用TSW-40型土壤真三轴仪,开展了不同试验条件下粉质黏土的真三轴试验。研究了真三维应力状态下粉质黏土应力-应变曲线变化规律,探讨了含水率、固结围压和中主应力比等试验条件对粉质黏土强度的影响规律,对比分析了真三维应力状态下 平面上粉质黏土试验值与Mohr-Coulomb（简称M-C）强度准则和双剪强度准则的差异。研究结果表明：真三维应力状态下粉质黏土的广义剪应力-剪应变曲线主要表现为应变硬化型。粉质黏土强度随中主应力比和固结围压的增大而增大,随着含水率的增大而减小。 平面上粉质黏土的强度包络线随静水压力 的增大向外逐渐扩展,随含水率 的增大向内逐渐收缩。 平面上M-C强度极限线和双剪强度屈服极限线与试验值对比分析表明,M-C强度准则描述土体强度时结果偏小,采用双剪强度准则时结果偏大。     

NGI‐ADP: Anisotropic shear strength model for clay

Many geotechnical problems involve undrained behavior of clay and the capacity in undrained loading. Most constitutive models used today are effective stress based and only indirectly obtain values for the undrained shear strength. To match the design profiles of undrained shear strengths, in active (A), direct simple shear (D) and passive (P) modes of loading are complicated. This paper presents the elastoplastic constitutive model NGI‐ADP which is based on the undrained shear strength approach with direct input of shear strengths. Consequently, exact match with design undrained shear strengths profiles is obtained and the well‐known anisotropy of undrained shear strength and stiffness is accounted for in the constitutive model. A non‐linear stress path‐dependent hardening relationship is used, defined from direct input of failure strains in the three directions of shearing represented by triaxial compression, direct simple shear and triaxial extension. With its clear input parameters the model has significant advantages for design analysis of undrained problems. The constitutive model is implemented, into finite element codes, with an implicit integration scheme. Its performance is demonstrated by a finite element analysis of a bearing capacity problem. Copyright © 2011 John Wiley & Sons, Ltd.     

软黏土蠕变特性试验研究：回顾与发展

大量的室内和现场试验表明软黏土存在长期变形难以收敛的问题,为了更深入地认识软黏土的蠕变特性,首先,从微观物理化学角度探讨蠕变机制;然后,深入分析一维应力条件下次固结系数的演变规律以及确定方法和三轴排水、不排水蠕变特性及黏土的长期抗剪强度、室内旁压试验和现场试验等复杂应力下的黏土蠕变特性等,并且讨论黏土蠕变特性与应变速率和应变松弛特性在一维和三轴条件下的关联性;最后,从黏土次固结特性与微观结构的相关性、如何准确描述次固结系数的非线性及蠕变过程中应力剪胀剪缩关系发展等三方面更深入地探讨软黏土蠕变特性。     

饱和软黏土动力学特性循环扭剪试验研究

以饱和软黏土的循环扭剪试验研究为基础,阐明不固结不排水条件下饱和软黏土的动力学特性。在循环扭剪和循环三轴两种不同试验应力状态下,通过研究不同围压和不同静、动应力组合下饱和软黏土的应力等效破坏关系和应变等效破坏关系,提出在循环荷载作用下饱和软黏土的循环破坏同样遵循Mises屈服准则,且在Mises屈服准则下,饱和软黏土的循环强度趋于不变量。循环破坏的过程可以等效为一种拟静力弹塑性循环蠕变,建立了饱和软黏土循环累积变形随静荷载、循环荷载以及循环破坏振次之间的关系式。上述结论分别从应力和变形两个方面阐述饱和软黏土的动力学特性,与试验应力状态和围压无关,可以推广到一般应力状态下。     

黄河三角洲粉土循环剪切特性的试验研究

利用土工静力-动力液压三轴-扭转多功能剪切仪,对黄河三角洲埕北海域原状粉土试样进行了循环三轴和循环扭剪试验,以此探讨了原状海底粉土的动力变形与强度特性,确定了最大动剪切模量及动力变形参数随剪应变的变化规律。基于循环三轴试验结果,建立了残余孔隙水压力发展变化规律的双曲线经验模式,对动强度关于循环次数的依赖性提出了经验关系,确定了有关试验参数及其关于固结应力比的依赖性。与循环三轴试验结果相对比,认为在土工抗震设计中采用循环三轴试验结果稍偏于保守。     

循环荷载对粉质黏土力学性质的影响

通过等压固结静、动三轴试验,研究了渤海湾粉质黏土在循环荷载作用下的动力性质和循环荷载后不排水静力性质。试验结果表明,循环应力幅值比越大,平均轴向应变和轴向应变幅值越大;循环应力幅值比达到0.4时,平均轴向应变和轴向应变幅值随着循环周数增加迅速增大;循环应力幅值比相同,固结应力越大,轴向应变幅值越大,而平均轴向应变越小。在较大的循环应力幅值比下,平均孔压比值和孔压幅值比值随着循环周数的增大会达到稳定;循环应力幅值比越大平均孔压比值和孔压幅值比值均越大;相同循环应力幅值,固结应力越大平均孔压比值越小,而固结应力对孔压幅值比值影响较小。循环荷载的作用会导致循环荷载后不排水剪在q-p’平面上有效应力路径和孔压发展表现出超固结土的性质。     

连云港软黏土的不排水强度试验研究

对连云港软黏土不扰动样和重塑样进行三轴CU对比试验研究。结果表明,不扰动样的不排水强度包线是由以屈服压力为界的双折线组成,而重塑样的不排水强度包线可以用通过圆点的直线表示。屈服前阶段,不扰动样由于受结构性影响产生的抵抗作用,其不排水强度包线高于重塑样;屈服后阶段,不扰动样的不排水强度包线介于具有不同初始含水率的重塑样的包线之间。通过对比分析,发现不扰动样与重塑样的不排水强度包线存在差异的主要原因是由于试样的不同初始含水率引起的。不扰动样在屈服后阶段的含水率变化规律与不排水强度增长规律与重塑土相一致,采用Burland提出的孔隙指数Iv,可将不同初始含水率重塑样及屈服后阶段不扰动样的压缩性状较好地归一化。     

上海软土动剪切模量衰减规律的试验研究

为了获得上海软土在循环动力荷载作用下的动力特性,特别是动剪切模量的衰减规律,采用原状第④层上海软土进行等向固结不排水动三轴试验,分析了加载频率和动应力比 对动剪切模量的影响。试验结果表明,对于上海第④层软土,动应力比达到0.2时,土样才发生破坏;土样的剪切模量-累积轴应变曲线不受加载频率和动应力比的影响;剪切模量的衰减下降主要发生在轴应变0～1%区间。最后提出了一个描述剪切模量衰减规律的数学模型。     

Lower bound limit analysis of an anisotropic undrained strength criterion using second‐order cone programming

In this paper, the formulation of the lower bound limit analysis of an anisotropic undrained strength criterion using second‐order cone programming is described. The finite element concept was used to discretize the soil mass into 3‐noded triangular elements. The stress field was modeled using a linear interpolation within the elements while stress discontinuities were permitted to occur at the shared edges of adjacent elements. An elliptical yield criterion was adopted to model the anisotropic undrained strength of the clay. A statically admissible stress field was defined by enforcing the equilibrium equations within all triangular elements and along all shared edges of adjacent elements, stress boundary conditions, and no stress violation of the anisotropic strength envelope cast in the form of a conic quadratic constraint. The lower bound solution of the proposed formulation was solved by second‐order cone programming. The proposed formulation of the anisotropic undrained strength criterion was validated through comparison of the model's predictions with the known exact solutions of strip footings, and was applied to solve undrained stability of a shallow unlined square tunnel. Computational performance between the proposed approach of second‐order cone programming and linear programming was examined and discussed.     

Stress-strain behaviour of reconstituted illitic clay at different temperatures

Tests on specimens of reconstituted illitic clay have examined the influence of temperature on the mechanical behaviour of clay soils. The program involved consolidation to effective confining pressures up to 1.5 MPa, heating to 100°C, and tests on normally consolidated and overconsolidated specimens with OCR = 2. The tests included isotropic consolidation, undrained triaxial compression with pore water pressure measurement, drained tests along controlled stress paths to investigate yielding behaviour, and undrained tests which involved heating and measurement of the resulting induced pore water pressures. The large strain strength envelope is independent of temperature. However, peak undrained strengths increase with temperature because smaller pore water pressures are generated during shearing. An important contribution from the study is a series of results for the yielding of illitic clay at three different temperatures. For the first time, there is clear evidence of yield loci decreasing in size with increasing temperature. An associated flow rule can be assumed without serious error. The results contribute to the confirmation of a thermal elastic-plastic soil model developed by the authors from cam clay following the addition of a small number of extra assumptions. Depending on the initial stress state, heating under undrained conditions may produce shear failure.     

德州松软土蠕变特性及其模型研究

利用TSS10型三轴流变试验机对德州松软土开展不同围压下的三轴固结不排水蠕变试验,获得变形、应变率与应力和时间的变化关系,分析蠕变性状和特征,建立其分数阶类Maxwell蠕变模型。结果表明:德州松软土具有明显非线性蠕变特性,高应力水平的蠕变变形大于低应力下的蠕变变形;介绍分数阶蠕变模型原理,引用了分数阶类Maxwell蠕变模型,并求取德州松软土的类Maxwell蠕变模型参数和验证其描述德州松软土蠕变变形的适宜性。     

超固结黏土单调和耦合循环的剪切特性研究

针对超固结黏土空心试样,利用土工静力-动力液压三轴-扭转多功能剪切仪,在均等固结条件下进行了单调扭剪和三轴-扭转耦合循环剪切试验。试验结果表明：不同超固结比黏土的单调扭剪强度也可由正常固结黏土的单调扭剪强度得到,得到了不同超固结比下饱和黏土的强度及模量的退化规律;随着超固结比的增大,相同破坏循环次数的动应力比和临界循环应力比均线性增大;超固结比对耦合循环剪切的孔隙水压力的发展模式影响显著。参考Yasuhara的测量方法,采用荷载停止后继续采集孔压的方法可以更好地反映黏土在耦合循环荷载下产生的真实孔压和孔压的增长情况。提出的综合应变式同时考虑了剪切变化和正向偏差变形的共同效应,适合作为主应力连续旋转的耦合循环剪切试验的破坏标准。     

A viscoplastic subloading soil model for rate‐dependent cyclic anisotropic structured behaviour

This paper presents a new purely viscoplastic soil model based on the subloading surface concept with a mobile centre of homothety, enabling the occurrence of viscoplastic strains inside the yield surface and avoiding the abrupt change in stiffness of the traditional overstress viscoplastic models. This is required for overconsolidated soils. The model is formulated to reproduce the soil rate‐dependent behaviour under cyclic loading (changes in loading direction) and incorporates both initial and induced anisotropy, as well as destructuring. The model shows good qualitative response to some imposed three‐dimensional stress paths under quasi‐inviscid (elastoplastic) behaviour. Some of the main time‐dependent aspects of soil behaviour that the model is capable of reproducing were also illustrated. The capability of the model to adequately reproduce the results from an undrained triaxial test performed on stiff overconsolidated clays from the Lisbon region (Formação de Benfica), with an unloading–reloading deviatoric stress cycle at constant mean stress, that incorporates a series of staggered fast loading and creep stages, was evaluated. The model was able to reproduce well the main observed aspects of the time‐dependent stress–strain response and pore pressure evolution of a stiff overconsolidated clay under complex loading. The revised and generalised viscoplastic subloading surface concept is viable and can be applied to a consistent extension to viscoplasticity, including in the interior of the yield surface, of existing elastoplastic models formulated for soils and other materials. Copyright © 2016 John Wiley & Sons, Ltd.