CN113404499B - Real-time control method for attitude of shield type heading machine in complex stratum - Google Patents

Real-time control method for attitude of shield type heading machine in complex stratum Download PDF

Info

Publication number
CN113404499B
CN113404499B CN202110785850.7A CN202110785850A CN113404499B CN 113404499 B CN113404499 B CN 113404499B CN 202110785850 A CN202110785850 A CN 202110785850A CN 113404499 B CN113404499 B CN 113404499B
Authority
CN
China
Prior art keywords
shield
tunneling machine
type tunneling
deviation
vertical
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110785850.7A
Other languages
Chinese (zh)
Other versions
CN113404499A (en
Inventor
唐亮
孔祥勋
凌贤长
唐文冲
张艺帆
丛晟亦
胡心雨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Chongqing Research Institute of Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Chongqing Research Institute of Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology, Chongqing Research Institute of Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202110785850.7A priority Critical patent/CN113404499B/en
Publication of CN113404499A publication Critical patent/CN113404499A/en
Application granted granted Critical
Publication of CN113404499B publication Critical patent/CN113404499B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21DSHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D9/00Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
    • E21D9/003Arrangement of measuring or indicating devices for use during driving of tunnels, e.g. for guiding machines
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21DSHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D9/00Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
    • E21D9/06Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
    • E21D9/093Control of the driving shield, e.g. of the hydraulic advancing cylinders

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Environmental & Geological Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

A method for controlling the posture of a shield-type tunneling machine in a complex stratum in real time belongs to the technical field of tunnel shield-type tunneling machine construction, and the specific scheme comprises the following steps: measuring the posture and the position of a shield type tunneling machine, and judging whether the shaft line offset s of the shield type tunneling machine reaches an allowable offset [ s ], namely judging whether posture correction is needed; secondly, when deviation correction is needed, planning a deviation correction path according to the offset s and the posture of the shaft line of the shield type tunneling machine at the current position, determining single-ring deviation correction amount and single-ring deviation correction angle, and calculating needed deviation correction torque; and step three, correspondingly adjusting the hydraulic propulsion system, and performing deviation correction work on the shield type tunneling machine. The invention adopts the shield-type heading machine posture real-time control method with theoretical basis to correct the deviation of the shield-type heading machine, can quantitatively determine the deviation correcting path and the deviation correcting moment, and effectively avoids the problems of under-correction and over-correction, thereby reducing the problems of segment damage, path deviation and potential operation.

Description

Real-time control method for attitude of shield type heading machine in complex stratum
Technical Field
The invention belongs to the technical field of tunnel shield type tunneling machine construction, and particularly relates to a method for controlling the posture of a shield type tunneling machine in a complex stratum in real time.
Background
In the tunneling process of the shield type tunneling machine, the moment of the shield type tunneling machine is easy to be unbalanced due to reasons such as stratum change and the like, so that the axis of the shield type tunneling machine deviates. In the prior art, according to attitude information of a shield-type tunneling machine obtained by an automatic measuring device, a driver of the shield-type tunneling machine corrects the deviation of the shield-type tunneling machine by virtue of personal operation experience, and the deviation correcting effect depends on the experience and decision level of an operator. Different technicians have different knowledge and understanding of the deviation correcting target value and the deviation correcting load of the shield type tunneling machine. An operator qualitatively predicts a deviation correcting target corner based on personal experience and determines the problems of under-correction and over-correction possibly caused by deviation correcting load: deviation quality problems can lead to segment damage, route drift, and potential operational problems. Therefore, it is necessary to provide a method for real-time quantitative control of the attitude of the shield-type heading machine based on the mechanical balance theory.
Disclosure of Invention
The invention aims to solve the problem that the existing shield-type tunneling machine posture adjustment is carried out by depending on manual experience, and provides a real-time control method for the posture of a shield-type tunneling machine in a complex stratum.
In order to achieve the purpose, the technical scheme adopted by the invention is as follows:
a method for controlling the posture of a shield-type heading machine in a complex stratum in real time comprises the following steps:
measuring the posture and the position of a shield type tunneling machine, and judging whether the shaft line offset s of the shield type tunneling machine reaches an allowable offset [ s ], namely judging whether posture correction is needed;
secondly, when deviation correction is needed, planning a deviation correction path according to the offset s and the posture of the shaft line of the shield type tunneling machine at the current position, determining single-ring deviation correction amount and single-ring deviation correction angle, and calculating needed deviation correction torque;
and step three, correspondingly adjusting the hydraulic propulsion system, and performing deviation correction work on the shield type tunneling machine.
Further, in the first step, the distance s of the central point of the shield type tunneling machine deviating from the design axis on the vertical plane and the horizontal plane of the current position is obtained by utilizing the automatic attitude measuring device of the shield type tunneling machine i And an included angle eta between the axis of the shield type tunneling machine and the designed axis, wherein the vertical plane i is equal to v, and the horizontal plane i is equal to h.
Further, in the second step, a double-circular-arc interpolation method is adopted to determine the deviation rectifying path, and the deviation rectifying path is determined according to the formula (1)Fixed deviation-correcting radius R rec (ii) a Calculating the deviation-correcting mileage L and the central angle theta of the double arcs according to the equation set (2) 1 And theta 2 (ii) a Determining the number k of segment rings required by rectification according to a formula (3); determining a curve equation of the double-arc deviation rectifying path according to a formula (4), so as to obtain the coordinates of the shield type tunneling machine at each ring position on the deviation rectifying path, wherein the difference of longitudinal coordinates at two adjacent ring positions is the deviation of each ring, and the included angle between the connecting line at the two adjacent ring positions and the design axis is the deviation rectifying angle of each ring;
Figure BDA0003159193790000021
Figure BDA0003159193790000022
k=[L/g]+1 (3)
Figure BDA0003159193790000023
in the formula, l is the length of the shield type tunneling machine, and d is the wedge-shaped amount of the duct piece; eta is the included angle between the axis of the shield type tunneling machine and the design axis, and g is the width of the single-ring duct piece; s i For the shaft line offset of the shield type tunneling machine, a vertical plane i is equal to v, and a horizontal plane i is equal to h; (x) 1 ,z 1 )
And (y) 1 ,z 1 ) Is the center O of a first section of circular arc on the horizontal and vertical surfaces 1 And determining a deviation rectifying path consisting of two sections of circular arcs by a double circular arc interpolation method, wherein the deviation rectifying path is respectively represented by AQ and QB.
Further, in the second step, the calculation method of the correction torque is as follows:
1) calculating the section displacement of the shield-type tunneling machine and the load around the shield shell according to the single-ring deviation correction angle beta of the shield-type tunneling machine in the vertical plane, the single-ring deviation correction angle alpha in the horizontal plane and the initial single-ring deviation correction quantity delta, substituting the formula (6), and judging whether the force balance condition is met;
2) if so, determining the correction torque according to the formula (7);
3) if not, after correcting the deviation correction angle of the shield type tunneling machine, repeating the steps 1) and 2) until the deviation correction torque is calculated
Figure BDA0003159193790000024
Figure BDA0003159193790000025
In the formula (I), the compound is shown in the specification,
Figure BDA0003159193790000026
a vertical adjusting moment for controlling the vertical pitch angle generated by the jack thrust difference,
Figure BDA0003159193790000027
a horizontal adjusting moment for controlling a horizontal swinging angle generated by the difference of jack thrust; f 1v Is the gravity of a shield type tunneling machine; f 5v Vertical load resultant force around the shield shell is caused by attitude change; f 5h The resultant force of horizontal load around the shield shell is caused by the change of the posture; m 1v The moment is caused by the eccentric distance between the gravity center and the centroid of the shield type tunneling machine; m 4v The moment is a vertical moment caused by the front load of a cutter head of the shield type tunneling machine; m 5v Vertical resultant moment generated by load around the shield shell due to attitude change; m 5h The horizontal resultant moment generated by the load around the shield shell is caused by the change of the posture.
Further, solve for F 5v 、F 5h 、M 1v 、M 4v 、M 5h And M 5v The calculation steps are as follows:
shield shell section centroid displacement delta s caused by attitude change v And Δ s h Calculated according to equation (8):
Figure BDA0003159193790000031
in the formula,. DELTA.s v Is the vertical displacement of the centroid of the section of the shield shell, and delta s is the vertical displacement of the centroid caused by the gravity action of the shield-type tunneling machine v Adding vertical displacement to the centroid caused by the action of attitude adjusting moment, wherein z is a coordinate along the length direction of the shield type tunneling machine, L is the length of the shield type tunneling machine, and delta s h Is the horizontal displacement of the centroid of the shield shell section, delta h Adding horizontal displacement to a centroid caused by the action of the attitude adjusting moment, wherein beta is a single-ring deviation rectifying angle of the shield type tunneling machine in a vertical plane, and alpha is a single-ring deviation rectifying angle of the shield type tunneling machine in a horizontal plane;
the curve equation of any point on the section of the shield shell is shown as the formula (9):
Figure BDA0003159193790000032
converting equation (9) to polar equation (10):
Figure BDA0003159193790000033
in the formula, R c The mean radius of the shield shell is shown, and theta is an included angle between a connecting line of a point on the shield shell and the circle center and a polar axis;
displacement of any point on the shield section
Figure BDA0003159193790000034
Can be determined by equations (11) and (12):
Figure BDA0003159193790000035
Figure BDA0003159193790000036
wherein r is the coordinate of the polar axis;
determining the full-line stratum distribution according to the geological survey report, and obtainingDetermining the parameters a and K in the foundation reaction force curve (5) of each layer of soil according to the engineering mechanical characteristics of each geotechnical layer 0 ,K min And K max
Figure BDA0003159193790000037
Wherein K is the soil pressure coefficient, i is v or h, which respectively represent the vertical direction and the horizontal direction, K i0 Is the coefficient of static soil pressure, K imin Is the coefficient of active earth pressure, a i Is a foundation reaction coefficient, U i For displacement of the soil mass in the direction i, K imax Is the passive soil pressure coefficient;
and (3) loading around the shield shell: dispersing the shield shell, wherein the shield shell respectively comprises m units and n units along the length direction and the circumferential direction,
the earth pressure at any point on the shield section can be determined by equation (13):
Figure BDA0003159193790000041
wherein p represents the p-th cell in the longitudinal direction, q represents the q-th cell in the circumferential direction, σ v,pq To node vertical pressure, K v,pq Is the vertical soil pressure coefficient, sigma v0,pq To initial vertical pressure, K h,pq Is the horizontal soil pressure coefficient, sigma h,pq To node horizontal pressure, σ h0,pq Is an initial horizontal pressure;
solving the attitude adjusting moment: the gravity of the shield type tunneling machine is G, and the eccentricity is l s And then:
M 1v =G·l s (14)
the formation reaction force and the moment thereof are as follows:
Figure BDA0003159193790000042
Figure BDA0003159193790000043
the cutter head load moment is as follows:
M 4v =FM 1 +FM 2 (17)
in the formula FM 1 Moment caused by the penetration resistance of the cutter head; FM 2 Moment caused by stratum lateral soil pressure on a cutter head panel; d is the average diameter of the shield shell.
Compared with the prior art, the invention has the beneficial effects that: the method can effectively avoid the problems of under-correction and over-correction compared with the traditional correction strategy of the shield-type tunneling machine based on manual experience, thereby reducing the damage of the duct piece, the deviation of the route and the potential operation problem.
Drawings
FIG. 1 is a flow chart of the operation of a method for controlling the attitude of a shield-type heading machine in real time;
FIG. 2 is a schematic diagram of a target deviation rectifying path;
FIG. 3 is a schematic diagram of a reaction force curve of a stratum foundation;
FIG. 4 is a schematic diagram of the displacement of the shield shell section caused by the attitude change;
FIG. 5 is a flow chart of the attitude deviation rectifying moment solution.
Detailed Description
The present invention will be described in detail with reference to the accompanying drawings.
The first embodiment is as follows:
a real-time control method for the attitude of a shield-type heading machine in a complex stratum is shown in the general route of figure 1. And obtaining the attitude and the position of the shield-type tunneling machine according to the automatic attitude measuring device of the shield-type tunneling machine, and judging whether the shaft line offset s of the shield-type tunneling machine reaches the allowable offset [ s ], namely judging whether attitude correction is needed. When the correction is needed, planning a correction path according to the offset s and the posture of the shaft line of the shield type tunneling machine at the current position, determining the single-ring correction amount and the single-ring correction angle, and calculating the needed correction torque; and the hydraulic propulsion system makes corresponding adjustment to perform deviation correction work of the shield type tunneling machine.
The method comprises the following specific steps:
(1) automatic measurement device of the posture of the shield-type heading machine: the existing automatic total station and gyroscope combined guide system or other more advanced and accurate shield type heading machine attitude information measuring devices can be adopted. The automatic attitude measuring device of the shield type tunneling machine can acquire the distance s of the center point of the shield type tunneling machine deviating from the design axis on the vertical plane and the horizontal plane of the current position i And an included angle eta between the axis of the shield type heading machine and the designed axis, wherein the vertical plane i is equal to v, and the horizontal plane i is equal to h.
(2) Deviation rectifying path planning system: determining the deviation rectifying path by adopting a double-arc interpolation method, as shown in figure 2, determining the deviation rectifying radius R according to the formula (1) rec (ii) a Calculating the deviation-correcting mileage L and the central angle theta of the double arcs according to the equation set (2) 1 And theta 2 (ii) a Determining according to a formula (3) to obtain the number k of the segment rings required by rectification; and (4) determining a curve equation of the double arcs according to the formula (4) so as to obtain the coordinates of the shield type tunneling machine at each ring position on the deviation rectifying path, wherein the difference of longitudinal coordinates at the positions of two adjacent rings is the deviation of each ring, and the included angle between the connecting line at the positions of two adjacent rings and the design axis is the deviation rectifying angle of each ring.
Figure BDA0003159193790000051
Figure BDA0003159193790000052
k=[L/g]+1 (3)
Figure BDA0003159193790000061
In the formula, l is the length of the shield type tunneling machine, and d is the wedge-shaped amount of the duct piece; g is the width of the single-ring pipe piece; s i For the shaft line offset of the shield type tunneling machine, a vertical plane i is equal to v, and a horizontal plane i is equal to h; (x) 1 ,z 1 ) And (y) 1 ,z 1 ) And determining a deviation rectifying path consisting of two sections of circular arcs by adopting a double circular arc interpolation method for the circle center coordinates of the first section of circular arc on the horizontal and vertical surfaces, wherein the deviation rectifying path is respectively represented by AQ and QB. (3) A stratum model: determining the all-line stratum distribution according to the geological survey report, acquiring the engineering mechanical characteristics of each geotechnical layer, and determining
Parameters a, K in the foundation reaction force curve (5) of the layers of soil in FIG. 3 0 ,K min And K max
Figure BDA0003159193790000062
Wherein K is the soil pressure coefficient; i is v or h, which respectively represent the vertical direction and the horizontal direction; k i0 Is the coefficient of static soil pressure; k imin Is the active soil pressure coefficient; a is i Is the foundation reaction coefficient; u shape i The displacement of the soil body in the direction i; k imax Is the passive soil pressure coefficient.
(4) Deviation rectifying moment model
The shield type heading machine satisfies the force balance condition and the moment balance condition shown in the formulas (6) and (7):
Figure BDA0003159193790000063
Figure BDA0003159193790000064
in the formula F 1v Is the gravity of a shield type tunneling machine; f 5v Vertical load resultant force around the shield shell is caused by attitude change; f 5h The resultant force of horizontal load around the shield shell is caused by the change of the posture; m 1v The moment is caused by the eccentric distance between the gravity center and the centroid of the shield type tunneling machine; m 4v The moment is a vertical moment caused by the front load of a cutter head of the shield type tunneling machine; m 5v Vertical resultant moment generated by load around the shield shell due to attitude change;
Figure BDA0003159193790000065
a vertical adjusting moment for controlling a vertical pitch angle generated by the difference of jack thrust; m 5h The horizontal resultant moment generated by the load around the shield shell caused by the attitude change;
Figure BDA0003159193790000066
the horizontal adjusting moment for controlling the horizontal swing angle generated by the jack thrust difference.
Shield shell section displacement caused by attitude change
As shown in fig. 4, the centroid displacement of the shield shell section caused by the attitude change is calculated according to equation (8):
Figure BDA0003159193790000071
in the formula,. DELTA.s v Vertical displacement of the section centroid of the shield type tunneling machine; delta s is the vertical displacement of the centroid caused by the gravity action of the shield type tunneling machine; delta v Adding vertical displacement to the centroid caused by the action of the attitude adjusting moment; z is a coordinate along the length direction of the shield type tunneling machine; l is the length of the shield type tunneling machine; Δ s h The horizontal displacement of the section centroid of the shield type tunneling machine is adopted; delta h Adding horizontal displacement to the centroid caused by the action of the attitude adjusting moment; beta is a single-ring deviation rectifying angle of the shield type tunneling machine in a vertical plane; alpha is a single-ring deviation rectifying angle of the shield type tunneling machine in the horizontal plane.
The curve equation of any point on the section of the shield shell is shown as the formula (9):
Figure BDA0003159193790000072
conversion to polar equation (10):
Figure BDA0003159193790000073
in the formula, R c Is a shield shellThe radius is uniform, and theta is an included angle between a connecting line of a point on the shield shell and the circle center and the polar axis;
displacement of any point on the shield section
Figure BDA0003159193790000074
Can be determined by equations (11) and (12):
Figure BDA0003159193790000075
Figure BDA0003159193790000076
wherein r is the coordinate of the polar axis;
② load around shield shell
Dispersing the shield shell, wherein the shield shell comprises m units and n units along the length direction and the circumferential direction respectively, and the soil pressure applied to any point on the cross section of the shield shell can be determined by the formula (13):
Figure BDA0003159193790000077
wherein p represents the p-th cell in the longitudinal direction, q represents the q-th cell in the circumferential direction, σ v,pq To node vertical pressure, K v,pq Is the vertical soil pressure coefficient, sigma v0,pq To initial vertical pressure, K h,pq Is the horizontal soil pressure coefficient, sigma h,pq To node horizontal pressure, σ h0,pq Is an initial horizontal pressure;
solving for attitude adjusting moment
The gravity of the shield type tunneling machine is G, and the eccentricity is l s And then:
M 1v =G·l s (14)
the formation reaction force and the moment thereof are as follows:
Figure BDA0003159193790000081
Figure BDA0003159193790000082
the cutter head load moment is as follows:
M 4v =FM 1 +FM 2 (17)
in the formula FM 1 Moment caused by the penetration resistance of the cutter head; FM 2 The moment caused by the stratum lateral soil pressure on the cutter head panel is shown as D, and the average diameter of the shield shell is shown as D.
The corrective moment solving process is shown in FIG. 5. And obtaining a deviation correcting angle according to the planned deviation correcting path, calculating the displacement of each point on the shield type tunneling machine according to a formula (12), calculating the load around the shield shell according to a formula (13), and judging whether the load meets a force balance equation (6). If the correction torque is satisfied, determining the correction torque according to the torque balance equation (7). If not, correcting the deviation correcting angle of the shield type tunneling machine, and then carrying out the process again until the deviation correcting moment is calculated.
(5) A hydraulic propulsion system of a shield type heading machine.
And setting the hydraulic propulsion system according to the deviation correcting moment to finish the propulsion and deviation correction of the shield type tunneling machine.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (4)

1. A method for controlling the posture of a shield-type heading machine in a complex stratum in real time is characterized by comprising the following steps:
measuring the posture and the position of a shield type tunneling machine, and judging whether the shaft line offset s of the shield type tunneling machine reaches an allowable offset [ s ], namely judging whether posture correction is needed;
step two, when deviation correction is needed, according to the current position protectionPlanning a deviation rectifying path by using the shaft line offset s and the posture of the shield type tunneling machine, determining single-ring deviation rectifying amount and single-ring deviation rectifying angle, and calculating required deviation rectifying torque; determining a deviation rectifying path by adopting a double-arc interpolation method, and determining a deviation rectifying radius R according to a formula (1) rec (ii) a Calculating the deviation-correcting mileage L and the central angle theta of the double arcs according to the equation set (2) 1 And theta 2 (ii) a Determining the number k of segment rings required by rectification according to a formula (3); determining a curve equation of the double-arc deviation rectifying path according to a formula (4), so as to obtain the coordinates of the shield type tunneling machine at each ring position on the deviation rectifying path, wherein the difference of longitudinal coordinates at two adjacent ring positions is the deviation of each ring, and the included angle between the connecting line at the two adjacent ring positions and the design axis is the deviation rectifying angle of each ring;
Figure FDA0003684765470000011
Figure FDA0003684765470000012
k=[L/g]+1 (3)
Figure FDA0003684765470000013
in the formula, l is the length of the shield type tunneling machine, and d is the wedge-shaped amount of the duct piece; eta is the included angle between the axis of the shield type tunneling machine and the design axis, and g is the width of the single-ring duct piece; s i For the shaft line offset of the shield type tunneling machine, a vertical plane i is equal to v, and a horizontal plane i is equal to h; (x) 1 ,z 1 ) And (y) 1 ,z 1 ) Is the center O of a first section of circular arc on the horizontal and vertical surfaces 1 Determining a deviation rectifying path consisting of two sections of circular arcs by adopting a double circular arc interpolation method, and respectively representing the deviation rectifying path by AQ and QB;
and step three, correspondingly adjusting the hydraulic propulsion system, and performing deviation correction work on the shield type tunneling machine.
2. The method for controlling the attitude of the shield-type heading machine for the complex stratum according to claim 1, which is characterized in that:
in the first step, the distance s of the central point of the shield type tunneling machine deviating from the design axis on the vertical plane and the horizontal plane of the current position is obtained by utilizing the automatic attitude measuring device of the shield type tunneling machine i And an included angle eta between the axis of the shield type tunneling machine and the designed axis, wherein the vertical plane i is equal to v, and the horizontal plane i is equal to h.
3. The method for controlling the attitude of the shield-type heading machine for the complex stratum according to claim 1, which is characterized in that:
in the second step, the calculation method of the deviation rectifying moment comprises the following steps:
1) calculating the section displacement of the shield-type tunneling machine and the load around the shield shell according to the single-ring deviation correction angle beta of the shield-type tunneling machine in the vertical plane, the single-ring deviation correction angle alpha in the horizontal plane and the initial single-ring deviation correction quantity delta, substituting the formula (6), and judging whether the force balance condition is met;
2) if so, determining the correction torque according to the formula (7);
3) if not, after correcting the deviation correction angle of the shield type tunneling machine, repeating the steps 1) and 2) until the deviation correction torque is calculated
Figure FDA0003684765470000021
Figure FDA0003684765470000022
In the formula (I), the compound is shown in the specification,
Figure FDA0003684765470000023
a vertical adjusting moment for controlling the vertical pitch angle generated by the jack thrust difference,
Figure FDA0003684765470000024
is a lifting jackA horizontal adjusting moment generated by the difference of the pushing force and used for controlling the horizontal swinging angle; f 1v Is the gravity of a shield type tunneling machine; f 5v Vertical load resultant force around the shield shell is caused by attitude change; f 5h The resultant force of horizontal load around the shield shell is caused by the change of the posture; m 1v The moment is caused by the eccentric distance between the gravity center and the centroid of the shield type tunneling machine; m 4v The moment is a vertical moment caused by the front load of a cutter head of the shield type tunneling machine; m 5v Vertical resultant moment generated by load around the shield shell due to attitude change; m 5h The horizontal resultant moment generated by the load around the shield shell is caused by the change of the posture.
4. The method for controlling the attitude of the shield-type heading machine for the complex stratum according to claim 3, which is characterized in that:
solving for F 5v 、F 5h 、M 1v 、M 4v 、M 5h And M 5v The calculation steps are as follows:
shield shell section centroid displacement delta s caused by attitude change v And Δ s h Calculated according to equation (8):
Figure FDA0003684765470000025
in the formula,. DELTA.s v Is the vertical displacement of the centroid of the section of the shield shell, and delta s is the vertical displacement of the centroid caused by the gravity action of the shield-type tunneling machine v Adding vertical displacement to the centroid caused by the action of attitude adjusting moment, wherein z is a coordinate along the length direction of the shield type tunneling machine, L is the length of the shield type tunneling machine, and delta s h Is the horizontal displacement of the centroid of the shield shell section, delta h Adding horizontal displacement to a centroid caused by the action of the attitude adjusting moment, wherein beta is a single-ring deviation rectifying angle of the shield type tunneling machine in a vertical plane, and alpha is a single-ring deviation rectifying angle of the shield type tunneling machine in a horizontal plane;
the curve equation of any point on the section of the shield shell is shown as the formula (9):
Figure FDA0003684765470000026
converting equation (9) to polar equation (10):
Figure FDA0003684765470000031
in the formula, R c The mean radius of the shield shell is shown, and theta is an included angle between a connecting line of a point on the shield shell and the circle center and a polar axis;
displacement of any point on the shield section
Figure FDA0003684765470000032
Can be determined by equations (11) and (12):
Figure FDA0003684765470000033
Figure FDA0003684765470000034
wherein r is the coordinate of the polar axis;
determining all-line stratum distribution according to the geological survey report, acquiring engineering mechanical characteristics of each geotechnical layer, and determining parameters a and K in a foundation reaction force curve (5) of each layer of soil 0 ,K min And K max
Figure FDA0003684765470000035
Wherein K is the soil pressure coefficient, i is v or h, which respectively represent the vertical direction and the horizontal direction, K i0 Is the coefficient of static soil pressure, K imin Is the coefficient of active earth pressure, a i Is a foundation reaction coefficient, U i For displacement of the soil mass in the direction i, K imax Is the passive soil pressure coefficient;
and (3) loading around the shield shell: dispersing the shield shell, wherein the shield shell respectively comprises m units and n units along the length direction and the circumferential direction,
the earth pressure at any point on the shield section can be determined by equation (13):
Figure FDA0003684765470000036
wherein p represents the p-th cell in the longitudinal direction, q represents the q-th cell in the circumferential direction, σ v,pq To node vertical pressure, K v,pq Is the vertical soil pressure coefficient, sigma v0,pq To initial vertical pressure, K h,pq Is the horizontal soil pressure coefficient, sigma h,pq To node horizontal pressure, σ h0,pq Is an initial horizontal pressure;
solving the attitude adjusting moment: the gravity of the shield type tunneling machine is G, and the eccentricity is l s And then:
M 1v =G·l s (14)
the formation reaction force and the moment thereof are as follows:
Figure FDA0003684765470000041
Figure FDA0003684765470000042
the cutter head load moment is as follows:
M 4v =FM 1 +FM 2 (17)
in the formula FM 1 Moment caused by the penetration resistance of the cutter head; FM 2 Moment caused by stratum lateral soil pressure on a cutter head panel; d is the average diameter of the shield shell.
CN202110785850.7A 2021-07-12 2021-07-12 Real-time control method for attitude of shield type heading machine in complex stratum Active CN113404499B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110785850.7A CN113404499B (en) 2021-07-12 2021-07-12 Real-time control method for attitude of shield type heading machine in complex stratum

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110785850.7A CN113404499B (en) 2021-07-12 2021-07-12 Real-time control method for attitude of shield type heading machine in complex stratum

Publications (2)

Publication Number Publication Date
CN113404499A CN113404499A (en) 2021-09-17
CN113404499B true CN113404499B (en) 2022-08-02

Family

ID=77686024

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110785850.7A Active CN113404499B (en) 2021-07-12 2021-07-12 Real-time control method for attitude of shield type heading machine in complex stratum

Country Status (1)

Country Link
CN (1) CN113404499B (en)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4421027B2 (en) * 1999-09-22 2010-02-24 株式会社 渡守建設 Measuring device, measuring method, propulsion trajectory management device, and propulsion trajectory management method for propulsion trajectory and propulsion posture in propulsion shield method
CN106246186B (en) * 2016-08-26 2018-08-03 北京中煤矿山工程有限公司 A kind of shaft excavation machine is oriented to control and method of adjustment
JP6316493B1 (en) * 2017-09-26 2018-04-25 大成建設株式会社 Simultaneous drilling assembly control device for shield tunneling machine, simultaneous drilling assembly control system, and simultaneous drilling assembly control method
CN108868807B (en) * 2018-09-07 2019-11-08 上海隧道工程有限公司 The intelligent control method of shield driving correction
CN110067566B (en) * 2019-05-30 2020-06-05 上海隧道工程有限公司 Method and system for predicting shield deviation rectifying moment
CN111271071A (en) * 2020-01-19 2020-06-12 浙江中创天成科技有限公司 Shield tunneling machine attitude control method based on fuzzy adaptive neural network

Also Published As

Publication number Publication date
CN113404499A (en) 2021-09-17

Similar Documents

Publication Publication Date Title
CN111271071A (en) Shield tunneling machine attitude control method based on fuzzy adaptive neural network
CN105909270B (en) A kind of shield machine axis control system
CN109779649A (en) Shield driving axis real-time deviation correcting system and method based on big data
CN110067568B (en) Self-adaptive control method and system for shield deviation-correcting oil pressure output
CN110067566B (en) Method and system for predicting shield deviation rectifying moment
CN112066955B (en) Method and system for measuring pose parameters of body of underground dynamic heading machine
JPH041396A (en) Operating method of small aperture pipe underground-excavator
CN113738390B (en) Time-sharing migration space attitude fine-adjustment deviation-rectifying method of rectangular pipe jacking machine
EP4354247A1 (en) Path planning method for wall-climbing robot
CN110039548B (en) Control method, device and equipment for assembling machine
CN116066123A (en) Automatic tracking control method for shield tunneling track based on model predictive control
CN109766610B (en) Cutter head design method for full-section profiling excavation of rounded rectangular tunnel
CN113404499B (en) Real-time control method for attitude of shield type heading machine in complex stratum
CN112101621A (en) Discrete variable-based soft soil stratum shield construction surface deformation prediction method
CN114991800A (en) Deviation correction control method and system for heading machine, heading machine and storage medium
CN112325887A (en) Method and system for verifying cutting track of boom-type heading machine
CN106703823A (en) Posture error correction system and method of large driving equipment
CN105571811B (en) The method for measuring aerocraft real angle of attack value in wind tunnel experiment
CN111101955A (en) Construction method for ultra-large diameter shield to penetrate through small-radius curve tunnel section
CN111460614A (en) Underground-moon L2 point transfer orbit midway correction method
CN114819256A (en) Continuous real-time trajectory planning method for backhoe excavator
CN103837115A (en) Three-dimensional attitude measurement method and device
CN110188947B (en) Method and system for predicting current ring target in shield deviation correction
CN116108701B (en) FAST novel feed cabin mechanism kinematics positive solution solving and control method
CN115562275A (en) Intelligent navigation control method for coal mine crawler-type heading machine based on MLRNN-PID algorithm

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant