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.

Numerical modelling of resonance mechanisms in jointed rocks using transfer functions

Resonance effects in parallel jointed rocks subject to stress waves are investigated using transfer functions,derived from signals generated through numerical modelling.Resonance is important for a range of engineering situations as it identifies the frequency of waves which will be favourably transmitted.Two different numerical methods are used for this study,adopting the finite difference method and the combined discrete element-finite difference method.The numerical models are validated by replicating results from previous studies.The two methods are found to behave similarly and show the same resonance effects;one operating at low frequency and the other operating at relatively high frequency.These resonance effects are interpreted in terms of simple physical systems and analytical equations are derived to predict the resonant frequencies of complex rock masses.Low frequency resonance is shown to be generated by a system synonymous with masses between springs,described as spring resonance,with an equal number of resonant frequencies as the number of blocks.High frequency resonance is generated through superposition of multiple reflected waves developing standing waves within intact blocks,described as superposition resonance.While resonance through superposition has previously been identified,resonance based on masses between springs has not been previously identified in jointed rocks.The findings of this study have implications for future analysis of multiple jointed rock masses,showing that a wave travelling through such materials can induce other modes of propagation of waves,i.e.spring resonance.

A 3D microseismic data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions and its application

Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,it is always challenging to accurately measure rock mass properties.In this study,a three-dimensional(3D)microseismic(MS)data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions is proposed.It is a 3D finite element model that takes seepage,damage and stress field effects into account jointly.Multiple factors(i.e.joints,water and microseismicity)are used to optimize the rock mass mechanical parameters at different scales.The model is applied in Shirengou iron mine to study the damage evolution of rock mass and assess the crown pillar stability during the transition from open-pit to underground mining.It is found that the damage pattern is mostly controlled by the structure,water and rock mass parameters.The damage pattern is evidently different from the two-dimensional result and is more consistent with the field observations.This difference is caused by the MS-derived damage acting on the rock mass.MS data are responsible for gradually correcting the damage zone,changing the direction in which it expands,and promoting it to evolve close to reality.For the crown pillar,the proposed model yields a more trustworthy safety factor.In order to guarantee the stability of the pillar,it is suggested to take waterproof and reinforcement measures in areas with a high degree of damage.

Assessment of cyclic deformation and critical stress amplitude of jointed rocks via cyclic triaxial testing

Jointed rock specimens with a natural replicated joint surface oriented at a mean dip angle of 60were prepared,and a series of cyclic triaxial tests was performed at different confining pressures and cyclic deviatoric stress amplitudes.The samples were subjected to 10,000 loading-unloading cycles with a frequency of 8 Hz.At each level of confining pressure,the applied cyclic deviatoric stress amplitude was increased incrementally until excessive deformation of the jointed rock specimen was observed.Analysis of the test results indicated that there existed a critical cyclic deviatoric stress amplitude(i.e.critical dynamic deviatoric stress)beyond which the jointed rock specimens yielded.The measured critical dynamic deviatoric stress was less than the corresponding static deviatoric stress.At cyclic deviatoric stress amplitudes less than the critical dynamic deviatoric stress,minor cumulative residual axial strains were observed,resulting in hysteretic damping.However,for cyclic deviatoric stresses beyond the critical dynamic deviatoric stress,the plastic strains increased promptly,and the resilient moduli degraded rapidly during the initial loading cycles.Cyclic triaxial test results showed that at higher confining pressures,the ultimate residual axial strain attained by the jointed rock specimen decreased,the steadystate dissipated energy density and steady-state damping ratio per load cycle decreased,while steadystate resilient moduli increased.

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.

Fracture behaviors of columnar jointed rock mass using interface mechanics theorem

For a special geological structure of columnar jointed rock mass(CJRM),its mechanical properties are strongly affected by the columnar joints.To describe the fracture behaviors of CJRM using the basic theories of interface mechanics for composite materials,the interface stresses of the vertical and horizontal joints,which are the two primary joints in the CJRM under triaxial compression,are studied,and their mathematical expressions are derived based on the superposition principle.Based on the obtained interface stresses of the vertical and horizontal joints in the CJRM,the crack initiation of the joint interface in the CJRM is studied using the maximum circumferential stress theory in fracture mechanics.Moreover,based on this investigation,the fracture behaviors of CJRM are analyzed.According to the results of similar material physical model tests for the CJRM,the theoretical study is verified.Finally,the influence of the mechanical parameters of the CJRM on the joint interface stress is discussed comprehensively.

Particle flow study on strength and meso-mechanism of Brazilian splitting test for jointed rock mass

A discrete element method （DEM） called particle flow code （PFC2D） was used to construct a model for Brazilian disc splitting test in the present study. Based on the experimental results of intact Brazilian disc of rock-like material, a set of micro-parameters in PFC2D that reflected the macro-mechanical behavior of rock-like materials were obtained. And then PFC2D was used to simulate Brazilian splitting test for jointed rock mass specimens and specimen containing a central straight notch. The effect of joint angle and notch angle on the tensile strength and failure mode of jointed rock specimens was detailed analyzed. In order to reveal the meso-mechanical mechanism of crack coalescence, displacement trend lines were applied to analyze the displacement evolution during the crack initiation and propagation. The investigated conclusions can be described as follows. （1） The tensile strength of jointed rock mass disc specimen is dependent to the joint angle. As the joint angle increases, the tensile strength of jointed rock specimen takes on a nonlinear variance. （2） The tensile strength of jointed rock mass disc specimen containing a central straight notch distributes as a function of both joint angle and notch angle. （3） Three major failure modes, i.e., pure tensile failure, shear failure and mixed tension and shear failure mode are observed in jointed rock mass disc specimens under Brazilian test. （4） The notch angle roles on crack initiation and and joint angle play important propagation characteristics of jointed rock mass disc specimen containing a central straight notch under Brazilian test.

Responses of jointed rock masses subjected to impact loading

Impact-induced damage to jointed rock masses has important consequences in various mining and civil engineering applications. This paper reports a numerical investigation to address the responses of jointed rock masses subjected to impact loading. It also focuses on the static and dynamic properties of an intact rock derived from a series of laboratory tests on meta-sandstone samples from a quarry in Nova Scotia, Canada. A distinct element code(PFC2D) was used to generate a bonded particle model(BPM) to simulate both the static and dynamic properties of the intact rock. The calibrated BPM was then used to construct large-scale jointed rock mass samples by incorporating discrete joint networks of multiple joint intensities into the intact rock matrix represented by the BPM. Finally, the impact-induced damage inflicted by a rigid projectile particle on the jointed rock mass samples was determined through the use of the numerical model. The simulation results show that joints play an important role in the impactinduced rock mass damage where higher joint intensity results in more damage to the rock mass. This is mainly attributed to variations of stress wave propagation in jointed rock masses as compared to intact rock devoid of joints.

Determination of geological strength index of jointed rock mass based on image processing

The geological strength index(GSI) system,widely used for the design and practice of mining process,is a unique rock mass classification system related to the rock mass strength and deformation parameters based on the generalized Hoek-Brown and Mohr-Coulomb failure criteria.The GSI can be estimated using standard chart and field observations of rock mass blockiness and discontinuity surface conditions.The GSI value gives a numerical representation of the overall geotechnical quality of the rock mass.In this study,we propose a method to determine the GSI quantitatively using photographic images of in situ jointed rock mass with image processing technology,fractal theory and artificial neural network(ANN).We employ the GSI system to characterize the jointed rock mass around the working in a coal mine.The relative error between the proposed value and the given value in the GSI chart is less than 3.6%.

Seismic stability of jointed rock slopes under obliquely incident earthquake waves

Seismic stability of slopes has been traditionally analyzed with vertically propagated earthquake waves.However,for rock slopes,the earthquake waves might approach the outcrop still with a evidently oblique direction.To investigate the impact of obliquely incident earthquake excitations,the input method for SV and P waves with arbitrary incident angles is conducted,respectively,by adopting the equivalent nodal force method together with a viscous-spring boundary.Then,the input method is introduced within the framework of ABAQUS software and verified by a numerical example.Both SV and P waves input are considered herein for a 2 D jointed rock slope.For the jointed rock mass,the jointed material model in ABAQUS software is employed to simulate its behavior as a continuum.Results of the study show that the earthquake incident angles have significance on the seismic stability of jointed rock slopes.The larger the incident angle,the greater the risk of slope instability.Furthermore,the stability of the jointed rock slopes also is affected by wave types of earthquakes heavily.P waves induce weaker responses and SV waves are shown to be more critical.

Rate-dependent mechanical behavior of jointed rock with an impersistent joint under different infill conditions

Transition in the rate-dependent mechanical response of rock was investigated due to the presence of impersistent joint with different infill conditions.Four types of samples,i.e.intact,jointed with no grouting,jointed and grouted with cement,and jointed and grouted with epoxy,were fabricated using model material.A series of dynamic split Hopkinson pressure bar(SHPB) tests was conducted on prepared samples with strain rates varying between 53-130 salong with static uniaxial compression tests(10s).Progression of fracture/failure along samples was monitored using high-speed imaging and digital image correlation(DIC).Strength was observed to be significantly lower for jointed samples as compared to intact samples.However,the increasing trend of strength with strain rates remained similar for all types of samples.Epoxy was observed to be a better grout due to greater improvement in the strength of epoxy grouted jointed samples than cement grouted samples under both static and dynamic conditions.Significant changes were observed in fracture behavior(initiation,pattern and mechanism)with strain rate for intact and jointed unfilled/grouted samples.Fracturing was dominated by shear and tensile cracks at high strain rates compared to tensile cracks at low strain rates in all samples.Unlike static loading conditions,the location of cracks initiation shifts away from joint tips with increasing strain rate and depending upon existing infill conditions(unfilled/grouted).

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.

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.

A new estimation method and an anisotropy index for the deformation modulus of jointed rock masses

The deformation modulus of a rock mass is an important parameter to describe its mechanical behavior.In this study,an analytical method is developed to determine the deformation modulus of jointed rock masses,which considers the mechanical properties of intact rocks and joints based on the superposition principle.Due to incorporating the variations in the orientations and sizes of joint sets,the proposed method is applicable to the rock mass with persistent and parallel joints as well as that with nonpersistent and nonparallel joints.In addition,an anisotropy index AIdmfor the deformation modulus is defined to quantitatively describe the anisotropy of rock masses.The range of AIdmis from 0 to 1,and the more anisotropic the rock mass is,the larger the value of AIdmwill be.To evaluate the proposed method,20 groups of numerical experiments are conducted with the universal distinct element code(UDEC).For each experimental group,the deformation modulus in 24 directions are obtained by UDEC(numerical value)and the proposed method(predicted value),and then the mean error rates are calculated.Note that the mean error rate is the mean value of the error rates of the deformation modulus in 24 directions,where for each direction,the error rate is equal to the ratio of numerical value minus predicted value to the numerical value.The results show that(i)for different experimental groups,the mean error rates vary between 5.06%and 22.03%;(ii)the error rates for the discrete fracture networks(DFNs)with two sets of joints are at the same level as those with one set of joints;and(iii)therefore,the proposed method for estimating the deformation modulus of jointed rock masses is valid.

Shear behaviour of anchored jointed rock mass under different engineering conditions considering damage and joint surface morphology

According to different geological conditions and engineering disturbance,three dimensional numerical models of anchored jointed rock mass with two kinds of boundary conditions of constant normal load(CNL)and constant normal stiffness(CNS)were constructed considering the ductility damage of rockbolt,the stiffness degradation of grouting layer and the joint surface roughness.The effects of anchorage angle,joint surface morphology,and boundary conditions on the shear performance of anchorage system were analyzed,and the failure characteristics under different working conditions were revealed.Finally,the analytical solution of shear strength of anchored system was established.Results show that the larger the anchorage angle is,the more serious the necking phenomenon of rockbolt will be.The damage degree of the bonding layer and the horizontal displacement of the bedding rock mass decrease with the increase of the joint surface roughness.The CNL condition is to instantaneously apply high normal stress,and the CNS condition is to gradually form a high normal stress environment through the superposition of increment on the basis of the initial value,which can resist greater transverse load.This is equivalent to enhancing the ductility of the rockbolt.The shear strength of the system increases with the increase of normal stress and normal stiffness.Ignoring the normal stiffness will underestimate the shear strength.

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.

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.

Characterization on jointed rock masses based on PFC2D

Geometrical parameters of discontinuities,such as spacing,length,dip and fault throw between joints have a great influence on the mechanical behavior of jointed rock masses.Accurate characterization for discontinuities is important for investigate the stability of rock masses.In this paper,the PFC2D is combined with joint network generation method to examine the mechanical behaviors of jointed mass.Taking Miaogou Open-pit Mine as an example,the information and statistical distributions of discontinuities of the slope rock masses are measured by ShapeMetriX3D measuring tool.Then,the automatic generation algorithm of random joints network based on the Monte-Carlo method is proposed using the programming language(FISH)embedded within PFC2D.This algorithm could represent the discontinuities compared with the geological surveys.In simulating the compression test of a jointed rock sample,the mechanical behavior and crack propagation were investigated.The results reveal that the failure mode and crack propagation of the jointed rock are dominated by the distribution of joints in addition to the intact rock properties.The simulation result shows the feasibility of the joints generating method in application to jointed rock mass.

STRENGTH PROPERTIES OF JOINTED ROCK MASSES BASED ON THE HOMOGENIZATION METHOD

This paper aims to determine the strength properties of jointed rock masses by means of the homogenization method. To reflect the microstructure of jointed rock masses, a representative element volume （REV） is selected. Assuming being rigid and perfectly plastic and utilizing the Mohr-Coulomb yield criterion for rock and joints, an upper bound limit analysis method is proposed based on the homogenization theory. By using the finite element method and Goodman joint element, a nonlinear mathematical programming with equality constraints is formulated, which can be solved by a direct iterative algorithm. Numerical results show the strength of jointed rock masses behaves anisotropy with different joint directions. This method presents an effective tool for the strength analysis of jointed rock masses.

Simulation of blast induced crater in jointed rock mass by discontinuous deformation analysis method

Rock blasting is a dynamic process accom panied with the propagations of shock waves and the dispersion of the explosion gas.This paper adopts the discontinuous deformation analysis(DDA)method to simulate the rock blasting process.A dynamic parameter adjustment and the non-reflecting boundary condition are implemented in the DDA method.The sub-block DDA method to simulate fracture problems is used.The blasting process in jointed rock mass is simulated by application of the explosion gas pressure on the expanding borehole walls and induced connected fracture surfaces around the boreholes.The blast craters with different overburdens are derived.The whole process including the explosion gas dispersion,borehole expansion,rock mass failure and cast,and the formation of the final blasting piles in rock blasting are well reproduced numerically.Parametric study for different overburdens is carried out,and the results are analyzed and discussed.