Abstract

In this paper, a three-dimensional numerical parametric study was conducted to predict the deformation mechanism of a 20-storey building sitting on a (4 × 4) piled raft to an adjacent 25 m deep basement excavation. The influences of different excavation depths were investigated. An advanced hypoplastic sand model (which is capable of taking small-strain stiffness into account) was adopted to capture soil behaviour. The computed results revealed that excavation adjacent to a building resting on the piled raft caused significant settlement, differential settlement, lateral deflection, and interstorey drift in the building. With settlement due to working load (i.e., 4.8dp%), the total settlements of the building (7.8dp%) exceed the maximum allowable foundation settlement (i.e., 50 mm). In addition, substantial bending moment, shear forces, and changes in axial load distribution along pile length were induced. The findings from this study revealed that the building and pile responses are significantly influenced by the excavation depth.

1. Introduction

Underground space has become increasingly important in relieving the shortage of lands in megacities with rapid economic development and dramatic increase of population [1]. Deep basements for underground facilities such as subway stations are usually excavated near existing buildings and bridges. Basement excavation-induced stress relief and ground movements inevitably cause additional deformation and stress on adjacent pile foundations [2, 3]. It is highly likely that excavation-induced ground movement and stress relief may lead to severe damage and might result in collapse of structure. Thus, it is a major concern for designers and engineers to evaluate adverse effects of basement excavation on buildings. To ensure the safety and serviceability of existing pile foundations, several researchers have studied this complex problem by means of different approaches. Rocha [4] systematically described scale modeling for problems in soil mechanics in a 1 g gravitation field. Moncarz and Krawinkler [5] considered that if an adequate model correctly scales the primary features of the problem, the scaling relations between the prototype and model are not significantly affected. In the previous research study, Sawwaf and Nazir [6] analyzed the results of laboratory model tests on the influence of deep excavation-induced lateral soil movements on the behaviour of a model strip footing adjacent to the excavation and discussed the variation of the footing measured vertical settlements with different parameters. Fang et al. [7] conducted physical model tests of highway tunnel construction and examined the stability of the surrounding rocks regarding different caved zone-tunnel distances and dip angles of the coal seam. Yang et al. [8] investigated the mechanical behaviour of a typical jointed rock block located adjacent to an underground excavation by comparing the failure modes of the numerical simulations with the corresponding experimental physical jointed block test results. Shahin et al. [911] studied about the influence of existing building load on the deformation and earth pressure of the ground in shallow tunnelling by two-dimensional model tests and numerical analysis. Kusui et al. [12] conducted a large number of scaled-down tunnel experiments to investigate the response of unsupported walls to an increased stress field. Xu et al. [13] carried out a series of three-dimensional shaking table tests to investigate the mechanism and effect of seismic measures of the mountain tunnel using a scaled model based on a real tunnel. Furthermore, the centrifuge test was often used to study the effect of excavation construction, and Liyanapathirana and Nishanthan [14] and Zhang et al. [15] analyzed the excavation-induced pile behaviour using the finite element method and the centrifuge test. In terms of influence parameters of excavation construction, Yashiro et al. [16] conducted parameter analysis focusing on the overburden and the stiffness coefficient of the ground by anecdotal surveying and numerical analysis. Boone et al. [17] provided a modified approach for estimating potential damage and compared to case histories using construction data from a large braced excavation. Hsiao et al. [18, 19] considered the calculated settlement as the load in the context of reliability analysis and proposed a simplified model for evaluating the damage potential of a building adjacent to a braced excavation. Long [20] summarized the general trends and patterns based on some 300 case histories of wall and ground movements due to deep excavations worldwide. Among those numerical studies, however, responses of the framed building to excavation were reported in a limited number of studies. Most of these aforementioned research studies only examined the tunnelling effects on buildings founded on shallow footings [2124]. Therefore, there is a gap of systematic research on the responses of a high-rise building subjected to live load and resting on the piled raft to advancement of adjacent excavation. The paper aims to understand the responses of a 20-storey building to ground movement due to excavation-induced stress relief. To achieve this, three-dimensional finite element analysis was conducted. The building deformation (i.e., differential settlement, deflection, and interstorey drift) and piled raft (i.e., induced bending moment and shear forces in the piles) were presented and discussed.

2. Development of Three-Dimensional Finite Element Model

2.1. General

With the prime objective of investigating the excavation building interaction, three-dimensional finite element analysis was adopted in this study. An excavation was carried out in the close proximity of a 20-storey high rising building which is resting on a (4 × 4) piled raft foundation in sand. Figure 1(a) illustrates the general setup of the 20-storey building founded on a piled raft and a basement. In practice, the final excavation depth of the basement commonly varied from 12 to 26 m, and the thickness of the retaining wall varied from 0.3 to 1.0 m [1, 6, 25]. To simulate an actual construction in practice, the basement with a final excavation depth of 25 m (staged excavation) and a wall thickness of 0.3 m was simulated in this study. The depth of first and eighth (i.e., final) stages was 2 m and the depth of each of remaining stages (i.e., second to seven) was taken as 3 m. The ratio of wall penetration depth to excavation depth is typically 0.5–2 in engineering practice [1, 6, 26], and thus a value of 0.6 is adopted in this study. The retaining wall was supported by eight levels of props with vertical and horizontal spacing of 3.0 m and 5 m, respectively. The first level of props was installed at 2.0 m below the ground surface. The props are modeled as soft with axial rigidity of 81 × 103 kNm [2126]. A 60 m high, 15 m wide 20-storey, concrete building with three spans in each direction was selected in the study. The building was founded on a square raft, which was 15 m wide and 1.5 m thick supported by a group of piles in a 4 × 4 configuration with a centre-to-centre distance of 3.75 m (see Figure 1(b)). The diameter (dp) and length (Lp) of each pile were 1 m and 20 m, respectively. The clear distance between the closest pile and the diaphragm wall was 2.0 m (2dp). Since the prime objective of this parametric study is to assess the impact of excavation on the 20-storey building under dead and live loads, the live load of 5 kPa [27] was adopted in this study. With this adopted piled raft system, the numerical prediction showed that settlement of the building due to dead and live load was 45.2 mm (4.52%dp) which is within the limit of the allowable foundation settlement of 50 mm given by Skempton and Macdonald [28] and O’Brien [29].

2.2. Features of the Building

Figure 2 illustrates the building and its sections of structural components. The structural sections were adopted from Rasouli and Fatahi [30]. They designed the sections by conducting a routine design procedure by conducting analysis and design using SAP2000 [31] based on AS3600 [32], AS1170.1 [27], and AS1170.4 [33]. The details of all sections are summarized in Table 1. The specific compressive strength () and mass density of the concrete members in the building were considered to be 32 MPa and 2400 kg/m3, respectively, and the modulus of elasticity of concrete was estimated to be equal to 30.1 GPa.

2.3. Characteristics of the Numerical Model

Figure 3 shows the finite element mesh developed in Abaqus software. The size of the mesh numerical run was taken as 50 m × 50 m  × 60 m. The different structural components of the building (i.e., columns and slabs) were modeled by using two-node linear beam elements (B31) and four-node shell elements (S4). The mechanical behaviour of structural elements of the building was modeled using an elastic-perfectly plastic constitutive model to capture any possible inelastic behaviour of the structural elements; these structural elements could behave elastically until they reached yielding stress that was equal to the compressive strength of concrete (i.e.,  = 32 MPa) taken from the study of Shing and Tanabe [25]. Eight-noded hexahedral brick elements are used to model the soil, the pile, and the diaphragm wall, while two-noded truss elements are adopted to model the props.

Since the stress-strain relationship of soils is highly nonlinear even at very small strain and the stiffness of soil depends on the recent stress or strain history of the soil [3438], an advanced hypoplastic model was used to simulate the behaviour of sand in this study. The basic hypoplastic soil model requires eight material parameters to describe the nonlinear behaviour of sand (, hs, n, ed0, ec0, ei0, α, and β). Parameter is the soil frictional angle at a critical state. Parameters hs and n are used to describe the shape of limiting void ratio lines. Parameters ed0, ec0, and ei0 are reference void ratios of isotropic normal compression line, critical state line, and minimum void ratio line, respectively. Effects of soil relative density on peak frictional angle and shear stiffness are characterized by parameters of α and β. To consider the small-strain stiffness and stress path-dependent soil behaviour, the concept of intergranular strain is incorporated into the basic model [39], which includes additional five parameters (mR, mT, R, βr, and χ). By using the concept of intergranular strain, the modified hypoplastic sand model has the ability to capture effects of shear strain and stress path on soil stiffness [40]. These five parameters were obtained by fitting the stiffness degradation curves of Toyoura sand (see Figure 4) obtained from the stress-path triaxial tests carried out by Ng et al. [41]. Figure 4 shows the measured secant modulus degradation curves corresponding to four varying recent stress histories (i.e., θ = 0, 90, −90, and 180°), based on the GDS and CKC triaxial automated testing systems. It can be seen that the measured results in the tests for u = 180° by the GDS and CKC triaxial testing systems had fairly similar trends. Quantitatively, the percentage differences between the secant modulus measured by the two triaxial testing systems were less than 18%. The model parameters were taken from a study by Shi et al. [21]. They calibrated and validated all the parameters of the hypoplastic sand model against their centrifuge test results (which was performed to simulate excavation in sand). These parameters are summarized in Table 2. The lateral coefficient of earth pressure (K0) is estimated by the empirical equation proposed by Jáky [42], where is internal frictional of soil at the critical state. Since the internal friction angle of Toyoura sand at the critical state is 31°, the K0 value is estimated as 0.5.

2.4. Soil Structure Interactions and Boundary Conditions

While modeling excavation-pile-soil problem, one of the important aspects of modeling soil structure interaction (SSI) is to establish interaction between pile and surrounding soil. To incorporate the interactions between pile-soil and raft-soil, the surface-to-surface contact technique provided in the Abaqus software package was used [43]. In this approach, two surfaces in contact were assigned as the master and slave surfaces. The penalty approach was used for tangential contact, and the normal behaviour is modeled as hard contact with no normal relative displacement between the pile and surrounding soil. The interface was modeled by the Coulomb friction law, in which the interface friction coefficient (μ) and limiting displacement lim) are required as input parameters. A limiting shear displacement of 5 mm was assumed to achieve full mobilization of the interface friction equal to μ×, where is the normal effective stress between two contact surfaces, and a typical value 0.35 of μ for a bored pile was used in all analyses [35]. 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 excavation process will be simulated by deactivating soil elements inside the excavation zone. In the meantime, the truss elements representing the props will be activated.

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 unit weight of the concrete was assumed to be 24 kN/m3. The parameters for the piles and the diaphragm wall are summarized in Table 3.

2.5. Numerical Simulation

The numerical simulations of all the three cases were carried out in the following steps.(a)Step 1: initial geostatic stresses are generated in the mesh by applying gravity load and the coefficient of lateral earth pressure of 0.5.(b)Step 2: the piles are constructed, and then a raft is placed on the top of the pile group and soil deposit.(c)Step 3: the 20-storey building is constructed on the top of the piled raft. A live load of 5 kPa is applied on the floor slab of each storey.(d)Step 4: the brick elements representing the diaphragm wall are activated.(e)Step 5: staged multipropped excavation is simulated as described in Section 2.1. After excavating to 2 m depth, the first level of props is installed at 2 m below the ground surface.(f)Step 6: repeat step 6 to excavate the next stages and install props until the last stage of excavation (i.e., He = 25 m) is completed.

2.6. Validation of Numerical Modeling against Centrifuge Test Result

Before conducting any numerical parametric study, the measured ground surface settlement in centrifuge test results reported by Shi et al. [21] was used to validate the soil constitutive model and model parameters adopted in the finite element analyses. Figure 5 shows the comparison between measured and computed ground surface settlement trough developed on completion of the excavation. Both the ground surface settlement (S) and the distance behind the wall in the transverse direction are normalised by final excavation depth (He). It can be observed from the figure that the measured and computed settlement troughs show reasonable agreement. In this study, multipropped basement excavation is simulated. Thus, concave-shaped profile of ground settlement behind the wall is observed, as expected. The measured maximum soil settlement is observed at a distance of 0.375He, while the induced settlement is negligible when the distance from the wall exceeds 2.5He. When the normalised distance from the basement is 0.375He, the computed ground surface settlement at the end of basement excavation is 0.17%He, which is 18% larger than the measured one (0.143%He).

3. Interpretation of Computed Results

3.1. Settlement of Piled Raft

Settlement of the piled raft at its centre (Sp) induced due to different excavation stages is plotted in Figure 6. The piled settlement is normalised by the pile diameter. Excavation depths are indicated by h which is normalised by final excavation depth (He). It can be seen from the figure that bilinear increment in Sp was observed with increasing excavation depth. Sp increased linearly when excavation depth reached h/He = 0.68; however, the rate of induced settlement further increased until final excavation depth (h/He = 1.00). This was because of degradation of sand stiffness due to excavation-induced stress and the plastic strain generated surrounding the piles due to excavation-induced stress release. Since the initial excavation stages (0 < h/He < 0.68) are adjacent to the pile shaft, the shafts of the piles were affected by excavation-induced plastic strain. However, the toes of the piles were resting on the intact (from plastic strain) stiffer layer, whereas toes of the piles were severely influenced by the plastic strain as excavation stages passed beyond the pile toes (h/He > 1.00). Consequently, the piled raft experienced larger settlement due to excavation stages beyond h/He > 0.68 as compared to the initial excavation stages. On completion of the excavation, the piled raft settlement was 30 mm (3.0dp%). With settlement due to working load (i.e., 4.8dp%), the total settlements of the building (7.8dp%) exceed the maximum allowable foundation settlement (i.e., 50 mm) according to Skempton and Macdonald [28] and O’Brien [29]. Hence, the numerical prediction implies that the serviceability limit state of a building founded on the pile can be affected when excavation is carried out adjacent to the building.

3.2. Differential Settlement of Piled Raft

In addition to the discussion of piled raft settlement in the previous section, differential settlement induced in the building is an important parameter to investigate. During the excavation of nearby piled raft foundation, as expected, piles closest to the excavation are subjected to larger ground movement and stress release than those which are farthest from the excavation. As a result, differential settlement is induced in the piles. Figure 7 shows settlements of the piled raft induced at various positions along the centreline of the raft (see inset in the figure) at different excavation stages (i.e., h/He = 0.08, 0.32, 0.56, 0.80, and 1.00). The settlement before excavation is added in the figure for reference. It can be seen that after completion of the building construction (before excavation), the settlement occurred in the building was uniform through the centreline of the raft. However, as excavation was carried out, the raft corner closest to the diaphragm wall settled larger than that furthest from the wall. After the completion of the first excavation stage (i.e., h/He = 0.08), negligible differential settlement was predicted. However, as excavation depth increases, the differential settlement increases. This observation is attributed to excavation-induced stress release which is larger and smaller around the pile closest to and farthest from the excavation. The induced differential settlement caused the building permanent interstorey drifts (discussed in Section 3.3). The differential settlement in the raft of 3.6, 7.4, and 16.7 mm was induced on completion of excavation stages h/He = 0.56, 0.80, and 1.00, respectively.

3.3. Flooring Lateral Displacement and Interstorey Drifts

The numerical predictions in the previous section showed that the excavation resulted in differential settlement of the raft. This led to displace the 20-storey building laterally towards the excavation. Figure 8(a) illustrates the induced lateral deflection of the building along each storey at different excavation stages (i.e., h/He = 0.08, 0.32, 0.56, 0.80, and 1.00). It can be seen from the figure that lateral deflection increased with the height of the building. Negligible lateral deflection was computed after the first excavation stage (i.e., h/He = 0.08). However, as the excavation continues, the substantial lateral deflection was induced in the building. This observation is consistent with induced settlement in the piled raft. The maximum lateral deflection occurred at the top of the building (the roof level of 20th storey). The general trend of the lateral deflection profile is consistent with differential settlement. Since interstorey drift may induce distress in structural component of the building, it is necessary to predict the drift induced in the building due to excavation at different stages. Figure 8(b) illustrates interstorey drift of the building during excavation. The interstorey drift of the building can be defined as ratio of difference of deflections at two storeys to storey height. It is noted that interstorey drift shows similar trend as that of lateral deflection. The drift of 0.024%, 0.05%, and 0.11% was predicted on completion of excavation stages h/He = 0.56, 0.80, and 1.00, respectively.

3.4. Computed Ground Deformation Mechanism and Deviatoric Strain

Figure 9 shows computed incremental displacement vectors due to excavation. In addition, computed incremental shear strain of the excavation-induced stress release 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 45° drawn from the wall toe. In addition to ground movement, significant shear strain (up to 1.0%) is also developed extending nearly at 45° from the wall toe. A similar pattern and magnitude of excavation-induced shear strain were calculated from the measured centrifuge test in sand [21]. This indicates that the finite element analysis results are appropriate to analyze the excavation pile interaction problem.

It can be seen that all the piles are located at the zone of ground movement. Moreover, the piles closest to the diaphragm wall are subjected to larger ground movement as compared to farthest piles. This caused the differential settlement (see Figure 7) in the piled raft. Because the entire pile length is affected by the ground movement due to excavation-induced stress release, shaft resistance decreased. To maintain equilibrium, the pile had to settle which resulted in increment in end bearing (explained in Section 3.6). On the other hand, the toe of the closest pile is located within the region where shear strain is developed. As illustrated, the presence of the pile further increases the shear strain of soil near the pile toe, as more end bearing has to be mobilized to support vertical load transferred from the pile shaft to the pile toe after excavation.

3.5. Changes in Load Taken by the Raft during Excavation

After construction of the building resting on a piled raft system, some of the building load is sustained by the raft and remaining load is transferred to the piles. However, the load taken by the raft can be altered when excavation is performed close to piles. Figure 10 represents the change in load sharing by the raft at different depths of excavations. It can be seen that before excavation (after application of working load), about 5.9% of the building load (i.e., 57.6 MN) was carried by the raft and the rest of the load was transferred to sixteen piles. It can be observed that the load taken by the raft kept decreasing during excavation. As the excavation depth increases, rate of reduction of the load is the highest. This indicates that the building load was transferred to piles on completion of excavation. This can be attributed to the raft settlement and induced ground settlement. On completion of the excavation, the load resisted by the raft reduced to 4.85% of the total load.

3.6. Changes in Axial Load Distribution along Piles

As discussed in the previous section, the load taken by the raft altered substantially due to the excavation. It suggests that changes in axial load distribution along 16 piles can be induced. Since arrangement of the 16 piles in the piled raft system is in a square pattern (i.e., 4 × 4 piles), the piles at critical positions are selected for discussion in this section. The piles are designated as P1, P2, P3, and P4 (see the inset in Figure 11). Figure 11(a) compares axial load distribution along the lengths of piles P1 and P4. For comparison, load distribution along the length of each pile before the excavation is included in each figure. It can be observed from the figure that the load taken by the piles positioned at the corner of the raft (i.e., P1 and P4) is same, whereas piles in the middle of the raft shared equal load before excavation. Before excavation, the load taken by the piles P1 and P4 is 3517 kN. Both piles resisted the load by shaft resistance (36% of the load), and the remainder resisted by end bearing. On completion of the excavation, the axial load along the entire pile length of piles P1 and P4 increased. This indicates that shaft resistance decreased along the length of both the piles. This observation can be attributed to excavation-induced shear stain due to stress release at the midportion of the pile (discussed in Section 3.4). To support the load, shaft resistance was mobilized at the lower portion of the pile. Figure 11(b) compares axial load distribution along the length of piles P2 and P3. For comparison, load distribution along the length of each pile before the excavation is included in each figure. It can be seen from the figure that before excavation, the load transferred to each of the piles P2 and P3 was 3160 kN. The piles resisted the load by mobilizing shaft resistance (40% of load) along the lower portion the pile length, and remainder 60% of the load was resisted by end bearing. It can be seen from the figure that similar to the changes in axial load distribution along the pile length of P1 and P4, the axial load along the entire pile length of P2 and P4 increased on completion of the excavation. This load redistribution is ascribed to excavation-induced stress release and load transfer from the raft to piles. The end bearing of both the piles P2 and P3 increased by 10% to carry the transferred load from the raft.

3.7. Structural Response of Pile (Induced Bending Moment and Shear Forces along Piles)

Before excavation, the piles are horizontally confined by the surrounding soil and are subjected to negligible shear forces and bending moments. However, when excavation is carried out near piles and induces stress release in the ground, the substantial shear forces and bending moment can be produced in the piles. Thus, the shear forces and bending moment induced along piles due to excavation are important parameters to be predicted. Figure 12 illustrates the induced shear forces along normalised depth of four selected piles (i.e., P1, P2, P3, and P4) on completion of excavation. A highly nonlinear shear force distribution was predicted along each pile in each case. Due to rigid connection between the piles and raft, the maximum shear force generated is at the head of the piles. However, the magnitude of the maximum shear force was dependent on the location of the piles in the raft. The piles at the corners of the raft (P1 and P4) and the middle of the raft (P2 and P3) are subjected to positive and negative shear forces, respectively. The maximum shear forces developed at pile heads P1, P2, P3, and P4 were 730 kN, 600 kN, 550 KN, and 250 kN, respectively. This is because of excavation-induced stress release and ground movement towards excavation.

Since the piles are mostly designed to carry vertical load, the bending moment induced in the piles is another important parameter to investigate due to adjacent excavation. Figure 13 shows the induced bending moments along the normalised depth of four selected piles (i.e., P1, P2, P3, and P4) on completion of excavation. It is clear evidence from the figures that the maximum bending moment is induced at the head of each pile. This is because the piles are rigidly connected to the raft. The maximum negative bending moment (300 kNm) was predicted at the head of pile P1. To counterbalance the negative bending moment, positive bending moment was induced at the lower portion of the pile (Z/Lp > 0.2). The remaining piles (P2, P3, and P4) were subjected to positive bending moment at their heads. Since the pile toes were free to move, negligible bending moment was induced at the toes of each piles.

4. Conclusions

Based upon the analysis and results, the following specific conclusions are drawn:(1)Owing to degradation of soil stiffness, bilinear increment in Sp was observed with increasing of excavation depth. Sp increased linearly when excavation depth reached h/He = 0.68; however, the rate of induced settlement further increased until final excavation depth (h/He = 1.00). With settlement due to working load (i.e., 4.8dp%), the total settlements of the building (7.8dp%) exceed the maximum allowable foundation settlement (i.e., 50 mm) according to Skempton and Macdonald and O’Brien.(2)After the completion of the first excavation stage (i.e., h/He = 0.08), negligible differential settlement was induced. However, as excavation depth increases, the differential settlement increases. This observation is attributed to excavation-induced stress release which is larger and smaller around the pile closest to and farthest from the excavation. The differential settlement in the raft of 3.6, 7.4, and 16.7 mm was induced on completion of excavation stages h/He = 0.56, 0.80, and 1.00, respectively.(3)The excavation resulted in differential settlement of the raft. This led to displace the 20-storey building laterally towards the excavation which caused the building permanent interstorey drifts. The interstorey drift may induce distress in structural component of the building.(4)The load taken by the raft (about 5.9% of the building loads, i.e., 57.6 MN before excavation) can be altered when excavation is performed close to piles. On completion of the excavation, the load resisted by the raft reduced to 4.85% of the total load.(5)Apart from induced settlement, differential settlement, lateral movement, and interstorey drift in the building due to excavation, the piles are also subjected to substantial shear forces and bending moment along pile lengths.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge the financial support provided by the Quaid-e-Awam University of Engineering, Science & Technology, Sindh, Pakistan.