岩石内变量蠕变模型研究

通过一系列砂岩、页岩和辉绿岩室内蠕变试验,获得了不同应力和初始应变状态下岩石蠕变特性。通过误差分析讨论了计算蠕变速率时如何选择时间增量。讨论了复杂加、卸载条件下时间相关蠕变模型的局限性,提出采用不可恢复应变作为内变量来描述复杂条件下岩石蠕变性质,并据此提出了岩石的内变量蠕变模型。试验数据与模型拟合结果对比表明,内变量蠕变模型较好地反映了3种岩石的蠕变行为。此外,还初步讨论了模型在岩石率相关性质研究中的适用性,并利用相关试验数据进行了验证。     

Numerical modelling for clays with internal discontinuity surfaces

A numerical model capable of performing deformation analysis of a medium containing discontinuity surfaces is presented. The discontinuity can be either a crack, which can be open or closed, or a shear band. The model consists of two separate numerical algorithms, which are coupled together by means of the principle of superposition. In particular, an integral equation scheme based on the theory of dislocations is adopted for modelling the discontinuity, while a finite element discretization is used for the continuous medium. In this paper the discontinuity modelling is illustrated in detail, together with the specific formulation of the principle of superposition adopted, and some simple examples of application are presented. The well-known modelling approach based on Fracture Mechanics theory is also briefly discussed. The two models are compared and some advantages and drawbacks of each are pointed out, comments are made regarding their applicability in the specific case of soil mechanics, and conclusions are drawn as regards the conditions under which one or the other is appropriate. Finally, a full-scale example of deformation analysis using the proposed model is presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Disturbed state model for sand–geosynthetic interfaces and application to pull-out tests

Successful numerical simulation of geosynthetic-reinforced earth structures depends on selecting proper constitutive models for soils, geosynthetics and soil–geosynthetic interfaces. Many constitutive models are available for modelling soils and geosynthetics. However, constitutive models for soil–geosynthetic interfaces which can capture most of the important characteristics of interface response are not readily available. In this paper, an elasto-plastic constitutive model based on the disturbed state concept (DSC) for geosynthetic–soil interfaces has been presented. The proposed model is capable of capturing most of the important characteristics of interface response, such as dilation, hardening and softening. The behaviour of interfaces under the direct shear test has been predicted by the model. The present model has been implemented in the finite element procedure in association with the thin-layer element. Five pull-out tests with two different geogrids have been simulated numerically using FEM. For the calibration of the constitutive models used in FEM, the standard laboratory tests used are: (1) triaxial tests for the sand, (2) direct shear tests for the interfaces and (3) axial tension tests for the geogrids. The results of the finite element simulations of pull-out tests agree well with the test data. The proposed model can be used for the stress-deformation study of geosynthetic-reinforced embankments through numerical simulation. Copyright © 1999 John Wiley & Sons, Ltd.     

非饱和黄土结构性定量试验研究

非饱和黄土是典型的结构性土,对4种含水率的原状非饱和黄土和重塑非饱和黄土进行了3种围压下的三轴剪切试验。基于二元介质理论,通过原状非饱和黄土和重塑非饱和黄土的应力-应变曲线对比,定义了结构性参数结构应力分担比。考察了不同围压和含水率对结构应力分担比随应变变化的影响,并给出了其定量表达式。试验结果表明,围压和含水率是影响结构性的两个主要因素;随着应变增长,结构应力分担比分先线性增长后指数形式减低,其最大值可以表示成围压和含水率的函数。最后给出了考虑结构性的原状非饱和黄土本构模型的建模思路。     

Degradation of rock fracture asperities in unloading,reloading, and reversal

Among the constitutive models for rock fractures developed over the years, Barton's empirical model has been widely used. Although Barton's failure criterion predicts peak shear strength of rock fractures with acceptable precision, it has some limitations in estimating the peak shear displacement, post‐peak shear strength, dilation, and surface degradation. The first author modified Barton's original model in order to address these limitations. Barton proposed his model for degradation of fracture asperities in unloading, reloading, and shear displacement reversal based on just one cyclic direct shear test. In this study, a database of results of 18 cyclic direct shear tests available in the literature was collected and analyzed. Modifications were made to Barton's original model (in terms of fracture cyclic shearing) to make it consistent with the modified model proposed by the first author. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Effect of suction on the mechanical behaviour of iron ore rock

The effect of suction on the behaviour of iron ore has been studied from both physical and mechanical points of view. The porosity and the suction phenomena have been analysed using different experimental techniques. Uniaxial compressive tests on partially saturated samples have shown that the suction is responsible for strength and cohesion improvement. Considering the theory of partially saturated porous soils of Coussy and Dangla (Mécanique des sols non saturés (2002 edn). Hermès Science: 2002; 390), we have proposed a constitutive law for partially saturated iron ore. The real increase in the apparent cohesion due to the capillary attraction forces is overestimated if the yield function is written in terms of effective stresses. The effect of the capillary cohesion has been modelled with a function in the expression of the apparent cohesion of the yield function. The effect of suction on the mechanical behaviour has been represented in the effective stresses space and in the total stresses space like the Alonso model (Géotechnique 1990; 40 :405–430). Copyright © 2005 John Wiley & Sons, Ltd.     

基于Burgers模型考虑损伤的非定常岩石蠕变模型

Burgers模型只能描述岩石蠕变的前两个阶段,为了描述蠕变全过程,考虑蠕变参数的时间效应及损伤带来的影响,采用非定常黏性元件取代Burgers模型中串联的定常黏性元件,使其能描述加速蠕变阶段。根据损伤变量的变化特征,假定了一个函数,经过Lemaitre应变等效原理,代入Kelvin模型可得到能描述蠕变衰减阶段的模型。试验数据拟合结果显示,改进的Burgers模型能很好地描述加速蠕变阶段,模型拟合曲线与试验曲线基本吻合,相关系数高,参数取值也在合理范围内;非定常Burgers模型更适合描述不同应力下的岩石蠕变特征曲线。     

Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model

A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd.     

New mineral activity–composition relations for thermodynamic calculations in metapelitic systems

New activity–composition (a–x) relations for minerals commonly occurring in metapelites are presented for use with the internally consistent thermodynamic dataset of Holland & Powell ( 2011 , Journal of Metamorphic Geology, 29 , 333–383). The a–x relations include a broader consideration of Fe₂O₃ in minerals, changes to the formalism of several phases and order–disorder in all ferromagnesian minerals where Fe–Mg mixing occurs on multiple sites. The a–x relations for chlorite, biotite, garnet, chloritoid, staurolite, cordierite, orthopyroxene, muscovite, paragonite and margarite have been substantially reparameterized using the approach outlined in the companion paper in this issue. For the first time, the entire set of a–x relations for the common ferromagnesian minerals in metapelitic rocks is parameterized simultaneously, with attention paid to ensuring that they can be used together to calculate phase diagrams of geologically appropriate topology. The a–x relations developed are for use in the Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCKFMASHTO) system for both subsolidus and suprasolidus conditions. Petrogenetic grids in KFMASH and KFMASHTO are similar in topology to those produced with earlier end‐member datasets and a–x relations, but with some notable differences. In particular, in subsolidus equilibria, the FeO/(FeO + MgO) of garnet is now greater than in coexisting staurolite, bringing a number of key staurolite‐bearing equilibria into better agreement with inferences from field and petrographic observations. Furthermore, the addition of Fe³⁺ and Ti to a number of silicate phases allows more plausible equilibria to be calculated in relevant systems. Pseudosections calculated with the new a–x relations are also topologically similar to equivalent diagrams using earlier a–x relations, although with many low variance fields shifting in P–T space to somewhat lower pressure conditions.