Historical Development of Dam Analysis: I.I. Earth and Rockfill Dams
Historical Development of Dam Analysis: I.I. Earth and Rockfill Dams
Historical Development of Dam Analysis: I.I. Earth and Rockfill Dams
Historical development of
dam analysis
Section as designed
Drainage
^^Z^^^y^ channel
Ballast
Fig. 2. Non-circular slip of the downstream slope of Lea River Dam near Chingford,
s
Essex, England, in July 1937
Fig. 6. Results oj model tests and finite clement analyses for the 220 m high
Dabuklamm Dam in the Darter Valley. Eastern Tyrol1''2"
1.2. G r a v i t y d a m s
(373-10) (372-95)
Fig. 8. Cross-section of Bouzey gravity dam near Epinal, France, built from 1877 to
21
1882 and collapsed in 1895
2%
larger t h a n the vertical compressive stress. I n 1898 Lieckfeldt for
m u l a t e d a p r o o f of stresses considering j o i n t w a t e r p r e s s u r e . T h e p r o b
lem of p e r c o l a t i o n considering d a m s as p o r o u s m e d i a w a s
29
t a k e n u p by Fillunger in 1913. H i s theoretical views were vigorously
30
contested by Terzaghi a n d in 1934 he c o m m i t t e d suicide. I n 1938
31
Heinrich showed t h a t the correct solution of the p e r c o l a t i o n p r o b
32
lem m u s t start from a t w o - p h a s e system. In 1936 Brahtz reported
the application of b o t h Airy's stress function a n d photoelasticity
t o the stress analysis of gravity d a m s , considering different elastic
properties for b o t h d a m a n d u n d e r g r o u n d . Different Young's m o d u l i
of elasticity for d a m a n d u n d e r g r o u n d were t a k e n i n t o a c c o u n t in
33
the m a t h e m a t i c a l t h e o r y of functions by Tolke in 1938 (Fig. 10).
T h e e n o r m o u s a m o u n t of c o m p u t a t i o n involved w a s a v o i d e d in 1947
34
by Zienkiewicz w h o solved the system of difference e q u a t i o n s with
35
the help of the m e t h o d of relaxation, devised by Southwell in 1940.
W h e n designing the 285 m high gravity d a m G r a n d e D i x e n c e in
36
Switzerland, Stucky studied the influence of r o c k elasticity with
the help of m o d e l tests as h a d been d o n e for the 221 m high H o o v e r
D a m in the U S A a b o u t 20 years previously. O f p a r t i c u l a r i m p o r t a n c e
to the p r o o f of stability of gravity d a m s is the a s s u m e d d i s t r i b u t i o n of
3 7
w a t e r pressure at the d a m base (Fig. I I ) w h i c h d e p e n d s b o t h o n the
position of the g r o u t c u r t a i n a n d o n the efficiency of the d a m base
3S
d r a i n a g e . In 1965 Rescher r e p o r t e d a detailed analysis of the m a i n
1 ft = 0-305 m
1 25
Fig. 9. Stresses in a gravity dam according to tests with rubber models * and to finite
3
difference analysis *
n=E /E
f b =2-0 =1-0 =0-5 n=E /Ei b =2-0 =1-0 =0-5
Fig. 10. Influence of the ratio of Young's moduli of elasticity of concrete and bedrock
33
on stresses in gravity dams
1 ft = 0-305 m
Fig. 11. Predicted and measured uplift distribution below the 89 m high Hiwassee
31
Dam of the Tennessee Valley Authority, USA
Fig. 12. Influence of the ratio of crest length to dam height on the base stresses of
1
gravity dams subjected to deadload and hydrostatic pressure*
1.3. A r c h d a m s
Analyses of a r c h d a m s were p e r f o r m e d until the 1920s (e.g. in 1925
for the 113 m high P a c o i m a D a m in California) with help of
Mariotte's (1620 - 1686) well k n o w n ring or t u b e formula. In 1889,
43
24 years before the w o r k of Ritter w h o is considered as the orig
44
i n a t o r of the m e t h o d in E u r o p e , Vischer a n d Wagoner described
the a d j u s t m e n t of r a d i a l displacements of a n a r c h d a m by assuming
several h o r i z o n t a l arches a n d only one (crown) cantilever. This
m e t h o d w a s first applied in 1904 t o the analysis of the 67 m high
Pathfinder D a m a n d the 99 m high Buffalo Bill D a m (Fig. 14) in
\e = H/Ebt !'4g /3
0 0
P
l l l l l l l
Fig. 13. Horizontal section of a gravity dam in a narrow valley acting as a tension-free
42
beam with the strain distribution in the pressure line of an equivalent arch
4 5
the U S A . In o r d e r t o utilize better the e m p l o y e d c o n c r e t e in his
design of the 52 m high S a l m o n C r e e k D a m in A l a s k a (Fig. 15) in
46
1 9 1 3 - 1 9 1 4 , Jorgensen following a n earlier p r o p o s a l by
4 7
Eastwood in 1910, d e p a r t e d from the c o n s t a n t r a d i u s type of a r c h
d a m a n d realized the first c o n s t a n t angle d a m . Since t h e n c o n s t a n t
angle d a m s h a v e been the d o m i n a n t d a m s h a p e in n o t t o o wide valleys.
Nevertheless, there h a v e been a few retrogressions: t h e 114 m high
Spitallamm D a m in Switzerland built in 1926 - 1932, the 180 m high
Tignes D a m in F r a n c e in 1949 - 1953 a n d the 245 m high
S a y a n o s h u s h e n s k D a m in Siberia in 1970 - 1988. T h e 55 m high
M o n t s a l v e n s D a m built in Switzerland ( C a n t o n F r i b o u r g ) in
4S
1919 - 1920 w a s analysed by Stucky using a r a d i a l a d j u s t m e n t
of four arches a n d nine cantilevers (Fig. 16). This grid m e t h o d
was extended in the U S A from 1923 with t a n g e n t i a l a n d twist
49
adjustments, a n d p r e s e n t e d in 1929 by Howell a n d Jaquith as
the 'trial l o a d m e t h o d ' , because the division of l o a d s h a d t o be found
50
by successive trials. By 1925 Vogt h a d s h o w n t h a t the displacements
of f o u n d a t i o n rock could be predicted a p p r o x i m a t e l y using three
51
formulae derived from Boussinesq's t h e o r y of the elastic s e m i - s p a c e .
T h e trial load m e t h o d h a d already been checked in 1926 w i t h the help
of the Stevenson C r e e k Test D a m (Fig. 17), a n 18 m high a n d only
5 2
61 cm thick experimental a r c h d a m n e a r F r e s n o in California.
T h e trial load m e t h o d u n d e r w e n t its first a p p l i c a t i o n o n a truly g r a n d
scale in H o u k ' s analysis of the 221 m high H o o v e r D a m (Fig. 18)
5 3
in the U S A . T h e results of the e x t r a o r d i n a r i l y v o l u m i n o u s c o m p u t
ations were checked with a m o d e l test (scale 1 in 240) at the University
54
of C o l o r a d o at Boulder (Fig. 1 9 ) .
1 ft = 0-305 m
Fig. 14. The 99 m high constant radius arch dam Buffalo Bill near Cody, Wyoming
45
USA
re
a
a
A Height in feet (1 ft 0-305 m) above river bed
B Height in feet
C Upstream radius in feet
D Ring stress in psi (1 psi = 0-007 MN/m2)
Fig. 15. The 52 m high constant angle arch dam on the Salmon Creek in Alaska**'
Portion of load
Fig. 16. Radial displacements and load distribution of the 55 m high constant angle
4
arch dam Montsalvens on the Jogne River in Switzerland *
Fig. 17. Radial displacements of the 18-3 m high constant radius arch dam on the
52
Stevenson Creek near Fresno in California
Fig. 18. The 221 m high Hoover Dam on the Colorado River in the Black Canyon
53
between Arizona and Nevada, USA
Fig. 20. Crest displacements of the 90 m high and very thin Tolla Dam in Corsica after
57
in situ measurement and prediction by six adjustments
Fig. 21. Radial, tangential and vertical displacements by shell theory of 130 m high
63
constant angle dam Punt dal Gall in the Livigno Valley, Italy
64
I n 1992 Herzog p r e s e n t e d a n a p p r o x i m a t e m e t h o d for prelimi
n a r y design a n d checking p u r p o s e s . T h e m e t h o d w a s shown to be
entirely a d e q u a t e by the p o s t e r i o r analysis of some 60 a r c h d a m s .
It is b a s e d o n solution of the t a n k e q u a t i o n by the classic
e d g e - p e r t u r b a t i o n m e t h o d a n d yields, with the help of closed form
ulae, a conclusive picture of the stresses a n d displacements of a n a r c h
d a m subjected t o w a t e r pressure a n d t e m p e r a t u r e changes after only
a n h o u r of m a n u a l calculation.
3 9
T h e finite element m e t h o d h a d to wait for d e v e l o p m e n t of the
6 5
m a t r i x displacement m e t h o d before it b e c a m e possible to analyse
69
Fig, 22. Arch dam in fissured bedrock with percolation
6 6 67
a r c h d a m s with the help of i s o p a r a m e t r i c e l e m e n t s . A t present
6 8 6 9
b o t h the f o u n d a t i o n r o c k a n d p e r c o l a t i o n of the d a m (Fig. 22)
can be a c c o u n t e d for. T h e newest d e v e l o p m e n t s consider the o p e n i n g
70
of c o n s t r u c t i o n j o i n t s a n d the f o r m a t i o n of c r a c k s . I n general, analy
7 1
sis of a r c h d a m s using the finite element m e t h o d requires very
n
powerful c o m p u t e r s . H o w e v e r , even as early as 1986 Stevenson et a l
r e p o r t e d the finite element analysis of a n a r c h d a m o n a p e r s o n a l c o m
puter.