Subject-specific loads on the lumbar spine in detailed finite element models scaled geometrically and kinematic-driven by radiography images

Effect of asymptomatic intervertebral flexion patterns on lumbar disc pressure: A finite element analysis study

Movement patterns may be a factor for manipulating the lumbar load, although little information is yet available in the literature about the relationship between this variable and intervertebral disc pressure (IDP). A finite element model of the lumbar spine (49-year-old asymptomatic female) was used to simulate intervertebral movements (L2-L5) of 127 asymptomatic participants. The data from participants that at least completed a simulation of lumbar vertebral movement during the first 53% of a movement cycle (flexion phase) were used for further analyses. Then, for each vertebral angular motion curve with constant spatial peaks, different temporal patterns were simulated in two stages: (1) in lumbar pattern exchange (LPE), each vertebral angle was simulated by the corresponding vertebrae of other participants data; (2) in vertebral pattern exchange (VPE), vertebral angles were simulated by each other. The k-mean algorithm was used to cluster two groups of variables; peak and cumulative IDP, in both stages of simulations (i.e., LPE and VPE). In the second stage of the simulation (VPE), Kendall's tau was utilized to consider the relationship between different temporal patterns and IDPs for each individual lumbar level. Cluster analyses showed that the temporal movement pattern did not exhibit any effect on the peak IDP while the cumulative IDP changed significantly for some patterns. Earlier involvement in lumbar motion at any level led to higher IDP in the majority of simulations. There is therefore a possibility of manipulating lumbar IDP by changing the temporal pattern with the same ROM, in which optimal distribution of the loads among lumbar levels may be applied as preventive or treatment interventions. Evaluating load benefits, such as load, on biomechanically relevant lumbar levels, dynamically measured by quantitative fluoroscopy, may help inform interventional exercises.

Effect of low back pain on the kinetics and kinematics of the lumbar spine - a combined in vivo and in silico investigation

Lifting is a significant risk factor for low back pain (LBP). Different biomechanical factors including spinal loads, kinematics, and muscle electromyography (EMG) activities have previously been investigated during lifting activities in LBP patients and asymptomatic individuals to identify their association with LBP. However, the findings were contradictory and inconclusive. Accurate and subject-specific prediction of spinal loads is crucial for understanding, diagnosing, planning tailored treatments, and preventing recurrent pain in LBP patients. Therefore, the present study aimed to estimate the L5-S1 compressive and resultant shear loads in 19 healthy and 17 non-specific chronic LBP individuals during various static load-holding tasks (holding a 10 kg box at hip, chest, and head height) using full-body and personalized musculoskeletal models driven by subject-specific in vivo kinematic/kinetic, EMG, and physiological cross-sectional areas (PCSAs) data. These biomechanical characteristics were concurrently analyzed to identify potential differences between the two groups. Statistical analyses showed that LBP had almost no significant effect on the range of motion (trunk, lumbar, pelvis), PCSA, and EMG. There were no significant differences (p > 0.05) in the predicted L5-S1 loads. However, as the task became more demanding, by elevating the hand-load from hip to head, LBP patients experienced significant increases in both compressive (33 %, p = 0.00) and shear (25 %, p = 0.02) loads, while asymptomatic individuals showed significant increases only in compressive loads (30 %, p = 0.01). This suggests that engaging in more challenging activities could potentially magnify the effect of LBP on the biomechanical factors and increase their discrimination capacity between LBP and asymptomatic individuals.

The 6 degrees-of-freedom range of motion of the L1-S1 vertebrae in young and middle-aged asymptomatic people

STUDY DESIGN

Controlled laboratory study.

OBJECTIVE

To determine the 6 degrees of freedom of lumbar vertebra in vivo during different functional activities in young and middle-aged asymptomatic subjects.

METHODS

A total of 26 asymptomatic subjects (M/F, 15/11; age, 20-55 years) were recruited in this study. They were divided into two groups: young group (number: 14; age: 20-30 years old) and middle-aged group (number: 12; age: 45-55 years old). The lumbar segment of each subject was scanned by computed tomography for the construction of three-dimensional (3D) models of the vertebra from L1 to S1. The lumbar spine was imaged by using a dual fluoroscopic system when the subjects performed different trunk postures. The 3D models of vertebrae were matched to two fluoroscopic images simultaneously in software. The range of motion (ROM) of vertebrae in the young and middle-aged groups was compared by using multiway analysis of variance, respectively.

RESULTS

During the supine to the upright posture, vertebral rotation of the L1-S1 occurred mainly around the mediolateral axis (mean: 3.9 ± 2.9°). Along the mediolateral axis, vertebral translation was significantly lower at L1-2 (7.7 ± 2.4 mm) and L2-3 (8.0 ± 3.5 mm) than at L3-4 (1.6 ± 1.2 mm), L4-5 (3.3 ± 2.6 mm), and L5-S1 (2.6 ± 1.9 mm). At the L4-5 level, the young group had a higher rotational ROM than the middle-aged group around all three axes during left-right bending. Along the anteroposterior axis, the young group had a lower translational ROM at L4-5 than the middle-aged group during left-right bending (4.6 ± 3.3 vs. 7.6 ± 4.8 mm; P < 0.05). At L5-S1, the young group had a lower translational ROM than the middle-aged group during flexion-extension, left-right bending, and left-right torsion.

CONCLUSION

This study explored the lumbar vertebral ROM at L1-S1 during different functional postures in both young and middle-aged volunteers. There were higher coupled translations at L3-4 and L4-5 than at the upper lumbar segments during supine to upright. The vertebral rotation decreased with age. In addition, the older subjects had a higher anteroposterior translation at the L4-5 segment and higher mediolateral translation at the L5-S1 segment than the young group. These data might provide basic data to be compared with spinal pathology.

The influence of the viscoelastic property of polycarbonate urethane as an artificial disc core material under various physiological motions at the L4-L5 level

Intervertebral disc (IVD) degeneration is one of the musculoskeletal disorders due to the Degenerative Disc Disease (DDD), that cause low back pain (LBP) and leads to a reduced range of motion. Spinal fusion and arthroplasty are the other surgical procedures that could replace the disc affected by DDD against artificial disc replacement (ADR). This study aims to analyse the biomechanical behaviour of proposed core material as Polycarbonate Urethane (PCU) in the L4-L5 lumbar segment for ADR with Ti-6Al-4V and Co-28Cr-6M as endplate materials and compare it to the performance of an ultra-high molecular weight polyethylene (UHMWPE) core. Finite element methods have been approached to measure the overall stress distribution along with other physiological motions like Flexion (FLEX), Extension (EXT), Axial rotation (AR) and Lateral bending (LB), respectively. Preload of 450 N compressive load, 8 N-m for Flex, 6 N-m for EXT, 6 N-m for AR and 4 N-m for LB are applied. It could be concluded that Ti-6Al-4V - PCU and Co-28Cr-6M - PCU is the best composition for the ADR for the L4-L5 level.

Novel force-displacement control passive finite element models of the spine to simulate intact and pathological conditions; comparisons with traditional passive and detailed musculoskeletal models

Passive finite element (FE) models of the spine are commonly used to simulate intact and various pre- and postoperative pathological conditions. Being devoid of muscles, these traditional models are driven by simplistic loading scenarios, e.g., a constant moment and compressive follower load (FL) that do not properly mimic the complex in vivo loading condition under muscle exertions. We aim to develop novel passive FE models that are driven by more realistic yet simple loading scenarios, i.e., in vivo vertebral rotations and pathological-condition dependent FLs (estimated based on detailed musculoskeletal finite element (MS-FE) models). In these novel force-displacement control FE models, unlike the traditional passive FE models, FLs vary not only at different spine segments (T12-S1) but between intact, pre- and postoperative conditions. Intact, preoperative degenerated, and postoperative fused conditions at the L4-L5 segment for five static in vivo activities in upright and flexed postures were simulated by the traditional passive FE, novel force-displacement control FE, and gold-standard detailed MS-FE spine models. Our findings indicate that, when compared to the MS-FE models, the force-displacement control passive FE models could accurately predict the magnitude of disc compression force, intradiscal pressure, annulus maximal von Mises stress, and vector sum of all ligament forces at adjacent segments (L3-L4 and L5-S1) but failed to predict disc shear and facet joint forces. In this regard, the force-displacement control passive FE models were much more accurate than the traditional passive FE models. Clinical recommendations made based on traditional passive FE models should, therefore, be interpreted with caution.

Biomechanical Properties of Paraspinal Muscles Influence Spinal Loading—A Musculoskeletal Simulation Study

Paraspinal muscles are vital to the functioning of the spine. Changes in muscle physiological cross-sectional area significantly affect spinal loading, but the importance of other muscle biomechanical properties remains unclear. This study explored the changes in spinal loading due to variation in five muscle biomechanical properties: passive stiffness, slack sarcomere length (SSL), in situ sarcomere length, specific tension, and pennation angle. An enhanced version of a musculoskeletal simulation model of the thoracolumbar spine with 210 muscle fascicles was used for this study and its predictions were validated for several tasks and multiple postures. Ranges of physiologically realistic values were selected for all five muscle parameters and their influence on L4-L5 intradiscal pressure (IDP) was investigated in standing and 36° flexion. We observed large changes in IDP due to changes in passive stiffness, SSL, in situ sarcomere length, and specific tension, often with interesting interplays between the parameters. For example, for upright standing, a change in stiffness value from one tenth to 10 times the baseline value increased the IDP only by 91% for the baseline model but by 945% when SSL was 0.4 μm shorter. Shorter SSL values and higher stiffnesses led to the largest increases in IDP. More changes were evident in flexion, as sarcomere lengths were longer in that posture and thus the passive curve is more influential. Our results highlight the importance of the muscle force-length curve and the parameters associated with it and motivate further experimental studies on in vivo measurement of those properties.

Computational modelling of the scoliotic spine: A literature review

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

Stress analysis of intervertebral disc during occupational activities

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

A novel coupled musculoskeletal finite element model of the spine - Critical evaluation of trunk models in some tasks

Spine musculoskeletal (MS) models make simplifying assumptions on the intervertebral joint degrees-of-freedom (rotational and/or translational), representation (spherical or beam-like joints), and properties (linear or nonlinear). They also generally neglect the realistic structure of the joints with disc nuclei/annuli, facets, and ligaments. We aim to develop a novel MS model where trunk muscles are incorporated into a detailed finite element (FE) model of the ligamentous T12-S1 spine thus constructing a gold standard coupled MS-FE model. Model predictions are compared under some tasks with those of our earlier spherical joints, beam joints, and hybrid (uncoupled) MS-FE models. The coupled model predicted L4-L5 intradiscal pressures (R² ≅ 0.97, RMSE ≅ 0.27 MPa) and L1-S1 centers of rotation (CoRs) in agreement to in vivo data. Differences in model predictions grew at larger trunk flexion angles; at the peak (80°) flexion the coupled model predicted, compared to the hybrid model, much smaller global/local muscle forces (~38%), segmental (~44%) and disc (~22%) compression forces but larger segmental (~9%) and disc (~17%) shear loads, ligament forces at the lower lumbar levels (by up to 57%) and facet forces at all levels. The spherical/beam joints models predicted much greater muscle forces and segmental loads under larger flexion angles. Unlike the spherical joints model with fixed CoRs, the beam joints model predicted CoRs closer (RMSE = 2.3 mm in flexion tasks) to those of the coupled model. The coupled model offers a great potential for future studies towards improvement of surgical techniques, management of musculoskeletal injuries and subject-specific simulations.

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

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