The actual structural systems have many failure modes. Due to the same random sources owned by the performance functions of these failure modes, there usually exist some nonlinear correlations between the various failure modes. How to handle the nonlinear correlations is one of the main scientific problems in the field of structural system reliability. In this paper, for the two-component systems and multiple-component systems with multiple failure modes, the mixed copula models for time-independent reliability analysis of series systems, parallel systems, series-parallel systems, and parallel-series systems are presented. These obtained mixed copula models, considering the nonlinear correlation between failure modes, are obtained with the chosen optimal copula functions with the Bayesian selection criteria and Monte Carlo Sampling (MCS) method. And a numerical example is provided to illustrate the feasibility and application of the built mixed models for structural system reliability.
1. Introduction
Today’s structural systems are becoming more complex and more sophisticated. Therefore, the evaluation of structural system reliability is becoming harder. This means that the derivations based on classical assumptions are no longer satisfactory for the analysis of systems in terms of reliability. The evaluation of today’s real life systems needs more detailed and complicated statistical analysis. Dependence between the components is one of the intractable realistic assumptions that need to be carefully considered [1–4].
Especially for the actual bridge system, there exist many failure modes such as flexural failure and shearing failure. As the limit state functions (performance function) of these failure modes may have the same random sources, it is not mutually exclusive yet among failure modes [5–7]. Naturally, how to model the correlation among failure modes in aspect of the reliability analysis of structural systems is one of the most significant topics. The classic Pearson correlation coefficient is mainly used for characterizing the correlation among failure modes of structural system before; however, it has a few disadvantages. The copulas, unlike the Pearson correlation coefficient only applied for describing linear correlation, offer a flexible tool for deriving nonlinear dependence, especially tail dependence among failure modes. The aim of this study is to introduce the copula function as a useful tool for modeling the dependence among failure modes of bridge system. More recently, the copula theory has been primarily used in mechanical engineering [8] and hydraulic engineering [9] but little used in bridge engineering.
In this paper, firstly, several commonly used elliptical copulas and Archimedean copulas were introduced, and their application in the correlation analysis was also described in detail. And then based on the introduced copulas, with the aid of copula Bayesian selection criteria [7, 10, 11], a flexible mixed copula model is constructed, by means of linear weighted model, to model the nonlinear dependence among failure modes. For the two-component systems and multiple-component systems with multiple failure modes, the performance function value of failure modes is chosen as the copula functions’ analytical variable to construct mixed copula model. With Monte Carlo Sampling method, the unknown parameters of copula functions can be approximately determined; therefore, the mixed copula function, considering the correlation between failure modes, is built. And then use the built mixed copula model to solve the structural system’s reliability. Finally, a numerical example is provided to illustrate the feasibility and application of the built mixed copula models.
2. Mixed Copula Models
The mixed copula models are commonly built through different combinations of several copula functions, which respectively possess various characteristics in aspect of the dependence modeling. For the complex and ever-changing correlations among random variables, they are not adequate for modeling the dependence by means of only a single copula function. Therefore, it is necessary to build a flexible mixed copula model to describe the complex dependence structures. Moreover, it is essential to select several appropriate copula functions from the existing copulas before building the mixed copula model. The Bayesian selection criteria are just chosen for constructing the mixed copula model, as this method is independent of the parameter estimation and may be applied to all known copula classes [7, 10].
2.1. Basic Theory about Copula Functions
Sklar’s theory [12] clearly indicated that, given random variables X=X1,X2,…,Xn with continuous marginal distribution functions FX1,FX2,…,FXn and n-dimensional joint distribution function F, there exists a unique copula function C such that, for all x=x1,x2,…,xn∈Rn, (1)Fx1,x2,…,xn=CFX1x1,FX2x2,…,FXnxn.
According to (1), the joint probability density function of the random variable X=X1,X2,…,Xn is given as follows:(2)px1,x2,…,xn=cFX1x1,FX2x2,…,FXnxn∏i=1npXixi,where pXixi is the marginal probability density function (PDF) and c(u) is the joint PDF of the copula function denoted as follows:(3)cu1,…,un=∂nC∂u1⋯∂unu1,…,un.
If FXi for i=1,2,…,n have generalized inverse functions FXi-1 for i=1,2,…,n, then, (1) can be rewritten as follows:(4)FFX1-1u1,FX2-1u2,…,FXn-1un=Cu1,u2,…,un,where ui=FXi(i=1,2,…,n), u1,u2,…,un=u.
Equation (4) reveals how to construct the copula function of a multivariate distribution with given marginal distributions. It follows from the probability integral transform that the random variables Ui≡FXiXi,i=1,2,…,n, are all uniformly distributed on (0, 1). Conversely, Xi=FXi-1Ui,i=1,2,…,n, U=U1,U2,…,Un; with this in mind, (4) is expressed as(5)PX1≤x1,X2≤x2,…,Xn≤xn=PX1≤FX1-1u1,X2≤FX2-1u2,…,Xn≤FXn-1un=PU1≤u1,U2≤u2,…,Un≤un=Cu1,u2,…,un.
According to (1)–(5), the copula function organically combines each marginal distribution function with the multivariate joint distribution functions; therefore, it not only considers the dependence among random variables but also simplifies the probability modeling process for multivariate random variables.
2.2. Bayesian Selection Criteria of Copula Functions
Huard et al. [10] proposed the Bayesian copula selection method applied to the bivariate copula functions. Bayesian copula selection method was deduced by the Bayesian hypothesis testing to choose the best copula. However, we apply this method to choose not just the best one but also the better copulas as the candidate copulas for constructing the mixed copula model.
For the Bayesian copula selection criteria, it is essential to propose the optical copula function selection criteria described by(6)Copt=argClmaxPHt∣D,I,where Copt is the optical copula function, Ht is denoted as the hypotheses that the right copula is copula Cl(l=1,2,…,L), as the candidate copulas, Cl is denoted by Cl=Clui,vi;θi(i=1,2,…,n), and n pair of quantiles ui,vi are supposed as n mutually independent pairs composing the data set D. Besides, I stands for some a priori information. According to Bayesian formula, it can be obtained that(7)PHt∣D,I=PD∣Ht,IPHt∣IPD∣I,where PD∣Ht,I is the likelihood, PHt∣I is the a priori copula family, and PD∣I is the normalization constant. The rank correlation coefficient Kendall’s τ, denoted by τ=gl(θ) to be the correlation parameters for all the selected copulas, is introduced in order to give the formula for PHt∣D,I:(8)PHt∣D,I=∫τ∈ΩPD∣Ht,τ,IPHt∣τ,IPτ∣IPD∣Idτ,where PD∣Ht,τ,I is the likelihood with respect to τ, because of the n mutually independent pairs, and is also described as the copula density:(9)PD∣Ht,τ,I=∏l=1nClui,vi;gl-1τ.
Substituting (9) into (8), we obtain(10)PHt∣D,I=1PD∣I∫τ∈Ω∏l=1nClui,vi;gl-1τPHt∣τ,IPτ∣Idτ,where P(D∣I) is the normalization constant and PHt∣τ,IP(τ∣I)=1/λ(Λ), where λ(Λ) is denoted as the Lebesgue measure of Λ=[-1,1] in this paper; therefore, (10) can be also expressed as follows:(11)PHt∣D,I∝1λΛ∫Ωl∩Λ∏l=1nClui,vi;gl-1τdτ.
Obviously, the optimal copula function selection criteria can be described as the computation of the highest weight Wl:(12)Wl=1λΛ∫Ωl∩Λ∏i=1nclui,vi;gl-1τdτ.
In especial, the “right” copulas include not just the best copula with the highest weight but also the better copula with higher weight for the mixed copula modeling in this paper.
Two particular classes of copulas that proved to be useful in dependence modeling are the elliptical and Archimedean classes [7, 13]. The four copulas (one of the elliptical copulas, Gaussian copula, and three of the Archimedean copulas, Gumbel copula, Clayton copula, and Frank copula) are employed as the candidates of the Bayesian copula selection criteria. Finally, the chosen ones are used for constructing the mixed copula model. Among the chosen copulas, Gaussian copulas and Frank copulas are able to capture symmetric dependence structures among random variables. Different from them, Gumbel copulas and Clayton copulas exhibit asymmetric dependence and Gumbel copulas are especially employed for describing upper tail dependence structures, while Clayton copulas are employed for that of the lower tail. And, then, the probability density function graphs and their contour plots are depicted as in Figures 1–4.
PDF and contour plots of Gaussian copula function.
PDF and contour plots of Gumbel copula function.
PDF and contour plots of Clayton copula function.
PDF and contour plots of Frank copula function.
2.3. Modeling Method of Mixed Copula Models
There exist complex nonlinear correlations among failure modes; therefore, just one copula is not enough to characterize dependence structures among failure modes. With the aid of a weighted combination of the selected copula functions, a mixed copula is approximately constructed in order to describe the complex dependence structures among failure modes [7, 14]; namely,(13)Cu,v;a,θ=∑i=1NaiCiu,v;θiN≤4,where Ci(i=1,2,…,N) are, respectively, denoted as one of the chosen copulas among the four copulas and a=a1,a2,…,aN is the weighting coefficient vector along with the inner parameter vector θ=θ1,θ2,…,θN of the corresponding copulas, and u,ν are variables of the chosen copula functions, and what is more, a1+a2+⋯+aN=1.
The copula function and Monte Carlo Simulation (MCS) method are contributed to the mixed copula model analysis. The specific steps can be stated as follows.
2.3.1. Monte Carlo Sampling (MCS)
According to the distribution types of random variables for the performance functions of failure modes, the samples for the random variables can be obtained with MCS method. After substituting the samples into the corresponding performance functions, we can get the random sequence gji,i=1,2,…,n,j=1,2,…,N, for each of the performance functions and obtain the sequences Fji=Fgji,i=1,2,…,n,j=1,2,…,N, of the corresponding empirical distribution function by mean of Matlab software, where g is the performance function and F is the empirical distribution function.
2.3.2. Copula Function Selection
With the scatter plot of samples, we can judge approximate distribution features of the samples. Furthermore, it is extremely vital to determine which classes of copula functions to choose as the candidate copulas with Bayesian copula selection criteria. In this paper, the copula classes we chose are Gaussian copulas, Gumbel copulas, and Clayton copulas as well as Frank copulas, whose parameter domain along with the definition and the domain of Kendall’s τ are listed in Table 1. The next is to build the mixed copula function through the selected copulas and (13).
Basic cases of the parameters of copula functions.
Copula
τ=g(θ)
τ∈Ω
θ∈Ω
Gaussian
2Xπ2πarcsinθ
-1,1
-1,1
Gumbel
1-θ-1
0,1
1,+∞
Clayton
1-22+θ
0,1
0,+∞
Frank
1-4θ1-1θ∫0θtet-1dt
-1,0∪0,1
-∞,0∪0,+∞
2.3.3. Parameter Estimation for the Mixed Copula
Applying the least residual error quadratic sum method OLS, one of the curve fitting criteria, we can obtain parameters of the mixed copula function. And the formula of the OLS is represented as follows:(14)FOLS=1n∑i=1nFF1i,F2i,…,FNi-CF1i,F2i,…,FNi2,where F(·) is the joint empirical distribution function of the empirical distribution function sequences for each performance function and C(F1)i,(F2)i,…,(FN)i is the mixed copula function value. The unknown parameter value for the mixed copula function can be obtained by the optimization computation with the rule of OLS, while the optimized parameters must be determined to make sure that the value of FOLS is the minimum.
2.3.4. Goodness-of-Fit Estimation for the Mixed Copula
Through MCS method, the values of joint empirical distribution functions and the corresponding mixed copula function are obtained, and, then, the scatter plot between them is drawn. Finally, the goodness of fit for the mixed copula can be determined through the scatter plot.
3. Mixed Copula Model Expressions about Joint Failure Probability of Structural System 3.1. Mixed Copula Model Expression about Joint Failure Probability of Two-Component Series System
For the two-component series system which is shown in Figure 5, suppose that the performance function of the component failure mode is(15)giX=giX1,X2,…,Xn,i=1,2.
Two-component series system.
With (5), as a check, the probability that both failure modes occur is denoted as(16)Pg1X≤0,g2X≤0=PFg1g1X≤Fg10,Fg2g2X≤Fg20=PU1≤Fg10,U2≤Fg20=CFg10,Fg20=Cpfg1,pfg2.
Fg(0,0)=CFg1(0),Fg2(0) is stated in the fundamental theorem of Sklar. Therefore, the failure probability of two-component series system (at least one failure mode of the two components occurred) can be solved as(17)pf=Pg1X≤0∪g2X≤0=Pg1X≤0+Pg2X≤0-Pg1X≤0,g2X≤0=pfg1+pfg2-Cpfg1,pfg2,where pfg1,pfg2, respectively, denote the failure probability of the two failure modes and C· is two-component series system copula function.
3.2. Mixed Copula Model Expression about Joint Failure Probability of Two-Component Parallel System
For the two-component parallel system which is shown in Figure 6, suppose that the performance function of the component failure mode is(18)giX=giX1,X2,…,Xn,i=1,2.
Two-component parallel system.
With (5) and (16), the failure probability of two-component parallel system (two failure modes of the two components meantime occurred) can be obtained:(19)pf=Pg1X≤0∩g2X≤0=Cpfg1,pfg2,where pfg1,pfg2, respectively, denote the failure probability of the two failure modes and C· is two-component parallel system copula function.
3.3. Mixed Copula Model Expression about Joint Failure Probability of Multiple-Component Series System
For the multiple-component series system which is shown in Figure 7, suppose that the performance function of the component failure mode is (20)giX=giX1,X2,…,Xn,i=1,2,…,N.
Multiple-component series system.
With (5), as a check, the probability that all the failure modes meantime occur is denoted as(21)Pg1X≤0,g2X≤0,…,gNX≤0=PFg1g1X≤Fg10,Fg2g2X≤Fg20,…,FgNgNX≤FgN0=PU1≤Fg10,U2≤Fg20,…,UN≤FgN0=CFg10,Fg20,…,FgN0=Cpfg1,pfg2,…,pfgN.
Fg(0,0,…,0)=CFg10,Fg20,…,FgN(0) is stated in the fundamental theorem of Sklar. Therefore, the failure probability of multiple-component series system (at least one failure mode of the multiple components occurred) can be solved as(22)pf=Pg1X≤0∪g2X≤0∪⋯∪gNX≤0=∑i=1NPgiX≤0-∑i1<i2Pgi1X≤0,gi2X≤0+⋯+-1n+1∑i1<i2<⋯<inPgi1X≤0,gi2X≤0,…,ginX≤0+⋯+-1N+1Pg1X≤0,g2X≤0,…,gNX≤0=∑i=1npfgi-∑i1<i2Cpfgi1,pfgi2+⋯+-1n+1∑i1<i2<⋯<inCpfgi1,pfgi2,…,pfgin+⋯+-1N+1Cpfg1,pfg2,…,pfgN,where pfg1,pfg2,…,pfgN, respectively, denote the failure probability of the multiple failure modes and C· is multiple-component series system copula function.
3.4. Mixed Copula Model Expression about Joint Failure Probability of Multiple-Component Parallel System
For the multiple-component parallel system which is shown in Figure 8, suppose that the performance function of the component failure mode is(23)giX=giX1,X2,…,Xn,i=1,2,…,N.
Multiple-component parallel system.
With (5), the failure probability of two-component parallel system (all the failure modes of all the components meantime occurred) can be obtained:(24)pf=Pg1X≤0,g2X≤0,…,gNX≤0=PFg1g1X≤Fg10,Fg2g2X≤Fg20,…,FgNgNX≤FgN0=PU1≤Fg10,U2≤Fg20,…,UN≤FgN0=CFg10,Fg20,…,FgN0=Cpfg1,pfg2,…,pfgN,where pfg1,pfg2,…,pfgN, respectively, denote the failure probability of the N failure modes and C· is multiple-component parallel system copula function.
3.5. Mixed Copula Model Expression about Joint Failure Probability of Series-Parallel System
For series-parallel system shown in Figure 9, in this paper, only the correlation between internal components of each subparallel system is considered, while the correlation between subparallel systems is not considered and considered to be mutually independent. Therefore, with (24), all the copula models of all the subparallel systems’ failure probability can be obtained, and then, with (25), the failure probability of the series-parallel system can be solved:(25)pfs-p=P⋃i=1NFi=P⋃i=1N⋂jNiFij=∑i=1NPFi-∑i1<i2PFi1PFi2+⋯+-1n+1∑i1<i2<⋯<in∏k=i1inPFk+⋯+-1N+1∏i=1NPFi,where N is the total number of the subparallel systems; Ni is the total number of the components in the ith subparallel system; Fi=⋂j=1NiFij is the failure probability of the ith subparallel system; Fij=gij(X)≤0 is the failure probability of the jth component in the ith subparallel system; PFi,i=1,2,…,N, can be solved with (24), which are the failure probability of the subparallel system considering correlation between internal components of each subparallel system.
Series-parallel system.
3.6. Mixed Copula Model Expression about Joint Failure Probability of Parallel-Series System
For parallel-series system shown in Figure 10, in this paper, only the correlation between internal components of each subseries system is considered, while the correlation between subseries systems is not considered and considered to be mutually independent. Therefore, with (22), all the copula models of all the subseries systems’ failure probability can be obtained, and then, with (26), the failure probability of the parallel-series system can be solved as(26)pfp-s=P⋂i=1NFi=∏i=1NPFi=∏i=1NP⋃j=1NiFij,where N is the total number of the subseries systems; Ni is the total number of the components in the ith subseries system; Fi=⋂j=1NiFij is the failure probability of the ith subseries system; Fij=gij(X)≤0 is the failure probability of the jth component in the ith subseries system; PFi,i=1,2,…,N, can be solved with (22), which are the failure probability of the subseries system considering correlation between internal components of each subseries system.
Parallel-series system.
In this paper, firstly with First-Order Reliability Method (FORM), the structural reliability index and the corresponding failure probability of each failure mode can be solved, and then with the constructed mix copula functions, such as (17), (19), (22), (24), (25), and (26), the failure probability of structural system considering correlation between failure modes can be obtained.
4. Numerical Example: System Reliability Analysis of Simply Supported Cored Slab Bridge
For the simply supported cored slab bridge shown in Figure 11, the total span is 13 m, the computed span is 12.6 m, the clear width of bridge deck is 7 m, the width of footway on both sides of bridge deck is 1 m, and the whole bridge is composed of nine concrete cored slabs [15]. The design reference period of this bridge is 100 years. And this bridge has been served for 32 years. At the 32nd year, the resistance of each girder follows normal distribution; the distribution parameters are, respectively, 796.04 kN·m (mean value) and 91.783 kN·m (standard deviation).
Cross section of bridge (unit: cm).
Based on Figure 11 and reference [16], the failure criterion of bridge system is as follows: if any two adjacent girders both failed, then the whole bridge system failed. According to the failure criterion, the bridge system is a series-parallel system, which is shown in Figure 12.
Series-parallel system.
The performance function of each girder’s failure mode is(27)Zi=gXi1,Xi2,Xi3=gRi,Gi,Qi=Ri-Gi-Qi,i=1,2,…,9,where Ri is the resistance of the ith girder, Gi is the dead load effect of the ith girder, and Qi is the live load effect of the ith girder. At the 32nd year, the distribution parameters about dead load effects are listed in Table 2, where the standard deviation is not changed, because, for load effect, the variation of variables is very small. And the distribution parameters about vehicle load effects are listed in Table 3, which occur due to symmetrical variable load and is not applicated to the unsymmetrical variable load.
Distribution parameters of beams’ dead load effects in 32nd year.
Girder number
1
2
3
4
5
6
7
8
9
Mean value/kN⋅m
212.86
221.79
223.91
230.74
260.69
230.74
223.91
221.79
212.86
Standard deviation/kN⋅m
9.68
9.68
9.68
9.68
9.68
9.68
9.68
9.68
9.68
Distribution parameters of beams’ maximum live load effects in 32nd year.
Girder number
1
2
3
4
5
6
7
8
9
Mean value/kN⋅m
85.42
96.42
116.8
144.4
180.9
144.4
116.8
96.42
85.42
Standard deviation/kN⋅m
12.71
14.35
17.38
21.49
26.92
21.49
17.38
14.35
12.71
Based on the reliability analysis method of two-component parallel system considering the correlation between failures modes described in Section 3.2, the reliability analysis processes of bridge system are as follows.
Because the bridge system is symmetrical, the subparallel systems (1#-2#, 2#-3#, 3#-4#, and 4#-5#) are used to analyze the reliability and failure probability of the bridge system.
Based on Tables 2 and 3 and (25), the scatter plot between two random sampling sequences of the corresponding limit state functions for the two failure modes of each subparallel system (1#-2#, 2#-3#, 3#-4#, and 4#-5#) can be obtained. Then, according to the characteristics of the obtained scatter plots, the candidate copula functions, which can approximately describe distribution features of the samples, can be selected, and then with Bayesian selection criteria described in Section 2.2, the suitable copula functions, from the selected candidate copula functions, can be obtained which can be used to build the mixed copula model. The parameters of the built mixed copula function for each subparallel system are listed in Table 4. Finally, the mixed copula modes, for each subparallel system, are built. Further, PDF plots and contour plots for the mixed copulas, the scatter plots between empirical distributions and mixed copula functions, and scatter plots for two limit state functions are presented, respectively, in Figures 13–16.
Parameters of each of the mixed copulas for four subsystems.
a
b
1-a-b
ρ
α
θ
1#-2#
0.7438
0.1958
0.0604
−0.1868
2.4872
4.4278
2#-3#
0.8521
0.1479
0
0.1014
2.4207
4.4078
3#-4#
0.9381
0.0619
0
0.1502
2.4611
4.3984
4#-5#
0.7690
0.2310
0
−0.3490
2.4814
4.4192
The subsystem 1#-2#.
PDF plot for the mixed copula
Contour plot
Empirical distribution versus mixed copula
Scatter plot for two limit state functions
The subsystem 2#-3#.
PDF plot for the mixed copula
Contour plot
Empirical distribution versus mixed copula
Scatter plot for two limit state functions
The subsystem 3#-4#.
PDF plot for the mixed copula
Contour plot
Empirical distribution versus mixed copula
Scatter plot for two limit state functions
The subsystem 4#-5#.
PDF plot for the mixed copula
Contour plot
Empirical distribution versus mixed copula
Scatter plot for two limit state functions
With FORM, the corresponding reliability indices to each girder’s failure mode are, respectively,(28)β1=5.3429,β2=5.1159,β3=4.8484,β4=4.4417,β5=3.6869.
Then, with the equation pf=Φ(-β), the corresponding failure probability to each girder’s failure mode is, respectively,(29)pfg1=4.57×10-8,pfg2=1.56×10-7,pfg3=6.22×10-7,pfg4=4.46×10-6,pfg5=1.14×10-4.
With (19) and Table 4, the following can be obtained:(30)pf12=Cpf1,pf2=8.61×10-11,pf23=Cpf2,pf3=3.12×10-10,pf34=Cpf3,pf4=1.32×10-9,pf45=Cpf4,pf5=1.31×10-7,where pf12 is the failure probability when girder 1# and girder 2# meantime failed, pf23 is the failure probability when girder 2# and girder 3# meantime failed, pf34 is the failure probability when girder 3# and girder 4# meantime failed, and pf45 is the failure probability when girder 4# and girder 5# meantime failed.
Suppose that failure modes between subsystems are mutually independent, and then the failure probability of structural system is approximately (31)pf=-12∑j=i+1,i=18pfij+-13∑t=k+1,j=i+1,k>i,i=17pfijpfkt+⋯+-19pf12pf23pf34pf45pf56pf67pf78pf89≈∑j=i+1,i=18pfij-∑t=k+1,j=i+1,k>i,i=17pfijpfkt=2pf12+pf23+pf34+pf45-4pf12pf23+pf12pf34+pf12pf45+pf23pf34+pf23pf45+pf34pf45-pf122-pf232-pf342-pf452=2.65×10-7,where, because the structural system is symmetric, the failure probability of the subsystems at symmetric positions is the same. Namely, pf12=pf89, pf23=pf78,…,pf45=pf56.
For each subparallel system, the failure probability without considering the correlation between the two components is(32)hf12=pf1pf2=7.13×10-15,hf23=pf2pf3=9.70×10-14,hf34=pf3pf4=2.77×10-12,hf45=pf4pf5=5.08×10-10.
Suppose that failure modes between subsystems are mutually independent, and the failure modes between the two components of each subsystem are also mutually independent; then, the failure probability of structural system [17] is approximately(33)pf=-12∑j=i+1,i=18hfij+-13∑t=k+1,j=i+1,k>i,i=17hfijhfkt+⋯+-19hf12hf23hf34hf45hf56hf67hf78hf89≈∑j=i+1,i=18hfij-∑t=k+1,j=i+1,k>i,i=17hfijhfkt=2hf12+hf23+hf34+hf45-4hf12hf23+hf12hf34+hf12hf45+hf23hf34+hf23hf45+hf34hf45-hf122-hf232-hf342-hf452=1.02×10-9.
From the above solved system’s failure probability shown in (31) and (33), it can be seen that the failure probability of the series-parallel system considering the correlation between two adjacent girders is larger than the failure probability without considering the correlation between two adjacent girders, which showed that the series-parallel system considering the correlation between two adjacent girders more easily failed. Further, it is illustrated that considering the correlation between two adjacent girders of each subparallel system is essential and applicable for solving the reliability of the series-parallel system.
5. Conclusions
For the two-component systems and multiple-component systems with multiple failure modes, this paper presents the mixed copula models for reliability analysis of series systems, parallel systems, series-parallel systems, and parallel-series systems. The mixed copula model is obtained with the chosen optimal copula functions with the Bayesian method. Through a numerical example, it is illustrated that the calculated failure probability when considering the correlation between failure modes is larger than that without considering the correlation between failure modes. It is verified that the solved failure probability is conservative without considering the correlation between failure modes.
This paper provided a new method for characterizing the correlation between failure modes and solving the reliability of the system considering correlation between failure modes.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the Fundamental Research Funds for the Central Universities (lzujbky-2015-300, lzujbky-2015-301).
EryilmazS.Multivariate copula based dynamic reliability modeling with application to weighted-k-out-of-n systems of dependent components201451232810.1016/j.strusafe.2014.05.0042-s2.0-84902962790BaluA. S.RaoB. N.Multicut-high dimensional model representation for structural reliability bounds estimation under mixed uncertainties201227641943810.1111/j.1467-8667.2012.00763.x2-s2.0-84862154355Dueñas-OsorioL.RojoJ.Reliability assessment of lifeline systems with radial topology201126211112810.1111/j.1467-8667.2010.00661.x2-s2.0-78651415479ValD. V.StewartM. G.MelchersR. E.Life-cycle performance of RC bridges: probabilistic approach2000151142510.1111/0885-9507.001672-s2.0-0034125478LiG.ChengG.-D.Correlation between structural failure modes and calculation of system reliability under hazard loads2001183192-s2.0-0035354243LiD.-Q.WuS.-B.System reliability analysis on main girder of plane gates considering multiple correlated failure modes20094078708772-s2.0-69249130088LiuY.LuD.Reliability analysis of two-dimensional series portal-framed bridge system based on mixed copula functions20145749510510.4028/www.scientific.net/KEM.574.952-s2.0-84884881954HanW. Q.ZhouJ. Y.SunK. Z.Copulas in the mechanical components reliability analysisReliability Technology Academic Exchange Meeting and the 4th National Machinery Industry Reliability Engineering Branch of the Second Plebary Committee Conference ProceedingsXiangtan, China2010 (Chinese)HouY. Y.SongS. B.Research of the joint distribution for flood peak and volume based on copula function201030113944HuardD.ÉvinG.FavreA.-C.Bayesian copula selection200651280982210.1016/j.csda.2005.08.0102-s2.0-33750412434FanN.HeX. L.ZhaoQ.Optimal copula function selection based on Bayesian methodology2012294516522MR3013379SklarA.Fonctions de répartition à n dimensions et leurs marges19598229231NelsenR. B.2006New York, NY, USASpringerMR2197664HuL.2002New Haven, Conn, USAYale UniversityFanX. P.2010Harbin, ChinaSchool of Civil Engineering, Harbin Institute of TechnologyGharaibehE. S.1999Denver, Colo, USAUniversity of ColoradoMelchersR. E.1987Chichester, UKEllis Horwood