Hierarchical Bayesian constitutive model selection for high-strain-rate soft material characterization

TL;DR

This paper introduces a hierarchical Bayesian approach combined with microcavitation techniques to accurately select and characterize constitutive models of soft materials under ultra-high-strain-rate conditions, accounting for uncertainty and model complexity.

Contribution

It develops a novel hierarchical Bayesian model selection framework for high-strain-rate soft material characterization using microcavitation data, improving model credibility and parameter estimation.

Findings

Synthetic tests recover true models and parameters.

Experimental data fits credible models well.

Consistent model selection across institutions.

Abstract

The high-fidelity characterization of soft, tissue-like materials under ultra-high-strain-rate conditions is critical in engineering and medicine. Still, it remains challenging due to limited optical access, sensitivity to initial conditions, and experimental variability. Microcavitation techniques (e.g., laser-induced microcavitation) have emerged as a viable method for determining the mechanical properties of soft materials in the ultra-high-strain-rate regime (higher than 10^3 s^{-1}); however, they are limited by measurement noise and uncertainty in parameter estimation. A hierarchical Bayesian model selection method is employed using the Inertial Microcavitation Rheometry (IMR) technique to address these limitations. With this method, the parameter space of different constitutive models is explored to determine the most credible constitutive model that describes laser-induced microcavitation bubble oscillations in soft, viscoelastic, transparent hydrogels. The target data/evidence is computed using a weighted Gaussian likelihood with a hierarchical noise scale, which enables the quantification of uncertainty in model plausibility. Physically informed priors, including range-invariant, stress-based parameter priors, a model-redundancy prior, and a Bayesian Information Criterion motivated model prior, penalize complex models to enforce Occam's razor. Using a precomputed grid of simulations, the probabilistic model selection process enables an initial guess for the Maximum A Posteriori (MAP) material parameter values. Synthetic tests recover the ground-truth models and expected parameters. Using experimental data for gelatin, fibrin, polyacrylamide, and agarose, MAP simulations of credible models reproduce the data. Moreover, a cross-institutional comparison of 10% gelatin indicates consistent constitutive model selection.