Compressive strength of interbody cages in the lumbar spine: the effect of cage shape, posterior instrumentation and bone density

One goal of interbody fusion is to increase the height of the degenerated disc space. Interbody cages in particular have been promoted with the claim that they can maintain the disc space better than other methods. There are many factors that can affect the disc height maintenance, including graft or cage design, the quality of the surrounding bone and the presence of supplementary posterior fixation. The present study is an in vitro biomechanical investigation of the compressive behaviour of three different interbody cage designs in a human cadaveric model. The effect of bone density and posterior instrumentation were assessed. Thirty-six lumbar functional spinal units were instrumented with one of three interbody cages: (1) a porous titanium implant with endplate fit (Stratec), (2) a porous, rectangular carbon-fibre implant (Brantigan) and (3) a porous, cylindrical threaded implant (Ray). Posterior instrumentation (USS) was applied to half of the specimens. All specimens were subjected to axial compression displacement until failure. Correlations between both the failure load and the load at 3 mm displacement with the bone density measurements were observed. Neither the cage design nor the presence of posterior instrumentation had a significant effect on the failure load. The loads at 3 mm were slightly less for the Stratec cage, implying lower axial stiffness, but were not different with posterior instrumentation. The large range of observed failure loads overlaps the potential in vivo compressive loads, implying that failure of the bone-implant interface may occur clinically. Preoperative measurements of bone density may be an effective tool to predict settling around interbody cages.

Compressive strength of elderly vertebrae is reduced by disc degeneration and additional flexion

Computer tomography (CT)-based finite element (FE) models assess vertebral strength better than dual energy X-ray absorptiometry. Osteoporotic vertebrae are usually loaded via degenerated intervertebral discs (IVD) and potentially at higher risk under forward bending, but the influences of the IVD and loading conditions are generally overlooked. Accordingly, magnetic resonance imaging was performed on 14 lumbar discs to generate FE models for the healthiest and most degenerated specimens. Compression, torsion, bending, flexion and extension conducted experimentally were used to calibrate both models. They were combined with CT-based FE models of 12 lumbar vertebral bodies to evaluate the effect of disc degeneration compared to a loading via endplates embedded in a stiff resin, the usual experimental paradigm. Compression and lifting were simulated, load and damage pattern were evaluated at failure. Adding flexion to the compression (lifting) and higher disc degeneration reduces the failure load (8–14%, 5–7%) and increases damage in the vertebrae. Under both loading scenarios, decreasing the disc height slightly increases the failure load; embedding and degenerated IVD provides respectively the highest and lowest failure load. Embedded vertebrae are more brittle, but failure loads induced via IVDs correlate highly with vertebral strength. In conclusion, osteoporotic vertebrae with degenerated IVDs are consistently weaker—especially under lifting, but clinical assessment of their strength is possible via FE analysis without extensive disc modelling, by extrapolating measures from the embedded situation.

Clinical versus pre-clinical FE models for vertebral body strength predictions

The finite element analysis is an accepted method to predict vertebral body compressive strength. This study compares measurements obtained from in vitro tests with the ones from two different simulation models: clinical quantitative computer tomography (QCT) based homogenized finite element (hFE) models and pre-clinical high-resolution peripheral QCT-based (HR-pQCT) hFE models. About 37 vertebral body sections were prepared by removing end-plates and posterior elements, scanned with QCT (390/450μm voxel size) as well as HR-pQCT (82μm voxel size), and tested in compression up to failure. Non-linear viscous damage hFE models were created from QCT/HT-pQCT images and compared to experimental results based on stiffness and ultimate load. As expected, the predictability of QCT/HR-pQCT-based hFE models for both apparent stiffness (r2=0.685/0.801r2=0.685/0.801) and strength (r2=0.774/0.924r2=0.774/0.924) increased if a better image resolution was used. An analysis of the damage distribution showed similar damage locations for all cases. In conclusion, HR-pQCT-based hFE models increased the predictability considerably and do not need any tuning of input parameters. In contrast, QCT-based hFE models usually need some tuning but are clinically the only possible choice at the moment.

Gravitational slopes, geomorphology, and material strengths of the nucleus of comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Aims. We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material (tensile, shear, and compressive strengths). Methods. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface (Imhotep, Ash, Seth, Hathor, and Agilkia), using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies (overhangs, collapsed structures, boulders, cliffs, and Philae's footprint) and mechanical considerations. Results. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 degrees) are covered by a fine material and contain a few large (>10 m) and isolated boulders; ii) intermediate-slope terrains (20-45 degrees) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from <1m to 10m for the majority of them; and iii) high-slope terrains (45-90 degrees) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

Thermal and mechanical properties of the near-surface layers of comet 67P/Churyumov-Gerasimenko

Thermal and mechanical material properties determine comet evolution and even solar system formation because comets are considered remnant volatile-rich planetesimals. Using data from the Multipurpose Sensors for Surface and Sub-Surface Science (MUPUS) instrument package gathered at the Philae landing site Abydos on comet 67P/Churyumov-Gerasimenko, we found the diurnal temperature to vary between 90 and 130 K. The surface emissivity was 0.97, and the local thermal inertia was 85 +/- 35 J m(-2) K(-1)s(-1/2). The MUPUS thermal probe did not fully penetrate the near-surface layers, suggesting a local resistance of the ground to penetration of >4 megapascals, equivalent to >2 megapascal uniaxial compressive strength. A sintered near-surface microporous dust-ice layer with a porosity of 30 to 65% is consistent with the data.

Bone mineral density (BMD) and computer tomographic measurements of the equine proximal phalanx in correlation with breaking strength

Despite the fact that bone mineral density (BMD) is an important fracture risk predictor in human medicine, studies in equine orthopedic research are still lacking. We hypothesized that BMD correlates with bone failure and fatigue fractures of this bone. Thus, the objectives of this study were to measure the structural and mechanical properties of the proximal phalanx with dual energy X-ray absorptiometry (DXA), to correlate the data obtained from DXA and computer tomography (CT) measurements to those obtained by loading pressure examination and to establish representative region of interest (ROI) for in vitro BMD measurements of the equine proximal phalanx for predicting bone failure force. DXA was used to measure the whole bone BMD and additional three ROI sites in 14 equine proximal phalanges. Following evaluation of the bone density, whole bone, cortical width and area in the mid-diaphyseal plane were measured on CT images. Bones were broken using a manually controlled universal bone crusher to measure bone failure force and reevaluated for the site of fractures on follow-up CT images. Compressive load was applied at a constant displacement rate of 2 mm/min until failure, defined as the first clear drop in the load measurement. The lowest BMD was measured at the trabecular region (mean +/- SD: 1.52 +/- 0.12 g/cm2; median: 1.48 g/cm2; range: 1.38-1.83 g/cm2). There was a significant positive linear correlation between trabelcular BMD and the breaking strength (P = 0.023, r = 0.62). The trabecular region of the proximal phalanx appears to be the only significant indicator of failure of strength in vitro. This finding should be reassessed to further reveal the prognostic value of trabecular BMD in an in vivo fracture risk model.

Experimental Validation of Finite Element Analysis of Human Vertebral Collapse Under Large Compressive Strains

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towardsthe prediction of extended vertebral collapse may help in assessing fracture stability in future work.

The Initial Slope of the Variogram, Foundation of the Trabecular Bone Score, Is Not or Poorly Associated With Vertebral Strength

Trabecular bone score (TBS) rests on the textural analysis of DXA to reflect the decay in trabecular structure characterising osteoporosis. Yet, its discriminative power in fracture studies remains incomprehensible as prior biomechanical tests found no correlation with vertebral strength. To verify this result possibly due to an unrealistic set-up and to cover a wide range of loading scenarios, the data from three previous biomechanical studies using different experimental settings was used. They involved the compressive failure of 62 human lumbar vertebrae loaded 1) via intervertebral discs to mimic the in vivo situation (“full vertebra”), 2) via the classical endplate embedding (“vertebral body”) or 3) via a ball joint to induce anterior wedge failure (“vertebral section”). HR-pQCT scans acquired prior testing were used to simulate anterior-posterior DXA from which areal bone mineral density (aBMD) and the initial slope of the variogram (ISV), the early definition of TBS, were evaluated. Finally, the relation of aBMD and ISV with failure load (Fexp) and apparent failure stress (σexp) was assessed and their relative contribution to a multi-linear model was quantified via ANOVA. We found that, unlike aBMD, ISV did not significantly correlate with Fexp and σexp, except for the “vertebral body” case (r2 = 0.396, p = 0.028). Aside from the “vertebra section” set-up where it explained only 6.4% of σexp (p = 0.037), it brought no significant improvement to aBMD. These results indicate that ISV, a replica of TBS, is a poor surrogate for vertebral strength no matter the testing set-up, which supports the prior observations and raises a fortiori the question of the deterministic factors underlying the statistical relationship between TBS and vertebral fracture risk.