Simulation calculation method and device for anchor bolt parameters in tunnel face based on extrusion deformation

Abstract

The present disclosure provides a tunnel face anchor parameter simulation calculation method and device based on extrusion deformation, belonging to the field of data processing technology, the method includes: obtaining tunnel site geological data, tunnel cross-section data and anchor parameters; performing anchor parameter update operations multiple times until the first stop condition is met, and obtaining tunnel face anchor parameter simulation calculation results; anchor parameter update operations include: calculating the second extrusion deformation based on the tunnel site geological data, tunnel cross-section data and anchor parameters; calculating the second safety factor according to the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face anchor reinforcement; if the second safety factor is less than or equal to the set threshold, adjust the anchor parameters; the first stop condition is: the second safety factor is greater than the set threshold. The tunnel face anchor parameter simulation calculation method and device 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 anchor parameters of a tunnel face based on extrusion deformation. Background Art

[0002] Due to the high-energy geological environment and the particular mechanical properties of soft rocks in high-geo-stress soft rock tunnel projects, the stability control of the tunnel face in high-geo-stress soft rock tunnels is a major engineering problem that tunnel practitioners and scientific researchers urgently need to solve.

[0003] As a commonly used face stability reinforcement measure, glass fiber anchor bolts (hereinafter referred to as face anchor bolts) have been increasingly used in actual projects of soft rock tunnels with high geostress. However, the current application and promotion of face anchor bolt reinforcement measures have the following two shortcomings: (1) Most of the face anchor bolt reinforcement parameters are determined by engineering experience or engineering analogy, which makes it impossible to reasonably quantitatively evaluate the face anchor bolt reinforcement effect; (2) Most of the existing face anchor bolt reinforcement stability evaluation methods are based on static equilibrium theory and are mainly applicable to soil or broken surrounding rock strata, but have limited applicability to soft rock tunnels with high geostress. Summary of the invention

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

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

[0006] Obtain geological data of the tunnel site, tunnel section data and anchor parameters;

[0007] The anchor parameter update operation is performed multiple times until the first stop condition is met, and the simulation calculation result of the anchor parameter of the tunnel face is obtained;

[0008] The anchor parameter update operation includes:

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

[0010] Calculating a second safety factor according to the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face anchor reinforcement;

[0011] If the second safety factor is less than or equal to a set threshold, adjusting the anchor parameters;

[0012] The first stopping condition is: the second safety factor is greater than a set threshold.

[0013] A second aspect of the embodiments of the present disclosure provides a tunnel face anchor parameter simulation calculation device based on extrusion deformation, comprising:

[0014] Data acquisition module, used to obtain geological data of tunnel site, tunnel section data and anchor parameters;

[0015] A parameter calculation module is used to perform the anchor parameter update operation multiple times until the first stop condition is met, and obtain the tunnel face anchor parameter simulation calculation result;

[0016] The anchor parameter update operation includes:

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

[0018] Calculating a second safety factor according to the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face anchor reinforcement;

[0019] If the second safety factor is less than or equal to a set threshold, adjusting the anchor parameters;

[0020] The first stopping condition is: the second safety factor is greater than a set threshold.

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

[0022] In the disclosed embodiment, the stability of the tunnel face under the condition of face anchor reinforcement is judged based on the second safety factor. When the second safety factor is less than or equal to a set threshold, the anchor 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 anchor 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.

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

[0024] 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.

[0025] Figure 1 A schematic diagram of a flow chart of a method for simulating and calculating parameters of anchor bolts at a tunnel face based on extrusion deformation provided in an embodiment of the present disclosure;

[0026] Figure 2 A schematic diagram of a flow chart of a method for simulating and calculating parameters of anchor bolts at a tunnel face based on extrusion deformation provided by another embodiment of the present disclosure;

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

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

[0029] Figure 5 A schematic diagram of the relationship between the anchor rod length and the surrounding rock plastic zone provided in an embodiment of the present disclosure; wherein, Figure 5 (a) is a schematic diagram of anchor length = plastic zone radius. Figure 5 (b) is a schematic diagram of anchor length > plastic zone radius. Figure 5 (c) is a schematic diagram of anchor length < plastic zone radius;

[0030] Figure 6 This is a structural block diagram of a tunnel face anchor parameter simulation calculation device based on extrusion deformation provided by an embodiment of the present disclosure. DETAILED DESCRIPTION

[0031] 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.

[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] The following is an introduction to the technical terms involved in this embodiment.

[0034] 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.

[0035] 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.

[0036] 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.

[0037] Please refer to Figure 1-Figure 2 , Figure 1 A schematic flow chart of a method for simulating and calculating parameters of anchor bolts at a tunnel face based on extrusion deformation provided in an embodiment of the present disclosure, the method comprising:

[0038] S101: Obtaining geological data of the tunnel site, tunnel section data and anchor parameters.

[0039] 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 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.). Anchor rod parameters include anchor rod length, anchor rod diameter, anchor rod spacing, etc. The geological data of the tunnel site area and the tunnel section data can be obtained through actual measurement, and the anchor rod parameters can be parameters pre-selected based on experience.

[0040] S102: performing the anchor parameter update operation multiple times until the first stop condition is met, and obtaining the tunnel face anchor parameter simulation calculation result.

[0041] The anchor parameter update operation includes:

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

[0043] The second safety factor is calculated according to the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face anchor reinforcement;

[0044] If the second safety factor is less than or equal to the set threshold, adjust the anchor parameters;

[0045] The first stopping condition is: the second safety factor is greater than a set threshold.

[0046] In this embodiment, 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 sphere theory, elastic-plastic theory and finite difference theory. This method is subject to the following assumptions:

[0047] (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 (also known as the second extrusion deformation) is transformed into the calculation of the radial deformation of the hollow sphere.

[0048] (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:

[0049]

[0050] 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, Represents cohesion, represents the cohesion of soft rock in the elastic stage under uniaxial conditions, Represents exponential operation with constant e as 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.

[0051] The soft rock mechanical parameter evolution model used 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. 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.

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

[0053]

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

[0055] Considering that different tunnel faces can withstand different ultimate extrusion deformations, the stability of the tunnel face cannot be directly characterized by the second extrusion deformation alone. Therefore, based on the second extrusion deformation, the ultimate extrusion deformation of the tunnel face ( ) and the second extrusion deformation ( ) as the second safety factor , the specific formula is as follows:

[0056]

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

[0058] 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).

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

[0060]

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

[0062]

[0063] By comparing the second safety factor under anchor reinforcement conditions with the set threshold of the face safety factor, the safety and stability of the face of a soft rock tunnel with high ground stress after anchor reinforcement can be evaluated. If the second safety factor is greater than the set threshold, the face is safe and stable, and the parameter calculation is completed; if the second safety factor is less than or equal to the set threshold, the face is unstable, and the face anchor reinforcement parameters need to be adjusted until the safety factor requirements are met.

[0064] 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.

[0065] From the above, it can be concluded that this embodiment judges the stability of the tunnel face under the condition of face anchor reinforcement based on the second safety factor. When the second safety factor is less than or equal to the set threshold, the anchor parameters are adjusted, and the second extrusion deformation is recalculated based on the tunnel site geological data, tunnel section data and anchor 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.

[0066] The disclosed embodiment takes into account that the evolution mechanism of deformation, instability and failure of the tunnel face in high-ground stress soft rock tunnels is significantly different from the evolution mechanism of deformation, instability and failure of the tunnel face in soil strata or broken strata. The extrusion deformation of the tunnel face is used as the basis for judging the stability of the tunnel face in high-ground stress soft rock tunnels. A method for designing anchor reinforcement parameters for the tunnel face in high-ground stress soft rock tunnels based on the extrusion deformation of the tunnel face is proposed, which realizes quantitative design and evaluation of anchor reinforcement parameters for the tunnel face, provides theoretical guidance and suggestions for actual tunnel engineering, and helps to improve the stability of the tunnel face and construction safety.

[0067] In one embodiment of the present disclosure, before performing the anchor parameter update operation multiple times, the tunnel face anchor parameter simulation calculation method based on extrusion deformation further includes:

[0068] 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;

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

[0070] If the first safety factor is less than or equal to the set threshold, the anchor parameter update operation is performed.

[0071] In this embodiment, before calculating the reinforcement system 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 grouting reinforcement is judged based on the first safety factor. When the first safety factor is less than or equal to the set threshold, the reinforcement system parameter calculation operation is performed to determine the appropriate advanced reinforcement system parameters, and the tunnel face is grouting 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 advanced reinforcement system parameter calculation operation to avoid unnecessary calculation operations.

[0072] 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:

[0073] .

[0074] After calculating the equivalent circle radius of the non-circular section, the tunnel surrounding rock is differentiated and radius normalized: the surrounding rock is divided into multiple first concentric rings from the elastic-plastic boundary to the tunnel excavation radius, and the confining pressure (i.e. radial stress) and plastic deviatoric strain of each first 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. Figure 4 As shown, Figure 4 Reference numeral 41 indicates a tunnel face, and reference numeral 42 indicates an elastic-plastic boundary.

[0075] Then, the radius of the surrounding rock ring is normalized, assuming that the radius of the plastic zone is , the expression of the normalized radius of the ring is as follows:

[0076]

[0077] 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 .

[0078] On this basis, the first extrusion deformation is calculated based on the geological data of the tunnel site and the tunnel section data, including:

[0079] Divide the tunnel surrounding rock into a plurality of first concentric circular rings from the first elastic-plastic boundary inwards to the tunnel excavation radius;

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

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

[0082]

[0083] 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;

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

[0085]

[0086] 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;

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

[0088]

[0089]

[0090]

[0091]

[0092]

[0093] 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 is Indicates The cohesion of the first concentric rings;

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

[0095]

[0096] 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;

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

[0098] 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.

[0099] In this embodiment, the calculation expressions of the first plastic zone radius and the first extrusion deformation are as follows:

[0100]

[0101]

[0102] 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 .

[0103] On the basis of obtaining the first extrusion deformation, the upper limit of the extrusion deformation (i.e., the tunnel face limit extrusion deformation ( )) and the first extrusion deformation ( ) as the first safety factor , the specific formula is as follows:

[0104] .

[0105] 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 anchor parameters, including:

[0106] If the bolt length is equal to the radius of the plastic zone, the surrounding rock is divided into A second concentric ring, until the radius of the tunnel excavation;

[0107] Starting from the first second concentric ring, the stress information, strain information and anchor force information of each second concentric ring are calculated step by step until the first second concentric ring. A second concentric ring;

[0108] 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.

[0109] In this embodiment, if Figure 5 As shown in the figure, due to the different lengths of anchor rods, according to the relationship between anchor rod length and plastic zone of surrounding rock, it can be divided into three cases: (1) The first case: anchor rod length = plastic zone radius; (2) The second case: anchor rod length > plastic zone radius; (3) The third case: anchor rod length < plastic zone radius. Therefore, the situation after anchor rod reinforcement should be judged first. It is assumed that the plastic zone radius of surrounding rock after anchor rod reinforcement is equal to the anchor rod length ( ), perform differential calculation based on the anchor length and solve the radial stress That is, starting from the first second concentric ring, the anchor length is used as the radius of the plastic zone, and the stress information, strain information and anchor force information of each second concentric ring are calculated step by step until the first By comparing the The radial stress of the second concentric ring and support reaction , determine which of the above three situations the anchor reinforcement belongs to, the judgment method is as follows:

[0110] First case: If ,but ; Anchor length = plastic zone radius;

[0111] Second case: If ,but ; Anchor length > plastic zone radius;

[0112] The third case: If ,but ; Anchor rod length < plastic zone radius.

[0113] Among them, for a circular tunnel under uniform surrounding rock pressure, the support reaction force can be obtained by combining the axisymmetric elastic theory with the tunnel section data. Figure 5 Reference numeral 41 indicates a tunnel face, reference numeral 42 indicates an elastic-plastic boundary, and reference numeral 43 indicates an anchor bolt.

[0114] For the first case: anchor length = plastic zone radius, since it already satisfies The condition is calculated to When a second concentric ring is formed, The inner diameter of the second concentric ring and the The second extrusion deformation is calculated by the total tangential strain of the second concentric ring. The specific calculation formula is:

[0115] .

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

[0117] In this embodiment, for the first case: anchor rod length = plastic zone radius, the calculation process of the stress information, strain information and anchor rod force information of each second concentric ring is as follows:

[0118] (S11) Differential processing within the plastic zone of the surrounding rock. The surrounding rock is divided into From the second concentric ring to the outer radius of the tunnel, the mechanical parameters of the surrounding rock of each ring obey the soft rock mechanical parameter evolution model; and, based on 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:

[0119]

[0120] in, represents the width of the second concentric ring, represents the excavation radius of the tunnel face, Indicates the length of the anchor rod.

[0121] (S12) Solve the stress information of the first ring. According to the elastic-plastic boundary conditions, solve the radial stress and tangential stress at the elastic-plastic boundary (that is, the first second concentric ring):

[0122]

[0123] 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.

[0124] (S13) Solve the strain information of the first ring. Since it is at the elastic-plastic boundary, based on the plastic increment theory, the strain information generated by plasticity is 0 at this time, so the strain information at the elastic-plastic boundary is as follows:

[0125]

[0126] 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.

[0127] (S14) Solve the first ring anchor force information. Assuming that the anchor and the surrounding rock deform cooperatively, the anchor force expression at the elastic-plastic boundary is as follows:

[0128]

[0129] in, Indicates the anchor force information of the first second concentric ring. represents the cross-sectional area of ​​the anchor rod, represents the elastic modulus of the anchor, represents the total radial strain of the first second concentric ring, .

[0130] (S15) Any ring in the middle ( According to the equilibrium equation and yield criterion, any middle ring ( The expression of the ring stress information is as follows:

[0131]

[0132]

[0133] in,

[0134]

[0135]

[0136]

[0137] in, represents the width of the second concentric ring, , represents the excavation radius of the tunnel face, Indicates the anchor length, Indicates The average radius of the second concentric rings, Indicates The strain increment of the second concentric ring is, For the The friction angle of the second concentric ring is Indicates The cohesion of the second concentric ring, Indicates the vertical spacing of anchor rods (m), Indicates the horizontal spacing of anchor rods (m).

[0138] (S16) Any ring in the middle ( The strain information of any ring in the middle ( The expressions of strain increment caused by elasticity, strain increment caused by plasticity and total strain increment are as follows:

[0139]

[0140] 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.

[0141] (S17) Any ring in the middle ( The anchor force information of any intermediate ring is as follows:

[0142] .

[0143] Repeat steps (S15) to (S17) to solve the stress and strain information of the intermediate rings ring by ring until the final ring (the ring).

[0144] 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 anchor parameters, including:

[0145] If the bolt length is greater than the plastic zone radius, the surrounding rock is divided into A second concentric ring, until the radius of the tunnel excavation;

[0146] The set boundary radial stress is determined as the radial stress of the first second concentric ring. Starting from the first second concentric ring, the stress information, strain information and anchor force information of each second concentric ring are calculated step by step until the second concentric ring is reached. A second concentric ring;

[0147] Jordi The radial stress of the second concentric ring is less than the support reaction force, the boundary radial stress is increased, and the process returns to the step of calculating the stress information, strain information, and anchor force information of each second concentric ring step by step from the first second concentric ring inward;

[0148] Jordi The radial stress of the second concentric ring is greater than the support reaction force, the boundary radial stress is reduced, and the process returns to the step of calculating the stress information, strain information, and anchor force information of each second concentric ring step by step from the first second concentric ring inward;

[0149] Jordi The radial stress of the second concentric ring is equal to the support reaction force. 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.

[0150] In this embodiment, for the second case: anchor rod length > plastic zone radius, due to , multiple iterations are required to continuously adjust the boundary radial stress until it meets Specifically, when calculating to the When there is a second concentric ring, it can be determined and ,like Greater than , then after increasing the boundary radial stress, start the calculation again from the first second concentric ring; if Less than , then after reducing the boundary radial stress, the calculation starts again from the first second concentric ring until , and the operation ends.

[0151] At this time, the calculation expression of the extrusion deformation of the reinforced tunnel face (that is, the second extrusion deformation) is as follows:

[0152]

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

[0154] In this embodiment, for the second case: anchor rod length>plastic zone radius, the calculation process of the stress information, strain information and anchor rod force information of each second concentric ring is as follows:

[0155] (S21) Surrounding rock differential processing. The surrounding rock within the length of the anchor is differentially processed, and the surrounding rock is divided from the inner end of the anchor into From the second concentric ring to the outer radius of the tunnel, the mechanical parameters of the surrounding rock of each ring obey the soft rock mechanical parameter evolution model; and, based on 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:

[0156]

[0157] in, represents the width of the second concentric ring, represents the excavation radius of the tunnel face, Indicates the length of the anchor rod.

[0158] (S22) The first ring stress information is solved. Assuming the radial stress at the boundary , and solve for the tangential stress according to the following formula:

[0159]

[0160] (S23) Solve the strain information of the first ring. Since it is in an elastic state, based on the plastic increment theory, the strain information generated by plasticity is 0 at this time, so the strain expression at the elastic-plastic boundary is as follows:

[0161]

[0162] 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.

[0163] (S24) Solve the first ring anchor force information. Assuming that the anchor and the surrounding rock deform cooperatively, the anchor force expression at the elastic-plastic boundary is as follows:

[0164]

[0165] in, Indicates the anchor force information of the first second concentric ring. represents the cross-sectional area of ​​the anchor rod, represents the elastic modulus of the anchor, represents the total radial strain of the first second concentric ring, .

[0166] (S25) Any ring in the middle ( According to the equilibrium equation and yield criterion, any middle ring ( The expression of the ring stress information is as follows:

[0167]

[0168] in,

[0169]

[0170] in, represents the width of the second concentric ring, , represents the excavation radius of the tunnel face, Indicates the anchor length, Indicates The average radius of the second concentric rings, Indicates The strain increment of the second concentric ring is, Indicates the vertical spacing of anchor rods (m), Indicates the horizontal spacing of anchors (m);

[0171] Discrimination and Relationship, determine whether the surrounding rock is in an elastic or plastic state at this time.

[0172] like Greater than It is still in elastic state. Less than The plastic state.

[0173] in, is the reference value of radial stress, represents the reference value of tangential stress, and The specific value of is based on the stress equilibrium boundary condition, which is:

[0174]

[0175] If it is elastic,

[0176] .

[0177] Any ring in the middle ( The calculation formula for the strain information of the ring) is the same as the calculation formula in step (S16) in the first case.

[0178] If it is a plastic state, let = , update the current ring radius differential expression as follows:

[0179]

[0180]

[0181] in, Indicates The average tangential stress of the second concentric ring, i.e. The average value of the tangential stress of the inner and outer rings of the second concentric ring, Indicates The average radial stress of the second concentric ring, i.e. The average value of the radial stress of the inner and outer rings of the second concentric ring.

[0182] Re-determine and :

[0183]

[0184] .

[0185] Then, using the calculation formula of step (S16) in the first case, calculate the middle arbitrary ring (the ring) strain information.

[0186] (S26) Repeat step (S25) to solve the stress and strain information of the intermediate rings ring by ring until the final ring (the ring).

[0187] 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 anchor parameters, including:

[0188] If the bolt length is less than the radius of the plastic zone, the surrounding rock is divided into A second concentric ring, until the radius of the tunnel excavation;

[0189] Starting from the first second concentric ring, the stress information, strain information and anchor force information of each second concentric ring are calculated step by step inward based on the set radius of the second plastic zone until the first second concentric ring is reached. A second concentric ring;

[0190] Jordi The radial stress of the second concentric ring is less than the support reaction force, the radius of the second plastic zone is increased, and the process returns to the step of calculating the stress information, strain information, and anchor force information of each second concentric ring gradually inward starting from the first second concentric ring;

[0191] Jordi The radial stress of the second concentric ring is greater than the support reaction force, the radius of the second plastic zone is reduced, and the process returns to the step of calculating the stress information, strain information, and anchor force information of each second concentric ring step by step from the first second concentric ring inward;

[0192] Jordi The radial stress of the second concentric ring is equal to the support reaction force. 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.

[0193] In this embodiment, for the third case: anchor rod length < plastic zone radius, due to , multiple iterations are required to continuously adjust the radius of the second plastic zone until Specifically, when calculating to the When there is a second concentric ring, it can be determined and ,like Greater than , then increase Then, start the calculation again from the first second concentric ring; if Less than , then reduce After that, start calculating again from the first second concentric ring until , and the operation ends.

[0194] At this time, the calculation expression of the extrusion deformation of the reinforced tunnel face (that is, the second extrusion deformation) is as follows:

[0195]

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

[0197] In this embodiment, for the third case: anchor rod length < plastic zone radius, the calculation process of the stress information, strain information and anchor rod force information of each second concentric ring is as follows:

[0198] (S31) Differential processing within the plastic zone of the surrounding rock. Assume that the radius of the plastic zone of the surrounding rock is , the surrounding rock in the plastic zone is differentially processed, and the surrounding rock is divided into From the concentric rings to the outer radius of the tunnel, the mechanical parameters of the surrounding rock of each ring obey the soft rock mechanical parameter evolution model; and, based on 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:

[0199]

[0200] in, represents the width of the second concentric ring, Indicates the excavation radius of the tunnel face.

[0201] (S32) Solve the stress information of the first ring. According to the elastic-plastic boundary conditions, solve the radial stress and tangential stress at the elastic-plastic boundary:

[0202]

[0203] 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.

[0204] (S33) Solve the strain information of the first ring. Since it is at the elastic-plastic boundary, based on the plastic increment theory, the strain information generated by plasticity is 0 at this time, so the strain information at the elastic-plastic boundary is as follows:

[0205]

[0206] 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.

[0207] (S34) Solve the first ring anchor force information. Assuming that the anchor and the surrounding rock deform in coordination and have not yet entered the anchor support range, the anchor force expression is as follows:

[0208] .

[0209] (S35) Any ring in the middle ( Ring) radius information is solved. and The expression is as follows:

[0210]

[0211] .

[0212] (S36) Determine whether anchor rods are involved.

[0213] like > Then the anchor has not yet participated;

[0214] like < Then the anchor rod is already involved.

[0215] (S37) If the anchor rod is not involved, determine the stress information differential expression:

[0216]

[0217]

[0218]

[0219]

[0220]

[0221] Using the calculation formula of step (S16) in the first case, calculate the middle arbitrary ring ( ring) strain information.

[0222] (S38) If the anchor rod is already involved, The inner diameter of the second concentric ring Updated to anchor length:

[0223]

[0224]

[0225] Determine Stress information for the second concentric ring:

[0226]

[0227]

[0228]

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

[0230] Using the calculation formula of step (S16) in the first case, calculate the middle arbitrary ring ( ring) strain information.

[0231] Re-determine and :

[0232]

[0233]

[0234] Based on the re-determination and Calculate the The stress and strain information of the second concentric ring.

[0235] (S39) Repeat steps (S35) to (S38) to solve the stress and strain information of the intermediate rings ring by ring until the final ring (the ring).

[0236] Corresponding to the tunnel face anchor parameter simulation calculation method based on extrusion deformation in the above embodiment, Figure 6 This is a structural block diagram of a tunnel face anchor 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 6 The tunnel face anchor parameter simulation calculation device 20 based on extrusion deformation includes: a data acquisition module 21 and a parameter calculation module 22.

[0237] The data acquisition module 21 is used to acquire the geological data of the tunnel site, the tunnel section data and the anchor parameters;

[0238] The parameter calculation module 22 is used to perform the anchor parameter update operation multiple times until the first stop condition is met, and obtain the tunnel face anchor parameter simulation calculation result;

[0239] The anchor parameter update operation includes:

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

[0241] The second safety factor is calculated according to the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face anchor reinforcement;

[0242] If the second safety factor is less than or equal to the set threshold, adjust the anchor parameters;

[0243] The first stopping condition is: the second safety factor is greater than a set threshold.

[0244] In one embodiment of the present disclosure, the parameter calculation module 22 is specifically used to:

[0245] 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;

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

[0247] If the first safety factor is less than or equal to the set threshold, the anchor parameter update operation is performed.

[0248] In one embodiment of the present disclosure, the parameter calculation module 22 is specifically used to:

[0249] Divide the tunnel surrounding rock into a plurality of first concentric circular rings from the first elastic-plastic boundary inwards to the tunnel excavation radius;

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

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

[0252]

[0253] 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;

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

[0255]

[0256] 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;

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

[0258]

[0259]

[0260]

[0261]

[0262]

[0263] 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;

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

[0265]

[0266] 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;

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

[0268] 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.

[0269] In one embodiment of the present disclosure, the parameter calculation module 22 is specifically used to:

[0270] If the bolt length is equal to the radius of the plastic zone, the surrounding rock is divided into A second concentric ring, until the radius of the tunnel excavation;

[0271] Starting from the first second concentric ring, the stress information, strain information and anchor force information of each second concentric ring are calculated step by step until the first second concentric ring. A second concentric ring;

[0272] 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.

[0273] In one embodiment of the present disclosure, the parameter calculation module 22 is specifically used to:

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

[0275]

[0276] 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;

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

[0278]

[0279] 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;

[0280] The anchor force information of the first second concentric ring is:

[0281]

[0282] in, Indicates the anchor force information of the first second concentric ring. represents the cross-sectional area of ​​the anchor rod, represents the elastic modulus of the anchor, represents the total radial strain of the first second concentric ring, ;

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

[0284]

[0285]

[0286] in,

[0287]

[0288]

[0289]

[0290] in, represents the width of the second concentric ring, , represents the excavation radius of the tunnel face, Indicates the anchor length, Indicates The average radius of the second concentric rings, Indicates The strain increment of the second concentric ring is, For the The friction angle of the second concentric ring is, Indicates The cohesion of the second concentric ring, Indicates the vertical spacing of anchor rods, Indicates the horizontal spacing of anchor rods;

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

[0292]

[0293] 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;

[0294] No. The anchor force information of the second concentric ring is:

[0295] .

[0296] 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 anchor parameter simulation calculation method based on extrusion deformation, characterized in that: include: Obtain geological data of the tunnel site, tunnel section data and anchor parameters; The anchor parameter update operation is performed multiple times until the first stop condition is met, and the simulation calculation result of the anchor parameter of the tunnel face is obtained; The anchor parameter update operation includes: The second extrusion deformation is calculated based on the geological data of the tunnel site, the tunnel section data and the anchor parameters; the second extrusion deformation is the extrusion deformation of the tunnel face under the condition of anchor reinforcement of the tunnel face; Calculating a second safety factor according to the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face anchor reinforcement; If the second safety factor is less than or equal to a set threshold, adjusting the anchor parameters; The first stop condition is: the second safety factor is greater than a set threshold; Before performing the anchor parameter update operation multiple times, the tunnel face anchor parameter simulation calculation method based on extrusion deformation further 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, executing the anchor parameter update operation; The first extrusion deformation amount is calculated based on the geological data of the tunnel site and the tunnel section data, including: Divide the tunnel surrounding rock into a plurality of first concentric circular rings from the first elastic-plastic boundary inwards to the tunnel excavation radius; 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 second 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 is 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 second 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 anchor parameter simulation calculation method based on extrusion deformation according to claim 1 is characterized in that: The calculating of the second extrusion deformation based on the tunnel site geological data, tunnel section data and anchor rod parameters includes: If the bolt length is equal to the radius of the plastic zone, the 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 information, strain information and anchor force information of each second concentric ring are calculated step by step until the first second concentric ring. A second concentric ring; 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.

3. The tunnel face anchor parameter simulation calculation method based on extrusion deformation according to claim 2 is characterized in that: The method starts from the first second concentric ring and gradually calculates the stress information, strain information and anchor force information of each second concentric ring inward, including: 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 caused by elasticity, represents the elastic modulus of the first second concentric ring, represents the Poisson's ratio of the first second concentric ring; The anchor force information of the first second concentric ring is: in, Indicates the anchor force information of the first second concentric ring. represents the cross-sectional area of ​​the anchor rod, represents the elastic modulus of the anchor, represents the total radial strain of the first second concentric ring, ; No. The stress information of the second concentric ring is: in, in, represents the width of the second concentric ring, , represents the excavation radius of the tunnel face, Indicates the anchor length, Indicates The average radius of the second concentric rings, Indicates The strain increment of the second concentric ring is, For the The friction angle of the second concentric ring is, Indicates The cohesion of the second concentric ring, Indicates the vertical spacing of anchor rods, Indicates the horizontal spacing of anchor rods; 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, 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; No. The anchor force information of the second concentric ring is: 。 4. The tunnel face anchor parameter simulation calculation method based on extrusion deformation according to claim 1 is characterized in that: The calculating of the second extrusion deformation based on the tunnel site geological data, tunnel section data and anchor rod parameters includes: If the bolt length is greater than the plastic zone radius, the surrounding rock is divided into A second concentric ring, until the radius of the tunnel excavation; The set boundary radial stress is determined as the radial stress of the first second concentric ring. Starting from the first second concentric ring, the stress information, strain information and anchor force information of each second concentric ring are calculated step by step until the second concentric ring is reached. A second concentric ring; Jordi The radial stress of the second concentric ring is less than the support reaction force, the boundary radial stress is increased, and the process returns to the step of calculating the stress information, strain information, and anchor force information of each second concentric ring step by step from the first second concentric ring inward; Jordi The radial stress of the second concentric ring is greater than the support reaction force, the boundary radial stress is reduced, and the process returns to the step of calculating the stress information, strain information, and anchor force information of each second concentric ring step by step from the first second concentric ring inward; Jordi The radial stress of the second concentric ring is equal to the support reaction force. 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.

5. The tunnel face anchor parameter simulation calculation method based on extrusion deformation according to claim 4 is characterized in that: The method starts from the first second concentric ring and gradually calculates the stress information, strain information and anchor force information of each second concentric ring inward, including: The tangential stress of the first second concentric ring is: in, represents the radial stress of the first second concentric ring, represents the tangential stress 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 caused by elasticity, represents the elastic modulus of the first second concentric ring, represents the Poisson's ratio of the first second concentric ring; The anchor force information of the first second concentric ring is: in, Indicates the anchor force information of the first second concentric ring. represents the cross-sectional area of ​​the anchor rod, represents the elastic modulus of the anchor, represents the total radial strain of the first second concentric ring, ; No. The stress information of the second concentric ring is: in, in, represents the width of the second concentric ring, , represents the excavation radius of the tunnel face, Indicates the anchor length, Indicates The average radius of the second concentric rings, Indicates The strain increment of the second concentric ring is, Indicates the vertical spacing of anchor rods, Indicates the horizontal spacing of anchor rods; like Greater than , ; like Less than , = , ; Update by the first formula and ;in, Indicates the reference value of radial stress; 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, 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; No. The anchor force information of the second concentric ring is: ； The first formula is: ; Indicates The average tangential stress of the second concentric ring is, Indicates The average radial stress of the second concentric ring.

6. The tunnel face anchor parameter simulation calculation method based on extrusion deformation according to claim 1 is characterized in that: The calculating of the second extrusion deformation based on the tunnel site geological data, tunnel section data and anchor rod parameters includes: If the bolt length is less than the radius of the plastic zone, the 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 information, strain information and anchor force information of each second concentric ring are calculated step by step inward based on the set radius of the second plastic zone until the first second concentric ring is reached. A second concentric ring; Jordi The radial stress of the second concentric ring is less than the support reaction force, the radius of the second plastic zone is increased, and the process returns to the step of calculating the stress information, strain information, and anchor force information of each second concentric ring gradually inward starting from the first second concentric ring; Jordi The radial stress of the second concentric ring is greater than the support reaction force, the radius of the second plastic zone is reduced, and the process returns to the step of calculating the stress information, strain information, and anchor force information of each second concentric ring step by step from the first second concentric ring inward; Jordi The radial stress of the second concentric ring is equal to the support reaction force. 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.

7. The method for simulating and calculating parameters of anchor bolts at a tunnel face based on extrusion deformation according to claim 6, characterized in that: The method starts from the first second concentric ring and gradually calculates the stress information, strain information and anchor force information of each second concentric ring inward, including: 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 caused by elasticity, represents the elastic modulus of the first second concentric ring, represents the Poisson's ratio of the first second concentric ring; The anchor force information of the first second concentric ring is: ； Jordi The inner diameter of the second concentric ring is larger than the length of the anchor rod. 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 second plastic radius, Indicates The average radius of the second concentric rings, For the The friction angle of the second concentric ring is, Indicates The cohesion of a second concentric ring; Jordi The inner diameter of the second concentric ring is smaller than the length of the anchor rod. The inner diameter of the second concentric ring is updated to the anchor length, The stress information of the second concentric ring is: represents the cross-sectional area of ​​the anchor rod, represents the elastic modulus of the anchor, Indicates the vertical spacing of anchor rods, Indicates the horizontal spacing of anchor rods; 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, 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.

8. A tunnel face anchor parameter simulation calculation device based on extrusion deformation, characterized in that: include: Data acquisition module, used to obtain tunnel site geological data, tunnel section data and anchor parameters; A parameter calculation module is used to perform the anchor parameter update operation multiple times until the first stop condition is met, and obtain the tunnel face anchor parameter simulation calculation result; The anchor parameter update operation includes: The second extrusion deformation is calculated based on the geological data of the tunnel site, the tunnel section data and the anchor parameters; the second extrusion deformation is the extrusion deformation of the tunnel face under the condition of anchor reinforcement of the tunnel face; Calculating a second safety factor according to the second extrusion deformation; the second safety factor is used to characterize the stability of the tunnel face under the condition of face anchor reinforcement; If the second safety factor is less than or equal to a set threshold, adjusting the anchor parameters; The first stop condition is: the second safety factor is greater than a set threshold; Before performing anchor parameter update operations multiple times, the parameter calculation module is specifically 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, executing the anchor parameter update operation; The first extrusion deformation amount is calculated based on the geological data of the tunnel site and the tunnel section data, including: Divide the tunnel surrounding rock into a plurality of first concentric circular rings from the first elastic-plastic boundary inwards to the tunnel excavation radius; 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 second 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 is 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 second 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