Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

Mechanical influence of facet tropism in patients with chronic discogenic pain disorder

Aims

The presence of facet tropism has been correlated with an elevated susceptibility to lumbar disc pathology. Our objective was to evaluate the impact of facet tropism on chronic lumbosacral discogenic pain through the analysis of clinical data and finite element modelling (FEM).

Methods

Retrospective analysis was conducted on clinical data, with a specific focus on the spinal units displaying facet tropism, utilizing FEM analysis for motion simulation. We studied 318 intervertebral levels in 156 patients who had undergone provocation discography. Significant predictors of clinical findings were identified by univariate and multivariate analyses. Loading conditions were applied in FEM simulations to mimic biomechanical effects on intervertebral discs, focusing on maximal displacement and intradiscal pressures, gauged through alterations in disc morphology and physical stress.

Results

A total of 144 discs were categorized as 'positive' and 174 discs as 'negative' by the results of provocation discography. The presence of defined facet tropism (OR 3.451, 95% CI 1.944 to 6.126) and higher Adams classification (OR 2.172, 95% CI 1.523 to 3.097) were important predictive parameters for discography-'positive' discs. FEM simulations showcased uneven stress distribution and significant disc displacement in tropism-affected discs, where loading exacerbated stress on facets with greater angles. During varied positions, notably increased stress and displacement were observed in discs with tropism compared to those with normal facet structure.

Conclusion

Our findings indicate that facet tropism can contribute to disc herniation and changes in intradiscal pressure, potentially exacerbating disc degeneration due to altered force distribution and increased mechanical stress.

Comparative FEM study on intervertebral disc modeling: Holzapfel-Gasser-Ogden vs. structural rebars

Introduction: Numerical modeling of the intervertebral disc (IVD) is challenging due to its complex and heterogeneous structure, requiring careful selection of constitutive models and material properties. A critical aspect of such modeling is the representation of annulus fibers, which significantly impact IVD biomechanics. This study presents a comparative analysis of different methods for fiber reinforcement in the annulus fibrosus of a finite element (FE) model of the human IVD. Methods: We utilized a reconstructed L4-L5 IVD geometry to compare three fiber modeling approaches: the anisotropic Holzapfel-Gasser-Ogden (HGO) model (HGO fiber model) and two sets of structural rebar elements with linear-elastic (linear rebar model) and hyperelastic (nonlinear rebar model) material definitions, respectively. Prior to calibration, we conducted a sensitivity analysis to identify the most important model parameters to be calibrated and improve the efficiency of the calibration. Calibration was performed using a genetic algorithm and in vitro range of motion (RoM) data from a published study with eight specimens tested under four loading scenarios. For validation, intradiscal pressure (IDP) measurements from the same study were used, along with additional RoM data from a separate publication involving five specimens subjected to four different loading conditions. Results: The sensitivity analysis revealed that most parameters, except for the Poisson ratio of the annulus fibers and C01 from the nucleus, significantly affected the RoM and IDP outcomes. Upon calibration, the HGO fiber model demonstrated the highest accuracy (R2 = 0.95), followed by the linear (R2 = 0.89) and nonlinear rebar models (R2 = 0.87). During the validation phase, the HGO fiber model maintained its high accuracy (RoM R2 = 0.85; IDP R2 = 0.87), while the linear and nonlinear rebar models had lower validation scores (RoM R2 = 0.71 and 0.69; IDP R2 = 0.86 and 0.8, respectively). Discussion: The results of the study demonstrate a successful calibration process that established good agreement with experimental data. Based on our findings, the HGO fiber model appears to be a more suitable option for accurate IVD FE modeling considering its higher fidelity in simulation results and computational efficiency.

Effect of patient-specific factors on regeneration in lumbar spine at healthy disc and total disc replacement. Computer simulation

BACKGROUND AND OBJECTIVE

Degenerative diseases of the spine have a negative impact on the quality of life of patients. This study presents the results of numerical modelling of the mechanical behaviour of the lumbar spine with patient-specific conditions at physiological loads. This paper aims to numerically study the influence of degenerative changes in the spine and the presence of an endoprosthesis on the creation of conditions for tissue regeneration.

METHODS

A numerical model of the mechanical behaviour of lumbar spine at healthy and after total disc replacement under low-energy impacts equivalent to physiological loads is presented. The model is based on the movable cellular automaton method (discrete elements), where the mechanical behaviour of bone tissue is described using the Biot poroelasticity accounting for the presence and transfer of interstitial biological fluid. The nutritional pathways of the intervertebral disc in cases of healthy and osteoporotic bone tissues were predicted based on the analysis of the simulation results according to the mechanobiological principles.

RESULTS

Simulation of total disc replacement showed that osseointegration of the artificial disc plates occurs only in healthy bone tissue. With total disc replacement in a patient with osteoporosis, there is an area of increased risk of bone resorption in the near-contact area, approximately 1 mm wide, around the fixators. Dynamic loads may improve the osseointegration of the implant in pathological conditions of the bone tissue.

CONCLUSIONS

The results obtained in the case of healthy spine and osteoporotic bone tissues correspond to the experimental data on biomechanics and possible methods of IVD regeneration from the position of mechanobiological principles. The results obtained with an artificial disc (with keel-type fixation) showed that the use of this type of endoprosthesis in healthy bone tissues allows to reproduce the function of the natural intervertebral disc and does not contribute to the development of neoplastic processes. In the case of an artificial disc with osteoporosis of bone tissues, there is a zone with increased risk of tissue resorption and development of neoplastic processes in the area near the contact of the implant attachment. This circumstance can be compensated by increasing the loading level.

Effect of novel polyethylene insert configurations on bone-implant micromotion and contact stresses in total ankle replacement prostheses: a finite element analysis

Purpose

This study investigates the impact of elastic improvements to the artificial ankle joint insert on prosthesis biomechanics to reduce the risk of prosthesis loosening in TAR patients.

Methods

CT data of the right ankle was collected from one elderly female volunteer. An original TAR model (Model A) was developed from CT images and the INBONE II implant system. The development of the new inserts adopts an elastic improvement design approach, where different geometric configurations of flexible layers are inserted into the traditional insert. The structure can be divided into continuous flexible layers and intermittent flexible layers. The flexible layers aim to improve the elasticity of the component by absorbing and dispersing more kinetic energy. The newly designed inserts are used to replace the original insert in Model A, resulting in the development of Models B-D. A finite element model of gait analysis was based by gait parameters. Discrepancies in micromotion and contact behaviour were analysed during the gait cycle, along with interface fretting and articular surface stress at 50% of the gait cycle.

Results

In terms of micromotion, the improved elastic models showed reduced micromotion at the tibial-implant interfaces compared to the original model. The peak average micromotion decreased by 12.1%, 13.1%, and 14.5% in Models B, C, and D, respectively. The micromotion distribution also improved in the improved models, especially in Model D. Regarding contact areas, all models showed increased contact areas of articular surfaces with axial load, with Models B, C, and D increasing by 26.8%, 23.9%, and 24.4%, respectively. Contact stress on articular surfaces increased with axial load, reaching peak stress during the late stance phase. Models with continuous flexible layer designs exhibited lower stress levels. The insert and the talar prosthetic articular surfaces showed more uniform stress distribution in the improved models.

Conclusion

Improving the elasticity of the insert can enhance component flexibility, absorb impact forces, reduce micromotion, and improve contact behavior. The design scheme of continuous flexible layers is more advantageous in transmitting and dispersing stress, providing reference value for insert improvement.

Development and evaluation of sex-specific thoracolumbar spine finite element models to study spine biomechanics

Musculoskeletal disorders and low back pain (LBP) are common global afflictions, with a higher prevalence observed in females. However, the cause of many LBP cases continues to elude researchers. Current approaches seldom consider differences in male and female spines. Thus, this study aimed to compare the load distribution between male and female spines through finite element modeling. Two finite element models of the spine, one male and one female, were developed, inclusive of sex-specific geometry and material properties. The models consisted of the vertebrae, intervertebral discs (IVD), tendons, surrounding spinal muscles, and thoracolumbar fascia and were subjected to loading conditions simulating flexion and extension. Following extensive validation against published literature, intersegmental rotation, IVD stress, and vertebral body stress were evaluated. The female model demonstrated increased magnitudes for rotation and stresses when compared to the male model. Results suggest that the augmented stresses in the female model indicate an increased load distribution throughout the spine compared to the male model. These findings may corroborate the higher prevalence of LBP in females. This study highlights the importance of using patient- and sex-specific models for patient analyses and care.

Validation and Estimation of Obesity-Induced Intervertebral Disc Degeneration through Subject-Specific Finite Element Modelling of Functional Spinal Units

(1) Background: Intervertebral disc degeneration has been linked to obesity; its potential mechanical effects on the intervertebral disc remain unknown. This study aimed to develop and validate a patient-specific model of L3-L4 vertebrae and then use the model to estimate the impact of increasing body weight on disc degeneration. (2) Methods: A three-dimensional model of the functional spinal unit of L3-L4 vertebrae and its components were developed and validated. Validation was achieved by comparing the range of motions (RoM) and intradiscal pressures with the previous literature. Subsequently, the validated model was loaded according to the body mass index and estimated stress, deformation, and RoM to assess disc degeneration. (3) Results: During validation, L3-L4 RoM and intradiscal pressures: flexion 5.17° and 1.04 MPa, extension 1.54° and 0.22 MPa, lateral bending 3.36° and 0.54 MPa, axial rotation 1.14° and 0.52 MPa, respectively. When investigating the impact of weight on disc degeneration, escalating from normal weight to obesity reveals an increased RoM, by 3.44% during flexion, 22.7% during extension, 29.71% during lateral bending, and 33.2% during axial rotation, respectively. Also, stress and disc deformation elevated with increasing weight across all RoM. (4) Conclusions: The predicted mechanical responses of the developed model closely matched the validation dataset. The validated model predicts disc degeneration under increased weight and could lay the foundation for future recommendations aimed at identifying predictors of lower back pain due to disc degeneration.

Hybrid cortical bone trajectory and modified cortical bone trajectory techniques in transforaminal lumbar interbody fusion at L4-L5 segment: A finite element analysis

Background

The academia has increasingly acknowledged the superior biomechanical performance of the hybrid fixation technique in recent years. However, there is a lack of research on the hybrid fixation technique using BCS (Bilateral Cortical Screws) and BMCS (Bilateral Modified Cortical Screws). This study aims to investigate the biomechanical performance of the BCS and BMCS hybrid fixation technique in transforaminal lumbar interbody fusion (TLIF) at the L4-L5 segment in a complete lumbar-sacral finite element model.

Methods

Three cadaver specimens are used to construct three lumbar-sacral finite element models. The biomechanical properties of various fixation technologies (BCS-BCS, BMCS-BMCS, BMCS-BCS, and BCS-BMCS) are evaluated at the L4-5 segment with a TLIF procedure conducted, including the range of motion (ROM) of the L4-5 segment, as well as the stress experienced by the cage, screws, and rods. The testing is conducted under specific loading conditions, including a compressive load of 400 N and a torque of 7.5Nm, subjecting the model to simulate flexion, extension, lateral bending, and rotation.

Results

No significant variations are seen in the ROM at the L4-5 segment when comparing the four fixation procedures during flexion and extension. However, when it comes to lateral bending and rotation, the ROM is ordered in descending order as BCS-BCS, BCS-BMCS, BMCS-BMCS, and BMCS-BCS. The maximum stress experienced by the cage is observed to be highest within the BMCS-BCS technique during movements including flexion, extension, and lateral bending. Conversely, the BMCS-BMCS technique exhibits the highest cage stress levels during rotational movements. The stress applies to the screws and rods order the sequence of BCS-BCS, BCS-BMCS, BMCS-BCS, and BMCS-BMCS throughout all four working conditions.

Conclusion

The BMCS-BCS technique shows better biomechanical performance with less ROM and lower stress on the internal fixation system compared to other fixation techniques. BMCS-BMCS technology has similar mechanical performance to BMCS-BCS but has more contact area between screws and cortical bone, making it better for patients with severe osteoporosis.

Preliminary exploration of the biomechanical properties of three novel cervical porous fusion cages using a finite element study

BACKGROUND

Porous cages are considered a promising alternative to high-density cages because their interconnectivity favours bony ingrowth and appropriate stiffness tuning reduces stress shielding and the risk of cage subsidence.

METHODS

This study proposes three approaches that combine macroscopic topology optimization and micropore design to establish three new types of porous cages by integrating lattices (gyroid, Schwarz, body-centred cubic) with the optimized cage frame. Using these three porous cages along with traditional high-density cages, four ACDF surgical models were developed to compare the mechanical properties of facet articular cartilage, discs, cortical bone, and cages under specific loads.

RESULTS

The facet joints in the porous cage groups had lower contact forces than those in the high-density cage group. The intervertebral discs in all models experienced maximum stress at the C5/6 segment. The stress distribution on the cortical bone surface was more uniform in the porous cage groups, leading to increased average stress values. The gyroid, Schwarz, and BCC cage groups showed higher average stress on the C5 cortical bone. The average stress on the surface of porous cages was higher than that on the surface of high-density cages, with the greatest difference observed under the lateral bending condition. The BCC cage demonstrated favourable mechanical stability.

CONCLUSION

The new porous cervical cages satifies requirements of low rigidity and serve as a favourable biological scaffold for bone ingrowth. This study provides valuable insights for the development of next-generation orthopaedic medical devices.

Bibliometric and Visualization Analysis of Biomechanical Research on Lumbar Intervertebral Disc

Background

Biomechanical research on the lumbar intervertebral disc (IVD) provides valuable information for the diagnosis, treatment, and prevention of related diseases, and has received increasing attention. Using bibliometric methods and visualization techniques, this study investigates for the first time the research status and development trends in this field, with the aim of providing guidance and support for subsequent research.

Methods

The Science Citation Index Expanded (SCI-Expanded) within the Web of Science Core Collection (WoSCC) database was used as the data source to select literature published from 2003 to 2022 related to biomechanical research on lumbar IVD. VOSviewer 1.6.19 and CiteSpace 6.2.R2 visualization software, as well as the online analysis platform of literature metrology, were utilized to generate scientific knowledge maps for visual display and data analysis.

Results

The United States is the most productive country in this field, with the Ulm University making the largest contribution. Wilke HJ is both the most prolific author and one of the highly cited authors, while Adams MA is the most cited author. Spine, J Biomech, Eur Spine J, Spine J, and Clin Biomech are not only the journals with the highest number of publications, but also highly cited journals. The main research topics in this field include constructing and validating three-dimensional (3D) finite element model (FEM) of lumbar spine, measuring intradiscal pressure, exploring the biomechanical effects and related risk factors of lumbar disc degeneration, studying the mechanical responses to different torque load combinations, and classifying lumbar disc degeneration based on magnetic resonance images (MRI), which are also the hot research themes in recent years.

Conclusion

This study systematically reviews the knowledge system and development trends in the field of biomechanics of lumbar IVD, providing valuable references for further research.

Muscle-driven forward dynamic active hybrid model of the lumbosacral spine: combined FEM and multibody simulation

Most spine models belong to either the musculoskeletal multibody (MB) or finite element (FE) method. Recently, coupling of MB and FE models has increasingly been used to combine advantages of both methods. Active hybrid FE-MB models, still rarely used in spine research, avoid the interface and convergence problems associated with model coupling. They provide the inherent ability to account for the full interplay of passive and active mechanisms for spinal stability. In this paper, we developed and validated a novel muscle-driven forward dynamic active hybrid FE-MB model of the lumbosacral spine (LSS) in ArtiSynth to simultaneously calculate muscle activation patterns, vertebral movements, and internal mechanical loads. The model consisted of the rigid vertebrae L1-S1 interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, facet joints, and force actuators representing the muscles. Morphological muscle data were implemented via a semi-automated registration procedure. Four auxiliary bodies were utilized to describe non-linear muscle paths by wrapping and attaching the anterior abdominal muscles. This included an abdominal plate whose kinematics was optimized using motion capture data from upper body movements. Intra-abdominal pressure was calculated from the forces of the abdominal muscles compressing the abdominal cavity. For the muscle-driven approach, forward dynamics assisted data tracking was used to predict muscle activation patterns that generate spinal postures and balance the spine without prescribing accurate spinal kinematics. During calibration, the maximum specific muscle tension and spinal rhythms resulting from the model dynamics were evaluated. To validate the model, load cases were simulated from -10° extension to +30° flexion with weights up to 20 kg in both hands. The biomechanical model responses were compared with in vivo literature data of intradiscal pressures, intra-abdominal pressures, and muscle activities. The results demonstrated high agreement with this data and highlight the advantages of active hybrid modeling for the LSS. Overall, this new self-contained tool provides a robust and efficient estimation of LSS biomechanical responses under in vivo similar loads, for example, to improve pain treatment by spinal stabilization therapies.

Influence of Simulated State of Disc Degeneration and Axial Stiffness of Coupler in a Hybrid Performance Stabilisation System on the Biomechanics of a Spine Segment Model

Spinal fusion surgery leads to the restriction of mobility in the vertebral segments postoperatively, thereby causing stress to rise at the adjacent levels, resulting in early degeneration and a high risk of adjacent vertebral fractures. Thus, to address this issue, non-fusion surgery applies some pedicle screw-based dynamic stabilisation systems to provide stability and micromotion, thereby reducing stress in the fusion segments. Among these systems, the hybrid performance stabilisation system (HPSS) combines a rigid rod, transfer screw, and coupler design to offer a semi-rigid fixation method that preserves some mobility near the fusion site and reduces the adjacent segment compensatory effects. However, further research and confirmation are needed regarding the biomechanical effects of the dynamic coupler stiffness of the HPSS on the intrinsic degenerated adjacent segment. Therefore, this study utilised the finite element method to investigate the impact of the coupler stiffness of the HPSS on the mobility of the lumbar vertebral segments and the stress distribution in the intervertebral discs under flexion, extension, and lateral bending, as well as the clinical applicability of the HPSS on the discs with intrinsic moderate and severe degeneration at the adjacent level. The analytical results indicated that, regardless of the degree of disc degeneration, the use of a dynamic coupler stiffness of 57 N/mm in the HPSS may reduce the stress concentrations at the adjacent levels. However, for severely degenerated discs, the postoperative stress on the adjacent segments with the HPSS was still higher compared with that of the discs with moderate degeneration. We conclude that, when the discs had moderate degeneration, increasing the coupler stiffness led to a decrease in disc mobility. In the case of severe disc degeneration, the effect on disc mobility by coupler stiffness was less pronounced. Increasing the coupler stiffness ked to higher stress on intervertebral discs with moderate degeneration, while its effect on stress was less pronounced for discs with severe degeneration. It is recommended that patients with severe degeneration who undergo spinal dynamic stabilisation should remain mindful of the risk of accelerated adjacent segment degeneration.

Development and validation of lumbar spine finite element model

The functional biomechanics of the lumbar spine have been better understood by finite element method (FEM) simulations. However, there are still areas where the behavior of soft tissues can be better modeled or described in a different way. The purpose of this research is to develop and validate a lumbar spine section intended for biomechanical research. A FE model of the 50th percentile adult male (AM) Total Human Model for Safety (THUMS) v6.1 was used to implement the modifications. The main modifications were to apply orthotropic material properties and nonlinear stress-strain behavior for ligaments, hyperelastic material properties for annulus fibrosus and nucleus pulposus, and the specific content of collagenous fibers in the annulus fibrosus ground substance. Additionally, a separation of the nucleus pulposus from surrounding bones and tissues was implemented. The FE model was subjected to different loading modes, in which intervertebral rotations and disc pressures were calculated. Loading modes contained different forces and moments acting on the lumbar section: axial forces (compression and tension), shear forces, pure moments, and combined loading modes of axial forces and pure moments. The obtained ranges of motion from the modified numerical model agreed with experimental data for all loading modes. Moreover, intradiscal pressure validation for the modified model presented a good agreement with the data available from the literature. This study demonstrated the modifications of the THUMS v6.1 model and validated the obtained numerical results with existing literature in the sub-injurious range. By applying the proposed changes, it is possible to better model the behavior of the human lumbar section under various loads and moments.

Biomechanical evaluation of the hybrid pedicle screw-cortical bone trajectory technique in transforaminal lumbar interbody fusion to adjacent segment degeneration-finite element analysis

BACKGROUND

Transforaminal lumbar interbody fusion is an effective surgical treatment of intervertebral disk herniation. However, its clinical efficacy for adjacent segment disk degeneration (ASDD) after hybrid bilateral pedicle screw - bilateral cortical screw (pedicle screw at L4 and cortical bone trajectory screw at L5) and hybrid bilateral cortical screw - bilateral pedicle screw (bilateral cortical screw at L4 and bilateral pedicle screw at L5) remains undiscovered. Therefore, the aim of this study is to evaluate the effect of the hybrid bilateral pedicle screw - bilateral cortical screw and hybrid bilateral cortical screw - bilateral pedicle screw on the adjacent segment via a 3-dimensional (3D) finite element (FE) analysis.

METHODS

Four human cadaveric lumbar spine specimens were provided by the anatomy teaching and research department of Xinjiang Medical University. Four finite element models of L1-S1 lumbar spine segment were generated. For each of these, four lumbar transforaminal lumbar interbody fusion models at L4-L5 segment with the following instruments were created: hybrid bilateral pedicle screw - bilateral cortical screw, bilateral cortical screw - bilateral cortical screw (bilateral cortical screw at both L4 and L5 segments), bilateral pedicle screw - bilateral pedicle screw (bilateral pedicle screw at both L4 and L5 segments), and hybrid bilateral cortical screw - bilateral pedicle screw. A 400-N compressive load with 7.5 Nm moments was applied for the simulation of flexion, extension, lateral bending, and rotation. The range of motion of L3-L4 and L5-S1 segments and von Mises stress of the intervertebral disc at the adjacent segment were compared.

RESULTS

Hybrid bilateral pedicle screw - bilateral cortical screw has the lowest range of motion at L3-L4 segment in flexion, extension, and lateral bending, and the highest disc stress in all motions, while the range of motion at L5-S1 segment and disc stress was lower than bilateral pedicle screw - bilateral pedicle screw in flexion, extension, and lateral bending, and higher than bilateral cortical screw - bilateral cortical screw in all motions. The range of motion of hybrid bilateral cortical screw - bilateral pedicle screw at L3-L4 segment was lower than bilateral pedicle screw - bilateral pedicle screw and higher than bilateral cortical screw - bilateral cortical screw in flexion, extension, and lateral bending, and the range of motion at L5-S1 segment was higher than bilateral pedicle screw - bilateral pedicle screw in flexion, lateral bending, and axial rotation. The disc stress at L3-L4 segment was lowest and more dispersed in all motions, and the disc stress at L5-S1 segment was higher than bilateral pedicle screw - bilateral pedicle screw in lateral bending and axial rotation, but more dispersed.

CONCLUSION

Hybrid bilateral cortical screw - bilateral pedicle screw decreases the impact on adjacent segments after spinal fusion, reduces the iatrogenic injury to the paravertebral tissues, and provides throughout decompression of the lateral recess.

Hybrid pedicle screw and modified cortical bone trajectory technique in transforaminal lumbar interbody fusion at L4-L5 segment: finite element analysis

BACKGROUND

Investigate the biomechanical properties of the hybrid fixation technique with bilateral pedicle screw (BPS) and bilateral modified cortical bone trajectory screw (BMCS) in L4-L5 transforaminal lumbar interbody fusion (TLIF).

METHODS

 Three finite element (FE) models of the L1-S1 lumbar spine were established according to the three human cadaveric lumbar specimens. BPS-BMCS (BPS at L4 and BMCS at L5), BMCS-BPS (BMCS at L4 and BPS at L5), BPS-BPS (BPS at L4 and L5), and BMCS-BMCS (BMCS at L4 and L5) were implanted into the L4-L5 segment of each FE model. The range of motion (ROM) of the L4-L5 segment, von Mises stress of the fixation, intervertebral cage, and rod were compared under a 400-N compressive load with 7.5 Nm moments in flexion, extension, bending, and rotation.

RESULTS

 BPS-BMCS technique has the lowest ROM in extension and rotation, and BMCS-BMCS technique has the lowest ROM in flexion and lateral bending. The BMCS-BMCS technique showed maximal cage stress in flexion and lateral bending, and the BPS-BPS technique in extension and rotation. Compared to the BPS-BPS and BMCS-BMCS technique, BPS-BMCS technique presented a lower risk of screw breakage and BMCS-BPS technique presented a lower risk of rod breakage.

CONCLUSION

 The results of this study support that the use of the BPS-BMCS and BMCS-BPS techniques in TLIF surgery for offering the superior stability and a lower risk of cage subsidence and instrument-related complication.

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.

Biomechanical performance of the novel assembled uncovertebral joint fusion cage in single-level anterior cervical discectomy and fusion: A finite element analysis

Introduction: Anterior cervical discectomy and fusion (ACDF) is widely accepted as the gold standard surgical procedure for treating cervical radiculopathy and myelopathy. However, there is concern about the low fusion rate in the early period after ACDF surgery using the Zero-P fusion cage. We creatively designed an assembled uncoupled joint fusion device to improve the fusion rate and solve the implantation difficulties. This study aimed to assess the biomechanical performance of the assembled uncovertebral joint fusion cage in single-level ACDF and compare it with the Zero-P device.

Methods: A three-dimensional finite element (FE) of a healthy cervical spine (C2−C7) was constructed and validated. In the one-level surgery model, either an assembled uncovertebral joint fusion cage or a zero-profile device was implanted at the C5–C6 segment of the model. A pure moment of 1.0 Nm combined with a follower load of 75 N was imposed at C2 to determine flexion, extension, lateral bending, and axial rotation. The segmental range of motion (ROM), facet contact force (FCF), maximum intradiscal pressure (IDP), and screw−bone stress were determined and compared with those of the zero-profile device.

Results: The results showed that the ROMs of the fused levels in both models were nearly zero, while the motions of the unfused segments were unevenly increased. The FCF at adjacent segments in the assembled uncovertebral joint fusion cage group was less than that that of the Zero-P group. The IDP at the adjacent segments and screw–bone stress were slightly higher in the assembled uncovertebral joint fusion cage group than in those of the Zero-P group. Stress on the cage was mainly concentrated on both sides of the wings, reaching 13.4–20.4 Mpa in the assembled uncovertebral joint fusion cage group.

Conclusion: The assembled uncovertebral joint fusion cage provided strong immobilization, similar to the Zero-P device. When compared with the Zero-P group, the assembled uncovertebral joint fusion cage achieved similar resultant values regarding FCF, IDP, and screw–bone stress. Moreover, the assembled uncovertebral joint fusion cage effectively achieved early bone formation and fusion, probably due to proper stress distributions in the wings of both sides.

Biomechanical and clinical studies on lumbar spine fusion surgery: a review

Low back pain is associated with degenerative disc diseases of the spine. Surgical treatment includes fusion and non-fusion types. The gold standard is fusion surgery, wherein the affected vertebral segment is fused. The common complication of fusion surgery is adjacent segment degeneration (ASD). The ASD often leads to revision surgery, calling for a further fusion of adjacent segments. The existing designs of nonfusion type implants are associated with clinical problems such as subsidence, difficulty in implantation, and the requirement of revision surgeries. Various surgical approaches have been adopted by the surgeons to insert the spinal implants into the affected segment. Over the years, extensive biomechanical investigations have been reported on various surgical approaches and prostheses to predict the outcomes of lumbar spine implantations. Computer models have been proven to be very effective in identifying the best prosthesis and surgical procedure. The objective of the study was to review the literature on biomechanical studies for the treatment of lumbar spinal degenerative diseases. A critical review of the clinical and biomechanical studies on fusion spine surgeries was undertaken. The important modeling parameters, challenges, and limitations of the current studies were identified, showing the future research directions.

Biomechanical consequences of the intervertebral disc centre of rotation kinematics during lateral bending and axial rotation

The location of the instantaneous centre of rotation (ICR) of a lumbar unit has a considerable clinical importance as a spinal health estimator. Consequently, many studies have been conducted to measure or estimate the ICR during rotations in the three anatomical planes; however the results reported are widely scattered. Even if some inter-subjects variability is to be expected, such inconsistencies are likely explained by the differences in methods and experiments. Therefore, in this paper we seek to model three behaviours of the ICR during lateral bending and axial rotation based on results published in the literature. In order to assess the metabolic and mechanical sensibility to the assumption made on the ICR kinematics, we used a previously validated three dimensional non-linear poroelastic model of a porcine intervertebral disc to simulate physiological lateral and axial rotations. The impact of the geometry was also briefly investigated by considering a 11[Formula: see text] wedge angle. From our simulations, it appears that the hypothesis made on the ICR location does not significantly affect the critical nutrients concentrations but gives disparate predictions of the intradiscal pressure at the centre of the disc (variation up to 0.7 MPa) and of the displacement fields (variation up to 0.4 mm). On the contrary, the wedge angle does not influence the estimated intradiscal pressure but leads to minimal oxygen concentration decreased up to 33% and increased maximal lactate concentration up to 13%. While we can not settle on which definition of the ICR is more accurate, this work suggests that patient-specific modeling of the ICR is required and brings new insights that can be useful for the development of new tools or the design of surgical material such as total lumbar disc prostheses.

Biomechanical analysis of customized cage conforming to the endplate morphology in anterior cervical discectomy fusion: A finite element analysis

In anterior cervical discectomy and fusion (ACDF), an interbody fusion device is an essential implant. An unsuitable interbody fusion device can cause postoperative complications, including subsidence and nonunion. We designed a customized intervertebral fusion device to reduce postoperative complications and validated it by finite element analysis. Herein, we built a non-homogeneous model of the C3-7 cervical spine. Three implant models (customized cage, commercial cage, and bone graft cage) were constructed and placed in the C45 cervical segment after ACDF surgery. The simulated range of motion (ROM), stress at the cage-bone interface, and stress on the cage and implants were compared under different conditions. The commercial cage showed maximum stress peaks at 40.3 MPa and 43.2 MPa in the inferior endplate of C4 and superior endplate of C5 under rotational conditions, higher compared to 29.7 MPa and 26.4 MPa, respectively, in the customized cage. The ROM was not significantly different between the three cages placed after ACDF. The stresses on the commercial cage were higher compared to the other two cages under all conditions. The bone graft in the customized cage was subject to higher stress than the commercial cage under all conditions, particularly lateral bending, wherein the maximum stress was 5.5 MPa. These results showed that a customized cage that better conformed to the vertebral anatomy was promising for reducing the risk of stress shielding and the occurrence of subsidence.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

Yeditepe spine mesh: Finite element modeling and validation of a parametric CAD model of lumbar spine

Finite element analysis is a powerful tool that is often used to study the biomechanical response of the spine. The primary objective of this study was to illustrate the mechanical behavior of a previously proposed parametric CAD spine model in comparison with a segmented FSU model and the literature. In this study, two finite element models of the L4-L5 spinal level were developed from the same patient's CT scan data. The first was developed using well-known segmentation methods, whereas the second was developed from the new by using a novel parametric CAD model. Both models were subjected to the same loading and boundary conditions to perform flexion, extension, lateral bending and axial rotation motions. The segmented finite element model was observed to be in good agreement with the literature. The parametric finite element model results were also observed to be in good agreement with the segmented finite element model and with the literature except under extension.

Biomechanical Evaluation of the Cross-link Usage and Position in the Single and Multiple Segment Posterior Lumbar Interbody Fusion

OBJECTIVE

Previous studies have neither explored the usage of cross-links nor investigated the optimal position of the cross-links in posterior lumbar interbody fusion (PLIF). This study evaluates biomechanical properties of cross-links in terms of different fixation segments and optimal position in single- and multi-segment posterior lumbar interbody fusion.

METHODS

Two finite element (FE) models of instrumented lumbosacral spine with single-(L4/5) and multi-segment (L3-S1) PLIF surgery were simulated. On the basis of the two models, the benefits of the usage of cross-links were assessed and compared with the status of no application of cross-links. Moreover, the effects of position of cross-links on multi-segment PLIF surgery were studied in Upper, Middle, and Lower positions.

RESULTS

No significant difference was found in the range of motion (ROM), intersegmental rotational angle (IRA) of adjacent segments, and intradiscal pressure (IDP) regardless of the usage of cross-links in the single-segment PLIF surgery, while the cross-link increased the maximum von Mises stress in the fixation (MSF) under the axial rotation (53.65 MPa vs 41.42 MPa). In the multi-segment PLIF surgery, the usage of cross-links showed anti-rotational advantages indicated by ROM (Without Cross-link 2.35^o , Upper, 2.24^o ; Middle, 2.26^o ; Lower, 2.30^o ) and IRA (Without Cross-link 1.19^o , Upper, 1.08^o ; Middle, 1.09^o ; Lower, 1.13^o ). The greatest values of MSF were found in without cross-link case under the flexion, lateral bending, and axial rotation (37.48, 62.61, and 86.73 MPa). The application of cross-links at the Middle and Lower positions had lower values of MSF (48.79 and 69.62 MPa) under the lateral bending and axial rotation, respectively.

CONCLUSION

The application of cross-links was not beneficial for the single-segment PLIF, while it was found highly advantageous for the multi-segment PLIF. Moreover, the usage of cross-links at the Middle or Lower positions resulted in a better biomechanical stability.

Patient-Specific Finite Element Modeling of the Whole Lumbar Spine Using Clinical Routine Multi-Detector Computed Tomography (MDCT) Data-A Pilot Study

(1) Background: To study the feasibility of developing finite element (FE) models of the whole lumbar spine using clinical routine multi-detector computed tomography (MDCT) scans to predict failure load (FL) and range of motion (ROM) parameters. (2) Methods: MDCT scans of 12 subjects (6 healthy controls (HC), mean age ± standard deviation (SD): 62.16 ± 10.24 years, and 6 osteoporotic patients (OP), mean age ± SD: 65.83 ± 11.19 years) were included in the current study. Comprehensive FE models of the lumbar spine (5 vertebrae + 4 intervertebral discs (IVDs) + ligaments) were generated (L1−L5) and simulated. The coefficients of correlation (ρ) were calculated to investigate the relationship between FE-based FL and ROM parameters and bone mineral density (BMD) values of L1−L3 derived from MDCT (BMDQCT-L1-3). Finally, Mann−Whitney U tests were performed to analyze differences in FL and ROM parameters between HC and OP cohorts. (3) Results: Mean FE-based FL value of the HC cohort was significantly higher than that of the OP cohort (1471.50 ± 275.69 N (HC) vs. 763.33 ± 166.70 N (OP), p < 0.01). A strong correlation of 0.8 (p < 0.01) was observed between FE-based FL and BMDQCT-L1-L3 values. However, no significant differences were observed between ROM parameters of HC and OP cohorts (p = 0.69 for flexion; p = 0.69 for extension; p = 0.47 for lateral bending; p = 0.13 for twisting). In addition, no statistically significant correlations were observed between ROM parameters and BMDQCT- L1-3. (4) Conclusions: Clinical routine MDCT data can be used for patient-specific FE modeling of the whole lumbar spine. ROM parameters do not seem to be significantly altered between HC and OP. In contrast, FE-derived FL may help identify patients with increased osteoporotic fracture risk in the future.

Biomechanical Effects of a Novel Anatomic Titanium Mesh Cage for Single-Level Anterior Cervical Corpectomy and Fusion: A Finite Element Analysis

Background: The traditional titanium mesh cage (TTMC) has become common as a classical instrument for Anterior Cervical Corpectomy and Fusion (ACCF), but a series of complications such as cage subsidence, adjacent segment degeneration (ASD), and implant-related complications by using the TTMC have often been reported in the previous literature. The aim of this study was to assess whether a novel anatomic titanium mesh cage (NTMC) could improve the biomechanical condition after surgery.

Methods: The NTMC model consists of two spacers located on both sides of the TTMC which match the anatomic structure between the endplates by measuring patient preoperative cervical computed tomography (CT) data. The ranges of motion (ROMs) of the surgical segments and the stress peaks in the C6 superior endplates, titanium mesh cage (TMC), screw–bone interface, anterior titanium plate, and adjacent intervertebral disc were compared.

Results: Compared with the TTMC, the NTMC reduced the surgical segmental ROMs by 89.4% postoperatively. The C6 superior endplate stress peaks were higher in the TTMC (4.473–23.890 MPa), followed by the NTMC (1.923–5.035 MPa). The stress peaks on the TMC were higher in the TTMC (47.896–349.525 MPa), and the stress peaks on the TMC were lower in the NTMC (17.907–92.799 MPa). TTMC induced higher stress peaks in the screw–bone interface (40.0–153.2 MPa), followed by the NTMC (14.8–67.8 MPa). About the stress peaks on the anterior titanium plate, the stress of TTMC is from 16.499 to 58.432 MPa, and that of the NTMC is from 12.456 to 34.607 MPa. Moreover, the TTMC induced higher stress peaks in the C3/4 and C6/7 intervertebral disc (0.201–6.691 MPa and 0.248–4.735 MPa, respectively), followed by the NTMC (0.227–3.690 MPa and 0.174–3.521 MPa, respectively).

Conclusion: First, the application of the NTMC can effectively decrease the risks of TMC subsidence after surgery. Second, in the NTMC, the stresses at the anterior screw-plate, bone–screw, and TMC interface are much less than in the TTMC, which decreased the risks of instrument-related complications after surgery. Finally, increases in IDP at adjacent levels are associated with the internal stresses of adjacent discs which may lead to ASD; therefore, the NTMC can effectively decrease the risks of ASD.

Effect of the In Situ Screw Implantation Region and Angle on the Stability of Lateral Lumbar Interbody Fusion: A Finite Element Study

Objective

To investigate the effect of the in situ screw implantation region and angle on the stability of lateral lumbar interbody fusion (LLIF) from a biomechanical perspective.

Methods

A validated L2‐4 finite element (FE) model was modified for simulation. The L3‐4 fused segment undergoing LLIF surgery was modeled. The area between the superior and inferior edges and the anterior and posterior edges of the vertebral body (VB) is divided into four zones by three parallel lines in coronal and horizontal planes. In situ screw implantation methods with different angles based on the three parallel lines in coronal plane were applied in Models A, B, and C (A: parallel to inferior line; B: from inferior line to midline; C: from inferior line to superior line). In addition, four implantation methods with different regions based on the three parallel lines in horizontal plane were simulated as types 1–2, 1–3, 2–2, and 2–3 (1–2: from anterior line to midline; 1–3: from anterior line to posterior line; 2–2: parallel to midline; 2–3: from midline to posterior line). L3‐4 ROM, interbody cage stress, screw‐bone interface stress, and L4 superior endplate stress were tracked and calculated for comparisons among these models.

Results

The L3‐4 ROM of Models A, B, and C decreased with the extent ranging from 47.9% (flexion‐extension) to 62.4% (lateral bending) with no significant differences under any loading condition. Types 2–2 and 2–3 had 45% restriction, while types 1–2 and 1–3 had 51% restriction in ROM under flexion‐extension conditions. Under lateral bending, types 2–2 and 2–3 had 70.6% restriction, while types 1–2 and 1–3 had 61.2% restriction in ROM. Under axial rotation, types 2–2 and 2–3 had 65.2% restriction, while types 1–2 and 1–3 had 59.3% restriction in ROM. The stress of the cage in types 2–2 and 2–3 was approximately 20% lower than that in types 1–2 and 1–3 under all loading conditions in all models. The peak stresses at the screw‐bone interface in types 2–2 and 2–3 were much lower (approximately 35%) than those in types 1–2 and 1–3 under lateral bending, while no significant differences were observed under flexion‐extension and axial rotation. The peak stress on the L4 superior endplate was approximately 30 MPa and was not significantly different in all models under any loading condition.

Conclusions

Different regions of entry‐exit screws induced multiple screw trajectories and influenced the stability and mechanical responses. However, different implantation angles did not. Considering the difficulty of implantation, the ipsilateral‐contralateral trajectory in the lateral middle region of the VB can be optimal for in situ screw implantation in LLIF surgery.

Comparison of the susceptibility to implant failure in the lateral, posterior, and transforaminal lumbar interbody fusion: A finite element analysis

OBJECTIVE

There are several techniques for lumbar interbody fusion, and implant failure following lumbar interbody fusion can be troublesome. This study aimed to compare the stress in posterior implant and peri-screw vertebral bodies among lateral lumbar interbody fusion (LLIF), posterior lumbar interbody fusion (PLIF), and transforaminal lumbar interbody fusion (TLIF) and to select the technique that is least likely to cause implant failure.

METHODS

We created an intact L3-L5 model and simulated the LLIF, PLIF, and TLIF techniques at L4-L5 using finite element methods. All models at the lower portion of L5 were fixed and imposed a preload of 400 N and a moment of 7.5 Nm on the upper portion of L3 to simulate flexion, extension, lateral bending, and axial rotation. We investigated the peak stresses and stress concentration in the posterior implant and peri-screw vertebral bodies for the LLIF, PLIF, and TLIF techniques.

RESULTS

The extension, flexion, bending, and rotation peak stresses and stress concentration in the posterior implant, and the peri-screw vertebral bodies, were the lowest in LLIF, followed by PLIF and TLIF, respectively.

CONCLUSIONS

It was found that implant failure was least likely to occur in LLIF, followed by PLIF and TLIF, respectively. Hence, surgeons should be aware of these factors when selecting an appropriate surgical technique and be careful for implant failure during postoperative follow-up.

Development of a Finite Element Model of the Pediatric Thoracic and Lumbar Spine, Ribcage, and Pelvis with Orthotropic Region-Specific Vertebral Growth

Finite element (FE) modeling of the spine has increasingly been applied in orthopedic precision-medicine approaches. Previously published FE models of the pediatric spine growth have made simplifications in geometry of anatomical structures, material properties, and representation of vertebral growth. To address those limitations, a comprehensive FE model of a pediatric (10-year-old) osteo-ligamentous thoracic and lumbar spine (T1-L5 with intervertebral discs (IVDs) and ligaments), ribcage, and pelvis with age- and level-specific ligament properties and orthotropic region-specific vertebral growth was developed and validated. Range of motion (ROM) measures, namely lateral bending, flexion-extension, and axial rotation, of the current 10 YO FE model were generally within reported ranges of scaled in vitro adult ROM data. Changes in T1-L5 spine height, as well as kyphosis (T2-T12) and lordosis (L1-L5) angles in the current FE model for two years of growth (from ages 10 to 12 years) were within ranges reported from corresponding pediatric clinical data. The use of such comprehensive pediatric FE models can provide clinically relevant insights into normative and pathological biomechanical responses of the spine, and also contribute to the development and optimization of clinical interventions for spine deformities.

Optimization of a lumbar interspinous fixation device for the lumbar spine with degenerative disc disease

Interspinous spacer devices used in interspinous fixation surgery remove soft tissues in the lumbar spine, such as ligaments and muscles and may cause degenerative diseases in adjacent segments its stiffness is higher than that of the lumbar spine. Therefore, this study aimed to structurally and kinematically optimize a lumbar interspinous fixation device (LIFD) using a full lumbar finite element model that allows for minimally invasive surgery, after which the normal behavior of the lumbar spine is not affected. The proposed healthy and degenerative lumbar spine models reflect the physiological characteristics of the lumbar spine in the human body. The optimum number of spring turns and spring wire diameter in the LIFD were selected as 3 mm and 2 turns, respectively—from a dynamic range of motion (ROM) perspective rather than a structural maximum stress perspective—by applying a 7.5 N∙m extension moment and 500 N follower load to the LIFD-inserted lumbar spine model. As the spring wire diameter in the LIFD increased, the maximum stress generated in the LIFD increased, and the ROM decreased. Further, as the number of spring turns decreased, both the maximum stress and ROM of the LIFD increased. When the optimized LIFD was inserted into a degenerative lumbar spine model with a degenerative disc, the facet joint force of the L3-L4 lumbar segment was reduced by 56%–98% in extension, lateral bending, and axial rotation. These results suggest that the optimized device can strengthen the stability of the lumbar spine that has undergone interspinous fixation surgery and reduce the risk of degenerative diseases at the adjacent lumbar segments.

A biomechanical study of ligament tethers strengthening for the prevention of proximal junctional kyphosis after posterior long-segment spinal fusion

BACKGROUND

Proximal junctional kyphosis is a known complication of posterior long-segment thoracolumbar fusion. Here, the biomechanical effectiveness of ligament tethers strengthening and vertebral body augmentation, in proximal junctional kyphosis prevention was explored using the finite element analysis.

METHODS

Based on a validated model of T1-L5 with the pedicle screw system instrumented T8-L5, strengthening models with different strategies were created to assess the range of motion in proximal vertebrae, vertebrae stress, pedicle screw stress, and pressure on intervertebral discs during extension, flexion, lateral bending, and axial rotation motions. Strengthening strategies included two- and three-level posterior ligament tethers (TE-T7-T9 and TE-T6-T9), and tethers with T7 &T8 vertebral body augmentation (TECE-T7-T9 and TECE-T6-T9).

FINDINGS

Compared to the spinal fusion model, the ligament tethers strengthening significantly reduced the flexion-extension range of motion difference among the proximal vertebrae. During the flexion-extension motion, the T8 vertebra stresses in the TE-T7-T9, TE-T6-T9, TECE-T7-T9, and TECE-T6-T9 models were distinctively reduced, the values decreased by 26.8%, 28.3%, 28.8%, and 9.6%, respectively, during flexion, and by 21.9%, 35.2%, 23%, and 18.6%, respectively, during extension. In the strengthening models, the maximum stresses on the T7/T8 intervertebral disc in the TE-T6-T9 model were reduced by 13.8% during flexion and by 14.7% during extension.

INTERPRETATION

Based on our results, the ligament strengthening configuration of the three-level posterior tethers produced a more gradual transition in range of motion, vertebrae stresses, and intervertebral discs stress between the fused and non-fused segments, especially during flexion-extension, which may significantly decrease the proximal junctional kyphosis biomechanical risk.

Biomechanical and Clinical Study of Rod Curvature in Single-Segment Posterior Lumbar Interbody Fusion

Objective: Pedicle screw fixation is a common technique used in posterior lumbar interbody fusion (PLIF) surgery for lumbar disorders. During operation, rod contouring is often subjective and not satisfactory, but only few studies focused on the rod-contouring issue previously. The aim of the study was to explore the effect of the rod contouring on the single-segment PLIF by the finite element (FE) method and retrospective study.

Methods: A FE model of the lumbosacral vertebrae was first reconstructed, and subsequently single-segmental (L4/5) PLIF surgeries with four rod curvatures (RCs) were simulated. Herein, three RCs were designed by referring to centroid, Cobb, and posterior tangent methods applied in the lumbar lordosis measurement, and zero RC indicating straight rods was included as well. Clinical data of patients subjected to L4/5 segmental PLIF were also analyzed to verify the correlation between RCs and clinical outcome.

Results: No difference was observed among the four RC models in the range of motion (ROM), intersegmental rotation angle (IRA), and intradiscal pressure (IDP) under four actions. The posterior tangent model had less maximum stress in fixation (MSF) in flexion, extension, and axial rotation than the other RC models. Patients with favorable prognosis had larger RC and positive RC minus posterior tangent angle (RC-PTA) of fused segments with respect to those who had poor prognosis and received revision surgery.

Conclusion: All RC models had similar biomechanical behaviors under four actions. The posterior tangent-based RC model was superior in fixation stress distribution compared to centroid, Cobb, and straight models. The retrospective study demonstrated that moderate RC and positive RC-PTA were associated with better postoperative results.

Low-Modulus PMMA Has the Potential to Reduce Stresses on Endplates after Cement Discoplasty

Cement discoplasty has been developed to treat patients with advanced intervertebral disc degeneration. In discoplasty, poly(methylmethacrylate) (PMMA) bone cement is injected into the disc, leading to reduced pain and certain spinal alignment correction. Standard PMMA-cements have much higher elastic modulus than the surrounding vertebral bone, which may lead to a propensity for adjacent fractures. A PMMA-cement with lower modulus might be biomechanically beneficial. In this study, PMMA-cements with lower modulus were obtained using previously established methods. A commercial PMMA-cement (V-steady^®, G21 srl) was used as control, and as base cement. The low-modulus PMMA-cements were modified by 12 vol% (LA12), 16 vol% (LA16) and 20 vol% (LA20) linoleic acid (LA). After storage in 37 °C PBS from 24 h up to 8 weeks, specimens were tested in compression to obtain the material properties. A lower E-modulus was obtained with increasing amount of LA. However, with storage time, the E-modulus increased. Standard and low-modulus PMMA discoplasty were compared in a previously developed and validated computational lumbar spine model. All discoplasty models showed the same trend, namely a substantial reduction in range of motion (ROM), compared to the healthy model. The V-steady model had the largest ROM-reduction (77%), and the LA20 model had the smallest (45%). The average stress at the endplate was higher for all discoplasty models than for the healthy model, but the stresses were reduced for cements with higher amounts of LA. The study indicates that low-modulus PMMA is promising for discoplasty from a mechanical viewpoint. However, validation experiments are needed, and the clinical setting needs to be further considered.

FEBio finite element model of a pediatric cervical spine

OBJECTIVE

The underlying biomechanical differences between the pediatric and adult cervical spine are incompletely understood. Computational spine modeling can address that knowledge gap. Using a computational method known as finite element modeling, the authors describe the creation and evaluation of a complete pediatric cervical spine model.

METHODS

Using a thin-slice CT scan of the cervical spine from a 5-year-old boy, a 3D model was created for finite element analysis. The material properties and boundary and loading conditions were created and model analysis performed using open-source software. Because the precise material properties of the pediatric cervical spine are not known, a published parametric approach of scaling adult properties by 50%, 25%, and 10% was used. Each scaled finite element model (FEM) underwent two types of simulations for pediatric cadaver testing (axial tension and cardinal ranges of motion [ROMs]) to assess axial stiffness, ROM, and facet joint force (FJF). The authors evaluated the axial stiffness and flexion-extension ROM predicted by the model using previously published experimental measurements obtained from pediatric cadaveric tissues.

RESULTS

In the axial tension simulation, the model with 50% adult ligamentous and annulus material properties predicted an axial stiffness of 49 N/mm, which corresponded with previously published data from similarly aged cadavers (46.1 ± 9.6 N/mm). In the flexion-extension simulation, the same 50% model predicted an ROM that was within the range of the similarly aged cohort of cadavers. The subaxial FJFs predicted by the model in extension, lateral bending, and axial rotation were in the range of 1-4 N and, as expected, tended to increase as the ligament and disc material properties decreased.

CONCLUSIONS

A pediatric cervical spine FEM was created that accurately predicts axial tension and flexion-extension ROM when ligamentous and annulus material properties are reduced to 50% of published adult properties. This model shows promise for use in surgical simulation procedures and as a normal comparison for disease-specific FEMs.

MDCT-Based Finite Element Analysis for the Prediction of Functional Spine Unit Strength-An In Vitro Study

(1) Objective: This study aimed to analyze the effect of ligaments on the strength of functional spine unit (FSU) assessed by finite element (FE) analysis of anatomical models developed from multi-detector computed tomography (MDCT) data. (2) Methods: MDCT scans for cadaveric specimens were acquired from 16 donors (7 males, mean age of 84.29 ± 6.06 years and 9 females, mean age of 81.00 ± 11.52 years). Two sets of FSU models (three vertebrae + two disks), one with and another without (w/o) ligaments, were generated. The vertebrae were segmented semi-automatically, intervertebral disks (IVD) were generated manually, and ligaments were modeled based on the anatomical location. FE-predicted failure loads of FSU models (with and w/o ligaments) were compared with the experimental failure loads obtained from the uniaxial biomechanical test of specimens. (3) Results: The mean and standard deviation of the experimental failure load of FSU specimens was 3513 ± 1029 N, whereas of FE-based failure loads were 2942 ± 943 N and 2537 ± 929 N for FSU models with ligaments and without ligament attachments, respectively. A good correlation (ρ = 0.79, and ρ = 0.75) was observed between the experimental and FE-based failure loads for the FSU model with and with ligaments, respectively. (4) Conclusions: The FE-based FSU model can be used to determine bone strength, and the ligaments seem to have an effect on the model accuracy for the failure load calculation; further studies are needed to understand the contribution of ligaments.

Biomechanical modeling of spinal ligaments: finite element analysis of L4-L5 spinal segment

The complex mechanical structure of spine is usually simplified in finite element (FE) modes. In this study, different 3D models of L4-L5 spinal segment distinguished by their ligament modelling were developed (1D truss, 2D shell and 3D space truss elements). All models could be considered validated with respect to range of motion and intradiscal pressure, although their ligament stresses/forces were substantially different. The models with 2D shell and 3D space truss ligaments showed the stress distribution and identified the potential failure/injury locations in ligaments. The model with 3D space truss ligaments showed the stress/force direction (representing collagen fiber directions).

Investigation of Reaction Forces in the Thoracolumbar Fascia during Different Activities: A Mechanistic Numerical Study

Spinal instability remains a complex phenomenon to study while the cause of low back pain continues to challenge researchers. The role of fascia in biomechanics adds to the complexity of spine biomechanics but offers a new window from which to investigate our spines. Specifically, the thoracolumbar fascia may have an important role in spine biomechanics, and thus the purpose of this study was to access the mechanical influence of the thoracolumbar fascia on spine biomechanics during different simulated activities. A numerical finite element model of the lumbar spine inclusive of the intra-abdominal and intra-muscular regions as well as the thoracolumbar fascia was constructed and validated. Four different loading scenarios were simulated while deformation, stress, pressure, and reaction forces between the thoracolumbar fascia and spine were measured. Model validation was accomplished through comparison to in vivo and ex vivo published studies. Force transmission between the thoracolumbar fascia and the spine increased 40% comparing kyphotic and squatting lifting patterns. Further, the importance of reciprocating paraspinal and intra-abdominal pressures was demonstrated. It was also found that tension in the thoracolumbar fascia remains even in a simulated prone position. This numerical analysis allowed for an objective interpretation of the loads conveyed through the thoracolumbar fascia in different positional or lifting scenarios. Based on validation studies, it would appear to be a viable experimental platform from which insight can be derived. The loads in the thoracolumbar fascia vary considerably based on simulated tasks and are linked to the pressures in the paraspinal and intra-abdominal regions.

Sensitivity of Intervertebral Disc Finite Element Models to Internal Geometric and Non-geometric Parameters

Finite element models are useful for investigating internal intervertebral disc (IVD) behaviours without using disruptive experimental techniques. Simplified geometries are commonly used to reduce computational time or because internal geometries cannot be acquired from CT scans. This study aimed to (1) investigate the effect of altered geometries both at endplates and the nucleus-anulus boundary on model response, and (2) to investigate model sensitivity to material and geometric inputs, and different modelling approaches (graduated or consistent fibre bundle angles and glued or cohesive inter-lamellar contact). Six models were developed from 9.4 T MRIs of bovine IVDs. Models had two variations of endplate geometry (a simple curved profile from the centre of the disc to the periphery, and precise geometry segmented from MRIs), and three variations of NP-AF boundary (linear, curved, and segmented). Models were subjected to axial compressive loading (to 0.86 mm at a strain rate of 0.1/s) and the effect on stiffness and strain distributions, and the sensitivity to modelling approaches was investigated. The model with the most complex geometry (segmented endplates, curved NP-AF boundary) was 3.1 times stiffer than the model with the simplest geometry (curved endplates, linear NP-AF boundary), although this difference may be exaggerated since segmenting the endplates in the complex geometry models resulted in a shorter average disc height. Peak strains were close to the endplates at locations of high curvature in the segmented endplate models which were not captured in the curved endplate models. Differences were also seen in sensitivity to material properties, graduated fibre angles, cohesive rather than glued inter-lamellar contact, and NP:AF ratios. These results show that FE modellers must take care to ensure geometries are realistic so that load is distributed and passes through IVDs accurately.

Topological design and biomechanical evaluation for 3D printed multi-segment artificial vertebral implants

Customized spinal implants fabricated by additive manufacturing have been increasingly used clinically to restore the physiological functions. However, the mechanisms and methods about the design for the spinal implants are not clear, especially for the reconstruction of multi-segment vertebral. This study aims to develop a novel multi-objective optimization methodology based on various normal spinal activities, to design the artificial vertebral implant (AVI) with lightweight, high-strength and high-stability. The biomechanical performance for two types of AVI was analyzed and compared under different loading conditions by finite element method. These implants were manufactured via selective laser melting technology and evaluated via compressive testing. Results showed the maximum Mises stress of the optimized implant under various load cases were about 41.5% of that of the trussed implant, and below fatigue strength of 3D printed titanium materials. The optimized implant was about 2 times to trussed implant in term of the maximum compression load and compression stiffness to per unit mass, which indicated the optimized implant can meet the safety requirement. Finally, the optimized implant has been used in clinical practice and good short-term clinical outcomes were achieved. Therefore, the novel developed method provides a favorable guarantee for the design of 3D printed multi-segment artificial vertebral implants.

Experimental and Computational Comparison of Intervertebral Disc Bulge for Specimen-Specific Model Evaluation Based on Imaging

Finite element modelling of the spinal unit is a promising preclinical tool to assess the biomechanical outcome of emerging interventions. Currently, most models are calibrated and validated against range of motion and rarely directly against soft-tissue deformation. The aim of this contribution was to develop an in vitro methodology to measure disc bulge and assess the ability of different specimen-specific modelling approaches to predict disc bulge. Bovine bone-disc-bone sections (N = 6) were prepared with 40 glass markers on the intervertebral disc surface. These were initially magnetic resonance (MR)-imaged and then sequentially imaged using peripheral-qCT under axial compression of 1 mm increments. Specimen-specific finite-element models were developed from the CT data, using three different methods to represent the nucleus pulposus geometry with and without complementary use of the MR images. Both calibrated specimen-specific and averaged compressive material properties for the disc tissues were investigated. A successful methodology was developed to quantify the disc bulge in vitro, enabling observation of surface displacement on qCT. From the finite element model results, no clear advantage was found in using geometrical information from the MR images in terms of the models' ability to predict stiffness or disc bulge for bovine intervertebral disc.

Biomechanical Effects of a Cross Connector in Sacral Fractures - A Finite Element Analysis

Background: Spinopelvic fractures and approaches of operative stabilization have been a source of controversial discussion. Biomechanical data support the benefit of a spinopelvic stabilization and minimally invasive procedures help to reduce the dissatisfying complication rate. The role of a cross connector within spinopelvic devices remains inconclusive. We aimed to analyze the effect of a cross connector in a finite element model (FE model). Study Design: A FE model of the L1-L5 spine segment with pelvis and a spinopelvic stabilization was reconstructed from patient-specific CT images. The biomechanical relevance of a cross connector in a Denis zone I (AO: 61-B2) sacrum fracture was assessed in the FE model by applying bending and twisting forces with and without a cross connector. Biomechanical outcomes from the numerical model were investigated also considering uncertainties in material properties and levels of osseointegration. Results: The designed FE model showed comparable values in range-of-motion (ROM) and stresses with reference to the literature. The superiority of the spinopelvic stabilization (L5/Os ilium) ± cross connector compared to a non-operative procedure was confirmed in all analyzed loading conditions by reduced ROM and principal stresses in the disk L5/S1, vertebral body L5 and the fracture area. By considering the combination of all loading cases, the presence of a cross connector reduced the maximum stresses in the fracture area of around 10%. This difference has been statistically validated (p < 0.0001). Conclusion: The implementation of a spinopelvic stabilization (L5/Os ilium) in sacrum fractures sustained the fracture and led to enhanced biomechanical properties compared to a non-reductive procedure. While the additional cross connector did not alter the resulting ROM in L4/L5 or L5/sacrum, the reduction of the maximum stresses in the fracture area was significant.

Maintaining Bone Health in the Lumbar Spine: Routine Activities Alone Are Not Enough

Public health organisations typically recommend a minimum amount of moderate intensity activities such as walking or cycling for two and a half hours a week, combined with some more demanding physical activity on at least 2 days a week to maintain a healthy musculoskeletal condition. For populations at risk of bone loss in the lumbar spine, these guidelines are particularly relevant. However, an understanding of how these different activities are influential in maintaining vertebral bone health is lacking. A predictive structural finite element modelling approach using a strain-driven algorithm was developed to study mechanical stimulus and bone adaptation in the lumbar spine under various physiological loading conditions. These loading conditions were obtained with a previously developed full-body musculoskeletal model for a range of daily living activities representative of a healthy lifestyle. Activities of interest for the simulations include moderate intensity activities involving limited spine movements in all directions such as, walking, stair ascent and descent, sitting down and standing up, and more demanding activities with large spine movements during reaching and lifting tasks. For a combination of moderate and more demanding activities, the finite element model predicted a trabecular and cortical bone architecture representative of a healthy vertebra. When more demanding activities were removed from the simulations, areas at risk of bone degradation were observed at all lumbar levels in the anterior part of the vertebral body, the transverse processes and the spinous process. Moderate intensity activities alone were found to be insufficient in providing a mechanical stimulus to prevent bone degradation. More demanding physical activities are essential to maintain bone health in the lumbar spine.

Screws Fixation for Oblique Lateral Lumbar Interbody Fusion (OL-LIF): A Finite Element Study

Background

The combination of screw fixation and cage can provide stability in lumbar interbody fusion (LIF), which is an important technique to treat lumbar degeneration diseases. As the narrow surface cage is developed in oblique lateral lumbar interbody fusion (OL-LIF), screw fixation should be improved at the same time. We used the finite element (FE) method to investigate the biomechanics response by three different ways of screw fixation in OL-LIF.

Methods

Using a validated FE model, OL-LIF with 3 different screw fixations was simulated, including percutaneous transverterbral screw (PTVS) fixation, percutaneous cortical bone trajectory screw (PCBTS) fixation, and percutaneous transpedical screw (PPS) fixation. Range of motion (ROM), vertebral body displacement, cage displacement, cage stress, cortical bone stress, and screw stress were compared.

Results

ROM in FE models significantly decreased by 84-89% in flexion, 91-93% in extension, 78-89% in right and left lateral bending, and 73-82% in right and left axial rotation compared to the original model. The maximum displacement of the vertebral body and the cage in six motions except for the extension of model PTVS was the smallest among models. Meanwhile, the model PTVS had the higher stress of screw-rods system and also the lowest stress of cage. In all moments, the maximum stresses of the cages were lower than their yield stress.

Conclusions

Three screw fixations can highly restrict the surgical functional spinal unit (FSU). PTVS provided the better stability than the other two screw fixations. It may be a good choice for OL-LIF.

Biomechanical comparison of four types of instrumentation constructs for revision surgery in lumbar adjacent segment disease: A finite element study

BACKGROUND

Different constructs are applied in revision surgery (RS) for adjacent segment disease (ASD) aiming to further decompress and fixate the affected segment(s) in two ways: replacing or preserving the primary implants. This study aimed to compare the biomechanical properties of four constructs with different configurations.

METHODS

An T12-L5 finite element (FE) model was constructed and validated. Primary surgery was performed at L4-L5 and instrumented from L3 to L5. Thereafter, RS was undertook by decompressing L2-L3 and fixated with implant-replacing construct A, or implant-preserving construct B, C or D. Range of motion (ROM) and intervertebral disc pressure (IDP) were compared. Maximum von Mises stress on the rods between Construct A and B was evaluated.

RESULTS

An obvious reduction of ROM was observed when the FE model was instrumented with four constructs respectively. The overall changing characteristics of ROM were approximately identical among four constructs. The changing characteristic of IDP among four constructs was similar. The degree of IDP reduction of Construct B was comparable to Construct A, while that of Construct C was comparable to Construct D. Maximum von Mises stress on the rods between Construct A and B indicated that no stress concentration was recorded at the locking part of the connector rod.

CONCLUSIONS

The biomechanics of implant-preserving constructs were comparable to the traditional implant-replacing construct. The location of side-by-side connector could not affect the stability of Construct C and D. Construct B might be an optimal choice in RS for less dissection, less complication and more convenience in manipulation.

Biomechanical analysis of lumbar interbody fusion supplemented with various posterior stabilization systems

PURPOSE

Biomechanical comparison between rigid and non-rigid posterior stabilization systems following lumbar interbody fusion has been conducted in several studies. However, most of these previous studies mainly focused on investigating biomechanics of adjacent spinal segments or spine stability. The objective of the present study was to compare biomechanical responses of the fusion devices when using different posterior instrumentations.

METHODS

Finite-element model of the intact human lumbar spine (L1-sacrum) was modified to simulate implantation of the fusion cage at L4-L5 level supplemented with different posterior stabilization systems including (i) pedicle screw-based fixation using rigid connecting rods (titanium rods), (ii) pedicle screw-based fixation using flexible connecting rods (PEEK rods) and (iii) dynamic interspinous spacer (DIAM). Stress responses were compared among these various models under bending moments.

RESULTS

The highest and lowest stresses in endplate, fusion cage and bone graft were found at the fused L4-L5 level with DIAM and titanium rod stabilization systems, respectively. When using PEEK rod for the pedicle screw fixation, peak stress in the pedicle screw was lower but the ratio of peak stress in the rods to yield stress of the rod material was higher than using titanium rod.

CONCLUSIONS

Compared with conventional rigid posterior stabilization system, the use of non-rigid stabilization system (i.e., the PEEK rod system and DIAM system) following lumbar interbody fusion might increase the risks of cage subsidence and cage damage, but promote bony fusion due to higher stress in the bone graft. For the pedicle screw-based rod stabilization system, using PEEK rod might reduce the risk of screw breakage but increased breakage risk of the rod itself.

Biomechanical Effect of C5 /C6 Intervertebral Reconstructive Height on Adjacent Segments in Anterior Cervical Discectomy and Fusion - A Finite Element Analysis

OBJECTIVE

To investigate the biomechanical effect of different intervertebral reconstructive heights on adjacent segments following C₅ /C₆ anterior cervical discectomy and fusion (ACDF) through finite element analysis.

METHODS

A finite element model of intact C₄ -C₇ segments was developed and validated for the present study. Five additional C₄ -C₇ postoperative models were constructed with 100%, 125%, 150%, 175%, and 200% of the benchmark height of C₅ /C₆ on the basis of the intact model. The changes in intradiscal pressure (IDP) and range of motion (ROM) of adjacent segments before and after reconstruction of C₅ /C₆ were analyzed.

RESULTS

For the upper adjacent segment (C₄ /C₅ ), the IDPs under the different loading conditions all increased after reconstruction. The maximum IDPs were 0.387, 0.489, 0.491, and 0.472 MPa under flexion, extension, axial rotation, and lateral bending, respectively, observed at the reconstructive height of 200%. The minimum IDPs were observed at 150% reconstructive height under all loading conditions except extension, and were 57, 86 and 81% of the maximum IDPs under flexion, axial rotation, and lateral bending, respectively. The minimum IDP under extension occurred when the reconstructive height is 125% of the benchmark height. For the lower adjacent segment (C₆ /C₇ ), the IDPs of postoperative models under all loading conditions also increased compared to the preoperative model. The maximum IDPs after reconstruction under flexion, extension, axial rotation, and lateral bending were 0.402, 0.411, 0.461, and 0.497 MPa, respectively, when the height of the reconstruction was 200% of the benchmark. The minimum IDPs were observed after a reconstruction at 150% of the benchmark, and were 59%, 85%, 82%, and 81% of the maximum IDPs under flexion, extension, axial rotation, and lateral bending loading conditions.

CONCLUSIONS

The reconstructive height is an important factor affecting the IDP and the ROM of adjacent segments after ACDF. To delay the adjacent segment disease, an intervertebral reconstructive height of 150% is an appropriate height in C₅ /C₆ ACDF.

Calibration and validation of a novel hybrid model of the lumbosacral spine in ArtiSynth-The passive structures

In computational biomechanics, two separate types of models have been used predominantly to enhance the understanding of the mechanisms of action of the lumbosacral spine (LSS): Finite element (FE) and musculoskeletal multibody (MB) models. To combine advantages of both models, hybrid FE-MB models are an increasingly used alternative. The aim of this paper is to develop, calibrate, and validate a novel passive hybrid FE-MB open-access simulation model of a ligamentous LSS using ArtiSynth. Based on anatomical data from the Male Visible Human Project, the LSS model is constructed from the L1-S1 rigid vertebrae interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, and facet joints. A mesh convergence study, sensitivity analyses, and systematic calibration were conducted with the hybrid functional spinal unit (FSU) L4/5. The predicted mechanical responses of the FSU L4/5, the lumbar spine (L1-L5), and the LSS were validated against literature data from in vivo and in vitro measurements and in silico models. Spinal mechanical responses considered when loaded with pure moments and combined loading modes were total and intervertebral range of motions, instantaneous axes and centers of rotation, facet joint contact forces, intradiscal pressures, disc bulges, and stiffnesses. Undesirable correlations with the FE mesh were minimized, the number of crisscrossed collagen fiber rings was reduced to five, and the individual influences of specific anatomical structures were adjusted to in vitro range of motions. Including intervertebral motion couplings for axial rotation and nonlinear stiffening under increasing axial compression, the predicted kinematic and structural mechanics responses were consistent with the comparative data. The results demonstrate that the hybrid simulation model is robust and efficient in reproducing valid mechanical responses to provide a starting point for upcoming optimizations and extensions, such as with active skeletal muscles.

Numerical analysis of the mechanical behaviour of intact and implanted lumbar functional spinal units: Effects of loading and boundary conditions

The objective of this study was to develop an improved finite element (FE) model of a lumbar functional spinal unit (FSU) and to subsequently analyse the deviations in load transfer owing to implantation. The effects of loading and boundary conditions on load transfer in intact and implanted FSUs and its relationship with the potential risk of vertebral fracture were investigated. The FE models of L1-L5 and L3-L4 FSUs, intact and implanted, were developed using patient-specific CT-scan dataset and segmentation of cortical and cancellous bone regions. The effect of submodelling technique, as compared to artificial boundary conditions, on the elastic behaviour of lumbar spine was examined. Applied forces and moments, corresponding to physiologic movements, were used as loading conditions. Results indicated that the loading and boundary conditions considerably affect stress-strain distributions within a FSU. This study, based on an improved FE model of a vertebra, highlights the importance of using the submodelling technique to adequately evaluate the mechanical behaviour of a FSU. In the intact FSU, strains of 200-400 µε were observed in the cancellous bone of vertebral body and pedicles. High equivalent stresses of 10-25 MPa and 1-5 MPa were generated around the pars interarticularis for cortical and cancellous regions, respectively. Implantation caused reductions of 85%-92% in the range of motion for all movements. Insertion of the intervertebral cage resulted in major deviations in load transfer across a FSU for all movements. The cancellous bone around cage experienced pronounced increase in stresses of 10-15 MPa, which indicated potential risk of failure initiation in the vertebra.

Biomechanical effect of different plate-to-disc distance on surgical and adjacent segment in anterior cervical discectomy and fusion - a finite element analysis

Background

The plate-to-disc distance (PDD) is an important factor affecting the degeneration of adjacent segments after anterior cervical discectomy and fusion (ACDF). However, the most suitable PDD is controversial. This study examined the adjacent intervertebral disc stress, bone graft stress, titanium plate stress and screw stress to evaluate the biomechanical effect of different PDD on surgical segment and adjacent segment following C5/C6 ACDF.

Methods

We constructed 10 preoperative finite element models of intact C4–C7 segments and validated them in the present study. We simulated ACDF surgery based on the 10 intact models in software. We designed three different distance of plate-to-disc titanium plates: long PDD (10 mm), medium PDD (5 mm) and short PDD (0 mm). The changes in C4/C5 and C6/C7 intervertebral disc stress, bone graft stress, titanium plate stress and screw stress were analyzed.

Results

The von Mises stress of C4/C5 and C6/C7 intervertebral discs had no significant differences (P > 0.05) in three different PDD groups. Titanium plate stress increased as the PDD decreased. The bone graft stress and screws stress decreased as the PDD decreased. The maximum stress of each part occurred was mostly in the conditions of rotation and lateral bending.

Conclusions

The PDD has no effect on adjacent intervertebral disc stress, but it is an important factor that affecting the bone graft stress, titanium plate stress and screws stress after ACDF. Shorter PDD plate can provide better stability to reduce stress on screws and bone graft, which may be helpful to prevent cage subsidence, pseudarthrosis and instrument failure. This can serve as a reference for clinical choice of plate.

MDCT-Based Finite Element Analyses: Are Measurements at the Lumbar Spine Associated with the Biomechanical Strength of Functional Spinal Units of Incidental Osteoporotic Fractures along the Thoracolumbar Spine?

Assessment of osteoporosis-associated fracture risk during clinical routine is based on the evaluation of clinical risk factors and T-scores, as derived from measurements of areal bone mineral density (aBMD). However, these parameters are limited in their ability to identify patients at high fracture risk. Finite element models (FEMs) have shown to improve bone strength prediction beyond aBMD. This study aims to investigate whether FEM measurements at the lumbar spine can predict the biomechanical strength of functional spinal units (FSUs) with incidental osteoporotic vertebral fractures (VFs) along the thoracolumbar spine. Multi-detector computed tomography (MDCT) data of 11 patients (5 females and 6 males, median age: 67 years) who underwent MDCT twice (median interval between baseline and follow-up MDCT: 18 months) and sustained an incidental osteoporotic VF between baseline and follow-up scanning were used. Based on baseline MDCT data, two FSUs consisting of vertebral bodies and intervertebral discs (IVDs) were modeled: one standardly capturing L1-IVD-L2-IVD-L3 (FSU_L1-L3) and one modeling the incidentally fractured vertebral body at the center of the FSU (FSU_F). Furthermore, volumetric BMD (vBMD) derived from MDCT, FEM-based displacement, and FEM-based load of the single vertebrae L1 to L3 were determined. Statistically significant correlations (adjusted for a BMD ratio of fracture/L1-L3 segments) were revealed between the FSU_F and mean load of L1 to L3 (r = 0.814, p = 0.004) and the mean vBMD of L1 to L3 (r = 0.745, p = 0.013), whereas there was no statistically significant association between the FSU_F and FSU_L1-L3 or between FSU_F and the mean displacement of L1 to L3 (p > 0.05). In conclusion, FEM measurements of single vertebrae at the lumbar spine may be able to predict the biomechanical strength of incidentally fractured vertebral segments along the thoracolumbar spine, while FSUs seem to predict only segment-specific fracture risk.

Biomechanical Evaluation of Oblique Lumbar Interbody Fusion with Various Fixation Options: A Finite Element Analysis

Objective

The aim of the present study was to clarify the biomechanical properties of oblique lumbar interbody fusion (OLIF) using different fixation methods in normal and osteoporosis spines.

Methods

Normal and osteoporosis intact finite element models of L₁–S₁ were established based on CT images of a healthy male volunteer. Group A was the normal models and group B was the osteoporosis model. Each group included four subgroups: (i) intact; (ii) stand‐alone cage (Cage); (iii) cage with lateral plate and two lateral screws (LP); and (iv) cage with bilateral pedicle screws and rods (BPSR). The L₃–L₄ level was defined as the surgical segment. After validating the normal intact model, compressive load of 400 N and torsional moment of 10 Nm were applied to the superior surface of L₂ to simulate flexion, extension, left bending, right bending, left rotation, and right rotation motions. Surgical segmental range of motion (ROM), cage stress, endplate stress, supplemental fixation stress, and stress distribution were analyzed in each group.

Results

Cage provided the minimal reduction of ROM among all motions (normal, 82.30%–98.81%; osteoporosis, 92.04%–97.29% of intact model). BPSR demonstrated the maximum reduction of ROM (normal, 43.94%–61.13%; osteoporosis, 45.61%–62.27% of intact model). The ROM of LP was between that of Cage and BPSR (normal, 63.25%–79.72%; osteoporosis, 70%–87.15% of intact model). Cage had the minimal cage stress and endplate stress. With the help of LP and BPSR fixation, cage stress and endplate stress were significantly reduced in all motions, both in normal and osteoporosis finite element models. However, BPSR had more advantages. For cage stress, BPSR was at least 75.73% less than that of Cage in the normal model, and it was at least 80.10% less than that of Cage in the osteoporosis model. For endplate stress, BPSR was at least 75.98% less than that of Cage in the normal model, and it was at least 78.06% less than that of Cage in the osteoporosis model. For supplemental fixation stress, BPSR and LP were much less than the yield strength in all motions in the two groups. In addition, the comparison between the two groups showed that the ROM, cage stress, endplate stress, and supplemental fixation stress in the normal model were less than in the osteoporosis model when using the same fixation option of OLIF.

Conclusion

Oblique lumbar interbody fusion with BPSR provided the best biomechanical stability both in normal and osteoporosis spines. The biomechanical properties of the normal spine were better than those of the osteoporosis spine when using the same fixation option of OLIF.