Effect of two-level decompressive procedures on the biomechanics of the lumbo-sacral spine: an ex vivo study

Hemilaminectomy and laminectomy are decompressive procedures commonly used in case of lumbar spinal stenosis, which involve the removal of the posterior elements of the spine. These procedures may compromise the stability of the spine segment and create critical strains in the intervertebral discs. Thus, this study aimed to investigate if decompressive procedures could alter the biomechanics of the lumbar spine. The focus was on the changes in the range of motion and strain distribution of the discs after two-level hemilaminectomy and laminectomy. Twelve L2-S1 cadaver specimens were prepared and mechanically tested in flexion, extension and both left and right lateral bending, in the intact condition, after a two-level hemilaminectomy on L4 and L5 vertebrae, and a full laminectomy. The range of motion (ROM) of the entire segment was assessed in all the conditions and loading configurations. In addition, Digital Image Correlation was used to measure the strain distribution on the surface of each specimen during the mechanical tests, focusing on the disc between the two decompressed vertebrae and in the two adjacent discs. Hemilaminectomy did not significantly affect the ROM, nor the strain on the discs. Laminectomy significantly increased the ROM in flexion, compared to the intact state. Laminectomy significantly increased the tensile strains on both L3-L4 and L4-L5 disc (p = 0.028 and p = 0.014) in ipsilateral bending, and the compressive strains on L4-L5 intervertebral disc, in both ipsilateral and contralateral bending (p = 0.014 and p = 0.0066), with respect to the intact condition. In conclusion, this study found out that hemilaminectomy did not significantly impact the biomechanics of the lumbar spine. Conversely, after the full laminectomy, flexion significantly increased the range of motion and lateral bending was the most critical configuration for largest principal strain.

Stereo Camera Setup for 360° Digital Image Correlation to Reveal Smart Structures of Hakea Fruits

About forty years after its first application, digital image correlation (DIC) has become an established method for measuring surface displacements and deformations of objects under stress. To date, DIC has been used in a variety of in vitro and in vivo studies to biomechanically characterise biological samples in order to reveal biomimetic principles. However, when surfaces of samples strongly deform or twist, they cannot be thoroughly traced. To overcome this challenge, different DIC setups have been developed to provide additional sensor perspectives and, thus, capture larger parts of an object's surface. Herein, we discuss current solutions for this multi-perspective DIC, and we present our own approach to a 360° DIC system based on a single stereo-camera setup. Using this setup, we are able to characterise the desiccation-driven opening mechanism of two woody Hakea fruits over their entire surfaces. Both the breaking mechanism and the actuation of the two valves in predominantly dead plant material are models for smart materials. Based on these results, an evaluation of the setup for 360° DIC regarding its use in deducing biomimetic principles is given. Furthermore, we propose a way to improve and apply the method for future measurements.

Applying a homogeneous pressure distribution to the upper vertebral endplate: Validation of a new loading system, pilot application to human vertebral bodies, and finite element predictions of DIC measured displacements and strains

Image-based personalized Finite Element Models (pFEM) could detect alterations in physiological deformation of human vertebral bodies, but their accuracy has been seldom reported. Meaningful validation experiments should allow vertebral endplate deformability and ensure well-controlled boundary conditions. This study aimed to (i) validate a new loading system to apply a homogeneous pressure on the vertebral endplate during vertebral body compression regardless of endplate deformation; (ii) perform a pilot study on human vertebral bodies measuring surface displacements and strains with Digital Image Correlation (DIC); (iii) determine the accuracy of pFEM of the vertebral bodies. Homogeneous pressure application was achieved by pressurizing a fluid silicone encased in a rubber silicone film acting on the cranial endplate. The loading system was validated by comparing DIC-measured longitudinal strains and lower-end contact pressures, measured on three homogeneous pseudovertebrae of constant transversal section at 2.0 kN, against theoretically calculated values. Longitudinal strains and contact pressures were rather homogeneous, and their mean values close to theoretical calculations (5% underestimation). DIC measurements of surface longitudinal and circumferential displacements and strains were obtained on three human vertebral bodies at 2.0 kN. Complete displacement and strain maps were achieved for anterolateral aspects with random errors ≤0.2 μm and ≤30 μstrain, respectively. Venous plexus and double curvatures limited the completeness and accuracy of DIC data in posterior aspects. pFEM of vertebral bodies, including cortical bone mapping, were built from computed tomography images. In anterolateral aspects, pFEM accuracy of the three vertebrae was: (i) comparable to literature in terms of longitudinal displacements (R2>0.8); (ii) extended to circumferential displacements (pooled data: R2>0.9) and longitudinal strains (zero median error, 95% error: <27%). Circumferential strains were overestimated (median error: 39%). The new methods presented may permit to study how physiological and pathologic conditions influence the ability of vertebral endplates/bodies to sustain loads.

Assessment of foraminal decompression following discoplasty using a combination of ex vivo testing and numerical tools

Percutaneous Cement Discoplasty (PCD) is a minimally invasive surgical technique to treat degenerated intervertebral discs. When the disc is severely degenerated, the vacuum observed in place of the nucleus pulposus can be filled with bone cement to restore the disc height, open the foramen space, and relieve pain. This study aimed to evaluate the foramen geometry change due to PCD, in the loaded spine. Cadaveric spines (n = 25) were tested in flexion and extension while Digital Image Correlation (DIC) measured displacements and deformations. Tests were performed on simulated pre-operative condition (nucleotomy) and after PCD. Registering DIC images and the 3D specimen geometry from CT scans, a 3D model of the specimens aligned in the experimental pose was obtained for nucleotomy and PCD. Foramen space volume was geometrically measured for both conditions. The volume of cement injected was measured to explore correlation with the change of foramen space. PCD induced a significant overall foraminal decompression in both flexion (foramen space increased by 835 ± 1289 mm³, p = 0.001) and extension (1205 ± 1106 mm³, p < 0.001), confirming that the expected improvements of PCD show also during spine motion. Furthermore, in extension when the foramen is the most challenged, the impact of PCD on the foramen correlated with the injected cement volume.

Biomechanical consequences of cement discoplasty: An in vitro study on thoraco-lumbar human spines

With the ageing of the population, there is an increasing need for minimally invasive spine surgeries to relieve pain and improve quality of life. Percutaneous Cement Discoplasty is a minimally invasive technique to treat advanced disc degeneration, including vacuum phenomenon. The present study aimed to develop an in vitro model of percutaneous cement discoplasty to investigate its consequences on the spine biomechanics in comparison with the degenerated condition. Human spinal segments (n = 27) were tested at 50% body weight in flexion and extension. Posterior disc height, range of motion, segment stiffness, and strains were measured using Digital Image Correlation. The cement distribution was also studied on CT scans. As main result, percutaneous cement discoplasty restored the posterior disc height by 41% for flexion and 35% for extension. Range of motion was significantly reduced only in flexion by 27%, and stiffness increased accordingly. The injected cement volume was 4.56 ± 1.78 ml (mean ± SD). Some specimens (n = 7) exhibited cement perforation of one endplate. The thickness of the cement mass moderately correlated with the posterior disc height and range of motion with different trends for flexions vs. extension. Finally, extreme strains on the discs were reduced by percutaneous cement discoplasty, with modified patterns of the distribution. To conclude, this study supported clinical observations in term of recovered disc height close to the foramen, while percutaneous cement discoplasty helped stabilize the spine in flexion and did not increase the risk of tissue damage in the annulus.

Subtle segmental angle changes of single-level lumbar fusions and adjacent-level biomechanics: cadaveric study of optically measured disc strain

OBJECTIVE

Changes to segmental lordosis at a single level may affect adjacent-level biomechanics and overall spinal alignment with an iatrogenic domino effect commonly seen in adult spinal deformity. This study investigated the effects of different segmental angles of single-level lumbar fixation on stability and principal strain across the surface of the adjacent-level disc.

METHODS

Seven human cadaveric L3-S1 specimens were instrumented at L4-5 and tested in 3 conditions: 1) neutral native angle ("neutral"), 2) increasing angle by 5° of lordosis ("lordosis"), and 3) decreasing angle by 5° of kyphosis ("kyphosis"). Pure moment loads (7.5 Nm) were applied in flexion, extension, lateral bending, and axial rotation, followed by 400 N of axial compression alone and together with pure moments. Range of motion (ROM), principal maximum strain (E1), and principal minimum strain (E2) across different surface subregions of the upper adjacent-level disc (L3-4) were optically assessed. Larger magnitudes of either E1 or E2 indicate larger tissue deformations and represent indirect measures of increased stress.

RESULTS

At the superior adjacent level, a significant increase in ROM was observed in kyphosis and lordosis versus neutral in flexion (p ≤ 0.001) and extension (p ≤ 0.02). ROM was increased in lordosis versus neutral (p = 0.03) and kyphosis (p = 0.004) during compression. ROM increased in kyphosis versus neutral and lordosis (both p = 0.03) in compression plus extension. Lordosis resulted in increased E1 across the midposterior subregion of the disc (Q3) versus neutral during right lateral bending (p = 0.04); lordosis and kyphosis resulted in decreased E1 in Q3 versus neutral with compression (p ≤ 0.03). Lordosis decreased E1 in Q3 versus neutral during compression plus flexion (p = 0.01), whereas kyphosis increased E1 in all quartiles and increased E2 in the midanterior subregion versus lordosis in compression plus flexion (p ≤ 0.047). Kyphosis decreased E1 in Q3 (p = 0.02) and E2 in the anterior-most subregion of the disc (Q1) (p = 0.006) versus neutral, whereas lordosis decreased E1 in Q3 (p = 0.008) versus neutral in compression plus extension.

CONCLUSIONS

Lumbar spine monosegmental fixation with 5° offset from the neutral individual segmental angle altered the motion and principal strain magnitudes at the upper adjacent disc, with induced kyphosis resulting in larger principal strains compared with lordosis. Segmental alignment of single-level fusion influences adjacent-segment biomechanics, and suboptimal alignment may play a role in the clinical development of adjacent-segment disease.

Experimental validation of a subject-specific finite element model of lumbar spine segment using digital image correlation

Pathologies such as cancer metastasis and osteoporosis strongly affect the mechanical properties of the vertebral bone and increase the risk of fragility fractures. The prediction of the fracture risk with a patient-specific model, directly generated from the diagnostic images of the patient, could help the clinician in the choice of the correct therapy to follow. But before such models can be used to support any clinical decision, their credibility must be demonstrated through verification, validation, and uncertainty quantification. In this study we describe a procedure for the generation of such patient-specific finite element models and present a first validation of the kinematics of the spine segment. Quantitative computed tomography images of a cadaveric lumbar spine segment presenting vertebral metastatic lesions were used to generate the model. The applied boundary conditions replicated a specific experimental test where the spine segment was loaded in compression-flexion. Model predictions in terms of vertebral surface displacements were compared against the full-field experimental displacements measured with Digital Image Correlation. A good agreement was obtained from the local comparison between experimental data and simulation results (R² > 0.9 and RMSE% <8%). In conclusion, this work demonstrates the possibility to apply the developed modelling pipeline to predict the displacement field of human spine segment under physiological loading conditions, which is a first fundamental step in the credibility assessment of these clinical decision-support technology.

Use of digital imaging correlation techniques for full-field strain distribution analysis of implantable devices and tissue in spinal biomechanics research

Few studies have used optical full-field surface strain mapping to study spinal biomechanics. We used a commercial digital imaging correlation (DIC) system to (1) compare posterior surface strains on spinal rods with those obtained from conventional foil strain gauges, (2) quantify bony vertebral body and intervertebral disc (IVD) surface strains on 3 L3-S cadaveric spines during gold-standard flexibility tests (7.5-Nm flexion-extension and 400-N compression), and (3) report our experience with the application and feasibility of DIC to comprehensively map strain in spinal biomechanics. Spinal rods were tested under zero load and using ASTM F1717 standard. For rod strain measures, the largest mean bias offset and baseline noise standard deviation under zero load for DIC were 7.6 με and 33.7 με, respectively. For tissue measures, the largest mean bias offset was 8 με for ε1 and -55 με for ε2 with baseline noise standard deviations of 19 με and 26 με, respectively. On average, DIC rod strain measurements were 5.3% less than strain gauge measurements throughout the load range. Principal IVD and bony surface strains were consistently measurable and showed marked regional differences in strain patterns under different load conditions. Strains measured on spinal rods using DIC techniques reasonably agreed with standard strain gauge measurements. Subregional strain analyses on soft and hard spinal tissues during standard flexibility tests were feasible. Optical strain mapping is a viable, accurate, and promising measurement technique for novel spinal biomechanical studies.

Global and local characterization explains the different mechanisms of failure of the human ribs

Knowledge of the mechanics and mechanistic reasons inducing rib fracture is fundamental for forensic investigations and for the design of implants and cardiopulmonary resuscitation devices. A mechanical rationale to explain the different rib mechanisms of failure is still a challenge. The aim of this work was to experimentally characterize human ribs to test the hypothesis that a correlation exists between the ribs properties and the mechanism of failure. 89 ribs were tested in antero-posterior compression. The full-field strain distribution was measured through Digital Image Correlation. The fracture load ranged 7-132 N. Two main different mechanisms of failure were observed: brittle and buckling. The strain analysis showed that the direction of principal strains was either aligned with the ribs, or oblique, around 45°, with a rather uniform direction in the most strained area. The maximum principal strains were in the range between 1000 and 30000 microstrain and the minimum principal strain between -30000 and -800 microstrain. The ribs undergoing brittle fracture had significantly thicker cortical bone than those undergoing buckling. Also, larger tensile strains were observed in the specimens with brittle fracture than in the buckling ones. These findings support the focus of cortical thickness modelling which could help in sharpening computational models for the aforesaid purposes.

Type, size, and position of metastatic lesions explain the deformation of the vertebrae under complex loading conditions

BACKGROUND

Bone metastases may lead to spine instability and increase the risk of fracture. Scoring systems are available to assess critical metastases, but they lack specificity, and provide uncertain indications over a wide range, where most cases fall. The aim of this work was to use a novel biomechanical approach to evaluate the effect of lesion type, size, and location on the deformation of the metastatic vertebra.

METHOD

Vertebrae with metastases were identified from 16 human spines from a donation programme. The size and position of the metastases, and the Spine Instability Neoplastic Score (SINS) were evaluated from clinical Quantitative Computed Tomography images. Thirty-five spine segments consisting of metastatic vertebrae and adjacent healthy controls were biomechanically tested in four different loading conditions. The strain distribution over the entire vertebral bodies was measured with Digital Image Correlation. Correlations between the features of the metastasis (type, size, position and SINS) and the deformation of the metastatic vertebrae were statistically explored.

RESULTS

The metastatic type (lytic, blastic, mixed) characterizes the vertebral behaviour (Kruskal-Wallis, p = 0.04). In fact, the lytic metastases showed more critical deformation compared to the control vertebrae (average: 2-fold increase, with peaks of 14-fold increase). By contrast, the vertebrae with mixed or blastic metastases did not show a clear trend, with deformations similar or lower than the controls. Once the position of the lytic lesion with respect to the loading direction was taken into account, the size of the lesion was significantly correlated with the perturbation to the strain distribution (r² = 0.72, p < 0.001). Conversely, the SINS poorly correlated with the mechanical evidence, and only in case of lytic lesions (r² = 0.25, p < 0.0001).

CONCLUSION

These results highlight the relevance of the size and location of the lytic lesion, which are marginally considered in the current clinical scoring systems, in driving the spinal biomechanical instability. The strong correlation with the biomechanical evidence indicates that these parameters are representative of the mechanical competence of the vertebra. The improved explanatory power compared to the SINS suggests including them in future guidelines for the clinical practice.

A novel specimen shape for measurement of linear strain fields by means of digital image correlation

Strains on the surface of engineering structures or biological tissues are non-homogeneous. These strain fields can be captured by means of Digital Image Correlation (DIC). However, DIC strain field measurements are prone to noise and filtering of these fields influences measured strain gradients. This study aims to design a novel tensile test specimen showing two linear gradients, to measure full-field linear strain measurements on the surface of test specimens, and to investigate the accuracy of DIC strain measurements globally (full-field) and locally (strain gauges' positions), with and without filtering of the DIC strain fields. Three materials were employed for this study: aluminium, polymer, and bovine bone. Normalized strain gradients were introduced that are load independent and evaluated at two local positions showing 3.6 and 6.9% strain change per mm. Such levels are typically found in human bones. At these two positions, two strain gauges were applied to check the experimental strain magnitudes. A third strain gauge was applied to measure the strain in a neutral position showing no gradient. The accuracy of the DIC field measurement was evaluated at two deformation stages (at [Formula: see text] 500 and 1750 μstrain) using the root mean square error (RMSE). The RMSE over the two linear strain fields was less than 500 μstrain for both deformation stages and all materials. Gaussian low-pass filter (LPF) reduced the DIC noise between 25% and 64% on average. As well, filtering improved the accuracy of the local normalized strain gradients measurements with relative difference less than 20% and 12% for the high- and low-gradient, respectively. In summary, a novel specimen shape and methodological approach are presented which are useful for evaluating and improving the accuracy of the DIC measurement where non-homogeneous strain fields are expected such as on bone tissue due to their hierarchical structure.

Digital image correlation in dental materials and related research: A review

OBJECTIVE

Digital image correlation (DIC) is a non-contact image processing technique for full-field strain measurement. Although DIC has been widely used in engineering and biomechanical fields, it is in the spotlight only recently in dental materials. Therefore, the purpose of this review paper is introducing the working principle of the DIC technique with some modifications and providing further potential applications in various dental materials and related fields.

METHODS

The accuracy of the algorithm depending on the environmental characteristics of the DIC technique, as well as the advantages and disadvantages of strain measurement using optical measurements, have been elaborated in dental materials and related fields. Applications to those researches have been classified into the following categories: shrinkage behavior of light-cured resin composite, resin-tooth interface, mechanical properties of tooth structure, crack extension and elastic properties of dental materials, and deformation of dental restoration and prosthesis. This classification and discussion were performed using literature survey and review based on numerous papers in the international journals published over the past 20 years. The future directions for predicting the precise deformation of dental materials under various environments, as well as limitations of the DIC technique, was presented in this review.

RESULTS

The DIC technique was demonstrated as a more effective tool to measure full-field polymerization shrinkage of composite resin, even in a simulated clinical condition over the existing methods. Moreover, the DIC combined with other technologies can be useful to evaluate the mechanical behavior of material-tooth interface, dentine structure and restorative and prosthetic materials with high accuracy. Three-dimensional DIC using two cameras extended the measurement range in-plane to out-of-plane, enabling measure of the strain directly on the surface of dental restorations or prosthesis.

SIGNIFICANCE

DIC technique is a potential tool for measuring and predicting the full-field deformation/strain of dental materials and actual prostheses in diverse clinical conditions. The versatility of DIC can replace the existing complex sensor devices in those studies.

Investigation by Digital Image Correlation of Mixed-Mode I and II Fracture Behavior of Polymeric IASCB Specimens with Additive Manufactured Crack-Like Notch

In this work, the fracture mechanics properties of polyamide (PA) specimens manufactured by the selective laser sintering (SLS) technology are investigated, in which an embedded crack-like notch was inserted in the design and produced during the additive manufacturing (AM) phase. To cover a wide variety of mode I/II mixity levels, the inclined asymmetrical semicircular specimen subjected to three points loading (IASCB) was employed. The investigation was carried out by analyzing the full displacement field in the proximity of the crack tip by means of the digital image correlation (DIC) technique. To characterize the material, which exhibits a marked elastic-plastic behavior, the quantity J-integral was evaluated by two different methods: the first one exploits the full fields measured by the DIC, whereas the second one exploits the experimental load-displacement curves along with FEM analysis. The DIC methodology was experimentally validated and proposed as an alternative method to evaluate the J-integral. It is especially suited for conditions in which it is not possible to use the conventional LDC method due to complex and possibly unknown loading conditions. Furthermore, results showed that the AM technique could be used effectively to induce cracks in this type of material. These two aspects together can lead to both a simplification of the fracture characterization process and to the possibility of dealing with a wider number of practical, real-world scenarios. Indeed, because of the nature of the additive manufacturing process, AM crack-like notches can be sintered even having complex geometry, being three-dimensional and/or inside the tested structure.

Load-sharing biomechanics of lumbar fixation and fusion with pedicle subtraction osteotomy

Pedicle subtraction osteotomy (PSO) is an invasive surgical technique allowing the restoration of a well-balanced sagittal profile, however, the risks of pseudarthrosis and instrumentation breakage are still high. Literature studied primary stability and posterior instrumentation loads, neglecting the load shared by the anterior column, which is fundamental to promote fusion early after surgery. The study aimed at quantifying the load-sharing occurring after PSO procedure across the ventral spinal structures and the posterior instrumentation, as affected by simple bilateral fixation alone, with interbody cages adjacent to PSO level and supplementary accessory rods. Lumbar spine segments were loaded in vitro under flexion–extension, lateral bending, and torsion using an established spine tester. Digital image correlation (DIC) and strain-gauge (SG) analyses measured, respectively, the full-field strain distribution on the ventral surface of the spine and the local strain on posterior primary rods. Ventral strains considerably decreased following PSO and instrumentation, confirming the effectiveness of posterior load-sharing. Supplemental accessory rods considerably reduced the posterior rod strains only with interbody cages, but the ventral strains were unaffected: this indicates that the load transfer across the osteotomy could be promoted, thus explaining the higher fusion rate with decreased rod fracture risk reported in clinical literature.

The Effect of the Loading Rate on the Full-Field Strain Distribution on the Surface on the Intervertebral Discs

Contrasting results are reported when the spine is tested at different strain rates. Tissue specimens from the ligaments or the intervertebral discs (IVD, including annulus fibrosus and nucleus pulposus) exhibit higher stiffness and lower dissipation at high strain rates. Counterintuitively, when spine segments are tested at high rates, the hysteresis area and loop width increase. It is unclear how the load is shared between the different structures at different loading rates. The hypotheses of this study were: (i) As the IVD stiffens at higher loading rates, the strain distribution around the disc would be different depending on the loading rate; (ii) Preconditioning attenuates the strain-rate dependency of the IVD, thus making differences in strain distribution smaller at the different rates. Six segments of three vertebrae (L4-L6) were extracted from porcine spines and tested in presso-flexion at different loading rates (reaching full load in 0.67, 6.7, and 67 s). The full-field strain maps were measured using digital image correlation on the surface of the IVDs from lateral. The posterior-to-anterior trends of the strain were computed in detail for each IVD, and compared between loading rates. The values and the direction of principal strain on the surface of the IVDs, vertebrae, and endplates remained unchanged at different rates. In the transition zone between IVD and vertebra, only slight differences due to the loading rate appeared but with no statistical significance. These findings will allow better understanding of the rate-dependent behavior and failure of the IVD.

Properties of PMMA end cap holders affect FE stiffness predictions of vertebral specimens

Bone cement is often used, in experimental biomechanics, as a potting agent for vertebral bodies (VB). As a consequence, it is usually included in finite element (FE) models to improve accuracy in boundary condition settings. However, bone cement material properties are typically assigned to these models based on literature data obtained from specimens created under conditions which often differ from those employed for cement end caps. These discrepancies can result in solids with different material properties from those reported. Therefore, this study aimed to analyse the effect of assigning different mechanical properties to bone cement in FE vertebral models. A porcine C2 vertebral body was potted in bone cement end caps, μCT scanned, and tested in compression. DIC was performed on the anterior surface of the specimen to monitor the displacement. Specimen stiffness was calculated from the load-displacement output of the materials testing machine and from the machine load output and average displacement measured by DIC. Fifteen bone cement cylinders with dimensions similar to the cement end caps were produced and subjected to the same compression protocol as the vertebral specimen and average stiffness and Young moduli were estimated. Two geometrically identical vertebral body FE models were created from the μCT images, the only difference residing in the values assigned to bone cement material properties: in one model these were obtained from the literature and in the other from the cylindrical cement samples previously tested. The average Youngs modulus of the bone cement cylindrical specimens was 1177 ± 3 MPa, considerably lower than the values reported in the literature. With this value, the FE model predicted a vertebral specimen stiffness 3% lower than that measured experimentally, while when using the value most commonly reported in similar studies, specimen stiffness was overestimated by 150%.

Body Anthropometry and Bone Strength Conjointly Determine the Risk of Hip Fracture in a Sideways Fall

We hypothesize that variations of body anthropometry, conjointly with the bone strength, determine the risk of hip fracture. To test the hypothesis, we compared, in a simulated sideways fall, the hip impact energy to the energy needed to fracture the femur. Ten femurs from elderly donors were tested using a novel drop-tower protocol for replicating the hip fracture dynamics during a fall on the side. The impact energy was varied for each femur according to the donor’s body weight, height and soft-tissue thickness, by adjusting the drop height and mass. The fracture pattern, force, energy, strain in the superior femoral neck, bone morphology and microarchitecture were evaluated. Fracture patterns were consistent with clinically relevant hip fractures, and the superior neck strains and timings were comparable with the literature. The hip impact energy (11 – 95 J) and the fracture energy (11 – 39 J) ranges overlapped and showed comparable variance (CV = 69 and 61%, respectively). The aBMD-based definition of osteoporosis correctly classified 7 (70%) fracture/non-fracture cases. The incorrectly classified cases presented large impact energy variations, morphology variations and large subcortical voids as seen in microcomputed tomography. In conclusion, the risk of osteoporotic hip fracture in a sideways fall depends on both body anthropometry and bone strength.

Testing the impact of discoplasty on the biomechanics of the intervertebral disc with simulated degeneration: An in vitro study

Percutaneous Cement Discoplasty has recently been developed to relieve pain in highly degenerated intervertebral discs presenting a vacuum phenomenon in patients that cannot undergo major surgery. Little is currently known about the biomechanical effects of discoplasty. This study aimed at investigating the feasibility of modelling empty discs and subsequent discoplasty surgery and measuring their impact over the specimen geometry and mechanical behaviour. Ten porcine lumbar spine segments were tested in flexion, extension, and lateral bending under 5.4 Nm (with a 200 N compressive force and a 27 mm offset). Tests were performed in three conditions for each specimen: with intact disc, after nucleotomy and after discoplasty. A 3D Digital Image Correlation (DIC) system was used to measure the surface displacements and strains. The posterior disc height, range of motion (ROM), and stiffness were measured at the peak load. CT scans were performed to confirm that the cement distribution was acceptable. Discoplasty recovered the height loss caused by nucleotomy (p = 0.04) with respect to the intact condition, but it did not impact significantly either the ROM or the stiffness. The strains over the disc surface increased after nucleotomy, while discoplasty concentrated the strains on the endplates. In conclusion, this preliminary study has shown that discoplasty recovered the intervertebral posterior height, opening the neuroforamen as clinically observed, but it did not influence the spine mobility or stiffness. This study confirms that this in vitro approach can be used to investigate discoplasty.

Assessing the Mechanical Weakness of Vertebrae Affected by Primary Tumors: A Feasibility Study

Patients spend months between the primary spinal tumor diagnosis and the surgical treatment, due to the need for performing chemotherapy and/or radiotherapy. During this period, they are exposed to an unknown risk of fracture. The aim of this study was to assess if it is possible to measure the mechanical strain in vertebrae affected by primary tumors, so as to open the way to an evidence-based scoring or prediction tool. We performed biomechanical tests on three vertebrae with bone tumor removed from patients. The tests were designed so as not to compromise the standard surgical and diagnostic procedures. Non-destructive mechanical tests in combination with state-of-the-art digital image correlation allowed to measure the distribution of strain on the surface of the vertebra. Our study has shown that the strains in the tumor region is circa 3 times higher than in the healthy bones, with principal strain peaks of 40,000/-20,000 microstrain, indicating a stress concentration potentially triggering vertebral fracture. This study has proven it is possible to analyze the mechanical behavior of primary tumor vertebrae as part of the clinical treatment protocol. This will allow building a tool for quantifying the risk of fracture and improving decision making in spine tumors.

Investigation by Digital Image Correlation of Mixed Mode I and II Fracture Behavior of Metallic IASCB Specimens with Additive Manufactured Crack-Like Notch

This work investigates the fracture behavior of maraging steel specimens manufactured by the selective laser sintering (SLS) technology, in which a crack-like notch (sharp notch) was directly produced during the additive manufacturing (AM) process. For the evaluation of the fracture toughness, the inclined asymmetrical semi-circular specimen subjected to three points loading (IASCB) was used, allowing to cover a wide variety of Mode I and II combinations. The effectiveness of manufacturing crack-like notches via the SLS technique in metals was evaluated by comparing the obtained experimental results with the ones obtained with pre-cracks induced by fatigue loading. The investigation was carried out by using the digital image correlation (DIC) technique, that allowed the evaluation of the full displacement fields around the crack tip. The displacement field was then used to compute the stress intensity factors (SIFs) for various combinations of Mode I and II, via a fitting technique which relies on the Williams’ model for the displacement. The SIFs obtained in this way were compared to the results obtained with the conventional critical load method. The results showed that the discrepancy between the two methods reduces by ranging from Mode I to Mode II loading condition. Finally, the experimental SIFs obtained by the two methods were described by the mixed mode local stress criterium.

Digital Image Correlation (DIC) Assessment of the Non-Linear Response of the Anterior Longitudinal Ligament of the Spine during Flexion and Extension

While the non-linear behavior of spine segments has been extensively investigated in the past, the behavior of the Anterior Longitudinal Ligament (ALL) and its contribution during flexion and extension has never been studied considering the spine as a whole. The aims of the present study were to exploit Digital Image Correlation (DIC) to: (I) characterize the strain distribution on the ALL during flexion-extension, (II) compare the strain on specific regions of interest (ROI) of the ALL in front of the vertebra and of the intervertebral disc, (III) analyze the non-linear relationship between the surface strain and the imposed rotation and the resultant moment. Three specimens consisting of 6 functional spinal units (FSUs) were tested in flexion-extension. The full-field strain maps were measured on the surface of the ALL, and the most strained areas were investigated in detail. The DIC-measured strains showed different values of peak strain in correspondence with the vertebra and the disc but the average over the ROIs was of the same order of magnitude. The strain-moment curves showed a non-linear response like the moment-angle curves: in flexion the slope of the strain-moment curve was greater than in extension and with a more abrupt change of slope. To the authors’ knowledge, this is the first study addressing, by means of a full-field strain measurement, the non-linear contribution of the ALL to spine biomechanics. This study was limited to only three specimens; hence the results must be taken with caution. This information could be used in the future to build more realistic numerical models of the spine.

The strain distribution in the lumbar anterior longitudinal ligament is affected by the loading condition and bony features: An in vitro full-field analysis

The role of the ligaments is fundamental in determining the spine biomechanics in physiological and pathological conditions. The anterior longitudinal ligament (ALL) is fundamental in constraining motions especially in the sagittal plane. The ALL also confines the intervertebral discs, preventing herniation. The specific contribution of the ALL has indirectly been investigated in the past as a part of whole spine segments where the structural flexibility was measured. The mechanical properties of isolated ALL have been measured as well. The strain distribution in the ALL has never been measured under pseudo-physiological conditions, as part of multi-vertebra spine segments. This would help elucidate the biomechanical function of the ALL. The aim of this study was to investigate in depth the biomechanical function of the ALL in front of the lumbar vertebrae and of the intervertebral disc. Five lumbar cadaveric spine specimens were subjected to different loading scenarios (flexion-extension, lateral bending, axial torsion) using a state-of-the-art spine tester. The full-field strain distribution on the anterior surface was measured using digital image correlation (DIC) adapted and validated for application to spine segments. The measured strain maps were highly inhomogeneous: the ALL was generally more strained in front of the discs than in front of the vertebrae, with some locally higher strains both imputable to ligament fibers and related to local bony defects. The strain distributions were significantly different among the loading configurations, but also between opposite directions of loading (flexion vs. extension, right vs. left lateral bending, clockwise vs. counterclockwise torsion). This study allowed for the first time to assess the biomechanical behaviour of the anterior longitudinal ligament for the different loading of the spine. We were able to identify both the average trends, and the local effects related to osteophytes, a key feature indicative of spine degeneration.

Effect of different motor tasks on hip cup primary stability and on the strains in the periacetabular bone: An in vitro study

BACKGROUND

Excessive prosthesis/bone motions and the bone strains around the acetabulum may prevent osteointegration and lead to cup loosening. These two factors depend on post-operative joint loading. We investigated how Walking (which is often simulated) and Standing-Up from seated (possibly more critical) influence the cup primary stability and periacetabular strains.

METHODS

Twelve composite hemipelvises were used in two test campaigns. Simplified loading conditions were adopted to simulate Walking and Standing-Up. For each motor task, a single-direction force was applied in load packages of increasing amplitude. Stable and unstable uncemented cups were implanted. Digital image correlation was used to measure implant/bone motions (three-dimensional translations and rotations, both permanent and inducible), and the strain distribution around the acetabulum.

FINDINGS

When stable implants were tested, higher permanent cranial translations were found during Walking (however the resultant migrations were comparable with Standing-Up); higher rotations were found for Standing-Up. When unstable implants were tested, motions were 1-2 order of magnitude higher. Strains increased significantly from stable to unstable implants. The peak strains were in the superior aspect of the acetabulum during Walking and in the superior-posterior aspect of the acetabulum and at the bottom of the posterior column during Standing-Up.

INTERPRETATION

Different cup migration trends were caused by simulated Walking and Standing-Up, both similar to those observed clinically. The cup mobilization pattern depended on the different simulated motor tasks. Pre-clinical testing of new uncemented cups could include simulation of both motor tasks. Our study could also translate to indication of what tasks should be avoided.

The Size of Simulated Lytic Metastases Affects the Strain Distribution on the Anterior Surface of the Vertebra

Metastatic lesions of the vertebra are associated with risk of fracture, which can be disabling and life-threatening. In the literature, attempts are found to identify the parameters that reduce the strength of a metastatic vertebra leading to spine instability. However, a number of controversial issues remain. Our aim was to quantify how the strain distribution in the vertebral body is affected by the presence and by the size of a simulated metastatic defect. Five cadaveric thoracic spine segments were subjected to non-destructive presso-flexion while intact, and after simulation of metastases of increasing size. For the largest defect, the specimens were eventually tested to failure. The full-field strain distribution in the elastic range was measured with digital image correlation (DIC) on the anterior surface of the vertebral body. The mean strain in the vertebra remained similar to the intact when the defects were smaller than 30% of the vertebral volume. The mean strains became significantly larger than in the intact for larger defects. The map of strain and its statistical distribution indicated a rather uniform condition in the intact vertebra and with defects smaller than 30%. Conversely, the strain distribution became significantly different from the intact for defects larger than 30%. A strain peak appeared in the region of the simulated metastasis, where fracture initiated during the final destructive test. This is a first step in understanding how the features of metastasis influence the vertebral strain and for the construction of a mechanistic model to predicted fracture.

High resolution three-dimensional strain mapping of bioprosthetic heart valves using digital image correlation

Transcatheter aortic valve replacement (TAVR) is a safe and effective treatment option for patients deemed at high and intermediate risk for surgical aortic valve replacement. Similar to surgical aortic valves (SAVs), transcatheter aortic valves (TAVs) undergo calcification and mechanical wear over time. However, to date, there have been limited publications on the long-term durability of TAV devices. To assess longevity and mechanical strength of TAVs in comparison to surgical bioprosthetic valves, three-dimensional deformation analysis and strain measurement of the leaflets become an inevitable part of the evaluation. The goal of this study was to measure and compare leaflet displacement and strain of two commonly used TAVs in a side-by-side comparison with a commonly used SAV using a high-resolution digital image correlation (DIC) system. 26-mm Edwards SAPIEN 3, 26-mm Medtronic CoreValve, and 25-mm Carpentier-Edwards PERIMOUNT Magna surgical bioprosthesis were examined in a custom-made valve testing apparatus. A time-varying, spatially uniform pressure was applied to the leaflets at different loading rates. GOM ARAMIS® software was used to map leaflet displacement and strain fields during loading and unloading. High displacement regions were found to be at the leaflet belly region of the three bioprosthetic valves. In addition, the frame of the surgical bioprosthesis was found to be remarkably flexible, in contrary to CoreValve and SAPIEN 3 in which the stent was nearly rigid under a similar loading condition. The experimental DIC measurements can be used to characterize the anisotropic materiel behavior of the bioprosthetic heart valve leaflets and validate heart valve computational simulations.

Neighborhood binary speckle pattern for deformation measurements insensitive to local illumination variation by digital image correlation

Speckle pattern-based characteristics of digital image correlation (DIC) restrict its application in engineering fields and nonlaboratory environments, since serious decorrelation effect occurs due to localized sudden illumination variation. A simple and efficient speckle pattern adjusting and optimizing approach presented in this paper is aimed at providing a novel speckle pattern robust enough to resist local illumination variation. The new speckle pattern, called neighborhood binary speckle pattern, derived from original speckle pattern, is obtained by means of thresholding the pixels of a neighborhood at its central pixel value and considering the result as a binary number. The efficiency of the proposed speckle pattern is evaluated in six experimental scenarios. Experiment results indicate that the DIC measurements based on neighborhood binary speckle pattern are able to provide reliable and accurate results, even though local brightness and contrast of the deformed images have been seriously changed. It is expected that the new speckle pattern will have more potential value in engineering applications.

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.