Studying slippage on pushing applications with snake robots
Abstract
In this paper, a framework for analyzing the motion resulting from the interaction between a snake robot and an object is shown. Metrics are derived to study the motion of the object and robot, showing that the addition of passive wheels to the snake robot helps to minimize slippage. However, the passive wheels do not have a significant impact on the force exerted onto the object. This puts snake robots in a similar framework as robotic arms, while considering special properties exclusive to snake robots (e.g., lack of a fixedbase, interaction with the environment through friction). It is also shown that the configuration (shape) of the snake robot, parameterized with the polar coordinates of the robot’s COM, plays an important role in the interaction with the object. Two examples, a snake robot with two joints and another with three joints, are studied to show the applicability of the model.
Keywords
Background
Robots that are capable of locomotion in unstructured conditions are necessary for realistic applications. However, locomotion alone may not be sufficient when more dexterous interaction with the environment is needed. Therefore, robotic systems with capability to locomote and also interact dexterously with their surroundings are desirable, and indeed a natural extension of robotics research.
Snake robots have shown promise regarding locomotion [1]. Locomotion in planar environments has been probably the main topic of research for snake robots [2, 3, 4] and has been extended to motion in planar slopes [5, 6], motion in 3Dspace [7, 8], and more broad studies on locomotion [9]. An interesting idea that combines locomotion and interaction with the environment, called obstacleaided locomotion (OAL), has been proposed in [10] where obstacles in the environment are used as auxiliary sources for propulsion or to avoid jamming.
Although snake robots could excel in locomotion, it is not clear if they can be used to interact with the environment (or an object) dexterously. Its structure resembles a robotic manipulator, but there are key differences that have not been fully addressed in previous research (c.f. Fig. 1).
The lack of a fixedbase makes it difficult for a snake robot to manipulate an object as dexterously as a robotic arm. Another difference is that a snake robot has contact with the environment through friction at several points of its body. Additionally, mass becomes a very important parameter to study. Unlike research regarding robotic arms where it is assumed that the arm can lift the object and it is a matter of choosing an optimal input, snake robots may not be able to move the object due to its inertial properties.
Because the kinematic structure of a snake robot resembles a robotic arm, papers that deal with similar (but not exactly the same) situations can be found in existing literature. In [11] a hyperredundant serial robot was considered and both locomotion and manipulation of an object were considered. However, the analysis was purely kinematic while assuming a fixedbase robotic system. In other words, there was not force analysis showing the conditions for feasibility of the problem. In [12], the duality between locomotion and manipulation of a snake robot was considered under the assumption that the snake robot can be treated similarly to a robotic arm with a fixedbase when manipulating an object. This was achieved by making the first link of the snake robot behave similarly to fixedbase (due to its shape and mass), but the results cannot be extended to the case of a general snake robot. The problem of analyzing and controlling a snake robot under these conditions has been reported in [13], where it is shown that several assumptions made in previous published literature are not enough to guarantee accurate control of a planar snake robot with frictional contacts with the ground.
The main objective of this paper is to study the resulting interaction between a snake robot and an object, when the task is to push the object. We consider this to be a prelude to more interesting interactions like grasping or dexterous manipulation. However, it is important to understand the basics of the interaction first. This is an extension of previously published work [14, 15] where a more complete mathematical modeling of the problem has been presented. In [15], the optimal configurations of the snake robot to maximize the force exerted onto the object have been presented.
This paper focuses on the motion of the system, rather than forces. The main motivation for this study is that, as presented in [13], calculating an optimal input instantaneously (i.e., at one instant of time) is not enough to accurately control the system. We conjecture that understanding how the system will behave over time is also important. In other words, if the task is to push an object, then the motion of the object must be maximized, while the motion of the snake robot minimized.
Throughout this paper, there are several assumptions that have to be made because of the complexity of the problem:
 1.
All bodies in the system are rigid.
 2.
Contacts between the snake robot and objects or the environment (except the ground) are considered frictionless point contacts.
 3.
The snake robot has passive wheels or any other mechanical means to achieve anisotropic friction between the robot’s belly and the ground.
 4.
There is only one constraint per link of the snake robot.
 5.
The operational space is a plane (embedded in a full 3D space).
 6.
The snake robot has only one contact point with the object to be manipulated.
Assumption 1 allows to have a clear mathematical model of the problem without making assumptions about the compliance of the bodies which may be unrealistic to know a priori in reallife situations. Although this assumption may be relaxed by considering some sort of virtualcompliance at the contacts (e.g., [16, 17]), it does not necessarily imply more realistic or correct results. In particular, it may lead to a stiff system of differential equations and several other problems. As presented in [14], a snake robot may have too many contacts with the environment leading to a statically indeterminate system [18]. In order to ensure Assumption 4, we assume that the constraints from passive wheels are removed when that link is contacting an object or a wall, for example. This can be done by lifting the links [19] or with retractable passive wheels, for example. Assumptions 1 through 4 allow to consider that all constraint forces are linearly independent and a unique solution can be found. Assumption 5 is made in order to limit the number of parameters and get clear and meaningful results which may be difficult for 3dimensional space, since the problem presented in this paper is still very broad and has not been clearly defined as it has been discussed in this section. However, the models presented in this paper are based on spatial vectors [20] which are trivial to extend from 2D to 3D, so in the future more results can be obtained for more specific tasks. Assumption 6 is made because it is not the intent of this paper to study any type of grasp closure or dexterous manipulation, but to understand the interaction itself first.
The paper is organized as follows. In “Mathematical background” section, the necessary mathematical background to understand this paper is presented along references necessary to develop the concepts further. “Motion of the system due to the interaction” section is the main body of the paper; the modeling of the system is presented, and metrics and quantities mentioned in this section are derived. In “Results” section, a specific example is studied to show the application of the proposed metrics. “Discussion and future applications” section includes several comments regarding the scope and limitations of the results presented in this paper. The paper concludes with some remarks in “Conclusion” section.
Mathematical background
In this section, we give a very brief introduction to the mathematical topics necessary to understand this paper. We recommend [20, 21, 22] for a more detailed treatment. In particular, the foundations of the model used in this paper have been presented in [15]; readers are encouraged to read this reference for a more detailed treatment of snake robots in the framework of articulatedbodies. As stressed in previous research, it is important to guarantee invariance of metrics in order for the results to be meaningful. Not only to avoid inconsistency of units, but also for the metric to be invariant to a change of coordinates. To derivate the analysis of the system to lead to meaningful metrics, we employ dual vectors [23] and basic differential geometry [21].
Differential geometry: twists, wrenches, and metrics
A twist (concatenation of linear and angular velocity) υυ→∈Mnυυ→∈Mn can be expressed w.r.t. a covariant basis ee=[ee→1,…,ee→n]Tee=[ee→1,…,ee→n]T as υυ→=eeTυυυυ→=eeTυυ. The element υυ∈Rnυυ∈Rn can be interpreted as the (vector of) contravariant components of υυ→υυ→. A wrench (concatenation of linear force and torque) ff→∈ F nff→∈ F n can be expressed w.r.t. a contravariant basis ee∗=[ee→∗1,…,ee→∗n]Tee∗=[ee→1∗,…,ee→n∗]T as ff→=ee∗Tff∗ff→=ee∗Tff∗and ff∗∈Rnff∗∈Rn can be interpreted as the covariant components of ff→ff→.
It is important to notice that both bases eeee and ee∗ee∗ may not be orthogonal, so the common definition of inner product (e.g., υυ⋅υυ=υυTυυυυ⋅υυ=υυTυυ) would give incorrect results. Let us denote the metric tensor of the covariant basis as II=eeeeTII=eeeeT and its inverse by II−1II−1. The (squared) length of a twist and wrench is an invariant quantity and can be obtained using the scalar product {∘}{∘} while taking into account the metric tensor as
Unconstrained model of the snake robot
The snake robot can be modeled as a series of rigid links connected by revolute joints. All joints have their axes parallel to each other; therefore, the snake robot is constrained to move on a plane (but is unconstrained in any other way). The kinematic model of a snake robot is similar to an openchain robotic manipulator (c.f. Fig. 2). The model has been previously studied in [13, 15, 24, 25].
The snake robot has a total of ns∈Nns∈N degrees of freedom (DOFs), and its generalized coordinates are encapsulated in the vector qqs(t)∈Rnsqqs(t)∈Rns. The snake robot has nℓ=na+1nℓ=na+1 links each with mass mimi.
The Jacobian for the ith link is a mapping from the vector of generalized velocities qq˙sqq˙s to the twist υυi∈R3υυi∈R3 of the link and is denoted as JJi∈R3×nsJJi∈R3×ns
The equations of motion of the snake robot can be presented in the canonical form
where MMs(qqs)∈Rns×nsMMs(qqs)∈Rns×ns is the inertia matrix of the snake robot (a symmetric positive definite (PD) matrix), hh∗s(qqs,qq˙s)∈Rns×1hhs∗(qqs,qq˙s)∈Rns×1 contains Coriolis and centripetal effects, and ττ∗ext(qq,qq˙)∈Rns×1ττext∗(qq,qq˙)∈Rns×1is a vector of torques produced by external forces (e.g., kinetic friction). The matrix BB∈Rns×naBB∈Rns×nadefined as
is a matrix that projects the vector of input forces ττ∗actττact∗ into the space of generalized forces. The matrix 1111 denotes the identity matrix of appropriate dimensions.
Unconstrained model of an object
A rigid body is able to move in its operational space with dimensions nop,nop, and its equations of motion can be compactly written as
where aaobjaaobj, pp∗objppobj∗, and ff∗obj∈Rnopffobj∗∈Rnop denote the acceleration, velocityproduced terms, and total wrench acting on the body, respectively. If the body is constrained to move in a plane (but unconstrained in any other way), it will have three DOFs (i.e., nop=3nop=3). IIobj∈ R nop×nopIIobj∈ R nop×nop denotes the inertia tensor of the rigid body. The mass of the object mobjmobj will be denoted as a multiple of the mass of a link of the snake robot as mobj=κmimobj=κmi. In other words, κκ is a proportionality coefficient relating the masses of interest.
If all links of the snake robot have the same mass m, then the inertia matrix of the snake robot can be factored as MMs:=mMM¯MMs:=mMM¯, and the inertia matrix of an object as IIobj=mobjII¯objIIobj=mobjII¯obj, where the new inertia matrices MM¯MM¯ and II¯objII¯obj correspond to inertia matrices with unitary mass.
Summary of constraints
The interaction between a snake robot and an external object creates a set of forces between them that avoid penetration (also called kinematic constraints or nonpenetrability constraints [21, 26, 27]). Additionally, the (static) friction forces between the belly of the robot and the ground can also be modeled as constraint forces (bounded by their friction limit). Assuming there are nc∈Nnc∈N constraint forces in total, the constraint forces ff∗c∈Rncffc∗∈Rnc span the constrained subspace
where the matrix TT∈nopnc×ncTT∈nopnc×nc is a matrix spanning the constraint forces on the operational space [21], and λλ∗∈Rncλλ∗∈Rnc contains the magnitude of the constraint forces (in the context of optimization this vector is usually called the Lagrangian multipliers [21, 26]).
To facilitate the coupling between the snake robot and environment/object(s), it is useful to put together all the constraints in vector/matrix form. All the constraints can be put together into the following form
where the AA∈Rns×nAA∈Rns×n is called the constraint matrix and takes the following form
where JJs∈Rnc×nsJJs∈Rnc×ns is called the robot Jacobian (also called hand Jacobian [26, 27]) which projects the vector of generalized velocities of the snake robot onto the constrained subspace
where TTkTTk spans the constrained space for the kth constraint and JJ∗JJ∗ denotes the Jacobian corresponding to the link under that constraint, without any specific ordering. The matrix GG∈Rnop×ncGG∈Rnop×nc is usually referred to as Grasp Matrix and its transpose is a mapping from the motion space of the object to the constrained subspace; it can be constructed in a similar manner to the robot Jacobian. The constraint forces projected back onto the snake robot and object are
where ττ∗c∈Rnsττc∗∈Rns is simply the projection of the constraint reaction force −ff∗c−ffc∗ onto the space of generalized forces of the snake robot.
Motion of the system due to the interaction
Now that it is assumed that the snake robot is touching at least one object, the new coupled equations of motion can be written as
along the constraints
where equality holds for constraints imposed by friction. Equation (13) is the derivative of (8). As discussed in [21], both holonomic and nonholonomic constraints can be obtained in this uniform manner at the acceleration level.
The change from the contravariant basis ee∗ee∗ to the basis of the constrained space ee∗ceec∗ can be obtained by using the projector ee∗cΦΦee∗:Rn→Rnceec∗ΦΦee∗:Rn→Rnc [15, 21] defined as
The positive semidefinite (PSD) matrix GGc:=(AAII−1AAT)GGc:=(AAII−1AAT) represents the metric tensor for the basis ee∗ceec∗ and has full rank if all constraints are linearly independent. (Additional comments regarding the rank of this metric tensor are located in “Discussion and future applications” section). The mapping (14) can be interpreted as the left pseudo inverse of the matrix AATAAT as
The constraint forces λλ→=ee∗Tcλλ∗λλ→=eec∗Tλλ∗ can be obtained by projecting the equations of motion (12) onto the constrained subspaces using the projector (14) and taking into account the constraint (13) as
The righthand side (RHS) of (16) has two terms. The first term depends purely on the set of forces exerted onto the system, either by the actuators of the snake robot or an external wrench exerted onto the object. The second term includes terms produced by velocity and will vanish if the system starts from an equilibrium configuration. This (affine) system of equations is usually interpreted as a force ellipsoid [27, 28, 29], and it maps a quadratic region in the input space ff∗ff∗to an ellipsoid in the output space of constraint forces λλ∗λλ∗, while the velocityproduced terms will shift the origin of such ellipsoid. In this paper, it is assumed the system starts from equilibrium so that the following linear mapping can be defined
The obtained contrained forces (due to the inputs in the system) can be substituted back onto the equations of motion (12). The resulting acceleration of the system is
where the new projector ee∗⊥cΦΦee∗:Rn→Rneec∗⊥ΦΦee∗:Rn→Rn defined as
is a projector from input wrenches ee∗Tff∗ee∗Tff∗ to the space orthogonal to the constrained space CC, but with coordinates expressed with respect to the original basis ee∗ee∗. The extra inverse inertia II−1II−1transforms to coordinates in eeee basis (i.e., transforms from wrenches to twists).
By solving for qq¨sqq¨s and aaobjaaobj, the motion of the snake robot and object can be obtained as
where the auxiliary mappings qq¨sΦΦττ∗act:Rna→Rnsqq¨sΦΦττact∗:Rna→Rns, qq¨sΦΦff∗obj:Rnop→Rnsqq¨sΦΦffobj∗:Rnop→Rns,
aaobjΦΦττ∗act:Rna→RnopaaobjΦΦττact∗:Rna→Rnop, and aaobjΦΦff∗obj:Rnop→RnopaaobjΦΦffobj∗:Rnop→Rnop can be defined as
The (squared) length of the accelerations of the snake robot qq¨→s2qq¨→s2 or object aa→obj2aa→obj2 can be obtained in an invariant way by considering its metric tensors, where the total expression can be divided into three terms as
We will concentrate on the contributions of the inputs of the snake robot. It can be verified that the auxiliary mappings ΞΞττ∗actΞΞττact∗ and ΩΩττ∗actΩΩττact∗, after some manipulation, can be defined as
where the auxiliary term
has been introduced for a more compact notation and any further simplification has been omitted for simplicity’s sake. However, the linear relationship w.r.t. the masses of the system becomes evident.
Slippage ratio
As stated in “Background” section, it is an important problem to predict motion and not only forces, in order to understand the interaction between the snake robot and object and try to accomplish a task. If the task is to manipulate an object, then it is desirable to maximize the motion of the object aa→obj2aa→obj2 while minimizing the slippage of the snake robot qq¨→s2qq¨→s2. On the other hand, a snake robot could locomote using the environment as a source of propulsive forces or as a support, similar to the idea of climbing [30, 31, 32]. This case resembles more a walking robot where the contact with the environment is necessary for the robot to move. To the best of our knowledge, this distinction has not been studied with snake robots. To analyze this, we propose the ratio of accelerations
and call it slippage ratio which is a dimensionless scalar quantity bounded as sr∈[0,1]sr∈[0,1]. Using this ratio, we can analyze the following three general situations:

sr→1sr→1 which implies that the acceleration of the snake robot is minimal (qq¨→s2≪aa→obj2qq¨→s2≪aa→obj2 or qq¨→s2≈0qq¨→s2≈0).

sr≈0.5sr≈0.5 which implies a similar magnitude of acceleration for the two subsystems (qq¨→s2≈aa→obj2qq¨→s2≈aa→obj2).

sr→0sr→0 which implies that the magnitude of the acceleration of the object is minimal (qq¨→s2≫aa→obj2qq¨→s2≫aa→obj2 or aa→obj2≈0aa→obj2≈0).
This quantity can be seen as the ratio between a desired output and the total output. By analyzing the slippage ratio sr, given a configuration and input, we can understand better the behavior of the system.
Polar coordinates of the COM of the snake robot
In order to compare snake robots with different number of joints, it is necessary to parameterize the configuration of the robot with a set of parameters in common. The idea of using the polar coordinates of the COM of the snake robot w.r.t. the contact point with the object has been introduced in [15] and further explored in [24]. The (unsigned) distance from the COM of the snake robot to the contact point is denoted by COM, and the angle between this vector and the link contacting the object is denoted as ∠COM∠COM. These quantities can be seen in Fig. 3.
Results
To study the interaction between snake robot and object, we can apply the framework proposed in this paper while changing the number and type of constraints and studying the resulting accelerations of the system. In general, we propose three different scenarios depending on the type of constraints present on the system as follows:

Scenario 1: The snake robot is in contact with an object but unconstrained in any other way. The friction between the snake robot and ground is negligible.

Scenario 2: The snake robot contacts one object and has passive wheels in all other links. The friction between the passive wheels and ground is bounded by its limit surface.

Scenario 3: The snake robot is contacting one object and has passive wheels in all other links. The passive wheels impose (unbounded and bilateral) nonholonomic constraints.
In other words, we will change the properties of the interaction between the snake robot and the environment (ground) and then analyze the resulting acceleration of the object aa→obj2aa→obj2, of the snake robot qq¨→s2qq¨→s2 and the slippage ratio sr as a result. Scenario 1 allows us to consider only the inertial properties of the system. Scenario 2, on the other hand, allows us to study the effect that passive wheels have on the system, but with a bounded friction coefficient μsμs. Scenario 3 considers ideal passive wheels and could be considered as the extreme case when μs→∞μs→∞. This is the most common model used for studying locomotion of snake robots.
The norms (26), (27) and slippage ratio (30) have to be studied over regions of the input space on a specific kth configuration. The inputs are restricted to the quadratic region of the input
which will be, in general, an ellipsoid and not a unitary sphere as is usually considered (i.e., ττ→act2⩽1ττ→act2⩽1 is not the same as ττ∗Tactττ∗act⩽1ττact∗Tττact∗⩽1).
Case study 1: Snake robot with two joints
In order to show more specific and qualitative results, we apply the mappings and study a snake robot with two joints (c.f. Fig. 3). The small number of joints allows us to show graphically the magnitude of the studied norms as a function of the joint torques. First, we study a snake robot with the parameters described in Table 1. The snake robot is contacting an object with its tail (first link) and the contact occurs at the middle of the link (c.f. Fig. 3). The angles of the joints are varied in the range [− 135°, 135°] every 10° (784 configurations in total), and the metrics (26), (27), and (30) are calculated within the quadratic region (31).
Table 1
Parameters of the simulation for case study 1
Symbol 
Value 
Unit 
Description 

n 
5 
Number of DOFs of the system 

nana 
2 
Number of actuated joints 

mimi 
1 
(kg) 
Mass of linkilinki, i=1,…,nℓi=1,…,nℓ 
ℓiℓi 
0.15 
(m) 
Length of linkilinki, 1=1,…,nℓ1=1,…,nℓ 
Icom,iIcom,i 
0.002 
(kg m^{2}) 
Rotational inertia for the ith link 
ττact=[τa1,τa2]Tττact=[τa1,τa2]T 
(N m) 
Input joint torques 

μsμs 
0.1 
Coefficient of (static) friction used for Scenario 2 
One example configuration can be seen in Fig. 4, where it is assumed that the object has a hundred times the mass of a link of the robot (i.e., κ=100κ=100). The first, second, and third columns represent the three scenarios depicted in Fig. 3, respectively. A lighter color represents a higher value of the depicted norm. Figure 4a shows the magnitude of the acceleration of the object aa→obj2aa→obj2. It can be seen that it barely changes regardless of the scenario (i.e., independently of the fact that the snake robot has or has not passive wheels, the object will accelerate the same given the same input). Figure 4b shows the magnitude of the acceleration of snake robot qq¨→s2qq¨→s2. This shows clearly that, even if the object’s acceleration is similar for all three scenarios, the behavior of the snake robot changes. The addition of passive wheels (second and third columns) increases the area where the snake robot’s slippage is minimal. Without passive wheels, the snake robot will slip in almost any direction of the input space.
The slippage ratio (30) gives quantitative information about the movement of the system and can be studied in the same way as the previous norms. Figure 4c shows the value of the slippage ratio for all three scenarios with several values for κκ for one configuration. It can be seen that sr→0sr→0 in the region where there is no contact with the object (i.e., the snake robot can move freely and therefore aa→obj2→0aa→obj2→0). It is interesting to see that in all three scenarios it is always possible to make the object move. However, Scenario 2 (nonideal passive wheels) has a limited region where the slippage ratio is high, compared to Scenario 3, where a whole region seems to give high values of slippage ratio. These regions in the input space are highlighted in Fig. 4b, c. Regions where the slippage of the snake robot are minimized tend to have higher slippage ratio.
Case study 2: Snake robot with three joints
The proposed framework and metrics can be applied to a snake robot with any number of joints. In this section, a snake robot with three joints is studied. However, studying the threedimensional input space could be cumbersome. Instead, the norms λλ→2λλ→2, aa→obj2aa→obj2 and slippage ratio (30) are studied as a function of the polar coordinates of the COM of the snake robot (COM,∠COM)(COM,∠COM) w.r.t. the contact point, as discussed in previous sections.
Figure 5a reports the result for the norm of constraint forces λλ→2λλ→2. The polar plots show the results for scenario 1, 2, and 3, respectively. The higher the value, the bigger the constraint forces. It can be seen that in scenario 1 (negligible friction) there is a clear trend for configurations with the COM of the snake robot at angles 90∘90∘ and −90∘−90∘ to have a higher impact on the wrench applied to the object. Although scenarios 2 and 3 report a higher norm of the constraint forces, this is due to the addition of passive wheels. From this figure alone, it is not possible to ascertain the impact on the object.
Figure 5b reports the result for the norm of the object’s acceleration aa→obj2aa→obj2. The polar plots report the results for scenarios 1, 2, and 3, respectively. It can be seen that the addition of passive wheels (even ideal ones) have little impact on the acceleration of the object. However, the configuration of the snake robot, parametrized with the polar coordinates of its COM, has a clear and meaningful impact on the acceleration of the object.
Although basic intuition would tell that the addition of constraints (i.e., passive wheels) should have an impact on the force applied to the object (through λλ→λλ→), and consequently on its acceleration aa→objaa→obj, this study shows that is not the case (at least, not that simply).
An important addition of this paper w.r.t. [15, 24] is the study of the slippage ratio sr. By studying the relationship between motions of both systems (snake robot and object), we can understand how the additional constraints have an impact on the system. A snake robot without passive wheels will slip as it pushes the object. Therefore, minimizing this motion while keeping a steady force on the object (and therefore producing an acceleration) is desirable. Figure 5c shows the result of the slippage ratio (30). It can be seen that in Scenario 1 (without passive wheels) the same trend as with the object’s acceleration appears. However, passive wheels (even nonideal ones) have a big impact on the slippage ratio (take notice of the change of scale).
To show more clearly the impact of the configuration of the snake robot on the acceleration and slippage of the system, Fig. 6 shows the best and worst configurations for the acceleration of the object Fig. 6a and for slippage ratio Fig. 6b. The results are summarized in Table 2. It can be seen that passive wheels (Scenario 2 and 3) have little impact on the acceleration, but a significant one on decreasing the slippage of the snake robot (sr→1sr→1).
Table 2
Results of the norms
Concept 
Scenario 1 
Scenario 2 
Scenario 3 

Worst aa→obj2aa→obj2 
0.000763215 
0.00107454 
0.00107454 
Best aa→obj2aa→obj2 
0.00785302 
0.00861168 
0.00861168 
Worst sr 
0.00138484 
0.00185713 
0.00185098 
Best sr 
0.0267348 
0.988638 
0.988643 
Discussion and future applications
Several assumptions have been made in this line of research, especially the number of contacts considered between objects, and their rigidity. This is because we are interested in giving a solution that is mathematically rigorous while guaranteeing uniqueness of solution. To include more contact points means to loose this in favor of robustness. For example, a penalty method (aka. virtual springs) or barrier functions may be considered, which is common for wholearm body manipulation (WAM) tasks. Although our assumptions are restrictive, it allows us to give a solid foundation for the research. Other models or considerations can be used for more realistic scenarios, but the rigidbody assumption used here allows to have a clear basis for comparison. Considering the gaps in knowledge regarding snake robots (as highlighted in “Background” section), we consider the model and results presented to be useful for moving research forward.
The holonomic constraints (e.g., constraints due to joints of the snake robot) are already encoded in the kinematic model of the snake robot presented in “Mathematical background” section. These holonomic constraints are described by the use of the robot’s Geometric Jacobians. Further distinction between holonomic and nonholonomic constraints is not necessary, as both can be expressed in the same unified manner (13), as mentioned in [21, 33].
The metrics and general framework presented can be used to analyze more complex tasks involving snake robots (or similar robotic systems) interacting with an object or the environment. This gives an opportunity to study both snake robots and robotic arms in the same framework, since the analysis is similar to the oftenused force/manipulability ellipsoids used to study robotic arms or hands [34].
However, the metrics presented in previous research do not consider the motion of the robot itself, since a fixedbase was always assumed. The framework presented in this paper can be applied to other mobile systems in a more complete manner than reported in the literature [28, 29]. More specifically, the analysis presented extends the concept of force or manipulability ellipsoids [26, 34], from the case of a fixedbase robot with endeffector, to a mobile robot without an endeffector. For a given task and configuration, an analysis can be carried out to find the optimal input (vector of joint torques) to minimize or maximize slippage of the system.
A few conclusions can be drawn from analyzing norms (26) and (27) on the input space (c.f. Fig. 4). In all scenarios, τ2τ2 which is the joint furthest away from the object has almost no effect on the acceleration of the object. However, the addition of passive wheels helps to anchor the snake robot and couples the effect of τ2τ2 on the system.
Conclusion
In this paper, a modeling and analysis framework for snake robots in contact with an external body has been presented. Results show that the addition of passive wheels has little effect on the wrench applied to the object and therefore, little change in its acceleration. However, the passive wheels do have an effect on the motion of the robot itself. In other words, under certain conditions the slippage of the robot can be minimized while pushing the object. This could be beneficial for pushing or manipulation tasks. To the best of our knowledge, this problem has not been fully studied with snake robots.
Declarations
Authors’ contributions
FR conducted mathematical analysis, programming to perform the simulations, and wrote the manuscript. SM supervised the research. Both authors read and approved the final manuscript.
Acknowledgements
This study was in part supported by “Strategic Research Foundation Grantaided Project for Private Universities (2013–2017) from Ministry of Education, Culture, Sports, Science and Technology, Japan, and RGIRO (Ritsumeikan Global Innovation Research Organization). We acknowledge and thank their support.
Competing interests
The authors declare that they have no competing interests.
Funding
This research did not receive funding from any particular organization.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Authors’ Affiliations
(1)