Presence of intervertebral discs alters observed stiffness and failure mechanisms in the vertebra

Subsidence of a partially porous titanium lumbar cage produced by electron beam melting technology

The lumbar intervertebral devices are widely used in the surgical treatment of lumbar diseases. The subsidence represents a serious clinical issue during the healing process, mainly when the interfaces between the implant and the vertebral bodies are not well designed. The aim of this study is the evaluation of subsidence risk for two different devices. The devices have the same shape, but one of them includes a filling micro lattice structure. The effect of the micro lattice structure on the subsidence behavior of the implant was evaluated by means of both experimental tests and finite element analyses. Compressive tests were carried out by using blocks made of grade 15 polyurethane, which simulate the vertebral bone. Non-linear, quasi-static finite element analyses were performed to simulate experimental and physiologic conditions. The experimental tests and the FE analyses showed that the subsidence risk is higher for the device without micro lattice structure, due to the smaller contact surface. Moreover, an overload in the central zone of the contact surface was detected in the same device and it could cause the implant failure. Thus, the micro lattice structure allows a homogenous pressure distribution at the implant-bone interface.

Measurement of bone damage caused by quasi-static compressive loading-unloading to explore dental implants stability: Simultaneous use of in-vitro tests, μ-CT images, and digital volume correlation

Primary stability of dental implants is the initial mechanical engagement of the implant with its adjacent bone. Implantation and the subsequent loading may cause mechanical damage in the peripheral bone, which ultimately reduces the stability of the implant. This study aimed at evaluating primary stability of dental implants through applying stepwise compressive displacement-controlled, loading-unloading cycles to obtain overall stiffness and dissipated energy of the bone-implant structure; and quantifying induced plastic strains in surrounding bone using digital volume correlation (DVC) method, through comparing μCT images in different loading steps. To this end, dental implants were inserted into the cylindrical trabecular bones, then the bone-implant structure was undergone step-wise loading-unloading cycles, and μCT images were taken in some particular steps, then comparison was made between undeformed and deformed configurations using DVC to quantify plastic strain within the trabecular bone. Comparing stiffness reduction and dissipated energy values in different loading steps, obtained from the force-displacement curve in each loading step, revealed that the maximum displacement of 0.16 mm can be deemed as a safe threshold above which damages in peri-implant bone started to increase considerably (p < 0.05). In addition, it was found here that peri-implant bone strain linearly increased with decreasing bone-implant stiffness (p < 0.05). Moreover, strain concentration in peri-implant bone region showed that the plastic strain in trabecular bone spread up to a distance of about 2.5 mm away from the implant surface. Research of this kind can be used to optimize the design of dental implants, with the ultimate goal of improving their stability, also to validate in-silico models, e.g., micro-finite element models, which can help gain a deeper understanding of bone-implant construct behavior.

MicroFE models of porcine vertebrae with induced bone focal lesions: Validation of predicted displacements with digital volume correlation

The evaluation of the local mechanical behavior as a result of metastatic lesions is fundamental for the characterization of the mechanical competence of metastatic vertebrae. Micro finite element (microFE) models have the potential of addressing this challenge through laboratory studies but their predictions of local deformation due to the complexity of the bone structure compromized by the lesion must be validated against experiments. In this study, the displacements predicted by homogeneous, linear and isotropic microFE models of vertebrae were validated against experimental Digital Volume Correlation (DVC) measurements. Porcine spine segments, with and without mechanically induced focal lesions, were tested in compression within a micro computed tomography (microCT) scanner. The displacement within the bone were measured with an optimized global DVC approach (BoneDVC). MicroFE models of the intact and lesioned vertebrae, including or excluding the growth plates, were developed from the microCT images. The microFE and DVC boundary conditions were matched. The displacements measured by the DVC and predicted by the microFE along each Cartesian direction were compared. The results showed an excellent agreement between the measured and predicted displacements, both for intact and metastatic vertebrae, in the middle of the vertebra, in those cases where the structure was not loaded beyond yield (0.69 < R² < 1.00). Models with growth plates showed the worst correlations (0.02 < R² < 0.99), while a clear improvement was observed if the growth plates were excluded (0.56 < R² < 1.00). In conclusion, these simplified models can predict complex displacement fields in the elastic regime with high reliability, more complex non-linear models should be implemented to predict regions with high deformation, when the bone is loaded beyond yield.

The opening size of the laminoplasty is dependent on the groove size: A numerical study

BACKGROUND

The expansion of the cervical vertebrae lamina appears to be crucial to related surgical procedures. The dimensions of the groove influence the strain concentration within the lamina of the vertebra and, thus, the potential success or failure of respective surgical procedure. The aim of this computational study is to clarify both the role of the size of the groove with concern to both the open door and the double door laminoplasty techniques.

METHODS

Finite element models were created via computer tomography with varying lamina groove dimensions. Displacements were applied to the models at the open side of the vertebral arch and the vertebral body was constrained prior to movement along all the axes. The maximal opening size measured on the inner side of the lamina and the percentage increase in the initial spinal areas were subsequently analyzed.

FINDINGS

The elastic strain concentration value was observed for the groove in all cases, while the maximal principal elastic strain concentration value was observed at the opposite side to the groove cut into the lamina, also in all cases. The maximal area increase related to the 4 mm groove accompanied by the preservation of the ventral cortex of the bone.

INTERPRETATION

The study suggested three conclusions a) the wider the groove, the greater is the opening potential, b) the maximal opening size following laminoplasty is not dependent on the depth of the bone cut for this type of groove, c) no benefit accrues in terms of the opening size following the cutting of a supplementary groove at the beginning of the lamina.

Assessment of Intravertebral Mechanical Strains and Cancellous Bone Texture Under Load Using a Clinically Available Digital Tomosynthesis Modality

Vertebral fractures are the most common osteoporotic fractures, but clinical means for assessment of vertebral bone integrity are limited in accuracy, as they typically use surrogate measures that are indirectly related to mechanics. The objective of this study was to examine the extent to which intravertebral strain distributions and changes in cancellous bone texture generated by a load of physiological magnitude can be characterized using a clinically available imaging modality. We hypothesized that digital tomosynthesis-based digital volume correlation (DTS-DVC) and image texture-based metrics of cancellous bone microstructure can detect development of mechanical strains under load. Isolated cadaveric T11 vertebrae and L2-L4 vertebral segments were DTS imaged in a nonloaded state and under physiological load levels. Axial strain, maximum principal strain, maximum compressive and tensile principal strains, and von Mises equivalent strain were calculated using the DVC technique. The change in textural parameters (line fraction deviation, anisotropy, and fractal parameters) under load was calculated within the cancellous centrum. The effect of load on measured strains and texture variables was tested using mixed model analysis of variance, and relationships of strain and texture variables with donor age, bone density parameters, and bone size were examined using regression models. Magnitudes and heterogeneity of intravertebral strain measures correlated with applied loading and were significantly different from background noise. Image texture parameters were found to change with applied loading, but these changes were not observed in the second experiment testing L2-L4 segments. DTS-DVC-derived strains correlated with age more strongly than did bone mineral density (BMD) for T11.

A novel approach to evaluate the effects of artificial bone focal lesion on the three-dimensional strain distributions within the vertebral body

The spine is the first site for incidence of bone metastasis. Thus, the vertebrae have a high potential risk of being weakened by metastatic tissues. The evaluation of strength of the bone affected by the presence of metastases is fundamental to assess the fracture risk. This work proposes a robust method to evaluate the variations of strain distributions due to artificial lesions within the vertebral body, based on in situ mechanical testing and digital volume correlation. Five porcine vertebrae were tested in compression up to 6500N inside a micro computed tomography scanner. For each specimen, images were acquired before and after the application of the load, before and after the introduction of the artificial lesions. Principal strains were computed within the bone by means of digital volume correlation (DVC). All intact specimens showed a consistent strain distribution, with peak minimum principal strain in the range -1.8% to -0.7% in the middle of the vertebra, demonstrating the robustness of the method. Similar distributions of strains were found for the intact vertebrae in the different regions. The artificial lesion generally doubled the strain in the middle portion of the specimen, probably due to stress concentrations close to the defect. In conclusion, a robust method to evaluate the redistribution of the strain due to artificial lesions within the vertebral body was developed and will be used in the future to improve current clinical assessment of fracture risk in metastatic spines.

Effect of fabric on the accuracy of computed tomography-based finite element analyses of the vertebra

Quantitative computed tomography (QCT)-based finite element (FE) models of the vertebra are widely used in studying spine biomechanics and mechanobiology, but their accuracy has not been fully established. Although the models typically assign material properties based only on local bone mineral density (BMD), the mechanical behavior of trabecular bone also depends on fabric. The goal of this study was to determine the effect of incorporating measurements of fabric on the accuracy of FE predictions of vertebral deformation. Accuracy was assessed by using displacement fields measured via digital volume correlation-applied to time-lapse microcomputed tomography (μCT)-as the gold standard. Two QCT-based FE models were generated from human L1 vertebrae (n = 11): the entire vertebral body and a cuboid-shaped portion of the trabecular centrum [dimensions: (20-30) × (15-20) × (15-20) mm3]. For axial compression boundary conditions, there was no difference (p = 0.40) in the accuracy of the FE-computed displacements for models using material properties based on local values of BMD versus those using material properties based on local values of fabric and volume fraction. However, when using BMD-based material properties, errors were higher for the vertebral-body models (8.4-50.1%) than cuboid models (1.5-19.6%), suggesting that these properties are inaccurate in the peripheral regions of the centrum. Errors also increased when assuming that the cuboid region experienced uniaxial loading during axial compression of the vertebra. These findings indicate that a BMD-based constitutive model is not sufficient for the peripheral region of the vertebral body when seeking accurate QCT-based FE modeling of the vertebra.

Effect of SR-microCT radiation on the mechanical integrity of trabecular bone using in situ mechanical testing and digital volume correlation

The use of synchrotron radiation micro-computed tomography (SR-microCT) is becoming increasingly popular for studying the relationship between microstructure and bone mechanics subjected to in situ mechanical testing. However, it is well known that the effect of SR X-ray radiation can considerably alter the mechanical properties of bone tissue. Digital volume correlation (DVC) has been extensively used to compute full-field strain distributions in bone specimens subjected to step-wise mechanical loading, but tissue damage from sequential SR-microCT scans has not been previously addressed. Therefore, the aim of this study is to examine the influence of SR irradiation-induced microdamage on the apparent elastic properties of trabecular bone using DVC applied to in situ SR-microCT tomograms obtained with different exposure times. Results showed how DVC was able to identify high local strain levels (> 10,000 µε) corresponding to visible microcracks at high irradiation doses (~ 230 kGy), despite the apparent elastic properties remained unaltered. Microcracks were not detected and bone plasticity was preserved for low irradiation doses (~ 33 kGy), although image quality and consequently, DVC performance were reduced. DVC results suggested some local deterioration of tissue that might have resulted from mechanical strain concentration further enhanced by some level of local irradiation even for low accumulated dose.

Femoral fracture type can be predicted from femoral structure: A finite element study validated by digital volume correlation experiments

Proximal femoral fractures can be categorized into two main types: Neck and intertrochanteric fractures accounting for 53% and 43% of all proximal femoral fractures, respectively. The possibility to predict the type of fracture a specific patient is predisposed to would allow drug and exercise therapies, hip protector design, and prophylactic surgery to be better targeted for this patient rendering fracture preventing strategies more effective. This study hypothesized that the type of fracture is closely related to the patient-specific femoral structure and predictable by finite element (FE) methods. Fourteen femora were DXA scanned, CT scanned, and mechanically tested to fracture. FE-predicted fracture patterns were compared to experimentally observed fracture patterns. Measurements of strain patterns to explain neck and intertrochanteric fracture patterns were performed using a digital volume correlation (DVC) technique and compared to FE-predicted strains and experimentally observed fracture patterns. Although loaded identically, the femora exhibited different fracture types (six neck and eight intertrochanteric fractures). CT-based FE models matched the experimental observations well (86%) demonstrating that the fracture type can be predicted. DVC-measured and FE-predicted strains showed obvious consistency. Neither DXA-based BMD nor any morphologic characteristics such as neck diameter, femoral neck length, or neck shaft angle were associated with fracture type. In conclusion, patient-specific femoral structure correlates with fracture type and FE analyses were able to predict these fracture types. Also, the demonstration of FE and DVC as metrics of the strains in bones may be of substantial clinical value, informing treatment strategies and device selection and design. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:993-1001, 2018.

Effect of the In Vitro Boundary Conditions on the Surface Strain Experienced by the Vertebral Body in the Elastic Regime

The vertebral strength and strain can be assessed in vitro by both using isolated vertebrae and sets of three adjacent vertebrae (the central one is loaded through the disks). Our goal was to elucidate if testing single-vertebra-specimens in the elastic regime provides different surface strains to three-vertebrae-segments. Twelve three-vertebrae sets were extracted from thoracolumbar human spines. To measure the principal strains, the central vertebra of each segment was prepared with eight strain-gauges. The sets were tested mechanically, allowing comparison of the surface strains between the two boundary conditions: first when the same vertebra was loaded through the disks (three-vertebrae-segment) and then with the endplates embedded in cement (single-vertebra). They were all subjected to four nondestructive tests (compression, traction, torsion clockwise, and counterclockwise). The magnitude of principal strains differed significantly between the two boundary conditions. For axial loading, the largest principal strains (along vertebral axis) were significantly higher when the same vertebra was tested isolated compared to the three-vertebrae-segment. Conversely, circumferential strains decreased significantly in the single vertebrae compared to the three-vertebrae-segment, with some variations exceeding 100% of the strain magnitude, including changes from tension to compression. For torsion, the differences between boundary conditions were smaller. This study shows that, in the elastic regime, when the vertebra is loaded through a cement pot, the surface strains differ from when it is loaded through the disks. Therefore, when single vertebrae are tested, surface strain should be taken with caution.

Application of digital volume correlation to study the efficacy of prophylactic vertebral augmentation

BACKGROUND

Prophylactic augmentation is meant to reinforce the vertebral body, but in some cases it is suspected to actually weaken it. Past studies only investigated structural failure and the surface strain distribution. To elucidate the failure mechanism of the augmented vertebra, more information is needed about the internal strain distribution. This study aims to measure, for the first time, the full-field three-dimensional strain distribution inside augmented vertebrae in the elastic regime and to failure.

METHODS

Eight porcine vertebrae were prophylactically-augmented using two augmentation materials. They were scanned with a micro-computed tomography scanner (38.8μm voxel resolution) while undeformed, and loaded at 5%, 10%, and 15% compressions. Internal strains (axial, antero-posterior and lateral-lateral components) were computed using digital volume correlation.

FINDINGS

For both augmentation materials, the highest strains were measured in the regions adjacent to the injected cement mass, whereas the cement-interdigitated-bone was less strained. While this was already visible in the elastic regime (5%), it was a predictor of the localization of failure, which became visible at higher degrees of compression (10% and 15%), when failure propagated across the trabecular bone. Localization of high strains and failure was consistent between specimens, but different between the cement types.

INTERPRETATION

This study indicated the potential of digital volume correlation in measuring the internal strain (elastic regime) and failure in augmented vertebrae. While the cement-interdigitated region becomes stiffer (less strained), the adjacent non-augmented trabecular bone is affected by the stress concentration induced by the cement mass. This approach can help establish better criteria to improve vertebroplasty.

Bone mechanical properties and changes with osteoporosis

This review will define the role of collagen and within-bone heterogeneity and elaborate the importance of trabecular and cortical architecture with regard to their effect on the mechanical strength of bone. For each of these factors, the changes seen with osteoporosis and ageing will be described and how they can compromise strength and eventually lead to bone fragility.

Three-Dimensional Local Measurements of Bone Strain and Displacement: Comparison of Three Digital Volume Correlation Approaches

Different digital volume correlation (DVC) approaches are currently available or under development for bone tissue micromechanics. The aim of this study was to compare accuracy and precision errors of three DVC approaches for a particular three-dimensional (3D) zero-strain condition. Trabecular and cortical bone specimens were repeatedly scanned with a micro-computed tomography (CT). The errors affecting computed displacements and strains were extracted for a known virtual translation, as well as for repeated scans. Three DVC strategies were tested: two local approaches, based on fast-Fourier-transform (DaVis-FFT) or direct-correlation (DaVis-DC), and a global approach based on elastic registration and a finite element (FE) solver (ShIRT-FE). Different computation subvolume sizes were tested. Much larger errors were found for the repeated scans than for the virtual translation test. For each algorithm, errors decreased asymptotically for larger subvolume sizes in the range explored. Considering this particular set of images, ShIRT-FE showed an overall better accuracy and precision (a few hundreds microstrain for a subvolume of 50 voxels). When the largest subvolume (50–52 voxels) was applied to cortical bone, the accuracy error obtained for repeated scans with ShIRT-FE was approximately half of that for the best local approach (DaVis-DC). The difference was lower (250 microstrain) in the case of trabecular bone. In terms of precision, the errors shown by DaVis-DC were closer to the ones computed by ShIRT-FE (differences of 131 microstrain and 157 microstrain for cortical and trabecular bone, respectively). The multipass computation available for DaVis software improved the accuracy and precision only for the DaVis-FFT in the virtual translation, particularly for trabecular bone. The better accuracy and precision of ShIRT-FE, followed by DaVis-DC, were obtained with a higher computational cost when compared to DaVis-FFT. The results underline the importance of performing a quantitative comparison of DVC methods on the same set of samples by using also repeated scans, other than virtual translation tests only. ShIRT-FE provides the most accurate and precise results for this set of images. However, both DaVis approaches show reasonable results for large nodal spacing, particularly for trabecular bone. Finally, this study highlights the importance of using sufficiently large subvolumes, in order to achieve better accuracy and precision.

Augmentation of failed human vertebrae with critical un-contained lytic defect restores their structural competence under functional loading: An experimental study

BACKGROUND

Lytic spinal lesions reduce vertebral strength and may result in their fracture. Vertebral augmentation is employed clinically to provide mechanical stability and pain relief for vertebrae with lytic lesions. However, little is known about its efficacy in strengthening fractured vertebrae containing lytic metastasis.

METHODS

Eighteen unembalmed human lumbar vertebrae, having simulated uncontained lytic defects and tested to failure in a prior study, were augmented using a transpedicular approach and re-tested to failure using a wedge fracture model. Axial and moment based strength and stiffness parameters were used to quantify the effect of augmentation on the structural response of the failed vertebrae. Effects of cement volume, bone mineral density and vertebral geometry on the change in structural response were investigated.

FINDINGS

Augmentation increased the failed lytic vertebral strength [compression: 85% (P<0.001), flexion: 80% (P<0.001), anterior-posterior shear: 95%, P<0.001)] and stiffness [(40% (P<0.05), 53% (P<0.05), 45% (P<0.05)]. Cement volume correlated with the compressive strength (r(2)=0.47, P<0.05) and anterior-posterior shear strength (r(2)=0.52, P<0.05) and stiffness (r(2)=0.45, P<0.05). Neither the geometry of the failed vertebrae nor its pre-fracture bone mineral density correlated with the volume of cement.

INTERPRETATION

Vertebral augmentation is effective in bolstering the failed lytic vertebrae compressive and axial structural competence, showing strength estimates up to 50-90% of historical values of osteoporotic vertebrae without lytic defects. This modest increase suggests that lytic vertebrae undergo a high degree of structural damage at failure, with strength only partially restored by vertebral augmentation. The positive effect of cement volume is self-limiting due to extravasation.

The role of patient-mode high-resolution peripheral quantitative computed tomography indices in the prediction of failure strength of the elderly women's thoracic vertebral body

The correlations between the failure load of 20 T12 vertebral bodies, their patient-mode high-resolution peripheral quantitative computed tomography (HR-pQCT) indices, and the L1 areal bone mineral density (aBMD) were investigated. For the prediction of the T12 vertebral failure load, the T12 HR-pQCT microarchitectural parameters added significant information to that of L1 aBMD and to that of cortical BMD, but not to that of T12 vertebral BMD and not to that of T12 trabecular BMD.

INTRODUCTION

HR-pQCT is a new in vivo imaging technique for assessing the three-dimensional microarchitecture of cortical and trabecular bone at the distal radius and tibia. But little is known about this technique in the direct measurement of vertebral body.

METHODS

Twenty female donors with the mean age of 80.1 (7.6) years were included in the study. Dual X-ray absorptiometry of the lumbar spine and femur was performed. The spinal specimens (T11/T12/L1) were dissected, scanned using HR-pQCT scanner, and mechanically tested under 4° wedge compression. The L1 aBMD, T12 patient-mode HR-pQCT indices, and T12 vertebral failure loads were analyzed.

RESULTS

For the prediction of vertebral failure load, the inclusion of BV/TV into L1 aBMD was the best model (R (2) = 0.52), Tb.N and Tb.Sp added significant information to the L1 aBMD and to the cortical BMD, but none of the vertebral microarchitectural parameters yielded additional significant information to the trabecular BMD (or BV/TV) and to the vertebral BMD.

CONCLUSION

Vertebral microarchitectural parameters obtained from the patient-mode HR-pQCT analysis provide significant information on bone strength complementary to that of aBMD and to that of cortical BMD, but not to that of vertebral BMD and not to that of trabecular BMD.