Head-Neck Osteoplasty has Minor Effect on the Strength of an Ovine Cam-FAI Model: In Vitro and Finite Element Analyses

Background Osteochondroplasty of the head-neck region is performed on patients with cam femoroacetabular impingement (FAI) without fully understanding its repercussion on the integrity of the femur. Cam-type FAI can be surgically and reproducibly induced in the ovine femur, which makes it suitable for studying corrective surgery in a consistent way. Finite element models built on quantitative CT (QCT) are computer tools that can be used to predict femoral strength and evaluate the mechanical effect of surgical correction. Questions/purposes We asked: (1) What is the effect of a resection of the superolateral aspect of the ovine femoral head-neck junction on failure load? (2) How does the failure load after osteochondroplasty compare with reported forces from activities of daily living in sheep? (3) How do failure loads and failure locations from the computer simulations compare with the experiments? Methods Osteochondroplasties (3, 6, 9 mm) were performed on one side of 18 ovine femoral pairs with the contralateral intact side as a control. The 36 femurs were scanned via QCT from which specimen-specific computer models were built. Destructive compression tests then were conducted experimentally using a servohydraulic testing system and numerically via the computer models. Safety factors were calculated as the ratio of the maximal force measured in vivo by telemeterized hip implants during the sheep’s walking and running activities to the failure load. The simulated failure loads and failure locations from the computer models were compared with the experimental results. Results Failure loads were reduced by 5% (95% CI, 2%–8%) for the 3-mm group (p = 0.0089), 10% (95% CI, 6%–14%) for the 6-mm group (p = 0.0015), and 19% (95% CI, 13%–26%) for the 9-mm group (p = 0.0097) compared with the controls. Yet, the weakest specimen still supported more than 2.4 times the peak load during running. Strong correspondence was found between the simulated and experimental failure loads (R2 = 0.83; p < 0.001) and failure locations. Conclusions The resistance of ovine femurs to fracture decreased with deeper resections. However, under in vitro testing conditions, the effect on femoral strength remains small even after 9 mm correction, suggesting that femoral head-neck osteochondroplasty could be done safely on the ovine femur. QCT-based finite element models were able to predict weakening of the femur resulting from the osteochondroplasty. Clinical Relevance The ovine femur provides a seemingly safe platform for scientific evaluation of FAI. It also appears that computer models based on preoperative CT scans may have the potential to provide patient-specific guidelines for preventing overcorrection of cam FAI

Image-based biomechanical assessment of vertebral body and intervertebral disc in the human lumbar spine

Life expectancy continuously increases but our society faces age-related conditions. Among musculoskeletal diseases, osteoporosis associated with risk of vertebral fracture and degenerative intervertebral disc (IVD) are painful pathologies responsible for tremendous healthcare costs. Hence, reliable diagnostic tools are necessary to plan a treatment or follow up its efficacy. Yet, radiographic and MRI techniques, respectively clinical standards for evaluation of bone strength and IVD degeneration, are unspecific and not objective. Increasingly used in biomedical engineering, CT-based finite element (FE) models constitute the state-of-art for vertebral strength prediction. However, as non-invasive biomechanical evaluation and personalised FE models of the IVD are not available, rigid boundary conditions (BCs) are applied on the FE models to avoid uncertainties of disc degeneration that might bias the predictions. Moreover, considering the impact of low back pain, the biomechanical status of the IVD is needed as a criterion for early disc degeneration. Thus, the first FE study focuses on two rigid BCs applied on the vertebral bodies during compression test of cadaver vertebral bodies, vertebral sections and PMMA embedding. The second FE study highlights the large influence of the intervertebral disc’s compliance on the vertebral strength, damage distribution and its initiation. The third study introduces a new protocol for normalisation of the IVD stiffness in compression, torsion and bending using MRI-based data to account for its morphology. In the last study, a new criterion (Otsu threshold) for disc degeneration based on quantitative MRI data (axial T2 map) is proposed. The results show that vertebral strength and damage distribution computed with rigid BCs are identical. Yet, large discrepancies in strength and damage localisation were observed when the vertebral bodies were loaded via IVDs. The normalisation protocol attenuated the effect of geometry on the IVD stiffnesses without complete suppression. Finally, the Otsu threshold computed in the posterior part of annulus fibrosus was related to the disc biomechanics and meet objectivity and simplicity required for a clinical application. In conclusion, the stiffness normalisation protocol necessary for consistent IVD comparisons and the relation found between degeneration, mechanical response of the IVD and Otsu threshold lead the way for non-invasive evaluation biomechanical status of the IVD. As the FE prediction of vertebral strength is largely influenced by the IVD conditions, this data could also improve the future FE models of osteoporotic vertebra

The effective elastic properties of human trabecular bone may be approximated using micro-finite element analyses of embedded volume elements

Boundary conditions (BCs) and sample size affect the measured elastic properties of cancellous bone. Samples too small to be representative appear stiffer under kinematic uniform BCs (KUBCs) than under periodicity-compatible mixed uniform BCs (PMUBCs). To avoid those effects, we propose to determine the effective properties of trabecular bone using an embedded configuration. Cubic samples of various sizes (2.63, 5.29, 7.96, 10.58 and 15.87 mm) were cropped from μCT scans of femoral heads and vertebral bodies. They were converted into μFE models and their stiffness tensor was established via six uniaxial and shear load cases. PMUBCs- and KUBCs-based tensors were determined for each sample. “In situ” stiffness tensors were also evaluated for the embedded configuration, i.e. when the loads were transmitted to the samples via a layer of trabecular bone. The Zysset–Curnier model accounting for bone volume fraction and fabric anisotropy was fitted to those stiffness tensors, and model parameters ν0 Poisson’s ratio) E 0 and μ0 (elastic and shear moduli) were compared between sizes. BCs and sample size had little impact on ν0. However, KUBCs- and PMUBCs-based E 0 and μ0 , respectively, decreased and increased with growing size, though convergence was not reached even for our largest samples. Both BCs produced upper and lower bounds for the in situ values that were almost constant across samples dimensions, thus appearing as an approximation of the effective properties. PMUBCs seem also appropriate for mimicking the trabecular core, but they still underestimate its elastic properties (especially in shear) even for nearly orthotropic samples

Embedding of human vertebral bodies leads to higher ultimate load and altered damage localisation under axial compression

Computer tomography (CT)-based finite element (FE) models of vertebral bodies assess fracture load in vitro better than dual energy X-ray absorptiometry, but boundary conditions affect stress distribution under the endplates that may influence ultimate load and damage localisation under post-yield strains. Therefore, HRpQCT-based homogenised FE models of 12 vertebral bodies were subjected to axial compression with two distinct boundary conditions: embedding in polymethylmethalcrylate (PMMA) and bonding to a healthy intervertebral disc (IVD) with distinct hyperelastic properties for nucleus and annulus. Bone volume fraction and fabric assessed from HRpQCT data were used to determine the elastic, plastic and damage behaviour of bone. Ultimate forces obtained with PMMA were 22% higher than with IVD but correlated highly (R2 = 0.99). At ultimate force, distinct fractions of damage were computed in the endplates (PMMA: 6%, IVD: 70%), cortex and trabecular sub-regions, which confirms previous observations that in contrast to PMMA embedding, failure initiated underneath the nuclei in healthy IVDs. In conclusion, axial loading of vertebral bodies via PMMA embedding versus healthy IVD overestimates ultimate load and leads to distinct damage localisation and failure pattern

μCT-based trabecular anisotropy can be reproducibly computed from HR-pQCT scans using the triangulated bone surface

The trabecular structure can be assessed at the wrist or tibia via high-resolution peripheral quantitative computed tomography (HR-pQCT). Yet on this modality, the performance of the existing methods, evaluating trabecular anisotropy is usually overlooked, especially in terms of reproducibility. We thus proposed to compare the TRI routine used by SCANCO Medical AG (Brüttisellen, Switzerland), the classical mean intercept length (MIL), and the grey-level structure tensor (GST) to the mean surface length (MSL), a new method for evaluating a second-order fabric tensor based on the triangulation of the bone surface. The distal radius of 24 fresh-frozen human forearms was scanned three times via HR-pQCT protocols (61 μm, 82 μm nominal voxel size), dissected, and imaged via micro computed tomography (μCT) at 16 μm nominal voxel size. After registering the scans, we compared for each resolution the fabric tensors, determined by the mentioned techniques for 182 trabecular regions of interest. We then evaluated the reproducibility of the fabric information measured by HR-pQCT via precision errors. On μCT, TRI and GST were respectively the best and worst surrogates for MILμCT (MIL computed on μCT) in terms of eigenvalues and main direction of anisotropy. On HR-pQCT, however, MSL provided the best approximation of MILμCT. Surprisingly, surface-based approaches (TRI, MSL) also proved to be more precise than both MIL and GST. Our findings confirm that MSL can reproducibly estimate MILμCT, the current gold standard. MSL thus enables the direct mapping of the fabric-dependent material properties required in homogenised HR-pQCT-based finite element models