Continuous dynamic spinal motion analysis

Musculoskeletal spine modeling in large patient cohorts: how morphological individualization affects lumbar load estimation

Introduction: Achieving an adequate level of detail is a crucial part of any modeling process. Thus, oversimplification of complex systems can lead to overestimation, underestimation, and general bias of effects, while elaborate models run the risk of losing validity due to the uncontrolled interaction of multiple influencing factors and error propagation. Methods: We used a validated pipeline for the automated generation of multi-body models of the trunk to create 279 models based on CT data from 93 patients to investigate how different degrees of individualization affect the observed effects of different morphological characteristics on lumbar loads. Specifically, individual parameters related to spinal morphology (thoracic kyphosis (TK), lumbar lordosis (LL), and torso height (TH)), as well as torso weight (TW) and distribution, were fully or partly considered in the respective models according to their degree of individualization, and the effect strengths of these parameters on spinal loading were compared between semi- and highly individualized models. T-distributed stochastic neighbor embedding (T-SNE) analysis was performed for overarching pattern recognition and multiple regression analyses to evaluate changes in occurring effects and significance. Results: We were able to identify significant effects (p < 0.05) of various morphological parameters on lumbar loads in models with different degrees of individualization. Torso weight and lumbar lordosis showed the strongest effects on compression (β ≈ 0.9) and anterior-posterior shear forces (β ≈ 0.7), respectively. We could further show that the effect strength of individual parameters tended to decrease if more individual characteristics were included in the models. Discussion: The induced variability due to model individualization could only partly be explained by simple morphological parameters. Our study shows that model simplification can lead to an emphasis on individual effects, which needs to be critically assessed with regard to in vivo complexity. At the same time, we demonstrated that individualized models representing a population-based cohort are still able to identify relevant influences on spinal loading while considering a variety of influencing factors and their interactions.

Evaluation of trunk muscle coactivation predictions in multi-body models

Musculoskeletal simulations with muscle optimization aim to minimize muscle effort, hence are considered unable to predict the activation of antagonistic muscles. However, activation of antagonistic muscles might be necessary to satisfy the dynamic equilibrium. This study aims to elucidate under which conditions coactivation can be predicted, to evaluate factors modulating it, and to compare the antagonistic activations predicted by the lumbar spine model with literature data. Simple 2D and 3D models, comprising of 2 or 3 rigid bodies, with simple or multi-joint muscles, were created to study conditions under which muscle coactivity is predicted. An existing musculoskeletal model of the lumbar spine developed in AnyBody was used to investigate the effects of modeling intra-abdominal pressure (IAP), linear/cubic and load/activity-based muscle recruitment criterion on predicted coactivation during forward flexion and lateral bending. The predicted antagonist activations were compared to reported EMG data. Muscle coactivity was predicted with simplified models when multi-joint muscles were present or the model was three-dimensional. During forward flexion and lateral bending, the coactivation ratio predicted by the model showed good agreement with experimental values. Predicted coactivation was negligibly influenced by IAP but substantially reduced with a force-based recruitment criterion. The conditions needed in multi-body models to predict coactivity are: three-dimensionality or multi-joint muscles, unless perfect antagonists. The antagonist activations are required to balance 3D moments but do not reflect other physiological phenomena, which might explain the discrepancies between model predictions and experimental data. Nevertheless, the findings confirm the ability of the multi-body trunk models to predict muscle coactivity and suggest their overall validity.

Biomechanical analysis of complications following T10-Pelvis spinal fusion: A population based computational study

Proximal junctional kyphosis (PJK) and proximal junctional failure (PJF) are challenging complications of long fusion constructs for the treatment of adult spinal deformity. The objective of this study is to understand the biomechanical stresses proximal to the upper instrumentation of a T10-pelvis fusion in a large patient cohort. The pre-fusion models were subject-specific thoracolumbar spine models that incorporate the height, weight, spine curvature, and muscle morphology of 250 individuals from the Framingham Heart Study Multidetector CT Study. To create post-fusion models, the subject-specific models were further modified to eliminate motion between the intervertebral joints from T10 to the pelvis. OpenSim analysis tools were used to calculate the medial lateral shear force, anterior posterior shear force, and compressive force on the T9 vertebra during the static postures. Differences between pre-fusion and post-fusion T9 biomechanics were consistent between increased segmental mobility and unchanged segmental mobility conditions. For all static postures, compression decreased (p < 0. 0005). Anterior-posterior shear force significantly increased (p < 0. 0005) during axial twist and significantly increased (p < 0. 0005) during trunk flexion. Medial lateral shear force significantly increased (p < 0. 0005) during axial twist. This computational study provided the first use of subject-specific models to investigate the biomechanics of long spinal fusions. Patients undergoing T10-Pelvis fusion were predicted to have increased shear forces and decreased compressive force at the T9 vertebra, independent of change in segmental mobility. The computational model shows potential for the investigation of spinal fusion biomechanics to reduce the risk of PJK or PJF.

An Enhanced Spine Model Validated for Simulating Dynamic Lifting Tasks in OpenSim

A fully articulated thoracolumbar spine model had been previously developed in OpenSim and had been extensively validated against experimental data during various static tasks. In the present study, we enhanced this detailed musculoskeletal model by adding the role of passive structures and adding kinematic constraints to make it suitable for dynamic tasks. We validated the spinal forces estimated by this enhanced model during nine dynamic lifting/lowering tasks. Moreover, we recently developed and evaluated five approaches in OpenSim to model the external loads applied to the hands during lifting/lowering tasks, and in the present study, we assessed which approach results in more accurate spinal forces. Regardless of the external load modeling approach, the maximum forces predicted by our enhanced spine model across all tasks, as well as the pattern of estimated spinal forces within each task, showed strong correlations (r-values and cross-correlation coefficients > 0.9) with experimental data. Given the biofidelity of our enhanced model, its accessibility via the open-source OpenSim software, and the extent to which this model has been validated, we recommend it for applications requiring estimation of spinal forces during lifting/lowering tasks using multibody-based models and inverse dynamic analyses.

Prediction of trunk muscle activation and spinal forces in adolescent idiopathic scoliosis during simulated trunk motion: A musculoskeletal modelling study

Due to lack of reference validation data, the common strategy in characterizing adolescent idiopathic scoliosis (AIS) by musculoskeletal modelling approach consists in adapting structure and parameters of validated body models of adult individuals with physiological alignments. Until now, only static postures have been replicated and investigated in AIS subjects. When aiming to simulate trunk motion, two critical factors need consideration: how distributing movement along the vertebral motion levels (lumbar spine rhythm), and if neglecting or accounting for the contribution of the stiffness of the motion segments (disc stiffness). The present study investigates the effect of three different lumbar spine rhythms and absence/presence of disc stiffness on trunk muscle imbalance in the lumbar region and on intervertebral lateral shear at different levels of the thoracolumbar/lumbar scoliotic curve, during simulated trunk motions in the three anatomical planes (flexion/extension, lateral bending, and axial rotation). A spine model with articulated ribcage previously developed in AnyBody software and adapted to replicate the spinal alignment in AIS subjects is employed. An existing dataset of 100 subjects with mild and moderate scoliosis is exploited. The results pointed out the significant impact of lumbar spine rhythm configuration and disc stiffness on changes in the evaluated outputs, as well as a relationship with scoliosis severity. Unfortunately, no optimal settings can be identified due to lack of reference validation data. According to that, extreme caution is recommended when aiming to adapt models of adult individuals with physiological alignments to adolescent subjects with scoliotic deformity.

Effect of obesity on spinal loads during load-reaching activities: A subject- and kinematics-specific musculoskeletal modeling approach

Obesity has been associated to increase the risk of low back disorders. Previous musculoskeletal models simulating the effect of body weight on intervertebral joint loads have assumed identical body postures for obese and normal-weight individuals during a given physical activity. Our recent kinematic-measurement studies, however, indicate that obese individuals adapt different body postures (segmental orientations) than normal-weight ones when performing load-reaching activities. The present study, therefore, used a subject- and kinematics-specific musculoskeletal modeling approach to compare spinal loads of nine normal-weight and nine obese individuals each performing twelve static two-handed load-reaching activities at different hand heights, anterior distances, and asymmetry angles (total of 12 tasks × 18 subjects = 216 model simulations). Each model incorporated personalized muscle architectures, body mass distributions, and full-body kinematics for each subject and task. Results indicated that even when accounting for subject-specific body kinematics obese individuals experienced significantly larger (by ∼38% in average) L5-S1 compression (2305 ± 468 N versus 1674 ± 337 N) and shear (508 ± 111 N versus 705 ± 150 N) loads during all reaching activities (p < 0.05 for all hand positions). This average difference of ∼38% was similar to the results obtained from previous modeling investigations that neglected kinematics differences between the two weight groups. Moreover, there was no significant interaction effect between body weight and hand position on the spinal loads; indicating that the effect of body weight on L5-S1 loads was not dependent on the position of hands. Postural differences alone appear, hence, ineffective in compensating the greater spinal loads that obese people experience during reaching activities.

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

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

Muscle-driven forward dynamic active hybrid model of the lumbosacral spine: combined FEM and multibody simulation

Most spine models belong to either the musculoskeletal multibody (MB) or finite element (FE) method. Recently, coupling of MB and FE models has increasingly been used to combine advantages of both methods. Active hybrid FE-MB models, still rarely used in spine research, avoid the interface and convergence problems associated with model coupling. They provide the inherent ability to account for the full interplay of passive and active mechanisms for spinal stability. In this paper, we developed and validated a novel muscle-driven forward dynamic active hybrid FE-MB model of the lumbosacral spine (LSS) in ArtiSynth to simultaneously calculate muscle activation patterns, vertebral movements, and internal mechanical loads. The model consisted of the rigid vertebrae L1-S1 interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, facet joints, and force actuators representing the muscles. Morphological muscle data were implemented via a semi-automated registration procedure. Four auxiliary bodies were utilized to describe non-linear muscle paths by wrapping and attaching the anterior abdominal muscles. This included an abdominal plate whose kinematics was optimized using motion capture data from upper body movements. Intra-abdominal pressure was calculated from the forces of the abdominal muscles compressing the abdominal cavity. For the muscle-driven approach, forward dynamics assisted data tracking was used to predict muscle activation patterns that generate spinal postures and balance the spine without prescribing accurate spinal kinematics. During calibration, the maximum specific muscle tension and spinal rhythms resulting from the model dynamics were evaluated. To validate the model, load cases were simulated from -10° extension to +30° flexion with weights up to 20 kg in both hands. The biomechanical model responses were compared with in vivo literature data of intradiscal pressures, intra-abdominal pressures, and muscle activities. The results demonstrated high agreement with this data and highlight the advantages of active hybrid modeling for the LSS. Overall, this new self-contained tool provides a robust and efficient estimation of LSS biomechanical responses under in vivo similar loads, for example, to improve pain treatment by spinal stabilization therapies.

The Variability of Lumbar Sequential Motion Patterns: Observational Study

BACKGROUND

Physiological motion of the lumbar spine is a topic of interest for musculoskeletal health care professionals since abnormal motion is believed to be related to lumbar complaints. Many researchers have described ranges of motion for the lumbar spine, but only few have mentioned specific motion patterns of each individual segment during flexion and extension, mostly comprising the sequence of segmental initiation in sagittal rotation. However, an adequate definition of physiological motion is still lacking. For the lower cervical spine, a consistent pattern of segmental contributions in a flexion-extension movement in young healthy individuals was described, resulting in a definition of physiological motion of the cervical spine.

OBJECTIVE

This study aimed to define the lumbar spines' physiological motion pattern by determining the sequence of segmental contribution in sagittal rotation of each vertebra during maximum flexion and extension in healthy male participants.

METHODS

Cinematographic recordings were performed twice in 11 healthy male participants, aged 18-25 years, without a history of spine problems, with a 2-week interval (time point T1 and T2). Image recognition software was used to identify specific patterns in the sequence of segmental contributions per individual by plotting segmental rotation of each individual segment against the cumulative rotation of segments L1 to S1. Intraindividual variability was determined by testing T1 against T2. Intraclass correlation coefficients were tested by reevaluation of 30 intervertebral sequences by a second researcher.

RESULTS

No consistent pattern was found when studying the graphs of the cinematographic recordings during flexion. A much more consistent pattern was found during extension, especially in the last phase. It consisted of a peak in rotation in L3L4, followed by a peak in L2L3, and finally, in L1L2. This pattern was present in 71% (15/21) of all recordings; 64% (7/11) of the participants had a consistent pattern at both time points. Sequence of segmental contribution was less consistent in the lumbar spine than the cervical spine, possibly caused by differences in facet orientation, intervertebral discs, overprojection of the pelvis, and muscle recruitment.

CONCLUSIONS

In 64% (7/11) of the recordings, a consistent motion pattern was found in the upper lumbar spine during the last phase of extension in asymptomatic young male participants. Physiological motion of the lumbar spine is a broad concept, influenced by multiple factors, which cannot be captured in a firm definition yet.

TRIAL REGISTRATION

ClinicalTrials.gov NCT03737227; https://clinicaltrials.gov/ct2/show/NCT03737227.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

RR2-10.2196/14741.

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

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULT

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

CONCLUSION

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

Multibody Models of the Thoracolumbar Spine: A Review on Applications, Limitations, and Challenges

Numerical models of the musculoskeletal system as investigative tools are an integral part of biomechanical and clinical research. While finite element modeling is primarily suitable for the examination of deformation states and internal stresses in flexible bodies, multibody modeling is based on the assumption of rigid bodies, that are connected via joints and flexible elements. This simplification allows the consideration of biomechanical systems from a holistic perspective and thus takes into account multiple influencing factors of mechanical loads. Being the source of major health issues worldwide, the human spine is subject to a variety of studies using these models to investigate and understand healthy and pathological biomechanics of the upper body. In this review, we summarize the current state-of-the-art literature on multibody models of the thoracolumbar spine and identify limitations and challenges related to current modeling approaches.

Lumbar spinal loads and lumbar muscle forces evaluation with various lumbar supports and backrest inclination angles in driving posture

PURPOSE

The low back pain of professional drivers could be linked to excessive lumbar load. This study aims at developing a musculoskeletal model to study the lumbar spinal loads and lumbar muscle forces of the human body in driving posture, so as to contribute to a better understanding of low back pain and to improve the design of vehicle seats.

METHODS

A standing musculoskeletal model, including limbs, head and neck, that can reflect several activities of daily living was established based on the Christophy spine model. The model was then validated by comparing the calculated lumbar loads and muscle forces to the experimental data in the previous studies. Referring to radiology studies, the musculoskeletal model was adjusted into different driving postures with several different lumbar supports (0, 2 and 4 cm) and inclinations of the backrest (from 23° to 33°, by 2° intervals). The lumbar biomechanical load with various lumbar supports and backrest inclination angles was calculated.

RESULTS

The results showed that the overall lumbar spinal load and lumbar muscle force with 4 cm lumbar support were reduced by 11.30 and 26.24%. The lumbar spinal loads and lumbar muscle forces increased first and then decreased with the increase in backrest inclination angles from 23° to 33°. The lumbar biomechanical load varied slightly with the backrest inclination angles from 29° to 33°.

CONCLUSIONS

There are two findings: (i) the lumbar spinal loads at the L3-L4, L4-L5 and L5-S1, and lumbar muscle forces decreased obviously with the 4 cm lumbar support, while the seat cushion inclination angle was set to 10°. (ii) The recommended backrest inclination angles are 29° to 33° with a 10° seat cushion to the horizontal, which can keep a low level of the lumbar spinal loads and lumbar muscle forces. This study could be used to explain the association between drivers' sitting posture and the lumbar load change, and provide a reference for the prevention of low back pain.

Evaluation of spinal force normalization techniques

Division normalization is commonly used in biomechanics studies to remove the effect of anthropometric differences (e.g., body weight) on kinetic variables, facilitating comparison across a population. In spine biomechanics, spinal forces are commonly divided by the body weight or the intervertebral load during a standing posture. However, it has been suggested that offset and power curve normalization are more appropriate than division normalization for normalizing kinetic variables such as ground reaction forces during walking and running. The present study investigated, for the first time, the effectiveness of four techniques for normalizing spinal forces to remove the effect of body weight. Spinal forces at all lumbar levels were estimated using a detailed OpenSim musculoskeletal model of the spine for 11 scaled models (50-100 kg) and during 13 trunk flexion tasks. Pearson correlations of raw and normalized forces against body weight were used to assess the effectiveness of each normalization technique. Body weight and standing division normalization could only successfully normalize L₄L₅ spinal forces in three tasks, and L₅S₁ loads in five and three tasks, respectively; however, offset and power curve normalization techniques were successful across all lumbar spine levels and all tasks. Offset normalization successfully removed the effect of body weight and maintained the influence of flexion angle on spinal forces. Thus, we recommend offset normalization to account for anthropometric differences in studies of spinal forces.

Validation of a Patient-Specific Musculoskeletal Model for Lumbar Load Estimation Generated by an Automated Pipeline From Whole Body CT

Background: Chronic back pain is a major health problem worldwide. Although its causes can be diverse, biomechanical factors leading to spinal degeneration are considered a central issue. Numerical biomechanical models can identify critical factors and, thus, help predict impending spinal degeneration. However, spinal biomechanics are subject to significant interindividual variations. Therefore, in order to achieve meaningful findings on potential pathologies, predictive models have to take into account individual characteristics. To make these highly individualized models suitable for systematic studies on spinal biomechanics and clinical practice, the automation of data processing and modeling itself is inevitable. The purpose of this study was to validate an automatically generated patient-specific musculoskeletal model of the spine simulating static loading tasks.

Methods: CT imaging data from two patients with non-degenerative spines were processed using an automated deep learning-based segmentation pipeline. In a semi-automated process with minimal user interaction, we generated patient-specific musculoskeletal models and simulated various static loading tasks. To validate the model, calculated vertebral loadings of the lumbar spine and muscle forces were compared with in vivo data from the literature. Finally, results from both models were compared to assess the potential of our process for interindividual analysis.

Results: Calculated vertebral loads and muscle activation overall stood in close correlation with data from the literature. Compression forces normalized to upright standing deviated by a maximum of 16% for flexion and 33% for lifting tasks. Interindividual comparison of compression, as well as lateral and anterior–posterior shear forces, could be linked plausibly to individual spinal alignment and bodyweight.

Conclusion: We developed a method to generate patient-specific musculoskeletal models of the lumbar spine. The models were able to calculate loads of the lumbar spine for static activities with respect to individual biomechanical properties, such as spinal alignment, bodyweight distribution, and ligament and muscle insertion points. The process is automated to a large extent, which makes it suitable for systematic investigation of spinal biomechanics in large datasets.

Kinematics of the Cervical Spine Under Healthy and Degenerative Conditions: A Systematic Review

Knowledge of spinal kinematics is essential for the diagnosis and management of spinal diseases. Distinguishing between physiological and pathological motion patterns can help diagnose these diseases, plan surgical interventions and improve relevant tools and software. During the last decades, numerous studies based on diverse methodologies attempted to elucidate spinal mobility in different planes of motion. The authors aimed to summarize and compare the evidence about cervical spine kinematics under healthy and degenerative conditions. This includes an illustrated description of the spectrum of physiological cervical spine kinematics, followed by a comparable presentation of kinematics of the degenerative cervical spine. Data was obtained through a systematic MEDLINE search including studies on angular/translational segmental motion contribution, range of motion, coupling and center of rotation. As far as the degenerative conditions are concerned, kinematic data regarding disc degeneration and spondylolisthesis were available. Although the majority of the studies identified repeating motion patterns for most motion planes, discrepancies associated with limited sample sizes and different imaging techniques and/or spine configurations, were noted. Among healthy/asymptomatic individuals, flexion extension (FE) and lateral bending (LB) are mainly facilitated by the subaxial cervical spine. C4-C5 and C5-C6 were the major FE contributors in the reported studies, exceeding the motion contribution of sub-adjacent segments. Axial rotation (AR) greatly depends on C1-C2. FE range of motion (ROM) is distributed between the atlantoaxial and subaxial segments, while AR ROM stems mainly from the former and LB ROM from the latter. In coupled motion rotation is quantitatively predominant over translation. Motion migrates caudally from C1-C2 and the center of rotation (COR) translocates anteriorly and superiorly for each successive subaxial segment. In degenerative settings, concurrent or subsequent lesions render the association between diseases and mobility alterations challenging. The affected segments seem to maintain translational and angular motion in early and moderate degeneration. However, the progression of degeneration restrains mobility, which seems to be maintained or compensated by adjacent non-affected segments. While the kinematics of the healthy cervical spine have been addressed by multiple studies, the entire nosological and kinematic spectrum of cervical spine degeneration is partially addressed. Large-scale in vivo studies can complement the existing evidence, cover the gaps and pave the way to technological and clinical breakthroughs.

CORRELATION BETWEEN THE MODIC I SIGN AND IMAGES OF VERTEBRAL INSTABILITY

ABSTRACT The lumbar region is highly affected by degenerative diseases and can be symptomatic as a result of inflammatory processes of the disc or segmental vertebral instability (SVI), according to concepts by White and Panjabi. Magnetic resonance imaging (MRI) and dynamic radiographs can be used to evaluate it. Objective: To determine whether images of Modic type I lesions are compatible with radiological vertebral instability. Methods: The MRIs and dynamic radiographs of 100 patients with a mean age of 46.94 years, 65 of whom were women and 35 of whom were men, were studied to evaluate the reaction of endplates according to Modic and SVI at levels L2-L3, L3-L4, L4-L5 and L5-S1. They were divided into 3 groups: A- Modic I and SVI; B- Modic I without SVI; and C- Without Modic I and SVI. Sex, age and type of disc were evaluated. Result: 46 Modic I lesions and 24 cases of SVI were recognized. There were 5 patients in Group A (OR 0.23, p=0.006), 38 in Group B; and 19 in Group C. No direct relationship between age and Modic I and/or SVI was observed. Pfirrmann grade 5 discs are 4.6 times more likely to present a Modic I signal. The most affected disc was L3-L4 for the Modic signal and L4-L5 with respect to the SVI, translational instability being more frequent. Conclusions: Modic type changes, identified by magnetic resonance, are clinically relevant in relation to low back pain, but their presence does not confirm either translational or angular instability. Level of evidence III; Cross-sectional cohort.

A Reference Database of Standardised Continuous Lumbar Intervertebral Motion Analysis for Conducting Patient-Specific Comparisons

Lumbar instability has long been thought of as the failure of lumbar vertebrae to maintain their normal patterns of displacement. However, it is unknown what these patterns consist of. Research using quantitative fluoroscopy (QF) has shown that continuous lumbar intervertebral patterns of rotational displacement can be reliably measured during standing flexion and return motion using standardised protocols and can be used to assess patients with suspected lumbar spine motion disorders. However, normative values are needed to make individualised comparisons. One hundred and thirty-one healthy asymptomatic participants were recruited and performed guided flexion and return motion by following the rotating arm of an upright motion frame. Fluoroscopic image acquisition at 15fps was performed and individual intervertebral levels from L2-3 to L5-S1 were tracked and analysed during separate outward flexion and return phases. Results were presented as proportional intervertebral motion representing these phases using continuous means and 95%CIs, followed by verification of the differences between levels using Statistical Parametric Mapping (SPM). A secondary analysis of 8 control participants matched to 8 patients with chronic, non-specific low back pain (CNSLBP) was performed for comparison. One hundred and twenty-seven asymptomatic participants’ data were analysed. Their ages ranged from 18 to 70 years (mean 38.6) with mean body mass index 23.8 kg/m² 48.8% were female. Both the flexion and return phases for each level evidenced continuous change in mean proportional motion share, with narrow confidence intervals, highly significant differences and discrete motion paths between levels as confirmed by SPM. Patients in the secondary analysis evidenced significantly less L5-S1 motion than controls (p < 0.05). A reference database of spinal displacement patterns during lumbar (L2-S1) intersegmental flexion and return motion using a standardised motion protocol using fluoroscopy is presented. Spinal displacement patterns in asymptomatic individuals were found to be distinctive and consistent for each intervertebral level, and to continuously change during bending and return. This database may be used to allow continuous intervertebral kinematics to drive dynamic models of joint and muscular forces as well as reference values against which to make patient-specific comparisons in suspected cases of lumbar spine motion disorders.

Comparing the risk of low-back injury using model-based optimization: Improved technique versus exoskeleton assistance

Although wearable robotic systems are designed to reduce the risk of low-back injury, it is unclear how effective assistance is, compared to improvements in lifting technique. We use a two-factor block study design to simulate how effective exoskeleton assistance and technical improvements are at reducing the risk of low-back injury when compared to a typical adult lifting a box. The effects of assistance are examined by simulating two different models: a model of just the human participant, and a model of the human participant wearing the SPEXOR exoskeleton. The effects of lifting technique are investigated by formulating two different types of optimal control problems: a least-squares problem which tracks the human participant's lifting technique, and a minimization problem where the model is free to use a different movement. Different lifting techniques are considered using three different cost functions related to risk factors for low-back injury: cumulative low-back load (CLBL), peak low-back load (PLBL), and a combination of both CLBL and PLBL (HYB). The results of our simulations indicate that an exoskeleton alone can make modest reductions in both CLBL and PLBL. In contrast, technical improvements alone are effective at reducing CLBL, but not PLBL. The largest reductions in both CLBL and PLBL occur when both an exoskeleton and technical improvements are used. While all three of the lifting technique cost functions reduce both CLBL and PLBL, the HYB cost function offers the most balanced reduction in both CLBL and PLBL.

The Influence of Kinematic Constraints on Model Performance During Inverse Kinematics Analysis of the Thoracolumbar Spine

Motion analysis is increasingly applied to spine musculoskeletal models using kinematic constraints to estimate individual intervertebral joint movements, which cannot be directly measured from the skin surface markers. Traditionally, kinematic constraints have allowed a single spinal degree of freedom (DOF) in each direction, and there has been little examination of how different kinematic constraints affect evaluations of spine motion. Thus, the objective of this study was to evaluate the performance of different kinematic constraints for inverse kinematics analysis. We collected motion analysis marker data in seven healthy participants (4F, 3M, aged 27–67) during flexion–extension, lateral bending, and axial rotation tasks. Inverse kinematics analyses were performed on subject-specific models with 17 thoracolumbar joints allowing 51 rotational DOF (51DOF) and corresponding models including seven sets of kinematic constraints that limited spine motion from 3 to 9DOF. Outcomes included: (1) root mean square (RMS) error of spine markers (measured vs. model); (2) lag-one autocorrelation coefficients to assess smoothness of angular motions; (3) maximum range of motion (ROM) of intervertebral joints in three directions of motion (FE, LB, AR) to assess whether they are physiologically reasonable; and (4) segmental spine angles in static ROM trials. We found that RMS error of spine markers was higher with constraints than without (p < 0.0001) but did not notably improve kinematic constraints above 6DOF. Compared to segmental angles calculated directly from spine markers, models with kinematic constraints had moderate to good intraclass correlation coefficients (ICCs) for flexion–extension and lateral bending, though weak to moderate ICCs for axial rotation. Adding more DOF to kinematic constraints did not improve performance in matching segmental angles. Kinematic constraints with 4–6DOF produced similar levels of smoothness across all tasks and generally improved smoothness compared to 9DOF or unconstrained (51DOF) models. Our results also revealed that the maximum joint ROMs predicted using 4–6DOF constraints were largely within physiologically acceptable ranges throughout the spine and in all directions of motions. We conclude that a kinematic constraint with 5DOF can produce smooth spine motions with physiologically reasonable joint ROMs and relatively low marker error.

An open-source musculoskeletal model of the lumbar spine and lower limbs: a validation for movements of the lumbar spine

Musculoskeletal models of the lumbar spine have been developed with varying levels of detail for a wide range of clinical applications. Providing consistency is ensured throughout the modelling approach, these models can be combined with other computational models and be used in predictive modelling studies to investigate bone health deterioration and the associated fracture risk. To provide precise physiological loading conditions for such predictive modelling studies, a new full-body musculoskeletal model including a detailed and consistent representation of the lower limbs and the lumbar spine was developed. The model was assessed against in vivo measurements from the literature for a range of spine movements representative of daily living activities. Comparison between model estimations and electromyography recordings was also made for a range of lifting tasks. This new musculoskeletal model will provide a comprehensive physiological mechanical environment for future predictive finite element modelling studies on bone structural adaptation. It is freely available on https://simtk.org/projects/llsm/.

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.

A novel method for prediction of postoperative global sagittal alignment based on full-body musculoskeletal modeling and posture optimization

Associations between spinal sagittal balance and pain and disability are well documented. Reciprocal changes after spinal surgery might be critical for the outcomes, but assessing their extent remains a challenge. This paper proposes a method for predicting full-body sagittal alignment including reciprocal changes in response to spinal fusion, based on musculoskeletal modeling and inverse-inverse dynamics approach. An established body model (AnyBody) was used, with fused segments modeled as rigid. Posture was optimized based on muscle expenditure minimization, following the concept of the cone of economy. The data of adult spinal fusion patients were obtained retrospectively from an ongoing clinical study. Patient spino-pelvic alignment, body weight and height, age- and pathology-related muscle deterioration, and underwent treatment details were represented in the model. Predicted postural changes were compared to follow-up radiographs to evaluate method validity. Twenty-one cases were analyzed in this preliminary study (age range = 48-74; number of fused segments 1-14). The model predictions correlated well with the radiographic measures at follow-up: TPA, r = 0.83; ΔPILL, r = 0.90; LL, r = 0.90; TK, r = 0.77. The model demonstrated high accuracy in predicting sagittal imbalance (positive predictive value = 1.00, negative predictive value = 0.75). The presented method for patient- and treatment-specific postoperative posture prediction can be used to guide preoperative planning of spinal fusion, but more extensive validation is needed.

Lumbar Intervertebral Motion Analysis During Flexion and Extension Cinematographic Recordings in Healthy Male Participants: Protocol

Background

Physiological motion of the lumbar spine is a subject of interest for musculoskeletal health care professionals, as abnormal motion is believed to be related to lumbar conditions and complaints. Many researchers have described ranges of motion for the lumbar spine, but only a few have mentioned specific motion patterns of each individual segment during flexion and extension. These motion patterns mostly comprise the sequence of segmental initiation in sagittal rotation. However, an adequate definition of physiological motion of the lumbar spine is still lacking. The reason for this is the reporting of different ranges of motion and sequences of segmental initiation in previous studies. Furthermore, due to insufficient fields of view, none of these papers have reported on maximum flexion and extension motion patterns of L1 to S1. In the lower cervical spine, a consistent pattern of segmental contributions was recently described. In order to understand physiological motion of the lumbar spine, it is necessary to systematically study motion patterns, including the sequence of segmental contribution, of vertebrae L1 to S1 in healthy individuals during maximum flexion and extension.

Objective

This study aims to define the lumbar spines’ physiological motion pattern of vertebrae L1, L2, L3, L4, L5, and S1 by determining the sequence of segmental contribution and the sequence of segmental initiation of motion in sagittal rotation of each vertebra during maximum flexion and extension. The secondary endpoint will be exploring the possibility of analyzing the intervertebral horizontal and vertical translation of each vertebra during maximum flexion and extension.

Methods

Cinematographic recordings will be performed on 11 healthy male participants, aged 18-25 years, without a history of spine problems. Cinematographic flexion and extension recordings will be made at two time points with a minimum 2-week interval in between.

Results

The study has been approved by the local institutional medical ethical committee (Medical Research Ethics Committee of Zuyderland and Zuyd University of Applied Sciences) on September 24, 2018. Inclusion of participants will be completed in 2020.

Conclusions

If successful, these physiological motion patterns can be compared with motion patterns of patients with lumbar conditions before or after surgery. Ultimately, researchers may be able to determine differences in biomechanics that can potentially be linked to physical complaints like low back pain.

Trial Registration

ClinicalTrials.gov NCT03737227; https://clinicaltrials.gov/ct2/show/NCT03737227

International Registered Report Identifier (IRRID)

DERR1-10.2196/14741

Twente Spine Model: A thorough investigation of the spinal loads in a complete and coherent musculoskeletal model of the human spine

Although in vivospinal loads have been previously measured, existing data are limited to certain lumbar and thoracic levels. A detailed investigation of spinal loads would assist with injury prevention and implant design but is unavailable. In this study, we developed a complete and coherent musculoskeletal model of the entire human spine and studied the intervertebral disc compression forces for physiological movements on three anatomical planes. This model incorporates the individual vertebrae at the cervical, thoracic, and lumbar regions, a flexible ribcage, and complete muscle anatomy. Intradiscal pressures were estimated from predicted compressive forces, and these were generally in close agreement with previously measured data. We found that compressive forces at the trunk discs increased during trunk lateral bending and axial rotation of the trunk. During flexion, compressive forces increased in the thoracolumbar and lumbar regions and slightly decreased at the middle thoracic discs. In extension, the forces generally decreased at the thoracolumbar and lumbar discs whereas they slightly increased at the upper and middle thoracic discs. Furthermore, similar to a previous biomechanical model of the cervical spine, our model predicted increased compression forces in neck flexion, lateral bending, and axial rotation, and decreased forces in neck extension.

Have Studies that Measure Lumbar Kinematics and Muscle Activity Concurrently during Sagittal Bending Improved Understanding of Spinal Stability and Sub-System Interactions? A Systematic Review

In order to improve understanding of the complex interactions between spinal sub-systems (i.e., the passive (ligaments, discs, fascia and bones), the active (muscles and tendons) and the neural control systems), it is necessary to take a dynamic approach that incorporates the measurement of multiple systems concurrently. There are currently no reviews of studies that have investigated dynamic sagittal bending movements using a combination of electromyography (EMG) and lumbar kinematic measurements. As such it is not clear how understanding of spinal stability concepts has advanced with regards to this functional movement of the spine. The primary aim of this review was therefore to evaluate how such studies have contributed to improved understanding of lumbar spinal stability mechanisms. PubMed and Cochrane databases were searched using combinations of the keywords related to spinal stability and sagittal bending tasks, using strict inclusion and exclusion criteria and adhering to PRISMA guidelines. Whilst examples of the interactions between the passive and active sub-systems were shown, typically small sample sizes meant that results were not generalizable. The majority of studies used regional kinematic measurements, and whilst this was appropriate in terms of individual study aims, the studies could not provide insight into sub-system interaction at the level of the spinal motion segment. In addition, the heterogeneity in methodologies made comparison between studies difficult. The review suggests that since Panjabi’s seminal spinal control papers, only limited advancement in the understanding of these theories has been provided by the studies under review, particularly at an inter-segmental level. This lack of progression indicates a requirement for new research approaches that incorporate multiple system measurements at a motion segment level.

Motion analysis of the cervical spine during extension and flexion: Reliability of the vertebral marking procedure

Cervical spine motion analysis using videofluoroscopy is currently a technique without a gold standard. We demonstrate the reliability of a rigid and reliable analysis methodology for cervical motion using videofluoroscopic images, representing the entire range of motion during flexion and extension, from the neutral position to the end-range in the sagittal plane. Two researchers with radiography and vertebral marking expertise, and two inexperienced researchers with 10 hours of training manually marked anatomical structures on fluoroscopic images in a procedure designed to control for vertebral rotation around the mid-plane axis. The average marking error across examiners and images was -0.12^∘ (standard deviation: 0.88°), and the intraexaminer error ranged from -1.00^∘ to 1.61° (standard deviation range: 0.27°-1.19°). Our method demonstrated lower errors compared to the higher resolution X-ray studies, and proved that vertebral marking can be performed by persons with no experience in radiographic image analysis.

Digital tracking algorithm reveals the influence of structural irregularities on joint movements in the human cervical spine

BACKGROUND

Disc height loss and osteophytes change the local mechanical environment in the spine; while previous research has examined kinematic dysfunction under degenerative change, none has looked at the influence of disc height loss and osteophytes throughout movement.

METHODS

Twenty patients with pain related to the head, neck or shoulders were imaged via videofluoroscopy as they underwent sagittal-plane flexion and extension. A clinician graded disc height loss and osteophytes as "severe/moderate", "mild", or "none". A novel tracking algorithm quantified motions of each vertebra. This information was used to calculate intervertebral angular and shear displacements. The digital algorithm made it practical to track individual vertebrae in multiple patients through hundreds of images without bias.

FINDINGS

Cases without height loss/osteophytes had a consistent increase in intervertebral angular displacement from C2/C3 to C5/C6, like that of healthy individuals, and mild height losses did not produce aberrations that were systematic or necessarily discernable. However, joints with moderate to severe disc height loss and osteophytes exhibited reduced range of motion compared to adjacent unaffected joints in that patient and corresponding joints in patients without structural irregularities.

INTERPRETATION

Digitally-obtained motion histories of individual joints allowed anatomical joint changes to be linked with changes in joint movement patterns. Specifically, disc height loss and osteophytes were found to influence cervical spine movement in the sagittal plane, reducing angular motions at affected joints by approximately 10% between those with and without height loss and osteophytes. Further, these joint changes were associated with perturbed intervertebral angular and shear movements.

Continuous lumbar spine rhythms during level walking, stair climbing and trunk flexion in people with and without lumbar disc herniation

Low back pain(LBP) is one of the most prevalent diseases afflicting people today. Abnormal musculoskeletal loadings during activities of daily living (ADLs) have been deemed to be associated with spine rhythm. But no studies have reported abnormal continuous spine rhythms during ADLs in LBP patients. Therefore, the objective of this study was to investigate the continuous lumbar spine rhythms and their difference between people with and without lumbar disc herniation (LDH). Twenty-six healthy people and seven patients with LDH were recruited in this study. They performed level walking, stair climbing, and trunk flexion. Active optical markers placed on the landmark of the spinous process and pelvis were captured using motion analysis system to drive a musculoskeletal model to calculate the continuous lumbar spine rhythms. It was found that the lumbar spine rhythm was roughly constant throughout the analyzed cycle in both healthy people and LDH patients during trunk flexion. LDH patients displayed fluctuant lumbar spine rhythms during level walking and stair climbing and significantly higher segmental contributions of the lumbar segments in the lower lumbar region during stair climbing and trunk flexion. In conclusion, there were different compensatory responses to LDH in the continuous lumbar spine rhythms during different ADLs. This study provides a new insight into the abnormal spinal motion in LDH patients.

Regarding loads after spinal fusion, every level should be seen separately: a musculoskeletal analysis

INTRODUCTION

The number of spinal fusion surgeries is steadily increasing and biomechanical consequences are still in debate. The aim of this study is to provide biomechanical insights into the sagittal balance of the spine and to compare spinal load before and after spinal fusion.

METHOD

The joint reaction forces of 52 patients were analyzed in proximo-distal and antero-posterior direction from the levels T12-L1 to L5-S1 using musculoskeletal simulations.

RESULTS

In 104 simulations, pre-surgical forces were equal to post-surgical. The levels L4-L5 and T12-L1, however, showed increased spinal forces compression forces with higher sagittal displacement. Improved restauration of sagittal balance was accompanied by lower spinal load. AP shear stress, interestingly decreased with sagittal imbalance.

CONCLUSION

Imbalanced spines have a risk of increased compression forces at Th12-L1. L4-L5 always has increased spinal loads. These slides can be retrieved under Electronic Supplementary Material.

Biomechanical effects of hybrid stabilization on the risk of proximal adjacent-segment degeneration following lumbar spinal fusion using an interspinous device or a pedicle screw–based dynamic fixator

OBJECTIVE

Pedicle screw-rod–based hybrid stabilization (PH) and interspinous device–based hybrid stabilization (IH) have been proposed to prevent adjacent-segment degeneration (ASD) and their effectiveness has been reported. However, a comparative study based on sound biomechanical proof has not yet been reported. The aim of this study was to compare the biomechanical effects of IH and PH on the transition and adjacent segments.

METHODS

A validated finite element model of the normal lumbosacral spine was used. Based on the normal model, a rigid fusion model was immobilized at the L4–5 level by a rigid fixator. The DIAM or NFlex model was added on the L3–4 segment of the fusion model to construct the IH and PH models, respectively. The developed models simulated 4 different loading directions using the hybrid loading protocol.

RESULTS

Compared with the intact case, fusion on L4–5 produced 18.8%, 9.3%, 11.7%, and 13.7% increments in motion at L3–4 under flexion, extension, lateral bending, and axial rotation, respectively. Additional instrumentation at L3–4 (transition segment) in hybrid models reduced motion changes at this level. The IH model showed 8.4%, −33.9%, 6.9%, and 2.0% change in motion at the segment, whereas the PH model showed −30.4%, −26.7%, −23.0%, and 12.9%. At L2–3 (adjacent segment), the PH model showed 14.3%, 3.4%, 15.0%, and 0.8% of motion increment compared with the motion in the IH model. Both hybrid models showed decreased intradiscal pressure (IDP) at the transition segment compared with the fusion model, but the pressure at L2–3 (adjacent segment) increased in all loading directions except under extension.

CONCLUSIONS

Both IH and PH models limited excessive motion and IDP at the transition segment compared with the fusion model. At the segment adjacent to the transition level, PH induced higher stress than IH model. Such differences may eventually influence the likelihood of ASD.

Is L5-S1 motion segment different from the rest? A radiographic kinematic assessment of 72 patients with chronic low back pain

PURPOSE

The relationship between biomechanical instability and degenerative changes in the lumbar spine in chronic low back pain (CLBP) patients remains controversial. The main objective of this retrospective radiographical study was to evaluate changes in kinematics at different lumbar levels (in particular the L5-S1 level) with progressive grades of disc degeneration and facet joint osteoarthritis in CLBP patients.

METHODS

Using standing neutral and dynamic flexion/extension (Fx/Ex) radiographs of the lumbar spine, in vivo segmental kinematics at L1-L2 through L5-S1 were evaluated in 72 consecutive CLBP patients. Disc degeneration was quantified using changes in signal intensity and central disc height on mid-sagittal T2-weighted magnetic resonance (MR) scans. Additionally, the presence or absence of facet joint osteoarthritis was noted on T2-weighted axial MR scans.

RESULTS

Disc degeneration and facet joint osteoarthritis occurred independent of each other at the L5-S1 level (p = 0.188), but an association was observed between the two at L4-L5 (p < 0.001) and L3-L4 (p < 0.05) levels. In the absence of facet joint osteoarthritis, the L5-S1 segment showed a greater range of motion (ROM) in Ex (3.3° ± 3.6°) and a smaller ROM in Fx (0.6° ± 4.2°) compared with the upper lumbar levels (p < 0.05), but the differences diminished in the presence of it. In the absence of facet joint osteoarthritis, no change in L5-S1 kinematics was observed with progressive disc degeneration, but in its presence, restabilisation of the L5-S1 segment was observed between mild and severe disc degeneration states.

CONCLUSION

The L5-S1 motion segment exhibited unique degenerative and kinematic characteristics compared with the upper lumbar motion segments. Disc degeneration and facet joint osteoarthritis occurred independent of each other at the L5-S1 level, but not at the other lumbar levels. Severe disc degeneration in the presence of facet joint osteoarthritis biomechanically restabilised the L5-S1 motion segment.

Uneven intervertebral motion sharing is related to disc degeneration and is greater in patients with chronic, non-specific low back pain: an in vivo, cross-sectional cohort comparison of intervertebral dynamics using quantitative fluoroscopy

PURPOSE

Evidence of intervertebral mechanical markers in chronic, non-specific low back pain (CNSLBP) is lacking. This research used dynamic fluoroscopic studies to compare intervertebral angular motion sharing inequality and variability (MSI and MSV) during continuous lumbar motion in CNSLBP patients and controls. Passive recumbent and active standing protocols were used and the relationships of these variables to age and disc degeneration were assessed.

METHODS

Twenty patients with CNSLBP and 20 matched controls received quantitative fluoroscopic lumbar spine examinations using a standardised protocol for data collection and image analysis. Composite disc degeneration (CDD) scores comprising the sum of Kellgren and Lawrence grades from L2-S1 were obtained. Indices of intervertebral motion sharing inequality (MSI) and variability (MSV) were derived and expressed in units of proportion of lumbar range of motion from outward and return motion sequences during lying (passive) and standing (active) lumbar bending and compared between patients and controls. Relationships between MSI, MSV, age and CDD were assessed by linear correlation.

RESULTS

MSI was significantly greater in the patients throughout the intervertebral motion sequences of recumbent flexion (0.29 vs. 0.22, p = 0.02) and when flexion, extension, left and right motion were combined to give a composite measure (1.40 vs. 0.92, p = 0.04). MSI correlated substantially with age (R = 0.85, p = 0.004) and CDD (R = 0.70, p = 0.03) in lying passive investigations in patients and not in controls. There were also substantial correlations between MSV and age (R = 0.77, p = 0.01) and CDD (R = 0.85, p = 0.004) in standing flexion in patients and not in controls.

CONCLUSION

Greater inequality and variability of motion sharing was found in patients with CNSLBP than in controls, confirming previous studies and suggesting a biomechanical marker for the disorder at intervertebral level. The relationship between disc degeneration and MSI was augmented in patients, but not in controls during passive motion and similarly for MSV during active motion, suggesting links between in vivo disc mechanics and pain generation.

Validation of the AnyBody full body musculoskeletal model in computing lumbar spine loads at L4L5 level

In the panorama of available musculoskeletal modeling software, AnyBody software is a commercial tool that provides a full body musculoskeletal model which is increasingly exploited by numerous researchers worldwide. In this regard, model validation becomes essential to guarantee the suitability of the model in representing the simulated system. When focusing on lumbar spine, the previous works aimed at validating the AnyBody model in computing the intervertebral loads held several limitations, and a comprehensive validation is to be considered as lacking. The present study was aimed at extensively validating the suitability of the AnyBody model in computing lumbar spine loads at L4L5 level. The intersegmental loads were calculated during twelve specific exercise tasks designed to accurately replicate the conditions during which Wilke et al. (2001) measured in vivo the L4L5 intradiscal pressure. Motion capture data of one volunteer subject were acquired during the execution of the tasks and then imported into AnyBody to set model kinematics. Two different approaches in computing intradiscal pressure from the intersegmental load were evaluated. Lumbopelvic rhythm was compared with reference in vivo measurements to assess the accuracy of the lumbopelvic kinematics. Positive agreement was confirmed between the calculated pressures and the in vivo measurements, thus demonstrating the suitability of the AnyBody model. Specific caution needs to be taken only when considering postures characterized by large lateral displacements. Minor discrepancy was found assessing lumbopelvic rhythm. The present findings promote the AnyBody model as an appropriate tool to non-invasively evaluate the lumbar loads at L4L5 in physiological activities.

A videofluoroscopy-based tracking algorithm for quantifying the time course of human intervertebral displacements

The motions of individual intervertebral joints can affect spine motion, injury risk, deterioration, pain, treatment strategies, and clinical outcomes. Since standard kinematic methods do not provide precise time-course details about individual vertebrae and intervertebral motions, information that could be useful for scientific advancement and clinical assessment, we developed an iterative template matching algorithm to obtain this data from videofluoroscopy images. To assess the bias of our approach, vertebrae in an intact porcine spine were tracked and compared to the motions of high-contrast markers. To estimate precision under clinical conditions, motions of three human cervical spines were tracked independently ten times and vertebral and intervertebral motions associated with individual trials were compared to corresponding averages. Both tests produced errors in intervertebral angular and shear displacements no greater than 0.4° and 0.055 mm, respectively. When applied to two patient cases, aberrant intervertebral motions in the cervical spine were typically found to correlate with patient-specific anatomical features such as disc height loss and osteophytes. The case studies suggest that intervertebral kinematic time-course data could have value in clinical assessments, lead to broader understanding of how specific anatomical features influence joint motions, and in due course inform clinical treatments.

Incorporating Six Degree-of-Freedom Intervertebral Joint Stiffness in a Lumbar Spine Musculoskeletal Model-Method and Performance in Flexed Postures

Intervertebral translations and rotations are likely dependent on intervertebral stiffness properties. The objective of this study was to incorporate realistic intervertebral stiffnesses in a musculoskeletal model of the lumbar spine using a novel force-dependent kinematics approach, and examine the effects on vertebral compressive loading and intervertebral motions. Predicted vertebral loading and intervertebral motions were compared to previously reported in vivo measurements. Intervertebral joint reaction forces and motions were strongly affected by flexion stiffness, as well as force-motion coupling of the intervertebral stiffness. Better understanding of intervertebral stiffness and force-motion coupling could improve musculoskeletal modeling, implant design, and surgical planning.

Attainment rate as a surrogate indicator of the intervertebral neutral zone length in lateral bending: an in vitro proof of concept study

Background

Lumbar segmental instability is often considered to be a cause of chronic low back pain. However, defining its measurement has been largely limited to laboratory studies. These have characterised segmental stability as the intrinsic resistance of spine specimens to initial bending moments by quantifying the dynamic neutral zone. However these measurements have been impossible to obtain in vivo without invasive procedures, preventing the assessment of intervertebral stability in patients. Quantitative fluoroscopy (QF), measures the initial velocity of the attainment of intervertebral rotational motion in patients, which may to some extent be representative of the dynamic neutral zone. This study sought to explore the possible relationship between the dynamic neutral zone and intervertebral rotational attainment rate as measured with (QF) in an in vitro preparation. The purpose was to find out if further work into this concept is worth pursuing.

Method

This study used passive recumbent QF in a multi-segmental porcine model. This assessed the intrinsic intervertebral responses to a minimal coronal plane bending moment as measured with a digital force guage. Bending moments about each intervertebral joint were calculated and correlated with the rate at which global motion was attained at each intervertebral segment in the first 10° of global motion where the intervertebral joint was rotating.

Results

Unlike previous studies of single segment specimens, a neutral zone was found to exist during lateral bending. The initial attainment rates for left and right lateral flexion were comparable to previously published in vivo values for healthy controls. Substantial and highly significant levels of correlation between initial attainment rate and neutral zone were found for left (Rho = 0.75, P = 0.0002) and combined left-right bending (Rho = 0.72, P = 0.0001) and moderate ones for right alone (Rho = 0.55, P = 0.0012).

Conclusions

This study found good correlation between the initial intervertebral attainment rate and the dynamic neutral zone, thereby opening the possibility to detect segmental instability from clinical studies. However the results must be treated with caution. Further studies with multiple specimens and adding sagittal plane motion are warranted.

Development and Validation of a Musculoskeletal Model of the Fully Articulated Thoracolumbar Spine and Rib Cage

We developed and validated a fully articulated model of the thoracolumbar spine in opensim that includes the individual vertebrae, ribs, and sternum. To ensure trunk muscles in the model accurately represent muscles in vivo, we used a novel approach to adjust muscle cross-sectional area (CSA) and position using computed tomography (CT) scans of the trunk sampled from a community-based cohort. Model predictions of vertebral compressive loading and trunk muscle tension were highly correlated to previous in vivo measures of intradiscal pressure (IDP), vertebral loading from telemeterized implants and trunk muscle myoelectric activity recorded by electromyography (EMG).

Examination of a lumbar spine biomechanical model for assessing axial compression, shear, and bending moment using selected Olympic lifts

BACKGROUND/AIMS

Loading during concurrent bending and compression associated with deadlift, hang clean and hang snatch lifts carries the potential for injury to the intervertebral discs, muscles and ligaments. This study examined the capacity of a newly developed spinal model to compute shear and compressive forces, and bending moments in lumbar spine for each lift.

METHODS

Five male subjects participated in the study. The spine was modeled as a chain of rigid bodies (vertebrae) connected via the intervertebral discs. Each vertebral reference frame was centered in the center of mass of the vertebral body, and its principal directions were axial, anterior-posterior, and medial-lateral.

RESULTS

The results demonstrated the capacity of this spinal model to assess forces and bending moments at and about the lumbar vertebrae by showing the variations among these variables with different lifting techniques.

CONCLUSION

These results show the model's potential as a diagnostic tool.