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

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³ 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.