# The Method of Solving Structural Reliability with Multiparameter Correlation Problem

Authors: Juan Du; Haibin Li; Yun He
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

Correlation among variables must be considered to accurately reflect the level of structure reliability. This problem has referential value to engineering practice and has attracted attention from relevant scholars and industries. In this paper, Copula function was used to build the joint probability density function among all variables. The key is to describe the correlation among variables, solve the correlation parameterθ of Copula function, and select the type of correlation structure among variables. The correlation parameter θ of Copula function was solved using Pearson linear correlation coefficient and maximum likelihood estimation. Based on the Akaike information criteria (AIC) and Bayesian information criteria (BIC), the optimal Copula function was selected, and the correlation structure among variables was determined. Monte Carlo method, which is based on Nataf inverse transformation, was introduced and used to evaluate the reliability of the correlated variable. Finally, this paper proposed the reliability calculation method based on dual neural network and direct integration by establishing the dual neural network of original and integrand functions. Compared with the Monte Carlo method, the proposed method can be utilized to efficiently and precisely calculate the structure reliability of multiple correlated random variables.


---

## Body

## 1. Introduction

In structural reliability analysis, variables are often correlated, for instance, the positive correlation between the seismic peak displacement and permanent displacement [1] and among the stresses of the weld fatigue damage [2] and the significant negative correlation between the shear strength parameters of rock and soil and the parameters of pile-load-displacement curve [3, 4]. Therefore, the correlation among variables must be considered to rationally analyze the structural reliability. The traditional calculation method for structural reliability includes first-order [5] and second-order moments [6], which are only limited to the linear correlation of variables. A 2D or multidimensional distribution model is often used to characterize the correlation among variables, such as 2D lognormal distribution [7] and multidimensional normal distribution [8]. However, one of the drawbacks of these models is the requirement of same-edge distribution variables, which greatly limit their application in reliability analysis.The emergence and development of Copula theory solve the above problems and provide a new way to construct joint probability density function of correlated variables [9]. Copula theory was first proposed by Sklar in 1959 [10]. Under the framework of Copula function, variables could follow arbitrary edge distributions, and the linear and nonlinear correlations among variables could be defined [11]. Copula functions were introduced in computational structural reliability because of its capacity to handle the arbitrary correlations of variables. Li et al. used the 2D Frank Copula function to construct the joint probability distribution function of two seismic attenuation models [12]. Based on the Copula function, Tang et al. presented the load-displacement hyperbolic probabilistic analysis of pile foundation and used the Copula function to establish the joint probability distribution function of hyperbolic parameters [13]. Huang et al. established a Copula function model of rock mass shear strength parameters based on measured data [14]. Liu et al. used Copulas to develop a reliability model for systems with s-dependent degradation processes. The proposed model accommodates the assumptions of s-dependence among the degradation processes and allows for different marginal distributions [15]. For two-component and multiple-component systems with multiple failure modes, Liu and Fan established the mixed Copula models for time-independent reliability analysis of series systems, parallel systems, series-parallel systems, and parallel-series systems [16]. Based on the Copula function, Xu et al. constructed a joint probability density function, which was integrated on the failure field to calculate the failure probability of the structure [17]. In reliability-based design optimization (RBDO) of the structure, the exact joint probability density function of the input relevant variables is necessary to obtain the optimal design [18]. Lee et al. used the Copula function to consider the correlation among variables and solve the problem of RBDO [19, 20] and found that the correlation significantly affected RBDO results [21–25]. They also analyzed the influence of the confidence level for Copula function on RBDO results [26] and used the Bayesian and Markov chain Monte Carlo methods to select the Copula function [27]. For RBDO with varying standard deviations (STDs), Cho et al. used the Copula function for the design of sensitivity and then improved the efficiency of the calculation [28].The calculation methods for structural reliability are classified in three categories. The first method solves the structural response or the probability feature quantity of performance function [29–31]; however, the computational accuracy in this method is heavily dependent on the form of performance functions. The second method was the Monte Carlo method, which directly applies sampling statistics [32–34]. The disadvantage of this method is the large computational complexity. The third method was direct integration [35], which has encountered mathematical difficulties in the multiple integral and quantitative description of computational accuracy.Neural networks can approximate any functions and thus could be properly applied to structural reliability. Papadrakakis et al. combined neural network with Monte Carlo simulation to analyze the reliability of elastoplastic structure; the time-consuming calculation of the Monte Carlo simulation was reduced, and the efficiency of the calculation was improved [36]. Lopes et al. analyzed reliability using neural network instead of finite element analysis. The neural network had advantages in computational efficiency when compared with Monte Carlo method [37]. Zuo et al. used neural network to fit the performance function of the structure. The values of performance function and partial derivatives were obtained at the point of mean values. Hence, the moments of performance function were calculated based on the moments of random variables [38]. Cheng and Li used the neural network to simulate the limit state equation of long-span bridge. The genetic algorithm was used to train the network, and the failure probability of the structure was obtained [39]. Meng et al. utilized BP neural network for nonlinear mapping function trains to obtain the explicit expression of stress in response to the random variable. A study analyzed reliability and the sensitivity of metal structure by combining random perturbation theory and first-order second-moment method [40]. Elhewy et al. proposed a response surface method based on artificial neural networks and consequently reduced the computational complexity for reliability analysis [41]. Li et al. proposed a structural reliability method, which is based on neural network and direct integration, to solve the reliability of independent variables [42]. Dai and Cao developed a new neural network model based on wavelet support vector machine for reliability analysis; this model has extended the application of wavelet neural network to a high dimension [43].This paper proposed a direct integration method that is based on dual neural network and could be applied for reliability calculation. The proposed method is composed of two neural networks with similar structure, multiple inputs, single output, and single hidden. By designing the function relation between the weights of two neural networks, one neural network is able to approximate integrands, whereas the other approximates original function. Therefore, the above networks were called integrand and original function neural networks. We only need to train the integrand neural network. Thus, the weights of original function neural network were provided directly by the function relation between the weights of the two neural networks. We subsequently used the original function neural network to calculate the multiple integral. In the proposed method, the integrand could easily obtain the sample data that were directly trained, and the integral computational accuracy would be greatly improved. Therefore, the proposed method can efficiently and accurately solve structural reliability problems.First, we introduced the method by constructing the joint probability density function of correlation variable based on Copula function. This method includes the solution of correlation parameterθ and the selection of the structural type of the Copula function. Second, we introduced Nataf-Monte Carlo method. Third, we proposed direct integration based on dual neural network, which is in turn based on the integral form of reliability computation and the normalization method of the integral area. Fourth, we compared the proposed technique with Monte Carlo method and verified its effectiveness in simulation. Finally, the full-text conclusion and future prospects were presented.


## 2. Construction of Joint Probability Density Function Based on Copula Function for Correlated Variables

### 2.1. Copula Function

Copula theory was first proposed by Sklar in 1959. Sklar pointed out that any multidimensional joint distribution function could be decomposed into a corresponding edge distribution function and a Copula function. The Copula function determines the correlation among variables, including the size of the correlation coefficient and the type of the correlation structure [10]. According to its strict definition stated by Nelsen, the Copula function is a function that associates the joint distribution function of the variable with its edge distribution function. In essence, Copula function is also a joint distribution function. For n-dimensional cases, the Copula function was defined as the n-dimensional joint distribution function. The edge distribution was [0,1] in the [0,1]n space [9].According to Sklar theorem [10], the joint distribution function Fx1,x2,…,xn of variables X1, X2,…,Xn could be expressed as(1)Fx1,x2,…,xn=CF1x1,F2x2,…,Fnxn;θ=Cu1,u2,…,un;θ,where ui=Fixi is the edge distribution function of the variable Xi: i=1,2,…,n. Cu1,u2,…,un;θ is the Copula function, and θ is the correlation parameter of the Copula function.From (1), the joint probability density function fx1,x2,…,xn could be obtained as(2)fx1,x2,…,xn=f1x1f2x2⋯fnxncF1x1,F2x2,…,Fnxn;θ=cu1,u2,…,un;θ∏i=1nfixi,where fixi is the edge probability density function of variables X1,X2,…,Xn and cu1,u2,…,un;θ=∂nCu1,u2,…,un;θ/∂u1∂u2⋯∂un is the density function for Cu1,u2,…,un;θ.If the edge distribution function of variablesX1,X2,…,Xn and the Copula function were known, then the multidimensional distribution model of the variables could be established using (1) and (2).


### 2.2. Correlation Parameterθ of Copula Function

The correlation parameterθ of the Copula function characterizes the correlation among variables. The method of solving the correlation parameter θ is different for different types of Copula functions.When the Copula function of multidimensional variableX1,X2,…,Xn was an Ellipse Copula function, the number of correlation parameters θij was the same as that of the correlation coefficients among variables and had one-to-one correspondence. The correlation parameter of any two variables Xi and Xj was θij. The relationship between θij and Pearson linear correlation coefficient ρij was [9](3)ρij=∫-∞∞∫-∞∞xi-μiσixj-μjσjfixifjxjcui,uj;θijdxidxj.When the Copula function of the multidimensional variableX1,X2,…,Xn was the Archimedean Copula function, the multidimensional Archimedean Copula function had only one correlation parameter θ because only one generator existed. This parameter describes the overall correlation between the multidimensional variables X1,X2,…,Xn. Maximum likelihood estimation is generally used to solve the correlation parameter θ [44].(4)Lθ=∑i=1nln⁡cu1,u2,…,un;θ.


### 2.3. Selection of the Optimal Copula Function

Different Copula functions could describe different correlations among variables. Table1 lists various Copula function types [9]. Therefore, selecting the optimal Copula function, which could best fit the correlation among variables, is necessary when constructing the joint probability density function among variables. In this paper, we used Akaike information criterion (AIC) and Bayesian information criterion (BIC) to select the optimal Copula function. The best Copula function has the smallest AIC or BIC values and is well fitted for the correlated structure among variables. For multidimensional distribution models, the AIC and BIC values were expressed as follows [45, 46]: (5)AIC=2k-2∑i=1nln⁡cu1,u2,…,un;θBIC=kln⁡n-2∑i=1nln⁡cu1,u2,…,un;θ,where k is the number of correlation parameters θ. In the Ellipse Copula function, k=0.5n(n-1), whereas in the Archimedean Copula function, k=1.Table 1
Copula function types.
Copula types
Copula distribution functionCu1,u2,…,un;θ
θ range of values
Gaussian
ϕ
n
ϕ
-
1
u
1
,
ϕ
-
1
u
2
,
…
,
ϕ
-
1
u
n
;
θ
-
1,1
Gumbel
e
x
p
-
∑
i
=
1
n
-
ln
⁡
u
i
θ
1
/
θ
1
,
∞
Clayton
∑
i
=
1
n
u
i
-
θ
-
n
+
1
-
1
/
θ
0
,
∞
Frank
-
1
θ
l
n
1
+
∏
i
=
1
n
e
-
θ
u
i
-
1
e
-
θ
-
1
n
-
1
n
=
2
,
θ
∈
-
∞
,
∞
∩
θ
≠
0 n≥3,θ∈0,∞
Number 16
1
2
A
+
A
2
+
4
θ A=∑i=1nui-n+1-θ∑i=1n1ui-n+1
n
=
2
,
θ
∈
0
,
∞ n≥3,θ∈0.1863,∞


## 2.1. Copula Function

Copula theory was first proposed by Sklar in 1959. Sklar pointed out that any multidimensional joint distribution function could be decomposed into a corresponding edge distribution function and a Copula function. The Copula function determines the correlation among variables, including the size of the correlation coefficient and the type of the correlation structure [10]. According to its strict definition stated by Nelsen, the Copula function is a function that associates the joint distribution function of the variable with its edge distribution function. In essence, Copula function is also a joint distribution function. For n-dimensional cases, the Copula function was defined as the n-dimensional joint distribution function. The edge distribution was [0,1] in the [0,1]n space [9].According to Sklar theorem [10], the joint distribution function Fx1,x2,…,xn of variables X1, X2,…,Xn could be expressed as(1)Fx1,x2,…,xn=CF1x1,F2x2,…,Fnxn;θ=Cu1,u2,…,un;θ,where ui=Fixi is the edge distribution function of the variable Xi: i=1,2,…,n. Cu1,u2,…,un;θ is the Copula function, and θ is the correlation parameter of the Copula function.From (1), the joint probability density function fx1,x2,…,xn could be obtained as(2)fx1,x2,…,xn=f1x1f2x2⋯fnxncF1x1,F2x2,…,Fnxn;θ=cu1,u2,…,un;θ∏i=1nfixi,where fixi is the edge probability density function of variables X1,X2,…,Xn and cu1,u2,…,un;θ=∂nCu1,u2,…,un;θ/∂u1∂u2⋯∂un is the density function for Cu1,u2,…,un;θ.If the edge distribution function of variablesX1,X2,…,Xn and the Copula function were known, then the multidimensional distribution model of the variables could be established using (1) and (2).


## 2.2. Correlation Parameterθ of Copula Function

The correlation parameterθ of the Copula function characterizes the correlation among variables. The method of solving the correlation parameter θ is different for different types of Copula functions.When the Copula function of multidimensional variableX1,X2,…,Xn was an Ellipse Copula function, the number of correlation parameters θij was the same as that of the correlation coefficients among variables and had one-to-one correspondence. The correlation parameter of any two variables Xi and Xj was θij. The relationship between θij and Pearson linear correlation coefficient ρij was [9](3)ρij=∫-∞∞∫-∞∞xi-μiσixj-μjσjfixifjxjcui,uj;θijdxidxj.When the Copula function of the multidimensional variableX1,X2,…,Xn was the Archimedean Copula function, the multidimensional Archimedean Copula function had only one correlation parameter θ because only one generator existed. This parameter describes the overall correlation between the multidimensional variables X1,X2,…,Xn. Maximum likelihood estimation is generally used to solve the correlation parameter θ [44].(4)Lθ=∑i=1nln⁡cu1,u2,…,un;θ.


## 2.3. Selection of the Optimal Copula Function

Different Copula functions could describe different correlations among variables. Table1 lists various Copula function types [9]. Therefore, selecting the optimal Copula function, which could best fit the correlation among variables, is necessary when constructing the joint probability density function among variables. In this paper, we used Akaike information criterion (AIC) and Bayesian information criterion (BIC) to select the optimal Copula function. The best Copula function has the smallest AIC or BIC values and is well fitted for the correlated structure among variables. For multidimensional distribution models, the AIC and BIC values were expressed as follows [45, 46]: (5)AIC=2k-2∑i=1nln⁡cu1,u2,…,un;θBIC=kln⁡n-2∑i=1nln⁡cu1,u2,…,un;θ,where k is the number of correlation parameters θ. In the Ellipse Copula function, k=0.5n(n-1), whereas in the Archimedean Copula function, k=1.Table 1
Copula function types.
Copula types
Copula distribution functionCu1,u2,…,un;θ
θ range of values
Gaussian
ϕ
n
ϕ
-
1
u
1
,
ϕ
-
1
u
2
,
…
,
ϕ
-
1
u
n
;
θ
-
1,1
Gumbel
e
x
p
-
∑
i
=
1
n
-
ln
⁡
u
i
θ
1
/
θ
1
,
∞
Clayton
∑
i
=
1
n
u
i
-
θ
-
n
+
1
-
1
/
θ
0
,
∞
Frank
-
1
θ
l
n
1
+
∏
i
=
1
n
e
-
θ
u
i
-
1
e
-
θ
-
1
n
-
1
n
=
2
,
θ
∈
-
∞
,
∞
∩
θ
≠
0 n≥3,θ∈0,∞
Number 16
1
2
A
+
A
2
+
4
θ A=∑i=1nui-n+1-θ∑i=1n1ui-n+1
n
=
2
,
θ
∈
0
,
∞ n≥3,θ∈0.1863,∞


## 3. Monte Carlo Method Based on Nataf Inverse Transformation

Monte Carlo is a numerical simulation method that is widely used in structural reliability engineering [47, 48]. Based on statistical sampling theory, Monte Carlo method solves structural reliability by statistical sampling or random simulation of random variables. This method could obtain the true failure probability if the distribution type of the variable is known, and the number of sampling times is sufficient. In reliability calculation, Monte Carlo method is considered a precise calculation method [49]. In this section, the Nataf-Monte Carlo method was introduced to calculate the structural reliability of correlated variables. This method was defined as Nataf-MCS.The Nataf distribution model was proposed by Nataf in 1962 [50]. The famous Nataf transform is a representation form of the Nataf distribution model, which includes positive and inverse transforms. Nataf transforms could convert correlated variables into independent variables. Nataf inverse transforms convert the independent variables into correlated variables.If the cumulative distribution function of the variablesX=X1,X2,…,Xn were F1x1,F2x2,…,Fnxn, then the correlation coefficient matrix was ρ. The concrete steps were given by the Nataf-Monte Carlo method to calculate the structural reliability of the relevant variables as follows [49–51].(1) The correlation coefficient matrixρ0 of the standard normal distribution variables Z=Z1,Z2,…,Zn was solved.(6)ρij=∫-∞∞∫-∞∞xi-μiσixj-μjσjfixifjxjφziφzjφzi,zj;ρ0ijdxidxj.Ang et al. established the following empirical equation [49] for various distribution types under normal circumstances ρij≤ρ0ij:(7)ρ0ij=ξρij,where the coefficient ξ is a function of the correlation coefficient ρij and the edge distribution Fixi. The equation for the coefficient ξ was provided in literature [51].(2) The standard normal distribution variablesZ=Z1,Z2,…,Zn were correlated, and their correlation coefficient matrices were ρ0. ρ0 was decomposed by Gorensky to obtain the lower triangular matrix L.(8)ρ0=LLT.(3) The traditional random number generation method was applied to produce independent variablesY=Y1,Y2,…,Yn, which all followed the standard normal distribution.(4) Based on the lower triangular matrixL, the independent normal random variables Y=Y1,Y2,…,Yn were transformed into standard normal distribution variables Z=Z1,Z2,…,Zn.(9)Z=LY.(5) The relevant standard normal distribution variablesZ=Z1,Z2,…,Zn were mapped to the original spatial variables X=X1,X2,…,Xn using the principle of equal probability transformation.(10)X=F-1ϕZ.(6) The original spatial variablesX=X1,X2,…,Xn were substituted into the functional functions of the reliability analysis, and the structure response was calculated: Gi=GXi. The number of times Gi<0 was the number of failure sample points M. The unbiased estimate of the structural failure probability was [49] as follows:(11)Pf=MN.


## 4. Structural Reliability Calculations Based on Dual Neural Network and Direct Integration for Correlated Variables

Structural reliability could be expressed as the probability of functionG>0. G could be expressed as a function of basic variables X=[x1,x2,…,xn]; that is, G=gx1,x2,…,xn. The joint probability density function of the correlated variables was fx1,x2,…,xn, introducing a weight function: (12)hx=1,x>00,x≤0,(13)yx1,x2,…,xn=hgx1,x2,…,xnfx1,x2,…,xn=hgx1,x2,…,xnfx1fx2⋯fxn·cx1,x2,…,xn;θ.Structural reliability could be expressed as multiple integrals as follows:(14)P=∫μn-bnσnμn+bnσn⋯∫μ1-a1σ1μ1+a1σ1yx1,x2,…,xndx1dx2⋯dxn=∫μn-bnσnμn+bnσn⋯∫μ1-a1σ1μ1+a1σ1hgx1,x2,…,xnfx1fx2⋯fxn·cx1,x2,…,xndx1dx2⋯dxn.In (16), μi and σi were the mean and standard deviation, respectively, of random variable i. ai and bi were positive integers, which were determined by accuracy requirements. In general, ai≥3 and bi≥3.The construction of dual neural network for multiple integral calculation was discussed in the following section.We constructed a single hidden layer BP neural network to establish the mapping relationship between the input variableX and the original function Y. The network structure is shown in Figure 1.Figure 1
Original function network structure.The relationships between the output variables and input variables could be expressed as(15)Y=∑j=1mfs∑i=1nwji1xi+ϑjwj2+b.The derivation of (17) could be given by(16)y=∂nY∂x1∂x2⋯∂xn=∑j=1mfsn∑i=1nwji1xi+ϑjwj11wj21⋯wjn1wj2,where Wj=wj11wj21⋯wjn1wj2. Equation (16) could be rewritten as a function of the relationship between the output and input variables as follows:(17)y=∑j=1mfsn∑i=1nwji1xi+ϑjWj.The neural network completed the mapping relationship between the input variableX and the integrand y. The network structure is shown in Figure 2.Figure 2
Integrand network structure.Let the original function of the integrandy(x1,x2,…,xn) be Y(x1,x2,…,xn). According to multiple integral theory, function Y is the weighted algebraic sum of each vertex in a hypercube D=x1,x2,…,xn∣x11≤x1≤x12,x21≤x2≤x22,…,xn1≤xn≤xn2 and could be expressed as(18)J=∑m1=12∑m2=12⋯∑mn=12-1∑k=1nmkYx1m1,x2m2,…,xnmn.Thus, the original function and integrand neural networks were combined as a dual neural network. The two networks, namely, integrand network and original function network, had three layers,n inputs, single output, and m unit in hidden layer. In the integrand network, the connection weights and threshold of input layer to hidden layer unit were wji1 and ϑj, respectively; the connection weights and threshold of hidden layer to output layer were Wj and 0, respectively; and the activation function of the hidden layer units was fsn. In the original function network, the connection weights and threshold of input layer to hidden layer unit were wji1 and ϑj, respectively; the connection weights and threshold of hidden layer to output layer were Wj/∏i=1nwji1 and b, respectively; and the activation function of the hidden layer units was fs. When the integrand network approximated the integrand in the integral, the original function network simultaneously approximated the original function.


## 5. Numerical Examples

### 5.1. Example 1

The performance function of a system was assumed asgX=x12-3x2, where x1 is the lognormal distribution random variables with mean value μ1=12.42, standard deviation σ1=2.484, and coefficient of variation η1=0.20. x2 was an extreme value distribution of type I with mean value μ2=17.75, standard deviation σ2=5, and coefficient of variation η2=0.2817. The correlation coefficient of x1 and x2 was ρ=-0.85. The structural reliability was calculated.Based on the Copula function, the 2D joint probability density function of the random variablesx1 and x2 was established and could be written as(19)fx1,x2=f1x1f2x2cF1x1,F2x2;θ,where(20)f1x1=12πx1σYexp⁡-12ln⁡x1-μYσY2,f2x2=αexp⁡-x2-δα-exp⁡-x2-δα,where μY=2.5, σY=0.198, α=0.246, and δ=15.4.Calculated by AIC and BIC criteria, in example 1, Plackett Copula, Frank Copula, number 16 Copula, and Gaussian Copula were selected to construct the 2D joint probability density function of random variablesx1 and x2. The values of the correlation parameter θ of each Copula function were obtained from (3) and are shown in Table 2.Table 2
Correlation parameterθ in example 1.
Copula function
Plackett Copula
Frank Copula
Number 16 Copula
Gaussian Copula
Correlation parameterθ
0.0229
-
11.3478
0.0038
-
0.89f
x
1
,
x
2 and gX were substituted for (13) to obtain the integrand.The integrand neural network had 2 inputs, single output, and 20 units in hidden layer. The transfer functions of hidden layer and output layer werefs=ex and Purelin function.The ranges ofx1 and x2μ1-4σ1,μ1+5σ1 and μ2-10σ2,μ2+6σ2 were divided into 50 parts. Two variables crossed each other to form the network input samples. The network output value y of the corresponding sample points was calculated by the integrand. The training sample set is shown in Table 3.Table 3
Training sample set of integrand network in example 1.
x
2
Y
x
1
2.4920
2.9479
3.4038
⋯
24.3741
24.8300
-
32.2500
0
0
0
⋯
0
0
-
30.6173
0
0
0
⋯
0
0
⋯
⋯
⋯
⋯
⋯
⋯
⋯
42.8520
1.188
0
e
-
30
3.306
8
e
-
19
6.750
0
e
-
16
⋯
0
0
44.4847
8.072
8
e
-
31
2.783
9
e
-
19
5.861
2
e
-
16
⋯
0
0
46.1173
5.483
2
e
-
31
2.367
5
e
-
19
5.014
2
e
-
16
⋯
0
0
47.7500
3.726
7
e
-
31
2.032
9
e
-
19
4.162
2
e
-
16
⋯
0
0This example used Levenberg-Marquardt training algorithm. Error convergence curve of training integrand network after 1000 steps is shown in Figure3.Figure 3
Training error curve of integrand network in example 1.The relationship of the dual neural network was used to construct the original function network. The original function network structure was the same as the integrand network. The integrand network training sample set of each vertex in a hypercube was calculated and simulated using the original function network. The sample set is shown in Table4. The results were substituted into (18) to obtain the result of reliability.Table 4
Simulative sample set of original function network in example 1.
x
1
24.8300
24.8300
2.4920
2.4920
x
2
47.7500
-
32.2500
47.7500
-
32.2500For comparison, we provide the results from Nataf-Monte Carlo method and the proposed method. The results calculated using Nataf-Monte Carlo method with sampling 1,000,000 were viewed as accurate solutions The relative errors are shown in Table5.Table 5
Example 1 reliability calculation results.
Method
Nataf-MCS
Direct integral method based on dual neural network
Plackett Copula
Frank Copula
Number 16 Copula
Gaussian Copula
Reliability
0.9609
0.9654
0.9528
0.9648
0.9799
Relative error (%)
0
0.47
0.84
0.41
2.00The ratio of the calculation workloads of the Nataf-MCS and the proposed method can be obtained from the following equation:(21)γ=W1NmW1Nn1+W2Nn2,where W1 is a single workload for structural analysis. Nm is the calculation times of the Nataf-MCS. Nn1 is the calculation times required to construct the network training sample. W2 is a single workload for neural network training. Nn2 is the neural network training times.We can compare the computational efficiency of the two methods using the ratioγ of calculation workload. In this paper, W1=10W2. From the calculation in this example, γ=80.From Table5, the accuracy of the proposed method was close to that of the Nataf-MCS. The maximum relative error was only 2.00%. Based on ratio γ=80, the calculation workload of Nataf-MCS was 80 times higher than that of the proposed method. With regard to computational efficiency, the proposed method saves calculation time and greatly improves the computational efficiency of reliability.


### 5.2. Example 2

In the cantilever beam shown in Figure4,Yis the horizontal load on the beam, Z is the lateral load on the beam, and w represents the beam width. The performance function of structure gY,Z,w=45-43Y/w-160Z/w2. Considering the correlation between variables, the reliability of the cantilever beam was analyzed [52].Figure 4
Cantilever structure.Based on the experimental data in [52], the mean, standard deviation, and edge distribution functions of variables Y, Z, and w were obtained and are listed in Table 6.Table 6
Statistical parameters and distribution types of variables in example 2.
Project
Y (lb)
Z (lb)
w (in)
Meanμ
988.72
506.28
1.98
Standard deviationσ
120.78
132.48
0.27
Edge distribution function
Normal
Normal
NormalCalculated by AIC and BIC criteria, in example 2, Frank Copula, Gumbel Copula, and Clayton Copula were selected to construct the 3D joint probability density function of random variablesY,Z,and w.(22)fY,Z,w=fYfZfw·cFY,FZ,Fw;θ.Based on the maximum likelihood estimation in (4), the values of each Copula function’s correlation parameter θ were calculated and are listed in Table 7.Table 7
Values of correlation parameterθ in example 2.
Copula function
Frank Copula
Gumbel Copula
Clayton Copula
Correlation parameterθ
32.47
7.66
9.05The 3D joint probability density functionfY,Z,w and the performance function gY,Z,w in example 2 were substituted into (13) to calculate the integrand.The structure of the integrand neural network could also be constructed. The integrand neural network had three inputs and single output. After adjustment, the number of hidden layer units ism=10. The transfer functions of hidden layer and output layer were fs=ex and Purelin function, respectively.VariablesY, Z, and w were separately divided into 30 parts in the range of μi-3σi,μi+3σi. The three variables intersect each other to form the integrand network input sample point. The network output value y of the corresponding sample point was calculated by the integrand function, and the training sample point set is shown in Table 8.Table 8
Training sample set of integrand network in example 2.
Y
Z
Y
w
1.17
1.23
⋯
2.74
2.80
626.40
108.84
3.740
6
e
-
019
7.247
9
e
-
019
⋯
2.056
6
e
-
017
9.277
7
e
-
019
626.40
136.25
6.176
9
e
-
019
1.192
6
e
-
018
⋯
3.327
8
e
-
017
1.472
1
e
-
018
⋯
⋯
⋯
⋯
⋯
⋯
⋯
626.40
903.70
1.158
1
e
-
020
2.159
2
e
-
020
⋯
4.997
3
e
-
019
1.670
4
e
-
020
⋯
⋯
⋯
⋯
⋯
⋯
⋯
1351.10
108.84
7.819
5
e
-
021
1.456
8
e
-
020
⋯
3.356
0
e
-
019
1.111
9
e
-
020
1351.10
136.25
1.238
5
e
-
020
2.306
3
e
-
020
⋯
5.299
6
e
-
019
1.747
2
e
-
020
⋯
⋯
⋯
⋯
⋯
⋯
⋯
1351.10
903.70
1.365
4
e
-
022
2.524
1
e
-
022
⋯
5.542
1
e
-
021
1.662
5
e
-
022We used the Levenberg-Marquardt training algorithm. When the training reached 100 steps, the error of the integrand network reached the convergence state. The error convergence curve is shown in Figure5.Figure 5
Training error curve of integrand network in example 2.The integrand network training sample set of each vertex in a hypercube was calculated and simulated with the original function network. The sample set is shown in Table9. The results were substituted into (18) to obtain the result of reliability.Table 9
Simulative sample set of original function network in example 2.
Y
626.40
1351.10
626.40
1351.10
626.40
1351.10
626.40
1351.10
Z
108.84
903.70
108.84
108.84
903.70
108.84
903.70
903.70
w
1.17
2.80
2.80
2.80
1.17
1.17
2.80
1.17Nataf-MCS was introduced in Section3. The sampling was 1,000,000 times, and the reliability of example 2 was 0.9998. In this example, γ=35.71. Table 10 shows the reliability values calculated by the two methods.Table 10
Example 2 reliability calculation results.
Method
Nataf-MCS
Direct integral method based on dual neural network
Frank Copula
Gumbel Copula
Clayton Copula
Reliability
0.9998
0.9960
0.9913
0.9936
Relative error (%)
0
0.38
0.85
0.62In this example, three kinds of Archimedean Copula functions were selected to construct the joint probability density function of variablesY, Z, and w. The reliability of the proposed method was 0.85%, which was higher than that of Nataf-MCS. The ratio γ of calculation workload was 35.71, indicating that the proposed method has superior computational efficiency and can reduce the calculation time.


### 5.3. Example 3

The soil parameters of rock and soil include liquidity index (LI), undrained shear strengthsu, remolded undrained shear strength sure, preconsolidation stress σp′, and vertical effective stress σv′. The experimental data in literature [53] were used as an example to illustrate the calculation steps of the proposed method. The performance function of the soil parameters gLI,su,sure,σp′,σv′=σp′sure+LIσp′-σv′-su2 was considered, and the reliability of rock and soil was analyzed under normal use limit state. Test data for soil parameters are shown in Table 11.Table 11
Test data for soil parameters.
Number
LI
s
u (kPa)
s
u
re (kPa)
σ
p
′ (kPa)
σ
v
′ (kPa)
1
1.28
2.58
0.22
7.94
9.56
2
1.27
4.74
0.16
29.41
12.82
3
1.45
5.42
0.39
17.14
16.44
4
1.20
5.82
0.59
27.63
20.06
5
1.26
6.97
0.65
23.37
24.04
⋯
⋯
⋯
⋯
⋯
⋯
343
2.11
12.26
1.53
80.24
61.90
344
2.05
14.10
1.53
86.36
68.02
345
1.57
15.20
1.84
94.76
75.66The meanμ, standard deviation σ, and edge distribution functions type of each variable were obtained by analyzing the experimental data as shown in Table 12.Table 12
Statistical parameters and distribution types of variables in example 3.
Items
LI
s
u (kPa)
s
u
re (kPa)
σ
p
′ (kPa)
σ
v
′ (kPa)
Meanμ
1.25
31.01
2.51
105.82
66.63
Standard deviationσ
0.61
29.48
3.81
103.21
53.48
Edge distributionfunction
Lognormal
Lognormal
Lognormal
Lognormal
LognormalFrom the calculation, we obtained the correlation coefficientρij of any two variables to form the matrix ρ.(23)ρ=LIsuresureσp′σv′LIsusureσp′σv′10.05-0.50-0.06-0.210.0510.170.840.57-0.500.1710.300.42-0.060.840.3010.73-0.210.570.420.731.Given the positive and negative correlations among variablesLI, su, sure, σp′, and σv′, we need to select the appropriate Copula function from Elliptical Copula function to construct the joint probability density function fLI,su,sure,σp′,σv′. Using AIC and BIC criteria, Gaussian Copula was selected for example 2. Afterwards, (3) was used to obtain any two variables related parameters θij and to form matrix θ.(24)θ=LIsuresureσp′σv′LIsusureσp′σv′10.07-0.91-0.08-0.260.0710.250.880.64-0.910.2510.410.53-0.080.880.4110.78-0.260.670.530.781.The 5D joint probability density function of variablesLI, su, sure, σp′, and σv′ can be expressed as the following equation:(25)fLI,su,sure,σp′,σv′=fLIfsufsurefσp′fσv′·cFLI,Fsu,Fsure,Fσp′,Fσv′;θ.The 5D joint probability density functionfLI,su,sure,σp′,σv′ and the performance function gLI,su,sure,σp′,σv′ were substituted into (13) to calculate the integrand.The structure of the integrand neural network could also be constructed. The integrand neural network had five inputs and single output. After adjustment, the number of hidden layer units wasm=5. The transfer functions of hidden layer and output layer were fs=ex and Purelin function, respectively.VariablesLI, su, sure, σp′, and σv′ are divided into 10 parts in the range of μi-4σi,μi+5σi. The five variables intersect each other to form the input sample points of the integrand network. The network output value y of the corresponding sample points was calculated by the integrand function. Given that each variable follows the lognormal distribution, we get μi-4σi=0.1 when the left interval of the variable μi-4σi≤0. The training sample set is shown in Table 13.Table 13
Training sample set of integrand network in example 3.
LI
s
u
s
u
re
y
σ
p
′
0.10
0.10
⋯
415.45
415.45
σ
v
′
0.10
25.32
201.85
227.07
0.10
0.10
0.10
1.072
8
e
-
006
5.355
0
e
-
007
⋯
4.194
5
e
-
011
2.767
2
e
-
011
0.10
0.10
1.64
3.355
9
e
-
007
1.681
9
e
-
007
⋯
1.389
7
e
-
011
9.205
3
e
-
012
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
0.10
0.10
13.94
3.854
6
e
-
009
1.994
8
e
-
009
⋯
2.527
4
e
-
013
1.728
5
e
-
013
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
3.08
119.45
0.10
1.330
3
e
-
010
6.429
6
e
-
011
⋯
5.826
7
e
-
014
3.721
8
e
-
014
3.08
119.45
1.64
3.859
7
e
-
011
1.872
9
e
-
011
⋯
1.790
5
e
-
014
1.148
3
e
-
014
⋯
⋯
⋯
⋯
⋯
⋯
⋯
3.08
119.45
13.94
2.426
6
e
-
013
1.215
9
e
-
013
⋯
1.782
3
e
-
016
1.180
2
e
-
016We used the Levenberg-Marquardt training algorithm. When the training reached 100 steps, the error of the integrand network reached the convergence state. The error convergence curve is shown in Figure6.Figure 6
Training error curve of integrand network in example 3.We used the same method for the previous example. The sample points of the original function network were obtained and are shown in Table14. The simulation values were substituted into (18) to obtain the reliability of the structure as shown in Table 15.Table 14
Simulative sample set of original function network in example 3.
LI
0.10
3.08
0.10
3.08
⋯
0.10
3.08
s
u
0.10
119.45
119.45
0.10
⋯
0.10
119.45
s
u
re
0.10
13.94
13.94
0.10
⋯
13.94
0.10
σ
p
′
0.10
415.45
415.45
0.10
⋯
0.10
415.45
σ
v
′
0.10
227.07
227.07
0.10
⋯
227.07
0.10Table 15
Example 3 reliability calculation results.
Method
Nataf-MCS
Direct integral method based on dual neural network
Gaussian Copula
Reliability
0.9594
0.9483
Relative error (%)
0
1.16The calculation workload of the two methods was substituted into (21), thus obtaining γ=9.90.The reliability of 5D correlated variables is analyzed using example 3. The reliability of the proposed method was extremely close to that of Nataf-MCS, and the relative error was only 1.16%. Fromγ=9.90, for the same problem, the workload of the two methods varied greatly. The neural network method has advantages in computational efficiency and reduces the calculation time for structural reliability analysis. This example shows that the proposed method can efficiently solve the reliability problem of multidimensional correlated variables.


## 5.1. Example 1

The performance function of a system was assumed asgX=x12-3x2, where x1 is the lognormal distribution random variables with mean value μ1=12.42, standard deviation σ1=2.484, and coefficient of variation η1=0.20. x2 was an extreme value distribution of type I with mean value μ2=17.75, standard deviation σ2=5, and coefficient of variation η2=0.2817. The correlation coefficient of x1 and x2 was ρ=-0.85. The structural reliability was calculated.Based on the Copula function, the 2D joint probability density function of the random variablesx1 and x2 was established and could be written as(19)fx1,x2=f1x1f2x2cF1x1,F2x2;θ,where(20)f1x1=12πx1σYexp⁡-12ln⁡x1-μYσY2,f2x2=αexp⁡-x2-δα-exp⁡-x2-δα,where μY=2.5, σY=0.198, α=0.246, and δ=15.4.Calculated by AIC and BIC criteria, in example 1, Plackett Copula, Frank Copula, number 16 Copula, and Gaussian Copula were selected to construct the 2D joint probability density function of random variablesx1 and x2. The values of the correlation parameter θ of each Copula function were obtained from (3) and are shown in Table 2.Table 2
Correlation parameterθ in example 1.
Copula function
Plackett Copula
Frank Copula
Number 16 Copula
Gaussian Copula
Correlation parameterθ
0.0229
-
11.3478
0.0038
-
0.89f
x
1
,
x
2 and gX were substituted for (13) to obtain the integrand.The integrand neural network had 2 inputs, single output, and 20 units in hidden layer. The transfer functions of hidden layer and output layer werefs=ex and Purelin function.The ranges ofx1 and x2μ1-4σ1,μ1+5σ1 and μ2-10σ2,μ2+6σ2 were divided into 50 parts. Two variables crossed each other to form the network input samples. The network output value y of the corresponding sample points was calculated by the integrand. The training sample set is shown in Table 3.Table 3
Training sample set of integrand network in example 1.
x
2
Y
x
1
2.4920
2.9479
3.4038
⋯
24.3741
24.8300
-
32.2500
0
0
0
⋯
0
0
-
30.6173
0
0
0
⋯
0
0
⋯
⋯
⋯
⋯
⋯
⋯
⋯
42.8520
1.188
0
e
-
30
3.306
8
e
-
19
6.750
0
e
-
16
⋯
0
0
44.4847
8.072
8
e
-
31
2.783
9
e
-
19
5.861
2
e
-
16
⋯
0
0
46.1173
5.483
2
e
-
31
2.367
5
e
-
19
5.014
2
e
-
16
⋯
0
0
47.7500
3.726
7
e
-
31
2.032
9
e
-
19
4.162
2
e
-
16
⋯
0
0This example used Levenberg-Marquardt training algorithm. Error convergence curve of training integrand network after 1000 steps is shown in Figure3.Figure 3
Training error curve of integrand network in example 1.The relationship of the dual neural network was used to construct the original function network. The original function network structure was the same as the integrand network. The integrand network training sample set of each vertex in a hypercube was calculated and simulated using the original function network. The sample set is shown in Table4. The results were substituted into (18) to obtain the result of reliability.Table 4
Simulative sample set of original function network in example 1.
x
1
24.8300
24.8300
2.4920
2.4920
x
2
47.7500
-
32.2500
47.7500
-
32.2500For comparison, we provide the results from Nataf-Monte Carlo method and the proposed method. The results calculated using Nataf-Monte Carlo method with sampling 1,000,000 were viewed as accurate solutions The relative errors are shown in Table5.Table 5
Example 1 reliability calculation results.
Method
Nataf-MCS
Direct integral method based on dual neural network
Plackett Copula
Frank Copula
Number 16 Copula
Gaussian Copula
Reliability
0.9609
0.9654
0.9528
0.9648
0.9799
Relative error (%)
0
0.47
0.84
0.41
2.00The ratio of the calculation workloads of the Nataf-MCS and the proposed method can be obtained from the following equation:(21)γ=W1NmW1Nn1+W2Nn2,where W1 is a single workload for structural analysis. Nm is the calculation times of the Nataf-MCS. Nn1 is the calculation times required to construct the network training sample. W2 is a single workload for neural network training. Nn2 is the neural network training times.We can compare the computational efficiency of the two methods using the ratioγ of calculation workload. In this paper, W1=10W2. From the calculation in this example, γ=80.From Table5, the accuracy of the proposed method was close to that of the Nataf-MCS. The maximum relative error was only 2.00%. Based on ratio γ=80, the calculation workload of Nataf-MCS was 80 times higher than that of the proposed method. With regard to computational efficiency, the proposed method saves calculation time and greatly improves the computational efficiency of reliability.


## 5.2. Example 2

In the cantilever beam shown in Figure4,Yis the horizontal load on the beam, Z is the lateral load on the beam, and w represents the beam width. The performance function of structure gY,Z,w=45-43Y/w-160Z/w2. Considering the correlation between variables, the reliability of the cantilever beam was analyzed [52].Figure 4
Cantilever structure.Based on the experimental data in [52], the mean, standard deviation, and edge distribution functions of variables Y, Z, and w were obtained and are listed in Table 6.Table 6
Statistical parameters and distribution types of variables in example 2.
Project
Y (lb)
Z (lb)
w (in)
Meanμ
988.72
506.28
1.98
Standard deviationσ
120.78
132.48
0.27
Edge distribution function
Normal
Normal
NormalCalculated by AIC and BIC criteria, in example 2, Frank Copula, Gumbel Copula, and Clayton Copula were selected to construct the 3D joint probability density function of random variablesY,Z,and w.(22)fY,Z,w=fYfZfw·cFY,FZ,Fw;θ.Based on the maximum likelihood estimation in (4), the values of each Copula function’s correlation parameter θ were calculated and are listed in Table 7.Table 7
Values of correlation parameterθ in example 2.
Copula function
Frank Copula
Gumbel Copula
Clayton Copula
Correlation parameterθ
32.47
7.66
9.05The 3D joint probability density functionfY,Z,w and the performance function gY,Z,w in example 2 were substituted into (13) to calculate the integrand.The structure of the integrand neural network could also be constructed. The integrand neural network had three inputs and single output. After adjustment, the number of hidden layer units ism=10. The transfer functions of hidden layer and output layer were fs=ex and Purelin function, respectively.VariablesY, Z, and w were separately divided into 30 parts in the range of μi-3σi,μi+3σi. The three variables intersect each other to form the integrand network input sample point. The network output value y of the corresponding sample point was calculated by the integrand function, and the training sample point set is shown in Table 8.Table 8
Training sample set of integrand network in example 2.
Y
Z
Y
w
1.17
1.23
⋯
2.74
2.80
626.40
108.84
3.740
6
e
-
019
7.247
9
e
-
019
⋯
2.056
6
e
-
017
9.277
7
e
-
019
626.40
136.25
6.176
9
e
-
019
1.192
6
e
-
018
⋯
3.327
8
e
-
017
1.472
1
e
-
018
⋯
⋯
⋯
⋯
⋯
⋯
⋯
626.40
903.70
1.158
1
e
-
020
2.159
2
e
-
020
⋯
4.997
3
e
-
019
1.670
4
e
-
020
⋯
⋯
⋯
⋯
⋯
⋯
⋯
1351.10
108.84
7.819
5
e
-
021
1.456
8
e
-
020
⋯
3.356
0
e
-
019
1.111
9
e
-
020
1351.10
136.25
1.238
5
e
-
020
2.306
3
e
-
020
⋯
5.299
6
e
-
019
1.747
2
e
-
020
⋯
⋯
⋯
⋯
⋯
⋯
⋯
1351.10
903.70
1.365
4
e
-
022
2.524
1
e
-
022
⋯
5.542
1
e
-
021
1.662
5
e
-
022We used the Levenberg-Marquardt training algorithm. When the training reached 100 steps, the error of the integrand network reached the convergence state. The error convergence curve is shown in Figure5.Figure 5
Training error curve of integrand network in example 2.The integrand network training sample set of each vertex in a hypercube was calculated and simulated with the original function network. The sample set is shown in Table9. The results were substituted into (18) to obtain the result of reliability.Table 9
Simulative sample set of original function network in example 2.
Y
626.40
1351.10
626.40
1351.10
626.40
1351.10
626.40
1351.10
Z
108.84
903.70
108.84
108.84
903.70
108.84
903.70
903.70
w
1.17
2.80
2.80
2.80
1.17
1.17
2.80
1.17Nataf-MCS was introduced in Section3. The sampling was 1,000,000 times, and the reliability of example 2 was 0.9998. In this example, γ=35.71. Table 10 shows the reliability values calculated by the two methods.Table 10
Example 2 reliability calculation results.
Method
Nataf-MCS
Direct integral method based on dual neural network
Frank Copula
Gumbel Copula
Clayton Copula
Reliability
0.9998
0.9960
0.9913
0.9936
Relative error (%)
0
0.38
0.85
0.62In this example, three kinds of Archimedean Copula functions were selected to construct the joint probability density function of variablesY, Z, and w. The reliability of the proposed method was 0.85%, which was higher than that of Nataf-MCS. The ratio γ of calculation workload was 35.71, indicating that the proposed method has superior computational efficiency and can reduce the calculation time.


## 5.3. Example 3

The soil parameters of rock and soil include liquidity index (LI), undrained shear strengthsu, remolded undrained shear strength sure, preconsolidation stress σp′, and vertical effective stress σv′. The experimental data in literature [53] were used as an example to illustrate the calculation steps of the proposed method. The performance function of the soil parameters gLI,su,sure,σp′,σv′=σp′sure+LIσp′-σv′-su2 was considered, and the reliability of rock and soil was analyzed under normal use limit state. Test data for soil parameters are shown in Table 11.Table 11
Test data for soil parameters.
Number
LI
s
u (kPa)
s
u
re (kPa)
σ
p
′ (kPa)
σ
v
′ (kPa)
1
1.28
2.58
0.22
7.94
9.56
2
1.27
4.74
0.16
29.41
12.82
3
1.45
5.42
0.39
17.14
16.44
4
1.20
5.82
0.59
27.63
20.06
5
1.26
6.97
0.65
23.37
24.04
⋯
⋯
⋯
⋯
⋯
⋯
343
2.11
12.26
1.53
80.24
61.90
344
2.05
14.10
1.53
86.36
68.02
345
1.57
15.20
1.84
94.76
75.66The meanμ, standard deviation σ, and edge distribution functions type of each variable were obtained by analyzing the experimental data as shown in Table 12.Table 12
Statistical parameters and distribution types of variables in example 3.
Items
LI
s
u (kPa)
s
u
re (kPa)
σ
p
′ (kPa)
σ
v
′ (kPa)
Meanμ
1.25
31.01
2.51
105.82
66.63
Standard deviationσ
0.61
29.48
3.81
103.21
53.48
Edge distributionfunction
Lognormal
Lognormal
Lognormal
Lognormal
LognormalFrom the calculation, we obtained the correlation coefficientρij of any two variables to form the matrix ρ.(23)ρ=LIsuresureσp′σv′LIsusureσp′σv′10.05-0.50-0.06-0.210.0510.170.840.57-0.500.1710.300.42-0.060.840.3010.73-0.210.570.420.731.Given the positive and negative correlations among variablesLI, su, sure, σp′, and σv′, we need to select the appropriate Copula function from Elliptical Copula function to construct the joint probability density function fLI,su,sure,σp′,σv′. Using AIC and BIC criteria, Gaussian Copula was selected for example 2. Afterwards, (3) was used to obtain any two variables related parameters θij and to form matrix θ.(24)θ=LIsuresureσp′σv′LIsusureσp′σv′10.07-0.91-0.08-0.260.0710.250.880.64-0.910.2510.410.53-0.080.880.4110.78-0.260.670.530.781.The 5D joint probability density function of variablesLI, su, sure, σp′, and σv′ can be expressed as the following equation:(25)fLI,su,sure,σp′,σv′=fLIfsufsurefσp′fσv′·cFLI,Fsu,Fsure,Fσp′,Fσv′;θ.The 5D joint probability density functionfLI,su,sure,σp′,σv′ and the performance function gLI,su,sure,σp′,σv′ were substituted into (13) to calculate the integrand.The structure of the integrand neural network could also be constructed. The integrand neural network had five inputs and single output. After adjustment, the number of hidden layer units wasm=5. The transfer functions of hidden layer and output layer were fs=ex and Purelin function, respectively.VariablesLI, su, sure, σp′, and σv′ are divided into 10 parts in the range of μi-4σi,μi+5σi. The five variables intersect each other to form the input sample points of the integrand network. The network output value y of the corresponding sample points was calculated by the integrand function. Given that each variable follows the lognormal distribution, we get μi-4σi=0.1 when the left interval of the variable μi-4σi≤0. The training sample set is shown in Table 13.Table 13
Training sample set of integrand network in example 3.
LI
s
u
s
u
re
y
σ
p
′
0.10
0.10
⋯
415.45
415.45
σ
v
′
0.10
25.32
201.85
227.07
0.10
0.10
0.10
1.072
8
e
-
006
5.355
0
e
-
007
⋯
4.194
5
e
-
011
2.767
2
e
-
011
0.10
0.10
1.64
3.355
9
e
-
007
1.681
9
e
-
007
⋯
1.389
7
e
-
011
9.205
3
e
-
012
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
0.10
0.10
13.94
3.854
6
e
-
009
1.994
8
e
-
009
⋯
2.527
4
e
-
013
1.728
5
e
-
013
⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯
3.08
119.45
0.10
1.330
3
e
-
010
6.429
6
e
-
011
⋯
5.826
7
e
-
014
3.721
8
e
-
014
3.08
119.45
1.64
3.859
7
e
-
011
1.872
9
e
-
011
⋯
1.790
5
e
-
014
1.148
3
e
-
014
⋯
⋯
⋯
⋯
⋯
⋯
⋯
3.08
119.45
13.94
2.426
6
e
-
013
1.215
9
e
-
013
⋯
1.782
3
e
-
016
1.180
2
e
-
016We used the Levenberg-Marquardt training algorithm. When the training reached 100 steps, the error of the integrand network reached the convergence state. The error convergence curve is shown in Figure6.Figure 6
Training error curve of integrand network in example 3.We used the same method for the previous example. The sample points of the original function network were obtained and are shown in Table14. The simulation values were substituted into (18) to obtain the reliability of the structure as shown in Table 15.Table 14
Simulative sample set of original function network in example 3.
LI
0.10
3.08
0.10
3.08
⋯
0.10
3.08
s
u
0.10
119.45
119.45
0.10
⋯
0.10
119.45
s
u
re
0.10
13.94
13.94
0.10
⋯
13.94
0.10
σ
p
′
0.10
415.45
415.45
0.10
⋯
0.10
415.45
σ
v
′
0.10
227.07
227.07
0.10
⋯
227.07
0.10Table 15
Example 3 reliability calculation results.
Method
Nataf-MCS
Direct integral method based on dual neural network
Gaussian Copula
Reliability
0.9594
0.9483
Relative error (%)
0
1.16The calculation workload of the two methods was substituted into (21), thus obtaining γ=9.90.The reliability of 5D correlated variables is analyzed using example 3. The reliability of the proposed method was extremely close to that of Nataf-MCS, and the relative error was only 1.16%. Fromγ=9.90, for the same problem, the workload of the two methods varied greatly. The neural network method has advantages in computational efficiency and reduces the calculation time for structural reliability analysis. This example shows that the proposed method can efficiently solve the reliability problem of multidimensional correlated variables.


## 6. Conclusion

Identifying the joint probability density function among variables is necessary to accurately calculate the structural reliability. This paper used the Copula function to construct the joint probability density function among correlated variables. According to Copula theory, the estimation of the edge probability density function among variables is independent of the selection of Copula function. The joint distribution function with arbitrary edge distribution and correlation structure could be constructed. This method was suitable for constructing joint probability density functions of 2D and multidimensional correlation problems. In addition, this paper proposed a Monte Carlo method based on Nataf inverse transformation. Nataf-MCS and the proposed method were used to solve the structural reliability among correlated variables. The proposed method could overcome the difficulty of solving multiple integral problems and achieve accurate calculation of reliability. However, we did not establish the relationship of the original function and every random variable because the training sample set originated from the integrand. We conveniently selected the sample set to improve the accuracy of calculating the reliability. Numerical examples showed that the proposed method is highly precise and efficient in calculating the reliability of multidimensional correlated variables.Copula function is often used to construct the static variable correlation model. In many applications, the correlation coefficient among variables has the characteristics of being dynamic and time-varying. Thus, a Copula model must be constructed to describe the effect of these factors on structural reliability. This issue must be further explored.In future work, the proposed method can also be extended for the reliability analysis of large and complex structural problems. In these cases, the workload of computing for dual neural network was negligible as compared with that of constructing network training samples. In addition, (21) can be simplified to γ≈Nm/Nn. Hence, the proposed method shows efficiency and advantages in analyzing large and complex structural problems.



---

*Source: 6976301-2017-12-11.xml*