Piled Raft Behaviour
Piled Raft Behaviour
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
Behaviour of Piled Rafts with Piles of Different Lengths and Diameters under
Vertical Loading
Abstract
The use of piled raft foundations has become popular in recent years as the system
can act as a settlement reducer. When the raft is subjected to non-uniform loadings,
piles of different diameters and lengths may be required to prevent the foundation
from rotatingand to reduce the differential settlement . This paper presents a finite
layer method used for the analysis of piled rafts with piles of different lengths and
diameters. The soil is divided into horizontal layers with different material properties
and only vertical loads may be applied to the raft. Interactions between raft-soil-pile
are computed. Results from the present method are shown to agree closely with those
from the finite element method.
Introduction
The use of a piled raft as the foundation for buildings has proven to be an effective
and economic way to control total and differential settlements as well as improving
bearing capacities. In the conventional foundation design approach, piles are
designed to carry the majority of the load, therefore, the behaviour of the foundation
system is governed by the pile group. In the design of a piled raft, the load is shared
between the raft and piles and it is necessary to take the complex soil-structure
interaction effects into consideration.
The behaviour of piled raft foundations has been studied by many researchers and
different techniques have been developed for the analyses in past decades. The
boundary element method which treated the raft as a plate, the piles as a series of
springs and the soil as an elastic continuum was used by Brown and Weisner (1975),
Chow (1986), Kuwabara (1989), Poulos (1994) and Sinha (1997). A combination of
finite element and boundary element methods for the analyses were used by Franke
et al. (1994), Russo and Viggiani (1998) and Mendonca and Paiva (2003). In this
approach, the raft is treated as a bending plate modeled by the finite element method
and the soil and pile are modeled by the boundary element method. Three-
dimensional finite element techniques have been used by Zhang et al. (1991),
Katzenbach and Reul (1997), and Reul and Randolph (2003 and 2004).
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
Finite Layer methods (Small and Booker 1986) were developed to compute the
deflections of horizontally layered soil. This technique can be used for the analysis of
piled rafts as it is capable of taking into account the raft-soil-pile interactions and has
been demonstrated for vertical loadings (Ta and Small 1996), horizontal loadings
(Zhang and Small 2000) and for rafts supported by different lengths of piles (Small
et al. 2004).
Method of Analysis
The analysis of a piled raft is carried out by a combination of finite layer and finite
element methods. The foundation system can be separated into two free body
diagrams (an isolated raft and a pile group) as shown in Figure 1. The raft is
subjected to external forces and contact stresses between the piles or the soil and the
raft. The pile group is subjected to the interface forces transferred from the raft to the
pile. The raft is divided into a number of rectangular elements with each pile head
fitted within one of the elements and the pile is divided into a number of rod
elements.
To carry out the analysis, four different interactions have to be considered. Figure 3
shows the following interactions:
(i) Pile-to-Pile – deflection of an unloaded pile due to an adjacent loaded
pile
{Pr}
Interface forces between the Interface forces between the
raft and piles raft and the soil
Interface forces transferred from raft to pile heads and soil surface
{Psp}
x
z
Soil Ring loads
acting on soil
nodes
Pile
Circular loads acting
on pile base
Figure 1. A raft and a pile group subjected to external forces and interface
forces in all directions (the y direction is not shown)
Material 1
Material 2
The pile-to-pile interaction can be obtained by using the stiffness approach (Zhang
and Small, 2000). A pair of piles being considered in the calculation can have
different lengths and diameters. Firstly, the influence matrix for the soil has to be
computed by applying a unit ring load at each soil layer along the pile shaft and a
unit uniform load at the pile base. Deflections computed at the locations of these ring
and base loads form the columns of the soil influence matrix [Isoil]. Deflections of the
soil at the locations of the ring or uniform loads can be expressed as
where fs are the forces acting at the soil interfaces, and s are the deflections of the
soil at the locations of ring or uniform loads. The soil stiffness matrix [Ksoil] can be
obtained by inverting the soil influence matrix.
Once the soil stiffness matrix has been constructed, we can assemble the pile
stiffness matrix [Kpile] for the two piles under consideration. The pile stiffness matrix
is simply the sum of individual element stiffness matrices.
where p are the deflections at the nodes of the pile, fp are the forces applied at the
nodes and P is the load vector of the applied load at the pile heads. By assuming that
no slip occurs along the pile shaft, the deflection of the soil and deflection at the pile
shaft are the same and the forces acting on the soil and along the pile shaft have
opposite sign. The pile-soil stiffness can be formed by adding the pile stiffness
matrix to the soil stiffness matrix.
The deflection of an unloaded pile due to the adjacent loaded pile can then be
computed by applying a unit load at the pile head.
The pile-to-soil interaction can be obtained by using the same stiffness approach as
for the pile-to-pile interaction. However, the stiffness matrix is set up for a single pile
only. By solving the equations, the magnitudes of the ring loads along the pile shaft
and circular load at the pile base can be computed. Then these loads are applied back
to the soil to solve for the surface deflections at the desired positions.
The interaction between the uniformly loaded soil surface and an unloaded pile (soil-
to-pile) can be computed by the use of the interaction method. Considering the soil
first, the deflection of the soil consists of two components – (i) due to the ring loads
along the pile shaft and (ii) the surface load applied to the soil. The deflection due to
the ring load can be computed by Equation 1, whilst the deflection due to the surface
load can be computed simply by the finite layer method.
where l are the deflections of the soil due to the surface load. The pile influence
matrix can be obtained by applying unit loads at each node to a pinned pile (pinned
at the pile head). The deflection of the pile can be written as
{ p} = [Ipile]{fp} + a (6)
where is the vertical translation for the pinned raft due to the surface load applied
on the soil and a = (1,1,1,…,1)T. The deflection of the pile can be computed by
combining these two equations with the assumption that no slip occurs along the pile
shaft (i.e. s = p and fs = - fp) and writing
I soil + I pile a f
l
= (7)
0 aT 0
The interaction between a uniformly loaded surface and the unloaded soil surface
(soil-to-soil) can be computed by applying a uniform rectangular load on the surface
to calculate the deflection of the soil at desired positions.
By applying rectangular unit pressures onto the soil surface or a pile head
corresponding to the raft elements, deflections at the centre of each element as a
result of different interactions can be obtained. These deflections form the columns
of the influence matrix for the soil and piles [Isp]. The deflection at the centres of the
elements due to the contact pressures between the raft and piles or soil can be
expressed as
{ sp} = [Isp]{Fsp} (8)
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
where sp are the deflections of the soil at the interface between the raft and soil and
raft and pile, Fsp are the contact pressures at the interface.
Analysis of Raft
The raft is analysed by the use of eight noded isoparametric elements with forty-eight
degrees of freedom. The contact stresses between the raft and soil surface are
assumed to be a uniform rectangular block of pressure acting vertically on each raft
element. The actual displacement at the centre of each raft element (Zhang and Small,
1991) can be computed and expressed as
To obtain the influence matrix of the raft [Ir], the raft must be pinned at one node to
restrain it from undergoing rigid body rotations and translations. Deflections at the
centre of all elements can be calculated by applying a unit pressure to any single
element. These deflections form the columns of the influence matrix of the raft.
where Ptot are the total loads applied to the raft and Mx and My are the total moments
due to the applied load about the pin in the x and y directions respectively. A1 to An
correspond to the areas of each raft element if the raft consists of n elements.
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
[I sp + Ir ] a b c F a
T
0 0 0 Ptot
T
= (11)
0 0 0 x Mx
T
0 0 0 y My
Results
Figure 5 shows the comparison between the results of the present program (finite
layer (FL)) and a finite element (FE) program for displacement of the piled raft along
the centre line (the x-axis in Figure 4). The overall depth of the clay layer is assumed
to be 28 m below the raft. It can be seen that the finite element results are slightly
smaller than the finite layer results. This may be due to several assumptions made in
the finite layer analyses (i) soil layers are of infinite lateral extent, (ii) assumption of
a perfectly smooth base of the raft and pile heads are not rigidly attached to the raft,
(iii) horizontal movements are neglected. Computational times were measured as 4.5
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
mins and 17.8 mins for the present program and the finite element program
respectively.
y
0.3w 0.2w
0.2w
0.3w
w x
Uniform load
L/2
Figure 6 shows the vertical displacements for the case where the piles are of equal
length (Case 1) and longer and larger beneath the heavy load (Case 2). The overall
depth of the clay layer used for the analysis is assumed to be 40 m below the raft. As
was expected, the raft tilts towards the more heavily loaded area of the raft. It is
clearly shown that when the raft is subjected to non-uniform loads, increasing the
pile length and diameter can help to minimize the total and differential displacements,
and the tilt.
0.010
Displacement (m)
0.020
0.030
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
0.040
0.050
0.060
FL FE
Figure 5. Comparison of vertical displacement along the centre line of the raft
(x-axis) with equal pile lengths
X-coord. (m)
0 5 10 15 20 25
0.000
0.010
Displacement (m)
0.020
0.030
0.040
0.050
0.060
CASE 1 CA SE 2
Case History
In order to reduce settlements and excessive tilt of the foundation, a piled-raft was
designed as the foundation for the building such that part of the structural load could
be transferred to the deeper soil layers (Katzenbach et al 2000). The raft is 58.8 m x
58.8 m supported by 64 bored piles. The raft is founded at a depth of 14 m below the
ground level and the thickness varies from 3 m at the edge to 6 m at the centre. The
piles have a diameter of 1.3m and are arranged in three concentric rings. The pile
length varies from 26.9 m in the outer ring through 30.9 m in the middle ring and to
34.9 m in the inner ring. The effective load applied to the raft was 1570 MN. Figure
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
7 shows the layout of the piled raft foundation. The piles in the middle and outer
rings are designed to carry the outer load bearing walls of the building to minimize
tilting and the piles in the inner circles are designed to carry the loads of the centre
core of the building (Franke, 1991).
The piled raft was installed with extensive instrumentation to monitor the behaviour.
These measurements include settlement of the foundation, loads at pile heads and
along the pile shafts and contact pressures on the raft.
Exterior corner
Core column
14 m
6m
34.9 m
58.8 m
28 piles 26.9 m long
20 piles 30.9 m long
16 piles 34.9 m long
The measured settlement at the centre of the raft upon completion in 1990 was 85
mm (Katzenbach et al 2000), increasing to 115 mm in 1995 about 5 years after
completion (Tamaro, 1996), and reaching 144 mm in 1998 (Reul and Randolph,
2003).
Analysis of the piled raft was carried out by the present program and comparisons
were made with the measurements. The thicknesses of Frankfurt Clay and Frankfurt
Limestone for the anaylsis have been assumed to be 74.8 m and 55.2 m respectively
(Reul and Randolph, 2003). Young’s modulus of Frankfurt clay can be calculated
using the following expression based on back-analysis (Reul, 2000):
) z 30 &
E = 45 + tanh' $ + 1 × 0.7 z (12)
( 15 %
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
where E is Young’s modulus (MPa) and z is the depth (m) below the surface. Table 2
summarizes the parameters used for the analysis.
Table 2: Properties of materials used in the finite layer analyses (Reul and Randolph,
2003)
The analysis of the foundation was carried out by the present program based on the
finite layer method. To save computation time, one-quarter of the foundation has
been modeled as shown in Figure 8 for the analysis. At the axes of symmetry, the raft
is restrained from rotations in the appropriate directions. As symmetry was assumed
in the analysis, tilting of the foundation has been neglected.
Figure 9 shows the settlement contours of the foundation. The computed settlement
at the centre of the raft amounts to 139 mm which shows good agreement with the
measurement of 144 mm in 1998 ((Reul and Randolph, 2003). As there is a lack of
support underneath the central region of the foundation, some small differential
settlements were observed from centre to edge in the finite layer analysis.
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
30
12.2 cm
12.4 cm
25
12.6 cm
20 12.8 cm
15
10
-5 0 5 10 15 20 25 30 35
Conclusion
The Finite Layer Method can be used for the analysis of piled rafts supported by
piles of different lengths and diameters and subjected to non-uniform loadings.
Comparisons of actual measurements and the results of finite element calculations
have shown that the finite layer method can provide reasonably accurate solutions. It
is very easy to prepare data for the finite layer analysis as only the pile locations and
properties need to be specified, and the piles discretised into linear elements. This
method is more efficient than the conventional three-dimensional finite element
analysis. The computational time for the finite layer analysis is about four times less
than for the finite element analysis in the example presented here.
References
Chow, Y. K. (1986). “Analysis of vertically loaded pile groups.” Int. Jl. for
Numerical and Analytical Methods in Geomechanics, 10, 59-72.
Downloaded from ascelibrary.org by EGE UNIVERSITY on 12/22/14. Copyright ASCE. For personal use only; all rights reserved.
Franke, E. (1991). “Measurements beneath piled rafts.” ENPC Conf. on Deep Found.,
Paris.
Katzenbach, R., Arslan, U. and Moormann C. (2000). “Chapter 13, Piled raft
foundation projects in Germany.” in Design applications of raft foundations, J. A.
Hemsley, Editor, Thomas Telford, London, pp.323-391
Reul, O. and Randolph, M. F. (2004). “Design Strategies for Piled Rafts Subjected to
Nonuniform Vertical Loading.” Jl. of Geotechnical and Geoenvironmental Eng.,
ASCE, 130(1)1-11
Sinha, J. (1997) Piled raft foundations subjected to swelling and shrinking soils. PhD
thesis, University of Sydney, Australia
Small, J. C. and Zhang, H. H. and Chow, H. (2004). “Behaviour of piled rafts with
piles of different lengths and diameters.” Proceedings of the 9th Australia New
Zealand Conference on Geomechanics, Auckland, Vol. 1,123-129.
Ta, L. D. and Small, J. C. (1996). “Analysis of piled raft systems in layered soils.”
Int. Jl. for Numerical and Analytical Methods in Geomechanics, 20,57-72.
Zhang, H. H. and Small, J. C. (2000). “Analysis of axially and laterally loaded pile
groups embedded in layered soils.” Proceedings of the 8th Australia New Zealand
Conference on Geomechanics, Hobart, Vol. 1, 475-483.
KEYWORDS: