Based on M-Copula Reliability Analysis of Random Load Correlation

Load is one of the main causes of structural failure, and the correlation among loads would affect the evaluation results of structural performance. &e purpose of this paper is to analyze the influence of the correlation among multiple loads on the structural reliability. In this paper, the nonparametric kernel density estimation (NKDE) method is used to estimate the probability density function (PDF) of related loads. In addition, themixed copula (M-Copula) model is proposed, which combines Gumbel copula, Frank copula, Clayton copula, and weight coefficient, and the model parameters are fitted by MATLAB software to get the correlation of related loads. &e reliability based on the related load combination is calculated according to the constructedmodel. After analyzing three numerical cases, the results show that the probability characteristics of NKDE estimation are very close to the actual conditions, and the reliability calculated by theM-Copula model is larger than those calculated by JCSS, Turkstra, and Gong methods. Using the M-Copula model for load correlation would avoid underestimating the reliability of the structure, which is conducive to structural economic development.


Introduction
With the diversified development of building structure forms, the functions of buildings are becoming more and more comprehensive, which leads to the diversified and complicated load forms in the service cycle of the structure [1][2][3]. For example, besides common loads, industrial buildings also bear fatigue, impact, additional loads for changing purposes and stacking loads for goods. [4][5][6]. According these conditions, it is necessary to put forward higher requirements for the load combination method and to consider the complex correlation among various loads in order to ensure economic design of structures and accurately evaluate reliability and safety of structures during service. In the past, random load combination methods mainly included Wen Yi-Kwei [7] coincidence method, upcrossing theory [8], and Kellishen, Ferry-Borges, JCSS rule, and Turkstra combination method commonly used [9]. ese methods or theories hold that all loads are independent and unrelated to each other. However, under influence of a certain condition, many loads may change significantly, and in general, there may be no clear correlation between the peak values of load components, but there are synchronous internal relations within the spatial neighborhood. e calculation based on independent variables is obviously unreasonable, which makes calculation results of structural reliability far from reality. At present, many scholars pay more and more attention to the study of load process correlation, such as Gong [9], Li [10], and Yao [11], who have made great contributions to development of load correlation. However, these methods or theories of Gong, Yao, and Li were all based on stationary binomial distribution or decompose loads into independent load pulses. And they emphasize correlation between adjacent loads and consider that correlation between random loads farther apart can be neglected. is idealized treatment is unreasonable, which will affect the structure of structural reliability analysis. e following is a brief introduction of the commonly used Turkstra method. e core of this method is to find out all variable load effects of structures and set them as the maximum value of the design reference period one by one.
Other components are combined according to random time value, and results of all the combinations are arranged to find the minimum value, which is taken as the design value of load effect [9]. Assuming that there are n variable loads participating in combination, the stochastic function of effect can be expressed as S i (t), t ∈ [0, T] , (i � 1,2, . . ., n). So, they have n maximum combinations and can be expressed as [12] S M1 � max 0≤t≤T S 1 (t) + S 2 t 0 + · · · + S n t 0 , S M2 � S 1 t 0 + max 0≤t≤T S 2 (t) + · · · + S n t 0 , where S Mi is maximum value of the i combination.
is combination method is simple in construction, which considered that variable loads are independent of each other. It only needs to list all the combination forms one by one and calculate the maximum value of each combination separately. However, this method is not necessarily the same as actual load combination, such as twin load [9,13]. e results of this calculation method are quite erroneous.
In recent studies, copula function has been widely used in constructing multidimensional joint distribution functions. In view of the problems of load correlations, this paper proposes a M-Copula model of multivariate complex loads.
e advantage of this model is that it can flexibly establish the joint distribution function of multivariate random loads and is not affected by the probability density function of single load.

Load Effect Combination Method Based on Copula Function
2.1. eory of Copula Function. After Sklar [14] proposed copula function in 1959, it was successfully applied in fields of statistical economics and finance [15], hydrology [16][17][18], and structural engineering [19][20][21][22][23][24], and a lot of research papers were published. e advantage of copula function is that it can directly construct joint distribution function for many related random variables or decouple a joint distribution function into multiple marginal distribution functions and a unique copula function [25]. ere are more than twenty copula functions commonly used, such as Gauss copula, Gumbel copula, Clayton copula, and Frank copula. e copula functions selected and constructed are different, and joint distribution functions are also different. According to Sklar theorem [14], for n consecutive marginal cumulative distribution functions F 1 , F 2 , . . . , F n , there is only one ndimensional copula function C(·), so that the following equation holds F x 1 , x 2 , . . . , x n � C F 1 x 1 , F 2 x 2 , . . . , F n x n , (2) where F(x 1 , x 2 , . . . , x n ) is the joint distribution function of the random vector x � x 1 , x 2 , . . . , x n , C(·)is the copula function, and F i (x i )is marginal distribution function of x i . erefore, seemingly unrelated single-load random variables can establish some relationship according to equation (2), and any change of random variables may affect results of combination. e joint probability density function of the random vector x � x 1 , x 2 , . . . , x n can be obtained by calculating the derivative of equation (2). e process is as follows: Equation (3) is simplified to where c(u 1 , u 2 , . .., u n ) � (z n C(u 1 , u 2 , ... , u n )/zu 1 zu 2 . ..z n ), In order to study relationship between loads, it is necessary to select a copula function.
e Gaussian copula function is a copula function of multivariate normal distribution. Its advantage is that it can construct correlation function between multiple random variables [14] and representation as where Φ R (·) is standard Gaussian distribution, R is correlation coefficient matrix, and φ − 1 (·) is inverse function of standard Gaussian distribution. When calculating random variables with the degree of freedom n, there exists n-fold integral in Gaussian copula function, which is inconvenient for digital calculation. Considering feasibility of calculation, Gumbel copula function, Frank copula function, and Clayton copula function belong to the Archimedean copulas family can be used. ey have obvious computing advantages: calculation process is relatively simple and high accuracy, so they are widely used in practical engineering. If there is a generator ϕ α (u) � − (− ln u) α , the Gumbel copula function can be deduced as C Gumbel α (u) or denoted as C Gumbel (u, α), and the calculation method be expressed as where α is the parameter of Gumbel copula function and α range of which is [1, ∞), and the change of α values will affect the correlation degree of random variables. If the value of α is 1, then the random variables u 1 , u 2 , . . . , u n are independent of each other; if the value range of α approaches +∞, then the random variables u 1 , u 2 , . . . , u n are completely correlated [14]. If there is a generator ϕ β (u) � − ln[(exp(− βu) − 1)/(exp(− β) − 1)], the Frank copula function can be deduced as C Frank β (u) or denoted as C Frank (u; β), and the calculation method be expressed as where β is the parameter of Frank copula function and β range of which is (− ∞, 0)∪(0, +∞), and change of β values will affect correlation degree of random variables. Similarly, if there is a generator ϕ c (u) � (1/c)(u − c − 1) and the range of c is (0, +∞), the Clayton copula function can be deduced as C Clayton c (u) or denoted as C Clayton (u; c), and the calculation method be expressed as Assuming the two-dimensional random variables (U 1 , U 2 ), the values of (α, β, c) are 1, 5, and 10, respectively, to generate 2000 random numbers. e scatter plot is shown in Figure 1. Figure 1 shows that different copula functions have their own distinct characteristics. Using a single copula function cannot accurately describe random characteristics of complex loads, so it is necessary to synthesize a variety of copula functions for analysis.

Copula Model of Load Combination.
Load action is an unstable process that changes with time, but the law of development is not uncorrelated [9][10][11]26]. e model of load effect combination can well reflect the correlation between loads. Total load effect is results of multiple loads. Most of loads have strong or weak correlation. Although the load combination model has obtained a lot of research results, most of them focus on the study of correlation of twodimensional loads. e combination model of three-dimensional loads and more than three-dimensional loads is not perfect, which is not conducive to development of load correlation research. e M-Copula function presented below can be used to establish a multiload correlation model, which provides a feasible method for the study of the correlation between multiloads.

Creating M-Copula Model.
According to the basic process of the Turkstra method, in the service period [0, T], load stochastic process is divided into k segments on average, and each calculation time is often τ; the length of j segment is τ j ; in any period of time interval τ j , the probability of load occurrence is p j , and then q j � 1 − p j . If the average score of load stochastic process is enough, then the probability density function of loads in [τ j − , τ j + ] time is consistent; that is, f In addition, another important assumption of the Turkstra method is that load effects in each time period are independent of each other. However, the load effect is a continuous function, and it is the result of the combination of various random loads, so the results of the Turkstra combination method are not accurate enough.
Copula function can solve problems of load dependence in each time period. Similarly, the discretization method is used. If there are n loads and they are correlated with each other, then the total load effect S can be expressed as Load S i (i � 1, 2, . . . , n) is also decomposed into k segments, and each calculation time is τ j . e specific process can be shown in Figure 2.
When load effect is discretized, the length of the j-th period is τ j , and the load effect is S i(τ j ) . en, load effect is composed of k load effect, which can be expressed as a function of k-dimensional random variables: Assuming that probability density function of the i load effect is f(s i ), C(u 1 , u 2 , . . . , u n ) can be calculated according to equation (4), but because of complexity of load correlation and diversity of load probability distribution, it is impossible to accurately simulate the correlation with a single copula function, so a hybrid copula function modeling method is proposed. In this paper, a new M-Copula load combination method is constructed by Gumbel copula function, Frank copula function, Clayton copula function, and some weights, which can be expressed by where C mix (·) is mixed copula function, (C G , C F , C C ) are expressed as Gumbel copula function, Frank copula function, and Clayton copula function, respectively, ω i is weight coefficient, and n i�1 ω i � 1, (α, β, c) are correlation degree of random variable u 1 , u 2 , . . . , u n . By establishing the M-Copula correlation model, the results of load effect combination can fully reflect relationship between loads and obtain more accurate reliability calculation results.

Nonparametric Kernel Density Function Estimation.
Load effect varies with time, and it is also a comprehensive response under various loads. If structures encounter ultimate loads during its service period, structures will be in a dangerous ultimate state. In this case, it is impossible to estimate probability density function of a single load. In other words, probability density function of a single load decoupled from presupposed probability density function under combined action of multiple loads will no longer be accurate. In this case, the nonparametric kernel density function estimation method [27][28][29][30] should be used to reconstruct the edge probability density function of a certain load. In order to analyze correlation of load random variables, the n-dimensional random variables (u 1 , u 2 , . . . , u n ) Mathematical Problems in Engineering 3 are analyzed, and combination equations (4) and (11) show that and If there are variable matrix u, when the number of n and h n is large enough, n and h n are correlated with each other, and h n > 0, then the kernel density estimate of u can be expressed as follows: where K(·) is kernel function, h n is bandwidth, and u � (u 1 , u 2 , . . . , u i ). Set u i is divided into k time periods and the j-th interval is τ j , then probability density function of the scorable quantity u i can be expressed as where Kernel function [K u 1 (·), K u 2 (·), . . . , K u n (·)] can generally be determined by the empirical method, and bandwidth (τ u 1j , τ u 2j , . . . , τ u nj ) can take different values.
. In order to ensure accuracy of results, the number of simulations k is as much as possible. In the case of the large sample, Gauss distribution can be chosen as kernel function of the kernel density estimation function and can be expressed as [31] where i � 1, 2, . . . , n and j � 1, 2, . . . , k. After obtaining kernel function, marginal probability density function . . , f u n (u n )] and marginal distribution function [F u 1 (u 1 ), F u 2 (u 2 ), . . . , F u n (u n )] of each variable of matrix variable u � (u 1 , u 2 , ..., u i ) can be obtained according to equation (15).  , x n,2 ), . . . , (x 1,k , x 2,k , . . . , x n,k )] are randomly generated samples, the joint distribution function of variables X can be expressed as H(X 1 , X 2 , . . . , X n ), and the corresponding copula function is C(X 1 , X 2 , . . . , X n ), then Kendall's rank correlation coefficient ρ of variable (X 1 , X 2 , . . . , X n ) can be expressed as [31] ρ

Multivariate Sequence
where ρ(X 1 , X 2 , . . . , X n ) are Kendall's rank correlation coefficient of the whole sample of variables X, ρ(X i , X j )are the Kendall's rank correlation coefficient of the sample of variables (X i , X j ), n is the degree of freedom of random variables, and k is the number of the random sample. According to the simple algorithm, the calculation method of Kendall's rank correlation coefficient ρ(X i , X k ) of any two random variables can be as follows: where ρ(X i , X υ ) are Kendall's rank correlation coefficient of any two variables, and then 1 ≤ i ≤ υ ≤ n and 1 ≤ ξ ≤ η ≤ k; among them, e Kendall's rank correlation coefficients of any two variables can be obtained by equation (18), and sample estimates of Kendall's rank correlation coefficients ρ of population variables can be calculated by equation (16), from which ρ(X 1 , X 2 , . . . , X n ) � ρ(X 1 , X 2 , . . . , X n ) can be obtained approximately. According to Table 1, the Kendall's rank correlation coefficient and copula function parameters can be transformed. e weight coefficient ω i is estimated by the maximum likelihood estimation method [24]. Based on the M-Copula load correlation analysis model, the specific process is divided into the following steps: Mathematical Problems in Engineering 5 (4) By using the maximum likelihood estimation method, weight parameters (ω 1 , ω 2 , ω 3 ) are estimated, and load correlation analysis model C(α, β, c; ω 1 , ω 2 , ω 3 ) based on M-Copula function is obtained. erefore, it is expressed in Figure 3.

Reliability Calculation Model
Based on M-Copula e traditional reliability calculation method is probability expression of resistance and load, which is independent when dealing with loads and does not consider the complex correlation between loads, so failure probability obtained by the general calculation method is inconsistent with actual situation. Considering the correlation of various loads will ultimately affect reliable calculation results. Different results of load combination will affect results of reliability calculation. e concept of reliability is probability of completing a predetermined function within a specified time under specified conditions. In order to obtain reliable probability, it is necessary to calculate load effect and resistance. If the structure is in a good working environment, it can be approximated that resistance of the structure will not change during service, but this assumption is unreasonable, because structures will suffer uncertain damages during service, which will greatly affect the performance of the structure.
Similarly to the principle of Figure 2, the resistance of structures can be decomposed into k time intervals during the whole service period [0, T], and resistance value in the jth period is R j , as shown in Figure 4.
If the number of τ j or k is large enough, it can be considered that R j(τ i ) does not change in the time range of τ j , and then the reliability model can be expressed as where Z is performance function, S τ i is load effect in the time range of τ i , and min(·) is small value function. P f represents failure probability of structures. For n-related loads, the failure probability of structures can be expressed as where (α, β, c) are parameters of copula function and (ω 1 , ω 2 , ω 3 ) are weight coefficient.

Case Analysis
e actual loads suffered by structures usually have many forms, such as Gauss distribution (GD) and Gumbel distribution (Gumbel Ι). In order to reflect correlation analysis of multidimensional correlated loads, the three-dimensional load correlation is selected to study. en, the M-Copula model can be expressed as Next, through analysis of specific statistical dates, the confidence level is 95%. Figure 5 shows sample dates of actual statistical results of loads (S 1 , S 2 , S 3 ) and results of NKDE fitting.
From the results of Figure 5, it is found that the three load effects are not completely random but conform to some form of distribution. e fitting distribution parameters are shown in Table 2.
By comparing the results of the above NKDE estimation method with actual data, it is found that probability density function obtained by the nonparametric kernel density estimation method has a good accuracy, which will be conducive to calculation of reliability. Common variable loads are usually distributed in extreme form. Load effects are combined on basis of Ferry-Borges, JCSS, and Turkstra combination methods commonly used in the past. e basic core of these methods is to determine maximum value in the time interval and then take the basic time point value as combination value or use each time interval. e maximum value in the interval is combined, which is obviously too safe. At the same time, it will lead to the high reliability design value of structure, which is not conducive to economic design of structures.
However, probabilistic density function of nonparametric kernel density estimation can weaken this effect from each component of load effect combination, and calculation results are reliable. According to the same method, other two groups of data samples are analyzed, and the distribution parameters of three groups of data are finally expressed in Table 3.
According to equation (16) and Table 1, the parameters of Gumbel copula, Frank copula, and Clayton copula functions can be estimated by using MATLAB R2014b programming.
e estimated results are expressed in Table 4.
According to the results of parameter estimation, the M-Copula model established in this paper is used to calculate structural reliability index R el , and the results are compared with those calculated by JCSS, Turkstra, and Gong methods. Comparison results are shown in Table 5. Calculation results in Table 5 show that correlation between complex loads will affect the calculation results of reliability and even change evaluation of structural performance. It is pointed out in the table that traditional JCSS and Turkstra methods are obviously conservative in calculating reliability evaluation indicators, which will inevitably lead to uneconomical treatment schemes. e M-Copula function model presented in this paper can well describe relationship between multiple complex loads, accurately calculate reliability index, and provide an accurate and scientific method for structural performance evaluation. 6 Mathematical Problems in Engineering

Results and Discussion
ere are many kinds of loads that contain relevance is supported during structural service, in which the load correlation would affect the results of structural reliability evaluation. e purpose of this study is to analyze the influence of load correlation on reliability. In the past research, the assumption is that the correlation among loads is always ignored, or the correlation between the two loads is considered, and it would underestimate the result of reliability evaluation and is not conducive to the balance between reliability and economy. In this paper, firstly, the probability characteristics of three groups of loads (S 1 , S 2 , S 3 ) are estimated by the NKDE method (the confidence level is 95%). Figure 5 shows that the accuracy of NKDE fitting is very high according to the results of estimation and statistics. Secondly, the M-Copula model is proposed to build the correlation of multiple loads, and MATLAB software is used to fit the model parameters. In three numerical case studies, the M-Copula model is used to calculate reliability of different load correlation, and the results are compared with those of the JCSS, Turkstra, and Gong method. Table 5 shows that the reliability calculated by the model is the highest of the four methods. It is also found in the comparative analysis that the reliability index calculated by the M-Copula model is 24.79%, 23.81%, and 24.05% higher than that calculated by the JCSS method in the international standard.
It can be seen that the past research methods underestimate the actual reliability, resulting in a large amount of waste of resources and adverse economic development of the structure.
is study analyzes the internal relationship among loads and provides a feasible mathematical model for the economic development of reliability design. e model can be used to analyze the correlation of loads greater than 3. However, a large number of load correlation analysis must involve multiple integration, which requires computer-aided calculation.
In this paper, the static load correlation analysis is solved. Firstly, the method of nonparametric kernel density estimation (NKDE) is used to estimate the probability density function (PDF), and a mix copula (M-Copula) function model is proposed. Finally, the correlation of related loads is obtained through MATLAB model parameter fitting. However, the dynamic load analysis of structural dynamics has some limitations, so we will carry out the study of dynamic load correlation in the future.

Conclusions
e M-Copula function model established in this paper and analysis results of specific cases can be summarized as follows: (1) Load effect results from the combined action of multiple complex loads, which are not completely independent, and most of the loads are correlated. Previous studies have only considered that loads are independent of each other, not related to each other, or that correlation between two loads is considered; obviously, these are unreasonable.
(2) Loads are random variables, and in fact they cannot be completely decoupled and independent. e NKDE method can be used to fit three kinds of related loads, and results show that the method has high accuracy.
(3) Establishing the M-Copula function model can accurately describe the relationship between loads and calculate correlation between loads through Gumbel copula, Frank copula, Clayton copula function, and comprehensive model of weight coefficient. By analyzing actual data of three kinds of related loads, the results show that M-Copula function has a high consistency in estimating correlation between multivariate loads. (4) Comparing the results of the reliability index calculated by the M-Copula model with those calculated by the traditional JCSS method and Turkstra method, it is found that the calculation method used Note: S i is the load effect value, R is the bearing capacity, μ and σ are mean and variance, respectively, and δ is coefficient of variation.  in the past is too conservative, which is not conducive to the result of structural performance evaluation and economy of later structural treatment, while the M-Copula model accurately calculates load effects, and current performance should be accurately evaluated by combining the values and structural reliability indicators.

Data Availability
All data generated or analyzed during this study are included in this article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.