Japanese Geotechnical Society Special Publication The 15th Asian Regional Conference on

Soil Mechanics and Geotechnical Engineering

Impacts from three-dimensional effect on the wall deflection induced by a deep excavation in
Kaohsiung, Taiwan

Bin-Chen Benson Hsiung i) and Sy-Dan Dao ii)

i) Associate Professor, Department of Civil Engineering, National Kaohsiung University of Applied Sciences, 415Chien-Kung Road,
Kaohsiung City, 807, Taiwan.
ii) PhD Candidate, Department of Civil Engineering, National Kaohsiung University of Applied Sciences, 415Chien-Kung Road,
Kaohsiung City, 807, Taiwan.

ABSTRACT

In this paper, it aims to examine impacts from three-dimensional effect on the wall deflection induced by a deep
excavation in Kaohsiung city, Taiwan. The commercial software PLAXIS 3D was used as a numerical tool for 3D
finite element analyses in this study. First, a benchmark analysis was performed to simulate a case history of deep
excavation in thick layers of sand (Case A) to verify the performance of 3D numerical analysis model in predicting
the wall displacements. It is aware that there is a little limitation in prediction of the wall movements by using a
constitutive soil model having single elastic modulus. Next, a series of parametric studies that uses the same input
parameters as the benchmark analysis was conducted to model the excavation of Case A with various values of
excavation length (L) and width (B). From these parametric studies, plane strain ratio (PSR), which is the ratio of the
maximum wall deflection of a certain section to the maximum wall deflection of the section under the plane strain
condition, was determined with various values of distance from evaluated section to the excavation corner (d), length
(L) and width (B) of excavation. A relationship between PSR, d and ratio of B/L was thus interpreted. It is
summarized that PSR is smaller than in sand rather than in clay for B/L more than 1.0, but it is larger than in sand
rather than in clay for B/L less than 0.5. Further verification may have to be delivered later in order to explore the
reason for the difference.

Keywords: deep excavation, wall deflection, three-dimensional effect, plane strain ratio, PLAXIS 3D.

1 INTRODUCTION was first proposed, which is the ratio of the maximum

wall deflection of a certain section along the wall to the
The wall deflection induced by a deep excavation
maximum wall deflection of the section under the plane
depends on many factors such as conditions of soil and
strain condition. The PSR values were also determined
groundwater, excavation geometry, surcharge load,
and evaluated for a typical excavation in condition of
existence of adjacent structures, construction method,
clayey soils in Taipei, Taiwan.
stiffness and penetration of the retaining wall, type and
In this paper, the PSR values will be determined for
installation method of struts, spacing and stiffness of
a deep excavation in condition of sandy soils in
struts, and ground improvement. In addition, in a
Kaohsiung, Taiwan. The PSR values will be interpreted
certain excavation, the wall deflection also depends on
using a series of parametric studies with various
position of evaluated section or distance from evaluated
excavation lengths and widths. It is expected that this
section to the excavation corner. Effect of the position
study could be useful to predict the 3D maximum wall
of evaluated section on the wall deflection is called to
deflections from 2D analyses, which is commonly
be the three-dimensional effect or corner effect on the
adopted in engineering practice. The commercial
wall deflection.
software PLAXIS 3D, version 2013, is selected as a
The wall deflections have been studied by many
tool for numerical analyses in this study.
researches such as Clough and O'Rourke (1990), Ou et
al. (1996), Ou et al. (2000), Ou (2006), Kung et al.
2 BENCHMARK ANALYSIS
(2007), Lin et al. (2007), Hsiung (2009), Wang et al.
(2010), Likitlersuang et al. (2013) and Khoiri and Ou In this section, a case history of deep excavation in
(2013). However, most of these researches mainly Kaohsiung, Taiwan, namely Case A, was used as a
focused excavations in clays, not in sands. In the study basic to verify the validity and performance of the
of Ou et al. (1996), a concept of plane strain ratio (PSR) benchmark analysis, which will be applied as a base of

http://doi.org/10.3208/jgssp.TWN-15 1602
later numerical analyses. Case A was located in the dominantly affect the excavation behavior.
central area of Kaohsiung city, next to the O7 station, According to previous researches, such as Hsieh et
which is on the Orange line of Kaohsiung MRT system. al. (2003), Ou (2006), Kung et al. (2007) and
The shape of Case A was rectangular with 70 m in Schweiger (2009), the constitutive soil model adopted
length and 20 m in width. The excavation was carried in numerical analysis has a very little influence on
out by the bottom-up construction method and was predicting the wall deflection induced by deep
retained by a diaphragm wall that is 0.9 m thick and 32 excavations. As a result, the linear elastic-perfectly
m deep. It was excavated by five-stage excavation with plastic Mohr-Coulomb model (MC model), a basic
the maximum excavation depth of 16.8 m. The model of soil, was adopted to simulate soils in the
retaining wall was propped by steel struts at four levels, model of benchmark analysis. The distance from the
and the horizontal spacing of the struts was average lateral boundaries of this model to the retaining wall
about 5.5 m. Fig. 1 below shows the cross section and was taken to be four times excavation depth. This
ground condition of Case A. distance is considered because the maximum length of
influence zone of settlement behind the retaining wall
Ground level: 0.0 m SID was four times excavation depth according to previous
CL, Su = 28MPa 1.50 m (1H400x400x13x21)
studies of Clough and O'Rourke (1990), Kung et al.
2.0 m
1 2.50 m (2007), Wang et al. (2010), and Ou and Hsieh (2011).
In the benchmark analysis, the soil layers of CL type
SM, N = 5-11 (clayey soils) were simulated with MC model in a total
6.5 m 6.30 m (2H400x400x13x21)
CL, Su = 21MPa 2 7.30 m
8.0 m stress undrained analysis, and the soil layers of SM type
9.65 m (2H400x400x13x21)
3 10.65 m
(sandy soils) were modeled with MC model in an
effective stress drained analysis. The strength
13.00 m (2H400x400x13x21) parameters of soils were directly taken from laboratory
4 14.00 m
tests. Poisson's ratio was assumed as 0.35 for sand
SM, N = 5-17 5 16.80 m
17.0 m layers and 0.495 ( 0.5) for clay layers. For stiffness
parameters, as reported by Hsiung (2009), effective
Young's modulus of sandy silts and silty sands could be
23.5 m
SM, N = 5-17 obtained by the following equation:
0.9 m diaphram wall
E '  2000 N (kPa) (1)
SM, N = 5-17
28.5 m in which, N is blow counts in Standard Penetration
CL, Su = 84MPa
30.5 m
32.0 m
Tests (SPT).
According to the studies of Bowles (1996), Lim et al.
(2010), Likitlersuang et al. (2013), and Khoiri and Ou
(2013), undrained Young's modulas of clayey silts and
42.0 m
SM, N = 18-26 silty clays could be estimated as follows:
Eu  500 Su (2)

in which, Su is undrained shear strength of clays.

The diaphragm wall was simulated by plate
SM, N = 28-42 elements, and the steel struts were simulated by
60.0 m
Mudstone elements of node-to-node anchor. The linear elastic
model was adopted to simulate both the diaphragm wall
Fig. 1. Cross section and ground condition of Case A. and steel struts. Poisson's ratio was taken to be 0.2 for
both the diaphragm wall and steel struts. The Young's
A three-dimensional finite element analysis will be modulus of the diaphragm wall was calculated by the
conducted to simulate the excavation of Case A, which equation of AASHTO (1998) as follows:
is called "benchmark analysis". The commercial For concrete with normal density:
software PLAXIS 3D, version 2013, is selected as a
numerical tool for 3D analyses in this study. PLAXIS Ec  4800 f 'c ( MPa ) (3)
3D is a three-dimensional finite element program,
which is developed at Deft University of Technology in in which, f'c (MPa) is the standard compressive strength
the Netherlands and is made commercially available by of the wall concrete. The Young’s modulus of steel
Plaxis Bv, Amsterdam, the Netherlands. struts was taken to be 2.1x105 MPa. The Young’s
As be seen from Fig. 1, the influences of clay layers moduli of both the diaphragm wall and steel struts were
on the excavation behavior are not significant due to reduced by 30% to consider the cracks in the diaphragm
their very thin thicknesses. In contrast, the sand layers wall due to bending moments and to consider the

1603
repeated uses and improper installation of steel struts The main reason can be possibly related to this fact that
(see Ou, 2006). the MC model dose not consider the strain-dependent
Fig. 2 below shows results obtained from the stiffness behavior or the small train characteristics that
benchmark analysis. The field measurements are also invole high stiffness modulus at small strain levels of soil.
shown in this figure for comparison, and only wall Only a single Young’s modulus of soil can be used for
deflection is demonstrated in this paper. MC model and such model also dose not distinguish
As can be noted from Fig. 2, the walls behave as a between loading and unloading/reloading stiffness. It is
cantilever at the first stage of excavation because the thus concluded that Young's modulus of sand adopted in
the MC model was thus underestimated at the earlier
struts at the first level have not yet installed and
stages of excavation due to wider range of small strain soil
preloaded. The walls then become prop-mode at the
area at these stages, and it thus leads to over-prediction of
subsequent stages of excavation. wall displacements herein.
The predicted wall deflections at the earlier stages of
Furthermore, it is aware that the heave of excavation
excavation (Stages 1, 2 and 3) in general are all slightly
bottom is significantly over-predicted because the MC
larger than the field measurements, respectively, no matter
model does not consider the higher stiffness of ground
in long side or short side. Once excavation goes deeper,
below the excavation level that is unloaded during
the predicted wall deflections at the later stages (Stages 4
excavation process. The over-prediction of the heave of
and 5) are slightly smaller than the field measurements at
excavation bottom then causes larger displacements at the
upper parts of the wall but comparatively larger than the
lower wall parts and smaller displacements at the upper
field measurements at lower parts of the wall, respectively.
wall parts.
However, the maximum lateral wall displcamenets from
both prediction and measurements are close to each other.

Displacement (mm) Displacement (mm)

-30 0 30 60 90 -30 0 30 60 90
0 0
The 1st strut: The 1st strut:
1.5m 1.5m

5 5

10 10
Depth (m)

The 5th
Depth (m)

15 stage: 16.8m 15 The 5th stage:

16.8m

Stage 1-Field Stage 1-Field

20 Stage 2-Field 20 Stage 2-Field

Stage 3-Field Stage 3-Field
Stage 4-Field Stage 4-Field
25 Stage 5-Field 25 Stage 5-Field

Stage 1-Analysis Stage 1-Analysis

Stage 2-Analysis Stage 2-Analysis

30 Stage 3-Analysis 30 Stage 3-Analysis

Stage 4-Analysis Stage 4-Analysis

Stage 5-Analysis Stage 5-Analysis

35 35
a) At the center of long side b) At the center of short side
Fig. 2. Comparisons of measured and computed wall deflections from the benchmark analysis of Case A.

3 THREE-DIMENSIONAL INFLUENCE values of excavation length and width. First, 3D

analyses were conducted to simulate Case A with
To evaluate the three-dimensional (3D) effect on the
various values of excavation length and width. Next,
wall deflection, a series of parametric studies was
analyses in plane strain condition (2D analyses) were
conducted on the excavation of Case A with various
also carried out with various values of excavation width

1604
to obtain 2D results for comparison with 3D results. is defined to further evaluate three-dimensional effect
The excavation length (L) was varied in values of 20 m, on wall deflection as follows:
40 m, 60 m, 80 m and 100 m. The excavation width (B)
A
was varied in values of 10 m, 20 m, 40 m, 60 m, 80 m PSR  (4)
and 100 m.  plane
In order to ignore influences from complex
characteristics of structure and soil, input parameters of in which, δ is the maximum wall deflection at a certain
structure and ground of said parametric studies were the section along the wall, and plane is the maximum wall
same as those used in the benchmark analysis. Thus, deflection of the section under plane strain condition.
only geometries of excavation (width and length) were Fig. 3 below presents results obtained from the
varied in the parametric studies. parametric studies.
As indicated previously, the Plane Strain Ratio (PSR)

Distance from the corner (m) Distance from the corner (m)
0 10 20 30 40 50 0 10 20 30 40 50
0 0
L=20 m L=20 m
Plane strain ratio (PSR)

Plane strain ratio (PSR)

0.2 L=40 m 0.2 L=40 m
L=60 m L=60 m
0.4 L=80 m 0.4 L=80 m
L=100 m L=100 m
0.6 0.6

0.8 0.8
B=10 m B=20 m
1 1
Distance from the corner (m) Distance from the corner (m)
0 10 20 30 40 50 0 10 20 30 40 50
0 0
L=20 m L=20 m
Plane strain ratio (PSR)

Plane strain ratio (PSR)

0.2 L=40 m 0.2 L=40 m

L=60 m L=60 m
0.4 L=80 m 0.4 L=80 m
L=100 m L=100 m
0.6 0.6

0.8 0.8
B=40 m B=60 m
1 1
Distance from the corner (m) Distance from the corner (m)
0 10 20 30 40 50 0 10 20 30 40 50
0 0
L=20 m L=20 m
Plane strain ratio (PSR)

Plane strain ratio (PSR)

0.2 L=40 m 0.2 L=40 m

L=60 m L=60 m
0.4 L=80 m 0.4 L=80 m
L=100 m L=100 m
0.6 0.6

0.8 0.8
B=80 m B=100 m
1 1
Fig. 3. Variation of the PSR with distance from the corner for constant excavation widths and various excavation lengths of Case A.

As shown in Fig. 3, PSR values corresponding to for the case having L of 20 m. For case having B of
various values of distance from the corner (d), length (L) 20 m, wall displacements are under plane strain
and width (B) of excavation are plotted. Once PSR condition for the place at 25 m from the corner, except
value is greater, it means displacement at such location cases of L of 20 m and 40 m. However, the maximum
is closer to plane strain condition. For the case of B of PSR value cannot reach 1.0 for B larger than 20 m. It
10 m, wall deformations all reach plane strain condition can only be 0.78 for case of B of 100 m. It is because
once it is approximately 20 m from the corner, except the arching effect on the wall displacements is inversely

1605
proportional to B/L. For case of B of 100 m, the wall 4 APPLICATION OF PSR
deflections cannot reach the plane strain condition
The relationship of Fig. 4 can to be used to estimate
because the minimum ratio of B/L is equal to 1.0.
3D maximum wall deflections by using 2D finite
Considering results shown from Fig. 3 above, a
element analyses for the excavation of Case A and
relationship between B/L and distance from the corner
other similar cases. Fig. 6 below shows comparison of
(d) for various PSR values is established by the linear
the maximum wall deflections at the section on the long
regression method as presented in Fig. 4. Results made
side of Case A that is 10 m away from the corner.
by Ou et al. (1996) based on excavations in clay are
It is clear that the 2D finite element analysis (2D
shown in Fig. 5 for comparison purpose.
FEA) is not able to accurately predict the maximum
The PSR values interpreted from the same d and
wall deflections at the section on the long side of Case
B/L in Fig. 4 and Fig. 5 are compared. It is clearly seen
A that is 10 m away from the corner, which are heavily
that for B/L more than 1.0, PSR is smaller than in sand
influenced from the three-dimensional effect. The PSR
rather than in clay with the same value of d. For
value (Ou et al., 1996), which is based on excavations
example, for B/L of 2.0 and d of 10 m, PSR is equal to
in clays, is also not able to well estimate the maximum
0.28 for sand and 0.33 for clay. In contract, for B/L less
wall deflections. The maximum wall deflections
than 0.5, PSR is larger than in sand rather than in clay
predicted from the PSR value (this study) are very close
with the same value of d. For example, for B/L of 0.25
to those obtained from 3D finite element analysis (3D
and d of 10 m, PSR is equal to 0.77 for sand and 0.53
FEA).
for clay. Further verification may have to be delivered
later in order to explore the reason for the difference.
120

5.0 2D FEA
PSR=0.1 100
0.1 0.2 0.3 0.4 0.5
Max. wall deflection (mm)

4.5
PSR=0.2
4.0 80 3D FEA
PSR=0.3
3.5
PSR=0.4 60
3.0 0.6
PSR=0.5 PSR value
B/L

2.5 40 (this study)

PSR=0.6
2.0
0.7 0.9 PSR=0.7 20 PSR value
1.5
(Ou et al., 1996)
1.0 0.8 PSR=0.8
0.9 0
0.5 PSR=0.9 0 5 10 15 20
1.0
PSR=1.0
Excavation depth (m)
0.0
0 10 20 30 40 50
Fig. 6. Comparison of the maximum wall deflections at the
Distance from the corner, d (m)
section on the long side of Case A that is 10 m away from the
Fig. 4. Relationship between B/L and distance from the corner corner.
(d) for various PSR values.
5 CONSLUSION
The following are conclusions drawn from this
study:
1) A restriction of the benchmark analysis in
prediction of the wall lateral displacements is seen. It is
because the benchmark analysis uses a basic soil model,
the MC model, which only adopts a single Young's
modulus of soil, dose not distinguish between loading
and unloading/reloading stiffness, and also does not
consider the strain-dependent stiffness behavior of soil.
2) The PSR increases gradually with increase of
distance from the corner (d). The PSR increases quickly
when the d is still small and then reaches to a certainly
constant value. For excavations of 10 m and 20 m
Fig. 5. Relationship between B/L and distance from the corner widths, the PSR reaches to 1.0 (i.e. wall displacements
(d) for various PSR values of excavations in clays (Ou et al., are under plane strain condition) as the d is larger than
1996). 20 m and 25 m, respectively. However, the maximum
PSR value cannot reach 1.0 for B larger than 20 m. It
can only be 0.78 for case of B of 100 m.

1606
 3) The PSR at a certain d value increases 10) Lin, D. G. and Woo, S. M. (2007): Three-dimensional
gradually with decrease of B/L. This result is similar to analyses of deep excavation in Taipei 101 construction
project, Journal of GeoEngineering, TGS, Vol. 2, No. 2, pp.
the previous study of Ou et al. (1996) for excavations in
29-41.
clays. However, for B/L more than 1.0, PSR is smaller 11) Ou, C. Y. (2006): Deep excavation: Theory and Practice,
than in sand rather than in clay. In contract, for B/L less Taylor & Francis, Netherlands.
than 0.5, PSR is larger than in sand rather than in clay. 12) Ou, C. Y. and Shiou, B. Y. (1998): Analysis of the corner
Further verification may have to be delivered later in effect on excavation behaviors, Canadian Geotechnical
order to explore the reason for the difference. Journal, Vol. 35, pp. 532-540.
13) Ou, C. Y., Chiou, D. C., and Wu, T. S. (1996):
4) A simple relationship for predicting PSR
Three-dimensional finite element analysis of deep
values or 3D maximum wall deflections by using 2D excavations, Journal of Geotechnical Engineering, ASCE,
finite element analyses is proposed as Fig. 4 for the Vol. 122, No. 5, pp. 337-345.
excavation of Case A. It is understood that such 14) Ou, C. Y., Shiau, B. Y., and Wang, I. W. (2000):
relationship does not consider all factors that affect the Three-dimensional deformation behavior of the Taipei
excavation behavior such as excavation sequence, national enterprise center (TNEC) excavation case history,
Canadian Geotechnical Journal, Vol. 37, pp. 438-448.
excavation duration, surcharge condition, existence of
15) Ou, C.Y. and Hsieh, P. G. (2011): A simplified method for
adjacent structures and construction method of the wall. predicting ground settlement profiles induced by excavation
However, the relationship developed in this study can in soft clay, Computers and Geotechnics, Vol. 38, pp.
still be applicable for other excavations with the same 987-997.
conditions. 16) PLAXIS (2013): Reference Manual, Plaxis BV, Amsterdam,
the Netherlands.
17) Schweiger, H. F. (2009): Influence of constitutive model and
ACKNOWLEDGEMENTS
EC7 design approach in FEM analysis of deep excavations,
The authors are grateful to Mr. Hsin-Nan Huang and Proceedings of ISSMGE international seminar on deep
Mr. Wei-Ya Song, graduate students of Department of excavations and retaining structures, Budapest, pp. 99-114.
18) Wang, J. H, Xu, Z.H, and Wang, W. D. (2010): Wall and
Civil Engineering, National Kaohsiung University of ground movements due to deep Excavations in Shanghai soft
Applied and Sciences, Kaohsiung, Taiwan, for helping soils, Journal of Geotechnical and Geoenvironmental
to collect the field data used in the study. Engineering, ASCE, Vol. 136, No. 7, pp. 985-994.

REFERENCES
1) AASHTO (1998): AASHTO LRFD Bridge design
specifications, American Association of Stage Highway and
Transportation Officials, Washington, US.
2) Bowles, J. E. (1996): Foundation analysis and design, 5th
edition, McGraw Hill Book Company, New York, USA.
3) Clough, G. W. and O’Rourke, T. D. (1990):
Construction-induced movements of in situ walls, Design
and Performance of Earth Retaining Structures, ASCE
Special Publication, No. 25, pp. 439-470.
4) Hsieh, P. G., Kung, T. C., Ou, C. Y., and Tang, Y. G. (2003):
Deep excavation analysis with consideration of small strain
modulus and its degradation behavior of clay, Proceedings of
12th Asian Regional Conference on Soil Mechanics and
Geotechnical Engineering, Singapore, Vol. 1, pp. 785-788.
5) Hsiung, B. B. C. (2009): A case study on the behavior of a
deep excavation in sand, Computers and Geotechnics, Vol.
36, pp. 665-675.
6) Khoiri, M., and Ou, C. Y. (2013): Evaluation of deformation
parameter for deep excavation in sand through case histories,
Computers and Geotechnics, Vol. 47, pp. 57-67.
7) Kung, G. T. C, Juang, C. H., Hsiao, E. C. L., and Hashash, Y.
M. A. (2007): Simplified model for wall deflection and
ground-surface settlement caused by braced excavation in
clays, Journal of Geotechnical and Geoenvironmental
Engineeing, ASCE, Vol. 133, No. 6, pp. 731-747.
8) Likitlersuang, S., Surarak, C., Wanatowski, D., Oh, E., and
Balasubramaniam, A. (2013): Finite element analysis of a
deep excavation: A case study from the Bangkok MRT, Soils
and Foundations, Vol. 53, No. 4, pp. 498-509.
9) Lim, A., Ou, C. Y., and Hsieh, P. G. (2010): Evaluation of
clay constitutive models for analysis of deep excavation
under undrained conditions, Journal of GeoEngineering,
TGS, Vol. 5, No. 1, pp. 9-20.

1607