堆石料的剪胀特性与广义塑性本构模型

以三轴试验成果为基础,考虑颗粒破碎引起堆石料剪胀比与应力比之间的非线性关系,提出了能够反映堆石料低围压剪胀、高围压剪缩特性的剪胀方程。在广义塑性理论框架内构造堆石料的塑性流动方向向量和加载方向向量,引入依赖于密实度与平均应力的压缩参数,构造随平均应力、剪应力比和密实度变化的塑性模量,建立了一个考虑颗粒破碎的堆石料弹塑性本构模型。阐述了该模型10个参数的确定方法,并通过模拟不同围压和不同应力路径下堆石料的三轴压缩试验资料验证了模型与参数的合理性。     

堆石料状态相关弹塑性本构模型

堆石料的强度变形特性与初始孔隙及应力状态等因素相关。建立了能够预测不同初始孔隙与初始围压影响的堆石料弹塑性本构模型。在剑桥类本构模型框架内,模型能够反映随着孔隙与围压的增大,变形特性由剪胀趋于剪缩的规律。模型采用了基于颗粒体材料细观结构变化的屈服函数和非关联流动准则,提出了能够反映堆石料正常固结线不唯一的硬化参数。为了反映状态相关性,假定堆石料存在唯一的临界状态面,探讨了考虑状态相关性需要满足的数学条件,从而对剪胀方程与硬化参数进行了修正。提出了基于粒子群优化算法的模型参数快速确定方法,将某筑坝堆石料不同初始孔隙比与围压条件下模型预测结果与三轴试验结果对比,验证了模型的合理性。     

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

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

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

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

不同加载方向条件下堆石料剪胀特性试验研究

剪胀方程描述的是塑性应变增量比例大小的变化规律,是构建弹塑性模型的重要组成部分。然而目前堆石料的剪胀方程大都是根据等围压 路径得到的,这些剪胀方程是否能反映不同应力路径条件下堆石料的剪胀规律还尚不明了。开展了等围压 、等轴向力 、等平均主应力p、等主应力比 共4种应力路径下的堆石料大型三轴试验,研究了应力路径对堆石料剪胀规律的影响。试验表明：不同路径的剪胀比 和应力比 均呈近似线性关系,但剪胀比与应力比不存在一致的关系,表现出明显的应力路径相关性;在 - 空间中,剪胀线随着应力增量方向的变化会发生平移和转动,其斜率 和截距参数 均随 ( )的增大而增大;目前的剪胀理论不能较好地描述剪胀的应力路径相关性,提出了一个考虑应力增量方向影响的修正剪胀方程并对其进行了试验验证。     

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

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

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

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

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

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

堆石料Gudehus-Bauer亚塑性本构模型改进及参数确定方法

研究了Gudehus-Bauer亚塑性本构模型和模型参数的求取方法。采用侧限压缩试验曲线求取模型参数颗粒硬度hs和指数n。根据模型方程的推导,建立了拟合指数 和 与围压之间的关系,并提出了新的拟合参数。考虑到堆石料具有明显的剪胀、剪缩性,在Gudehus-Bauer模型线性项中增加了主要控制体积应变项 ,以改善模型对堆石料体积应变曲线的描述。采用堆石料大型侧限压缩试验、常规三轴试验分别验证了新的拟合参数和改进后的Gudehus-Bauer亚塑性本构模型。与堆石料试验成果比较,提出的新拟合参数与改进后的Gudehus-Bauer亚塑性本构模型可以较好地模拟堆石料的应力-应变特性,并较好地改善了堆石料体变曲线的模拟结果。对改进后的模型作了常规三轴加、卸载模拟,模拟结果反映了改进的Gudehus-Bauer亚塑性本构模型具有一定的卸载适应性。     

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

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

一个用于面板坝流变分析的堆石流变模型

在前人工作的基础上, 从宏观上建立了一个用于面板坝流变分析的堆石流变模型: 堆石的瞬时弹塑性变形由椭圆-抛物双屈服面模型确定,并计入时间效应, 由双曲线函数经验公式计算堆石与时间有关的粘塑性变形。通过现场实测结果的反分析验证了模型的可靠性。     

掺复合浆液堆石料压缩特性试验研究

针对高堆石坝坝体变形大且持续时间长的问题,提出了一种提高堆石体压缩模量、减缓流变变形的方法,即利用易流动、后期硬化的高掺粉煤灰水泥砂浆（简称复合浆液）在堆石体摊铺过程中充填堆石体孔隙,在施工碾压的基础上进一步减小堆石体孔隙率。首先,通过观测复合浆液在透明玻璃球中的流动性态,研究不同配比与掺量对复合浆液流动特性的影响;然后,通过侧限压缩和流变试验,研究掺复合浆液堆石料的压缩模量、流变变形的变化规律。试验结果表明：复合浆液在颗粒体孔隙中的流动性能与粉煤灰掺比、水胶比、砂粒最大粒径以及掺量有关;复合浆液的掺入能有效提高堆石料压缩模量,减少流变变形,且对于软岩料、级配不良料效果更为明显;复合浆液在堆石料中主要起到了充填、胶结与润滑作用。     

An elastoplastic constitutive model for rockfills incorporating energy dissipation of nonlinear friction and particle breakage

A new elastoplastic model is developed for rockfills within the general critical state framework incorporating the state parameter. Two state functions are proposed to characterize the evolution of volume dilation and strain softening of rockfills, and a modified breakage index based on the concept of Hardin's relative breakage is defined to describe the progressive crushing of rockfills. The nonassociated plastic flow rule is derived from a state dependent dilatancy equation, and it incorporates energy dissipation due to intrinsic nonlinear friction and particle breakage upon shearing. Thus, their couple effect on the plastic deviatoric and volumetric deformation of rockfills is taken into account in the current model. The numerical analyses are carried out for a series of drained triaxial tests on the modeled rockfills at various consolidation pressures and stress paths. The volume dilation/contraction and strain softening/hardening of rockfills are accurately predicted by the proposed model, and the particle breakage and nonlinear critical state shear strength of rockfills are also well captured. The research findings indicate that the current model is applicable to represent the complex stress–strain–volume change behavior of rockfills in general. Copyright © 2013 John Wiley & Sons, Ltd.     

Bounding surface plasticity model for stress-strain and grain-crushing behaviors of rockfill materials

Crushing of grains can greatly influence the strength,dilatancy,and stress-strain relationship of rockfill materials.The critical state line(CSL)in the void ratio versus mean effective stress plane was extended to the breakage critical state plane(BCSP).A state void-ratio-pressure index that incorporated the effect of grain crushing was proposed according to the BCSP.Rowe’s stress-dilatancy equation was modified by adding the breakage voidratio-pressure index,which was also incorporated into the formulations of the bounding stress ratio and plastic modulus.A BCSP-based bounding surface plasticity model was proposed to describe the state-dependent stressstrain behaviors and the evolution of grain crushing during shearing process of rockfill materials,and was shown to sufficiently capture the breakage phenomenon.     

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.

     

坝体填筑料压缩特性及影响因素分析

堆石料的压缩性直接关系工程的安全问题,对某一抽水蓄能电站坝体填筑料开展多组大型侧限压缩试验,结果表明,堆石料的压缩模量较大、压缩指数很小。只有在应力超过颗粒破碎应力以后堆石料压缩曲线才会发生较大程度的下降,且认为颗粒破碎应力是堆石料发生弹塑性变形的转折点,对于不同初始孔隙比的同一种堆石料,在极高的应力状态下它们的压缩曲线会聚到唯一的极限压缩曲线LCC。压缩曲线中堆石料的卸荷回弹量非常小,表明其是弹性变形。不论何种加载路径材料的初始级配决定了堆石料的最终压缩状态,且随着母岩强度的增加、初始孔隙比或者泥岩含量的减小堆石料的压缩性减小。     

Unified 3D critical state bounding-surface plasticity model for soils incorporating continuous plastic loading under cyclic paths. Part I: Constitutive relations

The unified three-dimensional (3D) critical state bounding-surface plasticity model gUTS enables clays, silts and sands to be treated within a single framework. Furthermore, loose and dense states of a particular soil subjected to a wide range of confinements are viewed as a single material defined by the same set of constants. The model is able to handle both monotonic and complex cyclic paths including those involving a rotation of the principal stress directions. The model incorporates the following features: combined use of radial and deviatoric mapping rules and the use of an apparent normal consolidation line for sands; use of a non-associated flow rule where the ratio of the rates of volumetric plastic strain to deviatoric plastic strain is a function only of the ratio of deviatoric to mean effective stresses and the Lode angle; adoption of a bi-linear critical state line projected onto the plane of the void ratio versus logarithm of mean effective stress; inclusion of a sub-elliptic, or super-elliptic, segment in the plastic dilatancy surface for stress ratios less than critical; use of elliptic segments in the deviatoric planes; movement of the projection centre in the deviatoric mapping region and incorporation of a plastic stiffening effect for cyclic paths which repeatedly load in the same deviatoric direction.     

A reliability-based approach to evaluating the stability of high rockfill dams using a nonlinear shear strength criterion

In recent years, many high rockfill dams have been constructed in China for the purpose of hydropower generation. One of the critical aspects of rockfill dam design is the analysis of slope stability. Triaxial compression tests show that the failure envelopes of rockfills are curved and that nonlinear shear strength criteria yield better predictions of the shear strengths of rockfills than the Mohr–Coulomb criterion. Because the determination of shear strength parameters involves uncertainties, a reliability-based methodology was developed for use in evaluating the failure probability of rockfill dam slopes by integrating Bishop’s simplified method, Duncan’s nonlinear strength criterion and the first-order second-moment reliability method presented in this paper. A computer program, SCU-SLIDE, was developed and its outputs validated by comparison with Monte Carlo simulation results. The approach developed was used to study the stability of the Shuang Jiang Kou rockfill embankment dam, which when completed will be the tallest dam in the world. The results of the stability analysis are discussed and conclusions are presented in this paper.