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

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

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

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

单调荷载下砂土变形过程数值模拟及细观机制研究

采用离散单元法的颗粒流理论,模拟了松砂和密砂在单调荷载作用下的变形过程,研究了砂土渐进破坏过程中的宏观力学行为和细观组构参量的演化规律。采用PFC的FISH语言开发了细观组构统计程序,通过记录加载不同时刻试样的细观参量,如配位数、接触法向分布、粒间法向接触力、切向接触力等的演化,分析了砂土变形过程中细观组构变化与宏观力学响应之间的内在联系。应用表征上述量的组构参数研究了砂土的诱发各向异性,探讨了松砂剪缩、密砂剪胀的细观机制。研究成果对于揭示砂土变形的细观机制以及建立砂土的细观力学模型具有重要意义。     

基于相变状态的砂土本构模型的研究

为了很好地描述砂土的应力应变特征,通常利用以临界状态孔隙比为基础的状态参数建立弹塑性本构关系,但是对于密砂排水试验,临界状态参数很难测到。相变状态作为一特征状态,相关参数较容易测得。因此,本文定义了以相变状态孔隙比为基础的状态参数,并引入砂土的剪胀和边界应力比表达式,建立砂土的弹塑性本构模型。最后,利用本文提出的模型参数得到的密砂模拟结果与试验结果相吻合,较好地反映了密砂的应变强化和软化力学特征。     

主应力方向对砂土剪胀性的影响分析

剪胀性是土特有的一种材料属性,而准确地描述砂土的剪胀性则是建立砂土本构模型的重要基础。大量常规三轴试验发现,在以相同加载条件下剪切时密砂和松砂会表现出完全不同的剪胀性和应力-应变关系特性,说明砂土的剪胀性不仅与其所处的应力状态有关,也与其物理状态相关。状态参量理论很好地解释了砂土所处应力状态和物理状态对剪胀性的共同作用。空心扭剪三轴试验仪可以实现不同主应力方向的单调剪切试验。试验结果表明,当砂土以不同主应力方向单调剪切时,即使处于相同初始应力条件和物理状态,砂土也会表现出不同的剪胀性,说明了主应力方向也是决定砂土剪胀性的重要条件。本文通过分析试验中主应力方向对砂土剪胀性的影响,提出了一个含有主应力方向的状态参量,并建立了相应的剪胀方程。通过与试验数据的对比,验证了该方法的正确性和准确性。     

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

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

等幅剪应变下砂土循环单剪行为的离散元模拟

循环动载下土体变形呈现的复杂各向异性,本质上依赖于土体微观组构特征的演化。为了揭示非比例循环动载下土体变形的微观机制,采用离散元方法模拟砂土的循环单剪行为。应用等幅剪应变的往复加载实现循环单剪应力路径,模拟得到了砂土的循环弱化、剪胀性、非共轴性及微观组构的演化规律。微观模拟表明,在循环剪切过程,土体表现为循环弱化行为,并最终趋于塑性安定状态。而试样组构主方向倾角和组构各向异性系数也会不断增大直至一个稳定值,同时主应变率、主应力轴和组构主方向旋转角度会存在一定的差异,即非共轴现象,随着单调剪切的进行,三者会逐渐趋于一致。土体的剪胀行为表现为循环压密性,而非共轴性呈现逐渐增强的演化规律。循环剪切过程中微观组构的主方向与应力主方向逐渐趋于一致,而微观组构各向异性呈现减弱的演化趋势。     

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

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

宏-细观结合的砂土单剪试验非共轴特性分析

针对传统本构理论无法描述土体单剪试验非共轴变形的不足,采用非共轴修正模型进行改进。模型基于材料状态相关临界状态理论,采用宏-细观结合的方法,将1个新的各向异性状态变量引入本构模型来描述砂土的各向异性。考虑细观组构张量和应力张量的几何关系的变化,模型可以描述砂土在主应力轴旋转条件下材料状态的变化,材料状态变化直接导致模型的硬化规律和剪胀性发生变化,因此,模型可以描述该条件下原生向异性对砂土变形的影响。引入非共轴理论对本构模型进行修正,建立了三维非共轴各向异性模型。单剪试验的加载条件会造成主应力轴相对土体沉积面发生旋转,修正模型不但能够描述砂土在主应力轴旋转条件下其原生各向异性对变形的影响,而且可以描述主应力轴旋转造成的应力诱发各向异性对土体变形的影响,因此,该模型能够对整个单剪试验的变形规律进行描述,而且物理意义清晰。通过铝棒堆积体和Toyoura砂单剪试验验证表明,非共轴修正各向异性模型能对单剪试验的整个变形过程进行较好的模拟。     

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

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

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.     

Constitutive model for cyclic behaviour of cohesionless sands

This paper presents a constitutive model for describing the stress-strain response of sands under cyclic loading. The model, formulated using the critical state theory within the bounding surface plasticity framework, is an upgraded version of an existing model developed for monotonic behaviour of cohesionless sands. With modification of the hardening law, plastic volumetric strain increment and unloading plastic modulus, the original model was modified to simulate cyclic loading. The proposed model was validated against triaxial cyclic loading tests for Fuji River sand, Toyoura sand and Nigata sand. Comparison between the measured and predicted results suggests that the proposed modified model can capture the main features of cohesionless sands under drained and undrained cyclic loading.     

A modification to dense sand dynamic simulation capability of Pastor–Zienkiewicz–Chan model

A modification to the nonlinear Pastor–Zienkiewicz–Chan (PZC) constitutive model without any change in the number of model parameters is introduced in order to simulate stiffness degradation of dense sands at dynamic loading. The PZC model is based on generalized plasticity and was verified by good prediction of liquefaction and undrained behavior of saturated sand. The PZC is a robust model that can predict drained dynamic behavior of sands, especially stiffness increase in loose sand at reloading of dynamic loading. Yet, this model does not show stiffness degradation of dense sand at reloading. The modification is made through modifying the stress memory factor, H _DM, which is multiplied by the plastic modulus, H _L. This modification does not influence reloading behavior of loose sand. The modified PZC model is verified via results of drained cyclic tests. Two cyclic triaxial tests on loose and dense specimens, along with two cyclic plane strain tests on dense sand are utilized for validation. The model simulation shows that the modified PZC model is able to predict the stiffness degradation of dense sand at reloading well.     

Modeling cyclic shearing of sands in the semifluidized state

Liquefaction is associated with the loss of mean effective stress and increase of the pore water pressure in saturated granular materials due to their contractive tendency under cyclic shear loading. The loss of mean effective stress is linked to loss of grain contacts, bringing the granular material to a “semifluidized state” and leading to development and accumulation of large cyclic shear strains. Constitutive modeling of the cyclic stress-strain response in earthquake-induced liquefaction and post-liquefaction is complex and yet very important for stress-deformation and performance-based analysis of sand deposits. A new state internal variable named strain liquefaction factor is introduced that evolves at low mean effective stresses, and its constitutive role is to reduce the plastic shear stiffness and dilatancy while maintaining the same plastic volumetric strain rate in the semifluidized state. This new constitutive ingredient is added to an existing critical state compatible, bounding surface plasticity reference model, that is well established for constitutive modeling of cyclic response of sands in the pre-liquefaction state. The roles of the key components of the proposed formulation are examined in a series of sensitivity analyses. Their combined effects in improving the performance of the reference model are examined by simulating undrained cyclic simple shear tests on Ottawa sand, with focus on reproducing the increasing shear strain amplitude as well as its saturation in the post-liquefaction response.     

A middle surface concept (MSC) model for saturated sands in general stress space

An elastoplastic constitutive model is proposed for saturated sands in general stress space using the middle surface concept (MSC). In MSC, different features of stress–strain response of a material are divided into different pseudo‐yield surfaces. The true‐yield surface representing the true response is established by using various links between the yield surfaces. In this MSC sand model, several well‐known features of sand response are represented by three different pseudo‐yield surfaces, which are developed in a simple and straightforward way. These features include the critical state behaviour, the effects of state parameter, unloading and reloading plastic deformation, the influence of fabric anisotropy, and phase transformation line related behaviour. Finally, the model predictions and test results are compared for two different types of sands under a variety of loading conditions and good comparisons are obtained. The application of MSC to saturated sand modelling shows the versatility of MSC as a general concept for modelling stress–strain response of materials. Copyright © 2006 John Wiley & Sons, Ltd.     

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

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