Dynamic Response Characteristics of the Structures under Excavation by High-Efficiency Blasting in the Urban Ultradeep Foundation Pit

To ensure the safe and eﬃcient excavation of ultradeep foundation pits in a complex urban environment, the ultradeep foundation pit excavation project of Liuguangmen Urban Complex in Guiyang City was taken as the study point. A high-eﬃciency blasting method was proposed for the excavation of vibration-isolating slot and electronic detonators, and a three-dimensional spatial calculation model of the foundation pit structure was established. In addition, the ﬁeld test and numerical simulation of the blasting vibration were developed. The feasibility of the high-eﬃciency blasting method was demonstrated, and the propagation law of the dynamic response characteristic parameters of the structure was explored using an electronic detonator and vibration-isolating slot. The results show that the electronic detonator carried out peak shaving within the group, unloading waves between the groups, and the blocking eﬀect of the vibration-isolating slot, and all these eﬀectively reduced the peak stress and peak particle vibration velocity, evenly distributed the vibration velocity, and guaranteed the safety of the surrounding existing buildings, proving the feasibility of the eﬃcient blasting method. There is a “wall eﬀect” at the supporting pile, whose propagation velocity is lower than that inside the rock mass, and the PPV is radial X > tangential Y > vertical Z . The peak values of stress, strain, and PPV all decrease with increasing supporting pile height, and there is an inﬂection point at 20m. The internal nodes of the pile reciprocate without sharp change and failure in displacement. The peak displacement of diﬀerent nodes is in the order X > Y > Z with increasing pile height, and the maximum peak value is 0.08mm.


Introduction
e safe and efficient excavation of foundation pits in a complex urban environment has been investigated in detail. For example, Li et al. [1] adopted a vibration-isolating hole combined with long and short time delays of an electronic detonator on the basis of controlled dosage to solve the problem of excavation of large quantity and deep rock foundation pit with complex strata in downtown area. Gao et al. [2] analyzed the dynamic response law and vibration velocity safety threshold of an adjacent sewage pipeline under the blasting action of foundation pit excavation. Ma et al. [3] carried out a shaking table model test to study the seismic stability of the retaining structure of underground continuous wall during the construction of deep foundation pit. Zhang et al. [4] studied the effect of different traffic load coefficients and supporting structure parameters on the stability of the supporting structure of deep foundation pit under construction. Xue et al. [5] discussed the vibration response of train load to the retaining structure of foundation pit and proposed the measures of velocity limit and vibration reduction. Xie et al. [6] conducted vibrational, acoustic, and numerical simulation tests on the vibration velocity threshold of foundation pit blasting. Dowding et al. [7] studied the strain response characteristics of the urban structure under blasting excavation. Yankelevsky et al. [8] studied the blasting resistance and response of reinforced concrete building structures under the terrorist threats of blasting in Israeli urban areas. Meins et al. [9] summarized the structural response of high-rise buildings under shortdistance blasting in New York City, providing a reference model for the development of the standard for short-distance blasting in high-rise buildings. Xie et al. [10] optimized the shallow-buried tunnel blasting scheme based on the complex urban environment and controlled the vibration velocity below 1.5 cm/s. Wang et al. [11] investigated ways to reduce vibration and monitor the foundation pit blasting to protect surrounding buildings. Zhou et al. [12] studied the dynamic response, damage evolution, and failure forms of BFRP bars in shallow-buried urban tunnels under repeated explosions. Zhang et al. [13] analyzed the spatial distribution characteristics of blasting vibration velocity and dynamic stress of high-rise buildings through dynamic numerical analysis. Wang et al. [14] captured the acceleration spectrum of Wenchuan earthquake to establish a reasonable threedimensional seismic response model, simulating the dynamic stability, instability abrupt transition process, and time effect characteristics of the foundation pit and support structure under strong earthquake. Wei et al. [15] simulated the vibration velocity response of different structures by the dynamic analysis method and verified that the vibrationisolating slot is an effective damping method. Sun et al. [16] experimentally investigated the damping effect of the vibration damping holes in the foundation pit, analyzed the vibration velocity attenuation law and the vibrationdamping rate of the single-and double-row vibrationdamping holes, and calculated and verified them by numerical simulation. However, there are a few reports on the evolution law of the dynamic response characteristic parameters under high-efficiency blasting effect of excavation vibration-isolating slot and electronic detonator initiation.
In view of this, based on the excavation of ultradeep foundation pits (the excavation depth exceeds 15 m) in the complicated environment around Liuguangmen in Guiyang city, combined with the vibration reduction excavation technology of electronic detonator blasting and the characteristics of the vibration-isolating slot, in this study, a high-efficiency blasting excavation and slope support scheme for foundation pits and blasting vibration monitoring scheme were formulated. A digital network vibration meter (CBSD-VM-M01) was used to conduct on-site blasting vibration test and analyze the vibration response characteristics of existing buildings according to the test results. en, based on the actual blasting site, a three-dimensional spatial calculation model of the foundation pit structure was established to study the dynamic response law of the supporting structure under the blasting excavation of the ultradeep foundation pit.

Test Site.
e experimental site is the foundation of an urban complex engineering foundation pit, located in the political and cultural center of Guiyang city, close to the provincial government building. e width of the foundation pit is about 120 m, the maximum section length is more than 260 m, and the maximum excavation depth is up to 42 m. e underground part of the proposed building consists of two floors of sewage treatment structure, three floors of the underground parking lot, and one floor of commercial district from bottom to top. Figure 1 shows the complex surrounding environment of the foundation pit excavation. Buildings (structures) close to the foundation pit include residential buildings at 22 m in the north, light rail stations under construction at 20 m in the west, and specialized hospitals at 55 m.

Blasting Excavation Scheme.
e upper layer of the rock mass is weathered, while the lower layer is relatively hard, and fissure karst is developed. Blasting was used to break the rock in the area where conditions permit. e slope of the foundation pit is vertical, and an anchor cable row pile with a diameter of 1800 mm and a pile center distance of 1500 mm apart was used for support. Part of the adjacent pile rock mass was retained during the earth-rock excavation. After the pile body reaches the design strength and the anchor cable is gradually applied to prestress, the retained rock mass was broken and cleared mechanically.
In terms of the adverse effects of blasting of this project, vibration, the main control object, was mainly controlled by the single-hole charge, the initiation network, and the vibration propagation path. According to Sadov's empirical formula for vibration calculation [17] and similar construction experience, the blasting area and the mechanical broken area (rock mass within 15 m of the adjacent pile) are divided and the single-hole charge is strictly controlled. Table 1 shows the hole pattern parameters of the blasting area. e effect of the Wenchuan earthquake acceleration on steel-supported foundation pit was studied following literature procedure [18]. e instability process is divided into three stages: reciprocating motion and a sharp increase in displacement and instability with the corresponding periods of t ≤ 10 s, 10 ≤ t ≤ 20 s, and t ≥ 20 s, respectively. e design of the detonation network ensures that the blasting time of the hole is independent of each other, avoids vibration superposition, and controls the single duration of vibration. Based on the convenience of digital electronic detonator delay setting, the detonation network of the electronic detonator was designed to have a long and short delay, namely, the blastholes in the network are grouped (one group has 5-6 holes), the holes were set as short delay (45 ms) within the group and long delay (1500 ms) between the groups. e delay network is shown in Figure 2. A vibration-isolating slot is a shallow ditch dug in the vicinity of buildings and other facilities to prevent vibration harm. To ensure the safety of the surrounding structures, the slot was set with a depth of 4 m, a width at the top of 4 m, and a width at the bottom of 3 m.

Monitoring Design.
A huge amount of energy is released instantly in the process of blasting as a result of the detonation of the explosive. A part of this energy is used for the crushing of the mining rock and the compaction of the soil, rock, and mud, while the other part is absorbed by the surrounding rock or soil rock and mud. In addition, the released energy spreads in all the directions in the form of an elastic vibration wave, causing vibration in the adjacent supporting structure and the existing building particles. To reveal the vibration response characteristics of the structure, a digital network vibration meter (CBSD-VM-M01) was used for the field test. A total of two vibration monitoring points were arranged, in which the measuring points 1 and 2 were arranged in the existing residence and the support structure, respectively, as shown in Figure 3.

Shock and Vibration
According to the monitoring scheme, the dirt on the residential foundation surface and the supporting pile wall surface where the measuring point is located was cleaned.
en, adjust the X direction (radial direction) of the intelligent sensor of the vibration meter to keep consistent with the direction of the explosion source, connect the sensor with the surface of the ground (or support pile wall) with gypsum (or super glue and fixture), and make the instrument level center, while ensuring the fastness of connection. Finally, set the parameters of the vibration meter through the control analyzer and press "Collect" on the intelligent sensor to conduct the vibration data monitoring test. Figure 4 shows the test system for the blasting vibration. e system mainly consists of a control analyzer, an intelligent sensor, a router, and a computer. When blasting was carried out, the intelligent sensor was fixed together with the existing building. e blasting action causes the surrounding building particles to vibrate, and the sensor housing vibrates accordingly. Under the action of inertia, the internal coil, damping ring, and core rod move relative to the housing, thus the coil generates induced electromotive force, which is converted into vibrational velocity. e control analyzer connects the intelligent sensor through the router to set vibration parameters and obtain vibration data and finally uploads the computer terminal for analysis and processing. According to the site environmental conditions and blasting vibration characteristics, the trigger level of the test device was set as 20%, the sampling frequency as 2 kHz, the recording duration as 6 s, and the delay as 0.2 s.

Analysis of Test Results of the Structural
Vibration Response Characteristics Figure 5 shows the characteristic curve of blasting vibration response of seismic wave extending to the existing buildings under the action of vibration-isolating slot and electronic detonator initiation. e vibration waveform of the first group (0-180 ms), the second group (1680-1905 ms), the third group (3405-3630 ms), and the fourth group (5130-5355 ms) was monitored during 0-1 s, 1.5-2.5 s, 3-4 s, and 5-6 s, respectively. In all the four groups, the PPV at the monitoring point of the existing houses was 0.210 cm/s, and the supporting pile was 0.195 cm/s. e duration of each group was roughly equal. e PPV and main frequency of each group are shown in Table 2. Figure 5 shows the tested vibration waveform of the residence and the supporting pile, indicating that the stress transfer process caused by the blasting of the hole between groups is basically independent of each other. e vibrational velocity of each group was relatively uniform, the wave peak and trough of the vibration velocity between the groups were superimposed to reduce, and the staggering peak reduced the peak of the vibration velocity. ere is a long time interval between the groups. After blasting the last hole of the previous group, the energy is not supplemented and the vibration velocity gradually decreases. In this way, the peak shaving within the group and unloading waves between the groups reduce the vibration.
Vibration is one of the adverse effects of blasting and is the main factor affecting the safety of construction and the most common cause of "disturbance in the surrounding structures and disturbing the people". erefore, vibration must be checked and analyzed strictly. In Table 2 of the tested signal data, the PPV and main frequency of the residence and the supporting piles as a whole show the radial direction X > tangential Y > vertical Z with their PPV and main frequency corresponding to (0.210 cm/s, 21.27 Hz) and (0.195 cm/s, 18.35 Hz), respectively. Referring to relevant national regulations and studies on the vibration velocity threshold for disturbing residents [6], the threshold of this type of house was taken as 1.5 cm/s to minimize disturbance to the residents. e impact of this blasting stress disturbance on the residential building in this study met relevant regulations and is lower than the threshold of the nuisance security in relevant studies and regulations. e PPV of radial, tangential, and vertical at the monitoring point (in the support pile) 15 m from the horizontal distance is 0.195, 0.148, and 0.029 cm/s, which are lower than that at the measuring point (in the residence), 35 m (0.210 cm/s in radial, 0.203 cm/s in tangential, and 0.041 cm/ s in vertical), showing abnormalities. e reason is that the tested data of the residence are the internal propagation of the rock mass, while there is a wall surface between the supporting pile and the rock mass, resulting in the "wall effect," with propagation velocity lower than the internal velocity of the rock mass. In other words, the stress wave propagates, reflects, and stretches on the wall surface of the supporting pile, which leads to energy dissipation and thus reduces the propagation velocity.

Numerical Analysis of Dynamic Response
Characteristics of the Structure  foundation pit blasting excavation. e fluid-solid coupling model was selected, and the Lagrange algorithm was used for rock and reinforced concrete. Considering that the material under the impact load of the explosion is extremely deformed, the ALE algorithm is set as the explosive. e numerical model adopts the SOLID164 solid element and was modeled in proportion according to the design of the blasting scheme. To reduce the calculation workload, 11 support piles were selected for modeling, with a pile diameter of 1.8 m and a pile center distance of 1.5 m, and the crown beam structure was established at the top. e reserved rock mass is 8 m wide at the top and 10 m wide at the bottom, and the height difference from the blasting excavation surface is 8 m. e vibration-isolating slot is 4 m deep, 4 m wide at the top, and 3 m wide at the bottom. Because of the independent stress transfer process between the groups in the detonating network, a set of holes (the number of holes is 6) in the network was selected for modeling, and the parameters used are as follows: (i) the difference between the holes was set as 45 ms, (ii) the equivalent volume in the holes was set as no. 2 rock emulsion explosive of 6 kg, (iii) the hole spacing was 3.0 m, (iv) the row spacing was 2.5 m, and (v) the termination time was set as 450 ms. e relative position and mesh division of the model are shown in Figure 6.

Model Parameters.
e model materials are pile and crown beam reinforced concrete materials, geotechnical materials, explosive materials, stemming materials, and air materials.  ere are some difficulties in modeling the processing of reinforced concrete materials for supporting piles. To facilitate modeling and simplify the calculation, the equivalent modulus treatment method was selected for modeling of reinforced concrete materials by referring to relevant studies on elastoplasticity constitutive models of reinforced concrete [19,20]; that is, according to the fraction of reinforcement in section and concrete strength, the reinforcement in reinforced concrete materials was dispersed among the materials for homogeneity.
According to relevant specifications, the concrete material is C30 with a Young's modulus of 30 GPa and a Poisson's ratio of 0.2, the fraction of reinforcement in the section is 2%, and Young's modulus of reinforcement is 200 GPa.
During the modeling process, the * MAT_BRITTLE_DAMAGE_(TITLE) (096) model was selected, and relevant parameter settings are listed in Table 3.
(2) Geotechnical Material Model. e geotechnical material itself is nonhomogeneous and is also filled with structural planes such as joints and cracks. To simplify the calculation, it was idealized into homogeneous and continuous elasticplastic materials in the simulation process. In this study, the * MAT_PLASTIC_KINEMATIC_(TITLE) material model was used to simulate the geotechnical material, and its mechanical parameters are listed in Table 4.
(3) Explosive Material Model. e no. 2 rock emulsion explosive was used in blasting construction, and the maximum single-hole charge in the simulation process is 6 kg. e * MAT_HIGH_EXPLOVE_BURN explosive material model with a mass density of 1200 kg/m was used. Its detonation velocity is 4500 m/s. e detailed parameters are shown in Table 5. According to the nature of rock mass and its stress state in the process of the explosion, the JWL equation of state was used to simulate the relationship between the pressure and specific volume in the process of explosive detonation, namely, where A, B, R 1 , R 2 , and ω are the material constants; P is the pressure from the equation of state; V is the relative volume; E 0 is the internal energy density per unit initial volume.
(4) Stemming Material Model. e rock powder and mud from the drilling hole were used for on-site stemming. e material model of * MAT_SOIL_AND_FOAM was selected. e main parameters are shown in Table 6.
Because of the complexity of the building structure, the structure was no longer modeled separately to simplify the numerical calculation process, and the vibration velocity waveform of the bottom node near the foundation pit side of residential building (actual monitoring point of the project) was taken for analysis.

Verification of the Numerical Results.
In the numerical modeling, referring to the actual engineering environment, the residential building on the north side of the adjacent foundation pit was taken as the research object, and the vibration waveform of the bottom node (the actual monitoring point) of the residential building near the foundation pit was verified. e tested and the simulated vibration waveform are plotted in Figure 7, indicating that the trend of the curve for the tested and simulated vibration velocity versus time basically coincides, and the peak value of the simulated vibration velocity is 0.350 cm/s, while the tested peak value is 0.277 cm/s. e value of numerical simulation is high, and the relative error with the tested value is 20.8% because the modeling process idealizes the size structure of the vibration-isolating slot, geotechnical medium, surrounding topography, and explosive state equation. Although there is a certain difference between the simulated and the tested value, the numerical result data within the difference range are reasonable and feasible.

e Numerical Effect.
According to the delay time of blasthole initiation, the von Mises stress cloud diagram after the blasting of no. 1, no. 3, no. 4, and no. 6 blastholes was taken, as shown in Figure 8, indicating that, in the blasting process of a single group of blastholes, as the number of blastholes increases, the stress concentration range in the blasting area increases, and the stress concentration in the blasting area gradually increases, reflecting the rockbreaking process to a certain extent. According to the stress nephogram, the rock mass in the left side of the vibration-isolating slot is basically blue with small stress, indicating that the slot plays a blocking role in the transfer of blasting stress and effectively protects the foundation pit supporting system and existing buildings on the left side of the slot. To quantify the barrier effect of the vibration-isolating slot on the blasting stress, the stress versus time curve of the characteristic points 10 m from the center of the detonation zone to the direction of the vibration-isolating slot and the opposite direction was taken, as shown in Figure 9. e peak stress of the characteristic points on the side without the vibration-isolating slot and that on the side with the slot were extracted, and the corresponding peak stress was found as 8.45 and 1.48 MPa, respectively. e peak stress of the characteristic points on the side of the vibration-isolating slot was only 17.6% of that on the other side, proving that the vibration-isolating slot is an economical and efficient engineering measure for the protection of the structure near the area.

Dynamic Response Characteristic Analysis.
e most unfavorable factor of blasting excavation for the foundation pit is the disturbance caused by blasting stress to the supporting structure. e safety and stability of the structure depend on its dynamic response under blasting disturbance. erefore, the dynamic response law of displacement, stress, strain, and velocity was investigated to provide reference for vibration control and prediction. e response characteristics of displacement, velocity, stress, and strain of the supporting pile were studied starting from the center of the top of the pile and taking the nodes and corresponding units of the supporting pile axis vertically downward every 10 m. e node numbers are 846585, 848037, 848169, 849554, and      Shock and Vibration 7 849915, and the unit numbers are H1545060, H1545492, H1546434, H1546467, and H1547250, as shown in Figure 10.

Displacement Response Characteristics of the Sup-
porting Pile. Figure 11 shows the displacement versus time response curves of the nodes at different positions of the supporting pile axis. Figure 11 shows that the displacement of the axial node of the supporting pile is large in the X direction, followed by in the Y direction, and smallest in the Z direction and is related to the position applied by blasting stress. e overall displacement of the node reflects the movement of the pile. When the displacement is about 300 ms, the node itself shows reciprocating movement, and the displacement does not change dramatically during the reciprocating process. e node in the numerical calculation process did not show any failure, indicating that there is no damage inside the pile body under the stress disturbance.

Shock and Vibration
In the calculated step size, the displacement of no. 849915 node, which is closest to the explosion zone, is the maximum, reaching 0.081 mm in the X direction. e displacement of the bottom three characteristic points (node no. 849915, no. 849554, and no. 848169) in the X direction showed a recovery trend. Each node displacement in the X direction versus time response curves is shown in Figure 11(f ). e bottom node of the supporting pile moves directly in the negative X direction after being disturbed by the blasting stress, while the top node first moves slightly in the positive X direction and then moves in the negative X direction.
is reflects the movement process of the rod body after being stressed and is following the motion law of the elastic rod body.

Velocity Response Characteristics of the Supporting Pile.
e response trend of vertical, horizontal radial, and tangential PPV with increasing supporting pile height under the action of the electronic detonator initiation and the vibration-isolating slot is statistically shown in Figure 12, indicating that the PPV is radial X > tangential Y > vertical Z and is consistent with the experimental outcome. At the bottom of the pile, the radial, tangential, and vertical PPVs are the largest, with 4.55, 2.83, and 0.56 cm/s, respectively, whereas at the top, they are small with the values of 0.17, 0.11, and 0.03 cm/s, respectively. e PPV of the blasting response in each direction attenuates strictly with increasing supporting pile height, and there is no "elevation amplification effect." It can be considered that the existence of the vibration-  isolating slot has an obvious effect on seismic wave propagation. e radial and the tangential peak vibration velocity attenuate sharply along with the pile in the range of 0-20 m.
e inflection point appears at 20 m, and the attenuation tends to be gentle in the range of 20-40 m. e vertical PPV is flat with increasing pile height and is generally gentle. Figure 13 shows the response trend of tensile stress, effective stress, and shear stress with increasing supporting pile height under the action of the electronic detonator initiation and the vibration-isolating slot. Clearly, the variation trend of tensile stress, effective stress, and shear stress peak is consistent and decreases with increasing pile height, and the peak stress at the bottom of the pile reaches the maximum value of 3.58, 3.48, and 1.95 MPa, respectively, and the peak stress is the smallest at the top with 0.12 MPa, 0.14 MPa, and 0.08 MPa, respectively. According to the strength theory, all the stresses are less than 4 MPa and the supporting piles are in a safe state. e height along the pile is in the range of 0-20 m in the following order: peak tensile stress > effective stress > shear stress. e height along the pile is in the range of 20-40 m, in the following order: peak effective stress > tensile stress > shear stress. Figure 14 shows the strain response trend of vertical, horizontal radial, and tangential directions with increasing supporting pile height under the action of the electronic detonator initiation and the vibration-isolating slot, indicating that the variation trend of the vertical, tangential, and radial strain peaks is consistent, and they all decrease with increasing pile height. e maximum strain peaks at the bottom of the pile are 2.5, 4.8, and 8.4, respectively, and at the top is the smallest, at 0.08, 0.3, and 0.1, respectively. e height along the pile is in the range of 0-20 m, and the strain follows the order horizontal > vertical. e height along the pile is in the range of 20-40 m with the strain value in the following order: tangential > vertical > radial.

Conclusion
Based on the ultradeep foundation pit excavation in the complex environment around Liuguangmen, Guiyang city, an effective blasting model of vibration-isolating slot and electronic detonator was successfully established, the vibration response characteristic test and dynamic response numerical test were carried out, and the dynamic response characteristics of existing structures under the blasting action of 6 kg single-hole charge were studied. e outcome of this study is as follows: (1) e vibration velocity and stress of the structure are in a safe range under the action of excavation vibration-isolating slot and electronic detonator initiation. In addition, the efficient blasting excavation scheme was found to be feasible, ensuring the safety of the surrounding existing buildings. (2) e detonation of the electronic detonator carries out peak shaving within the group, unloading waves between the groups, and the blocking effect of the vibration-isolating slot, effectively reducing the peak stress and PPV, evenly distributing the vibration velocity, and ensuring the safety and efficiency of ultradeep foundation pit excavation in the complex urban environment. ere is a contact surface between the supporting pile and the rock mass, and the stress wave in the surface is reflected and stretched, resulting in the "wall effect", and the propagation velocity is lower than that inside the rock mass. e PPV is in the following order: radial X > tangential Y > vertical Z.  Shock and Vibration 11 (3) Under the action of the electronic detonator initiation and vibration-isolating slot, the peak values of stress, strain, and PPV all show a decreasing trend with increasing supporting pile height, and there is an inflection point at 20 m. In the 0-20 m range, the peak stress is tensile > effective > shear, the peak strain is horizontal > vertical, and the PPV attenuates rapidly. In 20-40 m range, the peak stress value is effective > tensile > shear, the peak strain value is in the order tangential > vertical > radial, and the PPV attenuates gently. (4) e internal nodes of the pile reciprocate without sharp change and failure in displacement. e peak displacement of different nodes with increasing pile height is in the following order: X > Y > Z. e maximum displacement value of the axial node of the pile within the calculated step size (45 ms) is 0.08 mm, which is located at the bottom of the pile, and its motion direction is opposite to the explosion zone with a trend of recovery, while the displacement is located at the top of the pile and its motion direction is first toward the explosion zone and then the reverse zone, following the motion law of the elastic rod body.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare no conflicts of interest.