Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2485-2492 2485 www.etasr.com Soomro et al.: Single Pile Settlement and Load Transfer Mechanism due to Excavation in Silty Clay Single Pile Settlement and Load Transfer Mechanism due to Excavation in Silty Clay Mukhtiar Ali Soomro Kiran Fatima Memon Mohsin Ali Soomro Department of Civil Engineering Quaid-e-Awam University of Engineering, Science & Technology Nawabshah, Pakistan eng.soomro@gmail.com Department of Civil Engineering Quaid-e-Awam University of Engineering, Science & Technology Nawabshah, Pakistan kiranf893@gmail.com Department of Civil Engineering Quaid-e-Awam University of Engineering, Science & Technology Nawabshah, Pakistan drmohsin@quest.edu.pk Asad Memon Manthar Ali Keerio Department of Civil Engineering Quaid-e-Awam University of Engineering, Science & Technology Nawabshah, Pakistan asadmemon10@hotmail.com Department of Civil Engineering Quaid-e-Awam University Engineering, Science & Technology Larkana, Pakistan mantharali99@quest.edu.pk Abstract—In densely built areas, development of underground transportation system often involves excavations for basement construction and cut-and-cover tunnels which are sometimes inevitable to be constructed adjacent to existing piled foundations. In order to gain new insights into single pile responses (i.e. settlement and load transfer mechanism) to an adjacent excavation in saturated silty clay, a three-dimensional coupled- consolidation numerical analysis is conducted in this study. An advanced hypoplasticity (clay) constitutive model with small-strain stiffness was adopted. A linear increase in pile settlement was observed due to excavation-induced stress release. This is because part of the pile is placed within the boundaries of a major influence zone due to excavation-induced ground movement. Based on a settlement criterion, apparent loss of pile‘s capacity is 14%. A maximum bending moment of about 350 kNm is induced in the pile with the maximum deflection of 28 mm. In addition, mobilisation of shear strength at the pile-soil interface was found to be a key factor governing pile-soil-excavation interaction. During excavation, a downward load-transfer mechanism in the piles can be identified. Keywords-excavation; pile; silty clay; numerical modelling I. INTRODUCTION The development of underground space induces ground movements and stress release in the surrounding soil that may affect existing structures. Owing to the inherent lack of surface space in congested urban areas, deep excavation for basements and cut-and-cover tunnels are inevitably constructed adjacent to existing pile foundations often supporting high-rise buildings [1-2]. Thus, it is important for designers to estimate potential damage to existing piles resulting from nearby deep excavations. Most of the earlier researchers estimated the buildings settlement and tilting considering wall movements and ground surface settlement trough using empirical approaches. The performance of pile group depends on the stress state in the soil and surrounding sub-surface soil movements [3]. In addition, deep excavation in soft clay induces negative excess pore water pressure [4] which induces long term pile group settlement with the dissipation of excess pore water pressure. Therefore, it is essential to investigate the pile group response mechanism adjacent to deep excavation in soft soil. Authors in [1, 2] reported case studies in granular soil and alluvium residual soil respectively. They demonstrated that lateral soil movements due to excavation can be detrimental to nearby existing piles. In both reported case studies, pile toe level was much deeper than the excavation level and they reported only the lateral pile behavior. Apart from field monitoring, a number of centrifuge tests were also conducted to investigate the response of single pile [5] and pile group in soft kaolin clay [6]. They concluded that the induced bending moment and lateral pile deflection were highly influenced by the distance from the wall and pile head condition. In studies, lateral response of end bearing piles without initial applied load was reported. In reality, pile group in soft clay behave as floating pile group and subjected to initial applied load from superstructure. In the presence of initial applied load, soil surrounding the pile group experiences higher stress level before the commencement of adjacent excavation. Authors in [7-8] conducted centrifuge tests to investigate the effects of an unpropped excavation on the behaviours of nearby single piles and pile groups in dry dense Toyoura sand. It was found that the distance from the pile to the retaining wall and pile head Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2485-2492 2486 www.etasr.com Soomro et al.: Single Pile Settlement and Load Transfer Mechanism due to Excavation in Silty Clay conditions had a large influence on the induced pile bending moment and lateral deflection. Authors in [9] reported the results of three centrifuge tests which were carried out to study the effects of a multipropped deep excavation in-flight on the behaviour of single piles in dry Toyoura sand. Piles were laterally restrained in terms of rotation and deflection right at or above ground surface in the three different tests. It was concluded that lateral restraints imposed on the pile head have a significant influence on induced pile bending moment. Induced bending moment due to excavation can exceed the pile bending capacity. Authors in [10] developed design charts to compute the lateral behaviour of a single pile adjacent to deep excavation in soft ground. They performed two staged analysis considering plane strain conditions and linear elastic soil model. Similar work was conducted in [11] in which authors used finite element method. In both these studies lateral response of single pile was investigated. The settlement behavior of pile and development of excess pore water pressure and consolidation settlement were not investigated. Authors in [12] proposed an analytical method to investigate reduction of capacity and increase in settlement of nearby pile during excavation. It was reported that pile settlement due to excavation depends on the percentages of end bearing and shaft friction of the pile, the soil movement pattern, and the distribution of the maximum shaft friction with depth. However, shaft resistance in these methods is calculated on the assumption that horizontal stress acting on the piles does not change during excavation. This assumption may not be valid and the pile settlement may be underestimated using the preceding methods, leading to a non-conservative prediction. Most of the previous studies focused on the lateral response of single pile foundation. The vertical response of floating pile group combined with working load adjacent to deep excavation in soft clay has not been studied. Moreover, in soft clay progressive change in excess pore water pressure with time and consolidation settlement of pile group were not investigated. It is well recognised that the stress-strain relationship of soils is highly nonlinear even at very small strain. The stiffness of most soils decreases as strain increases and depends on the recent stress or strain history of the soil [13-14]. Owing to non- linear soil behaviour, an excavation can cause reduction to the ground stiffness. Therefore, it is vital to investigate the pile responses to adjacent excavation in silty clay. To obtain a satisfactory numerical model of the single pile responses to excavation-induced stress relief, the analysis must take account of the small strain non-linearity of soil. In view of the aforementioned issues, this study aims at investigating the settlement and load transfer mechanism of an existing single pile due to an adjacent excavation in saturated silty clay. To achieve these objectives, a three-dimensional coupled- consolidation numerical modelling was carried out. Settlement, axial load distribution along the pile, stress changes and excessive pore pressure generation during excavation are reported and discussed. II. THREE-DIMENSIONAL COUPLED CONSOLIDATION ANALYSIS To gain new insights into single pile responses to a nearby multipropped excavation in saturated silty clay, this study conducts a three-dimensional coupled consolidation numerical analysis. Figure 1(a) shows the elevation view of the configuration of numerical simulation in which a multipropped excavation was carried out adjacent to a single pile. The final depth of the excavation (He) was 10 m. The embedded length (Lp) and diameter (dp) of the pile were 18 m and 0.8 m, respectively. The modelled pile represents a cylindrical reinforced concrete (grade 40, reinforcement ratio = 1) with a bending moment capacity of 800 kNm. The clear distance between diaphragm wall and the pile was 3.0 m (0.3He). The excavation was supported by 0.6m thick diaphragm wall. The ratio of wall penetration depth to excavation depth is typically 0.5-2 in engineering practice [15-16], and thus a value of 0.5 was adopted in this study. The retaining wall was supported by three levels of props, at 1.0, 4.0 and 7.0 m depth. The props were modelled as soft with axial rigidity of 81 × 103 kNm [15]. Horizontal spacing of props was 10 m. Ground water table Working load (1073 kN) Existing single pile Silty clay Prop 1 (1.0 m) Prop 2 (4.0 m) Prop 3 (7.0 m) P la n e o f sy m m et ry Excavation level 1 (3.0 m) Excavation level 2 (7.0 m) Formation level (10.0 m) Embedded length = 18 m Diameter = 0.8 m Diaphragm wall Thickness = 0.6 m Depth = 14.0 m 6.0 m (a) Monitoring section 50m 20m P la ne o f sy m m et ry 0.6 m thick diaphragm wall Existing single pile Prop Prop 10m 6m (b) Fig. 1. Configuration of the numerical run (a) elevation, (b) plan view It is worth noting that, in reality, high-rise buildings are unlikely to be built on a single pile. This hypothesised study is a virtual case [3, 17-18]. This simplification is made to understand the settlement and load transfer mechanism clearly. Figure 1(b) illustrates the plan view of the configuration of the numerical simulation. The length of the excavation is 12 m. Due to symmetry, only half of the excavation was simulated. A monitoring section was selected at the transverse centreline of the excavation. In addition to this simulation, a pile load test (L) was conducted numerically in “greenfield” conditions (i.e., without excavation) to obtain the ultimate capacity of the pile in silty clay. Based on this, the working load was then calculated with a factor of safety of 3.0. The obtained working load was applied to the pile in the analysis simulating excavation. Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2485-2492 2487 www.etasr.com Soomro et al.: Single Pile Settlement and Load Transfer Mechanism due to Excavation in Silty Clay A. Finite Element Mesh and Boundary Conditions Figure 2 shows an isometric view of a finite element mesh used to analyze the soil-pile-excavation problem. The size of the mesh for each numerical runs is 50 m × 20 m × 40 m. These dimensions were sufficiently large to minimize boundary effects in the numerical simulation as further increment of the dimensions of the finite element mesh did not lead to any change in the computed results. Regarding the element size in the mesh, it is found that further halving of the adopted mesh size only leads to a change of computed results of no more than 0.2%, suggesting the mesh is sufficiently fine. Eight-noded hexahedral brick elements were used to model the soil, the pile and the diaphragm wall, while two-noded truss elements were adopted to model the props. Roller and pin supports were applied to the vertical sides and the base of the mesh, respectively. Therefore, movements normal to the vertical boundaries and in all directions of the base were restrained. The water table was assumed to be at ground surface. Initially, the pore water pressure distribution was assumed to be hydrostatic. Free drainage was allowed at the top boundary of the mesh. The excavation process was simulated by deactivating soil elements inside excavation zone. In the meantime, the truss elements representing the props were activated. Single p ile Diap hragm wall Props Pin support Roller sup port Fig. 2. Three dimensional finite element mesh and boundary conditions B. Constitutive Model and Model Parameters Used in Finite Element Analysis The basic hypoplastic model was developed to capture the non-linear behavior (upon monotonic loading at medium to large-strain levels) of granular materials [19-20]. The basic model consists of five parameters (N, λ*, κ*, φc and r) . Parameters N and λ* define the position and the slope of the isotropic normal compression line in the ln(1+e) versus lnp’ plane [21], e is the void ratio and p’ is mean effective stress. The parameter κ* defines the slope of the isotropic unloading line in the same plane. φc is the critical state friction angle and the parameter r controls the large strain shear modulus. To account for the strain-dependency and path-dependency of the soil stiffness (at small strains), authors in [22] further improved the basic hypoplastic model by incorporating the concept of intergranular strain. The intergranular strain concept requires five additional parameters (R, βr, χ, mT and mR): R controlling the size of the elastic range, βr and χ controlling the rate of stiffness degradation. The parameters mT and mR control the initial shear modulus upon 180° and 90° strain path reversal, respectively. In this hypothesised study, the parameters for silty clay were adopted the ones in [23]. All the model parameters for silty clay reported in [23] are summarised in Table I. The coefficient of lateral earth pressure at rest, Ko is estimated in [24] by the following equation:      sin0 sin1 OCRK (1) TABLE I. MODEL PARAMETERS OF SILTY CLAY [23] Description Value Effective angle of shearing resistance at critical state: ’ 33o Parameter controlling the slope of the isotropic normal compression line in the ln(1 + e) versus lnp plane, * 0.103 Parameter controlling the slope of the isotropic normal compression line in the ln(1 + e) versus lnp plane, * 0.015 Parameter controlling the position of the isotropic normal compression line in the ln(1 +e) versus lnp plane, N 1.31 Parameter controlling the shear stiffness at medium- to large- strain levels, r 0.3 Parameter controlling initial shear modulus upon 180 strain path reversal, mR 12 Parameter controlling initial shear modulus upon 90 strain path reversal, mT 12 Size of elastic range, R 2×10-5 Parameter controlling the rate of degradation of the stiffness with strainr 0.09 Parameter controlling degradation rate of stiffness with strain 0.7 Initial void ratio, e 0.7 Dry density (kN/m3) 1615 Coefficient of permeability, k (m/s) 1×10-9 The concrete pile, the diaphragm wall and the props were assumed to be linear elastic with Young's modulus of 35 GPa and Poisson's ratio of 0.25. The wall thickness was taken as 0.60 m. The concrete unit weight was taken as 24 kN/m3. The parameters for the piles and the diaphragm wall are summarized in Table II. TABLE II. CONCRETE PARAMETERS ADOPTED IN FINITE ELEMENT ANALYSIS Description Value Young's Modulus, E 35 GPa Poisson's ratio,  0.3 Density,  2400 kg/m3 C. Numerical Modelling Procedure Each numerical analysis is modeled according to the following steps: Step 1: Set up the initial boundary and initial stress conditions (i.e., static stress conditions with varying K0 with depth). Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2485-2492 2488 www.etasr.com Soomro et al.: Single Pile Settlement and Load Transfer Mechanism due to Excavation in Silty Clay Step 2: Activate the brick elements representing single pile (modeled as “wished-in-place”). Step 3: Apply the working load (determined from numerical pile load test) on the pile. Step 4: Allow excess pore pressure, which generated in result of the application of working load on the pile, to dissipate. Step 5: Activate the brick elements representing the diaphragm wall. Step 6: Staged multipropped excavation was simulated as described above. After excavating to 3 m depth, the first level of props was installed at 1 m below the ground surface. Soil was then excavated to 7 m below the ground surface, followed by installation of the second level of props at 4 m depth. Finally, excavation was extended to the target level of 10 m depth with installation of third level of props at 7 m. III. INTERPRETATION OF COMPUTED RESULTS A. Determination of Working Load for the Pile The objective of this parametric study was to investigate single pile load responses to an adjacent excavation. Prior to the simulation of the excavation, it was necessary to determine the ultimate axial load carrying capacity of the pile group. The working load can then be obtained using a factor of safety (FOS) of 3.0. Therefore, a numerical pile load test was carried out on a different finite element mesh to obtain the load settlement relationship and the capacity of the single pile without excavation. The load applied on the pile was gradually increased to 4 MN over a period of 24 h. The resulting pile load-displacement curve for the simulated pile is shown in Figure 3. The ultimate axial load capacity was determined based on a displacement-based failure criterion in [25]. This failure criterion is expressed as follows: ,max 1 0.045 2 h p ph p p p P L d A E    (2) where δph,max is the maximum pile head movement which defines the ultimate load, Ph is the pile head load, Lp is the pile length, Ep is the pile shaft elastic modulus, Ap is the cross- sectional area of the pile, and dp is the pile diameter. Based on the failure criterion, the ultimate bearing capacity of 3.22 MN was calculated. With a factor of safety (FOS) of 3.0, the working load was determined to be 1.07 MN. Owing to the applied working load, the pile settled by 0.7% dp (Figure 3). B. Progressive Pile Settlement and Apparent Loss of the Pile Capacity During Excavation Figure 4 shows the incremental settlement (Sp) of the pile with different excavation stages. Construction stages of the excavation are indicated by the depth h the ground surface. Sp and h are normalized by the pile diameter (dp) and the final excavation depth (He), respectively. The measured induced pile settlement due to excavation in centrifuge modeling reported in [9] is also shown for comparison. The centrifuge test was carried out to investigate the effects of a multipropped deep excavation (final depth of excavation= 8 m in prototype) in- flight on the behavior of single pile (diameter of the pile= 1.25 m in prototype) in dry Toyoura sand (i.e., Dr= 70%). The embedded depth of the pile and clear distance between the pile and model wall is 20 m and 3 m, respectively. A linear increase in Sp of the pile was observed with the first two excavation stages (h/He=0.7). However, the pile settled at increased rate due to subsequent excavation stage (h/He=1.0). This may be attributed to two reasons. First, owing to reduction of shaft resistance due to excavation induced stress release, the pile had to settle substantially to further mobilize end-bearing (explained below). Second, the most part of the pile is placed within the boundaries of a major influence zone due to excavation-induced ground movement. Similar settlement characteristics of a single pile due to excavation observed from the centrifuge test reported in [9]. As expected, induced settlement of the single pile in medium dense sand was less than that of single pile in stiff clay. This is because the sand stiffness is higher than that of clay. The induced pile settlement due to excavation was 9.5 mm (1.2% dp). Fig. 3. Computed load settlement curve from the pile load test without excavation Total pile settlement (including settlement due to working load and excavation) was 15.1 mm (1.9% pile diameter). This value satisfies a reliability-based serviceability criterion (56 mm), which was developed based on 95 composite (reinforced concrete and steel) buildings subjected to settlement. This conclusion may not be applicable to scenarios in which the ground conditions or excavation depth are different from those adopted in this study. Considering pile capacity is often interpreted as using settlement criteria, and pile settlement due to excavation may be correlated with an apparent loss of pile capacity (ALPC). Therefore, the excavation induced settlement can be regarded as an additional load on the pile head. In this study, settlement of the pile due to the applied working load (before excavation) was about 0.7% dp (5.6 mm) (see Figure 3). The excavation resulted in an additional pile settlement of 1.2% Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2485-2492 2489 www.etasr.com Soomro et al.: Single Pile Settlement and Load Transfer Mechanism due to Excavation in Silty Clay dp (9.5 mm). With this additional settlement, the pile behaved as if it were loaded by an additional axial force of 1.0 MN on the pile cap (based on load settlement curve shown in Figure 3). The additional load accounts for 31% of the ultimate pile capacity of the pile (based on the failure criterion presented in [26]). In other words, ALPC of 31% was due excavation. It should be noted that the reduction of pile capacity does not mean that the ultimate load capacity is physically reduced due to the excavation. It instead indicates the violation of the serviceability limit state of the pile as a result of the excavation-induced pile settlement and designers should take measures to control or mitigate this adverse effect. C. Computed Ground Deformation Mechanism and Deviatoric Strain Figure 5 shows computed incremental displacement vectors due to excavation. Computed incremental shear strain that the excavation induced stress released is also superimposed in the Figure. It can be observed that soil on the retained side moved towards the excavation, whereas the soil underneath the excavation heaved upwards due to vertical stress relief inside the exaction zone. Furthermore, it can be seen that excavation- induced ground movement is mainly bounded by a wedge at an angle of 450 drawn from the wall toe. In addition to ground movement, significant shear strain (up to 1.44%) is also developed extending nearly at 450 from the wall toe. It can be seen that two-thirds of the pile shaft (0 ≤Z/Lp ≤0.67) are located inside and pile toe is located outside the sliding wedge zone of ground movement. Because the two-thirds of the pile shaft are affected by the ground movement due to excavation induced stress release, shaft resistance decreased. To maintain equilibrium, the pile had to settle resulting in increment in end- bearing. On the other hand, the pile toe is located within the region where shear strain is developed. As illustrated, the presence of the pile increases further the soil shear strain near the pile toe, as more end-bearing has to be mobilised to support vertical load transferred from the pile shaft to the pile toe after excavation. D. Changes in Axial Load Distribution Figure 6 illustrates the axial force distribution along the pile with normalised depth (Z/Lp) below the ground surface during excavation. The axial load distribution before excavation (after applying the working load) is also included in the figure as a reference. Before excavation, the pile carried approximately 70% of the working load (1073 kN) with its shaft resistance and the remainder with its end-bearing resistance. It can be observed that excavation-induced reduction in normal stresses to the pile shaft and downward soil movement caused an increase in axial load along the entire length of the pile at different excavation stages. At the end of the first excavation stage, the maximum increment in the axial force (14% of that at working load) was computed at Z/Lp = 0.43, which is above the formation level. By inspecting the axial load distribution after the first excavation stage, it is observed that along the upper portion of the pile (Z/Lp ≥0.2), the shaft resistance decreases to zero. Consequently, the load transferred to the lower portion of the pile. To maintain equilibrium, the pile had to settle (see Figure 5) to further mobilise the end-bearing and shaft resistance along the lower portion (Z/Lp>0.2). As the excavation reached at final depth, it caused a further reduction of the normal stresses to the pile shaft. Consequently, soil settled more than the pile resulting in negative skin friction (NSF) along the upper portion of the pile (0 ≤ Z/Lp ≤ 0.6). This suggests that this portion of the pile is subjected to “dragload” by the surrounding soil. This caused the pile to settle more than that due to pervious excavation stages (see Figure 6). To maintain vertical equilibrium of the pile, the soil surrounding the lower part of the pile (Z/Lp>0.6) resisted its settlement by mobilising positive skin friction (PSF) at the pile-soil interface and end-bearing resistance at the toe of the pile. Owing to end of the excavation, the end-bearing increased to 74%. Fig. 4. Computed (this study) and measured ([9]) normalised pile head settlement during excavation % mm E x ca v at io n 45o Sliding wedge zone of ground movement Fig. 5. Computed incremental ground movement and shear strain at the excavation end. E. Mobilised Shaft Resistance During Excavation To substantiate the previous discussion, the mobilized shaft resistance in the pile after the application of the working load (before the excavation) and the excavation is shown in Figure 7. In the figure, the depth below the ground surface (Z) is Engineering, Technology & Applied Science Research Vol. 8, No. 1, 2018, 2485-2492 2490 www.etasr.com Soomro et al.: Single Pile Settlement and Load Transfer Mechanism due to Excavation in Silty Clay normalized by the pile length (Lp). The computed average mobilized unit shaft resistance f(Z) at various depths was calculated based on the following equation:    Z F f Z s Z    (3) where ΔF is the difference between the computed axial loads at two consecutive depths, ΔZ is the vertical distance between the two consecutive depths, and s is the perimeter of the pile. As expected, after the application of the working load, the positive shaft resistance (soil supporting the pile) is mobilized along the entire length of the pile. The mobilized positive shaft resistance increased along pile depth with maximum shaft resistance of 31 kPa near the pile toe. It can be seen from the Figure that as excavation depth increases, mobilized shaft resistance decreased above the formation level (Z/Lp<0.5). On the completion of the excavation (h/He=1.0), the mobilized shaft resistance at the upper portion of the pile (Z/Lp<0.5) decreased to negative and increased below formation level (0.5