AMathematical Model for Calculating the “Brittleness-Ductility” Drop Coefficient of Sandstone in Mining Zones

Sandstone is different from pure brittle rocks such as marble and granite. Its brittleness weakens and then increases with increasing confining pressure. *is paper establishes a calculation model of the “brittleness-ductility” energy drop coefficient that characterizes the whole process of sandstone failure from the perspective of energy. *e obtained result was compared with the result of traditional calculation model. *e results show that the established model is more consistent with the experimental phenomenon and can describe the whole process of sandstone failure more precisely at a macro level. Moreover, the model is more practical since it avoids the tedious integration calculation process. *e model and tests have extended the applicable range and applied rock types of the energy drop coefficient and enriched the knowledge system of elastoplastic mechanics. At the same time, they carry important significance for environmental restoration and geological disaster prevention and control in sedimentary rock mining zones, which can promote the sustainable development of related industries in the mining zones.


Introduction
China is one of the few countries in the world which use coal as the main energy source. Its coal output has exceeded onethird of the world's total [1,2]. Resource destruction, environmental damage, and industrial accidents related to coal production are very prominent [3,4]. In particular, during the construction of new mine and waste mine backfilling, frequent mine collapses threaten the safety of workers and their property, seriously damage the ecological environment and groundwater quality, or even trigger more severe geographic hazards [5][6][7]. e sustainable development of coal resources and related industries is facing challenges. Sandstone is the most common surrounding rock in coal mines in China [8]. Frequent landslides are often closely related to the failure of sandstone. erefore, a deep understanding of its failure mechanism can ensure the safety of mine construction and exploitation [9,10]. In addition, it helps to prevent the negative impact of mining industry on the surroundings and has far-reaching significance for environmental restoration and disaster prevention and control in the mining zones.
In traditional elastic-plastic mechanics theory, the stressstrain curve is usually used to describe the process of material deformation and failure, thus establishing the constitutive model of the material and the corresponding strength theory. For linear elastic materials, this method can be used to reflect mechanical characteristics of the material accurately. However, due to the anisotropy and inhomogeneity of the material, the stress-strain relationship of the rock is distinctively nonlinear [5,11]. As a result, for stones it is difficult to establish a reasonable strength criterion based on the stress-strain curve alone. In fact, the loading process is a process of exerting external energy input. Before its stress reaches the peak, the rock absorbs and stores the external energy taking advantage of its own elastoplasticity to maintain a stable state [12][13][14]. At the peak intensity, the rock is in a critical state of high energy instability [15]. With the continuous input of external energy, the rock finally releases the stored energy and loses its stability, and the stress quickly drops to the residual strength [16,17]. is means that the rock deformation and failure process is an irreversible process of energy release and dissipation. e ultimate failure of the rock is actually a sudden change in the process of energy dissipation. erefore, to describe the deformation and failure process of rocks from the perspective of energy can better reflect the physical nature of its failure mechanism.
In recent years, many researchers have studied the failure mechanism of brittle rocks based on energy analysis. Zhao et al. demonstrated in detail that the failure mode and energy dissipation process change with the anisotropy, inhomogeneity, and stress state of rocks and proposed the principle of minimum energy for dynamic failure based on the findings [18]. Yang et al. performed triaxial tests on marble rocks by using servo machines to investigate the deformation and failure process of rock samples and their energy characteristics [19]. Su et al. studied the patterns of plastic deformation, uniaxial compression strength, average modulus, energy dissipation, and failure characteristics of rock samples during triaxial tests [20]. Xie et al. proposed failure criterion for stratified rock based on distortion energy and generalized volume expansion potential energy, which was in good agreement with the results of the double pressure tests. e established rock damage evolution equation was verified through experiments, and related parameter indicators were determined [21]. Zuo et al. analyzed the shortcomings of the stress drop coefficient of brittle rocks and put forward a new characterization parameter of rock failure-the energy drop coefficient [22]. Song et al. analyzed the applicability and limitation of the energy drop coefficient and further considered the proportion of releasable elastic energy in the total energy before peak stage to present an improved energy drop coefficient [23]. Liu et al. presented valency-based topological indices of chemical networks [24,25]. By using some real world data, we performed certain comparative testing to investigate the performance of almost all well-known valencybased indices. In light of this, many scholars have rediscovered the changes in the stress-strain state of rocks under external loads from an energy perspective and analyzed the energy transfer and conversion in the process of rock deformation and failure in order to establish the strength theory using energy changes as failure criteria, which is a more realistic description of the failure mechanism of rocks.
However, there are some strong limitations in current researches. First, these energy threshold values reflecting rock strength usually only take into account the energy evolution characteristics of the postpeak stage and thus cannot reflect the entire process of rock deformation and failure [26]. Second, the calculation of the traditional energy drop coefficient of rocks mainly requires a large number of integral calculations, which is cumbersome and inconvenient. ird, previous studies have used magmatic rocks and metamorphic rocks such as marble and granite as research objects more, and few studies have been conducted on sedimentary rocks such as sandstone [27]. To address these research gaps, this paper obtains the whole stress-strain curves of sandstones under different confining pressures based on triaxial tests and proposes a calculation model of "brittleness-ductility" energy drop coefficient applicable to sandstone. Based on established model, the mechanism of sandstone failure is explained and influence of confining pressure on energy release and dissipation during rock deformation and failure is discussed. e calculation model of "brittleness-ductility" energy drop coefficient proposed in this paper can truly describe the no-linear deformation and failure process of sandstone. e model fully considers the two states before and after the failure of the rock samples and further enriches the physical meaning of the traditional stress-strain curve in describing material failure mechanism, which can provide references for studies on other rocks with similar mechanical characteristics. In addition, the model can efficiently calculate and characterize the mechanical properties of sandstone and provide a reference for the comprehensive development and management of mining zones, which is beneficial to the environmental restoration and geological disaster prevention in these zones. e study can promote the sustainable development of related industries as well.

Energy Drop Coefficient of Rock Based on
Triaxial Tests e rock triaxial test is essentially a simulation of the mechanical properties of the rock in a three-way compression environment. Compared to conventional test methods, triaxial tests are more intricately advanced tests, which can obtain compressive strength, shear strength, elastic modulus, Poisson's ratio, cohesion, and internal friction angle under different confining pressures (or depths) of rocks. Without changing the ambient temperature, the heat exchange between the test piece and the environment can be ignored. erefore, the energy conversion of the test piece during the tests is mainly mechanical energy exchange with the test machine.
In the test, the work done by the testing machine for the axial compression of test piece is where W z is the work done by the testing machine for the axial compression, and it is also the total energy input into the test piece. F z is the axial equivalent effect, d the diameter of the test piece, h the height, σ z the axial stress, and ε z the axial strain. For the triaxial test with constant confining pressure, the work done by the test piece in the circumferential direction on the hydraulic oil is where F r is the circumferential equivalent force, d the diameter of the test piece, h the height, σ z the radial stress, and ε z the radial strain.
From (1) and (2), it can be known that the actual mechanical energy absorbed by the sample W is the difference between W z and W r : 2 Discrete Dynamics in Nature and Society It is difficult to measure W z and W r directly from the tests. In fact, the deformation and failure of the pieces essentially result from the continuous dissipation and rapid release of energy. It means that part of the input mechanical energy is released as elastic energy and the other converted to dissipative energy. erefore, the actual mechanical energy W absorbed when the test piece reaches the peak pressure can be divided into dissipative energy and releasable elastic energy: where W p is dissipative energy and W e is releasable elastic energy. Figure 1 shows the typical stress-strain curve of rocks. W p corresponds to the area enclosed by the OAB and W e corresponds to the area enclosed by the BAC. When the axial strain is ε, according to (1)-(4), the energy absorbed by the test piece is W(ε) and is equal to the work performed by the test machine on the test piece: (5) where V is the volume of the test piece; σ z , ε z , σ r , and ε r are integral variables; ε is the axial strain; and μ is the Poisson's ratio.
When the confining pressure remains unchanged, W e , the elastic energy released by the rock, can be measured by the unloading test. If the elastic modulus of the test piece is E, where E is the elastic modulus of the test piece and ε e is recoverable elastic strain: From (6) and (7), it can be known that When the stress drops from point A to point D, the change of releasable elastic energy is According to (5), when confining pressure is constant, the stress drops from point A to point D, and the externally input mechanical energy ΔW is e ratio of ΔW and ΔW e is used to characterize the energy change during the stress drop when the test piece fails. is ratio H is defined as the energy drop coefficient of the rocks: where σ A and ε A are the peak stress and peak strain, respectively; μ is the Poisson's ratio; σ D and σ D are the residual stress and residual strain, respectively; and E A is the elastic modulus of unloading at the peak. e energy drop coefficient H expresses a proportional relationship between the mechanical energy input into the rock and the releasable elastic energy. It demonstrates that the deformation and failure of the brittle rock are the result of energy dissipation and sudden release. H has a certain correlation with the brittle failure of the rock. e larger the H value, the less obvious the brittle stress drop of the rock, and the stronger the ductility. erefore, the energy drop coefficient H has important reference value for characterizing the degree of rock brittle failure.

Triaxial Test of Sandstone
In order to find out the failure mechanism of sandstone from the energy perspective and to obtain the relevant mechanical parameters, triaxial tests were performed on typical sandstone samples collected from the Hulusu Mine in Inner Mongolia, China. e test pieces were standard cylinder with a diameter of 50 mm and a height of 100 mm. e tests were conducted on the MTS816 rock test system. e device is shown in Figure 2. is system uses a three-axis servo press, which can record the stress-strain curve and thus obtain the data of test piece including residual stress and permanent deformation. e 1-hydraulic system was adopted to load hydraulic oil into 6-triaxial cavity to simulate and control the confining pressure stress state of the test piece. e 3-control system was used to control and operate 9-loading device and to load the test piece [28]. e test adopted the force-displacement mixed control loading system. In the initial stage, different confining pressures were applied to the test piece through force control. Confining pressures were applied to predetermined   Figure 3 shows the triaxial compressive stress-strain curves of sandstone samples under confining pressures of 4, 6, 8, and 10 MPa, respectively.
It can be learned that before the stress reached its peak value, it increased linearly with the rise of the strain, showing relatively good elasticity. When it reached its peak value, the sample went through buckling failure rapidly, exhibiting obvious brittle stress drop phenomenon. After that, an obvious yielding plateau appeared, which indicated that the sandstones still maintained some ductility under the confining pressures. e carrying capacity increased with the rise of confining pressure. As the confining pressure increased from 4 MPa to 8 MPa, the magnitude of stress drop also gradually went up. However, as the confining pressure continued to do so, the magnitude dropped instead. It showed that as the confining pressure increased, the sandstone exhibited the mechanical characteristics of changing from brittleness to ductility and then back to brittleness. is is inconsistent with the characteristic that the brittleness of granite and marble gradually declines with increasing confining pressure. It is difficult to reveal its failure mechanism using traditional energy drop coefficient. Figure 4 shows the characteristics of the failure surfaces of sandstone samples under different confining pressures. It can be concluded that the samples under different confining pressures had only one shear failure surface which changed along with different confining pressures. Under 4 MPa, the macroscopic fracture surfaces of the samples were cracked from the middle, indicating that they had already been completely destroyed under the low confining pressure state, exhibiting distinguished brittleness. As the confining pressure increased, the angles θ1-θ3 between the normal direction of the sandstone failure surface and the horizontal line gradually increased, and the macroscopic failure surfaces became more and more flat, indicating enhanced ductility of the samples. At 10 MPa, the fracture surfaces corresponding to the confining pressures were obviously different from those under other confining pressure states, with θ 3 being the largest angle. e surrounding areas were severely fragmented, showing more severe brittleness. is again showed that sandstone has different mechanical properties from marble and granite rocks.

"Brittleness-Ductility" Energy
Drop Coefficient e traditional calculation method of energy drop coefficient has been well applied in rocks such as marble and granite. Based on the results of triaxial tests, it can be concluded that sandstone exhibits transformative characteristics of "brittleness-ductility-brittleness" with increasing confining pressure, which makes it difficult for the energy drop coefficient of brittle rock to fully reflect the deformation and failure process of sandstone. In addition, the traditional method requires the calculation of a large number of integrals. e complicated and cumbersome process is difficult to be applied in practical engineering. erefore, it is necessary to establish a calculation model that is easy to use and more suitable for sandstone failure characteristics. When the sandstone enters the ductility stage, the elastic modulus E A corresponding to the unloading at the peak in (11) cannot accurately reflect the stress-strain relationship of the samples at this time. When the effect of sandstone damage and defect is not considered, the rock damage factor D is considered to be 0. At this time, according to E D � (1 − D)E A , it can be replaced by the elastic modulus E D corresponding to the residual strength unloading. In Figure 1, AB and DE stand for E A and E D , respectively. In order to simplify the calculation method, the area of the stressstrain curve is used to represent the energy involved in the deformation and failure process. ε D ε A σ z dε z is represented by the area of trapezoid ACDF, which represents the external input of the mechanical energy into the samples after the peak point A. e releasable elastic energy of the samples can be expressed by the difference between SBAC (in red in Figure 1) and SEDF (in yellow in Figure 1). is value represents the relative magnitudes of the ductility and brittleness of the samples-the smaller the value, the higher the ductility. 2μ(ε D − ε A )σ r is represented by W C . erefore, the energy drop coefficient calculation model for sandstone can be expressed as He can better characterize the relative magnitudes of the ductility and brittleness of the rock, which is called the "brittleness-ductility" energy drop coefficient. e coefficient calculation model is more in line with the energy change of the actual failure process of sandstone and can better reflect the energy conversion of the collapse of sandstone at various depths.

Comparison Analysis
e residual energy, elastic energy, work done by rock samples to hydraulic oil under different confining pressures, and the drop coefficients H and He obtained by the two different calculation models according to (11) and (12) are calculated, and the results are shown in Table 1.
It can be learned that as the confining pressure increases, the mechanical energy input of the samples gradually increases; the energy loss caused by the work performed by the samples to the hydraulic oil rises as well; and the releasable elastic energy gradually decreases. e traditional energy drop coefficient of sandstone H and the "brittle-ductile energy drop coefficient" He have basically the same variation tendency. As the confining pressure increases, the energy drop coefficient also gradually goes up, the brittleness exhibited by the failure of the samples weakens, and the ductility enhances.
is proves the correctness of the "brittle-ductility" energy drop coefficient. However, when the confining pressure rises to 10 MPa, the H value of the sandstone continues to increase while He drops. e results show that when the confining pressure rises to 10 MPa, the brittleness enhances. is demonstrates that the calculation of the coefficient He proposed in this paper can better characterize the relative magnitudes of the ductility and brittleness of the rock and thus reflect the physical failure mechanism of sandstone failure more precisely.

Conclusions.
is paper establishes a calculation model of the "brittle-ductility" energy drop coefficient that characterizes the whole process of sandstone deformation and  iscrete Dynamics in Nature and Society 5 failure from the perspective of energy based on the stressstrain curves of sandstone samples obtained from triaxial tests. It also discusses the influence of confining pressure on energy drop during sandstone failure. e following conclusions are reached: (1) According to the fracture surface characteristics and stress-strain curves of the sandstone obtained from the triaxial tests, as the confining pressure increases, the brittleness of the failed samples diminishes and then gradually enhances after reaching the peak. is is different from other rocks such as granite and marble.
(2) Compared with traditional energy drop coefficients, the calculation model of "brittleness-ductility" energy drop coefficient proposed in this paper can accurately describe the whole process of sandstone deformation and failure from an energy perspective. At the same time, the improved model is easy to adopt because it can avoid tedious integration calculations.
is model extends the applicable rock types of the energy drop coefficient and enriches the knowledge system of elastoplastic mechanics. At the same time, a deep understanding of the sandstone failure mechanism can ensure the safety of construction and exploitation in mining zones while reducing the environmental damage. It is of great significance for environmental restoration and geological disaster prevention.

Recommendations.
For future research, the following research directions deserve scholarly attention: is paper only establishes a calculation model for the "brittle-ductility" energy drop coefficient suitable for characterizing the sandstone failure mechanism but does not explain why the tendency of sandstone brittleness to change with confining pressure differs from that of marble and granite. Future research can try to explain the mechanism from the composition and mechanical properties.
According to the genesis, rocks can be divided into magmatic rocks, sedimentary rocks, and metamorphic rocks. In addition to sandstone, there are other common surrounding rocks like mudstone in the mining zones. Difference between different types of rocks is actually very large. e new calculation model proved to be feasible and efficient to describe the failure of sandstone. However, whether it can suit other types of rocks is still unclear.
erefore, more samples could be collected to study the energy drop law. In light of this, the strength thresholds of rocks can be identified by applying correspondent energy coefficients according to rock types to better guide the safe and sustainable production in the mining zones.
Data Availability e data are generated from experiments.

Conflicts of Interest
e authors declare no conflicts of interest.