A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

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.

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.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

Effects of Artificially Induced Breast Augmentation on the Electromyographic Activity of Neck and Trunk Muscles during Common Daily Movements

A female breast can be a potential source of musculoskeletal problems, especially if it is disproportionately large. The purpose of the present study was to examine the effect of artificially induced breast volume on the EMG activity of neck and trunk musculature during common everyday movements. The EMG activity of the sternocleidomastoid (SCM), the upper trapezius (UT), and the thoracic and lumbar erector spinae (TES, LES) were recorded during 45° trunk inclination from the upright standing and sitting postures (TIST45°, TISI45°) as well as during stand-to-sit and sit-to-stand (STSI, SIST) in 24 healthy females with minimal and ideal breast volume (M-NBV, I-NBV). All movements were performed before and after increasing M-NBV and I-NBV by 1.5-, 3.0-, 4.5-, and 6-times using silicone-gel implants. Significantly higher EMG activity for TES and LES were found at 6.0- and ≥4.5-times increase the I-NBV, respectively, compared to smaller breast volumes during TIST45°. EMG activity of UT was higher, and TES was lower in M-NBV females compared to I-NBV females in all movements but were significantly different only during SIST. The female breast can affect the activity of neck and trunk muscles only when its volume increases above a certain limit, potentially contributing to muscle dysfunction.

How Can Biomechanical Multibody Models of Scoliosis Be Accurate in Simulating Spine Movement Behavior While Neglecting the Changes of Spinal Length?

This paper studies how biomechanical multibody models of scoliosis can neglect the changes of spinal length and yet be accurate in reconstructing spinal columns. As these models with fixed length comprise rigid links interconnected by rotary joints, they resemble polygonal chains that approximate spine curves with a finite number of line segments. In mathematics, using more segments with shorter lengths can result in more accurate curve approximations. This raises the question of whether more accurate spine curve approximations by increasing the number of links/joints can yield more accurate spinal column reconstructions. For this, the accuracy of spine curve approximation was improved consistently by increasing the number of links/joints, and its effects on the accuracy of spinal column reconstruction were assessed. Positive correlation was found between the accuracy of spine reconstruction and curve approximation. It was shown that while increasing the accuracy of curve approximations, the representation of scoliosis concavity and its side-to-side deviations were improved. Moreover, reconstruction errors of the spine regions separated by the inflection vertebrae had minimal impacts on each other. Overall, multibody scoliosis models with fixed spinal lengths can benefit from the extra rotational joints that contribute toward the accuracy of spine curve approximation. The outcome of this study leads to concurrent accuracy improvement and simplification of multibody models; joint-link configurations can be independently defined for the regions separated by the inflection vertebrae, enabling local optimization of the models for higher accuracy without unnecessary added complexity to the whole model.

Calibration and validation of a novel hybrid model of the lumbosacral spine in ArtiSynth-The passive structures

In computational biomechanics, two separate types of models have been used predominantly to enhance the understanding of the mechanisms of action of the lumbosacral spine (LSS): Finite element (FE) and musculoskeletal multibody (MB) models. To combine advantages of both models, hybrid FE-MB models are an increasingly used alternative. The aim of this paper is to develop, calibrate, and validate a novel passive hybrid FE-MB open-access simulation model of a ligamentous LSS using ArtiSynth. Based on anatomical data from the Male Visible Human Project, the LSS model is constructed from the L1-S1 rigid vertebrae interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, and facet joints. A mesh convergence study, sensitivity analyses, and systematic calibration were conducted with the hybrid functional spinal unit (FSU) L4/5. The predicted mechanical responses of the FSU L4/5, the lumbar spine (L1-L5), and the LSS were validated against literature data from in vivo and in vitro measurements and in silico models. Spinal mechanical responses considered when loaded with pure moments and combined loading modes were total and intervertebral range of motions, instantaneous axes and centers of rotation, facet joint contact forces, intradiscal pressures, disc bulges, and stiffnesses. Undesirable correlations with the FE mesh were minimized, the number of crisscrossed collagen fiber rings was reduced to five, and the individual influences of specific anatomical structures were adjusted to in vitro range of motions. Including intervertebral motion couplings for axial rotation and nonlinear stiffening under increasing axial compression, the predicted kinematic and structural mechanics responses were consistent with the comparative data. The results demonstrate that the hybrid simulation model is robust and efficient in reproducing valid mechanical responses to provide a starting point for upcoming optimizations and extensions, such as with active skeletal muscles.

A Comparison of Clinical Spinal Mobility Measures to Experimentally Derived Lumbar Spine Passive Stiffness

Spinal stiffness and mobility assessments vary between clinical and research settings, potentially hindering the understanding and treatment of low back pain. A total of 71 healthy participants were evaluated using 2 clinical assessments (posteroanterior spring and passive intervertebral motion) and 2 quantitative measures: lumped mechanical stiffness of the lumbar spine and local tissue stiffness (lumbar erector spinae and supraspinous ligament) measured via myotonometry. The authors hypothesized that clinical, mechanical, and local tissue measures would be correlated, that clinical tests would not alter mechanical stiffness, and that males would demonstrate greater lumbar stiffness than females. Clinical, lumped mechanical, and tissue stiffness were not correlated; however, gradings from the posteroanterior spring and passive intervertebral motion tests were positively correlated with each other. Clinical assessments had no effect on lumped mechanical stiffness. The males had greater lumped mechanical and lumbar erector spinae stiffness compared with the females. The lack of correlation between clinical, tissue, and lumped mechanical measures of spinal stiffness indicates that the use of the term "stiffness" by clinicians may require reevaluation; clinicians should be confident that they are not altering mechanical stiffness of the spine through segmental mobility assessments; and greater resting lumbar erector stiffness in males suggests that sex should be considered in the assessment and treatment of the low back.

An EMG-assisted modeling approach to assess passive lumbar tissue loading in vivo during trunk bending

Lower back pain (LBP) is a condition with high prevalence and high cost both in the United States and around the world. The magnitude of mechanical loading on spine is strongly associated with the occurrence of LBP. Previously, to assess spinal loading, biologically assisted biomechanical models were developed to estimate trunk muscle contraction forces. Loadings on lumbar passive tissues are estimated using anatomical models. However, despite the substantial individual variability in lumbar ligament geometry and viscoelastic properties, the existing anatomical models do not account for these differences. As such, the accuracy of model prediction is compromised especially when mid to full range of trunk motions are involved. This paper describes a new modeling approach to assess lumbar passive tissue loading with the consideration of individual differences in lumbar passive tissue properties. A data set that has trunk bending data from 13 human participants was analyzed; on average, lumbar passive tissue contributes to ∼89% of the total spinal compression force at fully flexed trunk postures; the estimated spinal tissue loadings were in feasible ranges as reported from previous cadaver studies; the estimated spinal loadings were also mostly in agreement with results from previous in vivo studies.

Semiautomated 3D Spine Reconstruction from Biplanar Radiographic Images: Prediction of Intervertebral Loading in Scoliotic Subjects

The present study proposes a semiautomatic software approach to reconstruct 3D subject-specific musculoskeletal model of thoracolumbar spine from radiographic digitized images acquired with EOS system. The approach is applied to evaluate the intervertebral loads in 38 standing adolescents with mild idiopathic scoliosis. For each vertebra, a set of landmarks was manually identified on radiographic images. The landmark coordinates were processed to calculate the following vertebral geometrical properties in the 3D space (i) location (ii) dimensions; and (iii) rotations. Spherical joints simulated disks, ligaments, and facet joints. Body weight distribution, muscles forces, and insertion points were placed according to physiological-anatomical values. Inverse static analysis, calculating joints' reactions in maintaining assigned spine configuration, was performed with AnyBody software. Reaction forces were computed to quantify intervertebral loads, and correlation with the patient anatomical parameters was then checked. Preliminary validation was performed comparing the model outcomes with that obtained from other authors in previous modeling works and from in vivo measurements. The comparison with previous modeling works and in vivo studies partially fulfilled the preliminary validation purpose. However, minor incongruities were pointed out that need further investigations. The subjects' intervertebral loads were found significantly correlated with the anatomical parameters in the sagittal and axial planes. Despite preliminary encouraging results that support model suitability, future investigations to consolidate the proposed approach are necessary. Nonetheless, the present method appears to be a promising tool that once fully validated could allow the subject-specific non-invasive evaluation of a deformed spine, providing supplementary information to the routine clinical examination and surgical intervention planning.

Timing and magnitude of lumbar spine contribution to trunk forward bending and backward return in patients with acute low back pain

Alterations in the lumbo-pelvic coordination denote changes in neuromuscular control of trunk motion as well as load sharing between passive and active tissues in the lower back. Differences in timing and magnitude aspects of lumbo-pelvic coordination between patients with chronic low back pain (LBP) and asymptomatic individuals have been reported; yet, the literature on lumbo-pelvic coordination in patients with acute LBP is scant. A case-control study was conducted to explore the differences in timing and magnitude aspects of lumbo-pelvic coordination between females with (n=19) and without (n=19) acute LBP. Participants in each group completed one experimental session wherein they performed trunk forward bending and backward return at preferred and fast paces. The amount of lumbar contribution to trunk motion (as the magnitude aspect) as well as the mean absolute relative phase (MARP) and deviation phase (DP) between thoracic and pelvic rotations (as the timing aspect) of lumbo-pelvic coordination were calculated. The lumbar contribution to trunk motion in the 2nd and the 3rd quarters of both forward bending and backward return phases was significantly smaller in the patient than the control group. The MARP and the DP were smaller in the patient vs. the control group during entire motion. The reduced lumbar contribution to trunk motion as well as the more in-phase and less variable lumbo-pelvic coordination in patients with acute LBP compared to the asymptomatic controls is likely the result of a neuromuscular adaptation to reduce painful deformation and to protect injured lower back tissues.

Fusion angle affects intervertebral adjacent spinal segment joint forces-Model-based analysis of patient specific alignment

This study addresses the hypothesis that adjacent segment intervertebral joint loads are sensitive to the degree of lordosis that is surgically imposed during vertebral fusion. Adjacent segment degeneration is often observed after lumbar fusion, but a causative mechanism is not yet clearly evident. Altered kinematics of the adjacent segments and potentially nonphysiological mechanical joint loads have been implicated in this process. However, little is known of how altered alignment and kinematics influence loading of the adjacent intervertebral joints under consideration of active muscle forces. This study investigated these effects by simulating L4/5 fusions using kinematics-driven musculoskeletal models of one generic and eight sagittal alignment-specific models. Models featured different spinopelvic configurations but were normalized by body height, masses, and muscle properties. Fusion of the L4/5 segment was implemented in an in situ (22°), hyperlordotic (32°), and hypolordotic (8°) fashion and kinematic input parameters were changed accordingly based on findings of an in vitro investigation. Bending motion from upright standing to 45° forward flexion and back was simulated for all models in intact and fused conditions. Joint loads at adjacent levels and moment arms of spinal muscles experienced changes after all types of fusion. Hypolordotic configuration led to an increase of adjacent segment (L3/4) shear forces of 29% on average, whereas hyperlordotic fusion reduced shear by 39%. Overall, L4/5 in situ fusion resulted in intervertebral joint forces closest to intact loading conditions. An artificial decrease in lumbar lordosis (minus 14° on average) caused by an L4/5 fusion lead to adverse loading conditions, particularly at the cranial adjacent levels, and altered muscle moment arms, in particular for muscles in the vicinity of the fusion. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:131-139, 2017.

Preserving Posterior Complex Can Prevent Adjacent Segment Disease following Posterior Lumbar Interbody Fusion Surgeries: A Finite Element Analysis

Objective

To investigate the biomechanical effects of the lumbar posterior complex on the adjacent segments after posterior lumbar interbody fusion (PLIF) surgeries.

Methods

A finite element model of the L1–S1 segment was modified to simulate PLIF with total laminectomy (PLIF-LAM) and PLIF with hemilaminectomy (PLIF-HEMI) procedures. The models were subjected to a 400N follower load with a 7.5-N.m moment of flexion, extension, torsion, and lateral bending. The range of motion (ROM), intradiscal pressure (IDP), and ligament force were compared.

Results

In Flexion, the ROM, IDP and ligament force of posterior longitudinal ligament, intertransverse ligament, and capsular ligament remarkably increased at the proximal adjacent segment in the PLIF-LAM model, and slightly increased in the PLIF-HEMI model. There was almost no difference for the ROM, IDP and ligament force at L5-S1 level between the two PLIF models although the ligament forces of ligamenta flava remarkably increased compared with the intact lumbar spine (INT) model. For the other loading conditions, these two models almost showed no difference in ROM, IDP and ligament force on the adjacent discs.

Conclusions

Preserved posterior complex acts as the posterior tension band during PLIF surgery and results in less ROM, IDP and ligament forces on the proximal adjacent segment in flexion. Preserving the posterior complex during decompression can be effective on preventing adjacent segment degeneration (ASD) following PLIF surgeries.

On the load-sharing along the ligamentous lumbosacral spine in flexed and extended postures: Finite element study

A harmonic synergy between the load-bearing and stabilizing components of the spine is necessary to maintain its normal function. This study aimed to investigate the load-sharing along the ligamentous lumbosacral spine under sagittal loading. A 3D nonlinear detailed Finite Element (FE) model of lumbosacral spine with realistic geometry was developed and validated using wide range of numerical and experimental (in-vivo and in-vitro) data. The model was subjected to 500 N compressive Follower Load (FL) combined with 7.5 Nm flexion (FLX) or extension (EXT) moments. Load-sharing was expressed as percentage of total internal force/moment developed along the spine that each spinal component carried. These internal forces and moments were determined at the discs centres and included the applied load and the resisting forces in the ligaments and facet joints. The contribution of the facet joints and ligaments in supporting bending moments produced additional forces and moments in the discs. The intervertebral discs carried up to 81% and 68% of the total internal force in case of FL combined with FLX and EXT, respectively. The ligaments withstood up to 67% and 81% of the total internal moment in cases of FL combined with EXT and FLX, respectively. Contribution of the facet joints in resisting internal force and moment was noticeable at levels L4-S1 only particularly in case of FL combined with EXT and reached up 29% and 52% of the internal moment and force, respectively. This study demonstrated that spinal load-sharing depended on applied load and varied along the spine.

Passive lumbar tissue loading during trunk bending at three speeds: An in vivo study

BACKGROUND

Low back disorders are closely related with the magnitude of mechanical loading on human spine. However, spinal loading contributed by the lumbar passive tissues is still not well understood. In this study, the effect of motion speed on lumbar passive moment output was investigated. In addition, the increase of lumbar passive moment during trunk bending was modeled.

METHODS

Twelve volunteers performed trunk-bending motions at three different speeds. Trunk kinematics and muscle activities were collected and used to estimate instantaneous spinal loading and the corresponding lumbar passive moment. The lumbar passive moments at different ranges of trunk motion were compared at different speed levels and the relationship between lumbar passive moment lumbar flexion was modeled.

FINDINGS

A non-linear, two-stage pattern of increase in lumbar passive moment was evident during trunk flexion. However, the effect of motion speed was not significant on lumbar passive moments or any of the model parameters.

INTERPRETATION

As reported previously, distinct lumbar ligaments may begin to generate tension at differing extents of trunk flexion, and this could be the cause of the observed two-stage increasing pattern of lumbar passive moment. The current results also suggest that changes in tissue strain rate may not have a significant impact on the total passive moment output at the relatively slow trunk motions examined here.

Capturing Three-Dimensional In Vivo Lumbar Intervertebral Joint Kinematics Using Dynamic Stereo-X-Ray Imaging

Availability of accurate three-dimensional (3D) kinematics of lumbar vertebrae is necessary to understand normal and pathological biomechanics of the lumbar spine. Due to the technical challenges of imaging the lumbar spine motion in vivo, it has been difficult to obtain comprehensive, 3D lumbar kinematics during dynamic functional tasks. The present study demonstrates a recently developed technique to acquire true 3D lumbar vertebral kinematics, in vivo, during a functional load-lifting task. The technique uses a high-speed dynamic stereo-radiography (DSX) system coupled with a volumetric model-based bone tracking procedure. Eight asymptomatic male participants performed weight-lifting tasks, while dynamic X-ray images of their lumbar spines were acquired at 30 fps. A custom-designed radiation attenuator reduced the radiation white-out effect and enhanced the image quality. High resolution CT scans of participants' lumbar spines were obtained to create 3D bone models, which were used to track the X-ray images via a volumetric bone tracking procedure. Continuous 3D intervertebral kinematics from the second lumbar vertebra (L2) to the sacrum (S1) were derived. Results revealed motions occurring simultaneously in all the segments. Differences in contributions to overall lumbar motion from individual segments, particularly L2–L3, L3–L4, and L4–L5, were not statistically significant. However, a reduced contribution from the L5–S1 segment was observed. Segmental extension was nominally linear in the middle range (20%–80%) of motion during the lifting task, but exhibited nonlinear behavior at the beginning and end of the motion. L5–S1 extension exhibited the greatest nonlinearity and variability across participants. Substantial AP translations occurred in all segments (5.0 ± 0.3 mm) and exhibited more scatter and deviation from a nominally linear path compared to segmental extension. Maximum out-of-plane rotations (<1.91 deg) and translations (<0.94 mm) were small compared to the dominant motion in the sagittal plane. The demonstrated success in capturing continuous 3D in vivo lumbar intervertebral kinematics during functional tasks affords the possibility to create a baseline data set for evaluating the lumbar spinal function. The technique can be used to address the gaps in knowledge of lumbar kinematics, to improve the accuracy of the kinematic input into biomechanical models, and to support development of new disk replacement designs more closely replicating the natural lumbar biomechanics.

Computational Biomechanical Modeling of Scoliotic Spine: Challenges and Opportunities

BACKGROUND

Biomechanical computer models of the spine have important roles in the treatment and correction of scoliosis by providing predictive information for surgeons and other clinicians.

OBJECTIVES

This article reviews computational models of intact and scoliotic spine and its components; vertebra, intervertebral disc, ligament, facet joints, and muscle. Several spine models, developed using multi-body modelling and finite element modelling schemes, and their pros and cons are discussed.

CONCLUSIONS

The review reveals that scoliosis modelling is performed for 3 main applications: 1) brace simulation; 2) analysis of surgical correction technique; and 3) patient positioning before surgical instrumentation. The models provide predictive information for a priori choice of brace configurations and mechanically effective surgical correction techniques and the expected degree of correction. However, they have many shortcomings: for instance, they do not fully reproduce the active behaviour of the spine and the models' properties are not personalized.