Abstract
Peristalsis, i.e., a motion pattern arising from the propagation of muscle contraction and expansion waves along the body, is a common locomotion strategy for limbless animals. Mimicking peristalsis in bio-inspired robots has attracted considerable attention in the literature. It has recently been observed that maximal velocity in a metameric earthworm-like robot is achieved by actuating the segments using a “phase coordination” principle. This paper shows that, in fact, peristalsis (which requires not only phase coordination, but also that all segments oscillate at same frequency and amplitude) emerges from optimization principles. More precisely, basing our analysis on the assumption of small deformations, we show that peristaltic waves provide the optimal actuation solution in the ideal case of a periodic infinite system, and that this is approximately true, modulo edge effects, for the real, finite length system. Therefore, this paper confirms the effectiveness of mimicking peristalsis in bio-inspired robots, at least in the small-deformation regime. Further research will be required to test the effectiveness of this strategy if large deformations are allowed.
1. Introduction
The study of self-propelled locomotors exploiting friction-induced traction as a result of body shape changes, is gaining attention because of the variety of physical systems which take advantage of such a locomotion strategy. One motivation is the desire to understand biological phenomena, such as cell migration on or within solid substrates, matrices, and tissues (Alberts et al., 2002). Another motivation is the attempt to replicate these mechanisms in robotics with the idea that biomimetic constructs may outperform traditional ones when confronted with unstructured and unpredictable environments.
In particular, robotic locomotion research has recently considered crawling and burrowing animals (e.g., earthworms, snakes, and caterpillars), whence an increasing number of research projects on bio-inspired metameric (soft) robots (Menciassi et al., 2006; Wang et al., 2009; Boxerbaum et al., 2012; Daltorio et al., 2013; Fang et al., 2015; Nemitz et al., 2016; Umedachi et al., 2016; Ge et al., 2017). As a matter of fact, many species such as earthworms, caterpillars, sea cucumbers and snails move using peristalsis which is a locomotion mechanism consisting of a series of wave-like muscle relaxation and contraction which propagate along the body (Quillin, 1999). One of the most studied biological species is Lumbricus terrestris (commonly known as nightcrawler) which is a kind of earthworm which uses peristalsis both for surface crawling and for burrowing. Each of its metameres (body segments) is endowed with longitudinal and circular muscles and can regulate frictional forces thanks to microscopic bristles called setae (Quillin, 1999). Understanding how relatively simple organisms are able to attain peristalsis and to which extent coordination is regulated by either the nervous system or spontaneous reflexes, are questions addressed by researchers for about a century and are still drawing attention (Garrey and Moore, 1915; Gray and Lissmann, 1938; Gardner, 1976; Quillin, 1999).
In the field of robotics, peristalsis has been mostly mimicked by a priori assignment of “gaits” defined by a few scalar parameters. Optimization of locomotion performances with respect to variations of these scalar parameters has been studied. Fang et al. (2015) consider harmonic deformations with a single, fixed, (time) frequency and amplitude, and determine the phase patterns of actuation maximizing the average velocity. Optimization leads to phase coordination, in the form of a pattern which is close to the identical-phase-difference (IPD) pattern corresponding to peristalsis. However, no rigorous proof of the connection between peristaltic waves and optimal actuation is given and, more importantly, the basic hypothesis of harmonic oscillations with a single fixed time frequency and amplitude is taken as an a priori assumption.
This paper aims to provide a deeper understanding of harmonic oscillations and peristalsis as result of an optimization problem rather than an a priori hypothesis. Indeed, we prove that - in the regime of small deformations - peristalsis is a symmetry property of the solution to an optimization problem. Symmetry of the solution comes from symmetry properties of operators in the equations governing the optimization problem, which are, in turn, the signature of geometric symmetries of the physical system.
The rest of the paper is organized in three main sections: material and methods, results and discussion. The first one, inspired by nightcrawlers' retractable setae, introduces a velocity-force law which is able to describe, for limit values of a single scalar parameter, the linear case (Newtonian) as well as the case of “free slip–perfect grip.” This model is presented in both continuous and discrete version and in the latter case we address some related optimal control problems. The second section illustrates the behavior of the continuous model by means of two examples and shows how peristalsis emerges, up to edge-effects, from optimization problems in the discrete framework. The last section presents a comparison with the work presented by Fang et al. (2015), states the main conclusions and provides directions for the future.
2. Materials and methods
2.1. Continuous self-propelled 1D crawlers
Model description and kinematics
Consider a 1D crawler moving along a straight line and assume its reference configuration is the segment
cf. Figure 1. We use the same formalism introduced by DeSimone and Tatone (2012) and DeSimone et al. (2013) so that X is the coordinate in the reference configuration while x(t) denotes the coordinate along the crawler's body in the current configuration (at time t). In particular x(t) is the image of X through the current transformation χ(·, t) which can be written in terms of the current distance s(·, t) from the left end, i.e.,
where x1(t): = χ(X1, t) and x2(t): = χ(X2, t). By definition, s(0, t)≡0 for all t and we assume that
in order to guarantee the monotonicity of χ(·, t) at any time t.
Figure 1
In what follows a prime will denote the derivative with respect to X while a superscript dot will denote the derivative with respect to t.
We define the displacement
so that
where we have implicitly defined the strain
in terms of which condition (1) reads
Finally, notice that, since the material (or Lagrangian) velocity is
the spatial (or Eulerian) velocity is given by
Equations of motion
Throughout this section we deal with the motility problem, namely, given a history of strain ϵ(X, t), the aim is to find x1(t) which determines the dynamics of the one-dimensional crawler.
Friction laws
The force at the interface between substrate and crawler is modelled through a force-velocity relationship. In particular, we write the density per unit current length of the tangential component of the friction force at time t, f(x, t), as a function of the Eulerian velocity v(x, t).
Several models for the resistance forces are conceivable such as a
Newtonian model, i.e., a linear viscous law
where μ>0 is a friction (or viscosity) coefficient;
or a more general “p-model”
for p∈[0, +∞). Parameter p in our force law (4) allows us to investigate different types of frictional behaviors. For p = 0, we obtain a force per unit current length that is a linear function of velocity alone, which reduces to the Newtonian model (3). For p>0, we obtain a friction law that is sensitive to the state of elongation of the segment, with force per unit length higher or lower than that of the Newtonian case depending on whether the element is contracted (λ < 1 or ϵ < 0) or extended (λ>1 or ϵ>0). In the limit p → ∞, this produces an idealized model for friction in which no force opposes slip when the segment is extended (free slip), while the segment can withstand any tangential force without sliding (perfect grip) when it is contracted. We call this idealized model “free slip–perfect grip.” Figure 2 displays the graphs of gp(ϵ) around ϵ = 0 for different values of p. In fact, our model is a continuous analog of the discrete model proposed by Fang et al. (2015) to mimic the behavior of earthworms' setae, which protrude when the body is axially contracted, resulting in an increment of the resistance (Edwards et al., 2018).
Figure 2
In what follows, we will use the p-model (4).
Force balance
The total friction is obtained by integrating the force per unit current length on the whole current domain, i.e.,
where
Then, by neglecting inertia, the force balance yields
The square bracket multiplying ẋ1(t) in the formula above is the drag for rigid motion at unit speed and fixed shape ϵ(X, t), while Fe(t) is an external force which, for instance, can take into account the gravity force acting on a crawler on an inclined plane. Solving for ẋ1(t), we obtain
which, in the case of zero external forces, is independent of μ.
Notice that, once the initial position x1(0) the strain ϵ(X, t) and the external force Fe(t) are provided, the whole dynamics x1(t) can be determined by integrating (5). Indeed, since
one has
and hence the right hand side of (5) is known once ϵ(X, t) is specified.
2.2. Discrete self-propelled 1D crawlers
We now move to a discrete model, directly inspired by studies on annelid worms (Quillin, 1999) and metameric robots (Menciassi et al., 2006; Daltorio et al., 2013; Fang et al., 2015).
Model description and kinematics
We model the crawler's body as made up of N segments of same length L in the reference configuration (cf. Figure 3)
for n = 1, …, N. Let x0(t) denote the current position of the left edge; any point X in the reference domain is mapped to a point x(t) in the current domain through the map
whence the definition of strain
Figure 3
Furthermore, we assume that each segment can be contracted or expanded according to a constant stretch so that the overall strain results to be a piecewise constant function of X (at any fixed time t), i.e.,
Consequently its time-derivative, , is piecewise constant and hence, from (7), ṡ(·, t) is piecewise affine, i.e.,
for X∈[Xn−1, Xn]. Note that in this new framework the monotonicity condition (2) reads
which is the only constraint for an admissible history of strains, the datum of our motility problem.
Equations of motion
Analogously to the previous case, the force balance yields
where, in view of (9),
and, in view of (10),
Solving for ẋ0(t) we obtain
which can be rewritten in the following vectorial form
where
Equation (11) fully describes the dynamics once x0(0) and ϵ(X, t) are provided. In particular, the displacement after T time units is given by
Relative displacement
We can rewrite everything in terms of a displacement relative to x0(t), defined as
in order to describe the displacement in a coordinate system which is “co-moving” with the left end x0(t).
In the discrete framework, the relative displacement turns out to be a piecewise affine function of X (at any fixed time t): if X∈[Xn−1, Xn],
Setting
we have
where u0≡0. Equivalently
where , and
Optimal control problems
In this section we address the problem of maximizing the net displacement Δx0 among periodic shape changes ϵ(X, t) with the same given energy cost.
We now describe the optimization problems with quadratic energy in the non-linear case first, and then in the small-deformation regime, for which general results can be established. We assume Fe≡0.
Feasible region
We assume that the shape function ϵ(t) is a C2 function defined from ℝ to ℝN.
In addition, we require ϵ(·) to be a time-periodic function. Finally we restrict our search to shape functions with a given cost per period, i.e., , where the energy functional is assumed to be of the following quadratic form (in both ϵ and )
where 𝔸 and 𝔹 are symmetric and positive definite N-dimensional matrices. Overall, the feasible region is
Optimization problem
The general (non-linear) optimization problem is
which is an isoperimetric problem (e.g., Van Brunt, 2004) involving N dependent variables ϵn. The corresponding Euler-Lagrange equations lead to a second order non-linear system of ODEs, i.e., for n = 1, …, N,
where (here λ denotes a Lagrange multiplier).
The small-deformation regime
We can focus on the small-deformation regime by expanding the objective function at the leading orders (about ϵ = 0), i.e.,
and, integrating by parts,
whence
In particular, it can be proved that the (skew-symmetric Toeplitz) matrix skw(vϵ(0)) =:V depends only on N, L and p. Indeed,
(see Appendix A in Supplementary Material).
Therefore, in the regime of small deformations, problem (18) can be replaced by the following linear problem
The corresponding Euler-Lagrange equations
where
lead to the following system of second order linear ODEs
In general, a solution to (21) might be difficult to determine due to the complexity of finding a common diagonalization of 𝔸 and 𝔹. However, following the procedure adopted by Wiezel et al. (2018), we can solve this problem when one of the two operators is null, say 𝔸≡0 (resp. 𝔹≡0), and the other one, 𝔹 (resp. 𝔸), is symmetric, positive definite and such that the eigenspaces associated with the maximum-modulus eigenvalues of (resp. ) have dimension 1. Indeed, as shown in sections 1 and 2 in Appendix B (Supplementary Material), it turns out that
for A = 0 and B symmetric and positive definite, up to a constant, a solution of (20) must be of the form
where α∈ℂ\{0} is a constant such that and is a suitable constant vector depending only on A and V.
for 𝔸 symmetric and positive definite and 𝔹 = 0, a solution of (20) with ϵ of unitary time frequency must be of the form
where α∈ℂ\{0} is a constant such that and is a suitable constant vector depending only on 𝔹 and 𝕍.
Both (22) and (23) have the form
i.e., they are circles in the plane (ℜ(e), ℑ(e)), regardless the number of links. Moreover, using the polar representations
we get, for n = 1, …, N,
i.e., the optimal gait depends only on the 2N+2 parameters {ϑn}n, {ϱn}n, ϑa and ϱa. Admittedly, since α is a constant with fixed modulus and free argument, we can always assume , i.e.,
thus reducing the number of parameters to 2N+1.
The problem for 𝔸 = 0 and 𝔹 = 𝕀N is essentially equivalent to the one for 𝔸 = 𝕀N and 𝔹 = 0 (provided that unitary time frequency of ϵ is prescribed). Indeed, if (25) is a solution to
then it is a solution also to
and vice versa.
In general the two problems, 𝔸 = 0 with 𝔹 symmetric positive definite and 𝔹 = 0 with 𝔸 symmetric positive definite, are not equivalent. In fact, constraining the norm induced by one operator does not determine the norm induced by the other one, but only provides a bound. Indeed, denoting by λmin(·) and λmax(·) the minimum and maximum eigenvalue respectively, observe that, for ϵ(t) like (25),
and, analogously,
3. Results
3.1. Continuous model: two examples of contraction waves
In this section we discuss two examples of contraction waves to illustrate the behavior of the p-model. We show that the parameter p determines the kind of motion: for p < 1 the motion is prograde (i.e., motion in the same direction as the one of the waves) while for p>1 the model reproduces an earthworm-like retrograde motion (i.e., motion in the opposite direction as the one of the waves).
For simplicity, in the following examples we neglect external forces, i.e., Fe≡0.
Smooth contraction wave
Consider a smooth traveling contraction wave by prescribing the strain along the body of the crawler as
or equivalently, in terms of the stretch,
where ϵ0 is the wave amplitude, L is the reference length of the crawler and c is a parameter which modulates time frequency and it is assumed to be strictly positive, i.e., the wave travels toward the right. By integrating over space,
and, by differentiating with respect to time,
Finally, in view of (5),
whence
The Newtonian case is recovered by setting p = 0, i.e.,
Figure 4 displays three numerical examples. For p < 1 and, in particular for p = 0, the case of Newtonian resistance, we always have prograde motion (i.e., motion in the same direction as the one of the waves). This is indeed observed for example in snails, although in this case the force-velocity laws that we use in this paper would not be fully adequate to capture the properties of the mucus present between the animal and the surface [non-Newtonian rheology, suction effects, see (Denny, 1980) and (DeSimone et al., 2013)]. For p>1 and, in particular, for the limit case p = ∞ describing the perfect-grip/free-slip ideal version of the modulated friction laws typical of animals with setae, the motion is retrograde (i.e., motion in the opposite direction as the one of the waves). This is the behavior typically observed for earthworms.
Figure 4
Square contraction wave
Consider the square contraction wave
or equivalently, in terms of the stretch,
where L is the reference length of the crawler, c is the wave speed, ξ is the measure of the interval where ϵ = δ and the subscript ~L denotes the “modulo L” operator (i.e., y~L stands for ymodL).
By integrating the stretch over space, we get
whence
Finally, in view of (5),
where
and
Defining and , it follows that
where {·} and ⌊·⌋ denote the fractional part and floor function respectively. The Newtonian case can be obtained as particular case by setting p = 0. Figure 5 shows three numerical examples and, as for the smooth contraction wave, the motion is prograde or retrograde whether p < 1 or p>1, respectively.
Figure 5
3.2. Discrete model: peristalsis as optimal gait
In the discrete framework, peristalsis is the result of phase coordination among the harmonic contractions of body segments, i.e., it has the form
where T is the period, ϱ is the amplitude and Δφ is the constant phase difference. As for the continuous case, discrete peristalsis produces prograde or retrograde motions according to the value of the parameter p in (4).
In this section we work out explicitly the problem of maximizing the displacement for a particular case from which peristalsis emerges, modulo an edge-effect.
Dissipation energy
Let us define an energy functional D:C2(ℝ, ℝN) → ℝ as
where
i.e., the energy cost is the time integral over a period of a dissipation rate which is sum of two terms: is times the energy expended to overcome the friction force and is the cost of control weighted by a scalar factor w. D[ϵ] is thus times the sum of the work due to the friction force plus the L2-norm of the controls suitably weighted to time the input direction.
Some calculations [see section 1 in Appendix C (Supplementary Material)] lead to
where 𝔻(ϵ)∈ℝN×N for any ϵ∈(−1, +∞]N, and
where 𝕀N is the N-dimensional identity matrix. Therefore the energy functional is
where 𝔾(ϵ): = 𝔻(ϵ)+wIN.
Non-linear optimal control problem
The non-linear optimization problem associated with energy functional (28) is
The Euler-Lagrange equations lead to a second order non-linear system of ODEs, i.e., for n = 1, …, N,
where , λ is the Lagrange multiplier.
The small-deformation regime
In the regime of small deformations we can expand the terms of problem (29) at the leading orders about ϵ = 0. As before, the net displacement per time period can be approximated by
and the energy functional by
where 𝔾: = 𝔾(0). Hence, in the small-deformation regime, the problem fits the form (17)-(18) for 𝔸 = 0 and 𝔹 = 𝔾. Moreover, G is bisymmetric (namely, symmetric about both of its diagonals) and depends only on N, L and w. Indeed
[see section 2 in Appendix C (Supplementary Material)]. Therefore a solution must be of the form (25)
The centrosymmetry of 𝔾 and the skew-centrosymmetry of 𝕍 imply a reflectional symmetry about the center [see section 3 in Appendix B (Supplementary Material)], i.e.,
the moduli of components of e are symmetric about the center (cf. Figure 6), i.e.,
phase differences between adjacent segments are symmetric about the center, i.e.,
so that the N-th phase differs from the (N−1)-th one by the same amount by which the second phase differs from the first one and so on (cf. Figure 6).
Figure 6
Equation (25) shows that the optimal gait requires a precise “phase coordination” of locomotion patterns among the segments, which is a common observation in Biology for several kinds of animals.
Numerical simulations show that the optimal solution turns out to be a discrete approximation of a traveling wave. In particular,
the moduli of en for n = 1, …, N can be approximated by a constant average value (cf. Figure 6), i.e.,
so that each segment undergoes a harmonic deformation with a certain initial phase;
phase differences between adjacent segments turn out to be almost constant, i.e., for a suitable φ0
holds true for n = 1, …, N (cf. Figure 6).
Therefore, in view of (32) and (33), the solution is a discrete approximation of a continuous traveling wave, i.e.,
where is the midpoint of the n-th segment and
whence the canonical form of traveling wave which describes peristalsis (cf. Figure 7)
where and .
Figure 7
The edge-effect
The symmetric structure of the optimal gait (in the small-deformation regime) arises from underlying physical symmetries which clearly stand out in the properties of the matrices 𝔾 and 𝕍. In particular, an “edge-effect” is apparent: the 1D crawler is symmetric about its geometric center and segments near the edges behave differently with respect to adjacent segments, but in the same way as their centrosymmetric counterparts.
As expected, this edge-effect vanishes when considering an “infinite” (periodic) 1D crawler because, due to the shift-invariance symmetry, each segment behaves as a “geometric center.”
To show this claim consider a 1D crawler made up of infinitely many segments and assume that it is a periodic structure of which each module consists of N components (cf. Figure 8).
Figure 8
At any time t, we already defined the relative displacement u(·, t) as the change of position of the material point X in the body's reference, i.e.,
The hypothesis of periodicity leads to
From (35) we obtain that the friction force is periodic and we can consider the force balance in a single module. Therefore, condition (35) reads
and, in the discrete framework (9), this leads to
The optimal control problem becomes
and it can be proved [see section 3 in Appendix C (Supplementary Material)] that its solutions need to be like (22), where the complex N-dimensional vector e has the form
for some k∈{1, …, N−1} and e1∈ℂ\{0}, cf. Figure 9. In particular,
- each component of e has modulus ϱ: = ||e1||;
- each component can be obtained from the previous one by a rotation of or, in other words, the phase difference between two consecutive components is constant, i.e., for n = 1, …, N
where and ;
Figure 9
whence the exact harmonic peristalsis.
Notice that problem (38) can be written in terms of relative displacements un through the periodic version of transformation (15), i.e.,
where
In particular, we get
where
are circulant matrices (namely, Toeplitz matrices where each row vector is rotated one element to the right relative to the preceding row vector), thus reflecting the geometric symmetry of the periodic structure, namely, the shift-invariance.
Considering the general (i.e., non-periodic) problem in terms of relative displacements yields
where
are “quasi-circulant” matrices indeed
where 𝔼V and 𝔼G are null a part from the last column and the last row, i.e.,
and
where a = 2L3(3−N)+12Nw and b = L3(9−4N)−12Nw.
Wavenumber
In the “periodic case” we can study the wavenumber (that is the number of waves travelling along the body of the crawler) of the optimal gait in relation to the number of metameres N and to the weight w. We fix the dissipation
and we let w vary from 0 to 102 for N∈[3, 250] (see Figure 10). As shown in Section 3 in Appendix C (Supplementary Material), the wavenumber of the optimal gait must be an integer close to the real number and hence, for any fixed N, it depends on the order of magnitude of the weight w and
for w → ∞, it tends to 1, corresponding to a single wave spanning the whole length NL;
for w = 0, it is close to , i.e., one full wave-length every three segments.
Figure 10
This behavior is qualitatively unaffected by the type of friction model which is adopted (i.e., by the choice of the parameter p).
4. Discussion
4.1. Comparison with previous studies
To put our study in perspective, we consider the discrete framework and we compare our results with the ones presented by Fang et al. (2015). Here the authors perform an optimization of the so-called “average steady-state velocity” us among harmonic shape functions having the form (in our notation)
where is the oscillation amplitude, T is the period and ηn is the actuation phase for the n-th segment (or actuator).
Since the average steady-state velocity is given by
the optimization problem reads
and in the small-deformation regime it can be replaced by
Denote the actuation phase differences between adjacent segments by
From observations of numerical simulations, Fang et al. (2015) report that “[…] the optimized phase-different patterns are always reflectionally symmetric [about the center, Ed.] regardless of the initial symmetry requirements […]” and of the number of segments. Thus, a solution to (45) fulfills
In fact these properties can be rigorously proved under the assumption that problem (45) admits a unique solution in [0, 2π)N (see Appendix D in Supplementary Material).
Furthermore, property (47) can be proved also for (46), assuming it admits a unique solution in [0, 2π)N. To this aim, denote the unique solution to (46) by
and consider the shape change associated with
where
Notice that for n = 1, …, N,
and hence, by exploiting the fact that V is skew-centrosymmetric (namely, 𝕂T𝕍𝕂 = −𝕍),
Thus
which leads to (47).
Problem (46) constrains the L2-norm of the time-derivatives, i.e., for strains having the form (43) we get
regardless of η. Therefore we can extend the maximization to the C2 periodic strains whose time derivative fulfills the same constraint, i.e.,
Since problem (21) reduces to (48) when 𝔸 = 0 and 𝔹 = 𝕀N, a solution to (48) must be of the form
where e = (en)n is a unit eigenvector associated with the maximum-modulus eigenvalue of V and ϑa is a constant. Notice that the reflectional symmetry about the center still holds true. As a matter of fact, (49) leads to a slight increment in the net displacement with respect to the solution to (46), cf. Figure 11.
Figure 11
4.2. Summary and outlook
Our analysis confirms the effectiveness of mimicking peristalsis in bio-inspired robots, at least in the small-deformation regime. This bio-inspired actuation strategy has been implemented on a trial-and-error basis many times in the robotics literature and, more recently, also proposed as optimal (in some suitably defined sense, and in some suitably defined class of actuation strategies). Our main result is a mathematically rigorous proof that, in the small deformation regime, actuation by peristaltic waves is an optimal control strategy emerging naturally from the geometric symmetry of the system, namely, the invariance under shifts along the body axis. This is true exactly in the periodic case, and approximately true in the case of finite length, modulo edge-effects. Our results is of theoretical nature. Nevertheless, we believe that it has important consequences for applications. For example, it helps us not to fall into the naive temptation to expect that peristaltic waves are always an optimal actuation strategy just because they are observed in nature, but to exercise critical judgment whenever the hypotheses on the geometric symmetry that are “responsible” for the optimality result in our case (invariance under shift of a homogeneous one-dimensional system) are false.
Actuation by phase coordination, optimal actuation by identical phase difference, and the connections between this and traveling waves have been already discussed in the literature (e.g., Fang et al., 2015), but never through a mathematically rigorous analysis of the optimal control problem, of the symmetry properties of the governing equations and operators, and of the relation between these and the geometric symmetries of the system. This is exactly what we do in this paper. The added value of this analysis is that we are able to show (for the first time, to the best of our knowledge, at least in the robotics literature) that peristaltic waves are the signature of the invariance with respect to shifts (a geometric symmetry) of a homogeneous one-dimensional system.
Further work will be needed to test the effectiveness of peristaltic waves as a locomotion strategies if large deformations are allowed. In addition, future work will explore the issue of how peristalsis is actually enforced in biological systems. Of particular interest is the dichotomy between the paradigm of actuations via a Central Pattern Generator (CPG), as opposed to local sensory and feedback mechanisms. The CPG paradigm is apparent in several different organisms (Marder and Bucher, 2001; Grillner, 2006) and has been employed in robotics with some success (Ijspeert, 2008; Boxerbaum et al., 2012). However, there is a growing awareness of the role played by proprioception, especially for lower organisms such as the nematode worm C. elegans (Boyle et al., 2012; Wen et al., 2012) and D. melanogaster larvae (Pehlevan et al., 2016).
Statements
Author contributions
AD and FA conceived research. DA executed research and performed numerical simulation. AD supervised research. FA contributed expertise on circulant matrices. All authors analyzed the data and wrote the manuscript.
Funding
We gratefully acknowledge the financial support of the European Research Council through the Advanced Grant 340685-MicroMotility.
Acknowledgments
This article was written while DA and AD where visiting Division C of the Engineering Department of the University of Cambridge. We thank Prof. Vikram Deshpande and Corpus Christi College for their warm hospitality.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The handling Editor declared a shared affiliation, though no other collaboration, with one of the authors, AD.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frobt.2018.00099/full#supplementary-material
References
1
AlbertsB.JohnsonA.LewisJ.RaffM. (2002). Molecular Biology of the Cell, 4th Edn.New York, NY: Garland Science.
2
BoxerbaumA.ShawK.ChielH.QuinnR. (2012). Continuous wave peristaltic motion in a robot. Int. J. Robot. Res.31, 302–318. 10.1177/0278364911432486
3
BoyleJ. H.BerriS.CohenN. (2012). Gait modulation in c. elegans: an integrated neuromechanical model. Front. Comput. Neurosci.6:10. 10.3389/fncom.2012.00010
4
CollarA. (1962). On centrosymmetric and centroskew matrices. Q. J. Mechan. Appl. Math.15, 265–281. 10.1093/qjmam/15.3.265
5
DaltorioK.BoxerbaumA.HorchlerA. S.ShawK. M.ChielH. J.QuinnR. D. (2013). Efficient worm-like locomotion: slip and control of soft-bodied peristaltic robots. Bioinspir. Biomimet.8:035003. 10.1088/1748-3182/8/3/035003
6
DennyM. (1980). The role of gastropod pedal mucus in locomotion. Nature285:160. 10.1038/285160a0
7
DeSimoneA.GuarnieriF.NoselliG.TatoneA. (2013). Crawlers in viscous environments: linear vs non-linear rheology. Int. J. Non Linear Mech.56, 142–147. 10.1016/j.ijnonlinmec.2013.02.007
8
DeSimoneA.TatoneA. (2012). Crawling motility through the analysis of model locomotors: two case studies. Eur. Phys. J. E Soft Matt. Biol. Phys.35, 1–8. 10.1140/epje/i2012-12085-x
9
EdwardsC.HendrixP.AranconN. (2018). Biology and Ecology of Earthworms. Springer. Available online at: https://www.springer.com/in/book/9780387749426
10
FangH.LiS.WangK. W.XuJ. (2015). Phase coordination and phase–velocity relationship in metameric robot locomotion. Bioinspir. Biomimet.10:066006. 10.1088/1748-3190/10/6/066006
11
GardnerC. (1976). The neuronal control of locomotion in the earthworm. Biol. Rev.51, 25–52. 10.1111/j.1469-185X.1976.tb01119.x
12
GarreyW.MooreA. (1915). Peristalsis and coordination in the earthworm. Am. J. Physiol. Legacy Cont.39, 139–148. 10.1152/ajplegacy.1915.39.2.139
13
GeJ.CalderónA.Pérez-ArancibiaN. (2017). An earthworm-inspired soft crawling robot controlled by friction. arXiv [preprint] arXiv:1707.04084.
14
GrayJ.LissmannH. (1938). Studies in animal locomotion: Vii. locomotory reflexes in the earthworm. J. Exp. Biol.15, 506–517.
15
GrillnerS. (2006). Biological pattern generation: the cellular and computational logic of networks in motion. Neuron52, 751–766. 10.1016/j.neuron.2006.11.008
16
IjspeertA. (2008). Central pattern generators for locomotion control in animals and robots: a review. Neural Netw.21, 642–653. 10.1016/j.neunet.2008.03.014
17
MarderE.BucherD. (2001). Central pattern generators and the control of rhythmic movements. Curr. Biol.11, R986–R996. 10.1016/S0960-9822(01)00581-4
18
McInerneyA. (2013). First Steps in Differential Geometry: Riemannian, Contact, Symplectic. Undergraduate Texts in Mathematics. New York, NY: Springer.
19
MenciassiA.AccotoD.GoriniS.DarioP. (2006). Development of a biomimetic miniature robotic crawler. Auton. Robots21, 155–163. 10.1007/s10514-006-7846-9
20
NemitzM. P.MihaylovP.BarracloughT. W.RossD.StokesA. A. (2016). Using voice coils to actuate modular soft robots: wormbot, an example. Soft Robot.3, 198–204. 10.1089/soro.2016.0009
21
PehlevanC.PaolettiP.MahadevanL. (2016). Integrative neuromechanics of crawling in D. melanogaster larvae. Elife5:23. 10.7554/eLife.11031
22
QuillinK. (1999). Kinematic scaling of locomotion by hydrostatic animals: ontogeny of peristaltic crawling by the earthworm lumbricus terrestris. J. Exp. Biol.202, 661–674.
23
UmedachiT.KanoT.IshiguroA.TrimmerB. (2016). Gait control in a soft robot by sensing interactions with the environment using self-deformation. R. Soc. Open Sci.3:160766. 10.1098/rsos.160766
24
Van BruntB. (2004). The Calculus of Variations. New York, NY: Springer.
25
WangK.YanG.MaG.YeD. (2009). An earthworm-like robotic endoscope system for human intestine: design, analysis, and experiment. Anna. Biomed. Eng.37, 210–221. 10.1007/s10439-008-9597-6
26
WenQ.PoM. D.HulmeE.ChenS.LiuX.KwokS. W.et al. (2012). Proprioceptive coupling within motor neurons drives c. elegans forward locomotion. Neuron76, 750–761. 10.1016/j.neuron.2012.08.039
27
WiezelO.GiraldiL.DeSimoneA.OrY.AlougesF. (2018). Energy-optimal small-amplitude strokes for multi-link microswimmers: Purcell's loops and taylor's waves reconciled. arXiv [preprint] arXiv:1801.04687.
Summary
Keywords
crawling motility, lumbricus terrestris, peristalsis, self-propulsion, metameric robots, biomimetic robots, soft robotics, optimization
Citation
Agostinelli D, Alouges F and DeSimone A (2018) Peristaltic Waves as Optimal Gaits in Metameric Bio-Inspired Robots. Front. Robot. AI 5:99. doi: 10.3389/frobt.2018.00099
Received
15 April 2018
Accepted
31 July 2018
Published
05 September 2018
Volume
5 - 2018
Edited by
Matteo Cianchetti, Scuola Sant'Anna di Studi Avanzati, Italy
Reviewed by
Shinya Aoi, Kyoto University, Japan; Barbara Mazzolai, Fondazione Istituto Italiano di Technologia, Italy
Updates
Copyright
© 2018 Agostinelli, Alouges and DeSimone.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Antonio DeSimone desimone@sissa.it
This article was submitted to Soft Robotics, a section of the journal Frontiers in Robotics and AI
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.