Mathematical Problems in Engineering
Volume 2017 (2017), Article ID 6976301, 12 pages
https://doi.org/10.1155/2017/6976301

Research Article

The Method of Solving Structural Reliability with Multiparameter Correlation Problem

Juan Du, Haibin Li, and Yun He

College of Science, Inner Mongolia University of Technology, Hohhot 010051, China

Correspondence should be addressed to Haibin Li; lhbnm2003@126.com

Received 1 July 2017; Revised 12 November 2017; Accepted 16 November 2017; Published 11 December 2017

Academic Editor: Fiorenzo A. Fazzolari

Copyright © 2017 Juan Du et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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 -dimensional cases, the Copula function was defined as the -dimensional joint distribution function. The edge distribution was in the n space [9].

According to Sklar theorem [10], the joint distribution function of variables , could be expressed as

where is the edge distribution function of the variable : . is the Copula function, and is the correlation parameter of the Copula function.

From (1), the joint probability density function could be obtained as

where is the edge probability density function of variables and is the density function for .

If the edge distribution function of variables 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 variable was an Ellipse Copula function, the number of correlation parameters was the same as that of the correlation coefficients among variables and had one-to-one correspondence. The correlation parameter of any two variables and was . The relationship between and Pearson linear correlation coefficient was [9]

When the Copula function of the multidimensional variable 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 . Maximum likelihood estimation is generally used to solve the correlation parameter [44].

2.3. Selection of the Optimal Copula Function

Different Copula functions could describe different correlations among variables. Table 1 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]:

where is the number of correlation parameters . In the Ellipse Copula function, , whereas in the Archimedean Copula function, .

Table 1: Copula function types.

Copula types
Copula distribution function

range of values

Gaussian

Gumbel

Clayton

Frank

Number 16

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 variables were , 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 of the standard normal distribution variables was solved.

Ang et al. established the following empirical equation [49] for various distribution types under normal circumstances :

where the coefficient is a function of the correlation coefficient and the edge distribution . The equation for the coefficient was provided in literature [51].

(2) The standard normal distribution variables were correlated, and their correlation coefficient matrices were . was decomposed by Gorensky to obtain the lower triangular matrix .

(3) The traditional random number generation method was applied to produce independent variables , which all followed the standard normal distribution.

(4) Based on the lower triangular matrix , the independent normal random variables were transformed into standard normal distribution variables .

(5) The relevant standard normal distribution variables were mapped to the original spatial variables using the principle of equal probability transformation.

(6) The original spatial variables were substituted into the functional functions of the reliability analysis, and the structure response was calculated: . The number of times was the number of failure sample points . The unbiased estimate of the structural failure probability was [49] as follows:

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

Structural reliability could be expressed as the probability of function . could be expressed as a function of basic variables ; that is, . The joint probability density function of the correlated variables was , introducing a weight function:

Structural reliability could be expressed as multiple integrals as follows:

In (16), and were the mean and standard deviation, respectively, of random variable . and were positive integers, which were determined by accuracy requirements. In general, and .

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 variable and the original function . 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

The derivation of (17) could be given by

where . Equation (16) could be rewritten as a function of the relationship between the output and input variables as follows:

The neural network completed the mapping relationship between the input variable and the integrand . The network structure is shown in Figure 2.

Figure 2: Integrand network structure.

Let the original function of the integrand be . According to multiple integral theory, function is the weighted algebraic sum of each vertex in a hypercube and could be expressed as

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, inputs, single output, and unit in hidden layer. In the integrand network, the connection weights and threshold of input layer to hidden layer unit were and , respectively; the connection weights and threshold of hidden layer to output layer were and 0, respectively; and the activation function of the hidden layer units was . In the original function network, the connection weights and threshold of input layer to hidden layer unit were and , respectively; the connection weights and threshold of hidden layer to output layer were and , respectively; and the activation function of the hidden layer units was . 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 as , where is the lognormal distribution random variables with mean value , standard deviation , and coefficient of variation . was an extreme value distribution of type I with mean value , standard deviation , and coefficient of variation . The correlation coefficient of and was . The structural reliability was calculated.

Based on the Copula function, the 2D joint probability density function of the random variables and was established and could be written as

where

where , , , and .

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 variables and . 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

and 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 were and Purelin function.

The ranges of and and were divided into 50 parts. Two variables crossed each other to form the network input samples. The network output value 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.

2.4920
2.9479
3.4038
⋯
24.3741
24.8300

0
0
0
⋯
0
0

0
0
0
⋯
0
0

⋯
⋯
⋯
⋯
⋯
⋯
⋯

42.8520

⋯
0
0

44.4847

⋯
0
0

46.1173

⋯
0
0

47.7500

⋯
0
0

This example used Levenberg-Marquardt training algorithm. Error convergence curve of training integrand network after 1000 steps is shown in Figure 3.

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 Table 4. The results were substituted into (18) to obtain the result of reliability.

Table 4: Simulative sample set of original function network in example 1.

24.8300
24.8300
2.4920
2.4920

47.7500

For 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 Table 5.

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.00

The ratio of the calculation workloads of the Nataf-MCS and the proposed method can be obtained from the following equation:

where is a single workload for structural analysis. is the calculation times of the Nataf-MCS. is the calculation times required to construct the network training sample. is a single workload for neural network training. 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, . From the calculation in this example, .

From Table 5, 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 , 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 Figure 4, Y is the horizontal load on the beam, is the lateral load on the beam, and represents the beam width. The performance function of structure . 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 , , and were obtained and are listed in Table 6.

Table 6: Statistical parameters and distribution types of variables in example 2.

Project
(lb)
(lb)
(in)

Mean
988.72
506.28
1.98

Standard deviation
120.78
132.48
0.27

Edge distribution function
Normal
Normal
Normal

Calculated 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 variables Y, Z, and .

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.05

The 3D joint probability density function and the performance function 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 is . The transfer functions of hidden layer and output layer were and Purelin function, respectively.

Variables , , and were separately divided into 30 parts in the range of . The three variables intersect each other to form the integrand network input sample point. The network output value 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.

1.17
1.23
⋯
2.74
2.80

626.40
108.84

⋯

626.40
136.25

⋯

⋯
⋯
⋯
⋯
⋯
⋯
⋯

626.40
903.70

⋯

⋯
⋯
⋯
⋯
⋯
⋯
⋯

1351.10
108.84

⋯

1351.10
136.25

⋯

⋯
⋯
⋯
⋯
⋯
⋯
⋯

1351.10
903.70

⋯

We 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 Figure 5.

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 Table 9. The results were substituted into (18) to obtain the result of reliability.

Table 9: Simulative sample set of original function network in example 2.

626.40
1351.10
626.40
1351.10
626.40
1351.10
626.40
1351.10

108.84
903.70
108.84
108.84
903.70
108.84
903.70
903.70

1.17
2.80
2.80
2.80
1.17
1.17
2.80
1.17

Nataf-MCS was introduced in Section 3. The sampling was 1,000,000 times, and the reliability of example 2 was 0.9998. In this example, . 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.62

In this example, three kinds of Archimedean Copula functions were selected to construct the joint probability density function of variables , , and . 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 strength , remolded undrained shear strength , preconsolidation stress , and vertical effective stress . 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 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
(kPa)
(kPa)
(kPa)
(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.66

The 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
(kPa)
(kPa)
(kPa)
(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 distribution
function
Lognormal
Lognormal
Lognormal
Lognormal
Lognormal

From the calculation, we obtained the correlation coefficient of any two variables to form the matrix .

Given the positive and negative correlations among variables , , , , and , we need to select the appropriate Copula function from Elliptical Copula function to construct the joint probability density function . Using AIC and BIC criteria, Gaussian Copula was selected for example 2. Afterwards, (3) was used to obtain any two variables related parameters and to form matrix .

The 5D joint probability density function of variables , , , , and can be expressed as the following equation:

The 5D joint probability density function and the performance function 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 was . The transfer functions of hidden layer and output layer were and Purelin function, respectively.

Variables , , , , and are divided into 10 parts in the range of . The five variables intersect each other to form the input sample points of the integrand network. The network output value of the corresponding sample points was calculated by the integrand function. Given that each variable follows the lognormal distribution, we get when the left interval of the variable . The training sample set is shown in Table 13.

Table 13: Training sample set of integrand network in example 3.

LI

0.10
0.10
⋯
415.45
415.45

0.10
25.32

201.85
227.07

0.10
0.10
0.10

⋯

0.10
0.10
1.64

⋯

⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯

0.10
0.10
13.94

⋯

⋯
⋯
⋯
⋯
⋯
⋯
⋯
⋯

3.08
119.45
0.10

⋯

3.08
119.45
1.64

⋯

⋯
⋯

⋯
⋯
⋯
⋯
⋯

3.08
119.45
13.94

⋯

We 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 Figure 6.

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 Table 14. 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

0.10
119.45
119.45
0.10
⋯
0.10
119.45

0.10
13.94
13.94
0.10
⋯
13.94
0.10

0.10
415.45
415.45
0.10
⋯
0.10
415.45

0.10
227.07
227.07
0.10
⋯
227.07
0.10

Table 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.16

The calculation workload of the two methods was substituted into (21), thus obtaining .

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 , 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 . Hence, the proposed method shows efficiency and advantages in analyzing large and complex structural problems.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant no. 11262014 and Graduate Student Scientific Research Innovation Project of Autonomous Region under Grant no. B20171012805.

References

K. Goda, “Statistical modeling of joint probability distribution using copula: Application to peak and permanent displacement seismic demands,” Structural Safety, vol. 32, no. 2, pp. 112–123, 2010.

B. J. Leira, “Probabilistic assessment of weld fatigue damage for a nonlinear combination of correlated stress components,” Probabilistic Engineering Mechanics, vol. 26, no. 3, pp. 492–500, 2011.

D. Li, Y. Chen, W. Lu, and C. Zhou, “Stochastic response surface method for reliability analysis of rock slopes involving correlated non-normal variables,” Computers & Geosciences, vol. 38, no. 1, pp. 58–68, 2011.

D.-Q. Li, X.-S. Tang, K.-K. Phoon, Y.-F. Chen, and C.-B. Zhou, “Bivariate simulation using copula and its application to probabilistic pile settlement analysis,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 37, no. 6, pp. 597–617, 2013.

H. F. Burcharth, J. D. Sorensen, and E. Christian, “On the evaluation of failure probability of monolithic vertical wall breakwaters,” Proceedings of Wave Barriers in Deep Waters, pp. 458–468, 1994.

A. Der Kiureghian, H.-Z. Lin, and S.-J. Hwang, “Second-order reliability approximations,” Journal of Engineering Mechanics, vol. 113, no. 8, pp. 1208–1225, 1987.

S. Yue, “The bivariate lognormal distribution for describing joint statistical properties of a multivariate storm event,” Environmetrics, vol. 13, no. 8, pp. 811–819, 2002.

A. C. Rencher, Methods of Multivariate Analysis, John Wiley & Sons, Inc., New York, NY, USA, 1995.

R. B. Nelsen, An Introduction to Copulas, Springer, New York, NY, USA, 2nd edition, 2006.

M. Sklar, “Fonctions de repartition an dimensions et leurs marges,” Publications de l'Institut de Statistique de l'Université de Paris, vol. 8, pp. 229–231, 1959.

X. Zeng, J. Ren, Z. Wang, S. Marshall, and T. Durrani, “Copulas for statistical signal processing (Part I): Extensions and generalization,” Signal Processing, vol. 94, no. 1, pp. 691–702, 2014.

Y. Li, B. Shi, and J. Zhang, “The copula joint function and its application in probabilistic seismic hazard analysis,” Earthquake Science, vol. 30, no. 3, pp. 292–301, 2008.

X.-S. Tang, D.-Q. Li, C.-B. Zhou, and K.-K. Phoon, “Probabilistic analysis of load-displacement hyperbolic curves of single pile using Copula,” Rock and Soil Mechanics, vol. 33, no. 1, pp. 171–178, 2012.

D. Huang, C. Yang, B. Zeng, and G. Fu, “A copula-based method for estimating shear strength parameters of rock mass,” Mathematical Problems in Engineering, vol. 2014, Article ID 693062, 11 pages, 2014.

Z. Liu, X. Ma, J. Yang, and Y. Zhao, “Reliability modeling for systems with multiple degradation processes using inverse Gaussian process and copulas,” Mathematical Problems in Engineering, vol. 2014, Article ID 829597, 2014.

Y. Liu and X. Fan, “Time-independent reliability analysis of bridge system based on mixed copula models,” Mathematical Problems in Engineering, vol. 2016, Article ID 2720614, 13 pages, 2016.

X. Xu, J. Li, J. Gong, H. Deng, and L. Wan, “Copula-Based Slope Reliability Analysis Using the Failure Domain Defined by the g-Line,” Mathematical Problems in Engineering, vol. 2016, Article ID 6141838, 11 pages, 2016.

I. Lee, K. K. Choi, and L. Zhao, “Sampling-based RBDO using the stochastic sensitivity analysis and dynamic kriging method,” Structural and Multidisciplinary Optimization, vol. 44, no. 3, pp. 299–317, 2011.

I. Lee, K. K. Choi, Y. Noh, L. Zhao, and D. Gorsich, “Sampling-based stochastic sensitivity analysis using score functions for RBDO problems with correlated random variables,” Journal of Mechanical Design, vol. 133, no. 2, Article ID 021003, 2011.

Y. Noh, K. K. Choi, and L. Du, “Reliability-based design optimization of problems with correlated input variables using a Gaussian Copula,” Structural and Multidisciplinary Optimization, vol. 38, no. 1, pp. 1–16, 2009.

Y. Noh, K. K. Choi, and I. Lee, “Identification of marginal and joint CDFs using Bayesian method for RBDO,” Structural and Multidisciplinary Optimization, vol. 40, no. 1-6, pp. 35–51, 2010.

Y. Noh, K. K. Choi, I. Lee, D. Gorsich, and D. Lamb, “Reliability-based design optimization with confidence level for non-Gaussian distributions using bootstrap method,” Journal of Mechanical Design, vol. 133, no. 9, Article ID 091001, 2011.

Y. Noh, Input Model Uncertainty and Reliability-Based Design Optimization with Associated Confidence Level, The University of Iowa, Iowa, USA, 2009.

H. Cho, K. K. Choi, and D. Lamb, “Sensitivity developments for RBDO with dependent input variable and varying input standard deviation,” Journal of Mechanical Design, vol. 139, no. 7, 9 pages, 2017.

Y. Noh, K. K. Choi, and I. Lee, “Reduction of ordering effect in reliability-based design optimization using dimension reduction method,” AIAA Journal, vol. 47, no. 4, pp. 994–1004, 2009.

Y. Noh, K. K. Choi, I. Lee, D. Gorsich, and D. Lamb, “Reliability-based design optimization with confidence level under input model uncertainty due to limited test data,” Structural and Multidisciplinary Optimization, vol. 43, no. 4, pp. 443–458, 2011.

Y. Noh, K. K. Choi, and I. Lee, “Comparison study between MCMC-based and weight-based Bayesian methods for identification of joint distribution,” Structural and Multidisciplinary Optimization, vol. 42, no. 6, pp. 823–833, 2010.

H. Cho, K. K. Choi, I. Lee, and D. Lamb, “Design Sensitivity Method for Sampling-Based RBDO with Varying Standard Deviation,” Journal of Mechanical Design, vol. 138, no. 1, Article ID 011405, 2016.

Y.-G. Zhao and T. Ono, “A general procedure for first/second-order reliability method (FORM/SORM),” Structural Safety, vol. 21, no. 2, pp. 95–112, 1999.

Y.-G. Zhao and T. Ono, “Third-moment standardization for structural reliability analysis,” Journal of Structural Engineering, vol. 126, no. 6, pp. 724–731, 2000.

Y.-G. Zhao and T. Ono, “Moment methods for structural reliability,” Structural Safety, vol. 23, no. 1, pp. 47–75, 2001.

Q. Gong, J. Zhang, C. Tan, and C. Wang, “Neural networks combined with importance sampling techniques for reliability evaluation of explosive initiating device,” Chinese Journal of Aeronautics, vol. 25, no. 2, pp. 208–215, 2012.

V. Papadopoulos, D. G. Giovanis, N. D. Lagaros, and M. Papadrakakis, “Accelerated subset simulation with neural networks for reliability analysis,” Computer Methods Applied Mechanics and Engineering, vol. 223/224, pp. 70–80, 2012.

R. Douc, F. Maire, and J. Olsson, “On the use of Markov chain MONte Carlo methods for the sampling of mixture models: a statistical perspective,” Statistics and Computing, vol. 25, no. 1, pp. 95–110, 2015.

A. Genz and A. Malik, “Remarks on algorithm 006: An adaptive algorithm for numerical integration over an N-dimensional rectangular region,” Journal of Computational and Applied Mathematics, vol. 6, no. 4, pp. 295–299, 1980.

M. Papadrakakis, V. Papadopoulos, and N. D. Lagaros, “Structural reliability analyis of elastic-plastic structures using neural networks and Monte Carlo simulation,” Computer Methods Applied Mechanics and Engineering, vol. 136, no. 1-2, pp. 145–163, 1996.

P. A. M. Lopes, H. M. Gomes, and A. M. Awruch, “Reliability analysis of laminated composite structures using finite elements and neural networks,” Composite Structures, vol. 92, no. 7, pp. 1603–1613, 2010.

Y.-L. Zuo, H.-H. Zhu, and X.-J. Li, “An ANN-based four order moments method for geotechnical engineering reliability analysis,” Rock and Soil Mechanics, vol. 34, no. 2, pp. 513–518, 2013.

J. Cheng and Q. S. Li, “Reliability analysis of structures using artificial neural network based genetic algorithms,” Computer Methods Applied Mechanics and Engineering, vol. 197, no. 45-48, pp. 3742–3750, 2008.

W. Meng, Z. Yang, X. Qi, and J. Cai, “Reliability analysis-based numerical calculation of metal structure of bridge crane,” Mathematical Problems in Engineering, vol. 2013, Article ID 260976, 5 pages, 2013.

A. H. Elhewy, E. Mesbahi, and Y. Pu, “Reliability analysis of structures using neural network method,” Probabilistic Engineering Mechanics, vol. 21, no. 1, pp. 44–53, 2006.

H. B. Li, Y. He, and X. B. Nie, “Structural reliability calculation method based on the dual neural network and direct integration method,” Neural Computing and Applications, pp. 1–9, 2016.

H. Dai and Z. Cao, “A Wavelet Support Vector Machine-Based Neural Network Metamodel for Structural Reliability Assessment,” Computer-Aided Civil and Infrastructure Engineering, vol. 32, no. 4, pp. 344–357, 2017.

U. Cherubini, E. Luciano, and W. Vecchiato, “Copula Methods in Finance,” Copula Methods in Finance, pp. 1–293, 2013.

H. Akaike, “A New Look at the Statistical Model Identification,” IEEE Transactions on Automatic Control, vol. 19, no. 6, pp. 716–723, 1974.

G. Schwarz, “Estimating the dimension of a model,” The Annals of Statistics, vol. 6, no. 2, pp. 461–464, 1978.

Y. Hao, X. Rong, L. Ma, P. Fan, and H. Lu, “Uncertainty analysis on risk assessment of water inrush in karst tunnels,” Mathematical Problems in Engineering, vol. 2016, Article ID 2947628, 2016.

S. Mishra, O. A. Vanli, B. P. Alduse, and S. Jung, “Hurricane loss estimation in wood-frame buildings using Bayesian model updating: Assessing uncertainty in fragility and reliability analyses,” Engineering Structures, vol. 135, pp. 81–94, 2017.

A. Ang, -S. Tang W, and -S. Tang, Probability Concepts in Engineering Planning and Design, John Wiley and Sons, New York, NY, USA, 1984.

A. Nataf, “Determination des distributions de probabilites dont les marges sont donnees,” Comptes Rendus De L'academie Des Sciences, vol. 255, pp. 42-43, 1962.

A. Der Kiureghian and P.-L. Liu, “Structural reliability under incomplete probability information,” Journal of Engineering Mechanics, vol. 112, no. 1, pp. 85–104, 1986.

C. Jiang, W. Zhang, B. Wang, and X. Han, “Structural reliability analysis using a copula-function-based evidence theory model,” Computers & Structures, vol. 143, pp. 19–31, 2014.

J. Ching and K. Phoon, “Corrigendum: Modeling parameters of structured clays as a multivariate normal distribution,” Canadian Geotechnical Journal, vol. 49, no. 12, pp. 1447–1450, 2012.