Next Article in Journal
Terrestrial Megafauna Response to Drone Noise Levels in Ex Situ Areas
Previous Article in Journal
Drones and Blockchain Integration to Manage Forest Fires in Remote Regions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Guidance, Navigation and Control System for Multi-Robot Network in Monitoring and Inspection Operations

by
Mohammad Hayajneh
1,* and
Ahmad Al Mahasneh
2
1
Department of Mechatronics Engineering, Faculty of Engineering, The Hashemite University, Zarqa 13133, Jordan
2
Department of Mechatronics Engineering, Faculty of Engineering, Philadelphia University, Amman 19392, Jordan
*
Author to whom correspondence should be addressed.
Submission received: 28 September 2022 / Revised: 22 October 2022 / Accepted: 27 October 2022 / Published: 30 October 2022

Abstract

:
This work focuses on the challenges associated with autonomous robot guidance, navigation, and control in multi-robot systems. This study provides an affordable solution by utilizing a group of small unmanned ground vehicles and quadrotors that collaborate on monitoring and inspection missions. The proposed system utilizes a potential fields path planning algorithm to allow a robot to track a moving target while avoiding obstacles in a dynamic environment. To achieve the required performance and provide robust tracking against wind disturbances, a backstepping controller is used to solve the essential stability problem and ensure that each robot follows the specified path asymptotically. Furthermore, the performance is also compared with a proportional-integral-derivative (PID) controller to ensure the superiority of the control system. The system combines a low-cost inertial measurement unit (IMU), a GNSS receiver, and a barometer for UAVs to generate a navigation solution (position, velocity, and attitude estimations), which is then used in the guidance and control algorithms. A similar solution is used for UGVs by integrating the IMU, a GNSS receiver, and encoders. Non-linear complementary filters integrate the measurements in the navigation system to produce high bandwidth estimates of the state of each robotic platform. Experimental results of several scenarios are discussed to prove the effectiveness of the approach.

1. Introduction

Today, unmanned vehicles (UVs) are used to carry out a number of tasks in order to automate processes, complete tasks faster and more accurately, and ensure the safety of humans [1]. Moreover, UVs are often used in smart cities [2,3]. Due to this, researchers have developed ground and aerial robots that were designed to operate under uncertain conditions. Often, it is a burden for a single robot to accomplish a task effectively [4]. For this reason, a group of small, inexpensive robots can work together in a network, each with its own capabilities, to accomplish a common goal [5].
The applications of robotic systems in clean energy systems are continuously developing since there is a global tendency to shift to renewable energy sources such as solar, wind, and hydro-energy systems [6]. However, for these systems to perform well, they must be inspected and monitored on a regular basis [7]. For instance, in [8], a drone with thermal imaging is used for solar panels defect detection. In addition, a climbing robot was utilized to inspect and maintain a solar power plant [9]. To maximize the efficiency of such operations, several studies used multi-robot networks to develop technology to analyze the normal operation and failure of solar modules by attaching optical and thermal infrared sensors to unmanned aerial vehicles (UAVs) [10,11]. A similar study used multi-drone equipped with vision and LiDAR sensors for global inspection, guiding climbing robots in analyzing structure parts using non-destructive inspection methods and locally repairing smaller defects in wind turbines [12]. With the advancement of IoT technology, the use of multi-robot systems to collect data from a network of wireless sensors has received increased attention and has grown dramatically in applications such as agriculture [4], railway track monitoring [13], and search and rescue [14].
To achieve a trusted level of autonomy in multi-robot system operations, it should include the following major components: guidance, navigation, and control [15]. The guidance system is in charge of providing the robot with a planned trajectory that allows it to travel from its starting point to the desired location while avoiding obstacles and tracking a target. In a multi-robot system, each agent must cooperate and coordinate with the others in order to achieve real-time cooperative navigation and free-collision path planning [5,16,17]. The navigation system is used to determine the location, orientation, and velocities of each robot in the multi-robot system at any desired time [18,19,20]. The control system provides the correct forces or torques to achieve the guidance goals. Backstepping (BS), sliding mode control (SMC), feedback linearization, PID, optimal and robust control, learning-based control, and other techniques were used to solve the stabilization control problem and trajectory control for a quadrotor and a mobile robot [21,22]. Backstepping control is considered one of the well-developed nonlinear control approaches that can stabilize dynamic systems while handling uncertainties to achieve the required performance [23,24].
In recent years, swarming robot dynamics and control have been active research topics. The difficulty of validating these algorithms in actual tests has been a barrier to more frequent and widespread use. Many researchers focused their efforts on developing a stable guidance, navigation and control (GNC) system to address this issue. In [25], robust navigation algorithms for multi-agent fixed-wing aircraft are presented. These algorithms are based on adaptive moving mesh partial differential equations controlled by the free energy heat flow equation, and they are experimentally validated using LQR controllers and multi-scale moving point guidance. A similar project used a PID controller and a waypoint-based guidance system to evaluate an INS/ GNSS navigation system [26]. An affordable GNC solution with a visual-based navigation system was tested on a group of small unmanned ground vehicles and quadrotors to achieve a common goal [27]. In a practical marine environment, a hybrid framework for guidance and navigation of swarms of unmanned surface vehicles (USVs) that combines two layers of offline planning and online planning was applied [28].
This work developed a complete system architecture for a guidance, navigation, and control solution to enable small UAVs and UGVs in monitoring and inspection applications. To solve the essential stability problem and ensure that each robot follows the specified path asymptotically, a backstepping controller is used and discussed in detail. A guidance algorithm for generating a flight trajectory based on a potential field method is also described, allowing a robot to track a moving target and avoid obstacles in a dynamic environment. To generate a navigation solution (position, velocity, and attitude (orientation) estimations), the system combines a low-cost inertial measurement unit (IMU), a GNSS receiver, and a barometer for UAVs, which is then used in the guidance and control algorithms. For UGVs, a similar solution is used by integrating IMU, a GNSS receiver, and encoders. The measurements in the navigation system are integrated by non-linear complementary filters that produce high-bandwidth estimates of the state of each robotic platform. What distinguishes the developed system is its direct implementation of its algorithms in embedded systems with minimal programming cost. This work provides users with a reliable and viable GNC structure, which is supported by several practical experiments conducted both internally and externally for life applications. The developed GNC system was tested using a true reference motion capture system to demonstrate the efficiency of the integrated system in terms of stability, accuracy, and maneuverability. Several outdoor tasks were carried out to evaluate the monitoring and inspection process, rather than a single robot. These operations demonstrate the navigation system’s ability to precisely locate the aerial and ground robots.
The current study is structured as follows: Section 2 provides an overview of the dynamic model of each robot. Section 3 discusses the adopted backstepping controller for a quadrotor and a ground robot. A motion planning algorithm is developed and discussed in Section 4. The navigation system for each robotic platform is adopted for accurate state estimations. Section 6 illustrates indoor and outdoor experiments of multi-robot missions. Furthermore, this section discusses the experimental results. Finally, concluding remarks and future research directions are summarized in Section 7.

2. System Modeling and Description

The kinematics and dynamics models of the quadrotor and wheeled robot, in addition to their relevance in the design of robust controllers, can be helpful in understanding the stability nature and operations under uncertain conditions.

2.1. Quadrotor Model

There has been much explanation and study of the mathematical model of the quadrotor in the literature [24,29,30], and we will summarize the general model here. It is established through:
v ˙ = 1 m ( T + f d + f w ) + g R T u 3 ω × v p ˙ = R v ω ˙ = J 1 [ τ S ( ω ) J ω ] R ˙ = R S ( ω )
As shown in Figure 1, the vectors of linear velocity and position are determined by v = [ u , v , w ] T and p = [ x , y , z ] T , respectively. Moreover, ω = [ ρ , q , r ] T is the angular velocity vector. Any vector z in body frame B can be translated into inertial frame I by a rotation matrix R , and any vector z can be represented by the skew symmetric matrix S ( z ) . The gravity constant g is represented in inertial frame by the vector u 3 = [ 0 , 0 , 1 ] T , and the moment of inertia is given by the matrix J and the mass by m. The terms ω × v , f d , and f w in Equation (1) are the gyroscopic effects, the drag forces, and the wind disturbances, respectively. Furthermore, the thrust force T and moments τ on each body axis are expressed as follows:
T = 0 0 c t ( ω m 1 2 + ω m 2 2 + ω m 3 2 + ω m 4 2 ) T τ = d c t ( ω m 1 2 + ω m 2 2 + ω m 3 2 ω m 4 2 ) d c t ( ω m 1 2 + ω m 2 2 ω m 3 2 ω m 4 2 ) c q ( ω m 1 2 ω m 2 2 + ω m 3 2 ω m 4 2 ) + J m q π 30 ( ω m 1 + ω m 2 ω m 3 + ω m 4 ) J m p π 30 ( ω m 1 ω m 2 + ω m 3 ω m 4 ) 0
The torque and thrust coefficients of the rotor (motor+propeller) are c q and c t , respectively; J m is the rotor’s inertia and d is the distance between the center of the rotor to the quadrotor’s center of mass.

2.2. Wheeled Robot Model

In the literature, the kinematic model of two-wheeled differential-drive robots is subjected to non-holonomic constraints, assuming that the wheels do not skid [24,31]. In this case, the robot motion is described as follows:
v = R ( ψ ) x ˙ y ˙ ψ ˙ T = r ω R + ω L 2 ω = r ω R ω L d R ( ψ ) = cos ψ sin ψ 0 sin ψ cos ψ 0 0 0 1
In the above equation, ψ and ω represent the angle of rotation and angular velocity around the robot’s z-axis, respectively, and v is the linear velocity; r and d are the radius of the robot’s wheel and the distance between the robot’s rear wheels, respectively; ω R and ω L are the angular speeds of the right and the left wheels, respectively. In Figure 1, the body frame of the robot moves relative to an inertial frame O. Therefore, the kinematic equations in Cartesian coordinates are given as follows:
x ˙ = v cos ψ y ˙ = v sin ψ ψ ˙ = ω

3. Backstepping Controller

3.1. Quadrotor Control

The backstepping (BS) controller is used in this work because of its capacity to solve the quadrotor’s essential stability problem. Most crucially, the control system ensures that the quadrotor’s position follows the specified path asymptotically [23,24]. The BS control system integrates both translation and rotation dynamics within a single control.
The adopted control law generates one thrust T and three torques τ = [ τ ϕ , τ θ , τ ψ ] T and a timing law for the path parameter γ ( t ) that guarantees convergence of the quadrotor’s position P ( t ) to the desired path P d ( γ ) , which is a function of a virtual arc γ . The control structures for the thrust and the torque vector are given by Equations (5) and (6), respectively:
T d = m | | k 1 2 k 2 e 2 + k 1 σ ˙ ( e 1 + g u 3 p ¨ d | | = m | | t d | | , T d > 0 r 3 d = t d | | t d | | T = r 3 d T r 3 T d
where k 1 and k 2 are positive gains, and m is the mass of the quadrotor; r 3 is the third column of the rotation matrix ( R ), and it indicates the direction of the quadrotor’s z body-axis. As a result, if r 3 equals the desired thrust direction ( r 3 d ), the thrust force ( T ) will be equal to the desired thrust force ( T d ). Based on the desired thrust value ( T d ) and thrust direction ( r 3 d ) in Equation (5), the thrust control law (T) is generated. Furthermore, the torque control law ( τ ( t ) ) is calculated in Equation (6) according to the angular acceleration ( ω ˙ c ). Moreover, an arbitrary function ( ω ˙ 3 c ) is used to define the dynamics of yaw angle.
ω ˙ 3 c = l 2 ( ω 3 ω 3 d + l 1 r 2 d T r 1 ) + ω ˙ 3 d l 1 d d t ( r 2 d T r 1 ) ω ˙ c = S ( u 3 ) ( k 4 e 4 R T r 3 d + k 3 S ( u 3 ) 2 ( R ˙ T r 3 d + R T r ˙ 3 d ) + 1 m k 1 S ( u 3 ) 2 ( T ˙ d R T e 2 + T d R T e ˙ 2 ) ) + R ˙ T R d ω d + R T R d ω ˙ d + [ 0 0 ω ˙ 3 c ] T τ = J ω ˙ c + S ( ω ) J ω
where l 1 , l 2 , k 3 , and k 4 are positive gains; r 1 and r 2 are the first and second columns of the rotation matrix ( R ), respectively. The skew symmetric matrix is represented by S ( . ) ; R d and ω d denote the desired rotation matrix and quadrotor’s desired angular speed, respectively. Based on the BS controller, four error vectors ( e 1 to e 4 ) are defined. Those error vectors and derivatives are expressed as follows:
e 1 = p p d ( γ ) , e ˙ 1 = R v p ˙ d , e ¨ 1 = T m R u 3 + g u 3 p ¨ d e 2 = σ ( e 1 ) + 1 k 1 e ˙ 1 , e ˙ 2 = σ ˙ ( e 1 ) + 1 k 1 e ¨ 1 e 3 = r 3 r 3 d e 4 = k 3 S ( u 3 ) 2 R T r 3 d + S ( u 3 ) ( ω R T R d ω d ) T d m k 1 S ( u 3 ) 2 R T e 2
where p d and v d are the desired positions and velocities of the quadrotor, respectively. A time-dependent position error vector ( e 1 ) can influence the initial position of the vehicle significantly. As a result, the control structure is distinguished by the fact that it is not directly related to the time factor associated with the first error. In other words, thrust (T) and torque vector ( τ ) are not directly linked to ( e 1 ). In addition to the position error, the other errors are regulated by the saturation function ( σ ) of ( e 1 ), which allows for bounded actuation. A detailed discussion of the stability proofs for these control laws can be found in [23].
It is essential for the control algorithm to keep a well-defined control law to ensure that the vehicle follows the path. In order to achieve this condition, we have a timing law, which is obtained by adding the second time-derivative of the path parameter γ :
γ ¨ = k γ σ ( γ ˙ γ ˙ d ) + u 1 T R T ( k 1 2 k 2 e 2 + k 1 σ ˙ ( e 1 ) γ ˙ p d ) p ¨ d = p d γ ˙ 2 + p d γ ¨
where p d and p d are the first and second partial derivatives of the desired path, respectively. The control law in Equation (8) drives γ ˙ to γ ˙ d , which is calculated based on the path’s geometry. Any vector in tangent frame { T } can be represented in inertial frame { I } by the rotation matrix R T at any points where p 0 . For ease of use, most variables, which were previously defined, are listed in Table 1.

3.2. Ground Robot Control

A backstepping controller is adopted to accurately regulate the actuators that drive the robot from its current position p = [ x y ψ ] T to the desired position p d = [ x d y d ψ d ] T . The aim of the control is to minimize the following error:
ϵ = R ( ψ ) ( p d p ) = ( x d x ) cos ψ + ( y d y ) sin ψ ( x d x ) sin ψ + ( y d y ) cos ψ ψ d ψ
The error derivatives after substitutions from Equations (4) and (9) can be simplified as follows:
ϵ ˙ = ω ϵ 2 v + v d cos ϵ 3 ω ϵ 1 + v d sin ϵ 3 ω d ω
As discussed in [24], two control inputs (i.e., v c and ω c ) are considered. Hence, the control law is implemented as follows:
v c ω c = k 1 ϵ 1 ( F v d ) + v d cos ϵ 3 ω d ( τ cos ϵ 3 )
where v d and ω d are the desired linear velocity and angular velocity of the robot, respectively. The variables F and τ are calculated as follows:
F = v d k 5 + 1 τ = k 1 ( δ sin ϵ 3 δ ϖ ) + η σ ϖ = η ( k 3 F ϵ 1 3 + k 2 v d ϵ 1 2 ϵ 2 + ( k 2 F 2 + k 5 2 F v d 2 ) ϵ 1 ) δ = v d + k 5 F v d η = v d v d + k 5 F σ = δ ( k 2 k 5 v d 3 k 3 F ) ϵ 1 3 + 2 δ k 2 k 5 v d ϵ 2 2 ϵ 1 + k 3 δ ( 2 F 3 k 5 F ) ϵ 1 2 ϵ 2 + δ k 2 v d 2 ϵ 1 2 sin ϵ 3 + δ k 3 ( F 3 v d 4 ) ϵ 1 + δ k 2 k 5 ( v d F + F 2 v d 3 ) ϵ 2
where k 1 , k 2 , k 3 , k 4 , and k 5 are the controller gains.

4. Motion Planning for Multi-Robot System

The purpose of this section is to introduce a new UV motion planning system based on potential forces. This system helps the robot (i.e., ground and aerial robots) follow a moving target within a dynamic environment containing obstacles. To achieve the desired target, an individual trajectory is generated for each robot, which has access to all the state information in one processor. In Section 3, the proposed control structure for each robot ensures that the tracking error will converge to zero, where tracking error of the ith robot is measured as:
e t i = q d i q i
where q d i is a vector of the robot desired trajectory, and q i is its position with respect the inertial frame. The robot will continue to move towards its destination as long as there are no obstacles in its path. Alternately, the robot’s predefined path must be modified by feeding small variations in its motions. Essentially, the proposed approach generates a dynamic trajectory that allows the robot to track the target and avoid obstacles by generating attractive and repulsive potential forces.

4.1. Potential Attractive Force

The attractive force model is defined by the relative position and velocity vectors of two robots as follows:
f A = A 1 ( p ) + A 2 ( v ) A 1 ( p ) = 2 k p δ p A 2 ( v ) = 2 k v δ v
The attractive force model includes two force components, A 1 and A 2 , as shown in Figure 2. The first component, A 1 , propels the robot toward the target, shortening the distance between them. The second component, A 2 , keeps the robot moves at the same velocity as the target. The relative position between the robot and the target is given by δ p , and the relative velocity is given by δ v . Each force component can be tuned correctly by the gains k p and k v . Figure 2 illustrates the relationship between the attractive force components and the relative position and velocity of the robot to the target in a 2-dimensional space. Normally, this system is applied to a ground robot that is following another robot, and the concept can also be extended to cover UAVs in the third dimension. The attractive force for each motion direction (i.e., x, y, and z) can be given as follows:
f x = 2 k p x δ x + 2 k v x δ v x f y = 2 k p y δ y + 2 k v y δ v y f z = 2 k p z δ z + 2 k v z δ v z

4.2. Potential Repulsive Force

When a robot encounters an obstacle, it will also be provided with an extended potential force to detour and avoid collision. Considering p o b s and v o b s are the position and the velocity of an obstacle which can be obtained on-line, then the relative velocity between the robot and the obstacle is ϵ v = v v o b s , and the relative position is ϵ p = p o b s p , as shown in Figure 3. Then, the value of ϵ v in the direction from the robot to the obstacle is given by:
v R O = ϵ v T n R O n R O = ϵ p | | ϵ p | |
where n R O is a unit vector pointing from the robot to the obstacle. As shown in Figure 3, in the case where the robot is moving away from the obstacle as v R O 0 , no avoidance action is needed from the robot. In contrast, if v R O > 0 , the robot requires a force to push it away from the obstacle. As a result, a model is used to generate appropriate repulsive forces. These forces are computed using the following function:
r = r 1 + r 2 . r 1 = η ( ρ s ρ m ) 2 ( 1 + v R O a m a x ) n R O r 2 = η v R O v R O ρ s a m a x ( ρ s ρ m ) 2 n R O
It can be seen in Figure 3 that r 1 is a force in the opposite direction of v R O n R O , which keeps the robot away from the obstacle. Moreover, r 2 is in the same direction of the perpendicular vector v R O n R O and acts as a steering force. Moreover, η is a positive constant, ρ s is the shortest distance between the center of the robot and the center of the obstacle, and ρ m is the distance that the robot travels before it approaches zero, which is a function of a maximum deceleration a m a x and is given as follows:
ρ m = v R O 2 2 a m a x .
In order to make the repulsive force model valid for multiple obstacles n o b s , the following relation is established:
r = n = 1 n o b s r i .
Using Equations (14) and (17), the robot is subjected to an overall virtual force equal to the sum of the calculated attractive and repulsive forces as follows:
f t o t = f + r .

5. Navigation System and State Estimations

In real-time platforms, an estimation system is essential to compensate for noise and biases in sensors, and, in addition, to provide accurate attitude, position, and velocity information to the aerial and ground robots. One of the navigation solutions can be implemented with a typical inertial measurement unit (IMU) and a low-cost GNSS receiver for outdoor use or tracking systems for indoor use. This section discusses the design and implementation of attitude, position, and velocity estimations for both the ground robot and the quadrotor, addresses the lateral drifting, and improves the vertical position of the quadrotor platform.

5.1. Quadrotor State Estimations

Using a complementary filter, the attitude measurement is performed by fusing high-frequency gyroscope measurements with low-frequency magnetometer and accelerometer measurements. For the sake of direct implementation in embedded systems, the discrete quaternion representation of the attitude estimation is given as follows [32]:
q ^ k + 1 b ^ k + 1 = I 4 T 2 Ξ ( q k ) 0 I 3 q ^ k b ^ k + 1 2 Ξ ( q k ) k 1 q k 2 q i = 1 n ( y i k × y ^ i k )
where T as the sampling time interval and the sub-index k abbreviates the time instant t = k T ; q ^ and b ^ are the estimated attitude in quaternion and the estimated gyro bias, respectively. In addition, the gains k 1 q and k 2 q are positive and tuned for best performance; I 3 and I 4 are 3 × 3 and 4 × 4 identity matrix, respectively. Furthermore, the term ( y i k × y ^ i k ) represents the variation between the estimated inertial vector and the measured one in the body frame. In practice, two vectors are utilized, which are the earth’s gravity vector and the earth magnetic field vector, and can be measured in body frame by an accelerometer ( a ) and a magnetometer ( m ). Hence the term i = 1 n ( y i k × y ^ i k ) can be implemented as follows:
i = 1 2 ( y i k × y ^ i k ) = k a ( a k × R ( q k ) k [ 0 0 g ] T ) + k m ( m k × R ( q k ) [ m x o m y o m z o ] T )
where the vectors [ 0 0 g ] T and [ m x o m y o m z o ] T are the earth’s gravity and magnetic field vectors, respectively, which are measured in inertial frame; R ( q k ) is the quaternion rotation matrix.
This work also implements another complementary filter for accurate position and velocity estimations. Due to the low sampling rate and lack of accuracy, standalone GNSS systems are not reliable when they are used in high-dynamic systems such as quadrotors. Thus, GNSS data is integrated with accelerometer measurements by a complementary filter to provide continuous position and velocity estimations. For more reliable estimations in altitude, a barometer is used to provide the measurements of height. The position and velocity complementary filter is provided in discrete-time form as follows:
p ^ k + 1 v ^ k + 1 b ^ a k + 1 = I T I T 2 2 R ( q k ) 0 I T R ( q k ) 0 0 I p ^ k v ^ k b ^ a k + T 2 2 I T I 0 R ( q k ) ( a m k + [ 0 0 g ] T ) + k p p ˜ k k v v ˜ k k b v ˜ k
where p ˜ and v ˜ are the errors in position and velocity, respectively, that represent the differences between the measured values (i.e., p and v ) and estimated ones (i.e., p ^ and v ^ ) at time point k; b ^ a is the estimated accelerometer bias, and k p , k v , and k b are positive gains.

5.2. Ground Robot State Estimation

In a mobile robot navigating through an indoor environment, only the state in plane (i.e., x, y, and ψ ) is needed assuming that the robot moves on a horizontal and flat surface. However, a robot’s motion could be significantly influenced by terrain characteristics when it is operating in an outdoor environment. Therefore, it can be oriented around the three axes of its reference frame in any way, resulting in six degrees of freedom (i.e., three positions and three orientations). Through a fusion of IMU, GNSS, and compass data, the work utilized the same technique as the quadrotor, which is discussed in the previous sub-section, to determine the localization and orientation of the ground robot. To enhance the estimations, the positions and orientation of the wheeled mobile robot is usually calculated based on an odometry according to measurements from motor encoders as follows [31]:
x k + 1 = x k + δ s k c o s ( ψ k + δ ψ k 2 ) y k + 1 = y k + δ s k s i n ( ψ k + δ ψ k 2 ) ψ k + 1 = ψ k + δ ψ k
At each sampling time k, δ s k and δ ψ k are the increments in distance and angle.

6. Experimental Setup and Results

In this section, the proposed guidance, navigation, and control (GNC) system is implemented and evaluated by a multi-robot system. The performance of the adopted GNC system was evaluated in two different environments. The first environment was indoor and included three robotics platforms (i.e., two quadrotors and one ground robot), as shown in Figure 4. The second environment was implemented outdoor with another three drones and two car-like robots, as shown in Figure 4. An indoor system demonstrated robotic platforms with a motion capture system that provided real-time position and orientation. Based on several retro-reflective markers mounted on each robot, the motion capture system ran a software that defined rigid bodies. Each quadrotor executed the control algorithms with an onboard Intel Aero computer and the ground robot with a Raspberry Pi 3. In this study, quadrotor drones and ground robots were used to evaluate the performance of the backstepping controllers with the presence of disturbances. In addition in this stage, the drone’s ability to track a ground robot and avoid obstacles was tested in several experiments. Moreover, multi-robot missions were performed outdoors in different scenarios for surveying and monitoring, as illustrated in Figure 5. A Pixhawk autopilot was used to test the proposed control and navigation algorithms. The guidance system and communication between robots were built on a Raspberry Pi using the Robotics Operating System (ROS).

6.1. Backstepping and PID Comparison

Two experiments were implemented on each robot platform (i.e., quadrotor and ground robot). Well tuned PID controllers were adopted in the first experiment and compared with a similar scenario in the second experiment in which backstepping controllers were implemented, as discussed in Section 3. The performance of each platform was analyzed by computing the mean squared error (MSE) in positions and velocities. The results are listed in Table 2 for a quadrotor and Table 3 for a ground robot. The tables show the superiority of the backstepping controlled quadrotor and ground robot over PID controllers, particularly in high-speed paths. Figure 6 illustrates how the BS-controlled quadrotor was correctly tracking the predefined circular-path, whereas the PID-controlled quadrotor was diverging from the high-speed path. A similar performance was demonstrated by a BS-controlled ground robot on the squared-path. Figure 7 and Figure 8 clearly show the better performance of the backstepping controller for the quadrotor and the ground robot, respectively.
Figure 9 shows the performance of quadrotors using both PID and BS controllers in an experiment on a squared-path with the addition of fan air flow as a disturbance, located at ( x = 0 , y = 1.5 m) and directed toward the y axis. The PID-controlled quadrotor’s highest deviation from the reference route as a result of the disturbance was 42.4 cm, whereas the BS-controlled quadrotor’s highest position divergence was 5 cm. Therefore, based on the BS controller’s good results, it was implemented in ground robots and quadrotors for the best stability and tracking in multi-robot missions.

6.2. GNC System Performance

The backstepping controller was used along with the navigation system in Section 5 and the guidance system in Section 4 for a full guidance, navigation, and control (GNC) structure to demonstrate the performance of the quadrotors and ground robot in a network of multi-robots.
An indoor experiment was conducted to demonstrate a mission with two drones and one ground robot in order to to validate the overall GNC system. In this mission, a main drone was assigned to follow a ground robot that was traveling in a 4 × 4 squared path and followed with a 1.2 radius circle. A traveling drone was used to travel to point (−1.5, −1.5, 1.2) m and hover for a while. This point is one of the squared path’s corners that the tracking drone will pass through. Figure 10 depicts how the tracking drone avoided colliding with the traveling drone. Furthermore, it avoided the traveling drone once more while the last was approaching it, as shown in Figure 10.
The position and velocity of the two drones and the ground robot are illustrated in Figure 11. The figure shows that the tracking drone was precisely following the ground robot (i.e., M S E x y = 0.173 m, M S E z = 0.046 m, and M S E v x y = 0.125 m/s), despite the fact that the obstacle avoidance task increased the error between them. The figure also shows how the tracking drone safely avoided the traveling drone on the 20th and 38th seconds and returned back to the desired path.

6.3. Outdoor Missions

Different missions were carried out in an outdoor environment to demonstrate the performance of a network of several drones and ground robots. The first mission scenario was carried out with three drones, one of which was assigned to lead the others on a predefined path, as illustrated in Figure 12. In this mission, three drones were launched from the same location, and the leading drone was tasked with following the predefined way-points while the other two drones followed it safely. Figure 13 illustrates the performance of the 140-s mission. The average lateral speed (i.e., in the xy direction) of each drone was approximately 10 m/s at different heights. The three drones completed the mission precisely and without any collisions, with the follower drones maintaining a nearly constant distance from the leading drone throughout the mission.
In the second scenario, the three drones were programmed to survey an area of approximately 87,000 m 2 at the same time. Each drone was assigned to survey a portion of the area and photograph various locations, as shown in Figure 14. The mission was completed in approximately 160 s, with 24 photos of the entire area being taken, as illustrated in Figure 15. For sake of comparison, the same task was assigned to a single drone. The mission took approximately 330 s to complete using one drone.
The final outdoor scenario involved two ground robots performing a monitoring task. As shown in Figure 16, each robot was assigned to follow a different path on campus at the same time. As shown in Figure 17, each robot was moving at a speed of 5 m/s along a path that was approximately 1770 m long for the first robot and 2390 m for the second robot. As a result, the first ground robot completed its assignment in 390 seconds, whereas the second ground robot completed the mission in 490 s, taking into account the time it took to launch to the first way-point and return to the starting point.
The goal of the last two tasks was to evaluate the monitoring and inspection process using a group of aerial and/or ground robots rather than a single robot. The previous operations demonstrated that using a group of robots speeds up and improves the efficiency of the work. The robots demonstrated the ability to track way-points without colliding or intersecting with one another. These two experiments also demonstrated the ability of the navigation system to precisely locate the aerial and ground robots.

7. Conclusions

This work demonstrated the practical implementation of light-computing nonlinear algorithms for a full guidance, navigation, and control system for a network of quadrotor drones and ground robots. For the best path tracking and position and orientation stability, a backstepping controller was used. In order to implement a robust navigation system for control and guidance systems, nonlinear complementary filters were used to accurately estimate attitudes, positions, and velocities. For optimal robot guidance, a potential field method was used for tracking and obstacle-free path planning. Several indoor and outdoor experiments were conducted to assess the effectiveness of the GNC system on multi-robot missions. The results demonstrated the advantages of the proposed system in monitoring and inspection applications. This work paves the way for future efforts to improve the navigation system for higher-level perception in order to recognize and classify detected objects/events and infer some of their attributes. Another interesting aspect of the research is the use of the proposed GNC on the scales of swarms of quadrotors and ground robots in a challenging environment.

Author Contributions

Methodology, M.H.; validation, M.H. and A.A.M.; resources, M.H.; writing—original draft preparation, A.A.M.; writing—review and editing, M.H. and A.A.M.; supervision, M.H.; funding acquisition, M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by funding from the Royal Academy of Engineering, Transforming Systems through Partnership TSP1040.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research is carried out in the context of the project “Drone-Assisted Micro Irrigation for Dry Lands in Jordan Using IoT Enabled Sensor Network (DAMIJO-IoT). All experimental results were obtained using all of the facilities provided by the Faculty of Engineering at the Hashemite University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. BaniHani, S.; Hayajneh, M.R.M.; Al-Jarrah, A.; Mutawe, S. New control approaches for trajectory tracking and motion planning of unmanned tracked robot. Adv. Electr. Electron. Eng. 2021, 19, 42–56. [Google Scholar] [CrossRef]
  2. Li, P.; Wu, X.; Shen, W.; Tong, W.; Guo, S. Collaboration of heterogeneous unmanned vehicles for smart cities. IEEE Netw. 2019, 33, 133–137. [Google Scholar] [CrossRef]
  3. Mohamed, N.; Al-Jaroodi, J.; Jawhar, I.; Idries, A.; Mohammed, F. Unmanned aerial vehicles applications in future smart cities. Technol. Forecast. Soc. Chang. 2020, 153, 119293. [Google Scholar] [CrossRef]
  4. Ju, C.; Son, H.I. Multiple UAV systems for agricultural applications: Control, implementation, and evaluation. Electronics 2018, 7, 162. [Google Scholar] [CrossRef] [Green Version]
  5. Mutawe, S.; Hayajneh, M.; BaniHani, S.; Al Qaderi, M. Simulation of Trajectory Tracking and Motion Coordination for Heterogeneous Multi-Robots System. Jordan J. Mech. Ind. Eng. 2021, 15, 337–345. [Google Scholar]
  6. Aminifar, F.; Rahmatian, F. Unmanned aerial vehicles in modern power systems: Technologies, use cases, outlooks, and challenges. IEEE Electrif. Mag. 2020, 8, 107–116. [Google Scholar] [CrossRef]
  7. Salahat, E.; Asselineau, C.A.; Coventry, J.; Mahony, R. Waypoint planning for autonomous aerial inspection of large-scale solar farms. In Proceedings of the IECON 2019-45th Annual Conference of the IEEE Industrial Electronics Society, Lisbon, Portugal, 14–17 October 2019; Volume 1, pp. 763–769. [Google Scholar]
  8. Ismail, H.; Chikte, R.; Bandyopadhyay, A.; Al Jasmi, N. Autonomous detection of PV panels using a drone. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, Salt Lake City, Utah, USA, 11–14 November 2019; American Society of Mechanical Engineers: New York, NY, USA, 2019; Volume 59414, p. V004T05A051. [Google Scholar]
  9. Felsch, T.; Strauss, G.; Perez, C.; Rego, J.M.; Maurtua, I.; Susperregi, L.; Rodríguez, J.R. Robotized Inspection of Vertical Structures of a Solar Power Plant Using NDT Techniques. Robotics 2015, 4, 103–119. [Google Scholar] [CrossRef] [Green Version]
  10. Lee, D.H.; Park, J.H. Developing inspection methodology of solar energy plants by thermal infrared sensor on board unmanned aerial vehicles. Energies 2019, 12, 2928. [Google Scholar] [CrossRef] [Green Version]
  11. Rezk, M.; Aljasmi, N.; Salim, R.; Ismail, H.; Nikolakakos, I. Autonomous PV Panel Inspection with Geotagging Capabilities Using Drone. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, Virtual, Online, 1–5 November 2021; American Society of Mechanical Engineers: New York, NY, USA, 2021; Volume 85611, p. V07AT07A040. [Google Scholar]
  12. Franko, J.; Du, S.; Kallweit, S.; Duelberg, E.; Engemann, H. Design of a multi-robot system for wind turbine maintenance. Energies 2020, 13, 2552. [Google Scholar] [CrossRef]
  13. Iyer, S.; Velmurugan, T.; Gandomi, A.H.; Noor Mohammed, V.; Saravanan, K.; Nandakumar, S. Structural health monitoring of railway tracks using IoT-based multi-robot system. Neural Comput. Appl. 2021, 33, 5897–5915. [Google Scholar] [CrossRef]
  14. Roy, S.; Vo, T.; Hernandez, S.; Lehrmann, A.; Ali, A.; Kalafatis, S. IoT Security and Computation Management on a Multi-Robot System for Rescue Operations Based on a Cloud Framework. Sensors 2022, 22, 5569. [Google Scholar] [CrossRef]
  15. Kendoul, F. Survey of advances in guidance, navigation, and control of unmanned rotorcraft systems. J. Field Robot. 2012, 29, 315–378. [Google Scholar] [CrossRef]
  16. Verma, J.K.; Ranga, V. Multi-robot coordination analysis, taxonomy, challenges and future scope. J. Intell. Robot. Syst. 2021, 102, 1–36. [Google Scholar] [CrossRef] [PubMed]
  17. Yan, Z.; Jouandeau, N.; Cherif, A.A. A survey and analysis of multi-robot coordination. Int. J. Adv. Robot. Syst. 2013, 10, 399. [Google Scholar] [CrossRef]
  18. Zhilenkov, A.A.; Chernyi, S.G.; Sokolov, S.S.; Nyrkov, A.P. Intelligent autonomous navigation system for UAV in randomly changing environmental conditions. J. Intell. Fuzzy Syst. 2020, 38, 6619–6625. [Google Scholar] [CrossRef]
  19. Huang, L.; Song, J.; Zhang, C.; Cai, G. Design and performance analysis of landmark-based INS/Vision Navigation System for UAV. Optik 2018, 172, 484–493. [Google Scholar] [CrossRef]
  20. Qiu, Z.; Lin, D.; Jin, R.; Lv, J.; Zheng, Z. A Global ArUco-Based Lidar Navigation System for UAV Navigation in GNSS-Denied Environments. Aerospace 2022, 9, 456. [Google Scholar] [CrossRef]
  21. Kim, J.; Kim, S.; Ju, C.; Son, H.I. Unmanned aerial vehicles in agriculture: A review of perspective of platform, control, and applications. IEEE Access 2019, 7, 105100–105115. [Google Scholar] [CrossRef]
  22. Tzafestas, S.G. Mobile robot control and navigation: A global overview. J. Intell. Robot. Syst. 2018, 91, 35–58. [Google Scholar] [CrossRef]
  23. Cabecinhas, D.; Cunha, R.; Silvestre, C. A globally stabilizing path following controller for rotorcraft with wind disturbance rejection. IEEE Trans. Control Syst. Technol. 2014, 23, 708–714. [Google Scholar] [CrossRef]
  24. Mutawe, S.; Hayajneh, M.; BaniHani, S. Robust Path Following Controllers for Quadrotor and Ground Robot. In Proceedings of the 2021 International Conference on Electrical, Communication, and Computer Engineering (ICECCE), Kuala Lumpur, Malaysia, 12–13 June 2021; pp. 1–6. [Google Scholar]
  25. Kim, A.R.; Keshmiri, S.; Blevins, A.; Shukla, D.; Huang, W. Control of multi-agent collaborative fixed-wing UASs in unstructured environment. J. Intell. Robot. Syst. 2020, 97, 205–225. [Google Scholar] [CrossRef]
  26. Elkaim, G.H.; Lie, F.A.P.; Gebre-Egziabher, D. Principles of guidance, navigation, and control of UAVs. In Handbook of Unmanned Aerial Vehicles; Springer: Berlin/Heidelberg, Germany, 2015; pp. 347–380. [Google Scholar]
  27. Edlerman, E.; Linker, R. Autonomous multi-robot system for use in vineyards and orchards. In Proceedings of the 2019 27th Mediterranean Conference on Control and Automation (MED), Akko, Israel, 1–4 July 2019; pp. 274–279. [Google Scholar]
  28. Singh, Y.; Bibuli, M.; Zereik, E.; Sharma, S.; Khan, A.; Sutton, R. A novel double layered hybrid multi-robot framework for guidance and navigation of unmanned surface vehicles in a practical maritime environment. J. Mar. Sci. Eng. 2020, 8, 624. [Google Scholar] [CrossRef]
  29. Al-Fetyani, M.; Hayajneh, M.; Alsharkawi, A. Design of an executable anfis-based control system to improve the attitude and altitude performances of a quadcopter drone. Int. J. Autom. Comput. 2021, 18, 124–140. [Google Scholar] [CrossRef]
  30. Kurak, S.; Hodzic, M. Control and estimation of a quadcopter dynamical model. Period. Eng. Nat. Sci. 2018, 6, 63–75. [Google Scholar] [CrossRef]
  31. Hayajneh, M. Experimental validation of integrated and robust control system for mobile robots. Int. J. Dyn. Control 2021, 9, 1491–1504. [Google Scholar] [CrossRef]
  32. Mutawe, S.; Hayajneh, M.; Al Momani, F. Accurate State Estimations and Velocity Drifting Compensations Using Complementary Filters for a Quadrotor in GPS-Drop Regions. Int. J. Eng. Appl. 2021, 9, 317–326. [Google Scholar] [CrossRef]
Figure 1. Multi-robot coordinates and control structure. Subscript “d” means “desired”.
Figure 1. Multi-robot coordinates and control structure. Subscript “d” means “desired”.
Drones 06 00332 g001
Figure 2. Attractive force in 2D space.
Figure 2. Attractive force in 2D space.
Drones 06 00332 g002
Figure 3. Repulsive force in 2D space.
Figure 3. Repulsive force in 2D space.
Drones 06 00332 g003
Figure 4. Multi-robot systems: (a) indoor and (b) outdoor.
Figure 4. Multi-robot systems: (a) indoor and (b) outdoor.
Drones 06 00332 g004
Figure 5. Surveying mission at the campus of the Hashemite University: (a) a photo from the drone’s camera and (b) planned automatic mission on a map using sequence of way-points.
Figure 5. Surveying mission at the campus of the Hashemite University: (a) a photo from the drone’s camera and (b) planned automatic mission on a map using sequence of way-points.
Drones 06 00332 g005
Figure 6. A comparison between PID and BS controllers for quadrotors and ground robots (speed = 3 m/s).
Figure 6. A comparison between PID and BS controllers for quadrotors and ground robots (speed = 3 m/s).
Drones 06 00332 g006
Figure 7. A comparison between PID and BS controllers for quadrotors.
Figure 7. A comparison between PID and BS controllers for quadrotors.
Drones 06 00332 g007
Figure 8. A comparison between PID and BS controllers for ground robots.
Figure 8. A comparison between PID and BS controllers for ground robots.
Drones 06 00332 g008
Figure 9. A comparison between PID and BS controllers due to disturbance.
Figure 9. A comparison between PID and BS controllers due to disturbance.
Drones 06 00332 g009
Figure 10. Multi-robot mission of two drones and one ground robot.
Figure 10. Multi-robot mission of two drones and one ground robot.
Drones 06 00332 g010
Figure 11. Positions and velocities of two drones and one ground robot mission.
Figure 11. Positions and velocities of two drones and one ground robot mission.
Drones 06 00332 g011
Figure 12. Three-drone mission with one leading the others on a predefined path: (a) four way-points on a map (b) the trajectories of the three drones.
Figure 12. Three-drone mission with one leading the others on a predefined path: (a) four way-points on a map (b) the trajectories of the three drones.
Drones 06 00332 g012
Figure 13. Positions and velocities for the three drone mission.
Figure 13. Positions and velocities for the three drone mission.
Drones 06 00332 g013
Figure 14. Three drones surveying a large area: (a) surveyed area on a map, and (b) the path of each drone with the position of each picture (P).
Figure 14. Three drones surveying a large area: (a) surveyed area on a map, and (b) the path of each drone with the position of each picture (P).
Drones 06 00332 g014aDrones 06 00332 g014b
Figure 15. Positions and velocities for the three drone mission.
Figure 15. Positions and velocities for the three drone mission.
Drones 06 00332 g015
Figure 16. Two ground robots in monitoring task: (a) robot’s path on a map using a sequence of way-points, and (b) the performed path of each robot.
Figure 16. Two ground robots in monitoring task: (a) robot’s path on a map using a sequence of way-points, and (b) the performed path of each robot.
Drones 06 00332 g016
Figure 17. Positions and velocities for the two ground robots mission.
Figure 17. Positions and velocities for the two ground robots mission.
Drones 06 00332 g017
Table 1. Parameters used in backstepping controller of a quadrotor. Subscript “d” means “desired”.
Table 1. Parameters used in backstepping controller of a quadrotor. Subscript “d” means “desired”.
ParameterDescriptionParameterDescription
k 1 , k 2 , k 3 , k 4 , l 1 , l 2 , k γ , P m a x controller gains σ ( x ) P m a x x 1 + | | x | | Saturation function
p [ x , y , z ] T p d Desired path position as a function of γ
v [ u , v , w ] T p d , p d First and second partial derivatives of p d
ω [ P , q , r ] T u 1 , u 2 , u 3 [ 1 , 0 , 0 ] T , [ 0 , 1 , 0 ] T , [ 0 , 0 , 1 ] T
R Rotation matrix of Euler angles ( ϕ , θ , ψ ) S ( x ) Skew symmetric matrix
Table 2. Mean squared errors (MSEs) for PID-controlled and BS-controlled quadrotors at different velocities.
Table 2. Mean squared errors (MSEs) for PID-controlled and BS-controlled quadrotors at different velocities.
Parameterv = 1 m/sv = 3 m/sv = 4 m/s
PIDBSPIDBSPIDBS
M S E x y 0.10740.09960.37120.10210.49000.1160
M S E v x y 0.18890.11890.45030.08391.17520.4881
M S E z 0.07070.02470.04750.02950.07540.0444
M S E v z 0.05440.04150.08650.08240.09440.0997
Table 3. Mean squared errors (MSEs) for PID-controlled and BS-controlled ground robots at different velocities.
Table 3. Mean squared errors (MSEs) for PID-controlled and BS-controlled ground robots at different velocities.
Parameterv = 1 m/sv = 3 m/sv = 4 m/s
PIDBSPIDBSPIDBS
M S E x y 0.07640.01980.26140.11060.43120.1995
M S E v x y 0.17960.08450.71090.31081.24500.6754
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hayajneh, M.; Al Mahasneh, A. Guidance, Navigation and Control System for Multi-Robot Network in Monitoring and Inspection Operations. Drones 2022, 6, 332. https://0-doi-org.brum.beds.ac.uk/10.3390/drones6110332

AMA Style

Hayajneh M, Al Mahasneh A. Guidance, Navigation and Control System for Multi-Robot Network in Monitoring and Inspection Operations. Drones. 2022; 6(11):332. https://0-doi-org.brum.beds.ac.uk/10.3390/drones6110332

Chicago/Turabian Style

Hayajneh, Mohammad, and Ahmad Al Mahasneh. 2022. "Guidance, Navigation and Control System for Multi-Robot Network in Monitoring and Inspection Operations" Drones 6, no. 11: 332. https://0-doi-org.brum.beds.ac.uk/10.3390/drones6110332

Article Metrics

Back to TopTop