基于离散元法的节理岩体边坡稳定性分析

节理岩体边坡的稳定性在很大程度上取决于节理的强度及其分布形式。由于节理岩体边坡的失稳破坏具有大变形和非连续的特点,因此,离散单元法成为研究节理岩体边坡破坏机制的最有效方法之一。通过采用离散元软件PFC2D进行数值模拟,对完整岩石及节理的力学性能进行研究,并建立含密集节理的岩体边坡模型,讨论了节理连通率对边坡破坏形式的影响。结果表明,节理岩体边坡的失稳破坏是一个渐进的过程;在多组节理密集分布的岩体边坡中,连通率越大,其稳定性越差;随着连通率的减小,边坡的破坏形式由大范围的滑坡转变为局部崩塌的形式     

某水电站地下厂房洞室群围岩稳定分析

采用有限元方法反演某水电站地下厂房厂区初始地应力场。采用Zienkiewicz-Pande的双曲线屈服准则,对某水电站地下厂房洞室群的开挖过程进行了三维弹塑性有限元模拟。为反映锚固支护对围岩的加固作用和评价地下厂房洞室群支护参数的合理性,提出了可以考虑抗剪作用的隐式锚杆单元和可以考虑预应力效果隐式锚索单元。计算结果表明,该洞室群在洞室布局、开挖顺序和支护参数等方面是合理的,各主要洞室群的稳定性总体较好。     

某水电站地下厂房围岩稳定性分析

对某水电站地下厂房采用有限元方法反演厂区初始地应力场。分别对地下厂房洞室群3个平行洞室的不同间距、洞室开挖的不同顺序以及施加锚杆和锚索后的方案与毛洞方案进行了分析比较。计算结果表明,经过适当优化后该洞室群在洞室布局、开挖顺序和支护参数等方面是安全可行的,各主要洞室群的稳定性总体较好。     

基于离散元的强度折减法分析岩质边坡稳定性

将通用离散元UDEC与强度折减法结合,对含多结构面的岩质边坡的稳定性进行了分析。通过对节理岩质边坡的UDEC模型中的可变形块体和节理单元的强度参数进行折减,使模型不能再达到平衡状态,此时的折减系数就是边坡的安全系数,另外,由对应的边坡块体的速度矢量可以确定滑动面和边坡的破坏形态。通过与传统的条分法的结果比较,表明基于UDEC的强度折减法是一种可靠、有效的方法,为复杂节理岩质边坡的滑动面确定与安全系数计算开辟了新的途径。     

Geoengineering techniques used in the construction of underground openings in jointed rocks

Aoki, K. and Shiogama, Y., 1993. Geoengineering techniques used in the construction of underground openings in jointed rocks. In: M. Langer, K. Hoshino and K. Aoki (Editors), Engineering Geology in the Utilization of Underground Space.Eng. Geol., 35: 167–173.

Examples of underground openings used for pumped storage hydroelectric power stations and storage facilities for crude oil in Japan are described. In these examples the observational construction technique is important to compensate for missing data in the preliminary surveys. Some new drilling and assessment techniques that have been developed are also introduced.     

Continuum modelling of jointed porous rock

Constitutive equations for the mechanical and hydraulic behaviour of saturated porous rock with joint sets of specified orientations are developed by superposing continuum representations for the mechanical and hydraulic properties of the intact rock and each of the joint sets. The resulting continuum theory allows for fluid diffusion through and between interconnected rock pores and joint sets of specified orientation, and also accounts for the anisotropy of the mechanical properties due to joint stiffnesses. The accuracy and reliability of this model are verified by finite element simulation of example problems. The first example considers joint orientation-dependent rock deformation in a hypothetical porous medium with one joint set of different dip angles. More realistic examples related to rock slope stability and reservoir-induced seismicity are also considered in which the constitutive law's utility for modelling time-dependent fluid pressures is illustrated.     

Undrained shear testing of jointed rock

SummaryUndrained Shear Testing of Jointed Rock Water pressures must change inside joints undergoing shear without perfect drainage. This paper describes a new direct shear machine in which jacketed samples with oriented joints can be tested under consolidated undrained conditions with pore pressure measurement. Triaxial compression techniques for such tests are also described and typical results with intact and jointed sandstone samples are compared with results from the direct shear tests. Whereas intact specimens displayed increasing pore pressure followed by dilatancy and pore pressure decline before peak loading, the pore pressure in jointed specimens continued increasing right up to the peak load.

---

ZusammenfassungScherversuche an klüftigen Gesteinen ohne Dränung Wasserdruck in Klüften muß sich verändern, wenn ohne perfekte Dränage Scherkraft ausgeübt wird. In der vorliegenden Arbeit wird ein neuer, direkt wirkender Scherapparat beschrieben, in welchem ummantelte Prüfkörper mit gerichteten Klüften unter Konsolidationsbedingungen ohne Dränage untersucht werden können, wobei der Porendruck gemessen wird. Es wird außerdem die Technik von Triaxialversuchen beschrieben und typische Ergebnisse mit ungestörtem und geklüftetem Sandstein werden mit Ergebnissen aus direkten Scherversuchen verglichen. Während ungestörte Proben steigenden Porendruck, gefolgt von Dilatanz, aufweisen und der Porendruck vor der Höchstbelastung nachläßt, nimmt derselbe in geklüfteten Proben bis zur Höchstbelastung ständig zu.

---

RésuméEssais de cisaillement non drainé sur des échantillons de roche contenant une discontinuité Les pressions interstitielles à l'intérieur des discontinuités doivent changer lorsqu'elles sont soumises au cisaillement et que le drainage n'est pas parfait. Cet article décrit une nouvelle machine pour essais de cisaillement direct permettant de tester, dans des conditions de consolidation non drainée, des échantillons de roche gainés comportant une discontinuité d'orientation déterminée. La machine permet la mesure de la pression interstitielle.L'article décrit également les techniques de compression triaxiale correspondantes et des résultats typiques obtenus sur des échantillons de grès intact ou traversé par un joint sont comparés à ceux provenant des essais de cisaillement direct. Alors que les échantillons de roche intacte montrent un accroissement de la pression interstitielle suivi de dilatance et d'une réduction de pression avant que la contrainte de cisaillement maximum ne soit atteinte, la pression interstitielle dans les échantillons contenant une discontinuité croît jusqu'au droit de la contrainte de cisaillement maximum.

---

---

With 11 Figures     

Stability assessment of jointed rock slope with different crack infillings under various thermomechanical loadings

Thermoelastic deformation of rock significantly affects the stability of rock slope because thermoelastic strains may cause fracture propagation under favorable condition of failure. Rock slope stability depends on the balance between shear stress and shear resistance along the plane of weakness. Due to warming of rock slopes by heat transfer phenomena, viz. conduction and convection, considerable change in induced stresses (normal and shear) and resistance takes place which further causes instability in rock slope. In this paper, a two-dimensional finite element model has been used to simulate the stability of jointed rock slope containing crack in its upper surface. Four different cases have been simulated on the basis of infilling material (air, water, ice, water and ice) in the crack. Stability of rock slope is examined in terms of shear displacement and factor of safety for different thermal conditions of slope surface. A comparative study has been done for the four cases of infilling material in the crack. The various affecting parameters, viz. shear displacement, factor of safety, shear strength along the joint, and different surface temperature conditions, are illustrated by means of graphs. It has been found that the values of horizontal and vertical displacements are in the range of millimeters. The maximum values of horizontal and vertical displacements are 2.17 mm. Moreover, the maximum values of vertical compressive and tensile stresses are 15.4 MPa and 4.45 MPa respectively for the said four cases. According to the infilling material in the crack, the stability of the rock slope for the given geometry of slope is found in the following order: crack filled with ice < crack filled with ice and water < crack filled with water < empty crack. Validations of numerical results have been done from previous studies, and it has been found that the trends of normal stress, shear strength, and shear displacement along the joint are well matched.     

Blast damage control in jointed rock mass

Highly jointed rocks often cause problems associated with blast damage and the stability of the back and/or walls of the excavation. A field study was performed to understand the role played by the joint parameters in inducing blast damage. The field work included blasting of small scale models, drift rounds and monitoring of blast damage at several operating mines. The damage was assessed by blast vibration monitoring, half cast factor, overbreak measurement and visual inspection.

The effect of spacing, orientation, aperture, condition, filling material and wall strength of joints on blast damage is described. The interaction between the joint planes and explosive energy has been discussed and the overbreak control measures have been suggested.     

The transmission of compressional waves in jointed rock

Experiments have been carried out to determine the effect of joints on the propagation of P-waves in specimens of both weak and strong rocks; in particular the effect of normal stress on the transmission of P-waves in jointed rock was studied. These experiments were supported by measurements of seismic wave velocity carried out in tunnels driven in Lower Chalk and Carboniferous strata. The laboratory tests showed that relatively small normal stresses are sufficient to close joints acoustically and this was confirmed by the tunnel site observations. The results show that the interpretation of seismic data for the exploration of engineering sites may be improved by allowing for effects which result in the acoustic closure of joints.     

节理岩体蠕变特性研究

在已知岩石和节理蠕变规律的前提下,推导了节理岩体蠕变模型的一般表达式。假定岩石体积变形和节理法向压缩变形为弹性变形,忽略节理的剪胀现象,认为只有岩石畸变和节理剪切滑移与时间有关,从而推导了含三组相交节理的岩体蠕变模型及其参数。根据反演出的岩石和节理蠕变模型,计算了含三组相交节理的岩体在单轴应力作用下一些节理参数对岩体单轴蠕变的影响。分析表明,节理间距、剪胀系数以及节理夹角都对岩体的单轴蠕变变形有明显影响。节理间距越大,剪胀系数越大,节理夹角越小,节理岩体的单轴蠕变柔量也就越小,岩体的蠕变变形也越小。     

High lateral strain ratio in jointed rock masses

Rock masses encountered in the field are jointed and highly anisotropic. An important manifestation of this characteristic of rock masses is that the blocks, when slipping, separate and create voids in the mass resulting in high lateral strain even under low axial stress. It has been shown in this study that the ratio of lateral to axial strain may be very high, especially, if the joints are critically oriented. The assumption of isotropic linearly elastic material is not applicable in such situations. This observation is based on the outcome of an extensive laboratory testing programme, in which a large number of specimens of a jointed rock mass with various joint configurations were tested under uniaxial loading conditions. The trends of experimental results for both lateral strain ratio and rock mass strength have also been verified through distinct element modelling. The reason for high lateral strains has been attributed to the creation of voids and also to the fact that permanent deformations due to slip commence along rock joints right from the start of loading process. A simple mechanistic model has also been suggested to explain the high values of lateral strain for rough and dilatant rock joints. The lateral strains in such situations are important in the design of rock bolts. An example from a recently completed hydroelectric project has also been discussed where high lateral strain was found to be responsible for excessive wall closure and failure of rock bolts.     

Progressive failure in stratified and jointed rock mass

Summary This paper describes a mathematical model for analyzing progressive failure in a stratified and jointed rock mass dipping parallel to the slope face. The relationship between the necessary support forces and permissible displacement in the rock mass is discussed. An analytical expression is provided for the minimum required support force corresponding to a critical displacement. This expression is used to evaluate the potential of unstable failure propagation. The effect of time on progressive failure is considered by using a simple rheological model for joints. A case study, involving an excavation in a stratified medium, is used to compare model prediction with monitored performance.     

节理岩体爆破的DDA方法模拟

在非连续变形分析（DDA）方法中,通过跟踪炮孔扩张和炮孔周边裂隙的发展贯通求得爆腔的即时体积,进而根据爆生压力状态方程计算爆腔即时压力,并将爆生压力载荷作用到主炮孔内壁和贯通裂隙面上,实现了爆生产物作用下节理岩体爆破的DDA方法模拟。采用DDA方法模拟了节理岩体中的水平柱状炮孔抛掷爆破问题,得到了爆腔的体积扩张和压力衰减时间曲线,模拟很好的再现了岩石爆破过程中的炮孔扩张、岩体破坏、块体抛掷和爆堆形成过程。     

Stability analysis of a jointed beam

A numerical step-by-step procedure, analogous to the ‘Initial Stress Method,’ is presented for the analysis of a single-layer jointed rock beam subjected to gravity loads and in-plane in situ formation pressure. The joints are permitted to open at locations where the flexural stresses exceed flexural strength. The material properties may be different for each rock block and joint. A detailed algorithm is given for the solution of the problem. The results of several analyses indicating the relative effects of initial formation pressure, transverse load, stiffness of the joint material, and joint spacing on the response of a jointed beam are presented. The convergence characteristics of the numerical procedure are included. The joint material is assumed to be ‘no-tension’ type. Both the geometric and material non-linear effects are considered.     

Analysis of excavations in jointed rock of infinite extent

The subject of this work is the development of a code to study the behaviour of stratified and jointed rock masses around underground excavations. The rock mass is divided into two types of regions, one which is supposed to exhibit linear elastic behaviour and which may extend to infinity, and the other which will include discontinuities that behave inelastically. The former has been simulated by a symmetric direct, boundary integral, plane strain, quadratic, orthotropic module, and the latter by quadratic plane strain, membrane and inelastic joint elements. The two modules are coupled in one program. Sequences of loading include static point, pressure, body and residual loads, construction and excavation sequences, and quasi-static earthquake load. The programe is interactive, with graphics. A numerical example is presented to illustrate the method.     

一种新的节理裂隙岩体弹塑性模型

提出一种节理裂隙岩体弹塑性模型,将结构划分为块体单元和缝单元,其中缝单元可以是实际节理裂隙,也可以是人为缝单元。以块体单元形心的刚体位移和块体的平均应变作为基本未知量。该模型能充分考虑节理裂隙材料和块体材料的本构关系,计算量不大,是一种位移不协调单元。该模型特别适合节理裂隙岩体的数值分析。     

Microplane damage model for jointed rock masses

The paper presents a new microplane constitutive model for the inelastic behavior of jointed rock masses that takes into account the mechanical behavior and geometric characteristics of cracks and joints. The basic idea is that the microplane modeling of rock masses under general triaxial loading, including compression, requires the isotropic rock matrix and the joints to be considered as two distinct phases coupled in parallel. A joint continuity factor is defined as a microplane damage variable to represent the stress‐carrying area fraction of the joint phase. Based on the assumption of parallel coupling between the rock joint and the rock matrix, the overall mechanical behavior of the rock is characterized by microplane constitutive laws for the rock matrix and for the rock joints, along with an evolution law for the microplane joint continuity factor. The inelastic response of the rock matrix and the rock joints is controlled on the microplane level by the stress–strain boundaries. Based on the arguments enunciated in developing the new microplane model M7 for concrete, the previously used volumetric–deviatoric splits of the elastic strains and of the tensile boundary are avoided. The boundaries are tensile normal, compressive normal, and shear. The numerical simulations demonstrate satisfactory fits of published triaxial test data on sandstone and on jointed plaster mortar, including quintessential features such as the strain softening and dilatancy under low confining pressure, as well as the brittle–ductile transition under higher confining pressure, and the decrease of jointed rock strength and Young's modulus with an increasing dip angle of the joint. Copyright © 2014 John Wiley & Sons, Ltd.