A new numerical procedure for elasto-plastic analysis of a circular opening excavated in a strain-softening rock mass

A simple numerical procedure for calculating the distribution of stresses and radial displacements around a circular tunnel excavated in a strain-softening Mohr–Coulomb or generalized Hoek–Brown rock mass is proposed. The problem is considered as axisymmetric, i.e. the initial stress state is assumed to be hydrostatic and the rock mass is said to be isotropic. By invoking the finite difference approximation of the equilibrium and compatibility equations, the increments of stresses and strains for each ring, starting from the outmost one for which boundary conditions are known a priori, are calculated in a successive manner. In the proposed approach, the potential plastic zone is divided into a finite number of concentric rings whose thicknesses are determined internally to satisfy the equilibrium equation. For the strain-softening behavior, it is assumed that all the strength parameters are a linear function of deviatoric plastic strain. Several illustrative examples are given to demonstrate the performance of the proposed method. For the brittle–plastic case, the results show a very good agreement with the closed-form solution. For strain-softening cases, the predictions by the proposed method are also in good agreement with the known rigorous numerical solutions. It is shown that the approximate solution converges to the exact solution when the increment of stress for each ring becomes smaller. The influence of the strength parameter ‘a’, appearing in the generalized Hoek–Brown criterion, on the elasto-plastic solutions is examined through the establishment of ground reaction curves and the discussion for the locations of the plastic radii.     

Linearization of the Hoek and Brown rock failure criterion for tunnelling in elasto-plastic rock masses

We present a novel methodology for estimation of equivalent Mohr–Coulomb strength parameters that can be used for design of supported tunnels in elasto-plastic rock masses satisfying the non-linear empirical Hoek–Brown failure criterion. We work with a general adimensional formulation of the Hoek–Brown failure criterion in the space of normalized Lambe's variables for plane stress, and we perform linearization considering the stress field in the plastic region around the tunnel. The procedure is validated using analytical solutions to a series of benchmark test cases. Numerical solutions are also employed to validate the procedure in cases for which analytical solutions are not available. Results indicate that the stress field in the plastic region around the tunnel, as well as the linearization method employed and the quality of the rock mass, has a significant impact on computed estimates of equivalent Mohr–Coulomb strength parameters. Results of numerical analyses also show that our proposed linearization method can be employed to estimate loads and moments on the tunnel support system. We recommend the equating model responses (EMR) method to compute equivalent Mohr–Coulomb strength parameters when the tunnel support pressure is accurately known, and we further show that our newly introduced linearization method can be employed as an alternative to the best fitting in the existing stress range (BFe) and best fitting in an artificial stress range (BFa) methods, providing performance estimates that are generally better than estimates of the BFe and BFa methods when differences with the response of the Hoek–Brown rock mass are of engineering significance (say more than 10%).     

An exact implementation of the Hoek–Brown criterion for elasto-plastic finite element calculations

A simple stress update algorithm for generalised Hoek–Brown plasticity is presented. It is intended for use in elasto-plastic finite element computations and utilises the return mapping concept for computing the stress increment belonging to a given increment in strain at a material point. In the algorithm all manipulations are carried out in principal stress space, where the Hoek–Brown failure criterion has a very simple form compared to its formulation in general stress space. In principal stress space it is also simple to determine whether the stress should be returned to one of the edges or to the apex of the yield surface and to form the constitutive matrices. As opposed to earlier finite element implementations of Hoek–Brown plasticity the exact criterion is used, i.e. no rounding of the yield surface corners or edges is attempted. Numerical examples and a comparison with an often used method for dealing with the corner singularities indicates the efficiency of the presented.     

The Hoek–Brown failure criterion and GSI – 2018 edition

The Hoek–Brown criterion was introduced in 1980 to provide input for the design of underground excavations in rock. The criterion now incorporates both intact rock and discontinuities, such as joints, characterized by the geological strength index (GSI), into a system designed to estimate the mechanical behaviour of typical rock masses encountered in tunnels, slopes and foundations. The strength and deformation properties of intact rock, derived from laboratory tests, are reduced based on the properties of discontinuities in the rock mass. The nonlinear Hoek–Brown criterion for rock masses is widely accepted and has been applied in many projects around the world. While, in general, it has been found to provide satisfactory estimates, there are several questions on the limits of its applicability and on the inaccuracies related to the quality of the input data. This paper introduces relatively few fundamental changes, but it does discuss many of the issues of utilization and presents case histories to demonstrate practical applications of the criterion and the GSI system.     

20世纪岩石强度理论的发展--纪念Mohr-Coulomb强度理论100周年

190 0年莫尔 (O. Mohr)教授建立了著名的莫尔 -库仑 (Mohr- Coulomb)强度理论。 10 0年来 ,岩石强度理论广泛吸引了工程师和物理学家 (包括工程地质专家、力学家、地球物理学家、材料科学家和土木、机械工程师等 )的注意 ,对岩石在复杂应力状态下的强度理论 (破坏准则 )以及它们的实验验证进行了大量研究 ,提出了众多岩石破坏准则 ,并进行了大量轴对称围压三轴试验和真三轴试验。对 10 0年来的岩石强度理论研究的发展和主要结果进行小结 ,特别是从理论上总结了从强度理论的下限 (莫尔 -库仑单剪强度理论 ,190 0 )到上限 (双剪强度理论 ,俞茂宏 ,1985 )到统一强度理论 (俞茂宏 ,1991)的发展。对 10 0年来岩石强度理论发展的回顾 ,将有助于我们对岩石强度理论的理性认识和合理选用     

Revisiting rock classification to estimate rock mass properties

This paper presents the results of ongoing research carried out by the author exploring methods to provide a more robust estimate of rock mass properties specifically for use in tunnel design. Data from various large-scale rock mass failures are introduced, including coal pillars. The damage-initiation,spalling-limit approach is compared to the coal pillar database. New comparisons of estimating the geological strength index(GSI) and relationships to estimate the Hoeke Brown failure criterion parameters, mb, s and a, are presented.     

Interaction of shotcrete with rock and rock bolts—A numerical study

The shotcrete–rock interaction is very complex and is influenced by a number of factors. The influence of the following factors was investigated by a series of numerical analyses: the surface roughness of the opening, the rock strength and Young's modulus, the discontinuities, the extent and properties of the excavated disturbed zone, the mechanical properties of the interface between shotcrete and rock, and the thickness of the shotcrete lining and the rock bolts. The study was carried out as a sensitivity analysis. The results showed that the rock strength and the surface roughness had significant impact on the number of failures at the rock–shotcrete interface and in the shotcrete lining. Furthermore, the behaviour of the lining is sensitive to small amplitudes of the surface roughness. In all the cases investigated, a high interface strength was favourable. The results indicate that if a thick shotcrete lining is dependent on the bond strength. The benefit of using a thicker lining can be doubtful. The analyses showed that for an uneven surface the extent of the EDZ had a minor effect on the behaviour of the shotcrete lining. Furthermore, if rock bolts were installed at the apex of the protrusion instead of at the depression, the number of failures decreased both at the interface and in the lining.     

Discussion on “Empirical methods for determining shaft bearing capacity of semi-deep foundations socketed in rocks”[J Rock Mech Geotech Eng 6(2017)1140-151]

A new comprehensive set of data(n = 178) is compiled by adding a data set(n = 72) collected by Arioglu et al.(2007) to the data set(n = 106) presented in Rezazadeh and Eslami(2017). Then, the compiled data set is evaluated regardless of the variation in lithology/strength. The proposed empirical equation in this study comprises a wider range of uniaxial compressive strength(UCS)(0.15 MPa σ_(rc) 156 MPa) and various rock types. Rock mass cuttability index(RMCI) is correlated with shaft resistance(r_s) to predict the shaft resistance of rock-socketed piles. The prediction capacity of the RMCI versus r_s equation is also found to be in a fair good agreement with the presented data in Rezazadeh and Eslami(2017). Since the RMCI is a promising parameter in the prediction of shaft resistance, the researchers in the rock-socketed pile design area should consider this parameter in the further investigations.     

Physical and mechanical properties of rock masses at Stromboli: a dataset for volcano instability evaluation

Stromboli island has a complex geological history with repeated changes in the volcanic activity alternating with destructive events, caldera collapses and flank landslides. The last activity resulted in the creation of the Sciara del Fuoco depression which was modified by the recent 2002–2003 landslide. The variation in lithology, degree of tectonization and disturbance has resulted in the presence of a wide spectrum of geotechnical materials. This paper summarises the physical and mechanical properties of Stromboli’s intact rocks, rock masses and loose deposits, based on field surveys and laboratory tests. A new classification of the rock succession is introduced and four lithotechnical units defined: Lava, Lava-Breccia, Breccia and Pyroclastic deposit. The range of variability in bulk volume, porosity, intact rock compressive strength and geological strength index is presented. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, tensile and compressive strength, global strength and modulus of deformation calculated in a specified stress range.      

A nonlinear strength criterion for rock-like materials based on fracture mechanics

A nonlinear strength criterion for rock-like materials is developed in this paper. Taking α as an angle of micro-failure orientation in rock-like materials, a formulation between α and load is derived from a mixed-mode fracture criterion based on linear elastic fracture mechanics. According to micro-failure experimental phenomena of rock-like materials, a failure characteristic parameter under triaxial compression condition is chosen, which is relevant to confining pressure and is an invariant. A theoretical nonlinear strength criterion is also derived, which is exactly in the same mathematical form as the original Hoek–Brown empirical strength criterion. In addition, it is also found that the coefficient m in the Hoek–Brown criterion has physical meaning which is related to the ratio between the uniaxial compressive strength and the uniaxial tensile strength.     

Overhanging rock slope by design: An integrated approach using rock mass strength characterisation,large-scale numerical modelling and limit equilibrium methods

Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.     

Effects of confinement on rock mass modulus: A synthetic rock mass modelling (SRM) study

The main objective of this paper is to examine the influence of the applied confining stress on the rock mass modulus of moderately jointed rocks (well interlocked undisturbed rock mass with blocks formed by three or less intersecting joints). A synthetic rock mass modelling (SRM) approach is employed to determine the mechanical properties of the rock mass. In this approach, the intact body of rock is represented by the discrete element method (DEM)-Voronoi grains with the ability of simulating the initiation and propagation of microcracks within the intact part of the model. The geometry of the pre-existing joints is generated by employing discrete fracture network (DFN) modelling based on field joint data collected from the Brockville Tunnel using LiDAR scanning. The geometrical characteristics of the simulated joints at a representative sample size are first validated against the field data, and then used to measure the rock quality designation (RQD), joint spacing, areal fracture intensity (P₂₁), and block volumes. These geometrical quantities are used to quantitatively determine a representative range of the geological strength index (GSI). The results show that estimating the GSI using the RQD tends to make a closer estimate of the degree of blockiness that leads to GSI values corresponding to those obtained from direct visual observations of the rock mass conditions in the field. The use of joint spacing and block volume in order to quantify the GSI value range for the studied rock mass suggests a lower range compared to that evaluated in situ. Based on numerical modelling results and laboratory data of rock testing reported in the literature, a semi-empirical equation is proposed that relates the rock mass modulus to confinement as a function of the areal fracture intensity and joint stiffness.     

Shear strength criteria for rock,rock joints,rockfill and rock masses: Problems and some solutions

Although many intact rock types can be very strong, a critical confining pressure can eventually be reached in triaxial testing, such that the Mohr shear strength envelope becomes horizontal. This critical state has recently been better defined, and correct curvature or correct deviation from linear Mohr–Coulomb (M-C) has finally been found. Standard shear testing procedures for rock joints, using multiple testing of the same sample, in case of insufficient samples, can be shown to exaggerate apparent cohesion. Even rough joints do not have any cohesion, but instead have very high friction angles at low stress, due to strong dilation. Rock masses, implying problems of large-scale interaction with engineering structures, may have both cohesive and frictional strength components. However, it is not correct to add these, following linear M-C or nonlinear Hoek–Brown (H-B) standard routines. Cohesion is broken at small strain, while friction is mobilized at larger strain and remains to the end of the shear deformation. The criterion ‘c then σ_n tan φ’ should replace ‘c plus σ_ntan φ’ for improved fit to reality. Transformation of principal stresses to a shear plane seems to ignore mobilized dilation, and caused great experimental difficulties until understood. There seems to be plenty of room for continued research, so that errors of judgement of the last 50 years can be corrected.     

Reliability-based design in rock engineering: Application of Bayesian regression methods to rock strength data

Reliability-based design(RBD) is being adopted by geotechnical design codes worldwide, and it is therefore necessary that rock engineering practice evolves to embrace RBD. This paper examines the Hoek-Brown(He B) strength criterion within the RBD framework, and presents three distinct analyses using a Bayesian approach. Firstly, a compilation of intact compressive strength test data for six rock types is used to examine uncertainty and variability in the estimated He B parameters m and sc, and corresponding predicted axial strength. The results suggest that within-and between-rock type variabilities are so large that these parameters need to be determined from rock testing campaigns, rather than reference values being used. The second analysis uses an extensive set of compressive and tensile(both direct and indirect) strength data for a granodiorite, together with a new Bayesian regression model, to develop joint probability distributions of m and scsuitable for use in RBD. This analysis also shows how compressive and indirect tensile strength data may be robustly used to fit an He B criterion.The third analysis uses the granodiorite data to investigate the important matter of developing characteristic strength criteria. Using definitions from Eurocode 7, a formal Bayesian interpretation of characteristic strength is proposed and used to analyse strength data to generate a characteristic criterion.These criteria are presented in terms of characteristic parameters mkand s ck, the values of which are shown to depend on the testing regime used to obtain the strength data. The paper confirms that careful use of appropriate Bayesian statistical analysis allows the He B criterion to be brought within the framework of RBD. It also reveals that testing guidelines such as the International Society for Rock Mechanics and Rock Engineering(ISRM) suggested methods will require modification in order to support RBD. Importantly, the need to fully understand the implications of uncertainty in nonlinear strength criteria is identified.     

一个用于预测金属板材断裂的损伤破坏准则

基于连续介质损伤力学 ,提出了描述金属板材塑性损伤的演变方程 ,并依此得到了一个用于预测金属板材断裂破坏的简单准则 .其中特别考虑了剪切机制和剪切应力在金属薄板损伤发展和最终破坏过程中的作用 ,同时也引入了材料微观损伤对材料机械性能的影响 .提出的损伤演变规律 ,可近似表征金属薄板在成型过程中损伤的发展变化 ,计算得出的损伤曲线与有关实验曲线比较吻合 .给出的损伤破坏准则 ,可直接用于预测金属板材的断裂极限应变     

Behavior of dam concrete under biaxial compression-tension and triaxial compression-compression-tension stresses

In order to meet the requirement for nonlinear analysis and design of mass concrete structures, the deformation behavior and strength of three-graded concrete specimens 250 mm × 250 mm × 400 mm with a maximum aggregate size of 80 mm and the corresponding wet-screened concrete specimens 150 mm × 150 mm × 300 mm with a maximum aggregate size of 40 mm were studied experimentally. Specimens subjected to biaxial compression-tension (C-T) and triaxial compression-compression-tension (C-C-T) stress states. Test data indicate that both the deformation and strength of the mass concrete specimens are lower than those of the corresponding wetscreened concrete small specimens, but the initial tangent modulus of the stress-strain curve of the former is greater than that of the latter. Test results show that the wetscreened effect and size effect of the specimens under complex stress states are obvious such that these should be considered in the design of mass concrete structures. In addition, respective failure criteria for mass concrete in principal stress space and octahedron stress space are proposed. __________ Translated from China Civil Engineering Journal, 2007, 40(7): 104–110 [译自: 土木工程学报]     

A simplified three-dimensional extension of Hoek-Brown strength criterion

The Hoek-Brown(HB) strength criterion has been applied widely in a large number of projects around the world.However,this criterion ignores the intermediate principal stress σ_2.Many evidences have demonstrated that the rock strength is dependent on σ_2. Thus it is necessary to extend the HB criterion into a three-dimensional(3D) form.In this study,the effect of σ_2 on the strength of rocks is identified by reviewing the true triaxial tests of various rock types reported in the literature.A simple 3D strength criterion is developed.The modified criterion is verified by the true triaxial tests of 13 rock types.The results indicate that the modified criterion can achieve a good fit to most of rock types.It can represent a series of criteria as b varies.For comparisons,several existing 3D versions of the HB criterion are selected to predict the strengths of these rock types.It is indicated that the proposed criterion works better than other criteria.A substantial relationship between parameter b and the unconfined compressive strength is established,which guarantees that the proposed criterion can still work well even in the absence of true triaxial test data.     

碾压混凝土抗弯拉强度特性的试验研究

通过室内小梁试验对用水量和水泥用量较大的普通碾压混凝土抗弯拉强度特性进行了研究,建立了各龄期抗弯拉强度与显著影响因素之间的经验回归公式,都具有较好的相关关系,并分析了各因素对抗弯拉强度的影响,试验研究对以碾压混凝土早期强度为设计标准的简易机场和公路结构设计具有指导意义.