Acessibilidade / Reportar erro

Proposal for considering the group effect in the prediction of settlements in pile groups through load transfer methods

Abstract

When designing a foundation project, it is necessary to ensure that all the elements meet both ultimate and serviceability limit states, which call for predictions of settlement and load capacity. The load transfer methods are a widely used alternative to estimate the load-settlement ratio of piles in the design of foundation projects. However, traditional load transfer methods do not consider the interactive effects between the elements in pile groups. This study proposes changes to the load transfer curves developed by Bohn et al. (2016), aiming to incorporate the group effect in the analysis of load-settlement relationships in pile groups. Comparisons between the predicted settlements obtained using the proposed method and the results of load tests performed by Dai et al. (2012) in Jiangsu, China, showed that the modifications proposed in this study agreed well with the experimental results for most of the analyzed groups.

Keywords
Pile settlements; Load transfer methods; Group effect; Hyperbolic model

1. Introduction

Settlement prediction is an important task that can be performed using different techniques, such as the load transfer methods, which use load transfer functions (called “t-z” and “q-z” curves) to determine the load-settlement relationship for a pile. These load transfer curves can be developed from theoretical solutions and empirical observations.

Load transfer functions were used to describe relationship between unit skin friction and pile settlement and the relationship between the pile tip resistance and the pile tip displacements. Some examples of this application are present in Zhang & Zhang (2012aZhang, Q.Q., & Zhang, Z.M. (2012a). A simplified nonlinear approach for single pile settlement analysis. Canadian Geotechnical Journal, 49(11), 11256-11266. http://dx.doi.org/10.1139/t11-110.
http://dx.doi.org/10.1139/t11-110...
, bZhang, Q.Q., & Zhang, Z.M. (2012b). Simplified calculation approach for settlement of single pile and pile groups. Journal of Computing in Civil Engineering, 26(6), 750-758. http://dx.doi.org/10.1061/(ASCE)CP.1943-5487.0000167.
http://dx.doi.org/10.1061/(ASCE)CP.1943-...
), Wang et al. (2012)Wang, Z.J., Xie, X.Y., & Wang, J.J.C. (2012). A new nonlinear method for vertical settlement prediction of a single pile and pile groups in layered soils. Computers and Geotechnics, 45(9), 118-126. http://dx.doi.org/10.1016/j.compgeo.2012.05.011.
http://dx.doi.org/10.1016/j.compgeo.2012...
and Lee et al. (2013)Lee, J., You, K., Jeong, S., & Kim, J. (2013). Proposed point bearing load transfer function in jointed rock-socketed drilled shafts. Soil and Foundation, 53(4), 596-606. http://dx.doi.org/10.1016/j.sandf.2013.06.010.
http://dx.doi.org/10.1016/j.sandf.2013.0...
.

Nowadays most multi-storey buildings use pile groups. Thus, in order to determine the settlements, the interactive effects between them (i.e., group effect) should be considered. However, most of the current methods to estimate settlements analyze the piles separately, without considering the influence of the neighboring ones on the field behavior of the pile group.

Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
used the load transfer curves to predict settlement in isolated piles. This study analyzed several load transfer curves, comparing the curves to the results collected in the field. Fellenius (2018)Fellenius, B.H. (2018). Discussion of “development of axial pile load transfer curves based on instrumented load tests”. Journal of Geotechnical and Geoenvironmental Engineering, 144(4), 07018005. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001867.
http://dx.doi.org/10.1061/(ASCE)GT.1943-...
discussed the analyses made by Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
and it was concluded that, among the several functions analyzed in the study, the hyperbolic load transfer curves obtained good agreement to the measured results.

The present study aims to propose changes to the load transfer functions for single piles proposed by Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
in order to incorporate the group effect in the analysis of load-settlement ratios of pile groups.

2. Development of a t-z load transfer function to predict settlement in pile groups

The hyperbolic t-z load transfer function proposed by Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
to calculate the settlement of a single pile can be expressed as shown in Equation 1.

τ s = τ s , u l t s s i M s D + s s i (1)

where τs is the mobilized skin friction; τs,ult is the ultimate unit resistance; ssi is the displacement of the analyzed element (the single pile); D is the pile diameter; and Ms is the deformation parameter, which can be obtained from Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
.

In order to consider the group effect in the function, the term ssi can be isolated, as shown in Equation 2.

s s i = τ s M s D τ s , u l t τ s (2)

Following the procedure suggested by Zhang et al. (2010Zhang, Q.Q., Zhang, Z.M., & He, J.Y. (2010). A simplified approach for settlement analysis of single pile and pile groups considering interaction between identical piles in multilayered soils. Computers and Geotechnics, 37(7-8), 969-976. http://dx.doi.org/10.1016/j.compgeo.2010.08.003.
http://dx.doi.org/10.1016/j.compgeo.2010...
, 2016Zhang, Q.Q., Liu, S.W., Zhang, S.M., Zhang, J., & Wang, K. (2016). Simplified non-linear approaches for response of a single pile and pile groups considering progressive deformation of pile-soil system. Soil and Foundation, 56(3), 473-484. http://dx.doi.org/10.1016/j.sandf.2016.04.013.
http://dx.doi.org/10.1016/j.sandf.2016.0...
), Zhang & Zhang (2012a)Zhang, Q.Q., & Zhang, Z.M. (2012a). A simplified nonlinear approach for single pile settlement analysis. Canadian Geotechnical Journal, 49(11), 11256-11266. http://dx.doi.org/10.1139/t11-110.
http://dx.doi.org/10.1139/t11-110...
and Pan et al. (2018)Pan, D.D., Zhang, Q.Q., Liu, S.W., & Zhang, M.S. (2018). Analysis on response of a single pile and pile groups based on the Runge-Kutta method. KSCE Journal of Civil Engineering, 22(1), 92-100. http://dx.doi.org/10.1007/s12205-017-0578-x.
http://dx.doi.org/10.1007/s12205-017-057...
, the term related to the displacement, which is induced by the stress qs from the neighboring piles, is added to Equation 2. In sequence, the additional settlement is calculated according to the formulations proposed by Randolph & Wroth (1979)Randolph, M.F., & Wroth, C.P. (1979). An analysis of the vertical deformation of pile groups. Geotechnique, 29(4), 423-439. http://dx.doi.org/10.1680/geot.1979.29.4.423.
http://dx.doi.org/10.1680/geot.1979.29.4...
and Mylonakis & Gazetas (1998)Mylonakis, G., & Gazetas, G. (1998). Settlement and additional internal forces of grouped piles in layered soil. Geotechnique, 48(1), 55-72. http://dx.doi.org/10.1680/geot.1998.48.1.55.
http://dx.doi.org/10.1680/geot.1998.48.1...
, as shown in Equations 3 and 4.

Δ S s = ζ s τ s (3)

where ΔSs is the increase in displacement caused by the group effect; ζs is the interaction factor; and τs is the mobilized lateral resistance.

ζ s = j 1, j = 1 n r 0 G s ln 2.5 l 1 υ r i j j i , j = 1 n r 0 2 G s r i j ln 2.5 l 1 υ r i j (4)

where r0 is the radius of analyzed pile; GS is the shear modulus of the soil layer where the analyzed element is located; l is the pile length; υ is the Poisson’s ratio of the soil layer where the analyzed element is located; and rij is the distance between axes (centers) of two considered piles.

The total displacement of the shaft that is caused by the applied loads on the analyzed pile plus the neighboring ones can thus be written as in Equation 5, where the additional portion of settlement caused by the group effect is added.

s s = τ s M s D τ s , u l t q s + ζ s τ s (5)

Equation 5 can be rewritten as shown in Equation 6.

τ s = M s D + τ s , u l t ζ s + s s M s D + τ s , u l t ζ s + s s 2 4 ζ s s s τ s , u l t 2 ζ s (6)

And Equation 6 can be simplified by defining the terms ags and bgs as in Equations 7 and 8.

a g s = M s D + τ s , u l t ζ s (7)
b g s = 4 ζ s τ s , u l t (8)

By substituting the values of ags and bgs, Equation 6 can be rewritten as Equation 9, thus obtaining the expression for the modified t-z load transfer function.

τ s = a g s + s s a g s + s s 2 + b g s s s 2 ζ s (9)

The development described above was implemented by Zhang et al. (2010)Zhang, Q.Q., Zhang, Z.M., & He, J.Y. (2010). A simplified approach for settlement analysis of single pile and pile groups considering interaction between identical piles in multilayered soils. Computers and Geotechnics, 37(7-8), 969-976. http://dx.doi.org/10.1016/j.compgeo.2010.08.003.
http://dx.doi.org/10.1016/j.compgeo.2010...
for load transfer curves of single piles. The present study implemented a similar procedure but modifying the load transfer functions proposed by Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
, which were developed based upon a greater number of load tests when compared to Zhang et al. (2010)Zhang, Q.Q., Zhang, Z.M., & He, J.Y. (2010). A simplified approach for settlement analysis of single pile and pile groups considering interaction between identical piles in multilayered soils. Computers and Geotechnics, 37(7-8), 969-976. http://dx.doi.org/10.1016/j.compgeo.2010.08.003.
http://dx.doi.org/10.1016/j.compgeo.2010...
and also considered load tests carried out in Brazilian soils.

3. Development of a q-z load transfer function to predict settlement in pile groups

The hyperbolic q-z load transfer function developed by Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
to calculate settlement of single piles can be written as Equation 10.

q b = q b , u l t s b i M b D + s b i (10)

The term sbi can be isolated, in a similar procedure to that used for the t-z curve, as shown in Equation 11.

s b i = q b M b D q b , u l t q b (11)

The additional settlement was calculated according to the formulations proposed by Randolph & Wroth (1979)Randolph, M.F., & Wroth, C.P. (1979). An analysis of the vertical deformation of pile groups. Geotechnique, 29(4), 423-439. http://dx.doi.org/10.1680/geot.1979.29.4.423.
http://dx.doi.org/10.1680/geot.1979.29.4...
and can be seen in Equations 12 and 13. The variable qb is the mobilized resistance, qb,ult is the ultimate resistance, sbi is the pile toe displacement, and Mb is the toe deformation parameter. Values for Mb can be obtained from Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
.

Δ S b = ζ b q b (12)
ζ b = j = 1, j i n 1 ν b 2 π G b r i j (13)

The pile tip displacement due to the loading on top of both analyzed pile and neighboring piles can then be expressed as shown in Equation 14.

s b = q b M b D q b , u l t q b + q b ζ b (14)

Isolating the term qb, the obtained solution (Equation 15) corresponds to the modified q-z load transfer function, with the terms agb and bgb being defined respectively as in Equations 16 and 17.

q b = a g b + s b a g b + s b 2 + b g b s b 2 ζ b (15)
a g b = M b D + q b , u l t ζ b (16)
b g b = 4 ζ b q b , u l t (17)

Thus, from the expressions proposed in Equations 9 and 15, the group effect can be considered in the prediction of settlements in pile groups.

4. Analyzed case - Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...

The experimental data used as input in this study were extracted from Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
and comprised load tests performed on 32 bored piles (single and groups of 2, 4, and 9 piles), in an experimental field located in Jiangsu, China. Table 1 presents a summary with the characteristics of single piles and pile groups.

Table 1
Characteristics of piles studied by Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
.

The compressive strength of the concrete used in pile shafts and caps was 25 MPa, with a Young’s Modulus of 29.2 GPa, obtained from compression tests performed in 6 prismatic specimens (100 mm × 100 mm × 300 mm), tested 28 days after casting. According to Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
, the pile caps could be considered as rigid.

A hole called BH was drilled down to a depth of 29.50 m, aiming to obtain soil samples and the surrounding stratigraphy. In addition, 4 cone penetration tests were carried out (CPT1, CPT2, CPT3, and CPT4). A layout of the experimental field displaying the location of piles, drilled hole BH (stratigraphy), and performed CPTs can be found in Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
.

The stratigraphy showed layers of clay and silt in the upper depths. A layer of soft clay was also identified between the depths of 17.0 m and 29.50 m and water level was located at 2.6 m below ground surface.

The load tests were performed until pile failure. The reaction system consisted of prismatic blocks of precast concrete and a platform supported by 12 m-long reaction piles, installed 5 m away from the center of tested groups. The maximum load applied was 1.2 times the estimated load capacity for each group.

4.1 Method

In order to obtain the load-settlement ratios for the pile groups analyzed by Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
, the modified t-z and q-z curves (Equations 9 and 15) were used. The pile shafts were subdivided into 1 m-long segments, each one of them corresponding to a t-z curve, determined according to the properties of the soil layer in which the segment was located. This same discretization was adopted by Zhang et al. (2016)Zhang, Q.Q., Liu, S.W., Zhang, S.M., Zhang, J., & Wang, K. (2016). Simplified non-linear approaches for response of a single pile and pile groups considering progressive deformation of pile-soil system. Soil and Foundation, 56(3), 473-484. http://dx.doi.org/10.1016/j.sandf.2016.04.013.
http://dx.doi.org/10.1016/j.sandf.2016.0...
and Pan et al. (2018)Pan, D.D., Zhang, Q.Q., Liu, S.W., & Zhang, M.S. (2018). Analysis on response of a single pile and pile groups based on the Runge-Kutta method. KSCE Journal of Civil Engineering, 22(1), 92-100. http://dx.doi.org/10.1007/s12205-017-0578-x.
http://dx.doi.org/10.1007/s12205-017-057...
.

Values for Ms were obtained based on the chart presented by Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
and Poisson's ratio was assumed as 0.45, according to Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
. Table 2 presents the values for the ultimate unitary resistance (qs,ult) and Ms per section of the 20-m-long piles (i.e. groups QZ2, QZ4, and QZ9).

Table 2
Values for qs,ult and Ms for the 20 m-long pile groups.

In this method, the values for ultimate unitary lateral resistances for pile groups – the input parameters – were those regarding single piles, as presented in Table 2. Thus, the values for ultimate unitary lateral resistances obtained for pile DZ1 were used for the 20 m-long pile groups, and the values for pile DZ1L, for the 24 m-long ones.

The t-z and q-z curves for single piles were converted into pile group curves by considering the interactive effects between piles (group effect), according to the development described in Sections 2 and 3 of this study. Table 3 shows the values for ultimate unitary resistance and Ms per section for groups QZ2L, QZ4L, and QZ9L (24 m-long piles).

Table 3
Values for qs,ult and Ms for the 24 m-long pile groups.

For both 20 m- and 24 m-long piles, the adopted Mb was 0.0088. The load capacity obtained by Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
for the single pile DZ1 (20 m) was 1,430 kN and for the single pile DZ1L (24 m), 1,540 kN. Both values were obtained through load tests.

Most of the applied load was absorbed by the shaft. The toe resistance mobilized in failure was 4.73 kN for DZ1 and 78.54 kN for DZ1L. These were the values used in q-z load transfer functions.

In the approach here described, the load-settlement relationships of the piles were obtained through an interactive procedure, similar to that proposed by Coyle & Reese (1966)Coyle, H.M., & Reese, L.C. (1966). Load transfer for axially loaded piles in clay. Journal of the Soil Mechanics and Foundations Division, 92(2), 1-26. http://dx.doi.org/10.1061/JSFEAQ.0000850.
http://dx.doi.org/10.1061/JSFEAQ.0000850...
, Zhang & Zhang (2012a)Zhang, Q.Q., & Zhang, Z.M. (2012a). A simplified nonlinear approach for single pile settlement analysis. Canadian Geotechnical Journal, 49(11), 11256-11266. http://dx.doi.org/10.1139/t11-110.
http://dx.doi.org/10.1139/t11-110...
and Zhang et al. (2016)Zhang, Q.Q., Liu, S.W., Zhang, S.M., Zhang, J., & Wang, K. (2016). Simplified non-linear approaches for response of a single pile and pile groups considering progressive deformation of pile-soil system. Soil and Foundation, 56(3), 473-484. http://dx.doi.org/10.1016/j.sandf.2016.04.013.
http://dx.doi.org/10.1016/j.sandf.2016.0...
.

  • [Step 1] The pile is subdivided into n segments. A displacement Sb is chosen for the pile toe (segment n). From this displacement and using the q-z curve (Equation 15), the mobilized load Pb at pile tip is calculated, multiplying toe stress qb by the tip area.

  • [Step 2] The displacement of the midpoint of analyzed segment (Ss) is estimated. For the first iteration, Ss = Sb.

  • [Step 3] Using the t-z curve (Equation 9) and the displacement Ss, the shear stress qs mobilized in the lateral area of analyzed segment can be calculated.

  • [Step 4] The force due to friction (ΔPn) acting on the side of analyzed segment is calculated (Equation 18):

    ΔPn=τsπLi(18)

where Li is the length of analyzed segment; and τs is shear stress mobilized in the shaft of analyzed segment.

  • [Step 5] The force acting at the top of analyzed segment (Ptn) is determined according to Equation 19.

    Ptn=Pb+ΔPn(19)

  • [Step 6] The elastic shortening of the lower half of analyzed segment is calculated according to Equation 21, using the axial force which acts at the midpoint of analyzed segment (Pn,med) and the force acting at its base (Pb).

    Pn,med=Ptn+Pb2(20)
    Δe/2=Pnmed+Pb2ΔLn/2EpAp(21)

where ΔLn is the variation in the length of analyzed element; Ep is Young’s Modulus of the pile; and Ap is pile cross-sectional area.

  • [Step 7] A new displacement for the midpoint of analyzed segment (Ss') is calculated, adding the elastic shortening (Δe/2) to the displacement of the base of segment Sb, according to Equation 22.

    Ss'=Sb+Δe/2(22)

  • [Step 8] The value of Ss' is compared with Ss. If they differ by more than 10-6 m, the procedure from steps 2 to 7 is repeated, assuming Ss' as the new value for Ss, until the convergence for analyzed segment is reached.

  • [Step 9] After the convergence is reached, the displacement at the top of analyzed segment (Stn) is calculated, according to Equation 23.

    Stn=Sb+Δe(23)

  • [Step 10] The displacement of the top of analyzed segment corresponds to the displacement of the base of the upper segment. Therefore, the procedure is repeated, until load and settlement are obtained for the top segment of pile.

The procedure described above is replicated for several values of Sb until the load-settlement ratio for the range of loads of interest is obtained.

Equations 9 and 15 require the soil shear modulus of each soil layer to be applied, which in this case was estimated for each stress level. For this purpose, the soil shear wave velocity (Vs) was determined using Equation 24 for granular soil layers, proposed by Baldi et al. (1989)Baldi, G., Bellotti, R., Gionna, V.N., Jamiolkowski, M., & Lo Presti, D.C.F. (1989). Modulus of sands from CPT’s and DMT’s. In Proceedings of the 12th International Conference on Soil Mechanics and Foundation Engineering (Vol. 1, pp. 1-6), Rio de Janeiro. Publications Committee of XII ICSMFE. and Equation 25 for fine soil layers, proposed by Mayne (1995)Mayne, P.W. (1995). Evaluation of foundation response by G/Gmax degradation relationships in soils (pp. 136-148). Taipei: National Science Council of ROC and NSF of U.S., which correlate CPT’s penetration resistance with Vs.

V s = 277 q c 0.13 σ v o ' 0.27 (24)
V s = 1.75 q c 0.627 (25)

where σ'vo is the effective geostatic stress of soil in the center of analyzed layer.

From the shear wave values for each layer, the maximum shear modulus (Gmax) can be calculated, according to Equation 26.

G m a x = γ n g V s 2 (26)

where g is the acceleration of gravity; and γn is the soil unit weight.

From Gmax, the shear modulus (G) is obtained for each analyzed load, using another iterative procedure, as described below (Steps A to E). This geotechnical parameter is obtained using Equations 24 and 25, for each layer, considering the thickness of each layer.

  • [Step A] A load at the top of pile (Pt') is estimated. From the load capacity of the single pile, a factor of safety (FS) is determined.

  • [Step B] Using Equation 27 (Fahey & Carter, 1993Fahey, M., & Carter, J.P. (1993). A finite element study of the pressuremeter test in sand using nonlinear elastic plastic model. Canadian Geotechnical Journal, 30(2), 348-362. http://dx.doi.org/10.1139/t93-029.
    http://dx.doi.org/10.1139/t93-029...
    ), the shear modulus (Gs) of each layer can be also determined.

    GsGmáx=11FS0.3(27)

  • [Step C] The iterative procedure described in Steps A and B is performed again, using the calculated values for G. Then, a new load Pt is obtained.

  • [Step D] The obtained Pt is compared with the estimated Pt. If the obtained value is smaller than the estimated one, a larger displacement Sb is used in the beginning of the procedure. If the obtained Pt is greater, a smaller Sb displacement should be used.

  • [Step E] The procedure is repeated until the difference between the estimated Pt' and the obtained Pt is not greater than 0.1 kN.

5. Results of predictions and analyses

5.1 Analyzed case study – Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...

The case study described by Dai et al. (2012)Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
was used in this study to evaluate the modified t-z and q-z curves and the proposed methodology. At the end of this section, some possible causes for differences between the predicted (calculated) and measured results will be presented.

Figure 1 shows both the predicted and the experimental (obtained by Dai et al., 2012Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
) load-settlement curves for groups QZ2 (2 piles) and QZ2L (2 piles). The method proposed in this study yielded results that agreed quite well with the experimental ones, especially for loads up to 2,500 kN. For higher loads, the predicted settlements were slightly higher than those experimentally measured.

Figure 1
Predicted and experimental load-settlement curves for (a) group QZ2 and (b) group QZ2L.

Considering a factor of safety (FS) of 2.0 (i.e. an adopted workload of 1,250 kN), a difference of only 1.5 mm (absolute value) was observed between the predicted and the experimentally measured settlements.

As to group QZ2L, the predicted settlements also agreed well with the experimental results for loads up to 1,500 kN. The predicted results were closest to those measured for the load of 1,100 kN. Between 1,500 kN and 2,800 kN, the predicted (calculated) results were slightly lower than those obtained experimentally. Above 2,800 kN, the method yielded higher values for settlement than those obtained experimentally.

Figure 2 shows the predicted and the experimental (Dai et al., 2012Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087.
http://dx.doi.org/10.1139/t2012-087...
) load-settlement curves for groups QZ4 (4 piles, 20-m-long piles) and QZ4L (4 piles, 24-m-long piles). Good agreement was observed between predicted and experimental results up to 1,500 kN. For higher loads, the predicted settlements were smaller than those obtained experimentally. These were the most discordant results (between experimental and predicted) considering all analyzed groups, which might have happened due to differences in stiffness of the soils around the single piles and the pile groups.

Figure 2
Predicted and experimental load-settlement curves for (a) group QZ4 and (b) group QZ4L.

Regarding group QZ4L, predicted results were remarkably close to experimental ones, notably for loads below 1,500 kN. Between 1,500 kN and 2,800 kN, experimental results were slightly higher than estimated ones. For load levels above 2,800 kN, the method provided settlements also slightly higher than those experimentally obtained.

Figure 3 shows the load-settlement curves for pile groups QZ9 and QZ9L. For group QZ9, excellent agreement between predicted and experimental results was obtained, especially for loads under 7,500 kN. For higher loads, the predicted results were slightly lower than those experimentally obtained. One of the possible causes for the divergence may be a reduced soil stiffness in the area of ​​installation of this group, in comparison to the area of ​​installation of isolated piles.

Figure 3
Predicted and experimental load-settlement curves for (a) group QZ9 and (b) group QZ9L.

For group QZ9L, predicted and experimental results agreed as well. The most consistent results were for loads between 0 and 3,000 kN. Around 6,000 kN, the predicted settlements were slightly higher than the experimental ones. For loads above 8,000 kN, the experimental results were slightly higher than those predicted.

In general, the method proposed in this study yielded settlements that reasonably agreed with experimental results, notably for groups QZ4L and QZ9. However, some divergences between experimental and predicted results were still observed, particularly for group QZ4.

These small divergences probably happened due to soil heterogeneity, for CPT results showed that soil stiffness was not the same in all assessed locations. Thus, as the used input parameters came both from load tests and CPTs carried out on single piles, it is possible that the stiffness in the regions where the load tests were performed was different from where CPTs were performed.

Another possible reason for the differences in the results reside in the choice for the parameters Ms, Mb, and shear modulus. Although estimated considering the results of 72 load tests, more accurate values for Ms and Mb could have been obtained if a greater number of load tests had been used for each pile and soil types.

As to shear modulus, although the correlations were calibrated considering a large amount of data and showed good adherence to them, it is possible that the values could be a little distant from field reality, which could also be a possible cause for the divergences between predicted and experimental settlements.

6. Conclusions

The present work proposed modifications to the t-z and q-z curves developed by Bohn et al. (2016)Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
https://doi.org/10.1061/(ASCE)GT.1943-56...
, in order to consider the group effect in the analysis of load-settlement ratios in pile groups. The modifications were carried out using the interaction factors proposed by Randolph & Wroth (1979)Randolph, M.F., & Wroth, C.P. (1979). An analysis of the vertical deformation of pile groups. Geotechnique, 29(4), 423-439. http://dx.doi.org/10.1680/geot.1979.29.4.423.
http://dx.doi.org/10.1680/geot.1979.29.4...
and Mylonakis & Gazetas (1998)Mylonakis, G., & Gazetas, G. (1998). Settlement and additional internal forces of grouped piles in layered soil. Geotechnique, 48(1), 55-72. http://dx.doi.org/10.1680/geot.1998.48.1.55.
http://dx.doi.org/10.1680/geot.1998.48.1...
.

Regarding groups QZ2 and QZ2L, the method proposed in this study yielded results with good agreement with experimental ones, especially for loads up to 2,500 kN. For loads above 2,500 kN, the predicted settlements were slightly higher than those experimentally measured.

For the group QZ4, good agreement was also observed between predicted and experimental results, mainly for loads up to 1,500 kN.

For the group QZ4L, predicted results were quite close to experimental ones, notably for loads below 750 kN. For load levels above 1,400 kN, the predicted settlements were slightly higher than those experimentally obtained.

Regarding groups QZ9 and QZ9L, good agreement between predicted and experimental results was also obtained, especially for loads under 7,500 kN. For higher loads, the predicted results were slightly lower than experimental ones. The geotechnical parameter Gmax were obtained by correlations with CPT.

Thus, the predictions for load-settlement ratios made according to the modifications proposed in this study for t-z and q-z curves agreed quite well with the experimental data for most analyzed pile groups, even though considerable divergences between experimental and predicted settlements were also obtained (block QZ4), most likely due to soil heterogeneity, evidenced by results from cone penetration tests (CPTs).

List of symbols

τsMobilized unit lateral resistance

τs,ultUltimate unit lateral resistance

DPile diameter

MsDeformation parameter for t-z curve

ssiDisplacement of the analyzed element

ζsLateral interaction factor

ΔSsIncrease in displacement of the analyzed element caused by the group effect

r0Radius of pile

GSShear modulus of soil layer where the element is located

υPoisson’s ratio of soil layer where the element is located

rijDistance between axes (centers) of two considered piles

qbMobilized unit base resistance

qb,ultUltimate unit base resistance

sbiBase displacement

MbDeformation parameter for q-z curve

ΔSsIncrease in base displacement caused by the group effect

ζbBase interaction factor

GbShear modulus of the base soil layer

υsPoisson’s ratio of the base soil layer

ΔPnForce due to friction

LiLength of analyzed segment

PtnForce acting at the top of analyzed element

PbForce acting at the base of analyzed element

Pn,medAxial force which acts at the midpoint of analyzed segment

ΔLn Variation in the length of analyzed element

ΔeElastic shortening

Ss'Displacement of the midpoint of analyzed segment

StnDisplacement at the top of analyzed segment

VsSoil shear wave velocity

qcCone penetration resistance

GmaxMaximum soil shear modulus

gAcceleration of gravity

γnSoil unit weight

FSSafety factor

  • Discussion open until August 31, 2021.

Acknowledgements

To the Post-Graduate Program POSDEHA. To the CNPQ - National Council for Scientific and Technological Development and to the Coordination for funding the research. To the CAPES Coordination for the Improvement of Higher Education Personnel.

References

  • Baldi, G., Bellotti, R., Gionna, V.N., Jamiolkowski, M., & Lo Presti, D.C.F. (1989). Modulus of sands from CPT’s and DMT’s. In Proceedings of the 12th International Conference on Soil Mechanics and Foundation Engineering (Vol. 1, pp. 1-6), Rio de Janeiro. Publications Committee of XII ICSMFE.
  • Bohn, C., Santos, A.L., & Frank, R. (2016). Development of axial pile load transfer curves based on instrumented load tests. Journal of Geotechnical and Geoenvironmental Engineering, 143, 2-6. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001579.
    » https://doi.org/10.1061/(ASCE)GT.1943-5606.0001579
  • Coyle, H.M., & Reese, L.C. (1966). Load transfer for axially loaded piles in clay. Journal of the Soil Mechanics and Foundations Division, 92(2), 1-26. http://dx.doi.org/10.1061/JSFEAQ.0000850
    » http://dx.doi.org/10.1061/JSFEAQ.0000850
  • Dai, G., Salgado, R., Gong, W., & Zhang, Y. (2012). Load tests on full-scale bored pile groups. Canadian Geotechnical Journal, 49(11), 1293-1308. http://dx.doi.org/10.1139/t2012-087
    » http://dx.doi.org/10.1139/t2012-087
  • Fahey, M., & Carter, J.P. (1993). A finite element study of the pressuremeter test in sand using nonlinear elastic plastic model. Canadian Geotechnical Journal, 30(2), 348-362. http://dx.doi.org/10.1139/t93-029
    » http://dx.doi.org/10.1139/t93-029
  • Fellenius, B.H. (2018). Discussion of “development of axial pile load transfer curves based on instrumented load tests”. Journal of Geotechnical and Geoenvironmental Engineering, 144(4), 07018005. http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001867
    » http://dx.doi.org/10.1061/(ASCE)GT.1943-5606.0001867
  • Lee, J., You, K., Jeong, S., & Kim, J. (2013). Proposed point bearing load transfer function in jointed rock-socketed drilled shafts. Soil and Foundation, 53(4), 596-606. http://dx.doi.org/10.1016/j.sandf.2013.06.010
    » http://dx.doi.org/10.1016/j.sandf.2013.06.010
  • Mayne, P.W. (1995). Evaluation of foundation response by G/Gmax degradation relationships in soils (pp. 136-148). Taipei: National Science Council of ROC and NSF of U.S.
  • Mylonakis, G., & Gazetas, G. (1998). Settlement and additional internal forces of grouped piles in layered soil. Geotechnique, 48(1), 55-72. http://dx.doi.org/10.1680/geot.1998.48.1.55
    » http://dx.doi.org/10.1680/geot.1998.48.1.55
  • Pan, D.D., Zhang, Q.Q., Liu, S.W., & Zhang, M.S. (2018). Analysis on response of a single pile and pile groups based on the Runge-Kutta method. KSCE Journal of Civil Engineering, 22(1), 92-100. http://dx.doi.org/10.1007/s12205-017-0578-x
    » http://dx.doi.org/10.1007/s12205-017-0578-x
  • Randolph, M.F., & Wroth, C.P. (1979). An analysis of the vertical deformation of pile groups. Geotechnique, 29(4), 423-439. http://dx.doi.org/10.1680/geot.1979.29.4.423
    » http://dx.doi.org/10.1680/geot.1979.29.4.423
  • Wang, Z.J., Xie, X.Y., & Wang, J.J.C. (2012). A new nonlinear method for vertical settlement prediction of a single pile and pile groups in layered soils. Computers and Geotechnics, 45(9), 118-126. http://dx.doi.org/10.1016/j.compgeo.2012.05.011
    » http://dx.doi.org/10.1016/j.compgeo.2012.05.011
  • Zhang, Q.Q., Zhang, Z.M., & He, J.Y. (2010). A simplified approach for settlement analysis of single pile and pile groups considering interaction between identical piles in multilayered soils. Computers and Geotechnics, 37(7-8), 969-976. http://dx.doi.org/10.1016/j.compgeo.2010.08.003
    » http://dx.doi.org/10.1016/j.compgeo.2010.08.003
  • Zhang, Q.Q., & Zhang, Z.M. (2012a). A simplified nonlinear approach for single pile settlement analysis. Canadian Geotechnical Journal, 49(11), 11256-11266. http://dx.doi.org/10.1139/t11-110
    » http://dx.doi.org/10.1139/t11-110
  • Zhang, Q.Q., & Zhang, Z.M. (2012b). Simplified calculation approach for settlement of single pile and pile groups. Journal of Computing in Civil Engineering, 26(6), 750-758. http://dx.doi.org/10.1061/(ASCE)CP.1943-5487.0000167
    » http://dx.doi.org/10.1061/(ASCE)CP.1943-5487.0000167
  • Zhang, Q.Q., Liu, S.W., Zhang, S.M., Zhang, J., & Wang, K. (2016). Simplified non-linear approaches for response of a single pile and pile groups considering progressive deformation of pile-soil system. Soil and Foundation, 56(3), 473-484. http://dx.doi.org/10.1016/j.sandf.2016.04.013
    » http://dx.doi.org/10.1016/j.sandf.2016.04.013

Publication Dates

  • Publication in this collection
    23 June 2021
  • Date of issue
    2021

History

  • Received
    13 Dec 2020
  • Accepted
    05 May 2021
Associação Brasileira de Mecânica dos Solos Av. Queiroz Filho, 1700 - Torre A, Sala 106, Cep: 05319-000, São Paulo - SP - Brasil, Tel: (11) 3833-0023 - São Paulo - SP - Brazil
E-mail: secretariat@soilsandrocks.com