Analysis on Deformation and Stress Characteristics of a Multibraced Pit-in-Pit Excavation in a Subway Transfer Station

As an important part of subway transfer station construction, pit-in-pit (PIP) excavation has always been a key object of monitoring and control. In this study, taking the PIP excavation project of Hefei metro line 4 and line 7 transfer station as the background, combined with on-site monitoring and numerical simulation, the variation law of lateral wall displacement, ground surface settlement, and strut force during the PIP excavation were analyzed. &e results showed that the maximum lateral deformation of the pile caused by the excavation of the external pit accounted for 80%–90% of the total deformation and the surface settlement accounted for 70% of the total settlement. &e excavation of the inner pit only made the maximum lateral wall displacement of the outer pit and the surface settlement behind the wall increase slightly, the growth rate tended to zero, and the maximum lateral deformation depth was all above the excavation surface, which indicated that, for this project, the inner excavation had little effect on the outer pit retaining structure. In addition, the increased exposure time of the soil at the bottom of the pit and the presence of the corner effect will cause the further development of the structure displacement and the surface settlement. Based on strut force measurements, using the tributary area load distribution procedure, it was obtained that 0.3cHe was used as the upper limit index of strut force in the region, where c is the weight of the soil and He is the excavation depth of the foundation pit. &e research conclusions were helpful for the construction and sustainable development of the PIP project.


Introduction
With the continuous development of the social economy and increasing urban traffic pressure, many large and medium-sized cities have to build underground transportation networks to alleviate surface traffic pressure. e subway transfer station serves as a node of the underground transportation network, and its special supporting structure in the PIP will increase the construction risk. e study of the degree of soil disturbance and the deformation law of the supporting structure in the PIP excavation is of great significance for the protection of the surrounding environment. e apparent pressure diagram (APD) proposed by Terzaghi [1] and Peck [2] is widely used in sand, hard-cracked clay, and soft-to-medium clay for braced excavations to estimate the magnitude and distribution of strut loads. Based on a finite element analysis, some scholars have studied the influence of the deep excavation on existing structures [3], the basal heave stability for excavation in spatial variable soils [4], and the distribution of axial forces under different wall stiffnesses in different soil layers [5,6]. Wang et al. [7] proposed a braced excavation back analysis method that considered the three-dimensional effect based on the plane-strain assumption. Mohamed et al. [8] used the finite difference code to simulate a three-dimensional excavation and compared it with measured data. In addition, Mohamed et al. also briefly remarked on the prediction of excavation behavior by means of 2D compared with 3D numerical analysis results. Rich data and achievements have also been provided by experimental studies [9,10], and field measurement analyses of excavation cases [11][12][13][14][15][16][17][18][19][20] also enable the further development of foundation pit excavation research.
e above studies are based on the fact that a foundation pit is excavated to a certain depth. Considering the requirements of building foundation structures, elevator shafts, basements, and subway transfer stations, the inside of the foundation pit often needs to be excavated twice or even multiple times, and its cross-sectional form is difficult to guarantee the square rule. e phenomenon of difference in depth inside the foundation pit is becoming more and more common. Combining the theory of ground loss with empirical equations, Sun et al. [21] proposed a prognoses method on curve and the maximum value of ground subsidence which is suitable for Beijing area. And, based on the hypothesis of theoretical excavation line and function mass，the mechanism of difference in surrounding strata deformation was proposed. Hou et al. [22] adopted the hardening-soil finite element model with unloaded soil to study the influence of the ratio β of internal depth to bottom width on the supporting structure and soil deformation. Based on the limited equilibrium theory and the hypothesis of planar sliding plane, Han et al. [23] deduced four formulas for calculating passive earth pressure with cohesion force and different shapes of sliding soil mass taken into account, and corresponding mathematical expressions of shear failure angles are given. e above studies on PIP excavation are based on specific assumptions using theoretical and numerical analysis methods. However, there are few studies on detailed field measurement data during the entire construction process. e research of the field measured data is often more representative than the theoretical and numerical analysis, and the research results are more instructive to the construction of the PIP project.
In this study, taking a PIP excavation in a subway transfer station in Hefei as the background, the lateral wall displacement, ground settlement, and strut force during the whole construction process were monitored in detail.
rough the analysis of monitoring data and comparison with numerical simulation results, the spatiotemporal evolution of the support structure deformation and ground settlement during PIP excavation is discussed. As the capital city of Anhui Province in southeastern China, Hefei has multiple subway lines under simultaneous construction. e impact assessment of subway construction on the surrounding environment is particularly important [24]. ere are few studies on the PIP excavation in this area, and it is hoped that the results of this study will provide guidance for the follow-up subway construction.

Site.
e excavation project used in this study was the transfer station of two orthogonally intersecting subway lines. Figure 1 shows a plan view of this site. e project for this station consisted of two parts. e outer excavation is the pit of the No. 4 line (first construction), the excavation plane is 261.8 m long × 21.9 m wide, and the excavation depth is 17.9 m. e inner excavation is the pit of the No. 7 line (second construction), the excavation plane is 176.9 m long × 23.1 m wide, and the excavation depth is 7.9 m. Because the south end of the No. 7 line intersects the standard section of the No. 4 line, a nonrectangular excavation surface is created at the north side of the intersection, which makes the design and construction of the PIP support structure more difficult.
e excavation examined in this study used the segmental earthwork transfer excavation method, and the excavation sequence is shown in Figures 1 and 2. e horizontal section length is based on the structural section and is controlled between 15 and 30 m. e vertical stratification is divided according to the mechanical excavation ability and the strut positions. e outer excavation is divided into four layers from top to bottom. e construction was conducted from the two ends to the middle, and the excavation was stopped when excavated to the intersection position. e inner excavation was divided into seven floors from top to bottom, and construction began from the intersection and moved from the two ends to the middle.

Geology.
Prior to the construction, the soil layer distribution of the site was obtained by drilling a pit location.
ere were six strata at the excavation site of the station: fill, medium clay, stiff clay, fully weathered sandy mudstone, strong weathered sandy mudstone, and moderately weathered sandy mudstone. e final excavation position of the outer pit was located in the stiff clay layer, and the final excavation position of the inner pit was located in the strong weathered sandy mudstone. e soil profile of the station intersection is shown in Figure 3. Before the PIP excavation, the geological survey personnel obtained the basic physical parameter indexes of the typical strata within the scope of the project through borehole sampling and laboratory tests, as shown in Table 1. e groundwater of this station consisted primarily of perched water and bedrock fissure water. e perched water was primarily distributed in the surface silty clay filling. e bedrock fissure water was primarily distributed in the fully weathered sandy mudstone and the strongly weathered sandy mudstone, and the single hole water inflow was generally less than 10 m 3 /d. e water richness was between extremely poor and poor. Figure 4 shows a plan view and two cross-sectional views of the north-south direction at the intersection of the two subway lines. Due to the use of the segmental earthwork transfer excavation method, in the excavation sequence shown in Figures 1 and 2, the outer excavation (line 4) had only the crossover portion of the soil during the excavation of the intersection, and the excavation of the east and west sides was completed, while the inner excavation (line 7) was not excavated. ree layers of concrete struts were used on the east and west sides of the intersection in the outer pits. Because the excavation depth of the southern end was different from that of the standard section, three concrete struts and two steel bracings with an external diameter of 609 mm were used on the south side of the intersection, and four concrete struts were used to support the north side of the intersection. e concrete strength grade was C30, and the steel was selected from        e support distribution, support material, and support section geometry parameters in the four directions are summarized in Table 2.

Advances in Civil Engineering
e excavation at the intersection was very deep, the shape of the plane was irregular, and the strut distribution was complicated. erefore, the construction of this portion of the structure was considered to be a risky link. e detailed construction stages of the intersections are summarized in Table 3. Because the excavation was conducted in stages, for the same layer of support at the intersection, the support in different directions differed in the installation time. For ease of analysis, Table 3 gives a uniform support installation time as a reference. e outer excavation was divided into east-west sides (first excavation) and intersection position (secondary excavation), which were divided into four layers of progressive construction according to the foregoing. To facilitate the study of the deformation characteristics of the excavation during different stages, the construction of outer excavation was divided into stage II and stage III, which represented the east-west excavation and the cross-position excavation, respectively.

On-Site Monitoring
3.1. Layout of Monitoring Points. By setting monitoring points at different positions in the four directions of the intersection, the lateral wall displacement, strut forces, and surface settlement during the different excavation stages were analyzed. To facilitate an analysis of the measurements, a uniform measuring point number was adopted, and the arrangement and numbering of the measuring points at the intersecting positions are as shown in Figure 5. e depth displacement of the pile body is reflected by the angle change of the inclinometer tube embedded in the internal portion of the pile. e buried depth of the inclinometer tube is equal to that of the retaining structure. It was fixed on the reinforcing cage and placed in the groove together and buried in the concrete. e axial force of the concrete support was calculated from the measured data of the steel bar stress meter embedded in the support section, assuming that the steel strain at the same section was equal to the concrete strain. Ground surface settlement is obtained by setting surface piles and observing with a leveling instrument. e ground settlement monitoring lines in different directions include 4 measuring points with a distance of 4 m. e accuracy of each instrument is shown in Table 4.

Measured Results of the Lateral Wall Displacement.
In this study, the representative diagonal positions H12 and H35 were selected at the intersection, and the actual displacement curves of the wall during the whole construction process are drawn, as shown in Figure 6. During the idealized excavation process, struts were installed during the initial stage of the excavation to reduce the cantilever rotation of the wall. However, due to the excavation sequence and site conditions, the struts needed to be installed layer by layer. e excavation process began from the cantilever state of the wall, and the support state transitions began from a single support to multiple supports. e first strut of the site was installed with a preexcavation groove. is was done, so there was no cantilever rotation of the wall during the initial stage of the excavation, which would have been in the form of a bulge.
From the measured results in Figure 6, it can be concluded that, during the excavation process, due to the lateral unloading, the outer retaining wall was laterally moved into the center of the excavation. In addition, the bulging of the external wall increased with the excavation depth, and the maximum displacement position (H m ) moved down with excavation depth. e wall on the north side had a large displacement at the H35 measuring point of approximately 18 mm, and the maximum position appeared at a depth of 12 m below ground surface. During the entire excavation process, the horizontal displacement at the top of the wall at the H35 measuring point was nearly zero, mainly because the pile at this position was constrained by the concrete support in two directions. In addition, it can be seen from the monitoring results that the displacement at the toe of the Advances in Civil Engineering pile was small, basically within 2 mm, indicating that the pile end was approximately fixed during the PIP excavation, and the pile has sufficient embedding depth. e wall displacement of the two measuring points increased slightly during the excavation of the inner pit. e maximum lateral wall displacement δ h−m of the H12 and H35 measuring points increased by 18.5% and 15%, respectively, when the inner excavation was completed.
is increased percentage was based on the completion of the outer excavation. In addition, H m remained basically unchanged during the inner excavation. It shows that the external excavation plays a leading role in the deformation of the outer pit retaining wall, and the inner excavation had little effect on the deformation of the outer wall. is is because the excavation position of the inner pit in this project is in the center area of the outer pit; the piles of the inner pit strengthen the soil at the bottom of the pit on the one hand and hinder the further development of the shear failure surface in the soil caused by the excavation of the outer pit on the other hand.

Measured Results of the Ground Surface Settlement.
According to the symmetry of the excavation plane, three lines S11, S12, and S21 were selected in this study. S11 was used to monitor the ground settlement behind the wall at the deepest excavation depth. S11-1 means the measuring point with a distance of 4 m behind the wall in No. 11 measuring point, which was the closest to the wall, and S11-4 was the farthest from the wall with a distance of 16 m. e measuring points S12 and S21 were located at the corners of the southeast and northeast sides, respectively. Figure 7 shows the on-site monitoring results of the ground settlement during the entire construction of the three measuring points.
It can be drawn from Figure 7 that surface settlement mainly occurs during the outer excavation. e maximum settlement of the S11, S12, and S21 survey lines at the completion of the outer excavation accounted for 65%, 72%, and 69% of the maximum settlement at the completion of the entire construction process, respectively. During the excavation of the inner pit, the settlement value of most of the measuring points tended to be stable, but the S21-3 further increased. e main reason was that the measuring point was affected by the simultaneous excavation of the standard section on the east side of line 4. e simultaneous excavation of two adjacent construction sections should be avoided, and the construction method of skip excavation should be adopted as much as possible.
In the final stage of the excavation of the inner pit (IV-5), the surface settlement of all the measuring points had  6 Advances in Civil Engineering further increased, because when the design elevation of the inner pit bottom was reached, the construction speed slows down, so that the time for the soil at the pit bottom to be exposed to the air increases. Due to the stress relaxation of the soil, the earth pressure in the active area behind the wall continues to increase with time, and the earth pressure in the passive area decreases continuously. e safety of the foundation pit gradually decreases, and the further lateral wall deformation causes the surface settlement behind the wall to increase. After excavation at the design elevation of the pit bottom, plain concrete cushions should be poured in a timely manner at certain intervals to reduce the exposure time and disturbance of the soil. Figure 8 shows the surface settlement profile of the different points at the completion of construction. It can be seen from the figure that the position of the maximum settlement value (V m ) of the surface behind the wall is different depending on the position of the measuring point. For the S11 measuring point with the largest excavation depth, the maximum settlement δ v−m � 11.89 mm. V m is located at S11-2, and V m is the closest to the wall of all the selected three measuring points. e S12 measuring point was affected by the excavation in the north and west S10-1 S10-2 S10-3 S10-4 S11-1 S11-2 S11-3 S11-4 S12-1 S12-2 S12-3 S12-4  Advances in Civil Engineering directions, and δ v−m � 13.84 mm and V m was located in S12-4. is demonstrates that the existence of the corner effect affected the deformation characteristics of the foundation pit excavation so that δ v−m increased and V m moved away from the wall. e surface settlement value of S21 was the largest at 15.9 mm. V m was between S11 and S12. is phenomenon occurred because there was a certain angle in space between the strut force and the earth pressure on the wall. Part of the strut force was used to balance the shearing force of the joint between the strut and the retaining wall, so that the support function of the struts was reduced.

Measured Results of the Strut Forces.
e west side data points SF17 and SF18 and the south side data points SF27, SF28, and SF29 were selected in the two orthogonal directions of the intersection. e monitoring results of the strut axial force during the PIP excavation are shown in Figure 9. e last digit of the monitoring point indicates the vertical brace position; for example, SF17-1 indicates the first brace of the SF17 measurement point.
For the two measuring points of the outer excavation, SF17 and SF18, the first strut wire was damaged during construction, and the second and third strut forces had an increasing trend throughout the construction process. e maximum strut force of the SF17 and SF18 data points was 2919 kN, both of which appeared in the second portions of the SF17 and SF18 data points, the SF17-2 and SF18-2 supports. For the inner excavation, in SF27, SF28, and SF29, the axial force of the first two supports was relatively gentle, while the axial force of the third support increased rapidly, especially during the inner excavation stage when the SF29-3 support axial force reached 3699 kN. Combined with the lateral displacement of the wall, it can be concluded that, during the excavation, the maximum lateral deformation position gradually moved down with an increase in the excavation depth, and the maximum displacement increased continuously. Under the coordinated action of structure and soil deformation, an increase in the excavation depth had more significant influence on the axial force of the bottom struts.

Details of the Numerical Simulation
Due to the limited measured data, the deformation characteristics of the entire excavation process could not be fully shown. Hence, the three-dimensional model of the finite element was used to further study the deformation laws at different positions in the deep excavation, as shown in Figure 10. Due to the low water content of the formation and the low groundwater level, the seepage effect was ignored in the simulation. By considering the calculation time and calculation accuracy, the mesh size of the excavation area was determined to be approximately 1.5 m, and the model boundary was a coarse grid with a size of 5 m. Because the focus of this study was on the intersection, the model reduced the number of grids by shortening the standard section distance for an excavation. Figure 10(a) shows a finite element model of a three-dimensional excavation. e overall model consists of 109,363 nodes and 220,964 elements. e soil was modeled using solid elements, and the support was modeled using node-coupled beam elements. e model considered the influence of the step method excavation on the deformation of the support system and excavated the soil on the east-west sides and the intersection position. Additionally, the soil interface at the intersection was in the form of three gradings (Figure 10(b)). e stability of the soil under self-weight   Advances in Civil Engineering was ensured, but it was closer to the actual construction conditions. Because the model could not simulate the retaining pile system in the actual project, the retaining pile was replaced with the diaphragm wall using the following equation that utilizes the principle of equal stiffness: where D is the pile diameter, t is the spacing of the center of the pile, and h is the converted wall thickness. e diaphragm wall was modeled using an 8-node quadrilateral plate element with six degrees of freedom per node. e parameters of the retaining pile and the thickness of the diaphragm wall after transformation are shown in Table 5. Figure 11 shows the distribution of the diaphragm wall after the equivalent of the retaining pile. e numerical simulation assumes that the soil follows the hardened soil (HS) model, because the model can simulate the hardening phenomenon where the increment of stress gradually decreases with strain. Geotechnical strength is described by triaxial test stiffness, triaxial unloading stiffness, and reloading strength, so it is more accurate than Mohr-Coulomb constitutive model to simulate the increase of foundation stiffness with compression. In addition, the model has functions of expansion limit to suppress excessive volume expansion and tensile strength limit to suppress excessive tension. e effective cohesion (C′), effective friction angle (φ ′ ), lateral pressure coefficient (K 0 ), and compression modulus (E S ) are shown in Table 1. Take on-site samples and use the GDS multifunctional triaxial apparatus ( Figure 12) [26] obtained the triaxial test stressstrain curves and consolidation test load-strain relationship curve of typical soil layers in Shanghai based on the triaxial consolidation drainage shear test, the triaxial consolidation drainage unloading-reloading test, and the standard consolidation test. e proportional relationship between these moduli is discussed based on the results of Brinkgreve's research. It can be expressed using the following equations: According to Janbu [27], the power exponent, m, for sand and silt associated with the modulus stress level, is generally 0.5. For cohesive soils, m ranges from 0.5 to 1. According to Bolton [28], the dilatancy angle (ψ) can be assumed to be (φ ′ − 30 ∘ ) for sand. For cohesive soils, ψ is generally taken to be 0. e parameters of the HS model used in the numerical simulation are shown in Table 6. Figure 13 shows the results of the three-dimensional numerical analysis of the deformation of the lateral wall during different excavation stages. e displacement of the wall at the H11 monitoring point with the largest excavation section was the largest, and the maximum lateral wall displacement (δ h−m ) of the outer excavation and the inner excavation was 32.77 mm and 54.34 mm, respectively, which were 0.19% and 0.32% of the excavation depth of the foundation pit, respectively. For the H12 monitoring point, when the excavation of the east-west soil was completed (stage II-7) δ h−m � 8 mm, the maximum deformation position, H m , of the lateral wall gradually moved downward with an increase in the excavation depth.

Lateral Wall Displacements.
is was consistent with the conclusions obtained from the measured data. When the outer excavation was completed (stage III-5), δ h−m of the H12 monitoring point was 19.1 mm, which is approximately 0.11% H e , less than 0.19% H e of the H11 measuring point. is indicated that an increase in the width of the excavation section will cause the displacement of the lateral wall to increase.
With the construction of the inner excavation, δ h−m and H m of each monitoring point increased, and H m was approximately equal to H e , which is consistent with the conclusions obtained by Ng et al. [17] based on the monitoring data. δ h−m at the completion of the inner excavation (stage IV-5) at the H12 and H35 points was increased by 12% Distance from the excavation (m) S11 S12 S21 Figure 8: Settlement profile of S11, S12, and S21 measuring points when the excavation was completed. and 7%, respectively, compared with the completion of the outer excavation (stage III-5). e significant increase in the displacement of the lateral wall of the H11 measuring point was primarily due to an increase in the excavation depth and a reduction in the strut stiffness. erefore, the H11 measuring point did not reflect the influence of the inner excavation on the deformation of the external excavation supporting system. It can be seen from the strut distribution in Figure 4(a) that the fourth and fifth layer supports of the S11 were both steel struts, and these were reduced in compressive strength and rigidity compared with the first three concrete supports. Specification [29] stipulates that the maximum horizontal displacement of the supporting structure is δ h−m ≤ 0.15% H e and δ h−m ≤ 30 mm. e numerical calculation results showed that when the outer excavation was completed and the construction of the inner excavation was conducted, the displacement of the wall at monitoring point H11 was greater than the control value.
Hence, a reinforcement plan should be developed in advance and on-site monitoring should be strengthened. en, when the displacement approaches the control value, needed reinforcement measures can be taken in time. Figure 13(d) shows the deformation curve of the H36 monitoring point at the inner pit retaining wall during different excavation stages. e excavation of the east-west soil did not affect the inner retaining structure, but the excavation of the soil at the upper part of the intersection caused a vertical unload that had a small effect on the   deformation of the inner excavation lateral wall. e installation of the 4th and 5th braces inhibited the development of a displacement of the wall at the inner excavation, and the maximum displacement of the lateral wall was close to the wall too. e deformation curve of the inner wall was different from the bulging of the outer wall, and it manifested more as tilting and deflection, but the maximum deformation was less than 10 mm, which was far less than the control value, and meets the requirements of the specification. Figure 14 shows the relationship between the maximum lateral deformation and the depth of excavation obtained  Advances in Civil Engineering from simulation results and measurements. It also summarizes the measured results of the deep excavation by Wang et al. [11]. It can be seen from the figure that both simulation and measurement in this study are included in δ h−m � (0.03% − 0.12%)H e . In other cases, as the depth of excavation increased, δ h−m grew faster and tended to 0.5%H e . However, the simulated and measured δ h−m growth rate of this site during the inner excavation was close to 0, indicating that the impact of the inner excavation on the diaphragm wall of the outer excavation was small.  Advances in Civil Engineering 13 Figure 15 shows the relationship between the maximum lateral deformation depth and the excavation depth obtained from the simulation results and measurements. It also summarizes the measured results of China [11]. e results show that, except for several data points during the initial stage, most of the data points in this study are below the line H m � H e , while the data from the Shanghai metro stations indicate that this occurs at a depth in the range of H m � H e and H m � H e + 10. e main reason for this difference was that the excavation in Shanghai uses prestressed props.

Ground Surface Settlement.
ere are generally two types of settlement profiles caused by excavation. One is the spandrel type, in which maximum surface settlement occurs very close to the wall when the excavation depth is small. e other is that as the depth of the excavation increases, the settlement profile becomes a concave type in which the maximum surface settlement occurs at a distance away from the supported wall.
Four surface settlement monitoring points corresponding to the lateral wall displacement monitoring points were selected to analyze the variations in surface settlement during different excavation stages. Figure 16 shows the simulation results for the four surface settlement monitoring points. e results show that, with excavation, the maximum settlement positions (V m ) of the different measuring points were not fixed but gradually moved away from the retaining wall, and the settlement profile changed from spandrel to concave. V m of S11 was the closest to the wall, followed by the S20 and S21 measuring points, and V m of S12 was the farthest from the wall. ese results are consistent with the conclusions obtained from the measured data. e settlement of S11 was the largest, and, since the depth of the excavation was the deepest h � 25.8 m, the maximum settlement value (δ v−m ) was 11.2 mm. e excavation of the east-west soil had little effect on the surface settlement after the intersection of the wall, and the surface settlement here was approximately 2 mm. Excluding the S11 measuring point, the internal excavation made δ v−m of S12, S20, and S21 increase by 59%, 17%, and 15%, respectively, compared with the completion of the external excavation. Because the S12 measuring point was close to the interface between the inner and outer excavations, the corner effect on the settlement of the ground surface behind the wall was significant. Hence, the settlement value fluctuated within 20 m behind the wall (Figure 16(b)), which was approximately 0.8 times the inner excavation depth. V m of S12 increased compared with other measuring points and was approximately 30 m behind the wall.
As the depth of excavation increased, the range of the settlement values continued to expand. e numerical simulation results were affected by the calculation accuracy. e area where the surface settlement curve reached a stable value away from the wall was considered an unaffected area. It can be concluded from Figure 16 that the outer excavation played a leading role in the surface settlement range and the inner excavation further increased δ v−m , which had no effect on the surface settlement range. e influence range of the excavation obtained from the numerical simulation was larger than 3.0H e , which was higher than that of Wang et al. [11] and Liu et al. [16], which were 2.0H e and 1.0H e , respectively.
Peck [2] reported a diagram of the normalized maximum settlement versus the distance behind the wall and divided it into three regions. Zone I was sand and soft-tohard clay, zone II was very soft-to-soft clay, and zone III was very soft-to-soft clay to a great depth. e simulated and measured data in this study and the three areas summarized by Peck are shown in Figure 17. It can be seen from the figure that the predicted and measured values found in this study are very small, and the data points all fall within zone I.

Apparent Earth
Pressure. To further analyze the stress level of the supporting axial force, the maximum load on each layer of strut was calculated using the apparent earth pressure diagram. In this study, the apparent earth pressure envelope during the excavation was estimated based on the maximum strut force of the on-site monitoring at the SF17 and SF18 measuring points of the outer excavation and the SF27, SF28, and SF29 of the inner excavation.
us, the relationship between the normalized apparent earth pressure (σ h /cH e ) and the normalized depth (z/H e ) was plotted, as shown in Figure 19. e limited measurement results in this study were compared with other cases in similar hard clays. Figure 19 includes the cases for Twine [19] and Li [30]. It can be seen from the figure that the apparent earth pressure is less than 0.5cH e , which is consistent with the characteristic line (red line) proposed by Twine and Roscoe. However, for most of the data in this study and in Li and Twine work, the apparent earth pressure value was less than 0.3cH e (blue line), with the exception of four values in Twine that exceeded the predicted range and one value that was equal to the predicted value. For the PIP excavation project in this region, it is more reasonable to use 0.3cH e as the upper limit index of the strut axial force. When this value  Figure 16: Numerical results of ground surface settlement at S11 (a), S12 (b), S20 (c), and H21 (d).
is exceeded, the engineers should be vigilant and review the measured data and consider taking reinforcement measures to guarantee construction safety.

Discussion
is study verifies the feasibility of the model in calculating the stress and deformation of stratum and structure during the PIP excavation through the comparison of the measurements and the numerical simulation results. It provides a guidance method for the construction of subsequent subway stations in this region. However, for the model and simulation method proposed in this study, there are still insufficient and required follow-up researches, summarized as follows: (1) In determining the soil parameters, this study was limited by the available undisturbed soil samples,   This site Twine and Roscoe [19] Li [30] The proposed and a semiexperimental and semiempirical method had to be used. In the succession of the studies, it is necessary to obtain accurate soil parameters through more undisturbed soil tests. In addition, the influence of parameter changes on the calculation results also needs further study. (2) In the foundation pit project, the steel brace has the advantages of light weight, low cost, and being conducive to construction compared with the concrete brace. Based on the calculation model in this study, one important future direction is to evaluate the feasibility of using steel braces with small spacing to replace concrete braces, which has an significant contribution to saving construction time and cost. (3) In view of uncertain external factors, such as the structural stress redistribution caused by local support fractures, the accuracy of the prediction of the method in this study is of great significance for the safe construction of the PIP project.

Conclusions
In this study, the lateral displacement, surface settlement, and strut force variations of retaining pile were analyzed based on measured field data. en, based on a triaxial consolidation draining shear test and previously reported research, the three-dimensional model parameters for simulating the deformation characteristics of a multibraced excavation were obtained. In addition, the deformation characteristics of different locations in the space were further studied. e measured results of the reported cases were also compared with the simulated and measured values from this study. e conclusions drawn from this research are as follows: (1) e excavation of the outer pit caused the maximum lateral deformation of the wall to gradually move down and the maximum displacement to increase. Both the measurement and the simulation results indicated that the maximum lateral deformation of the pile caused by the outer excavation accounts for 80%-90% of the total deformation. Affected by the location of the inner excavation, the increase in the displacement of the outer pit pile caused by the inner excavation in this project was small, within 20%. e simulated and measured maximum deformations in this study were included in (0.03%-0.12%)H e , the growth rate during inner excavation tended to 0, and the maximum lateral deformation depth was above the excavation surface, which indicated that, for this project, the inner excavation had little effect on the outer pit supporting structure. (2) e surface settlement mainly occurs during the outer excavation. e surface settlement caused by the outer excavation accounts for 70% of the total settlement. During the excavation of the inner pit, the settlement value of most of the measuring points tended to be stable. However, when the inner pit was excavated to the design elevation of the bottom of the pit, as the construction speed slowed down, the exposure time of the soil at the bottom of the pit increases, causing a further increase in surface settlement. During construction of the bottom of the pit, plain concrete cushions should be poured in time to reduce the exposure time and degree of disturbance of the soil. (3) Under the influence of the corner effect, the maximum settlement of the corner position was furthest away from the wall. e simulation results further reflected the existence of the corner effect, which caused fluctuations in the settlement value within 0.8 times of the excavation depth behind the wall. By plotting the normalized maximum settlement and the distance behind the wall, it can be concluded that the simulation results and the measured values were small. In addition, the data points fell within the range of zone I (sand and soft-to-hard clay), as summarized by Peck. (4) During the entire excavation process, the strut force of each layer increased with an increase in the excavation depth, and the axial force of the bottom support changed most significantly. It can be found by drawing the apparent earth pressure diagram that, for the PIP excavation in the local area, it is more reasonable to use 0.3cH e as the upper limit index of the strut axial force. When this value is exceeded, engineers should be vigilant and be prepared to take reinforcement measures.

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

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.