Study on Vertical Alignment Maintenance Technique using GNSS in Skyscraper

Similar documents
1314. Estimation of mode shapes expanded from incomplete measurements

INCREMENTAL DISPLACEMENT ESTIMATION METHOD FOR VISUALLY SERVOED PARIED STRUCTURED LIGHT SYSTEM (ViSP)

Lecture 3: Camera Calibration, DLT, SVD

THE POSITION AND ORIENTATION MEASUREMENT OF GONDOLA USING A VISUAL CAMERA

A STRUCTURAL OPTIMIZATION METHODOLOGY USING THE INDEPENDENCE AXIOM

CS231A Course Notes 4: Stereo Systems and Structure from Motion

Exterior Orientation Parameters

Euclidean Space. Definition 1 (Euclidean Space) A Euclidean space is a finite-dimensional vector space over the reals R, with an inner product,.

The Use of Biplot Analysis and Euclidean Distance with Procrustes Measure for Outliers Detection

Rectification and Distortion Correction

Designing a GPS Receiver Network with GNSS Algorithm for Accuracy and Safety

ctbuh.org/papers Structural Member Monitoring of High-Rise Buildings Using a 2D Laser Scanner Title:

Extension of the ABC-Procrustes algorithm for 3D affine coordinate transformation

calibrated coordinates Linear transformation pixel coordinates

A METHOD TO MODELIZE THE OVERALL STIFFNESS OF A BUILDING IN A STICK MODEL FITTED TO A 3D MODEL

Horizontal Positional Accuracy (HPA) How is it calculated? Gord Gamble and Peter Goodier, Practice Advisory Department

Transducer and LBL calibration - Integrated functions in HiPAP systems

Monocular Vision-based Displacement Measurement System Robust to Angle and Distance Using Homography

Korean Geodetic Datum 2002(KGD2002) : Nationwide GPS Network Densification. Mr. Gwang-Ho Jung, 3) Mr. Sang-Hun Cha. Kyungil University, Rep.

A Study on Development of Azimuth Angle Tracking Algorithm for Tracking-type Floating Photovoltaic System

Chapters 1 7: Overview

Presented at the FIG Working Week 2017, May 29 - June 2, 2017 in Helsinki, Finland

Reminder: Lecture 20: The Eight-Point Algorithm. Essential/Fundamental Matrix. E/F Matrix Summary. Computing F. Computing F from Point Matches

UNIVERSITY CALIFORNIA, RIVERSIDE AERIAL TARGET GROUND CONTROL SURVEY REPORT JOB # DATE: MARCH 2011

Efficient Stereo Image Rectification Method Using Horizontal Baseline

Aaronia GPS Logger Programming Guide

Serially-Linked Parallel Leg Design for Biped Robots

DIGITAL ROAD PROFILE USING KINEMATIC GPS

Epipolar Geometry and the Essential Matrix

Observed Differential Carrier Phase Float Accuracies on three Kinematic Surveys using a u-blox Antaris GPS Receiver under Open Sky Conditions

Subspace Clustering with Global Dimension Minimization And Application to Motion Segmentation

ξ ν ecliptic sun m λ s equator

GI-Eye II GPS/Inertial System For Target Geo-Location and Image Geo-Referencing

THE SIMPLE POLYNOMAIL TECHNIQUE TO COMPUTE THE ORTHOMETRIC HEIGHT IN EGYPT

Use of n-vector for Radar Applications

Navigational Aids 1 st Semester/2007/TF 7:30 PM -9:00 PM

Two-view geometry Computer Vision Spring 2018, Lecture 10

Geodetic Deformation Analysis as the main part of monitoring projects. Using Adjustments in a Monitoring Network. Dr. I Milev

Geometric Transformations

Lecture 5 2D Transformation

Geometric Approach to Position Determination in Space: Advantages and Limitations

ME 597: AUTONOMOUS MOBILE ROBOTICS SECTION 2 COORDINATE TRANSFORMS. Prof. Steven Waslander

LOCAL GEODETIC HORIZON COORDINATES

MATRIX REVIEW PROBLEMS: Our matrix test will be on Friday May 23rd. Here are some problems to help you review.

LIGHTWEIGHT DESIGN OF SEAT CUSHION EXTENSION MODULES USING THE PROPERTIES OF PLASTIC AND HCA-SIMP

Image Processing. Image Features

Chapter 14: Controlled Remote Tug Telemetry

Advance Convergence Characteristic Based on Recycling Buffer Structure in Adaptive Transversal Filter

CHAPTER 2 SENSOR DATA SIMULATION: A KINEMATIC APPROACH

High-Precision Positioning Unit 2.2 Student Exercise: Calculating Topographic Change

Computer Project #2 (Matrix Operations)

Rigid ICP registration with Kinect

Reconstruction of the Vibro-Acoustic Field on the Surface of the Refrigerator Compressor by Using the BEM-Based Acoustic Holography

Instantaneous Multi-Baseline Ambiguity Resolution with Constraints

Radial Basis Function Networks: Algorithms

Partial Sound Field Decomposition in Multi-Reference Nearfield Acoustical Holography by Using Optimally-Located Virtual References

Shape spaces. Shape usually defined explicitly to be residual: independent of size.

3D Reconstruction with two Calibrated Cameras

Chapter 2 File Management (updated September 5, 2009)

5 Classifications of Accuracy and Standards

CALCULATING TRANSFORMATIONS OF KINEMATIC CHAINS USING HOMOGENEOUS COORDINATES

Matching. Compare region of image to region of image. Today, simplest kind of matching. Intensities similar.

10. Cartesian Trajectory Planning for Robot Manipulators

Evaluating the Performance of a Vehicle Pose Measurement System

Strapdown system technology

The Determination of Telescope and Antenna Invariant Point (IVP)

TechnicalNotes. Trimble Total Control Software

Lightweight Corrugated Plate for Cargo box by using DOE and Topology Optimization

Wireless Sensor Networks Localization Methods: Multidimensional Scaling vs. Semidefinite Programming Approach

Graphics and Interaction Transformation geometry and homogeneous coordinates

MLEP Introduction to GPS Exercise Four Exercise requires pencil/pen to answer questions

COMP30019 Graphics and Interaction Transformation geometry and homogeneous coordinates

Basic Triangle Congruence Lesson Plan

CHARACTERIZATION AND CALIBRATION OF MEMS INERTIAL MEASUREMENT UNITS

An idea which can be used once is a trick. If it can be used more than once it becomes a method

Chapter 18. Geometric Operations

A Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation

A comparative study of the integer ambiguity validation procedures

Kinematics Analysis of Free-Floating Redundant Space Manipulator based on Momentum Conservation. Germany, ,

Principal Coordinate Clustering

Human pose estimation using Active Shape Models

Towards Real-time Geometric Orbit Determination of Low Earth Orbiters

Stereo and Epipolar geometry

Surveying. Session GPS Surveying 1. GPS Surveying. Carrier-Phase (RTK) Pseudo-Range (DGPS) Slide 1

BOWLING GREEN VERTICAL CONTROL NETWORK (BGVCN)

DIMENSIONAL SYNTHESIS OF SPATIAL RR ROBOTS

Relating Local Vision Measurements to Global Navigation Satellite Systems Using Waypoint Based Maps

Math background. 2D Geometric Transformations. Implicit representations. Explicit representations. Read: CS 4620 Lecture 6

Short on camera geometry and camera calibration

Matrices. Chapter Matrix A Mathematical Definition Matrix Dimensions and Notation

Convert Local Coordinate Systems to Standard Coordinate Systems

Mathematics Standards for High School Geometry

ifp Universität Stuttgart Performance of IGI AEROcontrol-IId GPS/Inertial System Final Report

GEOMETRY. Changes to the original 2010 COS is in red. If it is red and crossed out, it has been moved to another course.

RELIABILITY OF PARAMETRIC ERROR ON CALIBRATION OF CMM

Determination of the Parameter for Transformation of Local Geodetic System to the World Geodetic System using GNSS

General Instructions. Questions

2D/3D Geometric Transformations and Scene Graphs

Fig 1. Geometry of DGPS

A GENTLE INTRODUCTION TO THE BASIC CONCEPTS OF SHAPE SPACE AND SHAPE STATISTICS

Transcription:

Study on Vertical Alignment Maintenance echnique using GNSS in Skyscraper Eunchurn Park 1, Yu-Seung Kim 1, Joong-Yub Lee 1, Jun-Sung Choi 1* and Yeon-Back Jung 2, Won-Kyun Seok 2, Kwang-Soo Jung 2, Soon-Jeon Park 2, Joo-Ho Lee 2 1 Korea Maintenance Co., LD., Seoul, South Korea (el:82-2-830-7071, E-mail: eunchurn@kmbest.co.kr, yuseungkim@kmctech.co.kr, joongyup@kmctech.co.kr, eunchurn@kmbest.co.kr, kimchan@kmbest.co.kr, ceo@kmctech.co.kr * ) 2 Lotte Engineering & Construction, Seoul, South Korea (el:82-2-718-4688, E-mail: jung61@lottenc.com, archief@lottenc.com, ksjung716@lottenc.com, soon1026@lottenc.com, joo7777@lottenc.com) Abstract: In this study, vertical alignment maintenance technique of building form was implemented by using GNSS(global navigation satellite system) in skyscraper. An example test building was chosen as apartment building named as Lotte Castle Firenze which is under construction in Pusan, Korea. For the faster construction term of works, RK-DGPS was chosen and for insuring the data set accuracy of RK- DGPS measurements, the network adjustment method was used. his method is data processing which has compensating and differential process by using multi-reference point data sets except survey point. he standard deviation of survey point can be reduced by optimizing the reference network. he optimization method was chosen as the Procrustes analysis which involves matching shape configuration, through translation, rotations, and possibly scaling, to minimize the Euclidean distance between them. his technique effectively assumes that the coordinate have an isotropic covariance structure. he four point reference sites were chosen for Procrustes analysis. he network adjustment result showed that the standard deviation of survey point could be reasonably minimized. Keywords: skyscraper, real-time kinematic, DGPS, network adjustment, Procrustes analysis, form adjustment 1. Introduction In recent years there has been considerable interest in the construction of super high-rise buildings. o keep pace with the recent high-rise construction techniques, GNSS technology has been developed tremendously with effort in architectural and civil engineering measurement and survey field each year. From the prior art, various procedures and devices for surveys during and after the phase of erection of a high-rise building are known. High-rise buildings are subject to strong external tilt effects caused, for instance, by wind pressures, unilateral thermal effects by exposure to sunlight, and unilateral loads. Such effects are a particular challenge in the phase of construction of a high-rise building, inasmuch as the high-rise building under construction is also subject to tilt effects, and will at least temporarily lose its as a rule exactly vertical alignment. Yet construction should progress in such a way that the building is aligned as planned, and particularly so in the vertical, when returning into an untilted basic state. his paper has been implicated in a GNSS system into our skyscraper project that maintains a construction by a mechanized system in each development. With this research, a system that mechanically control and maintain the skyscraper s verticalalignment was developed and arranged. he GNSS system on the field was applied by creating an accurate vertical-alignment data and send to the central station. With this data set, the central station creates the system of structure and attitude and applies to the field condition to confirm the system s possibilities. here are several structures in architecture and civil engineering cases that apply a GNSS system, but it has limitation of use in measurement. But with advance development of GNSS technology and its accuracy, a measurement system has gotten close to more accurate and precise operating system. In addition, these types of measurement research are evaluated with the accelerometer in safety of structure and usability during wind load. But if we implicate a GNSS system in the on-going construction field, it would be difficult to measure the displacement at the present due to vibrations, a short measuring time and the other noise problems. o secure our accuracy of the system, three RK reading stations were used with real-time network adjustment reference coordinate, which lead to adjust the structure coordinate. 2. Elements of the theory In this chapter, the theoretical frameworks describing the processes were summarized. First, the post processing DGPS method which has taken a long time to calculate and convince coordinates. It is not proper to apply to construction processes of building, because reducing times of the each term of works is very important element in the construction processes of a highrise building. In order to overcome the accuracy and precision of RK-DGPS, the real-time network adjustment using multipoint controlling and displacement compensation using accelerometers were applied to proposed method. 2.1 Real-time Procrustes analysis network adjustment 2.1.1 Basic concept of network adjustment

Reference 2 F = 2A A 2A B + ( L + L ) = 0 (5) where ( A A ) and ( + ) L L are symmetric matrices. Multiplying Eq. (5) on the left side by, Survey point 1 Reference 3 Survey point 2 Reference 1 Figure 1. Basic conceptual diagram of network adjustment For insuring the data set accuracy of RK-DGPS measurements, the network adjustment method was used. his method is data processing which has compensating and differential process by using multi-reference point data sets except survey point. As the coordinates of reference points were measured by the postprocessing DGPS method, it could be assumed as convincible reference coordinates. he standard deviation of survey point can be reduced by optimizing the reference network. he optimization method was chosen as the Procrustes analysis method. In order to transform the WGS84 coordination to M coordination, the translation matrix which was Geodetic to ECEF and ECEF to ENU transformation was needed. 2.1.2 Applied Procrustes analysis We suppose that A is the set of floating points surveyed from RK-fixed mode data. B is the set of reference points surveyed from the post-processed DGPS data. A = {( x, y ), i = 1,, k}, B = {( x, y ), i = 1,, k} (1) i i i i Orthogonal Procrustes problem [Schoenemann, 1970] is the least squares solution of the problem that is the transformation of a given matrix A into a given matrix B by an orthogonal transformation matrix in such a way to minimize the sum of squares of the residual matrix E=A-B. Matrices A and B are (p k) dimensional, in which contain p corresponding points in the k-dimensional space. A least squares solution must satisfy the following condition, minimize to : tr { E E} = tr {( A B) ( A B )} he problem also has another condition, which is the orthogonal transformation matrix, = = I Both of the conditions can be combined in a Lagrangian function, F = tr{ } + tr{ ( )} {( ) ( )} { ( )} tr{ A A A B B A B B} tr { L( I) } E E L I (2) F = tr A B A B + tr L I (3) F = + (4) + where L is a matrix of Lagrangian multipliers, and tr{ } stands for trace of the matrix. he derivation of this function with respect to unknown matrix must be set to zero. ( L + L ) 2 L + L AA AB + = 0 (6) 2 ( ) ( ) ( ) = A B A B A A (7) ( L + L ) = 2 Since ( A B) is symmetric, ( ) symmetric. Remind that ( + ) A B must also be L L is also symmetric. herefore, the following condition must be satisfied. ( ) = ( ) A B A B (8) Multiplying Eq.(8) on the left side by, ( ) = ( ) A B A B (9) and on the right side by, ( ) = ( ) A B A B (10) Finally, we have the following equation using Eq.(9) and (10), ( )( ) = ( ) ( ) A B A B A B A B (11) A B A B are symmetric. Both of them have same eigenvalues. Matrices ( A B )( A B ) and ( ) ( ) svd {( )( ) } = svd ( ) ( ) { } A B A B A B A B (12) where svd { } stands for Singular Value Decomposition, namely Eckart-Young Decomposition. he result is, his means that, s = s VD V WD W (13) V = W (14) Finally, we can solve the unknown orthogonal transformation matrix. Obviously, a direct computation of the sum of squared distances between corresponding points of two configurations seems the simplest way of measuring the degree of resemblance. However, such a direct computation is not very meaningful due to the likely arbitrary location, orientation, and scale of one configuration relative to the other. As noted in the above section, it is more appropriate to compare the difference between

configurations after linear transformations have been performed on one configuration relative to the other [Gower & Dijksterhuis, 2004]. Geometrically, a linear transformation can be seen as a rigid movement of a network, where the network equals the internal relationships of reference points in a configuration. Initially, two types of movement that preserves distances among the points are performed, translation and rotation. When comparing shapes of configurations, the origins to which they refer are irrelevant [Gower & Dijksterhuis, 2004]. Hence, a shape-preserving translation is performed that centers the coordinates of one object to the other. Finally, the coordinate of survey points could be attained as timehistory data set which has minimized standard deviation. 3. est model and configuration of equipments 3.1 est site est site is new construction site named Lotte Castle Firenze in Pusan as shown in Figure 2. Figure 4. Measuring Control Point 2 able 1 and 2 shows the comparison of CP1 and CP2 measurement result. As observed in these tables, the each baseline of CP1 and CP2 have 6.757mm difference with both method. able 1. Conventional measurement result of CP1 and CP2 Conventional Survey Measurement E N B.M(Level) Point 206122.630 187997.490 EL=16.825 CP1 206196.050 188001.330 EL=19.760 CP2 73.52035092 (m) CP1-CP2 Baseline able 2. RK measurement result of CP1 and CP2 GPS Measurement (RK-fixed) Latitude (dms) Longitude (dms) 35 11 28.2444984 129 04 4.2990396 CP1 35 11 28.36783412046 129 04 7.200966627 CP2 73.51359392 (m) CP1-CP2 Baseline 3.2 Installation, Equipment and measurement software Figure 2. Gang form plan and form survey point Figure 3 shows the aerial photograph of the reference sites for network adjustment. As shown in figure, the survey point should be in the triangle of three reference site. he GPS antenna was installed in the survey point, gang form of under construction building. Figure 5 shows the GPS antenna of the gang form and the equipment list is able 3. LOERIA BUILDING Construction site Office 4F 1F Figure 5. GPS antenna of the survey point Figure 3. Aerial photograph of the reference point s network able 3. Equipment list of main measurement system Equipment Name Model Manufacture Items GPS AeroAntenna Septentrio PolaNt Antenna ech., Inc. GPS Receiver RF Modem Septentrio NV FreeWave PolarRx2@ PolarRx2 PolarRx2e HP-900RE IM-800X009 Septentrio FreeWave ech., Inc 1EA 1EA 2EA

19.63 19.62 19.61 E component (m) Geodetic and ENU coordinates transformation, NMEA ASCII data conversion and decoder, Real-time Procrustes Network Adjustment module and Optimization programming modules were in the integrated software as shown in Figure 6. 19.6 19.59 19.58 RK-fixed raw data :σ =7.2846(mm), m=19.5864592(m) 3P Network Adjustment : σ =5.9235(mm), m=19.5898911(m) 4P Network Adjustment : σ =3.6905(mm), m=19.5901835(m) 19.57 19.56 19.55 (a) East Component 64.35 N component (m) 64.3 64.25 64.2 64.15 RK GANG FORM : σ =20.6448(mm), m=64.1953661(m) 3P Network Adjustment : σ =15.4139(mm), m=64.1959707(m) 4P Network Adjustment : σ =9.6031(mm), m=64.1911750(m) 64.1 64.05 (b) North Component Figure 8. Network adjustment result of survey point able 4 is the result of standard deviation of one hour measurement data. As observed in the table, 4P PANA(Procrustes Analysis Network Adjustment) could be the more minimized. (a) Front panel (b) Block diagram Figure 6. Integrated measurement software 4. est Result 4.1 Network Adjustment Result -55.05 64.85 RK fixed raw data: σ=6.3442(mm), m=64.7939(m) 64.84 3P Network Adjustment : σ=2.9456(mm), m=64.7977(m) 4P Network Adjustment : σ=3.4005(mm), m=64.7988(m) 64.83 RK fixed raw data: σ=17.2193(mm), m=-55.1939(m) 3P Network Adjustment : σ=7.8272(mm), m=-55.1956(m) -55.1 ests of network adjustment which is 3-point (only reference points) and 4-point (including a survey point as reference) optimization were implemented. As observed in Figure 7, the each time history for 1hour result of reference points have the same differential components and the standard deviation of each data were reduced. 64.82 64.81 64.8 64.79 4P Network Adjustment : σ=7.9889(mm), m=-55.1983(m) -55.15 able 4. Standard deviation of one hour measurement data Standard deviation (σ) Site Method East (mm) North (mm) RK-fixed 6.3442 17.2193 1F 2.9456 7.8272 3.4005 7.9889 RK-fixed 9.0285 18.3166 4F 1.8578 3.5686 2.0758 3.5534 RK-fixed 6.1846 21.1810 Lotteria 2.8771 4.3686 3.4724 8.5527 RK-fixed 7.2846 20.6448 Gang 5.9235 15.4139 Form 3.6905 9.6031 4.2 Form Adjustment Result Figure 9 shows the result of form adjustment result corresponding to each 9, 13, 16 and 19 story. As shown in Figure, those optimized points were in the control line each direction with 10mm. Between 16 story and 19 story, gang form was adjusted. -55.2-55.25 20 64.78-55.3 (a) 1F East (b) 1F North 15 RK fixed raw data: σ=9.0285(mm), m=258.8741(m) 258.91 3P Network Adjustment : σ=1.8578(mm), m=258.8813(m) 4P Network Adjustment : σ=2.0758(mm), m=258.8801(m) 258.9-83.95 4F 258.92 9 Story 13 Story 16 Story 19 Story 258.89 258.88 258.87 RK fixed raw data: σ=18.3166(mm), m=-84.0455(m) 3P Network Adjustment : σ=3.5686(mm), m=-84.0485(m) 4P Network Adjustment : σ=3.5534(mm), m=-84.0485(m) -84 10-84.05-84.1 258.86-84.15 258.85 5 (c) 4F East RK fixed raw data: σ=6.1846(mm), m=-11.4954(m) 3P Network Adjustment : σ=2.8771(mm), m=-11.4925(m) RK fixed raw data: σ=21.181(mm), m=75.0402(m) 3P Network Adjustment : σ=4.3686(mm), m=75.0411(m) 75.1 4P Network Adjustment : σ=3.4724(mm), m=-11.4919(m) -11.46 (d) 4F North -11.48-11.5 Y (mm) 4P Network Adjustment : σ=8.5527(mm), m=75.0358(m) 9 Story 19 Story 0 13 Story 75.05 16 Story -5 75 74.95-11.52 74.9 (e) LOERIA East -10 (f) LOERIA North Figure 7. Network adjustment result of each reference points As observed in Figure 8, the each time history result of survey point has the same differential components and the standard deviation of data was reduced in the control value range, 10mm. -15-20 -20-15 -10-5 0 X (mm) 5 10 Figure 9. Form adjustment result 15 20

5. Concluding remarks he vertical alignment maintenance technique of building form was implemented by using GNSS in skyscraper. An example test building, RK-DGPS, Real-time network adjustment and the four point reference sites were chosen. he network adjustment result showed that the standard deviation of survey point could be reasonably minimized. For the further study, with this system, sets up the accuracy standard that compare the data set from the accelerometer by creating an algorithm that has a common system mode. the data set from the accelerometer and the GNSS was combined, and then a system transfer function model which defines the condition of status and its standard was created. With this result, next test could be planned with more upgrade one than previous experiments. Acknowledgement he work presented in this paper was supported by the Lotte Engineering & Construction. Reference 1. Gower, J.C., (1975) Generalized Procrustes analysis. Psychometrika, 40(1), pp. 33-51. 2. Luo, B., Hancock, W.R. (1999) Feature matching with Procrustes alignment and graph editing. 7th International Conference on Image Processing and its Applications. 3. Schoenemann, P.H., Carroll, R., (1970). Fitting one matrix to another under choice of a central dilation and a rigid motion. Psychometrika, 35(2), pp. 245-255. 4. Gower, John C. and Dijksterhuis, Garmt B. (2004) Procrustes Problems, Oxford University Press.