Next Article in Journal
Non-Invasive Identification of Vehicle Suspension Parameters: A Methodology Based on Synthetic Data Analysis
Previous Article in Journal
Convergence Criteria for Fixed Point Problems and Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Derivative-Variance Hybrid Global Sensitivity Measure with Optimal Sampling Method Selection

1
School of Mechanical and Electrical Engineering, Wuhan Institute of Technology, Wuhan 430205, China
2
Automotive Institute, Wuhan Technical College of Communication, Wuhan 430065, China
3
Wuhan Fiberhome Technical Services Co., Ltd., Wuhan 430205, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(3), 396; https://doi.org/10.3390/math12030396
Submission received: 20 October 2023 / Revised: 21 January 2024 / Accepted: 22 January 2024 / Published: 26 January 2024
(This article belongs to the Section Engineering Mathematics)

Abstract

:
This paper proposes a derivative-variance hybrid global sensitivity measure with optimal sampling method selection. The proposed sensitivity measure is as computationally efficient as the derivative-based global sensitivity measure, which also serves as the conservative estimation of the corresponding variance-based global sensitivity measure. Moreover, the optimal sampling method for the proposed sensitivity measure is studied. In search of the optimal sampling method, we investigated the performances of six widely used sampling methods, namely Monte Carlo sampling, Latin hypercube sampling, stratified sampling, Latinized stratified sampling, and quasi-Monte Carlo sampling using the Sobol and Halton sequences. In addition, the proposed sensitivity measure is validated through its application to a rural bridge.

1. Introduction

Sensitivity analysis serves as a fundamental instrument for importance ranking and the dimension reduction of input parameters [1], and its applications are spreading across diverse topics in civil engineering, such as strength analysis [2], failure analysis [3], rock mechanics [4], slope stability analysis [5], and structural reliability analysis [6].
Sensitivity analysis examines how the output, or model’s response, is affected by the inputs and therefore calculates their corresponding importance. Compared to qualitative or subjective sensitivity measures [7,8], quantitative sensitivity measures have been demonstrated to be much more trustworthy [9,10,11]. Quantitative sensitivity measures can be roughly classified into local sensitivity measures and global sensitivity measures [12].
The local sensitivity measures, also referred to as the elementary effect, assess how the input factor locally affects the model output by varying the input of interest while keeping the rest of the inputs constant, which is the one factor of a time method [13]. The main drawback of local sensitivity measures is that they are fixed to the corresponding nominal points in the input space, which means they cannot provide information on the inputs’ total effects or detect interactions among inputs. In contrast to local sensitivity measures, global sensitivity measures evaluate the inputs’ overall effects on the model output, which can also account for interactions between inputs. Obviously, global sensitivity measures offer a more comprehensive sensitivity analysis and are thus more widely accepted in practical engineering. Invented in the 1970s, the Fourier amplitude sensitivity test [14] is one of the quickest global sensitivity measures. It is, however, merely a first-order measure and can be insufficient if higher-order terms dominate the model output [9]. Another widely adopted global sensitivity measure is the variance-based global sensitivity measure, which was first proposed by Sobol and further developed by Saltelli and his colleagues [1,11]. The variance-based global sensitivity measure calculates the model output’s variances of different subsets of input variables, which usually requires a large number of model evaluations in order to achieve a reasonable convergence and then becomes computationally unaffordable. To this end, a number of alternative global sensitivity measures have been recently proposed, one of which is the Morris measure [7,15]. The Morris measure relies on averaging the absolute values of local sensitivity measures, which is considerably cheaper in terms of computation than the variance-based global sensitivity measure. Nevertheless, the Morris method uses local sensitivity measures of random points on a fixed grid for averaging, which cannot correctly reflect effects if the model’s characteristic dimensions are smaller than the grid size [16]. To solve this issue, the derivative-based global sensitivity measure was recently proposed by employing local sensitivity measures of randomly distributed points rather than a fixed grid to estimate the integral representing the global effects [17]. The derivative-based global sensitivity measure has been demonstrated to be good at quickly identifying unimportant input factors, while the variance-based global sensitivity measure is more comprehensive but also more computationally expensive. Another advantage of the variance-based global sensitivity measure over the derivative-based global sensitivity measure is that the variance-based global sensitivity measure provides absolute values for the sensitivity measures, while the derivative-based global sensitivity measure only provides relative values for sensitivity measures. In other words, in using a derivative-based global sensitivity measure for the sensitivity analysis, one has to compute the derivative-based global sensitivity measure of each input before determining their sensitivities or importance. Recently, a variety of innovative sensitivity measures have been proposed for various application scopes. A novel global sensitivity measure, called PAWN, was proposed to efficiently compute density-based sensitivity indices [18]. As for structural engineering, a novel vine copula-based approach for analyzing the variance-based global sensitivity measure is proposed [19]. In reliability engineering, a new reliability-oriented global sensitivity measure, which is obtained by replacing the variance with entropy, was proposed [20]. Connections between entropy and variance in reliability engineering have also been recently explored [21].
It should also be pointed out that the computational efficiency of global sensitivity measures heavily depends on different sampling methods [12,22,23]. Although recently developed meta-models appear to be more efficient in terms of sampling [24,25,26], the traditional quasi-Monte Carlo or Monte Carlo methods remain the primary sampling methods for practical sensitivity analysis. It is well known that quasi-Monte Carlo sampling should outperform Monte Carlo sampling in a sensitivity analysis for a sufficiently large sample size [23,27]. A board comparison including other popular sampling methods, such as Latin hypercube sampling, stratified sampling, and Latinized stratified sampling, has not been conducted yet. In addition, the impact of nonlinear dependency among random variables is also important in sensitivity analysis, which can be investigated through the copula theory [28].
Sensitivity measures that combine the advantages of both the variance-based global sensitivity measures and derivative-based global sensitivity measures would provide better tools for sensitivity analysis. In light of such an idea, a new global sensitivity measure based on derivative-integral and variance decomposition was recently proposed; however, this measure is too complicated for many engineering applications [29], and a simpler derivate-variance hybrid sensitivity analysis is needed. Therefore, this research presents a simple derivative-variance hybrid sensitivity measure, which we named the derivative-variance hybrid global sensitivity measure. The proposed derivative-variance hybrid global sensitivity measure is as computationally efficient as the derivative-based global sensitivity measure, which can also serve as a conservative approximation of the corresponding variance-based global sensitivity measure. Furthermore, in seeking the optimal sampling method, we investigated the performances of six frequently used sampling methods, namely Monte Carlo samplings, Latin hypercube samplings, stratified samplings, Latinized stratified samplings, and two quasi-Monte Carlo samplings using the Sobol sequence and the Halton sequence.
The rest of this paper is organized as follows: Section 2 reviews the necessary backgrounds of sensitivity measures; Section 3 proposes the derivative-variance hybrid global sensitivity measure; Section 4 determines the optimal sampling method via an investigation of the performances of six widely used sampling methods; and Section 5 validates the proposed sensitivity measure through its application to the sensitivity analysis of a rural bridge.

2. Fundamentals of Sensitivity Measures

Quantitative sensitivity measures are used to calculate the influences of inputs on the model output as certain numbers. In this section, we review three important quantitative sensitivity measures; we consider a model output represented by the function f ( x ) = f ( x 1 , , x n ) defined in the unit hypercube I n = [ 0 , 1 ] n , where f ( x ) is both differentiable and square-integrable.

2.1. Local Sensitivity Measure

For local sensitivity measures, the sensitivity of a variable is measured locally via a nominal point in the space of the inputs, and the value of this variable is represented by the corresponding partial derivative or the approximated partial difference, i.e.,
ε i ( x ) = f ( x ) x i f ( x + Δ i ) f ( x ) Δ i
where εi is the local sensitivity measure associated with the nominal point x0 and Δ i = [ 0 , , Δ , , 0 ] is a vector of length ∆ along the direction of xi.

2.2. Derivative-Based Global Sensitivity Measure

In light of the fact that a high value of the output’s derivative indicates a large variation in the model output, the derivative-based global sensitivity measures are achieved globally via stochastic integrals of the randomly distributed local sensitivity measures. There exist various definitions of the derivative-based global sensitivity measures. Here we present a widely used one:
v i = I n ε i ( x ) 2 d x 1 N k = 1 N ε i ( x k ) 2
where xk is obtained through random sampling, N is the total number of sample points, and vi is Sobol’s derivative-based global sensitivity measure [15]. It is noted that a derivative-based global sensitivity measure does not reflect the relative importance of xi. Thus, one has to rank all the inputs according to their derivative-based global sensitivity measures to obtain their relative importance, while the value of a derivative-based global sensitivity measure itself does not provide any useful information on sensitivity.

2.3. Variance-Based Global Sensitivity Measure

Compared to various versions of the derivative-based global sensitivity measures, the definition of variance-based global sensitivity measures is unique, and their core is the analysis of variance decomposition [30]. The analysis of variance decomposition of f ( x ) has the following form:
f ( x ) = f 0 + i = 0 n f i + 1 i < j n f i j ( x i , x j ) + + f 1 n ( x i , , x n )
which should satisfy
0 1 f j i ( x j , , x i ) d x k = 0
where 1 ≤ jkin.
Each term in Equation (3) can be calculated in the following recursive manner:
f 0 = I n f x d x
f i ( x i ) = I n f x k i d x k f 0
f i j ( x i , x j ) = I n f x k i , j d x k f 0 f i ( x i ) f j ( x j )
f i j ( x i , , x j ) = I n f x k i , j d x k f 0 k = i j f k ( x k ) k = i j l = i , l k j f k l ( x k , x k ) k = i j f i ( k 1 ) ( k + 1 ) j ( x i , , x k 1 , x k + 1 , , x j )
Let x A = x A 1 , , x A s denote an arbitrary subset of x = x 1 , , x n , where 1 ≤ A1 < Asn and 1 ≤ sn, and let the corresponding complement denote x B = x B 1 , , x B n s , which satisfies x = x A , x B and A 1 , , A s , B 1 , , B n s = 1 , 2 , , n . The variance with respect to xA is defined as
V x A = I n f i j ( x i , , x j ) 2 x i x j .
The total variance is defined by
V = s = 1 n A 1 < < A s s V A 1 A s .
Finally, the variance-based global sensitivity measure with respect to xA is defined as
S x A = V x A V .
The total variance-based global sensitivity measure with respect to xA is defined as
S x A T = 1 S x B .
The first-order variance-based global sensitivity measure with respect to xi is defined as
S i = V i V .
where Vi is referred to as the first-order variance.
And the first-order total variance-based global sensitivity measure with respect to xi is computed as
S i T = 1 V 1 ( i 1 ) ( i + 1 ) n V
Sobol proposed a Monte Carlo sampling-based computation method for the variance-based global sensitivity measure [11], which consists of the following four steps: Firstly, construct two input matrixes U and V, i.e.,
U = x 11 x 1 n x N 1 x N n
V = x 11 x 1 n x N 1 x N n
Secondly, according to the relationship of xA and xB to x, compute UA and UB from U, and compute VA and VB from V, i.e.,
U A = x 1 A 1 x 1 A s x N A 1 x N A s
U B = x 1 B 1 x 1 B n s x N B 1 x N B n s
V A = x 1 A 1 x 1 A s x N A 1 x N A s
V B = x 1 B 1 x 1 B n s x N B 1 x N B n s
Thirdly, composite the following two matrices: WA = [UA,VB] and WB = [VA,UB]. Lastly, compute the variance-based global sensitivity measure using Equations (9)–(14), and apply the following stochastic integrals:
V x A = I n f ( x ) f ( x A , x B ) d x d x B f 0 2 1 N i = 1 N f ( U A ) f ( W A ) f 0 2
V x B = I n f ( x ) f ( x A , x B ) d x d x A f 0 2 1 N i = 1 N f ( U A ) f ( W B ) f 0 2
V = I n f ( x ) 2 d x f 0 2 1 N i = 1 N f ( U A ) 2 f 0 2
f 0 = I n f ( x ) d x 1 N i = 1 N f ( U A )
where x’ is another input space that is independent from x.
If the aim is to compute the total of the variance-based global sensitivity measure, the computational cost is nN in terms of model evaluations, which linearly increases with both the sample size and the dimensions of the inputs. What hampers many engineering applications of the variance-based global sensitivity measure is the fact that the sample size N required is usually too large to afford, which is significantly greater than that of the derivative-based global sensitivity measures [31]. In addition, the uniformity and stochasticity of the samples can dramatically influence the convergence rate [32].

3. The Derivative-Variance Hybrid Global Sensitivity Measure

In this section, the derivative-variance hybrid global sensitivity measure is proposed, which is as computationally efficient as the derivative-based global sensitivity measure while providing as much meaningful absolute value for sensitivity analysis as the variance-based global sensitivity measure. Section 3.1 proposes the sensitivity measure; Section 3.2 presents a method for computing the proposed sensitivity measure; Section 3.3 compares the derivative-variance hybrid global sensitivity measure with the derivative-based global sensitivity measure and the variance-based global sensitivity measure using three benchmark examples; and Section 3.4 discusses the advantages and limitations of the proposed sensitivity measure.

3.1. Presentation of the Sensitivity Measure

The derivative-based global sensitivity measure and the variance-based global sensitivity measure can be linked using [17]:
π 2 V S i T v i
where S i T is the variance-based global sensitivity measure computed via Equation (14), V is the total variance computed via Equation (23), and vi is the derivative-based global sensitivity measure computed via Equation (2).
On the basis of Equation (25), the following sensitivity measure, referred to as the derivative-variance hybrid global sensitivity measure, is proposed:
ϒ i = ν i π 2 V S i T .
Evidently, the derivative-variance hybrid global sensitivity measure ϒi is obtained by dividing the derivative-based global sensitivity measure vi by π2V, so they should share a similar computational cost that is much lower than that of the variance-based global sensitivity measure. Derivative-variance hybrid global sensitivity measures can be treated as reduced derivative-based global sensitivity measures in which the reduction allows the value of the derivative-variance hybrid global sensitivity measure to serve as a conservative approximation of the corresponding variance-based global sensitivity measure.

3.2. Computational Method

As shown in Figure 1, the derivative-variance hybrid global sensitivity measure can be computed using the following steps: (1) sample the input matrix x; (2) compute the total variance V using Equations (23) and (24); (3) enumerate each input xi with Steps 4 and 5; (4) compute the derivative-based global sensitivity measure vi using Equation (2); and (5) compute the derivative-variance hybrid global sensitivity measure ϒi using Equation (26).

3.3. Comparison with Other Global Sensitivity Measures

In this subsection, the performances of the derivative-variance hybrid global sensitivity measure, the derivative-based global sensitivity measure, and the variance-based global sensitivity measure using three benchmark examples are compared and analyzed, and the corresponding results are obtained through Monte Carlo simulations, in which the sample size is 108.
Example 1. 
The analysis of variance decomposition is provided below [33]:
f ( x ) = i = 1 3 c i x i 1 2 + c 12 x 1 1 2 x 2 1 2 5
where 0 < xi < 1, c1 = 10, c2 = 1, c3 = 0.1, and c12 = 10. In this example, x1 is the most sensitive input, x3 is trivial, the importance of x2 is in-between, and the sensitivity of a highly non-linear term is controlled by c12. The corresponding results are shown in Table 1, from which it is seen that: (1) the derivative-based global sensitivity measures cannot directly indicate sensitivities of the corresponding inputs, while the variance-based global sensitivity measure and the derivative-variance hybrid global sensitivity measure can; (2) the derivative-variance hybrid global sensitivity measure can serve as a conservative estimation of the corresponding variance-based global sensitivity measure since the value of the derivative-variance hybrid global sensitivity measure is slightly larger than that of the corresponding variance-based global sensitivity measure; the value of derivative-variance hybrid global sensitivity measure can be greater than 1.
Example 2. 
The g-function is provided below [33]:
f ( x ) = i = 1 n 4 x i 2 + a i 1 + a i
where 0 < xi < 1, n = 3, a1 = 1, a2 = 10, and a3 = 100. In this example, a larger value of ai corresponds to a smaller importance of xi. The corresponding results are shown in Table 2. It is seen from Table 2 that the derivative-variance hybrid global sensitivity measure can still effectively identify insensitive input x3; however, the ratio between the derivative-variance hybrid global sensitivity measure and the corresponding variance-based global sensitivity measure can be significantly greater than 1.
Example 3. 
The Ishigami function is provided below [34]:
f ( x ) = sin 2 π x 1 π + a sin 2 π x 2 π 2 + b 2 π x 3 π 4 sin 2 π x 1 π
where 0 < xi < 1, i = 1, 2, or 3, and a and b in the following examples were set to 7 and 0.1, respectively. The corresponding results are shown in Table 3, from which it is seen that the derivative-variance hybrid global sensitivity measure may not be efficient with regard to the ranking of inputs of significant sensitivities.

3.4. Advantages and Limitations

Since the derivative-variance hybrid global sensitivity measure is the conservative approximation of the corresponding variance-based global sensitivity measure, the values of the variance-based global sensitivity measures can directly reflect the inputs’ sensitivities, which provides a solid argument for identifying insensitive inputs. For instance, if it is found that ϒi is negligibly small, then xi can be safely reduced without further analyzing the rest of the inputs. For the identification of insensitive inputs using the derivative-based global sensitivity measures, one has to exhaust each input and rank their derivative-based global sensitivity measures in ascending order. On the other hand, the computational cost of the derivative-variance hybrid global sensitivity measure is similar to that of the derivative-based global sensitivity measure, which is considerably lower than that of the variance-based global sensitivity measure. At this point, we list the advantages of the proposed derivative-variance hybrid global sensitivity measure as follows: (1) inputs with a small derivative-variance hybrid global sensitivity measure can be labeled as unimportant without investigating the rest of the inputs; (2) the computational cost of the derivative-variance hybrid global sensitivity measure is comparable to that of the derivative-based global sensitivity measure, which only requires the additional computation of the total variance V; (3) since the derivative-variance hybrid global sensitivity measure is the upper bound for the corresponding variance-based global sensitivity measure, the derivative-based global sensitivity measures can serve as a conservative approximation of the variance-based global sensitivity measure. Thus, the derivative-variance hybrid global sensitivity measure is recommended to replace derivative-based global sensitivity measure or the variance-based global sensitivity measure for sensitivity analysis.
On the other hand, since the computation of the derivative-variance hybrid global sensitivity measure requires the value of the corresponding derivative-based global sensitivity measure, it shares the same limitations as the derivative-based global sensitivity measure, such as: (1) it requires the output function to be continuous so that the corresponding derivatives are available; (2) its applications on highly nonlinear output functions may result in considerable errors; and (3) it cannot detect the interactive effect between different inputs.

4. Selection of the Optimal Sampling Method

In this section, six frequently used sampling methods are evaluated in terms of their performances of sensitivity analysis using the derivative-variance hybrid global sensitivity measure, namely: Monte Carlo sampling [23]; Latin hypercube sampling [12]; stratified sampling [22]; Latinized stratified sampling [22]; and two quasi-Monte Carlo samplings using the Sobol sequence and Halton sequence [23]. Figure 2 presents one realization of random points in domain I2 that are generated from the above six sampling methods. It is noted that, when the simple size N equals kn and k is a positive integer, stratified sampling and Latinized stratified sampling can guarantee that there is exactly one point contained in each stratum (a hypercube with a length of 1/k defined by a corresponding grid).
The performances of the derivative-variance hybrid global sensitivity measure of the aforementioned six sampling methods were investigated using the aforementioned three examples described by Equations (27)–(29). Since all these sampling methods produce unbiased estimations for random integrals, only the stochastic errors were studied, and we chose the standard deviation of the relative error as the index for the stochastic errors. At this point, the relative error of variable χ is defined as: r ( χ i ) = χ i χ 0 / χ 0 , where χ i is one realization of χ and χ 0 is the corresponding analytical solution. Figure 3 shows the standard deviations (STD) of the relative errors for the derivative-variance hybrid global sensitivity measure ϒi as functions of sample size N, in which N equals 103, 253, and 503. In the calculation, 107 sample points were divided into the corresponding number of sample sets according to a different sample size N, and each sample set yields one group of relative errors. As can be seen from Figure 3, Monte Carlo sampling and Latin hypercube sampling always have the largest stochastic errors among these six sampling methods, suggesting a ban of these two sampling methods on sensitivity analysis. Figure 3 presents the stochastic errors for derivative-variance hybrid global sensitivity measure evaluations. It is noted that, due to the derivative of x3 of Example 1 being a constant, the errors’ STDs obtained using different sampling methods are equal to zero. In Figure 3, it can be seen that the Latinized stratified sampling and the two quasi-Monte Carlo samplings produce the lowest stochastic error, and their relative performances may vary for different cases. In this instance, the two quasi-Monte Carlo samplings outperform the Latinized stratified sampling in the evaluations of ϒ1 and ϒ3, while the Latinized stratified sampling provides the lowest stochastic error for the evaluations of ϒ2. Thus, Latinized stratified sampling and quasi-Monte Carlo sampling are recommended.
An explanation for the above phenomenon can involve the sample’s stochasticity and uniformity. On the one hand, a random integral requires the sample to have good uniformity to cover sudden changes in small areas. On the other hand, reasonable stochasticity is also necessary to avoid periodic related issues. Nevertheless, the dilemma is that high uniformity and high stochasticity cannot be achieved simultaneously since the uniformity decreases as the stochasticity increases. Monte Carlo sampling and Latin hypercube sampling have poor uniformity, resulting in undesirable performances. In contrast, quasi-Monte Carlo samplings and Latinized stratified sampling have the right balance between uniformity and stochasticity, which contributes to their satisfactory accuracy and convergence rate in sensitivity analysis.

5. Numerical Application and Validation

In order to demonstrate its validity and practicability, the proposed sensitivity measure is applied to the sensitivity analysis of a rural bridge. The bridge, located in the Xinzhou district (Wuhan city, Hubei Province, China), is a solid plate bridge with a span-combination of 4 × 6 m, whose width equals 10 m (shown in Figure 4). Each span consists of 10 solid plates of 0.99 m wide.
The performance function [35] of the weakest plate is chosen as the model output, which can be expressed as
M = f s d A s h a s f s d A s 2 f c d b b h ρ 1 l 2 / 8 b h ρ 2 l 2 / 8 F α l / 4
in which the specific parameters of the bridge are obtained from in situ investigations and reasonable error analysis (shown in Table 4).
Firstly, the derivative-variance hybrid global sensitivity measure is applied to the sensitivity analysis of the rural bridge. In the analysis, the sampling method used was quasi-Monte Carlo sampling with Sobol sequence, and the convergence threshold in terms of absolute error was set to 0.01. The values of the derivative-based global sensitivity measures obtained from 108 samples are considered the target for error computation. As shown in Table 5, it was found that 1000 samples can already guarantee the above criterion of convergence, and such a good convergence using merely 1000 samples indicates the convergence rate of derivative-variance hybrid global sensitivity measure using quasi-Monte Carlo sampling is sufficiently high. In addition, it is seen from Table 5 that the six insensitive inputs, namely l, b, h , ρ 1 , ρ 2 , and f c d , are found (highlighted in gray), and their derivative-variance hybrid global sensitivity measures are less than the corresponding sensitivity threshold of 1.5%. At this juncture, it is noted that the relative errors for the different inputs’ derivative-variance hybrid global sensitivity measures are surprisingly close to each other, and the reason for this remains to be further studied.
On the other hand, in order to compare the derivative-variance hybrid global sensitivity measure and the variance-based global sensitivity measure in terms of computational efficiency, the variance-based global sensitivity measure was also imposed on the sensitivity analysis of the above example, and the identical convergence threshold of 0.01 was adopted. The target variance-based global sensitivity measure values for error computation were obtained from 108 samples. It is found that, if the sample size equals 240,000, the above criterion for convergence was just met, and the corresponding result is shown in Table 6. By comparing Table 5 and Table 6, it is evident that the computational efficiency of the derivative-variance hybrid global sensitivity measure is much higher than that of the variance-based global sensitivity measure, and the required sample sizes for the derivative-variance hybrid global sensitivity measure and the variance-based global sensitivity measure to meet the convergence threshold of 0.01 are 1000 and 240000, respectively.
Lastly, in order to demonstrate that the derivative-variance hybrid global sensitivity measures can serve as good approximations for the corresponding variance-based global sensitivity measures, the derivative-variance hybrid global sensitivity measures, the derivative-based global sensitivity measures, and the variance-based global sensitivity measures for the above example are computed and compared, and the sample sizes for computing the variance-based global sensitivity measures, the derivative-based global sensitivity measures, and the derivative-variance hybrid global sensitivity measures are 240,000, 1000, and 1000, respectively. Table 7 shows the comparison among the derivative-based global sensitivity measures, the variance-based global sensitivity measures, and derivative-variance hybrid global sensitivity measures. From Table 7, it is seen that: (1) the absolute values of the derivative-based global sensitivity measures are meaningless, so one must compute all the inputs’ derivative-based global sensitivity measures and rank them before determining their sensitivities, while the values of the derivative-variance hybrid global sensitivity measures and the variance-based global sensitivity measures directly reflect the sensitivity of the corresponding inputs; (2) for sensitive inputs, the derivative-variance hybrid global sensitivity measure is consistently greater than the corresponding variance-based global sensitivity measure; and (3) for sensitive inputs, the relative errors between the derivative-based global sensitivity measures and the corresponding variance-based global sensitivity measures are all less than 26%. Therefore, it is demonstrated that derivative-variance hybrid global sensitivity measures are good and conservative approximations for the variance-based global sensitivity measures.

6. Conclusions

This research proposes a new sensitivity measure for sensitivity analysis, which is referred to as the derivative-variance hybrid global sensitivity measure. Optimal sampling methods for computing the derivative-variance hybrid global sensitivity measure are investigated by evaluating the corresponding performances of six frequently used sampling methods. To demonstrate its validity and practicality, the proposed sensitivity measure is applied to the sensitivity analysis of a rural bridge, which is further compared to the derivative-based global sensitivity measure and the variance-based global sensitivity measure. We found that the derivative-variance hybrid global sensitivity measure can directly measure the sensitivities of the corresponding inputs of the variance-based global sensitivity measure, while the derivative-based global sensitivity measures cannot directly reflect the inputs’ sensitivities, and the computational efficiency of the derivative-variance hybrid global sensitivity measure is much higher than that of the variance-based global sensitivity measure. Thus, the derivative-variance hybrid global sensitivity measure has promising sensitivity analysis applications in a variety of engineering fields, such as strength analysis, failure analysis, rock mechanics, slope stability analysis, and structural reliability analysis. In addition, the following conclusions are drawn:
  • The derivative-variance hybrid global sensitivity measure and the derivative-based global sensitivity measure share similar computational costs, but the values of the derivative-variance hybrid global sensitivity measures directly imply the sensitivities of the corresponding inputs, while the values of the derivative-based global sensitivity measures are meaningless.
  • The derivative-variance hybrid global sensitivity measure can serve as a good and conservative approximation for the corresponding variance-based global sensitivity measure, while the computational cost of the derivative-variance hybrid global sensitivity measure is much lower than that of the variance-based global sensitivity measure.
  • Latinized stratified sampling and quasi-Monte Carlo sampling are recommended for computing the derivative-variance hybrid global sensitivity measure.

Author Contributions

Conceptualization, H.L.; Methodology, J.C.; Software, J.H.; Investigation, A.X.; Writing—original draft, J.L.; Writing—review & editing, C.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Gusu Innovation and Entrepreneurship Leading Talent Program (No. ZXL2022518) and Wuhan’s 2022 Provincial Service Industry Development Guidance Fund 282 Plan (Services of Wuhan Municipal Development and Reform Commission (2023) No. 120).

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors wish to acknowledge the financial support provided by the Gusu Innovation and Entrepreneurship Leading Talent Program (No. ZXL2022518) and Wuhan’s 2022 Provincial Service Industry Development Guidance Fund 282 Plan (Services of Wuhan Municipal Development and Reform Commission (2023) No. 120).

Conflicts of Interest

Authors Aibo Xu and Jiwei Hu were employed by the Wuhan Fiberhome Technical Services Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The Wuhan Fiberhome Technical Services Co., Ltd. had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Saltelli, A. Sensitivity analysis for importance assessment. Risk Anal. 2002, 22, 579–590. [Google Scholar] [CrossRef]
  2. Feng, D.-C.; Fu, B. Shear strength of internal reinforced concrete beam-column joints: Intelligent modeling approach and sensitivity analysis. Adv. Civ. Eng. 2020, 2020, 8850417. [Google Scholar] [CrossRef]
  3. Jolfaei, S.; Lakirouhani, A. Sensitivity Analysis of Effective Parameters in Borehole Failure, Using Neural Network. Adv. Civ. Eng. 2022, 2022, 4958004. [Google Scholar] [CrossRef]
  4. Wang, H.-S.; Shuang, H.-Q.; Li, L.; Xiao, S.-S. The Stability Factors’ Sensitivity Analysis of Key Rock B and Its Engineering Application of Gob-Side Entry Driving in Fully-Mechanized Caving Faces. Adv. Civ. Eng. 2021, 2021, 9963450. [Google Scholar] [CrossRef]
  5. Yuan, B.; Li, Z.; Su, Z.; Luo, Q.; Chen, M.; Zhao, Z. Sensitivity of multistage fill slope based on finite element model. Adv. Civ. Eng. 2021, 2021, 6622936. [Google Scholar] [CrossRef]
  6. Wei, P.; Tang, C.; Yang, Y. Structural reliability and reliability sensitivity analysis of extremely rare failure events by combining sampling and surrogate model methods. Proc. Inst. Mech. Eng. Part O J. Risk Reliab. 2019, 233, 943–957. [Google Scholar] [CrossRef]
  7. Campolongo, F.; Saltelli, A.; Cariboni, J. From screening to quantitative sensitivity analysis. A unified approach. Comput. Phys. Commun. 2011, 182, 978–988. [Google Scholar] [CrossRef]
  8. Hamby, D. A comparison of sensitivity analysis techniques. Health Phys. 1995, 68, 195–204. [Google Scholar] [CrossRef]
  9. Morio, J. Global and local sensitivity analysis methods for a physical system. Eur. J. Phys. 2011, 32, 1577–1583. [Google Scholar] [CrossRef]
  10. Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. Global Sensitivity Analysis: The Primer; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  11. Sobol, I.M. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Math. Comput. Simul. 2001, 55, 271–280. [Google Scholar] [CrossRef]
  12. 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]
  13. Saltelli, A.; Annoni, P. How to avoid a perfunctory sensitivity analysis. Environ. Model. Softw. 2010, 25, 1508–1517. [Google Scholar] [CrossRef]
  14. Cukier, R.; Levine, H.; Shuler, K. Nonlinear sensitivity analysis of multiparameter model systems. J. Comput. Phys. 1978, 26, 1–42. [Google Scholar] [CrossRef]
  15. Morris, M.D. Factorial Sampling Plans for Preliminary Computational Experiments. Technometrics 1991, 33, 161–174. [Google Scholar] [CrossRef]
  16. Kucherenko, S. Application of global sensitivity indices for measuring the effectiveness of quasi-monte carlo methods and parameter estimation. In Proceedings of the Fifth International Conference on Sensitivity Analysis of Model Output (SAMO 2007), Budapest, Hungary, 18–22 June 2007; pp. 35–36. [Google Scholar]
  17. Kucherenko, S.; Song, S. The derivative-based global sensitivity measures and Their Link with Sobol’ Sensitivity Indices. In Monte Carlo and Quasi-Monte Carlo Methods; Springer Proceedings in Mathematics & Statistics; Springer: Berlin/Heidelberg, Germany, 2016; pp. 455–469. [Google Scholar]
  18. Pianosi, F.; Wagener, T. A simple and efficient method for global sensitivity analysis based on cumulative distribution functions. Environ. Model. Softw. 2015, 67, 1–11. [Google Scholar] [CrossRef]
  19. Bai, Z.; Wei, H.; Xiao, Y.; Song, S.; Kucherenko, S. A vine copula-based global sensitivity analysis method for structures with multidimensional dependent variables. Mathematics 2021, 9, 2489. [Google Scholar] [CrossRef]
  20. Kala, Z. New importance measures based on failure probability in global sensitivity analysis of reliability. Mathematics 2021, 9, 2425. [Google Scholar] [CrossRef]
  21. Kala, Z. Quantification of Model Uncertainty Based on Variance and Entropy of Bernoulli Distribution. Mathematics 2022, 10, 3980. [Google Scholar] [CrossRef]
  22. Shields, M.D.; Zhang, J. The generalization of Latin hypercube sampling. Reliab. Eng. Syst. Saf. 2016, 148, 96–108. [Google Scholar] [CrossRef]
  23. Caflisch, R.E. Monte carlo and quasi-monte carlo methods. Acta Numer. 1998, 7, 1–49. [Google Scholar] [CrossRef]
  24. Sudret, B.; Mai, C.V. Computing the derivative-based global sensitivity measures using polynomial chaos expansions. Reliab. Eng. Syst. Saf. 2015, 134, 241–250. [Google Scholar] [CrossRef]
  25. Dubourg, V.; Sudret, B. Meta-model-based importance sampling for reliability sensitivity analysis. Struct. Saf. 2014, 49, 27–36. [Google Scholar] [CrossRef]
  26. Song, S.; Lu, Z.; Qiao, H. Subset simulation for structural reliability sensitivity analysis. Reliab. Eng. Syst. Saf. 2009, 94, 658–665. [Google Scholar] [CrossRef]
  27. Kucherenko, S.; Rodriguez-Fernandez, M.; Pantelides, C.; Shah, N. Monte Carlo evaluation of the derivative-based global sensitivity measures. Reliab. Eng. Syst. Saf. 2009, 94, 1135–1148. [Google Scholar] [CrossRef]
  28. 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]
  29. Liu, J.; Liu, Q.; Han, X.; Jiang, C.; Tao, Y. A new global sensitivity measure based on derivative-integral and variance decomposition and its application in structural crashworthiness. Struct. Multidiscip. Optim. 2019, 60, 2249–2264. [Google Scholar] [CrossRef]
  30. Sobol, I.M. Sensitivity analysis for non-linear mathematical models. Math. Model. Comput. Exp. 1993, 1, 407–414. [Google Scholar]
  31. Sobol, I.M.; Kucherenko, S. Derivative based global sensitivity measures. Procedia Soc. Behav. Sci. 2010, 2, 7745–7746. [Google Scholar] [CrossRef]
  32. De-xin, Z.; Yun-zheng, J.; Jian, C. Analytic sensitivity analysis for shape optimization. Appl. Math. Mech. 2001, 22, 1325–1332. [Google Scholar] [CrossRef]
  33. Sobol’, I.M.; Kucherenko, S. Derivative based global sensitivity measures and their link with global sensitivity indices. Math. Comput. Simul. 2009, 79, 3009–3017. [Google Scholar] [CrossRef]
  34. Kala, Z. Benchmark of goal oriented sensitivity analysis methods using Ishigami function. Int. J. Math. Comput. Methods 2018, 3. [Google Scholar]
  35. Valdebenito, M.A.; Schuëller, G.I. A survey on approaches for reliability-based optimization. Struct. Multidiscip. Optim. 2010, 42, 645–663. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the computation of derivative-variance hybrid global sensitivity measure.
Figure 1. Flowchart of the computation of derivative-variance hybrid global sensitivity measure.
Mathematics 12 00396 g001
Figure 2. One realization obtained using different sampling methods: (a) Monte Carlo sampling, (b) Latin hypercube sampling, (c) stratified sampling, (d) Latinized stratified sampling, (e) quasi-Monte Carlo sampling using Sobol sequence, and (f) quasi-Monte Carlo sampling using Halton sequence.
Figure 2. One realization obtained using different sampling methods: (a) Monte Carlo sampling, (b) Latin hypercube sampling, (c) stratified sampling, (d) Latinized stratified sampling, (e) quasi-Monte Carlo sampling using Sobol sequence, and (f) quasi-Monte Carlo sampling using Halton sequence.
Mathematics 12 00396 g002
Figure 3. Standard deviations of relative errors produced using different sampling methods: (a) ϒ1 of Example 1, (b) ϒ2 of Example 1, (c) ϒ3 of Example 1, (d) ϒ1 of Example 2, (e) ϒ2 of Example 2, (f) ϒ3 of Example 2, (g) ϒ1 of Example 3, (h) ϒ2 of Example 3, and (i) ϒ3 of Example 3.
Figure 3. Standard deviations of relative errors produced using different sampling methods: (a) ϒ1 of Example 1, (b) ϒ2 of Example 1, (c) ϒ3 of Example 1, (d) ϒ1 of Example 2, (e) ϒ2 of Example 2, (f) ϒ3 of Example 2, (g) ϒ1 of Example 3, (h) ϒ2 of Example 3, and (i) ϒ3 of Example 3.
Mathematics 12 00396 g003
Figure 4. Picture of the solid plate bridge.
Figure 4. Picture of the solid plate bridge.
Mathematics 12 00396 g004
Table 1. Results of the derivative-variance hybrid global sensitivity measure, the derivative-based global sensitivity measure, and the variance-based global sensitivity measure from Example 1.
Table 1. Results of the derivative-variance hybrid global sensitivity measure, the derivative-based global sensitivity measure, and the variance-based global sensitivity measure from Example 1.
Variance-Based Global Sensitivity MeasureDerivative-Based Global Sensitivity MeasureDerivative-Variance Hybrid Global Sensitivity Measure
x10.9901100.01.204
x20.01001.0910.0131
x39.657 × 10−50.0101.204 × 10−4
Table 2. Results of the derivative-variance hybrid global sensitivity measure, the derivative-based global sensitivity measure, and the variance-based global sensitivity measure from Example 2.
Table 2. Results of the derivative-variance hybrid global sensitivity measure, the derivative-based global sensitivity measure, and the variance-based global sensitivity measure from Example 2.
Variance-Based Global Sensitivity MeasureDerivative-Based Global Sensitivity MeasureDerivative-Variance Hybrid Global Sensitivity Measure
x10.96794.0114.7063
x20.03520.14330.1681
x34.125 × 10−40.00170.0020
Table 3. Results of the derivative-variance hybrid global sensitivity measure, the derivative-based global sensitivity measure, and the variance-based global sensitivity measure from Example 3.
Table 3. Results of the derivative-variance hybrid global sensitivity measure, the derivative-based global sensitivity measure, and the variance-based global sensitivity measure from Example 3.
Variance-Based Global Sensitivity MeasureDerivative-Based Global Sensitivity MeasureDerivative-Variance Hybrid Global Sensitivity Measure
x10.5574304.72.230
x20.4421967.27.078
x30.2433433.73.174
Table 4. List of the inputs.
Table 4. List of the inputs.
Input ParameterMedianDiscrepancyCoefficient of Variance
Tensile strength of reinforcements fsd (kPa)388,28026,6000.0685
Cross-sectional area of reinforcements As (m2)0.0028150.0000990.035
Thickness of reinforcement protection layer as (m)0.050.0030.06
Axial compressive strength of concrete fcd (kPa)31,20046800.15
Length of the solid plate l (m)5.60.0150.0027
Width of the solid plate b (m)0.990.010.01
Thickness of the solid plate h (m)0.320.010.031
Density of the solid plate ρ1 (kN/m3)2520.08
Thickness of the deck pavement h′ (m)0.10.010.1
Density of the deck pavement ρ2 (kN/m3)262.080.08
Vehicle load F (kN)300300.1
Transverse distribution coefficient α0.260.0260.1
Table 5. Sensitivity analysis of the bridge using the derivative-variance hybrid global sensitivity measure.
Table 5. Sensitivity analysis of the bridge using the derivative-variance hybrid global sensitivity measure.
fsdAsasfcdlbhρ1hρ2Fα
N = 10000.48510.13210.01740.01144.51 × 10−48.7 × 10−50.16090.00940.00160.00110.19290.1929
N = 1080.48940.13320.01760.01154.60 × 10−48.21 × 10−50.16240.00950.00170.00110.19470.1947
Absolute error0.00430.00110.00020.0001000.00150.00010.000100.00180.0018
Relative error0.88%0.88%0.88%0.89%0.89%0.89%0.88%0.88%0.88%0.88%0.89%0.89%
Table 6. Sensitivity analysis of the bridge using the variance-based global sensitivity measure.
Table 6. Sensitivity analysis of the bridge using the variance-based global sensitivity measure.
fsdAsasfcdlbhρ1hρ2Fα
N = 240,0000.39420.10790.0160.00856.04 × 10−42.24 × 10−50.12850.00840.00080.0010.16490.156
N = 1080.40210.11020.01470.00933.62 × 10−49.18 × 10−50.13230.00790.00140.00090.16020.1606
Absolute error0.00790.00230.00130.00082.42 × 10−46.94 × 10−50.00380.00050.00060.00010.00470.0046
Relative error1.96%2.09%2.09%8.60%66.90%75.60%2.87%6.33%42.80%11.10%2.93%2.86%
Table 7. Comparison between the derivative-based global sensitivity measures, the derivative-based global sensitivity measures, and the variance-based global sensitivity measures.
Table 7. Comparison between the derivative-based global sensitivity measures, the derivative-based global sensitivity measures, and the variance-based global sensitivity measures.
fsdAsasfcdlbhρ1hρ2Fα
Derivative-based global sensitivity measure3.3 × 1088.9 × 1071.2 × 1077.7 × 1073.1 × 1055.9 × 1049.2 × 1076.4 × 1061.1 × 1067.4 × 1051.3 × 1081.3 × 108
Variance-based global sensitivity measure0.39420.10790.0160.00846.04 × 10−42.2 × 10−50.12850.00840.00080.0010.16490.156
Derivative-variance hybrid global sensitivity measure0.48510.13210.01740.01144.51 × 10−48.7 × 10−50.16090.00940.00160.00110.19290.1929
Absolute error0.09090.02420.00140.0031.53 × 10−46.5 × 10−50.03240.0010.00080.00010.0280.0369
Relative error23.10%18.30%8.75%35.70%25.30%295%25.21%11.90%100%10.00%16.90%23.60%
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

Liu, J.; Liu, H.; Zhang, C.; Cao, J.; Xu, A.; Hu, J. Derivative-Variance Hybrid Global Sensitivity Measure with Optimal Sampling Method Selection. Mathematics 2024, 12, 396. https://doi.org/10.3390/math12030396

AMA Style

Liu J, Liu H, Zhang C, Cao J, Xu A, Hu J. Derivative-Variance Hybrid Global Sensitivity Measure with Optimal Sampling Method Selection. Mathematics. 2024; 12(3):396. https://doi.org/10.3390/math12030396

Chicago/Turabian Style

Liu, Jiacheng, Haiyun Liu, Cong Zhang, Jiyin Cao, Aibo Xu, and Jiwei Hu. 2024. "Derivative-Variance Hybrid Global Sensitivity Measure with Optimal Sampling Method Selection" Mathematics 12, no. 3: 396. https://doi.org/10.3390/math12030396

APA Style

Liu, J., Liu, H., Zhang, C., Cao, J., Xu, A., & Hu, J. (2024). Derivative-Variance Hybrid Global Sensitivity Measure with Optimal Sampling Method Selection. Mathematics, 12(3), 396. https://doi.org/10.3390/math12030396

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