Abstract

Stochastic nonlinear dependencies have been reported extensively between different uncertain parameters or in their time or spatial variance. However, the description of dependency is commonly not provided except a linear correlation. The structural reliability incorporating nonlinear dependencies thus needs to be addressed based on the linear correlations. This paper first demonstrates the capture of nonlinear dependency by fitting various bivariate non-Gaussian copulas to limited data samples of structural material properties. The vine copula model is used to enable a flexible modeling of multiple nonlinear dependencies by mapping the linear correlations into the non-Gaussian copula parameters. A sequential search strategy is applied to achieve the estimate of numerous copula parameters, and a simplified algorithm is further designed for reliability involving stationary stochastic processes. The subset simulation is then adopted to efficiently generate random variables from the corresponding distribution for high reliability evaluation. Two examples including a frame structure with different stochastic material properties and a cantilever beam with spatially variable stochastic modulus are investigated to discuss the possible effects of nonlinear dependency on structural reliability. Since the dependency can be determined qualitatively from limited data, the proposed method provides a feasible way for reliability evaluation with prescriptions on correlated stochastic parameters.

1. Introduction

Reliability analysis of structural systems involves the distribution modeling of a random vector or a stochastic process , which is used to represent the inherent uncertainties in material properties, geometrical dimensions, loading conditions, or their time/spatial variance in engineering practice [1]. The failure probability of a structural system is defined as the integral of the joint probability density function (PDF) of the random vector over the failure domain defined by the performance function :

Equation (1) implies the need of a joint probability distribution for estimating failure probability. If the random variables are mutually independent, the joint probability distribution will be a simple multiplication of its marginal distributions. However, correlations may exist between different stochastic parameters or in the variance of a stochastic parameter over time or space. For example, correlations are reported between wind speed and wind direction or between Young’s modulus and Poisson’s ratio [2, 3]. The time-dependent or spatially variable stochastic parameter (e.g., time series of external loads, material modulus varying along space, or variation of pore structure in radial distribution [4, 5]) that is often characterized by a stochastic process can also be viewed as correlated random variables with high dimensions [69]. Ignoring these correlations will lead to a biased estimate of the structure reliability.

In structural reliability engineering practice, the probabilistic description of uncertain parameters in terms of univariate distributions and pairwise correlations is ubiquitous. The dependency that shapes the scatter of data points is commonly assumed linear for convenience (e.g., [10, 11]). However, the validity of this assumption is questioned and the corresponding model approximation is coarse if the dependency is nonlinear [3]. In other words, reliability evaluation involving nonlinear data dependencies by linear transformations (e.g., the orthogonal and Nataf transformation) can be inaccurate [12].

The copula theory provides a feasible way of characterizing nonlinear dependencies [13]. It has been widely applied to reliability problems where the data scatter exhibits nonlinear characteristics. For instance, the torsional moment and external forces can be upper tail dependent that is better captured by a Gumbel copula [14]. The time-variant peaks of snow load can also show the Gumbel copula characteristics [15]. For a given copula, the distribution model is usually established by matching it to the correlations of uncertain parameters [16, 17]. There are three major measures of correlations, namely, Pearson’s rho, Spearman’s rho, and Kendall’s tau [18]. The first measure is known as the linear correlation while the other two are the rank correlations. Retrieving copula parameters from the rank correlations is relatively straightforward because the estimate is independent of the marginal distributions [18]. However, the use of linear correlation coefficient is a more common practice in the literature [15]. Compared to the rank correlations, the retrieval of copula parameters from the linear correlations is more complicated because the estimate depends on both the copula function and the marginal distributions [18]. Although some efficient quadrature techniques can be adopted to solve the mathematics [16, 17], the proposed approach to the retrieval of non-Gaussian copula parameters from linear correlation coefficients are restricted to the correlated bivariate cases [16, 17, 19]. The difficulty of generalization to the correlated multivariate cases is rooted in the inflexibility of multivariate nonelliptical copulas for representing multiple nonlinear dependencies [20].

The multiple nonlinear dependencies thus need to be addressed, particularly when the linear correlations of stochastic parameters are ubiquitous in the literature. In this paper, a vine copula approach is adopted for multiple dependence modeling. The copula parameters are retrieved from the univariate distributions and correlation matrix via a sequential search strategy (S3) [21, 22]. In [21], a simplified version of S3 was particularly designed for parameter retrieval in reliability involving stationary stochastic processes (i.e., spatial-variant reliability). Once the joint probability distribution is represented in terms of vine copulas, the iso-probabilistic transformation for nonlinear dependent random variable can be integrated into the reliability approach. Since the expected reliability for structures are often high, the subset simulation (SS) is adopted to estimate the failure probability. The main objectives are to demonstrate the possible nonlinear dependencies in practice and assess its impact on high structural reliability. To this end, Section 2 briefly reviews the copula theory and shows the possibly nonlinear dependencies underlying the stochastic structural material properties. The vine copula approach with the copula parameters retrieved from univariate distributions and correlation matrix by S3 is provided in Section 3. The SS is then extended to nonlinear dependences where the correlated Gaussian/non-Gaussian is targeted, as shown in Section 4. Section 5 uses two examples of finite element (FE) analysis—a frame structure with different stochastic material properties and a cantilever beam with spatially variable stochastic modulus—to discuss the possible effects of nonlinear dependence on structural reliability. Finally, conclusions are presented in Section 6.

2. Nonlinear Dependencies in Structural Material Properties

A copula is an n-dimensional cumulative distribution functions (CDF) such that each of its marginal CDFs is uniform on the interval [13]. Copulas can model the dependence between random variables with arbitrary marginal distributions and link multivariate distributions functions of random variables to their univariate marginal distributions. Let , denote the marginal CDFs of X, and the joint probability distribution of X can be expressed aswhere denotes a dependent, uniformly distributed vector of random variables and , denotes the set of parameters for the copula, and is the copula function, which can be derived aswhere is the inverse function of the corresponding marginal distributions.

In structural systems, the linear dependence characterized by a Gaussian copula is usually assumed for correlated random variables (e.g., the Nataf transformation). In this section, the validity of the Gaussian assumption is examined through experimental data from standard tensile tests of low carbon steels. The low carbon steels are widely used materials in structures and they have various types. Here, a batch of low carbon steels (grade Q235 in Chinese Standard GB/T 17107-1997, equivalently ISO Fe360 [23]) collected from the same construction site is used for experiments. It is known that there exists dependencies between diameter D and yield load YL and between yield load YL and extreme load EL. The best-fit copula to the measured dependence structure can be identified based on the minimum value of the Akaike information criterion (AIC) and the Bayesian information criterion (BIC) [16], which are calculated, respectively, aswhere N is the sample size, k is the number of copula parameters, is the density function of a bivariate copula, and , , and are the coordinates in the unit square space (Figures 1(b), 1(d), and 1(f), respectively) transformed from the measured data in the original space (Figures 1(a), 1(c), and 1(e)), which is defined aswhere denotes the rank of Di, YLi and ELi among the set , , and , respectively. For illustrative purposes, the Gaussian copula, Clayton copula, CClayton copula, and Frank copula, as listed in Table 1, are examined in this paper. The AIC and BIC values for the four copulas are presented in Table 2.

Although the materials are collected from the same batch, the uncertainty of material properties as a result of manufacturing error is observed. Results from Table 2 indicate that the Clayton copula is the best-fit copula for modeling the dependence between D and YL and between D and EL, while the Frank copula is the best-fit copula for modeling the dependence between YL and EL. In other words, none of the dependencies are linear. Thus, the incorporation of nonlinear dependences into structural reliability is necessary. However, the above copulas are fitted only for the bidimensional marginal distributions, while the reliability evaluation requires the joint probability distribution (i.e., the trivariate distribution of D, YL, and EL). The construction of the multivariate distribution that involves the above bidimensional marginal distributions requires a hierarchical nesting of the bivariate copulas, which will be presented in Section 3.

3. Pair-Copula Decomposition Approach to Multivariate Distribution Construction from Marginals and Linear Correlations

3.1. Multivariate Probability Distribution Construction by Pair-Copula Decomposition

A multivariate copula is needed for modeling the multivariate probability distribution. However, the construction of an n-dimensional copula while maintaining high flexibility is challenging due to the compatibility problem [13]. A promising way to solve this problem is the pair-copula decomposition approach [24].

According to Aas et al. [24], the n-dimensional joint PDF of X can be decomposed aswhere the conditional PDF in the form of can be represented aswhere is an arbitrarily selected element of while denotes the vector excluding and is the density function of . Using equation (7) recursively, a multivariate joint PDF can be represented by a product of pair-copula density functions, acting on several different conditional probability distributions. The pair-copula construction involves marginal conditional distributions in the form of , which can also be recursively expressed in terms of the partial derivative of a bivariate copula with respect to the second parameter [25]:

If is a univariate, equation (8) can be simplified as follows:

For high-dimensional distributions, there are generally many different representations of using the pair-copula constructions. The vines are commonly used to organize the decompositions [26]. The C-vine and D-vine [27] are two special cases of regular vines, as shown in Figure 2.

For an n-dimensional problem, both C-vine and D-vine consist of n − 1 trees . Tree Tj has nodes and edges. Each edge corresponds to a pair-copula density function (e.g., edge corresponds to the copula function ) and the edge labels the nodes to be drawn. The edges in tree are determined only by the nodes in tree . The way of connecting nodes in C-vines and D-vines is different. In a D-vine, each node is connected by no more than two edges. In a C-vine, each tree has a unique node that is connected by the other edges.

corresponding to a D-vine and C-vine is, respectively, given as

3.2. Retrieving Vine Copula Parameters from Incomplete Probability Information

The determination of the pair-copula parameters is critical to the modeling of multiple dependences. Given the marginal distributions and correlations, the copula parameter can be retrieved using the following relationship [28]:where and are the means and the standard deviations of , respectively, and is the bivariate joint PDF, which can be obtained from the multivariate joint PDF bywhere denotes the random vector excluding the ith and jth variables. By analogy with equation (6), the conditional PDF can also be written as a product of marginal conditional PDFs. For example,

In particular, for a bivariate distribution , there is only one unconditional bivariate copula function involved, and its parameter can be solely determined by

If is a multivariate distribution, conditional copulas must be involved in the decomposition. However, the associated conditional copula parameters cannot be analytically solved from one single integral function. For example, suppose a trivariate distribution , which can be represented by three copulas: , , and , and three marginal PDFs. To relate the conditional copula parameter with a correlation coefficient (such as ), the conditional distribution has to be represented by other conditional copulas (e.g., ); otherwise, the desired parameter to be estimated will be removed from equation (11) after representing the bidimensional marginal PDF by equation (12). In other words, the analytical approach to retrieving the desired conditional copula parameters requires a simultaneous estimate of irrelevant copula parameters (such as in the above trivariate example). To address the problem, a sequential search strategy (S3) and a simple simulation-based bisection search method proposed in [21, 22] were adopted to approximate the pair-copula parameters. The major advantage of the algorithm is that only one copula parameter is estimated per iteration and the estimate of irrelevant parameters is avoided. For stationary stochastic processes, the correlation coefficient is dependent on the relative separation such that the correlation coefficients between any two neighboring nodes in the D-vine structure are identical. Moreover, the marginal distribution of each node is identical if the stochastic process is stationary. Thus, the pair-copula parameters at each tree level are identical, and the efforts for parameter searching can be significantly reduced. For this reason, the D-vine structured model is particularly suitable for reliability problems involving stationary stochastic processes.

4. Subset Simulation considering Dependent Random Variables

4.1. Subset Simulation

The Monte Carlo simulation (MCS) is widely applied to estimate the failure probability given in equation (1) due to the simplicity and robustness of the algorithm [29]. However, the crude MCS can be unacceptable for its high computational demands when assessing small failure probabilities [30]. To address the limitation, various advanced simulation techniques have been proposed [31]. The SS uses a Markov Chain Monte Carlo (MCMC) sampling strategy, which is identified as one of the most efficient method for estimating small failure probabilities [3234]. Let F denote the failure event in the U-space. The basic idea of SS is to progressively approach the target failure region through a sequence of nested subsets of the random space:where is the certain event. By conditioning on the intermediate failure regions , the failure probability is estimated as a product of conditional probabilities:

The probability is equivalently computed by crude MCS. To estimate the conditional probabilities, samples from the conditional PDFs should be generated, which is

In equation (17), is the indicator function of . Generation of independent and identically distributed (i.i.d.) samples from can be achieved by application of the acceptance-rejection method.

Without loss of generality, suppose that the intermediate and target failure regions are defined as, where . The values of can be chosen adaptively so that the estimates of the conditional probabilities , equal to a constant value of P0.

For details of the implementation procedure of SS, readers are referred to [32].

4.2. Rosenblatt Transformation

The original SS assumes that the random variables are mutually independent (or have been transformed into independent ones). To incorporate the correlations into reliability evaluation, the iso-probabilistic transformation that maps the dependent random variables into independent ones is needed [29]. The probability of failure can thus be expressed in the transformed space aswhere , is the standard normal PDF, and is the limit-state function in the independent standard normal space (U-space), and is the indicator function, where if and , otherwise.

The Rosenblatt transformation can be used for simulating nonlinear dependent random variables [35]. The Rosenblatt transformation is defined aswhere is the inverse function of the CDF of a univariate standard Gaussian distribution and is the marginal conditional distributions. By using non-Gaussian copulas, the nonlinearity of dependency is expressed in the copula representation of the marginal conditional distributions . For simplicity, the excluding variable in the C-vine model is always chosen so thatand of the D-vine models is chosen so that

Then, the dependent random variables in the original space can be obtained from independent random variables using the inverse of the Rosenblatt transformation. In particular, for C-vine models,and for D-vine models,where is the inverse function of equation (8).

5. Numerical Examples

In this section, two examples of FE analysis are used to demonstrate the effect of nonlinear dependency on structural reliability. For both examples, the sample size at each SS simulation level is set to 10000 and the conditional failure probability is set to 0.1. The aim of using a relatively larger sample size is to provide a more robust estimate. A uniform proposal PDF is used for MCMC simulation. The spread of the uniform PDF is taken as twice of the standard deviation of each random variable at the current level.

5.1. Example 1: Three-Bay Five-Story Frame Structure

The three-bay five-story frame structure has been analyzed in many studies [3638]. The schematic diagram of the structure is shown in Figure 3. There are 21 random variables representing various loads and material properties. The 21 random variables are independent by three blocks. The first block involves 3 correlated random variables representing external loads (P1, P2, and P3 with pairwise correlation coefficients of 0.5). The second block involves 2 correlated random variables representing the modulus of beams and columns, respectively (E4 and E5 with a correlation coefficient of 0.9). The third block involves the other 16 correlated random variables representing material properties including cross-sectional area A and moment of inertia I. The cross-sectional area and the moment of inertia of the same element are highly correlated with a correlation coefficient of 0.95 while those of different elements are assumed a correlation coefficient of 0.13 (i.e., and ). The element properties and their statistics are listed in Tables 3 and 4, respectively.

The structure is considered failure if the horizontal displacement of the upper floor exceeds a threshold value of 0.061 m. Thus, the limit state function is

The bivariate Gaussian copula is first adopted for all pair copulas. The result corresponding to this linear dependency should be consistent with that based on the Nataf transformation [12]. Table 5 compares the failure probability by the extended SS with results in previous works. The consistent results indicate the accuracy of the proposed method under linear dependency assumption.

To investigate the impact of nonlinear dependency on the structural reliability, the correlated loads are assumed to have an upper tail correlation that can be captured by a Clayton copula, while the correlated material properties are assumed to have a lower tail correlation that is better modeled by a Clayton copula. The failure probability under this hybrid copulas is , which is 51% higher than the Gaussian copula case. This deviation may seem not significant compared to that caused by other uncertainties. However, the difference can enlarge as the expected failure probability decreases. As shown in Figure 4, denotes the ratio of the failure probability under the nonlinear dependency over that of the linear case, which increased as the tolerant displacement becomes larger. Note that structural systems are often designed towards high reliability, which implies that the effect of dependency can be significant.

To illustrate how reliability results are biased by nonlinear dependency, Figure 5 presents the MCMC simulation in successively generating one pair of correlated random variables (E4 and E5) under different copulas. With the use of a fixed seed, the random walks under different dependencies are exactly the same in the U-space. However, the values differed after mapping into the original space despite a similar route is observed. This effect is also true for multiple nonlinear dependencies. The reason is that the nonlinear characteristics for each pair of uncertain parameters can be maintained using the same non-Gaussian copula in a vine structure regardless of their conditioning order [39]. For example, compared to the Gaussian copula, the cross-sectional area and moment of inertia are more likely to have similar values if any one of these parameters has a small value under the Clayton copula.

5.2. Example 2: The Cantilever Beam with Spatially Variable Modulus

The cantilever beam has been analyzed by several researchers [40, 41]. Consider a cantilever beam with a length of L = 1 unit and a constant cross-sectional area A = 1 unit, as shown in Figure 6. The beam is subjected to a linearly increasing uniaxial load along the x-direction. The modulus of elasticity varies spatially along the x-direction thatwhere unit is the mean of elastic modulus and is a homogeneous Gaussian random field with mean zero and autocovariance function given bywhere x and are the coordinates of two sections of the beam, is the expectation operator, is the standard deviation of , and is the correlation length. Although the proposed method is illustrated with a Gaussian process, the application to non-Gaussian processes is straightforward since the choice of marginal distribution is independent.

Two different meshes were used for the discretization of the FE model and the random field, respectively. To provide a consistent comparison with previous works, the random field is discretized into 256 random variables and further averaged every 16 points. The resulting 16 values of material modulus are then mapped onto the 480 FE meshes. The FE mesh grid is sufficiently fine as the deterministic result converges. The limit state function is defined with respect to the axial displacement at the free end (point B):

Table 6 compares the result with that from previous works, where the random field was parameterized using the spectral decomposition method (i.e., ) is represented by a set of independent standard normal random variables as , where are the eigenvalue and eigenfunction of the covariance kernel satisfying ) [42]. The consistency in failure probability indicates that the proposed method is accurate. The slight deviation can be attributed to the variation of the SS estimator. It has been acknowledged that the spectral decomposition is exact if the random variables within a field are jointly Gaussian [43]. From the copula perspective, jointly Gaussian means Gaussian marginals plus Gaussian copulas [13], which implies a linear dependency is used. However, nonlinear dependencies of stochastic processes have been reported in the literature [15]. Therefore, the effect of nonlinear dependency in a random field is of interest.

To illustrate the impact of nonlinear dependency in a random field, the bivariate Clayton copula and the bivariate Frank copula are, respectively, taken as the pair copulas. Moreover, the effect at different correlation lengths is examined. Figure 7 shows the comparison of failure probabilities as a function of the correlation length. When b approaches to infinity or tends to be zero, the deviation is negligible. This result can be expected as the random variables which are equivalent to a single random variable when fully correlated (), while the zero correlations () lead to the degeneration of all pair copulas to independence copula; thus, the joint probability distribution of the random field is simply a multiplication of its marginals. However, the random field is seldom highly ragged or smooth. For example, the ratio of correlation length over field length is in the order of 10−1 for various properties of CFRP composite panels [44]. Refer to Figure 7; the failure probability from the linear dependency case can be increased or reduced by more than thousands of times if the dependency is nonlinear. Figure 8 presents the CDFs corresponding to different copulas when and , respectively. For highly smooth random field (), the deviation is relatively insignificant regardless of the expected failure probability; however, the deviation becomes more significant as the cantilever beam approaches to high reliability when the spatial variation is moderate (). Figure 9 shows the realizations of the stochastic process under different dependencies. Note the seed for simulation is identical. The value of modulus at each point is different despite the similar pattern of variance.

Figures 10(a)10(d) show the means of modulus at different sections as the simulation level increases from 1 to 4 when different copulas are used. The mean values are generally the same at level 1 (crude MCS). The coordinate corresponding to the minimum mean value of modulus, which indicates the probable failure section, shifted to the left side far away from the free end as the simulation level increased. However, the dependency can influence the minimum mean value at each simulation level. For example, the Clayton copula, which has a lower tail dependency of modulus, makes the targeted failure region reached fast as the minimum mean value of conditional samples decreases sharply. Figures 1112 further show the PDFs of conditional samples at significant and insignificant section, respectively. The great impact of dependency on distribution of conditional samples at significant section is obvious.

6. Concluding Remarks

The investigation to the uncertainties of low carbon steel properties demonstrates the possible nonlinear dependencies in practice. However, the incorporation of this nonlinearity into reliability analysis under the ubiquitous linear correlations is not addressed. This paper extends the subset simulation to structural reliability involving nonlinear dependencies and investigates its effect. The vine copula model is adopted to allow for modeling of multiple nonlinear dependencies, while the model parameters are estimated from linear correlations by a sequential search strategy. The algorithm for model parameter searching is optimized so that high-dimensional reliability problems involving time or spatial variance can be handled efficiently.

The deviation in failure probability due to different dependencies is generally trivial if the expected failure probability tends to be large. However, structural systems are usually designed towards high reliability, where the deviation can increase significantly. For spatially variable parameters, the deviation is more sensitive to the correlation length except for highly ragged or smooth random fields. The deviation in failure probabilities corresponding to different dependencies can reach to a thousand of times if the spatial variation is neither drastic nor smooth. Unfortunately, the correlation length in practice is commonly somewhere between the two extreme cases, which highlighted the importance of dependence modeling.

The accurate modeling of multivariate distributions requires statistical fitting from substantial experimental data, which is not always fulfilled in practice. Instead, the statistical moments such as means, variance, and covariance can be easily obtained from available information. Therefore, a feasible way is to determine the dependency qualitatively from the limited data. Then, the reliability can be quantitatively assessed by integrating the best-fit copulas with given statistics of each uncertain parameter.

Notation

:Cross-sectional area
:Correlation length
:Copula function
:Copula density function
:Diameter
:Modulus of elasticity
:Expectation operator
:Marginal PDF
:Marginal CDF
:Joint PDF
:Failure event
:Performance function
:Moment of inertia
:Failure probability
:Load
:Random variables uniformly distributed on [0, 1]
:Random variables in the U-space
:Random vector in the original space
:Standard Gaussian PDF
:Standard Gaussian CDF
:Parameter of bivariate copula function
:Linear correlation coefficient
:Mean
:Standard deviation
:Autocovariance function
:Displacement threshold value.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (NSFC), under Grant no. 51608399, and “Jiangsu Smart Factory Engineering Research Center (Huaiyin Institute of Technology),” under Grant no. JSFE1903. The authors gratefully acknowledge the assistance of Dr. Fan Wang for his support in the calculation.