Voluntary EMG-to-force estimation with a multi-scale physiological muscle model

BioMedical Engineering OnLine

Mitsuhiro Hayashibe¹  David Guiraud¹ 

[1] INRIA DEMAR Project and LIRMM, UMR5506 CNRS University of Montpellier, 161 Rue Ada, 34095 Montpellier, France

关键词: Multi-scale physiology;    Hill model;    Muscle force estimation;    EMG;    Muscle model;   

Others  :  797384
DOI  :  10.1186/1475-925X-12-86

 received in 2013-01-08, accepted in 2013-08-21,  发布年份 2013

【 摘 要 】

Background

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.

Methods

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.

Results

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.

Conclusions

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.

【 授权许可】

   
2013 Hayashibe and Guiraud; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

Figure 6.

【 参考文献 】

• [1]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.
• [2]Flash T, Hogan N: The coordination of arm movements: an experimentally confirmed mathematical model. J Neurosci 1985, 5(7):1688-1703.
• [3]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.
• [4]Smaga S: Tremor. Am Fam Physician 2003, 68:1545-1552.
• [5]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.
• [6]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.
• [7]Stegeman DF, Blok JH, Hermens HJ, Roeleveld K: Surface EMG models: properties and applications. J Electromyogr Kinesiol 2000, 10:313-326.
• [8]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.
• [9]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]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.
• [11]Hill AV: The heat of shortening and the dynamic constants in muscle. Proc R S Ser B; London 1938, 126:136-195.
• [12]Zajac FE: Muscle and tendon: properties, models, scaling and application to biomechanics and motor control. CRC Critic Rev Biomed Eng 1989, 17:359-411.
• [13]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.
• [14]Manal K, Buchanan TS: A one-parameter neural activation to muscle activation model. J Biomech 2003, 36:1197-1202.
• [15]Huxley AF: Muscle structure and theories of contraction. Prog Biophys Biophys Chem 1957, 7:255-318.
• [16]Leeuwen JL: Optimum power output and structural design of sarcomeres. J Theor Biol 1991, 19:229-256.
• [17]Zahalak GI: A distribution-moment approximation for kinetic theories of muscular contraction. Math Biosci 1981, 55:89-114.
• [18]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.
• [19]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.
• [20]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.
• [21]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.
• [22]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 webcite
• [23]Stroeve S: Learning combined feedback and feedforward control of a musculoskeletal system. Biol Cybern 1996, 75:73-83.
• [24]Stroeve S: Impedance characteristics of a neuro-musculoskeletal model of the human arm I: posture control. Biol Cybern 1999, 81:475-494.
• [25]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; 1990
• [26]Keener J, Sneyd J: Mathematical Physiology, II: Systems Physiology. Series 8, Vol. 2. New York: Springer; 2009.
• [27]Hatze H: A general myocybernetic control model of skeletal muscle. Biol Cybern 1978, 28:143-157.
• [28]Bestel J: Modele Differentiel de la Contraction Musculaire controlee application au system cardio-vasculaire. PhD thesis,Dissertation, Universite de Paris IX, France; 2000
• [29]Hawkins D, Hull M: A method for determining lower extremity muscle-tendon lengths during flexion/extension movements. J Biomech 1990, 23:487-494.
• [30]Cook C, McDonagh M: Measurement of muscle and tendon stiffness in man. Eur J Appl Physiol Occup Physiol 1996, 72(4):380-382.
• [31]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 610709
• [32]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.