Displacements in a Soil Mass Due to Pile Groups

H.G. Poulos and N.S. Mattes

List of Symbols

ApA_pCross-sectional area of pile.
DDDepth of finite layer.
EpE_pYoung’s modulus of pile material.
EsE_sModulus of elasticity of soil.
EsjE_{sj}Modulus of elasticity of soil layer jj.
EsE’_sModulus of elasticity of soil skeleton.
EuE_uUndrained modulus of elasticity of soil.
HHDepth of point below soil surface.
pI_pIInfluence factor for vertical displacement due to a vertical
point load.
IibI_{ib}Influence factor for displacement at ii due to the base stress Ρb\Rho_b.
IijI_{ij}Influence factor for displacement at ii due to uniform shear stress at jj.
Ii,jI_{i,j}Influence factor for displacement beneath the centre of the pile group for the level of the top of layer jj due to pile ii.
IjI_jDisplacement influence factor on the pile axis at the level of the top of layer jj.
IxjI_{xj}Displacement influence factor II for H/L=0H/L = 0 corresponding to the distance between the centre of pile jj and the point xx.
IxjI’_{xj}Value of II for H=DH = D corresponding to the distance between the centre of pile jj and the point xx.
IρI_{\rho}Influence factor for displacement due to a single pile.
IIInfluence factor for displacement of a pile in a layer of depth D.
KKPile stiffness factor, =(Ep/Es)RA= (E_p/E_s) \cdot R_A.
LLLength of pile.
LeL_eLength of equivalent pier.
PPTotal load on pile.
PavP_{av}Average pile load.
Pi,PjP_i, P_jLoad in pile ii and pile jj respectively.
RRSettlement ratio.
Z,Z1,R,R1,cZ, Z_1, R, R_1 , c Variables in Mindlin’s equation (Fig. 1 (b)).
ddDiameter of pile shaft.
dbd_bDiameter of pile base.
ded_eDiameter of equivalent pier.
i,ji, jIntegers.
mmTotal number of different layers.
nnNumber of pile elements.
ppUniformly distributed shear stress.
pbp_bUniform vertical stress on base.
pjp_jUniform shear stress on element jj.
rrRadial distance.
δ\deltaHeight of pile element =L/n= L/n.
ρ\rhoDisplacement.
ρGD\rho G DDisplacement of the pile group in a layer of depth DD and
modulus of elasticity E1E_1 which is underlain by a rigid base.
ρi\rho iDisplacement at ii due to all elements of the pile.
ρx\rho xDisplacement at point xx.
νs\nu_sPoisson’s ratio of soil.

1.—INTRODUCTION

The settlement behaviour of single piles and pile groups has recently been examined theoretically in a series of papers (Refs. 5, 9, 10 and 11). In these papers, solutions have been obtained for the settlement at the head of a pile or pile group situated in an ideal elastic homogeneous soil mass, and it has been found that such solutions for an ideal soil reproduce observed behaviour of piles in real soils with encouraging accuracy. In many cases in practice however, the soil profile is layered with compressible layers present below the piles, and hence the solutions previously obtained for the displacement of piles in a uniform soil may not be directly applicable. An outstanding example of such a case has been reported by Golder and Osler (Ref. 4) in which a pile foundation supporting a furnace was founded in compact sand which was underlain at a considerable depth by a deep deposit of compressible clay. It was found that the settlement of the piles in the sand was responsible for only about 10% of the final settlement, the remaining 90% being attributable to the underlying clay deposit. In view of the not-infrequent occurrence of such situations, it is desirable to develop a method by which the settlement of compressible soil layers underlying the piles may be determined. In some cases, it may also be of interest to estimate the settlement of the ground surface at some distance away from the piles, for example, in determining the additional settlement of an existing building due to a new structure.

In this paper, the above two problems are considered in detail. Numerical solutions are obtained firstly for the vertical displacement within a soil mass due to a single loaded pile, and using these solutions, methods are described for calculating the surface displacement profile around a pile group and the settlement of soil layers beneath a pile or pile group.

In obtaining the solutions for the displacements due to a single pile, it is assumed that no local yield occurs between the pile and the soil and that conditions within the soil remain purely elastic. Also, in considering the settlement of pile groups, the groups are assumed to be free-standing with no contact between the pile cap and the soil. Analysis has shown that, even if the pile cap is, in fact, in contact with the soil, it has little influence on the settlement of the group, provided the pile spacing is relatively close and elastic conditions prevail within the soil.

2.—ANALYSIS

A pile of length LL, shaft diameter dd, base diameter dbd_b, area ApA_p and Young’s modulus EpE_p , is situated in an ideal elastic homogeneous isotropic soil mass having parameters EsE_s and νs\nu_s. At any point ii within the soil mass, the vertical displacement due to the stress pjp_j on a cylindrical element jj may be obtained by double integration of the Mindlin equation for the vertical displacement due to a vertical point load within a semi-infinite elastic mass. Referring to Fig. 1 (a), the displacement influence factor for the point i due to loading on element jj is

Iij=2(N1)δNδ0πpId2 dθ dz(1)\begin{equation} I_{ij} = 2 \int_{(N-1)_\delta}^{N\delta} {\int^{\pi}_{0}} {_pI} \frac{d} {2} \ d \theta \ dz \label{eq:1} \end{equation} \ref{eq:1} (1)

where

δ=L/n\delta = L/n

pI=_pI = Influence factor for vertical displacement due to a vertical point load.

From Mindlin’s equation,

pI=(1+νs)8π(1νs)Es[z12R13+34νs+8νs2R+(34ν)z22cz+2c2R3+6cz2(zc)R5](2)\begin{equation} _pI = \frac{ (1 + \nu_s) }{ 8\pi(1 – \nu_s) E_s } \left[ \frac{z_1{^2}}{R_1{^3}} + \frac{3 – 4 \nu_s + 8 \nu_s{^2}}{R} + \frac{(3-4\nu) z^2 – 2cz + 2c^2}{R^3} + \frac{6cz^2(z-c)}{R^5} \right] \label{eq:2} \end{equation} \ref{eq:2} (2)

where z1,z,R1,R and cz_1, z, R_1, R \text { and} \ c are geometrical parameters defined in Fig. 1 (b).

Fig. 1.

As described by Poulos and Davis (Ref. 10), the integration with respect to zz is done analytically, while the integration with respect to θ\theta is done numerically.

Similarly, the influence factor for vertical displacement at ii due to the uniform base stress pbp_b is given by

Iib=02π0db/2pI r dr dθ(3)\begin{equation} I_{ib} = \int^{2 \pi}_0 {\int^{d_{b/2}}_0} {_pI} \ r \ dr \ d \theta \label{eq:3} \end{equation} \ref{eq:3} (3)

where pI_pI is given by Eq. (2).

The integration is done analytically with respect to r and numerically with respect to θ\theta.

The displacement at point ii due to all elements of the pile is then given by

ρi=j=1nIijpj+Iibpb(4)\begin{equation} \rho i = \overset{n}{ \underset{j = 1}{\sum}} I_{ij} p_j + I_{ib} p_b \label{eq:4} \end{equation} \ref{eq:4} (4)

Values of Pi and Pb have been obtained for various values of Lid by Poulos and Davis (Ref. 10) for incompressible piles, and Mattes and Poulos (Ref. 5) for compressible piles. In the latter paper, the relative compressibility of the pile is defined by the pile stiffness factor K=(Ep/Es)RAK = (E_p/E_s) \cdot R_A , where RA=Ap/(π d2/4)R_A = A_p/(\pi\ d^2/4). The larger KK, the less compressible is the pile. These values have been used with the values of Iii and lib calculated from Eqs. (1) and (3) to obtain displacement distributions for a wide range of cases. In evaluating IijI_{ij} and IibI_{ib}, intervals of π/20\pi / 20 have been used for the numerical integrations. In all cases computed, dbd_b has been taken equal to dd.

3.—DISPLACEMENT DISTRIBUTIONS DUE TO A SINGLE LOADED PILE

Influence factors IpI_p for the displacement due to a single loaded pile have been obtained for values of L/dL/d ranging from 1 to 100 and are shown in Figs. 2 to 7 for various dimensionless depths H/LH/L and radii r/Lr/L, and for νs=0.5\nu_s = 0.5. The actual displacement ρ\rho is given by

ρ=PL EsIρ(5)\begin{equation} \rho = \frac{P}{L \ E_s} \cdot I_\rho \label{eq:5} \end{equation} \ref{eq:5} (5)

In Figs. 2 to 6, K=1000K = 1000, while in Fig. 7, values of IρI_\rho are given for K=20000K = 20 000, for L/d=100L/d = 100. KK has a considerable influence on the displacements near the pile for L/d>10L/d \gt 10, IpI_p increasing as KK decreases, but only a small influence for L/d<10L/d \lt 10 or at points remote from the pile. The influence of νs\nu_s on IρI_\rho is relatively small, especially for L/d>10L/d \gt 10, so that for practical purposes, the values of IρI_\rho for νs=0.5\nu_s = 0.5 may be used for all values of νs\nu_s.

Fig. 2.—Displacement Influence Factors.
Fig. 3.—Displacement Influence Factors.
Fig. 4.—Displacement Influence Factors.
Fig. 5.—Displacement Influence Factors.
Fig. 6.—Displacement Influence Factors.
Fig. 7.—Displacement Influence Factors.

For points remote from the axis, it is foUnd that the pile may be replaced by a point load of magnitude equal to the load on the pile, and acting at a depth of 2/3L below the surface. The Mindlin equation for vertical displacement due to a point load may then be used, viz.,

ρ=PpI(6)\begin{equation} \rho = P \cdot _pI \label{eq:6} \end{equation} \ref{eq:6} (6)

where PP is the total load on the pile, and pI_pI is defined in Eq. (2).

On the soil surface (H=0H = 0), this equation reduces to the following simple form:

ρ=P(1+νs)2π Es[c2R3+2(1νs)R](7)\begin{equation} \rho = \frac{ P(1 + \nu_s) }{ 2 \pi \ E_s } \left[ \frac{c^2}{R^3} + \frac{ 2(1 – \nu_s) }{ R } \right] \label{eq:7} \end{equation} \ref{eq:7}(7)

where R=(c2+r2)12R = (c^2 + r^2)^{\frac{1}{2}}.

It is found that for r/L0.4r/L \geqslant 0.4, the surface displacement calculated from Eq. (7) gives a value which agrees to within ± 3% of the more correct value calculated from Eq. (4), provided that L/d10L/d \geqslant 10. Even for shallow piers (e.g., L/d=1L/d = 1), the use of Eq. (7) for surface displacements gives results almost identical with the more correct analysis for r/L0.75r/L \geqslant 0.75.

4.—DETERMINATION OF SURFACE DISPLACEMENT PROFILE DUE TO A PILE GROUP

The displacement distributions for a single pile in Figs. 2 to 7 may be used to determine the surface displacement profiles due to a loaded pile group by superposition of the displacements due to each individual pile in the group at the point in question. Thus, for a group of n piles situated in a deep homogeneous soil layer, the displacement ρx\rho_x at a point XX on the soil surface is

ρx=1L Es j=1n Pj Ixj(8)\begin{equation} \rho_x = \frac{1}{L \ E_s} \ \overset{n}{\underset{j = 1}{\sum}} \ P_j \ I_{xj} \label{eq:8} \end{equation} \ref{eq:8}(8)

where

PjP_j is the load on pile jj,

IxjI_{xj} is the displacement influence factor IρI_\rho, for H/L=0H/L = 0, corresponding to the distance between the centre of pile jj and the point XX.

It should be noted that in this method, no account can be taken of the reinforcing effect of any piles which may be present between the point XX and the pile jj. The presence of such piles will tend to decrease the influence of pile jj on the displacement at XX and hence the displacement will be over-estimated by using Eq. (8). Some indication of the extent to which the displacement may be overestimated may be inferred from Fig. 8 (see later) which suggests that the overestimate will be relatively small.

If the pile group is situated on a finite layer of depth DD, which is underlain by a rigid base, the Steinbrenner approximation (Ref. 13) may be used to estimate the displacements. In this case,

ρx=1L Es j=1n Pj(IxjIxj)(9)\begin{equation} \rho_x = \frac{1}{L \ E_s} \ \overset{n}{\underset{j=1}{\sum}} \ P_j (I_{xj} – I’_{xj}) \label{eq:9} \end{equation} \ref{eq:9}(9)

where IIxjI’_{xj} is the value of IρI_\rho for H=DH = D and for the value of rr equal to the distance between the centre of pile jj and the point XX.

The work of Davis and Taylor (Ref. 2) and Poulos (Ref. 8) in connection with surface footings suggests that the above approximation will be reasonably accurate in the vicinity of the pile but may be inaccurate remote from the pile (e.g., for r/D>1r/D \gt 1) or for very shallow layers.

If the pile group has a rigid cap, the distribution of loads Pi within the group will generally be non-uniform. This load distribution may be calculated by the method described by Poulos (Ref. 9) and Poulos and Mattes (Ref. 12). For square groups of incompressible piles, values of Pi are tabulated by Poulos (Ref. 9) for 323^2, 424^2 and 525^2 groups for L/d=25L/d = 25, while some load distributions for groups of compressible floating piles are presented by Poulos and Mattes (Ref. 12).

In cases in which the load distribution within a group is not available or where a rough estimate only of the surface displacement is required, it is more convenient to assume a uniform distribution of pile loads (i.e., Pj=P/nP_j = P/n) or else, if the group is large, to consider it as an equivalent single pier. In the latter case, the diameter of the pier should be taken such that the cross-sectional area of the equivalent pier equals the gross plan area of the group while the equivalent length of the pier may generally be taken in the range 0.6L to 0.9L (Ref. 9).

For a square group of 9 piles with a centre-to-centre spacing of 5 dia., Fig. 8 shows a typical comparison between the surface displacement profile for the group with a rigid cap, calculated from Eq. (9), that calculated by assuming the pile loads in the group are equal (i.e., that the pile cap is perfectly flexible) and that calculated by assuming the group to have a rigid cap and to be a single equivalent cylindrical pier of equivalent length Le=0.8LL_e = 0.8L. It will be seen that there is close correspondence between all three displacement profiles away from the immediate vicinity of the group. Underneath the pile group, the displacement distribution for the case of equal loads in all piles is not uniform and hence there is some difference between this case and the other two curves for a rigid pile cap. However, the differences are localised and for x/L0.3x/L \geqslant 0.3, the surface displacement profile for equal pile loads is virtually identical with that for the case of a rigid pile cap. Remote from the group, there is, however, a small discrepancy between the profile for the equivalent single pier and the profiles calculated for the actual group, the displacements in the latter case being greater than those for the equivalent pier. This discrepancy may possibly reflect the small order of overestimation of displacement involved in neglecting the reinforcing effect of intervening piles when using Eq. (9).

Fig. 8.—Comparison between Correct and Approximate Surface Displacement Profiiles.
Fig. 9.—Influence of K on Surface Displacement Profiles.

Typical surface displacement profiles are shown in Fig. 9 for a square four-pile group (222^2 group) in a semi-infinite mass. The increased displacements with decreasing values of KK are clearly shown. However, for x/L>1x/L \gt 1, the surface displacements are almost independent of KK.

The influence of the layer depth on the surface displacement profile is shown in Fig. 10 for a 222^2 group. As the layer depth decreases, the surface displacements are “damped” and, for relatively shallow layers (e.g., D/L=1.5D/L = 1.5), may be negligible at distances of one pile length or more from the centre of the group.

Fig. 10.—Influence of Layer Depth on Surface Displacement Profiles.

5.—DETERMINATION OF SETTLEMENT OF PILES AND PILE GROUPS UNDERLAIN BY LAYERED STRATA

The method outlined in the previous section for calculating surface displacement may be extended to the case where a pile or pile group is founded in a homogeneous layer which is underlain by further layers of different soils. In this case, the settlement may be calculated approximately by adding to the settlement of the pile or pile group in the founding layer, the contributions of the underlying layers.

Considering first the case of a single pile founded in a uniform layer of depth DD, a convenient method of calculating the settlement of the pile is as follows:

ρ=P I1L Es1+PL{j=2m1(IjIj+1)Esj+ImEsm}(10)\begin{equation} \rho = \frac{P\ I_1}{L\ E_{s^1}} + \frac{P}{L} \left\{ \overset{m – 1}{ \underset{j = 2}{\sum} } \frac{(I_j – I_{j+1})}{E_{sj}} + \frac{I_m}{E_{sm}} \right\} \label{eq:10} \end{equation} \ref{eq:10}(10)

where

I1=I_1 = displacement influence factor for a pile in a layer of depth DD (obtained from Poulos and Davis (Ref. 10) for incompressible piles, Mattes and Poulos (Ref. 5) for compressible piles or Poulos and Mattes (Ref. 11) for end-bearing piles),
Ij=I_j = displacement influence factor on the pile axis at the level of the top of layer jj (obtained from Figs. 2 to 7),
Esj=E_{sj} = Young’s modulus of layer jj,
m=m = total number of different layers (including the founding layer).

The above method of calculation is similar to that suggested by Nair (Ref. 7) for piles and that used by Egorov et al (Ref. 3) for surface foundations and, in effect, makes use of the Steinbrenner approximation for the calculation of the displacements of the underlying layers.

Before proceeding, it is useful to consider briefly the limitations of Eq. (10) to calculate settlements. The major limitation is that no account is taken of the influence of changes in stress distribution below the pile which may occur because of the layering of the soil profile. In particular, if a stiffer layer overlies a softer layer, some concentration of stress into the stiffer layer will occur. Although the influence of such stress concentration on pile behaviour cannot readily be determined at the present time, some idea of the errors involved in using Eq. (10) may be obtained by considering the case of a circular surface footing on a two-layer elastic system underlain by a rough rigid base. It is found that the displacements calculated from Eq. (10) are generally greater than the more correct values, (obtained either from Burmister (Ref. 1) or from Odemark’s method, quoted by Ueshita and Meyerhof (Ref. 14». The errors are greatest when a stiff layer of considerable thickness overlies a thick softer layer, but unless the ratio of the Young’s moduli of the upper and lower layers IS greater than about 10, the resulting overestimate is. not serious. Furthermore, if softer layers overlie stiffer layers, Eq. (10) gives a displacement which agrees closely with the correct value.

A further limitation of Eq. (10) is that it inherently assumes that the shear stress distribution along the pile is that for a pile in a semi-infinite homogeneous soil mass and is unaffected by any heterogeneity of the soil below the pile. For cases involving different soil layers at some depth below the pile tip, the above assumption should be satisfactory as it has been found by Poulos and Davis (Ref. 10) that the shear stress distribution remains substantially unchanged if a rigid base exists below the tip of the pile. However, the approximation may be in error if the pile bears directly on to a stratum of considerably greater stiffness than the upper layer. The distribution of shear stress along the pile in this case may be considerably different from that for a pile in a semi-infinite mass, especially if the pile is relatively incompressible (Poulos and Mattes, Ref. 11). In such a case, the displacement factor for the pile in the upper layer of depth HH should be calculated from the solutions for an end-bearing pile presented by Poulos and Mattes (Ref. 11). The displacement factors for the underlying layers are calculated as before.

It therefore appears that Eq. (10) will give a reasonable estimate of the settlement of the pile unless a very stiff layer of considerable thickness beneath the pile overlies a very soft layer. The error involved in using Eq. (10) will generally lead to an overestimate of the settlement.

The calculation of the settlement of a single pile from Eq. (10) may be extended to the case of a pile group. For a group of n piles situated in a homogeneous layer of depth DD, underlain by a series of other strata, the displacement at the centre of the pile group can be calculated by summation of the settlement of the underlying layers due to all piles in the group. Thus

ρ=ρGD+i=1nPiLj=2m1{(Ii,jIi,j+1)Esj+Ii,mEsm}(11)\begin{equation} \rho = \rho_{GD} + \overset{n}{\underset{i=1}{\sum}} \frac{P_i}{L} \overset{m-1}{\underset{j=2}{\sum}} \left\{ \frac{ (I_{i, j} – I_{i, j+1}) }{E_{sj}} + \frac{I_{i, m}}{E_{sm}} \right\} \label{eq:11} \end{equation} \ref{eq:11}(11)

where
ρGD=\rho_{GD} = displacement of the pile group in a layer of depth DD and Young’s modulus E1E_1 which is underlain by a rigid base,
Pi=P_i = load in pile ii,
Ii,j=I_{i,j} = displacement influence factor beneath the centre of the pile group for the level of the top of layer jj due to pile ii (obtained from Figs. 2 to 7).

The value of ρGD\rho_{GD} may be calculated as

ρGD=Rρ1(12)\begin{equation} \rho_{GD} = R \cdot \rho_1 \label{eq:12} \end{equation} \ref{eq:12}(12)

where
Rs=R_s = the settlement ratio of group settlement to the settlement of a single pile under the same average load as a pile in the group,
ρ1=\rho_1 = settlement of single pile at the average pile load.

For a floating pile group, ρ1\rho_1 may either be calculated from the influence factors presented by Poulos and Davis (Ref. 10) or Mattes and Poulos (Ref. 5), or from a field load test on a single pile. RsR_s may be calculated by the method described by Poulos (Ref. 9). For a group which rests on a firmer stratum below which a softer layer exists, ρ1\rho_1 may be calculated from the influence factors presented by Poulos and Mattes (Ref. 11), or again from a field load test, while RsR_s may be calculated, as if the group was end-bearing, by the method described by Poulos and Mattes (Ref. 12).

For a pile group with a rigid cap, the values of PiP_i in the group will not be uniform, as previously mentioned. The load distributions obtained by Poulos (Ref. 9) and Poulos and Mattes (Ref. 12) for a pile group in a homogeneous mass may be used to obtain values of PiP_i although some slight inaccuracy may result because of the influence of the underlying compressible layers in modifying the load distribution. However, as shown in the previous section, little error will ensue if the load distribution within the group is assumed uniform. Alternatively, and again more expediently, the pile group may be replaced by an equivalent single pier as previously outlined. A further advantage of this latter method is that, for the relatively small value of L/dL/d which the equivalent pier will have, the values of IρI_\rho will be almost independent of KK and thus it will generally not be necessary to consider different values of KK for different layers.

An example is given below of the calculation of the settlement of a floating pile group using both the correct . load distribution and a uniform load distribution in the group and also replacing the group by an equivalent single pier. An example is also given of the settlement calculation for a group which rests on a firmer stratum below which a softer layer exists.

6.—ILLUSTRATIVE EXAMPLES

6.1 Floating Pile Group:

The example considered is shown in Fig. 11. A 323^2 group of concrete piles with a rigid cap is founded in a soil layer which is underlain by a further three layers and it is desired to calculate the settlement of the pile group. The values of Young’s modulus EsE_s of the soil layers are given in Fig. 11 and are taken to be average values of the Young’s modulus of the soil skeleton so that the settlements calculated will be total final settlements.

Fig. 11.—Illustrative Example No.1.

The settlement of the pile group in layer 1 may be calculated from Eq. (12). For the piles in layer 1, assuming EE of the concrete to be 3x1063 x 10^6 lb./sq. in., K=(3×106)/(5×102)=6,000K = (3 \times 10^6)/(5 \times 10^2) = 6,000. For this value of KK, the piles may be considered incompressible (Mattes and Poulos, Ref. 5). From the influence factors given by Poulos and Davis (Ref. 10), for L/d=25L/d = 25 and H/L=1.5H/L = 1.5, I1=1.38I_1 = 1.38. For the pile spacing of s=5ds = 5d, the value of RsR_s is estimated from the results of Poulos (Ref. 9) to be 2.3. Hence the settlement of the pile group in layer 1=(PG/n)(Rs/LEs1)I1=3.78 in.1 = (P_G/n) \cdot (R_s / LE_{s1}) \cdot I_1 = 3.78\ \text{in.}, where PGP_G is the total load on the pile group.

In order to calculate the additional settlements due to the underlying layers, three alternative methods will be used. Firstly, the correct non-uniform load’ distribution in the group will be considered; secondly, it will be assumed that the load distribution in the group is uniform; thirdly, the group will be replaced by an equivalent single pier.

(i) Considering Non-Uniform Pile Loads

From the tabulated results presented by Poulos (Ref. 9), the load distribution in the group may be obtained. The average pile load Pav=720/9=80 tonsP_{av} = 720/9 = 80\ \text{tons}, and using the key to the piles shown in Fig. 11, the pile loads are as follows:

P1=1.32 Pav=106 tonsP2=0.84 Pav=67 tonsP3=0.34 Pav=28 tons\begin{equation} \begin{array}{l l l} P_1 & = 1.32\ P_{av} = & 106 \text{ tons} \\ P_2 & = 0.84\ P_{av} = & 67 \text{ tons} \\ P_3 & = 0.34\ P_{av} = & 28 \text{ tons} \\ \end{array} \end{equation}

Since the settlement beneath the centre of the group is required, only three piles of the group need be considered because of symmetry. Because of the relatively large values of H/LH/L involved in this problem, the displacement influence factors are unaffected by the value of KK. The calculations are shown in Table I. The total additional settlement due to the underlying layers is found to be 2.25 in. The overall settlement of the group is therefore 3.78 + 2.25 = 6.03 in.

(ii) Considering Uniform Pile Load Distribution

In this case, the load PiP_i on each pile is 80 tons. From Table I, the values of (IjIj+1)/Esj(I_j – I_{j+1})/E_{sj} may be extracted for each pile and each layer. The resulting calculations are shown in Table II. In this case, the settlement of the group is 3.78 + 2.29 = 6.07 in. It will be seen that the calculated settlement is almost identical with that calculated from Table I using the correct non-uniform pile load distribution.

(iii) Considering the Group as an Equivalent Single Pier

The gross plan area of the group = 11 x 11 = 121 sq. ft. Therefore the equivalent diameter of a single pier de = 12.4 ft. From Poulos (Ref. 9), the equivalent length Le of the pier may be taken as 0.8L = 20 ft.

 Le/de=20/12.4=1.6\begin{equation} \begin{array}{l r} \therefore &\ L_e / d_e = 20/12.4 = 1.6 \end{array} \end{equation}

Eq. (10) may now be used to determine the settlement of the underlying layers. Since values of Ip are only given for L/d = 1 and L/d = 2, the values of I for the required value of Lelde = 1.6 will be interpolated from the value/for L/d = 1 and 2. The calculations are shown in Table III. Because the settlements are calculated beneath the centre of the group, only values of I p for r / L = 0 are required. The settlement of the group in this case is 3.78 + 2.52 = 6.30 in.

The calculated additional settlement due to the underlying layers from Table III is about 10% greater than the values calculated by the two previous methods. Nevertheless, the agreement is sufficiently good to suggest that, for a group containing a large number of piles, the replacement of the group by an equivalent single pier provides a convenient and reasonably accurate method of calculating the settlement of the group.

TABLE I

Pile No. iiPile Load PiP_i (tons)riL\frac{r_i}{L}Layer 2
H2L\frac{H_2}{L}I2I_2H3L\frac{H_3}{L}I3I_3P1(I2I3)Es2(in.2)\frac{P_1(I_2 – I_3)}{E_{s2}} \\ (\text{in.}^2)
11060.281.50.40020.28524.6
2670.201.50.41020.28515.7
32801.50.43520.2857.3
Pile No. iiPile Load PiP_i (tons)riL\frac{r_i}{L}Layer 3
H3L\frac{H_3}{L}I3I_3H4L\frac{H_4}{L}I4I_4Pi(I3I4)Es3(in.2)\frac{P_i(I_3 – I_4)}{E_{s3}} \\ (\text{in.}^2)
11060.2820.28530.17574.5
2670.2020.28530.18045.0
328020.29530.18519.7
Pile No. iiPile Load PiP_i (tons)Layer 4
H4L\frac{H_4}{L}I4I_4Pi(I4)Es4(in.2)\frac{P_i(I_4)}{E_{s4}} \\ (\text{in.}^2)
110630.1752.1
26730.1801.4
32830.1850.6
Settlement of each layer = 4 X (settlement due to pile 1 + settlement due to pile 2 + settlement due to pile 3).
For layer 2, Settlement =125×12 [4×(24.6×15.7)×7.3]=0.56 in.= \frac{1}{25 \times 12}\ [ 4 \times (24.6 \times 15.7) \times 7.3 ] = 0.56 \text{ in.}
For layer 3, Settlement =125×12 [4×(74.5×45.0)×19.7]=1.64 in.= \frac{1}{25 \times 12}\ [ 4 \times (74.5 \times 45.0) \times 19.7 ] = 1.64 \text{ in.}
For layer 4, Settlement =125×12 [4×(2.1×1.4)×0.6]=0.05 in.= \frac{1}{25 \times 12}\ [ 4 \times (2.1 \times 1.4) \times 0.6 ] = 0.05 \text{ in.}
Total settlement of layers 2 to 4 =2.25 in.= 2.25 \text{ in.}

TABLE II

Pile No. iiValues of Pi (IjIj+1)Esj\frac{P_i\ (I_j – I_{j+1})}{E_{sj}}
Layer 2Layer 3Layer 4
118.856.21.6
218.854.71.7
320.856.31.7
For layer 2, settlement =125×12 [4×(18.8+18.8)+20.8]=0.57 in.= \frac{1}{25 \times 12}\ [ 4 \times (18.8 + 18.8) + 20.8 ] = 0.57 \text{ in.}
For layer 3, settlement =125×12 [4×(56.2+54.7)+56.3]=1.67 in.= \frac{1}{25 \times 12}\ [ 4 \times (56.2 + 54.7) + 56.3 ] = 1.67 \text{ in.}
For layer 4, settlement =125×12 [4×(1.6+1.7)+1.7]=0.05 in.= \frac{1}{25 \times 12}\ [ 4 \times (1.6 + 1.7) + 1.7 ] = 0.05 \text{ in.}
Total settlement of layers 2 to 4 =2.29 in.= 2.29 \text{ in.}

TABLE III

L/dL/dLayer 2
H2Le\frac{H_2}{L_e}I2I_2H3Le\frac{H_3}{L_e}I3I_3(I2I3)Es2(sq.in./ton)\frac{(I_2 – I_3)}{E_{s2}} \\ (\text{sq.in./ton})
11.880.1862.50.136
21.880.4862.50.300
1.6
(by linear interpolation)
1.880.3552.50.2350.223
L/dL/dLayer 3
H3Le\frac{H_3}{L_e}I3I_3H4Le\frac{H_4}{L_e}I4I_4(I2I3)Es3(sq.in./ton)\frac{(I_2 – I_3)}{E_{s3}} \\ (\text{sq.in./ton})
12.50.1363.750.082
22.50.3003.750.180
1.6
(by linear interpolation)
2.50.2353.750.1410.601
L/dL/dLayer 4
H4Le\frac{H_4}{L_e}I4I_4I4Es4 (sq.in./ton)\frac{I_4}{E_{s4}} \\ \text{ (sq.in./ton)}
13.750.082
23.750.180
1.6
(by linear interpolation)
3.750.1410.16
For layer 2, settlement =P12Le(I2I3)Es2=72020×12×0.223=0.67 in.= \frac{P \frac{1}{2}}{L_e} \cdot \frac{(I_2 – I_3)}{E_{s2}} = \frac{720}{20 \times 12} \times 0.223 = 0.67 \text{ in.}
For layer 3, settlement =P12Le(I2I3)Es2=72020×12×0.601=1.80 in.= \frac{P \frac{1}{2}}{L_e} \cdot \frac{(I_2 – I_3)}{E_{s2}} = \frac{720}{20 \times 12} \times 0.601 = 1.80 \text{ in.}
For layer 4, settlement =P12Le(I2I3)Es2=72020×12×0.016=0.05 in.= \frac{P \frac{1}{2}}{L_e} \cdot \frac{(I_2 – I_3)}{E_{s2}} = \frac{720}{20 \times 12} \times 0.016 = 0.05 \text{ in.}
Total settlement of layers 2 to 4 =2.52 in.= 2.52 \text{ in.}

The above example emphasizes that compressible layers well beneath the pile group may contribute significantly to the settlement of the group. In this case, the additional settlement due to the underlying layers comprises about 37% of the overall settlement. The example also highlights the possible consequence of incomplete site investigation.

6.2 End-Bearing Pile Group:

The problem is shown in Fig. 12. A 323^2 group of concrete piles is founded through alluvial clay on to a medium dense sand, which is underlain by a soft clay layer which rests on shale. Average values of both undrained Young’s modulus EuE_u and Young’s modulus of the soil skeleton EsE’_s are given and both the immediate and total final settlements are calculated. Because the displacement influence factors IρI_\rho are not seriously influenced by Poisson’s ratio νs\nu_s (see Section 3), it has been assumed that, under both undrained and drained conditions, IρI_\rho values are those for νs=0.5\nu_s = 0.5.

The group is considered first as end-bearing. The value of KK of the piles is about 6,000, so that from the results of Poulos and Mattes (Ref. 11), such piles are virtually incompressible, and hence the settlement of a single pile is merely the elastic settlement of the pile itself. Furthermore, from the results of Poulos and Mattes (Ref. 12), no interaction occurs between adjacent piles, i.e., Rs=1R_s = 1, and all piles are equally loaded**.

**Because the supporting sand layer is not perfectly rigid, some interaction may in fact occur, but it is unlikely to be significant.

Hence, the settlement of the piles, when considered as bearing on a rigid stratum, is

S1=P LAp Ep\begin{equation} S_1 = \frac{P\ L}{A_p\ E_p} \end{equation}

where

PP is the average pile load,
ApAp is the area of the pile,
EpEp is Young’s modulus of the pile.

In this case, it is found that, for Ep=3×106 Ib./sq. in.,S1=0.16 in.E_p = 3 \times 10^6 \text{ Ib./sq. in.}, S_1 = 0.16\text{ in.}

Fig. 12.—Illustrative Example No.2.

The calculation of the additional settlements due to the underlying layers is shown in Table IV. The settlement of the sand layer is assumed to occur immediately. The immediate settlement of the group is the sum of the elastic pile settlement and the immediate settlement of the underlying layers, i.e., 0.16 + 3.21 = 3.37 in. The total final settlement is 0.16 + 4.46 = 4.62 in. Thus, because of the underlying soft. clay layer, a considerable amount of consolidation settlement (1.25 in.) occurs in this case, whereas if the pile were founded on a truly rigid stratum, only an immediate settlement of 0.16 in. would occur.

7.—COMPARISON WITH FIELD CASE

As an indication of the applicability of the method described to field situations, an attempt was made to compare the measured settlements of the furnace foundation described by Golder and Osler (Ref. 4) with the calculated values. The foundation supported a total load of 2,400 tons and consisted of 32 Franki piles of 22-ft. length arranged in a generally rectangular pattern consisting of 3 interior rows of 8 piles each, the rows spaced at 10.5 ft. and the piles at 6 ft., and at each end, 4 piles spaced at 7-ft. centres. The soil profile could generally be characterised as about 95 ft. of sand interspersed with layers of silty clay and sandy clay, and underlain by a very deep deposit of a firm grey silty clay (Leda clay).

A single pile test was carried out and revealed that, under the average working load, the single pile settlement was about 0.04 in. An analysis carried out, as described by Poulos (Ref. 9), gave a value of settlement ratio RsR_s of the pile group in the sand of 8.5. Thus, the settlement of the piles in the sand (ρGD\rho_{GD}) is 8.5 x 0.04 = 0.34 in.

Considering now the settlement due to the Leda clay, only oedometer results are given for this clay and it is therefore necessary to use these results to estimate Young’s modulus EsE’_s of the clay. The data given is compression index Cc=0.33C_c = 0.33, and void ratio = 1.25. From these results, a value of the coefficient of volume decrease mvm_v may be derived for a typical depth of say 50 ft. below the sand-clay interface. This gives mv=6.3×106 ft.2/lb.m_v = 6.3 \times 10^{-6} \text{ ft.}^2 /\text{lb.} Using now the theory for an ideal elastic soil,

Es=(1+νs)(12νs)(1νs)mv\begin{equation} E’_s = \frac{(1+\nu’_s)(1-2\nu’_s)}{(1-\nu’_s) \cdot m_v} \end{equation}

Assuming a typical value of νs\nu’_s for the firm clay of 0.2, it is found that Es=1,000 lb./sq. in.E’_s = 1,000 \text{ lb.}/\text{sq. in.} Considering now the calculation of the settlement due to the clay, and assuming the clay to be semi-infinite, it is found that the settlement is 3.3 in. Thus, the total settlement of the group is 3.3 + 0.34 = 3.7 in.

The actual foundation was in fact demolished before the final settlement had been attained, at which time, the measured settlement had reached about 1.8 in. Measured values of the coefficient of consolidation cv vary too widely to be able to accurately estimate the degree of settlement at the time of demolition. However, if an average value of cv of 1×102 in.2/min1 \times 10^{-2} \text{ in.}^2/ \text{min}. is assumed, the dimensionless time factor Tv=(cvt)/a2T_v=(c_vt)/a^2 (assuming the foundation to be replaced by an equivalent circular foundation of radius aa) after 11 years is 0.91. From the three-dimensional rate of settlement analysis of Mandel (Ref. 6) for a circular foundation on a sand layer overlying a deep clay layer, it is found that the appropriate degree of settlement is about 0.38. Thus, the estimated final settlement (taking account of an initial sand settlement of 0.5 in.) = (1.3/0.38) + 0.5 = 3.9 in., which is in general agreement with the calculated value. In addition, the initial settlement of the group, of about 0.5 in., due largely to the sand, is consistent with the calculated value of ρGD\rho_{GD} of 0.34 in.

Although the settlements calculated by the method in this paper are very similar to those produced by Golder and Osler using one-dimensional theory, the present method would appear somewhat more logical as it does not require any arbitrary selection of rigid boundaries beneath the clay or reductions for foundation rigidity as is the case with the one-dimensional predictions.

8.—CONCLUSIONS

Distributions of vertical displacement within a semi-infinite mass due to a loaded pile have been obtained for various length-to-diameter ratios of the pile. It has been found that, even at points relatively close to the pile, the displacements due to a pile agree closely with those due to an equivalent single point load acting at a depth of two-thirds of the pile length. A method of calculating the surface displacement profile around a pile group and the settlement of a pile group underlain by layered strata, using the single pile displacement distributions, has been suggested. It has been found that, for a pile group having a rigid pile cap, the displacements may be calculated with sufficient accuracy by assuming all piles are equally-loaded and ignoring the correct non-uniform load distribution within the group. The replacement of the pile group by an equivalent single pier has been shown to lead to an estimate of the displacement which is quite adequate for practical purposes. Such a procedure appears to be most convenient when considering displacements due to a group containing 3 large number of piles.

ACKNOWLEDGMENTS

The work described in this paper forms part of a general programme of research into the settlement of all types of foundations, under the general direction of Professor E. H. Davis, whose comments are gratefully acknowledged. The work was carried out in the School of Civil Engineering at The University of Sydney, the Head of which is Professor J. W. Roderick, and was supported by a grant from the Australian Research Grants Committee. The computational work was carried out using the facilities of the Basser Computing Department of the School of Physics at the University of Sydney.

References

  1. BURMISTER, D. M.-The Theory of Stresses and Displacements in Layered Systems and Applications to the Design of Airport Runways. Proc. Highway Research Board, Vol. 23, 1943, pp. 126-48.
  2. DAVIS, E. H. and TAYLOR, H.-The Movement of Bridge Approaches and Abutments on Soft Foundation Soils. Proc. First Conf. Aust. Road Research Board, Canberra, 1962, Vol. I, Part 2, pp. 740-64.
  3. EGOROV, K. E., KUZMIN, P. G. and POPov, B. P.-The Observed Settlements of Buildings ‘as Compared with Preliminary Calculation. Proc. Fourth Int. Conf. Soil Mechanics and Foundation Engg., London, 12-24 August, 1957. London, Butterworth, 1957, Vol. I, pp. 291-6. 35
  4. GOLDER, H. Q. and OSLER, J. C.-Settlement of a Furnace Foundation, Sorel, Quebec. Canadian Geotech. Jour., Vol. 5, No. I, Feb., 1968, pp.46-56.
  5. MATTES, N. S. and POULOS, H. G.-The Settlement of a Single Compressible Pile. Proc. A.S.C.E., Jour, Soil Mechanics and Foundations, Div., Vol. 95, No. SM1, Jan., 1969, pp. 189-207.
  6. MANDEL, J.-Tassements produits par la consolidation d’une couche d’argile de grande epaisseur. Proc. Fifth Int. Conf. Soil Mechanics and Foundation Engg., Paris, 17-22 July, 1961, Vol. 1, pp. 733-6.
  7. NAIR, K.-Load Settlement and Load Transfer Characteristics of a Friction Pile Subject to a Vertical Load. Proc. Third Pan-American Conf. Soil Mechanics, Caracas, Venezuela, 1967, p. 565.
  8. POULOS, H. G.-Stresses and Displacements in an Elastic Layer Underlain by a Rough Rigid Base. Geotechnique, Vol. 17, No.4, 1967, pp. 378-410.
  9. POULOS, H. G.-Analysis of the Settlement of Pile Groups. Geotechnique, Vol. 18, No.4, 1968, pp. 449-71.
  10. POULOS, H. G. and DAVIS, E. H.-The Settlement Behaviour of Single Axially-Loaded Incompressible Piles and Piers. Geotechnique, Vol. 18, No.3, 1968, pp. 351-71.
  11. POULOS, H. G. and MATTES, N. S.-The Behaviour of Axially Loaded End-Bearing Piles. Geotechnique, Vol. 19, No.2, 1969, pp. 285-300.
  12. POULOS, H. G. and MATTES, N. S.-Settlement and Load Distribution Analysis of Pile Groups. Aust. Geomechanics Jour., Vol. G1. No.1, 1971, pp. 18-28.
  13. STEINBRENNER, W.-Tafeln zur Setzungberechnung. Die Strasse, Vol. 1, 1934, p. 12l.
  14. UESHITA, K. and MEYERHOF, G. G.-Deflection of Multilayer Soil Systems. Proc. A.S.C.E., Jour. Soil Mechanics and Foundations Div., Vol. 93, No. SM5, Sept., 1967, pp. 257-82.