Three-dimensional limit variation analysis on the ultimate pullout capacity of the anchor cables based on the Hoek-Brown failure criterion

Only simplified two-dimensional model and a single failure mode are adopted to calculate the ultimate pullout capacity(UPC)of anchor cables in most previous research.This study focuses on a more comprehensive combination failure mode that consists of bond failure of an anchorage body and failure of an anchored rock mass.The three-dimensional ultimate pullout capacity of the anchor cables is calculated based on the Hoek-Brown failure criterion and variation analysis method.The numerical solution for the curvilinear function in fracture plane is obtained based on the finite difference theory,which more accurately reflects the failure state of the anchor cable,as opposed to that being assumed in advance.The results reveal that relying solely on a single failure mode for UPC calculations has limitations,as changes in parameter values not only directly impact the UPC value but also can alter the failure model and thus the calculation method.

Energy analysis of rock plug thickness in karst tunnels based on non-associated flow rule and nonlinear failure criterion

The geological hazards, such as water inrush and mud outburst, are easily induced by the high water pressure caverns ahead of a karst tunnel face. Therefore, it is a pivotal issue to determine the reserved thickness of rock plug during the construction of tunnels. The limit analysis principle is employed to analyze the safe thickness from the point of energy dissipation, and the nonlinear and non-associated characteristics of geotechnical materials are both considered. On the basis of a plane failure pattern of rock plug, the expressions of detaching curve and rock plug thickness are derived. The effect of each parameter on the safe thickness of rock plug is discussed in detail, which interprets the corresponding failure scope of rock plug. The obtained results indicate that the thickness of rock plug is highly influenced by the nonlinear dilatancy coefficient and the nonlinear coefficient. The proposed method is validated by a comparison of the calculated results with those of the engineering project of the "526 karst cavern" of Yunwushan tunnel. This proposed method can provide reference basis for the design and excavation of karst tunnels in the future.

Validation of Hoek-Brown failure criterion charts for rock slopes

Although stability charts suggested by Hoek and Bray on the basis of Mohr-Coulomb criterion are used for rock slopes,but complete and precise recognition is required for distinguishing cohesive strength and Mohr-Coulomb equivalent internal friction angle for rock mass.The paper by Lia et al.[6]is the only one that introduced rock slope charts according to Hoek-Brown failure criterion.In this paper,at first,this type of charts is introduced.Then,Mohr-Coulomb failure criterion charts[2]are compared and validated with Hoek-Brown failure criterion ones[6].Next,Bishop method utilizing Slide software is compared with Hoek-Brown failure criterion stability charts.Average standard deviation(ASD),root mean square error(RMSE)and variance account for(VAF)were used for the comparison.According to the results,because of high distribution and very low correlation among the comparisons,Hoek-Brown failure criterion charts are not efficient.

Upper bound solutions for supporting pressures of shallow tunnels with nonlinear failure criterion

线性失败标准广泛地在对浅隧道起作用的地球压力的计算被使用。然而，试验性的证据证明那个非线性的失败标准能表示相当，很好，几乎所有的失败岩石打字。一个非线性的 Hoek 布郎失败标准被采用在限制分析框架估计浅隧道的支持的压力。二失败机制为计算外部力量和内部精力驱散的工作率被建议。到非线性的失败标准的一根正切的线被用来作为一个非线性的编程问题提出支持的压力问题。这样提出的客观功能关于失败机制和接触点的地点被最小化。以便估计有效性，为建议失败机制的支持的压力被计算，当非线性的标准被归结为线性标准时，与对方，和礼品相比，结果与以前出版的答案相比。协议支持建议失败机制的有效性。一个实验被进行在倒塌形状和浅隧道的支持的压力上调查非线性的标准的影响。

Kinematical analysis of highway tunnel collapse using nonlinear failure criterion

For different kinds of rocks,the collapse range of tunnel was studied in the previously published literature.However,some tunnels were buried in soils,and test data showed that the strength envelopes of the soils followed power-law failure criterion.In this work,deep buried highway tunnel with large section was taken as objective,and the basic expressions of collapse shape and region were deduced for the highway tunnels in soils,based on kinematical approach and power-law failure criterion.In order to see the effectiveness of the proposed expressions,the solutions presented in this work agree well with previous results if the nonlinear failure criterion is reduced to a linear Mohr-Coulomb failure criterion.The present results are compared with practical projects and tunnel design code.The numerical results show that the height and width of tunnel collapse are greatly affected by the nonlinear criterion for the tunnel in soil.

Upper bound analysis of slope stability with nonlinear failure criterion based on strength reduction technique

Based on the upper bound limit analysis theorem and the shear strength reduction technique,the equation for expressing critical limit-equilibrium state was employed to define the safety factor of a given slope and its corresponding critical failure mechanism by means of the kinematical approach of limit analysis theory.The nonlinear shear strength parameters were treated as variable parameters and a kinematically admissible failure mechanism was considered for calculation schemes.The iterative optimization method was adopted to obtain the safety factors.Case study and comparative analysis show that solutions presented here agree with available predictions when nonlinear criterion reduces to linear criterion,and the validity of present method could be illuminated.From the numerical results,it can also be seen that nonlinear parameter m,slope foot gradient β,height of slope H,slope top gradient α and soil bulk density γ have significant effects on the safety factor of the slope.

Safe retaining pressures for pressurized tunnel face using nonlinear failure criterion and reliability theory

Based on the active failure mechanism and passive failure mechanism for a pressurized tunnel face, the analytical solutions of the minimum collapse pressure and maximum blowout pressure that could maintain the stability of pressurized tunnel faces were deduced using limit analysis in conjunction with nonlinear failure criterion under the condition of pore water pressure. Due to the objective existence of the parameter randomness of soil, the statistical properties of random variables were determined by the maximum entropy principle, and the Monte Carlo method was employed to calculate the failure probability of a pressurized tunnel. The results show that the randomness of soil parameters exerts great influence on the stability of a pressurized tunnel, which indicates that the research should be done on the topic of determination of statistical distribution for geotechnical parameters and the level of variability. For the failure probability of a pressurized tunnel under multiple failure modes, the corresponding safe retaining pressures and optimal range of safe retaining pressures are calculated by introducing allowable failure probability and minimum allowable failure probability. The results can provide practical use in the pressurized tunnel engineering.

Energy analysis of stability of twin shallow tunnels based on nonlinear failure criterion

Based on nonlinear Mohr-Coulomb failure criterion, the analytical solutions of stability number and supporting force on twin shallow tunnels were derived using upper bound theorem of limit analysis. The optimized solutions were obtained by the technique of sequential quadratic programming. When nonlinear coefficient equals 1 and internal friction angle equals 0, the nonlinear Mohr-Coulomb failure criterion degenerates into linear failure criterion. The calculated results of stability number in this work were compared with previous results, and the agreement verifies the effectiveness of the present method. Under the condition of nonlinear Mohr-Coulomb failure criterion, the results show that the supporting force on twin shallow tunnels obviously increases when the nonlinear coefficient, burial depth, ground load or pore water pressure coefficients increase. When the clear distance is 0.5to 1.0 times the diameter of tunnel, the supporting force of twin shallow tunnels reaches its maximum value, which means that the tunnels are the easiest to collapse. While the clear distance increases to 3.5 times the diameter of tunnel, the calculation for twin shallow tunnels can be carried out by the method for independent single shallow tunnel. Therefore, 3.5 times the diameter of tunnel serves as a critical value to determine whether twin shallow tunnels influence each other. In designing twin shallow tunnels,appropriate clear distance value must be selected according to its change rules and actual topographic conditions, meanwhile, the influences of nonlinear failure criterion of soil materials and pore water must be completely considered. During the excavation process, supporting system should be intensified at the positions of larger burial depth or ground load to avoid collapses.

A generalized three-dimensional failure criterion for rock masses

The smooth convex generalized failure function, which represents 1/6 part of envelope in tile deviatoric plane, is proposed. The proposed function relies on four shape parameters （L, a, b and c）, in which two parameters （a and b） are dependent on the others. The parameter Ls is called extension ratio. The proposed failure function could be incorporated with any two-dimensional （2D） failure criteria to make it a three-dimensional （3D） version. In this paper, a mathematical formulation for incorporation of Hoek-Brown failure criterion with the proposed function is presented. The Hoek-Brown failure criterion is the most suited 2D failure criterion tbr geomaterials. Two types of analyses for best-fitting solution of published true tri-axial test data were made by considering （1） constant extension ratio and （2） variable extension ratio. The shape and strength parameters for different types of rocks have been determined by best-fitting the published true tri-axial test data for both the analyses. It is observed from the best-fitting solution by considering uniform extension ratio （L~） that shape constants have a correlation with Hoek-Brown strength parameters. Thus, only two parameters （c~. and m） are needed for representing the 3D failure criterion for intact rock. The statistical expression between shape and Hoek-Brown strength parameters is given. In the second analysis, when considering varying extension ratio, another parameterfis introduced. The modified extension ratio is related tofand extension ratio. The results at minimum mean misfit for all the nine rocks indicate that the range off varies from 0.7 to 1.0. It is found that mean misfit by considering varying extension ratio is lower than that in the first analysis. But it requires three parameters. A statistical expression betweenfand Hoek-Brown strength parameters has been established. Though coefficient of correlation is not reasonable, we may eliminate it as an extra parameter. At the end of the paper, a methodology has also been given for its application to isotropic jointed rock mass, so that it can be implemented in a numerical code for stability analysis of jointed rock mass structures.

A simplified approach to directly consider intact rock anisotropy in Hoek——Brown failure criterion

Many rock types have naturally occurring inherent anisotropic planes, such as bedding planes, foliation,or flow structures. Such characteristic induces directional features and anisotropy in rocks＇ strength anddeformational properties. The HoekeBrown （HeB） failure criterion is an empirical strength criterionwidely applied to rock mechanics and engineering. A direct modification to HeB failure criterion toaccount for rock anisotropy is considered as the base of the research. Such modification introduced a newdefinition of the anisotropy as direct parameter named the anisotropic parameter （Kb）. However, thecomputation of this parameter takes much experimental work and cannot be calculated in a simple way.The aim of this paper is to study the trend of the relation between the degree of anisotropy （Rc） and theminimum value of anisotropic parameter （Kmin）, and to predict the Kmin directly from the uniaxialcompression tests instead of triaxial tests, and also to decrease the amount of experimental work. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Retrospective and prospective review of the generalized nonlinear strength theory for geomaterials

Strength theory is the basic theory for calculating and designing the strength of engineering materials in civil,hydraulic,mechanical,aerospace,military,and other engineering disciplines.Therefore,the comprehensive study of the generalized nonlinear strength theory(GNST)of geomaterials has significance for the construction of engineering rock strength.This paper reviews the GNST of geomaterials to demonstrate the research status of nonlinear strength characteristics of geomaterials under complex stress paths.First,it systematically summarizes the research progress of GNST(classical and empirical criteria).Then,the latest research the authors conducted over the past five years on the GNST is introduced,and a generalized three-dimensional(3D)nonlinear Hoek‒Brown(HB)criterion(NGHB criterion)is proposed for practical applications.This criterion can be degenerated into the existing three modified HB criteria and has a better prediction performance.The strength prediction errors for six rocks and two in-situ rock masses are 2.0724%-3.5091%and 1.0144%-3.2321%,respectively.Finally,the development and outlook of the GNST are expounded,and a new topic about the building strength index of rock mass and determining the strength of in-situ engineering rock mass is proposed.The summarization of the GNST provides theoretical traceability and optimization for constructing in-situ engineering rock mass strength.

THREE-DIMENSIONAL SLOPE STABILITY ANALYSIS BASED ON NONLINEAR FAILURE ENVELOPE

The effects of nonlinearity of strength envelopes on 3D slope stability analysis are investigated.A power relation for the nonlinear envelope is employed to derive the 3D factor of safety equations of an extended Spencer method hich satisfies boty force equilibrium and moment equilibrium.Then,a search procedure is presented based on dynamic programming to determine the 3D critical slip surface for a general slope,Linear and nonlinear strength envelopes used for slope stability computations are obtained by fitting curves to the 103 strength data of consolidated-undrained(CU)triaxial compression tests for compacted Israeli clay.Results of a typical 3D problem show that a linear approximation of the nonlinear strength envelope may lead to a significant overestimation of calculated safety factors.

Calibration of Four Nonlinear Failure Envelopes from Triaxial Test Data and Influence of Nonlinearity on Geotechnical Computations

It is now recognized that many geomaterials have nonlinear failure envelopes. This non-linearity is most marked at lower stress levels, the failure envelope being of quasi-parabolic shape. It is not easy to calibrate these nonlinear failure envelopes from triaxial test data. Currently only the power-type failure envelope has been studied with an established formal procedure for its determination from triaxial test data. In this paper, a simplified procedure is evolved for the development of four different types of nonlinear envelopes. These are of invaluable assistance in the evaluation of true factors of safety in problems of slope stability and correct computation of lateral earth pressure and bearing capacity. The use of the Mohr-Coulomb failure envelopes leads to an overestimation of the factors of safety and other geotechnical quantities.

Slope reliability and back analysis of failure with geotechnical parameters estimated using Bayesian inference

A Bayesian approach is proposed for the inference of the geotechnical parameters used in slope design.The methodology involves the construction of posterior probability distributions that combine prior information on the parameter values with typical data from laboratory tests and site investigations used in design.The posterior distributions are often complex,multidimensional functions whose analysis requires the use of Markov chain Monte Carlo(MCMC)methods.These procedures are used to draw representative samples of the parameters investigated,providing information on their best estimate values,variability and correlations.The paper describes the methodology to define the posterior distributions of the input parameters for slope design and the use of these results for evaluation of the reliability of a slope with the first order reliability method(FORM).The reliability analysis corresponds to a forward stability analysis of the slope where the factor of safety(FS)is calculated with a surrogate model from the more likely values of the input parameters.The Bayesian model is also used to update the estimation of the input parameters based on the back analysis of slope failure.In this case,the condition FS?1 is treated as a data point that is compared with the model prediction of FS.The analysis requires a sufficient number of observations of failure to outbalance the effect of the initial input parameters.The parameters are updated according to their uncertainty,which is determined by the amount of data supporting them.The methodology is illustrated with an example of a rock slope characterised with a Hoek-Brown rock mass strength.The example is used to highlight the advantages of using Bayesian methods for the slope reliability analysis and to show the effects of data support on the results of the updating process from back analysis of failure.

Pseudo-static stability analysis of rock slopes reinforced by passive bolts using the generalized Hoek-Brown criterion

The stability analysis of passive bolt-reinforced rock slopes under seismic loads is investigated within the framework of the kinematic approach of limit analysis theory.A pseudo-static method is adopted to account for the inertial forces induced in the rock mass by seismic events.The strength properties of the rock material are described by a modified Hoek-Brown strength criterion,whereas the passive bolts are modeled as bar-like inclusions that exhibit only resistance to tensile-compressive forces.Taking advantage of the ability to compute closed-form expressions for the support functions associated with the modified Hoek-Brown strength criterion,a rotational failure mechanism is implemented to derive rigorous lower bound estimates for the amount of reinforcement strength to prevent slope failure.The approach is then applied to investigating the effects of relevant geometry,strength and loading parameters in light of a preliminary parametric study.The accuracy of the approach is assessed by comparison of the lower bound estimates with finite element limit analysis solutions,thus emphasizing the ability of the approach to properly predict the stability conditions and to capture the essential features of deformation localization pattern.Finally,the extension of the approach to account for slipping at the interface between reinforcements and surrounding rock mass is outlined.

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 s2.Many evidences have demonstrated that the rock strength is dependent on s2.Thus it is necessary to extend the HB criterion into a three-dimensional(3D)form.In this study,the effect of s2 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.

考虑水力效应的裂缝边坡稳定性非线性能耗分析

现有考虑水力效应的裂缝边坡稳定性分析大多基于线性破坏准则,而岩土体破坏往往呈现非线性特征,因此开展水力效应影响下的裂缝边坡稳定性非线性极限上限分析具有重要意义。基于极限分析上限定理及强度折减技术,结合“外切线法”非线性破坏准则,构建坡顶含竖直裂缝的边坡对数螺旋线破坏模式,根据虚功原理推导出裂缝边坡安全系数解析式,通过MATLAB优化计算,结合边坡工程实例探讨了典型因素对裂缝边坡稳定性、临界裂缝及滑动面位置的影响规律。研究结果表明:随着地下水位h的持续上升,边坡安全系数不断降低,临界裂缝深度逐渐增大且临界裂缝位置逐渐向坡顶缘偏移;非线性系数m明显影响边坡的稳定性,边坡安全系数随着非线性系数的增大显著减小,采用线性破坏准则会高估地下水位变化对边坡稳定性的影响;随着非线性系数的增大,临界裂缝深度随之增大,临界裂缝位置距离坡顶缘越来越远,滑坡体体积逐渐增大;裂缝与坡顶缘距离l m随着非线性系数m的增大而增大(当m增大0.2时,l m增大约1 m),且随着裂缝与坡顶缘距离的增大,边坡安全系数无明显变化,临界裂缝深度先逐渐减小后趋于稳定。

非线性破坏准则下浅埋隧道的稳定性分析

隧道围岩压力值是浅埋隧道设计计算的重要参数。为进一步提高围岩压力计算结果的可靠性,以传统浅埋隧道坍塌模型为基础,提出了圆弧滑移面破坏模式,同时基于非线性Mohr-Coulomb破坏准则,运用极限分析上限法推导了地震力作用下浅埋隧道围岩压力的计算公式,验证了计算结果的可靠性,并对影响因素进行了探讨。与实际工程和已有成果的对比分析证明了理论公式具有较高的可靠性和准确度。研究还发现:围岩压力受非线性系数和初始粘聚力的影响较为显著,非线性系数越大、初始粘聚力越小,围岩压力越大;随着待定参数K(水平围岩压力与竖向围岩压力的比值)的减小,水平围岩压力减小,竖向围岩压力增大;地震力对围岩压力的影响不容忽视,竖向地震力单独作用对围岩压力的影响最大,水平和竖向地震力共同作用的影响次之,而仅有水平地震力作用的影响最小。随着水平与竖向地震力系数的增大,围岩压力也在增大。分析结果对于大部分的隧道,尤其是易受地震影响的浅埋隧道研究具有较高的参考价值。

基于Hoek-Brown非线性破坏准则的软岩地基极限承载力数值模拟

基于有限差分分析(FLAC)的软质岩石地基承载力数值计算模型,分析中将软质岩石地基视为均质、各向同性的连续介质,采用理想线性弹塑性本构关系模拟和Hoek-Brown非线性破坏准则及非关联流动法则。为便于数值模拟,采用切线法,把Hoek-Brown非线性破坏准则转换成Mohr-Coulomb线性破坏准则,并依据等量原则推导了瞬时内摩擦角和瞬时黏聚力的计算公式,编制了基于FLAC内置FISH语言的子程序。通过施加竖向速度向量,分别获得了圆形基础下考虑和忽略岩体自重的极限承载力。与理论值和现场载荷板承载力试验结果进行的对比分析表明,Hoek-Brown解结果比Bell解结果更接近实测值,但都偏于保守,而数值计算结果与实测值非常接近,特别是在考虑自重的情况下,误差仅为2%。说明对于软质岩石地基的极限承载力,采用基于Hoek-Brown非线性破坏准则的FLAC数值模型进行分析是合适的。

基于广义Hoek-Brown非线性破坏准则的节理岩体地基承载力研究

针对节理岩体非线性破坏特征,发展基于Hoek-Brown破坏准则的临界滑动场理论进行节理岩体地基承载力的计算。首先将Hoek-Brown准则的剪切强度逐点等效到Mohr-Coulomb强度线上,求得每点的瞬时内摩擦角和瞬时黏聚力;在此基础上,改进基于Mohr-Coulomb准则的临界滑动场理论,建立新的迭代算法,将Hoek-Brown强度准则与临界滑动场理论结合起来,求解被动土压力和地基承载力。同时,综合分析GSI和mi值对地基承载力的影响。算例对比分析结果表明,该方法能迅速准确地确定节理岩体地基最危险破坏滑面并得到相应的地基极限承载力值。