Rock creep disturbance induced instability and post-disturbance residual life prediction method

Abstract

The invention discloses a rock creep disturbance induced instability and post-disturbance residual life prediction method. The method comprises the following steps: constructing a creep damage factor and a disturbance damage factor based on irreversible deformation; the coupling effect of creep damage and disturbance damage is comprehensively considered, a theoretical formula capable of analyzing rock creep response in a complex disturbance environment is obtained, and the problem of nonlinear response prediction of the rock creep process under the disturbance effect is creatively solved. According to the method, the disturbance aftereffect in the rock creep process can be scientifically represented, and disturbance-induced instant instability and lag instability can be predicted; according to the method, the combined elements with irreversible disturbance deformation are described in series, nonlinear correction is carried out, damage evolution and element parameters are dynamically coupled, a composite element model considering the historical disturbance cumulative effect is established, and the memory characteristic of rock mass deformation on the historical disturbance effect can be described.

Description

Technical Field

[0001] This invention belongs to the field of stability analysis of deep underground engineering tunnels, specifically involving a method for predicting rock creep disturbance-induced instability and post-disturbance remaining life. Background Technology

[0002] Deep underground engineering rock masses are often in a relatively stable time-deformation process due to significant stress. However, the time-deformation process can be accelerated or even lead to instability when subjected to various disturbances such as tectonic movements, engineering activities, and underground explosions. Revealing the development and evolution of time-deformation in engineering rock masses under weak disturbances and the physical mechanisms underlying disasters is crucial for optimizing time-delay disaster prevention and control technologies and improving the long-term stability and resilience of surrounding rock. Clarifying this process requires constructing theoretical models, decoupling the interaction between creep and disturbance factors on the stability of surrounding rock, determining relevant instability thresholds, and enabling early warning and proactive prevention of time-delay disasters.

[0003] Based on damage theory, current combined element models often characterize damage and describe accelerated creep stages by reducing model parameters over time or by connecting time-dependent nonlinear elements in series. However, two key scientific bottlenecks remain: First, damage variables and damage thresholds are highly correlated with time parameters and rely on empirical posterior calibration, failing to consider damage changes caused by various disturbances and their impact on subsequent creep processes. Second, creep damage and disturbance damage cannot be decoupled, making it difficult to quantify the contribution of disturbances to the long-term stability of engineering surrounding rock. These shortcomings significantly hinder the extension of such models to the study of the delayed instability process of engineering surrounding rock under weak disturbance and creep coupling. Summary of the Invention

[0004] The purpose of this invention is to provide a method for predicting rock creep disturbance-induced instability based on irreversible deformation. This method constructs creep damage factors and disturbance damage factors based on irreversible deformation, comprehensively considering the coupling effect of creep damage and disturbance damage, to obtain a theoretical formula that can analyze the rock creep response under complex disturbance environments. This innovatively solves the problem of predicting the nonlinear response of rock creep processes under disturbance. This invention can scientifically characterize the disturbance aftereffects of rock creep processes, that is, the creep process after disturbance exhibits a different state or behavior than the undisturbed state, and predict disturbance-induced immediate and delayed instability. This invention is modular and combinable, and can also degenerate to describe conventional creep instability.

[0005] The technical solution to achieve the objective of this invention is: a method for predicting rock creep disturbance-induced instability and post-disturbance remaining life, comprising the following steps:

[0006] S1: Based on the irreversible deformation evolution law of rocks, the plastic element parameters are regarded as variable parameters related to the irreversible creep deformation of rocks, and a rock creep damage function D is constructed.c ;

[0007] S2: Based on the creep deformation composition of rocks, a Burgers model containing elastic, viscoelastic, and plastic elements is constructed, and the rock creep damage function D is applied. c Substituting into the Burgers model, a rock creep equation ε is established, in which the parameters of the plastic elements are iteratively reduced as irreversible creep deformation increases. c ;

[0008] S3: Treating the component parameters as variable parameters related to the irreversible deformation of rock disturbance, constructing the rock disturbance damage function D. d ;

[0009] S4: The rock disturbance damage function D d Substituting into the Burgers model, a perturbation damage equation is established to describe the evolution of irreversible deformation during rock perturbation instability.

[0010] S5: Based on the deformation decomposition process of rock creep instability under disturbance, the disturbance damage equation is... The element describing the irreversible deformation of the initial disturbance is replaced by the rock creep equation ε. c The paper describes the damage plastic element of irreversible creep deformation and establishes the perturbation damage equation for the evolution of irreversible rock deformation coupled with creep load.

[0011] S6: Applying the rock creep equation ε c The damage plastic element describing irreversible creep deformation is replaced by a perturbation damage equation describing the evolution of irreversible rock deformation coupled with creep load. The creep-disturbance coupled nonlinear plastic body was characterized, and the model equation ε of the rock creep-disturbance hysteresis instability combined element was established.

[0012] S7: Based on the rock creep disturbance hysteresis instability combined element model equation ε, and combined with the rock instability irreversible deformation criterion, determine the rock instability type and state, and predict the rock's post-disturbance instability time t. f ;

[0013] S8: Based on the rock's post-disturbance instability time t f In conjunction with creep lifetime, the remaining lifetime of the rock after disturbance is calculated.

[0014] Furthermore, step S1 specifically includes:

[0015] S11: Based on the deformation composition and rock creep instability deformation threshold of the Burgers model, the allowable irreversible deformation ε of the rock is obtained. a for:

[0016]

[0017] In the formula, ε a For the permissible irreversible deformation of rocks, ε represents the total creep deformation when the rock completely breaks down. e For the instantaneous elastic deformation of rock, ε ve This refers to the viscoelastic deformation of the rock;

[0018] S12: Define the rock creep damage function D c Therefore, creep is an irreversible deformation. An exponential function of variables is expressed as:

[0019]

[0020] Ignoring the instantaneous irreversible deformation caused by the application of creep load, when At that time, D c =0; while when At that time, D c =1.

[0021] Furthermore, in step S2, the rock creep equation ε, in which the element parameters are iteratively reduced with increasing irreversible creep deformation, is further refined. c The construction method includes the following steps:

[0022] S21: In the Burgers model, the viscosity coefficient η1 of the plastic element decreases due to the continuous effect of damage, expressed as:

[0023] η1(D c )=(1-D c )η1;

[0024] In the formula, η1(D c () represents the viscosity coefficient considering creep damage reduction;

[0025] S22: Replace the viscosity coefficient η1 of the plastic element in the Burgers model with a viscosity coefficient η1(D) that takes into account creep damage reduction. c This yields the rock creep equation ε, which iteratively reduces the plastic element parameters as irreversible creep deformation increases. c The expression:

[0026]

[0027] Furthermore, in step S3, the rock disturbance damage function D is constructed. d The method includes the following steps:

[0028] S31: Defines the irreversible perturbation deformation ε when rock completely breaks down. b The expression is:

[0029]

[0030] In the formula, The total disturbance and irreversible deformation when the rock completely breaks down. The initial disturbance results in irreversible deformation. This represents the cumulative irreversible deformation increment during disturbance instability.

[0031] S32: Define the rock disturbance damage function D d Therefore, the disturbance causes irreversible deformation. Exponential function with variable:

[0032]

[0033] In the formula, β is the disturbance damage factor. The total irreversible deformation during the disturbance process is expressed as:

[0034]

[0035] In the formula, The initial disturbance results in irreversible deformation. This represents the cumulative irreversible deformation increment during the disturbance process. Before the disturbance, there is an initial irreversible deformation due to the disturbance. D d ≠0; if during the disturbance process

[0036] D d =1.

[0037] Furthermore, in step S4, the disturbance damage equation is derived from the evolution of irreversible deformation during the rock disturbance instability process. The construction method includes the following steps:

[0038] S41: Disturbance Damage Equation The irreversible deformation of the rock under disturbance is further decomposed, and all model elements are used to describe the irreversible deformation. When constructing the disturbance damage equation, except for the initial irreversible deformation described by the initial disturbance instantaneous modulus E3 which does not participate in nonlinear growth, all other model parameters are reduced, as expressed below:

[0039]

[0040] In the formula, η3(D d ), η4(D d E4(D) represents the viscosity coefficient considering disturbance damage reduction; d The viscoelastic modulus is calculated considering the reduction due to disturbance damage.

[0041] S42: Replace the viscosity coefficient η3 of the plastic element, the viscosity coefficient η4 of the viscoelastic element, and the viscoelastic modulus E4 in the Burgers model with the viscosity coefficient η3 (D) considering the reduction of disturbance damage. d ), η4(D d )and

[0042] E4(D d The disturbance damage equation for the irreversible deformation evolution during rock disturbance instability was obtained. The expression is:

[0043]

[0044] Furthermore, in step S5, the perturbation damage equation for the irreversible deformation evolution of rock coupled with creep load is... The construction method includes the following steps:

[0045] S51: Irreversible deformation at the interface of different stress states is continuous. The irreversible creep deformation that ultimately occurs during the creep process before disturbance is equal to the initial irreversible deformation during disturbance, expressed as:

[0046]

[0047] In the formula, For the irreversible creep deformation at the start time of the disturbance. The disturbance at the start time is an irreversible deformation;

[0048] Irreversible Deformation of Rock under Coupled Creep and Disturbance Loads It consists of two parts, and the expression is:

[0049]

[0050] After the disturbance ends, the rock continues to creep. The total irreversible deformation resulting from the disturbance process is equal to the initial irreversible creep deformation in the subsequent creep process, expressed as:

[0051]

[0052] In the formula, The total irreversible deformation of rock under the coupled action of creep load and disturbance load at the end time of the disturbance is defined as follows: The initial irreversible creep deformation at the end of the disturbance;

[0053] S52: Rock creep damage function D under coupled creep load and disturbance load c Disturbance damage function D d Total irreversible deformation of rock under the coupled action of creep load and disturbance load Control, the expression is:

[0054]

[0055] S53: Extending the cumulative irreversible deformation during the perturbation process to a continuous function in the time domain t allows general creep damage and perturbation damage to be coupled in the same model. The expression for the number of perturbations n is:

[0056]

[0057] In the formula, Δt is the duration of the disturbance load, and its value is (tt) x ), t x Let T be the disturbance initiation time, T be the disturbance period, and Δt / T be defined as the generalized disturbance number. Then, the irreversible deformation of the disturbance will continue to develop in the time domain t.

[0058] S54: Based on the continuity condition of irreversible deformation, the perturbation damage equation for the evolution of irreversible deformation of rock coupled with creep load. The expression is:

[0059]

[0060] Furthermore, in step S6, based on different rock loading states, the perturbation damage equation for the irreversible deformation evolution of rock coupled with creep load is... There are three modes, and the construction method includes the following steps:

[0061] S61: Mode 1: The expression for the rock creep disturbance hysteresis instability combined element model equation ε before the application of disturbance load is the same as the rock creep equation ε in step S2.2. c Same expression:

[0062]

[0063] S62: Mode 2: The irreversible creep deformation described in step S22 The damaged plastic element is replaced by the perturbation damage equation for the irreversible deformation evolution of rock under coupled creep load as described in step S54. The creep-disturbance coupled nonlinear plastic body was characterized, and the model equation ε of the rock creep-disturbance hysteresis instability combined element under the coupled action of creep load and disturbance load was obtained, and the expression is:

[0064]

[0065] S63: After the disturbance ends, the rock continues to creep. The expression for the cumulative irreversible deformation caused by the coupling effect of creep load and disturbance load is:

[0066]

[0067] S64: Mode 3: The perturbation damage equation for the irreversible deformation evolution of rock coupled with creep load in step S62. Replace with the cumulative irreversible deformation caused by the coupling effect of creep load and disturbance load in step S63. The equation ε of the rock creep disturbance hysteresis instability combined element model after the application of disturbance load is obtained, and its expression is:

[0068]

[0069] Compared with the prior art, the significant advantages of this invention are:

[0070] Based on the Burgers model, this invention constructs a time-varying damage function with cumulative irreversible deformation as the internal variable. Both the creep damage function and the disturbance damage function depend on the current irreversible deformation and the unstable irreversible deformation. By connecting the combined elements that describe the disturbance irreversible deformation in series and performing nonlinear correction, the damage evolution is dynamically coupled with the element parameters to establish a composite element model that considers the cumulative effect of disturbance history. This model can describe the "memory" characteristic of rock mass deformation to the effect of historical disturbance. Attached Figure Description

[0071] Figure 1 This is a flowchart of a method for predicting rock creep disturbance-induced instability and post-disturbance remaining life.

[0072] Figure 2 This is a schematic diagram of the Burgers model.

[0073] Figure 3 This is a schematic diagram of the rock creep composite element model proposed in this invention.

[0074] Figure 4 This is a schematic diagram of the deformation decomposition of rock creep instability under disturbance.

[0075] Figure 5 This is a schematic diagram of the combined element model for irreversible deformation and damage caused by rock disturbance proposed in this invention.

[0076] Figure 6 This is a schematic diagram of the rock irreversible deformation disturbance damage combination element model with coupled creep load proposed in this invention.

[0077] Figure 7 This is a schematic diagram of the rock creep disturbance hysteresis instability combined element model proposed in this invention.

[0078] Figure 8 This is a comparison diagram of the measured data of Embodiment 1 of the present invention and the combined element model proposed in the present invention.

[0079] Figure 9This is a prediction result diagram of Example 1 based on the combined element model proposed in this invention.

[0080] Figure 10 This is a comparison diagram of the measured data of Embodiment 2 of the present invention and the combined element model proposed in the present invention. Detailed Implementation

[0081] The present invention will now be described in further detail with reference to the accompanying drawings.

[0082] To enhance understanding of the present invention, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. The following embodiments are only used to more clearly illustrate the technical solutions of the present invention and should not be construed as limiting the scope of protection of the present invention.

[0083] This invention solves the problem that existing technologies are unable to describe the hysteretic instability of rocks under the coupled effects of disturbance and creep.

[0084] like Figures 1 to 7 As shown, this invention provides a method for predicting rock creep disturbance-induced instability and post-disturbance remaining life, specifically including the following steps:

[0085] S1. Based on the irreversible deformation evolution law of rocks, the parameters of plastic elements are regarded as variable parameters related to the irreversible creep deformation of rocks, and a rock creep damage function D is constructed. c ;

[0086] S2, the rock creep damage function D c Substituting into the Burgers model, a rock creep equation ε is established, in which the parameters of the plastic elements are iteratively reduced as irreversible creep deformation increases. c ;

[0087] S3. Treat the component parameters as variable parameters related to the irreversible deformation of rock disturbance, and construct the rock disturbance damage function D. d ;

[0088] S4. The rock disturbance damage function D d Substituting into the Burgers model, a perturbation damage equation is established to describe the evolution of irreversible deformation during rock perturbation instability.

[0089] S5. Based on the deformation decomposition process of rock creep instability under disturbance, the disturbance damage equation is... The element describing the irreversible deformation of the initial disturbance is replaced with the rock creep equation ε. c The paper describes the damage plastic element of irreversible creep deformation and establishes the perturbation damage equation for the evolution of irreversible rock deformation coupled with creep load.

[0090] S6. The rock creep equation ε c The damage plastic element describing irreversible creep deformation is replaced by the perturbation damage equation for the evolution of irreversible rock deformation coupled with creep load. The creep-disturbance coupled nonlinear plastic body was characterized, and the model equation ε of the rock creep-disturbance hysteresis instability combined element was established.

[0091] like Figure 2 As shown, the Burgers model consists of three deformed parts:

[0092]

[0093] The Burgers model can describe decelerating and steady creep of rocks well, but its components are all linear, and it cannot describe accelerated creep or nonlinear creep evolution under disturbance. Under continuous creep load, the accumulation and development of damage deteriorates the deformation parameters and eventually leads to accelerated creep and instability. Based on this, as... Figure 3 As shown, the parameters of the plastic elements describing irreversible deformation in the Burgers model are nonlinearly modified to construct a nonlinear Burgers volume that reflects accelerated rock creep. The nonlinear Burgers volume consists of three deformable parts:

[0094]

[0095] Constructed rock creep equation ε c Presented in implicit function form, it includes both the forward problem of predicting results from the model and the "inverse problem" of determining model parameters based on observations. Iterative calculations are performed using Matlab, and the expression is:

[0096]

[0097] In the formula, η is the viscosity coefficient considering creep damage reduction in the i-th time step. 1(i) (D c Irreversible deformation resulting solely from the previous time step Sure.

[0098] like Figure 4 As shown, for a specific disturbance, the deformation response caused by the disturbance can be divided into different components according to its nature. The deformation before the disturbance follows the general creep law, consisting of instantaneous elastic deformation ε. e Viscoelastic deformation ε ve and irreversible deformation of creep Composition; the deformation increment during the disturbance process can also be decomposed into instantaneous elastic deformation increment. Viscoelastic deformation increment and the cumulative irreversible deformation increment during the disturbance process After the disturbance ended, the rock returned to the creep process, but its irreversible deformation development was affected by the disturbance.

[0099] like Figure 5 As shown, the elastic element is used to describe the irreversible deformation of the rock under initial disturbance. Viscoelastic and plastic elements considering disturbance damage reduction are used to describe the cumulative irreversible deformation increment of rocks during disturbance. Disturbance damage equation for the evolution of irreversible deformation during rock disturbance instability It consists of two deformed parts:

[0100]

[0101] In Matlab, iterative calculations are performed using the following expression:

[0102]

[0103] like Figure 6 As shown, based on the continuity condition of irreversible deformation, the perturbation damage equation for the evolution of irreversible deformation of rock coupled with creep load is presented. It consists of two deformed parts:

[0104]

[0105] In Matlab, iterative calculations are performed using the following expression:

[0106]

[0107] like Figure 7 As shown, the damage plastic element describing irreversible creep deformation in the nonlinear Burgers volume is replaced with a creep-disturbance coupled nonlinear plastic body to establish a rock creep-disturbance hysteresis instability combined element model based on irreversible deformation.

[0108] Example 1

[0109] The rock creep disturbance test data of red sandstone in the literature "Wang Bo, et al. Experimental study on microscopic damage of rock mass in sensitive neighborhood under rheological disturbance conditions [J]. Chinese Journal of Rock Mechanics and Engineering, 2024, 43(S2):3820-3831" were used to apply the rock creep disturbance hysteresis instability combined element model proposed in this invention. The model parameters obtained by nonlinear least square fitting of the relevant experimental data are shown in Table 1.

[0110] Table 1. Parameters of the combined element model for rock creep disturbance hysteresis instability based on irreversible deformation.

[0111]

[0112] like Figure 8 As shown, the model proposed in this invention can well describe the deformation evolution process of rocks under creep-disturbance coupling, demonstrating the rationality and applicability of the model.

[0113] like Figure 9 As shown, by combining disturbance parameters and initial creep state, real-time prediction of creep instability time can be achieved, and the prediction results can be updated through in-situ monitoring data, overcoming the technical bottleneck of traditional time-dependent damage models that rely on posterior parameter calibration. This model innovatively separates the action mechanisms of creep and disturbance, establishing element modules to characterize creep damage and disturbance damage respectively. Furthermore, it achieves the coupling of creep and disturbance effects through cumulative irreversible deformation, overcoming the parameter calibration bias caused by excessive correlation between damage parameters and time parameters in traditional models.

[0114] Example 2

[0115] The key to the model proposed in this invention's ability to describe the perturbation-accelerated creep process lies in the "irreversibility" of the model parameter reduction. For intact rocks, this "irreversibility" is manifested in the process of internal microcrack initiation, development, and aggregation to form macroscopic fracture surfaces, which is closely related to the continuous development of damage. Similar patterns exist for rock structural surfaces. Therefore, the rock creep perturbation hysteresis instability combined element model based on irreversible deformation established in this invention is also applicable to rock structural surfaces.

[0116] The model proposed in this invention was based on impact creep test data of sandstone structural surfaces from the literature “Niu et al. Shear creep deformation of rock fracture distrubed by dynamic loading[J]. International Journal of Rock Mechanics and Mining Sciences, 2024, 183: 105943.”. The model parameters obtained by nonlinear least square fitting of the relevant test data are shown in Table 2.

[0117] Table 2. Parameters of the combined element model based on creep disturbance hysteresis instability of irreversibly deformable rock structural surfaces.

[0118]

[0119] like Figure 10As shown, the model parameters exhibit high accuracy through successful verification of creep impact tests on rock structural surfaces under different impact energy conditions. This model demonstrates strong application potential in engineering scenarios involving creep response, specifically covering key analytical needs such as creep remaining life prediction and disturbance parameter inversion. It addresses the time threshold failure problem of traditional models, achieving accurate assessment of the long-term stability of surrounding rock creep. This model decouples creep damage from disturbance damage through deformation decomposition, overcoming the time-dependent damage coupling of traditional models, and enabling disturbance source identification and instability process reproduction.

[0120] It should be noted that the above-described embodiments are illustrative of the technical solutions of the present invention and not limiting thereof. Equivalent substitutions or other modifications made by those skilled in the art based on the prior art, as long as they do not exceed the concept and scope of the technical solutions of the present invention, should be included within the scope of the claims of the present invention.

Claims

1. A method for predicting rock creep disturbance-induced instability and post-disturbance residual life, characterized in that, The method comprises the following steps: S1: Based on the irreversible deformation evolution law of rock, the plastic element parameters are regarded as variable parameters related to the irreversible deformation of rock creep, and the rock creep damage function D is constructed c ; S2: Based on the creep deformation composition of the rock, a Burgers model containing elastic elements, viscoelastic elements and plastic elements is constructed, and a rock creep damage function D c Substitute the Burgers model to establish a rock creep equation ε c ; S3: The element parameter is regarded as a variable parameter related to irreversible deformation of rock disturbance, and a rock disturbance damage function D is constructed d ; S4: The rock disturbance damage function D d The disturbance damage equation describing the evolution of the irreversible deformation of the rock during the disturbance instability process is established by substituting the Burgers model S5: Based on the decomposition process of rock creep instability deformation under disturbance, the disturbance damage equation The element describing the initial disturbance irreversible deformation in the equation ε c The damage plastic element describing the creep irreversible deformation, the disturbance damage equation coupling the creep load irreversible deformation evolution of rock is established S6: rock creep equation ε c The damage plastic element describing the irreversible deformation of creep is replaced by a perturbed damage equation coupled with the evolution of irreversible deformation of rock under creep load The creep-perturbation coupled nonlinear plastic body is characterized, and the rock creep-perturbation hysteresis instability combined element model equation ε is established; S7: According to the rock creep disturbance lag instability combined element model equation ε, combined with the rock instability irreversible deformation criterion to determine the instability type and state of the rock, and predict the post-disturbance instability time t of the rock f ; S8: Based on the rock's post-disturbance instability time t f In conjunction with creep lifetime, the remaining lifetime of the rock after disturbance is calculated.

2. The method of claim 1, wherein, The step S1 is specifically: S11: Based on the deformation composition of Burgers model and the rock creep instability deformation threshold, the allowable irreversible deformation of the rock ε is obtained a is: where ε a is the allowable irreversible deformation of the rock, is the total creep deformation of the rock at complete failure, ε e is the instantaneous elastic deformation of the rock, ε ve is the viscoelastic deformation of the rock; S12: Defining rock creep damage function D c is an exponential function of the irreversible deformation due to creep, expressed as: Ignoring the instantaneous irreversible deformation produced when the creep load is applied, when D c = 0; while when D c = 1.

3. The method of claim 2, wherein, In step S2, the rock creep equation ε of the element parameter is iteratively reduced with the increase of the irreversible deformation of creep c The construction method of the rock creep equation ε includes the following steps: S21: The viscous coefficient η1 of the plastic element in the Burgers model is reduced due to the continuous effect of damage, and is expressed as: η1(D c ) = (1 - D c ) η1; where η1(D c ) is the viscosity coefficient considering the reduction of creep damage; S22: Replace the viscous coefficient η1 of the plastic element in the Burgers model with the viscous coefficient η1(D c ) considering the creep damage reduction, to obtain the expression of the rock creep equation ε c whose plastic element parameters are iteratively reduced with the increase of the creep irreversible deformation.

4. The method of claim 3, wherein, The method of constructing the rock disturbance damage function D in step S3 d includes the following steps: S31: define the disturbed irreversible deformation ε of the rock at complete failure b The expression is: wherein is the total disturbed irreversible deformation at complete failure of the rock, is the initial disturbed irreversible deformation, is the accumulated irreversible deformation increment at the disturbed instability. S32: define rock disturbance damage function D d is an exponential function of the disturbance irreversible deformation as a variable: where β is the perturbation damage factor, is the total irreversible deformation during the perturbation process, which is expressed as wherein is the initial disturbance irreversible deformation, is the accumulated irreversible deformation increment during the disturbance, and D d ≠ 0; if during the disturbance D d = 1.

5. The method of claim 4, wherein, In step S4, the disturbance damage equation of the disturbance irreversible deformation evolution development in the rock disturbance instability process is constructed The construction method of the disturbance damage equation includes the following steps: S41: perturbation damage equation The perturbation irreversible deformation of the rock is further decomposed, and each model element is used to describe irreversible deformation. In constructing the perturbation damage equation, except that the initial perturbation irreversible deformation described by the initial perturbation instantaneous modulus E3 does not participate in nonlinear growth, the remaining model parameters are all reduced, and the expression is as follows: wherein η3(D d ) and η4(D d ) are the viscosity coefficients taking into account the damage reduction of the perturbation; E4(D d ) is the viscoelastic modulus taking into account the damage reduction of the perturbation; S42: replace the viscous coefficient η3 of the plastic element, the viscous coefficient η4 of the viscoelastic element and the viscoelastic modulus E4 in the Burgers model with the viscous coefficient η3(D d ), the viscous coefficient η4(D d ) and the viscoelastic modulus E4(D d ) considering the disturbance damage reduction, respectively, to obtain the expression of the disturbance damage equation of the disturbance irreversible deformation evolution development in the disturbance instability process of the rock .

6. The method of claim 5, wherein, In step S5, the perturbed damage equation of irreversible deformation evolution of rock coupled with creep load The construction method of the equation of state includes the following steps: S51: The irreversible deformation at the junction of different stress states has continuity, and the irreversible deformation of the creep finally generated before the disturbance process is equal to the initial irreversible deformation of the disturbance process, and is expressed as: wherein creep irreversible deformation at the start of the disturbance, disturbance irreversible deformation at the start of the disturbance; Irreversible deformation of rock under coupling of creep load and disturbance load The expression is composed of two parts: After the disturbance ends, the rock continues to creep, and the total irreversible deformation finally generated in the disturbance process is equal to the initial irreversible deformation of the subsequent creep process, and is expressed as: In the formula, is the total irreversible deformation of rock under the coupling action of creep load and disturbance load at the end of disturbance time, is the initial creep irreversible deformation at the end of disturbance time; S52: Rock creep damage function D under coupling action of creep load and disturbance load c , disturbance damage function D d Total irreversible deformation of rock under coupling action of creep load and disturbance load Control, the expression is: S53: Expanding the cumulative irreversible deformation in the disturbance process into a continuous function in the time domain t can couple the general creep damage and the disturbance damage in the same model, and the expression of the disturbance number n is: In the formula, △t is the time of disturbance load, and its value is (t-t x ) x The initial time of disturbance t, T is the disturbance period, and △t / T is defined as the generalized disturbance number, and the irreversible deformation of the disturbance continuously develops in the time domain t. S54: According to the continuity condition of irreversible deformation, the perturbed damage equation of the irreversible deformation evolution of the rock coupled with the creep load The expression is:

7. The method of claim 6, wherein, In step S6, according to different rock load states, the disturbance damage equation of the irreversible deformation evolution development of the rock coupled with the creep load The three modes are divided, and the construction method includes the following steps: S61: Mode 1: Rock creep disturbance lag instability combined element model equation ε expression before the action of disturbance load and rock creep equation ε in step S2.2 c The expressions are the same: S62: Mode 2: replace the damage plastic element described in step S22 with the perturbed damage equation of the irreversible deformation evolution of the rock coupled with the creep load described in step S54 The creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic body is characterized, and the creep-perturbation coupled nonlinear plastic​ S63: After the disturbance ends, the rock continues to creep, and the expression of the cumulative irreversible deformation generated by the coupling effect of the creep load and the disturbance load is: S64: Mode 3: Disturbance damage equation of the irreversible deformation evolution development of the rock coupled with the creep load in step S62 Cumulative disturbance irreversible deformation generated by the coupling of the creep load and the disturbance load in step S63 The rock creep disturbance hysteresis instability combined element model equation ε after the disturbance load is obtained, and the expression is:

Images