CN109556539B - Large-gradient tunnel axis detection method based on shield segment determination - Google Patents
Large-gradient tunnel axis detection method based on shield segment determination Download PDFInfo
- Publication number
- CN109556539B CN109556539B CN201811640233.2A CN201811640233A CN109556539B CN 109556539 B CN109556539 B CN 109556539B CN 201811640233 A CN201811640233 A CN 201811640233A CN 109556539 B CN109556539 B CN 109556539B
- Authority
- CN
- China
- Prior art keywords
- plane
- coordinate
- fitting
- coordinates
- formula
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/26—Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
- G01B11/27—Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes for testing the alignment of axes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
The invention relates to a method for detecting the axis of a large-gradient tunnel based on shield segment measurement. Compared with the prior art, the method provided by the invention fully considers the correction of the vertical deviation of the acquired section during the tunnel axis detection, can truly reflect the construction deviation of the tunnel axis, and is more scientific and feasible, small in measurement error and low in implementation cost.
Description
Technical Field
The invention relates to an axis detection method, in particular to a large-gradient tunnel axis detection method based on shield segment determination.
Background
In the process of tunneling the steep-slope tunnel, the front part of the shield tunneling machine has large mass, once attitude deviation is generated, the ring piece is corrected to cause great difficulty, the tunnel forming is not favorably controlled, and even the shield tunneling is influenced, so that the tunnel axis needs to be detected in real time in the tunnel construction process to ensure the accuracy of shield tunneling. The general axis detection method has two methods:
1. measuring the height of a horizontal ruler: the observation point position of the method is not the center of the actual ring piece, particularly when the gradient of the tunnel is large, the gradient cannot be judged by human eyes in the tunnel, and if the height of a vertical surface is used for replacing the radial length of the ring piece, namely the vertical surface of the tunnel, obvious deviation is generated, and the actual measurement precision of the axis is directly influenced.
2. Measuring the section of the total station: the method usually collects three-dimensional coordinates of ring points in a vertical plane, when the gradient of a tunnel is larger, larger deviation between an actual collection plane and a plane to be collected exists, even though some section collection software considers that the collection section direction is determined according to the design gradient of the tunnel axis, the deviation between the actual construction axis and the design axis needs to be measured just due to axis detection, obviously, when the construction axis and the design axis have larger deviation, the true deviation of the construction axis is difficult to find due to orientation according to the design axis.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a method for detecting the axis of a large-gradient tunnel based on shield segment measurement.
The purpose of the invention can be realized by the following technical scheme:
a method for detecting the axis of a large-gradient tunnel based on shield segment determination specifically comprises the following steps:
(1) the total station is oriented vertically to the longitudinal direction of the duct piece alongMeasuring three-dimensional coordinates of the edge of the ring piece at one side of the instrument at the joint of the ring pieces point by point, wherein the three-dimensional coordinates of the measuring points are (x)i,yi,hi)TFitting a plane;
(2) calculating the vertical distance from the measuring points to the fitting plane, eliminating rough difference points, and calculating the coordinates of the qualified measuring points on the fitting plane;
(3) fitting a circular curve by using two-dimensional points on a fitting plane, and calculating a circle center coordinate;
(4) and converting the coordinates of the circle center on the fitting plane into the coordinates of the original coordinate system, and comparing the coordinates with the designed coordinates to finish the detection.
Further, in the step (1): the measuring point is at an equal distance from the seam edge.
Further, in the step (1): the equation of the point fitting space plane o '-x' y 'z' is ax + by + ch + d is 0.
Further, the coordinate of the measuring point in the step (2) is (xp)i,ypi,hpi)TThe calculation equation is:
Further, in the step (3):
the equation for the circle is: (xp)i-x0)2+(ypi-y0)2=R2In the formula: (x)0,y0) As the center plane coordinate, R is the radius of the circle,
error equation:
in the formula:
further, in the step (4): coordinates (x) with the center of the circle in the fitting plane0,y00) to coordinates (x) in the original planex,yy,hh),
In the formula: (x)c,yc,hc) As translation amount of coordinate system, alpha, beta, gamma as transformation rotation angle of coordinate system, R1(α),R2(β),R3(γ)Is a rotation matrix, and the calculation formula is:
the method skillfully utilizes the ring piece joint to determine the tunnel section collecting plane, can objectively reflect the current tunnel construction situation compared with the traditional measuring scheme, establishes an independent coordinate system, fits a space circle by least square, and has the advantages of strict coordinate conversion model, high calculation precision and small error.
Compared with the prior art, the invention has the following beneficial effects:
1) the correction of vertical deviation of the collected section in tunnel axis detection is fully considered, and the construction deviation of the tunnel axis is truly reflected;
2) the field operation implementation scheme is simple and convenient, and the method is scientific and feasible;
3) the same station can detect a plurality of tunnel sections at one time;
4) the measurement precision can be effectively ensured, and the influence of measurement errors is small;
5) and an independent coordinate system is established, so that the calculation is more simplified, the calculation precision is high, and the error is small.
Drawings
FIG. 1 is a schematic diagram of the detection method of the present invention;
FIG. 2 is a schematic illustration of a fitted plane and an original plane;
FIG. 3 is a schematic diagram of an embodiment.
Detailed Description
The invention is described in detail below with reference to the figures and specific embodiments.
Examples
A method for detecting the axis of a large-gradient tunnel based on shield segment determination comprises the following specific scheme:
1) and (3) arranging measurement stations on the ground of the tunnel, and measuring the plane coordinates and the elevation of the measurement stations by using underground control points, as shown in figure 1.
2) And a high-precision prism-free total station is arranged at the measuring station, and the measuring instrument is high.
3) The total station is oriented in the longitudinal direction of the duct piece, and the three-dimensional coordinate (x) is measured point by point along the edge of the ring piece at one side of the instrument at the joint of the ring piecesi,yi,hi)TAnd the measured point is required to be separated from the edge of the seam by the same distance as much as possible (the human eye estimation error is within 5 cm), and the observation data is automatically recorded in the instrument memory.
4) Fitting the spatial plane o '-x' y 'z' with three-dimensional coordinates acquired by a station, as shown in FIG. 2:
ax+by+ch+d=0。
5) calculating the vertical distance from each measuring point to the plane point by point, eliminating the rough difference points, projecting all qualified points to the fitted plane, and calculating the coordinates of the points on the fitted plane: (xp)i,ypi,hpi)T
6) To fit two-dimensional points (xp) on a planei,ypi)TFitting a circular curve by using a least square method and calculating the three-dimensional center of a circleCoordinates are as follows:
the equation for the circle is: (xp)i-x0)2+(ypi-y0)2=R2In the formula: (x)0,y0) As plane coordinates of the center of a circle
Error equation:
in the formula:
7) and (3) coordinate conversion: coordinates (x) centered on the o '-x' y 'z' coordinate system0,y00) to coordinates (x) in the tunnel construction coordinate system o-xyhx,yy,hh)
In the above formula, α, β, γ are coordinate system conversion rotation angles, R1(α),R2(β),R3(γ)Is a rotation matrix, and the calculation formula is as follows:
as shown in fig. 3, the designed gradient of a certain tunnel section is 35 per mill, the inner diameter of the ring piece is 10.5m, two times of section data acquisition of the same measuring station are respectively carried out by adopting a conventional section acquisition method and the method, and the observation data and the calculation result are as follows:
comparing the calculated results without considering the influence of other measuring factors on the deviation, and obtaining the following conclusion:
the vertical deviation of the collected section causes the difference between the center points of the tunnel in the range of the axial line to be 183mm and the difference in the elevation of the axial line to be 7.5mm, and the section line is changed from a circle to an ellipse by collecting the section of the tunnel in the vertical plane, so that the difference of the calculation results of the fitting radius is 4.0 mm.
Usually the tunnel axis measurement error is required to be less than 25mm, and obviously, the above-mentioned multiple deviations will have obvious influence on the plane and elevation detection results of the tunnel axis.
The embodiments described above are described to facilitate an understanding and use of the invention by those skilled in the art. It will be readily apparent to those skilled in the art that various modifications to these embodiments may be made, and the generic principles described herein may be applied to other embodiments without the use of the inventive faculty. Therefore, the present invention is not limited to the above embodiments, and those skilled in the art should make improvements and modifications within the scope of the present invention based on the disclosure of the present invention.
Claims (2)
1. A large-gradient tunnel axis detection method based on shield segment determination is characterized by comprising the following specific steps:
(1) the total station is oriented in the longitudinal direction of the duct piece, and measures three-dimensional coordinates point by point along the edge of the duct piece at one side of the apparatus at the joint of the duct pieces, wherein the three-dimensional coordinates of the measuring points are (x)i,yi,hi)TFitting a plane;
(2) calculating the vertical distance from the measuring points to the fitting plane, eliminating rough difference points, and calculating the coordinates of the qualified measuring points on the fitting plane;
(3) fitting a circular curve by using two-dimensional points on a fitting plane, and calculating a circle center coordinate;
(4) converting the coordinate of the circle center on the fitting plane into the coordinate of the original coordinate system, and comparing the coordinate with the designed coordinate to finish detection;
in the step (1): the measuring points are away from the edge of the seam by the same distance;
the equation of the measuring point fitting space plane o '-x' y 'z' is ax + by + ch + d is 0;
the coordinate of the measuring point in the step (2) is (xp)i,ypi,hpi)TThe calculation equation is:
In the step (4): coordinates (x) with the center of the circle in the fitting plane0,y00) to coordinates (x) in the original planex,yy,hh),
In the formula: (x)c,yc,hc) As translation amount of coordinate system, alpha, beta, gamma as transformation rotation angle of coordinate system, R1(α),R2(β),R3(γ)Is a rotation matrix, and the calculation formula is:
2. the method for detecting the axis of the large-gradient tunnel based on shield segment determination as claimed in claim 1, wherein in the step (3):
the equation for the circle is: (xp)i-x0)2+(ypi-y0)2=R2In the formula: (x)0,y0) As the center plane coordinate, R is the radius of the circle,
error equation:
in the formula:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811640233.2A CN109556539B (en) | 2018-12-29 | 2018-12-29 | Large-gradient tunnel axis detection method based on shield segment determination |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811640233.2A CN109556539B (en) | 2018-12-29 | 2018-12-29 | Large-gradient tunnel axis detection method based on shield segment determination |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109556539A CN109556539A (en) | 2019-04-02 |
CN109556539B true CN109556539B (en) | 2021-07-02 |
Family
ID=65872010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811640233.2A Active CN109556539B (en) | 2018-12-29 | 2018-12-29 | Large-gradient tunnel axis detection method based on shield segment determination |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109556539B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111489390B (en) * | 2020-03-25 | 2023-10-13 | 山东科技大学 | Central axis extraction method for circular tunnel |
CN115248949B (en) * | 2022-09-22 | 2023-02-03 | 中交第一航务工程局有限公司 | End steel shell linearity evaluation method and system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101685012A (en) * | 2008-09-27 | 2010-03-31 | 上海市机械施工有限公司 | Measuring and calculating method of subway tunnel section |
CN103900528A (en) * | 2012-12-28 | 2014-07-02 | 朱志洁 | Detection method for three-dimensional space section roundness of large component |
CN104482918A (en) * | 2014-11-01 | 2015-04-01 | 中铁一局集团有限公司 | Metro shield tunnel forming tunnel portal measurement method |
CN106761797A (en) * | 2016-12-24 | 2017-05-31 | 中铁十八局集团第四工程有限公司 | A kind of heavy grade curve shield driving vertical missing control and the method for correction |
CN107167112A (en) * | 2017-05-24 | 2017-09-15 | 中国葛洲坝集团三峡建设工程有限公司 | A kind of device and method suitable for linear pattern canal axes positioning measurement |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6946834B2 (en) * | 2001-06-01 | 2005-09-20 | Koninklijke Philips Electronics N.V. | Method of orienting an axis of magnetization of a first magnetic element with respect to a second magnetic element, semimanufacture for obtaining a sensor, sensor for measuring a magnetic field |
DE102004010114A1 (en) * | 2004-02-27 | 2005-09-22 | Schlepütz, Roland | Laser-based measurement equipment for positioning e.g. pipe driving machine, has computer to automatically inquire angles and lengths of machine, where angles and lengths are used for positional computation of machine |
CN202382742U (en) * | 2011-11-09 | 2012-08-15 | 苏京平 | Prism-free measurement chip for total station |
CN104389615B (en) * | 2014-09-30 | 2016-08-24 | 中铁四局集团第二工程有限公司 | Littoral soft soil stratum sharp radius curve heavy grade shield construction method |
CN104537162B (en) * | 2014-12-16 | 2017-07-11 | 上海交通大学 | Determine shield tunnel liner interannular seam resistance faulting of slab ends and open deformability method |
JP6484140B2 (en) * | 2015-08-05 | 2019-03-13 | 本田技研工業株式会社 | Vehicle side collision detection assembly |
CN106949886B (en) * | 2017-04-28 | 2023-10-27 | 中国电建集团铁路建设有限公司 | Rapid positioning device and positioning method for subway tunnel engineering line measuring point positions |
CN207228180U (en) * | 2017-08-18 | 2018-04-13 | 上海山南勘测设计有限公司 | A kind of horizontal displacement measuring system based on collimation line method |
CN108614944A (en) * | 2018-05-10 | 2018-10-02 | 西安电子科技大学 | A kind of shield track axis correction parameter prediction technique |
CN109029277A (en) * | 2018-06-27 | 2018-12-18 | 常州沃翌智能科技有限公司 | A kind of tunnel deformation monitoring system and method |
-
2018
- 2018-12-29 CN CN201811640233.2A patent/CN109556539B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101685012A (en) * | 2008-09-27 | 2010-03-31 | 上海市机械施工有限公司 | Measuring and calculating method of subway tunnel section |
CN103900528A (en) * | 2012-12-28 | 2014-07-02 | 朱志洁 | Detection method for three-dimensional space section roundness of large component |
CN104482918A (en) * | 2014-11-01 | 2015-04-01 | 中铁一局集团有限公司 | Metro shield tunnel forming tunnel portal measurement method |
CN106761797A (en) * | 2016-12-24 | 2017-05-31 | 中铁十八局集团第四工程有限公司 | A kind of heavy grade curve shield driving vertical missing control and the method for correction |
CN107167112A (en) * | 2017-05-24 | 2017-09-15 | 中国葛洲坝集团三峡建设工程有限公司 | A kind of device and method suitable for linear pattern canal axes positioning measurement |
Also Published As
Publication number | Publication date |
---|---|
CN109556539A (en) | 2019-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107167790B (en) | A kind of two step scaling method of laser radar based on Calibration Field | |
US10422619B2 (en) | Identification of geometric deviations of a motion guide in a coordinate-measuring machine or in a machine tool | |
CN100562708C (en) | Six-freedom degree articulated flexible COMERO high precision calibration method | |
WO2015078874A1 (en) | Calibration apparatus and method for computed tomography | |
CN109186487A (en) | A kind of pipeline ovality automatic checkout equipment and its detection method | |
CN109612406B (en) | Random detection method for ring splicing quality of splicing pipe pieces of shield tunnel | |
CN109556539B (en) | Large-gradient tunnel axis detection method based on shield segment determination | |
CN108311545B (en) | Y-type rolling mill continuous rolling centering and hole pattern detection system and method | |
CN103386640B (en) | Large caliber reflecting mirror machining tool accuracy alignment method | |
CN102494657A (en) | Measuring head radius compensation method for curve surface profile measuring and detecting | |
CN104101296A (en) | Digital positioning detection method and system in large-scale structure precision assembly | |
CN104070418B (en) | Axial symmetry surface shape of optical aspheric surface On-line Measuring Method | |
CN103791868B (en) | A kind of space nominal volume based on virtual ball and scaling method thereof | |
CN103735282A (en) | Device and method for geometric correction of detector of cone-beam CT (computed tomography) system | |
CN106017319A (en) | Coordinate transformation tool and method of 3D laser scanning data based on high-precision point measurement | |
CN108317999B (en) | Communication iron tower inclination angle measuring method based on inertial sensor | |
CN103900490A (en) | Device and method for detecting ovality of tyre | |
CN110017823B (en) | Rail positioning coordinate calculation method and device based on total station instrument free-standing station | |
WO2018099051A1 (en) | Method for verifying accuracy of rtk tilt compensation measurement | |
CN105043381A (en) | Magnetic nail-based positioning method | |
CN202393362U (en) | Long axis straightness detection system | |
CN101013065A (en) | Pixel frequency based star sensor high accuracy calibration method | |
CN106052568A (en) | Compensation method for screw type surface detection device displacement sensor system mounting error | |
CN106989670A (en) | A kind of non-contact type high-precision large-scale workpiece tracking measurement method of robot collaboration | |
CN105559809B (en) | Scanning method and device |
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 |