堆石料的临界状态与考虑颗粒破碎的本构模型

高应力水平时堆石料的颗粒破碎对其强度和变形机制有重要影响。临界状态土力学理论对重塑土的应力-应变关系的描述较为成功,但目前颗粒破碎对堆石料的临界状态的影响及其数学描述鲜有研究。对堆石料进行了固结应力从0.4 MPa到4 MPa的18组固结排水和固结不排水常规三轴压缩试验,以及6组等向压缩试验。试验结果表明：在排水条件和不排水条件下,不同的固结应力试样都趋于临界状态;堆石料的临界状态在q-p′平面和e-lgp′平面均为非线性变化。基于此试验结果,通过引入状态参数,在广义塑性力学的理论框架下,建立了考虑颗粒破碎的堆石料本构模型,并给出了模型参数的确定方法。与长河坝料的试验进行了对比,结果表明所建议的本构模型可以较好地模拟堆石料从低围压0.4 MPa到高围压 3.5 MPa下的应力-应变特性     

循环荷载作用下考虑颗粒破碎的堆石体本构模型

中国西部兴建的很多200～300 m高的堆石坝处于高烈度地震区。应力水平高时堆石体的颗粒破碎对其在循环荷载作用下的应力、应变特性有重要的影响。基于广义塑性理论,通过引入状态参数,建立了循环荷载作用下考虑颗粒破碎的堆石体的本构模型,并给出了模型参数的确定方法。与堆石体在400、800、1 500、2 200 kPa围压作用时的试验结果对比,表明所提出的本构模型可以较好地模拟循环荷载作用下颗粒破碎时堆石体的动应力和动应变响应。     

Structural features and mechanical behaviour of a pyroclastic weak rock

The results of an extensive programme of laboratory testing on intact and reconstituted samples of a pyroclastic weak rock from the volcanic complex of the Colli Albani (Central Italy) are presented. The deposit is known as Pozzolana Nera and may be assimilated to a bonded coarse grained material. The nature of bonds and the micro‐structural features were examined by means of diffractometry, optical and electron microscopy. As bonds are made of the same constituents of grains and aggregates of grains, bond deterioration and particles breakage upon loading are indistinguishable features of the mechanical behaviour. The testing programme consisted mainly of one‐dimensional and drained and undrained triaxial compression tests in a wide range of confining pressures up to 58 MPa. As confining stress increases, the mechanical behaviour of the material changes from brittle and dilatant to ductile and contractant; for both brittle and ductile behaviour failure is associated with the formation of shear surfaces separating the sample in several parts at the end of test. The experimental stress–dilatancy relationships are compared with the classical stress–dilatancy theories for a purely frictional material and for a material with friction and cohesion between particles. The analysis of the data indicates that peak strength results from the interplay between degradation of inter‐particle bonds, increasing friction between particles and increasing rate of dilation. Copyright © 2001 John Wiley & Son, Ltd.     

A combined method to model grain crushing with DEM

This paper presents a combined method to model grain crushing effects with discrete element method.This method combines the two most commonly used concepts to model grain crushing in DEM,i.e.the replacement method and the agglomerate method,so that it is both accurate and efficient.The method can be easily implemented.The performance is shown by several DEM simulations of biaxial tests.Particles with different crush-abilities are modeled.DEM simulation results with and without grain crushing are compared and discussed.The change of grain size distribution due to grain crushing is also investigated.     

A double structure generalized plasticity model for expansive materials

The constitutive model presented in this work is built on a conceptual approach for unsaturated expansive soils in which the fundamental characteristic is the explicit consideration of two pore levels. The distinction between the macro‐ and microstructure provides the opportunity to take into account the dominant phenomena that affect the behaviour of each structural level and the main interactions between them. The microstructure is associated with the active clay minerals, while the macrostructure accounts for the larger‐scale structure of the material. The model has been formulated considering concepts of classical and generalized plasticity theories. The generalized stress–strain rate equations are derived within a framework of multidissipative materials, which provides a consistent and formal approach when there are several sources of energy dissipation. The model is formulated in the space of stresses, suction and temperature; and has been implemented in a finite element code. The approach has been applied to explaining and reproducing the behaviour of expansive soils in a variety of problems for which experimental data are available. Three application cases are presented in this paper. Of particular interest is the modelling of an accidental overheating, that took place in a large‐scale heating test. This test allows the capabilities of the model to be checked when a complex thermo‐hydro‐mechanical (THM) path is followed. Copyright © 2005 John Wiley & Sons, Ltd.     

关于火山碎屑岩熔结过程的讨论及国内外研究现状综述

对国内外关于火山碎屑物熔结过程的研究历史做出了简要回顾,对与熔结过程相关的术语进行了归纳和说明,进一步探讨了熔结过程中出现的结构构造.此外,分别从野外调查和室内实验两方面总结了火山碎屑的熔结机理,对其研究意义做出主要说明,以期引起国内对火山作用中的火山碎屑物及碎屑流熔结过程研究的重视,以便更好地理解火山的喷发作用和重建火山喷发历史等火山学研究中的重要问题.     

Stress–dilatancy based modelling of granular materials and extensions to soils with crushable grains

Stress–dilatancy relations have played a crucial role in the understanding of the mechanical behaviour of soils and in the development of realistic constitutive models for their response. Recent investigations on the mechanical behaviour of materials with crushable grains have called into question the validity of classical relations such as those used in critical state soil mechanics. In this paper, a method to construct thermodynamically consistent (isotropic, three‐invariant) elasto‐plastic models based on a given stress–dilatancy relation is discussed. Extensions to cover the case of granular materials with crushable grains are also presented, based on the interpretation of some classical model parameters (e.g. the stress ratio at critical state) as internal variables that evolve according to suitable hardening laws. Copyright © 2004 John Wiley & Sons, Ltd.     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.     

A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.     

Integration and calibration of a plasticity model for granular materials

A new macroscopic constitutive model for non‐cohesive granular materials, with the focus on coarse‐sized materials (railway ballast), is presented. The model is based on the concepts of rate‐independent isotropic plasticity. The Backward Euler rule is used for integrating the pertinent evolution equations. The resulting incremental relations are solved in the strain space that is extended with the internal (hardening) variables. The model is calibrated using data from Conventional Triaxial Compression (CTC) tests, carried out at the University of Colorado at Boulder. A function evaluation method is used for the optimization, whereby a ‘multi‐vector’ strategy for choosing the appropriate start vector is proposed. Copyright © 2002 John Wiley & Sons, Ltd.     

滨海软土非线性流变模型及其工程应用研究

基于室内一维固结蠕变试验,研究了天津滨海淤泥质软黏土的蠕变变形规律,建立了反映滨海软黏土的非线性流变本构模型,确定了模型中的各个参数。通过此本构模型推导出天津软黏土的沉降计算公式,利用该公式对天津滨海地区北塘水库大坝软土堤基进行了沉降计算和分析,计算值与实测值比较吻合,证明了流变模型的合理性,同时也说明在软土地基沉降计算中,考虑非线性流变性状的影响是必要的。并对天津滨海新区软土长期蠕变沉降作出预测。     

Application of unconstrained optimization and sensitivity analysis to calibration of a soil constitutive model

Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

软土流变模型及在建筑物纠偏工程中的应用

根据地基应力解除法纠偏方法原理和软土的变形机理,建立了土体非线性弹-粘弹流变本构模型,并将其成功运用于软粘土地基应力解除法纠偏工程实践.结合某6层住宅楼的倾斜扶正过程,运用所建的本构模型和有限元技术进行了全过程的数值模拟,揭示了建筑物倾斜过程和扶正过程中地基应力-应变的变化规律,并预测了纠偏施工后建筑物的附加沉降量,所得结果与实测结果有较好的吻合.研究表明,土的流变特性对软土地基建筑物的纠偏实践有明显的影响,在实际工程中应予以充分考虑.     

A destructuration theory and its application to SANICLAY model

Many natural clays have an undisturbed shear strength in excess of the remoulded strength. Destructuration modeling provides a means to account for such sensitivity in a constitutive model. This paper extends the SANICLAY model to include destructuration. Two distinct types of destructuration are considered: isotropic and frictional. The former is a concept already presented in relation to other models and in essence constitutes a mechanism of isotropic softening of the yield surface with destructuration. The latter refers to the reduction of the critical stress ratio reflecting the effect of destructuration on the friction angle, and is believed to be a novel proposition. Both the types depend on a measure of destructuration rate expressed in terms of combined plastic volumetric and deviatoric strain rates. The SANICLAY model itself is generalized from its previous form by additional dependence of the yield surface on the third isotropic stress invariant. Such a generalization allows to obtain as particular cases simplified model versions of lower complexity including one with a single surface and associative flow rule, by simply setting accordingly parameters of the generalized version. A detailed calibration procedure of the relatively few model constants is presented, and the performance of three versions of the model, in descending order of complexity, is validated by comparison of simulations to various data for oedometric consolidation followed by triaxial undrained compression and extension tests on two structured clays. Copyright © 2009 John Wiley & Sons, Ltd.     

The rate of sedimentation from turbulent suspension: An experimental model with application to pyroclastic density currents and discussion on the grain‐size dependence of flow runout

Large‐scale experiments generating ground‐hugging multiphase flows were carried out with the aim of modelling the rate of sedimentation, of pyroclastic density currents. The current was initiated by the impact on the ground of a dense gas‐particle fountain issuing from a vertical conduit. On impact, a thick massive deposit was formed. The grain size of the massive deposit was almost identical to that of the mixture feeding the fountain, suggesting that similar layers formed at the impact of a natural volcanic fountain should be representative of the parent grain‐size distribution of the eruption. The flow evolved laterally into a turbulent suspension current that sedimented a thin, tractive layer. A good correlation was found between the ratio of transported/sedimented load and the normalized Rouse number of the turbulent current. A model of the sedimentation rate was developed, which shows a relationship between grain size and flow runout. A current fed with coarser particles has a higher sedimentation rate, a larger grain‐size selectivity and runs shorter than a current fed with finer particles. Application of the model to pyroclastic deposits of Vesuvius and Campi Flegrei of Southern Italy resulted in sedimentation rates falling inside the range of experiments and allowed definition of the duration of pyroclastic density currents which add important information on the hazard of such dangerous flows. The model could possibly be extended, in the future, to other geological density currents as, for example, turbidity currents.     

基于非线性损伤理论的改进CVISC模型及其在FLAC3D中实现

滑坡形成是一个典型的岩土变形破坏时效过程,低速滑坡时效过程尤为显著。基于流变理论,建立反映滑坡变形破坏发展过程的本构模型、预测滑坡活动趋势,一直是国内外滑坡研究的基本途径和难点问题之一。然而,目前国内外已有的多数流变模型仅能反映滑坡岩土蠕变的第一、第二阶段,不能刻画滑坡岩土蠕变的第三阶段（加速蠕变阶段）。本文基于滑带在滑坡中的作用以及滑带岩土剪切蠕变发展的累进性和非线性特征,借助损伤理论,在FLAC3D内置的CVISC流变本构模型中引入非线性损伤黏塑性元件,构建了可描述滑坡加速蠕变过程的非线性损伤流变本构模型,依据类比原理建立了改进的CVISC三维差分模型,通过FLAC3D开放的用户接口实现了本构模型的二次开发,并将改进的CVISC模型用于长期缓慢滑移、伴随间歇性剧烈活动的甘肃舟曲泄流坡数值模拟中。模拟结果显示,该模型不仅呈现了滑坡的加速蠕变特征,而且揭示的滑坡活动特征与其曾经出现的活动特征基本一致,从而证实基于非线性损伤理论的改进CVISC模型具有较好的有效性。     

粗粒土三轴试验数值模拟与试样颗粒初始架构初探

试样颗粒初始架构是粗粒土室内三轴试验中一个难以人为控制的因素。基于元胞自动机方法,结合粗粒土的室内三轴试验开发了能制备不同颗粒初始架构粗粒土试样的HHC-CA模型,该模型制备的粗粒土试样表征了粗粒土各粒组分布的不均匀性和随机性,再借助FLAC3D进行粗粒土三轴数值模拟试验,并探讨了试样剪切带内的砾石含量和试样砾石含量对内摩擦角的影响。数值模拟结果表明,同一粒径级配下,粗粒土的内摩擦角值与剪切带内的砾石含量呈增函数关系;试样砾石含量的增加可能会出现内摩擦角减小的现象,但同一粒径级配下的内摩擦角均值随试样砾石含量的增加而明显增大。     

湖相软土流变模型识别及其工程应用分析

对苏州湖相软土进行了三轴流变试验,根据试验结果拟采用S-6、K-B、K-S 3个流变模型来描述土体的流变特征,在模型参数反演中发现,S-6模型对加载初期变形的拟合精度不够,K-B模型对加载中期变形的拟合精度不够,而K-S模型能够很好地描述苏州湖相软土的流变特征。在模型的实际应用中引入应力比的概念,发现应力比与模型参数成指数关系,以此可以根据不同的工程应力条件来选择模型参数,最后,将该模型用于实际工程的沉降计算,通过与实测数据的对比显示模型的拟合精度较高,可以用于指导该地区的工程建设。     

Induced fabric under cyclic and rotational loads in a strain space multiple mechanism model for granular materials

The strain space multiple mechanism model idealizes the behavior of granular materials on the basis of 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 dependent on the contact distribution function as well as the normal and tangential components of interparticle 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 contact distribution and interparticle contact forces) to the macroscopic stress. The stress space fabric produces macroscopic stress through the tensorial average. Thus, the stress space fabric characterizes the fundamental and higher modes of anisotropy induced in granular materials. Herein, the induced fabric is associated with monotonic and cyclic loadings, loading with the rotation of the principal stress, and general loading. Upon loading with the rotation of the principal stress axis, some of the virtual simple shear mechanisms undergo loading whereas others undergo unloading. This process of fabric evolution is the primary cause of noncoaxiality between the axes of principal stresses and strains. Although cyclic behavior and behavior under the rotation of the principal stress axis seem to originate from two distinct mechanisms, the strain space multiple mechanism model demonstrates that these behaviors are closely related through the hysteretic damping factor. Copyright © 2011 John Wiley & Sons, Ltd.