A Case Study on the Application of the Steel Tube Slab Structure in Construction of a Subway Station

: It is an effective approach to use Steel Tube Slab (STS) structure combined with the Pile-Beam-Arch (PBA) method to construct a large-space underground station. Traditional construction methods cannot meet the requirement of construction because of the complicated soil layers and high building densities in urban areas. The STS method can effectively increase the rigidity of the supporting system by using steel tubes. Firstly, the stress of bolts and steel tubes are investigated in the construction process based on the ﬁeld monitoring data. Subsequently, FLAC3D is used to establish a three-dimensional model, which is veriﬁed based on the in-situ monitoring data; the effect of excavation process on ground settlement, deformation of STS structure and bridge pile are studied by numerical results. Moreover, the key parameters such as welding of ﬂanges and the step length are studied. The results show that the stress of the steel tubes and ﬂanges does not exceed the designed strength during the construction process. Based on the numerical simulation data, it is indicated that the STS structure can be a very effective and dependable measure in controlling and reducing the surface settlement and the existing adjacent buildings. The numerical results can be used to guide the later construction.


Introduction
The stress redistribution induced by underground excavations causes soil movements around the excavation face in the earth mass, which is also translated to ground settlement. These soil movements can cause severe damage to adjacent buildings and underground facilities [1,2]. Control of the soil movements due to excavations is of utter importance in subway projects, because they are often planned in the populated areas of cities that are already well developed. Several methods for surface settlement prediction due to underground excavations are currently available including those proposed by Peck, Reilly and Novozhenin [3][4][5].
Underground excavations in urban area have some characteristic features such as large spans, ultra-shallow buried depth, complex environments, and heavy overburdens due to street traffic. Hence, excavation-induced instability and settlement has become a major concern in the design and construction of the new underground excavations [6].
The pipe roofing method evolved from pipe umbrella method adopted in shallow ground cover and poor ground conditions during underground construction. As for this method, steel pipes are jacked to the design location by pipe jacking machines to act as the initial support system. These adjacent pipes are connected together and form a supporting structure, filled with grout or concrete in the inside and outside of the steel pipes before tunnel excavation. The pipe roofing method was firstly to be used to construct an underpass railway in 1971 and a total of six such projects were completed in the next decade [7]. The Antwerp subway station in Belgium was an early engineering feature, constructed by applying pipe roofing method in Europe in 1979 [8,9]. It was first used in America in 1994 for a large shallow covered tunnel [10]. Subsequently, this method is applied to construct highways in Japan, and on the basis of the pipe roofing method, Endless Self-Advancing (ESA) and Front Jacking (FJ) were proposed [11,12]. Meanwhile, the methods were widely applied to certain engineering projects in Singapore and Malaysia [7,13]. A typical engineering project, 923 stations of the Seoul subway line 9 in Korea was constructed by using the pipe roofing method in the complex environments [14,15].
Based on some engineering examples of previous pipe roofing method in China, it was used to build an underpass in Hong Kong in 1984, and the Taipei Songshan airport underpass was constructed using the pipe roofing method in 1989. The pipe roofing method was employed to construct the vehicle tunnel, which passes the Taipei International Airport [16]. The Roof-Box Jacking Method (RBJ) was used for the underpass project of Shanghai Beihong street in 2004. This support structure formed a rectangular section and was made up of 80 steel pipes with a diameter of 0.97 m and longitudinal length of 125 m [17]. By an arc-shaped supporting system, Xinle Yizhi station of Shenyang subway line 2 was constructed by the new roof tubular of method with 19 steel pipes with diameter of approximately 2 m in 2011 [18]. The Gongbei tunnel, as a part of Hongkong-Zhuhai-Macao Bridge project, was constructed a combination technology of curved pipe jacking roof for ground support and ground freezing for waterproof. This pipe roofing structure is composed of 36 steel pipe strings of 1.620 m diameter with a 355-358 mm gap between adjacent strings [19,20].
In term of studies on the working mechanism of the pipe-roofing support systems, these have shown that the maximum bending moment of the steel grid is often much larger than that of the pipe-roofing system in the longitudinal direction [21]. Therefore, to ensure the stability of the pipe-roofing system during excavation, an intensive temporary supporting structure is often required. The tedious and expensive set up of the steel supporting structure is a big disadvantage of the traditional pipe-roofing method. And as for forming an arch structure, the space utilization is lower than the rectangular section.
To overcome the problem associated with supporting system requirements as mentioned above, a new pipe-roofing structure, also known as the Steel Tube Slab (STS) structure, has been proposed [22]. STS design considerably improves the transverse load-carrying capacity of the traditional pipe-roofing system by inclusion of high strength transverse bolts and flange plates between the tubes. Such a system is attractive as it can avoid the use of temporary support during excavation.
In the current study, the application of the proposed method in the design and construction of Northeastern Street station of line 10 of the Shenyang Subway is then thoroughly discussed. Data from monitoring stations in the construction site are presented and compared with the results of the numerical simulations. The effect excavation process on ground settlement, deformation of STS structure and bridge pile are studied; moreover the key parameter such as welding of flanges and the step length are investigated by conducting a suite of numerical simulations.

Support Structures in Underground Excavations
Different types of supporting structures have so far been used to control the soil movement during underground excavations. These methods are briefly discussed in the following sections.

Pile, Beam and Arch (PBA) Method
The concept of the traditional Pile, Beam and Arch (PBA) method is to develop a supporting structure in underground excavations to control soil movements and ground settlements. In this Appl. Sci. 2018, 8,1437 3 of 29 method, as shown in Figure 1, pilot tunnels are first constructed to obtain access for manual-excavation and pile construction. The top of the piles is then attached by a cap beam. The soil is then excavated by bench method to form arches whose ends are braced against the cap beams. The soil is then excavated in a top-down approach and medium level floors are built as necessary. Finally, the bottom floor is constructed after the last layer of the soil is excavated. The PBA method is very effective and has been used in many projects in the past [23,24].
Appl. Sci. 2017, 7, x 3 of 29 soil is then excavated by bench method to form arches whose ends are braced against the cap beams. The soil is then excavated in a top-down approach and medium level floors are built as necessary. Finally, the bottom floor is constructed after the last layer of the soil is excavated. The PBA method is very effective and has been used in many projects in the past [23,24].

The Traditional Pipe Roofing Structure
Another approach used for underground supporting system is traditional pipe roofing method. In this technique, adjacent steel tubes are connected by interlock to form the pipe roofing structure, and then temporary supports are established while the soil below the pipe roofing is excavated. Previous applications of this method show that it is effective in reducing surface deformation [25,26]. However, in this technique, adjacent steel tubes are only connected by interlock, resulting in low

The Traditional Pipe Roofing Structure
Another approach used for underground supporting system is traditional pipe roofing method. In this technique, adjacent steel tubes are connected by interlock to form the pipe roofing structure, and then temporary supports are established while the soil below the pipe roofing is excavated. Previous applications of this method show that it is effective in reducing surface deformation [25,26]. However, As for this new supporting structure, the authors have comprehensively investigated the performance of the new STS structure and the effect of different parameters on its stiffness and limit load using both numerical and experimental techniques. The results of these investigations will not be repeated here. Interested readers are referred to the authors [22,28,29] for more details.

Engineering Overview
While several stations of the Shenyang subway (line 10) are constructed using the traditional PBA method, for the construction of the Northeastern Street station (located west of the intersection of Beihai street and Northeastern Road) the traditional PBA method could not be used due to the complexity of the existing structures close to the station, and the heavy traffic load that had to be

The Proposed Supporting Structure
In this section, a new technique for underground excavations supporting system is presented. The technique proposed is based on the STS method, however, unlike the traditional pipe roofing method, it is proposed that the adjacent steel tubes are connected by flange plate, bolt and concrete [27], as shown in Figure 3. This will significantly increase the lateral stiffness and bearing capacity of the pipe roofing system. Therefore, the need for other temporary supports during excavation is circumvented, which shortens the construction period, reduces the complexity of the construction and increase the safety during construction. As will be shown in the remainder of this paper, the proposed technique is superior to the standard PBA and the traditional pipe roofing methods, particularly in large span and ultra-shallow underground excavations.
Appl. Sci. 2017, 7, x 4 of 29 lateral stiffness of the pipe-roof structure. As a result, other temporary supports must be established to maintain the stability of the whole structure when the soil is excavated, increasing the construction complexity and cost, as shown in Figure 2.

The Proposed Supporting Structure
In this section, a new technique for underground excavations supporting system is presented. The technique proposed is based on the STS method, however, unlike the traditional pipe roofing method, it is proposed that the adjacent steel tubes are connected by flange plate, bolt and concrete [27], as shown in Figure 3. This will significantly increase the lateral stiffness and bearing capacity of the pipe roofing system. Therefore, the need for other temporary supports during excavation is circumvented, which shortens the construction period, reduces the complexity of the construction and increase the safety during construction. As will be shown in the remainder of this paper, the proposed technique is superior to the standard PBA and the traditional pipe roofing methods, particularly in large span and ultra-shallow underground excavations. As for this new supporting structure, the authors have comprehensively investigated the performance of the new STS structure and the effect of different parameters on its stiffness and limit load using both numerical and experimental techniques. The results of these investigations will not be repeated here. Interested readers are referred to the authors [22,28,29] for more details.

Engineering Overview
While several stations of the Shenyang subway (line 10) are constructed using the traditional PBA method, for the construction of the Northeastern Street station (located west of the intersection of Beihai street and Northeastern Road) the traditional PBA method could not be used due to the complexity of the existing structures close to the station, and the heavy traffic load that had to be As for this new supporting structure, the authors have comprehensively investigated the performance of the new STS structure and the effect of different parameters on its stiffness and limit load using both numerical and experimental techniques. The results of these investigations will not be repeated here. Interested readers are referred to the authors [22,28,29] for more details.

Engineering Overview
While several stations of the Shenyang subway (line 10) are constructed using the traditional PBA method, for the construction of the Northeastern Street station (located west of the intersection of Beihai street and Northeastern Road) the traditional PBA method could not be used due to the complexity of the existing structures close to the station, and the heavy traffic load that had to be handled during the excavation. As a result, with considerations based on geographical location and geologic situation of the subway station, it was decided to use a combination of STS and PBA methods, as discussed in the previous section, for the design and the construction of the station.
The total length and standard cross section of the subway station is 225.95 m and 24.7 m, respectively. And the buried depth of the bottom floor of the platform layer is approximate 17.5 m, with an overburden layer of about 3.4 m thick. For the construction, the station is divided into three sections: north of the station for which covered excavation is used; south of the station for which open surface method is adopted; and the middle section of the station for which a combination of STS and PBA methods is adopted due to the intensity of the surrounding structures and heavy surface traffic load. The length of the southern, middle, and the northern segments of the station are 68.65 m, 114.4 m, 42.9 m, respectively. The engineering details of the Northeastern Street station is shown in Figure 4. handled during the excavation. As a result, with considerations based on geographical location and geologic situation of the subway station, it was decided to use a combination of STS and PBA methods, as discussed in the previous section, for the design and the construction of the station. The total length and standard cross section of the subway station is 225.95 m and 24.7 m, respectively. And the buried depth of the bottom floor of the platform layer is approximate 17.5 m, with an overburden layer of about 3.4 m thick. For the construction, the station is divided into three sections: north of the station for which covered excavation is used; south of the station for which open surface method is adopted; and the middle section of the station for which a combination of STS and PBA methods is adopted due to the intensity of the surrounding structures and heavy surface traffic load. The length of the southern, middle, and the northern segments of the station are 68.65 m, 114.4 m, 42.9 m, respectively. The engineering details of the Northeastern Street station is shown in Figure 4.   The soil stratigraphy at the station location was determined from geological and geotechnical investigations. Accordingly, the soil deposit can be divided into five layers based on the mechanical properties and formation lithology, as shown in Figure 5. And depth of the water table is from 5.2 m to 14.9 m in the construction area. In order to get the strength parameters of sandy soils, a series of super heavy Dynamic cone Penetration test (DPSH) were performed, as shown in Table 1.
Appl. Sci. 2017, 7, x 6 of 29 The soil stratigraphy at the station location was determined from geological and geotechnical investigations. Accordingly, the soil deposit can be divided into five layers based on the mechanical properties and formation lithology, as shown in Figure 5. And depth of the water table is from 5.2 m to 14.9 m in the construction area. In order to get the strength parameters of sandy soils, a series of super heavy Dynamic cone Penetration test (DPSH) were performed, as shown in Table 1.  Several steps are involved in the construction of the proposed supporting structure based on combination of STS and PBA methods. These steps are demonstrated in details in Figure 6.
(1) The steel tubes are pushed into the soil by the jacking machine, and then the soil inside the steel tubes are extracted using horizontal directional drilling (Figure 6a).
(2) Lateral connecting bolts are installed between adjacent steel tubes by artificial operations. Blocking plate at both ends of the steel tubes are then installed, and concrete is poured through grouting holes on the blocking plates ( Figure 6b).
(4) The inter grid of the top pilot tunnels is removed, and the waterproof layer at the top and the middle floor slab are constructed. When the strength of the floor slab reaches to the 80% of the design strength, the soil above the middle floor slab is excavated for the whole cross section of the station, followed by the construction of the middle longitudinal beams and side walls (Figure 6d).  Several steps are involved in the construction of the proposed supporting structure based on combination of STS and PBA methods. These steps are demonstrated in details in Figure 6.
(1) The steel tubes are pushed into the soil by the jacking machine, and then the soil inside the steel tubes are extracted using horizontal directional drilling (Figure 6a).
(2) Lateral connecting bolts are installed between adjacent steel tubes by artificial operations. Blocking plate at both ends of the steel tubes are then installed, and concrete is poured through grouting holes on the blocking plates ( Figure 6b).
(4) The inter grid of the top pilot tunnels is removed, and the waterproof layer at the top and the middle floor slab are constructed. When the strength of the floor slab reaches to the 80% of the design strength, the soil above the middle floor slab is excavated for the whole cross section of the station, followed by the construction of the middle longitudinal beams and side walls (Figure 6d).
(5) When the strength of the middle floor slabs and side walls of the station hall reach up to 80% of the design values, the soil from the middle span to the bottom slab is excavated. Followed by the construction of the waterproof layer and the bottom slab (Figure 6e).  (6) The rest of side walls and waterproof layers of top slabs are constructed to complete the construction of the subway station ( Figure 6).
Appl. Sci. 2017, 7, x 7 of 29 (5) When the strength of the middle floor slabs and side walls of the station hall reach up to 80% of the design values, the soil from the middle span to the bottom slab is excavated. Followed by the construction of the waterproof layer and the bottom slab ( Figure 6e).
The rest of side walls and waterproof layers of top slabs are constructed to complete the construction of the subway station.
(e) (f) Figure 6. Construction stages in the new support system. Note: Z1, Z2, Z3 and Z4 represent pilot tunnel; C1 and C2 represent cylindrical steel column. The characteristics of this construction method combining STS and PBA methods are as follows: (1) The support structure can be designed to have adequate bearing capacity to provide support for the overburden pressure when soil below the structure is excavated. So other temporary supports may not be required resulting in a flexibility of the method and potential for wide applications.
(2) Not only can the proposed supporting system effectively control the soil surface deformation, but also the intensive temporary supports during excavation stage can be avoided owing to support provided by the cylindrical columns involved in the proposed technique.
(3) The proposed method results in a rectangular cross section for the excavated area, which increases the utilization of the underground space compared to the traditional arch formation.

Connection of Adjacent Steel Tube
For the Northeastern Street station, the design specifications required the outer diameter and the thickness of the steel tubes of 0.9 m and 16 mm, respectively, and the center distance between two adjacent steel tubes of 1. The characteristics of this construction method combining STS and PBA methods are as follows: (1) The support structure can be designed to have adequate bearing capacity to provide support for the overburden pressure when soil below the structure is excavated. So other temporary supports may not be required resulting in a flexibility of the method and potential for wide applications.
(2) Not only can the proposed supporting system effectively control the soil surface deformation, but also the intensive temporary supports during excavation stage can be avoided owing to support provided by the cylindrical columns involved in the proposed technique.
(3) The proposed method results in a rectangular cross section for the excavated area, which increases the utilization of the underground space compared to the traditional arch formation.

Connection of Adjacent Steel Tube
For the Northeastern Street station, the design specifications required the outer diameter and the thickness of the steel tubes of 0.9 m and 16 mm, respectively, and the center distance between two adjacent steel tubes of 1.

Connection between Cylindrical Steel Columns and Longitudinal Beams
The bottom and top longitudinal beams are made of C40 concrete, which are directly poured in place, in connection with the cylindrical steel columns. Sand is poured into the space between the cylindrical steel columns and the surrounding soil to increase the stability of the cylindrical steel columns (see Figure 8a). The cylindrical steel columns and longitudinal beams are connected by a reinforcement cage, which is adopted in the cylindrical columns. Furthermore, flanges are used to provide further strength in the connection between longitudinal beam and cylindrical column, as shown in Figure 8b. Besides, anchoring bolts are lengthened and bent at the ends to increase anchorage force.

Connection between Cylindrical Steel Columns and Longitudinal Beams
The bottom and top longitudinal beams are made of C40 concrete, which are directly poured in place, in connection with the cylindrical steel columns. Sand is poured into the space between the cylindrical steel columns and the surrounding soil to increase the stability of the cylindrical steel columns (see Figure 8a). The cylindrical steel columns and longitudinal beams are connected by a reinforcement cage, which is adopted in the cylindrical columns. Furthermore, flanges are used to provide further strength in the connection between longitudinal beam and cylindrical column, as shown in Figure 8b. Besides, anchoring bolts are lengthened and bent at the ends to increase anchorage force.

Monitoring Scheme
JTM-V5000c strain meter made by Jiangsu gold civil construction group, China were used to monitor strain of steel tubes. The working range of the sensors is from 1200 με to 800 με, with the accuracy of 0.1 με. The dates were collected by a digital auto-reading instrument of DH3816N made by Donghua instrument group, China, and the stress values are calculated by the existing computing parameters. The installation of strain meter is shown in Figure 9. Considering circumferential strain meters merely measure the chord length corresponding to a circular arc of inner wall of steel tubes, in order to analyze the measurement error of circumferential strain, the feasibility and accuracy of monitoring schedule is investigated. Based on the installation in Figure 10, the chord length and in corresponding to ach length are D and l, respectively, and central angle and the inner radius are 2α and r. In the jacking stage, the strain of chord length is εD, and the circumferential deformation and strain are Δl and εr. So, the relative error is calculated based on the Formula 1 [30].

Monitoring Scheme
JTM-V5000c strain meter made by Jiangsu gold civil construction group, China were used to monitor strain of steel tubes. The working range of the sensors is from 1200 με to 800 με, with the accuracy of 0.1 με. The dates were collected by a digital auto-reading instrument of DH3816N made by Donghua instrument group, China, and the stress values are calculated by the existing computing parameters. The installation of strain meter is shown in Figure 9. Considering circumferential strain meters merely measure the chord length corresponding to a circular arc of inner wall of steel tubes, in order to analyze the measurement error of circumferential strain, the feasibility and accuracy of monitoring schedule is investigated. Based on the installation in Figure 10, the chord length and in corresponding to ach length are D and l, respectively, and central angle and the inner radius are 2α and r. In the jacking stage, the strain of chord length is εD, and the circumferential deformation and strain are Δl and εr. So, the relative error is calculated based on the Formula 1 [30]. Considering circumferential strain meters merely measure the chord length corresponding to a circular arc of inner wall of steel tubes, in order to analyze the measurement error of circumferential strain, the feasibility and accuracy of monitoring schedule is investigated. Based on the installation in Figure 10, the chord length and in corresponding to ach length are D and l, respectively, and central angle and the inner radius are 2α and r. In the jacking stage, the strain of chord length is ε D , and the circumferential deformation and strain are ∆ l and ε r . So, the relative error is calculated based on the Formula 1 [30]. Based on the results of in-situ monitoring results, the monitoring values are within the limits of safety values. Applying the length of strain meter, inner diameter of steel tube and the monitoring value of strain meter to the Formula (1), the relative error between the linear strain monitored by strain meter and the circumferential strain of steel tube is 0.215%, and it is permitted. It is therefore reasonable for the installation to monitor the circumferential stain.
In order to study the force behavior and the force differences between different jacking distance in the jacking construction, the two monitoring sections A and B are arranged in Steel tube 5#, and the distance from the initial end are 9 m and 15 m, respectively. Four monitoring points are prepared in every section, and each point is installed two strain meters, one longitudinal direction and one circumferential direction. The point above the steel tube is defined as 1#, and the others point are defined as 2#, 3# and 4# according to a clockwise direction, as shown in Figure 11.  Based on the results of in-situ monitoring results, the monitoring values are within the limits of safety values. Applying the length of strain meter, inner diameter of steel tube and the monitoring value of strain meter to the Formula (1), the relative error between the linear strain monitored by strain meter and the circumferential strain of steel tube is 0.215%, and it is permitted. It is therefore reasonable for the installation to monitor the circumferential stain.
In order to study the force behavior and the force differences between different jacking distance in the jacking construction, the two monitoring sections A and B are arranged in Steel tube 5#, and the distance from the initial end are 9 m and 15 m, respectively. Four monitoring points are prepared in every section, and each point is installed two strain meters, one longitudinal direction and one circumferential direction. The point above the steel tube is defined as 1#, and the others point are defined as 2#, 3# and 4# according to a clockwise direction, as shown in Figure 11.
Strain meter Locker r Steel tube Figure 10. Calculation of measurement error strain meter.
Based on the results of in-situ monitoring results, the monitoring values are within the limits of safety values. Applying the length of strain meter, inner diameter of steel tube and the monitoring value of strain meter to the Formula (1), the relative error between the linear strain monitored by strain meter and the circumferential strain of steel tube is 0.215%, and it is permitted. It is therefore reasonable for the installation to monitor the circumferential stain.
In order to study the force behavior and the force differences between different jacking distance in the jacking construction, the two monitoring sections A and B are arranged in Steel tube 5#, and the distance from the initial end are 9 m and 15 m, respectively. Four monitoring points are prepared in every section, and each point is installed two strain meters, one longitudinal direction and one circumferential direction. The point above the steel tube is defined as 1#, and the others point are defined as 2#, 3# and 4# according to a clockwise direction, as shown in Figure 11.

Relationship between Stress and Jacking Distance
Based on the monitoring results, the strain of steel tube is elastic in the jacking stage. So the relationship between the stress and strain is defined in the elasticity stage based on the Hooke law, as shown in Formula (2) [31].
However, the force mechanism is attributed to plane stress problem, σ z = 0, as shown in Formula (3) [31].
The formula can be transformed into another form, as shown in Formula (4) where, E is the elasticity modulus of steel tube; νis the poisson ratio; σ x , σ y and σ z are the longitudinal strain, the circumferential strain and the radius strain, respectively. Based on the in-situ monitoring data, the longitudinal and circumferential stain (ε x and ε y ) is obtained, and the stress values with the two direction are required using Formula (4).
The longitudinal and circumferential stress are shown in Figures 12 and 13, along with the increase of the jacking distance in the construction stage. Based on the monitoring results, the stress trend of section A is similar to section B. The results are thus analyzed according to section A. As point 3 is damaged in the construction stage, it is not analyzed.

Longitudinal Stress of Steel Tube at Section A
The longitudinal stress is shown in Figure 12. From the figure, the left and right side is the compressive stress, and the inner side is larger than the outer side because the steel tube turns right, which leads to compressive stress concentration at right side. Although the bottom of the steel tube performs the tensile strain, but the final longitudinal stress is tensile stress, this is because it is less than the circumferential stress and affected by Poisson effect. Therefore, the longitudinal stress is controlled by the lacking of force in the horizontal direction.

Circumferential Stress of Steel Tube at Section A
The circumferential stress is shown in Figure 13. As it is indicated that the bottom of the steel tube reflects the tensile stress, the left and right side is the compressive stress. Moreover, the value of the left and right side is less than the bottom, and the inner side is larger the outer side. This is because the circumferential tensile strain is larger than the longitudinal compressive strain, which decrease the influences of longitudinal compressive deformation. As for the inner stress is larger than the outer side, this is because the outer resistance of soil is usually larger than the inner side, and the right pressure decreases the compression deformation, which lead to the decrease of the compressive stress.

Monitoring Scheme
Since the proposed STS method had not been applied to construct a shallow buried subway station before, stress monitoring for the bolts connecting adjacent steel tubes was planned for the support structure to gain confidence in the theoretical and practical aspects of the design and performance. According to the previous experimental research about STS structure, when failure is

Circumferential Stress of Steel Tube at Section A
The circumferential stress is shown in Figure 13. As it is indicated that the bottom of the steel tube reflects the tensile stress, the left and right side is the compressive stress. Moreover, the value of the left and right side is less than the bottom, and the inner side is larger the outer side. This is because the circumferential tensile strain is larger than the longitudinal compressive strain, which decrease the influences of longitudinal compressive deformation. As for the inner stress is larger than the outer side, this is because the outer resistance of soil is usually larger than the inner side, and the right pressure decreases the compression deformation, which lead to the decrease of the compressive stress.

Circumferential Stress of Steel Tube at Section A
The circumferential stress is shown in Figure 13. As it is indicated that the bottom of the steel tube reflects the tensile stress, the left and right side is the compressive stress. Moreover, the value of the left and right side is less than the bottom, and the inner side is larger the outer side. This is because the circumferential tensile strain is larger than the longitudinal compressive strain, which decrease the influences of longitudinal compressive deformation. As for the inner stress is larger than the outer side, this is because the outer resistance of soil is usually larger than the inner side, and the right pressure decreases the compression deformation, which lead to the decrease of the compressive stress.

Monitoring Scheme
Since the proposed STS method had not been applied to construct a shallow buried subway station before, stress monitoring for the bolts connecting adjacent steel tubes was planned for the support structure to gain confidence in the theoretical and practical aspects of the design and performance. According to the previous experimental research about STS structure, when failure is

Monitoring Scheme
Since the proposed STS method had not been applied to construct a shallow buried subway station before, stress monitoring for the bolts connecting adjacent steel tubes was planned for the support structure to gain confidence in the theoretical and practical aspects of the design and performance. According to the previous experimental research about STS structure, when failure is going to happen for the component beams, the bolts are also damaged simultaneously [32]. So the construction risk can be predicted by the bolts of field observation.
JTM-V1000 vibrating strings (2014, Jiangsu Gold Civil Construction Group, Changzhou, China) were used to monitor Stress in the bolts connecting adjacent steel tubes. The working range of the strings is from −100 MPa to 200 MPa, with the accuracy of 0.1 MPa. The dates were collected by a digital auto-reading instrument of GPC-3 type made by Dandong three of instrument group, China. The monitoring point were selected based on the potential risk areas, as shown in Figure 14. The initial readings of the strings were recorded prior to the removal of the soil beneath the STS structure.
Considering Tube 9 and Tube 10 are over the top pilot tunnel (Z3), and are located in one of the most disadvantage region in the pilot construction stage, as shown in Figure 3. So 8 monitoring sensors were installed between Tube 9 and Tube 10 to investigate mechanical behavior of the bolts inside steel tubes when the soil was excavated. Each monitoring section was arranged with two force sensors to verify the correctness of the data. The detail of the installation of the monitoring sensors is shown in Figure 15. going to happen for the component beams, the bolts are also damaged simultaneously [32]. So the construction risk can be predicted by the bolts of field observation. JTM-V1000 vibrating strings (2014, Jiangsu Gold Civil Construction Group, Changzhou, China) were used to monitor Stress in the bolts connecting adjacent steel tubes. The working range of the strings is from −100 MPa to 200 MPa, with the accuracy of 0.1 MPa. The dates were collected by a digital auto-reading instrument of GPC-3 type made by Dandong three of instrument group, China. The monitoring point were selected based on the potential risk areas, as shown in Figure 14. The initial readings of the strings were recorded prior to the removal of the soil beneath the STS structure.
Considering Tube 9 and Tube 10 are over the top pilot tunnel (Z3), and are located in one of the most disadvantage region in the pilot construction stage, as shown in Figure 3. So 8 monitoring sensors were installed between Tube 9 and Tube 10 to investigate mechanical behavior of the bolts inside steel tubes when the soil was excavated. Each monitoring section was arranged with two force sensors to verify the correctness of the data. The detail of the installation of the monitoring sensors is shown in Figure 15.

Analysis of Monitoring Results
The stress on the bolts connecting Tube 9 and Tube 10 were monitored during the construction sequence of the station. The stresses recorded in the eight monitoring sensors at different construction stages are shown in Figure 16. The following observations can be made at different construction stages: (1) Pilot tunnels construction: The bolts started going in tension as the excavation started, and the tensile stresses gradually increased along with the soil excavation. When the pilot tunnels were excavated at the distance of 15m, the tensile stress of about 5 kN was recorded in the monitoring point 4#. The tensile stress in the bolts began to increase rapidly when the pilot tunnels were excavated further, and the peak value of 28.9 kN was recorded when the tunnel excavation reached to its final distance of 43.5 m.
As seen in Figure 16, with the tunnel excavation, the monitoring values increase gradually, and the increment shows up an increasingly large trend. The reasons were as follows: (a) The influence of steel pipe jacking in accuracy in the pre-construction stage. According to the construction records, going to happen for the component beams, the bolts are also damaged simultaneously [32]. So the construction risk can be predicted by the bolts of field observation. JTM-V1000 vibrating strings (2014, Jiangsu Gold Civil Construction Group, Changzhou, China) were used to monitor Stress in the bolts connecting adjacent steel tubes. The working range of the strings is from −100 MPa to 200 MPa, with the accuracy of 0.1 MPa. The dates were collected by a digital auto-reading instrument of GPC-3 type made by Dandong three of instrument group, China. The monitoring point were selected based on the potential risk areas, as shown in Figure 14. The initial readings of the strings were recorded prior to the removal of the soil beneath the STS structure.
Considering Tube 9 and Tube 10 are over the top pilot tunnel (Z3), and are located in one of the most disadvantage region in the pilot construction stage, as shown in Figure 3. So 8 monitoring sensors were installed between Tube 9 and Tube 10 to investigate mechanical behavior of the bolts inside steel tubes when the soil was excavated. Each monitoring section was arranged with two force sensors to verify the correctness of the data. The detail of the installation of the monitoring sensors is shown in Figure 15.

Analysis of Monitoring Results
The stress on the bolts connecting Tube 9 and Tube 10 were monitored during the construction sequence of the station. The stresses recorded in the eight monitoring sensors at different construction stages are shown in Figure 16. The following observations can be made at different construction stages: (1) Pilot tunnels construction: The bolts started going in tension as the excavation started, and the tensile stresses gradually increased along with the soil excavation. When the pilot tunnels were excavated at the distance of 15m, the tensile stress of about 5 kN was recorded in the monitoring point 4#. The tensile stress in the bolts began to increase rapidly when the pilot tunnels were excavated further, and the peak value of 28.9 kN was recorded when the tunnel excavation reached to its final distance of 43.5 m.
As seen in Figure 16, with the tunnel excavation, the monitoring values increase gradually, and the increment shows up an increasingly large trend. The reasons were as follows: (a) The influence of steel pipe jacking in accuracy in the pre-construction stage. According to the construction records,

Analysis of Monitoring Results
The stress on the bolts connecting Tube 9 and Tube 10 were monitored during the construction sequence of the station. The stresses recorded in the eight monitoring sensors at different construction stages are shown in Figure 16. The following observations can be made at different construction stages: (1) Pilot tunnels construction: The bolts started going in tension as the excavation started, and the tensile stresses gradually increased along with the soil excavation. When the pilot tunnels were excavated at the distance of 15m, the tensile stress of about 5 kN was recorded in the monitoring point 4#. The tensile stress in the bolts began to increase rapidly when the pilot tunnels were excavated further, and the peak value of 28.9 kN was recorded when the tunnel excavation reached to its final distance of 43.5 m.
As seen in Figure 16, with the tunnel excavation, the monitoring values increase gradually, and the increment shows up an increasingly large trend. The reasons were as follows: (a) The influence of steel pipe jacking in accuracy in the pre-construction stage. According to the construction records, the deviation of the early 25-m is smaller than the rest distance in the steel pipe jacking stage, which caused the quality of flange welding in the last half was inferior to the first half. So the main stress was borne by bolts in the last half of the whole support system; (b) The influence of construction progress. The excavation period of the early 25-m is less slow than the rest distance in excavation soil below STS structure, which causes the ground disturbance to be smaller.
(2) Beam-column supporting system construction: As the bottom longitudinal beam, cylindrical steel column and top longitudinal beam located inside the pilot tunnel are built, respectively. Tensile stresses in the bolts gradually reduced at this stage. The reason is that part of the overburden load is carried by the column, which releases some of the stress in the bolts.
(3) Construction of the station hall: The stresses in the bolts decreased during this stage, and even changed from tension to compression. The monitoring location produced negative moment along with soil excavation, and it happened because the stiffness of the STS structure supported by the steel columns was higher than the stiffness of the rest, so the settlement of the structure was larger on the sides compared to the supporting section.
(4) Construction of the platform layer and the later stage: stresses in the bolts remained almost constant during this stage because the whole supporting system had already reached equilibrium.
Appl. Sci. 2017, 7, x 14 of 29 the deviation of the early 25-m is smaller than the rest distance in the steel pipe jacking stage, which caused the quality of flange welding in the last half was inferior to the first half. So the main stress was borne by bolts in the last half of the whole support system; (b) The influence of construction progress. The excavation period of the early 25-m is less slow than the rest distance in excavation soil below STS structure, which causes the ground disturbance to be smaller.
(2) Beam-column supporting system construction: As the bottom longitudinal beam, cylindrical steel column and top longitudinal beam located inside the pilot tunnel are built, respectively. Tensile stresses in the bolts gradually reduced at this stage. The reason is that part of the overburden load is carried by the column, which releases some of the stress in the bolts.
(3) Construction of the station hall: The stresses in the bolts decreased during this stage, and even changed from tension to compression. The monitoring location produced negative moment along with soil excavation, and it happened because the stiffness of the STS structure supported by the steel columns was higher than the stiffness of the rest, so the settlement of the structure was larger on the sides compared to the supporting section.
(4) Construction of the platform layer and the later stage: stresses in the bolts remained almost constant during this stage because the whole supporting system had already reached equilibrium.

Simplification of STS Structure
Considering the complexity of the STS structure, accounting for all its components in a numerical model is extremely complicated and greatly increases computational expense. To avoid this problem, an equivalent structure with a similar behavior is first established for the use in the subsequent numerical simulations.
In order to obtain the mechanical behavior of the STS beam, a full-scale laboratory test is performed by the authors [33]. It was observed that the beam failed in flexure. In the early stage of the load application, the specimen was in an elastic state and no cracking was observed. Cracking started to occur along the interface between the steel tubes and the concrete in the mid-span at the level of about 7% of the measured peak load. As the applied load was increased, more flexural cracks appeared within both the pure bending and shear regions of the beam and the existing cracks were propagated toward the loading points with increasing width. The yielding of the bottom flange and connecting bolts began at the level of about 86% of the recorded peak load. During the loading, the crack propagation of the beam occurred slowly due to the constraining effect of the steel tube and flange plate. The flexural failure of the beam was initiated by the yielding of the bottom flange followed by the sudden rupture of the bottom flange and bolts, with crushing of the concrete. The cracks were developed mainly along the interface between the tubes and concrete in the region between the tubes. Several cracks were also formed along the interface between the tubes and the core concrete inside the tubes because of the tensile effect of the flange plate on the tubes.
Based on the full-scale laboratory test a general load-deformation behavior of an STS beam can be established. The comparison results between the laboratory and numerical calculation were shown in Figure 17, the trend of laboratory test was familiar with numerical calculation. It indicates that element, interaction and stress-strain applied the model were valid.

Simplification of STS Structure
Considering the complexity of the STS structure, accounting for all its components in a numerical model is extremely complicated and greatly increases computational expense. To avoid this problem, an equivalent structure with a similar behavior is first established for the use in the subsequent numerical simulations.
In order to obtain the mechanical behavior of the STS beam, a full-scale laboratory test is performed by the authors [33]. It was observed that the beam failed in flexure. In the early stage of the load application, the specimen was in an elastic state and no cracking was observed. Cracking started to occur along the interface between the steel tubes and the concrete in the mid-span at the level of about 7% of the measured peak load. As the applied load was increased, more flexural cracks appeared within both the pure bending and shear regions of the beam and the existing cracks were propagated toward the loading points with increasing width. The yielding of the bottom flange and connecting bolts began at the level of about 86% of the recorded peak load. During the loading, the crack propagation of the beam occurred slowly due to the constraining effect of the steel tube and flange plate. The flexural failure of the beam was initiated by the yielding of the bottom flange followed by the sudden rupture of the bottom flange and bolts, with crushing of the concrete. The cracks were developed mainly along the interface between the tubes and concrete in the region between the tubes. Several cracks were also formed along the interface between the tubes and the core concrete inside the tubes because of the tensile effect of the flange plate on the tubes.
Based on the full-scale laboratory test a general load-deformation behavior of an STS beam can be established. The comparison results between the laboratory and numerical calculation were shown in Figure 17, the trend of laboratory test was familiar with numerical calculation. It indicates that element, interaction and stress-strain applied the model were valid. According to the construction schedule, the maximum free surface of the subway station are 7 steel tubes, considering the influence size-effect on calculation results. The model was established to study the effect welding of flange on bearing capacity (see Figure 18). Due to the STS structure being in the elastics stage in the construction of subway station, and according to the deflection formula of simply-supported beam [34], the equivalent stiffness of welding of bottom flange and no welding of bottom flange were 1.91 × 108 N·m 2 and 1.05 × 108 N·m 2 , respectively (see Figure 19).  Simplification model of the STS structure was divided into two materials (material 1 and material 2) to simulate the weak part of adjacent steel tubes, as shown in Figure 20. Not only can it ensure the whole stiffness of STS structure, but also retain the characteristics of STS structure. Considering the longitudinal direction of STS structure has a high stiffness, the longitudinal stiffness of material 1 is equivalent to longitudinal stiffness of concrete filled steel tube, and the stiffness is 92 Figure 17.
Comparison of load-deformation curves between laboratory test results and numerical calculation.
According to the construction schedule, the maximum free surface of the subway station are 7 steel tubes, considering the influence size-effect on calculation results. The model was established to study the effect welding of flange on bearing capacity (see Figure 18). Due to the STS structure being in the elastics stage in the construction of subway station, and according to the deflection formula of simply-supported beam [34], the equivalent stiffness of welding of bottom flange and no welding of bottom flange were 1.91 × 108 N·m 2 and 1.05 × 108 N·m 2 , respectively (see Figure 19). According to the construction schedule, the maximum free surface of the subway station are 7 steel tubes, considering the influence size-effect on calculation results. The model was established to study the effect welding of flange on bearing capacity (see Figure 18). Due to the STS structure being in the elastics stage in the construction of subway station, and according to the deflection formula of simply-supported beam [34], the equivalent stiffness of welding of bottom flange and no welding of bottom flange were 1.91 × 108 N·m 2 and 1.05 × 108 N·m 2 , respectively (see Figure 19).  Simplification model of the STS structure was divided into two materials (material 1 and material 2) to simulate the weak part of adjacent steel tubes, as shown in Figure 20. Not only can it ensure the whole stiffness of STS structure, but also retain the characteristics of STS structure. Considering the longitudinal direction of STS structure has a high stiffness, the longitudinal stiffness of material 1 is equivalent to longitudinal stiffness of concrete filled steel tube, and the stiffness is 92 According to the construction schedule, the maximum free surface of the subway station are 7 steel tubes, considering the influence size-effect on calculation results. The model was established to study the effect welding of flange on bearing capacity (see Figure 18). Due to the STS structure being in the elastics stage in the construction of subway station, and according to the deflection formula of simply-supported beam [34], the equivalent stiffness of welding of bottom flange and no welding of bottom flange were 1.91 × 108 N·m 2 and 1.05 × 108 N·m 2 , respectively (see Figure 19).  Simplification model of the STS structure was divided into two materials (material 1 and material 2) to simulate the weak part of adjacent steel tubes, as shown in Figure 20. Not only can it ensure the whole stiffness of STS structure, but also retain the characteristics of STS structure. Considering the longitudinal direction of STS structure has a high stiffness, the longitudinal stiffness of material 1 is equivalent to longitudinal stiffness of concrete filled steel tube, and the stiffness is 92 Simplification model of the STS structure was divided into two materials (material 1 and material 2) to simulate the weak part of adjacent steel tubes, as shown in Figure 20. Not only can it ensure the whole stiffness of STS structure, but also retain the characteristics of STS structure. Considering the longitudinal direction of STS structure has a high stiffness, the longitudinal stiffness of material 1 is equivalent to longitudinal stiffness of concrete filled steel tube, and the stiffness is 92 GPa according to Formula (5). Considering the effect of the steel tube with the hole on the whole stiffness, the stiffness was reduced to 70%, that is 65 GPa [35].
where, E-the equivalent stiffness, I-the equivalent moment of inertia; E s -the stiffness of steel tube, I s -the moment of inertia of steel tube; E c -the stiffness of concrete, I c -the moment of inertia of concrete.
Appl. Sci. 2017, 7, x 17 of 29 GPa according to Formula (5). Considering the effect of the steel tube with the hole on the whole stiffness, the stiffness was reduced to 70%, that is 65 GPa [35].
where, E-the equivalent stiffness, I-the equivalent moment of inertia; Es-the stiffness of steel tube, Is-the moment of inertia of steel tube; Ec-the stiffness of concrete, Ic-the moment of inertia of concrete. The seven models considered different elastic modulus of material 2 were established to obtain load-deformation curve of the models by exerting concentrated load of 500 kN, and the simplified model occur elastics deformation from Figure 21. The equivalent lateral stiffness of the seven models was calculated by load-deformation curve, and the relationship between elastic modulus of material 2 and the lateral stiffness of the simplified model was shown Figure 22. So when the bottom flange was welded, the elastic modulus of the material 2 is 2.3 GPa, and when the bottom flange was not welded, the elastic modulus of the material 2 is 1.21 GPa.  7.0x10 8 1.4x10 9 2.1x10 9 2.8x10 9 3.5x10 9 3.0x10 7 6.0x10 The seven models considered different elastic modulus of material 2 were established to obtain load-deformation curve of the models by exerting concentrated load of 500 kN, and the simplified model occur elastics deformation from Figure 21. The equivalent lateral stiffness of the seven models was calculated by load-deformation curve, and the relationship between elastic modulus of material 2 and the lateral stiffness of the simplified model was shown Figure 22. So when the bottom flange was welded, the elastic modulus of the material 2 is 2.3 GPa, and when the bottom flange was not welded, the elastic modulus of the material 2 is 1.21 GPa. GPa according to Formula (5). Considering the effect of the steel tube with the hole on the whole stiffness, the stiffness was reduced to 70%, that is 65 GPa [35]. The seven models considered different elastic modulus of material 2 were established to obtain load-deformation curve of the models by exerting concentrated load of 500 kN, and the simplified model occur elastics deformation from Figure 21. The equivalent lateral stiffness of the seven models was calculated by load-deformation curve, and the relationship between elastic modulus of material 2 and the lateral stiffness of the simplified model was shown Figure 22. So when the bottom flange was welded, the elastic modulus of the material 2 is 2.3 GPa, and when the bottom flange was not welded, the elastic modulus of the material 2 is 1.21 GPa.  7.0x10 8 1.4x10 9 2.1x10 9 2.8x10 9 3.5x10 9 3.0x10 7 6.0x10 GPa according to Formula (5). Considering the effect of the steel tube with the hole on the whole stiffness, the stiffness was reduced to 70%, that is 65 GPa [35]. The seven models considered different elastic modulus of material 2 were established to obtain load-deformation curve of the models by exerting concentrated load of 500 kN, and the simplified model occur elastics deformation from Figure 21. The equivalent lateral stiffness of the seven models was calculated by load-deformation curve, and the relationship between elastic modulus of material 2 and the lateral stiffness of the simplified model was shown Figure 22. So when the bottom flange was welded, the elastic modulus of the material 2 is 2.3 GPa, and when the bottom flange was not welded, the elastic modulus of the material 2 is 1.21 GPa.  7.0x10 8 1.4x10 9 2.1x10 9 2.8x10 9 3.5x10 9 3.0x10 7 6.0x10

Model and Material Properties
A large span underground subway station with relatively limited size has three-dimensional characteristics. Therefore, to numerically study the general performance of the suggested support system, a full three-dimensional model of the station construction was established using FLAC3D 5.0 software (2012, Itasca Consulting Group, Minneapolis MN, USA). The dimension of the numerical model was selected 120 m × 39 m × 50 m to include most of the soil affected by the excavation, as shown in Figure 23. The bottom boundary of the model was entirely fixed in the whole directions. The horizontal displacements were fixed in the left and wright sides, while vertical movements were free. The top boundary was free to move in all directions. A surcharge of 20 kPa was applied on the top surface to simulate surface loads [36]. The diameter of the cylindrical steel columns is 0.8 m, the distance between cylindrical steel columns is 8.5 m, the length of the feet-lock bolt is 2.5 m, and the distance between the feet-lock bolts is 1.0 m in the model, in accordance to the real structure.

Model and Material Properties
A large span underground subway station with relatively limited size has three-dimensional characteristics. Therefore, to numerically study the general performance of the suggested support system, a full three-dimensional model of the station construction was established using FLAC3D 5.0 software (2012, Itasca Consulting Group, Minneapolis MN, USA). The dimension of the numerical model was selected 120 m × 39 m × 50 m to include most of the soil affected by the excavation, as shown in Figure 23. The bottom boundary of the model was entirely fixed in the whole directions. The horizontal displacements were fixed in the left and wright sides, while vertical movements were free. The top boundary was free to move in all directions. A surcharge of 20 kPa was applied on the top surface to simulate surface loads [36]. The diameter of the cylindrical steel columns is 0.8 m, the distance between cylindrical steel columns is 8.5 m, the length of the feet-lock bolt is 2.5 m, and the distance between the feet-lock bolts is 1.0 m in the model, in accordance to the real structure. Beam and pile elements are adopted to simulate the temporary steel supports and feet-lock bolt, respectively. Different materials were assigned to the different soil layers and different components of the STS structure as seen in Figure 20b. Line and shell elements were used to simulate the cylindrical steel column and the primary lining, respectively. In total, 315,055 zones and 323,540 nodes were used, meanwhile the maximum element size is 6 m × 4 m × 1 m, and the minimum element size is 0.25 m × 0.25 m × 1 m. For simplicity, the Mohr-Coulomb model was adopted to define the behavior of the natural and improved soil in the numerical analyses, the parameters of soil were shown in Table 2. The cohesions and friction angles of sandy soils were obtained according to the recommended relationship between the strength parameter and DPSH value from the scientific references. Besides, in the Shenyang area, the cohesions (c) are usually negligible based on engineering experience and some scientific references. But for the numerical simulation software, c = 0 may cause some options performing hardly, even calculation non-convergence. Therefore, the cohesions of sandy soils are generally assumed to be 1 kPa in the numerical analysis [37][38][39]. STS structure elements, primary lining, beam and columns are deemed to have a liner-elastic behavior. Material properties used for the numerical simulations are listed in Table 3.
As for the initial condition, the initial stress field of the soil is generated by balancing the gravity of the soil. In the initial stage of the model calculation, and the initial displacement is reset to zero. Meanwhile, based on the construction experience of Shenyang and the relevant research, it is proved that dewatering has little influence on soil deformation. So hydrostatic or geostatic condition is generally ignored.
The interface elements are used to simulate the interaction between soil and structure, and the interface elements are comprised of a series of the triangular elements contained three nodes, as is Beam and pile elements are adopted to simulate the temporary steel supports and feet-lock bolt, respectively. Different materials were assigned to the different soil layers and different components of the STS structure as seen in Figure 20b. Line and shell elements were used to simulate the cylindrical steel column and the primary lining, respectively. In total, 315,055 zones and 323,540 nodes were used, meanwhile the maximum element size is 6 m × 4 m × 1 m, and the minimum element size is 0.25 m × 0.25 m × 1 m. For simplicity, the Mohr-Coulomb model was adopted to define the behavior of the natural and improved soil in the numerical analyses, the parameters of soil were shown in Table 2. The cohesions and friction angles of sandy soils were obtained according to the recommended relationship between the strength parameter and DPSH value from the scientific references. Besides, in the Shenyang area, the cohesions (c) are usually negligible based on engineering experience and some scientific references. But for the numerical simulation software, c = 0 may cause some options performing hardly, even calculation non-convergence. Therefore, the cohesions of sandy soils are generally assumed to be 1 kPa in the numerical analysis [37][38][39]. STS structure elements, primary lining, beam and columns are deemed to have a liner-elastic behavior. Material properties used for the numerical simulations are listed in Table 3.
As for the initial condition, the initial stress field of the soil is generated by balancing the gravity of the soil. In the initial stage of the model calculation, and the initial displacement is reset to zero. Meanwhile, based on the construction experience of Shenyang and the relevant research, it is proved that dewatering has little influence on soil deformation. So hydrostatic or geostatic condition is generally ignored.
The interface elements are used to simulate the interaction between soil and structure, and the interface elements are comprised of a series of the triangular elements contained three nodes, as is shown in Figure 24. The normal and shear constitutive equation of interfaces elements are as follows [40]: where, F n (t+∆t) is the normal at time (t + ∆t); F si (t+∆t) is the shear force vector at time (t + ∆t); µ n is the absolute normal penetration of the interface node into the target face; µ si is the incremental relative shear displacement vector; σ n is the additional normal stress added due to interface stress initialization; σ si is the additional shear stress added due to interface stress initialization; k s is shear stiffness, k n is normal stiffness, and A is representative area associated with the interface node. shown in Figure 24. The normal and shear constitutive equation of interfaces elements are as follows [40]: where, Fn (t+Δt) is the normal at time (t + Δt); Fsi (t+Δt) is the shear force vector at time (t + Δt); μn is the absolute normal penetration of the interface node into the target face; μsi is the incremental relative shear displacement vector; σn is the additional normal stress added due to interface stress initialization; σsi is the additional shear stress added due to interface stress initialization; ks is shear stiffness, kn is normal stiffness, and A is representative area associated with the interface node. The absolute normal penetration of interface surfaces and contact nodes and shear relative velocity are calculated in every step. Subsequently, they are substituted into the constitutive equations of interface elements to obtain normal stress vector and shear stress vector.  The absolute normal penetration of interface surfaces and contact nodes and shear relative velocity are calculated in every step. Subsequently, they are substituted into the constitutive equations of interface elements to obtain normal stress vector and shear stress vector.

Surface Settlement
In accordance with the construction sequence of Northeastern Road subway station, the whole construction simulation was divided into 6 stages, as follows: (1) Stage 1: Construction of the side pilot tunnels Z5 and Z6; (2) Stage 2: Construction of the pilot tunnels Z1 and Z2 and cylindrical steel column C1; The soil excavation was simulated in the subway construction stage, which causes stress redistribution of in the surrounding soil and surface settlement. A key point on the centerline of the cross-section, which is 15 m away from the initial excavation section was selected as the typical monitoring point, and the behavior of surface settlement is studied in the construction procedure. Comparison of surface settlements from in-situ monitoring and the numerical simulations is shown in Figure 24. As can be seen from the figure, the maximum ground settlement recorded in the site is around 12 mm occurred at the end of stage 5 of the construction. The ground uplift of 1.5 mm was observed after the platform layer is constructed mainly due to the heave of the soil under the supporting structure. Figure 25 shows a reasonable agreement between the field data and the numerical simulations given the complexity of the nature of the problem.

Surface Settlement
In accordance with the construction sequence of Northeastern Road subway station, the whole construction simulation was divided into 6 stages, as follows: (1) Stage 1: Construction of the side pilot tunnels Z5 and Z6; (2) Stage 2: Construction of the pilot tunnels Z1 and Z2 and cylindrical steel column C1; The soil excavation was simulated in the subway construction stage, which causes stress redistribution of in the surrounding soil and surface settlement. A key point on the centerline of the cross-section, which is 15 m away from the initial excavation section was selected as the typical monitoring point, and the behavior of surface settlement is studied in the construction procedure. Comparison of surface settlements from in-situ monitoring and the numerical simulations is shown in Figure 24. As can be seen from the figure, the maximum ground settlement recorded in the site is around 12 mm occurred at the end of stage 5 of the construction. The ground uplift of 1.5 mm was observed after the platform layer is constructed mainly due to the heave of the soil under the supporting structure. Figure 25 shows a reasonable agreement between the field data and the numerical simulations given the complexity of the nature of the problem. Based on the finite difference model (FDM) of the construction, the ground settlement profile at the end of each construction stage was obtained, as shown in Figure 26. The following observations can be made: End of stage 1: The maximum monitored ground settlement is 2.2 mm, which accounts for around 20% of the maximum total ground settlement. Soil disturbance due to the construction of the side pilot tunnels are clearly seen.
End of stage 2: The total ground settlement is increased to 6.1 mm accounting for 33% of the maximum ground settlement. The maximum ground settlement is 2.5 mm due to the excavation of pilot tunnel Z2 and is 0.65 mm due to the excavation of pilot tunnel Z1. The maximum ground settlement increases 0.7 mm due to the construction of column C1 and excavation of soils surrounding it. The ground settlement due to the construction of pilot tunnel Z2 accounts for 65% of the ground settlement of this stage, indicating that the excavation of lower pilot tunnels results in a Based on the finite difference model (FDM) of the construction, the ground settlement profile at the end of each construction stage was obtained, as shown in Figure 26. The following observations can be made: End of stage 1: The maximum monitored ground settlement is 2.2 mm, which accounts for around 20% of the maximum total ground settlement. Soil disturbance due to the construction of the side pilot tunnels are clearly seen.
End of stage 2: The total ground settlement is increased to 6.1 mm accounting for 33% of the maximum ground settlement. The maximum ground settlement is 2.5 mm due to the excavation of pilot tunnel Z2 and is 0.65 mm due to the excavation of pilot tunnel Z1. The maximum ground settlement increases 0.7 mm due to the construction of column C1 and excavation of soils surrounding it. The ground settlement due to the construction of pilot tunnel Z2 accounts for 65% of the ground settlement of this stage, indicating that the excavation of lower pilot tunnels results in a considerable portion of the ground settlement when STS pipe-roofing supporting system is used. The initial lining acts as supporting points for pipe-roof structure, which reduces the span of horizontal pipe-roofing system during the excavation of upper pilot tunnels. The large stiffness of the STS support system prevents large ground settlements.
End of stage 3: The cumulative ground settlement is 8.3 mm. The ground settlement increases 2.2 mm at this stage which accounts for around 20% of the maximum ground settlement. The contribution to the final ground settlement reduces in this stage compared to stage 2, most probably because the left side of supporting system is already complete and settled during the construction of stage 3, which reduces the effective lateral span of the pipe-roofing system, therefore increasing its rigidity, resulting in less ground settlement.
End of stage 4: The cumulative ground settlement is 11 mm. The ground settlement increases 2.7 mm in this stage, which accounts for 23% of the maximum ground settlement. The vertical loading is transferred to the beam and pile system due to the removal of the linings on the pilot tunnels, which is the cause of the settlement in this stage.
End of stage 5: The cumulative ground settlement is 11.6 mm, and the ground settlement increases 0.52 mm which accounts for 5.16% of the maximum ground settlement. The whole force bearing system is completed in this step, thus, the construction of the station hall has a little impact on total ground settlement.
End of stage 6: The ground heaves 1.5 mm because it occurs stress release after soil located in the platform layer is excavated and induces the behave of the supporting system and reflects to surface. considerable portion of the ground settlement when STS pipe-roofing supporting system is used. The initial lining acts as supporting points for pipe-roof structure, which reduces the span of horizontal pipe-roofing system during the excavation of upper pilot tunnels. The large stiffness of the STS support system prevents large ground settlements. End of stage 3: The cumulative ground settlement is 8.3 mm. The ground settlement increases 2.2 mm at this stage which accounts for around 20% of the maximum ground settlement. The contribution to the final ground settlement reduces in this stage compared to stage 2, most probably because the left side of supporting system is already complete and settled during the construction of stage 3, which reduces the effective lateral span of the pipe-roofing system, therefore increasing its rigidity, resulting in less ground settlement.
End of stage 4: The cumulative ground settlement is 11 mm. The ground settlement increases 2.7 mm in this stage, which accounts for 23% of the maximum ground settlement. The vertical loading is transferred to the beam and pile system due to the removal of the linings on the pilot tunnels, which is the cause of the settlement in this stage.
End of stage 5: The cumulative ground settlement is 11.6 mm, and the ground settlement increases 0.52 mm which accounts for 5.16% of the maximum ground settlement. The whole force bearing system is completed in this step, thus, the construction of the station hall has a little impact on total ground settlement.
End of stage 6: The ground heaves 1.5 mm because it occurs stress release after soil located in the platform layer is excavated and induces the behave of the supporting system and reflects to surface.

Deformation of the STS Structure
The stability of the STS structure directly influences the safety of the whole construction. Therefore, deformation of the STS structure in different excavation stages is discussed in the following.

Vertical Deformation of the Horizontal Section of the STS Structure
The deformation profile of the horizontal section of the STS structure is shown in Figure 27. Initial ground deformation is incurred by the excavation of the side tunnels (Z5 and Z6), and the deformation is larger at the ends of the horizontal section of the STS structure. The STS structure above pilot tunnel Z1 shows the largest deformation of about 5.9 mm. This is because the STS structure is unsupported along the whole span prior to the construction of column C1, resulting in a comparably large deflection in the STS structure when tunnel Z1 is completed.
The largest deformation of the STS structure is at the top of pilot tunnel Z3 with a deformation of about 3.3 mm during the construction of pilot tunnel Z3, Z4 and steel column C2. The deflection in

Deformation of the STS Structure
The stability of the STS structure directly influences the safety of the whole construction. Therefore, deformation of the STS structure in different excavation stages is discussed in the following.

Vertical Deformation of the Horizontal Section of the STS Structure
The deformation profile of the horizontal section of the STS structure is shown in Figure 27. Initial ground deformation is incurred by the excavation of the side tunnels (Z5 and Z6), and the deformation is larger at the ends of the horizontal section of the STS structure. The STS structure above pilot tunnel Z1 shows the largest deformation of about 5.9 mm. This is because the STS structure is unsupported along the whole span prior to the construction of column C1, resulting in a comparably large deflection in the STS structure when tunnel Z1 is completed.
The largest deformation of the STS structure is at the top of pilot tunnel Z3 with a deformation of about 3.3 mm during the construction of pilot tunnel Z3, Z4 and steel column C2. The deflection in the STS structure after stage 3 of the construction is not as pronounced as stage 2, given that the structure is now partly supported by column C1, as discussed previously. The largest deflection in the STS structure, which is 12.7 mm, happens when the soil is excavated, and the station hall is constructed.
Appl. Sci. 2017, 7, x 22 of 29 the STS structure after stage 3 of the construction is not as pronounced as stage 2, given that the structure is now partly supported by column C1, as discussed previously. The largest deflection in the STS structure, which is 12.7 mm, happens when the soil is excavated, and the station hall is constructed. Horizontal Deformation of the Vertical Section of the STS Structure and the Side Walls Figure 28 shows the horizontal deformation of the vertical section of the STS structure and the side walls at y = 15 m. The deformation profiles for both the right and left sections are almost identical, however, the right section shows slightly more deformation due to unsymmetrical nature of the construction sequence. The largest deformation of around 5 mm appears at the intersection of the right vertical section of the STS structure and the right-side wall, where a weak joint may have resulted in the large deformation. The largest horizontal deformation of the bottom of the side walls is 2.4 mm. After construction of the station hall, the largest deformation of STS structure is 3.4 mm, which occurs at the temporary braces. The three steps excavation is applied during the construction of the station hall, and when the first step is excavated the temporary braces were not installed in time, which may cause a large releasing of stress.

Horizontal Deformation of Bridge Piles
As shown in Figure 29, the horizontal deformation profile for all four bridge piles are almost the same, which indicates that it is not sensitive for time effect of construction compared to surface Horizontal Deformation of the Vertical Section of the STS Structure and the Side Walls Figure 28 shows the horizontal deformation of the vertical section of the STS structure and the side walls at y = 15 m. The deformation profiles for both the right and left sections are almost identical, however, the right section shows slightly more deformation due to unsymmetrical nature of the construction sequence. The largest deformation of around 5 mm appears at the intersection of the right vertical section of the STS structure and the right-side wall, where a weak joint may have resulted in the large deformation. The largest horizontal deformation of the bottom of the side walls is 2.4 mm. After construction of the station hall, the largest deformation of STS structure is 3.4 mm, which occurs at the temporary braces. The three steps excavation is applied during the construction of the station hall, and when the first step is excavated the temporary braces were not installed in time, which may cause a large releasing of stress. the STS structure after stage 3 of the construction is not as pronounced as stage 2, given that the structure is now partly supported by column C1, as discussed previously. The largest deflection in the STS structure, which is 12.7 mm, happens when the soil is excavated, and the station hall is constructed. Horizontal Deformation of the Vertical Section of the STS Structure and the Side Walls Figure 28 shows the horizontal deformation of the vertical section of the STS structure and the side walls at y = 15 m. The deformation profiles for both the right and left sections are almost identical, however, the right section shows slightly more deformation due to unsymmetrical nature of the construction sequence. The largest deformation of around 5 mm appears at the intersection of the right vertical section of the STS structure and the right-side wall, where a weak joint may have resulted in the large deformation. The largest horizontal deformation of the bottom of the side walls is 2.4 mm. After construction of the station hall, the largest deformation of STS structure is 3.4 mm, which occurs at the temporary braces. The three steps excavation is applied during the construction of the station hall, and when the first step is excavated the temporary braces were not installed in time, which may cause a large releasing of stress. As shown in Figure 29, the horizontal deformation profile for all four bridge piles are almost the same, which indicates that it is not sensitive for time effect of construction compared to surface

Horizontal Deformation of Bridge Piles
As shown in Figure 29, the horizontal deformation profile for all four bridge piles are almost the same, which indicates that it is not sensitive for time effect of construction compared to surface settlement, and showing a maximum deformation of around 3 mm. The support provide by the vertical sections of the STS system largely decreases the impact of construction on bridge piles. The largest horizontal deformation occurs at the top of the piles (around 3 mm), and the horizontal displacements at the bottom of the piles are almost negligible. The horizontal displacement of the bottom half of the bridge piles mainly occurs in the last two stages of the construction. In their previous study, the authors have shown the importance of bottom flanges welding on the performance of STS structures through extensive experimental work [18,19]. In the construction of the Northwest Street station, the flanges at the free face of the STS structure were welded timely after soil excavation as shown in Figure 30a. In this section, the effect of the welding of the flanges are investigated on the performance of the STS structure in the station construction using the numerical simulations. As discussed before and shown in Figure 30b, the STS structure is simulated as a composite of two different materials. The STS structure with or without welding of the flanges can be simulated by appropriately adjusting the material properties of the composite system as discussed earlier in the paper. The impact of welding of the flanges in the STS system on the ground In their previous study, the authors have shown the importance of bottom flanges welding on the performance of STS structures through extensive experimental work [18,19]. In the construction of the Northwest Street station, the flanges at the free face of the STS structure were welded timely after soil excavation as shown in Figure 30a. In this section, the effect of the welding of the flanges are investigated on the performance of the STS structure in the station construction using the numerical simulations. As discussed before and shown in Figure 30b, the STS structure is simulated as a composite of two different materials. The STS structure with or without welding of the flanges can be simulated by appropriately adjusting the material properties of the composite system as discussed earlier in the paper. The impact of welding of the flanges in the STS system on the ground displacement is studied by simulating different welding conditions for the upper left pilot tunnel (Z1).   The settlement profile for the crown of Z1 and Z2 tunnels, and after the bottom (Z2) is completed, the ground with and without welding of the flanges are shown in Figures 31 and 32 for the section at y = 15 m. The ground settlement above the pilot tunnel Z1 increases around 1.3 mm with the welding of flanges during the excavation, and the ground settlement increases 1.35 mm without welding of the flange, increasing 5.47%. The crown settlement of the pilot tunnel Z1 also shows only a slight decrease (from 2.1 mm to 1.9 mm) when flanges are welded. It can be concluded that the welding of flanges during the excavation of the upper pilot tunnel Z1 has little impact on the settlement of the whole structure.   (1) Stability analysis of excavation face The stability of the excavation face must be considered in tunnel construction. Pilot construction supported by STS structure is different from other arching excavation face. Considering the bottom tunnel pilots (Z1 and Z3) satisfy the symmetry, the tunnel pilot (Z1) is listed as an example to investigate the stability of top pilot tunnel in this section.
The horizontal displacement of the excavation face in the pilot tunnel Z1 is shown in Figure 33. When the step length is 3 m, the maximum horizontal displacement is 3.4 mm which occurs slightly lower than the center of the upper step section. The displacements decrease gradually with the increase of step length. When the step length is larger than 3 m, the step length has a little effect on the displacement of the excavation face. The numerical results indicate that sliding between upper and lower steps may happen when the step length is less than 3 m. The numerical results highlight the importance of adopting an appropriate step length in the stability of the excavation face. Based on the results of the numerical analysis, it seems that one tunnel diameter, which is 3 m in the case of this simulation, is appropriate for this purpose. Although larger step lengths will also have similar face displacements, they may result in larger ground settlement due to the delay they cause in installation of the initial supports of the tunnel. Step location Figure 33. Horizontal displacement at the tunnel face for different step lengths. (1) Stability analysis of excavation face The stability of the excavation face must be considered in tunnel construction. Pilot construction supported by STS structure is different from other arching excavation face. Considering the bottom tunnel pilots (Z1 and Z3) satisfy the symmetry, the tunnel pilot (Z1) is listed as an example to investigate the stability of top pilot tunnel in this section.
The horizontal displacement of the excavation face in the pilot tunnel Z1 is shown in Figure 33. When the step length is 3 m, the maximum horizontal displacement is 3.4 mm which occurs slightly lower than the center of the upper step section. The displacements decrease gradually with the increase of step length. When the step length is larger than 3 m, the step length has a little effect on the displacement of the excavation face. The numerical results indicate that sliding between upper and lower steps may happen when the step length is less than 3 m. The numerical results highlight the importance of adopting an appropriate step length in the stability of the excavation face. Based on the results of the numerical analysis, it seems that one tunnel diameter, which is 3 m in the case of this simulation, is appropriate for this purpose. Although larger step lengths will also have similar face displacements, they may result in larger ground settlement due to the delay they cause in installation of the initial supports of the tunnel. (1) Stability analysis of excavation face The stability of the excavation face must be considered in tunnel construction. Pilot construction supported by STS structure is different from other arching excavation face. Considering the bottom tunnel pilots (Z1 and Z3) satisfy the symmetry, the tunnel pilot (Z1) is listed as an example to investigate the stability of top pilot tunnel in this section.
The horizontal displacement of the excavation face in the pilot tunnel Z1 is shown in Figure 33. When the step length is 3 m, the maximum horizontal displacement is 3.4 mm which occurs slightly lower than the center of the upper step section. The displacements decrease gradually with the increase of step length. When the step length is larger than 3 m, the step length has a little effect on the displacement of the excavation face. The numerical results indicate that sliding between upper and lower steps may happen when the step length is less than 3 m. The numerical results highlight the importance of adopting an appropriate step length in the stability of the excavation face. Based on the results of the numerical analysis, it seems that one tunnel diameter, which is 3 m in the case of this simulation, is appropriate for this purpose. Although larger step lengths will also have similar face displacements, they may result in larger ground settlement due to the delay they cause in installation of the initial supports of the tunnel. Step location Figure 33. Horizontal displacement at the tunnel face for different step lengths. (2) Ground Settlement Monitoring points were installed on the crown of the pilot tunnel Z1 in longitudinal direction at y = 20 m in order to investigate the impact of step length on ground settlement.
In this section, the effect of different step lengths on the ground settlement and crown settlement at cross section with y = 20 m is investigated through the numerical simulations, and the ground settlement profile and crown settlement profile are showed in Figures 34 and 35. (2) Ground Settlement Monitoring points were installed on the crown of the pilot tunnel Z1 in longitudinal direction at y = 20 m in order to investigate the impact of step length on ground settlement.
In this section, the effect of different step lengths on the ground settlement and crown settlement at cross section with y = 20 m is investigated through the numerical simulations, and the ground settlement profile and crown settlement profile are showed in Figures 34 and 35.  As seen from Figures 33 and 34, the step length has great influence on the ground and tunnel crown settlement above the excavation face. The ground settlement increases with the increase of step length. This is because when the step length is large, it is hard to effectively control the settlement because the ring supports can only be installed when both steps of the excavation are completed. The small local heaves observed in the ground settlement profiles are mainly due to the negative moment created at the end of the horizontal section of the STS system due to settlements of the middle section of the STS system. As shown in Figure 35, larger stage length can cause more crown settlement, appropriate step length must be selected so that the support rings can be quickly installed and the ground settlements can be controlled.

Longitudinal Displacement Profile of the Pilot Tunnel
The longitudinal displacement profile (LDT) of the pilot tunnel Z1 is investigated here. The step length is assumed 3 m in the numerical simulations in accordance to what was used in the site.
As for the ground and the crown of the pilot tunnel, the simulation and monitoring value show similar settlement profiles when pilot tunnel Z1 is excavated, and the largest difference of both is within 10%, as shown in Figure 36. The maximum settlement of ground is no more than 1.25 mm, and the crown settlement is less than 2.0 mm. Settlements are noticeable only for sections that are more than 15 m away from the initial face of the tunnel, and the settlement tends to be a constant (2) Ground Settlement Monitoring points were installed on the crown of the pilot tunnel Z1 in longitudinal direction at y = 20 m in order to investigate the impact of step length on ground settlement.
In this section, the effect of different step lengths on the ground settlement and crown settlement at cross section with y = 20 m is investigated through the numerical simulations, and the ground settlement profile and crown settlement profile are showed in Figures 34 and 35.  As seen from Figures 33 and 34, the step length has great influence on the ground and tunnel crown settlement above the excavation face. The ground settlement increases with the increase of step length. This is because when the step length is large, it is hard to effectively control the settlement because the ring supports can only be installed when both steps of the excavation are completed. The small local heaves observed in the ground settlement profiles are mainly due to the negative moment created at the end of the horizontal section of the STS system due to settlements of the middle section of the STS system. As shown in Figure 35, larger stage length can cause more crown settlement, appropriate step length must be selected so that the support rings can be quickly installed and the ground settlements can be controlled.

Longitudinal Displacement Profile of the Pilot Tunnel
The longitudinal displacement profile (LDT) of the pilot tunnel Z1 is investigated here. The step length is assumed 3 m in the numerical simulations in accordance to what was used in the site.
As for the ground and the crown of the pilot tunnel, the simulation and monitoring value show similar settlement profiles when pilot tunnel Z1 is excavated, and the largest difference of both is within 10%, as shown in Figure 36. The maximum settlement of ground is no more than 1.25 mm, and the crown settlement is less than 2.0 mm. Settlements are noticeable only for sections that are more than 15 m away from the initial face of the tunnel, and the settlement tends to be a constant As seen from Figures 33 and 34, the step length has great influence on the ground and tunnel crown settlement above the excavation face. The ground settlement increases with the increase of step length. This is because when the step length is large, it is hard to effectively control the settlement because the ring supports can only be installed when both steps of the excavation are completed. The small local heaves observed in the ground settlement profiles are mainly due to the negative moment created at the end of the horizontal section of the STS system due to settlements of the middle section of the STS system. As shown in Figure 35, larger stage length can cause more crown settlement, appropriate step length must be selected so that the support rings can be quickly installed and the ground settlements can be controlled.

Longitudinal Displacement Profile of the Pilot Tunnel
The longitudinal displacement profile (LDT) of the pilot tunnel Z1 is investigated here. The step length is assumed 3 m in the numerical simulations in accordance to what was used in the site.
As for the ground and the crown of the pilot tunnel, the simulation and monitoring value show similar settlement profiles when pilot tunnel Z1 is excavated, and the largest difference of both is within 10%, as shown in Figure 36. The maximum settlement of ground is no more than 1.25 mm, and the crown settlement is less than 2.0 mm. Settlements are noticeable only for sections that are more than 15 m away from the initial face of the tunnel, and the settlement tends to be a constant when the excavation face exceeds the longitudinal length of 30 m. So, the main influence range is before 5 m and after 10 m of the monitoring section in settlement.

Conclusion
A new pre-supporting system, Steel Tube Slab (STS) for construction of large-span underground spaces is presented. The conclusions are as follows.
(1) Steel Tube Slab (STS) structure used in large span underground excavation with shallow depth can effectively control and reduce the surface settlement and can prevent the impact of soil excavation on adjacent buildings. The STS structure acts as the initial supporting system during the excavation stages and can increase the stiffness of the whole construction system. So, it reduces the final ground surface settlement and the deformation of the existing building.
(2) Choosing the reasonable step length can be an effective way to control the ground surface settlement. It will not only control the slipping of excavation face and longitudinal displacement, but also can close into a ring as soon as possible, which will control the ground settlement effectively.
(3) STS structure is a permanent supporting system, so it is more dependable and serviceable during its service period. STS structure simplifies the force mechanism in the initial supporting system which causes less construction sequences and reduces the engineering cost.
(4) It can be concluded that STS structure can be a very dependable and effectively method in controlling and reducing the surface settlement and the influence of the construction sequence on the existing building, especially as for the large span and ultra-shallow underground structures.

Acknowledgments:
The authors would like to thank Northeastern University and The University of New South Wales in providing the access to the software and other facilities, and thank Arman Khoshghalb for his kind help in the research contents, analyses and the English language and style of the manuscript. The research described in this paper was supported by the Fundamental Research Funds for the Central Universities (Grant No. N160106006), the National Natural Science Foundation of China (Grant No. 51578116) and China Scholarship Council. These supports are gratefully acknowledged. And the data is available for the journal.

Conflicts of Interest:
The authors declare no conflict of interest.

Conclusion
A new pre-supporting system, Steel Tube Slab (STS) for construction of large-span underground spaces is presented. The conclusions are as follows.
(1) Steel Tube Slab (STS) structure used in large span underground excavation with shallow depth can effectively control and reduce the surface settlement and can prevent the impact of soil excavation on adjacent buildings. The STS structure acts as the initial supporting system during the excavation stages and can increase the stiffness of the whole construction system. So, it reduces the final ground surface settlement and the deformation of the existing building.
(2) Choosing the reasonable step length can be an effective way to control the ground surface settlement. It will not only control the slipping of excavation face and longitudinal displacement, but also can close into a ring as soon as possible, which will control the ground settlement effectively.
(3) STS structure is a permanent supporting system, so it is more dependable and serviceable during its service period. STS structure simplifies the force mechanism in the initial supporting system which causes less construction sequences and reduces the engineering cost.
(4) It can be concluded that STS structure can be a very dependable and effectively method in controlling and reducing the surface settlement and the influence of the construction sequence on the existing building, especially as for the large span and ultra-shallow underground structures.