Revising Seismic Vulnerability of Bridges Based on Bayesian Updating Method to Evaluate Tra ﬃ c Capacity of Bridges

: Analyzing the seismic vulnerability of bridge structures is of great signiﬁcance in assessing the sustainable capacity of bridges. However, previous seismic vulnerability analysis of bridge structures was based on theoretical simulation or experimental research, which couldn’t really reﬂect the actual seismic damage. Therefore, this paper used Bayesian updating method to obtain a corrected vulnerability curve of bridge structures considering the theoretical simulation and historical seismic data, which overcame the shortcomings that information couldn’t be updated. Speciﬁcally speaking, the seismic demand probability function of the bridge structure was obtained through numerical simulation, and a seismic vulnerability curve of a log-normal distribution with two parameters (median value c j and log-standard deviation ξ j ) was established, which was taken as the prior information. Then, combining the historical empirical damage data, the demand information is updated to obtain the posterior probability of the structural seismic demand. The possibility of exceeding the limit damage state of the bridge structure under di ﬀ erent ground motion intensity can be obtained by using the corrected seismic vulnerability curve. The results show that the corrected seismic vulnerability curves are closer to theoretical simulation results, which provides a new idea for the reliability analysis of the bridge structure. Moreover, based on the corrected seismic vulnerability curve of Guxigou Middle Bridge, the post-earthquake tra ﬃ c capacity of the bridge was evaluated by combining with tra ﬃ c ﬂow and tra ﬃ c time. The research results show that the proposed method can quickly evaluate the sustainable tra ﬃ c capacity of bridges after earthquakes, which has certain practicability and scientiﬁcity and provides the theoretical basis and practical guiding signiﬁcance for the relief work after earthquakes.


Introduction
In recent years, frequent earthquakes have resulted in huge economic loss and severe social impact on people's life. As an important part of the traffic lifeline project, bridges are widely distributed and are easily damaged by earthquakes. After the earthquake, many bridges will be directly damaged, which not only brings many problems to social life and economic operation, but also causes great difficulties in rescue and relief work after the earthquake. Therefore, analyzing the sustainable traffic capacity of bridges under earthquake disaster conditions is very important for emergency decision-making of traffic lifeline systems [1]. The Transportation Research Board (TRB) published the first edition of "Highway Capacity Manual (HCM)" [2] in 1950, which defined the concept of capacity. Most previous researches on the traffic capacity of bridges mainly relied on monitoring the health state of bridge structures. Cabboi et al. [3] proposed an automated procedure based on parametric identification methods for health monitoring for bridge structures. Chalouhi et al. [4] presented a novelty detection method for railway bridges, and based on a small a priori probability of anomaly in the bridge, the conditional/posterior probability of anomaly was computed by using the Bayesian inference. As we all know, the reliability of the bridge structures after earthquakes will directly affect the sustainable traffic capacity of bridges. However, effective disaster risk management strategies and practices depend on a rigorous understanding of the dimensions of exposure and vulnerability [5]. Barbat et al. [6] proposed seismic vulnerability and risk evaluation methods for urban areas considering the influence of multiple factors. Chieffo et al. [7] proposed that the seismic risk assessment is a multivariate problem based on the estimation of three major factors such as vulnerability, exposure, and hazard. In previous studies, seismic vulnerability analysis has been widely used in seismic reliability analysis of structures [8,9].
At present, the analysis of seismic performance of bridges under different working conditions usually takes the form of seismic vulnerability curves [10][11][12][13]. An empirical and theoretical method was used to construct the seismic vulnerability curves by researchers [14,15]. Laura et al. [16] used operational modal analysis and experimental modal analysis to compare with the numerical model, which was used to evaluate the load bearing capacity and structural behavior of a viaduct under traffic or seismic excitation. Zani et al. [17] built a detailed 3D finite-element (FE) model to analyze the seismic behavior of the Azzone Visconti Bridge by using nonlinear dynamic analysis method. According to Cornell's research [18], seismic demand and ground motion parameters-ground motion peak ground acceleration (PGA) followed the log-normal distribution. California Department of Transportation classified the damage states of bridge structures into five categories, including no damage (ND), slight damage (SD), moderate damage (MD), extensive damage (ED), and complete damage (CD) [19]. Shinozuka et al. [20] built the empirical vulnerability curve based on bridge damage data in the Kobe earthquake. Hwang et al. [21] used a displacement ductility ratio of piers as a damage degree to define the damage levels of the bridge, and the seismic vulnerability curves under four damage states were obtained. Choi et al. [22] established a theoretical seismic vulnerability curve for bridge structures in the mid-earthquake region. Duan et al. [23] proposed a modified coefficient model to evaluate the traffic capacity of the rescue road after the earthquake. Lan et al. [24] established the seismic vulnerability curves of beam bridges and arch bridges based on the damage data and proposed a method for evaluating the traffic capacity of common bridges after the earthquake. Liu [25] proposed a method for assessing the overall impact of the seismic vulnerability of urban bridges on the transportation network. However, all these methods have their deficiencies. Although the results of empirical analysis method are credible, these results are limited by lack of enough seismic information in most areas [26]. While the theoretical analysis method is used to develop the seismic vulnerability curve of bridge structures, it is conducted only through nonlinear analysis of the bridge structure. In addition, the information obtained through the above methods cannot be updated in a timely manner. Therefore, this study established firstly a seismic vulnerability curve of the logarithmic normal distribution with two parameters (median value c j and log-standard deviation ξ j ) based on the numerical simulation. Based on Bayesian updating method, combining the historical empirical damage data, the theoretical simulation results were taken as the prior information to update the seismic demand information of bridge structures. Finally, a corrected seismic vulnerability curve of bridges was drawn to evaluate the damage states of bridge structures under different ground motion intensity. Moreover, the sustainable traffic capacity of bridges after the earthquake was evaluated based on the corrected seismic vulnerability curve.

Seismic Vulnerability Analysis of Bridge
Based on existing research [18,[27][28][29][30][31], it is assumed that the probability distribution model of the seismic vulnerability curves can be expressed in the form of lognormal distribution function with two parameter (median value c j and log-standard deviation ξ j ), and the maximum likelihood method was used to estimate these two parameters.
where, F(a i ; c j , ξ j ) is the seismic vulnerability curve of the i-th bridge under the damage state j. a i is the peak ground acceleration (PGA) of the i-th bridge.
The likelihood function is expressed as follows: where, x i is a Boolean random variable of a given damage state of the i-th bridge at a i . The maximum likelihood estimations of two parameters can be obtained by solving the following equation [18,29,[32][33][34]: Two parameters of the seismic vulnerability curve calculated by this method are brought into Equation (1) to establish the seismic vulnerability curves under different damage states.

Bayesian Method Principle
According to probability theory and multiplication law [35], the probability of the product of two random variables is as follows: According to the total probability formula, if random variable ε is consisted of many unrelated events, there is: Bayesian theorem is obtained by bringing Equation (5) into Equation (4): where, p(θ/ε) is the posterior probability of parameter θ, p(θ) is the prior probability of parameter θ. The Bayesian theorem is expressed that the posterior probability of a random variable is the product of its prior probability, likelihood function, and normalization constant.

Parameter Estimation Method
Bayesian parameter estimation method uses prior distribution information (that is, numerical simulation) to make a preliminary estimation of the parameter distribution probability model. Then, based on the posterior data, which is historical empirical damage data of bridge structures, the probability distribution parameters are re-estimated and updated to make full use of existing data.
Due to the seismic demand and ground motion parameters-ground motion peak ground acceleration (PGA) follow the log-normal distribution [18], suppose that f (x) is a lognormal distribution density function of random variable x. µ and σ are the median value and standard deviation of data samples (lnx 1 , lnx 2 , . . . , ln x n ) based on seismic damage tests or theoretical simulation of n bridges, respectively. µ and σ are the median value and standard deviation of statistical samples (lnx 1 , lnx 2 , . . . , ln x m ) based on historical bridge earthquake damage of m bridges, respectively. If the prior distribution of µ is a lognormal function that is f (µ) ∼ LN(µ , σ ), the posterior distribution f (µ) of µ is expressed as follows [36][37][38]: According to Bayesian principle, two posterior lognormal distribution parameters (median value µ and standard deviation σ ) of the random variable x can be expressed as follows: So far, the median value µ and standard deviation σ based on Bayesian updating method are the two parameters (median value c j and log-standard deviation ξ j ) of Equations (1)-(3).

Seismic Vulnerability Curves Based on Historical Damage Data
At 14:18 on May 12, 2008, a magnitude 8.0M earthquake occurred Yingxiu, Wenchuan country, Sichuan province, China. Plenty of bridges suffered damage in the western area of Sichuan province. Lin [39] comprehensively investigated the seismic damage situation in the Wenchuan earthquake. By using these data and other supplementary data obtained from surveys, a total of 2150 bridges were finally sorted out, including 611 arc bridges and 1539 girder bridges. According to damage level, the damage states of these bridges were divided into four categories, including slight damage (SD), moderate damage (MD), extensive damage (ED), and complete damage (CD).
Due to the damage data include the structural vulnerability factors and uncertainty factors, the calculated seismic vulnerability curves ( Figure 1) are reasonable and reliable. Table 1 shows two parameters of the logarithmic normal distribution [39]. By counting the damage degree of different bridges during Wenchuan earthquake, the median value c j and log-standard deviation ξ j of the seismic vulnerability parameters can be obtained and used as the basis of empirical parameters. statistical samples ( , , … , ln ) based on historical bridge earthquake damage of bridges, respectively. If the prior distribution of μ is a lognormal function that is ( )~ ( , ), the posterior distribution ′′( ) of is expressed as follows [36][37][38]: According to Bayesian principle, two posterior lognormal distribution parameters (median value and standard deviation ) of the random variable x can be expressed as follows: So far, the median value and standard deviation based on Bayesian updating method are the two parameters (median value and log-standard deviation ) of Equations (1)~(3).

Seismic Vulnerability Curves Based on Historical Damage Data
At 14:18 on May 12, 2008, a magnitude 8.0M earthquake occurred Yingxiu, Wenchuan country, Sichuan province, China. Plenty of bridges suffered damage in the western area of Sichuan province. Lin [39] comprehensively investigated the seismic damage situation in the Wenchuan earthquake. By using these data and other supplementary data obtained from surveys, a total of 2150 bridges were finally sorted out, including 611 arc bridges and 1539 girder bridges. According to damage level, the damage states of these bridges were divided into four categories, including slight damage (SD), moderate damage (MD), extensive damage (ED), and complete damage (CD).
Due to the damage data include the structural vulnerability factors and uncertainty factors, the calculated seismic vulnerability curves ( Figure 1) are reasonable and reliable. Table 1 shows two parameters of the logarithmic normal distribution [39]. By counting the damage degree of different bridges during Wenchuan earthquake, the median value and log-standard deviation of the seismic vulnerability parameters can be obtained and used as the basis of empirical parameters.

Correcting the Theoretical Seismic Vulnerability Curve of Simply Supported Bridge
Zheng et al. [40] used the probabilistic seismic demand analysis approach to establish an analytical seismic vulnerability model of typically simply supported girder bridge (Guxigou Middle Bridge, see Figure 2) in Wenchuan area. Based on survey data of damaged bridges from the Wenchuan earthquake, the uncertainty of the bridge structure parameters was evaluated. The Latin Hypercube Sampling method was used to sample the uncertainty parameters in the model to generate 160 typical bridge analysis samples in turn. Nonlinear dynamic transient analysis was preformed using the Wenchuan Earthquake ground motions and the seismic peak responses of key bridge components were recorded. By using the regression analysis method, a relationship between the ground motion intensity and seismic demand of key bridge components was established. Figure 3 shows the probabilistic seismic demand models [33,40]. After determining the limit state of the bridge components corresponding to four damage states (SD, MD, ED, and CD) of the bridge, the seismic vulnerability curves corresponding to the different damage states were generated based on the logarithmic normal distribution assumption. Based on adopting the reliability theory, the upper and lower bounds of the seismic vulnerability of the bridge system were calculated. Finally, the least squares regression was used to obtain the median value and log-standard deviation of the seismic vulnerability function of the bridge system, as shown in Table 2 [40]. Figure 4 shows the seismic vulnerability curve of the Guxigou Middle Bridge based on the simulation.

Correcting the Theoretical Seismic Vulnerability Curve of Simply Supported Bridge
Zheng et al. [40] used the probabilistic seismic demand analysis approach to establish an analytical seismic vulnerability model of typically simply supported girder bridge (Guxigou Middle Bridge, see Figure 2) in Wenchuan area. Based on survey data of damaged bridges from the Wenchuan earthquake, the uncertainty of the bridge structure parameters was evaluated. The Latin Hypercube Sampling method was used to sample the uncertainty parameters in the model to generate 160 typical bridge analysis samples in turn. Nonlinear dynamic transient analysis was preformed using the Wenchuan Earthquake ground motions and the seismic peak responses of key bridge components were recorded. By using the regression analysis method, a relationship between the ground motion intensity and seismic demand of key bridge components was established. Figure  3 shows the probabilistic seismic demand models [33,40]. After determining the limit state of the bridge components corresponding to four damage states (SD, MD, ED, and CD) of the bridge, the seismic vulnerability curves corresponding to the different damage states were generated based on the logarithmic normal distribution assumption. Based on adopting the reliability theory, the upper and lower bounds of the seismic vulnerability of the bridge system were calculated. Finally, the least squares regression was used to obtain the median value and log-standard deviation of the seismic vulnerability function of the bridge system, as shown in Table 2 [40]. Figure 4 shows the seismic vulnerability curve of the Guxigou Middle Bridge based on the simulation.

Correcting the Theoretical Seismic Vulnerability Curve of Simply Supported Bridge
Zheng et al. [40] used the probabilistic seismic demand analysis approach to establish an analytical seismic vulnerability model of typically simply supported girder bridge (Guxigou Middle Bridge, see Figure 2) in Wenchuan area. Based on survey data of damaged bridges from the Wenchuan earthquake, the uncertainty of the bridge structure parameters was evaluated. The Latin Hypercube Sampling method was used to sample the uncertainty parameters in the model to generate 160 typical bridge analysis samples in turn. Nonlinear dynamic transient analysis was preformed using the Wenchuan Earthquake ground motions and the seismic peak responses of key bridge components were recorded. By using the regression analysis method, a relationship between the ground motion intensity and seismic demand of key bridge components was established. Figure  3 shows the probabilistic seismic demand models [33,40]. After determining the limit state of the bridge components corresponding to four damage states (SD, MD, ED, and CD) of the bridge, the seismic vulnerability curves corresponding to the different damage states were generated based on the logarithmic normal distribution assumption. Based on adopting the reliability theory, the upper and lower bounds of the seismic vulnerability of the bridge system were calculated. Finally, the least squares regression was used to obtain the median value and log-standard deviation of the seismic vulnerability function of the bridge system, as shown in Table 2 [40]. Figure 4 shows the seismic vulnerability curve of the Guxigou Middle Bridge based on the simulation.    In this paper, the seismic vulnerability function parameters of Guxigou Middle Bridge obtained based on the numerical simulation were used as the prior information of Bayesian correction. Based on Bayesian method, the structural seismic demand information was calculated and corrected by combining the damage survey data during Wenchuan earthquake. Based on Bayesian parameter estimation method, the posterior probability of seismic demand information was obtained by Equations (7)~(9). Figure 5 shows the contrast diagram of the corrected seismic vulnerability curve of Guxigou Middle Bridge based on numerical simulation and seismic survey data. Table 2 shows two parameters of the numerical simulation method and Bayesian updating method. It is found that the corrected seismic vulnerability curve is slightly lower than the seismic vulnerability curve of numerical simulation in Figure 5. The deviation of the corrected seismic vulnerability curve and the numerical simulation curve is caused by the uncertainty of the seismic wave. Therefore, the results of numerical simulation are safe to ensure the structural performance, safety, and stability under earthquakes. Figure 5 shows that the seismic reliability of the Guxigou In this paper, the seismic vulnerability function parameters of Guxigou Middle Bridge obtained based on the numerical simulation were used as the prior information of Bayesian correction. Based on Bayesian method, the structural seismic demand information was calculated and corrected by combining the damage survey data during Wenchuan earthquake. Based on Bayesian parameter estimation method, the posterior probability of seismic demand information was obtained by Equations (7)- (9). Figure 5 shows the contrast diagram of the corrected seismic vulnerability curve of Guxigou Middle Bridge based on numerical simulation and seismic survey data. Table 2 shows two parameters of the numerical simulation method and Bayesian updating method. In this paper, the seismic vulnerability function parameters of Guxigou Middle Bridge obtained based on the numerical simulation were used as the prior information of Bayesian correction. Based on Bayesian method, the structural seismic demand information was calculated and corrected by combining the damage survey data during Wenchuan earthquake. Based on Bayesian parameter estimation method, the posterior probability of seismic demand information was obtained by Equations (7)~(9). Figure 5 shows the contrast diagram of the corrected seismic vulnerability curve of Guxigou Middle Bridge based on numerical simulation and seismic survey data. Table 2 shows two parameters of the numerical simulation method and Bayesian updating method. It is found that the corrected seismic vulnerability curve is slightly lower than the seismic vulnerability curve of numerical simulation in Figure 5. The deviation of the corrected seismic vulnerability curve and the numerical simulation curve is caused by the uncertainty of the seismic wave. Therefore, the results of numerical simulation are safe to ensure the structural performance, safety, and stability under earthquakes. Figure 5 shows that the seismic reliability of the Guxigou It is found that the corrected seismic vulnerability curve is slightly lower than the seismic vulnerability curve of numerical simulation in Figure 5. The deviation of the corrected seismic vulnerability curve and the numerical simulation curve is caused by the uncertainty of the seismic wave. Therefore, the results of numerical simulation are safe to ensure the structural performance, safety, and stability under earthquakes. Figure 5 shows that the seismic reliability of the Guxigou Middle Bridge is higher when the earthquake intensity is lower, and its ability for maintaining the structural performance level is stronger. However, the seismic performance of this bridge is weaker when there are larger intensity earthquakes.

Correcting the Theoretical Seismic Vulnerability Curve of Continuous Girder Bridge
Wang et al. [41] selected a three-span pre-stressed reinforced concrete rigid frame bridge as the research object. The structural layout plan and elevation view are shown in Figure 6. Wang et al. used the structural nonlinear analysis software SAP2000 to analyze the dynamic characteristics and seismic response of the bridge and established a three-dimensional finite element model of the bridge. Based on the drift ratio of the bridge column, three limit states (ND, MD, and CD) of the bridge were taken as the limit states of the bridge, and the probability parameters of the three limit states were obtained. The incremental dynamic analysis method was used to calculate the probabilistic characteristics of the engineering demand parameters of the bridge at a given ground motion intensity. Based on the total probability theory, the probability was calculated when the bridge structure reached or exceeded the specified damage state at a given ground motion intensity. Finally, the seismic vulnerability curve of the bridge structure was obtained by fitting the discrete points. Table 3 lists the median value and log-standard deviation of the seismic vulnerability function of the bridge system based on numerical simulation [41]. Figure 7 shows the vulnerability curve of the bridge based on the simulation. Wang et al. [41] selected a three-span pre-stressed reinforced concrete rigid frame bridge as the research object. The structural layout plan and elevation view are shown in Figure 6. Wang et al. used the structural nonlinear analysis software SAP2000 to analyze the dynamic characteristics and seismic response of the bridge and established a three-dimensional finite element model of the bridge. Based on the drift ratio of the bridge column, three limit states (ND, MD, and CD) of the bridge were taken as the limit states of the bridge, and the probability parameters of the three limit states were obtained. The incremental dynamic analysis method was used to calculate the probabilistic characteristics of the engineering demand parameters of the bridge at a given ground motion intensity. Based on the total probability theory, the probability was calculated when the bridge structure reached or exceeded the specified damage state at a given ground motion intensity. Finally, the seismic vulnerability curve of the bridge structure was obtained by fitting the discrete points. Table 3 lists the median value and log-standard deviation of the seismic vulnerability function of the bridge system based on numerical simulation [41]. Figure 7 shows the vulnerability curve of the bridge based on the simulation.     In this paper, the seismic vulnerability function parameters (median value and log-standard deviation) of this bridge obtained based on the numerical simulation were used as the prior information of Bayesian correction. Based on the Bayesian parameter estimation method and combining survey damage data during the Wenchuan earthquake, the seismic demand information of the bridge structure was calculated according to Equations (7)~(9). The posterior probability of the post-earthquake demand of the bridge was obtained. Figure 8 shows a comparison of the seismic vulnerability curves of a three-span pre-stressed reinforced concrete continuous rigid frame bridge based on numerical simulation and seismic damage investigation data.  Figure 8 shows that the exceedance probability of the numerical simulation and the corrected seismic vulnerability curve based on seismic damage investigation data under the three working conditions is consistent and coincide with each other. Although there are some biases between the exceedance probability of the structural performance based on Bayesian updating method and the exceedance probability of the structural performance based on numerical simulation, the seismic performance of bridges is reliable, economical, and reasonable. In this paper, the seismic vulnerability function parameters (median value and log-standard deviation) of this bridge obtained based on the numerical simulation were used as the prior information of Bayesian correction. Based on the Bayesian parameter estimation method and combining survey damage data during the Wenchuan earthquake, the seismic demand information of the bridge structure was calculated according to Equations (7)- (9). The posterior probability of the post-earthquake demand of the bridge was obtained. Figure 8 shows a comparison of the seismic vulnerability curves of a three-span pre-stressed reinforced concrete continuous rigid frame bridge based on numerical simulation and seismic damage investigation data. In this paper, the seismic vulnerability function parameters (median value and log-standard deviation) of this bridge obtained based on the numerical simulation were used as the prior information of Bayesian correction. Based on the Bayesian parameter estimation method and combining survey damage data during the Wenchuan earthquake, the seismic demand information of the bridge structure was calculated according to Equations (7)~(9). The posterior probability of the post-earthquake demand of the bridge was obtained. Figure 8 shows a comparison of the seismic vulnerability curves of a three-span pre-stressed reinforced concrete continuous rigid frame bridge based on numerical simulation and seismic damage investigation data.  Figure 8 shows that the exceedance probability of the numerical simulation and the corrected seismic vulnerability curve based on seismic damage investigation data under the three working conditions is consistent and coincide with each other. Although there are some biases between the exceedance probability of the structural performance based on Bayesian updating method and the exceedance probability of the structural performance based on numerical simulation, the seismic performance of bridges is reliable, economical, and reasonable.  Figure 8 shows that the exceedance probability of the numerical simulation and the corrected seismic vulnerability curve based on seismic damage investigation data under the three working conditions is consistent and coincide with each other. Although there are some biases between the exceedance probability of the structural performance based on Bayesian updating method and the exceedance probability of the structural performance based on numerical simulation, the seismic performance of bridges is reliable, economical, and reasonable.

Post-Earthquake Traffic Capacity Assessment
As we all know, the transportation system in the earthquake area will suffer different degrees of damage after the earthquake. As an important part of the transportation system, the damaged degree and traffic capacity of the bridge system after the earthquake have a very important impact on the earthquake relief work. If the bridge structure is completely destroyed after the earthquake, its traffic is almost paralyzed, more air disaster relief forces must be considered. However, if the bridge structure is only slightly damaged after the earthquake, and ground transportation can still maintain a part of the operating state, then it may still be considered to rely mainly on the ground transportation system to complete the disaster relief transportation [24,42]. Therefore, it is of great significance to evaluate the post-earthquake state of bridges in the field of earthquake relief work.
In this paper, Guxigou Middle Bridge was taken as an example, which seismic vulnerability curve based on Bayesian updating method for each damage was used as the basis for assessing the traffic capacity after earthquake. Specifically, based on the corrected seismic vulnerability curves of the bridge, the percentages of various damage states under a certain ground motion were calculated by using the method provided in the reference [43]. Then, combining probabilistic results and expected values assigned for each damage level, the average damage level of the bridge under different ground motion intensity was calculated to obtain change of design traffic flow and variation curve of travel time under a certain traffic flow, which can provide a basis for the quick assessing of the traffic capacity of the bridge after an earthquake.

The Damage State Failure Probability of Guxigou Middle Bridge
Limit states (LS) describe the failure limit state of the structure, which is the limit of the adjacent damage state (DS), and the damage state is the interval between two adjacent limit states. Figure 9 shows the relationship between the limit state and damage state [34]. The limit state of Guxigou Middle Bridge structure was divided into four limit states, which includes SD (LS 1 ), MD (LS 2 ), ED (LS 3 ), and CD (LS 4 ). These four limit states divide the damage process of the bridge structure into five damaged states, which includes no damage (DS 1 ), slight damage (DS 1 ), moderate damage (DS 3 ), extensive damage (DS 4 ), and complete damage (DS 5 ) [34]. As we all know, the transportation system in the earthquake area will suffer different degrees of damage after the earthquake. As an important part of the transportation system, the damaged degree and traffic capacity of the bridge system after the earthquake have a very important impact on the earthquake relief work. If the bridge structure is completely destroyed after the earthquake, its traffic is almost paralyzed, more air disaster relief forces must be considered. However, if the bridge structure is only slightly damaged after the earthquake, and ground transportation can still maintain a part of the operating state, then it may still be considered to rely mainly on the ground transportation system to complete the disaster relief transportation [24,42]. Therefore, it is of great significance to evaluate the post-earthquake state of bridges in the field of earthquake relief work.
In this paper, Guxigou Middle Bridge was taken as an example, which seismic vulnerability curve based on Bayesian updating method for each damage was used as the basis for assessing the traffic capacity after earthquake. Specifically, based on the corrected seismic vulnerability curves of the bridge, the percentages of various damage states under a certain ground motion were calculated by using the method provided in the reference [43]. Then, combining probabilistic results and expected values assigned for each damage level, the average damage level of the bridge under different ground motion intensity was calculated to obtain change of design traffic flow and variation curve of travel time under a certain traffic flow, which can provide a basis for the quick assessing of the traffic capacity of the bridge after an earthquake.

The Damage State Failure Probability of Guxigou Middle Bridge
Limit states (LS) describe the failure limit state of the structure, which is the limit of the adjacent damage state (DS), and the damage state is the interval between two adjacent limit states. Figure 9 shows the relationship between the limit state and damage state [34]. The limit state of Guxigou Middle Bridge structure was divided into four limit states, which includes SD (LS1), MD (LS2), ED (LS3), and CD (LS4) . These four limit states divide the damage process of the bridge structure into five damaged states, which includes no damage (DS1), slight damage (DS1), moderate damage (DS3), extensive damage (DS4), and complete damage (DS5) [34]. Suppose that , ( ) represents the probability of exceeding the limit state LSi of the bridge structure when the ground motion intensity is x, which is the value ( ) of the seismic vulnerability function corresponding to LSi at x. , ( ) is called the limit state failure probability. Suppose that , ( ) represents the probability when the bridge structure is in the damage state DSi, which is the interval probability between the limit state LSi-1 and the limit state LSi of the bridge structure. , ( ) is called the damage state failure probability. Taking LS4 as an example, the limit state failure probability The seismic vulnerability curve (that is, the limit state failure probability curve) and correlations with the limit state and the damage state are shown in Figure 10. Suppose that P f ,LS i (x) represents the probability of exceeding the limit state LS i of the bridge structure when the ground motion intensity is x, which is the value F R i (x) of the seismic vulnerability function corresponding to LS i at x. P f ,LS i (x) is called the limit state failure probability. Suppose that P f ,DS i (x) represents the probability when the bridge structure is in the damage state DS i , which is the interval probability between the limit state LS i-1 and the limit state LS i of the bridge structure. P f ,DS i (x) is called the damage state failure probability. Taking LS 4 as an example, the limit state failure probability P f ,LS 3 (x) is equal to the sum of P f ,DS 4 (x) and P f ,DS 5 (x). The relationship between the limit state and the damage state can be expressed by Equation (10) [34].
The seismic vulnerability curve (that is, the limit state failure probability curve) and correlations with the limit state and the damage state are shown in Figure 10. According to the seismic vulnerability curves based on Bayesian updating method In Figure 5, the limit state failure probability of the bridge structure can be calculated by Equation (11). Table 4 gives the damage state failure probability of Guxigou Middle Bridge under different peak ground accelerations PGA. Figure 11 shows the failure probability curves of five damage states.   According to the seismic vulnerability curves based on Bayesian updating method In Figure 5, the limit state failure probability of the bridge structure can be calculated by Equation (11). Table 4 gives the damage state failure probability of Guxigou Middle Bridge under different peak ground accelerations PGA. Figure 11 shows the failure probability curves of five damage states.  According to the seismic vulnerability curves based on Bayesian updating method In Figure 5, the limit state failure probability of the bridge structure can be calculated by Equation (11). Table 4 gives the damage state failure probability of Guxigou Middle Bridge under different peak ground accelerations PGA. Figure 11 shows the failure probability curves of five damage states.

Design Traffic Flow Assessment
Based on existing research [43,44], the upper bound of the corrected seismic vulnerability curves of Guxigou Middle Bridge corresponding to each damage state were taken as the basis for calculating traffic flow. According to the references [42,43], the probabilistic results and expected values (the values ∈ [0, 4], where, 0 is no damage state and 4 is complete damage state) assigned for each damage level were combined to obtain the average damage level value l b of Guxigou Middle Bridge under different ground motion intensity, which can be calculated by Equation (12). l b = 0 * P f ,DS 1 + 1 * P f ,DS 2 + 2 * P f ,DS 3 + 3 * P f ,DS 4 + 4 * P f ,DS 5 (12) In this paper, the design speed v c of Guxigou Middle Bridge was 100 km/h, the number of lanes b was 3, and the minimum distance s min between the two vehicles at the design speed was 100 m. The designed critical flow f co was the critical value of traffic flow at the beginning of congested road conditions (that is the actual traffic capacity). According to Equation (13), the designed critical flow was 3000 vehicles/hour.
The change of traffic flow f c with the increase of PGA under the ground motion intensity can be calculated according to Equation (14) [43]. Table 5 shows the change of traffic flow f c changing with the peak ground acceleration PGA.
where, l b represents the average damage level of the bridge under different seismic intensities, a RC and β RC are related parameters, which generally take 0.03 and 4, respectively. PGA The change of traffic flow with the increase of the ground motion intensity were calculated according to Equation (14). Figure 12 shows that the normal performance of the bridge structure declines little, and the traffic flow does not change significantly when the earthquake acceleration is less than 0.10 g (that is, the earthquake intensity is small). When the seismic acceleration is more than 0.10 g, Guxigou Middle Bridge is greatly damaged and the bearing capacity of the bridge structure reduces with the increase of the ground motion intensity. Moreover, the main function and the traffic flow of the bridge rapidly decrease. When the earthquake acceleration reaches 0.30 g, the traffic flow reduces to 1441 vehicles/hour. As a result, the traffic flow of Guxigou Middle Bridge should be strictly controlled at this time, which provides a basis for the decision-makers to exercise the rescue and transportation after the earthquake.

Traffic Time Assessment
Traffic impedance refers to the relationship between road travel time and traffic load, intersection delay, and intersection load. Under earthquake action, the actual traffic capacity of each section reduces to a certain extent because the damage degree of each section is different. Therefore, the time required for the vehicles to pass through the bridges is also changed. Based on the calculation of the designed traffic flow, the time required for the vehicle to cross Guxigou Middle Bridge was calculated according to Equation (15) [42,43].
where, = v ⁄ is the time required for the vehicle to pass through the bridge during free flow. In the traffic flow distribution procedure of the US Highway Bureau, and are the blocking factors, which take 0.15 and 4, respectively. The impedance c on the section/bridge is the traffic time when the traffic flow is .
As the research object of this paper, the size of the main bridge of Guxigou Middle Bridge is 3m*30m. The time required for the vehicle to pass through the bridge under earthquake action is shown in Figure 13.

Traffic Time Assessment
Traffic impedance refers to the relationship between road travel time and traffic load, intersection delay, and intersection load. Under earthquake action, the actual traffic capacity of each section reduces to a certain extent because the damage degree of each section is different. Therefore, the time required for the vehicles to pass through the bridges is also changed. Based on the calculation of the designed traffic flow, the time required for the vehicle to cross Guxigou Middle Bridge was calculated according to Equation (15) [42,43].
where, c 0 = L/v c is the time required for the vehicle to pass through the bridge during free flow. In the traffic flow distribution procedure of the US Highway Bureau, α CF and β CF are the blocking factors, which take 0.15 and 4, respectively. The impedance c on the section/bridge is the traffic time when the traffic flow is f .
As the research object of this paper, the size of the main bridge of Guxigou Middle Bridge is 3 m × 30 m. The time required for the vehicle to pass through the bridge under earthquake action is shown in Figure 13.

Traffic Time Assessment
Traffic impedance refers to the relationship between road travel time and traffic load, intersection delay, and intersection load. Under earthquake action, the actual traffic capacity of each section reduces to a certain extent because the damage degree of each section is different. Therefore, the time required for the vehicles to pass through the bridges is also changed. Based on the calculation of the designed traffic flow, the time required for the vehicle to cross Guxigou Middle Bridge was calculated according to Equation (15) [42,43].
where, = v ⁄ is the time required for the vehicle to pass through the bridge during free flow. In the traffic flow distribution procedure of the US Highway Bureau, and are the blocking factors, which take 0.15 and 4, respectively. The impedance c on the section/bridge is the traffic time when the traffic flow is .
As the research object of this paper, the size of the main bridge of Guxigou Middle Bridge is 3m*30m. The time required for the vehicle to pass through the bridge under earthquake action is shown in Figure 13.
0% 20% 40% 60% 80% 100% 120% 140% 160% 180% 200%   Figure 13 shows that there is no fluctuation in the time of opening when the traffic flow is below the critical flow (that is the percentage of traffic flow is less than 100%). While the flow traffic exceeds the critical traffic flow, traffic time curve is increased with the increase in traffic flow. Road congestion situation intensifies, which affects the normal passage of vehicles. When the bridge is damaged during earthquakes, the structure is damaged, affecting the normal service performance and the normal traffic, the passing time of vehicle is increased. According to Equation (15), the traffic time changes with the earthquake acceleration under certain traffic flow is calculated, and the result is shown in Figure 14. Figure 13 shows that there is no fluctuation in the time of opening when the traffic flow is below the critical flow (that is the percentage of traffic flow is less than 100%). While the flow traffic exceeds the critical traffic flow, traffic time curve is increased with the increase in traffic flow. Road congestion situation intensifies, which affects the normal passage of vehicles. When the bridge is damaged during earthquakes, the structure is damaged, affecting the normal service performance and the normal traffic, the passing time of vehicle is increased. According to Equation (15), the traffic time changes with the earthquake acceleration under certain traffic flow is calculated, and the result is shown in Figure 14.  Figure 14 shows the traffic time is almost unchanged when the peak acceleration of earthquake is below 0.20g. At the same time, it does not affect the normal service performance of the bridge. When the peak acceleration of earthquake increases from 0.20g to 0.30g, the traffic time changes significantly, and the bridge structure is greatly damaged, and the normal passage of vehicles is blocked. When the rescue is needed after the earthquake, the traffic diversion control should be considered to avoid causing congestion in the prime time for rescue. When the peak acceleration of earthquake is greater than 0.30g, the traffic time is of rapid growth. At this time, the bridge structure suffered severe damage or even lost its normal using performance and the carrying capacity could not operate normally. Therefore, the earthquake relief should consider changing the rescue route or using air traffic force to avoid this traffic node.

Conclusion
(1) This study used Bayesian correction method to update the seismic vulnerability curves of bridges under several damage states. In the case that the experimental conditions are not sufficient for statistical probability analysis, this method makes full use of the existing historical earthquake damage experience and numerical simulation results to determine the probability distribution parameters of the bridge structure. The results show that the corrected seismic vulnerability curves are closer to the actual situation, which provides a new idea for the reliability analysis of bridge structure.
(2) Based on the corrected seismic vulnerability curves of Guxigou Middle Bridge, the traffic flow and traffic time of the bridge after earthquakes were quickly evaluated. The results show that the traffic flow decreases and the traffic time increases when the peak acceleration of earthquake increases. This method has certain practicability and scientificity, which is of great significance for evaluating the sustainable traffic capacity of bridges after earthquakes.   Figure 14 shows the traffic time is almost unchanged when the peak acceleration of earthquake is below 0.20 g. At the same time, it does not affect the normal service performance of the bridge. When the peak acceleration of earthquake increases from 0.20 g to 0.30 g, the traffic time changes significantly, and the bridge structure is greatly damaged, and the normal passage of vehicles is blocked. When the rescue is needed after the earthquake, the traffic diversion control should be considered to avoid causing congestion in the prime time for rescue. When the peak acceleration of earthquake is greater than 0.30 g, the traffic time is of rapid growth. At this time, the bridge structure suffered severe damage or even lost its normal using performance and the carrying capacity could not operate normally. Therefore, the earthquake relief should consider changing the rescue route or using air traffic force to avoid this traffic node.

Conclusion
(1) This study used Bayesian correction method to update the seismic vulnerability curves of bridges under several damage states. In the case that the experimental conditions are not sufficient for statistical probability analysis, this method makes full use of the existing historical earthquake damage experience and numerical simulation results to determine the probability distribution parameters of the bridge structure. The results show that the corrected seismic vulnerability curves are closer to the actual situation, which provides a new idea for the reliability analysis of bridge structure. (2) Based on the corrected seismic vulnerability curves of Guxigou Middle Bridge, the traffic flow and traffic time of the bridge after earthquakes were quickly evaluated. The results show that the traffic flow decreases and the traffic time increases when the peak acceleration of earthquake increases. This method has certain practicability and scientificity, which is of great significance for evaluating the sustainable traffic capacity of bridges after earthquakes.