11th ICSGE

17-19 May 2005

Ain Shams University

Cairo - Egypt

Eleventh International Colloquium on Structural and Geotechnical Engineering

ANALYSIS OF STABILIZING SLOPES USING VERTICAL PILES

GEHAN E. ABDELRAHMAN1 MAHMOUD S. ABDELBAKI2 YOUSSEF G. YOUSSEF3 ABSTRACT This parametric study reveals the important parameters that affect the behavior and factor of safety of piles in slopes. These parameters are; the pile position, pile diameter and depth, soil properties, soil layer thickness and surcharge load. The effect of different undrained shear strength of clayey soil on the pile bending moment and the slope factor of safety are also studied. The nonlinear elasto-plastic finite element program PLAXIS is used to analyze the response of vertical piles due to the lateral soil movement generated by road surcharge. In addition an X-STABLE program is used to calculate the critical slip surface for the different cases studied before pile installation to predict the required pile length for stabilization. A comparison is made for no stabilization, stabilization using one row of piles and stabilization using two rows of piles. The studied case is for soft clayey soil overlying sandy soil. The results show that piles designed to stabilize slopes, raise the factor of safety and the allowable road surcharge in case of using one or two pile rows, for that type of soil. A comparison is made between the results of the finite element analysis and case study results from Alsalam-Canal. In many cases the results are reasonably in good agreement. Keywords: safety factor, laterally loaded piles, PLAXIS, road surcharge, slope stability. 1 INTRODUCTION Piles may be designed to restrain soil movements when used to stabilize slopes or potential landslides. The lateral loads resulting from the soil movement induce bending moments and deflections in the pile which may lead to their structural failure. Many theoretical and empirical methods were introduced for solving certain types of problems. Stewart et al. [1], Chen and Poulos [2], Chow [3], Hassiotis et al. [4] introduced other methods for the design of slopes stabilized with a single row of piles. Chen and Poulos [2], took into account the pile group effect and the presence of a slope adjacent to two pile rows. The methods of the analysis of piles and pile groups subjected to lateral loading from horizontal soil movements, have been classified into four categories as described by Stewart et al. [5]. One of those method is the comprehensive finite-element analysis, where the piles are represented in the mesh and the overall soil-pile-embankment response is included. This method is adopted here using the nonlinear elasto-plastic finite element program PLAXIS.
1

Lecturer civil Engineering Departement, Cairo University-Fayoum Branch. Fayoum, Egypt [email protected] 2 Lecturer civil Engineering Departement, Cairo University-Fayoum Branch. Fayoum, Egypt 3 Assistant Lecturer civil Engineering Departement, Cairo University-Fayoum Branch. Fayoum, Egypt E05GE28 - 1

2 GEOTECHNICAL ANALYSIS Two slope profiles have been studied. The first consists of homogenous soft to medium silty clay, whereas the second is made of homogenous soft to medium silty clay followed by a sand layer as shown in Figure 1.

Fig. 1: Schematic drawing of the soil profile and pile rows (PLAXIS). The soft clay was assigned a rang of different values of undrained cohesion as shown in Table 1. These values are in the range of the soft to medium clay as stated in the Egyptian code of practice. The value of the Modulus of Elasticity is assumed to be 500 Cu. The piles are modelled as a beam element using the moment of inertia of a single pile and distributing it on every meter to get the moment of inertia/m as shown in Table 2. Since the analysis presented here is intended for addressing the slope stability problem of Alsalam-Canal. The water level is modelled to be at the water level in the Alsalam-Canal with a seepage line 15:1 towards the inland. A paved road exists on the right bank of the canal. The width of the canal is about 32m, while the width of the road is 12m. The only external load applied is the live load on the road. Two different approaches are adopted to study the canal bank stability. In the first, the vertical load intensity on the bank is increased up to failure. This method was used for three cases; without piles, with one row of piles, and with two rows of piles. In the second approach, the value of the cohesion is decreased till failure and the respective load at failure is determined.

E05GE28 - 2

Table 1: Soil properties Parameter Dry unit weight ( kN/m3 ) Wet unit weight ( kN/m3 ) Young's modulus E ( kN/m ) Poisson's ratio Cohesion C ( kN/m ) Friction angle Type of behavior
2 2

Clay 16 18 500 C 0.49 10,20,30,40 0 Undrained

Sand 17 20 3x104 0.3 0 31 Drained

Table 2: Pile properties Parameter Young's Modulus, E ( kN/m2 ) EA (kN/m) EI (kN.m2/m) Poisson's ratio, Pile Diameter, D (m) Pile Length, L (m) Pile 2.1x107 4.75x106 1.07x105 0.3 0.5 17

3 ANALYTICAL MODEL The slope stability analysis is performed using the specialized finite element program, PLAXIS. The program is an elasto-plastic, plane strain, finite element program for soil modeling. The soil is modeled using 6 nodded triangular elements. The soil model is MohrColumb method with nonlinear failure envelope. Both vertical and lateral soil movements may occur, but only the lateral soil movement is of significance in this case. The pile is modeled as a simple elastic beam with interface elements to model the soil structure interaction. The strength ratio of the interface elements, R, is assumed to be 0.7. The lateral displacement of each element of the pile can be related to the pile bending stiffness and the horizontal pile-soil interaction stresses. The lateral displacement of the corresponding soil elements is related to the soil modulus or stiffness, the pile-soil interaction stresses, and the free-field horizontal soil movements. A limiting lateral pile-soil stress can be specified so that local failure of the soil can be allowed for, thus allowing a nonlinear response to be obtained. The stability of a slope can be investigated by a number of limit equilibrium methods. The program X-STABLE is used to calculate the actual slip surface plane for the different shear strength cases, using Janbu slices method before pile installation. 4 FACTORS INFLUENCING THE LATERAL PILE RESPONSE The parametric study presented by Poulos [6] identified some of the most important factors influencing the behavior of piles within or near an embankment on clay. The most significant factors appear to include pile position, undrained shear strength of clay, soil layer thickness,

E05GE28 - 3

and road surcharge intensity. He also indicated that delayed installation of the pile had substantial benefit in reducing pile bending moments and lateral movements. Table (3) presents the results of the parametric study performed in the current study. Table 3: The effect of clay cohesion on pile bending moment at failure.
Cohesion Clay One Pile C kN/m 10 12 14 15 16 17 18 25 30 40
2

Clay underlain by Sand One Pile M1 kN.m 1300.0 1210.0 1040.0 947.1 904.0 845.1 792.7 544.6 417.7 308.0 M2 kN.M 23.0 24.6 48.2 52.2 59.1 68.5 77.5 106.6 125.1 159.0

Clay Two Piles Canal Pile Road Pile M M kN.m 719.8 537.2 453.0 443.0 456.7 443.0 376.0 374.7 416.8 419.6 kN.m 1010.0 814.6 820.0 828.6 733.7 719.4 700.0 647.3 694.1 696.7

Clay underlain by Sand Two Piles Canal Pile M1 M2 kN.m 310.0 312.0 312.0 312.7 306.7 297.3 284.7 190.8 128.3 30.4 kN.m 437.0 442.4 366.9 349.9 335.1 324.3 315.5 295.0 275.0 250.0 Two Piles Road Pile M1 M2 kN.m 344.1 338.7 342.8 312.0 277.7 259.0 242.9 169.9 98.9 54.4 kN.m 408.0 425.0 340.7 327.0 319.1 310.6 304.4 282.4 284.1 273.2

M kN.m 261.2 250.7 268.6 247.4 272.5 265.9 249.2 264.6 266.7 268.2

* M1 maximum negative moment * M2 maximum positive moment 4.1 Effect of soil layer thickness Two types of soils are used, The first type consisted of extended homogeneous clay, while the other consisted of homogeneous clay underlain by a bearing stratum of sand. A relation is obtained between pile bending moment and the undrained shear strength of clayey soil C for the two soil types. As shown in Figure 2 the maximum negative bending moment, M1, decrease as the undrained shear strength increases in case the pile extend in underlying bearing stratum. For the extended homogeneous clay, the pile bending moment is not affected by the increase in the undrained shear strength as shown in Figure 2.
1400 1200
(M Pile In Clay) (M1 Pile In Sand) (M2 Pile In Sand)

M om ent M , (kN .m )

1000 800 600 400 200 0 10 12 14 15 16 17 18

25

30

40

Undrained Cohesion C, (kN/m2)

Fig. 2: One pile row floating in clayey soil and another ending in sandy soil.
E05GE28 - 4

4.2 Pile arrangement (piles in one and two rows.) In the finite-element analysis, the pile-soil-pile system is simplified to a two-dimensional plain strain study, and is loaded by a soil movement due to road surcharge, which is applied incrementally until the system reaches an ultimate plastic state. Chen et al. [2] studied the group effect, fp, on the ultimate pile load, pu, in case of floating pile. In this study the transverse pile spacing as shown in Figure 1 Sh= 3D while in the longitudinal direction the pile spacing Sv = 6D, where D is the pile diameter. According to Chen et al, fp equals 1.1 in case of one long row. In case of two long rows, for the pile located close to the dredged side, fp, equals to 1.4, while for the pile located close to road side, fp, equals to 1.2. Consistent with the above-mentioned results, in case of two pile rows, the pile near the canal side has bending moment less than the pile that lies near the road side as shown in Figure 3. In case of one pile row the pile bending moment is less than that of two pile rows. This is attributed to the greater moment of inertia of two pile rows compared to the single row of piles.
1400 1200 Moment, (kN.m) 1000 800 600 400 200 0 10 12 14 15 16 17 18 25 Undrained Cohesion C, (kN/m2) 30 40

One Pile Canal Pile Road Pile

Fig. 3: One and two pile rows floating in clay. 4.3 Pile size and position In all analysis cases, the pile length has been chosen deeper than the slip failure surface. XSTABLE program was used to calculate the position of the slip surface before pile installation. This slip surface was calculated to be at 12.65m from the berm level. Thus the selected pile length of 17m, is suitable. Poulos [7] found that piles installed within the embankment are subjected to significantly larger bending moments than piles at the toe, due to both larger horizontal movements and also the effect of the embankment fill moving past the pile. The maximum moment appears to be developed for piles located near the middle of slope, and is more than twice the value for piles at the embankment toe. In this study the pile is locate at the berm level which is considered at the mid-slope according to Poulos [7]. Piles floating in clayey soil had maximum bending moment double that for those ending in the sand layer. The existence of underlying bearing stratum of sand decreases the bending moment on the pile as shown in Figure 4.

E05GE28 - 5

1400 1200

(M Piles In Clay) (M1 Piles In Sand) (M2 Piles In Sand)

Moment M, (kN .m)

1000 800 600 400 200 0 10 12 14 15 16 17 18

Undrained Cohesion C, (kN/m )

25

30

40

Fig. 4: Bending moment of pile near the canal side.

(M Piles In Clay) (M1 Piles In Sand) (M2 Piles In Sand)

M oment M , (kN .m)

1000 800 600 400 200 0 10 12 14

Undrained Cohesion C, (kN/m )

15

16

17

18

25

30

40

Fig. 5: Bending moment of pile near the road side. In case of two pile rows and according to group effect, the pile near canal side had bending moment less than the pile which lies near the road side as shown in Figures 4 and 5 in case of piles floating in clayey soil. In case of two pile rows ending in sandy layer, both piles had the same bending moment and equal value for both maximum negative bending moment, M1, and maximum positive bending moment, M2, along the pile length as shown in Figure 5. 5 MAXIMUM BENDING MOMENT WITH UNDRAINED SHEAR STRENGTH Figure 6 shows that the increase in undrained shear strength decreases the pile bending moment. This may be due to the smaller lateral movements which are developed in the stronger soils. In case of one pile row the maximum negative bending moment, M1, is much higher than maximum positive bending moment, M2, as shown in Figures 6and 7, while in case of two piles M1 and M2 had the same value for both canal and road piles. Pile position has no effect on the value of bending moment.

E05GE28 - 6

1400 1200

One Pile Canal Pile Road Pile

Moment M1,(kN.m)

1000 800 600 400 200 0 10 12 14 15 16 17 18

25

30

40

Undrained Cohesion C, (kN/m )

Fig. 6: Maximum negative pile moment, M1, in one and two pile rows socket in sandy soil.
1400 M oment M2, (kN .m) 1200 1000 800 600 400 200 0 10 12 14 15 16 17 18
2

One Pile Canal Pile Road Pile

25

30

40

Undrained Cohesion C, (kN/m )

Fig. 7: Maximum positive pile moment, M2, in one and two pile rows socket in sandy soil.

E05GE28 - 7

6 SURCHARGE LOAD The increase in failure (surcharge) load with the shear strength of the soil is shown in Figure 8. It is seen that the percentage of improvement increases as the shear strength of the soil decreases. At higher values of cohesion (more than 40 kN/m2), the increase in the percentage of improvement decreases significantly and the need to use piles for stabilization is thus diminished. It may be seen that the single row and two pile rows ending in sand layer increased the failure load on the road to a value of about 28kN/m2. One or two pile rows had the same effect either floated in clay or ended in sand on the surcharge causing failure as shown in Figure 8.
175 NO PILE (CLAY) 150 125 100 75 50 25 0 10 12 14 15 ONE PILE (CLAY) TWO PILES (CLAY) NO PILE (SAND) ONE PILE (SAND) TOW PILES (SAND)

Failure Surcharge Load Q, kN /m2

Undrained Cohesion, C, kN/m2

16

17

18

25

30

40

Fig. 8: Road surcharge with different undrained shear strength ( F.O.S. =1) 7 FACTORS INFLUENCING FACTOR OF SAFETY The geotechnical analysis is concerned with calculating the value of the safety factor for the different cases and the comparison between the performances of the system in the different cases. The safety factor may be calculated by two different methods. The first which used in this study, is the safety factor in terms of the surcharge load. This was calculated as the ratio between the value of the maximum load at failure for each case and the value of the applied load. The second method of safety factor calculation is in terms of the shear strength of the soil; thats the ratio between the actual value of shear strength and the value of shear strength at failure. Both the presence of sand layer underlying the clay strata, or the existence of one or two pile rows floating in clayey soil had no effect on the factor of safety. However with the increase of undrained shear strength, a slight increase occurs in the factor of safety, in case of one or two pile rows floating in clay. In case of piles ending in the sand layer, factor of safety increased as shown in Figure 9. The two pile rows give higher factor of safety than that for the one pile row.

E05GE28 - 8

Table 4: Calculated factor of safety in different cases

Clay No Pile 0.80 1.08 1.27 1.36 1.45 1.54 1.63 2.26 2.72 3.61 One Pile 1.06 1.13 1.35 1.42 1.54 1.63 1.70 2.40 2.88 3.84 Two Piles 1.05 1.15 1.34 1.43 1.53 1.63 1.72 2.40 2.87 3.83

Clay underlain by Sand No Pile 1.01 1.10 1.28 1.36 1.45 1.54 1.63 2.27 2.73 3.65 One Pile 1.47 1.56 1.77 1.87 1.95 2.04 2.13 2.81 3.28 4.23 Two Piles 1.65 1.76 1.94 2.04 2.14 2.23 2.33 3.02 3.53 4.54

5.0 4.5 4.0 3.5 NO PILES(CLAY) ONE PILE(CLAY) TWO PILES(CLAY) NO PILES (SAND) ONE PILE(SAND)

Factor of Safety FOS

3.0 2.5 2.0 1.5 1.0 0.5 0.0 10 12

14

15

16

17

18
2

25

30

40

Undrained Cohesion (C) kN/m

Fig. 9: Factor of safety with different undrained shear strength. 8 COMPARISON WITH FIELD MEASUREMENT The embankments of the Alsalam-Canal were made of un-compacted marine deposits, and was founded over soft clay soil in the north of Egypt. The stability of the embankment slopes was studied. In many sections of the canal, the factor of safety of the side slopes did not meet the requirements of the Egyptian code. The stabilizing system consisted of reinforced concrete piles with a diameter of 0.5 m spaced centerto-center at 1.25m parallel to the canal. the lateral movements, vertical movements, and the straining actions, were measured. The measurements prove the efficiency of the stabilizing system in reducing significantly the lateral movements of the embankment mass. The strains measured in the pile reinforcement bars show that the critical factor of safety of the pile is 1.63 during either the construction or loading process. Abdel-Motalab [8] concluded that the stabilizing system has efficiently controlled the sliding of the embankment side slopes under the most critical loading conditions. The maximum negative moment, M1, occurred near the bed level. Below the bed level, where the major resisting forces exist, the moment started to reverse to reach the maximum positive moment, M2, at the surface of the lower layer.
E05GE28 - 9

Furthermore, the bending moment diminished near the surface of the continuous sand layer where the pile is embedded for at least 2.0 m. The maximum bending moment was about 163.76 kN.m, which is about 61% of its maximum bending capacity 267 kN.m. That means that the factor of safety of the pile is 1.63 during both the construction and loading process. By comparing the finite element results in Figure 9 for the case of extended clay layer or clay underlain by sand without piles with the case of using two pile rows, the factor of safety increased from 1 to 1.6. Finite element analysis using PLAXIS program, shows good agreement with that measured in the field in the case when using two pile rows. 9 COCLUSIONS 1. The benefit of using piles in stabilizing clay slop is more apparent at lower values of cohesion, (C>40 kN/m2). 2. The presence either of one or two pile rows floating in clayey soil caused slight increases to failure surcharge on the road, and no effect on the factor of safety. 3. The presence of one pile or two piles rows in case of a firm sand layer below the soft clay, increased the failure load on the road. The safety factor also increased by 60%. Two double rows of piles is more effective in this case than one pile row. 4. Finite element analysis using PLAXIS program, shows good agreement with that measured in the field in case of using two pile rows. REFERENCES [1] Stewart, D. P. (1992), Analysis of Piles Subjected to Embankment Induced Lateral Soil Movements. Journal of Geotechnical and Geo-environmental Engineering, pp. 425-435. [2] Chen, L.T. and Poulos, H.G. (1996), Piles Subjected to Lateral Soil Movements. Journal of Geotechnical and Geo-environmental Engineering, pp. 802-811. [3] Chow, Y.K. (1996), Analysis of Piles Used for Slope Stabilization. International Journal for Numerical and Analytical Methods in Geomechanics, Vol.20, pp.635-646. [4] Hassiotis, S., Chameau, J.L., and Gunaratne, M. (1997), Design Method for Stabilization of Slopes with Piles. Journal of Geotechnical and Geo-environmental Engineering, ASCE, Vol.123, No. 4, pp. 314-323. [5] Stewart, D. P., Jewell, R. J. and Randolph, M.F. (1994), Design of Piled Bridge Abutments on Soft Clay for Loading From Lateral Soil Movements. Geotechnique, No.2, pp.277-296. [6] Poulos, H.G. (1994), Analysis and Design of Piles Through Embankments. Design and Construction of Deep Foundation, Vol. III, pp1403-1421. [7] Poulos, H.G. (1976), Behavior of Laterally Loaded Piles Near a Cut or Slope. Australian Geomechanics, July, Vol. C6, No.1, pp1-12. [8] Abdel-Motaleb, A.A. (2005), Monitoring of A Stabilizing Embankment Using Reinforced Concrete Piles. The Fifth International Geotechnical Engineering Conference-Cairo University, January, pp. 183-199.

E05GE28 - 10