Mathematical Problems in Engineering

Volume 2013, Article ID 364791, 7 pages

http://dx.doi.org/10.1155/2013/364791

## Dynamics-Based Stranded-Crowd Model for Evacuation in Building Bottlenecks

^{1}School of Architecture, Harbin Institute of Technology, Harbin, Heilongjiang 150001, China^{2}College of Civil and Architecture Engineering, Northeast Petroleum University, Daqing, Heilongjiang 163318, China

Received 7 December 2012; Revised 5 January 2013; Accepted 15 January 2013

Academic Editor: Sheng-Yong Chen

Copyright © 2013 Lidi Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In high-density public buildings, it is difficult to evacuate. So in this paper, we propose a novel quantitative evacuation model to insure people’s safety and reduce the risk of crowding. We analyze the mechanism of arch-like clogging phenomena during evacuation and the influencing factors in emergency situations at bottleneck passages; then we design a model based on crowd dynamics and apply the model to a stadium example. The example is used to compare evacuation results of crowd density with different egress widths in stranded zones. The results show this model proposed can guide the safe and dangerous egress widths in performance design and can help evacuation routes to be selected and optimized.

#### 1. Introduction

One of the most important essential requirements for safety design of a building is safe evacuation of people from buildings in emergency conditions. From the early 20th century, there have been a lot of disasters that occurred during emergency evacuation of crowds. For example, 173 people died in the 1943 stampede incident in the subway station in London, UK; in the 1989 Hillsborough English FA Cup Stampede, 96 people died, and 400 were injured [1]; 20 soccer fans were killed during a stampede at a soccer game which became Africa’s worst sports-related disaster in Accra, Ghana, in 2001; 21 children died, and 47 people were injured in a grandstand failure at a school in Xuzhou, China, in 2002, and so forth [2]. So it is very significant to study on safe design of large scale buildings and influence factors during the evacuation process, since these incidents often occurred in high-density crowds buildings.

There have been more than one hundred years for the research on human behavior during the evacuation process. With improvement of computer technology, the research has developed to computer simulation [3, 4] from the initial postdisaster investigations. Egress modeling is one of the most important means of egress investigation [5, 6]. Popular models, such as the social force model [7] and the discrete model [8], the latter including lattice gas model and cellular automata model, are able to successfully simulate some typical phenomena observed in pedestrian dynamics. Software based on the above models, EVACENT4, STEPS, building EXODUS, EXITT, Legion, and so forth, have been applied widely to simulate and evaluate very complex factors influencing people evacuation. Recently, some experimental results have also been presented to validate the models of pedestrian flow and evacuation [9–11].

However, it is very complicated to analyze human behavior in emergency conditions. Furthermore, systematic study on panic behavior and quantitative theories of predicting crowd dynamics are few. Computation for especially determination of people evacuation is so complex that their application is limited. Thus, several authors have devoted their work to the research of simple model and the calculation of evacuation time. Zhang et al. [12, 13] proposed an SCM to calculate the number of stranded crowd with different egress widths and modify the evacuation time; Liu et al. [14] developed an evacuation model of crowd streams at a stadium to calculate and analyze the total evacuation time.

Relieving panic and stress of people from building space and construction environment is the duty of architects. However, evacuation design research has not been properly established in architectural discipline. Thus, the current study hopes to serve such a need and is inspired by the study of Zhang et al. [12] and Liu et al. [14]. In this paper, a simple continuous-time model based on the dependency of physical parameters (density, velocity, and flow-rate) is set up to discuss the concepts of safe space design in terms of evacuation processes efficiency in high-density crowd situations and panicked behavior. This paper intends to introduce this model and apply the analysis results to architectures, to provide modifications in architectural layouts to optimize evacuation.

A brief description of the organization of this paper is as follows. The main calculation principles in this paper are stated in Section 3. The parameters analysis of Section 3 is described in Section 2. The application is given in Section 4. The results and analysis of the application are considered in Section 5. Finally, a conclusion is drawn to summarize the paper.

#### 2. Preliminaries

##### 2.1. Concept

The evacuation route of pedestrians in a stadium can be simply attributed as leaving from spectators’ stands to vertical passages, from passages to doors or stairs, and then to safe areas. Traffic nodes among the evacuation routes, such as changed building space, narrowed channels, incompletely opened exits, merging points of several crowd streams (Figure 1(a)), and stairs, are called “Building Bottlenecks,” and where potential hazards would break out in the case of high gathering crowd. In an urgent situation, crowd streams from spacious areas rushing here would hinder the positive streams flowing out, increase crowd density, and then form an irregular succession of arch-like blockings at exit bottlenecks (Figure 1(b)), which instantaneously causes a cease in evacuation flow, and when the arches break, pedestrians will suddenly leave just like avalanche bunches [15, 16]. More seriously, this phenomenon can cause people psychological panic and induce collective human behavior, a crowd stampede, which is one of the most disastrous forms, often leading to fatalities as people are crushed or trampled [17].

##### 2.2. Parameters Analysis

During the last few decades, pedestrians flow and evacuation have attracted researchers’ attention. They concluded that egress evacuation capacity is related to egress width, , m, egress stream flow-rate, , persons/(ms), crowd velocity, , m/s, and crowd density, , persons/m^{2}.

Egress stream flow-rate is the number of persons passing through per unit of distance width of evacuation passage or exit per unit of time, person/(ms), crowd velocity is moving distances per unit of time, m/s, and crowd density is the number of persons per unit area, persons/m^{2}.

It is essential to understand the main characteristics of their relationship among the above parameters. In the past, researches about these parameters were mainly through two kinds of means: field observation and video record. So far, Predteckenskii and Milinski [18], Fruin [19], Ando et al. [20], Paul [21], Nelson and Maclennan [22], Smith and Petly [23], and Lu [24] et al., have already accumulated a large amount of observation data, from which typical relationship curves between the crowd density and velocity was summarized by Thompson (Figure 2) [25]. In general, in different environments and with different observation methods, each result of fleeing velocity varies slightly because of all kinds of constraints, but the changing trend of these curves is similar; that is, an increase in crowd streams density makes the distance between people decrease, and makes personnel movement velocity slow; on the contrary, a decrease in the density makes personnel movement velocity increase.

Based on influence mechanism of personnel acceleration that is produced by people around in urgency evacuation from buildings, Lu et al. [24] proposed a velocity dynamics equation under crowding environment, which is a function of density: where is fleeing velocity, m/s, generally taken as 1.669 m/s, the range of is 0.25~0.44, is 0.014~0.088, and is 0.15~0.26.

Egress stream flow-rate depends on velocity and density, showed in (2): Then the number of people evacuating through exit or passage is expressed as where is the width of the evacuation exit or passage, m and is evacuation time, s.

Figure 3 shows the essential relationship of three parameters [24]. is free flow velocity, m/s, is the maximum flow-rate value, person/(ms), and and are critical density and velocity, respectively, persons/m^{2}, m/s. When crowd density is lower than the critical density, flow-rate increases as density increases, pedestrian flow in free flow state, but when crowd density is greater than the critical density, flow-rate decreases as density increases, pedestrian flow in blocked state, where flow-rate and velocity decrease and density increases. is the safe density; in this situation, crowd is gathering but slowly flowing, and is the maximum density; flow-rate and velocity all equal zero, in a completely jammed situation [26].

The Green Guide [26] defines the limit of 4 persons/m^{2} as a safe density for gathering crowd, that is, persons/m^{2}, and proposes persons/m^{2}. Under extreme conditions, density of 15 persons/m^{2} has been observed. According to the Chinese physical characteristic, we choose persons/m^{2} as a safe density for a moving crowd and persons/m^{2} as a maximum density [27] in the jammed situation.

#### 3. Calculation Principle

Streams of pedestrians evacuating towards a certain direction and the whole evacuation path inevitably encounter changed building space, such as path nodes, where streams merged into one, narrowed egresses, stairs, doors, steps, and so forth; these bottlenecks possibly become stranded zones due to clogging. Figure 4 shows a sketch of crowd flow. Assuming all leaving crowd will flow through the bottleneck into the next space, we define people who continue to flow in the zone as inflow crowd and out of the zone as outflow crowd; if jamming happens, crowd staying the zone as stranded crowd. Obviously, the balance principle of import and export applied and the number of stranded crowd should be the difference between numbers of inflow and outflow.

##### 3.1. The Number of Inflow Crowd

Assuming several pedestrians’ streams from internal passages merge into the bottleneck and then walk towards the next egress, the number of inflow crowd can be
where is the number of inflow crowd from branch entrances; and are outflow velocity and density of branch entrances No. at time , m/s and persons/m^{2}; is the width of entrances No. , m; is the number of branch entrances; is the evacuation time, s.

##### 3.2. The Number of Outflow Crowd

During the time , the evacuation crowd number of flowout of the bottleneck is
where is the number of outflow crowd; is the time that the first person reaches the bottleneck exit, s; is the outflow velocity of the exit at time , m/s; is the width of total exit, m; is the density of bottleneck zone in front of the exit at time , persons/m^{2}.

##### 3.3. The Number of Stranded Crowd

Till the evacuation ends, the number of stranded crowds in the bottleneck zone is

##### 3.4. The Total Evacuation Time

When crowd streams all flow out of the exit to the next space, the evacuation ends. The time is where is the number of total evacuation people; , are the outflow velocity and flow-rate of the exit at free flow, m/s, person/(ms).

The model is supposed that, when the evacuation velocity of the exit is equal to zero, no one is moving, and the exit is absolutely blocked, crowd density in the bottleneck is very high and arch-like clogging emerges, which must be prevented. This model also can calculate the smallest safe egress width based on the safe density and can calculate total evacuation time.

#### 4. Application Example

In some stands of Tianjin Olympic Center Stadium as an example, evacuation routes include three longitudinal passages, one transverse passage, the hall in front of the egress, and the egress. The hall in front of the egress is a typical building bottleneck and becomes our main study field. The passage exit structure is shown in Figure 5.

##### 4.1. Parameters Selection

Assuming the discussed stands have a seating capacity of 1400; that is,; all persons are evenly distributed in the stands. Three longitudinal passages have the same width , m, enough to allow up to two spectators to walk simultaneously. Transverse passage width is m, stranded-crowd zone width is m and egress width , is a varying value.

The distance from one longitudinal passage entrance to the stranded-crowd zone is , m. The length from the egress to the edge of transverse passage is , m. The area of stranded-crowd zone is , m^{2}.

According to the Design Code for Sports Building (JCJ31-2003) of China [28], assuming one stream of crowd flow is 40 persons/min and one crowd stream breadth is 0.55 m, then flow-rate is

In this paper, we ignore the time of spectators leaving seats to longitudinal passages and only consider the process of spectators beginning to flow into the transverse passage from stands entrances.

##### 4.2. The Number of Stranded Crowd

At the very beginning of the evacuation, because the number of people flowing into the bottleneck is small, they can freely flow out of the egress. As people numbers increase, stranded crowd appears.

The number of stranded crowd can be calculated as follows:

() where is the initial number of stranded crowd in the bottleneck during the time ; is the desired velocity during the evacuation, m/s; is the nearest distance from stands entrances to the bottleneck egress, ; is inflow rate, person/(ms) and ;

()

Because this model is based on crowd dynamics, the relationship between crowd density and crowd velocity of the egress of the bottleneck as shown in (1) is used, which is more suitable for the high-density crowd of buildings. We adopt the parameters from HMSO [24]; that is, , , and , and substituting (1) into (10), then

At time , the egress crowd velocity is equal to zero, which means no one can pass through the egress. is the number of outflow crowd during the time ; is the time step, because and are variables of time, and is the dynamic value with time, and , are the number of stranded and outflow crowd at the prior time, the time .

#### 5. Results and Analysis

##### 5.1. The Curves of the Number of Stranded Crowd and Outflow Crowd

From the above calculation, the curves of the number of stranded crowd and outflow crowd varying with time are obtained and compared with four different egress widths ( m, m, m, and m).

() m: Figure 6 shows that the changes of two curves are not linear, and their gradients are similar. In the initial time, stranded-crowd number in the bottleneck is greater than outflow number, but when the time s, they are equal; then outflow number is greater than stranded-crowd number. When (the velocity of outflow) approaches zero, stranded-crowd number increases quickly, and the density of the zone increases up to (p/m^{2}); at that time, , s, and the egress is jammed.

() m: Figure 7 shows that the gradients of two curves are apparently different. Outflow crowd number soon exceeds the stranded and becomes linear, but the number of stranded-crowd has been slow in growth until the time s, and the rate of stranded-crowd number growth begins to have a sensible rise. When the time s, , and (p/m^{2}), the egress is jammed.

() m: Figure 8 shows the curve of outflow crowd number is simply linear and increasing rapidly, and the curve of the stranded number is linear initially. However, when the time s, the curve is horizontal, which means the stranded number is constant and the density of bottleneck (p/m^{2}). In that situation, we consider the crowd is safe, and evacuation crowd could freely flow out without any obstacle; m/s.

() m: The egress width in a stadium should be greater than or equal to the sum of the widths of served stand passages or egresses [28]; so we choose m for calculation and that is equivalent to the egress completely opened. Figure 9 shows the calculation results are similar to the curves when the width m, except for when the density decreases, (p/m^{2}), m/s. All evacuation is not obstructed at all. In addition, the results are according to the regulation of the Design Codes for Sports Building (JCJ31-2003) of China [28], which further verifies this model’s effectiveness.

##### 5.2. The Relationship of Crowd Density and Egress Width in the Bottleneck

Comparing the density with the egress width from 1.8 m to 2.2 m in the bottleneck, we observed that it slips very fast (Table 1); so we consider there may be a critical width which is the threshold value of safe width. In order to solve this problem, the density with different egress widths is calculated, and the results are shown in Figure 10.

From the above analysis, , , the corresponding safe width , and dangerous width are labeled, respectively, in Figure 10.

When and m, the density of the crowd is extremely high, and the physical movement of pedestrians is almost impossible. This situation is terribly dangerous, because while pedestrians are still standing, when crushing occurs, the high pressure within the crowd compresses the lungs of some individuals in the crowd, producing fatalities from asphyxiation [17]. So the width of egress must be bigger than , and inspection of the bottleneck density is very significant.

When , m; although the width increases a little, the density decreases significantly, which seemingly small changes to the architecture of the surroundings can have such a large impact on the behavior of crowd, particularly when a crowd is panicked. After the density suddenly inclines up to persons/m^{2}, with the width increasing, density is smoothly decreasing; so the safe width, , is the threshold width value, which does not lead to high density in theory.

#### 6. Conclusion

The people evacuation calculation model, based on the crowd dynamics when the exit jammed, is presented in this paper. The main characteristics of the evacuation streams, people’s safety under emergency conditions, and comparison of the results using dynamics formula with different egress widths for the movement in a stadium example are all discussed.

Safe egress width is within a stone’s throw of dangerous egress width, when a crowd is panicked. Small changes to the architecture of the surroundings can have large effects on the behavior of crowd.

This stranded-crowd model can calculate total evacuation time, which will be carried out in future work.

#### References

- Y. Liu, D. M. Liu, B. Norman et al., “Analysis of evacuation performance of merging points in stadiums based on crowd simulation,” in
*Proceedings of the 12th Conference of the In-ternational Building Performance Simulation Association (IBPSA, BS'11)*, pp. 2651–2658, November 2011. - X. Pan, C. S. Han, K. Dauber, and K. H. Law, “Human and social behavior in computational modeling and analysis of egress,”
*Automation in Construction*, vol. 15, no. 4, pp. 448–461, 2006. View at Publisher · View at Google Scholar · View at Scopus - S. Chen, Y. Zheng, C. Cattani, and W. Wang, “Modeling of biological intelligence for SCM system optimization,”
*Computational and Mathematical Methods in Medicine*, vol. 2012, Article ID 769702, 10 pages, 2012. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Wen, W. Zheng, J. Zhu, X. Li, and S. Chen, “Elman fuzzy adaptive control for obstacle avoidance of mobile robots using hybrid force/position incorporation,”
*IEEE Transactions on Systems, Man and Cybernetics C*, vol. 42, no. 4, pp. 603–608, 2012. View at Publisher · View at Google Scholar · View at Scopus - K. J. Zhu and L. Z. Yang, “The effects of exit position and internal layout of classroom on evacuation efficiency,”
*Acta Physica Sinica*, vol. 59, no. 11, pp. 7701–7707, 2010. View at Google Scholar · View at Scopus - W. G. Song, Y. F. Yu, and T. Chen, “Influences of exit condition on pedestrian evacuation,”
*Fire Safety Science*, vol. 12, no. 2, pp. 100–105, 2003. View at Google Scholar - D. Helbing and P. Molnar, “Social force model for pedestrian dynamics,”
*Physical Review E*, vol. 51, no. 5, pp. 4282–4288, 1995. View at Publisher · View at Google Scholar - M. Isobe, T. Adachi, and T. Nagatani, “Experiment and simulation of pedestrian counter flow,”
*Journal of Fire Protection Engineering*, vol. 4, no. 3, pp. 81–92, 1992. View at Publisher · View at Google Scholar - V. Papinigis, E. Geda, and K. Lukošius, “Design of people evacuation from rooms and buildings,”
*Journal of Civil Engineering and Management*, vol. 16, no. 1, pp. 131–139, 2010. View at Publisher · View at Google Scholar · View at Scopus - S. Chen, W. Huang, C. Cattani, and G. Altieri, “Traffic dynamics on complex networks: a survey,”
*Mathematical Problems in Engineering*, vol. 2012, Article ID 732698, 23 pages, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - C. Cattani, R. Badea, S. Chen, and M. Crisan, “Biomedical signal processing and modeling complexity of living systems,”
*Computational and Mathematical Methods in Medicine*, vol. 2012, Article ID 298634, 2 pages, 2012. View at Publisher · View at Google Scholar - Q. S. Zhang, M. Liu, C. H. Wu, and G. Zhao, “A stranded-crowd model (SCM) for performance-based design of stadium egress,”
*Building and Environment*, vol. 42, no. 7, pp. 2630–2636, 2007. View at Publisher · View at Google Scholar · View at Scopus - Q. S. Zhang, M. Liu, J. Liu, and G. Zhao, “Modification of evacuation time computational model for stadium crowd risk analysis,”
*Process Safety and Environmental Protection*, vol. 85, no. 6, pp. 541–548, 2007. View at Publisher · View at Google Scholar · View at Scopus - Q. Liu, H. Yang, H. P. Lu, and S. Jing, “Study on the dispersal of people flow in the stadium,”
*China Civil Engineering Journal*, vol. 37, no. 10, pp. 92–98, 2004. View at Google Scholar - S. Y. Chen, Y. H. Wang, and C. Cattani, “Key issues in modeling of complex 3D structures from video sequences,”
*Mathematical Problems in Engineering*, vol. 2012, Article ID 856523, 17 pages, 2012. View at Publisher · View at Google Scholar - D. Helbing, I. Farkas, and T. Vicsek, “Simulating dynamical features of escape panic,”
*Nature*, vol. 407, no. 6803, pp. 487–490, 2000. View at Publisher · View at Google Scholar · View at Scopus - R. S. C. Lee and R. L. Hughes, “Exploring trampling and crushing in a crowd,”
*Journal of Transportation Engineering*, vol. 131, no. 8, pp. 575–582, 2005. View at Publisher · View at Google Scholar · View at Scopus - V. M. Predtechenskii and A. I. Milinski,
*Planning for Foot Traffic Flow in Buildings*, Stroiizdat, Moscow, Russia, 1969. - J. J. Fruin,
*Pedestrian Planning and Design*, Metropolitan Association of Urban Designers and Environmental Planners, New York, NY, USA, 1971. - K. Ando, H. Ota, and T. Oki, “Forecasting the flow of people,”
*Railway Research Review*, vol. 8, no. 45, pp. 8–14, 1988. View at Google Scholar - J. L. Paul, “Effective width model for evacuation flow in buildings,”
*Society of Fire Protection Engineers*, vol. 4, no. 6, pp. 215–232, 1980. View at Google Scholar - H. E. Nelson and H. A. Maclennan,
*Emergency Movement*, NFPA, Quincy, Mass, USA, 1995. - K. Smith and D. N. Petley,
*Environmental Hazards Assessing Risk and Reducing Disaster*, Routledge, London, UK, 5th edition, 2008. - J. A. Lu, Z. Fang, Z. M. Lu et al., “Mathematical model of evacuation speed for personnel in buildings,”
*Engineering Journal of Wuhan University*, vol. 35, no. 2, pp. 66–70, 2002. View at Google Scholar - P. A. Thompson and E. W. Marchant, “Computer and fluid modelling of evacuation,”
*Safety Science*, vol. 18, no. 4, pp. 277–289, 1995. View at Publisher · View at Google Scholar · View at Scopus *Guide to Safety at Sports Grounds (Green Guide)*, Football Licensing Authority Department for Culture, Media and Sport (DCMS), TSO (The Stationery Office), 5th edition, 2008.- J. Chen, P. J. Zhang, W. Yang, and Q. Li, “Personnel evacuation strategy for theater based on systematic dynamics model,”
*Journal of Natural Disasters*, vol. 14, no. 6, pp. 125–132, 2005. View at Google Scholar · View at Scopus *Design Code for Sports Building (JCJ31-2003) of China*, The Industry Standard of the People's Republic of China, Beijing, China, 2003.