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.

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.

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.

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.

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.

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.

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.

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.

Mechanical properties and reinforcement effect of jointed rock mass with pre-stressed bolt

Pre-stressed bolt anchorage is the key technology for jointed rock masses in rock tunnelling,slope treatment and mining engineering.To investigate the mechanical properties and reinforcement effect of jointed rock masses with pre-stressed bolts,in this study,uniaxial compression tests were conducted on specimens with different anchoring types and flaw inclination angles.ABAQUS software was used to verify and supplement the laboratory tests.The laws of the uniaxial compressive strength(UCS)obtained from the numerical simulations and laboratory tests were consistent.The results showed that under the same flaw angle,both the UCS and elastic modulus of the bolted specimens were improved compared with those of the specimens without bolts and the improvements increased with an increase in the bolt pre-stress.Under the same anchoring type,the UCS and elastic modulus of the jointed specimens increased with an increase in the flaw angle.The pre-stressed bolt could not only restrain the slip of the specimens along the flaw surface but also change the propagation mode of the secondary cracks and limit the initiation of cracks.In addition,the plot contours of the maximum principal strain and the Tresca stress of the numerical models were influenced by the anchoring type,flaw angle,anchoring angle and bolt position.

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

Simulation of failure process of jointed rock

A modified discontinuous deformation analysis (DDA) algorithm was proposed to simulate the failure behavior of jointed rock. In the proposed algorithm, by using the Monte-Carlo technique, random joint network was generated in the domain of interest. Based on the joint network, the triangular DDA block system was automatically generated by adopting the advanced front method. In the process of generating blocks, numerous artificial joints came into being, and once the stress states at some artificial joints satisfy the failure criterion given beforehand, artificial joints will turn into real joints. In this way, the whole fragmentation process of rock mass can be replicated. The algorithm logic was described in detail, and several numerical examples were carried out to obtain some insight into the failure behavior of rock mass containing random joints. From the numerical results, it can be found that the crack initiates from the crack tip, the growth direction of the crack depends upon the loading and constraint conditions, and the proposed method can reproduce some complicated phenomena in the whole process of rock failure.

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.

Elasto-plastic Analysis of Circular Tunnels in Jointed Rock Masses Satisfy the Hoek-Brown Failure Criterion

The relationship between the Hoek-Brown parameters and the mechanical response of circular tunnels is il-lustrated. Closed-form and approximate solutions are given for the extent of the plastic zone and the stress and dis-placement fields under axisymmetrical and asymmetric stress conditions. For the same rock masses and under axisym-metrical stress conditions,the radius of the plastic zone in terms of Hoek-Brown criterion is generally an approximation of the radius in terms of the Mohr-Coulomb criterion. The radius in terms of the Hoek-Brown criterion is larger under low stress conditions. For poor quality rock masses (GSI<25),measures (such as grouting,setting rock bolts,etc.) that improve the GSI of rock masses are effective in improving the stability of tunnels. It is not advisable to improve the sta-bility of the tunnels by providing a small support resistance p through shotcrete,except for very poor quality jointed rock masses. Without reference to the quality of the rock mass,the disturbance factor D should not less than 0.5. Meas-ures which disturb rock masses during tunnel construction should be taken carefully when the tunnel depth increases.

Experimental study on acoustic emission characteristics of jointed rock mass by double disc cutter

The characteristics of joints are crucial factors which influence the penetration efficiency of tunnel boring machine(TBM).Based on the theoretical study,numerical simulation and experimental research,many researchers have studied the interaction between TBM disc cutters and jointed rock mass.However,in most of these works,the effect of joint on rock fragmentation by double disc cutter has been scarcely investigated.Thus,the effects of joint orientation and joint space on rock fragmentation by double disc cutter are highlighted in this study.During the test,jointed concrete specimens are adopted to simulate jointed rock mass.Improved RYL-600rock shear rheological instrument was employed during the indentation process under disc cutters,and acoustic emission location system was used to analyze the rock damage and physical deterioration.The results show that there are four failure modes and three modes of crack initiation and propagation in jointed rock mass.It is concluded that the existing joint planes have obviously restrained the crack initiation and propagation during the rock fragmentation process.The results also indicate that samples are damaged most seriously when joint orientation equals60°,which is proved to be the optimum joint orientation in TBM penetration.

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.