Stiffness-distortion sarcomere model for muscle simulation

A relatively simple method is presented for incorporating cross-bridge mechanisms into a muscle model. The method is based on representing force in a half sarcomere as the product of the stiffness of all parallel cross bridges and their average distortion. Differential equations for sarcomeric stiffness are derived from a three-state kinetic scheme for the cross-bridge cycle. Differential equations for average distortion are derived from a distortional balance that accounts for distortion entering and leaving due to cross-bridge cycling and for distortion imposed by shearing motion between thick and thin filaments. The distortion equations are unique and enable sarcomere mechanodynamics to be described by only a few ordinary differential equations. Model predictions of small-amplitude step and sinusoidal responses agreed well with previously described experimental results and allowed unique interpretations to be made of various response components. Similarly good results were obtained for model reproductions of force-velocity and large-amplitude step and ramp responses. The model allowed reasonable predictions of contractile behavior by taking into account what is understood to be basic muscle contractile mechanisms.