- Department of Engineering ‘Enzo Ferrari’ (DIEF), University of Modena and Reggio Emilia, Modena, Italy
A reliable prediction of the human-induced vibrations of footbridges relies on an accurate representation of the pedestrian excitation for different loading scenario. Particularly, the modeling of crowd-induced dynamic loading is a critical issue for the serviceability assessment of footbridges. At the design stage, the modeling of crowd loading is often derived from single pedestrian models, neglecting the effect of the structural vibrations as well as the interactions among pedestrians. A detailed description of the crowd behavior can be achieved employing a social force model that describes the different influences affecting individual pedestrian motion. These models are widely adopted to describe the crowd behavior especially in the field of evacuation of public buildings, public safety and transport station management while applications in the serviceability assessment of footbridges are less common. To simulate unidirectional pedestrian flows on footbridges, this paper proposes a parameter calibration of the Helbing’s social force model performed adopting the response surface methodology. Parameters of the social force model are calibrated so as to represent the fundamental relation between mean walking speed and density of the pedestrian crowd. The crowd-induced vibrations are then simulated by modeling each pedestrian in the crowd as a vertical load that crosses the footbridge with time varying trajectory and velocity estimated from the calibrated social force model. Finally, results are compared to those obtained from a multiplication factor approach proposed in literature. This considers the crowd as a uniform distribution of pedestrians with constant speed and given synchronization level and the footbridge response is evaluated as the response to a single pedestrian scaled by a proper enhancement factor.
1 Introduction
The serviceability assessment under pedestrian-induced vibrations is a key aspect in the design of modern footbridges, which are usually slender and light structures. This can make them particularly vulnerable to human-induced vibrations. The great interest of the scientific community in the serviceability assessment of modern footbridges has been mainly motivated by swaying problems experienced by several long-span footbridges such as the Millenium Bridge in London (Dallard et al., 2001) and the Solferino footbridge in Paris (Danbon and Grillaud, 2005). However, this is not a new phenomenon and also earlier researches can be found in the literature (Matsumoto et al., 1978; Fujino et al., 1993). In the last 20 years, much efforts were devoted to the characterization of the load induced by a single pedestrian walking or running and the subsequent extrapolation of these findings to the case of a moving crowd (Racic et al., 2009; Caprani et al., 2012; Pan et al., 2017; Younis et al., 2017; Bassoli et al., 2018a; Wang et al., 2019). Some of these studies have contributed to the publication of several guidelines and standards for the design of slender footbridges under pedestrian action, such as Butz et al. (2008). However, codes of practice still present some shortages, especially when dealing with pedestrian crowds. In those cases, in fact, the simulated dynamic response often differs from the experimental evidence (Živanović et al., 2010; da Silva et al., 2013). This is mainly due to the simplified crowd models adopted, based on an enhancement factor related to the pedestrian number that multiplies the single person load.
Recently, new important aspects of the pedestrian behavior have been introduced and taken into account to improve the modeling of pedestrian crowds, which are the crowd behavior, the human-structure interaction and the inter- and intra-subject variability of the walking force (Venuti et al., 2016; Jiménez-Alonso and Sáez, 2017; Bassoli et al., 2018b; Moreu et al., 2020). The latter is usually simulated adopting a probabilistic approach, namely modeling the walking parameters as random variables (Živanović et al., 2007; Van Nimmen et al., 2014, 2020). Several models able to address the three above mentioned key aspects have been proposed, among which da Silva et al. (2013); Venuti et al. (2016); Jiménez-Alonso and Sáez (2017). They share a common pattern based on the coupling of two sub-models: a pedestrian-structure interaction model and a crowd model. The most common approach, but not the only one, to model the pedestrian-structure interaction is to represent each pedestrian as a single degree of freedom (SDOF) system crossing the structure in addition to a external force attached to its base. The external force represents the walking force caused on a rigid surface while the SDOF system contributes to modify the dynamic properties of the crowd-structure coupled system (Venuti et al., 2016). As far as the crowd behavior is concerned, two modeling approaches exist: macroscopic and microscopic models. The firsts are based on the modeling of the crowd behavior as a continuous flow of a fluid. The second models allow for a more detailed description of the crowd thanks to the evaluation of the time varying position and velocity of each pedestrian. Macroscopic models consider the behavior of the crowd as a whole and assume the continuity of flow, assumption that may not be satisfied for low and medium pedestrian densities (Carroll et al., 2012). Moreover, the variability in the walking parameters among the crowd cannot be accounted for. The above mentioned limitations can be overcome adopting the microscopic modeling, in which the behavior of each pedestrian is governed by the different motivations and influences that he/she experiences according to the equations of particle dynamics (Jiménez-Alonso et al., 2016). This approach was initially proposed by Helbing and Molnar (1995), who described the various stimuli experienced by each pedestrian in the crowd as social forces. This allows to consider that i) the pedestrian crowd may not be uniformly distributed, ii) each pedestrian can contribute differently to the modal excitation depending on his own pacing frequency and iii) pedestrians may react by stopping or slowing down depending on the perceived level of footbridge vibrations (Carroll et al., 2012).
Social force models (SFM) have been widely adopted to simulate the behavior of pedestrian crowds in the fields of transport station management, safety assurance of large pedestrian flow events and building evacuation (Chen et al., 2018). Their large diffusion depends on the good capability of describing movement processes using simple mathematical formulations. Starting from the original version proposed by Helbing and Molnar (1995), different developments have been proposed to adapt the model to the application field. A comprehensive review summarizing the existing social force models for pedestrian traffic is presented in Chen et al. (2018). Furthermore, some applications of social force models for the serviceability assessment of footbridges under pedestrian flows can be also found in literature, both with reference to lateral (Carroll et al., 2012) and vertical (Jiménez-Alonso et al., 2016; Venuti et al., 2016) loads. The crowd model adopted by Carroll et al. (2012) is adapted from the one proposed by Langston et al. (2006) and based on the findings of Helbing et al. (2000). In particular, Langston et al. (2006) studied a multi-circle pedestrian model to simulate a single enclosure entry scenario. However, the calibration of the parameter describing pedestrian rotational movements remains a challenging task. Moreover, the parameters adopted in Helbing et al. (2000) are chosen to simulate the mechanisms of panic and jamming by uncoordinated motion in crowds. Similarly, social force model parameters adopted in Jiménez-Alonso et al. (2016) are based on Carroll et al. (2012) and Helbing and Molnar (1995). In the latter case, simulations are performed to describe the self-organization of collective phenomena of pedestrian behavior considering two pedestrian groups trying to pass a narrow door and walking in opposite directions. Finally, Venuti et al. (2016) adopt a crowd model originally proposed by Cristiani et al. (2011) who described the interaction mechanisms by velocity terms.
This paper proposes a parameter calibration of the Helbing’s social force model (Helbing et al., 2005) for the vibration serviceability of footbridges. For a realistic simulation of the pedestrian traffic, model parameters are adjusted so that the simulated pedestrian flow fits the fundamental speed-density relation proposed by Weidmann (1993). The optimization process is carried out adopting the response surface methodology (RSM), based on the approximation of the objective function using simple analytical functions. The RSM is widely adopted in model optimization problems thanks to an adequate accuracy of results combined with computational efficiency (Vincenzi and Savoia, 2015; Myers et al., 2016). Once the social force model is properly calibrated, simulations are performed for different pedestrian flows allowing for the evaluation of the time-varying position and velocity of each pedestrian in the crowd. These are adopted as input data for detailed simulations of the pedestrian-induced loads on a step-by-step basis. Finally, simulation results are compared to those obtained from a multiplication factor approach proposed by Caprani et al. (2012). Based on the assumption that the crowd is composed of uniformly distributed pedestrians with constant speed and given synchronization level, they calculate enhancement factors to evaluate the footbridge response to crowd loads scaling the single pedestrian load.
The paper is organized as follows. The basics of the social force model are presented in section 2. Section 3 describes how the simulated speed-density relation is obtained starting from the results of the social force model. The optimal model parameters are evaluated in section 4 and the calibrated social force model is adopted in section 5 to simulate crowd-induced vibrations. Results are compared to those obtained adopting a multiplication factor approach in section 6. Finally, conclusions are drawn in section 7.
2 The Social Force Model
The social force model describes the behavior of individual pedestrians by a superposition of social forces that reflect motivations and environmental influences. The motivation to move toward a specific destination with a certain desired velocity is represented by a driving term, while repulsive forces describe the tendency to keep a certain distance from other pedestrians, borders, obstacles, and dangers. Finally, attractive forces express the tendency of a group or family members to stay together or to move towards window displays, sights or unusual events. The basics of the social force model are shown in the following, while the reader is referred to Helbing and Molnar (1995); Helbing et al. (2000, 2005) for a detailed description.
The time-varying velocity of a generic pedestrian α,
where
In Eq. 2,
where NP is the total number of pedestrians in the simulation. The social force, pedestrian positions and velocities are all function of time t. However, this dependency may be omitted in the following equations for the sake of simplicity.
The driving term represents i) the intention of each pedestrian to walk with a desired speed
where
where
The interaction repulsive forces describe that the pedestrian α tends to keep a situation-dependent distance from the other pedestrians β:
The first term of the right-hand side of Eq. 6 describes the territorial effect, namely the tendency to respect the private sphere of each pedestrian, and helps avoiding collisions if there are sudden velocity changes. The second term express physical interactions at high pedestrian densities and pushy crowds, when frictional effects are ignored. Note that other definitions of the repulsive force can be found in literature, such as the velocity dependent interaction force described in Johansson et al. (2008) and Helbing and Johansson (2009). However, Johansson et al. (2008) pointed out that the prediction of pedestrian motion does not significantly improve by including additional speed-dependent parameters, in light of an increasing computational time and complexity. The parameters Aα and Bα identify the repulsive interaction strength and range, respectively, and partially depend on cultural influences. Although the dependence on α would enable defining those parameters for each individual, it is usually assumed
FIGURE 1. Action area of the territorial effects for a pedestrian (black dot) moving to the right considering different
Furthermore, pedestrians tend to keep a certain distance from borders, such as guardrails or walls. The repulsive effect of borders is similar to the repulsion among pedestrians except for the anisotropic behavior. It is described by a force
Similarly to the repulsive interaction among pedestrians,
Pedestrian behavior can also be affected by the so-called attractive effects, which represent the tendency of pedestrian groups, such as family members or tourist groups, to walk together. These attractive forces can be modeled accordingly to the repulsive forces among pedestrians of Eq. 6 with longer interaction range and negative interaction strength.
The social force model of Eqs. (1–7) is described by
3 Evaluation of the Simulated SPEED-DENSITY Relation
This section shows how the speed-density relation can be evaluated from the results of the social force model. First, input parameters and assumptions made for the simulations of the pedestrian flow are listed in section 3.1, while an example of the simulated speed-density relation is presented in section 3.2.
3.1 Problem Formulation
The desired velocity of each pedestrian
As regards the repulsive force among pedestrians, the first term of the right-hand side of Eq. 6 is considered, which represents the tendency to respect the private pedestrian sphere and accounts for the anisotropic pedestrian behavior. It is assumed that the repulsive force between two pedestrians acts when
The model parameters to be calibrated are
Because the destination of a pedestrian crossing the footbridge is typically the opposite end and not a specific point, the desired destination is re-evaluated at each iteration. By doing this, in the event of interaction with other pedestrians or with the borders, the pedestrian is able to re-evaluate its desired destination as the point on the footbridge edge closest to its current position (Parry, 2007).
Finally, according to applications of the social force model found in literature, attractive effects
3.2 Simulated Speed-Density Relation
Simulations are performed considering a footbridge 10 m long and 3 m wide. The footbridge is completed by an access route 60 m long and a 300 m long way out (Figure 2). Access road and way out are necessary for the modeling of a roughly stationary pedestrian flow on the footbridge. Initial pedestrian positions are randomly distributed in the access route. Considering that the dimension of each pedestrian is 2rα, the uniform random distribution of initial positions is defined on the condition that pedestrians shall not overlap. An example of the random initial positions of a pedestrian group is presented in Figure 2. Finally, the way out needs to be long enough so as not to stop the pedestrian flow in the central part.
FIGURE 2. Geometry of the footbridge (gray) and of the access route and way out (white). Black dots: example of random initial positions of pedestrians. Dimensions are in meters.
To define the simulated speed-density relation, simulations are repeated for different-sized pedestrian groups. In particular, groups of 25 to 500 pedestrians are considered. As results depend on the random initial position and velocity of each pedestrians, for each group dimension the simulation is repeated 10 times to increase the number of results (namely NA = 10). The total analysis duration ranges from 125 to 195 s depending on the group dimension. The analysis duration is chosen so as to allow approximately all pedestrians to cross the footbridge.
The social force model allows for the evaluation of the positions and velocities (both in X and Y direction) of each pedestrian during the whole analysis. As the parameters of interest are limited to the footbridge (highlighted in gray in Figure 2), results related to positions in X direction between 60 and 70 m are extracted and analyzed in the rest of the paper. Figure 3A shows an example of the number of pedestrians present on the footbridge at each time instant for the case of 350 pedestrians. It can be observed that the pedestrian number gradually increases at the beginning of the analysis. Once the majority of pedestrians arrived on the footbridge, the pedestrian number remains roughly constant as long as they start leaving the footbridge. Density and velocity of people walking are evaluated according to two different criteria:
FIGURE 3. (A) Number of pedestrians on the footbridge and (B) instantaneous crowd density (black dots) and velocity (red dots) for the case of 350 pedestrians.
• Criterion 1: crowd density and mean pedestrian velocity are evaluated at each time instant. An example of the instantaneous velocity and density obtained for the case of 350 pedestrians is shown in Figure 3B. The analysis can be ideally divided in three time windows: 0 − T1, T1 − T2, and T2 − Tend. T1 and T2 represent the first and the last time instant when the 80% of NP,max pedestrians is on the footbridge, being NP,max the maximum number of pedestrians simultaneously present on the footbridge. As regards the example of Figure 3B, T1 = 15.7 s, T2 = 97.5 s and Tend = 125 s. In the first time window, the mean pedestrian velocity decreases as the pedestrian density increases, while in the central part of the analysis they remain roughly stationary. Finally, from T2 to Tend the pedestrian density decreases since pedestrians start leaving the footbridge. In this case, a decreasing pedestrian density does not correspond to an increasing mean velocity as pedestrians are not free to move but they are slowed down by pedestrians that have already left the footbridge. To avoid meaningless results, instantaneous densities and velocities are evaluated up to T2. Instantaneous values can be represented in a speed-density diagram, as shown in Figure 4A, where the relation between increasing density and decreasing mean velocity can be clearly observed.
FIGURE 4. (A) Example of the crowd density and velocity from criterion 1 (blue dots) and criterion 2 (green dot) for the case of 350 pedestrians (B) example of the crowd density and velocity from criterion 1 (blue dots) and criterion 2 (green dots) and simulated speed-density relation (red dots) for group dimensions from 25 to 500 pedestrians.
• Criterion 2: mean pedestrian density and velocity are evaluated over the central time window of the analysis (T1 − T2), namely when the pedestrian flow is roughly stationary. With the reference to the previous example, this results in the green dot of Figure 4A.
Figure 4B presents an example of the speed-density diagram obtained considering results from criterion 1 and criterion 2 for different analyses, namely 11 pedestrian group sizes and 10 analyses for each case. Finally, the simulated speed-density relation represented by red dots is obtained from criterion 1 by evaluating the mean crowd velocity in correspondence of each pedestrian density.
4 Model Parameter Calibration
The social force model described above is calibrated so that the simulated pedestrian flow fits the fundamental speed-density relation presented by Venuti and Bruno (2007) in the form:
where
The model is calibrated through the response surface methodology, a widely adopted mathematical and statistical method for modeling and analyzing processes in which the response of interest is affected by various variables. The aim of the RSM is to find the values of input variables that produce the best value of the response (Myers et al., 2016). This study means to find the values of the model parameters to be calibrated, i.e.,
4.1 Response Surface Methodology
The basic idea of the RSM is to approximate the objective function using simple and explicit interpolation functions. Originally proposed by Box and Winson (1951) to optimize chemical processes, the application of RSM has been subsequently extended to other fields such as engineering problems involving complex and time consuming analyses in order to reduce the computational effort (Vincenzi and Savoia, 2015).
According to the RSM, the objective function H can be approximated by an analytical estimation function
with p the D-dimensional vector collecting the model parameters to be calibrated and
where Q is a D × D matrix gathering the quadratic terms, L is a D-dimension vector collecting the linear coefficient and β0 a constant.
Khuri and Cornell (1996) proposed the procedure described in the following to analytically correlate values of identification parameters and
where βi are the unknown coefficients of the RS. Considering NS observations, namely NS evaluation of
where
In Eq. 13, all observations have the same weight. To increase the accuracy of the RS close to the optimal solution, Kaymaz and McMahon (2005) and Myers et al. (2016) proposed the weighted regression method where weights of sampling points
where w is a NS × NS diagonal matrix of weight coefficients. The weight coefficients wi can be evaluated as (Vincenzi and Savoia, 2015):
where
Once the parameters β of the response surface are evaluated, the optimal parameter vector p* minimizing
Finally, the estimation of the optimal solution is further improved by iteratively re-calibrating the response surface close to p*. A detailed flowchart of the updating process is presented in Figure 5.
4.2 Results and Discussion
Results of the social force model calibration are presented in Table 2 and Figure 6. Table 2 lists the ranges of variation and the optimal values of the updating parameters, while Figure 6 shows the speed-density relation obtained from the calibrated social force model together with the theoretical curve. In this case, the simulated speed-density relation is calculated considering six pedestrian group dimensions (from 250 to 500 pedestrians) and 10 analyses for each case. A pretty good agreement between the theoretical and simulated curve can be observed, even though slightly higher discrepancies are found for very low densities (ρ< 0.3 ped/m2).
FIGURE 6. Crowd density and velocity from criterion 1 (blue dots) and criterion 2 (green dots) and simulated speed-density relation (red dots) obtained from the calibrated model. Black curve: theoretical speed-density relation.
As an example,Figures 7A,B show, respectively, the effect of the model parameter
FIGURE 7. Crowd density and velocity from criterion 1 (blue dots) and criterion 2 (green dots) and simulated speed-density relation (red dots) obtained considering (A)
To assess the reliability of results, the optimal simulated speed-density relation is also compared to experimental measurements of pedestrian flows reported in literature. In particular, results presented by Oeding (1963), Older (1968), Mōri and Tsukaguchi (1987), Fruin (1987), Weidmann (1993) and Seyfried et al. (2005) are considered. Note that the results presented by Weidmann (1993) represent a fitting of the experiments and they are the basis from which the analytical model of Eq. 8 was obtained. Calibration results are in line with the experiments although the high variability of these latter.
Finally, the effect of the footbridge geometry on the simulated speed-density relation is evaluated. Results obtained considering a footbridge width of 2.5, 4.5 and 7.0 m are compared to those obtained from the reference width of 3.0 m in Figure 8B. It is firstly observed that an increase in the footbridge width does not affect the results while a narrower footbridge implies a slightly lower curve. This is because in the second case the repulsive effects of borders further limit the pedestrian density on the footbridge. On the contrary, increasing the footbridge geometry implies that the same mean pedestrian densities and velocities can be reached. Accordingly, the footbridge length does not affect the simulated speed-density relation, unless very short and unusual footbridges are considered. However, for each geometry variation, the pedestrian group size needed to reach a certain pedestrian density has to be evaluated as well as an adequate length of the access route and way out.
FIGURE 8. (A) Comparison between the simulated speed-density relation (red curve) and experimental measurements of pedestrian flows. Black dots: Weidmann (1993), green dots: Oeding (1963), grey dots: Mōri and Tsukaguchi (1987), orange crosses: Older (1968), light blue dots: Seyfried et al. (2005), blue circles: Fruin (1987)(B) Simulated speed-density relations obtained considering different footbridge width: 3.0 m (red curve), 2.5 m (black curve), 4.5 m (blue curve) and 7.0 m (magenta curve).
5 Pedestrian Flow Simulations
This section presents the footbridge vibrations induced by pedestrian flows simulated through the calibrated social force model.
5.1 Footbridge Parameters
The footbridge considered in the simulations is a simply supported beam 10 m long and 3 m wide having a linear dynamic behavior. It is characterized by a modal mass of 25 × 103 kg, a damping ratio of 0.5% and a natural frequency varying from 0.5 to 3.0 Hz. Only the contribution of the fundamental mode, namely the first bending mode having a half-sine mode shape, is taken into account.
5.2 Crowd Loading
The crowd loading is given by the superposition of the walking forces induced by each pedestrian in the crowd. The single pedestrian walking force is obtained as a series of successive footfall forces, each one described as a Fourier series according to the model proposed by Li et al. (2010). In particular the jth footstep force of a generic pedestrian α is expressed:
where
The next scheme outlines how the walking force of a generic pedestrian α in the crowd is evaluated:
•The pedestrian weight
•The trajectory of pedestrian α is evaluated from the results of the SFM, as the one shown in Figure 9 (black line).
•Application points of the single foot forces are identified on the pedestrian trajectory placing them at a distance of
•As the contribution of the first bending mode is taken into account, the 2-D problem can be reduced to a 1-D problem by projecting the foot standing points on the footbridge centerline (green dots in Figure 9). Note that this simplification would not be appropriate if a torsional mode contributed to the structural response.
•Instants of application of the single foot forces are evaluated as the instants when pedestrian α occupies the positions identified by red dots in Figure 9. Note that the initial time when the pedestrian starts to walk on the footbridge depends on its initial position as well as the path followed in the 60 m preceding the footbridge (Figure 2). Moreover, this method allows accounting for the fact that a pedestrian may stop walking and then restart. In this case, in fact, a larger time gap between two successive steps would be observed.
•The velocity in correspondence of each step is obtained by combining the instantaneous velocities in X and Y directions. Figure 10A presents the time varying position (in X direction) and velocity of the same pedestrian of Figure 9. Blue and red dots identify velocities and positions of the considered footsteps along with the corresponding instants of application.
•For each footstep j, the velocity magnitude is converted into a pacing frequency through the following relation proposed by Bruno and Venuti (2008) based on the experimental data presented in Bertram and Ruina (2001):
•where
•The period
•Starting from the pedestrian weight
•Finally, the single pedestrian walking force is obtained by adding the contribution of all footstep forces, each one applied in the position and at the time instant previously evaluated.
FIGURE 9. Example of pedestrian trajectory (black line), foot positions on the trajectory (red dots) and projection of the foot positions on the bridge centerline (green dots).
FIGURE 10. Example of (A) pedestrian velocities (blue line and dots) and positions in X direction (red line and dots) vs. time and (B) consecutive footfall forces.
Figure 10B shows the temporal sequence of footstep forces for the pedestrian α, which have to be added to obtain the single pedestrian walking force. This procedure is repeated for each pedestrian in the crowd and, finally, the crowd loading is obtained by the superposition of the individual pedestrian forces.
5.3 Footbridge Response Simulation
The footbridge response is evaluated by numerically integrating the equation of motion in the modal space considering the footbridge as a single degree of freedom system and time increments of 0.001 s (see, for instance, Bassoli et al. (2018a)). To calculate the modal force, the crowd loading of Section 5.2 is weighted by the amplitude of the mode shape in correspondence of the footfall positions. The variable representative of the response is chosen as the mid-span acceleration. The vibration response is assessed using a 1-s root-mean-square (RMS) moving average value from the acceleration time history.
Finally, the effects of the crowd-bridge interaction in terms of added mass are evaluated. In particular, the mean number of pedestrians simultaneously present on the footbridge during the time interval T1 − T2 (section 3.2) is added to the footbridge mass and the corresponding reduction in the footbridge natural frequency is estimated.
5.4 Results and Discussion
In the following, results obtained for pedestrian groups of 100 and 350 pedestrians are presented. Table 3 lists the values of the mean density and velocity obtained from criterion 2 (see section 3.2) for each of the 10 analyses performed, together with the mean walking frequency estimated from Eq. 19. Simulations with 100 pedestrians lead to mean pedestrian densities ranging from 0.51 to 0.84 ped/m2 and corresponding velocities from 1.20 to 1.02 m/s2, while 350 pedestrians involve densities in the range [1.61–1.90] ped/m2 and velocities from 0.85 to 0.78 m/s2. Velocities from criterion 2 give an indication of the pedestrian velocity. However, in the crowd load simulations, the velocity of each pedestrian is evaluated on a step-by-step basis as detailed in section 5.2. Note that group dimensions of 100 and 350 pedestrians are chosen as representative of correlated and very dense traffic, respectively. In the first case pedestrian movements are not completely free but partially constrained by the presence of other pedestrians, while in the second case they are highly constrained by other pedestrians. Moreover, it has to be stressed that the Fourier coefficients of Eq. 18 are defined for step frequencies in the range [1.60–2.40] Hz. Nevertheless, numerical simulations are performed considering also frequencies out of this range, as a consequence of the velocities involved.
The maximum RMS accelerations obtained for pedestrian groups of 100 and 350 pedestrians and for footbridge natural frequencies in the range [0.5–3.0] Hz are shown in Figure 11. In particular, for each group size and for each natural frequency the footbridge response is evaluated 10 times starting from the results of the 10 SFM simulations. These results are indicated with black dots in Figure 11, while mean values are highlighted in red. As expected, accelerations caused by 100 pedestrians are lower than those due to 350 pedestrians. The main amplifications of the structural response are obtained for frequencies between 1.5 and 2 Hz for the case of 100 pedestrians, while they are observed also for lower structural frequencies in the second case. This is an effect of the lower pedestrian velocities (and consequently step frequencies) characterizing higher pedestrian densities. Moreover, it is observed that, in both cases, amplifications of the structural response are obtained for a wide range of footbridge frequencies because each pedestrian in the crowd contributes differently to the modal excitation depending on his own pacing frequency, which changes step-by-step.
FIGURE 11. Maximum RMS accelerations caused by (A) 100 and (B) 350 pedestrians for different footbridge natural frequencies: results for each of the 10 analyses (black dots) and mean values (red line). Green dotted line: mean accelerations obtained accounting for the added pedestrian mass.
Finally, the effects of the pedestrian mass are evaluated. The mean number of pedestrians simultaneously present on the footbridge during the time interval
6 Comparison
Results of section 5 are compared with those obtained from the multiplication factor approach proposed by Caprani et al. (2012). In particular, they evaluate a set of enhancement factors for predicting the response due to a crowd based on the predicted accelerations of a single pedestrian. Enhancement factors are evaluated for crowd densities ranging from 0.44 ped/m2 to 2.11 ped/m2 (namely 0.44, 0.55, 0.75, 1.5 and 2.11 ped/m2) and synchronization proportions of 0,
In this study, the single pedestrian load is based on the single step load model of Li et al. (2010) (Eq. 18). In contrast to the crowd loading of section 5.2, in this case it is assumed that the pedestrian produces the same footfall force at each step leading to a periodic walking force. In line with Caprani et al. (2012), the pedestrian weight G and step length
Multiplication factors are plotted in Figure 12 together with the values of m proposed by Caprani et al. (2012). The multiplication factor obtained from each simulation is associated to the corresponding mean crowd density listed in Table 3. Results from 350 pedestrians, representing crowd densities typical of constrained traffic, are mainly consistent with those proposed by Caprani et al. (2012) for synchronization levels from 0.5N to N. On the contrary, results from 100 pedestrians, describing densities typical of correlated but not constrained traffic, are in quite good agreement with those of Caprani et al. (2012) for synchronization proportion lower than 0.2N. This findings are in line with the fact that pedestrians tend to synchronize more and more with increasing crowd density.
FIGURE 12. Multiplication factors m obtained from 100 pedestrians (black markers) and 350 pedestrians (blue markers) for footbridge natural frequencies of 1.94 Hz (.), 2.0 Hz (+) and 2.1 Hz (×). Colored lines: multiplication factors of Caprani et al. (2012) for different synchronization proportions.
7 Conclusions
This paper proposes a parameter calibration of the Helbing’s social force model (Helbing et al., 2005) performed adopting the response surface methodology. Model parameters are calibrated so as to describe the fundamental relation between mean walking speed and density of pedestrian crowds. The calibrated social force model enables a detailed simulation of unidirectional pedestrian flows on footbridges suitable for vibration assessment purposes. The updated model is adopted in this paper to simulate pedestrian flows and results are compared with the multiplication factor approach proposed by Caprani et al. (2012).
The speed-density relation obtained from the calibrated social force model shows a pretty good match with the theoretical curve, with larger discrepancies observed for ρ < 0.3 ped/m2. However, the simulated curve obtained considering
The social force model allows for the evaluation of the instantaneous position and velocity of each pedestrian in the crowd, assumed as the input of a step-by-step simulation of the pedestrian loads. Simulations presented in the paper highlight the advantages of a discrete crowd model, that is the direct possibility of modeling the variability of the crowd loading and the suitability for low and medium traffic densities and discontinuous pedestrian distribution. being the traffic continuity not required.
The main drawback of the presented crowd model is that analyses can be time-consuming, especially for high pedestrian densities. Moreover, the random nature of the crowd loading requires a statistical characterization of the response, incrementing the computing time.
It is worth noting that the proposed calibrated parameters are adequate to simulate pedestrian flows on footbridge with different dimensions than those adopted in this study, with exception of width lower than 3.0 m. In these case, in fact, the simulated speed-density relation is slightly lower. However, for each geometry variation, the pedestrian group size needed to reach a certain pedestrian density has to be evaluated as well as an adequate length of the access route and way out.
The comparison with the multiplication factors proposed by Caprani et al. (2012) shows a good agreement with the simulation results, even though they are scattered. In particular, the coherence of multiplication factors obtained for high pedestrians densities with those of Caprani et al. (2012) for high levels of synchronization, and vice versa for low pedestrian densities, demonstrates the ability of the model to describe the crowd behavior without the need to define in advance the synchronization level.
Data Availability Statement
The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.
Author Contributions
Conceptualization, EB and LV; Data collection, EB and LV; Analysis and interpretation of results, EB and LV; Supervision, LV; Writing, review and editing, EB and LV.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Barela, A. M. F., and Duarte, M. (2008). Biomechanical Characteristics of Elderly Individuals Walking on Land and in Water. J. Electromyogr. Kinesiol. 18, 446–454. doi:10.1016/j.jelekin.2006.10.008
Bassoli, E., Gambarelli, P., and Vincenzi, L. (2018a). Human-induced Vibrations of a Curved Cable-Stayed Footbridge. J. Constructional Steel Res. 146, 84–96. doi:10.1016/j.jcsr.2018.02.001
Bassoli, E., Van Nimmen, K., Vincenzi, L., and Van den Broeck, P. (2018b). A Spectral Load Model for Pedestrian Excitation Including Vertical Human-Structure Interaction. Eng. Structures 156, 537–547. doi:10.1016/j.engstruct.2017.11.050
Bertram, J. E. A., and Ruina, A. (2001). Multiple Walking Speed-Frequency Relations Are Predicted by Constrained Optimization. J. Theor. Biol. 209, 445–453. doi:10.1006/jtbi.2001.2279
Box, G. E. P., and Wilson, K. B. (1951). On the Experimental Attainment of Optimum Conditions. J. R. Stat. Soc. Ser. B (Methodological). 13, 1–38. doi:10.1111/j.2517-6161.1951.tb00067.x
Bruno, L., and Venuti, F. (2008). The Pedestrian Speed-Density Relation: Modelling and Application. in Proceedings Of the 3rd International Footbridge Conference, Porto, Portugal
Buchmüller, S., and Weidmann, U. (2006). Parameters of Pedestrians, Pedestrian Traffic and Walking Facilities. Tech. rep., ETH Zürich, IVT. 132, 1–58. doi:10.3929/ethz-b-000047950
Butz, C., Feldmann, M., Heinemeyer, C., and Sedlacek, G. (2008). Advanced Load Models for Synchronous Pedestrian Excitation and Optimised Design Guidelines for Steel Footbridges. Tech. rep., Res. Fund Coal Steel
Caprani, C. C., Keogh, J., Archbold, P., and Fanning, P. (2012). Enhancement Factors for the Vertical Response of Footbridges Subjected to Stochastic Crowd Loading. Comput. Structures. 102-103–103, 87–96. doi:10.1016/j.compstruc.2012.03.006
Carroll, S. P., Owen, J. S., and Hussein, M. F. M. (2012). Modelling Crowd-Bridge Dynamic Interaction with a Discretely Defined Crowd. J. Sound Vibration 331, 2685–2709. doi:10.1016/j.jsv.2012.01.025
Chen, X., Treiber, M., Kanagaraj, V., and Li, H. (2018). Social Force Models for Pedestrian Traffic - State of the Art. Transport Rev. 38, 625–653. doi:10.1080/01441647.2017.1396265
Cristiani, E., Piccoli, B., and Tosin, A. (2011). Multiscale Modeling of Granular Flows with Application to Crowd Dynamics. Multiscale Model. Simul. 9, 155–182. doi:10.1137/100797515
da Silva, F. T., Brito, H. M. B. F., and Pimentel, R. L. (2013). Modeling of Crowd Load in Vertical Direction Using Biodynamic Model for Pedestrians Crossing Footbridges. Can. J. Civ. Eng. 40, 1196–1204. doi:10.1139/cjce-2011-0587
Dallard, P., Fitzpatrick, T., Bourva, S. L., Low, A., Smith, R. R., and Willford, M. (2001). The London Millennium Footbridge. The Struct. Engineer. 79, 17–33.
Danbon, F., and Grillaud, G. (2005). Dynamic Behaviour of a Steel Footbridge. Characterization and Modelling of the Dynamic Loading Induced by a Moving Crowd on the Solferino Footbridge in Paris. in Proceedings Of the 2nd International Footbridge Conference (Venice, Italy)
Ferrarotti, A., and Tubino, F. (2015). Equivalent Spectral Model for Pedestrian-Induced Forces on Footbridges: a Generalized Formulation. in Proceedings Of the 5th ECCOMAS Thematic Conference On Computational Methods In Structural Dynamics And Earthquake Engineering (Crete Island, Greece). doi:10.7712/120115.3551.800
Fujino, Y., Pacheco, B. M., Nakamura, S.-I., and Warnitchai, P. (1993). Synchronization of Human Walking Observed during Lateral Vibration of a Congested Pedestrian Bridge. Earthquake Engng. Struct. Dyn. 22, 741–758. doi:10.1002/eqe.4290220902
Hariri-Ardebili, M. A., Seyed-Kolbadi, S. M., and Noori, M. (2018). Response Surface Method for Material Uncertainty Quantification of Infrastructures. Shock and Vibration. 2018, 1, 14. doi:10.1155/2018/1784203
Helbing, D., Buzna, L., Johansson, A., and Werner, T. (2005). Self-organized Pedestrian Crowd Dynamics: Experiments, Simulations, and Design Solutions. Transportation Sci. 39, 1–24. doi:10.1287/trsc.1040.0108
Helbing, D., Farkas, I., and Vicsek, T. (2000). Simulating Dynamical Features of Escape Panic. Nature. 407, 487–490. doi:10.1038/35035023
Helbing, D., and Johansson, A. (2009). Pedestrian, Crowd and Evacuation Dynamics. New York, NY: Springer New York, 6476–6495. doi:10.1007/978-0-387-30440-3_382
Helbing, D., and Molnár, P. (1995). Social Force Model for Pedestrian Dynamics. Phys. Rev. E. 51, 4282–4286. doi:10.1103/PhysRevE.51.4282
Jiménez-Alonso, J. F., Sáez, A., Caetano, E., and Magalhães, F. (2016). Vertical Crowd-Structure Interaction Model to Analyze the Change of the Modal Properties of a Footbridge. J. Bridge Eng. 21. doi:10.1061/(ASCE)BE.1943-5592.0000828
Jiménez-Alonso, J. F., and Sáez, A. (2017). Recent Advances in the Serviceability Assessment of Footbridges under Pedestrian-Induced Vibrations. Bridge Eng. Chapter 5. doi:10.5772/intechopen.71888
Johansson, A., Helbing, D., and Shukla, P. K. (2007). Specification of the Social Force Pedestrian Model by Evolutionary Adjustment to Video Tracking Data. Advs. Complex Syst. 10, 271–288. doi:10.1142/S0219525907001355
Kaymaz, I., and McMahon, C. A. (2005). A Response Surface Method Based on Weighted Regression for Structural Reliability Analysis. Probabilistic Eng. Mech. 20, 11–17. doi:10.1016/j.probengmech.2004.05.005
Khuri, A. I., and Cornell, J. A. (1996). Response Surfaces: Designs and Analyses, 2nd Edition (Boca Raton: CRC Press). doi:10.1201/9780203740774
Langston, P. A., Masling, R., and Asmar, B. N. (2006). Crowd Dynamics Discrete Element Multi-Circle Model. Saf. Sci. 44, 395–417. doi:10.1016/j.ssci.2005.11.007
Li, Q., Fan, J., Nie, J., Li, Q., and Chen, Y. (2010). Crowd-induced Random Vibration of Footbridge and Vibration Control Using Multiple Tuned Mass Dampers. J. Sound Vibration. 329, 4068–4092. doi:10.1016/j.jsv.2010.04.013
Mōri, M., and Tsukaguchi, H. (1987). A New Method for Evaluation of Level of Service in Pedestrian Facilities. Transportation Res. A: Gen. 21, 223–234. doi:10.1016/0191-2607(87)90016-1
Matsumoto, Y., Nishioka, T., Shiojiri, H., and Matsuzaki, K. (1978). Dynamic Design of Footbridges. IABSE Proc. 1–15. doi:10.5169/seals-33221
Moreu, F., Maharjan, D., Zhu, C., and Wyckoff, E. (2020). Monitoring Human Induced Floor Vibrations for Quantifying Dance Moves: A Study of Human-Structure Interaction. Front. Built Environ. 6, 36. doi:10.3389/fbuil.2020.00036
Myers, R. H., Montgomery, D. C., and Anderson-Cook, C. M. (2016). Response Surface Methodology: Process and Product Optimization Using Designed Experiments. 4th Edition. Wiley Series in Probability and Statistics (Wiley)
Oeding, D. (1963). Verkehrsbelastung und Dimensionierung von Gehwegen und anderen Anlagen des Fußgängerverkehrs (Strassenbau und Strassenterkehrstechnik). Bonn: German
Older, S. (1968). Movement of Pedestrians on Footways in Shopping Streets. Traffic Eng. Control. 10, 160–163.
Pan, S., Xu, S., Mirshekari, M., Zhang, P., and Noh, H. Y. (2017). Collaboratively Adaptive Vibration Sensing System for High-Fidelity Monitoring of Structural Responses Induced by Pedestrians. Front. Built Environ. 3, 28. doi:10.3389/fbuil.2017.00028
Parry, G. W. (2007). The Dynamics Of Crowds. Master’s thesis. Dordrecht: Department of Mathematical Sciences - University of Bath. doi:10.1007/1-4020-5312-6
Portier, K., Keith Tolson, J., and Roberts, S. M. (2007). Body Weight Distributions for Risk Assessment. Risk Anal. 27 (1), 11–26. doi:10.1111/j.1539-6924.2006.00856.x
Pula, W., and Bauer, J. (2007). Probabilistic Methods in Geotechnical Engineering (Vienna: Springer), Vol. 491 of CISM Courses And Lectures, Chap. Application of the Response Surface Method. doi:10.1007/978-3-211-73366-0_6
Racic, V., Pavic, A., and Brownjohn, J. M. W. (2009). Experimental Identification and Analytical Modelling of Human Walking Forces: Literature Review. J. Sound Vibration. 326, 1–49. doi:10.1016/j.jsv.2009.04.020
Seyfried, A., Steffen, B., Klingsch, W., and Boltes, M. (2005). The Fundamental Dyagram of Pedestrian Movements Revisited. J. Stat. Mech. 10, P10002. doi:10.1088/1742-5468/2005/10/p10002
Van Nimmen, K., Lombaert, G., Jonkers, I., De Roeck, G., and Van den Broeck, P. (2014). Characterisation of Walking Loads by 3D Inertial Motion Tracking. J. Sound Vibration. 333, 5212–5226. doi:10.1016/j.jsv.2014.05.022
Van Nimmen, K., Van den Broeck, P., Lombaert, G., and Tubino, F. (2020). Pedestrian-induced Vibrations of Footbridges: An Extended Spectral Approach. J. Bridge Eng. 25, 04020058. doi:10.1061/(ASCE)BE.1943-5592.0001582
Venuti, F., and Bruno, L. (2007). An Interpretative Model of the Pedestrian Fundamental Relation. Comptes Rendus Mécanique. 335, 194–200. doi:10.1016/j.crme.2007.03.008
Venuti, F., Racic, V., and Corbetta, A. (2016). Modelling Framework for Dynamic Interaction between Multiple Pedestrians and Vertical Vibrations of Footbridges. J. Sound Vibration 379, 245–263. doi:10.1016/j.jsv.2016.05.047
Vincenzi, L., and Savoia, M. (2015). Coupling Response Surface and Differential Evolution for Parameter Identification Problems. Computer-Aided Civil Infrastructure Eng. 30, 376–393. doi:10.1111/mice.12124
Wang, Y., Brownjohn, J., Dai, K., and Patel, M. (2019). An Estimation of Pedestrian Action on Footbridges Using Computer Vision Approaches. Front. Built Environ. 5, 133. doi:10.3389/fbuil.2019.00133
Wei, X., Broeck, P. V. d., Roeck, G. D., and Nimmen, K. V. (2017). A Simplified Method to Account for the Effect of Human-Human Interaction on the Pedestrian-Induced Vibrations of Footbridges. Proced. Eng. 199, 2907–2912. doi:10.1016/j.proeng.2017.09.331International Conference on Structural Dynamics, EURODYN 2017
Younis, A., Avci, O., Hussein, M., Davis, B., and Reynolds, P. (2017). Dynamic Forces Induced by a Single Pedestrian: A Literature Review. Appl. Mech. Rev. 69. doi:10.1115/1.4036327
Živanović, S., Pavić, A., and Ingólfsson, E. (2010). Modeling Spatially Unrestricted Pedestrian Traffic on Footbridges. J. Struct. Eng. ASCE 136, 1296–1308. doi:10.1061/(ASCE)ST.1943-541X.0000226
Keywords: social force model, crowd loading, model updating, response surface, footbridges
Citation: Bassoli E and Vincenzi L (2021) Parameter Calibration of a Social Force Model for the Crowd-Induced Vibrations of Footbridges. Front. Built Environ. 7:656799. doi: 10.3389/fbuil.2021.656799
Received: 21 January 2021; Accepted: 28 April 2021;
Published: 17 May 2021.
Edited by:
Onur Avci, Iowa State University, United StatesReviewed by:
Chuanzhi Dong, University of Central Florida, United StatesSin-Chi Kuok, University of Cambridge, United Kingdom
Copyright © 2021 Bassoli and Vincenzi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Loris Vincenzi, loris.vincenzi@unimore.it