Sensitivity Analysis for Pedestrian-Induced Vibration in Footbridges

: This paper aims to provide a novel insight into the inﬂuence of uncertainties in system-and pedestrian-induced load parameters on the vibration response of footbridges. The study begins with a sensitivity analysis for the vertical vibration response of a representative footbridge to two loading cases: a single pedestrian and a crowd. Two methods are utilized: the Sobol’-based global sensitivity analysis method and the local sensitivity analysis method. Uncertainties in all model parameters (which include bridge and human body dynamics in a walking posture, as well as dynamic force generated by humans) are considered in stochastic response estimation. Parametric analysis is then performed to investigate the inﬂuence of the variation of the mean values of the bridge modal mass, damping ratio, and natural frequency on the results of global and local sensitivity analysis. Systematic comparison of the results of global and local sensitivity analysis is performed to identify their similarities and differences. It has been found that the sensitive parameters and their importance ranking strongly depend on bridge modal properties and loading scenarios (i


Introduction
The application of the lightweight high-strength construction materials enabled the design of slender and aesthetically remarkable footbridges [1,2].Such structures are prone to excessive vibration under the excitation of walking pedestrian(s) [3,4].Therefore, vibration serviceability assessment plays a key role in the design of modern footbridges [5].
For the vibration serviceability assessment of a slender footbridge with well-separated vibration modes, design standards routinely require the estimation of the vertical resonant response of the footbridge induced by an average pedestrian traversing the bridge.In such an analysis, the bridge is modeled, in the modal domain, as a single degree of freedom (SDoF) system representing a relevant mode.The walking force is assumed to be a harmonic force having the frequency that matches the natural frequency of the relevant vibration mode to create resonance.More advanced design guidelines also require estimation of the vibration response induced by a crowd [6,7].In this case, the crowd-induced response is calculated by multiplying the resonant response to an average pedestrian by a factor that is a function of the number of pedestrians present on the bridge at any one time and, in some cases, the damping ratio of the bridge.The human-structure interaction (HSI) is usually excluded from the analysis, apart from a basic consideration that pedestrians contribute to the modal mass of the structure.
The dynamic properties and load parameters of a human-structure system are assumed to be deterministic in the analysis above.This leads to a deterministic estimation of vibration response.However, there are inherent uncertainties in load parameters that are due to the inter-and intra-subject variabilities of the pedestrian action.In addition, there unavoidably exist uncertainties in dynamic properties of the structure owing to the inherent uncertainties in the construction process of the structure.As a consequence, the resultant vibration response is difficult to predict accurately.Recently, probability-based human-induced vibration estimation in footbridges by taking into account uncertainties in load parameters has drawn increasing attention [8][9][10][11][12], as it has become clear that stochastic, rather than deterministic, assessment of human-induced vibration response is more suitable for vibration serviceability assessment.
In the last four decades, great efforts have been devoted to experimental characterization of the uncertainties in walking-force parameters [1,[13][14][15][16][17][18][19][20][21][22][23][24][25].The mean value and coefficient of variation (COV) for step frequency, step length, and pedestrian weight obtained by various researchers are summarized in Table 1.The dynamic load factor was found to correlate with step frequency in the references [13,26], while in the study [5] the dynamic load factor was assumed as an independent variable following the normal distribution with a mean value of 0.35 and a COV of 20%.The normal distribution parameters expressed using the two parameters will be presented in short form hereafter, in this case being N (0.35, 20%).
Given that there are inherent uncertainties in the construction process of the structure, the natural frequencies, modal mass (corresponding to modal shapes normalized to the largest component equal to unity), and damping ratios of a bridge are also random parameters.The COVs of the fundamental natural frequency and modal mass are usually no more than 10% [5], while the COV for the damping ratio could be relatively large, ranging from 4% to 25% [43].Note that experiments for stochastic characterization of uncertain parameters in a human-structure system usually require a significant number of test subjects, a large number of trials, and expensive test facilities, making them time consuming and expensive.
Recently, increasing attention has been paid to the sensitivity of human-induced vibration response of a footbridge to uncertainties in the system and load parameters.Sensitivity analysis aims to study how variations in model outputs may be attributed to model inputs.It informs selection of influential parameters required for analyzing uncertainty propagation in stochastic response estimation, as well as for experimental characterization of randomness.This analysis therefore results in significant savings in computational effort in stochastic response estimation, as well as in reducing the cost of experimental characterization.Wei et al. [44] used a polynomial chaos-expansion-based global sensitivity analysis (GSA) method to determine the influence of uncertainties in the mechanical properties of fiber-reinforced polymer (FRP) components on the vibration response of an FRP footbridge exposed to a single pedestrian crossing.They found that the influence of the uncertainties could be neglected on the investigated footbridge.Van Nimmen et al. [45] focused on the influence of variations in a footbridge's natural frequency on the estimated response caused by a crowd crossing.The vibration response was estimated according to design guides SE-TRA [6] and HIVOSS [46], in which the walking force is treated as a deterministic, uniformly distributed load.It was found that for bridges with a frequency of no more than the second harmonic of the walking force, a COV of 10% for natural frequency may lead to a big scatter of estimated peak acceleration responses, e.g., ranging from 1.4 m/s 2 to 5.2 m/s 2 .Caprani et al. [47] evaluated the stochastic response of a virtual footbridge with dynamic properties treated as deterministic parameters while a single pedestrian loading was modeled as a stochastic force model.The uncertainties in step frequency, step length, and weight of pedestrian were considered.They found the estimated footbridge acceleration response is most sensitive to the step frequency.Similarly, Pedersen and Frier [8] studied the sensitivity of footbridge response to step frequency, walking speed, pedestrian weight, and dynamic load factor characterizing a single pedestrian and concluded that the step frequency is the most influential parameter.Shahabpoor et al. [29] examined the impact of the mass, natural frequency, and damping ratio of the human body and a pedestrian's walking speed, as well as the pedestrians' arrival rate on the multi-pedestrian-induced footbridge response experienced by walking pedestrians rather than at a fixed location.The HSI was considered by modeling each pedestrian as an SDoF system.The results showed that the footbridge acceleration response is most sensitive to a human body's frequency in cases when this frequency is close to the bridge frequency.The sensitivity of the response to the other four parameters was significantly lower.More specifically, a 30% variation in one of these four parameters changed the response up to only 10%, whereas a 10% variation in a human body's frequency can achieve the same effect.
The existing literature provides an incomplete picture of the sensitive parameters, i.e., parameters that strongly influence footbridge vibration response.This is mainly due to the fact that only partial uncertainties were examined, i.e., either uncertainties in the parameters of the structure, or human-induced dynamic force, or human body dynamics.Another possible reason is the choice of the local sensitivity analysis (LSA) in the existing research.LSA is a one-at-a-time technique, evaluating variations in model outputs with respect to one input, with the remaining parameters fixed [48].For an explicit input-output mapping relationship, the local sensitivity of the output may be defined as the partial derivative of the output with respect to one input of interest.A local sensitivity index is only valid in the vicinity of the base point where it is evaluated and may vary with the location of the base point [49].Hence, LSA provides only a limited insight into model sensitivity and may lead to inconsistent sensitivity analysis results.A third possible reason is that the results from sensitivity analysis may vary with the input spaces of uncertain parameters.In existing literature, different stochastic distributions or intervals were employed for the same uncertain parameters by different researchers.
The contribution of this paper is to provide a more conclusive insight into the sensitivity of a pedestrian-induced vibration response.More specifically, stochastic vibration responses of a representative footbridge induced by a single pedestrian or a crowd walking are first estimated.Uncertainties in the dynamic properties of the structure and a human body, as well as parameters of human-induced load, are all included, which belong to aleatory uncertainty.The epistemic uncertainty, such as the uncertainty in distribution types of parameters, are out of the scope of the paper.The HSI is taken into account by modeling the human body as an SDoF system in conjunction with a vertical ground reaction force crossing the bridge at a constant speed.The crowd is assumed to constitute spatially unrestricted pedestrian traffic (characterized by a relatively low traffic density in which each person is free to choose a walking speed independently from the other pedestrians).GSA is then carried out using the Sobol' method to investigate the influence of uncertainties in the system and load parameters.In this analysis, all parameters are varied simultaneously over the specified input space.As a result, the contributions of each individual parameter and the interactions between parameters to model output are evaluated simultaneously.GSA provides more reliable results than LSA by characterizing the response sensitivity across the entire input space.Then, parametric analysis is conducted to study the influence of the input spaces of bridge modal mass, damping ratio, and natural frequency on the results of GSA.In addition, LSAs are also performed in all cases so that the similarities and differences between GSA and LSA can be observed and evaluated.
The layout of this article is as follows.This introductory section is followed by Section 2, which briefly introduces the Sobol' method.Section 3 presents the estimation methods for the vertical vibration response of a footbridge induced by a single pedestrian and a crowd.In Section 4, sensitivity analysis is carried out to reveal the global and local sensitivity of the acceleration response of a representative footbridge to the two load cases.Parametric analysis is performed in Section 5 to investigate how the variations in the input spaces of modal mass, damping ratio, and natural frequency of the footbridge affect the results of GSA and LSA.The discussion is presented in Section 6.The concluding remarks are given in Section 7.

The Sobol' Method
This section presents a brief introduction of the Sobol' method.A detailed description of the method can be found elsewhere [50].

Formulation of the Sobol' Method
The model under investigation may be described by the following function: where Y is a scalar model output and x = (x 1 , x 2 , . . . ,x k ) is a vector of k model inputs.
It is assumed that f (x) is square integrable.In addition, the input parameters are assumed to be independent, random, and within the unit hypercube (i.e., x ∈ I k and I is the unit interval [0, 1]).This incurs no loss of generality because any space of a random variable can be transformed onto this unit hypercube.
The right-hand side of Equation ( 1) can be decomposed into the form of ANOVArepresentation [51]: It follows that the summands in Equation ( 2) are orthogonal and can be expressed as the integrals of f (x).The physical meanings of the summands are: f 0 represents the average of the model f (x), f i (x i ) describes the main effect of input x i , f ij x i , x j represents the effect of the interaction between x i and x j , and so on.
Squaring Equation (2) and integrating over I k , we can deduce that the variance of model output can be decomposed into fractions which are contributed by the sets of inputs: where D is the variance of model output.D i is the variance corresponding to input x i , D ij is the variance contributed by both inputs x i and x j , and so on.The Sobol' indices can be defined as: where S i is the first-order sensitivity index (FSI) that measures the contribution of x i alone to the output variance D. S ij is the second-order sensitivity index that measures the contribution of the terms involving the interaction between x i and x j to the total variance D, and so on.Each sensitivity index is in the range from 0 to 1.The sum of all of sensitivity indices is equal to 1.
In addition, a total sensitivity index (TSI), measuring the contribution of the input x i and all the possible joint terms between x i and all the other inputs, may be given by: where S T i is TSI for input x i .D −i represents the total variance uncorrelated to input x i .Usually, the larger the TSI of an input, the more influence it has on the model output.In addition, the difference between the FSI and TSI of a parameter is accounted by its interaction with other parameters.
Let y = x l 1 , . . ., x l m , 1 ≤ l 1 < • • • < l m ≤ k, be an arbitrary subset of x and z is the complementary set.Thus, x = (y, z).Let L = (l 1 , • • • , l m ).The variance corresponding to input y can be defined as [50]: Similarly, D z can be defined.The FSI S y and TSI S T y of subset y then can be expressed as: Equations ( 7) and (8) become Equation (4) (the first equation) and Equation ( 5), respectively, when y = x i .

Implementation of the Sobol' Method
The variances D, D y , and D z can be calculated using a Monte Carlo algorithm as [50]: Buildings 2022, 12, 883 6 of 21 where N is the number of Monte Carlo trails, which should be large enough to ensure convergence.To perform Equations ( 10) and (11), y (z) needs to be sampled twice for each simulation, such that y d and y d (z d and z d ) are independent from each other.y d represents the sample of y in the first sampling for the dth simulation and y d represents the sample of y in the second sampling for the dth simulation.Similarly, z d and z d can be defined.

Pedestrian-Induced Vibration Analysis for Footbridges
This section describes the approaches for estimating the vibration response of a footbridge induced by a single pedestrian and a crowd.The footbridge is modeled as an SDoF system associated with a relevant vibration mode, and the influence of a pedestrian walking on a footbridge is modeled as another SDoF system in conjunction with a ground reaction force crossing the bridge at a constant speed.

Estimation of the Vibration Response Induced by a Single Pedestrian
The equations of motion of a footbridge with a single pedestrian crossing may be written as [4]: where: m b , ξ b , and f b are the modal mass (kg), damping ratio (%), and natural frequency (Hz) of the footbridge, respectively, while m p , ξ p , and f p are the mass (kg), damping ratio (%) and natural frequency (Hz) of the pedestrian, respectively.x b (t) and x p (t) are the modal displacements of the footbridge and pedestrian, respectively.F(t) = Wα cos(2π f s t) is the walking force.W, α, and f s are the weight (N), dynamic load factor, and step frequency (Hz) of the pedestrian, respectively.Φ(vt) is the mode shape of the bridge and v is the walking speed, which is the product of step frequency f s and step length l s .( .. ) and ( . ) are the second-order and first-order derivatives of () with respect to time t.

Estimation of the Vibration Response Induced by a Crowd
Consider a crowd of r pedestrians crossing a footbridge in one direction.Pedestrian arrival is assumed to be a Poisson process.This means that the time lag (τ) between pedestrians follows an exponential distribution [52].The equations of motion of a footbridge with a crowd crossing are [4]: where: The submatrices in Equation ( 13) are defined as follows: Buildings 2022, 12, 883 7 of 21 are the mass, damping ratio, natural frequency, walking speed, modal displacement, weight, dynamic load factor, step frequency, and phase angle for the n-th pedestrian, respectively.The arrival time of the n-th pedestrian on the bridge can be expressed as t on,n = n ∑ i=1 τ i .The departure time of the n-th pedestrian off the bridge can be expressed as t o f f ,n = t on,n + L/( f s,n l s,n ).L, l s,n and τ i are the bridge length, step length of the n-th pedestrian, and time lag between the i-th pedestrian and the (i − 1)-th pedestrian, respectively.The definitions of the other variables have been given in Section 3.1.

Sensitivity Analysis of the Vibration Response of a Representative Footbridge
This section presents the GSA and LSA of the vibration response of a representative footbridge induced by a single pedestrian and a crowd.The representative footbridge is an idealized bridge, having a span of 50 m and a deck width of 2 m, which is typically employed in numerical examples in the literature related to human-induced vibration [4,5].

Representative Values for System and Load Parameters
The first vertical mode of the bridge is considered only.Thus, the mode shape is a halfsine function.According to the database of measured dynamic properties of 138 footbridges summarized by Wei et al. [53], the mean value for the modal mass, damping ratio, and fundamental frequency of the representative footbridge can be assumed to take values of 40,000 kg (based on a representative mass per deck area of 800 kg/m 2 ), 1% and 2 Hz, respectively.The COVs for the modal mass and fundamental frequency of the footbridge are both assumed to be 10%, while the COV for the damping ratio is doubled because of the usually larger uncertainty in damping identification.
The pedestrian-generated dynamic load from Section 3 is a harmonic function characterized by force amplitude (α times W), step frequency f s , and step length l s in the case of a single pedestrian.There are two additional parameters, phase angle ϕ and time-lag τ, in the case of a crowd.The dynamic model for a walking pedestrian is characterized by the mass m p (W divided by gravitational acceleration), natural frequency f p , and damping ratio ξ p .The mean and COV values for step frequency, step length, weight, damping ratio, and natural frequency of a single pedestrian are assigned the average measured values reported in the literature and summarized in Section 1.They are listed in Table 2.The pedestrian density is set to be 0.3 pedestrians/m 2 of deck area, which corresponds to free (i.e., spatially unrestricted) walking [54].All the parameters characterizing dynamic properties and dynamic load are assumed to follow the normal distribution, except the phase angle and time lag between pedestrians.The phase angle is assumed to obey the uniform distribution between 0 and 2 π [4].The time lag between pedestrians follows an exponential distribution with a mean value of 1.2 s, which is the reciprocal of mean arrival rate of 0.84 pedestrians/s.Such an arrival rate corresponds to an average pedestrian density of 0.3 pedestrians/m 2 for a crowd walking, on average, at a walking speed of 1.4 m/s (which is equal to the mean step frequency of 2 Hz multiplied by the mean step length of 0.7 m).In addition, the mean value and COV of the dynamic load factor are in agreement with those provided by Tubino et al. [5], in which the dynamic load factor was assumed to independent of step frequency.All parameters in Table 2 are assumed to be mutually independent random variables.Since the fundamental frequency of the bridge is in the range of the step frequency of a normal-walking pedestrian, only the first-order harmonic of the walking force is considered when estimating human-induced vibration.In the case of a single pedestrian, the input parameters for sensitivity analysis are the step frequency f s , step length l s , dynamic load factor α, weight W, damping ratio ξ p , and natural frequency f p of the pedestrian, and the fundamental frequency f b , damping ratio ξ b , and modal mass m b of the footbridge.For the case of a crowd, there are two additional inputs: phase angle ϕ and time lag τ.In both cases, the model output for sensitivity analysis is the maximum of one-second root-mean-squared (RMS) acceleration at the mid-span of the footbridge.Based on the input and output parameters mentioned above, the sensitivity analysis of the human-induced footbridge response can be carried out using the Sobol' method and the local sensitivity analysis method.

Vibration Response Induced by a Single Pedestrian
The vibration response of the footbridge induced by a single pedestrian is estimated by solving Equation ( 12) using the Newmark-β numerical integration method [55].Due to the step frequency and step length of pedestrian are random variables, the time to cross a 50 m long bridge is different for each pedestrian.Therefore, the time duration for each simulation is different, which is set to be 110% of the time for crossing the bridge to ensure the time history experiencing peak response and then decaying to a low response level at the end [4].Monte Carlo simulations based on low-discrepancy sequences are conducted.It is found that the results of the sensitivity analysis are almost unchanged with a simulation size larger than 5000.Therefore, 10,000 Monte Carlo simulations are performed for reliable results.

Global Sensitivity Analysis
In this paper, the aim of the sensitivity analysis is to identify and rank the parameters that influence the vibration response most.Both FSI and TSI from GSA could be employed for distinguishing between sensitive and insensitive parameters, but a literature search does not provide consistent thresholds for either.For example, Tang et al. [56], Coppitters et al. [57], and Hsieh et al. [58] considered parameters of an FSI below 0.01, 0.02, and 0.05, respectively, to be insensitive.By contrast, Zhang et al. [59] and Chan et al. [60] employed a TSI with threshold values of 0.05 and 0.3, respectively.Although FSIs are frequently utilized for parameter selection and ranking, they are insufficient to characterize parameter importance when there are strong parameter interactions [50].Hence, the selection and ranking of sensitive parameters in this paper are decided using a TSI threshold value of 0.07.However, both FSI and TSI will be calculated and presented.
The FSIs and TSIs of the maximum of one-second RMS acceleration at the mid-span to all the parameters are calculated using the Sobol' method described in Section 2. The 9 of 21 results from GSA are depicted in Figure 1a, where S 1 and S T denote the FSI and TSI, respectively.According to the TSIs, the bridge's natural frequency is the most important parameter, followed closely by step frequency and then dynamic load factor, while the TSI for other parameters is less than 0.07 and is of little importance.The interactions between the uncertain parameters can be detected according to the difference between TSI and FSI, owing to the fact that TSI can measure the contributions of an uncertain parameter, as well as of the interactions between this parameter and all other parameters to the variation of vibration response, but FSI measures the contribution of this parameter only.It can be found that the TSIs for a bridge's natural frequency and step frequency are significantly larger than their first-order counterparts.This indicates a strong interaction between these two parameters, which is accounted for by the high probability of occurrence of resonance with a step frequency close to the bridge's natural frequency.
a TSI with threshold values of 0.05 and 0.3, respectively.Although FSIs are frequently utilized for parameter selection and ranking, they are insufficient to characterize parameter importance when there are strong parameter interactions [50].Hence, the selection and ranking of sensitive parameters in this paper are decided using a TSI threshold value of 0.07.However, both FSI and TSI will be calculated and presented.
The FSIs and TSIs of the maximum of one-second RMS acceleration at the mid-span to all the parameters are calculated using the Sobol' method described in Section 2. The results from GSA are depicted in Figure 1a, where  1 and   denote the FSI and TSI, respectively.According to the TSIs, the bridge's natural frequency is the most important parameter, followed closely by step frequency and then dynamic load factor, while the TSI for other parameters is less than 0.07 and is of little importance.The interactions between the uncertain parameters can be detected according to the difference between TSI and FSI, owing to the fact that TSI can measure the contributions of an uncertain parameter, as well as of the interactions between this parameter and all other parameters to the variation of vibration response, but FSI measures the contribution of this parameter only.It can be found that the TSIs for a bridge's natural frequency and step frequency are significantly larger than their first-order counterparts.This indicates a strong interaction between these two parameters, which is accounted for by the high probability of occurrence of resonance with a step frequency close to the bridge's natural frequency.

Local Sensitivity Analysis
In this section, LSA for vibration response of the representative footbridge under the excitation of a single pedestrian is carried out.Specifically, one input parameter varies at a time according to its statistical characteristics and the other inputs are kept at their nominal values.The corresponding COV for the vibration response of the bridge, estimated using the Monte Carlo simulation method, is then calculated.Note that the LSA can only produce qualitative results [48], i.e., ranking the input parameters.In other words, it is

Local Sensitivity Analysis
In this section, LSA for vibration response of the representative footbridge under the excitation of a single pedestrian is carried out.Specifically, one input parameter varies at a time according to its statistical characteristics and the other inputs are kept at their nominal values.The corresponding COV for the vibration response of the bridge, estimated using the Monte Carlo simulation method, is then calculated.Note that the LSA can only produce qualitative results [48], i.e., ranking the input parameters.In other words, it is impossible to establish a quantitative criterion to distinguish important from unimportant parameters.For the sake of comparing the results of GSA and LSA, a top few parameters ranked by LSA are considered as the sensitive parameters as much as those given by GSA.
Figure 1b presents the results of LSA for vibration response induced by a single pedestrian.By comparing with Figure 1a, it can be found that the result of LSA is similar to the TSIs, rather than the FSIs, from the GSA.This is mainly due to the fact that for an uncertain parameter, both the TSI and the result of LSA simultaneously measure the contributions of this parameter as well as of the interactions between this parameter and all other parameters to the variation of vibration response, while the FSI only measure the contribution of this single parameter.The results of LSA show that the bridge natural frequency is the most important parameter, followed by step frequency and then dynamic load factor.This agrees with the conclusion from GSA.However, LSA cannot identify the interactions between the uncertain parameters as the same as the GSA.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 4.2 is in the range of 0.16 m/s 2 to 0.23 m/s 2 .

Vibration Response Induced by a Crowd
The vibration response of the footbridge induced by a crowd is computed using the same numerical integration method to solve Equation (13).To ensure a long-enough duration of steady-state pedestrian flow, i.e., the period after the departure time of the first pedestrian off the bridge and before the arrival time of the last pedestrian on the bridge, the number of people in a crowd is set to be 150.In this case, the maximum of one-second RMS acceleration exhibits no obvious variation during this period.Due to the randomness in the pedestrian step frequency, step length and time-lag between pedestrians, the time required for a crowd to cross the bridge is different each time.Therefore, the time duration for each simulation is set to be the time for a crowd to cross the bridge completely.Again, 10,000 Monte Carlo simulations are conducted.

Global Sensitivity Analysis
The FSIs and TSIs for the system and load parameters are shown in Figure 2a.The sensitive parameters, in descending order (of TSI), are the bridge natural frequency, step frequency, phase angle, time lag, step length, and bridge modal mass.Hereafter, the sensitive parameters are listed in descending order by default.Compared with the case of a single pedestrian, the dynamic load factor loses on significance.In addition, there are big differences between the TSIs and FSIs for the bridge's natural frequency, step frequency, phase angle, and time lag because of the interactions between these four parameters.It should be noted that for a system with nonlinearity like the one in this paper, TSI is recommended to measure the sensitivity of uncertain parameters instead of FSI.This is because the FSI cannot evaluate interactions from nonlinear effects, which may lead to the misleading result.
The vibration response of the footbridge induced by a crowd is computed using the same numerical integration method to solve Equation (13).To ensure a long-enough duration of steady-state pedestrian flow, i.e., the period after the departure time of the first pedestrian off the bridge and before the arrival time of the last pedestrian on the bridge, the number of people in a crowd is set to be 150.In this case, the maximum of one-second RMS acceleration exhibits no obvious variation during this period.Due to the randomness in the pedestrian step frequency, step length and time-lag between pedestrians, the time required for a crowd to cross the bridge is different each time.Therefore, the time duration for each simulation is set to be the time for a crowd to cross the bridge completely.Again, 10,000 Monte Carlo simulations are conducted.

Global Sensitivity Analysis
The FSIs and TSIs for the system and load parameters are shown in Figure 2a.The sensitive parameters, in descending order (of TSI), are the bridge natural frequency, step frequency, phase angle, time lag, step length, and bridge modal mass.Hereafter, the sensitive parameters are listed in descending order by default.Compared with the case of a single pedestrian, the dynamic load factor loses on significance.In addition, there are big differences between the TSIs and FSIs for the bridge's natural frequency, step frequency, phase angle, and time lag because of the interactions between these four parameters.It should be noted that for a system with nonlinearity like the one in this paper, TSI is recommended to measure the sensitivity of uncertain parameters instead of FSI.This is because the FSI cannot evaluate interactions from nonlinear effects, which may lead to the misleading result.

Local Sensitivity Analysis
The results from LSA are shown in Figure 2b.The most influential parameter is the bridge's natural frequency, followed by phase-angle and then step frequency.The bridge damping ratio, time lag, and bridge modal mass are the least sensitive parameters.The other parameters are of little influence.Comparison of the results in Figure 2a,b shows that LSA and GSA differ in ranking the parameters.This inconsistency may be attributed to the fact that a local sensitivity index is estimated only in the vicinity of the nominal point, leading to the incomplete investigation of the full input space, while GSA produces an estimation of model sensitivity over the full input space including the nominal point.As a consequence, GSA, rather than LSA, is recommended as a choice for the sensitivity analysis of the pedestrian-induced vibration of footbridges.This is owing to the fact that the results of LSA may not include the effects of interactions between uncertain parameters in the unexplored input space, while GSA does not have this problem and therefore leads to more reliable results.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 4.3 is in the range of 0.61 m/s 2 to 2.64 m/s 2 .

Parametric Analysis of the Input Spaces of Modal Mass, Damping Ratio, and Natural Frequency
The GSA and LSA presented in Section 4 are for a particular footbridge.To provide more robust insight, the sensitivity analysis for different footbridges is performed in this section.It includes both GSA and LSA for the vibration response of footbridges with different mean values of modal mass, damping ratio, and natural frequency.

Natural Frequency
The footbridges with a natural frequency up to 3 Hz may be strongly excited by the first-order harmonic of a pedestrian's walking force.Thus, six footbridges with different mean natural frequencies from 0.5 Hz to 3 Hz (with a step size of 0.5 Hz) are considered in the parametric analysis.The natural frequency of these six footbridges is described by normal distribution: N (0.5-3 HZ, 10%).The statistical characteristics for the other parameters are taken from Table 2.These six footbridges are referred to as Bridge Series A.

Vibration Response Induced by a Single Pedestrian Global Sensitivity Analysis
The generated TSIs for Bridge Series A are presented in Figure 3a, leading to the following observations:

•
For the bridges with the mean natural frequency of 0.5 Hz, the vibration responses are most sensitive to dynamic load factor, followed by pedestrian weight and then bridge modal mass.For the bridges with the mean natural frequency of 1 Hz, there are two additional sensitive parameters step frequency and bridge natural frequency.

•
For the bridges with the mean natural frequency no smaller than 1.5 Hz, the variations in step frequency and bridge natural frequency have most significant impact on the vibration response.The dynamic load factor becomes the third sensitive parameter for the bridges with the mean natural frequency of 2 Hz.The results of LSAs for Bridge Series A are presented in Figure 3b.The sensitivity of parameters and their ranking agree with those from GSA.
The results from GSA and LSA suggest that the parameters that influence vibration response most strongly depend on the natural frequency of the bridge, or more precisely, its closeness to the step frequency and therefore likelihood of exciting the resonance.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.

Global Sensitivity Analysis
The TSIs from the GSAs for Bridge Series A are presented in Figure 4a, which suggests that:

•
For the bridges with mean natural frequency of 0.5 Hz, the step frequency and phase angle are the first two contributors to the response's variation, followed by the time lag and then the bridge modal mass.The step length and dynamic load factor are the

Local Sensitivity Analysis
The results of LSAs for Bridge Series A are presented in Figure 3b.The sensitivity of parameters and their ranking agree with those from GSA.
The results from GSA and LSA suggest that the parameters that influence vibration response most strongly depend on the natural frequency of the bridge, or more precisely, its closeness to the step frequency and therefore likelihood of exciting the resonance.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.1.1 is in the range of 0.01 m/s 2 to 0.23 m/s 2 .

Vibration Response Induced by a Crowd Global Sensitivity Analysis
The TSIs from the GSAs for Bridge Series A are presented in Figure 4a, which suggests that:

•
For the bridges with mean natural frequency of 0.5 Hz, the step frequency and phase angle are the first two contributors to the response's variation, followed by the time lag and then the bridge modal mass.The step length and dynamic load factor are the remaining two sensitive parameters.

•
For the bridges with the mean natural frequency of 1 Hz, an additional sensitive parameter is the bridge natural frequency.

•
For the bridges with the mean natural frequencies of 1.5 Hz and 2.5 Hz, the most sensitive parameter is the bridge frequency, followed by step frequency and then phase angle.For the bridges with the mean natural frequency at 2 Hz, there are three additional sensitive parameters: time lag, step length, and bridge modal mass.For the bridges with the mean natural frequency of 3 Hz, only the first two parameters are sensitive.
Buildings 2022, 12, x FOR PEER REVIEW 13 of 22 • For the bridges with the mean natural frequency of 0.5 Hz, there are six sensitive parameters, i.e., the phase angle, step frequency, time lag, bridge modal mass, dynamic load factor, and step length.Apart from these six parameters, an additional parameter, the bridge natural frequency, becomes sensitive for the bridges with the mean natural frequency of 1 Hz.

•
For the bridges with the mean natural frequency of 1.5 Hz and 2.5 Hz, the most sensitive parameter is the bridge natural frequency, followed by step frequency and then phase angle.When the bridges have the mean natural frequency of 3 Hz, the step frequency becomes insensitive.

•
For the bridges with the mean natural frequency of 2 Hz, the variation in bridge natural frequency impacts most the vibration response, followed by phase angle and then step frequency.The bridge damping ratio, time lag, and bridge modal mass are the last three sensitive parameters.By comparing Figure 4a,b, it is found that GSA and LSA generate the same sensitive parameters and their ranking for the bridges with the mean natural frequency of 1.5 Hz and 2.5 Hz, different sensitive parameters when the mean natural frequency is taken to be 2 Hz and 3 Hz, and the same parameters but different ranking when the mean natural frequency no larger than 1 Hz.Again, both the results from GSA and LSA vary with the input space of bridge natural frequency.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.1.2is in the range of 0.04 m/s 2 to 2.64 m/s 2 .

Local Sensitivity Analysis
The results of LSAs for Bridge Series A induced by a crowd are presented in Figure 4b.

•
For the bridges with the mean natural frequency of 0.5 Hz, there are six sensitive parameters, i.e., the phase angle, step frequency, time lag, bridge modal mass, dynamic load factor, and step length.Apart from these six parameters, an additional parameter, the bridge natural frequency, becomes sensitive for the bridges with the mean natural frequency of 1 Hz.

•
For the bridges with the mean natural frequency of 1.5 Hz and 2.5 Hz, the most sensitive parameter is the bridge natural frequency, followed by step frequency and then phase angle.When the bridges have the mean natural frequency of 3 Hz, the step frequency becomes insensitive.

•
For the bridges with the mean natural frequency of 2 Hz, the variation in bridge natural frequency impacts most the vibration response, followed by phase angle and then step frequency.The bridge damping ratio, time lag, and bridge modal mass are the last three sensitive parameters.
By comparing Figure 4a,b, it is found that GSA and LSA generate the same sensitive parameters and their ranking for the bridges with the mean natural frequency of 1.5 Hz and 2.5 Hz, different sensitive parameters when the mean natural frequency is taken to be 2 Hz and 3 Hz, and the same parameters but different ranking when the mean natural frequency no larger than 1 Hz.Again, both the results from GSA and LSA vary with the input space of bridge natural frequency.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.1.2is in the range of 0.04 m/s 2 to 2.64 m/s 2 .

Modal Mass
Similarly, parametric analysis is carried out to investigate the robustness of results of sensitivity analysis to the variation of the input space of modal mass.According to [53], the physical mass per square meter of deck area for a footbridge ranges from 52 kg/m 2 to 1955 kg/m 2 .The footbridge with low modal mass is prone to vibration problems.Here, seven footbridges, the mean physical mass density ranging from 50 kg/m 2 to 1600 kg/m 2 , are analyzed.Their mean modal mass is taken to be 2500 kg, 5000 kg, 7500 kg, 10,000 kg, 20,000 kg, 40,000 kg, and 80,000 kg, respectively.The modal mass of these seven footbridges is described by normal distribution: N (5000-80,000 kg, 10%).The statistical characteristics of the remaining parameters are taken from Table 2.These seven footbridges are designated as Bridge Series B.

Vibration Response Induced by a Single Pedestrian Global Sensitivity Analysis
Figure 5a shows the TSIs from the GSAs for Bridge Series B, which indicate that the bridge frequency is always the main contributor, with its TSI slightly larger than that of step frequency.The dynamic load factor is the third sensitive parameter with a TSI of about 0.07.1955 kg/m 2 .The footbridge with low modal mass is prone to vibration problems.Here, seven footbridges, the mean physical mass density ranging from 50 kg/m 2 to 1600 kg/m 2 , are analyzed.Their mean modal mass is taken to be 2500 kg, 5000 kg, 7500 kg, 10,000 kg, 20,000 kg, 40,000 kg, and 80,000 kg, respectively.The modal mass of these seven footbridges is described by normal distribution: N (5000-80,000 kg, 10%).The statistical characteristics of the remaining parameters are taken from Table 2.These seven footbridges are designated as Bridge Series B.

Vibration Response Induced by a Single Pedestrian Global Sensitivity Analysis
Figure 5a shows the TSIs from the GSAs for Bridge Series B, which indicate that the bridge frequency is always the main contributor, with its TSI slightly larger than that of step frequency.The dynamic load factor is the third sensitive parameter with a TSI of about 0.07.

Local Sensitivity Analysis
The results of LSAs are shown in Figure 5b.The influential parameters are the bridge frequency, step frequency, and dynamic load factor, and they are independent from modal mass.
The comparison between Figure 5a,b indicates that LSA and GSA produce consistent results in terms of sensitive parameters and their ranking.In addition, it can be concluded that the variation of the input space of modal mass will not affect the results of GSA and LSA under the excitation of a single pedestrian.

Local Sensitivity Analysis
The results of LSAs are shown in Figure 5b.The influential parameters are the bridge frequency, step frequency, and dynamic load factor, and they are independent from modal mass.
The comparison between Figure 5a,b indicates that LSA and GSA produce consistent results in terms of sensitive parameters and their ranking.In addition, it can be concluded that the variation of the input space of modal mass will not affect the results of GSA and LSA under the excitation of a single pedestrian.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.2.1 is in the range of 0.08 m/s 2 to 2.17 m/s 2 .

Vibration Response Induced by a Crowd Global Sensitivity Analysis
The TSIs from GSAs for Bridge Series B are obtained and presented in Figure 6a.For all seven bridges, there are six influential parameters, i.e., the bridge natural frequency, step frequency, phase angle, time lag, step length, and bridge modal mass.For the bridges with the mean modal mass no larger than 10,000 kg, there is an additional sensitive parameter, the dynamic load factor.The TSIs from GSAs for Bridge Series B are obtained and presented in Figure 6a.For all seven bridges, there are six influential parameters, i.e., the bridge natural frequency, step frequency, phase angle, time lag, step length, and bridge modal mass.For the bridges with the mean modal mass no larger than 10,000 kg, there is an additional sensitive parameter, the dynamic load factor.

Local Sensitivity Analysis
For the scenario of a crowd, the results from LSA are presented in Figure 6b.For the bridges with the mean modal mass no smaller than 20,000 kg, there are six influential parameters, which are the bridge's natural frequency, phase angle, step frequency, bridge damping ratio, time lag, and bridge modal mass.For the bridges with the mean modal mass of 2500 kg, there is an additional sensitive parameter, the step length, and for the bridges with the mean modal mass ranging from 5000 kg to 10,000 kg, the seventh sensitive parameter is the dynamic load factor.Compared with the results of GSA, LSA may lead to different influential parameters.Based the analysis above, it can be concluded that the results of GSA and LSA may vary with the input space of modal mass.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.2.2 is in the range of 0.35 m/s 2 to 8.4 m/s 2 .

Damping Ratio
This section examines the influence of the variation of the input space of the damping ratio on the results of sensitivity analysis.The footbridge with a low damping ratio, usually no more than 2%, is prone to vibration problems.Here, nine footbridges are analyzed.Their mean damping ratios are taken to be 0.2%, 0.3%, 0.4%, 0.5%, 0.75%, 1%, 1.25%, 1.5%, and 2%, respectively.The damping ratio of these nine footbridges is described by normal distribution: N (0.2-2%, 20%).The statistical characteristics of the remaining parameters are taken from Table 2.These nine footbridges are denoted as Bridge Series C.

Local Sensitivity Analysis
For the scenario of a crowd, the results from LSA are presented in Figure 6b.For the bridges with the mean modal mass no smaller than 20,000 kg, there are six influential parameters, which are the bridge's natural frequency, phase angle, step frequency, bridge damping ratio, time lag, and bridge modal mass.For the bridges with the mean modal mass of 2500 kg, there is an additional sensitive parameter, the step length, and for the bridges with the mean modal mass ranging from 5000 kg to 10,000 kg, the seventh sensitive parameter is the dynamic load factor.
Compared with the results of GSA, LSA may lead to different influential parameters.Based the analysis above, it can be concluded that the results of GSA and LSA may vary with the input space of modal mass.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.2.2 is in the range of 0.35 m/s 2 to 8.4 m/s 2 .

Damping Ratio
This section examines the influence of the variation of the input space of the damping ratio on the results of sensitivity analysis.The footbridge with a low damping ratio, usually no more than 2%, is prone to vibration problems.Here, nine footbridges are analyzed.Their mean damping ratios are taken to be 0.2%, 0.3%, 0.4%, 0.5%, 0.75%, 1%, 1.25%, 1.5%, and 2%, respectively.The damping ratio of these nine footbridges is described by normal distribution: N (0.2-2%, 20%).The statistical characteristics of the remaining parameters are taken from Table 2.These nine footbridges are denoted as Bridge Series C.

Vibration Response Induced by a Single Pedestrian Global Sensitivity Analysis
The generated TSIs for Bridge Series C are displayed in Figure 7a.For all nine bridges, there are two sensitive parameters, among which the step frequency is the most important with the TSI slightly larger than that of the bridge's natural frequency.The dynamic load factor becomes the third sensitive parameter for the bridges with the mean damping ratio no smaller than 1%.The generated TSIs for Bridge Series C are displayed in Figure 7a.For all nine bridges, there are two sensitive parameters, among which the step frequency is the most important with the TSI slightly larger than that of the bridge's natural frequency.The dynamic load factor becomes the third sensitive parameter for the bridges with the mean damping ratio no smaller than 1%.

Local Sensitivity Analysis
The LSAs for Bridge Series C in the single pedestrian scenario are presented in Figure 7b.For all nine bridges, the vibration responses are most sensitive to bridge natural frequency and step frequency.For the bridges with the mean damping ratio no smaller than 1%, there is an additional sensitive parameter, the dynamic load factor.The comparison between the results in Figure 7a,b indicates that GSA and LSA produce the same sensitive parameters but different rankings for all nine bridges.It can be found that the results from GSA and LSA may change with the input space of the damping ratio.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.3.1 is in the range of 0.1 m/s 2 to 0.5 m/s 2 .

Global Sensitivity Analysis
The TSIs from GSAs for Bridge Series C are shown in Figure 8a.For all nine bridges, the sensitive parameters are the bridge natural frequency, step frequency, phase angle, time lag, step length, and bridge modal mass.For the bridges with the mean damping ratio of 2%, there is an additional sensitive parameter, the bridge damping ratio.

Local Sensitivity Analysis
The LSAs for Bridge Series C in the single pedestrian scenario are presented in Figure 7b.For all nine bridges, the vibration responses are most sensitive to bridge natural frequency and step frequency.For the bridges with the mean damping ratio no smaller than 1%, there is an additional sensitive parameter, the dynamic load factor.
The comparison between the results in Figure 7a,b indicates that GSA and LSA produce the same sensitive parameters but different rankings for all nine bridges.It can be found that the results from GSA and LSA may change with the input space of the damping ratio.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.3.1 is in the range of 0.1 m/s 2 to 0.5 m/s 2 .

Vibration Response Induced by a Crowd Global Sensitivity Analysis
The TSIs from GSAs for Bridge Series C are shown in Figure 8a.For all nine bridges, the sensitive parameters are the bridge natural frequency, step frequency, phase angle, time lag, step length, and bridge modal mass.For the bridges with the mean damping ratio of 2%, there is an additional sensitive parameter, the bridge damping ratio.
bridge damping ratio, time lag, and dynamic load factor.

•
For the bridges with the mean damping ratio ranging from 0.4% to 1.5%, the first five sensitive parameters are the same as those for the bridges with the mean value of damping ratio no bigger than 0.3%, as well as the ranking, but the sixth sensitive parameter is replaced by the bridge modal mass.Apart from these six parameters, an additional parameter, the dynamic load factor, becomes sensitive for the bridges with the mean damping ratio of 2%.
In summary, the results of GSA and LSA sensitive to the input space of damping ratio.Compared with GSA, LSA produces the different sensitive parameters for each damping ratio case.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.

Discussion
The results of GSA and LSA for the vibration responses of the representative bridge, Bridge Series A, Bridge Series B, and Bridge Series C to the two load cases are summarized in Table 3.In the table,    denotes the local sensitivity index of parameter  and    denotes the TSI obtained from the Sobol' method.It can be concluded that:

•
In the case of a single pedestrian crossing, both GSA and LSA select the same sensitive parameters and produce the same ranking for the representative bridge, Bridge Series A and Bridge Series B. However, for the bridges in Bridge Series C, the same sensitive parameters are selected but ordered differently.

•
In the case of the excitation of a crowd, both GSA and LSA select the same influential parameters and their ranking for the bridges in Bridge Series A with the mean natural

Local Sensitivity Analysis
The results of LSAs for Bridge Series C caused by the excitation of a crowd are shown in Figure 8b, which shows:

•
When the bridges have the mean damping ratio no bigger than 0.3%, there are six sensitive parameters, i.e., the bridge natural frequency, phase angle, step frequency, bridge damping ratio, time lag, and dynamic load factor.

•
For the bridges with the mean damping ratio ranging from 0.4% to 1.5%, the first five sensitive parameters are the same as those for the bridges with the mean value of damping ratio no bigger than 0.3%, as well as the ranking, but the sixth sensitive parameter is replaced by the bridge modal mass.Apart from these six parameters, an additional parameter, the dynamic load factor, becomes sensitive for the bridges with the mean damping ratio of 2%.
In summary, the results of GSA and LSA sensitive to the input space of damping ratio.Compared with GSA, LSA produces the different sensitive parameters for each damping ratio case.
Note that a 95th percentile of the peak of one-second RMS acceleration of the 10,000 cases for each analysis in the GSA and LSA in Section 5.3.2 is in the range of 0.47 m/s 2 to 4.31 m/s 2 .

Discussion
The results of GSA and LSA for the vibration responses of the representative bridge, Bridge Series A, Bridge Series B, and Bridge Series C to the two load cases are summarized in Table 3.In the table, S L x denotes the local sensitivity index of parameter x and S T x denotes the TSI obtained from the Sobol' method.It can be concluded that:

•
In the case of a single pedestrian crossing, both GSA and LSA select the same sensitive parameters and produce the same ranking for the representative bridge, Bridge Series A and Bridge Series B. However, for the bridges in Bridge Series C, the same sensitive parameters are selected but ordered differently.

•
In the case of the excitation of a crowd, both GSA and LSA select the same influential parameters and their ranking for the bridges in Bridge Series A with the mean natural frequency of 1.5 Hz and 2.5 Hz.The same sensitive parameters but ordered differently are generated for bridges in Bridge Series A with the mean natural frequency no larger than 1 Hz.However, for the representative bridge, bridges in Bridge Series A the mean natural frequency of 2 Hz and 3 Hz, Bridge Series B and Bridge Series C, GSA and LSA generate different sensitive parameters.

•
The results of GSA and LSA for all kinds of bridges in the two load cases are generally sensitive to the variation of input space of modal mass, natural frequency, and damping ratio of the bridge.The exception is the results of GSA and LSA for Bridge Series B under the excitation of a single pedestrian.

•
In addition, when the pedestrian step frequency is close to the bridge frequency, the TSIs of important parameters in all cases are found to be much greater than their firstorder counterparts, which indicates strong interaction between important parameters.This agrees with the fact that large human-induced vibration response only occurs when these two frequencies are close to each other.It is strong interaction between these two frequencies that leads to significant discrepancies between the TSIs and FSIs.

•
The damping ratio and natural frequency of a human body are found to be the only two insensitive parameters in all the cases.Additional numerical results, which are not presented in the paper, show that the human body's frequency may become sensitive only in the case where a single pedestrian crosses a bridge with a mean modal mass no larger than 350 kg.However, there are rarely real footbridges with modal masses lower than 350 kg.Therefore, there is no need to quantify the uncertainty in these two parameters in practice.

•
Sensitive parameters from both the GSA and LSA may vary with the load case.For example, for a footbridge in Bridge Series A with the mean natural frequency of 1.5 Hz, the step frequency is the most influential parameter in the case of a single pedestrian crossing, while the bridge's natural frequency becomes the most influential parameter in the case of a crowd crossing.This is explained by the fact that the model outputs for sensitivity analysis are different between the situations of a single pedestrian and a crowd.
The analysis presented in this paper indicates why there is inconsistency in the results of sensitivity analysis for the pedestrian-induced vibrations of footbridges in existing literature.These reasons are given as follows.

•
The results of sensitivity analysis may vary with the input space of parameters.In other words, for different footbridges and different load cases, different sensitive parameters may be produced by GSA or LSA.Therefore, the analysis in the paper is limited to the parameter space covered in the paper.

•
The inclusion of partial uncertain parameters in sensitivity analysis may lead to a misleading conclusion.Take the representative footbridge considered in Section 4 for example: additional analysis shows that the vibration response to the two load cases is most sensitive to step frequency if uncertainties in the step frequency, step length, and pedestrian weight are considered only.However, the bridge's natural frequency becomes most influential if all uncertain parameters are considered, as shown in Section 4.

•
The results of LSA may be inconsistent with those of GSA.This is attributed to the fact that LSA evaluates the sensitivity in the vicinity of a particular point and provides only a limited view of model sensitivity, which may exclude the effects of interactions between uncertain parameters in the unexplored input space and lead to misleading results.However, GSA provides more reliable results by characterizing the sensitivity across all input spaces.Therefore, GSA is suggested for the sensitive analysis of pedestrian-induced vibration of footbridges.

Conclusions
This paper investigates the impact of uncertainty in system and load parameters on vibration responses induced by walking pedestrian(s) using the approaches of LSA and Sobol'-based GSA.First, the LSA and GSA are conducted for the vibration responses of a representative bridge under two load cases, i.e., a single pedestrian and a crowd.Parametric analysis is then carried out to study the sensitivity of the results of LSA and GSA to the input spaces of the modal mass, natural frequency, and damping ratio of a footbridge.It is concluded that in the two load cases, GSA and LSA may produce consistent results in terms of sensitive parameters and their ranking, different sensitive parameters, or the same sensitive parameters but ordered differently, depending upon the bridge's modal properties.The results of GSA and LSA for all kinds of bridges in the two load cases are generally sensitive to the variation of input space of modal mass, natural frequency, and damping ratio of the bridge.The exception is the result for Bridge Series B under the excitation of a single pedestrian.The damping ratio and natural frequency of human body are found to be the only two insensitive parameters.Hence, there is no need to quantify experimentally or numerically the uncertainty in these two parameters in vibration response analysis for footbridges.Compared with the LSA, GSA provides more reliable results by characterizing the sensitivity across entire input spaces; therefore, GSA is suggested for the sensitive analysis of the pedestrian-induced vibration of footbridges.This study provides a comprehensive and conclusive insight into the sensitivity of a pedestrian-induced vibration response and explains why there are inconsistent results in existing literature.It is beneficial for selecting sensitive parameters required for analyzing uncertainty propagation in stochastic response assessment, as well as for the experimental characterization of randomness.

Funding:
The authors acknowledge the support of the National Natural Science Foundation of China (grant number 52178181) and the Fundamental Research Funds for the Central Universities, China (grant number 2020CX035).
Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Figure 1 .
Figure 1.The results from GSA and LSA of the vibration response of a representative footbridge induced by a single pedestrian: (a) the Sobol' indices for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 1 .
Figure 1.The results from GSA and LSA of the vibration response of a representative footbridge induced by a single pedestrian: (a) the Sobol' indices for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 2 .
Figure 2. The results from GSA and LSA of the vibration response of a representative footbridge induced by a crowd: (a) the Sobol' indices for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 3 .
Figure 3.The results from GSAs and LSAs of the vibration response of Bridge Series A with a single pedestrian crossing: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 3 .
Figure 3.The results from GSAs and LSAs of the vibration response of Bridge Series A with a single pedestrian crossing: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 4 .
Figure 4.The results from GSAs and LSAs of the vibration response of Bridge Series A induced by a crowd: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 4 .
Figure 4.The results from GSAs and LSAs of the vibration response of Bridge Series A induced by a crowd: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 5 .
Figure 5.The results from GSAs and LSAs of the vibration response of Bridge Series B with a single pedestrian crossing: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 5 .
Figure 5.The results from GSAs and LSAs of the vibration response of Bridge Series B with a single pedestrian crossing: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 6 .
Figure 6.The results from GSAs and LSAs of the vibration response of Bridge Series B induced by a crowd: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 6 .
Figure 6.The results from GSAs and LSAs of the vibration response of Bridge Series B induced by a crowd: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 7 .
Figure 7.The results from GSAs and LSAs of the vibration response of Bridge Series C with a single pedestrian crossing: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 7 .
Figure 7.The results from GSAs and LSAs of the vibration response of Bridge Series C with a single pedestrian crossing: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 8 .
Figure 8.The results from GSAs and LSAs of the vibration response of Bridge Series C induced by a crowd: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.

Figure 8 .
Figure 8.The results from GSAs and LSAs of the vibration response of Bridge Series C induced by a crowd: (a) the TSIs for system and load parameters; and (b) the COV of vibration response when one of the system and load parameters varies.
α > S L ls ; f b = 1.0 Hz, S L ϕ > S L fs > S L τ > α > S L ls ; f b = 1.5 Hz, S L f b > S L fs > S L ϕ ; f b = 2.0 Hz, S L f b > S L ϕ > S L fs > S L ξ b > S L τ > S L m b ; f b = 2.5 Hz, S L f b > S L fs > S L ϕ ; f b = 3.0 Hz, S L f b > S L ϕ ; f b = 0.5 Hz, S T fs > S T ϕ > S T τ > S T m b > S T ls > S T α ; f b = 1.0 Hz, S T fs > S T ϕ > S T m b > S T τ > S T f b > S T ls > S T α ; f b = 1.5 Hz, S T f b > S T fs > S T ϕ ; f b = 2.0 Hz, S T f b > S T fs > S T ϕ > S T τ > S T ls > S T m b ; f b = 2.5 Hz, S T f b > S T fs > S T ϕ ; f b = 3.0 Hz, ≤ ξ b ≤ 1.5%,

Table 1 .
Stochastic properties of load parameters.

Table 2 .
Statistical characteristics of system and load parameters.
1: The range for uniform distribution.

Table 3 .
The results from GSA and LSA of vibration responses of bridges.