Effects of eight different ligament property datasets on biomechanics of a lumbar L4-L5 finite element model

Comparative finite element analysis of posterior short segment fixation constructs with or without intermediate screws in the fractured vertebrae for the treatment of type a thoracolumbar fracture

Six-screw short-segment posterior fixation for thoracolumbar fractures, which involves intermediate screws at the fractured vertebrae has been proposed to reduce the rates of kyphosis recurrence and implant failure. Yet, little is known about the mechanisms and biomechanical responses by which intermediate screws at the fracture vertebrae enhance fixation strength. The objective of this study was to investigate the biomechanical properties that are associated with the augmentation of intermediate screws in relation to the severity of type A thoracolumbar fracture using finite element analysis. Short-segment stabilization models with or without augmentation screws at fractured vertebrae were established based on finite element model of moderate compressive fractures, severe compressive fractures and burst fractures. The spinal stiffness, stresses at the implanted hardware, and axial displacement of the bony defect were measured and compared under mechanical loading conditions. All six-screw stabilization showed a decreased range of motion in extension, lateral bending, and axial rotation compared to the traditional four-screw fixation models. Burst thoracolumbar fracture benefited more from augmentation of intermediate screws at the fracture vertebrae. The stress of the rod in six-screw models increased while decreased that of pedicle screws. Our results suggested that patients with more unstable fractures might achieve greater benefits from augmentation of intermediate screws at the fracture vertebrae. Augmentation of intermediate screws at the fracture vertebrae is recommended for patients with higher wedge-shaped or burst fractures to reduce the risk of hardware failure and postoperative re-collapse of injured vertebrae.

Cumulative creep response of viscoelastic lumbar tissue as a function of work-rest schedule

We explore the effect of stress-recovery schedule on the cumulative creep response of lumbar tissues. Twelve participants performed a 48-minute protocol that consisted of 12 min of full trunk flexion and 36 min of upright standing. Two stress-recovery (work-rest) schedules were considered: a) three minutes of full trunk flexion followed by twelve minutes of upright standing (3:12), and b) one minute of full trunk flexion followed by four minutes of upright standing (1:4). Lumbar kinematics and EMG activity of erector spinae muscles were collected. Cumulative creep deformation was explored by considering the changes in peak lumbar flexion angles during full flexion and changes in the angles of flexion-relaxation (EMG-off) of the lumbar extensor musculature after the 48-minute protocol. The results of time-dependent lumbar flexion angle during full flexion revealed a noticeable creep response in both work-rest schedules, but the cumulative creep response was significantly greater in the 3:12 schedule (Δ3.5°) than in the 1:4 schedule (Δ1.6°). Similarly, the change in the EMG-off lumbar flexion angle in the 3:12 schedule was significantly greater than in the 1:4 schedule (Δ2.5° vs -Δ0.2°, respectively). These results indicate that the passive lumbar tissues recover their force producing capability more rapidly with shorter cycle times.

Biomechanical evaluation of different oblique lumbar interbody fusion constructs: a finite element analysis

BACKGROUND

Finite element analysis (FEA) was performed to investigate the biomechanical differences between different adjunct fixation methods for oblique lumbar interbody fusion (OLIF) and to further analyze its effect on adjacent segmental degeneration.

METHODS

We built a single-segment (Si-segment) finite element model (FEM) for L4-5 and a double-segment (Do-segment) FEM for L3-5. Each complete FEM was supplemented and modified, and both developed two surgical models of OLIF with assisted internal fixation. They were OLIF with posterior bilateral percutaneous pedicle screw (TINA system) fixation (OLIF + BPS) and OLIF with lateral plate system (OLIF + LPS). The range of motion (ROM) and displacement of the vertebral body, cage stress, adjacent segment disc stress, and spinal ligament tension were recorded for the four models during flexion/extension, right/left bending, and right/left rotation by applying follower load.

RESULTS

For the BPS and LPS systems in the six postures of flexion, extension, right/left bending, and right/left rotation, the ROM of L4 in the Si-segment FEM were 0.32°/1.83°, 0.33°/1.34°, 0.23°/0.47°, 0.24°/0.45°, 0.33°/0.79°, and 0.34°/0.62°; the ROM of L4 in the Do-segment FEM were 0.39°/2.00°, 0.37°/1.38°, 0.23°/0.47°, 0.21°/0.44°, 0.33°/0.57°, and 0.31°/0.62°, and the ROM of L3 in the Do-segment FEM were 6.03°/7.31°, 2.52°/3.50°, 4.21°/4.38°, 4.21°/4.42°, 2.09°/2.32°, and 2.07°/2.43°. BPS system had less vertebral displacement, less cage maximum stress, and less spinal ligament tension in Si/Do-segment FEM relative to the LPS system. BPS system had a smaller upper adjacent vertebral ROM, greater intervertebral disc stress in terms of left and right bending as well as left and right rotation compared to the LPS system in the L3-4 of the Do-segment FEM. There was little biomechanical difference between the same fixation system in the Si/Do-segment FEM.

CONCLUSIONS

Our finite element analysis showed that compared to OLIF + LPS, OLIF + BPS (TINA) is more effective in reducing interbody stress and spinal ligament tension, and it better maintains the stability of the target segment and provides a better fusion environment to resist cage subsidence. However, OLIF + BPS (TINA) may be more likely to cause adjacent segment degeneration than OLIF + LPS.

Finite element modelling and biodynamic response prediction of the seated human body exposed to whole-body vibration

Biodynamic modelling of seat-occupant systems can assist in seat comfort design. A finite element (FE) model of the seated human body, including detailed modelling of the lumbar spine, was established to reflect the human response to vibration and biodynamic response of the lumbar spine under whole-body vibration (WBV). The lumbar spine model was established and validated against the in-vitro results and calculated data. The posture of the lumbar spine was adjusted according to the radiological research results, and the adjusted model was combined to establish a FE model of the seated human body. The present seated human model with backrest inclination angles of 10, 20, and 30°, validated by comparing the measured apparent mass and seat-to-lumbar spine transmissibility, was used to calculate the biodynamic response of the lumbar spine with three inclined backrests under WBV. The results showed that the model could characterise the apparent mass, seat-to-lumbar spine transmissibility, and the biodynamic response of the lumbar spine. Practitioner summary: Biodynamic models can represent dynamic characteristics of the human body exposed to vibration and assist in seat comfort design. The three-dimensional FE model of the human body can be used to explore the human response to vibration and the biodynamic response of the lumbar spine under WBV.

The Impact of Spine Pathology on Posterior Ligamentous Complex Structure and Function

PURPOSE OF REVIEW

Spinal ligament is an important component of the spinal column in mitigating biomechanical stress. Particularly the posterior ligamentous complex, which is composed of the ligamentum flavum, interspinous, and supraspinous ligaments. However, research characterizing the biomechanics and role of ligament health in spinal pathology and clinical context are scarce. This article provides a comprehensive review of the implications of spinal pathology on the structure, function, and biomechanical properties of the posterior ligamentous complex.

RECENT FINDINGS

Current research characterizing biomechanical properties of the posterior ligamentous complex is primarily composed of cadaveric studies and finite element modeling, and more recently incorporating patient-specific anatomy into finite element models. The ultimate goal of current research is to understand the relative contributions of these ligamentous structures in healthy and pathological spine, and whether preserving ligaments may play an important role in spinal surgical techniques. At baseline, posterior ligamentous complex structures account for 30-40% of spinal stability, which is highly dependent on the intrinsic biomechanical properties of each ligament. Biomechanics vary widely with pathology and following rigid surgical fixation techniques and are generally maladaptive. Often secondary to morphological changes in the setting of spinal pathology, but morphological changes in ligament may also serve as a primary pathology. Biomechanical maladaptations of the spinal ligament adversely influence overall spinal column integrity and ultimately predispose to increased risk for surgical failure and poor clinical outcomes. Future research is needed, particularly in living subjects, to better characterize adaptations in ligaments that can provide targets for improved treatment of spinal pathology.

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.

A comparison of intervertebral ligament properties utilized in a thoracic spine functional unit through kinematic evaluation

Ligament properties in the literature are variable, yet scarce, but needed to calibrate computational models for spine clinical research applications. A comparison of ligament stiffness properties and their effect on the kinematic behavior of a thoracic functional spinal unit (FSU) is examined in this paper. Six unique ligament property sets were utilized within a volumetric T7-T8 finite element (FE) model developed using computer-aided design (CAD) spinal geometry. A 7.5 Nm moment was applied along three anatomical planes both with and without costovertebral (CV) joints present. Range of Motion (RoM) was assessed for each property set and compared to published experimental data. Intact and serial ligament removal procedures were implemented in accordance with experimental protocol. The variance in both kinematic behavior and comparability with experimental data among property sets emphasizes the role nonlinear characterization plays in determining proper kinematic behavior in spinal FE models. Additionally, a decrease in RoM variation among property sets was exhibited when the model setup incorporated the CV joint. With proper assessment of the source and size of each ligament, the material properties considered here could be expanded and justified for implementation into thoracic spine clinical studies.

A neuromuscular human body model for lumbar injury risk analysis in a vibration loading environment

BACKGROUND AND OBJECTIVE

Long-term intensive exposure to whole-body vibration substantially increases the risk of low back pain and degenerative diseases in special occupational groups, like motor vehicle drivers, military vehicle occupants, aircraft pilots, etc. This study aims to establish and validate a neuromuscular human body model focusing on improvement of the detailed description of anatomic structures and neural reflex control, for lumbar injury analysis in vibration loading environments.

METHODS

A whole-body musculoskeletal in Opensim codes was first improved by including a detailed anatomic description of spinal ligaments, non-linear intervertebral disc, and lumbar facet joints, and coupling a proprioceptive feedback closed-loop control strategy with GTOs and muscle spindles modeling in Python codes. Then, the established neuromuscular model was multi-levelly validated from sub-segments to the whole model, from regular movements to dynamic responses to vibration loadings. Finally, the neuromuscular model was combined with a dynamic model of an armored vehicle to analyze occupant lumbar injury risk in vibration loadings due to different road conditions and traveling velocities.

RESULT

Based on a series of biomechanical indexes, including lumbar joint rotation angles, the lumbar intervertebral pressures, the displacement of the lumbar segments, and the lumbar muscle activities, the validation results show that the present neuromuscular model is available and feasible in predicting lumbar biomechanical responses in normal daily movement and vibration loading environments. Furthermore, the combined analysis with the armored vehicle model predicted similar lumbar injury risk to the experimental or epidemiologic studies. The preliminary analysis results also showed that road types and travelling velocities have substantial combined effects on lumbar muscle activities, and indicated that intervertebral joint pressure and muscle activity indexes can need to be jointly considered for lumbar injury risk evaluation.

CONCLUSION

In conclusion, the established neuromuscular model is an effective tool to evaluate vibration loading effects on injury risk of the human body and assist vehicle design vibration comfort by directly concerning the human body injury itself.

Implications of range of motion requirements for the laxity of ligaments in a lumbar finite element model

Finite element models of the lumbar spine often adopt ligament properties from tensile tests without accounting for possible differences between testing and in situ initial ligament length. Such differences could result in laxities or preloads at the beginning of a simulation that would affect the ligament forces, tangent stiffness, and the posture at which they fail. In vivo and in vitro human experimental data reported laxities or preloads. However, laxities or preloads, which could also result from postural differences, are often neglected in simulation studies. This study proposes a numerical methodology to identify ranges of ligament laxities or preloads compatible with the selected tensile ligament properties, the model, and the range of motion (RoM) the model aims to simulate. The approach assumes that ligaments should remain in a safe elongation range for the complete RoM, and that each ligament should play a significant mechanical role in at least one load case. The methodology was applied to the functional spinal unit (FSU) models using the RoM from healthy subjects and ligament properties from the literature. Without laxity, some ligaments reached their elongation at failure within the RoM. Laxity ranges varied considerably (from -9.2 mm preload to 10.7 mm laxity) and flexion was the most critical load case to determine them. Their effect on the mobility response was also assessed. The effect on the mobility of a FSU was also assessed. While the proposed method cannot determine an exact laxity value, it is simple and it can be applied to any model to identify a plausible range of ligament initial length.

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.

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.

The effects of pre-stressed rods contoured by different bending techniques on posterior instrumentation of L4-L5 lumbar spine segment: A finite element study

Posterior pedicle screw instrumentation (PPSI) is a well-known method in lumbar spine surgery. Understanding how PPSI affects the biomechanics of the lumbar spine is an important issue. In particular, during PPSI operations, surgeons bend rods according to the patients’ spinal curvatures based on their own experiences. As a result, residual stresses develop on the rods due to this bending. Although many finite element-based biomechanical studies have been performed for PPSI, studies comparing the effects of residual stresses arising from rod contouring on the construct stresses are lacking. Thus, this study aimed to investigate the effects of residual stress in PPSI using rods contoured with a French bender and an in-situ bender, as well as comparing the differences in stress increment with straight and contoured rods for different physiological motions. Accordingly, a finite element (FE) model of the L4-L5 lumbar spine segment was developed for PPSI and the effects of residual stresses on rods were investigated by using different bending methods. In the simulations, it was found that rods contoured with a French bender with residual stress resulted in significantly more increased stress in PPSI compared to those contoured with an in-situ bender. The results of this study emphasize that increased stress in PPSI due to the residual stresses for different physiological motions may increase the risk of PPSI failures. Additionally, the finite element modeling approach employed here could be used as a fundamental tool for further investigations of topics such as PPSI fatigue life and failure studies.

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.

Effects of Revision Rod Position on Spinal Construct Stability in Lumbar Revision Surgery: A Finite Element Study

Revision surgery (RS) is a necessary surgical intervention in clinical practice to treat spinal instrumentation–related symptomatic complications. Three constructs with different configurations have been applied in RS. One distinguishing characteristic of these configurations is that the revision rods connecting previous segments and revision segments are placed alongside, outside, or inside the previous rods at the level of facetectomy. Whether the position of the revision rod could generate mechanical disparities in revision constructs is unknown. The objective of this study was to assess the influence of the revision rod position on the construct after RS. A validated spinal finite element (FE) model was developed to simulate RS after previous instrumented fusion using a modified dual-rod construct (DRCm), satellite-rod construct (SRC), and cortical bone trajectory construct (CBTC). Thereafter, maximum von Mises stress (VMS) on the annulus fibrosus and cages and the ligament force of the interspinous ligament, supraspinous ligament, and ligamentum flavum under a pure moment load and a follower load in six directions were applied to assess the influence of the revision rod position on the revision construct. An approximately identical overall reducing tendency of VMS was observed among the three constructs. The changing tendency of the maximum VMS on the cages placed at L4-L5 was nearly equal among the three constructs. However, the changing tendency of the maximum VMS on the cage placed at L2-L3 was notable, especially in the CBTC under right bending and left axial rotation. The overall changing tendency of the ligament force in the DRCm, SRC, and CBTC was also approximately equal, while the ligament force of the CBTC was found to be significantly greater than that of the DRCm and SRC at L1-L2. The results indicated that the stiffness associated with the CBTC might be lower than that associated with the DRCm and SRC in RS. The results of the present study indicated that the DRCm, SRC, and CBTC could provide sufficient stabilization in RS. The CBTC was a less rigid construct. Rather than the revision rod position, the method of constructing spinal instrumentation played a role in influencing the biomechanics of revision.

Development of a finite element lumbar spine model to predict intervertebral disc herniation risk factors

The objective of the current study was to develop a lumbar motion segment FE model that predicts disc herniation risk. The posterolateral nucleus extrusion force and disc pressure increased as the amount of flexion and magnitude of compression was increased in all loading scenarios. The nucleus extrusion force and posterior stress in the annulus both increased when exposed to a combination of compression and flexion. Results of the current study confirmed the authors hypothesis that the model would accurately predict herniation risk when exposed to a biomechanical environment known to cause herniations.

Computational modeling of lumbar disc degeneration before and after spinal fusion

BACKGROUND

Advancing age and degeneration frequently lead to low back pain, which is the most prevalent musculoskeletal disorder worldwide. Degenerative changes in intervertebral discs and musculo-ligamentous incapacity to compensate sagittal imbalance are typically amongst the sources of instability, with spinal fusion techniques being the main treatment options to relieve pain. The aims of this work were to: (i) assess the link between ligament degeneration and spinal instability by determining the role of each ligament per movement, (ii) evaluate the impact of disc height reduction in degenerative changes, and (iii) unveil the most advantageous type of posterior fixation in Oblique Lumbar Interbody Fusion to prevent adjacent disc degeneration.

METHODS

Two L3-L5 finite element models were developed, being the first in healthy condition and the second having reduced L4-L5 height. Different degrees of degeneration were tested, combined with different fixation configurations for Oblique Lumbar Interbody Fusion.

FINDINGS

Facet capsular ligament and anterior longitudinal ligament were the most influential ligaments for spinal stability, particularly with increasing degeneration and disc height reduction. Pre-existent degeneration had lower influence than the fusion procedure for the risk of adjacent disc degeneration, being the highest stability and minimal degeneration achieved with bilateral fixation. Right unilateral fixation was more suited to reduce disc stress than left unilateral fixation.

INTERPRETATION

Bilateral fixation is the best option to stabilize the spinal segment, but unilateral right fixation may suffice. This has direct implications for clinical practice, and the extension to a population-based study will allow for more efficient fusion surgeries.

Influence of posterior pedicle screw fixation at L4-L5 level on biomechanics of the lumbar spine with and without fusion: a finite element method

Background

Posterior pedicle screw (PS) fixation, a common treatment method for widespread low-back pain problems, has many uncertain aspects including stress concentration levels, effects on adjacent segments, and relationships with physiological motions. A better understanding of how posterior PS fixation affects the biomechanics of the lumbar spine is needed. For this purpose, a finite element (FE) model of a lumbar spine with posterior PS fixation at the L4–L5 segment level was developed by partially removing facet joints (FJs) to imitate an actual surgical procedure. This FE study aimed to investigate the influence of the posterior PS fixation system on the biomechanics of the lumbar spine before and after fusion by determining which physiological motions have the most increase in posterior instrumentation (PI) stresses and FJ loading.

Results

It was determined that posterior PS fixation increased FJ loading by approximately 35% and 23% at the L3–L4 adjacent level with extension and lateral bending motion, respectively. This increase in FJ loading at the adjacent level could point to the possibility that adjacent segment disease has developed or progressed after posterior lumbar interbody fusion. Furthermore, analyses of peak von Mises stresses on PI showed that the maximum PI stresses of 272.1 MPa and 263.7 MPa occurred in lateral bending and flexion motion before fusion, respectively.

Conclusions

The effects of a posterior PS fixation system on the biomechanics of the lumbar spine before and after fusion were investigated for all physiological motions. This model could be used as a fundamental tool for further studies, providing a better understanding of the effects of posterior PS fixation by clearing up uncertain aspects.

Computational modelling of the scoliotic spine: A literature review

Scoliosis is a deformity of the spine that in severe cases requires surgical treatment. There is still disagreement among clinicians as to what the aim of such treatment is as well as the optimal surgical technique. Numerical models can aid clinical decision-making by estimating the outcome of a given surgical intervention. This paper provided some background information on the modelling of the healthy spine and a review of the literature on scoliotic spine models, their validation, and their application. An overview of the methods and techniques used to construct scoliotic finite element and multibody models was given as well as the boundary conditions used in the simulations. The current limitations of the models were discussed as well as how such limitations are addressed in non-scoliotic spine models. Finally, future directions for the numerical modelling of scoliosis were addressed.

An Electromyographically Driven Cervical Spine Model in OpenSim

Relatively few biomechanical models exist aimed at quantifying the mechanical risk factors associated with neck pain. In addition, there is a need to validate spinal-rhythm techniques for inverse dynamics spine models. Therefore, the present investigation was 3-fold: (1) the development of a cervical spine model in OpenSim, (2) a test of a novel spinal-rhythm technique based on minimizing the potential energy in the passive tissues, and (3) comparison of an electromyographically driven approach to estimating compression and shear to other cervical spine models. The authors developed ligament force-deflection and intervertebral joint moment-angle curves from published data. The 218 Hill-type muscle elements, representing 58 muscles, were included and their passive forces validated against in vivo data. Our novel spinal-rhythm technique, based on minimizing the potential energy in the passive tissues, disproportionately assigned motion to the upper cervical spine that was not physiological. Finally, using kinematics and electromyography collected from 8 healthy male volunteers, the authors calculated the compression at C7-T1 as a function of the head-trunk Euler angles. Differences from other models varied from 25.5 to 368.1 N. These differences in forces may result in differences in model geometry, passive components, number of degrees of freedom, or objective functions.

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).

Development of a flexible instrumented lumbar spine finite element model and comparison with in-vitro experiments

Surgical corrections of degenerative lumbar scoliosis and sagittal malalignment are associated with significant complications, such as rod fractures and pseudarthrosis, particularly in the lumbosacral junction. Finite element studies can provide relevant insights to improve performance of spinal implants. The aim of the present study was to present the development of non-instrumented and instrumented Finite Element Models (FEMs) of the lumbopelvic spine and to compare numerical results with experimental data available in the literature. The lumbo-pelvic spine FEM was based on a CT-scan from an asymptomatic volunteer representing the 50th percentile male. In a first step a calibration of mechanical properties was performed in order to obtain a quantitative agreement between numerical results and experimental data for defect stages of spinal segments. Then, FEM results were compared in terms of range of motion and strains in rods to in-vitro experimental data from the literature for flexible non-instrumented and instrumented lumbar spines. Numerical results from the calibration process were consistent with experimental data, especially in flexion. A positive agreement was obtained between FEM and experimental results for the lumbar and sacroiliac segments. Instrumented FEMs predicted the same trends as experimental in-vitro studies. The instrumentation configuration consisting of double rods and an interbody cage at L5-S1 maximally reduced range of motion and strains in main rods and thus had the lowest risk of pseudarthrosis and rod fracture. The developed FEMs were found to be consistent with published experimental results; therefore they can be used for further post-operative complication investigations.

Construction of Vertebral Body Tracking Algorithm Based on Dynamic Imaging Parameter Measurement and Its Application in the Treatment of Lumbar Instability

Static imaging measurements could not truly reflect the dynamic panorama of the lumbar movement process, and the abnormal activities between the lumbar vertebrae and their dynamic balance could not be observed, resulting in difficulties in the mechanism analysis of lumbar instability and the efficacy evaluation of manipulation therapy. Therefore, this paper constructed a vertebral tracking algorithm based on dynamic imaging parameter measurement through imaging parameter measurement and calculation. According to the imaging data obtained by vertebral body tracking algorithm, the corresponding statistical methods were used to compare the functional scores before and after manipulation and the changes of imaging data, so as to evaluate the therapeutic effect of manipulation on lumbar instability. Through the clinical observation and imaging analysis of 15 patients with lumbar instability before and after manipulation treatment, it is verified that the vertebra tracking algorithm is effective in the vertebra tracking and plays a positive role in the treatment of lumbar instability.

Distraction forces on the spine in early-onset scoliosis: A systematic review and meta-analysis of clinical and biomechanical literature

Distraction-based growing rods are frequently used to treat Early-Onset Scoliosis. These use intermittent spinal distractions to maintain correction and allow for growth. It is unknown how much spinal distraction can be applied safely. We performed a systematic review and meta-analysis of clinical and biomechanical literature to identify such safety limits for the pediatric spine. This systematic review and meta-analysis was performed according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement. Three systematic searches were performed including in-vivo, ex-vivo and in-silico literature. Study quality was assessed in all studies and data including patient- or specimen characteristics, distraction magnitude and spinal failure location and ultimate force at failure were collected. Twelve studies were included, 6 in-vivo, 4 ex-vivo and 2 in-silico studies. Mean in-vivo distraction forces ranged between 242 and 621 N with maxima of 422-981 N, without structural failures when using pedicle screw constructs. In the ex-vivo studies (only cervical spines), segment C0-C2 was strongest, with decreasing strength in more distal segments. Meta-regression analysis demonstrated that ultimate force at birth is 300-350 N, which increases approximately 100 N each year until adulthood. Ex-vivo and in-silico studies showed that yielding occurs at 70-90% of ultimate force, failure starts at the junction between endplate and intervertebral disc, after which the posterior- and anterior long ligament rupture. While data on safety of distraction forces is limited, this systematic review and meta-analysis may aid in the development of guidelines on spinal distraction and may benefit the development and optimization of contemporary and future distraction-based technologies.

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.

A comprehensive approach for the validation of lumbar spine finite element models investigating post-fusion adjacent segment effects

Spinal fusion surgery is usually followed by accelerated degenerative changes in the unfused segments above and below the treated segment(s), i.e., adjacent segment disease (ASD). While a number of risk factors for ASD have been suggested, its exact pathogenesis remains to be identified. Finite element (FE) models are indispensable tools to investigate mechanical effects of fusion surgeries on post-fusion changes in the adjacent segment kinematics and kinetics. Existing modeling studies validate only their intact FE model against in vitro data and subsequently simulate post-fusion in vivo conditions. The present study provides a novel approach for the comprehensive validation of a lumbar (T12-S1) FE model in post-fusion conditions. Sixteen simulated fusion surgeries, performed on cadaveric specimens using various testing and loading conditions, were modeled by this FE model. Predictions for adjacent segment range of motion (RoM) and intradiscal pressure (IDP) were compared with those obtained from the corresponding in vitro tests. Overall, 70% of the predicted adjacent segment RoMs were within the range of in vitro data for both intact and post-fusion conditions. Correlation (r) values between model and in vitro findings for the adjacent segment RoMs were positive and greater than 0.84. Most of the predicted IDPs were, however, out of the narrow range of in vitro IDPs at the adjacent segments but with great positive correlations (r ≥ 0.89). FE modeling studies investigating the effect of fusion surgery on in vivo adjacent segment biomechanics are encouraged to use post-surgery in vitro data to validate their FE model.

Finite element analysis of lumbar spine with different backpack positions in parachuting landing

The purpose of this study was to investigate the lumbar spine stress with different backpack positions in parachuting landing using a finite element model of lumbar vertebra 1-5. The backpack gravity center was set at three positions (posterior-high (case PH), posterior-low (case PL), and anterior-low (case AL)) respectively. In results, the peak Von-Mises stresses of the matrix, nucleus, fibers, endplate and ligament in case AL were 2.765 MPa, 0.534 MPa, 6.561 MPa, 4.045 MPa and 1.790 MPa respectively, lower than those in case PL (6.913 MPa, 1.316 MPa, 20.716 MPa, 10.917 MPa and 5.147 MPa respectively) and case PH (7.328 MPa, 1.394 MPa, 22.147 MPa, 11.617 MPa and 5.464 MPa respectively). In conclusion, setting the gravity center of backpack at anterior-low position would reduce lumbar spine stress and reduce lumbar spine injuries.

The importance of intervertebral disc material model on the prediction of mechanical function of the cervical spine

BACKGROUND

Linear elastic, hyperelastic, and multiphasic material constitutive models are frequently used for spinal intervertebral disc simulations. While the characteristics of each model are known, their effect on spine mechanical response requires a careful investigation. The use of advanced material models may not be applicable when material constants are not available, model convergence is unlikely, and computational time is a concern. On the other hand, poor estimations of tissue's mechanical response are likely if the spine model is oversimplified. In this study, discrepancies in load response introduced by material models will be investigated.

METHODS

Three fiber-reinforced C2-C3 disc models were developed with linear elastic, hyperelastic, and biphasic behaviors. Three different loading modes were investigated: compression, flexion and extension in quasi-static and dynamic conditions. The deformed disc height, disc fluid pressure, range of motion, and stresses were compared.

RESULTS

Results indicated that the intervertebral disc material model has a strong effect on load-sharing and disc height change when compression and flexion were applied. The predicted mechanical response of three models under extension had less discrepancy than its counterparts under flexion and compression. The fluid-solid interaction showed more relevance in dynamic than quasi-static loading conditions. The fiber-reinforced linear elastic and hyperelastic material models underestimated the load-sharing of the intervertebral disc annular collagen fibers.

CONCLUSION

This study confirmed the central role of the disc fluid pressure in spinal load-sharing and highlighted loading conditions where linear elastic and hyperelastic models predicted energy distribution different than that of the biphasic model.

Patterns of Load-to-Strength Ratios Along the Spine in a Population-Based Cohort to Evaluate the Contribution of Spinal Loading to Vertebral Fractures

Vertebral fractures (VFx) are common among older adults. Epidemiological studies report high occurrence of VFx at mid-thoracic and thoracolumbar regions of the spine; however, reasons for this observation remain poorly understood. Prior reports of high ratios of spinal loading to vertebral strength in the thoracolumbar region suggest a possible biomechanical explanation. However, no studies have evaluated load-to-strength ratios (LSRs) throughout the spine for a large number of activities in a sizeable cohort. Thus, we performed a cross-sectional study in a sample of adult men and women from a population-based cohort to: 1) determine which activities cause the largest vertebral LSRs, and 2) examine patterns of LSRs along the spine for these high-load activities. We used subject-specific musculoskeletal models of the trunk to determine vertebral compressive loads for 109 activities in 250 individuals (aged 41 to 90 years, 50% women) from the Framingham Heart Study. Vertebral compressive strengths from T₄ to L₄ were calculated from computed tomography-based vertebral size and bone density measurements. We determined which activities caused maximum LSRs at each of these spinal levels. We identified nine activities that accounted for >95% of the maximum LSRs overall and at least 89.6% at each spinal level. The activity with the highest LSR varied by spinal level, and three distinct spinal regions could be identified by the activity producing maximum LSRs: lateral bending with a weight in one hand (upper thoracic), holding weights with elbows flexed (lower thoracic), and forward flexion with weight (lumbar). This study highlights the need to consider a range of lifting, holding, and non-symmetric activities when evaluating vertebral LSRs. Moreover, we identified key activities that produce higher loading in multiple regions of the spine. These results provide the first guidance on what activities to consider when evaluating vertebral load-to-strength ratios in future studies, including those examining dynamic motions and the biomechanics of VFx. © 2020 American Society for Bone and Mineral Research (ASBMR).

Computational Challenges in Tissue Engineering for the Spine

This paper deals with a brief review of the recent developments in computational modelling applied to innovative treatments of spine diseases. Additionally, it provides a perspective on the research directions expected for the forthcoming years. The spine is composed of distinct and complex tissues that require specific modelling approaches. With the advent of additive manufacturing and increasing computational power, patient-specific treatments have moved from being a research trend to a reality in clinical practice, but there are many issues to be addressed before such approaches become universal. Here, it is identified that the major setback resides in validation of these computational techniques prior to approval by regulatory agencies. Nevertheless, there are very promising indicators in terms of optimised scaffold modelling for both disc arthroplasty and vertebroplasty, powered by a decisive contribution from imaging methods.

The importance of curve severity, type and instrumentation strategy in the surgical correction of adolescent idiopathic scoliosis: an in silico clinical trial on 64 cases

Adolescent idiopathic scoliosis is a three-dimensional deformity of the spine which is frequently corrected with the implantation of instrumentation with generally good or excellent clinical results; mechanical post-operative complications such as implant loosening and breakage are however relatively frequent. The rate of complications is associated with a lack of consensus about the surgical decision-making process; choices about the instrumentation length, the anchoring implants and the degree of correction are indeed mostly based on personal views and previous experience of the surgeon. In this work, we performed an in silico clinical trial on a large number of subjects in order to clarify which factors have the highest importance in determining the risk of complications by quantitatively analysing the mechanical stresses and loads in the instrumentation after the correction maneuvers. The results of the simulations highlighted the fundamental role of the curve severity, also in its three-dimensional aspect, and of the instrumentation strategy, whereas the length of the fixation had a lower importance.

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.

Effect of changes in the lumbar posture in lifting on trunk muscle and spinal loads: A combined in vivo, musculoskeletal, and finite element model study

Irrespective of the lifting technique (squat or stoop), the lumbar spine posture (more kyphotic versus more lordotic) adopted during lifting activities is an important parameter affecting the active-passive spinal load distribution. The advantages in either posture while lifting remains, however, a matter of debate. To comprehensively investigate the role on the trunk biomechanics of changes in the lumbar posture (lordotic, free or kyphotic) during forward trunk flexion, validated musculoskeletal and finite element models, driven by in vivo kinematics data, were used to estimate detailed internal tissue stresses-forces in and load-sharing among various joint active-passive tissues. Findings indicated that the lordotic posture, as compared to the kyphotic one, resulted in marked increases in back global muscle activities (~14-19%), overall segmental compression (~7.5-46.1%) and shear (~5.4-47.5%) forces, and L5-S1 facet joint forces (by up to 80 N). At the L5-S1 level, the lordotic lumbar posture caused considerable decreases in the moment resisted by passive structures (spine and musculature, ~14-27%), negligible reductions in the maximum disc fiber strains (by ~0.4-4.7%) and small increases in intradiscal pressure (~1.8-3.4%). Collectively and with due consideration of the risk of fatigue and viscoelastic creep especially under repetitive lifts, current results support a free posture (in between the extreme kyphotic and lordotic postures) with moderate contributions from both active and passive structures during lifting activities involving trunk forward flexion.

Lumbar spinal ligament characteristics extracted from stepwise reduction experiments allow for preciser modeling than literature data

Lumbar ligaments play a key role in stabilizing the spine, particularly assisting muscles at wide-range movements. Hence, valid ligament force–strain data are required to generate physiological model predictions. These data have been obtained by experiments on single ligaments or functional units throughout the literature. However, contrary to detailed spine geometries, gained, for instance, from CT data, ligament characteristics are often inattentively transferred to multi-body system (MBS) or finite element models. In this paper, we use an elaborated MBS model of the lumbar spine to demonstrate how individualized ligament characteristics can be obtained by reversely reenacting stepwise reduction experiments, where the range of motion (ROM) was measured. We additionally validated the extracted characteristics with physiological experiments on intradiscal pressure (IDP). Our results on a total of in each case 160 ROM and 49 IDP simulations indicated superiority of our procedure (seven and eight outliers) toward the incorporation of classical literature data (on average 71 and 31 outliers).

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 role of the facet capsular ligament in providing spinal stability

Low back pain (LBP) is the most common type of pain in America, and spinal instability is a primary cause. The facet capsular ligament (FCL) encloses the articulating joints of the spine and is of particular interest due to its high innervation - as instability ensues, high stretch values likely are a cause of this pain. Therefore, this work investigated the FCL's role in providing stability to the lumbar spine. A previously validated finite element model of the L4-L5 spinal motion segment was used to simulate pure moment bending in multiple planes. FCL failure was simulated and the following outcome measures were calculated: helical axes of motion, range of motion (ROM), bending stiffness, facet joint space, and FCL stretch. ROM increased, bending stiffness decreased, and altered helical axis patterns were observed with the removal of the FCL. Additionally, a large increase in FCL stretch was measured with diminished FCL mechanical competency, providing support that the FCL plays an important role in spinal stability.