Modelling torque generation by the mero-carpopodite joint of the american lobster and the snow crab

The torque generated by a rotating joint comprises the useful force exerted by the joint on the external environment, and both the magnitude and distribution of torque through the step cycle during walking are important variables in understanding the mechanics of walking. The mechanics of the American lobster (Homarus americanus) and snow crab (Chionoecetes opilio) during walking were modelled to examine the relative roles of flexor versus extensor apodeme-muscle complexes, investigate which legs of these decapods likely contribute the greatest to locomotion, determine scaling effects of torque generation, and assess the relative roles of various model variables on torque production. Force generated along the length of the apodeme by the muscle was modelled based on apodeme surface area, muscle stress, and muscle fibre pinnation angle. Torque was then calculated from this estimated force and the corresponding moment arm. The flexor apodeme-muscle complex is calculated to generate consistently greater forces than the extensor, and generally this results in flexor torque being larger than extensor, though the snow crab does illustrate the opposite in two of its legs. This greater torque generation in flexion suggests that, in addition to the pushing of the trailing legs, the pulling action of the leading legs may play a significant role, at least during lateral walking. Leg 4 of both species appears to generate greater torques and thus provide the greatest forces for locomotion. Torque generation as a function of body size shows a second order response due to the increase in apodeme surface area. The pinnation angle of the muscle fibre is found to be insignificant in force generation, apodeme surface area (representing muscle cross sectional area) likely plays the most influential role in total force production, and moment arm controls the distribution of this force through the step cycle. Muscle stress remain a largely unknown quantity however, and may significantly affect both magnitude and distribution through step cycle of forces, and thus torque. Despite the uncertainty associated with the muscle stress parameter, the modelled results fit well with previously published force measurements.     

A dynamic simulator to evaluate distal radio-ulnar joint kinematics

In order to perform cadaveric biomechanical studies of the human forearm and distal radio-ulnar joint, a dynamic simulator has been constructed. The device is based upon a Plexiglas frame, to which the ulna is secured in a vertical orientation and the humerus in a horizontal orientation. The hand is secured in a sliding bar linkage to a stepper-motor that is used to rotate the forearm. The tendons to be loaded are connected to pneumatic actuators that provide agonist and antagonist muscle loading resulting in torque along the forearm axis. The muscle loading profiles and magnitudes are programmable as a function of the pronation-supination position and direction. A magnetic tracking system is used to collect three-dimensional kinematics data of up to four segments, in conjunction with the muscle tendon loads, forearm torque and other prescribed experimental measures. All functions are under PC control using custom software written with LabVIEW (National Instruments, Austin, TX). For the DRUJ testing, the validity of the tendon loading protocol to produce physiologic torque/rotation patterns was verified using in vivo data. The relationship of individual muscle forces to forearm torque was determined by a cadaveric study.     

Task Dependency of Grip Stiffness—A Study of Human Grip Force and Grip Stiffness Dependency during Two Different Tasks with Same Grip Forces

It is widely known that the pinch-grip forces of the human hand are linearly related to the weight of the grasped object. Less is known about the relationship between grip force and grip stiffness. We set out to determine variations to these dependencies in different tasks with and without visual feedback. In two different settings, subjects were asked to (a) grasp and hold a stiffness-measuring manipulandum with a predefined grip force, differing from experiment to experiment, or (b) grasp and hold this manipulandum of which we varied the weight between trials in a more natural task. Both situations led to grip forces in comparable ranges. As the measured grip stiffness is the result of muscle and tendon properties, and since muscle/tendon stiffness increases more-or-less linearly as a function of muscle force, we found, as might be predicted, a linear relationship between grip force and grip stiffness. However, the measured stiffness ranges and the increase of stiffness with grip force varied significantly between the two tasks. Furthermore, we found a strong correlation between regression slope and mean stiffness for the force task which we ascribe to a force stiffness curve going through the origin. Based on a biomechanical model, we attributed the difference between both tasks to changes in wrist configuration, rather than to changes in cocontraction. In a new set of experiments where we prevent the wrist from moving by fixing it and resting it on a pedestal, we found subjects exhibiting similar stiffness/force characteristics in both tasks.     

Evaluation by Expert Dancers of a Robot That Performs Partnered Stepping via Haptic Interaction

Our long-term goal is to enable a robot to engage in partner dance for use in rehabilitation therapy, assessment, diagnosis, and scientific investigations of two-person whole-body motor coordination. Partner dance has been shown to improve balance and gait in people with Parkinson''s disease and in older adults, which motivates our work. During partner dance, dance couples rely heavily on haptic interaction to convey motor intent such as speed and direction. In this paper, we investigate the potential for a wheeled mobile robot with a human-like upper-body to perform partnered stepping with people based on the forces applied to its end effectors. Blindfolded expert dancers (N=10) performed a forward/backward walking step to a recorded drum beat while holding the robot''s end effectors. We varied the admittance gain of the robot''s mobile base controller and the stiffness of the robot''s arms. The robot followed the participants with low lag (M=224, SD=194 ms) across all trials. High admittance gain and high arm stiffness conditions resulted in significantly improved performance with respect to subjective and objective measures. Biomechanical measures such as the human hand to human sternum distance, center-of-mass of leader to center-of-mass of follower (CoM-CoM) distance, and interaction forces correlated with the expert dancers'' subjective ratings of their interactions with the robot, which were internally consistent (Cronbach''s α=0.92). In response to a final questionnaire, 1/10 expert dancers strongly agreed, 5/10 agreed, and 1/10 disagreed with the statement "The robot was a good follower." 2/10 strongly agreed, 3/10 agreed, and 2/10 disagreed with the statement "The robot was fun to dance with." The remaining participants were neutral with respect to these two questions.     

Effect of muscle model parameter scaling for isometric plantar flexion torque prediction

This paper uses a EMG-driven Hill-type muscle model to estimate individual muscle forces of the triceps surae in isometric plantar flexion contractions. A uniform group of 20 young physical-active adult males was instructed to follow a specific contraction protocol with low (20%MVC) and medium-high (60%MVC) contractions, separated by relaxing intervals. The torque calculated by summing the individual muscle forces multiplied by the respective moment arms was compared to the torque measured by a dynamometer. Musculoskeletal parameters from the literature were used. Then, three different “correction factors” or bias have been applied on some of the muscle model parameters. These factors were based on anthropometric and dynamometric measurements: moment arm scaled by bimalleolar diameter, tendon slack length by leg length and optimal force by the maximum torque. Model torque agreement with dynamometer was recalculated with the parameter scales. It was observed that the relative torque estimation error decreased slightly but significantly when all factors were applied simultaneously (12.92±4.94% without scaling to 10.12±1.73%), which resulted mainly from the correction of the maximal muscle force parameter.     

Human plantarflexor stiffness to multiple single-stretch trials.

The purpose of this investigation was to determine the influence of different stretch velocities, different rates of pre-stretch force development, and different pre-stretch muscle lengths on the intrinsic stiffness exhibited by the quasi-statically contracting active human plantarflexors during multiple single-stretch trials at 20-60% of maximum isometric contraction. Subjects were positioned prone, with the knee flexed 1.57 rad(90 degrees), shank stabilized, and foot secured in a hard plastic orthotic. Slowly increasing isometric plantarflexion force was produced until the plantarflexors were stretched by a rapid 0.2 rad (12 degrees) dorsiflexion movement. Plantarflexion forces and ankle positions were determined during these stretches as well as during resting stretches when the muscle was inactive. Resting forces were subtracted from the active trials, forces converted to torques, and stiffnesses determined for the first 62 ms of the stretch. The slope of the stiffness vs pre-stretch torque relationship averaged 4.30 +/- 0.34 Nm rad-1 Nm-1. Little difference was found between stiffness determined through the single-stretch method and the results of previous studies employing different mechanical inputs. Differences in stiffnesses with different stretching velocities were caused by computational artifact rather than by differences in intrinsic muscular reaction. Faster rates of pre-stretch force increase prior to the stretch resulted in slightly lower stiffnesses. Different pre-stretch muscle lengths apparently did not result in different stiffnesses. The shape of the torque vs displacement curve was remarkably insensitive to the planned manipulations of the testing conditions, responding in a stereotypical manner.     

Evaluation of a subject-specific finite-element model of the equine metacarpophalangeal joint under physiological load

The equine metacarpophalangeal (MCP) joint is frequently injured, especially by racehorses in training. Most injuries result from repetitive loading of the subchondral bone and articular cartilage rather than from acute events. The likelihood of injury is multi-factorial but the magnitude of mechanical loading and the number of loading cycles are believed to play an important role. Therefore, an important step in understanding injury is to determine the distribution of load across the articular surface during normal locomotion. A subject-specific finite-element model of the MCP joint was developed (including deformable cartilage, elastic ligaments, muscle forces and rigid representations of bone), evaluated against measurements obtained from cadaver experiments, and then loaded using data from gait experiments. The sensitivity of the model to force inputs, cartilage stiffness, and cartilage geometry was studied. The FE model predicted MCP joint torque and sesamoid bone flexion angles within 5% of experimental measurements. Muscle–tendon forces, joint loads and cartilage stresses all increased as locomotion speed increased from walking to trotting and finally cantering. Perturbations to muscle–tendon forces resulted in small changes in articular cartilage stresses, whereas variations in joint torque, cartilage geometry and stiffness produced much larger effects. Non-subject-specific cartilage geometry changed the magnitude and distribution of pressure and the von Mises stress markedly. The mean and peak cartilage stresses generally increased with an increase in cartilage stiffness. Areas of peak stress correlated qualitatively with sites of common injury, suggesting that further modelling work may elucidate the types of loading that precede joint injury and may assist in the development of techniques for injury mitigation.     

Straightening Out the Elasticity of Myosin Cross-Bridges

In a contracting muscle, myosin cross-bridges extending from thick filaments pull the interdigitating thin (actin-containing) filaments during cyclical ATP-driven interactions toward the center of the sarcomere, the structural unit of striated muscle. Cross-bridge attachments in the sarcomere have been reported to exhibit a similar stiffness under both positive and negative forces. However, in vitro measurements on filaments with a sparse complement of heads detected a decrease of the cross-bridge stiffness at negative forces attributed to the buckling of the subfragment 2 tail portion. Here, we review some old and new data that confirm that cross-bridge stiffness is nearly linear in the muscle filament lattice. The implications of high myosin stiffness at positive and negative strains are considered in muscle fibers and in nonmuscle intracellular cargo transport.     

A mathematical analysis of the force-stiffness characteristics of muscles in control of a single joint system

Feldman (1966) has proposed that a muscle endowed with its spinal reflex system behaves as a non-linear spring with an adjustable resting length. In contrast, because of the length-tension properties of muscles, many researchers have modeled them as non-linear springs with adjustable stiffness. Here we test the merits of each approach: Initially, it is proven that the adjustable stiffness model predicts that isometric muscle force and stiffness are linearly related. We show that this prediction is not supported by data on the static stiffness-force characteristics of reflexive muscles, where stiffness grows non-linearly with force. Therefore, an intact muscle-reflex system does not behave as a non-linear spring with an adjustable stiffness. However, when the same muscle is devoid of its reflexes, the data shows that stiffness grows linearly with force. We aim to understand the functional advantage of the non-linear stiffness-force relationship present in the reflexive muscle. Control of an inverted pendulum with a pair of antagonist muscles is considered. Using an active-state muscle model we describe force development in an areflexive muscle. From the data on the relationship of stiffness and force in the intact muscle we derive the length-tension properties of a reflexive muscle. It is shown that a muscle under the control of its spinal reflexes resembles a non-linear spring with an adjustable resting length. This provides independent evidence in support of the Feldman hypothesis of an adjustable resting length as the control parameter of a reflexive muscle, but it disagrees with his particular formulation. In order to maintain stability of the single joint system, we prove that a necessary condition is that muscle stiffness must grow at least linearly with force at isometric conditions. This shows that co-contraction of antagonist muscles may actually destabilize the limb if the slope of this stiffness-force relationship is less than an amount specified by the change in the moment arm of the muscle as a function of joint configuration. In a reflexive muscle where stiffness grows faster than linearly with force, co-contraction will always lead to an increase in stiffness. Furthermore, with the reflexive muscles, the same level of joint stiffness can be produced by much smaller muscle forces because of the non-linear stiffness-force relationship. This allows the joint to remain stable at a fraction of the metabolic energy cost associated with maintaining stability with areflexive muscles.This work was supported in part by grant no. 1R01 NS 24926 from the NIH (Michael Arbib, PI). R.S. was supported by an IBM Graduate Fellowship in Computer Science     

Finger force vectors in multi-finger prehension

In a majority of studies on grasp, only normal forces were measured and only when a zero torque was exerted on a hand-held object. This study concerns finger force vectors during the torque production tasks. Subjects (n=8) stabilized a handle with an attachment that allowed for change of external torque from -1.5 to 1.5 Nm. Forces and moments exerted by the digit tips on the object were recorded. At the large (>-0.375 Nm) supination torques the index/middle and ring/little pairs of fingers generated oppositely directed tangential forces. The index and middle finger produced forces in a downward direction and therefore did not support the load. At a zero torque and pronation torques, the middle, ring and little fingers produced forces along nearly the same direction. The vector of the index finger force was always directed differently from the vectors of other finger forces, the angles ranged from 19 degrees 30' to 47 degrees 40'. The points of force application were systematically displaced with the torque, with the exception of the little finger. Tangential finger forces contributed substantially to the total torque exerted on the hand-held object.     

Comparison of global and joint-to-joint methods for estimating the hip joint load and the muscle forces during walking

A three-dimensional musculoskeletal model of the lower limb was developed to study the influence of biarticular muscles on the muscle force distribution and joint loads during walking. A complete walking cycle was recorded for 9 healthy subjects using the standard optoelectronic motion tracking system. Ground contact forces were also measured using a 6-axes force plate. Inverse dynamics was used to compute net joint reactions (forces and torques) in the lower limb. A static optimization method was then used to estimate muscle forces. Two different approaches were used: in the first one named global method, the biarticular muscles exerted a torque on the two joints they spanned at the same time, and in the second one called joint-by-joint method, these biarticular muscles were divided into two mono-articular muscles with geometrical (insertion, origin, via points) and physiological properties remained unchanged. The hip joint load during the gait cycle was then calculated taking into account the effect of muscle contractions. The two approaches resulted in different muscle force repartition: the biarticular muscles were favoured over any set of single-joint muscles with the same physiological function when using the global method. While the two approaches yielded only little difference in the resultant hip load, the examination of muscle power showed that biarticular muscles could produce positive work at one joint and negative work at the other, transferring energy between body segments and thus decreasing the metabolic cost of movement.     

Magnetic torque tweezers: measuring torsional stiffness in DNA and RecA-DNA filaments

We introduce magnetic torque tweezers, which enable direct single-molecule measurements of torque. Our measurements of the effective torsional stiffness C of dsDNA indicated a substantial force dependence, with C = approximately 40 nm at low forces up to C = approximately 100 nm at high forces. The initial torsional stiffness of RecA filaments was nearly twofold larger than that for dsDNA, yet at moderate torques further build-up of torsional strain was prevented.     

Spinal stiffness increases with axial load: another stabilizing consequence of muscle action.

This paper addresses the role of lumbar spinal motion segment stiffness in spinal stability. The stability of the lumbar spine was modelled with loadings of 30 Nm or 60 Nm efforts about each of the three principal axes, together with the partial body weight above the lumbar spine. Two assumptions about motion segment stiffness were made: first the stiffness was represented by an 'equivalent beam' with constant stiffness properties; second the stiffness was updated based on the motion segment axial loading using a relationship determined experimentally from human lumbar spinal specimens tested with 0, 250 and 500 N of axial compressive preload. Two physiologically plausible muscle activation strategies were used in turn for calculating the muscle forces required for equilibrium. Stability analyses provided estimates of the minimum muscle stiffness required for stability. These critical muscle stiffness values decreased when preload effects were used in estimating spinal stiffness in all cases of loadings and muscle activation strategies, indicating that stability increased. These analytical findings emphasize that the spinal stiffness (as well as muscular stiffness) is important in maintaining spinal stability, and that the stiffness-increasing effect of 'preloading' should be taken into account in stability analyses.     

Estimation of muscle forces in gait using a simulation of the electromyographic activity and numerical optimization

Clinical gait analysis provides great contributions to the understanding of gait patterns. However, a complete distribution of muscle forces throughout the gait cycle is a current challenge for many researchers. Two techniques are often used to estimate muscle forces: inverse dynamics with static optimization and computer muscle control that uses forward dynamics to minimize tracking. The first method often involves limitations due to changing muscle dynamics and possible signal artefacts that depend on day-to-day variation in the position of electromyographic (EMG) electrodes. Nevertheless, in clinical gait analysis, the method of inverse dynamics is a fundamental and commonly used computational procedure to calculate the force and torque reactions at various body joints. Our aim was to develop a generic musculoskeletal model that could be able to be applied in the clinical setting. The musculoskeletal model of the lower limb presents a simulation for the EMG data to address the common limitations of these techniques. This model presents a new point of view from the inverse dynamics used on clinical gait analysis, including the EMG information, and shows a similar performance to another model available in the OpenSim software. The main problem of these methods to achieve a correct muscle coordination is the lack of complete EMG data for all muscles modelled. We present a technique that simulates the EMG activity and presents a good correlation with the muscle forces throughout the gait cycle. Also, this method showed great similarities whit the real EMG data recorded from the subjects doing the same movement.     

Effects of force detecting sense organs on muscle synergies are correlated with their response properties

Sense organs that monitor forces in legs can contribute to activation of muscles as synergist groups. Previous studies in cockroaches and stick insects showed that campaniform sensilla, receptors that encode forces via exoskeletal strains, enhance muscle synergies in substrate grip. However synergist activation was mediated by different groups of receptors in cockroaches (trochanteral sensilla) and stick insects (femoral sensilla). The factors underlying the differential effects are unclear as the responses of femoral campaniform sensilla have not previously been characterized. The present study characterized the structure and response properties (via extracellular recording) of the femoral sensilla in both insects. The cockroach trochantero-femoral (TrF) joint is mobile and the joint membrane acts as an elastic antagonist to the reductor muscle. Cockroach femoral campaniform sensilla show weak discharges to forces in the coxo-trochanteral (CTr) joint plane (in which forces are generated by coxal muscles) but instead encode forces directed posteriorly (TrF joint plane). In stick insects, the TrF joint is fused and femoral campaniform sensilla discharge both to forces directed posteriorly and forces in the CTr joint plane. These findings support the idea that receptors that enhance synergies encode forces in the plane of action of leg muscles used in support and propulsion.     

Novel bioinspired control approaches to increase the stiffness variability in multi-muscle driven joints

Antagonistic muscle pairs pulling on a joint are in general able to modulate stiffness through co-activation. Closer analysis of the stiffness, however, shows that, depending on the muscle and joint parameters, domains might occur in joint angle space for which stiffness variation is limited (low stiffness variability) or even impossible (stiffness nodes). As a consequence, stiffness control utilizing pure co-activation might fail. This work presents novel strategies for simultaneous control of torque and stiffness in a hinge joint actuated by two antagonistic muscle pairs. One strategy handles stiffness nodes by shifting them away from the current joint position and thus regaining stiffness controllability. To prevent domains of low stiffness variation, an optimal muscle configuration is sought and finally defined which allows for a maximal stiffness variation across a wide joint angle range. Based on this optimal configuration, four additional control strategies are proposed and tested which deliver stiffnesses and torques comparable to those obtained in the optimal case. The strategies combine torque control and stiffness control by co-activation with novel ideas like activation overflow and an inverse model approach. All strategies are tested in simulation and the results are compared with those of the optimal setup.     

Impedance characteristics of a neuromusculoskeletal model of the human arm I. Posture control

The mechanical impedance of neuromusculoskeletal models of the human arm is studied in this paper. The model analysis provides a better understanding of the contributions of possible intrinsic and reflexive components of arm impedance, makes clear the limitations of second-order mass-viscosity-stiffness models and reveals possible task effects on the impedance. The musculoskeletal model describes planar movements of the upper arm and forearm, which are moved by six lumped muscles with nonlinear dynamics. The motor control system is represented by a neural network which combines feedforward and feedback control. It is optimized for the control of movements or for posture control in the presence of external forces. The achieved impedance characteristics depend on the conditions during the learning process. In particular, the impedance is adapted in a suitable way to the frequency content and direction of external forces acting on the hand during an isometric task. The impedance characteristics of a model, which is optimized for movement control, are similar to experimental data in the literature. The achieved stiffness is, to a large extent, reflexively determined whereas the approximated viscosity is primarily due to intrinsic attributes. It is argued that usually applied Hill-type muscle models do not properly represent intrinsic muscle stiffness. Received: 14 October 1997 / Accepted in revised form: 18 May 1999     

Single-molecule measurement of the stiffness of the rigor myosin head

The force-extension curve of single myosin subfragment-1 molecules, interacting in the rigor state with an actin filament, has been investigated at low [ATP] by applying a slow triangle-wave movement to the optical traps holding a bead-actin-bead dumbbell. In combination with a measurement of the overall stiffness of the dumbbell, this allowed characterization of the three extensible elements, the actin-bead links and the myosin. Simultaneously, another method, based on an analysis of bead position covariance, gave satisfactory agreement. The mean covariance-based estimate for the myosin stiffness was 1.79 pN/nm (SD = 0.7 pN/nm; SE = 0.06 pN/nm (n = 166 myosin molecules)), consistent with a recent report (1.7 pN/nm) from rabbit muscle fibers. In the triangle-wave protocol, the motion of the trapped beads during interactions was linear within experimental error over the physiological range of force applied to myosin (±10 pN), consistent with a Hookean model; any nonlinear terms could not be characterized. Bound states subjected to forces that resisted the working stroke (i.e., positive forces) detached at a significantly lower force than when subjected to negative forces, which is indicative of a strain-dependent dissociation rate.     

重复经颅磁刺激联合等速肌力训练对不完全性脊髓损伤患者神经电生理指标、下肢肌力和脊髓功能独立性的影响

摘要 目的：探讨重复经颅磁刺激（TMS）联合等速肌力训练对不完全性脊髓损伤（SCI）患者神经电生理指标、下肢肌力和脊髓功能独立性的影响。方法：选取2018年3月~2019年12月期间我院收治87例不完全性SCI患者,根据入院奇偶顺序分为观察组（n=44）和对照组（n=43）,两组均给予常规康复训练,对照组在此基础上联合等速肌力训练,观察组在对照组基础上联合TMS,对比两组神经电生理指标[静息运动阈值（RMT）和运动诱发电位（MEP）]、下肢肌力指标[屈、伸肌群的峰力矩（PT）、力矩加速能（TAE）以及胭绳肌与股四头肌肌力比率（H/Q）]、功能独立性评定(FIM)量表、疼痛简化McGill疼痛问卷(SF-MPQ）、Barthel指数评定量表（BI）评分。结果：治疗4周后,观察组RMT较治疗前降低,且低于对照组（P<0.05）;MEP较治疗前升高,且高于对照组（P<0.05）。治疗4周后,两组屈肌群PT、屈肌群TAE、伸肌群PT、伸肌群TAE、H/Q、FIM、BI评分均较治疗前升高,且观察组高于对照组（P<0.05）;两组SF-MPQ评分均较治疗前下降（P<0.05）,且观察组低于对照组（P<0.05）。结论：TMS联合等速肌力训练治疗不完全性SCI患者可刺激患者神经功能恢复,提高患者脊髓功能独立性,改善下肢肌力,减轻患者的神经性疼痛。     

Evidence of changes in load sharing during isometric elbow flexion with ramped torque

This study aimed to: (1) test the repeatability of Supersonic Shear Imaging measures of muscle shear elastic modulus of four elbow flexor muscles during isometric elbow flexion with ramped torque; (2) determine the relationship between muscle shear elastic modulus and elbow torque for the four elbow flexor muscles, and (3) investigate changes in load sharing between synergist elbow flexor muscles with increases in elbow flexor torque. Ten subjects performed ten isometric elbow flexions consisting of linear torque ramps of 30-s from 0 to 40% of maximal voluntary contraction. The shear elastic modulus of each elbow flexor muscle (biceps brachii long head [BB(LH)], biceps brachii short head [BB(SH)], brachialis [BA], and brachoradialis [BR]) and of triceps brachii long head [TB] was measured twice with individual muscles recorded in separate trials in random order. A good repeatability of the shape of the changes in shear elastic modulus as a function of torque was found for each elbow flexor muscle (r-values: 0.85 to 0.94). Relationships between the shear elastic modulus and torque were best explained by a second order polynomial, except BA where a higher polynomial was required. Statistical analysis showed that BB(SH) and BB(LH) had an initial slow change at low torques followed by an increasing rate of increase in modulus with higher torques. In contrast, the BA shear elastic modulus increased rapidly at low forces, but plateaued at higher forces. These results suggest that changes in load sharing between synergist elbow flexors could partly explain the non-linear EMG-torque relationship classically reported for BB during isometric efforts.