Effect of internal damping on locomotion in frictional environments

The gaits of undulating animals arise from a complex interaction of their central nervous system, muscle, connective tissue, bone, and environment. As a simplifying assumption, many previous studies have often assumed that sufficient internal force is available to produce observed kinematics, thus not focusing on quantifying the interconnection between muscle effort, body shape, and external reaction forces. This interplay, however, is critical to locomotion performance in crawling animals, especially when accompanied by body viscoelasticity. Moreover, in bioinspired robotic applications, the body's internal damping is indeed a parameter that the designer can tune. Still, the effect of internal damping is not well understood. This study explores how internal damping affects the locomotion performance of a crawler with a continuous, viscoelastic, nonlinear beam model. Crawler muscle actuation is modeled as a traveling wave of bending moment propagating posteriorly along the body. Consistent with the friction properties of the scales of snakes and limbless lizards, environmental forces are modeled using anisotropic Coulomb friction. It is found that by varying the crawler body's internal damping, the crawler's performance can be altered, and distinct gaits could be achieved, including changing the net locomotion direction from forward to back. We will discuss this forward and backward control and identify the optimal internal damping for peak crawling speed.

Efficient bending and lifting patterns in snake locomotion

We optimize three-dimensional snake kinematics for locomotor efficiency. We assume a general space-curve representation of the snake backbone with small-to-moderate lifting off the ground and negligible body inertia. The cost of locomotion includes work against friction and internal viscous dissipation. When restricted to planar kinematics, our population-based optimization method finds the same types of optima as a previous Newton-based method. With lifting, a few types of optimal motions prevail. We have an s-shaped body with alternating lifting of the middle and ends at small-to-moderate transverse friction. With large transverse friction, curling and sliding motions are typical at small viscous dissipation, replaced by large-amplitude bending at large viscous dissipation. With small viscous dissipation, we find local optima that resemble sidewinding motions across friction coefficient space. They are always suboptimal to alternating lifting motions, with average input power 10–100% higher.

Proprioceptive Mechanism for Bioinspired Fish Swimming

In this Letter, we propose a mechanism for driving bioinspired fish swimming locomotion based on proprioceptive sensing. Proprioception provides information about and representation of a body's position, motion, and acceleration in addition to the usual five senses. We hypothesize that a feedback loop based on this "sixth" sense results in an instability, driving the locomotion. In order to test our assumptions, we use a biomimetic robot and compare the experimental results to a simple yet generic model with excellent agreement.

Continuous body 3-D reconstruction of limbless animals

Limbless animals such as snakes, limbless lizards, worms, eels and lampreys move their slender, long bodies in three dimensions to traverse diverse environments. Accurately quantifying their continuous body's 3-D shape and motion is important for understanding body-environment interactions in complex terrain, but this is difficult to achieve (especially for local orientation and rotation). Here, we describe an interpolation method to quantify continuous body 3-D position and orientation. We simplify the body as an elastic rod and apply a backbone optimization method to interpolate continuous body shape between end constraints imposed by tracked markers. Despite over-simplifying the biomechanics, our method achieves a higher interpolation accuracy (∼50% error) in both 3-D position and orientation compared with the widely used cubic B-spline interpolation method. Beyond snakes traversing large obstacles as demonstrated, our method applies to other long, slender, limbless animals and continuum robots. We provide codes and demo files for easy application of our method.

Experimental-numerical method for calculating bending moments in swimming fish shows that fish larvae control undulatory swimming with simple actuation

Most fish swim with body undulations that result from fluid-structure interactions between the fish's internal tissues and the surrounding water. Gaining insight into these complex fluid-structure interactions is essential to understand how fish swim. To this end, we developed a dedicated experimental-numerical inverse dynamics approach to calculate the lateral bending moment distributions for a large-amplitude undulatory swimmer that moves freely in three-dimensional space. We combined automated motion tracking from multiple synchronised high-speed video sequences, computation of fluid dynamic stresses on the swimmer's body from computational fluid dynamics, and bending moment calculations using these stresses as input for a novel beam model of the body. The bending moment, which represent the system's net actuation, varies over time and along the fish's central axis due to muscle actions, passive tissues, inertia, and fluid dynamics. Our three-dimensional analysis of 113 swimming events of zebrafish larvae ranging in age from 3 to 12 days after fertilisation shows that these bending moment patterns are not only relatively simple but also strikingly similar throughout early development and from fast starts to periodic swimming. This suggests that fish larvae may produce and adjust swimming movements relatively simply, yet effectively, while restructuring their neuromuscular control system throughout their rapid development.

Energetics of optimal undulatory swimming organisms

Energy consumption is one of the primary considerations in animal locomotion. In swimming locomotion, a number of questions related to swimming energetics of an organism and how the energetic quantities scale with body size remain open, largely due to the difficulties with modeling and measuring the power production and consumption. Based on a comprehensive theoretical framework that incorporates cyclic muscle behavior, structural dynamics and swimming hydrodynamics, we perform extensive computational simulations and show that many of the outstanding problems in swimming energetics can be explained by considering the coupling between hydrodynamics and muscle contraction characteristics, as well as the trade-offs between the conflicting performance goals of sustained swimming speed U and cost of transport COT. Our results lead to three main conclusions: (1) in contrast to previous hypotheses, achieving optimal values of U and COT is independent of producing maximal power or efficiency; (2) muscle efficiency in swimming, in contrast to that in flying or running, decreases with increasing body size, consistent with muscle contraction characteristics; (3) the long-standing problem of two disparate patterns of longitudinal power output distributions in swimming fish can be reconciled by relating the two patterns to U-optimal or COT-optimal swimmers, respectively. We also provide further evidence that the use of tendons in caudal regions is beneficial from an energetic perspective. Our conclusions explain and unify many existing observations and are supported by computational data covering nine orders of magnitude in body size.

3D computational models explain muscle activation patterns and energetic functions of internal structures in fish swimming

How muscles are used is a key to understanding the internal driving of fish swimming. However, the underlying mechanisms of some features of the muscle activation patterns and their differential appearance in different species are still obscure. In this study, we explain the muscle activation patterns by using 3D computational fluid dynamics models coupled to the motion of fish with prescribed deformation and examining the torque and power required along the fish body with two primary swimming modes. We find that the torque required by the hydrodynamic forces and body inertia exhibits a wave pattern that travels faster than the curvature wave in both anguilliform and carangiform swimmers, which can explain the traveling wave speeds of the muscle activations. Notably, intermittent negative power (i.e., power delivered by the fluid to the body) on the posterior part, along with a timely transfer of torque and energy by tendons, explains the decrease in the duration of muscle activation towards the tail. The torque contribution from the body elasticity further clarifies the wave speed increase or the reverse of the wave direction of the muscle activation on the posterior part of a carangiform swimmer. For anguilliform swimmers, the absence of the aforementioned changes in the muscle activation on the posterior part is consistent with our torque prediction and the absence of long tendons from experimental observations. These results provide novel insights into the functions of muscles and tendons as an integral part of the internal driving system, especially from an energy perspective, and they highlight the differences in the internal driving systems between the two primary swimming modes.

For undulatory swimming, fish form posteriorly traveling waves of body bending by activating their muscles sequentially along the body. However, experimental observations have shown that the muscle activation wave does not simply match the bending wave. Researchers have previously computed the torque required for muscles along the body based on classic hydrodynamic theories and explained the higher wave speed of the muscle activation compared to the curvature wave. However, the origins of other features of the muscle activation pattern and their variation among different species are still obscure after decades of research. In this study, we use 3D computational fluid dynamics models to compute the spatiotemporal distributions of both the torque and power required for eel-like and mackerel-like swimming. By examining both the torque and power patterns and considering the energy transfer, storage, and release by tendons and body viscoelasticity, we can explain not only the features and variations in the muscle activation patterns as observed from fish experiments but also how tendons and body elasticity save energy. We provide a mechanical picture in which the body shape, body movement, muscles, tendons, and body elasticity of a mackerel (or similar) orchestrate to make swimming efficient.

Snakes partition their body to traverse large steps stably

Many snakes live in deserts, forests and river valleys and traverse challenging 3-D terrain such as rocks, felled trees and rubble, with obstacles as large as themselves and variable surface properties. By contrast, apart from branch cantilevering, burrowing, swimming and gliding, laboratory studies of snake locomotion have focused on locomotion on simple flat surfaces. Here, to begin to understand snake locomotion in complex 3-D terrain, we studied how the variable kingsnake, a terrestrial generalist, traversed a large step of variable surface friction and step height (up to 30% snout-vent length). The snake traversed by partitioning its body into three sections with distinct functions. Body sections below and above the step oscillated laterally on horizontal surfaces for propulsion, whereas the body section in between cantilevered in a vertical plane to bridge the large height increase. As the animal progressed, these three sections traveled down its body, conforming overall body shape to the step. In addition, the snake adjusted the partitioned gait in response to increase in step height and decrease in surface friction, at the cost of reduced speed. As surface friction decreased, body movement below and above the step changed from a continuous lateral undulation with little slip to an intermittent oscillatory movement with much slip, and initial head lift-off became closer to the step. Given these adjustments, body partitioning allowed the snake to be always stable, even when initially cantilevering but before reaching the surface above. Such a partitioned gait may be generally useful for diverse, complex 3-D terrain.

Efficient immersed-boundary lattice Boltzmann scheme for fluid-structure interaction problems involving large solid deformation

A hybrid numerical method which couples the immersed-boundary lattice Boltzmann method with the smoothed point interpolation method (S-PIM) is presented in this paper for the fluid-structure interaction problems involving large solid deformation. In the method, the lattice Boltzmann method is adopted for its advantages in modeling complex fluid flow, the S-PIM is coupled for its robustness in dealing with large solid deformation, and the immersed-boundary method is used for its efficiency in handling the interaction of fluid and solid. In the fluid-solid coupling procedure, a force correction technique based direct-forcing scheme is introduced to enforce nonslip boundary condition with high accuracy, and an averaged dual time stepping scheme is proposed to get stronger robustness of the present method. Numerical experiments are carefully carried out from benchmark problems such as cylinder Couette flow and a beam in a fluid tunnel to more challenging problems such as a flexible beam in the wake of a cylinder and the swimming of a two-dimensional fishlike body. Comparisons of the numerical results with the referenced solutions show that all desirable features of these coupled methods are inherited in the present coupling scheme, and the efficiency of the present method to model such complex problems is verified.

Transition and formation of the torque pattern of undulatory locomotion in resistive force dominated media

In undulatory locomotion, torques along the body are required to overcome external forces from the environment and bend the body. These torques are usually generated by muscles in animals and closely related to muscle activations. In previous studies, researchers observed a single traveling wave pattern of the torque or muscle activation, but the formation of the torque pattern is still not well understood. To elucidate the formation of the torque pattern required by external resistive forces and the transition as kinematic parameters vary, we use simplistic resistive force theory models of self-propelled, steady undulatory locomotors and examine the spatio-temporal variation of the internal torque. We find that the internal torque has a traveling wave pattern with a decreasing speed normalized by the curvature speed as the wave number (the number of wavelengths on the locomotor's body) increases from 0.5 to 1.8. As the wave number increases to 2 and greater values, the torque transitions into a two-wave-like pattern and complex patterns. Using phasor diagram analysis, we reveal that the formation and transitions of the pattern are consequences of the integration and cancellation of force phasors.

Modeling and simulation of fish swimming with active muscles

Our goal is to reproduce the key features of carangiform swimming by modeling muscle functioning using the notion of active distortions, thus emphasizing the kinematical role of muscle, the generation of movement, rather than the dynamical one, the production of force. This approach, already proposed to model the action of muscles in different contexts, is here tested again for the problem of developing an effective and reliable framework to model and simulate swimming. A proper undulatory movement of a fish-like body is reproduced by defining a pattern of distortions, tuned in both space and time, meant to model the muscles activation which produce the flexural motion of body fish; eventually, interactions with the surrounding water yields the desired thrust. Carangiform swimmers have a relatively inflexible anterior body section and a generally flat, flexible posterior section. Because of this configuration, undulations sent rearward along the body attain a significant amplitude only in the posterior section. We compare the performances of different swimming gaits, and we are able to find some important relations between key parameters such as frequencies, wavelength, tail amplitude, and the achieved swim velocity, or the generated thrust, which summarize the swimming performance. In particular, an interesting relation is found between the Strouhal number and the wavelength of muscles activation. We highlight the muscle function during fish locomotion describing the activation of muscles and the relation between the force production and the shortening-lengthening cycle of muscle. We found a great accordance between results and empirical relations, giving an implicit validation of our models.

Automated Planar Tracking the Waving Bodies of Multiple Zebrafish Swimming in Shallow Water

Zebrafish (Danio rerio) is one of the most widely used model organisms in collective behavior research. Multi-object tracking with high speed camera is currently the most feasible way to accurately measure their motion states for quantitative study of their collective behavior. However, due to difficulties such as their similar appearance, complex body deformation and frequent occlusions, it is a big challenge for an automated system to be able to reliably track the body geometry of each individual fish. To accomplish this task, we propose a novel fish body model that uses a chain of rectangles to represent fish body. Then in detection stage, the point of maximum curvature along fish boundary is detected and set as fish nose point. Afterwards, in tracking stage, we firstly apply Kalman filter to track fish head, then use rectangle chain fitting to fit fish body, which at the same time further judge the head tracking results and remove the incorrect ones. At last, a tracklets relinking stage further solves trajectory fragmentation due to occlusion. Experiment results show that the proposed tracking system can track a group of zebrafish with their body geometry accurately even when occlusion occurs from time to time.

Scaling in Free-Swimming Fish and Implications for Measuring Size-at-Time in the Wild

This study was motivated by the need to measure size-at-age, and thus growth rate, in fish in the wild. We postulated that this could be achieved using accelerometer tags based first on early isometric scaling models that hypothesize that similar animals should move at the same speed with a stroke frequency that scales with length-1, and second on observations that the speed of primarily air-breathing free-swimming animals, presumably swimming 'efficiently', is independent of size, confirming that stroke frequency scales as length-1. However, such scaling relations between size and swimming parameters for fish remain mostly theoretical. Based on free-swimming saithe and sturgeon tagged with accelerometers, we introduce a species-specific scaling relationship between dominant tail beat frequency (TBF) and fork length. Dominant TBF was proportional to length-1 (r2 = 0.73, n = 40), and estimated swimming speed within species was independent of length. Similar scaling relations accrued in relation to body mass-0.29. We demonstrate that the dominant TBF can be used to estimate size-at-time and that accelerometer tags with onboard processing may be able to provide size-at-time estimates among free-swimming fish and thus the estimation of growth rate (change in size-at-time) in the wild.

Gait and speed selection in slender inertial swimmers

Inertial swimmers use flexural movements to push water and generate thrust. We quantify this dynamical process for a slender body in a fluid by accounting for passive elasticity and hydrodynamics and active muscular force generation and proprioception. Our coupled elastohydrodynamic model takes the form of a nonlinear eigenvalue problem for the swimming speed and locomotion gait. The solution of this problem shows that swimmers use quantized resonant interactions with the fluid environment to enhance speed and efficiency. Thus, a fish is like an optimized diode that converts a prescribed alternating transverse motion to forward motion. Our results also allow for a broad comparative view of swimming locomotion and provide a mechanistic basis for the empirical relation linking the swimmer's speed U, length L, and tail beat frequency f, given by U/L ~ f [Bainbridge R (1958) J Exp Biol 35:109-133]. Furthermore, we show that a simple form of proprioceptive sensory feedback, wherein local muscle activation is function of body curvature, suffices to drive elastic instabilities associated with thrust production and leads to a spontaneous swimming gait without the need for a central pattern generator. Taken together, our results provide a simple mechanistic view of swimming consistent with natural observations and suggest ways to engineer artificial swimmers for optimal performance.

X-ray motion analysis of the vertebral column during the startle response in striped bass, Morone saxatilis

Whole-body stiffness has a substantial impact on propulsive wave speed during axial undulatory locomotion in fishes. The connective tissues of the vertebral column may contribute to body stiffness, but without mechanical and kinematic analysis it is unclear whether the in vivo range of motion of intervertebral joints (IVJs) is great enough to stress IVJ tissues, thus generating stiffness. The present study used 2D videoradiography and 3D X-ray reconstruction of moving morphology (XROMM) to quantify vertebral kinematics during the startle response in striped bass (Morone saxatilis). X-ray video revealed two distinct patterns of bending: pattern I begins in the abdominal region and then proceeds to maximum IVJ angles in the caudal region, whereas pattern II begins in the cervical region and proceeds to maximum IVJ angles in the abdominal and then the caudal joints. In pattern II bends, the cervical joints exhibit a greater in vivo range of motion than previously reported in other species. XROMM analysis of caudal IVJs suggests primarily lateral bending: mean axial and dorsoventral rotations were less than 2 deg and inconsistent across 51 sequences analyzed from five individuals, whereas mean maximum lateral bending angles were 10.4±3.57 deg. These angles, combined with previous investigations of mechanical properties, reveal that the maximum angles all occur within the neutral zone of bending, indicating that little stress is experienced about the joint. This suggests that the IVJs of striped bass are quite compliant and likely do not contribute significantly to whole-body stiffness or elastic recoil during swimming in vivo.

Analytical insights into optimality and resonance in fish swimming

This paper provides analytical insights into the hypothesis that fish exploit resonance to reduce the mechanical cost of swimming. A simple body-fluid fish model, representing carangiform locomotion, is developed. Steady swimming at various speeds is analysed using optimal gait theory by minimizing bending moment over tail movements and stiffness, and the results are shown to match with data from observed swimming. Our analysis indicates the following: thrust-drag balance leads to the Strouhal number being predetermined based on the drag coefficient and the ratio of wetted body area to cross-sectional area of accelerated fluid. Muscle tension is reduced when undulation frequency matches resonance frequency, which maximizes the ratio of tail-tip velocity to bending moment. Finally, hydrodynamic resonance determines tail-beat frequency, whereas muscle stiffness is actively adjusted, so that overall body-fluid resonance is exploited.

A kinematic model of Kármán gaiting in rainbow trout

A mechanistic understanding of how fishes swim in unsteady flows is challenging despite its prevalence in nature. Previous kinematic studies of fish Kármán gaiting in a vortex street behind a cylinder only report time-averaged measurements, precluding our ability to formally describe motions on a cycle-by-cycle basis. Here we present the first analytical model that describes the swimming kinematics of Kármán gaiting trout with 70-90% accuracy. We found that body bending kinematics can be modelled with a travelling wave equation, which has also been shown to accurately model free-stream swimming kinematics. However, free-stream swimming and Kármán gaiting are separated in the parameter space; the amplitude, wavelength and frequency values of the traveling wave equation are substantially different for each behavior. During Kármán gaiting, the wave is initiated at the body center, which is 0.2L (where L is total body length) further down the body compared with the initiation point in free-stream swimming. The wave travels with a constant speed, which is higher than the nominal flow speed just as in free-stream swimming. In addition to undulation, we observed that Kármán gaiting fish also exhibit substantial lateral translations and body rotations, which can constitute up to 75% of the behavior. These motions are periodic and their frequencies also match the vortex shedding frequency. There is an inverse correlation between head angle and body angle: when the body rotates in one direction, the head of the fish turns into the opposite direction. Our kinematic model mathematically describes how fish swim in vortical flows in real time and provides a platform to better understand the effects of flow variations as well as the contribution of muscle activity during corrective motions.

Energy harvesting from the tail beating of a carangiform swimmer using ionic polymer-metal composites

In this paper, we study energy harvesting from the beating of a biomimetic fish tail using ionic polymer-metal composites. The design of the biomimetic tail is based on carangiform swimmers and is specifically inspired by the morphology of the heterocercal tail of thresher sharks. The tail is constituted of a soft silicone matrix molded in the form of the heterocercal tail and reinforced by a steel beam of rectangular cross section. We propose a modeling framework for the underwater vibration of the biomimetic tail, wherein the tail is assimilated to a cantilever beam with rectangular cross section and heterogeneous physical properties. We focus on base excitation in the form of a superimposed rotation about a fixed axis and we consider the regime of moderately large-amplitude vibrations. In this context, the effect of the encompassing fluid is described through a hydrodynamic function, which accounts for inertial, viscous and convective phenomena. The model is validated through experiments in which the base excitation is systematically varied and the motion of selected points on the biomimetic tail tracked in time. The feasibility of harvesting energy from an ionic polymer-metal composite attached to the vibrating structure is experimentally and theoretically assessed. The response of the transducer is described using a black-box model, where the voltage output is controlled by the rate of change of the mean curvature. Experiments are performed to elucidate the impact of the shunting resistance, the frequency of the base excitation and the placement of the ionic polymer-metal composite on energy harvesting from the considered biomimetic tail.

Emergence of the advancing neuromechanical phase in a resistive force dominated medium

Undulatory locomotion, a gait in which thrust is produced in the opposite direction of a traveling wave of body bending, is a common mode of propulsion used by animals in fluids, on land, and even within sand. As such, it has been an excellent system for discovery of neuromechanical principles of movement. In nearly all animals studied, the wave of muscle activation progresses faster than the wave of body bending, leading to an advancing phase of activation relative to the curvature toward the tail. This is referred to as "neuromechanical phase lags" (NPL). Several multiparameter neuromechanical models have reproduced this phenomenon, but due to model complexity, the origin of the NPL has proved difficult to identify. Here, we use perhaps the simplest model of undulatory swimming to predict the NPL accurately during sand-swimming by the sandfish lizard, with no fitting parameters. The sinusoidal wave used in sandfish locomotion, the friction-dominated and noninertial granular resistive force environment, and the simplicity of the model allow detailed analysis, and reveal the fundamental mechanism responsible for the phenomenon: the combination of synchronized torques from distant points on the body and local traveling torques. This general mechanism should help explain the NPL in organisms in other environments; we therefore propose that sand-swimming could be an excellent system with which to generate and test other neuromechanical models of movement quantitatively. Such a system can also provide guidance for the design and control of robotic undulatory locomotors in complex environments.

Optimal shape and motion of undulatory swimming organisms

Undulatory swimming animals exhibit diverse ranges of body shapes and motion patterns and are often considered as having superior locomotory performance. The extent to which morphological traits of swimming animals have evolved owing to primarily locomotion considerations is, however, not clear. To shed some light on that question, we present here the optimal shape and motion of undulatory swimming organisms obtained by optimizing locomotive performance measures within the framework of a combined hydrodynamical, structural and novel muscular model. We develop a muscular model for periodic muscle contraction which provides relevant kinematic and energetic quantities required to describe swimming. Using an evolutionary algorithm, we performed a multi-objective optimization for achieving maximum sustained swimming speed U and minimum cost of transport (COT)--two conflicting locomotive performance measures that have been conjectured as likely to increase fitness for survival. Starting from an initial population of random characteristics, our results show that, for a range of size scales, fish-like body shapes and motion indeed emerge when U and COT are optimized. Inherent boundary-layer-dependent allometric scaling between body mass and kinematic and energetic quantities of the optimal populations is observed. The trade-off between U and COT affects the geometry, kinematics and energetics of swimming organisms. Our results are corroborated by empirical data from swimming animals over nine orders of magnitude in size, supporting the notion that optimizing U and COT could be the driving force of evolution in many species.

Mechanisms underlying rhythmic locomotion: interactions between activation, tension and body curvature waves

Undulatory animal locomotion arises from three closely related propagating waves that sweep rostrocaudally along the body: activation of segmental muscles by motoneurons (MNs), strain of the body wall, and muscle tension induced by activation and strain. Neuromechanical models that predict the relative propagation speeds of neural/muscle activation, muscle tension and body curvature can reveal crucial underlying control features of the central nervous system and the power-generating mechanisms of the muscle. We provide an analytical explanation of the relative speeds of these three waves based on a model of neuromuscular activation and a model of the body-fluid interactions for leech anguilliform-like swimming. First, we deduced the motoneuron spike frequencies that activate the muscle and the resulting muscle tension during swimming in intact leeches from muscle bending moments. Muscle bending moments were derived from our video-recorded kinematic motion data by our body-fluid interaction model. The phase relationships of neural activation and muscle tension in the strain cycle were then calculated. Our study predicts that the MN activation and body curvature waves have roughly the same speed (the ratio of curvature to MN activation speed ≈0.84), whereas the tension wave travels about twice as fast. The high speed of the tension wave resulting from slow MN activation is explained by the multiplicative effects of MN activation and muscle strain on tension development. That is, the product of two slower waves (activation and strain) with appropriate amplitude, bias and phase can generate a tension wave with twice the propagation speed of the factors. Our study predicts that (1) the bending moment required for swimming is achieved by minimal MN spike frequency, rather than by minimal muscle tension; (2) MN activity is greater in the mid-body than in the head and tail regions; (3) inhibitory MNs not only accelerate the muscle relaxation but also reduce the intrinsic tonus tension during one sector of the swim cycle; and (4) movements of the caudal end are passive during swimming. These predictions await verification or rejection through further experiments on swimming animals.

Mechanisms underlying rhythmic locomotion: body-fluid interaction in undulatory swimming

Swimming of fish and other animals results from interactions of rhythmic body movements with the surrounding fluid. This paper develops a model for the body-fluid interaction in undulatory swimming of leeches, where the body is represented by a chain of rigid links and the hydrodynamic force model is based on resistive and reactive force theories. The drag and added-mass coefficients for the fluid force model were determined from experimental data of kinematic variables during intact swimming, measured through video recording and image processing. Parameter optimizations to minimize errors in simulated model behaviors revealed that the resistive force is dominant, and a simple static function of relative velocity captures the essence of hydrodynamic forces acting on the body. The model thus developed, together with the experimental kinematic data, allows us to investigate temporal and spatial (along the body) distributions of muscle actuation, body curvature, hydrodynamic thrust and drag, muscle power supply and energy dissipation into the fluid. We have found that: (1) thrust is generated continuously along the body with increasing magnitude toward the tail, (2) drag is nearly constant along the body, (3) muscle actuation waves travel two or three times faster than the body curvature waves and (4) energy for swimming is supplied primarily by the mid-body muscles, transmitted through the body in the form of elastic energy, and dissipated into the water near the tail.

Numerical analysis of a unique mode of locomotion: vertical climbing by Pacific lamprey

Pacific lampreys are capable of climbing vertical wetted surfaces through a two-phase (bending and stretching) locomotion mode using the oral disc for adherence. We investigate the physical mechanism and performance of this process by using a continuous beam model. Two mechanisms, one akin to the jumping process and the other related to the fast stretching of the body, have been identified. This locomotion mode may inspire biomimetic designs of anguilliform swimming devices capable of overcoming steep obstacles. By using a genetic algorithm simulation we identify the combination of kinematic parameters corresponding to optimal efficiency (defined as the gravitational potential energy gained in each climbing step divided by the energy spent to activate the motion). These parameters are similar to laboratory observations of lamprey motion, suggesting that this type of locomotion has been optimized for maximum efficiency through evolution.

Interactions between internal forces, body stiffness, and fluid environment in a neuromechanical model of lamprey swimming

Animal movements result from a complex balance of many different forces. Muscles produce force to move the body; the body has inertial, elastic, and damping properties that may aid or oppose the muscle force; and the environment produces reaction forces back on the body. The actual motion is an emergent property of these interactions. To examine the roles of body stiffness, muscle activation, and fluid environment for swimming animals, a computational model of a lamprey was developed. The model uses an immersed boundary framework that fully couples the Navier-Stokes equations of fluid dynamics with an actuated, elastic body model. This is the first model at a Reynolds number appropriate for a swimming fish that captures the complete fluid-structure interaction, in which the body deforms according to both internal muscular forces and external fluid forces. Results indicate that identical muscle activation patterns can produce different kinematics depending on body stiffness, and the optimal value of stiffness for maximum acceleration is different from that for maximum steady swimming speed. Additionally, negative muscle work, observed in many fishes, emerges at higher tail beat frequencies without sensory input and may contribute to energy efficiency. Swimming fishes that can tune their body stiffness by appropriately timed muscle contractions may therefore be able to optimize the passive dynamics of their bodies to maximize peak acceleration or swimming speed.

Propulsion performance of a skeleton-strengthened fin

We examine numerically the performance of a thin foil reinforced by embedded rays resembling the caudal fins of many fishes. In our study, the supporting rays are depicted as nonlinear Euler–Bernoulli beams with three-dimensional deformability. This structural model is then incorporated into a boundary-element hydrodynamic model to achieve coupled fluid–structure interaction simulation. Kinematically, we incorporate both a homocercal mode with dorso-ventral symmetry and a heterocercal mode with dorso-ventral asymmetry. Using the homocercal mode, our results demonstrate that the anisotropic deformability of the ray-reinforced fin significantly increases its capacity of force generation. This performance enhancement manifests as increased propulsion efficiency, reduced transverse force and reduced sensitivity to kinematic parameters. Further reduction in transverse force is observed by using the heterocercal mode. In the heterocercal model, the fin also generates a small lifting force, which may be important in vertical maneuvers. Via three-dimensional flow visualization, a chain of vortex rings is observed in the wake. Detailed features of the wake, e.g. the orientation of the vortex rings in the heterocercal mode, agree with predictions based upon particle image velocimetry (PIV) measurements of flow around live fish.

Automated visual tracking for studying the ontogeny of zebrafish swimming

The zebrafish Danio rerio is a widely used model organism in studies of genetics, developmental biology, and recently, biomechanics. In order to quantify changes in swimming during all stages of development, we have developed a visual tracking system that estimates the posture of fish. Our current approach assumes planar motion of the fish, given image sequences taken from a top view. An accurate geometric fish model is automatically designed and fit to the images at each time frame. Our approach works across a range of fish shapes and sizes and is therefore well suited for studying the ontogeny of fish swimming, while also being robust to common environmental occlusions. Our current analysis focuses on measuring the influence of vertebra development on the swimming capabilities of zebrafish. We examine wild-type zebrafish and mutants with stiff vertebrae (stocksteif) and quantify their body kinematics as a function of their development from larvae to adult (mutants made available by the Hubrecht laboratory, The Netherlands). By tracking the fish, we are able to measure the curvature and net acceleration along the body that result from the fish's body wave. Here, we demonstrate the capabilities of the tracking system for the escape response of wild-type zebrafish and stocksteif mutant zebrafish. The response was filmed with a digital high-speed camera at 1500 frames s–1. Our approach enables biomechanists and ethologists to process much larger datasets than possible at present. Our automated tracking scheme can therefore accelerate insight in the swimming behavior of many species of (developing)fish.

Limbless undulatory propulsion on land

We analyze the lateral undulatory motion of a natural or artificial snake or other slender organism that "swims" on land by propagating retrograde flexural waves. The governing equations for the planar lateral undulation of a thin filament that interacts frictionally with its environment lead to an incomplete system. Closures accounting for the forces generated by the internal muscles and the interaction of the filament with its environment lead to a nonlinear boundary value problem, which we solve using a combination of analytical and numerical methods. We find that the primary determinant of the shape of the organism is its interaction with the external environment, whereas the speed of the organism is determined primarily by the internal muscular forces, consistent with prior qualitative observations. Our model also allows us to pose and solve a variety of optimization problems such as those associated with maximum speed and mechanical efficiency, thus defining the performance envelope of this mode of locomotion.

An elastic rod model for anguilliform swimming

We develop a model for anguilliform (eel-like) swimming as an elastic rod actuated via time-dependent intrinsic curvature and subject to hydrodynamic drag forces, the latter as proposed by Taylor (in Proc Roy Proc Lond A 214:158-183, 1952). We employ a eometrically exact theory and discretize the resulting nonlinear partial differential evolution both to perform numerical simulations, and to compare with previous models consisting of chains of rigid links or masses connected by springs, dampers, and prescribed force generators representing muscles. We show that muscle activations driven by motoneuronal spike trains via calcium dynamics produce intrinsic curvatures corresponding to near-sinusoidal body shapes in longitudinally-uniform rods, but that passive elasticity causes Taylor's assumption of prescribed shape to fail, leading to time-periodic motions and lower speeds than those predicted Taylor (in Proc Roy Proc Lond A 214:158-183, 1952). We investigate the effects of bending stiffness, body geometry, and activation patterns on swimming speed, turning behavior, and acceleration to steady swimming. We show that laterally-uniform activation yields stable straight swimming and laterally differential activation levels lead to stable turns, and we argue that tapered bodies with reduced caudal (tail-end) activation (to produce uniform intrinsic curvature) swim faster than ones with uniform activation.

Conservation rules, their breakdown, and optimality in Caenorhabditis sinusoidal locomotion

Undulatory locomotion is common to nematodes as well as to limbless vertebrates, but its control is not understood in spite of the identification of hundred of genes involved in Caenorhabditis elegans locomotion. To reveal the mechanisms of nematode undulatory locomotion, we quantitatively analysed the movement of C. elegans with genetic perturbations to neurons, muscles, and skeleton (cuticle). We also compared locomotion of different Caenorhabditis species. We constructed a theoretical model that combines mechanics and biophysics, and that is constrained by the observations of propulsion and muscular velocities, as well as wavelength and amplitude of undulations. We find that normalized wavelength is a conserved quantity among wild-type C. elegans individuals, across mutants, and across different species. The velocity of forward propulsion scales linearly with the velocity of the muscular wave and the corresponding slope is also a conserved quantity and almost optimal; the exceptions are in some mutants affecting cuticle structure. In theoretical terms, the optimality of the slope is equivalent to the exact balance between muscular and visco-elastic body reaction bending moments. We find that the amplitude and frequency of undulations are inversely correlated and provide a theoretical explanation for this fact. These experimental results are valid both for young adults and for all larval stages of wild-type C. elegans. In particular, during development, the amplitude scales linearly with the wavelength, consistent with our theory. We also investigated the influence of substrate firmness on motion parameters, and found that it does not affect the above invariants. In general, our biomechanical model can explain the observed robustness of the mechanisms controlling nematode undulatory locomotion.

Adaptive bone formation in acellular vertebrae of sea bass(Dicentrarchus labrax L.)

Mammalian bone is an active tissue in which osteoblasts and osteoclasts balance bone mass. This process of adaptive modelling and remodelling is probably regulated by strain-sensing osteocytes. Bone of advanced teleosts is acellular yet, despite the lack of osteocytes, it is capable of an adaptive response to physical stimuli. Strenuous exercise is known to induce lordosis. Lordosis is a ventrad curvature of the vertebral column, and the affected vertebrae show an increase in bone formation.

The effects of lordosis on the strain distribution in sea bass(Dicentrarchus labrax L.) vertebrae are assessed using finite element modelling. The response of the local tissue is analyzed spatially and ontogenetically in terms of bone volume.

Lordotic vertebrae show a significantly increased strain energy due to the increased load compared with normal vertebrae when loaded in compression. High strain regions are found in the vertebral centrum and parasagittal ridges. The increase in strain energy is attenuated by a change in architecture due to the increased bone formation. The increased bone formation is seen mainly at the articular surfaces of the vertebrae, although some extra bone is formed in the vertebral centrum.

Regions in which the highest strains are found do not spatially correlate with regions in which the most extensive bone apposition occurs in lordotic vertebrae of sea bass. Mammalian-like strain-regulated bone modelling is probably not the guiding mechanism in adaptive bone modelling of acellular sea bass vertebrae. Chondroidal ossification is found at the articular surfaces where it mediates a rapid adaptive response, potentially attenuating high stresses on the dorsal zygapophyses.

Force transmission via axial tendons in undulating fish: a dynamic analysis

Sonomicrometrics of in vivo axial strain of muscle has shown that the swimming fish body bends like a homogenous, continuous beam in all species except tuna. This simple beam-like behavior is surprising because the underlying tendon structure, muscle structure and behavior are complex. Given this incongruence, our goal was to understand the mechanical role of various myoseptal tendons. We modeled a pumpkinseed sunfish, Lepomis gibbosus, using experimentally-derived physical and mechanical attributes, swimming from rest with steady muscle activity. Axially oriented muscle-tendons, transverse and axial myoseptal tendons, as suggested by current morphological knowledge, interacted to replicate the force and moment distribution. Dynamic stiffness and damping associated with muscle activation, realistic muscle force generation, and force distribution following tendon geometry were incorporated. The vertebral column consisted of 11 rigid vertebrae connected by joints that restricted bending to the lateral plane and endowed the body with its passive viscoelasticity. In reaction to the acceleration of the body in an inviscid fluid and its internal transmission of moment via the vertebral column, the model predicted the kinematic response. Varying only tendon geometry and stiffness, four different simulations were run. Simulations with only intrasegmental tendons produced unstable axial and lateral tail forces and body motions. Only the simulation that included both intra- and intersegmental tendons, muscle-enhanced segment stiffness, and a stiffened caudal joint produced stable and large lateral and axial forces at the tail. Thus this model predicts that axial tendons function within a myomere to (1) convert axial force to moment (moment transduction), (2) transmit axial forces between adjacent myosepta (segment coupling), and, intersegmentally, to (3) distribute axial forces (force entrainment), and (4) stiffen joints in bending (flexural stiffening). The fact that all four functions are needed to produce the most realistic swimming motions suggests that axial tendons are essential to the simple beam-like behavior of fish.

Computational hydrodynamics of animal swimming: boundary element method and three-dimensional vortex wake structure

The slender body theory, lifting surface theories, and more recently panel methods and Navier-Stokes solvers have been used to study the hydrodynamics of fish swimming. This paper presents progress on swimming hydrodynamics using a boundary integral equation method (or boundary element method) based on potential flow model. The unsteady three-dimensional BEM code 3DynaFS that we developed and used is able to model realistic body geometries, arbitrary movements, and resulting wake evolution. Pressure distribution over the body surface, vorticity in the wake, and the velocity field around the body can be computed. The structure and dynamic behavior of the vortex wakes generated by the swimming body are responsible for the underlying fluid dynamic mechanisms to realize the high-efficiency propulsion and high-agility maneuvering. Three-dimensional vortex wake structures are not well known, although two-dimensional structures termed 'reverse Karman Vortex Street' have been observed and studied. In this paper, simulations about a swimming saithe (Pollachius virens) using our BEM code have demonstrated that undulatory swimming reduces three-dimensional effects due to substantially weakened tail tip vortex, resulting in a reverse Karman Vortex Street as the major flow pattern in the three-dimensional wake of an undulating swimming fish.

Muscle dynamics in skipjack tuna: timing of red muscle shortening in relation to activation and body curvature during steady swimming

Cyclic length changes in the internal red muscle of skipjack tuna (Katsuwonus pelamis) were measured using sonomicrometry while the fish swam in a water tunnel at steady speeds of 1.1-2.3 L s(−)(1), where L is fork length. These data were coupled with simultaneous electromyographic (EMG) recordings. The onset of EMG activity occurred at virtually the same phase of the strain cycle for muscle at axial locations between approximately 0.4L and 0.74L, where the majority of the internal red muscle is located. Furthermore, EMG activity always began during muscle lengthening, 40–50 prior to peak length, suggesting that force enhancement by stretching and net positive work probably occur in red muscle all along the body. Our results support the idea that positive contractile power is derived from all the aerobic swimming muscle in tunas, while force transmission is provided primarily by connective tissue structures, such as skin and tendons, rather than by muscles performing negative work. We also compared measured muscle length changes with midline curvature (as a potential index of muscle strain) calculated from synchronised video image analysis. Unlike contraction of the superficial red muscle in other fish, the shortening of internal red muscle in skipjack tuna substantially lags behind changes in the local midline curvature. The temporal separation of red muscle shortening and local curvature is so pronounced that, in the mid-body region, muscle shortening at each location is synchronous with midline curvature at locations that are 7–8 cm (i.e. 8–10 vertebral segments) more posterior. These results suggest that contraction of the internal red muscle causes deformation of the body at more posterior locations, rather than locally. This situation represents a unique departure from the model of a homogeneous bending beam, which describes red muscle strain in other fish during steady swimming, but is consistent with the idea that tunas produce thrust by motion of the caudal fin rather than by undulation of segments along the body.