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
Bayesian inference and uncertainty propagation using efficient fractional-order viscoelastic models for dielectric elastomers
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
- Award ID(s):
- 1745654
- NSF-PAR ID:
- 10281338
- Date Published:
- Journal Name:
- Journal of Intelligent Material Systems and Structures
- Volume:
- 32
- Issue:
- 4
- ISSN:
- 1045-389X
- Page Range / eLocation ID:
- 486 to 496
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We introduce a data-driven fractional modeling framework for complex materials, and particularly bio-tissues. From multi-step relaxation experiments of distinct anatomical locations of porcine urinary bladder, we identify an anomalous relaxation character, with two power-law-like behaviors for short/long long times, and nonlinearity for strains greater than 25%. The first component of our framework is an existence study, to determine admissible fractional viscoelastic models that qualitatively describe linear relaxation. After the linear viscoelastic model is selected, the second stage adds large-strain effects to the framework through a fractional quasi-linear viscoelastic approach for the nonlinear elastic response of the bio-tissue of interest. From single-step relaxation data of the urinary bladder, a fractional Maxwell model captures both short/long-term behaviors with two fractional orders, being the most suitable model for small strains at the first stage. For the second stage, multi-step relaxation data under large strains were employed to calibrate a four-parameter fractional quasi-linear viscoelastic model, that combines a Scott-Blair relaxation function and an exponential instantaneous stress response, to describe the elastin/collagen phases of bladder rheology. Our obtained results demonstrate that the employed fractional quasi-linear model, with a single fractional order in the range α = 0.25–0.30, is suitable for the porcine urinary bladder, producing errors below 2% without need for recalibration over subsequent applied strains. We conclude that fractional models are attractive tools to capture the bladder tissue behavior under small-to-large strains and multiple time scales, therefore being potential alternatives to describe multiple stages of bladder functionality.more » « less
-
Despite the success of multiscale modeling in science and engineering, embedding molecular-level information into nonlinear reactor design and control optimization problems remains challenging. In this work, we propose a computationally tractable scale-bridging approach that incorporates information from multi-product microkinetic (MK) models with thousands of rates and chemical species into nonlinear reactor design optimization problems. We demonstrate reduced-order kinetic (ROK) modeling approaches for catalytic oligomerization in shale gas processing. We assemble a library of six candidate ROK models based on literature and MK model structure. We find that three metrics—quality of fit (e.g., mean squared logarithmic error), thermodynamic consistency (e.g., low conversion of exothermic reactions at high temperatures), and model identifiability—are all necessary to train and select ROK models. The ROK models that closely mimic the structure of the MK model offer the best compromise to emulate the product distribution. Using the four best ROK models, we optimize the temperature profiles in staged reactors to maximize conversions to heavier oligomerization products. The optimal temperature starts at 630–900K and monotonically decreases to approximately 560 K in the final stage, depending on the choice of ROK model. For all models, staging increases heavier olefin production by 2.5% and there is minimal benefit to more than four stages. The choice of ROK model, i.e., model-form uncertainty, results in a 22% difference in the objective function, which is twice the impact of parametric uncertainty; we demonstrate sequential eigendecomposition of the Fisher information matrix to identify and fix sloppy model parameters, which allows for more reliable estimation of the covariance of the identifiable calibrated model parameters. First-order uncertainty propagation determines this parametric uncertainty induces less than a 10% variability in the reactor optimization objective function. This result highlights the importance of quantifying model-form uncertainty, in addition to parametric uncertainty, in multi-scale reactor and process design and optimization. Moreover, the fast dynamic optimization solution times suggest the ROK strategy is suitable for incorporating molecular information in sequential modular or equation-oriented process simulation and optimization frameworks.more » « less
-
more » « less
SUMMARY Full-waveform inversion (FWI) methods rely on accurate numerical simulation of wave propagation in the analysed medium. Acoustic or elastic wave equations are often used to model seismic wave propagation. These types of simulations do not account for intrinsic attenuation effects due to material anelasticity, and thus correction techniques have been utilized in practice to partially compensate the anelasticity. These techniques often only consider the waveform amplitude correction based on averaging of overall amplitude response over the entire data set, and ignore the phase correction. Viscoelastic wave equations account for the anelastic response in both waveform amplitude and phase, and are therefore a more suitable alternative. In this study, we present a novel 3-D Gauss–Newton viscoelastic FWI (3-D GN-VFWI) method. To address the main challenge of the Gauss–Newton optimization, we develop formulas to compute the Jacobian efficiently by the convolution of virtual sources and backward wavefields. The virtual sources are obtained by directly differentiating the viscoelastic wave equations with respect to model parameters. In order to resolve complex 3-D structures with reasonable computational effort, a homogeneous attenuation (Q factor) is used throughout the analysis to model the anelastic effects. Synthetic and field experiments are performed to demonstrate the utility of the method. The synthetic results clearly demonstrate the ability of the method in characterizing a challenging velocity profile, including voids and reverse velocity layers. The field experimental results show that method successfully characterizes the complex substructure with two voids and undulating limestone bedrock, which are confirmed by invasive tests. Compared to 3-D elastic FWI results, the presented viscoelastic method produces more accurate results regarding depths of the voids and bedrock. This study suggests that the improvement of imaging accuracy would warrant the widespread use of viscoelastic wave equations in FWI problems. To our best knowledge, this is the first reported study on 3-D GN-VFWI at any scale. This study provides the new theory and formulation for the use of Gauss–Newton optimization on the 3-D viscoelastic problem.
-
more » « less
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.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1177/1045389X20969847
