Paper 2005
Paper 2005
net/publication/280147380
CITATIONS READS
16 791
3 authors, including:
Ashraf A. A. Beshr
Mansoura University
68 PUBLICATIONS 485 CITATIONS
SEE PROFILE
All content following this page was uploaded by Ashraf A. A. Beshr on 31 August 2017.
1 2
Dr. Zaki M. Ziedan Dr. Hisham M. Abou Halima
3
Eng. Ashraf A. A. Beshr
1
Associate prof., Public Works Dept., Mansoura University, EGYPT
2
Associate prof., Public Works Dept., Mansoura University, EGYPT
3
Assistant lecturer, Public Works Dept., Mansoura University, EGYPT
Abstract
The safety concepts form the basis of modern structures design and assessment codes. The detailed
information about the structural deformations can help to determine the health of these structures, as
well as to evaluate whether such deformations are the same as those the structure was designed to
tolerate. This paper investigates an integrated monitoring system for the estimation of the
deformation behavior of structural members. Three different surveying techniques (precise leveling,
one and two total stations measurement techniques) are presented to evaluate the deformation
behavior of structural members. The comparison study between the surveying and structural
techniques for computation the structural deformation is introduced and discussed. These
techniques are applied to measure the structural deformation of reinforced concrete beams. Finally,
the three used surveying techniques with a special mathematical and adjustment models can be used
in monitoring the structural deformations with high accuracy. Moreover, the surveying techniques
can be motorized to give the continuous monitoring the structural deformations. The results of the
practical measurements, calculations and analysis of the interesting deformation using least squares
theory and computer programs are presented.
Keywords: Deformation, Monitoring, Precise leveling, Total Station, High Strength Concrete,
Beams.
1. Introduction
The security of civil engineering works demands a periodical monitoring of the structures. In many
civil structures like bridges, tunnels and dams, the deformations are the most relevant parameters to
be monitored. So monitoring the structural deformation and dynamic response to the large variety
of external loadings has a great importance for maintaining structures safety and economical design
of man-made structures. Dial gauge, accelerometer, Tiltmeter, etc. are traditional tools and methods
to measure structure displacement, rotation and together with temperature, wind speed and direction
allow the comprehensive investigation of structure dynamics behaviors [4,6]. These tools must be
installed, maintained, and frequently recalibrated to produce reliable results. The collected data
from these tools need to be interpreted to obtain direct geometric results which in many cases is
very complicated procedure and out of the control of the general structural engineers [5, 11]. Hence,
a flexible surveying technique is needed to overcome these obstacles, and make the process of
measurements easier and more accurate.
When the double observing procedure with double scales (left and right scale) rods is used, four
readings are taken at each setup, and the height difference at each setup is calculated by:
Where BLi, FLi and BRi, FRi are the backsight and foresight readings respectively. The distances
between the instrument and expected rod position are not equal, so the systematic errors must be
taken into consideration by applying the necessary corrections to the raw data [10]. The variance of
∆ hi is then determined by the error propagation method as following:
2 1 2 1 2 1 2 1 2
σ ∆ hi = σ B Li + σ F Li + σ B Ri + σ F Ri ………….(2)
4 4 4 4
2.2 One total station technique:
L
M o n ito re d S tru c tu re
Z
H
B
N S
Y
α γ Z B
A XB
Z A
XA
YB
YA
X
P o in t ( A ) is th e k n o w n c o o r d in a te s p o in t ( o c c u p ie d S ta tio n )
P o in t ( B ) is th e m o n ito r e d p o in t ( S h e e t P ris m )
From figure (2), a local three-dimensional coordinates system is needed to calculate the spatial
coordinates of any target point. The X-axis is chosen arbitrary as a horizontal line in the direction of
the base of the monitoring building, where the Y-axis is a horizontal line perpendicular to the
building base direction and positive in the direction towards the monitoring object, and the Z- axis
is a vertical line determined by the vertical axis of the instrument at occupied station. There is a
known coordinates point (A), and these coordinates are (XA, YA, ZA). From this point, we can
monitor the movements of any point (B) in space in order to determine its local coordinates (XB,
YB, XB).
X =X + ∆X
B A A−B
X =X + S cosγ .sinα
B A
Y =Y + ∆Y
B A A−B ………….. (3)
Y =Y + S.cosγ.cosα
B A
Z =Z + ∆Z
B A A−B
Z =Z + S.sinγ
B A
Where:
(S) is the slope distance between the instrument and the monitoring point, (α) is the horizontal
angle, and (γ) is the vertical angle. From figure (2), we have three unknown parameters (XB, YB,
XB) i.e. U=3, and there is no redundancy (r = 0). Then this has a unique solution, so the multivariate
propagation technique will be used [9]. By using the covariance law:
C X = J .C l .J T
……………(4)
( 3 ,3 ) ( 3 ,3 ) (3,3) ( 3,3 )
Where:
(CX) is the variance covariance matrix of unknowns, (J) is the Jacobian matrix, and (CL) is the
variance covariance matrix of observations. We can put equations (3) in the matrix form as follows:
∂X ∂X ∂X ∂X ∂Y ∂Z
2
∂S ∂α ∂ γ σ l 0 . 0 0 . 0 ∂ S ∂S ∂S
σ X 2 σ XY σ XZ
2
∂Y ∂Y ∂Y ∂X ∂Y ∂Z
σ YX σ Y σ YZ = . 0 . 0 σ α 2 0 . 0 .
∂S ∂α ∂γ ∂α ∂α ∂α
2
σ ZX σ ZY σ Z 2
∂Z ∂Z ∂ Z 0 .0 0 .0 σ γ ∂ X ∂Y ∂Z
∂S
∂α ∂ γ ∂γ ∂γ ∂ γ
Then, we can deduce that:
∂X 2 2 ∂X 2 2 ∂X 2 2
σ 2X = ( ) σS +( ) σα + ( ) σγ
∂S ∂α ∂γ
……(5)
σ 2 X = (cos γ sin α ) 2 σ S 2 + ( S cos γ cos α ) 2 σ α 2 + ( − S sin γ sin α ) 2 σ γ 2
∂Y 2 2 ∂Y 2 2 ∂Y 2 2
σ 2Y = ( ) σ s + ( ) σα + ( ) σγ
∂S ∂α ∂γ
σ 2Y = (cos α cos γ )2σ S 2 + (−S cos γ sin α )2σ α 2 + (−S sin γ cos α )2σ γ 2 …. (6)
2 ∂Z 2 2 ∂Z 2 2 ∂Z 2 2
σ Z = ( ) σ + ( ) σ α + ( ) σ γ
∂α ∂γ
s
∂S
2
……(7)
σ 2
Z = (sin γ ) 2 σ 2
s + ( S cos γ ) 2 σ γ
L
M onitored Structure
Z
H
B
γ2
γ1
α2
α1
From figure (3), there are three unknowns (XB, YB, ZB) and six observations (two slope distances
S1, S2, two horizontal angles α 1, α 2, and two vertical angles γ1, γ2). Then the least squares
adjustment technique will be used to get the most probable value of unknowns. The observation
equation technique will be used. In this model of adjustment, the number of equations is equal to
the number of observations (n = 6).
The two lengths of the lines (S1, S2) in the space can be written as:
S1 = (X B − X A ) 2 + (Y B − Y A ) 2 + ( Z B − Z A )2
…………(8)
S2 = (X B − X C ) 2 + (Y B − Y C ) 2 + ( Z B − ZC )2
The two lines (S1, S2) can be resolved into two components in the horizontal projection as:
AB = S1 cos γ1
………….(9)
CB = S2 cos γ2
By using the coordinate’s formulae, the two lines AB and CD can be written as follows:
AB = ( X B − X A ) 2 + (Y B − Y A ) 2
CB = (X B − X C ) 2 + (Y B − Y C ) 2 ………(10)
From figure (3), the two horizontal angles (α 1 and α2) can be calculated as follows:
AB 2 + AC 2 − CB 2
α 1 = arccos ( )
2. AB. AC
………(11)
AC 2 + BC 2 − AB 2
α 2 = arccos ( )
2. AC .BC
By using the coordinates formulae, we can write equation (11) as:
From figure (3), the two vertical angles (γ1 and γ2) can be calculated as following:
ZB − Z A
γ 1 = arctan [ ]
2 2
( X B − X A ) + (YB − YA )
…………(13)
ZB − ZC
γ 2 = arctan1[ ]
( X B − X C )2 + (YB − YC )2
The equations (8), (12) and (13) are the six observation equations; these equations are nonlinear
functions of parameters. They can be treated by least squares adjustment technique. The first step in
the solution is finding the approximate values of unknowns. Referring to figure (3), the coordinates
of point (B) can be computed according to the following formula [1,7]:
0 AC . cot α 1
X =
cot α 1 + cot α 2
B
AC …………..(14)
Y B0 =
cot α 1 + cot α 2
Z B0 = Z A + S 1 . sin γ 1
By substitute these approximate values in the six observation equations, the approximate values of
observations (L0) can be computed, then the steps of least squares technique using the observations
equation can be done to calculate the coordinates of point (B) and their accuracy.
Each movement vector has magnitude and direction expressed as point displacement coordinate
differences. These vectors describe the displacement field over a given time interval. Comparison of
the magnitude of the calculated displacement and its associated accuracy indicates whether the
reported movement is more likely due to observations error [6].
│dn│< (en)
Where:
│dn│ is the magnitude of the displacement for point n. It can be calculated as:
dn = ( ∆ X ) 2 + ( ∆ Y ) 2 + ( ∆Z ) 2 …………….(16)
and (en) is the maximum dimension of combined 95% confidence ellipse for point (n), it can be
calculated as following [6]:
e n = 1 . 96 . (σ ) 2 + (σ i ) 2 …………….(17)
f
Where:
σf is the standard error in position for the (final) or most recent survey,
σi is the standard error in position for the (initial) or reference survey.
Then
│dn│< (en) the point isn’t moved.
│dn│> (en) the point is moved.
4. Experimental program
The structural application consists of four reinforced concrete beams, which are tested in Steel
laboratory in faculty of Engineering - Mansoura University to estimate the deformation of these
beams subjected to specified loads [3]. The four tested beams have the same section (225 cm*20
cm*12 cm), but differ in reinforcement. They have the same number of bars at upper and lower
reinforcement. Two of them have 2Φ12 and the others have 2Φ16. The steel used is high mild steel.
The beams also have 5Ø6/m/ as stirrups. High Strength Concrete (HSC) mix is used [3].
Ordinary Portland cement and natural sand with high fineness modulus of 2.65 and Coarse
aggregate (natural gravel) with a maximum size of 12 mm are used. Powder silica fume with SiO2
of 92%, specific gravity of 2.2 and specific surface area of 16.8 m2/g is used. High Range Water
Reducers (HRWR) superplasticizers with trade name (Conplast 430) are used to improve both fresh
and hardened concrete properties. The High Reinforced Concrete is applied for this paper.
The automatic precise level NI007, two 3m double scale invar staves, two total stations (DTM 850
and SET300), sheet prisms of diameter 1 cm and calibrated dial and strain gauges are used in the
field measurements. The accuracy of all instruments and effect of the systematic errors are taken
into considerations during the practical measurements.
The searching for the “best” deformation model is based on either a priori knowledge of the
expected deformations or a qualitative analysis of the deformation trend deduced from all the
observations taken together. So, investigation and studies on the beam deformation by the available
structural software are done to predict the expected deflections values. Structural Analysis Program
(SAP) is a software program using the finite element system to find the straining actions at each
nodal of the element. The model geometry is entered in terms of features, which are subdivided,
into finite elements in order to perform the analysis but the analysis is done in the elastic stage.
7. Discussion on results
Actually, the full results are shown in M. Sc. thesis [3]. The important and vital results are
displayed in the fowolling:
7.1 Analysis of precise leveling observations
Precise leveling tested the first beam; L= 225 cm
the beam face was divided into four L/2 =112.5cm L/2 =112.5cm
Table (2) Comparison between deflections from precise leveling and dial
gauge readings
It is obvious that the resulted deflection values from the precise leveling analysis are very close to
those from dial gauge readings. The differences between the two techniques are too small. The
relationship between the acting forces and the resulting deflections deduced from precise leveling
for section (3) is illustrated in figure (5). As indicated in table (2), the deflection values at section
(2) are close to the values at section (4) because the two sections are at equal distance from support.
The maximum deflection values at section (3), because this section is at mid span of the beam. The
values of deflection at section (1) – at support- due to the rotation of the beam resulting from the
loading
4.50
3.75
3.00
Load
2.25
(ton)
1.50
0.75
0.00
0 5 10 15 20 25 30 35 40
Deflection (mm)
Figure (5) The relationship between the acting loads and the
resulting deflection from leveling at Sec. (3)
7.2 Analysis of one total station observations
L= 225 cm
L0/5 = 41 cm L0/5 = 41 cm
Load cell
2 cm 2 13
10 9 8 7 6
16 cm d = 20 cm 5 6\m
1 2 3 4 5
2 cm 2 13
Support Support b = 12 cm
L0 = 205 cm
(b) Cross Section
10 cm (a) Vertical plane with monitored points 10 cm
of beam
A comparison between the magnitude of the calculated coordinate differences and their associated
accuracy in Z direction for loads 0.35, 0.7 and 1.05 ton as a sample output is shown in table (3). The
values of ∆Z and enV can be calculated as:
2 2
e nv = 1.96 * σ ∆h i + σ ∆h0 …………(18)
and
∆zi = zi − z 0 ………….(19)
Table (3) Comparison of the magnitude of the calculated coordinate differences and their
associated accuracy for loads (0.35 ton, 0.7 ton, 1.05 ton)
The adjusted vertical displacements of ten monitoring points under all cases of loading are
calculated. Since the displacement of each point was known, the contour lines of the movements in
vertical direction (Z-direction) can be graphically illustrated by using Surfer program as shown in
figure (8).
The plan spacing between the contour lines indicates the steepness of slopes. From figure (8), it is
obvious that the maximum displacement occurs at mid span and the value decreases towards the
support. A Comparison between the deflection values from one total station technique and dial
gauge readings is illustrated in table (4).
The resulted deflection values from the one total station analysis are very close to those obtained
from dial gauge readings. The differences between the two techniques are too small.. By using the
same structural analysis technique, the adjusted displacements in X- direction can be calculated.
The displacements in X-direction for upper raw of monitoring points (6, 7, 8, 9 and 10) are greater
than the lower raw of monitoring points (1, 2, 3, 4 and 5). The maximum displacement value at
point 10 and load 4.2 ton, and this value is 6.93mm. The directions of point’s displacements in X-
direction at failure load can be graphically shown in figure (9).
Load cell
10 9 8 7 6
1 2 3 4 5
Z
Y
The contour lines of the movements in horizontal displacement (X-direction) at failure load can be
graphically illustrated by using Surfer Program as shown in figure (10).
The contour lines are closely spaced at support, which indicates a steep gradient. This means that
the displacements in X-direction at support
10 9 8 7 6
are strongly varied. But the contour lines / /
10 6
9/
8/
are widely spaced at mid span, which 7/
L= 225 cm
The last beam is tested by using the two
L/2 =112.5cm L/2 =112.5cm
total stations technique. The beam face is
L0/5 = 41 cm L0/5 = 41 cm
divided into five monitoring points, the Load cell
A comparison between the deflection values obtained from the two total stations technique and dial
gauge readings is illustrated in table (5).
Table (5) Comparison the vertical displacement between two total stations
analysis and dial gauge readings
It is obvious from table (5), that the deflection values from the two total stations technique are very
close to dial gauge readings from p= 0.35 ton to load p=3.85 ton. After load p=3.85 ton, there is
obvious difference because of the vibrations of dial gauge during loading especially the dial gauges
are placed under the tested beam. This dial gauge vibration is one of the disadvantages of using dial
gauge in structural members monitoring.
8. Conclusions
The results of experimental work lead to the fowolling conclusions:
1. The three used surveying techniques (precise leveling, one total station and two total stations) can
provide valuable data on the deflection of the structural members and movement of buildings
because the resulted deflection values from surveying techniques with the discussed adjustment
techniques are very close to the values from dial gauge readings.
2. Precise leveling is valid as a monitoring tool. It can read up to (0.05mm) with deflection accuracy
up to (±0.03 mm).
3. One total station has the potential to be used in monitoring the deformations of the structural
members because it can record the distance to 0.1mm with accuracy (±0.05 mm) in X-direction,
(±0.8 mm) in Y-direction and (±0.14 mm) in Z-direction.
4. The accuracy of the monitoring target coordinates is improved if the two total stations are set in
the site at their best locations instead of using one total station.
5. The accuracy of the monitoring points from one total station technique depends strongly on
horizontal angle, vertical angle and the instrument position distance. The best parameters are:
Horizontal angle (α) = 0 and vertical angle (γ) = 0
For distance, there is no optimum distance. But when the distance (S) between the instrument and
the monitoring point increases, the accuracy will decrease. The effect of horizontal angle variation
on the accuracy in X- direction is more than the effect of variation of the instrument position
distance. For Y-direction, the effect of the instrument position variation is more than horizontal and
vertical angle variation. In Z-direction, the accuracy depends mainly on the vertical angle value.
6. The upper surface displacement in X-direction of the beam is more than lower surface due to the
rotation of the beam under loading.
7. Achieving the required accuracy for surveying monitoring technique is based on the following
factors:
a- The used instruments specifications (instrument resolution, data collection options and the proper
operating instructions).
b- The field observing and modeling procedures.
ACKNOWLEDGMENT
The experimental work was carried out in the structure and steel laboratory of civil engineering
department of Mansoura University, EGYPT. All the academic and technical staff at the public
works and structural departments is gratefully acknowledged.
REFERENCES
1- Allan A.N., 1996, “ Surveying Building Surface by Theodolite Intersection” Survey Review.
2- A. M. Behairy, 1991, " Application of First Order Design Problem to Building Construction
Surveying Networks”, Faculty of engineering, Shobra, Zagazig University. Paper, CERM-Vol.
4- Brown C.J., Karuna R., Ashkenazi V., Roberts G.W. and Evans R.a, 1999, “ Monitoring of
structures using the Global Position System”, Institute of Engineering Surveying and Space
Geodesy, U.K.
5- G.W. Roberts, A.H. Dodson and V. Ashkenazi, 2000," Comparison of GPS Measurements and
Finite Element Modeling for Deformation Measurements of the Humber Bridge”, Institute of
Engineering Surveying and Space Geodesy, U.K.
6- Joseph Schroedel, June 2002," Engineering and Design- Structural Deformation surveying ",
U.S Army Corps of Engineering, U.S.A
7- Khalil M.M., 1991, “Best Location of Theodolite stations for Monitoring Deformations of
buildings”. M.Sc. Thesis, Cairo University, 84 pp
8- Mikhail, E.M. and Gracie, 1981, "Analysis and Adjustment of Survey Measurements", Van
Nostrand Reinhold Company, New York. U.S.