A Bonded Particle Model Simulation of Shear Strength and Asperity Degradation for Rough Rock Fractures

Different failure modes during fracture shearing have been introduced including normal dilation or sliding, asperity cut-off and degradation. Attempts have been made to study these mechanisms using analytical, experimental and numerical methods. However, the majority of the existing models simplify the problem, which leads to unrealistic results. With this in mind, the aim of this paper is to simulate the mechanical behaviour of synthetic and rock fracture profiles during direct shear tests by using the two-dimensional particle flow computer code PFC2D. Correlations between the simulated peak shear strength and the fracture roughness parameter D _R1 recently proposed by Rasouli and Harrison (2010) are developed. Shear test simulations are carried out with PFC2D and the effects of the geometrical features as well as the model micro-properties on the fracture shear behaviour are studied. The shear strength and asperity degradation processes of synthetic profiles including triangular, sinusoidal and randomly generated profiles are analysed. Different failure modes including asperity sliding, cut-off, and asperity degradation are explicitly observed and compared with the available models. The D _R1 parameter is applied to the analysis of synthetic and rock fracture profiles. Accordingly, correlations are developed between D _R1 and the peak shear strength obtained from simulations and by using analytical solutions. The results are shown to be in good agreement with the basic understanding of rock fracture shear behaviour and asperity contact degradation.     

Numerical Investigation of the Dynamic Compressive Behaviour of Rock Materials at High Strain Rate

The dynamic compressive strength of rock materials increases with the strain rate. They are usually obtained by conducting laboratory tests such as split Hopkinson pressure bar (SHPB) test or drop-weight test. It is commonly agreed now that the dynamic increase factor (DIF) obtained from impact test is affected by lateral inertia confinement, friction confinement between the specimen and impact materials and the specimen sizes and geometries. Therefore, those derived directly from testing data do not necessarily reflect the true dynamic material properties. The influences of these parameters, however, are not straightforward to be quantified in laboratory tests. Therefore, the empirical DIF relations of rock materials obtained directly from impact tests consist of contributions from lateral inertia and end friction confinements, which need be eliminated to reflect the true dynamic material properties. Moreover, different rocks, such as granite, limestone and tuff have different material parameters, e.g., equation of state (EOS) and strength, which may also affect the DIF of materials but are not explicitly studied in the open literature. In the present study, numerical models of granite, limestone and tuff materials with different EOS and strength under impact loads are developed to simulate SHPB tests and to study the influences of EOS and strength, lateral inertia confinement and end friction confinement effects on their respective DIFs in the strain rate range between 1 and 1,000 s^?1. The commercial software AUTODYN with user-provided subroutines is used to perform the numerical simulations of SHPB tests. Numerical simulation results indicate that the lateral inertia confinement, friction confinement and specimen aspect (L/D) ratio significantly influence DIF obtained from impact tests and the inertia confinement effect is different for different rocks. Based on the numerical results, quantifications on the relative contributions from the lateral inertia confinement and the material strain rate effect on DIF of granite, limestone and tuff material compressive strength are made. The effects of friction coefficient, L/D ratio and rock type on DIF are discussed. Empirical relations of DIF with strain rate for the three rock materials representing the true material strain rate effect are also proposed.     

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.     

Evidence for a Long-Term Strength Threshold in Crystalline Rock

The mechanical response of brittle rock to long-duration compression loading is of particular concern in underground disposal of nuclear waste, where radionuclides must be isolated from the biosphere for periods of the order of a million years. Does the strength decrease without limit over such time, or is there, for some rock types, a lower “threshold” strength below which the rock will cease to deform? This paper examines the possibility of such a threshold in silicate crystalline rocks from several perspectives, including: (1) interpretation of the results of short-term creep tests on rock; (2) numerical analysis of the effect of decrease in fracture toughness due to stress corrosion on the strength of a crystalline rock; and (3) evidence from plate tectonics, and observations of in situ rock stress in granite quarries. The study concludes that there is clear evidence of threshold strength. The threshold is of the order of 40% of the unconfined compressive strength or higher for laboratory specimens under unconfined compressive loading, and increases rapidly in absolute value with confinement. Field evidence also leads to the conclusion that the long-term strength of crystalline rock in situ is of comparable magnitude to the laboratory value.     

交叉裂隙岩体裂纹扩展试验及混合有限-离散元数值模拟研究

对预制交叉裂隙岩石试样进行裂纹扩展试验,研究了裂纹萌生、扩展、聚合过程,分析了主裂隙和轴向载荷夹角及主次裂隙夹角对裂纹起裂应力和聚合应力的影响,并用混合有限元-离散元程序,即图形处理器并行化的3D Y-HFDEM代码对试验进行了仿真计算,实现了岩石破坏从连续介质向非连续介质的过渡,对裂纹类及损伤破坏模式进行了识别,捕捉到了试验中难以发现的现象。研究表明：随主裂隙与轴向载荷夹角增加,裂纹聚合区的拉伸裂纹数量增加;裂纹起裂和聚合应力与主裂隙与轴向载荷夹角成正比;主次裂隙夹角增加,岩石的破坏模式由拉伸破坏转为剪切破坏,交叉裂隙加剧岩石破碎程度;主裂隙尖端萌生扩展的拉伸-剪切混合裂缝引起的破坏在岩石破坏中占主导地位,是导致岩体失去承载能力的主控裂纹;混合有限元-离散元仿真软件GPGPU并行化的3D Y-HFDEM IDE在岩石裂纹扩展研究中具有优势,可以捕捉实验室难以发现的损伤断裂类型,可以作为岩石裂纹扩展研究的有力工具。     

Influence of inclination angles and confining pressures on mechanical behavior of rock materials containing a preexisting crack

To deeply understand the cracking mechanical behavior of brittle rock materials, numerical simulations of a rock specimen containing a single preexisting crack were carried out by the expanded distinct element method (EDEM). Based on the analysis of crack tips and a comparison between stress- and strain-based methods, the strain strength criterion was adopted in the numerical models to simulate the crack initiation and propagation processes under uniaxial and biaxial compression. The simulation results indicated that the crack inclination angle and confining pressure had a great influence on the tensile and shear properties, peak strength, and failure behaviors, which also showed a good agreement with the experimental results. If the specimen was under uniaxial compression, it was found that the initiation stress and peak strength first decreased and then increased with an increasing inclination angle α. Regardless of the size of α, tensile cracks initiated prior to shear cracks. If α was small (such as α ≤ 30°), the tensile cracks dominated the specimen failure, the wing cracks propagated towards the direction of uniaxial compression, and the propagation of shear cracks was inhibited by the high concentration of tensile stress. In contrast, if α was large (such as α ≥ 45°), mixed cracks dominated the specimen failure, and the external loading favored the further propagation of shear cracks. Analyzing the numerical results of the specimen with a 45° inclination angle under biaxial compression, it was revealed that lateral confinement had a significant influence on the initiation sequence and the mechanical properties of new cracks.     

FDEM-TH3D: A three-dimensional coupled hydrothermal model for fractured rock

This paper proposes a three-dimensional coupled hydrothermal model for fractured rock based on the finite-discrete element method to simulate fluid flow and heat transport. The 3D coupled hydrothermal model is composed of three main parts: a heat conduction model for the rock matrix, a heat transfer model for the fluid in the fractures (including heat conduction and heat convection), and a heat exchange model between the rock matrix and the fluid in the fractures. Four examples with analytical solutions are provided to verify the model. A heat exchange experiment of circulating water in a cylindrical granite sample with one fracture is simulated. The simulation results agree well with the experimental results. The effects of the fracture aperture, fluid viscosity, and pressure difference on the heat exchange between the fluid and rock are studied. Finally, an application concerned with heat transport and fluid flow in fractured rock is presented. The simulation results indicate that the 3D fully coupled hydrothermal model can capture the fluid flow and temperature evolution of rocks and fluids.     

Polyaxial stress-dependent permeability of a three-dimensional fractured rock layer

A study about the influence of polyaxial (true-triaxial) stresses on the permeability of a three-dimensional (3D) fractured rock layer is presented. The 3D fracture system is constructed by extruding a two-dimensional (2D) outcrop pattern of a limestone bed that exhibits a ladder structure consisting of a “through-going” joint set abutted by later-stage short fractures. Geomechanical behaviour of the 3D fractured rock in response to in-situ stresses is modelled by the finite-discrete element method, which can capture the deformation of matrix blocks, variation of stress fields, reactivation of pre-existing rough fractures and propagation of new cracks. A series of numerical simulations is designed to load the fractured rock using various polyaxial in-situ stresses and the stress-dependent flow properties are further calculated. The fractured layer tends to exhibit stronger flow localisation and higher equivalent permeability as the far-field stress ratio is increased and the stress field is rotated such that fractures are preferentially oriented for shearing. The shear dilation of pre-existing fractures has dominant effects on flow localisation in the system, while the propagation of new fractures has minor impacts. The role of the overburden stress suggests that the conventional 2D analysis that neglects the effect of the out-of-plane stress (perpendicular to the bedding interface) may provide indicative approximations but not fully capture the polyaxial stress-dependent fracture network behaviour. The results of this study have important implications for understanding the heterogeneous flow of geological fluids (e.g. groundwater, petroleum) in subsurface and upscaling permeability for large-scale assessments.     

不同倾角节理组和锚固效应对岩体特性的影响

首先采用DDARF(discontinuous deformation analysis for rock failure)分析方法对双裂隙岩块进行单轴和双轴压缩模拟试验,研究了裂隙角度和侧向应力大小对岩块特性的影响,得到了裂隙岩块在这两种加载试验中的破坏过程、应力-应变曲线以及岩块中裂隙的起裂应力和岩块的峰值强度。在双轴压缩模拟试验中绘制了裂隙角度为45°的岩块在不同侧向压力下的强度包络线。其次,采用DDARF分析方法模拟劈裂试验中含裂隙试块的锚固效果,得到了4种不同锚固角度试块的轴向荷载–位移变化曲线和裂隙扩展规律。模拟结果与前人的类似条件下的试验结果相符良好。随后又将双裂隙试块双轴压缩模拟试验中得到的参数运用到一个地下洞室的工程实例中,用等效力学特性的方法分析对比了完整岩体和节理岩体洞室开挖完成后的破损状态的差异。最后运用DDARF分析方法研究了随机生成4组节理岩体的地下洞室的稳定性,得到了洞室节理围岩的裂隙扩展过程。同时通过对关键点位移的监测分析了锚杆的锚固效应。     

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 the Anisotropic Mechanical Behaviour of Opalinus Clay at the Laboratory-Scale Using FEM/DEM

The Opalinus Clay (OPA) is an argillaceous rock formation selected to host a deep geologic repository for high-level nuclear waste in Switzerland. It has been shown that the excavation damaged zone (EDZ) in this formation is heavily affected by the anisotropic mechanical response of the material related to the presence of bedding planes. In this context, the purpose of this study is twofold: (i) to illustrate the new developments that have been introduced into the combined finite-discrete element method (FEM/DEM) to model layered materials and (ii) to demonstrate the effectiveness of this new modelling approach in simulating the short-term mechanical response of OPA at the laboratory-scale. A transversely isotropic elastic constitutive law is implemented to account for the anisotropic elastic modulus, while a procedure to incorporate a distribution of preferentially oriented defects is devised to capture the anisotropic strength. Laboratory results of indirect tensile tests and uniaxial compression tests are used to calibrate the numerical model. Emergent strength and deformation properties, together with the simulated damage mechanisms, are shown to be in strong agreement with experimental observations. Subsequently, the calibrated model is validated by investigating the effect of confinement and the influence of the loading angle with respect to the specimen anisotropy. Simulated fracture patterns are discussed in the context of the theory of brittle rock failure and analyzed with reference to the EDZ formation mechanisms observed at the Mont Terri Underground Research Laboratory.     

Simulations of hydro-fracking in rock mass at meso-scale using fully coupled DEM/CFD approach

The paper deals with two-dimensional (2D) numerical modelling of hydro-fracking (hydraulic fracturing) in rocks at the meso-scale. A numerical model was developed to characterize the properties of fluid-driven fractures in rocks by combining the discrete element method (DEM) with computational fluid dynamics (CFD). The mechanical behaviour of the rock matrix was simulated with DEM and the behaviour of the fracturing fluid flow in newly developed and pre-existing fractures with CFD. The changes in the void geometry in the rock matrix were taken into account. The initial 2D hydro-fracking simulation tests were carried out for a rock segment under biaxial compression with one injection slot in order to validate the numerical model. The qualitative effect of several parameters on the propagation of a hydraulic fracture was studied: initial porosity of the rock matrix, dynamic viscosity of the fracking fluid, rock strength and pre-existing fracture. The characteristic features of a fractured rock mass due to a high-pressure injection of fluid were realistically modelled by the proposed coupled approach.

     

Granite rock fragmentation at percussive drilling – experimental and numerical investigation

The aim of this study is to numerically model the fracture system at percussive drilling. Because of the complex behavior of rock materials, a continuum approach is employed relying upon a plasticity model with yield surface locus as a quadratic function of the mean pressure in the principal stress space coupled with an anisotropic damage model. In particular, Bohus granite rock is investigated, and the material parameters are defined based on previous experiments. This includes different tests such as direct tension and compression, three‐point bending, and quasi‐oedometric tests to investigate the material behavior at both tension and confined compression stress states. The equation of motion is discretized using a finite element approach, and the explicit time integration method is employed. Edge‐on impact tests are performed, and the results are used to validate the numerical model. The percussive drilling problem is then modeled in 3D, and the bit‐rock interaction is considered using contact mechanics. The fracture mechanism in the rock and the bit penetration‐ resisting force response are realistically captured by the numerical model. Copyright © 2013 John Wiley & Sons, Ltd.     

高温作用后花岗岩三轴压缩试验研究

为综合考察温度、围压对花岗岩力学性质及破坏方式的影响,在高温（25℃～1 000 ℃）作用后,利用MTS815.02电液伺服材料试验系统对花岗岩岩样进行不同围压作用下的三轴压缩试验。研究结果表明,（1）围压一定时,经历不同高温作用后花岗岩三轴压缩全应力-应变曲线经历了压密、弹性、屈服、破坏、塑性流动5个阶段;（2）经历不同高温作用后岩样三轴抗压强度与围压呈非线性二次多项式增长关系,围压为40 MPa时的抗压强度比单轴抗压强度提高了382.30%;常规三轴压缩条件下,400 ℃是花岗岩力学参数的阀值温度;（3）经历高温作用后,岩样弹性模量随围压升高呈增大趋势,围压为40 MPa时的弹性模量比单轴时提高了90.26%;随温度升高呈二次非线性减小,1 000 ℃时的弹性模量比25℃时降低了57.16%;（4）花岗岩的失稳型式同时取决于围压和温度。单轴压缩状态下,随着温度的升高,岩样变形破坏型式由脆性破裂向塑性变形过渡,失稳型式在低温时为突发失稳、中高温为准突发失稳,温度高于800 ℃为渐进破坏;三轴压缩状态下,随着围压的增大,岩样破裂型式由脆性张拉破裂逐渐向剪切破裂过渡,岩样的失稳型式以突发失稳为主。在试验温压范围内,影响花岗岩力学性质的首要因素是温度,其次是围压。     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

双轴应力下非贯通节理岩体压缩损伤本构模型

针对目前节理岩体损伤变量定义中大多仅考虑节理长度、倾角等几何性质,而未考虑节理抗剪强度等力学性质的不足,基于断裂力学中的由于单个节理存在引起的附加应变能增量与损伤力学中的损伤应变能释放量相关联的观点,推导出了在双轴应力下含单条非贯通闭合节理岩体的损伤变量计算公式,并根据断裂力学理论对双轴压缩荷载下的单个节理尖端应力强度因子计算方法进行了研究,得出了应力强度因子KⅠ、KⅡ的计算公式;同时考虑多节理间的相互作用,给出了单组单排及单组多排非贯通节理尖端应力强度因子计算公式,由此建立了相应的节理岩体双轴压缩损伤本构模型,并利用该模型进行了相应的算例分析。结果表明：对含单条非贯通闭合节理的岩体而言,当节理倾角小于其内摩擦角时,岩体强度与完整岩石相同,岩体损伤为0,而后随着节理倾角增加,岩体强度、损伤随节理倾角的变化分别呈开口向上及向下的抛物线,当节理倾角约为60°时,岩体损伤最大,强度最低。随着节理长度增加,岩体损伤增加,而随着节理内摩擦角的增加,岩体损伤则减小;对含单组单排非贯通闭合节理的岩体而言,当节理总长度一定时,随着单条节理长度的减小及节理条数的增加,岩体损伤则逐渐减小,但其减小幅度与节理条数并不呈线性关系。     

An Empirical Failure Criterion for Intact Rocks

The parameter m _i is an important rock property parameter required for use of the Hoek–Brown failure criterion. The conventional method for determining m _i is to fit a series of triaxial compression test data. In the absence of laboratory test data, guideline charts have been provided by Hoek to estimate the m _i value. In the conventional Hoek–Brown failure criterion, the m _i value is a constant for a given rock. It is observed that using a constant m _i may not fit the triaxial compression test data well for some rocks. In this paper, a negative exponent empirical model is proposed to express m _i as a function of confinement, and this exercise leads us to a new empirical failure criterion for intact rocks. Triaxial compression test data of various rocks are used to fit parameters of this model. It is seen that the new empirical failure criterion fits the test data better than the conventional Hoek–Brown failure criterion for intact rocks. The conventional Hoek–Brown criterion fits the test data well in the high-confinement region but fails to match data well in the low-confinement and tension regions. In particular, it overestimates the uniaxial compressive strength (UCS) and the uniaxial tensile strength of rocks. On the other hand, curves fitted by the proposed empirical failure criterion match test data very well, and the estimated UCS and tensile strength agree well with test data.     

Depth Variability of Compressive Strength Test Results of Toki Granite from Shobasama and Mizunami Construction Sites,Japan

This paper studies the depth variability of uniaxial compressive laboratory test results on intact Toki granite (i.e., sound rock without macroscopic fractures) from the Shobasama and Mizunami Construction Sites, Japan. Some of the depth variability observed in the laboratory results can be indirectly attributed to the high fracture frequency of the “upper highly fractured rock domain” from which some of the samples were taken. For samples taken from the “lower sparsely fractured rock domain,” however, the uniaxial compressive strength of the granite seems to be very strongly correlated to the level of in situ rock stress (i.e., maximum shear stress) determined by measurement results obtained from hydrofracturing tests. The correlation between the laboratory results and the level of in situ stress is explained by the damage due to the complex stress path that the cores undergo during drilling, besides the stress concentrations at the drill-bit/rock contact, which can also affect the microcracking of the samples. An attempt to adjust laboratory test results to estimate the in situ intact rock strength of Toki granite based on its correlation with in situ stresses was carried out.     

A 2D coupled hydro-thermal model for the combined finite-discrete element method

Based on the combined finite-discrete element method (FDEM), a two-dimensional coupled hydro-thermal model is proposed. This model can simulate fluid flow and heat transfer in rock masses with arbitrary complex fracture networks. The model consists of three parts: a heat conduction model of the rock matrix, a heat-transfer model of the fluid in the fracture (including the heat conduction and convection of fluid), and a heat exchange model between the fluid and rock at the fracture surface. Three examples with analytical solutions are given to verify the correctness of the coupled model. Finally, the coupled model is applied to hydro-thermal coupling simulations of a rock mass with a fracture network. The temperature field evolution, the effect of thermal conductivity of the rock matrix thermal conductivity and the fracture aperture on the outlet temperature are studied. The coupled model presented in this paper will enable the application of FDEM to study rock rupture driven by the effect of hydro-thermo-mechanical coupling in geomaterials such as in geothermal systems, petroleum engineering, environmental engineering and nuclear waste geological storage.

     

Correlation of production rate of ornamental stone with rock brittleness indexes

The prediction of production rate in ornamental stones sawing is very important in cost estimation and process planning of the rock sawing plants. The main aim of this paper is finding a mathematical correlation between production rate and rock brittleness indexes. The utilized data have been collected from several stone factories in Iran. Seventeen different granite and carbonate rocks have been experienced sawing conditions with large-diameter circular saws. The laboratory tests such as uniaxial compressive strength and tensile strength, were carried out on the rock samples which were collected from these factories. The ratio of compressive strength to tensile strength (B ₁), the ratio of compressive strength minus tensile strength to compressive strength plus tensile strength (B ₂), and the area under the line of compressive strength and the line of tensile strength (B ₃) were used as brittleness indexes in this study. Then, the relationships between these brittleness indexes and production rate were investigated by using a regression analysis. The results show that, there is not a good relationship between the production rate and B ₁ and B ₂. However, the production rate correlates with B ₃ very well. Therefore, there is a reliable prediction for ornamental stones production based on B ₃ as a brittleness indexes.