CC BY
5Russian Journal of Nonlinear Dynamics, 2019, vol. 15, no. 4, pp. 623-632. Full-texts are available at http://nd.ics.org.ru DOI: 10.20537/nd190421
MSC 2010: 93B18, 93B52
An Investigation of Motion of a Crawling Robot with Supports with Controllable Friction
L. Yu. Vorochaeva, S. I. Savin, A. S. Yatsun
This paper deals with the gait of a robot crawling on a horizontal rough surface with dry friction. A distinctive feature of the device is the presence of two supporting elements with a controlled coefficient of friction, which allows for alternately fixing the supports on the surface. As a result of numerical simulation, the patterns (laws) of influence on the motion characteristics of the mass-dimensional and control parameters of the robot, as well as the parameters of the supporting surface, are obtained, so that it is possible to find application in the design of the specified devices.
Keywords: crawling robot, supporting element, controlled friction coefficient, dry friction model
1. Introduction
Currently, scientists from different countries are conducting research aimed at creating crawling robots to carry out rescue, search, and reconnaissance operations in completely or partially destroyed buildings. The choice of exactly these crawling devices as such robots is due to their higher maneuverability and terrain crossing capacity over a non-deterministic surface in the presence of obstacles, as well as the ability to penetrate through narrow restricted spaces (slots, pipes, ventilation shafts) that are most often observed after natural and man-inflicted disasters.
Received May 31, 2019 Accepted July 27, 2019
This work was carried out within the framework of the Presidential Grant MK-200.2019.1.
Ludmila Yu. Vorochaeva mila180888@yandex.ru Andrey S. Yatsun ayatsun@yandex.ru
Southwest State University
ul. 50 let Oktryabrya 94, Kursk, 305040 Russia
Sergey I. Savin s.savin@innopolis.ru
Innopolis University
ul. Universitetskaya 1, Innopolis, 420500 Russia
In this connection, the articles [1-3] address issues related to the development of gaits of crawling robots, as well as issues of movement on surfaces with obstacles. For instance, in [1] the development of such an algorithm (of the gait) for the movement of a crawling snake-like robot on the surface with obstacles is described, so that the trajectory of movement of the robot is smooth. The study [2] presents a method for designing two gaits of a crawling robot moving through pipes and over a rough terrain, where the required configuration of the robot is formed by connecting curve segments. Spatial movement of the crawling snake-robot in the presence of obstacles is considered in [3], overcoming obstacles is provided by vertical and horizontal deformations of the robot's body in the presence of anisotropic friction between its links and the surface. The papers [4-6] present movement issues of the crawling robots in limited spaces, including pipes and along cables.
In this paper, we consider one of the simplest design schemes of a crawling robot, which consists of a link made in the form of a rod and supporting elements installed at its ends, which provides for switching between two values of the friction coefficient, one of which allows the support to be fixed on the surface, and the other, to slide over it. The aim of this work is to identify numerically the patterns of movement of the robot on a rough surface, which can later be used in the design of such devices.
2. Description of the crawling robot
In this paper we consider a crawling robot consisting of one link and two supporting elements with a controlled friction coefficient, moving along a horizontal rough surface in the plane Oxy (Fig. 1a).
The link of the robot is a rod with length l and mass m. The center of mass of the link is located at the point of symmetry of the link, denoted as point C. Let the position of the point O1 be set by the radius vector:
roi = (xoi yo1) > (2-1)
where xO1, yO1 denote the coordinates of the point O1 in the projections on the axis Ox and Oy, and the angle of rotation of the link relative to the positive direction of the axis Ox is equal to p. Then the position of the center of mass of the robot is determined by the formula
rC = (xO1 + l cos p/2 yO1 + l sin . (2.2)
O
o
yA
Mi ,
02
F/Pf
o
(a)
(b)
(c)
Fig. 1. a) A scheme of a crawling robot. b, c) Initial positions of the robot's motion stages: b) the first stage, c) the second stage.
Supporting elements are set at points Ox and O2 and allow for switching between two possible values of the friction coefficient: maximal fmax and minimal fmin [7]. We assume that the value fmax is sufficient for the support to be fixed on the surface, while when the friction coefficient fmin is observed, the contact element slides along the surface. At the support points, the drives are installed which create the moments Ml and M2 and allow the robot link to rotate in a horizontal plane relative to the fixed support. In the second support, at the same moment, the dry friction force Ffr arises.
3. Friction force models
In general terms, the dry friction force formula is written as
if £ = 0,
Ffr =
'Fd F
F
per
if (£ = 0) A (|Fo| < Fpred), if (£ = 0) A (|Fo| > Fpred),
(3-1)
where Fd stands for the friction force in dynamics, Fs denotes the static friction force, Fper means the transient friction force, Fpred is the ultimate friction force, Fo is the resultant of all the forces except for the force of dry friction applied to the contact point of the robot, and £ denotes the angular velocity of the point of the robot in which the friction force is applied.
In this paper, we confine ourselves to three options for the force of dry friction: the force of symmetric dry friction (Fig. 2a), the friction force with a reduced dynamic component (Fig. 2b) and the friction force with a continuously drooping characteristic (Fig. 2c). Table 1 shows formulas for each of the friction forces.
kFfr
iFfr o i
>Ffr 0 v
<P
(a)
(b)
(c)
Fig. 2. Models of the dry friction force: a) symmetric, b) with a reduced dynamic component, c) with a continuously drooping characteristic.
Table 1. Formulas of models of the dry friction force
Friction forces Symmetric With a reduced dynamic component With a continuously drooping characteristic
Fs -F0
F per -fminNsgn(F0) -f^inNSgn(F0) -fminNsgn(F0)
Fd -/min^Sgn^) -fiinNs gn(¥>) -fmmNe-a^ 1 sgn(^)
F pred f ■ N J mm P- N J mm f ■ N J mm
F0 Mlj2//
Table 1 uses the following notation: fmin stands for the dry friction coefficient; /min, /¿m denote the coefficients of the static friction and the slip, respectively, N = mg/2 is the normal robot support reaction, and a is the parameter determining the properties of a rough surface.
4. Gait of the crawling robot
By alternately controlling the friction coefficient in the contact elements, as well as intermittently feeding torques to the drives and turning them off, you can implement the gait of this robot, as shown in Figs. 1b, 1c. The robot's locomotion includes two stages. At the first stage (n = 1) the point O1 is fixed. The drive set in this support generates torque M1, which is formed until the equality p = p01 is true, where p01 indicates some given value of the rotation angle of the robot link, corresponding to the completion of the first stage. After that, the second stage of the movement begins (n = 2), point O2 is fixed on the surface, the drive installed in it generates the torque M2. The stage will be terminated when the condition p = —p02 is satisfied, where the angle value p02 corresponds to the completion of the second stage. Then the stages are repeated. It should be noted that the values p01 and p02 are determined by the robot control system based on the trajectory along which the object must move, for example, when p01 = p02 the movement will occur in a straight line, at p01 > p02 it will be counterclockwise rotation, and at p01 < p02, clockwise rotation. More explicitly the gait of the robot is illustrated in Table 2.
Table 2. Description of the crawling robot gait
Motion 01 Motion 02 fo 1 fo2 Termination condition
stage 1 - + fmax f . J mm
stage 2 + - f . J mm fmax
The differential equation of motion of the robot during each of the stages is written as
J(p = Mlt2 - Ffrl,
where J stands for the moment of inertia relative to the fixed support point. Let the moment of the drive in support 1 be changed according to the law
(4.1)
' M0 sin(t - Ti) if [M0 sin(t - Tï) < Mc] A [n = 1], Mi = \mc if [(Mo sin(t - TÎ) ^ Mc) V (Mi = Mc)] A [n = 1], ,0 if [n = 2],
(4.2)
« — 1
« —9
0 1 2 3 At
Fig. 3. Dependence graphs of Ml(t), M2(t).
where M0 denotes the moment amplitude, Tf is the moment of time at which the equality p = p01 holds, Mc stands for the constant drive torque for which the condition Mc < M0 holds
The moment of the drive in support 2 is described by the formula ( 0 if [n = 1],
M2 = i - M0 sin(t - T22) if [-M0 sin(t - T22) > -Mc] A [n = 2], (4.3)
[ - Mc if [(-Mo sin(t - T2) < -Mc) V (M2 = -Mc)] A [n = 2],
where Tf denotes the moment of time at which the equality p = p02 holds. The drive torque graphs are shown in Fig. 3.
5. Simulation results
Let us consider the rectilinear motion of the robot with the introduced models of friction forces with dimensionless values of the parameters: m = 1, l = 1, M0 = 1.5, Mc = 1.45,
fmin = fmn = °.2, fLn = a = °.4, p01 = p02 = P0 = °.5.
Simulation of the robot's motion allows for the identification of the influence on the characteristics of movement of the mass-dimensional (the length and the mass of the link) and control parameters (the value of the control moments, the laws of their changes, the angle of rotation of the link), as well as the parameters of the supporting surface (the value of the friction coefficient, the model of the friction force). Figure 4 shows the time dependence of the change in the coordinates of the extreme points of the robot and its center of mass in the projections on the axes Ox and Oy, as well as the trajectory graphs constructed for the symmetric dry friction model.
Fig. 4. Dependence graphs: a) x01(t), x02(t), xc(t), b) yQ1(i), yo2(t), yc(t), c) yoi(xoi), yo2(xo2) yc (xc).
According to the graphs, it can be seen that during the first stage of the motion the points O2 and C dislocate, and during the second the points O1 and C dislocate, moreover, the movement occurs along the arcs of circles, the radius of which for the extreme points of the robot is twice larger than the radius for its center of mass. During each of the stages there are areas corresponding to the stationary state of the robot. This is due to the laws of the moment variation in the drives used in the work and the model of friction force (dry friction). It should be noted that the qualitative character of the graphs of Fig. 4 will not change with the variation of the model of dry friction force, and the quantitative changes will be due to the values of the static friction coefficients fmin and the slip f^in for the friction force with a reduced
dynamic component and the coefficient value a for the friction force with a constantly drooping characteristic.
In Fig. 5 the graphs of time dependences of the friction forces are shown during one cycle of movement (of two stages), where the areas of the static state of the object are highlighted in color. At these time intervals, the graphs of the friction forces repeat the corresponding laws of change of the control moments generated by the drives. In dynamics, the values of friction forces are either constant (for the symmetric friction force and the friction force with a reduced dynamic component, and in the first case the value of the sliding friction force is equal to the largest value of the static friction force, and in the second case is fs/fd times less) either changed exponentially (for the friction force with a constantly drooping characteristic).
Ffr n= 1 n = 2
Ffr n= 1 n = 2
3 t
Ffr 71 = 1 , n = 2
Fig. 5. Dependence graphs of Ffr(t) in one cycle of motion with the models of friction forces: a) symmetric, b) with a reduced dynamic component, c) with a constantly drooping characteristic.
Based on these graphs, it can be concluded that the areas of the static position of the robot do not depend on the model of the friction force, but are determined by the value of the friction coefficient fmin = fmin, the larger this value, the wider the static area ts on the time interval of each stage Tf = T2 = T*. In addition, the boundaries of this area depend on the weight-dimensional parameters of the robot (the mass m and the length of the link l), and also on the moment Mc, generated by the drives, as illustrated in the diagrams of Fig. 6.
/min — fu
/min — fn
/min — fn
Fig. 6. Section (area) diagrams of fmin, fSin(Mc): a) m, l = 0.8, b) m, l = 1, c) m, l = 1.2, 1 — the range of values when the robot's movement is possible, 2 — the range of values at which the robot is in a stationary position.
Figure 6 shows the ranges of the coefficients fmin = fms in, at which the movement of the robot is possible (area 1) and in which the object is stationary (area 2, in it ts = T*) depending on the moment Mc for three mass values m of the link at l = 1 or, on the contrary, for three lengths l of the link at m = 1. The given graphs indicate that, as the value Mc increases, the maximum permissible value of the friction coefficient increases, at which the robot will move
along the surface. Moreover, the dependence /min, /min(Mc) is linear, the slope of a straight line to the axis of moments decreases as the mass or the length of the link increases.
Further simulation was carried out for such values of friction coefficients at which ts <T*. Its goal is to establish the patterns of influence of the parameters of the friction model on the time intervals ts and T*. To assess the ratio of the time stage and the time during which the object is stationary, we introduce the ratio ts and T*. To assess the ratio of the time stage and the time during which the object is stationary, we introduce the ratio
P
î*. rp* '
(5.1)
The corresponding graphs are shown in Figs. 7-9. For the symmetric dry friction model, the friction coefficient fmin is considered as the parameter under study. According to the graphs in Fig. 7, it can be seen that the time of one stage of movement increases on the curve as the friction coefficient increases, and the greater the value of the drive torque Mc, the lower the graph of the corresponding curve is located.
Fig. 7. Dependence graphs for symmetric dry friction force: a) T*(/min), b) ts(/min), c) p(/min) at 1 M = 1.9, 2 — M = 2.45, 3 — M = 2.9.
The time of the stationary state of the device in this case increases according to the laws close to proportional, the angle of inclination of the straight line corresponding to the largest value Mc, to the horizontal axis of the graph is the smallest (Fig. 7b). The ratio p increases with the increase of the friction coefficient, reaches its highest value and then decreases, the nature of the curves is monotonic (Fig. 7c). With the increasing Mc value the highest parameter value p is observed at a higher value fmin and this value itself increases in number.
For the other two models, the corresponding graphs are drawn at the moment of the drive Mc = 1.9. For the model of friction force with a reduced dynamic component, the parameter under study is the coefficient of sliding friction fmin. The time of each stage increases with the increase of /min according to a curvilinear law, moreover, the greater the values of the coefficient of static friction /min, the higher is the corresponding curve (Fig. 8a). The time values ts do not depend on the coefficient /min, and are determined only by the coefficient /min, the greater the coefficient, the greater the rest time (Fig. 8b). The parameter p decreases with the increase of /min on the curve, the convexity of which is directed upwards, in this case there is no peak value observed (Fig. 8c).
The graphs in Fig. 10 for the friction force with a constantly drooping characteristic are built for the variable parameter a at three different values of the friction coefficient /min. The type of the dependences T*(a) differs from that given above in Figs. 7 and 8, herein, there is a decrease in the time of the stage with the increase of the parameter a on the curves, close to
/
1 N /
2
v3
3.5 3 2.5 2 1.5
1 0.2 0.4 0.6
fd■ J mm
(a)
Fig. 8. Graphs of dependences for the friction force with a reduced dynamic component: a) T*(f^nin), b)
ts am in), c) p(im in) at i — fm in = 0.57,2 — fm in = 0.5,3 — fm in = 0.45.
Fig. 9. Dependence graphs for the friction force with a constantly drooping characteristic: a) T*(a) b) ts(a), c) p(a) at 1 — ^ = 0.57, 2 — ^ = 0.5, 3 — ^ = 0.45.
the exponentials, and with the increase of the friction coefficient fmin the corresponding curve is located higher relative to the horizontal axis of the graph (Fig. 9a).
The time of the object at rest (stationary state) (Fig. 9b), as well as for the force with a reduced dynamic component, does not depend on the variable parameter (in this case, the parameter a). The ratios of p with the increasing parameter a are increasing in the monotonic type according to curvilinear laws, the convexities of which are directed upwards; the ratio p reaches the largest values at the highest of the considered parameter values of a (Fig. 9c). Based on the graphs of Figs. 8-10, it can be concluded that for the most rapid movement of the robot on the surface it is necessary to ensure the minimum possible values of the friction coefficient fmin for the symmetric friction force, the friction coefficients fmin and fmn for the friction force with a reduced dynamic component, the friction coefficient fmin and the parameter a for the friction force with a continuously drooping characteristic.
Another variable parameter affecting the motion characteristics of the device is the value of the angle . To assess its impact, we have drawn graphs of the time of motion, T, it takes the robot to travel the distance s = 1, in the event of variation of with the three models of the friction force (Fig. 10).
The nature of these dependences does not depend on the model of the friction force; they are due only to the numerical characteristics of the curves. The graphs can be divided into three sections: in the first of them (0 ^ < 0.5) there is a decrease in time T under the laws close to proportional, in the second (0.5 ^ < 1.4), the motion time does not depend on the angle ^>0, these areas feature horizontal lines, in the third (1.4 ^ ^ 1.57), a decrease in the motion time is observed again as the angle grows.
0.5 1 1.5 ip0
Fig. 10. Dependence graphs of T(y0) with the friction models: 1 — symmetric, 2 — with a reduced dynamic component, 3 — with a constantly drooping characteristic.
6. Conclusion
This paper considers the movement of a crawling robot formed by a link in the form of a rod and two supporting elements with a controlled friction coefficient due to the change of the supporting surfaces. The gait of the device is proposed, consisting of alternate fixation of the supports and the rotation of the link relative to the fixed support. A simulation has been carried out of the movement of the robot on a horizontal rough surface, the interaction of which with the supports of the robot is described by three models of dry friction force: symmetric, with a reduced dynamic component, and with a constantly drooping characteristic.
As a result of this research, the influence on the motion characteristics of the mass-dimensional parameters of the device, control values and the parameters of the supporting surface has been established. It was revealed that, to increase the velocity of the robot when moving over the surface, it is necessary to minimize its mass m and the length of the link l, as well as to use the largest possible moments Mc of the drives and move at each stage to the largest possible angles . In addition, it has been established that the time of the robot's motion to a specified distance is reduced while minimizing the friction coefficients fmin, fmin and f^in and the parameter a for the appropriate friction models. The diagrams of the static and dynamic areas of the robot are drawn depending on the mass-dimensional and control parameters, it is established that the qualitative and quantitative characteristics of these diagrams are determined by the value of the static friction coefficient fmin or fmin, but do not depend on the parameters of the model of the friction force, manifested in dynamics. The revealed regularities can be used in the development of the device design, the elaboration of the laws for controlling its drives, and in choosing materials for the supporting elements.
References
[1] Conkur, E.S. and Gurbuz, R., Path Planning Algorithm for Snake-Like Robots, Inform. Technol. Contr., 2008, vol. 37, no. 2. pp. 159-162.
[2] Takemori, T., Tanaka, M., and Matsuno, F., Gait Design for a Snake Robot by Connecting Curve Segments and Experimental Demonstration, IEEE Trans. on Robotics, 2018, vol. 34, no. 5, pp. 13841391.
[3] Fu, Q., Mitchel, Th., Yi, N., Gart, S., and Li, Ch., Snake Robot's Poor 3-D Obstacle Traversal Reveals Snake's Better Stability Mechanisms, APS March Meeting 2018, 2018, vol.63, no. 1.
[4] Douadi, L., Spinello, D., Gueaieb, W., and Sarfraz, H., Planar Kinematics Analysis of a Snake-Like Robot, Robotica, 2014, vol. 32, no. 5, pp. 659-675.
[5] Wei, W., Yu, J.X., Lin, W.H., Mai, R.Sh., Chen, Ch.X., He, Zh.Y., and Li, Zh., Structure Design of Climbing Snake-Like Robot for Detection of Cable-Stayed Bridge, Appl. Mech. Mater., 2014, vol.598, pp. 610-618.
[6] Huang, Ch.-W., Huang, Ch.-H., Hung, Y.-H., and Chang, Ch.-Y., Sensing Pipes of a Nuclear Power Mechanism Using Low-Cost Snake Robot, Adv. Mech. Eng., 2018, vol. 10, no. 6, 8 pp.
[7] Vorochaeva, L. Yu., Yatsun, A. S., and Yatsun, S. F. Simulation of the Motion of a Five-Link Crawling Robot with Controlled Friction on a Surface Having Obstacles, J. Comput. Syst. Sci. Int., 2017, vol. 56, no. 3, pp. 527-552.
