Parametric Sensitivity Analysis for Importance Measure on Failure Probability and Its Efficient Kriging Solution

Themoment-independent importance measure (IM) on the failure probability is important in system reliability engineering, and it is always influenced by the distribution parameters of inputs. For the purpose of identifying the influential distribution parameters, the parametric sensitivity of IM on the failure probability based on local and global sensitivity analysis technology is proposed. Then the definitions of the parametric sensitivities of IM on the failure probability are given, and their computational formulae are derived.The parametric sensitivity finds out how the IM can be changed by varying the distribution parameters, which provides an important reference to improve or modify the reliability properties. When the sensitivity indicator is larger, the basic distribution parameter becomes more important to the IM. Meanwhile, for the issue that the computational effort of the IM and its parametric sensitivity is usually too expensive, an active learning Kriging (ALK) solution is established in this study. Two numerical examples and two engineering examples are examined to demonstrate the significance of the proposed parametric sensitivity index, as well as the efficiency and precision of the calculation method.


Introduction
Uncertainties existing in engineering analysis and design are inherently unavoidable in nature associated with the manufacturing error, material property, loads, and so forth.Fortunately, reliability analysis and sensitivity analysis are now available to deal with the uncertainty existing in design variables to improve the performance of a mechanical or structural system [1][2][3].Reliability analysis aims at predicting the failure probability (or reliability) of the structure under the effects of random uncertainties.On the other hand, sensitivity analysis focuses on the contribution of each uncertainty or distribution parameters of the input variables [4][5][6][7].It is reasonable and practicable to obtain reliability sensitivity analysis for quantifying and ranking the effects of random uncertainties on the failure probability.In this paper, the reliability sensitivity analysis is concerned.
Sensitivity analysis (SA) is the study of how the output response of a model (numerical or otherwise) is affected by the input uncertainty, which can be classified into two groups: local SA and global SA [8].The local SA investigates how small variation of the distribution parameters near a reference point changes the output value.The classical local SA is defined as the partial derivative of the output with respect to the distribution parameters of inputs [9].The global SA, also named as importance measure (IM), gives consideration to measure the effect of the output uncertainty on the uncertainty of the input parameters, covering their variation range space as opposed to local SA using partial derivatives.
Saltelli and Marivoet [10] and Helton and Davis [11] proposed the nonparametric techniques (input-output correlation), but this method lacks model independence.With the advantage of "global, quantitative and model free, " the variance-based importance measures are gaining the increasing attention of practitioners and have been used extensively for quantitative analysis [12][13][14][15][16].However, Borgonovo addressed the following fact: "the premise of variance-based GSA technique that the variance is sufficient to identify the variability of model output is not always true" [17].The "moment-independent" importance measures have been presented [6,[16][17][18][19].They are also global, quantitative, 2 Mathematical Problems in Engineering model-free, and additionally moment-independent, thus attracting more and more attention of practitioners recently.Generally speaking, in reliability analysis, researchers often pay the most attention to the failure probability.With this respect, Cui et al. [6] introduced a moment-independent importance measure of the basic variable on the failure probability which was further developed by Li et al. [20].
This moment-independent importance measure on failure probability is applied to quantify the average effect of the basic variables on the reliability of the model and obtain the importance ranking.The IM on the failure probability can be used in the priordesign stage for variables screening when a reliability design solution is yet identified and the postdesign stage for uncertainty reduction after an optimal design has been determined.Uncertain inputs inherent in most engineering problems are assumed as random variables obeying probabilistic distributions.Obviously, system reliability and reliability IM on failure probability are decided by distribution parameters.One can directly change the input's IMs by controlling or modifying some input's distribution parameters; namely, changing the input's distribution parameters can also influence the failure probability, which would facilitate its use under various scenarios of design under uncertainty, for instance, in reliability-based design.It is necessary to further recognize effects of the distribution parameters within system reliability on the importance ranking.At present, Cui et al. [7] defined the parametric sensitivities to illustrate the influences of the distribution parameters on the importance measures.
Combined with the local SA technique of input parameters, the effects of the distribution parameters on the IM on failure probability can be introduced, by which IMs of the inputs can be controlled or modified by changing the distribution parameters.This can provide important guidance for robust design, reliability-based design, and reliability-based optimization in engineering.However, its solution still relies on the corresponding method for failure probability and the computation of the derivative operation on failure probability existing in the parametric sensitivity of IM.
The Monte Carlo simulation (MCS) procedure is easy to implement and is available for computing the parametric sensitivity of IM based purely on model evaluation [15,21], but it has to face the problem of "curse of computational cost" for the problem with small failure probability (10 −3 -10 −4 or smaller).Thus, to deal with this problem, the Kriging approach is widely used for deterministic optimization problems [22] and reliability analysis [23] has been intensively investigated.Furthermore, an advanced Kriging method, named active learning Kriging (ALK), has been proved to be highly efficient in reliability analysis problems [24][25][26].This work would employ the ALK method to compute the parametric sensitivity of IM on failure probability.In the ALK method, the Kriging model is updated by adding new training points to the design of experiment (DOE) in iterations by active learning until the Kriging model satisfies necessary accuracy.The computational efficiency of the Kriging method can be validated by several numerical and engineering examples.
The remainder of this work is organized as follows.Section 2 reviews the definition of the moment-independent importance measure of the basic variable on the failure probability.And the parametric sensitivity of IM on failure probability is firstly presented.In Section 3, the established ALK solution can effectively solve the problem that the computational cost of the parametric sensitivity of IM relies on small failure probability.The effectiveness of the proposed parametric sensitivity of IMs and efficiency of the ALK method are demonstrated by several examples in Section 4. The discussions and conclusions are given at the end of this paper.When the th input variable   is fixed at one given value, the conditional failure probability   |X  can be obtained.Based on the idea of the moment-independent importance analysis, the importance measure of basic variable X  on the failure probability is defined by Cui et al. [6] as

Definition of the Parametric Sensitivity of IM on Failure Probability
where X  represents a random basic variable   or a set of random basic variables is the operator of expectation.
As the absolute value in (1) is difficult to compute, it is transformed into square operation by Li et al. [20].The modified version of importance measures on the failure probability can be expressed as follows: In the reliability analysis, the failure domain of this structure system is defined as  = (x :  (x) < 0) . ( Suppose the indicator function of this failure domain is given as   (x); that is, Then, the unconditional failure probability and conditional failure probability on   can be expressed as Here,  X (x |   ) is the conditional joint PDF on   .According to the probability theory,  X (x |   ) is defined as where    (  ) is the PDF of   .To analyze the effect of changing the th parameter     ∈   −  of input on   X  , the derivative of the IM can be defined as

The
To compute the above formula, the derivatives of    and   |X  contribution with respect to the distribution parameter of input can be given by

Mathematical Problems in Engineering
Here, the effects of the distribution parameters on the IMs of input variables can be known, combining local and global sensitivity analysis technology.One can directly optimize the inputs' IMs by controlling or changing some inputs' distribution parameters.It can provide useful information for robust design, reliability-based design, and reliability-based optimization.
In computing the parametric sensitivity of IM, the derivative of the    and   |X  with respect to the distribution parameters must be calculated.The performance function expression is characterized by complex implicit limit state functions in most engineering problems, so the derivative of the    and   |X  is without a definite analytic expression.The Monte Carlo simulation (MCS) method can be used for computing the parametric sensitivity of IM.This method has the advantages of high-accuracy and usability based on model evaluation.However, for problems with small failure probability (10 −3 -10 −4 or smaller), only a very small portion of samples will drop into the failure domain.The computational cost is unacceptable for failure probability and the derivative of the    and   |X . For the latter problem, an active learning Kriging method combining MCS is proposed in the next section.

The Active Learning Kriging Method for the Parametric Sensitivity of IM on Failure Probability
Section 2 tells us that the key point in estimating the global reliability sensitivity indices is to provide the failure probability    , the unconditional failure probability   |X  , and the derivative of    and   |X  to the distribution parameters.In practice, there are many engineering problems characterized by complex implicit performance functions, so the computational cost of double-loop MCS method becomes larger especially for finite element (FE) model.To solve this problem, the solution of ALK method is employed to approximate the implicit performance function and improve the efficiency greatly using few training points.Generally speaking, Kriging metamodel is combined by a linear regression component and a stochastic process.It is an interpolation technique based on the statistical theory and has been used to construct the input and output (I/O) systems.The Kriging models are performed with the toolbox DACE, which is a MATLAB toolbox developed by Lophaven et al. [27].

Kriging Metamodel. The present Kriging model expresses the unknown function 𝐺(x) as [22]
where (x, ) is the deterministic part which gives an approximation of the response in mean.The second part () is the realization of stochastic process, which provides the approximation of the local fluctuation so that the whole model interpolates exactly.This part shows the following statistical characteristics: Var [ ()] =  2 , ( The mean of this component is zero, the variance is  2 , and ( 14) defines its covariance between two points in space.In (15),   and   denote two arbitrary points; (  ,   ) is the correlation function of   and   .
Several types of correlation function (  ,   ) can be used, and, in this study, the Gaussian function can be formulated as where  dv is the dimension of random input variables,    and    denote the th component of training samples   and   , respectively, and   is the correlation parameter to ensure that the metamodel is flexible enough to approximate the true function, which is usually obtained by an optimization process.
A set of   experimental samples is denoted by x = [x 1 , x 2 , . . ., x   ]  , where x  is the th training point, and responses Y = [ 1 ,  2 , . . .,    ]  is the corresponding response to x.The predicted value   (x) and predicted variance  2  (x) of the objective function (x) at an unknown point x are In (17),  denotes a  ×  correlation matrix by  , = (  ,   );   (x) = [(x, x 1 ), (x, x 2 ), . . ., (x, x   )]  is the correlative relations between unknown point x and sample points x = [x 1 , x 2 , . . ., x   ]  ; the unknown parameters β can be expressed as The Kriging variance  2  (x) is defined as the minimum square error between true response (x) and estimated value Ĝ(x) and is expressed as follows: where  =    −1 (x) − ; the unknown parameters  2  can be expressed as However, at a given unknown point X, the predicted value   (X) is not the true value of (X) and there exist some uncertainties. 2  (x) demonstrates uncertainty of the predictor value; it also provides an important index to adjudge the fitting accuracy and enables quantifying the uncertainty of predictions with an easy approach.

The ALK Method and the Solution for the Parametric Sensitivity of IM 𝛿 𝑁𝑃
.The ALK method has been applied to different fields in engineering, like efficient global optimization (EGO) [22], probabilistic analysis (PRA) [24][25][26], and reliability-based design optimization [28].In the ALK method, the learning function plays an important role in constructing the active learning Kriging model [26].A new Kriging model is updated by adding a new training point to the design of experiment (DOE) in subsequent iterations by active learning until the Kriging model satisfies necessary accuracy.According to the indicator function of this failure domain in (5), we only need to focus on the sign of the performance function in the reliability analysis.
For the uncertainty of predicted value   (X), there may exist some risk that the value of (X) is positive ((X) > 0) even if predicted value is negative (  (X) < 0).So, the points owing a high potential risk to cross the predicted separator (X) = 0 have to be added to the DOE and evaluated by the real performance function.These potentially "dangerous" points determining the precision of the failure probability are in the region: close to the limit state, have high Kriging variance or both.To identify them, a new learning function named the expected risk function (ERF) is proposed in our recent work [29].
First, if   (X) < 0, we define an indicator to measure such potential risk as As shown in Figure 2(a), (X) measures the shift of value of (X) to be positive and the larger the value of (X) is, the more risky the sign of (X) is to be changed from negative to positive.
Covering the range ((X) > 0) to calculate the mean of (X), the ERF for the case   (X) < 0 is obtained as follows: where   and   are the predicted value and variance in (17) and (19) and Φ(⋅) and (⋅) are the cumulative distribution function (CDF) and PDF of the standard normal distribution.
If   (X) > 0, as shown in Figure 2(b), we define an indicator to measure such potential risk as And the ERF for the case   (X) > 0 is obtained as follows: ) .
( The ERF is employed to measure the potential possibility that the sign of the limit state function in a point X is changed from positive to negative (or negative to positive) in the Kriging model.The point maximizing ERF should be added to the initial set of training points in DOE.
The ALK-based solution for the parametric sensitivity can be simply divided into five steps and provided as follows.
Step 1.For a model  = (x), generate  samples x  = ( 1 ,  2 ⋅ ⋅ ⋅   ) ( = 1, 2 ⋅ ⋅ ⋅ ) of the input variables by PDF  X (x) using Sobol's low-discrepancy samples.One can refer to Sobol [15] for more details.A  ×  dimension matrix X is formed and shown as where the th column of matrix X represents the generated random realizations for input variable   . Step At the same time,   |  /    can be computed combining the reliability sensitivity analysis method by (11).
Step 5.The IM     can be calculated according to the definition of ( 2 It can be seen that a large number of samples must be taken for providing precise estimates when calculating IM   and its parametric sensitivity with Monte Carlo method.In the second step of ALK procedure, one only needs to compute the limit state function values using the metamodel instead of original model; thus, this procedure can further improve the efficiency without loss of precision, especially for the finite element model (FEM).
In the next section, we introduce two numerical examples and two engineering examples for demonstrating the efficiency and precision of the calculation procedure and illustrating the engineering significance of the IM   and its parametric sensitivity.[25], which behaves nonlinearly around the limit state function.The nonlinear performance function is given as where the basic variables  1 and  2 are independent and follow normal distributions  1 ∼ (1.5, 1) and  2 ∼ (2.5, 1), respectively.
The estimates of the importance measure on the failure probability     computed by the MCS and ALK procedures and the error in each estimate are reported in Table 1.The numbers in the last column are the evaluation numbers of the performance function when calculating     .The results of their parametric sensitivities are computed using the same active learning Kriging model and listed in Table 2.
As revealed by Tables 1 and 2, the results computed by the ALK procedure are precise enough compared with those computed by the MCS procedure.The results of MCS procedure can be seen as the accurate results.Compared with the tremendously large computational cost (2 × 10 4 × 10 4 samples) of Monte Carlo method, the ALK method only begins with 12 training points, while 36 points whose sign of response has the largest potential risk to be wrongly predicted are added into the initial DOE.Additionally, Figure 3 gives iteration history of the true values of the 36 training points for Example 1 by the ALK method.The sign of the candidate points (the samples in MCS procedure) is estimated using the ALK model.Compared with the true sign of the candidate points using originality performance function, there are only 8 candidate points whose signs of responses are wrongly predicted.This can be linked to Figure 4 in which it shows the sign of the response at each Monte Carlo sampling point.It is seen that most of the added points are located in the vicinity of the limit state function.Consequently, this method certifies a correct approximation of the response using a minimum number of calls to the performance function.
From the results in Table 1, it is noted that there is no difference between the rankings of the basic variables on failure probability by the MCS and ALK procedures, namely  2 >  1 , which illustrates that the ALK method is meaningful and reasonable.Therefore,  2 should be paid more attention  to in the reliability analysis on failure probability.Results shown in Table 2 point out the influences of varying some variables' distribution parameters to the IMs on the failure probability.The uncertainty of IMs     on failure probability can be improved through modifying or controlling the distribution parameters of the variables indirectly.
Example 2 (Ishigami Function).Ishigami function [30] is a commonly used example in importance measure analysis [6,20].It can be written as follows: where the basic variables  1 ,  2 , and  3 are independent and uniformly distributed in [   ,    ] ( = 1, 2, 3),    = −, and    = .The values of constants are set as  = 5 and  = 0.1.The estimates of the importance measure on failure probability     and the results of their parametric sensitivities are listed in Tables 3 and 4, respectively.
As revealed by Tables 3 and 4, the computational cost of Monte Carlo method (3 × 10 8 samples) is tremendously large, which illustrates that our model is computationally challenging.The ALK method begins with 12 training points, while 44 points are added into the initial DOE to satisfy the accuracy.We plot the tendency of the true values of  the training points for Ishigami function computed by ALK method Figure 5. validates the high efficiency of the ALK method.The importance ranking of basic variables on failure probability obtained by     is  1 ,  2 , and  3 .Therefore,  1 should be paid more attention to in the reliability analysis on failure probability.Results shown in Table 4 point out the influences of varying some variables' distribution parameters to the the failure probability.
From Examples 1 and 2, can be found that the tance of a variable to failure probability is not just affected by the distribution of this variable but affected by those of the rest of input variables, and sometimes the latter may even play a more important role than the former.Surely, the engineering examples in Sections 4.2 and 4.3 are also following the fact.Accordingly, it is necessary to   find out all the factors that may affect the importance measures according to the failure probability, and parametric sensitivities on IMs provide a way to solve this problem.Except for that, Table 4 also shows an important property as all the variables of the model are uniform variables; namely, the numerical values of two parametric sensitivities of the IM are equal, and their signs are opposite.This property result of the parametric sensitivity of the IM is decided by ( 8) and (9).For the uniform variable, the signs of (1/ X (x))( X (x)/   ) are opposite.

Roof Truss.
In order to test the applicability of the proposed method for problems with more random variables and expressed in a more engineering way, the roof truss is selected as example.The truss is simply illustrated as in Figure 6.The top chord and the compression bars of the truss are reinforced by concrete and the bottom chord and the tension bars are all made of steel.Assume the truss bears uniformly distributed load , which can be transformed into nodal load  = /4.The perpendicular deflection of truss peak node C Δ  can be calculated using the knowledge of structural mechanics and where   and   are the cross-sectional areas of the reinforced concrete and steel bars, respectively,   and   are the corresponding elastic moduli of reinforced concrete and steel, and  is the length of the truss as Figure 6 shows.The distribution parameters of the independent normal random basic variables are given in Table 5.
Considering the safety of the truss, the perpendicular deflection Δ  should satisfy the constraint Δ  ≤ 0.03 m.Hence, the structural performance function can be given as follows: For this highly nonlinear example, the computational results of the importance measure on failure probability measures by ALK and MCS are listed in Tables 6 and 7.
In addition, the IMs on the failure probability by the state dependent parameter (SDP) method in literature [20] are also listed for comparison in Table 6.The MCS procedure takes 6 × 10 7 samples.The SDP method needs only 1024 model runs for calculating the importance measures.However, the Kriging method needs to call the performance function only 72 times.The ALK method begins with 12 training points, and 60 points are added into the initial DOE using ALK.The tendency of the true values of the training points for roof truss by ALK method is plotted in Figure 7. Compared with the SDP method, the proposed method has obtained a more accurate result and lower computational cost.
Additionally, it can be seen from Table 6 that the effects of the basic variable of  on the failure probability are notable and the influences of basic variables sectional area   and    on the failure probability are relatively smaller, whereas elastic moduli   and   and the length  are the least influential ones which can attract less attention.In engineering, parameters of the load  and the length  are more easily controlled or modified than other inputs, so the paper especially pays attention to their parametric sensitivities of IM     .Due to space limitation, Table 7 lists the parametric sensitivities of basic variables IM on the failure probability.Taking IM    for example, the distribution parameters of the load  are relatively higher than those of the length .Thus, it can be seen that changing those parameters with high parametric sensitivity has more influences on the IM of the corresponding variables than other parameters.The parametric sensitivity of     is valuable in the reliability engineering because it can provide indirect information for reliability design and reliability-based optimization., respectively.Joint 2 is also subjected to a horizontal load  3 = 10 KN.The input variables are all normally distributed, and their coefficient of variance is 0.05.We assume the displacement of node 3 in vertical direction not exceeding 3 mm as the constrain condition.The limit state function can be constructed,  = 3 mm − | 2 |, where  2 is an implicit function of the basic random variables.As shown in Figure 9, the finite element model can be obtained in Ansys 11.0.The results of the importance measures   of inputs are listed in Table 8.To identify the influential distribution parameters,  the computational results of parametric sensitivity are shown in Table 9.

A Planar 10-Bar
From Tables 8 and 9, it can be seen that the computational results of Monte Carlo method and ALK method are in good agreement.The MCS procedure takes 6 × 10 7 samples and its results can be seen as the accurate results, but the computational cost of Monte Carlo method is tremendously large especially for FEM.The ALK method only needs to call the finite element model 93 times to satisfy the accuracy.The efficient Kriging method begins to get convergence at 12 training points, and 81 points are added into the initial DOE using ALK.The true values of the added training points for roof truss by the ALK method are shown in Figure 10.The ALK method can also improve the computational efficiency impressively and ensure acceptable accuracy for the implicit finite element model.
As revealed in Table 8, the ranking of IM not to exceed as 3 mm the constrain condition   is as follows: the basic variables of , , and  on the failure probability are notable; the influences of basic variable  1 on the failure probability are less; and the influences of basic variables  2 and  3 on the failure probability are very small, which are even near zeros.From Table 9, the distribution parameters   ,   ,   , and   are the most influential ones on the IMs.Thus, it can be seen that changing those parameters with high parametric sensitivities has more influences on the IM of the corresponding variable than other parameters.In the sight of this structural design, we need to pay more attention to the distribution parameters with high parametric sensitivity.In order to obtain the IM on failure probability results, especially for the high ranking IM, the important distribution parameters to them must be given precisely.To do this, it is necessary to collect the information and improve the understanding of the distribution parameters.

Conclusions
This paper investigates the influence of the distribution parameters on the IM on failure probability.It is noted that the IM of basic variable not only is influenced by its distribution parameters but also is influenced by other basic variables' distribution parameters.By further developing the presented moment-independent IM on failure probability, the parametric sensitivity of IM is first presented according to the derivative theory; thus, how the influential distribution parameters influence the influential IM can be made clear.Meanwhile, we can decrease the variability of the IM on failure probability by collecting the information and improving the understanding of those most influential parameters.The parametric sensitivity of     is valuable in reliability engineering because it can provide direct   and useful information for reliability design and reliabilitybased optimization.
The computation of the IM on failure probability and its parametric sensitivity is often feasible by the MCS, but the computational cost of MCS method is tremendously large with small failure probability (10 −3 -10 −4 or smaller).For dealing with this problem, the ALK method is employed to calculate the IM and its parametric sensitivity.It can be seen by the numerical and engineering examples that the ALK method is more efficient than MC method.To ensure the computational accuracy, the large number of training points used in the traditional Kriging method is essential.Thanks to the existence of active learning process, the points which may greatly affect the metamodel's fitting accuracy can be precisely selected, which can make the Kriging metamodel more accurate and the additional computational cost is acceptable.It is noticed that a small quantity of points in the inter- The training point iteration that the proposed method is validated to be rational and efficient.

Figure 2 :
Figure 2: Risk of the sign of (X) to be wrongly predicted in a Kriging model: (a) the sign of (X) is negative; (b) the sign of (X) is positive.

Figure 3 :
Figure 3: Iteration history of the true value of added training points for Example 1 by ALK.

Figure 4 :
Figure 4: Sign of the response at each sample predicted for Example 1 by ALK.

Figure 5 :
Figure 5: Iteration history of the true value of added training points for Ishigami function by ALK.

Figure 6 :
Figure 6: Schematic diagram of a roof truss.
Structure.A planar 10-bar structure shown in Figure 8 is investigated.The horizontal bars have the same length  = 1 m.The diagonal bars and the vertical members have the same length √ 2.The cross-sectional area of 10 bars is denoted by  = 0.001 m 2 .The elastic modulus of all bars is  = 100 GPa. 1 = 80KN and  2 = 10 KN are the external loads subjected to joints 4 and 2

Figure 7 :
Figure 7: Iteration history of the true value of added training points for a roof truss by ALK.
esting region are added to construct the Kriging predictor model until the Kriging model satisfies necessary accuracy.The computational results of several examples demonstrate

Figure 10 :
Figure 10: Iteration histories of the true value of added training points for 10-bar structure.

The contributions of distribution parameter 𝜃 𝑗 𝑋 𝑖 on 𝑃 𝑓 𝑌 and 𝑃 𝑓 𝑌|X 𝑖 can be shown in Figure 1. To analyze the effect of changing the 𝑗th parameter 𝜃 𝑗 𝑋 𝑖 (𝑖 = 1, 2 ⋅ ⋅ ⋅ 𝑛) of input 𝑋 𝑖 , the sensitivity derivative of the IM on failure probability 𝛿 𝑃 X 𝑖 can be defined as
, but it still influences    and   |X  . −  is a vector containing all distribution parameters of the input variables but   .Influence of the distribution parameter on the    and   |X  .
(24)he initial step, only a dozen of points in the DOE are enough according to our experience.(b)ConstructtheKrigingmodel,compute the ERF of candidate points by (22) and(24), and judge whether it is smaller than the given tolerance.If so, the active learning process can stop and turn to Step 3 directly.If the stopping condition is not satisfied, it should go to Step 2(c) to add new training point.(c) Add the point with max value of ERF to the DOE and loop back to Step 2(b).Based on the active learning Kriging model, the corresponding  values of Kriging predictions y  =   (x  ) ( = 1, 2 ⋅ ⋅ ⋅ ) are obtained.Compute the values I   ( = 1, 2 ⋅ ⋅ ⋅ ) of the indicator function of this failure domain and    = ∑  =1    /.Using the same active learning Kriging model,    /    and    /    can be computed combining the reliability sensitivity analysis method by (10).By fixing the (, )th component ( , ) of X, generate  samples k    = ( 1 , . . .,  −1, ,   ,  +1, ⋅ ⋅ ⋅   ) ( = 1, 2 ⋅ ⋅ ⋅ ) according to the conditional PDF  X (x |   ), and the corresponding  values of Kriging predictions y   , =   (k   , )( = 1, 2 ⋅ ⋅ ⋅ ) can be obtained based on the same active learning Kriging model.With these  samples y    ( = 1, 2 ⋅ ⋅ ⋅ ), the    |  of conditional on   can be obtained as 2. Construct an active learning Kriging model as follows.(a)The samples x  in Step 1 are set as candidate points.Randomly choose some training points in the samples x  and evaluate the corresponding performance function.

Table 1 :
Computational results of the importance measure     in Example 1.

Table 2 :
The parametric sensitivities of the importance measure     in Example 1.

Table 3 :
Computational results of the importance measure     in Ishigami function.

Table 4 :
The parametric sensitivities of the importance measure     in Ishigami function.

Table 5 :
Statistical properties of random variables for roof truss.

Table 6 :
Computational results of the importance measure     for roof truss.

Table 7 :
The parametric sensitivities of the importance measure     for roof truss.

Table 8 :
Computational results of the importance measure     in 10-bar structure.

Table 9 :
The parametric sensitivities of the importance measure     in 10-bar structure.