Parameter Sensitivity of Shallow-Bias Tunnel with a Clear Distance Located in Rock

In order to obtain the seismic internal force response laws of a shallow-bias tunnel with a small clear distance, the reliability of the numerical simulation is verified by the shaking table model test. 0e parameter sensitivity of the tunnel is studied by using MIDAS-NX finite element software. 0e effects of seismic wave peak (0.1 g, 0.2 g, 0.3 g, 0.4 g, 0.5 g, and 0.6 g), existing slope angle (30°, 45°, 60°, and 90°), clear distance (1.0D, 1.5D, 2.0D, and 3.0D), and excitation mode (X direction, Z direction, XY direction, and XYZ direction) on the internal force response law of the tunnel are studied, respectively. 0e results show that (1) the shear force gradually increases with the increasing of seismic peak. 0e amplification is different with different measuring points. (2) Under different existing slope-angle conditions, the variation trend of shear force of the tunnel is similar, but the shear force is different.0e existing slope has significant effect on the shear force response of the tunnel, and the degree is different with different slope angles. (3) Under the conditions of 1.5D and 2.0D, the shear force response of the tunnel is stronger, but the response of other conditions is relatively weak. 0e tunnel with 1.5D to 2.0D clear distance should be avoided. Different excitation modes have a significant effect on the shear force response of the tunnel. (4) Under the same excitation mode, the different excitation directions also have a significant effect on the shear force response. (5) 0e shear force response of the tunnel crosssection shows nonlinear variation trend. 0e shear force response is strongest at the arch shoulder and arch foot of the tunnel. 0e research results provide a useful reference for the design of antishock and vibration resistance of the tunnel.


Introduction
Tunnels have been important components of the transportation in mountain areas, which is beneficial to optimize the direction of the line and save construction cost.Previously, many scholars believed that the earthquake-resistant behavior of tunnels were stronger than that of aboveground buildings.After the earthquake observation, however, it is found that under the action of strong ground motion, tunnels may be subjected to extensive deformation or even collapse [1][2][3][4].e surrounding rock has various pattern fractures, and a series of morphological changes occur in an earthquake [5][6][7][8].Moreover, the dynamic response of tunnel structure is highly distinct from that of superstructure [9][10][11].
e study, currently, on seismic force response of the new type tunnel (shallow-bias tunnel with a small distance) has a few been studied by scholars.e seismic force response of the new type tunnel is significantly different from that of the tunnel conventional type.e numerical methods (e.g., FEM, FDM, and BEM), as the most popular approaches, have been used to study the dynamic response of tunnels.A few experimental studies have been reported the acceleration response of the tunnel.Based on the comparison between the shaking table test and the numerical method, the results obtained by the numerical method were credible and reasonable. is study aims to perform a systematic experimental work followed by means of an extended numerical parametric study to elaborate some dimensions of the problem.e effects of various parameters, including seismic wave peak (0.1 g, 0.2 g, 0.3 g, 0.4 g, 0.5 g, and 0.6 g), seismic wave excitation mode (X direction, Z direction, XZ direction, and XYZ direction), clear distance (1.0 D, 1.5 D, 2.0 D, and 3.0 D), and existing slope angle (30 °, 45 °, 60 °, and 90 °) on the seismic force response, were investigated in the present study.

Numerical Simulation
2.1.Eigenvalue Analysis.In order to calculate the mode shape and the natural period of undamped free vibration in MIDAS-NX, the characteristic equation is as follows: where K is the stiffness matrix of the structure, M is the mass matrix of the structure, ω 2 n is the eigenvalue of the n mode, and Φ n is the mode vector of the n mode.
Eigenvalue analysis, also known as free vibration analysis, is used to analyze the inherent dynamic characteristics of structures.By eigenvalue analysis, the main dynamic characteristics, such as mode shape, natural vibration period, and vibration mode participation coefficient, are determined by the structure quality and stiffness.
In the calculation, the load and damping in the equation of motion of the SDOF system are assumed to be 0. en the equation degenerates into free vibration equation, which is as follows: where u is the displacement induced by vibration.If u � A cos ωt (A is related to the initial displacement constant), (2) can be changed into the following formula: e conditions of (2) are that the item in parentheses is 0. So the eigenvalues are as follows: where ω 2 is the eigenvalue, ω is the proper circular frequencies, T is the natural period, and f is the natural frequency.e vibration mode participation coefficient is calculated by the ratio of response between the mode shape and all involved vibration modes.e calculation formula is as follows: where M is the order of the mode shapes, τ m is the mode participation coefficient, M i is the quality of i nodes, and φ im is the m order vibration vector of i node position.
In order to ensure that the earthquake mainly contains enough vibration, modal effective mass M is greater than the sum of all the effective mass of 90% general provisions in the seismic design code.
where M m is the effective mass of various modes.

Computation Damping.
e Rayleigh damping type is used in the numerical simulation.In order to reduce the uncertainty of stiffness damping in the high modes, the sum of both mass proportional damping and stiffness damping is used as the damping matrix.
e proportional damping matrix C adopts the formula proposed by Caughey: where j and N are the nodal degrees of freedom.M −1 K can be calculated by the free vibration formula of the undamped system: Substituting formula (9) in (8), the following formula can be obtained: In M −1 K � ω 2 , ω 2 is equivalent to the number of modes, and it is expressed by ω 2 s .Substituting M −1 K in formula (7) and multiplying (7) by u s   T and u s on both sides, the following formula can be obtained: e damping constant h s of the s order vibration mode is expressed by the following formulas (12) and (13): e damping constant of the mass ratio and stiffness type is shown in the following formulas ( 14) and (15), respectively.

2
Advances in Civil Engineering e Rayleigh damping matrix is shown in the following formula (16).

Mass proportion type:
Stiffness proportional type: Rayleigh type:

Time History Analysis.
e structural dynamic time history analysis refers to the process of calculating the structural response (displacement, velocity, internal force, etc.) at any moment, when the structure is subjected to dynamic loads.e dynamic equilibrium equation is used in the dynamic time history analysis of MIDAS-NX, which is as follows: where [M] is the mass matrix, [C] is the damping matrix, [K] is the stiffness matrix, p(t) is the dynamic loading, € u(t) is the relative acceleration, _ u(t) is the relative velocity, and u(t) is the relative displacement.
e modal superposition method for structural dynamic analysis is used to carry out dynamic time history analysis in MIDAS-NX.e mode superposition method means that the displacement of the structure is solved by a linear combination of orthogonal displacement vectors. is method is better for linear dynamic analysis of large structures.
e premise of using this method is that the damping matrix can be represented by a linear combination of the mass matrix and the stiffness matrix.
e mode superposition method is one of the most widely used methods of structural analysis program.But in nonlinear dynamic analysis, this method has some limitations.In order to make up for this shortcoming, the nonlinear characteristics of stiffness and damping can be taken into account in MIDAS-NX.In nonlinear dynamic time history analysis, the direct integration method is adopted.e direct integral method is a method of time as the integral parameter solution of the dynamic equilibrium equation.Analysis of dynamic time history using the Newmark method with better convergence in MIDAS-NX is performed.e basic integral method is as follows: € u t+Δt can be obtained from (19).It is substituted in formula (18), and _ u t+Δt is calculated to obtain the relationship between displacement, velocity, and acceleration of the current stage and that of the last stage as follows: Substituting formula (20) in formula ( 21), the present displacement can be obtained.e present velocity and acceleration can be obtained by present and last displacements: where � δΔt, and α and β are the integral parameters of Newmark (α � 0.5, β � 0.25).Δt is the integral time interval.Midas-NX analysis software based on finite element theory can transform differential equations into linear algebraic equations to solve problems.It is applicable to anisotropic, nonlinear, and heterogeneous materials and has an effective applicability for complex boundary conditions.It can better reveal the dynamic response law of tunnel under earthquake loading.In this paper, Midas-NX finite element software is used to analyze the nonlinear dynamic response of the tunnel.In order to reduce the boundary effect in numerical simulation, the model size is more than 5 times the diameter of the tunnel.erefore, the length, width, and height of the numerical model are 60 m, 40 m, and 55 m, respectively.
Lysmer and Kuhlemeyer [36] showed that for accurate representation of wave transmission through a model, the element size Δ l must be smaller than approximately 1/10 to 1/8 of the wavelength associated with the highest frequency component of the input wave: where λ is the wavelength of the propagated wave in the model.According to the calculation results, the total number of nodes and units in the numerical calculation model are 9,548 Advances in Civil Engineering and 42,539, respectively.In the numerical simulation, the surrounding rock and the lining are simulated by a solid element, and the elastic-plastic constitutive model and Mohr-Coulomb yield criterion are used.In the computational model, the free eld boundary and Rayleigh damping are used, and the critical damping ratio of 5% is considered.
e speci c calculation model is shown in Figure 1. e Darui arti cial wave is used as the loading wave.e time history curve of acceleration of the Darui arti cial wave and its Fourier spectrum are shown in Figure 2. Five measuring points are respectively located at the arch foot, arch shoulder, and vault of each lining inside the double hole.LF and RF refer to the left arch foot and right arch foot, respectively.LS and RS refer to the left arch shoulder and right arch shoulder, respectively.VA refers to vault.e number of measuring points is 1∼5 and 6∼10 in the left hole and right hole, respectively.e detailed measuring points are shown in Figure 3.
e physical and mechanical parameters of the surrounding rock and lining model are given in Table 1.Test.In order to verify the reliability of numerical simulation, a series of large-scale shaking table physical model tests have been carried out at the National Engineering Laboratory of the high-speed railway construction technology in Central South University (Figure 4).In the model test, the static and dynamic parameters of the model and prototype must satisfy the similarity relation.e factors such as the size of the shaking table, the relevant parameters of the test instrument, the carrying capacity, and the model boundary e ect are comprehensively considered.Finally, the geometric similarity ratio of the model is 1 : 10; the acceleration similarity ratio is 1 : 1; and the density similarity ratio is 1 : 1.Other main similarity constants can be deduced from similarity theory and dimensional analysis.

Veri cation of Shaking Table
e rigid model box was used in the shaking table test, which was made of a steel plate, section steel, and plexiglass material (Figure 4).In order to eliminate the boundary e ect of the model box, the middle sand and crushed stone were used as the friction boundary at the bottom of the model box.e polystyrene foam was pasted on the inner wall of the model box.A layer of PVC lm was pasted on the polystyrene foam board at the left and right sides of the model box.e boundary of the model was dealt within literature [37,38].
e lining model was made of microconcrete which was similar to the physical properties of concrete.Steel bar is simulated with a galvanized iron wire.According to the similarity constant, the lining thickness is determined to be 4 cm.After several times of proportioning test, the optimum lining model material ratio is 1 : 6.9 : 1.3 (cement : sand : water).e compressive strength of lining is converted to 5 MPa in accordance with the stress similitude ratio.e surrounding rock of the tunnel is divided into three layers.
e rst layer is weakly weathered rock; second layer is weak rock; and third layer is hard rock.e mixed proportion of similar materials from top to bottom is 1 : 6.30 : 1.17 e location and number of measuring points in shaking table test is the same as that in numerical simulation.One measuring point was used as the reference point, which was located at the shaking table, and the number is 11.
e reliability of numerical simulation is veri ed by the acceleration response of the tunnel.e 0.4 g Darui arti cial wave (DR-XZ) curves of the shaking table test and numerical simulation are shown in Figure 5.In the 0.4 g condition, the acceleration response peak of each measuring point in horizontal and vertical directions is shown in Table 2.
As shown in Figure 4, the time history curves of the shaking table test and numerical simulation are similar, indicating that the numerical simulation is in good agreement with the shaking table test.As shown in Table 2, the acceleration response peak value obtained by vertical simulation calculation is basically consistent with the peak value obtained by the shaking table test.e error is within the reasonable range.e results show that the shaking table test results are reasonable and the numerical simulation results are reliable.

Analysis of Tunnel Shear Force with Di erent Slope
Angles.In this paper, shear force is taken as an index to investigate the di erence of dynamic response of the tunnel under di erent existing slope-angle conditions (30 °, 45 °, 60 °, and 90 °).Numerical calculation results show the variation trends of shear force, which are shown in Figure 6.
As shown in Figure 6(a), the shear force at the arch foot and vault of the tunnel is smaller, while the shear force at the arch shoulder is larger.e overall variation trend of each hole presents an "M" shape.At the same measuring point, 4 Advances in Civil Engineering the shear force varies greatly with di erent loading peaks.e shear force increases with the increase of the seismic wave peak.Under the seismic wave peak of 0.1 g, the maximum of shear force at arch shoulder is 30 kN.Under the seismic wave peak of 0.6 g, the maximum of shear force at spandrel is 210 kN, which increases by 6 times.erefore seismic wave peak has a signi cant in uence on the shear force of the tunnel crosssection.e shear forces at both left arch shoulder and right arch shoulder are the same in the left hole, and the result is 180 kN.e arch shoulders of the right hole have di erent shear force.e result of the left arch shoulder is 140 kN, but the result of the right arch shoulder is 210 kN.And the di erence is 70 kN.e reason is that the upper rock layer of the left hole is thick and has good stability, and it is far away from the existing slope and is less a ected by the existing slope.e right hole is close to the existing slope, and the slope has a great in uence on the seismic force response of the right tunnel under the earthquake dynamic action.e right arch shoulder of the right hole is close to the existing slope, and the upper rock layer has poor stability.
By comparing Figures 6(a)-6(d), it can be seen that the variation trend of shear force of the tunnel crosssection is similar under di erent existing slope-angle conditions, but the shear strength is di erent.It can be seen that the existing slope angle never changes the seismic force response characteristics of the tunnel but only a ects the response intensity of the seismic force.With the increase of the existing slope angle, the di erence of shear force between the two arches gradually increases in the left hole.e di erence of shear force between the two arches gradually decreases in the right hole.When the existing slope angle is 90 °, the shear strength of the left hole is the same as that of the right hole.e variation trend presents symmetrical distribution.It can be seen that the existing slope has signi cant in uence on the seismic force response of the tunnel.e degree of in uence is di erent under di erent existing slope angles.In actual engineering practice, the in uence of the existing slope angle on the dynamic stability of the tunnel should be fully considered.

Analysis of Tunnel Shear Force with Di erent Clear
Distances.In this paper, shear force is taken as an index to   Advances in Civil Engineering investigate the di erence of seismic force response of the tunnel under di erent clear distance conditions (1.0 D, 1.5 D, 2.0 D, and 2.5 D, where D is the tunnel diameter).e shear force is calculated, and variation trends of shear force are shown in Figure 7.By Figure 7(a), it can be seen that the shear force of the left hole rstly increases, then decreases, and nally increases.e variation trend of shear force of the right hole is complex.
ere are di erences in the variation trend between left and right holes.e main reason is the in uence of the existing slope on the right tunnel.e shear force at the arch shoulder and arch foot is larger.
e maximum is 330 kN in 0.6 g condition.At the vault, the shear force is small.e seismic force response of the left hole is larger than that of the right hole.e reason is that the overburden thickness of the left hole is larger than that of the right hole.e inertia force of the left hole is larger under the action of seismic wave.Under the condition of 1.5 D and 2.0 D, the maximum at the right arch foot of the left hole is 920 kN.Under the 3.0 D condition, the maximum at the right arch foot of the left hole is 750 kN.
By comparing the four conditions' (1.0 D, 1.5 D, 2.0 D and 3.0 D, where D is the tunnel diameter) variation trend of shear force response and response values, it is found that under 1.5 D and 2.0 D conditions the seismic force response of the tunnel crosssection is strong, while the 1.0 D and 3.0 D response is relatively weak.
is conclusion can provide reference for practical engineering examples of small clear distance tunnels and avoid constructing tunnels with 1.5 D to 2.0 D spacing as far as possible.e magnitude of the shear force response of the tunnel crosssection is a ected by the seismic wave peak.e seismic force response increases with the increase of the seismic wave peak.A linear variation is presented.

Analysis of Tunnel Shear Force with Di erent Loading
Modes.Under di erent loading modes, the shear force is taken as an index to investigate the di erence of seismic force response of the tunnel (X direction, Z direction, XZ direction, and XYZ direction).X is the horizontal direction, Y is the tunnel axis direction, and Z is the vertical direction.
e shear force is calculated, and variation trends of shear force are shown in Figure 8.
As shown in Figure 8(a), the shear force of the tunnel crosssection shows an "M" shape trend under the action of one-way loading wave (X direction).Under the same loading peak, the shear force reaches the maximum at the arch shoulder, and the maximum is 218 kN. e minimum is in the arch foot and vault.Because of the in uence of the existing slope on the shear force response of the tunnel, the variation trends of the shear force in the left and right tunnels never show a symmetrical trend.
As shown in Figure 8(b), the shear force of the tunnel crosssection shows a "V" shape trend under the action of one-way loading wave (Z direction).Under the same loading   Advances in Civil Engineering peak, the shear force reaches the maximum at the arch foot, and the maximum is 280 kN. e minimum is in the vault.Because of the in uence of the existing slope on the shear force response of the tunnel, the variation trends of the shear force of the left and the right holes never show a symmetrical trend, and the response of the left hole is more intense.e main reason is that in the vertical direction of loading wave (Z direction), the overlying strata of the left hole are thicker; the action force of the left hole is greater; and the "rise e ect" has a signi cant impact on the shear force response of the left hole.By comparing Figures 6(a) and 6(b), it can be seen that the loading direction of the single direction loading method is di erent, and the magnitude and variation trend of tunnel shear force are di erent.It shows that the loading direction of seismic wave has a signi cant in uence on the seismic force response of the tunnel.
From Figure 8(c), it can be seen that under the action of bidirectional wave loading (XZ), the variation trend of shear force is more complex.e left hole shows the trend of rst increasing, then decreasing, and nal increasing, and the maximum is located at the arch shoulder and arch foot.e shear force of the right tunnel is larger at the arch shoulder Advances in Civil Engineering and arch foot.e left half arch shear force is larger than the right half arch in the right hole.e reason may be that the response of the vertical direction is stronger than that of the horizontal direction under the action of the bidirectional wave.In the vertical direction, the response of the right arch is weaker than that of the left half arch, and it is a ected by the dynamic response of the middle rock column due to the close distance between the left arch and the middle rock column.
As shown in Figure 8(d), the variation trend of shear force is complex under the action of bidirectional wave loading (XYZ).e shear force at the arch foot and the arch shoulder is larger, and the vault is smaller.e maximum is located at the right arch foot of the left hole.Under 0.6 g loading peak, the maximum is 700 kN. e variation tends at both the left hole and the right hole are similar, but the numerical value of shear force is di erent.e reason may be that the seismic wave acts as the main factor under the action of the three directions of seismic waves, while the in uence of the existing slope weakens.
By comparing Figures 7(a)-7(d), it can be seen that the seismic wave peak has an e ect on the shear force of the  e transverse shear force of the small clear distance tunnel shows nonlinear variation trend.Di erent loading modes have a signi cant in uence on the shear force response.Under the same loading mode, the di erent loading directions also have a signi cant in uence on the shear force response.e maximum shear force occurs at the arch shoulder and arch foot of the tunnel.erefore, the seismic forti cation measures should be strengthened at the arch shoulder and arch foot of the tunnel.Advances in Civil Engineering MIDAS-NX finite element software.e tunnel shear responses under different seismic wave peaks (0.1 g, 0.2 g, 0.3 g, 0.4 g, 0.5 g, and 0.6 g), different loading modes (X direction, Z direction, XZ direction and XYZ direction), different existing slope angles (30,45,60  shows nonlinear variation trend.e shear force response is strongest at the arch shoulder and arch foot of the tunnel.erefore, it is necessary to strengthen the seismic fortification measures at the key point.

Conclusions
(cement : sand : water) of M7.5 mortar, 1 : 7.25 : 1.35 (cement : sand : water) of M5 mortar, and 1 : 5.58 : 1.04 (cement : sand : water) of M10 mortar, respectively.e buried depth of the tunnel is 0.9 m; the clear width of the tunnel is 0.7 m; and thickness of the middle partition wall is 0.4 m.e acceleration sensor was used to record the response of the measuring point in the shaking table test.e type, range, and sensitivity were 1221L-002, 1221L-002, ±20 m•s −2 , and 2000 mv/g, respectively.

Figure 2 :
Figure 2: (a) e time history curve of acceleration of the Darui arti cial wave and (b) its Fourier spectrum.

Figure 5 :
Figure 5: Time history curves of numerical simulation and shaking table test.(a) Horizontal direction.(b) Vertical direction.
and 90 °), and different clear distances (1.0 D, 1.5 D, 2.0 D and 2.5 D) are studied.e following conclusions are obtained: (1) e shear force gradually increases with the increasing of seismic peak.e amplification is different with different measuring points.Moreover, the variation trend is invariable, although the conditions (such as different loading modes, different existing slope angles, and different clear distances) are different.(2) Under different existing slope-angle conditions, the variation trend of shear force of tunnel crosssection is similar, but the shear strength is different.e existing slope angle never changes the seismic force response characteristics of the tunnel, but only affects the seismic force response.With the increase of the existing slope angle, the difference of shear strength of the arch is becoming larger in the left hole.However, the difference in the right hole decreases gradually.e existing slope has a significant influence on the seismic force response of the tunnel, and the influence degree is different under different existing slope angles.In actual engineering practice, the influence of the existing slope angle on the dynamic stability of the tunnel should be fully considered.(3) Under different clear distance conditions, the shear force at the arch shoulder and arch foot are larger, while the shear force at the vault is smaller.Under 1.5 D and 2.0 D conditions, the shear force response of the tunnel crosssection is stronger.e shear force response of 1.0 D and 3.0 D is relatively weak.In practical engineering examples, the tunnels with 1.5 D to 2.0 D clear distance should be avoided as far as possible.(4) Under different loading modes, the shear force response of tunnel is different.Under the same loading mode, the shear force response in different loading directions is different.e loading mode of the seismic wave has a significant influence on the shear force of the tunnel crosssection.(5) e shear force response of the tunnel crosssection

Table 1 :
Physical and mechanical parameters of the surrounding rock and lining model.
e reliability of the numerical simulation is veri ed by the shaking table model test.e parametric study of a shallowbias tunnel with a clear distance is carried out by using