A statistical damage-based constitutive model for shearing of rock joints in brittle drop mode

Some rock joints exhibit significant brittleness,characterized by a sharp decrease in shear stress upon reaching the peak strength.However,existing models often fail to accurately represent this behavior and are encumbered by numerous parameters lacking clear mechanical significance.This study presents a new statistical damage constitutive model rooted in both damage mechanics and statistics,containing only three model parameters.The proposed model encompasses all stages of joint shearing,including the compaction stage,linear stage,plastic yielding stage,drop stage,strain softening stage,and residual strength stage.To derive the analytical expression of the constitutive model,three boundary conditions are introduced.Experimental data from both natural and artificial rock joints is utilized to validate the model,resulting in average absolute relative errors ranging from 3%to 8%.Moreover,a comparative analysis with established models illustrates that the proposed model captures stress drop and post-peak strain softening more effectively,with model parameters possessing clearer mechanical interpretations.Furthermore,parameter analysis is conducted to investigate the impacts of model parameters on the curves and unveil the relationship between these parameters and the mechanical properties of rock joints.Importantly,the proposed model is straightforward in form,and all model parameters can be obtained from direct shear tests,thus facilitating the utilization in numerical simulations.

Heterogeneities of grain boundary contact for simulation of laboratoryscale mechanical behavior of granitic rocks

From a practical point of view,grain structure heterogeneities are key parameters that control the rock response and still remains a challenge to incorporate in a quantitative manner.One of the less discussed topics in the context of the grain-based model(GBM)in the particle flow code(PFC)is the contact heterogeneities and the appropriate contact model to mimic the grain boundary behavior.Generally,the smooth joint(SJ)model and linear parallel bond(LPB)model are used to simulate the grain boundary behavior.However,the literature does not document the suitability of different models for specific problems.Another challenge in implementing GBM in PFC is that only a single bonding parameter is used at the grain boundaries.The aim of this study is to investigate the responses of a laboratory-scale specimen with SJ and LPB models,considering grain boundary heterogeneous and homogeneous contact parameters.Uniaxial and biaxial compression tests are performed to calibrate the response of Creighton granite.The stressestrain curves,volumetric dilation,inter-crack(crack in the grain boundary),and intra-crack(crack within the grain)development,and failure patterns associated with different contact models are examined.It was found that both the SJ and LPB models can reproduce the pre-peak behavior observed for a granitic rock type.However,the LPB model is unable to reproduce the post-peak behavior.Due to the large interlocking effect originating from the balls in contact and the ball size in the LPB model,local dilation is induced at the grain boundaries.This overestimates the volumetric dilation and residual shear strength.The LPB model tends to result in discontinuous inter-cracks and stress localization in the rock specimen,resulting in fine fragments at the rock surface during failure.

Strength and deformation characteristics of irregular columnar jointed rock mass: A combined experimental and theoretical study

The irregularity of jointed network poses a challenge to the determination of field mechanical param-eters of columnar jointed rock mass(CJRM),and a reasonable prediction of deformation and strength characteristics of CJRM is important for engineering construction.The Voronoi diagram and three-dimensional printing technology were used to make an irregular columnar jointed mold,and the irregular CJRM(ICJRM)specimens with different dip directions and dip angles were prepared.Uniaxial compression tests were performed,and the anisotropic strength and deformation characteristics of ICJRM were described.The failure modes and mechanisms were revealed in accordance with the final appearances of the ICJRM specimens.Based on the model test results,the empirical correlations for determining the field deformation and strength parameters of CJRM were derived using the dip angle and modified joint factor.The proposed empirical equations were used in the Baihetan Project,and the calculated mechanical parameters were compared with the field test results and those obtained from the tunneling quality index method.Results showed that the deformation parameters determined by the two proposed methods are all consistent with the field test results,and these two methods can also estimate the strength parameters effectively.

A method to predict the peak shear strength of rock joints based on machine learning

In geotechnical and tunneling engineering,accurately determining the mechanical properties of jointed rock holds great significance for project safety assessments.Peak shear strength(PSS),being the paramount mechanical property of joints,has been a focal point in the research field.There are limitations in the current peak shear strength(PSS)prediction models for jointed rock:(i)the models do not comprehensively consider various influencing factors,and a PSS prediction model covering seven factors has not been established,including the sampling interval of the joints,the surface roughness of the joints,the normal stress,the basic friction angle,the uniaxial tensile strength,the uniaxial compressive strength,and the joint size for coupled joints;(ii)the datasets used to train the models are relatively limited;and(iii)there is a controversy regarding whether compressive or tensile strength should be used as the strength term among the influencing factors.To overcome these limitations,we developed four machine learning models covering these seven influencing factors,three relying on Support Vector Regression(SVR)with different kernel functions(linear,polynomial,and Radial Basis Function(RBF))and one using deep learning(DL).Based on these seven influencing factors,we compiled a dataset comprising the outcomes of 493 published direct shear tests for the training and validation of these four models.We compared the prediction performance of these four machine learning models with Tang’s and Tatone’s models.The prediction errors of Tang’s and Tatone’s models are 21.8%and 17.7%,respectively,while SVR_linear is at 16.6%,SVR_poly is at 14.0%,and SVR_RBF is at 12.1%.DL outperforms the two existing models with only an 8.5%error.Additionally,we performed shear tests on granite joints to validate the predictive capability of the DL-based model.With the DL approach,the results suggest that uniaxial tensile strength is recommended as the material strength term in the PSS model for more reliable outcomes.

Nonlinear shear behavior of rock joints using a linearized implementation of the Barton-Bandis model

Experiments on rock joint behaviors have shown that joint surface roughness is mobilized under shearing,inducing dilation and resulting in nonlinear joint shear strength and shear stress vs.shear displacement behaviors.The Barton-Bandis(B-B) joint model provides the most realistic prediction for the nonlinear shear behavior of rock joints.The B-B model accounts for asperity roughness and strength through the joint roughness coefficient(JRC) and joint wall compressive strength(JCS) parameters.Nevertheless,many computer codes for rock engineering analysis still use the constant shear strength parameters from the linear Mohr-Coulomb(M-C) model,which is only appropriate for smooth and non-dilatant joints.This limitation prevents fractured rock models from capturing the nonlinearity of joint shear behavior.To bridge the B-B and the M C models,this paper aims to provide a linearized implementation of the B-B model using a tangential technique to obtain the equivalent M-C parameters that can satisfy the nonlinear shear behavior of rock joints.These equivalent parameters,namely the equivalent peak cohesion,friction angle,and dilation angle,are then converted into their mobilized forms to account for the mobilization and degradation of JRC under shearing.The conversion is done by expressing JRC in the equivalent peak parameters as functions of joint shear displacement using proposed hyperbolic and logarithmic functions at the pre-and post-peak regions of shear displacement,respectively.Likewise,the pre-and post-peak joint shear stiffnesses are derived so that a complete shear stress-shear displacement relationship can be established.Verifications of the linearized implementation of the B-B model show that the shear stress-shear displacement curves,the dilation behavior,and the shear strength envelopes of rock joints are consistent with available experimental and numerical results.

Discrete modeling of rock joints with a smooth-joint contact model

Structural defects such as joints or faults are inherent to almost any rock mass.In many situations those defects have a major impact on slope stability as they can control the possible failure mechanisms.Having a good estimate of their strength then becomes crucial.The roughness of a structure is a major contributor to its strength through two different aspects,i.e.the morphology of the surface(or the shape)and the strength of the asperities(related to the strength of the rock).In the current state of practice,roughness is assessed through idealized descriptions(Patton strength criterion)or through empirical parameters(Barton JRC).In both cases,the multi-dimensionality of the roughness is ignored.In this study,we propose to take advantage of the latest developments in numerical techniques.With3D photogrammetry and/or laser mapping,practitioners have access to the real morphology of an exposed structure.The derived triangulated surface was introduced into the DEM(discrete element method)code PFC3D to create a synthetic rock joint.The interaction between particles on either side of the discontinuity was described by a smooth-joint model(SJM),hence suppressing the artificial roughness introduced by the particle discretization.Shear tests were then performed on the synthetic rock joint.A good correspondence between strengths predicted by the model and strengths derived from well-established techniques was obtained for thefirst time.Amongst the benefits of the methodology is the possibility offered by the model to be used in a quantitative way for shear strength estimates,to reproduce the progressive degradation of the asperities upon shearing and to analyze structures of different scales without introducing any empirical relation.

A five-parameter constitutive model for hysteresis shearing and energy dissipation of rock joints

Rock joints exhibit hysteresis shearing behavior and produce energy dissipation under shear cyclic loads,which however cannot be accurately depicted by existing constitutive models. This paper establishes a constitutive model for hysteresis shearing and associated energy dissipation of rock joints. Analytical expressions of the model during cyclic shearing processes are derived. Derivation of the model indicates no energy dissipation in the elastic stage. When the shear load exceeds elastic boundary, nonlinear energy dissipation takes place. Validations with experiments show that the proposed model provides good conformities with direct shear curves and hysteresis loops, and can predict the energy dissipation characteristics of rock joints under different working conditions. Compared to the constitutive models using Weibull's distribution, the proposed one is smooth at the elastic boundary and can accurately capture the maximum shear stress. Unlike the existing incremental-type models, the proposed one provides clear and direct analytical expressions for both shear stress and energy dissipation during the whole displacement domain, which is more convenient in application.

Analytical model of shear mechanical behaviour of bolted rock joints considering influence of normal stress on bolt guide rail effect

Rock bolts have been widely used in slopes as a reinforcement measure.Modelling the shear mechanical behaviours of bolted rock joints is very complicated due to the complex factors that affect the axial force and shear force on the bolts.Rock bolts under shear action exhibit the guide rail effect;that is,the rock mass slides along the rock bolt as if the rock bolt is a rail.The normal stress can inhibit the guide rail effect and reduce the axial force on bolts.However,this factor is not considered by the existing analysis models.Shear tests of bolted joints under different normal stresses were carried out in the laboratory.During the test,the axial force on each point monitored on the bolt was recorded by a strain gauge,and the attenuation trend of the strain was studied.An analytical model that considers the inhibition of the bolt rail effect due to an increase in the normal stress was proposed to predict the shear mechanical behaviour of rock bolted joints.The new model accommodates the bolt shear behaviours in the elastic stage and plastic stage,and the estimated values agree well with the results of the direct shear tests in the laboratory.The validation shows that the proposed model can effectively describe the deformation characteristics of the bolts in the shear tests.

An Anisotropic Geomechanical Model for Jointed Rock Masses at Taishan Nuclear Power Station, Southern China

An anisotropic geomechanical model for jointed rock mass is presented. Simultaneously with deriving the orthotropic anisotropy elastic parameters along the positive axis, the equivalent compliance matrix for the deflection axis orthotropic anisotropy was derived through a three- dimensional coordinate transformation. In addition, Singh＇s analysis of the stress concentration effects of intermittent joints was adopted, based on two groups of intermittent joints and a set of cross- cutting joints in the jointed rock mass. The stress concentration effects caused by intermittent joints and the coupling effect of cross-cutting joints along the deflection-axis are also considered. The proposed anisotropic mechanics parameters method is applied to determine the deformation parameters of jointed granite at the Taishan Nuclear Power Station. Combined with the deterministic mechanical parameters of rock blocks and joints, the deformation parameters and their variability in jointed rock masses are estimated quantitatively. The computed results show that jointed granite at the Taishan Nuclear Power Station exhibits typical anisotropic mechanical characteristics; the elastic moduli in the two horizontal directions were similar, but the elastic modulus in the vertical direction was much greater. Jointed rock elastic moduli in the two horizontal and vertical directions were respectively about 24% and 37% of the core of rock, showing weakly orthotropic anisotropy; the ratio of elastic moduli in the vertical and horizontal directions was 1.53, clearly indicating the transversely isotropic rock mass mechanical characteristics. The method can be popularized to solve other rock mechanics problems in nuclear power engineering.

Distinct element modelling of fracture plan control in continuum and jointed rock mass in presplitting method of surface mining

Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surface blast design. The purpose of presplitting is to form a fracture plane across which the radial cracks from the production blast cannot travel. The purpose of this study is to investigate of effect of presplitting on the generation of a smooth wall in continuum and jointed rock mass. The 2D distinct element code was used to simulate the presplitting in a rock slope. The blast load history as a function of time was applied to the inner wall of each blasthole. Important parameters that were considered in the analysis were stress tensor and fracturing pattern. The blast loading magnitude and blasthole spacing and jointing pattern were found to be very significant in the final results.

AN ANALYSIS MODEL OF THE MECHANICS OF JOINTED ROCK MASS

The creep fracture damage model of jointed rock mass is proposed, the basic methods and properties of engineering structure forecast analysis in split surrounding rock mass is discussed, numericalized way is analyzed, and the computing program of RCFD rock mechanics problems known is modified.Moreover, this program is used to make practical forecast analysis for the stability of surrounding rock mass of one large section (vertical)shaft which is buried in the earth 300 m deep.

Crack initiation stress and strain of jointed rock containing multi-cracks under uniaxial compressive loading: A particle flow code approach

The ratio of crack initiation stress to the uniaxial compressive strength(SCI,B/SUC,B) and the ratio of axial strain at the crack initiation stress to the axial strain at the uniaxial compressive strength(B,UCB,CI,A,A/SSSS) were studied by performing numerical stress analysis on blocks having multi flaws at close spacing's under uniaxial loading using PFC3 D. The following findings are obtained: SCI,B/SUC,B has an average value of about 0.5 with a variability of ± 0.1. This range agrees quite well with the values obtained by former research. For joint inclination angle, β=90°,B,UCB,CI,A,A/SSSS is found to be around 0.48 irrespective of the value of joint continuity factor, k. No particular relation is found betweenB,UCB,CI,A,A/SSSS and β; however, the average B,UCB,CI,A,A/SSSS seems to slightly decrease with increasing k. The variability ofB,UCB,CI,A,A/SSSS is found to increase with k.Based on the cases studied in this work,B,UCB,CI,A,A/SSSS ranges between 0.3 and 0.5. This range is quite close to the range of 0.4to 0.6 obtained for SCI,B/SUC,B. The highest variability of ± 0.12 forB,UCB,CI,A,A/SSSS is obtained for k=0.8. For the remaining k values the variability ofB,UCB,CI,A,A/SSSS can be expressed within ± 0.05. This finding is very similar to the finding obtained for the variability of SCI,B/SUC,B.

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（MC） 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 σn tan φ＇ 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.

Stresses induced by post-tensioned anchor in jointed rock mass

A new analytical study on stresses around a post-tensioned anchor in rocks with two perpendicular joint sets is presented. The assumptions of orthotropic elastic rock with plane strain conditions are made in derivation of the formulations. A tri-linear bond-slip constitutive law is used for modeling the tendon-grout interface behavior and debonding of this interface. The bearing plate width is also considered in the analysis. The obtained solutions are in the integral forms and numerical techniques that have been used for evaluation. In the illustrative example given, the major principal stress is compressive in the anchor free zone and compressive stress concentrations of 815 k Pa and 727 k Pa(for the anchor load of 300 k N) are observed under the bearing plate and the bond length proximal end, respectively. However, large values of tensile stresses with the maximum of-434 k Pa are formed at the bond length distal end. The results obtained using the proposed solution are compared very those of numerical method(FEM).

An extended displacement discontinuity method for analysis of stress wave propagation in viscoelastic rock mass

An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primary group with an average geometrical size larger than or in the same order of magnitude of wavelength of a concerned stress wave is defined as 'macro-joints',while the secondary group with a high density and relatively small geometrical size compared to the wavelength is known as 'micro-defects'.The rock mass with micro-defects is modeled as an equivalent viscoelastic medium while the macro-joints in the rock mass are modeled explicitly as physical discontinuities.Viscoelastic properties of a micro-defected sedimentary rock are obtained by longitudinally impacting a cored long sedimentary rod with a pendulum.Wave propagation coefficient and dynamic viscoelastic modulus are measured.The EDDM is then successfully employed to analyze the wave propagation across macro-joint in VRM.The effect of the rock viscosity on the stress wave propagation is evaluated by comparing the results of VRM from the presented EDDM with those of an elastic rock mass (ERM) from the conventional displacement discontinuity method (CDDM).The CDDM is a special case of the EDDM under the condition that the rock viscosity is ignored.Comparison of the reflected and transmitted waves shows that the essential rock viscosity has a significant effect on stress wave attenuation.When a short propagation distance of a stress wave is considered,the results obtained from the CDDM approximate to the EDDM solutions,however,when the propagation distance is sufficiently long relative to the wavelength,the effect of rock viscosity on the stress wave propagation cannot be ignored.

Determination of mechanical parameters for jointed rock masses

Combining with empirical method, laboratory test and numerical simulation, a comprehensive system was presented to determine the mechanical parameters of jointed rock masses. The system has the following four functions: (1) Based on the field investigation of joints, the system can consider rock mass structures, by using network simulation technology. (2) Rock samples are conducted by numerical simulation with the input engineering mechanical parameters of rocks and joints obtained from laboratory tests. (3) The whole stress-strain curve of jointed rock masses under certain normal stress can be plotted from numerical simulation, and then the shear strength parameters of jointed rock masses can be obtained from the whole stress-strain curves under different normal stresses. (4) The statistical values of mechanical parameters of jointed rock masses can be determined according to numerical simulation. Based on the statistical values, combining with engineering experiences and geological investigations, the comprehensive mechanical parameters of jointed rock masses can be achieved finally. Several cases are presented to prove the engineering feasibility and suitability of this system.

Multi-scale data joint inversion of minerals and porosity in altered igneous reservoirs—A case study in the South China Sea

There are abundant igneous gas reservoirs in the South China Sea with significant value of research,and lithology classification,mineral analysis and porosity inversion are important links in reservoir evaluation.However,affected by the diverse lithology,complicated mineral and widespread alteration,conventional logging lithology classification and mineral inversion become considerably difficult.At the same time,owing to the limitation of the wireline log response equation,the quantity and accuracy of minerals can hardly meet the exploration requirements of igneous formations.To overcome those issues,this study takes the South China Sea as an example,and combines multi-scale data such as micro rock slices,petrophysical experiments,wireline log and element cutting log to establish a set of joint inversion methods for minerals and porosity of altered igneous rocks.Specifically,we define the lithology and mineral characteristics through core slices and mineral data,and establish an igneous multi-mineral volumetric model.Then we determine element cutting log correction method based on core element data,and combine wireline log and corrected element cutting log to perform the lithology classification and joint inversion of minerals and porosity.However,it is always difficult to determine the elemental eigenvalues of different minerals in inversion.This paper uses multiple linear regression methods to solve this problem.Finally,an integrated inversion technique for altered igneous formations was developed.The results show that the corrected element cutting log are in good agreement with the core element data,and the mineral and porosity results obtained from the joint inversion based on the wireline log and corrected element cutting log are also in good agreement with the core data from X-ray diffraction.The results demonstrate that the inversion technique is applicable and this study provides a new direction for the mineral inversion research of altered igneous formations.

A CUSPC ATASTROPHIC MODEL AND TIME EFFECT OF ROCK BURSTS OF A

Ａｍｅｃｈａｎｉｃａｌｍｏｄｅｌｆｏｒｔｈｅｒｏｃｋｂｕｒｓｔｓｏｆａｎｉｓｏｌａｔｅｄｃｏａｌｐｉｌｌａｒｗａｓｄｅｖｅｌｏｐｅｄａｎｄｔｈｅｍｅｃｈａｎｉｓｍｏｆｔｈｅｉｎｓｔａｂｉｌ－ｉｔｙｏｆｔｈｅｒｏｃｋｂｕｒｓｔｓｗａｓｓｔｕｄｉｅｄｂｙａｐｐｌｙｉｎｇｃｕｓｐｃａｔａｓｔｒｏｐｈｉｃｔｈｅｏｒｙ．Ｉｔｗａｓｅｘｐｏｕｎｄｅｄｔｈａｔｔｈｅｏｃｃｕｒｒｅｎｃｅｓｃｏｎｄｉ－ｔｉｏｎｓｏｆｒｏｃｋｂｕｒｓｔｓｃａｎｂｅａｔｔｒｉｂｕｔｅｄｔｏｔｈｅｃｏｎｄｉｔｉｏｎｔｈａｔｔｈｅｓｙｓｔｅｍａｃｃｅｐｔｓｅｎｏｕｇｈｅｎｅｒｇｙｆｒｏｍｅｘｔｅｒｎａｌｅｎｖｉｒｏｎ－ｍｅｎｔａｎｄｔｈｅｃｏｎｄｉｔｉｏｎｔｈａｔｔｈｅｅｎｅｒｇｙａｃｃｕｍｕｌａｔｅｄｉｎｓｙｓｔｅｍｉｓｒｅｌｅａｓｅｄｓｕｄｄｅｎｌｙ．Ｉｔｉｓｔｈｅｃａｕｓｅｏｆｔｈｅｒｏｃｋｂｕｒｓｔｓｈｙｓｔｅｒｅｓｉｓｔｈａｔｃｏｎｔｒｏｌｖａｒｉａｂｌｅｓｄｏｎ’ｔｓａｔｉｓｆｙｔｈｅｅｑｕａｔｉｏｎｏｆｂｉｆｕｒｃａｔｉｏｎｓｅｔｕｎｌｅｓｓａｐｉｅｃｅｏｆｐｅｒｉｏｄｈａｓｅｌａｐｓｅｄ．Ｆｉｎａｌ－ｌｙ，ｔｈｅｉｎｆｌｕｅｎｃｅｏ?

Joints in unsaturated rocks:Thermo-hydro-mechanical formulation and constitutive behaviour

A formulation for the coupled analysis of thermo-hydro-mechanical （THM） problems in joints is first presented. The work involves the establishment of equilibrium and mass and energy balance equations. Balance equations were formulated taking into account two phases： water and air. The joint element developed was implemented in a general purpose finite element computer code for THM analysis of porous media （Code_Bright）. The program was then used to study a number of cases ranging from laboratory tests to large scale in situ tests. A numerical simulation of coupled hydraulic shear tests of rough granite joints is first presented. The tests as well as the model show the coupling between permeability and the deformation of thejoints. The experimental investigation was focused on the effects of suction on the mechanical behaviour of rock joints. Laboratory tests were performed in a direct shear cell equipped with suction control. Suction was imposed using a vapour forced convection circuit connected to the cell and controlled by an air pump. Artificial joints of Lilla claystone were prepared.Joint roughness of varying intensity was created by carving the surfaces in contact in such a manner that rock ridges of different tip angles were formed. These angles ranged from 0° （smooth joint） to 45° （very rough joint profile）. The geometric profiles of the two surfaces in contact were initially positioned in a ＂matching＂ situation. Several tests were performed for different values of suctions （200, 100, and 20 MPa） and for different values of vertical stresses （30, 60, and 150 kPa）. A constitutive model including the effects of suction and joint roughness is proposed to simulate the unsaturated behaviour of rock joints. The new constitutive law was incorporated in the code and experimental results were numerically simulated.

Advances in statistical mechanics of rock masses and its engineering applications

To efficiently link the continuum mechanics for rocks with the structural statistics of rock masses,a theoretical and methodological system called the statistical mechanics of rock masses(SMRM)was developed in the past three decades.In SMRM,equivalent continuum models of stressestrain relationship,strength and failure probability for jointed rock masses were established,which were based on the geometric probability models characterising the rock mass structure.This follows the statistical physics,the continuum mechanics,the fracture mechanics and the weakest link hypothesis.A general constitutive model and complete stressestrain models under compressive and shear conditions were also developed as the derivatives of the SMRM theory.An SMRM calculation system was then developed to provide fast and precise solutions for parameter estimations of rock masses,such as full-direction rock quality designation(RQD),elastic modulus,Coulomb compressive strength,rock mass quality rating,and Poisson’s ratio and shear strength.The constitutive equations involved in SMRM were integrated into a FLAC3D based numerical module to apply for engineering rock masses.It is also capable of analysing the complete deformation of rock masses and active reinforcement of engineering rock masses.Examples of engineering applications of SMRM were presented,including a rock mass at QBT hydropower station in northwestern China,a dam slope of Zongo II hydropower station in D.R.Congo,an open-pit mine in Dexing,China,an underground powerhouse of Jinping I hydropower station in southwestern China,and a typical circular tunnel in Lanzhou-Chongqing railway,China.These applications verified the reliability of the SMRM and demonstrated its applicability to broad engineering issues associated with jointed rock masses.