Abstract

This study introduces a method to obtain accurate results regarding bridge seismic cost and reduce structural seismic losses. A three-span reinforced concrete (RC) continuous bridges model consistent with the reality was established, and seismic vulnerability of the RC bridge was calculated via incremental dynamic analysis combined with local earthquake disaster data. Direct losses from earthquakes are calculated based on local earthquake damage condition, specific maintenance reinforcement methods, and analysis of seismic vulnerability. The indirect economic losses caused by inspection and maintenance are calculated based on local traffic volume. To reduce the effects of bridge aging and deterioration, this paper examines periodic inspection of bridges after earthquakes. By calculating the life cycle cost of structure under different inspection frequencies, the optimal time interval was determined. The results presented for this study can serve as references for the calculation of structural seismic loss and thereby provide the basis for strategies of reducing seismic loss.

1. Introduction

Bridges are often subjected to various events such as earthquakes, fires, and traffic accidents [1]. Earthquakes are common and can cause serious economic losses; hence, seismic risk assessment in transportation networks has become a popular research area [2, 3]. To reduce earthquake-related economic losses for bridges and to identify reasonable inspection and maintenance intervals, it is necessary to evaluate seismic loss precisely.

Studies related to earthquake loss assessment are generally concentrated in the field of architecture. Probability analysis can reflect the variability of ground motion parameters [46]; seismic loss assessment is thus evolving from deterministic parsing to probabilistic analysis. Luo et al. [7] provided a seismic loss calculation method using an RC frame structure with different seismic fortification levels based on PEER-PBEE to explore ground motion intensity, the applicable vulnerability function, and the loss function, all of which should be considered in earthquake loss evaluation. Kiremidjian et al. [8] conducted a study on direct losses resulting from bridge failure and driving delays under specific earthquake conditions. In addition, Chang et al. [9] developed a method to identify a series of seismic scenarios to estimate regional seismic intensity. Wang and Sun [10] analyzed an existing bridge’s mechanical behavior and vulnerability. Zhou et al. [11] selected common earthquake damage factors and obtained the corresponding relationships between the damage grade and earthquake damage index via multiple linear regression and the least squares method. Li et al. [12] investigated and calculated cost data for a city bridge project. Mahmoud et al. [13] introduced a new double-stage static method to estimate capacity curves of MR frames; in addition to enhanced accuracy, the N1-N2 method was found to be exceptionally efficient. Only recently have the effects of bridge columns’ cumulative seismic vulnerability been considered in bridges’ life cycle cost analysis [1416]. The loss ratio obtained by a special study was used to calculate earthquake losses of railway bridges; however, it is not accurate to calculate direct structural losses from earthquakes via the loss ratio. Establishing a model that corresponds to reality to quantify bridge-related seismic loss remains challenging. Generally, special inspections are only carried out for bridge structures after moderate earthquakes, although earthquakes can happen at any time. It is therefore necessary to identify reasonable inspection intervals, which can be added as special inspections and significantly reduce earthquake losses.

In this paper, an accurate model and structural vulnerability curve obtained via incremental dynamic analysis (IDA) is presented. The results of direct economic structural loss after earthquakes are calculated based on seismic data at the case site, statistical analysis of various maintenance and reinforcement measures taken after the earthquake, city bridges’ engineering cost data, and statistical analysis of local traffic information to calculate indirect economic losses after earthquakes. Assumptions for regular inspection and maintenance activities include 1-, 2-, 3-, 4-, and 5-year time intervals and calculating the cost of a bridge at every interval in its service life to determine the optimal inspection interval.

2. Methods

The first step in assessing seismic loss to determine the optimal inspection frequency is to identify the seismic parameters, source characteristics of earthquakes, and types of earthquakes throughout history in the region. Seismic hazard analysis follows from identifying the magnitude and frequency of earthquakes in a region over a certain period of time [4, 9, 17]. Then, fragility curves can be obtained by analyzing seismic vulnerability of structure [1820]. For a specific earthquake intensity, the probability of a structure being in various states of damage after an earthquake can be determined using fragility curves. The second step is statistical analysis of repair strategies, preferably those associated with different damage states, such as strategies’ unit cost, and traffic information (average daily traffic, downtime associated with different damage levels, and so on). An increase in the damage state can reduce the link traffic capacity and speed limit, resulting in additional travel time. The amount of traffic volume that exceeds the capacity of the damaged link has to follow the detour [21]. Next, it is necessary to compute the direct and indirect losses of the affected structure associated with the bridge’s damage state to calculate the cost of routine maintenance and inspection at certain intervals in the bridge’s service life. Finally, to determine the optimal inspection interval, overall performance requirements must be satisfied while the total resource costs remain stable. A flow chart summarizing this methodology is presented in Figure 1.

3. Bridge Modeling

The bridge in this study is designed for continuous T girder with a total length of 3 × 25 m, four traffic lanes, and a bridge width of 21 m. The main beam is composed of seven T beams made of concrete C40 (see Figure 2). Static and dynamic analyses of the structure were carried out in OpenSees (http://opensees.berkeley.edu/) to evaluate the comprehensive performance of the postearthquake bridge. The bent cap should be designed as a capability protection component [22], and the T girder and bent cap are modeled using linear elastic beam-column elements.

The pier is a double-column pier with a diameter of 1.4 m, a height of 10 m, and C30 material (see Figure 3), which is modeled by a nonlinear three-dimensional beam-column element. The concrete cross section is discretized into several fibers. Each longitudinal reinforcement is discretized into a fiber. Thirty steel fibers of unconfined concrete are radially divided into 2 layers; there are 16 portions in the ring for a total of 32 nonconstrained concrete fibers. Constrained concrete is radially divided into 10 layers, and the ring is further divided into 16 portions (cf. Figure 2). The moment-curvature relationship of the column cross section is also represented in Figure 2. The columns on pile cap are supported by driven pile foundations with a length of 10 m. The pier-to-deck connection is a plate-type elastomeric bearing, which is modeled by an elastomeric bearing element [23].

The abutment is supported by the pile foundation. The total height of the abutment is 9.7 m, the height of the back wall is 2.2 m, the depth of the reinforced concrete rock-socketed pile is 8.5 m, and the pile diameter is 1.5 m. Interaction between abutment and backfill is considered via the hyperbolic force-displacement (HFD) model [24]. The HFD abutment model is simulated by Hyperbolic Gap Material in OpenSees. And the parameters of the HFD abutment model are based on the study of Wilson [25]. The piles are modeled by a linear elastic beam-column element, and the soil-structure interaction is simulated via equivalent soil springs with linear behavior [26]. The abutment is connected to the deck by a lead rubber bearing [23]. The piers are most likely to be damaged under earthquake excitation. In this paper, the seismic losses of the piers after an earthquake are calculated accurately. Other components can be calculated with reference to the piers.

4. Seismic Vulnerability Analysis

4.1. Seismic Hazard Curve

Seismic hazard is the size and frequency of seismic action to which a region may be subjected for a certain period of time. The curve reflects the seismic parameters, source characteristics, and types of earthquakes in the region. The results of the seismic hazard analysis are shown in Table 1 and Figure 4.

When selecting different functions in MATLAB to fit the data in Table 1, the error of logistic regression in Table 1 was smaller, and the fitting formula was simple. The curve and function expression of the earthquake hazard are shown in Figure 4.

4.2. Probabilistic Seismic Vulnerability Analysis

When conducting a vulnerability analysis of the pier under earthquake conditions, the column’s curvature ductility is used as the damage index of the structure during an earthquake based on the procedure given by Aryan and Ghassemieh [27]. The curvature ductility is defined as the ratio of the maximum column curvature recorded from a nonlinear time-history analysis to the column yield curvature obtained from moment-curvature analysis (Figure 5) [28]. The current paper defines four damage states based on the following [29]: slight damage, s = 1; moderate damage, s = 2; major damage, s = 3; and complete damage, s = 4. The damage limit states are assumed to equal the ductility of 1.0, 2.0, 4.0, and 7.0 for slight, moderate, major, and complete damage states, respectively [28]. Selecting 50 seismic waves from Appendix A of the literature [30], when scaling the PGA of each seismic wave to 0∼1.5 g, the step length is 0.1 g [3133], and the bridge structure is analyzed via the IDA method. The time-history response of the pier curvature is recorded. For every working condition, the pier’s top section has 50 time responses corresponding to the seismic waves; it only draws the curvature response diagram of the pier’s top section under the action of two seismic waves that are limited in space (cf. Figure 5). Based on the statistical characteristics of the dynamic response of the structure under different seismic waves, the curvature distribution of piers at the corresponding peak acceleration is obtained. Here, we give the case of PGA = 0.5 g, 1.0 g, and 1.5 g (Figure 5). Then, the damage to the structure is judged according to the defined damage level. If the total number of seismic waves is , and the number of seismic waves to which the curvature of the pier is subjected exceeding the specified damage index is , then the failure probability under the specified damage level can be obtained as follows [34]:

The vulnerability curve can be obtained after obtaining the failure probability of each damage level, but the vulnerability curve is a wave line that is not smooth. Fitting correction for the vulnerability curve via logarithmic normal cumulative distribution (Formula (2)) and seismic vulnerability curves for specific peak ground acceleration [35] is shown in Figure 5.where and are the median value and logarithmic standard deviation of PGA associated with a damage state, respectively, and is the standard normal cumulative distribution function. To avoid any intersection, the fragility curves in different damage states result in different deviation values, and Zhang and Huo [36] suggested considering a common deviation value for all damage states. In this study, the median value of log-standard deviations (equal to 0.58) is selected as the identical log-standard deviation as shown in Table 2.

5. Earthquake Loss Assessment

5.1. Calculation of Earthquake Loss

Assuming that the life of a bridge is 100 years, the frequency of inspection and maintenance is 1, 2, 3, 4, and 5 years, and the life cycle cost of bridges under earthquake conditions is calculated, respectively. The life cycle cost of the bridge includes routine inspection and maintenance costs at certain time intervals and the costs associated with different damage states of the bridge subjected to earthquakes . The mentioned costs represent the major seismic costs of the bridge from project completion to the end of its service lifetime and can be calculated as [28]where is inspection cost and is maintenance cost. is indirect cost of the maintenance procedure, including detour cost and a time loss for users and goods when traveling and is direct postearthquake cost. is indirect postearthquake cost, including detour cost and a time loss for users and goods when traveling. The recurrent costs are combined by weighting them using a discount factor, , that takes into account the time value of money [28], and r represents the after-inflation risk-free discount ratio to convert future costs into present values; in this article, r assumes to be 5%.

If the time interval of regular inspection and maintenance activities is 1 year, then the process of calculating structural earthquake loss is shown in Figure 6. The direct cost of regular inspection and maintenance calculated by the function in Figure 6 is ¥560160. Traffic restrictions are required in the process of bridge maintenance and reinforcement, as users are forced to detour [21]. The detour cost can be calculated as follows (based on Stein et al. [6]):where and are the average costs for running cars and trucks per unit length, respectively, is the daily traffic volume of the bridge in the interval, is the average daily truck traffic ratio in the interval, and is the length of the detour. is the downtime caused by maintenance, and is index of usage disruption. In the current study, only 1/4 of a bridge was closed each time for 1 week [28]. The monetary value of time loss for users and goods while traveling can be expressed as (based on You et al. [21])where is the average wage per hour, is the average total compensation per hour, is the time value of the goods, and are the average vehicle occupancy for cars and trucks, respectively, is the traffic speed on the undamaged bridge, is the traffic speed on the roundabout line, and is the length of the bridge.

The direct loss in postearthquake costs is related to the maintenance measures taken after the earthquake. Assume that a bridge will return to its original state after each inspection and maintenance procedure without considering the effects of structural deterioration. Considering the different maintenance and reinforcement schemes, assume the seismic loss of different structural damage states is calculated accurately, and each result is multiplied by the appropriate weighting factor. The general formula for direct losses in postearthquake costs is shown in Figure 6, where represents the overall performance of the bridge under a specified seismic hazard risk and can be determined based on the annual probability of exceeding a given damage state [28], expressed as

Convolved with the bridge’s site-specific hazard curve, the exceedance probability of different damage states of the pier is obtained by Formula (6):

Through the above calculation, when the time interval of inspection and maintenance is 1 year, the cost of regular inspection and maintenance is ¥7005360, the postearthquake loss is ¥3911367, and the life cycle cost is ¥10916727. When the regular interval of inspection and maintenance activities changes from every year to every 5 years, the results are as shown in Table 3.

5.2. Basic Parameters

The preferred repair strategies [3742] for a pier subjected to an earthquake and the repair percentages of different strategies are shown in Table 4.

The unit costs of different repair items can be obtained based on bridge construction costs [29], the market price of building materials, and the project quota of different repair strategies (Table 5) [3]. In other cases, the cost of maintenance and reinforcement should be addressed in light of actual local conditions.

The damage ratio index is assumed to be 0.03, 0.08, 0.25, and 1.0 for the slight, moderate, major, and complete damage states, respectively [3942]. For example, when the damage ratio index is 0.03 for the slight damage state, the length of the crack in the pier to be repaired is 0.03 × 10 = 0.3 m. The traffic information calculated in this region is shown in Tables 68.

5.3. Determination of Optimal Inspection Interval

The frequency of inspection and maintenance is 1 year, and the various loss costs of bridges under earthquake conditions are shown in Figure 7. The results of the direct and indirect structural costs in the bridge service life at different time intervals are shown in Figure 8.

As indicated, the indirect postearthquake cost and indirect maintenance cost are higher than the direct cost. The ratio of indirect postearthquake cost to direct postearthquake cost is 9.4, and the ratio of indirect cost of maintenance to the direct cost is 11.5, 10.4, 9.6, 8.6, and 8.2 at inspection and maintenance intervals of 1, 2, 3, 4, and 5 years, respectively. Routine inspection and maintenance costs are higher than postearthquake costs.

Although the bridge length is short, the route is the key traffic line, and the traffic volume is large. When calculating direct postearthquake costs, the unit costs of different repair items used in maintenance represent the base price; thus, the indirect costs of maintenance are higher than the direct.

As shown above (Figures 910), the shorter the time interval, the higher the regular inspection cost of the structure. When the time interval changes from 1 to 2 years, the total cost is reduced significantly; when the time interval changes from 4 to 5 years, the change in total cost is not as apparent. In this case, the optimal time interval is 5 years.

6. Conclusions

This article presents a computational framework for the seismic vulnerability, seismic loss assessment, and frequency analysis of optimal routine bridge inspection. The optimal frequency for regular inspection is obtained based on economic metrics. Assuming the bridge will return to its original state after each maintenance procedure, seismic loss estimation is conducted via a framework based on seismic vulnerability analysis. The proposed approach also accounts for inspection and maintenance costs, uncertainty in repair procedures, and indirect losses. This methodology can be effectively used to arrive at more accurate estimates of economic losses for bridges. The seismic vulnerability of the structure is obtained by establishing an accurate model that accounts for the probability of a set of seismic scenario events that reflect earthquake activity in the region. The following conclusions are drawn:(1)Routine inspection and maintenance are introduced to minimize loss of life and property caused by earthquakes and to reduce the effects of aging, deterioration, and corrosion. A computational framework is presented for life cycle costs at different time intervals of routine inspection, and this paper identifies an optimal frequency of regular inspection by evaluating the associated monetary values.(2)Comparing calculated life cycle costs, it is found that as the time interval becomes larger, the total cost tends to remain stable. The most reasonable time routine inspection and maintenance interval is the minimum value of the stable value corresponding to the time interval.(3)The information obtained from this approach can be used in maintenance and inspection procedures. The ratio of the indirect maintenance costs to direct costs is approximately 9.2; more attention can hence be paid to indirect costs.(4)The life cycle cost is sensitive to the parameters used in earthquake loss assessment; therefore, the parameters involved in cost analysis should be estimated carefully.(5)As the time interval grows, the effects of aging and deterioration on structures’ seismic loss assessment cannot be neglected. Future research will concentrate on time-dependent seismic vulnerability in bridge seismic loss estimation.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (nos. 51578184 and 51408167) and the Postdoctoral Scientific Research Developmental Fund of Heilongjiang Province (no. LBH-Q15058).