Equivalent Dynamic Load Factor of Different Non-Exceedance Probability for Crowd Jumping Loads

: Existing dynamic load factors (DLF) for crowd jumping loads are modeled by identifying the peaks and energies from the corresponding Fourier amplitude spectrum, which may underestimate and overestimate structural responses, respectively. Based on the principle of equal structural responses, this study herein develops an equivalent DLF, using the frequency response function to weight and integrate the power spectral density (PSD) matrix of crowd jumping loads. Firstly, massive PSD matrices of different crowd sizes and different metronome frequencies are simulated using a random field model of crowd jumping loads. Thereafter, the equivalent DLF of different structural damping ratios, crowd sizes, metronome frequencies, and non-exceedance probabilities are investigated, and a model of the equivalent DLF is established accordingly. It is believed that this model could provide an efficient and accurate way to conduct reliability design for structures subjected to crowd jumping loads.


Introduction
Currently, the modern assembly structures (e.g., large-span floor slabs, footbridges, grandstands, etc.), gradually exhibit the characteristics of low natural frequencies and low damping ratios, due to the extensive use of lightweight, high-strength building materials [1][2][3][4][5]. Accordingly, structural vibration serviceability issues [6,7] as well as safety risks [8,9] caused by crowd activities have been increasingly emerging, which may further lead to huge economic losses and even heavy casualties in extreme cases. One typical example of this was the London Millennium Bridge in 2000, in which severe lateral vibrations of the bridge were induced by pedestrian use [10]. The bridge was therefore closed to traffic for three days after the grand opening, and it cost GBP 5 million to solve the vibration problem. Additionally, at a sports concert given in Vigo, Spain in 2018, the wooden stands collapsed due to the audience's rhythmic activities, with over 300 people injured [11]. Jumping is usually considered to generate the most significant vertical loads in contrast to other common activities [12,13], in particular when the crowd is coordinated by an external beat. This will be the focus of this research.
Load models are necessary, both in the design stage and the as-built stage to evaluate the vibration performance of structures subjected to crowd jumping activities [14]. Considering that jumping activities are approximate periodic, Fourier series models are mostly used to represent jumping loads [13,[15][16][17][18][19], where the dynamic load factor (DLF) is the key parameter to characterize the load magnitude and crowd synchronization. When crowds jump to an external beat, it is unlikely to achieve perfect synchronization because of both inter-and intra-subject variability [20], and the crowd behave differently to the beat of different dominant frequencies [17]. Therefore, the DLF is a stochastic variable and related to the crowd size and guidance frequency, where the guidance frequency denotes the dominant frequency of the metronome or music that guides the crowd to jump. Massive investigations on jumping loads have been conducted [13][14][15][16][17][18][19], and DLF models have been established accordingly. The DLF of these models is identified by the peaks [16] or energies [17] from the corresponding Fourier amplitude spectrum. The intra-subject variability of a jumping activity caused the load energy to spread around the jumping frequency and its multiplier [21], where the jumping frequency denotes the dominant frequency of a person's jumping load time history. When a crowd jumps at a given guidance frequency, most people can follow the beat, and their jumping frequency is equal to the guidance frequency.
As for a few people who cannot follow the beat, their jumping frequency is unequal to the guidance frequency. The DLF obtained by identifying the peaks can underestimate structural responses because it overlooks the energy near the jumping frequency and its multiplier [22]. The DLF obtained by identifying the energy can overestimate structural responses because it assumes that the energy is all concentrated at the jumping frequency and its multiplier. When using the resonance assumption to calculate structural responses, the energy at the jumping frequency and its multiplier contributes more to structural responses than the energy around them. Moreover, most of the DLF models lack the data support of large crowd jumping tests. For one thing, the DLF model is developed based on the individual jumping load records aligned by the time instant of the beat [17,19], in which the influence of human-human interaction is not involved. For another, the DLF is modeled through crowd jumping loads identified by structural responses [15], in which the jumping load of each person could not be recorded. Moreover, owing to the limitation of test conditions, the model ignores the change of the DLF with the guidance frequency. In addition, only the mean value or a certain quantile value of the DLF is provided, which is not applicable for the structural design of different reliabilities.
The accurate calculation of structural responses is the main purpose of load modeling. Hence, the modeling parameters could be obtained based on the principle of equal structural responses. This idea has been applied to the modeling of crowd walking loads, in which the number of equivalent synchronized pedestrians is obtained using this principle [23,24]. Similarly, the DLF can be identified with the equal structural response, hereafter termed as the equivalent DLF. This is the main highlight of this study. It can overcome the shortcomings of the existing DLF modeled by identifying the peaks or energies from the corresponding Fourier amplitude spectrum. Moreover, structural response could be calculated in the form of Fourier series using the equivalent DLF, which is further more convenient for application to structural design and evaluation than other complicated stochastic models of jumping loads [12,[25][26][27][28][29]. In addition, to further apply the model to structural reliability design, the equivalent DLF of different non-exceedance probability is established. This is another main highlight of this study. It can overcome the shortcomings of the existing DLF, which only provides its mean value or a certain quantile value. The paper begins with the formula of the equivalent DLF calculation, which adopts the frequency response function (FRF) to weight and integrate the power spectral density (PSD) matrix of crowd jumping loads. This approach makes the root mean square value of structural acceleration responses equal. The formula of simulating PSD matrices of crowd jumping loads is then provided. This is followed by the modeling of the equivalent DLF with four parameters, i.e., structural damping ratio, guidance frequency, crowd size, and non-exceedance probabilities. The Engineering application of the equivalent DLF is followed. Afterwards, discussions on this work are provided.

Formula of the Equivalent Dynamic Load Factor
The Fourier series-based model for crowd jumping loads, F(t), is usually expressed as [15]: where W is the total body weight of the jumping crowd, r j and ϕ j are the DLF and phase angle of the jth harmonic, f gu is the guidance frequency, and n h is the number of load harmonics considered. The values of r j for jumping loads are commonly determined by two methods. In the first method, hereafter termed as the peak method, the maximum value of the jth harmonic region in the Fourier amplitude spectrum is adopted as r j , as depicted in Figure 1, and the values are listed in Table 1 for further comparison. In the above, the jth harmonic region is ((j − 0.5)f gu , (j + 0.5)f gu ). As for the second method, hereafter termed as the energy method, the total energy of the jth harmonic region in the Fourier amplitude spectrum is adopted as r j , as calculated by: where A jh is the amplitude of the hth sinusoid in the jth harmonic region and n s is the number of the sinusoids. The values of r j calculated by the energy method for the Fourier amplitude spectrum in Figure 1 are provided in Table 1. The peak method only considers the energy at the jumping frequency and its multiplier, and the energy in between is ignored, which can underestimate structural responses. The energy method assumes that all the energy of a jumping load is concentrated at the jumping frequency and its multiplier, which overestimates structural responses when the resonance assumption is made.
as [15]: where W is the total body weight of the jumping crowd, rj and φj are the DLF and phase angle of the jth harmonic, fgu is the guidance frequency, and nh is the number of load harmonics considered. The values of rj for jumping loads are commonly determined by two methods. In the first method, hereafter termed as the peak method, the maximum value of the jth harmonic region in the Fourier amplitude spectrum is adopted as rj, as depicted in Figure 1, and the values are listed in Table 1 for further comparison. In the above, the jth harmonic region is ((j − 0.5)fgu, (j + 0.5)fgu). As for the second method, hereafter termed as the energy method, the total energy of the jth harmonic region in the Fourier amplitude spectrum is adopted as rj, as calculated by: where Ajh is the amplitude of the hth sinusoid in the jth harmonic region and ns is the number of the sinusoids. The values of rj calculated by the energy method for the Fourier amplitude spectrum in Figure 1 are provided in Table 1. The peak method only considers the energy at the jumping frequency and its multiplier, and the energy in between is ignored, which can underestimate structural responses. The energy method assumes that all the energy of a jumping load is concentrated at the jumping frequency and its multiplier, which overestimates structural responses when the resonance assumption is made.  A novel DLF calculating method that can accurately estimate the root mean square value of structural acceleration responses is therefore proposed. Considering a single degree-of-freedom system with a unit mass subjected to crowd jumping loads, as shown in Figure 2, the equation of motion of the system is expressed as:  A novel DLF calculating method that can accurately estimate the root mean square value of structural acceleration responses is therefore proposed. Considering a single degree-of-freedom system with a unit mass subjected to crowd jumping loads, as shown in Figure 2, the equation of motion of the system is expressed as: ..
where ζ s , ω s , and u s (t) are the damping ratio, circular frequency, and displacement of the system, respectively, W l is the body weight of the lth jumping person, and x l (t) is the ground reaction force time history of the lth jumping person divided by W l , which is dimensionless. Therefore, x l (t) reflects the change of center when an individual jumps, and it is the most direct variable to reflect the jumping pattern. The FRF H(f ; f s , ζ s ) related to the input, where f s = 2πω s is the cyclic frequency of the structure and i is the imaginary unit. Then, the auto-PSD [21] G(f ) of .. u s (t) can be obtained by the stochastic vibration theory [30]: where G lk (f ) is the lth row and kth column element of the PSD matrix of crowd jumping loads. If l = k, G lk (f ) is the cross-PSD of x l (t) and x k (t). If l = k, G ll (f ) is the auto-PSD of x l (t).
Integrating G(f ) within the frequency range yields the root mean square value a rms of .. u s (t): dimensionless. Therefore, xl(t) reflects the change of center when an individual jumps, and it is the most direct variable to reflect the jumping pattern.  i where fs = 2πωs is the cyclic frequency of the structure and i is the imaginary unit. Then, the auto-PSD [21] G(f) of ( )  s u t can be obtained by the stochastic vibration theory [30]: where Glk(f) is the lth row and kth column element of the PSD matrix of crowd jumping loads. If l ≠ k, Glk(f) is the cross-PSD of xl(t) and xk(t). If l = k, Gll(f) is the auto-PSD of xl(t).
Integrating G(f) within the frequency range yields the root mean square value arms of ( ) If the crowd jumping load, i.e., the right term of Equation (3), is modeled by the Fourier series, as expressed by Equation (1), the auto-PSD GF(f) of ( ) F t is: where ( ) δ is the Dirac delta function. Then, as the derivation procedure of Equations If the crowd jumping load, i.e., the right term of Equation (3), is modeled by the Fourier series, as expressed by Equation (1), the auto-PSD G F (f ) of F(t) is: where δ() is the Dirac delta function. Then, as the derivation procedure of Equations (3)-(5), a rms of .. u s (t) can be calculated by: The resonance scenario is commonly assumed for structural response estimation to consider the severest case [31], where structural frequency is assumed to be equal to the guidance frequency or its multiplier. Consequently, structural responses induced by the jth harmonic of jumping load is much larger than the rest of the harmonics: Making a rms calculated by Equations (6) and (9) equal and assuming a resonance scenario in Equation (6) yields: If the weight of every jumping person is assumed to be equal, Equation (10) will be further simplified as: In the above, H(jf gu ; jf gu , ζ s ) is the maximum of H(f ; jf gu , ζ s ), so it is indicated that the equivalent DLF for crowd jumping loads can be identified using the FRF to weight and integrate the PSD matrix of crowd jumping loads, hereafter termed as the FRF-weighting method. The equivalent DLF is related to the structural damping ratio because Equation (11) is derived from the principle of equal structural responses. Then, the equivalent DLF of the Fourier amplitude spectrum in Figure 1 calculated by Equation (11) is provided in Table 1. It is observed from Table 1 that the DLF obtained by the peak method and energy method are quite different, especially for high-order harmonics, because the degree of energy diffusion increases as the harmonic order increases [21]. The maximum-normalized FRF is adopted in the FRF-weighting method to weight the contribution of different sinusoids to structural responses, and the equivalent DLF obtained by this method is between the DLF obtained by the peak method and the DLF obtained by the energy method. Moreover, the equivalent DLF obtained by the FRF-weighting method is strongly related to structural damping ratio. When the structural damping ratio is small, the bandwidth of the FRF is narrow, so the energy around the jumping frequency and its multiplier contributes less to structural responses. Therefore, the equivalent DLF is small and close to the DLF obtained by the peak method. When the structural damping ratio is large, the bandwidth of the FRF is wide and the energy around the jumping frequency and its multiplier contributes more to structural responses. Therefore, the equivalent DLF is large and close to the DLF obtained by the energy method.

Simulation of the Equivalent Dynamic Load Factor
The PSD matrix for crowd jumping loads is fundamental to calculating the equivalent DLF according to Equation (11). In this section, a random field model for crowd jumping loads proposed by the authors [32] is adopted for the simulation of the PSD matrix.

Random Field Model for Crowd Jumping Loads
The auto-PSD of individual jumping loads, i.e., the diagonal element of the random field model, is established according to the individual jumping load records collected by force plates. The cross-PSD of crowd jumping load, i.e., the off-diagonal element of the random field model, is established according to the crowd jumping records collected by three-dimensional motion capture technology. It is proved by spectral analysis that the energy after the second harmonic of the load is weak, and the test error is large for high-frequency components due to the high-frequency trembling of the skin when threedimensional motion capture technology is used. Therefore, the first two harmonics of crowd jumping loads are considered in the random field model. Spectral analysis of the individual load records demonstrates that the energy of the jumping load is mainly distributed around the fundamental jumping frequency and its multiple integers (harmonics), and the energy in between is almost equal to zero. Therefore, the auto-PSD can be modeled separately for each harmonic. The auto-PSD sample is then normalized by the jumping frequency and the harmonic order, and a two-term Gaussian function is used to model the shape of the normalized auto-PSD. Finally, the lth diagonal element of the PSD matrix, i.e., the auto-PSD G ll (f ) of the lth person's jumping load, is modeled by: where G ll,j (f ; f ju ) is the jth harmonic of the auto-PSD and s j (f ju ) is the energy for the jth harmonic, which can be calculated by Equation (13); the coefficients identified by the test records are listed in Table 2. ρ(f ju ) is used to reflect the energy truncation effect and can be calculated by Equation (14). G no (f /jf ju ) is used to describe the shape of the normalized-PSD, which can be calculated by Equation (15), and the coefficients are provided in Table 2. The coherence function is the key to the modeling of the cross-PSD, because the auto-PSD has been well modeled by Equation (12). First, the coherence function is modeled by filtering the records with a dominant frequency equal to the metronome frequency. The magnitude of the coherence function is then simplified as an energy scaling factor, because the real part and imaginary part of cross-PSD have the same curve shape, as auto-PSD. The phase angle of the coherence function is simplified as a set of time lag shifts, because different load records have the same dominant frequency. Finally, the off-diagonal element of the PSD matrix, i.e., the cross-PSD G lk (f ) of the lth person's jumping load and kth person's jumping load is modeled by: where f ju,l is the lth person's jumping frequency, G ll,j (f ; f ju,l ) is defined by Equation (12), and |γ lk (j; f ju,l )| is the magnitude of the coherence function at the jth harmonic, as provided in Table 3. ∆t l (f ju,l ) denotes the time lag shift, which follows a normal distribution N(0, σ 2 ) [s], and the value of σ is provided in Table 3. When a crowd jumps at a given guidance frequency, some people can follow the beat while others cannot. Therefore, the jumping frequency of each person, i.e., f ju,l , is modelled as a random variable, and Table 4 in the literature [32] shows its probability distribution.

Simulation Procedure
Because the first two harmonics of crowd jumping loads are considered in the random field model, the equivalent DLF for the first two harmonics is simulated in this section. As detailed in Table 1, the equivalent DLF is relevant to the structural damping ratio, and the structural damping ratio of common structures subjected to crowd jumping loads, such as gymnasia, grandstands, and concert halls, is in the range of 0.005-0.05 [31,33]. Therefore, the structural damping ratio is taken from 0.005 to 0.05 with an interval of 0.005. Moreover, a crowd reacts differently to different guidance frequencies, so it is crucial to investigate the relationship between the equivalent DLF and the guidance frequency. The guidance frequency is taken from 1.5 to 3.5 Hz, with an interval of 0.1 Hz, because it is hard for a crowd to keep jumping outside this range [32]. The DLF and crowd size exhibit a nonlinear relationship because of the inter-subject variability. In addition, crowd coordination tends to be a constant when the crowd size exceeds fifty [15,18], so the crowd size is taken from 1 to 50 with an interval of 1. To further apply the model to structural reliability design, 1000 Monte Carlo simulations are performed for each working condition; to investigate the equivalent DLF of different non-exceedance probability, the adequacy of the number of Monte Carlo simulations will be discussed in Section 4.4. Details of the simulation algorithm can be seen in Algorithm 1.

Modeling of the Equivalent Dynamic Load Factor
The equivalent DLF is related to the structural damping ratio, guidance frequency, crowd size, and non-exceedance probability, so the modeling of the equivalent DLF is a multivariable regression problem. The working condition, where the crowd size is equal to 50, the structural damping ratio is equal to 0.005, and the non-exceedance probability is equal to 0.5, is defined as the standard working condition. First, the equivalent DLF model of the standard working condition is built. Afterwards, the conversion models for other crowd size, structural damping ratio, and non-exceedance probability are established.

Modeling of the Equivalent DLF for the Standard Working Condition
According to the definition of the standard working condition, the median of 1000 simulation data with a structural damping ratio equal to 0.005 and crowd size equal to 50 are adopted to build the model. The crowd reacts differently to different beats. This is also proven by the probability distribution of the jumping frequency, the probability distribution of the time lag shift, and the values of the coherency function magnitude, all of which are related to the guidance frequency. Therefore, the equivalent DLF of the standard working condition is related to the guidance frequency, as depicted in Figure 3, from which it is observed that the equivalent DLF increases first and then decreases with the guidance frequency. This is because it is hard for a crowd to keep up with the beat, and the jumping activity is exhausted when the guidance frequency is high or low. Particularly, the equivalent DLF at the low guidance frequency is only about one-third of the maximum value at the moderate guidance frequency. A fourth-order polynomial is used to model the equivalent DLF for the standard working condition, as depicted in Figure 3, in which Equation (17) is used for the first harmonic and Equation (18)  the jumping activity is exhausted when the guidance frequency is high or low. Particularly, the equivalent DLF at the low guidance frequency is only about one-third of the maximum value at the moderate guidance frequency. A fourth-order polynomial is used to model the equivalent DLF for the standard working condition, as depicted in Figure  3, in which Equation (17) is used for the first harmonic and Equation (18)

Conversion Model for Structural Damping Ratio
According to Equation (11), the structural damping ratio affects the equivalent DLF by changing the magnitude of the FRF. The bandwidth of the FRF increases as the structural damping ratio increases, so the contribution for the energy between different harmonics to structural responses increases. As a result, the equivalent DLF will increase with the structural damping ratio. Moreover, the influence degree of the structural damping ratio on the equivalent DLF is related to the PSD matrix. In particular, the impact degree increases as the bandwidth of the sum of the PSD matrix's elements increases, because more energies are distributed between the jumping frequency and its multiplier, and these energies are weighted by the FRF to contribute for structural responses. First, the bandwidth increases as the jumping frequency increases, which is illustrated by the expression of the normalized PSD, i.e., Equation (15). Moreover, this bandwidth increases as the dispersion degree of the jumping frequency increases according to Equation (16). The dispersion degree of the jumping frequency is large at the low guidance frequency [32], because many people cannot keep up with the beat. Consequently, the influence degree of the structural damping ratio on the equivalent DLF for different guidance frequencies is complicated. To further yield the influence degree, the median

Conversion Model for Structural Damping Ratio
According to Equation (11), the structural damping ratio affects the equivalent DLF by changing the magnitude of the FRF. The bandwidth of the FRF increases as the structural damping ratio increases, so the contribution for the energy between different harmonics to structural responses increases. As a result, the equivalent DLF will increase with the structural damping ratio. Moreover, the influence degree of the structural damping ratio on the equivalent DLF is related to the PSD matrix. In particular, the impact degree increases as the bandwidth of the sum of the PSD matrix's elements increases, because more energies are distributed between the jumping frequency and its multiplier, and these energies are weighted by the FRF to contribute for structural responses. First, the bandwidth increases as the jumping frequency increases, which is illustrated by the expression of the normalized PSD, i.e., Equation (15). Moreover, this bandwidth increases as the dispersion degree of the jumping frequency increases according to Equation (16). The dispersion degree of the jumping frequency is large at the low guidance frequency [32], because many people cannot keep up with the beat. Consequently, the influence degree of the structural damping ratio on the equivalent DLF for different guidance frequencies is complicated. To further yield the influence degree, the median of the 1000 simulations with a crowd size equal to 50 is used to plot Figure 4, where the equivalent DLF is normalized by the equivalent DLF with a structural damping ration equal to 0.005.  It is observed from Figure 4 that the influence of the structural damping ratio on the equivalent DLF is not greatly affected by the guidance frequency in general. Only at the low guidance frequency, such as 1.5 and 1.6 Hz, the influence is different and is slightly greater than the other guidance frequencies. Therefore, the influence of the guidance It is observed from Figure 4 that the influence of the structural damping ratio on the equivalent DLF is not greatly affected by the guidance frequency in general. Only at the low guidance frequency, such as 1.5 and 1.6 Hz, the influence is different and is slightly greater than the other guidance frequencies. Therefore, the influence of the guidance frequency is overlooked in the conversion model of Equations (19) and (20) for the structural damping ratio. In addition, the crowd size and non-exceedance probability are also overlooked in the conversion model of Equations (19) and (20), because these two variables are not related to the influence of the structural damping ratio on the equivalent DLF according to Equation (11). Then, by fitting the average curve of different guidance frequencies in Figure 4, power functions, as depicted in Figure 5, are adopted to establish the conversion model, as shown in Equations (19) and (20).
where c 1 (ζ s ) is the conversion model for the first harmonic and c 2 (ζ s ) is the conversion model for the second c 1 (ζ s ) = −0.5241ζ s −0.2253 + 2.728 harmonic. frequency is overlooked in the conversion model of Equations (19) and (20) for the structural damping ratio. In addition, the crowd size and non-exceedance probability are also overlooked in the conversion model of Equations (19) and (20), because these two variables are not related to the influence of the structural damping ratio on the equivalent DLF according to Equation (11). Then, by fitting the average curve of different guidance frequencies in Figure 4, power functions, as depicted in Figure 5, are adopted to establish the conversion model, as shown in Equations (19) and (20).

Conversion Model for Crowd Size
Owing to the inter-subject variability, the equivalent DLF is not proportional to the crowd size. In existing models, the crowd DLF divided by the crowd size is defined as the coordination factor to quantify the crowd synchronization [15], because the coordination factor tends to be a constant as the crowd size increases, which makes the modeling of the crowd DLF easy. Learning from this idea, the coordination factor, coi, in this section is defined as: where the subscript j denotes the harmonic order, cj is the conversion model for the crowd size, and rj(np) is the equivalent DLF with crowd size np. To further investigate the law of the coordination factor, the median of the 1000 simulations with structural damping ratio equal to 0.005 is used to plot Figure 6, where the ordinate is obtained by Equation (21).

Conversion Model for Crowd Size
Owing to the inter-subject variability, the equivalent DLF is not proportional to the crowd size. In existing models, the crowd DLF divided by the crowd size is defined as the coordination factor to quantify the crowd synchronization [15], because the coordination factor tends to be a constant as the crowd size increases, which makes the modeling of the crowd DLF easy. Learning from this idea, the coordination factor, co i , in this section is defined as: co j n p = c j n p n p = r j n p n p r j (50) (21) where the subscript j denotes the harmonic order, c j is the conversion model for the crowd size, and r j (n p ) is the equivalent DLF with crowd size n p . To further investigate the law of the coordination factor, the median of the 1000 simulations with structural damping ratio equal to 0.005 is used to plot Figure 6, where the ordinate is obtained by Equation (21). It is observed from Figure 6 that the influence of the crowd size on the coordination factor is affected by the guidance frequency when the crowd size is small. Nevertheless, as the crowd size increases, the relationships between the coordination factor and crowd size for different guidance frequencies tend to coincide, and practical structures are commonly subjected to jumping loads with large crowd. Therefore, the influence of the guidance frequency is overlooked in the modeling of the coordination factor. In addition, the structural damping ratio and non-exceedance probability are also overlooked in the coordination factor model, because these two variables are not related to crowd synchronization. Then, by fitting the average curve of different guidance frequencies in Figure 6, power functions, as depicted in Figure 7, are adopted to establish the model of coordination factor: co 1 n p = 0.01993n p −1.033 + 0.01973 (22) co 2 n p = 0.1012n p −0.6797 + 0.01308 (23) Buildings 2022, 12, x FOR PEER REVIEW (a) (b) Figure 6. Relationship between coordination factor and crowd size.
It is observed from Figure 6 that the influence of the crowd size on the coord factor is affected by the guidance frequency when the crowd size is small. Never as the crowd size increases, the relationships between the coordination factor and size for different guidance frequencies tend to coincide, and practical structu commonly subjected to jumping loads with large crowd. Therefore, the influenc

Conversion Model for Non-Exceedance Probability
Due to the randomness of crowd jumping loads, the equivalent DLF is a random variable. Therefore, the probability distribution of the equivalent DLF is the key to building the conversion model for non-exceedance probability. Some histograms with crowd size equal to 50 and structural damping ratio equal to 0.005 are provided in Figures 8  and 9. It is observed that the equivalent DLF for both the first harmonic and second harmonic follows a normal distribution, and similar histograms are found for other crowd sizes and structural damping ratios.  The equivalent DLF rj of a certain non-exceedance probability p can be deduced by: where F −1 ( ) is the inverse function of the cumulative distribution function for the standard normal distribution, j r is the average of the equivalent DLF, and cvj is the coefficient of variation of the equivalent DLF. According to Equations (17)- (20) and Equations (24) The equivalent DLF r j of a certain non-exceedance probability p can be deduced by: where F −1 ( ) is the inverse function of the cumulative distribution function for the standard normal distribution, r j is the average of the equivalent DLF, and c vj is the coefficient of variation of the equivalent DLF. According to Equations (17)- (20) and Equations (24) and (25), the average of the equivalent DLF can be calculated by: Introducing Equation (27) into Equation 26) and reorganizing the equation yields: Therefore, the coefficient of variation of the equivalent DLF is the key to building the conversion model for the non-exceedance probability. First, the simulations with crowd size equal to 50 and structural damping ratio equal to 0.005 are taken to investigate the c vj of different guidance frequencies. To further determine the number of Monte Carlo simulations, the value of c vj of a different number of Monte Carlo simulations is calculated, as shown in Tables 4 and 5. It can be observed that when the number of Monte Carlo simulations exceeds 600, c vj tends to be stabilized. Therefore, 1000 is adopted as the number of Monte Carlo simulation in Section 3.2, and the simulated results are depicted in Figure 10. It is observed from Figure 10a that the trend of the coefficient of variation for the first harmonic is inverse to the trend in Figure 3a, because the randomness of the crowd jumping load is weak when the crowd synchronization is good. A fourth-order polynomial is then used to fit the coefficient of variation for the first harmonic of different guidance frequencies: As for the second harmonic in Figure 10b, the difference of the coefficient of variation for different guidance frequency is small, so the average of the coefficient of variation for different guidance frequency is adopted: Simulations with a structural damping ratio equal to 0.005 are taken to investigate the relationship between the coefficient of variation and crowd size, and the relationships of some guidance frequencies are illustrated in Figure 11. It is observed from Figure 11 that the coefficient of variation tends to be a constant as the crowd size increases. Because practical structures are commonly subjected to large crowd jumping loads, the influence of the crowd size on the coefficient of variation is overlooked. Then, the simulations with a crowd size equal to 50 are taken to investigate the relationship between coefficient of variation and structural damping ratio, and the relationships of some guidance frequencies are illustrated in Figure 12. It is observed from Figure 12 that the coefficient of variation is only slightly affected by the structural damping ratio, so the influence of the structural damping ratio on the coefficient of variation is ignored.   Simulations with a structural damping ratio equal to 0.005 are taken to investigate the relationship between the coefficient of variation and crowd size, and the relationships of some guidance frequencies are illustrated in Figure 11. It is observed from Figure 11 that the coefficient of variation tends to be a constant as the crowd size increases. Because practical structures are commonly subjected to large crowd jumping loads, the influence of the crowd size on the coefficient of variation is overlooked. Then, the simulations with a crowd size equal to 50 are taken to investigate the relationship between coefficient of variation and structural damping ratio, and the relationships of some guidance frequencies are illustrated in Figure 12. It is observed from Figure 12 that the coefficient of variation is only slightly affected by the structural damping ratio, so the influence of the structural damping ratio on the coefficient of variation is ignored. Simulations with a structural damping ratio equal to 0.005 are taken to investigate the relationship between the coefficient of variation and crowd size, and the relationships of some guidance frequencies are illustrated in Figure 11. It is observed from Figure 11 that the coefficient of variation tends to be a constant as the crowd size increases. Because practical structures are commonly subjected to large crowd jumping loads, the influence of the crowd size on the coefficient of variation is overlooked. Then, the simulations with a crowd size equal to 50 are taken to investigate the relationship between coefficient of variation and structural damping ratio, and the relationships of some guidance frequencies are illustrated in Figure 12. It is observed from Figure 12 that the coefficient of variation is only slightly affected by the structural damping ratio, so the influence of the structural damping ratio on the coefficient of variation is ignored.

Model of the Equivalent DLF
To date, the model of the equivalent DLF for the standard working condition (Equations (17) and (18)) and the conversion model for structural damping ratio (Equations (19) and (20)), crowd size (Equations (24) and (25)), and non-exceedance probability (Equations (28)-(30)) have been established. Consequently, the equivalent DLF can be calculated by:

Model of the Equivalent DLF
To date, the model of the equivalent DLF for the standard working condition (Equations (17) and (18)) and the conversion model for structural damping ratio (Equations (19) and (20)), crowd size (Equations (24) and (25)), and non-exceedance probability (Equations (28)-(30)) have been established. Consequently, the equivalent DLF can be calculated by:

Engineering Application of the Equivalent DLF
The obtained equivalent DLF r j can be used to estimate structural responses under crowd jumping activities in the design stage or as-built stage. First, structural dynamic properties can be obtained by dynamic tests or finite element models. Normally, structural responses are controlled by the resonant mode, where it is assumed the sth vibration mode is resonant, i.e., f s = jf gu . The crowd size, n p , can be determined by structural functional requirements, and the non-exceedance probability, p, can be determined by structural reliability requirements. Then, introducing the values of j, f gu , n p , and p and the damping ratio ζ s , r j can be calculated by Equations (31) and (32). The root mean square value of structural acceleration a rms can then be calculated by: where M s is the modal mass of the resonant mode, normalized by the maximum modal coordinate.

Discussions
The principal contribution of this paper is to develop an equivalent DLF model for crowd jumping loads. The FRF is applied to weight and integrate the PSD matrix to calculate the equivalent DLF. The structural damping ratio is introduced as an independent variable in the model, so that the root mean square of the structural responses can be accurately calculated, which overcomes the shortcomings of the existing DLF [13,[15][16][17][18][19] overestimating or underestimating structural responses. Moreover, the Fourier series model based on the developed equivalent DLF is more convenient and efficient to calculate structural responses than complicated stochastic models [12,[25][26][27][28][29]. In addition, by introducing the non-exceedance probability as an independent variable, a normal distribution-based model is built to obtain the equivalent DLF of different reliabilities, which could avoid using massive Monte Carlo simulations for structural reliability assessment and design. However, to better predict structural responses in real scenes, some limitations and extensions can be described as follows: (1) Practical engineering structures, such as concert halls and grandstands, are commonly occupied by crowds of a large size. A crowd that is in constant contact with the structure, such as a standing crowd or a sitting crowd, can influence structural dynamic properties [34][35][36], causing an increase of the structural damping ratio in particular. Therefore, when Equations (31) and (32) are adopted to obtain the equivalent DLF, the influence of crowd on the structural damping ratio should be considered.
(2) Crowd jumping loads are multi-point excitations, because different people are at different positions in the structure. The Fourier series model based on the developed equivalent DLF simplifies multi-point excitation to single-point excitation, which cannot consider the modal value of each person. As a result, the modal value of each jumping person is assumed to be 1 to predict structural responses, as shown in Equation (33), which indicates that every jumping person is located at the maximum vibration mode, and it can overestimate structural responses. To reduce the errors in structural response calculation, a representative value to characterize each jumping person's modal value needs to be further investigated.

Conclusions
This study builds on an equivalent DLF model for crowd jumping loads based on the massive simulated PSD matrices using a random field model. The equivalent DLF is identified from the principle of equaling the root mean square of structural responses. Four independent variables, the structural damping ratio, guidance frequency, crowd size, and non-exceedance probability, are introduced in the equivalent DLF model. A fourth-order polynomial is used to fit the equivalent DLF of different guidance frequencies. Power functions are used to convert the equivalent DLF of different damping ratios and crowd sizes. A normal distribution-based model is used to convert the equivalent DLF of different non-exceedance probabilities.