Hill-type models are ubiquitous in the field of biomechanics, providing estimates of a muscle’s force as a function of its activation state and its assumed force-length and force-velocity properties. However, despite their routine use, the accuracy with which Hill-type models predict the forces generated by muscles during submaximal, dynamic tasks remains largely unknown. This study compared human gastrocnemii forces predicted by Hill-type models to the forces estimated from ultrasound-based measures of tendon length changes and stiffness during cycling, over a range of loads and cadences. We tested both a traditional model, with one contractile element, and a differential model, with two contractile elements that accounted for independent contributions of slow and fast muscle fibres. Both models were driven by subject-specific, ultrasound-based measures of fascicle lengths, velocities, and pennation angles and by activation patterns of slow and fast muscle fibres derived from surface electromyographic recordings. The models predicted on average, 54 % the time-varying gastrocnemii forces estimated from the ultrasound-based methods. However, differences between predicted and estimated forces were smaller under low speed-high activation conditions, with models able to predict nearly 80 % of the gastrocnemii force over a complete pedal cycle. Additionally, the predictions from the Hill-type muscle models tested here showed that a similar pattern of force production could be achieved for most conditions with and without accounting for the independent contributions of different muscle fibre types.
Flying animals must successfully contend with obstacles in their natural environments. Inspired by the robust manoeuvring abilities of flying animals, unmanned aerial systems are being developed and tested to improve flight control through cluttered environments. We previously examined steering strategies that pigeons adopt to fly through an array of vertical obstacles (VOs). Modelling VO flight guidance revealed that pigeons steer towards larger visual gaps when making fast steering decisions. In the present experiments, we recorded three-dimensional flight kinematics of pigeons as they flew through randomized arrays of horizontal obstacles (HOs). We found that pigeons still decelerated upon approach but flew faster through a denser array of HOs compared with the VO array previously tested. Pigeons exhibited limited steering and chose gaps between obstacles most aligned to their immediate flight direction, in contrast to VO navigation that favoured widest gap steering. In addition, pigeons navigated past the HOs with more variable and decreased wing stroke span and adjusted their wing stroke plane to reduce contact with the obstacles. Variability in wing extension, stroke plane and wing stroke path was greater during HO flight. Pigeons also exhibited pronounced head movements when negotiating HOs, which potentially serve a visual function. These head-bobbing-like movements were most pronounced in the horizontal (flight direction) and vertical directions, consistent with engaging motion vision mechanisms for obstacle detection. These results show that pigeons exhibit a keen kinesthetic sense of their body and wings in relation to obstacles. Together with aerodynamic flapping flight mechanics that favours vertical manoeuvring, pigeons are able to navigate HOs using simple rules, with remarkable success.
Numerous historical descriptions of the Lesser Egyptian jerboa, Jaculus jaculus, a small bipedal mammal with elongate hindlimbs, make special note of their extraordinary leaping ability. We observed jerboa locomotion in a laboratory setting and performed inverse dynamics analysis to understand how this small rodent generates such impressive leaps. We combined kinematic data from video, kinetic data from a force platform, and morphometric data from dissections to calculate the relative contributions of each hindlimb muscle and tendon to the total movement.
Jerboas leapt in excess of 10 times their hip height. At the maximum recorded leap height (not the maximum observed leap height), peak moments for metatarso-phalangeal, ankle, knee, and hip joints were 13.1, 58.4, 65.1, and 66.9 Nmm, respectively. Muscles acting at the ankle joint contributed the most work (mean 231.6 mJ / kg Body Mass) to produce the energy of vertical leaping, while muscles acting at the metatarso-phalangeal joint produced the most stress (peak 317.1 kPa). The plantaris, digital flexors, and gastrocnemius tendons encountered peak stresses of 25.6, 19.1, and 6.0 MPa, respectively, transmitting the forces of their corresponding muscles (peak force 3.3, 2.0, and 3.8 N, respectively). Notably, we found that the mean elastic energy recovered in the primary tendons of both hindlimbs comprised on average only 4.4% of the energy of the associated leap.
The limited use of tendon elastic energy storage in the jerboa parallels the morphologically similar heteromyid kangaroo rat, Dipodomys spectabilis. When compared to larger saltatory kangaroos and wallabies that sustain hopping over longer periods of time, these small saltatory rodents store and recover less elastic strain energy in their tendons. The large contribution of muscle work, rather than elastic strain energy, to the vertical leap suggests that the fitness benefit of rapid acceleration for predator avoidance dominated over the need to enhance locomotor economy in the evolutionary history of jerboas.
Electronic supplementary material
The online version of this article (doi:10.1186/s12983-017-0215-z) contains supplementary material, which is available to authorized users.
Mechanistically linking movement behaviors and ecology is key to understanding the adaptive evolution of locomotion. Predator evasion, a behavior that enhances fitness, may depend upon short bursts or complex patterns of locomotion. However, such movements are poorly characterized by existing biomechanical metrics. We present methods based on the entropy measure of randomness from Information Theory to quantitatively characterize the unpredictability of non-steady-state locomotion. We then apply the method by examining sympatric rodent species whose escape trajectories differ in dimensionality. Unlike the speed-regulated gait use of cursorial animals to enhance locomotor economy, bipedal jerboa (family Dipodidae) gait transitions likely enhance maneuverability. In field-based observations, jerboa trajectories are significantly less predictable than those of quadrupedal rodents, likely increasing predator evasion ability. Consistent with this hypothesis, jerboas exhibit lower anxiety in open fields than quadrupedal rodents, a behavior that varies inversely with predator evasion ability. Our unpredictability metric expands the scope of quantitative biomechanical studies to include non-steady-state locomotion in a variety of evolutionary and ecologically significant contexts.
While prior studies of swallow manoeuvering have focused on slow-speed flight and obstacle avoidance in still air, swallows survive by foraging at high speeds in windy environments. Recent advances in field-portable, high-speed video systems, coupled with precise anemometry, permit measures of high-speed aerial performance of birds in a natural state. We undertook the present study to test: (i) the manner in which barn swallows (Hirundo rustica) may exploit wind dynamics and ground effect while foraging and (ii) the relative importance of flapping versus gliding for accomplishing high-speed manoeuvers. Using multi-camera videography synchronized with wind-velocity measurements, we tracked coursing manoeuvers in pursuit of prey. Wind speed averaged 1.3–2.0 m s−1 across the atmospheric boundary layer, exhibiting a shear gradient greater than expected, with instantaneous speeds of 0.02–6.1 m s−1. While barn swallows tended to flap throughout turns, they exhibited reduced wingbeat frequency, relying on glides and partial bounds during maximal manoeuvers. Further, the birds capitalized on the near-earth wind speed gradient to gain kinetic and potential energy during both flapping and gliding turns; providing evidence that such behaviour is not limited to large, fixed-wing soaring seabirds and that exploitation of wind gradients by small aerial insectivores may be a significant aspect of their aeroecology.
Flying birds rely on visual cues for retinal image stabilization by negating rotation-induced optic flow, the motion of the visual panorama across the retina, through corrective eye and head movements. In combination with vestibular and proprioceptive feedback, birds may also use visual cues to stabilize their body during flight. Here, we test whether artificially induced wide-field motion generated through projected visual patterns elicits maneuvers in body orientation and flight position, in addition to stabilizing vision. To test this hypothesis, we present hummingbirds flying freely within a 1.2 m cylindrical visual arena with a virtual surround rotated at different speeds about its vertical axis. The birds responded robustly to these visual perturbations by rotating their heads and bodies with the moving visual surround, and by adjusting their flight trajectories, following the surround. Thus, similar to insects, hummingbirds appear to use optic flow cues to control flight maneuvers as well as to stabilize their visual inputs.
Skeletal muscles share many common, highly conserved features of organization at the molecular and myofilament levels, giving skeletal muscle fibers generally similar and characteristic mechanical and energetic properties; properties well described by classical studies of muscle mechanics and energetics. However, skeletal muscles can differ considerably in architectural design (fiber length, pinnation, and connective tissue organization), as well as fiber type, and how they contract in relation to the timing of neuromotor activation and in vivo length change. The in vivo dynamics of muscle contraction is, therefore, crucial to assessing muscle design and the roles that muscles play in animal movement. Architectural differences in muscle–tendon organization combined with differences in the phase of activation and resulting fiber length changes greatly affect the time-varying force and work that muscles produce, as well as the energetic cost of force generation. Intrinsic force–length and force–velocity properties of muscles, together with their architecture, also play important roles in the control of movement, facilitating rapid adjustments to changing motor demands. Such adjustments complement slower, reflex-mediated neural feedback control of motor recruitment. Understanding how individual fiber forces are integrated to whole-muscle forces, which are transmitted to the skeleton for producing and controlling locomotor movement, is therefore essential for assessing muscle design in relation to the dynamics of movement.
Airflow conditions close to the Earth's surface are often complex, posing challenges to flight stability and control for volant taxa. Relatively little is known about how well flying animals can contend with complex, adverse air flows, or about the flight-control mechanisms employed by animals to mitigate wind disturbances. Several recent studies have examined flight in the unsteady von Karman vortex streets that form behind cylinders, generating flow disturbances that are predictable in space and time; these structures are relatively rare in nature, as they occur only the immediate, downstream vicinity of an object. In contrast, freestream turbulence is characterized by rapid, unpredictable flow disturbances across a wide range of spatial and temporal scales, and is nearly ubiquitous in natural habitats. Hummingbirds are ideal organisms for studying the influence of freestream turbulence on flight, as they forage in a variety of aerial conditions and are powerful flyers. We filmed ruby-throated hummingbirds (A. colubris) maintaining position at a feeder in laminar and strongly turbulent (intensity approximately 15%) airflow environments within a wind tunnel, and compared their mean head, body, tail and wing kinematics, as well as variability in these parameters. Hummingbirds exhibited remarkably stable head position and orientation in both smooth and turbulent flow while maintaining position at the feeder. However, the hummingbird's body was less stable in turbulent flow and appeared to be most sensitive to disturbances along the mediolateral axis, displaying large lateral accelerations, translations, and rolling motions during flight. The hummingbirds mitigated these disturbances by increasing mean wing stroke amplitude and stroke plane angle, and by varying these parameters asymmetrically between the wings, and from one stroke to the next. They also actively varied the orientation and fan angle of the tail, maintaining a larger mean fan angle when flying in turbulent flow; this may improve their passive stability, but likely incurs an energetic cost due to increased drag. Overall, we observed many of the same kinematic changes noted previously for hummingbirds flying in a von Karman vortex street, but we also observed kinematic changes associated with high force production, similar to those seen during load-lifting or high-speed flight. These findings suggest that flight may be particularly costly in fully mixed, freestream turbulence, the flow condition that hummingbirds are likely to encounter most frequently in natural habitats.
The complexity of low speed maneuvering flight is apparent from the combination of two critical aspects of this behavior: high power and precise control. To understand how such control is achieved, we examined the underlying kinematics and resulting aerodynamic mechanisms of low speed turning flight in the pigeon (Columba livia). Three birds were trained to perform 90 deg level turns in a stereotypical fashion and detailed three-dimensional (3D) kinematics were recorded at high speeds. Applying the angular momentum principle, we used mechanical modeling based on time-varying 3D inertia properties of individual sections of the pigeon's body to separate angular accelerations of the torso based on aerodynamics from those based on inertial effects. Directly measured angular accelerations of the torso were predicted by aerodynamic torques, justifying inferences of aerodynamic torque generation based on inside wing versus outside wing kinematics. Surprisingly, contralateral asymmetries in wing speed did not appear to underlie the 90 deg aerial turns, nor did contralateral differences in wing area, angle of attack, wingbeat amplitude or timing. Instead, torso angular accelerations into the turn were associated with the outside wing sweeping more anteriorly compared with a more laterally directed inside wing. In addition to moving through a relatively more retracted path, the inside wing was also more strongly pronated about its long axis compared with the outside wing, offsetting any difference in aerodynamic angle of attack that might arise from the observed asymmetry in wing trajectories. Therefore, to generate roll and pitch torques into the turn, pigeons simply reorient their wing trajectories toward the desired flight direction. As a result, by acting above the center of mass, the net aerodynamic force produced by the wings is directed inward, generating the necessary torques for turning.
Individuals traversing challenging obstacles are faced with a decision: they can adopt traversal strategies that minimally disrupt their normal locomotion patterns or they can adopt strategies that substantially alter their gait, conferring new advantages and disadvantages. We flew pigeons (Columba livia) through an array of vertical obstacles in a flight arena, presenting them with this choice. The pigeons selected either a strategy involving only a slight pause in the normal wing beat cycle, or a wings-folded posture granting reduced efficiency but greater stability should a misjudgment lead to collision. The more stable but less efficient flight strategy was not used to traverse easy obstacles with wide gaps for passage but came to dominate the postures used as obstacle challenge increased with narrower gaps and there was a greater chance of a collision. These results indicate that birds weigh potential obstacle negotiation strategies and estimate task difficulty during locomotor pattern selection.
In this paper we propose a constitutive model to analyze in-plane extension of goat fascia lata. We first perform a histological analysis of the fascia that shows a well-organized bi-layered arrangement of undulated collagen fascicles oriented along two well defined directions. To develop a model consistent with the tissue structure we identify the absolute and relative thickness of each layer and the orientation of the preferred directions. New data are presented showing the mechanical response in uniaxial and planar biaxial extension. The paper proposes a constitutive relation to describe the mechanical response. We provide a summary of the main ingredients of the nonlinear theory of elasticity and introduce a suitable strain-energy function to describe the anisotropic response of the fascia. We validate the model by showing good fit of the numerical results and the experimental data. Comments are included about differences and analogies between goat fascia lata and the human iliotibial band.
Various flight navigation strategies for birds have been identified at the large spatial scales of migratory and homing behaviours. However, relatively little is known about close-range obstacle negotiation through cluttered environments. To examine obstacle flight guidance, we tracked pigeons (Columba livia) flying through an artificial forest of vertical poles. Interestingly, pigeons adjusted their flight path only approximately 1.5 m from the forest entry, suggesting a reactive mode of path planning. Combining flight trajectories with obstacle pole positions, we reconstructed the visual experience of the pigeons throughout obstacle flights. Assuming proportional-derivative control with a constant delay, we searched the relevant parameter space of steering gains and visuomotor delays that best explained the observed steering. We found that a pigeon's steering resembles proportional control driven by the error angle between the flight direction and the desired opening, or gap, between obstacles. Using this pigeon steering controller, we simulated obstacle flights and showed that pigeons do not simply steer to the nearest opening in the direction of flight or destination. Pigeons bias their flight direction towards larger visual gaps when making fast steering decisions. The proposed behavioural modelling method converts the obstacle avoidance behaviour into a (piecewise) target-aiming behaviour, which is better defined and understood. This study demonstrates how such an approach decomposes open-loop free-flight behaviours into components that can be independently evaluated.
The output of skeletal muscle can be varied by selectively recruiting different motor units. However, our knowledge of muscle function is largely derived from muscle in which all motor units are activated. This discrepancy may limit our understanding of in vivo muscle function. Hence, this study aimed to characterize the mechanical properties of muscle with different motor unit activation. We determined the isometric properties and isotonic force-velocity relationship of rat plantaris muscles in situ with all of the muscle active, 30% of the muscle containing predominately slower motor units active or 20% of the muscle containing predominately faster motor units active. There was a significant effect of active motor unit type on isometric force rise time (p < 0.001) and the force-velocity relationship (p < 0.001). Surprisingly, force rise time was longer and maximum shortening velocity higher when all motor units were active than when either fast or slow motor units were selectively activated. We propose this is due to the greater relative effects of factors such as series compliance and muscle resistance to shortening during sub-maximal contractions. The findings presented here suggest that recruitment according to the size principle, where slow motor units are activated first and faster ones recruited as demand increases, may not pose a mechanical paradox, as has been previously suggested.
Trotting, bipedal running, and especially hopping have long been considered the principal bouncing gaits of legged animals. We use the radial-leg spring constant [Formula: see text] to quantify the stiffness of the physical leg during bouncing gaits. The radial-leg is modeled as an extensible strut between the hip and the ground and [Formula: see text] is determined from the force and deflection of this strut in each instance of stance. A Hookean spring is modeled in-series with a linear actuator and the stiffness of this spring [Formula: see text] is determined by minimizing the work of the actuator while reproducing the measured force-deflection dynamics of an individual leg during trotting or running, and of the paired legs during hopping. Prior studies have estimated leg stiffness using [Formula: see text], a metric that imagines a virtual-leg connected to the center of mass. While [Formula: see text] has been applied extensively in human and comparative biomechanics, we show that [Formula: see text] more accurately models the spring in the leg when actuation is allowed, as is the case in biological and robotic systems. Our allometric analysis of [Formula: see text] in the kangaroo rat, tammar wallaby, dog, goat, and human during hopping, trotting, or running show that [Formula: see text] scales as body mass to the two-third power, which is consistent with the predictions of dynamic similarity and with the scaling of [Formula: see text]. Hence, two-third scaling of locomotor spring constants among mammals is supported by both the radial-leg and virtual-leg models, yet the scaling of [Formula: see text] emerges from work-minimization in the radial-leg model instead of being a defacto result of the ratio of force to length used to compute [Formula: see text]. Another key distinction between the virtual-leg and radial-leg is that [Formula: see text] is substantially greater than [Formula: see text], as indicated by a 30-37% greater scaling coefficient for [Formula: see text]. We also show that the legs of goats are on average twice as stiff as those of dogs of the same mass and that goats increase the stiffness of their legs, in part, by more nearly aligning their distal limb-joints with the ground reaction force vector. This study is the first allometric analysis of leg spring constants in two decades. By means of an independent model, our findings reinforce the two-third scaling of spring constants with body mass, while showing that springs in-series with actuators are stiffer than those predicted by energy-conservative models of the leg.
Fascia is a highly organized collagenous tissue that is ubiquitous in the body, but whose function is not well understood. Because fascia has a sheet-like structure attaching to muscles and bones at multiple sites, it is exposed to different states of multi- or biaxial strain. In order to measure how biaxial strain affects fascia material behavior, planar biaxial tests with strain control were performed on longitudinal and transversely oriented samples of goat fascia lata (FL). Cruciform samples were cycled to multiple strain levels while the perpendicular direction was held at a constant strain. Structural differences among FL layers were examined using histology and SEM. Results show that FL stiffness, hysteresis, and strain energy density are greater in the longitudinal vs. transverse direction. Increased stiffness in the longitudinal layer is likely due to its greater thickness and greater average fibril diameter compared to the transverse layer(s). Perpendicular strain did not affect FL material behavior. Differential loading in the longitudinal vs. transverse directions may lead to structural changes, enhancing the ability of the longitudinal FL to transmit force, store energy, or stabilize the limb during locomotion. The relative compliance of the transverse fibers may allow expansion of underlying muscles when they contract.
We review here the use and reliability of Hill-type muscle models to predict muscle performance under varying conditions, ranging from in situ production of isometric force to in vivo dynamics of muscle length change and force in response to activation. Muscle models are frequently used in musculoskeletal simulations of movement, particularly when applied to studies of human motor performance in which surgically implanted transducers have limited use. Musculoskeletal simulations of different animal species also are being developed to evaluate comparative and evolutionary aspects of locomotor performance. However, such models are rarely validated against direct measures of fascicle strain or recordings of muscle-tendon force. Historically, Hill-type models simplify properties of whole muscle by scaling salient properties of single fibers to whole muscles, typically accounting for a muscle's architecture and series elasticity. Activation of the model's single contractile element (assigned the properties of homogenous fibers) is also simplified and is often based on temporal features of myoelectric (EMG) activation recorded from the muscle. Comparison of standard one-element models with a novel two-element model and with in situ and in vivo measures of EMG, fascicle strain, and force recorded from the gastrocnemius muscles of goats shows that a two-element Hill-type model, which allows independent recruitment of slow and fast units, better predicts temporal patterns of in situ and in vivo force. Recruitment patterns of slow/fast units based on wavelet decomposition of EMG activity in frequency-time space are generally correlated with the intensity spectra of the EMG signals, the strain rates of the fascicles, and the muscle-tendon forces measured in vivo, with faster units linked to greater strain rates and to more rapid forces. Using direct measures of muscle performance to further test Hill-type models, whether traditional or more complex, remains critical for establishing their accuracy and essential for verifying their applicability to scientific and clinical studies of musculoskeletal function.
Animals modulate the power output needed for different locomotor tasks by changing muscle forces and fascicle strain rates. To generate the necessary forces, appropriate motor units must be recruited. Faster motor units have faster activation-deactivation rates than slower motor units, and they contract at higher strain rates; therefore, recruitment of faster motor units may be advantageous for tasks that involve rapid movements or high rates of work. This study identified motor unit recruitment patterns in the gastrocnemii muscles of goats and examined whether faster motor units are recruited when locomotor speed is increased. The study also examined whether locomotor tasks that elicit faster (or slower) motor units are associated with increased (or decreased) in vivo tendon forces, force rise and relaxation rates, fascicle strains and/or strain rates. Electromyography (EMG), sonomicrometry and muscle-tendon force data were collected from the lateral and medial gastrocnemius muscles of goats during level walking, trotting and galloping and during inclined walking and trotting. EMG signals were analyzed using wavelet and principal component analyses to quantify changes in the EMG frequency spectra across the different locomotor conditions. Fascicle strain and strain rate were calculated from the sonomicrometric data, and force rise and relaxation rates were determined from the tendon force data. The results of this study showed that faster motor units were recruited as goats increased their locomotor speeds from level walking to galloping. Slow inclined walking elicited EMG intensities similar to those of fast level galloping but different EMG frequency spectra, indicating that recruitment of the different motor unit types depended, in part, on characteristics of the task. For the locomotor tasks and muscles analyzed here, recruitment patterns were generally associated with in vivo fascicle strain rates, EMG intensity and tendon force. Together, these data provide new evidence that changes in motor unit recruitment have an underlying mechanical basis, at least for certain locomotor tasks.
Hill-type models are commonly used to estimate muscle forces during human and animal movement-yet the accuracy of the forces estimated during walking, running, and other tasks remains largely unknown. Further, most Hill-type models assume a single contractile element, despite evidence that faster and slower motor units, which have different activation-deactivation dynamics, may be independently or collectively excited. This study evaluated a novel, two-element Hill-type model with "differential" activation of fast and slow contractile elements. Model performance was assessed using a comprehensive data set (including measures of EMG intensity, fascicle length, and tendon force) collected from the gastrocnemius muscles of goats during locomotor experiments. Muscle forces predicted by the new two-element model were compared to the forces estimated using traditional one-element models and to the forces measured in vivo using tendon buckle transducers. Overall, the two-element model resulted in the best predictions of in vivo gastrocnemius force. The coefficient of determination, r(2), was up to 26.9% higher and the root mean square error, RMSE, was up to 37.4% lower for the two-element model than for the one-element models tested. All models captured salient features of the measured muscle force during walking, trotting, and galloping (r(2)=0.26-0.51), and all exhibited some errors (RMSE=9.63-32.2% of the maximum in vivo force). These comparisons provide important insight into the accuracy of Hill-type models. The results also show that incorporation of fast and slow contractile elements within muscle models can improve estimates of time-varying, whole muscle force during locomotor tasks.
Goats and other quadrupeds must modulate the work output of their muscles to accommodate the changing mechanical demands associated with locomotion in their natural environments. This study examined which hindlimb joint moments goats use to generate and absorb mechanical energy on level and sloped surfaces over a range of locomotor speeds. Ground reaction forces and the three-dimensional locations of joint markers were recorded as goats walked, trotted and galloped over 0, +15 and -15 deg sloped surfaces. Net joint moments, powers and work were estimated at the goats' hip, knee, ankle and metatarsophalangeal joints throughout the stance phase via inverse dynamics calculations. Differences in locomotor speed on the level, inclined and declined surfaces were characterized and accounted for by fitting regression equations to the joint moment, power and work data plotted versus non-dimensionalized speed. During level locomotion, the net work generated by moments at each of the hindlimb joints was small (less than 0.1 J kg(-1) body mass) and did not vary substantially with gait or locomotor speed. During uphill running, by contrast, mechanical energy was generated at the hip, knee and ankle, and the net work at each of these joints increased dramatically with speed (P<0.05). The greatest increases in positive joint work occurred at the hip and ankle. During downhill running, mechanical energy was decreased in two main ways: goats generated larger knee extension moments in the first half of stance, absorbing energy as the knee flexed, and goats generated smaller ankle extension moments in the second half of stance, delivering less energy. The goats' hip extension moment in mid-stance was also diminished, contributing to the decrease in energy. These analyses offer new insight into quadrupedal locomotion, clarifying how the moments generated by hindlimb muscles modulate mechanical energy at different locomotor speeds and grades, as needed to accommodate the demands of variable terrain.
Biomechanical models of whole muscles commonly used in simulations of musculoskeletal function and movement typically assume that the muscle generates force as a scaled-up muscle fiber. However, muscles are comprised of motor units that have different intrinsic properties and that can be activated at different times. This study tested whether a muscle model comprised of motor units that could be independently activated resulted in more accurate predictions of force than traditional Hill-type models. Forces predicted by the models were evaluated by direct comparison with the muscle forces measured in situ from the gastrocnemii in goats. The muscle was stimulated tetanically at a range of frequencies, muscle fiber strains were measured using sonomicrometry, and the activation patterns of the different types of motor unit were calculated from electromyographic recordings. Activation patterns were input into five different muscle models. Four models were traditional Hill-type models with different intrinsic speeds and fiber-type properties. The fifth model incorporated differential groups of fast and slow motor units. For all goats, muscles and stimulation frequencies the differential model resulted in the best predictions of muscle force. The in situ muscle output was shown to depend on the recruitment of different motor units within the muscle.