CN117848276A - Large deformation monitoring method suitable for rough surface - Google Patents
Large deformation monitoring method suitable for rough surface Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000012544 monitoring process Methods 0.000 title claims abstract description 15
- 238000005192 partition Methods 0.000 claims abstract description 21
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 238000013461 design Methods 0.000 claims description 27
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 238000000638 solvent extraction Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 description 7
- 238000007405 data analysis Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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
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.
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) |
-
2023
- 2023-12-26 CN CN202311805655.1A patent/CN117848276A/en active Pending
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 |