Applying a follower load delivers realistic results for simulating standing

Anatomically and mechanically conforming patient-specific spinal fusion cages designed by full-scale topology optimization

Cage subsidence after instrumented lumbar spinal fusion surgery remains a significant cause of treatment failure, specifically for posterior or transforaminal lumbar interbody fusion. Recent advancements in computational techniques and additive manufacturing, have enabled the development of patient-specific implants and implant optimization to specific functional targets. This study aimed to introduce a novel full-scale topology optimization formulation that takes the structural response of the adjacent bone structures into account in the optimization process. The formulation includes maximum and minimum principal strain constraints that lower strain concentrations in the adjacent vertebrae. This optimization approach resulted in anatomically and mechanically conforming spinal fusion cages. Subsidence risk was quantified in a commercial finite element solver for off-the-shelf, anatomically conforming and the optimized cages, in two representative patients. We demonstrated that the anatomically and mechanically conforming cages reduced subsidence risk by 91% compared to an off-the-shelf implant with the same footprint for a patient with normal bone quality and 54% for a patient with osteopenia. Prototypes of the optimized cage were additively manufactured and mechanically tested to evaluate the manufacturability and integrity of the design and to validate the finite element model.

Primary Stability of Kyphoplasty in Incomplete Vertebral Body Burst Fractures in Osteoporosis: A Biomechanical Investigation

Background: The objective of our study was to biomechanically evaluate the use of kyphoplasty to stabilize post-traumatic segmental instability in incomplete burst fractures of the vertebrae. Methods: The study was performed on 14 osteoporotic spine postmortem samples (Th11-L3). First, acquisition of the native multisegmental kinematics in our robot-based spine tester with three-dimensional motion analysis was set as a baseline for each sample. Then, an incomplete burst fracture was generated in the vertebral body L1 with renewed kinematic testing. After subsequent kyphoplasty was performed on the fractured vertebral body, primary stability was examined again. Results: Initially, a significant increase in the range of motion after incomplete burst fracture generation in all three directions of motion (extension-flexion, lateral tilt, axial rotation) was detected as proof of post-traumatic instability. There were no significant changes to the native state in the adjacent segments. Radiologically, a significant loss of height in the fractured vertebral body was also shown. Traumatic instability was significantly reduced by kyphoplasty. However, native kinematics were not restored. Conclusions: Although post-traumatic segmental instability was significantly reduced by kyphoplasty in our in vitro model, native kinematics could not be reconstructed, and significant instability remained.

Effect of compressive and tensile forces on glucose concentration and cell viability within the intervertebral disc: A finite element study

Understanding the role of mechanical force on tissue nutrient transport is essential, as sustained force may affect nutrient levels within the disc and initiate disc degeneration. This study aims to evaluate the time-dependent effects of different compressive force amplitudes as well as tensile force on glucose concentration and cell viability within the disc. Based on the mechano-electrochemical mixture theory, a multiphasic finite element model of the lumbar intervertebral disc was developed. The minimum glucose concentration and minimum cell density in both normal and degenerated discs were predicted for different compressive force amplitudes, tensile force, and corresponding creep time. Under high compressive force, the minimum glucose concentration exhibited an increasing and then decreasing trend with creep time in the normal disc, whereas that of the degenerated disc increased, then decreased, and finally increased again. At steady state, a higher compressive force was accompanied by a lower glucose concentration distribution. In the degenerated disc, the minimum cell density was negatively correlated with creep time, with a greater range of affected tissue under a higher compressive force. For tensile force, the minimum glucose concentration of the degenerated disc raised over time. This study highlighted the importance of creep time, force magnitude, and force type in affecting nutrient concentration and cell viability. Sustained weight-bearing activities could deteriorate the nutrient environment of the degenerated disc, while tensile force might have a nonnegligible role in effectively improving nutrient levels within the degenerated disc.

The effect of polymethylmethacrylate augmentation on the primary stability of stand-alone implant construct versus posterior stabilization in oblique lumbar interbody fusion with osteoporotic bone quality- a finite element study

BACKGROUND CONTEXT

Oblique lumbar interbody fusion (OLIF) can provide an ideal minimally invasive solution for achieving spinal fusion in an older, more frail population where decreased bone quality can be a limiting factor. Stabilization can be achieved with bilateral pedicle screws (BPS), which require additional incisions and longer operative time. Alternatively, a novel self-anchoring stand-alone lateral plate system (SSA) can be used, where no additional incisions are required. Based on the relevant literature, BPS constructs provide greater primary biomechanical stability compared to lateral plate constructs, including SSA. This difference is further increased by osteoporosis. Screw augmentation in spinal fusion surgeries is commonly used; however, in the case of OLIF, it is a fairly new concept, lacking a consensus-based guideline.

PURPOSE

This comparative finite element (FE) study aimed to investigate the effect of PMMA screw augmentation on the primary stability of a stand-alone implant construct versus posterior stabilization in OLIF with osteoporotic bone quality.

STUDY DESIGN

The biomechanical effect of screw augmentation was studied inside an in-silico environment using computer-aided FE analysis.

METHODS

A previously validated and published L2-L4 FE model with normal and osteoporotic bone material properties was used. Geometries based on the OLIF implants (BPS, SSA) were created and placed inside the L3-L4 motion segment with increasing volumes (1-6 cm3) of PMMA augmentation. A follower load of 400 N and 10 Nm bending moment (in the three anatomical planes) were applied to the surgical FE models with different bone material properties. The operated L3-L4 segmental range of motion (ROM), the inserted cage's maximal caudal displacements, and L4 cranial bony endplate principal stress values were measured.

RESULTS

The nonaugmented values for the BPS construct were generally lower compared to SSA, and the difference was increased by osteoporosis. In osteoporotic bone, PMMA augmentation gradually decreased the investigated parameters and the difference between the two constructs as well. Between 3 cm3 and 4 cm3 of injected PMMA volume per screw, the difference between augmented SSA and standard BPS became comparable.

CONCLUSIONS

Based on this study, augmentation can enhance the primary stability of the constructs and decrease the difference between them. Considering leakage as a possible complication, between 3 cm3 and 4 cm3 of injected PMMA per screw can be an adequate amount for SSA augmentation. However, further in silico, and possibly in vitro and clinical testing is required to thoroughly understand the investigated biomechanical aspects.

CLINICAL SIGNIFICANCE

This study sheds light on the possible biomechanical advantage offered by augmented OLIF implants and provides a theoretical augmentation amount for the SSA construct. Based on the findings, the concept of an SSA device with PMMA augmentation capability is desirable.

Pediatric Solid-State 3D Models of Lumbar Vertebrae and Spine

Introduction While various 3D vertebral models have been utilized in numerous studies, there is a notable gap in the representation of pediatric lumbar vertebrae and spine. This study aimed to describe the changing shapes of lumbar vertebrae and spine with age and to develop precise 3D models. Materials and methods Solid-state 3D models of pediatric lumbar vertebrae and spine were created using SOLIDWORKS® Simulation software for five age groups: newborns, infants (ages 0-1), toddlers (ages 1-3), middle childhood (ages 4-7), and preadolescents (ages 8-12). Models were composed of components with varying biomechanical characteristics. Results Created 3D models replicate variations in the dimensions and configurations of vertebrae, taking into account osteometric analyses conducted on actual vertebral specimens. These models also include elements made of cartilage, representing various phases of vertebral growth during ontogeny. Additionally, through 3D parametric design, we developed comprehensive lumbar spine models, incorporating both the vertebrae and intervertebral disks. Conclusion Created pediatric solid-state vertebral 3D models can be utilized in developing virtual or augmented reality applications and for medical research. Users can interact with models, allowing virtual exploration and manipulation, enhancing learning experiences and facilitating a better understanding of spatial relationships. These solid-state 3D models allow finite element analysis and can be used for further research to calculate internal relative deformations and stress distribution under different conditions.

Finite element method-based study for spinal vibration characteristics of the scoliosis and kyphosis lumbar spine to whole-body vibration under a compressive follower preload

PURPOSE

To analyze the dynamic response of the lumbosacral vertebrae structure of a scoliosis spine and a kyphosis spine under whole-body vibration.

METHODS

Typical Lenke4 (kyphosis) and Lenke3 (scoliosis) spinal columns were used as research objects. A finite element model of the lumbosacral vertebrae segment was established and validated based on CT scanning images. Modal, harmonic response, and transient dynamic analyses were performed on the lumbar-sacral scoliosis model using the finite element software abaqus.

RESULTS

The first four resonance frequencies of kyphosis spine extracted from modal analysis were 0.86, 1.45, 8.51, and 55.71 Hz. The first four resonance frequencies of scoliosis spine extracted from modal analysis were 0.76, 1.45, 10.51, and 63.82 Hz. The scoliosis spine had the maximum resonance amplitude in the transverse direction, while the kyphosis spine had the maximum resonance amplitude in the anteroposterior direction. The dynamic response in transient analysis exhibited periodic response over time at all levels.

CONCLUSION

The scoliosis and kyphosis deformity of the spine significantly complicates the vibration response in the scoliosis and kyphosis areas at the top of the spine. Scoliosis and kyphosis patients are more likely to experience vibrational spinal diseases than healthy people. Besides, applying vertical cyclic loads on a malformed spine may cause further rotation of scoliosis and kyphosis deformities.

A finite element method study of the effect of vibration on the dynamic biomechanical response of the lumbar spine

BACKGROUND

Studies focusing on lumbar spine biomechanics are very limited, and the mechanism of the effect of vibration on lumbar spine biodynamics is unclear. To provide guidance and reference for lumbar spine biodynamics research and vibration safety assessment, this study aims to investigate the effects of different vibrations on lumbar spine biodynamics.

METHODS

A validated finite element model of the lumbosacral spine was utilized. The model incorporated a 40 kg mass on the upper side and a 400 N follower preload. As a comparison, another model without a coupled mass was also employed. A sinusoidal acceleration with an amplitude of 1 m/s2 and a frequency of 5 Hz was applied to the upper and lower sides of the model respectively.

FINDINGS

When the coupled mass point is not introduced: in the case of upper-side excitation, the lumbar spine shows a significantly larger response in the x-direction than in the z-direction, while in the case of lower-side excitation, the lumbar spine experiences rigid body displacement in the z-direction without any movement, deformation, rotation, or stress changes in the x-direction. When the coupled mass point is introduced: both upper and lower-side excitations result in significant differences in z-directional displacement, with relatively small differences in vertebral rotation angle, disc deformation, and stress. Under upper excitation, low-frequency oscillations occur in the x-direction. In both types of excitations, the anterior-posterior deformation of the L2-L3 and L4-L5 intervertebral discs is greater than the vertical deformation. The peak (maximum) disc stress exceeds the average stress and stress amplitude across the entire disc. Regardless of the excitation type, the stress distribution within the disc at the moment of peak displacement remains nearly identical, with the maximum stress consistently localized on the anterior side of the L4-L5 disc.

INTERPRETATION

Accurately simulating lumbar spine biodynamics requires the inclusion of the upper body mass in the lumbosacral spine model. The physiological curvature of the lumbar spine could escalate the risk of lumbar spine vibration injuries. It is more instructive to apply local high stress in the disc as a lumbar spine vibration safety evaluation parameter.

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.

Biomechanical response of lumbar intervertebral disc in daily sitting postures: a poroelastic finite element analysis

This study aims to establish and validate a poroelastic L4-L5 finite element model to evaluate the effect of different sitting postures and their durations on the mechanical responses of the disc. During the sustained loading conditions, the height loss, fluid loss and von-Mises stress gradually increased, but the intradiscal pressure decreased. The varying rates of aforementioned parameters were more significant at the initial loading stage and less so at the end. The predicted values in the flexed sitting posture were significantly greater than other postures. The extended sitting posture caused an obvious von-Mises stress concentration in the posterior region of the inter-lamellar matrix. From the biomechanical perspective, prolonged sitting may pose a high risk of lumbar disc degeneration, and therefore adjusting the posture properly in the early stage of sitting time may be useful to mitigate that. Additionally, upright sitting is a safer posture, while flexed sitting posture is more harmful.

Clinic choice of long or short segment pedicle screw-rod fixation in the treatment of thoracolumbar burst fracture: From scan data to numerical study

Based on computerized tomography scanning images of human lumbar vertebrae, finite element (FE) analysis is performed to predict the stress of pedicle screws, rods, and fractured vertebra as well as the displacement of fractured vertebra after internal fixation treatment of thoracolumbar burst fracture. A three-dimensional FE model of L1-L5 lumbar vertebrae with L3 burst fracture has been established and four fixation methods, namely, short segment cross- and trans-injured vertebrae, long segment cross- and trans-injured vertebrae fixations, have been adopted to perform posterior pedicle fixation. The stress distributions of the screws, rods, and fractured vertebra and the total deformation of the fractured vertebra are investigated under six different physiological motions. From the view of the stress on the screw-rod system and the deformation of the fractured vertebral body, the long segment cross-injured vertebra fixation has the best mechanical performance, followed by the long segment trans-injured vertebra fixation, and then the short segment fixation trans-injured vertebra. The short segment fixation cross-injured vertebra performs the worst. Among the six motions, the forward flexion movement has the greatest impact on the screw-rod system and the fractured vertebra. However, the rotation motion greatly affects the stress of the screw in the long segment fixation. This indicates that the longer the fixed segment is, the more susceptible it is to human rotation. Thus, for patients with severe fracture, the long segment cross-injured vertebra is preferred. On the contrary, the short segment trans-injured vertebra fixation is optimal.

Biomechanical evaluation of the thoracolumbar spine comparing healthy and irregular thoracic and lumbar curvatures

BACKGROUND

The geometric alignment of the spine plays an integral role in stability, biomechanical loading, and consequently, pain, and a range of healthy sagittal curvatures has been identified. Spinal biomechanics when sagittal curvature is outside the optimal range remains a debate and may provide insight into the load distribution throughout the spinal column.

METHOD

A thoracolumbar spine model (Healthy) was developed. Thoracic and lumbar curvatures were adjusted by 50% to create models with varying sagittal profiles: hypolordotic (HypoL), hyperlordotic (HyperL), hypokyphotic (HypoK), and hyperkyphotic (HyperK). In addition, lumbar spine models were constructed for the former three profiles. The models were subjected to loading conditions simulating flexion and extension. Following validation, intervertebral disc stresses, vertebral body stresses, disc heights, and intersegmental rotations were compared across all models.

RESULTS

Overall trends demonstrated that HyperL and HyperK models had a noticeable reduction in disc height and greater vertebral body stresses compared to the Healthy model. In comparison, the HypoL and HypoK models displayed opposite trends. Considering the lumbar models, the HypoL model had reduced disc stresses and flexibility, while the contrary was observed in the HyperL model. Results indicate that the models with excessive curvature may be subjected to greater stress magnitudes, while the straighter spine models may reduce these stresses.

CONCLUSIONS

Finite element modeling of spine biomechanics demonstrated that variations in sagittal profiles influence the load distribution and range of motion of the spine. Considering patient-specific sagittal profiles in finite element modeling may provide valuable insight for biomechanical analyses and targeted treatments.

Effect of different fixation methods on biomechanical property of cervical vertebral body replacement and fusion

BACKGROUND

The main purpose of this study was to examine the effect of different fixation methods (anterior fixation, self-stabilizing fixation and anterior-posterior fixation) on biomechanical property of vertebral body replacement and fusion.

METHODS

Three finite element models of cervical vertebral body replacement and fusion were established. The implanted models included artificial vertebral body and fixation system, and the loads imposed on the models included 75 N compression load and 1 Nm moment load.

FINDINGS

For anterior-posterior fixation, the cervical load was mainly transmitted by the posterior pedicle screw and rod (more than 50%), and the stress shielding problem was the most significant than the self-stabilizing and anterior fixation. Self-stabilizing fixation was more helpful to the fusion of implant and vertebrae, but the higher risk of vertebral body collapse was worthy of attention if the cervical spine with osteoporosis. The stress of bone was mainly concentrated around the screw hole. The maximum stress (20.03 MPa) was lower than the yield stress of cortical bone and the possibility of fracture around the fixation device of cervical spine was low. The anterior fixation could meet the requirement of vertebral body replacement and fusion, and the addition of posterior pedicle screws and rods might obtain better treatment in cases of severe spine injury or osteoporosis.

INTERPRETATION

The findings of this study may provide guidance on clinical treatments for choosing more appropriate fixation methods for different patients.

Biomechanical effect of age-related structural changes on cervical intervertebral disc: A finite element study

Previous literature has investigated the biomechanical response of healthy and degenerative discs, but the biomechanical response of suboptimal healthy intervertebral discs received less attention. The purpose was to compare the biomechanical responses and risk of herniation of young healthy, suboptimal healthy, and degenerative intervertebral discs. A cervical spine model was established and validated using the finite element method. Suboptimal healthy, mildly, moderately, and severely degenerative disc models were developed. Disc height deformation, range of motion, intradiscal pressure, and von Mises stress in annulus fibrosus were analyzed by applying a moment of 4 Nm in flexion, extension, lateral bending, and axial rotation with 100 N compressive loads. Disc height deformation in young healthy, suboptimal healthy, mildly, moderately, and severely degenerative discs was 40%, 37%, 21%, 12%, and 8%, respectively. The decreasing order of the range of motion was young healthy spine > suboptimal healthy spine > mildly degenerative spine > moderately degenerative spine > severely degenerative spine. The mean stress of annulus ground substance in the suboptimal healthy disc was higher than in the young healthy disc. The mean stress of inter-lamellar matrix and annulus ground substance in moderately and severely degenerative discs was higher than in other discs. Age-related structural changes and degenerative changes increased the stiffness and reduced the elastic deformation of intervertebral discs. Decreased range of motion due to the effects of aging or degeneration on the intervertebral disc, may cause compensation of adjacent segments and lead to progressive degeneration of multiple segments. The effect of aging on the intervertebral disc increased the risk of annulus fibrosus damage from the biomechanical point of view. Moderately and severely degenerative discs may have a higher risk of herniation due to the higher risk of damage and layers separation of annulus fibrosus caused by increased stress in the annulus ground substance and inter-lamellar matrix.

Kinematic and biomechanical responses of the spine to distraction surgery in children with early onset scoliosis: A 3-D finite element analysis

Periodical and consecutive distraction is an effective treatment for severe early onset scoliosis (EOS), which enables the spinal coronal and sagittal plane deformity correction. However, the rate of rod fractures and postoperative complications was still high mainly related to the distraction process. Previous studies have primarily investigated the maximum safe distraction force without a rod broken, neglecting the spinal re-imbalance and distraction energy consumption, which is equally vital to evaluate the operative value. This study aimed to reveal the kinematic and biomechanical responses occurring after spinal distraction surgery, which were affected by traditional bilateral fixation. The spinal models (C6-S1) before four distractions were reconstructed based on CT images and the growing rods were applied with the upward displacement load of 0–25 mm at an interval of 5 mm. Relationships between the distraction distance, the distraction force and the thoracic and lumbar Cobb angle were revealed, and the spinal displacement and rotation in three-dimensional directions were measured. The spinal overall imbalance would also happen during the distraction process even under the safe force, which was characterized by unexpected cervical lordosis and lateral displacement. Additionally, the law of diminishing return has been confirmed by comparing the distraction energy consumption in different distraction distances, which suggests that more attention paid to the spinal kinematic and biomechanical changes is better than to the distraction force. Notably, the selection of fixed segments significantly impacts the distraction force at the same distraction distance. Accordingly, some results could provide a better understanding of spinal distraction surgery.

A finite element analysis on different bone cement forms and injection volumes injected into lumbar vertebral body in percutaneous kyphoplasty

Background

To investigate the stress changes between different bone cement forms and injection volumes in adjacent vertebrae after percutaneous kyphoplasty (PKP) by establishing a three-dimensional finite element model of osteoporosis.

Methods

A male healthy volunteer was selected. CT of scans L1 to L3 vertebrae were imported into Mimics 21.0 software.The vertebral model of osteoporosiswas established based on previous literature reference. The models were divided into three groups: unilateral, bilateral integration and bilateral separation groups, with each group injecting 2 ml, 4,ml and 6 ml of bone cement, respectively. In all models, a vertical compressive load of 500 N, anterior flexion/posterior extension, left/right bending, and left/right rotation were applied with a moment of 7.5 N/m, of which 85% was applied to the anterior mid-column and 15% to the posterior column. The stress changes between adjacent vertebrae under different conditions were calculated.

Results

After percutaneous kyphoplasty was applied to the L2 vertebral body, some differences can be found between the effects of different cement injection volumes and cement morphology on adjacent structures. There was no major difference between the groups when the bone cement injection volume was 2 ml. When the amount of bone cement injected was 4 ml, the bone cement morphology of the bilateral integration group (BIG) produced less stress between adjacent vertebral bodies. The minimum stress was 14.95 MPa in the L3 vertebral body in posterior extension. Whereas the stress levels on adjacent intervertebral structures, BIG shaped bone cement shows some superiority. In addition, the adjacent vertebrae and intervertebral structures are subjected to less stress during left and right rotation.

Conclusions

The present finite element study suggested that bilateral integration bone cement is a suitable form of cement injection, and when the injection volume is 4 ml, reduces stress on adjacent segments by approximately 15% while maintaining the stability of the injected vertebral body.

Finite Element Analysis of a New Type of Spinal Protection Device for the Prevention and Treatment of Osteoporotic Vertebral Compression Fractures

OBJECTIVE

To study the effectiveness of a new spinal protection device for preventing and treating osteoporotic vertebral compression fractures (OVCFs) by finite element analysis (FEA).

METHODS

One healthy volunteer and one patient with 1-segment lumbar vertebral compression fractures were included in this experimental study. The DICOM files of two different lumbar spiral computed tomography (CT) scans were converted into STL files, and 3D finite element models of the lumbar spine were generated for normal and L1 vertebral fracture spines. A new type of spinal protection device was applied to reduce the stress on the anterior vertebral edge and direct the center of gravity posteriorly. The stress distribution characteristics of different finite element models of the lumbar spine were analyzed, revealing the characteristics of the stress distributed along the spine under the action of the new spinal protection device.

RESULTS

Under normal conditions, the stress was mainly distributed in the middle and posterior columns of the spine. When the anterior border of the L1 vertebral body was fractured and collapsed, the stress distribution shifted toward the anterior column due to the center of gravity being directed forward. According to finite element analysis of the spine with the new protection device, the stress in the middle and posterior columns tended to increase, and that in the anterior column decreased. After the new type of spinal fixation device was applied, the stress at the L1 and L2 vertebral endplates decreased to a certain extent, especially that at the L1 vertebral body. The maximum stress on the L1 vertebral body decreased by 20% after the auxiliary device was applied.

CONCLUSIONS

According to the FEA results, the new spinal protection device can effectively prevent and treat osteoporotic vertebral compression fractures (OVCFs), and can alter the stress distribution in the spine and reduce the stress in the anterior column of the vertebral body, especially in vertebral compression fractures.

Finite Element Analysis of Unilateral versus Bipedicular Bone-Filling Mesh Container for the Management of Osteoporotic Compression Fractures

The effect of unilateral and bilateral bone-filling mesh containers (BFC) on osteoporotic vertebral compression fracture (OVCF) was analyzed by the finite element method. The CT scan data of the T12-L2 vertebral body were obtained from a healthy female volunteer with no history of lumbar spine injury or obvious abnormality of vertebral body morphology. The normal finite element model of the T12-L2 vertebral body and the finite element model of osteoporosis were established, and the models were validated. The L1 in the normal model of the vertebral body was used to simulate the vertebral compression fracture, after which the unilateral and bilateral BFC were simulated to establish models representing the two surgical approaches. We analyzed changes in the deformation and von Mises stress in vertebral bodies and intervertebral discs in the two models under seven working conditions (axial direction, anteflexion, rear protraction, left-side bending, right-side bending, left rotation, and right rotation) and found that the unilateral and bilateral approaches are biomechanically comparable, with no statistical difference between the two overall models.

Effects of Growing Rod Technique with Different Surgical Modes and Growth Phases on the Treatment Outcome of Early Onset Scoliosis: A 3-D Finite Element Analysis

Early onset scoliosis (EOS) is emerging as a serious threat to children’s health and is the third largest threat to their health after myopia and obesity. At present, the growing rod technique (GRT), which allows patients to regain a well-balanced sagittal profile, is commonly considered as an invasive surgical procedure for the treatment of EOS. However, the risk of postoperative complications and instrumentation breakage remains high, which is mainly related to the choice of fixed mode. Several authors have studied primary stability and instrumentation loads, neglecting the mechanical transmission of the spinal long-segment model in different growth phases, which is fundamental to building a complete biomechanical environment. The present study aimed to investigate the kinematic and biomechanical properties that occur after GRT, across the long spinal structure and the posterior instrumentation, which are affected by unilateral or bilateral fixation. Accordingly, spinal segments (C6-S1) were loaded under flexion (Flex), extension (Ext), left lateral bending (LB), right lateral bending (RB), left torsion (LT), and right torsion (RT) using 11 established spinal models, which were from three growth phases. The stress distribution, spinal and intervertebral range of motion (ROM), counter torque of the vertebra, and bracing force on the rods were measured. The results showed that bilateral posterior fixation (BPF) is more stable than unilateral posterior fixation (UPF), at the expense of more compensations for the superior adjacent segment (SAS), especially when the superior fixed segment is closer to the head. Additionally, the bracing force of the instrumentation on the spine increases as the Cobb angle decreases. Accordingly, this biomechanical analysis provides theoretical suggestions for the selection of BPF or UPF and fixed segments in different growing phases.

The Simulation of Muscles Forces Increases the Stresses in Lumbar Fixation Implants with Respect to Pure Moment Loading

Simplified loading conditions such as pure moments are frequently used to compare different instrumentation techniques to treat spine disorders. The purpose of this study was to determine if the use of realistic loading conditions such as muscle forces can alter the stresses in the implants with respect to pure moment loading. A musculoskeletal model and a finite element model sharing the same anatomy were built and validated against in vitro data, and coupled in order to drive the finite element model with muscle forces calculated by the musculoskeletal one for a prescribed motion. Intact conditions as well as a L1-L5 posterior fixation with pedicle screws and rods were simulated in flexion-extension and lateral bending. The hardware stresses calculated with the finite element model with instrumentation under simplified and realistic loading conditions were compared. The ROM under simplified loading conditions showed good agreement with in vitro data. As expected, the ROMs between the two types of loading conditions showed relatively small differences. Realistic loading conditions increased the stresses in the pedicle screws and in the posterior rods with respect to simplified loading conditions; an increase of hardware stresses up to 40 MPa in extension for the posterior rods and 57 MPa in flexion for the pedicle screws were observed with respect to simplified loading conditions. This conclusion can be critical for the literature since it means that previous models which used pure moments may have underestimated the stresses in the implants in flexion-extension and in lateral bending.

Stability Evaluation of Different Oblique Lumbar Interbody Fusion Constructs in Normal and Osteoporotic Condition - A Finite Element Based Study

Introduction: In developed countries, the age structure of the population is currently undergoing an upward shift, resulting a decrease in general bone quality and surgical durability. Over the past decade, oblique lumbar interbody fusion (OLIF) has been globally accepted as a minimally invasive surgical technique. There are several stabilization options available for OLIF cage fixation such as self-anchored stand-alone (SSA), lateral plate-screw (LPS), and bilateral pedicle screw (BPS) systems. The constructs’ stability are crucial for the immediate and long-term success of the surgery. The aim of this study is to investigate the biomechanical effect of the aforementioned constructs, using finite element analysis with different bone qualities (osteoporotic and normal).

Method: A bi-segmental (L2–L4) finite element (FE) model was created, using a CT scan of a 24-year-old healthy male. After the FE model validation, CAD geometries of the implants were inserted into the L3–L4 motion segment during a virtual surgery. For the simulations, a 150 N follower load was applied on the models, then 10 Nm of torque was used in six general directions (flexion, extension, right/left bending, and right/left rotation), with different bone material properties.

Results: The smallest segmental (L3–L4) ROM (range of motion) was observed in the BPS system, except for right bending. Osteoporosis increased ROMs in all constructs, especially in the LPS system (right bending increase: 140.26%). Osteoporosis also increased the caudal displacement of the implanted cage in all models (healthy bone: 0.06 ± 0.03 mm, osteoporosis: 0.106 ± 0.07 mm), particularly with right bending, where the displacement doubled in SSA and LPS constructs. The displacement of the screws inside the L4 vertebra increased by 59% on average (59.33 ± 21.53%) due to osteoporosis (100% in LPS, rotation). BPS-L4 screw displacements were the least affected by osteoporosis.

Conclusions: The investigated constructs provide different levels of stability to the spine depending on the quality of the bone, which can affect the outcome of the surgery. In our model, the BPS system was found to be the most stable construct in osteoporosis. The presented model, after further development, has the potential to help the surgeon in planning a particular spinal surgery by adjusting the stabilization type to the patient’s bone quality.

Innovative sacropelvic fixation using iliac screws and triangular titanium implants

PURPOSE

Sacropelvic fixation is frequently used in combination with thoracolumbar instrumentation for the correction of severe spinal deformities. The purpose of this study was to explore the effects of the triangular titanium implants on the iliac screw fixation. Our hypothesis was that the use of triangular titanium implants can increase the stability of the iliac screw fixation.

METHODS

Three T10-pelvis instrumented models were created: pedicle screws and rods in T10-S1, and bilateral iliac screws (IL); posterior fixation and bilateral iliac screws and triangular implants inserted bilaterally in a sacro-alar-iliac trajectory (IL-Tri-SAI); posterior fixation and bilateral iliac screws and two bilateral triangular titanium implants inserted in a lateral trajectory (IL-Tri-Lat). Outputs of these models, such as hardware stresses, were compared against a model with pedicle screws and rods in T10-S1 (PED).

RESULTS

Sacropelvic fixation decreased the L5-S1 motion by 75-90%. The motion of the SIJ was reduced by 55-80% after iliac fixation; the addition of triangular titanium implants further reduced it. IL, IL-Tri-SAI and IL-Tri-Lat demonstrated lower S1 pedicle stresses with respect to PED. Triangular implants had a protective effect on the iliac screw stresses.

CONCLUSION

Sacropelvic fixation decreased L5-S1 range of motion suggesting increased stability of the joint. The combination of triangular titanium implants and iliac screws reduced the residual flexibility of the sacroiliac joint, and resulted in a protective effect on the S1 pedicle screws and iliac screws themselves. Clinical studies may be performed to demonstrate applicability of these FEA results to patient outcomes.

Subject-Specific Alignment and Mass Distribution in Musculoskeletal Models of the Lumbar Spine

Musculoskeletal modeling is a well-established method in spine biomechanics and generally employed for investigations concerning both the healthy and the pathological spine. It commonly involves inverse kinematics and optimization of muscle activity and provides detailed insight into joint loading. The aim of the present work was to develop and validate a procedure for the automatized generation of semi-subject-specific multi-rigid body models with an articulated lumbar spine. Individualization of the models was achieved with a novel approach incorporating information from annotated EOS images. The size and alignment of bony structures, as well as specific body weight distribution along the spine segments, were accurately reproduced in the 3D models. To ensure the pipeline’s robustness, models based on 145 EOS images of subjects with various weight distributions and spinopelvic parameters were generated. For validation, we performed kinematics-dependent and segment-dependent comparisons of the average joint loads obtained for our cohort with the outcome of various published in vivo and in situ studies. Overall, our results agreed well with literature data. The here described method is a promising tool for studying a variety of clinical questions, ranging from the evaluation of the effects of alignment variation on joint loading to the assessment of possible pathomechanisms involved in adjacent segment disease.

Lag-Screw Osteosynthesis in Thoracolumbar Pincer Fractures

STUDY DESIGN

Biomechanical.

OBJECTIVE

This study evaluates the biomechanical properties of lag-screws used in vertebral pincer fractures at the thoracolumbar junction.

METHODS

Pincer fractures were created in 18 bisegmental human specimens. The specimens were assigned to three groups depending on their treatment perspective, either bolted, with the thread positioned in the cortical or cancellous bone, or control. The specimens were mounted in a servo-hydraulic testing machine and loaded with a 500 N follower load. They were consecutively tested in 3 different conditions: intact, fractured, and bolted/control. For each condition 10 cycles in extension/flexion, torsion, and lateral bending were applied. After each tested condition, a computed tomography (CT) scan was performed. Finally, an extension/flexion fatigue loading was applied to all specimens.

RESULTS

Biomechanical results revealed a nonsignificant increase in stiffness in extension/flexion of the fractured specimens compared with the intact ones. For lateral bending and torsion, the stiffness was significantly lower. Compared with the fractured specimens, no changes in stiffness due to bolting were discovered. CT scans showed an increasing fracture gap during axial loading both in extension/flexion, torsion, and lateral bending in the control specimens. In bolted specimens, the anterior fragment was approximated, and the fracture gap nullified. This refers to both the cortical and the cancellous thread positions.

CONCLUSION

The results of this study concerning the effect of lag-screws on pincer fractures appear promising. Though there was little effect on stiffness, CT scans reveal a bony contact in the bolted specimens, which is a requirement for bony healing.

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.

Modified poly(methyl methacrylate) bone cement in the treatment of Kümmell disease

Kümmell disease (KD) causes serious vertebral body collapse in patients. However, only a few case reports have been conducted and the number of patients with KD investigated was limited. Additionally, the frequently used poly(methyl methacrylate) (PMMA) bone cement for KD is limited by excessive modulus and poor biocompatibility. Herein, we aimed to modify PMMA bone cement with mineralized collagen (MC), and compare the clinical effects, image performance and finite element analysis between the modified bone cement and PMMA bone cement for the treatment of phase I and II KD. Thirty-nine KD patients treated with PMMA bone cement and 40 KD patients treated with MC-modified PMMA bone cement from June 2015 to March 2017 were retrospectively analyzed. The surgical procedure, intraoperative blood loss, hospital stay and complications were compared between different groups. Visual analog scale, Oswestry disability index, anterior vertebral height, posterior vertebral height, computed tomography value, adjacent vertebral re-fracture, Cobb angle and wedge-shaped correction angle were evaluated. Additionally, the representative sample was selected for finite element analysis. We found that the MC-modified PMMA bone cement could achieve the same effect as that of PMMA bone cement and was associated with better vertebral height restoration in the long term.

Sensitivity analysis of biomechanical effect in vertebral body of two different augmenters

BACKGROUND

Transvertebral Bone Graft and Augmentation (TBGA) has achieved good clinical effects in the treatment of osteoporotic vertebral compression fractures (OVCFs). This study aimed to investigate the postoperatively biomechanical effects of TBGA and compare the biomechanical sensitivity of two different augmenters: a cylindrical enhancement device (CED) and bone cement.

METHODS

Finite element models of the spine segment T11-L3 were created, including one model based on normal segment and the other three with L1 augmentation for pathological conditions. Three treatments were simulated including CED implant treatment A, CED implant treatment B, and bone cement treatment. The stress distribution and maximum displacement of the four models under different treatments were analyzed. A method of linear fitting of dummy variables was used to analyze the sensitivity of biomechanical parameters to the degree of osteoporosis (DO) and load.

FINDINGS

The reduction of stress with increasing DO in augmented and adjacent vertebral bodies under bone cement augmentation was less than that under CED augmentation. The stress of augmented vertebral body and the adjacent vertebral body was most sensitive to extension and rotation loading conditions. As DO increasing, the bone cement augmentation significantly increased the stress level on the upper and lower endplates.

INTERPRETATION

When the degree of osteoporosis increased, CED outperforms bone cement in terms of the stress reduction in augmented vertebral and adjacent vertebral, which could be beneficial for avoiding re-fracture. Using TBGA to treat OVCFs, especially with Plan B method, the condition of the pathological spine is closer to the original status in terms of the sensitivity to stress and the spinal range of motion. The TBGA treatment is sensitive to lateral bending and torsion, therefore patients should be advised to avoid high-risk motions like lateral bending and rotation.

Assessment of trunk muscle activation and intervertebral load in adolescent idiopathic scoliosis by musculoskeletal modelling approach

Adolescent idiopathic scoliosis (AIS) is a three-dimensional deformity of the spine, the aetiology and pathogenesis of which are poorly understood. Unfortunately, biomechanical data describing trunk muscle activation and intervertebral load, which can contribute to understanding the pathomechanics of the AIS spine, cannot be measured in vivo due to the invasiveness of the procedures. The present study provides the biomechanical characterization of the spinal loads in scoliotic subjects by exploiting musculoskeletal modelling approach, allowing for calculating biomechanical measures in an assigned posture. A spine model with articulated ribcage previously developed in AnyBody software was applied. The predicted outcomes were evaluated in the upright posture, depending on scoliosis severity and curve type, in a population of 132 scoliotic subjects with mild, moderate, and severe scoliosis. Radiographic-based three dimensional reconstruction of vertebral orientations and scaling of body segments and trunk muscle cross-section area guaranteed geometrical subject-specificity. Validation analysis supporting the application of the model was performed. Trunk muscles were found more activated in the convex side of the scoliotic curve, in agreement with reference in vivo measurements, with progressive increase with scoliosis severity. The intervertebral lateral shear was found positively correlated with the severity of the scoliosis, demonstrating that the transferred load is not a priori orthogonal to vertebral endplate in the frontal plane, and thus questioning the assumption of the 'follower load' approach in case of experimental or computational study on the scoliotic spine. The study opens the way for the subject-specific characterization of scoliosis in assigned loading and motion conditions.

Development of a multiscale model of the human lumbar spine for investigation of tissue loads in people with and without a transtibial amputation during sit-to-stand

Quantification of lumbar spine load transfer is important for understanding low back pain, especially among persons with a lower limb amputation. Computational modeling provides a helpful solution for obtaining estimates of in vivo loads. A multiscale model was constructed by combining musculoskeletal and finite element (FE) models of the lumbar spine to determine tissue loading during daily activities. Three-dimensional kinematic and ground reaction force data were collected from participants with ([Formula: see text]) and without ([Formula: see text]) a unilateral transtibial amputation (TTA) during 5 sit-to-stand trials. We estimated tissue-level load transfer from the multiscale model by controlling the FE model with intervertebral kinematics and muscle forces predicted by the musculoskeletal model. Annulus fibrosis stress, intradiscal pressure (IDP), and facet contact forces were calculated using the FE model. Differences in whole-body kinematics, muscle forces, and tissue-level loads were found between participant groups. Notably, participants with TTA had greater axial rotation toward their intact limb ([Formula: see text]), greater abdominal muscle activity ([Formula: see text]), and greater overall tissue loading throughout sit-to-stand ([Formula: see text]) compared to able-bodied participants. Both normalized (to upright standing) and absolute estimates of L4-L5 IDP were close to in vivo values reported in the literature. The multiscale model can be used to estimate the distribution of loads within different lumbar spine tissue structures and can be adapted for use with different activities, populations, and spinal geometries.

The use of triangular implants to enhance sacropelvic fixation: a finite element investigation

BACKGROUND CONTEXT

Long thoracolumbar fixation and fusion have become a consolidated treatment for severe spinal disorders. Concomitant sacropelvic fixation with S2 alar-iliac (S2AI) screws is frequently performed to limit instrumentation failure and pseudarthrosis at the lumbosacral junction.

PURPOSE

This study explored the use of triangular titanium implants in different configurations in which the implants supplemented standard sacropelvic fixation with S2AI screws in order to further increase the stability of S2AI fixation.

STUDY DESIGN

Finite element study.

METHODS

Four T10-pelvis instrumented models were built: pedicle screws and rods in T10-S1 (PED); pedicle screws and rods in T10-S1, and bilateral S2 alar-iliac screws (S2AI); pedicle screws and rods in T10-S1, bilateral S2AI screws, and triangular implants inserted bilaterally in a sacral alar-iliac trajectory (Tri-SAI); pedicle screws and rods in T10-S1, bilateral S2AI screws and two bilateral triangular titanium implants inserted in a lateral trajectory (Tri-Lat). The models were tested under pure moments of 7.5 Nm in flexion-extension, lateral bending and axial rotation.

RESULTS

SIJ motion was reduced by 50% to 66% after S2AI fixation; the addition of triangular titanium implants in either a SAI or a lateral trajectory further reduced it. S2AI, Tri-SAI, and Tri-Lat resulted in significantly lower stresses in S1 pedicle screws when compared to PED. Triangular implants had a protective effect on the maximal stresses in S2AI screws, especially when placed in the SAI trajectory. Sacropelvic fixation did not have any protective effect on the posterior rods.

CONCLUSIONS

Supplementing S2AI screws with triangular implants had a protective effect on the S2AI screws themselves, as well as the S1 pedicle screws, in the tested model.

CLINICAL SIGNIFICANCE

Triangular implants can substantially reduce the residual flexibility of the SIJ with respect to S2AI fixation alone, suggesting a possible role in patients needing reinforced fixation. In vivo investigation is needed to determine if these in vitro effects translate into clinically important differences.

Muscle-driven and torque-driven centrodes during modeled flexion of individual lumbar spines are disparate

Lumbar spine biomechanics during the forward-bending of the upper body (flexion) are well investigated by both in vivo and in vitro experiments. In both cases, the experimentally observed relative motion of vertebral bodies can be used to calculate the instantaneous center of rotation (ICR). The timely evolution of the ICR, the centrode, is widely utilized for validating computer models and is thought to serve as a criterion for distinguishing healthy and degenerative motion patterns. While in vivo motion can be induced by physiological active structures (muscles), in vitro spinal segments have to be driven by external torque-applying equipment such as spine testers. It is implicitly assumed that muscle-driven and torque-driven centrodes are similar. Here, however, we show that centrodes qualitatively depend on the impetus. Distinction is achieved by introducing confidence regions (ellipses) that comprise centrodes of seven individual multi-body simulation models, performing flexion with and without preload. Muscle-driven centrodes were generally directed superior–anterior and tail-shaped, while torque-driven centrodes were located in a comparably narrow region close to the center of mass of the caudal vertebrae. We thus argue that centrodes resulting from different experimental conditions ought to be compared with caution. Finally, the applicability of our method regarding the analysis of clinical syndromes and the assessment of surgical methods is discussed.

Design of Patient Specific Spinal Implant (Pedicle Screw Fixation) using FE Analysis and Soft Computing Techniques

BACKGROUND

This work uses genetic algorithm (GA) for optimum design of patient specific spinal implants (pedicle screw) with varying implant diameter and bone condition. The optimum pedicle screw fixation in terms of implant diameter is on the basis of minimum strain difference from intact (natural) to implantation at peri-prosthetic bone for the considered six different peri-implant positions.

METHODS

This design problem is expressed as an optimization problem using the desirability function, where the data generated by finite element analysis is converted into an artificial neural network (ANN) model. The finite element model is generated from CT scan data. Thereafter all the ANN predictions of the microstrain in six positions are converted to unitless desirability value varying between 0 and 1, which is then combined to form the composite desirability. Maximization of the composite desirability is done using GA where composite desirability should be made to go up as close as possible to 1. If the composite desirability is 1, then all 'strain difference values in 6 positions' are 0.

RESULTS

The optimum solutions obtained can easily be used for making patient-specific spinal implants.

A COMPARATIVE ANALYSIS OF THE EFFECTS OF TWO AUGMENTERS ON THE SENSITIVITY OF VERTEBRAL BIOMECHANICAL BEHAVIOR IN VERTEBROPLASTY

The treatment of osteoporotic vertebral compression fractures (OVCFS) by transvertebral bone graft and augmentation (TBGA) has achieved satisfactory clinical results, but its biomechanical effects are not clear. The purpose of this study was to investigate the biomechanical effects of TBGA and compare the biomechanical sensitivity of the augmenter used in TBGA — a cylindrical enhancement device (CED) with bone cement. The finite element (FE) model of healthy segments T11-L3 (M1) was built, and two other models with L1 augmentation (M2, M3) were established to simulate CED and bone cement treatment, respectively. The stress and displacement distribution of the three models under five physiological loads were calculated and analyzed by the FE method. Based on the results, the sensitivities of biomechanical parameters to the degree of osteoporosis (DO) and loads were analyzed by linear fitting method using dummy variables. With the increase of DO, the CED is superior to bone cement in preventing the fractures of the augmented vertebral and the adjacent vertebral under the set loading conditions. Simulating TBGA method, the model 2 with L1 reconstructed was closer to the normal T11-L3 model in terms of sensitivity of stress and displacement under different loading conditions.

Biomechanical analysis of sagittal correction parameters for surgical instrumentation with pedicle subtraction osteotomy in adult spinal deformity

BACKGROUND

High complication rate has been documented in spinal instrumentation with pedicle subtraction osteotomy. The objective was to analyze the biomechanics of spinal instrumentation with pedicle subtraction osteotomy as functions of three instrumentation parameters.

METHODS

Patient-specific biomechanical computer models of 3 adult patients who were instrumented with pedicle subtraction osteotomy were used to simulate the actual instrumentations and alternative instrumentations derived by varying instrumentation parameters: pedicle subtraction osteotomy wedge angle and rod contouring angle, location of the osteotomy, and number of rods (2 vs. 4).

FINDINGS

A change of the sagittal correction by +7.5° (-7.5°) resulted in a change in the screw axial forces and rods bending loads around the osteotomy by +38% (-19%) and +28% (-11%), respectively. The bending moments in the rods were 31% lower at the osteotomy site when it was located at one level above, and 20% higher when it was located at a level below. Additional rods allowed the rod bending and screw axial loads to be reduced by 24% and 22%, respectively.

INTERPRETATION

The amount of sagittal correction was positively correlated with loads sustained by the screws and rods. Rods were subjected to higher bending loads at the pedicle subtraction osteotomy site when it was done at a lower level. A 4-rod construct is an effective way to reduce the risk of rod breakage by reducing the loads sustained by the rods.

Biomechanics of sacropelvic fixation: a comprehensive finite element comparison of three techniques

PURPOSE

Sacropelvic fixation is frequently used in combination with thoracolumbar instrumentation for complex deformity correction and is commonly associated with pseudoarthrosis, implant failure and loosening. This study compared pedicle screw fixation (PED) with three different sacropelvic fixation techniques, namely iliac screws (IL), S2 alar-iliac screws (S2AI) and laterally placed triangular titanium implants (SI), all in combination with lumbosacral instrumentation, accounting for implant micromotion.

METHODS

Existing finite element models of pelvis-L5 of three patients including lumbopelvic instrumentation were utilized. Moments of 7.5 Nm in the three directions combined with a 500 N compressive load were simulated. Measured metrics included flexibility, instrumentation stresses and bone-implant interface loads.

RESULTS

Fixation effectively reduced the sacroiliac flexibility. Compared to PED, IL and S2AI induced a reduction in peak stresses in the S1 pedicle screws. Rod stresses were mostly unaffected by S2AI and SI, but IL demonstrated a stress increase. In comparison with a previous work depicting full osteointegration, SI was found to have similar instrumentation stresses as those due to PED.

CONCLUSIONS

Fixation with triangular implants did not result in stress increase on the lumbosacral instrumentation, likely due to the lack of connection with the posterior rods. IL and S2AI had a mild protective effect on S1 pedicle screws in terms of stresses and bone-implant loads. IL resulted in an increase in the rod stresses. A comparison between this study and previous work incorporating full osteointegration demonstrates how these results may be applied clinically to better understand the effects of different treatments on patient outcomes. These slides can be retrieved under Electronic Supplementary Material.

A comparison of rigid, semi-rigid and flexible spinal stabilization devices: A finite element study

Pedicle-screw-based spinal fixation system has shown its success for treating degenerative disc disease-related back pain and spinal instability. In the last few decades, several non-fusion implants ( 'flexible' or 'dynamic' fixation) are developed for treating slight degenerate disc disease. The aim of this study is to characterize and compare the biomechanical responses of pedicle-screw-based fusion with various rod materials and a flexible spinal stabilization device on the lumbar spine (L3L5). Computed tomography scan-based finite element model and pedicle screw fixation with rigid rod material stainless steel rod, semi-rigid rod material poly-ether-ether-ketone and flexible rod device made of stainless steel are used in this study. Intact model of the lumbar spine and treated with all the different implants are simulated under typical physiological loading conditions. Compared with the intact model, pedicle screw with the stainless steel rod fixation system is found to offer very less range of motion. Poly-ether-ether-ketone rod system increased range of motion 3.8, 7 and 1.8 times for axial rotation, lateral bending and flexion-extension, respectively, compared to the stainless steel rod system. The flexible rod device rod system is found to reduce stress on vertebral body, carrying out more loads as compared to poly-ether-ether-ketone rods. In the case of stainless steel rods, range of motion is almost restricted on the fusion zone, which is overcome by the poly-ether-ether-ketone rod system to some extent and farther improved by the flexible rod device rod system. So, the poly-ether-ether-ketone rod and flexible rod device rod systems may be implemented for better clinical results after succeeding experimental validation and clinical trial.

Finite element analysis of a ball-and-socket artificial disc design to suppress excessive loading on facet joints: A comparative study with ProDisc

Facet arthrosis at surgical level was identified as major complication after total disc replacement (TDR). One of the reasons for facet arthrosis after TDR has been speculated to be the hypermobility of artificial discs. Accordingly, the artificial disc that can constrain the hypermobility of ball-and-socket type artificial discs and reduce loading on facet joints is demanded. The proposed artificial disc, which is named as NewPro, was constructed based on the FDA-approved ProDisc but contained an interlocking system consisting of additional bars and grooves to control the range of motion (ROM) of lumbar spine in all anatomical planes. The three-dimensional finite element model of L1 to L5 was developed first, and the biomechanical effects were compared between ProDisc and NewPro. The ROM and facet contact force of NewPro were significantly decreased by 42.7% and 14% in bending and by 45.6% and 34.4% in torsion, respectively, compared with the values of ProDisc, thanks to the interlocking system. In addition, the ROM and facet contact force could be selectively constrained by modifying the location of the bars. The proposed artificial disc with the interlocking system was able to constrain the intersegmental rotation effectively and reduce excessive loading on facet joints, although wear and strength tests would be needed prior to clinical applications.

The Spine: A Strong, Stable, and Flexible Structure with Biomimetics Potential

From its first appearance in early vertebrates, the spine evolved the function of protecting the spinal cord, avoiding excessive straining during body motion. Its stiffness and strength provided the basis for the development of the axial skeleton as the mechanical support of later animals, especially those which moved to the terrestrial environment where gravity loads are not alleviated by the buoyant force of water. In tetrapods, the functions of the spine can be summarized as follows: protecting the spinal cord; supporting the weight of the body, transmitting it to the ground through the limbs; allowing the motion of the trunk, through to its flexibility; providing robust origins and insertions to the muscles of trunk and limbs. This narrative review provides a brief perspective on the development of the spine in vertebrates, first from an evolutionary, and then from an embryological point of view. The paper describes functions and the shape of the spine throughout the whole evolution of vertebrates and vertebrate embryos, from primordial jawless fish to extant animals such as birds and humans, highlighting its fundamental features such as strength, stability, and flexibility, which gives it huge potential as a basis for bio-inspired technologies.

Differential response to vibration of three forms of scoliosis during axial cyclic loading: a finite element study

Background

Scoliosis is a serious disease that can affect all segments of society. Few studies have investigated the response to vibration of differing sinusoidal axial cyclic loading frequencies for different forms of scoliosis in the lumbar spine.

Methods

In this study, four finite element models, comprising a healthy spine, Lenke-A, Lenke-B and Lenke-C scoliosis of the lumbar S1-L1 region were developed. Modal analysis extracted resonant frequencies of the FE models with an upper body mass of 40 kg and 400 N preload. A transient dynamic analysis was performed to obtain the response to vibration of models under a sinusoidal axial loading of ± 40N at frequencies of 3, 5, 7, 9, 11 and 13 Hz using an upper body mass of 40 kg and 400 N preload.

Results

The first-order resonant frequencies of healthy, Lenke-A, Lenke-B and Lenke-C spines were 9.2, 3.9, 4.6 and 5.7 Hz, respectively. A Lenke-A lumbar spine was more likely to deform at a lower vibration frequency and Lenke-C deformed more easily at a higher vibration frequency. Furthermore, the vibration amplitude in the Y-direction (left-right) was greatest and least in the Z-direction (top-bottom). The frequency of cyclic loading closest to the resonant frequency resulted in a maximum value of peak-to-peak vibrational displacement. Furthermore, the vibrational amplitudes in patients with scoliosis were larger than they were in healthy subjects. In addition, axial displacement of the vertebrae in the healthy spine changed steadily whereas fluctuations in the scoliotic vertebrae in scoliosis patients were greater than that of other vertebrae.

Conclusions

Different forms of scoliosis may have different vibrational characteristics, the scoliotic vertebrae being the weak link in scoliosis under loading condition of whole body vibration. Scoliosis was more sensitive to this form of vibration. Where the frequency of axial cyclic vibrational loading of the lumbar spine was closer to its resonant frequency, the vibrational amplitude was larger. These results suggest that vibration will exacerbate the degree of scoliosis and so such patients should reduce their exposure to vibration. Clinical treatment should pay attention to the scoliotic vertebrae and reduce their vibration. These findings may assist in the clinical prevention and treatment of scoliosis.

Development of a Biconcave Mobile-Bearing Lumbar Total Disc Arthroplasty Concept Using Finite Element Analysis and Design Optimization

Total disc arthroplasty (TDA) is a motion-preserving surgical treatment for spinal disorders. However, adverse surgical outcomes, such as abnormal kinematics, facet joint (FJ) overloading, and polyethylene (PE) failures, have limited wide application of lumbar TDAs. The objectives of this computational study were to elucidate how implant design and FJ articulation both influence metal-on-polyethylene (MoP) motion and contact mechanics, as well as to propose and refine a new mobile-bearing TDA concept which enhanced postoperative performance. Simulation results show that abnormal motions (lift-off and/or unsymmetrical motion) are alleviated in fixed-/mobile-bearing TDA-treated segments, as the FJ gap increases. It clearly demonstrates that FJ articulation guides segmental motion and interferes with intended MoP articulation. For an existing biconvex mobile-bearing design, component impingement leads to a peak PE stress of 20.8 MPa (yield stress: 13 MPa), indicating a high risk of PE creep/fracture. Therefore, we proposed a new TDA concept featuring a biconcave PE core with a smooth shape, in order to strengthen the PE rim and mitigate edge-loading. Furthermore, the biconcave-core TDA was optimally designed to promote normal segmental range of motion (ROM), or to minimize polyethylene contact pressure (PCP). In extension (the severest loading scenario), the biconvex-core TDA design caused a ROM 3.6° (+88%) greater than the intact segment and a peak PCP of 116.5 MPa. In contrast, ROM-optimal or PCP-optimal biconcave-core TDA designs decreased the ROM difference to 0.0° or the peak PCP to 24.3 MPa. Therefore, this new TDA design can potentially reduce the incidence of hypermotion and PE damage. © 2019 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:1805-1816, 2019.

Effects of standing on lumbar spine alignment and intervertebral disc geometry in young, healthy individuals determined by positional magnetic resonance imaging

BACKGROUND

Most diagnostic imaging of the spine is performed in supine, a relatively unloaded position. Although the spine is subjected to functional loading that changes the spinal alignment and intervertebral disc geometry, little data exists on how healthy spines adapt to standing. This study seeks to quantify the changes of the lumbar spine from supine to standing in young, back-healthy individuals using a positional magnetic resonance imaging system.

METHODS

This is an observational study that examined the changes in the lumbar spine alignment and intervertebral disc geometry between supine and standing of forty participants (19 males/21 females) without a history of low back pain. The regional lumbar spinal alignment was measured by the sagittal Cobb angle. Segmental intervertebral disc measurements included the segmental Cobb angle, anterior-to-posterior height ratio, and intervertebral disc width measured at L1/L2 - L5/S1 levels. Intra-class correlation was performed for intra- and inter-observer measurements.

FINDINGS

The intra-observer intra-class correlation consistency model ranged from 0.76 to 0.98 with the inter-observer correlation ranging from 0.68 to 0.99. The Cobb angle decreased in standing. The L5/S1 segmental Cobb angle decreased in standing. The L2/L3 and L3/L4 anterior-to-posterior height ratios increased and the L5/S1 anterior-to-posterior height ratio decreased in standing. No difference in intervertebral disc widths was observed from supine to standing.

INTERPRETATIONS

We established normative data for a back-healthy population, using a positional magnetic resonance imaging system, that could inform future investigations that examine the standing-induced adaptations of the lumbar spine in individuals with spinal or intervertebral disc pathologies.

Development and validation of a geometrically personalized finite element model of the lower ligamentous cervical spine for clinical applications

Epidemiological and clinical studies show that the magnitude and scope of cervical disease are on the rise, along with the world's rising aging population. From a biomechanical perspective, the cervical spine presents a wide inter-individual variability, where its motion patterns and load sharing strongly depend on the anatomy. This study aimed to first develop and validate a geometrically patient-specific model of the lower cervical spine for clinical applications, and secondly to use the model to investigate the spinal biomechanics associated with typical cervical disorders. Based on measurements of 30 parameters from X-ray radiographs, the 3D geometry of the vertebrae and intervertebral discs (IVDs) were developed, and detailed finite element models (FEMs) of the lower ligamentous cervical spine for 6 subjects were constructed and simulated. The models were then used for the investigation of different grades of IVD alteration. The multi directional range of motion (ROM) results were in alignment with the in-vitro and in-Silico studies confirming the validity of the model. Severe disc alteration (Grade 3) presented a significant decrease in the ROM and intradiscal pressure (flexion, extension, and axial rotation) on the C5-C6 and slightly increase on the adjacent levels. Maximum stress in Annulus Fibrosus (AF) and facet joint forces increased for Grade 3 for both altered and adjacent levels. The novel validated geometrically-personalized FEM presented in this study potentially offers the clinical community a valuable quantitative tool for the noninvasive analyses of the biomechanical alterations associated with cervical spine disease towards improved surgical planning and enhanced clinical outcomes.

Is the sheep a suitable model to study the mechanical alterations of disc degeneration in humans? A probabilistic finite element model study

Intervertebral disc degeneration is one major source of low back pain, which because of its complex multifactorial nature renders the treatment challenging and thus necessitates extensive research. Experimental animal models have proven valuable in improving our understanding of degenerative processes and potentially promising therapies. Currently, the sheep is the most frequently used large animal in vivo model in intervertebral disc research. However, despite its undoubted value for investigations of the complex biological and cellular aspects, to date, it is unclear whether the sheep is also suited to study the mechanical aspects of disc degeneration in humans. A parametric finite element (FE) model of the L4-5 spinal motion segment was developed. Using this model, the geometry and the material properties of both the human and the ovine spinal segment as well as different appearances of disc degeneration can be depicted. Under pure and combined loads, it was investigated whether degenerative changes to both the human and the ovine model equivalent caused the same mechanical response. Different patterns of degeneration resulted in large variations in the ranges of motion, intradiscal pressure, ligament and facet loads. In the human, but not in the ovine model, all these results differed significantly between different degrees of degeneration. This FE model study highlighted possible differences in the mechanical response to disc degeneration between human and ovine intervertebral discs and indicates the necessity of further, more detailed, investigations.

The Path to Deliver the Most Realistic Follower Load for a Lumbar Spine in Standing Posture: A Finite Element Study

A spine is proven to be subjected to a follower load which is a compressive load of physiologic magnitude acting on the whole spine. The path of the follower load approximates the tangent to the curve of the spine in in vivo neutral standing posture. However, the specific path location of the follower load is still unclear. The aim of this study is to find out the most realistic location of the follower load path (FLP) for a lumbar spine in standing. A three-dimensional (3D) nonlinear finite element model (FEM) of lumbosacral vertebrae (L1-S1) with consideration of the calibrated material properties was established and validated by comparing with the experimental data. We show that the shape of the lumbosacral spine is strongly affected by the location of FLP. An evident nonlinear relationship between the FLP location and the kinematic response of the L1-S1 lumbosacral spine exists. The FLP at about 4 and 3 mm posterior to the curve connecting the center of the vertebral bodies delivers the most realistic location in standing for healthy people and patients having low back pains (LPBs), respectively. Moreover, the “sweeping” method introduced in this study can be applicable to all individualized FEM to determine the location of FLP.

Biomechanical effects of USS fixation with different screw insertion depths on the vertebrae stiffness and screw stress for the treatment of the L1 fracture

OBJECTIVE

To evaluate the biomechanical effects of internal fixation with different screw insertion depths on vertebrae stiffness and screw stress for L1 fracture.

METHODS

The established L1 fracture was fixed with 10 different depths of screw insertion: 10-100% screw-path length (SPL). Loading on the T12 endplate was simulated.

RESULTS

Screws inserted to 60-100% depths has a higher axial displacement of screw against injured vertebrae and maximum stress of screws compared to those of screws inserted to 30-50% depths and 10-20% (P< 0.05). No significant difference was noted among 60-100% SPL groups. Under single loading condition, the incidence rate of maximum stress of each screw ranged from 16.7-50.0%. Chi-square test showed superior screw has a higher incidence rate of maximum stress than inferior screw (P< 0.05).

CONCLUSIONS

Screws inserted to 60% depth or more can achieve effective strength to withstand the postoperative height correction loss of the L1 vertebrae fracture. However, continuous prolonged depth of screw insertion did not significantly increase the effective strength of the screw against injured vertebrae and maximum equivalent stress of screws. The incidence rate of the maximum stress of each screw in correlated with position of screw insertion but not associated with the screw insertion depth.

Impact of material properties of intervertebral disc on dynamic response of the human lumbar spine to vertical vibration: a finite element sensitivity study

This study aimed to determine the effect of variations in material properties of the intervertebral disc on dynamic response of the human lumbar spine to vertical vibration using a finite element model of the lumbar L1-S1 motion segment. The present material sensitivity study was conducted by varying elastic moduli for annulus ground substance (AGS), annulus fibers (AF), and nucleus pulposus (NP) in the disc. Transient dynamic analysis was performed initially on the model with basic material property under a sinusoidal vertical vibration load. Subsequently, the same analysis was done for each of the three disc components corresponding to high and low material property cases. The computed results were plotted as a function of time and compared. The AGS property displayed a larger impact on vertebral axial displacement and von Mises stress in AGS, and the AF property displayed a larger impact on disc bulge. In contrast, the NP property had little effect on all the response parameters. Additionally, the intradiscal pressure was found to be not sensitive to any of the disc properties. These findings may be helpful in adoption of appropriate material parameters for the intervertebral disc in finite element model of the lumbar spine used for vibration analysis. Graphical abstract Material property sensitivity analysis on vibration characteristics of the human lumbar spine.

Biomechanical analysis of spino-pelvic postural configurations in spondylolysis subjected to various sport-related dynamic loading conditions

PURPOSE

To study the risks of spondylolysis due to extrinsic loading conditions related to sports activities and intrinsic spino-pelvic postural parameters [pelvic incidence (PI) and sacral slope (SS)].

METHODS

A comprehensive osseo-disco-ligamentous L4-S1 finite element model was built for three cases with spondylolysis representing three different spino-pelvic angular configurations (SS = 32°, 47°, 59° and PI = 49°, 58°, 72°, respectively). After simulating the standing posture, 16 dynamic loading conditions were computationally tested for each configuration by combining four sport-related loads (compression, sagittal and lateral bending and axial torque). For each simulation, the Von Mises stress, L5-S1 facet contact force and resultant internal loads at the sacral endplate were computed. Significant effects were determined with an ANOVA.

RESULTS

The maximal stress and volume of cancellous bone in the pars with stress higher than 75% of the ultimate stress were higher with 900 N simulated compression (2.2 MPa and 145 mm³) compared to only the body weight (1.36 MPa and 20.9 mm³) (p < 0.001). Combined compression with 10 Nm of flexion and an axial torque of 6 Nm generated the highest stress conditions (up to 2.7 MPa), and L5-S1 facet contact force (up to 430 N). The maximal stress was on average 17% higher for the case with the highest SS compared to the one with lowest SS for the 16 tested conditions (p = 0.0028).

CONCLUSIONS

Combined flexion and axial rotation with compression generated the highest stress conditions related to risks of spondylolysis. The stress conditions intensify in patients with higher PI and SS. These slides can be retrieved under Electronic Supplementary Material.