Professional Documents
Culture Documents
2020-09-25 Prof Harry Woh Hup Lecture On 25 Sep 2020
2020-09-25 Prof Harry Woh Hup Lecture On 25 Sep 2020
Lecture Outline
• Background to lecture
• Part 1: Use of Practical Advanced Soil Models eg Hardening
Soil (HS) with Small Strain Stiffness (HSS) and Generalized HS
(GHS)
• Three Real Cases and Measured Results
• Part 2: Large Deep Shaft at Changi Water Reclamation
Plant in DTSS1
• Part 3: Re-visit Nicoll Highway Collapse tunnel at Zone M3
• Part 4: Hybrid Raft Foundation for Wisteria Condo at
Yishun Avenue 4
• Summary 2
1
9/24/2020
4
HS, HSS & GHS Practical Advanced Models
1. Introduction
2
9/24/2020
1. Introduction 5
General soil behavior
𝑑𝜀 𝑝 𝑑𝜀 𝑒
3
9/24/2020
9
PMT
8 FEM
7
6
P' (kg/cm2)
5
4
3
2
1
0
3.95 4 4.05 4.1 4.15 4.2 4.25
R0 (cm)
εyy
0.04
settlement troughs behind ERSS or above tunneling
0.05
works. 0.06
7
▪ FE model can be boundary independent 0.07
• Observe unrealistic
heaving in MC models
4
9/24/2020
10
5
9/24/2020
11
Generalized
Hardening Soil Model
• User defined soil model
• Alternative stress dependency formulation
• Modular version of the Hardening Soil model
11
11
12
12
6
9/24/2020
13
Before
After
excavation
excavation
14
7
9/24/2020
Parameters Meaning
1 E50ref Secant stiffness in standard drained triaxial test
2 Eoedref Tangent stiffness for primary oedometer loading
Deformation
3 Eurref Unloading / reloading stiffness
4 m Power for stress-level dependency of stiffness
5 c'ref (Effective) cohesion
6 φ (Effective) angle of internal friction Relevant tests:
Strength
7 ψ Angle of dilatancy • Drained triaxial test
8 σt Tension cut-off and tensile strength • 1D consolidation test
9 νur Poisson's ratio for unloading-reloading
10 pref Reference stress for stiffnesses • Pressuremeter test
Advanced
11 K0 nc K0 value for normal consolidation (Default=1-sinφ)
• CPTu test
12 Rf Failure ratio
13 G0ref Reference shear modulus at very small strains
Small strain • Resonant column / Bender element
14 γ0.7 Threshold shear strain at which Gs = 0.722G0
• Seismic CPT
15
Soil stress history is input by OCR or POP for Ko Initial Conditions • Geophysical surveys
15
16
16
8
9/24/2020
Interface elements
Prescribed displacement
Or Line load
▪ Full (depth) numerical model allows for a direct input of soil stress
history
17
17
1
2
3
▪ The FEM simulation will be able to capture the pressuremeter’s primary loading,
unloading/reloading behavior once initial contact stresses with BH have been established.
18
18
9
9/24/2020
19
4. Conclusion 20
Some Concluding remarks
20
20
10
9/24/2020
2. Pore pressure measured at IP2 (under base center of shaft) is key data to understand soil
unloading behavior history
3. Study Consolidation models using Plaxis Classical Consolidation analysis (Excess EPP), and
more advanced Fully Coupled Flow-Deformation Consolidation analysis (Total TPP)
4. Compare Consolidation predictions as bounded by UnDrained and Drained analysis, using the
same sets of Effective Stress parameters
21
21
Deep Shafts
22
22
11
9/24/2020
BACKGROUND
Excavation from Feb to Oct 2002 for IPS2 Shaft
Dia 42.6m
Depth 70m
Piezo PZ-2 (2m below FML)
Dia 42.6m
Depth 70m
Dia 35.8m
Depth 69.5m
23
23
24
12
9/24/2020
25
Soil Permeability
Function:
Kr vs 𝜓
26
13
9/24/2020
UnDrained is
most stable
27
27
PLAXIS CE 2020
PLAXIS COUPLED FLOW-
63 m depth 70 m depth PLAXIS CONSOLIDATION DEFORMATION TPP CONSD AT
DRAINED AT RL39.3m (failed) UNDRAINED AT RL33.0m(FML) EPP CONSD AT RL33.0m(FML) RL33.0m (FML)
28
28
14
9/24/2020
29
29
• UnDrained Model gave unrealistic high excess PP soil suctions below the base (+500 kPa)
• Classical EPP Model gave larger excess PP suctions compared to Coupled TPP model below FEL level (+150 kPa)
30
• Coupled TPP Model gave most realistic estimates of remaining excess PP suction below FEL level (+120 kPa)
30
15
9/24/2020
• COUPLED TPP
CONSD GIVE
REALISTIC LOWEST
ESTIMATE OF BASE
FOS
• BUT ALL THREE
HAVE HIGH FOS
AGAINST BASE
UPLIFT FAILURE
31
31
MEASURED PZ-2 is
yellow dash line
DRAINED
FAILURE
INITIATED
Fully Coupled Consolidation Model gave most realistic response of PZ2 predictions – YELLOW DASH LINE
32
32
16
9/24/2020
ACTIVE PP IN OA at 2m below
FML shifted to match Total Mean
Stress at start of Excavation
33
34
34
17
9/24/2020
90 90 90 90
80 80
80 80
70 70
70 70
Reduced Level (m)
50 50
50 50
40 40
40 40
30 30
30 30
20 Measured at I1 Measured at I1
20
Measured at I3 20 Measured at S4 Measured at S4
Measured at I3 20
FEM, With T emperature FEM, With T emperature Measured at S12 Measured at S12
10 10 Measured at S20
FEM. With T emperature (Adjusted) FEM. With T emperature (Adjusted) Measured at S20
10 10
FEM, With T emperature FEM, With T emperature
0 0
0 0
35
35
36
18
9/24/2020
37
37
38
38
19
9/24/2020
39
39
40
40
20
9/24/2020
Importance of Correct
UnDrained Shear
Strengths in Soft Clays
to assess Short-term
Safety of Cut slopes and
excavations
41
41
• UnDrained strengths is
stress-path dependent
• Different soil models will
produce different effective
stress paths
• UnDrained “B” is strictly
applicable to Mohr-
Coulomb model, and not
applicable to most
Advanced soil models
• For Advanced soil models
need to check UnDrained
Strengths profiles against
measured values for each
stage of excavations
42
42
21
9/24/2020
43
43
44
22
9/24/2020
45
45
FULL FEM MODEL SHOWED THAT WEAKER SIDE IS SOUTH SIDE (BH31 OR BH84)
46
46
23
9/24/2020
REDO WITH CE-2020 YIELDED STATE AT RL72.5m (After removal of Upper JGP)
47
48
48
24
9/24/2020
• Is Linear trend
extrapolation of PP
in OA soils below
JGP valid and
reasonable?
49
49
Ideal UnDrained
MEASURED PP AT GW-v24 IN
OA COMPARE WITH PLAXIS
COUPLED FLOW RESULTS
50
25
9/24/2020
51
51
UNDRAINED INDICATED BY
• Based on mean effective stress, OA
Sig_P’ is CONSTANT FOR 180 DAYS
below FML is UnDrained for about
P’ Effective
180 days (up to about 15m depths)
Mean Stress
and partially drained up to end of
280 days at about 30m depths
52
52
26
9/24/2020
PREDICTION OF PWP AT GW-v25-1 (-12,72.5) in Lower Marine Clay (LMC) using FULLY
COUPLED FLOW
NHW MEASURED PWP
MEASURED PP AT GW-v25-1 IN
LMC COMPARE WITH PLAXIS
COUPLED FLOW RESULTS
53
UNDRAINED
Sig_P’= CONSTANT (280 days) P’ Effective Mean Stress
• Based on Mean Effective Stresses, LMC is
essentially UnDrained throughout the
excavation process, as we have correctly
assumed in our typical design practice.
54
54
27
9/24/2020
PREDICTION OF PWP at GW-v25-2 (-12,85.5) in Upper Marine Clay (UMC) using FULLY
COUPLED FLOW
MEASURED PP AT GW-v25-2
IN UMC COMPARE WITH
PLAXIS COUPLED FLOW
RESULTS
55
56
56
28
9/24/2020
57
57
58
58
29
9/24/2020
59
59
• 3 tower blocks
• Core walls
• 8.4m x 8.4m Columns
• Podium
• 8.4m x 8.4m Columns
60
60
30
9/24/2020
61
61
Mixed Foundation
• Raft Foundation
• 1.7m thick Raft on
• 3 x rigid inclusion
zones (blue)
• Base slab overhang
• Bored Pile Zones (red)
• 600mm to
1000mm diameter
62
62
31
9/24/2020
• THUS, THIS GROUND CONDITION WILL ALLOW FOR RAFT FOUNDATION (WITH PILE
INCLUSIONS USED AS GI UNDER HEAVY TOWERS, TO ACT AS SETTLEMENT REDUCERS)
63
63
64
64
32
9/24/2020
65
65
Extensive SI work has been conducted in July 2015 with 11 boreholes PBH-1 to PBH-11, coupled with
detail lab Triaxial CID tests and site pressuremeter tests to ascertain soil strength parameters and
stiffness parameters for HS model. There are also 4 preliminary borehole BH-1 to BH-4 conducted in
Aug 2014 shown below:
PBH-6 PBH-5
PBH-9 PBH-8
PBH-7
PBH-2
PBH-11
BH-2
PBH-10 BH-3
PBH-3
66
66
33
9/24/2020
• The predominant geological feature of the site is the Bukit Timah Formation
with typically about 3m top fill material followed by the competent Bukit
Timah residual soil with SPT N increasing with depth.
• It is noted that ALL the boreholes consistently show that there exists NO
trace of Kallang formation soft soils at site which is very favorable factor
for Hybrid Raft foundation.
• As will be demonstrated in the next slides, the Raft foundation with base
resting on an elevation of about 6.6m below GL, the founding ground
condition is predominantly stiff Residual Soil with typical SPT N value in the
range of 10~20 with average SPT N value of about 15 which is deemed to
be sufficiently strong to support the Raft foundation with tolerable
settlements that can be managed upon application of building loads.
67
67
68
68
34
9/24/2020
PBH1 – Top fill 3m followed by competent Bukit Timah Residual Soil with SPT N
increasing steadily with depth
69
69
PBH5 – Top fill 3m followed by competent Bukit Timah Residual Soil with SPT N
increasing steadily with depth
70
70
35
9/24/2020
PBH10 – Top fill 3m followed by competent Bukit Timah Residual Soil with SPT N
increasing steadily with depth
71
71
BH-1 – Top fill 1m followed by competent Bukit Timah Residual Soil with SPT N
increasing steadily with depth
72
72
36
9/24/2020
BH-4 – Top fill 1m followed by competent Bukit Timah Residual Soil with SPT N
increasing steadily with depth
73
73
Just for information: The underlying rock elevation profile (below existing
GL) at site was shown below which is observed to be quite variable, from
shallow elevation of about -12 ~ -14m bgl around BH-1 and PBH-9 to
about -32~-34m bgl around BH-2 and PBH-10
74
74
37
9/24/2020
Two rock core boxes shown below as an example to illustrate very good
quality of the underlying granite rocks with typically high TCR and RQD
75
75
76
76
38
9/24/2020
Deployed boreholes
Top fill
Firm RS (N=7~10)
Stiff RS (N=10~20)
Very stiff RS
(N=20~40)
Hard RS (N=40~80)
Underlying rocks
77
77
• Hardening Soil Model (HS) was adopted for analysis, as it is very suitable for
construction involving intermediate unloading (1-level basement excavation) and
subsequently column/wall loading.
• Pressuremeter tests was conducted at site with results summarized below. It can be seen
that the 1st time loading modulus Em typically hovers around 1 times SPT N value (in
MPa) except for PBH8 which is 3.5 times SPT N value. It is noted that OYO Elastmeter
was lowered into “Pre-drilled” borehole instead of “Self-bored” method, which could
explain certain soil disturbance and relatively lower 1st time loading modulus obtained.
• On the other hand, Eur/E50 was observed to be very high in the range of 7.4~12.0 as
tabulated below:
Depth of
pressurement Test 1st loading modulus, Unload-reload
Borehole (m) SPT N value Em (kPa) Em vs SPT N modulus, Er (kPa) Er / Em
PBH7 20~20.5m 79 79,494 kPa 1.01N (in MPa) 700,462 kPa 8.8
PBH8 22~22.5m 100 351,124 kPa 3.51N (in MPa) 2,594,818 kPa 7.4
PBH9 10.5~11m 40 25,368 kPa .63N (in MPa) 222,214 kPa 8.8
PBH10 25~25.5m 64 51,974 kPa .81N (in MPa) 624,937 kPa 12.0
PBH11 14~14.5m 17 19,203 kPa 1.13N (in MPa) 184,764 kPa 9.6
78
78
39
9/24/2020
• Multiple lab Triaxial CID tests have been conducted with the soil strength
parameters summarized below. It is noted that generally very high soil cohesion
was derived with correspondingly lower friction angles.
• It is assessed that for such Bukit Timah residual soil with typically sandy silt
material, a lower cohesion coupled with correspondingly higher friction angle
would be more reasonable as recommended in the last column of the summary
table below, and was adopted in the present geotechnical 3D FEM analysis.
Sample Sample SPT N Lab cohesion Lab friction Recommended Recommended friction
Category Borehole Sample Depth (m) value c' (kPa) angle f ' () c' (kPa) angle f ' ()
Firm Soil
(N<10) PBH11 TW1 3~3.8m 9 11 kPa 32 degree 5 kPa 30 degree
79
Determination of Soil Stiffness Parameters for Bt Timah Residual Soils at Wisteria Site
Borehole Depth m SPT N Em MPA Em/N Eur MPA Eur/N Engineering Judgment on Parameters to
PBH7 20-22.5 79 79.5 1.01 700.5 8.87
PBH8 22-22.5 100 351.1 3.51 2594.8 25.95
Use:
PBH9 10.5-11 40 25.4 0.64 222.2 5.56 • Correlations to SPT N values
PBH10 22-22.5 64 52 0.81 624.9 9.76
PBH11 14-14.5 17 19.2 1.13 184.8 10.87 • E50_ref range from 0.64N to 3.51N
MEAN= 1.42 12.20 • More realistic range from 0.64 to 1.5N
STDDEV= 1.18 7.94
CHAR VALUE 0.83 8.23 • Characteristic value is 0.79 to 0.83N
Borehole Depth m SPT N Em MPA Em/N Eur MPA Eur/N
• We use E_50ref=1N MPA, with N as
PBH7 20-22.5 79 79.5 1.01 700.5 8.87 average values for each distinct soil
PBH9 10.5-11 40 25.4 0.64 222.2 5.56
PBH10 22-22.5 64 52 0.81 624.9 9.76 layers
PBH11 14-14.5 17 19.2 1.13 184.8 10.87
MEAN= 0.90 8.76
STDDEV= 0.22 2.29 • Eur_ref range from 5.56 to 25.95N
CHAR VALUE 0.79 7.62
• More realistic is 5.56 to 10.87N
• Characteristic value is 7.62 to 8.23N
• We use Eur_ref=5N MPA, with N as
average values for each distinct soil
layers 80
80
40
9/24/2020
81
81
82
82
41
9/24/2020
EC7 (annex H) explicitly states that for normal structures with isolated
foundations, total settlements up to 50mm are often acceptable; while
maximum relative rotation is 1/500 (see below)
83
Column loads simulated as “point” loads in 3D FEM mesh (some top soil
elements hidden to display loadings below GL)
84
84
42
9/24/2020
Wall loads simulated as “line” loads in 3D FEM mesh (some top soil elements
hidden to display loadings below GL)
Line loads
on walls
85
85
There are also some patch surface loads for backfills (94kPa) and the Detention
tanks/Sprinkler tank (50kpa) and some local areas (16kPa) as well as general
basement surcharge of DL+SLL = 10.1kPa
Some patch
surface loads for General basement surcharge
backfill, DL+SLL = 10.1kPa
Detention
tanks/sprinkler
tank and some
local areas
(94/50/16kPa)
86
86
43
9/24/2020
Loads on Raft
Loads on Raft within basement
outside basement
Loads on bored
piles
87
87
The calculated Final raft settlement contour is shown below with Max
settlement of about 42mm around the 3 Lift corewall areas, which exceeds
the Max allowable settlement of 40mm, thus the Original proposed Raft
scheme has very tight Serviceability Issue.
88
88
44
9/24/2020
89
89
90
90
45
9/24/2020
91
91
GI work using H-piles as Settlement Reducers (SR) to embed into the hardcore
layer by 200mm
1.7m Raft
92
H-piles driving or Jackin until reaching REFUSAL of underlying hard soil with SPT N=60~100
92
46
9/24/2020
GI work by H-piles
(400x400x232kg/m) in
2m x 2m grid pattern Raft edge outline
to cover the footprint
of the heavily loaded
lift corewell area #1
Local soil profile (PBH-9):
0~3m below GL: Fill
3~9m below GL: Stiff residual
soil (N=13~18)
9~13m below GL: Very stiff
residual soil (N=40~41)
13m downwards: Underlying
rocks
(Estimated H-pile toe elevation = 93
-13m below existing GL)
93
94
94
47
9/24/2020
95
Adjacent piled
foundations
outside basement
96
48
9/24/2020
97
98
98
49
9/24/2020
With the above Settlement Reducer (SR) piles, Max settlement of the raft was
effectively reduced to about 25.1mm as shown below, much less than the 40mm
allowable settlement → OK!
99
99
For the continuous Raft foundation, the differential settlement can be visualized by plotting a cross
section across the Max settlement trough as marked in Section 1-1 below and presented in next slide.
100
100
50
9/24/2020
-10
-20
Settlement (mm)
-30
-40
-50
-60
-70
Along Max gradient segment (green line):
-80 Max differential settlement = 21.8mm - 11.1mm = 10.1mm
Over the span of 15.9m - 2.8m = 13.1m
-90 So, Max differential distortion ratio = 10.1mm/13,100mm = 1: 1290 -- > OK!
-100
101
101
Likewise, for the settlement profile along the cross section 2-2 as marked below and
presented in next slide.
102
102
51
9/24/2020
0 10 20 30 40 50 60 70 80 90
0
-10
-20
-30
Settlement (mm)
-40
-50
-60
-100
Along 2-2
103
103
104
104
52
9/24/2020
105
105
Final measured Max Raft settlement of about 20~23mm around the 3 BLOCKS areas as
illustrated below, matching the predicted Max raft settlement of 20~25mm very well:
106
106
53
9/24/2020
107
107
On the other hand, the measured PODIUM settlement is much less of around
2~10mm, which also matches the 3D FEM contour pattern reasonably well as
illustrated in next slide:
108
108
54
9/24/2020
109
109
110
110
55
9/24/2020
Geotechnical 3D FEM prediction of 6.5mm bored pile settlement (see below), compared
reasonably well with the final settlement marker BSM-21 measurement of 6.0mm:
4W3 on bored
piles
4W2 on Raft
Differential Column
Column settlement settlement spacing Differential
Columns (mm) (mm) (m) settlement
111
112
112
56
9/24/2020
113
113
114
114
57
116
115
Relative movement (mm) Relative movement (mm)
-5
-25
-20
-15
-10
0
5
10
-5
-25
-20
-15
-10
0
5
10
Relative movement (mm) Relative movement (mm)
-5
-25
-20
-15
-10
0
5
10
-25
-20
-15
-10
-5
0
5
10
22-Feb-17 22-Feb-17
22-Feb-17 22-Feb-17
24-Apr-17 24-Apr-17
24-Apr-17 24-Apr-17
24-Jun-17 24-Jun-17
24-Jun-17 24-Jun-17
24-Aug-17 24-Aug-17
BSM 007
BSM 005
24-Aug-17 24-Aug-17
BSM 003
BSM 001
24-Oct-17 24-Oct-17
24-Oct-17 24-Oct-17
24-Dec-17 24-Dec-17
24-Dec-17 24-Dec-17
23-Feb-18 23-Feb-18
23-Feb-18 23-Feb-18
Relative movement (mm) Relative movement (mm)
-5
-25
-20
-15
-10
0
5
10
Relative movement (mm) Relative movement (mm)
-25
-20
-15
-10
-5
0
5
10
-5
-25
-20
-15
-10
0
5
10
-5
-25
-20
-15
-10
0
5
10
22-Feb-17 22-Feb-17
22-Feb-17 22-Feb-17
24-Apr-17 24-Apr-17
24-Apr-17 24-Apr-17
24-Jun-17 24-Jun-17
24-Jun-17 24-Jun-17
24-Aug-17 24-Aug-17
BSM 008
BSM 006
24-Aug-17 24-Aug-17
BSM 004
BSM 002
24-Oct-17 24-Oct-17
24-Oct-17 24-Oct-17
24-Dec-17 24-Dec-17
24-Dec-17 24-Dec-17
23-Feb-18 23-Feb-18
23-Feb-18 23-Feb-18
Comparison: Individual BSM
116
115
9/24/2020
58
118
117
Relative movement (mm) Relative movement (mm) Relative movement (mm) Relative movement (mm)
-25
-20
-15
-10
-5
0
5
10
-5
-25
-20
-15
-10
0
5
10
-5
-5
-25
-20
-15
-10
0
5
-25
-20
-15
-10
0
5
10
10
22-Feb-17 22-Feb-17 22-Feb-17 22-Feb-17
BSM 011
BSM 009
BSM 015
BSM 013
24-Oct-17 24-Oct-17 24-Oct-17 24-Oct-17
-25
-20
-15
-10
-5
0
5
10
-25
-20
-15
-10
-5
0
5
10
-25
-20
-15
-10
-5
0
5
10
22-Feb-17 22-Feb-17 22-Feb-17 22-Feb-17
BSM 016
BSM 014
BSM 012
BSM 010
not clear
adjusted but stage of
*BSM 014 & 015 should be
construction at installation is
118
117
9/24/2020
59
9/24/2020
SUMMARY
• Basis of Calibration of HS, HSS and GHS type models for practical
applications was presented
• Use of these models in Deep Shafts, Cut and Cover tunnel with pore
pressures response data showed good predictions capabilities when used
in Fully Coupled consolidation analysis (Cases of DTSS1 Very Deep Shaft
and NHW C824 Section M3 failure)
• Good measurements of soil stiffnesses with PMT and laboratory CID plus
Oedometer tests, allow for reliable foundation settlements performance
predictions with 3D-FEM analysis in stiff to hard residual soils (Case
Wisteria project in Bt Timah Soils)
119
119
60