Research Article  Open Access
Adaptive Path Following and Locomotion Optimization of SnakeLike Robot Controlled by the Central Pattern Generator
Abstract
This work investigates locomotion efficiency optimization and adaptive path following of snakelike robots in a complex environment. To optimize the locomotion efficiency, it takes energy consumption and forward velocity into account to investigate the optimal locomotion parameters of snakelike robots controlled by a central pattern generator (CPG) controller. A cuckoo search (CS) algorithm is applied to optimize locomotion parameters of the robot for environments with variable fractions and obstacle distribution. An adaptive path following method is proposed to steer the snakelike robot forward and along a desired path. The efficiency and accuracy of the proposed path following method is researched. In addition, a control framework that includes a CPG network, a locomotion efficiency optimization algorithm, and an adaptive path following method is designed to control snakelike robots move in different environments. Simulation and experimental results are presented to illustrate the performance of the proposed locomotion optimization method and adaptive path following controller for snakelike robots in complexity terrains.
1. Introduction
Inspired by biological snakes, snakelike robots have received significant attention to take the place of human work in some special situations such as exploration, monitoring, and surveillance tasks [1]. Similar to the stability and flexibility of biological snakes, snakelike robots can be used in challenging environments [2].
Since the first snakelike robot was developed by Hirose in 1972, many researchers have developed many snakelike robots with different mechanisms. Some snakelike robots are configured with active joints and passive wheels such as robots in [3, 4]. Some snake robots are connected by some joints that can bend, lengthen, and shorten, such as robots in [5, 6]. Some snakelike robots are made up by active joints and active wheels, such as robots in [7, 8]. Some snake robots are consisted of active wheels and passive joints [9]. Some robots are combined by joints and crawlers [10].
Several control approaches for robots with many degrees of freedom (DOFs), such as quadruped robot and snakelike robot, have been proposed in literatures [11, 12]. In the past years, researchers find that the central pattern generator (CPG) inside the spinal cord could generate rhythmic locomotion without sensory feedback or regulation signals from brain [13]. Matsuoka was a pioneer of CPG research and he proposed the first practical model in 1987 [14]. In the past years, CPG was successfully used to control robots, especially in control for bionic robot such as snakelike robots [15]. Rhythmic locomotion, such as lateral undulation and sidewinding, can be controlled by a CPG network connected by neuron oscillators [16, 17]. In [18], a doublechain Hopf CPG network is used to control serpentine and sidewinding gaits of a snake robot. In [19], a unified neural oscillator network is researched to control three different types of locomotion by adjusting several parameters. Smooth transition and autonomous locomotion control methods are investigated in [20–22]. However, some important performances such as locomotion efficiency and path following accuracy of snakelike robots are not considered in these studies.
There are many studies aiming at analysing and optimizing the locomotion efficiency of snakelike robots. The extremely important problem that should be solved is the locomotion efficiency of robots. In [12], three different gaits of a snakelike robot are analysed based on its dynamic model. The tradeoff between locomotion speed and energy efficiency is studied for several gaits. In [23], the relationships among the gait parameters, energy consumption, and forward velocity of different gaits for underwater snake robots are investigated, and the most efficient motion pattern is selected by empirical rules according to simulations.
In [24], gait parameters are optimized to stabilize the head joint and speed of a snakelike robot. However, the energy efficiency is not taken into account. In [25], a genetic algorithm is adopted to select parameters of serpentine gait parameters for a snake robot, but other performances are not considered except the speed of the robot. In [26], forward speeds of amphibious snakelike robots are optimized by the Powell’s method. However, the frequency of locomotion is not taken into account in the optimization algorithm [26]. A framework is researched to optimize locomotion of underwater snake robots and their energy efficiencies are investigated in [27]. To the best of our knowledge, there is no general optimization framework for snakelike robots in a random environment with different fractions and obstacles.
To realize the practicability of snakelike robots in various tasks, their path following in unknown terrains has been researched. In [28], a dynamic feedback controller is proposed to control the orientation of a snake robot to a desired angle. The theoretical approaches are validated through experimental results in [29]. In [30], a planar snakelike robot can move forward and traverse the reference path with a desired velocity controlled by a feedback control strategy that enforces virtual constraints. In [31], path following of snake robots is researched by using approximate dynamic programming and neural networks.
Inspired by the path following control of vessels on the sea in [32, 33], an integral line of sight (LOS) path following controller is designed to control an underwater snakelike robot go along a straight line based on a dynamic model in [34]. However, some of its key parameters are chosen arbitrarily. It is thus difficult to realize efficient and accurate path following due to the mechanism complicacy and high degree of freedoms of snakelike robots.
In this work, a doublechain Hopf GPG is used to replace the sine function to generate joint reference functions to avoid sudden transition of locomotion. To reduce energy consumption and extend autonomous working hours of snakelike robots, optimization of locomotion efficiency is considered. The energy consumption and forward speed are taken into account to investigate the optimal gait parameters. To explore the effect factors of locomotion efficiency, we have conducted various simulations. The simulation results show that the frequency and amplitude of a CPG controller are relevant to the efficiency. To find the most efficient gait parameters combination of the Hopf CPG, a CS algorithm is used to find the optimal gait parameters for different environment. An adaptive path following approach is proposed to improve the accuracy and the efficiency of path following. The approach is used to control a snakelike robot following curves and straight lines in narrow spaces with disorderly obstacles. Position, direction and forward velocity, efficiency, and error accumulation of snakelike robots are taken into account in the control of path following. In addition, a control framework includes a Hopf CPG network, an exponentially stable controller, a locomotion efficiency optimization algorithm, and an adaptive path following method. This work also includes more accurate and detailed simulation and experiment results about efficiency optimization and adaptive path following in complex environments.
This paper is organized as follows. The mathematical model of a snakelike robot and a CPG network used in this paper are presented in Section 2. Relations among locomotion efficiency, gait parameters, and environment factors are studied and the locomotion efficiency is optimized by the CS algorithm in Section 3. An adaptive LOS path following controller is proposed in Section 4. Extensive simulations are conducted to verify the efficiency and accuracy of the proposed path following method in Section 5. Finally, Section 6 concludes the paper.
2. Mathematical Model
Mathematical model of snakelike robot has attracted much attention [34, 36, 37]. In this paper, we briefly present the mathematical model of a snakelike robot that will be used in efficiency optimization and path following control. Our control design is based on the model presented in [36].
2.1. Mathematical Model of SnakeLike Robot
As shown in Figure 1, the snakelike robot is connected by links. The head position is , orientation of link is , and pitch angle of joint is . The mass of a link is set as and the moment of inertia of a link is defined as . To simplify computation, the center of mass (CM) is assumed to locate at its center point. Other notations and equations of the mathematical model of the robot in this paper are the same as those in [36].
In this dynamic model, viscous friction is assumed to act only on the CM of a link. The friction act of link is calculated aswhere and are viscous fraction coefficients between links and the ground in the tangential and normal directions, respectively. In this paper, the snakelike robot is connected by 10 links, the length of a link is , and the mass of a link is . The ground is flat and fraction coefficients are and . The rotation matrix of link is defined asThe model can be written in the standard form of a controlaffine system by defining , , , , and . Thus we havewhere and are computed byand matrices and are defined in [36].
The joint dynamics of a snakelike robot is given by . The new control input of the model is defined according to the exponentially stable joint controller: Reference joint angle of joint is defined aswhere is amplitude, is frequency of reference joint functions, and is the phase offset that affects the shape of a snakelike robot by changing the number of shapes. The relation between the shape and can be given aswhere is the number of shape in the horizontal plane and is the number of links.
However, actuators will be damaged if the snakelike robot is controlled by (7) directly when the reference function transits suddenly. To avoid jerky movements and physical damages of actuators, the CPG controller is applied to generate reference functions of joints.
2.2. Dynamic Model and Topological Structure of CPG Network
Some neural oscillator models were used to control the snakelike robots in recent years, such as Kuramoto oscillator model [38], Van Der Pol oscillator [16], and Hopf oscillator model [39]. However, the Kuramoto model needs big calculation and the Van Der Pol model converges slowly. Hence, the Hopf neural model is used to control a snakelike robot in this work due to its stability and fast convergence. The dynamics of a Hopf oscillator is given aswhere , is the radius of a limit cycle, is the convergence rate to the limit cycle, and decide the center of the limit cycle of the Hopf oscillator, and oscillate around and respectively, is a bifurcation parameter of the Hopf oscillator and in this paper, is accepted as the output signal, and the offset of the signals of is decided by . Solutions of and are sinusoidal in steady states. In this paper, notations and parameters of dynamics of the Hopf oscillator follow those in [35].
In this paper, the Hopf CPG network is connected by two series of Hopf oscillators. As shown in Figure 2, the left column (blue) and the right column (red) oscillators control horizontal and vertical movement snakelike robot, respectively. Parameters , , , and shown in Figure 2(a) are connection weights among neighboring oscillators, and in Figure 2(b) is the phase offset between two neighboring oscillators. The phase difference between horizontal and vertical oscillator is set as . Dynamics of the Hopf CPG network are given aswhere and are outputs of horizontal and vertical oscillators, respectively. and control the frequency and amplitude of the limit cycle. The rotation matrix is computed bywhere and are reference values of the yaw actuators and pitch actuators of joint . In this paper, reference angles of joint are computed by
(a)
(b)
As known from previous literatures, CPG signals change suddenly and discontinuously when parameters change abruptly. That leads to unstable locomotion and damage to motors of joints. Hence, we use the smooth gait transit method proposed in our previous work [35] to avoid sudden transition. The smooth transit function is defined aswhere amplitude and offset are computed aswhere and are predifined and desired values of , , and are the start and the end time of parameter changing.
3. Efficiency Optimization for SnakeLike Robots in Different Environments
As known from previous literatures [16–19], the Hopf CPG network is capable of controlling the rhythmic locomotion of snakelike robots. Nature snakes can always adjust their body shapes to move efficiently in different terrains. In order improve the adaptability of snakelike robots to complex environments, the lateral undulation locomotion with several combinations of CPG parameters and fractions is discussed. Herein, eight representative environments are investigated to obtain the efficient locomotion. Snakelike robots controlled by different CPG parameters are simulated to find out the connection between CPG parameters and robot locomotion efficiency.
3.1. Locomotion Efficiency Gait Parameters
To evaluate the locomotion efficiency of snakelike robots, we take their speed and energy consumption into account. To simplify the optimization objective, we compute locomotion efficiency as follows: where is the torque of the actuator of joint and is the distance of the snakelike robot moving in one period. is arbitrary time and , where is a period of locomotion. is the sum of squares of torques of all joints from to . It is the energy consumption in one period.
In this paper, different gaits are controlled by a Hopf CPG controller. The controller combined with the efficiency critic function (20) can generate uniform outputs with desired amplitudes and phase differences without additional adjustments.
To explore the relation between locomotion efficiency and parameters, locomotion controlled by different gait parameters in the same flat ground with tangential fraction and normal fraction is simulated. We simulate 12 proportional combinations of frequency and amplitude , and other parameters are constants. From the simulations results of Table 1, locomotion efficiency is affected by and of the CPG network. Thus, the relevant amplitude and frequency can be adjusted to obtain the efficient locomotion.

3.2. Locomotion Efficiency of SnakeLike Robots in Different Environments
Frictions in normal and tangential directions are important for lateral undulation of snakelike robots. To find out the relation between frictions and locomotion efficiency, we simulate eight different situations with the same locomotion parameters. In the eight situations, fractions in tangential and normal directions are set as (0.1, 0.3), (0.1, 0.4), (0.1, 0.5), (0.1, 1), (1, 3), (1, 4), (1, 5), and (1, 10), respectively. The first four situations are set with same friction coefficient in the tangential direction and increasing friction coefficient in the normal direction. The last four situations are set with ten times friction coefficient of the first four situations in tangential and normal directions.
As shown in Table 2, locomotion efficiencies are different in eight situations. It can be found that the efficiency is proportional to the fraction ratio . The locomotion efficiencies are similar to the same fraction .

Except for the fraction, a snakelike robot demonstrates different performances in different environments with different slop and different passable widths. As shown in Figure 3, a snakelike robot should transit body shape and locomotion flexibly to cross the cluttered environment covered with obstacles. The optimized efficiency locomotion parameters are selected from computing a large number of combinations of locomotion parameters. However, the optimized velocity is slow and the computation is too complex. It is necessary and significant to take advantage of an intelligent optimization algorithm to optimize parameters fast when we face a complex environment.
3.3. Parameter Optimization Using CS Algorithm
Traditional parameter selection method selects desired parameters by comparing all combinations of the amplitude and the frequency in the allowed ranges. To optimize the locomotion parameters fast and effectively, the CS algorithm [40–42] is used to explore the suitable optimal amplitude and frequency while the snakelike robot is moving in a new environment. The optimization algorithm used in this research is inspired by the obligate brood parasitic behaviour of some cuckoo species in combination with the flights behaviour of some fruit flies. The optimization objection is maximized locomotion efficiency in (20), to given ranges of these parameters.
As shown in Figure 3, the width of the wholeshape of a snakelike robot should be less than the width of a passable path. The amplitude of the wholeshape is decided by the width of a path, and is relevant to in (7). The ranges of should be derived according to the width of passable path to guarantee that snakelike robot to crosses a passable path without colliding obstacles on both sides of the latter.
In order to guarantee the constringency of searching, the optimal individual retention method is accepted in the CS algorithm. New values of our interested parameters are searched by using the CS algorithm aswhere and is the step size, which is related to the scales of the optimization problems and should be greater than zero.
In this paper, is set just according to the scale of amplitude to reduce the computation complexity since the values of and are approximate. is computed by . The product means the entry wise multiplication of two vectors. The step length of Ley flight is computed by a Ley distribution: while the Ley distribution has an infinite variance with an infinite mean.
4. Adaptive LOS Path Following Control
To optimize the accuracy and efficiency of path following of a snakelike robot executing special tasks in a narrow space, we propose an adaptive path following control scheme. As shown in Figure 4, an adaptive control frame is formed by a CPG locomotion generator, a direction controller based on the adaptive LOS guidance law, locomotion regulator, and an energy efficiency critic component. A closedloop control scheme based on the LOS law and the efficiency critic element is summarized to adjust the locomotion of a snakelike robot adaptively according to energy efficiency in different terrains. In the control scheme, information of environment including fractions and slop is measured by sensors, and appropriate locomotion parameters including and are selected by using the CS algorithm according to efficiency. The selected locomotion parameters are transported to the CPG network that generates reference swing angles of joints, and torques of actuators of the robot are computed by the exponentially stable controller as (6). The modules of the control scheme are presented in above sections.
4.1. LOS Guidance Law
In this paper, a LOS guidance law is applied to the curve path following of a snakelike robot. As shown in Figure 5, the tangent line of a reference path at point which is the nearest point to the CM of a snakelike robot is the desired path. The distance from the CM of a robot to is defined as , which corresponds to the cross track error, and the direction angle of the snakelike robot is represented by .
The path following control objective is to control a snakelike robot to converge to a reference path and subsequently move along the path at nonzero forward velocity. The objectives of path following can be formalized aswhere is an analytic function of a reference path. Snakelike robots cannot track a reference path flawlessly due to the oscillatory gait pattern. Our control objectives imply that the CM snakelike robot can have steady state oscillation around the reference path.
In this section we study LOS pathfollowing control for snakelike robots as shown in Figure 5. Their reference orientation is computed using the LOS guidance law aswhere is the lookahead distance that influences the convergence rate to the desired path, is the position of the CM of the robot, and is the nearest point in the reference path to the CM of the robot. As known from Section 2, the direction of the robot can be controlled by phase offset . Therefore, we can steer the heading direction according to the LOS angle in (26), and the joint angle offset is computed aswhere controller gain and heading direction are
4.2. Adaptive Path Following Based on LOS
Previous studies on path following of snakelike robots have focused on undulation locomotion with constant gait parameters in a stable environment, and gain in (29) and lookhead distance in (26) are fixed. The error accumulation and the energy efficiency of path following are not to be considered. Although snakelike robots can follow reference paths in a stable environment, it is hard to guarantee their path following performances of in complex environments especially when they face some uncertain factors such as fractions and distribution of obstacles.
To adapt to complex environments, this work proposes an adaptive path following approach for snakelike robots. In this path follow approach, gain and lookahead distance are adjusted adaptively according to the relation between their position and orientation. To control them to follow desired paths accurately and efficiently in different situations, we analyse their characteristics when they move in eight situations. Taking an example of following reference path , it can be divided into eight typical different situations according to the relative position of their CM to the reference path and the difference value between and as shown in Figure 6.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
After analysing these situations, we find that several factors are influential to values of and including the difference between the CM of a robot and the reference path and difference between and . Taking these influential factors into account, adaptive parameters and are computed bywhere and are initial values of and and and are constants that influence the amplitude of oscillation of a snakelike robot around the reference path. These parameters can be adjusted according to width and the requirements of following accuracy. In this paper, parameters are set as , , , , , , , , , , and .
Energy efficiency and error accumulation are taken into account beside the control objectives (26)(28). The oscillation of a snakelike robot around the reference path is also considered. The error accumulation can be computed as
5. Evaluation
To test the performance of the locomotion efficiency optimization method proposed in Section 3 and the adaptive LOS path following control scheme described in Section 4, we have conducted extensive simulations. The dynamic model and parameters of a locomotion controller of a snakelike robot are set in Section 2, and the reference angles of active joints are generated by the CPG network as described in Section 3.
5.1. Locomotion Efficiency Optimization in Complex Environment
In this simulation, the locomotion parameters of a snakelike robot in a complex environment are optimized by the CS algorithm. The parameters of the robot and control parameters are the same as before. When the robot moves in a new environment, the ground fraction is measured and the scales of amplitude and frequency are derived according to the space width. Then, the optimal combination of and is selected via the CS algorithm. The step size is and in the CS algorithm. As shown in Table 3, gait parameters of the robot are optimized by the CS algorithm in eight situations with different fractions. As shown in Figure 7, the gait parameters of the robot are selected by the combination method that chooses combinations of amplitude and frequency at random, where and . Compared with random combination method, the CS algorithm improves the accuracy and the speed of parameter optimization.

5.2. Adaptive Path Following of the SnakeLike Robot
The simulation is actualized in Matlab R2016b in a computer running Windows 10. The dynamics model is computed by the ode45 solver with a relative and absolute error tolerance of . The model is simulated with the controller given by (6), and controller gains are and . Reference angles of joints are computed by (7) and gait parameters are set as , , and .
Using the proposed path following method in Section 4, a snakelike robot is controlled to follow the reference path including a horizontal straight line and a sine curve. As shown in Figure 8, snakelike robot controlled by our approach moves forward and along with the horizontal line more accurately compared with using the LOS path following method. As shown in the Figure 8(b), the error accumulation is reduced with the adaptive path following method proposed in this paper. As shown in Figure 8(c), the efficiency of the robot controlled by proposed the method is predominant compared with the traditional LOS path following method. As shown in Figure 9, snakelike robot controlled by our approach moves forward and along with the sine curve path more accurately compared with using the classical LOS path following method. As shown in the Figure 9(b), the error accumulation is reduced with the adaptive path following method proposed in this paper. As shown in Figure 9(c), the efficiency of the robot controlled by proposed the method is predominant compared with the traditional LOS path following method. The simulation results are proposed to show that the position of the robot converges nicely to the reference horizontal straight line and sine curve paths.
(a)
(b)
(c)
(a)
(b)
(c)
To verify the proposed method, a simulator is implemented in Virtual Robot Experimentation Platform (VREP EDU Version), as shown in Figure 10. In the simulation, we modify the structure and control method based on the model of University of Castilla la Mancha. We equip a pair of passive wheels for each link to realize anisotropic friction between the robot and the ground. The actuators are installed on the horizontal and vertical joints to imitate natural snakes swinging from side to side. Inside VREP, we choose the bullet as the physics engines and set the time step at 50ms. The locomotion parameter optimization and adaptive path following methods are used in this simulation. The sine transition method is used to avoid sudden changes while locomotion parameters adjusting. As shown in Figure 10, there are some irregular obstacles in the environment and the robot through differentwidth paths adaptively by adjusting locomotion parameters. As shown in Figure 10(a), the passable path is narrow and the amplitude of the robot is small. The amplitude of the robot is transited adaptively according to the width of path, as shown in Figures 10(b) and 10(c).
(a)
(b)
(c)
(d)
Extensive prototype experiments are conducted to demonstrate that our method can adjust locomotion parameters adaptively and steer a snakelike robot toward and along the desired path. Our snakelike robot is connected by five pitch modules and five yaw modules. Each module is equipped with a steering gear (SC5115). The steering gear can rotate between and the angles of steering gears can be read by their potentiometers. Modules communicate with each other via RS485 bus. To increase surface friction, our snakelike robot is covered with plastic sheeting. As shown in Figure 11, the snakelike robot adjusts its body shape adaptively while moving across the narrow path with some irregular obstacles on both sides. As shown in Figure 11(a), the path in front of the robot is complex and distributed with some irregular obstacles on sides, and widths of path is different. As shown in Figure 11(b), the path is narrow and the robot moves toward with a small amplitude. As shown in Figure 11(c), the path is wider and the robot adjusts amplitude according to the path and selects corresponding frequency optimized energy efficiency. As shown in Figure 11(d), the path becomes narrow and the robot reduces the amplitude and selects efficient frequency. In this experiment, the robot can be controlled to go toward and along the desired horizontal line in the middle of the road and keep efficient locomotion using parameters according to paths and optimized by the CS algorithm.
(a)
(b)
(c)
(d)
6. Conclusion and Future Works
In this paper, a Hopf CPG network is used to replace the serpentine function to generate reference functions of joints of a snakelike robot to improve the smooth degree of control signal. To improve the energy efficiency of such robots in different environment, a CS algorithm is used to optimize locomotion parameters of the Hopf CPG network. An adaptive path following controller based on the improved LOS guidance law is designed to improve the adaptability of the robot to a complex environment with variable fractions and path widths. The adaptive LOS path following law is combined with a locomotion efficiency optimization system and CPG controller to steer the robot to the reference path accurately and efficiently. Simulation and experimental results show that the proposed methods optimize robot locomotion efficiency and steer the robot moving along with desired paths including straight and curve paths.
As for future work, we will equip more sensors on the snakelike robot and intend to study adaptive control and multiobjective locomotion optimization of the robot for more complex environment. In particular, the proposed adaptive efficiency optimization and path following control approach will be used to other 3D gaits of snakelike robots.
Data Availability
Some key codes and dates of my paper can be accessed by https://github.com/zhang3dong/Adaptivepathfollowingofsnakelikerobots.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work is supported by the National Natural Science Foundation of China (Grant nos. 51575034, 91748102, and 91848103), State Key Laboratory of Robotics (2018O15), and State Key Laboratory of Digital Manufacturing Equipment and Technology (DMETKF2018021).
Supplementary Materials
The supplementary video includes simulation and prototype experiment of our snakelike robot moves in a complex environment controlled by methods proposed in this paper. (Supplementary Materials)
References
 J. K. Hopkins and S. K. Gupta, “Design and modeling of a new drive system and exaggerated rectilineargait for a snakeinspired robot,” Journal of Mechanisms and Robotics, vol. 6, no. 2, p. 021001, 2014. View at: Google Scholar
 S. Hirose, Biologically Inspired Robots: SnakeLike Locomotors and Manipulators, Oxford. Univ. Press, London, UK, 1993.
 M. Moris and S. Hirose, “Threedimensional serpentine motion and lateral rolling by active cord mechanism ACMR3,” in Proceedings of the 2002 IEEE/RSJ International Conference on Intelligent Robots and Systems, pp. 829–834, Switzerland, October 2002. View at: Google Scholar
 A. Crespi, A. Badertscher, A. Guignard, and A. J. Ijspeert, “AmphiBot I: An amphibious snakelike robot,” Robotics and Autonomous Systems, vol. 50, no. 4, pp. 163–175, 2005. View at: Publisher Site  Google Scholar
 H. Ohno and S. Hirose, “Design of slim slime robot and its gait of locomotion,” in Proceedings of the 2001 IEEE/RSJ International Conference on Intelligent Robots and Systems, pp. 707–715, USA, November 2001. View at: Google Scholar
 S. Sugita, K. Ogami, M. Guarnieri et al., “A study on the mechanism and locomotion strategy for new snakelike robot active cord mechanism Slime model 1 ACMS1,” Journal of Robotics Mechatronics, vol. 20, no. 2, pp. 302–310, 2008. View at: Publisher Site  Google Scholar
 H. Yamada and S. Hirose, “Development of practical 3dimensional active cord mechanism ACMR4,” Journal of Robotics & Mechatronics, vol. 18, no. 3, pp. 305–311, 2006. View at: Publisher Site  Google Scholar
 H. Yamada, S. Chigisaki, M. Mori et al., “Development of amphibious snakelike robot ACMR5,” in Proceedings of the International Symposium in Robotics, pp. 617–624, 2005. View at: Google Scholar
 H. Kimura and S. Hirose, “Development of Genbu: Active wheel passive joint articulated mobile robot,” in Proceedings of the 2002 IEEE/RSJ International Conference on Intelligent Robots and Systems, pp. 823–828, Switzerland, October 2002. View at: Google Scholar
 M. Arai, Y. Tanaka, S. Hirose, H. Kuwahara, and S. Tsukui, “Development of SouryuIV and SouryuV: serially connected crawler vehicles for inrubble searching operations,” Journal of Field Robotics, vol. 25, no. 12, pp. 31–65, 2008. View at: Publisher Site  Google Scholar
 D. Gong, P. Wang, S. Zhao, L. Du, and Y. Duan, “Bionic quadruped robot dynamic gait control strategy based on twenty degrees of freedom,” IEEE/CAA Journal of Automatica Sinica, vol. 5, no. 1, pp. 382–388, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 R. Ariizumi and F. Matsuno, “Dynamic Analysis of Three Snake Robot Gaits,” IEEE Transactions on Robotics, vol. 33, no. 5, pp. 1075–1087, 2017. View at: Publisher Site  Google Scholar
 S. Grillner, T. Deliagina, A. El Manira et al., “Neural networks that coordinate locomotion and body orientation in lamprey,” Trends in Neurosciences, vol. 18, no. 6, pp. 270–279, 1995. View at: Publisher Site  Google Scholar
 K. Matsuoka, “Mechanisms of frequency and pattern control in the neural rhythm generators,” Biological Cybernetics, vol. 56, no. 56, pp. 345–353, 1987. View at: Publisher Site  Google Scholar
 L. Wen, T. Wang, G. Wu, J. Liang, and C. Wang, “Novel method for the modeling and control investigation of efficient swimming for robotic fish,” IEEE Transactions on Industrial Electronics, vol. 59, no. 8, pp. 3176–3188, 2012. View at: Publisher Site  Google Scholar
 H. Yu, W. Guo, J. Deng, M. Li, and H. Cai, “A CPGbased locomotion control architecture for hexapod robot,” in Proceedings of the 2013 26th IEEE/RSJ International Conference on Intelligent Robots and Systems: New Horizon, IROS 2013, pp. 5615–5621, Japan, November 2013. View at: Google Scholar
 N. M. Nor and S. Ma, “Smooth transition for CPGbased body shape control of a snakelike robot,” Bioinspiration & Biomimetics, vol. 9, no. 1, p. 016003, 2014. View at: Google Scholar
 Z. L. Wang and Q. Gao, “CPGinspired locomotion control for a snake robot basing on nonlinear oscillators,” Journal of Intelligent Robotic Systems, vol. 85, no. 2, pp. 1–19, 2016. View at: Google Scholar
 S. Manzoor and Y. Choi, “A unified neural oscillator model for various rhythmic locomotions of snakelike robot,” Neurocomputing, vol. 173, pp. 1112–1123, 2016. View at: Publisher Site  Google Scholar
 Z. Bing, L. Cheng, G. Chen, F. Röhrbein, K. Huang, and A. Knoll, “Towards autonomous locomotion: CPGbased control of smooth 3D slithering gait transition of a snakelike robot,” Bioinspiration & Biomimetics, vol. 12, no. 3, p. 035001, 2017. View at: Google Scholar
 Z. Bing, L. Cheng, K. Huang, M. Zhou, and A. Knoll, “CPGbased control of smooth transition for body shape and locomotion speed of a snakelike robot,” in Proceedings of the 2017 IEEE International Conference on Robotics and Automation, ICRA 2017, pp. 4146–4153, Singapore, June 2017. View at: Google Scholar
 G. Chen, Z. S. Bing, F. Roehrbein et al., “Efficient braininspired learning with the neuromorphic snakelike robot and the neurorobotic platform,” IEEE Transactions on Cognitive Developmental Systems, vol. 99, pp. 1–1, 2017. View at: Google Scholar
 X. Wu and S. Ma, “Adaptive creeping locomotion of a CPGcontrolled snakelike robot to environment change,” Autonomous Robots, vol. 28, no. 3, pp. 283–294, 2010. View at: Publisher Site  Google Scholar
 M. Tesch, J. Schneider, and H. Choset, “Expensive multiobjective optimization and validation with a robotics application,” in Proceedings of the Neural Information Processing Systems: Workshop on Bayesian Optimization Decision Making, pp. 1–5, 2012. View at: Google Scholar
 S. Hasanzadeh and A. A. Tootoonchi, “Ground adaptive and optimized locomotion of snake robot moving with a novel gait,” Autonomous Robots, vol. 28, no. 4, pp. 457–470, 2010. View at: Publisher Site  Google Scholar
 A. Crespi and A. J. Ijspeert, “Online optimization of swimming and crawling in an amphibious snake robot,” IEEE Transactions on Robotics, vol. 24, no. 1, pp. 75–87, 2008. View at: Publisher Site  Google Scholar
 E. Kelasidi, M. Jesmani, K. Y. Pettersen, and J. T. Gravdahl, “Multiobjective optimization for efficient motion of underwater snake robots,” Artificial Life and Robotics, vol. 21, no. 4, pp. 1–12, 2016. View at: Publisher Site  Google Scholar
 E. Rezapour, K. Y. Pettersen, P. Liljeback, and J. T. Gravdahl, “Path following control of planar snake robots using virtual holonomic constraints,” in Proceedings of the 2013 IEEE International Conference on Robotics and Biomimetics, ROBIO 2013, pp. 530–537, China, December 2013. View at: Google Scholar
 E. Rezapour, K. Y. Pettersen, P. Liljebäck, J. T. Gravdahl, and E. Kelasidi, “Path following control of planar snake robots using virtual holonomic constraints: theory and experiments,” in Proceedings of the IEEE International Conference on Robotics and Biomimetics. in Robotics and Biomimetics, Springer, New York, NY, USA, 2014. View at: Publisher Site  Google Scholar
 A. Mohammadi, E. Rezapour, M. Maggiore, and K. Y. Pettersen, “Maneuvering Control of Planar Snake Robots Using Virtual Holonomic Constraints,” IEEE Transactions on Control Systems Technology, vol. 24, no. 3, pp. 884–899, 2016. View at: Publisher Site  Google Scholar
 Z. Cao, Q. Xiao, R. Huang, and M. Zhou, “Robust neurooptimal control of underactuated snake robots with experience replay,” IEEE Transactions on Neural Networks and Learning Systems, vol. 29, no. 1, pp. 208–217, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 E. Borhaug, A. Pavlov, and K. Pettersen, “Integral LOS control for path following of underactuated marine surface vessels in the presence of constant ocean currents,” in Proceedings of the IEEE Conference on Decision and Control, pp. 4984–4991, 2008. View at: Google Scholar
 W. Caharija, K. Y. Pettersen, A. J. Sørensen, M. Candeloro, and J. T. Gravdahl, “Relative velocity control and integral line of sight for path following of autonomous surface vessels: Merging intuition with theory,” Proceedings of the Institution of Mechanical Engineers, Part M: Journal of Engineering for the Maritime Environment, vol. 228, no. 2, pp. 180–191, 2014. View at: Publisher Site  Google Scholar
 E. Kelasidi, P. Liljebäck, K. Y. Pettersen, and J. T. Gravdahl, “Integral Lineofsight guidance for path following control of underwater snake robots: Theory and experiments,” IEEE Transactions on Robotics, vol. 33, no. 3, pp. 610–628, 2017. View at: Publisher Site  Google Scholar
 D. Zhang, Q. Xiao, Z. Cao, R. Huang, and Y. Fu, “Smooth transition of the CPGbased controller for snakelike robots,” in Proceedings of the 2017 IEEE International Conference on Robotics and Biomimetics (ROBIO), pp. 2716–2721, Macau, December 2017. View at: Publisher Site  Google Scholar
 K. Y. Pettersen and J. T Gravdahl, “Controllability analysis of planar snake robots influenced by viscous ground friction,” in Proceedings of the IEEE International Conference on Intelligent Robots and Systems, pp. 3615–3622, 2009. View at: Google Scholar
 P. Liljebäck, K. Y. Pettersen, Ø. Stavdahl, and J. T. Gravdahl, “A simplified model of planar snake robot locomotion,” in Proceedings of the 23rd IEEE/RSJ 2010 International Conference on Intelligent Robots and Systems, IROS 2010, pp. 2868–2875, Taiwan, October 2010. View at: Google Scholar
 A. J. Ijspeert, “Central pattern generators for locomotion control in animals and robots: a review,” Neural Networks, vol. 21, no. 4, pp. 642–653, 2008. View at: Publisher Site  Google Scholar
 K. Seo, S.J. Chung, and J.J. E. Slotine, “CPGbased control of a turtlelike underwater vehicle,” Autonomous Robots, vol. 28, no. 3, pp. 247–269, 2010. View at: Publisher Site  Google Scholar
 X. S. Yang and S. Deb, “Cuckoo search via levy flights,” Mathematics, pp. 210–214, 2010. View at: Google Scholar
 J. Zhao, S. Liu, M. Zhou, X. Guo, and L. Qi, “An Improved Binary Cuckoo Search Algorithm for Solving Unit Commitment Problems: Methodological Description,” IEEE Access, vol. 6, pp. 43535–43545, 2018. View at: Publisher Site  Google Scholar
 J. Zhao, S. Liu, M. Zhou, X. Guo, and L. Qi, “Modified cuckoo search algorithm to solve economic power dispatch optimization problems,” IEEE/CAA Journal of Automatica Sinica, vol. 5, no. 4, pp. 794–806, 2018. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2019 Zhengcai Cao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.