Seepage Analysis by FEM (20 SEP 2011)
Seepage Analysis by FEM (20 SEP 2011)
Seepage Analysis by FEM (20 SEP 2011)
Aug 2010
CE5101 Lecture 4
Seepage and FEM
by
Prof Harryy Tan
SEP 2011
1
Outline
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Seepage Analysis
Simple Flow nets
Laplace
L l
E
Equation
ti
q x k x
x
p
Darcys Law
Groundwater Head or Potential
2
q
k x 2 0
x
x
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Equipotentials are
perpendicular to slope;so
piezometer will only rise by
hcosb
tan '
u
F 1
2
tan
cos
F 1
2
cos
tan
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
F 1 w
z cos 2 tan
F 1 w
z cos 2 tan
F 1 w
z cos 2 tan
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
F
z
tan
(With c=0)
h
tan c 1 w tan '
z
10
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Seepage in FEM
2D Formulation in FEM
Material Model and Darcy Law
Validation with Standard Problems
Application to SICC slope failure
Application to excavation
11
12
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
13
Log_10 scale
Why do we need a
permeability function?
Can the problem be
solved
l d without
ith t ititerations?
ti
?
Normal scale
14
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
TRANSITION SATURATED/UNSATURATED
qy K r k y
y
qx K r k x
Kr 1
saturated zone
K 10
r
unsaturated zone
K r 104 h hk
log( K r )
4h
hk
hk 0.7m (PLAXIS)
15
Confined flow
Domain defined
Unconfined flow
16
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
PERMEABILITY
log(k/ko)
Flow vector
Permeability Matrix
18
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
19
Interfaces
ON means Seepage Cutoff
OFF means Seepage allowed through Interface
Drains Zero pore pressure condition
Wells Prescribed flow condition; Inflow ((Recharge)
g ) or
Outflow (Discharge-Well Pumps) Q
Boundary Conditions
Prescribed Heads
Closed BC No Flow Allowed
20
10
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
21
22
11
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
PLAXIS Results
23
24
12
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
H=15m
H=13m
25
Groundwater
Head
H=15m
H 13
H=13m
26
13
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
27
28
14
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
29
30
15
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
31
32
16
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
33
34
17
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
35
36
18
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
37
Elevation (m)
1.714
10
20
30
40
50
60
Distance (m)
38
19
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
1.022
144
142
140
138
Elevation (m)
136
134
132
130
128
126
124
122
120
118
116
114
112
110
0
10
20
30
40
50
60
Distance (m)
39
40
20
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
41
42
21
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
43
44
22
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
5m CUT Undra...
1.5
1.3
1.2
1.1
Displacement [m]
45
46
23
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
47
5m CUT Undra...
1.5
1.3
1.2
1.1
1
0
Displacement [m]
24
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Modeling of
Ground Water in Excavation
Analysis
49
25
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Possible GWT
Conditions in
Excavations
2ba
b
w
uC
2b a
uG
(2b c)a
w
2b c a
51
h=Hb(const)
Ha
Hb
26
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
h H (
h=Ha(const)
)
h=Hb(const)
Ha
Hb
h=Ha(const)
Ha
h=Hb(const)
Suppress uplift
pressures
Hb
27
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
55
56
28
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
57
58
29
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
8m
110.000
GW(S)18
105.000
GW(S)19
16m
GW(S)17
100.000
10-Jan-06
10-Feb-06
10-Dec-05
10-Oct-05
10-Nov-05
10-Sep-05
10-Jul-05
10-Aug-05
10-Jun-05
10-Apr-05
10-May-05
10-Mar-05
10-Jan-05
10-Feb-05
10-Dec-04
10-Oct-04
10-Nov-04
10-Jul-04
10-Jun-04
10-Apr-04
10-May-04
10-Mar-04
10-Jan-04
10-Feb-04
10-Dec-03
10-Oct-03
10-Nov-03
90.000
95.000
10-Aug-04
115.000
59
Date
20
80
100
120
0.00
5.00
10.00
15.00
Depth (m)
undrained analysis
50 mm
drained analysis
97 mm
actual 85 mm
20.00
25.00
Drained Analysis
30.00
Drained
Undrained
I19
35.00
40.00
60
30
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
cv t
kE
h w
andd cv oed so, t
2
h
w
kEoed
Stage
1
2
3
4
5
Date
15-Nov-03
27-Nov-03
15-Jan-04
16-Feb-04
29-Mar-04
12 J l 04
12-Jul-04
7
8
13-Sep-04
18-Dec-04
9 & 10 18-Dec-04
11
12
13
14
3-Mar-05
26-Apr-05
1-Jun-05
1-Jun-05
125 T anchor
75 T anchor
150 T anchor
Construction Activity
Install 1.8m diameter CBP wall GL at RL117.0m
Excavate trench toRL115.0m to cast capping beam
Install Raker Anchor with 80% of 150T preload
Excavate to RL110.5m and install 1st row drain at RL112.5m
Excavate to RL102.5m and
Install drains at RL108.5, 106.5 and 104.5m
F S
For
Sectt 1 andd 22, excavate
t tto RL 96.5m
96 5
Install drains at RL102.5m and 75 T anchors at RL96.5m
For Sect 3, excavate to RL98.5m
Install drains at RL102.5m and 100T anchors at RL100.5m
Excavate to formation level at RL95.9m
Cut small rock berms to RL86.0m; gunnite exposed rock slope
Cast basement wall to RL95.9m and CD slab at RL86.0m
For Sect 1 and 2, remove 75T anchors
Cast basement wall to RL102.5m and slab at RL98.05m
For Sect 3, remove 100T anchors
Cast basement wall and slab at RL105.0m
Cast basement wall and slab at RL115.0m
Excavate to capping beam and remove raker anchors
Backfill to GL at RL117.0m
62
31
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Measured
112.5
Reduce
ed Level (m)
Preloading force
multiplier of 1.4 best
reflects the actual
deflected shape
117.5
102.5
97.5
92.5
87.5
82.5
77.5
0.00
50.00
100.00
150.00
Deflection (mm)
63
4 drains
2 drains
6 drains
64
32
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
117.5m
4 drains
103m
2 drains
6 drains
108m
100m
65
When no drains
installed, max. wall
deflection is greater
300mm
Collapse of wall
Wall deformation
increase with level of
drains which determine
height of water level
behind the wall
102.5
no drains
6 drains
97.5
4 drains
2 drains
92.5
87.5
82 5
82.5
77.5
0
100
200
300
400
Deflection (mm)
66
33
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
GWT
67
GWT
Wall deflect = 83 mm
68
34
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
CBP Elastic,
Elastic
Failure with no
Plastic Hinge,
FOS=1.75
CBP Elasto
ElastoPlastic Failure
with Plastic
Hinge, FOS=1.40
El ti DWall
Elastic
DW ll FOS=1.75
FOS 1 75
Plastic DWall FOS=1.40
70
35
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
117.500
112.500
112.500
107.500
107.500
Section 1 - Stage 5
(after excavate to RL102.5m and installation of
first 2 drains)
117.500
112.500
97.500
Calculated
92.500
102.500
Measured
97.500
Calculated
92.500
87.500
87.500
82.500
82.500
Measured
Reduced Lev
vel (m)
107.500
102.500
102.500
Measured
97.500
Calculated
92.500
87.500
77.500
-50.00
0.00
77.500
-50.00
0.00
50.00
82.500
50.00
100.00
150.00
77.500
0.00
200.00
50.00
Deflection (mm)
Deflection (mm)
150.00
200.00
S e c tio n 1 - S ta g e 7
(a fte r c u t b e rm )
Section 1 - Stage 6
(after excavate to berm top and installing of last
2 drains and anchors)
100.00
Deflection (mm)
1 1 7.5
117.5
117.5
1 1 2.5
112.5
112.5
1 0 7.5
107.5
107.5
Calculated
Me a s u re d
9 7.5
C a lc u l a te d
9 2.5
Measured
97.5
1 0 2.5
102.5
102.5
Measured
97.5
Calculated
92.5
92.5
8 7.5
87.5
87.5
8 2.5
82.5
82.5
7 7.5
0.00
77.5
0.00
50.00
100.00
Deflection (m m)
150.00
200.00
5 0.00
1 0 0.00
1 50 .0 0
D e fle c tio n (m m )
20 0 .0 0
77.5
0
50
100
150
200
Deflection (mm)
71
72
36
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
GWT is
i nearly
l
proportional to
excavation depth
73
Lower 5.6m
Lower 5.6m
37
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
75
H=-12.0m
Closed
Closed
Run SS Seepage analysis with Head settings and closed
flow BCs as in the figure
76
38
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
cluster
DRY
H=-12.0m
Interpolate
Run Phreatic level analysis with General Phreatic level, Cluster Dry in
Excavated Zone, Interpolate below excavated zone, results as above
77
Excavate 5m in 30 days.
Excavate 5m, k=1e-9 m/s
Lower 0.3m
39
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Lower 3.0m
Lower 1.8m
40
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
81
Governing Equations
Steady-state
continuity condition
82
41
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Governing Equations
84
42
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
85
86
43
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
AEV defines the suction value that must be exceeded before air enters the soil pore
Clays have very high AEV compared to Sands
ga is inversely related to AEV
87
88
44
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
89
90
45
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
91
Boundary Conditions
1
1. Water Table
Pw hp w
h y
Pw
2. Inflow
qx nx q y n y qexternal
4. Close boundary
q x nx q y n y 0
3. Outflow
qx nx q y n y qexternal
92
46
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Boundary Conditions
5. Prescribed heads
h h1 , h h2
6. Well/Drain
Q Q
6
7. Free Seepage
h y
8. Screen
q x nx q y n y 0
93
Boundary Conditions
Precipitation
h y hmax
if
Ponding
q x nx q y n y qrain
i
if
h y hmin
if
No infiltration
94
47
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
Boundary Conditions
95
Boundary Conditions
Eg Zone A and B
Eg. Zone C
96
48
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
h(t)
B
9
h(t)
h(t)
h(t)
97
F
From
H=25m
H 25 tto H
H=5m
5 iin 50 d
days
H=25m
H=5m
98
49
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
F
From
H=25m
H 25 tto H
H=5m
5 iin 50 d
days
H=25m
H=5m
99
y
From H=25m to H=5m in 50 days
Potential Slip Surface by c/phi reduction for the
Case of Slow DD in 50 days
100
50
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
WL at 25m FOS=1.74
WL at 5m Very Slow DD FOS=1.63
1.6
1.2
0.2
0.4
0.6
0.8
|U| [m]
101
102
51
AP Harry Tan
CE5101 Seepage FEM
Aug 2010
103
52