Research Papers

A Highly Efficient Semiphenomenological Model of a Half-Sarcomere for Real-Time Prediction of Mechanical Behavior

[+] Author and Article Information
Xing Chen

Robotics Institute,
State Key Laboratory of Mechanical System and Vibration,
Shanghai Jiao Tong University,
R. 301, Mechanical Building B,
800 Dong Chuan Road,
Shanghai 200240, China
e-mail: sing@sjtu.edu.cn

Yue Hong Yin

Robotics Institute,
State Key Laboratory of Mechanical System and Vibration,
Shanghai Jiao Tong University,
R. 914, Mechanical Building A,
800 Dong Chuan Road,
Shanghai 200240, China
e-mail: yhyin@sjtu.edu.cn

1Corresponding author.

Manuscript received April 22, 2014; final manuscript received September 2, 2014; accepted manuscript posted September 11, 2014; published online October 15, 2014. Assoc. Editor: Tammy L. Haut Donahue.

J Biomech Eng 136(12), 121001 (Oct 15, 2014) (9 pages) Paper No: BIO-14-1177; doi: 10.1115/1.4028536 History: Received April 22, 2014; Revised September 02, 2014; Accepted September 11, 2014

With existent biomechanical models of skeletal muscle, challenges still exist in implementing real-time predictions for contraction statuses that are particularly significant to biomechanical and biomedical engineering. Because of this difficulty, this paper proposed a decoupled scheme of the links involved in the working process of a sarcomere and established a semiphenomenological model integrating both linear and nonlinear frames of no higher than a second-order system. In order to facilitate engineering application and cybernetics, the proposed model contains a reduced number of parameters and no partial differential equation, making it highly concise and computationally efficient. Through the simulations of various contraction modes, including isometric, isotonic, successive stretch and release, and cyclic contractions, the correctness and efficiency of the model, are validated. Although this study targets half-sarcomeres, the proposed model can be easily extended to describe the larger-scale mechanical behavior of a muscle fiber or a whole muscle.

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


Hill, A., 1938, “The Heat of Shortening and the Dynamic Constants of Muscle,” Proc. R. Soc. Lond. Ser. B, 126(843), pp. 136–195. [CrossRef]
Huxley, A. F., 1957, “Muscle Structure and Theories of Contraction,” Prog. Biophys. Biophys. Chem., 7, pp. 255–318. [PubMed]
Julicher, F., Ajdari, A., and Prost, J., 1997, “Modeling Molecular Motors,” Rev. Mod. Phys., 69(4), pp. 1269–1281. [CrossRef]
Yin, Y., and Guo, Z., 2011, “Collective Mechanism of Molecular Motors and a Dynamic Mechanical Model for Sarcomere,” Sci. China Technol. Sci., 54(8), pp. 2130–2137. [CrossRef]
Fan, Y., and Yin, Y., 2013, “Active and Progressive Exoskeleton Rehabilitation Using Multi-Source Information Fusion From EMG and Force & Position-EPP,” IEEE Trans. Biomed. Eng., 60(12), pp. 3314–3321. [CrossRef]
Neptune, R. R., Burnfield, J. M., and Mulroy, S. J., 2007, “The Neuromuscular Demands of Toe Walking: A Forward Dynamics Simulation Analysis,” J. Biomech., 40(6), pp. 1293–1300. [CrossRef] [PubMed]
Yin, Y. H., Fan, Y. J., and Xu, L. D., 2012, “EMG and EPP-Integrated Human–Machine Interface Between the Paralyzed and Rehabilitation Exoskeleton,” IEEE Trans. Inf. Technol. Biomed., 16(4), pp. 542–549. [CrossRef] [PubMed]
Chen, X., and Yin, Y., 2013, “A Dynamical System-Markov Model for Active Postsynaptic Responses of Muscle Spindle Afferent Nerve,” Chin. Sci. Bull., 58(6), pp. 603–612. [CrossRef]
Destexhe, A., Mainen, Z. F., and Sejnowski, T. J., 1994, “Synthesis of Models for Excitable Membranes, Synaptic Transmission and Neuromodulation Using a Common Kinetic Formalism,” J. Comput. Neurosci., 1(3), pp. 195–230. [CrossRef] [PubMed]
Baylor, S. M., and Hollingworth, S., 2011, “Calcium Indicators and Calcium Signalling in Skeletal Muscle Fibres During Excitation–Contraction Coupling,” Prog. Biophys. Mol. Biol., 105(3), pp. 162–179. [CrossRef] [PubMed]
Stern, M. D., Pizarro, G., and Rios, E., 1997, “Local Control Model of Excitation–Contraction Coupling in Skeletal Muscle,” J. Gen. Physiol., 110(4), pp. 415–440. [CrossRef] [PubMed]
Endo, M., 2009, “Calcium-Induced Calcium Release in Skeletal Muscle,” Physiol. Rev., 89(4), pp. 1153–1176. [CrossRef] [PubMed]
Palmer, B. M., Suzuki, T., Wang, Y., Barnes, W. D., Miller, M. S., and Maughan, D. W., 2007, “Two-State Model of Acto-Myosin Attachment-Detachment Predicts C-Process of Sinusoidal Analysis,” Biophys. J., 93(3), pp. 760–769. [CrossRef] [PubMed]
Piazzesi, G., and Lombardi, V., 1995, “A Cross-Bridge Model That Is Able to Explain Mechanical and Energetic Properties of Shortening Muscle,” Biophys. J., 68(5), pp. 1966–1979. [CrossRef] [PubMed]
Propp, M. B., 1986, “A Model of Muscle Contraction Based Upon Component Studies,” Lect. Math. Life Sci., 16, pp. 61–119.
Zahalak, G. I., 2000, “The Two-State Cross-Bridge Model of Muscle is an Asymptotic Limit of Multi-State Models,” J. Theor. Biol., 204(1), pp. 67–82. [CrossRef] [PubMed]
Fan, Y., Guo, Z., and Yin, Y., 2011, “sEMG-Based Neuro-Fuzzy Controller for a Parallel Ankle Exoskeleton With Proprioception,” Int. J. Rob. Autom., 26(4), pp. 450–460. [CrossRef]
Tsianos, G. A., Rustin, C., and Loeb, G. E., 2012, “Mammalian Muscle Model for Predicting Force and Energetics During Physiological Behaviors,” IEEE Trans. Neural Syst. Rehab. Eng., 20(2), pp. 117–133. [CrossRef]
Iqbal, K., and Roy, A., 2004, “Stabilizing Pid Controllers for a Single-Link Biomechanical Model With Position, Velocity, and Force Feedback,” ASME J. Biomech. Eng., 126(6), pp. 838–843. [CrossRef]
Gollee, H., Murray-Smith, D. J., and Jarvis, J. C., 2001, “A Nonlinear Approach to Modeling of Electrically Stimulated Skeletal Muscle,” IEEE Trans. Biomed. Eng., 48(4), pp. 406–415. [CrossRef] [PubMed]
Bobet, J., Stein, R. B., and Oguztoreli, M. N., 1993, “A Linear Time-Varying Model of Force Generation in Skeletal Muscle,” IEEE Trans. Biomed. Eng., 40(10), pp. 1000–1006. [CrossRef] [PubMed]
Perreault, E. J., Heckman, C. J., and Sandercock, T. G., 2003, “Hill Muscle Model Errors During Movement Are Greatest Within the Physiologically Relevant Range of Motor Unit Firing Rates,” J. Biomech., 36(2), pp. 211–218. [CrossRef] [PubMed]
Dorgan, S. J., and O'malley, M. J., 1998, “A Mathematical Model for Skeletal Muscle Activated by N-Let Pulse Trains,” IEEE Trans. Rehab. Eng., 6(3), pp. 286–299. [CrossRef]
Kellermayer, M. S., Smith, S. B., Granzier, H. L., and Bustamante, C., 1997, “Folding–Unfolding Transitions in Single Titin Molecules Characterized With Laser Tweezers,” Science, 276(5315), pp. 1112–1116. [CrossRef] [PubMed]
Hatze, H., 1977, “A Myocybernetic Control Model of Skeletal Muscle,” Biol. Cybern., 25(2), pp. 103–119. [CrossRef] [PubMed]
Stoecker, U., Telley, I. A., Stüssi, E., and Denoth, J., 2009, “A Multisegmental Cross-Bridge Kinetics Model of the Myofibril,” J. Theor. Biol., 259(4), pp. 714–726. [CrossRef] [PubMed]
Huxley, A. F., 1974, “Muscular Contraction,” J. Physiol., 243(1), pp. 1–43. [PubMed]
Denoth, J., Stüssi, E., Csucs, G., and Danuser, G., 2002, “Single Muscle Fiber Contraction Is Dictated by Inter-Sarcomere Dynamics,” J. Theor. Biol., 216(1), pp. 101–122. [CrossRef] [PubMed]
Ramírez, A., Grasa, J., Alonso, A., Soteras, F., Osta, R., Muñoz, M., and Calvo, B., 2010, “Active Response of Skeletal Muscle: In Vivo Experimental Results and Model Formulation,” J. Theor. Biol., 267(4), pp. 546–553. [CrossRef] [PubMed]
Chen, X., Yin, Y., and Fan, Y., 2014, “EMG Oscillator Model-Based Energy Kernel Method for Characterizing Muscle Intrinsic Property Under Isometric Contraction,” Chin. Sci. Bull., 59(14), pp. 1556–1567. [CrossRef]
Ishii, Y., Nishiyama, M., and Yanagida, T., 2004, “Mechano-Chemical Coupling of Molecular Motors Revealed by Single Molecule Measurements,” Curr. Protein Pept. Sci., 5(2), pp. 81–87. [CrossRef] [PubMed]
Lymn, R. W., and Taylor, E. W., 1971, “Mechanism of Adenosine Triphosphate Hydrolysis by Actomyosin,” Biochemistry, 10(25), pp. 4617–4624. [CrossRef] [PubMed]
Huxley, A. F., and Simmons, R. M., 1971, “Proposed Mechanism of Force Generation in Striated Muscle,” Nature, 233(5321), pp. 533–538. [CrossRef] [PubMed]
De Tombe, P. P., and Ter Keurs, H., 1992, “An Internal Viscous Element Limits Unloaded Velocity of Sarcomere Shortening in Rat Myocardium,” J. Physiol., 454(1), pp. 619–642. [PubMed]
Roots, H., Offer, G., and Ranatunga, K., 2007, “Comparison of the Tension Responses to Ramp Shortening and Lengthening in Intact Mammalian Muscle Fibres: Crossbridge and Non-Crossbridge Contributions,” J. Muscle Res. Cell Motil., 28(2–3), pp. 123–139. [CrossRef] [PubMed]
Julian, F. J., and Moss, R. L., 1981, “Effects of Calcium and Ionic Strength on Shortening Velocity and Tension Development in Frog Skinned Muscle Fibres,” J. Physiol., 311(1), pp. 179–199. [PubMed]
Julian, F., 1971, “The Effect of Calcium on the Force–Velocity Relation of Briefly Glycerinated Frog Muscle Fibres,” J. Physiol., 218(1), pp. 117–145. [PubMed]
Meyer, G., Lieber, R., and Mcculloch, A., 2011, “A Nonlinear Model of Passive Muscle Viscosity,” ASME J. Biomech. Eng., 133(9), p. 091007. [CrossRef]
Sun, Y. B., Hilber, K., and Irving, M., 2001, “Effect of Active Shortening on the Rate of ATP Utilisation by Rabbit Psoas Muscle Fibres,” J. Physiol., 531(3), pp. 781–791. [CrossRef] [PubMed]
Josephson, R., 1993, “Contraction Dynamics and Power Output of Skeletal Muscle,” Annu. Rev. Physiol., 55(1), pp. 527–546. [CrossRef] [PubMed]
Josephson, R. K., 1985, “Mechanical Power Output From Striated Muscle During Cyclic Contraction,” J. Exp. Biol., 114(1), pp. 493–512.
Ko, C.-Y., Chang, Y., Kim, S.-B., Kim, S., Kim, G., Ryu, J., and Mun, M., 2014, “Linear- and Nonlinear-Electromyographic Analysis of Supracutaneous Vibration Stimuli of the Forearm Using Diverse Frequencies and Considering Skin Physiological Properties,” ASME J. Biomech. Eng., 136(1), p. 011008. [CrossRef]
Gordon, A. M., Huxley, A. F., and Julian, F. J., 1966, “The Variation in Isometric Tension With Sarcomere Length in Vertebrate Muscle Fibres,” J. Physiol., 184(1), pp. 170–192. [PubMed]


Grahic Jump Location
Fig. 1

The mechanical configuration of a half-sarcomere in the proposed model

Grahic Jump Location
Fig. 2

Activation kinetics represented by the isometric force of muscle. (a) The tension development and relax transients of a muscle under FES (functional electrical stimulation) of various frequencies. The dashed curves represent the smoothed trends. (b) The relationship between the stimulus frequency and the steady state isometric force. Here, the fitting parameter p is 5.77. The plots of all the experimental data are modified from Ref. [23].

Grahic Jump Location
Fig. 3

Schematic representation of the second-order system characterizing the activation kinetics

Grahic Jump Location
Fig. 4

(a) Schematic representation of the two-state cross-bridge model; (b) the bias of the distribution of P(x, t) and the mean deformation of a myosin motor under various sliding velocities

Grahic Jump Location
Fig. 5

Schematic representation of the fluctuation kinetics of the AE's force

Grahic Jump Location
Fig. 6

(a) The FaV relationships and FextV relationships under various activation degrees. The values of the parameters describing the FextV relationship at β = 1 are: a = 0.333, b = −0.333, A = 0.75, q = 15, and bc = 0.1; (b) The viscoelastic force of P1 under the sinusoidal length change as indicated by the inset. The parameters involved are: c1 = 0.04, c2 = 2.4, c3 = 2, and L1 = 2.6.

Grahic Jump Location
Fig. 7

Simulation of the force transients under the isometric stimulations of various intensities. All the stimuli are applied at the start and stopped at 0.4 s.

Grahic Jump Location
Fig. 8

(a) The length transients of the half-sarcomere under the isotonic contractions with various magnitudes of force steps; (b) the comparison of the FextV relationships between the simulated results and the experimental curve

Grahic Jump Location
Fig. 9

The force response of the half-sarcomere under the successive stretch and release at constant velocities during activation. The arrows indicate the time when the length ramps are applied or finished. The plot of experimental data is modified from Ref. [25].

Grahic Jump Location
Fig. 10

(a) The transients of f, β, Fext, and L during five cycles of stimulus with the phase of 10.4%; (b) the force–length orbits under the stimulus phase of 10.4%, 33%, −24%, and −4%, respectively. (c) The experimental force–length orbits under the stimulus phase of 18%, 33%, −24%, and −4%. The plots are modified from [41].

Grahic Jump Location
Fig. 11

The cybernetic architecture of the proposed model for the force-control mode, wherein the insets (a)–(d) represent the force–velocity relationship, the force–length relationship, the passive force of titin, and the response of drag force, respectively

Grahic Jump Location
Fig. 12

The form of G(L) that characterizes the isometric force–length relationship of a half-sarcomere



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