1 Introduction

As social productivity improves, people’s demand for mineral resources increases, shallow mineral resources are gradually exploited and exhausted, and deep mineral resource mining gradually becomes more common. The deformation and damage caused by cyclic loading of rocks under deep burial or high in-situ stress is becoming more critical (Taheri et al. 2016; Wang et al. 2021; Erarslan 2016; Hoek and Martin 2014; Erarslan 2013). The superimposed evolution of micro-cracks in the rock mass results in macroscopic damage to the surrounding rock mass in deep underground engineering. Therefore, it is crucial to understand the evolution of cracks in rock to ensure the stable state of the surrounding rock mass, especially during cyclic loading.

Recently, experts and scholars in the field of rock mechanics have intensively studied the evolution characteristics of rock fracture under cyclic loading (Ma et al. 2013; Liu et al. 2022; Qiu et al. 2014; Liu et al. 2022). Meng studied the deformation and volume expansion characteristics of rock under different confining pressures as well as cyclic loading and unloading conditions, and the author found that confining pressure significantly affects the expansion of micro-cracks in rock (Meng et al. 2021). Xiao studied the internal fracture and deformation law of rock under cyclic loading and unloading conditions using the rock residual strain method (Xiao et al. 2009). Momeni analyzed the stress–strain relationship of granite under different peak loads and different loading amplitudes under cyclic loading. The author observed three-stage characteristics in the deformation and failure process of the rock (Momeni et al. 2015). Ghamgosar studied the development behavior of cracks in rock under split conditions and found that more high-density cracks can be generated in the rock under the action of cyclic loading (Ghamgosar and Erarslan 2016). Xu propose a two-dimensional meso-scale elasto-plastic damage model to simulate the fatigue behavior of rock during cyclic loading. And they found that specimen failure occurs once the cumulative axial deformation under cyclic loading reaches the axial strain that corresponds to the post-peak phase of the stress–strain curve from a conventional uniaxial test simula-tion (Xu et al. 2021). Zhou used computed tomography imaging technology and digital image correlation testing technology to analyze the expansion and evolution behavior of the internal fractures of a rock under the action of load. The author found that the deformation and failure process of the rock can be divided into five stages (Zhou et al. 2020). Zhang used acoustic emission and high-speed camera technology to find that the initiation and development of rock cracks were consistent with the change in the crack inclination angle (Zhang et al. 2019). Ma systematically examine the influence of fissure angles, lengths and positions on the mechanical response. They found that the burst tendency of composite specimens can be effectively reduced by a fissure in either the weak (lower, SWL) or strong (upper, SWU) layers. However, the effectiveness of the fissure depends on the angle and length of the fissure (Ma et al. 2023). Ghasemi used the fluorescent labeling method to identify the initial fractures in a rock as intergranular and transgranular fractures, and the intergranular fracture expansion under load eventually caused rock failure (Ghasemi et al. 2020). Liu used the dual-damage constitutive model to study the damage and fracture mechanism of deep surrounding rock under high confining pressure (Liu et al. 2020). For the research on the characteristics of acoustic emission signals during the rock cyclic loading and unloading stages, Ji found that under the triaxial test of cyclic loading and unloading of monzonitic granite rock samples, the distribution of the characteristic parameters of granite rock acoustic emission signals has four variation characteristics stages (Ji et al. 2012). Harte analyzed the relationship between information entropy and system stability and used information entropy for earthquake prediction (Harte and Vere-Jones 2005). Bressan claimed that the evolution of entropy value is a description of the expansion process of a thermodynamic system from disorder to order (Bressan et al. 2017). Ju used industrial CT to achieve fracture description in rock deformation and failure engineering under different confining pressures (Ju et al. 2018). Cai found that the development process of rock fracture can be divided into three stages: fracture initiation, fracture expansion and fracture acceleration (Cai et al. 2004). Duan studied the fracture volume variation of Longmaxi rock under uniaxial compression and found that the fracture volume variation presents two stages of steady growth and accelerated jet expansion (Duan et al. 2019). Chen's study found that the higher the confining pressure, the larger the volume strain of fracture during rock fracture (Chen et al. 2020). Zhang found that the variation of acoustic emission energy rate could reflect the development of rock cracks, and when the AE energy rate was high, the cracks inside the sample gradually increased (Zhang et al. 2006). Akdag found that the AE energy rate corresponded well with the stage of rock compression and deformation, and the AE energy rate could reflect the fracture characteristics during rock deformation and failure (Akdag et al. 2021).

At present, most of the research results focus on the change law of single parameter of rock specimen surface and fracture volume strain, lacking of considering the influence of rock volume change on fracture volume development. In fact, deep stratum rocks are affected by excavation, rock stress field is in the process of dynamic adjustment, and rock volume change has a direct control effect on fracture expansion.

Therefore, the research focus of this paper is to investigate the fracture volume evolution law of monzogranite under cyclic loading, and then to describe the development of cracks in rock samples by introducing fracture volume strain ratio coefficient. This parameter can reflect the variation of fracture volume strain and its relation to the total volume strain of rock samples during loading failure process, and establish AE entropy to describe the development characteristics of fractures in rock. The results are of great significance for rock engineering stability analysis.

2 Experimental scheme design

The experimental rock samples were obtained from the monzogranite at a depth of -1600 m in a gold mine located in the Jiaodong Peninsula. The core column of the collected samples was processed into a cylindrical standard specimen of 50 mm and 100 mm, the allowable deviation of the height and diameter of the rock samples was 3 mm, and the allowable deviation of the unevenness of the two ends of the rock samples was 0.05 mm.

The experimental mechanical loading system was a TAW-2000 rock triaxial rigidity testing machine and a PCI-2 acoustic emission instrument located at the University of Science and Technology Beijing. The experimental instruments are shown as Fig. 1. The resonant frequencies of the acoustic emission probes used for the experiment were 55 kHz and 150 kHz, respectively. Both the pre-gain and the threshold value were 40 dB.

Fig. 1
figure 1

Schematic of experimental apparatus

Loading scheme design: At the beginning of the experiment, the testing machine was preloaded by axial pressure at 500 N/s until 5 MPa. It should be noted that the unloading stress is set to 5 MPa to ensure that the indenter of the testing machine is in a close contact with the rock sample. After stabilizing the data, they were cleared and changed to deformation control data, and the deformation rate was 0.01 mm/min. Rocks from a gold mine at a depth of − 1600 m are characterized as high-strength and high-hardness. Consequently, the rock samples obtained from the mine had high strength and hardness, and the uniaxial compressive strength reached 220 MPa. The first axial loading was at 100 MPa. The initial load is set to 100 MPa to prevent fatigue damage due to excessive cyclic loading under low stress, while also avoiding significant acoustic emission. To ensure that the experiment can be conducted between seven and nine cycles before the rock sample failure, the loading rate was kept constant, and the stress increase of cyclic loading was 20 MPa before the rock sample failure. Three groups of cyclic loading and unloading experiments were set under different confining pressures (10 MPa, 20 MPa, and 30 MPa). When the target loading stress of each cycle was reached, the samples were unloaded until the stress decrease to 5 MPa, then the next cycle of loading was performed till the specimen failed. The corresponding stress–strain curves and fracture volume strain curves are shown in Figs. 2 and 3.

Fig. 2
figure 2

The monzogranite’s stress–strain curves under cyclic loading and unloading with different confining pressures

Fig. 3
figure 3

Volume strain curves of monzogranite fractures under different confining pressures

3 Deformation characteristics of monzogranite

3.1 Evolution of fissure volume strain ratio coefficient

The initiation and evolution of fractures in the rock changed the rock physical and mechanical properties, and the volume strain of the fracture is an important parameter for describing the fracture evolution during cyclic loading and unloading. Total rock volume strain is equal to the elastic volume strain plus fracture volume strain (Erarslan 2016; Liu and Dai 2021), as follows:

$$\varepsilon_{v} = \varepsilon_{ve} + \varepsilon_{vc} ,$$
(1)

where \(\varepsilon_{v}\) is the total volume strain, \(\varepsilon_{ve}\) is the elastic volume strain, and \(\varepsilon_{vc}\) is the crack volume strain.

The total volume strain of the rock sample can be determined using Eq. (2):

$$\varepsilon_{v} = \varepsilon_{1} + 2\varepsilon_{3} ,$$
(2)

where \(\varepsilon_{1}\) is the axial strain, and \(\varepsilon_{3}\) is the hoop strain.

The volume strain of the rock sample elastic body can be determined using Eq. (3):

$$\varepsilon_{ve} = \frac{1 - 2\upsilon }{E}(\sigma_{1} + 2\sigma_{3} ),$$
(3)

where \(\upsilon\) is the Poisson’s ratio in the elastic stage, E is the elastic modulus,\(\sigma_{1}\) is the maximum principal stress, and \(\sigma_{3}\) is the minimum principal stress.

The proportion of the volume strain of the rock sample fractures can be calculated using Eq. (4):

$$F = \frac{{\varepsilon_{vc} }}{{\varepsilon_{v} }},$$
(4)

where \(F\) is the proportion coefficient of the fracture volume strain, which is the proportion of the fracture volume strain in the rock to the total volume strain and describes the relative scale of the fracture volume strain in the rock to the total volume strain.

The variation characteristics of the volume strain ratio of the monzogranite fissures under different confining pressures were obtained according to Eq. (4). The lower the confining pressure is, the weaker the lateral deformation constraint of the rock is. To better describe the fracture development behavior inside the rock, the volume strain curve of the rock sample fractures under 10-MPa confining pressure was used for analysis and investigation. After the initial cyclic unloading stage, the fracture volume strain curve of the rock sample had obvious residual fracture volume strain because the original fractures of the rock were compacted and closed during the initial cyclic loading stage. The fractures in the rock were completely compacted and closed, and the change rule of the fracture volume strain ratio curve in the later period was roughly the same. Therefore, the representative fracture volume strain ratio curve under the first, second, and fourth cycles was analyzed, as shown in Figs. 3 and 4.

Fig. 4
figure 4

Fracture volume strain ratio coefficient under the first, second, and fourth cyclic loads at a confining pressure of 10 MPa: a First cycle; b Second cycle; c Fourth cycle

The ratio of the fracture volume strain to the total volume strain was determined by the ratio of the fractured volume strain to the total volume strain, which can show the relative scale of the fractured volume strain to the total volume strain under different stress states. Relative stress is defined as the ratio of stress to peak stress at different states. Figures 4 and 5 show that under different stress states, fracture volume strain ratio coefficient curve of the monzogranite presented a logarithmic change. The unloading process has three stages: linear growth, slow decrease, and sharp decline.

Fig. 5
figure 5

Fracture volume strain ratio coefficient of monzogranite under different confining pressures: a Confining pressure: 10 MPa; b Confining pressure: 20 MPa; c Confining pressure: 30 MPa

In the cyclic loading stage of the rock sample, when the relative stress was between 0 and 0.05, the stage was the crack closure stage. The cracks and pores inside the rock were compressed and closed under the action of pressure, and the volume strain ratio of the crack was from 0.6 to 0.7, showing that the rock deformation in this stage was mainly fractured compaction deformation. As the stress level increased, the elastic deformation gradually increased, and the growth rate of the fracture volume strain ratio gradually decreased, stabilizing between 0.8 and 0.85. This behavior indicates that the fracture volumetric and elastic volume strains accounted for the relative volume strain during the loading process. The scale was stable and unchanged, but the fracture volume strain accounted for the main part. When the relative stress level was greater than 0.4, cracks were gradually formed in the rock. Correspondingly, the volume strain ratio of the cracks increased slowly, indicating that the cracks in the rock continued to develop. Meanwhile, the opening degree of the cracks in the rock sample increased under the load, and the volumetric and elastic strains of the cracks increased simultaneously. This stage is the fissure stable development stage.

During the cyclic unloading stage of the rock sample, the volume strain ratio of the fracture increased linearly as the stress level decreased, showing that the elastic deformation inside the rock recovered at first. Then, the volume strain of some uncompacted fractures recovered, as shown in Fig. 4a. When the stress level was unloaded to less than the crack initiation stress, the cracks were gradually closed, the volume strain ratio of the cracks decreased nonlinearly, and the relative scale of the total rock volume strain gradually decreased. Finally, the volume strain ratio of the cracks was 0.63, indicating that the sample contained volume strains of compacted fractures. This phenomenon showed that the internal fractures in the rock were compacted and closed after one cyclic loading, and the volume of the closed fractures accounted for 63% of the total volumetric deformation under the stress state.

As the number of cycles increased, the proportion coefficient of the fracture volume strain increased gradually under the same stress level, indicating that the relative scale of the fracture volume strain in the rock to the total volume strain continuously expanded. The fractures in the rock were constantly sprouting and developing, and the damage in the rock was constantly increasing. When the stress level reached the peak stress, the proportion of the fracture volume strain was close to 1, indicating that the deformation at this time was mainly caused by the fracture strain, and the rock was damaged.

3.2 Fracture morphology analysis

To characterize the fracture development inside granite rocks under different confining pressures, the failure map of monzogranite samples under different confining pressures and the scanning electron microscope images of the fracture surface of the rock samples are presented in Figs. 6 and 7.

Fig. 6
figure 6

Failure patterns of rocks under different confining pressures: a Confining pressure: 10 MPa; b Confining pressure: 20 MPa; c Confining pressure: 30 MPa

Fig. 7
figure 7

Cross-section morphology of rocks under different confining pressures: a Confining pressure: 10 MPa; b Confining pressure: 20 MPa; c Confining pressure: 30 MPa

Figures 6 and 7 show that the monzogranite samples mainly had two failure modes: shearing and tensioning. Under the confining pressure of 10 MPa, the rock sample was mainly damaged by tension: longitudinal tearing cracks were initiated, and the bottom of the rock sample peeled off. At this time, the tensile strength of the rock sample is higher, so it fails in tension. The micro-cracks were mainly concentrated on both sides of the main crack, and many cluster fractures were inside the rock sample. The inclusion distribution was uneven, and the number of ruptured inclusions increased considerably along the main crack direction. The failure mode of the rock sample under the confining pressure of 20 MPa was the combined tension and shear failure modes. A bulging crack occurred in the middle of the rock sample. The failure crack of the rock sample was mainly shear failure, and the number of secondary cracks was considerably reduced. Under the confining pressure of 30 MPa, the rock sample suffered shear damage where the internal intergranular cracks and pores developed in the rock. With the increase of confining pressure, the tensile strength of the rock sample decreases, while the shear strength increases. When the shear stress component exceeds the tensile stress component, the rock sample fails in shear. The volume strain ratio of the fractures at the upper limit stress of the rock sample under the 30 MPa confining pressure steadily increased, which resulted from the random development of new micro-cracks and pores inside the rock. The scale of developed fractures and pores did not increase until a certain range because of the high-confining-pressure constraint. The increase in cyclic stress generated sufficient new pores and cracks in the rock sample, making the crack volume grow steadily while an obvious directionality appeared along the shear failure surface. Under the low confining pressure, the confinement effect of the confining pressure was not obvious, the scale of randomly developed fractures and pores in the rock continued to develop, and the continuous friction between the fracture surfaces produced numerous fractured inclusions and micro-cracks; the volume strain ratio of the fractures fluctuated and did not substantially increase, as shown in Fig. 8.

Fig. 8
figure 8

Proportion coefficient of volume strain of monzogranite fractures under different confining pressures

3.3 Characteristics of AE entropy

Entropy is an indicator used for describing system complexity and disorder. The acoustic emission entropy value can be used to describe the distribution of acoustic emission amplitude and show crack development (Harte and Vere-Jones 2005; Bressan et al. 2017) since the process of crack development and expansion in rock mass is a nonlinear process from disorder to order. When fractures in rocks develop into fractures of the same scale, the acoustic emission amplitude distribution is typically relatively close, and the entropy value is large; when the internal fractures in rocks develop into fractures of different scales, the acoustic emission amplitude distribution is wider, and the entropy value is low. Currently, few studies are using the acoustic emission entropy theory to study crack evolution behavior during rock loading. According to the information entropy calculation method in probability statistics, it is assumed that all possible distributions of results satisfy the following probability conditions:

$$0 \le \delta_{i} \le 1\;\;\;\left( {i = 1,2, \ldots ,n} \right),$$
(6)
$$\sum\limits_{i = 1}^{n} {\delta_{i} = 1} ,$$
(7)
$$H\left( x \right) = H\left( {\delta_{1} ,\delta_{2} , \ldots ,\delta_{n} } \right) = - \sum\limits_{i = 1}^{n} {\delta_{i} \log \delta_{i} } .$$
(8)

Normalizing the entropy value,

$$H\left( t \right) = H\left( x \right)/\log N,$$
(9)

where \(H\left( x \right)\) represents the entropy value; \(N\) represents all grade states in the whole system, and every 5 dB is a level in the acoustic emission amplitude classification; \(\delta_{i}\) is the occurrence probability of this state, which is the ratio of the number of acoustic emission events to the total number of events within the same amplitude range; \(H\left( t \right)\) is the normalized time entropy.

The entropy value of the acoustic emission event was calculated as follows:

  1. (a)

    A Python program was written to divide the time calculation window. Every 1000 s was selected as the calculation window to count the number of acoustic emission events because the cycle elapsed for a long time.

  2. (b)

    The acoustic emission amplitude levels were classified in the time period, with a level of every 5 dB, and the distribution of the acoustic emission signals was counted in each level.

  3. (c)

    The probability of acoustic emission events was calculated in each window.

  4. (d)

    Equations (8) and (9) were used to calculate the entropy value of acoustic emission in the rock cycle loading and unloading stages.

Figure 9 illustrates the stress in the whole process of the cyclic loading and unloading deformation. The failure stress of monzogranite is captured under different confining pressures. The curve of acoustic emission cumulative energy and acoustic emission energy rate with time are correspondingly recorded.

Fig. 9
figure 9

Variation curves of acoustic emission cumulative energy and acoustic emission energy rate of rocks under different confining pressures: a Confining pressure: 10 MPa; b Confining pressure: 20 MPa; c Confining pressure: 30 MPa

According to the calculation method of the acoustic emission b value, the corresponding Python code was compiled to obtain the variation characteristics of the acoustic emission b value during the whole process of the monzogranite cyclic loading and unloading (Erarslan 2016; Liu and Dai 2021; Martin and Chandler 1994; Yu et al. 2021; Zhang and Zhang 2017; Liu et al. 2021; Wang et al. 2023; Yang et al. 2022; He 2022). The acoustic emission b value and entropy value curves of the whole deformation process and failure of the monzogranite under different confining pressure cyclic loads were obtained, as illustrated in Fig. 10.

Fig. 10
figure 10

Curves of b and entropy values of acoustic emission of rocks under different confining pressures: a Confining pressure: 10 MPa; b Confining pressure: 20 MPa; c Confining pressure: 30 MPa

Figures 9 and 10 show that the overall acoustic emission b value curve of the monzogranite failure process exhibited a downward trend while the acoustic emission entropy exhibited a fluctuating upward trend. The change trend of the acoustic emission b value and that of the acoustic emission entropy value had opposite trends. Compared with the rock samples under 30 MPa confining pressure, the rock samples under the confining pressure of 10 MPa and 20 MPa in the initial stage showed a lower the acoustic emission b value fluctuated, indicating the acoustic emission behavior inside the rock at this stage. The signal was weak, corresponding to the acoustic emission cumulative energy as illustrated in Fig. 9. The acoustic emission cumulative energy was at a low level, and the rock was in the stage of fissure closure, with mainly small fractures of the same scale. The acoustic emission signal amplitude was relatively close, and the entropy value was relatively high; As the cyclic upper limit stress increased, the rock entered the plastic stage, the acoustic emission b value began to decrease slowly, larger-scale fractures were gradually generated inside the rock sample, and the acoustic emission signal energy increased. At the internal accumulation stage, the acoustic emission entropy value had a short decline process, the amplitude signal was widely distributed, and the acoustic emission energy rate increased considerably. When the upper limit stress of the rock internal accumulation was close to the peak stress, the acoustic emission was in the quiet period just before the rock rupture. The acoustic emission entropy value decreased sharply, and the entropy value was at the lowest level. When the stress level was equal to the peak stress, the acoustic emission signal had a large, sudden increase in the acoustic emission cumulative energy. The corresponding acoustic emission entropy value rapidly increased to the maximum value. The crack was in an unstable development stage and eventually turned into an unstable failure. The acoustic emission entropy value before and after the rock rupture showed an obvious “V”-shaped change characteristic. The lowest point of the entropy value corresponded to the quiet period of the acoustic emission signal. When the confining pressure was 30 MPa, 20 MPa and 10 MPa, the lowest entropy value was ahead of the rock rupture point at 813.5 s 625.3 s and 343.5 s, respectively. Therefore, the higher the confining pressure was, the higher the acoustic emission entropy was. The timing point of the lowest value was earlier than the rock rupture point.

4 Discussion

Figure 11 shows the schematic diagram of the curve change of the monzogranite time-acoustic emission event number-entropy value-fracture volume strain ratio under different confining pressures. As the stress level increased, the number of the acoustic emission events and fracture volume strain ratio curve of the monzogranite showed an upward trend. This phenomenon showed that under the action of cyclic loading, the relative scale of the fracture volume strain to the total volume strain increased, and the fractures in the rock were constantly developing. The acoustic emission entropy value had an obvious stage division. The acoustic emission entropy value showed a rapid decrease in the crack closure stage, a slow rise in the crack stable development stage, an unstable crack development stage and a post-peak residual deformation stage, and a rapid upward trend. Overall, the acoustic emission entropy value showed that the more fractures developed in the rock, the larger the acoustic emission entropy value was. Therefore, the acoustic emission entropy value can better reflect the development of the internal fractures in the rock.

Fig. 11
figure 11

Time acoustic emission event number-entropy-volume strain ratio curve of different confining pressures of monzogranite: a Confining pressure: 10 MPa; b Confining pressure: 20 MPa; c Confining pressure: 30 MPa

The initial loading stage of the rock was the crack closure stage. At this stage, the rock was compressed and closed. The volume strain of the rock was mainly the volume strain of the crack. The acoustic emission entropy value was at a high level. As the stress level increased, the rock entered the elastic deformation stage, the number of acoustic emissions signal events was not obvious, and the entropy value was at a low level. As the stress level continued to increase, the rock entered the crack stable development stage, and the acoustic emission signal and amplitude were generated. The distribution range gradually expanded, and the acoustic emission entropy began to increase slowly. When the stress level was close to the peak strength, the volume strain ratio of the fracture increased rapidly, indicating that the fracture scale increased rapidly at this stage, and the stage was the unstable fracture expansion stage. The acoustic emission signal, amplitude, and entropy value increased gradually. When the stress reached the rock peak intensity, many acoustic emission signal events were generated; the large-scale signal accounted for the main part of the acoustic emission signal, and the acoustic emission entropy value was at a relatively high level. After reaching the peak intensity, the rock entered the post-peak residual deformation stage. A macroscopic fracture surface was formed inside the rock, and the fracture volume strain ratio was close to 1, indicating that the rock fracture volume accounted for the main part of the main total volume strain at this stage.

The acoustic emission entropy value before and after the rock rupture presented a distinct “V”-shaped change, and the lowest point of the “V”-shaped change of the acoustic emission entropy value corresponded to the quiet period of the acoustic emission. Moreover, the higher the confining pressure of the rock was, the greater the mutation degree of the acoustic emission entropy value was, the greater the entropy reduction degree was, and the earlier the rock rupture point was at the lowest point of the acoustic emission entropy value.

5 Conclusion

The fissure volume evolution and acoustic emission entropy characteristics of monzogranite under cyclic loading were obtained through triaxial cyclic loading–unloading tests and acoustic emission monitoring analysis for monzogranite taken from a gold mine located in the Jiaodong Peninsula at a depth of − 1600 m. The main conclusions of this study are as follows:

  1. (1)

    The fracture volume strain ratio factor is introduced to comprehensively describe the relative relationship between fracture volume strain and total volume strain. The experimental results demonstrate that the fracture volume strain ratio factor can effectively reflect the relative scale of fracture development in the rock. As the stress level increased, the fracture development scale gradually expanded, and the fracture volume strain ratio factor showed an upward trend.

  2. (2)

    Under different confining pressures, the volume strain ratios of fractures had different variation characteristics which were obviously determined by confining pressure. As the confining pressure increased, the volume strain ratio curve of the fracture showed three changing characteristics: stable fluctuation, slow rise, and finally rapid rise.

  3. (3)

    The acoustic emission entropy value effectively reflected the evolution behavior of crack propagation inside the loaded rock. As the number of the acoustic emission events increased, the acoustic emission entropy value showed a dominant upward trend. The acoustic emission entropy value before and after the rock failure presented a “V”-shaped change. The lowest point of the “V”-shaped change in the acoustic emission entropy value corresponded to the quiet period of the acoustic emission. A higher confining pressure results in a greater the mutation degree of the acoustic emission entropy value. Also, the higher the confining pressure was, the greater the decrease in the entropy value was, and the earlier the rock rupture point was reached.