A Bayesian constitutive model selection framework for biaxial mechanical testing of planar soft tissues: application to porcine aortic valves

TL;DR

This paper introduces a Bayesian framework for selecting the most appropriate constitutive models for biaxial testing data of soft tissues, demonstrated on porcine aortic valves, accounting for sample variability and tissue differences.

Contribution

It develops a Bayesian model selection method that evaluates the probability of different constitutive models based on experimental data, incorporating tissue variability and PCA analysis.

Findings

May--Newman model is most probable for porcine aortic valve data

Different cusp types show different most probable models

PCA captures significant variations in tissue mechanical properties

Abstract

A variety of constitutive models have been developed for soft tissue mechanics. However, there is no established criterion to select a suitable model for a specific application. Although the model that best fits the experimental data can be deemed the most suitable model, this practice often can be insufficient given the inter-sample variability of experimental observations. Herein, we present a Bayesian approach to calculate the relative probabilities of constitutive models based on biaxial mechanical testing of tissue samples. 46 samples of porcine aortic valve tissue were tested using a biaxial stretching setup. For each sample, seven ratios of stresses along and perpendicular to the fiber direction were applied. The probabilities of eight invariant-based constitutive models were calculated based on the experimental data using the proposed model selection framework. The calculated probabilities showed that, out of the considered models and based on the information available through the utilized experimental dataset, the May--Newman model was the most probable model for the porcine aortic valve data. When the samples were grouped into different cusp types, the May--Newman model remained the most probable for the left- and right-coronary cusps, whereas for non-coronary cusps two models were found to be equally probable: the Lee--Sacks model and the May--Newman model. This difference between cusp types was found to be associated with the first principal component analysis (PCA) mode, where this mode's amplitudes of the non-coronary and right-coronary cusps were found to be significantly different. Our results show that a PCA-based statistical model can capture significant variations in the mechanical properties of soft tissues. The presented framework is applicable to any tissue type, and has the potential to provide a structured and rational way of making simulations population-based.