1. Introduction
The current transportation market, driven by a steady increase in global, is growing rapidly driven by a steady increase in the global economy, especially in developing countries [
1]. There is an urgent demand for constructing highways and railways in mountainous and cross-sea areas, where long-span bridges are widely built or in construction [
2]. Cable-supported bridges, including cable-stayed bridges and suspension bridges, are widely used in highways crossing gorges, rivers, and gulfs, due to their superior structural mechanical property and beautiful appearance [
3]. A cable-supported bridge is usually the key junction of a highway or a railway demanding a higher safety margin. However, long-span bridges suffer from harsh environmental effects and complex loading conditions, such as heavy traffic loading, significant wind load effects, severe corrosion effects, and other natural disasters [
4,
5,
6]. Subsequently, these effects may result in changes in the structural-mechanical behavior, which may alter the dynamic characteristics and resistance of the bridge during its lifetime. Due to the structural flexibility of cable-supported bridges, wind load and seismic effect are critical factors affecting bridge safety [
7,
8]. Furthermore, the steadily growing traffic load and environmental corrosion effects can lead to bridge safety degradation [
9,
10]. Therefore, it’s critical to ensure the structural safety of cable-supported bridges in harsh environments.
Cable-supported bridges have one or more pylons and flexible decks supported by prestressed cables. The cables are critical components to support girder gravity and moving vehicle loads [
11]. However, cables are subjected to corrosion under the coupling effect of fatigue loading and environmental corrosion. As a result, cable corrosion can lead to strength degradation, which can result in cable rupture or even collapse of the bridge. A research conducted by Mehrabi et al. [
12] revealed that more than half of the stay cables of an in-service cable-stayed bridge were severely corroded, demanding prompt replacement. Cable-supported bridges are associated with numerous uncertainties, such as material performance behaviors, cross-sectional dimensions, and various external loads. These uncertainties make the safety analysis of cable-supported bridges more complex for two main reasons. First, the structural mechanical behavior is highly nonlinear. Second, the bridge structure is statically indeterminate, which may lead to a failure sequence of more complicated systems.
System reliability theory is a powerful tool for the estimation of the life-cycle safety of cable-supported bridges accounting for uncertainties [
13,
14]. For a statically indeterminate structure, the scenario of a component failure may not result in safety issues. For example, the rupture of a stay cable may not affect the safety and serviceability of a cable-stayed bridge. However, a progressive failure of critical components may result in safety problems or even a bridge collapse event. The system reliability theory can be employed to search for potential failure sequences and identify sensitive components or parameters, which is beneficial for design optimization. There are two critical problems in structural system reliability evaluation. First and foremost, long-span bridges have a higher indeterminate degree that makes the structural performance functions more complicated. Thus, it is important to accurately and efficiently approximate the limit state function by utilizing an adequate approach. In this regard, there are several commonly used approaches, such as the artificial neural network (ANN) and the respond surface method (RSM). With the rapid development of artificial intelligence (AI), the application of AI-based approaches in structural reliability analysis is becoming increasingly popular. Dai et al. [
15] proposed a least squares support vector regression (LS-SVR) for structural reliability analysis with high efficiency and accuracy. Lu et al. [
16] utilized a machine learning approach to estimate the system reliability of a cable-stayed bridge considering cable damage due to fatigue and corrosion effects. Deep learning is a recently developed ANN approach that is widely used in the fields of graphics and image recognition [
17]. A deep belief network (DBN) is a type of deep learning method that is commonly used for probability training [
18]. DBN is feasible for approximating the structural mechanical behavior, which is the basis for subsequent reliability evaluation. Compared to traditional AI-based approaches, DBN has the following advantages [
19,
20,
21]. First, it utilizes hidden layers which are more efficient compared to traditional multilayer perceptions, such as the back-propagation (BP) neural networks, and thus the training process is time-saving. Second, DBN has a deep architecture like a mammal brain, which is able to capture the depth nonlinearity of a structural complex system, while SVM is a shallow architecture benefiting from the nonlinear kernel function. Finally, DBN is more robust in the classification of various data with large-scale parameters.
The second critical problem is how to search for the failure sequence of a system efficiently. In this regard, the most commonly used approach is the β-bound approach that utilizes the reliability index range to determine the components that are likely to fail. Li and Song [
22] presented a branch-and-bound method together with the finite element (FE) method to update the limit state function. Kang et al. [
23] developed a matrix-based system reliability method to replace the failure tree of complex structures. Liu et al. [
24] proposed an adaptive SVR approach combined with the advanced β-bound approach to develop the fault tree of a pre-stressed concrete cable-stayed bridge.
This study proposes a novel and intelligent approach for system reliability evaluation of cable-supported bridges under stochastic traffic load by utilizing DBNs. Mathematical models for the system reliability of cable-supported bridges were derived taking into consideration structural nonlinearities and high-order statically indeterminate characteristics. A computational framework was developed to exhibit the procedures for system reliability evaluation using DBNs. In a case study, a prototype suspension bridge was selected to investigate the system reliability under stochastic traffic loading considering site-specific traffic monitoring data. Numerical results indicated the feasibility of utilizing DBN as a meta-model to estimate the structural failure probability. The dominant failure mode of the suspension bridge was determined, and the effect of cable degradation due to fatigue-corrosion damage on the system reliability of the suspension bridge was investigated.
3. A reliability Evaluation Framework Based on DBNs
3.1. Theoretical Basis of DBNs
In general, a DBN is a new type of generative neural networks which are constructed by stacking restricted Boltzmann machines (RBMs) that are usually trained based on the probability principle [
29]. In a DBN, the RBMs are connected between hidden and visible layers. In the visible layer, the nodes are connected without directions. In the pre-training stage, the main training strategy is unsupervised learning, where the fine-tuning stage is trained by a supervised learning method. The outputs of a visible layer are treated as inputs to the hidden layer. In general, DBNs are widely used in the field of image identification. Their application in structural reliability analysis is relatively insufficient. In this study, the DBN and BP networks were utilized to approximate the structural response surface functions. The framework of predicting the effect of structural load effect using the DBN is shown in
Figure 4. The DBN principle will be demonstrated below.
The energy function of a joint structure (
V,
H) of visible and hidden units can be written as:
where
V = (
v1,
v2,
v3,…,
vm) is a visible layer,
H = (
h1,
h2,
h3,…,
hn) is a hidden layer,
m and
n are the number of element in the visible and the hidden layer, respectively;
vi and
hj are the
ith and
jth neurons in the visible and hidden layers, respectively;
w = (
wij) is the weight of elements, and
ai and
bj are the bias in the visible layer and the hidden layer, respectively.
As it can be seen in
Figure 4, the parameters
ai, bi, and
wij are trained and optimized from bottom to top based on an unsupervised learning method. Subsequently, the parameters are adjusted from top to bottom. The training function is:
where
P is a probability function, ζ is the training rate which is in the range of 0.05 and 0.2, and
hj(k) is the
kth optimization of the
jth neuron in the hidden layer.
3.2. Proposed Computational Framework
In
Figure 5, a computational framework is presented to illustrate the approach of utilizing the DBN to estimate the reliability of a cable-supported bridge. Detailed descriptions of the followed computational steps are summarized below.
In the first step, the critical random variables, such as the cable and concrete elasticity modulus, and the cross-sectional inertia moments of girders and pylon elements, are selected. The training samples should be selected according to the probability distribution of these random variables, which will be used to train the mechanical behavior of the bridge structure.
In the second step, a parametric finite element model is developed using a commercial software, such as ANSYS or MIDAS. It is recommended to use ANSYS for the FE modeling, since its APDL code can be connected with MATLAB script files.
In the third step, the DBNs are trained the generated samples as the input data and the FE results as the output data. The number of hidden layers should be optimized based on the production accuracy of the DBNs. If the accuracy is not sufficient, the number of training samples will be increased. This will improve the prediction accuracy, but will increase the computational cost.
In the fourth step, a large number of samples are generated by using the Monte Carlo simulation (MCS) according to the probability distribution of the random variables. Subsequently, the load effect of the critical components is predicted taking into consideration the samples of the random variables. Therefore, structural reliability can be estimated by counting the number of failure samples.
In the fifth step, the load-carrying capability of the residual structure is checked with consideration of the potential failure of one of the critical components. The system failure criterion is defined as when the structural system has lost its load-carrying capacity, or when the structural deflection is so large enough that it cannot carry the external load. Such critical states are extremely dangerous for in-service bridges, since the bridge is possible to collapse at any time.
In the sixth step, the components that are likely to fail are deleted or replaced based on the criterion of plastic or brittle failure. For example, cable failure is brittle, and the corresponding updating approach is to delete the failed cables. On the other hand, the bending failure of a girder is plastic, and thus the corresponding updating approach is to add a plastic hinge at the failure point. Subsequently, the structural system is updated, and the process goes back to the second step.
Finally, based on the above steps, all failure sequences are screened out. The correlation coefficient matrix can be estimated for the system reliability evaluation, and the system event tree can be estimated considering the critical failure sequences. Eventually, the system reliability can be assessed based on the series-parallel connection criterion.
4. Case Study of a Suspension Bridge
4.1. Prototype Suspension Bridge Characteristics
The prototype bridge used in the present study is the Nanxi Yangtze River Bridge, which is a suspension bridge in the Luyu highway of China. The general arrangement and the serial number of each component are shown in
Figure 6. The main cable is composed of parallel wire strands, the stiffening girders are steel box girders, and the pylon has a concrete box section. According to the construction sequence, the suspenders, main cables, and girders are divided into 64 segments. Details of the bridge characteristics can be found in Lu et al. [
1].
The random variables of the suspension bridge were the elasticity modulus of the concrete, steel box girder, and cable; the cross-sectional areas of the main girder, pylon, and cable; the unit weights of steel and concrete; and the moving vehicle load. Statistical distributions of the above parameters are given in
Table 1 and are based on the cable-stayed bridge in Lu et al. [
24].
The upper and lower bound of the random variables were determined according to the 3
σ principle [
30], where the range of the samples is defined as “
μ ± 3
σ”. Due to that the wires have manufacturing errors and the number of wires is huge, in the present study, the coefficient of variation for cables and suspenders was assumed as 0.05.
4.2. Traffic Load Modeling Using Weigh-in-Motion Data
Weigh-in-motion (WIM) data collected from a highway bridge in China were selected for the probabilistic modeling of traffic loads. The WIM system is illustrated in
Figure 7, and more detailed information can be found in Lu et al. [
11]. The traffic parameters utilized in the present study were vehicle weights, axle weights, driving lanes, and vehicle spacing. The proportion of trucks and the ratio of truck overloading were assumed as 12% and 21%, respectively.
All vehicles were classified into six types, based on vehicle configurations and axle characteristics. The occupancy of the different vehicle types and their distributions in each lane are displayed in
Figure 8.
As it can be seen in
Figure 8a, light cars (V1) had the highest proportion of 35.64%, while two-axle (V2) and six-axle trucks (V6) had a higher proportion compared to other types of trucks. As it can be observed in
Figure 8b, there is a high probability that the light cars are driving in the fast lane, while heavy trucks more mostly likely to drive in the slow lane. This phenomenon is in accordance with the practical situation.
The V6 trucks were selected as a prototype to investigate the probability of vehicle loads.
Figure 9a shows histograms of the distribution of monitored empirical gross vehicle weights (GVWs). They were fitted using a Gaussian mixture model (GMM) that is a combination of several Gaussian distributions. The two peaks of the PDF in
Figure 9a represent the normal and overloaded trucks. It can be observed that the GVW followed a multimodal distribution, where the truck overloading effect was captured by the probabilistic model.
Figure 9b plots the relationship between axle weight (AW) and GVW, where the star, circle, and plus symbols represent the weights of the first, second, and fifth axle, respectively, and
AWij represents the AW of the
j-th axle of the
i-th type of vehicles. It is observed that the AW and GVW were almost linearly dependent.
Based on the estimated probability distribution models, a stochastic traffic load model was developed using MCS.
Figure 10 plots the simulated dense traffic flow model, where each vehicle was simulated as a label. The stochastic traffic load model contained the parameters of vehicle type, GVW, driving lane, and vehicle spacing.
4.3. Reliability Analysis Based on the DBN Approach
Since structural failure is determined by the internal force of the bridge components, an appropriate FE model is essential to conduct the structural mechanical analysis. An FE of the bridge was developed in the commercial software ANSYS (
Figure 11).
Girders and pylons were modeled with Beam188 elements, and the main cables and girders were modeled with Lin180 elements. The main cable and the pylon were connected with coupled degrees of freedom, while the bottom parts were fully constrained. The initial tension forces in the main cable and suspenders were introduced as pre-stressed tension.
Using the FE model, the structural response can be approximated through the DBN. Initially, 20 uniformly-distributed training samples were selected as the input data. Corresponding structural response including cable forces and bending moments, were computed using the FE model and were used as the output data. The DBN was treated following the framework presented in
Figure 3. As an example, the approximated cable force is presented in
Figure 12. It can be observed that the response surface has captured the nonlinearity of structural mechanics.
In general, there are three types of failure modes related to bridge components: suspender fracture, main cable rupture, and girder bending moment failure. The reliability index of each bridge component is shown in
Figure 13.
It can be observed that the suspender had a relatively lower reliability index, while that of the main cable was the highest. Therefore, the first component more likely to fail should be the suspender. The most important step for the system reliability evaluation is to develop the fault tree of the structural system. The fault tree is composed of failure sequences, which can be evaluated based on the β-bound method.
4.4. System Reliability Evaluation
According to the first failure searching procedure, the components that are likely to fail were H20 and H44. Since the failure mode of suspenders is brittle failure, the failed cable should be directly removed from the structural system. Subsequently, the structural system was updated and approximated by the DBNs. At the second stage, the components that were likely to fail were stiffening girders associated with G19, G21, G43, and G45. Since the failure mode of girders is plastic failure, a hinge should be added at the node between two girder elements. At the third failure stage, the main cable is the most likely to fail component associated with C18, C20, C22, C42, C44, and C46. Consequently, the three-level fault tree of the bridge structural system is constructed, as shown in
Figure 14, where
βc is the conditional reliability index of the component corresponding to each failure stage.
Subsequently, the system reliability of the bridge was estimated through the series-parallel connection criterion. The corresponding series-parallel relationship is presented in
Figure 15. It can be observed that the system reliability index was 9.63 by considering the three-level fault tree.
According to the standard for structural reliability of highway engineering in China, a bridge controlled by brittle failure should have a reliability index greater than 5.2. Therefore, the bridge used in the present study has an adequate safety factor during its lifetime. Nevertheless, fatigue-corrosion damage is the main factor that leads to of bridge safety degradation. In future studies, th cable strength should be determined and updated in the system reliability evaluation.
Cable damage due to the fatigue-corrosion effect is a critical factor that may lead to service life-shortening, safety threats, and increasing life-cycle cost. Thus, the effect of cable damage on the system reliability of suspension bridge deserves investigation. According to the cable strength degradation model shown in
Figure 3, the system reliability of a bridge is adaptively updated.
Figure 16 shows the degradation curve of the system reliability of the suspension bridge based on the service lifetime of the suspender.
It can be observed that cable degradation due to fatigue damage leads to a slow reduction of the system reliability index from 9.63 to 7.22. However, the combined effect of fatigue-corrosion damage leads to a relatively sharp decrease of the system reliability index from 9.63 to 4.41. Therefore, cable corrosion is the dominant factor affecting bridge safety. In addition, given the target reliability index of 5.2, the suspension bridge has sufficient safety capacity accounting for cable fatigue damage. On the other hand, considering the fatigue-corrosion damage of suspenders, the maintenance time for cable replacement should be 26 years.
5. Conclusions
An intelligent approach for structural system reliability evaluation based on DBNs was presented. Mathematical models for system reliability evaluation of cable-stayed bridges were derived takinginto consideration the structural nonlinearity and high-order statically indeterminate characteristics. The theoretical basis for utilizing DBNs to approximate the structural load was introduced. A computational framework was presented to illustrate the procedures followed to evaluate the bridge system reliability via DBNs. The feasibility of the computational framework was demonstrated in a case study using a prototype suspension bridge. The conclusions are summarized as follows:
- (1)
DBN provides an accurate approximation of the load effect of a cable-supported bridge, taking into consideration the structural nonlinearity in the behavior of the system. Thus, DBN can be utilized as a meta-model for a large number of simulations as a result of reliability analysis.
- (2)
The suspender was found to have a relatively low reliability index, while that of the main cable was the highest. Therefore, the first component that was mostly likely to fail should be the suspender. In addition, the degradation of suspenders due to fatigue-corrosion damage was found to have a significant effect on the system reliability of the cable-supported bridge.
- (3)
The system reliability index of the prototype suspension bridge is 9.63 considering the three-level system fault tree. Accounting for fatigue damage and fatigue-corrosion damage, the system reliability index decrease to 7.22 and 4.41, respectively.
- (4)
According to the design standard for the structural reliability of highway bridges in China, a bridge controlled by brittle failure should have a reliability index greater than 5.2. Cable replacement should take place after 26 years of service time accounting for the fatigue-corrosion damage of suspenders.
Since fatigue-corrosion damage is the main factor leading to bridge safety degradation, further studies should focus on the cable degradation behavior and its effect on the failure sequence of bridges. In addition, more site-specific inspection data of the cables are needed to formulate the cable replacement scheme. Wind loads, seismic loads, and temperature loads should be taken into consideration in the reliability model, since these parameters are more critical for flexible bridges. Finally, structural health monitoring data can be employed to update the reliability model.