# Global Sensitivity Analysis of Fuzzy Distribution Parameter  on Failure Probability and Its Single-Loop Estimation

Authors: Lei Cheng; Zhenzhou Lu; Luyi Li
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

An extending Borgonovo’s global sensitivity analysis is proposed to measure the influence of fuzzy distribution parameters on fuzzy failure probability by averaging the shift between the membership functions (MFs) of unconditional and conditional failure probability. The presented global sensitivity indices can reasonably reflect the influence of fuzzy-valued distribution parameters on the character of the failure probability, whereas solving the MFs of unconditional and conditional failure probability is time-consuming due to the involved multiple-loop sampling and optimization operators. To overcome the large computational cost, a single-loop simulation (SLS) is introduced to estimate the global sensitivity indices. By establishing a sampling probability density, only a set of samples of input variables are essential to evaluate the MFs of unconditional and conditional failure probability in the presented SLS method. Significance of the global sensitivity indices can be verified and demonstrated through several numerical and engineering examples.


---

## Body

## 1. Introduction

Two different uncertainty sources are involved in reliability engineering: aleatory uncertainty and epistemic uncertainty [1–3]. Aleatory uncertainty describes the inherent variability associated with a structural system, which is referred to as irreducible, objective uncertainty, while epistemic uncertainty results from lack of knowledge of fundamental phenomena and is related to our ability to understand, measure, and describe the system under study. Employing the probability theory, the aleatory uncertainty can be modeled as random variable with probability density function (PDF). When dealing with the epistemic uncertainty, more than one theory has been developed, such as the interval analysis [4–6], the evidence theory [7, 8], and the fuzzy set theory [9].Input variables can be modeled as random variables when pure aleatory uncertainty is present. Thus available sensitivity analysis (SA) techniques can be applied to study how uncertainty in the output of a model can be apportioned to different sources of uncertainty in the model input. SA is commonly divided into two groups: local SA (LSA) and global SA (GSA). LSA is defined as the partial derivatives of model output with respect to model input, which is generally computed at one given point. GSA, also known as importance analysis, focuses on determining which of the input variables affects output most in the whole uncertainty. With the ranking of the input variables resulting from importance analysis, one can pay more attention to variables with high importance and neglect variables with low importance during the process of design and optimization, thus providing useful information and guidance for engineers. At present, many GSA techniques and indices are available in the case of randomness, such as nonparametric techniques [10–12], variance-based importance measure indices [13–16], and moment-independent importance measures [17–19], among which moment independent sensitivity analysis is attracting growing attention among both academics and practitioners. The moment independent sensitivity index δ
i proposed by Borgonovo is widely extended in GSA: Castaings et al. [20], Plischke et al. [21], and Luo et al. [22] developed methods for moment independent indices recently; the moment independent model was applied to environmental models [23] and climate model [24] et al.; Wei et al. [25] and Borgonovo et al. [26] brought up new moment independent indices.When fuzzy set theory is employed, assigned membership function (MF) is used to model the epistemic uncertainty. Subsequently, methods for calculating the possibility distribution of output response are proposed [27–32]. Based on the fuzzy set theory which models the epistemic uncertainty as a fuzzy variable represented by MF, some methods have been developed for evaluating the importance of fuzzy variables [29–31, 33], among which Song et al. [33] proposed a generalized Borgonovo’s sensitivity indicator to analyze the effect of fuzzy-valued variables on the output response, and optimization techniques are employed to calculate the MF of output response [34]. When two kinds of uncertainties are concomitant in reliability problems, where aleatory uncertainty is described by PDF and epistemic uncertainty is modeled by fuzzy set theory, GSA for both fuzzy and random input variables is constructed by Zhangchun et al. in recent years [35–37].However, obviously, the importance of basic variable is affected by distribution parameters of all involved basic random variables which need to be designed for the required output response. Therefore, it is essential to analyze the importance of distribution parameter on probabilistic response in engineering applications. Recently, GSA is constructed for systems with both input variables and distribution parameters described by PDF [38–40], whereas in this paper, the input variables are described by PDF and the distribution parameters are modeled as fuzzy parameters with assigned MF [41–43].Then, a global sensitivity analysis is proposed to estimate the effect of fuzzy-valued distribution parameters on failure probability. Estimation of the established global sensitivity indices involves calculating the MFs of unconditional and conditional failure probability. For the problem that the computational cost for direct Monte-Carlo simulation (MCS) is rather too expensive to be accepted, a single-loop simulation (SLS) based on the concept of importance sampling is introduced to calculate the failure probability function (FPF), where a sampling PDF is introduced to generate a group of samples of the input variables. Thus, only one set of function calls is necessary to evaluate the sensitivity indices, and the efficiency has been greatly improved comparing with the existed direct MCS.The paper is organized as follows. In Section2, sensitivity indices of fuzzy-valued distribution parameters on failure probability are proposed and their properties are established. Section 3 gives the detailed procedures for direct MCS. Section 4 presents the SLS which can greatly save the computation cost. Applicability of sensitivity indices and accuracy and efficiency of the SLS are demonstrated by the numerical and engineering examples in Section 5. The paper ends with conclusions in Section 6.


## 2. Global Sensitivity Indices for Fuzzy Distribution Parameter Based on Failure Probability

### 2.1. Definitions of Global Sensitivity Indices

Considering structural systems with both epistemic and aleatory uncertainties, the output response can be given as(1)
Y
=
g
(
X
)
,
where Y is the output response. Let X denote the aleatory variables whose PDF is given by f
X
(
x
∣
θ
*
), and θ
=
(
θ
1
,
θ
2
,
…
,
θ
p
) refers to the fuzzy-valued distribution parameter vector, and its characteristic is represented by a family of MFs noted by μ
θ
i
(
θ
i
)
(
i
=
1,2
,
…
,
p
). The membership grade μ
θ
i
(
ω
) quantifies the degree of possibility that parameter θ
i takes value ω. It is possible to evaluate fuzzy parameters using λ-cuts or membership levels. Generally, the cut level λ takes value from the interval [
0,1
]. At each λ-cut level, the variation of the fuzzy variable is defined by a lower bound θ
i
L
(
λ
) and upper bound θ
i
U
(
λ
).Once the bounds of distribution parameters atλ-cut are obtained by the corresponding MF μ
θ
i
(
θ
i
), the bounds [
P
f
L
(
θ
(
λ
)
)
,
P
f
U
(
θ
(
λ
)
)
] of failure probability at various λ-cut can be calculated; that is, the MF μ
P
f of unconditional failure probability is obtained.While theith distribution parameter θ
i was fixed at a value θ
i
*, that is, θ
*
=
(
θ
1
,
…
,
θ
i
-
1
,
θ
i
*
,
θ
i
+
1
,
…
,
θ
p
), the MF μ
P
f
′
θ
i of conditional failure probability can also be evaluated. Obviously, the shift between μ
P
f and μ
P
f
′
θ
i measured by the area s
θ
i, which is the shadow area shown in Figure 1, can reflect the effect of the fixed parameter θ
i on failure probability:
(2)
s
θ
i
=
∫
0
1
|
μ
P
f
-
μ
P
f
′
θ
i
|
d
P
f
.Figure 1
Shift betweenμ
P
f and μ
P
f
′
θ
i.The larger the value ofs
θ
i is, the more important the fixed value θ
i for failure probability is. In average, the global sensitivity index of θ
i can be measured by
(3)
τ
θ
i
=
∫
-
∞
+
∞
μ
θ
i
s
θ
i
d
θ
i
.However, in the engineering, the fuzzy-valued distribution parameters may have different dimensions. In order to eliminate the infections of the units of distribution parameters, we normalize the distribution parameters by introducing the equivalent PDFμ
θ
i
(
e
)
=
μ
θ
i
/
∫
μ
θ
i
d
θ
i, which satisfies the property that ∫
μ
θ
i
(
e
)
d
θ
i
=
1. The equivalent PDF can eliminate the effect of the dimensions, while it does not change the distribution type of the MFs of distribution parameters [44] Thus, we have
(4)
Γ
θ
i
=
∫
-
∞
+
∞
μ
θ
i
∫
μ
θ
i
d
θ
i
s
θ
i
d
θ
i
=
∫
-
∞
+
∞
μ
θ
i
∫
μ
θ
i
d
θ
i
∫
0
1
|
μ
P
f
-
μ
P
f
′
θ
i
|
d
P
f
d
θ
i
.From the above definitions, we know thatΓ
θ
i represents the effect of the distribution parameter θ
i on failure probability of structural model.Generally, MFs of failure probability cannot be expressed by explicit expression, which are determined by the lower and upper bounds at variousλ-cuts. So the integral of (4) needs numerical integration to be evaluated; that is,
(5)
Γ
θ
i
=
∑
j
=
1
N
{
∑
j
=
1
N
1
Δ
θ
i
·
μ
θ
i
(
θ
i
,
j
)
∫
μ
θ
i
d
θ
i
h
h
h
h
h
×
∑
k
=
1
M
[
Δ
P
f
·
|
μ
P
f
(
P
f
(
k
)
)
-
μ
P
f
′
θ
i
(
P
f
(
k
)
)
|
]
μ
θ
i
(
θ
i
,
j
)
∫
μ
θ
i
d
θ
i
∑
j
=
1
N
1
}
,
where Δ
θ
i and Δ
P
f are the spans of isometric sampling, θ
i
,
j is the jth fixed value of the distribution parameter θ
i, and P
f
(
k
) is the kth value of unconditional failure probability P
f.Global sensitivity index for an individual distribution parameter can be easily extended to a group of distribution parameters. Then the sensitivity index for a group of parameters can be given as:(6)
Γ
θ
i
,
θ
j
=
∫
-
∞
+
∞
∫
-
∞
+
∞
μ
θ
i
∫
μ
θ
i
d
θ
i
μ
θ
j
∫
μ
θ
j
d
θ
j
h
h
h
h
h
h
h
h
h
×
∫
0
1
|
μ
P
f
-
μ
P
f
′
θ
i
,
θ
j
|
d
P
f
d
θ
i
d
θ
j
Γ
θ
i
∣
θ
j
=
∫
-
∞
+
∞
∫
-
∞
+
∞
μ
θ
i
∫
μ
θ
i
d
θ
i
μ
θ
j
∫
μ
θ
j
d
θ
j
h
h
h
h
h
h
h
h
h
×
∫
0
1
|
μ
P
f
′
θ
i
-
μ
P
f
′
θ
i
,
θ
j
|
d
P
f
d
θ
i
d
θ
j
,
where μ
P
f
′
θ
i
,
θ
j is the MF of conditional failure probability with distribution parameters θ
i and θ
j taken the fixed values θ
i
* and θ
j
*, respectively.Similarly, the numerical integrations for (6) are shown as follows:
(7)
Γ
θ
i
,
θ
j
=
∑
l
=
1
N
i
∑
m
=
1
N
j
Δ
θ
i
Δ
θ
j
·
μ
θ
i
(
θ
i
,
l
)
∫
μ
θ
i
d
θ
i
μ
θ
j
(
θ
j
,
m
)
∫
μ
θ
j
d
θ
j
h
h
h
h
h
h
h
×
{
∑
q
=
1
M
[
Δ
P
f
·
|
μ
P
f
(
P
f
(
k
)
)
-
μ
P
f
′
θ
i
,
θ
j
(
P
f
(
k
)
)
|
]
}
Γ
θ
i
∣
θ
j
=
∑
l
=
1
N
i
∑
m
=
1
N
j
Δ
θ
i
Δ
θ
j
·
μ
θ
i
(
θ
i
,
l
)
∫
μ
θ
i
d
θ
i
μ
θ
j
(
θ
j
,
m
)
∫
μ
θ
j
d
θ
j
h
h
h
h
h
h
h
×
{
∑
q
=
1
M
[
Δ
P
f
·
|
μ
P
f
′
θ
i
(
P
f
(
k
)
)
-
μ
P
f
′
θ
i
,
θ
j
(
P
f
(
k
)
)
|
]
}
.Mathematical properties of the sensitivity indices are discussed in Section2.2.


### 2.2. Mathematical Properties of the Proposed Sensitivity Indices

Similar to the sensitivity indices of the random input variables proposed by Song et al. [33], the mathematical properties of the proposed sensitivity indices for fuzzy distribution parameters can be explained as follows.Property 1.
Γ
θ
i
≥
0.Property 2.
If failure probabilityP
f is independent of θ
i, then Γ
θ
i
=
0.Property 3.
IfP
f is dependent of θ
i but independent of θ
j, then Γ
θ
i
θ
j
=
Γ
θ
i, Γ
θ
i
∣
θ
j
=
0, and Γ
θ
j
∣
θ
i
=
Γ
θ
i.Proofs for the above properties are similar with those in [19].Note that since the sensitivity indexδ
i proposed by Borgonovo et al. has the property of transformation invariance [45, 46], thus the algorithms of computing the delta index can converge even faster if a monotonic transformation is applied when dealing with large distribution region of the model output. From the definitions in (2) and (3), we can see that the inner statistic of Γ
θ
i and inner statistic of δ
i both are L
1-norm on densities; thus; the moment independent sensitivity index for fuzzy-valued parameters is also monotonic transformation invariant. For knowing about the detailed property of moment-independent sensitivity index, one could refer to [45, 46].


## 2.1. Definitions of Global Sensitivity Indices

Considering structural systems with both epistemic and aleatory uncertainties, the output response can be given as(1)
Y
=
g
(
X
)
,
where Y is the output response. Let X denote the aleatory variables whose PDF is given by f
X
(
x
∣
θ
*
), and θ
=
(
θ
1
,
θ
2
,
…
,
θ
p
) refers to the fuzzy-valued distribution parameter vector, and its characteristic is represented by a family of MFs noted by μ
θ
i
(
θ
i
)
(
i
=
1,2
,
…
,
p
). The membership grade μ
θ
i
(
ω
) quantifies the degree of possibility that parameter θ
i takes value ω. It is possible to evaluate fuzzy parameters using λ-cuts or membership levels. Generally, the cut level λ takes value from the interval [
0,1
]. At each λ-cut level, the variation of the fuzzy variable is defined by a lower bound θ
i
L
(
λ
) and upper bound θ
i
U
(
λ
).Once the bounds of distribution parameters atλ-cut are obtained by the corresponding MF μ
θ
i
(
θ
i
), the bounds [
P
f
L
(
θ
(
λ
)
)
,
P
f
U
(
θ
(
λ
)
)
] of failure probability at various λ-cut can be calculated; that is, the MF μ
P
f of unconditional failure probability is obtained.While theith distribution parameter θ
i was fixed at a value θ
i
*, that is, θ
*
=
(
θ
1
,
…
,
θ
i
-
1
,
θ
i
*
,
θ
i
+
1
,
…
,
θ
p
), the MF μ
P
f
′
θ
i of conditional failure probability can also be evaluated. Obviously, the shift between μ
P
f and μ
P
f
′
θ
i measured by the area s
θ
i, which is the shadow area shown in Figure 1, can reflect the effect of the fixed parameter θ
i on failure probability:
(2)
s
θ
i
=
∫
0
1
|
μ
P
f
-
μ
P
f
′
θ
i
|
d
P
f
.Figure 1
Shift betweenμ
P
f and μ
P
f
′
θ
i.The larger the value ofs
θ
i is, the more important the fixed value θ
i for failure probability is. In average, the global sensitivity index of θ
i can be measured by
(3)
τ
θ
i
=
∫
-
∞
+
∞
μ
θ
i
s
θ
i
d
θ
i
.However, in the engineering, the fuzzy-valued distribution parameters may have different dimensions. In order to eliminate the infections of the units of distribution parameters, we normalize the distribution parameters by introducing the equivalent PDFμ
θ
i
(
e
)
=
μ
θ
i
/
∫
μ
θ
i
d
θ
i, which satisfies the property that ∫
μ
θ
i
(
e
)
d
θ
i
=
1. The equivalent PDF can eliminate the effect of the dimensions, while it does not change the distribution type of the MFs of distribution parameters [44] Thus, we have
(4)
Γ
θ
i
=
∫
-
∞
+
∞
μ
θ
i
∫
μ
θ
i
d
θ
i
s
θ
i
d
θ
i
=
∫
-
∞
+
∞
μ
θ
i
∫
μ
θ
i
d
θ
i
∫
0
1
|
μ
P
f
-
μ
P
f
′
θ
i
|
d
P
f
d
θ
i
.From the above definitions, we know thatΓ
θ
i represents the effect of the distribution parameter θ
i on failure probability of structural model.Generally, MFs of failure probability cannot be expressed by explicit expression, which are determined by the lower and upper bounds at variousλ-cuts. So the integral of (4) needs numerical integration to be evaluated; that is,
(5)
Γ
θ
i
=
∑
j
=
1
N
{
∑
j
=
1
N
1
Δ
θ
i
·
μ
θ
i
(
θ
i
,
j
)
∫
μ
θ
i
d
θ
i
h
h
h
h
h
×
∑
k
=
1
M
[
Δ
P
f
·
|
μ
P
f
(
P
f
(
k
)
)
-
μ
P
f
′
θ
i
(
P
f
(
k
)
)
|
]
μ
θ
i
(
θ
i
,
j
)
∫
μ
θ
i
d
θ
i
∑
j
=
1
N
1
}
,
where Δ
θ
i and Δ
P
f are the spans of isometric sampling, θ
i
,
j is the jth fixed value of the distribution parameter θ
i, and P
f
(
k
) is the kth value of unconditional failure probability P
f.Global sensitivity index for an individual distribution parameter can be easily extended to a group of distribution parameters. Then the sensitivity index for a group of parameters can be given as:(6)
Γ
θ
i
,
θ
j
=
∫
-
∞
+
∞
∫
-
∞
+
∞
μ
θ
i
∫
μ
θ
i
d
θ
i
μ
θ
j
∫
μ
θ
j
d
θ
j
h
h
h
h
h
h
h
h
h
×
∫
0
1
|
μ
P
f
-
μ
P
f
′
θ
i
,
θ
j
|
d
P
f
d
θ
i
d
θ
j
Γ
θ
i
∣
θ
j
=
∫
-
∞
+
∞
∫
-
∞
+
∞
μ
θ
i
∫
μ
θ
i
d
θ
i
μ
θ
j
∫
μ
θ
j
d
θ
j
h
h
h
h
h
h
h
h
h
×
∫
0
1
|
μ
P
f
′
θ
i
-
μ
P
f
′
θ
i
,
θ
j
|
d
P
f
d
θ
i
d
θ
j
,
where μ
P
f
′
θ
i
,
θ
j is the MF of conditional failure probability with distribution parameters θ
i and θ
j taken the fixed values θ
i
* and θ
j
*, respectively.Similarly, the numerical integrations for (6) are shown as follows:
(7)
Γ
θ
i
,
θ
j
=
∑
l
=
1
N
i
∑
m
=
1
N
j
Δ
θ
i
Δ
θ
j
·
μ
θ
i
(
θ
i
,
l
)
∫
μ
θ
i
d
θ
i
μ
θ
j
(
θ
j
,
m
)
∫
μ
θ
j
d
θ
j
h
h
h
h
h
h
h
×
{
∑
q
=
1
M
[
Δ
P
f
·
|
μ
P
f
(
P
f
(
k
)
)
-
μ
P
f
′
θ
i
,
θ
j
(
P
f
(
k
)
)
|
]
}
Γ
θ
i
∣
θ
j
=
∑
l
=
1
N
i
∑
m
=
1
N
j
Δ
θ
i
Δ
θ
j
·
μ
θ
i
(
θ
i
,
l
)
∫
μ
θ
i
d
θ
i
μ
θ
j
(
θ
j
,
m
)
∫
μ
θ
j
d
θ
j
h
h
h
h
h
h
h
×
{
∑
q
=
1
M
[
Δ
P
f
·
|
μ
P
f
′
θ
i
(
P
f
(
k
)
)
-
μ
P
f
′
θ
i
,
θ
j
(
P
f
(
k
)
)
|
]
}
.Mathematical properties of the sensitivity indices are discussed in Section2.2.


## 2.2. Mathematical Properties of the Proposed Sensitivity Indices

Similar to the sensitivity indices of the random input variables proposed by Song et al. [33], the mathematical properties of the proposed sensitivity indices for fuzzy distribution parameters can be explained as follows.Property 1.
Γ
θ
i
≥
0.Property 2.
If failure probabilityP
f is independent of θ
i, then Γ
θ
i
=
0.Property 3.
IfP
f is dependent of θ
i but independent of θ
j, then Γ
θ
i
θ
j
=
Γ
θ
i, Γ
θ
i
∣
θ
j
=
0, and Γ
θ
j
∣
θ
i
=
Γ
θ
i.Proofs for the above properties are similar with those in [19].Note that since the sensitivity indexδ
i proposed by Borgonovo et al. has the property of transformation invariance [45, 46], thus the algorithms of computing the delta index can converge even faster if a monotonic transformation is applied when dealing with large distribution region of the model output. From the definitions in (2) and (3), we can see that the inner statistic of Γ
θ
i and inner statistic of δ
i both are L
1-norm on densities; thus; the moment independent sensitivity index for fuzzy-valued parameters is also monotonic transformation invariant. For knowing about the detailed property of moment-independent sensitivity index, one could refer to [45, 46].


## 3. MCS for Evaluating MFs of Failure Probability

Failure probabilityP
f is a function of the basic fuzzy-valued distribution parameter θ
=
(
θ
1
,
θ
2
,
…
,
θ
p
), given by P
f
=
ψ
(
θ
1
,
θ
2
,
…
,
θ
p
). Details of MCS for evaluating the MF of failure probability are as follows.(1)
Atλ
k-cut (
k
=
1,2
,
…
,
N
λ
), the lower bounds and upper bounds of the parameters θ
=
(
θ
1
,
θ
2
,
…
,
θ
p
) can be obtained as [
θ
i
L
(
λ
k
)
,
θ
i
U
(
λ
k
)
]
(
i
=
1,2
,
…
,
p
) according to the MF μ
θ
i
(
θ
i
) of basic fuzzy-valued parameter θ
i. Generate N samples θ
t
=
(
θ
1
t
,
θ
2
t
,
…
,
θ
p
t
)
(
t
=
1,2
,
…
,
N
) according to the intervals of distribution parameters at the given λ
k-cut.
(2)
For each sampleθ
t, generate N samples x
s
=
(
x
1
s
,
x
2
s
,
…
,
x
n
s
)
(
s
=
1,2
,
…
,
N
) of aleatory variables X, and then the corresponding N samples P
f
s
(
s
=
1,2
,
…
,
N
) of failure probability P
f can be evaluated. Thus, the lower bound P
f
L
(
λ
k
) and upper bound P
f
U
(
λ
k
) of unconditional failure probability at given λ
k-cut can be obtained.
(3)
Iterate (1) and (2) forN
λ times; then, the MF of unconditional failure probability μ
P
f is evaluated.
(4)
Change parameter matrixθ
t
=
(
θ
1
t
,
θ
2
t
,
…
,
θ
p
t
) to θ
t
*
=
(
θ
1
t
,
θ
2
t
,
…
,
θ
i
t
*
,
…
,
θ
p
t
) by fixing the ith distribution parameter θ
i at its fixed value θ
i
t
(
i
=
1,2
,
…
,
p
,
t
=
1,2
,
…
,
N
).
(5)
For each sampleθ
t
*, generate N samples x
s
*
=
(
x
1
s
,
x
2
s
,
…
,
x
n
s
)
(
s
=
1,2
,
…
,
N
) of aleatory variables X; then, the corresponding N samples P
f
s
*
(
s
=
1,2
,
…
,
M
) of failure probability P
f
* can be evaluated. Thus, the lower bound P
f
L
*
(
λ
k
) and upper bound P
f
U
*
(
λ
k
) of conditional failure probability at given λ
k-cut can be obtained.
(6)
Iterate (4) and (5) forN times; then, the N MFs μ
P
f
∣
θ
i
=
θ
i
t
(
t
=
1,2
,
…
,
N
) of conditional failure probability can be obtained.
(7)
Iterate (4), (5), and (6) forp times; all the MFs of conditional failure probability of distribution parameters μ
P
f
|
θ
j
=
θ
j
t
(
j
=
1,2
,
…
,
p
,
t
=
1,2
,
…
,
N
) can be evaluated.We can see from the procedures in Figure2 that using MCS combing optimization techniques N
×
N
×
N
λ and p
×
N
×
N
×
N
×
N
λ function calls are needed to calculate the MFs of unconditional and conditional failure probability, which is unacceptable especially in engineering problems. Thus, a single-loop method is proposed to evaluate the FPF of the distribution parameters in Section 4. (N is the sample number of distribution parameters and N
λ is the cut number of the membership functions.)Flow chart of MCS for evaluating the MFs of unconditional and conditional failure probability.
(a)
MCS for evaluating MF of unconditional failure probability
(b)
MCS for evaluating MFs of conditional failure probability


## 4. SLS for Evaluating MFs of Failure Probability

By definition, FPF is expressed as(8)
P
f
(
θ
)
=
∫
g
(
x
)
≤
0
I
F
(
x
)
f
(
x
θ
)
d
x
.Suppose now we have a joint PDFh
x
(
x
) of model input vector with constant distribution parameter vector θ
*
=
(
θ
1
*
,
θ
2
*
,
…
,
θ
m
*
). Then the FPF can be derived as
(9)
P
f
(
θ
)
=
∫
I
F
(
x
)
f
x
(
x
∣
θ
)
h
x
(
x
)
h
x
(
x
)
d
x
=
E
h
(
I
F
(
x
)
f
x
(
x
∣
θ
)
h
x
(
x
)
)
,
where h
x
(
x
) is an arbitrary selected joint PDF of X with constant distribution parameters θ
*. Equation (9) indicates that P
f
(
θ
) can be expressed as the expectation of the function I
F
(
x
)
f
x
(
x
∣
θ
)
/
h
x
(
x
). Given a set of model input vector samples x
(
k
)
=
(
x
1
(
k
)
,
x
2
(
k
)
,
…
,
x
n
(
k
)
)
(
k
=
1,2
,
…
,
N
) generated by the joint PDF h
x
(
x
), then P
f
(
θ
) can be estimated by
(10)
P
^
f
(
θ
)
=
1
N
∑
k
=
1
N
I
F
(
x
(
k
)
)
f
x
(
x
(
k
)
θ
)
h
x
(
x
(
k
)
)
.Detailed procedures for SLS are as follows.(1)
GenerateN samples x
s
=
(
x
1
s
,
x
2
s
,
…
,
x
n
s
)
(
s
=
1,2
,
…
,
N
) according to the joint PDF h
x
(
x
). For each x
s
(
s
=
1,2
,
…
,
N
), calculate the values of selected joint PDF h
x
(
x
s
) and failure indicator I
F
(
x
s
) as h
s and I
F
s.
(2)
Atλ
k-cut (
k
=
1,2
,
…
,
N
λ
), the lower and upper bounds of the parameters θ
=
(
θ
1
,
θ
2
,
…
,
θ
p
) can be obtained as [
θ
i
L
(
λ
k
)
,
θ
i
U
(
λ
k
)
]
(
i
=
1,2
,
…
,
p
) according to the MF μ
θ
i
(
θ
i
) of basic fuzzy-valued parameter. Generate N samples θ
t
=
(
θ
1
t
,
θ
2
t
,
…
,
θ
p
t
)
(
t
=
1,2
,
…
,
N
) according to the intervals of distribution parameters at the given λ
k-cut.
(3)
For eachθ
t
(
t
=
1,2
,
…
,
N
), calculate the values of the original joint PDF f
x
(
x
s
∣
θ
) as f
s. Then (10) can be calculated as (
1
/
N
)
∑
s
=
1
N
I
F
s
f
s
/
h
s; that is, P
^
f
(
θ
t
) is obtained.
(4)
Iterate (3) forN times; then, P
^
f
(
θ
) is calculated. Thus, the lower bound P
^
f
L and upper bound P
^
f
U of unconditional failure probability at given λ
k-cut are obtained.
(5)
Iterate (3) and (4) forN
λ times, and the MF μ
P
f of unconditional failure probability can be evaluated.
(6)
Change parameter matrixθ
t
=
(
θ
1
t
,
θ
2
t
,
…
,
θ
p
t
) to θ
t
*
=
(
θ
1
t
,
θ
2
t
,
…
,
θ
i
t
*
,
…
,
θ
p
t
) by fixing the ith distribution parameter θ
i at its fixed value θ
i
t
(
i
=
1,2
,
…
,
p
,
t
=
1,2
,
…
,
N
).
(7)
For eachθ
t
*
(
t
=
1,2
,
…
,
N
), calculate the values of the original joint PDF f
x
(
x
s
∣
θ
) as f
s
*. Then (10) can be calculated as (
1
/
N
)
∑
k
=
1
N
I
F
s
f
s
*
/
h
s; that is, P
^
f
(
θ
t
*
) is obtained.
(8)
Iterate (7) forN times; then, P
^
f
(
θ
*
) is calculated. Thus, the lower bound P
^
f
L
* and upper bound P
^
f
U
* of conditional failure probability at given λ
k-cut are obtained.
(9)
Iterate (7) and (8) forN
λ times, and the MF μ
P
f
∣
θ
i
=
θ
i
t of conditional failure probability can be evaluated.
(10)
Iterate (7), (8), and (9) forN
×
p times, and all the MFs μ
P
f
∣
θ
i
=
θ
i
t
(
i
=
1,2
,
…
,
p
,
t
=
1,2
,
…
,
N
) of conditional failure probability are evaluated.The procedure in Figure3 shows that given any λ
k-cut (
k
=
1,2
,
…
,
N
λ
), the interval of distribution parameters [
θ
L
(
λ
k
)
,
θ
U
(
λ
k
)
] can be obtained; thus, the failure probability can be estimated using (10); that is, (10) can be treated as the imitate explicit function of the distribution parameters. Then any arbitrary optimization method can be used to get the upper and lower bounds of failure probability [
P
f
L
(
λ
k
)
,
P
f
U
(
λ
k
)
] at given interval [
θ
L
(
λ
k
)
,
θ
U
(
λ
k
)
] of distribution parameters. Only one set of samples of the aleatory variables is used in the procedure at different λ
k-cut and only N function calls are needed to calculate the MF of unconditional failure probability. The procedures for calculating the MFs of conditional failure probability are similar, and the only difference is that the one or a group of distribution parameters are fixed at certain value.Flow chart of SLS for evaluating the MFs of unconditional and conditional failure probability.
(a)
SLS for evaluating MF of unconditional failure probability
(b)
SLS for evaluating MFs of conditional failure probabilityHence the above procedures for estimating the FPF need only one set of samples for implementing it; that is, the function calls for SLS isN; thus, it is computationally efficient.Considering that the computational cost is rather expensive, optimization techniques can be employed to calculate the minimum and maximum values of the failure probability within the specified bounds of the cuts of the fuzzy-valued parameters corresponding toλ-cut. Thus, combing MCS and SLS with local/global optimization techniques can be efficient. Accuracy and efficiency of the SLS can be established in Section 5.


## 5. Examples

A number of examples are used to demonstrate the availability and feasibility of the presented sensitivity indicator. These examples involve some cases met in engineering commonly, for example, normal, triangular, or trapezoidal types MFs, and MFs of failure probability are calculated by employing the optimization techniques.Example 1.
The response functionz
=
4
-
x
1
-
x
2 has two variables. Case 1: x
i
~
N
(
μ
i
,
2
), μ
i
(
i
=
1,2
) are fuzzy-valued distribution parameters. Case 2: x
i
~
N
(
1
,
σ
i
), σ
i
(
i
=
1,2
) are fuzzy-valued parameters. The MFs are listed as follows:
(11)
μ
μ
i
(
μ
i
)
=
{
4
-
μ
i
3
≤
μ
i
≤
4
μ
i
-
2
2
≤
μ
i
≤
3
,
μ
σ
i
(
σ
i
)
=
{
4
-
σ
i
3
≤
σ
i
≤
4
σ
i
-
2
2
≤
σ
i
≤
3
.The MFs of failure probability are drawn in Figures4 and 5, and the results of global sensitivity indices calculated by MCS and SLS combing with optimization techniques are listed in Table 1.Table 1
Results of global sensitivity indices for Example1, Case 1.
MCS
SLS
{
Γ
μ
1
,
Γ
μ
2
}
{0.0461, 0.0461}
{0.0458, 0.0463}
{
Γ
μ
1
μ
2
}
{0.0356}
{0.0362}
{
Γ
μ
1
∣
μ
2
,
Γ
μ
2
∣
μ
1
}
{0.0192, 0.0193}
{0.0187, 0.0189}
Function Calls
(
N
=
2000
,
N
λ
=
21
,
p
=
2
)
p
×
N
×
N
×
N
×
N
λ
+
N
×
N
×
N
λ
NFigure 4
MFs of failure probability for Example1, Case 1.Figure 5
MFs of failure probability for Example1, Case 2.Case 1.
See Figure4 and Table 1.Case 2.
See Figure5 and Table 2.Table 2
Results of global sensitivity indices for Example1, Case 2.
MCS
SLS
{
Γ
σ
1
,
Γ
σ
2
}
{0.0172, 0.0178}
{0.0175, 0.0173}
{
Γ
σ
1
σ
2
}
{0.0281}
{0.0279}
{
Γ
σ
1
∣
σ
2
,
Γ
σ
2
∣
σ
1
}
{0.0149, 0.0147}
{0.0151, 0.0146}Figures4 and 5 show that the MFs of failure probability obtained by the SLS are in good accordance with those obtained by MCS, which testifies the accuracy of SLS. And Tables 1 and 2 indicate that the sensitivity indices of the two distribution parameters are equal, which is reasonable in this example. In addition, the accuracy of SLS is proved again by the results estimated by the two methods in Tables 1 and 2. In addition, the function calls of the two methods are also listed in Table 1, proving that SLS is greatly efficient in evaluating the proposed sensitivity indices.To make the numerical results more robust, we carry out the calculations of sensitivity indices at increasing sample sizes to demonstrate convergence. Figures6 and 7 show the convergence plot for Case 1 using MCS and SLS, respectively. One can easily find that the results begin to converge at the level of 2
×
1
0
3. Thus we take N
=
2
×
1
0
3 as shown in Table 1.Figure 6
Convergence plot of MCS for Example1, Case 1.Figure 7
Convergence plot of SLS for Example1, Case 1.Example 2.
The response functiong
=
sin
x
1
+
7
si
n
2
x
2
+
0.1
x
3
4
sin
x
1 includes three variables, where x
i
~
N
(
0
,
σ
i
). And the MFs of distribution parameters σ
i are given as
(12)
μ
σ
i
(
σ
i
)
=
{
0
σ
i
≤
0.3
π
5
π
σ
i
-
1.5
0.3
π
≤
σ
i
≤
0.5
π
1
0.5
π
≤
σ
i
≤
0.7
π
8
-
10
σ
i
π
0.7
π
≤
σ
i
≤
0.8
π
0
σ
i
≥
0.8
π
.The MFs of failure probability are drawn in Figure8, and the results of sensitivity indices calculated by MCS and SLS combing with optimization techniques are listed in Table 3.Table 3
Results of global sensitivity indices for Example2.
MC
SL
{
Γ
σ
1
,
Γ
σ
2
,
Γ
σ
3
}
{0.0230, 0.0207, 0.0278}
{0.0228, 0.0201, 0.0291}
{
Γ
σ
1
σ
2
,
Γ
σ
1
σ
3
,
Γ
σ
2
σ
3
}
{0.0213, 0.0318, 0.0308}
{0.0221, 0.0336, 0.0324}
{
Γ
σ
1
∣
σ
2
,
Γ
σ
1
∣
σ
3
}
{0.0191, 0.0342}
{0.0195, 0.0361}
{
Γ
σ
2
∣
σ
1
,
Γ
σ
2
∣
σ
3
}
{0.0169, 0.0294}
{0.0187, 0.0302}
{
Γ
σ
3
∣
σ
1
,
Γ
σ
3
∣
σ
2
}
{0.0371, 0.0376}
{0.0369, 0.0381}Figure 8
MFs of failure probability for Example2.Results of sensitivity indices in Table3 indicate that the fuzzy-valued parameter σ
3 makes more contribution to failure probability of the structural model than the other parameters. And the results evaluated by SLS are in good agreement with those evaluated by MCS.Example 3.
A roof truss is shown in Figure9, and the top boom and the compression bars are reinforced by concrete, and the bottom boom and the tension bars are steel. Assume the uniformly distributed load q is applied on the roof truss, and the uniformly distributed load can be transformed into the nodal load P
=
q
l
/
4. Taking the safety and applicability into account, the perpendicular deflection Δ
C of the peak of structure node C not exceeding 2.8 cm is taken as the constraint condition, and the performance response function can be constructed by g
(
x
)
=
0.028
-
Δ
C, where Δ
C is the function of the basic random variables and Δ
C
=
(
q
l
2
/
2
)
(
3.81
/
A
C
E
C
+
1.13
/
A
S
E
S
), and A
C, A
S, E
C, E
S, and l, respectively, are sectional area, elastic modulus, length of the concrete, and steel bars; the distribution parameters of these independent normal basic random variables are listed in Table 4. It is assumed that the standard deviation of the variables is fuzzy parameters with membership function μ
σ
i
(
σ
i
)
=
exp
⁡
(
-
(
x
-
a
)
2
/
2
b
2
). And the coefficients are listed in Table 5.Table 4
Distribution parameters of the input variables of roof truss.
Random variables
q (N/m)
l (m)
A
S (m2)
A
C (m2)
E
S (N/m2)
E
C (N/m2)
Meanμ
x
20000
12
9.82
×
10
-
4
0.04
1
×
10
11
2
×
10
10
Standard deviationσ
x
σ
q
σ
l
σ
A
s
σ
A
c
σ
E
s
σ
E
cTable 5
Coefficients for membership function ofσ
x.
Standard deviationσ
x
σ
q
σ
l
σ
A
s
σ
A
c
σ
E
s
σ
E
c
Coefficienta
1400
0.12
5.892
×
10
-
5
0.0048
6
×
10
9
1.2
×
10
9
Coefficientb
700
0.006
2.95
×
10
-
6
0.0002
3
×
10
8
6
×
10
7Schematic diagram of a roof truss.
(a)
(b)The MFs of failure probability are drawn in Figure10, and the results of sensitivity indices calculated by MCS and SLS combing with optimization techniques are listed in Table 6.Table 6
Results of global sensitivity indices for Example3.
MCS
SLS
{
Γ
σ
q
,
Γ
σ
l
,
Γ
σ
A
s
,
Γ
σ
A
c
,
Γ
σ
E
s
,
Γ
σ
E
c
}
10
-
3
×
{
0.2553,0.4535,0.8478,0.7000,0.7010,0.4694
}
10
-
3
×
{
0.2385,0.4412,0.8324,0.7132,0.6829,0.4367
}
{
Γ
σ
q
σ
l
,
Γ
σ
q
σ
A
s
,
Γ
σ
q
σ
A
c
,
Γ
σ
q
σ
E
s
,
Γ
σ
q
σ
E
c
}
{0.0009, 0.0010, 0.0011, 0.0009, 0.0009}
{0.0007, 0.0008, 0.0009, 0.0007, 0.0007}
{
Γ
σ
l
σ
A
s
,
Γ
σ
l
σ
A
c
,
Γ
σ
l
σ
E
s
,
Γ
σ
l
σ
E
c
}
{0.0027, 0.0032, 0.0025, 0.0024}
{0.0023, 0.0028, 0.0022, 0.0021}
{
Γ
σ
A
s
σ
A
c
,
Γ
σ
A
s
σ
E
s
,
Γ
σ
A
s
σ
E
c
}
{0.0035, 0.0029, 0.0027}
{0.0028, 0.0025, 0.0023}
{
Γ
σ
A
c
σ
E
s
,
Γ
σ
A
c
σ
E
c
}
{0.0033, 0.0032}
{0.0027, 0.0026}
{
Γ
σ
E
s
σ
E
c
}
{0.0024}
{0.0021}
{
Γ
σ
q
∣
σ
l
,
Γ
σ
q
∣
σ
A
s
,
Γ
σ
q
∣
σ
A
c
,
Γ
σ
q
∣
σ
E
s
,
Γ
σ
q
∣
σ
E
c
}
{0.0024, 0.0023, 0.0029, 0.0024, 0.0024}
{0.0021, 0.0020, 0.0027, 0.0021, 0.0021}
{
Γ
σ
l
∣
σ
q
,
Γ
σ
l
∣
σ
A
s
,
Γ
σ
l
∣
σ
A
c
,
Γ
σ
l
∣
σ
E
s
,
Γ
σ
l
∣
σ
E
c
}
{0.0023, 0.0070, 0.0084, 0.0069, 0.0069}
{0.0020, 0.0064, 0.0078, 0.0064, 0.0064}
{
Γ
σ
A
s
∣
σ
q
,
Γ
σ
A
s
∣
σ
l
,
Γ
σ
A
s
∣
σ
A
c
,
Γ
σ
A
s
∣
σ
E
s
,
Γ
σ
A
s
∣
σ
E
c
}
{0.0023, 0.0072, 0.0087, 0.0073, 0.0072}
{0.0020, 0.0066, 0.0080, 0.0067, 0.0067}
{
Γ
σ
A
c
∣
σ
q
,
Γ
σ
A
c
∣
σ
l
,
Γ
σ
A
c
∣
σ
A
s
,
Γ
σ
A
c
∣
σ
E
s
,
Γ
σ
A
c
∣
σ
E
c
}
{0.0028, 0.0085, 0.0085, 0.0085, 0.0085}
{0.0024, 0.0080, 0.0080, 0.0080, 0.0080}
{
Γ
σ
E
s
∣
σ
q
,
Γ
σ
E
s
∣
σ
l
,
Γ
σ
E
s
∣
σ
A
s
,
Γ
σ
E
s
∣
σ
A
c
,
Γ
σ
E
s
∣
σ
E
c
}
{0.0024, 0.0072, 0.0073, 0.0087, 0.0070}
{0.0021, 0.0066, 0.0066, 0.0080, 0.0064}
{
Γ
σ
E
c
∣
σ
q
,
Γ
σ
E
c
∣
σ
l
,
Γ
σ
E
c
∣
σ
A
s
,
Γ
σ
E
c
∣
σ
A
c
,
Γ
σ
E
c
∣
σ
E
s
}
{0.0023, 0.0070, 0.0071, 0.0085, 0.0070}
{0.0020, 0.0064, 0.0064, 0.0067, 0.0064}Figure 10
MFs of failure probability for Example3.Results of sensitivity indices in Table6 show that the rank of fuzzy valued distribution parameter is σ
A
s
>
σ
E
s
>
σ
A
c
>
σ
E
c
>
σ
l
>
σ
q, among which σ
q has little effect on the failure probability of the fuzzy structural model, and the effect of σ
q with other parameters is the smallest among all the Γ
θ
i
θ
j. In addition, the sensitivity indices can be estimated by the SLS guaranteeing acceptable accuracy with only one set of samples, which can greatly improve the computational efficiency. As testified by the results of Examples 1–3, properties of the sensitivity indices are obviously satisfied.


## 6. Conclusion

A global sensitivity index of fuzzy-valued distribution parameter on the fuzzy failure probability is defined in this work. The key to calculate the sensitivity index is to evaluate the MFs of failure probability. For the problem that the computational cost of MCS is rather expensive, a SLS is introduced to estimate the proposed sensitivity indices, which can greatly improve the computational efficiency considerably with acceptable accuracy. Applications of the sensitivity indices proposed in this paper to the numerical and engineering examples show that fuzzy sensitivity indices can provide useful information about the sequence of the distribution parameters on the fuzzy uncertainty of failure probability, and the proposed sensitivity indices can identify the relative important fuzzy distribution parameters of the structural model, thus providing useful information for engineers in the process of reliability design and optimization.



---

*Source: 490718-2014-07-15.xml*