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.

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.

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 modeling of a tunnel in jointed rocks subjected to seismic loading

The studies on the performance of tunnels under static loads are reported extensively in the literature but their performances under dynamic loads are limited.The present study highlights some of the important aspects of jointed rock tunnels during seismic loading.The literature review provides a shake table experimental study of a jointed rock tunnel.A Universal Distinct Element Code(UDEC)model is developed from this shake table experiment.The model tunnel is subjected to a scaled input motion of the 1985 Mexico earthquake.The numerical results are validated systematically with the findings of the shake table experiment.Further,the developed numerical model is used to perform parametric studies to understand the effect of in-situ stress,joint angles,joint stiffness,and joint friction angle on the deformation and stability of the tunnel under the same earthquake input motion.It is observed that some joint angle combinations form a wedge that yields excessive deformation and subsequently a complete failure.An exponential decrease in deformation occurred in the tunnel as the joint stiffness increases.It is found that the shallow tunnels are more susceptible to damage under the action of earthquake loads.

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

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

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

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.

Analysis of the Elastic-Plastic Dynamic Finite Element on a Jointed Rock Mass

Aim To study the elastic plastic dynamical constitutive relations about a jointed rock mass under explosion load and its computer simulation. Methods\ Stress history is taken into account and stresses will follow changes in time during a period of explosion load. According to the principle of static force balance, the corresponding nodal concentrated force is calculated and the nodal displacement is counted. The elastic plastic dynamic finite element equations are thus obtained. Results\ A finite element method is given for a jointed rock mass under explosion load. Conclusion\ The problem of large plastic deformation for jointed rock mass on blasting was efficiently resolved through dynamic finite element analysis and the range of damages by blasting simulated, and this pushes forward the problem to engineering practice.

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

Analysis profile of the fully grouted rock bolt in jointed rock using analytical and numerical methods

The purpose of this study was to investigate the effect of bolt profile on load transfer mechanism of fully grouted bolts in jointed rocks using analytical and numerical methods. Based on the analytical method with development of methods, a new model is presented. To validate the analytical model, five different profiles modeled by ANSYS software. The profile of rock bolts T3 and T4with load transfer capacity,respectively 180 and 195 kN in the jointed rocks was selected as the optimum profiles. Finally, the selected profiles were examined in Tabas Coal Mine. FLAC analysis indicates that patterns 6+7 with2 NO flexi bolt 4 m better than other patterns within the faulted zone.

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.

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

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

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.

Optimizing profile line interval for enhanced accuracy in rock joint morphology and shear strength assessments

2D profile lines play a critical role in cost-effectively evaluating rock joint properties and shear strength.However, the interval(ΔI_(L)) between these lines significantly impacts roughness and shear strength assessments. A detailed study of 45 joint samples using four statistical measures across 500 different ΔI_(L)values identified a clear line interval effect with two stages: stable and fluctuation-discrete.Further statistical analysis showed a linear relationship between the error bounds of four parameters,shear strength evaluation, and their corresponding maximum ΔI_(L)values, where the gradient k of this linear relationship was influenced by the basic friction angle and normal stress. Accounting for these factors,lower-limit linear models were employed to determine the optimal ΔI_(L)values that met error tolerances(1%–10%) for all metrics and shear strength. The study also explored the consistent size effect on joints regardless of ΔI_(L)changes, revealing three types of size effects based on morphological heterogeneity.Notably, larger joints required generally higher ΔI_(L)to maintain the predefined error limits, suggesting an increased interval for large joint analyses. Consequently, this research provides a basis for determining the optimal ΔI_(L), improving accuracy in 2D profile line assessments of joint characteristics.

Energy evolution mechanism and failure criteria of jointed surrounding rock under uniaxial compression

The object of this article is to investigate the energy evolution mechanism and failure criteria of cross-jointed samples containing an opening during deformation and failure based on the uniaxial compression test and rock energy principle.The results show that the energy evolution characteristics of the samples correspond to a typical progressive damage mode.The peak total energy,peak elastic energy,and total input energy of the samples all first decrease and then increase with an increase of half of the included angle,reaching their minimum values when this angle is 45°,while the dissipated energy generally increases with this angle.The existence of the opening and cross joints can obviously weaken the energy storage capacity of the rock,and the change in the included angle of the cross joint has a great influence on the elastic energy ratio of the sample before the peak stress,which leads to some differences in the distribution laws of the input energy.The continuous change and the subsequent sharp change in the rate of change in the energy consumption ratio can be used as the criteria of the crack initiation and propagation and the unstable failure of the sample,respectively.

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

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

Macro and meso analysis of jointed rock mass triaxial compression test by using equivalent rock mass(ERM) technique

Methods that can efficiently model the effects of rock joints on rock mass behavior can be beneficial in rock engineering. The suitability of equivalent rock mass(ERM) technique based upon particle methods is investigated. The ERM methodology is first validated by comparing calculated and experimental data of lab triaxial compression test on a set of cylindrical rock mass samples, each containing a single joint oriented in various dip angles. The simulated results are then used to study the stress-strain nonlinearity and failure mechanism as a function of the joint dip angle and confining stress. The anisotropy and size effects are also investigated by using multi-scale cubic ERM models subjected to triaxial compression test. The deformation and failure behavior are found to be influenced by joint degradation, the micro-crack formation in the intact rock, the interaction between two joints, and the interactions of micro-cracks and joints.

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.

Dynamic-static coupling analysis on rockburst mechanism in jointed rock mass

A numerical code called RFPA-Dynamics was used to study the rockburst mechanism under dynamic load based on coupled static-dynamic analysis.The results show that dynamic disturbance has a very distinct triggering effect on rockburst.Under the dynamic load,rockburst is motivated by tensile stress formed by the overlapping of dynamic waves in the form of instantaneous open and cutting through of cracks in weak planes and pre-damaged areas.Meanwhile,the orientation of joint sets has an obvious leading effect on rockburst locations.Finally,a higher initial static stress state before dynamic loading can cause more pre-damaged area,thus leading to a larger rockburst scope.

Numerical investigation on the sensitivity of jointed rock mass strength to various factors

The mechanical properties of jointed rock masses, such as strength, deformation and the failure mechanism, can be understood only by studying the sensitivity of jointed rock mass strength (both the peak and residual strengths) to the factors that affect it. An orthogonal design of uniaxial compression tests was simulated on eighteen groups of jointed rock specimens having different geometric and mechanical properties using RFPA2D (Rock Failure Process Analysis) code. The results show that the peak strength is controlled by the geometric parameters of the joints, but that the residual strength is controlled by the mechanical prop- erties of the joint interfaces. The failure mode of jointed rock specimens is mainly shear failure. Joint quantity, or density, is the most important index that affects jointed rock mass strength and engineering quality.