Study of Obstacle-Crossing and Pitch Control Characteristic of a Novel Jumping Robot

In this study, we demonstrated a novel jumping robot that has the ability of accurate obstacle-crossing jumping and aerial pitch control. The novel robot can quickly leap high into the air with a powerful water jet thruster. The robot was designed to overcome multiple general obstacles via accurate jumping. Then a modified whale optimization algorithm (MWOA) was proposed to determine an optimized jumping trajectory according to the form of obstacles. By comparing with classical intelligent optimization algorithms, the MWOA revealed superiority in convergence rate and precision. Besides, the dynamics model of aerial pitch control was built and its effect was verified by the pitch control experiment. Lastly, the robot’s obstacle-crossing experiments were performed and the results validated the robot’s good ability of obstacle-crossing and aerial body righting. We believe the optimization of trajectory and the pitch control are of great help for the jumping robot’s complex jumping and obstacle-crossing performance.


Introduction
In nature, animals use various locomotion methods to travel across different environments, so that their energy consumptions can be effectively minimized [1]. For instance, the frogs usually apply jumping motion to prey or escape from predators, but sometimes they also crawl in a relatively leisurely situation. Thus, the animal's multiple locomotion abilities can motivate the innovative design of structure and locomotion mechanism of a multi-locomotion robot.
The combination of jumping movement and wheeled movement is common in multilocomotion mode robot, the typical robots include the UMN Scout robot [2], SANDIA wheeled hopping vehicle [3,4], the Boston Dynamics' sand flea robot [5], etc. The jumping movement enables the robot to overcome large-size obstacles and the wheeled movement makes the robot advance rapidly with low energy consumption. Thus, the wheeled jumping robot has better locomotion ability and obstacle-crossing ability compared with other robots.
Inspired by the terrestrial and marine animals, researchers have developed various mobile robots excelling at terrestrial jumping [6][7][8][9] or aquatic jumping [10][11][12][13]. Most terrestrial jumping robots rely on the surficial reactive force to realize jumping actions while the aquatic jumping robots usually accomplish aquatic jumping in virtue of the surface tension of water or the reactive force from the water jetting process. The gas-powered water jet propulsion scheme proposed in [12] is a proper way to achieve immediate takeoff and effective aquatic jumping on account of its high power density [14]. A well-designed water jet thruster is able to convert the internal high-pressure gas into its own kinetic energy, and the generated thrust is continuous in both terrestrial and aquatic environments since it depends on the surficial reactive force. In addition, the propellant water is environmentallyfriendly and easily obtainable. Therefore, the gas-powered water jet propulsion scheme is applicable for terrestrial jumping robots and aquatic jumping robots. The combination of flexible wheels and powerful water jet thrusters on a robot will definitely improve its terrestrial mobility and aerial payload, which is a significative area pending exploration.
Besides wheeled movement and jumping movement, the ability of aerial attitude control (including pitch control and roll control) is also quite important for the jumping robot. The rapid and effective aerial attitude control can increase the jumping robot's possibility of a safe landing. Some small animals, such as lizards [15][16][17], are good at controlling their body's attitude by swinging tails. Inspired by that, several robot prototypes were designed to research the tail's function in body attitude control. Johnson et al. achieved a legged robot's pitch angle control and safe landing by swaying its tail [18]. Chang-siu also experimented with wheeled robot's tail swinging and pitch control [19]. Kohut et al. researched the tail's effect on dynamic turning characteristics for a miniature legged robot [20].
We have designed a series of water jet thrusters which had a good load capacity and was capable of aquatic jumping and terrestrial jumping [21,22]. However, the water jet thruster can neither achieve aerial attitude control nor accurate trajectory control. Therefore, a novel jumping robot was devised with the proposed water jet thruster in [22], which is the first-ever robotic integration of a wheeled system and a water jet thruster. Moreover, a customized optimization algorithm and aerial pitch control model were proposed to optimize the jumping trajectory and help aerial body righting, respectively. We believe the optimization of trajectory and the pitch control are of great help for the jumping robot's complex jumping and obstacle-crossing performance. Specifically, the decision criteria for overcoming general obstacles was firstly proposed to generate the alternative jumping trajectory set. Then a modified whale optimization method aiming at optimizing the jumping trajectory was put forward, and its superiority was verified by comparing it with other different optimization algorithms. Afterwards, the dynamics model of aerial pitch control was built and its reasonability was verified by the jumping robot's free-falling experiment. Lastly, the jumping robot's jumping and obstacle-crossing experiments were carried out, and the experimental results showed that the proposed robot had the good ability of obstacle-crossing and aerial body righting.

Structural Design of the Jumping Robot
The prototype of the proposed jumping robot is shown in Figure 1, and its boundary dimension is 487 × 493 × 452 mm 3 . It mainly comprises the water jet thruster, shell, front legs, hind legs, assistant legs, battery and control electronics. The water jet thruster is used to accomplish terrestrial and aquatic jumping movement for the jumping robot. The front legs and hind legs were made of high-performance nylon which has good toughness. Then the dampers and foam wheels were installed on the front legs and hind legs, respectively. Hence the robot can achieve compliant landing by buffering of its front legs and hind legs. In addition, the driving wheels of the assistant legs and the driven wheels of the hind legs enables the jumping robot to execute sustainable wheeled movement on general road. Another function of the assistant legs is serving as the aerial attitude control device, which aiming at increasing the jumping robot's possibility of safe landing. It is noteworthy that only the robot's pitch angle can be controlled as the assistant legs have just one rotational degree of freedom. Therefore, the details of the jumping robot's components are listed in Table 1.  The jumping robot's control system consists of the control electronics and upper computer software. It is noteworthy that the microcontroller Arduino UNO R3 was applied as the lower computer, and the upper computer can send commands to the upper computer via wireless communication. The jumping robot used 6 channels of PWM signal on the Arduino UNO R3 to control 2 electromagnets (the electromagnets can activate the robot's jumping action and they are connected with electronic switches which also controlled by PWM signals), 2 servos of assistant legs and 2 servos of driving wheel, respectively. Then the attitude sensor MPU6050 module was used to measure the robot's accelerations and attitude angles on 3 axes. In addition, the Bluetooth module was applied to achieve the communication between Arduino and upper computer. Lastly, the working process of the jumping robot's control system are shown in Figure 2. The working process of jumping robot's control system. a set1 and a set2 are two set values of acceleration. e p and e set are the measured value and set value of pitch angle error, respectively, and 0 • is the expected value of pitch angle.

Decision Criteria of Overcoming General Obstacles
When the jumping robot jumps in the air, the pressure energy of high-pressure gas is converted into kinetic energy of the robot and the propellant water [22]. Supposing the bottle plug of the water jet thruster can be released anytime during the pressurization process, and the duration of the water jetting process is quite short, then the maximum pressure of pressurization P max determines the take-off velocity v 0 , and the jumping trajectory can be treated as parabola.
As shown in Figure 3, taking the robot's gravity center o as the original point and the central axis of water jet thruster as the y axis, the carrier coordinate system xoy can be established; taking the initial location of robot's gravity center O as the original point and the vertical direction as Y axis, the ground coordinate system XOY can be established. Then there are two controllable parameters while jumping: the take-off angle ϕ 0 and the take-off velocity v 0 , which are shown in Figure 3a.
It is obvious that different values of ϕ 0 and v 0 can generate different parabolic trajectory. Supposing the parabolic peak's coordinate in ground coordinate system is (x p , y p ), then it can be obtained that: As the parabola goes through original point, the parabolic equation can be written as: In actual situation, the jumping robot may mainly come across three types of obstacles: ditch, wall and platform. Therefore, there are four types of obstacle-crossing tasks: jumping over ditch, jumping over wall, jumping over platform and jumping onto platform ( Figure 3). It is noteworthy that the ditch and wall are special cases of platform, as the ditch and wall can be treated as platforms with little height and width respectively. Therefore, the discussions below are for overcoming the platform.
Supposing the platform's width and height as a and b respectively, the distance between its left side and the origin O is l c , the height of origin O from ground is l r . In consideration of the safety while landing, there should be a safety distance l s (Figure 3c). Therefore, the necessary conditions of jumping over the platform are: the coordinates of the platform's top vertexes should be located below the parabola; the minimum distance between the platform's top vertexes and the parabola should be no more than the safety distance l s . Thus, the decision criteria can be written as: Then the constraint conditions are: In the case of jumping onto the platform, the necessary condition is relatively simple: the coordinates of the platform's top left vertex should be located below the parabola; the minimum distance between the platform's top left vertex and the parabola should be no more than the safety distance l s . Thus, the decision criteria can be written as: Then the constraint conditions are: Hence the value ranges of ϕ 0 and v 0 can be determined according to Equations (3)-(6) in the case of jumping over the platform and jumping onto the platform.

Modified Whale Optimization Algorithm
After determining the value range of the jumping robot's ϕ 0 and v 0 , there will be an obstacle-crossing trajectory set for a particular obstacle. It is necessary to optimize the jumping trajectory according to certain indicators by proper optimization algorithms. The whale optimization algorithm (WOA) was proposed by Mirjalili and Lewis in 2016 [23]. It has good performance in terms of convergence rate and convergence precision. However, it also has the problem of outputting a local optimal solution when dealing with complex problems. Thus, a modified whale optimization algorithm was proposed and applied in the optimization of the jumping trajectory.

Standard Whale Optimization Algorithm
The whale optimization algorithm simulates the foraging behavior of humpback whales. The algorithm can be divided into three stages: encircling prey stage, bubbling-net attacking stage and searching for prey.
1. Encircling prey stage. Supposing the scale of the whale population is N, the dimension of the searching space is d, then the whale individual i can be expressed as X i = (x i 1 , x i 2 , . . . , x i d ), i = 1, 2, . . . , N. The prey's position corresponds to the global optimal solution of objective function. Therefore, the mathematic model of the encircling prey stage can be expressed as: where D is the distance vector between the current optimal solution and the searching individual, t is the current iterations, X (t) is the position vector of whale individual and X p (t) is the current optimal solution. The coefficient vector A and C can be calculated by where r 1 and r 2 are random vectors in [0,1], and a linearly decreases from 2 to 0 with t's growth: 2. Bubbling-net attacking stage. This stage stimulates the process that the whale individual moves spirally around the current prey to approach it. The mathematical model of the spiral motion is: where D' = |X p (t) − X(t)| is the distance vector between the current optimal solution and the searching individual, b is a constant of defining the shape of the logarithmic spiral and l is a random number in [−1,1]. In order to synchronize the behaviors of encircling prey and bubbling-net attack, a probability variable p in [0,1] is introduced, thus the mathematical model can be written as: 3. Searching for prey stage. The whale individual decides encircling for prey or searching for prey by the value of |A|. When |A| > 1, the whale cannot clearly determine the prey's position. Therefore, it is necessary to pick a random individual in place of the prey, that is: 2.3.2. Improvement Strategy 1. Tent mapping strategy. The standard WOA generates an initial population by random method, this may lead to its poor diversity. Therefore, a new chaotic mapping method named tent mapping was applied [24]. For a population with scale N and dimension d, the chaotic sequence y = {y d , d = 1, 2, . . . , d}, y d = {y id , i = 1, 2, . . . , N} can be firstly obtained by tent mapping with the equation Mapping the chaotic sequence to solution space, a population X = {X i , i = 1, 2, . . . , N}, X i = {X id , d = 1, 2, . . . , d} can be obtained, and where X max,d and X min,d are the upper and lower limit of X id respectively. Then the reverse population OX can be obtained according to X by Lastly, merge the population X and OX, and sort the new population according to the objective function's fitness, the top N individuals can be taken as the initial population.
2. Non-linear convergence factor and self-adaptive weight strategy. The linear decrease of convergence factor a in standard WOA cannot coordinate global searching and local exploitation. Therefore, it is better to make the descending speed of a decrease with the progress of iterations t. The proposed equation of non-linear convergence factor is where µ = 0.5 and θ =−π/2 are relevant parameters. Besides, the standard WOA does not take the prey' guiding force diversity for the whale's position update with the progress of iterations into consideration. Therefore, a self-adaptive weight can be introduced in the position update equation to utilize the current optimal solution effectively. The equation of self-adaptive weight is: Then Equations (8) and (12) can be modified as follows: 3. Diversity variation strategy. The population's diversity will gradually decrease with the iterative process in standard WOA, which may lead to the algorithm's premature convergence. Therefore, the variation operation for the population should be conducted when the diversity is below the threshold value. Supposing the fitness value of individual I is f i , f max , f min and f mean are the maximum, minimum and mean fitness values of the current population respectively. Then the variance of current population is: When the variance is below threshold value, the variation operation can be conducted as follows: where r 3 and r 4 are random numbers in [0,1], and X'(t) is a random individual in the population. Therefore, the flow diagram of MWOA is presented in Figure 4.

Performance Test
In order to test the optimization performance of MWOA, 10 benchmark functions were chosen for the performance test. The details of the benchmark functions are listed in Table A1 (Appendix A). It is noteworthy that f 1~f 4 are unimodal functions, f 5~f 8 are multimodal functions and f 9~f 10 are fixed dimension functions, respectively. The performance test was conducted in two aspects: 1. Test the performance improvement effect of each improvement strategy. 2. Compare the convergence rate and precision of MWOA with standard WOA and another two intelligent optimization algorithms, including DA (Dragonfly Algorithm) [25] and HHO (Harris Hawk Optimization) [26]. The population scale and maximum iterations in both performance tests were 30 and 500 respectively; each benchmark function was tested 30 times, and then the mean value and standard deviation were calculated.

Fitness Function
After determining the specific optimization algorithm, the conditions and objective of jumping trajectory optimization should also be determined. Thus, the condition was: the details of jumping robot's location, obstacle's location and size were known; the objective was: to optimize the two take-off parameters (take-off angle ϕ 0 and take-off velocity v 0 ) under the fitness function for a specific obstacle-crossing task.
1. Constraint condition of take-off parameters. The jumping trajectory is determined by the take-off parameters ϕ 0 (rad) and v 0 (m/s), while v 0 was determined by the pressurization pressure, jumping robot's total mass and the mass of propellant water inside the water jet thruster. The mapping relationship between v 0 and these variables can be obtained by the theoretical jumping model in [22]. Then it is noticeable that the water jet thruster's propulsive performance would not be so good if ϕ 0 takes a small value, because the gas-water interface inside the thruster would have a relatively small angle with the thruster's inner wall and the effect of the water-jetting process would be weaken. Thus, according to the actual situation, the constraint condition of ϕ 0 and v 0 are given to be: 2. Fitness function. The fitness function comprises three parts: horizontal jumping distance, vertical jumping distance and landing safety. If the robot confronts a ditch or a flat obstacle, the horizontal jumping distance is its priority. Thus, the fitness function of the horizontal jumping ability is where the denominator is the maximum value of the numerator in the constraint condition.
If the robot confronts a wall or a slim obstacle, the vertical jumping distance is its priority. Thus, the fitness function of the vertical jumping ability is where v 0max is the upper limit of v 0 . If the robot lands on hard ground with high speed, it can be impacted violently and toppled easily. Thus, the fitness function of landing safety is It is noteworthy that these three indicators restrict each other, and they cannot achieve the best value simultaneously. Therefore, each indicator can be weighted and the total fitness function is: where ω i is the weight coefficient and ω 1 + ω 2 + ω 3 = 1.

Dynamics Model of Aerial Pitch Control
The dynamics model for the jumping robot's aerial pitch control belongs to the dynamics of coupled rigid bodies while free falling, which were studied by relevant researchers. Yang et al. built the dynamics model of two rigid bodies connected by a spherical joint and then proposed the controllers by input-output linearization [27]. Agrawal et al. applied differential flatness to achieve the control of two coupled rigid bodies [28]. Chang-Siu et al. put forward the non-linear control model for a robot with a two-link tail and two degreeof-freedom actuations [29].
In order to achieve a reliable and effective aerial pitch control for the jumping robot, the dynamics model was first built. For simplification, the robot was seen as being composed of two parts: the main body and the assistant legs; then the robot only rotates in the pitching direction. Since the water jet thruster's propellant water is drained in a short period of time, the dynamics model assumes that the robotic system does not contain propellant water during the aerial pitch control process. As shown in Figure 5, the assistant leg's centroid is A, its center of rotation on the main body is C, the main body's centroid is B and the whole robot's centroid is O'. Then a dynamic coordinate system O'X'Y' was built where the X' axis and Y' axis are along the horizontal direction and the vertical direction, respectively. Then the robot system's dynamics model can be derived by the Euler-Lagrange method. Due to the establishment of O'X'Y', the robot's translational movement can be decoupled from its rotational movement [30]. As the robot's translational movement is a parabolic movement [31], only the rotational movement is discussed in this section. Therefore, the Lagrangian of the robot system can be written as [32]: where m b and m t are the mass of main body and assistant legs, respectively; I b and I t are the rotational inertia of main body and the assistant legs, respectively; θ b and θ t are the angle of main body and assistant leg relative to X' axis, respectively, and θ m = θ b − θ t is the angle of main body relative to assistant leg. Then applying the Euler-Lagrange method to the robot system, it can be obtained that: It is noteworthy that there is only one external input (the torque τ of assistant leg's servo), so the dynamics Equations can be obtained by substituting (30) into (31): For that dynamics equation, there are two controllable parameters, θ b and θ t , while there is only one input parameter, τ. However, only θ b should be precisely controlled, so the two dynamics equations in (32) can be properly transformed to one which only contains the parameter θ m = θ b − θ t . Firstly, the equations of θ b and θ t can be attained according to (32) as As ..
θ t , thus it can be derived from (33) that ..
Then the parameters . θ b and . θ t can be eliminated by the theorem of angular momentum conservation. Supposing the robotic system's total angular momentum is zero, then the robotic system's angular momentum is Thus, the equations of θ t can be derived as Taking the integral of the first Equation of (36), it can be obtained that where θ b (0) is the initial value of θ b of the aerial pitch control process. Then replacing the integration variable t with θ m , it can be obtained that: where θ m0 and θ mt are θ m 's values at the initial and end moment of the aerial pitch control process respectively. Then combining Equations (32) and (33), the antiderivative of Q/(R + Q) in (38) can be calculated by MATLAB 2016 [33], and the calculation result is Thus, Equation (38) can be transformed to be Therefore, once the initial value θ b (0), final value θ b (t) of θ b and the initial value θ m0 of θ m are determined, the final value θ mt of θ m can be calculated and then the necessary rotation angle of assistant leg ∆θ t = θ t (t) − θ t (0) can also be determined. In most cases, θ b (t) = −π/2, and the ranges of θ b (0) and θ m0 , were given to be [−π,0] and [−π/2, π/2], respectively.

Aerial Pitch Control Experiments of Jumping Robot
Since the dynamics model of aerial pitch control of the jumping robot were proposed, its correctness and reliability should be verified by experiments. The experimental method lets the robot fall freely from the given initial height and records the pitch angle's variation with or without the effect of aerial pitch control movement. The experiments were divided into two groups and each group contained five experiments. There was 580 mL propellant water in the water jet thruster's inner tank in the Group 2, while there was no water in Group 1. Thus, the experimental scheme of aerial pitch control was shown in Table 2. It is noteworthy that there were no aerial pitch control actions in experiment No. 5 and No. 10, which were meant to be the control experiments of experiment No. 4 and No. 9, respectively; besides, the robot's pitch angle θ p = θ b + π/2. In order to prevent obvious error in the robot's roll direction during the free falling process, the robot's initial roll angle was controlled to be zero and its center of gravity in the roll direction was also balanced. The assistant leg's theoretical rotation angle of each experiment can be calculated from Equation (40). For the jumping robot, the related parameters in (40) are: m t = 0.438 kg, m b_1 = 3.922 kg, m b_2 = 4.502 kg, l t = 0.146 m, l b = 0.041 m, I t = 0.00668 kg·m 2 , I b_1 = 0.085 kg·m 2 and I b_2 = 0.087 kg·m 2 . It is noteworthy that m b_1 and m b_2 are the main body's mass in Group 1 and Group 2; I t are the total rotational inertia of the robot's assistant legs as they were meant to rotate simultaneously; I b_1 and I b_2 are the main body's rotational inertia in Group 1 and Group 2. The expected value of the main body and the initial angle of the main body relative to the assistant leg were given to be θ b (t)= −π/2 and θ m0 = 0, respectively. Hence for each experiment in Group 1 and Group 2, the corresponding assistant leg's theoretical rotation angle ∆θ t = θ t (t) − θ t (0) can be calculated according to Equation (40), and the calculation results were listed in Table 2.
In addition, the falling and aerial pitch control process of the jumping robot was recorded by a high-speed camera (i-speed 221, iX Cameras, Britain), and the shooting frame rate was 400 Hz. The variation of the robot's pitch angle θ p was obtained from the analysis of experimental video in the motion analysis software ProAnalyst, and the sampling frequency was 400 Hz; the robot's height from the ground was measured by a telemeter rod.

Jumping Experiments of the Robot 2.6.1. Trajectory Accuracy Verification Experiments
In order to ensure the experimental jumping trajectory's accuracy, it is necessary to carry out the trajectory accuracy verification experiments, which contained four terrestrial jumping experiments of the jumping robot with different take-off parameters (take-off angle ϕ 0 and take-off velocity v 0 ). It is noteworthy that the terrestrial jumping movements were conducted by the water jetting process of the water jet thruster's inner tank. The takeoff velocity v 0 was positively related to the maximum pressurization pressure P max when the robot's total mass and propellant water's mass were determined, and the mapping relationship between v 0 and P max were given by the theoretical jumping model in [22]. Thus, the take-off parameters of four experiments were: a. ϕ 0 = 87 • , P max = 1.81 MPa; b. ϕ 0 = 87 • , P max = 2.23 MPa; c. ϕ 0 = 78 • , P max = 1.81 MPa; d. ϕ 0 = 78 • , P max = 2.23 MPa.
In addition, the volumes of propellant water in four experiments were all 580 mL, and the jumping robot's total mass was 4.36 kg. Thus, the corresponding theoretical take-off velocities with given maximum pressurization pressures of 1.81 MPa and 2.23 MPa can be calculated to be 5.36 m/s and 6.46 m/s, respectively.
In the robot's terrestrial jumping experiment, the propellant water was firstly injected into the water jet thruster. Then the robot's initial attitude was adjusted according to the given take-off angle. The robot started jumping movement a while after the liquid nitrogen was injected and the water jet thruster's inner tank was sealed. The jumping process was recorded by a digital camera (FDR-AX45, Sony, Tokyo) and the experimental jumping trajectory of robot's centroid was extracted from the experimental video by ProAnalyst.

Obstacle-Crossing Experiments
After the accuracy of robot's experimental jumping trajectory was verified, the obstaclecrossing experiments can be conducted with the third optimized take-off parameters whose priority was landing safety. These experiments were meant to verify the robot's abilities of obstacle-crossing and aerial pitch control. Similarly, the volume of propellant water inside the water jet thruster's inner tank was 580 mL, and the jumping robot's total mass was 4.36 kg. As the optimized take-off parameters of the experiments of jumping onto the platform and jumping over the platform were both ϕ 0 = 76.0 • and v 0 = 5.57 m/s, then the corresponding maximum pressurization pressure of the water jet thruster was 1.89 MPa.
As shown in Figure 3c,d, the experimental obstacle-crossing parameters were in accordance with those of the optimization process. Thus, for the experiment of jumping onto the platform, l c = 0.615 m and b = 1.04 m; and for the experiment of jumping over the platform, l c = 0.57 m, a = 0.105 m and b = 0.67 m. The experimental procedure was the same with that of the trajectory accuracy verification experiment. In addition, the robot's experimental jumping trajectories and pitch angle curves were obtained by analysis of experimental video in ProAnalyst and the onboard attitude sensor respectively.

Performance Improvement Effect of Each Improvement Strategy
In order to analyze the performance improvement effect of each improvement strategy, the 10 benchmark functions were used to test each improvement strategy. Therefore, the WOA modified by tent mapping strategy was named MWOA-I; the WOA modified by non-linear convergence factor and self-adaptive weight strategy was named MWOA-II; the WOA modified by diversity variation strategy was named MWOA-III. The comparison results were shown in Table 3. The bold column for each benchmark function in Table 3 is the best result. According to Table 3, the MWOA had best performances in 8 out of the 10 benchmark functions, then the MWOA-II and MWOA-III had 6 best results and 4 best results respectively. Therefore, it can be concluded that the non-linear convergence factor and self-adaptive weight strategy had a larger performance improvement effect than the other two strategies. However, the MWOA, which combined the three strategies still outnumber any other algorithms which utilized only one strategy in the best optimization results.

Performance Comparison with Other Intelligent Optimization Algorithms
In order to compare the MWOA's performance with other intelligent optimization algorithms, the 10 benchmark functions were used again to test the performance of WOA, DA, HHO and MWOA. The comparison results were shown in Table 4, and the convergence curves of these four algorithms for f 1 and f 2 were shown Figure 6 (the other convergence curves were presented in Figure A1). Similarly, the bold column for each benchmark function in Table 4 is the best result. According to Table 4, the MWOA had the 7 best results among the 10 benchmark functions, then the HHO and DA had the 5 best results and 1 of the best results, respectively. Thus, it can be concluded that the convergence precision and stability of the MWOA were the best according to the evaluation indicators of the mean value and standard deviation.
By analyzing the convergence curves of four intelligent optimization algorithms, it can be concluded that the MWOA had the fastest convergence speeds in 6 benchmark functions (f 1 , f 2 , f 3 , f 6 , f 7 and f 8 ), which surpassed the other algorithms. Compared with standard WOA, the MWOA had better convergence speeds in all 10 benchmark functions. For the optimization of f 6 , f 7 , f 8 and f 10 , the convergence curves of MWOA showed multiple inflection points while descending, which means its ability of escaping from local optimal solution was obviously improved. Thus, by performance comparison with other intelligent optimization algorithms, the reasonability and effectiveness of MWOA can be proved.

Optimization of Jumping Trajectory
After determining the total fitness function, the centroid jumping trajectory of jumping onto a platform and jumping over a platform can be optimized under different weight coefficient combinations. As shown in Figure 3c Therefore, the searching space was firstly determined by the decision criteria (3) and (5). Then the weight coefficient combinations were determined and the fitness function (29) was optimized by the MWOA. Lastly, the optimization results were listed in Table 5, where optimization I and II were for jumping onto a platform and jumping over a platform, respectively. The corresponding centroid jumping trajectories were shown in Figure 7.   Table 5 for jumping onto the platform and jumping over the platform, respectively. It can be seen that the optimum take-off parameters for jumping onto the platform and jumping over the platform are both ϕ 0 = 76.0 • , v 0 = 5.57 m/s in their separate searching spaces, thus the reasonability of optimization results of theoretical jumping trajectory can be verified.    Figure 10 shows the jumping robot's pitch angle curves of the aerial pitch control experiment, each curve started at the moment when the robot began to fall and ended at the moment when it contacted the ground. It can be seen from Figure 9 that the attitude sensor module detected the robot's acceleration variation when the robot started falling freely, thus the assistant legs were driven to rotate a theoretical angle ∆θ t = −48.5 • . During the rotation of the assistant legs, the main body's pitch angle also increased from the initial value −17.3 • . According to Figure 10a, the robot's pitch angle reached 0 • at 422 ms and then 1.81 • at 531 ms when the robot contacted the ground. After landing, the assistant legs reset to their initial position (parallel to the central axis of the water jet thruster), and the aerial pitch control process was completed.  As the expected value of pitch angle is 0 • , it can be known from Figure 10 5 and RMSE 9 < RMSE 10 , then the effectiveness of the mathematical model of aerial pitch control movement can be verified and it can be applied in the actual jumping and obstacle-crossing experiment for the jumping robot.

Trajectory Accuracy Verification Experiments
After the outdoor terrestrial jumping experiment, the four experimental jumping trajectories can be obtained and their comparisons with corresponding theoretical jumping trajectories were shown in Figure 11. Therefore, the RMSE of each experimental jumping trajectory can be calculated to evaluate its accuracy. The calculation results of RMSE for Figure 11a-d were: RMSE a = 0.057 m, RMSE b = 0.095 m, RMSE c = 0.074, RMSE d = 0.071 m, respectively. The maximum average degree of deviation in vertical displacement was 0.095 m, thus it can be concluded that the error between experimental jumping trajectory and theoretical jumping trajectory was small.

Obstacle-Crossing Experiments
The time sequence of jumping robot's process of jumping onto platform and jumping over platform are shown in Figures 12 and 13 respectively (The whole obstacle-crossing jumping processes are presented in https://www.youtube.com/watch?v=op6UnGVcYnw, (accessed on 30 March 2021) and in Video S1 (Supplementary Materials)). It can be known that the robot did successfully jump over and jump onto the platform with given take-off parameters. Although the robot did not compensate well in its roll direction due to the lack of attitude control action in the roll direction, the assistant legs' rotation did have an obvious adjusting effect on the robot main body's pitch angle. Then the jumping robot's jumping trajectory comparisons and pitch angle curves of the main body are shown in Figure 14.  The comparison of the jumping robot's centroid trajectory of the experiment of jumping onto the platform and jumping over the platform are shown in Figure 14a,b. It can be calculated that the RMSE values of these two experiments are RMSE a = 0.0537 m and RMSE b = 0.0731 m, respectively, which are both within the reasonable range. Thus, the accuracies of the robot's jumping trajectory are acceptable.
According to Figure 14c, the robot's pitch angle was stable at around −14 • in the initial stage. Then the main body's pitch angle started approaching the expected value of 0 • in the effect of the assistant legs' rotation after t = 550 ms. Lastly, the main body's pitch angle increased to −2.45 • at t = 890 ms when the robot landed on the platform. Thus, the duration and the absolute value of error of aerial pitch control was 340 ms and 2.45 • , respectively. Similarly, in Figure 14d, the robot's pitch angle was also stable at around −14 • initially. Then the main body's pitch angle started approaching the expected value of 0 • after t = 290 ms. However, the pitch angle did not change obviously after t = 820 ms, and the final pitch angle was −2.79 • at t = 1100 ms when the robot landed on the platform. Thus, the duration and the absolute value of error of aerial pitch control was 530 ms and 2.79 • respectively. Additionally, it can be calculated from Figure 14a that the initial experimental take-off angle and take-off velocity were ϕ 0e1 = 75.54 • and v 0e1 = 5.15 m/s respectively for the experiment of jumping onto the platform, and therefore the corresponding fitness value was f e1 = 0.8328. Similarly, the calculation results for the experiment of jumping over the platform were ϕ 0e2 = 75.53 • and v 0e2 = 5.24 m/s, respectively, therefore the corresponding fitness value was f e2 = 0.8337. The absolute value of relative error of fitness value for the experiment of jumping onto the platform and jumping over the platform were e f1 = 0.23%, e f2 = 0.12%. Therefore, the aerial pitch control model's reasonability and the jumping robot's accurate obstacle-crossing ability can be verified by these two obstacle-crossing experiments.

Conclusions
In this paper, a novel jumping robot was devised and the process of how it completed actual obstacle-crossing tasks and aerial pitch control movements are presented. Firstly, the decision criteria for overcoming general obstacles were proposed to obtain the feasible jumping trajectory set for the jumping robot. In order to optimize the jumping trajectory according to the fitness function value, a modified whale optimization algorithm (MWOA) was put forward and its performance superiority was verified by comparing it with other different intelligent optimization algorithms. Then the dynamics model of aerial pitch control was established and the jumping robot's aerial pitch control experiment was conducted. The experimental results demonstrated the effectiveness of the proposed dynamics model. In addition, the robot's trajectory accuracy verification experiment showed that the robot's jumping trajectory was accurate and reliable. Finally, the robot's obstaclecrossing experiments were conducted and the experimental results verified the robot's good ability of obstacle-crossing and aerial pitch control. Collectively, these achievements provide meaningful guidance for the jumping robot's performance in dealing with complex obstacle-crossing tasks and future research on jumping robot's reliability as well as environmental adaptability.