Abstract

The subway station is easy to be attacked by terrorist bombings, and it will cause heavy casualties. In this paper, a comprehensive casualty assessment method for personnel in the subway structure was established based on the existing personnel injury model. The spatial distribution characteristics of the shock wave suffered by the personnel in the subway platform were obtained. Combined with the comprehensive casualty assessment method, the personnel casualty area for the explosions in the subway platform was divided. The results show that for the same explosive charge, the maximum positive phase impulse generated by the explosion at the edge of the platform is smallest. The “notch effect” for the stair exit will increase the shock wave load. When the explosive is exploded in the center of the platform, the smaller the explosive charge is, the more obvious the “notch effect” is. When the explosive charge reaches 40 kg, the personnel safety area is reduced to a certain extent behind the stair except for the explosion happening at the stair. Also, the higher the shield door is, the larger the safety area behind the stair is.

1. Introduction

In recent years, the subway as a crowded public space has gradually become the preferred target for terrorist attacks. Because of the closed underground structure adopted in the subway system, the propagation law of the shock wave generated by the explosion in the subway structure is much more complicated than that in the free atmospheric space. Hence, the casualties caused by explosion are completely different from those in the unconstraint space. So far, many efforts have been devoted to the propagation law of explosive shock waves in the structure and the dynamic response of the structure for the explosion in the underground structures such as subway. Rigas and Sklavounos [1] calculated the propagation process of shock wave for the explosion in complex constrained space structures with CFX. The Italian European Structural Assessment Laboratory (ELSA) [2] used the finite element software to simulate the propagation of shock waves and the damage of carriage structure for the explosion in subway trains; Wu et al. [3] proposed a test procedure for the phenomenon of large-scale underground explosion and studied the propagation characteristics of stress wave in soil medium. Li et al. [4] compiled an analytical procedure for the dynamic response of double-track subway tunnel under unilateral tunnel explosion load by using FORTRAN program. Choi et al. [5] obtained the change of the effective strain of the concrete lining structure of the tunnel with the parameters such as explosive charge, explosion distance, and tunnel diameter by numerical simulation of the explosion.

In fact, the explosive charge used in the terrorist incidents in the subway is only few kilograms of TNT equivalent. The station structure cannot be completely damaged in this case but can cause the casualties. So, it is a key issue to reasonably assess the impact of personnel injury caused by explosions in the subway structure. Experimental and theoretical research has been conducted on the evaluation criteria of personnel injury by the effect of explosion. For example, “TM5-1300 Antiaccidental Explosive Design Manual” of the United States [6] proposed injury standard of shock wave with duration of 3∼5 ms to the human body. Li [7] gave the overpressure threshold corresponding to different casualty levels through a large number of animal tests combined with the relevant standards of the United States and the Soviet Union. Li et al. [8] analyzed the casualties of rabbits at different distances from the explosion point through the free-field explosion test and obtained the safety threshold of the shock wave for rabbits and the safety distance for different charges. Furthermore, Pietersen [9] obtained the harm equivalent and the formula of injury probability caused by the overpressure. Ferradás et al. [10] obtained the method of directly analyzing the impact of the shock wave to personnel injury on the existing characteristic curve of shock wave overpressure-impulse-distance.

Compared with many shock wave injury criteria, it is found that there are some differences in shock wave injury criteria proposed due to different test conditions. In this paper, the personnel injury model suitable for the explosion inside the subway structure is applied to conduct the in-depth study.

2. Method for Assessing Personnel Injury by Blast Loading in Subway

The severity of shock wave injury to personnel is usually measured by overpressure and positive phase impulse. The shock wave generated by the explosion in the subway structure is a complex shock wave [11], which can produce multiple reflection peaks. The maximum peak overpressure and positive phase impulse in the simulation are evaluated aswhere is the maximum peak overpressure, Pa; is the standard atmospheric pressure, Pa; and is the entire simulated time, s.where is the positive phase impulse of shock wave, Pa·s, and is the action time of positive pressure of shock wave, s.

The probit model is commonly used for the casualty assessment. There is a certain logarithmic relationship between probit units and reaction injury factors [10, 12], and the general formula is as follows:where is the probit units, is the injury variable of a certain injury event, and A and B are the parameters of the probit unit model.

Each individual injury caused by explosion can be evaluated by a corresponding probit unit model. The lung injury model is based on overpressure and impulse given by Bowen et al. [13], but this model, developed from free-field blasts, is not applicable to explosions in the subway. Therefore, in this paper, the adjusted severity of injury index (ASII) scoring method, applicable to complex shock wave, is used for lung injury assessment [14, 15], and the mathematical model is as follows:where is the maximum inward chest wall velocity, m/s, and is the equivalent triangular pulse duration, s, and calculated as

The relationship between injury degree and ASII score is shown in Table 1 [14].

The assessment model of the eardrum injury caused by the direct action of shock waves is shown in the following equation [11]:

The throwing injury caused by the indirect action of the shock wave can be divided into the whole body impact injury and the head impact injury. The injury assessment models are as shown in equations (7)∼(8) [16].

The whole body impact injury probability model is as follows:

The head impact injury probability model is as follows:

The probability of occurrence can be obtained from the probit unit and can calculated as [17]where P is the probability of occurrence; when Y is taken as 5, the probability of occurrence is 50%.

In this paper, the injury consequences of human body are divided into three types: death, serious injury, and light injury [18]. Combined with the results of individual injury assessment, the comprehensive personnel casualty assessment process for the explosion in the subway structure is obtained, as shown in Figure 1.

3. Finite Element Model of Explosion in Subway Station

3.1. Size of the Subway Platform

A typical station in Beijing was chosen as the research object, and it is an island-style and double-column structure. The effective length and the width of the platform are 120 m and 13 m, respectively. The height at the track and platform is 6 m and 4.8 m, respectively. The radius of platform column is 0.5 m. Longitudinal distance and horizontal distance between the columns are 6 m and 7 m, respectively. There is a staircase on the platform from 30 m to the center of the platform, respectively. Also, the width of each stair is 5 m, with the longitudinal bottom 13 m long and the short side 4 m long. There is an exit with 8 m long and 5 m wide at the top of the stair. The cross section of the subway platform is shown in Figure 2, and the half of longitudinal section of the subway platform is shown in Figure 3.

3.2. Structure of the Shield Door

At present, there are three main types of shield doors installed in the subway: full-height and fully closed shield doors, full-height and semiclosed shield doors, and half-height shield doors. The total height of the full-height and fully closed shield door is 3.2 m, whose shield door glass is 2 m high and the upper cover is 1.2 m high. The full-height and semiclosed shield door is 2.5 m high, whose shield door glass is 2 m high and the upper cover is 0.5 m high. The height of the semihigh shield door is 1.5 m, whose lower skirting board is 0.5 m high and the shielding door glass is 1 m high. The dimensions of the three doors in the longitudinal direction of 20 m (corresponding to one compartment) are shown in Figures 46, respectively.

3.3. Parameters of Explosion Source

Considering human bombs or luggage bombs as the two main sources for terrorist attack, the explosion height of all the explosion sources is determined as 1 m in this paper. According to the data of the explosive equivalent of various terrorist attacks given in No. 452 document issued by US Federal Emergency Management Agency [19], the maximum amount of belt and vest type human bomb is about 5 kg and 10 kg, respectively, and the average amount and maximum amount of ordinary bag type bombs are about 20 kg and 40 kg, respectively. Therefore, the four TNT equivalents (5 kg, 10 kg, 20 kg, and 40 kg) are selected in this paper.

Considering the actual cases of the explosion in the subway station and the simulation calculation effect, the center point of the platform is selected as location of the explosion source. And because the personnel at the edge of the platform and at the entrance to the bottom of the stairs are dense, these two positions are also selected. The specific locations of the explosion source are shown in Figure 7.

3.4. Calculation Conditions

Combined with the parameters of the explosion source in the subway platform and the characteristics of the subway platform structure, the calculation conditions of the explosion in the subway platform are shown in Table 2.

3.5. Numerical Model

The AUTODYN software is used to establish model for numerical simulation analysis. Also, in order to improve computational efficiency, the segmentation meshing is used, that is, within 0 m–7 m, 50 mm grid is used; for more than 7 m, 100 mm grid is used. With the consideration of the symmetry, 1/4 of the platform is selected to conduct modeling at the center of platform, and 1/2 of the platform is selected for explosion at the edge of the platform and the entrance to the bottom of the stair.

The platform board, wall, and column can be simplified as rigid materials. The boundary of the structure is set as the reflection boundary, and the transmission boundary is set at the exit of stair and the entrance of tunnel. The calculating result of shock wave propagation for the TNT exploding in one-dimension within 1 meter is remapped into the three-dimensional model by using REMAP technique. The numerical model is shown in Figure 8.

The arrangement of the monitoring points on the platform plane is as shown in Figure 9 (taking the 1/4 of the platform as an example). The monitoring points are set at intervals of 1 m. The vertical and horizontal intersections in the figure are the monitoring points on the platform plane, and the heights of the monitoring points are 1 m above the platform. At the stair, the monitoring points are set along the stairs slope where the person stands, and the heights of the monitoring points are 1.53 m, 2.07 m, 2.60 m, 3.13 m, 3.67 m, 4.20 m, 4.73 m, 4.75 m, and 4.79 m from 31 m to 39 m.

4. Analysis of Calculation Results

4.1. Shock Wave Load Distribution for Explosion of Subway Stations

According to equations (1) and (2), the values of overpressure and positive phase impulse can be obtained by extracting the pressure value at the position of the monitoring point for each numerical simulation, as shown in Figure 10.

The overpressure peak and positive phase impulse increase with the increase of TNT equivalent. The strength of the shock wave for the explosion on the platform under the different conditions is compared in Table 3.

When TNT equivalent is 40 kg, the maximum overpressure peak is the same because the overpressure at the blasting center is not affected by the surrounding environment. The shock wave continues to spread outwards, and the maximum positive phase impulses under different conditions vary due to the influence of the subway structure. When 40 kg TNT explodes at the stair entrance, the maximum positive phase impulse is the largest due to the “notch effect” of the stair exit, while the maximum positive phase impulse of the same equivalent exploding at the center of the platform under the other four working conditions (with shield door and unshield door) is relatively close, indicating that the shock waves are mainly reflected by the structural wall of the platform, and the shock wave positive phase duration is similar. The maximum positive phase impulse of the explosion at the edge of the platform is the smallest. That is because the shock wave expands to the train track of platform at the initial stage of propagation, and the shock wave energy is weakened by the effect of pressure relief of the train track.

When the explosive of different TNT equivalent explodes at the center of the platform without the shield door, the shock wave load on the stair slope is enhanced compared with that at other positions of the same section due to the “notch effect” of the stair. Also, the smaller the explosive charge is, the more obvious the “gap effect” is. Since the direction of the shock wave is perpendicular to the normal direction of the stair, the shock wave of large energy will continue to propagate forward due to inertia. The larger the explosive charge is, the greater the inertia is, and the less the pressure relief effect of the stair is. When explosive explodes at the center of the platform with the shield door, the shock wave overpressure and impulse increase near the shield door at the initial stage of explosion due to the reflection of the shield door.

The shock wave overpressure and impulse at the column and the stair will oscillate back and forth due to reflection and diffraction. When the shock wave propagates to the end wall of the platform, the shock wave overpressure and impulse near the end wall of the platform increase due to the reflection effect. The minimum values of the shock wave overpressure and impulse under different working conditions are all located at a certain position behind the stair due to the blocking effect of the stair on the shock wave.

4.2. Assessment for Personnel Casualty Area on the Platform

Using the comprehensive assessment method of personnel casualty obtained above, the casualties on the platform are divided into different levels, and different casualties are represented by different colors. Distribution maps of casualties under different conditions are obtained, as shown in Figure 11.

The distribution of casualties is irregular due to the influence of the column, stair, and shield door on the platform, but it presents radial expansion with the explosion source.

As the explosive charge increases, the safety area on the platform gradually decreases. When the explosive charge reaches 40 kg, in addition to the explosion that occurs at the stair, the safety area under other conditions is reduced to a certain extent behind the stair. That is because the stair and the column play a certain blocking effect on the shock wave. Moreover, when there is a shield door structure on the platform, the safety area on the platform is larger than other cases. The higher the shield door structure is, the larger the safety area behind the stair is. At the early stage of the explosion, the shock wave with large energy is blown through the shield door into the track area, resulting in a pressure relief effect. With the increase of the shock wave propagation distance, shield door on the platform plays a certain interception effect on the shock wave from the track area, which weakens the energy of the shock wave on the platform.

Due to the pressure relief effect of the stair, the shock wave overpressure and impulse are greatly increased, thereby increasing the degree of injury of the person at the position. When the explosive charge at the center of the platform is 5 kg and 10 kg, respectively, a separate light injury area appears at the slope of the stair, and other locations in the same section are safe. Similarly, when the explosives explode at the entrance to stair, the casualty at the slope of the stairs is also significantly higher than other areas of the same section.

Due to the reflection and diffraction of the shock wave caused by the platform column, the casualty of the explosion-facing surface of the column is greater than that of other locations in the same explosion distance and the casualty of the back-explosion surface of the column is relatively small compared with the surrounding. When the shock wave propagates to the end wall of the platform, the shock wave overpressure and impulse near the end wall increase due to the reflection of the end wall of the platform. So, when the explosives of 10 kg and 20 kg explode at center of the platform, there is a separate light injury area near the end wall.

The death and injury distance of personnel on the platform under different conditions (centered on the explosion source) is shown in Table 4.

According to Table 4, combined with dimensional analysis [20], the quantitative relationship between the explosive charge and the casualty distance for the condition of explosion at center of the platform with no shield door can be fitted, that is,where R is the casualty distance, m; a is a constant; and is the TNT equivalent, kg.

For the death distance, a is 1.8357, and the fitting correlation coefficient is 0.99798. For the distance of serious injury, a is 3.3693, and the fitting correlation coefficient is 0.99647.

5. Conclusion

According to the characteristics of shock wave propagation in the closed structure, ASII-based comprehensive casualty assessment method and process are established. Based on the investigation of the structure size of the subway platform and the cases of the explosion terrorist attack in the subway, the numerical models of the explosion in the subway platform under different conditions are established. The spatial distribution law of the overpressure and positive phase impulse of the shock wave in the platform is obtained. The ASII-based comprehensive assessment method for casualties is used to obtain the distribution map of casualties under different conditions. It deserves further study on how to develop an efficient and accurate calculation method to simulate the propagation of the explosion shock wave in the large-scale structure of the subway.

Data Availability

All data included in this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Key Research and Development Program of China (no. 2017YFC0803308-05) and National Social Science Foundation of China (grant no. 15BFX059). Special thanks are given to Ph.D Nafeng Wang for her helpful consultation in accomplishing this paper.