Module 3

Snake Locomotion

Without limbs, snakes deploy four distinct gaits — lateral undulation, concertina, rectilinear, and sidewinding — each suited to a different substrate. The underlying principle is anisotropic ventral-scale friction: the belly slides forward more easily than sideways or backward, turning body waves into net thrust. This module quantifies each gait and the biophysics that makes limbless locomotion competitive.

1. Four Gaits (Jayne 1986)

Lateral undulation: a sinusoidal wave travels head-to-tail; the body pushes against substrate asperities (rocks, grass stems) or uses friction anisotropy on flat ground. The commonest gait.
Concertina: the snake anchors posterior coils, extends the anterior body, anchors anteriorly, then hauls the posterior in. Used in tight burrows and up tree branches.
Rectilinear: large ventral scales walk the body forward via coordinated costocutaneous muscle contractions. Slow but used by heavy pythons and boas on open ground.
Sidewinding: only two body segments touch the ground at any instant; the snake rolls a static-contact pattern sideways, leaving characteristic J-shaped tracks. Enables movement on slip-prone dune sand (Crotalus cerastes, Cerastes cerastes).

Jayne’s 1986 electromyographic study of corn snakes (Pantherophis guttatus) resolved the muscle-firing patterns of each gait; Marvi 2014 (Science) characterised the sidewinder’s tangential-force redistribution on sand.

2. Friction Anisotropy

Hu 2009 (PNAS) measured the ventral-scale friction coefficients of a corn snake on fabric: μ_forward ≈ 0.11, μ_lateral≈ 0.38, μ_backward ≈ 0.41. The ratio ~4× is enough to convert a purely lateral body wave into net forward thrust on flat ground.

\[ \mu_{eff}(\theta) = \sqrt{\mu_f^2\cos^2\theta + \mu_s^2\sin^2\theta},\qquad v_{cm} \approx \frac{c}{2}\cdot\frac{\mu_s - \mu_f}{\mu_s + \mu_f} \]

where c is the body-wave speed. The microstructure that produces this anisotropy is a directional array of microscale tribological denticles on each ventral scale (Baum 2014); mimicking the pattern on synthetic surfaces yields low-friction anisotropic tapes used in soft robotics.

Simulation: Undulation & Friction Polar

Body midline snapshots of a travelling sinusoid plus the elliptical friction polar plot that explains why lateral undulation produces net forward motion.

Python

script.py55 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# Snake lateral undulation: body midline as travelling sinusoid
# s is arclength along body, t is time
L = 1.0                          # body length (m)
N = 200                          # number of segments
s = np.linspace(0, L, N)
k = 2 * np.pi / 0.3              # wavelength = 0.3 m
A = 0.08                         # amplitude (m)
omega = 2 * np.pi * 1.5          # 1.5 Hz

fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 6), facecolor='#0a0a1a')
for ax in (ax1, ax2):
    ax.set_facecolor('#111827')
    ax.tick_params(colors='#cbd5e1')
    for s_ in ax.spines.values(): s_.set_color('#334155')

# (1) Body shape snapshots over a cycle
for i, t in enumerate(np.linspace(0, 2*np.pi/omega, 6)):
    x = s
    y = A * np.sin(k*s - omega*t) + 0.06 * i
    ax1.plot(x, y, color=plt.cm.viridis(i/6), lw=2.2)
ax1.set_xlabel('X (m)', color='#cbd5e1')
ax1.set_ylabel('Y (offset snapshots)', color='#cbd5e1')
ax1.set_title('Lateral undulation: travelling wave', color='#6ee7b7', fontweight='bold')

# (2) Friction anisotropy: mu_forward << mu_side
mu_f = 0.11
mu_s = 0.38
mu_b = 0.41
angles = np.linspace(0, 2*np.pi, 200)
# Elliptical friction ~ Jayne 1986 / Hu 2009
mu_eff = np.sqrt((mu_f*np.cos(angles))**2 + (mu_s*np.sin(angles))**2)
ax2.plot(np.cos(angles)*mu_eff, np.sin(angles)*mu_eff, color='#34d399', lw=2.5)
for ang, lab, col in [(0, 'forward\nmu=0.11', '#fbbf24'),
                       (np.pi/2, 'lateral\nmu=0.38', '#f87171'),
                       (np.pi, 'back\nmu=0.41', '#60a5fa')]:
    ax2.plot(np.cos(ang)*mu_eff[int(ang/(2*np.pi)*200) % 200],
             np.sin(ang)*mu_eff[int(ang/(2*np.pi)*200) % 200],
             'o', color=col, markersize=12)
    ax2.annotate(lab, (np.cos(ang)*0.44, np.sin(ang)*0.44), color=col, ha='center', fontweight='bold')
ax2.set_aspect('equal')
ax2.set_title('Belly-scale friction polar plot\n(Hu 2009)', color='#6ee7b7', fontweight='bold')
ax2.set_xlim(-0.5, 0.5); ax2.set_ylim(-0.5, 0.5)
ax2.grid(True, color='#334155', alpha=0.3)

plt.tight_layout()
plt.savefig('output.png', dpi=120, bbox_inches='tight', facecolor='#0a0a1a')
speed = omega / k
print(f'Wave speed c = omega/k = {speed:.2f} m/s')
print(f'Ratio mu_side / mu_forward = {mu_s/mu_f:.2f}')
print(f'Corn-snake speed ~ 0.5 x wave speed = {0.5*speed:.2f} m/s (Hu 2009)')

Click Run to execute the Python code

Code will be executed with Python 3 on the server

3. Arboreal Cantilever Limits

Arboreal snakes cross gaps by cantilevering: extending the body as an unsupported beam until anchor contact. The maximum cantilever length scales with muscle cross-sectional area and body-wall second moment of area:

\[ L_{max} \propto \left(\frac{\sigma_{max}\,I}{\rho g A}\right)^{1/3} \]

Lillywhite 2000 and Jayne 2019 measured gap-crossing in Corallus, Boiga, and Chondropython; most species cantilever 30–50% of body length before tail prehensile anchoring becomes essential. Sidelong arboreal snakes (Chrysopelea) exploit passive gliding by flattening the body into a concave airfoil (Socha 2005, Nature).

4. Aquatic & Burrowing Modes

Sea snakes (Hydrophiinae) use anguilliform swimming with a flattened, oar-like tail. Propulsive efficiency is comparable to eels; Graham 1971 measured ~30–40% efficiency with Strouhal number St = 2fA/U ≈ 0.3, within the 0.2–0.4 optimum that Triantafyllou 1993 derived for oscillating-foil thrust. Burrowing sand boas and blind snakes use concertina within loose sand, generating tunnel-stabilising body-wall pressures far higher than surface sliders.

Key References

• Jayne, B. C. (1986). “Kinematics of terrestrial snake locomotion.” Copeia, 1986, 915–927.

• Hu, D. L., Nirody, J., Scott, T. & Shelley, M. J. (2009). “The mechanics of slithering locomotion.” Proc. Natl. Acad. Sci., 106, 10081–10085.

• Marvi, H. et al. (2014). “Sidewinding with minimal slip: snake and robot ascent of sandy slopes.” Science, 346, 224–229.

• Baum, M. J. et al. (2014). “Anisotropic friction of the ventral scales of the snake Lampropeltis getula californiae.” Tribol. Lett., 54, 139–150.

• Socha, J. J., O’Dempsey, T. & LaBarbera, M. (2005). “A 3-D kinematic analysis of gliding in a flying snake, Chrysopelea paradisi.” J. Exp. Biol., 208, 1817–1833.

• Lillywhite, H. B., LaFrentz, J. R., Lin, Y. C. & Tu, M. C. (2000). “The cantilever abilities of snakes.” J. Herpetol., 34, 523–528.

Share:X Reddit LinkedIn

← Module 2 Module 4: Sensing →