Motion Control of a Hybrid Quadruped-Quadrotor Robot

: Multimodal motion capability is an emerging topic in the robotics field, and this paper presents a hybrid robot system maneuvering in both terrestrial and aerial environments. Firstly, a micro quadruped–quadrotor robot with onboard sensing and computing is developed. This robot incorporates both the high mobility of unmanned aerial vehicles and the long endurance of mobile robots on the ground. A coordinated motion control scheme is then exploited for adaptive terrestrial– aerial motion transition. In this scheme, a bio-inspired terrestrial locomotion controller is proposed to generate various quadruped locomotions, and a model-based aerial locomotion controller is proposed to generate various quadrotor configurations. Then, an unified motion controller for the two subsystems which dynamically adjusts crawling and flying motion in a complicated environment is presented. Consequently, several practical trials are conducted to demonstrate the adaptability and the robustness of the proposed system.


Introduction
Robots with multimodal locomotion capabilities, which benefit from combinational advantages to overcome the individual limitations, have drawn increasingly great attention in robotics.In order to maneuver in versatile domains such as terrestrial, aerial, and aquatic, several state-of-the-art hybrid robots have been explored [1][2][3][4][5][6][7][8][9].Among these multimodal forms, legged systems with dexterous locomotion possess advanced potential in irregular terrain like walls [10], slopes [11], and uneven terrain [12], but the legged defect limits large-scale mobilities and tremendous obstacles.In nature, insects usually have evolved the ability of both crawling and flying over millions of years of natural selection.Much progress has been made in figuring out and replicating animal locomotion.As a result, robotics can benefit from these biological counterparts to develop corresponding hybrid behaviors.In order to enhance the adaptability of insect-like robots, several novel studies on leg-multirotor robots have attempted to adopt the ideas and mechanisms from nature: in [13], the authors investigate a land-air robot that performs structural separation with the four-legged robot and the quadrotor for multipurpose investigation, so the size of the quadrotor is fixed and only the ground motion is discussed; a novel air-ground amphibious quadruped robot called SPIDAR is developed to explore an integrated PIDbased control strategy and optimization-based control allocation in [14], but the computing cost and solving complexity are increasing for such a large 16-rotor and 8-joint robot in real-time scenarios; the authors of [15] presented a hexarotor with three-manipulator with the advantages of taking off and landing on unstable surface, but the leg-like manipulators are only used to stand rather than walking; a walking-and-flying robotic system based on the ANYmal quadruped robot is proposed for collaborative exploration in [16], but the system is divided into two subsystems to explore individually; the authors of [17] designed a multimodal locomotion robot called LEONARDO with two legs, but the bipedal motions require strong stability and delicate balance on uneven terrains via complicatedly controlling the rotary propellers of aerial locomotion.
Although many multimodal robot designs with both terrestrial locomotion and aerial locomotion have been widely investigated in recent years, only a few of these robots can deal with uncertainties in practical environments.The robot size, motion agility, and control complexity have a great effect on limited spaces and costs.In this work, a hybrid Quadruped-Quadrotor Robot, namely QQR, is customized, and either increases the mechanical complexity or decreases the controllability of the system.Owing to the long power endurance and great stability of quadrupeds as well as the excellent maneuverability and transformable flight of quadrotors [18], the proposed QQR may be applied for long-term missions and extreme terrains such as search and rescue and active detection, as well as by appropriately switching between available locomotion modes, while a single quadrotor with relatively short endurance can hardly complete the entire mission.Furthermore, the deformation characteristics of the QQR with several controllable degrees of freedom can modify any configuration during flight to undertake specific tasks [19] such as traversing a narrow gap to enter a collapsed building in Figure 1: firstly, the QQR flies fast with standard X configuration; then, the QQR switches to the H configuration for smaller size to pass through the narrow gap; finally, after entering the building, the QQR uses the T configuration for closer inspection, and the lack of flight condition will enforce the crawling modes of the QQR to localize the targets.In order to execute such collaborative missions, a hierarchical coordination control scheme is further proposed to improve the performance in practical environments.Several experiments and comparisons in actual environments are conducted to verify the effectiveness and reliability of the proposed hierarchical coordination control scheme.The results demonstrate that the proposed controllers can generate and transfer the locomotion of the robot with the changing environments.Contributions of this paper are as follows: (1) A micro, lightweight, and integrated robot prototype that enables both crawling motion and morphing flight.(2) A hierarchical coordination control scheme consisting of a low-level locomotion controller and a high-level motion controller that provides hybrid air-ground robots with multimodal locomotion capabilities.(3) Implementing foundational experiments and comparisons on such a hybrid robotic system to address the challenge of practical applications.The rest of this paper is organized as follows.The model of the quadruped-quadrotor robot is detailed in Section 2. Section 3 proposes the bio-inspired control scheme.The simulation and experimental results are carried out in Section 4. Finally, Section 5 gives the conclusion.

System Overview
In this section, the mechanical design of the QQR combining the peculiarity of both a common quadrotor and a bio-inspired quadruped is first presented.The drawbacks as short endurances in quadrotors and low rates in quadrupeds can be compensated by the adaptivity in our designed QQR.Then, the mathematical models are established separately for further control.

Robot Design
Figure 2 illustrates the mechanical configuration of the QQR prototype consisting of four legs and one body.Each leg with two degrees of freedom (DoF) is equipped with one brushless motor (connected to the propeller) and two servo motors (connected together to act as two actuated joints).The lightweight legs are significantly vital for the flight performance and morphing motion.A microprocessor (CUAV Nora+ autopilot) is located in the body to send control commands to brushless motors (TMOTOR F40) and servo motors (KINGMAX KM0940), as well as receiving decision commands from the upper intelligent computing module (Jetson Xavier NX).A stereo camera (RealSense D435), a GPS (NEO 3) and an integrated optical flow (HereFlow) are employed to sense environmental information.Besides that, a customized printed circuit board (PCB) is designed for the reliability of the electrical connections and a 4S battery is embedded under the body.Hence, the entire robotic system can be divided into legged module, flight module and control module.Moreover, the specifications of the QQR prototype are given in Table 1.

Modeling
In this section, the mathematical model of the QQR is divided into two parts: the multi-joint torque model for quadrupeds and the multi-rotor thrust model for quadrotors, as illustrated in Figure 3.

Terrestrial Dynamics of Multi-Joint Torque
In order to build the terrestrial dynamics of quadrupeds, the joint coordinates of the robot should be firstly defined, as depicted in Figure 4.The kinematic model is represented by Denavit-Hartenberg (DH) parameters for resolving inverse kinematic of each leg.According to these fixed frames, the transformation parameters and DH parameters are demonstrated in Table 2 and Table 3, respectively.
Table 2. Transformation parameters from the {O i0 } to the {O B }.The thigh joint position is defined as (d ix , d iy ).φ i denotes the direction angle in the body frame {O B }.Where l ix = l iy = 5 and Table 3. Denavit-Hartenberg parameters, where j = 1, 2 is the joint number from thigh joint to shin joint of each leg.And α ij is the link twist indicating the angle from Z i(k−1) to Z ik around X ik ; a ij is the joint distance meaning the distance from the intersection of Z i(k−1) and X ik to the origin of Z i(k−1) ; d ij is the link length representing the distance from the intersection of Z i(k−1) and X ik to the origin of X ik ; θ ij is the joint angle showing the angle from X i(k−1) to X ik around Z i(k−1) .

Joint j α
The relative translation and rotation between the (j − 1) th and the j th joint coordinates are computed by the transformation matrix (1): In addition, the transition matrix between the body coordinate {O B } and the hip joint coordinate {O i0 } is represented by (2): Consequently, the foot tip coordinate {O i2 } can be transformed into the body coordinate {O B } by multiplying the previous matrixes sequentially: Thus, the position of the foot tip with respect to the body coordinate {O B } can be derived from the below: where [p xi p yi p zi ] T is the position coordinate of the i th foot tip and θ ij is the joint angle of the joint j on the leg i.
The leg of the QQR is designed to connect the body trunk with the ground as a complex joint-link system, so closed kinematics chains exist in the QQR system.Since forces and moments propagate via the kinematics chains among different legs [20], the kinematics and dynamics of the robot are coupled.The dynamic model of such a coupled quadruped robot with eight actuators is derived from the Lagrangian-Euler method: where is the joint torque vector of the i th leg consisting of thigh joint torque τ i1 and shin joint torque τ i2 .
] T ∈ R 2 are the joint angle, joint angle velocity, and joint angle acceleration vector of the i th leg, respectively.
is the Coriolis force matrix of the i th leg.G i (θ) ∈ R 2 is the link gravitational force vector of the i th leg.
represents ground reaction forces of the i th foot tip with the coordinate {O i2 }, which are solved by the minimization of norm of feet forces approach [20].J i ∈ R 3×2 is the Jacobian matrix of the i th leg, computed by (6).To clarify the controller design, the position and velocity of the QQR in this work are transformed to the global coordinate {O G }.

Aerial Dynamics of Multi-Rotor Thrust
The aerial coordinates of each rotor i are shown in Figures 3 and 4. The six-DoF rigid body kinematic and dynamic equations are adopted to establish a quadrotor model with including the position p G ∈ R 3 , the linear velocity v G ∈ R 3 , the unit quaternion q G ∈ SO 3 , the body angular rate ω B ∈ R 3 , and the rotational speed vector of the four rotors where Ω Bi represents the rotational speed of the i th rotor [21].
The translational dynamics can be expressed by where f t and f d are the collective thrust force and the external drag force, respectively.M is the total mass of the QQR, g = [0, 0, −g] T is the Earth's gravitational vector, and z B is the Z axis of the body coordinate O B in the global coordinate O G .
The rotational dynamics can be written as where [•] * is the skew-symmetric matrix.τ and τ d are the total torque of the QQR and the body torque of model uncertainties, respectively.J is the diagonal inertial tensor matrix.
The rotational speed vector of the four rotors Ω B is set as a first-order system: where Ω d is the desired speed and k cons is the time constant.Hence, the collective thrust f t and torque τ produced by the rotors can be calculated as where f i (i = 1, 2, 3, 4) is the individual motor thrust on each rotor and f = [ H is the time-variant control allocation matrix, expressed by [22] where c t and c k declare the thrust coefficient and the torque coefficient, respectively.r = (r x , r y , r z ) is the position of the center of gravity (COG).It should be noted that all of the motors have limited thrust thresholds [ f min , f max ]:

Hierarchical Coordination Control Scheme
The hierarchical coordination control scheme by controlling rotors and articulated joints is depicted in Figure 5, which consists of a separate locomotion control for locomotion generation at low levels and an unified motion control for motion transition at high levels.

Low-Level Locomotion Control
In low-level locomotion control, to ensure the controllability and stability of flying and crawling simultaneously on such a hybrid QQR, a terrestrial locomotion controller and an aerial locomotion controller are proposed to produce active terrestrial locomotion (such as walk, trot, translation, and rotation) and aerial locomotion (such as X, O, H, and T configurations), respectively.

Terrestrial Locomotion Control
In nature, biological central pattern generations (CPGs) are mainly involved in the generation of coordinated and rhythmic movements for animals [23].Considering the similarity between legged animals and quadruped robots, a bio-inspired terrestrial locomotion controller based on central pattern generation (CPG) is adopted in our work.Moreover, the CPG-based model is able to couple dynamics of the QQR to produce various locomotion behaviors, as well as being modulated further via a high-level motion controller [24] in Figure 6.Similar to their biological counterparts, the QQR performs dexterous terrestrial locomotion capabilities in cluttered terrains, such as walk locomotion for slow movements, trot locomotion for fast running, and translation locomotion and rotation locomotion for agile turning [25].The aforementioned terrestrial dynamics model helps the QQR to learn various locomotion patterns.The terrestrial locomotions and the phase characteristics are illustrated in Figure 7.In every cycle, four legs in the walk locomotion are separated into four sets of {leg 1}, {leg 2}, {leg 3}, and {leg 4} to move successively, while in the trot locomotion, they are separated into two sets of {leg 1, leg 3} and {leg 2, leg 4}, moving alternately.Moreover, the translation locomotion and the rotation locomotion are supposed to accomplish gentle turns and sharp turns, respectively.In this work, the locomotion patterns can be altered in different situations adaptively.Owing to stable limit cycles and explainable parameters, Hopf oscillators are suitable for building CPG networks for robotic locomotion patterns [26].A traditional Hopf oscillator can be formulated by ( 16): where x and y are two state variables.µ and ω represent the amplitude and the frequency of the Hopf oscillator.α and β are constants determining the convergence rate of the limit cycle.In this work, x is defined as the output signal.
As the proposed QQR in this work is equipped with four 2-DoF legs, eight Hopf oscillators are utilized in the QQR's CPG network and need to be connected together to guarantee the performance of synchronization and coordination.In order to generate desired terrestrial locomotion, we propose a terrestrial CPG network based on the aforementioned multi-joint dynamic model, and the topology of the CPG network is illustrated as shown in Figure 8. Due to complicated couplings and high dimensions on the QQR, the proposed CPG-based locomotion control is decomposed into two layers: (1) The thigh layer includes four oscillators with bidirectional couplings.(2) The shin layer consists of four oscillators with unidirectional connection under the control of the corresponding oscillator in the thigh layer.The thigh layer provides priori knowledge to determine locomotion patterns of the QQR, while the shin layer fine-tunes locomotion performance for quick response to external contact environments.Motivated by the work presented by [27], the proposed CPG model with diffusive couplings has the advantage of dealing with kinematic constraints, which can be simplified and described by where x ij , y ij with the leg number i = 1, 2, 3, 4 and the joint number j = 1, 2 denote two state variables, and x ij is chosen as the output signal.The constant coupling strength k = 0.1 and the convergence coefficients α = β are set for all oscillators, which are determined through trial-and-error simulations on the stability of limit cycle in this work.The oscillator frequencies are unified as ω ij = ω for simplifying the further coordinated motion control.
Additionally, θ ij mn with m = 1, 2, 3, 4 and n = 1, 2 represents the phase difference between the joint ij and the joint mn.The associated 2D rotation matrix R(θ ij mn ) can be written as This proposed two-layer terrestrial CPG model not only introduce the coupling relations among various Hopf oscillators to achieve locomotion patterns of the QQR, but also can fine-tune locomotion performance via explicit CPG parameters to improve environmental adaptability.The proposed CPG-based architecture effectively reduces the control difficulty of the eight-joint QQR.Although there are only fewer parameters like amplitudes, frequencies, and phase differences in the proposed terrestrial locomotion controller, the parameter modulation process affects the interacting performance with environments that need to be fine-tuned by a high-level motion controller.While morphing in midair with the overlap between the propellers and the central body of the robot, the QQR will experience a sudden loss of thrust to deviate from the original position.To overcome the effect of fast changes in aerial configurations, a robust geometry-aware compensation controller aiming at modifying the motor speeds of the propeller is deployed as described in Figure 10.The individual thrust f i provided by a single rotor without morphing is given by Equation ( 13), where the thrust coefficient called c t is a positive constant.Considering the presence of a morphing case under the actuated thigh servo motors, we reformulate Equation ( 13) as [28]: where σ i is the morphing angle of the i th rotor depicting the overlap between the propellers and the central body of the robot.c(σ i ) is a related variable coefficient called the angledependent coefficient.Notably, c(σ i = 0 • ) in this case is equal to c t .Assuming that the tips of the rotor contribute the major thrust due to the high velocity airflow [29] and each arm with the morphing servo motor is capable of rotating in the limited range of θ i1 ∈ [0, π 2 ], the angular speeds of the rotors Ω Bi need to be increased for avoiding the QQR's thrust loss and the modified rotor speeds can be defined as where f di is the desired thrust of the i th rotor, and C is a scaling factor that depends on the angle-dependent coefficient c(σ i ) and the thrust coefficient c t .In this work, the angledependent coefficient function c(σ i ) is attained via an experimental parameter identification method.If there is no overlap between the propellers and the central body of the QQR, the thrust mapping maintains unchanged.Moreover, a thigh servo motor controller with regard to transformation flight can be approximated as a first-order system with time constant k d : where l re f is related to the QQR's position under commanded aerial configurations.

High-Level Motion Control
Owing to the efficient locomotion controller, the high-level control strategy can be defined by if-else switching law for smooth and fast motion transition in this work.Hence, a unified motion control based on a full finite-state machine with various locomotion patterns is demonstrated as Algorithm 1.In order to accomplish collaborative missions, various terrestrial-aerial motions of the QQR must be dynamically and automatically transformed according to the needs of actual situations.The locomotion switching control strategies is exactly assisted by onboard stereo camera and IMU/optical flow sensors to capture the terrains and judge the QQR's center of gravity, the terrain types allow the QQR to choose aerial or terrestrial modes and the center of gravity variations help the QQR to modify the locomotion patterns.
In the quadruped parameters modulation process, the reward function r is adopted and defined as where R v and R e are the positive weight coefficients.The first term drives the QQR to move forward fast and the second term penalizes the energy consumption of the robot.The average reward of the terrestrial controller is proved to converge to a stable value in different terrain types [30].
In the quadrotor locomotion compensation process, a simplified identification mechanism based on experimental approach is developed for any combination of the rotors morphing.In this work, we used the motion capture system to capture the QQR position and computed the range of the trust loss in different aerial configurations by the drop distance of the QQR, where f t_X , f t_H , f t_T and f t_O are the collective trust of X configuration, H configuration, T configuration and O configuration, respectively.Due to the a priori knowledge of the average loss rate of desired configurations, the morphology compensation in Algorithm 1 can be quantified effectively.Compared with the standard X configuration, the loss of collective trust in the H configuration and the T configuration are smaller than the O configuration and the results of Figure 11 in the QQR shows the loss rate are f t_H / f t_X = 30%, f t_T / f t_X = 20% and f t_O / f t_X = 40% (with the error [−10%, 30%]), respectively.The stability of the aerial controller is quantified by experimental identification data.

Algorithm 1
The full finite-state machine control algorithm for terrestrial-aerial motion transition.Require: The external environmental feedback from onboard sensors.Ensure: The corresponding results of terrestrial-aerial motion.Learn the contact surfaces with the onboard IMU and optical flow sensors.Identify the bounding of target obstacles with the camera inspection.Run the X-configuration via morphing joints; Run the H-configuration via morphing joints;  The proposed unified motion control method enables the stable crawling and morphing flight of the QQR.Multimodal motions are desired to extend potential practical applications in long-endurance and specific tasks.

Results
The proposed hierarchical coordination control scheme has been evaluated on the QQR platform and numerous trials in diverse challenging environments have been executed.

Simulations
Firstly, two simulations with the terrestrial locomotion controller and the aerial locomotion controller are conducted in MATLAB R2018 to test the effectiveness of locomotion generations.

Terrestrial CPG-Based Controller Verification
A single oscillator is usually adopted to produce rhythmic control signals for a single joint.The oscillator output model ( 16) is simulated by x curves and from the results of Figure 12, a single joint signal can be modulated by the typical Hopf oscillator parameters, such as the amplitude µ, the frequency ω and the convergence rate α, β.Moreover, as can be seen from Figure 13, the stable limit cycle property of the Hopf oscillator for leg 1 with the convergence coefficient α = β = 80 was plotted, which alleviates the partial effect of external disturbances and guarantees smooth transition among the QQR's terrestrial locomotion.These simulation results show that the proposed Hopf-based locomotion controller carries the potential of simple controllability and excellent robustness on complex terrestrial terrains.The output of a single oscillator should be coupled to fit the locomotion patterns with all joints of the QQR.A good terrestrial locomotion controller is expected to generate various locomotion patterns and allows smooth transition among these locomotion patterns.Taking the walk locomotion pattern as an instance with four Hopf oscillators in thigh layer are set as amplitude = √ 10 deg, f requency = 2 Hz, phase di f f erence = 90 deg, we demonstrated the locomotion generation capability of the proposed CPG-based terrestrial locomotion controller.The output signals in the thigh layer of the CPG network model (17) are presented in Figure 14 and the result shows that the thigh layer produce four coupled curves for four joints to constitute a corresponding walk locomotion pattern cycle by cycle.Meanwhile, the shin layer also receives reflex feedback from the external environment to improve the interactive performance of this walk locomotion pattern.

Aerial Compensation Controller Verification
The aerial compensation controller enhances position tracking performances in midair when morphing occurs.Taking three romdom curves (circle, sine and helix) as reference paths, we simulated trajectory tracking comparisons to emphasize the effectiveness of the proposed compensation locomotion controller in Matlab Simulink, as shown in Figure 15.As the comparisons show in Figure 16, the QQR follows a circular trajectory at a speed of 0.5 m/s on a circle of radius 1 m, at a constant height of 1 m.The QQR changes its shape during flight switching between the X and O configurations on and off the compensation mechanism with a varying error that ranges from a minimum 0.1 m to a maximum of 0.2 m.Moreover, we also summarize some quantitative results in Figures 17 and 18, in which the QQR transformed from X configuration to H configuration with compensation shows better tracking performance on different trajectories, and the de-stability error was less than 0.1m away from the desired path.The compensation controller improves the position tracking both in hovering and forward flight.

Experiments
Several experimental setups similar to the scenes in Figure 1 have been established for practical applications, and comparative experiments have been conducted.The QQR is shown running over different surfaces, including smooth ground, rough sand, and midair.Firstly, the aerial endurance test is set as hovering on the constant height of 2.5 m, and the terrestrial endurance test is set as spinning in circles with a second timer, as shown in Figure 19.The endurance time in the aerial case is approximately 10 min and the maximum flight distance is around 2000 m, while a longer terrestrial duration is up to 30 minutes and the maximum crawling distance on the ground is about 150 m.The QQR is able to address the trade-off between flight time and crawling time on different tasks.Then, to validate the effectiveness of the proposed hierarchical coordination control scheme, we decomposed the motion skills of the QQR against dynamic circumstances into three motion primitives: crawling on smooth and rough surfaces, crossing through big and small obstacles, adapting to terrestrial-to-aerial and aerial-to-terrestrial motion transition.

Crawling on Unknown Terrains
As stated before, various terrestrial locomotion patterns suit different terrains.To demonstrate the environmental adaptability of the proposed bio-inspired terrestrial controller, we evaluated the performance with ground friction comparisons ranging from smooth surfaces to rough surfaces.From the experimental results in Figures 20 and 21 (where the red circle on the ground is regarded as the object of motion reference), the trot locomotion pattern is usually chosen on a smooth case by the QQR automatically, which can bring the fastest motion at the forward speed of 10 cm/s.However, on a rough terrain, the QQR tends to the walk locomotion pattern to enhance crawling reliability and reduce energy consumption.Both trot and the walk locomotion were demonstrated as effective locomotion patterns.Moreover, the QQR actively selected various locomotion patterns adapting to environmental uncertainties via the terrestrial locomotion controller in Figure 22.It is worth highlighting that more terrestrial locomotions can be mimicked and developed via a CPG-based controller to close the gap between the QQR and biological counterparts.

Morphing and Trajectory Tracking in Midair
Morphing flight was assessed through a series of laboratory tests, the aerial configurations are described in detail in Section 3.1.2.From the test results in Figures 23-25, the aerial compensation controller was shown to be a useful strategy to make the QQR stable under morphing situations without requiring any symmetry of the morphology.Compared with common quadrotors, the actuated thigh joints along with scalable body sizes enable the QQR to adapt to aerial obstructive environments, especially narrow spaces.A big horizontal semicircle with a diameter of 60 cm and a small rectangle with the width of 20 cm are tailored to test adaptive transformations of the robot in midair.The initial size of the QQR with standard X configuration is aforementioned 26 cm and the robot cannot cross through holes directly where the robot needs to change its configuration.
As presented in Figures 26 and 27, we also offered contrast experiments to explain the superiority of the designed morphing characteristics in the presence of obstacles and holes.It can be found that if the target hole is big enough, the QQR holds on the standard X configuration for less energy cost.Nevertheless, the QQR shrinks actively from X configuration to H configuration while detecting dimensional risks with onboard perceptions.This morphing superiority stems from our aerial locomotion controller, which compensates in advance for this morphing configuration with less performance degradation.Then, to show the effectiveness of the proposed aerial compensation controller in a real test, the following comparison between the simulation and the experiment is designed.Consider an arbitrary ellipse as a reference based on the limited experimental field: The simulation result obtained from the Section 4.1.2and the real test result carried on Figure 28 are shown as Figures 29 and 30. Figure 29 shows the tracking performance of the proposed controller, and the maximum error is less than 0.1 m, even in real turning situations.Moreover, the tracking error on X axis, Y axis and Z axis are given in Figure 30, respectively.The results imply that even though the aerial configuration of the QQR are changed and the stochastic disturbances are considered, the proposed aerial compensation controller can ensure the stability of the QQR system.

Adaptive Terrestrial-Aerial Motion Transition
Based on an open-source simulation scenario in Figure 31, we quickly checked the proposed coordinated motion control scheme at an extreme environment of a whole linear obstacle.From motion transition simulations, the QQR tries to utilize terrestrial locomotion patterns for lower energy consumption until the robot cannot search feasible spaces on the ground.Then, the QQR transfers to aerial locomotion patterns for traverse maneuverability.
The feasibility of seamless motion transition has been further validated experimentally.As shown in Figure 32, adaptive motion transition from crawling mode to flight mode permits the QQR to traverse a complex occasion successfully: the QQR firstly passes through a board obstacle via a legged mode with a much better power utilization; then, the QQR automatically starts the propellers to cross a box obstacle at a higher rate.Hence, the QQR is able to overcome any obstacle using the transition between its terrestrial and aerial modes and the terrestrial-aerial adaptivity shows the superiority of optimal motion allocation in intricate missions.The terrestrial-aerial motion transition strategy endows robots with hybrid intelligence and brings tremendous application prospects on collaborative missions.
Using only the proposed control algorithm, numerous experimental results show that the QQR can adapt to various typical terrains based on the optimized mechanical design.Compared with other hybrid robots shown in Table 4, the QQR exhibits perfect multimodality and maneuverability under the limited mass and dimension in complicated terrains.

Conclusions
The proposed terrestrial-aerial QQR enables both a crawling motion and morphing flight with the compact and agile combination of the high maneuverability in quadrotors and the long endurance in quadrupeds.In this paper, based on the design of a micro, lightweight, and integrated quadruped-quadrotor robot, a hierarchical coordination control scheme consisting of a low-level locomotion control and a high-level motion control is explored.Hence, the adaptive robot system achieves both terrestrial locomotion and aerial locomotion as well as smooth terrestrial-to-aerial or aerial-to-terrestrial transition in unknown and uneven environments.Several crucial trials are performed to validate the proposed QQR's agility and adaptability.In the future, we will focus on enhancing autonomous navigation for intelligent motion planning with onboard sensing and computing in demanding terrains and real-world application scenarios [33][34][35][36][37].

Figure 1 .
Figure 1.Passing through a narrow gap to enter a collapsed building.

Figure 2 .
Figure 2. The mechanical configuration of the QQR prototype.

Figure 3 .
Figure 3. Illustration of the QQR.{O B } is the body coordinate.

Figure 4 .
Figure 4. Coordinates at different joints on each leg (such as leg 2).{O G } denotes the global coordinate, and {O ik } (i = 1, 2, 3, 4 and k = 0, 1, 2) represents the floating frame whose origin fixed on the joints or the foot tip.l ix , l iy are defined as the position of the thigh motors in the body frame O B and l ix = 5, l iy = 5 in the QQR case.l 1 , l 2 are the length of the leg links, respectively.θ i1 , θ i2 are the rotational joint angles around Z i0 , Z i1 axis of the leg and 2φ 2 = 180 − 2φ 1 .Because the symmetrical design of the QQR's thigh motors in the diagonal of the square body, φ 1 = φ 2 = 45 deg in this work.
the current position of the i th rotor in the body coordinate O B which changes in different aerial locomotion patterns and L = [L 1 , L 2 , L 3 , L 4 ].

Figure 5 .
Figure 5. Diagram of the proposed hierarchical coordination control scheme.

Figure 8 .
Figure 8.The topology network of the proposed terrestrial CPG.

3. 1 . 2 .
Aerial Locomotion Control The QQR utilizes four servo motors in the thigh layer for aerial morphing locomotion in midair, including X configuration for efficient flight, H configuration for narrow flight (with the angle change of 45 deg in four thigh joints conversely), T configuration for closer inspection (with the angle change of 45 deg in anterior two thigh joints), and O configuration for omnidirectional shrinking (with the angle change of 45 deg in four thigh joints uniformly), as shown in Figure 9.

Figure 10 .
Figure 10.The aerial locomotion control diagram for the QQR's morphing.

6 : 7 : 8 :
while locomotion_mode == mode::trans_Walk() do Run the walk locomotion pattern via CPG; Compute the cost function to modulate the terrestrial motion control performance; while locomotion_mode == mode::trans_Trot() do 11: Run the trot locomotion pattern via CPG; 12: Compute the cost function to modulate the terrestrial motion control performance; while locomotion_mode == mode::trans_Translation() do 15: Run the translation locomotion pattern via CPG; 16:Compute the cost function to modulate the terrestrial motion control performance;

Figure 11 .
Figure 11.The collective trust loss in different aerial configurations.

Figure 12 .
Figure 12.Effect of oscillator parameters on a single joint output.

Figure 13 .
Figure 13.The limit cycle property of the proposed Hopf oscillator.

4 Figure 14 .
Figure 14.The walk locomotion pattern generated by the proposed CPG-based terrestrial controller.

Figure 16 .Figure 17 .
Figure 16.The tracking performance of the proposed aerial compensation controller on a circle path.

Figure 18 .
Figure 18.The tracking performance of the proposed aerial compensation controller on a helix path.

Figure 19 .
Figure 19.The practical endurance tests on the QQR.

Figure 22 .
Figure 22.QQR crawls on complex terrains with the adaptive locomotion.

Figure 23 .
Figure 23.X configuration to H configuration for narrow spaces.

Figure 24 .
Figure 24.X configuration to T configuration for closer inspection.

Figure 25 .
Figure 25.X configuration to O configuration for omnidirectional shrinking.

Figure 26 .
Figure 26.A big target beyond the size of the QQR (without morphing).

Figure 27 .
Figure 27.A small target below the size of the QQR (with morphing).

Figure 31 .
Figure 31.Simulations in a cluttered case with terrestrial-aerial motion transition.

Figure 32 .
Figure 32.Experiments in a complex occasion with terrestrial-aerial motion transition.

Table 1 .
Technical specifications of the QQR prototype.
1: Unlock the QQR and checking the robot state for safe operation; 2: Detect control commands FSM_State based on the stereo camera, i.e.

Table 4 .
Characteristics of the hybrid robots.