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. The 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, the mixed 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. The reliability based on the related load combination is calculated according to the constructed model. 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 the M-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.
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–3]. For example, besides common loads, industrial buildings also bear fatigue, impact, additional loads for changing purposes and stacking loads for goods. [4–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  coincidence method, upcrossing theory , and Kellishen, Ferry-Borges, JCSS rule, and Turkstra combination method commonly used . These 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. The 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 , Li , and Yao , 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. This idealized treatment is unreasonable, which will affect the structure of structural reliability analysis.
The following is a brief introduction of the commonly used Turkstra method. The 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 . Assuming that there are n variable loads participating in combination, the stochastic function of effect can be expressed as , (i = 1,2, …, n). So, they have n maximum combinations and can be expressed as where is maximum value of the i combination. This 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]. The 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. The 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.
2. Load Effect Combination Method Based on Copula Function
2.1. Theory of Copula Function
After Sklar  proposed copula function in 1959, it was successfully applied in fields of statistical economics and finance , hydrology [16–18], and structural engineering [19–24], and a lot of research papers were published. The 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 . There are more than twenty copula functions commonly used, such as Gauss copula, Gumbel copula, Clayton copula, and Frank copula. The copula functions selected and constructed are different, and joint distribution functions are also different. According to Sklar theorem , for n consecutive marginal cumulative distribution functions , there is only one n-dimensional copula function , so that the following equation holdswhere is the joint distribution function of the random vector , is the copula function, and is marginal distribution function of . Therefore, 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. The joint probability density function of the random vector can be obtained by calculating the derivative of equation (2). The process is as follows:
Equation (3) is simplified towhere , , is marginal probability density function of , and is the joint probability density function of .
In order to study relationship between loads, it is necessary to select a copula function. The Gaussian copula function is a copula function of multivariate normal distribution. Its advantage is that it can construct correlation function between multiple random variables  and representation aswhere is standard Gaussian distribution, is correlation coefficient matrix, and 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. They 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 , the Gumbel copula function can be deduced as or denoted as , and the calculation method be expressed aswhere is the parameter of Gumbel copula function and range of which is , and the change of values will affect the correlation degree of random variables. If the value of is 1, then the random variables are independent of each other; if the value range of approaches , then the random variables are completely correlated .
If there is a generator , the Frank copula function can be deduced as or denoted as , and the calculation method be expressed aswhere is the parameter of Frank copula function and range of which is , and change of values will affect correlation degree of random variables. Similarly, if there is a generator and the range of is , the Clayton copula function can be deduced as or denoted as , and the calculation method be expressed as
Assuming the two-dimensional random variables , the values of are 1, 5, and 10, respectively, to generate 2000 random numbers. The 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.
2.2. Copula Model of Load Combination
Load action is an unstable process that changes with time, but the law of development is not uncorrelated [9–11, 26]. The 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 two-dimensional loads. The 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. The 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.
2.2.1. Creating M-Copula Model
According to the basic process of the Turkstra method, in the service period , load stochastic process is divided into k segments on average, and each calculation time is often ; the length of j segment is ; in any period of time interval , the probability of load occurrence is , and then . If the average score of load stochastic process is enough, then the probability density function of loads in time is consistent; that is, . 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 is also decomposed into k segments, and each calculation time is . The specific process can be shown in Figure 2.
When load effect is discretized, the length of the j-th period is , and the load effect is . Then, 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 , 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 bywhere is mixed copula function, are expressed as Gumbel copula function, Frank copula function, and Clayton copula function, respectively, is weight coefficient, and , are correlation degree of random variable . 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.
2.2.2. 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–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 are analyzed, and combination equations (4) and (11) show that and . If there are variable matrix , when the number of n and is large enough, n and are correlated with each other, and , then the kernel density estimate of can be expressed as follows:where is kernel function, is bandwidth, and . Set is divided into k time periods and the j-th interval is , then probability density function of the scorable quantity can be expressed aswhere Kernel function can generally be determined by the empirical method, and bandwidth can take different values. When or , . 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 where and . After obtaining kernel function, marginal probability density function and marginal distribution function of each variable of matrix variable can be obtained according to equation (15).
2.2.3. Multivariate Sequence Correlation Analysis
If there are n variables , then and are the i-1 and i components of variable X. Assuming that marginal distribution function of each component is unchanged, it can be expressed as . If are randomly generated samples, the joint distribution function of variables X can be expressed as , and the corresponding copula function is , then Kendall’s rank correlation coefficient of variable can be expressed as where are Kendall’s rank correlation coefficient of the whole sample of variables X, are the Kendall’s rank correlation coefficient of the sample of variables , 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 of any two random variables can be as follows:where are Kendall’s rank correlation coefficient of any two variables, and then and ; among them,
The 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 can be obtained approximately. According to Table 1, the Kendall’s rank correlation coefficient and copula function parameters can be transformed.
The weight coefficient is estimated by the maximum likelihood estimation method . Based on the M-Copula load correlation analysis model, the specific process is divided into the following steps:(1)Estimate the number of relevant loads of structures in the service period , divide it into n time intervals, the i time interval is , and .(2)According to the sample dates of random load in the interval of each period and denoted , .(3)The Kendall’s rank correlation coefficients of any two kinds of related loads are calculated, and then the Kendall’s rank correlation coefficients and the parameters of copula function are estimated.(4)By using the maximum likelihood estimation method, weight parameters are estimated, and load correlation analysis model based on M-Copula function is obtained. Therefore, it is expressed in Figure 3.
3. Reliability Calculation Model Based on M-Copula
The 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. The 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 , and resistance value in the j-th period is , as shown in Figure 4.
If the number of or k is large enough, it can be considered that does not change in the time range of , and then the reliability model can be expressed aswhere Z is performance function, is load effect in the time range of , and is small value function. represents failure probability of structures. For n-related loads, the failure probability of structures can be expressed aswhere are parameters of copula function and are weight coefficient.
4. Case Analysis
The 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. Then, 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 and results of NKDE fitting.
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. The 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. The 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. The 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 Rel, 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. The 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.
5. Results and Discussion
There 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. The 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 (S1, S2, S3) 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. This study analyzes the internal relationship among loads and provides a feasible mathematical model for the economic development of reliability design. The 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.
The 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. The 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 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.
|:||Joint distribution function|
|:||Joint probability density function|
|:||Standard Gauss distribution function|
|:||Generators of Gumbel, Frank, and Clayton function|
|:||Parameters of copula function|
|:||Copula function of Gumbel, Frank, and Clayton|
|:||Mix copula function|
|:||Kendall’s correlation coefficient|
|:||Coefficient of variation|
|Rel:||The reliability index.|
All data generated or analyzed during this study are included in this article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work was supported by the National Key Fund (51890903) and the National Key R&D Program of China Grant (grant numbers 2016YFC 0701301-01)—Probabilistic model and time-varying reliability of the complex environmental effects of industrial buildings.
The specimen statistics are provided as the supplementary materials. (Supplementary Materials)
K. M. C. Konthesingha, M. G. Stewart, P. Ryan, J. Ginger, and D. Henderson, “Reliability based vulnerability modelling of metal-clad industrial buildings to extreme wind loading for cyclonic regions,” Journal of Wind Engineering and Industrial Aerodynamics, vol. 147, pp. 176–185, 2015.View at: Publisher Site | Google Scholar
Y. K. Wen, “Statistical combination of extreme loads,” Journal of the Structural.Division, vol. 103, no. 5, pp. 1079–1093, 1977.View at: Google Scholar
J. X. Gong and G. F. Zhao, “Combination of relevant load effects and reliability calculation,” Journal of Engineering Mechanics, vol. 18, no. 4, pp. 1–6, 2001.View at: Google Scholar
Q. W. Li and C. Wang, “Effect of load random process relevance on structural time-varying reliability,” Journal of .Tsinghua University: Natural Science Edition, vol. 54, no. 10, pp. 1316–1320, 2014.View at: Google Scholar
J. T. Yao, “Time-interval analysis method for structural reliability,” Journal of Civil Engineering, vol. 10, no. 7, pp. 1–5, 2005.View at: Google Scholar
Z. M. Lin, Reliability Design and Evaluation of Engineering Structures, People’s Communications Publishing House, Beijing, China, 1990.
M. Dupraz, S. I. Rao, and H. Van Swygenhoven, “Large scale 3-dimensional atomistic simulations of screw dislocations interacting with coherent twin boundaries in Al, Cu and Ni under uniaxial and multiaxial loading conditions,” Acta Materialia, vol. 174, pp. 16–28, 2019.View at: Publisher Site | Google Scholar
A. Sklar, Fonctions de repartition a n dimensions et leurs marges, Publications De I’Institut De Statistique De L’Uinversite De Paris, Paris, France, 1959.
U. Cherubini, E. Luciano, and Vecchiato, Copula Methods in Finance, John Wiley & Sons, Chichester, England, 2004.
Z. B. Tu, M. F. Huang, and W. J. Lou, “Correlation of wind loads on building dynamic based on Copula function,” Journal of Zhejiang University-Science, vol. 48, no. 8, pp. 1370–1375, 2014.View at: Google Scholar
S. B. Hu, Z. L. Jiang, and X. R. Zhu, “Research on Parzen window based on improve Gaussian matrix in medical image registration,” Journal of Computational Information System, vol. 8, no. 12, pp. 5103–5110, 2012.View at: Google Scholar
A. Gramacki, Nonparametric Kernel Density Estimation and its Computational Aspects, Springer, New York, NY, USA, 1st edition, 2018.
A. M. Elgammal, Efficient Nonparametric Kernel Density Estimation for Real Time Computer Vision, University of Maryland at College Park, College Park, MD, USA, 2002.