# BDVCM about the Steel Box Girder Reliability Prediction Based on Monitoring Extreme Stresses

Authors: Yuefei Liu; Xueping Fan; Guanghong Yang; Zhipeng Shang; Xiaoxiong Zhao
Journal: Advances in Civil Engineering
 (2021)

---

## Abstract

To reasonably predict the steel box girder reliability considering the dynamic dependence among the performance functions corresponding to the failure modes of the multiple monitoring points, this paper firstly adopts the dynamic monitoring extreme stresses of the multiple control points to build the Bayesian Dynamic Vine Copula Model (BDVCM) taking into account the dynamic dependence of the multiple monitoring variables through combining the vine copula technique with Bayesian Dynamic Linear Models (BDLM); secondly, with first-order second-moment method and the built BDVCM, the steel box girder reliability, taking into account dynamic dependence among the performance functions corresponding to the failure modes of the multiple monitoring points, is predicted; finally, the monitoring data from the five sections of an existing steel box girder were provided to illustrate the proposed model and approach. The analytical results illustrated that the predicted results, without considering the dynamic nonlinear dependence among the failure modes of the multiple monitoring points, are conservative.


---

## Body

## 1. Introduction

Structural Health Monitoring (SHM) systems have accumulated a large amount of monitored data in the long-term service periods. How to reasonably predict bridge system reliability with these data has become one of the key scientific problems in SHM field. Because of the same input loads, the different monitoring data of the multiple monitoring points from bridge system show the dependence, randomness, and so on. Further, dynamic dependence among the performance functions corresponding to the multiple failure modes of the monitoring points can be shown. Therefore, based on the dependent monitoring data, the dependence models of the multiple failure modes can be presented; further, it is more reasonable to dynamically predict the bridge girder reliability.Until now, some studies about the dynamic reliability prediction of bridge systems, taking into account the dynamic dependence among the multiple failure modes, have been done based on copula functions; for example, a copula-based method is presented to investigate the impact of copulas characterizing the dependence structures among the random variables for modeling bivariate distributions, respectively, on time-independent and time-dependent system reliability under incomplete probability information [1, 2]; a new model named BDGCM for characterizing the dynamic dependence between only two monitoring points is firstly built [3], and, based on it, the bridge system reliability is dynamically predicted, where correlation coefficients are not accurately computed and the built BDLM does not take the nonlinear dependence between two variables into account; copula-functions-based dynamic dependence model between only two monitoring points of arch bridge girder is given [4], where the correlation coefficients are approximately obtained and the nonlinear correlation between two monitoring variables is not taken into account in the built BDLM; Gaussian copula-Bayesian Dynamic Linear Model-based dynamic reliability prediction for Yitong River Bridge (a flying swallow profiled concrete-filled steel tube arch bridge) taking into account the time-variant dependence between the performance functions of only one pair of failure modes is made [5, 6], where the correlation coefficients are not accurately solved and the built BDLM also does not consider the nonlinear dependence between a pair of variables. The above research studies show that the existing dynamic reliability prediction methods of bridge structures considered time-variant nonlinear dependence between one pair of performance functions, and the correlation coefficients are not accurately solved. Meanwhile, actual bridge systems commonly have multiple control monitoring points. Therefore, building more accurate dynamic nonlinear dependence models among the multiple failure modes corresponding to the control monitoring points and further predicting the dynamic reliability of bridge systems should be further studied and become the aim of the present research.In view of the above problems, this study adopts the steel box girder as the research object, takes a newly built Bayesian Dynamic Vine Copula Model (BDVCM) to characterize the dynamic nonlinear correlations among the performance functions of the failure modes at multiple control monitoring points of the steel box girder based on BHM data, and further dynamically predicts the steel box girder reliability. In this paper, monitoring extreme data is “the everyday extreme stress data.” For each section, the control monitoring point is “the monitoring point where the absolute value of the monitoring stress is maximal”; each failure mode refers to the fact that “the monitoring extreme stress is bigger than the allowable stress.” The detailed contents are as follows: firstly, based on monitoring extreme stress data, the BDVCM about the multiple monitoring variables is built, which can characterize the dynamic dependence among the monitoring variables; secondly, with first-order second-moment (FOSM) method and the built BDVCM, the dynamic reliability prediction method of the steel box girder, considering dynamic nonlinear correlations among the performance functions for multiple failure modes corresponding to the multiple control monitoring points, is given; thirdly, the monitoring extreme stress data of a steel box girder are provided to illustrate the feasibility and application of the proposed model and method; finally, some valuable conclusions are summarized.


## 2. BDVCM of Multiple Monitoring Variables

This section firstly adopts the dynamic monitoring extreme stress data to establish the corresponding Bayesian Dynamic Linear Models (BDLM) about the multiple monitoring variables of the steel box girder sections, and the BDLM can dynamically predict the extreme data of monitoring variables; then, the vine copula model, characterizing the nonlinear correlation among multiple monitoring variable about the steel box girder sections, is built; finally, the Bayesian Dynamic Vine Copula Model (BDVCM) about multiple monitoring variables is given through combining the built BDLM with vine copula. The detailed contents are as follows.


### 2.1. BDLM of Monitoring Variables and the Corresponding Recursive Processes

Bayesian Dynamic Linear Models (BDLM) include Dynamic Linear Models (DLM) and the corresponding recursive processes based on Bayes method. The detailed contents are as follows.


#### 2.1.1. Construction of DLM

The DLM includes the monitoring equation, state equation, and initial state information, where the state equation shows changes of the system with time and reflects inner dynamic changes of the system and random disturbances; the monitoring equation expresses the relationship between the monitoring data and the current state parameters of the system; the initial state information is the probability distribution information about the initial state.Suppose that all the dynamic state variablesθt,i,i∈1,2,…,n follow normal distribution, are all Markov chains, and are internally independent and mutually independent between each other. The state error ωi,t and the monitoring error vi,t all follow normal distribution and are internally independent and mutually independent between each other. DLM is defined in the following form.Monitoring equation:(1)yi,t=θi,t+vi,t,vi,t∼Ν0,Vi,t.State equation:(2)θi,t=θi,t−1+ωi,t,ωi,t∼Ν0,Wi,t.Initial state information:(3)θi,t−1|Di,t−1∼Νmi,t−1,Ci,t−1,where yi,t,vi,t, and ωi,t are, respectively, the monitoring extreme stress data, monitoring error, and state error about the DLM of the ith monitoring variable at time t; Ν⋅ is normal distribution function; Vi,t is a variance; Wi,t is an evolution variance indicating the model recurrence uncertainty from time t-1 to time t; Di,t−1 is the information set of the ith monitoring variable at and before time t-1, and Di,t=yi,t,Di,t−1 is the information set of the DLM about the ith monitoring variable at and before time t, which includes mi,t (mean) and Ci,t (variance). It is assumed that error sequences vi,t and ωi,t are internally independent, mutually independent, and independent of θi,t.


#### 2.1.2. Bayesian Probability Recursive Processes of DLM

In this section, the Bayesian probability recursive processes of DLM are provided, which are as follows:(1)
The a posteriori probability density function (PDF) of the state at timet–1:(4)θi,t−1|Di,t−1∼Νmi,t−1,Ci,t−1.(2)
The a priori distribution of the state at timet:(5)θi,t|Di,t−1∼Νai,t,Ri,t.whereai,t=mi,t−1 and Ri,t=Ci,t−1+Wi,t.(3)
One-step prediction distribution of monitoring variable at timet:(6)yi,t|Di,t−1∼Νμi,t,σi,t2.whereμi,t=Εyi,t|Di,t−1=ai,t, Qi,t=σi,t2=varyi,t|Di,t−1=Ri,t+Vi,t, and Qi,t−1 is the prediction precision of BDLM.According to the definition of Highest Posterior Density (HPD) region [7], the one-step prediction interval of the monitoring variable with a 95% guarantee rate at time t is(7)μi,t−1.645σi,t2,μi,t+1.645σi,t2,whereμi,t−1.645σi,t2 is the predicted lower limit data and μi,t+1.645σi,t2 is the predicted upper limit data.(4)
The a posteriori distribution of the state at timet:(8)θi,t|Di,t∼Νmi,t,Ci,t,where mi,t=Εθi,t|Di,t=ai,t+Ai,tei,t; Ci,t=varθi,t|Di,t=Ri,t−Ai,tAi,tTσi,t2; Ai,t=Ri,tσi,t2−1; ei,t=yi,t−fi,t; Ai,t is the adaptive coefficient.


#### 2.1.3. Determination of the Main Parameters about BDLM

In this paper, the updating time interval of BDLM is one day. The main probability distribution parameters about BDLM includeVi,t, Wi,t,mi,t−1, and Ci,t−1.The monitoring extreme data at and before timet are smoothly processed with cubical smoothing algorithm with five-point approximation [3, 5–7], and the processed data can be approximately taken as the initial state data [3, 7]. Through carrying out parameter estimation for the initial state data, mi,t−1 and Ci,t−1 can be obtained. The variance Vi,t of the monitoring errors can be approximately solved with the difference between monitoring extreme data and the initial state data; the discount factor δi is considered as an aid to choosing Wi,t. In the proposed BDLM, the predicting behavior is usually rapidly achieved. It is convenient and natural to adopt a constant rate of increase of uncertainty or decay of information. Thus, for a given discount factor, equation (9) is chosen to solve the variance Wi,t of the state error.(9)Wi,t+1=−Ci,t+Ci,tδ,where δi is the discount factor and δi=0.98 according to the research experience of the authors [3, 7].


### 2.2. Vine Copula Model of Five-Dimensional Monitoring Variables

In this study, the bridge girder has multiple monitoring points which are corresponding to multiple output variables. Because of the same input vehicle loads, environmental loads, and so on, there exists nonlinear correlation between these output variables. In this section, the vine copula model considering nonlinear correlation of the five output variables is built with pair-copula theory and bivariate Gaussian copula theory.


#### 2.2.1. Pair-Copula Theory

Pair-copula construction module based joint probability distribution function of multivariate random variables is proposed by Bedford [8, 9]. The pair-copula construction module only considers the nonlinear correlation of two random variables.Assume that there is a five-dimensional random variableX=x1,x2,…,x5, and, according to the conditional probability density function theory [10], copula model based joint probability density function fx1,x2,…,x5 can be written as(10)fx1,x2,…,x5=cF1x1,F2x2,…,F5x5∏i=15fixi,where c is copula density function and Fixi and fixi are, respectively, cumulative distribution function (CDF) and probability density function (PDF) of xi.With equation (10), the bivariate joint PDF can be obtained with(11)fxi,xj=cijFixi,Fjxjfixifjxj,and, further, conditional PDF can be obtained; namely,(12)fxi|xj=fxi,xjf2xj=cijFixi,Fjxjfixi,and, with equation (12), conditional PDF about each variable x of five-dimensional random variables can be reached; namely,(13)fx|u=cxui|u−iFx|u−i,Fui|u−ifx|u−i,where ui represents the ith variable of five-dimensional random variables and u−i represents the four-dimensional random variables obtained through removing the ith variable from five-dimensional random variables.


#### 2.2.2. Bivariate Gaussian Copula Model

This section adopts the bivariate Gaussian copula model [6, 11] to build the nonlinear correlation model between two random variables. The CDF and PDF about the bivariate Gaussian copula model are, respectively,(14)Cu1,u2;ρ=ΦGΦ−1u1,Φ−1u2;ρ=∫−∞Φ−1u1∫−∞Φ−1u212π1−ρ21/2exp−r2−2ρrs+s221−ρ2drds,cu1,u2;ρ=11−ρ2expΦ−1u12+Φ−1u22−2ρΦ−1u1Φ−1u221−ρ2exp−Φ−1u12Φ−1u222.where ui=Fixi,i=1,2,…,5, Fixi is the marginal CDF of xi; Φ is standard normal distribution function; and ρ restricted to the interval −1,1 is the relevant parameter of Gaussian copula function.In this paper, the relevant parameterρ of bivariate Gaussian copula function is determined by the Pearson linear correlation coefficient [4–6].Suppose that Pearson linear correlation coefficient of random variablesxi and xj is ρi,j. The relevant parameter ρ of bivariate Gaussian copula function can be obtained with the following equation:(15)ρ=ρΦ−1ui,Φ−1uj=ρΦ−1Fixi,Φ−1Fjxj=ρxi,xj=ρi,j.


#### 2.2.3. Vine Copula Model of Five-Dimensional Random Variables

This section adopts multiple pair-copula models to construct the multivariate joint PDF. From the definition of vine structures [8, 9], as can be known that each vine consists of multiple trees, each tree contains multiple nodes, the line connecting two nodes is termed as edge, and different edges are independent of each other. Different vine structures have different properties, and the widely used ones are regular vine structures [8, 9], which include many types. The widely applied regular vine structure types are C-vine and D-vine. This section takes five-dimensional random variables as the example to build C-vine structure and D-vine structure, respectively, as shown in Figures 1 and 2. Figure 1 shows that C-vine structure of five-dimensional random variables contains four trees, and every tree has one main node that is connected to the other nodes. Each line (edge) connecting main node to the other node is corresponding to one pair-copula model, and different lines (edges) are independent of each other. Figure 2 shows that D-vine structure of five-dimensional random variables also has four trees and ten edges (lines), and each edge is corresponding to one pair-copula model. Through Figure 1, C-vine is composed of ten pair-copula models, and any two pair-copula models are independent. Meanwhile, in Figure 2, D-vine is also composed of ten pair-copula models, and any two pair-copula models are also independent.Figure 1
Five-dimensional C-vine structure.Figure 2
Five-dimensional D-vine structure.


### 2.3. BDVCM of Five-Dimensional Monitoring Variables

In this section, the Bayesian Dynamic Vine Copula Model (BDVCM) is a fusion of BDLM and vine copula models. The multivariate joint distribution function is decomposed into a series of BDVCM about any two monitoring variables; further, the BDVCM of multiple monitoring variables can be obtained through connecting the BDVCM of any two monitoring variables with vine structure.With equation (6) and [4–6], the dynamic Pearson linear correlation coefficient between any two monitoring variables (xi, xj) can be expressed with(16)ρti,j≈covDi,t−1,Dj,t−1σi,tσj,t,and, with equations (15) and (16), the dynamic relevant parameter ρt of bivariate Gaussian copula function can be obtained with(17)ρt=ρt1,2.The BDVCM of bivariate monitoring variables can be obtained through combining the built BDLM about two monitoring variables and Gaussian copula function of bivariate monitoring variables. Further, the BDVCM of multiple monitoring variables can be obtained through connecting the BDVCM of any two monitoring variables with the vine structure. The detailed modeling processes are as follows.With equation (6), it is known that two one-step prediction variables (y1, y2) both follow normal distribution; namely,(18)y1,t+1∼Nμ1,t+1,σ1,t+12,y2,t+1∼Nμ2,t+1,σ2,t+12,F1y1,t+1=Φy1,t+1−μ1,t+1σ1,t+1,F2y2,t+1=Φy2,t+1−μ2,t+1σ2,t+1.Let(19)u1=F1y1,t+1,u2=F2y2,t+1.Further, BDVCM about bivariate monitoring variables can be obtained with(20)y1,t+1,y2,t+1∼CF1y1,t+1,F2y2,t+1=Cu1,u2;ρt+1=ΦGF1y1,t+1,F2y2,t+1;ρt+1,where ρt is the dynamic relevant parameter of bivariate Gaussian copula function, which can be computed with equation (17). ΦG⋅ is Gaussian copula function.With equation (20) and vine structures shown in Figures 1 and 2, the BDVCM of five-dimensional monitoring variables can be reached.


## 2.1. BDLM of Monitoring Variables and the Corresponding Recursive Processes

Bayesian Dynamic Linear Models (BDLM) include Dynamic Linear Models (DLM) and the corresponding recursive processes based on Bayes method. The detailed contents are as follows.


### 2.1.1. Construction of DLM

The DLM includes the monitoring equation, state equation, and initial state information, where the state equation shows changes of the system with time and reflects inner dynamic changes of the system and random disturbances; the monitoring equation expresses the relationship between the monitoring data and the current state parameters of the system; the initial state information is the probability distribution information about the initial state.Suppose that all the dynamic state variablesθt,i,i∈1,2,…,n follow normal distribution, are all Markov chains, and are internally independent and mutually independent between each other. The state error ωi,t and the monitoring error vi,t all follow normal distribution and are internally independent and mutually independent between each other. DLM is defined in the following form.Monitoring equation:(1)yi,t=θi,t+vi,t,vi,t∼Ν0,Vi,t.State equation:(2)θi,t=θi,t−1+ωi,t,ωi,t∼Ν0,Wi,t.Initial state information:(3)θi,t−1|Di,t−1∼Νmi,t−1,Ci,t−1,where yi,t,vi,t, and ωi,t are, respectively, the monitoring extreme stress data, monitoring error, and state error about the DLM of the ith monitoring variable at time t; Ν⋅ is normal distribution function; Vi,t is a variance; Wi,t is an evolution variance indicating the model recurrence uncertainty from time t-1 to time t; Di,t−1 is the information set of the ith monitoring variable at and before time t-1, and Di,t=yi,t,Di,t−1 is the information set of the DLM about the ith monitoring variable at and before time t, which includes mi,t (mean) and Ci,t (variance). It is assumed that error sequences vi,t and ωi,t are internally independent, mutually independent, and independent of θi,t.


### 2.1.2. Bayesian Probability Recursive Processes of DLM

In this section, the Bayesian probability recursive processes of DLM are provided, which are as follows:(1)
The a posteriori probability density function (PDF) of the state at timet–1:(4)θi,t−1|Di,t−1∼Νmi,t−1,Ci,t−1.(2)
The a priori distribution of the state at timet:(5)θi,t|Di,t−1∼Νai,t,Ri,t.whereai,t=mi,t−1 and Ri,t=Ci,t−1+Wi,t.(3)
One-step prediction distribution of monitoring variable at timet:(6)yi,t|Di,t−1∼Νμi,t,σi,t2.whereμi,t=Εyi,t|Di,t−1=ai,t, Qi,t=σi,t2=varyi,t|Di,t−1=Ri,t+Vi,t, and Qi,t−1 is the prediction precision of BDLM.According to the definition of Highest Posterior Density (HPD) region [7], the one-step prediction interval of the monitoring variable with a 95% guarantee rate at time t is(7)μi,t−1.645σi,t2,μi,t+1.645σi,t2,whereμi,t−1.645σi,t2 is the predicted lower limit data and μi,t+1.645σi,t2 is the predicted upper limit data.(4)
The a posteriori distribution of the state at timet:(8)θi,t|Di,t∼Νmi,t,Ci,t,where mi,t=Εθi,t|Di,t=ai,t+Ai,tei,t; Ci,t=varθi,t|Di,t=Ri,t−Ai,tAi,tTσi,t2; Ai,t=Ri,tσi,t2−1; ei,t=yi,t−fi,t; Ai,t is the adaptive coefficient.


### 2.1.3. Determination of the Main Parameters about BDLM

In this paper, the updating time interval of BDLM is one day. The main probability distribution parameters about BDLM includeVi,t, Wi,t,mi,t−1, and Ci,t−1.The monitoring extreme data at and before timet are smoothly processed with cubical smoothing algorithm with five-point approximation [3, 5–7], and the processed data can be approximately taken as the initial state data [3, 7]. Through carrying out parameter estimation for the initial state data, mi,t−1 and Ci,t−1 can be obtained. The variance Vi,t of the monitoring errors can be approximately solved with the difference between monitoring extreme data and the initial state data; the discount factor δi is considered as an aid to choosing Wi,t. In the proposed BDLM, the predicting behavior is usually rapidly achieved. It is convenient and natural to adopt a constant rate of increase of uncertainty or decay of information. Thus, for a given discount factor, equation (9) is chosen to solve the variance Wi,t of the state error.(9)Wi,t+1=−Ci,t+Ci,tδ,where δi is the discount factor and δi=0.98 according to the research experience of the authors [3, 7].


## 2.1.1. Construction of DLM

The DLM includes the monitoring equation, state equation, and initial state information, where the state equation shows changes of the system with time and reflects inner dynamic changes of the system and random disturbances; the monitoring equation expresses the relationship between the monitoring data and the current state parameters of the system; the initial state information is the probability distribution information about the initial state.Suppose that all the dynamic state variablesθt,i,i∈1,2,…,n follow normal distribution, are all Markov chains, and are internally independent and mutually independent between each other. The state error ωi,t and the monitoring error vi,t all follow normal distribution and are internally independent and mutually independent between each other. DLM is defined in the following form.Monitoring equation:(1)yi,t=θi,t+vi,t,vi,t∼Ν0,Vi,t.State equation:(2)θi,t=θi,t−1+ωi,t,ωi,t∼Ν0,Wi,t.Initial state information:(3)θi,t−1|Di,t−1∼Νmi,t−1,Ci,t−1,where yi,t,vi,t, and ωi,t are, respectively, the monitoring extreme stress data, monitoring error, and state error about the DLM of the ith monitoring variable at time t; Ν⋅ is normal distribution function; Vi,t is a variance; Wi,t is an evolution variance indicating the model recurrence uncertainty from time t-1 to time t; Di,t−1 is the information set of the ith monitoring variable at and before time t-1, and Di,t=yi,t,Di,t−1 is the information set of the DLM about the ith monitoring variable at and before time t, which includes mi,t (mean) and Ci,t (variance). It is assumed that error sequences vi,t and ωi,t are internally independent, mutually independent, and independent of θi,t.


## 2.1.2. Bayesian Probability Recursive Processes of DLM

In this section, the Bayesian probability recursive processes of DLM are provided, which are as follows:(1)
The a posteriori probability density function (PDF) of the state at timet–1:(4)θi,t−1|Di,t−1∼Νmi,t−1,Ci,t−1.(2)
The a priori distribution of the state at timet:(5)θi,t|Di,t−1∼Νai,t,Ri,t.whereai,t=mi,t−1 and Ri,t=Ci,t−1+Wi,t.(3)
One-step prediction distribution of monitoring variable at timet:(6)yi,t|Di,t−1∼Νμi,t,σi,t2.whereμi,t=Εyi,t|Di,t−1=ai,t, Qi,t=σi,t2=varyi,t|Di,t−1=Ri,t+Vi,t, and Qi,t−1 is the prediction precision of BDLM.According to the definition of Highest Posterior Density (HPD) region [7], the one-step prediction interval of the monitoring variable with a 95% guarantee rate at time t is(7)μi,t−1.645σi,t2,μi,t+1.645σi,t2,whereμi,t−1.645σi,t2 is the predicted lower limit data and μi,t+1.645σi,t2 is the predicted upper limit data.(4)
The a posteriori distribution of the state at timet:(8)θi,t|Di,t∼Νmi,t,Ci,t,where mi,t=Εθi,t|Di,t=ai,t+Ai,tei,t; Ci,t=varθi,t|Di,t=Ri,t−Ai,tAi,tTσi,t2; Ai,t=Ri,tσi,t2−1; ei,t=yi,t−fi,t; Ai,t is the adaptive coefficient.


## 2.1.3. Determination of the Main Parameters about BDLM

In this paper, the updating time interval of BDLM is one day. The main probability distribution parameters about BDLM includeVi,t, Wi,t,mi,t−1, and Ci,t−1.The monitoring extreme data at and before timet are smoothly processed with cubical smoothing algorithm with five-point approximation [3, 5–7], and the processed data can be approximately taken as the initial state data [3, 7]. Through carrying out parameter estimation for the initial state data, mi,t−1 and Ci,t−1 can be obtained. The variance Vi,t of the monitoring errors can be approximately solved with the difference between monitoring extreme data and the initial state data; the discount factor δi is considered as an aid to choosing Wi,t. In the proposed BDLM, the predicting behavior is usually rapidly achieved. It is convenient and natural to adopt a constant rate of increase of uncertainty or decay of information. Thus, for a given discount factor, equation (9) is chosen to solve the variance Wi,t of the state error.(9)Wi,t+1=−Ci,t+Ci,tδ,where δi is the discount factor and δi=0.98 according to the research experience of the authors [3, 7].


## 2.2. Vine Copula Model of Five-Dimensional Monitoring Variables

In this study, the bridge girder has multiple monitoring points which are corresponding to multiple output variables. Because of the same input vehicle loads, environmental loads, and so on, there exists nonlinear correlation between these output variables. In this section, the vine copula model considering nonlinear correlation of the five output variables is built with pair-copula theory and bivariate Gaussian copula theory.


### 2.2.1. Pair-Copula Theory

Pair-copula construction module based joint probability distribution function of multivariate random variables is proposed by Bedford [8, 9]. The pair-copula construction module only considers the nonlinear correlation of two random variables.Assume that there is a five-dimensional random variableX=x1,x2,…,x5, and, according to the conditional probability density function theory [10], copula model based joint probability density function fx1,x2,…,x5 can be written as(10)fx1,x2,…,x5=cF1x1,F2x2,…,F5x5∏i=15fixi,where c is copula density function and Fixi and fixi are, respectively, cumulative distribution function (CDF) and probability density function (PDF) of xi.With equation (10), the bivariate joint PDF can be obtained with(11)fxi,xj=cijFixi,Fjxjfixifjxj,and, further, conditional PDF can be obtained; namely,(12)fxi|xj=fxi,xjf2xj=cijFixi,Fjxjfixi,and, with equation (12), conditional PDF about each variable x of five-dimensional random variables can be reached; namely,(13)fx|u=cxui|u−iFx|u−i,Fui|u−ifx|u−i,where ui represents the ith variable of five-dimensional random variables and u−i represents the four-dimensional random variables obtained through removing the ith variable from five-dimensional random variables.


### 2.2.2. Bivariate Gaussian Copula Model

This section adopts the bivariate Gaussian copula model [6, 11] to build the nonlinear correlation model between two random variables. The CDF and PDF about the bivariate Gaussian copula model are, respectively,(14)Cu1,u2;ρ=ΦGΦ−1u1,Φ−1u2;ρ=∫−∞Φ−1u1∫−∞Φ−1u212π1−ρ21/2exp−r2−2ρrs+s221−ρ2drds,cu1,u2;ρ=11−ρ2expΦ−1u12+Φ−1u22−2ρΦ−1u1Φ−1u221−ρ2exp−Φ−1u12Φ−1u222.where ui=Fixi,i=1,2,…,5, Fixi is the marginal CDF of xi; Φ is standard normal distribution function; and ρ restricted to the interval −1,1 is the relevant parameter of Gaussian copula function.In this paper, the relevant parameterρ of bivariate Gaussian copula function is determined by the Pearson linear correlation coefficient [4–6].Suppose that Pearson linear correlation coefficient of random variablesxi and xj is ρi,j. The relevant parameter ρ of bivariate Gaussian copula function can be obtained with the following equation:(15)ρ=ρΦ−1ui,Φ−1uj=ρΦ−1Fixi,Φ−1Fjxj=ρxi,xj=ρi,j.


### 2.2.3. Vine Copula Model of Five-Dimensional Random Variables

This section adopts multiple pair-copula models to construct the multivariate joint PDF. From the definition of vine structures [8, 9], as can be known that each vine consists of multiple trees, each tree contains multiple nodes, the line connecting two nodes is termed as edge, and different edges are independent of each other. Different vine structures have different properties, and the widely used ones are regular vine structures [8, 9], which include many types. The widely applied regular vine structure types are C-vine and D-vine. This section takes five-dimensional random variables as the example to build C-vine structure and D-vine structure, respectively, as shown in Figures 1 and 2. Figure 1 shows that C-vine structure of five-dimensional random variables contains four trees, and every tree has one main node that is connected to the other nodes. Each line (edge) connecting main node to the other node is corresponding to one pair-copula model, and different lines (edges) are independent of each other. Figure 2 shows that D-vine structure of five-dimensional random variables also has four trees and ten edges (lines), and each edge is corresponding to one pair-copula model. Through Figure 1, C-vine is composed of ten pair-copula models, and any two pair-copula models are independent. Meanwhile, in Figure 2, D-vine is also composed of ten pair-copula models, and any two pair-copula models are also independent.Figure 1
Five-dimensional C-vine structure.Figure 2
Five-dimensional D-vine structure.


## 2.2.1. Pair-Copula Theory

Pair-copula construction module based joint probability distribution function of multivariate random variables is proposed by Bedford [8, 9]. The pair-copula construction module only considers the nonlinear correlation of two random variables.Assume that there is a five-dimensional random variableX=x1,x2,…,x5, and, according to the conditional probability density function theory [10], copula model based joint probability density function fx1,x2,…,x5 can be written as(10)fx1,x2,…,x5=cF1x1,F2x2,…,F5x5∏i=15fixi,where c is copula density function and Fixi and fixi are, respectively, cumulative distribution function (CDF) and probability density function (PDF) of xi.With equation (10), the bivariate joint PDF can be obtained with(11)fxi,xj=cijFixi,Fjxjfixifjxj,and, further, conditional PDF can be obtained; namely,(12)fxi|xj=fxi,xjf2xj=cijFixi,Fjxjfixi,and, with equation (12), conditional PDF about each variable x of five-dimensional random variables can be reached; namely,(13)fx|u=cxui|u−iFx|u−i,Fui|u−ifx|u−i,where ui represents the ith variable of five-dimensional random variables and u−i represents the four-dimensional random variables obtained through removing the ith variable from five-dimensional random variables.


## 2.2.2. Bivariate Gaussian Copula Model

This section adopts the bivariate Gaussian copula model [6, 11] to build the nonlinear correlation model between two random variables. The CDF and PDF about the bivariate Gaussian copula model are, respectively,(14)Cu1,u2;ρ=ΦGΦ−1u1,Φ−1u2;ρ=∫−∞Φ−1u1∫−∞Φ−1u212π1−ρ21/2exp−r2−2ρrs+s221−ρ2drds,cu1,u2;ρ=11−ρ2expΦ−1u12+Φ−1u22−2ρΦ−1u1Φ−1u221−ρ2exp−Φ−1u12Φ−1u222.where ui=Fixi,i=1,2,…,5, Fixi is the marginal CDF of xi; Φ is standard normal distribution function; and ρ restricted to the interval −1,1 is the relevant parameter of Gaussian copula function.In this paper, the relevant parameterρ of bivariate Gaussian copula function is determined by the Pearson linear correlation coefficient [4–6].Suppose that Pearson linear correlation coefficient of random variablesxi and xj is ρi,j. The relevant parameter ρ of bivariate Gaussian copula function can be obtained with the following equation:(15)ρ=ρΦ−1ui,Φ−1uj=ρΦ−1Fixi,Φ−1Fjxj=ρxi,xj=ρi,j.


## 2.2.3. Vine Copula Model of Five-Dimensional Random Variables

This section adopts multiple pair-copula models to construct the multivariate joint PDF. From the definition of vine structures [8, 9], as can be known that each vine consists of multiple trees, each tree contains multiple nodes, the line connecting two nodes is termed as edge, and different edges are independent of each other. Different vine structures have different properties, and the widely used ones are regular vine structures [8, 9], which include many types. The widely applied regular vine structure types are C-vine and D-vine. This section takes five-dimensional random variables as the example to build C-vine structure and D-vine structure, respectively, as shown in Figures 1 and 2. Figure 1 shows that C-vine structure of five-dimensional random variables contains four trees, and every tree has one main node that is connected to the other nodes. Each line (edge) connecting main node to the other node is corresponding to one pair-copula model, and different lines (edges) are independent of each other. Figure 2 shows that D-vine structure of five-dimensional random variables also has four trees and ten edges (lines), and each edge is corresponding to one pair-copula model. Through Figure 1, C-vine is composed of ten pair-copula models, and any two pair-copula models are independent. Meanwhile, in Figure 2, D-vine is also composed of ten pair-copula models, and any two pair-copula models are also independent.Figure 1
Five-dimensional C-vine structure.Figure 2
Five-dimensional D-vine structure.


## 2.3. BDVCM of Five-Dimensional Monitoring Variables

In this section, the Bayesian Dynamic Vine Copula Model (BDVCM) is a fusion of BDLM and vine copula models. The multivariate joint distribution function is decomposed into a series of BDVCM about any two monitoring variables; further, the BDVCM of multiple monitoring variables can be obtained through connecting the BDVCM of any two monitoring variables with vine structure.With equation (6) and [4–6], the dynamic Pearson linear correlation coefficient between any two monitoring variables (xi, xj) can be expressed with(16)ρti,j≈covDi,t−1,Dj,t−1σi,tσj,t,and, with equations (15) and (16), the dynamic relevant parameter ρt of bivariate Gaussian copula function can be obtained with(17)ρt=ρt1,2.The BDVCM of bivariate monitoring variables can be obtained through combining the built BDLM about two monitoring variables and Gaussian copula function of bivariate monitoring variables. Further, the BDVCM of multiple monitoring variables can be obtained through connecting the BDVCM of any two monitoring variables with the vine structure. The detailed modeling processes are as follows.With equation (6), it is known that two one-step prediction variables (y1, y2) both follow normal distribution; namely,(18)y1,t+1∼Nμ1,t+1,σ1,t+12,y2,t+1∼Nμ2,t+1,σ2,t+12,F1y1,t+1=Φy1,t+1−μ1,t+1σ1,t+1,F2y2,t+1=Φy2,t+1−μ2,t+1σ2,t+1.Let(19)u1=F1y1,t+1,u2=F2y2,t+1.Further, BDVCM about bivariate monitoring variables can be obtained with(20)y1,t+1,y2,t+1∼CF1y1,t+1,F2y2,t+1=Cu1,u2;ρt+1=ΦGF1y1,t+1,F2y2,t+1;ρt+1,where ρt is the dynamic relevant parameter of bivariate Gaussian copula function, which can be computed with equation (17). ΦG⋅ is Gaussian copula function.With equation (20) and vine structures shown in Figures 1 and 2, the BDVCM of five-dimensional monitoring variables can be reached.


## 3. Dynamic Reliability Prediction of the Steel Box Girder considering Nonlinear Correlation among Performance Functions about Failure Modes for Five Cross Sections Based on SHM Data

With equations (16)-(20) and Section 2.3, the BDVCM, considering nonlinear dependence among the five monitoring points of five cross sections, is provided. Further, bridge girder system reliability can be dynamically predicted. Reliability indices and failure probability about multiple monitoring points can be solved with the FOSM method (see Section 3.1). Failure probability of bridge girder system can be computed through Section 3.2.


### 3.1. Reliability Prediction Formula Based on the FOSM Method

In this study, the adopted bridge girder has five monitoring cross sections, and each cross section has a control monitoring point [12]. The five predicted performance functions about the five control monitoring points are(21)hiyi,t=σ−γPyi,t,i=1,2,…,5,where hi is the ith performance function, σ is the strength (allowable stress), yi,t is the monitoring or predicted everyday extreme stress of the ith control monitoring point at time t, and γP=1.15 is the correction factor of the sensors [7].Based on the FOSM method [13], the predicted reliability index βi,t and the corresponding failure probability pfi,t of the ith control monitoring point can be solved with(22)βi,t=μσ−γPμyi,tσσ2+γPσyi,t2,(23)pfi,t=Φ−βi,t.


### 3.2. The Steel Box Girder Reliability Prediction Based on BDVCM

This section supposes that any system composed of two control monitoring points is parallel. With equation (21), the failure probability of parallel system composed of two control monitoring points can be computed with the following equation:(24)gi,t+1=σ−γPyi,t+1,i=1,2,Pg1,t+1≤0,g2,t+1≤0;=PG1g1,t+1≤G10,G2g2,t+1≤G20=PF1,t+1≤G10,F2,t+1≤G20=CG10,G20;ρ12,t+1=Cpf1,t+1,pf2,t+1;ρt+1=ΦGpf1,t+1,pf2,t+1;ρt+1.where Fi,t+1=Φσ−γPyi,t+1−γPμi,t+1/γPσi,t+1,Gi0=Φ−μi,t+1/σi,t+1,i=1,2, pf1,t+1 and pf2,t+1 can be computed with equation (23). C is copula function, and ρt+1 is the time-variant relevant parameter obtained with equation (17). ΦG⋅ is Gaussian copula function.Because the multiple parallel systems composed of any two control monitoring points are serial, with equations (21)-(24), the predicted failure probability of the steel box girder considering the nonlinear correlation among performance functions of failure modes about multiple control monitoring points can be solved with the following equation:(25)Pfsystem,t+1=maxCpf1,t+1,pf2,t+1;ρ12,t+1,Cpf1,t+1,pf3,t+1;ρ13,t+1,Cpf1,t+1,pf4,t+1;ρ14,t+1,Cpf1,t+1,pf5,t+1;ρ15,t+1,Cpf2,t+1,pf3,t+1;ρ23,t+1,Cpf2,t+1,pf4,t+1;ρ24,t+1,Cpf2,t+1,pf5,t+1;ρ25,t+1,Cpf3,t+1,pf4,t+1;ρ34,t+1,Cpf3,t+1,pf5,t+1;ρ35,t+1,Cpf4,t+1,pf5,t+1;ρ45,t+1.


## 3.1. Reliability Prediction Formula Based on the FOSM Method

In this study, the adopted bridge girder has five monitoring cross sections, and each cross section has a control monitoring point [12]. The five predicted performance functions about the five control monitoring points are(21)hiyi,t=σ−γPyi,t,i=1,2,…,5,where hi is the ith performance function, σ is the strength (allowable stress), yi,t is the monitoring or predicted everyday extreme stress of the ith control monitoring point at time t, and γP=1.15 is the correction factor of the sensors [7].Based on the FOSM method [13], the predicted reliability index βi,t and the corresponding failure probability pfi,t of the ith control monitoring point can be solved with(22)βi,t=μσ−γPμyi,tσσ2+γPσyi,t2,(23)pfi,t=Φ−βi,t.


## 3.2. The Steel Box Girder Reliability Prediction Based on BDVCM

This section supposes that any system composed of two control monitoring points is parallel. With equation (21), the failure probability of parallel system composed of two control monitoring points can be computed with the following equation:(24)gi,t+1=σ−γPyi,t+1,i=1,2,Pg1,t+1≤0,g2,t+1≤0;=PG1g1,t+1≤G10,G2g2,t+1≤G20=PF1,t+1≤G10,F2,t+1≤G20=CG10,G20;ρ12,t+1=Cpf1,t+1,pf2,t+1;ρt+1=ΦGpf1,t+1,pf2,t+1;ρt+1.where Fi,t+1=Φσ−γPyi,t+1−γPμi,t+1/γPσi,t+1,Gi0=Φ−μi,t+1/σi,t+1,i=1,2, pf1,t+1 and pf2,t+1 can be computed with equation (23). C is copula function, and ρt+1 is the time-variant relevant parameter obtained with equation (17). ΦG⋅ is Gaussian copula function.Because the multiple parallel systems composed of any two control monitoring points are serial, with equations (21)-(24), the predicted failure probability of the steel box girder considering the nonlinear correlation among performance functions of failure modes about multiple control monitoring points can be solved with the following equation:(25)Pfsystem,t+1=maxCpf1,t+1,pf2,t+1;ρ12,t+1,Cpf1,t+1,pf3,t+1;ρ13,t+1,Cpf1,t+1,pf4,t+1;ρ14,t+1,Cpf1,t+1,pf5,t+1;ρ15,t+1,Cpf2,t+1,pf3,t+1;ρ23,t+1,Cpf2,t+1,pf4,t+1;ρ24,t+1,Cpf2,t+1,pf5,t+1;ρ25,t+1,Cpf3,t+1,pf4,t+1;ρ34,t+1,Cpf3,t+1,pf5,t+1;ρ35,t+1,Cpf4,t+1,pf5,t+1;ρ45,t+1.


## 4. Application to an Existing Bridge

This paper takes Tianjin Fumin Bridge from China as the analytical example [7, 14, 15]. The adopted girder has five monitoring cross sections which are sections A, B, C, D, and E shown in Figure 3. Each cross section has installed three stress sensors at three different monitoring points which are all Fiber Bragg Grating sensors [12]. The longitudinal stresses about these five sections are monitored through an array of stress sensors [12]. These stress sensors about the five monitoring cross sections are shown in Table 1 and Figure 4, and these cross sections are serial. The steel type of the steel box girder is Q345qD, and the mean value and coefficient of variance are, respectively, 411 MPa and 0.098 [12].Figure 3
(a) Fumin Bridge. (b) Five monitoring sections of Fumin Bridge girder system.
(a)(b)Table 1
Stress sensors of the five monitoring cross sections [12].
SectionsDifferent stress sensors at different monitoring pointsAFBG01074, FBG01081, FBG01078BFBG01947, FBG01949, FBG01046CFBG01059, FBG01071, FBG01073DFBG01012, FBG01015, FBG01005EFBG01029, FBG01021, FBG01028Figure 4
Stress sensors of the five monitoring cross sections. (a) Section A. (b) Section B. (c) Section C. (d) Section D. (e) Section E.
(a)(b)(c)(d)(e)


### 4.1. Monitoring Everyday Extreme Stress Data of the Five Monitoring Cross Sections

Sampling rate about the stress data is 2.5 Hz sampling frequency. The everyday extreme stresses of the five cross sections are, respectively, monitored for two hundred days. The monitoring data are shown in Figure5, which show that each section is monitored by three sensors at three different locations and the monitoring data are dynamic and random [12]. According to measurement program [7, 15], these data can reflect the safety-based dynamic reliability of girder system. The maximum absolute value of everyday extreme stresses at the three points of each section is defined as the everyday extreme stress of each section [12].Figure 5
Monitoring extreme stresses of the five sections. (a) Section A. (b) Section B. (c) Section C. (d) Section D. (e) Section E.
(a)(b)(c)(d)(e)


### 4.2. BDLM of Monitoring Everyday Extreme Stress Data of Each Monitoring Section

Figure5 shows that the sensors FBG01074, FBG01947, FBG01071, FBG01015, and FBG01028, respectively, monitored the everyday extreme stress about the control locations of sections A, B, C, D, and E. So, the monitoring data of the five sensors are applied to build the BDLM for predicting the corresponding everyday extreme stresses. The detailed contents are as follows.Firstly, the monitoring extreme data of the five sensors are, respectively, smoothly processed with cubical smoothing algorithm with five-point approximation [7], and the processed data can be approximately considered as the initial state information [12]. Through Kolmogorov-Smirnov (K-S) test, these initial state data about the five sensors are all approximately simulated with normal distribution. Secondly, based on equations (1)-(9) of Section 2.1 including Sections 2.1.1-2.1.3, the BDLM can be easily built. Finally, with the built BDLM, the everyday extreme data can be dynamically predicted. The predicted results are shown in Figure 6. As can be seen, the predicted extreme stress data and the predicted extreme stress ranges all fit the changing rules of monitoring extreme stress data.Figure 6
Predicted extreme stresses of the five sections. (a) Section A. (b) Section B. (c) Section C. (d) Section D. (e) Section E.
(a)(b)(c)(d)(e)


### 4.3. Dynamic Reliability Prediction of the Bridge Girder System

#### 4.3.1. Bridge Girder System Reliability Prediction without considering the Dependence among the Multiple Failure Modes

In this section, dynamic reliability index data at the control points in Tianjin Fumin Bridge girder system are predicted with equation (22) and the built BDLM. The allowable stress of bridge steel follows normal distribution with 411 MPa (mean value) and 0.098 (coefficient of variance).With equation (22), the computation formula of reliability indices can be obtained with(26)βi,P,t=411−1.15μi,t411×0.0982+1.15σi,t2,i=1,2,…,5,(27)pfi,P,t=Φ−βi,P,t,i=1,2,…,5,where βi,P,t is the predicted reliability index at the ith monitoring point at time t; μi,t,σi,t can be obtained with equation (6).With equation (26), the reliability indices of the five control monitored points can be dynamically predicted, and the predicted results are shown in Figure 7, from which, as can be seen, the predicted reliability indices and the predicted ranges all fit the changing rules of monitored reliability indices.Figure 7
Predicted reliability indices of the five sections. (a) Section A. (b) Section B. (c) Section C. (d) Section D. (e) Section E.
(a)(b)(c)(d)(e)The comparative analysis is made, and the minimum value is selected as bridge girder reliability index. The results are shown in Figure8, from which, as can be noticed, the reliability indices of the control monitored point at section B are the bridge girder reliability indices.Figure 8
Comparative analysis of predicted reliability indices.Based on these solved reliability indices and equation (27), failure probability of the bridge girder system can be calculated. The comparative analysis about the failure probability of the five control points is made, and the maximum value shown in Figure 9 is selected as the time-variant failure probability of the steel box girder. From Figure 9, as can be seen, the failure probability of the control monitoring point at section B is the bridge girder failure probability.Figure 9
Comparative analysis of predicted failure probability.


#### 4.3.2. Bridge Girder Reliability Prediction considering the Dependence among the Multiple Failure Modes

In this section, dynamic reliability of steel box girder considering correlation is predicted. Based on C-vine structure shown in Figure1, with equation (16), the time-variant correlation coefficients (e.g., ρAB,ρAC,ρAD,ρAE,ρBC|A,ρBD|A,ρBE|A,ρCD|AB,ρCE|AB,ρDE|ABC) between the performance functions of any two points can be obtained, where the time-variant correlation coefficients are equal to the relevant parameters of the corresponding copula functions. ρAB,ρAC,ρAD,ρAE, and ρBC|A are shown in Figure 10(a), and ρBD|A,ρBE|A,ρCD|AB,ρCE|AB, and ρDE|ABC are displayed in Figure 10(b).Figure 10
Time-dependent nonlinear correlation coefficients.
(a)(b)Based on Figure10, bivariant Gaussian-copula PDF functions and contour maps considering nonlinear dependence between any pair of failure modes about any pair of control monitoring points are analyzed. Conditional on sections A and B, the time-variant correlation between sections C and D was analyzed, and the analytical results of the 100th day are selected and shown in Figure 11, from which, as can be known, the dependence between any pair of failure modes is nonlinear.Figure 11
The time-variant correlation between sections C and D conditional on sections A and B.With equation (24), pfAB,pfAC,pfAD,pfAE, and pfBC|A can be computed and shown in Figure 12(a), and pfBD|A,pfBE|A,pfCD|AB,pfCE|AB, and pfDE|ABC can be solved and shown in Figure 12(b). Further, with equation (24), girder system failure probability considering nonlinear dependence among the failure modes can be obtained and shown in Figure 13, from which, as can be known, the bridge girder system failure probability taking into account the nonlinear dependence among the performance functions of failure modes is smaller than the one obtained without taking into account the nonlinear dependence, which proves that it is more reasonable and necessary to take into account the nonlinear dependence among the performance functions of failure modes for dynamic reliability prediction of structural system.Figure 12
Failure probability of bivariate vine structure.
(a)(b)Figure 13
Predicted bridge girder failure probability.


## 4.1. Monitoring Everyday Extreme Stress Data of the Five Monitoring Cross Sections

Sampling rate about the stress data is 2.5 Hz sampling frequency. The everyday extreme stresses of the five cross sections are, respectively, monitored for two hundred days. The monitoring data are shown in Figure5, which show that each section is monitored by three sensors at three different locations and the monitoring data are dynamic and random [12]. According to measurement program [7, 15], these data can reflect the safety-based dynamic reliability of girder system. The maximum absolute value of everyday extreme stresses at the three points of each section is defined as the everyday extreme stress of each section [12].Figure 5
Monitoring extreme stresses of the five sections. (a) Section A. (b) Section B. (c) Section C. (d) Section D. (e) Section E.
(a)(b)(c)(d)(e)


## 4.2. BDLM of Monitoring Everyday Extreme Stress Data of Each Monitoring Section

Figure5 shows that the sensors FBG01074, FBG01947, FBG01071, FBG01015, and FBG01028, respectively, monitored the everyday extreme stress about the control locations of sections A, B, C, D, and E. So, the monitoring data of the five sensors are applied to build the BDLM for predicting the corresponding everyday extreme stresses. The detailed contents are as follows.Firstly, the monitoring extreme data of the five sensors are, respectively, smoothly processed with cubical smoothing algorithm with five-point approximation [7], and the processed data can be approximately considered as the initial state information [12]. Through Kolmogorov-Smirnov (K-S) test, these initial state data about the five sensors are all approximately simulated with normal distribution. Secondly, based on equations (1)-(9) of Section 2.1 including Sections 2.1.1-2.1.3, the BDLM can be easily built. Finally, with the built BDLM, the everyday extreme data can be dynamically predicted. The predicted results are shown in Figure 6. As can be seen, the predicted extreme stress data and the predicted extreme stress ranges all fit the changing rules of monitoring extreme stress data.Figure 6
Predicted extreme stresses of the five sections. (a) Section A. (b) Section B. (c) Section C. (d) Section D. (e) Section E.
(a)(b)(c)(d)(e)


## 4.3. Dynamic Reliability Prediction of the Bridge Girder System

### 4.3.1. Bridge Girder System Reliability Prediction without considering the Dependence among the Multiple Failure Modes

In this section, dynamic reliability index data at the control points in Tianjin Fumin Bridge girder system are predicted with equation (22) and the built BDLM. The allowable stress of bridge steel follows normal distribution with 411 MPa (mean value) and 0.098 (coefficient of variance).With equation (22), the computation formula of reliability indices can be obtained with(26)βi,P,t=411−1.15μi,t411×0.0982+1.15σi,t2,i=1,2,…,5,(27)pfi,P,t=Φ−βi,P,t,i=1,2,…,5,where βi,P,t is the predicted reliability index at the ith monitoring point at time t; μi,t,σi,t can be obtained with equation (6).With equation (26), the reliability indices of the five control monitored points can be dynamically predicted, and the predicted results are shown in Figure 7, from which, as can be seen, the predicted reliability indices and the predicted ranges all fit the changing rules of monitored reliability indices.Figure 7
Predicted reliability indices of the five sections. (a) Section A. (b) Section B. (c) Section C. (d) Section D. (e) Section E.
(a)(b)(c)(d)(e)The comparative analysis is made, and the minimum value is selected as bridge girder reliability index. The results are shown in Figure8, from which, as can be noticed, the reliability indices of the control monitored point at section B are the bridge girder reliability indices.Figure 8
Comparative analysis of predicted reliability indices.Based on these solved reliability indices and equation (27), failure probability of the bridge girder system can be calculated. The comparative analysis about the failure probability of the five control points is made, and the maximum value shown in Figure 9 is selected as the time-variant failure probability of the steel box girder. From Figure 9, as can be seen, the failure probability of the control monitoring point at section B is the bridge girder failure probability.Figure 9
Comparative analysis of predicted failure probability.


### 4.3.2. Bridge Girder Reliability Prediction considering the Dependence among the Multiple Failure Modes

In this section, dynamic reliability of steel box girder considering correlation is predicted. Based on C-vine structure shown in Figure1, with equation (16), the time-variant correlation coefficients (e.g., ρAB,ρAC,ρAD,ρAE,ρBC|A,ρBD|A,ρBE|A,ρCD|AB,ρCE|AB,ρDE|ABC) between the performance functions of any two points can be obtained, where the time-variant correlation coefficients are equal to the relevant parameters of the corresponding copula functions. ρAB,ρAC,ρAD,ρAE, and ρBC|A are shown in Figure 10(a), and ρBD|A,ρBE|A,ρCD|AB,ρCE|AB, and ρDE|ABC are displayed in Figure 10(b).Figure 10
Time-dependent nonlinear correlation coefficients.
(a)(b)Based on Figure10, bivariant Gaussian-copula PDF functions and contour maps considering nonlinear dependence between any pair of failure modes about any pair of control monitoring points are analyzed. Conditional on sections A and B, the time-variant correlation between sections C and D was analyzed, and the analytical results of the 100th day are selected and shown in Figure 11, from which, as can be known, the dependence between any pair of failure modes is nonlinear.Figure 11
The time-variant correlation between sections C and D conditional on sections A and B.With equation (24), pfAB,pfAC,pfAD,pfAE, and pfBC|A can be computed and shown in Figure 12(a), and pfBD|A,pfBE|A,pfCD|AB,pfCE|AB, and pfDE|ABC can be solved and shown in Figure 12(b). Further, with equation (24), girder system failure probability considering nonlinear dependence among the failure modes can be obtained and shown in Figure 13, from which, as can be known, the bridge girder system failure probability taking into account the nonlinear dependence among the performance functions of failure modes is smaller than the one obtained without taking into account the nonlinear dependence, which proves that it is more reasonable and necessary to take into account the nonlinear dependence among the performance functions of failure modes for dynamic reliability prediction of structural system.Figure 12
Failure probability of bivariate vine structure.
(a)(b)Figure 13
Predicted bridge girder failure probability.


## 4.3.1. Bridge Girder System Reliability Prediction without considering the Dependence among the Multiple Failure Modes

In this section, dynamic reliability index data at the control points in Tianjin Fumin Bridge girder system are predicted with equation (22) and the built BDLM. The allowable stress of bridge steel follows normal distribution with 411 MPa (mean value) and 0.098 (coefficient of variance).With equation (22), the computation formula of reliability indices can be obtained with(26)βi,P,t=411−1.15μi,t411×0.0982+1.15σi,t2,i=1,2,…,5,(27)pfi,P,t=Φ−βi,P,t,i=1,2,…,5,where βi,P,t is the predicted reliability index at the ith monitoring point at time t; μi,t,σi,t can be obtained with equation (6).With equation (26), the reliability indices of the five control monitored points can be dynamically predicted, and the predicted results are shown in Figure 7, from which, as can be seen, the predicted reliability indices and the predicted ranges all fit the changing rules of monitored reliability indices.Figure 7
Predicted reliability indices of the five sections. (a) Section A. (b) Section B. (c) Section C. (d) Section D. (e) Section E.
(a)(b)(c)(d)(e)The comparative analysis is made, and the minimum value is selected as bridge girder reliability index. The results are shown in Figure8, from which, as can be noticed, the reliability indices of the control monitored point at section B are the bridge girder reliability indices.Figure 8
Comparative analysis of predicted reliability indices.Based on these solved reliability indices and equation (27), failure probability of the bridge girder system can be calculated. The comparative analysis about the failure probability of the five control points is made, and the maximum value shown in Figure 9 is selected as the time-variant failure probability of the steel box girder. From Figure 9, as can be seen, the failure probability of the control monitoring point at section B is the bridge girder failure probability.Figure 9
Comparative analysis of predicted failure probability.


## 4.3.2. Bridge Girder Reliability Prediction considering the Dependence among the Multiple Failure Modes

In this section, dynamic reliability of steel box girder considering correlation is predicted. Based on C-vine structure shown in Figure1, with equation (16), the time-variant correlation coefficients (e.g., ρAB,ρAC,ρAD,ρAE,ρBC|A,ρBD|A,ρBE|A,ρCD|AB,ρCE|AB,ρDE|ABC) between the performance functions of any two points can be obtained, where the time-variant correlation coefficients are equal to the relevant parameters of the corresponding copula functions. ρAB,ρAC,ρAD,ρAE, and ρBC|A are shown in Figure 10(a), and ρBD|A,ρBE|A,ρCD|AB,ρCE|AB, and ρDE|ABC are displayed in Figure 10(b).Figure 10
Time-dependent nonlinear correlation coefficients.
(a)(b)Based on Figure10, bivariant Gaussian-copula PDF functions and contour maps considering nonlinear dependence between any pair of failure modes about any pair of control monitoring points are analyzed. Conditional on sections A and B, the time-variant correlation between sections C and D was analyzed, and the analytical results of the 100th day are selected and shown in Figure 11, from which, as can be known, the dependence between any pair of failure modes is nonlinear.Figure 11
The time-variant correlation between sections C and D conditional on sections A and B.With equation (24), pfAB,pfAC,pfAD,pfAE, and pfBC|A can be computed and shown in Figure 12(a), and pfBD|A,pfBE|A,pfCD|AB,pfCE|AB, and pfDE|ABC can be solved and shown in Figure 12(b). Further, with equation (24), girder system failure probability considering nonlinear dependence among the failure modes can be obtained and shown in Figure 13, from which, as can be known, the bridge girder system failure probability taking into account the nonlinear dependence among the performance functions of failure modes is smaller than the one obtained without taking into account the nonlinear dependence, which proves that it is more reasonable and necessary to take into account the nonlinear dependence among the performance functions of failure modes for dynamic reliability prediction of structural system.Figure 12
Failure probability of bivariate vine structure.
(a)(b)Figure 13
Predicted bridge girder failure probability.


## 5. Conclusions

This paper presented a new dynamic reliability prediction method for the steel box girder. By combining Bayesian Dynamic Linear Model (BDLM) with vine copula functions, the Bayesian Dynamic Vine Copula Model (BDVCM) is provided. The predicted results illustrate that the dynamic failure probability taking into account the nonlinear dependence is lower than the one without taking into account the nonlinear dependence. It is indicated that the computed results without taking into account the dynamic nonlinear dependence are conservative, and the predicted results taking into account the dynamic nonlinear dependence are more reasonable and necessary.



---

*Source: 8816459-2021-02-20.xml*