Effects of Water Inrush from Tunnel Excavation Face on the Deformation and Mechanical Performance of Shield Tunnel Segment Joints

Water inrush from the excavation face often occurs in the current shield construction of metro tunnels. In this study, the discontinuity of shield tunnel lining and the interaction between the tunnel segments, the grouting layer, and the surrounding rock are considered. Based on the 3D nonlinear contact theory, a hybrid model of the shield tunnel is constructed. Considering the 8uid-solid coupling e9ect of water and soil, the in8uences of di9erent water head di9erences on the mechanical performance and deformation of segments and joints in the shield tunnel are studied. ,e water gushing from the excavation face leads to vertical convergence of the cross-sectional area of the shield tunnel, and joint opening and dislocation result in sharp decrease of the waterproof capacity of joints. Meanwhile, the stress in the vicinity of segment joints increases sharply, and local cracks occur in the segment lining. ,e axial force, shear force, and bending moment in the joint bolt are also signi;cantly increased. Based on the current metro standard and the computational results in this study, an emergency control criterion is put forward by means of controlling the discharge of water: the water head di9erence over the excavation face is required less than 4.6M.


Introduction
Elon Musk, the CEO of Tesla, announced lately that they will excavate underground tunnels to solve the problem of ground tra c congestion, and they started a company named Boring to develop a new type of shield construction machine for this purpose.Boring, led by Steve Davis, an engineer from SpaceX, is currently building an underground tunnel network in Los Angeles.In the tunnels, cars will be transported on an electrical skateboard, which could run at the maximum speed of 130 miles per hour, enough to deliver a car from Westwood to Los Angeles in 5 minutes [1].erefore, underground public transportation has great potential to alleviate urban tra c congestion.In recent years, metro construction has developed rapidly in China, and over the next ve years, new metro lines will be built in 42 cities, covering a total mileage of 4,448 km [2].In urban metro construction, shield tunneling is widely used, but at present, water inrush accidents occur frequently during the construction of the shield tunnel in China.Large amount of gushing water will inevitably cause unloading and deformation of the ground, very likely altering the stress and deformation of the shield tunnel and eventually a ecting the service performance of subway tunnels.Lee and Ishihara [3] reported that a section of tunnel in Kaohsiung collapsed because of uncontrolled leakage of groundwater spurted from the bottom of the operation face of shield tunneling as a result of pipeline or hydraulic fracturing.
In order to ensure the safety of metro tunnel shield construction, researchers have used various methods to study the stability of the excavation face during the shield construction of the metro tunnel.e related studies mainly include determining the limit of supporting pressure for the excavation face, the failure mode and mechanism of the excavation face, and the in uence of the surrounding environment on the stability of the excavation face.For example, Davis et al. [4] obtained the upper and lower bound stability solutions for heading collapse under undrained conditions, in consideration of three di erent shapes of shallow underground opening relevant to tunneling.Krause [5] calculated the limit support pressure for tunnel face failure by assuming that the failure zone was a half sphere, half circle, or quarter circle.Leca and Dormieux [6] assumed that the failure zone in front of the tunnel face consisted of a series of conical bodies, and the lower and upper bound limits were derived for the support pressure in a dry Mohr-Coulomb material.In addition, Borja [7] presented a nite element model for strain localization analysis of elastoplastic solids subjected to discontinuous displacement elds.Vermeer et al. [8] conducted a series of three-dimensional nite element simulations for tunnel face stability.
e results indicated that the failure zone in front of the tunnel face was closely related to the friction angle of sand, and the limit support pressure decreased with the friction angle of sand.Schuller and Schweiger [9] conducted a series of FEM simulations of tunnel excavation, and the results demonstrated that the development of plastic shear strains leads to a failure mechanism involving shear banding that can be captured by the model.Augarde et al. [10] investigated the stability of a plane strain heading under undrained soil conditions using nite element limit analysis methods based on the lower and upper bound theorems.Lee et al. [11] numerically analyzed the in uence of the seepage force generated by the ow of groundwater in the excavation face during tunnel construction and believed that the stability of the excavation face depends on the e ective supporting pressure acting on the soil skeleton and the seepage force induced by the water head di erence over the excavation face.e seepage force caused by groundwater ow constitutes a major part of the total supporting pressure.With the increase of groundwater table, the proportion of seepage force in the total supporting pressure tends to increase.Mollon et al. [12] analyzed the limit support pressure of the face in a circular tunnel with a kinematical approach based on a three-dimensional multiblock failure mechanism.Perazzelli et al. [13] numerically analyzed the stability of the tunnel face under seepage ow conditions with the so-called "method of slices," and the results indicated that the tensile failure of the wedge may be more critical than the shear failure if the gradient of the hydraulic head in the ground ahead of the face is high.
Most of the existing studies focus on the stability condition of the excavation face under drained or undrained conditions and the in uence of unstable excavation face on the overall stress and deformation of the shield tunnel; an equivalent model is usually used to calculate the results.However, few dealt with or considered the in uence of water inrush from the excavation surface of the shield tunnel (under drained condition) on segment joints when the excavation face has not become unstable and the ground surface has no large subsidence.e lining of the shield tunnel is composed of segments with a large number of joints, and these segment joints are the weakest part in the shield tunnel.erefore, the stress and deformation of the lining are di erent from those of the whole tunnel.e water inrush from the excavation face will not only lead to the overall deformation of tunnel but also induce the generation of joint opening and dislocation.e contact stress between the segment and the bolt at the joint will then increase sharply and weaken the joint, and the longterm safety of the tunnel will be at risk.erefore, in studying the in uence of water inrush from the excavation face on the segments of the shield tunnel, the variations of the stress and deformation in the vicinity of segment joints should be taken into account.
Several numerical models have been developed to simulate the segment joints of shield tunnels.e models can be classi ed according to the simpli cation of joints in the model as follows: simpli ed equivalent models, ring with multiple hinged joint models, and beam-spring or beamshell models.In a simpli ed equivalent model [14], the ratio of e ective lateral sti ness is introduced to reduce the overall exural sti ness of segment lining ring, so as to represent the reduction of the overall exural performance of the lining caused by the joint (including the segment joint and the longitudinal joint).In a multihinged ring model [15], the exural e ect of segment joints is simulated by a hinge structure, which may be one of two types: free hinge and elastic hinge.e in uence of girth joints is not considered but that of longitudinal joints is exaggerated.In beam-spring or beam-shell models [16][17][18], a three-dimensional shell or solid element is used to simulate the segment, and the mechanical e ects of segment and ring joints are simulated using rotating spring and shear spring.us, the planar problem is extended to a three-dimensional problem, and the longitudinal mechanical model and lateral mechanical model of the shield tunnel are uni ed; the requirement of ne analysis of the segment lining structure is also satis ed.However, it is extremely di cult to simulate the interaction between segment joints, especially to illustrate the cracking and dislocation of joints in case of water inrush from the excavation face.
is study takes as an example the water inrush from the excavation surface occurred during the shield construction in a section of railway from the Xiaodongmen station of Wuhan subway line 7 to the Wuchang station (Figure 1), to investigate the stress and deformation of metro tunnel segment joints when large amount of water gushes from the excavation face. is section from the Xiaodongmen to the Wuchang station (XWS) is set up in the north-south direction, with the left tunnel 1115.208m long and the right tunnel 1107.169m long, and mainly laid along the Zhongshan Road.e left tunnel is constructed with a Liaoning 33 shield machine, and the right tunnel is constructed with a Komatsu shield machine.An EPB (earth-pressure balance) shield tunneling machine is adopted for the excavation in the studied railway section, whose minimum radius of curvature in the horizontal plane is 400 mm, and the spacing between the left and right tunnels is 12.6-50.5m; the minimum 2 Advances in Civil Engineering radius of curvature in the longitudinal plane is 3000 m, with "V"-shaped slopes, and the thickness of the soil cover is 10.5-45.4 m. e tunnel structure is made of prefabricated reinforced concrete segments.e soil pro le in the water gushing section of the shield tunnel is shown in Figure 2. e rest of the paper is organized as follows: Section 2 introduces a 3D discontinuous contact model proposed in this study, the multiscale modeling approach, and computational parameters.Section 3 introduces and discusses the in uences of water inrush from the excavation face on the segment joints of the metro tunnel.Section 4 puts forward an emergency control criterion for shield construction under the condition of water inrush from the excavation face.Finally, Section 5 summarizes the conclusions of this study.

3D Discontinuous Contact Model.
e shield tunnel is formed by assembling segments interacting with each other at the joints.In a numerical model, the interaction between nodes is treated as a complex contact action, and the differences between a 3D discontinuous contact model and other computational models lie in the setup and description of these contact relationships.
e entire calculation model is shown in Figure 3.With the stratigraphic structure method, the soil from top to bottom is divided into 5 layers, which are miscellaneous lled soil, silty clay, strongly weathered mudstone, medium weathered mudstone, and hard mudstone, and the groundwater level is −1.20 to 2.50 m. e width, depth, and length are 100, 60, and 100 m, respectively.In this model, individual segments are simulated by solid elements, and segment caulking and ethylene propylene diene monomer (EPDM) elastic rubber-sealing gaskets are also considered in the model.
e EPDM elastomeric rubber gasket is simulated with hyperelastic solid elements, and segment bolts are simulated using beam elements with Φ30 mm circular cross-sectional area and the actual length of individual bolts.e grouting layer and surrounding rock mass are all simulated using solid elements.

Advances in Civil Engineering
According to the data of eld geological survey, the basic mechanical parameters of the soil in the water gushing area (Figure 1) during the construction of the shield tunnel are listed in Table 1.e perched water mainly occurs in the arti cially lled soil layer, receiving the atmospheric precipitation and the vertical in ltration of the living water of the surrounding residents.During the exploration, the measured initial water level of the perched water is 1.20-2.50m under the ground, and the static water level is 0.6 to 4.2 m under the ground.e pore con ned water occurs in the ancient Yangtze River layer (11-1) containing clay and ne sand, in the (11-2) layer containing clay and ne-medium sand, and in the (12) layer containing ne sand mixed with gravels and receives the osmotic recharge from the con ned water in the surrounding area.e karst ssure water mainly occurs in the karst ssures in limestone and argillaceous limestone, and it is generally distributed in karst areas and has weak local pressure-bearing capacity.Below the miscellaneous lled soil is saturated soil, indicating su cient water supply and constant water head, and it is easy to cause a large number of water gushing sections in the shield construction.
Each ring-shaped segment of subway tunnel is divided into a block K with back cover (21.5 °), two adjacent blocks B (68 °), and three standard blocks A (67.5 °).
e thickness of the segment is 350 mm, and the width of each ring is 1.5 m. e  Figure 4 shows the contact relationship among the grouting layer, segment, and the surrounding rocks in the calculation model.
e mechanical and deformation characteristics of segment joints and the interaction between the lining with back lling and the ground are simulated with the nonlinear contact theory.We used the master-slave algorithm for contact surface for contact searching and the nite sliding in the contact area for master-slave contact tracking, which is suitable for two known interface contact searches.A large number of master contact surface and corresponding slave surface need to be de ned in building a numerical model, to reduce the processing time for a large number of contacts in the model.e two ends of the connecting bolt are embedded into the corresponding segment elements to simulate the mechanical properties related to the tension, compression, and shear deformation of the joint, and the states of deformation of joints under di erent stress conditions can also be revealed.e contact property of the segment joint is "hard contact" in the normal direction, that is, the magnitude of contact stress between the contact surfaces is not limited [19].When the contact stress is zero or negative, the two contact surfaces are separated, and the contact constraints on the corresponding nodes are reduced; meanwhile, the Coulomb frictional force is applied in the tangent direction, according to the following formula: where τ crit is the critical shear stress, μ is the frictional coe cient of segment joints (i.e., 0.62), and p is the contact pressure in the normal direction.No relative sliding occurs between the friction surfaces until the shear force reaches the critical shear stress.
For the contact properties between the segment and grouting layer, the rigidity model with penalisation is used in the normal direction [19]: where p is the normal contact force, Δl is the embedded length, "kin" is the embedded rigidity with penalisation of the contact surface, and f is the penalisation function.e Coulomb frictional force is applied in the tangent direction, and the outer side of the grouting layer is tied with the soil mass.

Determination of Key Components with Multiscale Modeling Techniques.
e purpose of this work is to analyze  the in uence of water inrush from the excavation face on the internal force of the segment joints behind.Since the grid size and model scale are inversely proportional to the computational cost, large and full-size models are avoided in order to reduce the computing time.e multiscale modeling techniques [20,21] have been developed to provide an alternative to overcome this di culty.e macroscale for numerical calculation may cover a wide range, and the mesoscopic state of dangerous local structures can be observed through multiscale analysis.erefore, a hybrid model using the multiscale modeling techniques is established for numerical analysis.e simpli ed analysis method is applied to the segment model outside the drainage area, and ne-scale numerical simulation is carried out only for the drainage area.
e bolt-segment connections and mortise and tenon structures are also considered, and the other parts located certain distance away from these parts are treated by an equivalent simpli ed model with a relatively coarse mesh. is method enables the researchers to obtain the response of the tunnel and simultaneously obtain the distributions of deformation and stress in segment joints.
e equivalent model and the parts for ne analysis in the hybrid model are connected by setting the form of contact for the circumferential joint of segments, that is, the connection between the equivalent section and the part for ne analysis is considered as the circumferential joint between segments.
e constitutive model of the Coulomb frictional force is used in the tangent direction, and "hard contact" is used in the normal direction.Similarly, each segment is assembled by high-strength bolts in the circumferential direction.On the premise of accurate calculation, transition between two grids with di erent mesh sizes is avoided.As shown in Figure 5, in the drainage area, ne mesh is used for the segment bolt structure, where the internal forces of both the segment and the bolt can be obtained.e rst task is to identify the dangerous parts of the model.An equivalent model is used to analyze and determine the hazard zones.Under the in uence of water inrush from the excavation face, the deformation nephogram of the tunnel is shown in Figure 5(a).Large vertical displacement of the tunnel (Section D) is observed after excavation, and the location of this critical section is also shown in Figure 5.

Constitutive Model and Material Parameters.
According to the eld construction data, a numerical model is established, which consists of ve major parts: soil, shield machine, grouting layer, segment, and bolt.e shield body, grouting layer, and segment are assumed to be constituted by linear elastic materials.e hardening of owing slurry is simulated using equivalent layer, and segment joints and bolts are also considered; the diameter of each bending bolt is 30 mm. e speci c dimensions and basic mechanical parameters of the above components are shown in Table 2.
e shield body, grouting layer, and segment are assumed to be constituted by elastic materials, and the temporal e ect in the hardening of owing slurry is simulated with the equivalent layer.e grouting layer is divided into soft and hard grouting layers, and the thickness of the grouting layer is determined according to the actual conditions.Considering the construction parameters on site, the average grouting volume per ring is  According to the eld conditions, the thickness of the grouting layer in the construction model is assumed to be 0.14 m; therefore, according to (3) and the clearance volume of each ring of soil, the grouting ratio can be calculated to be 167%.At the same time, the grouting pressure is applied to the outer ring of the segment and the inner wall of the soil to simulate the grouting pressure at the shield tail: In addition, the model introduces a 0.04 m thick overexcavation layer on the outer wall of the shield machine, so as to simplify the in uence of the cone degree of shield machine and release the stress in the surrounding rock to some extent, making the numerical simulation closer to the actual construction process.e EPDM elastomeric gasket is modeled by the Mooney-Rivlin rubber model with two parameters: C1 and C2 are 0.7 and 0.35 MPa, respectively, and the density of gasket is 1.3 g/cm 3 .

Boundary Conditions De ned in the Numerical Model and Simulation Steps.
e seepage boundary is simulated by the uid-structure interaction of water and soil analyzed with "the transient analysis of soil" in ABAQUS.e miscellaneous back ll is dry soil, and the soil below is saturated soil.Because the speed of eld shield construction is fast and the segment has good water stopping e ect, the excavation process does not involve water loss.erefore, the pore pressure at the boundary of miscellaneous back ll is set to be zero, and except for the excavation surface, the other parts have undrained boundaries.e mechanical boundaries are set as follows: the bottom of the model is a xed boundary, only the normal displacements of the four lateral sides are restrained, and the top is a free surface.
In the process of excavation, the supporting force for the excavation face provided by the shield machine and the synchronous grouting pressure at the shield tail (as shown in Figure 6) are set as follows: the supporting force of the excavation face is 350 kPa; it increases linearly with the depth, and the gradient of increment is 20 kPa/m.e force is slightly larger than the pressure imposed by the lateral static soil, in order to prevent the ground subsidence induced by the collapse of soil above and in front of the excavation face.e grouting pressure is applied to the excavation soil and the outer wall of the segment to simulate the grouting at the shield tail.
e upper and lower grouting pressure is kept to be the same, and the value is 200 kPa.
A total of 57 analysis steps are set in the model, and each analysis step is 1-hour long.e tunneling proceeds step by step, and the width of each ring is 1.5 m. e speci c excavation steps are as follows.
Figure 7(a) is the rst step in the construction of a shield tunnel.All the shield machine body enters the soil; the shield body is 9 m long, so the soil in 9 m is inactivated simultaneously.
e supporting force for the excavation face is then added, and the outer ring (outside the overexcavation mass) of the shield machine is "tied" to the soil mass.
Figure 7(b) shows the second step: activate the front shield ring and kill the rear shield ring; at the same time, add the lining segment and the equivalent grouting layer and apply the grouting pressure to the inner wall of the soil and the outer ring of the lining; at the same time, Advances in Civil Engineering apply the supporting force for the excavation face to ensure the stability of the excavation face.In order to simplify the calculation, the Tie constraint is applied between the exterior of the slurry and the soil and between the inner part of the slurry and the outer ring of the segment.
Figure 7(c) shows that after the fth step, four rings have been supported by the segment grouting layer, and the lasting length of grouting pressure is 4 rings and 6 m.At this point, the slurry strength is 7 MPa, and the slurry behaves as a owing soft paste.
After the sixth step, the slurry at the end of the tunnel begins to harden, and the material parameters of slurry are changed using the " eld" variable in ABAQUS to increase the elastic modulus until 50 MPa; in the meanwhile, the grouting pressure behind the slurry is eliminated.Figure 7(d) shows the states of the shield machine body, lining segment, and soft-hard grouting layer after the ninth step; at this time, the segment has 8 rings, and behind the segment are 4 rings of soft paste and 4 rings of hard paste.is method is applied in sequence until the completion of 28 steps.
e water inrush from the excavation face in shield construction is a dynamic process, which gradually evolves from well-controlled drainage to small amount of drainage and then to large amount of drainage at some local segments of the excavation face.In this paper, the evolution of water inrush from the excavation face is determined by setting the water head di erence over the excavation face, and ne-mesh simulation of segments and bolts is carried out to monitor the in uence of water head di erence on the internal force of segments and bolts.Figure 8 shows the numerical model of water gushing from the local excavation face, which in combination with the condition of construction on site and the engineering geological parameters of the corresponding sections shown in Figure 3 and Table 1 suggests that when the shield tunneling initially proceeds in mudstone with good soil conditions, no water inrush occurs from the excavation face, and when passing through water-rich sand layer such as hidden river and the karst area, earth-pressure-based shield construction is easy to cause serious water inrush from the excavation face and further cause accidents such as dislocation and cracking of the rear segments.

Analysis of the Effects of Water Inrush from Shield Excavation Face
3.1.Overall Deformation of Tunnel.Figure 9 shows the distribution of deformation of tunnel lining under the condition of water inrush from the excavation face in shield tunneling.It can be seen from Figure 9 that the water inrush results in relatively large vertical displacement of segments.e maximum vertical deviation of the tunnel axis is 32 mm.e maximum vertical deviation angle appears on the tunnel axis in the dangerous Section D, and the hazard area of segments appears 2-3 rings behind the shield tail.e vertical deviation of adjacent segment vaults reaches 15 mm, that is, the maximum amount of dislocation is 15 mm.

Deformation of Joints after Water Inrush from the Excavation Face.
Water inrush from the excavation face will lead to the overall deformation of the existing lining  (1) After the occurrence of water inrush from the excavation face, joint dislocation and opening increase signi cantly in the critical Section D. is kind of joint deformation will weaken the intrinsic safety of the structure and the waterproof performance of the shield tunnel.
(2) After the occurrence of water inrush from the excavation face, dislocation usually occurs at the circumferential joints in the water gushing section, and no obvious dislocation is found at longitudinal joints.Opening occurs at longitudinal joints.e tunnel top has larger subsidence, and the volume of water gushing from the excavation face determines the size of segment dislocation and opening.Figure 10 shows the diagram of segment deformation when the excavation face is not drained and the pore water pressure of the excavation face is set to be 0, that is, the excavation face is fully drained.As can be seen from the diagram, when the excavation face is not drained, small extrusion of 1-2 mm appears at the top and bottom of the segment.However, when the excavation face is fully drained, as shown in Figure 10(b), very large compressive deformation with the maximum displacement of 32.12 mm is found at the top of the segment.Because shield tunneling works in hard mudstone, the rebound of the lower part of the tunnel is smaller, and the main deformation of segment comes from the subsidence of tunnel vault, resulting in serious change of segment internal force.In engineering practice, segment crushing failure may be caused by water inrush.
(3) e deformation of each ring of segment is shown in Figure 11.e upper part of the segment is subjected to greater extruding failure, which will cause large dislocation of circumferential segment joints and opening of longitudinal joints.Large dislocation and opening will further result in water inrush hazard of segment joints.3. e length of drainage is 6 rings, or 9 m, and the drainage time of each ring is 5 h, with reference to the undrained condition.

Variation of Joint Deformation with
Figure 12 shows the maximum deformation of segment (caused by circumferential joint dislocation and longitudinal joint opening) in the critical Section D, which changes with the progress of water inrush from the excavation face.e water head di erence over the excavation face is varied from 0 m to 18 m.e calculation results display the following variations: (1) When the water head di erence over the excavation face is 0-5 m, the deformation of the segment joint is smaller and the extents of joint dislocation and opening are smaller, but when the water head difference over the excavation face is 6-15 m, the deformation of joint increases remarkably.
(2) When the water head di erence over the excavation face reaches the limit value of 18 m, longitudinal joint opening at tunnel vault reaches the maximum value of 4 mm.Without the occurrence of water inrush, the  Advances in Civil Engineering longitudinal joint opening can be neglected.e size of joint opening exceeds the tolerance value for longitudinal joint speci ed in the metro design standard (3 mm), so water seepage may occur through longitudinal joints.
(3) When the water head di erence over the excavation face reaches the limit value of 18 m, the maximum dislocation of circumferential joint at tunnel vault reaches 15 mm.When water inrush from the excavation face does not occur, the dislocation of circumferential joint is 0.3 mm.As a result, remarkable dislocation of the circumferential joint at the tunnel vault occurs, nearly 50 times larger than that compared with the case without water inrush.When the dislocation exceeds the tolerance value of circumferential joint speci ed in the metro design standard (15 mm), the waterproof performance of the joint may be weakened.

Stress in Joint Concrete.
During shield construction of the metro tunnel, the water gushing from the excavation face will inevitably change the state of stress in segment Continuous variation of water head di erence over the excavation face results in stress concentration in segment joints.e maximum tensile stress in segment A1 is 3.1 MPa before the occurrence of water inrush and becomes 15.77 MPa after that, higher than the ultimate tensile strength of C50 concrete (5 MPa).In such cases, the local concrete will su er structural damage, and the maximum compressive stress will also change signicantly.e maximum compressive stress in segment B1 is 10.23 MPa before water inrush occurs and becomes 35.41 MPa after that.
e concrete in the vicinity of the bolt in segment joints is subjected to greater tension and shear force, because the displacement of the segment joint is larger after the occurrence of water inrush.At the same time, the maximum   Advances in Civil Engineering compressive stress in the concrete is 35.41MPa, which occurs under the maximum tensile stress near the same bolt.Special attention should be paid to the abnormal state of stress in the concrete in the vicinity of the bolt in segment joints after the occurrence of water inrush from the excavation face.e maximum tensile stress and maximum compressive stress in the concrete at segment joints change with the amount of water gushing from the excavation face, and the corresponding curves are shown in Figures 15 and 16.
(1) As shown in Figure 15, when the water head di erence over the excavation face is greater than 5 m, the stress in the concrete at the joints increases signi cantly and the maximum tensile stress in the segment will exceed the ultimate tensile strength.ere is a risk of cracking at this time.
(2) For the local compressive stress in the concrete at the joints, the in uence of the amount of water gushing from the excavation face on circumferential joints is greater than that on longitudinal joints.Considering the local tensile stress at the joints, the in uence of the water head di erence over the excavation face on circumferential joints is greater than that on longitudinal joints.

Internal Force of Coupling Bolt.
e nephograms of the tension, shear, and moment of the bolt are shown in Figures 17-19.e internal force of the bolt varies with the degree of opening and dislocation of segment joints.When the amount of water gushing from the excavation face is large, the internal force is large in high-strength 12 Advances in Civil Engineering bending bolt, and the tension-or shear-yielding damage and even fracture of bolts are easy to occur.For the common load-bearing carrier of bolts and segments, the destruction of one node will inevitably lead to progressive failure of the whole segment lining. is is extremely detrimental to the safety of the tunnel structure.erefore, researchers need to nd an approach to eliminate the inuence of the water head over the excavation face on the internal force of the segment and control the water head di erence over the excavation face to prevent the occurrence of hazards.
e variations of the maximum axial force, maximum shear force, and maximum bending moment with the progress of construction are shown in Figures 20-22.e strength and the distribution of axial force on the bolts of the segment joint change after the water inrush occurs from the excavation face.In the critical Section D, the axial force of the bolt increases sharply at the joint of the vault and arch of the tunnel.When the bolt is in a state with large deformation, its maximum tensile force reaches 316.9 kN, which is 16 times higher than that without drainage, and the maximum shear force is 561.2 kN, which is 14 times that before excavation.e maximum bending moment of bolt is 14.18 kN * m, nearly 13 times higher than that before drainage.

Criterion of Emergency Control of Water Inrush from the Excavation Face in Shield Tunneling
It is very di cult to monitor and measure the deformation and internal force of the segments and bolts in the vicinity of joints on the site of the metro tunnel.However, it is easier to monitor the amount of water gushing from the excavation face when it occurs.As shown in Figures 23-25, the in uence of di erent seepage velocities in the drainage area on the whole construction process is analyzed, and the following data are sorted out: the maximum deformation of the segment at the segment joint during the entire course of drainage in the tunnel, the maximum tensile stress in the segment, the maximum compressive stress in the segment, the maximum tensile stress in the bolt, the maximum shear stress in the bolt, and the trend of variation of the maximum bending moment with per unit Advances in Civil Engineering 13 water inrush at one step (mainly controlled by the water head di erence Δh).As can be seen from the diagram, when the water gushing from the excavation face increases, the internal force and deformation in both the segment and bolt show a trend of nonlinear increase.For a C50-reinforced concrete segment, when the stress in the segment exceeds the tensile 14 Advances in Civil Engineering strength, the segment joint may crack, causing more serious water inrush.For a high-strength M30 bolt, excessive stress will lead to the yielding failure and threaten the safety of engineering structure.erefore, in order to reduce the risk of segment fracture in case of water inrush from the excavation face, it is necessary to make emergency plans to control the water head di erence over the excavation face. is result may provide reference for the engineering practice.

Relationship between Segment Deformation and Water
Head Di erence over the Excavation Face.During water inrush from the excavation face, the in uence of water head di erence Δh on the circumferential joint is greater than that on the longitudinal joint, and the amount of dislocation caused by the deformation of the segment connecting to a joint is mainly re ected in the circumferential joint.e relationship between the maximum deformation of the segment connecting to circumferential joint and the water head di erence over the excavation face is shown in Figure 23.
e relationship between the deformation of the maximum segment joint μ and the water head di erence Δh can be expressed by a function: In order to ensure the waterproof performance of segment joint, the design speci cation of the joint deformation of the metro tunnel requires that the dislocation of circumferential joint should not be more than 15 mm.According to the calculation above, the corresponding water head di erence over tunnel lining should be less than 9.5 m. erefore, measures should be taken to ensure that the water head di erence shall not be greater than 9.5 m in case that the amount of water gushing from the excavation face continuously changed during shield construction.It helps to control the amount of circumferential joint dislocation caused by the deformation of the segments connecting to the joint and to ensure the waterproof performance of tunnel segment joints.

Relationship between Segment Internal Force and Water Head Di erence over the Excavation Face.
e relationship between the tensile stress and compressive stress in the concrete constituting segment joints and the water head di erence over the excavation face is shown in Figure 24, where σ t represents the maximum tensile stress in the concrete and σ p represents the maximum compressive stress in the concrete.
e relationship between the maximum stress in the concrete and the water head di erence index can be expressed by the following functions: e uniaxial compressive strength of the C50 concrete for subway segment is 50 MPa, and the corresponding water head di erence over the excavation face is less than 24.5 m.
ese values indicate that the concrete in the vicinity of the segment joint will not be crushed.
e splitting tensile strength test [22] is conducted for the C50 segment concrete used in the Wuhan no.7 subway line, and the obtained ultimate tensile strength is 5 MPa, and the corresponding water head di erence over the excavation face is less than 4.6 m. ese values indicate that the concrete in the vicinity of the segment joint will not be damaged by excessive tension.Hence, in order to ensure the safety of segment joints, the water head di erence over the excavation face should be less than 4.6 m.

Relationship between the Internal Force of Bolts and the
Water Head Di erence over the Excavation Face.When the amount of water gushing from the excavation face keeps varying in the process of shield construction, the relationship between the internal force of bolts connecting adjacent segments and the water head di erence over the excavation face is also shown in Figure 25, where N represents the maximum axial force acting on the bolt, Q represents the maximum shear force, and M represents the maximum bending moment.e relationship between e minimum tensile strength of the 8.8 grade bolt reaches 830 MPa, and the shear strength reaches 320 MPa.
erefore, the axial force and shear force of the M30 bolt are 466 and 215 kN, respectively, and the corresponding water head di erence is 22.9 m and 8.8 m, respectively.erefore, in order to avoid the yielding of bolts connecting adjacent tunnel segments, the water head di erence over the excavation face should be less than 8.8 m.
In a word, in order to guarantee the long-term safety of segments and joints of a metro shield tunnel, the water head di erence over the excavation face should not be greater than 4.6 m under the condition of water inrush from the excavation face.e water head di erence can be adjusted by controlling the discharge of water on site.

Conclusions
When water inrush occurs in the process of shield construction of the metro tunnel, the water head di erence over the excavation face varies.With the aim of studying the inuences of continuous water inrush from the excavation face on the mechanical properties and deformation of segment joints (including the deformation and internal force of the segments and bolts in the vicinity of joints), a multiscale hybrid model of the excavation face in a shield tunnel is constructed based on the three-dimensional nonlinear contact theory and applied to simulate a section of XWS. e results of this study may provide reference for the engineering practice of shield tunneling.e following conclusions are drawn: (1) e vertical convergence of the shield tunnel mainly re ects the overall deformation of the existing tunnel during the continuous change of water in ow at the excavation face.Dislocation of circumferential joints and opening of longitudinal joints are observed, and after excavation, the extent of both segment joint opening and dislocation is 16.5 times higher than that when the excavation face has no drainage.is situation directly results in a signi cant reduction in the e ectiveness of the waterproof function of joints.(2) e stress at the segment joints and connecting bolts changes sharply because of the constant change of water inrush at the excavation face.e local tensile stress of the segment exceeds the tensile strength limit of the C50 concrete; the axial force of the bolt is greater at the crown and the waist; in the critical position D, the shear force of the bolt on the circumferential joint increases signi cantly, and the maximum shear force of the bolt reaches 561.2 kN.
(3) e emergency control standard should be established based on the water head di erence Δh of the excavation face during the tunnel construction.When the water head di erence is over 4.6 m with continuously changing rate of water inrush from the excavation face, the limits of deformation and stress of segment joints are considered to have reached, and the risks of segment dislocation and extrusion cracking are increased.e research team will conduct further centrifuge simulation experiments to enhance the reliability of the calculation results.

Figure 1 :
Figure 1: Location of the XWS in the Wuhan metro system.

Figure 2 : 2 S tr o n g ly w e a th e r e d m u d s to n e 1 9 − 1 M o d e r a t e ly w e a t h e r e d m u d s t o n e 1 9 c − 2 H a r d m u d s t o n e 1 9 c − 1 DFigure 3 :
Figure 2: e water gushing section of XWS.

Figure 4 :
Figure 4: e contact relationship among the grouting layer, segment, and the surrounding rocks.(a) Model of the segmental lining.(b) 3D contact relationship.

DFigure 5 :
Figure 5: Multiscale hybrid model.(a) Critical section of the tunnel under water inrush.(b) Hybrid model of critical Section D.

Figure 6 :
Figure 6: e supporting force by shield machine and grouting pressure at the shield tail.(a) e supporting force at the excavation face.(b)e grouting pressure at the shield tail.
g segm ent linin g e leng th of shie ld ma chi ne = 9 m So ene d pas te = 6 m e shiel d body (c) e shiel d body So ene d pas te Ha rde ned pas te 8-rin g segm ent linin g (d)

Figure 7 :
Figure 7: e simulation of the excavation of the shield tunnel.(a) e rst step.(b) e second step.(c) e fth step.(d) e sixth step.
Di erent Amounts of Water Gushing.With the continuous change of the amount of water gushing from the excavation face, the deformation of segment joints change continuously.Seven values of pore water pressure of the excavation face ("por") are set (0, 20, 40, 70, 100, 130, and 160) to monitor the change of segment displacement and internal force.e pore water pressure values 0, 20, 40, 70, 100, 130, and 160 correspond to the reduction of water head (Δh) of 18, 16, 14, 11, 8, 5, and 2 m, respectively, as shown in Table

e dr ain ag e ar eaFigure 8 :
Figure 8: Local area of water inrush from the local excavation face.

Figure 9 :
Figure 9: Deformation distribution of tunnel lining after water inrush.

Figure 10 :Figure 11 :
Figure 10: e segment deformation under di erent water drains.(a) No draining at the excavation face.(b) Fully drained at the excavation face.

Figure 13 :
Figure 13: e e ect of water head di erence on the maximum tensile stress of segment.Water head di erence of (a) 18 m and (b) 0 m.

Figure 12 :
Figure 12: e maximum deformation of segment.

Figure 14 :Figure 15 :
Figure 14: e e ect of water head di erence on the maximum compressive stress of segment.Water head di erence of (a) 18 m and (b) 0 m.

Figure 25
Figure 24: e relationship between the tensile and compressive stress in the concrete constituting segment joints and the water head di erence over the excavation face.(a) e maximum tensile stress.(b) e maximum compressive stress.

Table 1 :
Basic parameters of soil mass. (b)
7 MPa; C2 0.35 MPa; q 1.3 g/cm26Advances in Civil Engineering set to be 7 m 3 ; considering the compressibility of the slurry at the end of shield and the volume loss caused by water loss, the grouting volume is larger than that of the clearance of soil over excavation, and the grouting volume is generally 150%-200% of the construction clearance.

Table 3 :
Pore pressure and the corresponding water head di erence at the excavation face.