A tunneling synchronous assembly control method

Abstract

The present invention provides a method for controlling synchronous segment assembly. According to the points and segments to be assembled, the cylinder numbers corresponding to the assembly area, transition area, and synchronization area are calculated, and the cylinders in the corresponding areas are specified to complete the corresponding actions according to the assembly steps. Since each segment point of the shield machine consists of multiple segments, and the number and corresponding propulsion cylinder numbers of each segment are different, when the segment assembler gives the segment point number and segment number, the synchronous segment assembly method of the present invention can quickly and accurately specify the action states for each propulsion cylinder.

Description

Technical Field

[0001] The present invention relates to a shield machine, and particularly to a control method for a shield machine. Background Art

[0002] Simultaneous tunneling and segment erection refers to a construction method in which, during tunnel construction by a shield machine, in order to improve the tunneling speed, segment erection is carried out while the shield machine is tunneling. By tunneling and erecting segments simultaneously, the shield construction efficiency is greatly improved. Especially in engineering projects such as long-distance tunnels and large-section tunnels, the construction period can be significantly shortened, the construction cost can be reduced, and huge economic and social benefits can be generated.

[0003] Since each position of the shield machine corresponds to multiple segments, and the cylinder numbers corresponding to different positions are also different, in order to quickly and accurately give the action state of each cylinder, this has become an urgent problem for us to solve. Summary of the Invention

[0004] Object of the Invention: A simultaneous tunneling and segment erection control method is provided, which judges the state of cylinders in the erection area and gives corresponding implementation steps, and gives the corresponding action states of the erection area, transition area and opposite area.

[0005] Technical Solution: A simultaneous tunneling and segment erection control method calculates the cylinder numbers corresponding to the erection area, transition area and synchronization area according to the position and segment to be erected, and designates the cylinders in the corresponding area to complete the corresponding actions according to the erection steps.

[0006] Further, the sequential number of each ring of segments is 1 to 10, represented by ; it is set that the number of the currently erected segment is Z, where the number of the closing segment is F, the number of the connecting segments is L1 and L2, and the number of the standard segments is B1, B2, B3, B4, B5, B6, B7;

[0007] When Z = F, = 1;

[0008] When Z = L2, = 2;

[0009] When Z = B7, = 3;

[0010] When Z = B6, = 4;

[0011] When Z = B5, = 5; [[ID=5l]]

[0012] When Z = B4, = 6;

[0013] When Z = B3, = 7;

[0014] When Z = B2, = 8;

[0015] When Z = B1, = 9;

[0016] When Z = L1, = 10.

[0017] Furthermore, it is known that the total number of propulsion cylinders is N, the number of propulsion cylinders in the assembly area is n, the number of propulsion cylinders in the transition area is s, and it is assumed that the number of propulsion cylinders in the synchronization area is m;

[0018] The number of propulsion cylinders corresponding to the closure segment is , and there are a total of 9 segments including the connecting segment and the standard segment. The number of propulsion cylinders corresponding to each segment is , and the segment point number Gt ∈ {1, 2,..., }.

[0019] Furthermore, calculate the cylinder numbers q1, q2,..., q j ,..., q n , where j ∈ {2, 3,..., n}:

[0020] When Z = F, , ;

[0021] If ≤ 0, ≤ 0; then = + N, = + N;

[0022] If ﹥ N, ﹥ N; then = - N = - N;

[0023] If N, ; then = , = ;

[0024] When Z = L1, L2, or B1 ~ B7, , ;

[0025] If ≤ 0, ≤ 0; then = +N, = +N;

[0026] If >N, >N; then = -N, = -N;

[0027] If N, ; then = , = .

[0028] In addition, calculate the cylinder numbers W1, W2,..., W k ,...W s , k ∈ {2, 3,..., s};

[0029] The arrangement order of each ring of segments is F-L2-B7-B6-B5-B4-B3-B2-B1-L1. According to this order, judge whether the currently assembled segment arrangement is on the left or right side of the previous segment;

[0030] ① There is no transition area for the assembly of the first segment;

[0031] ② Transition area of the second segment:

[0032] Ⅰ. If the currently assembled segment arrangement is on the left side of the previous segment, then W1 = -s + 1, W k = W1 + k - 1; where is the cylinder number of the nth cylinder in the assembly area of the previous segment, the ;

[0033] If W1 ≤ 0, then W1 = W1 + N;

[0034] If W k >N, then W k = W k -N;

[0035] If 0 < W1 ≤ N, 0 < W k ≤ N, then W1 = W 1， W k = W k ;

[0036] Ⅱ. If the currently assembled segment arrangement is on the right side of the previous segment, then W1 = , W k = W1 + k - 1; It is the number of the first oil cylinder in the segment erection area of the previous segment;

[0037] If W k > N, then W k = W k - N;

[0038] If 0 < W k ≤ N, then W k = W k ;

[0039] ③ For the remaining segment transition areas:

[0040] If the currently erected segment is arranged on the left side of the previous segment, and the previously erected segment is arranged on the right side of the segment before the previous one, then W1 = - s + 1, W k = W1 + k - 1; If the previously erected segment is arranged on the left side of the segment before the previous one, match in reverse order according to the order of the previously erected segments until a segment on the right side of its previous segment or the first segment is matched:

[0041] If W1 ≤ 0, then W1 = W1 + N;

[0042] If W k > N, then W k = W k - N;

[0043] If the currently erected segment is arranged on the right side of the previous segment, and the previous segment is arranged on the left side of the segment before the previous one, then W1 = q1, W k = W1 + k - 1; If the previous segment is arranged on the right side of the segment before the previous one, then until a segment arranged on the left side or the first segment is matched:

[0044] If W k > N, then W k = W k – N;

[0045] If 0 < W k ≤ N, then W k [[ID=6②]]= W k .

[0046] In addition, calculate the number m of the propulsion oil cylinders in the homology area, and take values of the number m of the oil cylinders in the homology area in the set in turn, where Round up and take the even number of 22% of the total number of propulsion cylinders. After each value is taken, calculate the thrust TQi of all propulsion cylinders and take its maximum value maxTQ, that is, the thrust peak value. Calculate a total of M / 2 + 1 times; compare the M / 2 + 1 thrust peak values maxTQ, and take the number of cylinders in the homology region with the smallest thrust peak value as the number of cylinders in the homology region m of the current segment.

[0047] Calculate the numbers P1, P2,…, P of the cylinders in the homology region f ,…P m , f ∈ {2, 3,…, m}:

[0048] ; ;

[0049] If P1 ≤ 0, P f ≤ 0, then P1 = P1 + N, P f = P f + N;

[0050] If P1 > N, P f > N, then P1 = P1 - N, P f ]>= P f - N;

[0051] If 0 < P1 ≤ N, 0 < P f ≤ N; then P1 = P1, P f = P f .

[0052] Furthermore, calculate the total thrust of the shield machine , with the unit of ;

[0053] (1)

[0054] In the formula, is the historical propulsion cylinder pressure, with the unit of MPa; is the ratio of cylinder thrust to pressure;

[0055] In the tunneling and synchronous segment erection mode, the thrust of the propulsion cylinders needs to be distributed according to a gradient to facilitate the control of cylinder pressure. Calculate the distribution value Qi of the force of each propulsion cylinder:

[0056] (2)

[0057] In the formula, represents the angle of the propulsion cylinder in the polar coordinates of the shield body, with the unit of rad, is the angle of the center point of the currently erected segment in the polar coordinates of the shield body, with the unit of rad, is the thrust operation intensity for erecting the current segment;

[0058] Calculate the thrust TQi of each propulsion cylinder:

[0059] (3)

[0060]

[0061] Divide the N propulsion cylinders into two groups, A and B, from the center line of the current segment erection area. Group A, with a deviation coefficient of , Group B, with a deviation coefficient of , there is:

[0062] (4).

[0063] Furthermore, calculate , and :

[0064] (5)

[0065] (6)

[0066] Re - distribute the historical thrust Fi according to formula (2) to obtain the thrust distribution value Pi and obtain the thrust TPi of each cylinder:

[0067] (7)

[0068] (8)

[0069] (9)

[0070] (10)

[0071] In the formula, is the angle at which the historical thrust is located, is the historical thrust operation intensity, is the radius of the circle where the cylinder is located, unit m; is the moment of the total thrust of the historical propulsion cylinders in the X - axis direction; is the moment of the total thrust of the historical propulsion cylinders in the Y - axis direction; is the moment of the total thrust of the propulsion cylinders during the current segment erection in the X - axis direction; is the moment of the total thrust of the propulsion cylinders during the current segment erection in the Y - axis direction;

[0072] In order to achieve thrust stability, that is, the moments in the X - axis and Y - axis directions remain unchanged during segment erection;

[0073] (11)

[0074] (12)

[0075] Solving the simultaneous equations (4), (11) and (12) gives 、 、 :

[0076] Let

[0077] ;

[0078] ;

[0079] ;

[0080] ;

[0081] ;

[0082] ;

[0083] ;

[0084] ;

[0085] ;

[0086] ;

[0087] ;

[0088] (13)

[0089] [[]]END]]is the moment in the X direction of the advancing cylinder in the homology region: [[]]END]]

[0090] [[]]END]] [[]]END]], [[]]END]](14) [[]]END]]

[0091] [[]]END]] [[]]END]]is the moment of the advancing cylinder in the homology region [[]]END]]direction: [[]]END]]

[0092] [[]]END]] [[]]END]](15) [[]]END]]

[0093] [[]]END]] [[]]END]]is the number of cylinders participating in tunneling excluding the cylinders in the homology region, F A [[]]END]]is the total thrust of the homology cylinders; from the above equations, we can get: [[]]END]]

[0094] (16)

[0095] Let again:

[0096] ;

[0097] ;

[0098] ;

[0099] ;

[0100] (17)

[0101] Similarly, let:

[0102] ;

[0103] ;

[0104] ;

[0105] ;

[0106] (18)

[0107] From equation it can be obtained that

[0108]

[0109] .

[0110] Furthermore,

[0111] At the start of assembly:

[0112] The state of the cylinders in the assembly area is initially assembly retraction. After retraction is in place, the state of the cylinders in this area changes to stop;

[0113] The state of the cylinders in the transition area is stop;

[0114] The state of the cylinders in the homology area is follow-up;

[0115] The state of the remaining advancing cylinders is tunneling extension;

[0116] When the segment is placed in place:

[0117] The state of the cylinders in the assembly area is assembly extension;

[0118] The state of the cylinders in the transition area is stop;

[0119] The cylinder state in the homology area is follow-up;

[0120] The states of the remaining propulsion cylinders are extended during tunneling;

[0121] When the assembly is completed:

[0122] The cylinder state in the assembly area is extended during tunneling;

[0123] The cylinder state in the transition area is stopped;

[0124] The cylinder state in the homology area is extended during tunneling;

[0125] The states of the remaining propulsion cylinders are extended during tunneling;

[0126] Finally, the action states of each propulsion cylinder are output.

[0127] Beneficial effects: Since each segment position of the shield machine consists of multiple segments, and the propulsion cylinder numbers and quantities corresponding to each segment are different, the segment position number and segment number are given by the erector. The tunneling synchronous assembly method of the present invention can quickly and accurately specify the action states for each propulsion cylinder. Brief Description of the Drawings

[0128] Figure 1 It is a flow chart of the method of the present invention;

[0129] Figure 2 It is a schematic diagram of cylinders in the segment assembly area and the transition area;

[0130] Figure 3 It is a schematic diagram of cylinders in the homology area;

[0131] Figure 4 It is a flow chart for calculating the peak thrust;

[0132] Figure 5 It is a schematic diagram of the thrust distribution of the tunneling synchronous assembly propulsion system;

[0133] Figure 6 For The calculation flow chart. Detailed Embodiment

[0134] Taking 56 cylinders as an example, the position where the F segment is located is the position number. The assembly sequence of the segments is different from the arrangement sequence of the segments, and this sequence is different due to different position numbers of the segments.

[0135] For the assembly of the B5 segment at the 1st position, the B5 segment is the segment that needs to be assembled at this position currently. The previous segment B4 has been assembled, and B3 is the next segment to be assembled. The cylinders corresponding to each area are artificially specified according to the segments in the assembly area, and the cylinder state is determined according to the steps of segment assembly.

[0136] Point number:

[0137] A tunneling synchronous assembly control method, as Figure 1 shown, the segment number is Z, Z ∈ {F, L1, L2, B1, B2, B3, B4, B5, B6, B7}, where the closing segment number is F, the connecting segment numbers are L1 and L2, and the standard segment numbers are B1, B2, B3, B4, B5, B6, B7. The arrangement order of each ring of segments is the same, and the clockwise sorting is F - L2 - B7 - B6 - B5 - B4 - B3 - B2 - B1 - L1. Let be the sequence number (1 - 10) of a ring of segments:

[0138] When Z = F, = 1;

[0139] When Z = L2, = 2;

[0140] When Z = B7, = 3;

[0141] When Z = B6, = 4;

[0142] When Z = B5, = 5;

[0143] When Z = B4, = 6;

[0144] When Z = B3, = 7;

[0145] When Z = B2, = 8;

[0146] When Z = B1, = 9;

[0147] When Z = L1, = 10.

[0148] It is known that the total number of propulsion cylinders is N, the number of propulsion cylinders in the assembly area is n, the number of propulsion cylinders in the transition area is s, and it is assumed that the number of propulsion cylinders in the synchronization area is m; the number of propulsion cylinders corresponding to the closing segment is , and there are 9 segments in total for the connecting segments and the standard segments, and the number of propulsion cylinders corresponding to each segment is , and the segment point number Gt ∈ {1, 2,..., }. Among them, the number of propulsion cylinders in the area for assembling the closing segment is , and the number of propulsion cylinders in the area of the connecting segments and the standard segments is .

[0149] 1. AsFigure 2 As shown, the segment erection area refers to the area where the propulsion cylinders corresponding to the segment placement area and the cylinders on both sides are located, leaving space for the erector to place the segments; calculate the cylinder numbers q1, q2…, q j ,…q n , j ∈ {2, 3,…, n}:

[0150] When Z = F, , , ;

[0151] If ≤ 0, ≤ 0; then = + N, = + N;

[0152] If ﹥ N, ﹥ N; then = , = - N;

[0153] If N, ; then = , = ;

[0154] When Z = L1, L2, or B1~B7, ,

[0155] , ;

[0156] If ≤ 0, ≤ 0; then = q1 + N, = + N;

[0157] If ﹥ N, ﹥ N; then = q1, = - N.

[0158] If N, ; then = , = .

[0159] 2. Since space for segment assembly needs to be reserved, the same hydraulic cylinders will stop when assembling the next ring of segments. These hydraulic cylinders are called the transition area. The hydraulic cylinders in the transition area are in a stopped state, while the remaining hydraulic cylinders continue to be in the tunneling state.

[0160] Therefore, the cylinder numbers that are the same between the cylinder numbers of the segment assembly area for assembling the previous segment and the cylinder numbers of the segment assembly area for assembling the next segment are the cylinder numbers of the transition area. If there are no identical cylinder numbers in the segment assembly areas of the previous segment and the next segment, then compare them with the cylinder numbers of the segment two segments before or two segments after in sequence until the same propulsion cylinder numbers are found.

[0161] The arrangement order of each ring of segments is F - L2 - B7 - B6 - B5 - B4 - B3 - B2 - B1 - L1. Judge whether the currently assembled segment is on the left or right side of the previous segment according to this order; if the currently assembled segment is the ηth segment to be assembled, where η ∈ {2, 10}, then the previous segment refers to the (η - 1)th assembled segment, the segment two segments before refers to the (η - 2)th assembled segment, and the first assembled segment refers to the 1st assembled segment.

[0162] Calculate the cylinder numbers W1, W2…, W k ,…W s , k ∈ {2, 3,…, s};

[0163] ① There is no transition area when assembling the first segment;

[0164] ② Transition area of the second segment:

[0165] Ⅰ. If the currently assembled segment is on the left side of the previous segment, then W1 = -s + 1, W k = W1 + k - 1; where is the nth cylinder number in the segment assembly area of the previous segment, ;

[0166] If W1 ≤ 0, then W1 = W1 + N;

[0167] If W k ﹥ N, then W k [ = W k - N;

[0168] If 0﹤W1 ≤ N, 0﹤W k ≤ N, then W1 = W 1， W k = W k ;

[0169] Ⅱ. If the currently assembled segment is on the right side of the previous segment, then W1 = ,Wk = W1 + k - 1; is the number of the first cylinder in the assembly area of the previous segment;

[0170] If W k > N, then W k = W k - N;

[0171] If 0 < W k ≤ N, then W k = W k ;

[0172] ③ For the remaining segment transition areas:

[0173] If the currently assembled segment is arranged on the left side of the previous segment, and the previously assembled segment is arranged on the right side of the segment before the previous one, then W1 = - s + 1, W k = W1 + k - 1; When the previously assembled segment is arranged on the left side of the segment before the previous one, match in reverse order according to the order of the previously assembled segments until a segment that is on the right side of its previous segment or the first segment is matched. At this time is the number of the nth cylinder in the assembly area of the segment that is matched to be on the right side or the first segment:

[0174] If W1 ≤ 0, then W1 = W1 + N;

[0175] If W k > N, then W k = W k - N;

[0176] If the currently assembled segment is arranged on the right side of the previous segment, and the previous segment is arranged on the left side of the segment before the previous one, then W1 = q1, W k = W1 + k - 1, is the number of the first cylinder in the assembly area of the previous segment; When the previous segment is arranged on the right side of the segment before the previous one, then match until a segment that is arranged on the left side or the first segment is matched. At this time is the number of the first cylinder in the assembly area of the segment that is matched to be on the left side or the first segment:

[0177] If W k > N, then W k = W k – N;

[0178] If 0 < W k ≤ N, then W k = W k .

[0179] 3. The cylinders in the homology region will have a pressure of 0 during the synchronous assembly in tunneling. To maintain the stability of the force application point during the propulsion process, the cylinders in the diagonal side area of the segment need to be adjusted to homology, that is, the pressure of the cylinders in this area is reduced to 3 Mpa to maintain the stability of the force application point, and its total thrust is F A .

[0180] As Figure 3 shown, during assembly, in order to maintain the force balance, the cylinders in the homology region should be on the diagonal side of the assembly region. For the stability of the system, the number m of cylinders in the homology region should be less than the maximum allowable number M of missing propulsion cylinders. M is the integer and even number obtained by taking 22% of the total number of propulsion cylinders. The diagonal side cylinders are as shown in the figure

[0181] As Figure 4 shown, calculate the number m of propulsion cylinders in the homology region, and take values of the number m of cylinders in the homology region in the set in turn, where is the integer and even number obtained by taking 22% of the total number of propulsion cylinders. After each value is taken, calculate the thrust TQi of all propulsion cylinders and take its maximum value maxTQ, that is, the thrust peak value, and calculate a total of M / 2 + 1 times; compare the M / 2 + 1 thrust peak values maxTQ, and take the number of cylinders in the homology region with the smallest thrust peak value as the number m of cylinders in the homology region of the current segment

[0182] Calculate the cylinder numbers P1, P2…, P f ,…P m , f ∈ {2, 3,…, m}: [[]]

[0183] [[]] [[]]; [[]]; [[]]

[0184] [[]]If P1 ≤ 0, P f ≤ 0, then P1 = P1 + N, P f = P f + N; [[]]

[0185] [[]]If P1 ﹥ N, P f ﹥ N, then P1 = P1 - N, P f = P f - N; [[]]

[0186] [[]]If 0 ﹤ P1 ≤ N, 0 ﹤ P f ≤ N; then P1 = P1, P f = P f . [[]]

[0187] [[]]As Figure 5 shown, first, calculate the total thrust of the shield machine [[]], in the unit of [[]]; [[]]

[0188] (1)

[0189] Wherein, is the pressure of the historical propulsion cylinder, with the unit of MPa; is the ratio of the cylinder thrust to the pressure;

[0190] In the tunneling synchronous assembly mode, the thrust of the propulsion cylinder needs to be distributed according to a certain gradient to facilitate the control of the cylinder pressure. Calculate the distribution value Qi of the force of each propulsion cylinder:

[0191] (2)

[0192] Wherein, represents the angle of the propulsion cylinder in the polar coordinates of the shield body, with the unit of , is the angle of the center point of the currently assembled segment in the polar coordinates of the shield body, with the unit of , is the thrust operation intensity for assembling the currently assembled segment;

[0193] Calculate the thrust TQi of each propulsion cylinder:

[0194] (3)

[0195] is , ;

[0196] As Figure 6 shown, divide the N propulsion cylinders into two groups A and B from the center line of the currently assembled segment area, is group A, and its deviation coefficient is , is group B, and its deviation coefficient is , there is:

[0197] (4)

[0198] Calculate

[0199] (5)

[0200] (6)

[0201] According to formula (2), redistribute the historical thrust Fi to obtain the thrust distribution value Pi and obtain the thrust TPi of each cylinder:

[0202] (7)

[0203] (8)

[0204] (9)

[0205] (10)

[0206] Wherein, is the angle where the historical thrust is located, is the historical thrust operation intensity, is the radius of the circle where the oil cylinder is located, unit m; is the torque of the total thrust of the historical propulsion oil cylinder in the X-axis direction; is the torque of the total thrust of the historical propulsion oil cylinder in the Y-axis direction; is the torque of the total thrust of the propulsion oil cylinder in the X-axis direction during the current assembly; is the torque of the total thrust of the propulsion oil cylinder in the Y-axis direction during the current assembly;

[0207] In order to achieve thrust stability, that is, the torques in the X-axis direction and Y-axis direction remain unchanged during assembly.

[0208] (11)

[0209] (12)

[0210] By solving the simultaneous equations (4), (11), and (12), we get :

[0211] Let:

[0212] ;

[0213] ;

[0214] ;

[0215] ;

[0216] ;

[0217] ;

[0218] ;

[0219] ;

[0220] ;

[0221] ;

[0222] ;

[0223] (13)

[0224] is the co - homology region propulsion oil cylinder Direction torque:

[0225] , (14)

[0226] MQaY is the co - homology region propulsion oil cylinder Direction torque:

[0227] (15)

[0228] is the number of oil cylinders participating in tunneling except for the co - homology region oil cylinders, F A is the total thrust of the co - homology oil cylinders; From the above equations, we can get:

[0229] (16)

[0230] Let again:

[0231] ;

[0232] ;

[0233] ;

[0234] ;

[0235] (17)

[0236] Similarly, let:

[0237] ;

[0238] ;

[0239] ;

[0240] ;

[0241] (18)

[0242] From equation we can get

[0243]

[0244] 。

[0245] 4. Output the status value Yi of the propulsion cylinder, where i = 1 to N.

[0246] The status of the propulsion cylinder is divided into the forward cylinder stopped (with the value a), the propulsion cylinder tunneling extended (with the value b), the propulsion cylinder following (with the value c), the propulsion cylinder assembling extended (with the value d), and the propulsion cylinder assembling retracted (with the value e). The action status values of the same propulsion cylinder in different assembling steps are different.

[0247] When starting to assemble:

[0248] The status of the cylinders in the assembling area is first assembling retracted. After retracting in place, the status of the cylinders in this area changes to stopped.

[0249] The status of the cylinders in the transition area is stopped.

[0250] The status of the cylinders in the homology area is following.

[0251] The status of the remaining propulsion cylinders is tunneling extended.

[0252] When the segment is placed in place:

[0253] The status of the cylinders in the assembling area is assembling extended.

[0254] The status of the cylinders in the transition area is stopped.

[0255] The status of the cylinders in the homology area is following.

[0256] The status of the remaining propulsion cylinders is tunneling extended.

[0257] When the assembly is completed:

[0258] The status of the cylinders in the assembling area is tunneling extended.

[0259] The status of the cylinders in the transition area is stopped.

[0260] The status of the cylinders in the homology area is tunneling extended.

[0261] The status of the remaining propulsion cylinders is tunneling extended.

[0262] Finally, output the action status Yi = (a, b, c, d, e) of each propulsion cylinder.

Claims

1. A method for controlling synchronous assembly of tunneling, characterized by: According to the points and segments to be assembled, the cylinder numbers corresponding to the assembly area, transition area, and synchronization area are calculated, and the cylinders in the corresponding areas are assigned to complete the corresponding actions according to the assembly steps; The sequence of each ring segment is numbered from 1 to 10. Indicates: Set the number of the currently assembled segment to Z, the capping block to F, the connecting blocks to L1 and L2, and the standard blocks to B1, B2, B3, B4, B5, B6, and B7; When Z=F, =1; When Z=L2, =2; When Z=B7, =3; When Z=B6, =4; When Z=B5, =5; When Z=B4, =6; When Z=B3, =7; When Z=B2, =8; When Z=B1, =9; When Z=L1, =10; Given that the total number of propulsion cylinders is N, the number of propulsion cylinders in the assembly area is n, the number of propulsion cylinders in the transition area is s, and the number of propulsion cylinders in the coherence area is m; The number of thrust cylinders corresponding to the capping block segment is , there are 9 segments in total, including connecting segments and standard segments. The number of propulsion cylinders corresponding to each segment is , segment point number Gt∈{1,2,…, }; Calculate the cylinder numbers q1, q2…, q in the assembly area j ,…q n , j∈{2, 3, …, n}: When Z=F, , ; like ≤0, ≤0; then = +N, = +N; like ﹥N, >N; then = -N = -N; like N, ;but = , = ; When Z=L1, L2, or B1~B7, , ; like ≤0, ≤0; then = +N, = +N; like ﹥N, >N; then = -N, = -N; like N, ;but = , = .

2. The tunneling synchronous assembly control method according to claim 1, characterized in that: Calculate the transition area cylinder numbers W1, W2…, W k ,…W s , k∈{2, 3, …, s}; The arrangement order of each ring of segments is F-L2-B7-B6-B5-B4-B3-B2-B1-L1. According to this order, the current segment is arranged on the left or right side of the previous segment. ① There is no transition area when assembling the first segment; ② Transition area of the second segment: Ⅰ. If the current segment is arranged on the left side of the previous segment, then W1= -s+1,W k =W1+k-1; where Number the nth cylinder in the assembly area of the previous segment. ; If W1≤0, then W1=W1+N; If W k ﹥N, then W k =W k -N; If 0﹤W1≤N,0﹤W k ≤N, then W1=W 1， W k =W k ; Ⅱ. If the current segment is arranged on the right side of the previous segment, then W1= , W k =W1+k-1; Number the first oil cylinder in the assembly area of the previous segment; If W k ﹥N, then W k =W k -N; If 0﹤W k ≤N, then W k =W k ; ③ Transition area of other segments: If the currently assembled segment is arranged on the left side of the previous segment, and the previous segment is arranged on the right side of the segment before that, then W1= -s+1,W k =W1+k-1; the previously assembled segment is arranged on the left side of the segment above it, and the segments are matched in reverse order according to the sequence of the assembled segments until the segment to the right of the segment above it or the first segment is matched: If W1≤0, then W1=W1+N; If W k ﹥N, then W k = W k -N; If the current segment is arranged on the right side of the previous segment, and the previous segment is arranged on the left side of the segment before that, then W1=q1, W k =W1+k-1; if the previous segment is arranged on the right side of the previous segment, then until it matches the segment arranged on the left side or the first segment: If W k ﹥N, then W k =W k –N; If 0﹤W k ≤N, then W k =W k .

3. The tunneling synchronous assembly control method according to claim 1, characterized in that: Calculate the number of propulsion cylinders m in the coherent region and put the number of cylinders m in the coherent region into the set The values are taken in sequence, where Round up to 22% of the total number of thrust cylinders. After each value is taken, calculate the thrust TQi of all thrust cylinders and take out its maximum value maxTQ, i.e., the thrust peak value. Calculate this value M / 2+1 times in total. Compare the M / 2+1 thrust peak values maxTQ. The number of cylinders in the coherent region with the smallest thrust peak value is taken as the number m of cylinders in the coherent region of the current segment. Calculate the cylinder numbers P1, P2…, P in the synchronization area f ,…P m , f∈{2, 3, …, m}: ； ； If P1≤0, P f ≤0, then P1=P1+N, P f =P f +N; If P1>N, P f ﹥N, then P1=P1-N, P f =P f -N; If 0﹤P1≤N, 0﹤P f ≤N; then P1=P1, P f =P f .

4. The tunneling synchronous assembly control method according to claim 3, characterized in that: Calculate the total thrust of the shield machine , the unit is ; （1） Where, is the historical propulsion cylinder pressure, unit: MPa; is the ratio of cylinder thrust to pressure; The synchronous excavation assembly mode requires that the thrust of the propulsion cylinder be distributed according to a gradient to facilitate cylinder pressure control and calculate the distribution value Qi of each propulsion cylinder force: （2） Where, Indicates the angle of the propulsion cylinder in the shield polar coordinates, in rad. The angle of the center point of the current assembled segment in the shield polar coordinates, in rad. The thrust operating strength for assembling the currently assembled segments; Calculate the thrust TQi of each propulsion cylinder: （3） for or ; Divide the N thrust cylinders into two groups, A and B, based on the centerline of the current segment assembly area. is group A, and its coefficient of variation is αA. is group B, and its coefficient of variation is ,have: （4）。 5. The tunneling synchronous assembly control method according to claim 4, characterized in that: calculate 、 and : （5） （6） According to formula (2), the historical thrust Fi is redistributed to obtain the thrust distribution value Pi And the thrust TPi of each cylinder is obtained: （7） （8） （9） （10） Where, From the perspective of historical thrust, is the historical thrust operating strength, is the radius of the circle where the oil cylinder is located, in m; is the total thrust torque of the historical propulsion cylinder in the X-axis direction; is the total thrust moment of the historical propulsion cylinder in the Y-axis direction; It is the total thrust torque of the propulsion cylinder in the X-axis direction during the current assembly; It is the total thrust moment of the propulsion cylinder in the Y-axis direction during the current assembly; In order to achieve thrust stability, that is, the torque in the X-axis and Y-axis directions remains unchanged during assembly; （11） （12） The simultaneous equations (4), (11), and (12) yield : make: ; ； ； ； ； ； ； ； ； ； ； ； （13） The X-direction torque of the thrust cylinder in the coherent area is: ,make (14) The Y-direction torque of the thrust cylinder in the coherent area is: （15） is the number of cylinders involved in tunneling except the cylinders in the coherent area, F A is the total thrust of the synchronous cylinder; from the above equation, we can get: （16） Then order: ； ； ； ； （17） Similarly, let: ； ； ； ； （18） From formula (17) (18) we can get 。 6. The tunneling synchronous assembly control method according to claim 1, characterized in that: When you start assembling: The state of the oil cylinder in the assembly area is first assembly retraction, and after retraction into place, the state of the oil cylinder in this area turns to stop; The state of the cylinder in the transition area is stopped; The state of the cylinder in the synchronization area is follow-up; The status of the remaining propulsion cylinders is excavation extension; When the segments are in place: The status of the oil cylinder in the assembly area is assembly extension; The state of the cylinder in the transition area is stopped; The state of the cylinder in the synchronization area is follow-up; The status of the other thrust cylinders is excavation extension; When assembly is complete: The status of the oil cylinder in the assembly area is excavation extension; The state of the cylinder in the transition area is stopped; The status of the cylinder in the synchronization area is advancing and extending; The status of the remaining propulsion cylinders is excavation extension; Finally, the action status of each propulsion cylinder is output.

Images