Moment-rotation behavior of intervertebral joints in flexion-extension, lateral bending, and axial rotation at all levels of the human spine: A structured review and meta-regression analysis

A new method to design energy-conserving surrogate models for the coupled, nonlinear responses of intervertebral discs

The aim of this study was to design physics-preserving and precise surrogate models of the nonlinear elastic behaviour of an intervertebral disc (IVD). Based on artificial force-displacement data sets from detailed finite element (FE) disc models, we used greedy kernel and polynomial approximations of second, third and fourth order to train surrogate models for the scalar force-torque -potential. Doing so, the resulting models of the elastic IVD responses ensured the conservation of mechanical energy through their structure. At the same time, they were capable of predicting disc forces in a physiological range of motion and for the coupling of all six degrees of freedom of an intervertebral joint. The performance of all surrogate models for a subject-specific L4 | 5 disc geometry was evaluated both on training and test data obtained from uncoupled (one-dimensional), weakly coupled (two-dimensional), and random movement trajectories in the entire six-dimensional (6d) physiological displacement range, as well as on synthetic kinematic data. We observed highest precisions for the kernel surrogate followed by the fourth-order polynomial model. Both clearly outperformed the second-order polynomial model which is equivalent to the commonly used stiffness matrix in neuro-musculoskeletal simulations. Hence, the proposed model architectures have the potential to improve the accuracy and, therewith, validity of load predictions in neuro-musculoskeletal spine models.

Variability of intervertebral joint stiffness between specimens and spine levels

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

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.

Diabetes-related Screw Loosening: The Distinction of Surgical Sites and the Relationship among Diabetes, Implant Stabilization and Clinical Outcomes

OBJECTIVES

Diabetes mellitus (DM) is correlated with poor clinical outcomes in spinal surgery. However, the effect of it on screw stabilization has not been investigated. The aim of this study was to evaluate the screw loosening rate and postoperative outcomes in diabetic patients and to identify potential risk factors associated with loosening.

METHODS

This was a retrospective study. Two hundred and forty-three patients who received cervical or lumbar internal fixation between 2015 and 2019 were enrolled. Screw loosening was assessed on radiography, and clinical outcomes were evaluated by the improvement of visual analogue scale (VAS), Oswestry disability index (ODI) or Japanese Orthopaedic Association (JOA) scores. The relationship of DM, screw loosening and clinical outcomes were analyzed with chi-square tests and regression analyses.

RESULTS

One hundred and twenty-two patients (50.2%) with diabetes were included in this study. Diabetes led to the increase of the rate of screw loosening in the lumbar spine, while the loosening rate did not vary significantly in the cervical spine. The occurrence of screw loosening in the lumbar spine was more likely to be associated with clinical outcomes for motor performance including walking and sitting. However, no significant effect on JOA and VAS scores in the cervical spine of screw loosening was found. Moreover, the history of DM affected the outcomes of the patients who underwent spinal surgery.

CONCLUSION

DM had an adverse effect on screw stabilization. The impaired improvement of clinical outcomes in diabetics after spinal surgery was related to screw loosening. In addition to the direct effects on operative wounds and neural function, the impact on the screws due to DM was also worth noting.

A novel artificial vertebral implant with Gyroid porous structures for reducing the subsidence and mechanical failure rate after vertebral body replacement

BACKGROUND

Prosthesis subsidence and mechanical failure were considered significant threats after vertebral body replacement during the long-term follow-up. Therefore, improving and optimizing the structure of vertebral substitutes for exceptional performance has become a pivotal challenge in spinal reconstruction.

METHODS

The study aimed to develop a novel artificial vertebral implant (AVI) with triply periodic minimal surface Gyroid porous structures to enhance the safety and stability of prostheses. The biomechanical performance of AVIs under different loading conditions was analyzed using the finite element method. These implants were fabricated using selective laser melting technology and evaluated through static compression and subsidence experiments.

RESULTS

The results demonstrated that the peak stress in the Gyroid porous AVI was consistently lower than that in the traditional porous AVI under all loading conditions, with a maximum reduction of 73.4%. Additionally, it effectively reduced peak stress at the bone-implant interface of the vertebrae. Static compression experiments demonstrated that the Gyroid porous AVI was about 1.63 times to traditional porous AVI in terms of the maximum compression load, indicating that Gyroid porous AVI could meet the safety requirement. Furthermore, static subsidence experiments revealed that the subsidence tendency of Gyroid porous AVI in polyurethane foam (simulated cancellous bone) was approximately 15.7% lower than that of traditional porous AVI.

CONCLUSIONS

The Gyroid porous AVI exhibited higher compressive strength and lower subsidence tendency than the strut-based traditional porous AVI, indicating it may be a promising substitute for spinal reconstruction.

Off-Axis Loading Fixture for Spine Biomechanics: Combined Compression and Bending

The spine is a multi-tissue musculoskeletal system that supports large multi-axial loads and motions during physiological activities. The healthy and pathological biomechanical function of the spine and its subtissues are generally studied using cadaveric specimens that often require multi-axis biomechanical test systems to mimic the complex loading environment of the spine. Unfortunately, an off-the-shelf device can easily exceed 200,000 USD, while a custom device requires extensive time and experience in mechatronics. Our goal was to develop a cost-appropriate compression and bending (flexion-extension and lateral bending) spine testing system that requires little time and minimal technical knowledge. Our solution was an off-axis loading fixture (OLaF) that mounts to an existing uni-axial test frame and requires no additional actuators. OLaF requires little machining, with most components purchased off-the-shelf, and costs less than 10,000 USD. The only external transducer required is a six-axis load cell. Furthermore, OLaF is controlled using the existing uni-axial test frame's software, while the load data is collected using the software included with the six-axis load cell. Here we provide the design rationale for how OLaF develops primary motions and loads and minimizes off-axis secondary constraints, verify the primary kinematics using motion capture, and demonstrate that the system is capable of applying physiologically relevant, noninjurious, axial compression and bending. While OLaF is limited to compression and bending studies it produces repeatable physiologically relevant biomechanics, with high quality data, and minimal startup costs.

Development and Validation of a Novel In Vitro Joint Testing System for Reproduction of In Vivo Dynamic Muscle Force

In vitro biomechanical experiments utilizing cadaveric specimens are one of the most effective methods for rehearsing surgical procedures, testing implants, and guiding postoperative rehabilitation. Applying dynamic physiological muscle force to the specimens is a challenge to reconstructing the environment of bionic mechanics in vivo, which is often ignored in the in vitro experiment. The current work aims to establish a hardware platform and numerical computation methods to reproduce dynamic muscle forces that can be applied to mechanical testing on in vitro specimens. Dynamic muscle loading is simulated through numerical computation, and the inputs of the platform will be derived. Then, the accuracy and robustness of the platform will be evaluated through actual muscle loading tests in vitro. The tests were run on three muscles (gastrocnemius lateralis, the rectus femoris, and the semitendinosus) around the knee joint and the results showed that the platform can accurately reproduce the magnitude of muscle strength (errors range from -6.2% to 1.81%) and changing pattern (goodness-of-fit range coefficient ranges from 0.00 to 0.06) of target muscle forces. The robustness of the platform is mainly manifested in that the platform can still accurately reproduce muscle force after changing the hardware combination. Additionally, the standard deviation of repeated test results is very small (standard ranges of hardware combination 1: 0.34 N~2.79 N vs. hardware combination 2: 0.68 N~2.93 N). Thus, the platform can stably and accurately reproduce muscle forces in vitro, and it has great potential to be applied in the future musculoskeletal loading system.

Effects of occipital-atlas stabilization on the upper cervical spine rotation combinations: an in vitro study

The purpose of this study is to compare axial rotation range of motion for the upper cervical spine during three movements: axial rotation, rotation + flexion + ipsilateral lateral bending and rotation + extension + contralateral lateral bending before and after occiput-atlas (C0-C1) stabilization. Ten cryopreserved C0-C2 specimens (mean age 74 years, range 63-85 years) were manually mobilized in 1. axial rotation, 2. rotation + flexion + ipsilateral lateral bending and 3. rotation + extension + contralateral lateral bending without and with a screw stabilization of C0-C1. Upper cervical range of motion and the force used to generate the motion were measured using an optical motion system and a load cell respectively. The range of motion (ROM) without C0-C1 stabilization was 9.8° ± 3.9° in right rotation + flexion + ipsilateral lateral bending and 15.5° ± 5.9° in left rotation + flexion + ipsilateral lateral bending. With stabilization, the ROM was 6.7° ± 4.3° and 13.6° ± 5.3°, respectively. The ROM without C0-C1 stabilization was 35.1° ± 6.0° in right rotation + extension + contralateral lateral bending and 29.0° ± 6.5° in left rotation + extension + contralateral lateral bending. With stabilization, the ROM was 25.7° ± 6.4° (p = 0.007) and 25.3° ± 7.1°, respectively. Neither rotation + flexion + ipsilateral lateral bending (left or right) or left rotation + extension + contralateral lateral bending reached statistical significance. ROM without C0-C1 stabilization was 33.9° ± 6.7° in right rotation and 28.0° ± 6.9° in left rotation. With stabilization, the ROM was 28.5° ± 7.0° (p = 0.005) and 23.7° ± 8.5° (p = 0.013) respectively. The stabilization of C0-C1 reduced the upper cervical axial rotation in right rotation + extension + contralateral lateral bending and right and left axial rotations; however, this reduction was not present in left rotation + extension + contralateral lateral bending or both combinations of rotation + flexion + ipsilateral lateral bending.

In Silico Meta-Analysis of Boundary Conditions for Experimental Tests on the Lumbar Spine

The study of the spine range of motion under given external load has been the object of many studies in literature, finalised to a better understanding of the spine biomechanics, its physiology, eventual pathologic conditions and possible rehabilitation strategies. However, the huge amount of experimental work performed so far cannot be straightforwardly analysed due to significant differences among loading set-ups. This work performs a meta-analysis of various boundary conditions in literature, focusing on the flexion/extension behaviour of the lumbar spine. The comparison among range of motions is performed virtually through a validated multibody model. Results clearly illustrated the effect of various boundary conditions which can be met in literature, so justifying differences of biomechanical behaviours reported by authors implementing different set-up: for example, a higher value of the follower load can indeed result in a stiffer behaviour; the application of force producing spurious moments results in an apparently more deformable behaviour, however the respective effects change at various segments along the spine due to its natural curvature. These outcomes are reported not only in qualitative, but also in quantitative terms. The numerical approach here followed to perform the meta-analysis is original and it proved to be effective thanks to the bypass of the natural variability among specimens which might completely or partially hinder the effect of some boundary conditions. In addition, it can provide very complete information since the behaviour of each functional spinal unit can be recorded. On the whole, the work provided an extensive review of lumbar spine loading in flexion/extension.

A finite element study on the effects of follower load on the continuous biomechanical responses of subaxial cervical spine

In spine biomechanics, follower loads are used to mimic the in vivo muscle forces acting on a human spine. However, the effects of the follower load on the continuous biomechanical responses of the subaxial cervical spines (C2-T1) have not been systematically clarified. This study aims at investigating the follower load effects on the continuous biomechanical responses of C2-T1. A nonlinear finite element model is reconstructed and validated for C2-T1. Six levels follower loads are considered along the follower load path that is optimized through a novel range of motion-based method. A moment up to 2 Nm is subsequently superimposed to produce motions in three anatomical planes. The continuous biomechanical responses, including the range of motion, facet joint force, intradiscal pressure and flexibility are evaluated for each motion segment. In the sagittal plane, the change of the overall range of motion arising from the follower loads is less than 6%. In the other two anatomical planes, both the magnitude and shape of the rotation-moment curves change with follower loads. At the neutral position, over 50% decrease in flexibility occurs as the follower load increases from zero to 250 N. In all three anatomical planes, over 50% and 30% decreases in flexibility occur in the first 0.5 Nm for small (≤100 N) and large (≥150 N) follower loads, respectively. Moreover, follower loads tend to increase both the facet joint forces and the intradiscal pressures. The shape of the intradiscal pressure-moment curves changes from nonlinear to roughly linear with increased follower load, especially in the coronal and transverse planes. The results obtained in this work provide a comprehensive understanding on the effects of follower load on the continuous biomechanical responses of the C2-T1.

Influence of Axial Load and a 45-Degree Flexion Head Position on Cervical Spinal Stiffness in Healthy Young Adults

Background: Neck pain is a major cause of disability worldwide. Poor neck posture such as using a smartphone or work-related additional cervical axial load, such headgear of aviators, can cause neck pain. This study aimed at investigating the role of head posture or additional axial load on spinal stiffness, a proxy measure to assess cervical motor control. Methods: The posterior-to-anterior cervical spinal stiffness of 49 young healthy male military employees [mean (SD) age 20 ± 1 years] was measured in two head positions: neutral and 45-degree flexed head position and two loading conditions: with and without additional 3 kg axial load. Each test condition comprised three trials. Measurements were taken at three cervical locations, i.e., spinous processes C2 and C7 and mid-cervical (MC). Results: Cervical spinal stiffness measurements showed good reliability in all test conditions. There was a significant three-way interaction between location × head position × load [F(2, 576) = 9.305, p < 0.001]. Significant two-way interactions were found between measurement locations × loading [F(2, 576) = 15.688, p < 0.001] and measurement locations × head position [F(2, 576) = 9.263, p < 0.001]. There was no significant interaction between loading × head position [F(1, 576) = 0.692, p = 0.406]. Post hoc analysis showed reduction of stiffness in all three measurement locations in flexion position. There was a decrease in stiffness in C2 with loading, increase in stiffness in C7 and no change in MC. Discussion: A flexed head posture leading to decreased stiffness of the cervical spine might contribute to neck pain, especially if the posture is prolonged and static, such as is the case with smartphone users. Regarding the additional load, stiffness decreased high cervical and increased low cervical. There was no change mid cervical. The lower spinal stiffness at the high cervical spine might be caused by capsular ligament laxity due to the buckling effect. At the lower cervical spine, the buckling effect seems to be less dominant, because the proximity to the ribs and sternum provide additional stiffness.

Intersegmental Kinematics of the Upper Cervical Spine: Normal Range of Motion and Its Alteration After Alar Ligament Transection

STUDY DESIGN

Biomechanical study using cadaveric cervical spines.

OBJECTIVE

To evaluate joint mobility and stiffness at the craniovertebral junction.

SUMMARY OF BACKGROUND DATA

Data on the intersegmental kinematics of the craniovertebral joints are available in the literature with a widespread range of values. The effect that alar ligament injuries have on intersegmental kinematics remains unclear and requires further biomechanical investigation.

METHODS

Ten occipito-atlanto-axial (C0-C1-C2) human specimens were articulated to flexion, extension, bilateral lateral bending, and bilateral axial rotation. The moment-rotation response was continuously tracked through the entire range of motion before and after unilateral alar ligament transection of the right side.

RESULTS

The intersegmental (C0-C1/C1-C2) moment-rotation response was continuously quantified in full flexion (7.2 ± 6.6°/12.1 ± 5.8°), extension (11.1 ± 6.4°/3.0 ± 2.8°), lateral bending to the right (3.1 ± 2.2°/1.6 ± 1.2°) and left sides (3.3 ± 1.6°/2.1 ± 1.5°), and axial rotation to the right (1.2 ± 3.5°/32.3 ± 9.3°) and left sides (2.7 ± 2.6°/25.3 ± 8.3°). Unilateral alar ligament transection increased the range of motion of C0-C2 in the three planes of movement; however, intersegmental motion alterations were not always observed.

CONCLUSION

Increases in the range of extension and lateral bending at C0-C1, which had not been reported previously, were observed. Further, the range of rotation on the right and left sides increased, in conjunction with the increased ranges at C0-C1 and C1-C2.Level of Evidence: N/A.

Biomechanical effects of lumbar fusion surgery on adjacent segments using musculoskeletal models of the intact, degenerated and fused spine

Adjacent segment disorders are prevalent in patients following a spinal fusion surgery. Postoperative alterations in the adjacent segment biomechanics play a role in the etiology of these conditions. While experimental approaches fail to directly quantify spinal loads, previous modeling studies have numerous shortcomings when simulating the complex structures of the spine and the pre/postoperative mechanobiology of the patient. The biomechanical effects of the L4–L5 fusion surgery on muscle forces and adjacent segment kinetics (compression, shear, and moment) were investigated using a validated musculoskeletal model. The model was driven by in vivo kinematics for both preoperative (intact or severely degenerated L4–L5) and postoperative conditions while accounting for muscle atrophies. Results indicated marked changes in the kinetics of adjacent L3–L4 and L5–S1 segments (e.g., by up to 115% and 73% in shear loads and passive moments, respectively) that depended on the preoperative L4–L5 disc condition, postoperative lumbopelvic kinematics and, to a lesser extent, postoperative changes in the L4–L5 segmental lordosis and muscle injuries. Upper adjacent segment was more affected post-fusion than the lower one. While these findings identify risk factors for adjacent segment disorders, they indicate that surgical and postoperative rehabilitation interventions should focus on the preservation/restoration of patient’s normal segmental kinematics.

Topological design and biomechanical evaluation for 3D printed multi-segment artificial vertebral implants

Customized spinal implants fabricated by additive manufacturing have been increasingly used clinically to restore the physiological functions. However, the mechanisms and methods about the design for the spinal implants are not clear, especially for the reconstruction of multi-segment vertebral. This study aims to develop a novel multi-objective optimization methodology based on various normal spinal activities, to design the artificial vertebral implant (AVI) with lightweight, high-strength and high-stability. The biomechanical performance for two types of AVI was analyzed and compared under different loading conditions by finite element method. These implants were manufactured via selective laser melting technology and evaluated via compressive testing. Results showed the maximum Mises stress of the optimized implant under various load cases were about 41.5% of that of the trussed implant, and below fatigue strength of 3D printed titanium materials. The optimized implant was about 2 times to trussed implant in term of the maximum compression load and compression stiffness to per unit mass, which indicated the optimized implant can meet the safety requirement. Finally, the optimized implant has been used in clinical practice and good short-term clinical outcomes were achieved. Therefore, the novel developed method provides a favorable guarantee for the design of 3D printed multi-segment artificial vertebral implants.

In vivo measurements of spinal stiffness according to a stepwise increase of axial load

Background

The spine has a complex motor control. Its different stabilization mechanisms through passive, active, and neurological subsystems may result in spinal stiffness. To better understand lumbar spinal motor control, this study aimed to measure the effects of increasing the axial load on spinal stiffness.

Methods

A total of 19 healthy young participants (mean age, 24 ± 2.1 years; 8 males and 11 females) were assessed in an upright standing position. Under different axial loads, the posterior-to-anterior spinal stiffness of the thoracic and lumbar spine was measured. Loads were 0%, 10%, 45%, and 80% of the participant’s body weight.

Results

Data were normally distributed and showed excellent reliability. A repeated-measures analysis of variance with a Greenhouse–Geisser correction showed an effect of the loading condition on the mean spinal stiffness [F (2.6, 744) = 3.456, p < 0.001]. Vertebrae and loading had no interaction [F (2.6, 741) = 0.656, p = 0.559]. Post hoc tests using Bonferroni correction revealed no changes with 10% loading (p = 1.000), and with every additional step of loading, spinal stiffness decreased: 0% or 10–45% loading (p < 0.001), 0% or 10–80% loading (p < 0.001), and 45–80% (p < 0.001).

Conclusion

We conclude that a load of ≥ 45% of the participant’s body weight can lead to changes in the spinal motor control. An axial load of 10% showed no significant changes. Rehabilitation should include high-axial-load exercise if needed in everyday living.

Patterns of Load-to-Strength Ratios Along the Spine in a Population-Based Cohort to Evaluate the Contribution of Spinal Loading to Vertebral Fractures

Vertebral fractures (VFx) are common among older adults. Epidemiological studies report high occurrence of VFx at mid-thoracic and thoracolumbar regions of the spine; however, reasons for this observation remain poorly understood. Prior reports of high ratios of spinal loading to vertebral strength in the thoracolumbar region suggest a possible biomechanical explanation. However, no studies have evaluated load-to-strength ratios (LSRs) throughout the spine for a large number of activities in a sizeable cohort. Thus, we performed a cross-sectional study in a sample of adult men and women from a population-based cohort to: 1) determine which activities cause the largest vertebral LSRs, and 2) examine patterns of LSRs along the spine for these high-load activities. We used subject-specific musculoskeletal models of the trunk to determine vertebral compressive loads for 109 activities in 250 individuals (aged 41 to 90 years, 50% women) from the Framingham Heart Study. Vertebral compressive strengths from T₄ to L₄ were calculated from computed tomography-based vertebral size and bone density measurements. We determined which activities caused maximum LSRs at each of these spinal levels. We identified nine activities that accounted for >95% of the maximum LSRs overall and at least 89.6% at each spinal level. The activity with the highest LSR varied by spinal level, and three distinct spinal regions could be identified by the activity producing maximum LSRs: lateral bending with a weight in one hand (upper thoracic), holding weights with elbows flexed (lower thoracic), and forward flexion with weight (lumbar). This study highlights the need to consider a range of lifting, holding, and non-symmetric activities when evaluating vertebral LSRs. Moreover, we identified key activities that produce higher loading in multiple regions of the spine. These results provide the first guidance on what activities to consider when evaluating vertebral load-to-strength ratios in future studies, including those examining dynamic motions and the biomechanics of VFx. © 2020 American Society for Bone and Mineral Research (ASBMR).

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

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

Selecting the appropriate input variables in a regression approach to estimate actively generated muscle moments around L5/S1 for exoskeleton control

Back support exoskeletons are designed to prevent work-related low-back pain by reducing mechanical loading. For actuated exoskeletons, support based on moments actively produced by the trunk muscles appears a viable approach. The moment can be estimated by a biomechanical model. However, one of the main challenges here is the feasibility of recording the required input variables (kinematics, EMG data, ground reaction forces) to run the model. The aim of this study was to evaluate how accurate different selections of input variables can estimate actively generated moments around L5/S1. Different multivariate regression analyses were performed using a dataset consisting of spinal load, body kinematics and trunk muscle activation levels during different lifting conditions with and without an exoskeleton. The accuracy of the resulting models depended on the number and type of input variables and the regression model order. The current study suggests that third-order polynomial regression of EMG signals of one or two bilateral back muscle pairs together with exoskeleton trunk and hip angle suffices to accurately estimate the actively generated muscle moment around L5/S1, and thereby design a proper control system for back support exoskeletons.

Spinal Compressive Forces in Adolescent Idiopathic Scoliosis With and Without Carrying Loads: A Musculoskeletal Modeling Study

The pathomechanisms of curve progression in adolescent idiopathic scoliosis (AIS) remain poorly understood and biomechanical data are limited. A deeper insight into spinal loading could provide valuable information toward the improvement of current treatment strategies. This work therefore aimed at using subject-specific musculoskeletal full-body models of patients with AIS to predict segmental compressive forces around the curve apex and to investigate how these forces are affected by simulated load carrying. Models were created based on spatially calibrated biplanar radiographic images from 24 patients with mild to moderate AIS and validated by comparing predictions of paravertebral muscle activity with reported values from in vivo studies. Spinal compressive forces were predicted during unloaded upright standing as well as standing with external loads of 10, 15, and 20% of body weight (BW) applied to the scapulae to simulate carrying a backpack in the regular way on the back as well as in front of the body and over the shoulder on the concave and convex sides of the scoliotic curve. The predicted muscle activities around the curve apex were higher on the convex side for the erector spinae (ES) and multifidi (MF) muscles, which was comparable to the EMG-based in vivo measurements from the literature. In terms of spinal loading, the implementation of spinal deformity resulted in a 10% increase of compressive force at the curve apex during unloaded upright standing. Apical compressive forces further increased by 50–62% for a simulated 10% BW load and by 77–94% and 103–128% for 15% and 20% BW loads, respectively. Moreover, load-dependent compressive force increases were the lowest in the regular backpack and the highest in the frontpack and convex conditions, with concave side-carrying forces in between. The predictions indicated increased segmental compressive forces during unloaded upright standing, which could be ascribed to the scoliotic deformation. When carrying loads, compressive forces further increased depending on the carrying mode and the weight of the load. These results can be used as a basis for further studies investigating segmental loading in AIS patients during functional activities. Models can thereby be created using the same approach as proposed in this study.