On the calibration of a shear stress criterion for rock joints to represent the full stress-strain profile

Conventional numerical solutions developed to describe the geomechanical behavior of rock interfaces subjected to differential load emphasize peak and residual shear strengths.The detailed analysis of preand post-peak shear stress-displacement behavior is central to various time-dependent and dynamic rock mechanic problems such as rockbursts and structural instabilities in highly stressed conditions.The complete stress-displacement surface(CSDS)model was developed to describe analytically the pre-and post-peak behavior of rock interfaces under differential loads.Original formulations of the CSDS model required extensive curve-fitting iterations which limited its practical applicability and transparent integration into engineering tools.The present work proposes modifications to the CSDS model aimed at developing a comprehensive and modern calibration protocol to describe the complete shear stressdisplacement behavior of rock interfaces under differential loads.The proposed update to the CSDS model incorporates the concept of mobilized shear strength to enhance the post-peak formulations.Barton’s concepts of joint roughness coefficient(JRC)and joint compressive strength(JCS)are incorporated to facilitate empirical estimations for peak shear stress and normal closure relations.Triaxial/uniaxial compression test and direct shear test results are used to validate the updated model and exemplify the proposed calibration method.The results illustrate that the revised model successfully predicts the post-peak and complete axial stressestrain and shear stressedisplacement curves for rock joints.

Effect of roughness on the shear behavior of rock joints subjected to impact loading

The shear behavior is regarded as the dominant property of rock joints and is dramatically affected by the joint surface roughness.To date,the effect of surface roughness on the shear behavior of rock joints under static or cyclic loading conditions has been extensively studied,but such effect under impact loading conditions keeps unclear.To address this issue,a series of impact shear tests was performed using a novel-designed dynamic experimental system combined with the digital image correlation(DIC)technique.The dynamic shear strength,deformability and failure mode of the jointed specimens with various joint roughness coefficients(JRC)are comprehensively analyzed.Results show that the shear strength and shear displacement characteristics of the rock joint under the impact loading keep consistent with those under static loading conditions.However,the temporal variations of shear stress,slip displacement and normal displacement under the impact loading conditions show obviously different behaviors.An elastic rebound of the slip displacement occurs during the impact shearing and its value increases with increasing joint roughness.Two identifiable stages(i.e.compression and dilation)are observed in the normal displacement curves for the rougher rock joints,whereas the joints with small roughness only manifest normal compression displacement.Besides,as the roughness increases,the maximum compression tends to decrease,while the maximum dilation gradually increases.More-over,the microstructural analysis based on scanning electron microscope(SEM)suggests that the roughness significantly affects the characteristics of the shear fractured zone enclosing the joint surface.

Experimental study on the shear performance of quasi-NPR steel bolted rock joints

Quasi-NPR(negative Poisson’s ratio)steel is a new type of super bolt material with high strength,high ductility,and a micro-negative Poisson’s effect.This material overcomes the contrasting characteristics of the high strength and high ductility of steel and it has significant energy-absorbing characteristics,which is of high value in deep rock and soil support engineering.However,research on the shear resistance of quasi-NPR steel has not been carried out.To study the shear performance of quasi-NPR steel bolted rock joints,indoor shear tests of bolted rock joints under different normal stress conditions were carried out.Q235 steel and#45 steel,two representative ordinary bolt steels,were set up as a control group for comparative tests to compare and analyze the shear strength,deformation and instability mode,shear energy absorption characteristics,and bolting contribution of different types of bolts.The results show that the jointed rock masses without bolt reinforcement undergo brittle failure under shear load,while the bolted jointed rock masses show obvious ductile failure characteristics.The shear deformation ca-pacity of quasi-NPR steel is more than 3.5 times that of Q235 steel and#45 steel.No fracture occurs in the quasi-NPR steel during large shear deformation and it can provide stable shear resistance.However,the other two types of control bolts become fractured under the same conditions.Quasi-NPR steel has significant energy-absorbing characteristics under shear load and has obvious advantages in terms of absorbing the energy released by shear deformation of jointed rock masses as compared with ordinary steel.In particular,the shear force plays a major role in resisting the shear deformation of Q235 steel and#45 steel,therefore,fracture failure occurs under small bolt deformation.However,the axial force of quasi-NPR steel can be fully exerted when resisting joint shear deformation;the steel itself does not break when large shear deformation occurs,and the supporting effect of the jointed rock mass is effectively guaranteed.

Investigation of long-wavelength elastic wave propagation through wet bentonite-filled rock joints

The saturation of the compacted bentonite buffer in the deep geological repository can cause bentonite swelling,intrusion into rock fractures,and erosion.Inevitably,erosion and subsequent bentonite mass loss due to groundwater inflow can aggravate the overall integrity of the engineered barrier system.Therefore,the coupled hydro-mechanical interaction between the buffer and rock during groundwater inflow and bentonite intrusion should be evaluated to guarantee the long-term safety of deep geological disposal.This study investigated the effect of bentonite erosion and intrusion on the elastic wave propagation characteristics in jointed rocks using a quasi-static resonant column test.Jointed rock specimens with different joint conditions(i.e.joint surface saturation and bentonite filling)were prepared using granite rock discs sampled from the Korea Underground Research Tunnel(KURT)and Gyeongju bentonite.The long-wavelength longitudinal and shear wave velocities were measured under different normal stress levels.A Hertzian-type power model was used to fit the wave velocities,and the relationship between the two fitted parameters provided the trend of joint conditions.Numerical simulations using three-dimensional distinct element code(3DEC)were conducted to better understand how the long-wavelength wave propagates through wet bentonite-filled rock joints.

Shear behavior of two-order morphology in rock joints

Two-order morphology of rock joints named as waviness and unevenness can be separated by morphology classification method,which plays a decisive role in the evolution of shear stress during the shear test.The joint morphology is obtained by using 3D printing and 3D laser scanning techniques and the joint model samples in two-order morphology are produced by cement mortar.Then,shear tests are performed under different normal loads.Results shows that the waviness is dominant in the total morphology during the shear test,and the shear contribution of unevenness mainly occurs in the climbing phase of shearing process.Comparing the failure modes of two-order morphology,waviness mainly embodies shear dilation characteristics and unevenness mainly shows shear wear characteristics.Based on this,a quantitative parameter is proposed to represent the ratio of the peak shear strength of the two-order morphology to that of total morphology.The functional relationship between the peak shear strength of total and two-order morphologies is determined,providing a theoretical method for further in-depth study on the shear strength of the interaction with two-order morphology of rock joints.

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.

Shear behaviour of a rock bridge sandwiched between incipient joints under the influence of hydraulic pressures

The rock bridges sandwiched in incipiently jointed rock mass were considered as barriers that block the fluid seepage,and provide certain shear strength reservation.For better revealing the influence of hydraulic pressure on the failure behaviour of rock bridges,direct shear tests were carried out through a newly proposed method on rock samples that contain two parallel incipient joints.By developing the gypsum-silicone pad coupling samples,a conventional triaxial test system was qualified to implement direct shear tests with satisfied sealing capability.The results showed that the rock bridges could be failed through the tensile failure,shear failure and mixed failure mechanism.The hydraulic pressure would facilitate the tensile failure mechanism and induce rougher fracture surfaces;while the normal stress would facilitate the shear failure mechanism and induce less rough fracture.The hydraulic pressure reduced the global shear strength of the rock block through reducing the efficient normal stress applied on the rock bridge area,which was highly dependent on the joint persistence,k.Moreover,because of the iterating occurrence of the hydraulic pressure lag with the fracture propagation,the rock bridge failure stage in the shear stress-shear displacement curves displayed a fluctuation trend.

Method for visualizing the shear process of rock joints using 3D laser scanning and 3D printing techniques

This study presents a visualized approach for tracking joint surface morphology.Three-dimensional laser scanning(3DLS)and 3D printing(3DP)techniques are adopted to record progressive failure during rock joint shearing.The 3DP resin is used to create transparent specimens to reproduce the surface morphology of a natural joint precisely.The freezing method is employed to enhance the mechanical properties of the 3DP specimens to reproduce the properties of hard rock more accurately.A video camera containing a charge-coupled device(CCD)camera is utilized to record the evolution of damaged area of joint surface during the direct shear test.The optimal shooting distance and shooting angle are recommended to be 800 mm and 40?,respectively.The images captured by the CCD camera are corrected to quantitatively describe the damaged area on the joint surface.Verification indicates that this method can accurately describe the total sheared areas at different shear stages.These findings may contribute to elucidating the shear behavior of rock joints.

Effect of rock joints on lined rock caverns subjected to high internal gaspressure

The storage of hydrogen gas in lined rock caverns(LRCs)may enable the implementation of the firstlarge-scale fossil-free steelmaking process in Sweden,but filling such storage causes joints in the rockmass to open,concentrating strains in the lining.The structural interaction between the LRC componentsmust be able to reduce the strain concentration in the sealing steel lining;however,this interaction iscomplex and difficult to predict with analytical methods.In this paper,the strain concentration in LRCsfrom the opening of rock joints is studied using finite element(FE)analyses,where the large-and small-scale deformation behaviors of the LRC are coupled.The model also includes concrete crack initiation anddevelopment with increasing gas pressure and rock joint width.The interaction between the jointed rockmass and the reinforced concrete,the sliding layer,and the steel lining is demonstrated.The results showthat the rock mass quality and the spacing of the rock joints have the greatest influence on the straindistributions in the steel lining.The largest effect of rock joints on the maximum strains in the steellining was observed for geological conditions of“good”quality rock masses.

Stability analysis of intermittently jointed rock slopes based on the stepped failure mode

In practical engineering,due to the noncontinuity characteristics of joints in rock slopes,in addition to plane failure,stepped sliding failure may occur for intermittently jointed rock slopes.Especially for intermittently bedding jointed rock slopes,the correlation and difference in strength parameters between joints and rock bridges,along with the various failure modes and intermittency of rock bridges,contribute to the complexity of stepped failure modes and the unpredictability of failure regions.Based on the upper-bound limit analysis method and multi-sliders step-path failure mode,considering the shear and tensile failure of rock bridges and the weakened relationship between the strength parameters of rock bridges and jointed surfaces,by introducing the modified M-C failure criterion and the formula for calculating the energy consumption of tensile failure of rock bridges,two failure mechanisms are constructed to obtain the safety factor(F_(s))of intermittently jointed rock slopes.The sequential quadratic programming method is used to obtain the optimal upper-bound solution for F_(s).The influence of multiple key parameters(slope height H,horizontal distance L,Slope angleβ,shear strength parameters of the rock bridgeφr and cr,Dimensionless parameter u,weakening coefficients of the internal friction angle and cohesion between the rock bridges and joint surfaces Kφand Kc)on the stability analysis of intermittently jointed rock slopes under the shear failure mode of rock bridges as well as under the tensile failure mode is also explored.The reliability of the failure mechanisms is verified by comparative analysis with theoretical results,numerical results,and landslide cases,and the variation rules of F_(s)with each key parameter are obtained.The results show that F_(s) varies linearly withφr and cr of the rock bridge and with K_(φ)and K_(c),whereas F_(s)changes nonlinearly with H and L.In particular,with the increase in Kφand Kc,Fs increases by approximately 52.78%and 171.02%on average,respectively.For rock bridge tensile failure,F_(s) shows a nonlinearly positive correlation withφr,cr,Kφand Kc.In particular,with the increase in Kφand Kc,Fs increases by approximately 13%and 61.69%on average,respectively.Fs decreases rapidly with increasing slope gradientβand decreasing dimensionless parameterμ.When Kφand Kc are both less than 1.0,the stepped sliding surface occurs more easily than the plane failure surface,especially in the case of tensile failure of the rock bridge.In addition,rock slopes with higher strength parameters,taller heights,and greater weakening coefficients are prone to rock bridge tension failure with lower Fs,and more attention should be given to the occurrence of such accidents in actual engineering.

A photogrammetric approach for quantifying the evolution of rock joint void geometry under varying contact states

Accurate measurement of the evolution of rock joint void geometry is essential for comprehending the distribution characteristics of asperities responsible for shear and seepage behaviors.However,existing techniques often require specialized equipment and skilled operators,posing practical challenges.In this study,a cost-effective photogrammetric approach is proposed.Particularly,local coordinate systems are established to facilitate the alignment and precise quantification of the relative position between two halves of a rock joint.Push/pull tests are conducted on rock joints with varying roughness levels to induce different contact states.A high-precision laser scanner serves as a benchmark for evaluating the photogrammetry method.Despite certain deviations exist,the measured evolution of void geometry is generally consistent with the qualitative findings of previous studies.The photogrammetric measurements yield comparable accuracy to laser scanning,with maximum errors of 13.2%for aperture and 14.4%for void volume.Most joint matching coefficient(JMC)measurement errors are below 20%.Larger measurement errors occur primarily in highly mismatched rock joints with JMC values below 0.2,but even in cases where measurement errors exceed 80%,the maximum JMC error is only 0.0434.Thus,the proposed photogrammetric approach holds promise for widespread application in void geometry measurements in rock joints.

Effect of layered joints on rockburst in deep tunnels

The existence of joints in the surrounding rock mass has a considerable efect on tunnel rockbursts.Herein,we studied the efect of layered joints with diferent inclination angles and spacings on rockburst in deep tunnels and investigated the failure area,deformation process of the surrounding rock mass,stress change inside the surrounding rock mass,velocity of the failed rock,and the kinetic energy of the failure.The failure type of the surrounding rock mass can thus be determined.The results showed that the intensity of rockburst increases as rock quality designation(RQD)decreases,while the deformation rate of the surrounding rock mass frst increases and then decreases.The deformation rate exhibits a turning point between RQD=50 and 70,below which the deformation rate of the surrounding rock mass gradually decreases,ultimately ceasing to be a rockburst.Rockburst always occurs perpendicular to the direction of the joint.Whenσ_(x)=σ_(y),as the joint inclination angle changes from 45°to 90°,the intensity of a rockburst frst decreases(from 45°to 60°),and then increases(from 60°to 90°).When combined with the evolution law of stress and strain energy,the rockburst process can be divided into four stages.

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.

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.

A roughness parameter considering joint material properties and peak shear strength model for rock joints

This study aims at proposing a reasonable roughness parameter that can reflect the peak shear strength(PSS)of rock joints.Firstly,the contribution of the asperities with different apparent dip angles to shear strength is studied.Then the shear strength of the entire joint asperities is derived.The results showed that the PSS of the entire joint asperities is proportional to a key parameter hs,which is related to the geometric character of the joint surface and the joint material properties.The parameter hsis taken as the new roughness parameter,and it is reasonable to associate the PSS with the geometric characteristics of the joint surface.Based on the new roughness parameter and shear test results of 20 sets of joint specimens,a new PSS model for rock joints is proposed.The new model is validated with the artificial joints in this paper and real rock joints in published studies.Results showed that it is suitable for different types of rock joints except for gneiss joints.The new model has the form of the Mohr-Coulomb model,which can directly reflect the relationship between the 3 D roughness parameters and the peak dilation angle.

Research on Complexity of Surface Undulating Shapes of Rock Joints

The surface undulating shapes of rock joints have been described qualitatively or experimental quantitatively for a long time. The non determined describing method can not fit quantitative evaluation of mechanical parameters of rock joints in engineering. In this paper, relative amplitude ( R A) is chosen as a quantitative describing index of surface measurement of 1 023 surface undulating curves which conducted by profile curve device(PCD). We discuss the nonuniformity,anisotropy and unhomogeneity of surface undulating shapes of joints. A new method that analyzes the complexity of surface undulating shapes of rock joints directional statistically in various rock joints is also put forward.

Laboratory-scale investigation of response characteristics of liquid-filled rock joints with different joint inclinations under dynamic loading

In underground rock engineering,water-bearing faults may be subjected to dynamic loading,resulting in the coupling of hydraulic and dynamic hazards.Understanding the interaction mechanism between the stress waves induced by dynamic loadings and liquid-filled rock joints is therefore crucial.In this study,an auxiliary device for simulating the liquid-filled layer was developed to analyze the dynamic response characteristics of liquid-filled rock joints in laboratory.Granite and polymethyl methacrylate(PMMA)specimens were chosen for testing,and high-amplitude shock waves induced by a split Hopkinson pressure bar(SHPB)were used to produce dynamic loadings.Impact loading tests were conducted on liquid-filled rock joints with different joint inclinations.The energy propagation coefficient and peak liquid pressure were proposed to investigate the energy propagation and attenuation of waves propagating across the joints,as well as the dynamic response characteristics of the liquid in the liquid-filled rock joints.For the inclination angle range considered herein,the experimental results showed that the energy propagation coefficient gently diminished with increasing joint inclination,and smaller coefficient values were obtained for granite specimens compared with PMMA specimens.The peak liquid pressure exhibited a gradually decreasing trend with increasing joint inclination,and the peak pressure for granite specimens was slightly higher than that for PMMA specimens.Overall,this paper may provide a considerably better method for studying liquid-filled rock joints at the laboratory scale,and serves as a guide for interpreting the underlying mechanisms for interactions between stress waves and liquid-filled rock joints.

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.

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.

An approach to measure infill matric suction of irregular infilled rock joints under constant normal stiffness shearing

Rock joints infilled with sediments can strongly influence the strength of rock mass. As infilled joints often exist under unsaturated condition, this study investigated the influence of matric suction of infill on the overall joint shear strength. A novel technique that allows direct measurement of matric suction of infill using high capacity tensiometers(HCTs) during direct shear of infilled joints under constant normal stiffness(CNS) is described. The CNS apparatus was modified to accommodate the HCT and the procedure is explained in detail. Joint specimens were simulated by gypsum plaster using threedimensional(3D) printed surface moulds, and filled with kaolin and sand mixture prepared at different water contents. Shear behaviours of both planar infilled joints and rough joints having joint roughness coefficients(JRCs) of 8-10 and 18-20 with the ratios of infill thickness to asperity height(t/a)equal to 0.5 were investigated. Matric suction shows predominantly unimodal behaviour during shearing of both planar and rough joints, which is closely associated with the variation of unloading rate and volumetric changes of the infill material. As expected, two-peak behaviour was observed for the rough joints and both peaks increased with the increase of infill matric suction. The results suggest that the contribution of matric suction of infill on the joint peak normalised shear stress is relatively independent of the joint roughness.