Influence of passive elements on prediction of intradiscal pressure and muscle activation in lumbar musculoskeletal models

Is the disappearance of the cervical flexion-relaxation phenomenon associated with cervical degeneration in healthy people?

PURPOSE

This study aims to explore the differences in cervical degeneration between healthy people with and without cervical flexion-relaxation phenomenon (FRP) and to identify whether the disappearance of cervical FRP is related to cervical degeneration.

METHODS

According to the flexion relaxation ratio (FRR), healthy subjects were divided into the normal FRP group and the abnormal FRP group. Besides, MRI was used to evaluate the degeneration of the passive subsystem (vertebral body, intervertebral disc, cervical sagittal balance, etc.) and the active subsystem (deep flexors [DEs], deep extensors [DFs], and superficial extensors [SEs]). In addition, the correlation of the FRR with the cervical degeneration score, C2-7Cobb, Borden method, relative total cross-sectional area (rTCSA), relative functional cross-sectional area (rFCSA), and fatty infiltration ratio (FIR) was analyzed.

RESULTS

A total of 128 healthy subjects were divided into the normal FRP group (n=52, 40.63%) and the abnormal FRP group (n=76, 59.38%). There were significant differences between the normal FRP group and the abnormal FRP group in the cervical degeneration score (z=-6.819, P<0.001), C2-7Cobb (t=2.994, P=0.004), Borden method (t=2.811, P=0.006), and FIR of DEs (t=-4.322, P<0.001). The FRR was significantly correlated with the cervical degeneration score (r=-0.457, P<0.001), C2-7Cobb (r=0.228, P=0.010), Borden method (r=0.197, P=0.026), and FIR of DEs (r=-0.253, P=0.004).

CONCLUSION

The disappearance of cervical FRP is related to cervical degeneration. A new hypothesis mechanism for FRP is proposed. The cervical FRP test is an effective and noninvasive examination for the differential diagnosis of healthy people, people with potential NSNP, and patients with NSNP.

Variability of intervertebral joint stiffness between specimens and spine levels

Introduction: Musculoskeletal multibody models of the spine can be used to investigate the biomechanical behaviour of the spine. In this context, a correct characterisation of the passive mechanical properties of the intervertebral joint is crucial. The intervertebral joint stiffness, in particular, is typically derived from the literature, and the differences between individuals and spine levels are often disregarded. Methods: This study tested if an optimisation method of personalising the intervertebral joint stiffnesses was able to capture expected stiffness variation between specimens and between spine levels and if the variation between spine levels could be accurately captured using a generic scaling ratio. Multibody models of six T12 to sacrum spine specimens were created from computed tomography data. For each specimen, two models were created: one with uniform stiffnesses across spine levels, and one accounting for level dependency. Three loading conditions were simulated. The initial stiffness values were optimised to minimize the kinematic error. Results: There was a range of optimised stiffnesses across the specimens and the models with level dependent stiffnesses were less accurate than the models without. Using an optimised stiffness substantially reduced prediction errors. Discussion: The optimisation captured the expected variation between specimens, and the prediction errors demonstrated the importance of accounting for level dependency. The inaccuracy of the predicted kinematics for the level-dependent models indicated that a generic scaling ratio is not a suitable method to account for the level dependency. The variation in the optimised stiffnesses for the different loading conditions indicates personalised stiffnesses should also be considered load-specific.

Skin marker-based subject-specific spinal alignment modeling: A feasibility study

Musculoskeletal models have the potential to improve diagnosis and optimize clinical treatment by predicting accurate outcomes on an individual basis. However, the subject-specific modeling of spinal alignment is often strongly simplified or is based on radiographic assessments, exposing subjects to unnecessary radiation. We therefore developed and introduced a novel skin marker-based approach for modeling subject-specific spinal alignment and evaluated its feasibility by comparing the predicted L1/L2 spinal loads during various functional activities with the loads predicted by the generically scaled models as well as with in vivo measured data obtained from the OrthoLoad database. Spinal loading simulations resulted in considerably higher compressive forces for both scaling approaches over all simulated activities, and AP shear forces that were closer or similar to the in vivo data for the subject-specific approach during upright standing activities and for the generic approach during activities that involved large flexions. These results underline the feasibility of the proposed method and associated workflow for inter- and intra-subject investigations using musculoskeletal simulations. When implemented into standard model scaling workflows, it is expected to improve the accuracy of muscle activity and joint loading simulations, which is crucial for investigations of treatment effects or pathology-dependent deviations.

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.

Prediction of the natural frequencies of different degrees of degenerated human lumbar segments L2-L3 using dynamic finite element analysis

BACKGROUND AND OBJECTIVE

Chronic exposure to resonant environment may cause more serious injuries to human lumbar spine than other vibrations. On the condition that the natural frequency of human lumbar spine is known, excitation frequency from an external vibration source can be optimized to keep away from the natural frequency and thus avoid lumbar resonance. Therefore, this study aimed to present an approach to predict the natural frequency of the human lumbar spine.

METHODS

Four poroelastic finite element models of human L2-L3 spinal motion segments with different degrees of degeneration were established. Dynamic finite element analyses of these models during 1 h of vibration were then conducted. The mechanical parameters of these models under vibrations at different excitation frequencies were predicted. The excitation frequencies that resulted in the greatest changes in the lumbar mechanical parameters were identified as the natural frequencies of the established L2-L3 spinal motion segments.

RESULTS

Simulation results showed that the natural frequencies of the healthy and mildly degenerated L2-L3 spinal motion segments, moderately degenerated L2-L3 spinal motion segments, and seriously degenerated L2-L3 spinal motion segments were in the range of 5-7, 3-5, and 1-3 Hz, respectively.

CONCLUSIONS

The predicted results indicated that the natural frequencies of the human L2-L3 spinal motion segments gradually decreased with the severity of degeneration. These phenomena may be related to changes in the lumbar structures and materials because of degeneration. This study provided a feasible method to predict the lumbar natural frequencies for different populations, which may be helpful in optimizing external vibration sources to avoid lumbar resonance.

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.

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.