Evaluation of a constitutive model for clays and sands: Part I – sand behaviour

This paper evaluates the performance of a generalized effective stress soil model for predicting the rate independent behaviour of freshly deposited sands, while a companion paper describes model capabilities for clays and silts. Most material parameters can be obtained from standard laboratory data, including hydrostatic or one‐dimensional compression, drained and undrained triaxial shear testing. A compilation of data on compression behaviour allows for estimation of compression parameters when this type of data is not available. Extensive comparisons of model predictions with measured data from undrained triaxial shear tests shows that the model gives excellent predictions of the transition from dilative to contractive shear response as the confining pressure and/or the initial formation void ratio increases. A parametric study of drained response shows that the model describes realistically the variation of peak friction angle and dilation rate as a function of confining pressure and density when compared with an empirical correlation valid for many sands. The proposed formulation predicts a unique critical state locus for both drained and undrained triaxial testing which is non‐linear over the entire range of stresses and is in excellent agreement with recent experimental data. Overall, the model provides excellent predictions of the stress–strain–strength relationships over a wide range of confining pressures and formation densities. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of Temperature on the Strength Parameters at the Plastic Domain for Unconsolidated Sandstones

Heavy and extra heavy oil reservoirs are subjected to high temperatures during thermal recovery processes like SAGD and CSS. In-situ temperature state changes throughout these operations may generate variations in the mechanical properties or strength parameters of the rock such as Young and Bulk’s modulus, internal friction angle, cohesion and dilation angle. This paper presents results about the variation of the elastic modulus and mechanical properties, in both elastic and plastic region due to changes in temperature and effective confining stress for reconstituted samples of Colombian unconsolidated sands. In order to study these changes, several drained triaxial tests are performed in a range of temperature between 50 and 230 °C and a range of effective confining stress between 0.4 and 8.2 MPa, these tests are carried out on reconstituted samples of oil sands recovered in an outcrop in the Magdalena middle valley of Colombia. The experimental results are analyzed with an analytical model based on the constitutive model of Mohr–Coulomb. The analysis is divided in two ways: First, the elastic region is analyzed calculating the variation of the elastic modulus (Young and bulk’s modulus) and the mechanical parameters (cohesion and internal friction angle) with temperature and effective confining stress, these calculus are done in the linear zone of the stress–strain curve. Then, the plastic region is analyzed using Mohr–Coulomb constitutive model with a non-associative flow rule in order to represent in a better way the dilatant/contractive behavior of the material, from this analysis is possible to obtain the behavior of the mechanical properties during plastic strains. The stress–strain curves obtained from the triaxial compression tests show that for the highest and medium effective confining stress (8.2 and 4.0 MPa), there is a hardening strain behavior after the yielding point that implies a loose sand behavior under confining conditions. Meanwhile, for lowest effective confining stress (0.4 MPa) the samples show a softening strain, a typical behavior for dense sands under confining conditions.     

Formulation of a unified constitutive model for clays and sands

This paper presents a new generalized effective stress model, referred to as MIT-S1, which is capable of predicting the rate independent, effective stress–strain–strength behaviour of uncemented soils over a wide range of confining pressures and densities. Freshly deposited sand specimens compressed from different initial formation densities approach a unique condition at high stress levels, referred to as the limiting compression curve (LCC), which is linear in a double logarithmic void ratio, e, mean effective stress space, p′. The model describes irrecoverable, plastic strains which develop throughout first loading using a simple four-parameter elasto-plastic model. The shear stiffness and strength properties of sands in the LCC regime can be normalized by the effective confining pressure and hence can be unified qualitatively, with the well-known behaviour of clays that are normally consolidated from a slurry condition along the virgin consolidation line (VCL). At lower confining pressures, the model characterizes the effects of formation density and fabric on the shear behaviour of sands through a number of key features: (a) void ratio is treated as a separate state variable in the incrementally linearized elasto-plastic formulation: (b) kinematic hardening describing the evolution of anisotropic stress–strain properties: (c) an aperture hardening function controls dilation as a function of ‘formation density’; and (d) the use of a single lemniscate-shaped yield surface with non-associated flow. These features enable the model to describe characteristic transitions from dilative to contractive shear response of sands as the confining pressure increases. This paper summarizes the procedures used to select input parameters for clays and sands, while a companion paper compares model predictions with measured data to illustrate the model capability for describing the shear behaviour of clays and sands. Copyright © 1999 John Wiley & Sons, Ltd.     

基于南水模型的钙质砂应力−应变关系模拟

钙质砂是一种海洋生物成因的特殊岩土材料,具有颗粒性状不规则、易破碎的特点,表现出与石英砂不同的力学性质。为模拟钙质砂在不同应力水平下的应力?应变关系,首先分析南水模型的不足：即无法描述应力不变情况下应变无限增长的特点和不能描述钙质砂的剪胀性;然后通过在切线模量中引入应力比与峰值应力比的比值,合理考虑了颗粒破碎的影响;提出了切线模量和切线体积比与应力比间的函数关系式,使得改进后的模型能较好地描述应变?体变关系。同时,将孔隙压力系数引入到孔隙压力表达式中,较好地模拟了三轴不排水剪切试验结果。与现有的模型相比,改进后的模型表达式简单易懂,参数获取更为方便。     

Acoustic emission characteristics of a dry sandy soil subjected to drained triaxial compression

Acoustic emission (AE) technique that is capable of diagnosing the failure process of stressed materials has rarely reported its application to sandy soils subjected to triaxial compression. In this paper, drained triaxial compression tests incorporating with a high-performance AE measurement system were conducted for dry sands with different confining stresses and initial relative densities. Generally, an increased confining stress or initial relative density generates more acoustic emissions, while there also exist exceptions due to different failure patterns. A good resemblance between stress–strain and AE hit rate–strain relations was observed, and power functions between the mechanical parameters and AE hit rate were well established regardless of different confining stresses and initial relative densities. Besides, the behavioral state of yield and peak during compression could be also evaluated by AE hit rate, compared with conventional stress–strain determination. Particularly, the peak AE hit rate is found not always synchronous to but fluctuating at around the peak stress depending on different failure patterns, which might provide beneficial insights into the incipient failure of sands. The present good consistencies suggest that AE characteristics could be used as alternative parameters to evaluate and even predict the mechanical behavior of dry sands.

     

A simple three‐dimensional distinct element modeling of the mechanical behavior of bonded sands

This paper presents a simple three‐dimensional (3D) Distinct Element Method (DEM) for numerical simulation of the mechanical behavior of bonded sands. First, a series of micro‐mechanical tests on a pair of aluminum rods glued together by cement with different bond sizes were performed to obtain the contact mechanical responses of ideally bonded granular material. Second, a 3D bond contact model, which takes into account the influences of bond sizes, was established by extending the obtained 2D experimental results to 3D case. Then, a DEM incorporating the new contact model was employed to perform a set of drained triaxial compression tests on the DEM bonded specimens with different cement contents under different confining pressures. Finally, the mechanical behavior of the bonded specimens was compared with the available experimental results. The results show that the DEM incorporating the simple 3D bond contact model is able to capture the main mechanical behavior of bonded sands. The bonded specimen with higher cement content under lower confining pressure exhibits more pronounced strain softening and shear dilatancy. The peak and residual strengths, the apparent cohesion and peak/residual friction angles, and the position and slope of the critical state line increase with increase in cement content. Microscopically, bond breakage starts when the system starts to dilate and the maximum rate of bond breakage coincides with the maximum rate of dilation. Bond breakage is primarily due to tension‐shear failure and the percentage of such failures is independent of both confining pressure and cement content. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

Pastor-Zienkiewicz状态相关本构模型及其参数确定方法研究

当考虑无黏性土初始密实度及围压时,Pastor-Zienkiewicz本构模型需采用多组参数模拟同种材料,给其在实际中的应用带来了困难。针对原始Pastor-Zienkiewicz本构模型存在的这一问题,将新近发展的状态相关参数引入剪胀方程、塑性势方向矢量、加（卸）载方向矢量及塑性模量方程中,提高了模型在考虑初始密实度及初始围压时模拟无黏性土力学特性的能力;详细探讨了修正模型参数的拟定方法,明确了各参数的物理意义;通过试验资料对修正模型的预测能力进行了验证,预测结果与试验资料吻合良好,表明了修正模型的可行性和合理性。     

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

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

A bonded-particle model for cemented sand

We propose an extension of the Discrete Element Method for the numerical simulation of cemented sands, in which spherical particles are bonded together by elastic beams connecting the centers of the spheres. The parameters of this model are the strengths and stiffnesses of the bonds and particles. For small strains, the elasticity of the bond element is equal to the well-known linear finite-element Timoshenko beam element with reduced integration. The finite rotations are represented by unit quaternions. An efficient way to compute relative rotations and to decompose them into their components is presented.The results of triaxial compression tests on artificially cemented sands are used to verify that the model can capture the macroscopic behavior of such materials. The results show that peak stress mainly depends on the strength of the bonds and the number of initially bonded particles in the material. Results of triaxial tests with different cement contents are reproduced by the analysis. An important parameter of the model is the strength difference between tension and compression of the bond element. This property controls the influence of the confining pressure on peak strength. In the future, the model could be adapted to other types of bonded materials like asphalt or rock.     

基于临界状态模型的砂土非共轴本构模拟

传统的砂土本构理论隐含了应力和塑性应变率的共轴条件,无法客观描述主应力轴旋转试验中的非共轴现象,并且当密度和围压变化较大时也不适用。基于材料状态相关砂土临界状态概念,将Pietruszczak和Stolle所提出的砂土本构模型进行了改进,并在模型中引入非共轴塑性流动理论来描述非共轴现象。通过对单剪试验和空心圆柱试验进行数值模拟,表明基于临界状态理论的非共轴模型能够合理描述主应力轴旋转过程中砂土的非共轴变形特性     

饱和土体静态液化失稳理论预测

通过在屈服准则和剪胀性方程中引入材料状态参数,建立了一个与材料状态相关的本构模型,用于模拟不同初始孔隙比和围压下土体的应力-应变关系。基于二阶功理论,建立了饱和土体静态液化失稳理论判别准则。通过预测一系列饱和松砂三轴不排水试验结果,验证了所建立的本构模型及判别预测准则的合理性。分析结果表明,静态液化发生于土体硬化阶段,静态液化触发伴随着剪应力达到峰值,其后剪应力降低且孔隙水压力持续增长。进一步地,模拟了充分降解的城市固体废弃物在不排水条件下的应力-应变特性,并预测到了潜在失稳线及静态液化触发点。     

土体应力-应变关系转型问题分析

通过对土体应力应变试验结果的分析,明确指出了土体普遍存在的应力-应变关系转型现象。提出了按后破坏模量和临界围压两个参数判别土体应力-应变关系转型的方法,对影响土体应力-应变关系转型的因素包括土的类型、围压、初始密度或含水率、排水条件及应力路径等进行了简要分析。最后,对砂土初步提出了一个能反映转型现象的组合应力-应变关系模型的理论框架。     

基于Hoek-Brown准则的岩石应变软化模型研究

基于Hoek-Brown（H-B）准则，根据大量试验数据，研究提出一种具有广泛适用性的岩石应变软化模型。利用40种硬岩和软岩三轴试验数据，验证了3参数幂函数型围压函数广泛适用于硬岩和软岩峰后临界塑性应变与围压的关系描述，并将其引入塑性内变量中；其次，比较分析当前基于H-B准则的岩石应变软化模型中软化参数的合理性，明确以GSI作为软化参数的方法最优，拟合结果正确且参数最少。然后，分析岩石峰后GSI随塑性内变量的演化规律，并建立一种普适的非线性演化模型。与其他演化模型相比，所建模型可更准确地描述不同种类岩石峰后H-B强度参数的非线性演化规律，同时参数较少。接着，采用非关联流动法则，将所建GSI演化模型嵌入有限元软件ABAQUS中，实现软化模型的数值求解。最后，对多种岩石的室内试验模拟结果表明，所建模型可正确描述不同围压下硬岩和软岩的峰后应变软化变形及剪胀变形，具备良好的适用性。与定临界塑性应变或定剪胀角的应变软化模型模拟结果相比，所建模型模拟曲线可一致性地与试验曲线接近。     

基于新型损伤定义的岩石损伤统计本构模型探讨

针对几何损伤理论中损伤定义的局限性与不足,对岩石损伤定义进行了改进;在此基础上,将在应力作用下的岩石材料抽象为破坏与未破坏两部分,并根据这两部分不同的受力情况,建立出新型损伤模型;其次,利用统计损伤理论,建立出岩石统计损伤演化方程,进而建立起基于特定围压下的岩石损伤统计本构模型;然后,通过探讨模型参数与围压的经验关系,对模型进行合理修正,从而建立了反映不同围压条件下的岩石损伤统计本构模型,并通过探讨模型参数的物理意义,初步建立了反映岩石软化与硬化特性相互转化的临界物理力学参数的确定方法。该模型能够反映岩石软化与硬化特性相互转化的特征,与试验结果进行比较分析表明,该模型与实测结果吻合良好。     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

Experimental and theoretical investigation of the high‐pressure,undrained response of a cohesionless sand

Key in predicting stability in sands during dynamic events is gaining a fundamental understanding of the physics of its deformation and failure at high pressures. In this paper are reported results of an experimental investigation into the high‐pressure (up to 700 MPa) mechanical response of Quikrete sand. During all triaxial compression tests, the material exhibited hardening up to failure while both compressibility and dilatancy regimes of the volumetric response were observed. Furthermore, the transition from compressibility to dilatancy was found to be highly dependent on confining pressure. By performing triaxial compression tests with several creep stages followed by unloading–reloading cycles, the time influence of the overall response was detected. Using the experimental data, a new model that captures both compressibility and dilatancy has been developed. Comparison between model predictions and data showed that the proposed model describes well the main characteristics of the high‐pressure response of sand. Copyright © 2012 John Wiley & Sons, Ltd.     

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

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