A method for generating a prestressed anchoring member for end anchoring and lengthening anchoring

Abstract

The application discloses a prestressed anchoring component generation method for end anchoring and lengthening anchoring, and belongs to the field of numerical calculation of roadway support. A grid model before roadway excavation is created by using FLAC3D; a constitutive model and material parameters are set for the grid model; boundary conditions are set for the grid model; gravity acceleration is set for the grid model; initial conditions are set for the grid model; a large deformation calculation mode is set; automatic solution is carried out on the grid model; part of the roadway in the grid model is excavated; anchoring component parameters are valued; the first row of anchoring components in the roadway is traversed; connection operation is carried out on the anchor unit nodes with ID 1; five material parameters are set for all the anchor units; two material parameters are set for the anchor units with ID 1 and ID 2; and automatic solution is carried out on the grid model again. The application has important significance for improving anchoring component generation efficiency and revealing the mechanical action principle between the anchoring component and the surrounding rock of the roadway.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of numerical calculation of roadway support, and particularly relates to a prestressed anchoring component generation method for end anchoring and lengthened anchoring. BACKGROUND

[0002] The anchoring component is commonly used in the mine roadway surrounding rock control and plays a crucial role in ensuring the stability of the mine roadway and chamber surrounding rock. According to the geometric shape, the commonly used anchoring components in mines are divided into two types of anchor rods and anchor cables. In order to effectively reveal the anchoring mechanism of the anchoring component, researchers and engineers often simulate the anchoring component by numerical calculation and analyze the reinforcement effect of the anchoring component on the roadway surrounding rock in engineering practice.

[0003] At present, in the analysis of mining engineering geotechnical, the commonly used element that can simulate the anchoring component is the anchor unit developed by Itasca Company. The anchor unit has simple structure and corresponding material parameters are designed, which can not only simulate the anchoring component of laboratory scale, but also simulate the anchoring component of engineering scale, and is widely used in the analysis of roadway support. However, there are still some problems in simulating the prestressed anchoring component of end anchoring and lengthened anchoring by using the anchor unit:

[0004] (1) It is difficult to strictly match the length of the anchoring segment with the length of the anchoring segment required by the user when simulating the anchoring component

[0005] The anchoring component of end anchoring and lengthened anchoring has anchoring segment and free segment. Ignoring the exposed segment outside the borehole, the sum of the length of the anchoring segment and the length of the free segment is the full length of the anchoring component in the borehole. When the user simulates the anchoring component by using the anchor unit, the anchor unit is often generated between the end of the anchoring segment and the end of the free segment, and then the anchor unit is used to simulate the anchoring component. In order to simulate the anchoring segment and the free segment respectively, the user often needs to divide the anchor unit into multiple segments, and uses the nodes between the segments to simulate the connection points between the anchoring segment and the free segment. However, due to the matching relationship between the length of the anchoring component and the number of segments, after the anchor unit is divided into multiple segments, the nodes between the segments are difficult to coincide with the connection points between the anchoring segment and the free segment, which leads to the difficulty in strictly matching the length of the anchoring segment with the length of the anchoring segment required by the user during simulation.

[0006] (2) The prestress is difficult to be accurately installed when the extension direction of the anchoring component is not parallel to the coordinate axis in the coordinate system

[0007] When modeling a roadway, for the convenience of excavating the roadway, the direction of the roadway is often parallel to a coordinate axis in a coordinate system, such as the X axis, the Y axis or the Z axis. However, the anchoring member is not necessarily perpendicular to the rock wall of the roadway when installed. For example, at the position of the corner between the roof of the roadway and the side of the roadway, the anchoring member is often installed in the rock wall of the roadway in an inclined manner at a certain angle. This causes part of the anchoring member to be non-parallel to any coordinate axis. This increases the difficulty of installing the prestress of the anchoring member. When the anchoring member is parallel to the direction of the coordinate axis, the user only needs to apply the prestress to the anchoring member along the direction of the coordinate axis. However, when the anchoring member is non-parallel to the direction of the coordinate axis, it is difficult to install the prestress along the extension direction of the anchoring member.

[0008] (3) The code is relatively redundant when installing the anchoring members in the roadway in rows

[0009] In order to ensure the stability of the surrounding rock of the roadway, the anchoring members in the roadway are often arranged in rows along the extension direction of the roadway. Under general conditions, the row spacing of the anchoring members along the extension direction of the roadway is a fixed value. When installing the anchoring members in rows by using the anchor unit, the user often needs to use a large number of commands to generate and install each row of anchoring members, which causes the code to be relatively redundant and causes inconvenience to subsequent calculations.

[0010] In view of the above problems, the present application provides a prestressed anchoring member generation method for end anchoring and lengthening anchoring, which aims to provide a quick and efficient generation method for anchoring members for end anchoring and lengthening anchoring, and has important significance for improving the generation efficiency of anchoring members and revealing the mechanical action principle between the anchoring members and the surrounding rock of the roadway. SUMMARY

[0011] The purpose of the present application is to provide a prestressed anchoring member generation method for end anchoring and lengthening anchoring. This method overcomes the defect that the user cannot accurately and efficiently simulate the anchoring member by relying on the traditional anchor unit, and can quickly and efficiently simulate the anchoring member in the roadway support.

[0012] The application adopts the technical scheme and provides a prestressed anchoring component generation method for end anchoring and lengthened anchoring, comprising the following steps: creating a grid model before roadway excavation by using FLAC3D; setting a constitutive model and material parameters for the grid model; setting boundary conditions for the grid model; setting gravity acceleration for the grid model; setting initial conditions for the grid model; setting a large deformation calculation mode; automatically solving the grid model; excavating part of the roadway in the grid model; assigning values to anchoring component parameters; traversing the first row of anchoring components in the roadway; connecting the anchor unit nodes with ID 1; setting five material parameters, namely young, cross-sectional-area, yield-tension, grout-stiffness and grout-perimeter, for all anchor units; setting two material parameters, namely grout-cohesion and grout-friction, for the anchor units with ID 1; setting two material parameters, namely grout-cohesion and grout-friction, for the anchor units with ID 2, and setting both to zero; and automatically solving the grid model again.

[0013] Further description of the above technical scheme is as follows:

[0014] The extension direction of the roadway in the grid model needs to be parallel to the direction of any coordinate axis.

[0015] Further description of the above technical scheme is as follows:

[0016] The large deformation calculation mode is false.

[0017] Further description of the above technical scheme is as follows:

[0018] The anchoring component parameters include five parameters, namely parameter , parameter , parameter , parameter and parameter ; parameter represents the anchoring segment length; parameter represents the anchoring component row distance; parameter represents the anchoring component row number; parameter represents the prestress size; and parameter represents the anchoring component row arrangement direction flag.

[0019] Further description of the above technical scheme is as follows:

[0020] When the anchoring component row arrangement direction is parallel to the X axis, the parameter is set to 1; when the arrangement direction of the anchor members in the row is parallel to the Y axis, the parameter is set to 2; when the arrangement direction of the anchor members in the row is parallel to the Z axis, the parameter is set to 3.

[0021] As a further description of the above technical solution:

[0022] When traversing the first row of anchor members in the roadway, after traversing any anchor member, the inner end head coordinates of the anchor member are first extracted and assigned to the parameters in vector form, at this time , , respectively represent the X coordinate, Y coordinate and Z coordinate of the inner end head of the anchor member; then the outer end head coordinates of the anchor member are extracted and assigned to the parameters in vector form, at this time , , respectively represent the X coordinate, Y coordinate and Z coordinate of the outer end head of the anchor member; finally, the FISH function cable_install is called.

[0023] As a further description of the above technical solution:

[0024] The logical structure of the FISH function cable_install is as follows: the parameter is set to 1; the parameter is set to 1; the parameter is set to 1; the parameter is set to 1; a first formula is executed; a second formula is executed; a third formula is executed; a fourth formula is executed; a loop operation is executed.

[0025] As a further description of the above technical solution:

[0026] The first formula is , in which is a parameter; the second formula is , in which is a parameter; the third formula is , in which is a parameter; the fourth formula is , in which is a parameter.

[0027] As a further description of the above technical solution:

[0028] The logical structure of the loop operation is as follows: it is judged whether the parameter is less than or equal to the parameter If not, exit the loop; if yes, execute the fifth, sixth, seventh, eighth, ninth, tenth, eleventh, twelfth, and thirteenth formulas in sequence; , , Starting from, , , Create an anchor cell for the endpoint and set its ID to 1; , , Starting from, , , Create an anchor element for the endpoint and set its ID to 2; execute formulas fourteen, fifteen, and sixteen; for coordinates... , , The anchor element node is applied along the X-axis direction The tension; for coordinates as , , The anchor unit node is applied along the Y-axis direction The tension; for coordinates as , , The anchor element node is applied along the Z-axis direction The tensile force; determine the parameters If the value is 1, exit the judgment; if it is, return the parameter. Increment the value by 1; check the parameter If the value is 2, exit the judgment; if it is, return the parameter. Increment the value by 1; check the parameter If the value is 3, exit the judgment; if it is, return the parameter. Increment the value by 1; parameter Increment the value by 1 and continue the loop.

[0029] As a further description of the above technical solution:

[0030] The fifth formula is: , in the formula The parameter is; the sixth formula is , in the formula The parameter is; the seventh formula is , in the formula The parameter is; the eighth formula is , in the formula The parameter is; the ninth formula is , in the formula The parameter is; the tenth formula is , in the formula is a parameter; the eleventh formula is , in the formula is a parameter; the twelfth formula is , in the formula is a parameter; the thirteenth formula is , in the formula is a parameter; the fourteenth formula is , in the formula is a parameter; the fifteenth formula is , in the formula is a parameter; the sixteenth formula is , in the formula is a parameter.

[0031] Advantages

[0032] The advantages of the present application mainly include:

[0033] (1) The present application designs a FISH function cable_install. When the user uses the FISH function, the FISH function will automatically calculate the connection point coordinates of the anchoring segments and free segments of all anchoring members according to the anchoring segment length provided by the user, and use anchor units with different IDs to simulate the anchoring segments and free segments respectively, so as to ensure that the length of the simulated anchoring segments in the anchoring member is strictly matched with the anchoring segment length required by the user. This solves the problem that in the conventional simulation, the nodes between the segments are difficult to coincide with the connection points between the anchoring segments and free segments when the anchoring member is divided into segments, and further solves the problem that in the conventional simulation, the length of the simulated anchoring segments is difficult to strictly match with the anchoring segment length required by the user. In addition, if the connection point coordinates of the anchoring segments and free segments are manually calculated and the anchoring segments and free segments are simulated by segments, a large amount of manual calculation and operation is required and a lot of time is consumed. By using the FISH function, the user does not need to manually calculate, and after providing the required anchoring segment length, the FISH function can quickly traverse all anchoring members, automatically calculate the connection point coordinates of the anchoring segments and free segments in all anchoring members, and then automatically simulate the anchoring segments and free segments by using anchor units, greatly simplifying the anchoring member generation process.

[0034] (2) When the anchoring member is in any direction, the FISH function cable_install designed by the application can calculate the extension direction components of the anchoring member extension direction along the X-axis direction, Y-axis direction and Z-axis direction, respectively. According to the calculated extension direction components, the FISH function will calculate the prestress components of each anchoring member total prestress in the X-axis direction, Y-axis direction and Z-axis direction, respectively. Then, the FISH function will automatically capture each anchoring member outer end node and install the calculated prestress components in each anchoring member outer end node. In actual calculation, the prestress components of each anchoring member outer end node along the X-axis direction, Y-axis direction and Z-axis direction will be combined into a prestress, and the prestress action direction is parallel to the anchoring member extension direction and points to the tunnel direction. This solves the problem that the prestress is difficult to be parallel to the anchoring member extension direction when the prestress is applied to the inclined anchoring member.

[0035] (3) The application sets the parameter in the FISH function cable_install to represent the arrangement direction of the anchoring members. When the user provides the anchoring member row distance and the anchoring member arrangement direction, the FISH function can generate multiple rows of anchoring members along the X-axis direction, Y-axis direction or Z-axis direction according to the anchoring member row distance and the anchoring member arrangement direction. When the number of anchoring member rows is equal to the number of rows specified by the user, the FISH function will stop generating anchoring members. Therefore, by using the FISH function, the redundant process of manually calculating the coordinates of each row of anchoring members can be greatly simplified, and multiple rows of anchoring members can be quickly and efficiently generated along the coordinate axis direction. BRIEF DESCRIPTION OF DRAWINGS

[0036] The accompanying drawings, which form a part of this application, are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification. In the drawings:

[0037] Figure 1 is an operation flow chart of the prestressed anchoring member generation method for end anchoring and lengthening anchoring according to the application.

[0038] Figure 2 is a logic structure diagram of the FISH function cable_install according to the application.

[0039] Figure 3 is a loop operation logic structure diagram according to the application.

[0040] Figure 4 is a prestress action direction display diagram after the anchoring member is applied with prestress by using the application. DETAILED DESCRIPTION

[0041] As Figure 1As shown, the present application provides a prestressed anchoring member generation method for end anchoring and lengthening anchoring, comprising: creating a grid model before roadway excavation by FLAC3D; setting a constitutive model and material parameters for the grid model; setting boundary conditions for the grid model; setting gravity acceleration for the grid model; setting initial conditions for the grid model; setting a large deformation calculation mode; automatically solving the grid model; excavating part of the roadway in the grid model; assigning values to anchoring member parameters; traversing the first row of anchoring members in the roadway; connecting the anchor unit nodes with ID 1; setting five material parameters of young, cross-sectional-area, yield-tension, grout-stiffness and grout-perimeter for all anchor units; setting two material parameters of grout-cohesion and grout-friction for the anchor unit with ID 1; setting two material parameters of grout-cohesion and grout-friction for the anchor unit with ID 2, and setting both to zero; automatically solving the grid model again.

[0042] In a specific embodiment:

[0043] The extension direction of the roadway in the grid model needs to be parallel to the direction of any coordinate axis.

[0044] In a specific embodiment:

[0045] The large deformation calculation mode is false.

[0046] In a specific embodiment:

[0047] The anchoring member parameters include five parameters, which are parameters , parameters , parameters , parameters and parameters ; parameter represents the anchoring segment length; parameter represents the anchoring member row distance; parameter represents the anchoring member row number; parameter represents the prestress size; and parameter represents the anchoring member row arrangement direction flag.

[0048] In a specific embodiment:

[0049] When the anchoring member row arrangement direction is parallel to the X axis, the parameter is set to 1; and when the anchoring member row arrangement direction is parallel to the Y axis, the parameter is set to 2; when the arrangement direction of the anchor members in the row is parallel to the Z axis, the parameter is set to 3.

[0050] In a specific embodiment:

[0051] When traversing the first row of anchor members in the roadway, after traversing any anchor member, the inner end head coordinates of the anchor member are first extracted and assigned to the parameters in vector form, at this time , , respectively represent the X coordinate, Y coordinate and Z coordinate of the inner end head of the anchor member; then the outer end head coordinates of the anchor member are extracted and assigned to the parameters in vector form, at this time , , respectively represent the X coordinate, Y coordinate and Z coordinate of the outer end head of the anchor member; finally, the FISH function cable_install is called.

[0052] In a specific embodiment:

[0053] The logical structure of the FISH function cable_install is as follows: the parameter is set to 1; the parameter is set to 1; the parameter is set to 1; the parameter is set to 1; a first formula is executed; a second formula is executed; a third formula is executed; a fourth formula is executed; a loop operation is executed.

[0054] In a specific embodiment:

[0055] The first formula is , in which is a parameter; the second formula is , in which is a parameter; the third formula is , in which is a parameter; the fourth formula is , in which is a parameter.

[0056] In a specific embodiment:

[0057] The logical structure of the loop operation is as follows: it is determined whether the parameter is less than or equal to the parameter , if not, the loop is exited; if yes, a fifth formula, a sixth formula, a seventh formula, an eighth formula, a ninth formula, a tenth formula, an eleventh formula, a twelfth formula and a thirteenth formula are executed in sequence; and , 、 , 、 、 , 、 、 , 、 、 , 、 、 , , 、 、 , , 、 、 , , , , , , , , ,

[0058] In a specific embodiment:

[0059] the fifth formula is , wherein is a parameter; the sixth formula is , wherein is a parameter; the seventh formula is , wherein is a parameter; the eighth formula is , wherein is a parameter; the ninth formula is , wherein is a parameter; the tenth formula is , wherein is a parameter; the eleventh formula is , wherein is a parameter; the twelfth formula is , wherein is a parameter; the thirteenth formula is , in the formula is a parameter; the fourteenth formula is , in the formula is a parameter; the fifteenth formula is , in the formula is a parameter; the sixteenth formula is , in the formula is a parameter.

[0060] To verify the effectiveness of the present application, a group of calculation cases are set and analyzed. A hexahedral mesh model is created in FLAC3D, and the sizes of the model along the X-axis direction, the Y-axis direction and the Z-axis direction are 30 m, 5 m and 30 m respectively. The grid model is set to the Mohr-Coulomb constitutive model, and the material parameters include Young's modulus, Poisson's ratio, cohesion, internal friction angle, tensile strength and density. The boundary conditions of the grid model are set, specifically, a compressive stress of 15 MPa is applied to the top boundary, and the other boundaries are supported by a roller. The gravity acceleration of the grid model is set, specifically, the gravity acceleration along the Z-axis direction is vertically downward and the size is 10 m / s 2 . The initial conditions of the grid model are set, specifically, the initial stress matches the top boundary compressive stress and the gravity acceleration. The large deformation calculation mode is set to false. The grid model is automatically solved. The part of the roadway in the grid model is excavated, specifically, the grid with X coordinate falling into the interval [12.5 17.5] and Z coordinate falling into the interval [13 17].

[0061] The anchor member parameters are assigned values, specifically, the parameter is 1.2 m, the parameter is 0.9 m, the parameter is 6, and the parameter is 150 kN. In this embodiment, the roadway extension direction is parallel to the Y-axis direction, so the anchor member row arrangement direction is parallel to the Y-axis, and therefore the parameter is set to 2.

[0062] The first row of anchor members in the roadway is traversed, and when any anchor member is traversed, the inner end head coordinates of the anchor member are first extracted and assigned to the parameter in the form of a vector, at this time , , represent the X coordinate, the Y coordinate and the Z coordinate of the inner end head of the anchor member respectively; then the outer end head coordinates of the anchor member are extracted and assigned to the parameter in the form of a vector, at this time , , These represent the X, Y, and Z coordinates of the outer end of the anchoring component, respectively; finally, the FISH function cable_install is called.

[0063] The logical structure of the FISH function cable_install is as follows: Figure 2 As shown, set the parameters Set the parameter to 1. Set the parameter to 1. Set the parameter to 1. Set the value to 1; execute the first formula; execute the second formula; execute the third formula; execute the fourth formula; execute the loop operation.

[0064] The first formula is , in the formula The parameter is; the second formula is , in the formula The parameter is; the third formula is , in the formula The parameter is; the fourth formula is , in the formula For parameters.

[0065] The logical structure of the loop operation is as follows: Figure 3 As shown, the parameters are determined. Is it less than or equal to the parameter? If not, exit the loop; if yes, execute the fifth, sixth, seventh, eighth, ninth, tenth, eleventh, twelfth, and thirteenth formulas in sequence; , , Starting from, , , Create an anchor cell for the endpoint and set its ID to 1; , , Starting from, , , Create an anchor element for the endpoint and set its ID to 2; execute formulas fourteen, fifteen, and sixteen; for coordinates... , , The anchor unit node is applied along the X-axis direction The tension; for coordinates as , , The anchor unit node is applied along the Y-axis direction The tension; for coordinates as , , The anchor unit nodes of the first group apply pulling force along the Z-axis direction ; determine whether the parameter is equal to 1, if not, exit the determination, if yes, increase the parameter value by 1 ; determine whether the parameter is equal to 2, if not, exit the determination, if yes, increase the parameter value by 1 ; determine whether the parameter is equal to 3, if not, exit the determination, if yes, increase the parameter value by 1 ; determine whether the parameter is equal to 4, if not, exit the determination, if yes, increase the parameter value by 1 ; determine whether the parameter is equal to 5, if not, exit the determination, if yes, increase the parameter value by 1 ; determine whether the parameter is equal to 6, if not, exit the determination, if yes, increase the parameter value by 1 ; determine whether the parameter is equal to 7, if not, exit the determination, if yes, increase the parameter value by 1 ; increase the parameter value by 1, and continue the loop.

[0066] The fifth formula is , wherein is a parameter; the sixth formula is , wherein is a parameter; the seventh formula is , wherein is a parameter; the eighth formula is , wherein is a parameter; the ninth formula is , wherein is a parameter; the tenth formula is , wherein is a parameter; the eleventh formula is , wherein is a parameter; the twelfth formula is , wherein is a parameter; the thirteenth formula is , wherein is a parameter; the fourteenth formula is , wherein is a parameter; the fifteenth formula is , wherein is a parameter; the sixteenth formula is , wherein is a parameter.

[0067] Connect the anchor unit nodes with ID 1; set young of all anchor units to 200 GPa, and cross-sectional-area to 380.133 mm 2, yield-tension of 186 kN, grout-stiffness of 100 MPa, grout-perimeter of 69.115 mm; grout-cohesion of 100 kN / m and grout-friction of 32° are set for the anchor unit with ID of 1 ; grout-cohesion of 0 and grout-friction of 0 are set for the anchor unit with ID of 2. The grid model is solved again automatically.

[0068] After the calculation, the coordinates of the connection points of the anchoring segments and the free segments of all the anchoring members are output. In order to verify the accuracy of the calculation results of the present application, all the anchoring members are plotted in AutoCAD according to the inner end coordinates and the outer end coordinates. According to the anchoring segment length of 1.2 m input by the user in the present embodiment, all the connection points of the anchoring segments and the free segments of the anchoring members are circled. The coordinates of all the connection points of the anchoring segments and the free segments of the anchoring members calculated by the present application and the coordinates of all the connection points of the anchoring segments and the free segments of the anchoring members measured by AutoCAD are compared, and it is found that they are consistent, which indicates the accuracy of the coordinates of the connection points of the anchoring segments and the free segments of the anchoring members calculated by the present application. This solves the problem that the simulated anchoring segment length is difficult to strictly match the anchoring segment length required by the user in the conventional simulation. In addition, if the user uses FLAC3D in combination with AutoCAD, measures the connection point coordinates of the anchoring segments and the free segments by AutoCAD, and manually segments and generates the anchoring members of the anchoring segments and the free segments in FLAC3D, a large amount of point drawing, point taking and anchoring member generation operations will be performed, which is extremely redundant and tedious. By using the present application, after the user inputs the anchoring segment length, the present application automatically calculates the connection point coordinates of the anchoring segments and the free segments, and the user does not need to manually measure a large number of connection point coordinates, which greatly simplifies the operation and improves the efficiency. The inner axial forces of the anchoring members are displayed, and it is found that the inner axial forces of all the anchoring members gradually increase from the inner end to the direction of the roadway in the anchoring segments, and the axial force is the largest at the connection point of the anchoring segment and the free segment. The axial force of the anchoring member in the free segment remains constant. In addition, the surface bonding stress of the anchoring member is displayed, and it is found that there is a bonding stress on the surface of the anchoring member in the anchoring segment, and the surface bonding stress of the anchoring member in the free segment is zero. The above-mentioned axial force law of the anchoring member and the surface bonding stress law of the anchoring member are consistent with the actual working conditions, which verifies the effectiveness of the anchoring members generated by the present application.

[0069] After the calculation, the prestress action directions of all the anchoring members are displayed as shown in Fig. 8. Figure 4As shown, the prestressing of the anchoring members arranged horizontally is all along the horizontal direction and points towards the excavation direction at the center of the roadway; the prestressing of the anchoring members arranged vertically is all along the vertical direction and points towards the excavation direction at the center of the roadway; and the prestressing of the anchoring members arranged inclinedly is all along the extension direction of the anchoring member and points towards the excavation direction at the center of the roadway. Therefore, this invention overcomes the problem of difficulty in installing prestressing parallel to the inclined direction of the anchoring member when using traditional simulation methods.

[0070] Finally, in this embodiment, the tunnel extension direction is parallel to the Y-axis direction. To support the tunnel, the anchoring components are arranged in rows along the Y-axis direction within the tunnel. In this embodiment, the user-defined row spacing of the anchoring components is 0.9 m, and the row arrangement direction indicator is 2. Therefore, this invention automatically installs the first row of anchoring components at equal intervals along the tunnel extension direction, i.e., the Y-axis direction, based on the user-input row spacing. In this embodiment, the user-defined number of anchoring component rows is 6. Therefore, when six rows of anchoring components are generated, this invention stops generating anchoring components. After the calculation is completed, the row spacing of the anchoring components (0.9 m) and the number of rows (6) are checked and found to be the same as the user-defined values, verifying the reliability of this invention. Compared to users installing anchoring components row by row through commands, this invention automatically installs anchoring components that meet user requirements by relying on the user-defined row spacing, number of rows, and row arrangement direction of the anchoring components. This greatly simplifies the redundant process of manually calculating the coordinates of each row of anchoring components and manually installing them, and significantly improves the installation efficiency of anchoring components.

[0071] This invention is not limited to the preferred embodiments described above. Anyone can derive other forms of products under the guidance of this invention. However, regardless of any changes made in their shape or structure, any technical solution that is the same as or similar to this application falls within the protection scope of this invention.

Claims

1. A method of generating a prestressed anchoring member for end anchoring and lengthening anchoring, characterized by, The method comprises the following steps: creating a grid model before roadway excavation by using FLAC3D; setting constitutive model and material parameters for the grid model; setting boundary conditions for the grid model; setting gravity acceleration for the grid model; setting initial conditions for the grid model; setting large deformation calculation mode; and automatically solving the grid model; The anchoring member parameters are assigned values, the anchoring member parameters including five parameters, respectively being a parameter , a parameter , a parameter , a parameter , and a parameter ; a parameter represents an anchoring segment length; a parameter represents an anchoring member row spacing; a parameter represents an anchoring member row number; a parameter represents a prestress size; and a parameter represents an anchoring member row arrangement direction flag. The first row of anchoring members in the roadway is traversed; when any anchoring member is traversed, the inner end head coordinates of the anchoring member are first extracted and are assigned to parameters in vector form At this time , , respectively represent the X coordinate, the Y coordinate and the Z coordinate of the inner end head of the anchoring member; then the outer end head coordinates of the anchoring member are extracted and are assigned to parameters in vector form At this time , , respectively represent the X coordinate, the Y coordinate and the Z coordinate of the outer end head of the anchoring member; finally, the FISH function cable_install is called. The FISH function cable_install logic structure is as follows: set parameter to 1; set parameter to 1; set parameter to 1; set parameter to 1; execute first formula; execute second formula; execute third formula; execute fourth formula; excavating part of the roadway in the grid model; The first formula is , where is a parameter; the second formula is , where is a parameter; the third formula is , where is a parameter; and the fourth formula is , where is a parameter. The loop operation logic structure is as follows: judging whether the parameter is less than or equal to the parameter , and if not, exiting the loop; If yes, execute the fifth formula, the sixth formula, the seventh formula, the eighth formula, the ninth formula, the tenth formula, the eleventh formula, the twelfth formula and the thirteenth formula in turn; create anchor units with the starting point 、 、 and the end point 、 、 and set the ID as 1; create anchor units with the starting point 、 、 and the end point 、 、 and set the ID as 2; execute the fourteenth formula, the fifteenth formula and the sixteenth formula; apply a pulling force of to the anchor unit node with the coordinates 、 、 in the X-axis direction; apply a pulling force of to the anchor unit node with the coordinates 、 、 in the Y-axis direction; apply a pulling force of to the anchor unit node with the coordinates 、 、 in the Z-axis direction; judge whether the parameter is equal to 1, if not, exit the judgment, if yes, add 1 to the parameter ; judge whether the parameter is equal to 2, if not, exit the judgment, if yes, add 1 to the parameter ; judge whether the parameter is equal to 3, if not, exit the judgment, if yes, add 1 to the parameter ; add 1 to the parameter and continue the loop; The fifth formula is , in which is a parameter; the sixth formula is , in which is a parameter; the seventh formula is , in which is a parameter; the eighth formula is , in which is a parameter; the ninth formula is , in which is a parameter; the tenth formula is , in which is a parameter; the eleventh formula is , in which is a parameter; the twelfth formula is , in which is a parameter; the thirteenth formula is , in which is a parameter; the fourteenth formula is , in which is a parameter; the fifteenth formula is , in which is a parameter; the sixteenth formula is , in which is a parameter; performing loop operation; 2. The method of claim 1, wherein the method is characterized by: connecting anchor unit nodes with ID 1; setting five material parameters, i.e. young, cross-sectional-area, yield-tension, grout-stiffness and grout-perimeter, for all anchor units; setting two material parameters, i.e. grout-cohesion and grout-friction, for anchor units with ID 1; setting two material parameters, i.e. grout-cohesion and grout-friction, for anchor units with ID 2, and setting both to zero; and automatically solving the grid model again.

3. The method of claim 1, wherein the method is characterized by: The roadway in the grid model extends in parallel with any coordinate axis. The large deformation calculation mode is false.

4. The method of claim 1, wherein the method is characterized by: when the anchoring members are arranged in a row parallel to the X axis, the parameter is set to 1 ; when the anchoring members are arranged in a row parallel to the Y axis, the parameter is set to 2; when the anchoring members are arranged in a row parallel to the Z axis, the parameter is set to 3.

Images