Two-Stage Analysis Method for the Mechanical Response of Adjacent Existing Tunnels Caused by Foundation Pit Excavation

: With the advancement of urban underground space networks, there has been a rise in foundation pit projects near existing tunnels. The construction of these foundation pits adjacent to existing tunnels can result in soil disturbance and stress redistribution, leading to additional deformation and internal force within the tunnels. This paper delves into the two-stage analysis method, outlining the calculation of additional stress in the initial stage considering various engineering factors and the methods for determining tunnel displacement and internal force in the subsequent stage. Through an engineering example and numerical simulations, the theoretical calculations were validated. The maximum displacement generated by the tunnel is − 4.85 mm and − 5.10 mm, respectively. The maximum error is only 5.9%, which confirms the validity of the theoretical approach. The analysis demonstrates that incorporating the unloading model of the bottom and surrounding side walls of the foundation pit is essential when calculating additional stress in the first stage. Moreover, the presence of engineering dewatering and double-hole tunnels can counterbalance the additional stress, with deviations of only 4.4% and 2.5%, respectively. In the second stage, factoring in the shear action and lateral soil action in the foundation and tunnel model enhances the accuracy of stress mode representation (accuracy increased by 18.8% and 29.3%, respectively). Additionally, accounting for the buried depth effect of the tunnel, soil non-uniformity, and foundation nonlinearity helps prevent excessive foundation reactions.


Introduction
The development of urban underground space has evolved from single, isolated structures to more complex and interconnected systems.As construction near existing underground structures increases, new foundation pits are becoming more common.Excavating these pits can impact the stress, deformation, and functionality of existing underground structures, potentially leading to engineering accidents.For instance, the collapse of a municipal channel foundation pit in Wuhan [1] resulted in multiple fractures of underground pipelines.Another example is the Nicoll Highway foundation pit accident in the Singapore Metro [2], which caused damage to a 100 m open-cut tunnel, damage to four fatalities, urban pipeline damage, and the collapse of the Nicoll fast track.To minimize these risks and prevent similar accidents, it is essential to conduct thorough research on how existing structures respond to foundation pit construction.
Only up to 2020, research methods for this issue, both domestically and internationally, could be classified into four main categories.The first category involves the numerical calculation method, where researchers establish models using numerical software, select relevant parameters, and conduct numerical simulation analyses of foundation pit construction [3][4][5][6][7].The second category is the field measurement method, which includes the field testing of existing engineering projects and the analysis of the collected data to understand the influence law [8][9][10][11][12].The third category is the model test method, where researchers conduct scale tests using similarity relationships, monitor real-time data during the test process, and analyze the data obtained through monitoring [13][14][15][16][17].The fourth category involves theoretical analysis, where researchers analyze previous data, deduce relevant formulas, or simplify tunnel analysis using relevant theories based on basic assumptions [18][19][20][21][22]. Theoretical analysis is considered simpler and faster compared to other methods, making it more suitable for preliminary design and guiding engineering construction.The two-stage analysis method is commonly used within theoretical analysis.Numerous scholars both domestically and internationally have conducted research on two-stage analysis methods.For example, Liu et al. [23] established a theoretical calculation model for the uplift deformation of the underlying tunnel caused by the excavation of foundation pits and found that the predicted results were basically consistent with the on-site measured results; Zhou et al. [24] used the incremental method to analyze the interaction between the support structure and the soil and found that the impact range of foundation pit construction on the underlying tunnel is about six times the excavation length; Zhou Zelin et al. [25] considered the effect of engineering precipitation on the deformation of existing tunnels, and the research results showed that excavation-related precipitation is a favorable factor in controlling the uplift of existing underlying tunnels.Although some research has been conducted through theoretical analysis methods, most of them are focused on underground tunnels, and there is relatively little theoretical research on the impact of excavation on side tunnels.
This paper analyzes and summarizes the calculation of additional stress in the first stage under various engineering factors, as well as the calculation method and research conclusions on the displacement and internal force of tunnel structures in the second stage.The vertical displacement distribution law of the side tunnel structure caused by the excavation of foundation pits was analyzed, and it was verified by using numerical simulation methods.By comparing research results from the two-stage analysis method, the paper identifies some advantages and disadvantages of the current approach, offering valuable insights for future scholars.

Two-Stage Analysis Method
The current method for theoretically analyzing the mechanical response of an existing tunnel structure due to the excavation of a foundation pit typically involves a two-stage analysis.In the first stage, the stress field theory of existing soil is used to calculate the vertical additional stress at the axis of the tunnel, assuming the tunnel does not exist.In the second stage, the tunnel structure is treated as an infinite beam on an elastic foundation, and the vertical additional stress calculated in the first stage is applied to the tunnel to solve for its deformation and internal forces analytically.A diagram illustrating the process of this two-stage analysis method can be seen in Figure 1.
The method for addressing additional stress in the first stage of the two-stage analysis primarily relies on the unloading model of the bottom and side walls of the foundation pit, following established stress solutions.In the second stage, addressing tunnel displacement and internal force is predominantly based on the tunnel-foundation interaction model, with the corresponding mechanical equilibrium differential equation being established for calculation.Various scholars have incorporated multiple engineering factors into the calculation process, enhancing the accuracy of determining additional stress [26], displacement [27], and internal force [28].The method for addressing additional stress in the first stage of the two-stage analysis primarily relies on the unloading model of the bottom and side walls of the foundation pit, following established stress solutions.In the second stage, addressing tunnel displacement and internal force is predominantly based on the tunnel-foundation interaction model, with the corresponding mechanical equilibrium differential equation being established for calculation.Various scholars have incorporated multiple engineering factors into the calculation process, enhancing the accuracy of determining additional stress [26], displacement [27], and internal force [28].

Additional Stress Calculation of the First Stage
The mechanical calculation model illustrating the impact of foundation pit excavation on adjacent tunnels is depicted in Figure 2. In the initial stage, the additional load at the tunnel position due to the unloading of foundation pit excavation is analyzed by considering the unloading effects of the foundation pit bottom and side walls (①~④), utilizing the Mindlin stress solution [29].A point (x1, L0, z0) on the tunnel axis experiences a unit force γddξdη from a point (ξ, η) in the rectangular uniform load at the pit bottom, resulting in vertical additional stress ( d z σ ) in the tunnel axis direction, as per Equation (1).Similarly, under a unit force K0γτdηdτ from a point (η, τ) in the triangular distributed load on the side wall (①), the vertical additional stress ( 1 c z σ ) in the tunnel axis direction is calculated using Equation (2).This process is repeated for side walls (②~④) through the superposition principle to determine the total additional stress on the tunnel induced by the soil unloading from the foundation pit bottom and walls post-excavation.Many scholars have further refined this calculation method by considering various factors to align the additional load calculation more closely with engineering practicality.

Additional Stress Calculation of the First Stage
The mechanical calculation model illustrating the impact of foundation pit excavation on adjacent tunnels is depicted in Figure 2. In the initial stage, the additional load at the tunnel position due to the unloading of foundation pit excavation is analyzed by considering the unloading effects of the foundation pit bottom and side walls ( 1⃝~4 ⃝), utilizing the Mindlin stress solution [29].A point (x 1 , L 0 , z 0 ) on the tunnel axis experiences a unit force γddξdη from a point (ξ, η) in the rectangular uniform load at the pit bottom, resulting in vertical additional stress (σ d z ) in the tunnel axis direction, as per Equation (1).Similarly, under a unit force K0γτdηdτ from a point (η, τ) in the triangular distributed load on the side wall ( 1⃝), the vertical additional stress (σ z ) in the tunnel axis direction is calculated using Equation (2).This process is repeated for side walls ( 2⃝~4 ⃝) through the superposition principle to determine the total additional stress on the tunnel induced by the soil unloading from the foundation pit bottom and walls post-excavation.Many scholars have further refined this calculation method by considering various factors to align the additional load calculation more closely with engineering practicality.

∫∫ ∫∫
In the formula, Ω represents the integral area at the bottom of the pit; Г represents the integral area on the side wall area (①); and ν represents the Poisson's ratio of the soil.The calculation methods of R1, R2, T1, and T2 are shown in Equations ( 3) and (4).

Base and Foundation Pit Side Wall Unloading
Chen et al. [30], Wang et al. [31], and Zhang et al. [32] have focused on the longitudinal additional stress induced at the existing tunnel position due to unloading at the foundation pit bottom.The analytical formula-derived tunnel uplift value is slightly lower than the field-measured value.This suggests that neglecting the unloading effect of the foundation pit's side wall impacts the displacement calculation of the tunnel structure.
Jiang et al. [33] analyzed the impact of unloading the pit wall and bottom soil on a tunnel, developing an analytical formula for the tunnel's longitudinal displacement due to foundation pit excavation.Building on this work, Zhang et al. [34,35] focused on the In the formula, Ω represents the integral area at the bottom of the pit; Γ represents the integral area on the side wall area ( 1⃝); and ν represents the Poisson's ratio of the soil.The calculation methods of R 1 , R 2 , T 1 , and T 2 are shown in Equations ( 3) and (4).

Base and Foundation Pit Side Wall Unloading
Chen et al. [30], Wang et al. [31], and Zhang et al. [32] have focused on the longitudinal additional stress induced at the existing tunnel position due to unloading at the foundation pit bottom.The analytical formula-derived tunnel uplift value is slightly lower than the field-measured value.This suggests that neglecting the unloading effect of the foundation pit's side wall impacts the displacement calculation of the tunnel structure.
Jiang et al. [33] analyzed the impact of unloading the pit wall and bottom soil on a tunnel, developing an analytical formula for the tunnel's longitudinal displacement due to foundation pit excavation.Building on this work, Zhang et al. [34,35] focused on the additional stress from unloading soil around the pit and obtained theoretical solutions for the tunnel displacement and bending moment.While existing studies address unloading effects on the foundation pit, they do not account for the protective role of the retaining structure during excavation, neglecting its impact on additional loads.

Considering the Influence of the Enclosure Structure
Wei et al. [36,37] suggested that the stress released at the bottom of a pit could be influenced by the side friction resistance generated by the retaining structure.They modified the formula for side friction resistance, as presented in Equation (5), where c represents the cohesion of soil, γ is the unit weight of soil, d is the depth of the foundation pit, H is the depth of the retaining structure of the foundation pit, K 0 is the lateral pressure coefficient, and φ is the internal friction angle of soil.Zhou et al. [24] utilized the incremental method [38,39] to analyze the interaction between the supporting structure and the soil.They derived an analytical solution for the vertical displacement of the tunnel, revealing that the impact of foundation pit construction on the tunnel was approximately six times the excavation length.

Considering the Precipitation Effect
Zhou et al. [25], Zhang et al. [40], Ou et al. [41], and Li [42] investigated the impact of engineering dewatering on the deformation of pre-existing tunnels, as illustrated in Figure 3.They utilized the effective stress principle and seepage theory to determine the extra effective stress induced in the soil due to precipitation, as detailed in Equations ( 6) and (7).Their findings indicate that precipitation raises the effective stress of the soil, counteracting the additional stress from excavating the foundation pit.The precipitation associated with excavation serves as a beneficial factor in managing the uplift of the existing underlying tunnel.
effective stress induced in the soil due to precipitation, as detailed in Equations (7).Their findings indicate that precipitation raises the effective stress of the soil, c acting the additional stress from excavating the foundation pit.The precipitation ated with excavation serves as a beneficial factor in managing the uplift of the underlying tunnel.( ) The formula is as follows: h1 represents the groundwater level before precip h2 represents the groundwater level after precipitation, and γw represents the w water.The hydraulic gradient, denoted as i, is positive when the tunnel is locate active area and negative otherwise.

Considering the Double-Hole Effect
Urban subways are commonly designed with double-line parallel configura cases where the clear distance is limited, a 'double-hole effect' may occur.When c ing the total additional stress on tunnel a, it is important to account for the ad stress caused by the contact pressure F between tunnel b and the soil.The addition ( F σ ) generated along the axis of tunnel a can be seen in Figure 4. Zhou et al. [43] and Bu et al. [44,45] have provided an analytical solution for gitudinal additional stress and deformation of a double-hole tunnel's underlying s due to the excavation of a foundation pit, considering the 'double-hole effect'.E The formula is as follows: h 1 represents the groundwater level before precipitation, h 2 represents the groundwater level after precipitation, and γ w represents the weight of water.The hydraulic gradient, denoted as i, is positive when the tunnel is located in the active area and negative otherwise.

Considering the Double-Hole Effect
Urban subways are commonly designed with double-line parallel configurations.In cases where the clear distance is limited, a 'double-hole effect' may occur.When calculating the total additional stress on tunnel a, it is important to account for the additional stress caused by the contact pressure F between tunnel b and the soil.The additional stress (σ F ) generated along the axis of tunnel a can be seen in Figure 4. ( ) The formula is as follows: h1 represents the groundwat h2 represents the groundwater level after precipitation, and water.The hydraulic gradient, denoted as i, is positive when active area and negative otherwise.

Considering the Double-Hole Effect
Urban subways are commonly designed with double-lin cases where the clear distance is limited, a 'double-hole effec ing the total additional stress on tunnel a, it is important t stress caused by the contact pressure F between tunnel b and t ( F σ ) generated along the axis of tunnel a can be seen in Figu Zhou et al. [43] and Bu et al. [44,45] have provided an an gitudinal additional stress and deformation of a double-hole t due to the excavation of a foundation pit, considering the 'd Zhou et al. [43] and Bu et al. [44,45] have provided an analytical solution for the longitudinal additional stress and deformation of a double-hole tunnel's underlying structure due to the excavation of a foundation pit, considering the 'double-hole effect'.Equation ( 8) illustrates the vertical additional stress resulting from the contact pressure between tunnel b (x 2 , y, z 2 ) and the soil on the axis of tunnel a (x 1 , y, z 1 ).The findings indicate that the calculated additional stress of a double-hole tunnel caused by foundation pit excavation is lower compared to that of a single-hole tunnel.This is attributed to the fact that the presence of an adjacent double-hole tunnel effectively increases the overall stiffness of the stratum.
The formula represents the contact pressure between tunnel b and the soil, where F is the contact pressure and L m is the integral range along the y-axis of the tunnel.The longitudinal influence range of the excavation load of the foundation pit is generally between −6L and 6L.H stands for the excavation depth of the foundation pit, with 9) and (10).

Considering the Residual Stress of Soil and Soil Nonlinearity
Zhang [46], Zhang et al. [47], and Zhang et al. [48] discussed soil behavior as viscoelastic and applied a viscoelastic transformation to the Mindlin elastic half-space stress solution.This approach allowed them to derive an analytical solution for the additional stress near a tunnel induced by excavation unloading in viscoelastic soil.
Liu [49] and Xu et al. [50] introduced the residual unloading stress coefficient α in their research.Equation ( 11) defines α 0 as a function of 0.3 and h, where h represents the thickness of the upper soil at the calculation point and h r is calculated as H/(0.612H+ 0.19), with H being the excavation depth of the foundation pit.By considering the residual stress of the soil, Equation ( 12) calculates the unloading stress, with γ representing the weight of the soil.Finally, utilizing the Mindlin solution, the additional load imposed on the tunnel structure due to the excavation of the foundation pit is determined.
Based on the scholarly research conducted by the aforementioned experts, it has been observed that the calculation method for determining the additional stress at the tunnel axis due to unloading at the bottom of the pit is evidently flawed.While taking into account the unloading of the side walls surrounding the foundation pit helps address these issues, practical engineering factors have not been fully considered.During the construction of foundation pits, excavation and support operations often occur simultaneously, necessitating the evaluation of the enclosure structure's impact on reducing additional stress.Factors such as project settlement, the presence of double-hole tunnels, and soil characteristics play significant roles in offsetting additional stress and should not be overlooked in the calculation process.
While the calculation method for the additional stress at the tunnel axis caused by the excavation of the foundation pit is relatively comprehensive, it fails to consider the impact of time-space effects.Factors such as the excavation of multiple foundation pits, varying shapes of foundation pits, and different excavation methods should be thoroughly analyzed for their influence on the calculation of additional stress.

The Second Stage of the Tunnel Displacement and Internal Force Calculation
The current predominant interaction model between the foundation and tunnel structures is the Winkler model [51], which is primarily defined by the foundation bed coefficient, k.However, this model may lead to some discrepancies between the calculated results and actual outcomes [52].To address this issue, some researchers have introduced twoparameter models like the Pasternak model [53] and Vlasov model [54].Additionally, a three-parameter Kerr model was proposed by another scholar [55].The interaction model of foundation and tunnel structures is illustrated in Figure 5, with the corresponding equilibrium differential equation presented in Equation ( 13).varying shapes of foundation pits, and different excavation methods should be thoroughly analyzed for their influence on the calculation of additional stress.

The Second Stage of the Tunnel Displacement and Internal Force Calculation
The current predominant interaction model between the foundation and tunnel structures is the Winkler model [51], which is primarily defined by the foundation bed coefficient, k.However, this model may lead to some discrepancies between the calculated results and actual outcomes [52].To address this issue, some researchers have introduced two-parameter models like the Pasternak model [53] and Vlasov model [54].Additionally, a three-parameter Kerr model was proposed by another scholar [55].The interaction model of foundation and tunnel structures is illustrated in Figure 5, with the corresponding equilibrium differential equation presented in Equation ( 13).
Buildings 2024, 14, 2246 9 of 22 In the formula, q(x) represents the distributed load on the tunnel, while w(x) represents the deflection of the tunnel.The variables k and c refer to the foundation bed coefficient, and G and 2t denote the shear stiffness of the tunnel.
In the second stage, the vertical displacement and internal force of the tunnel are calculated using an equilibrium differential equation based on the foundation-tunnel interaction model.By considering the boundary conditions, the vertical displacement of the tunnel due to additional stress is determined, and the internal force of the tunnel is calculated using the differential relationship.
The Winkler model is utilized in this study.Initially, the equilibrium differential equation (Equation ( 14)) is derived.Subsequently, based on the general solution and the boundary conditions, it is determined that at any point ξ along the tunnel structure, the concentrated load is q(ξ)dξ.This load results in the displacement of the point on the tunnel axis, as depicted in Equation (15).
The calculation method of λ is illustrated in Equation ( 16).
Similar to the Pasternak model, the equilibrium differential equation is first established, as depicted in Equation ( 17).Subsequently, the numerical solution can be derived using the finite difference method.The tunnel of length L is divided into n units, with each unit having a length of l = L/n.In this case, the tunnel is partitioned into n + 5 units (including four virtual units at both ends), as depicted in Figure 6.The differential results are then obtained, as shown in Equation (18).These results are further transformed into a stiffness matrix, as shown in Equation (19), to solve for the displacement and internal force of the tunnel based on the differential relationship.
In the formula, q(x) represents the distributed load on the tunnel, while w(x) repre sents the deflection of the tunnel.The variables k and c refer to the foundation bed coeffi cient, and G and 2t denote the shear stiffness of the tunnel.
In the second stage, the vertical displacement and internal force of the tunnel ar calculated using an equilibrium differential equation based on the foundation-tunnel in teraction model.By considering the boundary conditions, the vertical displacement of th tunnel due to additional stress is determined, and the internal force of the tunnel is calcu lated using the differential relationship.
The Winkler model is utilized in this study.Initially, the equilibrium differentia equation (Equation ( 14)) is derived.Subsequently, based on the general solution and th boundary conditions, it is determined that at any point ξ along the tunnel structure, th concentrated load is q(ξ)dξ.This load results in the displacement of the point on the tunne axis, as depicted in Equation (15).
The calculation method of λ is illustrated in Equation ( 16).17).Subsequently, the numerical solution can be derive using the finite difference method.The tunnel of length L is divided into n units, with eac unit having a length of l = L/n.In this case, the tunnel is partitioned into n + 5 units (in cluding four virtual units at both ends), as depicted in Figure 6.The differential result are then obtained, as shown in Equation (18).These results are further transformed into stiffness matrix, as shown in Equation (19), to solve for the displacement and internal forc of the tunnel based on the differential relationship.
In the formula, [K1], [K2], and [G0] represent the tunnel bending stiffness matrix, th foundation elastic stiffness matrix, and the foundation shear stiffness matrix, respectively In the formula, [K 1 ], [K 2 ], and [G 0 ] represent the tunnel bending stiffness matrix, the foundation elastic stiffness matrix, and the foundation shear stiffness matrix, respectively.{w} and {q} denote the vertical displacement and additional stress vector of the tunnel, respectively.
Similar to the Kerr model, this study first establishes the mechanical equilibrium equation (Equation ( 20)).Subsequently, the solutions for ψ 1 , ψ 2 , and ψ 3 are obtained using the general solution presented in Equation (21).By combining Equation (20) and Equation (21), the solution for the equation in Equation ( 22) is derived.Furthermore, the general solution for the shear layer displacement is determined based on its characteristic equation, as shown in Equation (23).The subsequent solution process follows a similar approach to the Winkler model, which will not be further elaborated in this context.

Considering the Shear Action
Current theoretical analysis methods often model tunnels as Euler-Bernoulli or elastic straight beams, focusing solely on bending effects and neglecting shear effects.While this approach is suitable for underground tunnels, for subway tunnels with discontinuous integral structures, it is more appropriate to analyze beam deformation under shear using the theory introduced by Timoshenko et al. [56].Figure 7 illustrates the Euler-Bernoulli beam and Timoshenko beam models.
tion, as shown in Equation (23).The subsequent solution process follows a similar approach to the Winkler model, which will not be further elaborated in this context.

Considering the Shear Action
Current theoretical analysis methods often model tunnels as Euler-Bernoulli or elastic straight beams, focusing solely on bending effects and neglecting shear effects.While this approach is suitable for underground tunnels, for subway tunnels with discontinuous integral structures, it is more appropriate to analyze beam deformation under shear using the theory introduced by Timoshenko et al. [56].Zong [57] and Liang et al. [58,59] both presented analytical solutions for the mechanical response of the Euler-Bernoulli beam and Timoshenko beam on the Kerr foundation under additional stress.The findings indicate that the displacement and internal force law derived from the Timoshenko beam model is more appropriate.Moreover, as the tunnel shear stiffness G increases, the maximum longitudinal deformation (w max ) gradually decreases.It is observed that when G reaches 9 × 10 9 N, the longitudinal deformation of the tunnel remains relatively constant, as illustrated in Figure 8. Zong [57] and Liang et al. [58,59] both presented analytical solutions fo ical response of the Euler-Bernoulli beam and Timoshenko beam on the K under additional stress.The findings indicate that the displacement and int derived from the Timoshenko beam model is more appropriate.Moreover shear stiffness G increases, the maximum longitudinal deformation (wmax) creases.It is observed that when G reaches 9 × 10 9 N, the longitudinal defo tunnel remains relatively constant, as illustrated in Figure 8. Feng et al. [60,61] simplified the tunnel as an infinite Timoshenko beam its behavior on a three-parameter Kerr foundation model to obtain an anal for longitudinal deformation.Their findings indicate that tunnel deforma with a higher shear stiffness, while the bending moment and shear force i tionally, an increasing foundation modulus and tunnel depth lead to a r deformation, bending moment, and shear force.
Liu et al. [62,63] simplified the shield tunnel as a Timoshenko beam on dation to analyze the uplift deformation of the underlying tunnel due to excavation.By optimizing the construction parameters of the new fou Guimiao Road in Shenzhen City based on this model, the study found tha results aligned well with field measurements, validating the accuracy and e the approach.

Considering the Effect of Lateral Soil
In practical engineering, the additional load from excavating the upp pit not only impacts the tunnel and transfers to the lower foundation but a lateral stratum.The deformation of the lateral soil differs from that of the t pressure on the tunnel.It is hypothesized that the lateral soil forces T1 and T are transferred through the shear layer and exerted on both sides of the illustrated in Equation ( 24) and Figure 9.

( ) T T tkw x = =
In the formula, 2t and k represent the shear stiffness and the foundat cient, respectively.Feng et al. [60,61] simplified the tunnel as an infinite Timoshenko beam and analyzed its behavior on a three-parameter Kerr foundation model to obtain an analytical solution for longitudinal deformation.Their findings indicate that tunnel deformation decreases with a higher shear stiffness, while the bending moment and shear force increase.Additionally, an increasing foundation modulus and tunnel depth lead to a reduced tunnel deformation, bending moment, and shear force.
Liu et al. [62,63] simplified the shield tunnel as a Timoshenko beam on a Vlazov foundation to analyze the uplift deformation of the underlying tunnel due to foundation pit excavation.By optimizing the construction parameters of the new foundation pit at Guimiao Road in Shenzhen City based on this model, the study found that the predicted results aligned well with field measurements, validating the accuracy and effectiveness of the approach.

Considering the Effect of Lateral Soil
In practical engineering, the additional load from excavating the upper foundation pit not only impacts the tunnel and transfers to the lower foundation but also affects the lateral stratum.The deformation of the lateral soil differs from that of the tunnel, causing pressure on the tunnel.It is hypothesized that the lateral soil forces T 1 and T 2 on the tunnel are transferred through the shear layer and exerted on both sides of the tunnel [64], as illustrated in Equation ( 24) and Figure 9.
In the formula, 2t and k represent the shear stiffness and the foundation bed coefficient, respectively.
Huang et al. [64] derived an analytical solution for the longitudinal deformation of a tunnel beneath lateral soil during the unloading conditions of foundation pit excavation.They recommended reinforcing the soil around the tunnel's lateral sides to minimize the impact of upper unloading on the tunnel.Tang [65] simplified the tunnel as an infinite Euler-Bernoulli beam on a Vlasov foundation model, accounting for the lateral soil effects and constraints at both tunnel ends.Their analytical solution for vertical tunnel deformation closely matched field monitoring data.Huang et al. [64] derived an analytical solution for the longitudinal deformat tunnel beneath lateral soil during the unloading conditions of foundation pit exc They recommended reinforcing the soil around the tunnel's lateral sides to minim impact of upper unloading on the tunnel.Tang [65] simplified the tunnel as an Euler-Bernoulli beam on a Vlasov foundation model, accounting for the lateral soi and constraints at both tunnel ends.Their analytical solution for vertical tunne mation closely matched field monitoring data.

Considering the Buried Depth Effect of the Tunnel
The foundation bed coefficient is calculated using the formula proposed by V based on long beam load tests on the surface.Attewell et al. [67] later suggeste twice the Vesić subgrade coefficient to estimate the subgrade coefficient of a tunn infinite depth, which could lead to an overestimation of soil stiffness.To addr many scholars have introduced the coefficient η proposed by Yu et al. [68], whic into account the influence of the tunnel buried depth.This coefficient considers th of the tunnel buried depth, as shown in Equations ( 25) and ( 26) (where h is the depth of the tunnel, D is the outer diameter of the tunnel, Es is the elastic modulu foundation soil, ν is the Poisson's ratio of the foundation soil, and EI is the equ bending stiffness of the tunnel).
In light of the impact of tunnel depth, Liang et al. [69] and Xu et al. [70] intro revised formula for the foundation bed coefficient and derived a mechanical resp lution for tunnels under foundation pit unloading, taking into account the influ tunnel depth.Building upon the foundation bed coefficient put forth by Cheng et Ying et al. [72], and Huang et al. [73], which incorporates the buried depth effec

Considering the Buried Depth Effect of the Tunnel
The foundation bed coefficient is calculated using the formula proposed by Vesić [66] based on long beam load tests on the surface.Attewell et al. [67] later suggested using twice the Vesić subgrade coefficient to estimate the subgrade coefficient of a tunnel with infinite depth, which could lead to an overestimation of soil stiffness.To address this, many scholars have introduced the coefficient η proposed by Yu et al. [68], which takes into account the influence of the tunnel buried depth.This coefficient considers the effect of the tunnel buried depth, as shown in Equations ( 25) and ( 26) (where h is the buried depth of the tunnel, D is the outer diameter of the tunnel, E s is the elastic modulus of the foundation soil, ν is the Poisson's ratio of the foundation soil, and EI is the equivalent bending stiffness of the tunnel).
In light of the impact of tunnel depth, Liang et al. [69] and Xu et al. [70] introduced a revised formula for the foundation bed coefficient and derived a mechanical response solution for tunnels under foundation pit unloading, taking into account the influence of tunnel depth.Building upon the foundation bed coefficient put forth by Cheng et al. [71], Ying et al. [72], and Huang et al. [73], which incorporates the buried depth effect of the tunnel, a calculation method for assessing the response of adjacent existing tunnels due to foundation pit excavation is proposed.The findings indicate that the maximum vertical displacement of the tunnel is significantly influenced by its depth, showing an approximately linear relationship.

Considering the Non-Uniformity of Soil
Jiang et al. [74] analyzed the variability of foundation soil by dividing it into n layers based on the Burmister elastic layered theory [75][76][77].Figure 10 illustrates the action diagram of the axisymmetric vertical load and presents the derived analytical solution for the longitudinal deformation of the tunnel in a layered foundation.The analytical formula accounts for variations in the elastic modulus and Poisson's ratio with soil stratification while also addressing the issue of excessive stress diffusion in a homogeneous foundation.
based on the Burmister elastic layered theory [75][76][77].Figure 10 illustrates the ac gram of the axisymmetric vertical load and presents the derived analytical solutio longitudinal deformation of the tunnel in a layered foundation.The analytical accounts for variations in the elastic modulus and Poisson's ratio with soil strat while also addressing the issue of excessive stress diffusion in a homogeneous fou

Considering the Effect of the Segment Joint and the Constraint Effect of the Tunnel E
Previous methods have treated the tunnel as a complete ring without joints assumed a fixed stiffness for the joints.In contrast, shield tunnels are constructed ing segments with bolts.This construction method assumes the shield tunnel is a uous structure but may not effectively account for the impact of joints on the tunn gitudinal deformation.
Zhou et al.
[78], Wei et al. [79], and Chen et al. [80] investigated the 'shear dis effect' of shield tunnels, as discussed by Shen et al. [81].Additionally, Wei et al.Zhang [85], and Yu [86] proposed that the longitudinal deformation of the tunne from the combined effects of shear dislocation and rigid body rotation, as indic Wang [87].They derived formulas for calculating longitudinal deformation, the tion amount between rings, the rotation angle between rings, and the shear force rings induced by foundation pit excavation.The study findings reveal that the dislocation deformation at the maximum vertical and horizontal displacement p the tunnel.The maximum dislocation between rings occurs at the inflection poi Gaussian curve, as illustrated in Figure 11.

Considering the Effect of the Segment Joint and the Constraint Effect of the Tunnel End
Previous methods have treated the tunnel as a complete ring without joints or have assumed a fixed stiffness for the joints.In contrast, shield tunnels are constructed by splicing segments with bolts.This construction method assumes the shield tunnel is a continuous structure but may not effectively account for the impact of joints on the tunnel's longitudinal deformation.
Zhou et al.
[78], Wei et al. [79], and Chen et al. [80] investigated the 'shear dislocation effect' of shield tunnels, as discussed by Shen et al. [81].Additionally, Wei et al. [82][83][84], Zhang [85], and Yu [86] proposed that the longitudinal deformation of the tunnel results from the combined effects of shear dislocation and rigid body rotation, as indicated by Wang [87].They derived formulas for calculating longitudinal deformation, the dislocation amount between rings, the rotation angle between rings, and the shear force between rings induced by foundation pit excavation.The study findings reveal that there is no dislocation deformation at the maximum vertical and horizontal displacement points of the tunnel.The maximum dislocation between rings occurs at the inflection point of the Gaussian curve, as illustrated in Figure 11.
accounts for variations in the elastic modulus and Poisson's ratio with soil strat while also addressing the issue of excessive stress diffusion in a homogeneous fou

Considering the Effect of the Segment Joint and the Constraint Effect of the Tunnel E
Previous methods have treated the tunnel as a complete ring without joints assumed a fixed stiffness for the joints.In contrast, shield tunnels are constructed ing segments with bolts.This construction method assumes the shield tunnel is a uous structure but may not effectively account for the impact of joints on the tunn gitudinal deformation.
Zhou et al.
[78], Wei et al. [79], and Chen et al. [80] investigated the 'shear dis effect' of shield tunnels, as discussed by Shen et al. [81].Additionally, Wei et al.Zhang [85], and Yu [86] proposed that the longitudinal deformation of the tunne from the combined effects of shear dislocation and rigid body rotation, as indic Wang [87].They derived formulas for calculating longitudinal deformation, the tion amount between rings, the rotation angle between rings, and the shear force rings induced by foundation pit excavation.The study findings reveal that the dislocation deformation at the maximum vertical and horizontal displacement p the tunnel.The maximum dislocation between rings occurs at the inflection poin Gaussian curve, as illustrated in Figure 11.Wang et al. [88] examined the impact of the station (working shaft) constraint on the end of the shield tunnel.The connection node was modeled as a rotating spring with a rotational stiffness of K θ and a vertical link, as depicted in Figure 12.The analytical solution for the longitudinal response of the adjacent shield tunnel due to the excavation of the foundation pit was derived.This study offers insights for addressing similar challenges within intricate boundary conditions.
Buildings 2024, 14, x FOR PEER REVIEW Wang et al. [88] examined the impact of the station (working shaft) constrain end of the shield tunnel.The connection node was modeled as a rotating spring rotational stiffness of Kθ and a vertical link, as depicted in Figure 12.The analyti tion for the longitudinal response of the adjacent shield tunnel due to the excav the foundation pit was derived.This study offers insights for addressing simil lenges within intricate boundary conditions.

Considering the Nonlinearity of the Segment Joint and Foundation
Foundation models such as Winkler, Pasternak, Vlasov, or Kerr all assume foundation reaction increases linearly with deformation.However, foundation so mation exhibits nonlinear characteristics.As depicted in Figure 13, the foundati tion force does not increase infinitely with soil displacement but rather approache cific critical value.Hence, it is essential to account for the nonlinear deformation teristics of the foundation when calculating displacements in tunnels under fou pit unloading.Zhang et al. [89] investigated the nonlinear variation characteristics of the stiffness of the lining joint to solve the analytical solution of the transverse intern and deformation of tunnels during foundation pit excavation.Kang et al. [90] an et al. [91] developed an analytical solution for the longitudinal deformation of shi nels during the excavation and unloading of foundation pits, based on the nonlin ternak foundation model introduced by Liang [92] (refer to Equation (27), where resents the ratio of the ultimate foundation reaction force qu to the required fou soil deformation).The findings suggest that the calculation method utilizing the n Pasternak foundation model is more dependable.

Considering the Nonlinearity of the Segment Joint and Foundation
Foundation models such as Winkler, Pasternak, Vlasov, or Kerr all assume that the foundation reaction increases linearly with deformation.However, foundation soil deformation exhibits nonlinear characteristics.As depicted in Figure 13, the foundation reaction force does not increase infinitely with soil displacement but rather approaches a specific critical value.Hence, it is essential to account for the nonlinear deformation characteristics of the foundation when calculating displacements in tunnels under foundation pit unloading.Wang et al. [88] examined the impact of the station (working shaft) constrain end of the shield tunnel.The connection node was modeled as a rotating sprin rotational stiffness of Kθ and a vertical link, as depicted in Figure 12.The analyti tion for the longitudinal response of the adjacent shield tunnel due to the excav the foundation pit was derived.This study offers insights for addressing simi lenges within intricate boundary conditions.

Considering the Nonlinearity of the Segment Joint and Foundation
Foundation models such as Winkler, Pasternak, Vlasov, or Kerr all assume foundation reaction increases linearly with deformation.However, foundation so mation exhibits nonlinear characteristics.As depicted in Figure 13, the foundat tion force does not increase infinitely with soil displacement but rather approach cific critical value.Hence, it is essential to account for the nonlinear deformation teristics of the foundation when calculating displacements in tunnels under fou pit unloading.Zhang et al. [89] investigated the nonlinear variation characteristics of the stiffness of the lining joint to solve the analytical solution of the transverse inter and deformation of tunnels during foundation pit excavation.Kang et al. [90] an et al. [91] developed an analytical solution for the longitudinal deformation of sh nels during the excavation and unloading of foundation pits, based on the nonlin ternak foundation model introduced by Liang [92] (refer to Equation (27), wher resents the ratio of the ultimate foundation reaction force qu to the required fou soil deformation).The findings suggest that the calculation method utilizing the n Pasternak foundation model is more dependable.
Based on the research conducted by the aforementioned scholars, it has bee mined that the foundation and tunnel model, which takes into account shear ac Zhang et al. [89] investigated the nonlinear variation characteristics of the bending stiffness of the lining joint to solve the analytical solution of the transverse internal force and deformation of tunnels during foundation pit excavation.Kang et al. [90] and Wang et al. [91] developed an analytical solution for the longitudinal deformation of shield tunnels during the excavation and unloading of foundation pits, based on the nonlinear Pasternak foundation model introduced by Liang [92] (refer to Equation (27), where k u represents the ratio of the ultimate foundation reaction force q u to the required foundation soil deformation).The findings suggest that the calculation method utilizing the nonlinear Pasternak foundation model is more dependable.
Based on the research conducted by the aforementioned scholars, it has been determined that the foundation and tunnel model, which takes into account shear action and lateral soil forces, offers a more accurate representation of the actual stress distribution.Additionally, incorporating a calculation model for segment joints enables a more precise analysis of additional displacement and internal forces in shield tunnels.By considering factors such as the tunnel depth, soil non-uniformity, and foundation non-linearity, potential drawbacks of excessive foundation reactions can be effectively mitigated.Moreover, the integration of rotating springs and vertical links presents a novel approach for calculating tunnel displacements under complex boundary conditions.
While the scholars have successfully addressed the issue of displacement and internal force at the tunnel axis resulting from foundation pit excavation, their calculation outcome only provides the overall displacement of the tunnel along the axis direction.It does not offer the displacement of specific positions within the tunnel, such as the vault or inverted arch.This limitation stems from simplifying the tunnel as a beam, and the further analysis of these specific areas could be explored in future research.

Engineering Case Verification
A new research project focusing on the transfer channel of the Beijing subway rail transit system is being conducted.The foundation pit for excavation measures 20 m in length, 20 m in width, and 30 m in depth.To support the pit, a combination of a diaphragm wall and steel support is being utilized.The No. 2 subway line is situated 6 m deep and 12 m away from the foundation pit. Figure 14 illustrates the spatial relationship between the foundation pit and the existing tunnel structure, while Table 1 provides the physical and mechanical parameters of the stratum and structure.
Buildings 2024, 14, x FOR PEER REVIEW lateral soil forces, offers a more accurate representation of the actual stress distr Additionally, incorporating a calculation model for segment joints enables a more analysis of additional displacement and internal forces in shield tunnels.By cons factors such as the tunnel depth, soil non-uniformity, and foundation non-linearity tial drawbacks of excessive foundation reactions can be effectively mitigated.Mo the integration of rotating springs and vertical links presents a novel approach fo lating tunnel displacements under complex boundary conditions.
While the scholars have successfully addressed the issue of displacement an nal force at the tunnel axis resulting from foundation pit excavation, their calculat come only provides the overall displacement of the tunnel along the axis direction not offer the displacement of specific positions within the tunnel, such as the vau verted arch.This limitation stems from simplifying the tunnel as a beam, and the analysis of these specific areas could be explored in future research.

Engineering Case Verification
A new research project focusing on the transfer channel of the Beijing subw transit system is being conducted.The foundation pit for excavation measures length, 20 m in width, and 30 m in depth.To support the pit, a combination of a dia wall and steel support is being utilized.The No. 2 subway line is situated 6 m d 12 m away from the foundation pit. Figure 14 illustrates the spatial relationship b the foundation pit and the existing tunnel structure, while Table 1 provides the p and mechanical parameters of the stratum and structure.3D 6.0 finite difference software is utilized to develop a three-dimensio merical model (Figure 15).The stratum is modeled using the Mohr-Coulomb cons  FLAC 3D 6.0 finite difference software is utilized to develop a three-dimensional numerical model (Figure 15).The stratum is modeled using the Mohr-Coulomb constitutive model, while the tunnel structure, diaphragm wall, and steel support are modeled using the elastic constitutive model.The top surface of the model is a free surface, and its front, rear, left, right, and bottom surfaces are all constrained to the displacement.The physical and mechanical parameters can be found in Table 1.The construction process of the numerical model consists of two main parts.The first part involves constructing the existing tunnel section and diaphragm wall.Upon completion, displacements are cleared, assuming soil consolidation settlement is complete, leaving only the stress field.This state is considered the initial state for foundation pit construction  In theoretical calculations, it is assumed that the soil is homogeneous within an elastic half space, and the tunnel structure is an infinitely long homogeneous elastic body, without considering the influence of the tunnel's existence on the calculation of additional stress in the soil.Engineering factors such as the base unloading, sidewall unloading, and burial depth are considered for the corresponding theoretical analysis.Moreover, the reduction effect caused by the unloading of the foundation and side walls of the foundation pit, as well as the support structure (using a reduction coefficient of 0.75), is considered, along with other relevant parameters detailed in Table 1, assuming a tunnel influence range of 30 m. Comparing the calculated results with the numerical simulation results depicted in Figure 16, it is evident that both patterns are largely consistent.The maximum vertical displacement occurs at the center of the tunnel and gradually decreases as the distance from the center increases.The maximum error observed is only 5.9%, thus validating the rationality of the theoretical calculation.In theoretical calculations, it is assumed that the soil is homogeneous within an elastic half space, and the tunnel structure is an infinitely long homogeneous elastic body, without considering the influence of the tunnel's existence on the calculation of additional stress in the soil.Engineering factors such as the base unloading, sidewall unloading, and burial depth are considered for the corresponding theoretical analysis.Moreover, the reduction effect caused by the unloading of the foundation and side walls of the foundation pit, as well as the support structure (using a reduction coefficient of 0.75), is considered, along with other relevant parameters detailed in Table 1, assuming a tunnel influence range of 30 m. Comparing the calculated results with the numerical simulation results depicted in Figure 16, it is evident that both patterns are largely consistent.The maximum vertical displacement occurs at the center of the tunnel and gradually decreases as the distance from the center increases.The maximum error observed is only 5.9%, thus validating the rationality of the theoretical calculation.To further investigate the influence of the excavation of the foundation pit longitudinal displacement of adjacent tunnels at different distances between the and foundation pit, numerical simulation methods were used.The depth to 6 m trolled, with distances of 3, 6, 12, 18, and 30 m, respectively.The remaining parame shown in Table 1, and the boundary conditions and construction process of the mo as described earlier.Figure 17 shows the variation curve of the vertical displace the tunnel structure caused by the excavation of the foundation pit at different di It can be observed that as the distance increases, the vertical displacement of the structure tends to be smaller on both sides and larger in the middle.Moreover, th mum vertical displacement of the tunnel structure shows a pattern of first increas then decreasing with the increase in the distance.There are too many parameters in the paper.In order to express them more c all parameters have been summarized in Table 2.

Parameters
Definition Parameters Definition The additional stress on the tunnel caused by To further investigate the influence of the excavation of the foundation pit on the longitudinal displacement of adjacent tunnels at different distances between the tunnel and foundation pit, numerical simulation methods were used.The depth to 6 m is controlled, with distances of 3, 6, 12, 18, and 30 m, respectively.The remaining parameters are shown in Table 1, and the boundary conditions and construction process of the model are as described earlier.Figure 17 shows the variation curve of the vertical displacement of the tunnel structure caused by the excavation of the foundation pit at different distances.It can be observed that as the distance increases, the vertical displacement of the tunnel structure tends to be smaller on both sides and larger in the middle.Moreover, the maximum vertical displacement of the tunnel structure shows a pattern of first increasing and then decreasing with the increase in the distance.To further investigate the influence of the excavation of the foundation pit longitudinal displacement of adjacent tunnels at different distances between the and foundation pit, numerical simulation methods were used.The depth to 6 m trolled, with distances of 3, 6, 12, 18, and 30 m, respectively.The remaining parame shown in Table 1, and the boundary conditions and construction process of the m as described earlier.Figure 17 shows the variation curve of the vertical displace the tunnel structure caused by the excavation of the foundation pit at different di It can be observed that as the distance increases, the vertical displacement of the structure tends to be smaller on both sides and larger in the middle.Moreover, th mum vertical displacement of the tunnel structure shows a pattern of first increas then decreasing with the increase in the distance.There are too many parameters in the paper.In order to express them more c all parameters have been summarized in Table 2.There are too many parameters in the paper.In order to express them more clearly, all parameters have been summarized in Table 2.The ratio of the ultimate foundation reaction q u to the required deformation of the foundation soil

Conclusions
The theoretical analysis method stands out among other research methods for its ability to uncover the mechanical mechanism of foundation pits near existing tunnels.This paper consolidates the research findings of numerous scholars worldwide on this topic, highlighting the strengths and weaknesses of the current theoretical analysis.It offers valuable insights and guidance for future scholars in this field.
(1) In the first stage, the calculation of additional stress takes into account factors such as the unloading model of the foundation pit bottom and side wall, the influence of the retaining structure, the dewatering effect, the double hole effect, the soil material, and residual stress.This approach aligns more closely with engineering practice and yields more accurate results.
(2) In the second stage, the displacement and internal force of the tunnel are calculated using the foundation-tunnel interaction model.The theoretical deduction takes into account various factors such as the shear effect, lateral soil effect, tunnel burial depth effect, soil non-uniformity, tunnel end constraint effect, segment joint dislocation deformation, and nonlinear influence of the tunnel and foundation.
(3) The existing tunnel will have a blocking effect on the surrounding soil.The blocking effect can reduce the soil pressure on the underground continuous wall after the excavation of the foundation pit, which may have an impact on the displacement and internal force of the existing tunnel structure.Further research can be conducted to investigate the impact of this blocking effect on the deformation and internal forces of the tunnel.
(4) Currently, the tunnel structure section is predominantly assumed to be circular.However, in reality, tunnel sections are often non-circular.Further theoretical research on the mechanical response of non-circular section tunnels during foundation pit excavation can be explored.It is recommended to combine the principle of conformal transformation in complex functions to process non-circular cross-sections in order to make the calculation results more in line with engineering practice.
(5) While numerous engineering factors are taken into account during theoretical analysis, the resulting analytical solutions often prove to be overly complex.These solutions can be simplified by expressing them as reduction or amplification factors, which are then combined with field-measured data.This simplified approach can be utilized with the basic formula to streamline calculations.(6) The validity of theoretical calculations was confirmed through engineering examples and numerical simulations, with a maximum error of only 5.9%, demonstrating good consistency between theoretical calculations and numerical simulations.The two-stage analysis method can be used to analyze the mechanical response of side tunnels caused by the excavation of foundation pits.

Figure 1 .
Figure 1.The schematic diagram of the solution process of the two-stage analysis method.

Figure 1 .
Figure 1.The schematic diagram of the solution process of the two-stage analysis method.

Figure 2 .
Figure 2. Mechanical calculation model of the influence of foundation pit excavation on the adjacent tunnel.

Figure 2 .
Figure 2. Mechanical calculation model of the influence of foundation pit excavation on the adjacent tunnel.

Figure 3 .
Figure 3.The influence diagram of engineering dewatering on the additional stress of soil.

Figure 4 .
Figure 4. Calculation diagram of additional stress considering the double-hole effect.

Figure 3 .
Figure 3.The influence diagram of engineering dewatering on the additional stress of soil.

Figure 3 .
Figure 3.The influence diagram of engineering dewatering on the a

Figure 4 .
Figure 4. Calculation diagram of additional stress considering the d

Figure 4 .
Figure 4. Calculation diagram of additional stress considering the double-hole effect.
(a) Winkler foundation model (b) Pasternak foundation model (c) Vlasov foundation model (d) Kerr foundation model

Figure 5 .
Figure 5. Interaction model between the foundation and tunnel.The advantage of the Winkler foundation model is that it has a simple form and is mathematically easy to handle, but it does not consider the continuity and shear action of the soil.The Pasternak and Vlasov foundation models overcome the shortcomings of the Winkler foundation model by introducing shear stiffness parameters and can more accurately simulate soil structure interactions, especially in cases of high soil shear stiffness.The Kerr foundation model adds an upper layer spring on the basis of the Pasternak model to simulate the vertical stress distribution of the soil, thereby more comprehensively reflecting the mechanical behavior of the soil.However, this also means that the calculation of the Kerr model may be more complex.

Figure 5 .
Figure 5. Interaction model between the foundation and tunnel.The advantage of the Winkler foundation model is that it has a simple form and is mathematically easy to handle, but it does not consider the continuity and shear action of the soil.The Pasternak and Vlasov foundation models overcome the shortcomings of the Winkler foundation model by introducing shear stiffness parameters and can more accurately simulate soil structure interactions, especially in cases of high soil shear stiffness.The Kerr foundation model adds an upper layer spring on the basis of the Pasternak model to simulate the vertical stress distribution of the soil, thereby more comprehensively reflecting the mechanical behavior of the soil.However, this also means that the calculation of the Kerr model may be more complex.
to the Pasternak model, the equilibrium differential equation is first estab lished, as depicted in Equation (

Figure 8 .
Figure 8. Vertical displacement of a tunnel under different shear stiffnesses.

Figure 8 .
Figure 8. Vertical displacement of a tunnel under different shear stiffnesses.

Figure 10 .
Figure 10.A schematic diagram of an axisymmetric load applied to the surface of a multila dation.

Figure 10 .
Figure 10.A schematic diagram of an axisymmetric load applied to the surface of a multilayer foundation.

Figure 10 .
Figure 10.A schematic diagram of an axisymmetric load applied to the surface of a multila dation.

Figure 11 .
Figure 11.Dislocation platform and corner of the shield tunnel.

Figure 11 .
Figure 11.Dislocation platform and corner of the shield tunnel.

Figure 12 .
Figure 12.Mechanical calculation model considering the end constraint effect.

Figure 13 .
Figure 13.Relationship between foundation reaction and displacement.

Figure 12 .
Figure 12.Mechanical calculation model considering the end constraint effect.

Figure 11 .
Figure 11.Dislocation platform and corner of the shield tunnel.

Figure 12 .
Figure 12.Mechanical calculation model considering the end constraint effect.

Figure 13 .
Figure 13.Relationship between foundation reaction and displacement.

Figure 13 .
Figure 13.Relationship between foundation reaction and displacement.

Figure 14 .
Figure 14.Location relationship between the foundation pit and tunnel (m).

Figure 14 .
Figure 14.Location relationship between the foundation pit and tunnel (m).
. The second part involves excavating the foundation pit in 12 layers, each 2.5 m thick, with a total of three steel supports installed at depths of −3 m, −12 m, and −21 m below the surface.PEER REVIEW 16 of 22 rear, left, right, and bottom surfaces are all constrained to the displacement.The physical and mechanical parameters can be found in Table 1.The construction process of the numerical model consists of two main parts.The first part involves constructing the existing tunnel section and diaphragm wall.Upon completion, displacements are cleared, assuming soil consolidation settlement is complete, leaving only the stress field.This state is considered the initial state for foundation pit construction.The second part involves excavating the foundation pit in 12 layers, each 2.5 m thick, with a total of three steel supports installed at depths of −3 m, −12 m, and −21 m below the surface.(a) Overall computational model (b) Existing structures and new structures

Buildings 2024 ,Figure 16 .
Figure 16.Comparative analysis of tunnel displacement caused by the excavation of the fo tion pit.

Figure 17 .
Figure 17.Vertical displacement of the tunnel structure caused by the excavation of the fou pit.

Figure 16 .
Figure 16.Comparative analysis of tunnel displacement caused by the excavation of the foundation pit.

Buildings 2024 ,Figure 16 .
Figure 16.Comparative analysis of tunnel displacement caused by the excavation of the fo tion pit.

Figure 17 .
Figure 17.Vertical displacement of the tunnel structure caused by the excavation of the fou pit.

σ
The additional stress on the tunnel caused by the unloading of the soil at the bottom of the foundation pit after excavation h1 Water level before precipitat Additional stress on the tunnel caused by

Figure 17 .
Figure 17.Vertical displacement of the tunnel structure caused by the excavation of the foundation pit.

Table 2 .z 1 Water level before precipitation σ c 1 zz 2 R 2 = 1 T 1 = (x 1 − L 2 ) 2 + 2 T 2 = (x 1 − L 2 ) 2 +
Parameters Summary.The additional stress on the tunnel caused by the unloading of the soil at the bottom of the foundation pit after excavation h Additional stress on the tunnel caused by unloading the soil on one side of the pit wall after excavation of the foundation pit h 2 Water level after precipitation σ z The total additional stress on the tunnel caused by the unloading of the soil at the bottom and walls of the foundation pit after excavation γ w Volume weight of water x 1 The horizontal coordinate of a point on the tunnel axis i Water surface slope or gradient L 0 The horizontal distance between the center of the foundation pit and the center of the tunnel σ w Additional effective stress in soil caused by precipitation z 0 Depth of tunnel b Contact pressure between tunnel b and soil (ξ, η) A point in the uniformly distributed rectangular load at the bottom of the pit L m Integration range along the b-axis of the tunnel (η, τ) A point in the triangular distribution load on the sidewall α Residual unloading stress coefficient γ Volume weight of soil h When calculating residual stress, calculate the thickness of the soil above the point K 0 Coefficient of lateral pressure q(x) Distributed loads on tunnels Ω Integral region on the bottom area of the pit w(x) The deflection of the tunnel Γ Integral region on the sidewall area k, c The foundation bed coefficient νPoisson's ratio of soil G, 2t Shear stiffness of tunnelsR 1 R 1 = (x 1 − ξ) 2 + (L 0 − η) 2 + (z 0 − d) 2 DThe outer diameter of the tunnelR (x 1 − ξ) 2 + (L 0 − η) 2 + (z 0 + d) 2 EIBending stiffness of tunnelsT (L 0 − η) 2 + (z 0 − τ) 2 T 1 ,T 2 Lateral soil force on tunnel T (L 0 − η) 2 + (z 0 + τ) 2 δ Coefficient that can consider the influence of tunnel depth L The length of the foundation pit E s The elastic modulus of soil B The width of the foundation pit ∆h i The thickness of the i-th layer of soil d The depth of the foundation pit K θ Rotational stiffness c Cohesive of soil q u Ultimate foundation reaction force φ Angle of internal friction of soil k u

Table 1 .
Physical and mechanical parameters of the strata and structure.

Table 1 .
Physical and mechanical parameters of the strata and structure.