We developed a three-dimensional, computational biomechanical model of a juvenile Nile crocodile (Crocodylus niloticus) pelvis and hindlimb, composed of 47 pelvic limb muscles, to investigate muscle function. We tested whether crocodiles, which are known to use a variety of limb postures during movement, use limb orientations (joint angles) that optimise the moment arms (leverages) or moment-generating capacities of their muscles during different limb postures ranging from a high walk to a sprawling motion. We also describe the three-dimensional (3D) kinematics of the crocodylian hindlimb during terrestrial locomotion across an instrumented walkway and a treadmill captured via X-ray Reconstruction of Moving Morphology (biplanar fluoroscopy; 'XROMM'). We reconstructed the 3D positions and orientations of each of the hindlimb bones and used dissection data for muscle lines of action to reconstruct a focal, subject-specific 3D musculoskeletal model. Motion data for different styles of walking (a high, crouched, bended and two types of sprawling motion) were fed into the 3D model to identify whether any joints adopted near-optimal poses for leverage across each of the behaviours. We found that (1) the hip adductors and knee extensors had their largest leverages during sprawling postures and (2) more erect postures typically involved greater peak moment arms about the hip (flexion-extension), knee (flexion) and metatarsophalangeal (flexion) joints. The results did not fully support the hypothesis that optimal poses are present during different locomotory behaviours because the peak capacities were not always reached around mid-stance phase. Furthermore, we obtained few clear trends for isometric moment-generating capacities. Therefore, perhaps peak muscular leverage in Nile crocodiles is instead reached either in early/late stance or possibly during swing phase or other locomotory behaviours that were not studied here, such as non-terrestrial movement. Alternatively, our findings could reflect a trade-off between having to execute different postures, meaning that hindlimb muscle leverage is not optimised for any singular posture or behaviour. Our model, however, provides a comprehensive set of 3D estimates of This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

The arrangement and physiology of muscle fibres can strongly influence musculoskeletal function and whole-organismal performance. However, experimental investigation of muscle function during in vivo activity is typically limited to relatively few muscles in a given system. Computational models and simulations of the musculoskeletal system can partly overcome these limitations, by exploring the dynamics of muscles, tendons and other tissues in a robust and quantitative fashion. Here, a high-fidelity, 26-degree-of-freedom musculoskeletal model was developed of the hindlimb of a small ground bird, the elegant-crested tinamou (Eudromia elegans, ~550 g), including all the major muscles of the limb (36 actuators per leg). The model was integrated with biplanar fluoroscopy (XROMM) and forceplate data for walking and running, where dynamic optimization was used to estimate muscle excitations and fibre length changes throughout both gaits. Following this, a series of static simulations over the total range of physiological limb postures were performed, to circumscribe the bounds of possible variation in fibre length. During gait, fibre lengths for all muscles remained between 0.5 to 1.21 times optimal fibre length, but operated mostly on the ascending limb and plateau of the active force-length curve, a result that parallels previous experimental findings for birds, humans and other species. However, the ranges of fibre length varied considerably among individual muscles, especially when considered across the total possible range of joint excursion. Net length change of muscle–tendon units was mostly less than optimal fibre length, sometimes markedly so, suggesting that approaches that use muscle–tendon length change to estimate optimal fibre length in extinct species are likely underestimating this important parameter for many muscles. The results of this study clarify and broaden understanding of muscle function in extant animals, and can help refine approaches used to study extinct species.

To capture and swallow food on land, a sticky tongue supported by the hyoid and gill arch skeleton has evolved in land vertebrates from aquatic ancestors that used mouth-cavity-expanding actions of the hyoid to suck food into the mouth. However, the evolutionary pathway bridging this drastic shift in feeding mechanism and associated hyoid motions remains unknown. Modern fish that feed on land may help to unravel the physical constraints and biomechanical solutions that led to terrestrialization of fish-feeding systems. Here, we show that the mudskipper emerges onto land with its mouth cavity filled with water, which it uses as a protruding and retracting 'hydrodynamic tongue' during the initial capture and subsequent intra-oral transport of food. Our analyses link this hydrodynamic action of the intra-oral water to a sequence of compressive and expansive cranial motions that diverge from the general pattern known for suction feeding in fishes. However, the hyoid motion pattern showed a remarkable resemblance to newts during tongue prehension. Consequently, although alternative scenarios cannot be excluded, hydrodynamic tongue usage may be a transitional step onto which the evolution of adhesive mucosa and intrinsic lingual muscles can be added to gain further independence from water for terrestrial foraging.

Many animals generate propulsive forces by coordinating legs, which contact and push against the surroundings, with bending of the body, which can only indirectly influence these forces. Such body–leg coordination is not commonly employed in quadrupedal robotic systems. To elucidate the role of back bending during quadrupedal locomotion, we study a model system: the salamander, a sprawled-posture quadruped that uses lateral bending of the elongate back in conjunction with stepping of the limbs during locomotion. We develop a geometric approach that yields a low-dimensional representation of the body and limb contributions to the locomotor performance quantified by stride displacement. For systems where the damping forces dominate inertial forces, our approach offers insight into appropriate coordination patterns, and improves the computational efficiency of optimization techniques. In particular, we demonstrate effect of the lateral undulation coordinated with leg movement in the forward, rotational, and lateral directions of the robot motion. We validate the theoretical results using numerical simulations, and then successfully test these approaches using robophysical experiments on granular media, a model deformable, frictional substrate. Although our focus lies primarily on robotics, we also demonstrate that our tools can accurately predict optimal body bending of a living salamander Salamandra salamandra.

Many quadrupedal animals use lateral degrees of freedom in their backs to assist locomotion. This paper seeks to use a robotic model to demonstrate that back bending assists not only forward motion, but also lateral and turning motions. We present a simple planner that uses geometric mechanics to prescribe gaits that coordinate both leg movements and back bending motion. Using these geometric tools, we show that back-bending can improve stride displacement in the forward, rotational, and lateral directions. Interestingly, we observe that the animal's (salamander, Salamandra salamandra) back-bending is close to our calculated back-bending to improve forward displacement. In addition to locomotion performance improvement, back bending can also expand the target position space a robot can reach within one gait cycle. Our results are verified by conducting experiments with a robot moving on granular materials.