基于离散裂隙网络模型的节理岩体渗透张量及特性分析

节理岩体几何结构非常复杂,研究其渗流特性对于指导含水岩层稳定性分析具有重要价值。应用离散裂隙网络模型DFN方法,基于VC++6.0软件平台,建立了平面渗流分析方法,分析了节理岩体不同几何分布情况下的渗透率张量特征,通过定义渗流定向性系数对岩体渗流的定向性特征进行了定量分析。结果表明：单组节理岩体渗流具有明显的各向异性特征,渗流定向性随着节理角度变化显著;节理随着节理贯通性增加,节理渗透率呈现对数增加趋势;两组节理情况下,各向异性特征随着节理组间夹角变化;两组节理岩体渗流特征研究中,正交分布下,岩体仍存在各向异性,但渗流定向性系数较低;当节理倾角服从正态分布时,随着节理倾角标准差增大,渗透率增加;两组节理夹角不同时,节理渗透主方向倾角随着夹角增大而相应增大,基本沿两组节理夹角方向的角平分线方向。     

Interpretation of torsional shear results for nonlinear stress–strain relationship

Torsional tests are frequently used in material testing to determine shear stress–strain relationships. In the tests, torque is applied to a solid or hollow cylindrical specimen, whose twist angle is recorded as the response. While the raw test results are readily available in terms of torque–twist sets, converting them into intrinsic stress–strain relationships is not always straightforward, owing to the nonlinear stress–strain behavior and the non‐uniform radial distribution of the stress and strain fields. It follows that the approaches currently used to interpret torsional data are essentially empirical based. The present paper addresses this data interpretation issue from a more rigorous perspective. The paper includes two parts: the first part concerns hollow cylindrical specimens under torsion, where an analytical‐based recursive procedure is developed; and the second part deals with solid torsional specimens, for which an analytical solution given by Nadai is directly employed. Examples pertaining to these two cases are examined and the resulting shear stress–strain curves are compared with those based on a number of popular empirical methods. The proposed analytical approaches can serve as a benchmark for checking the accuracy and validity of those empirical propositions. It is worth noting that, for brittle materials exhibiting softening behavior, only the analytical method gives the correct interpretation. Copyright © 2007 John Wiley & Sons, Ltd.     

Characterisation of confinement effect on jointed rock pillars using a Synthetic Rock Mass approach

Confinement effect on jointed rock pillars is numerically characterised in this research using a Synthetic Rock Mass (SRM) approach. The SRM is an integrated model incorporating a discrete fracture network within a Particle Flow Code 3D particle assembly. In this paper, the confinement effect on a 3D jointed pillar SRM model is investigated in a series of simulations, including biaxial compression tests and true and conventional triaxial compression tests. The numerical results suggest that the applied confining stresses generally result in higher pillar strengths and ductile post‐peak responses. More brittle post‐peak behaviour is simulated in the biaxial and true triaxial tests when the pillar is confined by a high stress in one lateral direction and by a zero/low stress in the other lateral direction. This phenomenon is attributed to significant lateral pillar dilation in the less confined direction. Detailed pillar failure modes are monitored in the uniaxial and triaxial tests. Axial splitting fractures and long shear zones cutting through the pillar are simulated when the pillar is able to dilate in the direction of least confinement. Localised shearing along joints and failed rock blocks is the dominant failure mode when the pillar dilation is resisted by the applied confining stresses. The pillar remains relatively intact with limited cracking in the pillar core in the highly confined triaxial tests. Copyright © 2016 John Wiley & Sons, Ltd.     

充水节理岩体中大断面隧洞开挖的数值模型研究

采用非饱和的渗流应力耦合模型分析了糯扎渡水电站2#导流洞的开挖过程,导流洞部分洞身穿过节理带,洞顶水头较高。详细地讨论了高水头和节理带对施工过程的影响。节理带采用描述高密度平行节理组的各向异性节理本构来逼近;考虑开挖引起的介质变形对渗透系数的影响;考虑排水引起的饱和度变化对渗透系数的影响。所用的非饱和瞬态耦合模型可以模拟出开挖引起的EDZ区域孔隙水压力急剧升高、有效应力减小、渗透系数动态的变化以及排水对洞室稳定性的提高。数值模拟的计算结果与国外类似试验的一般性观测结论相吻合,因此,可以用来评价水位以下隧洞施工方法和施工速度的合理性和经济性。     

Fluid permeability of sedimentary rocks in a complete stress–strain process

The permeability of sedimentary rocks during triaxial compression tests was investigated to relate it to the complete strain–stress process. It was found that the permeability was not constant, but varied with the stress and strain states in the rocks. Prior to the peak strength, the permeability decreases with increasing load. A dramatic increase in permeability occurs during the strain softening period. In the present study, in situ measurements of fluid flow and pressure in floor strata was carried out in a double longwall mining face in the Yangzhuang colliery. These measurements show that both the strata pressure and the position with respect to the mining face influence the hydrogeologic properties. The permeability increased in the floor strata behind the mining face because those mining induced fractures opened as the strata pressure decreased. To better understand this change in hydraulic behavior around the mining faces, 3-D numerical modeling was carried out. The model provides the general picture of the stress distribution and failure zone both in the floor and roof strata. The field and model results demonstrate the importance of changes in the stress and strain states on the hydrogeology of a site.     

Investigation on strength and stability of jointed rock mass using three‐dimensional numerical manifold method

This paper proposes a numerical model for jointed rock masses within the 3‐D numerical manifold method (NMM) framework equipped with a customized contact algorithm. The strength of rock sample containing a few sets of discontinuities is first investigated. The results of models with simple geometries are compared with the available analytical solutions to verify the developed computer code, whereas models with complex geometries are simulated to better understand the fundamental behavior and failure mechanism of jointed rock mass. Furthermore, the stability of jointed rock mass in an underground excavation is studied, where rock failure process is determined by the 3‐D NMM simulation. The simulation results provide valuable guidance on excavation process design and stabilization design in rock engineering practice. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

3D numerical modelling of bit–rock fracture mechanisms in percussive drilling with a multiple‐button bit

This paper considers numerical modelling of rock fracture induced by dynamic bit–rock interaction in percussive drilling. The work presented here extends the author's earlier research on the topic from the axisymmetric case to 3D case. The numerical method for modelling rock fracture includes a constitutive model for rock and a contact mechanics‐based technique to simulate the bit–rock interaction. The constitutive model is based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. This model is improved here from its earlier state by calibrating the softening laws using fracture energies G_Ic and G_IIc in tension and compression, respectively. Moreover, the viscosity modulus in tension is calibrated based on the dynamic Brazilian disc test. With these enhancements, the developed method is applied to 3D case of the bit–rock interaction problem assuming one symmetry plane. Single impact with single and multiple‐button bits is simulated. In the latter case, an initial borehole is modelled in order to simulate the usual in‐situ drilling conditions. The different failure types observed in the experiments as well as the interaction between the buttons resulting in chipping are realistically captured in the simulations. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical modelling of bit–rock fracture mechanisms in percussive drilling with a continuum approach

This paper presents a numerical method for continuum modelling of the dynamic bit–rock interaction process in percussive drilling. The method includes a constitutive model based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. The interaction between the drill bit and rock is modelled using contact mechanics by treating the bit as a rigid body. As the bit–rock interaction in percussive drilling is a transient event, the method is implemented in explicit dynamics FEM. The rock strength heterogeneity is characterized at the mesoscopic level statistically using the Weibull distribution. The bit–rock interaction is simulated under axisymmetric conditions using cylindrical and hemispherical buttons. The choice of the quite complex constitutive model accounting, e.g. for plastic compaction, viscoplastic shear and tensile failure along with induced damage and rate dependency is justified by numerical simulations. Moreover, the quasi‐static and dynamic cases are compared in plane strain simulations. Finally, some results clarifying the discrepancy of opinions found in the literature concerning the side (lateral) crack formation are obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

基于分形几何理论的DFN模型构建方法研究

以重庆梁-忠（梁平县-忠县）高速公路礼让隧道为工程研究背景,通过测线法调查现场节理,获得了节理产状分布概率密度函数,并从分形几何学的角度分析了节理间距及迹长的分形分布规律,推导出能反映节理间距及迹长分布状态的分形维数D及分形分布概率密度函数。在该基础上采用Matlab软件以及Monte-Carlo随机分析方法,产生节理参数随机数,结合3DEC中最新模块离散裂隙网络(DFN)技术,建立了能反映节理裂隙分布特征的离散裂隙网络模型并验证了模型的有效性,结果表明,分形分布比负指数分布包含更多的间距、迹长分布信息,更接近于实际分布;分形维数D反映了节理间距、迹长在其变化范围内的分布特征,分形维数的大小取决于小间距、小迹长部分数量在总节理数量中的比例,为节理裂隙岩体网络模型构建提供了一种新方法。     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

节理岩体动态破坏的SHPB相似材料试验研究

采用相似材料模型试验对不同节理倾角、节理贯通度、节理条数、载荷应变率、节理充填物厚度、节理充填物类型及试件长径比等7种工况下的节理岩体动态强度及破坏模式进行了分离式Hopkinson压杆(SHPB)试验研究。结果表明：节理岩体动态破坏模式及强度与节理构造形态密切相关。对于单节理岩体,其强度及破坏特征在很大程度上受节理倾角控制,节理倾角0°、90°试件动强度分别为完整试件的90%和71%,且其破坏形式均为张拉破坏;倾角60°试件动强度几乎为0;倾角30°、45°试件的动强度分别为完整试件的50%和18%,且其破坏以剪切破坏为主,兼有张拉破坏。中心1/4、1/2、4/5及全贯通节理试件的峰值强度分别为完整试件的95%、74% 、28%和17%,即随节理贯通度增加,试件动强度逐渐降低。含1~3条节理的试件动强度分别为完整试件的54%、23%和10%,即随节理条数增加,试件动强度随之有较大幅度降低,但节理条数的增加并没有改变其破坏模式。随着节理充填物厚度增加及节理充填物强度降低,试件强度依次递减,但破坏模式并没有改变。完整试件和节理试件的动强度均随着载荷应变率的增加而变大,且前者对载荷应变率的敏感性要远远高于后者,相应地试件的破坏模式也变得更加复杂。两类试件的动强度均随着试件长径比的增加先增大后减小,即存在一个最佳长径比。     

An examination of the validation of a model of the hydro/thermo/mechanical behaviour of engineered clay barriers

This paper focuses attention on the development of a numerical model of the hydro/thermo/mechanical behaviour of unsaturated clay and its consequent verification and validation. The work presented describes on-going collaboration between the Cardiff School of Engineering and Atomic Energy of Canada. The model development, which was carried out at Cardiff, can be described as being based on a mechanistic approach to coupled heat, moisture and air flow. This is then linked to a deformation analysis of the material within a ‘consolidation’ type of model. The whole is solved via the finite element method to yield a computer software code named COMPASS (COde for Modelling PArtly Saturated Soil). Some aspects of verification and validation of the model have been addressed in-house. However, the purpose of current AECL work is to provide an independent, rigorous, structured programme of validation and the paper will also explore the further validation of COMPASS within this context. © 1998 by John Wiley & Sons, Ltd.     

非常低延展性裂隙岩体REV存在性研究

裂隙岩体渗流问题一直是水文地质与工程地质的前沿课题之一,等效多孔介质模型是裂隙岩体渗流研究和工程应用的主要方法,此方法的重点在于裂隙岩体渗流模型即典型单元体(REV)的确定。本文采用作者自行开发的离散元软件FractureToKarst,根据国际岩石力学协会对岩石的分类,讨论了非常低延展性裂隙岩体典型单元体(REV)的存在性及大小,并给出确定其存在性的一般方法。通过研究可知,非常低延展性裂隙岩体REV存在的条件就是裂隙平均间距在0.2m以下,即极密间距、很密间距和密间距裂隙岩体存在REV。裂隙平均间距大于0.2 m,即中等间距、宽间距、很宽间距和极宽间距裂隙岩体不存在REV。     

Characterization of clastic rock using a biconcave bond model of DEM

This paper presents a biconcave bond model to investigate the effect of the cementation between grains on the mechanical behavior of rock. The proposed model considers the shape of the bonds among particles that have a biconcave cement form, based on observations of microscopic rock images. The general equations of the proposed model are based on Dvorkin theory. The accuracy and efficiency of the bond model is improved in three ways. After the biconcave bond model is implemented in the discrete element method software Particle Flow Code in 2 Dimensions, a series of numerical uniaxial compression tests were performed to investigate the relationships between the micro‐ to macro‐parameters. The simulations revealed that the biconcave bond model reflects the effect of micro‐parameters, such as the elastic modulus and Poisson's ratio of the cement, on the macroscopic deformation of cemented granular material. Variations in the bond geometry caused extremely diverse macro‐mechanical behaviors. Experimental results concerning rock demonstrate that the biconcave bond model accurately captures the mechanical behavior of intact rock and supports an innovative method for investigating the relationships between the micro‐ and macro‐parameters of cemented granular material. Copyright © 2016 John Wiley & Sons, Ltd.     

A new potential function for the calculation of contact forces in the combined finite–discrete element method

Although the potential contact force proposed by Munjiza overcomes the difficulties inherent in the traditional discrete element methods, the physical meaning of the potential is not clear and the contact force derived from the original potential function is strongly dependent on the mesh configuration. In this study, we redefine a potential function and propose a new contact force calculation method based on a unified standard. Moreover, the new potential function retains all the advantages of the original potential function but has less mesh dependency. Copyright © 2016 John Wiley & Sons, Ltd.     

基于非均质系数的裂隙岩体表征单元体研究

表征单元体（REV）是岩体力学中的一个基础性概念,其存在性是确定岩体等效参数和应用连续介质方法进行研究的前提条件。为综合反映裂隙大小、产状和密度对岩体REV的影响,本文提出采用岩体非均质系数HI来确定岩体的REV。首先,详细介绍了非均质系数的概念和含义。利用General Block软件建立了中等间距-中等延展性（MS1-MP1）裂隙的网络模型,并完成20次随机实现。从各模型中选取10个不同尺寸的岩体模型计算非均质系数,结果表明：非均质系数具有明显的尺寸效应,当岩体尺寸不小于8 m时,非均质系数的平均值和标准差分别为0.5和0.14,综合确定该岩体的REV为8 m。采用HI对三峡地下电站厂房围岩的研究表明,该岩体的REV为60 m,可作为非连续介质方法适用性的判断标准。HI是从岩体结构角度提出的新指标,适用于确定岩体的REV和统计范围,对岩体的统计分析和等效参数研究具有重要意义。     

Seismic bearing capacity of shallow foundations near rock slopes using the generalized Hoek–Brown criterion

The seismic bearing capacity of shallow foundations resting on a modified Hoek–Brown rock mass is investigated within the framework of the kinematic approach of limit analysis theory. The analysis focuses on evaluating the reduction in bearing capacity induced by seismic loading and by the proximity of a rock slope. A pseudo‐static approach is adopted to account for the earthquake effects for the seismic bearing capacity evaluations. At the rock material level, the closed‐form expressions previously obtained for the support functions of the rock failure criterion allow the implementation of different failure mechanisms families, and thus to derive rigorous upper bounds estimates of the load‐bearing capacity in both static and seismic conditions. The effects of geometrical, strength and loading parameters are assessed through a large number of parametric computations. Finally, design tables are presented for practical use in rock engineering. Copyright © 2010 John Wiley & Sons, Ltd.     

Determining the equivalent permeability tensor for fractured rock masses using a stochastic REV approach: Method and application to the field data from Sellafield,UK

A numerical procedure to determine the equivalent permeability tensor of a fractured rock is presented, using a stochastic REV (Representative Elementary Volume) concept that uses multiple realizations of stochastic DFN (Discrete Fracture Network) models. Ten square DFN models are generated using the Monte Carlo simulations of the fracture system based on the data obtained from a site characterization program at Sellafield, Cumbria, UK. Smaller models with varying sizes of from 0.25 m×0.25 m to 10 m×10 m are extracted from the generated DFN models and are used as two-dimensional geometrical models for calculation of equivalent permeability tensor. The DFN models are also rotated in 30º intervals to evaluate the tensor characteristics of calculated directional permeability. Results show that the variance of the calculated permeability values decreases significantly as the side lengths of the DFN models increase, which justifies the existence of a REV. The REV side length found in this analysis is about 5 m and 8 m with 20% and 10% acceptable variations, respectively. The calculated directional permeability values at the REV size have tensor characteristic that is confirmed by a close approximation of an ellipse in a polar plot of the reciprocal of square roots of the directional permeability.

基于量化GSI系统的破碎岩体应变软化行为研究

松软破碎岩体的复杂力学特性给深部开采工程带来一系列棘手的问题,深入研究破碎岩体峰后力学行为,对于深部巷道围岩稳定性分析及支护结构设计具有重要意义。基于量化GSI围岩评级系统及连续介质理论,建立了考虑围压影响的破碎岩体峰后应变软化力学模型。通过工程应用实例验证了数值模型的可靠性。研究结果表明：通过计算围岩与支护结构相互作用关系可知,理想弹塑性模型和应变软化模型所得计算结果差异较大,采用理想弹塑性模型进行支护设计,可能导致支护结构的安全性偏低,甚至失稳。现场应用结果表明：数值计算结果与现场实际情况吻合;提出的深部破碎岩体应变软化模型能够较为真实地反映破碎岩体的非线性破坏行为,可为硐室围岩稳定性分析及支护结构设计等提供新的思路。