Surface EMG-force modelling for the biceps brachii and its experimental evaluation during isometric isotonic contractions

The aim of this study was to evaluate a surface electromyography (sEMG) signal and force model for the biceps brachii muscle during isotonic isometric contractions for an experimental set-up as well as for a simulation. The proposed model includes a new rate coding scheme and a new analytical formulation of the muscle force generation. The proposed rate coding scheme supposes varying minimum and peak firing frequencies according to motor unit (MU) type (I or II). Practically, the proposed analytical mechanogram allows us to tune the force contribution of each active MU according to its type and instantaneous firing rate. A subsequent sensitivity analysis using a Monte Carlo simulation allows deducing optimised input parameter ranges that guarantee a realistic behaviour of the proposed model according to two existing criteria and an additional one. In fact, this proposed new criterion evaluates the force generation efficiency according to neural intent. Experiments and simulations, at varying force levels and using the optimised parameter ranges, were performed to evaluate the proposed model. As a result, our study showed that the proposed sEMG–force modelling can emulate the biceps brachii behaviour during isotonic isometric contractions.     

A Finite Element Model Approach to Determine the Influence of Electrode Design and Muscle Architecture on Myoelectric Signal Properties

Introduction

Surface electromyography (sEMG) is the measurement of the electrical activity of the skeletal muscle tissue detected at the skin’s surface. Typically, a bipolar electrode configuration is used. Most muscles have pennate and/or curved fibres, meaning it is not always feasible to align the bipolar electrodes along the fibres direction. Hence, there is a need to explore how different electrode designs can affect sEMG measurements.

Method

A three layer finite element (skin, fat, muscle) muscle model was used to explore different electrode designs. The implemented model used as source signal an experimentally recorded intramuscular EMG taken from the biceps brachii muscle of one healthy male. A wavelet based intensity analysis of the simulated sEMG signal was performed to analyze the power of the signal in the time and frequency domain.

Results

The model showed muscle tissue causing a bandwidth reduction (to 20-92- Hz). The inter-electrode distance (IED) and the electrode orientation relative to the fibres affected the total power but not the frequency filtering response. The effect of significant misalignment between the electrodes and the fibres (60°- 90°) could be reduced by increasing the IED (25–30 mm), which attenuates signal cancellation. When modelling pennated fibres, the muscle tissue started to act as a low pass filter. The effect of different IED seems to be enhanced in the pennated model, while the filtering response is changed considerably only when the electrodes are close to the signal termination within the model. For pennation angle greater than 20°, more than 50% of the source signal was attenuated, which can be compensated by increasing the IED to 25 mm.

Conclusion

Differences in tissue filtering properties, shown in our model, indicates that different electrode designs should be considered for muscle with different geometric properties (i.e. pennated muscles).     

静态负荷诱发肌肉疲劳后恢复期sEMG信号变化规律

目的:探讨肌肉疲劳过程中sEMG功率谱变化与H 的关系以及可能存在的其它影响因素.方法:利用肌肉进行疲劳收缩结束后,短时间内肌肉pH值尚无明显改变的特性,观察恢复期30 s内s EMG功率谱的变化规律.八名男性受试者,以肱二头肌为目标肌肉,负荷强度为60%MVC,静态持续负荷至疲劳点后,在恢复期以同样负荷分别观察2 s、4 s、6 s、8 s、10 s、20 s、30 s时的sEMG信号特征.结果:肱二头肌在以60%MVC静态疲劳负荷过程中MPF呈线性下降.在疲劳负荷后的恢复期,MPF恢复极其迅速,运动结束后仅2 s,MPF已恢复到整个下降范围的26.5%;至30 s,MPF已恢复到整个下降范围的87.7%.结论:由[H ]增加引起的肌纤维动作电位传导速度下降不是决定sEMG功率谱左移的唯一因素,提示sEMG功率谱左移可能与神经源性的中枢机制的作用有关.     

Evaluation of upper limb muscle fatigue based on surface electromyography

Fatigue is believed to be a major contributory factor to occupational injuries in machine operators. The development of accurate and usable techniques to measure operator fatigue is therefore important. In this study, we used a novel method based on surface electromyography (sEMG) of the biceps brachii and the Borg scale to evaluate local muscle fatigue in the upper limb after isometric muscle action. Thirteen young males performed isometric actions with the upper limb at different force levels. sEMG activities of the biceps brachii were recorded during the actions. Borg scales were used to evaluate the subjective sensation of local fatigue of the biceps brachii after the actions. sEMG activities were analyzed using the one-third band octave method, and an equation to determine the degree of fatigue was derived based on the relationship between the variable and the Borg scale. The results showed that the relationship could be expressed by a conic curve, and could be used to evaluate muscle fatigue during machine operation.     

High frequency vibration conditioning stimulation centrally reduces myoelectrical manifestation of fatigue in healthy subjects

ObjectiveVibration conditioning has been adopted as a tool to improve muscle force and reduce fatigue onset in various rehabilitation settings. This study was designed to asses if high frequency vibration can induce some conditioning effects detectable in surface EMG (sEMG) signal; and whether these effects are central or peripheral in origin.Design300 Hz vibration was applied for 30 min during 5 consecutive days, to the right biceps brachii muscle of 10 healthy males aged from 25 to 50 years. sEMG was recorded with a 16 electrode linear array placed on the skin overlying the vibrated muscle. The test protocol consisted of 30% and 60% maximal voluntary contraction (MVC) as well as involuntary (electrically elicited) contractions before and after treatment.ResultsNo statistically significant differences were found between PRE and POST vibration conditioning when involuntary stimulus-evoked contraction and 30% MVC were used. Significant differences in the initial values and rates of change of muscle fibre conduction velocity were found only at 60% MVC.Conclusions300 Hz vibration did not induce any peripheral changes as demonstrated by the lack of differences when fatigue was electrically induced. Differences were found only when the muscle was voluntarily fatigued at 60% MVC suggesting a modification in the centrally driven motor unit recruitment order, and interpreted as an adaptive response to the reiteration of the vibratory conditioning.     

局部肌肉疲劳的表面肌电信号复杂度和熵变化

目的 在于探讨静态和动态疲劳性运动过程中肱二头肌和腰部脊竖肌表面肌电(surface electromyography,sEMG)信号的Lempel-Ziv复杂度和Kolmogorov熵的变化规律。18名男性大学生志愿者被随机分为肱二头肌和腰部脊竖肌运动负荷组,分别完成静态和动态疲劳运动负荷试验。运动负荷期间连续记录sEMG信号,在对运动负荷时间和重复次数进行标准化处理后,截取相应时段的sEMG信号,计算Lempel-Ziv复杂度和Kolmogorov熵,观察它们随肌肉疲劳发展的变化规律。研究结果表明,无论是静态还是动态疲劳运动条件下,被检肌肉sEMG信号的复杂度和熵均随着运动负荷时间呈现明显的单调递减型变化。该变化可能与神经系统渐进性协调众多运动单位同步收缩的‘协同效应”有关。     

Comparison of surface electromyography and myotonometric measurements during voluntary isometric contractions

Objectives: Muscle stiffness increases during muscle contraction. The purpose of this study was to determine the strength of the correlation between myotonometric measurements of muscle stiffness and surface electromyography (sEMG) measurements during various levels of voluntary isometric contractions of the biceps brachii muscle. Subjects: Eight subjects (four female; four male), with mean age of 30.6±8.23 years, volunteered to participate in this study. Methods: Myotonometer and sEMG measurements were taken simultaneously from the right biceps brachii muscle. Data were obtained: (1) at rest, (2) while the subject held a 15 lb (6.8 kg) weight isometrically and, (3) during a maximal voluntary isometric contraction. Myotonometer force–displacement curves (amount of tissue displacement to a given unit of force applied perpendicular to the muscle) were compared with sEMG measurements using Pearson’s product–moment correlation coefficients. Results: Myotonometer and sEMG measurement correlations ranged from −0.70 to −0.90. The strongest correlations to sEMG were from Myotonometer force measurements between 1.00 and 2.00 kg. Conclusions: Myotonometer and sEMG measurements were highly correlated. Tissue stiffness, as measured by the Myotonometer, appears capable of assessing changes in muscle activation levels.     

Effects of mechanical assistance on muscle activity and motor performance during isometric elbow flexion

Mechanical assistance on joint movement is generally beneficial; however, its effects on cooperative performance and muscle activity needs to be further explored. This study examined how motor performance and muscle activity are altered when mechanical assistance is provided during isometric force control of ramp-down and hold phases. Thirteen right-handed participants (age: 24.7 ± 1.8 years) performed trajectory tracking tasks. Participants were asked to maintain the reference magnitude of 47 N (REF) during isometric elbow flexion. The force was released to a step-down magnitude of either 75% REF or 50% REF and maintained, with and without mechanical assistance. The ramp-down durations of force release were set to 0.5, 2.5, or 5.0 s. Throughout the experiment, we measured the following: (1) the force output using load cells to compute force variability and overshoot ratio; (2) peak perturbation on the elbow movement using an accelerometer; (3) the surface electromyography (sEMG) from biceps brachii and triceps brachii muscles; and (4) EMG oscillation from the biceps brachii muscle in the bandwidth of 15–45 Hz. Our results indicated that mechanical assistance, which involved greater peak perturbation, demonstrated lower force variability than non-assistance (p < 0.01), while EMG oscillation in the biceps brachii muscle from 15 to 45 Hz was increased (p < 0.05). These findings imply that if assistive force is provided during isometric force control, the central nervous system actively tries to stabilize motor performance by controlling specific motor unit activity in the agonist muscle.     

A Wavelet-Based Method to Predict Muscle Forces From Surface Electromyography Signals in Weightlifting

The purpose of this study was to develop a wavelet-based method to predict muscle forces from surface electromyography (EMG) signals in vivo.The weightlifting motor task was implemented as the case study.EMG signals of biceps brachii,triceps brachii and deltoid muscles were recorded when the subject carried out a standard weightlifting motor task.The wavelet-based algorithm was used to process raw EMG signals and extract features which could be input to the Hill-type muscle force models to predict muscle forces.At the same time,the musculoskeletal model of subject's weightlifting motor task was built and simulated using the Computed Muscle Control (CMC) method via a motion capture experiment.The results of CMC were compared with the muscle force predictions by the proposed method.The correlation coefficient between two results was 0.99(p＜0.01).However,the proposed method was easier and more efficiency than the CMC method.It has potential to be used clinically to predict muscle forces in vivo.     

Propagation direction of natural mechanical oscillations in the biceps brachii muscle during voluntary contraction

The aim of the study was to determine the directionality of the coupling of mechanical vibrations across the biceps brachii muscle at different frequencies of interest during voluntary contraction. The vibrations that are naturally generated by skeletal muscles were recorded by a two-dimensional array of skin mounted accelerometers over the biceps brachii muscle (surface mechanomyogram, S-MMG) during voluntary isometric contractions in ten healthy young men. As a measure of the similarity of vibration between a given pair of accelerometers, the spatial coherence of S-MMG at low (f < 25 Hz) and high (f > 25 Hz) frequency bands were investigated to determine if the coupling of the natural mechanical vibrations were due to the different physiological muscle activity at low and high frequencies. In both frequency bands, spatial coherence values for sensor pairs aligned longitudinally along the proximal to distal ends of the biceps were significantly higher compared with those for the sensor pairs oriented perpendicular to the muscle fibers. This difference was more evident at the higher frequency band. The findings indicated that coherent mechanical oscillations mainly propagated along the longitudinal direction of the biceps brachii muscle fibers at high frequencies (f > 25 Hz).     

Elimination of electrocardiogram contamination from electromyogram signals: An evaluation of currently used removal techniques.

Trunk electromyographic signals (EMG) are often contaminated with heart muscle electrical activity (ECG) due to the proximity of the collection sites to the heart and the volume conduction characteristics of the ECG through the torso. Few studies have quantified ECG removal techniques relative to an uncontaminated EMG signal (gold standard or criterion measure), or made direct comparisons between different methods for a given set of data. Understanding the impacts of both untreated contaminated EMG and ECG elimination techniques on the amplitude and frequency parameters is vital given the widespread use of EMG. The purpose of this study was to evaluate four groups of current and commonly used techniques for the removal of ECG contamination from EMG signals. ECG recordings at two intensity levels (rest and 50% maximum predicted heart rate) were superimposed on 11 uncontaminated biceps brachii EMG signals (rest, 7 isometric and 3 isoinertial levels). The 23 removal methods used were high pass digital filtering (finite impulse response (FIR) using a Hamming window, and fourth-order Butterworth (BW) filter) at five cutoff frequencies (20, 30, 40, 50, and 60 Hz), template techniques (template subtraction and an amplitude gating template), combinations of the subtraction template and high pass digital filtering, and a frequency subtraction/signal reconstruction method. For muscle activation levels between 10% and 25% of maximum voluntary contraction, the template subtraction and BW filter with a 30 Hz cutoff were the two best methods for maximal ECG removal with minimal EMG distortion. The BW filter with a 30 Hz cutoff provided the optimal balance between ease of implementation, time investment, and performance across all contractions and heart rate levels for the EMG levels evaluated in this study.     

Facilitation of triceps brachii muscle contraction by tendon vibration after chronic cervical spinal cord injury.

One way to improve the weak triceps brachii voluntary forces of people with chronic cervical spinal cord injury may be to excite the paralyzed or submaximally activated fraction of muscle. Here we examined whether elbow extensor force was enhanced by vibration (80 Hz) of the triceps or biceps brachii tendons at rest and during maximum isometric voluntary contractions (MVCs) of the elbow extensors performed by spinal cord-injured subjects. The mean +/- SE elbow extensor MVC force was 22 +/- 17.5 N (range: 0-23% control force, n = 11 muscles). Supramaximal radial nerve stimuli delivered during elbow extensor MVCs evoked force in six muscles that could be stimulated selectively, suggesting potential for force improvement. Biceps vibration at rest always evoked a tonic vibration reflex in biceps, but extension force did not improve with biceps vibration during triceps MVCs. Triceps vibration induced a tonic vibration reflex at rest in one-half of the triceps muscles tested. Elbow extensor MVC force (when >1% of control force) was enhanced by vibration of the triceps tendon in one-half of the muscles. Thus triceps, but not biceps, brachii tendon vibration increases the contraction strength of some partially paralyzed triceps brachii muscles.     

Exercise induced muscle damage and recovery assessed by means of linear and non-linear sEMG analysis and ultrasonography.

This study was aimed at investigating the time-course and recovery from eccentric (EC) exercise induced muscle damage by means of surface electromyography (sEMG), ultrasonography (US), and blood enzymes. Five subjects (EC Group) performed two bouts of 35 EC maximum contractions with the biceps brachii of their non dominant arm, five subjects were tested without performing EC (Control Group: CNT). The maximal isometric force (MVC) was measured. Force and sEMG signals were recorded during 80% MVC isometric contractions. In EC and CNT subjects US assessment on non-dominant biceps brachii was performed; creatin kinase (CK) and lactic dehydrogenasis (LDH) plasma levels were also assessed. Force, sEMG and CK-LDH measurements were performed before EC and after it periodically for 4 weeks. The sEMG was analysed in time and frequency domains; a non-linear analysis (Lyapunov 1st exponent, L1) of sEMG was also performed. After EC, the MVC was reduced by 40% on average with respect to the pre-EC values. A significant decrease in the initial frequency content, and in the MDF and L1 decay (13-42% less than the pre-EC values, respectively) was also observed. The sEMG amplitude (Root Mean Square, RMS) was unchanged after EC. The US revealed an increase in muscle belly thickness and in local muscle blood flow after EC. A complete recovery of all the considered parameters was achieved in two weeks. In conclusion sEMG analysis was confirmed as an early indicator of muscle damage. Muscle recovery from damage is followed by both sEMG and US and this may have useful clinical implications. Non linear analysis (L1) was revealed to be sensitive to early sEMG modifications induced by EC as well as able to follow the post EC changes in the sEMG.     

Whole-body vibration increases upper and lower body muscle activity in older adults: potential use of vibration accessories

The current study examined the effects of whole-body vibration (WBV) on upper and lower body muscle activity during static muscle contractions (squat and bicep curls). The use of WBV accessories such as hand straps attached to the platform and a soft surface mat were also evaluated. Surface electromyography (sEMG) was measured for the medial gastrocnemius (MG), vastus lateralis (VL), and biceps brachii (BB) muscles in fourteen healthy older adults (74.8±4.5 years; mean±SD) with a WBV stimulus at an acceleration of 40 m s(-2) (30 Hz High, 2.5 mm or 46 Hz Low, 1.1 mm). WBV increased lower body (VL and MG) sEMG vs baseline (no WBV) though this was decreased with the use of the soft mat. The addition of the bicep curl with hand straps had no effect on lower body sEMG. WBV also increased BB sEMG vs baseline which was further increased when using the hand straps. There was no upper body effect of the soft mat. This study demonstrates WBV increases both lower and upper body muscle activity in healthy older adults. Moreover, WBV accessories such as hand straps attached to the platform or a soft surface mat may be used to alter exercise intensity.     

Neuromuscular fatigue during maximal concurrent hand grip and elbow flexion or extension.

The objective was to investigate muscle fatigue measuring changes in force output and force tremor and electromyographic activity (EMG) during two sustained maximal isometric contractions for 60s: (1) concurrent hand grip and elbow flexion (HG and EF); or (2) hand grip and elbow extension (HG and EE). Each force tremor amplitude was decomposed into four frequency bands (1-3, 4-10, 11-20, and 21-50Hz). Surface EMGs were recorded from the flexor digitorum superficialis (FDS), extensor digitorum (ED), biceps brachii (BB) and lateral head of triceps brachii (TB). The HG and EF forces for the HG and EF and the HG force for the HG and EE declined rapidly, whereas the EE force remained almost constant near to the initial value for the first 40s and then declined. The decrease in EMG amplitude was observed not for the FDS muscle but for the ED muscle. The HG tremor amplitude for each frequency band showed similar decreasing rate, whereas the decreases in EF and EE tremor amplitudes for the lower band (below 10Hz) were slower than those for the higher band (above 11Hz). The neuromuscular mechanisms underlying muscle fatigue during sustained maximal concurrent contractions of hand grip and elbow flexion or extension are discussed.     

Spectral analysis of the surface electromyogram as a tool for studying rate modulation: A comparison between theory,simulation, and experiment

Theoretical work suggests that if the interpulse intervals (IPIs) of motor unit action potential trains (MUAPTs) are independently and normally distributed, then spectral analysis of the electromyogram could be a useful tool for studying rate modulation by virtue of the presence of a peak in the power spectrum at the average firing frequency of all active motor units. It is shown in this paper that IPIs need not be normally distributed, specifically that the results are very much the same if the IPIs are distributed according to a Gamma probability density function (PDF). Simulation of the electromyogram based on this theory proved the applicability of the method. Experimental results obtained for the masseter, biceps brachii and first dorsal interosseus (FDI) muscles, however, were in disagreement with both theory and simulation except for the biceps muscle at force levels up to 20% of the maximal force and for the masseter and FDI muscles in 1 out of 5 subjects. This indicates that the models for MUAPTs hitherto used might not be generally correct. Apart from this discrepancy, our results reveal differences between masseter and FDI muscles on the one hand and the biceps brachii on the other, which indicate that motor unit synchronisation is much more pronounced in the latter muscle.     

Distribution of motor unit potential velocities in the biceps brachii muscle of sprinters and endurance athletes during short static contractions at low force levels

In surface electromyography (sEMG), the distribution of motor unit potential (MUP) velocities has been shown to reflect the proportion of faster and slower propagating MUPs. This study investigated whether the distribution of MUP velocities could distinguish between sprinters and endurance athletes in not-specifically trained muscle (biceps brachii). sEMG results were acquired from 15 sprinters and 18 endurance athletes during short static contractions (3.8 s) at three force levels: unloaded, 10% and 20% of maximum voluntary contraction. The features extracted from the sEMG were: the mean muscle conduction velocity (CV) – estimated using the inter-peak latency and the cross-correlation methods, the within-subject skewness of MUP velocities (expressing the relative proportions of faster and slower propagating MUPs), and the within-subject standard deviation of MUP velocities. Sprinters had a higher CV than endurance athletes using both methods. Sprinters also demonstrated a greater proportion of fast propagating MUPs, as indicated by the skewness. Thus, the distribution of MUP velocities was able to demonstrate physiological differences between sprinters and endurance athletes during short contractions at low forces. The findings can be extrapolated to the motor unit level. Since the investigated muscle was not involved in specific training, the differences seem to reflect inherited properties.     

Repeatability of surface EMG parameters at various isometric contraction levels and during fatigue using bipolar and Laplacian electrode configurations.

The objective of this work was to assess the repeatability of two surface electromyographic (sEMG) recording techniques, the classical bipolar configuration and a Laplacian configuration to document their ability to provide reliable information during follow-up studies. The signals were recorded on 10 healthy subjects during voluntary isometric contractions of the biceps brachii muscle at different constant contraction levels. Slopes, area ratios (at 60% of the maximal voluntary contraction (MVC)) and initial values (at 20%, 40%, 60%, 80% and 100% MVC) of the root mean square (RMS), the mean power frequency (MPF) and the muscle fibre conduction velocity (CV) were estimated. Experimental sessions were repeated on three different days with both electrode sets to evaluate the repeatability of sEMG parameter estimates. Classical results were observed, such as an increase in the RMS and the CV with the contraction level. Only initial values of RMS and MPF were shown to be dependent on electrode type. These two parameters presented intra-class correlation coefficient values higher than .80 for high contraction levels. On the whole, the repeatability of the measures was good; however it was better for all sEMG parameter estimates with bipolar electrodes than Laplacian electrodes. Because a bipolar configuration is less selective than a Laplacian one, it provides a global view of muscular activity, which is more repeatable, hence more suitable for follow-up studies.     

Force-length, torque-angle and EMG-joint angle relationships of the human in vivo biceps brachii

The relationships of EMG and muscle force with elbow joint angle were investigated for muscle modelling purposes. Eight subjects had their arms fixed in an isometric elbow jig where the biceps brachii was electrically stimulated (30 Hz) and also in maximum voluntary contraction (MVC). Biceps EMG and elbow torque transduced at the wrist were recorded at 0.175 rad intervals through 1.75 rad of elbow extension. The results revealed that while the torque-length relationship displayed the classic inverted U pattern in both evoked and MVC conditions, the force-length relationship displayed a monotonically increasing pattern. Analyses of variance of the EMG data showed that there were no significant changes in the EMG amplitudes for the different joint angles during evoked or voluntary contractions. The result also showed that electrical stimulation can effectively isolated the torque-angle and force-length relationships of the biceps brachii and that the myoelectric signal during isometric contraction is uniform regardless of the length of the muscle or the joint angle.     

Effects of EMG processing on biomechanical models of muscle joint systems: sensitivity of trunk muscle moments, spinal forces, and stability

Biomechanical models are in use to estimate parameters such as contact forces and stability at various joints. In one class of these models, surface electromyography (EMG) is used to address the problem of mechanical indeterminacy such that individual muscle activation patterns are accounted for. Unfortunately, because of the stochastical properties of EMG signals, EMG based estimates of muscle force suffer from substantial estimation errors. Recent studies have shown that improvements in muscle force estimation can be achieved through adequate EMG processing, specifically whitening and high-pass (HP) filtering of the signals. The aim of this paper is to determine the effect of such processing on outcomes of a biomechanical model of the lumbosacral joint and surrounding musculature. Goodness of fit of estimated muscle moments to net moments and also estimated joint stability significantly increased with increasing cut-off frequencies in HP filtering, whereas no effect on joint contact forces was found. Whitening resulted in moment estimations comparable to those obtained from optimal HP filtering with cut-off frequencies over 250 Hz. Moreover, compared to HP filtering, whitening led to a further increase in estimated joint-stability. Based on theoretical models and on our experimental results, we hypothesize that the processing leads to an increase in pick-up area. This then would explain the improvements from a better balance between deep and superficial motor unit contributions to the signal.