Simulation calculation method and device for tunnel face grouting parameters based on extrusion deformation

Abstract

The present disclosure provides a method and device for simulating and calculating grouting parameters of a tunnel face based on extrusion deformation, which belongs to the field of data processing technology. The method includes: calculating a second extrusion deformation based on geological data of the tunnel site area, tunnel cross-section data and grouting parameters; the second extrusion deformation is the extrusion deformation of the tunnel face under the condition of grouting reinforcement of the face; calculating a second safety factor based on the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of grouting reinforcement of the face; if the second safety factor is less than or equal to a set threshold, adjusting the grouting parameters, and returning to execute the step of calculating the second extrusion deformation based on geological data of the tunnel site area, tunnel cross-section data and grouting parameters; if the second safety factor is greater than the set threshold, outputting the grouting parameter calculation result. The method and device for simulating and calculating grouting parameters of a tunnel face based on extrusion deformation provided by the present disclosure can improve the stability of the tunnel face and the safety of construction.

Description

Technical Field

[0001] The present disclosure belongs to the technical field of data processing, and more specifically, relates to a method and device for simulating and calculating grouting parameters of a tunnel face based on extrusion deformation. Background Art

[0002] The control of tunnel face stability during the construction of high-in-situ stress soft rock tunnels is the key to whether the tunnel construction can proceed smoothly. However, due to the high-energy geological environment and the particularity of soft rock mechanical properties in high-in-situ stress soft rock tunnels, the control of tunnel face stability is still a major engineering problem that needs to be solved. As a common and commonly used measure for strengthening the stability of tunnel faces, tunnel face grouting has been increasingly used in actual projects of high-in-situ stress soft rock tunnels. However, the current application and promotion of tunnel face grouting reinforcement measures have the following two shortcomings: (1) The determination of tunnel face grouting reinforcement parameters is mostly determined by engineering experience or engineering analogy, which makes it impossible to reasonably quantitatively evaluate the effect of tunnel face grouting reinforcement; (2) Most of the existing tunnel face grouting reinforcement stability evaluation methods are based on static equilibrium theory and are mainly applicable to soil or broken surrounding rock strata, but not to high-in-situ stress soft rock tunnels.

[0003] Therefore, it is urgent to propose a tunnel face grouting reinforcement parameter design method suitable for high ground stress soft rock, to achieve quantitative design and evaluation of tunnel face grouting reinforcement parameters, to provide theoretical guidance and suggestions for actual tunnel engineering, and to help improve tunnel face stability and construction safety. Summary of the invention

[0004] The purpose of the present invention is to provide a method and device for simulating and calculating tunnel face grouting parameters based on extrusion deformation, so as to realize the quantitative design and evaluation of tunnel face grouting reinforcement parameters, provide theoretical guidance and suggestions for actual tunnel engineering, and help improve the stability of tunnel face and construction safety.

[0005] A first aspect of the embodiments of the present disclosure provides a method for simulating and calculating grouting parameters of a tunnel face based on extrusion deformation, comprising:

[0006] The second extrusion deformation is calculated based on the geological data of the tunnel site, the tunnel cross-section data and the grouting parameters; the second extrusion deformation is the extrusion deformation of the tunnel face under the condition of grouting reinforcement of the tunnel face;

[0007] Calculating a second safety factor based on the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face grouting reinforcement;

[0008] If the second safety factor is less than or equal to the set threshold, adjust the grouting parameters, and return to the step of calculating the second extrusion deformation based on the tunnel site geological data, tunnel section data and the grouting parameters;

[0009] If the second safety factor is greater than a set threshold, the grouting parameter calculation result is output.

[0010] A second aspect of the disclosed embodiment provides a tunnel face grouting parameter simulation calculation device based on extrusion deformation, comprising:

[0011] The first calculation module is used to calculate the second extrusion deformation based on the geological data of the tunnel site, the tunnel section data and the grouting parameters; the second extrusion deformation is the extrusion deformation of the tunnel face under the condition of grouting reinforcement of the tunnel face;

[0012] A second calculation module is used to calculate a second safety factor based on the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face grouting reinforcement;

[0013] A parameter output module, if the second safety factor is less than or equal to a set threshold, adjusts the grouting parameters, and returns to execute the step of calculating the second extrusion deformation based on the geological data of the tunnel site, the tunnel section data and the grouting parameters;

[0014] If the second safety factor is greater than a set threshold, the grouting parameter calculation result is output.

[0015] According to a third aspect of an embodiment of the present disclosure, an electronic device is provided, comprising a memory, a processor, and a computer program stored in the memory and running on the processor, wherein when the processor executes the computer program, the steps of the above-mentioned method for simulating and calculating grouting parameters of a tunnel face based on extrusion deformation are implemented.

[0016] According to a fourth aspect of the embodiments of the present disclosure, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the above-mentioned method for simulating and calculating the grouting parameters of a tunnel face based on extrusion deformation are implemented.

[0017] The beneficial effects of the tunnel face grouting parameter simulation calculation method and device based on extrusion deformation provided by the embodiments of the present disclosure are:

[0018] In the disclosed embodiment, the stability of the tunnel face under the condition of face grouting reinforcement is judged based on the second safety factor. When the second safety factor is less than or equal to a set threshold, the grouting parameters are adjusted, and the second extrusion deformation is recalculated based on the geological data of the tunnel site, the tunnel cross-section data and the grouting parameters. Then, the second safety factor is calculated based on the second extrusion deformation until the second safety factor is greater than the set threshold, that is, the stability of the tunnel face meets the set requirements.

[0019] The disclosed embodiments take into account that the deformation, instability and failure evolution mechanism of the heading face of high-ground stress soft rock tunnels is significantly different from that of the heading face of tunnels in soil strata or broken strata. Based on the extrusion deformation of the heading face as the basis for judging the stability of the heading face of high-ground stress soft rock tunnels, a grouting reinforcement parameter design method for the heading face of high-ground stress soft rock tunnels based on the extrusion deformation of the heading face is proposed, which realizes the quantitative design of the heading face grouting parameters, provides theoretical guidance and suggestions for actual tunnel engineering, and helps to improve the stability of the tunnel heading face and the construction safety. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] In order to more clearly illustrate the technical solutions in the embodiments of the present disclosure, 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 disclosure. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0021] Figure 1 A schematic diagram of a flow chart of a tunnel face grouting parameter simulation calculation method based on extrusion deformation provided in an embodiment of the present disclosure;

[0022] Figure 2 A schematic flow chart of a tunnel face grouting parameter simulation calculation method based on extrusion deformation provided by another embodiment of the present disclosure;

[0023] Figure 3 A schematic diagram of calculating the equivalent circle of a horseshoe section provided by an embodiment of the present disclosure;

[0024] Figure 4 A schematic diagram of tunnel surrounding rock differentials provided by an embodiment of the present disclosure;

[0025] Figure 5 A structural block diagram of a tunnel face grouting parameter simulation calculation device based on extrusion deformation provided by an embodiment of the present disclosure;

[0026] Figure 6 A schematic block diagram of an electronic device provided by an embodiment of the present disclosure. DETAILED DESCRIPTION

[0027] In the following description, specific details such as specific system structures and technologies are provided for the purpose of illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present disclosure. However, it should be clear to those skilled in the art that the present disclosure may be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to avoid obstructing the description of the present disclosure with unnecessary details.

[0028] The following is an introduction to the technical terms involved in this embodiment.

[0029] Elastic-plastic boundary: In tunnel excavation, foundation treatment and other projects, the elastic-plastic properties of rock and soil need to be considered to evaluate the stability and safety of the project. Therefore, the concept of elastic-plastic boundary is often used. The elastic-plastic boundary refers to the boundary position where the material changes from the elastic state to the plastic state. At the elastic-plastic boundary, the stress of the material reaches the yield stress (that is, the critical stress value for the material to change from the elastic state to the plastic state). At the elastic-plastic boundary, the strain of the material will also change. In the elastic stage, the strain of the material is reversible. When unloaded, the material can return to its original shape and size; in the plastic stage, the strain of the material is irreversible. Even after unloading, the material will retain a certain permanent deformation.

[0030] Excavation radius: The excavation radius is the radius of the circular or arc area formed by excavation work in the construction of tunnels, mines, foundation projects, etc. In tunnel excavation, the excavation radius usually refers to the radius of the tunnel cross section.

[0031] Plastic zone radius: In geotechnical engineering, tunnel excavation will cause stress redistribution in the surrounding rock and soil. When the stress exceeds the yield strength of the rock and soil, a plastic zone will be formed around the tunnel. The plastic zone radius refers to the size of the area where the material undergoes plastic deformation under the action of force.

[0032] In order to make the purpose, technical solutions and advantages of the present disclosure more clear, specific embodiments will be described below in conjunction with the accompanying drawings.

[0033] Please refer to Figure 1 , Figure 1 A schematic flow chart of a tunnel face grouting parameter simulation calculation method based on extrusion deformation provided in an embodiment of the present disclosure, the method comprising:

[0034] S101: Calculate the second extrusion deformation based on the geological data of the tunnel site, the tunnel cross-section data and the grouting parameters; the second extrusion deformation is the extrusion deformation of the tunnel face under the condition of grouting reinforcement of the tunnel face.

[0035] In this embodiment, the geological data of the tunnel site area include ground stress, mechanical parameters of surrounding rock (such as elastic modulus, Poisson's ratio, cohesion, internal friction angle, etc.), which determine the deformation and failure characteristics of the surrounding rock under stress. Tunnel section data include tunnel section shape (such as circular, elliptical or horseshoe), section size (including specific values ​​such as tunnel height, width, radius, etc.). Grouting parameters include grouting material type (such as cement slurry, chemical slurry), material strength (including compressive strength, tensile strength, etc.), grouting pressure, grouting volume, etc. The geological data of the tunnel site area and the tunnel section data can be obtained through actual measurement, and the grouting parameters can be parameters pre-selected based on experience.

[0036] Considering the characteristics of the evolution of soft rock mechanical parameters, a method for calculating the extrusion deformation of the tunnel face under high ground stress soft rock under face grouting reinforcement is proposed based on the sphere theory, elastic-plastic theory and finite difference theory. This method is subject to the following assumptions:

[0037] (1) Based on the sphere theory, the tunnel face excavation process is considered as a hollow sphere, where the sphere radius is the tunnel face excavation radius. According to existing research results, the radial deformation of the hollow sphere is approximately equal to the extrusion deformation of the tunnel face. Therefore, the calculation of the extrusion deformation of the tunnel face is transformed into the calculation of the radial deformation of the hollow sphere.

[0038] (2) The tunnel surrounding rock is regarded as an isotropic and homogeneous material, and its mechanical parameters (including elastic modulus, Poisson's ratio, friction angle and cohesion) obey the soft rock mechanical parameter evolution model. The soft rock mechanical parameter evolution model is as follows:

[0039]

[0040] in, represents the plastic deviatoric strain, represents the strength stress ratio ( ), represents the confining pressure, Indicates the uniaxial compressive strength of the rock sample; represents the plastic deviatoric strain ratio, represents the critical plastic deviatoric strain; represents the dimensionless value of the elastic modulus, Indicates the real-time value of the elastic modulus, represents the elastic modulus of soft rock in the elastic stage under uniaxial conditions, represents the dimensionless value of Poisson's ratio, represents the real-time value of Poisson's ratio, represents the Poisson's ratio of soft rock in elastic stage under uniaxial condition, represents the dimensionless value of the friction angle, Indicates the real-time value of the friction angle, represents the friction angle of soft rock in elastic stage under uniaxial condition, is the dimensionless value of cohesion, Indicates cohesion, represents the cohesion of soft rock in the elastic stage under uniaxial conditions, Indicates a constant e Exponential operation with base ; is a proportional parameter, which is obtained by fitting the triaxial compression test data of soft rock under cyclic loading and unloading conditions under different confining pressures.

[0041] The soft rock mechanical parameter evolution model adopted in this embodiment can fully characterize the change process of mechanical parameters of soft rock in different deformation stages (elastic stage, softening stage, residual stage) with confining pressure and plastic deviatoric strain. Among them, the surrounding rock mechanical parameters include elastic modulus, Poisson's ratio, friction angle and cohesion. In addition, the above-mentioned soft rock mechanical parameter evolution model fully considers the coupling effect of confining pressure and plastic strain, and can accurately characterize the evolution process of soft rock mechanical parameters under the coupling effect of confining pressure and plastic strain, making the numerical simulation calculation results more accurate.

[0042] (3) The surrounding rock obeys the Mohr-Coulomb yield criterion, which is expressed as follows.

[0043]

[0044] in, is the radial stress, is the tangential stress, is the friction angle, For cohesion.

[0045] S102: Calculate a second safety factor based on the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face grouting reinforcement.

[0046] In this embodiment, considering that different tunnel faces can withstand different ultimate extrusion deformations, the stability of the tunnel face cannot be directly represented by the second extrusion deformation alone. Therefore, the ultimate extrusion deformation of the tunnel face ( ) and the second extrusion deformation ( ) as the second safety factor , the specific formula is as follows:

[0047]

[0048] Second safety factor It can intuitively characterize the stability of the tunnel face under the condition of face grouting reinforcement.

[0049] S103: If the second safety factor is less than or equal to the set threshold, adjust the grouting parameters, and return to the step of calculating the second extrusion deformation based on the tunnel site geological data, tunnel section data and grouting parameters;

[0050] If the second safety factor is greater than the set threshold, the grouting parameter calculation results are output.

[0051] In this embodiment, by comparing the second safety factor under grouting reinforcement conditions with the set threshold of the tunnel face safety factor, the safety and stability of the tunnel face of a high-stress soft rock tunnel after grouting reinforcement can be evaluated. If the second safety factor is greater than the set threshold, the tunnel face is safe and stable, and the design is completed; if the second safety factor is less than or equal to the set threshold, the tunnel face is unstable, and the tunnel face grouting reinforcement parameters need to be adjusted, such as increasing the grouting volume or grouting pressure, until the safety factor requirements are met.

[0052] Considering the complex geological conditions and construction factors, it is recommended to set the threshold to 1.5. If the surrounding rock conditions are poor and the ground stress level is high, the selection can be enlarged according to the actual situation to ensure a more adequate safety reserve.

[0053] From the above, it can be concluded that this embodiment judges the stability of the tunnel face under the condition of face grouting reinforcement based on the second safety factor. When the second safety factor is less than or equal to the set threshold, the grouting parameters are adjusted, and the second extrusion deformation is recalculated based on the geological data of the tunnel site, the tunnel section data and the grouting parameters, and then the second safety factor is calculated based on the second extrusion deformation, until the second safety factor is greater than the set threshold, that is, the stability of the tunnel face meets the set requirements.

[0054] The disclosed embodiment takes into account that the deformation, instability and failure evolution mechanism of the heading face of high-ground stress soft rock tunnels is significantly different from that of the heading face of tunnels in soil strata or broken strata. The heading face extrusion deformation is used as the basis for judging the stability of the heading face of high-ground stress soft rock tunnels. A grouting reinforcement parameter design method for the heading face of high-ground stress soft rock tunnels based on the heading face extrusion deformation is proposed, which realizes the quantitative design of the heading face grouting parameters, which is beneficial to improving the grouting reinforcement effect of the heading face.

[0055] Please refer to Figure 2 In one embodiment of the present disclosure, before the step of calculating the second extrusion deformation based on the tunnel site geological data, tunnel section data and grouting parameters, the tunnel face grouting parameter simulation calculation method based on extrusion deformation further includes:

[0056] The first extrusion deformation is calculated based on the geological data of the tunnel site and the tunnel section data; the first extrusion deformation is the extrusion deformation of the tunnel face without reinforcement;

[0057] Calculating a first safety factor based on a first extrusion deformation;

[0058] If the first safety factor is less than or equal to the set threshold, the step of calculating the second extrusion deformation based on the geological data of the tunnel site, the tunnel section data and the grouting parameters is executed.

[0059] In this embodiment, before calculating the grouting parameters, the first extrusion deformation of the tunnel face without reinforcement can be calculated based on the geological data of the tunnel site and the tunnel section data, and then the first safety factor can be calculated based on the first extrusion deformation. The stability of the tunnel face without reinforcement is judged based on the first safety factor. When the first safety factor is less than or equal to the set threshold, the grouting parameter calculation operation is performed to determine the appropriate grouting parameters, and the tunnel face is grout-reinforced to improve the stability of the tunnel face. When the first safety factor is greater than the set threshold, there is no need to perform the grouting parameter calculation operation to avoid unnecessary calculation operations.

[0060] In one embodiment of the present disclosure, the first extrusion deformation amount is calculated based on the geological data of the tunnel site and the tunnel section data, including:

[0061] The tunnel surrounding rock is divided into a plurality of first concentric rings from the first elastic-plastic boundary inward to the tunnel excavation radius; the first elastic-plastic boundary is the elastic-plastic boundary of the tunnel surrounding rock without reinforcement;

[0062] Starting from the first concentric ring, the stress information and strain information of each first concentric ring are calculated step by step until the first stop condition is met;

[0063] The stress information of the first concentric ring satisfies the stress equilibrium boundary condition, specifically:

[0064]

[0065] in, represents the radial stress of the first concentric ring, represents the tangential stress of the first concentric ring, represents the initial ground stress, represents the yield function, represents the friction angle of the first concentric ring, represents the cohesion of the first concentric ring;

[0066] The strain information of the first concentric ring is:

[0067]

[0068] in, represents the radial strain of the first concentric ring caused by elasticity, represents the tangential strain of the first concentric ring generated by elasticity, with the superscript Represents the elastic part, represents the elastic modulus of the first concentric ring, represents the Poisson's ratio of the first concentric ring;

[0069] No. The stress information of the first concentric ring is:

[0070]

[0071]

[0072]

[0073]

[0074]

[0075] in, is a set constant, indicating the width of the first concentric ring. Indicates The average radius of the first concentric rings, For the The friction angle of the first concentric ring, Indicates The cohesion of the first concentric rings;

[0076] No. The strain information of the first concentric ring is:

[0077]

[0078] in, The elasticity The radial strain increment of the first concentric ring is, The elasticity The tangential strain increment of the first concentric ring is, Indicates the first The tangential strain increment of the first concentric ring is, Represents the shear expansion coefficient, the specific value can be -1, Indicates the first The radial strain increment of the first concentric ring is, Indicates The total radial strain of the first concentric ring is, Indicates The total tangential strain of the first concentric ring;

[0079] The first stop condition is: The radial stress of the first concentric ring is equal to the support reaction force;

[0080] Based on The radius of the first plastic zone is calculated based on the inner diameter of the first concentric ring. The inner diameter of the first concentric ring and the The first extrusion deformation is calculated based on the total tangential strain of the first concentric rings.

[0081] In this embodiment, please refer to Figure 3 Before calculating the first extrusion deformation, the equivalent circle radius of the non-circular section is calculated based on the cross-sectional shape of the soft rock tunnel with high ground stress. The cross-sectional shapes of soft rock tunnels are mainly horseshoe-shaped and circular. The calculation step of the equivalent circle radius of the non-circular section mainly performs equivalent circle transformation on the non-circular section. If the cross-sectional shape of the soft rock tunnel is circular, skip this step. The equivalent circle radius obtained by the calculation method of the horseshoe equivalent circle radius ( ) is the height of the horseshoe section ( ) and span ( ) is 1 / 4 of the sum of . The specific calculation formula is as follows:

[0082] .

[0083] On this basis, the calculation process of the first extrusion deformation is as follows:

[0084] (S1) Perform differential and radius normalization processing on the tunnel surrounding rock: divide the surrounding rock from the elastic-plastic boundary inward into multiple first concentric rings to the tunnel excavation radius, substitute the confining pressure (i.e. radial stress) and plastic deviatoric strain of each first concentric ring into the soft rock mechanical parameter evolution model to calculate the corresponding mechanical parameters; and, according to symmetry, the mechanical parameters at different angles of the same ring are the same. Figure 4 As shown, Figure 4 In the figure, reference numeral 41 indicates a tunnel face, and reference numeral 42 indicates an elastic-plastic boundary.

[0085] The radius of the surrounding rock ring is normalized and the radius of the plastic zone is assumed to be , the expression of the normalized radius of the ring is as follows:

[0086]

[0087] in, Indicates The normalized radius of the first concentric ring, Indicates The inner diameter of the first concentric ring; the normalized radius of the first concentric ring (that is, the first ring) is .

[0088] (S2) Solve the stress and strain information of the first ring. According to the stress equilibrium boundary condition, the stress information of the outer boundary of the first ring can be solved. The expression of the stress equilibrium boundary condition is as follows:

[0089]

[0090] According to the theory of elastic-plastic mechanics, the expression of the strain information of the outer boundary of the first ring is as follows:

[0091]

[0092] Among them, since this ring is the outer boundary of the first ring (at the elastic-plastic boundary), according to the plastic increment theory, the radial strain and tangential strain generated by plasticity are both 0.

[0093] (S3) Any ring in the middle ( According to the yield criterion and equilibrium equation, any middle ring ( The expressions of radial stress and tangential stress of the ring are as follows:

[0094]

[0095]

[0096] in,

[0097]

[0098]

[0099]

[0100] represents the width of the first concentric ring, that is, the difference between the inner diameter and the outer diameter of the first concentric ring, which is set to:

[0101] .

[0102] (S4) Any ring in the middle ( According to the yield criterion, strain compatibility equation and plastic increment theory, the strain information of any middle ring ( The expressions of strain increment caused by elasticity, strain increment caused by plasticity and total strain increment are as follows:

[0103] .

[0104] (S5) Repeat steps (S3) and (S4) to solve the stress and strain information of the intermediate rings ring by ring until the first stop condition is met. The first stop condition is as follows:

[0105]

[0106] in, For the support reaction force, for a circular tunnel under uniform surrounding rock pressure, the axisymmetric elastic theory can be used to calculate the support reaction force.

[0107] (S6) Calculate the radius of the first plastic zone and the first extrusion deformation. The calculation expressions of the radius of the first plastic zone and the first extrusion deformation are as follows:

[0108]

[0109]

[0110] Among them, when When the radial stress of the first concentric ring meets the first stop condition, The inner diameter of the first concentric ring is the tunnel excavation radius .

[0111] It can be concluded from the above that, based on the sphere theory, this embodiment transforms the calculation of the extrusion deformation of the tunnel face into the calculation of the radial deformation of the hollow sphere, thereby realizing the quantitative calculation of the first extrusion deformation.

[0112] In one embodiment of the present disclosure, calculating the first safety factor based on the first extrusion deformation includes:

[0113] The ratio of the upper limit value of the extrusion deformation to the first extrusion deformation is determined as the first safety factor.

[0114] In this embodiment, the upper limit of the extrusion deformation (that is, the tunnel face limit extrusion deformation ( )) and the first extrusion deformation ( ) as the first safety factor , the specific formula is as follows:

[0115] .

[0116] Among them, the calculation of the extreme extrusion deformation of the tunnel face of high-in-situ stress soft rock tunnel can be carried out by numerical simulation method, and numerical simulation analysis can be carried out according to the specific high-in-situ stress soft rock tunnel engineering example to obtain the extreme extrusion deformation of the tunnel face corresponding to this engineering example. If the above numerical simulation method is not used, combined with the existing research results, it can be selected according to the span of the high-in-situ stress soft rock tunnel (as shown in Table 1) and the extrusion deformation control standard of the railway tunnel face in my country (as shown in Table 2).

[0117] Table 1 Classification standards for railway tunnel spans in my country

[0118]

[0119] Table 2 my country's railway tunnel face extrusion deformation control standards

[0120]

[0121] In one embodiment of the present disclosure, the second extrusion deformation is calculated based on the geological data of the tunnel site, the tunnel section data and the grouting parameters, including:

[0122] Performing multiple stress-strain calculation operations until a second stop condition is met;

[0123] Stress-strain calculation operations include:

[0124] Determine the second elastic-plastic boundary based on the set second plastic zone radius, where the second plastic zone radius is the plastic zone radius under the grouting reinforcement condition;

[0125] The tunnel surrounding rock is divided into A second concentric ring, until the radius of the tunnel excavation;

[0126] Starting from the first second concentric ring, the stress and strain information of each second concentric ring is calculated step by step until the first The second concentric ring, Radial stress of the second concentric ring;

[0127] Jordi The radial stress of the second concentric ring is less than the support reaction force, which reduces the radius of the second plastic zone;

[0128] The second stop condition is: The radial stress of the second concentric ring is equal to the support reaction force;

[0129] Based on The inner diameter of the second concentric ring and the The second extrusion deformation is calculated based on the total tangential strain of the second concentric rings.

[0130] In this embodiment, the method for calculating the extrusion deformation of the tunnel face under high ground stress soft rock tunnel grouting reinforcement (i.e., the second extrusion deformation) is also subject to the relevant assumptions in step S101, which will not be repeated here. The specific calculation steps of the second extrusion deformation are as follows:

[0131] (ST1) According to existing research, face grouting reinforcement measures can effectively reduce the range of the plastic zone. First, the radius of the second plastic zone under grouting reinforcement is set to .

[0132] (ST2) Differential processing within the plastic zone of surrounding rock. Differential processing is performed on the surrounding rock within the plastic zone of surrounding rock, and the surrounding rock is divided into The confining pressure (i.e. radial stress) and plastic deviatoric strain of each second concentric ring are substituted into the soft rock mechanical parameter evolution model to calculate the corresponding mechanical parameters; and according to symmetry, the mechanical parameters at different angles of the same ring are the same. It can be taken as 500. The differential step length expression is as follows:

[0133] .

[0134] (ST3) Starting from the first second concentric ring, the stress information and strain information of each second concentric ring is calculated step by step inward.

[0135] (ST4) By comparing The radial stress of the second concentric ring and support reaction ,like ,but , the calculation ends. ,but , reduce Then repeat steps (ST2) to (ST3) until ,at this time , the calculation is finished, As the final value of the radius of the second plastic zone.

[0136] (ST5) Calculate the extrusion deformation of the tunnel face after grouting reinforcement. The calculation expression of the extrusion deformation of the tunnel face after reinforcement is as shown in the following formula:

[0137]

[0138] in, is the second extrusion deformation, For the The tangential total strain of the second concentric ring is, when the second stop condition is met, The inner diameter of the second concentric ring is the tunnel face excavation radius .

[0139] In one embodiment of the present disclosure, starting from the first second concentric ring, gradually calculating the stress information and strain information of each second concentric ring inward includes:

[0140] The stress of the first second concentric ring satisfies the stress equilibrium boundary condition, specifically:

[0141]

[0142] in, represents the radial stress of the first second concentric ring, represents the tangential stress of the first second concentric ring, represents the initial ground stress, represents the yield function, represents the friction angle of the first second concentric ring, represents the cohesion of the first second concentric ring;

[0143] The strain information of the first second concentric ring is:

[0144]

[0145] in, represents the radial strain of the first second concentric ring caused by elasticity, represents the tangential strain of the first second concentric ring generated by elasticity, with the superscript Represents the elastic part, represents the elastic modulus of the first second concentric ring, represents the Poisson's ratio of the first second concentric ring;

[0146] Jordi The inner diameter of the second concentric ring is less than or equal to the grouting depth. The first cohesion is calculated based on the radial stress and soft rock mechanical parameter evolution model of the second concentric ring; the first adjustment parameter is calculated based on the first cohesion and the grouting depth; the first cohesion is adjusted based on the first adjustment parameter to obtain the first The cohesion of a second concentric ring;

[0147] No. The stress information of the second concentric ring is:

[0148]

[0149]

[0150]

[0151]

[0152]

[0153] in, represents the width of the second concentric ring, , represents the excavation radius of the tunnel face, Indicates the set radius of the second plastic zone, Indicates The average radius of the second concentric rings, Indicates The friction angle of the second concentric ring is Indicates The cohesion of a second concentric ring;

[0154] No. The strain information of the second concentric ring is:

[0155]

[0156] in, The elasticity The radial strain increment of the second concentric ring is, The elasticity The tangential strain increment of the second concentric ring is, Indicates the first The tangential strain increment of the second concentric ring is, Represents the shear expansion coefficient, the specific value can be -1, Indicates the first The radial strain increment of the second concentric ring is, Indicates The total radial strain of the second concentric ring is, Indicates The total tangential strain of the second concentric ring.

[0157] In this embodiment, step ST3 specifically includes:

[0158] (ST31) Middle arbitrary ring ( First, by comparing and face reinforcement depth ( That is, the grouting depth) size, judge Whether the grouting reinforcement scope has been entered.

[0159] like Greater than , that is, at this time, the grouting reinforcement range has not yet been entered, and any ring in the middle (the The expressions of radial stress and tangential stress of the ring are the same as the calculation formula of step S3 in the above embodiment.

[0160] like Less than or equal to , that is, at this time, the grouting reinforcement range has been entered, and the influence of grouting reinforcement on the tunnel face needs to be considered. The specific calculation method is: first, The radial stress of the second concentric ring The first cohesion (without considering the influence of grouting reinforcement) can be obtained by inputting the confining pressure into the soft rock mechanical parameter evolution model; then, it is multiplied by the cohesion increase coefficient (that is, the first adjustment parameter), and get The cohesion of the second concentric ring (considering the influence of grouting reinforcement) is calculated on this basis using the same calculation formula as step S3 in the above embodiment to obtain the cohesion of any middle ring (the Ring) radial stress and tangential stress expressions, cohesion enhancement coefficient The calculation formula is as follows:

[0161]

[0162] in, is the cohesion of the grouting body (Pa), is the grouting filling rate, is the grouting depth (m), is the tunnel face excavation diameter, is the failure angle of the tunnel face (°), which can be calculated by Perform calculations.

[0163] (ST32) Middle arbitrary ring ( The strain information of any middle ring ( The expressions for the strain increment caused by elasticity, the strain increment caused by plasticity and the total strain increment are the same as the calculation formula of step S4 in the above embodiment.

[0164] (ST33) Repeat steps (ST31) and (ST32) to solve the stress and strain information of the intermediate rings ring by ring until the final ring (the ring).

[0165] In summary, this embodiment takes into account that the deformation, instability and destruction evolution mechanism of the tunnel face of high-in-situ stress soft rock tunnel is significantly different from that of the tunnel face of soil strata or broken strata. Therefore, a stability evaluation model of the tunnel face of high-in-situ stress soft rock tunnel based on the extrusion deformation of the tunnel face is proposed, and a grouting reinforcement parameter design method of the tunnel face of high-in-situ stress soft rock tunnel based on the extrusion deformation of the tunnel face is established, realizing the quantitative design of the grouting parameters of the tunnel face. At the same time, the tunnel face extrusion deformation calculation method proposed in this embodiment takes into account the evolution process of soft rock mechanical parameters under high-in-situ stress conditions, making the calculation results more accurate and reasonable, and in line with the actual high-in-situ stress soft rock tunnel engineering.

[0166] Corresponding to the tunnel face grouting parameter simulation calculation method based on extrusion deformation in the above embodiment, Figure 5 This is a structural block diagram of a tunnel face grouting parameter simulation calculation device based on extrusion deformation provided by an embodiment of the present disclosure. For ease of explanation, only the parts related to the embodiment of the present disclosure are shown. Figure 5 The tunnel face grouting parameter simulation calculation device 20 based on extrusion deformation includes: a first calculation module 21, a second calculation module 22 and a parameter output module 23.

[0167] The first calculation module 21 is used to calculate the second extrusion deformation based on the geological data of the tunnel site, the tunnel section data and the grouting parameters; the second extrusion deformation is the extrusion deformation of the tunnel face under the condition of grouting reinforcement of the tunnel face;

[0168] A second calculation module 22 is used to calculate a second safety factor based on the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face grouting reinforcement;

[0169] The parameter output module 23 adjusts the grouting parameters if the second safety factor is less than or equal to the set threshold value, and returns to the step of calculating the second extrusion deformation based on the tunnel site geological data, the tunnel section data and the grouting parameters;

[0170] If the second safety factor is greater than the set threshold, the grouting parameter calculation results are output.

[0171] In one embodiment of the present disclosure, before the step of calculating the second extrusion deformation based on the tunnel site geological data, tunnel section data and grouting parameters, the first calculation module 21 is further used to:

[0172] The first extrusion deformation is calculated based on the geological data of the tunnel site and the tunnel section data; the first extrusion deformation is the extrusion deformation of the tunnel face without reinforcement;

[0173] Calculating a first safety factor based on a first extrusion deformation;

[0174] If the first safety factor is less than or equal to the set threshold, the step of calculating the second extrusion deformation based on the geological data of the tunnel site, the tunnel section data and the grouting parameters is executed.

[0175] In one embodiment of the present disclosure, the first calculation module 21 is further used to:

[0176] The tunnel surrounding rock is divided into a plurality of first concentric rings from the first elastic-plastic boundary inward to the tunnel excavation radius; the first elastic-plastic boundary is the elastic-plastic boundary of the tunnel surrounding rock without reinforcement;

[0177] Starting from the first concentric ring, the stress information and strain information of each first concentric ring are calculated step by step until the first stop condition is met;

[0178] The stress information of the first concentric ring satisfies the stress equilibrium boundary condition, specifically:

[0179]

[0180] in, represents the radial stress of the first concentric ring, represents the tangential stress of the first concentric ring, represents the initial ground stress, represents the yield function, represents the friction angle of the first concentric ring, represents the cohesion of the first concentric ring;

[0181] The strain information of the first concentric ring is:

[0182]

[0183] in, represents the radial strain of the first concentric ring caused by elasticity, represents the tangential strain of the first concentric ring generated by elasticity, with the superscript Represents the elastic part, represents the elastic modulus of the first concentric ring, represents the Poisson's ratio of the first concentric ring;

[0184] No. The stress information of the first concentric ring is:

[0185]

[0186]

[0187]

[0188]

[0189]

[0190] in, is a set constant, indicating the width of the first concentric ring. Indicates The average radius of the first concentric rings, For the The friction angle of the first concentric ring, Indicates The cohesion of the first concentric rings;

[0191] No. The strain information of the first concentric ring is:

[0192]

[0193] in, The elasticity The radial strain increment of the first concentric ring is, The elasticity The tangential strain increment of the first concentric ring is, Indicates the first The tangential strain increment of the first concentric ring is, Represents the shear expansion coefficient, the specific value can be -1, Indicates the first The radial strain increment of the first concentric ring is, Indicates The total radial strain of the first concentric ring is, Indicates The total tangential strain of the first concentric ring;

[0194] The first stop condition is: The radial stress of the first concentric ring is equal to the support reaction force;

[0195] Based on The radius of the first plastic zone is calculated based on the inner diameter of the first concentric ring. The inner diameter of the first concentric ring and the The first extrusion deformation is calculated based on the total tangential strain of the first concentric rings.

[0196] In one embodiment of the present disclosure, the first calculation module 21 is further configured to:

[0197] The ratio of the upper limit value of the extrusion deformation to the first extrusion deformation is determined as the first safety factor.

[0198] In one embodiment of the present disclosure, the first calculation module 21 is specifically used for:

[0199] Performing multiple stress-strain calculation operations until a second stop condition is met;

[0200] Stress-strain calculation operations include:

[0201] Determine the second elastic-plastic boundary based on the set second plastic zone radius, where the second plastic zone radius is the plastic zone radius under the grouting reinforcement condition;

[0202] The tunnel surrounding rock is divided into A second concentric ring, until the radius of the tunnel excavation;

[0203] Starting from the first second concentric ring, the stress and strain information of each second concentric ring is calculated step by step until the first The second concentric ring, Radial stress of the second concentric ring;

[0204] Jordi The radial stress of the second concentric ring is less than the support reaction force, which reduces the radius of the second plastic zone;

[0205] The second stop condition is: The radial stress of the second concentric ring is equal to the support reaction force;

[0206] Based on The inner diameter of the second concentric ring and the The second extrusion deformation is calculated based on the total tangential strain of the second concentric rings.

[0207] In one embodiment of the present disclosure, the second calculation module 21 is further configured to:

[0208] The stress of the first second concentric ring satisfies the stress equilibrium boundary condition, specifically:

[0209]

[0210] in, represents the radial stress of the first second concentric ring, represents the tangential stress of the first second concentric ring, represents the initial ground stress, represents the yield function, represents the friction angle of the first second concentric ring, represents the cohesion of the first second concentric ring;

[0211] The strain information of the first second concentric ring is:

[0212]

[0213] in, represents the radial strain of the first second concentric ring caused by elasticity, represents the tangential strain of the first second concentric ring generated by elasticity, with the superscript Represents the elastic part, represents the elastic modulus of the first second concentric ring, represents the Poisson's ratio of the first second concentric ring;

[0214] Jordi The inner diameter of the second concentric ring is less than or equal to the grouting depth. The first cohesion is calculated based on the radial stress and soft rock mechanical parameter evolution model of the second concentric ring; the first adjustment parameter is calculated based on the first cohesion and the grouting depth; the first cohesion is adjusted based on the first adjustment parameter to obtain the first The cohesion of a second concentric ring;

[0215] No. The stress information of the second concentric ring is:

[0216]

[0217]

[0218]

[0219]

[0220]

[0221] in, represents the width of the second concentric ring, , represents the excavation radius of the tunnel face, Indicates the set radius of the second plastic zone, Indicates The average radius of the second concentric rings, Indicates The friction angle of the second concentric ring is Indicates The cohesion of a second concentric ring;

[0222] No. The strain information of the second concentric ring is:

[0223]

[0224] in, The elasticity The radial strain increment of the second concentric ring is, The elasticity The tangential strain increment of the second concentric ring is, Indicates the first The tangential strain increment of the second concentric ring is, represents the shear dilatancy coefficient, Indicates the first The radial strain increment of the second concentric ring is, Indicates The total radial strain of the second concentric ring is, Indicates The total tangential strain of the second concentric ring.

[0225] See also Figure 6 , Figure 6 A schematic block diagram of an electronic device provided by an embodiment of the present disclosure. Figure 6 The electronic device 300 in the embodiment shown may include: one or more processors 301, one or more input devices 302, one or more output devices 303 and one or more memories 304. The processors 301, input devices 302, output devices 303 and memories 304 communicate with each other via a communication bus 305. The memory 304 is used to store computer programs, which include program instructions. The processor 301 is used to execute the program instructions stored in the memory 304. The processor 301 is configured to call the program instructions to execute the functions of each module / unit in the above-mentioned device embodiments, such as Figure 5 The functions of modules 21 to 23 are shown.

[0226] It should be understood that in the embodiment of the present disclosure, the processor 301 may be a central processing unit (CPU), and the processor may also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0227] The input device 302 may include a touch panel, a fingerprint collection sensor (for collecting the user's fingerprint information and fingerprint direction information), a microphone, etc., and the output device 303 may include a display (LCD, etc.), a speaker, etc.

[0228] The memory 304 may include a read-only memory and a random access memory, and provide instructions and data to the processor 301. A portion of the memory 304 may also include a non-volatile random access memory. For example, the memory 304 may also store information about the device type.

[0229] In a specific implementation, the processor 301, input device 302, and output device 303 described in the embodiments of the present disclosure can execute the second embodiment of the tunnel face grouting parameter simulation calculation method based on extrusion deformation provided in the embodiments of the present disclosure and the implementation method described in the second embodiment, and can also execute the implementation method of the electronic device described in the embodiments of the present disclosure, which will not be repeated here.

[0230] In another embodiment of the present disclosure, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program, wherein the computer program includes program instructions, and when the program instructions are executed by the processor, all or part of the processes in the above-mentioned embodiment method are implemented, and the computer program can also be completed by instructing the relevant hardware through the computer program. The computer program can be stored in a computer-readable storage medium, and when the computer program is executed by the processor, the steps of each of the above-mentioned method embodiments can be implemented. Among them, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form, etc. The computer-readable medium can include: any entity or device capable of carrying computer program code, recording medium, U disk, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium, etc.

[0231] The computer-readable storage medium may be an internal storage unit of the electronic device of any of the aforementioned embodiments, such as a hard disk or memory of the electronic device. The computer-readable storage medium may also be an external storage device of the electronic device, such as a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (SecureDigital, SD) card, a flash card (Flash Card), etc. equipped on the electronic device. Furthermore, the computer-readable storage medium may also include both an internal storage unit of the electronic device and an external storage device. The computer-readable storage medium is used to store computer programs and other programs and data required by the electronic device. The computer-readable storage medium may also be used to temporarily store data that has been output or is to be output.

[0232] Those of ordinary skill in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this disclosure.

[0233] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the electronic devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0234] In the several embodiments provided in the present application, it should be understood that the disclosed electronic devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic, for example, the division of units is only a logical function division, and there may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces or units, or it can be an electrical, mechanical or other form of connection.

[0235] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the embodiments of the present disclosure.

[0236] In addition, each functional unit in each embodiment of the present disclosure may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0237] The above are only specific embodiments of the present disclosure, but the protection scope of the present disclosure is not limited thereto. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed in the present disclosure, and these modifications or replacements should be included in the protection scope of the present disclosure. Therefore, the protection scope of the present disclosure should be based on the protection scope of the claims.

Claims

1. A tunnel face grouting parameter simulation calculation method based on extrusion deformation, characterized in that: include: The second extrusion deformation is calculated based on the geological data of the tunnel site, the tunnel cross-section data and the grouting parameters; the second extrusion deformation is the extrusion deformation of the tunnel face under the condition of grouting reinforcement of the tunnel face; Calculating a second safety factor based on the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face grouting reinforcement; If the second safety factor is less than or equal to the set threshold, the grouting parameters are adjusted, and the process returns to the step of calculating the second extrusion deformation based on the tunnel site geological data, the tunnel section data and the grouting parameters; If the second safety factor is greater than a set threshold, output the grouting parameter calculation result; Before the step of calculating the second extrusion deformation based on the tunnel site geological data, tunnel section data and grouting parameters, the tunnel face grouting parameter simulation calculation method based on extrusion deformation also includes: Calculate the first extrusion deformation based on the geological data of the tunnel site and the tunnel section data; the first extrusion deformation is the extrusion deformation of the tunnel face without reinforcement; Calculating a first safety factor based on the first extrusion deformation; If the first safety factor is less than or equal to a set threshold, the step of calculating the second extrusion deformation based on the tunnel site geological data, tunnel section data and grouting parameters is performed; The first extrusion deformation amount is calculated based on the geological data of the tunnel site and the tunnel section data, including: Dividing the tunnel surrounding rock into a plurality of first concentric circular rings from a first elastic-plastic boundary inwards to the tunnel excavation radius; the first elastic-plastic boundary is the elastic-plastic boundary of the tunnel surrounding rock without reinforcement; Starting from the first of the first concentric rings, gradually calculating the stress information and strain information of each of the first concentric rings inward until the first stop condition is met; The stress information of the first concentric ring satisfies the stress equilibrium boundary condition, specifically: in, represents the radial stress of the first concentric ring, represents the tangential stress of the first concentric ring, represents the initial ground stress, represents the yield function, represents the friction angle of the first concentric ring, represents the cohesion of the first concentric ring; The strain information of the first concentric ring is: in, represents the radial strain of the first concentric ring caused by elasticity, represents the tangential strain of the first concentric ring generated by elasticity, with the superscript Represents the elastic part, represents the elastic modulus of the first concentric ring, represents the Poisson's ratio of the first concentric ring; No. The stress information of the first concentric ring is: in, is a set constant, indicating the width of the first concentric ring. Indicates The average radius of the first concentric rings, For the The friction angle of the first concentric ring, Indicates The cohesion of the first concentric rings; No. The strain information of the first concentric ring is: in, The elasticity The radial strain increment of the first concentric ring is, The elasticity The tangential strain increment of the first concentric ring is, Indicates the first The tangential strain increment of the first concentric ring is, represents the shear dilatancy coefficient, Indicates the first The radial strain increment of the first concentric ring is, Indicates The total radial strain of the first concentric ring is, Indicates The total tangential strain of the first concentric ring; The first stop condition is: The radial stress of the first concentric ring is equal to the support reaction force; Based on The radius of the first plastic zone is calculated based on the inner diameter of the first concentric ring. The inner diameter of the first concentric ring and the The first extrusion deformation is calculated based on the total tangential strain of the first concentric rings.

2. The tunnel face grouting parameter simulation calculation method based on extrusion deformation according to claim 1 is characterized in that: The calculating the first safety factor based on the first extrusion deformation comprises: The ratio of the upper limit value of the extrusion deformation to the first extrusion deformation is determined as the first safety factor.

3. The tunnel face grouting parameter simulation calculation method based on extrusion deformation according to claim 1 is characterized in that: The second extrusion deformation amount is calculated based on the geological data of the tunnel site, the tunnel section data and the grouting parameters, including: Performing multiple stress-strain calculation operations until a second stop condition is met; The stress-strain calculation operation includes: Determining a second elastic-plastic boundary based on a set second plastic zone radius, wherein the second plastic zone radius is a plastic zone radius under grouting reinforcement conditions; The tunnel surrounding rock is divided into A second concentric ring, until the radius of the tunnel excavation; Starting from the first second concentric ring, the stress and strain information of each second concentric ring is calculated step by step until the first The second concentric ring, The radial stress and tangential total strain of the second concentric ring; Jordi The radial stress of the second concentric ring is less than the support reaction force, which reduces the radius of the second plastic zone; The second stop condition is: The radial stress of the second concentric ring is equal to the support reaction force; Based on The inner diameter of the second concentric ring and the The second extrusion deformation is calculated based on the total tangential strain of the second concentric rings.

4. The tunnel face grouting parameter simulation calculation method based on extrusion deformation according to claim 3 is characterized in that: The step of calculating the stress information and strain information of each second concentric ring gradually inward starting from the first second concentric ring includes: The stress of the first second concentric ring satisfies the stress equilibrium boundary condition, specifically: in, represents the radial stress of the first second concentric ring, represents the tangential stress of the first second concentric ring, represents the friction angle of the first second concentric ring, represents the cohesion of the first second concentric ring; The strain information of the first second concentric ring is: in, represents the radial strain of the first second concentric ring caused by elasticity, represents the tangential strain of the first second concentric ring generated by elasticity, with the superscript Represents the elastic part, represents the elastic modulus of the first second concentric ring, represents the Poisson's ratio of the first second concentric ring; Jordi The inner diameter of the second concentric ring is less than or equal to the grouting depth. The first cohesion is calculated based on the radial stress and soft rock mechanical parameter evolution model of the second concentric ring; the first adjustment parameter is calculated based on the first cohesion and the grouting depth; the first cohesion is adjusted based on the first adjustment parameter to obtain the first The cohesion of a second concentric ring; No. The stress information of the second concentric ring is: in, represents the width of the second concentric ring, , represents the excavation radius of the tunnel face, Indicates the set radius of the second plastic zone, Indicates The average radius of the second concentric rings, Indicates The friction angle of the second concentric ring is Indicates The cohesion of a second concentric ring; No. The strain information of the second concentric ring is: in, The elasticity The radial strain increment of the second concentric ring is, The elasticity The tangential strain increment of the second concentric ring is, Indicates the first The tangential strain increment of the second concentric ring is, represents the shear dilatancy coefficient, Indicates the first The radial strain increment of the second concentric ring is, Indicates The total radial strain of the second concentric ring is, Indicates The total tangential strain of the second concentric ring.

5. The tunnel face grouting parameter simulation calculation method based on extrusion deformation according to claim 1 is characterized in that: The calculating the second safety factor based on the second extrusion deformation comprises: The ratio of the upper limit value of the extrusion deformation to the second extrusion deformation is determined as the second safety factor.

6. A tunnel face grouting parameter simulation calculation device based on extrusion deformation, characterized in that: include: A first calculation module is used to calculate the second extrusion deformation based on the geological data of the tunnel site, the tunnel section data and the grouting parameters; The second extrusion deformation is the extrusion deformation of the tunnel face under the condition of face grouting reinforcement; A second calculation module, used for calculating a second safety factor based on the second extrusion deformation; The second safety factor is used to characterize the stability of the tunnel face under the condition of face grouting reinforcement; A parameter output module, if the second safety factor is less than or equal to a set threshold, adjusts the grouting parameters, and returns to execute the step of calculating the second extrusion deformation based on the geological data of the tunnel site, the tunnel section data and the grouting parameters; If the second safety factor is greater than a set threshold, output the grouting parameter calculation result; Before the step of calculating the second extrusion deformation based on the tunnel site geological data, tunnel section data and grouting parameters, the first calculation module is further used to: Calculate the first extrusion deformation based on the geological data of the tunnel site and the tunnel section data; the first extrusion deformation is the extrusion deformation of the tunnel face without reinforcement; Calculating a first safety factor based on the first extrusion deformation; If the first safety factor is less than or equal to a set threshold, the step of calculating the second extrusion deformation based on the tunnel site geological data, tunnel section data and grouting parameters is performed; The first calculation module is also used for: Dividing the tunnel surrounding rock into a plurality of first concentric circular rings from a first elastic-plastic boundary inwards to the tunnel excavation radius; the first elastic-plastic boundary is the elastic-plastic boundary of the tunnel surrounding rock without reinforcement; Starting from the first of the first concentric rings, gradually calculating the stress information and strain information of each of the first concentric rings inward until the first stop condition is met; The stress information of the first concentric ring satisfies the stress equilibrium boundary condition, specifically: in, represents the radial stress of the first concentric ring, represents the tangential stress of the first concentric ring, represents the initial ground stress, represents the yield function, represents the friction angle of the first concentric ring, represents the cohesion of the first concentric ring; The strain information of the first concentric ring is: in, represents the radial strain of the first concentric ring caused by elasticity, represents the tangential strain of the first concentric ring generated by elasticity, with the superscript Represents the elastic part, represents the elastic modulus of the first concentric ring, represents the Poisson's ratio of the first concentric ring; No. The stress information of the first concentric ring is: in, is a set constant, indicating the width of the first concentric ring. Indicates The average radius of the first concentric rings, For the The friction angle of the first concentric ring, Indicates The cohesion of the first concentric rings; No. The strain information of the first concentric ring is: in, The elasticity The radial strain increment of the first concentric ring is, The elasticity The tangential strain increment of the first concentric ring is, Indicates the first The tangential strain increment of the first concentric ring is, represents the shear dilatancy coefficient, Indicates the first The radial strain increment of the first concentric ring is, Indicates The total radial strain of the first concentric ring is, Indicates The total tangential strain of the first concentric ring; The first stop condition is: The radial stress of the first concentric ring is equal to the support reaction force; Based on The radius of the first plastic zone is calculated based on the inner diameter of the first concentric ring. The inner diameter of the first concentric ring and the The first extrusion deformation is calculated based on the total tangential strain of the first concentric rings.

Images