Model Test Study on the Inﬂuence of Ground Surcharges on the Deformation of Shield Tunnels

: Aiming at studying the harm caused by sudden ground loadings on existing shield tunnels, a indoor scaled model test with a geometric similarity ratio of 1:15.5 was adopted. Considering the inﬂuencing factors such as ground loading, burial depth of the shield tunnel, loading position and soil properties, tunnel convergence deformation, tunnel settlement and deep settlement of soil caused by sudden ground loadings are studied. A three-dimensional ﬁnite element simulation is carried out using the Midas software, and deep settlement of soil is calculated by a theoretical method. The purpose of this model test is to further understand the inﬂuence of ground surcharges on shield tunnel deformation. The results show that the greater the ground surcharge, the greater the settlement and vertical convergence deformation of the shield tunnel; The further away from the ground surcharge, the smaller the settlement, vertical convergence deformation and lateral convergence deformation of the tunnel. When the pile load size is constant, the greater the burial depth of the tunnel, the smaller the vertical convergence deformation and settlement of the tunnel; the maximum value of deep settlement of the soil always remains at the closest point to the ground surcharge; compared with the use of dry sand, the vertical convergence deformation and settlement of the tunnel are signiﬁcantly reduced when using wet sand. Both the theoretical calculation results and the numerical simulation results are in good agreement with the indoor model test results.


Introduction
Sudden ground surcharges near tunnels are an important factor threatening the safety of shield tunnel structures. Many tunnel structural damage accidents caused by sudden surcharges have occurred at home and abroad for many years [1][2][3][4]. When a sudden ground surcharge occurs around a tunnel, the ground surcharge will produce additional loads on the shield segment, which will cause a certain deformation (convergence deformation and longitudinal uneven settlement) of the segment. When the deformation is too large, the shield segment will be damaged, the gap between the segments will increase, and the connecting bolts will fail, etc. [5][6][7][8][9]. In order to ensure the safety of existing shield tunnel structures, it is therefore of great significance to study the deformation laws of shield tunnels under sudden ground surcharge conditions. At present, scholars at home and abroad have done a lot of research on the influence of ground surcharges on existing shield tunnels, including full-scale model tests [10], scaled model tests [11][12][13], field measurements [14,15], finite element simulations [14], theoretical calculations [16][17][18][19], etc. Xian et al. [10] used a full-scale model test to study the structural bearing capacity of shield tunnels under ground surcharge conditions; Wu et al. [11] studied the influence of the burial depth of the tunnel and the position of the ground surcharge on the deformation of existing shield tunnels by using the scaled indoor model test; Huang et al. [12] studied the deformation of the tunnel structure and the change of the earth pressure around the tunnel under the action of surface overloads by using a reduced scale indoor model test; Atkinson and Pott [13] studied the stability of tunnel structure underground loads by indoor model tests, and obtained the internal force distribution of tunnel lining structures, but the deformation of these tunnel lining structures could not be determined; Based on an actual project, Ma et al. [14] analyzed the field monitoring data of tunnel convergence deformation, established a numerical calculation model, and simulated the whole process from loading before tunnel construction to unloading after construction; Based on the Tianjin Binhai Metro Line Z2 project, Lu et al. [15] analyzed the influence of different factors on tunnel deformations under surcharge loads; [16], aiming at the tunnel damage caused by eccentric ground surcharges in the tunnel neighborhood, considering the influence of eccentric ground surcharges and the force between shield lining rings, based on the modified conventional method, deduced the calculation formulas of lining confining pressure and internal force (bending moment, axial force and shear force); Yamamoto et al. [17,18] studied and analyzed the stability of one tunnel and two tunnels in cohesive soil under the action of ground surcharges, compared the results with the predicted values and obtained relevant results; According to various engineering cases, Huang et al. [19] established a model to evaluate the resilience of shield tunnels after settlement, and compared it with the measured values. Up to now, when scholars use indoor model tests to study tunnel deformations, the tunnel model is generally simplified, and the connecting bolts between tunnel segments are simply simulated by a reduction method, and no one has studied the degree of influence of convergence deformation, settlement and deep soil settlement of shield tunnels under different ground surcharge conditions [20][21][22][23][24][25]. Therefore, it is necessary to use indoor model test methods for further research.
In this paper, the indoor scaled model test method, supplemented by a threedimensional finite element simulation and theoretical calculation method, are used to study the lateral and vertical convergence deformation, tunnel settlement and deep soil settlement of an existing shield tunnel caused by a sudden ground surcharge on the ground, and to analyze the laws that influence ground surcharge, tunnel burial depth, ground surcharge position and soil properties. At the same time, the finite element simulation results and theoretical calculation results are compared with the indoor model test results.

Introduction of Test Model
In this test, a total of 23 segments of a shield tunnel were selected as the research object. The outer diameter of the shield tunnel is 6.2 m, the ring width of the segment is 1.2 m, and the segment thickness is 0.348 m, with a total length of 27.6 m. The geometric similarity ratio of indoor model test is 1:15.5.
(1) Shield tunnel model In order to truly simulate the actual shield tunnel, the segments of the tunnel model are made of plexiglass with E = 2.06 GPa. Each ring tunnel consists of five segments with 67.5 • and one segment with 22.5 • , and the segments of the tunnel model are connected by bolts. The tunnel model has an outer diameter (D 0 ) of 0.4 m and a length of 1.78 m. See Figure 1 for the overall tunnel model. According to the second similarity theorem, the similarity constants of each quantity are deduced as shown in Table 1. Table 2 lists the geometric parameters terial characteristics of tunnel model. Table 3 contains the geometric parameters terial characteristics of the tunnel connecting bolts [26].   According to the second similarity theorem, the similarity constants of each physical quantity are deduced as shown in Table 1. Table 2 lists the geometric parameters and material characteristics of tunnel model. Table 3 contains the geometric parameters and material characteristics of the tunnel connecting bolts [26].  (2) Indoor model box In this indoor test, the size of model box is 1.8 m × 1.8 m × 1.5 m, the width of model box is 1.8 m, which is about 4.5 D 0 , and the height is 1.5 m, which is 3.75 D 0 . See Figure 2 for the model box used in our indoor model tests.

Test Soil Material
Dried sea sand was used as the soil material in this study. In order to avoid the influ ence of impurities in the sea sand on the test, a sieve with an aperture of 18 meshes (abou 1 mm) is used to screen out the large particles in the sea sand. See Table 4 for the physica and mechanical indexes of the dry sand.

Test Conditions
According to four test control variables, i.e., the size of ground surcharge, the buria depth of the shield tunnel (the distance from the top of tunnel to the ground), the position of ground surcharge and the physical properties of soil, nine sets of model tests were car ried out to study the transverse and vertical convergence deformation, tunnel settlemen and deep settlement of soil caused by sudden ground surcharges. Table 5 lists nine groups of test conditions, in which cases 1-4 mainly measure th transverse convergence deformation of the shield tunnel and deep settlement of soil unde the action of ground surcharges, while working conditions 5-9 mainly measure the verti cal convergence deformation and tunnel settlement of the shield tunnel under the action of ground surcharges. Because the tunnel is small, the convergence deformation can only monitor one lin at a time, so the test parameters of condition 1 and condition 6 are the same. As a standard working condition, two sets of tests are conducted, which respectively test the lateral and

Test Soil Material
Dried sea sand was used as the soil material in this study. In order to avoid the influence of impurities in the sea sand on the test, a sieve with an aperture of 18 meshes (about 1 mm) is used to screen out the large particles in the sea sand. See Table 4 for the physical and mechanical indexes of the dry sand.

Test Conditions
According to four test control variables, i.e., the size of ground surcharge, the burial depth of the shield tunnel (the distance from the top of tunnel to the ground), the position of ground surcharge and the physical properties of soil, nine sets of model tests were carried out to study the transverse and vertical convergence deformation, tunnel settlement and deep settlement of soil caused by sudden ground surcharges. Table 5 lists nine groups of test conditions, in which cases 1-4 mainly measure the transverse convergence deformation of the shield tunnel and deep settlement of soil under the action of ground surcharges, while working conditions 5-9 mainly measure the vertical convergence deformation and tunnel settlement of the shield tunnel under the action of ground surcharges. Because the tunnel is small, the convergence deformation can only monitor one line at a time, so the test parameters of condition 1 and condition 6 are the same. As a standard working condition, two sets of tests are conducted, which respectively test the lateral and vertical convergence of the tunnel, the tunnel settlement and the deep settlement of soil.
The layout of displacement measuring points for vertical convergence deformation and tunnel settlement of shield tunnel is shown in Figures 3 and 4. A total of five groups of measuring points are arranged. The steel rod used to fix the displacement meter in the middle of Figure 4 is fixed on the model box and is immobile, so the difference between the data of measurement point 5 and measurement point 10 is the vertical convergence value of the tunnel, and the data of measurement point 5 is the settlement value of the tunnel. The layout of displacement measuring points for vertical convergence deformation and tunnel settlement of shield tunnel is shown in Figures 3 and 4. A total of five groups of measuring points are arranged. The steel rod used to fix the displacement meter in the middle of Figure 4 is fixed on the model box and is immobile, so the difference between the data of measurement point 5 and measurement point 10 is the vertical convergence value of the tunnel, and the data of measurement point 5 is the settlement value of the tunnel.    The arrangement of displacement measuring points for transverse convergence deformation of the shield tunnel is shown in Figures 5 and 6, with a total of five groups of measuring points. The difference between the data of measuring point 15 and measuring point 20 in Figure 6 is the horizontal convergence value of the tunnel.  The layout of displacement measuring points for vertical convergence deform and tunnel settlement of shield tunnel is shown in Figures 3 and 4. A total of five of measuring points are arranged. The steel rod used to fix the displacement meter middle of Figure 4 is fixed on the model box and is immobile, so the difference be the data of measurement point 5 and measurement point 10 is the vertical conve value of the tunnel, and the data of measurement point 5 is the settlement value tunnel.    The arrangement of displacement measuring points for transverse convergen formation of the shield tunnel is shown in Figures 5 and 6, with a total of five gro measuring points. The difference between the data of measuring point 15 and mea point 20 in Figure 6 is the horizontal convergence value of the tunnel.  The arrangement of displacement measuring points for transverse convergence deformation of the shield tunnel is shown in Figures 5 and 6, with a total of five groups of measuring points. The difference between the data of measuring point 15 and measuring point 20 in Figure 6 is the horizontal convergence value of the tunnel.   The arrangement of displacement measuring points for deep settlement of shown in Figure 7, with a total of nine measuring points, all located 0.2 m belo ground.     The arrangement of displacement measuring points for deep settlement of soil shown in Figure 7, with a total of nine measuring points, all located 0.2 m below t ground.  The arrangement of displacement measuring points for deep settlement of soil is shown in Figure 7, with a total of nine measuring points, all located 0.2 m below the ground. The arrangement of displacement measuring points for deep settlem shown in Figure 7, with a total of nine measuring points, all located 0.2 m ground.

Test Procedure
In this model test, a jack is used to press the pressure-bearing steel plate to simulate the ground surcharge. The forces exerted by the jack are 34.4, 68.8, 103.2 and 137.6 kg, respectively, and the pressure of each stage is applied at an interval of 1 h. The size of the bearing plate is 0.4 m × 0.4 m and the actual ground surcharge simulated in the test are 36, 72, 108 and 144 kPa, respectively. According to the soil gravity of 18 kN/m 3 , the corresponding stacking height is 2, 4 m, 6 and 8 m, respectively, and the corresponding actual stacking range is 6.2 m × 6.2 m. The ground surcharge moves in the direction perpendicular to the shield tunnel, and the center line coincides with the center point of the shield tunnel. See Figure 8 for the model loading process.

Test Procedure
In this model test, a jack is used to press the pressure-bearing steel plate to simulate the ground surcharge. The forces exerted by the jack are 34.4, 68.8, 103.2 and 137.6 kg, respectively, and the pressure of each stage is applied at an interval of 1 h. The size of the bearing plate is 0.4 m × 0.4 m and the actual ground surcharge simulated in the test are 36, 72, 108 and 144 kPa, respectively. According to the soil gravity of 18 kN/m 3 , the corresponding stacking height is 2, 4 m, 6 and 8 m, respectively, and the corresponding actual stacking range is 6.2 m × 6.2 m. The ground surcharge moves in the direction perpendicular to the shield tunnel, and the center line coincides with the center point of the shield tunnel. See Figure 8 for the model loading process. During the model tests, the tunnel model should be hoisted to the designated position first, and then sand should be added step by step. First, 20 cm of sand should be added to make the sand reach the bottom of the tunnel model. After standing for 24 h, the remaining height of sand should be evenly added in two times, that is, the height of the two times is the same, and each time after adding, it needs to stand for 24 h. After each group of tests, the sand in the model box is completely discharged. After that, the next group of model tests were carried out.

Analysis of Indoor Model Test Results
Note: The values in the analysis of test results are actually measured data. During the model tests, the tunnel model should be hoisted to the designated position first, and then sand should be added step by step. First, 20 cm of sand should be added to make the sand reach the bottom of the tunnel model. After standing for 24 h, the remaining height of sand should be evenly added in two times, that is, the height of the two times is the same, and each time after adding, it needs to stand for 24 h. After each group of tests, the sand in the model box is completely discharged. After that, the next group of model tests were carried out.

Analysis of Indoor Model Test Results
Note: The values in the analysis of test results are actually measured data.

Analysis of Test Data under Standard Operating Condition
Let the test parameters of working condition 1 and working condition 6 be the same, and the ground surcharge is just above the tunnel (symmetrical loading) as the standard working conditions. When the burial depth of the tunnel is 0.6 m and the central ground surcharge is 137.6 kg, the tunnel settlement, vertical convergence deformation, transverse convergence deformation, deep soil settlement curve (parallel to the tunnel direction) and deep soil settlement curve (perpendicular to the tunnel direction) caused by ground surcharge are shown in Figures 9-13, respectively, and the negative value in the tunnel settlement curve represents the downward settlement of the tunnel. The negative value in the vertical convergence deformation curve represents the vertical shrinkage of the tunnel; The positive value in the curve of tunnel transverse convergence deformation represents the tunnel transverse enlargement; The negative value in the deep settlement curve (including the direction parallel to and perpendicular to the tunnel) represents the decrease of soil.   It can be seen from Figures 9-11 that: (1) Under the action of ground surch to the increase of vertical earth pressure, the tunnel has a downward settlemen transverse convergence deformation law of the tunnel is basically the same as t deformation of the tunnel, but the numerical value is the opposite, which ind the vertical diameter of the tunnel decreases and the transverse diameter increas ing in a transverse elliptical deformation. (2) The deformation of the tunnel is t just below the ground surcharge. With the distance from the ground surcharge settlement, vertical convergence deformation and transverse convergence defor It can be seen from Figures 9-11 that: (1) Under the action of ground surcharges, due to the increase of vertical earth pressure, the tunnel has a downward settlement, and the transverse convergence deformation law of the tunnel is basically the same as the vertical deformation of the tunnel, but the numerical value is the opposite, which indicates that the vertical diameter of the tunnel decreases and the transverse diameter increases, resulting in a transverse elliptical deformation. (2) The deformation of the tunnel is the largest just below the ground surcharge. With the distance from the ground surcharge center, the settlement, vertical convergence deformation and transverse convergence deformation of the tunnel are all smaller, and the variation law basically accords with a logistic function. The maximum settlement of tunnel is −0.341 mm, which is converted into actual engineering data by the similarity ratio to be −5.286 mm; The maximum vertical convergence deformation of the tunnel is −0.128 mm, which is converted into actual engineering data of −1.984 mm; The maximum transverse convergence deformation value of the tunnel is 0.118 mm, which is converted into the actual engineering data of 1.829 mm.
It can be seen from Figures 12 and 13 that: (1) the deep settlement of soil is related to the distance of ground surcharge, and the deep settlement of soil on both sides of ground surcharge is basically symmetrically distributed, and the maximum deep settlement of soil parallel to the tunnel direction is −1.852 mm, which is converted into actual engineering data by similarity ratio as −28.706 mm; The maximum deep settlement of soil perpendicular to the tunnel direction is −1.980 mm, which is converted into actual engineering data of −30.690 mm. (2) When the distance from the center of the tunnel is 0.2-0.4 m, the deep settlement of soil decreases by 43.5% parallel to the tunnel and 40.0% perpendicular to the tunnel.  It can be seen from Figures 14 and 15 that the settlement and vertical conve deformation of the tunnel basically increase linearly with the increase of groun It can be seen from Figures 14 and 15 that the settlement and vertical convergence deformation of the tunnel basically increase linearly with the increase of ground surcharge. The reason is that the ground surcharge value is not large, the convergence deformation and settlement value of the tunnel do not change much under this working condition, and the properties of the tunnel and sand do not change much before and after loading, so the settlement and vertical convergence deformation of the tunnel basically increase linearly during loading.  It can be seen from Figures 14 and 15 that the settlement and vertical conve deformation of the tunnel basically increase linearly with the increase of groun charge. The reason is that the ground surcharge value is not large, the convergence mation and settlement value of the tunnel do not change much under this workin dition, and the properties of the tunnel and sand do not change much before an loading, so the settlement and vertical convergence deformation of the tunnel ba increase linearly during loading.  As can be seen in Figures 16 and 17, the magnitude of the stacking load remai stant and as the tunnel burial depth increases and there is a significant reduction vertical convergence deformation and settlement of the tunnel. Taking the tunnel se at the center of the tunnel as an example, when the burial depth of the tunnel c from 0.5 m to 0.8 m, the vertical convergence deformation decreases from 0.145 0.094 mm, with a decrease of 35.2%; Settlement decreased from 0.392 mm to 0.24 with a decrease of 38.0%. Therefore, the influence of ground surcharge on the v convergence deformation and settlement of tunnel decreases obviously with the in of buried depth. The test results are basically consistent with the test results in ref [11], but the change trend of the test results is more obvious, which may be related As can be seen in Figures 16 and 17, the magnitude of the stacking load remains constant and as the tunnel burial depth increases and there is a significant reduction in the vertical convergence deformation and settlement of the tunnel. Taking the tunnel segment at the center of the tunnel as an example, when the burial depth of the tunnel changes from 0.5 m to 0.8 m, the vertical convergence deformation decreases from 0.145 mm to 0.094 mm, with a decrease of 35.2%; Settlement decreased from 0.392 mm to 0.243 mm, with a decrease of 38.0%. Therefore, the influence of ground surcharge on the vertical convergence deformation and settlement of tunnel decreases obviously with the increase of buried depth. The test results are basically consistent with the test results in reference [11], but the change trend of the test results is more obvious, which may be related to the different model tunnels used in the test. Wu et al. [11] directly used an iron sheet tube as their tunnel model, while each ring of the tunnel model used in this test is composed of five 67.5 • and one 22.5 • segment, which is more in line with the actual working conditions.

The Influence of Ground Surcharge Position on Tunnel Transverse Convergence Deformation and Deep Settlement of Soil (1) Analysis of model test results
According to working conditions 1, 2, 3 and 4, when the burial depth of the tunnel is 0.6 m, the ground surcharge is 137.6 kg, and the eccentric distance of ground surcharge is 0, 0.2, 0.4 and 0.6 m, respectively, the transverse convergence deformation of the center point of the tunnel, the transverse convergence deformation of each measuring point of the tunnel segment, the deep settlement of soil parallel to and perpendicular to the tunnel direction are shown in Figures 18-21. with a decrease of 38.0%. Therefore, the influence of ground surcharge on the v convergence deformation and settlement of tunnel decreases obviously with the in of buried depth. The test results are basically consistent with the test results in re [11], but the change trend of the test results is more obvious, which may be related different model tunnels used in the test. Wu et al. [11] directly used an iron sheet their tunnel model, while each ring of the tunnel model used in this test is comp five 67.5° and one 22.5° segment, which is more in line with the actual working cond

The Influence of Ground Surcharge Position on Tunnel Transverse Converge Deformation and Deep Settlement of Soil (1) Analysis of model test results
According to working conditions 1, 2, 3 and 4, when the burial depth of the tu 0.6 m, the ground surcharge is 137.6 kg, and the eccentric distance of ground surch 0, 0.2, 0.4 and 0.6 m, respectively, the transverse convergence deformation of the point of the tunnel, the transverse convergence deformation of each measuring p the tunnel segment, the deep settlement of soil parallel to and perpendicular to the direction are shown in Figures 18-21. It can be seen from Figure 18 that the transverse convergence deformation of t nel shows a linear downward trend with the increase of the eccentricity distance pile load. This is because the distance from the pile load to the measuring point in during the drift of the pile load, and the influence of the pile load on the tunnel be smaller and smaller.
It can be seen from Figure 19 that near the center of the tunnel, the surcharge tricity has a great influence on the transverse convergence deformation of each mea point of the tunnel segment. At the distance away from the center of the tunnel, the verse convergence deformation of the tunnel segment is very close at different sur eccentricity. When the surcharge eccentricity is 0-0.2 m, the surcharge eccentric tance has the greatest influence on the transverse convergence deformation of each uring point of the tunnel segment, and then the influence decreases.  It can be seen from Figure 21 that with the movement of the pile load, the dee tlement curve of the soil body approximately produces the whole translation in the direction, and the maximum value is always kept closest to the pile load center, an settlement on both sides of the pile load center is symmetrically distributed. (2) Verification of the Boussinesq solution According to the calculation formula of vertical displacement of soil calcula the Boussinesq solution, combined with uniformly distributed ground load, the ca tion formula of deep settlement of soil is deduced as follows: where: is vertical settlement; μ is poisson's ratio of soil, taking 0.35; E is the modulus of soil, taking 19 × 10 6 pa; P is the vertical force of ground surcharge on so Using Equation (1)  It can be seen from Figure 21 that with the movement of the pile load, the dee tlement curve of the soil body approximately produces the whole translation in the direction, and the maximum value is always kept closest to the pile load center, a settlement on both sides of the pile load center is symmetrically distributed. (2) Verification of the Boussinesq solution According to the calculation formula of vertical displacement of soil calcula the Boussinesq solution, combined with uniformly distributed ground load, the ca tion formula of deep settlement of soil is deduced as follows: where: is vertical settlement; μ is poisson's ratio of soil, taking 0.35; E is the modulus of soil, taking 19 × 10 6 pa; P is the vertical force of ground surcharge on so Using Equation (  It can be seen from Figure 18 that the transverse convergence deformation of the tunnel shows a linear downward trend with the increase of the eccentricity distance of the pile load. This is because the distance from the pile load to the measuring point increases during the drift of the pile load, and the influence of the pile load on the tunnel becomes smaller and smaller.
It can be seen from Figure 19 that near the center of the tunnel, the surcharge eccentricity has a great influence on the transverse convergence deformation of each measuring point of the tunnel segment. At the distance away from the center of the tunnel, the transverse convergence deformation of the tunnel segment is very close at different surcharge eccentricity. When the surcharge eccentricity is 0-0.2 m, the surcharge eccentricity distance has the greatest influence on the transverse convergence deformation of each measuring point of the tunnel segment, and then the influence decreases.
It can be seen from Figure 20 that deep soil settlement occurs under the action of ground surcharge, and the deep settlement of soil decreases with the increase of ground surcharge eccentricity. The initial stage of ground surcharge eccentricity has a great influence on the change value of deep settlement of soil, and with the increase of ground surcharge eccentricity, it has less and less influence on the change value of deep settlement of soil, but the curves of deep settlement of soil caused by each eccentric ground surcharge are similar in shape; The settlement of soil below the pile load is the largest, and the relationship between the settlement of soil at different measuring points presents a semi-normal distribution, which is related to the additional stress distribution of soil caused by the pile load.
It can be seen from Figure 21 that with the movement of the pile load, the deep settlement curve of the soil body approximately produces the whole translation in the same direction, and the maximum value is always kept closest to the pile load center, and the settlement on both sides of the pile load center is symmetrically distributed.
(2) Verification of the Boussinesq solution According to the calculation formula of vertical displacement of soil calculated by the Boussinesq solution, combined with uniformly distributed ground load, the calculation formula of deep settlement of soil is deduced as follows: where: ω is vertical settlement; µ is poisson's ratio of soil, taking 0.35; E is the elastic modulus of soil, taking 19 × 10 6 pa; P is the vertical force of ground surcharge on soil.
Using Equation (1) Comparing Figure 20 with Figure 22 and Figure 21 with Figure 23, it can be f that the change trend of the model test results is quite consistent with the calculatio sults of the Boussinesqs solution, which indicates that the deep settlement results o Comparing Figure 20 with Figure 22 and Figure 21 with Figure 23, it can be that the change trend of the model test results is quite consistent with the calculati sults of the Boussinesqs solution, which indicates that the deep settlement results obtained from the test are reliable. Comparing Figure 20 with Figures 21 and 22 with Figure 23, it can be found that the change trend of the model test results is quite consistent with the calculation results of the Boussinesqs solution, which indicates that the deep settlement results of soil obtained from the test are reliable.
Generally, the Poisson's ratio µ of sandy soil is taken in the range of 0.3-0.35, and the elastic modulus E is taken in the range of 10-25 × 10 6 pa. Taking measurement point 27 as an example, the interaction effect of Poisson's ratio µ and elastic modulus E on settlement ω is shown in Figure 24. Comparing Figure 20 with Figure 22 and Figure 21 with Figure 23, it can be fou that the change trend of the model test results is quite consistent with the calculation sults of the Boussinesqs solution, which indicates that the deep settlement results of obtained from the test are reliable.
Generally, the Poisson's ratio μ of sandy soil is taken in the range of 0.3-0.35, and elastic modulus E is taken in the range of 10-25 × 10 6 pa. Taking measurement point 27 an example, the interaction effect of Poisson's ratio μ and elastic modulus E on settlem ω is shown in Figure 24.  From Figure 24, it can be seen that with the increase of E, the deep settlement of the soil decreases significantly; with the increase of µ, the deep settlement of the soil increases; compared with both, the effect of E on the deep settlement of the soil is greater than that of µ.

Influence of Different Soil Properties on Vertical Convergence Deformation and Settlement of Tunnel
According to working conditions 6 and 9, when the buried depth of the tunnel is 0.6 m, the influence of different soil quality on the vertical convergence deformation and settlement of the tunnel is shown in Figures 25 and 26. The variation of stress σ with soil layer thickness Z underground surcharge is shown in Figure 27, and the stress σ isolines of different soils are shown in Figure 28.
Symmetry 2021, 13, x FOR PEER REVIEW From Figure 24, it can be seen that with the increase of E, the deep settlement soil decreases significantly; with the increase of μ, the deep settlement of the s creases; compared with both, the effect of E on the deep settlement of the soil is than that of μ.

Influence of Different Soil Properties on Vertical Convergence Deformation a Settlement of Tunnel
According to working conditions 6 and 9, when the buried depth of the tunne m, the influence of different soil quality on the vertical convergence deformation a tlement of the tunnel is shown in Figures 25 and 26. The variation of stress w layer thickness Z underground surcharge is shown in Figure 27, and the stress i of different soils are shown in Figure 28.  It can be seen from Figures 25 and 26 that compared with dry sand, when the test soil used is wet sand, the vertical convergence deformation and settlement of the tunnel are obviously reduced. Taking the vertical deformation at the center of the tunnel as an example, when the test soil is dry sand, the vertical convergence deformation of the tunnel is 0.128 mm, and when the test sand is wet sand, the vertical convergence deformation of the tunnel is 0.053 mm, with a decrease of 58.6%.The tunnel settlement decreased from It can be seen from Figures 25 and 26 that compared with dry sand, when the test so used is wet sand, the vertical convergence deformation and settlement of the tunnel a obviously reduced. Taking the vertical deformation at the center of the tunnel as an e ample, when the test soil is dry sand, the vertical convergence deformation of the tunn It can be seen from Figures 25 and 26 that compared with dry sand, when the test soil used is wet sand, the vertical convergence deformation and settlement of the tunnel are obviously reduced. Taking the vertical deformation at the center of the tunnel as an example, when the test soil is dry sand, the vertical convergence deformation of the tunnel is 0.128 mm, and when the test sand is wet sand, the vertical convergence deformation of the tunnel is 0.053 mm, with a decrease of 58.6%. The tunnel settlement decreased from 0.512 mm to 0.197 mm, with a decrease of 61.5%. The reasons for the above changes may be as follows: (1) When the test soil is wet sand, the sand in the model box is surrounded by an impermeable boundary, and there is excess pore water pressure in the sand, and there is still excess pore water pressure when the soil is damaged. When the excess pore water pressure rises in case of sudden ground surcharge, the effective stress of the soil decreases, and the sand is in a flowing state, which leads to the reduction of additional load on the tunnel, resulting in the reduction of settlement and convergence deformation of the tunnel; (2) As shown in Figures 27 and 28, under the action of ground surcharge, due to the different nature of the soil, the stress σ contours are flat when the soil is wet sand, and the stress σ decreases rapidly with the increase of depth Z. When the soil is dry sand, the stress σ contours are relatively long and lean, and the decrease of stress σ with the increase of depth Z is relatively small, which leads to the settlement and deformation of the tunnel is significantly smaller when the soil is wet sand.

Establishment of 3D Finite Element Model
In this paper, the Midas software is used to establish a finite element three-dimensional numerical model of the influence of ground surcharge on the deformation of adjacent existing shield tunnels. The tunnel segment and soil in the three-dimensional model adopt three-dimensional solid, in which the segment model adopts isotropic elastic model and the soil model adopts isotropic Mohr-Coulomb model. In order to reduce the influence of boundary effect on the calculation results, the horizontal plane size of the model is selected as 40 m × 57.6 m. The depth dimension of the model varies with the tunnel depth, ranging from 28.45 m to 33.1 m. The upper part of the model simulates the free boundary of the surface, the lower part imposes normal and tangential constraints, and the front, back, left and right boundaries impose normal constraints. The working conditions of three-dimensional finite element numerical model are shown in Table 6, and the finite element model with buried depth of 9.3 m is shown in Figure 29.

Model Calculation Parameters
The sand model parameters are weight 16.2 kN/m 3 , Young's modulus 19 MPa, Poisson's ratio 0.29, and internal friction angle 29 • . The parameters of shield segment are set according to the geometric parameters of tunnel prototype in Table 2

Model Calculation Parameters
The sand model parameters are weight 16.2 kN/m 3 , Young's modulus 19 MPa, son's ratio 0.29, and internal friction angle 29°. The parameters of shield segment ar according to the geometric parameters of tunnel prototype in Table 2 except elastic m ulus. Considering the influence of actual segment joints and bolts, the elastic modul segment is reduced by 90% of the values in the table.

Comparative Analysis of Simulation Results and Test Results
The overall settlement trend of the finite element tunnel model is shown in Figure 30, which is basically consistent with the test results.
Symmetry 2021, 13, x FOR PEER REVIEW

Comparative Analysis of Simulation Results and Test Results
The overall settlement trend of the finite element tunnel model is shown i which is basically consistent with the test results. The data measured in laboratory tests are converted into actual engineer using the similarity ratio, and the obtained data are compared with the finite e ulation results. Although there is a certain gap between the two results, they a of the same order of magnitude. Because the finite element method is often us itative research, a dimensionless comparative study is carried out in this pape See Figures 31 and 32   The data measured in laboratory tests are converted into actual engineering data by using the similarity ratio, and the obtained data are compared with the finite element simulation results. Although there is a certain gap between the two results, they are certainly of the same order of magnitude. Because the finite element method is often used for qualitative research, a dimensionless comparative study is carried out in this paper.
See       It can be seen from Figure 31 that the tunnel settlement test results and element simulation results are very consistent with each other under differe surcharge conditions applied at the center of the tunnel.  It can be seen from Figure 31 that the tunnel settlement test results and element simulation results are very consistent with each other under differe surcharge conditions applied at the center of the tunnel.
It can be seen from Figures 32, 33 and 35 that the test results of vertical co It can be seen from Figures 34 and 36 that the deep settlement of soil parallel to the tunnel direction and the vertical convergence deformation test of tunnel under different burial depths are in good agreement with the finite element simulation, but there is a certain difference between the ratio of test results and the ratio of finite element simulation. Possible reasons are as follows: (1) Due to the site, similarity ratio and economic reasons, the size of the model box should not be too large in the indoor test, which leads to the influence of the boundary effect of the model box on the test, and some errors will be produced when the model test results are placed directly according to the similarity ratio.
(2) Although the sand used in the test has been static, the static time is short and the sand is still loose. (3) In the finite element simulation, the segments are continuous without considering the bolts and joints, which is different from the tunnel model used in the indoor model test.

Conclusions
(1) The maximum values of tunnel settlement, vertical convergence deformation, transverse convergence deformation and deep settlement of soil all occur at the nearest place to the center of ground surcharge, which is the weakest point of the tunnel under ground surcharge conditions. (2) When the ground surcharge is the same, the vertical convergence deformation and settlement of the shield tunnel are obviously reduced with the increase of the burial depth of the shield tunnel, which indicates that increasing the burial depth of the shield tunnel can reduce the harm caused by ground surcharges. (3) Compared with the dry sand, when the test soil is wet sand, the vertical convergence deformation and settlement of the tunnel are significantly reduced, indicating that the nature of the soil has a great impact on the vertical convergence deformation and settlement of the tunnel, and the wet sand is more resistant to the vertical convergence deformation and settlement of the tunnel caused by the ground surcharge than the dry sand. (4) In calculating the deep settlement of soil caused by surface loading using the Boussinesq solution, the Poisson's ratio µ and the modulus of elasticity E of the soil have a greater influence on the calculation results of the deep settlement of the soil, and the influence of E on the deep settlement of the soil is greater than that of when comparing the two. (5) The finite element method can simulate well the variation laws of transverse and vertical convergence deformation, tunnel settlement and deep settlement of soil caused by sudden surcharges on the ground. Using less model tests to determine the model and parameters of finite element, and then using finite element to carry out a large number of simulations can improve the research efficiency while ensuring the research accuracy.
Horizontal movement of tunnel was not tested in this model test. There are few kinds of soil used in studying the influence of different soil properties on the vertical convergence deformation and settlement of tunnels. These aspects still need further improvement in future experiments and research.