Dynamic Response and Failure Mechanism of Deep-Buried Tunnel with Small Net Distance under Blasting Load

: Under blasting load, a series of safety problems, such as lining cracking and surrounding rock instability, are prone to occur in deep-buried tunnels with a small net distance. It is signiﬁcant to understand the dynamic response and failure mechanism of tunnels under blasting. The blasting attenuation formula is optimized through theoretical analysis and ﬁeld experiments. The measuring point vibration is monitored in real time and the tunnel blasting model is established by ANSYS/LS-DYNA software. The model was set as having no reﬂective boundary and an uncoupled charge structure was used. The attenuation law of blasting seismic waves is studied from the adjacent tunnel lining and the direction of the tunnel cross-section and length. The inner and outer sides of the tunnel lining are investigated, respectively. The displacement and acceleration of lining measuring point are also analyzed. The dynamic response of the tunnel lining under blasting excavation is analyzed from multiple angles. The results show that the arch foot on the inner side of the lining (the side in contact with the tunnel headroom) is the ﬁrst to generate vibration. On the outside of the lining (the side in contact with the rock),the peak vibration velocity is reached after blasting load unloading. There is little difference in the vibration velocity at different positions of the transverse section, but great difference in the vibration velocity of the longitudinal section. The inﬂuence of the horizontal displacement was greater than that of the vertical displacement. The vibration acceleration of the measuring point at the arch foot of the section is the largest and the detonation is also the largest.


Introduction
The depth of a tunnel has a great impact on its stability and safety.Excessive depth will lead to a high-stress field, which may lead to rockburst and large deformation [1].The rock clamped in the small-distance tunnel is subjected to severe disturbance by excavation on both sides, and is in an unfavorable state of unidirectional or bidirectional force [2].Especially under the action of blasting load, the instability probability of a deep-buried small-distance tunnel is obviously increased, so it is necessary to explore the dynamic response of deep-buried tunnels with a small net distance under blasting load.
The secondary lining concrete in the tunnel is close to the work face and is susceptible to blasting vibration during construction [3].A large number of scholars have carried out research on the stress characteristics of the surrounding rock of tunnels with a small net distance.Li et al. [4] monitored the surrounding rock pressure behind the initial support of the section and found that the tunnel excavation changed the maximum principal stress of the surrounding rock at the arches and inverts into tensile stress.Pan et al. [5] explored the influence of tunnel blasting on the secondary lining of an intermediate railway tunnel.Through numerical calculation and comparative analysis, it was revealed that the area where the plastic zone is concentrated and the left and right hances produce large deformation.Wu et al. [6] used shaking table model tests and damage theory to study the change of shear strength and the mechanical cumulative damage model of a fault fracture zone under the blasting vibration cyclic load.In order to explore the classification standard of large compressive deformation of deep-buried tunnels, Zhu et al. [7] found by numerical simulation that compressive stress was the main force on the initial support structure under blasting load, and tensile stress was distributed at the arch bottom and arch foot.The formulas for calculating the surrounding rock stress are also different, and different tunnels focus on different factors.Based on the upper limit method of limit analysis, An et al. [8] derived the formula of surrounding rock pressure in deep tunnels with a small net distance, and analyzed the characteristics of surrounding rock pressure by using the control variable theory.Wang et al. [9] use a power function to express the relationship between the main vibration frequency and the distance from the blasting source.Wu et al. [10], aiming at the surrounding rock pressure problem of an asymmetric three-hole tunnel, comprehensively considered such factors as rock thickness, compressive strength and internal friction angle on the basis of the general theory.The calculation method of the surrounding rock pressure in an asymmetric three span tunnel with a small net distance was obtained.
Numerical simulation is an effective way to analyze tunnel dynamic response.Zhu [11] established a three-dimensional finite element model of blasting vibration to explore the dynamic response law of rock piles tunnel under blasting load.The comparison between the simulated velocity and the measured velocity at the tunnel inlet shows that the finite element software can accurately predict the dynamic response at different positions of the tunnel.Hao [12] aimed to explore the dynamic response and velocity variation rule of tunnels under blasting load under different net distances and advance distances.Based on the finite element software LS-DYNA, the model was established and the engineering example was combined to analyze the dynamic response of the tunnel to the stress and vibration velocity.Ma et al. [13] established a 4D-LSM numerical calculation model in order to study the failure mechanism of tunnels under the action of close-range explosions.The dynamic response and failure characteristics of the tunnel under different blast pressures and distances were analyzed.The vibration velocity is an important index to measure the dynamic response of a tunnel, and the propagation law of vibration velocity also reflects the failure mechanism of a tunnel under blasting.In order to explore the blasting vibration response characteristics of the advance tunnel with a small net distance and large cross-section, Ma et al. [14] monitored and analyzed the distribution law of blasting vibration velocity in the longitudinal and cross-sectional sections of the advance tunnel.It was found that the radial vibration velocity of the side wall of the tunnel was the largest, and the vibration velocity of the side facing the explosion was obviously higher than that of the side facing the explosion.There are also laws that reflect the characteristics of the dynamic response: under the action of multi-frequency blasting, the acoustic velocity of the surrounding rock fluctuates from the inside of the middle rock to the tunnel contour [15], the blasting vibration speed in the middle of the floor attenuates the fastest and the damage range of surrounding rock is the largest [16], and the blasting vibration has great influence on the horizontal direction of the tunnel [17].
To sum up, the dynamic response of a tunnel to blasting load is special and different.The current domestic and foreign studies on deep-buried small-distance tunnels mainly focus on the stress characteristics of the surrounding rock and the propagation law of different tunnel sections.Therefore, against the background of the blasting excavation of three-way separation tunnel, the theory of data analysis is firstly expounded.Then the field blasting experiment is used to explore the law of seismic wave propagation in the tunnel.At the same time, the dynamic response of the tunnel is further explored with the help of numerical simulation software.This paper mainly explores the dynamic response inside and outside the lining and the transverse and longitudinal vibration velocity propagation law of the tunnel.It also analyzes the horizontal and vertical displacement and acceleration changes of the tunnel lining.Finally, the dynamic response and damage mechanism of a deep-buried tunnel with a small net distance under blasting load are obtained.

Methodology
In the experimental data processing, the linear regression and nonlinear regression are used to compare and analyze the Sadowski formula.
At present, the Sadowski empirical formula is widely recognized at home and abroad for the prediction of blasting vibration speed, and the Sadowski formula has been incorporated into China's Blasting Safety Regulations (GB6722-2014).The Sadowski empirical formula: V is peak velocity of blasting vibration; Q is blasting charge; R is blasting distance; K is the coefficient related to medium and blasting conditions; ∝ is blasting vibration attenuation index.
(1) Linear regression analysis Take the log of both sides of this equation: (2) Nonlinear regression analysis [18] Let Q 1/3 /R = q, then V = Kq α The sum of squares of nonlinear residuals is M 2 gets the minimum value, and the necessary condition of the function above is construct nonlinear equation about α after solving for α, K can be obtained from Equation ( 10)

Project Summary
Badaling Great Wall Station is located in the new Badaling Tunnel.Its location is shown in the Figure 1.The total length of the station is 470 m.The underground construction area of the station is 36,000 square meters and the depth between rail surface and ground is 102 m.With many stations and a large number of variable types of caverns, it is a complex underground cave group of stations in China.The underground structure of the station has three floors, as shown in Figure 2, which are platform floor, inbound channel floor, outbound channel floor and equipment cavern.The levels of surrounding rock are mainly III and V.The horizontal spacing of the standard tunnel section is 2.27~6 m.The vertical spacing between the platform layer and the passageway layer is 4.55 m.The clear passage distance of the escalator in and out of the station building is 4.1~3.78m.

∑
V q K ∑ q = ∑ V q ln q K ∑ q ln q = 1 (10) construct nonlinear equation about α after solving for α, K can be obtained from Equation ( 10)

Project Summary
Badaling Great Wall Station is located in the new Badaling Tunnel.Its location is shown in the Figure 1.The total length of the station is 470 m.The underground construction area of the station is 36,000 square meters and the depth between rail surface and ground is 102 m.With many stations and a large number of variable types of caverns, it is a complex underground cave group of stations in China.The underground structure of the station has three floors, as shown in Figure 2   The overall construction plan is to excavate the left line first, followed by the right line and the middle line.This is to select a reasonable tunneling scheme according to the parameters involved [19].The surrounding rock of grade Ⅲ is constructed by bench blasting.The upper and lower steps are staggered by 5~8 m, and the second lining is 51~80 m away from the tunnel face.The surrounding rock of Grade Ⅴ adopts the temporary transverse support method of three steps, the construction interval of the upper, middle and lower steps is 3~5 m, and the second lining is not more than 70 m away from the tunnel face.Emulsified explosives and plastic detonators are selected as blasting equipment.Millisecond differential initiation is used for blasting.The section size and spacing of the three holes are shown in Figure 3.

Design of Blasting Parameters for the Standard Section of Three Holes Separation
Bench blasting is adopted for the three-hole separation of the grade Ⅲ surrounding rock section, the step length is 5~8 m, and the blasting construction is carried out in the order of up and down [20].The borehole is drilled by hand-held pneumatic drill, the diameter of the borehole is 42 mm, and the blasting explosive is Φ32 mm rock emulsion explosive.In terms of charging structure, layered charging [21] and uncoupled charging are adopted.Taking the third-level surrounding rock of the middle hole in the three-hole separation section as an example, the middle hole in the cut blasting plays the role of the free surface [22], and its hole layout is shown in Figures 4 and 5.The overall construction plan is to excavate the left line first, followed by the right line and the middle line.This is to select a reasonable tunneling scheme according to the parameters involved [19].The surrounding rock of grade III is constructed by bench blasting.The upper and lower steps are staggered by 5~8 m, and the second lining is 51~80 m away from the tunnel face.The surrounding rock of Grade V adopts the temporary transverse support method of three steps, the construction interval of the upper, middle and lower steps is 3~5 m, and the second lining is not more than 70 m away from the tunnel face.Emulsified explosives and plastic detonators are selected as blasting equipment.Millisecond differential initiation is used for blasting.The section size and spacing of the three holes are shown in Figure 3.The overall construction plan is to excavate the left line first, followed by the right line and the middle line.This is to select a reasonable tunneling scheme according to the parameters involved [19].The surrounding rock of grade Ⅲ is constructed by bench blasting.The upper and lower steps are staggered by 5~8 m, and the second lining is 51~80 m away from the tunnel face.The surrounding rock of Grade Ⅴ adopts the temporary transverse support method of three steps, the construction interval of the upper, middle and lower steps is 3~5 m, and the second lining is not more than 70 m away from the tunnel face.Emulsified explosives and plastic detonators are selected as blasting equipment.Millisecond differential initiation is used for blasting.The section size and spacing of the three holes are shown in Figure 3.

Design of Blasting Parameters for the Standard Section of Three Holes Separation
Bench blasting is adopted for the three-hole separation of the grade Ⅲ surrounding rock section, the step length is 5~8 m, and the blasting construction is carried out in the order of up and down [20].The borehole is drilled by hand-held pneumatic drill, the diameter of the borehole is 42 mm, and the blasting explosive is Φ32 mm rock emulsion explosive.In terms of charging structure, layered charging [21] and uncoupled charging are adopted.Taking the third-level surrounding rock of the middle hole in the three-hole separation section as an example, the middle hole in the cut blasting plays the role of the free surface [22], and its hole layout is shown in Figures 4 and 5.  Bench blasting is adopted for the three-hole separation of the grade III surrounding rock section, the step length is 5~8 m, and the blasting construction is carried out in the order of up and down [20].The borehole is drilled by hand-held pneumatic drill, the diameter of the borehole is 42 mm, and the blasting explosive is Φ32 mm rock emulsion explosive.In terms of charging structure, layered charging [21] and uncoupled charging are adopted.Taking the third-level surrounding rock of the middle hole in the three-hole separation section as an example, the middle hole in the cut blasting plays the role of the free surface [22], and its hole layout is shown in Figures 4 and 5.

Vibration Monitoring Scheme
A total of seven measuring points are arranged for the vibration test, and their positions are shown in Figure 6.The mileage corresponding to the face of the middle hole is DK68+158, and the thickness of the wall is 6 m.Measuring point 2 is on the same plane as the charging section on the tunnel face.Measuring point 1 is set in the excavated stage of the middle tunnel, 5 m away from measuring point 2. Measuring points 3, 4, and 5 are set in the unexploded area of the middle tunnel.The distances from measuring point 2 are 5 m, 13 m and 23 m.Measuring points 6 and 7 are on the same plane with measuring point 2, but 0.8 m lower than that of measuring point 2. By monitoring the vibration velocity of the measuring point in the same plane and a measuring line, the propagation law of blasting vibrations in different positions is analyzed.

Monitoring Results and Analysis of Blasting Vibration 3.3.1. Vibration Monitoring Results
Tunnel blasting vibration signals have chaotic characteristics.As the distance between the blasting operation surface and the detonation center of existing tunnels decreases, the chaotic characteristics of the blasting vibration response are enhanced.Therefore, drilling and blasting parameter optimization and real-time monitoring of blasting vibration should be strengthened to ensure construction quality and blasting safety [23].In order to study the vibration response of the existing tunnel, the velocity vibration waveform of measuring point 2 was selected for analysis, as shown in Figure 7.
decreases, the chaotic characteristics of the blasting vibration response are enha Therefore, drilling and blasting parameter optimization and real-time monitori blasting vibration should be strengthened to ensure construction quality and bla safety [23].In order to study the vibration response of the existing tunnel, the ve vibration waveform of measuring point 2 was selected for analysis, as shown in Fig   The field measured data of each measurement point are shown in Table 1.
Table 1.Field measured data graph.The field measured data of each measurement point are shown in Table 1.As can be seen from Table 1 and Figure 7, the maximum vibration velocity of tunnel blasting appears in cut blasting, because cut blasting has only one free surface and it is under a lot of pressure.A large part of the explosion energy is directly transformed into a vibration wave and transmitted, so the blasting vibration speed is relatively high.After the cutting surface is formed, the free surface is provided for the subsequent blasting.Therefore, although the total amount of subsequent explosions is large, the vibration generated gradually decreases, and the impact on the existing tunnel is smaller.According to the above analysis, it is necessary to start with the cut hole vibration to solve the problem of excessive vibration velocity.At the same time, the comparison of the three-way vibration waveform shows that the radial peak vibration velocity has the largest fluctuation range, and the tangential vibration velocity has the longest fluctuation time.

Vibration
Due to the small distance between the two tunnels, the maximum vibration velocity of measuring point 2, which is directly opposite to the middle of the charge position, reaches 27.95 cm/s, but the attenuation speed is also relatively fast.When the horizontal distance from the tunnel face is 13 m, the vibration velocity of measuring point 4 has dropped to 5.22 cm/s.As can be seen from Table 1, the vibration velocity of measuring point 1 is 16.45 cm/s, while that of measuring point 3 is 17.22 cm/s.The former is smaller than the latter, indicating that although the propagation distance of the vibration wave is the same, the rock mass on the excavated side has been disturbed.The surrounding rock is damaged, and the vibration velocity transmitted to the measuring point is reduced.

Analysis of Measured Data
The measured data were analyzed by linear regression according to the least-squares method.The linear correlation between the blasting vibration velocity and the distance to the center of the blast is established, and the corresponding correlation coefficient (R 2 ) is given.The results are shown in Table 2.The Sadowski formula suitable for the project is obtained from the above regression fitting results, as shown in Table 3.

Formula of Vibration Velocity Attenuation
Correlation Coefficient The results of the nonlinear regression and Sadowski attenuation formula suitable for this project are shown in Tables 4 and 5.

Formula of Vibration Velocity Attenuation Correlation Coefficient
From the above two analysis results, it can be seen that there is a certain gap between the linear regression results and the nonlinear regression.This is caused by the fact that the sum of squares of the residuals in the linear regression cannot obtain the minimum value [18], so the nonlinear results are taken as the benchmark.The cumulative damage of the tunnel under blasting load increases with the increase in blasting times.However, the damage increments decrease gradually [24].Both changes show a nonlinear relationship.
The velocity changes in each direction of each measuring point are plotted according to Table 1, as shown in Figure 8.
As can be seen from Figure 8, the variation law of vibration velocity in x direction and resultant velocity is roughly consistent.The velocity in x direction is the highest, followed by the magnitude of the velocity in the y direction and the velocity in the z direction.The distance between measuring point 1 and the explosion source is 5 m; measuring point 3 has the same distance.The vibration velocity in x direction of measuring point 1 is slightly less than that of measuring point 3 due to the damage of the excavated surrounding rock.Measuring point 2 and the explosion source are in the same cross-section, and the distance is the closest, so the vibration velocity is the maximum.Measuring points 3 to 5 are distributed longitudinally along the tunnel, and the vibration velocity decreases with the increase in the distance from the detonation source.Measuring points 6 and 7 are in the same cross-section as measuring point 2. It can be seen from the figure that the vibration velocity is inversely proportional to the detonation source distance.
value [18], so the nonlinear results are taken as the benchmark.The cumulative damage of the tunnel under blasting load increases with the increase in blasting times.However, the damage increments decrease gradually [24].Both changes show a nonlinear relationship.
The velocity changes in each direction of each measuring point are plotted according to Table 1, as shown in Figure 8.As can be seen from Figure 8, the variation law of vibration velocity in x direction and resultant velocity is roughly consistent.The velocity in x direction is the highest, followed by the magnitude of the velocity in the y direction and the velocity in the z direction.The distance between measuring point 1 and the explosion source is 5 m; measuring point 3 has the same distance.The vibration velocity in x direction of measuring point 1 is slightly less than that of measuring point 3 due to the damage of the excavated surrounding rock.Measuring point 2 and the explosion source are in the same cross-section, and the distance is the closest, so the vibration velocity is the maximum.Measuring points 3 to 5 are distributed longitudinally along the tunnel, and the vibration velocity decreases with the increase in the distance from the detonation source.Measuring points 6 and 7 are in the same cross-section as measuring point 2. It can be seen from the figure that the vibration velocity is inversely proportional to the detonation source distance.

Numerical Simulation of Dynamic Response of Lining 4.1. Parameter Selection
The third section mainly carries on the field test and data analysis according to the general situation of the project.At the same time, the field blasting vibration propagation law is predicted by the traditional blasting vibration attenuation formula.The distribution law of the vibration velocity along different directions is obtained.Also, it is necessary to use numerical simulation software to study the blasting vibration propagation law more accurately.The help of numerical simulation software can be a good complement to the field experiment.The advantages of the numerical simulation method in instantaneous force, material non-linearity, and other complex mathematical problems are incomparable to the traditional experience.
This paper is based on the ANSYS/LS-DYNA finite element program to simulate the dynamic response of existing tunnels to blasting.In the separation section of three holes in this paper, the main tunnel is 11.7 m high and 14.0 m wide, and the tunnels on both sides are 11.1 m high and 12.6 m wide.The net distance of the tunnel is 6 m.Therefore, the selected boundary size of the model is 200 m horizontally, 196 m vertically and 100 m vertically, that is, the model size is 200 m × 196 m × 100 m, as shown in Figure 9.The uncoupled charging structure is adopted.At the same time, the 3D model is set to no reflective boundary, so that the influence of reflection on the solution can be avoided when the wave reaches the interface.Finally, non-reflective boundary conditions are applied to the front, back, left, right, and lower boundaries of the model.
the selected boundary size of the model is 200 m horizontally, 196 m vertically and 100 m vertically, that is, the model size is 200 m × 196 m × 100 m, as shown in Figure 9.The uncoupled charging structure is adopted.At the same time, the 3D model is set to no reflective boundary, so that the influence of reflection on the solution can be avoided when the wave reaches the interface.Finally, non-reflective boundary conditions are applied to the front, back, left, right, and lower boundaries of the model.The algorithm adopted in this simulation is the shared node method, and the selected unit is the 3D solid164 element [25].
There is no material model of explosives included in ANSYS pre-processing.The linear elastic material model is used in the modeling of explosives.After the K file is created, the linear elastic material model parameters are changed to *MAT_HIGH_EXPLOSIVE_BURN.The LS-DYNA 3D program describes the pressurevolume relationship of the detonation products of high energy explosives using the JWL equation of state.The unit pressure of the detonation product of high energy explosives can be obtained by the equation of state.The relationship of the JWL equation of state is as follows: V is relative volume; E is initial internal energy density; A, B, R , R , ω are constants determined by experiment.The constitutive model of surrounding rock is the HJC (Holmquist-Johnson-Cook) material model [26,27].The model is mainly applied to the simulation of concrete and rock with high strain rate and large deformation.The HJC model is defined by *MAT_JOHNSON_HOLMGUIST_CONCRETE.The parameter The algorithm adopted in this simulation is the shared node method, and the selected unit is the 3D solid164 element [25].
There is no material model of explosives included in ANSYS pre-processing.The linear elastic material model is used in the modeling of explosives.After the K file is created, the linear elastic material model parameters are changed to *MAT_HIGH_EXPLOSIVE_BURN.The LS-DYNA 3D program describes the pressure-volume relationship of the detonation products of high energy explosives using the JWL equation of state.The unit pressure of the detonation product of high energy explosives can be obtained by the equation of state.The relationship of the JWL equation of state is as follows: V is relative volume; E 0 is initial internal energy density; A, B, R 1 , R 2 , ω are constants determined by experiment.The constitutive model of surrounding rock is the HJC (Holmquist-Johnson-Cook) material model [26,27].The model is mainly applied to the simulation of concrete and rock with high strain rate and large deformation.The HJC model is defined by *MAT_JOHNSON_HOLMGUIST_CONCRETE.The parameter selection of III level surrounding rock is shown in Table 6.ρ is density; e is elasticity modulus; ν is Poisson's ratio; γ is unit weight; c is cohesive force; σ is tensile limit; ϕ is internal friction angle.The secondary lining and invert filling were simulated using elastic solid elements.Shell elements were used to simulate the initial shotcrete, and parameters were selected as shown in Table 7.That is: In Equation ( 19), the column coefficient can be determined by the following equation: ω min is the minimum center frequency, which is related to both the natural frequency and the input load frequency of the system under study, and the frequency is 50 Hertz in the calculation; ζ min is the minimum critical damping ratio, which is related to the material properties of rock mass, and the ratio is 0.5% in the calculation.

Simulated Result
This paper mainly studies the effect of the blasting construction of an excavated tunnel on existing tunnels.The dynamic response rules of blasting excavation at different positions and fault faces on the same section were studied, taking a cross-section in the same plane with the tunnel face, which is called section 0, and four cross-sections behind the tunnel face.These sections have the same longitudinal monitoring points.They are section 5, section 13 and section 23, as shown in Figure 10. Figure 11 shows the measuring point of section 0, and other sections are consistent with it.The comparison between the numerical simulation and measured vibration velocity is shown in Table 8.As can be seen from Table 8, the vibration velocity monitored on site is slightly smaller than the numerical calculation data, but the two are numerically similar.Yong [28], Li [29], and F. Kirzhneret et al. [30] all analyzed the dynamic response of the tunnel with the help of finite element software.This shows that the model established by ANSYS/LS-DYNA can better reflect the actual tunnel blasting construction process.

Analysis of Tunnel Force under the Influence of Dead Weight
After the completion of excavation, the buried tunnel will be affected by the gravity of its top rock mass, and the tunnel will produce certain displacement and stress concentrations.The numerical simulation of the tunnel under the action of gravity on the upper rock mass is carried out to find out the stress concentration and the most unfavorable position of possible displacement.The cloud picture of the stress and displacement of the three-way tunnel under the influence of dead weight is shown in Figures 12-16.Through the simulation of the tunnel's dynamic response under the action of dead weight, the conclusions are as follows: (1) As can be seen Figures 12-14, in the X and Z directions, the tunnel is mainly affected by the stress released by the surrounding rock mass due to excavation, and there is no large stress concentration around the tunnel in the X direction.In the Z direction there is a certain stress concentration in the vault and the hance.In the Y direction i.e., the direction of the dead weight of rock mass, the stress concentration is greate at the vault and above, while the stress is smaller at the arch foot and tunnel bottom (2) In Figure 15, under the action of X, Y and Z directions, the stress concentration is the largest at the hance, followed by the vault and the arch foot.The surrounding rock within a certain range will produce a certain stress.(3) As can be seen from Figure 16, the rock mass around the tunnel does not produce displacement under the action of rock mass gravity.This is because the tunnel is deep underground and gravity cancels out under the action of ground stress and surrounding rock stress.
The location of stress concentration in different directions is different.According to the magnitude of influence, we should focus on the influence of self-weight in Y direction and stress concentration at the hance.Through the simulation of the tunnel's dynamic response under the action of dead weight, the conclusions are as follows: (1) As can be seen Figures 12-14, in the X and Z directions, the tunnel is mainly affected by the stress released by the surrounding rock mass due to excavation, and there i no large stress concentration around the tunnel in the X direction.In the Z direction there is a certain stress concentration in the vault and the hance.In the Y direction i.e., the direction of the dead weight of rock mass, the stress concentration is greate at the vault and above, while the stress is smaller at the arch foot and tunnel bottom (2) In Figure 15, under the action of X, Y and Z directions, the stress concentration is the largest at the hance, followed by the vault and the arch foot.The surrounding rock within a certain range will produce a certain stress.(3) As can be seen from Figure 16, the rock mass around the tunnel does not produc displacement under the action of rock mass gravity.This is because the tunnel is deep underground and gravity cancels out under the action of ground stress and surrounding rock stress.
The location of stress concentration in different directions is different.According to the magnitude of influence, we should focus on the influence of self-weight in Y direction and stress concentration at the hance.Through the simulation of the tunnel's dynamic response under the action of dead the conclusions are as follows: (1) As can be seen Figures 12-14, in the X and Z directions, the tunnel is mainly affected by the stress released by the surrounding rock mass due to excavation, and there is no large stress concentration around the tunnel in the X direction.In the Z direction, there is a certain stress concentration in the vault and the hance.In the Y direction, i.e., the direction of the dead weight of rock mass, the stress concentration is greater at the vault and above, while the stress is smaller at the arch foot and tunnel bottom.(2) In Figure 15, under the action of X, Y and Z directions, the stress concentration is the largest at the hance, followed by the vault and the arch foot.The surrounding rock within a certain range will produce a certain stress.(3) As can be seen from Figure 16, the rock mass around the tunnel does not produce displacement under the action of rock mass gravity.This is because the tunnel is deep underground and gravity cancels out under the action of ground stress and surrounding rock stress.
The location of stress concentration in different directions is different.According to the magnitude of influence, we should focus on the influence of self-weight in Y direction and stress concentration at the hance.

Vibration Velocity Analysis of the Left Tunnel Lining of a Deep-Buried Tunnel with Small Interval
The inner and outer sides of the lining were analyzed for their respective dynamic response.The inner side was the side in contact with the tunnel headroom, and the outer side was the side in contact with the surrounding rock.
(1) Analysis of vibration velocity distribution on inner lining Figure 17 shows the nephogram of the vibration velocity corresponding to 1 ms, 3 ms, 5 ms, and 7 ms on the inner side of the lining.
response.The inner side was the side in contact with the tunnel headroom, and the outer side was the side in contact with the surrounding rock.
(1) Analysis of vibration velocity distribution on inner lining As can be seen from Figure 17: (1) As shown in Figure 17a, on the inside of the lining, the blasting vibration velocity first occurs at the arch foot, and its vibration velocity is also small, about 7.52 cm/s; (2) As can be seen from Figure 17b, the vibration velocity at the arch foot increases to 15.25 cm/s in a very short time and extends to the tunnel hance.It can be found in Figure 17c that the peak vibration velocity at the hance reaches 23.25 cm/s and the vibration velocity spreads in the upper part and the bottom part.(3) As can be seen from Figure 17d, the vibration velocity continues to spread along the tunnel lining towards the vault and tunnel bottom.In the whole blasting process, the maximum vibration velocity is at the arch foot followed by the tunnel hance, and the minimum velocity is at the vault.After that the blasting load gradually decreases, and the vibration velocity decreases rapidly until it approaches 0.
(2) Analysis of vibration velocity distribution on external lining Figure 18 shows the nephogram of vibration velocity corresponding to 1 ms, 3 ms, 5 ms, and 7 ms on the external side of the lining.As can be seen from Figure 17: As shown in Figure 17a, on the inside of the lining, the blasting vibration velocity first occurs at the arch foot, and its vibration velocity is also small, about 7.52 cm/s; (2) As can be seen from Figure 17b, the vibration velocity at the arch foot increases to 15.25 cm/s in a very short time and extends to the tunnel hance.It can be found in Figure 17c that the peak vibration velocity at the hance reaches 23.25 cm/s, and the vibration velocity spreads in the upper part and the bottom part.(3) As can be seen from Figure 17d, the vibration velocity continues to spread along the tunnel lining towards the vault and tunnel bottom.
In the whole blasting process, the maximum vibration velocity is at the arch foot, followed by the tunnel hance, and the minimum velocity is at the vault.After that, the blasting load gradually decreases, and the vibration velocity decreases rapidly until it approaches 0.
(2) Analysis of vibration velocity distribution on external lining  As can be seen from Figure 18: (1) At the initial stage of blasting load, the lining of the blasting side of the existing left tunnel is affected by vibration first, and a vibration velocity circle is formed The vibration velocity of the vibration velocity circle is relatively small, and the scope of action is not large, as shown in Figure 18a,b.The vibration velocity o the arch foot in the vibration velocity circle is the largest.(2) When the blasting load continues to increase, the vibration velocity circle increases rapidly and so does the vibration velocity, but it does not reach the maximum value, as shown in Figure 18c.(3) With the increase in time, the blasting stress wave continues to propagate and reflect in the medium, and reaches a maximum value of 30.55 cm/s in Figure 18d.At this time, the blasting load has been unloaded and then vibrates.The peak value of dynamic velocity decreases gradually until it approaches 0. Compared with the inner lining, the velocity propagation of the outer lining shows that the peak vibration velocity at every moment is greater than that of the inner lining, and the peak vibration velocity at 7 ms appears at the hance.This is because the outer lining is closer to the detonation source and the detonation wave propagates more rapidly. (

3) Analysis of transverse vibration velocity distribution
The vibration velocities of the measured points at each monitoring section are shown in Table 9, and the variation of peak vibration velocity is shown in Figure 19.As can be seen from Figure 18: (1) At the initial stage of blasting load, the lining of the blasting side of the existing left tunnel is affected by vibration first, and a vibration velocity circle is formed.The vibration velocity of the vibration velocity circle is relatively small, and the scope of action is not large, as shown in Figure 18a,b.The vibration velocity of the arch foot in the vibration velocity circle is the largest.(2) When the blasting load continues to increase, the vibration velocity circle increases rapidly and so does the vibration velocity, but it does not reach the maximum value, as shown in Figure 18c.(3) With the increase in time, the blasting stress wave continues to propagate and reflect in the medium, and reaches a maximum value of 30.55 cm/s in Figure 18d.At this time, the blasting load has been unloaded and then vibrates.
The peak value of dynamic velocity decreases gradually until it approaches 0.
Compared with the inner lining, the velocity propagation of the outer lining shows that the peak vibration velocity at every moment is greater than that of the inner lining, and the peak vibration velocity at 7 ms appears at the hance.This is because the outer lining is closer to the detonation source and the detonation wave propagates more rapidly. (

3) Analysis of transverse vibration velocity distribution
The vibration velocities of the measured points at each monitoring section are shown in Table 9, and the variation of peak vibration velocity is shown in Figure 19.  Figure 20 shows the waveform of vibration velocity at the measured points of each section.It can be seen from Figures 19 and 20 that: (1) The vibration velocity is the largest at the arch foot of the blasting side lining of the existing tunnel.The velocity of hance is slightly smaller than that of arch foot but larger than that of vault.As section 0 is very close to the detonation source, Figure 20 shows the waveform of vibration velocity at the measured points of each section.Figure 20 shows the waveform of vibration velocity at the measured points of each section.It can be seen from Figures 19 and 20 that: (1) The vibration velocity is the largest at the arch foot of the blasting side lining of the existing tunnel.The velocity of hance is slightly smaller than that of arch foot but larger than that of vault.As section 0 is very close to the detonation source, It can be seen from Figures 19 and 20 that: (1) The vibration velocity is the largest at the arch foot of the blasting side lining of the existing tunnel.The velocity of hance is slightly smaller than that of arch foot but larger than that of vault.As section 0 is very close to the detonation source, the vibration speed reaches 30.55 cm/s, but it rapidly decreases to 6.55 cm/s at section 23, as shown in Figure 19.(2) Section 0 first reaches the peak value of blasting vibration velocity, and its blasting vibration velocity rapidly decreases after reaching the maximum value.
The blasting vibration wave attenuates the fastest.Section 5 is the second, and its blasting vibration attenuation rate is less than that of section 0. The attenuation velocity of section 13 is obviously less than that of sections 0 and 5 after it reaches the peak velocity of blasting vibration.The blasting vibration velocity of section 23 reaches the peak time the latest, and the attenuation speed is also the smallest.
It can be seen that with the increase in detonation source distance, the time for the measuring point to reach the blasting vibration peak velocity also increases.After reaching the peak vibration velocity, the attenuation velocity also decreases with the increase in detonation source distance.The vibration velocity of the cross-section has little change.However, with the increase in longitudinal distance, there is a big difference in the vibration velocity of different sections.the vibration speed reaches 30.55 cm/s, but it rapidly decreases to 6.55 cm/s at section 23, as shown in Figure 19.
(2) Section 0 first reaches the peak value of blasting vibration velocity, and its blasting vibration velocity rapidly decreases after reaching the maximum value.
The blasting vibration wave attenuates the fastest.Section 5 is the second, and its blasting vibration attenuation rate is less than that of section 0. The attenuation velocity of section 13 is obviously less than that of sections 0 and 5 after it reaches the peak velocity of blasting vibration.The blasting vibration velocity of section 23 reaches the peak time the latest, and the attenuation speed is the smallest.It can be seen that with the increase in detonation source distance, the time for the measuring point to reach the blasting vibration peak velocity also increases.After reaching the peak vibration velocity, the attenuation velocity also decreases with the increase in detonation source distance.The vibration velocity of the cross-section has little change.However, with the increase in longitudinal distance, there is a big difference in the vibration velocity of different sections.the vibration speed reaches 30.55 cm/s, but it rapidly decreases to 6.55 cm/s at section 23, as shown in Figure 19.
(2) Section 0 first reaches the peak value of blasting vibration velocity, and its blasting vibration velocity rapidly decreases after reaching the maximum value.
The blasting vibration wave attenuates the fastest.Section 5 is the second, and its blasting vibration attenuation rate is less than that of section 0. The attenuation velocity of section 13 is obviously less than that of sections 0 and 5 after it reaches the peak velocity of blasting vibration.The blasting vibration velocity of section 23 reaches the peak time the latest, and the attenuation speed is also the smallest.It can be seen that with the increase in detonation source distance, the time for the measuring point to reach the blasting vibration peak velocity also increases.After reaching the peak vibration velocity, the attenuation velocity also decreases with the increase in detonation source distance.The vibration velocity of the cross-section has little change.However, with the increase in longitudinal distance, there is a big difference in the vibration velocity of different sections.As can be seen from Figures 21 and 22: (1) The blasting vibration velocity of section 0 and section 5 (i.e., n source) attenuates the fastest after reaching the peak.However, th rate of the vibration velocity of section 13 and section 23 is signif than that of the former, and their attenuation velocity roughly th (2) In the longitudinal direction of the existing tunnel (i.e., in the measuring points 2, 3, 4 and 5), the vibration velocity of the blasting existing tunnel decreases gradually.This indicates that the vibratio the measuring point decreases rapidly with the increase of the di blasting source.(3) With the increase in the detonation source distance, the moment w sections reach the peak vibration velocity is different.The mome measuring points of section 0 and section 5 reach the peak vibrati almost the same, while 13 and 23 reach the peak vibration veloci indicates that the surrounding rock vibration frequencies are simil of section 0 to section 5.The vibration velocity at the arch foot of se 23 is 0 in 0-2 ms.The sections close to the blasting source (section 0, section 5) have simil vibration velocity and is greatly affected.The two sections far aw detonation source (section 13 and section 23) have similar changes velocity.The main vibration band is within 0 to 2 ms.

Analysis of Displacement of Left Tunnel
A basic model diagram of tunnel nodes is shown in Figure 23, with nodes and 24-45 external nodes.As can be seen from Figures 21 and 22: (1) The blasting vibration velocity of section 0 and section 5 (i.e., near the blast source) attenuates the fastest after reaching the peak.However, the attenuation rate of the vibration velocity of section 13 and section 23 is significantly lower than that of the former, and their attenuation velocity is roughly the same.
In the longitudinal direction of the existing tunnel (i.e., in the directions of measuring points 2, 3, 4 and 5), the vibration velocity of the blasting lining of the existing tunnel decreases gradually.This indicates that the vibration velocity of the measuring point decreases rapidly with the increase of the distance to the blasting source.
With the increase in the detonation source distance, the moment when different sections reach the peak vibration velocity is different.The moment when the measuring points of section 0 and section 5 reach the peak vibration velocity is almost the same, while 13 and 23 reach the peak vibration velocity later.This indicates that the surrounding rock vibration frequencies are similar in the area of section 0 to section 5.The vibration velocity at the arch foot of sections 13 and 23 is 0 in 0-2 ms.
The sections close to the blasting source (section 0, section 5) have similar changes in vibration velocity and is greatly affected.The two sections far away from the detonation source (section 13 and section 23) have similar changes in vibration velocity.The main vibration band is within 0 to 2 ms.

Analysis of Displacement of Left Tunnel
A basic model diagram of tunnel nodes is shown in Figure 23, with 1-24 internal nodes and 24-45 external nodes.
An overall lining deformation diagram, the X and Y direction and vector synthesis displacement nephograms, are shown in Figures 24-27.
According to Figures 24-26, it can be seen that the maximum horizontal displacement of the tunnel occurs on the inner node, namely, the inner side of the lining, while the maximum vertical displacement of the tunnel occurs on the outer node, namely, the outer side of the lining.Compared with the horizontal displacement, the vertical displacement fluctuates greatly, but it still has relevant rules.Different from the horizontal displacement, the vertical displacement fluctuates greatly on the power exponential curve.This shows that its changes are more complex and more random.Of course, the overall vertical displacement of the tunnel is smaller than that of the horizontal displacement.It can be seen that under the condition of the same distance in most cases, the horizontal displacement should be paid more attention than the vertical displacement [31].According to Figures 24-26, it can be seen that the maximum horizontal displacement of the tunnel occurs on the inner node, namely, the inner side of the lining, while the maximum vertical displacement of the tunnel occurs on the outer node, namely, the outer side of the lining.Compared with the horizontal displacement, the vertical displacement fluctuates greatly, but it still has relevant rules.Different from the horizontal displacement, the vertical displacement fluctuates greatly on the power exponential curve.This shows that its changes are more complex and more random.Of course, the overall vertical displacement of the tunnel is smaller than that of the horizontal displacement.It can be seen that under the condition of the same distance in most cases, the horizontal displacement should be paid more attention than the vertical displacement [31].
It can be seen from Figure 27 that the influence range of displacement vector synthesis is mainly the inner side of the lining, which also confirms that more attention should be paid to the horizontal displacement mentioned above.The maximum displacement occurs at the tunnel vault, followed by the displacement at the hance.The variation of displacement with numbering is shown in Figure 28.It can be seen from Figure 27 that the influence range of displacement vector synthesis is mainly the inner side of the lining, which also confirms that more attention should be paid to the horizontal displacement mentioned above.The maximum displacement occurs at the tunnel vault, followed by the displacement at the hance.The variation of displacement with numbering is shown in Figure 28.According to Figures 24-26, it can be seen that the maximum horizontal displacement of the tunnel occurs on the inner node, namely, the inner side of the lining, while the maximum vertical displacement of the tunnel occurs on the outer node, namely, the outer side of the lining.Compared with the horizontal displacement, the vertical displacement fluctuates greatly, but it still has relevant rules.Different from the horizontal displacement, the vertical displacement fluctuates greatly on the power exponential curve.This shows that its changes are more complex and more random.Of course, the overall vertical displacement of the tunnel is smaller than that of the horizontal displacement.It can be seen that under the condition of the same distance in most cases, the horizontal displacement should be paid more attention than the vertical displacement [31].
It can be seen from Figure 27 that the influence range of displacement vector synthesis is mainly the inner side of the lining, which also confirms that more attention should be paid to the horizontal displacement mentioned above.The maximum displacement occurs at the tunnel vault, followed by the displacement at the hance.The variation of displacement with numbering is shown in Figure 28.

Acceleration Analysis of Measured Points
Figure 29 shows the time-history curve of blasting vibration acceleration at each measuring point at section 0.

Acceleration Analysis of Measured Points
Figure 29 shows the time-history curve of blasting vibration acceleration at each measuring point at section 0.  It can be seen from Figure 29 that the vibration acceleration of the measuring point at the arch foot of section 0 on the blasting side of the left tunnel is the largest, followed by the measuring point at the tunnel hance, and the peak value of vibration acceleration of the vault is the smallest.By analyzing the change of the blasting vibration acceleration, the resultant force change law of the blasting side can be directly observed: the blasting vibration acceleration at the arch foot is the largest, indicating that the outside lining of the blasting side has the largest detonation force.

Conclusions
Through the combination of field test and numerical simulation, this paper analyzes the dynamic response of a deep-buried tunnel with small interval under blasting load, and the specific conclusions are as follows: 1. Through the spectrum analysis of the field measured waveform and the nonlinear regression of the data, the empirical formula conforming to the attenuation law of the Sadowski vibration velocity of the project is obtained.2. On the inside of the lining (the side contacting the tunnel headroom), the vibration continues from the arch foot along the tunnel lining towards the vault and tunnel bottom.On the outside of the lining (the side contacting the surrounding rock), the lining on the blasting side is affected by vibration before the back side at the initial stage of blasting.The vibration velocity in the order of large to small is arch foot, tunnel hance, and vault.3.After the blasting action is completed, the vibration velocity of the lining structure gradually decreases but will not disappear, and will only approach 0. The vibration It can be seen from Figure 29 that the vibration acceleration of the measuring point at the arch foot of section 0 on the blasting side of the left tunnel is the largest, followed by the measuring point at the tunnel hance, and the peak value of vibration acceleration of the vault is the smallest.By analyzing the change of the blasting vibration acceleration, the resultant force change law of the blasting side can be directly observed: the blasting vibration acceleration at the arch foot is the largest, indicating that the outside lining of the blasting side has the largest detonation force.

Conclusions
Through the combination of field test and numerical simulation, this paper analyzes the dynamic response of a deep-buried tunnel with small interval under blasting load, and the specific conclusions are as follows: 1.
Through the spectrum analysis of the field measured waveform and the nonlinear regression of the data, the empirical formula conforming to the attenuation law of the Sadowski vibration velocity of the project is obtained.2.
On the inside of the lining (the side contacting the tunnel headroom), the vibration continues from the arch foot along the tunnel lining towards the vault and tunnel bottom.On the outside of the lining (the side contacting the surrounding rock), the lining on the blasting side is affected by vibration before the back side at the initial stage of blasting.The vibration velocity in the order of large to small is arch foot, tunnel hance, and vault.

3.
After the blasting action is completed, the vibration velocity of the lining structure gradually decreases but will not disappear, and will only approach 0. The vibration velocity of the cross-section shows little change.The vibration velocity of the section near the detonation source attenuates faster, while the vibration velocity of the section far away from the detonation source attenuates slower and there is little difference in the attenuation velocity.

4.
The influence of the horizontal displacement of the tunnel lining is greater than that of the vertical displacement, and the randomness of the vertical displacement is larger.The influence of the inner displacement of the lining is greater than that of the outer displacement.5.
The vibration acceleration at the arch foot of the section is the largest, followed by that at the tunnel hance, and the peak value of vibration acceleration at the vault is the smallest.

Figure 1 .
Figure 1.Location of Badaling Great Wall Station.(a) Location in China.(b) Surrounding environment.(c) The relative location of the Great Wall Station.

Figure 1 .
Figure 1.Location of Badaling Great Wall Station.(a) Location in China.(b) Surrounding environment.(c) The relative location of the Great Wall Station.

Figure 3 .
Figure 3. Section size and spacing of separation section of three holes (unit: cm).

Buildings 2023 , 26 Figure 1 .
Figure 1.Location of Badaling Great Wall Station.(a) Location in China.(b) Surrounding environment.(c) The relative location of the Great Wall Station.

Figure 3 .
Figure 3. Section size and spacing of separation section of three holes (unit: cm).

Figure 3 .
Figure 3. Section size and spacing of separation section of three holes (unit: cm).

3. 2 .
Experimental Scheme 3.2.1.Design of Blasting Parameters for the Standard Section of Three Holes Separation

Figure 4 .
Figure 4.The layout of the upper step gun hole of grade Ⅲ surrounding rock in the middle of the three-hole separation section (unit: mm).

Figure 4 .
Figure 4.The layout of the upper step gun hole of grade III surrounding rock in the middle of the three-hole separation section (unit: mm).

Figure 4 .
Figure 4.The layout of the upper step gun hole of grade Ⅲ surrounding rock in the middle of the three-hole separation section (unit: mm).

Figure 5 .
Figure 5.The layout of the lower step gun hole of grade Ⅲ surrounding rock in the middle of the three-hole separation section (unit: mm).

Figure 6 .
Figure 6.Diagram of arrangement of vibration measuring points.

Figure 5 .
Figure 5.The layout of the lower step gun hole of grade III surrounding rock in the middle of the three-hole separation section (unit: mm).

3. 2 . 2 .
Vibration Monitoring Scheme A total of seven measuring points are arranged for the vibration test, and their positions are shown in Figure 6.The mileage corresponding to the face of the middle hole is DK68+158, and the thickness of the wall is 6 m.Measuring point 2 is on the same plane as the charging section on the tunnel face.Measuring point 1 is set in the excavated stage of the middle tunnel, 5 m away from measuring point 2. Measuring points 3, 4, and 5 are set in the unexploded area of the middle tunnel.The distances from measuring point 2 are 5 m, 13 m and 23 m.Measuring points 6 and 7 are on the same plane with measuring point 2, but 0.8 m lower than that of measuring point 2. By monitoring the vibration velocity of the measuring point in the same plane and a measuring line, the propagation law of blasting vibrations in different positions is analyzed.

Figure 4 .
Figure 4.The layout of the upper step gun hole of grade Ⅲ surrounding rock in the middle of the three-hole separation section (unit: mm).

Figure 5 .
Figure 5.The layout of the lower step gun hole of grade Ⅲ surrounding rock in the middle of the three-hole separation section (unit: mm).

3. 2 . 2 .
Vibration Monitoring Scheme A total of seven measuring points are arranged for the vibration test, and their positions are shown in Figure 6.The mileage corresponding to the face of the middle hole is DK68+158, and the thickness of the wall is 6 m.Measuring point 2 is on the same plane as the charging section on the tunnel face.Measuring point 1 is set in the excavated stage of the middle tunnel, 5 m away from measuring point 2. Measuring points 3, 4, and 5 are set in the unexploded area of the middle tunnel.The distances from measuring point 2 are 5 m, 13 m and 23 m.Measuring points 6 and 7 are on the same plane with measuring point 2, but 0.8 m lower than that of measuring point 2. By monitoring the vibration velocity of the measuring point in the same plane and a measuring line, the propagation law of blasting vibrations in different positions is analyzed.

Figure 6 .
Figure 6.Diagram of arrangement of vibration measuring points.

Figure 6 .
Figure 6.Diagram of arrangement of vibration measuring points.

Figure 7 .
Figure 7. Velocity vibration waveform of measurement point 2. (a) Vertical velocity vibr waveform of measuring point 2. (b) Radial velocity vibration waveform of measuring point Tangential velocity vibration waveform of measuring point 2.

Figure 7 .
Figure 7. Velocity vibration waveform of measurement point 2. (a) Vertical velocity vibration waveform of measuring point 2. (b) Radial velocity vibration waveform of measuring point 2. (c) Tangential velocity vibration waveform of measuring point 2.

Buildings 2023 , 2 Figure 12 .
Figure 12.Stress distribution in the X direction of the three-way tunnel under the influence of dead weight.

Figure 13 .
Figure 13.Stress distribution in the Y direction of the three-way tunnel under the influence of dead weight.

Figure 14 .
Figure 14.Stress distribution in the Z direction of the three-way tunnel under the influence of dead weight.

Figure 12 . 26 Figure 12 .
Figure 12.Stress distribution in the X direction of the three-way tunnel under the influence of dead weight.

Figure 13 .
Figure 13.Stress distribution in the Y direction of the three-way tunnel under the influence of dead weight.

Figure 14 .
Figure 14.Stress distribution in the Z direction of the three-way tunnel under the influence of dead weight.

Figure 13 . 26 Figure 12 .
Figure 13.Stress distribution in the Y direction of the three-way tunnel under the influence of dead weight.

Figure 13 .
Figure 13.Stress distribution in the Y direction of the three-way tunnel under the influence of dead weight.

Figure 14 .
Figure 14.Stress distribution in the Z direction of the three-way tunnel under the influence of dead weight.

Figure 14 .
Figure 14.Stress distribution in the Z direction of the three-way tunnel under the influence of dead weight.

Figure 15 .
Figure 15.Stress intensity distribution of three-way tunnel under the influence of dead weight.

Figure 16 .
Figure 16.Y-direction displacement distribution cloud map of the three-way tunnel under the influence of dead weight.

Figure 15 . 2 Figure 15 .
Figure 15.Stress intensity distribution of three-way tunnel under the influence of dead weight.

Figure 16 .
Figure 16.Y-direction displacement distribution cloud map of the three-way tunnel under th influence of dead weight.

Figure 16 .
Figure 16.Y-direction displacement distribution cloud map of the three-way tunnel under the influence of dead weight.

Figure 17 Figure 17 .
Figure17shows the nephogram of the vibration velocity corresponding to 1 ms, 3 ms, 5 ms, and 7 ms on the inner side of the lining.

Figure 18
Figure18shows the nephogram of vibration velocity corresponding to 1 ms, 3 ms, 5 ms, and 7 ms on the external side of the lining.

( 4 )
Analysis of longitudinal vibration velocity distributionThe comparison of the vibration velocity and vibration waveform of each longitudinal measurement point in the left tunnel is shown in Figures21 and 22 .

( 4 )
Analysis of longitudinal vibration velocity distributionThe comparison of the vibration velocity and vibration waveform of each longitudinal measurement point in the left tunnel is shown in Figures21 and 22 .

Figure 21 .
Figure 21.Diagram of longitudinal vibration velocity change of measuring points.

Figure 21 .
Figure 21.Diagram of longitudinal vibration velocity change of measuring points.

( 4 )
Analysis of longitudinal vibration velocity distributionThe comparison of the vibration velocity and vibration waveform of each longitudinal measurement point in the left tunnel is shown in Figures21 and 22 .

Figure 21 .
Figure 21.Diagram of longitudinal vibration velocity change of measuring points.

Table 1 .
Field measured data graph.

Table 6 .
III level surrounding rock parameters.

Table 7 .
Concrete parameter table.In terms of the damping of rock mass, Rayleigh damping is the most widely used in dynamic analysis and practical engineering.The linear combination of global mass matrix [M] and global stiffness matrix [K] represent global damping matrix [C].

Table 8 .
Comparison between the calculated vibration velocity and the measured vibration velocity.

Table 8 .
Comparison between the calculated vibration velocity and the measured vibration velocity.

Table 8 .
Comparison between the calculated vibration velocity and the measured vibration velocity.

Table 9 .
Vibration velocity of measuring point in section.

Table 9 .
Vibration velocity of measuring point in section.