Next Article in Journal
Interfacial Microstructure Analysis of AZ31 Magnesium Alloy during Plastic Deformation Bonding
Next Article in Special Issue
Sliding Mode Observer-Based Fault Detection and Isolation Approach for a Wind Turbine Benchmark
Previous Article in Journal
Thermoeconomic Evaluation and Optimization of Using Different Environmentally Friendly Refrigerant Pairs for a Dual-Evaporator Cascade Refrigeration System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modified Dimension Reduction-Based Polynomial Chaos Expansion for Nonstandard Uncertainty Propagation and Its Application in Reliability Analysis

Department of Chemical & Biomolecular Engineering, Clarkson University, Potsdam, NY 13699, USA
*
Author to whom correspondence should be addressed.
Processes 2021, 9(10), 1856; https://doi.org/10.3390/pr9101856
Submission received: 23 September 2021 / Revised: 11 October 2021 / Accepted: 15 October 2021 / Published: 19 October 2021

Abstract

:
This paper presents an algorithm for efficient uncertainty quantification (UQ) in the presence of many uncertainties that follow a nonstandard distribution (e.g., lognormal). Using the polynomial chaos expansion (PCE), the algorithm builds surrogate models of uncertainty as functions of a standard distribution (e.g., Gaussian variables). The key to build these surrogate models is to calculate PCE coefficients of model outputs, which is computationally challenging, especially when dealing with models defined by complex functions (e.g., nonpolynomial terms) under many uncertainties. To address this issue, an algorithm that integrates the PCE with the generalized dimension reduction method (gDRM) is utilized to convert the high-dimensional integrals, required to calculate the PCE coefficients of model predictions, into several lower-dimensional ones that can be rapidly solved with quadrature rules. The accuracy of the algorithm is validated with four examples in structural reliability analysis and compared to other existing techniques, such as Monte Carlo simulations and the least angle regression-based PCE. Our results show our algorithm provides accurate UQ results and is computationally efficient when dealing with many uncertainties, thus laying the foundation to address UQ in complex control systems.

1. Introduction

Uncertainty, originating from the inherent randomness of a complex system, is common in first-principle models widely used in various engineering problems. Uncertainty quantification (UQ), which quantitatively studies the impact of uncertainty on a system’s performance, is important since uncertainty in model parameters (e.g., external loads, geometry, and material properties in a complex system) can degrade the accuracy of model predictions, thus affecting prediction-based control design and analysis. Many UQ approaches have been developed and applied to engineering problems to improve a system’s performance [1,2,3,4,5,6,7,8]. For example, reliability assessment has been explored for the controlled structures of viscoelastic dampers [3], and the stochastic responses of bridge structures under uncertainty have been studied [4]. Further, reliability-based design optimization (RBDO), used to consider uncertainty in process design, has been developed for structural analysis, including the sequential approximate optimization to deal with multimodal random variables [6] and the adaptive surrogate-based algorithm to replace expensive-to-evaluate models [7]. A hybrid reliability-based optimization method was also proposed for thermal structure design [9]. However, several challenges remain unsolved. For example, a major challenge of existing UQ tools is how to quantify the effect of a large number of uncertainties on model responses in an accurate and computationally efficient manner.
There are several popular UQ methods, including Monte Carlo (MC) simulations [10,11], Taylor expansion-based methods [12,13], and polynomial chaos expansion (PCE) [14,15]. Among these, MC has been widely used since it is easy to implement. It only requires random generation of samples from the distribution of uncertainty and repetitive execution of the first principle models with each sample. Despite the apparent simplicity, MC and its associated approaches are insufficient because they often require a large number of samples to ensure UQ accuracy, thus increasing the computational cost. To reduce the computational cost, different variants of MC were developed, which include an active learning approach to combine Kriging and MC simulation (i.e., AK-MCS) [8]. Further, the Taylor expansion-based method is another way for UQ, which is effective for less complicated models [16,17]. This includes the first-order reliability method (FORM) and second-order reliability method (SORM) [12,18] in structural reliability. However, when models involve complex nonpolynomial terms, the UQ accuracy of the Taylor expansion-based method cannot be guaranteed [12,13,18].
Recently, PCE has been used in different areas including the communities of structural reliability and control theory [4,19]. This approach was proposed by Wiener [14] to quantify uncertainty that follows a normal distribution and later extended to uncertainty with several standard distributions (e.g., uniform) by Xiu [20,21]. The PCE-based methods approximate uncertainty and its effect on model responses with a spectral expansion as a function of PCE coefficients. The PCE coefficients of uncertainty in models are given by a prior probability density function (PDF), but coefficients of model outputs are unknown. Thus, the key of PCE is to solve the PCE coefficients of model responses. Once the coefficients of model responses (e.g., performance function) are obtained, they can be used to approximate the PDF of performance function [19,20] and to evaluate failure probability in reliability analysis [18,22,23]. Since each model response can be represented explicitly with random variables of input uncertainties, it is useful to obtain the statistics of the response. Depending on how the PCE coefficients of model responses are calculated, it is categorized into intrusive and nonintrusive methods [20,24,25].
For the intrusive method, the coefficients of model outputs are obtained by projecting the governing equations into PCE, which requires calculating the inner product between a basis polynomial and the equations described with a truncated PCE in a stochastic Galerkin approach. The inner product generates a family of nested deterministic models to represent the original stochastic models. Calculating the coefficients with the inner product can be computationally prohibitive when models have complex functions (e.g., nonpolynomial functions) and many uncertainties. Nonintrusive PCE is another way, which has been used in the reliability analysis of bridge, truss, and frame structures [4,19,26,27]. Similar to MC, the convergence and accuracy of nonintrusive PCE depend on samples used for UQ. These samples are often referred to as collocation points. While practical, the convergence and accuracy of the nonintrusive PCE-based UQ can be affected by the total number of collocation points and by how the collocation points are generated [28,29]. The nonintrusive methods can be further categorized as a regression-based approach and projection-based approach, depending on how the PCE expressions of output responses are calculated. The former estimates the PCE coefficients by minimizing the mean square error of the approximation of model responses, while each PCE coefficient for the latter is defined with a multidimensional integral. It should be noted that as the number of uncertainties increase, both the size of dimension in the integral and the total number of PCE terms increase. Thus, the projection-based approach can be computationally demanding.
To address these issues, our previous works [29,30] have presented algorithms to combine the PCE with the generalized dimension reduction method (gDRM) or M-variate DRM [31]. This approach approximates high-dimensional integrals in the spectral projection (SP) [28] with a few low-dimensional ones. These low-dimensional integrals involve at most M integration variables [29] that are much smaller than the number of uncertainties. To further reduce the computational cost, the lower-dimensional integrals resulting from the gDRM were calculated with a sampling-based method (i.e., Gaussian quadrature rules) [30]. However, our previous works only deal with uncertainty that follows a standard distribution as done in the classic PCE theory.
Uncertainty may have nonstandard distributions. For example, lognormal distribution describes several natural phenomena, e.g., the bubble size distributions in the gas-liquid and gas-liquid-solid systems [32]. Because there is an explicit relationship between the lognormal and normal variables, we capitalize on the relationship following the work by Ghanem in [33] to efficiently deal with the lognormal uncertainty. Several tools were developed to approximate lognormal random variables [20,24,34,35,36,37,38,39], which include the isoprobabilistic transformation, such as the Rosenblatt or Nataf transformation [37,38], that transforms the original random variable into independent normal variables. Besides, the lognormal distribution was described in [39] by constructing orthogonal polynomials using the Stieltjes–Wigert polynomials. Compared to existing techniques, our algorithm approximates a lognormal distribution with a normal random variable based on the explicit relationship [33] and couples the PCE-based approximation with the gDRM to quantify the effect of many uncertainties on model responses in a computationally efficient way.
In summary, this work expands our previous work to couple the gDRM with the PCE to quantify the effect of nonstandard uncertainties on model outputs. For algorithm validation, lognormal and Weibull distributions are chosen as the testbed, which are approximated with standard random variables [20,33]. In addition, four examples in structural reliability analysis were chosen to show the efficiency of the algorithm and to discuss the UQ accuracy by comparing the results with recent works in the literature.
This paper is organized as follows. A brief description of the PCE-based UQ is given in Section 2. Section 3 shows procedures to approximate nonstandard uncertainty, followed by the gDRM-based PCE to deal with high-dimensional UQ problems. For algorithm verification, four numerical examples in structural reliability analysis are presented in Section 4, and a brief conclusion is given in Section 5.

2. Background of Polynomial Chaos Expansion (PCE)

Suppose that a stochastic process can be defined as u = K ( g ) , where g = ( g 1 ,   ,   g n ) R n is a Gaussian vector involving n uncertain parameters ( n 1 ) and K is the function to describe the relationship between a model response u and g . It is assumed each parameter in g is independent, i.e., any correlation among uncertain parameters is not considered. In this work, uncertainty in model parameters, i.e., input random variable of the models, is hereafter referred to as parametric uncertainty.
To obtain a PCE expression of u , the first step is to rewrite each parametric uncertainty g i as a function of a Gaussian variable x i as in [20]:
g i = g i ( x i ) = k = 0 p g ¯ i , k H k ( x i )
where { g ¯ i , k } are the PCE coefficients to estimate the ith parametric uncertainty, which is defined by the ith standard normal distribution, i.e., x i ~ N ( 0 , 1 ) . These coefficients { g ¯ i , k } are often assumed to be a given a priori or can be estimated with parameter estimation techniques [40]. In (1), H k is the one-dimensional Hermite polynomial basis function, and the polynomial order p defines the number of terms to approximate g i . Uncertainty in model parameters will introduce uncertainty in model responses. Thus, once the PCE coefficients of uncertainty in (1) are available, the PCE expansion of model response u can be subsequently defined with random variables as [20]:
u ( x ) = K ( g ) = j = 0 m 1 u ¯ j H j , n ( x )
where x = ( x 1 , , x i , , x n ) and { H j , n ( x ) } are the jth  n -dimensional orthogonal polynomials defined by basis functions { H k ( x i ) } for each uncertainty. Explicit expressions of { H j , n ( x ) } can be found in [15]. In (2), { u ¯ j } are the deterministic PCE coefficients of model response u , which are unknown and have to be calculated. Further, m is the total number of terms to approximate u ( x ) , which can be calculated with the number of parametric uncertainty n and the polynomial order p as in [20]:
m = ( n + p ) ! p ! n !
Unlike { g ¯ i , k } , PCE coefficients of model response, { u ¯ j } , are unknown and have to be calculated by projecting (2) onto each polynomial basis function H j , n and by using the spectral projection (SP), which can be defined as in [20]:
u ¯ j = u ( x ) H j , n ( x ) H j , n 2 ( x )   ( j { 0 ,   ,   m 1 } )
where the angle brackets · represent the inner product operator for the n -dimensional random space R n defined by x . For example, the inner product of a function f ( x ) is defined as:
f ( x ) = R n f ( x ) W ( x )   d x
where W ( x ) is the weighted function defined by the joint PDFs of random variables x . Note that the inner product of f ( x ) in (5) can be further specified by the expectation operator E [ · ] , i.e., E [ f ( x ) ] = f ( x ) . Accordingly, the PCE coefficient u ¯ j in (4) can be further written as in [20]:
u ¯ j = 1 γ j R n u ( x ) H j , n ( x ) W ( x )   d x
where γ j = E [ H j , n 2 ] is the normalization factor. It is worth mentioning that the jth normalization factor of the Hermite polynomial basis can be easily calculated as γ j = j ! [20]. In addition, the multidimensional integral in (6) is solved over the entire random domain R n of x, and W ( x ) here is the Gaussian weighted function. By substituting the PCE coefficients in (6) into (2), uncertainty in model responses can be rapidly estimated. For example, the statistical central moments μ u and σ u 2 , i.e., mean and variance of model response u , can be analytically calculated as in [20]:
μ u = E [ u ] = E [ j = 0 m 1 u ¯ j H j , n ] = u ¯ 0 E [ H 0 , n ] + E [ j = 1 m 1 u ¯ j H j , n ] = u ¯ 0
σ u 2 = E [ { u E [ u ] } 2 ] = E [ { j = 0 m 1 u ¯ j H j , n u ¯ j , 0 } 2 ] = E [ { j = 1 m 1 u ¯ j H j , n } 2 ] = j = 1 m 1 { u ¯ j } 2 E [ H j , n 2 ]
where H 0 , n 1 is the constant polynomial. As seen in (7) and (8), the first PCE coefficient u ¯ 0 can be used to estimate the mean, while the rest of the coefficients ( u ¯ j 0 ) and the expectation of the squared basis functions ( E [ H j 0 , n 2 ] ) can be used to estimate the uncertainty around the mean of u [20]. Similarly, other higher-order statistical central moments, such as skewness s u and kurtosis κ u , can be rapidly calculated with the high-order PCE coefficients as in [41].
s u = E [ { u E [ u ] } 3 ] 1 σ u 3 = E [ { j = 0 m 1 u ¯ j H j , n u ¯ j , 0 } 3 ] 1 σ u 3 = E [ { j = 1 m 1 u ¯ j H j , n } 3 ] 1 σ u 3
κ u = E [ { u E [ u ] } 4 ] 1 σ u 4 = E [ { j = 0 m 1 u ¯ j H j , n u ¯ j , 0 } 4 ] 1 σ u 4 = E [ { j = 1 m 1 u ¯ j H j , n } 4 ] 1 σ u 4
Details to calculate these statistical central moments can be found in [41,42,43]. These central moments are used to evaluate the reliability index and failure probability for structural reliability analysis in moment methods [18,22]. Thus, they are chosen in this work to discuss the performance of different UQ algorithms in terms of UQ accuracy and computational efficiency in Section 4.

3. Efficient UQ for Nonstandard Uncertainties

3.1. Approximation of Lognormal Uncertainties

Since lognormal uncertainty has an explicit relationship with Gaussian uncertainty, it will be transformed and estimated with the Gaussian variable as in [33]. For clarity, only a parametric uncertainty, which refers to an input random variable of the model, is considered, thus the subscript i in (1) will not be used.
Suppose a normally distributed uncertainty g ( x ) is approximated with a random variable x; and a lognormal uncertainty l ( x ) can be calculated by taking the exponential operator e x p [ · ] of g ( x ) , which results in a relationship as in [33]:
g ( x ) = ln l ( x ) ,   l ( x ) = e x p [ g ( x ) ]
In (11), g ( x ) will be referred to as the normalized Gaussian distribution, and the mean and the variance of g ( x ) are defined as μ g and σ g 2 , respectively. As done for the Gaussian uncertainty in (1), let one assume that l ( x ) can be approximated with a PCE expansion as in:
l ( x ) = k = 0 p l ¯ k H k ( x )
where { l ¯ k } are the PCE coefficients, describing the original lognormal uncertainty, which will be derived based on normalized Gaussian variables. As defined in (1), H k is the Hermite polynomial basis function. Since the first coefficient of PCE in (12) represents the mean value as described in (7), l ¯ 0 can be easily given as μ l = e μ g + ( σ g 2 / 2 ) [33]. For other higher-order terms, we apply the SP to calculate PCE coefficients as a function of the normalized Gaussian random variables. Following the similar procedures to calculate the coefficients of model responses in (4), the resulting PCE expansion of the lognormal uncertainty in (12) can be represented as in [20,33]:
l ( x ) = k = 0 p l ¯ k H k ( x ) = l ¯ 0 k = 0 p σ g k k ! H k ( x ) = l ¯ 0 ( 1 + σ g x + 1 2 ! σ g 2 ( x 2 1 ) + 1 3 ! σ g 3 ( x 3 3 x ) + )
Once the PCE coefficients of a lognormal uncertainty are available, the resulting uncertainty in model response as in (2) can be approximated following the similar procedures as in Section 2.
Compared to the PDF of the Gaussian random variable that has a symmetrical shape, the required polynomial order p in (13) for approximating the original lognormal uncertainty needs to be carefully selected to capture its asymmetric property, such as the long tail. As an example, Figure 1 shows the simulation results of a lognormal uncertainty, where three different polynomial orders (i.e., p = 2, 3, and 4) were used. Note that the PDFs shown in Figure 1 are estimated with the kernel density estimation function in MATLAB. The accuracy to approximate a lognormal distribution with different polynomial orders is compared to MC with 107 samples.
As seen in Figure 1, the accuracy to approximate a lognormal uncertainty can be affected by p and the magnitude of uncertainty. When uncertainty is small, the approximation of a lognormal uncertainty with different p values exhibits almost identical results as in Figure 1a. Whereas, when uncertainty is large, there is a noticeable difference between the PCE-based approximation and MC. However, when p was increased (>2), as seen in Figure 1b, the difference between the PCE-based approximation and MC is insignificant. Thus, we focus on two different values of p (2 vs. 3) in this work to show the performance of the UQ method, which will be discussed in Section 4.

3.2. Approximation of Other Nonstandard Uncertainties

Compared to approximating a lognormal uncertainty with a PCE expression that has an explicit correlation to a Gaussian random variable, the approximation of other nonstandard uncertainties can be made with its probability distribution function [20]. To estimate a nonstandard uncertainty in PCE, suppose Y is an uncertainty with nonstandard distributions (e.g., a Weibull distribution), which has a PCE expression as:
Y ( x ) = k = 0 p Y ¯ k Ψ k ( x )
where { Y ¯ k } are PCE coefficients, and { Ψ k } are polynomial basis functions. Considering Y follows a nonstandard distribution, the calculation of { Y ¯ k } is difficult since the explicit correlation between Y and x is unknown [20]. In this work, a technique in [20] is used to estimate an uncertainty that follows a nonstandard distribution. The resulting PCE coefficients in (14) can be calculated as:
Y ¯ k = 1 γ k E [ F Y 1 ( F x ( x ) ) Ψ k ( x ) ] = 1 γ k R n F Y 1 ( F x ( x ) ) Ψ k ( x ) W ( x ) d x
where F x is the cumulative distribution function (CDF) defined as F x ( x ) = x W ( t ) d t , and F Y 1 is the inverse of a cumulative distribution function of Y , i.e., F Y . For example, let one suppose that a Gaussian random variable x is used to approximate a Weibull-distributed uncertainty Y . Then, W ( x ) is the probability distribution function of x ; Ψ k is the Hermite polynomial basis function; F x is the CDF of the Gaussian random variable x ; and F Y 1 is the inverse of the CDF of Weibull-distributed uncertainty Y . Details about the derivation of (15) and its approximation are not given for brevity and can be found in [20].

3.3. Modified gDRM-Based PCE Using Quadrature Rules

When the PCE coefficients of uncertainty are available, such as (1) for the normal and (12) for lognormal uncertainty, SP is used to calculate the PCE coefficients of model response in (2), which requires calculating multivariate integrals as in (6). As in (3), the total number of terms ( m ) to accurately approximate uncertainty in model response in (2) is a function of the polynomial order ( p ) and the total number of uncertainties ( n ). When p is large (e.g., >2) and when the number of uncertainties increases, the number of terms in (2) can be greatly increased. This can increase the computational cost to calculate the PCE coefficients of model outputs. To reduce the computational cost, the gDRM will be used with quadrature rules to quickly calculate the PCE coefficients as in [30]. This approach is hereafter referred to as the modified gDRM-based PCE (or mgDRM-PCE).
Suppose that a continuous and differentiable function can be defined as f ( x ) = u ( x ) H j , n ( x ) , which is a portion of the integrand in (6). Then, the multivariate integral in (6) is rewritten as:
E [ f ( x ) ] = R n f ( x ) W ( x )   d x
where W ( x ) represents the weighted function described by the joint PDFs of x . Note that the function f ( x ) is defined to easily calculate a specific PCE coefficient u ¯ j in (6), and this expression can also be reused for different j = 0 ,   1 ,   ,   m 1 .
Following the definition of gDRM [31] and the representation in (16), each of the unknown PCE coefficients of the model response in (6) can be calculated as in [29]:
u ¯ j ( x ) = 1 γ j E [ f ( x ) ] 1 γ j r = 0 M ( 1 ) r ( n M + r 1 r ) E [ f M r ]
where f M r in the expectation operator E [ · ] is a function of ( M r ) random variables selected from x . In this way, several lower-dimensional functions, involving up to M random variables of x , can be used to approximate the original function f, which is written as in [31]:
f M r = d 1 < d 2 < < d M r f ( 0 , , 0 , x d 1 , 0 , , 0 , x d 2 , 0 , , 0 , x d M r , 0 )
where d 1 ,   d 2 ,   ,   d M r { 1 ,   2 , , n } . Since the gDRM considers combinations of the random variables in x , the total number of terms in (18) can be calculated as ( n M r ) . Notably, the last term in (18), when M = r , can be calculated as f 0 = f ( 0 ) , which can be solved by setting the mean values of all random variables x to 0. For example, when a bi-variate dimension reduction method (BiDRM) is used, i.e., M = 2, the n -variate function f ( x ) will be approximated with ( n 2 ) two-variate functions, ( n 1 ) one-variate functions, and a constant term f 0 . Thus, the expectation of f M r , i.e., E [ f M r ] in (17), results in a ( M r ) -dimensional integral for each combination of ( M r ) -random variables selected from x, which gives ( n M r ) ( M r ) -dimensional integrals in total.
As noted in [29,30], the accuracy to approximate a high-dimensional integral with several lower-dimensional ones increases as M increases. For example, compared to the BiDRM, a tri-variate DRM (TriDRM) approximates a high-dimensional integral in (16) with ( n 3 ) three-dimensional integrals, ( n 2 ) two-dimensional integrals, ( n 1 ) one-dimensional integrals, and a constant term f 0 . As such, extra effort is required to calculate ( n 3 ) additional three-dimensional integrals. Since the total number of lower-dimensional integrals increases, the computational time to approximate (16) may increase. To address this issue, we will estimate the resulting lower-dimensional integrals with quadrature rules [44] to reduce the computational cost [30,31]. It should be noted that this approach reduces the computational difficulty to evaluate high-dimensional integrals using sampling methods, which appears to be methodologically similar to the anchored analysis-of-variance (ANOVA) expansion in [45]. However, the anchored ANOVA expansion is dependent on the choice of the anchor point, which impacts the accuracy of the expansion and the truncation dimension [45], while the mgDRM does not require any strategy for the decomposition procedure in (18). A detailed description on the ANOVA expansion can be found in [45,46].
To calculate the expectation of one-variate function E [ f 1 ] —an integral only involves a single random variable, e.g., x d 1 of f 1 in (18), a one-dimensional quadrature rule can be defined as [30]:
f ( 0 , , 0 , x d 1 , 0 ,   , 0 ) W ( x d 1 ) d x d 1 q 1 = 1 θ 1 f ( 0 , , 0 , x d 1 q 1 , 0 ,   , 0 ) · α d 1 q 1
where { x d 1 q 1 ,   α d 1 q 1 } q 1 = 1 θ 1 is a set of quadrature points ( x d 1 q 1 ) and their corresponding weights ( α d 1 q 1 ) used to compute the one-dimensional integral resulting from the gDRM. In this work, Gauss–Hermite quadrature rules [44] are used to estimate low-dimensional integrals in (17) in the gDRM step.
To generate quadrature points to approximate multidimensional integrals in (17), a full tensor product grid can be constructed based on the one-dimensional quadrature rules. Note that these lower-dimensional integrals resulting from the gDRM only involve a few integration variables. Thus, the approximation with the full tensor product grid can be finished in real time.
To better illustrate the quadrature points-based approximation, suppose the PCE coefficients of model outputs in (17) are calculated with the TriDRM. This implies that these lower-dimensional integrals to approximate the multidimensional integrals in SP involve at most three random variables in x (or M = 3). That is, a multidimensional integral ( E [ f ( x ) ] ) is estimated with ( n 1 ) one-dimensional, ( n 2 ) two-dimensional, and ( n 3 ) three-dimensional integrals in total, which can be computed with the quadrature rules as in [30]:
E [ f 1 ] d 1 { q 1 = 1 θ 1 f ( 0 , , 0 , x d 1 q 1 , 0 ,   , 0 ) · α d 1 q 1 }
E [ f 2 ]   d 1 < d 2 { q 1 = 1 θ 1 q 2 = 1 θ 2 f ( 0 , , 0 , x d 1 q 1 , 0 , , 0 , x d 2 q 2 , 0 , , 0 ) · ( α d 1 q 1 α d 2 q 2 ) }
E [ f 3 ]   d 1 < d 2 < d 3 { q 1 = 1 θ 1 q 2 = 1 θ 2 q 3 = 1 θ 3 f ( 0 , , 0 , x d 1 q 1 , 0 , , 0 , x d 2 q 2 , 0 , , 0 , x d 3 q 3 , 0 ) · ( α d 1 q 1 α d 2 q 2 α d 3 q 3 ) }
where θ i ( i = 1, 2, 3) is the number of quadrature points,   x d i q i , for each integration variable; α d i q i is the weight of each quadrature point; and means the tensor product. In summary, the total number of quadrature points for low-dimensional integrals can be defined as Q = i = 1 M r θ i = θ i M r [44]. It is important to note that as the dimension of a random space increases, the number of integrals resulting from gDRM increases. In this case, the total number of evaluations of these integrals with quadrature points is also increased, which may increase the computational cost. This issue will be discussed with benchmark examples in Section 4.

3.4. Summary of the UQ Algorithm

Our algorithm quantifies the effect of nonstandard uncertainty on model responses by coupling PCE with the gDRM. A flowchart summarizing the procedures of the modified gDRM-based PCE to deal with nonstandard random variables is shown in Figure 2.
As in Figure 2, nonstandard uncertainties are firstly approximated with standard random variables. Once the PCE-based surrogate model of parametric uncertainty, i.e., an input random variable of the models is available, model responses can be approximated with coupled PCE coefficients as in (2). The PCE coefficients of model responses are calculated with (4). This produces high-dimensional integrals, which are not trivial to solve, especially when the number of uncertainties is large and when the model involves nonpolynomial terms. To address this, the gDRM, such as the BiDRM and TriDRM noted in Section 3.3, is used to convert the high-dimensional integral into few low-dimensional ones. To further decrease the computational cost, quadrature points generated with Gaussian-quadrature rules are used to numerically solve these low-dimensional integrals. The approximation of lower-dimensional integrals is highlighted in green boxes in Figure 2. Once the PCE coefficients of model responses are available, the statistical moments of model outputs (e.g., mean, variance, skewness, and kurtosis) can be quickly calculated with (7)~(10).

4. Benchmark Examples in Structural Reliability Analysis

Four examples in structural reliability analysis are chosen to illustrate the performance of the UQ algorithm to deal with nonstandard uncertainties. The correlation among uncertainties is not considered, since our objective is to study the applicability of the algorithm to tackle many uncertainties. When uncertainties are correlated, the number of random variables to approximate uncertainties will be decreased, thus reducing the number of dimensions. To evaluate the UQ accuracy, statistical moments (mean, standard deviation, skewness, and kurtosis) are calculated and compared with other techniques, which include MC simulations and the least angle regression-based PCE.

4.1. Example 1: Linear Performance Function

As given in Figure 3, a linear performance function is used to study the plastic collapse mechanisms of a one-bay frame, which is mathematically described as in [47,48]:
Z = u ( X ) = X 1 + 2 X 2 + 2 X 3 + X 4 5 X 5 5 X 6
Following [47,48], X 1 ~ X 6 in (23) are assumed to be independent and lognormally distributed; the mean and standard deviation of each variable are listed in Table 1. In this case study, all parameters are considered as parametric uncertainties, resulting in a six-dimensional random space. Since uncertainty follows a lognormal distribution, Hermite polynomials were used as the basis functions to build the PCE-based surrogate model of the response Z in (23).
The first step to build a surrogate model of Z in (23) is to formulate the PCE expressions of lognormally distributed uncertainties, X 1 ~ X 6 . In this case study, (13) was used to build a PCE model for each uncertainty with normalized Gaussian variables that are related to the lognormal uncertainty as in Section 3.1. The second coefficients of the normalized Gaussian variables (i.e., σ g for X 1 ~ X 6 ) were calculated as 0.0998 for X 1 ~ X 4 and 0.2936 for X 5 and X 6 , respectively. These were also used to derive the explicit PCE expressions of lognormal uncertainties as in (13). Once the PCE expressions of uncertainties ( X 1 ~ X 6 ) were available, the mgDRM-PCE in Section 3.3 was used to compute the PCE coefficients of Z in (23).
In this case study, different polynomial orders of X 1 ~ X 6 ( p = 2 vs. 3 ) and different numbers of random variables in the gDRM ( M = 2 vs. 3 ) were used. The UQ accuracy with different combinations of p and M was compared to the results of MC by calculating the relative error ( ϵ R ), which is defined as in:
ϵ R = | e M C e   e M C |
where e M C is a reference of a specific statistical moment calculated with MC, and e represents the corresponding statistical moment calculated with other UQ methods. The simulation results are summarized in Table 2 and compared with existing algorithms: a recent work in [48] that integrates the BiDRM with a high-order unscented transformation [49] and the least angle regression-based PCE (LAR-based PCE) [50,51]. For MC, 107 samples were used for each uncertainty to ensure UQ accuracy.
In Table 2, mBiDRM-PCE means that the quadrature rules were used to estimate integrals from the BiDRM, which involve at most two integration variables, while mTriDRM-PCE means that the quadrature rules were applied to the TriDRM, which has at most three integration variables. The BiDRM converted a six-dimensional integral into 15 two- and 6 one-dimensional ones, while 20 three-, 15 two-, and 6 one-dimensional integrals were used in the TriDRM. In addition, five quadrature points of each dimension were generated following the Gauss–Hermite quadrature rules (i.e., θ i = 5), for which a full-tensor product grid was built to calculate the low-dimensional integrals from the BiDRM and the TriDRM.
In addition, for the LAR-based PCE, the LAR procedure was implemented to select a set of basis polynomials to build a sparse PCE and estimate the unknown PCE coefficients of the response Z in (23) with the least square regression. It is important to note the LAR-based PCE in this work does not involve the adaptive solver as in [52] because our objective is to assess the accuracy of the surrogate model derived by the gDRM-based PCE. To obtain model evaluations of Z , 1000 pairs of random samples of parametric uncertainties were used. Thus, the total number of model runs is 1000. Further, the LAR method was used to identify the sensitive basis polynomials of output Z until the maximum correlation between the residual and the basis polynomials is less than a preselected effective value, which was set to 1.0 × 10−6. In this case study, two different polynomial orders ( p = 2 and 3) were studied. Thus, the number of PCE terms of LAR-based PCE was 13 and 19, respectively. In contrast, the number of PCE terms for our algorithm was 28 and 84, respectively.
As seen in Table 2, the difference in the mean value of Z ( μ Z ) calculated with the mBiDRM-PCE and mTriDRM-PCE is insignificant, and their results are almost identical to MC. For the other higher-order statistical moments, it was found that the relative error is relatively larger when p was set to 2, as compared to the results when p was set to 3. This indicates that more polynomial terms to approximate a lognormal uncertainty can improve the UQ accuracy. Besides, compared to the LAR-based PCE, our algorithm shows identical results for the four statistical moments, which demonstrates the accuracy of the gDRM-based PCE. It was also found that, when p was set to 3, the mBiDRM-PCE outperforms the results in a recent work [48], for which only the BiDRM was used. This shows the advantage of combining the PCE with the gDRM, because PCE can not only provide an analytical expression for model responses as a function of random variables, but also explicitly quantify the impact of approximation on UQ accuracy.
For comparison, the simulation results of mBiDRM-PCE and mTriDRM-PCE with different polynomial orders and MC are shown in Figure 4, where the cumulative density function (CDF) of each UQ method was approximated with the built-in kernel density estimation function in MATLAB. As seen, the CDF of Z in (23) is almost identical to the results obtained with MC, when p was 3. This shows our algorithm can accurately quantify the effect of lognormal uncertainties on model responses.

4.2. Example 2: Roof Structure

The performance function of a roof truss structure as illustrated in Figure 5 was used to study the UQ accuracy of our algorithm to deal with different nonstandard uncertainties, when the system involves nonlinear polynomial functions. The performance function is described as in [19,53]:
Z = u ( X ) = 0.03 0.5 q l 2 ( 3.81 A 2 E 2 + 1.13 A 1 E 1 )
where q and l represent a vertical load and roof span that are used to calculate the nodal load P = q l / 4 , and A and E are the cross-sectional area and elastic modulus, respectively. In this case study, all parameters in the performance function in (25) are uncertain with predefined PDFs. This yields a six-dimensional random space ( n = 6 ) . Details about the statistical description of uncertainties are given in Table 3.
Like Example 1, Hermite polynomials were used as the basis functions, and the PCE expressions of lognormal uncertainties were constructed with normalized Gaussian variables. In addition, it is assumed that l follows a Weibull distribution, for which the PCE expression was approximated with standard random variables using its PDF as described in Section 3.2. Additionally, to compute the unknown PCE coefficients of the performance function Z in (25), the mBiDRM and mTriDRM were used to reduce the computational cost by converting the high-dimensional integrals into low-dimensional ones, which were further approximated with quadrature points constructed by the Gaussian–Hermite quadrature rules. Note that the number of quadrature points for each dimension θ i was set to 5 as done in Example 1. To quantify the UQ accuracy, the relative errors ϵ R defined in (24) were calculated by referring to the results of MC. For MC, 107 samples were used to ensure UQ accuracy.
The results of Example 2 are summarized in Table 4. It is worth mentioning that the results from previous works (e.g., [19,53]) were not given. This is because the uncertainty in [19,53] follows a normal distribution and we intentionally introduced lognormal and Weibull distributions to test the performance of our method to deal with different types of nonstandard uncertainties. Following Example 1 in the previous section, however, the LAR-based PCE was chosen to compare the UQ accuracy with two different polynomial orders ( p = 2 and 3). The simulation results are given in Table 4. The total number of PCE terms for our algorithm is 28 and 84, when the polynomial order for each uncertain parameter was set to 2 and 3, respectively, while the number of PCE terms of the LAR-based PCE is 28 and 77 for each polynomial order, respectively.
Similar to Example 1, it was found that the UQ accuracy can be affected by the polynomial order ( p ) and the number of random variables in the gDRM ( M ). It is important to note that the mean values ( μ Z ) of different PCE-based gDRM methods in Table 4 reached the same result (0.0064), since only two significant figures were shown. Specifically, as the number of p and M increases, the UQ accuracy can be improved. For example, when M was set to 3 for the mTriDRM and when p was set to 3, the relative error ϵ R of κ Z is one order of magnitude smaller than other methods (~0.0202%). It was also found that when the polynomial order p was set to 2, the UQ accuracy of LAR-based PCE is close to our method (mBiDRM and mTriDRM) for all four statistical moments. When p was set to 3, the UQ accuracy of mTriDRM-PCE has the same order of magnitude as the LAR-based PCE. This clearly shows the algorithm in this work can deal with nonlinear problems with nonstandard uncertainties and provide accurate results as one of the most popular techniques in the literature.
For comparison, Figure 6 shows the CDFs of the performance function Z in (25), with different combinations of p and M . As seen, the UQ results of our algorithm converge to MC, as the number of p and M increases. It is important to note that the accuracy of MC is related to the total number of samples used for UQ. When the model is highly nonlinear and when the number of uncertainties is large, a large number of samples is required, thereby increasing the computational cost for MC. This is further discussed with two more complicated examples below.

4.3. Example 3: Truss Structure with 13 Members

A 13-bar truss structure [48] was chosen to study the performance of the UQ algorithm for dealing with a combination of different types of uncertainties. A schematic of the truss structure is shown in Figure 7, which involves 8 nodes and 13 members.
In this example, it is assumed that there are three external loads, i.e., P 1 , P 2 , and P 3 , imposed on nodes 6, 7, and 8, respectively, which follow a normal distribution. Additionally, the elastance E and the sectional area A for each member are assumed to be independent and lognormally distributed. The statistical description of each parameter is listed in Table 5. Finite element analysis was used to quantify uncertainty in the performance function Z of the 13-bar truss structure, which is mathematically defined as in [48]:
Z = u ( X ) = h m a x h
where h m a x is the threshold describing the maximum allowable deflection, i.e., displacement, and h is the vertical displacement on node 3. In reliability analysis, when the performance function Z exceeds a limit, e.g., zero in (26), it would be considered as a failure event. Additionally, the probability of failure events is referred to as the failure probability, which is often used for structural reliability analysis [54]. Thus, it is essential to calculate the statistics of the prediction h in a precise and computationally efficient way. Here we focus on evaluating the displacement h with different UQ methods. The proposed methods in this work were integrated with finite element analysis to assess the precise prediction h, and the performance was validated with other UQ techniques, such as MC, LAR-based PCE, and the method in [48].
As seen in Table 5, there are two lognormal uncertainties and three normal uncertainties. For the lognormal uncertainty, the normalized Gaussian variables and Hermite basis functions were used to build the PCE surrogate models in (13), while for normal uncertainty, the formulation of PCE models is defined in (1). Once these PCE models of uncertainties are available, the resulting surrogate model of h in (26) is described with unknown PCE coefficients that can be solved with the mgDRM involving five quadrature points in each dimension (i.e., θ i = 5). Two different polynomial orders (2 vs. 3) and two different values of M (2 vs. 3) were considered. Table 6 summarizes the simulation results. It is important to note that the total number of samples for each parametric uncertainty in MC was set to 107 to ensure UQ accuracy.
For the implementation of LAR-based PCE, 1000 model runs were used to evaluate the truss structure model as in Figure 7. Further, the sparse PCE terms were determined with the LAR algorithm to find the sensitive basis polynomials with a predefined value (i.e., 1.0 × 10−6) as done in previous examples. Note that the value was used to terminate the LAR algorithm when the maximum correlation between the residual and the basis polynomials was below the value. Details of the LAR procedures can be found in [50,51]. In addition, the resulting number of PCE terms for the LAR-based PCE was 21 and 40, when the polynomial order p was 2 and 3, respectively. In contrast, the total numbers of PCE coefficients of our method were 21 and 56, which were determined using (3).
As seen in Table 6, it was found that UQ accuracy can be affected by the polynomial order ( p ) and the total number of integration variables ( M ) used in the gDRM. For example, when the mTriDRM-PCE was used ( p = 3 and M = 3), the relative error ϵ R of κ h is ~0.1628% in Table 6. This is smaller than the results (~0.4273%) in [48]. As compared to the LAR-based PCE, our algorithm provides comparable results. This shows the potential of the gDRM-based PCE to deal with more complicated problems, especially when estimating the higher-order statistical moments since it was previously considered challenging [55]. To graphically study the UQ accuracy with different combinations of p and M , Figure 8 shows the estimated CDFs of the displacement h . As seen, the mTriDRM-PCE, when p was 3, can provide almost identical results, as compared to MC.
As compared to the first two examples, UQ in this case study requires integrating the finite element analysis with the PCE, which may increase the computational burden for UQ. We further studied the computational time to calculate the displacement h with an office desktop (Core i5-8400 central processing unit (CPU) at 2.80 GHz). Using MC with 107 samples for each uncertainty, it was found that ~31.03 min were required to simulate all nodes in Figure 7. For the mBiDRM-PCE, it took ~21.36 and 48.06 s to simulate all nodes, when the polynomial order ( p ) was set to 2 and 3, respectively. In addition, for the mTriDRM-PCE, it was found that ~26.28 and 67.31 s were required, when p was 2 and 3, respectively. The mTriDRM-PCE took longer than mBiDRM-PCE due to the larger number of integrals that need to be solved. Further, depending on the polynomial order p, the requisite number of terms in PCE to approximate the displacement h differs as shown in (3), which leads to the difference in computational times of each method (e.g., mBiDRM- and mTriDRM-PCE).

4.4. Example 4: Planar Truss Structure with 23 Members

A planar truss structure with 23 members as in Figure 9 is considered in this case study. It is assumed that all parameters cannot be known with certainty, resulting in a nine-dimensional random space. This allows us to validate the UQ accuracy of the proposed method for dealing with many uncertainties. Details of model parameters are given in Table 7.
In this case study, E is the elastic modulus; A 1 and A 2 are the cross-sectional area corresponding to the specific member; and P 1 to P 6 represent external loads imposed on nodes from 8 to 13. As done in Example 3, the finite-element analysis was used to calculate the performance function Z as in [19,56]:
Z = u ( X ) = v m a x v
where v m a x represents a specific threshold, i.e., the maximum deflection, and v is the vertical displacement on node 4 that will be approximated with different UQ methods. Our objective in this example is to accurately approximate uncertainty in v such that the failure probability can be quantified for structural reliability analysis.
As seen in Table 7, three parameters ( E , A 1 , and A 2 ) follow a lognormal distribution, while the rest of the parameters follow a Weibull distribution. Thus, the resulting random space is nine-dimensional, i.e., n = 9. To derive a PCE expression of v and to compute the statistical moments, the finite element analysis was coupled with the proposed methods, i.e., mBiDRM- and mTriDRM-PCE, and two different polynomial orders, i.e., p = 2 and 3, were investigated. Note that in these methods, the BiDRM approximates a nine-dimensional integral with 36 two- and 9 one-dimensional integrals, while the TriDRM requires the evaluations of 84 three-, 36 two-, and 9 one-dimensional ones. These integrals were calculated with the Gauss–Hermite quadrature rule [44], where the number of quadrature points for each random variable was set to 5. In addition, the number of samples for MC was set to 107 in order to compute the relative errors ϵ R in (24) for comparison purposes. The UQ results and relative errors ( ϵ R ) are shown in Table 8. Since the distribution of uncertainties is different as compared to [19,56], the results in these works were not given. However, LAR-based PCE was simulated for two different polynomial orders ( p = 2 and 3) for comparison. The simulation results are summarized in Table 8. For the LAR-based PCE, the total number of PCE terms was 55 and 212, when p was set to 2 and 3, respectively. In constrast, 55 and 220 PCE coefficients were used in our algorithm for the output response, when p was set to 2 and 3, respectively.
Similar to previous examples, our method provides accurate results, as the number of polynomial order ( p ) and the integration variables ( M ) in gDRM increase. For example, the mTriDRM outperforms others when p was 3, since the relative errors are smaller as in Table 8. In addition, it was found that the mTriDRM-PCE can provide comparable results to the LAR-based PCE for all statistical moments, thus confirming the accuracy of the mTriDRM-PCE-based algorithm in this work. For comparison, the CDF of the displacement v on node 4 was approximated with the proposed method and is shown in Figure 10. As compared to MC, it was found that mTriDRM-PCE provided the most accurate result, when p was set to 3. This shows the proposed algorithm can deal with complicated problems involving many uncertainties.
We further studied the computational time of the proposed method and compared the results to MC. Similar to the previous case study, it was found that the computational time for MC is larger, as compared to the gDRM-based method. For example, ~55.54 min were required in this case study, when 107 samples were used for each parametric uncertainty for MC. However, the computational time for the mTriDRM-PCE was found to be ~11.52 min, when p was set to 3, which is about 80% lower than MC. It is also important to note that the gDRM-PCE based UQ method, as compared to sampling-based MC, can provide an analytical expression of the displacement v as a function of uncertainties. This provides mathematical explanations to gain a deep understanding of the problem for improved reliability analysis. The analytical expression, for instance, can be combined with sensitivity analysis techniques to find the most sensitive random variable that can significantly affect the system’s performance.

5. Conclusions

Using the polynomial chaos expansion (PCE), an uncertainty quantification (UQ) algorithm is presented to deal with many uncertainties that follow nonstandard distributions. To build PCE-based surrogate models, standard random variables are used to identify the relationship between nonstandard and standard distributions. To reduce the computational cost for calculating the PCE coefficients of model outputs, a generalized dimension reduction method is used to transform a high-dimensional integral resulting from the spectral projection (SP) into a few lower-dimensional integrals, which can be rapidly solved with quadrature rules in real-time. To show the UQ accuracy of our algorithms, four examples of structural reliability analysis were used. As compared to Monte Carlo simulations and other works in the literature, our results show the superior performance of the algorithms in terms of UQ accuracy and computational time. This shows the potential of the algorithm to tackle UQ in more complicated engineering problems that require consideration of many uncertainties. However, when uncertainty has a nonstandard distribution with a large variance, more PCE coefficients of uncertainty might be needed in the proposed approach to ensure UQ accuracy. In this case, it can lead to intensive computational burden in the presence of many uncertainties. Future work will improve the proposed algorithm to address this issue by identifying orthogonal polynomial functions only for a given uncertainty.

Author Contributions

Conceptualization, Y.D.; methodology, Y.D. and J.S.; validation, J.S.; investigation, J.S.; writing—original draft preparation, J.S.; writing—review and editing, Y.D.; supervision, Y.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the National Science Foundation, Division of Civil, Mechanical and Manufacturing Innovation (CMMI), under the award No. 1727487.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANOVAAnalysis of variance
BiDRMBivariate dimension reduction method
CDFCumulative distribution function
gDRMGeneralized dimension reduction method
gDRM-PCEgDRM-based PCE
LARLeast angle regression
MCMote Carlo
mgDRM-PCEModified gDRM-based PCE
PCEPolynomial chaos expansion
PDFProbability density function
SPSpectral projection
TriDRMTrivariate dimension reduction method
UQUncertainty quantification

References

  1. Coit, D.W.; Jin, T.; Wattanapongsakorn, N. System optimization with component reliability estimation uncertainty: A multi-criteria approach. IEEE Trans. Reliab. 2004, 53, 369–380. [Google Scholar] [CrossRef]
  2. Zaman, K.; Mahadevan, S. Reliability-based design optimization of multidisciplinary system under aleatory and epistemic uncertainty. Struct. Multidiscip. Optim. 2017, 55, 681–699. [Google Scholar] [CrossRef]
  3. Xu, J.; Li, J. Stochastic dynamic response and reliability assessment of controlled structures with fractional derivative model of viscoelastic dampers. Mech. Syst. Signal Process. 2016, 72–73, 865–896. [Google Scholar] [CrossRef]
  4. Ni, P.; Xia, Y.; Li, J.; Hao, H. Using polynomial chaos expansion for uncertainty and sensitivity analysis of bridge structures. Mech. Syst. Signal Process. 2019, 119, 293–311. [Google Scholar] [CrossRef]
  5. Kwon, K.; Ryu, N.; Seo, M.; Kim, S.; Lee, T.H.; Min, S. Efficient uncertainty quantification for integrated performance of complex vehicle system. Mech. Syst. Signal Process. 2020, 139, 106601. [Google Scholar] [CrossRef]
  6. Zhang, Z.; Deng, W.; Jiang, C. Sequential approximate reliability-based design optimization for structures with multimodal random variables. Struct. Multidiscip. Optim. 2020, 62, 511–528. [Google Scholar] [CrossRef]
  7. Dubourg, V.; Sudret, B.; Bourinet, J.-M. Reliability-based design optimization using kriging surrogates and subset simulation. Struct. Multidiscip. Optim. 2011, 44, 673–690. [Google Scholar] [CrossRef] [Green Version]
  8. Echard, B.; Gayton, N.; Lemaire, M. AK-MCS: An active learning reliability method combining kriging and monte Carlo simulation. Struct. Saf. 2011, 33, 145–154. [Google Scholar] [CrossRef]
  9. Wang, C.; Qiu, Z.; Xu, M.; Li, Y. Novel reliability-based optimization method for thermal structure with hybrid random, interval and fuzzy parameters. Appl. Math. Model. 2017, 47, 573–586. [Google Scholar] [CrossRef]
  10. Fishman, G.S. Monte Carlo: Concepts, Algorithms, and Applications; Springer: New York, NY, USA, 1996. [Google Scholar]
  11. Schueller, G.I. Efficient Monte Carlo simulation procedures in structural uncertainty and reliability analysis—Recent advances. Struct. Eng. Mech. 2009, 32, 1–20. [Google Scholar] [CrossRef]
  12. Zhang, Z.; Jiang, C.; Wang, G.G.; Han, X. First and second order approximate reliability analysis methods using evidence theory. Reliab. Eng. Syst. Saf. 2015, 137, 40–49. [Google Scholar] [CrossRef]
  13. Smith, C.L. Uncertainty propagation using taylor series expansion and a spreadsheet. J. Ida. Acad. Sci. 1994, 30, 93–105. [Google Scholar]
  14. Wiener, N. The homogeneous chaos. Am. J. Math. 1938, 60, 897–936. [Google Scholar] [CrossRef]
  15. Ghanem, R.G.; Spanos, P.D. Stochastic Finite Elements: A Spectral Approach; Springer Science and Business Media LLC: New York, NY, USA, 1991. [Google Scholar]
  16. Hong, J.; Shaked, S.; Rosenbaum, R.K.; Jolliet, O. Analytical uncertainty propagation in life cycle inventory and impact assessment: Application to an automobile front panel. Int. J. Life Cycle Assess. 2010, 15, 499–510. [Google Scholar] [CrossRef]
  17. MacLeod, M.; Fraser, A.J.; Mackay, D. Evaluating and expressing the propagation of uncertainty in chemical fate and bioaccumulation models. Environ. Toxicol. Chem. 2002, 21, 700–709. [Google Scholar] [CrossRef] [PubMed]
  18. Zhao, Y.-G.; Ono, T. Moment methods for structural reliability. Struct. Saf. 2001, 23, 47–75. [Google Scholar] [CrossRef]
  19. Xu, J.; Kong, F. A cubature collocation based sparse polynomial chaos expansion for efficient structural reliability analysis. Struct. Saf. 2018, 74, 24–31. [Google Scholar] [CrossRef]
  20. Xiu, D. Numerical Methods for Stochastic Computations: A Spectral Method Approach; Princeton University Press: Princeton, NJ, USA, 2010. [Google Scholar]
  21. Xiu, D.; Karniadakis, G.E. The Wiener–Askey polynomial chaos for stochastic differential equations. SIAM J. Sci. Comput. 2002, 24, 619–644. [Google Scholar] [CrossRef]
  22. Wang, C.; Zhang, H.; Li, Q. Moment-based evaluation of structural reliability. Reliab. Eng. Syst. Saf. 2019, 181, 38–45. [Google Scholar] [CrossRef]
  23. Pulch, R. Polynomial Chaos for the Computation of Failure Probabilities in Periodic Problems. In Scientific Computing in Electrical Engineering SCEE 2008; Roos, J., Costa, L.R.J., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 191–198. [Google Scholar]
  24. Le Maître, O.P.; Knio, O.M. Spectral Methods for Uncertainty Quantification: With Applications to Computational Fluid Dynamics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  25. Najm, H.N. Uncertainty quantification and polynomial chaos techniques in computational fluid dynamics. Annu. Rev. Fluid Mech. 2009, 41, 35–52. [Google Scholar] [CrossRef]
  26. Wan, H.-P.; Ren, W.-X.; Todd, M.D. Arbitrary polynomial chaos expansion method for uncertainty quantification and global sensitivity analysis in structural dynamics. Mech. Syst. Signal Process. 2020, 142, 106732. [Google Scholar] [CrossRef]
  27. Wang, C.; Matthies, H.G.; Xu, M.; Li, Y. Dual interval-and-fuzzy analysis method for temperature prediction with hybrid epistemic uncertainties via polynomial chaos expansion. Comput. Methods Appl. Mech. Eng. 2018, 336, 171–186. [Google Scholar] [CrossRef]
  28. Eldred, M.; Burkardt, J. Comparison of Non-Intrusive Polynomial Chaos and Stochastic Collocation Methods for Uncertainty Quantification. In Proceedings of the 47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 5–8 January 2009. [Google Scholar]
  29. Son, J.; Du, Y. Comparison of intrusive and non-intrusive polynomial chaos expansion-based approaches for high dimensional parametric uncertainty quantification and propagation. Comput. Chem. Eng. 2020, 134, 106685. [Google Scholar] [CrossRef]
  30. Son, J.; Du, D.; Du, Y. Modified polynomial chaos expansion for efficient uncertainty quantification in biological systems. Appl. Mech. 2020, 1, 153–173. [Google Scholar] [CrossRef]
  31. Xu, H.; Rahman, S. A generalized dimension-reduction method for multidimensional integration in stochastic mechanics. Int. J. Numer. Methods Eng. 2004, 61, 1992–2019. [Google Scholar] [CrossRef] [Green Version]
  32. Zhang, L.; Li, T.; Ying, W.; Fang, D. Rising and descending bubble size distributions in gas–liquid and gas–liquid–solid slurry bubble column reactor. Chem. Eng. Res. Des. 2008, 86, 1143–1154. [Google Scholar] [CrossRef]
  33. Ghanem, R. The nonlinear gaussian spectrum of log-normal stochastic processes and variables. J. Appl. Mech. 1999, 66, 964–973. [Google Scholar] [CrossRef] [Green Version]
  34. Kundu, A.; Adhikari, S.; Friswell, M.I. Stochastic finite elements of discretely parameterized random systems on domains with boundary uncertainty. Int. J. Numer. Methods Eng. 2014, 100, 183–221. [Google Scholar] [CrossRef]
  35. Mohan, P.S.; Nair, P.B.; Keane, A.J. Stochastic projection schemes for deterministic linear elliptic partial differential equations on random domains. Int. J. Numer. Methods Eng. 2011, 85, 874–895. [Google Scholar] [CrossRef]
  36. Witteveen, J.A.; Bijl, H. Modeling Arbitrary Uncertainties Using Gram-Schmidt Polynomial Chaos. In Proceedings of the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 9–12 January 2006. [Google Scholar]
  37. Rosenblatt, M. Remarks on a multivariate transformation. Ann. Math. Stat. 1952, 23, 470–472. [Google Scholar] [CrossRef]
  38. Nataf, A. Détermination des distributions de probabilité dont les marges sont données. Comptes Rendus Acad. Sci. 1962, 225, 42–43. [Google Scholar]
  39. Ernst, O.G.; Mugler, A.; Starkloff, H.-J.; Ullmann, E. On the convergence of generalized polynomial chaos expansions. ESAIM Math. Model. Numer. Anal. 2012, 46, 317–339. [Google Scholar] [CrossRef] [Green Version]
  40. Du, Y.; Budman, H.; Duever, T. Parameter estimation for an inverse nonlinear stochastic problem: Reactivity ratio studies in copolymerization. Macromol. Theory Simul. 2017, 26, 1600095. [Google Scholar] [CrossRef]
  41. Lasota, R.; Stocki, R.; Tauzowski, P.; Szolc, T. Polynomial chaos expansion method in estimating probability distribution of rotor-shaft dynamic responses. Bull. Pol. Acad. Sci. Tech. Sci. 2015, 63, 413–422. [Google Scholar] [CrossRef] [Green Version]
  42. Krishnamoorthy, K. Handbook of Statistical Distributions with Applications; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  43. Lomax, R.G.; Hahs-Vaughn, D.L. An Introduction to Statistical Concepts, 3rd ed.; Routledge: New York, NY, USA, 2013. [Google Scholar]
  44. McClarren, R.G. Gauss Quadrature and Multi-dimensional Integrals. In Computational Nuclear Engineering and Radiological Science Using Python; Academic Press: Cambridge, MA, USA, 2018; pp. 287–299. [Google Scholar]
  45. Gao, Z.; Hesthaven, J.S. On ANOVA expansions and strategies for choosing the anchor point. Appl. Math. Comput. 2010, 217, 3274–3285. [Google Scholar] [CrossRef] [Green Version]
  46. Cao, Y.; Chen, Z.; Gunzburger, M. ANOVA expansions and efficient sampling methods for parameter dependent nonlinear PDEs. Int. J. Numer. Anal. Model. 2009, 6, 256–273. [Google Scholar]
  47. Der Kiureghian, A.; Lin, H.-Z.; Hwang, S.-J. Second-order reliability approximations. J. Eng. Mech. 1987, 113, 1208–1225. [Google Scholar] [CrossRef]
  48. Xu, J.; Dang, C. A new bivariate dimension reduction method for efficient structural reliability analysis. Mech. Syst. Signal Process. 2019, 115, 281–300. [Google Scholar] [CrossRef]
  49. Zhang, Y.-G.; Huang, Y.-L.; Wu, Z.-M.; Li, N. A high order unscented Kalman filtering method. Acta Autom. Sin. 2014, 40, 838–848. [Google Scholar]
  50. Chen, X. Welcome to Xiaohui Chen’s Webpage, LARS: Least Angle Regression (LARS). Available online: https://publish.illinois.edu/xiaohuichen/code/lars/ (accessed on 28 March 2020).
  51. Efron, B.; Hastie, T.; Johnstone, I.; Tibshirani, R. Least angle regression. Ann Stat. 2004, 32, 407–499. [Google Scholar] [CrossRef] [Green Version]
  52. Blatman, G.; Sudret, B. Adaptive sparse polynomial chaos expansion based on least angle regression. J. Comput. Phys. 2011, 230, 2345–2367. [Google Scholar] [CrossRef]
  53. Song, S.; Lu, Z.; Qiao, H. Subset simulation for structural reliability sensitivity analysis. Reliab. Eng. Syst. Saf. 2009, 94, 658–665. [Google Scholar] [CrossRef]
  54. Ayyub, B.M.; McCuen, R.H. Probability, Statistics, and Reliability for Engineers and Scientists; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  55. Grigoriu, M. Stochastic Calculus: Applications in Science and Engineering; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  56. Lee, S.H.; Kwak, B.M. Response surface augmented moment method for efficient reliability analysis. Struct. Saf. 2006, 28, 261–272. [Google Scholar] [CrossRef]
Figure 1. Illustration of the effect of polynomial order on the approximation accuracy of a lognormal uncertainty. The PCE-based approximation of a lognormal uncertainty is built with the normalized Gaussian variable with different polynomial orders and compared to the MC-based method. The mean value of l was set to 1 in both graphs, but two different standard deviations were used: 0.1 was used in (a) and 0.3 was used in (b).
Figure 1. Illustration of the effect of polynomial order on the approximation accuracy of a lognormal uncertainty. The PCE-based approximation of a lognormal uncertainty is built with the normalized Gaussian variable with different polynomial orders and compared to the MC-based method. The mean value of l was set to 1 in both graphs, but two different standard deviations were used: 0.1 was used in (a) and 0.3 was used in (b).
Processes 09 01856 g001
Figure 2. Schematic of the modified gDRM-based PCE for UQ under many uncertainties that follow nonstandard distributions.
Figure 2. Schematic of the modified gDRM-based PCE for UQ under many uncertainties that follow nonstandard distributions.
Processes 09 01856 g002
Figure 3. Illustration of the one-bay frame for a plastic collapse mechanism [48]. Reprinted from [48], Copyright (2019), with permission from Elsevier.
Figure 3. Illustration of the one-bay frame for a plastic collapse mechanism [48]. Reprinted from [48], Copyright (2019), with permission from Elsevier.
Processes 09 01856 g003
Figure 4. Comparison of the cumulative density functions of Z calculated with the algorithm in this work and MC in Example 1. (a) mBiDRM-PCE and (b) mTriDRM-PCE.
Figure 4. Comparison of the cumulative density functions of Z calculated with the algorithm in this work and MC in Example 1. (a) mBiDRM-PCE and (b) mTriDRM-PCE.
Processes 09 01856 g004
Figure 5. Illustration of a roof truss structure [19,53]. Reprinted from [19], Copyright (2018), with permission from Elsevier.
Figure 5. Illustration of a roof truss structure [19,53]. Reprinted from [19], Copyright (2018), with permission from Elsevier.
Processes 09 01856 g005
Figure 6. Comparison of the cumulative density functions of Z calculated with the algorithm in this work and with MC in Example 2. (a) mBiDRM-PCE and (b) mTriDRM-PCE.
Figure 6. Comparison of the cumulative density functions of Z calculated with the algorithm in this work and with MC in Example 2. (a) mBiDRM-PCE and (b) mTriDRM-PCE.
Processes 09 01856 g006
Figure 7. Illustration of a 13-bar truss structure [48]. Reprinted from [48], Copyright (2019), with permission from Elsevier.
Figure 7. Illustration of a 13-bar truss structure [48]. Reprinted from [48], Copyright (2019), with permission from Elsevier.
Processes 09 01856 g007
Figure 8. Comparison of the cumulative density function of the displacement h calculated with the algorithm in this work and with MC in Example 3. (a) mBiDRM-PCE and (b) mTriDRM-PCE.
Figure 8. Comparison of the cumulative density function of the displacement h calculated with the algorithm in this work and with MC in Example 3. (a) mBiDRM-PCE and (b) mTriDRM-PCE.
Processes 09 01856 g008
Figure 9. Illustration of a planar truss structure with 23 members [19,56]. Reprinted from [56], Copyright (2006), with permission from Elsevier.
Figure 9. Illustration of a planar truss structure with 23 members [19,56]. Reprinted from [56], Copyright (2006), with permission from Elsevier.
Processes 09 01856 g009
Figure 10. Comparison of the cumulative density functions of the displacement v calculated with the algorithm in this work and with MC in Example 4. (a) mBiDRM-PCE and (b) mTriDRM-PCE.
Figure 10. Comparison of the cumulative density functions of the displacement v calculated with the algorithm in this work and with MC in Example 4. (a) mBiDRM-PCE and (b) mTriDRM-PCE.
Processes 09 01856 g010
Table 1. Details of the uncertain variables in Example 1 [47,48]. Reprinted from [48], Copyright (2019), with permission from Elsevier.
Table 1. Details of the uncertain variables in Example 1 [47,48]. Reprinted from [48], Copyright (2019), with permission from Elsevier.
VariableDistribution TypeMeanStandard Deviation
X 1 Lognormal12012
X 2 Lognormal12012
X 3 Lognormal12012
X 4 Lognormal12012
X 5 Lognormal5015
X 6 Lognormal4012
Table 2. Summary of statistical moments and their relative errors in Example 1.
Table 2. Summary of statistical moments and their relative errors in Example 1.
UQ Methods μ Z   [ ϵ R ] σ Z   [ ϵ R ] s Z   [ ϵ R ] κ Z   [ ϵ R ]
MC269.9226103.2812 0.52763.6135
mBiDRM-PCE ( p = 2)270.0000
[0.0287%]
103.2173
[0.0619%]
0.4840
[8.2758%]
3.3813
[6.4250%]
mBiDRM-PCE ( p = 3)270.0000
[0.0287%]
103.2703
[0.0105%]
0.5261
[0.2911%]
3.5897
[0.6577%]
mTriDRM-PCE ( p = 2)270.0000
[0.0287%]
103.2173
[0.0619%]
0.4840
[8.2758%]
3.3813
[6.4250%]
mTriDRM-PCE ( p = 3)270.0000
[0.0287%]
103.2703
[0.0105%]
0.5261
[0.2911%]
3.5897
[0.6577%]
LAR-based PCE ( p = 2)270.0000
[0.0287%]
103.2173
[0.0619%]
0.4840
[8.2758%]
3.3813
[6.4250%]
LAR-based PCE ( p = 3)270.0000
[0.0287%]
103.2703
[0.0105%]
0.5261
[0.2911%]
3.5897
[0.6577%]
Reference [48]270.0010
[0.0291%]
103.2146
[0.0645%]
0.5165
[2.1124%]
3.6980
[2.3383%]
Table 3. Details of the uncertain variables in Example 2 [19,53]. Reprinted from [19], Copyright (2018), with permission from Elsevier.
Table 3. Details of the uncertain variables in Example 2 [19,53]. Reprinted from [19], Copyright (2018), with permission from Elsevier.
VariableDistribution TypeMeanStandard DeviationUnit
q Lognormal20,0001400N/m
l Weibull120.12m
A 1 Lognormal9.8 × 10−45.89 × 10−5m2
A 2 Lognormal400 × 10−448 × 10−4m2
E 1 Lognormal1 × 10116 × 109N/m2
E 2 Lognormal2 × 10101.2 × 109N/m2
Table 4. Summary of statistical moments and their relative errors in Example 2.
Table 4. Summary of statistical moments and their relative errors in Example 2.
UQ Methods μ Z   [ ϵ R ] σ Z   [ ϵ R ] s Z   [ ϵ R ] κ Z   [ ϵ R ]
MC0.00640.0024 0.30043.1614
mBiDRM-PCE ( p = 2)0.0064
[0.0065%]
0.0024
[0.0134%]
0.2955
[1.6162%]
3.1189
[1.3427%]
mBiDRM-PCE ( p = 3)0.0064
[0.0060%]
0.0024
[0.0030%]
0.3000
[0.1340%]
3.1485
[0.4067%]
mTriDRM-PCE ( p = 2)0.0064
[0.0065%]
0.0024
[0.0156%]
0.2952
[1.7095%]
3.1186
[1.3512%]
mTriDRM-PCE ( p = 3)0.0064
[0.0060%]
0.0024
[0.0011%]
0.3000
[0.1374%]
3.1620
[0.0202%]
LAR-based PCE ( p = 2)0.0064
[0.0106%]
0.0024
[0.0231%]
0.2933
[2.3358%]
3.1172
[1.3961%]
LAR-based PCE ( p = 3)0.0064
[0.0058%]
0.0024
[0.0033%]
0.2995
[0.2836%]
3.1610
[0.0120%]
Table 5. Details of the uncertain variables used in Example 3 [48]. Reprinted from [48], Copyright (2019), with permission from Elsevier.
Table 5. Details of the uncertain variables used in Example 3 [48]. Reprinted from [48], Copyright (2019), with permission from Elsevier.
VariableDistribution TypeMeanStandard DeviationUnit
E Lognormal20620.6Gpa
A Lognormal50050mm2
P 1 Normal203kN
P 2 Normal203kN
P 3 Normal203kN
Table 6. Summary of statistical moments and their relative errors in Example 3.
Table 6. Summary of statistical moments and their relative errors in Example 3.
UQ Methods μ Z   [ ϵ R ] σ Z   [ ϵ R ] s Z   [ ϵ R ] κ Z   [ ϵ R ]
MC10.68661.78080.46613.3907
mBiDRM-PCE ( p = 2)10.6864
[0.0021%]
1.7802
[0.0306%]
0.4507
[3.3024%]
3.2762
[3.3761%]
mBiDRM-PCE ( p = 3)10.6864
[0.0019%]
1.7807
[0.0030%]
0.4590
[1.5131%]
3.3243
[1.9594%]
mTriDRM-PCE ( p = 2)10.6864
[0.0021%]
1.7803
[0.0273%]
0.4518
[3.0748%]
3.2774
[3.3414%]
mTriDRM-PCE ( p = 3)10.6864
[0.0019%]
1.7808
[0.0000%]
0.4655
[0.1377%]
3.3852
[0.1628%]
LAR-based PCE ( p = 2)10.6863
[0.0032%]
1.7802
[0.0302%]
0.4545
[2.4969%]
3.2810
[3.2371%]
LAR-based PCE ( p = 3)10.6865
[0.0016%]
1.7808
[0.0004%]
0.4652
[0.1860%]
3.3841
[0.1954%]
Reference [48]10.6865
[0.0013%]
1.7809
[0.0083%]
0.4685
[0.5153%]
3.4052
[0.4273%]
Table 7. Details of the uncertain variables in Example 4 [19,56]. Reprinted from [56], Copyright (2006), with permission from Elsevier.
Table 7. Details of the uncertain variables in Example 4 [19,56]. Reprinted from [56], Copyright (2006), with permission from Elsevier.
VariableDistribution TypeMeanStandard DeviationUnit
E Lognormal2.10 × 10112.10 × 1010Pa
A 1 Lognormal2.0 × 10−32.0 × 10−4m2
A 2 Lognormal1.0 × 10−31.0 × 10−4m2
P 1 P 6 Weibull5.0 × 1047.5 × 103N
Table 8. Summary of statistical moments and their relative errors in Example 4.
Table 8. Summary of statistical moments and their relative errors in Example 4.
UQ Methods μ Z   [ ϵ R ] σ Z   [ ϵ R ] s Z   [ ϵ R ] κ Z   [ ϵ R ]
MC0.07940.01180.40673.3010
mBiDRM-PCE ( p = 2)0.0794
[0.0002%]
0.0118
[0.0224%]
0.3976
[2.2522%]
3.2234
[2.3512%]
mBiDRM-PCE ( p = 3)0.0794
[0.0000%]
0.0118
[0.0181%]
0.4131
[1.5571%]
3.3429
[1.2683%]
mTriDRM-PCE ( p = 2)0.0794
[0.0002%]
0.0118
[0.0282%]
0.3951
[2.8698%]
3.2205
[2.4406%]
mTriDRM-PCE ( p = 3)0.0794
[0.0000%]
0.0118
[0.0027%]
0.4086
[0.4477%]
3.3039
[0.0863%]
LAR-based PCE ( p = 2)0.0794
[0.0006%]
0.0118
[0.0162%]
0.3961
[2.6238%]
3.2210
[2.4256%]
LAR-based PCE ( p = 3)0.0794
[0.0008%]
0.0118
[0.0080%]
0.4084
[0.4159%]
3.3032
[0.0667%]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Son, J.; Du, Y. Modified Dimension Reduction-Based Polynomial Chaos Expansion for Nonstandard Uncertainty Propagation and Its Application in Reliability Analysis. Processes 2021, 9, 1856. https://doi.org/10.3390/pr9101856

AMA Style

Son J, Du Y. Modified Dimension Reduction-Based Polynomial Chaos Expansion for Nonstandard Uncertainty Propagation and Its Application in Reliability Analysis. Processes. 2021; 9(10):1856. https://doi.org/10.3390/pr9101856

Chicago/Turabian Style

Son, Jeongeun, and Yuncheng Du. 2021. "Modified Dimension Reduction-Based Polynomial Chaos Expansion for Nonstandard Uncertainty Propagation and Its Application in Reliability Analysis" Processes 9, no. 10: 1856. https://doi.org/10.3390/pr9101856

APA Style

Son, J., & Du, Y. (2021). Modified Dimension Reduction-Based Polynomial Chaos Expansion for Nonstandard Uncertainty Propagation and Its Application in Reliability Analysis. Processes, 9(10), 1856. https://doi.org/10.3390/pr9101856

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop