Synergetic Theory of Information Entropy Based on Failure Approach Index for Stability Analysis of Surrounding Rock System

It is generally acknowledged that the stability evaluation of surrounding rock denotes nonlinear complex system engineering. In order to accurately and quantitatively assess the safety states of surrounding rock and provide a scientific basis for the prevention and control of surrounding rock stability, the analysis method of the synergetic theory of information entropy using the failure approach index has been proposed. By means of deriving the general relationship between the total two-dimensional plastic shear strain and the total three-dimensional plastic shear strain and obtaining the numerical limit analysis step of the plastic shear strain, the threshold value of the ultimate plastic shear strain can be determined, which has provided the key criterion for the calculation of the information entropy based on the failure approach index. In addition, combining with the synergetic theory of the principle of maximum information entropy, the evolution equation of the excavation step and information entropy based on the failure approach index of the surrounding rock system in underground mining space are established, and the equations of the general solution and particular solution as well as the expression of the destabilizing excavation step are given. To account for this, the method is applied to analyze the failure states of the floor surrounding rock after the mining of the 71 coal seam in Xutuan Coal Mine and involve the disturbance effect and stability control method of the underlying 72 coal seam roof from the macroscopic and microscopic aspects. Consequently, the validity of the analysis method of synergetic theory of information entropy based on the failure approach index has been verified, which presents an updated approach for the stability evaluation of surrounding rock systems that is of satisfactory capability and value in engineering applications.


Introduction
In recent years, subway tunnels, hydraulic tunnels, urban underground space, coal underground mining, and other projects have increased, and the stability of geotechnical bodies in underground engineering needs to be solved urgently.Therefore, the evaluation of the stability of the surrounding rock in underground engineering has been a research hotspot in the discipline of rock mechanics and engineering and a key technical difficulty that needs to be solved urgently in practical engineering [1][2][3][4].
In terms of the stability analysis of rock and soil mass, many scholars use various information entropy models to analyze and evaluate the stability of surrounding rock.Xu et al. [5] applied nonlinear scientific theories, such as the information entropy of plastic deformation, the mutation theory, etc., to investigate the nonlinear stability problems of rock mass so that the instability process and failure mechanism of rock mass could be Entropy 2023, 25 revealed and a failure analysis method and instability criterion of rock mass that conform to reality could be deduced.According to the information entropy theory, Liu et al. [6] proposed a method for extracting the vibration characteristics of coal gangue based on the information entropy of the Hilbert spectrum and attained the outcome that the information entropy of the Hilbert spectrum when top coal falls performs better than that when mixing the coal gangue.In the light of the uncertainty measure theory, the information entropy theory, the reliability identification criteria, and the ranking criteria, He et al. [7] established a classification and ranking model of surrounding rock stability, providing a novel clue for the grading evaluation of tunnel surrounding rock stability.Li et al. [8] selected the formation lithology, the groundwater level, the rock occurrence, and the surrounding rock grade as risk evaluation indexes of water and mud inrush with respect to the geological and engineering factors of the tunnel, determined the index weights of each factor by virtue of the information entropy theory, and established an unascertained measurement model for risk evaluation of water and mud inrush.Zhao et al. [9] adopted the information entropy theory to determine the weight of each index and created a multi-index comprehensive measure evaluation model.Yang et al. [10] employed the stress information entropy principle to study the evolution law of mining-induced overburden stress, creating conditions for the establishment of the relationship between the change of stress information entropy and the evolution of mining-induced fractures, thereby laying a foundation for the formation of water channels and the identification of water hazard sources.In order to explore the damage evolution process of granite under triaxial compression, Ji et al. [11] introduced the AE amplitude information entropy variation coefficient and correlation dimension to identify granite failure precursors, presenting a fresh idea for rock instability prediction.With regard to the gray correlation analysis method and the entropy weight method, Zhang et al. [12] established a tunnel deformation analysis model for canopy settlement and lateral wall horizontal shrinkage during tunnel excavation and drew the conclusion that the main reason for the collapse accident of the Hilbert spectrum could be attributed to the tunnel overdigging.Moreover, Gao et al. [13] proposed a thermal infrared analysis method of exponential infrared image entropy to study the infrared radiation characteristics and fracture precursors during the process of rock fracture.Zhang et al. [14] put forward the concept of failure approach index on the basis of yield approach index and quantitatively revealed the degree to which rock and soil mass approaches or reaches the yield or failure state in the most unfavorable way under loads.Hence, the calculation formula can be expressed as: where FAI represents the failure approach index; ω refers to the complementary parameter of yield approach index (YAI), ω = 1 − YAI; and FD signifies the failure degree: where γ p denotes the plastic shear strain; γ p = e When 0 ≤ FAI < 1.0, the rock and soil mass is in the elastic stress state, and the greater the FAI value in this interval, the closer to the yield state it will become; when 1.0 ≤ FAI < 2.0, it is in the plastic yield state; at the time of FAI ≥ 2.0, it is in the failure state.
In simpler terms, the stability states of surrounding rock in engineering directly affect the stability control of surrounding rock.In the meantime, the plastic shear strain values and the ultimate plastic shear strain value of the surrounding rock are not prone to being quantitatively expressed through numerical calculation, and the surrounding rock has complicated efforts to determine the stability states.Accordingly, the analysis method of the synergetic theory of information entropy based on the failure approach index is introduced in this study to analyze this kind of problem.Under the condition of the principal strain order ε 1 ≤ ε 2 ≤ ε 3 , the general relationship between the total twodimensional plastic shear strain and the total three-dimensional plastic shear strain has been derived, and the numerical limit analysis method of plastic shear strain has been proposed so that the correctness and reliability of the method can be adequately validated.Therefore, the threshold of the ultimate plastic shear strain will be determined, which provides the essential criterion for the calculation of information entropy based on the failure approach index.In relation to the synergetic theory of the maximum information entropy principle, the evolution equation of the information entropy function based on the failure approach index and the excavation step in the surrounding rock system of underground mining space has been established.In addition, the general and solution formulae and the expression of the unstable excavation step are obtained, which provides a theoretical basis for the stability analysis of the surrounding rock system.Ultimately, the stability status of the 7 1 coal seam floor in Xutuan Coal Mine (Figure 1) has been analyzed from the macroscopic and microscopic aspects.The research results will extend the failure approach index theory and present a method and technical reference for the stability evaluation of similar rock and soil masses.to being quantitatively expressed through numerical calculation, and the surrounding rock has complicated efforts to determine the stability states.Accordingly, the analysis method of the synergetic theory of information entropy based on the failure approach index is introduced in this study to analyze this kind of problem.Under the condition of the principal strain order ε 1 ≤ ε 2 ≤ ε 3 , the general relationship between the total two-dimensional plastic shear strain and the total three-dimensional plastic shear strain has been derived, and the numerical limit analysis method of plastic shear strain has been proposed so that the correctness and reliability of the method can be adequately validated.Therefore, the threshold of the ultimate plastic shear strain will be determined, which provides the essential criterion for the calculation of information entropy based on the failure approach index.In relation to the synergetic theory of the maximum information entropy principle, the evolution equation of the information entropy function based on the failure approach index and the excavation step in the surrounding rock system of underground mining space has been established.In addition, the general and solution formulae and the expression of the unstable excavation step are obtained, which provides a theoretical basis for the stability analysis of the surrounding rock system.Ultimately, the stability status of the 71 coal seam floor in Xutuan Coal Mine (Figure 1) has been analyzed from the macroscopic and microscopic aspects.The research results will extend the failure approach index theory and present a method and technical reference for the stability evaluation of similar rock and soil masses.

Relationship between Two-Dimensional Plastic Shear Strain and Three-Dimensional Plastic Shear Strain
In the numerical simulation software, the plastic shear strain adopts the three-dimensional plastic shear strain increment form.Nevertheless, the more concise two-dimensional plastic shear strain full form is often employed in the derivation and calculation of the theoretical model (for example, the full form   is used in the expression of the information entropy based on the failure approach index) [15].In order to facilitate calculation, it is necessary to establish the relationship between the full form   and the incremental form ∆ε ps .
Under the condition of the principal strain order ε 1 ≤ ε 2 ≤ ε 3 , the expression of the full form   can be represented as:

Relationship between Two-Dimensional Plastic Shear Strain and Three-Dimensional Plastic Shear Strain
In the numerical simulation software, the plastic shear strain adopts the three-dimensional plastic shear strain increment form.Nevertheless, the more concise two-dimensional plastic shear strain full form is often employed in the derivation and calculation of the theoretical model (for example, the full form γ p is used in the expression of the information entropy based on the failure approach index) [15].In order to facilitate calculation, it is necessary to establish the relationship between the full form γ p and the incremental form ∆ε ps .
Under the condition of the principal strain order ε 1 ≤ ε 2 ≤ ε 3 , the expression of the full form γ p can be represented as: Entropy 2023, 25, 1237 where ε ps 3 and ε ps 1 signify the maximum and the minimum plastic principal strain, respectively.
For the purpose of facilitating the subsequent derivation, γ p is expressed in the incremental form: In the numerical simulation software, the strain-softened plasticity model adopts the incremental form ∆ε ps , whose expression can be determined as: where ∆ε Under the condition of the principal stress sequence σ 1 ≤ σ 2 ≤ σ 3 , the expression of the shear plastic potential function g s [16] is obtained as: where where ψ indicates the dilation angle.
With respect to the shear failure, the form of the non-associated flow rule is denoted as follows: where λ s depicts the scaling factor.According to Equation (9), the three plastic principal strain increment expressions can be attained as: In relation to Equation (10), the relationship between the principal strain increments ∆ε ps 3 and ∆ε ps 1 is established, and the expression can be signified as: In accordance with Equations ( 6) and (11), Equation ( 5) can be rewritten as: Based on Equation (11), Equation ( 12) can be denoted as: Combined with Equations ( 4), ( 12) and ( 13), the incremental form ∆ε ps can be expressed as: Entropy 2023, 25, 1237 With regard to Equation ( 14), ∆ε ps and ∆γ p are affected by the dilation angle of rock and soil material ψ, and 0 • ≤ ψ < 90 • .The ratio range of both can be determined as follows: Under the assumption of small deformation, no significant difference between the incremental form and the full form of plastic shear strain can be examined without considering the rotation of the stress spindle [17,18].Therefore, the relationship between the full amount of the two-dimensional plastic shear strain γ p and the three-dimensional plastic shear strain ε ps can be written as: Under the limit strain state, Equation ( 16) can be represented as where γ p f and ε ps f depict the two-dimensional and the three-dimensional ultimate plastic shear strain, respectively.

Calculation Process of the Ultimate Plastic Shear Strain
With the development of numerical methods for rock and soil mechanics, the numerical limit analysis method has gradually emerged, which has both comprehensive applicability and feasible practicability [19][20][21].In this research, the ultimate load of rock and soil mass is determined by the overload method under the strain-softening constitutive model, which presents the prerequisite to defining the ultimate plastic shear strain of rock and soil mass, as displayed in Figure 2.With regard to Equation ( 14), ∆ε ps and ∆  are affected by the dilation angle of rock and soil material , and 0°≤  < 90°.The ratio range of both can be determined as follows: Under the assumption of small deformation, no significant difference between the incremental form and the full form of plastic shear strain can be examined without considering the rotation of the stress spindle [17,18].Therefore, the relationship between the full amount of the two-dimensional plastic shear strain   and the three-dimensional plastic shear strain ε ps can be written as: Under the limit strain state, Equation ( 16) can be represented as where    and    depict the two-dimensional and the three-dimensional ultimate plastic shear strain, respectively.

Calculation Process of the Ultimate Plastic Shear Strain
With the development of numerical methods for rock and soil mechanics, the numerical limit analysis method has gradually emerged, which has both comprehensive applicability and feasible practicability [19][20][21].In this research, the ultimate load of rock and soil mass is determined by the overload method under the strain-softening constitutive model, which presents the prerequisite to defining the ultimate plastic shear strain of rock and soil mass, as displayed in Figure 2.This study predominantly focuses on the ultimate plastic shear strain under the ideal plastic model and employs the numerical limit analysis method to determine the ultimate plastic shear strain of rock and soil materials.In the first place, with respect to the specific engineering background, the basic procedures, such as creating a geometric model, defining the strain softening constitutive model and parameters, setting boundary conditions, applying loads, etc., are established.In addition, the monitoring points for the plastic shear strain are set in the calculation model.Afterwards, the extraction and output programs of the plastic shear strain are self-compiled by the Python language and embedded in the numerical simulation program to output the plastic shear strain values of the monitoring points.Ultimately, the 'axial load-plastic shear strain' curve is drawn to compare and analyze the variation trend of the plastic shear strain at each monitoring point under each level of load and determine the ultimate plastic shear strain value of rock and soil materials.
The key points of structural deformation weakness for the ultimate strain analysis of geotechnical materials have been normally selected.For instance, through the ultimate strain analysis of concrete blocks [22], 12 representative key nodes (elements) were chosen, and the load of each stage as well as the corresponding plastic shear strain values This study predominantly focuses on the ultimate plastic shear strain under the ideal plastic model and employs the numerical limit analysis method to determine the ultimate plastic shear strain of rock and soil materials.In the first place, with respect to the specific engineering background, the basic procedures, such as creating a geometric model, defining the strain softening constitutive model and parameters, setting boundary conditions, applying loads, etc., are established.In addition, the monitoring points for the plastic shear strain are set in the calculation model.Afterwards, the extraction and output programs of the plastic shear strain are self-compiled by the Python language and embedded in the numerical simulation program to output the plastic shear strain values of the monitoring points.Ultimately, the 'axial load-plastic shear strain' curve is drawn to compare and analyze the variation trend of the plastic shear strain at each monitoring point under each level of load and determine the ultimate plastic shear strain value of rock and soil materials.
The key points of structural deformation weakness for the ultimate strain analysis of geotechnical materials have been normally selected.For instance, through the ultimate strain analysis of concrete blocks [22], 12 representative key nodes (elements) were chosen, and the load of each stage as well as the corresponding plastic shear strain values were extracted and output by means of the self-written Python program.Meanwhile, by comparing and analyzing the variation trends of the plastic shear strain of 12 key elements with the increase in axial load, the threshold value of the ultimate plastic shear strain of concrete material was finally determined.Therefore, the ultimate plastic shear strain in geotechnical engineering can also be obtained by the numerical limit analysis method, and the specific flow chart is depicted in Figure 3. were extracted and output by means of the self-written Python program.Meanwhile, by comparing and analyzing the variation trends of the plastic shear strain of 12 key elements with the increase in axial load, the threshold value of the ultimate plastic shear strain of concrete material was finally determined.Therefore, the ultimate plastic shear strain in geotechnical engineering can also be obtained by the numerical limit analysis method, and the specific flow chart is depicted in Figure 3.

Solution and Verification of Ultimate Plastic Shear Strain
Given that there are numerous studies on the ultimate strain of concrete [23,24], the numerical limit analysis method is used to calculate and analyze the plastic shear strain of a C25 ordinary concrete cube specimen in this study.Hence, the ultimate plastic shear strain value of C25 ordinary concrete is obtained, and a comparison is made with the results of previous studies to verify the reliability of the numerical limit analysis method of plastic shear strain of C25 ordinary concrete.
In this research, the finite difference software FLAC3D of version 6.0 is adopted to numerically simulate the C25 concrete cube specimen.First, the ultimate load of C25 concrete is determined by the overload method.The calculation model is a cubic specimen with a side length of 150 mm, and each side is divided into 20 grids, as shown in Figure 4a.Afterwards, the strain-softened plasticity model is employed, whose material attribute parameters are defined.Hence, the physical and mechanical parameters of C25 concrete are exhibited in Table 1.Furthermore, the bottom surface of the calculation model is fixed, while the top surface is subjected to a downward vertical plane load without considering the friction, as displayed in Figure 4b.According to the failure morphology of the cubic concrete specimen in Figure 4c,d, 12 key recording joints (elements) were set in the computational model (Figure 4b).[25].

Solution and Verification of Ultimate Plastic Shear Strain
Given that there are numerous studies on the ultimate strain of concrete [23,24], the numerical limit analysis method is used to calculate and analyze the plastic shear strain of a C25 ordinary concrete cube specimen in this study.Hence, the ultimate plastic shear strain value of C25 ordinary concrete is obtained, and a comparison is made with the results of previous studies to verify the reliability of the numerical limit analysis method of plastic shear strain of C25 ordinary concrete.
In this research, the finite difference software FLAC3D of version 6.0 is adopted to numerically simulate the C25 concrete cube specimen.First, the ultimate load of C25 concrete is determined by the overload method.The calculation model is a cubic specimen with a side length of 150 mm, and each side is divided into 20 grids, as shown in Figure 4a.Afterwards, the strain-softened plasticity model is employed, whose material attribute parameters are defined.Hence, the physical and mechanical parameters of C25 concrete are exhibited in Table 1.Furthermore, the bottom surface of the calculation model is fixed, while the top surface is subjected to a downward vertical plane load without considering the friction, as displayed in Figure 4b.According to the failure morphology of the cubic concrete specimen in Figure 4c,d, 12 key recording joints (elements) were set in the computational model (Figure 4b).By means of the numerical simulation of the finite difference method, the main criteria to judge the failure include plastic zone penetration, the strength reduction method, the iteration non-convergence, etc. [26,27].In this study, the failure criterion, whether the average force ratio can reach the default value of the convergence standard, is adopted.In the meantime, the axial load is imposed on the calculation model step by step until it achieves the failure state.Currently, the prior-level load refers to the ultimate load of the model.
During the process of numerical simulation, the maximum unbalance force of the C25 concrete cube specimen model was sampled and recorded.Figure 5a manifests the maximum unbalance force curve of 24.772 MPa for the C25 concrete cube specimen.In this case, there is a trend of an exponential decline and a rapid convergence in the curve of the middle section and rear section, which reveals that the model is apparently undamaged.Simultaneously, Figure 5b shows the maximum unbalance force curve that rises sharply at the last stage and keeps fluctuating up and down of 24.773 MPa for a C25 concrete cube specimen, which demonstrates that the calculation result displays non-convergence and the model is damaged.Therefore, the ultimate load of the specimen, 24.772 MPa, is basically in accordance with the strength of C25 concrete, which verifies the accuracy and reliability of the physical and mechanical parameters of C25 concrete.By means of the numerical simulation of the finite difference method, the main criteria to judge the failure include plastic zone penetration, the strength reduction method, the iteration non-convergence, etc. [26,27].In this study, the failure criterion, whether the average force ratio can reach the default value of the convergence standard, is adopted.In the meantime, the axial load is imposed on the calculation model step by step until it achieves the failure state.Currently, the prior-level load refers to the ultimate load of the model.During the process of numerical simulation, the maximum unbalance force of the C25 concrete cube specimen model was sampled and recorded.Figure 5a manifests the maximum unbalance force curve of 24.772 MPa for the C25 concrete cube specimen.In this case, there is a trend of an exponential decline and a rapid convergence in the curve of the middle section and rear section, which reveals that the model is apparently undamaged.Simultaneously, Figure 5b shows the maximum unbalance force curve that rises sharply at the last stage and keeps fluctuating up and down of 24.773 MPa for a C25 concrete cube specimen, which demonstrates that the calculation result displays non-convergence and the model is damaged.Therefore, the ultimate load of the specimen, 24.772 MPa, is basically in accordance with the strength of C25 concrete, which verifies the accuracy and reliability of the physical and mechanical parameters of C25 concrete.In the finite difference software FLAC3D of version 6.0, the numerical limit analysis method is adopted to determine the ultimate plastic shear strain of C25 ordinary concrete.At that point, the basic steps can be presented as follows: (1) Carry the load to the ultimate one by using the function 'numpy.arange()' in the Python computer language and forming an arithmetic array; the function of the ultimate load, the maximum value of the plastic shear strain occurs at the upper part or corners of the model [22], and the damage to the upper side elements will inevitably lead to the overall failure of the specimen model.Therefore, the maximum value of the plastic shear strain is taken as the ultimate plastic shear strain of geotechnical materials.
Through the aforementioned basic steps, the relationship curve between the axial load and the plastic shear strain of 12 key elements in the C25 ordinary concrete cube specimen model is depicted in Figure 6.Among them, the plastic shear strains of units 9-12 become extremely small, and the plastic shear strain of unit 4 takes up the largest proportion, and failure occurs.Since unit 4 is located on the upper side of the calculation model, the damage to the element will inevitably contribute to the failure of the entire model.Thus, the plastic shear strain of this element signifies the ultimate plastic shear In the finite difference software FLAC3D of version 6.0, the numerical limit analysis method is adopted to determine the ultimate plastic shear strain of C25 ordinary concrete.At that point, the basic steps can be presented as follows: (1) Carry the load to the ultimate one by using the function 'numpy.arange()' in the Python computer language and forming an arithmetic array; (2) A cube specimen model with a size of 150 mm × 150 mm × 150 mm is constructed, whose sides are divided into 20 grids; (3) Assign the strain-softened plasticity model and set the corresponding material attribute parameters (including softening parameters); (4) Apply the displacement limits in all directions at the bottom of the model and impose the uniform load at the top side; (5) Extract the load from the arithmetic array and solve it to the equilibrium state; (6) According to the failure morphology of an ordinary concrete cube specimen, 12 key nodes (units) are set in the calculation model as monitoring points of plastic shear strain.Meanwhile, the coordinates of key nodes are loaded with the function 'numpy.array()' in the Python computer language.In addition, the function 'griddata( )' is used to extract the plastic shear strain values of 12 key nodes under each level of load.Finally, the function 'open( )' is adopted to output the values of each level of load and the plastic shear strain values of 12 key nodes; (7) Repeat the steps (1-6) and cycle to the end of the ultimate load; (8) Draw the 'axial load-plastic shear strain' curve of 12 key nodes (elements).Under the function of the ultimate load, the maximum value of the plastic shear strain occurs at the upper part or corners of the model [22], and the damage to the upper side elements will inevitably lead to the overall failure of the specimen model.Therefore, the maximum value of the plastic shear strain is taken as the ultimate plastic shear strain of geotechnical materials.
Through the aforementioned basic steps, the relationship curve between the axial load and the plastic shear strain of 12 key elements in the C25 ordinary concrete cube specimen model is depicted in Figure 6.Among them, the plastic shear strains of units 9-12 become extremely small, and the plastic shear strain of unit 4 takes up the largest proportion, and failure occurs.Since unit 4 is located on the upper side of the calculation model, the damage to the element will inevitably contribute to the failure of the entire model.Thus, the plastic shear strain of this element signifies the ultimate plastic shear strain, which equals 0.53‰.With respect to Equation ( 17), the two-dimensional ultimate plastic shear strain γ p f of an ordinary C25 concrete specimen is equivalent to 1.06‰.strain, which equals 0.53‰.With respect to Equation ( 17), the two-dimensional ultimate plastic shear strain    of an ordinary C25 concrete specimen is equivalent to 1.06‰.Under the condition of the principal strain sequence ε 1 ≤ ε 2 ≤ ε 3 , the numerical limit analysis method of plastic shear strain is employed to obtain the ultimate plastic shear strain value of C25 common concrete as 1.06‰.In addition, the ultimate plastic shear strain value of C25 common concrete, as given in the literature [22], equals 1.08‰.Thus, the results of these two studies remain approximately identical, which has verified the reliability of the numerical limit analysis method of plastic shear strain.

Definition of the Information Entropy Based on the Failure Approach Index
Assuming that the rock mass engineering system is equipped with n elements of given topological forms, boundary conditions, and applied loads, the failure approach index of the i-th element is   .Hence, the total failure approach index of the rock mass engineering system  can be expressed as: constructs: Therefore, there are: ,   ≥ 0 ( = 1,2, ⋯ , ).
Remarkably, the mechanical meaning of the complete and non-negative physical quantity   represents the proportion of the failure approach index of the i-th element in the total one; that is, the distribution of the failure approach index in the rock mass en- Under the condition of the principal strain sequence ε 1 ≤ ε 2 ≤ ε 3 , the numerical limit analysis method of plastic shear strain is employed to obtain the ultimate plastic shear strain value of C25 common concrete as 1.06‰.In addition, the ultimate plastic shear strain value of C25 common concrete, as given in the literature [22], equals 1.08‰.Thus, the results of these two studies remain approximately identical, which has verified the reliability of the numerical limit analysis method of plastic shear strain.

Definition of the Information Entropy Based on the Failure Approach Index
Assuming that the rock mass engineering system is equipped with n elements of given topological forms, boundary conditions, and applied loads, the failure approach index of the i-th element is q i .Hence, the total failure approach index of the rock mass engineering system Q can be expressed as: constructs: Therefore, there are: Remarkably, the mechanical meaning of the complete and non-negative physical quantity λ i represents the proportion of the failure approach index of the i-th element in the total one; that is, the distribution of the failure approach index in the rock mass engineering system can be described by λ i .In order to comprehensively reflect the disorder indicated by the local state variables of the whole system, the information entropy function based on failure approach index s of the rock mass engineering system is defined as: where φ symbolizes a positive constant and normally equals 1.

Evolution Equation of Synergetic Theory of Information Entropy Based on Failure Approach Index
In order to facilitate the calculation, two variables (s, x) are selected in this study to describe the stope evolution process of the surrounding rock system in underground mining space.In accordance with the synergetic theory, the evolution equation of the system can generally be written in the form of the Langevin equation [28]: where s signifies the slow variable, x denotes the fast variable, and K(s, x) depicts a nonlinear function containing the slow variable and the fast variable; l represents the excavation step, and F(l) refers to the fluctuation force.Equation ( 21) demonstrates that the evolution process of any nonlinear system is controlled by the internal factors of the system, which can be expressed by the nonlinear function K(s, x).On the other hand, it is also affected by external random factors, which can be replaced by the fluctuation force F(l).
Indeed, system evolution lies at the root of the internal factors, whereas the role of external factors is mainly manifested in promoting variations in the internal factors and triggering them at critical points of qualitative change.To be more specific, the fluctuation force F(l) belongs to an external factor that has no decisive effect on the stability of the surrounding rock system in underground mining space, so the fluctuation term F(l) in Equation ( 21) can be ignored.Regarding the two-dimensional systems, the nonlinear function in Equation ( 21) can be represented as: According to Equations ( 21) and ( 22), it can be realized as: .
By eliminating the fast variable x from Equation ( 23), the following integral expression can be obtained: The integration by parts method can transform x(l) in Equation ( 24) into a function of s(l), namely: when s changes slowly, it can be treated as a small quantity.By ignoring the integral term in Equation ( 25), we can obtain: where β signifies the growth rate of a fast variable, b = 1/β.
By substituting Equation ( 26) into Equation (23), the evolution equation of information entropy based on failure approach index and stope excavation step of surrounding rock systems in underground mining space can be attained, namely: The general solution of Equation ( 27) can be derived as follows: where a and b represent the undetermined coefficients, C indicates any constant, and l depicts the excavation step.
Since the influence of the fluctuation term F(l) is ignored, there are some errors in Equation ( 27).In order to reduce the effect of such errors, the gray system method is employed to accumulatively process the original sequence data of information entropy based on the failure approach index and dilute the impact of random factors (the fluctuation terms) on the original sequence data of information entropy [29].Let s be the non-negative sequence of the original information entropy based on the failure approach index, and the new sequence generated after one-time accumulation refers to s , that is: where n denotes the number of excavation steps.The expression for the new sequence generated by a single summation is: New sequence data created by mean values can be calculated as follows: Hence, Equation ( 27) can be written as: The values of a and b in Equation ( 33) can be obtained by using the least squares estimation method: where: The specific solution of Equation ( 33) can be attained as follows: In relation to Equation (33), the right-side term actually represents the change rate of information entropy.Therefore, the excavation step corresponding to the point of maximum change rate is taken as the unstable excavation step of the surrounding rock system in underground mining space, and construction: When dP ds = 0, at this point, P can get a maximum value, and the information entropy corresponding to its maximum value is: By combining Equations ( 37) and ( 39), the excavation step corresponding to the instability of the surrounding rock system in underground mining space can be acquired, namely: where l 1 represents the initial sequence number of excavation steps, which generally equals one and stays constant.
The expression of the safety factor of the surrounding rock system in underground mining space can be expressed as follows: where η signifies the safety factor; l denotes the unstable excavation step in the synergetic theory of information entropy based on the failure approach index; and L symbolizes the actual excavation step of the problem.If η > 1, the surrounding rock system of underground mining space is in a stable state; if η ≤ 1, it is in an unstable state.

Project Overview
Indeed, the main mining coal seams of Xutuan Coal Mine, the 7 1 coal seam and the 7 2 coal seam, generally exhibit the relationship of bifurcation in the superficial parts and merger in the deep regions, as displayed in Figure 7.In this case, there is an extensive variation between the spacing of the two coal layers, and the morphological characteristics of bifurcation and merging are extremely irregular.Occasionally, the conditions of bifurcation and merger will appear in the two ends and center sections of a working face, respectively.Therefore, the thickness of the coal seam varies tremendously, which makes the selection of the stoping technology exceedingly difficult.
The basic roof of 7 1 coal is fine sandstone with a fine grain structure, mainly quartz minerals.The direct roof of 7 1 coal is gray or dark gray mudstone.7 1 coal is observed to be black and lumpy.The direct floor of 7 1 coal (and the direct roof of 7 2 coal) is gray or dark gray mudstone, containing fossil plant roots.The coal of 7 2 is explored to be black and massive with endogenous fracture development, which generally contains a layer of gangue in the coal seam.The direct floor of 7 2 coal is gray or dark gray mudstone, and a coal streak is developed in the lower part.The basic floor of 7 2 coal is fine sandstone with gentle wavy horizontal bedding, primarily composed of quartz minerals.
Additionally, the average thicknesses of 7 1 and 7 2 coal refer to 1.9 m and 2.5 m, respectively.For the combined area of the coal seams 7 1 and 7 2 , that is, the rock strata with the thickness of gangue, which is less than 1 m, the coal mining technology of full-seam mining is adopted.However, in terms of the area where the gangue is thick, full-seam mining is not suitable as a consequence of the support failure, which needs to be mined by layers.In this case, the roof of the 7 2 coal seam is abnormally broken due to the impact of mining the 7 1 coal seam on the floor.During the early stage of stoping 7 2 14 fully mechanized mining faces on site are affected by the mining of 7 1 coal seams and the erosion of water, and the strength and bearing capacity of the rock mass of the roof of 7 2 coal seams are decreased, which will lead to a broken roof and poor integrity of the working face through the process of stoping.Meanwhile, the high fraction of coal wall caving and serious roof collapse of the end face during the mining will arise, and finally, the working face is obliged to be closed in advance, which has drastically influenced the normal replacement of mines.Therefore, how to ensure the safe and efficient mining of the underlying layered 7 2 coal in the case of slicing mining has become the major problem that the mine is facing.In this study, taking the 7 1 18 goaf of Xutuan Coal Mine as the engineering background, the influence of 7 1 coal mining on the stability of the floor surrounding the rock system is analyzed by means of the synergetic theory of information entropy based on the failure approach index, and the grouting reinforcement range of the local damaged area of the underlying 7 2 coal seam roof is determined, providing a guarantee for the safe and efficient mining of 7 2 coal seam.
where  1 represents the initial sequence number of excavation steps, which generally equals one and stays constant.
The expression of the safety factor of the surrounding rock system in underground mining space can be expressed as follows: where  signifies the safety factor;  ′ denotes the unstable excavation step in the synergetic theory of information entropy based on the failure approach index; and  symbolizes the actual excavation step of the problem.If  > 1, the surrounding rock system of underground mining space is in a stable state; if  ≤ 1, it is in an unstable state.

Project Overview
Indeed, the main mining coal seams of Xutuan Coal Mine, the 71 coal seam and the 72 coal seam, generally exhibit the relationship of bifurcation in the superficial parts and merger in the deep regions, as displayed in Figure 7.In this case, there is an extensive variation between the spacing of the two coal layers, and the morphological characteristics of bifurcation and merging are extremely irregular.Occasionally, the conditions of bifurcation and merger will appear in the two ends and center sections of a working face, respectively.Therefore, the thickness of the coal seam varies tremendously, which makes the selection of the stoping technology exceedingly difficult.The basic roof of 71 coal is fine sandstone with a fine grain structure, mainly quartz minerals.The direct roof of 71 coal is gray or dark gray mudstone.71 coal is observed to be black and lumpy.The direct floor of 71 coal (and the direct roof of 72 coal) is gray or dark gray mudstone, containing fossil plant roots.The coal of 72 is explored to be black and massive with endogenous fracture development, which generally contains a layer of gangue in the coal seam.The direct floor of 72 coal is gray or dark gray mudstone, and a coal streak is developed in the lower part.The basic floor of 72 coal is fine sandstone with gentle wavy horizontal bedding, primarily composed of quartz minerals.Additionally, the average thicknesses of 71 and 72 coal refer to 1.9 m and 2.5 m, respectively.For the combined area of the coal seams 71 and 72, that is, the rock strata with the thickness of gangue, which is less than 1 m, the coal mining technology of full-seam mining is adopted.However, in terms of the area where the gangue is thick, full-seam mining is not suitable as a consequence of the support failure, which needs to be mined by layers.In this case, the roof of the 72 coal seam is abnormally broken due to the impact of mining the 71 coal seam on the floor.During the early stage of stoping 7214 fully mechanized mining faces on site are affected by the mining of 71 coal seams and the erosion of water, and the strength and bearing capacity of the rock mass of the roof of 72 coal seams are decreased, which will lead to a broken roof and poor integrity of the working face through the process of stoping.Meanwhile, the high fraction of coal wall caving and serious roof collapse of the end face during the mining will arise, and finally, the working face is obliged to be closed in advance, which has drastically influenced the normal replacement of mines.Therefore, how to ensure the safe and efficient mining of the underlying layered 72 coal in the case of slicing mining has become the major problem that the mine is facing.In this study, taking the 7118 goaf of Xutuan Coal Mine as the engineering background, the influence of 71 coal mining on the stability of the floor surrounding the rock system is analyzed by means of the synergetic theory of information entropy based on the failure approach index, and the grouting reinforcement range of the local damaged area of the underlying 72 coal seam roof is determined, providing a guarantee for the safe and efficient mining of 72 coal seam.To illustrate, the coal and rock sampling sites are distributed on the roof of the air roadway in the 7225 working face, the falling gangue in the goaf, and the 71 working face.The gangue was cored, sliced, and polished into standard cylindrical specimens through the necessary machining procedures, and the mechanical tests, such as the compressive and tensile strengths of coal and rock, were carried out by the RMT rock mechanics test system, as displayed in Figure 9. Hence, the physical and mechanical parameters of coal and rock have been obtained, as listed in Table 2. To illustrate, the coal and rock sampling sites are distributed on the roof of the air roadway in the 7 2 25 working face, the falling gangue in the goaf, and the 7 1 working face.The gangue was cored, sliced, and polished into standard cylindrical specimens through the necessary machining procedures, and the mechanical tests, such as the compressive and tensile strengths of coal and rock, were carried out by the RMT rock mechanics test system, as displayed in Figure 9. Hence, the physical and mechanical parameters of coal and rock have been obtained, as listed in Table 2.

Stability Analysis of the Floor Surrounding Rock System of 71 Coal Seam
In light of the aforementioned numerical limit analysis method of plastic shear strain, the ultimate plastic shear strain value of mudstone in the floor of No. 71 coal seam in Xutuan Coal Mine is determined to be 0.14‰, which provides a crucial criterion for the calculation of information entropy based on the failure approach index.According to the literature [30], the roof pressure step distance of No. 71 coal seam mining in Xutuan Coal Mine is within the range of 8-12 m, so the model excavation range equals 15 m.Simultaneously, considering that the 7118 working face is affected by a periodic weighting, the synergetic theory of information entropy based on the failure approach index is adopted to analyze the stability of the floor surrounding the rock structure after stoping on the 7118 working face.
With the stoping of the working face 7118, the values of the failure approach index of the floor surrounding rock unit below the goafs are extracted 1 m per advance, and the information entropy based on the failure approach index of the floor surrounding rock system is calculated by Equation (20).The sequence data is listed in Table 3.In light of the aforementioned numerical limit analysis method of plastic shear strain, the ultimate plastic shear strain value of mudstone in the floor of No. 7 1 coal seam in Xutuan Coal Mine is determined to be 0.14‰, which provides a crucial criterion for the calculation of information entropy based on the failure approach index.According to the literature [30], the roof pressure step distance of No. 7 1 coal seam mining in Xutuan Coal Mine is within the range of 8-12 m, so the model excavation range equals 15 m.Simultaneously, considering that the 7 1 18 working face is affected by a periodic weighting, the synergetic theory of information entropy based on the failure approach index is adopted to analyze the stability of the floor surrounding the rock structure after stoping on the 7 1 18 working face.
With the stoping of the working face 7 1 18, the values of the failure approach index of the floor surrounding rock unit below the goafs are extracted 1 m per advance, and the information entropy based on the failure approach index of the floor surrounding rock system is calculated by Equation (20).The sequence data is listed in Table 3.In virtue of the sequence data of the original information entropy based on the failure approach index in Table 3, a = 0.1887 and b = 0.00023744 are calculated in accordance with Equations (34)-(36).Afterwards, l = 15.6922 is computed in relation to Equation (40).In this way, the safety factor of the floor surrounding the rock structure equals 1.0461.Therefore, the floor surrounding the rock system of the 7 1 coal seam proves to be firm.
Although the synergetic theory of information entropy based on the failure approach index has calculated that the floor surrounding rock system of 7 1 coal seams performs steady from a macroscopic perspective, the local area of floor surrounding rock is damaged with respect to the microscopic rock mechanics, as displayed in Figure 10.According to the partition significance of failure approach index, when FAI ≥ 2.0, the rock and soil mass is in a state of failure [14].When advancing 0-6 m above the 7 1 18 working face, there is no damage to the floor surrounding the rock.Nonetheless, when the working face is excavated to 6 m, the local failure of the floor surrounding the rock will emerge, and the maximum failure depth will reach 2.3 m.With the advance of the 7 1 18 coal mining work face, the abutment pressure of the stoping work face in front and behind becomes increasingly evident.When the working face is excavated to 12 m, the local failure depth of the floor surrounding the rock will reach 5.6 m.At this moment, the failure depth of the floor surrounding the rock grows to a peak of 5.6 m and tends to be stable, as exhibited in Figure 11.Consequently, it is consistent with the reference [31] that the depth of floor failure caused by the mining of 7 1 coal is no more than 6 m.In view of the local failure area's enormous stability influences, advanced grouting reinforcement measures should be taken to ensure the safe and efficient stoping of 7 2 coal seams.In virtue of the sequence data of the original information entropy based on the failure approach index in Table 3, a = 0.1887 and b = 0.00023744 are calculated in accordance with Equations (34), (35), and (36).Afterwards,  ′ = 15.6922 is computed in relation to Equation (40).In this way, the safety factor of the floor surrounding the rock structure equals 1.0461.Therefore, the floor surrounding the rock system of the 71 coal seam proves to be firm.
Although the synergetic theory of information entropy based on the failure approach index has calculated that the floor surrounding rock system of 71 coal seams performs steady from a macroscopic perspective, the local area of floor surrounding rock is damaged with respect to the microscopic rock mechanics, as displayed in Figure 10.According to the partition significance of failure approach index, when FAI ≥ 2.0, the rock and soil mass is in a state of failure [14].When advancing 0-6 m above the 7118 working face, there is no damage to the floor surrounding the rock.Nonetheless, when the working face is excavated to 6 m, the local failure of the floor surrounding the rock will emerge, and the maximum failure depth will reach 2.3 m.With the advance of the 7118 coal mining work face, the abutment pressure of the stoping work face in front and behind becomes increasingly evident.When the working face is excavated to 12 m, the local failure depth of the floor surrounding the rock will reach 5.6 m.At this moment, the failure depth of the floor surrounding the rock grows to a peak of 5.6 m and tends to be stable, as exhibited in Figure 11.Consequently, it is consistent with the reference [31] that the depth of floor failure caused by the mining of 71 coal is no more than 6 m.In view of the local failure area's enormous stability influences, advanced grouting reinforcement measures should be taken to ensure the safe and efficient stoping of 72 coal seams.

Conclusions
(1) Based on the strain-softened plasticity model, the general relationship between the total two-dimensional plastic shear strain and the total three-dimensional plastic shear strain has been deduced, and the numerical limit analysis method of plastic shear strain has been proposed to obtain the threshold of the ultimate plastic shear strain, which has provided the crucial criterion for the calculation of information entropy based on the failure approach index; (2) The numerical limit analysis method was used to calculate and analyze the plastic shear strain of C25 ordinary concrete specimens, and the ultimate plastic shear strain value obtained was consistent with the result of previous studies, which verified the correctness and reliability of the numerical limit analysis method of plastic shear strain; (3) With regard to the synergetic theory of the principle of maximum information entropy, the evolution equation of the information entropy function based on the failure approach index and excavation step of the surrounding rock system in underground mining space is established.Meanwhile, the general and solution formulae and the expression of the unstable excavation step are attained, which has presented a theoretical analysis approach for the stability analysis of the surrounding rock system;

Conclusions
(1) Based on the strain-softened plasticity model, the general relationship between the total two-dimensional plastic shear strain and the total three-dimensional plastic shear strain has been deduced, and the numerical limit analysis method of plastic shear strain has been proposed to obtain the threshold of the ultimate plastic shear strain, which has provided the crucial criterion for the calculation of information entropy based on the failure approach index; (2) The numerical limit analysis method was used to calculate and analyze the plastic shear strain of C25 ordinary concrete specimens, and the ultimate plastic shear strain value obtained was consistent with the result of previous studies, which verified the correctness and reliability of the numerical limit analysis method of plastic shear strain; (3) With regard to the synergetic theory of the principle of maximum information entropy, the evolution equation of the information entropy function based on the failure approach index and excavation step of the surrounding rock system in underground mining space is established.Meanwhile, the general and solution formulae and the expression of the unstable excavation step are attained, which has presented a theoretical analysis approach for the stability analysis of the surrounding rock system;

Conclusions
(1) Based on the strain-softened plasticity model, the general relationship between the total two-dimensional plastic shear strain and the total three-dimensional plastic shear strain has been deduced, and the numerical limit analysis method of plastic shear strain has been proposed to obtain the threshold of the ultimate plastic shear strain, which has provided the crucial criterion for the calculation of information entropy based on the failure approach index; (2) The numerical limit analysis method was used to calculate and analyze the plastic shear strain of C25 ordinary concrete specimens, and the ultimate plastic shear strain value obtained was consistent with the result of previous studies, which verified the correctness and reliability of the numerical limit analysis method of plastic shear strain; (3) With regard to the synergetic theory of the principle of maximum information entropy, the evolution equation of the information entropy function based on the failure approach index and excavation step of the surrounding rock system in underground mining space is established.Meanwhile, the general and solution formulae and the expression of the unstable excavation step are attained, which has presented a theoretical analysis approach for the stability analysis of the surrounding rock system; (4) Taking the close coal seam mining project in Xutuan Coal Mine as an example, the synergetic theory of information entropy based on the failure approach index was used to analyze whether the surrounding rock system of the 7 1 coal bottom plate with a layer thickness of 6 m is stable.Considering the periodic weighting of 7 1 coal roofs, the local damage depth of the mining floor reaches 5.6 m.In addition, the appropriate technical measures for controlling the stability of the underlying 7 2 coal mining roofs have been put forward to optimize the integrity of the stope roof and ensure the safety and efficiency of coal extraction.
ij ; and γ p f stands for the ultimate plastic shear strain.

Figure 1 .
Figure 1.Geographical position of Xutuan Coal Mine and column diagram of coal and rock layers.

Figure 1 .
Figure 1.Geographical position of Xutuan Coal Mine and column diagram of coal and rock layers.

Figure 2 .
Figure 2. Numerical simulation flow chart of the ultimate load.

Figure 2 .
Figure 2. Numerical simulation flow chart of the ultimate load.

Figure 3 .
Figure 3. Numerical calculation flow chart of the ultimate plastic shear strain.

Figure 3 .
Figure 3. Numerical calculation flow chart of the ultimate plastic shear strain.

Figure 4 .
Figure 4. Calculation model: (a) geometric model; (b) boundary conditions and key nodes (units): key node (unit) 1~12 represent the monitoring node (unit) of plastic shear strain ; (c) plane figure of failure morphology; and (d) real picture of failure morphology.

Figure 4 .
Figure 4. Calculation model: (a) geometric model; (b) boundary conditions and key nodes (units): key node (unit) 1~12 represent the monitoring node (unit) of plastic shear strain; (c) plane figure of failure morphology; and (d) real picture of failure morphology.

Figure 5 .
Figure 5. Curve of maximum unbalanced force: (a) area load on top surface with 24.772 MPa; (b) area load on top surface with 24.773 MPa.
(2) A cube specimen model with a size of 150 mm × 150 mm × 150 mm is constructed, whose sides are divided into 20 grids; (3) Assign the strain-softened plasticity model and set the corresponding material attribute parameters (including softening parameters); (4) Apply the displacement limits in all directions at the bottom of the model and impose the uniform load at the top side; (5) Extract the load from the arithmetic array and solve it to the equilibrium state; (6) According to the failure morphology of an ordinary concrete cube specimen, 12 key nodes (units) are set in the calculation model as monitoring points of plastic shear strain.Meanwhile, the coordinates of key nodes are loaded with the function 'numpy.array()' in the Python computer language.In addition, the function 'griddata( )' is used to extract the plastic shear strain values of 12 key nodes under each level of load.Finally, the function 'open( )' is adopted to output the values of each level of load and the plastic shear strain values of 12 key nodes; (7) Repeat the steps (1-6) and cycle to the end of the ultimate load; (8) Draw the 'axial load-plastic shear strain' curve of 12 key nodes (elements).Under

Figure 5 .
Figure 5. Curve of maximum unbalanced force: (a) area load on top surface with 24.772 MPa; (b) area load on top surface with 24.773 MPa.

Figure 6 .
Figure 6.Relational graph between the axial load and the plastic shear strain.

Figure 6 .
Figure 6.Relational graph between the axial load and the plastic shear strain.

Figure 7 .
Figure 7. Schematic diagram of the bifurcation and merge of 71 and 72 coal in the Xutuan Coal Mine.

Figure 7 .
Figure 7. Schematic diagram of the bifurcation and merge of 7 1 and 7 2 coal in the Xutuan Coal Mine.Furthermore, the 7 2 coal seam belongs to the lower Shihezi Formation of the lower Permian System with a thickness of 1.7 m-3.7 m and an average thickness of 2.5 m, which reflects the stable occurrence and has a spacing of 0.8 m to 9.3 m from the upper 7 1 coal seam.The spatial relationship between the 7 1 18 goaf and the 7 2 25 working face discussed in this study is depicted in Figure8.The elevation of the 7 1 18 goaf is −447.9m~−480.6 m, and the ground elevation is +25.0 m.

Furthermore, the
72 coal seam belongs to the lower Shihezi Formation of the lower Permian System with a thickness of 1.7 m-3.7 m and an average thickness of 2.5 m, which reflects the stable occurrence and has a spacing of 0.8 m to 9.3 m from the upper 71 coal seam.The spatial relationship between the 7118 goaf and the 7225 working face discussed in this study is depicted in Figure 8.The elevation of the 7118 goaf is −447.9m~−480.6 m, and the ground elevation is +25.0 m.

Figure 8 .
Figure 8.The spatial position relationship between the 7118 goaf and the 7225 working face.

Figure 11 .
Figure 11.The curve of excavation distance from the 7118 working face and failure depth from the local area of the floor surrounding the rock.

Figure 11 .
Figure 11.The curve of excavation distance from the 7118 working face and failure depth from the local area of the floor surrounding the rock.

Figure 11 .
Figure 11.The curve of excavation distance from the 7 1 18 working face and failure depth from the local area of the floor surrounding the rock.

Table 2 .
Physical and mechanical parameters of each rock layer.

Table 3 .
Sequence data of information based on the failure approach index.

Table 2 .
Physical and mechanical parameters of each rock layer.

Table 3 .
Sequence data of information entropy based on the failure approach index.