1. Introduction
The evaluation of rock brittleness holds significant research significance in the field of academic studies [
1,
2]. Understanding the brittleness characteristics of rocks is crucial for various applications, such as geotechnical engineering, rock mechanics, and petroleum reservoir characterization [
3,
4]. Evaluating rock brittleness helps in predicting rock failure behavior, designing safe and stable structures, optimizing drilling and hydraulic fracturing operations, and assessing the potential for induced seismicity [
5,
6,
7].
Currently, numerous studies have been conducted to investigate rock brittleness evaluation methods [
2]. These studies have proposed various approaches, including empirical indexes, laboratory testing, and numerical modeling techniques [
2,
3,
8]. Empirical indexes often rely on simple parameters, such as the ratio of compressive to tensile strength or the strength ratio index, to estimate rock brittleness [
9,
10]. Laboratory testing involves conducting experiments, such as uniaxial or triaxial compression tests, to analyze the stress-strain behavior and failure patterns of rocks [
11,
12,
13]. Numerical modeling techniques utilize advanced computational methods, such as finite element analysis or discrete element modeling, to simulate rock behavior and predict brittleness. Brittleness indexes based on mineral composition vary with the classification criteria of brittle minerals and do not consider the influence of rock diagenesis and strength components on rock brittleness [
1,
3,
14]. Although brittleness indexes based on mechanical properties take into account the influence of rock mass strength characteristics, they are not suitable for evaluating the brittleness characteristics of different types of rocks and complex stress conditions in the rock mass [
5,
6]. Brittleness indexes based on the characteristics of the complete stress-strain curve can reflect the influence of crack propagation on rock failure deformation during the rock’s damage process but fail to explain the strength evolution process of the rock [
7,
12]. Brittleness indexes based on hardness testing are primarily used for the brittleness evaluation of homogeneous materials, and the heterogeneity and anisotropy of rocks inevitably result in significant deviations in brittleness evaluation in such cases [
7]. Despite the progress made in rock brittleness evaluation, challenges still exist. The complex nature of rock behavior, including heterogeneity, anisotropy, and nonlinear stress-strain relationships, poses difficulties in accurately assessing brittleness [
15,
16]. Additionally, the lack of standardized brittleness evaluation criteria and the need for integrating multiple parameters for comprehensive assessments remain areas of active research.
In light of these research gaps, this study aims to contribute to the field of rock brittleness evaluation by proposing a novel method that incorporates energy considerations into the assessment. By analyzing the stress-strain curve characteristics before and after the peak stress, this method provides a more physically meaningful approach to quantifying rock brittleness. The findings of this study will enhance our understanding of rock failure mechanisms and provide valuable insights for engineering and geoscience applications.
2. The Existing Methods for Determining Rock Brittleness
Brittleness evaluation plays a significant role in assessing the mechanical behavior of rocks for various engineering applications [
2]. This section aims to discuss the current state of the research of rock brittleness assessment and critically analyze the advantages and disadvantages of different approaches based on mineral composition, strain, stress, elastic parameter, and energy.
Table 1 shows the existing methods for determining the rock brittleness.
- (a)
Definition of Brittleness Index based on Stress
The stress-based approach defines rock brittleness in terms of stress levels at failure or the stress intensity required for crack propagation. This method allows for the assessment of the rock’s resistance to fracturing under different loading conditions. The advantage of stress-based definitions is their direct relationship with fracture initiation and propagation. However, accurately determining stress levels in rocks can be challenging, particularly in heterogeneous formations.
- (b)
Definition of Brittleness Index based on Strain
Another approach to evaluating rock brittleness is through strain-based definitions. This method focuses on quantifying the deformation characteristics of rocks under stress. By considering the strain at failure or the strain energy release, this approach provides insights into the rock’s ability to undergo plastic deformation before failure. The advantage of strain-based definitions is their ability to capture the ductility or plasticity of rocks [
17,
18]. However, their drawback lies in the complexity of strain measurements and the need for advanced testing equipment.
- (c)
Definition of Brittleness Index Based on Elastic Parameters
Another approach involves defining brittleness indexes using elastic deformation parameters, such as Young’s modulus and Poisson’s ratio. These parameters reflect the rock’s ability to undergo elastic deformation without fracturing. This approach provides a more direct connection to rock mechanics properties and can be determined through laboratory testing or empirical correlations. However, it neglects the influence of plastic deformation and post-peak behavior, limiting its ability to capture the complete brittleness spectrum.
- (d)
Definition of Brittleness Index based on Energy
The energy-based definition of rock brittleness focuses on quantifying the fracture energy or the energy required to propagate cracks within rocks [
19,
20,
21,
22]. This approach considers the relationship between the energy dissipated during fracture and the rock’s brittleness characteristics. The advantage of energy-based definitions is their ability to capture the rock’s fracture toughness and resistance to crack propagation. However, accurately measuring fracture energy necessitates sophisticated experimental setups and can be time-consuming.
- (e)
Definition of Brittleness Index based on Rock Mineral Composition
One approach to defining rock brittleness is by considering the mineral composition. This method involves characterizing the rock’s brittleness index based on the proportions and mechanical properties of its constituent minerals. The advantage of this approach is its ability to capture the influence of mineralogy on rock behavior. However, its major limitation lies in the difficulty of accurately quantifying the mineral composition and establishing consistent correlations between mineralogy and brittleness.
Based on
Figure 1, the shortcomings of the current brittleness evaluation method based on stress-strain curves can be identified. For lines 1 and 2, they exhibit the same stress drop after the peak but have different slopes before the peak. If only the post-peak portion of the stress-strain curve is considered, it would inevitably lead to the calculation of the same brittleness index. Clearly, this type of brittleness index is not suitable. For lines 1 and 3, they have the same slope before the peak but different peak strengths and different stress drops after the peak. If the brittleness index is calculated based on the pre-peak elastic modulus, it neglects the influence of the post-peak curve on rock brittleness. It is evident that this type of brittleness index is not appropriate. For lines 2 and 4, they have the same slope before the peak, the same stress drop after the peak, but different peak strengths. If only the post-peak portion of the stress-strain curve and the pre-peak elastic deformation are considered, it overlooks the influence of the rock’s inherent peak strength. Clearly, this type of brittleness index is not suitable. For lines 3 and 4, they have different slopes before the peak, different stress drops after the peak, but the same peak strength and the same post-peak slope. Most brittleness indexes also fail to reflect this characteristic. Therefore, it is necessary to propose a brittleness index that can simultaneously reflect the pre-peak and post-peak curve characteristics and capture the elastic and plastic deformation features during the rock failure process.
3. A Method to Evaluate the Brittleness in Consideration of Plastic Deformation in Pre-Peak and Failure Energy in Post-Peak
In the evaluation of rock brittleness, it is essential to consider both pre-peak and post-peak stress-strain curves. This is because the behavior of rocks during the pre-peak and post-peak stages provides valuable insights into their brittleness characteristics.
During the pre-peak stage, the rock undergoes elastic deformation and initial crack propagation. By analyzing the stress-strain curve in this stage, we can determine the elastic modulus, which represents the rock’s stiffness and ability to withstand deformation without fracturing. The elastic modulus is a crucial parameter in defining the rock’s brittleness, as a higher modulus indicates a more brittle behavior. The inclusion of pre-peak plastic deformation in the brittle evaluation method is necessitated by the recognition that rocks often exhibit significant ductility before reaching their peak strength. By quantifying the extent of plastic deformation prior to failure, this method captures the rock’s ability to deform without fracturing. This information is essential for engineering design, as it helps identify potential failure modes and select appropriate mitigation measures.
After reaching the peak stress, the rock experiences post-peak behavior characterized by progressive failure and plastic deformation. By examining the stress-strain curve in this stage, we can assess the extent of plastic deformation, energy absorption, and the ability to undergo ductile deformation. These parameters influence the rock’s brittleness, as a higher degree of plastic deformation and energy absorption indicates a more ductile behavior. The consideration of post-peak failure energy in the brittle evaluation method is crucial as it provides insights into the rock’s fracture toughness and energy dissipation capacity. By quantifying the energy required for crack propagation and failure, this approach can assess the rock’s resistance to catastrophic failure and its ability to absorb and dissipate energy. Such information is vital for optimizing engineering designs and ensuring the stability of rock structures.
Considering both pre-peak and post-peak stress-strain curves allows for a comprehensive evaluation of rock brittleness. It provides a holistic understanding of the rock’s mechanical response to stress, including its initial stiffness, crack propagation, and ability to undergo plastic deformation. By analyzing the entire stress-strain curve, we can capture the complete spectrum of brittleness and understand the rock’s failure mechanisms.
Based on the definitions provided by existing scholars, this study considers brittleness and ductility as opposing characteristics [
17,
18]. The strain occurring in the rock mass prior to failure can be divided into recoverable elastic strain and irrecoverable plastic strain induced by crack propagation. The higher the degree of plastic deformation, the lower the brittleness of the rock mass, while the brittleness is independent of the amount of elastic deformation during the failure process. Therefore, if the amount of plastic strain occurring prior to failure is minimal, the rock mass is considered more brittle. As shown in
Figure 2, the brown area represents the stored elastic strain energy prior to failure, and the slope of the line dividing the yellow area and the brown area corresponds to the elastic modulus. Hence, 2
We/
E can be used to obtain the elastic strain occurring prior to failure, while the deformation prior to failure consists of plastic deformation and elastic deformation. Therefore, subtracting the elastic deformation from the total strain prior to failure provides the plastic deformation prior to failure, thus evaluating the brittleness value of the rock mass. The formula for the pre-peak brittleness index is given as follows:
where
is the strain before peak strength,
is the elastic strain energy, and
E is the elastic modulus. Herein, the log symbol represents the log for the decadic logarithm.
Analyzing the brittleness of rock masses solely from an energy perspective prior to the peak does not allow for distinguishing between the energy required to overcome cohesive forces and the energy needed for friction and heat generation. Therefore, in this study, the brittleness of the rock mass after the peak is assessed based on plastic strain. This approach provides a comprehensive evaluation considering the specific energy consumed during crack propagation and accounts for the rock mass’s post-peak behavior.
The post-peak brittleness index represents the total strain energy required for failure after reaching the peak stress. The difference in crack propagation between brittle and non-brittle rocks during the failure process using a schematic diagram is demonstrated. It can be observed that brittle materials require less energy for crack propagation, while non-brittle materials require energy to sustain crack growth. Therefore, the lower the energy required for crack propagation in the post-peak stage, the higher the brittleness of the rock mass. As shown in
Figure 2, the energy absorbed by the rock mass prior to the peak is divided into dissipation energy (yellow area) used for crack propagation, overcoming cohesive forces, friction, and heat dissipation, and the stored elastic strain energy (brown area) of the rock mass. Some of the stored elastic energy is utilized for crack propagation after the peak, while the remaining portion is residual elastic energy (green area). The energy for rock mass failure after the peak includes not only the stored elastic energy prior to the peak but also the external energy input after the peak (wall shape). Therefore, the total energy for post-peak rock mass failure is given by (brown area + green area − green area), which can be expressed as
We +
Wout −
Wre. However,
We +
Wout −
Wre is inversely proportional to the brittleness of the rock mass. Therefore, the reciprocal form is used to express the brittleness of the rock mass.
where
is the stored elastic strain energy before the peak,
is the external energy input after the peak, and
is the residual elastic strain energy.
In this study, a logarithmic function is employed to reduce the numerical disparity between the two components of the brittleness index, thus avoiding the influence of one component’s magnitude on the magnitude of the other component. For the coupling of the two brittleness index components, a multiplication form is selected, following the same principle as mentioned earlier. Additionally, the two components of the brittleness index describe the characteristics of stress-strain curve at different stages. Therefore, based on the energy parameters at the pre-peak and post-peak stages, the proposed brittleness index in this study is presented as follows:
5. Comparison of Brittleness Index with Other Existing Methods
As shown in
Figure 7, the changes in brittleness under different confining pressures for red-bed sandstone under saturated conditions and in its natural state were evaluated using B7, B9, and B17, which were compared to the brittleness evaluation method proposed in this study. Brittleness indexes B7 and B9 based on stress-strain parameters are selected to evaluate rock brittleness, while B17 is a brittleness index proposed based on energy parameters. The above three indexes have been validated by numerous cases, demonstrating their ability to evaluate rock brittleness [
2,
10]. Therefore, the rationality of the proposed method is verified through various different definitions of brittleness indexes. Due to the different orders of magnitude of the various brittleness indexes, the brittleness indexes of each method were standardized to facilitate the presentation of their respective trends. In
Figure 7, under natural conditions, the proposed brittleness index in this study exhibits the maximum values at all confining pressures, with the lowest value observed for B7. Furthermore, the magnitude of decrease in values for B7 and B9 diminishes as the confining pressure increases. The four brittleness indexes consistently characterize the brittleness of the red-bed sandstone under natural conditions. Under saturated conditions, the brittleness indexes exhibit different patterns of variation. The values of the proposed brittleness index, in comparison with the other indexes, fall within the intermediate range, and their trends of change are generally consistent.
Compared to B7 and B9, the advantages of the brittleness evaluation proposed in this study are as follows: the proposed method describes the stress-strain curve characteristics before and after the peak stress, and compared to B7 and B9, which use the stress increment after the peak stress and the slope of the curve to describe the brittleness of rock masses, the use of rock mass energy in this study for evaluating brittleness is more physically meaningful. This study argues that the stress growth rate and stress-strain rate only indicate the strength characteristics of the rock mass, without a direct relationship with brittleness. However, analyzing the characteristics of rock mass failure from an energy perspective can better reflect the brittleness characteristics of the rock mass. Compared to the brittleness evaluation proposed in this study, B17 lacks the characteristics after the peak stress and uses the proportion of elastic energy to evaluate the brittleness of the rock mass, which would give an advantage to highly elastic rock masses. Therefore, the study suggests that if energy is used to evaluate the brittleness of rock masses, it should be based on the amount of energy consumed by the rock mass in overcoming friction and heat release.
As shown in
Figure 8, the changes in brittleness under different confining pressures for cyan sandstone under freezing-thawing conditions and in its natural state were evaluated using B7, B9, and B17, which were compared to the brittleness evaluation method proposed in this study. Due to the different orders of magnitude of the various brittleness indexes, the brittleness indexes of each method were standardized to facilitate the presentation of their respective trends. Under natural conditions, the three brittleness indexes, B7, B9, and Bnew, exhibit a similar trend in characterizing the brittleness of the cyan sandstone. They gradually decrease with increasing confining pressure, and the value of Bnew is larger than B9 and B7. However, the trend of B17 is essentially larger than the previous three indexes, as it decreases gradually with increasing confining pressure at a relatively uniform rate. Under freezing-thawing conditions, the variation patterns of the three brittleness indexes, B7, B9, and Bnew, are similar, and the value of Bnew also falls between B9 and B7. At lower confining pressures, the three brittleness indexes show a greater rate of decrease, and as the confining pressure increases, their rate of decrease tends to stabilize. Similarly, the trend of B17 is essentially opposite to the previous three indexes.
In
Figure 9, the
x-axis represents the normalized brittleness indexes, where 1 represents the normalized brittleness index B7′, 2 represents B9′, 3 represents B17′, and 4 represents Bnew’. Comparing the normalized brittleness indexes reflects the variations in rock brittleness under different conditions. From the graph, it can be observed that, for the red-bed sandstone, under natural conditions, the brittleness gradually weakens with increasing confining pressure, indicating that confining pressure has an inhibitory effect on the brittleness of red-bed sandstone. However, under the influence of water–rock interaction, the brittleness initially increases and then decreases with increasing confining pressure, suggesting that the brittleness of red-bed sandstone changes after water–rock interaction. This change is more pronounced under isotropic conditions, and the brittleness enhances under anisotropic conditions.
As for the cyan sandstone, it is evident that the change trend of the brittleness index Bnew’ is consitent with the other three indexes. For example, at a confining pressure of 3 MPa, after freezing-thawing action, the radii of the spheres representing B7′, B9′, and Bnew’ are significantly smaller than those under natural conditions, indicating that freezing-thawing action accelerates the rate of brittle-ductile transition in the cyan sandstone.
6. Discussion
6.1. Rock Failure Mechanism from the Energy Evolution
The transformation of rock brittleness and ductility depends on factors such as stress level, physical properties of the rock (including crack density and rock toughness), mechanical properties of the rock (including tensile strength and shear strength), and loading rate. These factors collectively determine the rock’s failure behavior and mechanical response under stress. According to the proposed brittleness index in this study, during the evolution from absolute ductility to absolute brittleness, both pre-peak dissipated energy and post-peak fracture energy gradually decrease. As shown in
Figure 10, two characteristic points exist in this process: the ductility characteristic point (E > −M = Q) and the moderate brittleness characteristic point (E = −M > Q). Based on these two characteristic points, the brittleness characteristics of the rock during the entire failure process can be classified as follows: absolute ductility (E > Q = M = 0), weak brittleness (E > −M > Q), strong brittleness (−M > E > Q), and extremely strong brittleness (E > M > Q).
The understanding of rock failure mechanisms is fundamental to evaluating rock brittleness and implementing appropriate engineering interventions [
23,
24]. This section focuses on elucidating the mechanisms of rock failure from the perspective of energy evolution. By examining energy redistribution, dissipation, and release during the progressive failure process, a comprehensive understanding of rock brittleness can be achieved. The discussion emphasizes the significance of energy-based analysis and its implications for engineering applications.
- (a)
Energy Redistribution during Rock Failure
Rock failure is often accompanied by energy redistribution within the rock mass [
25]. This redistribution occurs due to stress concentration and the initiation and propagation of cracks. Energy-based analysis allows for quantifying the redistribution of potential and kinetic energy, providing insights into the spatial and temporal evolution of failure. By understanding energy redistribution, engineers can anticipate potential failure zones, optimize support systems, and design effective reinforcement measures.
- (b)
Energy Dissipation during Rock Failure
Energy dissipation plays a crucial role in rock failure processes. As cracks propagate and interact, energy is dissipated through various mechanisms, such as friction, microcracking, and fracture surface roughness. Energy-based analysis enables the quantification of dissipated energy, helping to assess the brittleness of rocks. The measurement of energy dissipation provides valuable information for designing energy-absorbing structures and selecting suitable materials to enhance rock stability.
- (c)
Energy Release during Catastrophic Rock Failure
Catastrophic failure events, such as rockbursts and collapses, involve the sudden release of stored energy. Energy-based analysis allows for understanding the magnitude and spatial distribution of released energy during such events. This information is crucial for assessing the potential impact on surrounding structures and developing effective hazard mitigation strategies. By quantifying energy release, engineers can optimize support systems, implement real-time monitoring, and design robust containment measures.
6.2. Application for the Rock Brittleness
To further the rationality of the proposed brittleness indexes for different rock types, the datasets of uniaxial and triaxial tests on two rock types, including granite and coal specimen, are collected from the literature [
26,
27]. The brittleness for the two rock types is calculated using the proposed formula. Herein, the brittleness index B7 is selected to verify the accuracy of the proposed brittleness index, as shown in
Figure 11, both B7 and Bnew have the same variation trend with the confining pressure, characterizing by gradual decline. Evidently, the Bnew conforms to the law of exponential decline, which better reflects the brittleness characteristics of the two rocks. Therefore, the proposed brittleness index can be used to evaluate the other rock types.
A brittleness index is proposed in this study, which describes the characteristics of stress-strain curves at different stages by considering the pre-peak plastic deformation and post-peak failure energy of rock mass. This leads to the development of a new brittleness evaluation method that offers the following advantages: (a) Consideration of pre-peak plastic deformation: Traditional brittleness evaluation indices often overlook the pre-peak plastic deformation of rock mass. However, this method introduces a pre-peak brittleness index factor that fully takes into account the plastic deformation characteristics of rocks during the loading process. This aids in a more accurate assessment of rock brittleness. (b) Consideration of post-peak failure energy: The method also introduces a post-peak brittleness index factor, considering the total energy consumed by the rock mass during the post-peak failure process. This reflects the failure characteristics of rocks in the post-peak stage, providing more comprehensive information for evaluating rock brittleness. (c) Clear physical significance: The method is based on the analysis of stress-strain curve characteristics and provides a more physically meaningful approach to quantify rock brittleness through energy considerations. Compared to traditional empirical indices or simple parameters, this energy-based evaluation method better reflects the mechanical properties and engineering behavior of rocks.
However, the method also has limitations and drawbacks, including requirement for experimental data support. The application of this method requires experimental testing to obtain relevant parameters such as stress-strain curves and failure energy of rocks. This may entail significant experimental workload and costs, and could be subject to limitations imposed by experimental conditions.
Rock brittleness evaluation methods are essential tools in rock engineering, offering valuable information for a wide range of engineering applications [
28]. This section aims to provide an overview of the engineering applications of these methods, emphasizing their importance in optimizing engineering designs, assessing stability, and predicting rock failure. The discussion also addresses the challenges associated with implementing these methods in practical engineering scenarios.
- (a)
Optimization of Engineering Designs
Rock brittleness evaluation methods are instrumental in the optimization of engineering designs, particularly in areas such as tunneling, mining, and geotechnical engineering. By assessing the brittleness characteristics of rocks, engineers can make informed decisions about support systems, excavation methods, and blasting techniques. This leads to more efficient and cost-effective designs that minimize the risk of rock failures and ensure the long-term stability of structures.
- (b)
Stability Analysis of Rock Structures
The evaluation of rock brittleness is crucial for the stability analysis of rock structures, including slopes, foundations, and underground excavations. By considering the brittleness indexes, engineers can assess the potential for brittle failure and identify areas prone to rock fracturing. This information helps in designing appropriate reinforcement measures and slope stabilization techniques, minimizing the risk of catastrophic failures and ensuring the safety of infrastructure [
29].
- (c)
Prediction of Rock Failure
Rock brittleness evaluation methods aid in predicting rock failure, providing early warnings and allowing for proactive measures to mitigate potential hazards. By quantifying the brittleness characteristics of rocks, engineers can assess their susceptibility to failure under different loading conditions. This information is vital for designing monitoring systems, implementing preventive measures, and developing reliable risk management strategies [
30].
- (d)
Challenges and Future Directions
The practical implementation of rock brittleness evaluation methods faces several challenges, including the accurate quantification of brittleness indexes, the heterogeneity of rock formations, and the need for advanced testing equipment. Future research should focus on developing standardized testing protocols, integrating multiple evaluation methods, and improving the applicability of these methods to diverse geological conditions. Additionally, advancements in non-destructive testing techniques and numerical modeling can enhance the accuracy and efficiency of rock brittleness assessments.