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.

Shear behaviours and roughness degeneration based on a quantified rock joint surface description

The asperity wear of rock joints significantly affects their shear behaviour.This study discusses the wear damage of the asperities on the joint surface,highlighting the roughness degradation characteristics during the shear process.The direct shear experiment of artificial specimens containing rock joints was conducted under different normal stresses based on three-dimensional scanning technology.These experimental results showed the contribution of joint wear to roughness degeneration,such as the height,zone,and volume of asperity degeneration.The wear coefficient of the rock joint was obtained based on the volume wear of asperities in the laboratory experiment.The functional relationship between the friction coefficient and wear coefficient is subsequently determined.To quantitatively analyse the wear damage of a joint surface,a calculation method for determining the wear depth of the rock joint after shearing was proposed based on wear theory.The relationship between the ultimate dilation and wear depth was analysed.A coefficient m,which can describe the damage degree of the joint surface,and a prediction method of joint surface roughness after shearing are established.Good agreement between analytical predictions and measured values demonstrates the capability of the developed model.Lastly,the sensitivity factors on the wear depth are explored.

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.

A comparative study for determining rock joint normal stiffness with destructive uniaxial compression and nondestructive ultrasonic wave testing

Rock joints are one of the vital discontinuities in a natural rock mass.How to accurately and conveniently determine joint normal stiffness is therefore significant in rock mechanics.Here,first,seven existing methods for determining joint normal stiffness were introduced and reviewed,among which MethodⅠ(the indirect measurement method),MethodⅡ(the direct determination method),MethodⅢ(the across-joint strain gauge measurement method)and MethodⅣ(the deformation measuring ring method)are via destructive uniaxial compression testing,while MethodⅤ(the best fitting method),MethodⅥ(the rapid evaluation method)and MethodⅦ(the effective modulus method)are through wave propagation principles and nondestructive ultrasonic testing.Subsequently,laboratory tests of intact and jointed sandstone specimens were conducted following the testing requirements and pro-cedures of those seven methods.A comparison among those methods was then performed.The results show that Method I,i.e.the benchmark method,is reliable and stable.MethodⅡhas a conceptual drawback,and its accuracy is acceptable at only very low stress levels.Relative errors in the results from MethodⅢare very large.With MethodⅣ,the testing results are sufficiently accurate despite the strict testing environment and complicated testing procedures.The results from MethodⅤare greatly unstable and significantly dependent on the natural frequency of the transducers.The joint normal stiffness determined with MethodⅥis stable and accurate,although data processing is complex.MethodⅦcould be adopted to determine the joint normal stiffness corresponding to the rock elastic deformation phase only.Consequently,it is suggested that MethodsⅠ,ⅣandⅥshould be adopted for the mea-surement of joint normal stiffness.The findings could be helpful in selecting an appropriate method to determine joint normal stiffness and,hence,to better solve discontinuous rock mass problems.

A novel method for geometric quality assurance of rock joint replicas in direct shear testing-Part 1:Derivation of quality assurance parameters and geometric reproducibility

Since each rock joint is unique by nature,the utilization of replicas in direct shear testing is required to carry out experimental parameter studies.However,information about the ability of the replicas to simulate the shear mechanical behavior of the rock joint and their dispersion in direct shear testing is lacking.With the aim to facilitate generation of high-quality direct shear test data from replicas,a novel component in the testing procedure is introduced by presenting two parameters for geometric quality assurance.The parameters are derived from surface comparisons of three-dimensional(3D)scanning data of the rock joint and its replicas.The first parameter,smf,captures morphological deviations between the replica and the rock joint surfaces.smf is derived as the standard deviation of the deviations between the coordinate points of the replica and the rock joint.Four sources of errors introduced in the replica manufacturing process employed in this study could be identified.These errors could be minimized,yielding replicas with smf
0.06 mm.The second parameter is a vector,VHp100,which describes deviations with respect to the shear direction.It is the projection of the 100 mm long normal vector of the best-fit plane of the replica joint surface to the corresponding plane of the rock joint.

VHp100

was found to be less than or equal to 0.36 mm in this study.Application of these two geometric quality assurance parameters demonstrates that it is possible to manufacture replicas with high geometric similarity to the rock joint.In a subsequent paper(part 2),smf and VHp100 are incorporated in a novel quality assurance method,in which the parameters shall be evaluated prior to direct shear testing.Replicas having parameter values below established thresholds shall have a known and narrow dispersion and imitate the shear mechanical behavior of the rock joint.

A novel method for geometric quality assurance of rock joint replicas in direct shear testing-Part 2:Validation and mechanical replicability

Each rock joint is unique by nature which means that utilization of replicas in direct shear tests is required in experimental parameter studies.However,a method to acquire knowledge about the ability of the replicas to imitate the shear mechanical behavior of the rock joint and their dispersion in direct shear testing is lacking.In this study,a novel method is presented for geometric quality assurance of replicas.The aim is to facilitate generation of high-quality direct shear testing data as a prerequisite for reliable subsequent analyses of the results.In Part 1 of this study,two quality assurance parameters,smf and V_(Hp100),are derived and their usefulness for evaluation of geometric deviations,i.e.geometric reproducibility,is shown.In Part 2,the parameters are validated by showing a correlation between the parameters and the shear mechanical behavior,which qualifies the parameters for usage in the quality assurance method.Unique results from direct shear tests presenting comparisons between replicas and the rock joint show that replicas fulfilling proposed threshold values of σ_(mf)<0.06 mm and
|V_(Hp100)|
<0.2 mm have a narrow dispersion and imitate the shear mechanical behavior of the rock joint in all aspects apart from having a slightly lower peak shear strength.The wear in these replicas,which have similar morphology as the rock joint,is in the same areas as in the rock joint.The wear is slightly larger in the rock joint and therefore the discrepancy in peak shear strength derives from differences in material properties,possibly from differences in toughness.It is shown by application of the suggested method that the quality assured replicas manufactured following the process employed in this study phenomenologically capture the shear strength characteristics,which makes them useful in parameter studies.

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.

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.

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.

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.

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 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.

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.

Modeling of Mechanical Behavior of Fractal Rock Joint

The present study shows that naturally developed fracture surfaces in rocks display the properties of self-affine fractals. Surface roughness can be quantitatively characterized by fractal dimension D and the intercept A on the log-log plot of variance: the former describes the irregularity and the later is statistically analogues to the slopes of asperities. In order to confirm the effects of these fractalparameters on the properties and mechanical behavior of rock joints, which have been observed in experiments under both normal andshear loadings, a theoretic model of rock joint is proposed on the basis of contact mechanics. The shape of asperity at contact is assumed to have a sinusoidal form in its representative scale r, with fractal dimension D and the intercept A. The model considers different local contact mechanisms, such as elastic deformation, frictional sliding and tensile fracture of the asperity. The empirical evolution law of surface damage developed in experiment is implemented into the model to up-date geometry of asperity in loading history. The effects of surface roughness characterized by D, A and re on normal and shear deformation of rock joint have been elaborated.

Stochastic seismic wave interaction with a rock joint having Coulomb-slip behavior

Seismic wave interaction with a slippery rock joint with an arbitrary impinging angle is analytically studied based on the conservation of momentum on the wave fronts. Based on the displacement discontinuity method, the wave propagation equations are derived for incident P- and S-waves. By comparison, the calculated transmission and reflection coefficients for normal incident waves are the same as the existing results, which proves the wave propagation equation obtained in the paper is correct. The wave propagation derived in the context can be applied to incident waves with different waveforms. Stochastic seismic waves are then used to analyze the seismic wave interaction with the slippery rock joint, where the stochastic seismic waves are generated from frequency spectra. The parametric studies are carried out to investigate the effect of type, intensity and impinging angle of the incident seismic waves on the wave propagation across the slippery rock joint.

Rock joint coefficients and their computerized classification

A computerized method for determining rock joint coefficients is presented.Two relative similarity indicators are introduced to classify surface morphology of rock joints.The classification enables to compare investigated and database rock joints.Such a comparison aims at finding the couple of surfaces that are distinguished by the highest dynamical conformity.The first absolute indicator results from the Fourier matrix and evaluates wavy shapes of surfaces.The second absolute indicator quantifies the heights of surface reliefs and is defined as the root mean square height of the surface outline.Numerical reliability of these indicators is tested within the surface analysis of a series of limestone specimens.Besides the computerized assessment,25 people have performed visual assessment of these limestone specimens.The results of visual assessments have been statistically processed and compared to the results received from the computerized procedure.The newly introduced absolute indicators have proved to be prospective numerical tools for evaluating joint rock coefficients.