Simulation and flight experiments of a quadrotor tail-sitter vertical take-off and landing unmanned aerial vehicle with wide flight envelope

Abstract

This paper presents the modeling, simulation, and control of a small-scale electric powered quadrotor tail-sitter vertical take-off and landing unmanned aerial vehicle. In the modeling part, a full attitude wind tunnel test is performed on the full-scale unmanned aerial vehicle to capture its aerodynamics over the flight envelope. To accurately capture the degradation of motor thrust and torque at the presence of the forward speed, a wind tunnel test on the motor and propeller is also carried out. The extensive wind tunnel tests, when combined with the unmanned aerial vehicle kinematics model, dynamics model and other practical constraints such as motor saturation and delay, lead to a complete flight simulator that can accurately reveal the actual aircraft dynamics as verified by actual flight experiments. Based on the developed model, a unified attitude controller and a stable transition controller are designed and verified. Both simulation and experiments show that the developed attitude controller can stabilize the unmanned aerial vehicle attitude over the entire flight envelope and the transition controller can successfully transit the unmanned aerial vehicle from vertical flight to level flight with negligible altitude dropping, a common and fundamental challenge for tail-sitter vertical take-off and landing aircrafts. Finally, when supplied with the designed controller, the tail-sitter unmanned aerial vehicle can achieve a wide flight speed envelope ranging from stationary hovering to fast level flight. This feature dramatically distinguishes our aircraft from conventional fixed-wing airplanes.

Introduction

Vertical take-off and landing (VTOL) vehicles which possess VTOL ability, high maneuverability during vertical flight, and high efficiency in level flight have been investigated for a long time.^1–3 Especially for the manned Vertical/Short Take-Off and Landing (V/STOL), many attempts have become mature and been serving for the military until now such as the Harrier GR7, V-22 Osprey, and F-35 Lightning. Small-scale electric powered VTOL unmanned aerial vehicles (UAVs), although in the infancy phase, also have attracted lots of recent attention. ⁴ They have great potential in many industrial applications: aerial photography, precise agriculture, intelligent surveying and monitoring, parcel delivery, border patrol, etc. Comparing with the widely used rotary-wing UAVs which have high maneuverability but low power efficiency and fixed-wing UAVs which have high cruise power efficiency but need extra facilities (e.g. runway, catapult) for take-off and landing, a VTOL UAV combining the advantages of rotary-wing vehicles and fixed-wing vehicles is a promising category of UAV which provides greater flexibility and capability (e.g. flight range, endurance, maneuverability) in many of the aforementioned applications.

There are mainly four categories of VTOL UAVs^4,5: dual-system (e.g. a quadplane), tilt-rotor, tilt-wing, ⁶ and tail-sitter. ⁷ Among them, the tail-sitter VTOL aircrafts are perhaps the simplest implementation as they do not require extra transition actuators. These are useful in saving weight and reducing manufacturing complexity especially for lightweight UAVs with small size. Currently, lots of research have been done on the small-size tail-sitter VTOL UAV. Bapst et al. ⁸ proposed a dual rotor tail-sitter VTOL which consists of a flying wing with elevons and two rotors. An airflow model is proposed to estimate the reference airspeed (i.e. the combination of vehicle airspeed and propeller-induced airspeed) according to the actuator power consumption. A velocity controller is designed to calculate the desired attitude and thrust. Due to the poor measurement of the airspeed sensor on a large angle of attack, this controller is only verified with a small angle of attack. Verling et al. ⁹ proposed a full attitude controller. A half-model wind tunnel test is conducted to collect the aerodynamic force and moment at different angles of attack, elevon deflections, rotor speeds, and airspeeds. The collected data are used to approximate an aerodynamic model proposed by Anderson. ¹⁰ Then, the approximated aerodynamic model is used in the attitude controller. The outdoor experiment shows that this designed controller can work in both vertical flight and level flight. Ritz and D’Andrea ¹¹ implemented a global controller for a dual rotor tail-sitter. A parameter learning scheme is used to estimate the pitching moment and aerodynamic lift and drag force of the vehicle. They use the on-board sensor to estimate the aerodynamics in real time. However, this may lead to inaccurate estimation if the vehicle enters a new flight regime. Similarly, Bronz et al. ¹² developed a dual rotor tail-sitter named “cyclone". The aerodynamic model comes from Jameson ¹³ and Bronz and Drouin. ¹⁴ Since the accurate knowledge of airspeed and angle of attack is difficult to obtain. They apply the incremental nonlinear dynamic inversion controller to the vehicle, which needs no modeling of the vehicle dynamics and is very strong at disturbance rejection. The flight test shows that the developed vehicle can accomplish hovering flight and the transition flight between vertical flight and level flight.

According to Saeed, et al., ⁴ the tail-sitter can be roughly divided into two different categories: control surface transitioning tail-sitter (CSTT) and differential thrust transitioning tail-sitter (DTTT). The research work mentioned above belongs to CSTT. Comparing with a CSTT, a DTTT can easily produce a rapid and strong control moment during vertical flight, which makes it easy to control and has a better anti-wind performance during the vertical flight. Sinha et al. ¹⁵ proposed a DTTT UAV named “Quadshot”u The control system has two modes: hover mode and forward flight mode. The forward flight mode can accomplish the forward transition by smoothly varying the pitch command. No aerodynamics is considered in the control system. Oosedo et al. ¹⁶ developed an asterisk-type quadrotor tail-sitter VTOL UAV, which only use the differential thrust from the four rotors to produce the control force and moment. An optimal transition strategy ¹⁷ is developed to accomplish the forward transition (no backward transition results are reported). In the modeling process, the authors only consider the aerodynamics in the longitudinal direction. Theys et al. ¹⁸ developed a velocity controller based on the full envelop aerodynamics of the vehicle. The authors only use NACA 0012 airfoil aerodynamic coefficients to model the aircraft aerodynamics, which will lead to modeling error. The flight results show that the modeling error degrades the controller performance.

For a DTTT tail-sitter that is concerned in this paper (see Figure 1 and Table 1), it needs to perform transition maneuver between vertical and level flight by varying its attitude (i.e. pitch angle). During this process, the velocity and angle of attack change dramatically, leading to complicated aerodynamics on the vehicle. The wide flight envelope presents a great challenge in modeling and controlling of such UAVs. To handle this problem, building a complete model which can accurately predict the vehicleel motion is crucially important. Unlike most of the quadrotor modeling process, which is regarded as a simple rigid body with a motor collective thrust and differential torque that solely depend on the input pulse-width modulation (PWM) signal,^19,20 the tail-sitter poses complicated aerodynamic forces and moments that dramatically deviate the UAV motion from a quadrotor. In addition, when the tail-sitter is in the level flight, a significant drop in propeller thrust and torque is encountered due to the presence of forward speed.

Figure 1.

The quadrotor tail-sitter with frames defined in the Modeling section and specifications shown in Table 1.

Table 1.

The tail-sitter VTOL UAV specification.

Component	Specification
Wing span	1.01 m
Surface area	0.24 m2
Gross weight	1.4 kg
Auto-pilot	Pixhawk 4 autopilot
Actuator	DJI E305
Battery	4 cell, 4480 mAh

In our previous work, ²¹ we have presented a complete model of a quadrotor tail-sitter, including full envelope aerodynamic model, propeller aerodynamics, motor dynamics, UAV kinematics and dynamics with all minor effects such as gyroscopic effects, motor acceleration/deceleration effects, and implemented the corresponding simulation. In this paper, wind tunnel tests are performed on the full-scale vehicle and the propulsion system, respectively, to capture the vehicle three-dimensional (3D) aerodynamic parameters and propeller aerodynamic coefficients to replace the original one (the corresponding test results are submitted along with this paper). The comparison of simulation results and experimental results show that the developed model is very accurate. A unified attitude controller and a stable transition controller are designed and the experiment results show that the developed attitude controller can stabilize the UAV attitude over the entire flight envelope which enables the vehicle to fly at a wide flight speed envelope ranging from stationary hovering to fast level flight, and the transition controller can successfully transit the UAV between vertical flight and level flight with negligible altitude dropping/gaining.

Modeling

Coordinate systems

Four coordinate systems are used in this paper. As shown in Figure 1, they are inertial frame, body frame, stability frame, and rotor frame. The origin of the inertial frame (shown as O_i) is set at the take-off point with x_i, y_i, and z_i, respectively, being aligned with North, East, and Downside. The body frame (denoted by x_b, y_b, and z_b) is defined to be the same as that of conventional fixed-wing aircrafts. To ease the aerodynamics analysis, the stability frame is commonly used. The stability frame is defined by rotating the body frame along y_b by $-\alpha$ such that the x axis of the stability frame is parallel to the longitudinal airspeed. ²² In addition to the inertial frame, body frame, and stability frame, we define a rotor frame for each of the four rotors to account for the fact that the rotors are usually inclined along the motor diagonal line by angle β_m to improve the roll (along x_b in Figure 1) moment during vertical flight. ²³ The origin of the rotor frame is set at the rotor (rotor here represents the rotation parts of the motor and the attached propeller) gravity center. $x_m^{\left(i\right)}$ , denoting the x axis of the ith motor frame, coincides with the ith motor rotation axis; $y_m^{\left(i\right)}$ , denoting the y axis of the ith motor, is aligned with the motor diagonal line.

Kinematics

The kinematics of the vehicle can be divided into two parts: translation and rotation. The translation part describes the displacement of the center of gravity (CG) of the aircraft which can be written as

$$\dot{\mathit{p}}=\mathit{v}$$

(1)where the p and v are, respectively, the position and velocity described in the inertial frame. The rotation part describes the attitude of the vehicle. There are three methods to parameterize the vehicleer attitude: rotation matrix, three-parameters method, and quaternion. The attitude kinematics of the vehicle using rotation matrix $\mathit{R}\in {\mathbb{ℝ}}^{3\times 3}$ is written as

$$\dot{\mathit{R}}=\mathit{R}\widehat{\mathit{\omega }}$$

(2)where the hat map $\widehat{\dot{}}$ is an operation that converts body angular velocity $\mathit{\omega }={\left[\begin{array}{ccc}{\omega }_x& {\omega }_y& {\omega }_z\end{array}\right]}^T\in {\mathbb{ℝ}}^3$ to a skew symmetric matrix which can be written as

$$\widehat{\mathit{\omega }}=\left[\begin{array}{ccc}0& -{\omega }_z& {\omega }_y\\ {\omega }_z& 0& -{\omega }_x\\ -{\omega }_y& {\omega }_x& 0\end{array}\right]$$

(3)

According to Schaub and Junkins, ²⁴ a generally rigid body only needs three independent parameters to parameterize the attitude. However, the rotation matrix has nine parameters which are highly redundant and will increase the computing complexity. To reduce such redundancy, a number of three-parameters methods have been developed, such as Euler angles or Tait-Bryan angles through three successive rotations, parameterization methods developed by Tsiotras and Longuski ²⁵ through two successive rotations, and angle-axis method through one single rotation. Unfortunately, it has been proven that none of the three-parameters methods can be singularity free.^24,26 A good balance between the singularity free and parameter redundancy is the quaternion-based parameterization method, which is singularity free while keeps only one redundant parameter. The attitude propagation using quaternion representation is written as ²²

$$\dot{\mathit{q}}=\frac{1}{2}\left[\begin{array}{cc}0& -{\mathit{\omega }}^T\\ \mathit{\omega }& -\widehat{\mathit{\omega }}\end{array}\right]\mathit{q}$$

(4)where $q={\left[\begin{array}{cc}\eta & {ϵ}^T\end{array}\right]}^T\in {\mathbb{ℝ}}^4$ , $\eta \in \mathbb{ℝ}$ and $\mathit{ϵ}\in {\mathbb{ℝ}}^3$ are, respectively, the scalar part and the vector part of the quaternion q . In our simulation, the attitude is updated using equation (4) and converted to rotation matrix for further usage through

$$\mathit{R}=\mathit{I}+2\eta \widehat{ϵ}+2{\widehat{ϵ}}^2$$

(5)

The quaternion is ideal to update the attitude in the simulation as it is singularity free and has few redundant parameters. However, it is hard for visualization and human interpretation. In conventional aircraft analysis, ZYX Tait-Brant angles are mostly used. ²² It is human intuitive but has singularity at pitch angle of ±90°. For the tail-sitter UAV, the vehicle is designed to transit from level to vertical flight by tilting its pitch angle nearly 90°, resulting in numerical instability for ZYX Tait-Brant angles. To avoid such singularity, ²¹ we choose ZXY Tait-Brant angles to represent the attitude that are used to monitor the aircraft attitude or accept human pilottd command. With the ZXY Tait-Brant angles, the singularity point is successfully shifted from $\theta =\pm {90}^{°}$ to $\varphi =\pm {90}^{°}$ which is seldom reached for a tail-sitter UAV. The ZXY Tait-Brant angles can be computed from the rotation matrix by using the following formulas

$$\begin{array}{l}\theta =\tan -1\left(\frac{-r_{31}}{r_{33}}\right)\\ \varphi =\sin -1\left(r_{32}\right)\\ \psi =\tan -1\left(\frac{-r_{12}}{r_{22}}\right)\end{array}$$

(6)where $r_{ij}\left(i,j=1,2,3\right)$ are the entries of the rotation matrix R .

Rigid body dynamics

We consider the vehicle to be a rigid body when building its dynamic model. Using the Newton–Euler equation, we can derive the rotation and translation dynamic of the vehicle as follows

$$m\dot{\mathit{v}}={\mathit{F}}_t$$

(7)

$${\mathit{J}}_b\dot{\mathit{\omega }}+\mathit{\omega }\times \left({\mathit{J}}_b\mathit{\omega }\right)={\mathit{M}}_t$$

(8)where m is the mass of the vehicle; ${\mathit{J}}_b$ is the inertial matrix of the vehicle; ${\mathit{F}}_t$ and ${\mathit{M}}_t$ are, respectively, the total force and moment applied on the vehicle, which can be, respectively, written as

$$\begin{array}{l}{\mathit{F}}_t=\mathit{R}\left({\displaystyle \sum {\mathit{R}}_m^{\left(i\right)}\left[\begin{array}{c}{T}_i\\ 0\\ 0\end{array}\right]+{\mathit{f}}_a}+\left[\begin{array}{c}0\\ 0\\ mg\end{array}\right]\right)\\ {\mathit{M}}_t={\displaystyle \sum {\widehat{r}}_i{\mathit{R}}_m^{\left(i\right)}\left[\begin{array}{c}{T}_i\\ 0\\ 0\end{array}\right]}+{\displaystyle \sum {\mathit{R}}_m^{\left(i\right)}}\left[\begin{array}{c}-s_i{Q}_i\\ 0\\ 0\end{array}\right]\\ +{\displaystyle \sum {\mathit{R}}_m^{\left(i\right)}}\left[\begin{array}{c}-s_i{J}_r^{\left(i\right)}{\dot{\sigma }}_i\\ 0\\ 0\end{array}\right]+\mathit{\omega }\times {\displaystyle \sum {\mathit{R}}_{\mathit{m}}^{\left(\mathit{i}\right)}\left[\begin{array}{c}-{\mathit{s}}_i{\mathit{J}}_{\mathit{r}}^{\left(\mathit{i}\right)}{\sigma }_{\mathit{i}}\\ 0\\ 0\end{array}\right]}+{\mathit{M}}_{\mathit{a}}\end{array}$$

(9)where T_i is the thrust produced by the ith rotor; ${\mathit{f}}_a$ is the aerodynamic force represented in body frame; R is the rotation matrix that converts a vector represented in the body frame to its representation in the inertial frame; ${\mathit{r}}_i$ is the location of the ith rotor in body frame; ${\mathit{R}}_m^{\left(i\right)}$ is the rotation matrix that converts a vector represented in the motor frame to its representation in the body frame; s_i indicates the rotating direction of the ith motor (s_i = 1 means that the rotation direction is same with $x_m^{\left(i\right)}$ ; $s_i=-1$ is the opposing situation); Q_i is the torque produced by the ith motor; $J_r^{\left(i\right)}$ is the inertial matrix of the ith rotating part of the motor and the attached propeller; σ_i is the rotation speed of the ith motor.

From equation (9), it can be seen that the total force ${\mathit{F}}_t$ consists of rotor thrust, aerodynamic force, and gravity. Similarly, the total moment ${\mathit{M}}_t$ comprises rotor thrust and torque, rotor acceleration and deceleration, gyroscopic effect of rotating rotors, and aerodynamic moment. The rotor related force and moment will be further discussed in next section and the airframe aerodynamic force and moment will be discussed in the Airframe aerodynamics section.

Actuator dynamics

According to Brandt and Selig, ²⁷ the propeller thrust and torque can be parameterized as follows

$$\begin{array}{c}T=C_T\left(J,\sigma \right)\rho {\sigma }^2{D}^4\\ Q=C_Q\left(J,\sigma \right)\rho {\sigma }^2{D}^5\end{array}$$

(10)where C_T and C_Q are, respectively, the thrust and torque coefficients; ρ is the air density; D is the propeller diameter; σ is the rotor speed in revolutions per minute (RPM) of the rotor; J is the advanced ratio which can be written as

$$J=\frac{V_f}{\sigma D}$$

(11)where V_f is the propeller forward speed (the magnitude of inflow air velocity along $x_m^{\left(i\right)}$ ).

To parameterize C_T and C_Q, a wind tunnel test was conducted with the test platform shown in Figure 2. The motor–propeller pair is directly installed on a six-axis force sensor to measure its thrust and torque. A constant voltage source is used to provide and measure the power (i.e. voltage and current) for the tested motor and propeller. Through a specially designed strut, the whole test platform is connected to a rotation table, which is used to adjust the orientation of the inflow air. The rotation axis of the rotor x_m is designed to be parallel to the ground. The wind speed generated by the wind tunnel is adjusted by tuning the fan voltage. A pitot tube (shown later in Figure 6) is used to measure the airspeed (i.e. the forward speed V_f for the propeller) of the wind tunnel. Through changing the PWM of the electric speed controller (ESC), we can adjust the speed of the rotor. The relationship among rotor speed, PWM, and forward speed V_f is shown in Figure 3. The data measurement process is as follows: (1) given a preset fan voltage and PWM, we wait for 30 s to make sure that the wind in the wind tunnel is constant and the propeller is rotating at a steady speed and (2) record the corresponding data with 10 s to get the mean value.

Figure 2.

The propulsion system wind tunnel test.

Figure 3.

The relationship among forward speed V_f, PWM, and rotor speed.

Similarly, we can obtain the relationship of thrust, torque, and power (the total power consumption of the propulsion system, gotten by multiplying the current and voltage), respectively, versus rotor speed and forward speed as shown in Figure 4(a), (c), and (e) by varying rotor speed and forward speed. Figure 4(a) and (c) is a two-dimensional graph with a legend of different constant wind velocities. Figure 4(e) is a 3D graph with the raw data and fitted curve. Using the results from Figure 4(a) and (c), we can further parameterize C_T and C_Q with respect to advance ratio J and rotor speed. MATLAB $fit$ () function is used to parameterize the data as a polynomial curve. The parameterized results of C_T and C_Q are represented as

$$\begin{array}{l}{C}_T\left(J,\sigma \right)=t_{00}+t_{10}J+t_{01}\sigma +t_{20}{J}^2+t_{11}J·\sigma +t_{02}{\sigma }^2\\ +t_{30}{J}^3+t_{21}{J}^2·\sigma +t_{12}J·{\sigma }^2+t_{03}{\sigma }^3\\ t_{00}=0.04438,t_{10}=-0.08226,t_{01}=0.001466,\\ t_{20}=-0.01371,t_{11}=-0.003201,t_{02}=-0.004358,\\ t_{30}=0.005713,t_{21}=-0.00722,t_{12}=-0.002443,\\ t_{03}=0.002003\end{array}$$

(12)

$$\begin{array}{l}{C}_Q\left(J,\sigma \right)=q_{00}+q_{10}J+q_{01}\sigma +q_{20}{J}^2+q_{11}J·\sigma \\ +q_{02}{\sigma }^2+q_{30}{J}^3+q_{21}{J}^2·\sigma +q_{12}J·{\sigma }^2\\ q_{00}=0.005461,q_{10}=-0.003079,q_{01}=-0.000154\\ q_{20}=-0.00083,q_{11}=-1.464e-06,q_{02}=-8.608e-05\\ q_{30}=0.000395,q_{21}=1.511e-05,q_{12}=-2.085e-05\end{array}$$

(13)which will be used in the simulation implementation. The corresponding fitted surfaces are shown in Figure 4(b) and (d), which will be used in the Flight envelope simulation section to estimate the total power consumption of our vehicle in the simulation.

Figure 4.

Propulsion system test data and fitting results. (a) Relationship among thrust, rotor speed, and forward speed, (b) Relationship among C_T, rotor speed, and advance ratio J, (c) Relationship among torque, rotor speed, and forward speed, (d) Relationship among C_Q, rotor speed, and advance ratio J, and (e) Relationship among power, rotor speed, and forward speed.

According to Theys et al., ²⁸ the inflow angle will also affect the thrust and moment of the propeller. In the worst case, the error is approximately 10%. Currently, we only consider the inflow along x_m axis, which is the principal factor to the thrust and moment of the propeller.

Airframe aerodynamics

The aerodynamic force ${\mathit{f}}_a$ and moment ${\mathit{M}}_a$ acting on the vehicle is shown in Figure 5(a), which can be represented as follows

$$\begin{array}{c}{\mathit{f}}_a=\left[\begin{array}{ccc}-\cos \alpha & 0& \sin \alpha \\ 0& 1& 0\\ -\sin \alpha & 0& -\cos \alpha \end{array}\right]\left[\begin{array}{c}\mathcal{D}\\ \mathcal{Y}\\ \mathcal{L}\end{array}\right]\\ {\mathit{M}}_a={\left[\begin{array}{lll}L\hfill & M\hfill & N\hfill \end{array}\right]}^T\end{array}$$

(14)where α is the angle of attack; $\mathcal{L},\mathcal{D}$ , and $\mathcal{Y}$ are, respectively, lift, drag, and side force represented in stability frame; L, M, and N are, respectively, aerodynamic roll, pitch, and yaw moment along body axis. According to Etkin and Reid, ²⁹ the aerodynamic force $\mathcal{L},\mathcal{D},\mathcal{Y}$ and aerodynamic moment L, M, N can be, respectively, parameterized as follows

$$\begin{array}{c}\mathcal{L}=\frac{1}{2}\rho V^{2}S{C}_L\left(\alpha ,\beta \right)\\ \mathcal{D}=\frac{1}{2}\rho V^{2}S{C}_D\left(\alpha ,\beta \right)\\ \mathcal{Y}=\frac{1}{2}\rho V^{2}S{C}_Y\left(\alpha ,\beta \right)\\ L=\frac{1}{2}\rho V^{2}S\overline{c}{C}_l\left(\alpha ,\beta \right)\\ M=\frac{1}{2}\rho V^{2}S\overline{c}{C}_m\left(\alpha ,\beta \right)\end{array}\begin{array}{c}N=\frac{1}{2}\rho V^{2}S\overline{c}{C}_n\left(\alpha ,\beta \right)\end{array}$$

(15)where ρ is the air density; V is the magnitude of the airspeed vector u ; S is the wing area; $\overline{c}$ is the vehiclede mean aerodynamic chord; C_L, C_D, and C_Y are, respectively, the lift, drag, and side force coefficients; and C_l, C_m, and C_n are, respectively, the roll, pitch, and yaw coefficients. The aerodynamic coefficients are all related to the angle of attack α and sideslip angle β. The airspeed vector u represented in body frame is the combination of ground speed v and wind speed w (both are represented in inertial frame), which can be written as

$$\mathit{u}={\mathit{R}}^T\left(\mathit{v}-\mathit{w}\right)$$

(16)

Figure 5.

Aerodynamic frame, force, and moment. (a) Aerodynamic frame, force, and moment 3D and (b) Aerodynamic frame, force, and moment in longitudinal direction.

The angle of attack α and sideslip angle β can be calculated as follows

$$\begin{array}{c}\mathit{u}={\left[\begin{array}{lll}{u}_x,\hfill & u_y,\hfill & u_z\hfill \end{array}\right]}^{T}V=\Vert \mathit{u}\Vert \\ \alpha =\tan -1\left(\frac{u_z}{u_x}\right)\beta =\sin -1\left(\frac{u_y}{V}\right)\end{array}$$

(17)

In order to calculate the aerodynamic force and moment at all possible angles of attack and sideslip angles that are likely to encounter in real flight, a full-scale vehicle (without propellers) wind tunnel test is conducted to parameterize the aerodynamic coefficients, as shown in Figure 6. The vehicle is fixed at the center of the wind tunnel on a rotation table that is used to adjust the vehicle yaw angle (sideslip angle). The UAV pitch angle can be adjusted by a rotation mechanism installed right below the UAV (see Figure 6 full-scale wind tunnel test). The wind speed level is adjusted by the voltage input to the wind tunnel fan. A six-axis force sensor is installed at the center of the rotation table to measure the aerodynamic force and moments. And a pitot tube is placed at the UAV shoulder to measure the actual airspeed. Similar to the data gathering method for the propeller, we collect the mean value of the measured data, including the wind speed measured by the pitot tube, UAV pitch angle read from the rotation mechanism, UAV yaw angle read from the rotation table, and the aerodynamic forces and moments measured by the six-axis force sensor, after the wind speed produced in the wind tunnel is steady. With the UAV pitch and yaw angle, we compute its angle of attack and sideslip angle. Meanwhile, we compute the aerodynamic lift, drag, side force, and moments by projecting the measured forces and moments to the corresponding directions. With these aerodynamic forces, moments, wind speed, and UAV wing area as the reference area, we compute the aerodynamic force and moment coefficients by equation (15). In the experiments, we adjust the rotation mechanism and also the rotation table such that the angle of attack reaches the full range of ( $-{180}^{°}\sim {180}^{°}$ ) and sideslip angle reaches the full range of 0 to 90° (0 to –90° is symmetric). The wind speed is from 2.9 m/s to 18.9 m/s, sufficient to cover the full operating speed of the UAV. The measured data is then interpolated by fitting it to a group of carefully chosen basis functions. The fitted aerodynamic force and moment coefficients are shown in Figure 7 which will be used in the simulation platform in the future. Since the transition flight mainly takes place in longitudinal direction, we particularly investigate the longitudinal aerodynamic coefficients (see Figure 8) for further use in the controller design in the Controller design section. The corresponding aerodynamic force in the longitudinal direction can be represented as

$${\mathit{f}}_{al}=\left[\begin{array}{ccc}-\cos \alpha & 0& \sin \alpha \\ 0& 1& 0\\ -\sin \alpha & 0& -\cos \alpha \end{array}\right]\left[\begin{array}{c}\mathcal{D}\\ 0\\ \mathcal{L}\end{array}\right]$$

(18)

Figure 6.

Full-scale wind tunnel test.

Figure 7.

The interpolated aerodynamic coefficients. (a) Lift coefficient C_L, (b) Roll coefficient C_l, (c) Drag coefficient C_D, (d) Pitch coefficient C_m, (e) Side force coefficient C_Y, and (f) Yaw coefficient C_n.

Figure 8.

Longitudinal aerodynamic coefficients (a) Lift coefficient C_L. and (b) Drag coefficient C_D.

The comprehensive modeling and wind tunnel test lead to a complete and accurate flight simulator. We will defer the simulator development in the Simulation implementation section and look into the controller design in the sequel discussion.

Controller design

Attitude and transition controller design

Our controller consists of two parts: an attitude controller and a transition controller. The attitude controller is used to track the attitude to the desired one in all flight envelopes. As shown in Figure 9, the attitude controller is a cascaded controller based on SO(3). The outer loop and the inner loop are, respectively, angular loop and angular rate loop. The outer loop is a proportional controller, the desired angular rate ${\mathit{\omega }}_d$ is

$${\mathit{\omega }}_d=-{\mathit{k}}_p{\left(\log \left({\mathit{R}}_e\right)\right)}^{\text{∨}},{\mathit{k}}_p>0$$

(19)where ${\mathit{k}}_p$ is the positive definite diagonal matrix, the vee map ${(\cdot )}^{\text{∨}}$ is the inverse operation of the hat map $\widehat{\dot{}},\log \left({\mathit{R}}_e\right)$ can be calculated from

$$\begin{array}{l}\log \left({\mathit{R}}_e\right)=\frac{\varphi }{2\sin \varphi }\left({\mathit{R}}_e-{\mathit{R}}_e^T\right)\\ \cos \varphi =\frac{1}{2}\left(tr\left(\mathit{R}\right)-1\right),\left|\varphi \right|<\pi \\ {\mathit{R}}_e={\mathit{R}}_d^T\mathit{R}\end{array}$$

(20)According to Bullo and Murray, ³⁰ this attitude control algorithm will be exponential stable if the angular velocity can be controlled to the desired value instantly. The angular rate loop is a PID controller with feed forward terms compensating the aerodynamic moment, which can be written as

$${\mathit{M}}_T={\mathit{K}}_p{\mathit{\omega }}_e+{\mathit{K}}_i{\displaystyle \int {\mathit{\omega }}_e\text{d}t+{\mathit{K}}_d\frac{\text{d}{\omega }_e}{\text{d}t}+\widehat{\mathit{\omega }}{\mathit{J}}_b\mathit{\omega }-{\mathit{M}}_a}$$

(21)where ${\mathit{K}}_p,{\mathit{K}}_i,{\mathit{K}}_d$ are positive diagonal gains for the PID controller, ${\mathit{\omega }}_e={\mathit{\omega }}_d-\mathit{\omega }$ is the angular velocity error. The Ziegler–Nichols method is used to tune the controller gain, respectively, in the roll, pitch, and yaw directions (i.e. the first parameter of ${\mathit{K}}_p,{\mathit{K}}_i,\text{and}{\mathit{K}}_d$ for roll direction, the second and third parameters of the matrices are, respectively, corresponding to the pitch and yaw directions).

Figure 9.

Attitude controller diagram.

The transition controller is developed to hold the UAVdl altitude during the transition process, which is particularly useful for accurate and stable transition flight. The altitude dynamics can be represented as follows

$$\begin{array}{l}{\ddot{h}}_z=\frac{1}{m}{\mathit{r}}_3^T\left({\mathit{f}}_{al}+{\mathit{f}}_T\right)+g\\ \mathit{R}=\left[\begin{array}{ccc}{\mathit{r}}_1& {\mathit{r}}_2& {\mathit{r}}_3\end{array}\right]\end{array}$$

(22)where ${\ddot{h}}_z{\hspace{0.17em}}_z$ is vehicle altitude represented in the inertial frame. Assuming the desired altitude is $h_d$ , the altitude error can be represented as $h_e=h_d-h_z$ . A PID controller is used to compute the altitude control actuation as follows

$${\mathit{r}}_3^T\left({\mathit{f}}_{al}+{\mathit{f}}_T\right)+mg=k_p{h}_e+k_i{\displaystyle \int h_e}\text{d}t+k_d\frac{\text{d}{h}_e}{\text{d}t}$$

(23)

The resulting thrust is therefore

$$f_{Tx}=\frac{k_p{h}_e+k_i{\displaystyle \int h_e}\text{d}t+k_d\frac{\text{d}{h}_e}{\text{d}t}-mg-{\mathit{r}}_3^T{\mathit{f}}_{al}}{r_{31}}$$

(24)where $r_{31}$ is the first entry of vector ${\mathit{r}}_3\in {\mathbb{ℝ}}^3,{\mathit{f}}_{al}$ is the aerodynamic force in the longitudinal direction.

Mixer

After the control thrust f_Tx and moment M _T are calculated, a mixer is used to allocate the rotor thrust. Considering the motor inclination design of our vehicle, the rotor thrust can be written as

$$\begin{array}{l}\left[\begin{array}{c}{T}_1\\ T_2\\ T_3\\ T_4\end{array}\right]=\frac{1}{4}\left[\begin{array}{cccc}1& -1& 1& -1\\ 1& 1& 1& 1\\ 1& -1& -1& 1\\ 1& 1& -1& -1\end{array}\right]\left[\begin{array}{c}\frac{f_{Tx}}{c\left({\beta }_m\right)}\\ \frac{M_{Tx}}{s\left({\beta }_m\right)d+\kappa c\left({\beta }_m\right)}\\ \frac{M_{Ty}}{\left(c\left({\beta }_m\right)d-\kappa s\left({\beta }_m\right)\right)s\left({\beta }_l\right)}\\ \frac{M_{Tz}}{\left(c\left({\beta }_m\right)d-\kappa s\left({\beta }_m\right)\right)c\left({\beta }_l\right)}\end{array}\right]\\ \kappa =\frac{C_Q\left(J,\sigma \right)}{C_T\left(J,\sigma \right)}\end{array}$$

(25)where $c(·)$ and $s(·)$ are, respectively, short for $\cos (·)$ and $\sin (·)$ , d is the arm length of each rotor, which is half the length of motor diagonal as shown in Figure 1. κ is the ratio between the torque and trust coefficient, which varies with J and rotor speed as shown in Figure 10. According to equation (25), the varying κ affects the calculation of rotor thrust. However, since the motor speed is not measurable due to the commercial-off-the-shell motor and ESC being used, we temporarily fix κ at a constant value corresponding to J = 0 and σ = 4500 (i.e. the hovering situation) in actual implementation. After getting the rotor thrust, equation (10) is used to get the rotor speed with J = 0, and the relationship between rotation speed and PWM with V_f = 0 (shown in Figure 11) is used to calculate the PWM command.

Figure 10.

The relationship among $\frac{C_Q}{C_T}$ , J, and rotor speed.

Figure 11.

The relationship between rotor speed and PWM with V_f = 0.

Simulation implementation

In the Modeling section, we have built the completed model of our vehicle and obtained all the required parameters through corresponding experiments. Combining the kinematics and dynamics model, the designed controllers, and the mixer, we implement a flight simulator in MATLAB Simulink as shown in Figure 12. The simulation is highly modularized, including Joystick, Transition Controller, Attitude Controller, VTOL Mixer, VTOL Dynamics, and Visualization module. The Joystick module is used to send the human pilot commands including attitude command, altitude command, and transition command. The Transition Controller module is used to calculate the thrust and send attitude command during the transition process. The Attitude Controller module receives the attitude command to calculate the desired control moment. The desired thrust and moment are then sent to the VTOL Mixer module to allocate the control action. The PWM commands are finally sent to VTOL Dynamics module which has been detailed in the Modeling section. The Visualization module can simultaneously show the attitude and position of the vehicle as shown in Figure 13, which demonstrates an example transition from vertical flight to level flight and then back to vertical flight. After the simulation is completed, all the simulation results are saved as a flight log for further analysis.

Figure 12.

The VTOL UAV simulation platform implemented in Simulink.

Figure 13.

The simulation visualization results.

Transition flight simulation

As we mentioned earlier, the vehicle is capable of both level flight with better power efficiency and vertical flight for hovering. To achieve this, a transition maneuver between these two distinctive flight modes is required. We show in the simulation that the previously designed controller can achieve this goal very well. In the simulation, the attitude controller is fed with a linear decreasing pitch angle command for the forward transition flight (from vertical flight to level flight) from the transition controller. Similarly, a command with a linear increasing pitch angle is used for the backward transition (from level flight to vertical flight). The simulation results are shown in Figure 14(a). During the transition process, the commands of roll, yaw, and altitude are all set to be constant. We can see from the simulation results that the pitch angle can be tracked very well except the backward transition process. However, the roll and yaw angles can be well tracked with the maximum absolute error less than $2^{°}$ and ${1.5}^{°}$ , respectively. Looking into the altitude response, we can see that the tracking error is rather small (i.e. within 0.15 m) at the forward transition flight. During the backward transition flight, the vehicle altitude increases by around 2 m, due to the excessive lift force created when the aircraft pitches up.

Figure 14.

Transition between level flight and vertical flight. (a) Simulation results and (b) Experimental results.

Flight envelope simulation

In the simulation, we set a series of commanded pitch angles with a constant altitude command. With the designed attitude and altitude controller designed in the Controller design section, the aircraft can track the desired pitch angles quite well without any significant altitude variation. Once reaching a steady state, the aircraft flight speed can be recorded. This flight speed is equal to the airspeed as no wind is simulated. As a consequence, the steady flight speed versus different pitch angle (i.e. the angle of attack) is obtained. The simulation results are shown in Figure 14, which shows that our vehicle can fly at any angle of attack between ${°}^5$ and ${90}^{°}$ . The simulation results also show that our vehicle can continuously reach cruise speed between 0 m/s to 15 m/s, which makes our vehicle suitable for different mission tasks.

Figure 15.

PWM of the motors for Figure 14(b).

Using the relationship between thrust and airspeed at different angles of attack as shown in Figure 16 and the fitted rotor parameters in Figure 4(e), the power consumption of our vehicle is estimated as follows:

First, we obtain the required total thrust at steady flight from Figure 16. And divide this total thrust by 4 to obtain the required thrust of each motor. By doing this, we are assuming that the vehicle flies at the trim condition which means the aerodynamic moment equals to zero. Thus, four rotors are considered to be working in the same condition.

From the same Figure 16, the steady level flight speed is read, based on which we calculate the forward speed of each rotors as $V_f=V\cos \left(\alpha \right)$ .

Using the relationship among the forward velocity, thrust, and rotation speed in Figure 4(a), we estimate the speed for each rotor.

Provided with rotation speed and forward speed in previous steps, we finally estimate the total power consumption using the fitted power consumption shown in Figure 4(e).

The estimated power consumption is shown in Figure 17, where we can see that the power decreases from 132.5 W to 76.02 W while the vehicle pitch forward from ${90}^{°}$ to ${16}^{°}$ , the corresponding airspeed is 9.5 m/s. When the angle of attack continues to decrease, the power increases rapidly (164.1 W at 7° (the corresponding airspeed is 12.7 m/s) which is more power consuming than hovering flight). This is due to the fact that we are using a propulsion system designed for multi-rotors. Thus, the rotor efficiency degrades dramatically with high forward speed V_f, which will cause the increase of power consumption. At the same time, as shown in Figure 18 the lift–drag ratio decreases dramatically while the angle of attack is below 9°. Thus, the rotors need to provide larger thrust to overcome the increasing drag as the decrease of angle of attack is below 9°, which will also lead to high power consumption.

Figure 16.

The relationship between angle of attack, velocity, and thrust.

Figure 17.

The relationship between estimated power and angle of attack/airspeed in simulation.

Flight experiments

Transition flight verification

During the transition flight verification process, a decreasing and increasing pitch command was, respectively, used to achieve forward and backward transition, similar to that used in the simulation. The experimental results are shown in Figure 14(b), which shows that our vehicle can accomplish the transition flight without gaining/dropping significant altitude (gains 0.38 m in the forward transition flight, drops 1.25 m in the level flight, and gains 1.85 m in the backward transition flight) and the attitude can be well tracked with maximum roll error less than 5° and yaw error less than 3.5°.

Comparing Figure 14(a) and (b), we can find out that the commanded attitude and altitude in simulation and experiment are the same. The transition flight verification shows that our simulation platform can predict the flight status of our vehicle quite well from the following aspects: (1) The pitch tracking performance is very similar for simulation and experiments, for both forward and backward transitions. For example, both the simulation and experiments show accurate pitch tracking during the forward flight, but a considerable tracking delay during backward transition. (2) The simulation shows that the vehicle gains around 2 m during the backward transition, which is also observed in the real flight. (3) Both roll and yaw directions are tracked quite well during the whole process.

Figure 15 shows the commanded PWM of the outdoor flight test for Figure 14(b). The red line and yellow line are, respectively, the upper bound and lower bound of the commanded PWM. As we can see, during the backward transition, the altitude of the vehicle has increased to around 2 m. To hold the desired altitude, the commanded PWMs are quite small (from 12 s to 14 s), which will decrease the control moment. However, the aerodynamic moments in pitch direction are quite significant due to the high airspeed at the beginning of the backward transition. The large aerodynamic moments exceeds the control moment provided by the motor–propeller pair leads to the pitch tracking error.

Figure 18.

The relationship between lift-drag ratio and angle of attack.

Flight envelope verification

Similar to the flight envelope simulation detailed in the Flight envelope simulation section, we conduct several flight experiments and set different pitch commands as the angle of attack in each experiment. The corresponding angle of attack and steady cruise speed are recorded during the flight experiment. A pitot tube is rigidly installed along x_b axis to measure the airspeed. However, this will lead to inaccurate and sometimes false airspeed measurement at high angle of attack. To solve this problem, we conducted our experiment in a clear day with very low wind. In this way, the vehicle ground speed, which can be measured by its global positioning system (GPS), can be used as the airspeed. To maximally cancel out the effect of residual wind, we conduct our experiments as that shown in Figure 19. We first measure the wind direction at any wind speed and then command the UAV to fly along and against the wind reversely to get the mean cruise speed. The results are summarized in Figure 20. It can be seen that the actual flight results agree with our simulation surprisingly well, which implies the effectiveness of our simulation model. The slight residual error may be caused by the uncompensated residual wind and ground speed measurement errors. Figure 21 plots the log data of some sample flights with ${60}^{°},{30}^{°}$ , and ${10}^{°}$ pitch commands, respectively.

Figure 19.

Outdoor test methodology.

Figure 20.

The relationship between angle of attack and airspeed.

Figure 21.

Vehicle ground speed versus different pitch angles during outdoor test. (a) Ground speed versus pitch angle 60°, (b) Ground speed versus pitch angle 30°, and (c) Ground speed versus pitch angle 10°.

Conclusion

In this paper, we build a detailed model of our quadrotor tail-sitter, including the kinematics, rigid body dynamics, actuator dynamics, and airframe aerodynamics. Extensive wind tunnel tests are conducted to obtain the coefficients of actuator dynamics and airframe aerodynamics. A transition controller and a unified attitude controller are designed to accomplish the vertical flight, level flight, and transition flight. A Simulink simulation platform is implemented and the corresponding transition flight simulation and flight envelope simulation are conducted. Simulation results that our vehicle can accomplish transition between hovering flight and level flight with insignificant altitude error and the ability to fly at a wide flight envelope from ${°}^5$ to ${90}^{°}$ , covering a wide flight speed envelope from stationary hovering to fast level flight. The corresponding experimental verification supports the simulation results and further prove the accuracy of our simulation platform. Some transition flight tests are shown in https://youtu.be/S3BVDLEFjGM. Future works will focus on the wind estimation and the design of the anti-wind controller, especially when hovering where the hovering accuracy is particularly important for the accurate landing.