Solvin 1
Solvin 1
Solvin 1
er
;Li > ;, 4
i CF \ ( ; ie So eS we 7
if X\VVJ o
Las
a -
AE G39
d oe J ; <¢ 7
Pa ; c
a
<—. Seth. Ae ee
— = is
+ *
.
;
-
ia - - Sf Fe
°
'
@2oe ana aaaaa
TRE,
FE EF
Sea e2e20
0 =
"s
Bry
il
ll iilulla
DATE DUE FOR RETURN
UNIVERSITY LIBRARY
31 may 2011
UNIVERSITY LIBRARY |Con~
GGL 27
f ./ SEP 2005
une Glas?
UNIVERSITY NBRap
Te 2 DM
{SR
*
2007
Thic book ay—be-recalled before the above date.
OE
gat Se ee 4
| 90014
TEEEEEEEIIIEE EEE
a
3 LONGMAN
Pearson Education Limited
Edinburgh Gate, Harlow
Essex CM20 2JE, England
and Associated Companies throughout the world
Bannister, A. (Arthur)
Solving problems in surveying / A. Bannister, R. Baker. — 2nd ed.
Pacino:
Includes index.
ISBN 0-470-23424-5 (Wiley)
1. Surveying—Problems, exercises, etc. I. Baker, R. (Raymond),
1956- . II. Title.
TA537.B36 1994 yo
ISBN 0-582-23644-4 ,
Preface vii
Levelling 1
Distance measurement 32
Theodolite and traverse surveying 70
Triangulation and the National Grid 103
Areas and volumes 143
Curve ranging 190
Setting out and point location 242
WW Observations
WH
Oth
Gh
aon and adjustments 288
The computer programs 342
Index 344
Preface
This book, first published in 1989, was written to assist students who are
preparing for examinations in surveying. The same intention applies to this
revised edition in which some changes have been made to reflect advances
in modern instrumentation and techniques, and to take advantage of the
comments fed back by various lecturers using the book. The overall format
has been retained in that, although the book is devoted primarily to the solution
of problems typically encountered in such examinations, each chapter
commences with some basic theory. Frequently, the solutions themselves are
introduced with mention to the salient features involved, before being developed
in accordance with the theoretical principles, and are supported in some cases
by simple computer programs in BASIC.
An appreciable number of the worked examples and problems have been
taken from past examination papers, and their sources have been duly indicated.
It must be emphasized that, although permission has been given for the
reproduction of questions as set, the solutions are ours and responsibility for
their validity rests entirely on our shoulders. Thanks are due to the Senates
of the University of London, University of Bradford and University of Salford,
to Professor Cusens of the University of Leeds, the Council of Engineering
Institutions and the Engineering Council and Councils of the Institution of Civil
Engineers, Structural Engineers and the Royal Institution of Chartered
Surveyors for their kind permission to reprint questions set in their examination
papers.
A. Bannister
R. Baker November 1993
|
Levelling
Equipment The basic items of equipment are the optical level and a graduated staff,
although when long sight distances are involved a theodolite and target can
be used. This latter method is termed trigonometrical levelling.
Surveyor’s telescope A telescope is an essential feature of both the level and theodolite since its
overall magnification facilitates the pointing on to a distant target under
satisfactory atmospheric conditions. The optics include a diaphragm or reticule,
often referred to as the crosshairs, which gives an internally fixed aiming mark.
The line of sight of the telescope is given by the line containing the intersection
of the vertical and horizontal crosshairs and the optical centre of the object
lens of the telescope. The image of the target has to lie in the plane of the
crosshairs, and in modern instruments this is effected by means of a movable
lens (Fig. 1.1).
Figure 1.1
LEVELLING 1
Level line : A level in proper adjustment, and correctly set up, will produce a horizontal
sight line which is at right angles to the direction of gravity. The sight line
is tangential to the level line at instrument height; this line follows a constant
height above mean sea level and hence is a curved line, as shown in Fig. 1.2.
Over short distances, such as those met in civil engineering work, the two
lines can be taken to coincide.
Staff held
i vertically
Horizontal line
Level line
*
\ |
Figure 1.2
Long sight lines Over long distances, a correction is required to reduce staff readings given
by the horizontal sight line to the level line equivalent. Refraction of the sight
line has also to be taken into account. The technique of reciprocal levelling
can be adopted to eliminate the need for the corrections when long sights are
being used.
1.1 Booking and A properly adjusted tilting level was set up at a point P and the
reducing levels following consecutive readings were taken on a staff positioned at points
A, B and C, respectively: 0.663, —0.841, 0.939.
The level was then moved to point Q and further readings at C and
D, respectively, were taken as follows: 1.198, 1.100.
Use this example to explain what is meant by the terms backsight,
foresight, intersight and change point. Book, reduce and check the levels
using both standard methods, given that the reduced level of A was
94.115 m AOD. (Note that the reading at B was taken on an inverted
staff and has therefore been recorded with a negative sign.)
What are the advantages and disadvantages that you would associate
with each method of booking? [Salford]
Solution Figure 1.3 shows the readings listed in the question and illustrates
the fact that the levelling involves the interconnection of separate groups, i.e.
AC and CD.
Backsight The first reading when the level has been set up (at P) is known as a backsight
(BS). The staff is held vertically at a point (A) of known reduced level; this
could be a bench mark but it need not be. When the instrument position has
to be changed, the first sight taken in the next section is also a backsight.
B — — — —
0.841 1.100
0.663---— -—:-+-— -—: »----—-—: —}1.198 QO
f 0.939 D
A P C a
Inverted Change
staff point
Figure 1.3
Foresight The last staff reading from a level station on to a staff held vertically at a point
is termed a foresight (FS). It is thus the last reading taken within a section
of levels before moving the instrument, and also is the last sighting made over
the whole series of levels. It is good practice to finish at a point of known
reduced level or to return to the start point so that an appraisal of the overall
accuracy of the work can be assessed from any misclosure.
Intersight The term ‘intersight’ covers all sightings and consequent staff readings made
between the backsight and foresight within each section. In Fig. 1.3 the inverted
staff reading is an intersight.
Change point To allow the levelling to cover an appreciable distance or to cater for difficult
terrain, the level position may have to be changed frequently. In order to relate
the various sections it is essential that a foresight and backsight be made on
a staff held at a change point. Point C in Fig. 1.3 fulfils this purpose. Change
points need to be firm features so that there is no relative displacement of the
staff when its face is turned towards the new position of the instrument. When
selecting such points it is good practice to try to arrange the lengths of the
backsight and foresight to be equal in order to eliminate collimation error and
any error due to curvature and refraction. Although modern levels allow
estimation of the staff reading to 1 mm over distances upwards from 100 m,
it is commonplace to have maximum sight distances of the order of 50 m.
Booking and reducing The level book will have columns headed ‘Backsight’, ‘Intersight’ and
‘Foresight’, together with others covering “Reduced Level’, ‘Distance’ and
‘Remarks’. It will also have provision for one of the two possible systems by
which the levels of specific points are calculated, namely the rise and fall
method and the height of collimation method.
Rise and fall method Note that, except for the change point, each staff reading is written on a separate
line so that each staff position has its unique reduced level. This remains true
at the change point since the staff does not move and the backsight from Q
is thereby related to the reduced level derived from the foresight from P. Hence
LEVELLING 3
Table 1.1
1.780
1.602
0.178
the two readings are entered on the same line. Successive staff readings are
used in calculating rise or fall between the points. From the booking we have:
A to B 0.663._— (—0.841) = + 1.504, 1. a.Fise,
B to C —0.841 — 0.939 = —1.780, i.e. a fall,
C to D 1.198 — 1.100 = +0.098, i.e. a rise.
Notice that a decrease in staff reading implies a rise from the first point to
the second, i.e. C to D.
The arithmetic involved in the reduction of the levels can be checked using
the following rule:
(ZFS — LBS) = (Lfall — Lrise) = (final level — first level).
Height of collimation —_Inso far as the booking of backsights, etc., is concerned the two systems are
method identical but the method of the reduction of levels is different. First, the height
of the collimation of the instrument for each of its positions is determined,
i.e. for the setting at P:
Table 1.2
Comparison of Less arithmetic is required in the reduction of levels with the collimation method
methods than with the rise and fall method, in particular when large numbers of
intermediate sights are involved. Moreover, the rule quoted for the rise and
fall method gives an arithmetical check on all the levels reduced, whereas only
change points are checked under the collimation method. However, a rule can
be established to give an arithmetic check for all reduced levels obtained by
the collimation method. Consider the levelling between stations A and B shown
in Table 1.3.
Table 1.3
b+x Station A
f Station B
The sum of all reduced levels except the known level at station A
=(b+x-f)+L0+x —- i)
=(n+1)60+2%)
—li -—f
in which n is the number of intermediate sights between A and B, and (b+ x)
is the height of collimation of the level.
This analysis can be readily extended to cover a number of consecutive
sections to produce the general rule:
LEVELLING 5
LD All reduced levels except the first = £ (Each height of collimation
times the number of intermediate and foresights related thereto) — LIS
— LFS.
1.2 Apportioning Reduce the levels below by the collimation method, apply appropriate
closing errors and checks. I is an OBM of specified level 56.174 m AOD. Adjust the
calculating gradients calculated levels so that the results at I agree with the specified level.
If the distance from A to H is 220 m calculate the mean gradient between
those points. Outline the practical precautions which must be taken for
accurate levelling.
Table 1.4
Remarks
G
H
|
[Bradford]
E G
Figure 1.4
Table 1.5
LEVELLING 7
Table 1.6
Oo <<
—-IrOoO7nmMoooaworD
1.3 Two-peg test (a) Explain how you would check a level for adjustment using the
‘classical’ two-peg test.
(b) Four stations P, Q, R and S were set out in a straight line such
that PQ = QR = RS = 30 m. A tilting level was set up at P and readings
of 2.148 m and 1.836 m were observed on a staff held vertically at Q
and R, respectively. The level was then set up at S and readings of
2.013 mand 1.755 m were observed on the staff held vertically in turn
at stations Q and R. The bubble was adjusted to be at the centre of its
run for each reading.
Check the overall adjustment of the level.
Solution (a) Figure 1.5 shows the method of conducting the two-peg test for
checking the adjustment of a level. Two rigid points A and B are marked on
the ground (the pegs) and the instrument is set up exactly between them at
point C. Readings are taken on to the staff held at A and B, and the difference
between them gives the difference in level of the pegs. The equality in length
of the backsights and foresights ensures that any instrumental error is equal
on both readings. The instrument is then moved so that it is outside the line
of the pegs and close to one peg. Readings are again taken on to the staff held
at points A and B, any discrepancy between the level difference given by the
first readings compared to that from the second readings is due to the instrument
being out of adjustment.
Error equal on
[= +}. ---| fase
both readings
—..
Figure 1 13
Figure 1.6
(b) This example shows an alternative method of carrying out the test. First
determine the apparent differences in level between Q and R. From Fig. 1.6,
(a) 2.148 — 1.836 = 0.312 m
(b) 2.013 — 1.755 = 0.258 m.
Since the two differences do not agree the horizonal axis of the telescope, called
its line of collimation, is not parallel to the bubble tube axis. To determine
the inclination of the line of collimation it is assumed that the line of collimation
is directed upwards, by 6, from the horizontal so for a horizontal sight line
the corrected staff readings are given in Table 1.7.
The two corrected differences in level are therefore
Table 1.7
2.148-—300 1.836—606
2.013—606 1.755—300
LEVELLING 9
The horizontal crosshair must be positioned by means of the tilting screw
to give readings of 2.067 m and 1.782 m, respectively, on the staff when held
in turn at Q and R. Actuating the tilting screw will displace the bubble which
must now be returned to its central position by means of the relevant adjusting
screws. The test would be repeated to ensure that readings conform to say
0.002 m.
1.735 *688
she
4=
PG RES. # ley Corrected
staff reading
Figure 1.7
Difference
1.5 Corrections for (a) What corrections should be applied to staff readings when level-
long sights ling over long distances? Derive an expression giving the combined
correction and determine the distance for which its value is 5 mm.
(b) In levelling across a wide river reciprocal observations gave the
following results for staffs held vertically at points P and Q from level
stations A and B respectively. A and P were near to each other on one
bank, whilst B and Q were similarly situated on the other bank.
Solution (a) In Figure 1.8 the level has been set up at MSL and sighted
towards a staff held vertically at a distance IH away. IH is the line containing
the instrument axis, but the line of sight has been refracted to give J as the
staff reading. The two corrections which have to be combined to give the
required correction, JK, are the correction for curvature, HK, and the correction
for refraction of the line of sight, HJ.
LEVELLING 11
Figure 1.8 .
IH
HJ = —
2R,
2 2 2
Combined correction JK = gbtsoed «Es = ae Pa Rg
2R 2R, 2R R,
IH? IH? 1000
Now = (km) = ——— IH? (m)
2R 12 740 12 740
= 0.078 IH? (m) with IH expressed in
kilometres.
JK = 0.078 IH* x $
= 0.067 IH? (m)
= 0.067 d* (m)
where d IH (km).
12 SOLVING PROBLEMS IN SURVEYING
If JK =5mm = 0.005m
0.005 = 0.067d?
Therefore
d = 0.273 km = 273 m.
Refracted
line of sight Level line
Figure 1.9
(b) Figure 1.9 shows the level stations and staff stations on each bank. The
line of sight and the level line can be considered the same between A and P,
B and Q, respectively, since the points are close together. Assuming no change
in refraction conditions when observing from A and then from B with the same
level, the level lines will give staff readings y below the actual staff readings
at Q and P, respectively.
The apparent differences in level are (1.495—1.180) m and
(1.730 —1.405) m, respectively. The corrected differences (H) in level are:
[1.495 — (1.180 —y)] and [1.730 —(1.405+y)], respectively. These differences
must be equal, and so
H = 1.495—1.180+y = 1.730—1.405—y.
Therefore 2H = (1.495 —1.180)+(1.730— 1.405)
or
LEVELLING 13
1.6 Trigonometrical Define a.
levelling (i) The coefficient of refraction (K,) in terms of the angle of
refraction and the angle subtended at the centre of the spheroid by the
arc joining the two survey stations.
(ii) The coefficient of refraction (Kz) in terms of the local radius of
the earth and the radius of curvature of the ray path. Prove from first
principles that K, = 2K, and quote an average value for the coefficient
of refraction.
Simultaneous reciprocal vertical angles have been observed, and the
recorded mean values are +06° 30’ 59.7” and —06° 31’ 37.5”, the
horizontal distance between the two observation stations being
1389.396 m. Assuming the local mean radius of the earth to be
6383.393 km, compute the coefficient of refraction, the corrections to
be applied for earth curvature and refraction and the corrected mean
vertical angle. [CEI]
Horizontal
at B
Vertical __,
at A
Vertical
ia di a WEP Be eal =/C atB
Horizontal
at A
Figure 1.10
a
sin (a+c/2—r)
BC
= AC
sin [180 —(90 +c/2) —(a+c/2—r)]
sin (a+c/2—r)
AC
cos (a+c—r)
As will be seen from Fig. 1.10 the correction for curvature and refraction
at A and B is (c/2 —r) and this refers angles a and 8 to chords AC and BD,
respectively.
Thus elevation angle BAC = at+c/2—r
and depression angle DBA = B—c/2+r
but BAC = DBA since AC and BD are parallel.
a“
Therefore BAC = =ac B= Dita.
A
= AC tan sesalal
2
+
i
Put K, aes
c
R
and K, =
R,
in which r is the angle of refraction, R is the mean radius of the earth, R;
is the mean radius of the sight line, and c is the angle contained between the
two verticals.
LEVELLING 15
In Example 1.5 the combined correction JK was evaluated as
pO2R GeeoRR 2
d2
eae
aR Kk2)
in which IH = d. Now from Fig. 1.8 in which JK = JJ x JiK, in which JiK
is the angle between the level line and the refracted line of sight.
A d
Then JIK = atl Ks) ‘
2R
From above c/2—r = c/2—K,c
= c/2 (1—2K))
d
—— (1-—2K,) since c = d/R.
2R
Therefore K, = 2K,.
K, depends upon temperature gradients, and in certain cases the line of sight
can curve upwards rather than downwards; an average value is 0.07 overland.
It is the coefficient used in practice, i.e. the coefficient of refraction is defined
as the ratio between the angle of refraction and the angle at the centre of the
earth.
Now determine the mean vertical angle.
Also chord AC = 2(R+AH4) sin c/2 in triangle ABC (Fig. 1.10). But for
practical purposes we can say
R+
chord AG: = arc AC = d BODINE
c |_ 1389,396 1
n— = ——— x —____
2 2 6 383 393
ie. <= 39,4",
2
Variables
Al = Elevation angle at A H1 = Height of instrument at A
A2 = Depression angle at B H2 = Height of instrument at B
A3 = Angle subtended at the H3 = Height of signal at A
centre of the earth H4 = Height of signal at B
A4 = Mean corrected angle L1 = Horizontal distance AB
Cl = Correction to angle Al M_ = Input/output, minutes
C2 = Correction to angle A2 R = Radius of the earth
C3 = Correction for curvature S = Input/output, seconds
C4 = Correction for refraction V1 = Ground level at A
C5 = Coefficient of refraction V2 = Ground level at B
D = Input/output, degrees
10 REM TRIGONOMETRIC LEVELLING
20 INPUT"HEIGHT OF INSTRUMENT AT A (M) ";H1
30 INPUT"HEIGHT OF INSTRUMENT AT B (M) ";H2
40 INPUT"HEIGHT OF SIGNAL AT A (M) "5H3
50 INPUT"HEIGHT OF SIGNAL AT B (M) ";H4
60 INPUT"ELEVATION ANGLE RECORDED AT A (D,M,S) ";D,M,S
70 Al=((3600*D) +(60*M) +S) /206264.8
80 INPUT"DEPRESSION ANGLE RECORDED AT B (D,M,S) ";D,M,S
90 A2=((3600*D) +(60*M) +S) /206264.8
100 INPUT"LEVEL OF GROUND AT A ABOVE MSL. (M) ";V1
110 INPUT"MEAN RADIUS OF THE EARTH (KM) ";R
120 R=R*1000
130 INPUT"HORIZONTAL DISTANCE AB (M) ";L1
140 A3=L1/(R+V1)
150 C1l=(H4-H1)/L1
160 C2=(H3-H2)/L1
170 A1l=A1-Cl
180 A2=A2+C2
190 A4=(A1+A2)/2
LEVELLING 17
200 V2=V1+INT((L1*TAN(A4) )*1000)/1000
210 C3=INT(A3*2062648!/2)/10 si:
220 C4=INT((A1+A3/2-A4) *206264.8%*10) /10
230 C5=ABS (INT (C4*1000/(2*C3)) /1000)
240 A4=A4*206264.8
250 D=INT(A4/3600)
260 M=INT((A4-(D*3600) ) /60)
270 S=INT((A4-(D*3600) -(M*60) ) *10) /10
280 PRINT"LEVEL AT B ="3V2;"M"
290 PRINT"MEAN OBSERVED ANGLE =";D;M;S
300 PRINT"CORRECTION FOR CURVATURE =";C3;"SEC"
310 PRINT"CORRECTION FOR REFRACTION =";C4;"SEC"
320 PRINT"COEFFICIENT OF REFRACTION =";C5
330 END
1.7 Levelling over Two survey stations A and B have been established on opposite sides
long distances of a large lake. The level of the station marker at A above water level
in the lake is 49.652 m and at B is 176.268 m. The distance between
A and B measured at water level is 67 256.4 m.
A theodolite with its trunnion axis 1.475 m above the station at A is
used to observe a target mounted 6.237 m above the station at B.
Determine the minimum clearance of the line of sight above water level.
How far from A does this occur?
Calculate the anticipated vertical angle at A. Will this be elevation
or depression?
Take the radius of the earth to be 6.382 48 x 10°m and the
coefficient of atmospheric refraction to be 0.082. [Bradford]
Vertical
at B
Vertical _ Horizontal
\ at A at A /
i Tangent to /
\ ae / refracted
i n|o line of sight
ee
67 256.4
ers tae. ee”
clearance /
\
Figure 1.11
Solution
Height of instrument at A above water level
= 1.475 + 49.652 = 51.127 m.
Height of target at B above water level
= 6.237 + 176.268 = 182.505 m.
131.378
Sean 206,265
67 256.4
ll 402.9”
; ee
Over distance x from A curvature correction
2
x | 206 265
2 ~ 6 382 480
= 0.016 159 x x”.
Refraction correction (r)
= 0.082 X c,
Il 0.082 x 2 x 0.016 159 x x
= 0.002 650 x x”.
Therefore combined correction
= 0.01351 x x”.
When x 67 256.4 m
cil2—r = 908.6”.
a + (c/2—r) = 402.9
a + 908.6 = 402.9.
Therefore a = —505.7”
—08’ 25.7” (i.e. depression).
(—505.7+0.013 51 Xx)x
Elevation of line of sight x from A =
206 265
For a minimum (by differentiation)
505.7 = 0.02702 x x
Therefore x = 18 716 m.
LEVELLING 19
at B is a depression. Such an eventuality can occur when the height of B above
A is small relative to distance AB. The expression
+
Hg — Hy =d X tan ote
was derived on the assumption that a was elevation and 6 was depression,
so in this case it is essential that a negative sign is given to a to acknowledge
that both measured vertical angles are depressions.
The expression for height difference in this case is
Hy — Hy = d x tan (F%)
*
a + (c/2—r) = 402.9”
= B — (c/2-n)
= B — 908.6”.
Therefore @ = 1311.5”.
131.372 m.
The very small discrepancy between the given height difference and that
calculated is due to rounding off during the whole calculation.
1.8 Eye and object Two survey stations A and B are 5126.1 m apart. During the course
correction of taking reciprocal trigonometrical levels between these stations the
readings in Table 1.9 have been recorded.
Table 1.9
Figure 1.12
Solution Determine the ‘eye and object corrections’. It will be noted that the
height of the instrument at A is not the same as the height of the signal at
B. Thus the observed angle of elevation does not refer to the ground levels
at A and B. This difference in height causes the observed vertical angle to
be larger than that which would be noted if one could observe directly from
those points. A correction (e), termed the ‘eye and object’ correction, is applied
to reduce the observed value to the required value. In its simplest form it can
be written as (Zg—h,)/d radian. (See Fig. 1.12.)
Hence:
Height of target at B = 5.00 Height of target at A = 2.00
Height of instrument at A = 1.50 Height of instrument at B = 1.40
360 x 60 x 60
Therefore angle c subtended at earth’s centre by AB
= a =" 1665162
30.86
Refraction
€
Therefore, Sor —r yp Bege
LEVELLING 21
Angle of elevation (a,) from A (corrected for eye and object)
=a, F cde — Fr
Matis 3.2” + 1 112" 1” 1s" 4.
If 6 is the angle of depression at B, the relevant eye and object correction
must be added to it to give the equivalent ground to ground measurement since
(Z, — hg) is positive.
Since a, + cl2—r = B, — (cl2—n)
= 1° 13’ 14.4”.
Therefore 6, = 6+24.1 = 1° 13’ 14.4” + 1’ 11.2”
= 1 )4°,25.6"
and thus B= 1° 14’ 1.5”.
+
= 109.2 m.
Had the slant distance between A and B been measured (as in Problem 15)
it should be reduced to its equivalent chord length at the mean height of the
stations.
1.9 Parallel plate Precise levels are normally fitted with a parallel plate micrometer.
micrometer Derive an expression relating the displacement of the line of sight d,
the thickness of the glass plate ¢, the angle of rotation of the plate i and
the refractive index of the glass n.
Calculate the thickness of a parallel plate micrometer made of glass
having a refractive index of 1.7 if it is specified that a 30° rotation of
the plate causes a 5.00 mm displacement of the line of sight. State clearly
any assumptions made. Check the answer from first principles in order
to determine whether or not the assumptions are justified. [Bradford]
Introduction In Fig. 1.13 the parallel plate, fitted in front of the objective
lens, is shown rotated through angle i from its vertical position. When vertical,
the line of sight to the staff will not be displaced vertically, but when the parallel
plate is rotated by means of a micrometer ‘knob’ that line can be directed
upwards or downwards. By adjusting the line to the nearest graduation of the
levelling staff the micrometer effectively allows readings to 0.01 mm.
Solution The derivation of the formula is as follows.
22 SOLVING PROBLEMS IN SURVEYING
'
|
position
j Displaced
ia line of sight
ll
position
Figure 1.13
In triangle PQR
PQ = t/cos 0
RQ = d.
Therefore d = : sin (i — @)
cos 0
ae : 6
=f (sn — cosi where sin i = n sin 6
Cos
sin i
cos i n
=tsini| 1 —
sin 1 sin2i
1 —- ae
n
pean cos 1
ll feSitbob: |) me eee
Vn?
—sini
fue 1 —sin‘i
tsini{l — Ss YS:
n=sint
a = mynioats
n
LEVELLING 23
Using the rigorous formula
ane V1i- sin2i
t X
et Vn sini
— Ss
d=
V1—-0.5?
5.0 x05 (1-
V1.77-0.57 ]
Therefore ¢ = 21.41 mm.
Problems
1 Book and reduce the following levels and carry out the necessary
checks on the arithmetic. A pair of numbers indicate a change point,
the first number being the backsight.
Staff Reading (m)
1.263 OBM, reduced level 26.294 m.
3.279, 0.796 Change point.
0.376 —_— Road level under bridge.
1.627, 0.291 Change point.
—2.162 Soffit of bridge arch.
1.582, 3.526 Change point.
2.014 TBM, reduced level 27.42 m.
Could a lorry 4.1 m high pass under the bridge? [Salford]
Answer No, clearance under bridge = 3.874 m.
2 Levels are taken to determine the height of two pegs a and b, and
to determine the soffit level of an overbridge. Using the values of levels
indicated in Table 1.10, and given that the first backsight is taken on
a BM at a church (RL 60.270 m), and the final foresight is on a BM
at a school (RL 59.960 m), determine the closing error and the height
between the underside of the bridge and the ground immediately below
it. Use both the collimation level and the rise and fall methods and apply
the usual checks. [Salford]
Answer 4mm; 5.473 m.
Table 1.10
4 Complete the extract in Table 1.11 from a level book, applying the
usual arithmetic checks. [Leeds]
Table 1.11
2.316
—0.874 Inverted staff
225.05 OBM
Table 1.12
88.41
Last RL — First RL
= 8.80
LEVELLING 25
The proposed improvement involves regrading the road to a 1 in 20
gradient, rising from chainage 0 to 300, afd passing through the existing
surface at chainage 180 m. Determine the maximum depths of required
excavation and fill. [Salford]
Answer Cut 1.17 m; Fill 5.17 m.
6 Reduce the levels in Table 1.13 by the rise and fall method, applying
appropriate checks.
Calculate the headroom of the bridge at A + 80 m and the mean
gradient of the ground between the points A and A + 120 m.
[Bradford]
Answer 5.093 m, 1 in 35.2. .
Table 1.13
——— =
level
7 The values indicated in Table 1.14 refer to levels taken along the
centre line of a proposed road at 30 metre intervals.
Determine the reduced levels, and the depth of the earthworks at each
point, given that the reduced level at chainage 100 metres is 46.400 metres
and formation level is to-rise at a gradient of 1 in 50 from a reduced
level of 45.000 m at chainage 100 metres.
Answer 1.4, 2.16, 3.77, 1.72. -
Table 1.14
100 RL 46.400 m
130
160
Foresight at change point
Backsight at change point
190
Table 1.15
BM 112.309 AOD
Point X
30 m from X
60 m from X
Change point
90 m from X
120 m from X
Inverted staff on underside
of bridge 126.8 m from X
Change point
150 m from X
Point Y (170 m from X)
BM 107.895 AOD
9 A level is set up at point X and readings of 0.219 and 1.674 are taken
on to two bench marks A and B respectively. The height of A is 166.84 m
above datum and of B is 165.37 m above datum. If the distances XA
and XB are 87.6 m and 33.8 m respectively, calculate the collimation
error per 100 m. If a further reading of 2.121 is taken from X on to
a point Y, 71.6 m from X, calculate the height of Y. [Leeds]
Answer 27.9 mm; 164.934 m.
Table 1.16
0.654
Te [| oe
0.356
LEVELLING 27
1.826 m on X and 1.493 m on Y were recorded. Determine the true
level of the top station. al
Answer 89.886 m.
Determine the difference in level between the two stations and the
coefficient of atmospheric refraction. Assume that the radius of the earth
is 6.37 X 10°m. [ICE]
Answer 139.27 m; 0.071
Table 1.17
If the height of Aga is 521.46 m AOD and the radius of the earth is
6380 km, calculate the heights of Beetle and Citation. [Leeds]
Answer 376.24 m; 346.34 m.
Table 1.18
Instrument at A
Height of instrument above ground = 1.45 m
Height of target at B above ground = 3.10 m
Vertical angle + 00° 11’ 12”
LEVELLING 29
Instrument at B
Height of instrument above ground = 1.52 m
Height of target at A above ground = 3.86 m
Vertical angle — 00° 20’ 52”
If the mean radius of the earth is 6375 km, calculate the coefficient of
atmospheric refraction and the height of station B above datum.
[Leeds]
Answer 0.062; 401.67 m.
24 Given that the refractive index from air to glass is 1.6 determine
the angular rotation of a parallel plate, of thickness 10 mm, to give a
vertical displacement of the image of the staff through 0.1 mm.
Answer 01° 31’ 40”
LEVELLING 31
2
Distance measurement
By tape (or wire) For many years the very accurate measurement of distance depended upon
the careful use of steel tapes or wires. Nowadays, EDM is used almost
exclusively for accurate work but the steel tape still is of value for measuring
limited lengths and for setting out purposes.
Up to seven corrections may be applied to the measured length, to give the
true length, dependent upon the circumstances of the measurement.
Correction for A steel tape will normally be provided with standardizing data, for example
standard it may be designated as 30 m long under a tension of 5 kgf at a temperature
of 20 °C when laid on the flat. With use the tape may stretch and it is imperative
that the tape is regularly checked against a reference tape kept specifically for
this purpose.
Correction for tension _ If the tape is of correct length under a standard tension and it is used under
a different tension the correction which should be applied is
(P—P,)L
AE °
where P is the tension applied in the field, P, is the standard tension, A is
the cross-sectional area of the tape, E is Young’s modulus for the tape material,
and L is the observed length.
The sign of the correction takes that of quantity (P—P,). Obviously, when
P is made equal to P, the correction is zero and the arithmetic is simplified
but the proviso#mentioned in the correction for sag should be borne in mind.
Correction for sag For very accurate work the tape can be allowed to hang in catenary, free of
the ground, between suitable supports. In the case of a long tape intermediate
supports can be used to reduce the magnitude of the correction (see Fig. 2.1).
If the tape has been standardized on the flat the correction which should
be applied to reduce the curved length to the chord length is
ee
Figure 2.1
w-L3
24 p*’
where w is the weight of the tape per unit length (P should be larger than 20wZ).
If such a tape is used on a plane surface which can be considered flat then
no correction is applicable.
If the tape has been standardized in catenary the equivalent chord length
is given by the graduations, provided that the standardizing tension is applied.
If such a tape is used on the flat then the correction is added to the tape length.
Correction for slope In surveying it is essential that horizontal lengths are determined. Thus length
L measured on the slope must be reduced to its equivalent plan length L cos
6 (see Fig. 2.2). The correction to be applied is
L
h
L cos @
Figure 2.2
2 94 s
Now cos 6 = 1— — + — + ..., which is nearly equal to 1 — rat
2! 4! 2L
h einer bass h?
where 6 = ra, Thus a simplified form of the correction is — Sap pro-
Correction for If a tape is used at a field temperature different from the standardization
temperature temperature then the correction is
iveaeea)
where «a is the coefficient of thermal expansion of the tape material, T is the
DISTANCE MEASUREMENT 33
field temperature, and Ts is the standardization temperature. The sign of the
correction takes the sign of (T— Ts). :
Correction to mean In the case of long lines in triangulation surveys the relationship between the
sea level length measured on the ground and the equivalent length at mean sea level
has to be considered. This is also the case when any distance measurements
are being related to National Grid co-ordinates (see Ch. 4). If the measured
length is L,, and the height of the line above datum is H then, since both
lengths subtend the same angle at the earth’s centre (see Fig. 2.3),
Figure 2.3
pS Ur Oe 4
and
Lnsi = RO,
f= Digs = |
R H+R
and so
R
Lnsi = Lm War
The correction to be applied is
R H
| OY Skee 1)RO hE) Ee apn
mst — Lm 0 a NE es Race
pase
H
ay
Electromagnetic Two main types of EDM are encountered in land surveying, namely the
distance electronic or microwave systems and electro-optical instruments. The principle
measurement (EDM) of operation is that a transmitter at the master station sends a modulated
2D = nv + spi. d\ + a constant,
2a
where n is the number of complete wavelengths contained within the double
distance and will be unknown. The purpose of deploying different modulation
frequencies is to evaluate this number by comparing the phase differences of
the various outgoing and measuring signals. Note that
. ae >
nf
where Cg is the velocity of the electromagnetic wave in a vacuum, f is the
frequency, and n is the refractive index of the medium through which the wave
passes.
Nowadays, most local survey work and setting out for civil engineering works
will be carried out using infrared based EDM, which falls within the electro-
optical group. The infrared carrier wave is transmitted to a passive reflector,
usually a retrodioptive prism, from which it is returned to the master at the
other end of the line. Ranges of the order of 1—3 km are attainable by standard
instruments with an accuracy of +5 mm and many include slope reduction
and calculation functions as part of the electronics.
Electronic or microwave instruments are mainly used over long ranges (up
to 100 km). The remote instrument needs an operator acting to instructions
from the master at the other end of the line because the signal is transmitted
from the master station, received by the remote station and retransmitted to
the master station. Lengths of 100 km can be measured to an accuracy of
+50 mm.
There are some instruments available that send modulation pulses along the
carrier wave and measure the transit time, but this is not the typical mode
cf measurement.
Optical distance Tacheometry or tachymetry is the method of surveying by which distance and
measurement heighting information can be determined from theodolite observations on either
(tacheometry or a levelling staff (stadia system) or a horizontally mounted bar (subtense system).
tachymetry) In practice, measurement by EDM is replacing measurements by optical
methods but they are still useful techniques when only conventional equipment
is available.
Stadia tacheometry For stadia tacheometry vertical angles need to be measured when the
theodolite’s line of sight is inclined. Stadia marks on the instrument diaphragm
relate the staff intercept to the slope distance.
DISTANCE MEASUREMENT 35
Subtense For subtense tacheometry the horizontal angle subtended by targets on a special
tacheometry horizontal bar (subtense bar) has to be measured. For heighting, a vertical
angle is also required. In Fig. 2.4,
Subtense bar
length b
fs Keeressernney) B
b -
2
pe ye
H
Figure 2.4
b/2 b eat’,
H= = =
tan 6/2 2 tan 6/2 6
when 6 is very small. Example 8.6 illustrates the use of a horizontal subtense
bar.
An alternative to the horizontal bar is a vertical bar carrying two targets;
in this case vertical angles are measured and the method is termed the tangential
system.
2.1 Basic tape A survey line was measured with a tape, .believed to be 20 m long,
measurements and a length of 284.62 m resulted. On checking, the tape was found to
measure 19.95 m long:
(a) What was the correct length of the line?
(b) If the line lay on a slope of 1 in 20 what would be the reduced
horizontal length used in the plotting of the survey?
(c) What reading is required to produce a horizontal distance of
15.08 m between two site pegs, one being 0.66 m above the other?
Solution. (a) Each time the tape is tensioned under its standardized value
a length of 20 m would be booked overall. But actually only a length of 19.95 m
has been covered. Hence the tape is reading ‘high’, the error is positive and
any subsequent correction must have a negative effect.
19.95
Correct length of line xX 284.62
20.00
= 283.91 m.
(b) A slope of 1 in 20 implies that there is a change in height of 1 m over
each 20 m in length horizontally.
In Fig. 2.5(a) angle 6 = tan~! 1/20
22 S2iby avg
Cc 20 B
Figure 2.5(a)
Alternatively, for small angles, @ = 1/20 radian, which gives the same value
for 0.
Thus the required length AC = AB cos 2° 52’
283.91 x 0.998 75
= 283.56 m.
Note that horizontal distances are required when plotting surveys.
5.08 m
15. Q
R 0.66 m
Site pegs R
Figure 2.5(b)
(c) In Fig. 2.5(b) rather than estimating the position of Q above R by eye
or by spirit level, it is preferable to set out the equivalent slope length be-
tween points P and R on top of the site pegs.
Let the slope length required be r.
0.66
Now tana = :
15.08
Therefore a = 2° 30’
thus r= 15.08 sec 2° 30”
= 15.09 m
2.2 Tape corrections The information in Table 2.1 was obtained when measuring the length
of a line by tape suspended in catenary under a pull of 134 N, the mean
temperature being 16°C.
Table 2.1
DISTANCE MEASUREMENT 37
If the tape was standardized on the flat under a pull of 89 N at 20°C
how long is the line?
Cross-sectional area of tape = 3.24 mm?
Mass of tape = 0.026 kg/m
Coefficient of linear expansion = 0.000 0009 /°C
Young’s modulus = 155 000 MN/m?
Mean height of line above mean sea level = 53.78 m
Radius of earth = 6367 km
[Salford]
Introduction. Assuming that the four bays are collinear, five corrections,
sag, pull, slope, temperature and mean sea level, will have to be evaluated
to obtain the actual length of the line.
Solution. It is convenient to tabulate the field measurements as shown in
Table 2.2.
Table 2.2
26 725.54
26 865.22
26 685.33
w2 L?
Sag correction =—-
24 P?
_ (0.026 x 9.806)? x 107 009.67
24 x 134?
—0.0161 m
(P=P)L
Pull correction =
AE
(134—89) 119.632
3.24 x 15.5 x 104
0.0107 m
In this example the simplified slope correction equation is used. A more
rigorous method is presented in Examples‘2.1 and 2.3.
= —0.0010 m
The following computer program can be used to solve this problem. Take
care that the data is entered in the correct units.
Variables
A = Cross-sectional area of tape L Length of any bay
C = Coefficient of thermal L1 = Total length of line
expansion L2 = Corrected length of line
Cl = Correction for tension L3 = Sum of L?
Cz = Correction for sag M = Tape mass
C3 = Correction for temperature N Number of bays
C4 = Correction for slope r Standard pull
= Correction for sea level Pl = Pull used at time of survey
E = Young’s modulus for the S = Height above sea level
tape T = Standard temperature
H = Level differences across any Tl = Temperature at time of
bay survey
= Counter W = Weight of tape
DISTANCE MEASUREMENT 39
210 Cl = INT(((P1 - P) * L1 / (A * E)) * 100000! + .5) / 100
220 C2 = INT((L3 * W * W * (-1) / (24 * Pl **P1)) * 100000! + .5) / 100
230 C3 = INT((L1 * C * (T1 - T)) * 100000! + .5) / 100
240 C4 = INT(C4 * 100000! + .5) / 100
250 C5 = INT((L1 * S * (-1) / (6367000! + S)) * 100000! + .5) / 100
260 L2 = L1 + (C1 + C2 + C3 + C5 + C5) / 1000
270 PRINT "APPARENT LENGTH OF BASE LINE ="; L1; "M"
280 PRINT "CORRECTION FOR TENSION ="; Cl; "MM"
290 PRINT "CORRECTION FOR SAG ="; C2; "MM"
300 PRINT "CORRECTION FOR TEMPERATURE ="; C3; "MM"
310 PRINT "CORRECTION FOR SLOPE ="; C4; "MM"
320 PRINT "CORRECTION FOR SEA LEVEL ="; C5; "MM"
330 PRINT "ACTUAL LENGTH OF BASE LINE ="; L2; "M"
340 END
Introduction. The embankment has been measured along its slope and this
distance needs to be reduced to the horizontal after first correcting for the effect
of pull and temperature. The tape was used on the flat and so there is no sag
correction, and adjustment to mean sea level is not necessary in this case.
Solution. First correct the measured length.
(147
— 49) 29.271
6 xX 207 000
0.0023 m
Temperature correction alot NL
0.000 011 (25-20) 29.271
= 0.0016 m
Slope length, s = 29.271+0.0023+0.0016
= 29.275 m
In Example 2.2 the approximate slope correction £ h?/2L was used. This can
be used for small values of h but would induce a very significant error (236 mm)
on the steep slope that occurs here. The exact reduced horizontal distance can
In = 14.346 X 1.5
= 21.519
Additional road width = 25.519—21.519
= 4.0m
2.4 Incorrect Derive expressions giving the errors in the pull and sag corrections
measurement of field due to an error of +6P in the value of the applied tension P.
pull The length of a baseline was deduced as 1319.774 m when measured
by a tape of length 30 m suspended in catenary. Determine the corrected
length of the line if the actual field tension was 170 N instead of the
intended value of 178 N.
The tape, which was standardized in the flat under a pull of 89 N,
had a mass of 0.026 kg/m and a cross-sectional area of 3.25 mm”. Take
Young’s modulus as 155 000 MN/m? and the acceleration due to
gravity as 9.806 m/s?. [Salford]
Solution. Consider the effects of the error in field pull. Let the error in
the nominal applied tension P be +6P. In the case of the pull correction
P+6P—P,)L
—PJL
actual pull correction = (so
AE
‘ (P=—P)L
nominal pull correction = —————
AE
error = actual correction
— nominal correction
(P+dP—P)L (P—PL
AE AE
'p
= +6P xX ——
AE
DISTANCE MEASUREMENT 41
6P X nominal correction
Therefore the error in = + b
—w-L3 -wJ?
actual sag correction = RAP 2ePy? ~ me
24P? (1 a
FP
w7L3
nominal sag correction = — ———~ '
24 P
error in sag correction = actual correction — nominal correction
wL3 ( | w?L3
24 P? P ) 24 P?
6P
neglecting all other terms in gre
: 6P 4 :
Therefore error in sag = * —— X nominal sag correction,
correction P
_ (178 —89)30
3.25 x 155 000
0.0053 m per 30 m, say.
w2L?
Nominal sag correction =-
24 P?
_ 9.806)?
0.026 x 307
24 x 1783
= —0.0023 m per 30 m, say.
Finally compute the errors in the total nominal corrections. Any change to
the deduced length of the line will depend‘upon changes in just the pull and
Now 6P = —8N.
—6P ;
Error due to —6P = X nominal correction
P-P,
= 2x 8 (0.0023)
—0, m per b ba
178 P Y
I ats 3
= x —0.0023 x 44 for the whole line
178
= —0.0091 m.
2.5 Measuring Prove the following equation for the elongation of a steel measuring
depths down a shaft tape when used for measurements in a vertical plane:
( + 4m (21—x) — *)
g
is the elongation.
is the gravitational acceleration.
the length of suspended tape.
the cross-sectional area of the tape.
the modulus of elasticity.
the attached mass.
&Sym
Ho
3 the mass of the tape per unit length.
DISTANCE MEASUREMENT 43
1 the total length of the tape.
P, the standard tension.
(
1
'
'
Aa
Figure 2.6
Solution. First derive the equation. Figure 2.6 shows the actual suspended
length x of the tape, which hangs from a fixed point C. (/—x) is the additional
portion of the tape not required in the measurement but still contributing to
the loading of the tape. Alternatively, the tape could be considered to be fully
run out to length / and x to be a specific measurement or reading with respect
to C. The tension sustained by the vertical tape due to self-loading is a maximum
at C and varies with y, being a minimum at the longest point under
consideration. Thus the extensions induced in the small elements of length dy
are greater in magnitude in the upper regions of the hanging tape than in the
lower regions but, naturally, all contribute to the overall elongation.
Consider an element of length dy at a point y from the ‘free end’ A.
Load on that element = mg X y
7 -2 x] ae
I $2
(-xyL AE 2
_ mgx (21% )
AE 2 oe
A mass may be attached to lower end A to ensure verticality and to minimize
oscillation. It will have a uniform effect over the tape in so far as elongation
is concerned, and in this context is analogous to the load applied to a ‘horizontal’
tape suspended in catenary between supports.
B
Extension due to mass M over length x = Mg x
AE
The tape has been standardized under tension P, and this has to be allowed
for in the same way as the standard tension in the pull correction mentioned
in Example 2.2.
2.6 Principle of EDM Explain carefully the principles involved in measuring distances using
a sinusoidal wave form as the medium of measurement.
A line AB was measured using EDM. The instrument was set up at
O in line with AB and on the side of A remote from B. The wavelength
of frequency 1 is 10 metres exactly. Frequency 2 is 9/10 frequency 1
and frequency 3 is 99/100 frequency 1.
DISTANCE MEASUREMENT 45
It is known that AB is less than 200 m.
Calculate the accurate length of AB from the phase difference readings
given in metres in Table 2.3.
Table 2.3
: [Bradford]
Double distance = 2D
Outward Frequency, Wavelength
Reflector
Figure 2.7
Solution. First, discuss the basic principles. Figure 2.7 shows three different
frequencies with their wavelengths enclosed within a distance of 100 m.
1
Since \ « —
we have a peg
ae
Therefore \, = fiNy
f
a HY,
10 f, —
Bites oid
Gi-f) f,/l0
At any point within the 100 m length, or stage, the phase of the (ff —/)
frequency wave is equal to the difference in phases of the other two waves.
For example, at the 50 m point the phase of f; is 5 x 27 whilst that of the
jo frequency is 4.5 x 27, giving a difference of 7 which is the phase of the
(f, —fz) frequency. This relationship allows distance to be measured within
100 m stages.
Now, when the third frequency f, is considered we have f; = 99/100f,,
whence \3 = 10.101 m and the wavelength of frequency (f, —f;) = 1000 m.
The same statement in respect of phase differences applies here as well and
further frequencies could be applied to extend the measurements of distance
to cover 10 000 m, etc., without any ambiguity. The term ‘fine’ frequency
can be assigned to f; which will appear in all the frequency difference values
i.e. (f, —f2) whilst the other frequencies needed to make up the stages, or
measurements of distance 100m, 1000 m, etc., are termed ‘coarse’
frequencies. Lengths for 0 m to 10 m are covered by the f, phase difference
measured at the ‘master’ station. In modern EDM the whole procedure is
carried out automatically by the electronics.
Next, we determine the lengths OA and OB. The phase differences have
been measured in metres in this problem and, as shown above, we have
sufficient data to resolve distances from 0 m to 10 m, 0 mto 100 m, and 0 m
to 1000 m.
If n; is the number of whole wavelengths of frequency f;, and n is the
number of frequency f, within OA, then
DISTANCE MEASUREMENT 47
Refer now to the f; frequency, which when related to f, gives the 0 to
1000 m stage. a
Let n,’ be the number of complete waves of f, and n3 be that for frequency
fh, then
Thus there are five whole wavelengths of 10 m length contained within the
0 m to 100 m stage, implying that a distance of 57.386 m is involved. Similarly,
10n,’ + 7.386 = 10.101n3 + 9.911
and nz =n,’ — 1.
Therefore ny) =75)
Thus double length OB = 757.386 m.
757.386 — 374.
Therefore length AB = : = alee)
= 191.525 m.
whilst \ = ee
fn
It will be noted that \ is inversely proportional to fn, in which n is the refractive
index in the air space. This index is influenced by atmospheric conditions,
i.e. air temperature, humidity and pressure, in both microwave and electro-
optical systems. It is very important in respect of the longer lines usually
measured by the microwave instruments.
In the case of electro-optical instruments the dry-air index is of consequence
and this varies through the spectrum. Accordingly a group refractive index
n, is adopted and is then subjected to correction for atmospheric conditions
to give ng.
DISTANCE MEASUREMENT 49
0.000 294 1 760
he 1 +. x —
12 760
ee
( as
= 1,000 281 7.
Note that when f = 14.9854 MHz
cies 299 792 500
> 1.000 281 7 x 14 985 400
= 20.000 m.
At the combination of t = 15 °C and p = 796.50 mmHg
0.000 294 1 - 796.50
i= i+
15 760.00
i+
( aaa)
= 1.000 292 2.
C
Now D « \ « —® in which D is a measured distance. Therefore Dfn is a
constant, i.e. nf
D, fn, = Ds fn;,
where D, is the recorded shaft depth and D; is the depth measured. Therefore
Introduction. Slant distance L between the EDM and the reflector (centred
over their respective stations) has been measured and the reduced horizontal
length, /, is required. To obtain / we require the vertical angle @ which the
~_ 0.004 m
Figure 2.8
measured length makes with the horizontal. This incorporates measured angle
a together with two corrections 6 and y, defined as ‘eye and object’ corrections,
due to the differences in height above the ground stations of the measuring
devices and their targets. The ‘eye and object’ correction and curvature and
refraction have been covered in detail in Chapter 1.
It was also mentioned therein that the measured slope length L between two
stations can be reduced to its chord length at the mean height of the stations.
Accordingly, the equivalent chord length at mean sea level can be computed
and a correction then added to determine the spheroidal distance between the
stations. Assuming that the slope distance was established by EDM, the relevant
corrections are L3/43 R* and L3/33 R* for microwave and infra-red systems,
respectively. These include corrections for the curvature of the path of the
signal and are of particular importance for long lines. Note that the theoretical
difference between spheroidal distance (Fig. 1.10) and chord length is Rc—2R
sin c/2 = d°/24 R?.
Parallel to AB
Measured length 0.172 m
Figure 2.9
DISTANCE MEASUREMENT 51
0.004 x 206 265”
114.652
17%
1.922, — 1.750
From Fig. 2.9 Y= radian
114.652
If dJ/ = + 2mm
0.002 = L sin 0 dé
114.652 sin 4° 30’ 31” dé.
dé = 0.000 221 9 radian.
Therefore dé = + 46”
for d/ = + 2 mm.
Table 2.4
Line of
collimation
Objective
lens
Diaphragm
Figure 2.10
Stadia lines
Horizontal ‘hair’
Vertical ‘hair’
Figure 2.11
DISTANCE MEASUREMENT 53
Solution. Derive basic formulae. Figure 2.10 shows a simple telescope,
whose objective is of focal length f. The image of a staff has been presented
at the diaphragm in such a way that the staff intercept s is given by the difference
between the upper and lower stadia lines in Fig. 2.11. The staff is normal
to the line of sight of the telescope (given by the centre hair reading).
1 1 1
Thus + = —,
u v ‘i
where u is the distance from the staff to the objective and v is the distance
from the lens to the diaphragm.
u u S
most Sy Ute tem Pe My
“3 v i
Therefore u = f +
Ss
1
Figure 2.12
Figure 2.13
angle measured from a vertical pointing such that (on one face) a reading of
90° is given when the line of sight is horizontal. The reading given by the
stadia lines on the vertical staff is s, and in Fig. 2.13 we have the reading
5; which would arise on a staff which is normal to the line of sight. Note that
S; = scos 0.
Now H = D cos 6
me sin
20
+ (f+d) sin 0.
i
(b) Determine distances and height differences. Between A and /, given that
DISTANCE MEASUREMENT 55
fs
H = — cos’ = 100 (1.479—0.500) cos” 4° 37’.
i
Since 6 = 94° 37’ —90° 00’
H = 97.27 m = horizontal distance A—1
i 1.479 —0.500
Pres, Se te ) sin
9°14’
i u 2
= 7.85 m.
Therefore height difference A-—1 = V+h,,—h;, (see Fig. 2.11)
= 7.85+0.990— 1.57 .
= 7.27 m.
Between A and 2 6 = 119° 25’ —90° 00’ = 29° 25’.
Therefore H = 100 (1.633 —1.000) cos? 29° 25’
= 48.03 m = horizontal distance A—2
= 27.08 m.
Therefore height difference A—2 = 27.08+1.318—1.57
= 26.83 m.
The following computer program will carry out stadia tacheometric calculations
by either the staff normal or the staff vertical method. The vertical angle input
in the program is related to zero as the horizontal and the ZD angles given
in this example must have 90° subtracted from them. In lines 20 and 30 the
instrument constants are set to ffi = 100, f+d = 0. For non-standard
instruments these lines must be changed before the program is run.
Variables
A = Vertical angle M_ = Input/output, Minutes
B = Height of instrument P = = Uphill/downhill indicatcr
C = Theodolite multiplying Q$ = Staff normal/vertical
constant (set to 100) indicator
D = Input/output, Degrees S = Input/output, Seconds
H = Horizontal distance V = Vertical distance
I = Stadia intercept X = Bottom stadia reading
K = Theodolite additive constant Y = Mid-crosshair reading
(set to 0) Z = Top stadia reading
L = Difference in level
We now find the standard errors of distances and height differences. This
part of the question is covered by Chapter 8 but, for continuity, treatment is
given here.
Assuming that the multiplying constant will be unchanged, horizontal distance
H is influenced by staff reading errors (ds) and errors in vertical circle reading
(dé).
5
H = — cos’6,
i
dH
therefore —— = As cos?6
ds i
dH
and — = = fs sin 20
dé i
aH \* dH \?
Thus sy? = dH? = (——) s,? + {——) 5s)
ds dé
2 2
ll @ cos’) s+ & sin 28) Se.
i i
where sy, 5, and sg are the standard errors in H, s and @, respectively.
Similarly for V
dV eee sin 26
ds i 2
DISTANCE MEASUREMENT 57
2 dL \2
ry ppm m
2 2
then 5,7 = dL? = e oh ry x 5,7 + ( cos 28) x 59
where sy, is the standard error in height difference L. Now the standard error
of a single stadia reading is 1.5 mm and there are two such readings involved
D0 We .
dh
Sin? = 1D = 2.25 mm’,
Between A and 1
fs sin 20 (100
x 0.979) sin 9° 14’m = 15709 mm.
i
Therefore sy? = (99.352? x 4.5) + 15709?
x (4.3633 x 10~4)?
= 44465.7 mm?
Sy = 210.9 mm = 0.21 m.
pill Se wey
sin Se 1A). si
pope: 2 2
Sy = 0.19 m,
eS 0.09 m.
Solution. Consider the basic theory. In Fig. 2.14 the staff is inclined
through angle 6 towards the instrument and so intercept RS would be given
by the stadia lines rather than the correct value XY on the vertical staff.
Let WV and TU be perpendicular to the line of sight. Since WV and TU
are very near together in practice we can assume that WV = TU. Moreover,
WV makes an angle of (6+6) with RS.
Inclined staff
seats Vertical staff
6 ae
Figure 2.14
Therefore TU XY cos 0
II WV
RS cos (6+6).
Therefore XY true staff intercept
RS cos (6+6)
cos 6
The same relationship holds for elevation 6 when the staff is inclined away
from the instrument.
In Fig. 2.15 the staff is inclined away from the instrument and so WV makes
an angle of (@—6) with RS.
As above TU = XY cos 6 = WV
= RS cos (@—6).
6—6
Therefore XY = RS SAG 9),
cos 8
This relationship holds for elevation @ with the staff inclined towards the
instrument.
We now find the error in the horizontal distance A—1. In Fig. 2.14 we have
6 = 4° 37’, 6 = 1° 00’ and RS = (1.479 —0.500) m.
DISTANCE MEASUREMENT 59
Inclined staff
is Vertical staff
Figure 2.15 2
+6
= Lf x eae) cos? 4° 37’
i cos 6
5° 37’
100 x 0.979 <*>. & cos? 4° 37’
cos 4° 37’
97.12 m.
The apparent distance = 100 x 0.979 cos” 4° 37’
= 97.27 m.
er eT cos (6 —4) :
os? 6
i cos
28° 25’
100° 0.633 a con? 00 8
COSEZ29- 25
48.50 m.
The apparent distance 100 x 0.633 cos? 29° 25’
= 48.03 m.
v 2.00
a 58’ se
2 tan {|———————_-
2
= 28.783 m.
Therefore H = b/60
—b
dH
We
—H?
- dé,
b
i.e. increase in 6 gives decrease in H. Given that the accuracy of measurement
is + 1 second of arc,
1
d@ = ———— radian
206 265
ee 2
Therefore dH = pain Ss in
2 206 265
a O00 9b00
DISTANCE MEASUREMENT 61
Note that dH depends upon H? for a given value d6, and consequently
accuracy is reduced markedly with increasing length.
Problems
Table 2.5
Field temperature = 17 °C
Tape details
Coefficient of thermal expansion 0.000 011/°C
Young’s modulus 207 kN/mm?
Density 7700 kg/m?
Cross-sectional area 6 mm?
An EDM was set up 1.37 m above end A of the base line and a prism
was set up 1.56 m-above end B. The distance between them was read
as 147.859 + 5 mm. What is the maximum and minimum possible length
of the 30 m division on the tape? [Salford]
Answer 30.008; 30.006 m.
Table 2.6
DISTANCE MEASUREMENT 63
dx, are phase differences on frequencies 1 and 2, respectively.
An EDM instrument produces a wavelength of exactly 10 m when
transmitting on frequency 1. Frequency 2 = 9/10 frequency 1 and
frequency 3 = 99/100 frequency 1.
This instrument was set up at O in line with P and Q in order to
determine the length PQ. Station O is on the side of P remote from Q.
Determine the shortest possible value for the length of the line from the
observations given in Table 2.7.
Table 2.7
Explain how this method eliminates the effects of zero error in the
instrument. [Bradford]
Answer 281.601 m
Arg
the formula (x2 —)) where x, and x, are fractions of two dif-
1—A2
ferent measuring patterns with wave lengths \, and ), respectively.
Answer the following questions relating to this formula and its use.
(i) What distance is the formula intended to represent?
(ii) How is the fomula derived?
(iii) What relationship must exist between ), and ),?
(iv) When does the formula fail to produce the required distance
even though x, and x2 have been correctly measured and what must
then be done to obtain the distance? Explain why.
(b) If an instrument capable of displaying the fractions referred to in
part (a) to 0.001 cycle is to be designed, propose with reasons wave
lengths for its measuring patterns which will allow distances up to 10 km
to be calculated by elementary arithmetic, i.e. requiring no electronic
or mechanical aids, with a resolution of 0.01 m. [London]
DISTANCE MEASUREMENT 65
> =| carrier wavelength (0.860 ym),
Co = velocity of electromagnetic radiation in free space
(299 792.5.km/s),
ow Il modulation frequency (14.985 40 MHz). [Eng. Council]
Answer +0.054 m
10 The data in Table 2.8 has been abstracted from the field notes of
a mine surface trilateration survey. Calculate the horizontal distance
between BH and LB. [CET]
Answer 1251.963 m
Table 2.8 ‘
Measured
From| To (m) (m) | vertical angle | (m) (m) (m)
Table 2.9
Table 2.10
Horizontal Vertical
Station circle circle Stadia
14 Three survey stations, F, A and B have been set out in line on steeply
inclined ground. A target at B, 1.219 m above the ground, is sighted
from the two instrument stations A and F. The angles of elevation are
45° 30’ from A and 30° 20’ from F. The height of the instrument axis
at A above the ground is 1.554 m and at F, 1.451 m. The horizontal
distance from A to F is 60.961 m. A levelling staff is held at F, and
a reading of 2.585 m obtained from the instrument at A, the telescope
being set to the horizontal.
Assuming station F to be 24.384 m above Ordnance Datum, calculate
(i) the horizontal distance from A to B, and (ii) the reduced level of
station B to Ordnance Datum. [Eng. Council]
Answer (i) 79.852 m; (ii) 107.509 m
Table 2.11
Vertical Horizontal
circle circle
DISTANCE MEASUREMENT 67
(iii) the heights of points A and B above datum. [Leeds]
Answer (i) 236.02 m; (ii) 281° 06’ 58”; (iii) 226.06 m, 217.96 m
Table 2.13
Stadia | Centre
Instrument} Station hair hair Vertical Horizontal
station | sighted| readings] reading circle circle
The positive and negative signs in the vertical circle column denote
angles of elevation and depression respectively. Calculate:
(a) the horizontal distance YZ,
(b) the heights of points Y and Z above datum and
(c) the error (in mm) which would occur if the slope distance rather
than the horizontal distance of YZ was plotted on a 1/500 scale plan.
[Leeds]
68 SOLVING PROBLEMS IN SURVEYING
Answer (a) 211.79 m; (b) 286.19 m, 265.39 m; (c) 2.0 mm
20 The gradient of the line joining two stations A and B was known
to be 1 in 18. A tacheometer having a multiplying constant of 100 and
zero additive constant was set up at A and, with an angle of depression
of 3°, observations were taken on a staff held vertically at B. If the height
of the instrument axis was 1.190 m and the reading of the lower stadia
line was 1.000 m, estimate the other staff readings, and deduce the
horizontal distance between A and B.
What other method of staff positioning may be adopted? Compare the
two methods with particular reference to the effects of staff tilt.
[CEI]
Answer 1.511; 2.023; 102.05 m
2
21 Use the slope correction expression aa to show that the dif-
ference between the actual correction and the nominal correction has
a magnitude of
h
2 —~ X nominal correction
due to an error of 6h in h.
The slope length of a line was measured as 29.8984 m and / was
booked as 1.382 m instead of the measured height 1.392 m. What error
will arise in the corrected horizontal length of the line?
Answer error = +0.0005 m (correction = —0.0005 m)
DISTANCE MEASUREMENT 69
3
Theodolite and traverse surveying
The theodolite There are three important lines or axes in a theodolite,.namely the line of sight,
the horizontal axis and the vertical axis. The line of sight has to be perpendicular
to the horizontal axis (trunnion axis), and their point of intersection has to
lie on the vertical axis (see Fig. 3.1). The line of sight then coincides with
the line of collimation which, in a correctly adjusted theodolite, describes a
vertical plane when rotated about the horizontal axis. In use, the vertical axis
has to be centred as accurately as possible over the station at which angles
are being measured and that axis has to be truly vertical.
Line of Z
collimation
Circle reading
reflecting
the-direction Vertical circle
Trunnion of line of (face right)
axis collimation Horizontal
circle
Errors in horizontal Errors in horizontal circle readings due to certain maladjustments of the
circle reading theodolite are given below, h being the altitude of the signal observed.
Maladjustment Error
(a) Line of collimation making an angle (90—c) with the csech
trunnion axis, i.e. vertical hair of diaphragm displaced
laterally.
(b) Trunnion axis inclined at (90—() to the vertical axis, i tanh
i.e. tilted from the horizontal.
Errors in vertical It can be taken that errors in the vertical circle readings due to (a) and (b)
circle reading are negligible and that, in each case, the mean of horizontal circle observations
Co-ordinates Normally, plane rectangular co-ordinates are used to identify the stations on
a traverse. A specific point is defined by its perpendicular distances from each
of two co-ordinate axes which are based on north—south and east—west
directions. The former is the reference axis and it can be:
(a) true north;
(b) magnetic north;
(c) National Grid north; and
(d) a chosen arbitrary direction, which could be one of the traverse lines
if so wished.
The intersection of the axes gives the origin for the survey and usually it is
to the south and west to ensure that all points have positive co-ordinates.
Easting/northing The co-ordinates given by the perpendicular distances from the two main axes
are termed:
Figure 3.2
Bearings The position of a point may also be referenced by stating length XY and bear-
ing ¢ of line XY, and these are referred to as polar co-ordinates. The bearing
¢ is termed the whole circle bearing (WCB) of XY. It is measured clockwise
from 0° to 360° at X between the north—south reference direction and the
direction of Y from X. In Fig. 3.3 the whole circle bearing of YZ is 6, and
the whole circle bearing of ZY is (@— 180°).
Similarly for YX in Fig. 3.2 the bearing is (6+ 180°) and, in general, bearing
of line 1—2 = bearing of line 2—1 +180°, 1 and 2 being points within the
system.
WCB of YZ = @
Figure 3.3
Table 3.1
Included angle
Introduction. True north refers to the north geographical pole. The true
or geographical meridian through a point is the trace of the plane through the
north and south poles and the point in question.
Magnetic north does not coincide with geographical north: the magnetic
meridian is the direction revealed by a freely floating magnetic needle. The
angle between it and the true meridian is termed declination.
The whole circle bearing of a line has been defined previously as the angle,
lying between 0° and 360°, between the direction of north and the direction
of the line, measured clockwise.
north A
Figure 3.4
Figure 3.4 shows the traverse survey, the orientation of line BC being 45°
from the meridian. This form of traverse is known as a closed-loop traverse
since it begins and ends at the same point. The sum of the internal angles of
a polygon is (2n—4) right angles, where n is the number of angles.
Thus the sum of the six angles of this example must be eight right angles,
or 720° 00’ 00”, whereas by measurement it is 719° 57’ 00” (Table 3.2).
The total error is therefore —3’ 00” or — 180” , and hence a total correction
of +180” has to be applied. Note that this is of some magnitude and implies
a relatively low order of work, for example fourth order implies an error of
the order of 60VN sec, i.e. 60V6 = 147”.
The corrections can be applied equally to each angle on the assumption that
conditions were constant at the time of measurement and that the angles had
been measured with the same accuracy. Hence a correction of (+180”/6) =
+30” is given to each angle in this example. Next we calculate whole circle
bearings.
It is usual to proceed in an anti-clockwise manner round the traverse when_
internal angles have been measured. To determine the whole circle bearing
of the line to the forward station it is necessary to add the whole circle bearing
of the previous line, i.e. that from the back station, to the internal angle at
the station, and then to add or deduct 180° depending upon whether that sum
is less or greater than 180°. For instance, at A we require the whole circle
bearing of AF knowing that of BA.
N
WCB of BA
%, = Internal
angle at A
Forward
station F
Figure 3.5
Zoe OT 40"
Deduct 180° 180° 00’ 00”
ce Re ee)ten De
Deduct 180° 180° 00’ 00”
Table 3.3
Sexagesimal Centesimal
system system
3.2 Correcting a A survey was carried out on a closed loop traverse with six sides.
traverse by With the traverse labelled anti-clockwise as shown in Figure 3.6 the
Bowditch’s method data in Table 3.4 were obtained.
Table 3.4
Solution. The first stage in the process is to determine the angular error
and apply corrections as discussed in Example 3.1.
Table 3.5 shows the tabulated angular data which sums to 720° 00’ 12”,
Table 3.5
The next stage is to calculate the whole circle bearings. In Example 3.1 a
method of calculating whole circle bearings was presented adding or subtracting
180° to the bearing of the back station plus the internal angle at the station.
The same set of bearings can be obtained by approaching the problem
differently. The reader is advised to try both methods and adopt the one that
is found to be easiest; the two methods cannot be mixed in any one calculation.
Consider Fig. 3.7.
ry WCBg,
WCB,,
Figure 3.7
WCBz, = WCBap + 180° (to lie in range 0° < WCBga < 360°)
WCBgc = WCBg, + B (subtract 360° if greater than 360°)
WCBcg = WCBgc + 180° etc
Table 3.6
Next, the closing error is determined. The algebraic sums of the easting
differences and northing differences should be zero because the traverse starts
and ends at A, but in fact we have total errors in so far as this traverse is
concerned of
AE = 0.067 m
and AN = —0.007m
Closing error = V(AE* + AN?)
= V(0.0677 +0.0077)
= 0.067 m
Expressed fractionally in terms of the total length of the traverse the linear
Displacement —"
due to linear
error
%
4 / |cotectionto ANgc
Vi
“Ci
Correction to AEg,
Displacement
due to angular
error
Figure 3.9
dE x length of line BC
correction to easting difference AEgc . =
total length of traverse
dN X length of line BC
correction to northing difference ANgc =
total length of traverse
in which dE and dW are the total corrections required for the easting differences
and northing differences, respectively.
In this example dE = —0.067 m and dN = +0.07 m. Hence the corrections
for line BC are:
85.771
correction to easting diffe rence —0.067 x
324.572
= —0.018 m
85.771
correction to northing difference +0.007 x ————_
324.572
= +0.002 m
= tan~! ——
dN
This is the bearing of the closing error and the correction in bearing applies
throughout the traverse for all lines.
The corrections can be tabulated as in Table 3.7.
Table 3.7
Correction |Corrected
to AN AN
Co-ordinates
Station
A 1000.00 1000.00
(line AB) — 12.689 6.485
B 987.311 1006.485
(line BC) — 63.136 — 58.074
Cc 924.175 948.411
(line CD) 42.180 — 64.787
4D 966.355 883.624
‘(line DE) 28.019 3.331
E 994.374 886.955
(line EF) 20.730 48.881
F 1015.104 935.836
(line FA) — 15.104 64.164
A 1000.00 1000.00
Variables
A(i) = Angle i Pi — Actual angular misclosure
B = WCB of line 1 to N S = Input/output, Seconds
Bi) = WCB of line i to i+1 S1 = Minimum theodolite
C = 360° in seconds division (least count)
D = Input/output, Degrees T = Sum of measured angles
E(i) = Easting of station i Tl = Expected sum of angles
I = Loop counter U1 = Easting misclosure
K = = Traverse accuracy factor U2 = Northing misclosure
(1-3) U3 = Linear error
L = Sum of side lengths U4 = Fractional error
L(i) = Length of side i to i+1 US = Cumulative easting
M = Input/output, Minutes correction
N = Number of sides_ U6 = Cumulative northing
N(i) = Northing of station i correction
P = Allowable angular X =I+1
misclosure
3.3 Correcting a A closed-loop traverse survey ABCDEA, shown in Fig. 3.10, gave
traverse by the the information in Table 3.9.
Transit rule
Table 3.9
: dE x AE
correction to AEpzp = path ieee 82°)
DAE:
‘ AN
dN x IE
cork i
correction to ANpzr SER
= Rc
L AN
Solution. First calculate the values of AE and AN in exactly the same manner
to Example 3.2. The data are tabulated in Table 3.10.
Table 3.10
Line Length AN
(m) n) (m)
45° 10’ 10” + 207.99 + 206.76
72°04" 65" + 685.87 + 221.77
161° 51’ 45” + 154.75 — 472.42
228° 43’ 10” — 393.28 — 345.27
a ae ETS
300° 41’ 50” —655.11 + 388.94
Table 3.12
E N
(m) (m)
A 1200.00 1200.00
(line AB) 207.97 206.79
B 1407.97 1406.79
(line BC) 685.80 221.80
c 2093.77 1628.59
(line CD) 154.73 ~ 472.36
D 2248.50 1156.23
(line DE) ~ 393.32 ~ 345,22
E 1855.18 811.01
(line EA) ~ 655.18 + 388.99
A 1200.00 1200.00
3.4 Vector A five-sided loop traverse (whose angles have an accepted misclosure)
misclosure of a has been computed giving the co-ordinate differences in Table 3.13 for
traverse each leg.
Table 3.13
Table 3.14
Solution. (i) Calculate the vector misclosure of the traverse. The easting
misclosure is +0.94 m and the northing misclosure is —0.49 m, as indicated
in Fig. 3.11.
: _—pe_E
—0.49 = AN
[ee
0.94 = AE
Figure 3.11
AE 0.94
and the closing error ——- = ———-_ = _— 1,92.
AN —0.49
If we accept that the length of BC should have been booked as 77.54 m (given
that it has been subjected to an error of 1 m) and recalculate the magnitudes
of the misclosures.
3.5 Closed-link Measurements of the traverse ABCDE, as shown in Fig. 3.12, are
traverse given in Table 3.15.
Table 3.15
; \ 123° 16’ 6”
Clockwise angle | Length | .—_.4-)\.
(m)
782.820 mE
260° ° 31’ ‘ 18 ” 129.352 460.901 mN
123°°50" 42” 81.700
233° 00’ 06” ;
101.112
158° 22’ 48” 94.273
283° 00’ 18” ;
740.270 mE
84.679 mN
Figure 3.12
The measured angles are as shown in the figure. Keeping the bearings
XA and EY and also the co-ordinates of A and E fixed, obtain the adjusted
This has to be shared out to five angles, i.e. 21” to A and E and 20” to B,
C and D.
Next determine the easting and northing differences. Having adjusted the
observed angles and deduced the whole circle bearings of the lines, the easting
Line Length AN
(m) im) (m)
129.352 203° 47’ 45” — 52.191 — 118.356
81.700 147° 38’ 47” + 43.721 — 69,017
101.112 200° 39’ 13” — 35.664 — 94.614
94.273 179° 02°21" + 1.581 — 94.260
740.270 84.679
782.820 460.901
and northing differences are calculated for each line on the basis of
AE = | sin (WCB) and AN = | cos (WCB).
From Table 3.16 it will be seen that the total difference is — 42.553 mE and
— 376.247 mN. The fixed co-ordinates for A and E are as in Table 3.17. Thus
the easting differences and northing differences given by the actual
measurements are in error by
eastings —42.553 — (—42.550) — 0.003 m
northings —376.247 — (—376.222) —0.025 m.
Therefore corrections of +0.003 m and +0.025 m are required for the easting
differences and northing distances, respectively.
Make the Bowditch adjustment of the easting and northing differences. As
stated previously the Bowditch linear adjustment for a particular line is
length of line x dE
correction to easting difference
total length of traverse
length of line x dN
correction to northing difference
total length of traverse :
in which dE = +0.003 m and dN = +0.025 m for this traverse. Hence the
corrections to be applied are as in Tables 3.18 and 3.19.
The computation of co-ordinates is given in Table 3.20.
eeLo [oe |
Length (m) 129.352 101.112] 94.273 | 406.437 m
AE correction (m) | + 0.001 + 0.001 + 0.003 m
AN correction (m) | + 0.008 + 0.006 + 0.025 m
Table 3.19
— 52.190
+ 43.721
— 35.663
+ 1.582
Table 3.20
3.6 Errors The magnitude of the errors in centring a theodolite over a station by
introduced when various methods has been investigated and the following results quoted
centring the
(i) plumb bob +3 mm
theodolite
(ii) centring rod +2 mm
(iii) optical plummet +1 mm
(iv) constrained centring +0.1 mm.
Figure 3.13
Introduction. In Fig. 3.13 let the true centring position of the theodolite
be S. It is possible for the vertical axis of the theodolite to lie anywhere within
a circle of radius x from that point, x being one of the centring errors quoted
above. However, there will be two points on the perimeter of that circle at
which the true horizontal angle RST will be subtended. These are S, and S,
which lie on the circumference of the circle containing R, S, and T. Accord-
ingly RS,T = RST = RS,T because all three angles stand on chord RT.
Figure 3.14
Solution. First determine the maximum angular error due to a centring error
of +3 mm. In Fig. 3.14 the theodolite has been centred at U, distance x from
S. The measured horizontal angle RUT = y;,
whilst the correct angle RST = +
=§@+ a
Now yi =%@+90+a+68
x sin 0 sin a
Il
Sint” y Z
It will be noted that the absolute maximum error E is given when sin a and
sin 6 have their maximum values. This occurs when a = 90° and 0 = (y—a)
= 90°, i.e. RST = 180°. .
Note also that E increases as y and z decrease. However, in this case we
are given that y is of the order of 120°, and for the maximum error arrising
; : dE
in this case —— = 0.
a
Now
ee a x sin(y — @) a sin a
sin 1” y Zz
3.7 Errors due to Derive expressions for the error in horizontal circle readings taken
maladjustment of the using a theodolite having the following maladjustments:
th lite . en :
eau (i) the line of collimation not perpendicular to the trunnion axis by
a small amount c
(ii) the trunnion axis not perpendicular to the vertical axis by a small
amount i.
Table 3.21
[Bradford]
; Path of
Path of
line of sight 2, 2 Z, ri # vd
(face right) A (face left)
LJ
S, T,Q,T, S,
\
Vertical Trunnion
circle axis
Figure 3.15
Therefore sin Z = ee
sin ZQ
For small angles we can write
_ SQxsin S
~~ sin(90—h)
in which h = QQ, (the altitude of Q).
SQ sec h
c sec h,
writing SQ = c.
In Fig. 3.16 the left-hand support of the trunnion axis is higher than the right-
hand support and consequently the line of sight sweeps out Z,QS3, making
angle i with the vertical circle ZS3. Q appears to be on that circle but is in
fact on vertical circle ZQ,. Thus the error in the horizontal circle reading for
this particular case is S3Q, and it is negative.
Consider spherical triangle QQ,S; in which Q, = 90°
sin Q,S3 = tani X tan QQ).
Trunnion S Qi os
axis
Figure 3.16
Solution. Determine angle PRQ. Tabulating the errors we have the values
in Table 3.22.
Table 3.22
RP +23” sec 42° 15’ 12” = + 31.1” +15” tan 42° 15’ 12” + 13.6”
RQ +23” sec 28° 12’ 34” = +26.1” —15” tan 28° 12’ 34” —8.0”
Table 3.23
Therefore the corrected angle PRQ = 85° 12’ 53” to the nearest second.
Table 3.24
From the information given in Table 3.24 calculate the closing error
in eastings and northings using the Bowditch method and hence determine
the corrected values of station co-ordinates if the known co-ordinates
of station A are 1500 N and 650 E.
Assume that all the angles, as given in Table 3.24, are correct and
therefore do not require any adjustment. [Salford]
Answer Co-ordinates of C 1543.5 mE, 1927.8 mN
Table 3.25
(i) Compute the traverse, and say whether you would consider the
misclosure satisfactory.
(ii) If you consider the misclosure to be unsatisfactory, say where
the bad observation responsible is most likely to have occurred.
(iii) On the basis that your assumption is correct, obtain adjusted
co-ordinates for B, C, D on the basis of:
ee, ee
Adjust the closing error using Bowditch’s method and hence calculate
the adjusted co-ordinates of A, B, C and D. [Leeds]
Answer Co-ordinates of C 1634.67 mE, 2037.13 mN
Table 3.27
Table 3.28
~E (m)
The mean clockwise horizontal angles and grid distances for the traverse
are given in Table 3.29.
Table 3.29
Whole circle
bearing
Table 3.31
10 Horizontal angle POR was measured, face left and face right, R
being clockwise of P. The mean value of the angle was computed as
75° 30’ 30”, the vertical angles being booked as 22° 00’ elevation on
P and 35° 00’ depression on R.
When measuring the angle it was noted that the horizontal plate bubble
(mounted parallel to the trunnion axis) had moved off centre during the
observations, in Table 3.32.
What value should be assigned to angle POR? The positive signs imply
bubble movement towards the left.
Answer 75° 30’ 16”
P +20" +30”
R 0” +10”
11 List the factors which determine the magnitude of the angular error
due to defective centring of the theodolite.
The centring error in setting up a theodolite over a survey station is
1 mm. Compute the maximum and minimum errors in the measurement
of the clockwise angle ABC induced by this centring error if the
magnitude of the angle is approximately 120° and the lengths of the lines
AB and BC are approximately 5 m and 20 m respectively.
What conclusions can be drawn from this computation?
[Eng. Council]
Answer Maximum error = 47.3”; minimum error = 0”.
Table 3.33
Instrument} Pointing Vertical circle Horizontal circle
station
Face left Face right
Determine the corrected horizontal circle readings given that the line
of sight was not at right angles to the trunnion axis, this being the only
maladjustment.
Answer 27° 24’ 34”; 246° 18’ 33”; 119° 47’ 20”; 338° 41’ 19”
Table 3.34
Horizontal circle Vertical circle
reading (face right) reading
Table 3.35
Altitude level
eS
Left 25° 20’ 40” 25° 21’ 00” 3.5 div. 2.5 div.
Right 25° 21’ 00” 25° 21’ 20” 4.5 div. 1.5 div.
1920 90
2850 260
1830 100
430 390
the co-ordinates and surface level of which are also listed. All quantities
are in metres.
Calculate
(i) the direction and rate of full-dip of the ore-body which may be
assumed to be uniform
(ii) the borehole depth at which the ore-body would be intersected at
point D.
{Eng. Council]
Answer WCB = 174° 40’ 54”; 1 in 3.066; 1205.6 m
Uses of triangulation Triangulation surveys can be used for the accurate location of control points
for plane surveys, aerial surveys and geodetic surveys covering appreciable
areas. They can also be used for the setting out of civil engineering works
such as the piers and abutments of long span bridges, and the measurement
of the deformation of dams.
Control points The control points lie at the corners of geometrical figures such as triangles,
quadrilaterals with diagonals and polygons with central points. All angles of
these figures are measured and their most probable values are determined as
illustrated in Chapter 8.
Base line Providing one side length has been measured all others can be calculated by
trigonometry; this length is referred to as the base line. At least one other length
should be measured, particularly over large areas, so that comparison may
be made between their lengths, calculated and measured. The accuracy of the
work can then be assessed.
If EDM systems are available, and all the sides are measured, a trilateration
survey is established. However, angular measurements define the shape of the
survey network better than wholly linear measurements and so it is likely that
a number of angles will be included in a trilateration survey.
Survey of Great In the geodetic triangulation of Great Britain, commenced in 1936, the primary
Britain triangles had side lengths generally lying between 40 km and 70 km. These
triangles supported a secondary system which, in turn, divided into tertiary
and fourth-order systems. In the latter case stations were at 1 km to 2 km
intervals in urban areas and were often fixed by the techniques of intersection
or resection.
National Grid Triangulation stations in Great Britain have co-ordinates expressed to 0.001 m
within the National Grid. This is composed of lines parallel to a central meridian
(through 2°W) and lines perpendicular thereto, producing a square grid of
side 10 km. The grid was derived from a Transverse Mercator projection whose
origin lies on the central meridian at a latitude of 49°N. The scale is constant
Convergence The direction of Grid North is that of the central meridian through 2°W but
elsewhere a meridian does not align with Grid North. Thus, in general, the
grid bearing of a line will not equal the true bearing of that line if measured
at a station by an astronomical method or by gyro-theodolite. To convert the
former to the latter a convergence factor has to be applied.
4.1. Sine and cosine In a triangulation survey for a certain civil engineering project line
rule AC was measured and found to be 1210.46 m long. Two stations B and
D were established on opposite sides of AC and the following angles
were observed:
ABD = 44° 40’ 59”
DBC = 67° 43’ 55”
ADB = 63° 19’ 28”
BDC = 29° 38’ 50”
Calculate the length BD.
Figure 4.1
re eee tr eG tary
sin 82° 37’ 15” ee es
Finally, use cosine rule in triangle ABC to obtain BD.
Therefore AC? =
ll (0.939 588 x BD)? + (0.498 789 x BD)?
+ (0.357 408 8 x BD’)
1.489 025 BD?.
Since AC = 1210.46 m
ae 1210.46 a
1.489 025
= 991.97 m.
4.2 Satellite station Explain the term triangulation survey and state the factors to be
considered when selecting stations for such surveys.
In triangle ABC, C could not be occupied and a satellite station S was
established north of C. The angles in Table 4.1 were then registered
by a theodolite set up at S.
Table 4.1
A B
p q
Cc
Figure 4.2
First we determine distance SC. In Fig. 4.2 let angles CAS and CBS be
p and q, respectively. SC will be a short distance in practice (less than AC/1000,
preferably) and so p and q will be relatively small angles. From Fig. 4.2
SC AC
sin p sin ASC
SC BC
and ae
sin g sin BSC
Now ASC = 360°—212° 34’ 45”
= Pa 25 157
and BSC = 227° 18’ 12” —74° 30’ 35”
2 es oyi,
: SCam
and si q = —— ‘sin Ia2-47 7
BC
For small angles
4.3 Spherical In a geodetic survey the mean angles in Table 4.2 were observed in
triangles one triangle, each having been observed the same number of times under
similar conditions.
Table 4.2
Triangulation, the curvature of the earth means that such planes are not parallel
at the apices, Fig. 4.3. Accordingly, the three angles of a large triangle do
not total 180°, as is the case for a plane triangle, but to 180° + spherical
excess. The spherical excess depends upon the area of the triangle, and
triangular error is given by
X observed angles — (180° + spherical excess).
Solution. First determine the spherical excess. The sum of the three mean
angles at A, B and C is 180° 00’ 02.7”. In order to estimate the spherical
excess it is necessary to estimate the area of triangle ABC. For this purpose
it is sufficiently accurate to assume that the triangle is plane. The three angles
should thus sum to 180° and we can deduct 2.7”/3 = 0.9” from each to give
AWS 162°a 246 75 ¢
B = 64° 56’ 09.0”
Co 52395533554
37 269.280 “fi BC
sin 52° 39’ 33.5” sin G2 °924) 91725)”
Therefore BC = 41 544.469 m
Area of triangle
R* sin 1”
where R is the radius of the earth.
701.3
Whence the spherical excess of triangle ABC =
6383.3937 x sin 1”
= 3.6”.
Now calculate the triangular error and the’ angle corrections.
Table 4.3
4215.65
205° 36’ 12”
3778.46
144° 23’ 20”
5237.28
Figure 4.4
The bearing of XQ is 36° 40’ 18” and the latitude of X is 52° 20’
45” N. Determine the bearings of YX and YR at Y and estimate the
latitude of Y.
Length of 1” Length of 1”
of latitude of longitude
30.9022 m 18.9364 m
30.9107 m 18.9008 m
Ps
Figure 4.5
Yp, = 90-9,
Cos
(Pe) r Pn
tan a COL
eatey 2
Ccos| SS
2
COs
(90 — d2) — (90 — 4)
a+(180—a—6a) 2
tan ee X COL
2 :(90 — $1 0- $1) -
Ff os2
wa(m0—) =
= ——
- —— X cot LS
“: = as nN N
ar
da = 6 sin(**2)
= @sind
where ¢ is the mean latitude of X and Y.
Figure 4.6
dX X 6¢ = I cos (a+)
p x 60 II I sin (a+——
A}
and ba II 60 x sin @.
Note that had X and Y been widely separated, triangle Xp, Ycould be solved
for angles X and Y using the standard expressions for tan [(X+Y)/2] and tan
[(X = Y)/2).
Solution. First calculate the traverse details. In the first instance the
relationship betwen X and Y is established as if the surface of the earth is plane.
Table 4.6
Easting Northing
difference difference
(m) (m)
4215.65 36° 40’ 18” 2517.71 3381.25
3778.46 62°°16"30" 3344.66 1757.85
5237.28 26° 39’ 50” 2350.26 4680.32
(b) The easting and northing differences are given in Table 4.6.
Length of XY II= V(8212.637 + 9819.427)
12 801.11 m = /(Fig. 4.6).
es 8212.63
Bearing of XY
9819.42
39° 54’ 28.7” =a.
Next determine the convergence effect. Had the latitude and longitude of X
and Y been given, 5¢, 60 and @ would be known and the expressions derived
previously could have been used directly. However, the latitude of Y and the
longitude difference 60 are unknown and accordingly the mid-latitude of X
and Y has to be estimated and 6¢, 60 and da evaluated by successive
approximations.
For the first trial we assume that the mid-latitude of XY is that of X, i.e.
52° 20’ 45” N.
45
Therefore 30.9022 + 300 (30.9107 — 30.9022)
30.9035 m.
The approximate difference in latitude between Y and X
= (52.923:23.9"..
203.9
= II 18.9364 — x 0.0356
18.9122 m
Isina __ easting difference
‘anda 00) =
18.9122 18.9122
434.25”.
Therefore ba 60 xX sing
434-25 Sie 25 25.9.
344.01” = 5’ 44.01’.
For the second trial we take the mid-latitude to be 52° 23’ 23.9”
203.
A = 30.9022 + eae (30.9107 — 30.9022)
300
= 30.9080 m.
Now da = 344.01”
$e _ 172.00".
2:
l ba
Hence 6g = — cos [a + —
r 2
odWi
= ne se: cos 39° 57’ 20.7”
30.9080
a= SlyLahe.
317.48”
Revised mid-latitude = 52° 20’ 45” +
ll S2metige ls! ye
203.7
Wh ence yb = 18.9364— (300 x 0.0356 )
== PR) Pease
Tes: ba
Revised 66 = — sin {a + —
m 2
It will be seen that only small changes have occurred in the values of 6¢
(decrease), 50 (increase) and 6a (increase). A further trial is not warranted,
particularly in view of the fact that the data were expressed to an accuracy
of a single second.
x
Whence 6a = 344.35
= 05’ 44”, say.
Therefore bearing of RY at R = 26° 39’ 50” + 05’ 44”
= 26° 45’ 34”.
Therefore bearing of YR at Y = 26° 45’ 34” + 180°
= 206° 45’ 34”.
Similarly,
bearing of YX at Y 39° 54’ 29” +05’ 44” +180°
220° 00’ 13”.
latitude of Y = latitude of X + 6¢
= 52° 20’ 45” + 317.48”
=. 52° 26’ 02” N.
4.5 Co-ordinates Two survey stations A and B have the National Grid co-ordinates and
latitudes given in Table 4.7.
Table 4.7
Calculate
Figure 4.7
Figure 4.7 shows a spheroid with point A on its surface. The normal at A
meets the minor axis at C and AC represents the radius of curvature (r)
perpendicular to the meridian. @ is the geographical latitude of A, and AD
is the radius of the parallel of latitude through A, having a value of r cos ¢.
Hence the length of 1” of longitude at this latitude is r cos @ X sin 1”.
Table 4.8
er
ae ee B 460 476.691 350 340.727
A 460 257.664 350 258.130
tan QAB
Ngp—-Na
219.027
82.597
2.651 754 9.
Therefore aap = 69° 20’ 17”.
(iv, ind the convergence at station A. All north—south lines on the grid
are pa lel to the central meridian running through 2°W, the true origin lying
at 49°N thereon. However, to ensure that all eastings are positive and all north-
ings are less than 1000 km a false origin 400 km west of and 100 km north
of the true origin was adopted. Hence the grid easting of the central meridian
is 400 km, so that A is at a perpendicular distance of (460 257.664 —
400 000) m = 60 257.664 m from the central meridian. Referring to Fig. 4.8,
let L be the perpendicular distance from the central meridian to A.
ike Meridian
at B
Meridian
Central
atA
meridian
Figure 4.8
L
= 206 265 — tan o”.
ic
Since the north—south grid lines are parallel to the central meridian, 5a is
the convergence, C,, at A.
60 257.664
Therefore Cy, = 206 265 x ad tan 53° 02’ 45”
6 384 100
= 2587.9”
= 43’ 07.9”.
(v) Find the convergence at station B. The perpendicular distance from the
central meridian to B is 460 476.691 — 400 000 = 60 476.691 m. Therefore
Grid north
Meridian
atA Grid bearing BA
Central
Meridian
meridian».
at B
Figure 4.9
Note that when the station is to the west of the central meridian, i.e. its easting
co-ordinate is less than 400 000 m, its meridian will be inclined in the other
direction with respect to the central meridian. Convergence (C) takes a negative
sign and is deducted from grid bearing to give geographical bearing. A further
point in respect of long lines is discussed in Example 4.6.
4.6 (t—T) correction Calculate the (t—T7) direction corrections to obtain straight-line
directions at two stations A and B having the National Grid co-ordinates
A E 300 120 m N 385 920 m
B £E310970m N 380 230m
Aide-mémoire
Grid bearing
of BA
Projected
sight line
Figure 4.10
4.7 Local scale The measured slope distances from a survey station A situated on a
factor colliery headgear 101.15 m above Ordnance Datum, to two Ordnance
Survey stations B and C, are AB 1923.400 m and AC 1398.446 m. The
National Grid co-ordinates and levels of stations B and C in metres are
Horizontal
plane
Figure 4.11
AA, = V(1398.446?
— 158.847)
= 1389.396 m.
Note that since the height differences are appreciable when compared to the
slope distances Pythagoras’s theorem has been used to determine the horizontal
distances, rather than the slope correction —h7/2L referred to in Chapter 2.
Grid distances can now be determined. In order to calculate the correspond-
ing grid distances between A and B, and A and C we must initially establish
the distances at mean sea level. The horizontal distances AAj, etc., are
positioned at their mean heights, h,,, as shown by the dotted line in Fig. 4.12.
Vertical Vertical
at B
Figure 4.12
Sap = 1902.248 m
and Sac = 1389.357 m.
Now FS = s.
Therefore saz = 0.999 672 x 1902.248
= 1901.624 m
and Sac = 0.999 672 X 1389.357
= 1388.901 m.
The calculations for grid distances have been presented in stages for clarity,
but these distances could be readily determined directly from the horizontal
distances L by the expression
I SON sl
ame) : ;
Table 4.9
Station E N
(m) (m)
C 324 022.07 342 846.89
B 323 679.35 340 431.32
Figure 4.13
b? = a?+c?-2 ac cos B,
whence 2 ac cos B = a?+c?—b?
a* = 2439.767 = 5 952 428.858
c? = 1901.62? = 3 616 158.624
add = 9 568 587.482
b* = 1388.90? = 1 929 043.210
deduct = 7 639 544.272
7 639 544.272
Therefore cos B =
2 X 2439.76 x 1901.62
0.823 316 1.
Therefore B44" 34" 54.7".
The grid bearing of BA is 34° 34’ 54.7” + 08° 04’ 30.7” = 42° 39’ 25.4”.
Whence between B and A
Mid-position
a, | of oscillations
(Horizontal circle
reading = N)
as Turning point
(right)
a3
Plot of movement
Of gyro mark
Turning point
as (left)
Figure 4.14
151° 44 84
5 ; x ° 5 ’
Soe te
° ; ,
peal, £38? To ag
4
4.9 Gyro-theodolite Horizontal circle readings on a reference object were recorded as 297°
observations — 32.8’ and 117° 32.6’ face left and face right, respectively, and an
transit method approximate gyro orientation of 159° 37.0’ was then established.
Successive transit times and amplitudes were recorded as follows:
a
Determine the azimuth of the reference object given that the proportion-
ality factor was 0.046’/s and the instrument calibration factor was
—2.5’. The gyro attachment was fixed with the telescope in the face-
left position and the gyro scale is marked positive and negative left and
right respectively of the zero mark.
Figure 4.15
Figure 4.16
ON = ca’ 6t,
in which c is the proportionality factor, a’ is the amplitude measured in scale
units, and df is the time difference derived from the times of three successive
transits. Reference may be made to ‘Surveying’, Bannister, Raymond and Baker
(Longman) for the derivation of the above expressions. The circle reading,
N’, is corrected by 6N to achieve the true meridian pointing, N.
Table 4.10
03 min 21.1s
07 min 13.8 s
10 min 34.7 s
Nt
Transit
times
noted
Figure 4.17
This has been booked as a positive time of swing: the gyro mark was to the
left of the scale zero during this period of time.
(b) 07 min 13.8 s—03 min 21.1 s = 03 min 52.7 s = 232.7s.
The corresponding amplitude was negative, i.e. to the right of the scale zero
and hence this particular time of swing has a negative sign.
(c) 10 min 34.7 s—07 min 13.8 s = 03 min 20.9 s = 200.9 s.
As in (a) this is a positive time of swing.
The differences in successive swing times can now be calculated to give dt
values of
+201.1—232.7 = —31.6s
and +200.9—232.7 = —31.8s,
respectively.
Now ON =ca’ 6,
Note that since the gyro attachment was positioned with the telescope face
left, the horizontal circle reading for the fixed pointing N’ was observed on
that face and the corresponding face reading had to be taken for the mean
pointing on the reference object.
ee ae
awAtw ate apAte
Surface observations
Known azimuth 128° 17’ 52”
Underground observations
Observations west of gyro north
Horizontal-circle reading during transit observations: 10° 08’ 00”
Horizontal-circle reading to reference object: 49° 44’ 15”
Time oftransit: 0 m00.0s, 3m08.6s, 7m26.7s, 10m 36.3 s,
14 m 54.8 s
Average amplitude: 4.8
6N
therefore c = ———————_
aw Atw + ap Atg
From the surface observations west of gyro north we obtain the data in Table
4.11.
Table 4.11
0m 00.0s
2m 35.5s
7m 26.4s
10m 02.2s
14m 54.5s
Mean At = 135.67
0m 00.0s
2m51.3s
7m 27.38
10m 18.4s
14m 53.9s
Mean At = 104.67
0m 00.0s
3m 08.6s 69.5
7m 26.75 68.5
10 m 36.3 s 68.9
14m 54.8s
Mean At = 69.0
Table 4.14
Time difference
At (s)
0m 00.0s
3m 12.8s
7m 27.0s
10m 39.3s
14m 54.0s
Mean At = 62.0
azimuth = M—N+E
= 49° 44’ 15” —(10° 08’ 00” + 14’ 54”) + 5’ 13”
= 9992267 .247
Table 4.15
x 52° 00’ 49.24” N | 00° 46’ 29.65” W | 484 176.101 | 235 737.866 | 81.141
Y 52° 01’ 43.24” N | 00° 46’ 20.44” W | 484 323.566 | 237 408.841 | 76.879
Solution. The first step in the calculations is to determine the grid distance.
From Table 4.16,
147.465 1670.975
Next use the GPS data to calculate the lengths of 1” of latitude and 1” of
longitude at the mean latitude of the two stations. Referring to Fig. 4.6,
Mean latitude of X and Y 52°01 16.24”
ll S2e0TE: 24"
Difference in latitude 66 = 54.00”
Difference in longitude 69 = 9.21”
Referring to the equations derived in Example 4.4, the length of 1” of latitude
at the mean latitude is
= 30.9049 m
and the length of 1” of longitude at the mean latitude is
30.9049 x 54.00
Therefore, xy =
cos 06° 00’ 23.95”
ll 1678.078 m
= 1678.057 m
and grid distance, s = FS
= 0.999 69 x 1678.057
= 1677.537 m
This compares to the previously calculated distance of 1677.469 m, i.e. an
overestimate of 0.068 m.
In practice the software converting latitude and longitude to grid co-ordinates
will use a more rigorous mathematical approach, leading to a greater degree
of accuracy, than the method used here which assumes that the meridians at
X and Y are parallel.
Problems
Table 4.17
Table 4.18
[ Pointing Horizontal circle reading
A 9927.15"
B 71° 04’ 42”
Cc 307° 38’ 15”
Table 4.19
Latitude Length of 1’ Length of 1”
of latitude of longitude
30.9155 m 18.5065 m
30.9160 m 18.4704 m
[Salford]
Answer 137° 08’ 49”
14 Given that the co-ordinates of two survey stations A and B are ea,na
and eb,nb on a local arbitrary grid and EA,NA and EB,NB on the
National Grid, derive expressions for the co-ordinates of the original
of the local grid relative to National Grid origin and for the sine and
cosine of the angle between the two sets of co-ordinate axes.
If the co-ordinates of A and B relative to the two grids are as tabulated
below, derive the National Grid co-ordinates of point P from the values
of the local co-ordinates given.
Point A ea 374.62 na 615.88
EA 49 624.31 NA 21 315.74
Point B- eb 47.23 nb 566.37
EB 46 676.24 NB 20 988.72
Point P_ ep 373.19 np 127.38
[Bradford]
Answer 50 111.34 mE, 21 517.92 mN
Table 4.20
Easting Northing
(m) (m)
Castle Hill 2619.43 2803.03
Beckwith Mount 4142.49 2005.50
The approximate co-ordinates of the rig are scaled from the Ordnance
Survey sheet as 3662 mE, 4400 mN. The distance from station S to the
rig is carefully measured as 18.835 m and the following angles observed:
Beckwith Mount — Castle Hill - S 85° 03’ 11”
Castle Hill — Beckwith Mount - S 50° 36’ 32”
Rig — S — Castle Hill 108° 30’ 43”
Rig — S — Beckwith Mount 64° 10’ 26”
Table 4.21
E N
(m) (m)
| A | 449 674.24 321 315.74
449 624.31 320 988.72
Given that the local scale factor is 0.999 63, determine the spheroid
distance between A and B.
Answer 331.24 m
18 A straight line AB, 18.5 km long, was set out at 90° to the meridian
at A in latitude 52° 22’ N. Determine the bearing of AB at B if 1” of
longitude has a length of 18.9222 m.
Answer 90° 12’ 54”
Table 4.22
E N
(m) (m)
| A | 527 398.249 304 854.427
551 477.378 315 465.283
(i) the National Grid distances from station R to the other three stations
Table 4.23
Table 4.24
d*d"d"d
Figure 5.1
0, +0
Trapezoidal rule Area = d (27% + on + 04 ee + o,-1}
In the trapezoidal rule the area is divided into trapezoids, boundaries being
assumed to be straight between pairs of offsets. The areas of individual
trapezoids, i.e.
0; +02 of 07 +03 a
2 2
are added together to derive the whole area.
d
Simpson’s rule Area = ci (0) F°40y F203 A 20,4 FF 0,):
The areas of successive pairs are added together when formulating the rule.
Since pairs of intercepts are taken, it will be evident that an even number of
intercepts is required when using the rule to determine an area. If an odd number
of intercepts is present then the first or last intercept is treated as a trapezoid.
Areas of straight- In so far as straight-sided figures are concerned there are standard expressions
sided figures available, for example:
Triangle area = +ab sin C,
in which C is the angle included between sides a and b.
4 a+b
Trapezium area = 3 Pp;
Planimeter The planimeter, an integrating device, is available for the direct measurement
of all shapes, irregular and regular, and high accuracy can be attained.
(a) Level across (b) Two-level section
Original Tink
ground d d.
icici fee ey = 1 2 = is
| |
“ins
Formation -—*~ bib 1ins
level a2 Original
ground surface
ee | 4
b
2
ny|o
(c) Three-level section (7) Part cut—part fill
Figure 5.2
Rv 2,
ayt 2 k+s :
|): ie Pe?a
a51\9P + sh) (dq,1 + d 2) > |"
= h(b + sh)
which is the ‘level-across’ case in Fig. 5.2.
Having computed the cross-sections at given intervals of chainage along the
centre line by standard expressions, as above, or by planimeter, etc., volumes
of cut in the case of excavation or volumes of fill in the case of embankments
can be determined. The end-areas rule or the prismoidal rule are applicable,
and they are analogous to the trapezoidal rule and Simpson’s rule, respectively,
for areas.
We will consider cross-sections which are D apart.
This method gives the correct volume providing that the area of the cross-
section midway between two cross-sections D apart equals the mean of the
two cross-sectional areas.
2D
Prismoidal rule Volume = am [Age 4Ap ot 2Agechs .. + 24,2; +-A,].
The prismoidal rule is based on the assumption that the earth forms a prismoid
between the two cross-sections 2D apart. For this to apply the linear dimensions
of the section midway between them have to be the mean of the corresponding
dimensions at the outer sections. It is accepted that the prismoidal rule is the
more accurate, but in practice ‘end-areas’ is frequently adopted for estimation
purposes, since some irregularities are likely to occur between various cross-
sections. In addition, there is the problem of bulking and shrinking when soil
is removed from its virgin state and later replaced or placed elsewhere.
Quantities for the construction of mass-haul curves are normally established
using ‘end-areas’.
Contours and spot Contours are used in a manner similar to cross-sections when calculating
heights volumes, the contour intercept value for the scheme being adopted for D. Spot
heights have a role to play in the computation of volumes of excavations for
basements, etc. A network or grid of levels is taken over the site, the net or
grid being such that the ground surfaces contained within the grids can be taken
5.1 Scaling from (a) On a certain Ordnance Survey sheet an area of 8.965 hectares on
plans with an the ground covers 143.44 cm?. What’is the scale?
allowance for (b) The plan of an old chain survey, plotted to 1/500 scale on linen
shrinkage cloth, was found to have shrunk such that a line originally 100 mm long
was now 98 mm. In addition, a footnote on the plan stated the 20 m
chain used during the survey was later found to be 20.02 m long after
completion of the plot. If a certain area on the plan is measured by
planimeter to be 2143 mm? estimate the correct area on the ground.
Assume uniform shrinkage.
5.2 Using Simpson’s State Simpson’s rule and the trapezoidal rule for the determination of
rule and the areas. Figure 5.3 indicates a field with two straight boundaries AB and
trapezoidal rule BC and an irregular third boundary AC. The lengths of the straight lines
between the stations and the offsets from AC at defined chainages from
A are as follows: ‘
Figure 5.3
Solution. (a) First calculate the area of triangle ABC. A standard expression
for the area of a triangle whose sides are known is
area = VS(S—a) (S—b) (S—c)
a+b+c
in which a, b, c, are the lengths of sides and S$ = 5
(b) Next calculate the area of the irregular figure. Simpson’s rule requires
an even number of increments in its application, whereas the trapezoidal rule
Table 5.1
110 x 71.5
= 7865 m?.
Total area enclosed is:
5.3 Area inside a Show that the area, abcd, enclosed by a traverse having four stations
traverse a(x; ¥1), b(x2 y2), C(x3 y3) and d(x4 y4) is given by the equation
Figure 5.4
Solution. Establish the expression. Figure 5.4 shows the four stations abcd,
lettered in a clockwise manner. By inspection: Area abcda = Area (abfea +
bchfb — dchgd — adgea). Now
+
Area bchfb = (
--
Area adgea = ( )
= [Ot tyox3+¥3%4t+Y4r})
— (YpX4 + y2x1 +y3X2 + 4x3)].
Next determine the area ABCDEA. In this example the basic expression
derived above must be extended to include the fifth station, i.e.,
Area = 5[1 (x. —x5) +y2 (x31) +3 (X44 —Xp) +94 (Xs — 3)
+5 (x1 —X4)]
in which x and y refer to easting and northing co-ordinates, respectively. Figure
5.5 shows the five stations together with the co-ordinate axes x and y. It is
convenient to tabulate the data with respect to the above expression, as in Table
3.3:
Therefore area = +2 yp (X41 — Xn—1)
= 4 x 290 000
Table 5.3
0 (+150 0
300 (+450 +135 000
70 (+320 + 11900
—130(— 60 + 66 300
— 240 ( 0 + 76800
{ +290 000
Variables
A Traverse area
I Loop counter
N Number of sides
X(i) X co-ordinate of station i
Y(i) Y co-ordinate of station i
ef wilhwasEP |
350 500 900 650 470
400 650 650 150 270
[Bradford]
In this example there are two unknown sets of co-ordinates which depend
upon lengths r and s in Fig. 5.6. The solution is based, first, on the relationship
between r and s and then, second, using that relationship to evaluate the co-
ordinates of either P or Q, knowing the area PCDQP.
Figure 5.7
Table 5.5
Table 5.6
Now, area
5.5 Cross-sections Prove that on uniformly sloping ground the cross-section area of a
and volumes from cutting or an embankment is given by the following formula:
cutting and
embankment details
where d, and d, are the horizontal distances from the centre line to the
limits of the side slopes, b is the formation width of the cutting or
embankment and 1 vertical in s horizontal is the gradient of the side
slopes.
An embankment for a new road is to be 20 m wide at the top, with
side slopes of 1 vertical to 3 horizontal. The heights of fill at the centre
line of three successive cross-sections, 50 m apart, are 2.4 m, 3.2 m
and 4.0 m respectively and the existing ground has a uniform cross fall
of 1 in 12. Calculate the volume of fill required. [Salford]
bt Figure 5.8
1 ese
Introduction. In Fig. 5.8 the side slopes are shown meeting the original
ground surface at A and B so that the side widths are d, and d,, respectively,
and hy and h, denote the differences in height between these points and the
formation level. Points E, F and G, lie on the centre line of the embankment.
b
Now d, —> + sh 2
and
d,
hn = h— °% where h = FG.
Therefore d, =
|
ll
nv
|e
Se Fee’
aa
n|o
—~ Sa
| a|&|>
a x Se
7 ee
and hy T +.
|
b
Therefore d,; = & + sh met ;
2 k-s
ACDBA is known as a ‘two-level’ section because two observations of height
will define cross fall 1 in k.
and EF 2
2s
did,
S 4s
1 b?
— |did),
> (aie si
: b k
In the expressions ag =\(— + sh) |——
a (Soa)
2 k—s
and A=| So |
|e
o
eet a(S
‘a
Cross- — + sh
section
(m)
Ss A 103.24)
50
In this example
Diife ke
='s
Figure 5.9
Solution. First determine the side widths. The two relevant expressions
for side widths are
k
d, = at = sh :
2 k,+s
d by + sh
d, = (— ky
* ; (; ) (G25)
in which : 6.5 m, kj = 10, s = 2.5 andh = 1.0m
16
a) =(8.5 + 2.5 x1.0) (———+—) =913.04 m.
16—2.5
Next determine the area of section 1.
E :
BJ = BE + EJ = BE + s ii
S Le
= 44m
EJe— 34m:
2 x 20
Therefore Vp = ———— [18.8 +. (4 x 20.3) + 2 X 27)
+ (4 Xx 28.8) + 30.2]
Vp = 1925.3 m?, say 1925 m>.
A,+A
Ven = D(A2AS 4 ty + Ay + As)
20 =26 4b j30,
+ 20.3 + 21.7 + 28.8)
1906 m°.
1ins | Cc
i —_ eS in r
t d1 + d> —
Figure 5.10
Solution. First calculate the sidewidths and areas. Consider the difference
in level between formation ABCD and points E and F which define the limits
of the sidewidths.
: d,—kh
Rise from
C to E = hy = k
= Rise from D to E
z d,—b/2
ayer ;
d,—b/2
Therefore — r
Woolies
b k
d, = {|— —-rh) |——}.
b k
Therefore d, = + sh)yah
2 k-s
Area in cut 7x CD x hy
N|-
ua das
N|—
oes tae
omy (ert)
Area in fill N|—
x AC X h,
“(+a @)
N|-
Note that if h had been the depth of cut at the centre line rather than fill, the
insertion of —h for h in the above expressions would give the respective
formulae for side widths and areas.
Next, calculate the volumes of cut and fill. As mentioned previously the
relevant volumes will be determined initially by end-areas.
Table 5.8
Section h
(m)
A,+A
volume of cut between X and Y = D eras
: = 0. =)
~ 10 'e
= 9dfeof rhe
os (see)
volume of fill between X and Y
2
142.3 m?.
Establish the prismoidal correction. Consider two parallel cross-sections,
triangular in shape, D apart as in Fig. 5.11, their areas being A; and A). Also
consider the cross-section midway between them, formed by the straight
connecting them. :
A, = area of mid-section
1 (bith hy +h,
RS pas
The volumes contained between the outer cross-sections are: by the end-area
method
ies D |
(b,h byh
2 2 2
By the prismoidal rule
Vp = D byh, a 4 b, +b, h, +h, ‘s bshy '
hall eae’ Z a 4 2
Now, the prismoidal correction or prismoidal excess is Vg, — Vp
8 D [(buf Gp Dalia
bsh ie Digebytage is, ofBirtby s
2 2 2 6 # 2
hithy\ | bah
2 f,
D
= — (b,—b,) 2) (hy
(oe! (h,—h).
— hy)
D ( difference By paige |
12 \ base lengths heights
for the triangular section. Using the expressions previously derived for the
cut section
d)—b/2
and height h, =
r
, ; b k
in which d, = {— — rh ‘
2 k=r
Now, b/2 = 7.50 m, k = 10 and r = 1. Whence, on tabulating, we get the
values in Table 5.9.
Table 5.9
b b
Cross-section h = —kh zs—rh d, hy,
20
Therefore prismoidal correction = re (5.00 — 3.70) (0.56 — 0.41)
= 0.33 m’.
Therefore volume when treating the earth between the two sections as a
prismoid
= 21.5 — 0.3
= 21.2 m’.
Similarly for the fill
b
base length = AC = ry + kh
d,—b/2
height =h, = pe ie
S
b k
in which ad, = (— + sh| |——}.
e kas
With s = 2 we get the results in Table 5.10.
20
Therefore prismoidal correction = % (10.00 — 11.30) (1.25 — 1.41)
0.35 m?.
Therefore volume of fill Il 142.3 — 0.3
142.0 m°.
Cross-section h
(m)
X 0.25
~ 0.38
5.8 Calculating In a cutting the width of the formation is 8 m and the side slopes are
volumes along a 1 in 2. The depths at the centre line, which lies on a circular curve of
curved section of radius 150 m, at three cross-sections 20 m apart are 2.50 m, 3.10 m
road and 4.30 m, respectively. The cutting is to be widened by increasing
the formation width to 11 m, the excavation being on the outside of the
curve and retaining the original side slopes.
Calculate the volume of excavation between the cross-sections using
the end-areas method. Assume that the transverse slope of the ground
is 1 in 5 at each cross-section.
Solution. First amend the end-areas and prismoidal formulae. If the cross-
sections are unsymmetrical, e.g. have transverse ground slopes, but have the
same shape and dimensions, the path of the centroids will be parallel to the
centre line. When there is some variation, as in this example, the path is not
parallel to the centre line. Some close approximations are made when using
either the end-areas of prismoidal formulae. For example, consider the three
successive cross-sections, D apart on the centre line, in Fig. 5.12, the centroids
of the ‘widening’ areas being e;, e) and e3, respectively, from the centre
line. The mean distances of the centroids of the two pairs of sections are
(e,+e,)/2 and (e,+¢3)/2 from that centre line. These values are taken to
represent all respective centroid positions within the relevant cross-sections
and therefore define the path.
If 6 is the angle subtended at the centre of the curve by the pairs of equi-
distant cross-sections
me \ | peeks eee
| |
babi
ioe de) |
Figure 5.12
6=—.D
R
: e,t+e
Therefore the distance between centroids = 0 (k= s) and
Had the widening been in the inside of the curve the above would change to
Ot Sis id]
5 2. ;
Figure 5.13
>
Figure 5.14
Therefore a= =.
Sinilafyse ett se ee
5 9)
+2
Therefore = SS
Also 2h = d—-7.
10H+14
Therefore y= here
Also QRTS = 3h, m’, and so we get the values (in metres) given in Table
5.11. It will be noted that (d—w) is constant, and so triangle PQR has fixed
dimensions, its height PV being 5.0/k = 1.0 m, with QR = 3 m.
Table 5.11
2h
d- (4.0+ 1.5 + =a =(d—h, —5,5)m
from PW.
(1.5 x 4x 3.5) +3h; (d—h,—5.5)
Therefore X =
1.5 + 3h,
From this we get the values in Table 5.12. We can now calculate the volume
of excavation.
Table 5.12
2 2R Z 2R
(20
18.0+21.0
2
(1eameteree|
1.56+12.55
2 x 150
+
21.0427,7.0 saan : ;
A ink 12.55+14.54
a 2ix 150
= 944.7 m?, say 945 m°.
Alternatively, the areas can be amended as in Table 5.13.
Table 5.13
A(i + e/R)
(m*)
19.39
22.76
29.62
Had the influence of curvature been neglected the calculated volume would
be 870 m?.
5.9 Volumes from a Figure 5.15 shows the distribution of 12 spot heights with a regular
rectangular grid of 10 m spacing covering a rectangular area which is to be graded to form
spot heights a horizontal plane. The spot height data (in metres) are:
VLG; 2 aay 3 17.63 1
A ehhS7s to Dally! 6 17.96
Del S8ruley Bot BO bie O. 418,25
10 17.83; 11 18.19; 12 18.42
Figure 5.15
Introduction. The solution to this example will be based on the fact that
the volume of a triangular prism is equal to the area of its normal section
multiplied by the mean length of side. Dividing the square grid of Fig. 5.15
into triangles, as shown, produces a set of vertical triangular prisms, one plane
end being at formation level and the other being the ground plane contained
between the corners of each triangle. Squares could be used rather than
triangles, but the ground plane should be better represented by the latter unless
the square grid size is small.
A
The volume of excavation in triangle 124 will be i (d,+d,+d4) m in
which A is the plan area of the triangle and d,, d, and d, are the depths from
ground level to formation level at the corners. For triangle 245 the volume
A
of excavation is a (d,+d4+ds), and so on.
in which A is the plan area of a triangle, LD, is the sum of depths used once,
XD, is the sum of depths used twice, LD; is the sum of depths used three
times, etc.
Had squares (or rectangles) been used four depths would arise in each prism
Table 5.14
1 17.06-—x
2 52.44 —3x
3 35.26 — 2x
4 52.11—-—3x
5 106.20 — 6x
6 53.88 — 3x
7 52.74 —-—3x
8 108.06 — 6x
9 54.75 —3x
35.66 — 2x
54.57 — 3x
WNWAWWOWNW 18.42—x
|=
=
XD xX n = 641.15—36x
It will be noted that in Table 5.14 the multiplying factors 2, 3 and 6 have
been used directly, rather than adding the depths used twice, three times and
six times, respectively, and then multiplying by those factors.
10 x 10
Area A of each triangle = oe m?.
10 x 10
Therefore volume of excavation = jx 5 LDn.
The volumes for this type of problem can be calculated by the following
computer program. If the formation level is input as 10.0 m then the value
4686 m? is output as above. This particular problem would need a trial and
error solution, altering the formation level until the volume was zero. Lines
120 and 150 make use of the logic functions AND and OR; on computers that
do not support these functions line 120 must be replaced by four separate IF
... THEN GOTO 150 statements and line 150 by
5.10 Volume of a Excavation is required for the construction of a tank, the existing
truncated cone ground sloping at 1 in 5. The tank is to have a diameter of 12 m at
formation level, its centre being 3.50 m below ground level. If the side
slopes are at 2 vertical to 1 horizontal, determine the volume of excavation
to the nearest m?. *
Figure 5.16
--
Eb. Sap
isonet
Figure 5.17
Solution. First calculate the lengths of major and minor axes. From Fig.
5.17, in which the ground slope is 1 vertical in 5 horizontal,
Pepe oy 1 3.50—DD,
CF, +6.00 Sa 6.00+AD; |
But CF, = FF,/2 and AD, = DD,/2 since the side slopes are 2 vertical to
1 horizontal.
Nee 2.09
and D Fy = 12,00°+ ——'+ -—
D D.
=15:66m
Therefore DF = 2DG = 15.97 m
length of major axis.
G is the centre of the ellipse, and thus it is at a height of (FF; + DD,)/2 above
the formation level ABC, i.e. 3.66 m.
The length of the minor axis of the ellipse is equal to the length of the chord
of the circle, centre E,, described by the intersection of the side slopes with
the horizontal plane through G (Fig. 5.18).
78a
Figure 5.18
15.97
Length DG along the slope = 5 m
15.20 m.
15.
Sera eee
3 2 2
Hf ogee Gate
ae
537.72 m°>, say 538 m°,
5.11 Triangular Two stations were established on the top of a spoil heap. The following
ground surface model readings were recorded using a total station instrument sighting on to
points at the base of the heap. The height of instrument and height of
prism were equal for both sets of readings.
From A the values in Table 5.15 were obtained.
Table 5.15
Table 5.16
ie" 18°12"
181° 8’ 47”
236° 38’ 22”
313° 41’ 19”
52° 25’ 38”
What is the volume above a horizontal plane 5.8 m below station A?
Introduction. Figure 5.19 shows the spoil heap; the reader should note that
the top of the heap would be truncated if point 6 did not lie on the ridge joining
stations A and B. The volume can be calculated by establishing a ground surface
model consisting of a network of triangles. The volume under each triangle
is the plan area of the triangle multiplied by the mean height of the corners
above datum.
té
Figure 5.19
Solution. To calculate the plan area of the triangles we have the horizontal
length of two sides and the included angle, and
Area = ab sin C.
AREAS AND VOLUMES 173
For triangle A12: em
angle = ei” 21 wt eemeae 1 The
= 51° 10° 54”
area = 4 X 50.23 x 41.69 X sin 51° 10’ 54”
= 815.7914 m’.
Height of 1 = 5.8 — 5.24 = 0.56
Height of 2 = 5.8 — 4.93 = 0,87
Height of A = 5.8 = 5.8
7.23
+3 = 2.41. _
Volume = 815.7914 x 2.41 = 1966.06 m?.
The calculation is best computed as a table, noting that B is 0.38 m lower than
A when the level differences for Point 1 are considered.
Table 5.17
The following computer program can be used to solve this problem. The
DIM statement in line 20 is set for 20 readings at any one station. For computers
with limited memory this can be reduced, as appropriate; the value 6 is required
for this problem. The program is designed to work with ground and formation
levels. To solve this example input the formation level as zero, and the heights
as positive values above formation, i.e. level of station A = 5.8, B = 5.42,
point 1 = 0.56, etc.
Important: for any station, input should start with the point with the lowest
horizontal circle reading and proceed in increasing values of horizontal circle,
i.e. for A Point 4, 5, 6, ... 3; for B Point 8, 1, 6, ... 7. This allows the
instrument to be referenced on to an object or bearing that is not one of the
points, since line 270 calculates the oh a angle either side of the zero
reference line.
5.12 Calculating A digital terrain model is being established, the characteristic points
heights taking the form of a square grid of side L, shown in Fig. 5.20. Show
that the level, Zg, of Q may be expressed by
Z(L—x)\L—y) , Zx(L—Y)
Ze a L?2 L?
Figure 5.21
= Nga
Z,
L-x
+ A_—.
(|_|)
ie
In Fig. 5.20(d)
Zo = hy + (m — hh) tk
B
very) anual | be
L—x x y L-x Zox
ZA ir a3 Z +
(Tr) 47,
i L—x
pe er x
L- L- x
L
xy L—x\ y
PZ 78 PERE)
© 1? “Hae
ZL —x)(L—y) + Zp aL —y) te Zoxy
ibs EA [?
In triangle PQR, Fig. 5.21, with the origin for the grid co-ordinates at 0 we
get the results in Table 5.19.
From the results in Table 5.19 we get
PQ = V[(66.0 — 30.0)? + (93.0 — 44.0)”]
= 60.80 m,
PR = V[(93.0 — 30.0)? + (23.0 — 44.0)7]
= 66.41 m,
and QR = wv[(93.0 — 66.0)? + (23.0 — 93.0)7]
= 75.03 m.
X CO-ae 0 y co-ordinate
(m)
30.0 44.0
66.0 93.0
93.0 23.0
Now PQ + PR + QR = 2S = 202.24 m
and so S = l0i12 m >
S — PQ = 40.32 m
S — PR = 34.71 m
S — QR = 26.09 m.
Therefore area of triangle PQR = V(101.12 x 40.32 x 34.71 x 26.09)
1921.5 m?.
For the calculation of Zp, Zg and Zp we get the values in Table 3220:
Table 5.20
1921°5
X (3.94+3.44+1.9)
5892.6 m?
= 5893 m’, say.
r Figure 5.22
Referring to Fig. 5.22, which has been drawn using the data of the example,
it will be noted that:
(a) decreasing aggregate volumes, i.e. a to g, s to m, imply the formation
of an embankment;
(b) a minimum point occurs in the curve at the end of an embankment, e.g. r;
(c) when the curve rises cut is involved, i.e. from r towards b;
(d) a maximum ordinate will occur on the mass-haul curve at the end of
a cut, i.e. at s.
If a horizontal line is drawn, i.e. gh, the ordinates of the mass-haul curve
are equal at g and h, and so the volumes of cut and fill balance over that length.
When the curve lies below that trace line earth is moved to the left, i.e. h—r—
g, and, similarly, when the curve lies above, earth is moved to the right, i.e.
b—s—c. The length of such a line, termed a balancing line, indicates the maxi-
mum distance that earth will be transferred within the particular loop of the
diagram formed by that line. It will be seen that the base line gives continuous
balancing lines ab and bc, but continuity is not essential, the balancing lines
being arranged to ensure the most economical solution.
Haul is defined as the total of the products of increments of volumes of cut
and their distances of travel to the embankment, i.e. area grh gives the haul
in length gh. Two other terms are of importance: namely free haul and overhaul.
Their definitions are given later in the solution.
therefore to give equal borrow at chainage 2000 and 2600 the balancing line
must bisect that particular ordinate to give 1000 m* borrow at the two
chainages. This is satisfied by balancing lines passing through d, e and f, as
shown.
Finally we estimate the costs incurred. The rate for excavation generally
allows for transport up to a specified distance, known as the free-haul distance
or free haul. For transport beyond this distance we have overhaul and a different
unit rate applies. Free-haul distances can be plotted on the mass-haul diagram
as shown in Fig. 5.22, balancing lines gh and lm being 200 m long in this
case. The total volumes of excavation involved are given by the intercepts from
gh to r and from Im to s. Since cut balances fill over the length of a balancing
line earth would be carted a maximum distance of 200 m from h to g and |
to m, respectively.
The volumes remaining, gg, and In, in Scheme 1, and gj with In in Scheme
2, i.e. the ordinates from the respective balancing lines, are the volumes to
which the unit cost of overhaul will refer.
In Scheme 1, balancing line abc:
free-haul volumes (based on gh and Im) = (6800 —2100) +(4100 — 1700)
= 7100 m?,
overhaul volumes (intercepts gg, and In,) = 2100+1700
= 3800 m’,
borrow at chainage 2600 = 2000 m’.
Cost = (7100 x 1.20) + (3800 x 1.70) + (2000 x 2.40)
= £19 780.
In Scheme 2, balancing line def:
free-haul volumes (based on gh and Im) = (6800 —2100)+(4100— 1700)
= 7100 m’,
Widen 45)m'
2 .
173 250 x at
50
= £18 720.
volume of cut between the two sections by both the end-area and
prismoidal methods. [Bradford]
Answer 1951.6 m? (end areas); 1893.9 m? (prismoidal)
3 An embankment on a line due north and laid on level ground has
a uniform height of 16 m. The width at the top is 30 m and at the base
80 m. A road, formed at ground level, is to be cut through the
embankment with WCB of 60°. The width of the cutting at ground level
is to be 25 m, and the sides are to slope at a gradient of 1 vertical to
15 horizontal. Calculate the volume of material in m? to be removed
from the embankment in forming the cutting for the road. Check your
calculations by means of an accurately scaled plan and cross-sectional
drawings, and draw an isometric view of the cutting. [Eng. Council]
Answer 46 096 m?
E 168.83 N 146.88
E 306.09 S 14.10
E 266.71 S 38.26
E 177.74 S 105.56
| Line | AB | BC | CD DE | EA
AE 30 | —140 50 —20 | 80
—90 | - 20 60 | 70 | —20
Point | Es | N |
A 50.000 100.000
B 200.000 200.000
Cc 150.00 100.000
[Bradford]
Answer AP = 23.11 m; CQ = 14.33 m
7 (a) Give reasons why the estimation of areas and volumes is important
in most engineering schemes.
(b) Describe the planimeter and explain briefly how it is used to
measure areas from plans.
(c) State Simpson’s rule for the determination of areas and mention
the assumption which underlies the rule.
(d) In a survey of a field enclosed by a fence, offsets were taken to
the fence from a chain line as follows:
Chainage (m) 0 LO 20 30 40 50 60 70 80 90
Offset (m) 4.41 6.61 9.08 11.14 11.20 9.16 7.08 4.82 2.56 0
Use Simpson’s rule to determine the area between the fence and the chain
line. [Salford]
Answer 641.5 m?
14 The provisional centre line for a new road has been pegged out but
no further survey operations have yet taken place.
Table 5.26
— =
0 0
100 + 20.2
200 +36.2
300 +55.6
400 +61.7
500 +40.8
600 + 15.6
700 +) oO
800 —30.5
900 —58.7
1000 —72.3
1100 —45.6
1200 —30.0
1300 —22.5
1400 —13.0
1500 0
16 Estimate the cost arising in the previous question when the following
rates apply:
Excavate, cart and fill within a free haul distance of 500 m
£1.40/m?
Overhaul £0.70/m? x 100 m
Borrow and fill £2.50/m?
Answer £70 280 ’
L
Figure 5.23
(i) If the excavation has vertical sides, calculate the volume of material
to be removed. (Assume no bulking.)
(ii) It was subsequently discovered that the grid actually had been set
out at 20.1 m centres, and an error in the adopted temporary bench mark
had resulted in all the ground levels being shown 0.50 m too high.
Calculate the corrected volume of excavation. [Salford]
Answer (i) 17 820 m?; (ii) 16 989 m?
18 A hillside, which may be considered to be a plane surface, has a
slope of 1 vertical in 9 horizontal. A straight road is to be constructed
thereon at a gradient of 1 vertical in 12 horizontal such that it is wholly
in embankment. Determine the volume of fill required between two cross-
sections 50 m apart if a side slope of 1 vertical to 2 horizontal obtains
and the width of formation is 16 m.
Answer 551.6 m°
Table 5.27
Table 5.28
640 5.2
645 9.4
650 16.3
655 22.4
660 40.7
665 61.5
670 112.2
675 198.1
680 272.4
(a) If the 640 m level represents the level floor of the reservoir, use
the prismoidal formula to calculate the total volume of water impounded
when the water level reaches 680 m.
(b) Determine the level of water at which one-third of the total capacity
is stored in the reservoir.
(c) On checking the calculations, it was found that the original plan
from which the areas of contour loops had been measured had shrunk
evenly by approximately 1.2% of linear measurement. What was the
corrected total volume of water? [Salford]
Answer (a) 29.693 x 10° m?; (b) 669.8 m; (c) 30.406 x 10° m3
24 The plan of a very old chain survey, plotted to 1/500 scale on linen
cloth, was found to have shrunk so that a line originally 200 mm long
was only 197 mm. Furthermore, a note on the plan stated that the 30 m
chain used for the survey had been found to be 20 mm too long after
completion of the plot. If a certain area on the plan is measured by
planimeter as 0.225 m?, estimate the correct area on the ground in
hectares. Uniform shrinkage of the plan may be assumed. [Salford]
Answer 5.805 hectares
Table 5.29
Distance
from A (m)
Reduce the levelling by the Rise and Fall method of booking, applying
all arithmetical checks, and calculate the volume of the excavation in
cubic metres.
[Eng. Council]
Answer 1459 m?
Circular curves The circular curve, shown in Figs 6.1 and 6.2, connects two straights between
tangent points T, and T> at a constant radius R. The straights meet at the
intersection point I, crossing at angle 6 which is known either as the total angle
of deflection or the angle of deviation or the angle of intersection.
Intersection
point as
Tangent
point
Curve deflecting
right
T,P produced
Tangent at P
Figuré 6.2
190 SOLVING PROBLEMS IN SURVEYING
Formulae for the design and setting out of this curve include the following.
Tangent length: IT, = IT, = R tan (6/2).
Curve length: T,T, = R@.
Chord length: T,P = 2R sin a.
er:
Offset from tangent: y=
2R
Offset from long chord: , 2.
b = V(R*—a— ie 2_
[RK {1112\?
:
: TyP
Tangential angle: a = 1718.9 os (minutes), where
PP,
Offset from chord produced: P”P; = on (T,P + PP,).
Compound curves A compound curve has two or more curves contained between the two main
tangents or straights. The individual curves meet tangentially at their junction
points.
Reverse curves A reverse curve is one in which two curves change direction as shown in Fig.
6.3 so that I,I, is a common tangent to the curves.
Figure 6.3
Clothoid A clothoid is a curve whose radius decreases linearly from infinity to zero.
It is used in practice to ensure smooth progress between sections of different
radius, for example a straight and a circular curve or between two circular
curves. It will be apparent that some form of compound curve arises.
Tangent at T,’
1°
Radius = A
55 Radius = r
Radius = o
Figure 6.4
Transition curve
Figure 6.5
: : 1800 1? ;
Deflection angle for a specific chord: 6 = ————— (minutes).
tRL
1800 L
6, = (minutes):
aR
II
hy I?
_ + Bs
40 K? 3456 K4
13 2 4
x= = ae Se
14 440
13 7 jl
+ 2
6K 336 K? 42 240 K*
¢ is the angle between tangent to curve point and T,I as in Fig. 6.17.
3
¥ mal x 6 LR (cubic spiral).
piral)
3
% 6 LR (cubic Ppparabola).)
2
Shift: Sa - :
24R
Tangent length: IT, = (R+S) tan (6/2) + L/2
Rate of change of The centrifugal force P acting on a vehicle of weight W as it moves along
radial acceleration 2 P
a circular curve at velocity V is given by P = ey in which ~ is known
8g
as the centrifugal ratio. By lifting (super-elevating) the outer edge of the road
or rail the resultant force can be made to act perpendicularly to the running
surface, Fig. 6.6. In practice to avoid large super-elevations an allowance (fB)
y2
for friction is made. Radial acceleration is given by the expression Ro
(P cosa — Wsina) /
Figure 6.6
In the case of a transition curve it changes as the vehicle moves along the curve
because the radius is variable. For constant velocity V the rate of change of
radial acceleration (assumed uniform) is
Vertical curves Vertical curves are introduced at the intersection of two gradients, either as
summit curves or sag curves. Usually they are parabolic in form, and so in
Fig. 6.7
—(-q)|x 2
Vertical offset: y = Bl elles!)
400 L
- Gage:
400 L - ‘
Summit curve
Figure 6.7
In many published texts the term (p+q), or (p—q) in the case of gradients
of an equal sense, is replaced by one variable, G. For flat gradients it is normal
to assume that the length along the gradients equals the length of the curve
(2L) and that these in turn equal the long chord AC and its horizontal projection.
This assumption is not to be made in the design of horizontal curves.
Although the parabola is adopted for the vertical curve, it can be approxi-
mated by a circular curve in respect of radial (centrifugal) acceleration when
estimating the length of curve at the design stage. For the curve connecting
gradients of opposite sense, i.e. +p% and —q%
Curve length: 2L =
V*ptq) _ K (p+q).
100 f
Tables of K values are published by the Department of Transport for design
purposes, to satisfy different road conditions and velocities. Also for parabolic
summit curves:
+
960.89E
961.41N
G®
840.35E
940.22N 957.33E 1060.04E
eg 943.82N 924.52N
Figure 6.8
Solution
should preferably be of the order of R/20. Then the chord length is virtually
the arc length, and the sum of the individual chords is essentially the length
of the circular curve.
Notice (from Table 6.1) that curve points have been located in ‘running-
chainage’ form or ‘through-chainage’ form. This facilitates the levelling of
the longitudinal and cross sections over the complete length of road.
The following computer program will calculate data to set out a circular
curve from the tangent points. It is an extract from the complete program for
a compound curve listed with Example 6.5; the variables are defined with the
main program.
856.78 — 840.35
bearing FT, = tan! (
923.56 — ee
135° 23’ 54”
length FT; = V{[(856.78 — 840.35)? + (923.56 — 940.22)?]
= 23.40 m.
In respect of point P (chainage 1000.00) the bearing of T,P is
78° 36’ 30” + 05° 59’ 36” = 84° 36’ 06”.
From Fig. 6.9
T,P = 2R sin 05° 59’ 36”
= 73.09 m.
Note that the chord length is used, not the curve length. Co-ordinates of P are
easting 856.78 + 73.09 sin 84° 36’ 06” = 929.55 mE
northing 923.56 + 73.09 cos 84° 36’ 06” = 930.44 mN.
Bearing GP Il
_, 929.55 — 960.89
930.44 — 961.41
= 225°.20' 25”.
Figure 6.9
A clockwise angle of (225° 20’ 25” — 110° 24’ 30”) = 114° 55’ 55” could
be set off at G from GH to give the required direction. In practice, the surveyor
would probably sight H and set the horizontal circle reading on the theodolite
to 110° 24” 30”, acircle reading of 225° 20’ 25” then gives the pointing on P.
Length GP = V[(929.55 — 960.89)? + (930.44 — 961.41)7]
= 44.06 m.
1020.123 742.194
1135.706 851.208
1297.041 863.197
1415.313 797.071
1198.827 794.213
| Straight 2
Straight 1
Figure 6.10
m, Ey + cy = m Ey +
Thus E; = a ts
m,—m)
Now,
a 1588.373 —(—219.986)
0.9432 —(—0.5591)
= 1203.727 mE
and N = mE; i Cy
= 0.9432 x 1203.727—219.986
= 915.369 mN
The deflection angle of the curve,
6 = tan~!(m,)—tan~
!(m))
tan~ !(0.9432)
—tan~ '(—0.5591)
72° 32403"
Note that if the curve had been in a north/south orientation @ would have been
obtained from the inverse of m, and mp.
6 can be positive or negative and in conjunction with the orientation of the
curve this indicates whether the curve is left-handed or right-handed.
The length to the tangent point, T, is given by
P=" tan oe
2;
Se5tM) Sn Vaewao OS
wv
= 146.738 m
23
— 1203.727)
Ey = 1203.727 + uN eLe x (1020.1
252.389
$2. f2...03”
1118.776 + 200 x a ———— XK 7
180
1371.968 m
For running chainages at 20 m intervals the first curve point will be at chainage
1120.000. The curve length to this point is (1120.000—1118.776) = 1.224 m.
1.224
Angle subtended at O = = 0.00612 radians
200
= 21/2”
WCBo; = 316° 40’ 28” + 21’ 2”
=S17% 1" "50"
Ey Eo + R sin WCBo;
= ea. a0
The logic function AND is also used in lines 390 and 630 and if necessary
this should be replaced by:
390 IF E9>E7 THEN GOTO 391
391 IF N9>N7 THEN GOTO 400
392 A4 = A4 4+2*P
with a similar routine at 630.
Variables
A WCB from control station to ES, NS Co-ordinates of 2nd point on
a point straight 2
Al Deflection angle E6, N6 Co-ordinates of intersection
A4 WCB from centre of circle to point
a point E7, N7 Co-ordinates of first tangent
C Running chainage point
Cl Chainage of first tangent E8, N8 Co-ordinates of second
point tangent point
C2 Chainage of second tangent E9, N9 Co-ordinates of centre of
point circle
D Output angle, degrees Gl ‘Gradient’ of straight 1
El, N1 Co-ordinates of control point G2 ‘Gradient’ of straight 2
E2, N2 Co-ordinates of 1st point on G3 ‘Gradient’ from first tangent
straight 1 point to centre of circle
E3, N3 Co-ordinates of 2nd point on G4 ‘Gradient’ from second
straight 1 tangent point to centre of
E4, N4 Co-ordinates of 1st point on circle
straight 2 of Chainage interval
6.3 Reverse curve Two parallel railway lines on the surface at a mine are to be connected
by a reverse curve, each section having the same radius. If the centre
lines are 50 m apart and the maximum distance between the tangent points
is 200 m calculate the maximum allowable radius that can be used.
[CEI]
Common tangent
at U
Vi
Figure 6.11
Solution.
0 1
Therefore sin — = ——
oy Are
4 0
and so cos — = — and tan— = —
2 aia? >
But 2R sin 6/2 = 25V17 i
6.4 Setting out a Calculate the data needed to set out one of the curves of Example 6.3
circular curve from using offsets from the long chord.
the long chord
Figure 6.12
Table 6.4
212.13
Tangent
at A,
Circular
ea | curve
=
\ R= 1000m R= 1000m
\
Figure 6.13
Solution. First calculate the curve lengths and the tangent lengths. The rate
of change of radial acceleration (a) is related to the design velocity (V) by
v3
Of ees
LR .
where L is the length of the transition curve and R is the radius of the circular
curve (and the transition curve) at A.
Shift ee L2 = 95.097 2
24R 24 x 1000
= 0.377 m.
Tangent length IA = (R+S) tan 6/2 + L/2
= 1000 x —— x 26.9511
180
= 470.385 m.
(a) We can now determine the chainages of the tangent points:
Chainage of intersection point = 1350.468 m
Deduct tangent length = 338.185 m
(c) Determine data to set out the curve with theodolite and tape. For the
transition curve the deflection angle is
_ 1800 1?
= 0.006 025 J” min.
a LR
For the transition curve from A the data can be tabulated as in Table 6.6.
Table 6.6
Table 6.8
in question is also output so that setting out can alternatively be carried out
with an EDM.
Variables
A Angle for DMS subroutine I Setting out chainage interval
Al Total angle of deflection of the Ll Length of the transition curve
curve, 0 L2 Length of the circular curve
A2 Deflection of the transition curve, L3 Running curve length on transition:
> chord length on circular curve
A3 Deflection of the circular curve, Setting out polar length
9-29, Input/output angle, minutes
A4 Setting out angle, 6 Radius of the circular curve
(6 Chainage of the intersection point Input/output angle, seconds
Cl Chainage of the first tangent point Shift
C2 Chainage of the second tangent Tangent length (IT)
point Design speed
GS Chainage of the third tangent point x<4nnwer
Offset of the transition from the
C4 Chainage of the final tangent point tangent
C5 Running chainage Offset at the tangent point
D Input/output angle, degrees see Offset of the transition along the
F Variable with value zero for output tangent
a Offset at the tangent point
9,
Figure 6.14
Solution. Determine the existing tangent lengths. For the first circular curve
R, 6; = 1192.95 — 740.40 m
452.55 m
Since R; = 540m
6, = 0.838 055 6 rad
= 48° 01’ 01”
Therefore tangent length AA; = R, tan @);/2 = 240.52 m.
For the second circular curve
Next determine the distance from the first tangent point to the intersection point.
The length of the common tangent A,C; at the junction point of the two
existing circular curves
ll240.52 + 87.56
II328.08 m.
Also 6 = 0, + 0 = 70° 02’ 16”.
By the sine rule
1207
24 x 660
0.91 m.
|7oron48"
Therefore new tangent length = (660+0.91) tan + 120
523.10 m.
Therefore chainage of new tangent point = 1111.79 — 523.10
588.69 m.
Now ¢, =
Bog SP = 0.090 909 1 rad
2LR 2 x 120 x 660
mS oS deol.
59.6206°.
= 686.78 m.
Chainage of second tangent point Il 588.69 + 120.00 + 686.78
+ 120.00
1515.47 m.
Now calculate offsets to locate the quarter points on the transition curve.
The equation of the offsets to the transition curve is
13
~~ 6LR ‘
212 SOLVING PROBLEMS IN SURVEYING
6.7 Locating a A curve, wholly transitional, is to be set out in a built-up area. At
transition curve from the preliminary survey and design stage it is found that the intersection
a traverse point I is inaccessible. A and B, E and F lie on the two straights and
they are linked with a traverse survey as shown on Fig. 6.15. The data
in Table 6.10 were obtained.
Table 6.10
Horizontal Horizontal
distance clockwise angle
(m)
196° 11’ 40”
189° 45’ 20”
204° 23’ 20”
193° 19’ 40”
Assuming the design value for the rate of change of radial acceleration
to be 1/3 m/s? for a velocity of 50 km/h determine the positions of the
tangent points.
Figure 6.15
Table 6.11
94.21 m 151.96 m
For the transition curves. The curve is to be wholly transitional and therefore,
by symmetry, for each of the two transition curves
When the curve turns through an angle of this magnitude the first-order
equations used between ¢, and (63° 40’ 00”)/2 in previous examples are not
valid. In Fig. 6.17 we can write dx = dl sin ¢ and dy = dl cos ¢.
3 5 2 4
dx =-(¢-$-+4)a eer:
3! 5! !
214 SOLVING PROBLEMS IN SURVEYING
Figure 6.17
]2 1° j}0 14 18
= (—— — ah 5} dl Os, fcTyee, sees oh enor,
2K 48K 3840K 8K 384K
Integrate
13 ]7 pl
=— — + —————
6K 336K? 42240K°
and
E i}?
= — + ————.,
40K2 3456K4
There are no constants of integration since ¢ = 0 when / = O. If L is the
length of each transition curve and the minimum radius at junction point T
is R, then
o, = 0.555 5965 =
50 ma
y3 3600
also LR. a me Sn dee
Therefore LR = 8037.55 = K.
Since Tae 1 LI 193k
R = 85.05 m
and L = 94.50 m.
Common tangent
at S
Circular
curve
Transition curve
Circular
curve
Figure 6.18
Solution. First establish the length of transition curve ST and the angle
contained between the tangent at P and the common tangent at S. From the
properties of the cubic spiral we have
Ir = LR =K =1,R,
where L is the total length PT, and /, is the length PS, the radii there being
550 m and 950 m, respectively.
The rate of change of radial acceleration
V3
[0AN i Spear
K
Now100 km/h = 27.78 m/s,
AR Ee
thereh L _ Ley
ee yy O8p. AT 550
Therefore 1, = 67.70 m and L = 116.94 m.
Therefore ST = 116.94 — 67.70 = 49.24 m.
Also LR = K = 64 315.84.
l,? :
Therefore @ = —’— = 0.035 631 6 rad
Table 6.12
. @-1)
6K
Whence we get the values in Table 6.13.
Table 6.13
Station Chainage
(m)
2474.00
2480.00
2490.00
2500.00
2510.00
2520.00
2523.24
6.9 Replacing a A circular curve of radius 1600 m connects two straights having a total
circular curve with a deflection angle of 65° 30’. It is to be shifted to allow a cubic spiral
compound curve transition curve of length 125 m to be inserted at each end, the total route
length remaining unchanged.
Calculate the distance between the new and the previous tangent points.
Existing
tangent
point
New
tangent
point
R = 1600m
Existing curve
Figure 6.19
Centre of
Figure 6.20
oy) = ——
: radian,
di
pets
where R is the radius of the new circular curve.
L?2
Shitts: = ———,
24R
where L = 125.00 m. Therefore new tangent length
Ve
IT,’ = ze + (R+S) tan 32° 45’
5 ce
= — <2 |R + ——) tan 32° 45°.
2 24R
The total curve length T,’T,’ has to equal the length of the existing circular
curve plus 2T,’T}, since the total route length must remain unchanged. Thus
L
1829.11 + 2T,T,’ =2L + (1.143 1907 = ae
L £2
But “T|T,” = me + (x+ Smal tan 32° 45’ — 1029.15
125.007
and so 1829.11 + 125.00 + 2(R+ | X 0.643 221 6
La 125.007
Now shift S = = 0.41 m.
24R = 24 x -1596.24 m
Therefore (R+S) 1596.24 + 0.41 = 1596.65 m.
125.00
IT,’ = bate + 1596.65 tan 32° 45’
Chainage
(m)
4180.00 0’ 04”
4200.00 1’ 47”
4220.00 5’ 49”
4240.00 12’ 08”
4260.00 20’ 45”
4280.00 31° 41”
4300.00 44’ 52”
1089.50 m.
Therefore T,T,’distance between the new and existing tangent points
1089.50 m — 1029.15 m
60.35 m.
Chainage of T,’ 5264.50 — 1089.50
= 4175.00 m.
(a) To set out transition curve T,’T,” deflection angle (6) can be calculated
for chord lengths (/) using the expression
1800/7,
= min
a RL
0.002 871 5 77 min in this case.
(b) Calculate data to locate X, the centre point of the curve (Fig. 6.20).
i (Karo ) = : + 4 0.41
1596.24 = 1898.43m
cos 32° 45’ cos 32° 45’
6.10 Compound A new road with a 100 km/h design speed has two straights joined
curve set out from by a curve consisting of two transitions and a 500 m radius circular curve.
control points The rate of change of radial acceleration on the transitions is limited
to 0.3 m/s?. The co-ordinates in Table 6.15 have been fixed.
Calculate the setting out angles at A and B,related to line AB, for
the mid-point of the curve. [Salford]
Introduction. Most major roads are now set out from control points using
either two theodolites and intersecting rays or, more usually, an EDM with
a bearing and distance as indicated in Example 6.1. The computer program
at the end of this problem allows suitable data to be calculated for the case
of a compound curve.
27.78?
sO Li=
0.3 x 500
142.890 m.
LE 142.89
$)
2R 2.x 500
II 0.142 89 radians = 8° 11’ 13”.
40(LR)?
3
ee a
40 x 500?
= 142.598 m.
Calculate the co-ordinates of tangent point T,. For the transition curve at T,
6.7
6; = tan! eee
142.598
= 2° 43’ 43”.
and length TT, = V(6.7967 + 142.5987)
= 142.760.
From the co-ordinate data
1553°203' = 1298.179
WCB of TI = tan!
770.889 — 719.840
78° 40’ 50”.
WCB of TT, 78° 40’ 50” + 2° 43’ 43”
Bios ao. 5
E co-ordinate of T, 1298.179 + 142.760 sin 81° 24’ 33”
1439.338 m.
N co-ordinate of T, 719.840 + 142.760 cos 81° 24’ 33”
= 741.165 m.
Next find the deflection angle of the curve from the co-ordinate data.
meet 770.889 — 719.840
Deflection angle of curve
1553.203 — 1298.179
770.889 — 641.230
+ tan~!
1778.662 — 1553.203
II PLP 19S 10 Se 28 54 1G
= 41° 13’ 20”.
41° 13720"
Deflection angle of mid-point, M = 20° 36’ 40”.
2
Deflection angle of circular curve to M 9d ad | sh: Uae teeect IS I
weal IAN ee Wie Odi
le SR
WCB of T,M (Fig. 6.21) = 81° 24’ 33” + $,—6,+
2
= $1524 soa, ern il? 13” —2° 43’ 43" + 6° 12° 43”
i.e. WCB of T;|M = 93° 04’ 46”.
12°25. 27
2
Figure 6.21
Male 22) 27
= 29s A) Sit
= 108.209 m
Co-ordinates of M:
By cosine rule
AB? + AM2 — BM?
cos (
2 x AB x AM
)
164.9447 + 168.3627 — 98.0152
Il cos~! (
2 X 164.944 x 168.362
34° 11’ 02”
)
164.9442 + 98.0152 — 168.362
ost (
2 xX 164.944 x 98.015
74° 49’ 06”
The following computer program will calculate the setting-out data for a
compound curve located from two control stations. Co-ordinates of the control
stations, the intersection point of the straights and the tangent points with the
straights must be known. The data is output in a table and readers with
Variables
Al Deflection of the whole curve, 0 L6, L7 Length from control points to
A2 Deflection of the transition curves, points being calculated
> L8 Length of a chord to the transition
A3 Deflection of the circular curve, curve
6—2¢, Output angle, minutes
A4 WCB of first tangent/final tangent qT
AS WCB of chord to transition curve QHve
‘S’ in the triangle formed by
A6 Setting-out angle for the circular control points 1 and 2 and the
curve point being calculated
A7 WCB of chord to the circular 2 Constant term in the expression for
curve the angle
A8 , A9 Component of Al Radius of the circular curve
GC Chainage of the intersection point Output angle, seconds
Cl Chainage of the first tangent point —HnDX offset of the transition from the
C2 Chainage of the second tangent tangent
point Y offset of the transition along the
Chainage of the third tangent point tangent
Chainage of the final tangent point Road design speed
C6 Chainage of the current Sin~' of the angle at the control
position point
Output angle, degrees Xx X co-ordinate of current point
+ counter X1 X co-ordinate of control point 1
Setting out chainage interval X2 X co-ordinate of control point 2
L1 * R for the transition X3 X co-ordinate of intersection point
Length of transition x4 X co-ordinate of tangent point with
Length of circular curve first straight
Running curve length X5 X co-ordinate of tangent point with
Length of a chord to the circular exit straight
curve ye Y co-ordinate of current point
L5 Distance between control points Y1—Y5 Y co-ordinate as above
0, Otp, 10, 15, 18, 20, 23, 29, 37, 39, 46, 58tp, 62, 60, 61, 62, 62,
63tp, 59, 46, 36, 32, 21, 17, 15, 10, 6, tp, —1, 0,
where tp indicates a tangent point on the curve. Calculate the necessary
slew of the track using the moment method.
é chord
points
Versine at 4 amended
by half the slew at 5
Figure 6.22
where L is the transition length and / is the running transition length. Note
that the versine is a constant for the circular curve and for the cubic parabola
transition curve it increases at a uniform rate.
Calculate the versines for the new curve. It is standard railway practice to
ease the entry and exit from the transition by applying a versine at the tangent
to the straight equal to 1/6 of the increase and by reducing the value at the
tangent to the circular curve by the same amount.
To correct the data by the method of moments a tabular format is used where:
Difference Required — Actual versine;
L difference Sum of differences;
Moment Sum of 2 differences;
Slew —2 X moment (sign convention outwards +ve,
inwards —ve). !
1
2
3
4
tS)
6
7
8
9
Note that from Fig. 6.22 a slew movement at point 5 amends the versine at
points 4 and 6 by half of the slew, and the method of moments carries this
correction cumulatively around the curve.
The 2 mm slew required for the exit straight would in practice be run out
evenly at the end of the transition and the start of the straight.
6.12 Sag vertical Two straights AB and BC falling to the right at 1 in 10 and 1 in 20,
curve respectively, are to be connected by a parabolic vertical curve 200 m
long. Given that the chainage and reduced level of B are 3627.00 m and
84.64 m, respectively, design the vertical curve.
What is the sighting distance for a car whose headlights are 0.70 m
above road level and the beams are inclined upwards at an angle of 1°?
ee point
A
A1
P =-10%
Figure 6.23
The general expression for the offsets (y) from a tangent to locate the curve is
(p—q) x’
400 1
where x is the distance from a tangent point; and / is the length of a tangent
(half the length of the curve). Therefore for this vertical curve
_ (1-9)
400 x 100
5,
ip
400 x 100’
when we ignore the negative sign. It is usual to take the curve length to be
equal to the total length of the tangents and, for gradients of this order, to
assume that they are equal to the horizontal lengths. When x = 100 m
The offsets can be evaluated from both tangents or from either tangent,
produced as necessary. The second alternative will be adopted in this example
and tabulated in Table 6.17 using running chainages at 20 m intervals. Offsets
must then be added to the corresponding grade (tangent) levels.
With the car at tangent point A,, the headlight beams will strike the road
surface at a point where the offset is (0.70 + x tan 1°), x being the distance
from A,. Thus
5
0.70 +.0.0175x = ————— x?
400 x 100
and x = sighting distance in this case
= 172.47 m.
6.13 Summit vertical Two straights PB and BQ are to be connected at a summit by a parabolic
curve vertical curve. P, reduced level 108.25 m, lies on the gradient rising
to the right at 1 in 60 at chainage 1862.00. Q, reduced level 106.85 m,
lies on the gradient falling to the right at 1 in 50 at chainage 2174.00.
The vertical curve must pass through R, reduced level 109.68 m and
chainage 1986.00.
Design the curve and determine the sighting distance between two
points 1.05 m above road level.
Solution. First determine the length of the curve. Let horizontal distance
PB = x, in Fig. 6.24.
=l10. 32m:
Curve level at R = 109.68 m.
Therefore offset at R = 0.64 m.
Sidney if ay.
4001
in which p = 42% and q = “4%.
Therefore / = 84.72 m, say 85 m for design purposes, and so the chainage
of T is 1994.00 — 85.0 = 1909.00 m.
Table 6.18
(m)
1909.00
1920.00
1940.00
1960.00
1980.00
1994.00
2000.00
2020.00
2040.00
2060.00
2079.00
Pre ew?!
100 4001
’
ae nie a
100 4001
ae 2p |
max (p+q)
2 X 10/6 x 85
= = 77.27 m.
10/6 + 2
At this point the grade level is 110.32 m AOD
(10/6 + 2)
and the offset = x (77.27)? = 0.64 m.
400 x 85
Therefore the highest point has a reduced level of 109.68 m.
Finally calculate the sight distance.
The offset AB at B = 0.78 m
= ING?
The sight line can be taken as the tangent to the vertical curve at A (Fig.
pl ql
100 100
6.24). It is parallel to TU whose slope is rad.
21
The angle contained between the tangent TB and chord TU at T is
|
Toe AB ee es Ria ee
+
200
s; = 14.73 m.
Therefore total sight distance = 2 (85.00 + 14.73)
= 200 m, say.
Note, as mentioned previously, curve lengths can be established in practice
using an expression of the type:
L = 21 = K (p+q) = KG,
in which K depends upon the value of centrifugal acceleration chosen for a
60 50
Therefore K = 46.4.
Problems .
2 (a) Describe a method for setting out a circular curve of small radius,
such as a minor road kerb line or a boundary wall, using only a chain
and tape. You should derive any necessary formulae from first principles.
(b) Tabulate the data required to set out, by chain and 20” theodolite,
a circular curve of radius 500 m to connect two straights having an
intersection angle 15° 20’ 00”. You may assume that the chainage of
the intersection point is 2100.00 m and that the curve is to be defined
in 20 m intervals. [Salford]
Answer Chainages of tangent points 2032.69 m, 2166.50 m
6 (a) Name the purposes which a transition curve serves in road design
and discuss briefly the factors which influence the length selected for it.
Establish the formula
B E 500.00 N 550.00
CE 1300.00 N 700.00
Calculate the radius of the circular curve passing through the three stations
and the co-ordinates of the centre of curvature.
Describe briefly how such a circular curve could be set out below
ground. [CEI]
Answer 1223.67 m; 1112.67 mE, 509.25 mS
11 Two straights AI and IB have bearings of 37° 15’ 00” and 48°
19’ 40”, respectively, and are to be connected by means of a wholly
transitional curve with a rate of change of radial acceleration of
0.3 m/s*. If the design speed of the road is 100 km/h and the co-
ordinates of point I are 2664,83 m E, 2742.88 m N calculate
(a) the reduced levels and chainages of the tangent points and the lowest
point on the curve,
(b) the reduced levels of the first two pegs on the curve, the pegs being
set at the 20 m points of the through chainage. [ICE]
Answer (a) 2605.92 m, 281.13 m; 2651.39 m, 264.90 m; (b)
280.54 m, 280.07 m
22 The bearings of two straights AI and IB are 15° 26’ 12” and 36°
19’ 20” respectively. The co-ordinates of the intersection point I are
2440.43 mE, 2336.94 mN. The two straights are to be connected by
a wholly transitional curve with a rate of change of radial acceleration
of 0.15 m/s?. If the design speed of the road is 120 km/h, calculate the
minimum radius of the curve, the co-ordinates of the two tangent points
and the co-ordinates of the point on the curve 100 m from the tangent
point on AI. [Leeds]
Answer R = 823.021 m; co-ordinates of T, 2359.95 mE,
2045.48 mN; co-ordinates of T, 2619.53 mE, 2580.55 mN
Table 6.19
Point mE mN
1031172 610.129
924.031 809.142
956.138 1025.836
1129.776 1221.781
Setting out The procedure adopted to ensure that a specific design feature, i.e. a building,
a road, a pipeline, etc., is correctly positioned at the construction stage is known
as setting out. Control is provided by stations whose positions were fixed during
the original survey or by subsidiary stations located by ‘intersection’, or by
‘resection’ from such stations. Setting out is thus the reverse of detail surveying
in that the control stations are used to place points on the ground in their relative
positions. If the setting out is referred to a control system based on National
Grid co-ordinates the local scale factor should be considered.
Control points If two control points A and B are known, a third point can be located in a
number of ways (see Fig. 7.1).
B B B B
f)
ee C 0
Cc Cc
A A A A
(a) (6) (c) (d)
Figure 7.1
Intersection Figure 7.1(c) shows the method of intersection by which C can either be co-
ordinated by observing 6, and 6, or be located if its co-ordinates are known,
since 6, and 6, can be calculated.
Resection By the method of resection a station or point can be located after pointings
have been made on at least three known stations. This technique is very useful
7
Tunnelling — the Generally, plumb wires are used to transfer directions underground. Essentially,
transfer of alignment the plumb wires produce a vertical reference plane and on the surface the plane
below ground can be placed in the line of sight of a correctly oriented theodolite; below
ground, the line of sight can be directed into that plane. This is known as ‘co-
planing’ and the line of sight when established can be used to set up floor or
roof stations within the tunnel.
Weisbach triangle In an alternative approach the theodolite is positioned out of the vertical plane
so that, in plan, each wire and the vertical axis of the theodolite form a triangle
— the Weisbach triangle. The small angle subtended at the theodolite is
measured and its value is included with other measurements to determine either
the orientation of the vertical plane at ground level or the fixing of reference
points below ground. Instrument positions can be changed (also when co-
planing) to improve accuracy, and in deep shafts more than one pair of plumb
wires can be put down.
7.1 Setting out by The four corners A, B, C, D of a rectangular building are to be set out
bearing and distance from control station P by a total station instrument. Using the co-ordinates
in Table 7.1, what whole circle bearing and distance is required to
establish each of the corners?
Table 7.1
102.731
67.376
17.879
53.234
WOBpr=xtan.- (fa)
Na-N,p
- (102,731=95.150
= tan ——
81.707
—70.192
a (=e ) “i
= tan
11.515
= 32° 25’ 43”
PA. = V((E,—E,)” + (Na-N,)”)
= V(7.601? + 11.515’)
= 13.797 m
For B,
E,—E
WCBpz = tan7! (2-2)
Ng—N,
Se h(i SR
117.062 — 70.192
et (<2)
= tan ——
46.870
This is not the WCB of the angle because a WCB cannot be negative. A
calculator presents angles from trigonometric functions relative to the north—
south axis with positive and negative values dependent upon the quadrant
containing the angle. For the case of tangent values the first and third quadrants
are positive and second and fourth quadrants negative. In this case the northing
of B is greater than the northing of P so B is in the second quadrant relative to P.
B
>N
A
, —30° 37’ 54”
D
i
Figure 7.2
244 SOLVING PROBLEMS IN SURVEYING
From Fig. 7.2,
WCBp3 = 360°—30° 37’ 54”
= 329° 22’ 06”
PB = V((—27.754)? + 46.8702)
= 54.471 m
For C,
This is not the WCB of PC because both the easting and northing of C are
less than those of P as indicated by the negative signs on the top and bottom
of the calculation for tan~' (compare this to the two positive signs obtained
for corner A). The angle is thus positive but in the third quadrant.
B
D Figure 7.3
For D,
En—E
WCBpp PD = tan! (Fo)
Np—N,
es (B26a 0)
32.210 —70.192
abe (cis)
— 37.982
47° 48’ 18” (on a calculator)
Variables
A WCB in radians E2, N2 Co-ordinates of the point
A$ Control variable L Output length
D Output angle, degrees m Output angle, minutes
El, Nl Co-ordinates of the control p T
station S Output angle, seconds
7.2 Intersection A new control station C is to be located from two stations A and B,
the horizontal clockwise angles recorded being BAC = 49° 27’ 18”
and ABC = 322° 45’ 24”. Determine the co-ordinates of C given that
those of A and B are:
A 950.00 mE 1200.00 mN
B 983.50 mE 1340.00 mN
AD (FE, —Ea)
teers
Ne me (Ne
Ne ea) 0|A DD) ee, on
Variables
Al Angle from station A E1, N1 Co-ordinates of station A
A2 Angle from station B E2, N2 Co-ordinates of station B
Cl, C2 Cot of Al, Cot of A2 E3, N3 Co-ordinates
of the point
D Input angle, degrees M Input angle, minutes
S Input angle, seconds
10 REM INTERSECTION
20 INPUT"INPUT OBSERVATION ANGLE FROM STATION A IN DEG,MIN,SEC ",D,M,S
30 Al=((D*3600)+(M*60) +S) /206264.8
40 INPUT"INPUT OBSERVATION ANGLE FROM:STATION B IN DEG,MIN,SEC ",D,M,S
50 A2=((D*3600)+(M*60) +S) /206264.8
60 INPUT"INPUT CO-ORDINATES OF A ",E1,N1 -
70 INPUT"INPUT CO-ORDINATES OF B ",E2,N2
80 C1=ABS (1/TAN(A1) )
90 C2=ABS(1/TAN(A2) )
100 E3=INT((((E1*C2)+(E2*C1)+N2-N1) /(C1+C2) ) *100+0.5)/100
110 N3=INT((((N1*C2)+(N2*C1) +E1-E2) /(C1+C2) ) *100+0.5) /100
120 PRINT"COORDINATES ARE ",E3,N3
130 END
7.3 Setting out by The co-ordinates of two survey stations F and G are:
intersecting rays
Station F 3812.07 mE 1631.32 mN
Station G. 3669.35 mE 1746.89 mN
Figure 7.5:
From F to G
Table 7.3
From F to H
Easting difference = —112.07 m
Northing difference = +43.68 m
— 112.07
Therefore bearing of FH = tan! Gacane 2OVGtr 31.
Table 7.4
1120.351 cot 48° 31° 18" +1271.782. cot 35° 11° 23%+
861.592— 759.787
coLds” 31 4418* cot 35° 11° 23”
1251.477 mE
- Na cot o2+Npg cot @, —EgtE,
159.787 cot 48° 31” 18* +861.592 cot 55° 11’ 23” —
1271.782 + 1120.351
cords 311s" Fcor gs” 11’ 23”
708.722 mN
Next calculate the vertical position of C.
From A,
Hpc = V((1271.782
— 1251.477)2 + (861.592 —708.722)*)
=9154.213 m
Zc = Hg + height of instrument + Hgc tan 6)
= 131.152 + 0.403 + 154.213 tan 4° 54’ 42”
= 144.8074 m
Hence,
= V(8? + 127)
= 14mm
7.5 Setting out Two stations on a traverse conducted in an urban area have co-ordinates
using the of 950.00 mE, 1200.00 mN and 983.50 mE, 1340.00 mN, respectively.
perpendicular A manhole is to be located at a point whose co-ordinates are 975.00 mE,
1266.00 mN. It is known that neither direct sightings nor measurements
can be made from the two stations and accordingly the engineer decides
to locate the manhole by linear measurements along and perpendicular
to the traverse line. Determine those measurements.
AD
but cot om = —
CD
and cot dD = —
poy, in |SpFig.
CD te
g a7.2.
AD (AB— AD)
And so —cp (Ng-—Nc)
“8 c) + ———Cp (N,—-N,-)
(Na—Nc) == (Ern—E
(Ep—Ea)
AD CD
whence ——
Ge (Ng—Na)—
(Np — Na) —
; (Fp—E,)
(Eg — Ea) = No—N,.
fom is’
wy AD CD
Similarly a (Ep — E,)+ Fie (Ng—Na) = Ec—Esg.
Thus AD and CD can be estimated provided that the three points have been
co-ordinated. i
AB = V[(Eg—E,)?+(Np—Na)?]
143.95 m.
Substituting in the two expressions
140.00 CD :
33.50
= = 66.00
143.95 143.95
and so AD = 70.01 m
CD = 8.95 m.
By setting out a length of 70.01 m from A along AB, and then setting off a
perpendicular of length 8.95 m from that position the manhole is located.
An alternative approach is as follows:
E,—E
Bearing of AB = tan~! —B——4-
=113°.20 25.3.7.
Ec—E
Bearing of AC = tan~! —S—4
Nc—Na
= 20°44" 45.9".
7.6 Resection In order to set out the centre line of a bridge it is required to establish
the point D which lies 115.00 m due south of the peg A, across the river.
Stations B and C were established and referred to A by subtense
tacheometry. Field notes were as in Table 7.6.
Target bearings
3 ses Station
Sg ei
Possible ae eae
locations of
Pp
Figure 7.6
BC BP BA BP
—
sin 6,
= —
sin a
and —
sin 05
= ——.
sin B
BC sina _ BA sin B
Therefore BP = : = ;
sin 6; sin 0,
: BA sin 6, ;
si o = —-————— Xam6.
BC sin 0,
: BA ;
Therefore - ae x - = k, say,
sin 6 BC sin 65
but at+pB 360° — CBA — 6, — 0;
o, which is known.
BA sin 0;
Therefore sina ar oe m4 Se
X sin (6—a).
BC sin 65
This expression allows the determination of a, and hence 8. There is now
sufficient information to calculate lengths CP, BP and AP, allowing P to be
fixed in position.
Calculation of lengths CB, BA and angle CBA. The lengths have been
determined by subtense measurements. Assuming that a bar length 2 m was
used (Fig. 2.4 refers)
b 2
AB — =
Cc
~N
sS 208
B 60.306 4
Figure 7.8
AC CB
Also ———— :
sin CBA sin 14° 30’ 00”
In triangle BAP
7.7 Resection — Point P was established within the triangle formed by control stations
Tienstra’s method A (1020 mE, 2560 mN), B (1360 mE, 2520 mN), C (1150 mE,
2070 mN), such that angles APB = 128° 20’ 20” and BPC =
135° 47’ 40”. What is the distance of point X (1200 mE, 2450 mN)
from P and what is the WCB of XP? [Salford]
_,[506.9527 + 342.345? — a
‘O
m 2 x 506.952 x 342.345
68° 25’ 54”
SETTING OUT AND POINT LOCATION 259
: _, (342.345? + 496.5887 — 506.952?
CBA = cos
2 x 342.345 x 496:588
= 71°4I34"
‘ 496.5882 + 506.9527 — 342.3457
ACB = cos~! os
2 x 496.588 x 506.952
BOeso2 32
Check: 68°.257-54" + 71° 41° 34” + 39° 52’ 32” = 180° 00° 00”.
1
K, — A A
1
eS - = 0.702 537
(cot 68° 25’ 54” — cot 135° 47’ 40”)
1
a = 2.306 219
(cot 71° 41’ 34” — cot 95° 52’ 00”)
i 1
> (cot ACB — cot APB)
1
a = 0.503 049
(cot 39° 52’ 32” — cot 128° 20’ 20”)
Whence K, + Ky + K3 = 3.511 81
Therefore
me (0.702 537
x 1020) + (2.306 219x 1360) + (0.503 049x 1150)
3 3.511 81
= 1261.900 m 4
10 REM RESECTION
20 INPUT"INPUT ANGLE BETWEEN STATIONS A AND B ",D,M,S
30 Al=((D*3600)+(M*60) +S) /206264.8
40 C1=1/TAN(A1)
50 INPUT"INPUT ANGLE BETWEEN STATIONS B AND C ",D,M,S
60 A2=((D*3600)+(M*60)+S) /206264.8
70 C2=1/TAN(A2)
80 C3=1/TAN(3.14159-A1-A2)
90 INPUT"INPUT CO-ORDINATES OF A ",E1,N1
100 INPUT"INPUT CO-ORDINATES OF B ",E2,N2
110 INPUT"INPUT CO-ORDINATES OF C ",E3,N3
120 B1=ATN((E1-E3) /(N1-N3) )
130 B2=ATN((E2-E3) /(N2-N3) )
140 B3=ATN((E1-E2) /(N1-N2) )
150 A4=B2-Bl
155 C4=1/TAN(A4)
160 A6=B3-B2-3.14159
170 C6=1/TAN(A6)
180 C5=1/TAN(3.14159-A4-A6)
190 K1=1/(C6-C3)
200 K2=1/(C4-C1)
210 K3=1/(C5-C2)
220 K4=K1+K2+K3
230 E4=INT((((K1*E2)+(K2*E3)+(K3*E1) ) /K4) *100+0.5)/100
240 N4=INT((((K1*N2) +(K2*N3)+(K3*N1) ) /K4) *100+0.5) /100
250 PRINT"CO-ORDINATES ARE ",E4,N4
260 END
7.8 Weisbach Two surface stations P and Q, having co-ordinates 1250.00 mE,
triangle 1200.00 mN and 1300.00 mE, 1350.00 mN, were observed during the
installation of plumb wires in a shaft. The readings in Table 7.7 were
recorded by a theodolite set up at a surface station A near to line XY.
Table 7.7
P 273° 42308,”
Q 93° 43’ 54”
Plumbwire X 08"=00" 50"
Plumbwire Y G7° 585102
Figure 7.9
Station N
(m)
1250.00 1200.00
1300.00 1350.00
Q, station A tying the two Weisbach triangles together. In Fig. 7.9 we have
the values in Table 7.8.
Therefore length PQ (50.007 + 150.007)
158.114 m.
_, 50.00
Bearing of PQ = tan
150.00
18° 26’ 05.8”.
Also, from the circle readings at A we have the values in Table 7.9.
Table 7.9
Table 7.10
6 8.3 74 x 160”
XYA ae
5.945
= 225.4”.
Also, from the circle readings, we get the values in Table 7.11.
Table 7.11
X é 8° 00’ 50”
P 273° 42’ 08”
ae oT 29.3”
deduct 180° 00’ 00”
195°.53" 25.2%
Figure 7.10
Solution. First consider the surface control (Fig. 7.10). The co-ordinates
of control points A and B have been given, so a ground point corresponding
to C could be located by methods previously discussed. For AB (Fig. 7.10)
we have the values in Table 7.12.
Table 7.12
Station Ee N
(m) (m)
A 1850.75 1270.40
B 1990.24 1310.66
Table 7.13
e N
(m) (m)
A 1850.75 1270.40
Cc 1942.00 1281.80
Table 7.14
Station E N
(m) (m)
B 1990.24 1310.66
C 1942.00 1281.80
X and Y are to be brought into this line. The greatest of care has to be exercised
and observations should be made from both sides of the shaft to obtain the
best positions.
Next establish the underground control. Figure 7.11 shows the disposition
of the underground station U as it forms a Weisbach triangle with plumb wires
X and Y, now suspended within the shaft. The relevant horizontal circle
readings are given in Table 7.15.
Table 7.15
20-120 12%
2a” 15° 35°
In triangle YUX
sin YUX _ sin UYXx
XY Ux
YU sm” UY sin 1”
3.640 5.295
The co-ordinate values have been ‘rounded off’ to be in accordance with those
of the control stations.
Next locate the centre line stations at 12.00 m and 25.00 m from C. The
first of these will be 10.18 m from X along the centre line. In triangle PUX
UX = 5.295m
PX = 10.180 m
and UXP = 679.9”.
Assume that UP + UX = PX
Table 7.16
Difference
Hence length UR = V
(5.1237 + 0.2257)
5128 m:
$= 55.123
Bearing UR = tan
0.225
272° 30’ 53.3”
Since the bearing of XU = 255° 11’ 19.9”
clockwise angle XUR = 197° 19’ 33.4”
But the horizontal circle reading for pointing on X is 23° 20’ 12”. Therefore
the horizontal circle reading for pointing on R
7.10 Construction A rotating construction laser is being used during the setting out of
lasers a car park and a reading of 0.913 m was obtained on to a TBM at level
44.13 m.
If the car park has a crossfall of 1 in 1000 longitudinally and 1 in 200 |
Solution. The first task is to calculate the levels of the four corners of the
construction bay.
Transverse rise = ay = 0.150 m
Check calculation
Offsets from the beam are obtained by subtracting the required level of the
corners from the level of the beam.
Level of TBM = 44.130 m
Level of beam = 44.13 + 0.913
= 45.043 m
Offset to lowest corner = 45.043 — 44.32
= 0.723 m
Offset to longitudinal corner = 45.043 — 44.326
= 0.717 m
Offset to transverse corner = 45.043 — 44.47
= 0.573 m
Offsets to highest corner = 45.043 — 44.476
= 0.567 m
7.11 Sight rails for a A sewer is to be laid between two points P and Q, 80 m apart. Levels
sewer were taken to establish the longitudinal profile, and the bookings in Table
7.17 obtained. Invert level at P is to be 120.750 m AOD and the sewer
is to fall towards Q at a gradient of 1 in 120.
0.633
rs
At the setting-out stage the level was set up close to its previous position
and a backsight of 0.587 m was recorded on a staff held at the bench
mark.
(a) Suggest a suitable length for the ‘traveller’.
(b) Determine:
(i) the staff readings required for the fixing of sight rails at P and Q;
(ii) the heights of those rails above ground level at P and Q.
Introduction. The sight rails are positioned so that the line connecting their
upper edges reflects the gradient of the trench bottom or the pipe invert, as
applicable. A boning rod, or traveller, of correct length is held with the upper
edge of its horizontal sight bar just in the line of sight given by the sight rails;
in this position its lower end stands at the required level (Fig. 7.12). The
horizontal sight rails are nailed to stout uprights, firmly installed on alternate
sides of the trench. These uprights must be well clear of the sides of the trench.
Frequent checking of their integrity is essential.
Sight rail
Traveller
Sight rail
Boning
rod
or
traveller
Sight rail
Figure 7.12
Height of RL Distance
collimation (m)
125.453 124.820 BM
122.528 Point P
122,796
122.893
122.538
122.433 Point Q
0.633 124.820 BM
Solution. First determine the existing ground levels. The level bookings
may be reduced as follows in Table 7.18.
Next calculate the length of boning rod. Invert level at P is to be 120.750 m
and the gradient of PQ is to be 1 in 120 falling towards Q, 80 m from P.
Table 7.19
120.750 122.528
55 120.583 122.756
40 120.416 122.893
60 120.250 122.538
80 120.083 122.433
Inner upright
yr nN Outer upright
met t +,
ne i. ting
link
Introduction. The profile boards are nailed to two uprights which are firmly
driven into the ground near the toes, A and C, of the embankment (Fig. 7.13).
The inner uprights need to have clearances from the toes of the order of 1.0 m
to prevent disturbance. The inner and outer uprights can be spaced up to 1.0 m
apart, since the sloping boards reflecting the side gradients need to be of
reasonable length for sighting purposes. A traveller is used in conjunction with
the upper surface of the boards to achieve the gradients.
Various configurations are discussed in this example.
and d= ra + sh oS
2 k+s
where b is the formation width; h is the height of embankment at centre line;
s is the side slope of the embankment; and k is the transverse cross fall.
Now h = 267.50—261.98 = 5.52 m
d, = (10+11.04) (4) = dy = (10+ 11.04) 4)
= 27.05 m = 17.21m
The reduced level of B = 261.98 m AOD
Figure 7.14
A and C can be set out by direct measurement from B, either sloping lengths
of 27.22 m and 17.32 m, respectively, being measured or the above horizontal
distances being set off perpendicular to the centre line. Alternatively, a trial
method using a level can be adopted as in Fig. 7.15. Say that a reading of
2.44 m was obtained when the staff was held at B, the level being in any suitable
position, not necessarily on the line of the cross-section. The staff man is now
moved along BC until the reading of 0.53 m is attained, i.e. to realize the
level difference of 1.91 m previously calculated. This position fixes C.
Next determine the profile board configuration. The two major considerations
are the length of the traveller and the relationship between the profile board,
existing ground and proposed slopes. In the first instance assume a traveller
of length 1.25 m with the centre lines of the uprights at 1.00 m and 2.00 m,
respectively, from toes A and C.
At A the sight line must be at an altitude of 1.25 + 258.97 = 260.22 m.
Thus at the uprights the corresponding altitudes of the sight line will be
Table 7.20
Cc Inner Outer
upright upright
0.03 m at the outer upright is too small, but if the spacing of the uprights
is reduced to 0.60 the difference at the outer upright becomes 264.34 — 264.07
= 0.27 m, and this difference should be considered to be a minimum, allowing
for the width of the board itself. Alternatively, if the traveller is made 1.5 m
long, which is probably a maximum for convenience, the values in Table 7.21
apply when the uprights are 1.0 m apart.
Table 7.21
A Inner Outer
upright upright
At C C Inner Outer
upright upright
Problems
2 List the various types of lasers used for mine surveying purposes.
Describe the particular applications they can be used for, and what their
limitations are both economically and practically.
With the aid of a neat sketch, show a suitable arrangement for the
use of an alignment laser in tunnelling.
[Eng. Council]
3 Control station X lies inside triangle PQR such that PXQ = 128°
23’ 24” and QXR = 97° 10’ 11”. Calculate the co-ordinates of point
X if the co-ordinates of
P are 1798.2 mE 7643.2 mN
Q are 2534.8 mE 6327.1 mN
R are 1531.7 mE 6215.6 mN
[Salford]
Answer 1916.7 mE, 6701.6 mN
Table 7.22
7063.122 1110.126
7632.157 827.553
5 (a) During the setting out of a tunnel the observations in Table 7.23
were made by theodolite at a surface station I near to a vertical shaft.
Table 7.23
The plumb wires down the shaft were 5.345 m apart, P being the nearer
of the two to the theodolite at I. If the whole circle bearing of PQ was
307° 47’ 24” deduce that of AB.
(b) Explain briefly two methods of using the gyro-theodolite to
determine the azimuth of AB. [Salford]
Answer 214° 14’ 50”
(a) Calculate the data to set out the first three columns from T using
a theodolite and steel tape.
(b) If the theodolite axis was set up 0.90 m vertically above T and
the columns terminated at 0.6 m below the ramp surface, calculate the
vertical angles required to check the stopping-off heights of the three
column heads. [I. Struct. E.]
Answer Chords 14.84 m long, deflection angles 14° 19’ 26” , 28° 38’
52” and 42° 58’ 18”; vertical angles 00° 00’ 00” ,02° 59’ 06”, 04°
143"
Table 7.24
Table 7.25
Table 7.26
E ~ OW’
(m) (rm) :
x 1000.0 1000.0
4 1300.0 1500.0
A 1100.0 1050.0
B 1180.0 1200.0
Cc 1325.0 1400.0
11 A, Band C are three control points which lie in a straight line and
have co-ordinates given in Table 7.27. The horizontal circle readings
on a theodolite stationed at T,, a point to the east of the line AC, were:
Table 7.27
Station
Table 7.28
Table 7.29
14 A and B are points on the centre line of a level mine roadway and
C and D are points on the centre line of a lower roadway having a uniform
gradient between C and D. It is proposed to connect the roadways by
a drivage from point B on a bearing of 165° 35’. Given the data in Table
7.30, calculate
Easting Northing
(m) (m)
Table 7.33
Station E N Hy
(m) (m) (m AOD)
Table 7.34
Table 7.35
(a) Book and reduce the levels carrying out appropriate arithmetical
checks.
(b) Draw a longitudinal section of the ground.
(c) Calculate the invert level of MH ‘A’ and the gradients of sewer
runs from MH ‘A’ to MH ‘B’ and MH ‘B’ to MH ‘C’, and show this
information on the section. Determine the minimum cover to the sewer.
(d) A traveller 3.0 m long was used when the trench was excavated.
Table 7.36
BS IS FS Distance Remarks
(m)
1.480 50 BM RL 48.290
0.855 75 Sight rail at MH 12
2.105 40 Sight rail at MH 13
1.480 50 BM.-RL 48.290
eastwards of the plumb wires and with X the nearer in each case. The
plumb wires were 6.250 m apart and the whole circle bearing of PQ
was 210° 17’ 22”. Determine the bearing of SR. [Salford]
Answer 107° 57’ 08”
Table 7.38
Error types A constant error, such as an incorrectly standardized tape, has the same
influence throughout the programme of observations, whilst a systematic error,
like the influence of field temperature on the length of that measuring tape,
can be of varying magnitude and changing sense. Accidental errors, such as
instrument imperfections, are of a random nature and could be said to be second-
order errors when compared to the types just mentioned.
In this chapter a number of examples are solved on the assumption that the
observations are free of all errors other than accidental errors. In other words,
the work has been carried out carefully, with as many checks as possible to
eradicate mistakes and, in addition, corrections have been applied to eliminate
any constant errors and systematic errors.
Accidental errors Accidental errors are unavoidable in so far as the observer is concerned. They
usually occur with positive and negative values of the same size and frequency,
the larger errors occurring less frequently than the small errors. If a number
of measurements of the length of a line are repeatedly made under similar
conditions a spread of values will be obtained and the difference between the
largest and the smallest values is defined by the term ‘range’. The values can
be placed into groups or classes having specified limits, the number of measures
in each group being termed the class frequency. If these numbers are expressed
as a percentage of the total then proportional frequencies are obtained.
Alternatively, the numbers themselves, taken in conjunction with their
respective classes, give the frequency distribution of the whole set of
measurements.
Probability If a very large number of measurements has been taken the frequency
distribution distribution could be considered to be the probability distribution. The normal
distribution is usually assumed for survey observations, being expressed by
the equation,
1,
1
ie Lene e 7M —p)*/20? dx
“. o V(2n)
where dp is the probability that value x will lie between the limits of x, and
(x; +dx); p is the true mean of the population; and o is the standard deviation.
Standard deviation Standard deviation is a measure of the spread of a distribution and for the
population, assuming the readings are of equal reliability,
deg Ee] n
a [PP
population standard deviation is estimated as
n—1
in which n is the number of observations and
ne ax
n
Confidence limits Having established the sample mean as an estimate of the true value of the
quantity, the observer might wish to state the range of values within which
that true value should lie for a given probability. This range is termed the
confidence interval, its bounds being the confidence limits. Having calculated
the standard deviation for a set of observations, confidence limits can be
established for that stated probability, statistics tables being available for this
purpose. Frequently, a figure of 95% is chosen, and this implies that nineteen
times out of twenty the true value will lie within the computed limits. If normally
distributed small random errors obtain, the presence of a very large residual
suggests an occurrence to the contrary. Such an observation can be rejected
if its residual error is larger than three times the standard deviation.
Weight When the measurements are not of equal reliability weights (w) are applied
to the individual measures to reduce them to one standard. The most probable
value is then the weighted mean.
The relationship w, 5,7 = w, s.2> =... = s’ applies, in which w, and
s; are the weight and standard error for observations 1, w, and sy are the
weight and standard error for observations 2, whilst s is the standard error
for the observations having unit weight. Hence .
s? s? ‘
Wi = >) Wo Se CIC,
512 592
Wi sy”
and = > ete.
W2 Sy
Thus the weight of an observation can be taken to be inversely proportional
to the square of its standard error. The most probable value of the weighted
mean (%,,) of weighted observations is
“wy Xj
Ww = Ed
“wy
mee | ial
m=)
(b) Standard deviation of an observation of weight w,,
= ot
|See
Ew, (x—%,)?
Te GC
w,(n—1)
Dwi (x —X,)"
= + Se ai i
(n—1)<xw,
Method of least This is a widely used method of determining the most probable values of
squares observed quantities, assuming that only accidental errors are present. It states
that the sum of the weighted residuals squared will be a minimum,
i.e. L[weight x (residual)”] is to be a minimum.
For n observations
Wy Tye Wola? ke Wats” + — + Wh Tn?
is to be a minimum. Thus
pah4)”- 0 5 eee
or r'Wr is to be a minimum. Note that the alternative form of L(Vw, 1)?
can be minimized.
8.1 The distribution The following readings were repeated on the same angle
ofa set of
of readi
eadings 73° 40’ plus (22 152, 09”, +42. 10*. 1*" 16°. is
ED? ai Lo
Determine:
Vy }™
(a) the most probable value of the angle, Cm
—
~ SY
Table 8.1
ae Oe ae
n
, } \ . 2 :
) ee Je = re ~~ 84
rh a mS = 9
i = + 3.1”
S
(d) Standard error of the mean = +
Vn
or
an
=)
eha MME Ares
" 1 Ce} . ; t W) ch 34
vane, ai
Cr . tats
rie a
Most a ee
Proboble oe
Vlog ok/10 .
= + 1.0”.
(e) 95% confidence limits:
aie ix s
The lower confidence limit = % —
Vn
Pea bX S
The upper confidence limit = *+
Vn
t is selected from statistical tables for a given value of n. When n = 10,
t = 2.26 and so
be ak
| = 2.26 x ——
Vin V10
= 2.2”.
The 95% confidence limits are 73° 40’ 14.0” + 2.2”.
It is common practice in surveying to reject any observation that differs from
the most probable value by more than three times the standard deviation.
8.2 Relative Prove that the mean value of a set of observations satisfies the least
precision squares condition.
The length of a line was measured using two different EDM instruments
A.and B with the following results:
Table 8.2
EDM B
785.546 785.545
785.538 785.549
785.535 785.543
785.550 785.546
785.542 785.545
785.547 785.550
785.547
785.543
xx \2
nia — :
n
which is always positive. Hence the sum of the squared residuals based on
the arithmetic mean is always less (i.e. a minimum) than the equivalent based
on any other value.
Table 8.3
785.546
785.538 7ws i
785.535 —8 64
785.550 +7 49
785.542 ay ‘s
785.547
eee “be
+ 5.73 mm.
Ws
= Sj = beeet + 2.34 mm.
v6
Similarly standard deviation Sg,
785.545
785.549
785.543
785.546
785.545
785.550 ah,
785.547
785.543 —
o-n-O00O
Lix—#)" 46
Sp= + ee +
n—-1 g~1
+ 2.56 mm.
2.56
eae
= + 0.91 mm.
Next calculate the relative precision. The length of the line given by EDM
A is 785.543 m with a standard error of 2.34 mm, whilst by EDM B it is
785.546 m with a standard error of 0.91 mm. Now
i a sa? * 0.917
Wp sa’ 2.34?
frout |
6.6
This ratio is a measure of the relative precision of the instruments.
Finally calculate the most probable length of the line. The most probable
length of the line is the weighted mean of the two observed lengths. Now
es LW x
Lw ;
In this example
a Wa ON te WR Iz
ww ’
Wa + Wp
in which La, and Lg are the mean lengths recorded by EDM A and EDM B,
respectively.
WB
——«= W
6.6 4
785.5456 m.
Xy could be written as 785.546 m to the nearest millimetre, in accordance
with the observed measurements. In practice, the arithmetic mean and standard
deviation or error may be expressed to one decimal place further than that
to which the measurements had been made.
Table 8.5
233
= 73° 40’ + ——
A
73° 40’ 13.7”.
(b) Standard deviation (s) of an observation of unit weight
n—1
TIS ;
s=t (since n = 8)
(S=1)
= + 3.33”.
(c) Standard deviation of an observation of weight 3.
Therefore $3 = + =
= + 1.92”.
(d) Find the standard error of the weighted mean.
+
| SS
77.53 |
(8-1) x 17
= + 0.81’,
or, since We Sn eo"
S 3133
Sm = Rit = + orl
= + 0.81”.
8.4 Probability Taking the standard deviation o of a single measure in Example 8.1
distribution to be +3” calculate:
(a) the magnitude of the deviation likely to occur once in every two
measurements; i
dp =
1 e ~ @1-#)"/20? dx.
P wN2n
Bic X— pb ,
On writing u = the expression becomes
oO
dp = we en"2 dy,
V20
This is the standardized form of the above expression, and the relationship
between dp/du and uw is illustrated in Fig. 8.1.
Figure 8.1
The curve is symmetrical and its total area is 1, the two parts about u =
having areas of 0.5. The shaded area has the value
oe ——e du
_» V20
and it gives the probability that u lies between — oo and +1, i.e. less or equal
to u;. The unshaded area gives the probability that u will be larger than +1).
Since the curve is symmetrical, the probability that u takes up a value outside
the range +u, to —u, is given by the two areas indicated in Fig. 8.2.
The values of the ordinates of the standardized form, and the corresponding
definite integrals, have been determined for a wide range of u. These have
been tabulated and are available in various publications. Typical values used
in this example are given in Table 8.6.
Solution. (a) For a deviation to occur once in every two measurements
a probability of 50% is implied. Thus in Fig. 8.2 the two shaded parts have
areas of 0.25 each and the shaded area is 0.5.
Bearing in mind that |*%! is the shaded area in Fig. 8.1, we require a value
of u in the table such that
+u; +u
ee = 0.25, ice. | = 0.75.
—oc —o
By inspection, u lies between 0.6 and 0.7, for which the integral has values
of 0.7257 and 0.7580, respectively. Hence the value u for the integral to be
0.75 is 0.6745.
0.6745 {= |
n—-1l
+2.0
| = 0.9772.
( +2.0
Hessel | =-1-09772
= 0.0228.
For the deviation to lie at the limits of, or outside, the range +6” to —6”
the probability is 2 x 0.0228 = 0.0456, or 4.6%.
(c) The standard deviation of the mean of nine observations
ae O54 A ip 3.0
J V9 3
= + 1.0”.
x—p |
Va
—— = ——
o 1.0
= Si?
1.5
whence | == (9332
—o
1.5
l— | = 0.0668.
8.5 Combinations of The co-ordinates of two surveying stations A and B are given below,
errors together with their standard deviations. Calculate the length and the
standard deviation of AB determined from these values.
Station A Easting 456.961 m + 20 mm,
Northing 573.237 m + 30 mm
Station B Easting 724.616m + 40 mm,
Northing 702.443 m + 50 mm
If an independent measurement by taping gives the length of AB to be
297.426 m + 70 mm and separate determination by EDM to be 297.155
+ 15 mm, calculate the most probable length of the line using all the
available information. What is the standard deviation of this value?
[Bradford]
in which s,, 5p, S- and sg etc. are the standard deviations (or the standard
errors) of a, b, c and d, respectively. Similarly, ifa =b+c+dt+...
0
,° = sy? +8,” + sy + ..., since es etc. = l.
[8.2]
: 0
Since’ (=) = K. [8.3]
0b
In the above relationships it is assumed that b, c, d, ... are independent,
implying that the probability of any single one having a certain value does
not depend on the values observed for the others.
Table 8.7
Solution. First determine length AB from the co-ordinates (see Table 8.7).
Therefore length of AB = |(267.655* + 129.2067)
= 297.209 m.
Next calculate the standard deviations of the easting and northing differences.
The easting difference was calculated as (724.616 — 456.961) m. Thus the
standard deviation (s,) of this difference is given by
5,2 = sep? + Sea7, as in eq. [8.2]
s,° = 40? + 20?
SF = Sia" +5 Sin
= 50° + 30?
sy = 3400
= 58.3 mm.
Now calculate the standard deviation of the calculated length AB.
Lis = x+y"
Lap = (x? + y?).
Whence, from eq. [8.1],
SAB" = oa ae * Oa ; ae
Ox dy :
in which s, = 44.7 mm and s, = 58.3 mm.
Now —
AB
aL
ax
= 1 x ay(x? + 2)" 2 > _ = Lap3
Ir +»)
OL
and — = +x 2y(x* ty?)
7 = — ; : “Files z. :
dy VQ" + y*) Lap
267.655 \? 129.206 \?
Whence § sap? = -(————) x 44.77 + eee x 55.3°
297.209 297.209
= 1622 + 642 = 2264 mm’.
Therefore sap = 47.58 mm,
297.426 | 1/4900
297.155 1/ 225
297.209 1/2264
Gea: Her,
= 297.171 m.
1 ae
1 a.
1
Now Wy = Yw =
4900 Pe) 2264
= 0.005 09 mm~?.
Since weights were applied as the inverse of the squares of the standard
deviations
‘bel
WE 5, :
5 1
Therefore s 2S
0.005 09
and Ss = + 14.0 mm.
Standard deviation of length 297.171 m is + 14.0 mm.
9.375 m 9.493 m
H, H,
Figure 8.3
= ee
: =tagyew henen 66 is
i sma ll
Therefore,
—b —H?
dH = we dé = dé (dé in radians)
Writing Spr, Spa, Sgr, Sy and Sg as the respective standard errors we have
SpoPQ = — —B
Hj Sq
9.3757 ° 1
1 206 265
0.00426 m
qe
1 206 265
= 0.000437 m
8.7 Accuracy of Two vertical levelling staves were securely fixed in position at stations
levelling A and B, 93.00 m apart. Ten instrument stations were then located
between A and B as shown in Fig. 8.4 and the differences in height
between A and B were then measured by a tilting level set up in turn
at the ten stations (see Table 8.9).
Determine the most probable values of the height differences between
A and B and the inclination of the line of sight of the level when the
bubble is centred.
OOaARON
COON
=
P, 4
Line of sight
Figure 8.4
Difference
in height
0.150 y = 0.1238 — 0.000 3045 D
Figure 8.5
Solution. First derive a general expression. In Fig. 8.4 it has been assumed
that the line of sight is inclined upwards when the bubble is central.
Apparent difference in level between A and B = a, — by.
True difference in level between A and B =a — Dp.
If t be the best, i.e. most probable, value for (a—b) and @ be the best value
for the inclination of the line of sight, the best value for the difference in level
(Y) obtained by observations from any instrument station will be Y = tf +
D, 6, and this is the equation of the regression line to fit the measurements.
Thus for an observed or apparent difference of y,; a residual of value
y,—(t+D, 9) results. Applying the principle of least squares the values of ¢
and @ which minimize the sum of the square of the residuals (or deviations)
between the observed measurements y, and Y, are found. L[y, —(t+D,6)]’
is to be a minimum. Writing this sum as E
E = X(y,’ — 2y,(t+D, 6)+07+2t D, 0+D,? 67).
dE
For a minimum “le = 0 and — =0.
Ot 00
Therefore —L2y, + L2t + L2D, 6 =0
and —L2y, D; + L2tD, + L2D,? 6 = 0.
Table 8.10
—0.0013
+0.0016
+ 0.0009
+ 0.0006
—0.0012
—0.0003
— 0.0004
+ 0.0004
Ly,
Therefore t = (n = 10 in this example),
n
i.e. the best value for t is the mean of the ten differences measured and
j= Ly,D,
2 ea
The data can be tabulated as in Table 8.10.
Most probable difference in height
t=
ry;
10
1.238
10
0.1238 m, say 0.124 m.
Most probable value of the inclination of the line of sight
LyD,
g=—
ED,”
6.702
22010
— 0.000 304 5 radian, i.e. 0.000 30 radian
(inclined downwards).
8.8 Least squares — A series of level networks are run, and the results are as follows:
the normal equations Weight
Height of A — Height of B = 18.614 3
Height of B — Height of C 16.264 Z
Height of C — Height of D = 22.385 2
Height of A — Height of D = 57.247 1
Calculate the most probable differences in height between A and C and
between B and D. [Leeds]
Introduction. When using the method of least squares to solve such prob-
lems as this two approaches present themselves:
a
Datum
Figure 8.6
r, = 57.247 —a
ry = 18.614 — (a—b)
r3 16.264 — (b-—c)
rg = 22.385 — c.
Dr,” is to be a minimum, i.e.
R = (57.247—a)? + [18.614—(a—b)]* + [16.264—(b—c)]?
+ (22.385—c)?, is to be a minimum.
aR aR aR er
Hence —— = 0, —— = 0 and —— = 0, S>5 a
da db dc Rk + ng
OO
and so shin:Dee oe toa Recon Wy,
(b) Adopt the assigned weights. The residuals r,, etc., maintain the same
values but they do not have the same reliability or precision. Hence Lwyry-
is to be a minimum in which
ry = 57.247 — a (weight 1),
ry = 18.614 — (a—b) (weight 3),
r; = 16.264 — (b—c) (weight 2),
rg = 22.385 — ¢ (weight 2);
and R = (57.247—a)* + 3[18.614—(a—b)]* + 2[16.264—(b—o)]?
+ 2(22.385—c)?.
OR OR OR
Again
—— = 0, —— = 0 and = 0, and so with respect to a
da 0b dc
(—2 X 57.247) + 2a — (6 X 18.614) + 6a — 6b = 0
with respect to b
(+6 X 18.614) — 6a + 6b —(4 X 16.264) + 4b — 4c = 0
with respect to c
(+4 X 16.264) — 4b + 4c — (4 X 22.385) + 4c = 0.
Dividing throughout by 2 results in the normal equations
4a — 3b — 113.089 = 0
~3a +5b — 2c_+ 23.314 = 0
— 2b + 4c — 12.242 = 0.
Whence a = 57.254 m, b = 38.642 m, c = 22.381 m. Therefore the most
probable difference in height between A and C is 34.873 m and the most prob-
able difference in height between B and D is 38.642 m.
Table 8.11
+ 57.247
+ 18.614
+ 16.264
+ 22.385
for the first line, (w) (c,) (term a) + (w) (c,) (term N),
i.e. (1)((—1) (—Da + (1)(—1) (57.247);
(3)(—1)(—
Da + (3)(-1)(1b
+ + (3)(—1) (18.614).
These summarize as
a — 57.247 + 3a — 3b — 55.842,
and so we obtain the normal equation
4a — 3b — 113.089 = 0.
For b, which appears on two lines
(3)(+1)(—Da + (3)(+1)(+Db + ()(+ 118.614)
(2)(—1)(—1)b + (2)(—1)(+
De + (2)(— 116.264)
which give:
Further details on the matrix method of solving this type of problem are given
in Example 8.10.
The DIM statement in line 20 is set for a problem with 10 unknowns in
10 equations. For readers with computers with limited memory the minimum
values for the array variables are: A(N2, N1); B(N2); D(N1, 2* N1); E(N1);
T(N1, N2); W(N2).
50 FOR I = 1 TO N2
60 PRINT "FOR EQUATION"; I
70 FOR J = 1 TO Nl
80 PRINT "INPUT COEFFICIENT OF VARIABLE"; J
90 INPUT A(I, J)
100 NEXT J
110 INPUT "PRODUCT OF THIS EQUATION "; B(I)
120 INPUT "WEIGHT OF THIS EQUATION "; W(I)
130 NEXT I
140 FOR I = 1 TO N1
150 FOR J = 1 TO N2
160 T(I, J) = A(J, I) * W(d)
170 NEXT J
180 NEXT I
190 FOR K = 1 TO N1
200 FOR I = 1 TO Nl
210 FOR J = 1 TO N2
220 .D(Ky Ly =wD (Kea D) et CE (K, i), * Ado per))
230 NEXT J
240 NEXT I
250 FOR I = Nl + 1 TO 2 * N1
260 IF Nl + K = I THEN GOTO 290
270 D(K, I) = 0
280 GOTO 300
290 D(K, I) =1
300 NEXT I
310 FOR J = 1 TO N2
320 E(K) = E(K) + (T(K, J) * B(J))
330 NEXT J
340 NEXT K
350 FOR K = 1 TO N1
360 FOR J = 1 TO N1
370 IF J = K THEN 420
380 M = D(J, K)
390 FOR I = 1 TO 2 * N1
400 D(J, I) = D(J, I) * D(K, K) - D(K, I) *M
410 NEXT I
420 NEXT J
430 NEXT K
440 FOR J = 1 TO N1
450 X = D(J, 3)
460 F = 0
470 FOR I = Nl + 1 TO 2 * Nl
A480 Did, I) = D(d, 1). / X
490 F = F + D(J, I) * E(I - N1)
500 NEXT I
510 F = INT(F * 10000 + .5) / 10000
520 PRINT "VARIABLE"; J; "="; F
530 NEXT J
540 END
+ 42.285
— 12.016
— 20.240
~ 8.114
— 22.224
Figure 8.7
[Salford]
Solution. The network of levels together with the directions of the rises
and falls are shown in Fig. 8.7. A study of this figure shows that there are
three closed circuits of level runs, namely ABCDA, ABDA and BCDB. Two
out of the three are needed in the analysis, but whichever are chosen the
‘directions’ of levelling should be compatible. For instance, in Fig. 8.8,
ABCDA and ABDA are compatible. Choosing ABDA and BCDB requires
some adjustment for conformity along BD and the falls along BC and CD are
transformed into rises to satisfy the changes in direction in Fig. 8.9. Such
transformations should minimize any confusion in respect of signs. In each
of the closed circuits it is required that the sum of the ‘rises’ will equal the
sum of the ‘falls’. Accordingly there is one condition equation for each circuit.
Let the configurations of Fig. 8.9 be adopted in the analysis and let
corrections (equal and opposite to the corresponding errors) e;, €2, 3, €4 and
es be required for the corresponding adjustment.
In circuit ABDA the error+42.285 — (20.240 + 22.224)
—0.179 m.
In circuit BCDB the error +(12.016 + 8.114) — 20.240
= — 0.110 m.
G
— 20.240
— 22.224 , gue
D D
Figure 8.8
A 42.285 3
B B 5
12.016
—22.224 Saseay Cc
e, es
8.114
e3
D D
Figure 8.9
The corrections are to be such that the two condition equations will be satisfied
and these reduce to
é; + e4 + es = +0.179m
eo + e3 + es = +0.110 m.
Since the equations of condition must remain unaltered by increments de, 6e),
etc., ;
Two runs were made on line AB and the mean value of level difference has
been used above.
It has been mentioned in Example 8.2 that
ele
10
Therefore e, = 10x, e, = 10y, e; = 20y, eg = 30x and
es = 25x + 25y.
Substituting in the equations of condition, which have to be satisfied
LOxet) 30x (25x te 25y)' = 0.179
10y + 20y + (25x + 25y) = 0.111
andso 65x + 25y = +0.179
25x + S5y = +0.110,
whence y = 0.000 91, x = 0.002 41,
e, = 10x = +0.024 m
e> = 10y = +0.009 m
e3 = 20y = +0.018 m
eg = 30x = +0.072 m
es = 25x + 25y = +0.083 m
8.10 Least squares From an observation station O four trigometrical stations A, B, C and
— angle observations D can be seen. The readings in Table 8.13 were taken in order to
determine the most likely values of the included angles.
Table 8.13
26° 161.221
37° 15’ 33°
47° 22 15”
62° 32’ 05”
84° 37’ 41”
109° 54’ 21”
and so
xe Ts alIgpie yt+z i: = y
Hence, 2) = —— = ee ae = pene a,
Wi W2 W3 4 Ws
a4
gE.
W6
Inserting the values for the weights and incorporating the above in eqns [8.9],
[8.10] and [8.11] gives
z
ete) hi Gyre) +O rz) + 3 = +06”
Q5,Xa op Ky FF an: + ay x; = by
in which n>i. If x; x) ... x, are the best values, residuals can be written as
ry = —(@y% + ay. % + ... +.a5;
x) + Dy
fap NS —(a> X{ se az2 X2 ci ne ae Ay; Xi) 3 by
Fo (Oy eX ay Xe a en dlsKi). he
and expressed as r = —Ax + b.
Allowing for weighted observations the principle of least squares requires
(—Ax + b)’ W(—Ax + b) = x" A'WAx — b'WAx
— x'A'Wb + b'Wb
to be a minimum,
Now let AOB = a, BOC = band COD =. In the example six observation
equations can be established.
Hence
=1 10° 270 3
Ov O22 3
Oo O-t.09.3
A'WA =|6 5 3
585
356
=[42 5 41 eT eH
-2 +2 -3 11.541 9009
+3 -2 +2 9.535 509 ‘
Note: The reader will realize that having determined the six residuals from
the observations, i.e. 25° 16’ 17” — a, 37° 15’ 33” — b, etc., a solution
could be derived using the method given in Example 8.8. The computer
program listed with Example 8.8 will solve this problem.
8.11 Braced In order to determine the co-ordinates of the two ends A and B of the
quadrilateral by least centre-line of a bridge across a wide river, a braced quadrilateral ABCD
squares and equal was set out. The mean observed angles and their log sines are given in
shifts Table 8.14. Adjust the angles.
[ICE]
The corrections also influence the log sine values. Before adjustment the sums
are
log sin 1 = 1.840 244 2 log sin 2 = 1.795 251 1
3 = 1.746 232 2 4 = 1,952 506 5
5 = 1.880 041 1 6 = 1.737 067 2
7 = 1.885 386 7 8 = 1.866 872 9
Cp 2 + 63 + 4 + es + Cg + er + Cg = 32”.
The total correction is shared equally to each of the eight angles, i.e.
e= +4”,
Step 2: eqn [8.17] is stated as
ey te PE eee fn Oe
Le, = +32”.
At this stage the corrections in Table 8.15 will satisfy the requirements of eqns
[8.16], [8.17] and [8.18].
Step 4: the log sine condition determined in (a) stated
1.351 904 2 fe eV) = C3V3 Ap e5V5 aie e7V7
Table 8.15
—
OANoahWND
Table 8.17
8.12 Variation of In Fig. 8.11 two control stations A and B have the co-ordinates given
co-ordinates in Table 8.18.
Table 8.18
Ee
ln pn
A 5210.15 12267.92
B 8785.62 9686.45
y+
| 76° 6’ 29” 4663.08
Figure 8.11
XG—*A
tan Gag =
YG— YA
and the change in that bearing due to the changes dx,, dya, etc., is
_ Ga=
Ge i Xa) dye
la
Similar expressions can be derived for BG, whilst the change in angle BGA
can be related to da,g and dagg. In the method of variation of co-ordinates
provisional co-ordinates are allocated to points requiring adjustment. These
points are subjected to displacements dx and dy duly determined by a ‘least-
squares’ analysis. Residuals (7) can be derived in the form
r = O-C-dy,
where O is the measured value of a length, bearing or angle, C is the value
of that quantity as calculated from relevant co-ordinates, and dy is the change
in that quantity due to the displacements of the respective points. (C+dy) is
the best value for the quantity.
Solution. (a) Calculate values of (O—C) for bearings aag and agg.
4526.39
soiatiaa SE EW72
Therefore Qag = 76° 06’ 17.5” =C
‘Therefore (O—C) = 76° 06’ 29” -— 76° 06’ 17.5”
= +11.5”
Xg = 9736.54 yg = 13 387.68 (provisional)
Xp = 8785.62 yg = 9686.45 (fixed)
4526.39 1119.76
=————
4662.84 dig 9 + ————
4662.84 dg
= 0,970 736.7 dxg + 0.240 145 5 dy,
. Aree gE , 3701.23
a ee aida 3821.43 °°
0.248 838 8 drg + 0.968 545 8 dyg.
Also (daag— dage) = —0.000 201 9 dig — 0.000 143 1 dyg.
4.7
r= — —— — 0.000 253 4 dx, + 0.000 065 1 dyg,
206 265
For lac
Table 8.19
ce a}?
Bearing 51 566.25
4.0
2
Length (: ) 20.00
0.05
2
Angle = oe 36 833.04
5.6
(d) Normal equations can now be derived. The method discussed in Example
8.8 will be adopted to minimize L[Vw, r,]* with respect to the two unknowns
dxg and dyg. In Table 8.20 the five residuals for a,c, agg, lag, lgg and angle
BGA have been multiplied by the respective values of Vw to obtain the specific
values within the columns. The two normal equations are obtained from the
table as
634.799 172 6 dxg + 156.472 281 3 dyg 48.697 989 3
156.472 281 3 dxg + 552.614 645 3 dy —99.665 040 7
and these solve for a
8.13 Trilateration A total station instrument is set up at point P and horizontal distances
to three control stations A, B and C are recorded as follows:
PA = 63.177 m + 0.003 m ai
PB 97.710 m + 0.004 m
PC = 136.930 m + 0.005 m
For the third step the residules have to be determined. It is thus necessary
to establish the changes in the lengths due to the variation of the co-ordinates
of point P. In Example 8.12 we established the expression
dig = — (Gg —Xa) dx, — (YG —Ya) dy, + (XG —Xa) dx,
lac AG AG
4+ Oa=Ya) aye
lac
This can be altered, with reference to Fig. 8.11, to read
dlag = —sin bac dx, — cos OaG dy, +sin OaG dig
+ cos OaG dyg
Whence,
823.307 — 567.596
— 969.066 — 246.802
105.382 994.432
For e,
(823.307) (1/9) (823.307e —567.569n) + (—969.066) (1/16)
(—969.066e —246.802n+ 1) + (105.382) (1725)
(105.382e+994.432n—12) = 0
75314.935e —51922.862n + 58693 .057e+ 14947 .964n — 60.567 +
444 .219e+4191.849n — 50.583 = 0
134452.210e —32783.049n —111.151 = 0
For n,
Variables
A = WCB to set co-ordinates F(i) = Calculation matrix
Al = Angle in triangle to set I, J, K = Loop counter
co-ordinates 1b = Calculated distance to
A(i, j) = Coefficient matrix point P
B(i) = Product matrix L(i) = Measured distances to
Cl = Cosine of angle Al point P
DG, j) = Calculation matrix M = Line divisor in matrix
E, N = Co-ordinate of point P inversion
E(i), NG) = Co-ordinates of control T(i, j) = Transpose matrix
stations W(i) = Weight matrix
F = Output matrix element X = Divisor of each line in
matrix inversion
10 REM TRILATERATION
BOD LM GAS exe epi C3 ine DCD iOpen keKodiget Olu. US )4 NCS), 23, I, WES)
DowP = 3 1AL59
30 FOR I =.1 To 3
40 PRINT "INPUT CO-ORDINATES OF "; I
50 INPUT E(I), N(I)
60 NEXT I
70 FOR I = 1 TO 3
80 PRINT "INPUT DISTANCE TO STATION "; I
90 INPUT L(I)
100 INPUT "ACCURACY OF MEASUREMENT ", W(I)
110 W(I) = 1 / (W(I) * W(I))
118 NEXT I
11D) ASAIN CCE 2) SEL) G7 AN (2) = eN(2)
120 IF E(1) > E(2) AND N(2) > N(1) THEN A = A + 2 * P
121 IF N(1) > N(2) THEN A = A + P
123 Ll = SQR((E(2) - E(1)) * 2 + (N(2) - N(1)) * 2)
WRC = CCN All) Aa KEE e— (U2) * (2))) f (2 = bE eres)
125 Al = ATN(SQR(1 - Cl * C1) / C1)
126 A=A + Al
U2 7) EPVAL >ii2u* .P THEN AvH yA = 2k 5 P
128 E = E(1) + L(1) * SIN(A)
129 N = N(1) + L(1) * COS(A)
131 FOR I = 1 TO 3
132 L = SQR((E - E(I)) * 2 + (N - N(I)) * 2)
133 A = ATN((E - E(I)) / (N - N(I)))
134 IF E > E(I) AND N(I) > N THEN A =A + 2 * P
135 IF N > N(I) THENA =A +P
136 A(I, 1) = SIN(A)
137 A(I, 2) = COS(A)
138 B(I) = L - L(I)
139 NEXT I
149 FOR I = 1 TO 2
150 FOR J = 1 TO 3
T6OMT( Lp) =TA(T eye Wd)
170 NEXT J
180 NEXT I
190 FOR K = 1 TO 2
200 FOR I = 1 TO 2
210 FOR J = 1 TO 3
220) D(K, 2) = D(K, =) + (E(K, ol), * A(d, 2))
230 NEXT J
240 NEXT I
250 FOR I = 3 TO 4
260 IF 2 + K = I THEN GOTO 290
270 D(K, I) = 0
280 GOTO 300
290 D(K, I) =1
300 NEXT I
310 FOR J = 1 TO 3
320 F(K) = F(K) + (T(K, J) * B(J))
330 NEXT J
340 NEXT K
350 FOR K = 1 TO 2
Problems
3 (a) State the law of propagation of random error and name the
limitation on its application. Hence show that if o is the standard error
of each of a set of n observations, the standard error of the mean equals
alNn.
ToB
———
W,
Figure 8.13
(b) Figure 8.13 represents a Weisbach triangle and the first leg of an
underground traverse. W,W> was measured to be 2.965 m and AW,
2.097 m with standard errors 3 mm and 2 mm respectively. Angle
W,AB is known with a standard error of 4”. The approximate value
of angle W,AW, is 0° 06’ 50”.
If the standard error of a single measurement of angle W,AW, is
estimated to be 10”, calculate how many times this angle must be
measured in order that the standard error of the bearing AB is not to
exceed 5”. (Assume that bearing W,W> is known without error.)
[London]
Answer 6
Table 8.22
1
1
1
1
1
1
1
1
1
1
Calculate the three individual bay lengths, the total length of the survey
line and the associated standard deviations of the means. Compare the
subtense bar and electromagnetic distance measurements and state, with
5 The mean values of four rounds of each of the three angles of the
triangle ABC, using three different theodolites with the same observer,
are as follows:
Using a theodolite reading to 20 seconds:
ABC = 64° 13’ 50” BCA = 61° 39’ 25”. CAB = 54° 06’ 45”
Using a theodolite reading to 5 seconds:
ABC = 64° 13’ 43” BCA = 61° 39’ 27” CAB = 54° 06’ 47”
Using a theodolite reading to 1 second:
ABC = 64° 13’ 45” BCA = 61° 39’ 29” CAB = 54° 06’ 49”
If the co-ordinates of A and C are 24 261.46 mE, 19 015.83 mN and
26 472.88 mE, 20 838.05 mN respectively and B lies to the north of
the line AC, calculate the most probable value of the co-ordinates of
station B. [Leeds]
Answer 24 085.52 mE, 21 810.82 mN
Table 8.23
Instrument A Instrument B
1056.429 1056.435
1056.435 1056.435
1056.440 1056.432
1056.432 1056.436
1056.433 1056.435
1056.436 1056.436
1056.434
1056.435
13. An angle was measured 15 times and the following data recorded.
What is the most probable value of the angle?
18° 59% SUC OT OU 18 S38 *
19> OO POO. 155 Soon", 19° OO Or”
1S SS. SO en LD, oO 06 18 pei SYS
197008 OOS 185959 580 Fis” See 59"
19°00" (00. T8559 758 0 18°. SP 57." [Salford]
Surveyor 1 Surveyor 2
17 A line of length D was set out and then divided into n bays of length
d,, dy, ..., d,. Each individual length D, d, etc., was then measured
by an EDM instrument. Show that:
(a) the additive constant of the device can be calculated from the
expression
Dp
Ld,
n—-1
(b) the standard deviation of that constant is
: V(n+1)
n—-1
s being the standard deviation of a single measurement of distance.
Table 8.26
Measured angle Weighting
Table 8.27
Station mE mN
1010.124 680.127 |
1270.736 750.132
850.037 541.891
Table 8.28
Length | Error
PA 313.073 +0,005
PB 306.190 +0,005
PC 379.915 +0,008
of A, Band C, yielded the data in Table 8.28. What are the co-ordinates
of P?
Answer 1219.329 mE, 447.970 mN
Table 8.30
[=
(m)
1797.47
1606.21
1594.23
1735.00
Table 8.31
Determine the best values of the levels at all points B, C, D, etc., taking
that of A to be zero and all values above to be equally reliable.
[London]
Answer level of D = 22.30 m
Table 8.32
[Salford]
Answer B 4.727 m; C 6.309 m; D 9.848 m
Table 8.33
AtoB + 2.275 3
BtoC —1.216 2
CtoD —4.342 1
DtoB +5.508 1
DtoA + 3.263 1
[Salford]
Answer B 2.273 m; C 1.066 m; D —3.258 m
Introduction The computer programs have been written in standard BASIC which should
be suitable for most home computers. The exception to this is the omission
of LET from the assign statements, for example LET X=X +1 is listed as
X=X+1. With reference to the manual for your computer it is possible to
enhance the programs, especially in the area of data input and output.
Angles All angles are input and output in degrees, minutes, seconds. The program
logic uses the angle as one integer value in seconds, this minimizes the round-
ing error that occurs if degrees are used. BASIC’s trigonometric functions
process the angle in radians and the conversion factor
T a 1
(180 * 3600) 206 264.8
is used in the programs.
Formulae Assign statements in the program listing are laid out in the same format as
the formulae in the text, this can lead to unnecessary brackets and hence
inefficient programming. The experienced computer user can remove
unnecessary brackets when entering the program into the computer.
Arc sine Standard BASIC has no function to determine sin~! of an angle, this is a
problem when using the sine rule. It can be overcome using two trigonometric
rules
and tan) ¢) = ae
cos 0 ’
sin 0
so 6 = tan~! |———_ |,
Famer 0)
where tan~! can be represented by BASIC’s function ATN.
Rounding Output data is rounded using the integer routine X=INT ((X * 100) +0.5)/100
to give two decimal places. Most computers have inbuilt functions like ROUND
or PRINT USING which will perform this job automatically.
IF/THEN statements The format of IF/THEN GOTO statements varies widely from machine to
machine, some allow other commands to be nested on the same line often with
an ELSE statement. In the listings a basic form of IF ..., THEN GOTO
x has been used.
Logic functions AND Some of the programs use the logic functions AND and OR in IF/THEN GOTO
or OR statements. Some computers do not support these functions and alternative
program statements are given with the introductions to the appropriate
programs.
Data checking It is good programming practice to check the data when it is entered for non-
sensible or illegal values. The programs can be enhanced by the inclusion of
data checking subroutines such as the example below for checking angles.
INDEX 345
probability distribution, 288, sight rails, 271—4 closed link, 87—90, 213—6
296—9 Simpson’s rule, 143, 146—8 _ closed loop, 76—85
profile boards, 274 sine rule, 104—5 Transit rule, 83—5
proportional frequency, 288 spherical triangular prism, 146, 167, 173
proportionality factor, 124, excess, 108 triangulation, 103—6
128—130 triangle, 94, 107—9, 110 trigonometrical levelling, 1, 2,
primary triangle, 102, 107—9 spheroid, 13 14-22
Pythagoras theorem, 41 spot height, 145, 167—9 trilateration, 103, 243, 327—332
stadia line, 53 true mean, 288
radial acceleration, 193—4 standard deviation, 288—290, trunnion axis, 70, 93—5
radius of curvature, 116 291-9 tunnelling, 243, 262—270
range, 288, 291 standard error, 50, 52, 57—8, turning point, 123
reference 289, 291—303, 323-331 two level section, 144, 154—6
axis, 71 subtense system, 36, 61—2, two-peg test, 8-11
spheroid, 115 257-8, 302-3
tape, 32 super-elevation, 193 variation of co-ordinates, 323—S,
refractive index, 22—4, 35, 48—50 325-331
regression analysis, 305 tacheometry vector misclosure, 85—6
relative precision, 52, 292—5 errors arising in, 57, 59 versine, 227—9
resection, 242, 255—261 stadia, 35, 53—60 vertical
residual, 304, 307, 324, 329 subtense, 36, 61—2, 257-8, _ axis, 70, 93-5
reticule, 1 302-3 curves, 195—6, 229—234
reversal method, 122—4 tangential angle, 191 volumes, 144
rise and fall method, 2—6 tangent point, 190 by spot heights, 167—9
rotating laser, 270—1 tape corrections, 32—4, 36—43 curved section, 163—7
used vertically, 43—5 end-areas rule, 145, 155, 158,
satellites (GPS), 132 “telescope, 1, 53 160
satellite station, 105—7 theodolite, 70 >of cut, 159
scale, 146 maladjustment of, 93—5 of fill, 154, 159
Schuler mean, 122—4 three level section, 157—8 of truncated cone, 170—2
setting out three point problem, 104—5 prismoidal rule, 145, 156, 158,
building, 243-6 Tienstra’s coefficients, 259-261 161
by bearing and distance, 243—6 transit method, 124-131
circular curve, 195—205 Transit rule, 83—5 weight, 290, 295—6, 299-302,
compound curve, 206—213, transition curve connecting 308-9, 311, 314-7,
221-5 circular curves, 216—8 327-331
side widths, 274—7 Transverse Mercator projection, weighted mean, 295—6
sight rails, 271—4 103, 115 Weisbach triangle, 243, 262—270
transition curves, 206—213 trapezoidal rule, 143, 146-8 WGS84 co-ordinates, 132
tunnel, 243, 262—270 traveller, 272—4 . whole circle bearing, 72—9, 86,
vertical curve, 229—234 traverse, 71, 109, 213 88, 243-6
sexagesimal system, 75—6 angular error in, 72—6
shift, 193, 206 area of, 144, 148-151 Young’s modulus, 32
shrinkage of plan, 146 Bowditch’s method, 76—83,
side widths, 144, 274-7 87—90 zenith distance, 54
a
a.
_
2:
hs ale ;
ade (Tete o at
+ Mitte © » Sie
Met ISe*9 Hp
ete
we
iccPet
“as
CA s-
2
SrolNalare Problems in: Surveying
A. BANNISTER & R. RYAN Gad NYelerolate Edition.
Solving Problems in Surveying aims to help prepare students of surveying for their
examinations, by presenting a clear explanation of theory and application in the form of solutions
to typical examination and assignment type questions. Each chapter comprises a summary of
basic theory, a range of worked examples and a selection of problems, with answers, to
encourage further practice and consolidate understanding. Simple computer programs in BASIC
covering common topics encountered by land surveyors are also included. .
This second edition has been updated to reflect advances in modern instrumentation and
techniques, such as GPS, construction lasers, deformation surveys and trilateration.
Coverage includes:
° levelling * areas and volumes
¢ distance measurements ° curve ranging
e theodolite and transverse surveys ¢ setting out and point location
¢ triangulation and the National Grid * observations and adjustments
Arthur Bannister and Raymond Baker are authors of the successful text book Surveying, also
published by Longman, and now in its sixth edition. — }
Raymond Baker is Lecturer in land surveying at the University of Salford, D: === of Civil
Engineering and: Construction. )
A full list of titles in the series appears at the front of this book.
a —— Saas
ISBN 0-S58e@-23b44-4
E33LONGMAN in
780582