Global and local identifiability analysis of a nonlinear biphasic constitutive model in confined compression

Application of biomechanical models relies on model parameters estimated from experimental data. Parameter non-identifiability, when the same model output can be produced by many sets of parameter values, introduces severe errors yet has received relatively little attention in biomechanics and is subtle enough to remain unnoticed in the absence of deliberate verification. The present work develops a global identifiability analysis method in which cluster analysis and singular value decomposition are applied to vectors of parameter-output variable correlation coefficients. This method provides a visual representation of which specific experimental design elements are beneficial or harmful in terms of parameter identifiability, supporting the correction of deficiencies in the test protocol prior to testing physical specimens. The method was applied to a representative nonlinear biphasic model for cartilaginous tissue, demonstrating that confined compression data does not provide identifiability for the biphasic model parameters. This result was confirmed by two independent analyses: local analysis of the Hessian of a sum-of-squares error cost function and observation of the behaviour of two optimization algorithms. Therefore, confined compression data are insufficient for the calibration of general-purpose biphasic models. Identifiability analysis by these or other methods is strongly recommended when planning future experiments.