Next Article in Journal
A Multi-Scale Hybrid Attention Network for Sentence Segmentation Line Detection in Dongba Scripture
Previous Article in Journal
An Improved Incremental Procedure for the Ground Reaction Based on Hoek-Brown Failure Criterion in the Tunnel Convergence-Confinement Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model-Based Optimization Approach for PID Control of Pitch–Roll UAV Orientation

by
Orlando Arrieta
1,2,*,
Daniel Campos
1,
Javier Rico-Azagra
3,
Montserrat Gil-Martínez
3,*,
José D. Rojas
1,2 and
Ramon Vilanova
2,*
1
Instituto de Investigaciones en Ingeniería, Facultad de Ingeniería, Universidad de Costa Rica, San José 11501-2060, Costa Rica
2
Departament de Telecomunicació i d’Enginyeria de Sistemes, Escola d’Enginyeria, Universitat Autònoma de Barcelona, Bellaterra, 08193 Barcelona, Spain
3
Control Engineering Research Group, Electrical Engineering Department, University of La Rioja, 26004 Logroño, Spain
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(15), 3390; https://doi.org/10.3390/math11153390
Submission received: 3 May 2023 / Revised: 27 June 2023 / Accepted: 16 July 2023 / Published: 3 August 2023
(This article belongs to the Section Engineering Mathematics)

Abstract

:
This study focuses on the modeling and control optimization of the pitch–roll angles orientation in a quad-rotor UAV (unmanned aerial vehicle) drone system. To address the control challenge, PID (proportional–integral–derivative) controllers are used, which have been tuned through a model-based optimization approach. A non-linear model of the system is obtained and further validated experimentally. The model is implemented in Simulink to conduct model-based optimization. The controller design uses an IMC (internal model control) model-based approach to establish the initial parameters for a PID controller. In the subsequent design phase, a multiobjective approach is taken, incorporating weighted cost functions that account for IAE (integral absolute error) and ITAE (integral of time-weighted absolute error) forms, along with battery usage. The outcomes of this design are showcased through simulations and real-world platform testing.

1. Introduction

Unmanned aerial vehicles (UAVs) are currently becoming increasingly popular in a wide range applications [1]. Although they were developed mainly for military purposes, it has recently become obvious that there are many other areas where UAVs might prove useful. In the field of agriculture, for example [2], they may be used for field observations or for chemical distribution. They can patrol over wide forest areas to monitor forest health [3], or they can be used for traffic observation in cities. Their potential applications in innumerable fields have led to their widespread use, and as a result, new challenges have arisen [4]. Additionally, the technological maturity of the sector has made it possible to configure ad hoc UAV equipment at a moderate cost [5,6,7,8], which has favored its use in training and research activities.
From a scientific and technical point of view, UAVs are multidisciplinary systems of special interest in engineering and, specifically, for the practice of control engineering [9,10]. In particular, controlling a multi-rotor-type UAV is a non-linear problem with six degrees of freedom (6-DoF) and with underactuation, especially when the thrusters cannot be reoriented [11,12]. For this reason, to address this complexity, a hierarchical control structure is commonly employed. The upper layers focus on positioning, while the lower layers handle orientation [13,14]. In turn, adequate control of position [15] makes the adequate tracking of trajectories possible [16]. However, the success of positioning is subject to adequate orientation control [17,18]. It is from this perspective that the present work focuses on this orientation problem.
UAV systems can be controlled using different methodologies, such as data-driven methodologies [19,20,21] or optimization techniques applied to proportional–integral–derivative (PID) controllers [22,23].
We use the experimental setup in [24] as our starting point. The test bench includes a fully equipped quad-rotor that rotates without displacement; these platforms are very popular for eliminating risky flight testing during controller design [25,26,27]. Its flight controller can be programmed from tools developed ad hoc in MATLAB-Simulink. A simulation environment allows us to prove control structures on a model that faithfully reproduces the behavior of the real system in order to establish a benchmark on control engineering [28]. This has been the core of the competition promoted by the Control Engineering Group of the Spanish Automatic Committee (CEA) [29]. The control strategy presented in this paper was an entry in the competition; it achieved the best result and therefore was the winner of the competition [30].
Considering this background, the objectives of this paper are twofold. We aim, on the one hand, to present a methodological application of our model-based controller design approach. On the other hand, we aim to exemplify controller tuning by successive refinement and improvement based on the optimization of PID controllers. This method of optimization is based on linear model approximation in the first stage in order to obtain suitable initial conditions for the next stage, where the optimization is carried out by taking a non-linear dynamics model as the basis. The results are tested on the experimental setup, therefore showing the validity of the overall approach. Equations for the flight patterns of the UAVs and non-linear dynamic equations are provided, and experimental adjustments to the experimental platform are also conducted.
The document is organized as follows. The next section describes our UAV system with a non-linear model. Then, in Section 3, the general framework for the control system is stated. Section 4 describes the procedure for the design of our controller and also describes the multi-objective problem we aim to optimize. We describe our implementation and validation processes on a real UAV platform in Section 5, and the paper ends in Section 6 with some conclusions and future work.

2. System Framework

The experimental setup and non-linear model considered in this work have been taken from [24]. In this study, a fully equipped quadrotor is affixed to a structure that permits unrestricted rotation without translation. The flight controller programming and real-time duplex transmission of commands and flight states to the remote pilot are facilitated through a set of MATLAB-Simulink tools. Moreover, a representative simulator is provided for testing purposes, making the configuration highly suitable for orientation control. The experimental setup is illustrated in Figure 1. In this section, the fist step is the reviewing of the non-linear model, followed by an experiment to validate the model and compare it against data generated using the experimental setup.

2.1. Non-Linear Model of the UAV

As mentioned, the main objective is to control the orientation of a multirotor UAV, characterized by the Euler angles ϕ , θ , and ψ . To maintain a fixed position and prevent displacement, a frame and two concentric rings are used, rotating on axes oriented 90 , towards each other, allowing movement in ϕ and θ . Additionally, a mechanical locking system can be engaged to optionally restrict rotation in ψ . In Figure 2, the coordinate axes of the body system ( X B , Y B , and Z B ) and the terrestrial inertial reference system ( X G , Y G , and Z G ) are illustrated.
The system description is adapted from [24]. As expected, the propulsion forces align with the coordinates of the body system, leading to the multirotor UAV attaining speeds of p, q, and r along the X B , Y B , and Z B coordinates, respectively, following the expressions
p ˙ = M x ( I z I y ) q r I x q ˙ = M y ( I x I z ) p r I y r ˙ = M z ( I y I x ) p q I z
where the values I x , I y , and I z represent the inertial values with respect to the axis, and the variables M x , M y , and M z are the total moments about each axis, which are built as
M x = M Γ x + M r x + M μ x + M g x M y = M Γ y + M r y + M μ y + M g y M z = M δ z + M r z + M μ z + M g z
where M g is the moment due to the gravitational force and is represented as
M g x = m g d sin ( ϕ ) cos ( θ ) M g y = m g d sin ( θ ) M g z = 0
the m value denotes the mass of the multirotor UAV, d represents the distance from the center of gravity to the center of rotation { O } along the Z B axis, and g corresponds to the gravitational constant. At the same time, the values of M r represent the gyroscopic moments given by a rotating mass,
M r x = I r ( ω 1 ω 2 + ω 3 ω 4 ) q M r y = I r ( ω 1 ω 2 + ω 3 ω 4 ) p M r z = I r ( ω 1 ˙ ω 2 ˙ + ω 3 ˙ ω 4 ˙ )
where I r stands for the inertia of each rotor and the ω refers to their respective turning speeds, which are directly influenced by the control signal. On the other hand, M μ is equivalent to the moments due to static and viscous frictions in the mechanical couplings of the structure and is described as:
M μ x = μ C x p | p | μ x p M μ y = μ C y q | q | μ y q M μ z = μ C z r | r | μ z r
where the μ C variables are the Coulomb friction coefficients and the μ variables are the viscous friction coefficients. The coefficients are combined to yield the moment M Γ , which represents the total moment generated by the thrust forces of the propellers Γ . Additionally, the moment M δ corresponds to the rotation around Z B resulting from the algebraic combination of the drag moments δ . The mathematical expression describing these moments is shown as
M Γ x = l ( Γ 1 + Γ 2 + Γ 3 Γ 4 ) 2 2 M Γ y = l ( Γ 1 Γ 2 + Γ 3 + Γ 4 ) 2 2 M δ z = ( δ 1 + δ 2 δ 3 + δ 4 )
Table 1 shows the values of parameters in previous equations, which were obtained through measurements and experimentation. The variables Γ and δ , as well as the previously mentioned speeds ω , depend on the command signals of each motor in the form
Γ i = k Γ 2 ( M i 1000 ) 2 k i 2 10 6 + k Γ 1 ( M i 1000 ) k i 10 3 δ i = k δ 2 ( M i 1000 ) 2 k i 2 10 6 + k δ 1 ( M i 1000 ) k i 10 3 ω i = k ω 1 ( M i 1000 ) k i 10 3 + k ω 0
where M represents an RC PWM signal in between 1000 μ and 2000 μ . The coefficients k Γ , k δ , and k ω vary with the supply voltage to the motors, as shown in Table 2; the experimental setup allows the variable voltage to emulate the battery charge level during flight. The constants k 1 = k 2 = 1 , k 3 = 0.8 , and k 4 = 0.86 are due to the experimentally observed differences between rotors. Finally, the control signals of each motor are defined as
M 1 = u z u ϕ u θ u ψ M 2 = u z + u ϕ u θ + u ψ M 3 = u z + u ϕ + u θ u ψ M 4 = u z u ϕ + u θ + u ψ
where u ϕ , u θ , and u ψ represent the control actions for turning on each principal axis and u z represents the control action for maintaining the altitude in hypothetical free flight. Then, taking u z = 1500 µ leaves a symmetric range of ± 500 µ to control the spins up to the saturation of the motor command. A mechanical lock on the structure optionally prevents yaw rotation when, as in this case, only attitude control is being performed. Consequently, u ψ is forced to zero.
Having acquired this data, it becomes feasible to determine the velocities p, q, and r, concerning the control actions and the multirotor UAV position. The final step involves transforming the variables from the inertial reference system to the body system using quaternions. The rotational speed, expressed in quaternions in the reference system of the body corresponds to
q 1 ˙ = q 2 p q 3 q q 4 r 2 q 2 ˙ = q 1 p q 4 q + q 3 r 2 q 3 ˙ = q 4 p + q 1 q q 2 r 2 q 4 ˙ = q 3 p + q 2 q + q 1 r 2
from which we obtain q = [ q 1 q 2 q 3 q 4 ] T , a quaternion is employed to represent the orientation of the terrestrial reference system relative to the body’s reference system. In the given equations, it is essential for the quaternion q to be normalized in the form q = q / | q | so that the conjugate quaternion Q = [ Q 1 Q 2 Q 3 Q 4 ] T = [ q 1 q 2 q 3 q 4 ] T = q ¯ matches the inverse quaternion q 1 . Thus, the Euler angles are calculated as
θ = arctan 2 ( 2 Q 1 Q 2 + 2 Q 3 Q 4 , 1 2 Q 2 2 + 2 Q 3 2 ) ϕ = arcsin ( 2 Q 1 Q 3 2 Q 2 Q 4 ) ψ = arctan 2 ( 2 Q 1 Q 4 + 2 Q 2 Q 3 , 1 2 Q 3 2 + 2 Q 4 2 )

2.2. Simulation Model and Validation

The above equations are integrated into the UAV block of the simulation model shown in Figure 3. This block contains, in addition to the UAV dynamics, a mathematical model of the IMU (inertial measurement unit) integrated in the flight controller and the estimation system used to determine the UAV attitude. This provides a representative model of the real system. It should be noted that the IMU model incorporates digital low-pass filters with a cut-off frequency of 20 Hz. These filters can significantly affect the system’s behavior when the control bandwidth approaches the filter’s cutoff frequency.
The non-linear model of the UAV is integrated into a feedback control structure; the control algorithms, as the one described in this paper, are implemented within the control system block, and the set point and monitoring blocks complete the structure of Figure 3. Further details on this simulation model can be found in [24]. Figure 4 shows one validation experiment that compares the real system and simulated outputs. Ignoring high-frequency unmodeled dynamics and noise, both results match quite well.
The complexity of the control system to be implemented is limited by the low-cost flight controller in the real platform. Thus, the computational load is constrained by a sampling frequency of 50 Hz, the precision is bounded by an 8-bit microcontroller resolution, and the limited memory and processing power recommend light control laws. PID-type controllers are clear candidates for this situation. The relevance and main contribution of this paper will thus be their off-line synthesis and optimization, which must lead to adequate performance in the nonlinear simulator and, finally, in the real test bench.
The non-linear model and the simulator are used hereinafter to construct a linear model-based controller design. Once the non-linear model is validated, the simulator allows us to take an optimal tuning approach to the two-degree-of-freedom PID controller parameters. The optimization algorithm interacts with the simulator to test and evaluate the current controller performance. The next section presents our controller design based on this non-linear model. The resulting controller is finally implemented and tested on the real platform.

3. Control System

This section presents a description of the control system and describes how we approached the controller design process.
The considered control system is shown in Figure 5, whose variables are as below.
  • y ( s ) is the process-controlled variable.
  • r ( s ) is the set point for the process output.
  • u ( s ) is the controller’s output signal.
  • d ( s ) is the load disturbance of the system.
The transfer functions of the system are as follows. P ( s ) is the controlled process, and the controller is a two degrees of freedom (2DoF) PID controller, with parameters of p ¯ = K p , T i , T d , β T . The controller is implemented in terms of its reference and output aspects C y ( s ) and C r ( s ) as
C y ( s ) = K p + K i s + K d s 0.1 K d s + 1
C r ( s ) = β K p + K i s
where
u ( s ) = C r ( s ) r ( s ) C y ( s ) y ( s )
For the UAV system, there are two control loops, one for each variable of interest. It is possible to delineate a relationship between the above general control structures with respect to the specific case where ϕ and θ (the roll and pitch angles, respectively) are the controlled variables, u ϕ and u θ are the control signals, r ϕ and r θ are the reference angles, and d is any disturbance input that can cause the angle to deviate from its reference. In free flight, the main input disturbance is wind. Inside the structure that attaches the quadrotor, differences between the four rotors create important couplings between the roll and pitch control loops. These couplings are considered disturbance inputs for the control design.

Performance of the Control System

The effectiveness of a control system can be assessed by employing a cost function that relies on the error, defined as the difference between the desired value r (set point) and the current value of the controlled variable y (the system’s output). Certainly, when time increases, the error becomes larger and longer, and therefore, the system’s performance tends to degrade. As a measure of evaluating the controller’s performance, a widely used reference is a functional based on the integral of the error, like the integral absolute error ( I A E ) or the integral time-weighted absolute error ( I T A E ), which also considers the duration of the error over time.
There is another index that measures the effort of the control action, which is computed as the integral of the absolute value for the variation of the control signal ( I A V U ). For this index, when a higher value is obtained, it indicates that the control action varies more abruptly.
The formulation of the indices is stated as
I A E = 0 t f | r ( t ) y ( t ) | d t
I T A E = 0 t f t | r ( t ) y ( t ) | d t
I A V U = 0 t f d u ( t ) d t d t
where they can be measured for changes in the set point (servo control) or in the load disturbance (regulatory control).
The global evaluation of performance combines the previous indices as
J o = w 1 I A E + w 2 I T A E + w 3 I A V U
with the weighting factors w 1 = w 2 = 0.4 and w 3 = 0.2 , which evaluate the error indices ( I A E and I T A E ) and the control action ( I A U V ). These weighting factors are given by the benchmark from [24].
Then, the controller tuning design can be solved as an optimization problem that can be stated as
J o o = ˙ min K p , K i , K d , β J o
where J o o is the achieved optimal value for the index J o (17).

4. Controller Design

The design and tuning of the PID controller was carried out in three steps. First of all, a second-order linear model approximation for each channel (roll ( ϕ ) and pitch ( θ )) was identified using data obtained from the non-linear simulator. The control system’s design was established from those models using the internal model control (IMC) technique and then improved by optimization algorithms. Finally, the controllers were optimized using the non-linear simulation model.

4.1. Model Identification

The purpose of the identification of the obtained models was to help determine the initial controllers that were taken as initial conditions by the optimization algorithm.
Since the model identification was conducted through open-loop experiments, our first step was to determine the control actions that stabilized the UAV at roll ( ϕ ) and pitch ( θ ) angles of zero. The rotor asymmetries required the non-null actions of u ϕ = 9.2 µs and u θ = 47 µs.
Once a stable orientation was achieved, three tests were carried out in open-loop conditions in order to obtain the process responses to step inputs for the following conditions:
  • The response for ϕ to a step increment of 15 µs at u ϕ ;
  • The response for θ to a step increment of 15 µs at u θ .
Some measures were taken from each response, including the maximum value for the output ( y p ) and the time of its occurrence ( t p ), as well as the steady-state value for the system’s output ( y s s ). With all this information, the models were identified using the form
P ( s ) P m ( s ) = K ω n 2 s 2 + 2 ξ ω n s + ω n 2
where P m ( s ) is the identified model of P ( s ) , where K represents the static gain, ξ indicates the damping ratio, and ω n is the natural frequency.
From Model (19), we can see that it is well-known that the overshoot of the response ( M p ) can be calculated as
M p = y p y s s y s s
and using the following relationships, the model parameters are:
ξ = l n ( M p ) 2 π 2 + l n ( M p ) 2
ω n = π t p 1 ξ 2
The results of these performed identifications allow us to obtain the parameters of Equations (21) and (22), which are systematized in Table 3. The open-loop responses for each angle (for the real system and model) are shown in Figure 6.
Once the models were obtained, the controller was designed following the below steps:
  • The PID controller was designed using the internal model controller (IMC) approach in (23).
  • The PID-IMC’s tuning was improved by analyzing the multiobjective problem (18) on the basis of the identified linear models in Table 3.
  • A second refinement of the controller’s parameters was performed using the non-linear model presented in Section 2 as a simulation model for the optimization of Problem (18).
The final, third controller was tested on the real equipment. The following subsections will describe the detailed processes of each of the design steps.

4.2. Initial IMC Controller

An initial PID controller was designed from the identified second-order models using the IMC method. The expression for this controller is
C y ( s ) = 1 T c s 1 P m ( s )
where T c is the desired closed-loop time constant. Note that the resulting structure for Equation (23) is the same as for Equation (11).
Taking into consideration that the purpose of this IMC controller is to serve as the starting point to be improved by means of the optimization procedure in the next step, the tuning value of T c for each channel was selected in such a way that they provided the fastest closed-loop responses (in fact, T c is the desired closed-loop time constant) while preserving stability (or, in other words, without the system becoming unstable). The resulting values are T c = 0.16 for ϕ and T c = 0.06 for θ . This represents a closed-loop system bandwidth between 8 and 10 times that of the corresponding open loop.
The parameters of the tuned controllers C ϕ and C θ for each controlled variable, following Equation (13), are shown in Table 4. In this case, the set-point weighting factor was set to the default value of β = 1 .
In order to see the initial performance of the controllers using the IMC method at each channel (roll ( ϕ ) and pitch ( θ )), a test was conducted with controllers from Table 4. The results are shown in Figure 7, and the performance indices I A E , I T A E , and I A U V , as well as the global index performance J o , are displayed in Table 5. They were computed by applying Equation (14) to (17).
As can be seen, even the reference is followed with coarse precision. Thus, it is possible to take this controller (IMC) as initial conditions to be improved as one of the main purposes of the paper.

4.3. Improved Optimized Controller

The optimal controller’s parameters were determined using the MATLAB function fminsearch. The initial controllers were obtained through the IMC method (as shown in Table 4) and the optimization process is using the linear models provided in Table 3.
In this manner, the optimization process runs smoother and faster. The function to be optimized is J o (17). The parameters obtained through the optimization problem in Equation (18) are shown in Table 6.
The responses for the controlled-variable roll ( ϕ ) and pitch ( θ ) using the optimized PID controllers shown in Table 6 are shown in Figure 8, and our evaluation of their performance is shown in Table 7.
The proposed control strategy not only performs well in tracking reference r changes but also in correcting any other cause of output deviations, i.e., compensating for disturbance inputs. In the present case, the experimentally observed differences between the rotors generate couplings between the roll and pitch control loops, which act as disturbance inputs. These couplings have more effect on the axis with less inertia (Table 1 concludes I x I y ). Consequently, Figure 8 shows how the roll ( ϕ ) angle ( Y B -axis rotation) hardly deviates when the pitch reference ( r θ ) changes, while the dynamic recovery of the pitch ( θ ) angle ( X B -axis rotation) is fast but visible after a change in the roll reference ( r ϕ ).
In addition, it is possible to see that the optimized controllers achieved from the linear model provide a significant improvement in the performance evaluation by comparing the data from Table 5 and Table 7.

4.4. Optimized Controller for the Non-Linear Model

In this step, the tuning of the controller consisted of the optimization of Equation (17) by means of Equation (18) and using the MATLAB function fminsearch. However, we also considered the non-linear model (NL) plant described in Section 2 in this process. The parameters in Table 6 were used as initial controllers for optimization.
The controller parameters achieved from the above for the current optimization are shown in Table 8.
The control structure of Figure 5 and controller parameters of Table 8 were implemented inside the control system block of the non-linear simulator in Figure 3. Consequently, Figure 9 shows the reference tracking and coupling compensation capabilities for the roll ( ϕ ) and pitch ( θ ) angles. Additionally, Table 9 shows the values of the performance indices for this case.
It is important to highlight that conducting the optimization directly with the non-linear model of the system, rather than using the identified linear model, resulted in a substantial improvement in performance. This can be checked by comparing the values in Table 7 and Table 9.

5. Implementation and Validation Results

The final controllers obtained in Section 4.4 were implemented in the real platform of the UAV system. This was a part of a competition promoted by the Control Engineering Group of the Spanish Automatic Committee (CEA) [29]. The three finalist designs are labeled T1, T2 and T3, with the latter being the one proposed in this work.
As a validation test, the controllers were evaluated using different power supply voltages in the UAV system. Experiments were conducted for battery voltages of 10 V, 11 V, and 12 V. This was performed with the aim of verifying how the state of the battery charge affects the behavior of the UAV when different supply voltages are used.
Table 10 shows the value for the performance indices for the roll ( ϕ ) and pitch ( θ ) angles for each of the power supply voltages. Moreover, the control system responses for the case of 11V are shown in Figure 10, where, in addition to the signal for each angle ( ϕ and θ ), the speeds reached by the multi-rotor UAVs (p and q) are shown, as well as the control signals ( u ϕ and u θ ). A video of this experiment in the real platform can be seen in [30].
In order to have more usable information, it is possible to calculate an averages for the experiments for each control design from the data of Table 10, as follows:
J o ¯ = J o ϕ + J o θ 2
Moreover, a general global index J O ¯ can be computed as an average of each of the tests. All these calculations are provided by Table 11, where it is possible to see that the proposed multi-objective PID design is the best and gives the lowest value for J O ¯ .

6. Conclusions

This study focuses on addressing the pitch–roll angle orientation in a quadrotor UAV drone system. The control approach involves the use of PID controllers and tuning them through multi-objective optimization. The UAV system is represented by a non-linear model, wit its respective interactions. The experimental setup successfully validates the non-linear model, as it aligns with the conducted control experiment.
The proposed control design for the PID controllers shows good performance for the system, providing low values for the indices I A E , I T A E , I A U V , and the global J o and improving the results of the multi-objective optimization problem.
Three different control designs were implemented in the real UAV platform. Ultimately, the control method proposed in this paper offered the best performance for the system compared to the other two methods.
Our future work will involve the inclusion of a third variable to control, the angle ψ , which adds more complexity to the problem since the drag forces acting on the rotors are much smaller than their thrust forces. Other improvements will include a robustness analysis to take into account the battery discharge, which reduces rotor force. Finally, it will be beneficial to extend our validations to a free-flight experiment from the confinement platform once altitude control is included.

Author Contributions

Conceptualization, O.A., J.D.R. and R.V.; data curation, D.C., J.D.R. and J.R.-A.; formal analysis, O.A., J.D.R., M.G.-M. and R.V.; funding acquisition, O.A. and R.V.; investigation, O.A., D.C., J.R.-A. and R.V.; methodology, D.C., J.D.R. and R.V.; software, D.C. and J.R.-A.; project administration, O.A.; resources, O.A., J.D.R. and M.G.-M.; supervision, M.G.-M. and R.V.; validation, J.D.R., M.G.-M. and R.V.; visualization, J.R.-A.; writing—original draft, O.A.; writing—review and editing, M.G.-M. and R.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work has received support from the Catalan Government under Project 2021 SGR 00197 and from the Spanish Government under the MICINN projects PID2019-105434RB-C33 and TED2021-806 129134B-I00, co-funded with European Union ERDF funds. The financial support provided by the University of Costa Rica under the grant 731-B9-265 is also greatly appreciated. Finally, this work has also been partially supported by the University of La Rioja, through the grant REGI 22/40.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mohsan, S.A.H.; Khan, M.A.; Noor, F.; Ullah, I.; Alsharif, M.H. Towards the Unmanned Aerial Vehicles (UAVs): A Comprehensive Review. Drones 2022, 6, 147. [Google Scholar] [CrossRef]
  2. Rejeb, A.; Abdollahi, A.; Rejeb, K.; Treiblmaier, H. Drones in agriculture: A review and bibliometric analysis. Comput. Electron. Agric. 2022, 198, 107017. [Google Scholar] [CrossRef]
  3. Ecke, S.; Dempewolf, J.; Frey, J.; Schwaller, A.; Endres, E.; Klemmt, H.J.; Tiede, D.; Seifert, T. UAV-Based Forest Health Monitoring: A Systematic Review. Remote Sens. 2022, 14, 3205. [Google Scholar] [CrossRef]
  4. Shakhatreh, H.; Sawalmeh, A.H.; Al-Fuqaha, A.; Dou, Z.; Almaita, E.; Khalil, I.; Othman, N.S.; Khreishah, A.; Guizani, M. Unmanned Aerial Vehicles (UAVs): A Survey on Civil Applications and Key Research Challenges. IEEE Access 2019, 7, 48572–48634. [Google Scholar] [CrossRef]
  5. Madridano, A.; Campos, S.; Al-Kaff, A.; García, F.; Martín, D.; Escalera, A. Unmanned aerial vehicle for fire surveillance and monitoring. Rev. Iberoam. Autom. Inform. Ind. 2020, 17, 254–263. [Google Scholar] [CrossRef]
  6. Kangunde, V.; Jamisola, R.S.; Theophilus, E.K. A review on drones controlled in real-time. Int. J. Dyn. Control 2021, 9, 1832–1846. [Google Scholar] [CrossRef]
  7. Bigazzi, L.; Gherardini, S.; Innocenti, G.; Basso, M. Development of Non Expensive Technologies for Precise Maneuvering of Completely Autonomous Unmanned Aerial Vehicles. Sensors 2021, 21, 391. [Google Scholar] [CrossRef]
  8. Rico, R.; Rico-Azagra, J.; Gil-Martínez, M. Hardware and RTOS Design of a Flight Controller for Professional Applications. IEEE Access 2022, 10, 134870–134883. [Google Scholar] [CrossRef]
  9. Khan, S.; Jaffery, M.H.; Hanif, A.; Asif, M.R. Teaching Tool for a Control Systems Laboratory Using a Quadrotor as a Plant in MATLAB. IEEE Trans. Educ. 2017, 60, 249–256. [Google Scholar] [CrossRef]
  10. Gonzalez-Vargas, A.M.; Serna-Ramirez, J.M.; Fory-Aguirre, C.; Ojeda-Misses, A.; Cardona-Ordoñez, J.M.; Tomba-Andrade, J.; Soria-Lopez, A. A low-cost, free-software platform with hard real-time performance for control engineering education. Comput. Appl. Eng. Educ. 2019, 27, 406–418. [Google Scholar] [CrossRef]
  11. Zhang, X.; Li, X.; Wang, K.; Lu, Y. A Survey of Modelling and Identification of Quadrotor Robot. Abstr. Appl. Anal. 2014, 2014, 320526. [Google Scholar] [CrossRef]
  12. Shraim, H.; Awada, A.; Youness, R. A survey on quadrotors: Configurations, modeling and identification, control, collision avoidance, fault diagnosis and tolerant control. IEEE Aerosp. Electron. Syst. Mag. 2018, 33, 14–33. [Google Scholar] [CrossRef]
  13. Nascimento, T.P.; Saska, M. Position and attitude control of multi-rotor aerial vehicles: A survey. Annu. Rev. Control. 2019, 48, 129–146. [Google Scholar] [CrossRef]
  14. Marshall, J.A.; Sun, W.; L’Afflitto, A. A survey of guidance, navigation, and control systems for autonomous multi-rotor small unmanned aerial systems. Annu. Rev. Control. 2021, 52, 390–427. [Google Scholar] [CrossRef]
  15. Gil-Martinez, M.; Rico-Azagra, J. Multi-rotor robust trajectory tracking. In Proceedings of the 2015 23rd Mediterranean Conference on Control and Automation (MED), Torremolinos, Spain, 16–19 June 2015; pp. 865–870. [Google Scholar] [CrossRef]
  16. Rubí, B.; Pérez, R.; Morcego, B. A Survey of Path Following Control Strategies for UAVs Focused on Quadrotors. J. Intell. Robot. Syst. 2020, 98, 241–265. [Google Scholar] [CrossRef] [Green Version]
  17. Lotufo, M.A.; Colangelo, L.; Perez-Montenegro, C.; Canuto, E.; Novara, C. UAV quadrotor attitude control: An ADRC-EMC combined approach. Control. Eng. Pract. 2019, 84, 13–22. [Google Scholar] [CrossRef]
  18. Chen, Y.; Zhang, G.; Zhuang, Y.; Hu, H. Autonomous Flight Control for Multi-Rotor UAVs Flying at Low Altitude. IEEE Access 2019, 7, 42614–42625. [Google Scholar] [CrossRef]
  19. Invernizzi, D.; Panizza, P.; Riccardi, F.; Formentin, S.; Lovera, M. Data-driven attitude control law of a variable-pitch quadrotor: A comparison study. IFAC-PapersOnLine 2016, 49, 236–241. [Google Scholar] [CrossRef]
  20. Capocchiano, S.; Panizza, P.; Invernizzi, D.; Lovera, M. Closed-loop data-driven attitude control design for a multirotor UAV. In Proceedings of the IEEE Conference on Control Technology and Applications (CCTA), Copenhagen, Denmark, 21–24 August 2018; pp. 153–158. [Google Scholar]
  21. Zangarini, A.; Invernizzi, D.; Panizza, P.; Lovera, M. Closed-loop MIMO data-driven attitude control design for a multirotor UAV. CEAS Aeronaut. J. 2020, 11, 873–884. [Google Scholar] [CrossRef]
  22. Bo, G.; Xin, L.; Hui, Z.; Ling, W. Quadrotor helicopter attitude control using cascade PID. In Proceedings of the 28th Chinese Control and Decision Conference, CCDC 2016, Yinchuan, China, 28–30 May 2016; pp. 5158–5163. [Google Scholar] [CrossRef]
  23. Wang, P.; Man, Z.; Cao, Z.; Zheng, J.; Zhao, Y. Dynamics modelling and linear control of quadcopter. In Proceedings of the International Conference on Advanced Mechatronic Systems, ICAMechS, Melbourne, Australia, 30 November–3 December 2016; pp. 498–503. [Google Scholar] [CrossRef]
  24. Rico-Azagra, J.; Gil-Martínez, M.; Rico, R.; Nájera, S.; Elvira, C. A benchmark for orientation control of a multirotor in a three degrees-of-freedom rotation structure. Rev. Iberoam. Autom. Inform. Ind. 2021, 18, 265–276. [Google Scholar] [CrossRef]
  25. Hancer, M.; Bitirgen, R.; Bayezit, I. Designing 3-DOF Hardware-In-The-Loop Test Platform Controlling Multirotor Vehicles. IFAC-PapersOnLine 2018, 51, 119–124. [Google Scholar] [CrossRef]
  26. Bajelani, M.; Tayefi, M.; Zhu, M. A real-test and simulation combined platform for developing intelligent tracking control of multirotors. Int. J. Intell. Unmanned Syst. 2022, in press. [CrossRef]
  27. Chávez-Gudiño, M.A.; Concha-Sánchez, A.; Maciel-Barboza, F.M.; Gadi, S.K.; Thenozhi, S.; Jiménez Betancourt, R. Development and control of a low cost 2 DOF laboratory helicopter. RIAI Rev. Iberoam. Autom. Inform. Ind. 2023, in press. [Google Scholar] [CrossRef]
  28. Ebeid, E.; Skriver, M.; Terkildsen, K.; Jensen, K.; Schultz, U. A survey of open-source UAV flight controllers and flight simulators. Microprocess. Microsyst. 2018, 61, 11–20. [Google Scholar] [CrossRef]
  29. Gallarta-Sáenz, D.; Rico-Azagra, J.; Gil-Martínez, M. Learning Enhancement of Control Engineering: A Competition-Based Case. IEEE Access 2023, 11, 38240–38250. [Google Scholar] [CrossRef]
  30. Control Engineering—Spanish Automatic Committee (CEA) and University of La Rioja Research Group. Competition: Orientation Control of a Multirotor—CIC2021. 2021. Available online: https://www.unirioja.es/dptos/die/cic2021/resultados.shtml (accessed on 16 July 2023).
Figure 1. Experimental UAV setup.
Figure 1. Experimental UAV setup.
Mathematics 11 03390 g001
Figure 2. Reference systems and main variables.
Figure 2. Reference systems and main variables.
Mathematics 11 03390 g002
Figure 3. Simulation model developed in MATLAB-Simulink.
Figure 3. Simulation model developed in MATLAB-Simulink.
Mathematics 11 03390 g003
Figure 4. Non-linear model validation with the experimental responses.
Figure 4. Non-linear model validation with the experimental responses.
Mathematics 11 03390 g004
Figure 5. The considered control system structure.
Figure 5. The considered control system structure.
Mathematics 11 03390 g005
Figure 6. Real system and model responses for roll and pitch to a 15 µs step in their control actions.
Figure 6. Real system and model responses for roll and pitch to a 15 µs step in their control actions.
Mathematics 11 03390 g006
Figure 7. Roll and pitch responses in reference tracking using the IMC controller.
Figure 7. Roll and pitch responses in reference tracking using the IMC controller.
Mathematics 11 03390 g007
Figure 8. Roll and pitch responses in reference tracking using the optimized (L) controller.
Figure 8. Roll and pitch responses in reference tracking using the optimized (L) controller.
Mathematics 11 03390 g008
Figure 9. Simulated non-linear (NL) responses of the final optimized controller for roll ( ϕ ) and pitch ( θ ) rotations.
Figure 9. Simulated non-linear (NL) responses of the final optimized controller for roll ( ϕ ) and pitch ( θ ) rotations.
Mathematics 11 03390 g009
Figure 10. UAV control system responses for roll and pitch rotations with 11 V battery charge.
Figure 10. UAV control system responses for roll and pitch rotations with 11 V battery charge.
Mathematics 11 03390 g010
Table 1. Rotary body parameters.
Table 1. Rotary body parameters.
Value Value
l0.25 (m) μ C x 0.006 (Nm/rad)
d0.015 (m) μ C y 0.006 (Nm/rad)
m0.375 (kg) μ C z 0.004 (Nm/rad)
I x 0.00135 (kg2) μ x 0.0015 (Nms/rad)
I y 0.019 (kg2) μ y 0.0015 (Nms/rad)
I z 0.04 (kg2) μ z 0.05 (Nms/rad)
I r 1.05 × 10 6 (kg2)
Table 2. Propulsion system parameters.
Table 2. Propulsion system parameters.
9 V10 V11 V12 V
k Γ 2 1.77961.20822.01052.0593
k Γ 1 0.95541.11791.57071.9021
k δ 2 0.02380.02390.03180.0326
k δ 1 0.02020.02270.02660.0298
k ω 2 1445.81592.41762.91921.5
k ω 1 63.47270.50575.90784.729
Table 3. Parameters for the second-order model P m ( s ) .
Table 3. Parameters for the second-order model P m ( s ) .
P ϕ P θ
K14.40014.700
ω n 3.0002.000
ξ 0.4520.316
Table 4. Controller parameters for the IMC method.
Table 4. Controller parameters for the IMC method.
IMC
C ϕ C θ
K p 0.0390.358
K i 0.1291.134
K d 0.0140.283
β 11
Table 5. Performance indices for the IMC controllers.
Table 5. Performance indices for the IMC controllers.
IAEITAEIAVU J o
Roll ( ϕ )2.9963.3650.0592.556
Pitch ( θ )2.4373.9670.3332.628
Table 6. Controller parameters for optimization with the linear model (L).
Table 6. Controller parameters for optimization with the linear model (L).
Optimal (L)
C ϕ C θ
K p 0.1290.803
K i 0.0040.403
K d 0.0240.159
β 1.5190.984
Table 7. Performance indices for the optimized (L) controllers.
Table 7. Performance indices for the optimized (L) controllers.
IAEITAEIAVU J o
Roll ( ϕ )1.0630.7810.1740.772
Pitch ( θ )0.8940.5250.4100.650
Table 8. Controller parameters for optimization with the non-linear model (NL).
Table 8. Controller parameters for optimization with the non-linear model (NL).
Optimal (NL)
C ϕ C θ
K p 0.1370.936
K i 0.1490.803
K d 0.0280.234
β 1.2730.826
Table 9. Performance indices for the optimized (NL) controllers.
Table 9. Performance indices for the optimized (NL) controllers.
IAEITAEIAVU J o
Roll ( ϕ )0.9580.6090.2220.671
Pitch ( θ )0.6460.3170.4980.484
Table 10. Performance evaluation and comparison for the UAV system.
Table 10. Performance evaluation and comparison for the UAV system.
Control
Design
[V]Roll( ϕ )
IAEITAEIAVU J o ϕ
101.2681.0000.2000.947
T1111.3091.2610.2451.077
121.5912.0820.3891.547
101.5232.3270.0771.556
T2111.5261.7540.0921.327
121.4981.8970.1291.384
101.4051.2900.1791.114
T3111.4081.6470.2321.268
121.3842.2330.3801.523
Control
Design
[V]Pitch ( θ )
IAEITAEIAVU J o ϕ
101.2801.9370.3041.348
T1111.0941.3800.3091.056
120.9731.0490.5030.905
101.7402.2640.1181.625
T2111.6702.5650.1591.726
121.5111.9430.1851.418
100.8670.7350.4560.732
T3110.8180.6770.5860.715
120.7800.6410.7590.720
Table 11. Average values for the performance index using the UAV system.
Table 11. Average values for the performance index using the UAV system.
Control
Design
[V] J o ¯ = J o ϕ + J o θ 2 J O ¯
101.148
T1111.0661.147
121.228
101.590
T2111.5261.506
121.401
100.923
T3110.9921.012
121.121
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Arrieta, O.; Campos, D.; Rico-Azagra, J.; Gil-Martínez, M.; Rojas, J.D.; Vilanova, R. Model-Based Optimization Approach for PID Control of Pitch–Roll UAV Orientation. Mathematics 2023, 11, 3390. https://doi.org/10.3390/math11153390

AMA Style

Arrieta O, Campos D, Rico-Azagra J, Gil-Martínez M, Rojas JD, Vilanova R. Model-Based Optimization Approach for PID Control of Pitch–Roll UAV Orientation. Mathematics. 2023; 11(15):3390. https://doi.org/10.3390/math11153390

Chicago/Turabian Style

Arrieta, Orlando, Daniel Campos, Javier Rico-Azagra, Montserrat Gil-Martínez, José D. Rojas, and Ramon Vilanova. 2023. "Model-Based Optimization Approach for PID Control of Pitch–Roll UAV Orientation" Mathematics 11, no. 15: 3390. https://doi.org/10.3390/math11153390

APA Style

Arrieta, O., Campos, D., Rico-Azagra, J., Gil-Martínez, M., Rojas, J. D., & Vilanova, R. (2023). Model-Based Optimization Approach for PID Control of Pitch–Roll UAV Orientation. Mathematics, 11(15), 3390. https://doi.org/10.3390/math11153390

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