skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Global sensitivity analysis of fractional-order viscoelasticity model
In this paper, we investigate hyperelastic and viscoelastic model parameters using Global Sensitivity Analysis(GSA). These models are used to characterize the physical response of many soft-elastomers, which are used in a wide variety of smart material applications. Recent research has shown the effectiveness of using fractional-order calculus operators in modeling the viscoelastic response. The GSA is performed using parameter subset selection (PSS), which quantifies the relative parameter contributions to the linear and nonlinear, fractional-order viscoelastic models. Calibration has been performed to quantify the model parameter uncertainty; however, this analysis has led to questions regarding parameter sensitivity and whether or not the parameters can be uniquely identified given the available data. By performing GSA we can determine which parameters are most influential in the model, and fix non-influential parameters at a nominal value. The model calibration can then be performed to quantify the uncertainty of the influential parameters.  more » « less
Award ID(s):
1745654
PAR ID:
10105020
Author(s) / Creator(s):
Date Published:
Journal Name:
SPIE Smart Structures + Nondestructive Evaluation
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Dielectric elastomers are employed for a wide variety of adaptive structures. Many of these soft elastomers exhibit significant rate-dependencies in their response. Accurately quantifying this viscoelastic behavior is non-trivial and in many cases a nonlinear modeling framework is required. Fractional-order operators have been applied to modeling viscoelastic behavior for many years, and recent research has shown fractional-order methods to be effective for nonlinear frameworks. This implementation can become computationally expensive to achieve an accurate approximation of the fractional-order derivative. Accurate estimation of the elastomer’s viscoelastic behavior to quantify parameter uncertainty motivates the use of Markov Chain Monte Carlo (MCMC) methods. Since MCMC is a sampling based method, requiring many model evaluations, efficient estimation of the fractional derivative operator is crucial. In this paper, we demonstrate the effectiveness of using quadrature techniques to approximate the Riemann–Liouville definition for fractional derivatives in the context of estimating the uncertainty of a nonlinear viscoelastic model. We also demonstrate the use of parameter subset selection techniques to isolate parameters that are identifiable in the sense that they are uniquely determined by measured data. For those identifiable parameters, we employ Bayesian inference to compute posterior distributions for parameters. Finally, we propagate parameter uncertainties through the models to compute prediction intervals for quantities of interest. 
    more » « less
  2. Abstract Calibration of agent‐based models (ABMs) is a major challenge due to the complex nature of the systems being modeled, the heterogeneous nature of geographical regions, the varying effects of model inputs on the outputs, and computational intensity. Nevertheless, ABMs need to be carefully tuned to achieve the desirable goal of simulating spatiotemporal phenomena of interest, and a well‐calibrated model is expected to achieve an improved understanding of the phenomena. To address some of the above challenges, this article proposes an integrated framework of global sensitivity analysis (GSA) and calibration, called GSA‐CAL. Specifically, variance‐based GSA is applied to identify input parameters with less influence on differences between simulated outputs and observations. By dropping these less influential input parameters in the calibration process, this research reduces the computational intensity of calibration. Since GSA requires many simulation runs, due to ABMs' stochasticity, we leverage the high‐performance computing power provided by the advanced cyberinfrastructure. A spatially explicit ABM of influenza transmission is used as the case study to demonstrate the utility of the framework. Leveraging GSA, we were able to exclude less influential parameters in the model calibration process and demonstrate the importance of revising local settings for an epidemic pattern in an outbreak. 
    more » « less
  3. We perform parameter subset selection and uncertainty analysis for phase-field models that are applied to the ferroelectric material lead titanate. A motivating objective is to determine which parameters are influential in the sense that their uncertainties directly affect the uncertainty in the model response, and fix noninfluential parameters at nominal values for subsequent uncertainty propagation. We employ Bayesian inference to quantify the uncertainties of gradient exchange parameters governing 180° and 90° tetragonal phase domain wall energies. The uncertainties of influential parameters determined by parameter subset selection are then propagated through the models to obtain credible intervals when estimating energy densities quantifying polarization and strain across domain walls. The results illustrate various properties of Landau and electromechanical coupling parameters and their influence on domain wall interactions. We employ energy statistics, which quantify distances between statistical observations, to compare credible intervals constructed using a complete set of parameters against an influential subset of parameters. These intervals are obtained from the uncertainty propagation of the model input parameters on the domain wall energy densities. The investigation provides critical insight into the development of parameter subset selection, uncertainty quantification, and propagation methodologies for material modeling domain wall structure evolution, informed by density functional theory simulations. 
    more » « less
  4. Abstract Developing constitutive models that can describe a complex fluid’s response to an applied stimulus has been one of the critical pursuits of rheologists. The complexity of the models typically goes hand-in-hand with that of the observed behaviors and can quickly become prohibitive depending on the choice of materials and/or flow protocols. Therefore, reducing the number of fitting parameters by seeking compact representations of those constitutive models can obviate extra experimentation to confine the parameter space. To this end, fractional derivatives in which the differential response of matter accepts non-integer orders have shown promise. Here, we develop neural networks that are informed by a series of different fractional constitutive models. These fractional rheology-informed neural networks (RhINNs) are then used to recover the relevant parameters (fractional derivative orders) of three fractional viscoelastic constitutive models, i.e., fractional Maxwell, Kelvin-Voigt, and Zener models. We find that for all three studied models, RhINNs recover the observed behavior accurately, although in some cases, the fractional derivative order is recovered with significant deviations from what is known as ground truth. This suggests that extra fractional elements are redundant when the material response is relatively simple. Therefore, choosing a fractional constitutive model for a given material response is contingent upon the response complexity, as fractional elements embody a wide range of transient material behaviors. 
    more » « less
  5. Abstract. Spatially distributed hydrological models are commonly employed to optimize the locations of engineering control measures across a watershed. Yet, parameter screening exercises that aim to reduce the dimensionality of the calibration search space are typically completed only for gauged locations, like the watershed outlet, and use screening metrics that are relevant to calibration instead of explicitly describing the engineering decision objectives. Identifying parameters that describe physical processes in ungauged locations that affect decision objectives should lead to a better understanding of control measure effectiveness. This paper provides guidance on evaluating model parameter uncertainty at the spatial scales and flow magnitudes of interest for such decision-making problems. We use global sensitivity analysis to screen parameters for model calibration, and to subsequently evaluate the appropriateness of using multipliers to adjust the values of spatially distributed parameters to further reduce dimensionality. We evaluate six sensitivity metrics, four of which align with decision objectives and two of which consider model residual error that would be considered in spatial optimizations of engineering designs. We compare the resulting parameter selection for the basin outlet and each hillslope. We also compare basin outlet results for four calibration-relevant metrics. These methods were applied to a RHESSys ecohydrological model of an exurban forested watershed near Baltimore, MD, USA. Results show that (1) the set of parameters selected by calibration-relevant metrics does not include parameters that control decision-relevant high and low streamflows, (2) evaluating sensitivity metrics at the basin outlet misses many parameters that control streamflows in hillslopes, and (3) for some multipliers, calibrating all parameters in the set being adjusted may be preferable to using the multiplier if parameter sensitivities are significantly different, while for others, calibrating a subset of the parameters may be preferable if they are not all influential. Thus, we recommend that parameter screening exercises use decision-relevant metrics that are evaluated at the spatial scales appropriate to decision making. While including more parameters in calibration will exacerbate equifinality, the resulting parametric uncertainty should be important to consider in discovering control measures that are robust to it. 
    more » « less