Optimised in vitro applicable loads for the simulation of lateral bending in the lumbar spine

Load Changes on a Short-Segment Posterior Instrumentation After Transosseous Disruption of L3 Vertebra - A Biomechanical Human Cadaveric Study

STUDY DESIGN

Biomechanical Cadaveric Study.

OBJECTIVES

Following the successful use of a novel implantable sensor (Monitor) in evaluating the progression of fracture healing in long bones and posterolateral fusion of the spine based on implant load monitoring, the aim of this study was to investigate its potential to assess healing of transosseous fractures of a lumbar vertebra stabilized with a pedicle-screw-rod construct.

METHODS

Six human cadaveric spines were instrumented with pedicle screws and rods spanning L3 vertebra. The spine was loaded in Flexion-Extension (FE), Lateral-Bending (LB) and Axial-Rotation (AR) with an intact L3 vertebra and after its transosseous disruption, creating an AO B1 type fracture. The implant load was measured on the one rod using the Monitor and on the contralateral rod by strain gauges to validate the Monitor's measurements. In parallel, the range of motion (ROM) was assessed.

RESULTS

ROM increased significantly in all directions in the fractured model (P ≤ 0.049). The Monitor measured a significant increase in implant load in FE (P = 0.002) and LB (P = 0.045), however, not in AR. The strain gauge - aligned with the rod axis and glued onto its posterior side - detected an increased implant load not only in FE (P = 0.001) and LB (P = 0.016) but also in AR (P = 0.047).

CONCLUSION

After a complete transosseous disruption of L3 vertebra, the implant load on the rods was considerably higher vs the state with an intact vertebral body. Innovative implantable sensors could monitor those changes, allowing assessment of the healing progression based on quantifiable data.

The effect of large channel-based foraminoplasty on lumbar biomechanics in percutaneous endoscopic discectomy: a finite element analysis

BACKGROUND

This study aimed to evaluate the effect of foraminoplasty using large-channel endoscopy during TESSYS on the biomechanics of the lumbar spine.

METHODS

A complete lumbar spine model, M1, was built using 3D finite elements, and models M2 and M3 were constructed to simulate the intraoperative removal of the superior articular process of L5 using a trephine saw with diameters of 5 mm and 8.5 mm, respectively, and applying normal physiological loads on the different models to simulate six working conditions-anterior flexion, posterior extension, left-right lateral bending, and left-right rotation-to investigate the displacement and facet joint stress change of the surgical segment, and the disc stress change of the surgical and adjacent segments.

RESULTS

Compared with the M1 model, the M2 and M3 models showed decreased stress at the L4-5 left FJ and a significant increase in stress at the right FJ in forward flexion. In the M2 and M3 models, the L4-5 FJ stresses were significantly greater in left lateral bending or left rotation than in right lateral bending or right rotation. The right FJ stress in M3 was greater during left rotation than that in M2, and that in M2 was greater than that in M1. The L4-5disc stress in the M3 model was greater during posterior extension than that in the M1 and M2 models. The L4-5disc stress in the M3 model was greater in the right rotation than in the M2 model, and that in the M2 model was greater than that in the M1 model.

CONCLUSION

Foraminoplasty using large-channel endoscopy could increase the stress on the FJ and disc of the surgical segment, which suggested unnecessary and excessive resection should be avoided in PTED to minimize biomechanical disruption.

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

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

Anatomical parameters alter the biomechanical responses of adjacent segments following lumbar fusion surgery: Personalized poroelastic finite element modelling investigations

Introduction: While the short-term post-operative outcome of lumbar fusion is satisfying for most patients, adjacent segment disease (ASD) can be prevalent in long-term clinical observations. It might be valuable to investigate if inherent geometrical differences among patients can significantly alter the biomechanics of adjacent levels post-surgery. This study aimed to utilize a validated geometrically personalized poroelastic finite element (FE) modeling technique to evaluate the alteration of biomechanical response in adjacent segments post-fusion. Methods: Thirty patients were categorized for evaluation in this study into two distinct groups [i.e., 1) non-ASD and 2) ASD patients] based on other long-term clinical follow-up investigations. To evaluate the time-dependent responses of the models subjected to cyclic loading, a daily cyclic loading scenario was applied to the FE models. Different rotational movements in different planes were superimposed using a 10 Nm moment after daily loading to compare the rotational motions with those at the beginning of cyclic loading. The biomechanical responses of the lumbosacral FE spine models in both groups were analyzed and compared before and after daily loading. Results: The achieved comparative errors between the FE results and clinical images were on average below 20% and 25% for pre-op and post-op models, respectively, which confirms the applicability of this predictive algorithm for rough pre-planning estimations. The results showed that the disc height loss and fluid loss were increased for the adjacent discs in post-op models after 16 h of cyclic loading. In addition, significant differences in disc height loss and fluid loss were observed between the patients who were in the non-ASD and ASD groups. Similarly, the increased stress and fiber strain in the annulus fibrosus (AF) was higher in the adjacent level of post-op models. However, the calculated stress and fiber strain values were significantly higher for patients with ASD. Discussion: Evaluating the biomechanical response of pre-op and post-op modeling in the non-ASD and ASD groups showed that the inherent geometric differences among patients cause significant variations in the estimated mechanical response. In conclusion, the results of the current study highlighted the effect of geometrical parameters (which may refer to the anatomical conditions or the induced modifications regarding surgical techniques) on time-dependent responses of lumbar spine biomechanics.

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.

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

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

Biphasic Properties of PVAH (Polyvinyl Alcohol Hydrogel) Reflecting Biomechanical Behavior of the Nucleus Pulposus of the Human Intervertebral Disc

PVAH is a mixture of solid and fluid, but its mechanical behavior has usually been described using solid material models. The purpose of this study was to obtain material properties that can reflect the mechanical behavior of polyvinyl alcohol hydrogel (PVAH) using finite element analysis, a biphasic continuum model, and to optimize the composition ratio of PVAH to replace the nucleus pulposus (NP) of the human intervertebral disc. Six types of PVAH specimens (3, 5, 7, 10, 15, 20 wt%) were prepared, then unconfined compression experiments were performed to acquire their material properties using the Holmes–Mow biphasic model. With an increasing weight percentage of PVA in PVAH, the Young’s modulus increased while the permeability parameter decreased. The Young’s modulus and permeability parameter were similar to those of the NP at 15 wt% and 20 wt%. The range of motion, facet joint force, and NP pressures measured from dynamic motional analysis of the lumbar segments with the NP model also exhibited similar values to those with 15~20 wt% PVAH models. Considering the structural stability and pain of the lumbar segments, it appears that 20 wt% PVAH is most suitable for replacing the NP.

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.

Stress analysis of intervertebral disc during occupational activities

BACKGROUND AND OBJECTIVE

Manual material handling activities cause large compression of the intervertebral disc of the lumbar spine. Intradiscal pressure (IDP) has generally been employed to predict the risk of low back injury. As an alternative to in vivo measurements, either motion analysis or finite element (FE) analysis has been used to estimate IDP. The purpose of this study is to propose a new biomechanical method that integrates FE analysis with motion analysis, in order to estimate the stresses and deformations of the intervertebral disc of the lumbar spine during occupational activities.

METHODS

In the proposed method, motion analysis is performed first by using motion capture data, and the results are employed as input data to FE analysis at specific times of interest during motion. In this method, an in-house interface program is used to scale an initial reference FE model to the subject of experiment, and transformed to the corresponding posture at a specific time during motion. The muscle forces and GRF obtained from motion analysis are applied to FE analysis as boundary and loading conditions. For a total of eighteen occupational activities, the IDP, shear stress, and strain of the L4-L5 segment are estimated.

RESULTS

Under each in vivo activity, the predicted IDP was in overall agreement with the available in vivo data. For lifting activities according to lift origin position, the maximum IDP occurred in the far-knee position immediately after lifting. As the lift origin position moved away from the spine, the stresses and strains in the disc increased.

CONCLUSIONS

This new proposed method is expected to allow the estimation of the stresses and deformations in the intervertebral disc during various occupational activities.

Numerical Evaluation of Spinal Stability after Posterior Spinal Fusion with Various Fixation Segments and Screw Types in Patients with Osteoporotic Thoracolumbar Burst Fracture Using Finite Element Analysis

The aim of this study was to analyze the spinal stability and safety after posterior spinal fusion with various fixation segments and screw types in patients with an osteoporotic thoracolumbar burst fracture based on finite element analysis (FEA). To realize various osteoporotic vertebral fracture conditions on T12, seven cases of Young’s modulus, namely 0%, 1%, 5%, 10%, 25%, 50%, and 100% of the Young’s modulus, for vertebral bones under intact conditions were considered. Four types of fixation for thoracolumbar fracture on T12 (fixed with T11-L1, T10-T11-L1, T11-L1-L2, and T10-T11-L1-L2) were applied to the thoracolumbar fusion model. The following screw types were considered: pedicle screw (PS) and cortical screw (CS). Using FEA, four motions were performed on the fixed spine, and the stress applied to the screw, peri-implant bone (PIB), and intervertebral disc (IVD) and the range of motion (ROM) were calculated. The lowest ROM calculated corresponded to the T10-T11-L1-L2 model, while the closest to the intact situation was achieved in the T11-L1-L2 fixation model using PS. The lowest stress in the screw and PB was detected in the T10-T11-L1-L2 fixation model.

Biomechanical Investigation Between Rigid and Semirigid Posterolateral Fixation During Daily Activities: Geometrically Parametric Poroelastic Finite Element Analyses

While spinal fusion using rigid rods remains the gold standard treatment modality for various lumbar degenerative conditions, its adverse effects, including accelerated adjacent segment disease (ASD), are well known. In order to better understand the performance of semirigid constructs using polyetheretherketone (PEEK) in fixation surgeries, the objective of this study was to analyze the biomechanical performance of PEEK versus Ti rods using a geometrically patient-specific poroelastic finite element (FE) analyses. Ten subject-specific preoperative models were developed, and the validity of the models was evaluated with previous studies. Furthermore, FE models of those lumbar spines were regenerated based on postoperation images for posterolateral fixation at the L4–L5 level. Biomechanical responses for instrumented and adjacent intervertebral discs (IVDs) were analyzed and compared subjected to static and cyclic loading. The preoperative model results were well comparable with previous FE studies. The PEEK construct demonstrated a slightly increased range of motion (ROM) at the instrumented level, but decreased ROM at adjacent levels, as compared with the Ti. However, no significant changes were detected during axial rotation. During cyclic loading, disc height loss, fluid loss, axial stress, and collagen fiber strain in the adjacent IVDs were higher for the Ti construct when compared with the intact and PEEK models. Increased ROM, experienced stress in AF, and fiber strain at adjacent levels were observed for the Ti rod group compared with the intact and PEEK rod group, which can indicate the risk of ASD for rigid fixation. Similar to the aforementioned pattern, disc height loss and fluid loss were significantly higher at adjacent levels in the Ti rod group after cycling loading which alter the fluid–solid interaction of the adjacent IVDs. This phenomenon debilitates the damping quality, which results in disc disability in absorbing stress. Such finding may suggest the advantage of using a semirigid fixation system to decrease the chance of ASD.

Spine biomechanical testing methodologies: The controversy of consensus vs scientific evidence

Biomechanical testing methodologies for the spine have developed over the past 50 years. During that time, there have been several paradigm shifts with respect to techniques. These techniques evolved by incorporating state-of-the-art engineering principles, in vivo measurements, anatomical structure-function relationships, and the scientific method. Multiple parametric studies have focused on the effects that the experimental technique has on outcomes. As a result, testing methodologies have evolved, but there are no standard testing protocols, which makes the comparison of findings between experiments difficult and conclusions about in vivo performance challenging. In 2019, the international spine research community was surveyed to determine the consensus on spine biomechanical testing and if the consensus opinion was consistent with the scientific evidence. More than 80 responses to the survey were received. The findings of this survey confirmed that while some methods have been commonly adopted, not all are consistent with the scientific evidence. This review summarizes the scientific literature, the current consensus, and the authors' recommendations on best practices based on the compendium of available evidence.

Development of a novel geometrically-parametric patient-specific finite element model to investigate the effects of the lumbar lordosis angle on fusion surgery

The success of lumbar interbody fusion, the key surgical procedure for treating different pathologies of the lumbar spine, is highly dependent on determining the patient-specific lumbar lordosis (LL) and restoring sagittal balance. This study aimed to (1) develop a personalized finite element (FE) model that automatically updates spinal geometry for different patients; and (2) apply this technique to study the influence of LL on post-fusion spinal biomechanics. Using an X-Ray image-based algorithm, the geometry of the lumbar spine (L1-S1) was updated using independent parameters. Ten subject-specific nonlinear osteoligamentous FE models were developed based on pre-operative images of fusion surgery candidate patients. Post-operative FE models of the same patients were consequently created. Comparison of the obtained results from FE models with pre- and post-operation functional images demonstrated the potential value of this technique in clinical applications. A parametric study of the effect of LL was conducted for cases with zero LL angle, positive LL angles (+6° and +12°) and negative LL angles (-3° and -6°) on fused level (L4-L5), resulting in a total of 50 fusion simulation models. The average range of motion, intradiscal pressure, and fiber strain at adjacent levels were significantly higher with decreased LL during different directions except axial rotation. This study demonstrates that the LL alters both the intersegmental motion and load-sharing in fusion, which may influence the initiation and rate of adjacent level degeneration. This personalized FE platform provides a practical, clinically applicable approach for the analyses of the biomechanical changes associated with lumbar spine fusion.

Optimization of compressive loading parameters to mimic in vivo cervical spine kinematics in vitro

The human cervical spine supports substantial compressive load in vivo. However, the traditional in vitro testing methods rarely include compressive loads, especially in investigations of multi-segment cervical spine constructs. Previously, a systematic comparison was performed between the standard pure moment with no compressive loading and published compressive loading techniques (follower load - FL, axial load - AL, and combined load - CL). The systematic comparison was structured a priori using a statistical design of experiments and the desirability function approach, which was chosen based on the goal of determining the optimal compressive loading parameters necessary to mimic the segmental contribution patterns exhibited in vivo. The optimized set of compressive loading parameters resulted in in vitro segmental rotations that were within one standard deviation and 10% of average percent error of the in vivo mean throughout the entire motion path. As hypothesized, the values for the optimized independent variables of FL and AL varied dynamically throughout the motion path. FL was not necessary at the extremes of the flexion-extension (FE) motion path but peaked through the neutral position, whereas, a large negative value of AL was necessary in extension and increased linearly to a large positive value in flexion. Although further validation is required, the long-term goal is to develop a "physiologic" in vitro testing method, which will be valuable for evaluating adjacent segment effect following spinal fusion surgery, disc arthroplasty instrumentation testing and design, as well as mechanobiology experiments where correct kinematics and arthrokinematics are critical.

Influence of follower load application on moment-rotation parameters and intradiscal pressure in the cervical spine

The objective of this study was to implement a follower load (FL) device within a robotic (universal force-moment sensor) testing system and utilize the system to explore the effect of FL on multi-segment cervical spine moment-rotation parameters and intradiscal pressure (IDP) at C45 and C56. Twelve fresh-frozen human cervical specimens (C3-C7) were biomechanically tested in a robotic testing system to a pure moment target of 2.0 Nm for flexion and extension (FE) with no compression and with 100 N of FL. Application of FL was accomplished by loading the specimens with bilateral cables passing through cable guides inserted into the vertebral bodies and attached to load controlled linear actuators. FL significantly increased neutral zone (NZ) stiffness and NZ width but resulted in no change in the range of motion (ROM) or elastic zone stiffness. C45 and C56 IDP measured in the neutral position were significantly increased with application of FL. The change in IDP with increasing flexion rotation was not significantly affected by the application of FL, whereas the change in IDP with increasing extension rotation was significantly reduced by the application of FL. Application of FL did not appear to affect the specimen's quantity of motion (ROM) but did affect the quality (the shape of the curve). Regarding IDP, the effects of adding FL compression approximates the effect of the patient going from supine to a seated position (FL compression increased the IDP in the neutral position). The change in IDP with increasing flexion rotation was not affected by the application of FL, but the change in IDP with increasing extension rotation was, however, significantly reduced by the application of FL.

Effect of Graded Facetectomy on Lumbar Biomechanics

Facetectomy is an important intervention for spinal stenosis but may lead to spinal instability. Biomechanical knowledge for facetectomy can be beneficial when deciding whether fusion is necessary. Therefore, the aim of this study was to investigate the biomechanical effect of different grades of facetectomy. A three-dimensional nonlinear finite element model of L3–L5 was constructed. The mobility of the model and the intradiscal pressure (IDP) of L4-L5 for standing were inside the data from the literature. The effect of graded facetectomy on intervertebral rotation, IDP, facet joint forces, and maximum von Mises equivalent stresses in the annuli was analyzed under flexion, extension, left/right lateral bending, and left/right axial rotation. Compared with the intact model, under extension, unilateral facetectomy increased the range of intervertebral rotation (IVR) by 11.7% and IDP by 10.7%, while the bilateral facetectomy increased IVR by 40.7% and IDP by 23.6%. Under axial rotation, the unilateral facetectomy and the bilateral facetectomy increased the IVR by 101.3% and 354.3%, respectively, when turned to the right and by 1.1% and 265.3%, respectively, when turned to the left. The results conclude that, after unilateral and bilateral facetectomy, care must be taken when placing the spine into extension and axial rotation posture from the biomechanical point of view.

Sensitivity of lumbar spine response to follower load and flexion moment: finite element study

The follower load (FL) combined with moments is commonly used to approximate flexed/extended posture of the lumbar spine in absence of muscles in biomechanical studies. There is a lack of consensus as to what magnitudes simulate better the physiological conditions. Considering the in-vivo measured values of the intradiscal pressure (IDP), intervertebral rotations (IVRs) and the disc loads, sensitivity of these spinal responses to different FL and flexion moment magnitudes was investigated using a 3D nonlinear finite element (FE) model of ligamentous lumbosacral spine. Optimal magnitudes of FL and moment that minimize deviation of the model predictions from in-vivo data were determined. Results revealed that the spinal parameters i.e. the IVRs, disc moment, and the increase in disc force and moment from neutral to flexed posture were more sensitive to moment magnitude than FL magnitude in case of flexion. The disc force and IDP were more sensitive to the FL magnitude than moment magnitude. The optimal ranges of FL and flexion moment magnitudes were 900-1100 N and 9.9-11.2 Nm, respectively. The FL magnitude had reverse effect on the IDP and disc force. Thus, magnitude for FL or flexion that minimizes the deviation of all the spinal parameters together from the in-vivo data can vary. To obtain reasonable compromise between the IDP and disc force, our findings recommend that FL of low magnitude must be combined with flexion moment of high intensity and vice versa.

Comparison of eight published static finite element models of the intact lumbar spine: predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centers around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.

Parameters influencing the outcome after total disc replacement at the lumbosacral junction. Part 1: misalignment of the vertebrae adjacent to a total disc replacement affects the facet joint and facet capsule forces in a probabilistic finite element analysis

PURPOSE

After total disc replacement with a ball-and-socket joint, reduced range of motion and progression of facet joint degeneration at the index level have been described. The aim of the study was to test the hypothesis that misalignment of the vertebrae adjacent to the implant reduces range of motion and increases facet joint or capsule tensile forces.

METHODS

A probabilistic finite element analysis was performed using a lumbosacral spine model with an artificial disc at level L5/S1. Misalignment of the L5 vertebra, the gap size of the facet joints, the transection of the posterior longitudinal ligament, and the spinal shape were varied. The model was loaded with pure moments.

RESULTS

Misalignment of the L5 vertebra reduced the range of motion up to 2°. A 2-mm displacement of the L5 vertebra in the anterior direction already led to facet joint forces of approximately 240 N. Extension, lateral bending, and axial rotation caused maximum facet joint forces between 280 and 380 N, while flexion caused maximum forces of approximately 200 N. A 2-mm displacement in the posterior direction led to capsule forces of approximately 80 N. Additional moments increased the maximum facet capsule forces to values between 120 and 230 N.

CONCLUSIONS

Misalignment of the vertebrae adjacent to an artificial disc strongly increases facet joint or capsule forces. It might, therefore, be an important reason for unsatisfactory clinical results. In an associated clinical study (Part 2), these findings are validated.

Comparison of four reconstruction methods after total sacrectomy: a finite element study

BACKGROUND

After total sacrectomy, it is mandatory to reconstruct the continuity between the lumbar spine and the pelvis. Only few biomechanical analyses exist which compare different reconstructions. Therefore, the aim of this study was to compare the lumbo-pelvic motion and the relative risk of implant breakage for four different reconstructions after total sacrectomy.

METHOD

Finite element analyses were performed for four general different reconstructions after total sacrectomy: sacral-rod reconstruction, four-rod reconstruction, bilateral fibular flaps reconstruction, and improved compound reconstruction. The rotations between L5 vertebra and ilium, the L5 shift-down displacement, and the maximum von Mises stress in the implants were calculated and evaluated for flexion, extension, lateral bending and axial rotation.

FINDINGS

The decreasing order of the rotations between L5 vertebra and ilium as well as of the L5 shift-down displacement for the studied reconstruction methods was four-rod reconstruction>sacral-rod reconstruction>bilateral fibular flaps reconstruction>improved compound reconstruction. The decreasing order of the maximum von Mises stress in the implants was sacral-rod reconstruction>four-rod reconstruction>bilateral fibular flaps reconstruction>improved compound reconstruction.

INTERPRETATION

From the mechanical point of view, improved compound reconstruction is superior to the other methods studied here as it shows the highest stability and the lowest maximum von Mises stress. However, clinical aspects must also be regarded when choosing a reconstruction method for a specific patient.

Optimal stiffness of a pedicle-screw-based motion preservation implant for the lumbar spine

PURPOSE

Pedicle-screw-based dynamic implants are intended to preserve intervertebral mobility while releasing certain spinal structures. The aim of the study was to determine the as yet unknown optimal stiffness value of the longitudinal rods that fulfils best these opposing tasks.

METHODS

A finite element model of the lumbar spine was used which includes the posterior implant at level L4/5. More than 250 variations of this model were generated by varying the diameter of the longitudinal rods between 6 and 12 mm and their elastic modulus between 10 MPa and 200 MPa. The loading cases flexion, extension, lateral bending and axial rotation were simulated. Evaluated optimization criteria were the ranges of motion, forces in the facet joints, posterior bulgings of the intervertebral disc and the intradiscal pressures. Various objective functions were evaluated.

RESULTS

The results show that the objective values depend more on the axial stiffness of the rods than on bending and torsional stiffness, rod diameter and elastic modulus. The optimal stiffness value for most of the investigated objective functions is approximately 50 N/mm and is achieved, e.g. using a rod diameter of 6 mm and an elastic modulus of 50 MPa. The design with the least axial stiffness was the best one with regard to the mobility. The forces in the facet joints and the intradiscal pressures were reduced mostly by an implant with the highest axial stiffness. When minimal posterior disc bulging was the criterion, the optimal axial stiffness was also approximately 50 N/mm.

CONCLUSIONS

The optimal axial stiffness of a pedicle-screw-based motion preservation implant for the lumbar spine is approximately 50 N/mm.