Snakes' bodies are covered in scales that make it easier to slide in some directions than in others.This frictional anisotropy allows for sliding locomotion with an undulatory gait, one of the most common for snakes. Isotropic friction is a simpler situation (that arises with snake robots for example) but is less understood. In this work we regularize a model for sliding locomotion to allow for static friction. We then propose a robust iterative numerical method to study the efficiency of a wide range of motions under isotropic Coulomb friction. We find that simple undulatory motions give little net locomotion in the isotropic regime. We compute general time-harmonic motions of three-link bodies and find three local optima for efficiency. The top two involve static friction to some extent. We then propose a class of smooth body motions that have similarities to concertina locomotion (including the involvement of static friction) and can achieve optimal efficiency for both isotropic and anisotropic friction.
I. INTRODUCTIONSnake locomotion has attracted the interest of biologists and engineers for several decades [1][2][3][4][5][6]. Many locomoting animals use appendages such as legs, wings, or fins to exert a force on the substrate or surrounding fluid, and propel the rest of the body forward [7]. Snakes lack appendages, and thus it is less clear which parts of the snake body should exert propulsive forces, and at which instants during the motion, to move forward efficiently.A typical way to understand how organisms move is to study physical or computational models and compare their motions with those of the actual organisms [7][8][9][10][11]. One can take a step further and pose and solve optimization problems for the models. This can suggest locomotion strategies that are effective for man-made vehicles [12][13][14]. It can also help understand why organisms have evolved in particular ways under a multitude of constraints [15][16][17].Often what is optimized is a measure of the efficiency of locomotion. For example, one can maximize the average speed for a given time-averaged power expended by the organism. One can study the effects of physical parameters and constraints by varying them and studying how the optimal solutions change. Well-known examples are optimization studies of organisms moving in low-[18-24] and high-Reynolds-number fluid flows [25][26][27][28][29][30][31]. For locomotion in frictional (terrestrial or granular) media, frictional forces can result in distinctive modes of efficient (or optimal) locomotion [32,33].Snakes are limbless reptiles with elongated bodies, supported by a backbone with 100-500 bony segments (vertebrae) * Electronic address: alben@umich.edu arXiv:1902.03163v1 [physics.bio-ph] 8 Feb 2019 [34]. The vertebrae allow for high flexibility particularly in the lateral (side-to-side) direction, with less flexibility for vertical (dorso-ventral) bending or for torsion. Running along the backbone are muscles that attach to the sides of the vertebrae and cause bending. The snake bod...