Analysis of Field Performance of Embankments On Soft Clay Deposit With and Without PVD-improvement
Analysis of Field Performance of Embankments On Soft Clay Deposit With and Without PVD-improvement
Analysis of Field Performance of Embankments On Soft Clay Deposit With and Without PVD-improvement
Received 19 September 2004; received in revised form 21 May 2005; accepted 31 May 2005
Abstract
This paper presents a case history of the performance of two full-scale test embankments
constructed on soft clay deposit in the eastern coastal region of China. One embankment was
constructed on natural subsoil and the other was constructed on prefabricated vertical drain
(PVD) improved subsoil. The thickness of the soft clay deposit without PVD-improvement
was 19 m and with PVD-improved case was 23 m. The PVDs were installed to a depth of 19 m
with spacing of 1.5 m in a triangular pattern. Field performance of the two embankments was
analyzed using the nite element method. The following inuential factors: (i) hydraulic
conductivity of subsoil in eld and (ii) discharge capacity of PVDs were investigated
numerically. The back-analyzed results for the embankment on natural subsoil showed that
the hydraulic conductivity ratio (Cf) of eld to laboratory values is about 6. To analyse the
PVD-improved subsoil, a simple approach using the equivalent vertical hydraulic conductivity
of PVD-improved subsoil was employed. The analytical results show that PVDs increased the
bulk vertical hydraulic conductivity of soft subsoil by about 30 times compared to the original
0266-1144/$ - see front matter r 2005 Elsevier Ltd. All rights reserved.
doi:10.1016/j.geotexmem.2005.05.002
ARTICLE IN PRESS
464 S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485
Nomenclature
non-treated subsoil. The discharge capacity of PVDs in this eld case is 79100 m3/a, which is
consistent with the ndings from laboratory tests and other reported values in literature.
r 2005 Elsevier Ltd. All rights reserved.
1. Introduction
SPIA
Taihu
Lake Shanghai
N
Expressway
Bay
gz hou
Han
Shaoxing
NIA
0 15 30 km
Ningbo
Table 1
Test embankments with ground improvement descriptions
Note: H embankment height; Hps Height of preloading embankment; n side slope of embankment;
DT thickness of crust layer; Dclay thickness of clay layer; PVD plastic drain improved; t thickness
of PVD; P drain pitch; D depth of drain.
2. Soil conditions
The soil prole and soil properties of the soft clay deposits at the test site are
shown in Fig. 2. The soil prole consists of a thin weathered crust (TC) from about
11.5 m thick overlying a silty clay (SC1) deposit of approximately 4 m thick. The
third layer is very soft clay (MC, in China it is called mucky clay) with a thickness of
approximately 10 m for the PVD-improved case and 11.5 m for the unimproved case.
Below the clay is a silt clay layer called mucky-silty clay (MSC) approximately 4 m
thick for PVD-improved case and 2.7 m for the unimproved case followed by a
medium to stiff silty clay layer (SC2) from 3 to 5 m thickness extending down to 23 m
depth, which is in turn underlain by a layer of loose clayey sand (CS). The soil layer
separation is plotted in Fig. 2, in which solid line denotes the PVD-improved subsoil
and dashed line denotes the subsoil without PVD-improvement. The thickness of the
soft clay deposit without PVD-improvement is about 19 m and that with PVD-
improvement is about 23 m. The soft silt clay and very soft clay have water contents
greater than their liquid limits, low hydraulic conductivity and lower shear strength.
The parameters such as compression index (Cc) and hydraulic conductivity (k) were
determined from laboratory oedometer tests; the values of kh and kv were obtained
from oedometer test results using the samples oriented in horizontal and vertical
directions, respectively. Shear strength Su were determined from eld vane shear test.
Fig. 3 illustrates the cross-section and plan view of the test embankments,
arrangement of PVDs, and instrumentations. As seen in Fig. 3, a 0.5 m thick sand
ARTICLE IN PRESS
S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485 467
Soil 3
LL , PL (kN/m ) Cc Su(kPa)
-10
kh (10 m/s)
-10
kv (10 m/s)
Description wn (%)
Top weathered 20 30 40 50 18 20 0.2 0.4 0.6 0 40 80 120 10 30 50 10 20 30 40
0
crust
Silt clay 1
5
Very soft Selected Selected Selected
clay
10 Selected
Depth (m)
15 Soft mucky
Silt clay
Fig. 2. Soil prole and properties at the test site of thinner deposit section.
Note: LL-liquid limit, PL-plastic limit, wn-water content, g-unit weight, Cc-compressive index, Su-Field
vane shear strength, kh, kv-hydraulic conductivity in horizontal and vertical direction from laboratory test.
10 26 m
TC: top weathered crust
:5 SC: silty clay
1:1
H
SC1
CS: clayey sand
MC
-10 piezometers
surface settlement gauges
MSC
inclinometers
MSC
SC2 piezometers
-20 SC2
inclinometers
CS CS
1.5 m
100 mm
0.867 m
6 mm
1.3 m
PVDs
26 m
piezometer 0.75 m
Fig. 3. Cross-section and plan view of embankment, eld instrumentation, and arrangement of PVDs in
eld.
ARTICLE IN PRESS
468 S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485
6 120
without PVD-improved case
5 100
Height (m) 4 80
Load (kPa)
3 Field data 60
simulated
2 40
1 20
0 0
0 200 400 600 800 1000 1200
(a) Elapsed time (days)
6 120
4 80
Height (m)
Load (kPa)
3 Field data 60
simulated
2 40
1 20
0 0
0 200 400 600 800 1000 1200
(b) Elapsed time (days)
Fig. 4. Load-time history of the embankment on: (a) Natural subsoil; (b) PVD-improved subsoil.
mat (hydraulic conductivity: k40.001 m/s) was placed on the soft ground initially.
Then, decomposed granite was lled and compacted to a unit weight of 20 kN/m3.
The loading-time history relations for the two embankments are plotted in Fig. 4.
The height of the embankment on PVD-improved subsoil is 5.88 m and that on
natural subsoil is 4.66 m. For the improved case, PVDs were installed in a triangular
pattern with the spacing of 1.519 m deep. Table 2 lists the properties of PVDs used
in the test embankment. The discharge capacity, Qw, provided by manufacturer is
greater than 1500 m3/year.
Field behaviour of soft ground under embankment loading were observed through
monitoring the surface settlement, lateral displacements of the subsoil at the middle
point of the embankment side slope, and pore pressure variation by piezometers at
different depths. The layout of the eld instrumentation is shown in Fig. 3.
ARTICLE IN PRESS
S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485 469
Table 2
Size and physical properties of PVDs used in this eld project
Note: Qw discharge capacity of PVD, the value in table was provided by manufacturer.
Table 3
Cf value for some clay deposits
Bangkok clay (100 km from sea) 25 Back calculated Chai et al. (1995)
Bangkok clay (close to sea) 4 Back calculated Chai et al. (1996)
Malaysia Muar clay deposit 2 Back calculated Chai and Bergado (1993)
Ariake clay 4 Back calculated Chai and Miura (1999)
Louiseville (Canada) 1 Self-boring permeameter Tavenas et al. (1986)
St-Alban (Canada) 3 Self-boring permeameter Tavenas et al. (1986)
Very soft mucky clay (China) 6 Back calculated Present study
Note: Cf Ratio between eld horizontal hydraulic conductivity kh value to laboratory value.
The hydraulic conductivity of subsoil in the eld is different from the laboratory
measured value, esp. calculated from laboratory oedometer test result. Laboratory
test normally underestimates the eld values (Tavenas et al., 1986). The relation can
be expressed as following formula:
kh f C f kh l or kv f C f kv l , (1)
For this eld case Cf 6 obtained, and the details will be discussed in the following
context.
The parameters for PVDs are listed in Table 4 and the determination methods are
presented as follows.
Table 4
Parameters related to the behaviour of PVD
where ks hydraulic conductivity in smear zone; and the other symbols are the same
as presented before.
The hydraulic conductivity in smear zone varies with the radius from the centee of
PVD. From laboratory test on various clay samples, the average values of kh =ks l
varies from 1 to 5 (Hansbo, 1987). For Ariake clay it is about 2.5 (Madhav et al.,
1993); for preconlidated kaolin clay, the ratio is about 3 (Hird and Moseley, 2000).
However, for this particular soil in this study, there is no available laboratory test
value. Zhu et al. (1992) reported a case history of sand drain improved subsoil
under the runway of Ninbo International Airport (see Fig. 1), where is about 100 km
away from the test site in this study. They estimated the value of kh =ks f in eld
of about 1215. Thus, due to C f 6, according to Eq. (6), kh =ks l 222:5 can
be evaluated.
ARTICLE IN PRESS
472 S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485
Table 5
Model parameter for subsoil in the test site
Layer n k l M e0 G gt (kN/m3) kh (103 m/d) kv (103 m/d) kve (103 m/d) kve/kv
TC 0.30 0.008 0.08 1.0 0.81 1.98 19.3 2.72 2.72 34.3 12.6
SC1 0.35 0.016 0.16 1.0 1.07 2.22 18.5 0.56 0.22 7.3 33.2
MC 0.35 0.028 0.28 0.8 1.36 3.22 17.3 2.54 1.69 28.9 17.1
SMC 0.35 0.018 0.18 0.8 1.10 2.98 17.9 2.06 1.01 24.1 23.9
SC2 0.30 0.010 0.10 1.0 0.81 1.99 19.3 0.39 0.18 5.2 28.9
CS 0.25 19.5 25.90 25.90
Fill 0.20 20.0 193.5 193.5
Note: TC Thin weathered crust; SC Silty clay; MC Very soft clay; SMC soft silty clay;
CS clayey sand; gt unit weight; e0 initial void ratio; l virgin loading slope in e-ln(p) plot (p is
effective mean stress); k reloading/unloading slope in e-ln(p) plot; and M slope of failure line in p
versus q plot (q is deviator stress); n Poissons ratio.
1993). The large deformation phenomenon can be considered by updating the nodal
coordinates during the incremental analysis. Usually, this operation does not include
the elements above the current construction level. As a result, the nal applied
embankment ll thickness will be larger than pre-dened value. There are two
methods to treat this problem: (i) include all the embankment elements into all steps
of incremental analysis and apply the gravity force by percentage in each step; or
(ii) apply the embankment elements layer by layer and correct the coordinate of the
elements above the current construction level according to the settlement of the
current construction top surface (Chai and Bergado, 1993). This study adopted Chai
and Bergados (1993) proposal.
In FEM analysis, the plain strain condition was assumed. The modelled range in
vertical direction was 25 m deep for the case without PVD-improvement and 29 m
deep for the case with PVD-improvement; and horizontally 120 m away from the
embankment centerline. The displacement boundary conditions were as follows: at
bottom, both vertical and horizontal displacements were xed, and for left and right
vertical boundaries, the horizontal displacement was xed. The adopted drainage
boundary conditions were as follows: the ground surface and bottom line (weathered
rock) were drained. The left and right boundaries were undrained. Fig. 5 shows the
FEM mesh for the embankment on natural subsoil and PVD-improved subsoil. The
element type for subsoil and the embankment ll is quadrilateral. Each element has
9 Gauss integral points.
The mechanical behaviour of the soft clay layers was represented by the Modied
CamClay model (Roscoe and Burland, 1968) and the clayey sand layer, the sand
ARTICLE IN PRESS
474 S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485
(5.88)
4.66
Top crust
PVD-improved zone
29.66(34.88)
25(29)
Mucky silt clay
120
for the top crust. The lateral earth pressure coefcient was calculated using the
equation proposed by Mayne and Kulhawy (1982).
0
K 0 1 sin f0 OCRsin f , (8)
where f0 is the effective internal friction angle of soil, OCR is over-consolidation
ratio.
The parameters for clayey sand were assumed as Youngs modulus E
25; 000 kPa and Poissons ratio v 0:25. The mechanical properties of the ll
material were assumed as follows: E 30; 000 kPa and v 0:25. The unit weight of
the ll material was 20 kN/m3. The ground-water level was about 1.5 m below
ground surface.
The aforementioned parameters such as Cf and qw were determined from the back-
calculation method. The procedure of back-calculation is as follows: (1) Back-
calculation of Cf based on the measured settlement of the embankment on natural
subsoil; (2) select kh =ks l 222:5, using Eq. (5), the eld value of kh =ks f
C f kh =ks l 12215 is obtained; (3) using above values to evaluate the value of qw
in the eld from PVD-improved case.
7.1.1. Cf value
Cf was obtained by analyzing the embankment on natural subsoil. In the analysis,
the initial values of hydraulic conductivity in Fig. 2 were multiplied by a factor Cf,
larger than 1.0 and the calculated settlement was compared with the measured
data. The settlement-time relation for the embankment on natural subsoil is plotted
in Fig. 6. As shown in the gure, with C f 6 (solid line), numerical result agrees well
with the measured values. For comparison, the calculated settlement-time curves
with respect to C f 1 are also plotted using dashed line in the same gure.
0.0
Field measured
center
0.2
right shoulder
left shoulder
Settlement (m) 0.4
0.6
Calculated
0.8
with Cf = 1
with Cf = 6
1.0
without PVD-improvedcase
1.2
0 200 400 600 800 1000 1200
Elapsed time (days)
Fig. 6. Measured and calculated values of settlement for the embankments on the subsoil without PVD-
improvement.
0.0
Field measured
3
Calculated (m /yr)
0.5 qw= 0, without PVD
qw= 10
qw= 79
Settlement (m)
1.5
2.0
settlement, and the results were plotted in Fig. 8. It can be considered that since the
hydraulic conductivity of soft clay is relatively low, when qw is larger than a certain
value, the hydraulic conductivity of soft clay controls the consolidation process
rather than qw. Only with this information, it can be concluded that the eld qw value
was larger than 79 m3/year. However, referring to other reported eld values (Holtz
et al. 1991; Li and Rowe, 2001), it is considered that qw values in the led was about
79100 m3/year. The nal adopted values are kh =ks f 13:5 and qw 100 m3 =a.
The nal calculated settlement-time relation for the embankment on PVD-improved
subsoil is plotted in Fig. 9.
ARTICLE IN PRESS
S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485 477
100
PVD section
-10.0m, in mucky clay (MC)
80
Excess pore pressure u (kPa) Field measured
3
FEM calculated, qw(m /yr)
60 qw = 0, without PVD
qw = 10
qw = 100
qw = 1500
40
20
0
0 200 400 600 800 1000 1200
Elapsed time (day)
Fig. 8. Variation of calculated excess pore pressure with discharge capacity qw.
0.0
with PVD-improved case
Field measured
0.5 right shoulder
center
Settlement (m)
left shoulder
1.0 calculated
1.5
2.0
Fig. 9. Measured and calculated values of settlement for the embankments without PVD-improvement.
Fig. 10 shows the observed excess pore pressures during embankment construc-
tion, which are presented as a function of the applied stress of the embankment at
the ground surface. As shown in Fig. 10a, for the embankment on natural subsoil,
the generation of pore pressure supports the tendency found by Tavenas and
Leroueil (1980), that is, two phases of pore pressure response during construction
can be identied: the rst, during which the increases in pore pressure are low; the
ARTICLE IN PRESS
478 S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485
H
=
u
40
20
0
0 20 40 60 80 100 120
(a) Load, H (kPa)
PVD-improved case
80 2.0 m
6.0 m
Excess pore pressure, u (kPa)
12 m
14 m
60 16 m
H
=
u
40
20
0
0 20 40 60 80 100 120
(b) Load, H (kPa)
Fig. 10. Pore pressures measured under the centre of the embankment test section, as a function of the
embankment load: (a) without PVD-improved case; (b) with PVD-improved case.
second, during which the increase in pore pressures is approximately equal to the
increase in total stress. For the embankment constructed on PVD-improved subsoil,
the measured pore pressure at a depth of 2 and 6 m still follows the ndings of
Tavenas and Leroueil (1980). This is because subsoil is partially drained and clay
begins to consolidation at the very low pressure during construction (Leroueil et al.,
1990). After the embankment pressure is approximately equal to the consolidation
yield stress, sp, of the soil. The soil becomes normally consolidated, and now posses
a coefcient of consolidation much lower than in the overconsolidated region, and it
ARTICLE IN PRESS
S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485 479
behaves essentially as undrained soil (Leroueil et al., 1990). However, the measured
pore pressures at deeper place did not agree with this tendency (see Fig. 10b). The
detailed discussion on the change of coefcient of consolidation, Cv, with
consolidation pressure will be discussed in the later context (see Fig. 12).
This phenomenon may be related to the following factors: (i) disturbance of
subsoil during installing PVDs; (ii) reduction of consolidation yield stress of crust
layer; (iii) piezometers might be installed within the smear zone. In this site, there is a
weathered top crust layer, which has much higher strength (about 5 times the soft
layer) and stiffness, as shown in Fig. 2. The stiff top crust may tend to distribute the
loading from the embankment over a larger area and reduce the level of stress to the
underlying softer subsoils. Moreover, top crust layer has a very high over-
consolidated ratio with OCR 5 and has a relative high coefcient of consolidation
(Cv). However, installing PVDs will disturb the top crust layer, reduce its
consolidation yield stress, and thus reduce Cv value. Hence, the increase of excess
pore pressure for the improved case becomes larger than the unimproved case at the
early stage of embankment construction in the shallower depth, i.e. less than 6.0 m
(see Fig. 10).
The reason for higher excess pore pressure for PVD-improved case in the deeper
places H412 m is not clear yet. One of the reasons is that piezometers might be
installed in the smear zone. In the eld, it is difcult to keep the verticality during
PVD installation and boring hole for piezometers installation. Piezometers were
designed to be installed at the mid point of PVDs (see Fig. 3). The distance between
piezometer and PVD is less than 867 mm. If there is only 21 inclining from vertical
direction during installation, the piezometers at 12 m depth will be installed within
the smear zone. If piezometers were installed in the smear zone, higher excess pore
pressure would be measured.
The comparison of the calculated and measured values of excess pore pressures
between PVD-improved case and unimproved case is plotted in Fig. 11. During
embankment construction, for the unimproved case, the simulated excess pore
pressure agrees the measured value fairly well. However, for the PVD-improved case,
the calculated values are much lower than the measured ones. For the both cases,
dissipation of calculated excess pore pressures is much faster than that of the
measured values.
7.3. Reasons for discrepancies between predicted and measured excess pore pressure
The discrepancy of excess pore pressure between FEM calculated results and eld
measured data might be due to possible errors in measurement and limitation of
calculation methods.
(1) Possible errors of the measured data: Pore pressure is measured by piezomters.
Pore water seeps into water chamber through a lter of porous stone. Fine
soil particles might inltrate into the pore of lter to cause clogging
ARTICLE IN PRESS
480 S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485
80 80
without PVD section with PVD section
60 60 -2.0m, in silty clay (SC1)
-2.5m, in silty clay (SC1)
Calculated 40 Calculated
40 Measured Measured
20 20
Excess pore pressure u (kPa)
40 40
20 20
0 0
80 80
-16.05m, in mucky silty clay (MSC) PVD section
60 60 -16.05m, in mucky silty clay (MSC)
40 40
20 20
0 0
0 200 400 600 800 1000 1200 0 200 400 600 800 1000 1200
(a) Elapsed time (day) (b) Elapsed time (day)
Fig. 11. Measured and calculated excess pore pressure-time curves: (a) without PVD; (b) with PVD.
of the piezomter. Hence, a certain amount of excess pore pressures are locked
in the piezometer. Some measurements on the excess pore pressure showed
that similar behaviour of excess pore pressures was observed, that is, with
progressing of consolidation, there were little excess pore pressure dissipations
as reported by Chai and Miura (1999), Hird et al. (1995), and Indraratna
et al. (1994). For the test embankments on Muar clay, Malaysia, some
piezometers showed no excess pore pressure dissipation during the consoli-
dation period. After installing additional piezomters aside the existing points,
close to the predicted values of excess pore pressures were measured (MHA,
1989).
(2) Another reason might be that piezometers were installed in smear zone as
discussed previously. The modelling method can not simulate this aspect.
(1) Limitation of modelling method on PVD behaviour: In the eld, a single PVD
works in a close axisymmetric drainage condition (three-dimensional (3D)). In
the present analysis plain strain condition is assumed (two-dimensional (2D)
condition) and PVD behaviour is modelled by considering the increased vertical
hydraulic conductivity in a macro sense. The equivalent conductivity is obtained
ARTICLE IN PRESS
S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485 481
by the average degree of consolidation. The smear effect is taken into account
through parameter m (Eq. (3)) by considering the average hydraulic conductivity,
ks, in smear zone. However, the variation of hydraulic conductivity in smear zone
was not taken into account. As a result, local pore pressure distribution in
horizontal direction is different between 3D and 2D conditions.
(2) Reduction of qw value with elapsed time: For discharge capacity, qw 100 m3/
year was used. As discussed in the previous context, increasing qw will have no
effect on the dissipation of excess pore pressure, however, reducing qw will
certainly affect the consolidation rate as illustrated in Fig. 7.
(3) Variation of Cv with consolidation pressure: In this study, the behaviour of
soft clay was modelled by Modied Camclay model (Roscoe and Burland,
1968). Rising and dissipation of excess pore pressure is related to the coefcient
of consolidation Cv. For linear e-lnp relationship, Cv can be expressed as
follows:
1 ep0 k
Cv , (9)
lgw
where k the hydraulic conductivity of soil, gw the unit weight of water, and
l the slope of virgin compression curve in e-lnp plot and p is effective mean
stress.
As shown in Eq. (9), Cv, is related to the void ratio, compressive index, stress level,
and hydraulic conductivity. Moreover, hydraulic conductivity k varies with the void
ratio (Tavenas et al., 1986; Taylor, 1948). Fig. 12 shows a comparison of the
calculated values Cv using Eq. (9) and measured Cv for Ariake clay (Chai et al.,
2004b). As shown in the gure, even considering the variation of k with e, there are
still discrepancies between the calculated values and the measured data.
-5
10
Cv = coefficient of consolidation
p = yield stress
c
-6
10
Cv (m2/s)
Measured
-7
10 Calculated
k constant
k varies with e (eq.7)
-8 pc
10
Vertical effetive stress, logP (kPa)
Fig. 12. Variation of the coefcient of consolidation with vertical effective stress (Chai et al., 2004b).
ARTICLE IN PRESS
482 S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485
7.4. Settlement
The calculated and measured values for the lateral displacement are shown in
Fig. 13. As illustrated in Fig. 13, the FEM can simulate the prole of lateral
displacement very well for the unimproved case. However, there are some
discrepancies between measured and simulated ones. Although the lateral displace-
ment at the end of construction for the PVD-improved case is much larger than that
for unimproved case, for the same embankment thickness 4.66 m, the calculated
lateral displacements are not different much. For the PVD-improved case, at
embankment thickness of 5.18 m there was a leaving period of 82 days (see Fig. 4b)
and this period might increase the lateral displacement much. If checking Fig. 13
carefully, for PVD-improved subsoil, the measured results show that the place where
maximum lateral displacement occurred is about 4 m deeper than that for the
unimproved case, however, the calculated values are the same as the unimproved one
and about 6 m from the ground surface. The exact reason is not clear but is probably
due to the spatial variation of subsoil condition.
8. Conclusions
The following conclusions are based on the eld measured data and the associated
FEM results:
0
without PVD case
10
Depth (m)
End of construction
Measured
20 Calculated
160 days later
Measured
Calculated
30
0 100 200 300 400 500
(a) Lateral displacement (mm)
0
with PVD case
10
Height = 4.66m
Depth (m)
Calculated
End of construction
Calculated
20 Measured
After 8.5 months
Calculated
Measured
30
0 200 400 600 800 1000
(b) Lateral displacement (mm)
Fig. 13. Measured and calculated lateral displacement: (a) without PVD; (b) with PVD.
3. In the macro-sense, installation of PVD with the spacing of 1.5 m makes increased
the vertical bulk hydraulic conductivity of the soft subsoils by about 30 times
compared to that for the untreated soil.
4. During embankment construction, the excess pore pressure measured for PVD-
improved subsoil were quite different from those without PVD-improvement. The
rapid increase in the excess pore pressure at shallower depths for PVD-improved
subsoil may be due to disturbance arising from PVD installation.
ARTICLE IN PRESS
484 S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485
References
Asaoka, A., Nakano, M., Matsuo, M., 1992. Prediction of partially drained behavior of soft clays under
embankment loading. Soils and Foundations 32, 4158.
Balaam, N.P., Booker, J.R., 1985. Analysis of rigid rafts supported by granular piles. International
Journal for Numerical and Analytical Methods in Geomechanics 5, 379403.
Bergado, D.T., Anderson, L.R., Miura, N., Balasubramaniam, A.S., 1996. Soft Ground Improvement, in
Lowland and other Environments. ASCE Press, pp. 427.
Chai, J.C., Bergado, D.T., 1993. Some techniques for FE analysis of embankment on soft ground.
Canadian Geotechnical Journal 30 (4), 710719.
Chai, J.C., Miura, N., 1999. Investigation of factors affecting vertical drain behavior. Journal of
Geotechnical and Geoenvironmental Engineering 125 (3), 216226.
Chai, J.C., Miura, N., Sakajo, S., Bergado, D.T., 1995. Behaviour of vertical drain improved subsoil under
embankment loading. Soils and Foundations 35 (4), 4961.
Chai, J.C., Bergado, D.T., Miura, N., Sakajo, S., 1996. Back calculated eld effect of vertical drain.
Proceedings of the second International Conference Soft Soil Engineering. Hohai University, Nanjing,
vol. 1, pp. 270275.
Chai, J.C., Shen, S.L., Miura, N., Bergado, D.T., 2001. A simple method of modeling PVD improved
subsoil. Journal of Geotechnical and Geoenvironmental Engineering 127 (11), 965972.
Chai, J.C., Miura, N., Nomura, T., 2004a. Effect of hydraulic radius on long-term drainage capacity of
geosynthetics drains. Geotextiles and Geomembranes 22, 316.
Chai, J.C., Miura, N., Zhu, H.H., Yudhbir, 2004b. Compression and consolidation characteristics of
structured natural clay. Canadian Geotechnical Journal 41 (12), 12501258.
Hansbo, S., 1981. Consolidation of ne-grained soils by prefabricated drains. Proceedings of the
Tenth International Conference on Soil Mechanics and Foundation Engineering, Stockholm, 3,
pp. 677682.
Hansbo, S., 1987. Design aspects of vertical drains and lime column installations. Proceedings of Ninth
Southeast Asian Geotechnical Conference. Southeast Asian Geotechnical Society, Bangkok, vol. 2(8),
pp. 112.
Hird, C.C., Moseley, V.J., 2000. Model study of seepage in smear zones around vertical drains in layered
soil. Geotechnique 50 (1), 8997.
Hird, C.C., Pyrah, I.C., Russell, D., 1992. Finite element modelling of vertical drains beneath
embankments on soft ground. Geotechnique 42 (3), 499511.
Hird, C.C., Pyrah, I.C., Russell, D., Cinicioglu, F., 1995. Modelling the effect of vertical drains in two-
dimensional nite element analyses of embankments on soft ground. Canadian Geotechnical Journal
32, 795807.
Holtz, R.D., Lancellotta, R., Jamiolkowski, M., Pedroni, S., 1991. Laboratory testing of prefabricated
wick drains. Proceedings of Geo-Coast91, Yokohama, Japan, pp. 311316.
Holtz, R.D., Shang, J.Q., Bergado, D.T., 2001. Soil improvement. Geotechnical and Geoenvironmental
Engineering Handbook. Kluwer Academic Publishing, Norwell, USA, pp. 429462 (Chapter 15).
Indraratna, B., Balasubramaniam, A.S., Ratnayake, P., 1994. Performance of embankment stabilized with
vertical drains on soft clay. Journal Geotechnical Engineering 120 (2), 257273.
Jamiolkowski, M., Lancellotta, R., Wolski, W., 1983. Pre-compression and speeding up consolidation,
general report. Special Session 6, Proceedings of Eight European Conference on Soil Mechanics and
Foundation Engineering. Balkema, Rotterdam, pp. 12011226.
Leroueil, S., Magnan, J.P., Tavenas, F., 1990. Embankments on soft clays. English Edition, translated by
D.M. Wood. Ellis Horwood, New York, 360pp.
Li, A.L., Rowe, R.K., 2001. Combined effects of reinforcement and prefabricated vertical drains on
embankment performance. Canadian Geotechnical Journal 38 (1), 12661282.
Madhav, R., Park, Y.M., Miura, N., 1993. Modelling and study of smear zones around band shaped
drains. Soils and Foundations 33 (4), 135147.
Mayne, P.W., Kulhawy, F.H., 1982. K0-OCR relationships in soils. Journal of Geotechnical Engineering
108 (GT6), 851872.
ARTICLE IN PRESS
S.-L. Shen et al. / Geotextiles and Geomembranes 23 (2005) 463485 485