高低压下不同力学特性的砂土统一模型

为构建能够反映砂土高低压下不同力学特性的统一模型,分析了砂土在较大的压力范围内的力学试验结果并获取其强度、等向压缩以及临界状态特性参数。通过引入应力路径相关因子来修正塑性应变增量中与应力路径相关的部分,从而使得模型硬化参量能够反映密实砂土在常压下的剪胀特性。同时,基于砂土的临界状态特性提出能够与砂土内部状态相对应的潜在状态面概念,由屈服面与潜在状态面间的动态关系确定加载过程中的动态密实参数及潜在强度,进而使得硬化参量也能够反映密实砂土在常压下的软化特性及高压下的剪缩、硬化特性。分析模型屈服面及潜在状态面间的演化规律并对不同压力等级下的砂土受荷力学行为进行模拟预测,证实了该模型能够反映密实砂土常压下剪胀软化及高压下剪缩硬化的特性。     

砂土的物态本构模型

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

广义塑性力学中的屈服面与应力-应变关系

详细讨论了广义塑性力学中屈服面和塑性势面的对应关系以及岩土材料的三类屈服面(即体积屈服面与 q 方向上及方向上的剪切屈服面)的基本特征, 尤其是提出了能考虑剪胀与剪缩的体积屈服面和应力 Lode 角θσ方向的剪切屈服面。指出在广义塑性力学中不必采用硬化定律, 就能得出塑性应变增量与应力增量的关系, 给出了求弹塑性矩阵的方法。     

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

基于适用于饱和砂土循环动力分析的边界面模型,利用修正广义Mises方法将其推广至三维应力空间中,引入与状态相关的剪胀函数,采用历史最大加载面硬化准则,以反映饱和砂土排水条件下刚度变化、剪胀剪缩等特性。结合子增量步显式积分算法的思想,建立了适用于饱和砂土循环动力分析边界面模型的显式积分算法流程。借助有限元软件子程序接口,将该算法流程开发到有限元软件中,通过建立土体单元数值模型,对Toyoura砂在单调荷载和循环荷载下的三轴试验工况进行模拟。结果表明,利用子增量步显式积分算法对模型进行积分,能够准确有效地模拟砂土的应力-应变曲线、以及砂土在单调和循环荷载作用下均展现的剪胀剪缩现象,验证了显式积分算法的有效性。通过设置不同的增量步长,验证了子增量步显式积分算法对于较小应变增量步的高度稳定性和收敛性。     

Gudehus-Bauer亚塑性本构模型研究

无黏性土的应力-应变关系可以用Gudehus-Bauer亚塑性本构模型来模拟,该模型强调应力增量的大小和方向不仅与当前应力状态有关,而且还取决于当前应变增量的大小和方向。为分析其与传统弹塑性理论的不同之处,对Gudehus-Bauer理论的线性项和非线性项进行了研究,并对不同初始孔隙比下Gudehus-Bauer亚塑性模型的应力-应变关系进行了探讨。结果表明Gudehus-Bauer亚塑性模型不用把应变分为弹性和塑性部分就能考虑不可逆变形,并能体现密砂的剪胀特性和应变软化特性以及松砂的剪缩特性和应变硬化特性。     

主应力轴变化下非共轴对砂土剪胀特性影响

传统塑性剪胀模型在描述应力比和塑性应变增量关系时都是基于共轴塑性流动法则,从而认为土体的剪胀性仅与应力比有关。大量试验结果表明,在涉及主应力轴变化的复杂应力条件下塑性流动过程中应力-应变是非共轴的,因而在分析砂土剪胀特性时非共轴是不可忽视的因素。为了研究主应力轴变化的复杂应力条件下非共轴对砂土剪胀特性的影响,利用空心圆柱仪对饱和砂土进行了一系列定轴剪切试验、纯主应力轴旋转试验以及组合加载试验。试验结果表明,不同应力路径下应力-应变非共轴都会引起剪胀曲线偏离Rowe直线,通过Gutiereez提出的考虑非共轴因子的修正剪胀方程可以修正非共轴引起的偏差,从而使得Rowe剪胀方程适用于涉及主应力轴旋转等更加复杂的加载条件。     

剪胀性砂土边界面模型的研究

剪胀性对于砂土,尤其是中密以及密实砂土,是一个非常显著的特性。相变线是剪胀性砂土的特征曲线,能够反映砂土的围压以及初时孔隙比对变形特性的影响。本文在边界面塑性理论的框架内,把相变状态参量引入到剪胀方程以及塑性硬化模量中,建立了一个能够描述砂土剪胀性以及循环特性的本构模型。本模型采用一套参量可以模拟不同初时孔隙比、不同围压、排水(或不排水)条件下单调(或循环)加载的应力-应变特性。验证表明本模型数值计算与试验结果相吻合。     

一个考虑循环荷载作用的简化模型

基于塑性硬化模量场理论和多重屈服面模型,结合各向同性硬化准则和移动硬化准则,在新的应力空间建立了一个新型不排水循环荷载作用下的多屈服面模型,并推导了一个适合三轴试验的简化的多屈服面模型。在此基础上,结合一个循环荷载作用下的动孔压模型,进行了饱和软黏土的动三轴模拟试验。结果表明,文中建立的多屈服面模型能够较好地模拟循环三轴试验、直剪试验和平面应变条件下的试验。     

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

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

土的三重屈服面应力-应变模型与SMP准则的结合

土的三重屈服面应力应变模型在国内有一定的影响力,它把土的塑性应变分成3部分,每一部分对应一个屈服面,分别为压缩屈服面、剪切屈服面和剪胀屈服面,从而在一定程度上反应了土的基本特性。但它在由三轴压缩应力状态向其它应力状态转化时,塑性系数需要重新确定,比较繁琐。而SMP准则能合理反映土的破坏特性,在国际上有一定的影响力。借用应力变换三维化方法,通过把土的三重屈服面应力应变模型和SMP准则相结合,使得原有的模型在不做任何假设的条件下,采用统一的塑性系数,由三轴压缩应力状态简单地转化到一般应力状态,并且能够考虑应变的分叉特性,能合理地预测已有的试验数据。     

Prediction of flow liquefaction instability of clean and silty sands

Adding a small amount of non-plastic silt to clean sands may lead to dramatic loss of shear strength and a noteworthy tendency toward contraction when the mechanical behavior of the mixture is compared with that of the clean host sand. Thus, simulation of the behavior of silty sands with varying fines content is still a challenging subject in geomechanics. A unified constitutive model for clean and silty sands is presented in this paper. To eliminate the factitious decrease of void ratio associated with inactive silt particles in various silty sand mixtures, the concept of equivalent void ratio is used in the model formulation instead of the global void ratio. In addition, the instantaneous soil state is expressed in terms of intergranular state parameter taking into account the combined influence of intergranular void ratio, mean principal effective stress and fines content. Then, dilatancy and plastic hardening modulus are directly linked to the intergranular state parameter. To improve the model capacity in simulation of cyclic tests, new features are added to the plastic hardening modulus. It is shown that the proposed model can reasonably reproduce the mechanical behavior as well as the onset of flow liquefaction instability of clean and silty sands using a unique set of parameters.     

SANISAND-FN: An evolving fabric-based sand model accounting for stress principal axes rotation

SANISAND is the name of a family of bounding surface plasticity constitutive models for sand within the framework of critical state theory, which have been able to realistically simulate the sand behavior under conventional monotonic and cyclic loading paths. In order to incorporate the important role of evolving fabric anisotropy, one such model was modified within the framework of the new anisotropic critical state theory and named SANISAND-F model. Yet the response under continuous stress principal axes rotation requires further modification to account for the effect of ensuing noncoaxiality on the dilatancy and plastic modulus. This modification is simpler than what is often proposed in the literature, since it does not incorporate an additional plastic loading mechanism and/or multiple dilatancy and plastic modulus expressions. The new model named SANISAND-FN is presented herein and is validated against published data for loading that includes drained stress principal axes rotation on Toyoura sand.     

考虑组构变化效应的饱和密砂循环本构模型

为了真实地描述饱和密砂在循环加载过程中的变形行为,需要引入考虑剪胀阶段组构变化的宏观参量。在已有的基于状态参量的本构模型基础上,引入反映组构变化的剪胀内变量,简称组构-剪胀内变量z。以相变线PTL作为参考线,采用基于相变的状态参量判断砂土在初始时刻和任意时刻体积变形的变化趋势,并通过z对剪胀比d的影响,考虑反向加载过程中塑性变形的累积,建立了一个针对饱和密砂的循环加载的弹塑性本构模型。该模型根据试验现象将已有模型中的塑性剪切模量区分为首次加载模量与再加载模量,能较好地模拟排水情况下砂土循环加载的胀-缩变化过程。最后,针对密砂的三轴排水情况,利用文中模型进行预测,并把预测结果与试验结果进行比较,结果表明该模型能够总体反映砂土循环加载的变形行为。     

Constitutive modelling of granular materials using a contact normal-based fabric tensor

This paper presents a fabric tensor-based bounding surface model accounting for anisotropic behaviour (e.g. the dependency of peak strength on loading direction and non-coaxial deformation) of granular materials. This model is developed based on a well-calibrated isotropic bounding surface model. The yield surface is modified by incorporating the back stress which is proportional to a contact normal-based fabric tensor for characterising fabric anisotropy. The evolution law of the fabric tensor, which is dependent on both rates of the stress ratio and the plastic strain, rules that the material fabric tends to align with the loading direction and evolves towards a unique critical state fabric tensor under monotonic shearing. The incorporation of the evolution law leads to a rotational hardening of the yield surface. The anisotropic critical state is assumed to be independent of the initial values of void ratio and fabric tensor. The critical state fabric tensor has the same intermediate stress ratio (i.e. b value) and principal directions as the critical state stress tensor. A non-associated flow rule in the deviatoric plane is adopted, which is able to predict the non-coaxial flow naturally. The stress–strain relation and fabric evolution of model predictions show a satisfactory agreement with DEM simulation results under monotonic shearing with different loading directions. The model is also validated by comparing with laboratory test results of Leighton Buzzard sand and Toyoura sand under various loading paths. The comparison results demonstrate encouraging applicability of the model for predicting the anisotropic behaviour of granular materials.

     

A critical state model for sands dependent on stress and density

An elastoplastic model for sands is presented in this paper, which can describe stress–strain behaviour dependent on mean effective stress level and void ratio. The main features of the proposed model are: (a) a new state parameter, which is dependent on the initial void ratio and initial mean stress, is proposed and applied to the yield function in order to predict the plastic deformation for very loose sands; and (b) another new state parameter, which is used to determine the peak strength and describe the critical state behaviour of sands during shearing, is proposed in order to predict simply negative/positive dilatancy and the hardening/softening behaviour of medium or dense sands. In addition, the proposed model can also predict the stress–strain behaviour of sands under three-dimensional stress conditions by using a transformed stress tensor instead of ordinary stress tensor. Copyright © 2004 John Wiley & Sons, Ltd.