Biomechanical response of lumbar facet joints under follower preload: a finite element study

The effect of different fixation systems on oblique lumbar interbody fusion under vibration conditions

Despite the fact that lower back pain caused by degenerative lumbar spine pathologies seriously affects the quality of life, however, there is a paucity of research on the biomechanical properties of different auxiliary fixation systems for its primary treatment (oblique lumbar interbody fusion) under vibratory environments. In order to study the effects of different fixation systems of OLIF surgery on the vibration characteristics of the human lumbar spine under whole-body vibration (WBV), a finite element (FE) model of OLIF surgery with five different fixation systems was established by modifying a previously established model of the normal lumbar spine (L1-S1). In this study, a compressive follower load of 500 N and a sinusoidal axial vertical load of ±40 N at the frequency of 5 Hz with a duration of 0.6 s was applied. The results showed that the bilateral pedicle screw fixation model had the highest resistance to cage subsidence and maintenance of disc height under WBV. In contrast, the lateral plate fixation model exerted very high stresses on important tissues, which would be detrimental to the patient's late recovery and reduction of complications. Therefore, this study suggests that drivers and related practitioners who are often in vibrating environments should have bilateral pedicle screws for OLIF surgery, and side plates are not recommended to be used as a separate immobilization system. Additionally, the lateral plate is not recommended to be used as a separate fixation system.

The impact of disc degeneration on the dynamic characteristics of the lumbar spine: a finite element analysis

The dynamics of disc degeneration was analyzed to determine the effect of disc degeneration at the L4-L5 segment on the dynamic characteristics of the total lumbar spine. A three-dimensional nonlinear finite element model of the L1-S1 normal lumbar spine was constructed and validated. This normal model was then modified to construct two degeneration models with different degrees of degeneration (mild, moderate) at the L4-L5 level. Modal analysis, harmonic response analysis, and transient dynamics analysis were performed on the total lumbar spine when experiencing following compressive loading (500 N). As the degree of disc degeneration increased, the vibration patterns corresponding to the first three orders of the model's intrinsic frequency were basically unchanged, with the first order being in the left-right lateral bending direction, the second order being in the forward-flexion and backward-extension direction, and the third order being in the axial stretching direction. The nucleus pulposus pressure peaks corresponding to the first order intrinsic frequency for the harmonic response analysis are all on the right side of the model, with sizes of 0.053 MPa, 0.061 MPa, and 0.036 MPa, respectively; the nucleus pulposus pressure peaks corresponding to the second order intrinsic frequency are all at the rear of the model, with sizes of 0.13 MPa, 0.087 MPa, and 0.11 MPa, respectively; and the nucleus pulposus pressure peaks corresponding to the third order intrinsic frequency are all at the front of the model, with sizes of 0.19 MPa, 0.22 MPa, and 0.22 MPa, respectively. The results of the transient analysis indicated that over time, the response curves of the healthy model, the mild model, and the moderate model all exhibited cyclic response characteristics. Intervertebral disc degeneration did not adversely affect the vibration characteristics of the entire lumbar spine system. Intervertebral disc degeneration significantly altered the dynamics of the degenerative segments and their neighboring normal segments. The process of disc degeneration gradually shifted the load from the nucleus pulposus to the annulus fibrosus when the entire lumbar spine was subjected to the same vibratory environment.

Biomechanical response of decompression alone in lower grade lumbar degenerative spondylolisthesis--A finite element analysis

BACKGROUND

Previous studies have demonstrated the clinical efficacy of decompression alone in lower-grade spondylolisthesis. A higher rate of surgical revision and a lower rate of back pain relief was also observed. However, there is a lack of relevant biomechanical evidence after decompression alone for lower-grade spondylolisthesis.

PURPOSE

Evaluating the biomechanical characteristics of total laminectomy, hemilaminectomy, and facetectomy for lower-grade spondylolisthesis by analyzing the range of motion (ROM), intradiscal pressure (IDP), annulus fibrosus stress (AFS), facet joints contact force (FJCF), and isthmus stress (IS).

METHODS

Firstly, we utilized finite element tools to develop a normal lumbar model and subsequently constructed a spondylolisthesis model based on the normal model. We then performed total laminectomy, hemilaminectomy, and one-third facetectomy in the normal model and spondylolisthesis model, respectively. Finally, we analyzed parameters, such as ROM, IDP, AFS, FJCF, and IS, for all the models under the same concentrate force and moment.

RESULTS

The intact spondylolisthesis model showed a significant increase in the relative parameters, including ROM, AFS, FJCF, and IS, compared to the intact normal lumbar model. Hemilaminectomy and one-third facetectomy in both spondylolisthesis and normal lumbar models did not result in an obvious change in ROM, IDP, AFS, FJCF, and IS compared to the pre-operative state. Moreover, there was no significant difference in the degree of parameter changes between the spondylolisthesis and normal lumbar models after undergoing the same surgical procedures. However, total laminectomy significantly increased ROM, AFS, and IS and decreased the FJCF in both normal lumbar models and spondylolisthesis models.

CONCLUSION

Hemilaminectomy and one-third facetectomy did not have a significant impact on the segment stability of lower-grade spondylolisthesis; however, patients with LDS undergoing hemilaminectomy and one-third facetectomy may experience higher isthmus stress on the surgical side during rotation. In addition, total laminectomy changes the biomechanics in both normal lumbar models and spondylolisthesis models.

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.

Analyzing isolated degeneration of lumbar facet joints: implications for degenerative instability and lumbar biomechanics using finite element analysis

The facet joint contributes to lumbar spine stability as it supports the weight of body along with the intervertebral discs. However, most studies on the causes of degenerative lumbar diseases focus on the intervertebral discs and often overlook the facet joints. This study aimed to investigate the impact of facet joint degeneration on the degenerative changes and diseases of the lumbar spine. A finite element model of the lumbar spine (L1-S1) was fabricated and validated to study the biomechanical characteristics of the facet joints. To simulate degeneration of the facet joint, the model was divided into four grades based on the number of degenerative segments (L4-L5 or L4-S1) and the contact condition between the facet joint surfaces. Finite element analysis was performed on four spine motions: flexion, extension, lateral bending, and axial torsion, by applying a pure moment to the upper surface of L1. Important parameters that could be used to confirm the effect of facet joint degeneration on the lumbar spine were calculated, including the range of motion (ROM) of the lumbar segments, maximum von Mises stress on the intervertebral discs, and reaction force at the facet joint. Facet joint degeneration affected the biomechanical characteristics of the lumbar spine depending on the movements of the spine. When analyzed by dividing it into degenerative onset and onset-adjacent segments, lumbar ROM and the maximum von Mises stress of the intervertebral discs decreased as the degree of degeneration increased in the degenerative onset segments. The reaction force at the facet joint decreased with flexion and increased with lateral bending and axial torsion. In contrast, lumbar ROM of the onset-adjacent segments remained almost unchanged despite severe degeneration of the facet joint, and the maximum von Mises stress of the intervertebral discs increased with flexion and extension but decreased with lateral bending and axial torsion. Additionally, the facet joint reaction force increased with extension, lateral bending, and axial rotation. This analysis, which combined the ROM of the lumbar segment, maximum von Mises stress on the intervertebral disc, and facet joint reaction force, confirmed the biomechanical changes in the lumbar spine due to the degeneration of isolated facet joints under the load of spinal motion. In the degenerative onset segment, spinal instability decreased, whereas in the onset-adjacent segment, a greater load was applied than in the intact state. When conducting biomechanical studies on the lumbar spine, considering facet joint degeneration is important since it can lead to degenerative spinal diseases, including adjacent segment diseases.

IVD fibrosis and disc collapse comprehensively aggravate vertebral body disuse osteoporosis and zygapophyseal joint osteoarthritis by posteriorly shifting the load transmission pattern

BACKGROUND

Disuse is a typical phenotype of osteoporosis, but the underlying mechanism has yet to be identified in elderly patients. Disc collapse and intervertebral disc (IVD) fibrosis are two main pathological changes in IVD degeneration (IDD) progression, given that these changes affect load transmission patterns, which may lead to disuse osteoporosis of vertebral bodies and zygapophyseal joint (ZJ) osteoarthritis (ZJOA) biomechanically.

METHODS

Clinical data from 59 patients were collected retrospectively. Patient vertebral bony density, ZJOA grade, and disc collapse status were judged via CT. The IVD fibrosis grade was determined based on the FA measurements. Regression analyses identified potential independent risk factors for osteoporosis and ZJOA. L4-L5 numerical models with and without disc collapse and IVD fibrosis were constructed; stress distributions on the bony endplate (BEP) and zygapophyseal joint (ZJ) cartilages were computed in models with and without disc collapse and IVD fibrosis.

RESULTS

A significantly lower disc height ratio and significantly greater FA were recorded in patients with ZJOA. A significant correlation was observed between lower HU values and two parameters related to IDD progression. These factors were also proven to be independent risk factors for both osteoporosis and ZJOA. Correspondingly, compared to the intact model without IDD. Lower stress on vertebral bodies and greater stress on ZJOA can be simultaneously recorded in models of disc collapse and IVD fibrosis.

CONCLUSIONS

IVD fibrosis and disc collapse simultaneously aggravate vertebral body disuse osteoporosis and ZJOA by posteriorly shifting the load transmission pattern.

Biomechanical Evaluation of 2 Endoscopic Spine Surgery Methods for Treating Lumbar Disc Herniation: A Finite Element Study

OBJECTIVE

This study aimed to evaluate the effects of 2 endoscopic spine surgeries on the biomechanical properties of normal and osteoporotic spines.

METHODS

Based on computed tomography images of a healthy adult volunteer, 6 finite element models were created. After validating the normal intact model, a concentrated force of 400 N and a moment of 7.5 Nm were exerted on the upper surface of L3 to simulate 6 physiological activities of the spine. Five types of indices were used to assess the biomechanical properties of the 6 models, range of motion (ROM), maximum displacement value, intervertebral disc stress, maximum stress value, and articular protrusion stress, and by combining them with finite element stress cloud.

RESULTS

In normal and osteoporotic spines, there was no meaningful change in ROM or disc stress in the 2 surgical models for the 6 motion states. Model N1 (osteoporotic percutaneous transforaminal endoscopic discectomy model) showed a decrease in maximum displacement value of 20.28% in right lateral bending. Model M2 (unilateral biportal endoscopic model) increased maximum displacement values of 16.88% and 17.82% during left and right lateral bending, respectively. The maximum stress value of L4-5 increased by 11.72% for model M2 during left rotation. In addition, using the same surgical approach, ROM, maximum displacement values, disc stress, and maximum stress values were more significant in the osteoporotic model than in the normal model.

CONCLUSION

In both normal and osteoporotic spines, both surgical approaches were less disruptive to the physiologic structure of the spine. Furthermore, using the same endoscopic spine surgery, normal spine biomechanical properties are superior to osteoporotic spines.

Biomechanical evaluation of Percutaneous endoscopic posterior lumbar interbody fusion and minimally invasive transforaminal lumbar interbody fusion: a biomechanical analysis

In order to analyze and evaluate the stability of lumbar spine and the risk of cage subsidence after different minimally invasive fusion operations, two finite element models Percutaneous endoscopic posterior lumbar interbody fusion (PE-PLIF) and minimally invasive transforaminal lumbar interbody Fusion (MIS-TLIF) were established. The results showed that compared with MIS-TLIF, PE-PLIF had better segmental stability, lower pedicle screw rod system stress, and lower risk of cage subsidence. The results suggest that the cage with appropriate height should be selected to ensure the segmental stability and avoid the risk of the subsidence caused by the cage with large height.

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.

Morphological Evaluation of Lumbar Facet Joints in Professional Baseball Players

Background

Many professional baseball players experience low back pain, a major cause of which is lumbar facet joint arthropathy.

Purpose

To evaluate the relationship between the dominant hand side and facet joint morphology in baseball movement.

Study Design

Cross-sectional study; Level of evidence, 3.

Methods

Participants were 25 Japanese professional baseball players (11 pitchers and 14 fielders) with low back pain and lower limb symptoms. Player age, hand dominance, and length of professional playing experience were recorded, and the lateral diameter of all lumbar facet joints was determined from the axial computed tomography scans. We defined the facet joints ipsilateral and contralateral to the dominant hand as dominant and nondominant, respectively. The nondominant-to-dominant (N/D) ratio of the lateral diameter was calculated, and differences between the pitchers and fielders were analyzed using the unpaired t test.

Results

The average player age and length of professional playing experience were 26.9 years (range, 19-37 years) and 7.2 years (range, 1-15 years), respectively. The right hand was dominant in 9 pitchers and 5 fielders, while the left hand was dominant in 2 pitchers and 9 fielders. In pitchers, the average lateral facet joint diameter on the nondominant side was significantly larger than on the dominant side at all vertebral levels except L1 to L2 (P < .05 for all). The N/D ratio for each facet joint was 1.06 (L1-L2), 1.11 (L2-L3), 1.10 (L3-L4), 1.12 (L4-L5), and 1.12 (L5-S1). In fielders, the average lateral facet joint diameter on the dominant side was significantly larger than on the nondominant side at L3 to L4 (P < .05), with N/D ratios of 0.98 (L1-L2), 0.96 (L2-L3), 0.94 (L3-L4), 0.97 (L4-L5), and 0.98 (L5-S1). The N/D ratio was significantly larger in pitchers than in fielders at all levels (P < .05 for all).

Conclusion

The facet joints of professional baseball players were enlarged asymmetrically, with different tendencies observed between pitchers and fielders. Although pitching and batting are movements that transmit the rotation from the lower limbs to the upper limbs, the effects of rotation and lateral flexion were associated significantly with facet joint hypertrophy.

Biomechanical Study of Minimally Invasive Nonfusion Surgery for Treatment of Disc Herniation Associated with Adjacent Segment Disease: A Finite Element Analysis

OBJECTIVE

We explored the biomechanical changes of 2 conventional minimally invasive nonfusion surgical methods for treating disc herniation in adjacent segment disease using 3-dimensional finite element analysis.

METHODS

A model comprising L3 to the sacrum was validated and used to establish an L4-L5 fusion model, and an adjacent segment disease (ASD) model was developed by modifying the material properties of the intervertebral discs. The ASD model was used to simulate 2 conventional minimally invasive nonfusion surgical methods, which resulted in the creation of 2 postoperative models (M1 and M2). The range of motion and the equivalent stress for each model were recorded under 6 different working conditions. The data are descriptive and were analyzed comparatively under a normal load.

RESULTS

Compared with the ASD group, the range of motion of the adjacent segment in the M1 and M2 groups remained unaffected. However, significant Von-Mises stress changes were found in the annulus fibrosus and nucleus pulposus (NP), especially during extension, ipsilateral bending, and rotation. Stress in the NP also shifted toward the surgical incision in the annulus fibrosus during these movements. The maximum Von-Mises stress in the NP of the cephalic segment increased more than did that of the caudal segment.

CONCLUSIONS

Minimal nonfusion surgery for ASD might not affect adjacent segment stability significantly. Nonetheless, it can lead to segmental degeneration deterioration and postoperative recurrence. The cephalic segment is affected more than the caudal segment. Therefore, consideration of disc degeneration and appropriate selection of surgical methods for ASD are crucial.

Biomechanical effect of intervertebral disc degeneration on the lower lumbar spine

Lumbar intervertebral disc degeneration can induce bone hyperplasia, lumbar intervertebral disc herniation and other diseases, is one of the causes of low back pain, which seriously affects people's quality of life. And the causes of degeneration are very complex, so it is essential to understand the underlying mechanism of intervertebral disc degeneration and its influence. In this study, biomechanical effects of L4∼L5 lumbar degeneration with different degrees of degeneration were studied based on the numerical simulations. The three-dimensional finite element model of normal L2∼S1 lumbar vertebrae was established based on CT images of average adult male and verified. Several key parameters (intervertebral disc height, nucleus pulposus size, properties of different materials, etc.) of the model were modified to construct L4∼L5 models with different degrees of degeneration (grade 1, grade 2, grade 3, and grade 4). The range of motion (ROM), the intradiscal pressure of the nucleus, and the maximum Von Mises stress were determined by applying torques in different directions to simulate the four postures of flexion, extension, lateral bending, and axial rotation under compression load (500 N) to simulate the upper body weight of the human body. In different postures, with the increase of L4∼L5 degeneration degree, the ROM of the L4∼L5 degeneration segment showed a decreasing trend (Grade 4 had decrease of 41.9% to 65.2% compared to normal at different postures), while the ROM of its adjacent normal segments showed an increasing trend (L3∼L4: Grade 4 had increase of 21%-94% compared to normal at different postures; L5∼S1: Grade 4 had increase of 32%-66% compared to normal at different postures). With the increase in the degree of degeneration, nucleus pulposus pressure in the L4∼L5 degeneration segment decreased continuously under different postural conditions (Grade 4 had decrease of 25%-134.6% compared to normal at different postures), while the nucleus pulposus pressure in adjacent normal segments (L3∼L4 and L5∼S1) showed a gradually increasing trend. The maximum Von Mises stress of the three segments increased with the increasing degree of degeneration at different postures (L4∼L5: Grade 4 increased to 1.75 ∼ 4 times compared to normal at different postures). In four different models of lumbar disc degeneration, the adjacent normal segment of the disc compensates for the movement and loading pattern of the degenerated segment. At the same time, the load pattern inside the degenerated segment also changes.

Annulus Calibration Increases the Computational Accuracy of the Lumbar Finite Element Model

STUDY DESIGN

Mechanical simulations.

OBJECTIVE

Inadequate calibration of annuli negatively affects the computational accuracy of finite element (FE) models. Specifically, the definition of annulus average radius (AR) does not have uniformity standards. Differences between the elastic moduli in the different layers and parts of the annulus were not fully calibrated when a linear elastic material is used to define its material properties. This study aims to optimize the computational accuracy of the FE model by calibrating the annulus.

METHODS

We calibrated the annulus AR and elastic modulus in our anterior-constructed lumbar model by eliminating the difference between the computed range of motion and that measured by in vitro studies under a flexion-extension loading condition. Multi-indicator validation was performed by comparing the computed indicators with those measured in in vitro studies. The computation time required for the different models has also been recorded to evaluate the computational efficiency.

RESULTS

The difference between computed and measured ROMs was less than 1% when the annulus AR and elastic modulus were calibrated. In the model validation process, all the indicators computed by the calibrated FE model were within ±1 standard deviation of the average values obtained from in vitro studies. The maximum difference between the computed and measured values was less than 10% under nearly all loading conditions. There is no apparent variation tendency for the computational time associated with different models.

CONCLUSION

The FE model with calibrated annulus AR and regional elastic modulus has higher computational accuracy and can be used in subsequent mechanical studies.

Biomechanical analysis of functionally graded porous interbody cage for lumbar spinal fusion

Functionally graded porous (FGP) interbody cage might offer a trade-off between porosity-based reduction of stiffness and mechanical properties. Using finite element models of intact and implanted lumbar functional spinal unit (FSU), the study investigated the quantitative deviations in load transfer and adaptive changes in bone density distributions around FGP interbody cages. The cage had three graded porosities: FGP-A, -B, and -C corresponded to a maximum porosity levels of 48%, 65% and 78%, respectively. Efficacy of the FGP cages were evaluated by comparing the numerically predicted results of solid-Ti and uniformly porous 78% porosity (P78) cage. Variations in stiffness and interface condition affected the strain distribution and bone remodelling around the cages. Peak strains of 0.5-1% were observed in less number of peri-prosthetic bone elements for the FGP cages as compared to the solid-Ti cage. Strains and bone apposition were considerably higher for the bonded implant-bone interface condition than the debonded case. For the FGP-C with bonded interface condition, bone apposition of 11-20% was predicted in the L4 and L5 regions of interest (ROIs); whereas the debonded model exhibited 6-10% increase in bone density. The deviations in bone density change between FGP-C and P78 model were 3-8% for L4 and L5 ROIs. FGP resulted in a reduced average micromotion (∼70-106 μm) as compared to solid-Ti (116 μm), for all physiologic movements. Compared to solid-Ti and uniformly porous cages, the FGP cage seems to be a viable alternative considering the conflicting nature of strength and porosity.

Influence of weight-bearing on the 3D movement of lumbar facet joints in the sitting position

OBJECTIVE

To analyze the motion characteristics of lumbar facet joints and to observe the effect of weight-bearing on lumbar facet joints in the sitting position.

METHODS

Ten normal subjects (5 males and 5 females) were recruited and scanned by CT, and their lumbar 3D models were reconstructed by software. The images of flexion and extension of lumbar facet joints in the sitting position were collected without weight-bearing and weight-bearing 10 kg, and the 2D model was constructed by software. The 2D-3D model was matched to restore the flexion and extension motion changes of the subjects' lumbar spine in the sitting position. Coordinates were established in the middle of the vertebral body and copied to the facet joints. Measure and record the lumbar facet joint movement distance through coordinate system. The relevant data of facet joints were collected.

RESULTS

In the L3/4 segment, after weight loading, the displacement of the left facet joint in the X axis became larger, while that in the Y axis and Z axis decreased. The displacement of the right facet joint in the X axis and Y axis increased, and the Z axis displacement decreased. The rotation angle of the bilateral facet joints also decreased. In the L4/5 segment, after loading, the displacements of the X, Y, and Z axis displacements of both sides increase, while the rotation angles of α and β increase, while the rotation angle of γ decreases. In the L5/S1 segment, the displacements of the X, Y, and Z axes on the left side decrease. The displacement of the X and Y axes on the right side decreases, while the displacement on the Z axis increases. The rotation angles of α and γ increase, and the rotation angle of the β axis decreases.

CONCLUSION

When sitting, the flexion and extension distance and rotational displacement of lumbar facet joints are not affected by weight-bearing. In addition, there is asymmetry in the movement of the left and right facet joints, and weight bearing has no effect on the asymmetry of the motion.

Biomechanical study of two-level oblique lumbar interbody fusion with different types of lateral instrumentation: a finite element analysis

Objective

The aim of this study was to verify the biomechanical properties of a newly designed angulated lateral plate (mini-LP) suited for two-level oblique lumbar interbody fusion (OLIF). The mini-LP is placed through the lateral ante-psoas surgical corridor, which reduces the operative time and complications associated with prolonged anesthesia and placement in the prone position.

Methods

A three-dimensional nonlinear finite element (FE) model of an intact L1-L5 lumbar spine was constructed and validated. The intact model was modified to generate a two-level OLIF surgery model augmented with three types of lateral fixation (stand-alone, SA; lateral rod screw, LRS; miniature lateral plate, mini-LP); the operative segments were L2-L3 and L3-L4. By applying a 500 N follower load and 7.5 Nm directional moment (flexion-extension, lateral bending, and axial rotation), all models were used to simulate human spine movement. Then, we extracted the range of motion (ROM), peak contact force of the bony endplate (PCFBE), peak equivalent stress of the cage (PESC), peak equivalent stress of fixation (PESF), and stress contour plots.

Results

When compared with the intact model, the SA model achieved the least reduction in ROM to surgical segments in all motions. The ROM of the mini-LP model was slightly smaller than that of the LRS model. There were no significant differences in surgical segments (L1-L2, L4-L5) between all surgical models and the intact model. The PCFBE and PESC of the LRS and the mini-LP fixation models were lower than those of the SA model. However, the differences in PCFBE or PESC between the LRS- and mini-LP-based models were not significant. The fixation stress of the LRS- and mini-LP-based models was significantly lower than the yield strength under all loading conditions. In addition, the variances in the PESF in the LRS- and mini-LP-based models were not obvious.

Conclusion

Our biomechanical FE analysis indicated that LRS or mini-LP fixation can both provide adequate biomechanical stability for two-level OLIF through a single incision. The newly designed mini-LP model seemed to be superior in installation convenience, and equally good outcomes were achieved with both LRS and mini-LP for two-level OLIF.

Effects of osteoporosis on the biomechanics of various supplemental fixations co-applied with oblique lumbar interbody fusion (OLIF): a finite element analysis

Background

Oblique lumbar interbody fusion (OLIF) is an important surgical modality for the treatment of degenerative lumbar spine disease. Various supplemental fixations can be co-applied with OLIF, increasing OLIF stability and reducing complications. However, it is unclear whether osteoporosis affects the success of supplemental fixations; therefore, this study analyzed the effects of osteoporosis on various supplemental fixations co-applied with OLIF.

Methods

We developed and validated an L3-S1 finite element (FE) model; we assigned different material properties to each component and established models of the osteoporotic and normal bone lumbar spine. We explored the outcomes of OLIF combined with each of five supplemental fixations: standalone OLIF; OLIF with lateral plate fixation (OLIF + LPF); OLIF with translaminar facet joint fixation and unilateral pedicle screw fixation (OLIF + TFJF + UPSF); OLIF with unilateral pedicle screw fixation (OLIF + UPSF); and OLIF with bilateral pedicle screw fixation (OLIF + BPSF). Under the various working conditions, we calculated the ranges of motion (ROMs) of the normal bone and osteoporosis models, the maximum Mises stresses of the fixation instruments (MMSFIs), and the average Mises stresses on cancellous bone (AMSCBs).

Results

Compared with the normal bone OLIF model, no demonstrable change in any segmental ROM was apparent. The MMSFIs increased in all five osteoporotic OLIF models. In the OLIF + TFJF + UPSF model, the MMSFIs increased sharply in forward flexion and extension. The stress changes of the OLIF + UPSF, OLIF + BPSF, and OLIF + TFJF + UPSF models were similar; all stresses trended upward. The AMSCBs decreased in all five osteoporotic OLIF models during flexion, extension, lateral bending, and axial rotation. The average stress change of cancellous bone was most obvious under extension. The AMSCBs of the five OLIF models decreased by 14%, 23.44%, 21.97%, 40.56%, and 22.44% respectively.

Conclusions

For some supplemental fixations, the AMSCBs were all reduced and the MMSFIs were all increased in the osteoporotic model, compared with the OLIF model of normal bone. Therefore, the biomechanical performance of an osteoporotic model may be inferior to the biomechanical performance of a normal model for the same fixation method; in some instances, it may increase the risks of fracture and internal fixation failure.

Biomechanical effect of C5-C6 intervertebral disc degeneration on the human lower cervical spine (C3-C7): a finite element study

The biomechanical effects of intervertebral discs and facet joints degeneration on the cervical spine are essential to understanding the mechanisms of spinal disorders to improve pathological and clinical treatment. In this study, the biomechanical effects of a progressively degenerated C5-C6 segment on the human lower cervical spine are determined by a detailed simulation of intervertebral disc degeneration. A detailed asymmetric three-dimension intact finite element model was developed using computed tomography scan data of the human lower cervical spine (C3-C7). The intact finite element model was then modified at the C5-C6 segment to build three degenerated models, such as mild, moderate, and severe degeneration. The physiological compressive load 73.6 N, and moment 1 Nm were applied at the superior endplate of the vertebra C3, and the inferior endplate of the C7 vertebra was a constraint for all degrees of freedom. Range of motion, maximum von Mises stress in the annulus, intradiscal pressure, and facet joint force of the degenerated models were computed. With progressive degeneration in the C5-C6 segment, the range of motion of degenerated and normal segments decreases in all postures. Intradiscal pressure of the degenerated segment decreases but increases in normal segments of degenerated segment C5-C6, and facet joint forces increase at both degenerated and normal segments. This study emphasizes that the degenerated disc alters the degenerated and normal segments' motion and loading patterns. The abnormal increase in facet joint force in the degenerated models threatened to accelerate the degeneration in the normal segments.

What happens at the L5/S1 facet joint when implants are placed across the sacroiliac joint?

BACKGROUND

Injuries to the posterior pelvic ring are often stabilized with fixation across the sacroiliac joint (SIJ). However, the compensatory changes at the neighboring L5/S1 facet joint are unknown. The objective of this study was to determine the compensatory change in pelvic kinematics and contact forces at the L5/S1 facet joint after fixation across the sacroiliac joint (SIJ) using a cadaveric model.

METHODS

Five fresh-frozen cadaveric pelvis specimens were dissected to remove non-structural soft tissue. Retroreflective markers were fixed to the L5 body, S1 body and bilateral anterior superior iliac spines to represent the motion of L5, S1 and the ileum, respectively. Pressure sensors were inserted in both L5/S1 facet joints. Testing was performed using a robotic system that applied load to mimic ambulation. Testing was performed prior to SIJ fixation, after unilateral SIJ fixation and bilateral fixation.

RESULTS

Contact force at the L5/S1 facet joint significantly increased by 55% from 48.4 N to 75.2 N following unilateral fixation (p = 0.0161) and increased by 100% to 96.9 N after bilateral fixation (p = 0.0038). Unilateral SIJ fixation increased flexion of the ilium relative to L5 from 1.2° to 2.0° (p = 0.01) and increased axial rotation of L5 relative to S1 from 0.7° to 1.6° (p = 0.001). Bilateral fixation increased flexion of the ilium relative to L5 to 2.0° from 1.2° prior to fixation (p = 0.001), increased axial rotation of L5 relative to S1 to 1.2° from 0.7° prior to fixation (p = 0.002) and increased flexion of L5 relative to S1 to 2.4° from 1.5° prior to fixation (p = 0.04).

CONCLUSION

The L5/S1 facet joint experiences compensatory increased motion under increased contact force after unilateral and bilateral SIJ fixation, possibly predisposing it to adjacent segment arthritis.

LEVEL OF EVIDENCE

V, cadaveric study.

Biomechanical study of oblique lumbar interbody fusion (OLIF) augmented with different types of instrumentation: a finite element analysis

Background

To explore the biomechanical differences in oblique lumbar interbody fusion (OLIF) augmented by different types of instrumentation.

Methods

A three-dimensional nonlinear finite element (FE) model of an intact L3-S1 lumbar spine was built and validated. The intact model was modified to develop five OLIF surgery models (Stand-alone OLIF; OLIF with lateral plate fixation [OLIF + LPF]; OLIF with unilateral pedicle screws fixation [OLIF + UPSF]; OLIF with bilateral pedicle screws fixation [OLIF + BPSF]; OLIF with translaminar facet joint fixation + unilateral pedicle screws fixation [OLIF + TFJF + UPSF]) in which the surgical segment was L4–L5. Under a follower load of 500 N, a 7.5-Nm moment was applied to all lumbar spine models to calculate the range of motion (ROM), equivalent stress peak of fixation instruments (ESPFI), equivalent stress peak of cage (ESPC), equivalent stress peak of cortical endplate (ESPCE), and equivalent stress average value of cancellous bone (ESAVCB).

Results

Compared with the intact model, the ROM of the L4–L5 segment in each OLIF surgery model decreased by > 80%. The ROM values of adjacent segments were not significantly different. The ESPFI, ESPC, and ESPCE values of the OLIF + BPSF model were smaller than those of the other OLIF surgery models. The ESAVCB value of the normal lumbar model was less than the ESAVCB values of all OLIF surgical models. In most postures, the ESPFI, ESPCE, and ESAVCB values of the OLIF + LPF model were the largest. The ESPC was higher in the Stand-alone OLIF model than in the other OLIF models. The stresses of several important components of the OLIF + UPSF and OLIF + TFJF + UPSF models were between those of the OLIF + LPF and OLIF + BPSF models.

Conclusions

Our biomechanical FE analysis indicated the greater ability of OLIF + BPSF to retain lumbar stability, resist cage subsidence, and maintain disc height. Therefore, in the augmentation of OLIF, bilateral pedicle screws fixation may be the best approach.

The Mismatch Between Bony Endplates and Grafted Bone Increases Screw Loosening Risk for OLIF Patients With ALSR Fixation Biomechanically

The mismatch between bony endplates (BEPs) and grafted bone (GB) triggers several complications biomechanically. However, no published study has identified whether this factor increases the risk of screw loosening by deteriorating the local stress levels. This study aimed to illustrate the biomechanical effects of the mismatch between BEP and GB and the related risk of screw loosening. In this study, radiographic and demographic data of 56 patients treated by single segment oblique lumbar interbody fusion (OLIF) with anterior lateral single rod (ALSR) fixation were collected retrospectively, and the match sufficiency between BEP and GB was measured and presented as the grafted bony occupancy rate (GBOR). Data in patients with and without screw loosening were compared; regression analyses identified independent risk factors. OLIF with different GBORs was simulated in a previously constructed and validated lumbosacral model, and biomechanical indicators related to screw loosening were computed in surgical models. The radiographic review and numerical simulations showed that the coronal plane’s GBOR was significantly lower in screw loosening patients both in the cranial and caudal vertebral bodies; the decrease in the coronal plane’s GBOR has been proven to be an independent risk factor for screw loosening. In addition, numerical mechanical simulations showed that the poor match between BEP and GB will lead to stress concentration on both screws and bone-screw interfaces. Therefore, we can conclude that the mismatch between the BEP and GB will increase the risk of screw loosening by deteriorating local stress levels, and the increase in the GBOR by modifying the OLIF cage’s design may be an effective method to optimize the patient’s prognosis.

Biomechanical changes at the adjacent segments induced by a lordotic porous interbody fusion cage

Biomechanical changes at the adjacent segments after interbody fusion are common instigators of adjacent segment degeneration (ASD). This study aims to investigate how the presence of a lordotic porous cage affects the biomechanical performance of the adjacent segments. A finite element model (FEM) of a lumbar spine implanted with a lordotic cage at L3-L4 was validated by in-vitro testing. The stress distribution on the cage and range of motion (ROM) of L3-L4 were used to assess the stability of the implant. Three angles of cage (0° = non-restoration, 7° = normal restoration and 11° = over-restoration) were modelled with different porosities (0%, 30% and 60%) and evaluated in the motions of flexion, extension, lateral bending and rotation. The ROM, intervertebral disc pressure (IDP) and facet joint force (FJF) were used to evaluate biomechanical changes at the adjacent segments in each model. The results indicated that porous cages produced more uniform stress distribution, but cage porosity did not influence the ROM, IDP and FJF at L2-L3 and L4-L5. Increasing the cage lordotic angle acted to decrease the ROM and IDP, and increase the FJF of L4-L5, but did not alter the ROM of L2-L3. In conclusion, changes in ROM, IDP and FJF at the adjacent segments were mainly influenced by the lordotic angle of the cage and not by the porosity. A larger angle of lordotic cage was shown to reduce the ROM and IDP, and increase the FJF of the lower segment (L4-L5), but had little effect on the ROM of the upper segment (L2-L3).

Biomechanical comparison of vertebral augmentation and cement discoplasty for the treatment of symptomatic Schmorl's node: a finite element analysis

Percutaneous vertebral augmentation (PVA) and percutaneous cement discoplasty (PCD) are two relatively new minimally invasive surgeries for symptomatic Schmorl's reported in recent decade. However, the clinical evidence for the effectiveness of these two surgeries is insufficient. The purpose of this study was to compare the biomechanical benefits and risks of the two surgeries in order to analyze their biomechanical differences and effectiveness. We reconstructed Five lumbar finite element models via computed tomography data, including control model, PVA-ideal model, PVA-nonideal model, PCD-ideal model, and PCD-nonideal model. The stress and strain of Schmorl's nodes, bone marrow edema zone (BMEZ), affected endplate, and the overall stability of segment were analyzed and compared. The validity of our models was confirmed. As a result, the PVA-ideal model can significantly reduce the stress of Schmorl's node and the strain of BMEZ, while this effect is inappreciable in PVA-nonideal model. The PCD-ideal model significantly reduced the strain of Schmorl's nodes and BMEZ, and significantly improve segmental stability, but also resulted in a significant increase in the stress of Schmorl's nodes, BMEZ and endplates. The PCD-nonideal model not only lacks blocking effect, but also sharply increases the strain of Schmorl's nodes and BMEZ. Thus, We recommend that both PVA and PCD surgeries in ideal distribution facilitated a more stable paranodular biomechanical microenvironment. However, due to the possibility of poor biomechanical outcomes caused by the non-ideal cement distribution, the non-ideal distribution of bone cement needs to be remedied in practice.

Single-segment central lumbar spinal stenosis: Correlation with lumbar X-ray measurements

BACKGROUND

Lumbar X-rays are usually preferred in patients with lower back pain, but lumbar spinal stenosis (LSS) cannot be directly observed on lumbar X-ray films.

OBJECTIVE

The purpose of this study is to explore the correlation between the degree of single-segment central LSS and lumbar X-ray measurements.

METHODS

The data of 60 male patients aged 39-78 years with single-segment central LSS were analyzed. Linear correlation analysis was used to determine the correlation between the single-segment central LSS and the various measurement parameters. Multiple linear regression analysis was used to analyze the factors affecting single-segment central LSS.

RESULTS

There were significant differences in S1/S0, E, B, L1-5Cobb, and M among the three groups (p< 0.05). S1/S0 was positively correlated with E, B, L1-5Cobb, and M (p< 0.05), but was not correlated with D (p= 0.66). After multiple linear regression analysis, B, L1-5Cobb, and M were independently associated with S1/S0.

CONCLUSIONS

The B, L1-5Cobb, and M parameters were independently associated with single-stage central LSS, and would likely be of particular value in evaluating the degree of single-segment central LSS; B, L1-5Cobb, and M served as independent predictors of the degree of LSS. These findings will guide clinicians' decision-making in the future.

TELD with limited foraminoplasty has potential biomechanical advantages over TELD with large annuloplasty: an in-silico study

Background

Facetectomy, an important procedure in the in–out and out–in techniques of transforaminal endoscopic lumbar discectomy (TELD), is related to the deterioration of the postoperative biomechanical environment and poor prognosis. Facetectomy may be avoided in TELD with large annuloplasty, but iatrogenic injury of the annulus and a high grade of nucleotomy have been reported as risk factors influencing poor prognosis. These risk factors may be alleviated in TELD with limited foraminoplasty, and the grade of facetectomy in this surgery can be reduced by using an endoscopic dynamic drill.

Methods

An intact lumbo-sacral finite element (FE) model and the corresponding model with adjacent segment degeneration were constructed and validated to evaluate the risk of biomechanical deterioration and related postoperative complications of TELD with large annuloplasty and TELD with limited foraminoplasty. Changes in various biomechanical indicators were then computed to evaluate the risk of postoperative complications in the surgical segment.

Results

Compared with the intact FE models, the model of TELD with limited foraminoplasty demonstrated slight biomechanical deterioration, whereas the model of TELD with large annuloplasty revealed obvious biomechanical deterioration. Degenerative changes in adjacent segments magnified, rather than altered, the overall trends of biomechanical change.

Conclusions

TELD with limited foraminoplasty presents potential biomechanical advantages over TELD with large annuloplasty. Iatrogenic injury of the annulus and a high grade of nucleotomy are risk factors for postoperative biomechanical deterioration and complications of the surgical segment.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Misaligned spinal rods can induce high internal forces consistent with those observed to cause screw pullout and disc degeneration

BACKGROUND CONTEXT

Manual contouring of spinal rods is often required intraoperatively for proper alignment of the rods within the pedicle screw heads. Residual misalignments are frequently reduced by using dedicated reduction devices. The forces exerted by these devices, however, are uncontrolled and may lead to excessive reaction forces. As a consequence, screw pullout might be provoked and surrounding tissue may experience unfavorable biomechanical loads. The corresponding loads and induced tissue deformations are however not well identified. Additionally, whether the forced reduction alters the biomechanical behavior of the lumbar spine during physiological movements postoperatively, remains unexplored.

PURPOSE

To predict whether the reduction of misaligned posterior instrumentation might result in clinical complications directly after reduction and during a subsequent physiological flexion movement.

STUDY DESIGN

Finite element analysis.

METHODS

A patient-specific, total lumbar (L1-S1) spine finite element model was available from previous research. The model consists of poro-elastic intervertebral discs with Pfirrmann grade-dependent material parameters, with linear elastic bone tissue with stiffness values related to the local bone density, and with the seven major ligaments per spinal motion segment described as nonlinear materials. Titanium instrumentation was implemented in this model to simulate a L4, L5, and S1 posterolateral fusion. Next, coronal and sagittal misalignments of 6 mm each were introduced between the rod and the screw head at L4. These misalignments were computationally reduced and a physiological flexion movement of 15° was prescribed. Non-instrumented and well-aligned instrumented models were used as control groups.

RESULTS

Pulling forces up to 1.0 kN were required to correct the induced misalignments of 6 mm. These forces affected the posture of the total lumbar spine, as motion segments were predicted to rotate up to 3 degrees and rotations propagated proximally to and even affect the L1-2 level. The facet contact pressures in the corrected misaligned models were asymmetrical suggesting non-physiological joint loading in the misaligned models. In addition, the discs and vertebrae experienced abnormally high forces as a result of the correction procedure. These effects were more pronounced after a 15° flexion movement following forced reduction.

CONCLUSIONS

The results of this study indicate that the correction of misaligned posterior instrumentation can result in high forces at the screws consistent with those reported to cause screw pullout, and may cause high-tissue strains in adjacent and downstream spinal segments.

CLINICAL SIGNIFICANCE

Proper alignment of spinal posterior instrumentation may reduce clinical complications secondary to unfavorable biomechanics.

Automated Pipeline to Generate Anatomically Accurate Patient-Specific Biomechanical Models of Healthy and Pathological FSUs

State-of-the-art preoperative biomechanical analysis for the planning of spinal surgery not only requires the generation of three-dimensional patient-specific models but also the accurate biomechanical representation of vertebral joints. The benefits offered by computational models suitable for such purposes are still outweighed by the time and effort required for their generation, thus compromising their applicability in a clinical environment. In this work, we aim to ease the integration of computerized methods into patient-specific planning of spinal surgery. We present the first pipeline combining deep learning and finite element methods that allows a completely automated model generation of functional spine units (FSUs) of the lumbar spine for patient-specific FE simulations (FEBio). The pipeline consists of three steps: (a) multiclass segmentation of cropped 3D CT images containing lumbar vertebrae using the DenseVNet network, (b) automatic landmark-based mesh fitting of statistical shape models onto 3D semantic segmented meshes of the vertebral models, and (c) automatic generation of patient-specific FE models of lumbar segments for the simulation of flexion-extension, lateral bending, and axial rotation movements. The automatic segmentation of FSUs was evaluated against the gold standard (manual segmentation) using 10-fold cross-validation. The obtained Dice coefficient was 93.7% on average, with a mean surface distance of 0.88 mm and a mean Hausdorff distance of 11.16 mm (N = 150). Automatic generation of finite element models to simulate the range of motion (ROM) was successfully performed for five healthy and five pathological FSUs. The results of the simulations were evaluated against the literature and showed comparable ROMs in both healthy and pathological cases, including the alteration of ROM typically observed in severely degenerated FSUs. The major intent of this work is to automate the creation of anatomically accurate patient-specific models by a single pipeline allowing functional modeling of spinal motion in healthy and pathological FSUs. Our approach reduces manual efforts to a minimum and the execution of the entire pipeline including simulations takes approximately 2 h. The automation, time-efficiency and robustness level of the pipeline represents a first step toward its clinical integration.

Predicting Vertebral Bone Strength Using Finite Element Analysis for Opportunistic Osteoporosis Screening in Routine Multidetector Computed Tomography Scans-A Feasibility Study

Purpose

To investigate the feasibility of using routine clinical multidetector computed tomography (MDCT) scans for conducting finite element (FE) analysis to predict vertebral bone strength for opportunistic osteoporosis screening.

Methods

Routine abdominal MDCT with and without intravenous contrast medium (IVCM) of seven subjects (five male; two female; mean age: 71.86 ± 7.40 years) without any bone disease were used. FE analysis was performed on individual vertebrae (T11, T12, L1, and L2) including the posterior elements to investigate the effect of IVCM and slice thickness (1 and 3 mm) on vertebral bone strength. Another subset of data from subjects with vs. without osteoporotic vertebral fractures (n = 9 age and gender-matched pairs) was analyzed for investigating the ability of FE-analysis to differentiate the two cohorts. Bland-Altman plots, box plots, and coefficient of correlation (R2) were calculated to determine the variations in FE-predicted failure loads for different conditions.

Results

The FE-predicted failure loads obtained from routine MDCT scans were strongly correlated with those from without IVCM (R2 = 0.91 for 1mm; R2 = 0.92 for 3mm slice thickness, respectively) and different slice thicknesses (R2 = 0.93 for 1mm vs. 3mm with IVCM). Furthermore, a good correlation was observed for 3mm slice thickness with IVCM vs. 1mm without IVCM (R2 = 0.87). Significant difference between FE-predicted failure loads of healthy and fractured patients was observed (4,705 ± 1,238 vs. 4,010 ± 1,297 N; p=0.026).

Conclusion

Routine clinical MDCT scans could be reliably used for assessment of fracture risk based on FE analysis and may be beneficial for patients who are at increased risk for osteoporotic fractures.

Parametric study of anterior percutaneous endoscopic cervical discectomy (APECD)

Anterior percutaneous endoscopic cervical discectomy (APECD) is a common treatment for cervical spondylotic radiculopathy (CSR). In this study, the effects of various channel diameters and approach angles on cervical vertebrae on postoperative outcomes in APECD surgery were explored. A finite element model of intact cervical C₃-C₇ was constructed and then modified to obtain six surgical models. Range of motion (ROM) and intradiscal pressure (IDP) were calculated under different conditions of flexion (Fle), extension (Ext), lateral bending, and axial rotation. During Fle and bending to the left (LB), the ROM was closer to the intact model when the angle of approach was 90°. During bending to the left (LB) and rotation to the left (LR), the ROM changed considerably (43.2%, 33.7%, respectively) where the angle of approach was 45°. As the surgical channel diameter increased, the extent of the change in ROM compared with the intact model also increased. IDP decreased by 48% and 49%, respectively, compared with the intact model at the C₅-C₆ segment where the angle of approach was 45° and 60° during Fle, while it changed little at 90°, by less than 10%. The IDP was increased noticeably by 117.6%, 82.1%, and 105.8%, for channel diameters of 2, 3 and 4 mm, respectively. And declined noticeably during LB and LR (LB: 27.1%, 27.1%, 38.5%; LR: 37.4%, 35.5%, 48.7%). The results demonstrated that the shorter the surgical path, the smaller surgical diameter, the less the biomechanical influence on the cervical vertebra.

Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis

There is an increased interest in studying the biomechanics of the facet joints. For in silico studies, it is therefore important to understand the level of reliability of models for outputs of interest related to the facet joints. In this work, a systematic review of finite element models of multi-level spinal section with facet joints output of interest was performed. The review focused on the methodology used to model the facet joints and its associated validation. From the 110 papers analysed, 18 presented some validation of the facet joints outputs. Validation was done by comparing outputs to literature data, either computational or experimental values; with the major drawback that, when comparing to computational values, the baseline data was rarely validated. Analysis of the modelling methodology showed that there seems to be a compromise made between accuracy of the geometry and nonlinearity of the cartilage behaviour in compression. Most models either used a soft contact representation of the cartilage layer at the joint or included a cartilage layer which was linear elastic. Most concerning, soft contact models usually did not contain much information on the pressure-overclosure law. This review shows that to increase the reliability of in silico model of the spine for facet joints outputs, more needs to be done regarding the description of the methods used to model the facet joints, and the validation for specific outputs of interest needs to be more thorough, with recommendation to systematically share input and output data of validation studies.

Does oblique lumbar interbody fusion promote adjacent degeneration in degenerative disc disease: A finite element analysis

BACKGROUND

The number of oblique lumbar interbody fusion (OLIF) procedures has continued to rise over recent years. Adjacent segment degeneration (ASD) is a common complication following vertebral body fusion. Although the precise mechanism remains uncertain, ASD has gradually become more common in OLIF. Therefore, the present study analyzed the association between disc degeneration and OLIF to explore whether adjacent degeneration was promoted by OLIF in degenerative disc disease.

METHODS

A three-dimensional nonlinear finite element (FE) model of the L3-S1 lumbar spine was developed and validated. Three lumbar spine degeneration models with different degrees of degeneration (mild, moderate and severe) and a model of OLIF surgery were constructed at the L4-L5 level. When subjected to a follower compressive load (500 N), hybrid moment loading was applied to all models of the lumbar spine and the range of motion (ROM), intradiscal pressure (IDP), facet joint force (FJF), average mises stress in the annulus (AMSA), average tresca stress in the annulus (ATSA) and average endplate stress (AES) were measured.

RESULTS

Compared with the healthy lumbar spine model, the ROM, IDP, FJF, AMSA, ATSA and AES of the segments adjacent to the degenerated segment increased in each posture as the degree of disc degeneration increased. In different directions of motion, the ROM, IDP, FJF, AMSA, ATSA and AES in the OLIF model in the L3-L4 and L5-S1 segments were higher than those of the healthy model and each degenerated model. Compared with the healthy model, the largest relative increase in biomechanical parameters above (ROM, IDP, FJF, AMSA, ATSA or AES) was observed in the L3-L4 segment in the OLIF model, of 77.13%, 32.63%, 237.19%, 45.36%, 110.92% and 80.28%, respectively. In the L5-S1 segment the corresponding values were 68.88%, 36.12%, 147.24%, 46.00%, 45.88% and 51.29%, respectively.

CONCLUSIONS

Both degenerated discs and OLIF surgery modified the pattern of motion and load distribution of adjacent segments (L3-L4 and L5-S1 segments). The increases in the biomechanical parameters of segments adjacent to the surgical segment in the OLIF model were more apparent than those of the degenerated models. In summary, OLIF risked accelerating the degeneration of segments adjacent to those of a surgical segment.

The influence of the rib cage on the static and dynamic stability responses of the scoliotic spine

The thoracic cage plays an important role in maintaining the stability of the thoracolumbar spine. In this study, the influence of a rib cage on static and dynamic responses in normal and scoliotic spines was investigated. Four spinal finite element (FE) models (T1-S), representing a normal spine with rib cage (N1), normal spine without rib cage (N2), a scoliotic spine with rib cage (S1) and a scoliotic spine without rib cage (S2), were established based on computed tomography (CT) images, and static, modal, and steady-state analyses were conducted. In S2, the Von Mises stress (VMS) was clearly decreased compared to S1 for four bending loadings. N2 and N1 showed a similar VMS to each other, and there was a significant increase in axial compression in N2 and S2 compared to N1 and S1, respectively. The U magnitude values of N2 and S2 were higher than in N1 and S1 for five loadings, respectively. The resonant frequencies of N2 and S2 were lower than those in N1 and S1, respectively. In steady-state analysis, maximum amplitudes of vibration for N2 and S2 were significantly larger than N1 and S1, respectively. This study has revealed that the rib cage improves spinal stability in vibrating environments and contributes to stability in scoliotic spines under static and dynamic loadings.

Analysis and improvement of the three-column spinal theory

Background

Denis and Ferguson et al.’s three-column spinal theory has been widely accepted and applied. However, this three-column theory was proposed based solely on observation and experience without thorough documented data and analysis. The aim of this study was to analyze and improve Denis and Ferguson et al.’s three-column spinal theory to propose a novel three-column concept in epidemiology, morphology and biomechanics.

Methods

A retrospective analysis of the computed tomography imaging data of patients with a diagnosis of T11-L5 vertebral fractures was conducted between February 2010 and December 2018. Three-dimensional (3D) distribution maps of fracture lines of all subjects were obtained based on 3D mapping techniques. In addition, a 25-year-old health male volunteer was recruited for the vertebral finite element force analysis.

Results

The present study enrolled 459 patients (age: 48 ± 11.42 years), containing a total of 521 fractured vertebrae. The fracture lines peaked in the upper and the outer third sections of the vertebra, starting from the anterior part of the vertebral pedicles in 3-D maps. Regarding flexion and extension of the spine, the last third of the vertebral body in front of the spinal canal was one main stress center in the finite element analysis. The stress on the vertebral body was greater in front of the pedicles in the lateral bending.

Conclusion

The study reveals that the posterior one-third of the vertebral body in front of the spinal canal and the posterior one-third of the vertebral body in front of the pedicle are very different in terms of fracture characteristics and risks to spinal canal (3D maps and stress distributing graphs), therefore, they should be classified as different columns. We provide strong evidence that Su’s three-column theory complies with the characteristics of vertebral physiological structure, vertebral fracture, and vertebral biomechanics.

Graphical abstract

Effect of Percutaneous Endoscopic Lumbar Foraminoplasty of Different Facet Joint Portions on Lumbar Biomechanics: A Finite Element Analysis

Objective

To evaluate the influence of percutaneous endoscopic lumbar foraminoplasty of different facet joint portions on segmental range of motion (ROM) and intradiscal pressure (IDP) of L₃/L₄ and L₄/L₅ motion segments by establishing three dimensional finite element (FE) model.

Method

Computed tomography images of a male adult volunteer of appropriate age and in good condition both mentally and physically. Obtained data was used in this study from July 2020 to December 2020, and an intact L_3–5 three dimensional finite element model was successfully constructed using ANSYS and MIMICS software (model M1). The M1 was modified to simulate the foraminoplasty of different facet joint portions, with unilateral cylindrical excision (diameter = 0.75 cm) performed on the tip (model M2) and the base (model M3) of right L₅ superior facet elements along with surrounding capsular ligaments, respectively. Under the same loading conditions, the ROM and IDP of L₃/₄ and L₄/L₅ segments in states of forward flexion, backward extension, left lateral bending, right lateral bending, left axial rotation and right axial rotation were all compared.

Result

Compared with the intact model in backward extension, M2 increased the ROM of L_4/5 segment by 9.4% and IDP by 11.7%, while the ROM and IDP of M3 changed only slightly. In right axial rotation, M2 and M3 increased the ROM of L_4/5 segment by 17.9% and by 3.6%, respectively. In left axial rotation, M2 and M3 increased the ROM of L₄/L₅ segment by 7.14% and 3.6%, respectively. As for other states including forward flexion, left lateral bending, right lateral bending, the ROM and IDP were not significantly distinct between these two models. While focusing on L₃/L₄ segment, obviously changes in the ROM and IDP have not been presented and neither M2 nor M3 changed in any loading condition.

Conclusion

This study provides evidence that the base‐facet foraminoplasty of L₅ superior facet provided a higher segmental stability compared with the tip‐facet foraminoplasty in flexion and axial rotation. Meanwhile, it also shows the two types of foraminoplasty make few differences to the L_4/5 segmental biomechanics. Besides, it does not appear to impact the stability of L₃/L₄ in six states of forward flexion, backward extension, left lateral bending, right lateral bending, left axial rotation and right axial rotation when superior facet of L₅ was partially removed. These findings might be useful in understanding biomechanics of the lumbar spine after foraminoplasty performed on different portions of the facet, thus providing endoscopic surgeons a better reference for operational approach to maintain the function and mobility of the spine.

Biomechanical Effect of L4 -L5 Intervertebral Disc Degeneration on the Lower Lumbar Spine: A Finite Element Study

OBJECTIVE

To ascertain the biomechanical effects of a degenerated L₄ -L₅ segment on the lower lumbar spine through a comprehensive simulation of disc degeneration.

METHODS

A three-dimensional nonlinear finite element model of a normal L₃ -S₁ lumbar spine was constructed and validated. This normal model was then modified such that three degenerated models with different degrees of degeneration (mild, moderate, or severe) at the L₄ -L₅ level were constructed. While experiencing a follower compressive load (500 N), hybrid moment loads were applied to all models to determine range of motion (ROM), intradiscal pressure (IDP), maximum von Mises stress in the annulus, maximum shear stress in the annulus, and facet joint force.

RESULTS

As the degree of disc degeneration increased, the ROM of the L₄ -L₅ degenerated segment declined dramatically in all postures (flexion: 5.79°-1.91°; extension: 5.53°-2.62°; right lateral bending: 4.47°-1.46°; left lateral bending: 4.86°-1.61°; right axial rotation: 2.69°-0.74°; left axial rotation: 2.69°-0.74°), while the ROM in adjacent segments increased (1.88°-8.19°). The largest percent decrease in motion of the L₄ -L₅ segment due to disc degeneration was in right axial rotation (75%), left axial rotation (69%), flexion (67%), right lateral bending (67%), left lateral bending right (67%), and extension (53%). The change in the trend of the IDP was the same as that of the ROM. Specifically, the IDP decreased (flexion: 0.592-0.09 MPa; extension: 0.678-0.334 MPa; right lateral bending: 0.498-0.205 MPa; left lateral bending: 0.523-0.272 MPa; right axial rotation: 0.535-0.246 MPa; left axial rotation: 0.53-0.266 MPa) in the L₄ -L₅ segment, while the IDP in adjacent segments increased (0.511-0.789 MPa). The maximum von Mises stress and maximum shear stress of the annulus in whole lumbar spine segments increased (L₄ -L₅ segment: 0.413-2.626 MPa and 0.412-2.783 MPa, respectively; adjacent segment of L₄ -L₅ : 0.356-1.493 MPa and 0.359-1.718 MPa, respectively) as degeneration of the disc progressively increased. There was no apparent regularity in facet joint force in the degenerated segment as the degree of disc degeneration increased. Nevertheless, facet joint forces in adjacent healthy segments increased as the degree of disc degeneration increased (extension: 49.7-295.3 N; lateral bending: 3.5-171.2 N; axial rotation: 140.2-258.8 N).

CONCLUSION

Degenerated discs caused changes in the motion and loading pattern of the degenerated segments and adjacent normal segments. The abnormal load and motion in the degenerated models risked accelerating degeneration in the adjacent normal segments. In addition, accurate simulation of degenerated facet joints is essential for predicting changes in facet joint loads following disc degeneration.

Effect of Spiral Nucleus Implant Parameters on the Compressive Biomechanics of Lumbar Intervertebral Disc

OBJECTIVE

To determine the effect of spiral nucleus implant parameters on the biomechanical behavior of the lumbar intervertebral disc after nucleus replacement under compressive loading.

METHODS

A finite element (FE) model of nucleus replacement in the L4-5 intervertebral disc was constructed. The effects of a spiral implant parameters, such as elasticity, size, and friction property, on the biomechanical behavior of the disc under a compressive load were determined. The effect of an implant with a sharp edge on disc biomechanics was also examined. The stress distribution and contact pressure on the endplate and AF, axial stiffness of disc, and annular bulge of the nucleus replacement models were investigated.

RESULTS

Axial stiffness, annular bulge, and contact pressure were all insensitive to friction properties. Insertion of the spiral implant reversed the changes in the AF and endplates due to the removal of the nucleus. There was a positive correlation between axial stiffness and elasticity with implant size. Annular bulge was positively correlated with size but negatively correlated with elasticity. Compared with the base model, the implant with a sharp edge caused a decrease in disc axial stiffness but an increase in contact pressure on the AF in an annular bulge in the sagittal and coronal axis, respectively.

CONCLUSIONS

A spiral implant may provide similar biomechanical behavior as a normal disc during compressive loading, with an optimal modulus of approximately 7 MPa. The spiral implant should fully conform to the nucleus cavity during replacement for the best biomechanical results.

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

Background

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

Methods

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

Results

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

Conclusions

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

Facet joint hypertrophy is a misnomer: A retrospective study

One of the major causes of lumbar spinal canal stenosis (LSCS) has been considered facet joint hypertrophy (FJH). However, a previous study asserted that "FJH" is a misnomer because common facet joints are no smaller than degenerative facet joints; however, this hypothesis has not been effectively demonstrated. Therefore, in order to verify that FJH is a misnomer in patients with LSCS, we devised new morphological parameters that we called facet joint thickness (FJT) and facet joint cross-sectional area (FJA).We collected FJT and FJA data from 114 patients with LSCS. A total of 86 control subjects underwent lumbar magnetic resonance imaging (MRI) as part of routine medical examinations, and axial T2-weighted MRI images were obtained from all participants. We measured FJT by drawing a line along the facet area and then measuring the narrowest point at L4-L5. We measured FJA as the whole cross-sectional area of the facet joint at the stenotic L4-L5 level.The average FJT was 1.60 ± 0.36 mm in the control group and 1.11 ± 0.32 mm in the LSCS group. The average FJA was 14.46 ± 5.17 mm in the control group and 9.31 ± 3.47 mm in the LSCS group. Patients with LSCS had significantly lower FJTs (P < .001) and FJAs (P < .001).FJH, a misnomer, should be renamed facet joint area narrowing. Using this terminology would eliminate confusion in descriptions of the facet joint.

Dynamic Response of the Lumbar Spine to Whole-body Vibration Under a Compressive Follower Preload

STUDY DESIGN

A finite element study of dynamic response of the lumbar spine to whole-body vibration.

OBJECTIVE

The aim of this study was to develop and validate a finite element model for exploring the impact of whole-body vibration on the entire lumbar spine with a compressive follower preload applied.

SUMMARY OF BACKGROUND DATA

Several finite element studies have investigated the biodynamic characteristics of the human lumbar spine when exposed to whole-body vibration. However, very limited studies have been performed to quantitatively describe dynamic response in time domain of the entire lumbar spine to vibration loading under a compressive follower preload.

METHODS

A three-dimensional nonlinear finite element model of the human lumbar spine (L1-sacrum) subjected to the compressive follower preload was created. Transient dynamic analysis was conducted on the model to compute the spinal response to a sinusoidal vertical vibration load of ±40 N under a 400 N preload. The obtained dynamic response results at all spinal levels were collected and plotted as a function of time. As a comparison, the corresponding results for vertical static loads (-40 and 40 N) under the preload (400 N) were also computed.

RESULTS

Plots of the dynamic response at all levels showed a cyclic response with time, and their vibration amplitudes (peak-to-bottom variations) were markedly higher than the corresponding changing amplitudes of static load cases. The increasing effect of the vibration load reached 314.5%, 263.2%, 242.4%, and 232.7%, respectively, in axial displacement of vertebral center, disc bulge, intradiscal pressure, and annulus stress (von-Mises stress). In addition, increasing the compressive follower preload led to an increase in the dynamic response and a decrease in their vibration amplitudes.

CONCLUSION

This study may be useful to help quantify the effect of cyclic loading on the entire lumbar spine under physiologic compressive loading, and better understand vibration characteristics of the spine.

LEVEL OF EVIDENCE

5.

The Effects of Physiological Biomechanical Loading on Intradiscal Pressure and Annulus Stress in Lumbar Spine: A Finite Element Analysis

The present study was conducted to examine the effects of body weight on intradiscal pressure (IDP) and annulus stress of intervertebral discs at lumbar spine. Three-dimensional finite element model of osseoligamentous lumbar spine was developed subjected to follower load of 500 N, 800 N, and 1200 N which represent the loads for individuals who are normal and overweight with the pure moments at 7.5 Nm in flexion and extension motions. It was observed that the maximum IDP was 1.26 MPa at L1-L2 vertebral segment. However, the highest increment of IDP was found at L4-L5 segment where the IDP was increased to 30% in flexion and it was more severe at extension motion reaching to 80%. Furthermore, the maximum annulus stress also occurred at the L1-L2 segment with 3.9 MPa in extension motion. However, the highest increment was also found at L4-L5 where the annulus stress increased to 17% in extension motion. Based on these results, the increase of physiological loading could be an important factor to the increment of intradiscal pressure and annulus fibrosis stress at all intervertebral discs at the lumbar spine which may lead to early intervertebral disc damage.

Segmental variations in facet joint translations during in vivo lumbar extension

The lumbar facet joint (FJ) is often associated with pathogenesis in the spine, but quantification of normal FJ motion remains limited to in vitro studies or static imaging of non-functional poses. The purpose of this study was to quantify lumbar FJ kinematics in healthy individuals during functional activity with dynamic stereo radiography (DSX) imaging. Ten asymptomatic participants lifted three known weights starting from a trunk-flexed (∼75°) position to an upright position while being imaged within the DSX system. High resolution computed tomography (CT) scan-derived 3D models of their lumbar vertebrae (L2-S1) were registered to the biplane 2D radiographs using a markerless model-based tracking technique providing instantaneous 3D vertebral kinematics throughout the lifting tasks. Effects of segment level and weight lifted were assessed using mixed-effect repeated measures ANOVA. Superior-inferior (SI) translation dominated FJ translation, with L5S1 showing significantly less translation magnitudes (Median (Md) = 3.5 mm, p < 0.0001) than L2L3, L3L4, and L4L5 segments (Md = 5.9 mm, 6.3 mm and 6.6 mm respectively). Linear regression-based slopes of continuous facet translations revealed strong linearity for SI translation (r² > 0.94), reasonably high linearity for sideways sliding (Z-) (r² > 0.8), but much less linearity for facet gap change (X-) (r² ∼ 0.5). Caudal segments (L4-S1), particularly L5S1, displayed greater coupling compared to cranial (L2-L4) segments, revealing distinct differences overall in FJ translation trends at L5S1. No significant effect of weight lifted on FJ translations was detected. The study presents a hitherto unavailable and highly precise baseline dataset of facet translations measured during a functional, dynamic lifting task.

The Effect of Lumbar Disc Herniation on Spine Loading Characteristics during Trunk Flexion and Two Types of Picking Up Activities

The main purpose of this study was to investigate the compensatory response of the muscle activities of seventeen major muscle groups in the spinal region, intradiscal forces of the five lumbar motion segment units (MSUs), and facet forces acting on the ten lumbar facet joints in patients with lumbar disc herniation (LDH). Twenty-six healthy adults and seven LDH patients performed trunk flexion, ipsilateral picking up, and contralateral picking up in sequence. Eight optical markers were placed on the landmarks of the pelvis and spinal process. The coordinates of these markers were captured to drive a musculoskeletal model to calculate the muscle activities, intradiscal forces, and facet forces. The muscle activities of the majority of the seventeen major muscle groups were found increases in LDH patients. In addition, the LDH patients displayed larger compressive forces and anteroposterior forces on all the five lumbar MSUs and more lumbar facet inventions on most facet joints. These findings suggest that the LDH patients demonstrate compensatory increases in the most trunk muscle activities and all spinal loads. These negative compensatory responses increase the risk of the aggravation of disc herniation. Therefore, treatment should intervene as earlier as possible for the severe LDH patients.

The effect of screw tunnels on the biomechanical stability of vertebral body after pedicle screws removal: a finite element analysis

PURPOSE

Posterior reduction and pedicle screw fixation is a widely used procedure for thoracic and lumbar vertebrae fractures. Usually, the pedicle screws would be removed after the fracture healing and screw tunnels would be left. The aim of this study is to evaluate the effect of screw tunnels on the biomechanical stability of the lumbar vertebral body after pedicle screws removal by finite element analysis (FEA).

METHODS

First, the CT values of the screw tunnels wall in the fractured vertebral bodies were measured in patients whose pedicle screws were removed, and they were then compared with the values of vertebral cortical bone. Second, an adult patient was included and the CT images of the lumbar spine were harvested. Three dimensional finite element models of the L1 vertebra with unilateral or bilateral screw tunnels were created based on the CT images. Different compressive loads were vertically acted on the models. The maximum loads which the models sustained and the distribution of the force in the different parts of the models were recorded and compared with each other.

RESULTS

The CT values of the tunnels wall and vertebral cortical bone were 387.126±62.342 and 399.204±53.612, which were not statistically different (P=0.149). The models of three dimensional tetrahedral mesh finite element of normal lumbar 1 vertebra were established with good geometric similarity and realistic appearance. After given the compressive loads, the cortical bone was the first one to reach its ultimate stress. The maximum loads which the bilateral screw tunnels model, unilateral screw tunnel model, and normal vertebral model can sustain were 3.97 Mpa, 3.83 Mpa, and 3.78 Mpa, respectively. For the diameter of the screw tunnels, the model with a diameter of 6.5 mm could sustain the largest load. In addition, the stress distributing on the outside of the cortical bone gradually decreased as the thickness of the tunnel wall increased.

CONCLUSIONS

Based on the FEA, pedicle screw tunnels would not decrease the biomechanical stability and strength of the vertebral body. A large diameter of screw tunnel and thick tunnel wall were helpful for the biomechanical stability of the vertebral body.

Comparing the effects of different dynamic sitting strategies in wheelchair seating on lumbar-pelvic angle

BACKGROUND

Prolonged static sitting in a wheelchair is associated with an increased risk of lower back pain. The wheelchair seating system is a key factor of this risk because it affects spinal loading in the sitting position. In this study, 7 dynamic sitting strategies (DSSs) are examined: lumbar prominent dynamic sitting (LPDS), back reclined dynamic sitting (BRDS), femur upward dynamic sitting (FUDS), lumbar prominent with back reclined dynamic sitting (LBDS), lumbar prominent with femur upward dynamic sitting (LFDS), back reclined with femur upward dynamic sitting (BFDS), and lumbar prominent with back reclined with femur upward dynamic sitting (LBFDS). The objective of this study was to analyze the biomechanical effects of these sitting strategies on lumbar-pelvic angles.

METHODS

Twenty able-bodied participants were recruited for the study. All participants performed LPDS, BRDS, FUDS, LBDS, LFDS, BFDS, and LBFDS in a random order. All lumbar-pelvic angle parameters, including the static lumbar angle, static pelvic angle, lumbar range of motion, and pelvic range of motion were measured and compared.

RESULTS

Results show that LBDS and LBFDS enabled the most beneficial lumbar movements, although the difference between the 2 strategies was nonsignificant. BRDS and BFDS enabled the most beneficial pelvic movements, although the difference between the 2 strategies was nonsignificant. Among all the upright DSSs, LPDS and LFDS enabled the most beneficial lumbar and pelvic movements, although no significant difference was observed between these 2 strategies.

CONCLUSIONS

We identified the effects and differences among 7 DSSs on lumbar-pelvic angles. Wheelchair users can choose the most suitable DSS that meets their needs. These findings may serve as a reference for practicing physicians or wheelchair users to choose an appropriate dynamic wheelchair seating system.

TRIAL REGISTRATION

ISRCTN12389808 , 18th November 2016, retrospectively registered.