Using inverse finite element analysis to identify spinal tissue behaviour in situ

Experimental and Computational Comparison of Intervertebral Disc Bulge for Specimen-Specific Model Evaluation Based on Imaging

Finite element modelling of the spinal unit is a promising preclinical tool to assess the biomechanical outcome of emerging interventions. Currently, most models are calibrated and validated against range of motion and rarely directly against soft-tissue deformation. The aim of this contribution was to develop an in vitro methodology to measure disc bulge and assess the ability of different specimen-specific modelling approaches to predict disc bulge. Bovine bone-disc-bone sections (N = 6) were prepared with 40 glass markers on the intervertebral disc surface. These were initially magnetic resonance (MR)-imaged and then sequentially imaged using peripheral-qCT under axial compression of 1 mm increments. Specimen-specific finite-element models were developed from the CT data, using three different methods to represent the nucleus pulposus geometry with and without complementary use of the MR images. Both calibrated specimen-specific and averaged compressive material properties for the disc tissues were investigated. A successful methodology was developed to quantify the disc bulge in vitro, enabling observation of surface displacement on qCT. From the finite element model results, no clear advantage was found in using geometrical information from the MR images in terms of the models' ability to predict stiffness or disc bulge for bovine intervertebral disc.

Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis

There is an increased interest in studying the biomechanics of the facet joints. For in silico studies, it is therefore important to understand the level of reliability of models for outputs of interest related to the facet joints. In this work, a systematic review of finite element models of multi-level spinal section with facet joints output of interest was performed. The review focused on the methodology used to model the facet joints and its associated validation. From the 110 papers analysed, 18 presented some validation of the facet joints outputs. Validation was done by comparing outputs to literature data, either computational or experimental values; with the major drawback that, when comparing to computational values, the baseline data was rarely validated. Analysis of the modelling methodology showed that there seems to be a compromise made between accuracy of the geometry and nonlinearity of the cartilage behaviour in compression. Most models either used a soft contact representation of the cartilage layer at the joint or included a cartilage layer which was linear elastic. Most concerning, soft contact models usually did not contain much information on the pressure-overclosure law. This review shows that to increase the reliability of in silico model of the spine for facet joints outputs, more needs to be done regarding the description of the methods used to model the facet joints, and the validation for specific outputs of interest needs to be more thorough, with recommendation to systematically share input and output data of validation studies.

Estimating tissue-level properties of porcine talar subchondral bone

Tissue-level properties of bone play an important role when characterising apparent-level bone biomechanical behaviour and yet little is known about its effect at this hierarchical level. In combination with trabecular morphological data these properties can be used to predict bone strength, which becomes an invaluable tool for clinicians in patient treatment planning. This study developed specimen-specific micro-finite element (μFE) models using validated continuum-level models, containing grayscale-derived material properties, to indirectly establish tissue-level properties of porcine talar subchondral bone. Specimen-specific continuum finite element (hFE) models of subchondral trabecular bone were setup using μCT data of ten cylindrical specimens extracted from juvenile porcine tali. The models were validated using quasi-static uniaxial compression testing. Validated hFE models were used to calibrate the tissue modulus of corresponding μFE models by minimising the difference between the μFE and hFE stiffness values. Key trabecular morphological indices (BV/TV, DA, Conn.D, Tb.Th, EF) were evaluated. Good agreement was observed between hFE models and experiment (CCC = 0.66). Calibrated E_tiss was 504 ± 37.65 MPa. Average BV/TV and DA for μFE specimens were 0.37 ± 0.05 and 0.68 ± 0.11, respectively. BV/TV (r² = 0.667) correlated highly with μFE stiffness. The small intra-specimen variation to tissue-level properties suggests that variations to apparent-level stiffness originate from variations to microarchitecture rather than tissue mechanical properties.