Understanding lumbar facet joint involvement and biomechanical changes post spinal fusion is limited. This study aimed to establish an in vitro model assessing mechanical effects of fusion on human lumbar facet joints, employing synchronized motion, pressure, and stiffness analysis. Seven human lumbar spinal units (age 54 to 92, ethics 15/YH/0096) underwent fusion via a partial nucleotomy model mimicking a lateral cage approach with PMMA cement injection. Mechanical testing pre and post-fusion included measuring compressive displacement and load, local motion capture, and pressure mapping at the facet joints. pQCT imaging (82 microns isotropic) was carried out at each stage to assess the integrity of the vertebral endplates and quantify the amount of cement injected. Before fusion, relative facet joint displacement (6.5 ± 4.1 mm) at maximum load (1.1 kN) exceeded crosshead displacement (3.9 ± 1.5 mm), with loads transferred across both facet joints. After fusion, facet displacement (2.0 ± 1.2 mm) reduced compared to pre-fusion, as was the crosshead displacement (2.2 ± 0.6 mm). Post-fusion loads (71.4 ± 73.2 N) transferred were reduced compared to pre-fusion levels (194.5 ± 125.4 N). Analysis of CT images showed no endplate damage post-fusion, whilst the IVD tissue: cement volume ratio did not correlate with the post-fusion behaviour of the specimens.Objectives
Methods and Results
While spinal fusion is known to be associated with adjacent disc degeneration, little is known on the role of the facet joints in the process, and whether their altered biomechanics following fusion plays a role in further spinal degeneration. This work aimed to develop a model and method to sequentially measure the effects of spinal fusion on lumbar facet joints through synchronisation of both motion analysis, pressure mapping and mechanical analysis. Parallel measurements of mature ovine lumbar facet joints (∼8yr old, n=3) were carried out using synchronised load and displacement measurements, motion capture during loading and pressure mapping of the joint spaces during loading. Functional units were prepared and cemented in PMMA endcaps. Displacement-controlled compression measurements were carried out using a materials testing machine (3365, Instron, USA) at 1 mm/min up to 950 N with the samples in a neutral position, while motion capture of the facet joints during compression was carried out using orthogonal HD webcams (Logitech, Switzerland) to measure the displacement of key facet joint features. The pressure mapping of load transfer during displacement was carried out using a flexible pressure sensor (6900 series, Tekscan, USA). Each sample was imaged at an isotropic resolution of 82 microns using a μCT scanner (XtremeCT, Scanco, Switzerland) to quantify the curvature within the facet joints.Abstract
Objectives
Methods
Autologous osteochondral grafting has demonstrated positive outcomes for treating articular cartilage defects by replacing the damaged region with a cylindrical graft consisting of bone with a layer of cartilage, taken from a non-loadbearing region of the knee. Despite positive clinical use, factors that cause graft subsidence or poor integration are relatively unknown. The aim of this study was to develop finite element (FE) models of osteochondral grafts within a tibiofemoral joint and to investigate parameters affecting osteochondral graft stability. Initial experimental tests on cadaveric femurs were performed to calibrate the bone properties and graft-bone frictional forces for use in corresponding FE models, generated from µCT scan data. The effects of cartilage defects and osteochondral graft repair were measured by examining contact pressure changes using in vitro tests on a single cadaveric human tibiofemoral joint. Six defects were created in the femoral condyles which were subsequently treated with osteochondral autografts or metal pins. Matching µCT scan-based FE models were created, and the contact patches were compared. Sensitivity to graft bone properties was investigated. The bone material properties and graft-bone frictional forces were successfully calibrated from the initial tests with good resulting levels of agreement (CCC=0.87). The tibiofemoral joint experiment provided a range of cases to model. These cases were well captured experimentally and represented accurately in the FE models. Graft properties relative to host bone had large effects on immediate graft stability despite limited changes to resultant cartilage contact pressure. Model confidence was built through extensive validation and sensitivity testing, and demonstrated that specimen-specific properties were required to accurately represent graft behaviour. The results indicate that graft bone properties affect the immediate stability, which is important for the selection of allografts and design of future synthetic grafts. Supported by the EPSRC-EP/P001076.Acknowledgements
Osteoarthritis (OA) is one of the lead causes of pain and disability in adults. Bone marrow lesions (BMLs) are one feature of subchondral bone involvement in OA. MRI images suggest changes in tissue content and properties in the affected regions however, it is not known if this alters the mechanical behavior of the bone, which could in turn affect OA progression. The aim of this study was to characterize the mechanical properties of BMLs, using a combined experimental and computational approach. Six human cadaveric patellae from donors aged 56–76 were used in this study; all exhibited BML regions under MRI. Bone plugs were taken from non-BML (n = 6) and BML (n = 7) regions within the patellae, with guidance from the MRI. The plugs were imaged at 82µm resolution using micro computed tomography (µCT) and tested under uniaxial compression. Finite element (FE) models were created for each plug from the µCT scans and morphological properties such as bone volume fraction (BV/TV) were also determined. The relationship between bone volume fraction and apparent modulus was investigated for both sample groups.Abstract
Introduction
Methods
Osteoarthritis (OA) affects more than four million people in the UK alone. Bone marrow lesions (BMLs) are a common feature of subchondral bone pathology in OA. Both bone volume fraction and mineral density within the BML are abnormal. The aim of this study was to investigate the effect of a potential treatment (bone augmentation) for BMLs on the knee joint mechanics in cases with healthy and fully degenerated cartilage, using finite element (FE) models of the joint to study the effect of BML size. FE models of a human tibiofemoral joint were created based on models from the Open Knee project (simtk.org). Following initial mesh convergence studies, each model was manipulated in ScanIP (Synopsys-Simpleware, UK) to incorporate a BML 2mm below the surface of the tibial contact region. Models representing extreme cases (healthy cartilage, no cartilage; BML region as an empty cavity or filled with bone substitution material (200GPa)) were generated, each with different sizes of BML. Models were tested under a representative physiological load of 2kN.Abstract
Introduction
Methods
Flattening of the talar dome is observed clinically in haemarthropathy as structural and functional changes advance but has not been quantified yet. In order to confirm clinical observation, and assess the degree of change, morphological measurements were derived from MR images. Four measurements were taken, using ImageJ (1.52v), from sagittal MRI projections at three locations – medial, lateral and central: Trochlear Tali Arc Length (TaAL), Talar Height (TaH), Trochlear Tali Length (TaL), and Trochlear Tali Radius (TaR). These measurements were used to generate three ratios of interest: TaR:TaAL, TaAL:TaL, and TaL:TaH. With the hypothesis of a flattening of the talar dome with haemarthropathy, it was expected that TaR:TaAL and TaL:TaH should be greater for haemophilic ankles, and TaAL:TaL should be smaller. A total of 126 MR images (ethics: MEEC 18–022) were included to assess the difference in those ratios between non-diseased ankles (33 images from 11 volunteers) and haemophilic ankles (93 images from 8 patients’ ankles). Non-diseased control measurements were compared to literature to assess the capacity of doing measurements on MRI instead of radiographs or CT.Abstract
OBJECTIVE
METHODS
Develop a methodology to assess the long term mechanical behavior of intervertebral discs by utilizing novel sequential state testing. Bovine functional spinal units were sequentially mechanically tested in (1) native (n=8), (2) degenerated (n=4), and (3) treated states (n=4). At stage (2), artificial degeneration was created using rapid enzymatic degeneration, followed by a 24 hour hold period under static load at 42°C. At stage (3), nucleus augmentation treatments were injected with a hydrogel or a ‘sham’ (water, chondroitin sulfate) injection. The mechanical protocol employed applied a static load hold period followed by cyclic compressive loading between ∼350 and 750 N at 1 Hz. 1000 cycles were applied at each stage, and the final test on each specimen was extended up to 20000 cycles. To verify if test time can be reduced, functions were fitted using stiffness data up to 100, 1000, 2500, 5000, 10000 and 20000 cycles. Linear regression for the native specimens comparing the stiffness at various cycles to the stiffness at 20000 cycles was completed.Abstract
Objectives
Methods
Osteochondral grafting (OCG) is one treatment strategy for osteoarthritis with good clinical results. Decellularised tissues provide a promising alternative to standard autografts or allografts. This study aimed to compare the stability of traditional OCG and decellularised scaffolds upon initial implantation. Host cubes (N=16) were extracted from porcine femoral condyles around an artificial defect hole. Grafts (N=11) were harvested from the trochlear groove; porcine decellularized osteochondral scaffold (N=5) were prepared. Each host was secured in fixtures and submerged in PBS at 37 ºC. Each graft or scaffold was press fit into one of the hosts, then pushed in for 5 mm, using an indenter (Instron3365) and pushed out in the opposite direction for 10 mm. Parameters analysed were the force required to initiate movement (Dislodging Force) and the maximum force (Max Force).Abstract
Objectives
Methods
Over 1% of the global population suffers with ankle osteoarthritis (OA), yet there is limited knowledge on the changes to subchondral bone with OA. In other joints, it has been shown that bone becomes osteosclerotic, with fewer, thicker trabeculae that become hypomineralised, causing an increased apparent bone volume fraction (BV/TV). Microstructural alterations reduce overall joint strength, which may impact the success of late-stage surgical interventions, such as total ankle arthroplasty (TAA). Previous ankle studies have evaluated changes to cartilage, bone plate and bone morphology with OA, hence this study aimed to characterise changes to trabecular architecture. Three ankle joints were isolated from non-diseased cadaveric feet (three males: 43, 50 and 57 years, MEEC 18-027). Cylindrical subchondral bone specimens (N=6, 6.5 mm Ø) were extracted from the tibial plafond. Osteoarthritic bone samples (N=6, distal tibia) were sourced from local patients (three males: 65, 58 and 68 years, NREC 07/Q1205/27) undergoing TAA surgery. Specimens were imaged using µCT at a 16 µm isotropic resolution (µCT-100 ScanCo Medical). Virtual cores of bone (6.5 mm Ø) were extracted from the image data of the osteoarthritic specimens and trimmed to a height of 4 mm. BoneJ was used to evaluate key morphological indices: BV/TV; anisotropy (DA); trabecular thickness (Tb.Th); trabecular density (Conn.D) and ellipsoid factor (EF) which characterises rod/plate geometry. Differences between the two groups of specimens were evaluated using a t-test with Bonferroni correction.Abstract
Objectives
Methods
Osteoarthritis (OA) affects over 8.75 million people in the UK creating the need for early stage interventions. Osteochondral (OC) grafting has been used to repair full thickness lesions but the efficacy of this therapy is questionable. As a first step in developing a testing framework able to predict the potential and suitability of OC grafts for repair, here, we present experimental data to be used in informing boundary conditions, input parameters and testing sequences for developing and verifying an FE model of the interaction of OC grafts and surrounding host tissue. Ten OC cylindrical grafts (height: 10mm; diameter: n=5–6.5mm; n=5–8.5mm) were harvested from adult porcine femurs (60–70kg). Unconfined compression tests were conducted using an Instron3365 with a 500N load cell and a BioBath filled with PBS at 37ºC. The OC grafts (prior to separation of cartilage and bone) and cartilage underwent four 5% strain (of cartilage layer) steps with displacement rate of 0.005mm/sec, each followed by a 45-minute relaxation. Final strain was 20%. Bone underwent a single displacement of 20% strain of bone at same displacement rate. Young's moduli ranged from 6.2–42.0MPa, 0.7–3.9MPa, 46.8–123.7MPa for OC graft, cartilage and bone, respectively. The coefficient of variance between OC Grafts, cartilage and bone was 70.6%, 71.8%, and 25.2%, respectively. Dispersion between samples was high. This may be due to intrinsic tissue variability but also due to the testing protocol: for cartilage in particular, the load was at the low end of the load cell capacity and the sample aspect ratio was poor for compressive testing. This work provides insight in understanding the effect of individual patient's and/or individual grafts used during osteochondral grafting. The results compel the need to accurately model these tissues when developing specimen-specific FE models for OC grafting.
Osteochondral (OC) grafting is one available method currently used to repair full thickness cartilage lesions with good results clinically when grafting occurs in patients with specific positive prognostic factors. However, there is poor understanding of the effect of individual patient and surgical factors. With limited tissue availability, development of Finite Element (FE) models taking into account these variations is essential. The aim of this study was to evaluate the effect of altering the material properties of OC grafts and their host environment through computer simulation. A generic FE model (ABAQUS CAE 2017) of a push-out test was developed as a press-fit bone cylinder (graft) sliding inside a bone ring (host tissue). Press-fit fixation was simulated using an interference fit. Overlap between host and graft (0.01mm–0.05mm) and coefficient of friction (0.3–0.7) were varied sequentially. Bone Young's moduli (YM) were varied individually between graft and host within the range of otherwise derived tissue moduli (46MPa, 82MPa, 123MPa). Increasing both overlap and frictional coefficient increased peak dislodging force independently (overlap: 490% & frictional coefficient: 176% across range tested). Increasing bone modulus also increased dislodging force, with host bone modulus (107%, 128%, and 140% increase across range, when Graft YM = 123MPa, 82 MPa, and 46MPa, respectively) having a greater influence than graft modulus (28%, 19% and 10% increase across range, when Host YM = 123 MPa, 82MPa and 46MPa, respectively). As anticipated increasing overlap and friction caused an increase in force necessary to dislodge the graft. Importantly, differentially changing the graft and host material properties changed the dislodging force indicating that difference between graft and host may be an important factor in the success or failure clinically of osteochondral grafting.
Osteoarthritis is a debilitating disease affecting over 1.7 million people in the UK annually. Total ankle replacements are an increasingly sought option for repairing a late stage arthritic ankle, but result in the removal of significant portions of bone regardless of tissue quality. Hence, the mapping of bone quality would allow the use of targeted treatments at earlier stages of the disease. This study aims to develop characterisation methodologies using porcine tissue to investigate the mechanical properties of subchondral bone in the ankle. N=11 talar bone plugs (6mm diameter) were extracted from porcine ankles and embedded into Delrin endcaps using a thin layer of PMMA cement. These were scanned under micro-CT (16 microns) and subjected to quasi-static uniaxial compression to determine apparent stiffness for each specimen. Specimen-specific continuum FE models were developed, with material properties derived from the greyscale value of the underlying image. A python-based least squares regression (Opti4Abq, N=6) was used to minimise the difference between experimental and model stiffness values, to determine the coefficient linking greyscale and mechanical properties. Apparent stiffness, elastic modulus and compressive strength were compared to BV/TV, which was derived using BoneJ (a bone image plugin for the NIH ImageJ). The results show positive correlations between BV/TV and compressive strength, stiffness and Young's modulus. Average BV/TV across all samples was 0.45. Average experimental and computational stiffness were 986N/mm and 891 N/mm respectively. A 21.8% RMS error was found using the validation set (N=5), which was of similar order to the calibration set. Some specimens saw issues with misalignment of the specimen faces and the loading platens, likely causing overestimation of mechanical properties. This study has developed methods that can be translated for use with human ankle bone and will lead to the development of an accurate means of mapping arthritic bone in the ankle.
Intervertebral disc (IVD) degeneration is one of the major causes of back pain. A number of emerging treatments for the condition have failed during clinical trial due to the lack of robust biomechanical testing during product development. The aim of this work was to develop improved in-vitro testing methods to enable new therapeutic approaches to be examined pre-clinically. It forms part of a wider programme of research to develop a minimally invasive nucleus augmentation procedure using self-assembling hydrogels. Previous static testing on extracted IVDs have shown large inter-specimen variation in the measured stiffness when specimen hydration and fluid flow were not well controlled. In this work, a method of normalising the hydration state of IVDs prior-to and during compressive testing was developed. Excised adult bovine IVDs underwent water-pik treatment and a 24-hour agitated bath in monosodium citrate solution to maximise fluid mobility. Specimens were submerged in a saline bath and held under constant pressure for 24 hours, after which the rate of change of displacement was low. Specimens were then cyclically loaded, from which the normalised specimen stiffness was determined. A degenerate disc model was developed with the use of enzymatic degeneration, allowing specimens to be tested sequentially in a healthy, degenerate, and then treated state. Self-assembling peptide-GAG hydrogels were tested using the developed method and the effect of treatment on stiffness and disc height were assessed. Compared to previous static tests, the improved method reduced the variation in the normalised specimen stiffness. In addition, statistically significant differences were seen before and after enzymatic degradation to simulate degeneration, thus providing controls against which to evaluate treatments. The augmentation of the nucleus with the hydrogel intervention reduces the stiffness of the degenerate disc towards that of the healthy disc. This method is now being used to further investigate nucleus augmentation devices.
Acetabular tissue damage is implicated in osteoarthritis (OA) and investigation of in situ acetabular soft tissues behaviour will improve understanding of tissue properties and interconnections. The study aim was to visualise acetabular soft tissues under load and to quantify displacements using computed tomography (CT) scans (XtremeCT, Scano Medical). A CT scan (resolution 82 μm) was performed on the disarticulated, unloaded porcine acetabulum. The femoral head was soaked in Sodium Iodide (NaI) solution and cling film wrapped to prevent transfer to the acetabular side. The joint was realigned, compressed using cable ties and re-scanned. The two images were down-sampled to 0.3 mm. Acetabular bone and soft tissues were segmented. Bony features were used to register the two background images, using Simpleware ScanIP 7.0 (Synopsys), to the same position and orientation (volume difference < 5%). Acetabular soft tissues displacements were measured by tracking the same points at the tissue edges on the two acetabular masks, along with difference in bone position as an additional error assessment. The use of radiopaque solution provided a clear contrast allowing separation of the femoral and acetabular soft tissues in the loaded image. The image registration process resulted in a difference in bone position in the areas of interest equivalent to image resolution (0.3 mm, a mean of 3 repeats by same user). A labral tip displacement of 1.7 mm and a cartilage thickness change from 1.5 mm unloaded to 0.9 mm loaded, were recorded. The combination of contrast enhancement, registration and focused local measurement was precise enough to reduce bone alignment error to that of image resolution and reveal local soft tissue displacements. These measurement methods can be used to develop models of soft tissues properties and behaviour, and therapy for hip tissue damage at early stage may be reviewed and optimised.
Ligaments and tendons are connective tissues with a highly hierarchical structure, from collagen fibres, to fibrils and fascicules. Their intricate structural arrangement produces an anisotropic non-linear elastic mechanical behaviour and a complex damage pattern before failure. Recent constitutive models have been developed with all parameters describing the structure of the tissue, with the advantage that they can in theory be measured on the tissue rather than being phenomenologically-derived. This is an ideal framework to model damage as its onset and propagation can be associated to changes in the structure directly. In this preliminary study, the possibility to identify damage mechanisms in the tissue structure using in silico models was analysed for both the anterior cruciate ligament, with fascicules forming a helix with its longitudinal axis, and the patellar tendon, with fascicules co-aligned with its longitudinal axis. Tissues of interest were modelled as cylinders submitted to uniaxial tension. Damage was modelled as either a reduction of collagen volume fraction with increased strain, assuming the number of collagen fibres sustaining load decreases as fibres fail, or a reduction of the modulus of the fibres, assuming pre-failure damage of the fibres. Each damage mechanism was associated with a damage variable with different fibre stretch threshold for damage initiation and assuming linear variation of damage until an arbitrary failure point. The apparent behaviour of the modelled tissues was significantly different as damage thresholds, damage mechanisms, type of fascicules were varied. This preliminary work showed that using a structural constitutive model to describe occurrence and propagation of structural damage in an in silico model of hierarchical connective tissues is a framework that can clearly differentiate at a macroscopic level between different values of damage threshold and different damage mechanisms for tissue with co-aligned or helical fascicules.
The clinical uptake of minimally invasive interventions for intervertebral disc, such as nucleus augmentation, is currently hampered by the lack of robust pre-clinical testing methods that can take into account the large variation in the mechanical behaviour of the tissues. Using computational modelling to develop new interventions could be a way to test scenarios accounting for variation. However, such models need to have been validated for relevant mechanical function, e.g. compressive, torsional or flexional stiffness, and local disc deformations. The aim of this work was to use a novel in-vitro imaging method to assess the validity of computational models of the disc that employed different degrees of sophistication in the anatomical representation of the nucleus. Bovine caudal bone-disc-bone entities (N=6) were dissected and tested in uniaxial compression in a custom-made rig. Forty glass markers were placed on the surface of each disc. The specimens were scanned both with MRI and micro-CT before and during loading. Specimen-specific computational models were built from CT images to replicate the compression test. The anatomy of the nucleus was represented in three ways: assuming a standard diameter ratio, assuming a cylindrical shape with its volume matching that measured from MRI, and deriving the shape directly from MRI. The three types of models were calibrated for force-displacement. The radial displacement of the glass markers were then compared with their experimental displacement derived from CT images. For a similar accuracy in modelling overall force-displacement, the mean error on the surface displacement was 35% for standard ratio nucleus, 38% for image-based cylindrical nucleus, and 32% for MRI-based nucleus geometry. This work shows that, as long as consistency is kept to develop and calibrate image-based computational models, the complexity of the nucleus geometry does not influence the ability of a model to predict surface displacement in the intervertebral disc.
Geometric variations of the hip joint can give rise to repeated abnormal contact between the femur and acetabular rim, resulting in cartilage and labrum damage. Population-based geometric parameterisation can facilitate the flexible and automated Parametric geometries can be developed into finite element models to assess the potential effects of morphological variations in bone on soft tissue strains. The aim of this study was to demonstrate the capabilities of our parameterisation research tool by assessing impingement severity resulting from a range of parametrically varied femoral and acetabular geometries. Custom made MATLAB (MathWorks) and Python codes[1] were used to generate bone surfaces, which were developed into finite element models in Abaqus (SIMULIA). Parametric femoral surfaces were defined by a spherical proximal head and ellipse sections through the neck/cam region. This method produced surfaces that were well fitted to bone geometry segmented from CT scans of cam patients and capable of producing trends in results similar to those found using segmented models. A simplified spherical geometry, including the labrum and acetabular cartilage, represented the acetabulum. Femoral parameters were adjusted to define relevant variations in cam size and position. Two radii (small and large cams) and two positions (anterior and superior cams) were defined resulting in four models. Alpha angles of these parametric femurs were measured in an anterior-posterior view and a cross-table lateral view using ImageJ (NIH). A further model was developed using a femur with a medium cam size and position, and the level of acetabular coverage and labrum length were varied. Bones were modelled as rigid bodies and soft tissues were modelled as transversely isotropic linearly elastic materials. With the acetabulum fully constrained in all cases, the femurs were constrained in translation and rotated to simulate flexion followed by internal rotation to cause impingement against the labrum.Introduction
Methods