Abstract

In this work, we solve the uncertain unmanned aerial vehicle smooth landing problem over a moving platform, assuming that the aircraft position relative to the platform and its acceleration is always measurable. The landing task is carried out by an output-feedback robust controller, together with a repulsive force. The robust controller controls the nominal model, accomplishes the needed tracking trajectory, and counteracts the unknown uncertainties. To assure that the aircraft is always above the platform, we include a repulsive force that only works in a small vicinity of the platform. To estimate the relative aircraft velocity and platform acceleration, we use a supertwisting-based observer, assuring finite-time convergence of these signals. This fact allowed us to design the feedback state stabilizer independently of the observer design (in accordance with the separation principle). We confirmed the effectiveness of our control approach by convincing numerical simulations.

1. Introduction

The challenge of landing unmanned aerial vehicles on moving platforms has attracted the attention of several researchers in the control community. It is a problem that, if solved, has many applications, such as landing an unmanned aerial vehicle (UAV) on the deck of a ship in high seas. This control task consists of regulating to zero the relative distance between the aircraft and the ship deck, the latter being subject to external motions, including pitch or vertical movements. This regulation task has to be accomplished smoothly, safely, precisely, and fast, with the primary concern of safeguarding the physical integrity of both the aircraft and the ship. Usually, this control task has been accomplished using identification, adaptive algorithms, sliding-mode control, robust control, and optimal control [14]. In the literature, we can find several control schemes to address the UAV shipboard landing problem. In [2], an optical flow based control law for landing on a moving platform is proposed. In [5], a vision-based algorithm is presented, where the ground effect experienced during the maneuver task is compensated for using by a sliding adaptive control law. In [6], another vision-based landing algorithm for an autonomous helicopter is presented. The landing algorithm is integrated with algorithms for the visual acquisition of the target and navigation to the goal, from an arbitrary initial position and orientation. In [7], a control strategy for autonomously landing a helicopter on a rocking ship (due to rough seas) is presented. This control strategy is implemented by two controllers that separate the time scales of translation and rotation. Marconi et al. [1] present a control strategy for the autonomous landing of a vertical take-off and landing vehicle (VTOL) on a ship with a deck that oscillates vertically due to high seas. There, the authors design an internal model-based error feedback dynamic regulator that is robust with respect to some uncertainties about the mechanical parameters that characterize the model and secures global convergence. Yang et al. [8] present a procedure for predicting a ship’s pitch and heave motions. This procedure is based on a predictive controller whose inputs are measurements of relative distance and whose outputs are the predictions of the pitch movements. Other exciting works based on the predictive control approach to solving this control task are [912]. In [13], the UAV landing problem is addressed using a three-step procedure: motion estimation, trajectory generation, and tracking control, where the latter is accomplished as an optimal control. The optimal control approach for landing UAVs is also used in two recently published works [14, 15]. A full review of the literature related to this problem is beyond the scope of this work; however, we suggest to the interested reader the following works [1625].

In this work, we approach the problems of feedback regulation and output regulation to accomplish the landing control maneuver task of successfully landing a vertical take-off and landing unmanned aerial vehicle (VTOL-UAV) over a vertically moving platform, under the assumptions that the vehicle model is partially known, and its relative position with respect to the platform and the VTOL-UAV vertical acceleration are both available. Due to the nature of this maneuver, we can formulate it as if it were a synchronization of two partially known systems, where the platform is the master, and the VTOL-UAV is the slave. To this end, we use the aircraft simplified dynamic model. That is, we neglected the dynamics of the other coordinates. That is, we assume that the aircraft horizontal and angular displacements are conveniently small, and in an actual implementation, such displacements can be eliminated by suitable robust controllers. Then, to solve this problem, we designed a control strategy were we first proposed a nonlinear control, consisting of two actions, in conjunction with a repulsive force. The first action, based on saturation, controls the nominal model, while the other action, which consists of a nonlinear integrator, compensates for the partially known nonlinearities. The controller and the compensator work together with the well-known supertwisting-based observer (STBO), which is in charge of estimates on a finite-time vehicle’s relative velocity and moving platform’s acceleration. The complementary repulsive force helps to guarantee that the aircraft is always above the platform. We use the Lyapunov method to carry out the closed-loop stability analysis. The novelty of this study is that, as far as we know, the Bézier spline-based trajectory planning in conjunction with a repulsive force has not been used before to solve the aforementioned control problem.

We organized the remainder of the study as follows: in Section 2, we introduce the VTOL-UAV model and establish the control problem; we develop the corresponding control strategy in Section 3, where we implement the needed STBO; in Section 4, we present the outcomes of the numerical simulations, which allow us to claim that our control strategy effectively lands the VTOL-UAV over the moving platform; finally, Section 5 gives the conclusions and final remarks.

Notation: the symbols and sgn refer, respectively, to a linear saturation function and the signum function of a real number, that is,

2. Problem Statement

To formally establish the control problem, we introduce a simplified model of a VTOL-UAV, described by the following system of equations:where is the vertical aircraft position; is the actual vertical rotor’s thrust; is the relative distance between the VTOL-UAV position, , and the position of the vertical moving platform, ; is the VTOL-UAV mass; is the gravity force; and is a nonlinear function modeled as follows (see equation (31) in [2]):where the positive constants and , with , can be identified from the physical. This nonlinear function is the thrust that the VTOL rotors generate. When the VTOL is approaching to the ground, the rotors will generate more thrust for a given power. This phenomenon is known as the ground effect [26]. These parameters are related to the aircraft physical dimensions [27]. The platform position is modeled aswhere is the amplitude of the component of the sinusoidal motion, is the frequency, and is the phase. Loosely speaking, models the platform vertical oscillatory movement due to the ocean waves. Having introduced the simplified VTOL-UAV system model, we present the primary concern of this work.

2.1. Control Problem

Consider the problem (2), with measurable outputs and and constants and unknown. The control objective consists in designing a control such thatfor all and where and are sufficiently small. In other words, we desire to render the system variables to a neighborhood of zero. That is to say, we want states of the VTOL-UAV to converge with that of the platform , with sufficient accuracy for some , restricted to for all .

To solve the control problem mentioned above, we introduce the following assumptions:A1: , with and A2: the signals and are always measurable

Notice that we must assure that for all because we assume that the VTOL-UAV is always above the platform. This control problem resembles, for instance, the landing of a helicopter on a ship whose deck oscillates vertically due to sea waves, as depicted in Figure 1. We must note that the reference signal and its first and second time derivatives have to be bounded (see , with ). On the contrary, we deal with the ideal case when signals and are noise free. However, if these signals are not noise free, we can overcome this problem filtering them using, for instance, a low-pass filter, obtaining the signals average values, which are very close to their actual values.

We finish this section by introducing a useful lemma, used in the forthcoming developments.

Lemma 1. Consider the following uncertain second-order system:where , are the states, is a continuous bounded perturbation, and is the system input, defined aswhere is a constant and is the smooth reference signal with time derivatives and . Suppose that with , , and , as long as , withthen the tracking errors, and , globally and asymptotically converge to the origin. Notice that with are strictly positive and constant bounds.

The proof of this lemma can be found in Appendix.

Remark 1. We must emphasize that the signal (refST0) is a bounded control. Consequently, it can be substituted by any other bounded controller, like a nested saturation-based controller, or by a twisting controller. For instance, the twisting controllercan be a suitable candidate. However, we use (7) because its stability analysis is easy to carry out.

3. Control Strategy

Roughly speaking, we can see the control problem established above as an uncertain system synchronization problem, where the platform acts as the master and the VTOL-UAV as the slave. To this end, we propose an output-feedback robust controller, together with a repulsive field-based force. The robust controller is devoted to assuring tracking error convergence to zero of the nominal model, and it simultaneously compensates for the unmodeled nonlinearities. The repulsive force acts if , where is related to the constant defined buffer zone. As we only have measurements of the relative position and the VTOL-UAV vertical acceleration , we estimate the unknown variables and by using an STBO [28], assuming that is also partially known. Having obtained and , we proceed to design a robust output-feedback-based controller together with a trajectory planner. The planner allows us to render the aircraft to the buffer zone and lands it very smoothly on the platform. We accomplish the reference trajectory using a fifth-grade Bézier spline, allowing us to plan the landing time.

Comment 1. In our opinion, the advantage of our proposal is that we give a robust control-based solution, which considers the dynamic variations in function given in (3), based on saturation functions, together with a nonlinear integrator. This combination has the advantage of solving the problem in a precise manner when the external perturbations have slow variations. Contrary, for instance, to the adaptive and PI control approaches where we use an approximation to solve this kind of problem, we assume that the time derivative of is equal to zero. On the contrary, if we compare our strategy with other well-known strategies, it has the drawback of having a slower time response, and the control gain parameters are difficult to tune adequately.

3.1. Reconstruction of the Unknown Variables

Defining the synchronization errors as and , we have thatwhere is in fact a measurable signal. Then, for estimating the missing variables and , we use the following observer [28, 29]:where constants and must be selected, such thatfor some and . Defining the observation errors as and , we obtain, from the first two equations of (10) and (11), the following equations:

As constants and are selected according to (12), we can assure that and asymptotically converge to zero in a finite time (for details see Sections III and IV of [28]). On the contrary, if we fix the step size of the numeric integration method, , such that , then we can assure that in a finite time, as long as [29]. That is, there exists , such that and , for all.

The disadvantage of this sliding-mode observer is that we need to know, at least, an estimated value for the constant . However, we can overcome this inconvenience by selecting the constant large enough. Similarly, to reconstruct , we need to select and small enough, depending on the natural frequency of .

3.2. Robust Control Law

Since we can recover variables and in a finite time, we propose a suitably robust control scheme for assuring that and asymptotically converge to zero. Before proceeding, we introduce the following useful proposition.

Proposition 1. Let us consider the following uncertain second-order system:where and are the states, , is the system input, and and are partially unknown functions. Let as assume the following:P1: there exits and , such that , for P2: for all Consider the following controller:where , , and is the continuous and bounded reference trajectory, with the first and the second time derivatives being also continuous and bounded, and the auxiliary control is defined aswithwhere the parameters are strictly positive with . Then, the closed-loop systems (14) and (15) assure that the tracking errors, and , asymptotically converge to zero.

Proof. From (14) and (15), we can write the closed-loop system asNow, from the definition of (18), it is easy to see that we can express the equation above in a short way asBefore continuing with the stability analysis of the above system, we first show that converges to zero. To this end, we must note that the variable , given in (18), can also be expressed asIndeed, it can be justified using the simple integration of the second equation of (19) and the definition (18). From assumption P3, we must emphasize that can be computed online with equation (18). Having defined the new system representation (20), we proceed to show that the virtual controller (17) makes variables and to asymptotically converge to zero. To this end, we substitute (17) into the time derivative of , given in equation (15), leading towhereHence, to analyze the behavior of the auxiliary variable , we use the Lyapunov function , in which the time derivative around the trajectories of (22) can be upper bounded aswhere and . Hence, selecting and using P1, we obtain that the following inequalityholds. From the last differential inequality, it is easy to conclude that exponentially converges to zero (see Appendix). Hence, the variable also exponentially converges to zero. Now, as (16) and are bounded, then is bounded (see (22)), implying that is uniformly continuous with , as long as . Thus, according to the lemma of Barbalat, we can conclude that , as long as .
Finally, after substituting (16) into equation (20), we obtainwhere is bounded and converges to zero. Besides, as exponentially converges to zero, we have that . Consequently, according to Lemma 1, we have that , as long as .

Remark 2. If we substitute the auxiliary controller , defined in (17), bythe variable will converge in a finite time. On the contrary, the performance of system (14), in closed loop with (15) and (16), could be improved if we substituted controller with any other control based on the saturation function, as follows:orwhere and are both positive constants, conveniently tuned.

3.2.1. The Dynamical Synchronization Equation

Substituting the equation of the VTOL-UAV system (2) into equation (10), we obtain the following synchronization equation:where and . Now, to propose the control to assure , as long as , we assume that we know , , and such that

To accomplish this, we introduce the main proposition of this work.

Proposition 2. Consider system (30) in closed loop withwhere is a Bézier spline and and evolve according to the system equations in (11), restricted to (12), and is proposed aswith, and

If constants and satisfy , then for any initial condition , we can assure that errors and asymptotically converge to zero.

Proof. See Appendix.

3.2.2. Planification of the Reference Trajectory

To land the VTOL-UAV on the moving platform, we use a Bézier spline-based trajectory [30]. We most note that the Bézier spline allows us to design translation trajectories from one equilibrium point to another equilibrium point, in a previously fixed time. We propose this trajectory to bring the aircraft, from any initial condition , to the buffer zone with being sufficiently small and conveniently fixed. That is, we define aswhere is the selected Bézier spline, which satisfies the following conditions:where a finite number of the time derivatives of are set to zero at the initial and final maneuver times and (for details, please refer [30]). To assure that both for all times and that the controller turns off when the VTOL-UAV reaches the buffer zone, we fixed it, for practical purpose, as

It is worthy to mention that acts as decision function, which allows to decide when to deactivate the control action. That is, when the aircraft is above and close enough to the platform, the controller is set to zero; otherwise, it is on. Thus, we establish the practical control strategy as follows:where and are defined in Proposition 2 and is a function that acts as a switch and is defined asand is the repulsive force defined aswith its corresponding potential defined aswhere . We must note that the repulsive force assures that . We omitted the stability proof of system (2) in closed loop with control (39) because it mimics Proposition 1.

Comment 2. Please note that the slower the aircraft descends, the lesser the possibility to crash it with the platform. We can expect the same result if we increase either the size of the buffer zone or the magnitude of the repulsive force, by conveniently tuning, for instance, and being sufficiently large. However, the crashing possibility cannot be avoided by any methodology, if an unexpected and strong enough crosswind appears suddenly.

4. Numerical Simulations

To assess the effectiveness of the proposed control strategy, we design and run three numerical experiments. In the first experiment, we perturb the aircraft with a crosswind. The second experiment consists of comparison among our control strategy (OCS), and the continuous and smooth PID-based strategy (PIDS), and an adaptive strategy (AS); the latter two are, respectively, proposed in [4, 13]. Finally, the third experiment is also a performance comparison, among OCS and a version of the discontinuous twisting-based strategy (TBS), found in [31].

4.1. Experiment Setup

We fix the system physical parameters as , , and and the controller gains, given in (15) to (17), asand the STBO gains as

For the experiments two and three, we suppose that the system acceleration is noisily perturbed by a normally distributed random process with zero mean and whose standard deviation is 1, defined as where is the random noise.

4.2. First Experiment

For this experiment, we fix the initial condition as and the Bézier spline parameters as and simulate the platform hypothetical vertical movement, , and the crosswind effect, , as follows:

We show the outcome of this experiment in Figure 2, where we see that the VTOL-UAV reaches the buffer zone after 27 seconds and the estimated velocity matches the actual velocity after 5 seconds. Consequently, we can conclude that our control approach effectively lands the aircraft over the platform, with a good performance. In fact, the final error is in the order of meters. We must note that there is a small error between and . This error can be explained because control is based on a saturation function, which has the inconvenience of having a significant time of convergence. However, if time is sufficiently large, then we can improve the tracking error.

4.3. Second Experiment

For this experiment, we fix the initial condition as and the Bézier spline parameters as . We simulate the platform hypothetical vertical movement as in the first experiment, while we externally perturbed the system by a crosswind, defined as , and we fixed the OCS saturation function parameter as . As we already mentioned, the experiment goal is to make a comparison among the OCS, PIDS, and AS strategies and to have an idea of how good our solution is with respect to the other two well-established approaches. In order to accomplish a fair comparison, we assume that is measurable, and we add the repulsive force to the strategies PIDS and AS and use the same control gains for both strategies. Then, we propose the PIDS asand the AS aswhere , , and and is the gain parameter estimation of based on the projective adaptive control law proposed in [13] and the control gains fixed as , , and , with proposed as in (41). We show the obtained results from the numerical simulation in Figure 3. In this figure, we show the performance index of the three control strategies, which is defined asand their corresponding control actions. From the performance indexes, we can see that OCS has a better performance than PIDS, while AS exhibits the less efficient performance index. In the same figure, we can see that the OCS control action shows fewer overshoots; however, the chattering phenomenon is always present. Between the PIDS and AS control actions, we can say that the first one exhibits more overshoots. We can also see that the PIDS strategy exhibits chattering when it is close to the equilibrium point, that is, after 20 seconds. Finally, in Figure 4 we show the trajectory tracking position error for the three control strategies. From this figure, we can see that OCS has a slightly better performance. Please note that the overshoots that we can see after 27 seconds are due to the effect of the repulsive force.

4.4. Third Experiment

The control task is the same as in the second experiment, but we fix the Bézier parameters as . In this experiment, once again, we perturbed the acceleration equation by adding it to the noisy signal . Additionally, to make the experiment more realistic and illustrative, we added a noise signal to the output . That is, , where and are, respectively, the measurable and actual outputs, and is a normally distributed random process with zero mean; we assume that is not available and estimated using the STBO (11). Then, we use the following TBS:where the controller gains satisfy and the following conditions:whereand is any constant. We most note that, for this experiment, we fix . Hence, according to (51), we can select and . Then, to make the comparison, we use the following error index:

We show the comparison outcome in Figure 5, where we can see that, as we expected, the TBS outperforms OCS because the closed-loop response of the TBS is almost immune to external matching perturbation, having the inconvenience of exhibiting the chattering phenomena in the controller, as we can see it in Figure 5, preventing its actual implementation. In favor of OCS, we mention that the controller has a smooth response, which, in our opinion, is due to the integration action because it is averaging the noise effect over the system. This behavior can allow us to implement OCS for real applications.

Comment 3. We must mention that we implement the sign function aswith integration step fixed as .

5. Conclusions

In this work, we propose a control strategy for landing a VTOL-UAV smoothly on a vertically moving platform. We designed this control strategy assuming that the VTOL-UAV model was partially known, while its position relative to the platform and the platform acceleration were both available. Our landing control approach is based on designing an output-feedback robust controller, in conjunction with a repulsive force. The robust controller consists of two control actions: one based on saturation functions, and the other is a nonlinear integrator. The former controls the nominal model and assures the tracking trajectory task, and the latter compensates for the partially known uncertainties. Additionally, we introduce a repulsive force controller to guarantee that the aircraft is always above the platform, and it only works within a very small vicinity (buffer zone). We use an STBO to estimate the nonavailable aircraft relative velocity and the moving platform acceleration. We select the aircraft reference trajectory as a Bézier spline, devoted to bringing the system from any initial condition towards the buffer zone in a fixed time. To assess the effectiveness of our control method, we carry out two numerical simulations, the first of which performs the landing task when a crosswind perturbs the aircraft. The second experiment consists of a comparison among the strategies OCS and AS. From the results obtained in this simulation, shown in Figure 3, we can see that our control strategy has a better performance than the other two. Finally, we carry out the corresponding stability analysis based on the Lyapunov method.

Appendix

Proof. Proof of Inequality (25): first of all, multiplying both sides of (25) by , we obtainwhere . For simplicity, we suppose that . Now, integrating both sides of the above inequality, we obtainwhich leads us, after using some simple algebra, to the following inequality:From the above, we conclude that exponentially converges to zero. We can probe the case when is in a similar fashion.

Proof of Lemma 1. : after some manipulation, it is easy to see that system (6), in closed loop with (7), can be written aswhere and . Notice that the right-hand side of the equation above is Lipchitz, then states and remain bounded in a finite time. That is, the system is forward complete. Now, introducing the linear transformations and into system (A.4), we obtainNow, since , as long as , we can select , such that , for all . Thus, selecting the positive definite function , we obtain that its time derivative, around the trajectories of the second equation of (A.5), leads to the following inequality:Because system (A.5) is Lipchitz with respect to the states , they cannot exhibit finite time of escape. Notice that, if , then from (A.6), we have that . This implies that there exists , such that , for all . When this condition is satisfied, we evidently have that . Consequently, the first equation of (A.5) reads asfor all . Using the positive definite function and recalling that , for all , we have that the time derivative of , around the trajectories of the system above, leads to the following inequality:Observe that, if , then . Thus, there exists a finite time such that for all . After these conditions hold, system (A.5) reads aswhere andAccording to (A.9), we have thatUsing the integration by parts method into the equation above, we have thatBecause the matrix is Hurwitz, there exist and , such thatHence, the solution of (A.12) can be upper bounded aswhere , , and . Notice thatNow, using the Cauchy mean value theorem, we can upper bound the first term of the right-hand side of (A.15) as follows (the fact that and are bonded implies that , as long as, ):for some , leading us to conclude thatSimilarly, the second term of the right-hand side of (A.15) satisfies the following equality:for some , implying thatTherefore, from inequalities (A.14)–(A.19), we conclude that the state converges to zero for any initial condition. That is, system (A.4) converges globally and asymptotically to the origin.

Proof of Proposition 2. For simplicity, we carry out the proof assuming that . As and are estimated by the STBO (11), with its corresponding gains fulfilling condition (12), we can conclude that converges in finite time to . That is,with and bounded, such that and , as long as , with being finite. Under these conditions, it is easy to see that the trajectories of the closed-loop system (32) with (30) does not escape to infinity for all . Moreover, after , signals , , , and are given as follows:where is defined in (33). As the above signals coincide with the expression given in Proposition 1, with , then we can assure that and converge to zero, as long as

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This research was supported by Secretaria de Investigacion y Posgrado of the Instituto Politecnico Nacional, under research grant nos. 20200671, 20201623, and 20201829.