Behavior Of Rock Mass Research Articles

Numerical Investigation of the Slope Stability in the Waste Dumps of Romanian Lignite Open-Pit Mines Using the Shear Strength Reduction Method

Open-pit mining generates significant amounts of waste material, leading to the formation of large waste dumps that pose environmental risks such as land degradation and potential slope failures. The paper presents a stability analysis of waste dump slopes in open-pit mining, focusing on the Motru coalfield in Romania. To assess the stability of these dumps, the study employs the Shear Strength Reduction Method (SSRM) implemented in the COMSOL Multiphysics version 6 software, considering both associative and non-associative plasticity models. (1) Various slope angles were analyzed, and the Factor of Safety (FoS) was calculated, showing that the FoS decreases as the slope angle increases. (2) The study also demonstrates that the use of non-associative plasticity leads to lower FoS values compared to associative plasticity. (3) The results are visualized through 2D and 3D models, highlighting failure surfaces and displacement patterns, which offer insight into the rock mass behavior prior to failure. (4) The research also emphasizes the effectiveness of numerical modeling in geotechnical assessments of stability. (5) The results suggest that a non-associative flow rule should be adopted for slope stability analysis. (7) Quantitative results are obtained, with small variations compared to those obtained by LEM. (6) Dilatation angle, soil moduli, or domain changes cause differences of just a few percent and are not critical for the use of the SSRM in engineering.

Deformation behaviour of strain-softening rock mass in tunnels considering deterioration model of elastic modulus

By analyzing the composite effects of the plastic shear strain and confining stress according to the laboratory test results, a model adopting a deteriorating elastic modulus for the intact rock is proposed. Soft rock such as coal and the strong rock such as granite with different quality are investigated. Geological Strength Index (GSI) is adopted to represent the integrity of rock, to be specific, from intact rock to the rock mass. Specifically, for the strain-softening rock mass, the elastic modulus is gained by introducing GSI into the deterioration model. The strength parameters are obtained by fitting GSI into Hoek-Brown failure criterion. Then, the elastic modulus and strength parameters of the rock mass are employed in the numerical procedure for solving out the rock deformation, the radii of the plastic softening and residual zones of the tunnels. The rationality of the numerical procedure is verified with the existing semi-analytical and analytical procedures. In the discussion, the influences of the critical plastic parameters and GSI, on the elastic modulus, rock deformation, radii of the plastic softening and residual zones are investigated. Five models for the variation of the modulus including the proposed deterioration model are defined. The rock mass deformation behaviours of the five models are compared. The results indicate that, in the plastic softening zone that is far away from the tunnel periphery, the elastic modulus is more sensitive to the plastic shear strain, whereas near the tunnel periphery, the elastic modulus is primarily affected by the confining stress. Regarding a linear decrease of the elastic modulus versus the plastic shear strain overestimates the elastic modulus to some extent, especially for the soft rock mass.

Experimental study on interaction characteristics of explosive cracks under confining pressure and its comparison with free boundary results

The influence of confining pressure on the dynamic fracture mechanical behavior of deep rock masses is crucial in the exploitation of deep resources and the development of underground spaces. In this study, dynamic caustics experiments were conducted using transparent epoxy resin materials to investigate the crack propagation behavior and interaction mechanisms of double cracks under confinement. The experimental results indicate that biaxial confining pressure suppresses crack propagation, interaction between double cracks, and penetration. Without confining pressure, cracks do not penetrate until the relative vertical distance ΔH between crack tips is 20 mm, whereas with confining pressure, cracks cease to penetrate when the relative vertical distance between crack tips is 10 mm; however, interaction between crack tips still exists at this point. During the crack interaction stage, both fracture toughness and propagation velocity of crack tips exhibit significant fluctuations under no confining pressure conditions, whereas under confining pressure, both decrease rapidly. Additionally, under confining pressure conditions, the fluctuation range of crack propagation direction during crack interaction is relatively small. Confining pressure rapidly reduces the degree of stress concentration at crack tips, increases the rate of energy release at crack tips, and enhances the ability of dynamic cracks to resist the influence of other crack tips or crack surfaces. Under confining pressure, continuous Rayleigh waves are typically observed at crack tips due to an imbalance between stored strain energy and the energy release rate. These waves dissipate part of the strain energy, promoting crack deflection and branching.

Numerical simulation of failure and micro-cracking behavior of non-persistent rock joint under direct shear

Persistency, as a key geometric parameter of joints, significantly affects shear strength parameters of jointed rock mass. A good understanding of how persistency affects shear behavior of joint is therefore crucial for better evaluation of stability of rock slope. To investigate the failure and micro-cracking behavior of non-persistent rock joint under direct shear, a novel Voronoi generation algorithm is first used to establish an improved grain-based model (GBM) of granite which considers the shape of feldspar. The calibrated model is then used to simulate the direct shear test of numerical models possessing different joint persistency under various normal stresses. The results reveal that the developed micro-cracks generally increase rapidly when the shear strain reaches to a value approximately 50 % of the peak shear strain and the grain boundary tensile micro-crack is dominant among these initiated micro-cracks. Micro-cracks generally initiate at the ends of rock bridge, and then gradually propagate to the central of rock bridge, forming en-echelon fractures. The failure mode of numerical model is closely related to the generated en-echelon fractures. An increase in both joint persistency and normal stress can lead to a shear failure. The finding in this study provides an important basis for understanding the mechanical behavior and failure mechanism of jointed rock mass.

Experimental study on mechanical and acoustic emission characteristics of red sandstone containing combined flaws of hole-fissure under biaxial compression

To investigate the mechanical behaviour and failure mechanism of flawed rock masses under biaxial stress conditions, a series of biaxial compression experiments are conducted on red sandstone specimens containing combined flaws of the circular hole and perforated symmetric fissure combined with acoustic emission (AE) technology in this work. The results show a close association between the stress-strain curve morphology and confining stress, both the biaxial compression strength and elastic modulus exhibit an upward trend with increasing fissure angle and confining stress. The maximum AE energy and cumulative AE energy both increase with higher confining stress. The crack types in the specimens are identified by analyzing the AF/RA value distribution, revealing a gradual decrease in the percentage of tensile cracks with increasing fissure angle. An AE localization algorithm based on the least absolute value method is applied to pinpoint AE events during biaxial compression, and AE events distribution is analyzed through the Kernel density estimation (KDE). The crack extension behaviour is described based on the movement of the maximum kernel density points, which initially exhibits a progression from both ends of the specimen towards the central region and subsequently from the central fissure tip towards the two ends.

The field-enriched finite element method with gravity effects for simulating the cracking behaviors of large-scale engineering rock masses

The field-enriched finite element method uses a scalar field defined as a field variable to describe cracks and characterize their impact on the displacement field and stress field of the solution model. It is capable of avoiding remeshing and employing level set functions to describe cracks when simulating the propagation of cracks. In this work, a field-enriched finite element model with gravity effects is proposed to simulate the large-scale failure process of engineering rock masses, and several numerical cases of geotechnical engineering are successfully analyzed. First, by introducing the unified tensile fracture criterion into the numerical model, the large-scale failure process of the intact slope is simulated. Second, the sliding process of rock slopes containing en echelon joints is numerically investigated. Third, the cracking process of the concrete dam is analyzed. Finally, the effects of joint and bedding plane inclination angles on the stability of tunnel chamber in transversely isotropic rock mass are studied. The numerical results indicate that the numerical method proposed in this work can accurately solve the large-scale failure process of rock masses.

Contribution to the numerical modelling analysis of displacements caused by tunnel construction in discontinuous rock masses

This paper presents a comprehensive study of the stability conditions of a rock mass surrounding a tunnel, using numerical modelling analysis of displacements induced by the tunnel construction. The case study focuses on the Boukhadra iron ore mine in Algeria. This research employs empirical equations from literature based on geomechanical classifications to provide the most appropriate geotechnical input data for simulation. In addition to those equations, a novel correlation equation relating the rock mass ratio (RMR) and the rock quality index (Q) was developed. This equation gives a high regression coefficient, revealing a strong correlation (R2) of 0.878. It gives a better estimation of the rock mass characteristics: Young’s moduli (E), Poisson’s ratio (υ), cohesion (C), and friction angle (φ), compared to the existing literature-based correlation equations. The estimated parameters were used to pass from a discontinuous to a continuous equivalent medium, making the numerical modelling easier with the finite element method. Compared with intact rock and direct equations, the numerical simulation results based on the new equation fit well with the in-situ observations and provide a more reliable comprehension of rock mass behaviour. The new equation proposed in this paper provides a substantial evaluation of the rock mass parameters needed in the design process of underground structures once used with respect to their acceptable ranges.

Individual-portal support with a pre-tensioned canopy for a coal mine roadway junction

Abstract This paper presents an alternative model of individual-portal support structure with a pre-tensioned canopy for coal mine roadway junctions in the geological and mining conditions of underground mines in Vietnam. The aim of the research is to determine the bearing capacity of the individual-portal support structure with various configuration. For this purpose, an analysis of the rock mass behaviour around the three-way roadway junction was carried out using the difference element method-based programme (FLAC3D). Additionally, other analysis of the bearing capacity of the selected individual-portal support structure was also conducted using the finite element method-based programme (COSMOS/M). The geological and mining conditions of Vietnamese mines were taken into account during performing these analyses. Based on the results, the optimal individual-portal support structure was proposed for the roadway junctions driven in the geological and mining conditions of underground mines in Vietnam.

Analysis of failure properties of red sandstone with structural plane subjected to true triaxial stress paths

In order to study the mechanical behavior and failure characteristics of rock mass with the structural plane in complex deep stress environment, the QKX-YB200 true triaxial rockburst test system, acoustic emission, and digital video recorder were adopted to conduct true triaxial loading and unloading indoor tests on the red sandstone cube specimens. The results show that: (1) Under true triaxial loading conditions, all specimens propagate in the form of anti-wing crack, and the dip angle of the structural plane is linearly related to the average value of the anti-wing crack propagation angle; under true triaxial unloading conditions, the specimens with structural plane dip angles of 30°, 45°, and 60° propagate in the form of wing crack, while those with a 90° angle propagate anti-wing cracks, and the closure of structural plane is related to the the crack type at the tips of structural plane. (2) The unstable failure stages in the stress–strain curve under true triaxial unloading are longer and exhibit more pronounced fluctuations; the peak strength of the specimens under the conditions of true triaxial loading and unloading decrease first and then increase with the increase of the dip angle, reaching the lowest at 45°. (3) During the true triaxial unloading tests, the rockburst degree of the intact rock specimens is more violent than that of the specimens with a structural plane. Additionally, the dip angle of the structural plane influences the severity of rockbursts under the same stress path. (4) The crack propagation mechanism and the closure degree of structural planes under different true triaxial stress paths and dip angles are discussed. It is observed that the crack initiation mode around the structural plane tips and the crack mode near free surface are closely related to the stress environment in different areas of the specimens. (5) Based on the indoor test results and previous experimental data, the differences in peak strength characteristics between open structural planes and closed structural planes are preliminarily analyzed and compared.

Review on Constitutive Model for Simulation of Weak Rock Mass

Understanding the behavior of weak rock masses is important for predicting the stability of structures under different loading conditions. Traditional models such as the generalized Hoek–Brown and Coulomb weak plane are widely used; however, they often fail to capture the nonlinear and irreversible behavior of weak rock masses. This study offers a comprehensive overview of a critical analysis of constitutive models’ strengths and limitations for simulating weak rock masses. By comparing traditional and advanced novel approaches such as the strength degradation of rock (SDR) masses and continuous damage mechanics (CDM), this investigation shows that the new advanced methods significantly enhance the quality and accuracy of simulations. Moreover, SDR models address the limitations of classical plasticity models by incorporating nonlinear stress paths and irreversible stress changes, while CDM offers detailed insights into microstructural defect progression. These advancements allow for more accurate and practical predictions of long-term stability in geomechanical engineering tailored to specific requirements of each project.

Shear behavior of BFRP anchor-jointed rock mass considering inclination angle and pre-tension

In this study, the double-shear tests of anchor-jointed rock mass were carried out to investigate the mechanical behaviors of Basalt Fiber Reinforced Polymer (BFRP) anchor bolt in the joint surface. A total of 21 specimens were prepared and two main influence factors were considered: the inclination angle of the BFRP anchor bolt (0°, 5°, 10°, 15°) and the load level of pre-tension (36, 72, 108 kN). The parametric study and design recommendations for the shear capacity of the joint surface under the BFRP anchor bolt were carried out. The results showed that a larger inclination angle or pre-tension force can improve the shear capacity and reduce the relative dislocation at the joint surface. When the inclination angle was 15°, the ultimate bearing capacity was increased by 51 % compared with that when the inclination angle was 0°. When the pre-tension force was 108 kN, the ultimate bearing capacity was increased by 142 % compared with that without pre-tension force. The analytical study found that Ranjbarnia's formulations were reliable to calculate of the shear capacity at the joint surface. The predicted results showed good agreement with the experimental results, with the relative errors less than 5 % for all groups. Then, by using this analytical model, the design parameters of the BFRP anchor bolt were tentatively determined at demanded shear capacities of 200 kN and 350 kN, respectively.

Effects of joint persistence and testing conditions on cyclic shear behavior of en-echelon joints under CNS conditions

Cyclic shear tests on rock joints serve as a practical strategy for understanding the shear behavior of jointed rock masses under seismic conditions. We explored the cyclic shear behavior of en-echelon and how joint persistence and test conditions (initial normal stress, normal stiffness, shear velocity, and cyclic distance) influence it through cyclic shear tests under CNS conditions. The results revealed a through-going shear zone induced by cyclic loads, characterized by abrasive rupture surfaces and brecciated material. Key findings included that increased joint persistence enlarged and smoothened the shear zone, while increased initial normal stress and cyclic distance, and decreased normal stiffness and shear velocity, diminished and roughened the brecciated material. Shear strength decreased across shear cycles, with the most significant reduction in the initial shear cycle. After ten cycles, the shear strength damage factor D varied from 0.785 to 0.909. Shear strength degradation was particularly sensitive to normal stiffness and cyclic distance. Low joint persistence, high initial normal stress, high normal stiffness, slow shear velocity, and large cyclic distance were the most destabilizing combinations. Cyclic loads significantly compressed en-echelon joints, with compressibility highly dependent on normal stress and stiffness. The frictional coefficient initially declined and then increased under a rising cycle number. This work provides crucial insights for understanding and predicting the mechanical response of en-echelon joints under seismic conditions.

A fuzzy K-Means algorithm based on Fisher distribution for the identification of rock discontinuity sets

Rock discontinuities significantly impact the mechanical and hydraulic behavior of rock masses. A crucial aspect of rock engineering involves classifying discontinuities with similar orientations into groups. For this purpose, clustering algorithms, such as K-Means and Fuzzy K-Means (FKM), have been employed. However, the outcomes of these algorithms are influenced by the selection of initial cluster centers. This paper proposes an improved FKM algorithm to automatically identify rock discontinuity sets based on the Fisher distribution (FFKM). The algorithm uses the Fisher distribution to generate and select appropriate initial cluster centers. The performance of FFKM was initially validated using a published data set, and its results were compared with other clustering methods commonly used for grouping discontinuities. Results demonstrated the superior performance of FFKM over the FKM algorithm, comparable to other methods. Subsequently, the proposed algorithm was employed to analyze a fracture data set sampled at an open-pit iron mine. The FFKM facilitated identifying the correct number of sets and produced results consistent with the fracture sets observed in the field. Finally, the algorithm was verified using an artificial discontinuity data set, and the results demonstrated that the method correctly identified the number of sets and provided discontinuity sets similar to the original data set. The FFKM algorithm offers significant advantages: it maintains the essential characteristics of the FKM algorithm, effectively addresses the challenge of selecting suitable initial cluster centers, requires only the expected number of discontinuity sets as an input parameter, processes data within an acceptable computation time, serves as a tool for defining the number of discontinuity sets, and mitigates the drawbacks of the stereographic projection method.

Study on the macro-mechanical behavior and micro-structure evolution law of broken rock mass under triaxial compression

Study on the macro-mechanical behavior and micro-structure evolution law of broken rock mass under triaxial compression

Unraveling time-dependent roof stability dynamics in Iran's coal mines through laboratory-based rock displacement testing

Roof stability is a critical concern in coal mines, as the potential for roof collapse poses a significant risk to miners' safety and productivity. Roof stability is heavily influenced by the time-dependent properties of the rock mass above the workings. This study uses rock displacement testing to examine the effect of time-dependent properties on roof stability and resistance reduction in coal mines in Iran. The study employed laboratory-based rock displacement tests on samples collected from coal mines in Iran, subjected to varying stress levels over time, to simulate the gradual deterioration of rock mass strength in underground mining conditions. The samples were monitored for displacement under these conditions to quantify the reduction in roof resistance over time and assess its effect on roof stability. The study found that areas with high stress at equilibrium gradually fail with time, and the stress transfers from the failure zone into deeper solid rock. The results demonstrate that varying viscous parameters can lead to different relaxation behaviors and stress distribution. Furthermore, incorporating strength degradation into numerical simulation can capture the failure under creep conditions and improve the accuracy of predicting time-dependent roof failure. This research aims to enhance safety measures and reduce the risk of collapses by investigating the time-dependent properties of roof stability through rock displacement testing in Iran's coal mines. The study's innovative approach uses numerical simulation based on the viscoelastic-plastic model to simulate the time-dependent behavior of the rock and incorporate strength degradation into the simulation. The results provide valuable insights into the time-dependent behavior of rock mass in coal mines in Iran and contribute to developing strategies for improving roof stability and lessening the chance of roof collapses. The instantaneous elastic strain was 4.35 × 10–4, and creep simulation was activated to run for a time equivalent to 2 × 106 s.

Stress Characteristics and Mechanical Behavior of Rock Masses with an Opening under Complex Deep Underground Stress Conditions

In this study, the impact of principal stress states on the stress characteristics and initial failure of the rock mass surrounding a three-center arch opening was investigated using complex variable function methods and Discrete Element Method (DEM) numerical modeling. First, the mapping function of the opening was determined using the trigonometric interpolation method, and the influence of the number of terms in the mapping function on its accuracy was revealed. Based on this, the far-field stress state of the underground rock mass was characterized by the ratio of the minimum to maximum principal stress (λ) and the angle (β) between the principal stress and the vertical direction. This stress state was then converted into normal and shear stresses. Using complex variable function theory, the stress characteristics at the boundary of the opening under different stress states were analyzed. Finally, DEM numerical modeling was employed to study the initial failure characteristics at the boundary of the opening and its relationship with the stress distribution. The results indicate that the lateral pressure coefficient significantly affects the stability of the opening by influencing stress concentration around the surrounding rock. Low lateral pressure coefficients lead to tensile stress concentration at the boundary perpendicular to the maximum principal stress. As the coefficient increases, tensile stress decreases, and compressive stress areas expand. While the principal stress direction has a minor effect on stress concentration, it notably impacts stress distribution at the boundary. When λ < 1.0 and β = 45°, stress distribution asymmetry is most pronounced, with the highest compressive stress. The early failure distribution aligns with stress concentration areas, validating the use of stress analysis in predicting opening stability and failure characteristics.

Fractal evolution characteristics of fracture meso-damage in uniaxial compression rock masses using bonded block model

With regard to deep mining in metal mines, an investigation into the failure mode of deep fractured rock masses and their corresponding acoustic emission signal characteristics is conducted via uniaxial compression tests. Subsequently, a fractal damage renormalization group mechanical model is developed to explain the behavior of those fractured rock masses. Employing the bonded block model (BBM) numerical simulation method, fracture process in synthetic rock samples is analyzed, thereby validating the efficacy of the mechanical model. The numerical simulations highlight the critical role of fractures expansion in underlying the deterioration of rock mass strength. As the peak load decreases, the fracture fractal dimension increases, leading to a significant 14.2% reduction in compressive strength accompanied by an approximate 8.7% rise in average fracture fractal dimension. A comparative analysis of tetrahedral and voronoi block synthetic rock samples reveals the tetrahedral block samples exhibit a superior ability to depict the fracture behavior of fractured rock masses. Specifically, they offer a more accurate simulation of acoustic emission characteristics and failure modes. Furthermore, variations in the fracture fractal dimension with respect to the hole defect’s position are observed, with the maximum value occurring along the vertical axis of the hole defect. This observation underscores the potential utility of visually monitoring deep rock fracture dynamics as an effective mean for quantitatively evaluating fracture damage and strength degradation in deep rock formations.

Determination of deformation modulus and characterization of anisotropic behavior of blocky rock masses

Determination of deformation modulus and characterization of anisotropic behavior of blocky rock masses

Influence of 3D joint roughness on fracture behaviours of rock mass subjected to compression

Jointed rock mass with internal rough structural planes have a significant influence on the mechanical behaviour and damage mode of the rock mass. Based on the Weierstrass-Mandelbrot (W-M) fractal function, a numerical model containing a single rough fracture was established. Multi-parameters sensitivity analysis was carried out to determine the key parameters affecting the mechanical properties. By conducting experimental research on the single fracture inclination, fractal dimension, and surface precision of numerical model uniaxial compression, the influence laws of different parameters and the model cracking process were revealed. The results show that the cracks go through the process of sprouting from the single fracture and expanding to both sides, and the expansion of shear cracks in the structural plane reflects the roughness pattern of the structural plane. The higher the roughness was, the lower the peak uniaxial compression strength of the model was. The change of the parameters influenced the macroscopic cracking morphology of the model, and the model finally showed the destruction along the diagonal, the V-shaped destruction, the X-like destruction, and the X-type destruction. The results have an outstanding reference value for further researches on the application of three-dimensional rough discrete fracture networks.

Investigate on the mechanical properties and microscopic three-dimensional morphology of rock failure surfaces under different stress states

The macro and micro morphology of rock failure surfaces play crucial roles in determining the rock mechanical and seepage properties. The morphology of unloaded deep rock failure surfaces exhibits significant variability and complexity. Surface roughness is closely linked to both shear strength and crack seepage behavior. Understanding these morphology parameters is vital for comprehending the mechanical behavior and seepage characteristics of rock masses. In this study, three-dimensional optical scanning technology was employed to analyze the micromorphological properties of limestone and sandstone failure surfaces under varying stress conditions. Line and surface roughness characteristics of different rock failure surfaces were then determined. Our findings reveal a critical confining pressure value (12 MPa) that influences the damage features of Ordovician limestone failure surfaces. With increasing confining pressure, pore depth and crack formation connecting the pores also increase. Beyond the critical confining pressure, the mesoscopic roughness of the failure surface decreases, and the range of interval-distributed pore roughness diminishes. Additionally, we conducted a detailed investigation into the water conductivity properties of rocks under different stress states using Barton's joint roughness coefficient (JRC) index and rock fractal theory. The roughness features of rock failure surfaces were classified into three categories based on mesoscopic pore and crack undulation forms: straight, wavy, and jagged. We also observed significant confining pressure effects on limestone and sandstone, which exceeding the critical confining pressure led to increased water conductivity in both rocks, albeit through different mechanisms. While sandstone exhibits fissures running across it, limestone shows shear abrasion holes. Beyond the critical confining pressure, the rock failure surface becomes smoother, leading to decreased water flow blocking capacity. The fractal dimension of Ordovician limestone increases significantly under critical confining pressure, leading to a more complex mesoscopic crack extension route.