A method for selecting and predicting tunnel segment types

Abstract

This invention discloses a method for selecting and predicting tunnel segment types, including: S1, determining the assembly rules for each ring segment based on the construction site conditions; S2, dividing the tunnel design axis into multiple straight and curved segments, and calculating the alignment parameters of each straight and curved segment in order of mileage; S3, obtaining the alignment parameters corresponding to the current location of the tunnel boring machine on the tunnel design axis; S4, calculating all possible assembly points K for the current ring segment. i S5. Calculate the weighted average of the current ring segment's axis rotation angle, shield tail clearance, propulsion cylinder stroke, and tunnel boring machine trend to determine the K values ​​for each optional assembly point in the current ring segment. i Weighted score J i Select the weighted score J i The assembly point K corresponding to the maximum value i As the current ring segment assembly point; S6, repeat steps S4 and S5 to predict the next ring segment type and assembly point. This invention can realize the reasonable selection of the current ring segment for standard rings and left and right turning ring segments, as well as the prediction of the next ring segment selection.

Description

Technical Field

[0001] This invention relates to the field of shield tunnel engineering. More specifically, this invention relates to a method for selecting and predicting shield tunnel segments. Background Technology

[0002] Currently, tunnel lining segments are mainly divided into general wedge-shaped ring segments and standard ring + turning ring segments. Existing technologies use general-purpose segment selection methods. Currently, the selection of standard ring + turning ring segments relies primarily on manual methods. However, manual selection is highly subjective and limited by the experience and understanding of operators. During construction, unreasonable assembly point selection frequently occurs, resulting in mismatches between segment locations and the tunnel's design axis and shield posture. This leads to serious consequences such as shield posture exceeding limits, segment misalignment or cracking, and tail-end leakage.

[0003] On the other hand, unlike general-purpose tunnel segments, after the selection of the current ring segment is completed, the segment type for the next ring (standard ring, left-turn ring, right-turn ring) needs to be selected as soon as possible to facilitate the on-site organization of the hoisting of the next ring segment. Currently, the selection and prediction of the next ring segment also relies on manual methods. However, manual selection and prediction are highly subjective and constrained by the experience and knowledge limitations of the operators. During construction, situations often arise where the selection of the next ring segment type is unreasonable and does not match the tunnel design axis and shield posture, affecting the shield posture and the quality of the formed segments. Summary of the Invention

[0004] One object of the present invention is to solve at least the above-mentioned problems and to provide at least the advantages that will be described later.

[0005] To achieve these objectives and other advantages according to the present invention, a method for selecting and predicting tunnel segment types is provided, comprising the following steps:

[0006] S1. Determine the assembly rules for each ring segment based on the construction site conditions;

[0007] S2. Divide the tunnel design axis into multiple straight and curved segments, and calculate the alignment parameters of each straight and curved segment in order of mileage.

[0008] S3. Obtain the alignment parameters corresponding to the current tunnel design axis position of the tunnel boring machine based on the data calculated in step S2.

[0009] S4. Calculate all possible assembly points K for the current ring segments based on the assembly points and assembly rules of the previous ring segments. i ;

[0010] S5. Calculate the weighted average of four factors: the current ring segment's axial rotation angle, the shield tail clearance, the propulsion cylinder stroke, and the tunnel boring machine's trend, to determine the K values ​​for each selectable assembly point within the current ring segment.i Weighted score J i Select the weighted score J i The assembly point K corresponding to the maximum value i As the current assembly point for the ring tunnel segments;

[0011] S6. Repeat steps S4 and S5 to predict the next ring segment type and segment assembly location.

[0012] Preferably, step S1 specifically includes:

[0013] S1-1. Choose either staggered or continuous seam assembly;

[0014] S1-2, Select the assembly points that are prohibited from being selected.

[0015] Preferably, in step S2, the linear parameters of the straight line segment include the length L of the straight line segment. ZH The curve segment includes a circular curve and transition curves symmetrically positioned on both sides of the circular curve. The linear parameters of the curve segment include the curve length L and the transition curve length L. S The curve turning angle α and the radius R of the circular curve.

[0016] Preferably, step S4 specifically includes:

[0017] S4-1. Calculate the preliminary selectable assembly points for the current ring segments;

[0018] S4-2. Eliminate assembly points that do not conform to the current ring segment assembly rules to obtain all possible assembly points K for the current ring segment. i .

[0019] Preferably, the weighted score J in step S5 is... i Calculated using the following formula:

[0020] J i =V 1i I 1i +V 2i I 2i +V 3i I 3i +V 4i I 4i

[0021] In the formula I 1i I 2i I 3i I 4i The importance coefficients for four factors are: axis rotation angle, tail clearance, propulsion cylinder stroke, and tunnel boring machine trend; V 2i V 3i V 4i For each optional assembly point K iEffect scores of four factors including the axis rotation angle, tail gap, stroke of the propulsion cylinder, and the trend of the shield machine after simulated assembly.

[0022] Preferably, for each optional segment erection point K i The importance coefficient I of the corresponding axis rotation angle 1i And the effect score V 1i The calculation method is as follows:

[0023] S5-1-1. Calculate the segment rotation angle β i Of each optional segment erection point K i ;

[0024] S5-1-2. Calculate the cumulative rotation angle θ i Of each optional segment erection point K i = β1 - α1 + β2 - α2 + … + β i-1 - α i-1 + β i - α i ;

[0025] S5-1-3. Assign a value to the importance coefficient I of the axis rotation angle according to the magnitude of the cumulative rotation angle θ i-1 Of the previous ring. When θ 1i Is 0, I i-1 Is 1, and the value of |θ 1i | is proportional to the value of I i-1 ; 1i ;

[0026] S5-1-4. Assign a value to the effect score a1 according to the β - α value calculated after simulated assembly of each optional segment erection point K i Of the current ring. 0 < a1 < 1. When |β - α| = 0, a1 = 1, and the value of |β - α| is inversely proportional to the value of a1;

[0027] S5-1-5. Assign a value to the effect score a2 according to the cumulative rotation angle θ i Calculated after simulated assembly of each optional segment erection point K i Of the current ring. 0 < a2 < 1. When |θ i | = 0, a1 = 1, and the value of |θ i | is inversely proportional to the value of a2;

[0028] S5-1-6. The effect score V of the axis rotation angle after simulated assembly of each optional segment erection point K i = 0.4 × a1 + 0.6 × a2. 1i

[0029] Preferably, for each optional segment erection point K i The importance coefficient I of the corresponding tail gap2i and the effect score V 2i The calculation method is as follows:

[0030] S5-2-1. Based on the minimum shield tail gap value δ among the shield tail gap values ​​at the four positions (up, down, left, and right) before the current annular segment assembly,... min Importance coefficient I 2i Assignment, when δ min >δ lim At that time, I 2i =1; when δ min <δ lim At that time, I 2i >1, and I 2i Value and δ min The values ​​are inversely proportional; where δ lim This is the critical value for the shield tail gap;

[0031] S5-2-2, Calculate the selectable assembly points K for the current ring based on the shield tail gap values ​​at the four positions (up, down, left, and right) before the current ring segment assembly. i Simulate the shield tail gap values ​​at the four positions (top, bottom, left, and right) after assembly, and calculate based on the minimum shield tail gap value δ′. min Performance score V 2i Assignment, 0 <V 2i <1, and V 2i Value and δ′ min The value is directly proportional.

[0032] Preferably, each selectable assembly point K i The importance coefficient I of the corresponding hydraulic cylinder stroke 3i And the effect score V 3i The calculation method is as follows:

[0033] S5-3-1. Based on the cylinder stroke input values ​​at the four positions (up, down, left, and right) before the current ring segment assembly, calculate the stroke difference between the up and down cylinders and the stroke difference between the left and right cylinders, and select the larger cylinder stroke difference value △L. max Importance coefficient I 3i Assignment, when △L max >△L lim At that time, I 3i =1; when △L max <△L lim At that time, I 3i >1, and I 3i Value and △L max The value is directly proportional; where ΔL max This is the critical value for the difference in cylinder stroke.

[0034] S5-3-2. Based on the stroke difference between the upper and lower cylinders and the stroke difference between the left and right cylinders before the current ring segment assembly, calculate the possible assembly points K for each ring segment. iSimulated cylinder stroke difference variation after assembly; selectable assembly points K for each front ring. i The calculated value of the cylinder stroke difference after simulated assembly = cylinder stroke difference of the previous ring segment + cylinder stroke difference change value, and based on the maximum cylinder stroke difference value △L′. max Performance score V 3i Assignment, 0 <V 3i <1, and V 2i Value and δ′ min The value is inversely proportional.

[0035] Preferably, each selectable assembly point K i The importance coefficient I of the corresponding tunnel boring machine trend 4i And the effect score V 4i The calculation method is as follows:

[0036] S5-4-1. Calculate the tunnel boring machine trend importance coefficient I based on the current trend value of the tunnel boring machine before the assembly of the ring segments. 4i Assignment, I 4i The value is proportional to the absolute value of the shield machine trend value before the current ring segment assembly. When the current shield machine trend value before the current ring segment assembly is 0, I 4i =0;

[0037] S5-4-2. The effectiveness score of determining the shield tunneling trend based on the direction of the shield tunneling trend and the direction of the pipe ring capping block is V. 4i .

[0038] Preferably, step S6 specifically includes the following steps:

[0039] S6-1. Based on the selected assembly points and assembly rules of the current ring segment, predict the type of the next ring segment and calculate the possible assembly points K for the next ring segment. i+1 ;

[0040] S6-2. Calculate the possible assembly points K for each segment of the next ring. i+1 Simulated shield tail clearance and hydraulic cylinder stroke difference after assembly;

[0041] S6-3. Based on the selected assembly points of the current ring segment, determine the difference between the shield tail gap and the hydraulic cylinder stroke after the simulated assembly of the selected assembly points of the current ring segment, and calculate the change in shield tail gap and the change in hydraulic cylinder stroke when the next ring shield tunneling is completed.

[0042] S6-4. Calculate the predicted value of the shield tail gap and the predicted value of the hydraulic cylinder stroke difference when the next ring of shield tunneling is completed. Based on step 5, complete the prediction of the assembly point and segment type of the next ring of segments.

[0043] The present invention has at least the following beneficial effects:

[0044] The shield segment selection and prediction method provided by this invention comprehensively considers the influence of tunnel alignment factors, including the shield tail gap, cylinder stroke difference, and shield machine trend, to achieve reasonable selection of the current ring segment for standard ring + left and right turning ring segments. Based on the current ring segment selection method, and comprehensively considering the shield machine attitude, it calculates the shield tail gap and cylinder stroke after the next ring shield is excavated, to achieve reasonable prediction of the location and type of the next ring segment. This method can reduce or replace manual operations, and plays a role in standardizing the selection and prediction of standard ring + turning ring segment construction on site, optimizing on-site construction procedures, and improving shield attitude and the quality of formed segments.

[0045] Other advantages, objectives and features of the present invention will become apparent in part from the following description, and in part from those skilled in the art through study and practice of the invention. Attached Figure Description

[0046] Figure 1 This is a flowchart of the shield tunnel segment selection and prediction method described in this invention;

[0047] Figure 2 This is a structural schematic diagram of the tunnel design axis described in this invention;

[0048] Figure 3 This is a schematic diagram illustrating the evaluation of the shield tunneling trend effect score according to the present invention; Detailed Implementation

[0049] The present invention will now be described in further detail with reference to the accompanying drawings, so that those skilled in the art can implement it based on the description.

[0050] It should be noted that, unless otherwise specified, the experimental methods described in the following embodiments are all conventional methods, and the reagents and materials described are all commercially available unless otherwise specified. In the description of this invention, the terms "lateral", "longitudinal", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", and "outer" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are only for the convenience of describing this invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this invention.

[0051] like Figure 1 As shown, this invention provides a method for selecting and predicting tunnel lining segments, comprising the following steps:

[0052] S1. Determine the assembly rules for each ring segment based on the construction site conditions;

[0053] S2. Divide the tunnel design axis into multiple straight and curved segments, and calculate the alignment parameters of each straight and curved segment in order of mileage.

[0054] S3. Obtain the alignment parameters corresponding to the current tunnel design axis position of the tunnel boring machine based on the data calculated in step S2.

[0055] S4. Calculate all possible assembly points K for the current ring segments based on the assembly points and assembly rules of the previous ring segments. i ;

[0056] S5. Calculate the weighted average of four factors: the axis rotation angle of the current ring segment, the shield tail clearance, the propulsion cylinder stroke, and the tunnel boring machine trend, to determine the K of each optional assembly point in the current segment type. i Weighted score J i Select the weighted score J i The assembly point K corresponding to the maximum value i As the current assembly point for the ring tunnel segments;

[0057] S6. Repeat steps S4 and S5 to predict the next ring segment type and segment assembly location.

[0058] In this technical solution, the shield tail gap, cylinder stroke difference, and shield machine trend are taken into account to comprehensively consider the influence factors of tunnel alignment, so as to achieve reasonable selection of the current ring of standard ring + left and right turning ring segments. Based on the current ring segment selection method, the shield machine attitude is comprehensively considered to calculate the shield tail gap and cylinder stroke after the next ring shield is excavated, so as to achieve reasonable prediction of the location and type of the next ring segment.

[0059] Step S1 specifically includes:

[0060] S1-1. Choose either staggered or continuous seam assembly;

[0061] S1-2, Select the assembly points that are prohibited from being selected.

[0062] Based on the construction site conditions, first select staggered or continuous joint assembly, and then select prohibited assembly points, such as 1 point prohibited at the bottom, 3 points prohibited at the bottom, below the waist, and 1 point prohibited at the top.

[0063] In step S2, the tunnel design axis is divided into multiple straight segments and curved segments. The linear parameters of the straight segments include the length L of the straight segment. ZH The curve segment includes a circular curve and transition curves symmetrically positioned on both sides of the circular curve. The linear parameters of the curve segment include the curve length L and the transition curve length L. S The curve turning angle α and the radius R of the circular curve.

[0064] The steps for calculating the linear parameters of the curve segment are as follows: Refer to... Figure 2 First, calculate the length L of the circular curve. R =L-2L S Then calculate the circular curve angle α. R =L R ×180° / πR, then calculate the number of rings H of the circular curve segment. R =L R / W, where W is the segment ring width, then calculate the average rotation angle θ of the circular curve. R =α R / H R The number of rings H in the transition curve S =L S / W, transition curve angle α S =(α-α) R) / 2, average turning angle θ of the transition curve S =α S / H S .

[0065] When the starting or ending mileage of the tunnel design axis section is located within a transition curve or circular curve section, it is necessary to calculate the length L of the transition curve or circular curve within the mileage range of the section separately. S ′ or L R According to L S ′ or L R L S or L R Calculate the corresponding α based on the proportion. S ′ and α R ′.

[0066] In actual use, each segment is numbered from smallest to largest mileage to facilitate the storage of the corresponding alignment parameters.

[0067] In step S3, based on the current position of the tunnel boring machine on the tunnel design axis, a straight segment or a curved segment is selected. For the curved segment, a circular curve or a transition curve, as well as a left turn or a right turn, are further selected, and the corresponding line parameters are obtained. The straight segment corresponds to a standard ring segment, and the curved segment corresponds to a left-turn ring segment or a right-turn ring segment, thus completing the selection of the segment type.

[0068] Step S4 specifically includes:

[0069] S4-1. Calculate the initial selectable assembly points for the current ring segment K = M + 2 + 3N, where M is the assembly point of the previous ring segment, N = 0, 1, 2, ...; the maximum value of K does not exceed the maximum number of assembly points for a single segment.

[0070] S4-2. Eliminate assembly points that do not conform to the current ring segment assembly rules to obtain all possible assembly points K for the current ring segment. i .

[0071] The weighted score J mentioned in step S5 i Calculated using the following formula:

[0072] J i =V 1i I 1i +V 2i I 2i +V 3i I 3i +V 4i I 4i

[0073] In the formula I 1i I 2i I 3i I 4i The importance coefficients for four factors are: axis rotation angle, tail clearance, propulsion cylinder stroke, and tunnel boring machine trend; V 2i V 3i V 4i For each optional assembly point K i The simulation scores the effects of four factors: axis rotation angle, tail clearance, propulsion cylinder stroke, and tunnel boring machine trend after assembly.

[0074] Furthermore, each optional assembly point K i The importance coefficient I of the corresponding axis rotation angle 1i and the effect score V 1i The calculation method is as follows:

[0075] S5-1-1 Calculate the available assembly points K for each current ring segment. i Segment rotation angle β i ;

[0076] The optional assembly points K of the standard ring i All segment rotation angles are 0°. Calculation of segment rotation angles for left and right rotation rings.

[0077] Formulas are as follows: (1) and (2)

[0078]

[0079]

[0080] In the formula, s is the wedge shape of the tunnel segment, m is the location of the tunnel segment, n is the total number of locations of the tunnel segment, d is the outer diameter of the tunnel segment, the left turn angle is negative, and the right turn angle is positive; the curve angle α is the angle of the axis segment where the shield is currently located, which is the angle of the curve segment selected in step 3. If it is a straight segment, the curve angle is 0, the left turn angle is negative, and the right turn angle is positive.

[0081] S5-1-2 Calculate the available assembly points K for each ring segment.i The cumulative rotation angle θ i = β1 - α1 + β2 - α2 + … + β i-1 - α i-1 + β i - α i ;

[0082] β - α is the rotation angle change after the assembly of one ring of segment; the cumulative rotation angle θ is the cumulative amount of the rotation angle change after the sequential assembly of segments starting from the starting mileage of the interval for the axis segment where the current shield machine is located. The cumulative rotation angle θ i-1 of the previous ring is β1 - α1 + β2 - α2 + … + β i-1 - α i-1 , and the cumulative rotation angle θ i after the simulated assembly of the segment at the optional point position of the current ring is β1 - α1 + β2 - α2 + … + β i-1 - α i-1 + β i - α i .

[0083] S5 - 1 - 3. Assign a value to the importance coefficient I i-1 of the axis rotation angle according to the magnitude of the cumulative rotation angle θ 1i of the previous ring. When θ i-1 is 0, I 1i is 1, and the value of |θ i-1 | is proportional to the value of I 1i ;

[0084] The larger the value of |θ i-1 |, the larger the value of I 1i . When the designed axis position of the tunnel where the current shield machine is located is a left turn, and the absolute values are the same, when θ[[ID=5..​​​​​​​​​​​​​​​​​​​​​​​​The larger |β - α| is, the smaller a1 is and it gradually approaches 0. When |β - α| is the same, when the current tunnel design axis position of the shield machine is a right turn, the weight when β - α > 0 is equal to a1, and the weight when β - α < 0 is 0.85 × a1; on the contrary, when β - α < 0, the weight is equal to a1, and when β - α > 0, the weight is 0.85 × a1; when the current tunnel design axis position of the shield machine is a straight section, a1 = 0.

[0087] S5-1-5. According to each optional assembly point K of the current ring i The cumulative rotation angle θ calculated after simulated assembly i Assign the effectiveness score a2, where 0 < a2 < 1. When |θ i | = 0, a1 = 1, and the value of |θ i is inversely proportional to the value of a2; the larger |θ i is, the smaller a2 is and it gradually approaches 0.

[0088] S5-1-6. The effectiveness score V of the axis rotation angle after simulated assembly of each optional assembly point Ki 1i = 0.4 × a1 + 0.6 × a2.

[0089] Furthermore, the importance coefficient I i corresponding to the shield tail clearance of each optional assembly point K 2i and the calculation method of the effectiveness score V 2i are as follows:

[0090] S5-2-1. According to the minimum shield tail clearance value δ among the shield tail clearance values at the upper, lower, left, and right positions before the current ring segment is assembled min Assign the importance coefficient I 2i When δ min > δ lim , I 2i = 1; when δ min < δ lim , I 2i > 1, and the value of I 2i is inversely proportional to the value of δ min ; where δ lim is the shield tail clearance critical value;

[0091] I 2i increases as the shield tail clearance value decreases, indicating that the smaller the shield tail clearance before segment assembly, the worse the situation, and the higher the weight of considering the shield tail clearance during segment selection. Take a certain shield tail clearance value as the critical value δ lim When δ min > δ lim , take I 2i = 1, indicating that this shield tail clearance value is relatively safe and there is no need to focus on it when selecting the assembly point of the current ring segment. When δmin <δ lim When, take I 2i >1, and I 2i It increases as the input shield tail gap value decreases.

[0092] S5-2-2, Calculate the selectable assembly points K for the current ring based on the shield tail gap values ​​at the four positions (up, down, left, and right) before the current ring segment assembly. i Simulate the shield tail gap values ​​at the four positions (top, bottom, left, and right) after assembly, and calculate based on the minimum shield tail gap value δ′. min Performance score V 2i Assignment, 0 <V 2i <1, and V 2i Value and δ′ min The value is directly proportional.

[0093] The current ring has several selectable assembly points K. i The specific calculation method for the shield tail gap values ​​at the four positions (top, bottom, left, and right) after simulation assembly is as follows:

[0094] First, calculate the top, bottom, left, and right widths of the standard ring, right-turn ring, and left-turn ring segments respectively;

[0095] The horizontal and vertical projections of the standard ring are both square, therefore the widths of the upper, lower, left, and right segments are all the standard segment width w;

[0096] The widths of the right-turn ring segment at the top, bottom, left, and right sides are calculated as follows:

[0097] S 右·T =wl·Cos[(2m·π / n)+π / 2] / 2 (3)

[0098] S 右·B =wl·Cos[(2m·π / n)-π / 2] / 2 (4)

[0099] S 右·L =w+l·Cos(2m·π / n) / 2 (5)

[0100] S 右·R =wl·Cos(2m·π / n) / 2 (6)

[0101] The widths of the left-turn ring segment at the top, bottom, left, and right sides are calculated as follows:

[0102] S 左·T =wl·Cos[(2m·π / n)-π / 2] / 2 (7)

[0103] S 左·B =wl·Cos[(2m·π / n)+π / 2] / 2 (8)

[0104] S 左·L =wl·Cos(2m·π / n) / 2 (9)

[0105] S 左·R =w+l·Cos(2m·π / n) / 2 (10)

[0106] In the above formula, w is the standard width of the segment, l is the wedge shape of the segment, m is the segment location, and n is the number of segment locations.

[0107] The changes in the gaps between the upper, lower, left, and right shield tails after the simulated assembly of the right-turning ring are calculated as follows:

[0108] Δδ 右·T =S 右·T ·Sin[arctan(S 右·T / 2d-S 右·B / 2d)-arctan(o / qp / q)] (11)

[0109] Δδ 右·B =-S 右·B ·Sin[arctan(S 右·T / 2d-S 右·B / 2d)-arctan(o / qp / q)] (12)

[0110] Δδ 右·L =S 右·L ·Sin[arctan(S 右·L / 2d-S 右·R / 2d)-arctan(o / qp / q)] (13)

[0111] Δδ 右·R =-S 右·R ·Sin[arctan(S 右·L / 2d-S 右·R / 2d)-arctan(o / qp / q)] (14)

[0112] The changes in the gaps between the upper, lower, left, and right shield tails after the simulated assembly of the left-turning ring are calculated as follows:

[0113] Δδ 左·T =S 左·T ·Sin[arctan(S 左·T / 2d-S 左·B / 2d)-arctan(o / qp / q)] (15)

[0114] Δδ 左·B =-S 左·B ·Sin[arctan(S 左·T / 2d-S 左·B / 2d)-arctan(o / qp / q)] (16)

[0115] Δδ 左·L =S 左·L ·Sin[arctan(S 左·L / 2d-S 左·R / 2d)-arctan(o / qp / q)] (17)

[0116] Δδ 左·R =-S 左·R ·Sin[arctan(S 左·L / 2d-S 左·R / 2d)-arctan(o / qp / q)] (18)

[0117] The changes in the gaps between the upper, lower, left, and right shield tails after the standard ring simulation assembly are calculated as follows:

[0118] Δδ 直·T =w·Sin[arctan(p / qo / q)] (19)

[0119] Δδ 直·T =w·Sin[arctan(p / qo / q)] (20)

[0120] Δδ 直·T =w·Sin[arctan(p / qo / q)] (21)

[0121] Δδ 直·T =w·Sin[arctan(p / qo / q)] (22)

[0122] In the above formula, o is the stroke of the upper cylinder, p is the stroke of the lower cylinder, q is the cylinder mounting diameter, and d is the outer diameter of the tube segment.

[0123] The shield tail gap values ​​at the four positions (up, down, left, and right) before the current ring segment assembly are as follows:

[0124] δ i-1·T δ i-1·B δ i-1·L δ i-1·R

[0125] The calculated values ​​for the shield tail gap at the four positions (up, down, left, and right) after the standard ring simulation assembly are as follows:

[0126] δ 直·T =δ i-1·T +Δδ 直·T δ 直·B =δ i-1·B +Δδ 直·B δ直·L =δ i-1·L +Δδ 直·L δ 直·R =δ i-1·R +Δδ 直·R ;

[0127] The calculated values ​​for the shield tail gap at the four positions (up, down, left, and right) after the simulated assembly of the left-turning ring are:

[0128] δ 左·T =δ i-1·T +Δδ 左·T δ 左·B =δ i-1·B +Δδ 左·B δ 左·L =δ i-1·L +Δδ 左·L δ 左·R =δ i-1·R +Δδ 左·R ;

[0129] The calculated values ​​for the shield tail gap at the four positions (up, down, left, and right) after the simulated assembly of the right-turn ring are:

[0130] δ 右·T =δ i-1·T +Δδ 右·T δ 右·B =δ i-1·B +Δδ 右·B δ 右·L =δ i-1·L +Δδ 右·L δ 右·R =δ i-1·R +Δδ 右·R ;

[0131] The score for the shield tail gap effect is assigned a range of 0. <V 2i <1, which increases as the shield tail gap increases, gradually approaching 1, and decreases as the shield tail gap decreases, gradually approaching 0. This indicates that the smaller the shield tail gap after segment assembly, the worse the segment selection effect, and the lower the corresponding effect score; conversely, the larger the gap, the higher the score.

[0132] Furthermore, each optional assembly point K i The importance coefficient I of the corresponding hydraulic cylinder stroke 3i And the effect score V 3i The calculation method is as follows:

[0133] S5-3-1. Based on the cylinder stroke input values ​​at the four positions (up, down, left, and right) before the current ring segment assembly, calculate the stroke difference between the up and down cylinders and the stroke difference between the left and right cylinders, and select the larger cylinder stroke difference value △L. max Importance coefficient I 3i Assignment, when △Lmax >△L lim At that time, I 3i =1; when △L max <△L lim At that time, I 3i >1, and I 3i Value and △L max The value is directly proportional; where ΔL max This is the critical value for the difference in cylinder stroke.

[0134] I 3i The value increases with the increase of the cylinder stroke difference, indicating that the larger the cylinder stroke difference before segment assembly, the worse the situation, and the higher the weight given to the cylinder stroke difference when selecting segments. A certain cylinder stroke difference is taken as the critical value ΔL. lim When △L max >△L lim When, take I 3i =1 indicates that this cylinder stroke difference is beneficial for shield tunneling construction, and does not need to be emphasized when selecting points for segment assembly; when △L max <△L lim When, take I 3i >1, and I 3i It increases as the difference in cylinder stroke increases.

[0135] S5-3-2. Based on the stroke difference between the upper and lower cylinders and the stroke difference between the left and right cylinders before the current ring segment assembly, calculate the possible assembly points K for each ring segment. i Simulated cylinder stroke difference variation after assembly; selectable assembly points K for each front ring. i The calculated value of the cylinder stroke difference after simulated assembly = cylinder stroke difference of the previous ring segment + cylinder stroke difference change value, and based on the maximum cylinder stroke difference value △L′. max Performance score V 3i Assignment, 0 <V 3i <1, and V 2i Value and δ′ min The value is inversely proportional.

[0136] Specifically, given the known stroke difference between the upper and lower, and left and right hydraulic cylinders before the current ring segment assembly, ΔL i-1·H =L 左 -L 右 ΔL i-1·V =L 上 -L 下 Based on formulas (3) to (10), given the upper, lower, left, and right widths of the standard ring, right-turn ring, and left-turn ring segments, calculate the selectable assembly points K for each segment type. i The change in cylinder stroke difference, calculated after simulation assembly = input cylinder stroke difference before assembly + change in cylinder stroke difference:

[0137] The stroke difference between the upper and lower, and left and right cylinders after the simulated assembly of the right-turn ring is:

[0138] L 右·V =ΔL i-1·V -(S 右·T -S 右·B L 右·H =ΔL i-1·H -(S 右·L -S 右·R );

[0139] The stroke differences between the upper and lower, and left and right cylinders after the simulated assembly of the left-turn ring are:

[0140] L 左·V =ΔL i-1·V -(S 左·T -S 左·B L 左·H =ΔL i-1·H -(S 左·L -S 左·R );

[0141] The stroke difference between the upper and lower, left and right cylinders remains unchanged after the standard ring simulation assembly.

[0142] Based on the maximum value of the cylinder stroke difference △L′ max Performance score V 2i The score for the hydraulic cylinder stroke effect is assigned a value within the range of 0. <V 3i <1, which increases as the cylinder stroke difference decreases and gradually approaches 1, and decreases as the cylinder stroke difference increases and gradually approaches 0. This means that the larger the cylinder stroke difference after segment assembly, the worse the segment selection effect and the lower the corresponding effect score. Conversely, the smaller the difference, the higher the score.

[0143] Furthermore, each optional assembly point K i The importance coefficient I of the corresponding tunnel boring machine trend 4i And the effect score V 4i The calculation method is as follows:

[0144] S5-4-1. Calculate the tunnel boring machine trend importance coefficient I based on the current trend value of the tunnel boring machine before the assembly of the ring segments. 4i Assignment, I 4i The value is proportional to the absolute value of the shield machine trend value before the current ring segment assembly. When the current shield machine trend value before the current ring segment assembly is 0, I 4i =0;

[0145] Taking the horizontal direction as an example, a rightward trend of the tunnel boring machine is considered positive; when the tunnel boring machine trend is less than 0, I... 4i >0, tunnel boring trend >0, I 4i <0, and I 4iThe value of I increases as the absolute value of the shield tunneling trend increases, and vice versa. When the shield tunneling trend before segment assembly is 0, I... 4i =0 indicates that the selection of tunnel segments does not need to be considered at this time.

[0146] S5-4-2. The effectiveness score of determining the shield tunneling trend based on the direction of the shield tunneling trend and the direction of the pipe ring capping block is V. 4i .

[0147] The shield tunneling trend effect score V is determined based on the direction of the shield tunneling trend and the direction of the pipe ring capping block. 4i The upper limit, middle value, and lower limit critical values, such as a positive value for the horizontal trend of the shield tunnel, indicate that the shield tunnel is trending towards the right side of the tunnel design axis. Figure 3 As shown in (a), the correction effect of the pipe ring on the left side of the capping block is optimal, with an effect score of the upper limit critical value of 1; Figure 3 As shown in (b), the ring effect of the capping block on the top and bottom sides is 0, and the corresponding effect score is the median value of 0; Figure 3 As shown in (c), the pipe ring with the capping block on the right has the worst effect, with an effect score of -1 (the lower limit critical value); for other pipe rings with the capping block on the left, the effect score is 0. <V 4i <1, the capping block is on the right, effect score -1. <V 4i <0.

[0148] Step S6 specifically includes the following steps:

[0149] S6-1. Based on the selected assembly points and assembly rules of the current ring segment, predict the type of the next ring segment and calculate the possible assembly points K for the next ring segment. i+1 ;

[0150] S6-2. Calculate the possible assembly points K for each segment of the next ring. i+1 Simulated shield tail clearance and hydraulic cylinder stroke difference after assembly;

[0151] S6-3. Based on the selected assembly points of the current ring segment, determine the difference between the shield tail gap and the hydraulic cylinder stroke after the simulated assembly of the selected assembly points of the current ring segment, and calculate the change in shield tail gap and the change in hydraulic cylinder stroke when the next ring shield tunneling is completed.

[0152] S6-4. Calculate the predicted value of the shield tail gap and the predicted value of the hydraulic cylinder stroke difference when the next ring of shield tunneling is completed. Based on step 5, complete the prediction of the assembly point and segment type of the next ring of segments.

[0153] First, based on the selected segment locations and segment assembly rules of the current ring, calculate the possible assembly points K for the next ring segments according to step S4. i+1Secondly, the formulas in steps S5-2-2 and S5-3-2 calculate the selectable assembly points K for each pipe type (standardized, left-turn ring, right-turn ring). i+1 Simulate the difference between the shield tail gap and the cylinder stroke during assembly; determine the difference between the simulated shield tail gap and the cylinder stroke corresponding to the selected segment locations of the current ring; calculate the change in shield tail gap and the change in cylinder stroke when the next ring of shield tunneling is completed according to formulas (23) and (24):

[0154]

[0155] ΔL′=[w·(R ​​+ r / 2) / (R - r / 2) ]- w (24)

[0156] In the above formula, t is the tail shield length, r is the cylinder mounting diameter, w is the tunneling length (segment width), and R is the axis radius.

[0157] The shield tail gap value δ after the current ring segment simulation assembly is known. T δ B δ L δ R The difference between the cylinder stroke and L V L H The predicted tail clearance δ′ when the next ring of shield tunneling is completed is determined based on the direction of the current shield machine axis. T ,δ′ B ,δ′ L ,δ′ R Predicted value L′ of cylinder stroke difference V L′ H :

[0158] If the current axis direction is left turn:

[0159] δ′ T =δ T ,δ′ B =δ B ,δ′ L =δ L -Δδ′、δ′ R =δ R +Δδ′,L′ V =L V L′ H =L H -ΔL′;

[0160] If the current axis direction is right turn:

[0161] δ′ T =δ T ,δ′ B =δ B,δ′ L =δ L +Δδ′、δ′ R =δ R -Δδ′,L′ V =L V L′ H =L H +ΔL′;

[0162] If the current axis direction is a straight segment, the change in shield tail gap and cylinder stroke when the next ring of shield tunneling is completed will be 0. The predicted value is the difference between the shield tail gap and cylinder stroke in the simulated assembly of the current ring segment.

[0163] After determining the tail gap value and cylinder stroke difference prediction value after the next ring of shield tunneling is completed, repeat step 5 to select the segment assembly points and segment types after the next ring of shield tunneling is completed, and complete the prediction of the segment assembly points and types for the next ring.

[0164] Although embodiments of the present invention have been disclosed above, they are not limited to the applications listed in the specification and embodiments. They can be applied to various fields suitable for the present invention. For those skilled in the art, other modifications can be easily made. Therefore, without departing from the general concept defined by the claims and their equivalents, the present invention is not limited to the specific details and illustrations shown and described herein.

Claims

1. A method for shield segment selection and prediction, characterized in that, Includes the following steps: S1. Determine the assembly rules for each ring segment based on the construction site conditions; S2. Divide the tunnel design axis into multiple straight and curved segments, and calculate the alignment parameters of each straight and curved segment in order of mileage. S3. Obtain the alignment parameters corresponding to the current tunnel design axis position of the shield machine based on the data calculated in step S2. S4、According to the last ring pipe splicing point position and the pipe splicing rule, all the selectable splicing point positions K of the current ring pipe are calculated i ; S5, the axis angle of the current ring segment, the shield tail gap, the push cylinder stroke, and the shield trend are weighted and calculated to obtain the weighted score of each selectable assembly point in the current ring segment The assembly point corresponding to the maximum weighted score is selected as the assembly point of the current ring segment ​​​ The weighted score Calculated using the following formula: In the formula , , , The importance coefficients are given to four factors: axis rotation angle, tail shield clearance, propulsion cylinder stroke, and tunnel boring machine trend. , , , For each optional assembly point The simulation scores the effects of four factors: axis rotation angle, tail shield clearance, propulsion cylinder stroke, and tunnel boring machine trend after assembly. Each optional assembly point Importance coefficient of the corresponding axis rotation angle and performance score The calculation method is as follows: S5-1-1 Calculate the available assembly points for each ring segment. Segment rotation ; S5-1-2 Calculate the available assembly points for each ring segment. Cumulative turning angle The curve turning angle α is the turning angle of the current shield tunnel axis segment; S5-1-3, Based on the accumulated rotation angle of the previous ring The importance factor of the size of the axis rotation angle Assignment, When it is 0, It is 1, and Value and The value is directly proportional to the value of ; S5-1-4. Based on the available assembly points of the current ring. The effect score is calculated based on the β-α value after the simulation assembly. Assignment, 0 < <1, when |β-α|=0, =1, |β-α| and The value is inversely proportional to the value of; S5-1-5, Based on the available assembly points of the current ring. Cumulative rotation angle calculated after simulated assembly Values ​​are used to score the effect. Assignment, 0 < <1, when When =0, =1, and The value is inversely proportional to the value of; S5-1-6, Optional assembly points Score for the effect of simulating the rotation angle of the axis after assembly ; S6. Repeat steps S4 and S5 to predict the next ring segment type and segment assembly location.

2. The shield tunnel segment selection and prediction method as described in claim 1, characterized in that, Step S1 specifically includes: S1-1. Choose either staggered or continuous seam assembly; S1-2, Select the assembly points that are prohibited from being selected.

3. The shield tunnel segment selection and prediction method as described in claim 1, characterized in that, In step S2, the linear parameters of the line segment include the length of the line segment. ; The curve segment includes circular curves and transition curves symmetrically positioned on both sides of the circular curve. The linear parameters of the curve segment include the curve length L and the transition curve length. The curve turning angle α and the radius R of the circular curve.

4. The shield tunnel segment selection and prediction method as described in claim 3, characterized in that, Step S4 specifically includes: S4-1. Calculate the preliminary selectable assembly points for the current ring segments; S4-2, eliminate the assembling point position not meeting the current ring pipe sheet assembling rule to obtain all selectable assembling point positions K of the current ring pipe sheet i .

5. The shield tunnel segment selection and prediction method as described in claim 1, characterized in that, Each optional assembly point K i The importance factor of the corresponding shield tail gap and performance score The calculation method is as follows: S5-2-1. Based on the minimum shield tail gap value among the shield tail gap values ​​at the four positions (up, down, left, and right) before the current annular segment assembly. Importance coefficient Assignment, when > hour, =1; when < hour, >1, and Value and The values ​​are inversely proportional; where This is the critical value for the shield tail gap; S5-2-2, Calculate the available assembly points for the current ring based on the shield tail gap values ​​at the four positions (up, down, left, and right) before the current ring segment assembly. Simulate the shield tail gap values ​​at the four positions (top, bottom, left, and right) after assembly, and calculate based on the minimum shield tail gap value. Score the results Assignment, 0 < <1, and Value and The value is directly proportional.

6. The shield tunnel segment selection and prediction method as described in claim 1, characterized in that, Each optional assembly point Importance coefficient of the corresponding hydraulic cylinder stroke and effect score The calculation method is as follows: S5-3-1. Based on the cylinder stroke input values ​​at the four positions (up, down, left, and right) before the current ring segment assembly, calculate the stroke difference between the up and down cylinders and the stroke difference between the left and right cylinders, and select the larger cylinder stroke difference value △L. max Importance coefficient Assignment, when hour, =1; when hour, >1, and Value and The value is directly proportional; among which This is the critical value for the difference in cylinder stroke. S5-3-2. Based on the stroke difference between the upper and lower cylinders and the stroke difference between the left and right cylinders before the current ring segment assembly, calculate the available assembly points for the current ring. Simulated cylinder stroke difference change value after assembly, current ring selectable assembly points. The calculated value of the cylinder stroke difference after simulated assembly = cylinder stroke difference of the previous ring segment + cylinder stroke difference change value, and is based on the maximum cylinder stroke difference. Score the results Assignment, 0 < <1, and Value and The value is inversely proportional.

7. The shield tunnel segment selection and prediction method as described in claim 1, characterized in that, Each optional assembly point K i The importance coefficient of the corresponding tunnel boring machine trend and effect score The calculation method is as follows: S5-4-1. Calculate the importance coefficient of the tunnel boring machine trend based on the current trend value of the tunnel boring machine before the assembly of the ring segments. Assignment, The value is directly proportional to the absolute value of the shield machine's trend value before the current ring segment assembly. When the current shield machine trend value before the current ring segment assembly is 0, =0; S5-4-2. Determine the effectiveness score of shield tunneling trend based on the direction of shield tunneling trend and the direction of the pipe ring capping block. .

8. The shield tunnel segment selection and prediction method as described in claim 1, characterized in that, Step S6 specifically includes the following steps: S6-1. Based on the selected assembly points and assembly rules of the current ring segment, predict the type of the next ring segment and calculate the possible assembly points for the next ring segment. ; S6-2. Calculate the possible assembly points for the next ring segment. Simulated shield tail clearance and hydraulic cylinder stroke difference after assembly; S6-3. Based on the selected assembly points of the current ring segment, determine the difference between the shield tail gap and the hydraulic cylinder stroke after the simulated assembly of the selected assembly points of the current ring segment, and calculate the change in shield tail gap and the change in hydraulic cylinder stroke when the next ring shield tunneling is completed. S6-4. Calculate the predicted value of the shield tail gap and the predicted value of the hydraulic cylinder stroke difference when the next ring of shield tunneling is completed. Based on step 5, complete the prediction of the assembly point and segment type of the next ring of segments.

Images