Relationships between individual isometric muscle forces, EMG activity and joint torque in monkeys

The aim of the present work was to determine the EMG activity and the moment of force developed by the main elbow flexor muscles, and to establish on this basis the degree of their participation in isometric contractions performed at various positions of the elbow. This was achieved by recording the following biomechanical parameters: EMG and tensile stress (or force) from biceps brachii (BB) and brachioradialis (BR); EMG from brachialis; external resultant force (FE). There was: a linear or quadratic relationship between the integrated EMG from each muscle and FE; a linear relationship between the force produced by BB or BR and FE. The slope of these relationships depended on the elbow angle, except for that between BB force and FE. It is proposed that iEMG changes compensate for those of the force lever arm. It has been calculated that the contribution of BR to external torque decreased from the extension to flexion while that of BB increased from 70 degrees to 90 degrees and then decreased. How far these data can be extrapolated to man is a matter of discussion based on iEMG and anthropometrical data.     

Effect of intervertebral translational flexibilities on estimations of trunk muscle forces,kinematics, loads,and stability

Due to the complexity of the human spinal motion segments, the intervertebral joints are often simulated in the musculoskeletal trunk models as pivots thus allowing no translational degrees of freedom (DOFs). This work aims to investigate, for the first time, the effect of such widely used assumption on trunk muscle forces, spinal loads, kinematics, and stability during a number of static activities. To address this, the shear deformable beam elements used in our nonlinear finite element (OFE) musculoskeletal model of the trunk were either substantially stiffened in translational directions (SFE model) or replaced by hinge joints interconnected through rotational springs (HFE model). Results indicated that ignoring intervertebral translational DOFs had in general low to moderate impact on model predictions. Compared with the OFE model, the SFE and HFE models predicted generally larger L4–L5 and L5–S1 compression and shear loads, especially for tasks with greater trunk angles; differences reached ～15% for the L4–L5 compression, ～36% for the L4–L5 shear and ～18% for the L5–S1 shear loads. Such differences increased, as location of the hinge joints in the HFE model moved from the mid-disc height to either the lower or upper endplates. Stability analyses of these models for some select activities revealed small changes in predicted margin of stability. Model studies dealing exclusively with the estimation of spinal loads and/or stability may, hence with small loss of accuracy, neglect intervertebral translational DOFs at smaller trunk flexion angles for the sake of computational simplicity.     

Response of trunk muscle coactivation to changes in spinal stability

The goal of this effort was to assess the neuromuscular response to changes in spinal stability. Biomechanical models suggest that antagonistic co-contraction may be related to stability constraints during lifting exertions. A two-dimensional biomechanical model of spinal equilibrium and stability was developed to predict trunk muscle co-contraction as a function of lifting height and external load. The model predicted antagonistic co-contraction must increase with potential energy of the system even when the external moment was maintained at a constant value. Predicted trends were compared with measured electromyographic (EMG) data recorded during static trunk extension exertions wherein subjects held weighted barbells at specific horizontal and vertical locations relative to the lumbo-sacral spine junction. The task was designed to assure the applied moment was identical during each height condition, thereby changing potential energy without influencing moment. Measured EMG activity in the trunk flexors increased with height of the external load as predicted by the model. Gender difference in spinal stability were also noted. Results empirically demonstrate that the neuromuscular system responds to changes in spinal stability and provide insight into the recruitment of trunk muscle activity.     

Neuromusculoskeletal modeling: estimation of muscle forces and joint moments and movements from measurements of neural command

This paper provides an overview of forward dynamic neuromusculoskeletal modeling. The aim of such models is to estimate or predict muscle forces, joint moments, and/or joint kinematics from neural signals. This is a four-step process. In the first step, muscle activation dynamics govern the transformation from the neural signal to a measure of muscle activation-a time varying parameter between 0 and 1. In the second step, muscle contraction dynamics characterize how muscle activations are transformed into muscle forces. The third step requires a model of the musculoskeletal geometry to transform muscle forces to joint moments. Finally, the equations of motion allow joint moments to be transformed into joint movements. Each step involves complex nonlinear relationships. The focus of this paper is on the details involved in the first two steps, since these are the most challenging to the biomechanician. The global process is then explained through applications to the study of predicting isometric elbow moments and dynamic knee kinetics.     

The effect of perturbing body segment parameters on calculated joint moments and muscle forces during gait

This study examined the effect of body segment parameter (BSP) perturbations on joint moments calculated using an inverse dynamics procedure and muscle forces calculated using computed muscle control (CMC) during gait. BSP (i.e. segment mass, center of mass location (com) and inertia tensor) of the left thigh, shank and foot of a scaled musculoskeletal model were perturbed. These perturbations started from their nominal value and were adjusted to ±40% in steps of 10%, for both individual as well as combined perturbations in BSP. For all perturbations, an inverse dynamics procedure calculated the ankle, knee and hip moments based on an identical inverse kinematics solution. Furthermore, the effect of applying a residual reduction algorithm (RRA) was investigated. Muscle excitations and resulting muscle forces were calculated using CMC. The results show only a limited effect of an individual parameter perturbation on the calculated moments, where the largest effect is found when perturbing the shank com (MS_com,shank, the ratio of absolute difference in torque and relative parameter perturbation, is maximally −7.81 N m for hip flexion moment). The additional influence of perturbing two parameters simultaneously is small (MS_{mass+com,thigh} is maximally 15.2 N m for hip flexion moment). RRA made small changes to the model to increase the dynamic consistency of the simulation (after RRA MS_com,shank is maximally 5.01 N m). CMC results show large differences in muscle forces when BSP are perturbed. These result from the underlying forward integration of the dynamic equations.     

Total trunk muscle force and spinal compression are lower in asymmetric moments as compared to pure extension moments.

The aim of the present study was to test the assumption that asymmetric trunk loading requires a higher total muscle force and consequently entails a higher compression forces on the spine as compared to symmetric loading. When the trunk musculature is modelled in sufficient detail, optimisation shows that there is no mechanical necessity for an increase in total muscle force (or compression force) with task asymmetry. A physiologically based optimisation does also not predict an increase in total muscle force or spinal loading with asymmetry. EMG data on 14 trunk muscles collected in eight subjects showed antagonistic coactivity to be present in both conditions. However, estimates of total muscle force based on the EMG were lower when producing an asymmetric moment. In conclusion, producing an asymmetric moment appears to cause slightly lower forces on the lumbosacral joint as compared to a symmetric moment. Only lateral shear forces increase with asymmetry but these remain well below failure levels.     

An EMG-driven musculoskeletal model to estimate muscle forces and knee joint moments in vivo

This paper examined if an electromyography (EMG) driven musculoskeletal model of the human knee could be used to predict knee moments, calculated using inverse dynamics, across a varied range of dynamic contractile conditions. Muscle-tendon lengths and moment arms of 13 muscles crossing the knee joint were determined from joint kinematics using a three-dimensional anatomical model of the lower limb. Muscle activation was determined using a second-order discrete non-linear model using rectified and low-pass filtered EMG as input. A modified Hill-type muscle model was used to calculate individual muscle forces using activation and muscle tendon lengths as inputs. The model was calibrated to six individuals by altering a set of physiologically based parameters using mathematical optimisation to match the net flexion/extension (FE) muscle moment with those measured by inverse dynamics. The model was calibrated for each subject using 5 different tasks, including passive and active FE in an isokinetic dynamometer, running, and cutting manoeuvres recorded using three-dimensional motion analysis. Once calibrated, the model was used to predict the FE moments, estimated via inverse dynamics, from over 200 isokinetic dynamometer, running and sidestepping tasks. The inverse dynamics joint moments were predicted with an average R(2) of 0.91 and mean residual error of approximately 12 Nm. A re-calibration of only the EMG-to-activation parameters revealed FE moments prediction across weeks of similar accuracy. Changing the muscle model to one that is more physiologically correct produced better predictions. The modelling method presented represents a good way to estimate in vivo muscle forces during movement tasks.     

A biomechanical model on muscle forces in the transfer of spinal load to the pelvis and legs.

Based on musculoskeletal anatomy of the lower back, abdominal wall, pelvis and upper legs, a biomechanical model has been developed on forces in the load transfer through the pelvis. The aim of this model is to obtain a tool for analyzing the relations between forces in muscles, ligaments and joints in the transfer of gravitational and external load from the upper body via the sacroiliac joints to the legs in normal situations and pathology. The study of the relation between muscle coordination patterns and forces in pelvic structures, in particular the sacroiliac joints, is relevant for a better understanding of the aetiology of low back pain and pelvic pain. The model comprises 94 muscle parts, 6 ligaments and 6 joints. It enables the calculation of forces in pelvic structures in various postures. The calculations are based on a linear/non-linear optimization scheme. To gain a better understanding of the function of individual muscles and ligaments, deviant properties of these structures can be preset. The model is validated by comparing calculations with EMG data from the literature. For agonistic muscles, good agreement is found between model calculations and EMG data. Antagonistic muscle activity is underestimated by the model. Imposed activity of modelled antagonistic muscles has a minor effect on the mutual proportions of agonistic muscle activities. Simulation of asymmetric muscle weakness shows higher activity of especially abdominal muscles.     

EMG muscle scanning: stability of hand-held surface electrodes

Electromyographic (EMG) muscle scanning measures 2-second samples of integrated muscle action potentials from individual neck and back muscles using a hand-held scanner with post-style surface electrodes separated by a fixed distance. This "scanning" technique is widely used to expeditiously assess muscle activity in the diagnosis of musculoskeletal disorders. In order to determine if the 2-second sample is sufficiently representative of electrical activity at a specific muscle site, the stability of the signal received by the hand-held scanner was measured bilaterally at six neck and back muscle sites over 40 seconds (20 2-second integration periods) in five seated subjects. Taking the overall average EMG activity as the "true" value, the mean number of 2-second integration periods required to achieve less than 5% standard error was calculated to be 1.47 for the 60 muscles tested. Only three sites required more than five integration periods. The validity of EMG scanning as a diagnostic tool is enhanced by longer integration periods.     

Effects of trunk muscle fatigue and load timing on spinal responses during sudden hand loading

The purpose of this study was to investigate the responses of the spine during sudden loading in the presence of back and abdominal muscle fatigue, with a primary focus on the implications for spinal stability. Fifteen females were studied and each received sudden loads to the hands, at both known and unknown times. Participants received these loading trials (a) while rested, (b) with back muscle fatigue, and (c) with a combination of back and abdominal muscle fatigue. Measures were taken on the EMG activity of two trunk extensor and two abdominal muscles, and on the trunk angle and centre of pressure. A 3 × 2 Repeated Measures ANOVA was also performed. There were no preparations made prior to the perturbation even when it could be anticipated. However, the peak responses that followed were greater in the unexpected versus the expected condition. In addition, trunk muscle fatigue led to an increase in the baseline activity of the trunk muscles but no additional increase in activity just prior to loading. There was increased activation of both (opposing) muscle groups when only one muscle group was fatigued. Because the peak responses following the perturbation were enhanced in the unknown timing condition, preparations must have taken place prior to the anticipated perturbations, perhaps in other segments of the body that were not measured. Also, the load impact may not have been great enough to elicit large preparations. The heightened baseline activity with fatigue suggests that there may have been increased spinal stiffness whenever the spine was fatigued, and not just immediately prior to an impending perturbation. The increased activation of opposing muscle groups is evidence of increased cocontraction in response to fatigue, possibly to maintain stability with decreasing coordination.     

Effect of instruction,surface stability,and load intensity on trunk muscle activity

The aim of this study was to assess the effect of verbal instruction, surface stability, and load intensity on trunk muscle activity levels during the free weight squat exercise. Twelve trained males performed a free weight squat under four conditions: (1) standing on stable ground lifting 50% of their 1-repetition maximum (RM), (2) standing on a BOSU balance trainer lifting 50% of their 1-RM, (3) standing on stable ground lifting 75% of their 1-RM, and (4) receiving verbal instructions to activate the trunk muscles followed by lifting 50% of their 1-RM. Surface EMG activity from muscles rectus abdominis (RA), external oblique (EO), transversus abdominis/internal oblique (TA/IO), and erector spinae (ES) were recorded for each condition and normalized for comparisons. Muscles RA, EO, and TA/IO displayed greater peak activity (39–167%) during squats with instructions compared to the other squat conditions (P = 0.04–0.007). Peak EMG activity of muscle ES was greater for the 75% 1-RM condition than squats with instructions or lifting 50% of 1-RM (P = 0.04–0.02). The results indicate that if the goal is to enhance EMG activity of the abdominal muscles during a multi-joint squat exercise then verbal instructions may be more effective than increasing load intensity or lifting on an unstable surface. However, in light of other research, conscious co-activation of the trunk muscles during the squat exercise may lead to spinal instability and hazardous compression forces in the lumbar spine.     

Interphalangeal joint and tendon forces: normal model and biomechanical consequences of surgical reconstruction

Soft tissue reconstructive surgery for rheumatoid-related proximal interphalangeal joint deformities frequently fails to produce the long-term predicted results. Detailed information on the biomechanics of this joint, under both normal and pathological conditions, is required to assess the efficacy of such surgical intervention. A biomechanical model of the proximal interphalangeal joint has been developed to investigate tendon and joint loading during real life three-dimensional activities. Based on a rigid body mechanics approach, the model uses high resolution MRI scans to obtain anatomical tendon and bone geometries in conjunction with three-dimensional kinematic and loading data. The model incorporates an optimisation routine which minimises overall maximum tendon stress in the eight individual elements considered. Radial and ulnar joint force components are included at the proximal interphalangeal joint level. Two simulated pathological versions of the mathematical model are developed to accommodate the altered anatomic relationships after tendon reconstructive surgery. Joint forces of up to 450N and common usage of the extensor mechanism during normal pinching and grasping activities are predicted. The ulnar lateral bands of the extensor tendon are generally loaded to a greater extent than the radial bands. Extensor tendon and joint forces in the simulated pathological models are significantly higher than those in the normal model. Combined with the poor tendon quality of rheumatoid arthritis patients generally, these amplified internal forces may lead to further joint deformation.     

Controlling for out-of-plane lumbar moments during unidirectional trunk efforts: Learning and reliability issues related to trunk muscle activation estimates

To assess the electromyographic (EMG) activation of trunk muscle during exertions performed in one primary plane (sagittal, frontal, transverse), we previously proposed a protocol allowing minimizing out-of-plane efforts (coupled moments – CMs) with the use of a static dynamometer combined with a visual feedback system. The aims of this study were to go further by testing motor learning and reliability issues related to such a protocol. Three identical sessions were conducted, where maximal voluntary contractions and submaximal ramp contractions were performed in six different directions while standing in the dynamometer. Two feedback conditions were tested, the simple 1D-feedback in the primary plane and the full 3D-feedback in all planes simultaneously. Surface EMG signals were collected from back and abdominal muscles and EMG amplitude and CMs were computed during the ramp contractions. Providing a 3D feedback to minimize CMs did not improve EMG reliability or in other words, did not reduce the within-subject variability. Providing three assessment days had practically no effect (no learning) on CMs and EMG variables. Overall, the reliability of EMG was at best moderate. However, although this limits its use on an individual basis, it still allows within- and between-group comparisons for research applications.     

A comparative study of two trunk biomechanical models under symmetric and asymmetric loadings

Despite recent advances in modeling of the human spine, simplifying assumptions are still required to tackle complexities. Such assumptions need to be scrutinized to assess their likely impacts on predictions. A comprehensive comparison of muscle forces and spinal loads estimated by a single-joint (L5–S1) optimisation-assisted EMG-driven (EMGAO) and a multi-joint Kinematics-driven (KD) model of the spine under symmetric (symmetric trunk flexion from neutral upright to maximum forward flexion) and asymmetric (holding a load at various heights in the right hand) activities is carried out. Regardless of the task simulated, the KD model predicted greater activities in extensor muscles as compared to the EMGAO model. Such differences in the symmetric tasks was due mainly to the distinct approaches to resolve the redundancy while in the asymmetric tasks they were due also to the different methods used to estimate joint moments. Shear and compression forces were generally higher in the KD model. Differences in predictions between these modeling approaches varied depending on the task simulated and the joint considered in the single-joint EMGAO model. The EMGAO model should incorporate a multi-joint strategy to satisfy equilibrium at different levels while the KD model should benefit from recorded EMG activities of the antagonistic muscles to supplement input measured kinematics.     

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.     

Influence of the 3D inverse dynamic method on the joint forces and moments during gait

The joint forces and moments are commonly used in gait analysis. They can be computed by four different 3D inverse dynamic methods proposed in the literature, either based on vectors and Euler angles, wrenches and quaternions, homogeneous matrices, or generalized coordinates and forces. In order to analyze the influence of the inverse dynamic method, the joint forces and moments were computed during gait on nine healthy subjects. A ratio was computed between the relative dispersions (due to the method) and the absolute amplitudes of the gait curves. The influence of the inverse dynamic method was negligible at the ankle (2%) but major at the knee and the hip joints (40%). This influence seems to be due to the dynamic computation rather than the kinematic computation. Compared to the influence of the joint center location, the body segment inertial parameter estimation, and more, the influence of the inverse dynamic method is at least of equivalent importance. This point should be confirmed with other subjects, possibly pathologic, and other movements.     

Effects of external trunk loads on lumbar spine stability

Stability of the lumbar spine is an important factor in determining spinal response to sudden loading. Using two different methods, this study evaluated how various trunk load magnitudes and directions affect lumbar spine stability. The first method was a quick release procedure in which effective trunk stiffness and stability were calculated from trunk kinematic response to a resisted-force release. The second method combined trunk muscle EMG data with a biomechanical model to calculate lumbar spine stability. Twelve subjects were tested in trunk flexion, extension, and lateral bending under nine permutations of vertical and horizontal trunk loading. The vertical load values were set at 0, 20, and 40% of the subject's body weight (BW). The horizontal loads were 0, 10, and 20% of BW. Effective spine stability as obtained from quick release experimentation increased significantly (p<0.01) with increased vertical and horizontal loading. It ranged from 785 (S.D.=580) Nm/rad under no-load conditions to 2200 (S.D.=1015) Nm/rad when the maximum horizontal and vertical loads were applied to the trunk simultaneously. Stability of the lumbar spine achieved prior to force release and estimated from the biomechanical model explained approximately 50% of variance in the effective spine stability obtained from quick release trials in extension and lateral bending (0.53相似文献   

A biomechanical analysis of finger joint forces and stresses developed during common daily activities

The problem of modelling stresses incurred at the finger joints is critical to the design of durable joint replacements in the hand. The goal of this study was to characterise the forces and stresses at the finger and thumb joints occurring during activities such as typing at a keyboard, playing piano, gripping a pen, carrying a weight and opening a jar. The metacarpal and proximal phalanx were modelled using a COMSOL-based finite element analysis. Analysis of these activities indicates that joint forces in excess of 100 N may be common at the metacarpophalangeal joint (MCP) due to carrying objects such as groceries or while opening jars. The model predicted that stresses in excess of 2 MPa, similar to stresses at the hip, occur at the MCP with the properties of cancellous bone playing a significant role in the magnitude and distribution of stress.