Stress calculation method and device for tension-compression dual-modulus rock and electronic equipment

Abstract

The invention relates to the technical field of rock mechanics and engineering, and discloses a stress calculation method and device for tension-compression dual-modulus rock and electronic equipment. The stress calculation method for the tension-compression dual-modulus rock comprises the following steps: performing parameter sensitivity analysis and data fitting on key parameters of the rock based on a numerical model of a rock tension-compression dual-modulus hypothesis, and obtaining a parameter sensitivity analysis result and a data fitting result; and constructing a mathematical model of rock stress based on the parameter sensitivity analysis result and the data fitting result. By adopting the method, the stress distribution of the center point of the Brazilian disc can be determined, and the problem that the stress of the center point of the Brazilian disc of the tension-compression dual-modulus rock cannot be accurately determined by an existing calculation method is solved.

Description

Technical Field

[0001] The present invention relates to the field of rock mechanics and engineering technology, and in particular to a stress calculation method, a calculation device and an electronic device for tensile and compressive dual modulus rock. Background Art

[0002] The tensile strength of rock plays an important role in geotechnical engineering, geological engineering, mining, tunnel construction, slope stability, water conservancy and hydropower engineering, and underground space development. As a key parameter to measure the ability of rock to resist tensile damage, it is directly related to the stability and safety of engineering structures. The tensile strength of rock is not only the basis for determining design parameters and optimizing the design of support structures, but also an important basis for formulating geological disaster prevention and control strategies. Therefore, accurate assessment of the tensile strength of rock is crucial to effectively prevent potential risks such as collapse and landslides, and to ensure engineering safety and long-term stable operation.

[0003] The Brazilian disk test, as a classic method for measuring the tensile strength of rocks, has been widely respected since its inception. The experiment is based on the isotropic elastic assumption, which simplifies the difficulty of analysis, but it implicitly assumes that the elastic properties of rocks under tension and compression are the same. However, this is not the case in reality. Most rocks exhibit different elastic responses under tension and compression, which is called the tensile-compression dual modulus characteristic. Specifically, the slope of the stress-strain curve of the rock during compression (i.e., the compression modulus) E c ) is often greater than its slope during stretching (i.e., the tensile modulus E t ). Different rock types E c / E t The ratio varies, with dense rocks close to 1:1, some sandstones at 2:1, and weaker rocks as high as 10:1. Therefore, it is necessary to further study more precise measurement methods to more accurately calculate the tensile strength of rocks with both tension and compression moduli. In order to accurately measure the tensile strength of rocks, it is first necessary to obtain the stress distribution state at the center of the disk in the Brazilian disk test. This stress state does not exist in isolation, but is closely dependent on the inherent elastic modulus characteristics of the rock material, and the two form a mutually dependent relationship. Therefore, in the absence of sufficient prior knowledge of the elastic modulus of rock, it is challenging to directly derive an accurate analytical solution for the stress distribution at the center of the Brazilian disk. Summary of the invention

[0004] The present invention provides a stress calculation method, a calculation device and an electronic device for tensile and compressive dual modulus rock, which are used to solve the problem in the prior art that the stress distribution of the center point of a Brazilian disk cannot be accurately derived.

[0005] The present invention provides a stress calculation method for tensile and compressive dual modulus rock, comprising: Based on the numerical model of rock tension and compression dual modulus assumption, parameter sensitivity analysis and data fitting are performed on the key parameters of rock to obtain parameter sensitivity analysis results and data fitting results; Based on the parameter sensitivity analysis results and the data fitting results, a mathematical model of rock stress is constructed.

[0006] According to the stress calculation method for tensile and compressive dual modulus rock provided by the present invention, in the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis result and the data fitting result, the method specifically includes: Based on the parameter sensitivity analysis result, constructing a first function of rock stress and the key parameter; Based on the data fitting result, obtaining a fitting function of the key parameters; Based on the first function and the fitting function, a mathematical model of the rock stress is constructed.

[0007] According to the stress calculation method for tensile and compressive dual modulus rock provided by the present invention, the step of constructing a first function of rock stress and the key parameter based on the parameter sensitivity analysis result includes: Based on the parameter sensitivity analysis results, construct a second function of a plurality of the key parameters; constructing a third function of rock stress and one of the key parameters; Based on the second function and the third function, a first function of rock stress and the key parameter is constructed.

[0008] According to the stress calculation method of the tensile and compressive dual modulus rock provided by the present invention, the key parameters include the Poisson's ratio of the rock under compression state. n c , rock strain, rock compression modulus E c and the tensile modulus of rock E t .

[0009] According to the stress calculation method for tensile and compressive dual modulus rock provided by the present invention, in the step of constructing the second function of the plurality of key parameters based on the parameter sensitivity analysis result, the method specifically includes: Based on the results of parameter sensitivity analysis, the parameters k and the strain and Poisson's ratio of the rock under compression n c The second function, where the parameter k The compression modulus E cWith the tensile modulus E t The absolute value of the ratio.

[0010] According to the stress calculation method for tensile and compressive dual modulus rock provided by the present invention, in the step of performing parameter sensitivity analysis and data fitting on key parameters of the rock and obtaining parameter sensitivity analysis results and data fitting results, the method for performing the data fitting includes: A bivariate cubic polynomial function is used to fit the key parameters.

[0011] According to the stress calculation method for tensile and compressive dual modulus rock provided by the present invention, after the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis result and the data fitting result, the method further includes: Based on the horizontal stress function in the mathematical model, the tensile strength of the rock is determined.

[0012] According to the stress calculation method for tensile and compressive dual modulus rock provided by the present invention, the constitutive equation of the numerical model in the principal stress coordinate system is: Formula (1)

[0013] Formula (1) is the stress tensor, is the strain tensor, is the fourth-order elastic tensor in the principal stress coordinate system; Formula (2); Formula (3); In formula (2), E t represents the tensile modulus of rock, E c They represent the compression modulus of rock, n t represents the Poisson's ratio associated with the stretched state, n c represents the Poisson's ratio associated with the compressive state.

[0014] A second aspect of the present invention provides a stress calculation device for tensile and compressive dual modulus rocks, comprising: The first acquisition module is used to perform parameter sensitivity analysis and data fitting on key parameters of rock based on the numerical model of rock tension and compression dual modulus assumption, and obtain parameter sensitivity analysis results and data fitting results; The model building module is used to build a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results.

[0015] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the stress calculation method for tensile and compressive dual modulus rock as described in any one of the above is implemented.

[0016] The present invention provides a stress calculation method, calculation device and electronic device for tension-compression dual modulus rock. Based on the numerical model, parameter sensitivity analysis and data fitting are performed on the key parameters of the rock; then, based on the parameter sensitivity analysis results and data fitting results, a mathematical model of rock stress is constructed. Compared with the stress distribution directly obtained by using the existing calculation formula, the mathematical model of stress of the present invention is more in line with the actual stress condition of the rock, more accurate and true, and solves the problem in the prior art that the stress distribution of the center point of the Brazilian disk of tension-compression dual modulus rock cannot be accurately derived. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are 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.

[0018] Figure 1 It is one of the flow diagrams of the stress calculation method for tensile and compressive dual modulus rock provided by the present invention.

[0019] Figure 2 This is the second flow chart of the stress calculation method for tensile and compressive dual modulus rock provided by the present invention.

[0020] Figure 3 This is a schematic diagram of the existing Brazilian disk experiment.

[0021] Figure 4 is a stress and strain sensitivity analysis diagram of the center point of the Brazilian disk of Example 1 of the present invention; wherein (a) E c and k = E c / E t The impact of n c =0.35); (b) k = E c / E t and n c The impact of E c =10GPa).

[0022] Figure 5 is the stress and e y / e x and ν c relationship diagram.

[0023] Figure 6 is a stress and strain sensitivity analysis diagram of the center point of the Brazilian disk of Example 2 of the present invention; wherein (a) E c and k = E c / E t The impact of n c =0.35); (b) k = E c / E t and n c The impact of E c =10GPa).

[0024] Figure 7 is the stress and e y / e x and ν c relationship diagram.

[0025] Figure 8 It is a structural schematic diagram of stress calculation of tensile and compressive dual modulus rock provided by the present invention.

[0026] Figure 9 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0027] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0028] In the embodiments of the present invention, "at least one" refers to one or more, and "more than one" refers to two or more. "And / or" describes the association relationship of associated objects, indicating that three relationships may exist. For example, A and / or B may represent: A exists alone, A and B exist at the same time, and B exists alone, where A and B may be singular or plural. In the textual description of the present invention, the character " / " generally indicates that the previous and next associated objects are in an "or" relationship. In addition, it should be noted that the serial numbers themselves, such as "first", "second", etc., used to distinguish the objects described in the present invention are only used to distinguish the objects described, and do not have any order or technical meaning.

[0029] In the Brazilian disk test, a thin circular rock disk is placed between two rigid platens and compressed so that the applied load acts along the diameter of the specimen, e.g. Figure 3 This test method induces a tensile stress state in the center area of ​​the disc specimen, which eventually causes the specimen to break and fail due to tension. It is generally believed that the horizontal force state and vertical stress state at the center point of a uniform isotropic Brazilian disc are: Formula (6-1); Formula (6-2).

[0030] In formula (6-1) and formula (6-2), P is the external load, D is the disk diameter, and L is the disk thickness. When P is the maximum load when the specimen fails, is the tensile strength of rock.

[0031] Under the condition of keeping the applied load constant at 10kN, the strain value at the center of the disk is currently determined by the existing numerical simulation method. and And the horizontal numerical stress at the center of the disk and vertical numerical stress Taking five different cases as examples, all numerical simulation results and the calculation results of formula (6-1) and formula (6-2) are shown in Table 1. In addition, in order to quantitatively evaluate their differences, the relative errors of the numerical stress and the calculation results of formula (6-1) and formula (6-2) are defined as follows: Formula (7-1).

[0032] Formula (7-2).

[0033] Table 1 Calculation results of numerical stress and theoretical stress

[0034] From Table 1, we can find thatE c / E t =1, the error between the two is almost 0. E c / E t increase, The maximum error is as high as 54.95%, and the minimum is 18.81%. The error is mainly due to the fact that formula (6-1) is based on the assumption of a single tensile and compressive modulus (i.e. E t = / E c and n t = n c ).

[0035] From the above analysis, it can be seen that the current theoretical stress obtained by formula (6-1) cannot accurately reflect the stress state of the rock. Therefore, how to accurately determine the stress distribution at the center of the Brazilian disk is an urgent problem to be solved in the industry.

[0036] In order to solve the above technical problems, the present invention provides a stress calculation method, a calculation device and an electronic device for tensile and compressive dual modulus rocks, which can accurately predict the stress distribution at the center point of the Brazilian disk, thereby ensuring the accuracy of the calculation results of the tensile strength of the rock and the accuracy of the calculation results of the elastic modulus.

[0037] Combine the following Figure 1 to Figure 9 The stress calculation method, computing device and electronic device of the tensile and compressive dual modulus rock of the present invention are described in detail. Among them, the execution subject of the stress calculation method of the tensile and compressive dual modulus rock is an electronic device or a server. The electronic device can be a personal computer (PC), a portable device, a laptop, a smart phone, a tablet computer, a portable wearable device and other devices; the server can refer to a single server, or a server cluster composed of multiple servers, a cloud server, etc. The present invention does not specifically limit the specific form of the electronic device or the server. Furthermore, the stress calculation method of the tensile and compressive dual modulus rock can also be applied to a computing device set in an electronic device or a server, and the computing device can be implemented by software, hardware or a combination of the two. The following takes the execution subject of the stress calculation method of the tensile and compressive dual modulus rock as an example of a computing device to describe the stress calculation method of the tensile and compressive dual modulus rock of the present invention.

[0038] In order to facilitate understanding of the stress calculation method for tensile and compressive dual modulus rock provided by the embodiment of the present invention, the stress calculation method for tensile and compressive dual modulus rock provided by the present invention will be described in detail below through the following exemplary embodiments. It can be understood that the following exemplary embodiments can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments.

[0039] like Figure 1 As shown, the first aspect of the present invention provides a stress calculation method for a tensile-compression dual modulus rock. The stress calculation method for a tensile-compression dual modulus rock comprises: S100. Based on the numerical model of rock tension and compression dual modulus assumption, parameter sensitivity analysis and data fitting are performed on key parameters of rock to obtain parameter sensitivity analysis results and data fitting results; S200. Based on the parameter sensitivity analysis results and data fitting results, a mathematical model of rock stress is constructed.

[0040] In this embodiment, based on the numerical model, parameter sensitivity analysis and data fitting are performed on the key parameters of the rock; then, based on the parameter sensitivity analysis results and the data fitting results, a mathematical model of rock stress is constructed. Compared with the stress distribution obtained by using the traditional stress calculation formula, the mathematical model of stress in this embodiment is more in line with the actual stress conditions of the rock, more realistic, and more accurate, and solves the problem in the prior art that the stress distribution at the center of the Brazilian disk of the tensile and compressive dual modulus rock cannot be accurately determined.

[0041] It should be noted that the mathematical model of rock stress in this embodiment includes a mathematical model of horizontal stress and / or a mathematical model of vertical stress.

[0042] like Figure 2 As shown, further, the stress calculation method of the tensile and compressive dual modulus rock of this embodiment also includes: S300: Determine the rock tensile strength based on the mathematical model of horizontal stress in the mathematical model. It can be understood that when the external load of the rock is equal to the breaking load, the rock horizontal stress at this time is the rock tensile strength.

[0043] In some embodiments, the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results specifically includes: Based on the results of parameter sensitivity analysis, the first function of rock stress and key parameters is constructed; Based on the data fitting results, obtain the fitting function of key parameters; Based on the first function and the fitting function, a mathematical model of rock stress is constructed.

[0044] Exemplarily, the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results specifically includes: Based on the results of parameter sensitivity analysis, the first function of rock horizontal stress and key parameters is constructed; Based on the data fitting results, obtain the fitting function of key parameters; Based on the first function and the fitting function, a mathematical model of rock horizontal stress is constructed.

[0045] Exemplarily, the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results specifically includes: Based on the results of parameter sensitivity analysis, the first function of rock vertical stress and key parameters is constructed; Based on the data fitting results, obtain the fitting function of key parameters; Based on the first function and the fitting function, a mathematical model of rock vertical stress is constructed.

[0046] Furthermore, based on the parameter sensitivity analysis results, the step of constructing the first function of rock stress and key parameters includes: Based on the results of parameter sensitivity analysis, a second function of multiple key parameters is constructed; Construct a third function of rock stress and a key parameter; Based on the second function and the third function, the first function of rock stress and key parameters is constructed.

[0047] Preferably, the key parameters include Poisson's ratio of the rock under compression n c , rock strain, rock compression modulus E c and the tensile modulus of rock E t .

[0048] Furthermore, the step of constructing a second function of multiple key parameters based on the parameter sensitivity analysis results specifically includes: Based on the results of parameter sensitivity analysis, the parameters k Poisson's ratio under strain and rock compression n c The second function, where the parameter k is the compression modulus E c and tensile modulus E t The absolute value of the ratio.

[0049] Exemplarily, based on the parameter sensitivity analysis results, the step of constructing a first function of rock stress and key parameters includes: Based on the results of parameter sensitivity analysis, the parameters k Poisson's ratio under strain and rock compression n c The second function of Constructing rock horizontal stress and parameters k The third function of Based on the second function and the third function, the first function of rock horizontal stress and key parameters is constructed.

[0050] Exemplarily, based on the parameter sensitivity analysis results, the step of constructing a first function of rock stress and key parameters includes: Based on the results of parameter sensitivity analysis, the parameters k Poisson's ratio under strain and rock compression n c The second function of Constructing rock vertical stress and parameters k The third function of Based on the second function and the third function, the first function of rock vertical stress and key parameters is constructed.

[0051] In some embodiments, in the step of performing parameter sensitivity analysis and data fitting on key parameters of rock and obtaining parameter sensitivity analysis results and data fitting results, the method of performing data fitting includes: A bivariate cubic polynomial function was used to fit the key parameters of the data.

[0052] Specifically, a cubic polynomial function is used to calculate the strain and Poisson's ratio of rock under compression. n c Perform data fitting.

[0053] Example 1

[0054] The following is a detailed description of the stress calculation method of the tensile and compressive dual modulus rock of the present invention by taking the construction of a mathematical model of the horizontal stress of the rock as an example.

[0055] Step 1: The principle of the numerical model of the rock tension and compression dual modulus assumption in this embodiment is as follows: Based on the strain-driven constitutive model of rock tensile and compressive dual modulus assumption, the constitutive equation in the principal stress coordinate system takes the following form: Formula (1); Among them, i, j, k, l can be 1, 2 or 3 respectively.

[0056] Formula (1) is the stress tensor, is the strain tensor, is the fourth-order elastic tensor in the principal stress coordinate system; The value of the element is determined by the sign of the principal stress, specifically: Formula (2); Formula (3); In formula (2), E t represents the tensile modulus of rock, E c They represent the compression modulus of rock, n t represents the Poisson's ratio associated with the stretched state, n c represents the Poisson's ratio associated with the compressive state.

[0057] It is worth noting that according to the laws of thermodynamics, the following relationship can be derived: n t / E t = n c / E c .

[0058] Formula (1) represents the constitutive equation in the principal stress coordinate system. In order to transform this equation from the principal coordinate system Oabg to the global coordinate system Oxyz, the following relationship is applied: Formula (4); ε in formula (4) g represents the stress tensor in the global coordinate system, ε p represents the stress tensor in the principal stress coordinate system, and T represents the transformation matrix, which is defined as: Formula (5); All the above formulas are solved using the finite element method. The elastic modulus and Poisson's ratio at each point in the model are not only determined by the basic properties of the material, but also affected by the current stress state at that point. Correspondingly, the stress and strain at each point are also affected by the local elastic constants. Therefore, there is a cross-coupling effect between the elastic constants and the stress state of the material. In order to ensure the convergence and accuracy of the numerical solution, an iterative stress transfer algorithm is implemented in the calculation. The convergence criterion is set as the Euclidean norm of the global elastic matrix reaching a preset tolerance threshold.

[0059] Step 2: Conduct sensitivity analysis on key parameters of rock based on the above numerical model.

[0060] like Figure 4 The parameter sensitivity analysis results obtained by the above numerical model are given. The key parameters at this time include the Poisson's ratio under rock compression n c , the horizontal strain at the center of the rock e x , the vertical strain at the center of the rock e y , the compression modulus of rock E c and the tensile modulus of rock E t The horizontal stress and strain were normalized to make it easier to compare the effects of different parameters. k changes significantly with changes in E c and n c From a mathematical point of view, the horizontal stress can be expressed as a function of the parameter k The relevant functions are as follows: Formula (8).

[0061] also, e y / e x With parameters k and n c changes with changes in E c Therefore, e y / e x Can be expressed as a parameter k and n c Function: Formula (9).

[0062] It should be noted that for the parameters k and n c For every unique combination of e y / e x The value corresponds to it, and vice versa. That is, e y / e x and n cEach unique combination of uniquely determines the parameter k. It reflects a one-to-one correspondence. Therefore, the parameter k can be expressed as e y / e x and n c There is an inverse function: Formula (10).

[0063] Substituting formula (10) into formula (8), we obtain: Formula (11).

[0064] Among them, formula (11) is the function model of parameter sensitivity analysis results.

[0065] Step 3: Perform data fitting on the key parameters of rock based on the above numerical model.

[0066] Figure 5 Shows stress, e y / e x and n c 3D scatter plot between . In this figure, discrete data points are presented in the form of solid circles and are marked as "value points" in the legend. After a detailed comparison and analysis of various curve fitting methods, a bivariate cubic polynomial function was selected to fit these data, and the specific form is as follows: Formula (12); In formula (12) , , a i and b i are the fitting coefficients.

[0067] Determine the fitting coefficients through regression analysis based on discrete data a i and b i Formula (12) is the data fitting result.

[0068] Step 4: Based on the parameter sensitivity analysis results and data fitting results, a mathematical model of rock horizontal stress is constructed.

[0069] Specifically, by substituting formula (12) into formula (11), we can obtain the expression formula (13) for rock fitting horizontal stress: Formula (13).

[0070] In order to verify the accuracy of the calculation results obtained by formula (13), the fitting stress ( ) and the numerical stress obtained by numerical simulation ( ) are compared, see Table 2. It can be seen that the fitting stress ( ) in all cases is consistent with the numerical stress ( ). The accuracy of the current calculation formula (13) is significantly improved compared to the results of the previous calculation formula (6-1). In all cases, The deviation is always less than 1%. This deviation occurs due to the small difference between the fitted stress and the actual stress. These results show that the reliability of the method proposed in this embodiment is much higher than that of previous methods.

[0071] Table 2 Calculation results obtained using formula (13)

[0072] Furthermore, the stress calculation method of the tensile and compressive dual modulus rock of this embodiment may also include: Step 5: Determine the rock tensile strength based on the fitted horizontal stress formula (13). Because this embodiment can more accurately obtain the stress value of the center point of the Brazilian disk of the tensile and compressive dual modulus rock, it effectively solves the existing problem of large error in the measurement of rock tensile strength.

[0073] Example 2

[0074] The following is a detailed description of the stress calculation method for tensile and compressive dual modulus rock of the present invention by taking the simultaneous construction of a mathematical model of rock horizontal stress and a mathematical model of rock vertical stress as an example.

[0075] Step 1: The principle of the numerical model of the rock tension and compression dual modulus assumption in this embodiment is the same as that in Example 1 and will not be repeated here.

[0076] Step 2: Conduct sensitivity analysis on key parameters of rock based on the above numerical model.

[0077] like Figure 6 The parameter sensitivity analysis results obtained by the above numerical model are given. The key parameters at this time include the Poisson's ratio under rock compression n c , the horizontal strain at the center of the rock e x , the vertical strain at the center of the rock e y , the compression modulus of rock E c and the tensile modulus of rock Et . For stress (including horizontal stress and vertical stress ) and strain variables were normalized to make it easier to compare the effects of different parameters. The analysis results show that stress (including horizontal stress and vertical stress ) changes significantly with the parameter k, however, it is not affected by E c and n c From a mathematical point of view, the horizontal stress and vertical stress can be expressed as only with parameters k The relevant functions are as follows: Formula (8); Formula (14).

[0078] also, e y / e x With parameters k and n c changes with changes in E c Therefore, e y / e x Can be expressed as a parameter k and n c Function: Formula (9).

[0079] It should be noted that for the parameter k and n c For every unique combination of e y / e x The value corresponds to it, and vice versa. That is, e y / e x and n c Each unique combination of uniquely determines the parameters k .function It reflects a one-to-one correspondence. Therefore, the parameter k It can be expressed as ey / e x and n c There is an inverse function: Formula (10).

[0080] Substituting formula (10) into formula (8) and formula (14) respectively, we can obtain: Formula (15).

[0081] Among them, formula (15) is the function model of parameter sensitivity analysis results.

[0082] Step 3: Perform data fitting on the key parameters of rock based on the above numerical model.

[0083] Figure 7 Shows stress, e y / e x and n c 3D scatter plot between . In this figure, discrete data points are presented in the form of solid circles and are marked as "value points" in the legend. After a detailed comparison and analysis of various curve fitting methods, a bivariate cubic polynomial function was selected to fit these data, and the specific form is as follows: Formula (16).

[0084] In formula (16) , , a i and b i are the fitting coefficients.

[0085] Determine the fitting coefficients through regression analysis based on discrete data a i and b i Formula (16) is the data fitting result.

[0086] Step 4: Based on the results of parameter sensitivity analysis and data fitting, a mathematical model of rock horizontal stress and a mathematical model of rock vertical stress are constructed.

[0087] Specifically, formula (16) is substituted into formula (15) to obtain the horizontal fitting stress of the rock and the vertical fitting stress of the rock, as shown in the following formula (17).

[0088] Formula (17).

[0089] Table 3 gives the calculation results obtained using formula (17) and compares them with the numerical simulation results. The results show that the fitting stress errors are all less than 1%.

[0090] Table 3 Calculation results obtained using formula (17)

[0091] like Figure 8 As shown, a specific embodiment of the second aspect of the present invention provides a stress calculation device for a tensile-compressive dual modulus rock. The stress calculation device for a tensile-compressive dual modulus rock includes a first acquisition module and a model construction module; wherein the first acquisition module is used to perform parameter sensitivity analysis and data fitting on key parameters of the rock based on a numerical model of the tensile-compressive dual modulus assumption of the rock, and obtain parameter sensitivity analysis results and data fitting results. The model construction module is used to construct a mathematical model of rock stress based on the parameter sensitivity analysis results and data fitting results.

[0092] Figure 9 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 9 As shown, the electronic device may include: a processor 810, a communication interface 820, a memory 830 and a communication bus 840, wherein the processor 810, the communication interface 820 and the memory 830 communicate with each other through the communication bus 840. The processor 810 may call the logic instructions in the memory 830 to execute the stress calculation method of the tension-compression dual modulus rock, the method comprising: based on the numerical model of the tension-compression dual modulus assumption of the rock, performing parameter sensitivity analysis and data fitting on the key parameters of the rock, obtaining the parameter sensitivity analysis results and the data fitting results; based on the parameter sensitivity analysis results and the data fitting results, constructing a mathematical model of rock stress.

[0093] In addition, the logic instructions in the above-mentioned memory 830 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0094] On the other hand, the present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the stress calculation method for tensile and compressive dual modulus rock provided by the above-mentioned methods. The method includes: based on a numerical model of the tensile and compressive dual modulus assumption of rock, performing parameter sensitivity analysis and data fitting on key parameters of the rock to obtain parameter sensitivity analysis results and data fitting results; based on the parameter sensitivity analysis results and the data fitting results, constructing a mathematical model of rock stress.

[0095] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, it is implemented to execute the stress calculation method for tensile and compressive dual modulus rock provided by the above-mentioned methods. The method includes: based on a numerical model of the tensile and compressive dual modulus assumption of rock, performing parameter sensitivity analysis and data fitting on key parameters of the rock, and obtaining parameter sensitivity analysis results and data fitting results; based on the parameter sensitivity analysis results and the data fitting results, constructing a mathematical model of rock stress.

[0096] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0097] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0098] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A stress calculation method for tensile and compressive dual modulus rock, characterized in that: include: Based on the numerical model of rock tension and compression dual modulus assumption, parameter sensitivity analysis and data fitting are performed on the key parameters of rock to obtain parameter sensitivity analysis results and data fitting results; Based on the parameter sensitivity analysis results and the data fitting results, a mathematical model of rock stress is constructed.

2. The stress calculation method for tensile and compressive dual modulus rock according to claim 1 is characterized in that: The step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results specifically includes: Based on the parameter sensitivity analysis result, constructing a first function of rock stress and the key parameter; Based on the data fitting result, obtaining a fitting function of the key parameters; Based on the first function and the fitting function, a mathematical model of the rock stress is constructed.

3. The stress calculation method for tensile and compressive dual modulus rock according to claim 3 is characterized in that: The step of constructing a first function of rock stress and the key parameter based on the parameter sensitivity analysis result includes: Based on the parameter sensitivity analysis results, construct a second function of a plurality of the key parameters; constructing a third function of rock stress and one of the key parameters; Based on the second function and the third function, a first function of rock stress and the key parameter is constructed.

4. The stress calculation method for tensile and compressive dual modulus rock according to claim 3 is characterized in that: The key parameters include Poisson's ratio under rock compression ν c , rock strain, rock compression modulus E c and the tensile modulus of rock E t .

5. The stress calculation method for tensile and compressive dual modulus rock according to claim 4 is characterized in that: The step of constructing a second function of a plurality of the key parameters based on the parameter sensitivity analysis result specifically includes: Based on the results of parameter sensitivity analysis, the parameters k and the strain and Poisson's ratio of the rock under compression ν c The second function, where the parameter k The compression modulus E c With the tensile modulus E t The absolute value of the ratio.

6. The stress calculation method for tensile and compressive dual modulus rock according to claim 1 is characterized in that: In the step of performing parameter sensitivity analysis and data fitting on the key parameters of the rock and obtaining parameter sensitivity analysis results and data fitting results, the method for performing the data fitting includes: A bivariate cubic polynomial function is used to fit the key parameters.

7. The stress calculation method for tensile and compressive dual modulus rock according to claim 1 is characterized in that: After the step of constructing a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results, the method further includes: Based on the mathematical model of horizontal stress in the mathematical model, the rock tensile strength is determined.

8. The stress calculation method for tensile and compressive dual modulus rock according to any one of claims 1 to 7, characterized in that: The constitutive equation of the numerical model in the principal stress coordinate system is: Formula (1): Formula (1) is the stress tensor, is the strain tensor, is the fourth-order elastic tensor in the principal stress coordinate system; Formula (2): Formula (3): In formula (2), E t represents the tensile modulus of rock, E c They represent the compression modulus of rock, ν t represents the Poisson's ratio associated with the stretched state, ν c represents the Poisson's ratio associated with the compressive state.

9. A stress calculation device for tension and compression dual modulus rock, characterized in that: include: The first acquisition module is used to perform parameter sensitivity analysis and data fitting on key parameters of rock based on the numerical model of rock tension and compression dual modulus assumption, and obtain parameter sensitivity analysis results and data fitting results; The model building module is used to build a mathematical model of rock stress based on the parameter sensitivity analysis results and the data fitting results.

10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the stress calculation method for tensile and compressive dual modulus rock as described in any one of claims 1 to 8 is implemented.

Images