ABSTRACT
Muscles convert metabolic energy into mechanical work. A computational model of muscle would ideally compute both effects efficiently for the entire range of muscle activation and kinematic conditions (force and length). We have extended the original Virtual Muscle algorithm (Cheng , 2000) to predict energy consumption for both slow- and fast-twitch muscle fiber types, partitioned according to the activation process (E(a)), cross-bridge cycling (E(xb)) and ATP/PCr recovery (E(recovery)). Because the terms of these functions correspond to identifiable physiological processes, their coefficients can be estimated directly from the types of experiments that are usually performed and extrapolated to dynamic conditions of natural motor behaviors. We also implemented a new approach to lumped modeling of the gradually recruited and frequency modulated motor units comprising each fiber type, which greatly reduced computational time. The emergent behavior of the model has significant implications for studies of optimal motor control and development of rehabilitation strategies because its trends were quite different from traditional estimates of energy (e.g., activation, force, stress, work, etc.). The model system was scaled to represent three different human experimental paradigms in which muscle heat was measured during voluntary exercise; predicted and observed energy rate agreed well both qualitatively and quantitatively.
