Construction method and system of rock damage constitutive model and simulation analysis method and system of rock damage

Abstract

The invention discloses a construction method and system of a rock damage constitutive model and a simulation analysis method and system of rock damage. The construction method comprises the following steps: soaking a target rock mass sample in a target liquid environment for a cyclic load test; according to continuous damage mechanics, a corresponding constitutive relation is established, and the constitutive relation is the relation of stress data, elastic modulus data, strain data and coupling damage; establishing a damage variable relationship, wherein the damage variable relationship is used for indicating a relationship between fatigue damage caused based on cyclic loading and mesoscopic damage caused based on target liquid environment corrosion and corresponding coupling damage; and taking irreversible plastic strain as an internal variable, and establishing a fatigue damage constitutive model of the target liquid environment corroded rock mass under the cyclic load effect based on the damage variable relationship and the constitutive relationship. According to the invention, the seawater corrosion-fatigue coupling damage variable is introduced, so that the deterioration characteristics of the fatigue mechanical property of the rock under the seawater corrosion and cyclic load coupling effect can be explored.

Description

Technical Field

[0001] The present invention relates to the field of data processing, and in particular to a method and system for constructing a constitutive model of seawater corrosion rock damage under cyclic load, and a method and system for simulating and analyzing seawater corrosion rock damage under cyclic load. Background Art

[0002] In the current construction of offshore wind power projects, the supporting structure of offshore wind power must be embedded in the seabed rock foundation. Since the rock foundation itself is corroded by sulfate ions and chloride ions in seawater, its mechanical properties deteriorate. The cyclic load generated by the dynamic disturbance of drilling construction will inevitably have an adverse effect on the dynamic characteristics and deformation damage of the structure. Therefore, the study of the dynamic characteristics and constitutive model of seawater corroded rocks under cyclic dynamic disturbance is undoubtedly of great theoretical value and engineering significance for optimizing the safe construction of engineering structures and the long-term stability evaluation. Summary of the invention

[0003] In view of the shortcoming of the lack of research on seawater corrosion under cyclic loading in the prior art, the present invention provides a method and system for constructing a seawater corrosion rock damage constitutive model under cyclic loading, and also discloses a simulation analysis method and system for seawater corrosion rock damage under cyclic loading based on the damage constitutive model.

[0004] In order to solve the above technical problems, the present invention is solved by the following technical solutions:

[0005] A method for constructing a rock damage constitutive model comprises the following steps:

[0006] The target rock mass specimen is immersed in the target liquid environment to carry out a cyclic load test;

[0007] According to continuum damage mechanics, a corresponding constitutive relationship is established, wherein the constitutive relationship is a relationship between stress data, elastic modulus data, strain data and coupled damage;

[0008] Establishing a damage variable relationship, wherein the damage variable relationship is used to indicate the relationship between fatigue damage caused by cyclic loading and mesoscopic damage caused by corrosion in a target liquid environment and corresponding coupled damage;

[0009] Taking the irreversible plastic strain as an internal variable, a fatigue damage constitutive model of the target liquid environment corroded rock mass under cyclic load is established based on the damage variable relationship and the constitutive relationship.

[0010] As an implementable embodiment, the damage variable relationship is:

[0011]

[0012] in:

[0013] D m Indicates coupling damage;

[0014] D c Indicates microscopic damage;

[0015] D f Indicates fatigue damage.

[0016] As an implementable method, the irreversible plastic strain is used as an internal variable, and the elastic modulus corresponding to each cycle process is calculated as follows:

[0017]

[0018] in:

[0019] E i represents the elastic modulus corresponding to the i-th cycle;

[0020] σ im represents the maximum cyclic stress corresponding to the i-th cycle;

[0021] ε ir represents the cyclic load generated by the construction dynamic disturbance corresponding to the i-th cycle;

[0022] P, Q, and R are model parameters obtained by fitting and calculating the elastic modulus, maximum cyclic stress, and irreversible plastic deformation corresponding to each cycle.

[0023] As an implementable method, a complete rock fatigue damage constitutive model containing internal variables is constructed to obtain a first constitutive model, wherein the first constitutive model is:

[0024]

[0025] in:

[0026] σ i represents the uniaxial stress corresponding to the i-th cycle;

[0027] E i represents the elastic modulus corresponding to the i-th cycle;

[0028] ε i represents the uniaxial strain corresponding to the i-th cycle;

[0029] D im represents the coupling damage corresponding to the i-th cycle;

[0030] σ im represents the maximum cyclic stress corresponding to the i-th cycle;

[0031] ε irrepresents the irreversible plastic deformation corresponding to the i-th cycle;

[0032] D c represents the corresponding mesoscopic damage;

[0033] D if represents the fatigue damage corresponding to the i-th cycle;

[0034] σ i represents the uniaxial stress corresponding to the i-th cycle;

[0035] P, Q, and R are model parameters.

[0036] As an implementable method, the first constitutive model is modified based on a compaction coefficient to obtain a second constitutive model, wherein the compaction coefficient is a ratio of a secant modulus to a Young's modulus;

[0037] The second constitutive model is:

[0038]

[0039] Where Cj is the compression coefficient.

[0040] A system for constructing a rock damage constitutive model, used to execute any one of the construction methods described above.

[0041] A simulation analysis method for rock damage, based on a rock damage constitutive model constructed by any one of the construction methods described above, simulates the condition of a target rock in a target liquid environment when subjected to a cyclic load caused by construction dynamic disturbance, wherein the specific steps of obtaining the damage condition and stress-strain state of the target rock at a target time are as follows:

[0042] Obtaining corresponding mesoscopic damage based on a preset corrosion damage relationship, wherein the corrosion damage relationship is a relationship between the mesoscopic damage and the pH value and immersion time of the target liquid environment;

[0043] Determining fatigue damage corresponding to the cyclic load;

[0044] Corresponding coupling damage is obtained based on the microscopic damage and the fatigue damage calculation

[0045] Based on the stress data corresponding to the target time, the mesoscopic damage and the fatigue damage, the strain condition of the target rock is determined based on the rock damage constitutive model.

[0046] As an implementable method:

[0047] Pre-immerse the rock samples based on the target rock in solutions of several pH values, record the mesoscopic losses corresponding to different immersion times, and establish the corrosion damage relationship;

[0048] The method for calculating the mesoscopic damage is to calculate and obtain the current mesoscopic damage of the rock sample based on the current porosity and initial porosity of the corresponding rock sample.

[0049] As an implementable method:

[0050] The irreversible strain method is used to calculate the fatigue damage caused by cyclic loading. The calculation formula is:

[0051]

[0052] in:

[0053] represents the fatigue damage variable of the i-th cycle;

[0054] represents the irreversible strain of the jth cycle;

[0055] N i represents the i-th cycle;

[0056] N f Indicates the total number of cycles.

[0057] A rock damage simulation analysis system is used to execute any one of the above-mentioned simulation analysis methods.

[0058] The present invention has significant technical effects due to the adoption of the above technical solution:

[0059] The present invention can explore the degradation characteristics of the fatigue mechanical properties of rocks under the coupling effects of seawater corrosion and cyclic loads by introducing seawater corrosion-fatigue coupling damage variables. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0061] Figure 1 It is a comparative schematic diagram of cyclic stress-strain curves of rocks immersed for 60 days at pH=13 during cyclic loading and unloading at 60% fatigue life in Example 1;

[0062] Figure 2 This is the test curve of fatigue stress-strain of the rock in Example 2 at pH = 13 and soaked for 60 days for the first 15 cycles;

[0063] Figure 3 This is the theoretical fatigue stress-strain curve of the rock in Example 2 at pH = 13 and immersed for 60 days for the first 15 cycles. DETAILED DESCRIPTION

[0064] The present invention is further described in detail below in conjunction with embodiments. The following embodiments are for explanation of the present invention but the present invention is not limited to the following embodiments.

[0065] Example 1: A method for constructing a rock damage constitutive model, such as Figure 1 As shown, the following steps are included:

[0066] S100. According to continuum damage mechanics, the corresponding constitutive relationship is established;

[0067] The constitutive relationship is the relationship between stress data, elastic modulus data, strain data and coupled damage;

[0068] According to the Lemaitre strain equivalence principle, any strain constitutive equation established for damaged materials can be derived in the same way as for undamaged materials by replacing the corresponding nominal stress with the effective stress represented by the damage variable.

[0069] The constitutive relation of the damaged material is expressed as follows:

[0070] [σ * ]=[σ] / (1-D)=[E][ε] / (1-D)

[0071] Where: [E] is the elastic modulus matrix of rock material, [σ*] is the effective stress matrix, [σ] is the nominal stress matrix, [ε] is the strain matrix, and D is the coupled damage.

[0072] S200, establishing damage variable relationship;

[0073] The damage variable relationship is used to indicate the relationship between fatigue damage caused by cyclic loading and mesoscopic damage caused by corrosion in the target liquid environment and the corresponding coupling damage;

[0074] Since seawater corrosion of rocks under cyclic loading contains both mesoscopic damage caused by seawater corrosion and fatigue damage caused by cyclic loading, the seawater corrosion-fatigue coupled damage should be considered in the constitutive model of rocks. The damage of seawater corrosion of rocks under cyclic loading is regarded as the coupling of seawater corrosion and fatigue damage under the same stress state. According to the Lemaitre strain equivalence hypothesis, it can be obtained:

[0075] ε m =ε c +ε f -ε v

[0076] in:

[0077] ε is the uniaxial strain;

[0078] m refers to rock material in its natural state;

[0079] c refers to the microscopic defective rock material caused by seawater corrosion;

[0080] f refers to the microscopic defective rock material caused by seawater corrosion;

[0081] v is undamaged rock material;

[0082] ε J is the uniaxial strain of the rock material in its natural state;

[0083] ε c It is the uniaxial strain of the mesoscopically defective rock material caused by seawater corrosion;

[0084] ε f is the uniaxial strain of the macroscopically damaged rock material caused by fatigue loading;

[0085] ε v is the uniaxial strain of the undamaged rock material.

[0086] Under the action of uniaxial stress σ1, the formula can be obtained:

[0087]

[0088] in:

[0089] σ1 is the uniaxial stress to which the rock material is subjected, and E is the elastic modulus of each rock material;

[0090] E m is the elastic modulus of the rock material in its natural state;

[0091] E c It is an elastic model of rock materials with mesoscopic defects caused by seawater corrosion;

[0092] E f It is an elastic model of macroscopically damaged rock materials caused by fatigue loading;

[0093] E v It is an elastic model for undamaged rock materials.

[0094] Assuming that the damage variable caused by microscopic damage in the loading direction is Dc, the damage variable caused by load damage is Df, and their coupled damage is D, then based on continuum damage mechanics, we can get:

[0095]

[0096] Based on the above formula, the coupled damage variables of seawater corrosion rock under cyclic loading can be derived, that is, the damage variable relationship:

[0097]

[0098] in:

[0099] D m Indicates coupling damage;

[0100] D c Indicates microscopic damage;

[0101] D f Indicates fatigue damage.

[0102] S300, taking irreversible plastic strain as an internal variable, and establishing a fatigue damage constitutive model of a target liquid environment corroded rock mass under cyclic load based on the damage variable relationship and the constitutive relationship.

[0103] S310, taking the irreversible plastic strain as an internal variable, determining the calculation relationship between the elastic modulus corresponding to each cycle process and the maximum cyclic stress and irreversible plastic deformation;

[0104] In addition to determining the damage variable, determining the deformation modulus is another prerequisite for constructing the constitutive model.

[0105] As is known to all, in a uniaxial cyclic compression test, the fatigue deformation modulus of rock materials decreases nonlinearly with the increase in the number of cycles. This embodiment adopts the internal variable theory to describe the degradation behavior of the fatigue deformation modulus. Since the fatigue damage variable defined by irreversible strain can better reflect the fatigue degradation of rock materials, the irreversible plastic strain is selected as the internal variable of this embodiment, and the following formula is proposed:

[0106] The calculation formula of the elastic modulus corresponding to the i-th cycle is:

[0107]

[0108] in:

[0109] E i represents the elastic modulus corresponding to the i-th cycle;

[0110] σ im represents the maximum cyclic stress corresponding to the i-th cycle;

[0111] ε ir Represents the cyclic load generated by the construction dynamic disturbance corresponding to the i-th cycle.

[0112] In this embodiment, the target rock sample is immersed in the target liquid environment for a cyclic load test. Since the pH value of seawater is generally between 8.0 and 8.5, which is weakly alkaline, in order to accelerate the progress of chemical corrosion, this experiment adopts a chemical solution with a pH of 13 to simulate the seawater corrosion environment. Those skilled in the art can set the pH value of the target liquid environment according to actual needs.

[0113] Based on the elastic modulus, maximum cyclic stress, and irreversible plastic deformation corresponding to each cycle, the model parameters corresponding to P, Q, and R are obtained by fitting and calculating according to the above calculation formula; in this embodiment, P, Q, and R can be determined by performing a multi-parameter linear fitting method in the Matlab compiler.

[0114] S320, constructing a complete rock fatigue damage constitutive model containing internal variables to obtain a first constitutive model;

[0115] The first constitutive model is:

[0116]

[0117] in:

[0118] σ i represents the uniaxial stress corresponding to the i-th cycle;

[0119] E i represents the elastic modulus corresponding to the i-th cycle;

[0120] ε i represents the uniaxial strain corresponding to the i-th cycle;

[0121] D im represents the coupling damage corresponding to the i-th cycle;

[0122] σ im represents the maximum cyclic stress corresponding to the i-th cycle;

[0123] ε ir represents the irreversible plastic deformation corresponding to the i-th cycle;

[0124] D c represents the corresponding mesoscopic damage;

[0125] D if represents the fatigue damage corresponding to the i-th cycle;

[0126] σ i represents the uniaxial stress corresponding to the i-th cycle;

[0127] P, Q, and R are model parameters.

[0128] Note that the uniaxial strain ε iThe initial value of is defined as the ratio of the minimum cyclic stress of each cycle to the initial tangent modulus. The loading modulus E is used in the loading and unloading processes. l and unloading modulus E u .

[0129] Further:

[0130] Considering the compaction phenomenon of rock mass during the first loading process, the compaction coefficient C is defined as J The first constitutive model mentioned above is modified, and the compaction coefficient C J Defined as the ratio of the secant modulus to Young's modulus.

[0131] The damage constitutive model of seawater corrosion of rocks under uniaxial cyclic compression was finally established, namely, the second constitutive model:

[0132]

[0133] Where Cj is the compression coefficient.

[0134] The constructed second constitutive model is verified based on the test data, that is, the corresponding uniaxial strain is calculated based on the second constitutive model based on the uniaxial stress, coupled damage, maximum cyclic stress and irreversible plastic deformation corresponding to each cycle in the test process, so as to verify the accuracy of the second constitutive model;

[0135] During the test, the current mesoscopic damage of the rock sample is calculated based on the current porosity and initial porosity of the corresponding rock sample. The calculation formula is as follows:

[0136]

[0137] in:

[0138] D c is the initial damage variable caused by seawater corrosion;

[0139] n0 is the porosity of the sandstone sample in the initial state;

[0140] n i is the porosity after solution corrosion during the i-th cycle.

[0141] During the test, the irreversible strain method was used to calculate the fatigue damage caused by cyclic loading. The calculation formula is as follows:

[0142] The calculation formula is:

[0143]

[0144] in:

[0145] represents the fatigue damage variable of the i-th cycle;

[0146] represents the irreversible strain of the jth cycle;

[0147] N i represents the i-th cycle;

[0148] N f Indicates the total number of cycles.

[0149] like Figure 1 As shown, the fatigue stress-strain theoretical curve and experimental data comparison of sandstone immersed in pH=13 solution at 60% fatigue life for 60 days;

[0150] like Figure 2 As shown, the fatigue stress-strain test curves for the first 15 cycles.

[0151] like Figure 3 Shown are the theoretical fatigue stress-strain curves for the first 15 cycles.

[0152] based on Figures 1 to 3 It can be seen that the fatigue stress-strain damage constitutive model of sandstone corroded by alkaline solution is in good agreement with the experimental results in both single cycle and overall conditions. Therefore, the seawater corrosion-fatigue constitutive model established in this study with the irreversible plastic deformation of intact rock as the internal variable can accurately describe the hysteresis effect and cumulative effect of the deformation of the submarine rock foundation of offshore wind power under cyclic loads, which has important theoretical value and engineering significance for optimizing the safe construction of offshore wind power engineering structures and long-term stability evaluation.

[0153] Example 2: A system for constructing a rock damage constitutive model, used to execute the construction method described in Example 1.

[0154] Example 3: A simulation analysis method for rock damage, based on the rock damage constitutive model constructed by the construction method described in Example 1, simulates the condition of the target rock in a target liquid environment when subjected to cyclic loads caused by construction dynamic disturbances, wherein the specific steps for obtaining the damage condition and stress-strain state of the target rock at the target time are:

[0155] S100, obtaining corresponding mesoscopic damage based on a preset corrosion damage relationship, wherein the corrosion damage relationship is a relationship between the mesoscopic damage and the pH value and immersion time of the target liquid environment;

[0156] In this embodiment, a rock sample based on the target rock is pre-immersed in solutions of several pH values, and the mesoscopic losses corresponding to different immersion times are recorded to establish a corrosion damage relationship;

[0157] The method for calculating the mesoscopic damage is to calculate and obtain the current mesoscopic damage of the rock sample based on the current porosity and initial porosity of the corresponding rock sample.

[0158] In the experiment, the rock was chemically corroded under different alkaline conditions. By calculating the porosity of the corroded sample, the fine damage variable (D c ) with corrosion time (y) and pH value (x), R2 = 99766, the fitting effect is good, as shown below:

[0159] D c =-0.05569-5.17737×10 -5 x 2 +0.00497x-4.34089×10 -6 y 2 +7.69103×10 -4 y

[0160] In actual engineering, the measured seawater pH value and seawater corrosion time are substituted into the above formula to obtain the mesoscopic damage variable under this condition, so as to predict the seawater corrosion damage in actual engineering.

[0161] S200, determining fatigue damage corresponding to the cyclic load;

[0162] The irreversible strain method is used to calculate the fatigue damage caused by cyclic loading.

[0163] S300, obtaining corresponding coupling damage based on the microscopic damage and the fatigue damage calculation

[0164] S400: Based on the stress data corresponding to the target time, the mesoscopic damage and the fatigue damage, the strain condition of the target rock is determined based on the rock damage constitutive model.

[0165] Embodiment 4: A rock damage simulation analysis system, used for the simulation analysis method described in Embodiment 3.

[0166] As for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0167] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0168] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, devices, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0169] The present invention is described with reference to the flowcharts and / or block diagrams of the method, terminal device (system), and computer program product according to the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing terminal device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0170] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing terminal device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0171] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal device so that a series of operating steps are executed on the computer or other programmable terminal device to produce a computer-implemented process, thereby providing instructions for executing on the computer or other programmable terminal device to implement the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0172] It should be noted that:

[0173] The "one embodiment" or "embodiment" mentioned in the specification means that the specific features, structures or characteristics described in conjunction with the embodiment are included in at least one embodiment of the present invention. Therefore, the phrases "one embodiment" or "embodiment" appearing in various places throughout the specification do not necessarily refer to the same embodiment.

[0174] Although the preferred embodiments of the present invention have been described, those skilled in the art may make additional changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the present invention.

[0175] In addition, it should be noted that the shapes and names of the parts and components of the specific embodiments described in this specification may be different. Any equivalent or simple changes made based on the structure, features and principles described in the patent concept of the present invention are included in the protection scope of the patent of the present invention. The technicians in the technical field of the present invention can make various modifications or supplements to the specific embodiments described or replace them in a similar manner, as long as they do not deviate from the structure of the present invention or exceed the scope defined by the claims, they should all fall within the protection scope of the present invention.

Claims

1. A method for constructing a rock damage constitutive model, characterized in that The following steps are involved: The target rock mass specimen is immersed in the target liquid environment to carry out a cyclic load test; According to continuum damage mechanics, a corresponding constitutive relationship is established, wherein the constitutive relationship is a relationship between stress data, elastic modulus data, strain data and coupled damage; Establishing a damage variable relationship, wherein the damage variable relationship is used to indicate the relationship between fatigue damage caused by cyclic loading and mesoscopic damage caused by corrosion in a target liquid environment and corresponding coupled damage; Taking the irreversible plastic strain as an internal variable, a fatigue damage constitutive model of a target liquid environment corroded rock mass under cyclic load is established based on the damage variable relationship and the constitutive relationship.

2. The method for constructing a rock damage constitutive model according to claim 1, characterized in that: The damage variable relationship is: in: D m Indicates coupling damage; D c Indicates microscopic damage; D f Indicates fatigue damage.

3. The method for constructing a rock damage constitutive model according to claim 2, characterized in that: Taking the irreversible plastic strain as the internal variable, the calculation formula of the elastic modulus corresponding to each cycle process is: in: E i represents the elastic modulus corresponding to the i-th cycle; σ im represents the maximum cyclic stress corresponding to the i-th cycle; ε ir represents the cyclic load generated by the construction dynamic disturbance corresponding to the i-th cycle; P, Q, and R are model parameters obtained by fitting and calculating the elastic modulus, maximum cyclic stress, and irreversible plastic deformation corresponding to each cycle.

4. The method for constructing a rock damage constitutive model according to claim 3, characterized in that: A complete rock fatigue damage constitutive model containing internal variables is constructed to obtain the first constitutive model, which is: in: σ i represents the uniaxial stress corresponding to the i-th cycle; E i represents the elastic modulus corresponding to the i-th cycle; ε i represents the uniaxial strain corresponding to the i-th cycle; D im represents the coupling damage corresponding to the i-th cycle; σ im represents the maximum cyclic stress corresponding to the i-th cycle; ε ir represents the irreversible plastic deformation corresponding to the i-th cycle; D c Indicates the corresponding mesoscopic damage; D if represents the fatigue damage corresponding to the i-th cycle; σ i represents the uniaxial stress corresponding to the i-th cycle; P, Q, and R are model parameters.

5. The method for constructing a rock damage constitutive model according to claim 4, characterized in that: The first constitutive model is modified based on the compression coefficient to obtain a second constitutive model, wherein the compression coefficient is the ratio of the secant modulus to the Young's modulus; The second constitutive model is: Where Cj is the compression coefficient.

6. A system for constructing a rock damage constitutive model, characterized in that: Used to execute the construction method described in any one of claims 1 to 5.

7. A simulation analysis method for rock damage, characterized in that: The rock damage constitutive model constructed based on the construction method according to any one of claims 1 to 5 is used to simulate the condition of the target rock under cyclic load caused by construction dynamic disturbance in the target liquid environment, wherein the specific steps of obtaining the damage condition and stress-strain state of the target rock at the target time are as follows: Obtaining corresponding mesoscopic damage based on a preset corrosion damage relationship, wherein the corrosion damage relationship is a relationship between the mesoscopic damage and the pH value and immersion time of the target liquid environment; Determining fatigue damage corresponding to the cyclic load; Corresponding coupling damage is obtained based on the microscopic damage and the fatigue damage calculation Based on the stress data corresponding to the target time, the mesoscopic damage and the fatigue damage, the strain condition of the target rock is determined based on the rock damage constitutive model.

8. A rock damage simulation analysis method according to claim 7, characterized in that: Pre-immerse the rock samples of the target rock in solutions of several pH values, record the mesoscopic losses corresponding to different immersion times, and establish the corrosion damage relationship; The method for calculating the mesoscopic damage is to calculate and obtain the current mesoscopic damage of the rock sample based on the current porosity and initial porosity of the corresponding rock sample.

9. A rock damage simulation analysis method according to claim 7, characterized in that: The irreversible strain method is used to calculate the fatigue damage caused by cyclic loading. The calculation formula is: in: represents the fatigue damage variable of the i-th cycle; represents the irreversible strain of the jth cycle; N i represents the i-th cycle; N f Indicates the total number of cycles.

10. A rock damage simulation analysis system, characterized in that: Used to execute the simulation analysis method described in any one of claims 7 to 9.

Images