Voluntary EMG-to-force estimation with a multi-scale physiological muscle model
© Hayashibe and Guiraud; licensee BioMed Central Ltd. 2013
Received: 8 January 2013
Accepted: 21 August 2013
Published: 4 September 2013
EMG-to-force estimation based on muscle models, for voluntary contraction has many applications in human motion analysis. The so-called Hill model is recognized as a standard model for this practical use. However, it is a phenomenological model whereby muscle activation, force-length and force-velocity properties are considered independently. Perreault reported Hill modeling errors were large for different firing frequencies, level of activation and speed of contraction. It may be due to the lack of coupling between activation and force-velocity properties. In this paper, we discuss EMG-force estimation with a multi-scale physiology based model, which has a link to underlying crossbridge dynamics. Differently from the Hill model, the proposed method provides dual dynamics of recruitment and calcium activation.
The ankle torque was measured for the plantar flexion along with EMG measurements of the medial gastrocnemius (GAS) and soleus (SOL). In addition to Hill representation of the passive elements, three models of the contractile parts have been compared. Using common EMG signals during isometric contraction in four able-bodied subjects, torque was estimated by the linear Hill model, the nonlinear Hill model and the multi-scale physiological model that refers to Huxley theory. The comparison was made in normalized scale versus the case in maximum voluntary contraction.
The estimation results obtained with the multi-scale model showed the best performances both in fast-short and slow-long term contraction in randomized tests for all the four subjects. The RMS errors were improved with the nonlinear Hill model compared to linear Hill, however it showed limitations to account for the different speed of contractions. Average error was 16.9% with the linear Hill model, 9.3% with the modified Hill model. In contrast, the error in the multi-scale model was 6.1% while maintaining a uniform estimation performance in both fast and slow contractions schemes.
We introduced a novel approach that allows EMG-force estimation based on a multi-scale physiology model integrating Hill approach for the passive elements and microscopic cross-bridge representations for the contractile element. The experimental evaluation highlights estimation improvements especially a larger range of contraction conditions with integration of the neural activation frequency property and force-velocity relationship through cross-bridge dynamics consideration.
Any human movement is produced by muscular and skeletal systems controlled by the nervous system. Mechanism of the human body dynamics has been revealed by many biomechanical researches. The general neuromusculoskeletal model of the whole body and its dynamics computation package have been developed and reported. However, the detailed neuromuscular activation system should still be closely analyzed and modeled from microscopic to macroscopic scales. Macroscopic modeling, such as the Hill-type muscle model, is a phenomenological model. Conversely, microscopic modeling is based on muscle physiology, more specifically the dynamics of actin-myosin, allowing it to display a richer response.
Neuromuscular modeling is important for neuroscience to understand how limb movements are smoothly and effectively controlled[2, 3]. It is also valuable for clinical applications to be used for quantitative analysis on abnormal muscle activation patterns such as spasticity induced by stroke or cerebral palsy. Moreover, neuroprosthetic control such as functional electrical stimulation (FES) for paralyzed muscles[5, 6] can be improved through neuromuscular modeling to optimize the muscle activation control.
Surface electromyography (EMG) is widely used to assess the muscle status in a wide range of clinical fields. The signal generated during a contraction is the summation of the signals of all different recruited motor units (MU) at a given time. Each MU receives information from the central nervous system in the form of action potentials, which are transmitted from the alpha motor neuron to the muscle fibers via the neuromuscular junction. During voluntary contraction, the EMG signal can be seen as the interference pattern of all active MUs, with each MU firing at its own mean frequency with its own timing. The contributions of individual MU cannot be easily recovered from the EMG signal. However, EMG signals still contain precious information about muscle activation and could thus represents the entry of a muscle model. Indeed, the important advantage of EMG usage is that it can account for a subject’s individual activation pattern to estimate muscle force. Moreover, as the EMG signals can be considered for each individual groups of antagonist muscles with their own muscle models, co-contraction can be captured and deeply investigated through the computation of individual force and stiffness. On the contrary, numerical optimization techniques typically do not account for muscular co-contraction.
EMG-based models were already used in many studies to estimate torques around joints or in musculoskeletal models[9, 10]. Most of the muscle models used in such studies are based on phenomenological models derived from Hill’s study and summarized by Zajac. The Hill-type model is attractive in biomechanical simulations because of the computational simplicity and easy to calibrate using experimentally measurable variables. It tends to produce accurate estimation results, however in limited ranges of function.
However, recent work performed the validation of Hill model under functionally wider relevant conditions. The authors concluded that the modeling errors were substantial for different firing frequencies and greatest at low motor unit firing rates that are most relevant in normal movement conditions both in voluntary and FES muscle contractions. They pointed out that this could be due to the Hill model assumption whereby muscle activation, force-length and force-velocity properties are considered independently. In their discussion, it was suggested that more physiological coupling between activation and force-velocity properties can be demonstrated in crossbridge models incorporating dependence between physiologically based activation and cross-bridge attachment. The important finding is that Hill modeling errors were large for different firing frequencies and greatest at low motor unit firing rates. In, classical linear Hill model was used. Especially for solving the problem of the mismatch at low muscle activation level, a nonlinear conversion was proposed and used by other researchers. This can solve partially the problem, however, there is a still the problem to correspond to the estimations in different contraction speeds. Besides, in Hill-type modeling, the cut off frequency should be carefully chosen depending on the type of task because the envelope of the estimated force is almost defined by the low-pass filtering of the activation. Thus, we aim at proposing multiscale muscle model, which can have more physiological coupling between activation and force-velocity properties by incorporating microscopic cross-bride and dual dynamics pathway for the recruitment / activation level and the calcium triggering the contraction cycle. However, the whole passive structure is based on Hill model as only the contractile element differs.
The first microscopic muscle model was proposed by Huxley. He detailed the interaction of the cross-bridges in a sarcomere in order to explain the force generation. The distinctions between microscopic and macroscopic levels are not absolute; thus a sarcomere model can be used to represent a whole muscle, which is assumed to be a homogeneous assembly of identical sarcomeres. Conversely, the Hill model, which was originally developed for whole muscles, has been used to represent the dynamics of individual sarcomeres within a fiber. The distribution-moment approach of Zahalak is a model for sarcomeres or whole muscle, which is extracted via a formal mathematical Gaussian approximation from Huxley cross-bridge model. This model fills a gap between microscopic and macroscopic levels. Based on Huxley and Hill-type models, Bestel-Sorine proposed an explanation of how the beating of cardiac muscle may take place, through chemical control input connected to the calcium dynamics of cells in muscles. This chemical input triggers the contractile element model.
Complete Huxley type physiological model have a greater repertoire of response, but are computationally complex. Bestel-Sorine representation provides an efficient trade-off between complexity and power of the model expression. It consists in an explicit computation of the first two moments reprenting stiffness and force of the contractile element, triggered by the chemical input as regard the contraction start and stop. Starting with this concept, we first adapted it to the striated muscle model to represent muscle responses under FES. The multi-scale muscle model was applied for voluntary contraction in a preliminary work. In this study, we aim at developing EMG-based muscle model to study voluntary muscle contraction with this multiscale modeling approach compared to Hill approaches, both linear and nonlinear. Indeed, in Hill approach, the dependency on cut-off frequency choice prevents from good performances in a wide range of contraction types. The multiscale approach should improve the performances in particular regarding fast and slow contraction patterns. In addition, we consider that it is of interest to apply the same type of EMG-to-force estimation with a physiologically detailed model and not only with a phenomenological Hill model because the internal biophysical dynamics could shed new light on neuromuscular activation. Using common data sets of isometric muscle contractions, the force estimation results are compared between Hill model approaches and newly proposed approach.
Hill-type muscle model
The transformation from EMG to muscle activation is a dominant process because the estimated muscle force is assumed to be proportional to the muscle activation.
high-pass filtering of raw EMG using zero-lag 4th order Butterworth filter (30 Hz) to remove motion artifacts
low-pass filtering with a 2 Hz cut-off frequency using zero-lag 2nd order Butterworth filter
normalization with the peak of maximum voluntary contraction (MVC)
There is an alternative method instead of the low-pass filtering of the third step. It is referred to as the mean absolute value (MAV) method. A moving average filter is applied within a fixed time window, which is equivalent to low-pass filtering. The recommendations of the SENIAM project indicate a time window of 0.25 to 0.5 seconds. It corresponds to the choice of cut-off frequency of the low-pass filter of 2 to 4 Hz. This critical parameter choice depends on the task dynamics. Indeed, this choice of cut off frequency is a trade off between how fast the fluctuations in amplitude can be and how reliable the estimate of the amplitude is. A low cut off frequency gives a more reliable estimate of the amplitude, but cannot capture fast changes, while higher cut off frequency results in noisier output but captures fast contractions. The cut off frequency mentioned in the SENIAM recommendations is 2 Hz for slow motions, and 6 Hz for fast motions. In Hill-type model, this cut off frequency should be carefully set depending on the type of task because the envelope of the estimated force is almost determined by this low-pass filtering step. The normalized EMG is processed with the following activation dynamics, which mainly captures the delay from muscle activation to mechanical output.
where d is the electromechanical delay and γ, β1 and β2 are the coefficients that define the second-order dynamics. To obtain a positive stable solution, a set of constraints are employed, i.e. β1=C1+C2,β2=C1C2 where |C1|<1,|C2|<1. In addition, the gain of this filter should be maintained by ensuring γ−β1−β2=1 as described in.
Nonlinearization of muscle activation for nonlinear Hill model
Many researchers assume that the above p(t) is a reasonable approximation of muscle activation. However, a nonlinear relationship has been reported between individual muscle EMG and the joint moment for some muscles, especially at lower forces. In studies on single motor units, multiple APs cause succesive twitch responses. If the time between APs decreases, twitches start to merge and the muscle force increases steadily. However, at high frequency we reach maximal contraction, where force development saturates even if the frequency increases. This induces a nonlinear relationship between the frequency and force for single motor units.
where A is a constant parameter for the nonlinear shape factor allowed to vary between -3 and 0, with A=−3 being highly exponential and A=0 being linear.
Hill-type contraction dynamics
where ε c is the strain of the contractile element, f l (ε c ) and are the normalized force-length and force-velocity relationships, respectively. F m is the maximum isometric muscle force.
where Lc 0 is the length of contractile element at rest position and V sh is a constant parameter. The b and V sh parameters were obtained from[23, 24]. At each time step, the fiber velocity should be solved and the muscle fiber length is computed by forward integration using Runge-Kutta algorithm. Since the value of ε c changed, the calculation should continue iteratively until the end of the input time series of a(t).
Muscle model parameters
Lc 0(G A S)
Lt 0(G A S)
F m (GAS)
Lc 0(S O L)
Lt 0(S O L)
F m (SOL)
F m ×20
Physiological muscle model
We breifly describe the main points of the multi-scale physiological muscle model. It integrates the characteristics both of macroscopic and microscopic muscle dynamics. The details of this modeling are given in. This model was originally designed to represent the muscle response under FES and is applied to voluntary muscular force estimation based on EMG in this paper.
Generation of chemical input
Thus, the chemical threshold consists basically in extracting the contraction time duration of the muscle and to indicating if the muscle is in contraction or not. The level of activation itself is managed by recruitment rate which is classically used in EMG usage for Hill model. Thus, especially for the low activation state, even if the chemical input is inaccurately determined, its effect is low since the recruitment activation p(t) is very low during this period. And for the higher activation state, the contraction period is accurately captured by this method. As long as the threshold is taken little above the baseline, the threshold is less sensitive as long as it gives the picture of contraction/relaxation state. The parameters for chemical input generation was obtained from[18, 28].
All the sarcomeres are assumed to be identical, and the deformation of both sarcomere and muscle is proportional. If S is the sarcomere length, we can write (S−S0)/S0=(L c −Lc 0)/Lc 0=ε c .
Huxley proposed that a cross-bridge between actin filaments and myosin heads could exist in two biochemical states, i.e. attached and detached states as in Figure2(B). Filaments sliding is the result of interactions between the myosin cross-bridges and the thin actin filaments. The cross-bridges reversibly bind to actin and produce a mechanical impulse, which in turn results in force transmission along the filaments. M and A in the figure represents the myosin head and actin binding site respectively.
where h is the maximum elongation of the myosin spring, x is the distance and y the normalized distance between the actin binding site and the myosin head:. n(y,t) is a distribution function representing the fraction of attached cross bridges relative to the normalized relative position y. To ensure a one way displacement, this attachment is considered possible only when y is between 0 and 1. represents the velocity of the actin filament relative to the myosin filament. f and g denote the rate functions of attachment and detachment, respectively.
Contrary to Zahalak approximation, we do not make any assumption on the distribution, rather the chosen f and g functions allows for a straightforward computation of the stiffness and the force generated by the contractile element. It provides a computational effective model that however relies on Huxley theory.
Myofiber and muscle scale
where k m =S0N k0f l (ε c )/Lc 0,F m =N k0h f l (ε c )/2. k0 (Nm −1) is the maximum stiffness obtained when all the available bridges are attached.
Recruitment rate is thus mixed with the dynamics of the contraction relaxation cycle providing a wider range of mechanical responses than Hill approaches. Compared to Zahalak approximation, we provide a two input model with intricate dynamics due to both internal muscle dynamics and recruitment rate dynamics in a straightforward way.
Computation of the whole dynamics
For the macroscopic representation, the same configuration with Hill models as in Figure2(A) is used including the muscle tendon parameters. The contractile element is replaced with the above nonlinear differential equations.
is the sign of. From the condition: k t Lc 0/c o s ϕ+k c Lc 0−F c >0, can be obtained from the sign of these terms:. Then we can compute and with (14). The internal state vector of this system should be set as.
Experimental measurement and data processing
Four healthy volunteers (2 males and 2 females, age =30.8±1.3 years, mass =66.5±16.0 kg, and stature =1.7±0.13 m) participated in the study after signing an informed consent form. The study was approved by the Agence française de sécurité sanitaire des produits de santé (Afssaps) committee for persons’ protection managed by CHRU Montpellier (2011-A01033-38), dedicated to the study of the lower limb biomechanics within ANR SoHuSim project.
In this preliminary test, we focused on the isometric condition and tried to compare the output force on a normalized scale against the output during MVC. The common parameters between three approaches were set at the same value as summarized in Table1. In order to maintain similar conditions, the common macroscopic model as in Figure2(A) was used for all the three approaches. Only the dynamics representation of the contractile element was replaced. A was the only parameter identified for each subject. A was obtained by minimizing the root mean square difference between the measured and the model estimation at MVC. Other muscle parameters except A are taken from literatures as explained in the model section.
The estimated muscle force is multiplied by the moment arm in order to obtain the torque. The moment arm was estimated from the Hawkins data from the joint angle in the measured condition. The moment arms for MG and SOL are 0.0515 [m] and 0.0464 [m], respectively. The contribution ratio was calculated using the values reported by Delp. This was obtained as the maximum force by the moment arm considering the pennation angle. The resulting ratio is MG 0.41 vs SOL 0.59. The SUM in the estimation result is plotted using the ratio as the sum of the two muscles.
The generated input command u(t) and activation rate α were given to the contractile element of the physiological model and the active stiffness k c and the muscle force F c were computed using Eq. 14. The same value for F m for the three approaches were used. k m is proportional to F m from the condition in Eq. 14. It was multipied by 20 times, as shown in Table1, based on the knowledge that the muscle active stiffness reaches the tendon stiffness at approximately maximum contraction. However, this muscle stiffness parameter does not influence much the force output in isometric condition.
The predicted joint torque based on EMG signals was compared with the directly measured torque of the ankle joint for the plantar flexion. The comparison was made against the case in MVC in normalized scale. Thus, we avoided the situation effected by the subject-specific difference regarding muscle strength, which is significant for torque estimation in absolute scale.
RMS errors between the measured and estimated results in 30% and 70% of MVC
RMS errors between the measured and estimated results in random contraction
Peak error (fast)
Peak error (fast)
Peak error (fast)
Hill-type muscle is a phenomenological model based on experimental facts with no link to the microscopic physiology. First, even if the estimation performances are equivalent, it is meaningful to understand and capture the muscle dynamics with a more detailed representation. Indeed, the same kind of force estimation could be obtained with the newly proposed physiological model using voluntary EMG signals.
Second, previous work by Perreault reported that Hill modeling errors were large for different firing frequencies and greatest at low MU firing rates. The classical linear Hill model was used. They suggested that more physiological coupling between activation and force-velocity properties may help to deal with these issues. Thus, the aim of introducing new approach based on a physiologically detailed model was to develop EMG-force estimation with natural velocity dependency incorporating microscopic cross-bridge dynamics on top of recruitment dynamics.
In general, the estimation accuracy of the proposed method is better than that of the linear and nonlinear Hill-type approach as in the obtained result regarding RMS errors. In addition, uniform accuracy was obtained for different types of contractions in the proposed method. In particular, the proposed model is able to properly deal with both slow and fast contraction speeds. When considering the Hill model results, Figure6(A,B) shows largest errors for the estimation of fast-short contractions. The reason of this error could be explained as follows: the signal measured in EMG is the summation of the APs of all different MUs. Even in the fast-short contraction, the amplitude of EMG itself does not reflect differences except duration with the same level as for slow-long contraction. However, the resulting fast-short contraction force is actually much less than that of slow-long contraction. This means that there is sometimes hysteresis regarding the neural command. In Hill approach, muscle activation a(t) is proportional to the resulting force. It therefore does not include the effect of the time hysteresis in contraction, thus Hill model can not accurately represent both short and long-term contractions simultaneously with only one choice of cut-off frequency in EMG low-pass filtering. In the proposed approach, the derivative of the contraction force is directly given by the neural command and it brings time hysteresis in force generation. It is interesting to note that this kind of uniform accuracy for different types of contractions was appeared along with the introduction of cross-bridge dynamics.
In addition, it is known that there is a nonlinear relationship between frequency during contraction and force for single MUs. In modified Hill model, this frequency dependency is offset only by the nonlinear conversion. This nonlinearization was recently proposed to modify the classical Hill model. This process was originally not introduced in the so-called Hill model. In fact, this modification brings much better estimation, especially at lower forces. However, even with the modification, it is not a time function so it still can not correspond to the different muscle contraction speeds. In Hill-type approach, neural activation p(t) is dominantly imposed by low-pass filtering with a 2–10 Hz cut-off frequency. The choice of cut-off frequency is very sensitive to the dynamics of p(t). p(t) is more or less proportional to the resulting force. Therefore, it is advised that the choice of cut off frequency should be carefully selected depending on the type of task in Hill-type modeling. In this study, a 2 Hz cut-off frequency was used according to the SENIAM recommendation. Sometimes by changing the cut-off frequency for different types of contraction, better results can be obtained with Hill-type model. However, in any case, the expected motion by subject is unknown in advance for practical use of EMG-force estimation, so it is difficult to have an appropriate cut-off frequency in advance. Moreover it assumes that the subject stays in the same kind of motion dynamics, that is not always the case. In the proposed model, the nonlinear activation property coming from force-velocity dependency is internally integrated in the model dynamics.
Furthermore, the physiological model satisfies the well established properties observed for the muscle’s behavior: i) the force-length relation is included in the definition of Eq. 13, ii) the force-velocity relationship can be expressed in isotonic and tetanic conditions. Before isotonic contraction could occur, the muscle contracts in isometric conditions until the force generated by the muscle balances the imposed one. Isotonic contraction then becomes possible. Assuming in Eq. 14, the following equation may be formulated:.
This correspondence can be found in Eq. 5. It can be verified that our physiological muscle model integrates the force-velocity relation naturally when considering the actin-myosin cross bridge. Finally, this improvement in estimation accuracy is in line with the suggestion put forward in a previous work on introducing crossbridge model.
Future study will be focused on simulations with multiple sets of the proposed model corresponding to multiple MUs in one muscle considering size principal, fast/slow fibers. However, for the application of voluntary EMG-force estimation, we should keep the simplicity as Hill-type model. Future work will also focus on increasing the number of experimental tests and we should develop an algorithm to statistically determine if a muscle is turned on/off regarding the chemical thresholding. In addition, the further interpretation of neuromuscular system both in voluntary and FES activations should be pursued. We have already performed the study on the identification and validation of the physiological model in in-vivo rabbit experiment and in paraplegic subjects including non-isometric situation under FES. In order to be applied to a broader range of clinical situations, further investigations in isokinetic and isotonic cases should be carried out.
In this paper we presented a method that allows estimation of muscle force from EMG signals with a multi-scale physiology based model with a link to underlying microscopic filament dynamics. The experimental results highlight the feasibility of the torque estimation and its comparison with Hill-type models using the same EMG signals. This is the first report on EMG-force estimation based on a multi-scale physiology model integrating Hill-type and microscopic cross bridge representations. The proposed method features:
a novel physiologically detailed model for EMG-force estimation in the place of a phenomenological Hill-type muscle model,
the estimation improvement especially for different types of contraction incorporating dual dynamics of recruitment and microscopic cross-bridge formation toward coupled activation-velocity relationship.
This work was supported by French National ANR SoHuSim project (http://sohusim.gforge.inria.fr/).
- Delp SL, Anderson FC, al E: OpenSim: Open-source software to create and analyze dynamic simulations of movement. IEEE Trans Biomed Eng 2007, 54(11):1940–1950.View ArticleGoogle Scholar
- Flash T, Hogan N: The coordination of arm movements: an experimentally confirmed mathematical model. J Neurosci 1985, 5(7):1688–1703.Google Scholar
- Kawato M, Maeda Y, Uno Y, Suzuki R: Trajectory formation of arm movement by cascade neural network model based on minimum torque-change criterion. Biol Cybern 1990, 62(4):275–288. 10.1007/BF00201442View ArticleGoogle Scholar
- Smaga S: Tremor. Am Fam Physician 2003, 68: 1545–1552.Google Scholar
- Kobetic R, Triolo RJ, Uhlir JP, Bieri C, Wibowo M, Polando G, Marsolais EB, Davis JA: Implanted functional electrical stimulation system for mobility in paraplegia:a follow-up case report. IEEE Trans Rehabil Eng 1999, 7(4):390–398. 10.1109/86.808942View ArticleGoogle Scholar
- Guiraud D, Stieglitz T, Koch KP, Divoux JL, Rabischong P: An implantable neuroprosthesis for standing and walking in paraplegia: 5-year patient follow-up. J Neural Eng 2006, 3: 268–275. 10.1088/1741-2560/3/4/003View ArticleGoogle Scholar
- Stegeman DF, Blok JH, Hermens HJ, Roeleveld K: Surface EMG models: properties and applications. J Electromyogr Kinesiol 2000, 10: 313–326. 10.1016/S1050-6411(00)00023-7View ArticleGoogle Scholar
- Buchanan TS, Lloyd DG, Manal K, Bessier TF: Neuromusculoskeletal modeling: estimation of muscle forces and joint moments and movements from measurements of neural command. J Appl Biomech 2004, 20: 367–395.Google Scholar
- Lloyd DG, Bessier TF: An emg-driven musculoskeletal model to estimate muscle forces and knee joint moment in vivo. J Biomech 2003, 36: 765–776. 10.1016/S0021-9290(03)00010-1View ArticleGoogle Scholar
- Koo TKK, Mak AFT: Feasibility of using emg driven neuromusculoskeletal model for prediction of dynamic movement of the elbow. J Electromyogr Kinesiol 2005, 15: 12–26. 10.1016/j.jelekin.2004.06.007View ArticleGoogle Scholar
- Hill AV: The heat of shortening and the dynamic constants in muscle. Proc R S Ser B; London 1938, 126: 136–195. 10.1098/rspb.1938.0050View ArticleGoogle Scholar
- Zajac FE: Muscle and tendon: properties, models, scaling and application to biomechanics and motor control. CRC Critic Rev Biomed Eng 1989, 17: 359–411.Google Scholar
- Perreault E, Heckman C, Sandercock T: Hill muscle model errors during movement are greatest within the physiologically relevant range of motor unit firing rates. J Biomech 2003, 36: 211–218. 10.1016/S0021-9290(02)00332-9View ArticleGoogle Scholar
- Manal K, Buchanan TS: A one-parameter neural activation to muscle activation model. J Biomech 2003, 36: 1197–1202. 10.1016/S0021-9290(03)00152-0View ArticleGoogle Scholar
- Huxley AF: Muscle structure and theories of contraction. Prog Biophys Biophys Chem 1957, 7: 255–318.Google Scholar
- Leeuwen JL: Optimum power output and structural design of sarcomeres. J Theor Biol 1991, 19: 229–256.View ArticleGoogle Scholar
- Zahalak GI: A distribution-moment approximation for kinetic theories of muscular contraction. Math Biosci 1981, 55: 89–114. 10.1016/0025-5564(81)90014-6View ArticleGoogle Scholar
- Bestel J, Clement F, Sorine M: A biomechanical model of muscle contraction. Proc Int Conf Med Image Comput Comput-Assist Intervent MICCAI 2001, Lect Notes Comput Sci 2001, 2208: 1159–1161.Google Scholar
- Makssoud HE, Guiraud D, Poignet P, Hayashibe M, Wieber PB, Yoshida K, Azevedo C: Multiscale modeling of skeletal muscle properties-experimental validations in isometric conditions. Biol Cybern 2011, 105(2):121–138. 10.1007/s00422-011-0445-7View ArticleGoogle Scholar
- Hayashibe M, Guiraud D, Poignet P: EMG-to-force estimation with full-scale physiology based muscle model. Proc IEEE/Int Conf Intell Robots Syst 2009, 1621–1626.Google Scholar
- Hermens HJ, Freriks B, Merletti R, Stegeman D, Blok J, Rau G, Disselhorst-Klug C, Hagg G: European recommendations for surface electromyography, results of the seniam project. Enschede, The Netherlands: Roessingh Research and Development; 1999.Google Scholar
- Lichtwark GA, Wilson AM: Optimal muscle fascicle length and tendon stiffness for maximising gastrocnemius efficiency during human walking and running. J Theor Biol 2008, 252(4):662–673. http://www.ncbi.nlm.nih.gov/pubmed/18374362 10.1016/j.jtbi.2008.01.018MathSciNetView ArticleGoogle Scholar
- Stroeve S: Learning combined feedback and feedforward control of a musculoskeletal system. Biol Cybern 1996, 75: 73–83. 10.1007/BF00238741View ArticleGoogle Scholar
- Stroeve S: Impedance characteristics of a neuro-musculoskeletal model of the human arm I: posture control. Biol Cybern 1999, 81: 475–494. 10.1007/s004220050577View ArticleGoogle Scholar
- Delp SL: Surgery simulation: a computer graphics system to analyze and design musculoskeletal reconstructions of the lower limb. PhD thesis, Dissertation, Stanford University, CA, USA; 1990Google Scholar
- Keener J, Sneyd J: Mathematical Physiology, II: Systems Physiology. Series 8, Vol. 2. New York: Springer; 2009.Google Scholar
- Hatze H: A general myocybernetic control model of skeletal muscle. Biol Cybern 1978, 28: 143–157. 10.1007/BF00337136View ArticleGoogle Scholar
- Bestel J: Modele Differentiel de la Contraction Musculaire controlee application au system cardio-vasculaire. PhD thesis, Dissertation, Universite de Paris IX, France; 2000Google Scholar
- Hawkins D, Hull M: A method for determining lower extremity muscle-tendon lengths during flexion/extension movements. J Biomech 1990, 23: 487–494. 10.1016/0021-9290(90)90304-LView ArticleGoogle Scholar
- Cook C, McDonagh M: Measurement of muscle and tendon stiffness in man. Eur J Appl Physiol Occup Physiol 1996, 72(4):380–382. 10.1007/BF00599700View ArticleGoogle Scholar
- Hayashibe M, Guiraud D, Poignet P: In vivo identification of skeletal muscle dynamics with nonlinear Kalman filter: comparison between EKF and SPKF. ISRN Rehabil 2013. 2013: Article ID 610709Google Scholar
- Benoussaad M, Poignet P, Hayashibe M, Azevedo C, Fattal C, Guiraud D: Experimental parameter identification of a multi-scale musculoskeletal model controlled by electrical stimulation: application to patients with spinal cord injury. Med Biol Eng Comput 2013, 51(6):617–631. 10.1007/s11517-013-1032-yView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.