  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.
