干湿-冻融循环下延吉膨胀岩的力学特性及其应力-应变归一化

为了研究干湿、冻融和干湿-冻融循环作用对延吉膨胀岩的应力-应变关系、体变特性、抗剪强度的影响,开展了一系列的固结排水剪切试验。结果表明：未循环试样的应力-应变关系曲线表现为稳定型或弱应变软化型,3种环境作用后试样应力-应变关系曲线呈现出一定的应变软化特性,剪胀明显,而且围压越小、循环次数越高的试样软化和剪胀的程度越高;随着循环次数的增加,3种环境作用下试样的黏聚力急剧减小,而内摩擦角略有增加;干湿-冻融循环作用下试样的应力-应变关系软化程度最高,剪胀最显著,试样的黏聚力降低的幅度最大。基于Konder双曲线模型,对循环处理前后试样的应力-应变特性进行归一化分析,建立了考虑干湿-冻融循环次数和围压影响,能同时描述应变软化型和硬化型的应力-应变关系式,并对不同干湿-冻融循环次数试样的应力-应变曲线进行预测,预测值与实测值较为接近,预测效果较好。     

Identification of non‐linear stress–strain–time relationship of soils using genetic algorithm

Recognition of non‐linear constitutive rock/soil model from experimental results is often multi‐modal in the large parameter space. A genetic evolution algorithm is thus proposed for its recognition, including that of structure of the model and coefficients in the model. The structure of the model can be firstly determined according to mechanical mechanism if the mechanism is clearly understood or searched by using evolutionary algorithm. The coefficients to be determined are then searched in global optional space. With the new evolutionary algorithm, the non‐linear stress–strain–time constitutive law to describe strain softening behaviours of diatomaceous soil under consolidated and undrained state was recognized by learning stress–strain–time behaviour of an intact sample under consolidated pressure of σ_c=0.1 MPa and strain velocity of_a=0.175%/min. This model gave reasonable prediction for diatomaceous soils under varying consolidated pressures (0.1–3.5 MPa) and strain velocities (0.0044–1.75%/min). It indicates that the methodology proposed in this paper is robust enough and strongly attractive for recognition of non‐linear constitutive model of soil and rock materials. Copyright © 2002 John Wiley & Sons, Ltd.     

中主应力对硬岩破裂机制的影响

采用数值方法模拟硬岩的三轴压缩试验,应用4个本构模型,即Mohr-Coulomb模型、Drucker-Prager模型、应变软化模型及考虑变形模量劣化的应变软化模型,研究中主应力对均质及非均质硬岩破裂机制的影响。结果表明：不论是均质还是非均质岩体,中主应力对采用Drucker-Prager模型的岩体强度影响较大,而对采用其他本构模型的岩体强度影响不大;当非均质岩体采用Mohr-Coulomb模型或Drucker-Prager模型时,中主应力对岩体破坏过程影响不大,但对其破坏模式有较大影响;当非均质岩体采用应变软化模型时,中主应力对岩体破坏过程及其模式均有较大影响。针对工程算例,采用不同的本构模型获得的均质或非均质岩体强度相差很大,除Drucker-Prager模型外,同一本构模型的均质岩体强度远大于非均质岩体。实际工程岩体对中主应力的响应是不同的,故在地下结构设计过程中应根据岩体特性选择合理的本构模型以保证工程安全。     

Nonlinear Evolutionary Mechanisms of Instability of Plane-Shear Slope: Catastrophe, Bifurcation, Chaos and Physical Prediction

Summary. A cusp catastrophe model is presented and the necessary and sufficient conditions leading to landslides are discussed. The sliding surface is assumed to be planar and is a combination of two media: medium 1 is elastic-brittle or strain-hardening and medium 2 is strain-softening. The shear stress-strain constitutive model for the strain-softening medium is described by the Weibull’s distribution law. This paper is a generalization and extension of the paper by Qin et al. (2001b), in which the shear stress-strain constitutive model for medium 2 was described by a negative exponent distribution; a special case of the Weibull’s distribution law. It is found that the instability of the slope relies mainly on both the stiffness ratio of the media and the homogeneity index and that a new role of water is to enlarge the material homogeneity or brittleness and hence to reduce the stiffness ratio. A nonlinear dynamic model (also called a physical forecasting model), which is derived by considering the time-dependent behavior of the strain-softening medium, is used to study the time prediction of landslides. An algorithm of inversion on the nonlinear dynamic model is suggested for seeking the precursory abnormality and abstracting mechanical parameters from the observed series of a landslide. A case study of the Jimingsi landslide is analysed and its nonlinear dynamic model is established from the observation series of this landslide using the suggested model and the algorithm of inversion. It is found that the catastrophic characteristic index |D| shows a quick rise till reaching an extremely high peak value after the slope evolves into the tertiary creep, and subsequently approaches a zero value prior to instability, which can be regarded as an important precursory abnormality index. By taking into account the evolutionary characteristic of the slope being in the secondary creep, a simplified nonlinear dynamic model is proposed for studying the properties of bifurcation and chaos. It is shown that the emergence of chaos depends on the mechanical parameters of the sliding-surface media.     

Arbitrary Lagrangian Eulerian based finite element analysis of cone penetration in soft clay

This paper presents a numerical model for the analysis of cone penetration in soft clay based on the finite element method. The constitutive behaviour of the soil is modelled by modifying an elastic, perfectly-plastic soil model obeying Von-Mises yield criterion to take into account the strain-softening, rate dependent behaviour of soft clay. Since this is a problem involving large soil deformations, the analysis is carried out using an Arbitrary Lagrangian Eulerian method where the quality of the mesh is preserved during penetration. The variation of cone resistance is examined with various parameters such as rigidity index of the soil, in situ stress anisotropy and roughness at the cone–soil interface, which influence the penetration resistance of the cone. A theoretical correlation has been developed incorporating these parameters and the results have been compared with previous correlations based on the cavity expansion theory, finite element method and strain path method. With the increase in strain-softening, relative brittleness of the soil increases and the penetration resistance is significantly reduced. With the rising strain-rate dependency, penetration resistance increases but this increase is independent of the degree of brittleness of the soil.     

单轴拉伸岩样破坏过程及尺寸效应数值模拟

由于从实验及理论角度研究岩样单轴拉伸条件下的破坏全过程及尺寸效应难度都很大。因此采用拉格朗日元法来研究这些问题。在峰值强度之前后,岩石材料的本构模型分别取为线弹性及拉破坏线性应变软化模型。为了使拉伸塑性区不出现在试样的端部,在试样的两侧面中部预制了2个凹槽。数值模拟结果表明,全程拉应力-拉应变曲线分为峰前和峰后阶段。在接近峰值的峰前阶段,由于两凹槽附近具有明显的拉应力集中现象,拉伸塑性区最先出现在两凹槽附近。随着轴向拉应变的增加,发生拉伸破坏的单元的数目增加,新发生拉伸破坏的单元越来越接近试样的中心,直到两块拉伸塑性区在应变软化阶段贯通。两凹槽连线上各单元拉应力的分布呈现3个阶段,“澡盆型”（“U型”）阶段,“双峰型”（“M型”）阶段及“单峰型”（“Π型”）阶段。“澡盆型”阶段对应于全程拉应力-拉应变曲线的弹性阶段。“双峰型”阶段及“单峰型”阶段对应于全程拉应力-拉应变曲线的非弹性阶段（包括峰值强度之前的一小段,即应变硬化阶段及峰后的应变软化阶段）。增加试样高度及降低试样宽度,拉应力-拉应变曲线的软化段变得越来越陡峭,因而试样越容易发生失稳破坏。由于试样宽度较大时,试样内部的单元并非处于单向拉应力状态,因此,增加试样宽度,全程拉应力-拉应变曲线的峰值强度增加。当试样宽度较小时,从出现塑性区,到塑性区贯通所需要的时间步较小,或应变范围较窄。这说明试样的脆性较强,前兆不明显。前兆不明显的脆性破坏对应常见的是洞室岩爆、冲击地压及地震等灾害。     

Assessment of using borehole-closure data to determine the constitutive behavior of salt

Borehole closure data have often been used to determine the constitutive behaviour of rock salt. This is an inverse problem and analysis of the data is confused by the fact that borehole closure is a structural response that does not uniquely determine constitutive behaviour since both the stress and deformation depend on the constitutive behaviour. This application brief assesses the conclusion of a recent study that the transient power law, ? = K τ^mtⁿ, which includes no term to account for steady-state deformation, is a better constitutive model for salt than are models that account for steady-state straining, and that an empirical model developed from the borehole closure data could be extrapolated to give conservative predictions of long-term closure of single or widely spaced openings.     

应变软化岩体分析原理及其应用

应变软化是指应力-应变曲线中轴向应力随应变的增加而减小的现象,许多种类的岩土介质在工程扰动的作用下呈现应变软化的行为。在分析应变软化问题时,其应力-应变关系式中的切线刚度矩阵是非正定的,由此导致计算求解的困难。将岩体应变软化过程简化为一系列脆塑性过程,于是应变软化问题的求解归结为一系列脆塑性过程的分析。基于经典弹塑性力学理论,提出了应变软化过程模拟方法及其相应的有限元求解过程,编制了计算程序,研究了应变软化本构模型中不同强度弱化速率对圆形洞室围岩塑形区分布的影响,进一步分析了应变软化模型对应的隧道径向变形沿洞轴方向的分布特征,并与已有监测数据得到的分布规律进行了对比。初步的研究结果表明,应变软化模型得到的计算结果是比较合理的     

辽西地区黄土的强度与本构特性

采用应变控制式三轴仪对取自辽宁西部地区阜新-朝阳段黄土的力学性质进行了试验,获得了黄土的应力-应变关系和抗剪强度的变化规律;采用人工制备不同含水率试样的方法,测试并分析了含水率对黄土本构和强度特征的影响。结果表明：(1)试样的应力-应变关系在不排水剪条件下呈强软化型,表现为脆性破坏;在固结不排水剪条件下随初始含水率的增加由弱软化型向硬化型转变;而在排水剪条件下低含水率时呈软化型,随含水率增加,呈弱硬化型和硬化型。(2)饱水后试样的抗剪强度有很显著的降低:c_UU为13.7～363.3 kPa时,φ_UU为18.6°~40.2°;而c_UU为0.2～29.6 kPa时,φ_CU则为4°～23.4°。(3)土的有效抗剪强度指标均随含水率增加呈非线性减小的关系,粘聚力与含水率呈半对数函数关系,内摩擦角则与含水率呈幂函数关系。给出了黄土本构的数学模型和黄土抗剪强度指标的取值方法。     

Comparative study on three-dimensional statistical damage constitutive modified model of rock based on power function and Weibull distribution

According to the characteristics of the microscopic unit strength of rock with random distribution, the power function distribution, and the Weibull distribution which is widely used in the past as the distribution function of the strength of the rock. Based on the theory of rock damage and statistical strength theory, two damage constitutive models established under three different confining pressures are modified. Then, the damage constitutive modified model of two kinds of distribution is verified and compared with the existing three axis test data. The results show that: (1) In the stage of elastic deformation of rock, the two theoretical damage constitutive model curves are in high agreement with the three axis test curve; (2) The rock at the plastic stage (hardening and softening), damage constitutive model established by Weibull probability distribution and experimental curve fits better than the damage constitutive model established by power function distribution. Especially in high confining pressure, damage constitutive model based on Weibull distribution can well describe the rock deformation from brittle to ductile transition process and power function is not; (3) In the residual strength in rock, damage variable D of damage constitutive model based on power function distribution appeared more than 1, which is deviation from the actual one and damage constitutive model based on Weibull distribution is not deviating. To summarize, using Weibull distribution statistical probability model to describe the microscopic unit strength of rock is more reasonable.     

煤岩两体模型变形破坏数值模拟

采用拉格朗日元法,在弹性岩石与弹性-应变软化煤体所构成的平面应变两体模型的上、下端面上不存在水平方向摩擦力条件下,模拟了模型的破坏过程、岩石高度对模型及煤体全程应力-应变曲线、煤体变形速率、煤体破坏模式及剪切应变增量分布的影响。结果表明,当模型的全程应力-应变曲线达到峰值时煤体内部的剪切带图案已经十分明显,在模型的应变硬化阶段,煤体中的应变局部化可视为模型失稳破坏的前兆,随岩石高度的增加,模型应力-应变曲线的软化段变得陡峭,这与单轴压缩条件下的解析解在定性上是一致的;煤体应力-应变曲线的软化段变得平缓,煤体消耗能量的能力增强;弹性阶段煤体的变形速率降低;煤体内部的剪切应变增量增加。煤体应力-应变曲线的软化段的斜率、弹性阶段煤体的变形速率、煤体内部的剪切应变增量及塑性耗散能都受岩石高度的影响,说明了岩石几何尺寸对煤体的影响(煤岩相互作用)是不容忽视的。     

An extended Cam Clay model for soft anisotropic rocks

The mechanical behaviour of soft sedimentary rocks in triaxial compression is modelled by means of a constitutive law which has been formulated by extending the original Cam Clay model. The geological material is assumed to be linearly elastic-plastic, strain-hardening, with fixed transverse isotropy. Theoretical predictions are compared to observed behaviour at different confining pressures of specimens of diatomite with various bedding inclinations. Although quantitative agreement is not always satisfactory, the trend of the tests is qualitatively well matched, what proves the paradigmatic value of this type of model not only for soils but even for soft rocks.     

A hierarchical approach for constitutive modelling of geologic materials

A hierarchical concept is proposed for the development of constitutive models to account for various factors that influence behaviour of (geologic) materials. It permits evolution of models of progressively higher grades from the basic model representing isotropic hardening with associative behaviour. Factors such as non-associativeness and induced anisotropy due to friction and cyclic loading, and softening are introduced as corrections or perturbations to the basic model. The influence of these factors is captured through non-associativeness manifested by deviation from normality of the plastic strain increments to the yield surface, F. Details of four models: isotropic hardening with associative behaviour, isotropic hardening with non-associative behavioural anisotropic hardening and strain-softening with a damage variable are presented. They are verified with respect to laboratory multiaxial test data under various paths of loading, unloading and reloading for typical soils, rock and concrete. The proposed concept is general, yet sufficiently simplified in terms of physical understanding, number of constants and their physical meanings, determination of the constants and implementation.     

Stability analysis of slopes in soils with strain-softening behaviour

This paper presents a numerical approach to analyse the stability of slopes in soils with strain-softening behaviour. In these materials, a progressive failure can occur owing to a reduction of strength with increasing strain. Such a phenomenon can be analysed using methods that are able to simulate the formation and development of the shear zones in which strain localises. From a computational point of view, this presents many difficulties because the numerical procedures currently used are often affected by a lack of convergence, and the solution may depend strongly on the mesh adopted. In order to overcome these numerical drawbacks, in the present study use is made of a non-local elasto-viscoplastic constitutive model within the framework of the finite element method. The Mohr–Coulomb yield function is adopted, and the strain-softening behaviour of the soil is simulated by reducing the strength parameters with the increasing deviatoric plastic strain. To assess the reliability of the proposed approach, some comparisons with the results obtained using other constitutive models for soils with strain-softening behaviour are presented. Finally, a slope subjected to a prescribed process of weathering is considered, and the effects of this process on the slope stability are discussed.     

Measurement of parameters impacting methane recovery from coal seams

Summary This paper describes the behaviour of coalbeds as gas reservoirs and discusses the results of a study carried out to establish the effect of release of methane on gas flow behaviour of coal. Experimental work consisted of microscopy, establishing adsorption/desorption isotherms, and monitoring changes in the volume of coal matrix with increasing and decreasing gas pressure. Micrographs obtained using small pieces of coal indicated that coal is made up of blocks, containing matrix and pores, separated by microfractures. This confirms the dual porosity model of coal structure with a primary porosity, and a fracture/cleat porosity-physical model used in coalbed methane simulators developed recently. Isotherms suggested that for the samples tested, a major part of the gas is released only after pressure falls below 600 psi, and this is primarily due to desorbing gas. Results of the volumetric strain experiments indicated that there is an increase in matrix volume with increase in gas pressure, in spite of matrix compressibility. Adsorption, therefore, induces swelling of the matrix. With decrease in gas pressure from 1000 psi to atmospheric, the matrix volume shrunk by 0.5%. These experimental results were inputted in a reservoir model and simulation runs made to determine the effect of pore volume and matrix shrinkage compressibilities on gas production. Over a five year period 60% more gas was produced when matrix shrinkage was used as an input parameter.Editor's note: The units used in this paper are generally used by the Gas Research Institute and are found in most oil and gas publications. Conversions of the more important units are: 1 MMCFD 28 300 m³/day; 1 MSCFD 28.3 m³/day; 100 psi 0.68 MPa; 100 ft²/lb 46.87 m²/kg; 1 ft 0.3048 m; 1 acre 0.40 hectare.     

Testing and evaluation of strength and deformation behaviour of jointed rocks

The objective of the paper is to derive the strength and modulus properties of rockmass as a function of intact rock strength and joint factor. The joint factor reflects the combined effect of joint frequency, joint inclination and joint strength. A study for the strength and deformation characteristics of jointed rock is done by conducting standard laboratory tests on cylindrical specimens of plaster of Paris after introducing artificial joints. The specimens having one to four joints at different inclinations which vary from 0° to 90° were tested at different confining conditions. The test results were examined to understand the effect of joint frequency and joint inclination on the strength and deformation behaviour of rock mass. Empirical correlations were developed for prediction of the uniaxial compressive strength and elastic modulus of jointed rocks. Results are compared with the earlier work on jointed specimens covering a wide variety of rocks. So, knowing the intact rock properties and the joint factor, the jointed rock properties can be estimated. These relations can be used for developing an equivalent continuum model for rock mass for handling boundary value problems. A failure criterion as proposed by Ramamurthy (1993 Ramamurthy, T. 1993. “Strength and modulus response of anisotropic rocks”. In Comprehensive rock engineering, Edited by: Hudson, J.A. Vol. 1, 313–329. Oxford: Pergamon Press.  [Google Scholar]) has been validated from these experimental results.     

应变软化模拟与圆形隧道衬砌分析

在岩土工程中,有许多类岩体呈现应变软化特性。在应变软化过程中,如果多个材料强度参数分析发生变化,就很难给出软化问题的解析解。已有的文献研究表明,在经典的弹塑性理论框架下,模拟岩土介质的应变软化过程依然值得深入研究。利用应变软化材料强度参数随塑性应变增加而弱化的特点,提出了将应变软化过程简化为一系列的脆性应力跌落与塑性流动过程的思想,应变软化行为的模拟转化为一系列脆塑性过程的分析。将这种模拟应变软化过程的思想方法用于求解均匀初始应力场中应变软化介质内的圆形隧道开挖问题,研究了衬砌对减少岩土开挖扰动的作用以及考虑应变软化对衬砌内力的影响。数值结果表明,对于应变软化岩土材料,忽略其强度软化的行为是偏于危险的。     

Constitutive model for tailing soils subjected to freeze–thaw cycles based on meso-mechanics and homogenization theory

A constitutive model is proposed for tailing soils subjected to freeze–thaw cycles based on the meso-mechanics and homogenization theory. The evolution of meso-structure upon loading was analyzed within the framework of breakage mechanism. When the new model is formed, tailing soils are idealized as composite materials composed of bonded elements described by an elastic brittle model and frictional elements described by a double hardening model. Based on meso-mechanics and homogenization theory, the nonuniform distributions of stress and strain within the representative volume element are given by introducing a structure parameter of breakage ratio with the derivation of the strain coefficient tensor, which connects the strains of the bonded elements and the representative volume element. The methods for determining model parameters are suggested based on the available tested results. The model proposed here can predict the deformation properties of tailing soils experiencing freeze–thaw cycles with acceptable accuracy. The strain-hardening and post-peak strain-softening behaviors of tailing soils under various confining pressures as well as different numbers of freeze–thaw cycles are well captured, and the dilatancy and contraction features are also adequately represented.

     

Deformation and stability of an elasto-plastic softening pillar

Summary A model of rock pillar or coal seam is considered assuming linear elastic behaviour before reaching the maximum strength and post-peak behaviour characterized by the residual strength. The deformation and stress across the pillar height are assumed to be uniform and the interaction with overlying rock strata is treated assuming beam model of the strata. The elasto-plastic stress distribution within pillar and the onset of instability occurring for the critical opening span are determined. Comparison with a solution for a simplified spring model of pillar is also presented.