Testing and modeling of the state-dependent behaviors of rockfill material

Dense Tacheng rockfill material (TRM) exhibits strain softening and dilation during drained triaxial tests, and therefore, an adapted Rowe’s stress–dilatancy equation was proposed for TRM. This equation incorporates an internal state index related to the density and pressure, as well as the coefficient of particle breakage and rotation. The adapted Rowe’s stress–dilatancy equation indicates that the relationship between stress and dilatancy is not constant, but varies with density and pressure. This result is in agreement with TRM test data. A state-dependent model was established for TRM using generalized plasticity theory combined with the adapted Rowe’s stress–dilatancy equation. The model includes twelve constants calibrated using TRM test data from Group A, and this model was used to predict the strain softening and dilatancy behaviors of dense TRM. Furthermore, the model predictions were validated using test data from Group B. In summary, the model accurately represented the stress–strain and dilatancy behaviors of TRM over a wide range of densities and pressures.     

DEM study on the effect of particle breakage on the macro- and micro-behavior of rockfill sheared along different stress paths

The behavior of crushable rockfill sheared along different stress paths is studied using the discrete element method. Rockfill particles are modeled as breakable agglomerates, and reasonable consistency is found between the predicted and experimental results. The simulation highlights the influence of agglomerate breakage on both the macro- and micro-behavior of the assembly. The evolution of both the number and mode of agglomerate breakage during shearing is significantly influenced by the confining pressure, deviator stress ratio, and loading direction. The relationship between the deviator stress ratio on the macro-scale and the deviator fabric on the micro-scale along different stress paths is explored.     

复杂应力路径下堆石体本构模型比较验证

目前对堆石体应力-应变关系特性的研究主要建立在常规三轴试验的基础上,对复杂应力路径上堆石体本构模型的验证工作尚不充分。利用糯扎渡高心墙堆石坝主堆石料复杂应力路径大型三轴试验成果,对国内常用的堆石体本构模型--邓肯-张EB模型、清华非线性解耦KG模型和沈珠江双屈服面模型进行了比较和验证。结果表明,邓肯-张EB模型通常会夸大堆石体的体积压缩变形,其加卸载准则无法正确判别一些复杂应力路径的加卸载状况,清华KG模型和沈珠江双屈服面模型对堆石体复杂应力路径的适应性相对较好。     

应力路径和干湿状态对堆石料颗粒破碎的影响研究

堆石料颗粒破碎是引起高土石坝变形的重要因素。在大坝填筑、蓄水期,堆石料的应力路径和干湿状态均是变化的。通过大型三轴试验,系统地研究了应力路径和干湿状态对堆石料颗粒破碎规律的影响。试验结果表明：（1）相同初始条件下,按不同应力路径达到同一轴向应变停机时测定的颗粒破碎率是不同的,等围压?3试验产生的颗粒破碎最大,等平均主应力p试验的次之,等最大主应力?1试验的最小,但不同应力路径下的颗粒级配演化规律是一致的。（2）相同初始条件下,湿样的颗粒破碎率明显高于干样,且二者的差距随着围压的增大而增大,不同干湿条件下的颗粒级配演化规律同样是一致的。（3）建立的考虑母岩强度的颗粒破碎率与塑性功的关系可以较好地统一描述不同应力路径及干湿状态下的颗粒破碎。该研究成果可为建立复杂应力路径及干湿变化条件下考虑颗粒破碎的堆石料弹塑性本构模型提供依据。     

Modeling the particle breakage of rockfill materials with the cohesive crack model

A practical combined finite–discrete element method was developed to simulate the breakage of irregularly shaped particles in granular geomaterials, e.g., rockfill. Using this method, each particle is discretized into a finite element mesh. The potential fracture paths are represented by pre-inserted cohesive interface elements (CIEs) with a progressive damage model. The Mohr–Coulomb model with a tension cut-off is employed as the damage initiation criterion to rupture the predominant failure mode occurs at the particle scale. Two series of biaxial tests were simulated for both the breakable and unbreakable particle assemblies. The two assemblies have identical configurations, with the exception that the former is inserted with CIEs and is breakable. The simulated stress–strain–dilation responses obtained for both assemblies are in agreement with experimental observations. We present a comprehensive study of the role of particle breakage on the mechanical behavior of rockfill materials at both the macroscopic and microscopic scales. The underlying mechanism of particle breakage can be explained by the force chain in the assemblies.     

堆石料等应力比路径四模量增量非线性模型

已有研究表明,土石坝内堆石料在施工填筑期的应力路径可近似为主应力比为常数的路径。根据大型三轴仪上进行的两种应力路径排水试验,即等向压缩的固结试验和等应力比值的剪切试验,提出了一个四模量增量非线性模型。模型除给出了堆石料体积模量K和剪切模量G的表达式外,对剪应力产生的体积应变和平均主应力产生的剪切应变也分别采用了交叉模量J1、J2表达。通过与试验数据的对比表明,四模量增量非线性模型能够较好地预测堆石料在等应力比路径下的应力-应变关系。     

A two-parameter expression for failure surfaces

The Mohr–Coulomb, Tresca and Von Mises criteria are classical failure criteria, widely accepted and used for various materials. To take into account the influence of intermediate principal stress on the strength of soil, Bishop, Lade–Duncan and Matsuoka proposed criteria in terms of three principal stresses or three stress invariants. This note describes an expression with two parameters for modeling the shape of the failure surface in the octahedral plane. By studying the roles of these two parameters and comparing the new criterion with the aforementioned criteria, the advantage and flexibility of the proposed function is explored. Application of the function is demonstrated by fitting the new surface to experimental data for various soils.     

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.     

复杂应力路径下堆石料本构关系研究

已有研究表明,土石坝内堆石料在坝体填筑过程的应力路径可近似为等应力比的路径（q/p=常数）,水库蓄水时应力路径将发生转折,呈复杂的应力路径形态（dq/dp=常数）。在大型三轴仪上进行了两种应力路径的排水试验,即等应力比路径下的偏压试验和复杂应力路径下的剪切试验。根据试验结果提出了一个堆石料应力路径增量非线性弹性模型,模型采用三模量形式除可以描述堆石料等应力比路径的应力-应变特征外,通过转折后的路径特征构造合适的柔度矩阵,能够表达转折应力路径下的本构关系。对试验曲线进行拟合表明,应力路径模型能够较好地反映堆石料在复杂应力路径下的应力与变形特性。     

堆石料颗粒破碎特性试验研究

颗粒破碎是堆石料的一项基本特性,它对堆石体的变形和强度特性具有明显的影响。对于高堆石坝而言,在高应力场作用下堆石颗粒发生明显破碎,可导致坝体变形率增加。为了正确认识堆石体及堆石坝的变形特性和机理,研究了堆石颗粒破碎特性以及颗粒破碎的影响因素。采用大型三轴试验研究了堆石料的颗粒破碎特性,分析了堆石体干密度、级配特征、堆石颗粒强度等对颗粒破碎的影响,研究了应力状态对颗粒破碎率的影响,建立了颗粒破碎率的计算模型以及颗粒破碎引起的堆石体应变增量与颗粒破碎率的关系。     

Dilatancy of granular materials in a strain space multiple mechanism model

A granular material consists of an assemblage of particles with contacts newly formed or disappeared, changing the micromechanical structures during macroscopic deformation. These structures are idealized through a strain space multiple mechanism model as a twofold structure consisting of a multitude of virtual two‐dimensional mechanisms, each of which consists of a multitude of virtual simple shear mechanisms of one‐dimensional nature. In particular, a second‐order fabric tensor describes direct macroscopic stress–strain relationship, and a fourth‐order fabric tensor describes incremental relationship. In this framework of modeling, the mechanism of interlocking defined as the energy less component of macroscopic strain provides an appropriate bridge between micromechanical and macroscopic dilative component of dilatancy. Another bridge for contractive component of dilatancy is provided through an obvious hypothesis on micromechanical counterparts being associated with virtual simple shear strain. It is also postulated that the dilatancy along the stress path beyond a line slightly above the phase transformation line is only due to the mechanism of interlocking and increment in dilatancy due to this interlocking eventually vanishing for a large shear strain. These classic postulates form the basis for formulating the dilatancy in the strain space multiple mechanism model. The performance of the proposed model is demonstrated through simulation of undrained behavior of sand under monotonic and cyclic loading. Copyright © 2010 John Wiley & Sons, Ltd.     

Constitutive modeling of inherent anisotropy in a strain space multiple mechanism model for granular materials

Consideration of fabric anisotropy is crucial to gaining an improved understanding of the behavior of granular materials. This paper presents a constitutive model to describe the sand behavior associated with fabric anisotropy within a framework of a strain space multiple mechanism model. In the proposed model, a second-order fabric tensor is extended by incorporating a new function that represents the effect of inherent (or initial fabric) anisotropy, along with three additional parameters: two of them, a₁ and a₂ , control the degree of anisotropy, and the second mode of inherent anisotropy can be expressed by introducing the parameter a₂ as well as the first mode by the parameter a₁ . The third parameter, θ₀ , expresses the principal direction of inherent anisotropy (eg, the normal vector direction of bedding planes relative to horizontal axis). The formulation of the dilative component of dilatancy (ie, positive dilatancy) is also extended to consider the effect of inherent anisotropy based on the interlocking mechanism. Experimental data on the complex anisotropic responses of Fraser River sand and Toyoura sand under monotonic loading is used to validate this model. The proposed model is shown to successfully capture anisotropic responses, which become contractive or dilative depending on different principal-stress directions, with a single set of anisotropy parameters; thus, the model is considered to possess the capability to simulate the anisotropic behaviors of granular materials. In addition to different loadings on the same fabric, the effects of different fabric anisotropies upon the sand behavior under the same loadings are also investigated.     

Evolution of induced fabric in a strain space multiple mechanism model for granular materials

The strain space multiple mechanism model idealizes the behavior of granular materials based on a multitude of virtual simple shear mechanisms oriented in arbitrary directions. Within this modeling framework, the virtual simple shear stress is defined as a quantity that depends on the contact distribution function as well as the normal and tangential components of inter‐particle contact forces, which evolve independently during the loading process. In other terms, the virtual simple shear stress is an intermediate quantity in the upscaling process from the microscopic level (characterized by the contact distribution and inter‐particle contact forces). The stress space fabric (i.e. the orientation distribution of the virtual simple shear stress) produces macroscopic stress through the tensorial average. Thus, the stress space fabric characterizes the fundamental and higher modes of anisotropy induced in granular materials. Comparing an induced fabric associated with the biaxial shear of plane granular assemblies obtained via a simulation using Discrete Element Method to the strain space multiple mechanism model suggests that the strain space multiple mechanism model has the capability to capture the essential features in the evolution of an induced fabric in granular materials. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

三维应力状态下粗粒料强度特性试验研究

粗粒料所处的应力状态具有随时间和空间而变的特点,往往处于较复杂的三维应力状态下。通过粗粒料大型真三轴各向等压固结等比例加载试验,研究了不同中主应力系数条件下粗粒料的强度特性。试验结果表明：中主应力对粗粒料的强度有重要影响,三维应力状态下,粗粒料的强度比常规三轴应力状态下有较大提高,大小主应力之差与大主应变关系曲线也更加陡峭;中主应力系数b从0增大到0.25时,破坏时大小主应力之差增加了39%～50%;摩尔–库仑强度参数咬合力c和内摩擦角? 均随着中主应力系数b的增大而增大,其中咬合力c的增长较为显著,特别是b从0增大到0.25时;相同小主应力条件下,b = 0时,破坏偏应力与球应力之比最小,b = 0.25时,破坏应力比达到最大值,随着b的增大破坏应力比又有所减小;相同中主应力系数条件下,随着小主应力的增大,破坏应力比逐渐减小。     

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

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

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

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

UH模型对粗颗粒土适用性验证及土石坝工程应用

对3种粗颗粒土进行了常规三轴CD试验、K₀压缩试验和等应力比路径试验,根据三轴试验数据标定UH(unified hardening,统一硬化)模型的参数,然后模拟另外两种路径下的应力-应变关系,并与邓肯E-B模型进行比较。结果表明:UH模型能合理地描述不同应力路径下粗颗粒土的力学特性,而邓肯E-B模型的预测结果则较差,验证了UH模型对粗颗粒土的适用性。最后利用UH模型对某心墙堆石坝进行了应力变形三维有限元分析,并将计算结果与邓肯E-B模型及现场监测值进行比较。对比显示:在坝体变形和应力的定性分布规律上,UH模型的计算结果与邓肯E-B模型基本一致;而在定量上,UH模型得到的坝体变形总体上更接近现场监测值,验证了UH模型在土石坝工程中的适用性和合理性。