循环荷载下砂土的剪切硬化边界面本构模型

基于临界状态土力学框架,建立了一个适用于往返循环荷载作用的砂土边界面本构模型。采用无纯弹性域假设,认为受到反向荷载的瞬时土体就产生塑性变形,砂土的弹性区域退化为一个点。屈服面为倒子弹头型,由于砂土孔隙比与压力之间不存在惟一对应的关系,使得屈服面大小无法与体积应变直接耦合,故采用塑性偏应变而不是剑桥模型那种塑性体应变作为硬化参数。流动法则采用加入状态参数的修正的Rowe应力剪胀关系,体现了依赖状态的剪胀思想。屈服面大小的比值 反映了塑性模量的演化,并推导了 的表达式。只用1套参数,该模型就能合理地模拟砂土在不同密度和固结压力下循环荷载的应力-应变关系曲线。     

Development in modeling cyclic loading of sands based on kinematic hardening

In this paper, there is presented an elastoplastic constitutive model to predict sandy soils behavior under monotonic and cyclic loadings. This model is based on an existing model (Cambou‐Jafari‐Sidoroff) that takes into account deviatoric and isotropic mechanisms of plasticity. The flow rule used in the deviatoric mechanism is non‐associated and a mixed hardening law controls the evolution of the yield surface. In this research the critical state surface and history surface, which separates the virgin and cyclic states in stress space, are defined. Kinematic hardening modulus and stress–dilatancy law for monotonic and cyclic loadings are effectively modified. With taking hardening modulus as a function of deviatoric and volumetric plastic strain and with defining the history surface and stress reversal, the model has the ability to predict the sandy soils' behavior. All of the model parameters have clear physical meanings and can be determined from usual laboratory tests. In order to validate the model, the results of homogeneous tests on Hostun and Toyoura sands are used. The results of validation show a good capability of the proposed model. Copyright © 2009 John Wiley & Sons, Ltd.     

Rotational kinematic hardening model for three‐dimensional stress reversals in sand

A kinematic hardening mechanism has previously been proposed to capture the behavior of soil during large stress reversals in the triaxial plane. This mechanism is now extended to the principal stress space. It incorporates rotation and intersection of yield surfaces to achieve a consistent and physically rational fit with experimentally observed soil behavior during large three‐dimensional stress reversals. An existing elasto‐plastic model with isotropic hardening is used as the basic framework to which the rotational kinematic hardening mechanism has been added. The new combined model preserves the behavior of the isotropic hardening model under monotonic loading conditions, and the extension from isotropic to rotational kinematic hardening under three‐dimensional conditions is accomplished without introducing new material parameters. The framework of the model is presented here with some comparisons between theoretical and experimental directions of strain increment vectors to indicate the potential of the model. Copyright © 2008 John Wiley & Sons, Ltd.     

基于旋转运动硬化理论的修正剑桥模型的改进

修正剑桥模型是最早建立和得到广泛承认的经典土体弹塑性模型之一,但不能模拟应力路径转折时土体的应力-应变特性以及应力引起的各向异性。将旋转运动硬化理论引入到剑桥模型中,给出了椭圆屈服面的旋转运动硬化机制,在不增加任何模型参数的情况下,把等向硬化的修正剑桥模型扩展为旋转运动硬化模型。扩展的新模型既保留了单调加载时的等向硬化的特性,又能反映应力路径转折时土体的本构特性与应力诱发的各向异性,初步验证了模型的有效性。     

Implicit numerical integration for a kinematic hardening soil plasticity model

Soil models based on kinematic hardening together with elements of bounding surface plasticity, provide a means of introducing some memory of recent history and stiffness variation in the predicted response of soils. Such models provide an improvement on simple elasto‐plastic models in describing soil behaviour under non‐monotonic loading. Routine use of such models requires robust numerical integration schemes. Explicit integration of highly non‐linear models requires extremely small steps in order to guarantee convergence. Here, a fully implicit scheme is presented for a simple kinematic hardening extension of the Cam clay soil model. The algorithm is based on the operator split methodology and the implicit Euler backward integration scheme is proposed to integrate the rate form of the constitutive relations. This algorithm maintains a quadratic rate of asymptotic convergence when used with a Newton–Raphson iterative procedure. Various strain‐driven axisymmetric triaxial paths are simulated in order to demonstrate the efficiency and good performance of the proposed algorithm. Copyright © 2001 John Wiley & Sons, Ltd.     

混凝土混合硬化黏塑性损伤统本构模型

传统的黏塑性统一本构模型只能用于金属类材料的本构分析,而不适用于混凝土的本构分析.基于此,以不可逆热力学理论为基础,选定了涉及混合硬化参量的Helmholtz自由能函数,推导出混凝土的黏塑性损伤耦合本构模型.通过对混凝土硬化模型的分析,构造了含运动硬化和等向硬化内变量的非相关流动势函数,又推导出流动方程以及内变量演化方程.所建模型对混凝土试验曲线的数值模拟显示,其能够正确描述混凝土的率相关性、非弹性体积膨胀特性、加载过程的软硬化特性以及由断裂和损伤引起的应力软化现象.     

A two-scale constitutive framework for modelling localised deformation in saturated dilative hardening material

Undrained deformation of dilative sand generates negative excess pore pressure. It enhances the strength, which is called dilative hardening. This increased suction is not permanent. The heterogeneity at the grain scale triggers localisations causing local volume changes. The negative hydraulic gradient drives fluid into dilating shear zones. It loosens the soil and diminishes the shear strength. It is essential to understand the mechanism behind this internal drainage and to capture it numerically. The purpose of this paper is to develop a macroscopic constitutive relationship for the undrained deformation of saturated dense sand in the presence of a locally fully or partially drained shear band. Separate constitutive relations are generated for the band and intact material. Both time and scale dependence during pore fluid diffusion in saturated sand are captured, eliminating the mesh dependency for finite element implementations. The model is applied to the Gauss points that satisfy the bifurcation criterion. The proposed method is calibrated to recreate the undrained macroscopic response bestowed by an extra-small mesh. The microscopic behaviours inside and outside shear band predicted by this model are qualitatively in good agreement with individual material point behaviours inside and outside the shear band in the extra-small mesh. Depending on the loading rate and the shear band thickness, the response inside the band can be fully or partially drained, which governs the ultimate global strength. The calibrated model is exploited to simulate an upscaled biaxial compression test with semipermeable boundaries.     

基于临界状态模型的砂土非共轴本构模拟

传统的砂土本构理论隐含了应力和塑性应变率的共轴条件,无法客观描述主应力轴旋转试验中的非共轴现象,并且当密度和围压变化较大时也不适用。基于材料状态相关砂土临界状态概念,将Pietruszczak和Stolle所提出的砂土本构模型进行了改进,并在模型中引入非共轴塑性流动理论来描述非共轴现象。通过对单剪试验和空心圆柱试验进行数值模拟,表明基于临界状态理论的非共轴模型能够合理描述主应力轴旋转过程中砂土的非共轴变形特性     

剪胀性砂土本构模型的研究

对于中密砂和密砂,剪胀现象是一个比较显著的特征,相变线是描述剪胀现象的一条特征线,它比较容易确定。因此,以相变线作为状态参考线,提出并定义了一种状态参量,在已有弹塑性模型基础上,将状态参量引入剪胀方程和塑性硬化模量表达式中,建立了一个新的考虑砂土剪胀性的本构模型。该模型尤其适合于中密砂和密砂,对于松砂则退化为原弹塑性模型。模型形式较简单,包含9个材料参数,且比较容易确定。数值模拟结果与试验数据吻合较好。     

胶结结构性土统一硬化模型

结构性土颗粒间的胶结使试样剪切破坏最终应力比高于相应重塑土,也限制了试样剪切时体积应变的自由发挥。在考虑结构垮塌为主的结构性土统一硬化(UH)模型基础上,将应力空间中静止的临界状态线扩展为动态的移动临界状态线。据此,通过建立新的屈服面方程并修正剪胀方程,将结构性土统一硬化(UH)模型扩展为胶结结构性土统一硬化(UH)模型。相对于原模型,新模型增加了1个模型参数,即初始胶结应力,反映土颗粒之间的初始胶结作用。通过4种结构性土试验数据与模型预测对照表明：所提模型能够较合理地描述结构性土等向压缩、常规三轴排水与不排水剪切等特性。     

A kinematic hardening based model for unsaturated soils considering different hydraulic conditions

At present, several of the existing elastoplastic constitutive models are adapted for describing the stress–strain behavior of unsaturated soils. However, most of them present certain limitations in this field. These limitations can be related to the basic model and/or added unsaturated state variables and formulations. In this regard, inability to model the hydro‐mechanical behavior in constant water (CW) conditions is an example of these limitations. In this paper, an advanced version of CJS model is selected for adaptation to the unsaturated states. Adaptation to unsaturated states is achieved in the framework of effective stress approach. Effective stress equation and unsaturated state variables are selected based on the recent research existing in the literature. The developed model is capable of describing the complex behavior of unsaturated soil in the CW condition in addition to predicting the behavior at failure and post–failure, nonlinear elastoplastic behavior at low levels of stress and strain (by selecting a very small elastic domain), as well as wetting and collapse behaviors. In order to validate the model, results of triaxial tests in CD and CW conditions are used. The validation results indicate the good capability of the proposed model. Behavior of the unsaturated soils during wetting is an important issue. For this reason, the model is also evaluated based on the results of wetting and collapse triaxial tests. A comparison between the tests and simulation results shows that the model is able to predict the soil behavior under the wetting path. Copyright © 2016 John Wiley & Sons, Ltd.     

无状态变量的状态依赖剪胀方程及其本构模型

粗粒土的剪胀行为具有状态依赖特性。为了考虑这一特性,不同的状态依赖变量被唯像地提出,并被经验性地内嵌入已有剑桥、修正剑桥等剪胀方程中。基于分数阶梯度律,用理论推导出了分数阶状态依赖剪胀方程,并阐述了分数阶数的物理意义。所得剪胀比大小受3个因素影响：分数阶求导阶数、当前加载应力以及当前应力到临界状态应力的距离。当分数阶求导阶数从1开始增大时,分数阶剪胀曲线自修正剑桥剪胀曲线向剑桥剪胀曲线移动;而当求导阶数从1开始减小时,分数阶剪胀曲线逐渐远离修正剑桥剪胀曲线;当求导阶数等于1时,分数阶剪胀曲线与修正剑桥剪胀曲线重合。为验证所提出的状态依赖剪胀方程,基于该方程进一步建立了砂土的状态依赖分数阶塑性力学本构模型,并对砂土和堆石料的三轴排水与不排水试验结果进行了模拟。研究表明,基于状态依赖分数阶剪胀方程建立的本构模型,可以合理地描述砂土在不同初始状态及加载条件下的应力-应变行为。与砂土UH模型预测结果对比发现,UH模型预测较好。     

砂土的物态本构模型

基于砂土的压缩回胀性、剪切非线性及剪缩剪胀性的系统分析和包括循环荷载、主应力轴旋转及应力路径偏转等复杂应力条件下的复杂变形反应,得到了三类应力-应变基本关系。在剪缩剪胀应力-应变关系中,引入了由偏应变分量确定的应变路径长度变量,揭示了应力主轴旋转、应力路径偏转引起的剪缩剪胀性。将这些基本关系与循环荷载下砂土的物态变化相联系,建立了砂土的物态动本构关系。     

A constitutive model for unsaturated soil interfaces

A constitutive model based on the disturbed state concept is presented to describe the behavior of interfaces in unsaturated soil. The model is an extension of an existing model developed for a sand–steel interface. As opposed to the original model, the modified model incorporates two independent stress variables, which are the net normal stress and matric suction. The saturated and dry state of the interface can be modeled as a special case using the constitutive model presented in this paper. The modified model is capable of capturing the main features of unsaturated interfaces observed during laboratory testing, including increasing shear strength and strain softening with increasing suction and net normal stress and increasing dilatancy with increasing suction. Laboratory tests were carried out on unsaturated interfaces in a modified direct shear test apparatus. The observed behavior of interfaces between unsaturated soil and steel plates (rough and smooth) is presented in comparison with model predictions. Copyright © 2008 John Wiley & Sons, Ltd.     

A phenomenological constitutive model for rocks with shear failure mode

In this paper, a phenomenological constitutive model is proposed to simulate the stress–strain behaviours of intact rocks with shear failure mode. The model captures a wide range of behaviours of rock material such as elastic, plastic, strain softening, ‘Class II’, strain localization, elastic modulus degradation, etc. The sensitivity of the stress–strain relation on the parameters is also investigated. Typical results obtained by testing a number of granite and marble specimens are used to validate the proposed model. Copyright © 2001 John Wiley & Sons, Ltd.     

考虑颗粒破碎效应的堆石料静动力本构模型

为合理反映颗粒破碎对堆石料力学特性的影响,基于试验结果分析,得出了堆石料在压缩和剪切作用下的颗粒破碎特性规律。通过引入压缩破碎和剪切破碎的相关参数,借鉴已有本构模型的合理定义,吸收临界状态理论和边界面理论的优点,发展了考虑颗粒破碎和状态相关的堆石料静动力统一弹塑性本构模型,并阐述了模型参数的确定方法。该模型不仅能够反映堆石料在静力荷载作用下的低压剪胀、高压剪缩、应变软化和硬化等特性,还能够反映在循环荷载作用下应力-应变的滞回特性和残余变形的累积效应。最后为验证模型的合理性,分别对堆石料的静力三轴和循环三轴试验进行了数值模拟预测,结果表明：模型预测与试验数据吻合良好,所提出的本构模型能够合理地描述颗粒破碎对堆石料静动力变形特性的影响。     

考虑温度影响的冻土蠕变本构模型

冰作为冻土的基本组成部分,对冻土蠕变的加速蠕变阶段有重要影响。温度通过影响冻土中冰的冻结与融化过程,及其黏塑性流动,引起冻土结构的强化与弱化,从而成为了决定冻土蠕变力学行为的关键因素之一。同时,外部应力也造成冻土的强化与弱化,影响着冻土的蠕变。通过引入硬化因子与损伤因子来考虑温度、应力造成的冻土材料强化与弱化。硬化因子H代表了蠕变过程中强化效应的大小,而损伤因子D则代表了由弱化效应造成的冻土材料相关参数的折减比例,进而提出了适用于冻土的改进西原蠕变本构模型。该模型预测值与试验数据的比较表明：改进的模型不仅能较好地描述初始蠕变阶段、稳定蠕变阶段,而且相比于传统模型,能更好地描述冻土加速蠕变阶段,具有合理性与一定的实用性。     

考虑冻融循环的季冻区高等级公路路基填土本构模型研究

在季冻区高等级公路路基场地取土,制成重塑土试件进行冻融循环试验。而后对经历不同冻融循环次数的土体进行三轴固结排水（CD）剪切试验,测得应力、应变数据。在大量试验数据的基础上,遵循岩土本构变化规律,以岩土弹塑性力学为理论基础,在p-q空间,以椭圆方程拟合体变屈服面,以抛物线方程拟合剪切屈服面,采用塑性关联流动,得到考虑冻融循环条件下融土的弹塑性本构方程。通过编程计算后,对比计算与试验结果,证明该双屈服面本构模型能够较为准确地预估冻融循环土体的应力-应变关系,该本构模型对季冻区路基土体长期稳定性分析和工程预测具有重要参考意义。