循环荷载作用下超固结软黏土软化-孔压模型研究

随着循环次数的增加,循环荷载作用下正常固结软黏土孔隙水压力逐渐增加,而孔压的增加将导致土体发生软化现象。对于超固结软黏土,循环荷载作用下在循环初期将产生负孔压,而此时也发生循环软化现象,这显然与有效应力原理是相矛盾的。以往的研究往往通过建立软化指数与循环次数的关系来描述土体的循环软化特性,从而不能反映土体的循环过程中残余孔压变化对土体循环软化的影响。通过对杭州饱和软黏土进行应力控制的循环三轴试验,对循环荷载作用下超固结软黏土的软化特性及孔压发展规律及两者关系进行了研究。在试验的基础上建立了超固结软黏土循环软化-孔压模型,该模型反映了残余孔压增长对超固结软黏土循环软化特性的影响规律。     

循环荷载作用下饱和软黏土的损伤模型

基于各向同性弹塑性损伤和Prevost模型的基本理论,把弹塑性等向硬化、运动硬化和各向同性损伤结合起来,推导了循环荷载作用下不排水饱和软黏土的弹塑性动力损伤本构模型。由循环累积塑性偏应变建立损伤演化方程,以描述循环加载对软黏土结构的破坏作用,并通过对循环三轴试验的模拟验证了模型的有效性。结果表明,提出的模型能较好地描述饱和软黏土在循环荷载作用下的变形、孔压变化及模量损伤过程。     

一种适用于饱和黏土循环动力分析边界面塑性模型的隐式积分算法

修正了传统隐式回映算法,建立了适用于饱和黏土循环动力分析的边界面塑性模型的完全隐式积分格式。该模型基于无弹性域概念和临界状态理论,采用各向同性、运动硬化准则、旋转的边界面,并引入表征土体结构损伤和重塑程度的损伤变量以反映循环载荷作用下饱和黏土的各向异性、刚度、强度软化及塑性变形累积等特征。针对等压固结 和偏压固结 的饱和高岭黏土的不排水三轴试验进行模拟,采用不同的应变增量步长进行计算,并与试验数据对比,结果表明,修正隐式回映算法应用于该类边界面模型的合理性、积分格式的精确性和稳定性;另外,结合有限元软件自动时间步长的增量迭代解法,对饱和黏土应力控制的不排水动三轴试验进行预测,结果表明,修正的适用于该边界面的塑性模型隐式回映算法可以得到比较合理的数值分析结果,能够反映饱和黏土的循环刚度的退化和强度的弱化等动力特性。     

循环荷载下非饱和结构性黏土的弹塑性双面模型

非饱和黏土的结构性能够显著影响其力学特性。基于非饱和土经典模型BBM（Barcelona basic model）和一种可描述循环塑性的硬化法则,引入体积破损率的作为标准土体结构破损的参数,建立了一个描述常吸力下非饱和结构性黏土静态及动态力学特性的弹塑性双面模型。模型在应力空间中包含与重塑非饱和土屈服面几何相似的结构性边界面和加载面,采用径向映射法则和可移动的记忆中心原理,通过结构性边界面和加载面在应力空间中的演化来反映循环加载过程中非饱和结构性黏土的循环塑性特征和结构损伤过程。通过与相关非饱和黏土控制吸力试验数据的比较,表明该模型能够较好地反映静态加载下非饱和结构性黏土的力学特性,而模型预测的循环荷载下的应力?应变特征也具有一定的合理性。     

考虑循环软化特性的饱和软土弹塑性本构关系研究

将软化指数关系与非等向硬化模量场理论相结合,研究了可描述循环荷载作用下饱和软土软化特性的增量弹塑性模型。该模型借助硬化模量的插值和映射中心的移动,在偏应力空间中构造硬化模量的演化规则;通过在弹塑性模量插值函数中引入初始弹塑性模量软化系数,模拟循环荷载作用下软土的刚度软化特性;通过引入硬化模量调整系数,增强循环加载时应力-应变曲线的滞回特性;再通过引入反映应变累积速率和大小的模型参数,描述循环加载时软土的应变累积特性。利用Idriss提出的指数关系式近似拟合软化系数随应力循环次数的变化关系,并通过引入循环应力参数建立了循环软化系数与静应力水平和循环应力水平的关系。阐述了确定模型参数的方法,并利用模型预测了相关试验结果,通过预测结果与试验结果的对比,验证了该模型描述循环荷载作用下软土软化特性的可行性。     

部分排水时饱和粉质黏土变围压循环三轴试验研究

围压循环变化将引起土体孔压累积,加剧土体累积塑性应变,导致地基灾变事故。为研究围压循环变化对饱和粉质黏土变形特性的影响,利用GDS双向振动三轴仪,进行部分排水时饱和粉质黏土的变围压循环三轴试验,研究不同应力路径斜率和循环动应力比时饱和粉质黏土的孔压、轴向累积塑性应变和体应变发展规律,建立部分排水条件时考虑循环围压和循环动应力耦合作用的地基粉质黏土累积应变数学模型。试验结果表明：部分排水时饱和粉质黏土孔压、轴向累积塑性应变和体应变均随循环动应力比和应力路径斜率的增大而增大;动孔压比随着振动次数的增加明显分为急剧增加、快速下降和持续平稳3个阶段;轴向累积塑性应变和体应变均随振动次数的增加而增大。当加载次数超过2 500次时,粉质黏土孔压趋于平稳,变形速率略有降低,但变形持续增加。补充试验结果表明,所建立模型与试验结果具有较高的一致性。研究结果将为交通循环荷载导致的地基灾变控制技术提供理论依据。     

考虑饱和黏土初始各向异性的循环边界面塑性模型

天然土体的初始各向异性通常可对其后继循环特性产生显著影响。现有考虑循环载荷作用的土体弹塑性模型,往往采用类似修正剑桥模型的椭圆形屈服面,已有研究表明,该椭圆形屈服面因其拉伸弹性区域偏大,针对天然K0固结状态的土体,其计算精度较差。基于新近提出的广义各向同性硬化准则,在边界面方程中引入初始各向异性张量,并采用空间滑动面破坏准则（SMP）的变换应力法,建立了能考虑饱和黏土初始各向异性的循环边界面塑性模型。分别针对等压和偏压固结的饱和黏土静、动三轴试验进行模拟,结果表明,该模型能合理反映土体的初始各向异性及其后继循环动力特性。     

恒定主应力轴偏转循环荷载下饱和软黏土累积孔压特性分析

现有的饱和软黏土循环荷载累积孔压显式模型未能反映路基土体单元在土体自重作用下主应力轴的偏转现象。在确定饱和软黏土固结时间的基础上,在固结压力分别为100、150、200 kPa下,对软黏土进行主应力轴偏转角为15°、30°、45°的不排水循环加载试验,得到了固结围压、主应力轴旋转角及动应力比对循环累积孔压影响的规律。基于第一次循环加载累积孔压变化规律和修正动偏应力水平,建立了能考虑恒定主应力轴偏转角循环加载下累积孔压显式模型,并利用恒定主应力轴偏转角循环加载试验结果验证了模型的合理性。建立的循环荷载下饱和软黏土累积孔压显式模型参数易于确定,参数具有明确的物理意义且能考虑主应力轴的恒定偏转现象。     

主应力轴旋转下K0固结饱和粉质黏土孔压及变形特性试验研究

主应力旋转将加速孔压和塑性应变的累积,影响土体的力学特性。为研究主应力轴循环旋转对 固结饱和粉质黏土孔压和变形的影响,利用GDS-HCA试验系统开展了不同应力路径的 固结饱和粉质黏土不排水循环三轴试验和循环扭剪试验,研究了主应力轴循环旋转对 固结饱和粉质黏土累积塑性应变、孔压等动力特性的影响。结果表明：动剪应力幅值、心形应力路径所包围的面积、循环偏应力幅值等均随循环剪应力比和循环动应力比的增大而增大。试样未发生破坏时, 固结饱和粉质黏土的孔压比随着振动荷载作用次数的增加而增大;当试样发生破坏时,孔隙压力迅速消散,孔压比急剧下降。循环扭剪试验时试件产生的轴向累积塑性应变始终大于循环三轴试验产生的累积塑性应变,说明主应力轴循环旋转会加速试件轴向应变的累积。试件的累积塑性应变随循环剪应力比和循环动应力比的增大而增大。循环动应力比越大时,主应力轴旋转造成的轴向累积塑性应变差异越显著。在Monismith幂次函数模型的基础上,建立了主应力轴循环旋转条件时K0固结饱和粉质黏土累积塑性应变方程,并对模型的可靠性进行了分析和验证。     

部分排水时饱和粉质黏土变围压循环三轴试验研究

围压循环变化将引起土体孔压累积,加剧土体累积塑性应变,导致地基灾变事故。为研究围压循环变化对饱和粉质黏土变形特性的影响,利用GDS双向振动三轴仪,进行部分排水时饱和粉质黏土的变围压循环三轴试验,研究不同应力路径斜率和循环动应力比时饱和粉质黏土的孔压、轴向累积塑性应变和体应变发展规律,建立部分排水条件时考虑循环围压和循环动应力耦合作用的地基粉质黏土累积应变数学模型。试验结果表明:部分排水时饱和粉质黏土孔压、轴向累积塑性应变和体应变均随循环动应力比和应力路径斜率的增大而增大;动孔压比随着振动次数的增加明显分为急剧增加、快速下降和持续平稳三个阶段;轴向累积塑性应变和体应变均随振动次数的增加而增大。当加载次数超过2500次时,粉质黏土孔压趋于平稳,变形速率略有降低,但变形持续增加;补充试验结果表明,所建立模型与试验结果具有较高的一致性。研究结果将为交通循环荷载导致的地基灾变控制技术提供理论依据。     

Anisotropic bounding-surface plasticity model for the cyclic shakedown and degradation of saturated clay

The cyclic behaviours of embedded offshore structures under different cyclic loading levels are related to the cyclic shakedown and degradation of the surrounding soils. In the present study, a damage-dependent bounding-surface model based on a newly proposed hardening rule was developed to predict the cyclic shakedown and degradation of saturated clay and the effect of the initial anisotropic stress state. By extending the Masing’s rule to the bounding-surface plasticity theory, the stress reversal point is taken as the generalised homological centre of the bounding surface. With movement of the generalised homological centre, at lower stress amplitudes, the cyclic process ends at a steady state, and cyclic shakedown is reached. At higher stress amplitudes, a damage parameter related to the accumulated deviatoric plastic strain is incorporated into the form of the bounding surface, which is hence able to contract to model degradations in stiffness and strength. To take into account the effects of initial anisotropic conditions on the cyclic behaviour of soils, an initial anisotropic tensor is introduced in the bounding surface. The developed model is validated through undrained isotropic and anisotropic cyclic triaxial tests in normally consolidated and overconsolidated saturated clay under both one-way and two-way loadings. Both cyclic shakedown and degradation are well reproduced by the model, as is the anisotropy effect induced by the initial anisotropic consolidation process.     

循环荷载作用下饱和软土的各向异性边界面模型

在临界状态弹塑性力学的框架内,建立了可以考虑循环荷载作用下各向异性对饱和软土力学特性影响的边界面塑性模型。该模型采用非关联的流动法则,引入了反映土体各向异性的内变量,利用该内变量的初始值描述初始各向异性,采用一种理论更为严谨、模型参数确定更为恰当的旋转硬化法则描述循环加载过程中各向异性的演化,利用更新映射中心的径向映射法则和与塑性偏应变路径长度有关的塑性模量插值规律,保证模型能够模拟循环加载时应力-应变响应的非线性、滞回性、应变累积性等基本特性,解释了模型参数的物理意义和确定方法,特别是给出了一种合理确定描述土体初始各向异性状态变量方法。通过文献中等压固结和偏压固结饱和黏土的循环三轴试验结果与模型预测结果的对比验证了模型的合理性。     

适用于砂土循环加载分析的边界面塑性模型

基于临界状态土力学框架,建立了一个适用于砂土排水循环加载的边界面塑性模型。采用了考虑虚拟峰值应力比的偏应变硬化准则,初始加载阶段应力点位于边界面上,反向加载阶段以历史最大屈服面作为边界面,同时实现了对密砂软化现象的模拟和对历史所受最大应力的记忆。边界面采用修正的椭圆形,引入考虑密度与应力水平的状态相关剪胀函数,采用非相关联流动法则和以应力反向点作为映射中心的径向映射准则。模型仅有10个参数,通过常规三轴试验即可确定,并且使用一套参数可以模拟不同围压、密度的单调和循环加载情况。分别对饱和砂土的单调、循环排水三轴试验进行模拟,结果表明,该模型能够合理地反映饱和砂土排水条件下的应力-应变特性。     

Kinematic hardening of soft clay in simple shear

In a separate paper, the authors have proposed a normalized, non-degrading form of the shear stress–shear strain relationship for undrained, cyclic simple shear of soft clay. This relationship is described in the present paper, and it is seen to include a single fatigue parameter—the mean effective stress. Application of the relationship therefore requires knowledge of the history of the mean effective stress during any loading history. The present paper proposes an effective stress path model which may be used for prediction of this history. The model is developed within the framework of bounding surface kinematic and isotropic hardening plasticity. It incorporates an isotropic hardening bounding surface, and a kinematic hardening yield surface, in which the elastic region vanishes, and so the yield surface reduces to the stress point. The normalized shear stress–shear strain relationship, developed on the basis of Iwan's model, is used to establish the shape of the cap of the bounding surface. A new translation rule is also incorporated in the model, allowing improved prediction of stress path development within the bounding surface during regular or irregular cyclic loading. Use of the proposed model to simulate the behaviour of soft clay in laboratory undrained cyclic simple shear tests shows excellent qualitative agreement, with most of the major features of the actual behaviour being predicted.     

Implicit and explicit integration schemes in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay

Two integration algorithms, namely the implicit return mapping and explicit sub-stepping schemes, are adopted in the anisotropic bounding surface plasticity model for cyclic behaviours of saturated clay and are implemented into finite element code. The model is a representative of a series of bounding surface models that have typical characteristics, including isotropic and kinematic hardening rules and a rotational bounding surface to capture complex but important cyclic behaviours of soils, such as cyclic shakedown and degradation. However, there is no explicit current yield surface in the model to which the conventional implicit algorithm returns the stress state back or the sub-stepping integration corrects the drift of the stress state. Hence, necessary modifications have been made for both of the integration schemes. First, the image stress point is mapped or corrected to the bounding surface instead of mapping back or correcting the stress state to the yield surface. Second, the unloading–loading criterion is checked to determine the image stress point rather than checking the yield criterion after giving the trial stress state in a conventional way. Comparative studies on the accuracy, stability and efficiency of the two integration schemes are conducted not only at the element level but also in solving boundary value problems of monotonic and cyclic bearing behaviours of rigid footings on saturated clay. For smaller strain increments, there is no significant difference in the accuracy between the two integration schemes, but the explicit integration shows a higher efficiency and accuracy. For relatively larger increments, the implicit return mapping algorithm presents good accuracy and more robustness, while the sub-stepping algorithm shows deteriorating accuracy and suffers the convergence problem. With the tolerance used in the present model, the bearing capacity of the rigid footing predicted by the return mapping algorithm is closer to the available analytical and numerical solutions, while the bearing capacity predicted by the sub-stepping algorithm shows a marginal increase.     

循环荷载作用下接触面的边界面模型研究

基于边界面的概念,建立了循环荷载作用下砂与结构物接触面的边界面塑性模型,以描述其剪切应力与剪应变的滞回特性。该模型考虑了初始相对密度和法向应力的影响,同时在弹性模量中引入损伤因子,反映了接触面在循环荷载下逐渐硬化的性质。通过比较理论计算与试验结果,显示了该模型的合理性。