法向卸荷下贯通裂隙岩样压剪破坏特征研究

开展贯通裂隙岩样在法向卸荷下的压剪破坏研究,有助于更好地理解岩体裂隙在卸荷条件下的失稳特征。对单裂隙贯通岩样开展剪切应力峰前一系列法向卸荷破坏试验,分析了法向卸荷下试样裂隙的变形及能量演化规律。结果表明,法向卸荷下裂隙压剪失稳时的法向应力大于常规直剪试验的对应值,裂隙的抗剪能力减弱;剪切位移随法向应力减小而增大;卸荷变形比K随裂隙轮廓面积比Rs的增加而增大;卸荷过程中总变形能U0由减小转变为增加的现象可用于预测岩石裂隙卸荷失稳。研究结果对理解贯通裂隙岩体压剪卸荷失稳破坏具有一定参考意义。     

Analytical‐numerical solution for stress distribution around tunnel reinforced by radial fully grouted rockbolts

This paper presents an analytical‐numerical approach to obtain the distribution of stresses and deformations around a reinforced tunnel. The increase in the radial stress of the reinforced tunnel, based on the performance of a bolt, is modeled by a function, which its maximum value is in the vicinity of the bolt periphery and it exponentially decreases in the far distance from the bolt. On the basis of this approach, the shear stiffness between the bolt and the rock mass and the shear stress distribution around the bolt within the rock mass are also analytically obtained. The results are compared with those obtained by the assumption of ‘uniform increase of radial stress’ method, which is made by the previous studies. The analyses show when the bolts' spacing is large, the safety factor must be increased if the ‘uniform increase of radial stress’ method is used for the design. Finally, a procedure is introduced to calculate the non‐equal deformation of the rock mass between the bolts at any radius that can be useful to compute the bending moment in shotcrete layer in New Austrian Tunnelling Method (NATM) approach. Copyright © 2016 John Wiley & Sons, Ltd.     

非贯通类节理岩体剪切变形规律分析

通过直剪条件下的模型试验研究具有相同连通率、含不同起伏角的共面闭合节理在不同法向应力作用下的剪切变形特征。试验结果表明：在达到峰值剪切强度前,剪切应力-位移曲线具有明显的线性段与非线性段;岩桥在法向与剪切荷载的共同作用下会产生一定程度的弱化,积累到一定程度时宏观表现为剪切应力软化。详细阐述了试验过程中裂纹发展的4个阶段并描述了剪切应力-剪切变形曲线变化的特征,提出一个能反映剪切应力-剪切变形曲线全过程的经验本构模型。新模型采用分段函数描述曲线变化特征。直剪试验数据分析表明,该模型能够拟合剪切应力-剪切变形变化发展趋势且具有较高的拟合精度,验证了模型的正确性。最后,对其存在的不足亦进行了简要分析     

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

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

Gouge Particle Evolution in a Rock Fracture Undergoing Shear: a Microscopic DEM Study

The evolution of gouge materials in rock fractures or faults undergoing shear can change fracture properties in terms of shear strength and dilation, fluid transmissivity and retardation for contaminants. In order to conceptually understand gouge mechanical behaviors including movement, microcracking, abrasion and redistribution, particle mechanics models were used to simulate single- and multi-gouge particles in a rough fracture segment undergoing shear. The results show that gouge particles behave in two different ways under low and high normal stresses, respectively. Under low normal stress, gouge particles mainly roll with the moving fracture walls, with little surface damage and small dilation during the shear process. Under high normal stress, gouge particles can be crushed into a few major pieces and a large number of minor comminuted particles, accompanied by more severe damage (abrasion and microcracking) in fracture walls and continuous fracture closure. The modeling results were also compared with published experiments and used to explain the observed macroscopic behaviors of rock fracture undergoing shear. The effects of microparameters used in the particle mechanics models on the simulation of gouge behaviors were also investigated through sensitivity analysis.     

Contact mechanism of a rock fracture subjected to normal loading and its impact on fast closure behavior during initial stage of fluid flow experiment

Fast closure of rock fractures has been commonly observed in the initial stage of fluid flow experiments at environmental temperatures under low or moderate normal stresses. To fully understand the mechanisms that drive this fast closure, the evolution of local stresses acting on contacting asperities on the fracture surfaces prior to fluid flow tests needs to be evaluated. In this study, we modeled numerically the asperity deformation and failure processes during initial normal loading, by adopting both elastic and elastic–plastic deformation models for the asperities on a real rock fracture with measured surface topography data, and estimated their impact on initial conditions for fluid flow test performed under laboratory conditions. Compared with the previous models that simulate the normal contact of a fracture as the approach of two rigid surfaces without deformations, our models of deformable asperities yielded smaller contact areas and higher stresses on contacting asperities at a given normal stress or normal displacement. The results show that the calculated local stresses were concentrated on the contacts of a few major asperities, resulting in crushing of asperity tips. With these higher contact stresses, however, the predicted closure rates by pressure solution are still several orders of magnitude lower than that of the experimental measurements at the initial stage of fluid flow test. This indicates that single pressure solution may not likely to be the principal compaction mechanism for this fast closure, and that the damages on contacting asperities that occur during the initial normal loading stage may play an important role. Copyright © 2015 John Wiley & Sons, Ltd.     

非规则砂岩节理剪切变形本构关系试验研究

岩石节理剪切变形对岩体工程的安全性和稳定性具有重要影响。为研究常法向应力下岩石节理剪切变形本构关系,采用RDS-200型岩石直剪仪对非规则砂岩节理进行了不同法向应力下的直剪试验。根据岩石节理剪切应力?位移全程曲线形状特征,将其依次划分为峰前压密阶段、线性阶段、屈服阶段和峰后软化阶段;根据剪切应力在峰后软化阶段降低幅度和速率大小,将岩石节理剪切应力?位移曲线划分为3种类型：峰后平台型、峰后缓降型和峰后跌落型。基于岩石节理剪切应力?位移曲线各阶段剪切变形特征,采用分段函数建立了岩石节理剪切变形本构模型。与其他模型相比,新提出的岩石节理剪切变形本构模型对试验数据拟合精度更高,更好地描述了岩石节理剪切应力?位移全程曲线。另外,在通过岩石节理直剪试验由经验公式确定模型参数之后,所提出本构模型可在不同法向应力下实现对不同粗糙度岩石节理剪切应力?位移曲线的预测。研究结果对岩石节理剪切变形的数值模拟和工程估计具有一定的实用价值。     

岩石多功能剪切试验测试系统研制

详细介绍了所研制的岩石多功能剪切试验测试系统的主要功能、技术指标和仪器组成,并开展了一系列相关力学试验。该试验测试系统主要包括试验装置、测量系统和控制系统3部分,最大法向拉伸应力为40 MPa,最大法向压缩应力为120 MPa,最大水平剪切应力为120 MPa,试样尺寸为50 mm?50 mm?50 mm;可开展多种力学试验,包括直接拉伸试验、直接剪切试验、拉伸-剪切试验和压缩–剪切试验。利用该试验测试系统对花岗岩进行试验研究,研究结果表明：直接拉伸试验中,试样发生脆性破坏,声发射信号瞬间达到峰值,破坏断面表现出拉伸破坏特征;直接剪切试验中,试样发生多次破坏,破坏瞬间声发射信号均发生突增,破坏断面表现出剪切破坏特征;拉伸-剪切试验中,试样在拉应力作用下剪切强度显著降低,声发射信号在破坏阶段表现强烈,破坏断面既有拉伸破坏特征也有一定的剪切破坏特征。上述力学试验结果,表明了所研制的岩石多功能剪切试验测试系统能够开展多种力学试验,为进一步研究岩石的剪切力学特性提供新的测试手段。     

岭澳核电三期强风化角岩边坡岩体直剪试验研究

采用堆载法进行天然状态与饱和状态下强风化角岩边坡岩体的现场直接剪切试验,获得不同正应力水平下剪切应力-剪切变形关系曲线和剪切强度参数,对其剪切应力-剪切变形关系曲线特征和不同正应力作用下剪应力随正应力的变化规律以及水-力耦合作用对剪切强度与变形特性影响进行了分析。试验结果表明：岩体峰值剪切强度和屈服剪切强度均随正应力的增大而增大;天然状态和饱水状态下剪应力随不同正应力的变化趋势基本相同;岩体剪切强度随含水率的增大而减小,在低法向应力下尤其敏感;水对岩体强度参数中黏聚力c的弱化作用更加明显,同时加大了岩体变形量,延长了岩体变形过程;通过现场直剪试验测得的法向变形可以估算岩体的压缩模量,为边坡稳定性分析提供参数。     

Stress‐dependent elastic wave velocity of microfractured sandstone

A model for the stress‐dependent elastic wave velocity response of fractured rock mass is proposed based on experimental evidence of stress‐dependent fracture normal and shear stiffness. Previously proposed models and previous experimental studies on stress‐dependent fracture stiffness have been reviewed to provide a basis for the new model. Most of the existing stress‐dependent elastic wave velocity models are empirical, with model parameters that do not have clear physical meanings. To propose the new model, the rock mass is assumed to have randomly oriented microscopic fractures. In addition, the characteristic length of microfractures is assumed to be sufficiently short compared to the rock mass dimensions. The macroscopic stress‐dependent elastic wave velocity response is assumed to be attributed to the stress dependency of fracture stiffness. The stress‐dependent fracture normal stiffness is defined as a generalized power law function of effective normal stress, which is a modification of the Goodman's model. On the other hand, the stress dependency of fracture shear stiffness is modeled as a linear function of normal stress based on experimental data. Ultrasonic wave velocity responses of a dry core sample of Berea sandstone were tested at effective stresses ranging from 2 to 55 MPa. Visual observation of thin sections obtained from the Berea sandstone confirms that the assumptions made for microstructure of rock mass model are appropriate. It is shown that the model can describe the stress‐dependent ultrasonic wave velocity responses of dry Berea sandstone with a set of reasonable material parameter values. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

The role of non‐coaxiality in the simulation of strain localization based on classical and Cosserat continua

It is normally accepted that materials inside the shear band undergo severe rotation of the principal stress direction, which causes non‐coaxiality between the principal stress and principal plastic strain rate. However, classical plasticity flow theory implicitly assumes that the principal stress and the principal plastic strain rate are coaxial; thus, it may not correctly predict the onset of the shear band. In addition, classical continuum does not contain any internal length scales; as a result, it cannot provide a reasonable shear band thickness. In this study, the original vertex non‐coaxial plastic model based on the classical continuum is extended to the Cosserat continuum. The corresponding codes are implemented via the interface of the user defined element subroutine in ABAQUS. Through a simple shear test, the effectiveness of the user's codes is verified. Through a uniaxial compression test, the influence of non‐coaxiality on the onset, the orientation, and the thickness of the shear band is investigated. Results show that the onset of the shear localization is delayed, and the thickness of the shear band is widened when the non‐coaxial degree increases, while the orientation of the shear band is little affected by the non‐coaxial degree. In addition, it is found that the non‐coaxiality can weaken the micro‐polar effect to some extent; nonetheless, the Cosserat non‐coaxial model still has its advantage over the classical non‐coaxial model in capturing the pre‐bifurcation as well as the post‐bifurcation behaviors of strain localization. Copyright © 2016 John Wiley & Sons, Ltd.     

节理分布位置对岩体剪切破裂特征影响试验研究

节理分布位置影响岩体综合抗剪强度,并控制“锁固段”型岩质边坡启动和破裂演化。通过模型试件剪切试验,分析了节理分布位置对岩体综合抗剪强度和破裂特征的影响程度。试件模型以水泥砂浆为材料,在剪切面不同位置预制封闭节理,使其成为四周密封、内部节理贯通的模型。在不同法向应力下进行剪切试验,测定综合抗剪强度变化规律和剪切破裂特征。试验结果表明:试样破裂面可分为翼部拉张区、翼部压剪区、端部拉张区三部分,各区面积所占比例随节理位置变化而波动。剪切强度参数也随节理位置变化而规律性波动,且法向应力越大,变化幅度越大。     

基于PIV技术的贯通单裂隙岩体压-剪破坏特征试验

基于PIV技术,采用相似材料预制贯通单裂隙岩体试件并开展压-剪试验,探究了裂隙开度对贯通单裂隙试件的强度、变形及破坏过程的影响。结果表明:(1)试件发生初始断裂时不同裂隙开度试验组应力均明显下降,且随着开度的增大裂隙岩体胶结部位塑性变形增强,岩体试件的特征应力值减小;(2)裂隙开度影响到裂隙和胶结部位变形在试件总变形中所占比例,进而影响裂隙试件的轴向峰值应变和横向峰值应变,开度越大裂隙试件的轴向峰值应变与横向峰值应变越小;(3)随着开度增大试件初始断裂部位由中下部向中上部转移,同时试件逐渐由拉伸剪切复合破坏转变为单一的剪切破坏;(4)需达到位移阈值k试件表面才能产生裂隙,而预制裂隙开度增大会降低试件的位移阈值k。其结果为研究深部矿井裂隙沟通造成的突水溃砂机理奠定了一定的基础。     

Macro‐and Microstructural,Textural Fabrics and Deformation Mechanism of Calcite Mylonites from Xar Moron‐Changchun Dextral Shear Zone,Northeast China

The calcite mylonites in the Xar Moron-Changchun shear zone show a significance dextral shearing characteristics. The asymmetric(σ-structure) calcite/quartz grains or aggregates, asymmetry of calcite c-axes fabric diagrams and the oblique foliation of recrystallized calcite grains correspond to a top-to-E shearing. Mineral deformation behaviors, twin morphology, C-axis EBSD fabrics, and quartz grain size-frequency diagrams demonstrate that the ductile shear zone was developed under conditions of greenschist facies, with the range of deformation temperatures from 200 to 300°C. These subgrains of host grains and surrounding recrystallized grains, strong undulose extinction, and slightly curved grain boundaries are probably results of intracrystalline deformation and dynamic recrystallization implying that the deformation took place within the dislocation-creep regime at shallow crustal levels. The calculated paleo-strain rates are between 10~(–7.87)s~(–1) and 10~(–11.49)s~(–1) with differential stresses of 32.63–63.94 MPa lying at the higher bound of typical strain rates in shear zones at crustal levels, and may indicate a relatively rapid deformation. The S-L-calcite tectonites have undergone a component of uplift which led to subhorizontal lifting in an already non-coaxial compressional deformation regime with a bulk pure shear-dominated general shear. This E-W large-scale dextral strike-slip movement is a consequence of the eastward extrusion of the Xing'an-Mongolian Orogenic Belt, and results from far-field forces associated with Late Triassic convergence domains after the final closure of the Paleo-Asian Ocean.     

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.     

Discrete element modelling of material non‐coaxiality in simple shear flows

We investigate the quasi‐static simple shear flow of a two‐dimensional assembly of cohesionless particles using discrete element method (DEM) simulations. We focus on the unsteady flow regime where the solid would experience significant evolution of stresses, mobilised shear strength and dilation. We construct the DEM model using a discretised‐wall confined granular cell where the apparent boundary is allowed to dilate or contract synchronously with the confined solid. A rather uniform simple shear field is achieved across the whole assembly, which benefits rheological studies in generalising constitutive laws for continuum methods. We examine two aspects of the simple shear behaviour: macroscopic stress and strain rate evolution, particularly the non‐coaxiality between the principal directions of the two; and micromechanics such as evolution of fabric. For an initially anisotropic specimen sheared under constant normal pressure, the direction of principal stress rotates towards that of the principal strain rate, gradually reducing the degree of non‐coaxiality from about 45° to fluctuating around 0°. The rate in approaching coaxiality is slower in samples with larger initial porosity, stress ratio and mean stress. Generally, a faster rate in approaching coaxiality in simple shear is observed in a more dilatant sample, which often shows a larger degree of mobilised fabric anisotropy, suggesting the possible important role of instantaneous internal friction angle. The evolution of principal fabric direction resembles that of the principal stress direction. © 2013 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons, Ltd.     

Effects of Scale and Normal Stress on Shear Behavior of a Soft Fracture Replica

The effects of the fracture scale (size) and normal stress on the shear behavior of a fracture were investigated in direct shear rig using mortar fracture replicas of different sizes (100–200 mm), which were reproduced from a tensile fracture created in granite. Monotonic shear loading was applied at normal stresses of 0.3 and 1 MPa for shear displacements of 2 and 20 mm. The change in surface topography was measured by using a profile measurement system with a laser profilometer to determine the changes in surface damage and the aperture during shear. Additionally, the effect of gouge material on the non-linearity of the closure curve during shear was investigated. The results show that the standard deviation of the initial aperture during shear increases not only with scale but also with the normal stress, since the matedness decreases due to the fact that asperities are more damaged with normal stress, and accordingly, the non-linearity in the closure curve during shear increases with normal stress. The results also show that the non-linearity of the closure curve during shear increases with the presence of gouge material/debris that is produced in the shearing process. This is mainly because gouge is soft and easily deforms under the application of normal stress. Moreover, gouge material is distributed non-uniformly around the contact areas. Thus, the normal stiffness during shear decreases not only with a decrease in the matedness between the surfaces, but also with the presence of gouge material.     

Non‐coaxial elasto‐plasticity model and bifurcation prediction of shear banding in sands

Numerous constitutive models built on coaxial theory and validated under axi‐symmetric condition often describe the stress–stain relationships and predict the inceptions of shear banding in sands inaccurately under true triaxial condition. By adopting an elaborated Mohr–Coulomb yield function and using non‐coaxial non‐associated flow rule, a 3D non‐coaxial elasto‐plasticity model is proposed and validated by a series of true triaxial tests on loose sands. The bifurcation analysis of true triaxial tests on dense sands predicts the influence of the intermediate principal stress ratio on the onset of shear band accurately. The failure of soils is shown to be related to the formation of shear band under most intermediate principal stress ratio conditions except for those which are close to the axi‐symmetric compression condition. Copyright © 2009 John Wiley & Sons, Ltd.     

Effects of boundary conditions and partial drainage on cyclic simple shear test results—a numerical study

Owing to imperfect boundary conditions in laboratory soil tests and the possibility of water diffusion inside the soil specimen in undrained tests, the assumption of uniform stress/strain over the sample is not valid. This study presents a qualitative assessment of the effects of non‐uniformities in stresses and strains, as well as effects of water diffusion within the soil sample on the global results of undrained cyclic simple shear tests. The possible implications of those phenomena on the results of liquefaction strength assessment are also discussed. A state‐of‐the‐art finite element code for transient analysis of multi‐phase systems is used to compare results of the so‐called ‘element tests’ (numerical constitutive experiments assuming uniform stress/strain/pore pressure distribution throughout the sample) with results of actual simulations of undrained cyclic simple shear tests using a finite element mesh and realistic boundary conditions. The finite element simulations are performed under various conditions, covering the entire range of practical situations: (1) perfectly drained soil specimen with constant volume, (2) perfectly undrained specimen, and (3) undrained test with possibility of water diffusion within the sample. The results presented here are restricted to strain‐driven tests performed for a loose uniform fine sand with relative density D_r=40%. Effects of system compliance in undrained laboratory simple shear tests are not investigated here. Copyright © 2004 John Wiley & Sons, Ltd.     

Fracture in a heterogeneous medium: a network approach

Fracture in a heterogeneous medium is simulated by a network model. Heterogeneity is introduced by randomly distributing an initial population of cracks on a triangular network of bonds. Each ruptured bond represents a crack. External stresses are applied onto the network and the evolution of the crack population is analysed when the stresses are increased. A crack will propagate to an intact bond if the critical crack extension force is reached. This propagation leads to the coalescence of bonds: crack clusters are formed. When their shapes are complex we approximate them by simpler ones to compute a crack extension force. Crack interactions are introduced by adding for a given bond the contributions of all adjacent dusters. By using this kind of approach combining fracture mechanics and network modelling, we are able to simulate the rupture of a rock specimen under compression/ traction or shear stresses and to give some characteristics of the macroscopic fracture.