A general first-order global sensitivity analysis method

Fourier amplitude sensitivity test (FAST) is one of the most popular global sensitivity analysis techniques. The main mechanism of FAST is to assign each parameter with a characteristic frequency through a search function. Then, for a specific parameter, the variance contribution can be singled out of the model output by the characteristic frequency. Although FAST has been widely applied, there are two limitations: (1) the aliasing effect among parameters by using integer characteristic frequencies and (2) the suitability for only models with independent parameters. In this paper, we synthesize the improvement to overcome the aliasing effect limitation [Tarantola S, Gatelli D, Mara TA. Random balance designs for the estimation of first order global sensitivity indices. Reliab Eng Syst Safety 2006; 91(6):717–27] and the improvement to overcome the independence limitation [Xu C, Gertner G. Extending a global sensitivity analysis technique to models with correlated parameters. Comput Stat Data Anal 2007, accepted for publication]. In this way, FAST can be a general first-order global sensitivity analysis method for linear/nonlinear models with as many correlated/uncorrelated parameters as the user specifies. We apply the general FAST to four test cases with correlated parameters. The results show that the sensitivity indices derived by the general FAST are in good agreement with the sensitivity indices derived by the correlation ratio method, which is a non-parametric method for models with correlated parameters.     

Full-field sensitivity analysis through dimension reduction and probabilistic surrogate models

Computational mechanics models often are compromised by uncertainty in their governing parameters, especially when the operating environment is incompletely known. Computational sensitivity analysis of a spatially distributed process to its governing parameters therefore is an essential, but often costly, step in uncertainty quantification. A sensitivity analysis method is described which features probabilistic surrogate models developed through equitable sampling of the parameter space, proper orthogonal decomposition (POD) for compact representations of the process’ variability from an ensemble of realizations, and cluster-weighted models of the joint probability density function of each POD coefficient and the governing parameters. Full-field sensitivities, i.e. sensitivities at every point in the computational grid, are computed by analytically differentiating the conditional expected value function of each POD coefficient and projecting the sensitivities onto the POD basis. Statistics of the full-field sensitivities are estimated by sampling the surrogate model throughout the parameter space. Major benefits of this method are: (1) the sensitivities are computed analytically and efficiently from regularized surrogate models, and (2) the conditional variances of the surrogate models may be used to estimate the statistical uncertainty in the sensitivities, which provides a basis for pursuing more data to improve the model. Synthetic examples and a physical example involving near-ground sound propagation through a refracting atmosphere are presented to illustrate the properties of the surrogate models and how full-field sensitivities and their uncertainties are computed.     

Uncertainty-based sensitivity indices for imprecise probability distributions

An uncertainty-based sensitivity index represents the contribution that uncertainty in model input X_i makes to the uncertainty in model output Y. This paper addresses the situation where the uncertainties in the model inputs are expressed as closed convex sets of probability measures, a situation that exists when inputs are expressed as intervals or sets of intervals with no particular distribution specified over the intervals, or as probability distributions with interval-valued parameters. Three different approaches to measuring uncertainty, and hence uncertainty-based sensitivity, are explored. Variance-based sensitivity analysis (VBSA) estimates the contribution that each uncertain input, acting individually or in combination, makes to variance in the model output. The partial expected value of perfect information (partial EVPI), quantifies the (financial) value of learning the true numeric value of an input. For both of these sensitivity indices the generalization to closed convex sets of probability measures yields lower and upper sensitivity indices. Finally, the use of relative entropy as an uncertainty-based sensitivity index is introduced and extended to the imprecise setting, drawing upon recent work on entropy measures for imprecise information.     

A new computational method of a moment-independent uncertainty importance measure

For a risk assessment model, the uncertainty in input parameters is propagated through the model and leads to the uncertainty in the model output. The study of how the uncertainty in the output of a model can be apportioned to the uncertainty in the model inputs is the job of sensitivity analysis. Saltelli [Sensitivity analysis for importance assessment. Risk Analysis 2002;22(3):579-90] pointed out that a good sensitivity indicator should be global, quantitative and model free. Borgonovo [A new uncertainty importance measure. Reliability Engineering and System Safety 2007;92(6):771-84] further extended these three requirements by adding the fourth feature, moment-independence, and proposed a new sensitivity measure, δ_i. It evaluates the influence of the input uncertainty on the entire output distribution without reference to any specific moment of the model output. In this paper, a new computational method of δ_i is proposed. It is conceptually simple and easier to implement. The feasibility of this new method is proved by applying it to two examples.     

Principal components analysis in sensitivity studies of dynamic systems

Local sensitivity analysis is a modeling tool for determining the effects of single parameter variations on the output of a first order differential system. To determine the effects of multi-parameter variations, the local sensitivity matrix can be used in a first order Taylor series to approximately model the variance of the system output. For dynamic systems, a principal component analysis based on this time varying approximation reveals the evolution of the directions and magnitudes of greatest variation of the system output derived from input variability in the parameters. Such an analysis acts as a means of modeling the robustness of dynamic differential systems.     

Sensitivity analysis of a complex, proposed geologic waste disposal system using the Fourier Amplitude Sensitivity Test method

The Fourier Amplitude Sensitivity Test (FAST) method has been used to perform a sensitivity analysis of a computer model developed for conducting total system performance assessment of the proposed high-level nuclear waste repository at Yucca Mountain, Nevada, USA. The computer model has a large number of random input parameters with assigned probability density functions, which may or may not be uniform, for representing data uncertainty. The FAST method, which was previously applied to models with parameters represented by the uniform probability distribution function only, has been modified to be applied to models with nonuniform probability distribution functions. Using an example problem with a small input parameter set, several aspects of the FAST method, such as the effects of integer frequency sets and random phase shifts in the functional transformations, and the number of discrete sampling points (equivalent to the number of model executions) on the ranking of the input parameters have been investigated. Because the number of input parameters of the computer model under investigation is too large to be handled by the FAST method, less important input parameters were first screened out using the Morris method. The FAST method was then used to rank the remaining parameters. The validity of the parameter ranking by the FAST method was verified using the conditional complementary cumulative distribution function (CCDF) of the output. The CCDF results revealed that the introduction of random phase shifts into the functional transformations, proposed by previous investigators to disrupt the repetitiveness of search curves, does not necessarily improve the sensitivity analysis results because it destroys the orthogonality of the trigonometric functions, which is required for Fourier analysis.     

Multivariate sensitivity analysis to measure global contribution of input factors in dynamic models

Many dynamic models are used for risk assessment and decision support in ecology and crop science. Such models generate time-dependent model predictions, with time either discretised or continuous. Their global sensitivity analysis is usually applied separately on each time output, but Campbell et al. (2006 [1]) advocated global sensitivity analyses on the expansion of the dynamics in a well-chosen functional basis. This paper focuses on the particular case when principal components analysis is combined with analysis of variance. In addition to the indices associated with the principal components, generalised sensitivity indices are proposed to synthesize the influence of each parameter on the whole time series output. Index definitions are given when the uncertainty on the input factors is either discrete or continuous and when the dynamic model is either discrete or functional. A general estimation algorithm is proposed, based on classical methods of global sensitivity analysis.The method is applied to a dynamic wheat crop model with 13 uncertain parameters. Three methods of global sensitivity analysis are compared: the Sobol'-Saltelli method, the extended FAST method, and the fractional factorial design of resolution 6.     

模糊分布参数的全局灵敏度分析新方法

为合理度量随机输入变量分布参数的模糊性对输出性能统计特征的影响,提出了模糊分布参数的全局灵敏度效应指标,并研究了指标的高效求解方法。首先,分析了不确定性从模糊分布参数至模型输出响应统计特征的传递机理,以输出性能期望响应为例,利用输出均值的无条件隶属函数与给定模糊分布参数取值条件下的隶属函数的平均差异来度量模糊分布参数的影响,建立了模糊分布参数的全局灵敏度效应指标。其次,为减少所提指标的计算成本、提高计算效率,采用了扩展蒙特卡罗模拟法(EMCS)来估算输入变量分布参数与模型输出响应统计特征的函数关系。最后通过对算例的计算,验证该文所提方法的准确性和高效性。     

CDF sensitivity analysis technique for ranking influential parameters in the performance assessment of the proposed high-level waste repository at Yucca Mountain, Nevada, USA

A cumulative distribution function (CDF)-based method has been used to perform sensitivity analysis on a computer model that conducts total system performance assessment of the proposed high-level nuclear waste repository at Yucca Mountain, and to identify the most influential input parameters affecting the output of the model. The performance assessment computer model referred to as the TPA code, was recently developed by the US nuclear regulatory commission (NRC) and the center for nuclear waste regulatory analyses (CNWRA), to evaluate the performance assessments conducted by the US department of energy (DOE) in support of their license application. The model uses a probabilistic framework implemented through Monte Carlo or Latin hypercube sampling (LHS) to permit the propagation of uncertainties associated with model parameters, conceptual models, and future system states. The problem involves more than 246 uncertain parameters (also referred to as random variables) of which the ones that have significant influence on the response or the uncertainty of the response must be identified and ranked. The CDF-based approach identifies and ranks important parameters based on the sensitivity of the response CDF to the input parameter distributions. Based on a reliability sensitivity concept [AIAA Journal 32 (1994) 1717], the response CDF is defined as the integral of the joint probability-density-function of the input parameters, with a domain of integration that is defined by a subset of the samples. The sensitivity analysis does not require explicit knowledge of any specific relationship between the response and the input parameters, and the sensitivity is dependent upon the magnitude of the response. The method allows for calculating sensitivity over a wide range of the response and is not limited to the mean value.     

Analysis of the sensitivity properties of a model of vector-borne bubonic plague.

Model sensitivity is a key to evaluation of mathematical models in ecology and evolution, especially in complex models with numerous parameters. In this paper, we use some recently developed methods for sensitivity analysis to study the parameter sensitivity of a model of vector-borne bubonic plague in a rodent population proposed by Keeling & Gilligan. The new sensitivity tools are based on a variational analysis involving the adjoint equation. The new approach provides a relatively inexpensive way to obtain derivative information about model output with respect to parameters. We use this approach to determine the sensitivity of a quantity of interest (the force of infection from rats and their fleas to humans) to various model parameters, determine a region over which linearization at a specific parameter reference point is valid, develop a global picture of the output surface, and search for maxima and minima in a given region in the parameter space.     

A new uncertainty importance measure

Uncertainty in parameters is present in many risk assessment problems and leads to uncertainty in model predictions. In this work, we introduce a global sensitivity indicator which looks at the influence of input uncertainty on the entire output distribution without reference to a specific moment of the output (moment independence) and which can be defined also in the presence of correlations among the parameters. We discuss its mathematical properties and highlight the differences between the present indicator, variance-based uncertainty importance measures and a moment independent sensitivity indicator previously introduced in the literature. Numerical results are discussed with application to the probabilistic risk assessment model on which Iman [A matrix-based approach to uncertainty and sensitivity analysis for fault trees. Risk Anal 1987;7(1):22–33] first introduced uncertainty importance measures.     

Sensitivity analysis in presence of model uncertainty and correlated inputs

The first motivation of this work is to take into account model uncertainty in sensitivity analysis (SA). We present with some examples, a methodology to treat uncertainty due to a mutation of the studied model. Development of this methodology has highlighted an important problem, frequently encountered in SA: how to interpret sensitivity indices when random inputs are non-independent? This paper suggests a strategy for the problem of SA of models with non-independent random inputs. We propose a new application of the multidimensional generalization of classical sensitivity indices, resulting from group sensitivities (sensitivity of the output of the model to a group of inputs), and describe an estimation method based on Monte-Carlo simulations. Practical and theoretical applications illustrate the interest of this method.     

A fatigue model with local sensitivity analysis

The goal of this paper is two fold. First, it introduces a general parametric lifetime model for high‐cycle fatigue regime derived from physical, statistical, engineering and dimensional analysis considerations. The proposed model has two threshold parameters and three Weibull distribution parameters. A two‐step procedure is presented to estimate the parameters. In the first step, the two threshold parameters are estimated by minimizing a least squares regression function. In the second step, the parameters are estimated by the maximum likelihood method after pooling together the data from different stress levels. Since parameter estimation should always be accompanied by a sensitivity analysis of the fitted model, the second goal of this paper is to propose a method for sensitivity analysis for fatigue models. We show that the proposed sensitivity analysis methods are general and can be applied to any fatigue or lifetime model, not just to the one proposed here. Although several fatigue models have been proposed in the literature, to our knowledge this is the first attempt to produce methods for sensitivity analysis for fatigue models. The proposed method makes use of the well‐known duality property of mathematical programming, which states that the partial derivatives of the primal objective function with respect to the constraints right hand side parameters are the optimal values of the negative of the dual problem variables. For the parameters or data, for which sensitivities are sought, to appear on the right hand side, they are converted into artificial variables and set to their actual values, thus obtaining the desired constraints. Both the estimation and sensitivity analysis methods are illustrated by two examples, one application using real fatigue data and the other using simulated data. In addition, the sensitivity proposed method is also applied to an alternative fatigue model. Finally, some specific conclusions and recommendations are also given.     

Relative importance of uncertain structural parameters. Part I: algorithm

A novel procedure for estimating the relative importance of uncertain parameters of complex FE model is presented. The method is specifically directed toward problems involving high-dimensional input parameter spaces, as they are encountered during uncertainty analysis of large scale, refined FE models. In these cases one is commonly faced with thousands of uncertain parameters and traditional techniques, e.g. finite difference or direct differentiation methods become expensive. In contrast, the presented method quickly filters out the most influential variables. Hence, the main objective is not to compute the sensitivity but to identify those parameters whose random variations have the biggest influence on the response. This is achieved by generating a set of samples with direct Monte Carlo simulation, which are closely scattered around the point at which the relative importance measures are sought. From these samples, estimators of the relative importance are synthesized and the most important ones are refined with a method of choice. In this paper, the underlying theory as well as the resulting algorithm is presented.     

Bayesian analysis for estimating statistical parameter distributions of elasto-viscoplastic material models

High temperature design methods rely on constitutive models for inelastic deformation and failure typically calibrated against the mean of experimental data without considering the associated scatter. Variability may arise from the experimental data acquisition process, from heat-to-heat material property variations, or both and need to be accurately captured to predict parameter bounds leading to efficient component design. Applying the Bayesian Markov Chain Monte Carlo (MCMC) method to produce statistical models capturing the underlying uncertainty in the experimental data is an area of ongoing research interest. This work varies aspects of the Bayesian MCMC method and explores their effect on the posterior parameter distributions for a uniaxial elasto-viscoplastic damage model using synthetically generated reference data. From our analysis with the uniaxial inelastic model we determine that an informed prior distribution including different types of test conditions results in more accurate posterior parameter distributions. The parameter posterior distributions, however, do not improve when increasing the number of similar experimental data. Additionally, changing the amount of scatter in the data affects the quality of the posterior distributions, especially for the less sensitive model parameters. Moreover, we perform a sensitivity study of the model parameters against the likelihood function prior to the Bayesian analysis. The results of the sensitivity analysis help to determine the reliability of the posterior distributions and reduce the dimensionality of the problem by fixing the insensitive parameters. The comprehensive study described in this work demonstrates how to efficiently apply the Bayesian MCMC methodology to capture parameter uncertainties in high temperature inelastic material models. Quantifying these uncertainties in inelastic models will improve high temperature engineering design practices and lead to safer, more effective component designs.     

考虑参数不确定性的无筋砌体结构地震易损性分析

鉴于无筋砌体结构参数的高离散性及其地震响应的强非线性,在地震易损性分析中考虑其结构参数的不确定性显得尤为必要。以4幢不同层数的无筋砌体结构为研究对象,在OpenSees中建立等效框架有限元模型,采用增量动力法和一次二阶矩法考察了地震动及结构参数不确定性对其地震易损性的影响。分析结果表明:在无筋砌体结构的地震易损性分析中需考虑地震动和结构参数的不确定性影响,且结构地震破坏程度越高这种影响程度越大;结构参数的不确定性影响程度与地震动的影响程度相当;结构参数不确定性相比于地震动不确定性的影响随着结构层数的降低而变得更为明显;敏感性分析表明,结构阻尼比变化对易损性曲线地震强度中位值的影响可达到其他单个参数的4倍,其敏感性最高。     

Case studies in Gaussian process modelling of computer codes

In this paper we present a number of recent applications in which an emulator of a computer code is created using a Gaussian process model. Tools are then applied to the emulator to perform sensitivity analysis and uncertainty analysis. Sensitivity analysis is used both as an aid to model improvement and as a guide to how much the output uncertainty might be reduced by learning about specific inputs. Uncertainty analysis allows us to reflect output uncertainty due to unknown input parameters, when the finished code is used for prediction.The computer codes themselves are currently being developed within the UK Centre for Terrestrial Carbon Dynamics.     

Probabilistic analysis of truss structures with uncertain parameters (virtual distortion method approach)

A new approach for probabilistic characterization of linear elastic redundant trusses with uncertainty on the various members subjected to deterministic loads acting on the nodes of the structure is presented. The method is based on the simple observation that variations of structural parameters are equivalent to superimposed strains on a reference structure depending on the axial forces on the elastic modulus of the original structure as well as on the uncertainty (virtual distortion method approach). Superposition principle may be applied to separate contribution to mechanical response due to external loads and parameter variations. Statically determinate trusses dealt with the proposed method yields explicit analytic solution in terms of displacements while redundant trusses have been studied by means of an asymptotic expansion exhibiting explicit dependence on parameter fluctuations. Probabilistic characterization of the response may then be obtained both for statically determinate and statically indeterminate stochastic trusses.     

Uncertainty quantification for multiscale modeling of polymer nanocomposites with correlated parameters

We propose a stochastic multiscale method to quantify the correlated key-input parameters influencing the mechanical properties of polymer nanocomposites (PNCs). The variations of parameters at nano-, micro-, meso- and macro-scales are connected by a hierarchical multiscale approach. The first-order and total-effect sensitivity indices are determined first. The input parameters include the single-walled carbon nanotube (SWNT) length, the SWNT waviness, the agglomeration and volume fraction of SWNTs. Stochastic methods consistently predict that the key parameters for the Young’s modulus of the composite are the volume fraction followed by the averaged longitudinal modulus of equivalent fiber (EF), the SWNT length, and the averaged transverse modulus of the EF, respectively. The averaged longitudinal modulus of the EF is estimated to be the most important parameter with respect to the Poisson’s ratio followed by the volume fraction, the SWNT length, and the averaged transverse modulus of the EF, respectively. On the other hand, the agglomeration parameters have insignificant effect on both Young’s modulus and Poisson’s ratio compared to other parameters. The sensitivity analysis (SA) also reveals the correlation between the input parameters and its effect on the mechanical properties.