CN117848276A - Large deformation monitoring method suitable for rough surface - Google Patents

Large deformation monitoring method suitable for rough surface Download PDF

Info

Publication number
CN117848276A
CN117848276A CN202311805655.1A CN202311805655A CN117848276A CN 117848276 A CN117848276 A CN 117848276A CN 202311805655 A CN202311805655 A CN 202311805655A CN 117848276 A CN117848276 A CN 117848276A
Authority
CN
China
Prior art keywords
scanning
deformation
value
deltay
curved surface
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.)
Pending
Application number
CN202311805655.1A
Other languages
Chinese (zh)
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.)
Shanghai Construction Group Co Ltd
Original Assignee
Shanghai Construction Group Co Ltd
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 Shanghai Construction Group Co Ltd filed Critical Shanghai Construction Group Co Ltd
Priority to CN202311805655.1A priority Critical patent/CN117848276A/en
Publication of CN117848276A publication Critical patent/CN117848276A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

The invention relates to a large deformation monitoring method suitable for a rough surface. First, three-dimensional scanning: the reference setting is that four target balls are fixed near four corners of the ground connection wall body to be measured; and (3) station arrangement scanning: starting from the edge of one side of the wall body, arranging a three-dimensional scanning instrument at the position 5-8 m in front of the wall body for scanning, and after the secondary scanning is finished, carrying out station arrangement scanning at intervals of 8-10 m along the direction of the wall body until the other side of the wall body is reached; and (3) splicing: summarizing and splicing the multi-station scanning data through software self-matched with the system to form a current scanning point cloud model of the whole wall; secondly, in view of deviation of the instrument erection position during each scanning, the whole wall coordinate unification work is completed through unifying the target ball coordinates; then, carrying out gridding partition division on the point cloud data scanned at the time according to the corrected coordinates; then, the subareas perform surface fitting: performing surface fitting on each partition point cloud surface by adopting a quadric surface fitting method; finally, deformation analysis.

Description

Large deformation monitoring method suitable for rough surface
Technical Field
The invention belongs to the technical field of underground engineering informatization measurement, and particularly relates to a large deformation monitoring method suitable for a rough surface.
Background
At present, a prism or a target is arranged on a wall body at fixed points for identifying the deformation of the structure of the underground continuous wall, the position change of the prism or the target is judged through a total station, the deformation of the wall body is indirectly judged, and the method is often in a missing state and cannot find dangerous points in time. In order to solve the defects of few measuring points and low measuring efficiency of the traditional measuring method, the integrated measurement of the underground continuous wall by adopting three-dimensional scanning becomes the current development direction. Tens of millions of point cloud data can be obtained in a few minutes through three-dimensional scanning, and deformation information of the whole wall body can be rapidly and comprehensively obtained.
After the enclosure deformation data is acquired, the analysis of the deformation of the diaphragm wall is required. However, under the current construction process conditions, the surface of the underground continuous wall is often mottled and uneven, so that the measurement result of the adjacent point cloud data is discontinuous, and a large amount of deviation of the measurement result caused by different point taking positions of two measurements is very easy; the long-time open air environment can further cause the surface to be damaged or clay, the point cloud data analysis is directly carried out, a large number of large deformation points can be generated suddenly, unnecessary alarms are caused, and the ground wall deformation judgment analysis is seriously affected.
Therefore, how to improve the existing data analysis method and find a method for monitoring deformation of an open air environment, which is suitable for the rough surface of a foundation pit underground wall, is a technical problem to be solved by those skilled in the art.
Disclosure of Invention
The invention aims to provide a large deformation monitoring method suitable for a rough surface, which avoids measurement errors caused by discontinuous monitoring data through quadric surface fitting, improves the accuracy of the monitoring data, realizes the accurate calculation of the large deformation of the whole underground continuous wall, and timely and comprehensively discovers a large deformation area of a surrounding.
In order to solve the technical problems, the invention comprises the following technical scheme:
a method for monitoring large deformations suitable for rough surfaces, comprising the steps of:
step S1: three-dimensional scanning:
(1) The reference setting comprises the steps of fixing four target balls near four corners of a ground connection wall to be measured, wherein the target balls can be fixed by a magnetic attraction method; (2) station-laying scanning: starting from the edge of one side of the wall body, arranging a three-dimensional scanning instrument at the position 5-8 m in front of the wall body for scanning, and after the secondary scanning is finished, carrying out station arrangement scanning at intervals of 8-10 m along the direction of the wall body until the other side of the wall body is reached; and (3) splicing: summarizing and splicing the multi-station scanning data through software self-matched with the system to form a current scanning point cloud model of the whole wall;
step S2: unified coordinates:
(1) Coordinate system transformation parameter determination: in view of the deviation of the instrument erection position during each scanning, the coordinates of the same point cloud in each scanning model are different, the coordinate system needing to be scanned for multiple times is unified, the unified work of the whole wall body coordinate can be completed through unifying the target ball coordinate, and the coordinate conversion formula is as follows:
wherein R (α, β, γ) is a rotation matrix, and α, β, γ is rotation angles around the X axis, Y axis, Z axis; to simplify the later analysis, the unified coordinate system marks the target ball K with the lower left corner of the foundation pit 1 The extending direction of the wall body is an X axis, and the vertical direction of the wall body is a Z axis direction;
assuming the nth scan, target sphere K 1 Is (X) K1，n ，Y K1，n ,Z K1，n ) Target ball K for right lower corner of foundation pit 2 Is (X) K2，n ，Y K2，n ,Z K2，n ) Due to K 1 Is the origin, so
As the extending direction of the wall body is X-axis, the Z-axis is not changed in angle, so
(2) Calculating a single target coordinate deviation value of each point cloud model, and assuming that the m data of a certain point cloud in the nth scanning is (X) m，n ，Y m，n ，Z m，n ) The transformed coordinates are:
step S3: point cloud meshing partition attribution: grid partition division is carried out on the point cloud data scanned at the time according to the corrected coordinates, and the grid size is recommended to be 20cm;
step S4: and (5) carrying out surface fitting in a partition: performing surface fitting on each partitioned point cloud surface by adopting a quadric surface fitting method so as to obtain the nth S A×B Partitioning;
step S5: and (5) deformation analysis.
Further, the meshing partition division calculation method in step S3 includes:
first, reference value setting, assuming K 1 Is a partition benchmark;
secondly, partitioning: suppose that the m coordinate of a point cloud in the nth scan is corrected to (X) m，n ，Y m，n ，Z m，n ) The partition to which the m point belongs is S A×B Wherein a=int [ (X) m，n -X K1，n) /0.02]+1，B＝int[(Z m，n -Z K1，n )/0.02]+1
Thus, all the point clouds can be subjected to partition attribution.
Further, the quadric surface fitting method in step S4 includes:
F A×B，n (x，y，z)＝a n x 2 +b n y 2 +c n z 2 +d n xy+e n xz+f n yz+g n x+h n y+i n z+j n wherein a to j are constants, and are calculated by least squares fitting, and when x and z are determined, the cumulative deformation y value is determined in consideration of the fact that the scan point cloud curved surface is positioned on the XOZ plane, so F A×B (x, y, z) can be converted into a quadratic function y with x and z as arguments A×B ＝S A×B (x，z)。
Further, the step S5 includes:
first, cumulative deformation analysis: subtracting the deformed curved surface from the n-time scanning
Initial deformed surface:
S A×B (x，z)＝S A×B，n (x，z)-S A×B，1 (x，z)；
obtaining the middle extremum y of curved surface P A×B ：
Will S A×B (x, z) deriving x and z and letting the derivative be 0, i.e
The x and z values at the extremum of the curved surface can be obtained and substituted into y A×B ＝S A×B (x, z) determination
AB region cumulative extremum y P A×B The value of the sum of the values,
obtaining the edge extremum y of the curved surface B A×B : substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into S A×B (x，z)，
Order the
X and z are substituted into y when extremum occurs on four edges of curved surface A×B ＝S A×B (x, z) finding the extremum y of 4 edges B A×B The value of the sum of the values,
obtaining a curved surface endpoint value y E A×B ：
Substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into S A×B (x, z) substituting S A×B (x, z) to obtain four end point values y E A×B ；
Let y A×B ＝max(|y P A×B |，|y E A×B |，|y B A×B |)，
Will y A×B Value and maximum value Y of envelope deformation in design file max Comparing, and judging whether an alarm is required;
when |y A×B |≤|Y max The deformation of the AB area is satisfied with the requirements of design files,
when |y A×B |＞|Y max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed;
analyzing the deformation of other subareas by adopting the same mode, thereby obtaining an alarm area and a maximum deformation value of the whole enclosure wall;
second, the current variation analysis: subtracting the last deformed curved surface from the current deformed curved surface:
Δy A×B，n ＝Q A×B，n (x，z)＝S A×B，n (x，z)-S A×B，n-1 (x, z) obtaining the curved surface intermediate extremum deltay P A×B : will delta y A×B，n ＝Q A×B，n (x, z) deriving x and z and letting the derivative be 0, i.eThe x and z values at the extremum of the curved surface can be obtained and substituted into deltay A×B，n ＝Q A×B，n (x, z) obtaining the variation extreme value deltay of the nth scan of the AB region P A×B，n ；
Obtaining the edge extremum deltay of the curved surface B A×B ：
Substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into Q A×B，n (x，z)，
Order the
X and z substitution deltay when extremum occurs to four edges of curved surface can be obtained A×B，n ＝Q A×B，n (x, z) obtaining the extreme value deltay of 4 edges B A×B，n ；
Obtaining a curved surface endpoint value delta y E A×B ：
Let x=0.02A-0.02, x=0.02A, z=0.02a-0.02, z=0.02a substituted into S A×B (x, z) to obtain four end points Δy B A×B，n ；
Let Deltay A×B，n ＝max(|Δy P A×B，n |，|Δy E A×B，n |，|Δy B A×B，n |) is provided; scanning frequency is not suitable to be lower than 1 day/time; when the scanning frequency is 1 day/time;
will delta y A×B，n Comparing the value with a maximum value delta Ymax of the deformation of the enclosure day in the design file, and judging whether an alarm is required;
when |Deltay A×B，n |≤|ΔY max The deformation of the AB area enclosure meets the requirements of design files;
when |Deltay A×B，n |＞|ΔY max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed; when the scanning frequency is 1 day multiple times (assuming T times/day);
will |Deltay A×B，n +Δy A×B，n-1 +…+Δy A×B，n-T+1 Comparing the value with the maximum value delta Ymax of the deformation of the enclosure day in the design file, and judging whether an alarm is required;
when |Deltay A×B，n +Δy A×B，n-2 +…+Δy A×B，n-T+1 |≤|ΔY max I, AB area enclosing deformation meets the requirement of design files
When |Deltay A×B，n +Δy A×B，n-1 +…+Δy A×B，n-T+1 |＞|ΔY max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed;
and the deformation of other subareas is analyzed in the same way, so that the current day alarm area and the maximum daily deformation value are obtained.
Compared with the prior art, the invention has the beneficial technical effects that:
according to the large deformation monitoring method suitable for the rough surface, the quadric surface fitting is adopted, so that measurement errors caused by discontinuous monitoring data are avoided, the accuracy of the monitoring data is improved, the accurate calculation of the large deformation of the whole underground continuous wall body is realized, and the large deformation area of the enclosure is found timely and comprehensively.
Detailed Description
The method for monitoring large deformation of a rough surface according to the present invention will be described in further detail with reference to specific examples. The advantages and features of the present invention will become more apparent from the following description. It should be noted that the embodiments of the present invention are illustrated in a very simplified form and are presented in a non-exclusive scale for convenience and clarity.
Example 1
A method for monitoring large deformations suitable for rough surfaces, comprising the steps of: the method comprises the following steps:
step S1: three-dimensional scanning:
(1) The reference setting comprises the steps of fixing four target balls near four corners of a ground connection wall to be measured, wherein the target balls can be fixed by a magnetic attraction method; (2) station-laying scanning: starting from the edge of one side of the wall body, arranging a three-dimensional scanning instrument at the position 5-8 m in front of the wall body for scanning, and after the secondary scanning is finished, carrying out station arrangement scanning at intervals of 8-10 m along the direction of the wall body until the other side of the wall body is reached; and (3) splicing: summarizing and splicing the multi-station scanning data through software self-matched with the system to form a current scanning point cloud model of the whole wall;
step S2: unified coordinates:
(1) Coordinate system transformation parameter determination: in view of the deviation of the instrument erection position during each scanning, the coordinates of the same point cloud in each scanning model are different, the coordinate system needing to be scanned for multiple times is unified, the unified work of the whole wall body coordinate can be completed through unifying the target ball coordinate, and the coordinate conversion formula is as follows:
wherein R (α, β, γ) is a rotation matrix, and α, β, γ is rotation angles around the X axis, Y axis, Z axis; to simplify the later analysis, the unified coordinate system marks the target ball K with the lower left corner of the foundation pit 1 The extending direction of the wall body is an X axis, and the vertical direction of the wall body is a Z axis direction;
assuming the nth scan, target sphere K 1 Is used for sittingMarked as (X) K1，n ，Y K1，n ，Z K1，n ) Target ball K for right lower corner of foundation pit 2 Is (X) K2，n ，Y K2，n ，Z K2，n ) Due to K 1 Is the origin, soAs the extending direction of the wall body is X-axis, the Z-axis is not changed in angle, so +.>
(2) Calculating a single target coordinate deviation value of each point cloud model, and assuming that the m data of a certain point cloud in the nth scanning is (X) m，n ，Y m，n ，Z m，n ) The transformed coordinates are:
step S3: point cloud meshing partition attribution: grid partition division is carried out on the point cloud data scanned at the time according to the corrected coordinates, and the grid size is recommended to be 20cm;
step S4: and (5) carrying out surface fitting in a partition: performing surface fitting on each partitioned point cloud surface by adopting a quadric surface fitting method so as to obtain the nth S A×B Partitioning;
step S5: and (5) deformation analysis.
In this embodiment, more preferably, the meshing partition division calculation method in step S3 includes:
first, reference value setting, assuming K 1 Is a partition benchmark;
secondly, partitioning: suppose that the m coordinate of a point cloud in the nth scan is corrected to (X) m，n ，Y m，n ，Z m，n ) The partition to which the m point belongs is S A×B Wherein a=int [ (X) m，n -X K1，n )/0.02]+1，B＝int[(Z m，n -Z K1，n )/0.02]+1 so that all point clouds can be zone-assigned.
In this embodiment, more preferably, the quadric surface fitting method in step S4 includes:
F A×B，n (x，y，z)＝a n x 2 +b n y 2 +c n z 2 +d n xy+e n xz+f n yz+g n x+h n y+i n z+j n ，
wherein a to j are constants, and are calculated by least squares fitting, and when x and z are determined, the accumulated deformation y value is determined in consideration of the fact that the scan point cloud curved surface is positioned on the XOZ plane, so F is obtained A×B (x, y, z) can be converted into a quadratic function y with x and z as arguments A×B ＝S A×B (x，z)。
In this embodiment, more preferably, the step S5 includes:
first, cumulative deformation analysis: subtracting the initial deformed surface from the n-time scanned deformed surface:
S A×B (x，z)＝S A×B，n (x，z)-S A×B，1 (x, z); curved surface intermediate extremum y P A×B Obtaining:
will S A×B (x, z) deriving x and z and letting the derivative be 0, i.eThe x and z values at the extremum of the curved surface can be obtained and substituted into y A×B ＝S A×B (x, z) determining the cumulative extremum y of the AB region P A×B The value of the sum of the values,
obtaining the edge extremum y of the curved surface B A×B : substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into S A×B (x，z)，
Order the
X and z are substituted into y when extremum occurs on four edges of curved surface A×B ＝S A×B (x, z) obtaining extremum yP of 4 edges A×B The value of the sum of the values,
obtaining a curved surface endpoint value y E A×B ：
Substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into S A×B (x, z) substituting S A×B (x, z) to obtain four end point values y E A×B ；
Let y A×B ＝max(|y P A×B |，|y E A×B |，|y B A×B |)，
Will y A×B Value and maximum value Y of envelope deformation in design file max Comparing, and judging whether an alarm is required;
when |y A×B |≤|Y max The deformation of the AB area is satisfied with the requirements of design files,
when |y A×B |＞|Y max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed;
analyzing the deformation of other subareas by adopting the same mode, thereby obtaining an alarm area and a maximum deformation value of the whole enclosure wall;
second, the current variation analysis: subtracting the last deformed curved surface from the current deformed curved surface:
Δy A×B，n ＝Q A×B，n (x，z)＝S A×B，n (x，z)-S A×B，n-1 (x, z) obtaining the curved surface intermediate extremum deltay P A×B : will delta y A×B，n ＝Q A×B，n (x, z) deriving x and z and making the derivative 0, i.eThe x and z values at the extremum of the curved surface can be obtained and substituted into deltay A×B，n ＝Q A×B，n (x, z) obtaining the variation extreme value deltay of the nth scan of the AB region P A×B，n ；
Obtaining the edge extremum deltay of the curved surface B A×B
Substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into Q A×B，n (x，z)，
Order the
X and z substitution deltay when extremum occurs to four edges of curved surface can be obtained A×B，n ＝Q A×B，n (x, z) obtaining the extreme value deltay of 4 edges B A×B，n ；
Obtaining a curved surface endpoint value delta y E A×B ：
Substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into S A×B (x, z) to obtain four end points Δy B A×B，n ；
Let Deltay A×B，n ＝max(|Δy P A×B，n |，|Δy E A×B，n |，|Δy B A×B，n |)；
Scanning frequency is not suitable to be lower than 1 day/time;
when the scanning frequency is 1 day/time;
will delta y A×B，n Comparing the value with a maximum value delta Ymax of the deformation of the enclosure day in the design file, and judging whether an alarm is required;
when |Deltay A×B，n |≤|ΔY max The deformation of the AB area enclosure meets the requirements of design files;
when |Deltay A×B，n |＞|ΔY max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed;
when the scanning frequency is 1 day multiple times (assuming T times/day);
will |Deltay A×B，n +Δy A×B，n-1 +…+Δy A×B，n-T+1 Comparing the value with the maximum value delta Ymax of the deformation of the enclosure day in the design file, and judging whether an alarm is required;
when |Deltay A×B，n +Δy A×B，n-1 +…+Δy A×B，n-T+1 |≤|ΔY max I, AB area surroundingThe deformation protection meets the requirements of design files
When |Deltay A×B，n +Δy A×B，n-1 +…+Δy A×B，n-T+1 |＞|ΔY max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed;
and the deformation of other subareas is analyzed in the same way, so that the current day alarm area and the maximum daily deformation value are obtained.
The above examples are preferred embodiments of the present invention, but the embodiments of the present invention are not limited to the above examples. The above examples illustrate only a few embodiments of the invention, which are described in detail and are not to be construed as limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention. Accordingly, the scope of protection of the present invention is to be determined by the appended claims.

Claims (4)

1. A method for monitoring large deformations suitable for rough surfaces, comprising the steps of:
step S1: three-dimensional scanning:
(1) The reference setting comprises the steps of fixing four target balls near four corners of a ground connection wall to be measured, wherein the target balls can be fixed by a magnetic attraction method; (2) station-laying scanning: starting from the edge of one side of the wall body, arranging a three-dimensional scanning instrument at the position 5-8 m in front of the wall body for scanning, and after the secondary scanning is finished, carrying out station arrangement scanning at intervals of 8-10 m along the direction of the wall body until the other side of the wall body is reached; and (3) splicing: summarizing and splicing the multi-station scanning data through software self-matched with the system to form a current scanning point cloud model of the whole wall;
step S2: unified coordinates:
(1) Coordinate system transformation parameter determination: in view of the deviation of the instrument erection position during each scanning, the coordinates of the same point cloud in each scanning model are different, the coordinate system needing to be scanned for multiple times is unified, the unified work of the whole wall body coordinate can be completed through unifying the target ball coordinate, and the coordinate conversion formula is as follows:
wherein R (α, β, γ) is a rotation matrix, and α, β, γ is rotation angles around the X axis, Y axis, Z axis; to simplify the later analysis, the unified coordinate system marks the target ball K with the lower left corner of the foundation pit 1 The extending direction of the wall body is an X axis, and the vertical direction of the wall body is a Z axis direction;
assuming the nth scan, target sphere K 1 Is (X) K1，n ，Y K1，n ,Z K1，n ) Target ball K for right lower corner of foundation pit 2 Is (X) K2，n ，Y K2，n ,Z K2，n ) Due to K 1 Is the origin, so
As the extending direction of the wall body is X-axis, the Z-axis is not changed in angle, so
α n ＝β n ＝1
(2) Calculating a single target coordinate deviation value of each point cloud model, and assuming that the m data of a certain point cloud in the nth scanning is (X) m，n ，Y m，n ，Z m，n ) The transformed coordinates are:
step S3: point cloud meshing partition attribution: grid partition division is carried out on the point cloud data scanned at the time according to the corrected coordinates, and the grid size is recommended to be 20cm;
step S4: and (5) carrying out surface fitting in a partition: performing surface fitting on each partitioned point cloud surface by adopting a quadric surface fitting method so as to obtain the nth S A×B Partitioning;
step S5: and (5) deformation analysis.
2. The method according to claim 1, wherein the meshing partition division calculation method in step S3 includes:
first, reference value setting, assuming K 1 Is a partition benchmark;
secondly, partitioning: suppose that the m coordinate of a point cloud in the nth scan is corrected to (X) m,n ，Y m,n ，Z m,n ) The partition to which the m point belongs is S A×B Wherein a=int [ (X) m，n -X K1，n )/0.02]+1，B＝int[(Z m，n -Z K1，n )/0.02]+1 so that all point clouds can be zone-assigned.
3. The method according to claim 2, wherein the quadric surface fitting method in step S4 includes:
F A×B，n (x，y，z)＝a n x 2 +b n y 2 +c n z 2 +d n xy+e n xz+f n yz+g n x+h n y+i n z+j n wherein a to j are constants, and are calculated by least squares fitting, and when x and z are determined, the cumulative deformation y value is determined in consideration of the fact that the scan point cloud curved surface is positioned on the XOZ plane, so F A×B (x, y, z) can be converted into a quadratic function y with x and z as arguments A×B ＝S A×B (x，z)。
4. A method according to claim 3, wherein said step S5 comprises:
first, cumulative deformation analysis: subtracting the initial deformed surface from the n-time scanned deformed surface:
S A×B (x，z)＝S A×B，n (x，z)-S A×B，1 (x，z)；
obtaining the middle extremum y of curved surface P A×B ：
Will S A×B (x, z) deriving x and z and letting the derivative be 0, i.eThe x and z values at the extremum of the curved surface can be obtained and substituted into y A×B ＝S A×B (x, z) determining the cumulative extremum y of the AB region P A×B The value of the sum of the values,
obtaining the edge extremum y of the curved surface P A×B : substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into S A×B (x, z), let
X and z are substituted into y when extremum occurs on four edges of curved surface A×B ＝S A×B (x, z) finding the extremum y of 4 edges B A×B The value of the sum of the values,
obtaining a curved surface endpoint value y E A×B ：
Substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into S A×B (x, z) substituting S A×B (x, z) to obtain four end point values y E A×B ；
Let y A×B ＝max(|y P A×B |，|y E A×B |，|y B A×B |)，
Will y A×B Value and maximum value Y of envelope deformation in design file max Comparing, and judging whether an alarm is required; when |y A×B |≤|Y max The deformation of the AB area is satisfied with the requirements of design files,
when |y A×B |＞|Y max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed;
analyzing the deformation of other subareas by adopting the same mode, thereby obtaining an alarm area and a maximum deformation value of the whole enclosure wall;
second, the current variation analysis: subtracting the last deformed curved surface from the current deformed curved surface: Δy A×B，n ＝Q A×B，n (x，z)＝S A×B，n (xz)-S A×B，n-1 (x, z) obtaining the curved surface intermediate extremum deltay P A×B : will delta y A×B，n ＝Q A×B，n (x, z) deriving x and z and letting the derivative be 0, i.eThe x and z values at the extremum of the curved surface can be obtained and substituted into deltay A×B，n ＝Q A×B，n (x, z) obtaining the variation extreme value deltay of the nth scan of the AB region P A×B，n ；
Obtaining the edge extremum deltay of the curved surface B A×B ：
Substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into Q A×B，n (x，z)，
Order the
X and z substitution deltay when extremum occurs to four edges of curved surface can be obtained A×B，n ＝Q A×B，n (x, z) obtaining the extreme value deltay of 4 edges B A×B，n ；
Obtaining a curved surface endpoint value delta y E A×B ：
Substituting x=0.02a-0.02, x=0.02a, z=0.02a-0.02, z=0.02a into S A×B (x, z) to obtain four end points Δy B A×B，n ；
Let Deltay A×B，n +max(|Δy P A×B，n |，|Δy E A×B，n |，|Δy B A×B，n |) is provided; the scanning frequency is not lowAt 1 day/time; when the scanning frequency is 1 day/time;
will delta y A×B，n Comparing the value with a maximum value delta Ymax of the deformation of the enclosure day in the design file, and judging whether an alarm is required;
when |Deltay A×B，n |≤|ΔY max The deformation of the AB area enclosure meets the requirements of design files;
when |Deltay A×B，n |＞|ΔY max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed;
when the scanning frequency is 1 day multiple times (assuming T times/day);
will |Deltay A×B，n +Δy A×B，n-1 +…+Δy A×B，n-T+1 Comparing the value with the maximum value delta Ymax of the deformation of the enclosure day in the design file, and judging whether an alarm is required;
when |Deltay A×B，n +Δy A×B，n-1 +…+Δy A×B，n-T+1 |＞|ΔY max I, AB area enclosing deformation meets the requirement of design files
When |Deltay A×B，n +Δy A×B，n-1 +…+Δy A×B，n-T+1 |＞|ΔY max I, the AB area enclosure deformation does not meet the requirement of design files, and an alarm is needed;
and the deformation of other subareas is analyzed in the same way, so that the current day alarm area and the maximum daily deformation value are obtained.
CN202311805655.1A 2023-12-26 2023-12-26 Large deformation monitoring method suitable for rough surface Pending CN117848276A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311805655.1A CN117848276A (en) 2023-12-26 2023-12-26 Large deformation monitoring method suitable for rough surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311805655.1A CN117848276A (en) 2023-12-26 2023-12-26 Large deformation monitoring method suitable for rough surface

Publications (1)

Publication Number Publication Date
CN117848276A true CN117848276A (en) 2024-04-09

Family

ID=90533872

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311805655.1A Pending CN117848276A (en) 2023-12-26 2023-12-26 Large deformation monitoring method suitable for rough surface

Country Status (1)

Country Link
CN (1) CN117848276A (en)

Similar Documents

Publication Publication Date Title
US20220018960A1 (en) Method for coordinate error correction with a three-dimensional lidar scanner
CN110793508B (en) Method and device for processing observation data of control points in three-dimensional control network
CN113435757B (en) BIM-based building masonry design and construction method
WO2021098808A1 (en) Method and system for determining laser tracker station, electronic device, and medium
CN104729529B (en) The method and system that map surveying systematic error judges
CN110516350A (en) A kind of ERS point tolerance modification method based on anisotropy weighting
CN111046584B (en) Precise adjustment method for satellite-borne instrument and equipment
CN116713993A (en) Six-axis serial robot D-H parameter calibration method and system
CN111691691A (en) Method for installing, measuring and positioning box body bending and twisting component and drum node
CN111694009B (en) Positioning system, method and device
CN117848276A (en) Large deformation monitoring method suitable for rough surface
CN109855611B (en) PC wall body rapid measurement and calibration method based on total station
CN106323587B (en) The monocular video high-precision measuring method of wing model in wind tunnel flexible deformation
CN110572777B (en) Cellular area restriction-based auxiliary positioning method
CN115289999B (en) Self-calibration method of three-dimensional laser scanning system
CN114001682B (en) Flatness judging method for mounting surface of heavy door frame of electric room module
CN115930800A (en) Tunnel face displacement field monitoring method based on three-dimensional laser point cloud
CN116068535A (en) Multi-laser radar external parameter self-calibration method, system, terminal equipment and storage medium
CN113203369B (en) Method for measuring normal verticality of robot hole making
CN115950352A (en) Iron tower bolt looseness detection and identification method based on laser vibration measurement
CN115619969A (en) Building damage trace measuring method based on laser point cloud difference
CN110285764B (en) Display panel cutting precision measuring method and display panel
CN110866951B (en) Method for correcting optical axis inclination of monocular camera
US20240282111A1 (en) Electronic counter scale intelligent verification method based on deep learning detection and identification
Algadhi et al. Assessment of accuracy and performance of terrestrial laser scanner in monitoring of retaining walls

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