Moment arms of the human neck muscles in flexion, bending and rotation

There is a paucity of data available for the moment arms of the muscles of the human neck. The objective of the present study was to measure the moment arms of the major cervical spine muscles in vitro. Experiments were performed on five fresh-frozen human head-neck specimens using a custom-designed robotic spine testing apparatus. The testing apparatus replicated flexion-extension, lateral bending and axial rotation of each individual intervertebral joint in the cervical spine while all other joints were kept immobile. The tendon excursion method was used to measure the moment arms of 30 muscle sub-regions involving 13 major muscles of the neck about all three axes of rotation of each joint for the neutral position of the cervical spine. Significant differences in the moment arm were observed across sub-regions of individual muscles and across the intervertebral joints spanned by each muscle (p<0.05). Overall, muscle moment arms were larger in flexion-extension and lateral bending than in axial rotation, and most muscles had prominent moment arms in at least 2 out of the 3 joint motions investigated. This study emphasizes the importance of detailed representation of a muscle's architecture in prediction of its torque capacity about the individual joints of the cervical spine. The dataset produced may be useful in developing and validating computational models of the human neck.     

Sensitivity of muscle and intervertebral disc force computations to variations in muscle attachment sites

Abstract

The current paper aims at assessing the sensitivity of muscle and intervertebral disc force computations against potential errors in modeling muscle attachment sites. We perturbed each attachment location in a complete and coherent musculoskeletal model of the human spine and quantified the changes in muscle and disc forces during standing upright, flexion, lateral bending, and axial rotation of the trunk. Although the majority of the muscles caused minor changes (less than 5%) in the disc forces, certain muscle groups, for example, quadratus lumborum, altered the shear and compressive forces as high as 353% and 17%, respectively. Furthermore, percent changes were higher in the shear forces than in the compressive forces. Our analyses identified certain muscles in the rib cage (intercostales interni and intercostales externi) and lumbar spine (quadratus lumborum and longissimus thoracis) as being more influential for computing muscle and disc forces. Furthermore, the disc forces at the L4/L5 joint were the most sensitive against muscle attachment sites, followed by T6/T7 and T12/L1 joints. Presented findings suggest that modeling muscle attachment sites based on solely anatomical illustrations might lead to erroneous evaluation of internal forces and promote using anatomical datasets where these locations were accurately measured. When developing a personalized model of the spine, certain care should also be paid especially for the muscles indicated in this work.     

Muscle and external load contribution to knee joint contact loads during normal gait

Large knee adduction moments during gait have been implicated as a mechanical factor related to the progression and severity of tibiofemoral osteoarthritis and it has been proposed that these moments increase the load on the medial compartment of the knee joint. However, this mechanism cannot be validated without taking into account the internal forces and moments generated by the muscles and ligaments, which cannot be easily measured. Previous musculoskeletal models suggest that the medial compartment of the tibiofemoral joint bears the majority of the tibiofemoral load, with the lateral compartment unloaded at times during stance. Yet these models did not utilise explicitly measured muscle activation patterns and measurements from an instrumented prosthesis which do not portray lateral compartment unloading. This paper utilised an EMG-driven model to estimate muscle forces and knee joint contact forces during healthy gait. Results indicate that while the medial compartment does bear the majority of the load during stance, muscles provide sufficient stability to counter the tendency of the external adduction moment to unload the lateral compartment. This stability was predominantly provided by the quadriceps, hamstrings, and gastrocnemii muscles, although the contribution from the tensor fascia latae was also significant. Lateral compartment unloading was not predicted by the EMG-driven model, suggesting that muscle activity patterns provide useful input to estimate muscle and joint contact forces.     

Role of trunk muscles in generating follower load in the lumbar spine of neutral standing posture

Recently, experimental results have demonstrated that the load carrying capacity of the human spine substantially increases under the follower load condition. Thus, it is essential to prove that a follower load can be generated in vivo by activating the appropriate muscles in order to demonstrate the possibility that the stability of the spinal column could be maintained through a follower load mechanism. The aim of this study was to analyze the coordination of the trunk muscles in order to understand the role of the muscles in generating the follower load. A three-dimensional finite element model of the lumbar spine was developed from T12 to S1 and 117 pairs of trunk muscles (58 pairs of superficial muscles and 59 pairs of deep muscles) were considered. The follower load concept was mathematically represented as an optimization problem. The muscle forces required to generate the follower load were predicted by solving the optimization problem. The corresponding displacements and rotations at all nodes were estimated along with the follower forces, shear forces, and joint moments acting on those nodes. In addition, the muscle forces and the corresponding responses were investigated when the activations of the deep muscles or the superficial muscles were restricted to 75% of the maximum activation, respectively. Significantly larger numbers of deep muscles were involved in the generation of the follower load than the number of superficial muscles, regardless of the restriction on muscle activation. The shear force and the resultant joint moment are more influenced by the change in muscle activation in the superficial muscles. A larger number of deep trunk muscles were activated in order to maintain the spinal posture in the lumbar spine. In addition, the deep muscles have a larger capability to reduce the shear force and the resultant joint moment with respect to the perturbation of the external load or muscle fatigue compared to the superficial muscles.     

Validation of lumbar spine loading from a musculoskeletal model including the lower limbs and lumbar spine

Low back mechanics are important to quantify to study injury, pain and disability. As in vivo forces are difficult to measure directly, modeling approaches are commonly used to estimate these forces. Validation of model estimates is critical to gain confidence in modeling results across populations of interest, such as people with lower-limb amputation. Motion capture, ground reaction force and electromyographic data were collected from ten participants without an amputation (five male/five female) and five participants with a unilateral transtibial amputation (four male/one female) during trunk-pelvis range of motion trials in flexion/extension, lateral bending and axial rotation. A musculoskeletal model with a detailed lumbar spine and the legs including 294 muscles was used to predict L4-L5 loading and muscle activations using static optimization. Model estimates of L4-L5 intervertebral joint loading were compared to measured intradiscal pressures from the literature and muscle activations were compared to electromyographic signals. Model loading estimates were only significantly different from experimental measurements during trunk extension for males without an amputation and for people with an amputation, which may suggest a greater portion of L4-L5 axial load transfer through the facet joints, as facet loads are not captured by intradiscal pressure transducers. Pressure estimates between the model and previous work were not significantly different for flexion, lateral bending or axial rotation. Timing of model-estimated muscle activations compared well with electromyographic activity of the lumbar paraspinals and upper erector spinae. Validated estimates of low back loading can increase the applicability of musculoskeletal models to clinical diagnosis and treatment.     

Evidence for a role of antagonistic cocontraction in controlling trunk stiffness during lifting

Activity of the abdominal muscles during symmetric lifting has been a consistent finding in many studies. It has been hypothesized that this antagonistic coactivation increases trunk stiffness to provide stability to the spine. To test this, we investigated whether abdominal activity in lifting is increased in response to destabilizing conditions.

Ten healthy male subjects lifted 35 l containers containing 15 l of water (unstable condition), or ice (stable condition). 3D-kinematics, ground reaction forces, and EMG of selected trunk muscles were recorded. Euler angles of the thorax relative to the pelvis were determined. Inverse dynamics was used to calculate moments about L5S1. Averaged normalized abdominal EMG activity was calculated to express coactivation and an EMG-driven trunk muscle model was used to estimate the flexor moment produced by these muscles and to estimate the L5S1 compression force.

Abdominal coactivation was significantly higher when lifting the unstable load. This coincided with significant increases in estimated moments produced by the antagonist muscles and in estimated compression forces on the L5S1 disc, except at the instant of the peak moment about L5S1. The lifting style was not affected by load instability as evidenced by the absence of effects on moments about L5S1 and angles of the thorax relative to the pelvis. The data support the interpretation of abdominal cocontraction during lifting as subserving spinal stability. An alternative function of the increased trunk stiffness due to cocontraction might be to achieve more precise control over the trajectory of lifted weight in order to avoid sloshing of the water mass in the box and the consequent perturbations.     

A frontal plane model of the lumbar spine subjected to a follower load: implications for the role of muscles

Compression on the lumbar spine is 1000 N for standing and walking and is higher during lifting. Ex vivo experiments show it buckles under a vertical load of 80-100 N. Conversely, the whole lumbar spine can support physiologic compressive loads without large displacements when the load is applied along a follower path that approximates the tangent to the curve of the lumbar spine. This study utilized a two-dimensional beam-column model of the lumbar spine in the frontal plane under gravitational and active muscle loads to address the following question: Can trunk muscle activation cause the path of the internal force resultant to approximate the tangent to the spinal curve and allow the lumbar spine to support compressive loads of physiologic magnitudes? The study identified muscle activation patterns that maintained the lumbar spine model under compressive follower load, resulting in the minimization of internal shear forces and bending moments simultaneously at all lumbar levels. The internal force resultant was compressive, and the lumbar spine model, loaded in compression along the follower load path, supported compressive loads of physiologic magnitudes with minimal change in curvature in the frontal plane. Trunk muscles may coactivate to generate a follower load path and allow the ligamentous lumbar spine to support physiologic compressive loads.     

Effects of lumbo-pelvic rhythm on trunk muscle forces and disc loads during forward flexion: A combined musculoskeletal and finite element simulation study

Previous in-vivo studies suggest that the ratio of total lumbar rotation over pelvic rotation (lumbo-pelvic rhythm) during trunk sagittal movement is essential to evaluate spinal loads and discriminate between low back pain and asymptomatic population. Similarly, there is also evidence that the lumbo-pelvic rhythm is key for evaluation of realistic muscle and joint reaction forces and moments predicted by various computational musculoskeletal models. This study investigated the effects of three lumbo-pelvic rhythms defined based on in-vivo measurements on the spinal response during moderate forward flexion (60°) using a combined approach of musculoskeletal modeling of the upper body and finite element model of the lumbosacral spine. The muscle forces and joint loads predicted by the musculoskeletal model, together with the gravitational forces, were applied to the finite element model to compute the disc force and moment, intradiscal pressure, annular fibers strain, and load-sharing. The results revealed that a rhythm with high pelvic rotation and low lumbar flexion involves more global muscles and increases the role of the disc in resisting spinal loads, while its counterpart, with low pelvic rotation, recruits more local muscles and engages the ligaments to lower the disc loads. On the other hand, a normal rhythm that has balanced pelvic and lumbar rotations yields almost equal disc and ligament load-sharing and results in more balanced synergy between global and local muscles. The lumbo-pelvic rhythm has less effect on the intradiscal pressure and annular fibers strain. This work demonstrated that the spinal response during forward flexion is highly dependent on the lumbo-pelvic rhythm. It is therefore, essential to adapt this parameter instead of using the default values in musculoskeletal models for accurate prediction of muscle forces and joint reaction forces and moments. The findings provided by this work are expected to improve knowledge of spinal response during forward flexion, and are clinically relevant towards low back pain treatment and disc injury prevention.     

Development of a novel MATLAB-based framework for implementing mechanical joint stability constraints within OpenSim musculoskeletal models

The Static Optimization (SO) solver in OpenSim estimates muscle activations and forces that only equilibrate applied moments. In this study, SO was enhanced through an open-access MATLAB interface, where calculated muscle activations can additionally satisfy crucial mechanical stability requirements. This Stability-Constrained SO (SCSO) is applicable to many OpenSim models and can potentially produce more biofidelic results than SO alone, especially when antagonistic muscle co-contraction is required to stabilize body joints. This hypothesis was tested using existing models and experimental data in the literature. Muscle activations were calculated by SO and SCSO for a spine model during two series of static trials (i.e. simulation 1 and 2), and also for a lower limb model (supplementary material 2). In simulation 1, symmetric and asymmetric flexion postures were compared, while in simulation 2, various external load heights were compared, where increases in load height did not change the external lumbar flexion moment, but necessitated higher EMG activations. During the tasks in simulation 1, the predicted muscle activations by SCSO demonstrated less average deviation from the EMG data (6.8% −7.5%) compared to those from SO (10.2%). In simulation 2, SO predicts constant muscle activations and forces, while SCSO predicts increases in the average activations of back and abdominal muscles that better match experimental data. Although the SCSO results are sensitive to some parameters (e.g. musculotendon stiffness), when considering the strategy of the central nervous system in distributing muscle forces and in activating antagonistic muscles, the assigned activations by SCSO are more biofidelic than SO.     

Mechanically simulated muscle forces strongly stabilize intact and injured upper cervical spine specimens

Although muscles are assumed to be capable of stabilizing the spinal column in vivo, they have only rarely been simulated in vitro. Their effect might be of particular importance in unstable segments. The present study therefore tests the hypothesis that mechanically simulated muscle forces stabilize intact and injured cervical spine specimens. In the first step, six human occipito-cervical spine specimens were loaded intact in a spine tester with pure moments in lateral bending (+/- 1.5 N m), flexion-extension (+/- 1.5 N m) and axial rotation (+/- 0.5 N m). In the second step, identical flexibility tests were carried out during constant traction of three mechanically simulated muscle pairs: splenius capitits (5 N), semispinalis capitis (5 N) and longus colli (15 N). Both steps were repeated after unilateral and bilateral transection of the alar ligaments. The muscle forces strongly stabilized C0-C2 in all loading and injury states. This was most obvious in axial rotation, where a reduction of range of motion (ROM) and neutral zone to <50% (without muscles=100%) was observed. With increasing injury the normalized ROM (intact condition=100%) increased with and without muscles approximately to the same extend. With bilateral injury this increase was 125-132% in lateral bending, 112%-119% in flexion-extension and 103-116% in axial rotation. Mechanically simulated cervical spine muscles strongly stabilized intact and injured cervical spine specimens. Nevertheless, it could be shown that in vitro flexibility tests without muscle force simulation do not necessarily lead to an overestimation of spinal instability if the results are normalized to the intact state.