Research Papers

Critical Damping Conditions for Third Order Muscle Models: Implications for Force Control

[+] Author and Article Information
Davide Piovesan

Sensory Motor Performance Program (SMPP),
Rehabilitation Institute of Chicago,
Chicago, IL 60611
e-mail: d-piovesan@northwestern.edu

Alberto Pierobon

Ashton Graybiel Spatial Orientation Laboratory,
Brandeis University,
Waltham, MA 02454
e-mail: pierobon@brandeis.edu

Ferdinando A. Mussa Ivaldi

Sensory Motor Performance Program (SMPP),
Rehabilitation Institute of Chicago,
Chicago, IL 60611;
Department of Physiology,
Northwestern University,
Chicago, IL 60611
e-mail: sandro@nortwestern.edu

1Corresponding author.

Contributed by the Bioengineering Division of ASME for publication in the Journal of Biomechanical Engineering. Manuscript received January 29, 2013; final manuscript received July 10, 2013; accepted manuscript posted July 29, 2013; published online September 20, 2013. Assoc. Editor: Zong-Ming Li.

J Biomech Eng 135(10), 101010 (Sep 20, 2013) (8 pages) Paper No: BIO-13-1054; doi: 10.1115/1.4025110 History: Received January 29, 2013; Revised July 10, 2013; Accepted July 29, 2013

Experimental results presented in the literature suggest that humans use a position control strategy to indirectly control force rather than direct force control. Modeling the muscle-tendon system as a third-order linear model, we provide an explanation of why an indirect force control strategy is preferred. We analyzed a third-order muscle system and verified that it is required for a faithful representation of muscle-tendon mechanics, especially when investigating critical damping conditions. We provided numerical examples using biomechanical properties of muscles and tendons reported in the literature. We demonstrated that at maximum isotonic contraction, for muscle and tendon stiffness within physiologically compatible ranges, a third-order muscle-tendon system can be under-damped. Over-damping occurs for values of the damping coefficient included within a finite interval defined by two separate critical limits (such interval is a semi-infinite region in second-order models). An increase in damping beyond the larger critical value would lead the system to mechanical instability. We proved the existence of a theoretical threshold for the ratio between tendon and muscle stiffness above which critical damping can never be achieved; thus resulting in an oscillatory free response of the system, independently of the value of the damping. Under such condition, combined with high muscle activation, oscillation of the system can be compensated only by active control.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Fig. 1

Force control models of second-order mechanical systems. (a) Direct force control. (b) Indirect force control via position control. (c) Generalized impedance model of a second order mechanical system indicating the force fields that appear in the d'Alembert equation, Eq. (1).

Grahic Jump Location
Fig. 2

Third-order mechanical models of muscle-tendon systems. (a) Maxwell model (b) generalized impedance schematics of Maxwell model (c) Poynting-Thomson model (d) generalized impedance schematics of the PT model. (b) and (d) indicate the force fields that appear in the d'Alembert equation, Eq. (1). Each force field is a function of the mechanical paramenters of the elements that generate it.

Grahic Jump Location
Fig. 3

Discriminant of Poynting-Thomson model's characteristic polynomial shown as a function of muscle damping CθP for four values of tendon/muscle stiffness ratio κ (2.64, 4.85, 8, 10). The function is shown for both wrist (top) and elbow (bottom). For κ < 8 the discriminant is positive, independent of the value of CθP, which translate in a free oscillatory response of the system. If κ ≥ 8 a finite interval of damping values exists within which the system does not present an oscillatory free response.

Grahic Jump Location
Fig. 4

Roots of Poynting-Thomson model's characteristic polynomial shown as functions of the damping CθP for four values of tendon/muscle stiffness ratio κ (2.64, 4.85, 8, 10). The top row represents the real root; center and bottom rows represent the real and imaginary part of the two complex conjugate root, respectively.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In