Next Article in Journal
Min–Max Dynamic Programming Control for Systems with Uncertain Mathematical Models via Differential Neural Network Bellman’s Function Approximation
Next Article in Special Issue
Uncertainty Analysis and Optimization Modeling with Application to Supply Chain Management: A Systematic Review
Previous Article in Journal
Application of Advanced Optimized Soft Computing Models for Atmospheric Variable Forecasting
Previous Article in Special Issue
Sufficient Conditions for the Existence and Uniqueness of Minimizers for Variational Problems under Uncertainty
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solution-Space-Reduction-Based Evidence Theory Method for Stiffness Evaluation of Air Springs with Epistemic Uncertainty

Key Laboratory of Traffic Safety on Track, Ministry of Education, School of Traffic & Transportation Engineering, Central South University, Changsha 410075, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(5), 1214; https://doi.org/10.3390/math11051214
Submission received: 13 January 2023 / Revised: 15 February 2023 / Accepted: 24 February 2023 / Published: 1 March 2023
(This article belongs to the Special Issue Uncertainty Analysis, Decision Making and Optimization)

Abstract

:
In the Dempster–Shafer evidence theory framework, extremum analysis, which should be repeatedly executed for uncertainty quantification (UQ), produces a heavy computational burden, particularly for a high-dimensional uncertain system with multiple joint focal elements. Although the polynomial surrogate can be used to reduce computational expenses, the size of the solution space hampers the efficiency of extremum analysis. To address this, a solution-space-reduction-based evidence theory method (SSR-ETM) is proposed in this paper. The SSR-ETM invests minimal additional time for potentially high-efficiency returns in dealing with epistemic uncertainty. In the SSR-ETM, monotonicity analysis of the polynomial surrogate over the range of evidence variables is first performed. Thereafter, the solution space can be narrowed to a smaller size to accelerate extremum analysis if the surrogate model is at least monotonic in one dimension. Four simple functions and an air spring system with epistemic uncertainty demonstrated the efficacy of the SSR-ETM, indicating an apparent superiority over the conventional method.

1. Introduction

Air springs are core components of the secondary suspension system of high-speed trains, whose stiffness performance is closely related to the stability and safety of train operations. The stiffness evaluation of air springs is conventionally based on deterministic approaches [1,2]. However, air springs involve inevitable uncertainties, subject to various raw materials and complex processing technologies. Their stiffness performance is more elusive under the coupling effect of the atmosphere and chronic external excitation. Given the crucial role of air springs in train operation, it is imperative to predict their stiffness performance under uncertainty.
There are two categories of uncertainty: aleatory uncertainty and epistemic uncertainty [3]. Aleatory uncertainty is attributed to the inherent variability of the system or working environment [4]. Although probability-based aleatory uncertainty is broadly researched in engineering [5,6,7,8], it is inappropriate for limited data. Epistemic uncertainty results from an incomplete or inadequate knowledge of the physical system [9]. Thus far, non-probabilistic methods such as the Dempster–Shafer evidence theory [10,11], fuzzy set theory [12,13,14], interval analysis [15,16,17,18], and generalised p-boxes [19,20] have been widely applied to describe epistemic uncertainty. Among the non-probabilistic methods mentioned above, the Dempster–Shafer evidence theory supports a relatively flexible mathematical framework equivalent to interval analysis or probability-based methods under particular conditions [9]. Such characteristics promote the extensive application of the Dempster–Shafer evidence theory in reliability analysis and response analysis [21,22,23].
Recently, polynomial chaos expansion (PCE) [24,25,26,27], which approximates computational models through a linear combination of orthogonal polynomials, has been introduced to the Dempster–Shafer evidence theory framework for UQ. Harsheel et al. [28] introduced the Legendre polynomial surrogate to address engineering problems with hybrid uncertainties. Subsequently, Yin et al. [9,29] employed a Jacobi polynomial surrogate and an arbitrary PCE for the response analysis of acoustic systems with epistemic uncertainty. However, the required number of model evaluations for a PCE increases significantly with the size of the uncertain vector; this is also known as the ‘curse of dimensionality’. Consequently, sparse PCE [30,31,32,33], which uses the sparsity-of-effects principle and powerful sparse regression solvers to approximate high-dimensional computational models through a few model evaluations, has been used [34].
Although sparse PCE provides an inexpensive approximation to the computational model, the efficiency of extremum analysis is limited to the size of the solution space. In particular, extremum analysis must be executed in each joint focal element for an evidence-theory-based UQ. Many engineering problems involve a large number of uncertain parameters, and the interval of the uncertain parameters is usually divided into as many focal elements as possible to improve the reliability of UQ. The calculation time of the conventional method [9,29] carrying out extremum analysis in the original solution space is unacceptable. To narrow the size of the solution space, SSR-ETM, which takes advantage of the monotonicity of the PCE-based surrogate model, is proposed. The monotonicity over the range of evidence variables is available for all joint focal elements. If the surrogate model is at least monotonic in one dimension, the size of the solution spaces can be reduced to increase the computational efficiency.
The contributions of this paper are as follows: (1) the Dempster–Shafer evidence theory was employed to quantify epistemic uncertainty; (2) the sparse PCE was introduced to construct the surrogate model, and the input variables with arbitrary probability measures were considered; (3) a method, namely SSR-ETM, was proposed to accelerate evidence-theory-based UQ; (4) different types of numerical examples are designed to discuss the shortcomings and advantages of SSR-ETM, and good results are obtained in an air spring system with epistemic uncertainty.
The remainder of this paper is organised as follows. Section 2 presents the theoretical framework of the Dempster–Shafer evidence theory. Section 3 introduces the SSR-ETM framework and its core steps. In Section 4, four simple functions and an air spring system with epistemic uncertainty are used to verify the superiority of the SSR-ETM. Section 5 presents the concluding remarks.

2. Fundamentals of the Dempster–Shafer Evidence Theory

The Dempster–Shafer evidence theory, abbreviated as evidence theory, owes its name to the pioneering work by Dempster [35] and Shafer [36,37,38]; it is a generalisation of the Bayesian theory of subjective probability [39]. In evidence theory, the frame of discernment (FD) Θ is defined as a finite exhaustive set of w pairwise mutually exclusive propositions as follows:
Θ = { h 1 , h 2 , , h w } .
where h 1 , h 2 , , h w are the elementary propositions.
The combination of all subsets of Θ constitutes a power set 2 Θ that contains the following possible propositions:
2 Θ = , { h 1 } , { h 2 } , , { h w } , { h 1 , h 2 } , { h 1 , h 3 } , , { h w 1 , h w } , , Θ .
Basic probability assignment (BPA), an independent probability-like measure, is assigned to each possible proposition to indicate belief information [40]. The BPA m ( A ) of the subset A is
m ( A ) 0 , A 2 Θ , m ( ) = 0 , A 2 Θ m ( A ) = 1 ,
where A is defined as the focal element once m ( A ) > 0 .
The evidence theory model is an uncertainty model based on evidence theory that converts the uncertain vector x = [ x 1 , x 2 , , x d ] into an evidence vector ξ = [ ξ 1 , ξ 2 , , ξ d ] composed of d independent evidence variables. A multi-index set κ that contains all index combinations of the focal element is defined as
κ = β + d : β i N i , i = 1 , 2 , , d   with   card ( κ ) = N f = i = 1 d N i ,
where + is the set of positive integers, N i is the number of focal elements of ξ i , and β is a multi-index. The joint focal element ε β , which can be regarded as a solution space, is spanned by focal elements, given by
ε β = ε β 1 × ε β 2 × ε β d
where ε β i = [ ε _ β i , ε ¯ β i ]   is the β i - th focal element of ξ i . The BPA of ε β , denoted by m ( ε β ) , is defined as the product of the BPAs of focal elements, as follows:
m ( ε β ) = i = 1 d m ( ε β i )
In the evidence theory model, the lower and upper bounds of the true extent of a proposition B , belief (Bel) and plausibility (Pl), are defined as
Bel ( B ) = A B m ( A ) ,
Pl ( B ) = A B = m ( A ) ,
Based on Bel and Pl, the cumulative belief function (CBF) and cumulative plausibility function (CPF) of system response are defined as
CBF ( u ) = Bel y ( ξ ) u = ε β y ¯ β u , β κ m ( ε β ) ,
CPF ( u ) = Pl y ( ξ ) u = ε β y _ β u , β κ m ( ε β ) ,
where y _ β and y ¯ β are the minimum and maximum values, respectively, of the system response in ε β , with u is a constant. As the lower and upper boundaries of the cumulative distribution functions (CDFs), CBF and CPF enclose the real CDF of the system response.
Other statistical properties, such as the system response expectation and variance, can be obtained. Owing to the uncertain probability distributions of the evidence variables, the expectation and variance of the system response, denoted by μ I and var I , are expressed in the following variation ranges [41]:
μ I = β κ y β I m ( ε β ) ,
var I = β κ y β I μ I 2 m ( ε β ) ,
where y β I = [ y _ β , y ¯ β ] is the response interval determined via extremum analysis, which is essentially equivalent to solving the following optimisation problem:
y _ β = min   y β   or   y ¯ β = max   y β s . t .   ξ β ε β , β κ .

3. SSR-ETM for Evidence-Theory-Based UQ

3.1. Framework of the SSR-ETM

In high-dimensional engineering applications, extremum analysis requires a large amount of computation. The size of the solution space is one of the main factors. Thus, SSR-ETM is proposed to mitigate the restriction of the size of solution space on efficiency under the precondition that the surrogate model is monotonic in one or more dimensions. As shown in Figure 1, the domain consists of solution spaces. In the SSR-ETM, the three-dimensional original solution space shown in Figure 1 can be reduced to a low-dimensional solution space. Depending on the number of monotonic dimensions, the reduced solution space may be a two-dimensional plane, line segment, or point.
The framework of the SSR-ETM is shown in Figure 2, which also presents the framework of the conventional methods in [9,29], collectively known as the evidence theory method (ETM) here, for a clear comparison. The grey parts, the blue part, and the red part are performed in the domain, reduced solution spaces, and original solution spaces, respectively. In the SSR-ETM, monotonicity analysis is performed only once, and the extremum analysis must be executed N f times.

3.2. Establishment of Sparse Representation

Suppose x = [ x 1 , x 2 , , x d ] d is a d-dimensional random vector in a continuous domain Ω , with each random variable being independent. Accordingly, the computational model y ( x ) can be expressed using PCE as follows:
y ( x ) α ƛ f α ψ α ( x ) ,
where α = [ α 1 , α 2 , , α d ] ƛ is a multi-index indicating the orders of multivariate orthogonal polynomials { ψ α ( x ) , α ƛ } , ƛ is a multi-index set, and { f α , α ƛ } are multivariate orthogonal polynomial coefficients [42]. PCE is usually truncated to a finite number of terms. The total order expansion technique [43] employed here yields an updated multi-index set ƛ d , τ , as follows:
ƛ d , τ = { α d : α τ } = { α ( 1 ) , , α ( P ) }   with   card ( ƛ d , τ ) P = ( τ + d ) ! τ ! d ! .
ψ α ( x ) is defined as the tensor product of univariate orthogonal polynomials, i.e.,
ψ α ( x ) = i = 1 d ψ α i ( x i ) ,
where ψ α i ( x i ) is a α i -order univariate orthogonal polynomial related to the probability density function (PDF) ρ ( x i ) . In addition, ψ α ( x ) satisfies orthogonality:
Ω ψ α 1 ( x ) ψ α 2 ( x ) ρ ( x ) d x = δ α 1 , α 2 ,
with δ α 1 , α 2 as the Kronecker delta and ρ ( x ) = i = 1 d ρ i ( x ) as the joint PDF.
Families of orthogonal polynomials about continuous distributions including Gaussian, uniform, gamma, and beta distributions were elaborated in [44]. However, random variables are not always limited to these classical distributions. In arbitrary PCE [45], orthogonal polynomials associated with arbitrary distributions can be constructed with the raw moment of random variables. The k-th raw moment u k of a continuous random variable is defined as
u k = Ω i x k ρ ( x ) d x .
The coefficient to the q-th power of k-order univariate orthogonal polynomial, ϖ q ( q = 1 , 2 , , k ) , is calculated by
u 0 u 1 u k u 1 u 2 u k + 1 u k 1 u k u 2 k 1 0 0 1 ϖ 0 ϖ 1 ϖ k 1 ϖ k = 0 0 0 1 .
Then, the k-order univariate orthogonal polynomial is given by
ψ k ( x ) = q = 0 k ϖ q x q .
However, the Vandermonde matrix may be ill-conditioned for a large value of k [46]. The alternative scheme relies on the following Hankel matrix [47]:
H = u 0 u 1 u k u 1 u 2 u k + 1 u k u k + 1 u 2 k
Condition det ( H ) > 0 , if satisfied, supports the Cholesky decomposition of H = R T R , with R being an upper triangular matrix. The matrix R provides information on recursive coefficients [46]
a k = r k , k + 1 r k , k r k 1 , k r k 1 , k 1 , b k = r k + 1 , k + 1 r k , k
of the following recurrence formula:
ψ 1 ( x ) = 0 , ψ 0 ( x ) = 1 , ψ k + 1 ( x ) = ( x a k ) ψ k ( x ) b k ψ k 1 ( x ) , k = 0 , 1 , 2 ,   .
In Equation (22), r i , j ( i , j = 1 , 2 , ) denotes the entry in throw i and the column j of R , with r 0 , 0 = 1 and r 0 , 1 = 0 . a 0 = μ 1 is defined here.
For a given experimental design (ED) X = { x ( 1 ) , , x ( N ) } T and the corresponding model evaluations Y = { y ( 1 ) , , y ( N ) } T , the least-squares approach can be used to solve the polynomial coefficients η = { f α , α ƛ d , τ } T as follows [48]:
η = arg min α ƛ d , τ Y Φ η 2 ; η = Φ Τ Φ 1 Φ Τ Y .
In Equation (24), Φ N × P is a measurement matrix with entries Φ i j = ψ α ( j ) ( x ( i ) ) . In [34,49,50], N 2 P , 3 P is recommended for a robust and unique η . This implies that numerous model evaluations are required for a multi-dimensional computational model with high polynomial orders. To construct a polynomial surrogate using few model evaluations, Blatman and Sudret [51] successfully applied least angle regression (LAR) to a sparse PCE to find those features with the highest impact on system response. It should be noted that the “curse of dimensionality” can be alleviated to some extent through the truncation criteria (such as the total order expansion technique mentioned above) and LAR, but it cannot be completely overcome. Referring to their work, we use sparse arbitrary polynomial chaos expansion (SAPCE) to establish sparse representations, as follows:
y ( x ) C ( x ) = α ƛ min f α ψ α ( x ) .
where ƛ min is a multi-index selected from ƛ d , τ . In SAPCE, Latin hypercube sampling [52] is used as the ED generator. Then, one can rewrite Equation (13) as follows:
C _ ( ξ β ) / C ¯ ( ξ β ) = min / max   C ( ξ β ) s . t .   ξ β ε β , β κ ,

3.3. Monotonicity Analysis

On most occasions, the optimal solution to the optimisation problem is difficult to solve analytically. Heuristic algorithms [53,54,55] are designed for desirable solutions within acceptable computational expenses to substitute optimal solutions. However, the efficiency of heuristic algorithms is inevitably hindered by dimensionality. Favourably, numerous surrogate models are monotonic in certain dimensions. Monotonicity analysis can be employed to simplify the optimisation problem. Let us use SAPCE to approximate a differentiable computational model whose partial derivative C ( ξ ) ξ i ( i = 1 , 2 , , d ) is given as
C ( ξ ) ξ i = α ƛ min f α i = 1 d ψ α i ( ξ i ) ξ i = α ƛ min f α ψ α ( ξ ) ψ α t ( ξ i ) ψ α t ( ξ i ) ξ i ,
where ψ α t ( ξ i ) ξ i follows the three-term recurrence:
ψ 0 ( ξ i ) ξ i = 0 , ψ 1 ( ξ i ) ξ i = 1 , ψ k + 1 ( ξ i ) ξ i = ψ k ( ξ i ) + ( ξ i a k ) ψ k ( ξ i ) ξ i b k ψ k 1 ( ξ i ) ξ i , k = 1 , 2 ,   .
However, the PDF of evidence variables is undetermined. In other words, only the probability assigned to each focal element is available. To construct optimal orthogonal polynomials, the authors of [29] put forward a reasonable assumption that all focal elements are assigned a uniform distribution. We employ the assumption in [29] and it yields the following PDF of evidence variables:
ρ ( ξ i ) = β i = 1 N i δ β i ( ξ i ) m ( ε β i ) ε ¯ β i ε _ β i , i = 1 , 2 , , d ,
where δ β i ( ξ i ) is an indicator function, indicating δ β i ( ξ i ) = 1 if ξ i ε β i , or else δ β i ( ξ i ) = 0 .
PSO [53] can be used to determine whether the partial derivative satisfies
min C ( ξ ) ξ i 0   or   max C ( ξ ) ξ i 0
If satisfied, then the surrogate model is an increasing or decreasing function related to ξ i ; otherwise, it is a non-monotonic function. The procedure of monotonicity analysis is shown in Algorithm 1. These three cases are mapped to p i = 1 , 2 , 3 , respectively. Therefore, the set ϑ = 1 , 2 , , d , which contains indices of the dimension, can be correspondingly split into index sets p 1 = i ϑ p i = 1 , p 2 = i ϑ p i = 2 , and p 3 = i ϑ p i = 3 . Let ξ β ( Ξ ) ε β and ξ ( Ξ ) Ω be two evidence vectors under conditions Ξ . Given ξ i ( 1 ) ξ i ( 2 ) and ξ i ( 1 ) , ξ i ( 2 ) Ω i , the monotonicity of the surrogate model makes [56,57]
C ξ ( ξ i = ξ i ( 1 ) ) C ξ ( ξ i = ξ i ( 2 ) ) ,   p i = 1   C ξ ( ξ i = ξ i ( 2 ) ) C ξ ( ξ i = ξ i ( 1 ) ) ,   p i = 2   .
Thus,
C ξ β ( ξ i = ε _ β i ) C ( ξ β ) C ξ β ( ξ i = ε ¯ β i ) ,   p i = 1 C ξ β ( ξ i = ε _ β i ) C ( ξ β ) C ξ β ( ξ i = ε ¯ β i ) ,   p i = 2 .
Algorithm 1 Monotonicity analysis
Input: Surrogate model C ( ξ ) generated by SAPCE.
1.Calculate the partial derivatives of C ( ξ ) using Equations (27) and (28).
2.Use PSO to calculate the extremum of partial derivatives of each dimension.
3.Assess the monotonicity of the surrogate model using Equation (30).
Output: Monotonicity of the surrogate model.
In essence, monotonicity analysis uses PSO to solve the extremum of the partial derivative of the surrogate model d times. Because the extremum calculated by PSO must not be a global extremum, it may lead to a false monotonicity judgment when the minimum or maximum of the above partial derivative is very close to zero, resulting in a large error in the subsequent calculation. Generally speaking, however, the global optimization performance of PSO has been widely verified.

3.4. Solution Space Reduction for Extremum Analysis of Joint Focal Elements

The monotonicity of the surrogate model facilitates in narrowing the solution space. Let r 1 , r 2 , and r 3 be the number of elements in p 1 , p 2 , and p 3 , respectively. Assume p 1 , p 2 , and p 3 divide any vector or set into three parts. It can be deduced from Equation (32) that
C ξ β ( ξ β p 1 = ε _ β p 1 ) C ( ξ β ) C ξ β ( ξ β p 1 = ε ¯ β p 1 ) , C ξ β ( ξ β p 2 = ε _ β p 2 ) C ( ξ β ) C ξ β ( ξ β p 2 = ε ¯ β p 2 ) ,
where ξ β ( ξ β p 1 = ε _ β p 1 ) , as an example, represents
ξ β ( ξ β p 1 = ε _ β p 1 ) = ξ β ( ξ β p 1 ( 1 ) = ε _ β p 1 ( 1 ) , ξ β p 1 ( 2 ) = ε _ β p 1 ( 2 ) , , ξ β p 1 ( r 1 ) = ε _ β p 1 ( r 1 ) )
with β p 1 ( r 1 ) being the r 1 - th index of β p 1 . The combination of the two inequalities in Equation (33) yields
C ξ β ( ξ β p 1 = ε _ β p 1 , ξ β p 2 = ε ¯ β p 2 ) C ( ξ β ) C ξ β ( ξ β p 1 = ε ¯ β p 1 , ξ β p 2 = ε _ β p 2 ) .
Equation (35) indicates that the extremums of the optimization problems are located in the following reduced solution spaces:
ε β min = ε _ β p 1 × ε ¯ β p 2 × ε β p 3 , ε β max = ε ¯ β p 1 × ε _ β p 2 × ε β p 3 ,
where ε β min and ε β max denote the reduced solution spaces corresponding to the minimisation/maximisation problem, respectively. The evidence vectors ξ β min and ξ β max generated in ε β min and ε β max , respectively, are expressed as
ξ β min / ξ β max = [ ξ β 1 , ξ β 2 , , ξ β d ] ,   where   ξ β i ( i = 1 , 2 , , d ) = ε _ β i / ε ¯ β i ,   p i = 1   ε ¯ β i / ε _ β i , p i = 2 ξ β i ε β i ,   p i = 3 .
For instance, provided that a four-dimensional surrogate model is increasing in the first and fourth dimensions, decreasing in the second dimension, and non-monotonic in the third dimension, then, based on Equation (36), the original solution space ε β is reduced to ε β min = ε _ β 1 × ε ¯ β 2 × ε β 3 × ε _ β 4 or ε β max = ε ¯ β 1 × ε _ β 2 × ε β 3 × ε ¯ β 4 . The new solution spaces are composed of only the third dimension. ξ β min and ξ β max are given by ξ β min = [ ε _ β 1 , ε ¯ β 2 , ξ β 3 , ε _ β 4 ] and ξ β max = [ ε ¯ β 1 , ε _ β 2 , ξ β 3 , ε ¯ β 4 ] , respectively. In addition, we can find ξ β p 1 = ε ¯ β 2 , ξ β p 2 = [ ε _ β 1 , ε _ β 4 ] , and ξ β p 3 = ξ β 3 . Thus, ξ β p 1 and ξ β p 2 contain all constants, and ξ β p 3 contains all unknown evidence variables, implying ψ α p 1 ( ξ β p 1 ) and ψ α p 2 ( ξ β p 2 ) as constants, and ψ α p 3 ( ξ β p 3 ) as the part that should be solved iteratively. This example intuitively illustrates the reduction in solution space.
The reduction in solution space depends on the number of monotonic dimensions of the surrogate model. If r 3 = d , the SSR-ETM is transformed into ETM, implying that the solution space cannot be reduced. In this case, the SSR-ETM takes additional time for monotonicity analysis, which is also the limitation of the SSR-ETM. Nevertheless, the SSR-ETM effectively narrows the size of the solution space from d to r 3 as long as r 3 < d . Moreover, the solution space is reduced to a point when r 3 = 0 .

3.5. Evidence-Theory-Based Response Analysis Using the Solution Space Reduction Technique

The extremum of the system response in each joint focal element should be determined for evidence-theory-based UQ. Suppose the monotonicity of the surrogate model determined by PSO is accurate, then, based on the inference in Section 3.4 and the solution space reduction technique, the optimisation problem in Equation (26) is exactly equivalent to
C _ ( ξ β ) / C ¯ ( ξ β ) = min/max α ƛ min f α ψ α ( ξ β ) ,   r 3 = d α ƛ min f α ψ α p 1 ( ξ β p 1 ) ψ α p 2 ( ξ β p 2 ) ψ α p 3 ( ξ β p 3 ) ,   0 < r 3 < d α ƛ min f α ψ α p 1 ( ξ β p 1 ) ψ α p 2 ( ξ β p 2 )   r 3 = 0 s . t . ξ β ε β , β k
This optimisation problem can be solved using PSO. Compared with Equation (26), the surrogate model in Equation (38) is divided into three parts according to monotonicity. On the premise that the monotonicity judgment is correct, the theoretical results of these two equations are exactly the same, indicating the same calculation accuracy of the SSR-ETM and ETM. Considering the randomness of PSO and the change in solution space, there may be a slight deviation between the calculation results of these two methods. To sum up, the SSR-ETM relies on the solution space reduction technique to improve computational efficiency by reducing the solution space in extremum analysis, and the computational accuracy is almost the same as that of the ETM.
After calculating C _ ( ξ β ) and C ¯ ( ξ β ) , four indicators can be rewritten as
CBF ( u ) = ε β C ¯ ( ξ β ) u , β κ m ( ε β ) ,
CPF ( u ) = ε β C _ ( ξ β ) u , β κ m ( ε β ) ,
μ I = β κ C I ( ξ β ) m ( ε β ) ,
var I = β κ C I ( ξ β ) μ I 2 m ( ε β ) ,
where C I ( ξ β ) = C _ ( ξ β ) , C ¯ ( ξ β ) .
Algorithm 2 shows the pseudo-code of the SSR-ETM. In Algorithm 2, ψ α p 1 ( ξ β p 1 ) and ψ α p 2 ( ξ β p 2 ) remain unchanged in the entire loop, and the value of ψ α p 3 ( ξ β p 3 ) depends on ξ β p 3 generated by PSO.
Algorithm 2 Solution-space-reduction-based evidence theory method (SSR-ETM)
Input: Monotonicity of the surrogate model determined by Algorithm 1, joint focal elements (solution spaces) ε β , and BPAs m ( ε β ) , β κ .
1. For   β κ
2.       If  r 3 < d
3.       Use Equation (37) to generate ξ β p 1 and ξ β p 2
4.       Use Equation (16) to calculate ψ α p 1 ( ξ β p 1 ) and ψ α p 2 ( ξ β p 2 ) .
5. End
6.       While PSO has not searched the extremum of C ( ξ β )  and r 3 > 0
7.          Use PSO to generate ξ β p 3 .
8.           Calculate ψ α p 3 ( ξ β p 3 ) using Equation (16).
9.          Calculate C ( ξ β ) .
10.        End
11.       Store C _ ( ξ β ) / C ¯ ( ξ β ) .
12. End
13. Calculate the four indicators using Equations (39)–(42).
Output: CBF, CPF, μ I , and var I .
Essentially, both the SSR-ETM and ETM use PSO to calculate the extremum of the surrogate model 2 N f times. The difference is that the particle dimension of PSO in the SSR-ETM is r 3 , while the particle dimension in the ETM is d . Therefore, when N f is relatively large and r 3 < d , the computational cost of the SSR-ETM for monotonicity analysis is relatively small, thus it is more efficient than the ETM.

4. Numerical Examples

In this section, four simple functions and an air spring system with epistemic uncertainty are used to compare the SSR-ETM and the ETM. The differences between the two methods are discussed in Section 3. The results were acquired using MATLAB R2021b on a 3.70 GHz AMD Ryzen Threadripper 3970X 32-Core CPU. The PSO used the built-in function “Particleswarm” of MATLAB, and all other parameters were default values, except for keeping the “Vectorized” on. The modified “leave-one-out” error ε was used as the error metrics; more details can be found in [51]. We also used some open-source code from UQLab [58]. The necessary parameters will be illustrated in each numerical example.

4.1. Simple Functions

Let us consider the following function:
y ( ξ ) = I m i = min { 1 , m } m ( ξ i + 2 ) 4 + I n j = min { 1 , n } n ( ξ j 0.5 ) 4 ,
where ξ i , ξ j [ 1 , 1 ] are independent evidence variables; I m and I n are defined as
I m = 0 , m = 0 1 m , m 0 , I n = 0 , n = 0 1 n , n 0 .
The partial derivative of the above function is given as
y ξ l = 4 ( ξ l + 2 ) 3 m , l = i > 0 4 ( ξ l 0.5 ) 3 n , l = j > 0 .
Thus, the function is monotonic in m dimensions and non-monotonic in n dimensions. Four functions derived from the function in Equation (43) are designed to illustrate the superiority of the SSR-ETM. In these four functions, N BPA focal elements with the same BPAs are allocated to each evidence variable. For example, the focal elements of an evidence variable in N BPA = 2 are set to [−1,0] and [0,1], with BPAs of 50% and 50%, respectively. Information on the four functions is presented in Table 1. In the simple functions, the number of samples was N = 60 ( m + n ) , the “leave-one-out” error was ε = 1 × 10 20 , and the highest order of the polynomial was eight. We set ε so small because there are many items in the candidate base for high-dimensional functions. In practical engineering problems, 10 4 ~ 10 3 is appropriate because PCE usually contains very few lower order terms.
The results showed that the four indicators about these simple functions obtained by the two methods are almost identical to the analytical solutions, indicating the high computational accuracy of the two methods. Because the plotted curves are almost coincident (similar to Figures 8 and 9) and limited by the fact that CPF and CBF are not suitable for drawing in the form of relative error, these indicators of simple functions are not exhibited here to improve readability and to emphasize the advantages of the SSR-ETM in terms of computational efficiency.

4.1.1. Non-Monotonic Function

The non-monotonic function shown in Equation (46) was used to study the performance of the SSR-ETM when all dimensions are non-monotonic. In this case, the SSR-ETM transforms into the ETM. This implies that the solution space remains unchanged.
y ( ξ ) = 1 n i = 1 n ( ξ i 0.5 ) 4 , n = 1 , 2 , , 20
As depicted in Figure 3, the SSR-ETM requires additional time for monotonicity analysis and solution space reduction, yet it is entirely acceptable. The SSR-ETM takes only approximately 0.54 s longer than the ETM when the number of non-monotonic dimensions is 20. Therefore, the SSR-ETM features a low-risk cost. In addition, the computation time of both methods rapidly increased with the number of dimensions. This is the limitation of the ETM in dealing with high-dimensional uncertainty problems. The performance of the SSR-ETM is not optimistic in this special case, either. However, it should be pointed out that such a special case rarely occurs in high-dimensional engineering problems.

4.1.2. Monotonic Function

To assess the performance of the SSR-ETM under the assumption that all dimensions are monotonic, a monotonic function was designed as follows:
y ( ξ ) = 1 m i = 1 m ( ξ i + 2 ) 4 , m = 1 , 2 , , 20
In this case, the SSR-ETM precisely locates the minimum and maximum points via monotonicity analysis. Accordingly, we can see from Figure 4 that the efficiency of the SSR-ETM is drastically improved compared with that of the ETM when m assumes a large value. Although the SSR-ETM requires additional time for monotonicity analysis and solution space reduction, it is still faster than the ETM. The SSR-ETM requires no more than 580 s to deal with the uncertainty problem containing 2 20 joint focal elements, while the ETM requires more than 90,000 s. Moreover, the SSR-ETM remains powerful despite the increase in monotonic dimension, while the computation time of the ETM grows explosively. Because of the short calculation time required for UQ, a small turn appeared around m = 3 .

4.1.3. Fifteen-Dimensional Function

Equation (48) expresses the 15-dimensional function designed to investigate the influence of the number of monotonic dimensions on computing efficiency.
y ( ξ ) = I m i = I m m ( ξ i + 2 ) 4 + I n j = I n n ( ξ j 0.5 ) 4 , m + n = 15 , m 0 , n 0
As shown in Figure 5, the computation time of the SSR-ETM is notably shortened as the number of monotonic dimensions increases, owing to the reduction in the solution space. In contrast, the computation time of ETM fluctuates within a certain range. The SSR-ETM can save approximately 2175 s for a completely monotonic function compared with a completely non-monotonic function; this comparison result is valid only for a function with a dimension of 15 and a number of joint focal elements of 2 15 . Thus, the efficiency improvement realized by solution space reduction is significant. Although the time invested in monotonicity analysis and solution space reduction may add additional time overall, it is worthwhile.

4.1.4. Five-Dimensional Function

Equation (49) expresses the five-dimensional function designed to illustrate the impact of the number of joint focal elements on the computational efficiency advantage of the SSR-ETM.
y ( ξ ) = I m i = min { 1 , m } m ( ξ i + 2 ) 4 + I n j = min { 1 , n } n ( ξ j 0.5 ) 4 , m + n = 5 ,   m 1 ,   n 0
Figure 6 reports the computation time ratio of the SSR-ETM to the ETM for different m and N BPA . This illustrates that the computation time ratio of the two methods tends to be stable when the number of joint focal elements reaches a certain degree, implying that the efficiency advantage of the SSR-ETM is mainly affected by the number of monotonic dimensions rather than the number of joint focal elements. Owing to the stability of the SSR-ETM, it is suitable for dealing with epistemic uncertainty problems with numerous joint focal elements.

4.1.5. Summary of the Results of the Four Simple Functions

The results obtained for the above four simple functions support the following conclusions.
(i)
The SSR-ETM features high computational accuracy that is almost the same as that of the ETM because it only divides the surrogate model into monotonic and non-monotonic parts.
(ii)
Although the SSR-ETM requires a minuscule amount of time for monotonicity analysis and solution space reduction, it may render remarkable benefits, specifically when the surrogate model is monotonic in all dimensions. When the surrogate model is non-monotonic in all dimensions, the SSR-ETM maintains almost the same computing efficiency as the ETM.
(iii)
The SSR-ETM performs better when there are more monotonic dimensions.
(iv)
The number of joint focal elements exhibits a negligible effect on the calculation efficiency advantage of the SSR-ETM over the ETM when it reaches a certain degree.

4.2. Air Spring System with Epistemic Uncertainty

4.2.1. Finite Element Model (FEM) of an Air Spring System

Figure 7 shows the air spring finite element model (FEM). The air spring system comprises an upper cover, lower seat, capsule, and cord layers. A total of 119,360 elements were used in modelling, including 30,960 C3D8H elements for the capsule, 35,360 C3D8R elements for the upper cover and lower seat, 49,360 SFM3D4R elements for the cord layer modelled by the rebar layer, and 80 C3D6H elements and 3600 C3D6 elements for the junction between the upper cover and the capsule. The capsule is made of a hyper-elastic material whose constitutive relation is characterised by the Mooney–Rivlin model [59,60].

4.2.2. Air Springs with Uncertain Parameters

Air springs are involved in these uncertainties. The uncertain vector Z , which is composed of five uncertain parameters of air springs that are assumed to be independent, is expressed as follows:
Z = [ E , φ , l , S , δ ] T ,
where E , φ , l , and S are the Young’s modulus, angle, spacing, and cross-sectional area of the cords, respectively, and δ is the thickness of the capsule. The evidence theory model defines these uncertain parameters as evidence variables, whose focal elements and corresponding BPAs are listed in Table 2. The system response K , the vertical static stiffness of air springs, can be obtained by solving the following finite element equation:
K = M ( Z ) .

4.2.3. SSR-ETM for Stiffness Evaluation

For a given ED X = { Z ( 1 ) , , Z ( N ) } T and its corresponding model evaluations Y = { K ( 1 ) , , K ( N ) } T , one can use the two methods to calculate the bounds of expectation and variance of the vertical static stiffness at different initial pressures. CBF and CPF at different initial pressures are also available. Owing to fewer uncertain variables of the air spring system, we incorporated 10 5 random samples and the vertex method [61] to obtain the extremum in each joint focal element and the reference solutions of the four indicators. Note that the reference solutions were obtained based on the surrogate model rather than the computational model or FEM. In this numerical example, the number of samples was 200, the “leave-one-out” error was 1 × 10 3 , and the highest order of the polynomial was 12. The results of the stiffness evaluation are presented in Figure 8, Figure 9 and Figure 10.
We can observe from Figure 8 that the bounds of the expectation and variance of the vertical static stiffness increase with the initial pressure. This is because the stress on the air springs increases with the initial pressure, leading to a higher vertical static stiffness owing to the adjustable stiffness characteristics.
As depicted in Figure 9, the curves of CPF are close to the CBF. This is because the minimum in each joint focal element is close to that of the maximum. The curves of the four indicators obtained by the two methods are almost coincident with that of the reference solutions, proving both methods provide accurate results of stiffness evaluation at different initial pressures as long as the surrogate model is well established.
Accelerated by the reduction in solution space, the SSR-ETM spends less time than the ETM on UQ (see Figure 10). Particularly, the surrogate model of the air spring system is monotonic in four dimensions when the initial pressure is between 0.46 MPa and 0.54 MPa. The SSR-ETM saves about 1400 s for this uncertain system with only five uncertain variables. The SSR-ETM still performs well under 0.56–0.68 MPa as two or three dimensions are reduced. Thus, the SSR-ETM was proved to be valid.

5. Conclusions

In this study, we propose the SSR-ETM. Initially, SAPCE was used to approximate the computational model. Thereafter, PSO was employed for monotonicity analysis. Based on monotonicity analysis, the solution space can be reduced if the surrogate model is monotonic in at least one dimension. The reduced solution space is conducive to accelerating extremum analysis. Ultimately, four simple functions and an air-spring system with epistemic uncertainty successfully verified the validity of the SSR-ETM. Our results support the following conclusions:
(i)
The SSR-ETM demonstrates a high computing accuracy almost comparable to that of the ETM as long as the surrogate model is well established.
(ii)
Compared with the ETM, the SSR-ETM adds minimal additional time for monotonicity analysis and solution space reduction.
(iii)
More monotonic dimensions contribute to a higher efficiency advantage of the SSR-ETM. In particular, when all dimensions are monotone, the SSR-ETM exhibits a significant efficiency advantage over the ETM.
(iv)
The SSR-ETM performed better than the ETM in the stiffness evaluation of the air springs with epistemic uncertainty.
In conclusion, the SSR-ETM is a promising UQ method with low risk and high return, particularly suitable for engineering applications with prominent monotonicity. However, the SSR-ETM still does not perform well for the high-dimensional surrogate model without any monotonic dimension or with few monotonic dimensions. In addition, it requires very high accuracy of the surrogate model. Moreover, the number of extremum analyses in SSR-ETM multiplies with the number of focal elements contained in the number of uncertain variables. Although we have listed several numerical examples with 2 20 joint focal elements, the number of extremum analyses in SSR-ETM is the product of the number of focal elements in all variables, which makes it difficult for the SSR-ETM to deal with problems with many uncertainty variables. In future work, it is necessary to study further how to solve these problems.

Author Contributions

Conceptualization, S.Y. and K.J.; methodology, S.Y.; software, K.J.; validation, S.Y. and K.J.; formal analysis, K.J.; investigation, Y.B.; resources, S.Y., W.Z. and Z.W.; data curation, K.J. and Y.B.; writing—original draft preparation, K.J.; writing—review and editing, K.J.; visualization, K.J.; supervision, S.Y.; project administration, S.Y., W.Z. and Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Please send an email to [email protected] if any questions about the code associated with this article. The air spring model is confidential and will not be available.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, J.-J.; Yin, Z.-H.; Yuan, X.-J.; Qiu, G.-Q.; Guo, K.-H.; Wang, X.-L. A refined stiffness model of rolling lobe air spring with structural parameters and the stiffness characteristics of rubber bellows. Measurement 2021, 169, 108355. [Google Scholar] [CrossRef]
  2. Chen, J.-J.; Yin, Z.-H.; Rakheja, S.; He, J.-H.; Guo, K.-H. Theoretical modelling and experimental analysis of the vertical stiffness of a convoluted air spring including the effect of the stiffness of the bellows. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2018, 232, 547–561. [Google Scholar] [CrossRef]
  3. Hora, S.C. Aleatory and epistemic uncertainty in probability elicitation with an example from hazardous waste management. Reliab. Eng. Syst. Saf. 1996, 54, 217–223. [Google Scholar] [CrossRef]
  4. Hüllermeier, E.; Waegeman, W. Aleatoric and epistemic uncertainty in machine learning: An introduction to concepts and methods. Mach. Learn. 2021, 110, 457–506. [Google Scholar] [CrossRef]
  5. Salehi, S.; Raisee, M.; Cervantes, M.J.; Nourbakhsh, A. An efficient multifidelity ℓ1-minimization method for sparse polynomial chaos. Comput. Methods Appl. Mech. Eng. 2018, 334, 183–207. [Google Scholar] [CrossRef]
  6. Tarakanov, A.; Elsheikh, A.H. Regression-based sparse polynomial chaos for uncertainty quantification of subsurface flow models. J. Comput. Phys. 2019, 399, 108909. [Google Scholar] [CrossRef] [Green Version]
  7. Thapa, M.; Mulani, S.B.; Walters, R.W. A new non-intrusive polynomial chaos using higher order sensitivities. Comput. Methods Appl. Mech. Eng. 2018, 328, 594–611. [Google Scholar] [CrossRef]
  8. Torre, E.; Marelli, S.; Embrechts, P.; Sudret, B. Data-driven polynomial chaos expansion for machine learning regression. J. Comput. Phys. 2019, 388, 601–623. [Google Scholar] [CrossRef] [Green Version]
  9. Yin, S.; Yu, D.; Luo, Z.; Xia, B. An arbitrary polynomial chaos expansion approach for response analysis of acoustic systems with epistemic uncertainty. Comput. Methods Appl. Mech. Eng. 2018, 332, 280–302. [Google Scholar] [CrossRef]
  10. Yin, S.; Yu, D.; Yin, H.; Lü, H.; Xia, B. Hybrid evidence theory-based finite element/statistical energy analysis method for mid-frequency analysis of built-up systems with epistemic uncertainties. Mech. Syst. Signal Process. 2017, 93, 204–224. [Google Scholar] [CrossRef]
  11. Yin, S.; Yu, D.; Ma, Z.; Xia, B. A unified model approach for probability response analysis of structure-acoustic system with random and epistemic uncertainties. Mech. Syst. Signal Process. 2018, 111, 509–528. [Google Scholar] [CrossRef]
  12. Pantelides, C.P.; Ganzerli, S. Comparison of fuzzy set and convex model theories in structural design. Mech. Syst. Signal Process. 2001, 15, 499–511. [Google Scholar] [CrossRef]
  13. Yang, S.; Han, J.; Xia, L.; Chen, Y.-H. An optimal fuzzy-theoretic setting of adaptive robust control design for a lower limb exoskeleton robot system. Mech. Syst. Signal Process. 2020, 141, 106706. [Google Scholar] [CrossRef]
  14. Huang, F.; Sava, A.; Adjallah, K.H.; Wang, Z. Fuzzy model identification based on mixture distribution analysis for bearings remaining useful life estimation using small training data set. Mech. Syst. Signal Process. 2021, 148, 107173. [Google Scholar] [CrossRef]
  15. van Mierlo, C.; Burmberger, L.; Daub, M.; Duddeck, F.; Faes, M.G.; Moens, D. Interval methods for lack-of-knowledge uncertainty in crash analysis. Mech. Syst. Signal Process. 2022, 168, 108574. [Google Scholar] [CrossRef]
  16. Zhang, W.; Liu, J.; Cho, C.; Han, X. A hybrid parameter identification method based on Bayesian approach and interval analysis for uncertain structures. Mech. Syst. Signal Process. 2015, 60, 853–865. [Google Scholar] [CrossRef]
  17. Deng, Z.; Guo, Z.; Zhang, X. Interval model updating using perturbation method and radial basis function neural networks. Mech. Syst. Signal Process. 2017, 84, 699–716. [Google Scholar] [CrossRef]
  18. Fang, S.-E.; Zhang, Q.-H.; Ren, W.-X. An interval model updating strategy using interval response surface models. Mech. Syst. Signal Process. 2015, 60, 909–927. [Google Scholar] [CrossRef]
  19. Zhu, W.; Chen, N.; Liu, J.; Beer, M. A probability-box-based method for propagation of multiple types of epistemic uncertainties and its application on composite structural-acoustic system. Mech. Syst. Signal Process. 2021, 149, 107184. [Google Scholar] [CrossRef]
  20. Crespo, L.G.; Kenny, S.P.; Giesy, D.P. Reliability analysis of polynomial systems subject to p-box uncertainties. Mech. Syst. Signal Process. 2013, 37, 121–136. [Google Scholar] [CrossRef]
  21. Su, X.; Mahadevan, S.; Xu, P.; Deng, Y. Dependence Assessment in Human Reliability Analysis Using Evidence Theory and AHP. Risk Anal. 2015, 35, 1296–1316. [Google Scholar] [CrossRef] [PubMed]
  22. 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]
  23. Shamlou, A.; Reza Feyzi, M.; Behjat, V. Interpretation of frequency response analysis of power transformer based on evidence theory. IET Gener. Transm. Distrib. 2019, 13, 3879–3887. [Google Scholar] [CrossRef]
  24. Amini, A.; Abdollahi, A.; Hariri-Ardebili, M.A.; Lall, U. Copula-based reliability and sensitivity analysis of aging dams: Adaptive Kriging and polynomial chaos Kriging methods. Appl. Soft Comput. 2021, 109, 107524. [Google Scholar] [CrossRef]
  25. Hariri-Ardebili, M.A.; Mahdavi, G.; Abdollahi, A.; Amini, A. An RF-PCE hybrid surrogate model for sensitivity analysis of dams. Water 2021, 13, 302. [Google Scholar] [CrossRef]
  26. Hariri-Ardebili, M.A.; Sudret, B. Polynomial chaos expansion for uncertainty quantification of dam engineering problems. Eng. Struct. 2020, 203, 109631. [Google Scholar] [CrossRef]
  27. Li, Y.; Hariri-Ardebili, M.A.; Deng, T.; Wei, Q.; Cao, M. A surrogate-assisted stochastic optimization inversion algorithm: Parameter identification of dams. Adv. Eng. Inform. 2023, 55, 101853. [Google Scholar] [CrossRef]
  28. Shah, H.R.; Hosder, S.; Winter, T. Quantification of Margins and Mixed Uncertainties Using Evidence Theory and Stochastic Expansions. In Proceedings of the 16th AIAA Non-Deterministic Approaches Conference, National Harbor, MD, USA, 13–17 January 2014. [Google Scholar]
  29. Yin, S.; Yu, D.; Yin, H.; Xia, B. A new evidence-theory-based method for response analysis of acoustic system with epistemic uncertainty by using Jacobi expansion. Comput. Methods Appl. Mech. Eng. 2017, 322, 419–440. [Google Scholar] [CrossRef]
  30. Zhang, J.; Gong, W.; Yue, X.; Shi, M.; Chen, L. Efficient reliability analysis using prediction-oriented active sparse polynomial chaos expansion. Reliab. Eng. Syst. Saf. 2022, 228, 108749. [Google Scholar] [CrossRef]
  31. Zhang, B.-Y.; Ni, Y.-Q. A hybrid sequential sampling strategy for sparse polynomial chaos expansion based on compressive sampling and Bayesian experimental design. Comput. Methods Appl. Mech. Eng. 2021, 386, 114130. [Google Scholar] [CrossRef]
  32. Zhang, J.; Yue, X.; Qiu, J.; Zhuo, L.; Zhu, J. Sparse polynomial chaos expansion based on Bregman-iterative greedy coordinate descent for global sensitivity analysis. Mech. Syst. Signal Process. 2021, 157, 107727. [Google Scholar] [CrossRef]
  33. Zhou, Y.; Lu, Z.; Cheng, K.; Ling, C. An efficient and robust adaptive sampling method for polynomial chaos expansion in sparse Bayesian learning framework. Comput. Methods Appl. Mech. Eng. 2019, 352, 654–674. [Google Scholar] [CrossRef]
  34. Luthen, N.; Marelli, S.; Sudret, B. Sparse polynomial chaos expansions: Literature survey and benchmark. SIAM/ASA J. Uncertain. Quantif. 2021, 9, 593–649. [Google Scholar] [CrossRef]
  35. Dempster, A.P. A generalization of Bayesian inference. J. R. Stat. Soc. Ser. B 1968, 30, 205–232. [Google Scholar] [CrossRef]
  36. Shafer, G. A Mathematical Theory of Evidence; Princeton University Press: Princeton, NJ, USA, 1976; Volume 42. [Google Scholar]
  37. Shafer, G.; Pearl, J. Readings in Uncertain Reasoning; Morgan Kaufmann Publishers Inc.: Burlington, MA, USA, 1990. [Google Scholar]
  38. Shafer, G. Perspectives on the theory and practice of belief functions. Int. J. Approx. Reason. 1990, 4, 323–362. [Google Scholar] [CrossRef] [Green Version]
  39. Shafer, G. Dempster-shafer theory. Encycl. Artif. Intell. 1992, 1, 330–331. [Google Scholar]
  40. Wang, C.; Matthies, H.G. Evidence theory-based reliability optimization design using polynomial chaos expansion. Comput. Methods Appl. Mech. Eng. 2018, 341, 640–657. [Google Scholar] [CrossRef]
  41. Bai, Y.C.; Jiang, C.; Han, X.; Hu, D.A. Evidence-theory-based structural static and dynamic response analysis under epistemic uncertainties. Finite Elem. Anal. Des. 2013, 68, 52–62. [Google Scholar] [CrossRef]
  42. Zhou, Y.; Lu, Z.; Cheng, K. Adaboost-based ensemble of polynomial chaos expansion with adaptive sampling. Comput. Methods Appl. Mech. Eng. 2022, 388, 114238. [Google Scholar] [CrossRef]
  43. Sudret, B. Global sensitivity analysis using polynomial chaos expansions. Reliab. Eng. Syst. Saf. 2008, 93, 964–979. [Google Scholar] [CrossRef]
  44. 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]
  45. Oladyshkin, S.; Nowak, W. Data-driven uncertainty quantification using the arbitrary polynomial chaos expansion. Reliab. Eng. Syst. Saf. 2012, 106, 179–190. [Google Scholar] [CrossRef]
  46. Guo, L.; Liu, Y.; Zhou, T. Data-driven polynomial chaos expansions: A weighted least-square approximation. J. Comput. Phys. 2019, 381, 129–145. [Google Scholar] [CrossRef] [Green Version]
  47. Ahlfeld, R.; Belkouchi, B.; Montomoli, F. SAMBA: Sparse approximation of moment-based arbitrary polynomial chaos. J. Comput. Phys. 2016, 320, 1–16. [Google Scholar] [CrossRef] [Green Version]
  48. Thapa, M.; Mulani, S.B.; Walters, R.W. Adaptive weighted least-squares polynomial chaos expansion with basis adaptivity and sequential adaptive sampling. Comput. Methods Appl. Mech. Eng. 2020, 360, 112759. [Google Scholar] [CrossRef]
  49. Fajraoui, N.; Marelli, S.; Sudret, B. Sequential Design of Experiment for Sparse Polynomial Chaos Expansions. SIAM/ASA J. Uncertain. Quantif. 2017, 5, 1061–1085. [Google Scholar] [CrossRef]
  50. Hosder, S.; Walters, R.; Balch, M. Efficient sampling for non-intrusive polynomial chaos applications with multiple uncertain input variables. In Proceedings of the 48th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Honolulu, HI, USA, 23–26 April 2007; p. 1939. [Google Scholar]
  51. Blatman, G.; Sudret, B. Adaptive sparse polynomial chaos expansion based on least angle regression. J. Comput. Phys. 2011, 230, 2345–2367. [Google Scholar] [CrossRef]
  52. McKay, M.D.; Beckman, R.J.; Conover, W.J. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 2000, 42, 55–61. [Google Scholar] [CrossRef]
  53. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95-International Conference on Neural Networks; Perth, WA, Australia, 27 November–1 December 1995, pp. 1942–1948.
  54. Holland, J.H. Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence; MIT Press: Cambridge, MA, USA, 1992. [Google Scholar]
  55. Dorigo, M.; Birattari, M.; Stutzle, T. Ant colony optimization. IEEE Comput. Intell. Mag. 2006, 1, 28–39. [Google Scholar] [CrossRef]
  56. Kriz, I.; Pultr, A. Multivariable Differential Calculus. In Introduction to Mathematical Analysis; Springer: Basel, Switzerland, 2013; pp. 65–95. [Google Scholar] [CrossRef]
  57. O’Regan, G. Calculus. In Mathematics in Computing: An Accessible Guide to Historical, Foundational and Application Contexts; Springer: London, UK, 2013; pp. 247–265. [Google Scholar] [CrossRef]
  58. Marelli, S.; Sudret, B. UQLab: A Framework for Uncertainty Quantification in Matlab; American Society of Civil Engineers: Reston, VA, USA, 2014. [Google Scholar]
  59. Mooney, M. A Theory of Large Elastic Deformation. J. Appl. Phys. 1940, 11, 582–592. [Google Scholar] [CrossRef]
  60. Rivlin, R.S. Large elastic deformations of isotropic materials IV. Further developments of the general theory. Philos. Trans. R. Soc. London. Ser. A Math. Phys. Sci. 1948, 241, 379–397. [Google Scholar]
  61. Qiu, Z.; Xia, Y.; Yang, J. The static displacement and the stress analysis of structures with bounded uncertainties using the vertex solution theorem. Comput. Methods Appl. Mech. Eng. 2007, 196, 4965–4984. [Google Scholar] [CrossRef]
Figure 1. Domain, original solution space, and its reduced solution spaces.
Figure 1. Domain, original solution space, and its reduced solution spaces.
Mathematics 11 01214 g001
Figure 2. Framework of the SSR-ETM and ETM.
Figure 2. Framework of the SSR-ETM and ETM.
Mathematics 11 01214 g002
Figure 3. Computation time of the two methods for the non-monotonic function.
Figure 3. Computation time of the two methods for the non-monotonic function.
Mathematics 11 01214 g003
Figure 4. Computation time of the two methods for the monotonic function.
Figure 4. Computation time of the two methods for the monotonic function.
Mathematics 11 01214 g004
Figure 5. Computation time of the two methods for the 15-dimensional function.
Figure 5. Computation time of the two methods for the 15-dimensional function.
Mathematics 11 01214 g005
Figure 6. Computation time ratio of the SSR-ETM to the ETM for the five-dimensional function.
Figure 6. Computation time ratio of the SSR-ETM to the ETM for the five-dimensional function.
Mathematics 11 01214 g006
Figure 7. FEM of an air spring.
Figure 7. FEM of an air spring.
Mathematics 11 01214 g007
Figure 8. Bounds of the expectation (a) and variance (b) of the vertical static stiffness and their reference solutions at different initial pressures.
Figure 8. Bounds of the expectation (a) and variance (b) of the vertical static stiffness and their reference solutions at different initial pressures.
Mathematics 11 01214 g008
Figure 9. CBF and CPF of the two methods and their reference solutions for the air spring system.
Figure 9. CBF and CPF of the two methods and their reference solutions for the air spring system.
Mathematics 11 01214 g009
Figure 10. Computation time of the two methods at different initial pressures.
Figure 10. Computation time of the two methods at different initial pressures.
Mathematics 11 01214 g010
Table 1. Information on the four functions.
Table 1. Information on the four functions.
Function Values   of   m   and   n N BPA
Non-monotonic function m = 0 ,   n = 1 , 2 , , 20 2
Monotonic function m = 1 , 2 , , 20 ,   n = 0 2
15-dimensional function m + n = 15 ,   m 0 ,   n 0 2
5-dimensional function m + n = 5 ,   m 1 ,   n 0 2, 3, …, 8
Table 2. Information about the air spring system.
Table 2. Information about the air spring system.
E (MPa) φ (Degree) l (mm) S (mm2) δ (mm)
Focal ElementBPA (%)Focal ElementBPA (%)Focal ElementBPA (%)Focal ElementBPA (%)Focal ElementBPA (%)
[1657,1675]3.90[40.8,41.5]0.85[1.40,1.42]5.05[5.4,5.5]0.47[6.0,6.2]0.65
[1675,1700]7.64[41.5,42.0]1.07[1.42,1.44]7.28[5.5,5.6]0.99[6.2,6.4]2.12
[1700,1725]7.11[42.0,42.5]3.06[1.44,1.46]8.33[5.6,5.7]4.74[6.4,6.6]25.36
[1725,1750]5.62[42.5,43.0]4.13[1.46,1.48]17.95[5.7,5.8]9.67[6.6,6.8]20.36
[1750,1775]7.28[43.0,43.5]5.11[1.48,1.50]14.63[5.8,5.9]20.97[6.8,7.0]21.21
[1775,1800]7.67[43.5,44.0]9.67[1.50,1.52]19.51[5.9,6.0]25.98[7.0,7.2]18.54
[1800,1825]7.78[44.0,44.5]8.75[1.52,1.54]11.46[6.0,6.1]20.86[7.2,7.4]5.10
[1825,1850]9.04[44.5,45.0]14.68[1.54,1.56]8.40[6.1,6.2]9.62[7.4,7.6]4.33
[1850,1875]7.25[45.0,45.5]12.61[1.56,1.58]5.98[6.2,6.3]4.95[7.6,7.8]1.98
[1875,1900]5.38[45.5,46.0]8.92[1.58,1.60]1.41[6.3,6.4]1.33[7.8,8.0]0.35
[1900,1925]7.47[46.0,46.5]7.61 [6.4,6.5]0.42
[1925,1950]7.69[46.5,47.0]7.05
[1950,1975]7.82[47.0,47.5]9.11
[1975,1996]8.35[47.5,48.0]4.62
[48.0,48.5]1.73
[48.5,49.1]1.03
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yin, S.; Jin, K.; Bai, Y.; Zhou, W.; Wang, Z. Solution-Space-Reduction-Based Evidence Theory Method for Stiffness Evaluation of Air Springs with Epistemic Uncertainty. Mathematics 2023, 11, 1214. https://doi.org/10.3390/math11051214

AMA Style

Yin S, Jin K, Bai Y, Zhou W, Wang Z. Solution-Space-Reduction-Based Evidence Theory Method for Stiffness Evaluation of Air Springs with Epistemic Uncertainty. Mathematics. 2023; 11(5):1214. https://doi.org/10.3390/math11051214

Chicago/Turabian Style

Yin, Shengwen, Keliang Jin, Yu Bai, Wei Zhou, and Zhonggang Wang. 2023. "Solution-Space-Reduction-Based Evidence Theory Method for Stiffness Evaluation of Air Springs with Epistemic Uncertainty" Mathematics 11, no. 5: 1214. https://doi.org/10.3390/math11051214

APA Style

Yin, S., Jin, K., Bai, Y., Zhou, W., & Wang, Z. (2023). Solution-Space-Reduction-Based Evidence Theory Method for Stiffness Evaluation of Air Springs with Epistemic Uncertainty. Mathematics, 11(5), 1214. https://doi.org/10.3390/math11051214

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