Land–Air–Wall Cross-Domain Robot Based on Gecko Landing Bionic Behavior: System Design, Modeling, and Experiment

: Based on the bionic behavior of geckos, this paper presents a land–air–wall cross-domain robot which can ﬂy in air, run on the ground, and adhere to various wall surfaces. When geckos jump and adsorb to vertical surfaces such as trunks, they can still adsorb to the wall with a large contact speed. Inspired by this phenomenon, we analyze the mechanism, apply it to our robot, and optimize the design of the robot structure. In addition, geckos use their tails to adjust posture to achieve abdominal landing during the process of falling. Inspired by this phenomenon, based on the rotor lift/power curve, we optimize the center of gravity by controlling the servo angle. The initial center of gravity offset of the robot is estimated by the extended state observer. The method reduces the distance between the center of gravity and the geometric center, balances the load of each propeller, and ﬁnally reduces the total power. The experiment and simulation results validate the feasibility of the land–air–wall cross-domain robot and the bionic methods.


Introduction
In recent years, land-air amphibious robots have become a hot research issue.Different countries have developed land-air amphibious robots with varied morphology.Land-air amphibious robots integrate the specialty of both flying robots and ground robots, including rapid space maneuverability and long-time working capability.At the same time, robots that are able to work in special environments have also attracted widespread attention.The flying adsorption amphibious robots have the ability to work in air and can perch on the surface of objects for a long time [1].However, land-air amphibious robots cannot work on vertical surfaces for a long period of time, and the flying adsorption amphibious robots do not have ground motion capability.

Motivation
In order to combine the advantages of the land-air amphibious robot and the flying adsorption amphibious robot, this paper develops a land-air-wall cross-domain robot (LAWCDR) that can be widely used in disaster rescuing and battlefield reconnaissance environment monitoring.
By observing the behavior of geckos, the way they jump and attach to vertical surfaces caught our attention.When the front claws of the gecko come in contact with the vertical wall, the front claws are bounced off due to the excessive speed, and then the rear claws and tail remain attached to the wall.The geckos are able to adsorb to the surface with rear claws, tail, and spine.The mechanism of this phenomenon is applied to our robot (LAWCDR), in which two front wheels simulate the rear claws of geckos, and the suction cup and support bars simulate the head and front claws.After the wheels make contact with the wall, the suction cup (head) and support bars (front claws) are moved forward to Appl.Sci.2022, 12, 3988.https://doi.org/10.3390/app12083988https://www.mdpi.com/journal/applscimake contact with the wall through the spine structure simulated by the servo to realize adsorption.The shock absorber of the wheels extends the dwelling time on the wall.Unlike real geckos, since our front wheels (rear claws) do not have a suction module, the front wheels may come off the wall after the suction cup (head) and support bars (front claws) come into contact with the wall.
In addition, the non-coincidence of the center of gravity and geometric center of the robot (LAWCDR) requires greater hovering power consumption.By observing the phenomenon that the belly landing is achieved by the gecko adjusting the attitude during the aerial motion, we adjust the body center of gravity in the air by adjusting the servos angle and reduce the aerial power by reducing the center of gravity offset.

Related Work
There are many different types of land-air amphibious robots.In 2009, the Hybird Land-Air Robot, researched by the University of Minnesota, folded the rotor of the vehicle in the ground state with the area between the two wheels, driven by ground motors.When flying, the robot stood up and unfolded to form a co-axial twin-rotor helicopter [2][3][4].In 2013, the HyTAQ Robot was proposed by the Robotics Laboratory of the Illinois Institute of Technology.By installing an external roller cover that acts as a rotor shield and as a roller for ground motion, the ground part was powered by wind propulsion [5].The robot called F-Star developed by Ben-Gurion University [6] had an adjustable angle of the arm while traveling on the ground, allowing F-Star to effectively raise or lower its fuselage when encountering an obstacle.
Flying adsorption amphibious robots can perch on the surface of objects in different ways, such as negative pressure adsorption [7,8], dry adsorption [9,10], and electrostatic adsorption.The collision of flying adsorption amphibious robots with walls and recovery from their collision have likewise received extensive attention.For example, Ref [11] used a compliant arm design that allowed it to fly freely while allowing a passive degree of freedom to absorb the impact.Additionally, Ref [12] installed a protective buffer around the rotor and modeled the collision of the airframe with a vertical wall.The simulation and experimental preliminary results verified the possibility and conditions for recovery from the collision.In [13], a collision recovery method was proposed for a quadrotor based on a protective cushion.It was validated by a comprehensive Monte Carlo collision bumper simulation, showing the feasibility of recovery from challenging collision scenarios.Research in [14] investigated the control problem of recovering a quadrotor after a collision with a rod object.In [15], the authors proposed a fast IMU-based collision detection method and potential field-based theory for collision recovery.
Designing flying adsorption amphibious robots through bionics is another effective method.Researchers have observed the existing roosting mechanisms in nature and applied them to robot design.For example, in [16][17][18] based on bat/gecko mimics, flight adsorption amphibious robots attached adsorption devices to multi-rotor vehicles for flight adsorption.In [19][20][21], the flight adsorption method based on bird claws was introduced.They used mechanical claws to mimic the claws of birds and made them capable of perching on various tree trunks.Furthermore, [22] used hooked claws and fixed-wing vehicles to achieve perching by mimicking the landing process of birds.

Contribution
In this paper, we mimic the behavior of the gecko and develop a land-air-wall crossdomain robot.The robot is capable of aerial flight, ground movement, and wall perching.Our robots have better adaptability compared to previous robots.For example, in disaster relief scenarios, the robot can maneuver quickly by flying, reconnoiter on the ground for a long time, or perch on the wall for surveillance and communication relay missions.
The specific contributions are as follows.

1.
By imitating the behavior of the gecko touching the wall, the robot can control the two joint movements for adsorption when the front wheels contact the wall.This mimics the gecko controlling the spine and tendons to achieve the process of adsorption.
Compared with previous flying adsorption amphibious robots, our bionic adsorption method reduces the requirement for the position feedback so that it is more suitable for practical scenarios.This is achieved by giving a certain pitch angle to first make contact between the wheel and the wall, and then controlling the servos to perform the adsorption.Position or force feedback are not necessary for the perching process.

2.
We design the mechanical structure to protect the rotor from contact with the wall and crashing.The optimized design of the adsorption device is carried out through the force analysis of the robot during the adsorption process.

3.
We mimic the behavior of the gecko in regulating its attitude during landing.The initial center of gravity offset of the robot is observed by the extended state observer.
Based on the lift/power curve, we balance the load of the four propellers by adjusting the joint angle to reduce the hovering power and improve the endurance.

Outline
This paper is organized as follows.In Section 2, we introduce the overall structural design and optimization of the land-air-wall cross-domain robot (LAWCDR).In Section 3, we present the dynamics modeling and control of the robot during different motions.Simulation and experimental results are given in Section 4 to verify the effectiveness of the robot (LAWCDR).The paper ends with conclusions in Section 5.

Overall Structural Design of Robot (LAWCDR)
The overall design of the land-air-wall cross-domain is shown in Figure 1, which mainly includes the ground module, the flight module, and the adsorption module.The ground module consists of the front wheels with shock absorbers and an omni-directional rear wheel.The robot turns by the differential speed of the two front wheels.The shock absorber is mainly responsible for providing a certain cushioning force and correcting heading when the robot contacts the wall during the adsorption process.The shock absorber does not work when it is on the ground, and the spring stretches to its maximum stroke under gravity.
The specific contributions are as follows.
1.By imitating the behavior of the gecko touching the wall, the robot can control the two joint movements for adsorption when the front wheels contact the wall.This mimics the gecko controlling the spine and tendons to achieve the process of adsorption.Compared with previous flying adsorption amphibious robots, our bionic adsorption method reduces the requirement for the position feedback so that it is more suitable for practical scenarios.This is achieved by giving a certain pitch angle to first make contact between the wheel and the wall, and then controlling the servos to perform the adsorption.Position or force feedback are not necessary for the perching process.2. We design the mechanical structure to protect the rotor from contact with the wall and crashing.The optimized design of the adsorption device is carried out through the force analysis of the robot during the adsorption process.3. We mimic the behavior of the gecko in regulating its attitude during landing.The initial center of gravity offset of the robot is observed by the extended state observer.
Based on the lift/power curve, we balance the load of the four propellers by adjusting the joint angle to reduce the hovering power and improve the endurance.

Outline
This paper is organized as follows.In Section 2, we introduce the overall structural design and optimization of the land-air-wall cross-domain robot (LAWCDR).In Section 3, we present the dynamics modeling and control of the robot during different motions.Simulation and experimental results are given in Section 4 to verify the effectiveness of the robot (LAWCDR).The paper ends with conclusions in Section 5.

Overall Structural Design of Robot (LAWCDR)
The overall design of the land-air-wall cross-domain is shown in Figure 1, which mainly includes the ground module, the flight module, and the adsorption module.The ground module consists of the front wheels with shock absorbers and an omni-directional rear wheel.The robot turns by the differential speed of the two front wheels.The shock absorber is mainly responsible for providing a certain cushioning force and correcting heading when the robot contacts the wall during the adsorption process.The shock absorber does not work when it is on the ground, and the spring stretches to its maximum stroke under gravity.The flight module mainly includes four rotor motors and the corresponding propellers, which are responsible for flying in three-dimensional space.The adsorption module consists of a suction cup, a vacuum pump, support bars, carbon tubes and silicone tubes connecting the various parts, and two RS485 servos responsible for the movement of the suction cup.While perching on the wall, the support bars prevent the suction cup from peeling off the wall due to the rotational torque generated by gravity.In addition, the wheel and the support bars work together to prevent the propellers from contacting the wall.
The structure and coordinates of the robot are shown in Figure 2, and the inertial coordinate system is represented by O w , O b , O h , and O c respectively, where O w denotes the world coordinate system and O b denotes the body coordinate system.T i denotes the lift generated by the ith rotor.The angular bisector of the second and third motors is taken as the y-axis, and the angular bisector of the third and fourth motors is taken as the x-axis.O h and O c denote the coordinate systems represented by the first and second joints, and O c coincides with the suction cup end coordinate system.The x-axis distance from the rotor to the geometric center of the body is l, the distance from the first joint to the second joint is d 2 , and the first joint is fixedly connected to the robot.The distance between the second joint and the suction cup is d 1 .h denotes the z-axis distance of O h of the body coordinate system.
suction cup from peeling off the wall due to the rotational torque generated by gra addition, the wheel and the support bars work together to prevent the propeller contacting the wall.
The structure and coordinates of the robot are shown in Figure 2 x-axis distance from the rotor to the geometric center of the body is l , the distanc the first joint to the second joint is 2 d , and the first joint is fixedly connected to the The distance between the second joint and the suction cup is  As shown in Figure 3, the robot can generally be divided into four working (a) Flight state: when the servo is adjusted to the best position, the power consump reduced according to the mass distribution.(b) Preparation state: as a preparator for adsorption, the suction cup is adjusted to appropriate position to reduce the reaction time when the wheels are in contact with the wall.(c) Adsorption state gle-line LIDAR is used to detect the relative position and attitude between the rob the wall.The head and front claws of the robot are controlled to face the wall.Af robot detects its contact with the wall, the robot controls the servo to reach forw achieve adsorption.(d) Roost state: at this time, the propellers and servos are stop save power and reduce heat.The robot adsorbs on the wall under the action of t uum pump and support bars.Because the body is in a free state, the front whe compressed, but they still prevent propellers from contacting the wall.As shown in Figure 3, the robot can generally be divided into four working states.(a) Flight state: when the servo is adjusted to the best position, the power consumption is reduced according to the mass distribution.(b) Preparation state: as a preparatory state for adsorption, the suction cup is adjusted to appropriate position to reduce the robot's reaction time when the wheels are in contact with the wall.(c) Adsorption state: a single-line LIDAR is used to detect the relative position and attitude between the robot and the wall.The head and front claws of the robot are controlled to face the wall.After the robot detects its contact with the wall, the robot controls the servo to reach forward to achieve adsorption.(d) Roost state: at this time, the propellers and servos are stopped to save power and reduce heat.The robot adsorbs on the wall under the action of the vacuum pump and support bars.Because the body is in a free state, the front wheels are compressed, but they still prevent propellers from contacting the wall.

Design and Analysis of the Adsorption of Bionic Structures
The robot (LAWCDR) mainly mimics the process of a gecko pe surface even after the collision, rather than being bounced off.In [23], slow motion of an Asian flat-tailed gecko (Hemidactylus platyurus) j to trunk, which allows the gecko to grasp the trunk through the sticky paws when it is bounced back from a vertical surface, with the tail pro the gecko to counteract the moment of flipping backward and again paws to contact with the wall.The researchers designed a simulat with four Velcro feet and a tail to control the body flexion of the gec tendons.Experiments have shown that the simulated gecko with the crease the adsorption success rate from 15% to 55%.In nature, the increase its landing success rate from nearly 0% to 87%.
We mimic the spine of a gecko by shock absorbers and a two When the wheels come into contact with the wall, the moment cause and affects the forward reaching of the suction cup.The front claws in order to improve the success rate of the adsorption.The wheels (rea wall and provide guidance of the heading before the front claw stret contact with the wall.As shown in Figure 4

Design and Analysis of the Adsorption of Bionic Structures
The robot (LAWCDR) mainly mimics the process of a gecko perching on a vertical surface even after the collision, rather than being bounced off.In [23], Jusufi analyzed the slow motion of an Asian flat-tailed gecko (Hemidactylus platyurus) jumping from trunk to trunk, which allows the gecko to grasp the trunk through the sticky toe pads of its hind paws when it is bounced back from a vertical surface, with the tail providing leverage for the gecko to counteract the moment of flipping backward and again allowing the front paws to contact with the wall.The researchers designed a simulated gecko equipped with four Velcro feet and a tail to control the body flexion of the gecko through internal tendons.Experiments have shown that the simulated gecko with the tail was able to increase the adsorption success rate from 15% to 55%.In nature, the tail helps the gecko increase its landing success rate from nearly 0% to 87%.
We mimic the spine of a gecko by shock absorbers and a two servos mechanism.When the wheels come into contact with the wall, the moment causes the body to rotate and affects the forward reaching of the suction cup.The front claws keep stretching out in order to improve the success rate of the adsorption.The wheels (rear claws) stay on the wall and provide guidance of the heading before the front claw stretches out and makes contact with the wall.As shown in Figure 4 where M  and M  are the y-axis and z-axis moment.with the wall in the body coordinate system, and 3 l is the y-axis distance from the two front wheels to the center of the body.i  is the rotational speed of the ith rotor.Take the left wheel as an example.
where M θ and M ψ are the y-axis and z-axis moment.M G is the gravitational moment.F L and F R are the equivalent contact forces of the left and right wheel on the wall.d wl and d wr are the z-axis distances of the left and right wheel contact points with the wall in the body coordinate system, and l 3 is the y-axis distance from the two front wheels to the center of the body.ω i is the rotational speed of the ith rotor.Take the left wheel as an example.
where I w is the rotational inertia of the wheel assembly.k s is the coefficient of elasticity of the spring.r w is the rotation radius of the spring assembly.r q is the rotation radius of the wheel assembly.h 1 and h 2 are the x-axis distance and z-axis distance of the spring mount point in the body coordinate system, respectively.σ L is the angle between the carbon tube of the wheel assembly and the z-axis in the body coordinate system.h p is the height of the propeller in the body coordinate system.R is the radius of the propeller.r is the radius of the wheel.The length of the spring, h s , is represented in the following equation.

cos( ) sin
sin( ) where w I is the rotational inertia of the wheel assembly. is the angle between the carbon tube of the wheel assembly and the z-axis in the body coordinate system.p h is the height of the propeller in the body coordinate system.R is the radius of the propeller.r is the radius of the wheel.The length of the spring, s h , is repre- sented in the following equation.
It must be ensured that the propeller cannot exceed the plane formed by the two front wheels and the suction cup when the robot is in contact with the wall.The coordinate of the suction cup in the world coordinate system is expressed as follows: The coordinate of the contact point between the wheels and the wall is expressed as follows: It must be ensured that the propeller cannot exceed the plane formed by the two front wheels and the suction cup when the robot is in contact with the wall.The coordinate of the suction cup in the world coordinate system is expressed as follows: The coordinate of the contact point between the wheels and the wall is expressed as follows: The coordinate of the point on the propeller closest to the wall is expressed as follows: In order to protect the propeller, when the suction cup and the wheels are in contact with the wall, the propellers are prevented from obtaining contact with the wall.In order to simplify the model, when the wheel is at the maximum compression, p prop and p b are simultaneously on the same side of the line formed by the suction cup and the wheel.
In the x-z plane, the equation of the line passing through p wheel and p prop is shown in the following equation: Therefore, The constraint of η 1 and η 2 can be calculated by Equations ( 6)-( 8) and (10).
The force of the robot is shown in Figure 6a.If the wheel is in a free state, the friction between the wall and the front wheels belongs to rolling friction, which is much smaller than the sliding friction (ignoring the rolling friction).f 1 denotes the friction between the suction cup and the wall, and f 2 denotes the friction between the support bars and the wall.N 1 and N 2 are the support forces of the wall for the suction cup and the support bars, respectively.T and G denote the total lift and gravity force, respectively.
The force of the robot is shown in Figure 6a.If the wheel is in a free state, the friction between the wall and the front wheels belongs to rolling friction, which is much smaller than the sliding friction (ignoring the rolling friction).As shown in Figure 6b, after adsorption succeeds, the robot loses the lift force T generated by the propeller.Therefore, we need the frictional forces f 1 and f 2 provided by the wall to balance with the gravity force G.At this time, the two servos can be in a free state to avoid the heat generated by the long working time.
In order to increase the weight capacity of the adsorption module, the structure of the suction cup and the support bars is shown in Figure 7.The angle between the carbon fiber tube which is connected to the suction cup and the support bars is θ s .The angle between the carbon fiber tube attached to the suction cup and the horizontal plane is α.Ignoring the mass of the support bars and the carbon fiber tube, the forces in the horizontal and vertical directions are balanced.As shown in Figure 6b, after adsorption succeeds, the robot loses the lift force T generated by the propeller.Therefore, we need the frictional forces 1 f and 2 f pro- vided by the wall to balance with the gravity force G .At this time, the two servos can be in a free state to avoid the heat generated by the long working time.
In order to increase the weight capacity of the adsorption module, the structure of the suction cup and the support bars is shown in Figure 7.The angle between the carbon fiber tube which is connected to the suction cup and the support bars is s θ .The angle between the carbon fiber tube attached to the suction cup and the horizontal plane is  .
Ignoring the mass of the support bars and the carbon fiber tube, the forces in the horizontal and vertical directions are balanced.Additionally, where µ 1 is the coefficient of friction between the suction cup and the wall, and µ 2 is the coefficient of friction between the support bars and the wall.d L is the distance from the connection point of the support bars to the suction cup.m is the mass of the robot and g denotes the acceleration of gravity Because suction cup and support bars are not the same material, µ 1 and µ 2 are not consistent.In general, µ 1 < µ 2 , therefore it is necessary to raise N 2 to increase the weight that the robot can hold.
The influence of N 2 on α is shown in Figure 8. N 2 becomes larger as α increases.As α becomes larger, the robot assigns more F suction to N 2 .Assuming d 1 = 0.35 m,d L = 0.2 m, N 2 is approximately equal to 1.996 × mg when α = 10 • which is 1.14 times as large as N 2 = 1.75 × mg when α = 0 • .The increase of α may cause the suction cup to be unable to fit against the wall and form a negative pressure environment.Therefore, we increased the contact angle of the suction cup by adding a pivot and a limit block, as shown in Figure 7.The suction cup is connected to the robot by a pivot and is limited by a limit block.The suction cup and the carbon fiber tube are connected by a flexible silicone tube.Therefore, the suction cup can rotate [−5 • , 5 • ] on the limit block, and because the double-layer suction cup itself has [−6.75 • , 6.75 • ] angle range, the robot's suction cup has [−11.75  = .The increase of  may cause the suction cup to be unable to fit against the wall and form a negative pressure environment.Therefore, we increased the contact angle of the suction cup by adding a pivot and a limit block, as shown in Figure 7.The suction cup is connected to the robot by a pivot and is limited by a limit block.The suction cup and the carbon fiber tube are connected by a flexible silicone tube.Therefore, the suction cup can rotate [ 5 ,5 ] −   on the limit block, and because the double-layer suction cup itself has [ 6.75 , 6.Then, the maximum friction force that the robot can provide is about 69.37 N, which is two times greater than the gravity.As shown in Figure 9, we tested the variation of N 2 for different clamping angles of the support bars θ s .As the clamping angle θ s becomes larger, N 2 gradually becomes smaller.However, N 1 cannot be lower than a certain threshold, otherwise it will lead to an increase in the leakage rate of the suction cup and further lead to the disengagement of the suction cup.Therefore, we take θ s as 45 degrees.When the suction cup works on the surface of the painted wall, the pressure difference is 40 kpa, and the diameter suction cup is 0.03 m.F suction is about 110 N. The total mass of the robot is 3 kg.Then, the maximum friction force that the robot can provide is about 69.37 N, which is two times greater than the gravity.

Air Dynamics Modeling and Ground Dynamics Modeling
The position of the robot in the world coordinate system is The dynamics model of the robot can be represented using the Euler-Lagrange equations.() Mq is the inertia matrix, ( , ) C q q is the Coriolis matrix, () Gq is the gravity vector,  is the generalized moment input, and B is the input transformation matrix.


The inertia matrix () Mq is shown as follows [24]:

Air Dynamics Modeling and Ground Dynamics Modeling
The position of the robot in the world coordinate system is p w = (x, y, z) T and the attitude angle is Ω = (ϕ, θ, ψ) T .The joint angle of the two servos is η = (η 1 , η 2 ) T .The rotation matrix R w b between O w and O b is denoted as: where cθ denotes cos(θ) and sθ denotes sin(θ).cϕ denotes cos(ϕ) and sϕ denotes sin(ϕ).cψ denotes cos(ψ) and sψ denotes sin(ψ).
The dynamics model of the robot can be represented using the Euler-Lagrange equations.
M(q) is the inertia matrix, C(q, .q) is the Coriolis matrix, G(q) is the gravity vector, τ is the generalized moment input, and B is the input transformation matrix.The inertia matrix M(q) is shown as follows [24]: where q are the matrices that form a velocity relationship.
. p and ω are the translational and angular velocity of body in the inertial frame.The Corillis matrix C(q, .q) is shown as follows [24,25]: where c kj indicates the element of C(q, . q) at the kth row and jth column.m ab indicates the element of M(q) at the ath row and bth column.
The gravity vector G(q) is shown as follows: where e g is the unit vector [ 0 0 1 T .p b i denotes the position of the ith joint in the body coordinate system.U is the potential energy of the robot.
Although the Euler-Lagrange Equation ( 16) can describe the motion of the robot, the accumulated dynamic model is too complicated to design an efficient controller.Because the range of motion of the joints during adsorption is limited, the joint firstly needs to be adjusted to the preparation state and the mass of the joints is small compared to the rest of the robot.Thus, we decouple the joint and the robot body, use the motion at the end of the suction cup as a perturbation to the robot body, and improve the efficiency of the control.
The rotor lift is expressed by the following equation.
where T i is the ith rotor lift, k f is the lift coefficient, ω i is the rotational speed of the ith rotor.
T tot is the total lift of the robot.The dynamics model of robot is shown in the following equation, ignoring air resistance.
where A = η θ −η θ − θ and d 2 sin η θ is the relative forward extension of the suction cup.
m is the mass of the body, m s is the mass of the first joint, and m v is the mass of second joint.
In the actual robot adsorption process, the roll angle is always close to 0 degrees in order to maintain the robot balance.Therefore, M can be expressed as follows.
where k m is the drag coefficient.τ ext is the moment due to the external force F ext .L 1 is the y-axis moment of the body mass in the body coordinate system.
When the suction cup is not in stable contact with the wall, the suction cup makes contact with the wall at small speeds modeled by a spring damping model.
where x p is the distance between O b and the wall, K is the contact stiffness, and C is the damping coefficient.
The velocity v 0 and angular velocity w 0 of the ground motion of the robot are shown in the following equation.
where v L and v R are the left and right wheel speeds, respectively.According to Equation ( 26), we can obtain the position and heading angle in the world coordinate system.
where (x 0 , y 0 ) is the position of the robot on the ground and ψ 0 is the heading angle for the ground motion.Similarly, the dynamics modeling of the ground part can be represented by the Euler-Lagrange equation.

Flight Bionic Control
According to [26], scientists discovered that geckos are able to adjust their airborne posture through their tails in the air, so that they can land on their bellies, just like a flying wingsuit athlete adjusts its center of gravity and drag to land successfully.In the same motor and propeller configuration, the faster the speed is, the higher the lift is, but the less efficient it is.Inspired by this, when robots fly in the air without adsorption, we also need to adjust the servos angle to balance the load of the four rotors, improve control performance, reduce power consumption, and extend the working time.Assume that W is the power that a single rotor requires to generate lift F.
Therefore, if the robot center of gravity and geometric center do not coincide, the front two rotors need to be increased by 0.5 × ∆F for each rotor.The lift generated by the front two rotors is F 1 , and the lift generated by the rear two rotors F 2 .
The total power W tot is greater than the total power in the case of uniform mass distribution.
As shown in Figure 10a, we use the flight controller to rotate the propeller at different speeds.The relationship between the lift and power of different types of propellers is calculated by the push-pull sensor and the voltage-current sensor, and the curves are fitted.
The curves of different types of propellers show a quadratic curve relationship as shown in Figure 10b, which satisfies the condition of Equation (28).As shown in Figure 11, the mass of the robot (LAWCDR) can be divided into three parts, the body mass, the first joint mass, and the second joint mass (suction cup and support bars).s L denotes the y-axis moment of the first joint in the body coordinate system.v L denotes the y-axis moment of the second joint in the body coordinate sys- tem.During the flight phase, we optimize the servo angle by calculating the horizontal and vertical moments that the center of gravity offset causes.As shown in Figure 11, the mass of the robot (LAWCDR) can be divided into three parts, the body mass, the first joint mass, and the second joint mass (suction cup and support bars).L s denotes the y-axis moment of the first joint in the body coordinate system.L v denotes the y-axis moment of the second joint in the body coordinate system.During the flight phase, we optimize the servo angle by calculating the horizontal and vertical moments that the center of gravity offset causes.As shown in Figure 11, the mass of the robot (LAWCDR) can be divided in parts, the body mass, the first joint mass, and the second joint mass (suction c support bars).s L denotes the y-axis moment of the first joint in the body co system.v L denotes the y-axis moment of the second joint in the body coordin tem.During the flight phase, we optimize the servo angle by calculating the ho and vertical moments that the center of gravity offset causes.
where S 1 is the horizontal offset moment.S 2 is the vertical offset moment.
Since the robot (LAWCDR) is not completely centrosymmetric, the mass and position of optional components such as the battery, on-board computer, camera, and image transmission are not completely fixed.Therefore, there is an unknown initial moment.It is necessary to observe this moment by an extended state observer (ESO) [27].The robot adjusts the servo angle according to the performance function Equation (31).The design of the extended state observer is shown in the following equation. where , δ is an adjustable threshold.a 1 and a 2 are positive adjustable parameters.θ is the observation of the pitch angle.β 1 , β 2 , and β 3 are the gain of the extended state observer, respectively.b is the scaling factor.f al(e, a i , δ) is the nonlinear function of the extended state observer.The final pitch angular velocity control quantity u is shown in the following equation.
where u θ is the original output of the controller.
is the disturbance moment.δ(t) is the moment due to the deformation of the servo, according to Equation (23).j( .ϕ, .ψ) is the pitch channel coupling term.The convergence of the observer is proved as follows.We make the following assumptions [27][28][29][30][31] in building the observer error model based on the robot dynamics model.
where v is the lump disturbance.Defining ψ + d(t) + δ(t), the system can be rewritten as: The observer can be designed as where β 1 , β 2 , β 3 are positive adjustable parameters.
According to Equations ( 35) and ( 36), the error model of the observation can be expressed as . e 2 = e 3 − g 2 (e 1 ) .
where ξ is the observer residual.The value is very small.Define the Lyapunov function V(e 1 (t), e 2 (t), e 3 (t)), and take the time derivation of V(t) shown as Equation (38): .

Simulation of the Initial Moment Observed by the ESO
To verify the effectiveness of the proposed extended state observer, we conducted simulations in Matlab/Simulink.The simulation uses fixed-step (t = 0.001) and ODE3 solver.In the simulation, we adopted the classic cascade PID control in the position control loop and attitude control loop.It is assumed that the robot has an initial moment of 1.5 Nm due to the non-uniform mass of the ground module.
As shown in Figure 12, the robot adjusted the servo angle from the preparation state to the flight state at 6 s.We calculated the average value of the output, marked by the red line.The pitch channel output was reduced from 1.50 Nm to 0.56 Nm.

Simulation of the Initial Moment Observed by the ESO
To verify the effectiveness of the proposed extended state observer, we conducted simulations in Matlab/Simulink.The simulation uses fixed-step (t = 0.001) and ODE3 solver.In the simulation, we adopted the classic cascade PID control in the position control loop and attitude control loop.It is assumed that the robot has an initial moment of 1.5 Nm due to the non-uniform mass of the ground module.
As shown in Figure 12, the robot adjusted the servo angle from the preparation state to the flight state at 6 s.We calculated the average value of the output, marked by the red line.The pitch channel output was reduced from 1.50 Nm to 0.56 Nm.The output of the extended state observer is shown in Figure 13.The initial moment observed by the ESO is −1.5 Nm.After servos angle change at 6 s, it changes to 0.55 Nm, which is roughly the same amount as the change in U3.The ESO observer is able to converge within 1 s.

Simulation of the Initial Moment Observed by the ESO
To verify the effectiveness of the proposed extended state observer, we conducted simulations in Matlab/Simulink.The simulation uses fixed-step (t = 0.001) and ODE3 solver.In the simulation, we adopted the classic cascade PID control in the position control loop and attitude control loop.It is assumed that the robot has an initial moment of 1.5 Nm due to the non-uniform mass of the ground module.
As shown in Figure 12, the robot adjusted the servo angle from the preparation state to the flight state at 6 s.We calculated the average value of the output, marked by the red line.The pitch channel output was reduced from 1.50 Nm to 0.56 Nm.

Bionic Power Optimization Experiments
As shown in Figure 14, the robot (LAWCDR) includes a Robomaster M2006 front wheel motor, four JuFeng 2814 kv700 motors, four 40A ESCs, four APC1055 propellers, a 60 mm double-layer nitrile rubber suction cup, and two Feetech SM45BL RS485 servos.The flight controller is CUAV fmuv5 nano flight controller.The wheel diameter is 14 cm.The onboard computer is Nvidia Jetson Xavier NX, through the mavlink protocol and ROS for planning trajectory.

Bionic Power Optimization Experiments
As shown in Figure 14, the robot (LAWCDR) includes a Robomaster M2006 front wheel motor, four JuFeng 2814 kv700 motors, four 40A ESCs, four APC1055 propellers, a 60 mm double-layer nitrile rubber suction cup, and two Feetech SM45BL RS485 servos.The flight controller is CUAV fmuv5 nano flight controller.The wheel diameter is 14 cm.The onboard computer is Nvidia Jetson Xavier NX, through the mavlink protocol and ROS for planning trajectory.As shown in Figure 15, the robot position feedback indoors is output through the Optitrack motion capture system and sent to the robot through the Xbee communication module.The outdoor position feedback is output via differential GPS and the RTCM differential data from the base station is sent to the robot via Xbee as well.The robot performs position decomposition to obtain outdoor position feedback.As shown in Figure 15, the robot position feedback indoors is output through the Optitrack motion capture system and sent to the robot through the Xbee communication module.The outdoor position feedback is output via differential GPS and the RTCM differential data from the base station is sent to the robot via Xbee as well.The robot performs position decomposition to obtain outdoor position feedback.
As shown in Figure 14, the robot (LAWCDR) includes a Robomaster M2006 front wheel motor, four JuFeng 2814 kv700 motors, four 40A ESCs, four APC1055 propellers, a 60 mm double-layer nitrile rubber suction cup, and two Feetech SM45BL RS485 servos.The flight controller is CUAV fmuv5 nano flight controller.The wheel diameter is 14 cm.The onboard computer is Nvidia Jetson Xavier NX, through the mavlink protocol and ROS for planning trajectory.As shown in Figure 15, the robot position feedback indoors is output through the Optitrack motion capture system and sent to the robot through the Xbee communication module.The outdoor position feedback is output via differential GPS and the RTCM differential data from the base station is sent to the robot via Xbee as well.The robot performs position decomposition to obtain outdoor position feedback.To verify the power optimization effect of the robot, we tested the power and control amount of the robot in different operating states when the robot was hovering in the air.The angular velocity control output of the pitch channel in the preparation state and flight state is shown in Figure 16.To better illustrate the variation, mean filtering of n = 200 was performed (shown as the red line).The control of the pitch axis is normalized to [-1,1].In order to balance the non-coincidence of the center of gravity and geometric center, the robot needs to reduce the control margin of the pitch axis.From the preparation state to the flight state, the amount of control is reduced from 0.38 to 0.15, and the control margin is increased from 0.62 to 0.85.This causes an increase of 37.1% in the control margin.Compared to the preparation state, the flight state increases the control margin and improves the control performance of the robot.
tion state to the flight state, the amount of control is reduced from 0.38 to 0 control margin is increased from 0.62 to 0.85.This causes an increase of 3 control margin.Compared to the preparation state, the flight state increases margin and improves the control performance of the robot.The total robot power is shown in Figure 17.We also performed the me with n = 200 (as shown in the red line).The total robot hovering power is a when the robot is in the preparation state and about 780 W when the robot is state.The power difference is about 6%.If the robot works in flight state, t time of the robot can be extended.Taking a 6 s 5000 mah battery as an examp working time can be extended from 481 s to 512 s.The total robot power is shown in Figure 17.We also performed the mean filtering with n = 200 (as shown in the red line).The total robot hovering power is about 830 W when the robot is in the preparation state and about 780 W when the robot is in the flight state.The power difference is about 6%.If the robot works in the flight state, the working time of the robot can be extended.Taking a 6 s 5000 mah battery as an example, the robot working time can be extended from 481 s to 512 s.
margin and improves the control performance of the robot.The total robot power is shown in Figure 17.We also performed the m with n = 200 (as shown in the red line).The total robot hovering power is when the robot is in the preparation state and about 780 W when the robot i state.The power difference is about 6%.If the robot works in the flight state time of the robot can be extended.Taking a 6 s 5000 mah battery as an exam working time can be extended from 481 s to 512 s.

Actual Motion Process Experiments (Including Bionic Adsorption)
The main working state transitions of the robot (LAWCDR) are air-wall transformation and land-air transformation.The air-wall transformation process includes the flight adsorption process and the inverse process.The robot takes off from the ground, hovers at a height of 1 m, flies to the front of the wall, and performs flight adsorption.Firstly, the robot adjusts the servos to the preparation state to reduce the reaction time of performing the adsorption action.Secondly, the robot adjusts its heading to face the wall.Thirdly, the robot flies toward the wall at a slow, desired speed.When the wheels contact the wall, the shock absorbers are compressed to further adjust the robot's heading.Fourthly, the flight control detects the contact.The servos horizontally stretch out the suction cup to contact the wall and the robot finally completes the flight adsorption process.
The robot performs a fast contact detection through the IMU.Once the acceleration of the x-axis that exceeds a threshold momentarily is detected, the robot controls the servos to start moving from the preparation position and stretches out the suction cup to contact the wall.As shown in Figure 18, after 138 s, the robot contacts the wall, and the actual position of the robot reaches 2 m without further change (as shown in the red line).The x-axis acceleration (as shown in the yellow line) steps from 0 to 20 m/s/s, exceeding the threshold value (8 m/s/s).The main working state transitions of the robot (LAWCDR) are air-wall transformation and land-air transformation.The air-wall transformation process includes the flight adsorption process and the inverse process.The robot takes off from the ground, hovers at a height of 1 m, flies to the front of the wall, and performs flight adsorption.Firstly, the robot adjusts the servos to the preparation state to reduce the reaction time of performing the adsorption action.Secondly, the robot adjusts its heading to face the wall.Thirdly, the robot flies toward the wall at a slow, desired speed.When the wheels contact the wall, the shock absorbers are compressed to further adjust the robot's heading.Fourthly, the flight control detects the contact.The servos horizontally stretch out the suction cup to contact the wall and the robot finally completes the flight adsorption process.
The robot performs a fast contact detection through the IMU.Once the acceleration of the x-axis that exceeds a threshold momentarily is detected, the robot controls the servos to start moving from the preparation position and stretches out the suction cup to contact the wall.As shown in Figure 18   As shown in Figure 20, the robot adjusts the desired attitude angle from 0 degrees to −25 degrees after contact is detected at 138 s.The desired and actual roll angles remain at about 0 degrees.The actual pitch angle changes from 0 to −25 degrees, and a horizontal component force is applied to the wall.This assures that the wheels are attached to the wall before the suction cup contacts the wall.During the roost state, the desired pitch angle is 0 and the rotors stop working.The actual pitch angle of the robot is kept at −10 degrees due to gravity.As shown in Figure 20, the robot adjusts the desired attitude angle from 0 degrees to −25 degrees after contact is detected at 138 s.The desired and actual roll angles remain at about 0 degrees.The actual pitch angle changes from 0 to −25 degrees, and a horizontal component force is applied to the wall.This assures that the wheels are attached to the wall before the suction cup contacts the wall.During the roost state, the desired pitch angle is 0 and the rotors stop working.The actual pitch angle of the robot is kept at −10 degrees due to gravity.The angle curves of servo 1 and servo 2 during the adsorption process are shown in Figure 21.When the contact is detected, the angle of servo 1 changes from 0 degrees to 24.5 degrees.The angle of servo 2 changed from 93 degrees to 41 degrees.Therefore, the servos keep horizontally stretching out the suction cup to contact the wall.The angle curves of servo 1 and servo 2 during the adsorption process are shown in Figure 21.When the contact is detected, the angle of servo 1 changes from 0 degrees to 24.5 degrees.The angle of servo 2 changed from 93 degrees to 41 degrees.Therefore, the servos keep horizontally stretching out the suction cup to contact the wall.The angle curves of servo 1 and servo 2 during the adsorption process are shown in Figure 21.When the contact is detected, the angle of servo 1 changes from 0 degrees to 24.5 degrees.The angle of servo 2 changed from 93 degrees to 41 degrees.Therefore, the servos keep horizontally stretching out the suction cup to contact the wall.

Conclusions
This paper develops a land-air-wall cross-domain robot that combines the advantages of the land-air amphibious robot and the flying adsorption amphibious robot and imitates the process of geckos using their spines to complete adsorption again after collision with the wall.The robot imitates the bionic adsorption process and uses servos and shock absorbers to simulate the spine.In addition, inspired by the phenomenon of

O
denotes the body coordinate syst denotes the lift generated by the ith rotor.The angular bisector of the second an motors is taken as the y-axis, and the angular bisector of the third and fourth mo taken as the x-axis.h O and c O denote the coordinate systems represented by t and second joints, and c O coincides with the suction cup end coordinate syste

O
of the body coordinate system.
, (a) the wheels (rear cla the wall.(b) The wheels (rear claws) are compressed to produce a The robot uses the lift force generated by the propeller to attach the (d) The suction cup (head) and support bars (front claws) are driven
, (a) the wheels (rear claws) start to contact the wall.(b) The wheels (rear claws) are compressed to produce a tilting moment.(c) The robot uses the lift force generated by the propeller to attach the wheels to the wall.(d) The suction cup (head) and support bars (front claws) are driven towards the wall by adjusting the posture with the servos (spine).(e) The suction cup is attached to the wall, establishes a negative pressure environment, and makes the robot adsorb to the wall.The rotor stops working.Panel (f) indicates the bionic function of each component.

Figure 4 .
Figure 4. Adsorption process of the robot (LAWCDR) compared with the gecko (rotor lift keeps the wheels fixed on the wall).(a) The wheels start to contact the wall.(b) The wheels are compressed.(c) The robot uses the lift force generated by the propeller to attach the wheels to the wall.(d) The suction cup and support bars are driven towards the wall.(e) The suction cup makes the robot adsorb to the wall.(f) The correspondence between geckos and robots.The structural dimensions of the robot (LAWCDR) are shown in Figure 5.The model of the wheel-wall contact can be represented as 2 2 2 2 1 2 3 4 2 2 2 2 3 3 1 2 3 4

F
are the equivalent contact forces of the left and right wheel on the wall.wl d and wr d are the z-axis distances of the left and right wheel contact points

Figure 4 .
Figure 4. Adsorption process of the robot (LAWCDR) compared with the gecko (rotor lift keeps the wheels fixed on the wall).(a) The wheels start to contact the wall.(b) The wheels are compressed.(c) The robot uses the lift force generated by the propeller to attach the wheels to the wall.(d) The suction cup and support bars are driven towards the wall.(e) The suction cup makes the robot adsorb to the wall.(f) The correspondence between geckos and robots.The structural dimensions of the robot (LAWCDR) are shown in Figure 5.The model of the wheel-wall contact can be represented as

Figure 5 .
Figure 5. Diagram of the spring mechanism installed on the front wheels.The red circles represent rotating shaft, orange device is the shock absorber.

skand 2 h
is the coefficient of elas- ticity of the spring.w r is the rotation radius of the spring assembly.q r is the rotation radius of the wheel assembly.1 h are the x-axis distance and z-axis distance of the spring mount point in the body coordinate system, respectively.L

Figure 5 .
Figure 5. Diagram of the spring mechanism installed on the front wheels.The red circles represent rotating shaft, orange device is the shock absorber.

1 f
denotes the friction between the suction cup and the wall, and 2 f denotes the friction between the support bars and the wall.

1 N and 2 N
are the support forces of the wall for the suction cup and the support bars, respectively.

Figure 6 .
Figure 6.Force analysis of the robot.(a) Force analysis of the approaching state (propellers are working).(b) Force analysis of the adsorbing state (propellers are not rotating).Figure 6. Force analysis of the robot.(a) Force analysis of the approaching state (propellers are working).(b) Force analysis of the adsorbing state (propellers are not rotating).

Figure 6 .
Figure 6.Force analysis of the robot.(a) Force analysis of the approaching state (propellers are working).(b) Force analysis of the adsorbing state (propellers are not rotating).Figure 6. Force analysis of the robot.(a) Force analysis of the approaching state (propellers are working).(b) Force analysis of the adsorbing state (propellers are not rotating).

Figure 7 .Figure 7 .
Figure 7. Structure of the suction cup and the support bars.

Figure 8 . 2 N 2 N becomes larger, 2 N 1 N
Figure 8. Influence of angle  to the support force

Figure 8 .
Figure 8. Influence of angle α to the support force N 2 on absorbing module.

Figure 9 .
Figure 9. Influence of the support bar angles s  to the support force 2 N .
containing all the generalized coordinate variables.

Figure 9 .
Figure 9. Influence of the support bar angles θ s to the support force N 2 .

.
p i and ω i are the translational and angular velocity of ith joint in the inertial frame.m b and m i are the mass of the body and the ith joint.I b and I i are the rotational inertia of the body and the ith joint.R = R w b .R i denotes the rotation matrix from the body frame to the ith joint.

Figure 10 .
Figure 10.Propeller lift and actual power test: (a) lift/power test bench (b) fitting curve between lift and power consumption.

Figure 11 .
Figure 11.The robot balancing the center of gravity by adjusting the servo.

Figure 10 .
Figure 10.Propeller lift and actual power test: (a) lift/power test bench (b) fitting curve between lift and power consumption.

Figure 10 .
Figure 10.Propeller lift and actual power test: (a) lift/power test bench (b) fitting curve bet and power consumption.

Figure 11 .
Figure 11.The robot balancing the center of gravity by adjusting the servo.

Figure 11 .
Figure 11.The robot balancing the center of gravity by adjusting the servo.

Figure 12 .
Figure 12.Simulation of control quantity of robot pitch channel.The output of the extended state observer is shown in Figure 13.The initial moment observed by the ESO is −1.5 Nm.After servos angle change at 6 s, it changes to 0.55 Nm, which is roughly the same amount as the change in U3.The ESO observer is able to converge within 1 s.

Figure 13 .
Figure 13.Simulation on robot extended state observer.

Figure 12 .
Figure 12.Simulation of control quantity of robot pitch channel.

Figure 12 .
Figure 12.Simulation of control quantity of robot pitch channel.The output of the extended state observer is shown in Figure 13.The initial moment observed by the ESO is −1.5 Nm.After servos angle change at 6 s, it changes to 0.55 Nm, which is roughly the same amount as the change in U3.The ESO observer is able to converge within 1 s.

Figure 13 .
Figure 13.Simulation on robot extended state observer.Figure 13.Simulation on robot extended state observer.

Figure 13 .
Figure 13.Simulation on robot extended state observer.Figure 13.Simulation on robot extended state observer.

Figure 15 .Figure 14 .
Figure 15.Setup of the robot in the experiments.To verify the power optimization effect of the robot, we tested the power and control amount of the robot in different operating states when the robot was hovering in the air.The angular velocity control output of the pitch channel in the preparation state and flight state is shown in Figure16.To better illustrate the variation, mean filtering of n =

Figure 15 .Figure 15 .
Figure 15.Setup of the robot in the experiments.To verify the power optimization effect of the robot, we tested the power and control amount of the robot in different operating states when the robot was hovering in the air.The angular velocity control output of the pitch channel in the preparation state and flight state is shown in Figure16.To better illustrate the variation, mean filtering of n =

Figure 16 .
Figure 16.The output control quantity at preparation state and flight state.

Figure 17 .
Figure 17.Total power at preparation state and flight state.

Figure 16 .
Figure 16.The output control quantity at preparation state and flight state.

Figure 16 .
Figure 16.The output control quantity at preparation state and flight state.

Figure 17 .
Figure 17.Total power at preparation state and flight state.

Figure 17 .
Figure 17.Total power at preparation state and flight state.
, after 138 s, the robot contacts the wall, and the actual position of the robot reaches 2 m without further change (as shown in the red line).The x-axis acceleration (as shown in the yellow line) steps from 0 to 20 m/s/s, exceeding the threshold value (8 m/s/s).

Figure 18 .
Figure 18.Process of the IMU-detecting robot contacting the wall.The robot imitates a gecko for adsorption as shown in Figure19, where (a) the robot is hovering in the air and the suction cup moves to the preparation state for adsorption.(b) The robot starts to fly toward the wall.(c) The front wheels start to contact the wall (the contact point between the right wheel and the wall is the blue point).The shock absorbers start to work, adjusting the robot's heading.(d) The suction cup driven by the servos moves towards the painted wall.As shown in Figure19e,f, the suction cup contacts the wall (the contact point between the suction cup and the wall is the red point) and a negative pressure environment is established.

Figure 18 .
Figure 18.Process of the IMU-detecting robot contacting the wall.The robot imitates a gecko for adsorption as shown in Figure19, where (a) the robot is hovering in the air and the suction cup moves to the preparation state for adsorption.(b) The robot starts to fly toward the wall.(c) The front wheels start to contact the wall (the contact point between the right wheel and the wall is the blue point).The shock absorbers start to work, adjusting the robot's heading.(d) The suction cup driven by the servos moves towards the painted wall.As shown in Figure19e,f, the suction cup contacts the wall (the contact point between the suction cup and the wall is the red point) and a negative pressure environment is established.As shown in Figure20, the robot adjusts the desired attitude angle from 0 degrees to −25 degrees after contact is detected at 138 s.The desired and actual roll angles remain at about 0 degrees.The actual pitch angle changes from 0 to −25 degrees, and a horizontal component force is applied to the wall.This assures that the wheels are attached to the wall before the suction cup contacts the wall.During the roost state, the desired pitch angle is 0 and the rotors stop working.The actual pitch angle of the robot is kept at −10 degrees due to gravity.

Figure 19 .
Figure 19.Robot flight adsorption step diagram (the red point is the contact point between the suction cup and the wall, and the blue point is the contact point between the right wheel and the wall).

Figure 19 .
Figure 19.Robot flight adsorption step diagram (the red point is the contact point between the suction cup and the wall, and the blue point is the contact point between the right wheel and the wall).

Figure 20 .
Figure 20.Attitude curve of the robot.(The different working states are divided by the yellow dotted line).

Figure 20 .
Figure 20.Attitude curve of the robot.(The different working states are divided by the yellow dotted line).

Figure 20 .
Figure 20.Attitude curve of the robot.(The different working states are divided by the yellow dotted line).

Figure 21 .
Figure 21.Servo angle curve of the robot.(The different working states are divided by the yellow dotted line).The trajectory of the land and air transformation process of the robot is shown in Figure22.The experiment to verify the land-air motion movement capability of the robot was carried out.The robot follows a planned trajectory of ground motion, air flight, takeoff, and landing.The robot moves from point A to point B, takes off at point B and

Figure 21 .
Figure 21.Servo angle curve of the robot.(The different working states are divided by the yellow dotted line).The trajectory of the land and air transformation process of the robot is shown in Figure22.The experiment to verify the land-air motion movement capability of the robot was carried out.The robot follows a planned trajectory of ground motion, air flight, takeoff, and landing.The robot moves from point A to point B, takes off at point B and moves to air point C, then flies through points D and E to point F, lands at point G, and then moves to points H and I on the ground.The entire movement process is completely autonomous.
Appl.Sci.2022, 12, x FOR PEER REVIEW 24 of 26 moves to air point C, then flies through points D and E to point F, lands at point G, and then moves to points H and I on the ground.The entire movement process is completely autonomous.

Figure 22 .
Figure 22.Actual land-air trajectory of the robot.

Figure 22 .
Figure 22.Actual land-air trajectory of the robot.