Next Article in Journal
Systematic Design of a 3-DOF Dual-Segment Continuum Robot for In Situ Maintenance in Nuclear Power Plants
Next Article in Special Issue
Navigation of a Differential Wheeled Robot Based on a Type-2 Fuzzy Inference Tree
Previous Article in Journal
Preparation Method of Lunar Soil Simulant and Experimental Verification of the Performance of an Impact Penetrator for Lunar Soil Exploration
Previous Article in Special Issue
Virtual Neuromuscular Control for Robotic Ankle Exoskeleton Standing Balance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Flatness-Based Active Disturbance Rejection Control for a PVTOL Aircraft System with an Inverted Pendular Load

by
Cesar Alejandro Villaseñor Rios
1,†,
Alberto Luviano-Juárez
1,†,
Norma Beatriz Lozada-Castillo
1,†,
Blanca Esther Carvajal-Gámez
1,†,
Dante Mújica-Vargas
2,† and
Octavio Gutiérrez-Frías
1,*
1
Instituto Politécnico Nacional, Unidad Profesional Interdisciplinaria en Ingeniería y Tecnologías Avanzadas, Mexico City 07340, Mexico
2
TecNM-CENIDET, Cuernavaca 62490, Morelos, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Machines 2022, 10(7), 595; https://doi.org/10.3390/machines10070595
Submission received: 22 June 2022 / Revised: 15 July 2022 / Accepted: 19 July 2022 / Published: 21 July 2022
(This article belongs to the Special Issue Control of Robotic Systems)

Abstract

:
This paper presents a systematic procedure for the control scheme design for a PVTOL aircraft system with an inverted pendular load, which is a nonlinear underactuated system. The control scheme is based on the use of angular movement as an artificial control in order to propose new auxiliary control inputs. This is achieved by a linear extended state observer-based active disturbance rejection control to reject both nonmodeled dynamics and external disturbances. The flying planar inverted pendulum is then linearized around an unstable equilibrium point, and the resulting system is subdivided into two subsystems: (1) the height system, and (2) the horizontal pendulum system. For the height system, a linear extended state observer-based active disturbance rejection control is proposed in order to accomplish a take-off and landing task in the presence of external disturbances and non-linearities neglected in the linearization process. The flatness property in the horizontal-pendulum system is exploited in order to propose another active disturbance rejection control of linear nature. The flatness of the tangentially linearized model provides a unique structural property that results in an advantageous low-order cascade decomposition of the linear extended state observer design. Numerical simulations show the effectiveness of the proposed control scheme in trajectory tracking tasks in the presence of disturbances caused by crosswinds with random amplitudes.

1. Introduction

The analysis and control of underactuated systems (UAS) is an active and intricate research topic in automatic control theory [1] because control inputs cannot produce arbitrary accelerations in some of the degrees of freedom at every instant in time [2]. This fact has several untoward consequences for design control laws. For instance, some UAS are not linearizable by static or dynamic feedback, making most of the traditional control schemes inapplicable. In contrast, the need for fewer actuators in UAS produces some practical advantages like lighter structures, lower costs, and reduced dimensions [3]; these have permitted the practical use of UAS in real applications such as damping systems [4], aerial vehicles [5], ground vehicles [6], and flexible robots [7], among others.
A classical benchmark in automatic control theory is the inverted pendulum system (IPS), which includes the pendubot, acrobot, wheeled inverted pendulum, spherical inverted pendulum, Furuta pendulum and, in recent years, the flying inverted pendulum. All these systems are nonlinear, underactuated, and unstable mechanical systems and many control strategies have been developed for stabilization and tracking the use of linear and non-linear control. Furuta et al. [8] are concerned with attitude control of a triple inverted pendulum system in which the lower hinge is free for rotation and two upper hinges are controlled. The controller is designed by using a linearized model in a neighbourhood of the upright position; simulations and experimental results are shown. In [9], robustness improvement in a class of UAS is addressed by combining the so-called IDA-PBC controller with an adaptive control technique. This control scheme is used in the inertia wheel inverted pendulum with numerical simulations and real-time experiments. The proposal in [10] leads global stabilization for the cart-pendulum system by control strategy based on the use of saturation functions. Aguilar–Ibañez et al. [11] present a non-linear controller for the stabilization of the Furuta pendulum by using a partial feedback linearisation. In the first stage, only the actuated coordinate is linearised; then, the stabilizing feedback controller is achieved by applying the Lyapunov direct method. This paper proves local asymptotic stability and demonstrates that the closed-loop system has a large region of attraction.
In recent decades, some popular UAS include the unmanned aerial vehicle (UAV) because of its robust electro-mechanical design, high maneuverability, and low cost; moreover, their use can be appliled to various fields, such as environmental monitoring, terrain mapping, emergency response, and military use, among others [12]. This wide diversity of applications has motivated novel device configurations, instrumentation schemes, and research in control theory to apply linear and nonlinear controls. Classical linear control is based on an approximate linear model near the equilibrium point, and many examples of linear control can be found in the literature, including the proportional integral derivative (PID) control [13] and linear quadratic regulator (LQR) control [14,15] for stabilization and tracking problems. In addition, some interesting works can be found related to nonlinear control approaches for UAS stabilization and tracking. In Hernández et al. [16], a control strategy is developed for take-off and landing maneuvers in a Quadcopter system. The control strategy consists of a combination of controllers based on nested saturations and a generalized proportional integral (GPI) controller, and it is designed considering the presence of disturbances. Consequently, larger robustness is obtained and the algorithm convergence is proven by means of Lyapunov’s second method. Azinheira et al. [17] present a backstepping-based controller used on the UAV system with input saturation, and the stability of the control solution is verified. Another backstepping approach is showcased by Das et al. [18], in which the backstepping approach is used for the quadrotor controller on the Lagrangian form of the dynamics; in addition, two neural networks are introduced to estimate the aerodynamic components. The proposed controller shows robustness cope of unmodeled disturbances. Xiong et al. [19] apply a second-order sliding mode control approach to design a robust flight controller for the Quadcopter system, which is designed to provide a robust position and attitude tracking control with respect to model parameter uncertainties.
Planar vertical take-off and landing (PVTOL) has been used as a benchmark for many research works given that it represents a simplification of other UAV systems, such as Quadcopter and aircraft. PVTOL is a hardly nonlinear control problem because it is an underactuated, nonminimum phase system and, due to the moment-to-force coupling in the model, its zero dynamic is unstable, and the exact input–output linearization methodology produces undesirable results. Numerous works about stabilization control for the PVTOL system can be located in the literature. Hauser et al. [20] showed a model extension in the PVTOL system and proposed a state feedback linearization control for stabilization. In [21], the state feedback control technique in a PVTOL extended model is analyzed, and the conditions under which it is possible to determine local asymptotic stability by using a control based on feedback linearization is found. In [22], Fantoni et al. proposed a global stabilizing scheme for the PVTOL system considering PVTOL’s angle as a fictitious control input. The control scheme is formed by a nonlinear combinations of linear saturation functions bounding the thrust input and the rolling moment. Global convergence of the state is proved. An extension of this work is shown in [23] wherein robustness is proven with respect to lateral force coupling, which was neglected in the previous work.
The robust control design for UAV systems is an active and challenging problem for stabilization and tracking given that these kinds of systems are subjected to disturbance in outdoor applications due to wind; consequently, there is a vast amount of literature on the subject. In [24], a backstepping procedure is applied to control the height position together with the nested saturation functions controller for horizontal position and roll angle stabilization, and an extended state observer is used to estimate the wind disturbance. In [25], a robust control strategy to stabilize a PVTOL aircraft in the presence of crosswind is proposed by the use of robust control Lyapunov functions (RCLF) and Sontag’s universal stabilizing feedback.
In recent years, the flying inverted pendulum (FIP) system has attracted interest in automatic control research; this system consists of an inverted pendulum linked to a UAV system. The FIP system was developed by Hehn and D’Andrea [26] in 2011. The LQR control was used to stabilize the position of the pendulum and Quadcopter systems considering that the pendulum’s mass was less than 5% of the UAV’s mass. Some other works have been developed based on the FIP system by using different control ideas; a three-level cascade strategy is proposed in [27], wherein each level of the cascade system executes a control law designed through a backstepping approach and numerical simulations were carried out. Chen et al. [28] used an improved genetic algorithm to determine parameters in the LQR controller in order to get better stabilization in the FIP system. Numerical simulations were carried out to observe the performance of the control methodology proposed. In [29], a backstepping control law based on geometric principles to swing up the FIP system is proposed, and numerical experiments have shown control actions for aggressive maneuvers.
The active disturbance rejection control (ADRC) has led to a new paradigmatic control vision whereby external and internal disturbances are actively estimated and rejected from a simplified version of the system. The ADRC field exhibits an increasing number of experimental results in diverse applications. In [30], a robust trajectory tracking control for the PVTOL system under crosswind is discussed. The control proposed is a combination of input–-output feedback linearization and ADRC to estimate and compensate the crosswind effects. In [31], a three-loop cascade control strategy based on ADRC is proposed for the FIP system. Simulation results showed a comparison to an LQR controller to illustrate the robustness of the control proposal. A flatness-based approach is proposed for the linear ADRC stabilization of a nonlinear inertia wheel pendulum in [32], which is subject to unmodelled dynamics and external disturbances. The approach exploits the cascade structure of the tangential linearization of the underactuated system and design a high-gain linear cascaded extended state observer of the generalized proportional integral type. Experimental results are shown to demonstrate the effectiveness and feasibility of the proposed approach with respect to a classic control technique.
The aforementioned works use restriction in the pendulum’s mass considered in [26]. Consequently, the pendulum’s dynamics are neglected in the UAV’s dynamics and, in fact, the flying pendulum system is in the presence of external disturbances produced by wind, as reported in above. Considering works in the literature, this paper is focused on the trajectory-tracking problem and stabilization in the PVTOL aircraft system with an inverted pendular load under crosswind disturbances without restriction in the pendulum’s mass. This consideration causes coupling in the pendulum and UAV dynamics. This work explores the differential flatness property in the PVTOL aircraft system with an inverted pendular load in order to design a control scheme that consists of three closed-loop controls; in each one, a linear state feedback control law is proposed, derived from differential flatness system, and an extended state observer (ESO) to estimate and compensate the crosswind is applied. In addition, this proposal considers that only position is available in state measurement. Consequently, the control scheme proposed takes advantage of the ESO to estimate the states that are not directly available for measurements.
The main contributions of this study are as follows:
(i)
a control strategy that uses the differential flatness property in the PVTOL aircraft system with an inverted pendular load in order to design a control for the height, roll attitude, horizontal position, and roll angle simultaneously, even in the presence of a crosswind;
(ii)
a control algorithm, which is robust against the presence of external disturbances and whose performance is competitive with respect to another robust controller of discontinuous nature; and
(iii)
a set of convenient transformations (cascade structure [3]), in which the linearized high-order system can be expressed as a various tandem lower-order systems depending on measurable variables that allow a controller formed by the combination of linear extended state observer-based ADRCs.
The rest of this article is divided as follows. In Section 2 the PVTOL aircraft system with an inverted pendular load is described, and the mathematical model is derived by using Euler–Lagrange formalism. Section 3 is devoted to describing the proposed control approach in four subsections. Section 3.1 presents the control loop for roll angle and a normalized form of the system. Section 3.2 shows the height control, and Section 3.3 presents the control loop for horizontal displacement and pendulum’s angle. Section 4 is devoted to numerical simulations, and Section 5 presents the conclusions.

2. Mathematical Model

The PVTOL aircraft system with an inverted pendular load is a mechanical system that consists in a PVTOL vehicle and an inverted pendulum attached to the PVTOL’s center of gravity; this system is a two-dimensional version of the FIP system. PVTOL aircraft system with an inverted pendular load is shown in Figure 1, and its dynamical model is derived by Euler–Lagrange equations. The Lagrangian of the system is constructed by the difference of the kinetic and potential energies of the PVTOL and the inverted pendulum,
K v = m v ϵ ˙ v ϵ ˙ v 2 + i v α ˙ 2 2 K p = m p ϵ ˙ p ϵ ˙ p 2 + i p θ ˙ 2 2
P v = m v g y v P p = m p g ( l p cos ( θ ) + y v )
where m v is the PVTOL’s mass, m p is the pendulum’s mass, i v is the PVTOL’s inertia, i p is the pendulum’s inertia and l p is the pendulum’s length, ϵ v = ( x v , y v ) is the reference frame fixed to PVTOL’s center of gravity related to an inertial reference frame fixed in the ground. ϵ p is the reference frame fixed to pendulum’s center of gravity and it is related to ϵ v as
ϵ p = ϵ v + R θ 0 l p
R θ is the rotational matrix defined as
R θ = cos ( θ ) sin ( θ ) sin ( θ ) cos ( θ )
The Lagrangian of the PVTOL aircraft system with an inverted pendular load is defined as follows:
L = 2 g l p m p cos ( θ ) 2 g ( m v + m p ) y v + ( m v + m p ) ( x ˙ v 2 + y ˙ v 2 ) 2 + i v α ˙ 2 2 l p m p ( cos ( θ ) x ˙ v + sin ( θ ) y ˙ v ) θ ˙ + ( i p + l p 2 m p ) θ ˙ 2 2 .
By using Euler–Lagrange formalism, the dynamics can be calculated as
d d t L q ˙ L q = U + d ,
where q is the generalized coordinates vector of the system, q = x v y v α θ , U and d are the control inputs matrix and the external disturbance matrix defined as follows:
U = u 1 sin ( α ) u 1 cos ( α ) u 2 0 d = d x sin ( α ) d y cos ( α ) d α 0 .
Control inputs u 1 and u 2 are defined through the forces generated by the pairs propeller-motor of PVTOL, u 1 = f 1 + f 2 and u 2 = f 2 f 1 and disturbances are taken as crosswind inputs as shown in [30] with amplitudes d x , d y , and d α respectively.
By solving the Euler–Lagrange Equation (2), differential equations of the system can be expressed as M ( q ) q ¨ + C ( q , q ˙ ) q ˙ + G ( q ) = U + d , where matrix M ( q ) is symmetric and non-singular, C ( q , q ˙ ) is the Coriolis matrix, and G ( q ) is the gravity vector. Solving this expression for q ¨ , the dynamical system can be expressed as
q ¨ = M 1 ( q ) ( U C ( q , q ˙ ) q ˙ G ( q ) ) + M 1 ( d ) .
Without loss of generality, PVTOL’s structure is taken as a rigid body, the pendulum is assumed as a rigid massless bar with a point mass in its end, i p = 0 , and PVTOL’s inertia is i v = 1 . Solving the expression (3) for the Lagrangian described in (1) of the PVTOL aircraft system with an inverted pendular load, the following differential equations are obtained:
x ¨ v = 0.5 u 1 ( ( 2 m v + m p ) sin ( α ) + m p sin ( α ) + 2 θ ) + l p m v m p sin ( θ ) θ ˙ 2 m ( m v + m p ) + D x
y ¨ v = 2 g m v ( m v + m p ) + ( 2 m v + m p ) u 1 cos ( α ) m p u 1 cos ( α + 2 θ ) + 2 l p m v m p cos ( θ ) θ ˙ 2 2 m v ( m v + m p )
+ D y
α ¨ = u 2 + D α
θ ¨ = u 1 sin ( α + θ ) l p m v + D θ
with
D x = 0.5 d x ( 2 m v + m p + m p cos ( 2 θ ) ) sin ( θ ) + d y m p cos ( α ) cos ( θ ) sin ( θ ) m v ( m v + m p )
D y = d y cos ( α ) ( 2 m v + m p m p cos ( 2 θ ) ) + d x m p sin ( α ) sin ( 2 θ ) 2 m v ( m v + m p )
D α = d α
D θ = d x cos ( θ ) sin ( α ) + d y cos ( α ) sin ( θ ) l p m v .
Control problem. Consider the PVTOL aircraft system with an inverted pendular load dynamical system ((4)– (8)). The problem is to design a control scheme for a take-off and landing maneuver while ensuring the stabilization in the inverted pendulum position around the unstable equilibrium point, through the control inputs u 1 and u 2 in spite of the presence of disturbances caused by crosswind and nonmodeled dynamics.

3. Control Scheme

The control scheme is conformed by three control loops: The first loop is designed for the trajectory tracking control in α angle. This inner loop keeps the trajectory tracking error in a sufficiently small vicinity of the origin to consider the desired trajectory as an artificial control in the PVTOL aircraft system with an inverted pendular load, where new auxiliary controllers are proposed. The PVTOL aircraft system with an inverted pendular load is tangentially linearized around the unstable equilibrium point, and the linearized system is used to design the second and third control loops. The second loop is devoted to trajectory tracking for take-off and landing maneuvers, and the third loop explores the differential flatness property in the horizontal displacement and pendulum’s angle system in order to stabilize both coordinates at the same time. A diagram of the proposed control scheme is depicted in Figure 2.
The following assumptions are considered for the control design.
Assumption 1.
The time derivatives of the desired trajectory for α angle are not well known.
Assumption 2.
Only linear and angular positions are available for measurement.

3.1. Control for α Angle

Dynamics for the α angle are described by Equations () and (). Where disturbance d α is caused by crosswind and is uniformly absolutely bounded, it is assumed that there exists a finite constant such that d α < δ α . The control objective is to keep α close to the desired trajectory, α d , which will be calculated by the following control loops; consequently, only the desired trajectory is known, but its time derivatives are not. Defining the trajectory tracking error e α = α d α and expressing the dynamical trajectory tracking error in terms of state variables e 1 α = e α and e 2 α = e ˙ α , the dynamical equations of the trajectory tracking error are expressed as follows:
e ˙ 1 α = e 2 α
e ˙ 2 α = α ¨ d u 2 D α .
By Assumption 1 α ¨ d is not well known, so it will be taken as disturbance, D α d = α d ¨ D α . It is clear that the nondisturbed system, i.e., D α d = 0 , (13), (14) is a pure integrator chain.
In order to reject the disturbance D α d , and to achieve Assumption 2, the following ESO is proposed for simultaneous estimation of e 2 α and the disturbance signal D α d :
e ^ ˙ 1 α = e ^ 2 α + λ 1 α ( e 1 α e ^ 1 α ) e ^ ˙ 2 α = u 2 + λ 2 α ( e 1 α e ^ 1 α ) + D ^ 1 α d D ^ ˙ 1 α d = D ^ 2 α d + λ 3 α ( e 1 α e ^ 1 α ) D ^ ˙ 2 α d = λ 4 α ( e 1 α e ^ 1 α ) .
The observation error of the state variable O α = e 1 α e ^ 1 α is subject to the following linearly dominant reconstruction error dynamics:
O α ( 4 ) + λ 1 α O α ( 3 ) + λ 2 α O ¨ α + λ 3 α O ˙ α + λ 4 α O α = ξ α .
The term ξ α represents the total uniformly absolutely bounded disturbance due to external disturbances and internal perturbations, i.e., ξ α < δ 1 ξ . A necessary and sufficient condition for having the observation error O α ultimately, uniformly, convergent toward a sufficiently small neighborhood of the acceleration estimation error phase space [3,33], consists in choosing the observer control parameters, λ j α j = 1 , 2 , 3 , 4 , such that the characteristic polynomial associated with the linear dominant dynamics is Hurwitz, making the linear injection error dynamics stable. The controller is thus given by
u 2 = k α 1 e 1 α + k α 2 e ^ 2 α D ^ 1 α d ,
where e 1 α is measurable and, is used instead of its estimate. The trajectory tracking error dynamics of the closed loop is given by
e ¨ α + k α 2 e ˙ α + k α 1 e α = η α ,
where η α = D α d D ^ α d < K α . The control parameters k α 1 and k α 2 are chosen such that a stable second-order dynamics with characteristic polynomial s 2 + 2 ω 0 ξ 0 s + ω 0 2 , being ω n , ξ 0 R 2 , is matched. In this case, the control parameters are chosen as k α 1 = ω 0 2 and k α 2 = 2 ω 0 ξ 0 . This means that angle α tends to be close to the desired trajectory α d in a sufficiently small neighborhood after a finite time t α , so
lim t e α < K α .

Linearized PVTOL Aircraft System with an Inverted Pendular Load

According to [16], after applying control u 2 (16) and a later finite time t α , the trajectory tracking error, e α , is sufficiently close to zero; conseuently, it is possible to replace α for α d and use it as artificial control input [10,22,23] for the PVTOL aircraft system with an inverted pendular load, where the error between α and α d is taken as a disturbance input in the subsequent control loop. Proposing control inputs u 1 and α d as
u 1 = r 1 2 + r 2 2 α d = arctan ( r 2 / r 1 ) ,
the PVTOL aircraft system with an inverted pendular load is expressed as follows:
x ¨ v = r 2 ( 2 m v + m p + m p cos ( 2 θ ) + m p r 1 sin ( 2 θ ) 2 l p m v m p s i n ( θ ) θ ˙ 2 2 m v ( m v + m p ) + D x y ¨ v = 2 g m v ( m v + m p ) + ( 2 m v + m p ) r 1 m p r 1 cos ( 2 θ ) 2 m v ( m v + m p )
+ m p r 2 sin ( 2 θ ) + l p m v m p cos ( θ ) θ ˙ 2 2 m v ( m v + m p ) + D y
θ ¨ = r 2 cos ( θ ) + r 1 sin ( θ ) l p m v + D θ
where
D x = d x ( m v + m p ) r 2 + m p sin ( θ ) ( d y r 1 cos ( θ ) + d x r 2 sin ( θ ) ) m v ( m v + m p ) r 1 2 + r 2 2
D y = d y r 1 ( m v + m p m p cos ( θ ) 2 ) + d x m p r 2 cos ( θ ) sin ( θ ) m v ( m v + m p ) r 1 2 + r 2 2
D θ = d x r 2 cos ( θ ) + d y r 1 sin ( θ ) l p m v r 1 2 + r 2 2 .
Equations (19)–(21) represent the reduced dynamical model of the PVTOL aircraft system with an inverted pendular load with r 1 and r 2 as control inputs.
A desired equilibrium point is given by: x v = 0 , y v = 0 , x ˙ v = 0 , y ˙ v = 0 , θ = 0 , θ ˙ = 0 , r 1 = g ( m v + m p ) , r 2 = 0 ; consequently, defining state variables x 1 = x , x 2 = x ˙ , y 1 = y , y 2 = y ˙ , θ 1 = θ , θ 2 = y ˙ and performing a tangent linearization of the system (19)–(21) [3], the resulting dynamics, ignoring disturbances, can be described in following form:
x ˙ 1 δ x ˙ 2 δ y ˙ 1 δ y ˙ 2 δ θ ˙ 1 δ θ ˙ 2 δ = 0 1 0 0 0 0 0 0 0 0 g m p m v 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 g ( m v + m p ) l p m v 0 x 1 δ x 2 δ y 1 δ y 2 δ θ 1 δ θ 2 δ + 0 0 0 1 m v 0 0 1 m v + m p 0 0 0 0 1 l p m v r 1 δ r 2 δ ,
where x 1 δ x 2 δ y 1 δ y 2 δ θ 1 δ θ 2 δ are the incremental state variables of the linear system, and r 1 δ r 2 δ are the incremental control inputs. The linear system (25) can be separated into two parts: the height subsystem ( y 1 δ ) and the horizontal pendulum subsystem ( x 1 δ , θ 1 δ ) as
y ˙ 1 δ y ˙ 3 δ = 0 1 0 0 y ˙ 1 δ y ˙ 3 δ + 0 1 m v + m p r 1 δ
x ˙ 1 δ x ˙ 2 δ θ ˙ 1 δ θ ˙ 2 δ = 0 1 0 0 0 0 g m p m v 0 0 0 0 1 0 1 g ( m v + m p ) l p m v 0 x 1 δ x 2 δ θ 1 δ θ 2 δ + 0 1 m v 0 1 m v + m p r 2 δ .
The subsystem (26) is the linear dynamical equation of the PVTOL’s height with incremental control input r 1 δ . It is a fully actuated and controllable system; the control law will be proposed to accomplish a take-off and landing maneuver with an ADRC to compensate the nonlinearities neglected in the linearization process and external disturbances.
Similarly, the subsystem (27) is the PVTOL’s horizontal position and the pendulum’s angle linear dynamics with control input r 2 δ . It is an underactuated system. A control law by differential flatness approach will be developed for a rest-to-rest maneuver, and nonlinearities and external disturbances will be compensated by another ADRC. Both control proposals will be acting at same time to accomplish the control objectives.

3.2. Height Control

As mentioned earlier, system (26) is the linear dynamical equation of the PVTOL’s height. The control objective is to achieve a take-off and landing maneuver. Consequently, we define the trajectory tracking error e y = y d y 1 δ . Let us express the dynamical equation for trajectory tracking error in terms of state variables e 1 y = e y , and e 2 y = e ˙ y , and the following dynamical equations are obtained:
e ˙ 1 y = e 2 y
e ˙ 2 y = r 1 δ m v + m p D y d .
By Assumption 1, y ¨ d is not measurable, so it will be taken as disturbance. D y d represents the total disturbances by high-order terms neglected in linearization process, the possibly unmodeled dynamics and external disturbances by crosswind (23). D y d is considered uniformly absolutely bounded so D y d < δ y , with δ y > 0 .
In order to reject the disturbance D y d , and to achieve Assumption 2, the following ESO is proposed for a simultaneous estimation of e 2 y and disturbance signal D y d ,
e ^ ˙ 1 y = e ^ 2 y + λ 1 y ( e 1 y e ^ 1 y ) e ^ ˙ 2 y = u 2 + λ 2 y ( e 1 y e ^ 1 y ) + D ^ 1 y d D ^ ˙ 1 y d = D ^ 2 y d + λ 3 y ( e 1 y e ^ 1 y ) D ^ ˙ 2 y d = λ 4 y ( e 1 y e ^ 1 y ) ,
where the state e 1 y is measurable. The observation error of the state variable O y = e 1 y e ^ 1 y generates the following linear reconstruction error dynamics:
O y 4 + λ 1 y + O y 3 + O y 2 + O y = ξ y .
The term ξ y represents the total uniformly, absolutely bounded disturbance due to external and internal disturbances, i.e., ξ y < δ 2 ξ . To force the observation error O y to be ultimately, uniformly, convergent toward a sufficiently small neighborhood of the acceleration estimation error phase space, let us choose the observer control parameters, λ j y j = 1 , 2 , 3 , 4 , such that the dominant linear injected error dynamics becomes stable (Hurwitz). The controller is thus given by
r 1 δ = ( m v + m p ) ( k y 1 e 1 y k y 2 e ^ 2 y D ^ 1 y d ) .
The trajectory tracking error dynamics of the closed loop is given by
e ¨ y + k y 2 e ˙ y + k y 1 e y = η y .
where η y = D y d D ^ y d < K y . The control parameters k y 1 and k y 2 are chosen such that a stable second-order dynamics with characteristic polynomial s 2 + 2 ω 1 ξ 1 s + ω 1 2 , being ω 1 , ξ 1 R 2 , is matched. In this case, the control parameters are chosen as k y 1 = ω 1 2 and k y 2 = 2 ω 1 ξ 1 . This means that the height y tends to be close to the desired trajectory y d in a sufficiently small neighbourhood after a finite time t y , so
lim t e y < K y
lim t r 1 δ < K r 1 ,
where K r 1 > 0 .

3.3. Control for Horizontal Displacement and Pendulum’s Angle

Considering the linear subsystem ( x , θ ), described in Equation (27), can be expressed in a compact form as
X ˙ h p = A h p X h p + B h p r 2 δ
with
A h p = 0 1 0 0 0 0 g m p m v 0 0 0 0 1 0 1 g ( m v + m p ) l p m v 0 ; B h p = 0 1 m v 0 1 m v + m p .
Moreover, X h p = ( x 1 δ , x 2 δ ) and r 2 δ as control input and according to [34,35], if the Kalman controllability matrix C k = [ B h p A h p B h p A h p 2 B h p A h p 3 B h p ] is not singular, the pair ( A h p , B h p ) is controllable and hence the system is flat with flat output computed as
F = [ 0 0 0 1 ] C k 1 X h p .
In order to simplify the calculation of the flat output time derivatives, flat output can be selected as
F = Γ [ 0 0 0 1 ] C k 1 X h p
with Γ = g l p 2 m is a constant; defining the row vector
c f = Γ [ 0 0 0 1 ] C k 1
where flat output F and and a finite number of its time derivatives can be obtained with the following equation:
F F ˙ F ¨ F ( 3 ) = c f c f A c f A 2 c f A 3 X h p .
Calculating (35), the flat output and its time derivative are completely parametrized in terms of the system variables and its time derivatives as
F = θ 1 δ x 1 δ l p F ˙ = θ 2 δ x 2 δ l p F ¨ = g θ 1 δ l p F ( 3 ) = g θ 2 δ l p .
The flat output fourth-order time derivative is obtained as follows:
F ( 4 ) = g ( r 2 δ + g ( m + m p ) θ 1 δ ) l p 2 m v .
Notice that the second and third time derivatives of flat outputs are computed by linear and angular velocities of θ δ and x δ and under Assumption 2, the velocities are not available. This restriction will be compensated in the control design procedure.
In the same way, all state variables can be parameterized by differential functions of the flat output F and can be computed by using Equation (36), as follows:
x 1 δ = l p F + l p 2 F ¨ g x 2 δ = l p F ˙ + l p 2 F ( 3 ) g θ 1 δ = l p F ¨ g θ 2 δ = l p F ( 3 ) g r 2 δ = l p 2 m v F ( 4 ) + g ( m v + m p ) l p F ¨ g .
The linearized system is clearly equivalent to the following input–output model
F ( 4 ) = g l p F ¨ ( m v + m p ) g r 2 δ l p 2 m v .
From the differential parametrization (36), it is clear that the tangent system naturally decomposes into a cascade connection of two independent blocks. The first is controlled by the incremental input r 2 δ with the output given by flat output F ¨ which coincides with a constant factor g l p , with incremental angular position of the pendulum θ 1 δ , i.e., F ¨ = θ 1 δ g l p . Then, the signal θ 1 δ g l p acts as an auxiliary control input to the second block, which consists of a chain of two integrators rendering the differential flat outputs F ˙ and F. The last variable F is the output of the second block and the output to be controlled for the overall system. This cascading property simplifies and decouples the observer design task in the flatness-based ADRC scheme.
Considering the system (38), we have the following simplified perturbed model for nonlinear subsystem ( x , θ )
F ( 4 ) = g r 2 δ l p 2 m v + ξ x θ ,
where ξ x θ represents state-dependent expressions, all the higher-order terms (h.o.t) neglected by the linearization, the possibly unmodeled dynamics, and external unknown disturbances affecting the system. All these uncertain terms are lumped into a single time-varying function represented by ξ x θ that is defined as
ξ x θ = g l p F ¨ ( m v + m p ) l p 2 m v + h . o . t .
defining the auxiliary control input r 2 δ as follows:
r 2 δ = v f l p 2 m v g ,
where v f is a new control input. Let F i = F ( i ) , i = 1 , 2 , 3 , 4 , the flat perturbed state space model is given by
F ˙ 1 = F 2 F ˙ 2 = F 3 F ˙ 3 = F 4 F ˙ 4 = v f + ξ x θ .
At this point, the cascading property is used, which implies that we are to view the previous system as the connection of two subsystems (see Figure 3). Notice that F 3 corresponds to a measurable position, which can be measurable, and it is the known input to the second-order pure integration system
F ˙ 1 = F 2 F ˙ 2 = F 3
and the remainder system is given by
F ˙ 3 = F 4 F ˙ 4 = v f + ξ x θ .
The perturbation term ξ x θ can be expressed in terms of the input, the output, and their algebraic combination. Then it is algebraically observable [36]. Next, let us propose an instantaneous virtual evolution model of a time-polynomial nature for such a time-varying function ξ x θ , denoted by z, and adopt, say, the following third-order time polynomial model for ξ x θ , i.e., z ( 4 ) = 0 . Then, with z ( t ) = z 1 , the the flat perturbed state space model (42) can be rewritten as follows:
F ˙ i = F i + 1 , i = 1 , 2 , 3 , F ˙ 4 = v f + z 1 z ˙ j = z j + 1 , j = 1 , 2 , 3 , z ˙ 4 = 0 .
According to [33,37], a set of coupled, high-gain, extended linear Luenberger observers, for the simultaneous estimation of the phase variables associated with the flat output and the time-polynomial approximation variable, can be proposed as follows:
F ^ ˙ 1 = F 2 + k 1 ( F 1 F ^ 1 ) F ^ ˙ 2 = F 3 + k 2 ( F 1 F ^ 1 ) F ^ ˙ 3 = F 4 + λ 1 ( F 3 F ^ 3 ) F ^ ˙ 4 = v f + λ 2 ( F 3 F ^ 3 ) + z ^ 1 z ^ ˙ 1 = z ^ 2 + λ 3 ( F 3 F ^ 3 ) z ^ ˙ 2 = z ^ 2 + λ 4 ( F 3 F ^ 3 ) z ^ ˙ 3 = z ^ 2 + λ 5 ( F 3 F ^ 3 ) z ^ ˙ 4 = λ 6 ( F 3 F ^ 3 ) .
The observation error, O F 1 = F 1 F ^ 1 , of the incremental flat output, generates the following linear injected estimation error dynamics:
O ¨ F 1 + k 1 O ˙ F 1 + k 2 O F 1 = 0 .
An appropriate choice of observer parameters, k 1 , k 2 such that the characteristic polynomial s 2 + k 1 s + k 2 is Hurwitz renders an asymptotically, exponentially decreasing estimation error state. The tracking error velocity for the flat output O ˙ F 1 is, thus, accurately estimated for feedback purposes.
In the same manner, considering the observation error O F 3 = F 3 F ^ 3 of the flat output acceleration tracking error, it generates the following dominantly linear reconstruction error dynamics:
O F 3 ( 6 ) + λ 1 O F 3 ( 5 ) + λ 2 O F 3 ( 4 ) + λ 3 O F 3 ( 3 ) + λ 4 O ¨ F 3 + λ 5 O ˙ F 3 + λ 6 O F 3 = ξ x θ ( 4 ) .
As stated in the previous ADRC designs, a necessary and sufficient condition for having the incremental flat output acceleration estimation error O F 3 and its time derivatives ultimately, uniformly, converge toward a small-as-desired neighborhood of the acceleration estimation error phase space is that ξ x θ ( 4 ) be uniformly, absolutely bounded. An appropriate choice of observer gains coefficients: λ i , i = 1 , 2 , . . . , 6 , such that placing the poles of the associated linear homogeneous system sufficiently far into the left half of the complex plane, renders a uniformly, asymptotically convergent estimation error, O F 3 , toward an arbitrary small vicinity of the origin along with a finite number of its time derivatives.
The control input may then be readily synthesized with an active disturbance-cancelling strategy for the uncertain input ξ x θ , in terms of the estimates value z ^ 1 and, for feedback purposes, of the estimated time derivatives associated with the incremental flat outputs F 1 and F 3 , both measurable. Then, the following control input v f is proposed:
v f = K F 1 F ^ 2 F 3 F ^ 4 z ^ 1 ,
where F 1 and F 3 are used in a way that can be calculated through the variables x and θ , which are assumed to be measurable, F ^ 2 and F ^ 4 are estimates of F 2 and F 4 , and K R 1 x 4 is the gain matrix. Notice that the coefficients in K matrix must be chosen in accordance with the fact that, asymptotically, the flat output is being approximately governed by the differential equation
F 1 ( 4 ) + K 1 F 1 ( 3 ) + K 2 F ¨ 1 + K 3 F ˙ 1 + K 4 F 1 = ξ x θ z ^ 1 ,
where the set of coefficients K 1 , K 2 , K 3 , K 4 should render Hurwitz the following characteristic polynomial: s ( 4 ) + K 1 s ( 3 ) + K 2 s 2 + K 3 s + K 4 . Finally, the diagram of the entire control scheme can be seen in Figure 4.

4. Numerical Simulations

To test the controller’s performance, numerical simulations were carried out by using MATLAB–Simulink, and the results were obtained based on the numerical method of Runge–Kutta of the fourth order with a fixed step of 0.01 s. The physical parameters of the system are m v = 0.9 kg, m p = 0.3 kg, l p = 0.5 m, g = 9.81 m/s 2 . The desired trajectory in height position y d was implemented by the use of gauss2mf() in Matlab software with the following parameters: σ 1 = 1.5 , μ 1 = 10 , σ 1 = 1.5 , and μ 1 = 20 .
The control parameters for u 2 (16) are proposed as k α 1 = 3.16 and k α 2 = 2.7 , and the observers gain parameters (15) are calculated as follows:
λ 1 α = 4 ϵ ω n λ 2 α = 2 ω n 2 + 4 ϵ 2 ω n 2 λ 3 α = 4 ϵ ω n 3 λ 4 α = ω n 4
with ω n = 20 and ϵ = 0.7 .
For the case of the controller r 1 δ (31), the control parameters are set to be k y 1 = 3.16 and k y 2 = 2.7 , and the observers gain parameters (30) are calculated according to the following expression:
λ 1 y = 4 ϵ ω n λ 2 y = 2 ω n 2 + 4 ϵ 2 ω n 2 λ 3 y = 4 ϵ ω n 3 λ 4 y = ω n 4
with ω n = 10 and ϵ = 0.7 .
In the case of the controller v f (46), the gain matrix K is proposed as
K = 1 3 4.2 3 ,
and the observers gain coefficients are calculated with the following equations:
k 1 = 2 ϵ ω n k 2 = ω n 2
with ω n = 10 and ϵ = 0.7
λ 1 = 6 ϵ ω n λ 2 = 3 ω n 2 + 12 ϵ 2 ω n 2 λ 3 = 12 ϵ ω n 3 + 8 ϵ 3 ω n 3 λ 4 = 3 ω n 4 + 12 ϵ 2 ω n 4 λ 5 = 6 ϵ ω n 5 λ 6 = ω n 6
with ω n = 0.5 and ϵ = 2 .
To show the effectiveness of the control scheme proposed in this work, a sliding mode control law was implemented to a comparison test against the robust differential flatness approach control scheme. The design of the sliding mode control is based on the linear approximation of the whole system that has the form (see Appendix A)
x ˙ L = A L x L + B L u L ,
where A L and B L are the linear constant matrices of the tangent linearized system, x L = [ x v , x v ˙ , y v , y v ˙ , α , α ˙ , θ , θ ˙ ] are the state variables of the system, and u L = u 1 u 2 is the set of control inputs of the linearized system. The proposed sliding mode control is
u 1 = 5 sign ( K L 1 x L )
u 2 = 15 sign ( K L 2 x L ) ,
where K L 1 and K L 2 are gain matrices designed by LQR technique based on the linearized system and sign ( x ) is the sign function defined as follows:
sign ( x ) = 1 i f x > 0 1 i f x < 0 0 i f x = 0
and the matrices K L 1 and K L 2 are
K L 1 = 0 0 1 1.8 0 0 0 0
K L 2 = 119 15.5 0 0 1 2.3 141 28 .
Both control laws were implemented in order to contrast the performance of the proposed control scheme with external disturbances caused by crosswind [30]. The crosswind dynamics are governed by the values of d x , d y and d a . These values are proposed as random bounded numbers with a Gaussian distribution as follows:
d x = 0.5 rand ( t ) d y = rand ( t ) d α = 0.5 rand ( t ) ,
where rand ( ) is a function that provides a random number in the interval ( 1 , 1 ) with the Gaussian distribution.
Figure 5 shows the behavior of α with the proposed control scheme ( ) and with sliding mode control ( · ). As can be seen, at the beginning of the simulation, the proposed control scheme shows a lesser overshooting effect than the sliding mode control response, and, in general, the error amplitude in the proposed control scheme is more reduced than the sliding mode case. This fact is shown in Figure 6 where an integral square error performance index of α is calculated for both cases.
Figure 7 shows the control input u 2 in both cases, for the proposed control scheme ( ) and the sliding mode control ( · ). The chattering problem in sliding mode control does not appear, but a peaking phenomenon arises in the transient response of the sliding mode controller and a larger control amplitude of the proposal is noticed in the transient behavior due to the high gain nature of the controller. In Figure 7c, the energy consumption performance of u 2 , in both cases, is evaluated by calculating the integral square of u 2 . In this, it is evident that the behavior of the proposed controller u 2 demands less energy than the discontinuous proposal in steady state.
Figure 8 shows the behavior of y v with the proposed control scheme ( ) and with sliding mode control ( · ), including the reference trajectory y d (−). The proposal has a better transient behavior than the sliding mode, which can be seen in Figure 9. After a while, the tracking trajectory error in both cases are similar in amplitude, which shows that both controllers are suitable for the task. Controller u 1 ’s behavior with the proposed control scheme and with sliding mode are shown in Figure 10; in this case, the minimum and maximum amplitude of the proposed control is similar with respect to the sliding mode behavior. However, the discontinuous nature of the sliding mode control makes the integral square index of u 1 show less energy consumption of the proposal while ensuring competitive error results.
The behavior of the coordinate x v is shown in Figure 11 and the behavior of the pendulum’s angle θ is shown in Figure 12. Physically, controller v f moves the PVTOL system horizontally to keep the θ angle close to zero. The controller returns a horizontal displacement to zero at the same time to try to hold θ close to zero. In both coordinates, it is easy to see that the overshoot of the control proposal is less than the one exhibited by the sliding mode control but, after the overshoot phenomenon, the sliding mode control case has a lesser amplitude than the proposed control scheme.
The dynamics of the observation errors in the designed ESO ( O α , O y , O F 3 and O F 5 ), can be seen in Figure 13, Figure 14, Figure 15 and Figure 16 respectively. In all cases, it is clear that the dynamics of observation errors converge to a close region of zero as expected.

5. Conclusions

This work presents a control scheme for a PVTOL aircraft system with an inverted pendular load. This control scheme is divided into three control loops. The first closed-loop control is focused on a trajectory tracking for α angle in order to use the desired trajectory α d as an artificial control. An ADRC controller is proposed in order to reject internal and external disturbances and unknown derivatives as well. After defining auxiliary control inputs, the PVTOL aircraft system with an inverted pendular load is tangentially linearized around the unstable equilibrium point. The linearized PVTOL aircraft system with an inverted pendular load is divided into two subsystems: the height system and the horizontal pendulum system.
In the height system, an ADRC controller is proposed in terms of the trajectory tracking error to achieve take-off and landing maneuvers in the presence of external disturbances. The third closed-loop control is designed by means of exploiting the flatness associated with the linearized model and via an extended state observer-based linear ADRC.
This control scheme achieves a competitive robustness in the presence of external disturbances, due to crosswind, and to reject the nonlinearities of the PVTOL aircraft system with an inverted pendular load model, neglected in the linearization process, which means that the proposed controller behaves well even in this undesirable realistic scenario. It is important to note that the scheme achieves greatly competitive numerical simulation results when compared with another effective robust controller, the sliding mode control. Note that the performance index of both controllers was computed to compare them, and the outcome revealed that our strategy has a better performance than the sliding mode controller.

Author Contributions

Conceptualization, C.A.V.R., A.L.-J., N.B.L.-C., B.E.C.-G., D.M.-V., and O.G.-F.; methodology, C.A.V.R., A.L.-J., N.B.L.-C., B.E.C.-G., D.M.-V., and O.G.-F.; software, C.A.V.R., A.L.-J., N.B.L.-C., B.E.C.-G., D.M.-V., and O.G.-F.; validation, C.A.V.R., A.L.-J., N.B.L.-C., B.E.C.-G., D.M.-V., and O.G.-F.; writing—original draft preparation, C.A.V.R., A.L.-J., and O.G.-F.; writing—review and editing, C.A.V.R., A.L.-J., N.B.L.-C., B.E.C.-G., D.M.-V., and O.G.-F.; visualization, A.L.-J. and N.B.L.-C.; funding acquisition, A.L.-J., N.B.L.-C., B.E.C.-G., and O.G.-F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Secretaría de Investigación y Posgrado del Instituto Politécnico Nacional (SIP-IPN) under grant numbers 20220542, 20220632, 20220623 and 20220633.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

C. Alejandro Villaseñor Rios thanks the support from the CONACYT.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Design of the Sliding Mode Control

The design of sliding mode control is based on linear approximation of the whole system that has the form
x ˙ L = A L x L + B L u L
where A L and B L ate the linear constant matrices of the linearized system by Taylor series, x L = ( y δ , y ˙ δ , α δ , α ˙ δ , x δ , x ˙ δ , θ δ , θ ˙ δ ) are the incremental variables states of the linearized system and u L = ( u 1 , u 2 ) is the set of control inputs of the linearized system, and the matrices of the tangent linearized system (A1) are:
A L = 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 g ( m v + m p ) m 0 0 0 g ( m v + m p ) m v 0 0 0 0 0 0 0 0 1 0 0 g ( m v + m p ) l p m 0 0 0 g ( m v + m p ) l p m v 0 B L = 0 0 1 m v + m p 0 0 0 0 1 0 0 0 0 0 0 0 0
The height dynamics is decoupled of the rest of the system with u 1 as control input, so; for which it is possible to separate the dynamics in the linear system in the following way.
x L h ˙ L = A L h x L h + B L h u 1
A L h = 0 1 0 0 B L h = 0 1 m v + m p
x L θ ˙ L = A L θ x L θ + B L θ u 2
A L θ = 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 g ( m v + m p ) m 0 0 0 g ( m v + m p ) m 0 0 0 0 0 0 1 g ( m v + m p ) l p m 0 0 0 g ( m V + m p ) l p m 0 B L θ = 0 1 0 0 0 0
where x L h = ( y δ , y ˙ δ ) = ( y 1 , y 2 ) is the state variable for height dynamics with A L h and B L h linear matrices of the system and u 1 as control input; x L θ = ( α δ , α ˙ δ , x δ , x ˙ δ , θ δ , θ ˙ δ ) = ( α 1 , α 2 , x 1 , x 2 , θ 1 , θ 2 ) is the state variable for ( α δ ) , x δ , θ δ ) dynamics with A L θ and B L θ linear matrices and u 2 as control input. The sliding mode control proposed is
u 1 = v 1 s i g n ( K L 1 x L h ) + g ( m v + m p )
u 2 = v 2 s i g n ( K L 2 x l θ )
where v 1 and v 2 are control gains and K L 1 and K L 2 are gain matrices designed by the LQR technique based on the linearized system with lqr() MATLAB function as follows
K L 1 = l q r ( A L h , B L h , e y e ( 2 ) , 1 )
K L 2 = l q r ( A L θ , B L θ , e y e ( 6 ) , 1 )
s i g n ( x ) is the sign function defined as follows
s i g n ( x ) = 1 i f x > 0 1 i f x < 0 0 i f x = 0
The stability proof for whole system is divided in two steps: stability proof for (A2) system and stability proof for (A3) system.

Appendix A.1. Stability Proof for Height Dynamics

To proof the height dynamics (A2) stability in closed loop with the control u 2 (A5), let us take the following sliding mode surface
σ y = K L 1 x L h = K L 1 1 y 1 + K L 1 2 y 2
Defining the following Lyapunov function
V y = 1 2 σ y 2
taking the time derivative of V y yields
V ˙ y = σ y σ ˙ y = σ y ( K L 1 1 y 2 + K L 1 2 y ˙ 2 )
V ˙ y = σ y ( K L 1 1 y 2 + K L 1 2 ( u 1 δ m v + m p + ϵ y )
V ˙ y = σ y ( K L 1 1 y 2 + K L 1 2 ( v 1 s i g n ( σ y ) + g ( m v + m p ) m v + m p + ϵ y )
V ˙ y = σ y ( K L 1 1 y 2 + K L 1 2 ( C 1 s i g n ( σ y ) + g + ϵ y )
V ˙ y = σ y ( P y s i g n ( σ y ) + ξ y )
where ξ y represents the total disturbance in y v dynamics with other terms that may affect in stability, and must be bounded so ξ y L y ; P y is a new control parameter.
V ˙ y σ y ( P y s i g n ( σ y ) + L y
V ˙ y P y σ y + σ y L y P y σ y + σ y L y
if P y > L y then V ˙ y < 0 ; so, the system is asymptotically stable at origin.

Appendix A.2. Stability Proof for (α,x v, θ) Dynamics

To proof the ( α , x v , θ ) dynamics (A3) stability in closed-loop with control u 1 (A4), let us take the following sliding surface
σ = K L 2 x L 2 = K L 2 1 α 1 + K L 2 2 α 2 + K L 2 3 x 1 + K L 2 4 x 2 + K L 2 5 θ 1 + K L 2 6 θ 2
Defining the following Lyapunov function
V = 1 2 σ 2
taking the time derivative of V yields
V ˙ = σ σ ˙ = σ ( K L 2 1 α 2 + K L 2 2 α ˙ 2 + K L 2 3 x 2 + K L 2 4 x ˙ 2 + K L 2 5 θ 2 + K L 2 6 θ ˙ 2 )
V ˙ = σ ( K L 2 1 α 2 + K L 2 2 ( v 2 s i g n ( σ ) ) + K L 2 3 x 2 + K L 2 4 x ˙ 2 + K L 2 5 θ 2 + K L 2 6 θ ˙ 2 )
V ˙ = σ ( P θ s i g n ( σ ) + ξ θ
where ξ θ denotes the total disturbances in y v dynamics with other terms that may affect in stability; and must be bounded, ξ y L θ ; P y is a new control parameter.
V ˙ σ ( P θ s i g n ( σ ) + L θ
V ˙ P θ σ + σ L θ P θ σ + σ L θ
if P θ > L θ then V ˙ < 0 ; so, the equilibrium given by the surface is asymptotically stable at the origin.

References

  1. Boubaker, O. The inverted pendulum: A fundamental benchmark in control theory and robotics. In Proceedings of the International Conference on Education and e-Learning Innovations, Sousse, Tunisia, 1–3 July 2012; pp. 1–6. [Google Scholar]
  2. Shkolnik, A.; Tedrake, R. High-dimensional underactuated motion planning via task space control. In Proceedings of the 2008 IEEE/RSJ International Conference on Intelligent Robots and Systems, Nice, France, 22–26 September 2008; pp. 3762–3768. [Google Scholar]
  3. Ramírez-Neria, M.; Sira-Ramírez, H.; Garrido-Moctezuma, R.; Luviano-Juarez, A. Linear active disturbance rejection control of underactuated systems: The case of the Furuta pendulum. ISA Trans. 2014, 53, 920–928. [Google Scholar] [CrossRef] [PubMed]
  4. Gómez-Estern, F.; Van der Schaft, A.J. Physical damping in IDA-PBC controlled underactuated mechanical systems. Eur. J. Control. 2004, 10, 451–468. [Google Scholar] [CrossRef] [Green Version]
  5. Jabbari Asl, H.; Oriolo, G.; Bolandi, H. Output feedback image-based visual servoing control of an underactuated unmanned aerial vehicle. Proc. Inst. Mech. Eng. Part I J. Syst. Control. Eng. 2014, 228, 435–448. [Google Scholar] [CrossRef]
  6. Rodríguez-Seda, E.J.; Tang, C.; Spong, M.W.; Stipanović, D.M. Trajectory tracking with collision avoidance for nonholonomic vehicles with acceleration constraints and limited sensing. Int. J. Robot. Res. 2014, 33, 1569–1592. [Google Scholar] [CrossRef]
  7. Ramirez-Neria, M.; Ochoa-Ortega, G.; Lozada-Castillo, N.; Trujano-Cabrera, M.A.; Campos-Lopez, J.P.; Luviano-Juárez, A. On the robust trajectory tracking task for flexible-joint robotic arm with unmodeled dynamics. IEEE Access 2016, 4, 7816–7827. [Google Scholar] [CrossRef]
  8. Furut, K.; Ochiai, T.; Ono, N. Attitude control of a triple inverted pendulum. Int. J. Control 1984, 39, 1351–1365. [Google Scholar] [CrossRef]
  9. Haddad, N.K.; Chemori, A.; Belghith, S. Robustness enhancement of IDA-PBC controller in stabilising the inertia wheel inverted pendulum: Theory and real-time experiments. Int. J. Control 2018, 91, 2657–2672. [Google Scholar] [CrossRef] [Green Version]
  10. Fantoni, I.; Lozano, R. Global stabilization of the cart-pendulum system using saturation functions. In Proceedings of the 42nd IEEE International Conference on Decision and Control (IEEE Cat. No. 03CH37475), Maui, HI, USA, 9–12 December 2003; Volume 5, pp. 4393–4398. [Google Scholar]
  11. Aguilar-Ibañez, C.; Suárez-Castañón, M.S.; Gutiérres-Frias, O.O. The direct Lyapunov method for the stabilisation of the Furuta pendulum. Int. J. Control 2010, 83, 2285–2293. [Google Scholar] [CrossRef]
  12. Moeini, A.; Lynch, A.F.; Zhao, Q. A backstepping disturbance observer control for multirotor UAVs: Theory and experiment. Int. J. Control 2021. [Google Scholar] [CrossRef]
  13. Hua, M.D.; Hamel, T.; Morin, P.; Samson, C. Introduction to feedback control of underactuated VTOLvehicles: A review of basic control design ideas and principles. IEEE Control. Syst. Mag. 2013, 33, 61–75. [Google Scholar]
  14. Liu, C.; Pan, J.; Chang, Y. PID and LQR trajectory tracking control for an unmanned quadrotor helicopter: Experimental studies. In Proceedings of the 2016 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; pp. 10845–10850. [Google Scholar]
  15. Xiong, L.; Yu, Z.; Wang, Y.; Yang, C.; Meng, Y. Vehicle dynamics control of four in-wheel motor drive electric vehicle using gain scheduling based on tyre cornering stiffness estimation. Veh. Syst. Dyn. 2012, 50, 831–846. [Google Scholar] [CrossRef]
  16. Lozano Hernández, Y.; Gutiérrez Frías, O.; Lozada-Castillo, N.; Luviano Juárez, A. Control algorithm for taking off and landing manoeuvres of quadrotors in open navigation environments. Int. J. Control. Autom. Syst. 2019, 17, 2331–2342. [Google Scholar] [CrossRef]
  17. Azinheira, J.R.; Moutinho, A. Hover control of an UAV with backstepping design including input saturations. IEEE Trans. Control. Syst. Technol. 2008, 16, 517–526. [Google Scholar] [CrossRef]
  18. Das, A.; Lewis, F.; Subbarao, K. Backstepping approach for controlling a quadrotor using lagrange form dynamics. J. Intell. Robot. Syst. 2009, 56, 127–151. [Google Scholar] [CrossRef]
  19. Zheng, E.H.; Xiong, J.J.; Luo, J.L. Second order sliding mode control for a quadrotor UAV. ISA Trans. 2014, 53, 1350–1356. [Google Scholar] [CrossRef] [PubMed]
  20. Hauser, J.; Sastry, S.; Meyer, G. Nonlinear control design for slightly non-minimum phase systems: Application to V/STOL aircraft. Automatica 1992, 28, 665–679. [Google Scholar] [CrossRef]
  21. Escobar, J.C.; Lozano, R.; Bonilla Estrada, M. PVTOL control using feedback linearisation with dynamic extension. Int. J. Control 2021, 94, 1794–1803. [Google Scholar] [CrossRef]
  22. Zavala-Rio, A.; Fantoni, I.; Lozano, R. Global stabilization of a PVTOL aircraft model with bounded inputs. Int. J. Control 2003, 76, 1833–1844. [Google Scholar] [CrossRef]
  23. López-Araujo, D.J.; Zavala-Río, A.; Fantoni, I.; Salazar, S.; Lozano, R. Global stabilisation of the PVTOL aircraft with lateral force coupling and bounded inputs. Int. J. Control 2010, 83, 1427–1441. [Google Scholar] [CrossRef]
  24. Aguilar-Ibanez, C.; Suarez-Castanon, M.S.; Meda-Campaña, J.; Gutierrez-Frias, O.; Merlo-Zapata, C.; Martinez-Castro, J.A. A simple approach to regulate a pvtol system using matching conditions. J. Intell. Robot. Syst. 2020, 98, 511–524. [Google Scholar] [CrossRef]
  25. Munoz, L.E.; Santos, O.; Castillo, P. Robust nonlinear real-time control strategy to stabilize a PVTOL aircraft in crosswind. In Proceedings of the 2010 IEEE/RSJ International Conference on Intelligent Robots and Systems, Taipei, Taiwan, 18–22 October 2010; pp. 1606–1611. [Google Scholar]
  26. Hehn, M.; D’Andrea, R. A flying inverted pendulum. In Proceedings of the 2011 IEEE International Conference on Robotics and Automation, Shangai, China, 9–13 May 2011; pp. 763–770. [Google Scholar]
  27. Krafes, S.; Chalh, Z.; Saka, A. Vision-based control of a flying spherical inverted pendulum. In Proceedings of the 2018 4th International Conference on Optimization and Applications (ICOA), Mohammedia, Morocco, 26–27 April 2018; pp. 1–6. [Google Scholar]
  28. Chen, H.; Yang, Y.; Sun, J. Improved Genetic Algorithm Based Optimal Control for A Flying Inverted Pendulum. In Proceedings of the 2019 3rd International Conference on Electronic Information Technology and Computer Engineering (EITCE), Xiamen, China, 18–20 October 2019; pp. 1428–1432. [Google Scholar]
  29. Nayak, A.; Banavar, R.N.; Maithripala, D. Stabilizing a spherical pendulum on a quadrotor. Asian J. Control 2020, 24, 1112–1121. [Google Scholar] [CrossRef]
  30. Aguilar-Ibanez, C.; Sira-Ramirez, H.; Suarez-Castanon, M.S.; Garrido, R. Robust trajectory-tracking control of a PVTOL under crosswind. Asian J. Control 2019, 21, 1293–1306. [Google Scholar] [CrossRef]
  31. Zhang, C.; Hu, H.; Gu, D.; Wang, J. Cascaded control for balancing an inverted pendulum on a flying quadrotor. Robotica 2017, 35, 1263–1279. [Google Scholar] [CrossRef] [Green Version]
  32. Ramírez-Neria, M.; Sira-Ramírez, H.; Garrido-Moctezuma, R.; Luviano-Juárez, A. Active disturbance rejection control of the inertia wheel pendulum through a tangent linearization approach. Int. J. Control. Autom. Syst. 2019, 17, 18–28. [Google Scholar] [CrossRef]
  33. Ramírez-Neria, M.; Sira-Ramírez, H.; Garrido-Moctezuma, R.; Luviano-Juárez, A. On the linear control of underactuated nonlinear systems via tangent flatness and active disturbance rejection control: The case of the ball and beam system. J. Dyn. Syst. Meas. Control 2016, 138, 104501. [Google Scholar] [CrossRef]
  34. Sira-Ramirez, H.; Agrawal, S.K. Differentially Flat systems; Crc Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  35. Sira-Ramírez, H.; Luviano-Juárez, A.; Ramírez-Neria, M.; Zurita-Bustamante, E.W. Active Disturbance Rejection Control of Dynamic Systems: A Flatness Based Approach; Butterworth-Heinemann: Oxford, UK, 2018. [Google Scholar]
  36. Diop, S.; Fliess, M. Nonlinear observability, identifiability, and persistent trajectories. In Proceedings of the 30th IEEE Conference on Decision and Control, Brighton, UK, 11–13 December 1991; pp. 714–719. [Google Scholar]
  37. Ramírez-Neria, M.; Sira-Ramírez, H.; Garrido-Moctezuma, R.; Luviano-Juárez, A.; Gao, Z. Active Disturbance Rejection Control for Reference Trajectory Tracking Tasks in the Pendubot System. IEEE Access 2021, 9, 102663–102670. [Google Scholar] [CrossRef]
Figure 1. PVTOL aircraft system with an inverted pendular load.
Figure 1. PVTOL aircraft system with an inverted pendular load.
Machines 10 00595 g001
Figure 2. Schematics of the control scheme.
Figure 2. Schematics of the control scheme.
Machines 10 00595 g002
Figure 3. Cascade structure of ( x , θ ) sub−system.
Figure 3. Cascade structure of ( x , θ ) sub−system.
Machines 10 00595 g003
Figure 4. Diagram of the proposed control scheme.
Figure 4. Diagram of the proposed control scheme.
Machines 10 00595 g004
Figure 5. Behavior of α .
Figure 5. Behavior of α .
Machines 10 00595 g005
Figure 6. Integral square error index of α tracking.
Figure 6. Integral square error index of α tracking.
Machines 10 00595 g006
Figure 7. u 2 behavior results. (a) u 2 behavior (proposed control); (b) u 2 behavior (sliding mode control); (c) integral square of u 2 calculation.
Figure 7. u 2 behavior results. (a) u 2 behavior (proposed control); (b) u 2 behavior (sliding mode control); (c) integral square of u 2 calculation.
Machines 10 00595 g007
Figure 8. Behavior of y v with the proposed control scheme and with the sliding mode control in trajectory tracking task.
Figure 8. Behavior of y v with the proposed control scheme and with the sliding mode control in trajectory tracking task.
Machines 10 00595 g008
Figure 9. Tracking error behavior of e y with the proposed controller and with the sliding mode control.
Figure 9. Tracking error behavior of e y with the proposed controller and with the sliding mode control.
Machines 10 00595 g009
Figure 10. u 1 behavior for the implemented control schemes. (a) u 1 behavior (proposed control); (b) u 1 behavior (sliding mode control); (c) integral square index of u 1 .
Figure 10. u 1 behavior for the implemented control schemes. (a) u 1 behavior (proposed control); (b) u 1 behavior (sliding mode control); (c) integral square index of u 1 .
Machines 10 00595 g010
Figure 11. Behavior of x v for the implemented control schemes.
Figure 11. Behavior of x v for the implemented control schemes.
Machines 10 00595 g011
Figure 12. Behavior for θ for the implemented control schemes.
Figure 12. Behavior for θ for the implemented control schemes.
Machines 10 00595 g012
Figure 13. Observation error dynamics O α .
Figure 13. Observation error dynamics O α .
Machines 10 00595 g013
Figure 14. Observation error dynamics O y .
Figure 14. Observation error dynamics O y .
Machines 10 00595 g014
Figure 15. Observation error dynamics O F 3 .
Figure 15. Observation error dynamics O F 3 .
Machines 10 00595 g015
Figure 16. Observation error dynamics O F 5 .
Figure 16. Observation error dynamics O F 5 .
Machines 10 00595 g016
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Villaseñor Rios, C.A.; Luviano-Juárez, A.; Lozada-Castillo, N.B.; Carvajal-Gámez, B.E.; Mújica-Vargas, D.; Gutiérrez-Frías, O. Flatness-Based Active Disturbance Rejection Control for a PVTOL Aircraft System with an Inverted Pendular Load. Machines 2022, 10, 595. https://doi.org/10.3390/machines10070595

AMA Style

Villaseñor Rios CA, Luviano-Juárez A, Lozada-Castillo NB, Carvajal-Gámez BE, Mújica-Vargas D, Gutiérrez-Frías O. Flatness-Based Active Disturbance Rejection Control for a PVTOL Aircraft System with an Inverted Pendular Load. Machines. 2022; 10(7):595. https://doi.org/10.3390/machines10070595

Chicago/Turabian Style

Villaseñor Rios, Cesar Alejandro, Alberto Luviano-Juárez, Norma Beatriz Lozada-Castillo, Blanca Esther Carvajal-Gámez, Dante Mújica-Vargas, and Octavio Gutiérrez-Frías. 2022. "Flatness-Based Active Disturbance Rejection Control for a PVTOL Aircraft System with an Inverted Pendular Load" Machines 10, no. 7: 595. https://doi.org/10.3390/machines10070595

APA Style

Villaseñor Rios, C. A., Luviano-Juárez, A., Lozada-Castillo, N. B., Carvajal-Gámez, B. E., Mújica-Vargas, D., & Gutiérrez-Frías, O. (2022). Flatness-Based Active Disturbance Rejection Control for a PVTOL Aircraft System with an Inverted Pendular Load. Machines, 10(7), 595. https://doi.org/10.3390/machines10070595

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop