Next Article in Journal
Potential-Field-RRT: A Path-Planning Algorithm for UAVs Based on Potential-Field-Oriented Greedy Strategy to Extend Random Tree
Previous Article in Journal
Research on Environment Perception System of Quadruped Robots Based on LiDAR and Vision
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model, Control, and Realistic Visual 3D Simulation of VTOL Fixed-Wing Transition Flight Considering Ground Effect

1
Department of Computer Science and Electronics, Faculty of Mathematics and Natural Sciences, Universitas Gadjah Mada, Yogyakarta 55281, Indonesia
2
Department of Aerospace Engineering, Sekolah Tinggi Teknologi Kedirgantaraan, Yogyakarta 55188, Indonesia
*
Author to whom correspondence should be addressed.
Drones 2023, 7(5), 330; https://doi.org/10.3390/drones7050330
Submission received: 6 April 2023 / Revised: 3 May 2023 / Accepted: 19 May 2023 / Published: 20 May 2023
(This article belongs to the Section Drone Design and Development)

Abstract

:
The research topic of VTOL (vertical take-off and landing) fixed wing (VFW) is gaining significant attention, particularly in the transition phase from VTOL to fixed wing and vice versa. One of the latest and most challenging transition strategies is the bird take-off mode, where vertical and horizontal take-off is carried out simultaneously, mimicking the behavior of birds. The condition that is rarely considered when taking off is the ground effect. Under natural conditions, a ground effect is bound to occur, which can significantly impact the stability of the transition when the VFW is close to the ground. This paper addresses this issue by proposing a model and control strategy and conducting realistic visual 3D simulations of the VFW transition that incorporates ground effect using full complex aerodynamic parameters. This research represents a novel approach, using the robot operating system (ROS) and Gazebo to conduct realistic visual 3D simulations for VFW transition. The linear quadratic regulator (LQR) control method is used to manage the transitions and compensate for any disturbances. The flight tests demonstrate the effectiveness of the proposed model and controller in executing flight missions using the bird take-off mode transition. Moreover, the controller has demonstrated reliability and robustness in compensating for attitude errors induced by ground effects and external disturbances.

1. Introduction

Unmanned aerial vehicles (UAVs), also known as drones, have become increasingly prevalent in military and civilian applications due to their numerous advantages. These advantages include low risk to human life, relatively inexpensive operational costs, and ease of use. UAVs come in various categories, such as multirotor and fixed-wing. Multirotor UAVs have vertical take-off and landing (VTOL) capabilities and can hover in place, but their cruising range is relatively short [1,2]. In contrast, fixed-wing UAVs can travel long distances at high speeds but require a long runway for take-off and landing [3,4].
Hybrid UAVs combine the advantages of both multirotor and fixed-wing types. They can take off and land vertically like multirotor and travel long distances like fixed-wing UAVs. Depending on their propulsion type, hybrid UAVs can be further divided into two categories: those with the same propulsion for vertical and horizontal flight and those with separate propulsion for each flight mode. The former includes tilt-wing, tilt-rotor, and tail sitter UAVs, while the latter includes VTOL fixed wing (VFW). There is a complicated mechanism for the same propulsion and a complex control system, so it is not used as a design reference. Unfortunately, the VFW is simple to implement and can be further expanded with quadrotors and fixed-wing configuration [5,6].
The advantage of VFW is that their geometry remains unchanged during the transition from VTOL flight mode to fixed-wing mode and vice versa, making their mechanical motions simple and straightforward [5,6,7,8]. However, the weakness of the VFW is that it has more propulsion than the same propulsion. This condition causes the VFW to be heavier than the same propulsion type. In addition, the VFW requires control of two propulsion systems for vertical take-off and cruising simultaneously during the transition process.
In general, VFW has four stages of flight; take-off, transition, cruising, and landing. Figure 1 shows the conventional transition stages. The movement from point A to point B is a take-off process. Meanwhile, moving from point B to point C is a transition process. The transition is the most critical process [5] in VFW. Point C is the end of the transition process and the beginning of cruising. The transition ends when the airspeed has exceeded the fixed wing’s stall speed [9].
The strategy above is safe enough because the transition process begins when the VFW is at cruising altitude [9]. However, vertical take-off to reach cruising altitude requires significant power, so it will shorten the cruising range. Based on the flight stages, power consumption from the largest to the smallest is taking off, landing, climbing, and cruising, respectively. Vertical take-off on a VFW requires much power because it rotates the four rotors [10,11].
To improve power efficiency, developing effective transition strategies for VFW UAVs is crucial. One promising approach is to emulate a bird’s take-off, as demonstrated in Figure 2. This innovative strategy involves a simultaneous vertical take-off and transition from point A by activating all four quadrotor propulsions to generate lift and fixed-wing propulsion to provide forward thrust. By combining these actions, the UAVs can achieve a more efficient and effective transition, reducing power consumption and improving flight performance [12]. These findings have important implications for the future of VFW technology, as power efficiency is an essential aspect to consider in the design and development of VFW. With this promising approach to improve power efficiency, VFW can continue to push the boundaries of aerial capabilities and contribute to the advancement of the field of unmanned aerial vehicles. However, one important factor to consider is the ground effect because the transition starts when VFW is still close to the ground.
In the context of UAVs, a phenomenon known as the ground effect occurs when the rotor UAVs operate at low altitudes. This phenomenon is characterized by a downwash propeller effect reflected by the ground, resulting in increased lift force but compromised flight stability. Prior research [13] has shown that the ground effect affects the flight stability of rotor UAVs during take-off when they are close to the ground. As the UAV moves further away from the ground, the ground effect diminishes rapidly, eventually disappearing altogether, thus destabilizing flight. It is a challenge for the bird take-off mode transition to carry out the stable transition when the ground effect occurs.
The ground effect on the VFW has not yet received attention from previous studies, both at the simulation and experimental stages. The ground effect in helicopters has been well-researched and studied for take-off, landing, and hovering near the ground [14,15,16]. In contrast, the ground effect in multirotor has less attention research [17,18]. Furthermore, while recent research has evaluated ground effects in VTOL planes, these studies have only focused on conventional take-off and have not explored the impact on flight stability [19]. This gap in research motivates this study, which aims to explore and apply the ground effect as a critical factor influencing the stability of VFW systems during bird take-off mode transition. By exploring this previously unexplored area of research, it is possible to advance the understanding of the complex interplay between flight dynamics and environmental factors and contribute to the development of more robust and adaptable flight control systems.
During the transition of bird take-off mode, the VFW needs a controller to maintain flight stability. Even more, there is a ground effect that affects flight stability. In autonomous VFW, a controller is at least needed to control attitude and trajectory to achieve stable flight [20,21]. Based on several studies [22,23,24], linear quadratic regulator (LQR) control has robustness for multirotor control. In addition, LQR also has satisfactory performance for fixed-wing control [25,26,27]. In this context, LQR control has the advantage of being able to handle nonlinear dynamics and uncertainties. Moreover, the feedback loop inherent in LQR control can help stabilize the VFW during the bird take-off transition mode, where the VFW experiences a significant change in its flight dynamics. Therefore, the application of LQR control for VFW, particularly during the bird take-off transition mode, is a promising solution to address the challenges of maintaining stability in the presence of ground effect and external disturbance.
The application of aerodynamics parameters in Gazebo 3D visual simulation for flapping wings has been discussed [28]. However, it is just a brief analysis, not the simulation application stage, because the flapping and morphing wings have airfoils that move relative to global inertia and body inertia to produce aerodynamic forces. This phenomenon is the same as in the multirotor. Furthermore, previous studies [29,30,31,32] simulated UAVs using Gazebo but did not apply full aerodynamic parameters and did not yield satisfactory outcomes. Different from previous studies, this research applies complete aerodynamic parameters to generate force and torque so that the simulation is like real conditions by considering the ground effect.
Overall, the contribution of this paper is summarized as follows:
  • The paper presents a dynamic model of VFW using the Newton–Euler approach and bird take-off mode transition trajectory.
  • The paper also proposes a ground effect model that influences flight stability during the transition. This model is an important consideration that was not previously addressed in similar research.
  • To handle the bird take-off mode transition and compensate for disturbances, this study applies a robust LQR controller.
  • Finally, the paper conducts a realistic 3D simulation of VFW’s bird take-off mode transition using Gazebo/ROS. This simulation is an essential step in verifying the model and robustness of the controller in the presence of external disturbances and ground effects.
This paper is organized into five sections, as follows; Section 2 describes VFW mathematical model, the ground effect model, and the system dynamic and aerodynamic parameters; Section 3 presents the state space model and the controller; Section 4 presents the visual 3D simulation using bird take off mode transition in any conditions. Finally, the conclusions are given in Section 5.

2. System Modeling and Parameters

Section 2.1 presents VFW modeling using the Newton–Euler approach, Section 2.2 explains the ground effect model, and Section 2.3 describes the aerodynamics parameters for realistic 3D visual simulation in ROS/Gazebo.

2.1. Mathematical Model of VTOL Fixed Wing

VFW modeling refers to the Earth’s inertial frame G: { O G ; X G , Y G , Z G } and the body’s inertial frame B: { O B ; X B , Y B , Z B }. According to the rotation matrix equation, GRB will represent the orientation of the body frame B to the earth frame G. GRB ∈ SO(3) is the rotation matrix representing the orientation of the body frame on the earth frame is presented as
R G B = c ψ c θ s ψ c ϕ + c ψ s θ s ϕ s ψ s ϕ c ψ s θ c ϕ s ψ c θ c ψ c ϕ + s ψ s θ s ϕ c ψ s ϕ + s ψ s θ c ϕ s θ c θ s ϕ c ϕ c θ
where c is cos, s is sin, θ is pitch angle, ϕ is roll angle and ψ is yaw angle.
As shown in Figure 3, the VFW was built using a Skysurfer airframe and attaching an H-type quadrotor frame. Skysurfer is a slow flayer fixed wing with a 1400 mm wing span and 915 mm length. Skysurfer has ailerons, an elevator, a rudder, and one puller propulsion. The force acting on the VFW during the transition involved two propulsion systems:
F = F Q + F E .
The subscript Q denotes the quadrotor, and the subscript E denotes the fixed wing. The resultant force F is the sum of forces from the quadrotor and fixed wing, while the force on each axis is
F x F y F z = F Q x F Q y F Q z + F E x F E y F E z
Therefore, the working torque is the sum of two torques from the quadrotor and fixed-wing torque is defined as
M = M Q + M E
and the torque on each axis is
M x M y M z = M Q x M Q y M Q z + M E x M E y M E z
The VFW dynamics can be expressed in the following form [6]
F x = F Q x + F E x m g   s i n   θ = m u ˙ + q w r v
F y = F Q y + F E y m g   c o s   θ . cos   ϕ = m v ˙ + r u p w
F z = F Q z + F E z + m g   c o s   θ . sin   ϕ = m w ˙ + p v q u
M x = M Q x + M E x = I x x r ˙ I x y q ˙ p r + I z z I y y q r
M y = M Q y + M E y = I y y q ˙ I x y p ˙ q r + I x x I z z p r
M z = M Q z + M E z = I z z r ˙ I x y q 2 p 2 + I y y I x x p q
where F, M, m, and g denote force, torque, mass, and gravity. While u, v, and w represent angular velocity on the x-axis, y-axis, and z-axis, translation velocity along the x-axis, y-axis, and z-axis are p, q, and r. The inertial matrix of the vehicle’s rigid body is described as
I = I x x I y x I x z I x y I y y I y z I x z I y z I z z
When VFW is maneuvering, there will be a change in angular and translational positions. The rate of change of angular position is defined as [7]
d ϕ d t = p + r   cos   ϕ + r   s i n   ϕ   t a n   θ
d θ d t = q   c o s   ϕ r   s i n   ϕ
d ψ d t = q   s i n   ϕ + r   c o s   ϕ   s e c   θ
d X G d t = u   c o s   θ . c o s   ψ + v c o s   ϕ . s i n   θ s i n   ϕ . s i n   ψ + w c o s   ϕ . s i n   ψ + s i n   ϕ . c o s   ψ
d Y G d t = u   s i n   θ . c o s   ψ + v   c o s   ϕ . c o s   ψ w   s i n   ϕ . c o s   θ
d Z G d t = u   c o s   θ . c o s   ψ + v c o s   γ . s i n   θ . s i n   ψ s i n   ϕ . c o s   ψ + w c o s   ϕ . c o s   ψ s i n   ϕ . s i n   θ . s i n   ψ
where X G , Y G , and Z G represent its position on the global axis. The total airspeed of VFW:
V = u 2 + v 2 + w 2
is the resultant of the aircraft’s speed on the Earth’s x-axis, y-axis, and z-axis. X G , Y G , and Z G represent its position on the global axis, while α is the angle of attack and β is the sideslip angle of the fixed wing given by
β = arcsin ( w V )
α = arctan ( v w )
To generate lift, the aircraft needs to maintain a positive pitch and angle of attack relative to the wind axis. The origin of the wind axis is at the vehicle’s center of mass. The direction cosine matrix transforms the inertial frame into a wind frame as follows:
F x F y F z = D C M w b F x a F y a F z a ,   M x M y M z = D C M w b M x a M y a M z a
D C M w b = cos   α . cos   β sin   α cos   α . sin   β sin   α . cos   β cos   α sin   α . sin   β sin   β 0 cos   β
where F x a , F y a , F z a , M x a , M y a , and M z a are the forces and moments in the wind frame.
At the beginning of the bird take-off mode transition, the quadrotor propulsion will dominate in producing lift. As airspeed increases, the fixed wing will produce more lift until the end of the transition. The thrust of each quadrotor propeller F i is calculated as
F i = b ω i 2 ,   i = 1 ,   2 ,   3 ,   4
where the thrust coefficient of each propeller b can be obtained as
b = C T ρ A R 2
ω , C T , ρ , A a n d R represent the thrust coefficient propeller, angular velocity, dimensionless thrust constant, air density, propeller area, and propeller radius, respectively. The total thrust of the quadrotor is calculated as
F Q = F 1 + F 2 + F 3 + F 4 ,   i = 1 ,   2 ,   3 ,   4
The aerodynamic moment resulting from the rotation of each propeller on the quadrotor is described as
M i = k ω i 2 ,   i = 1 ,   2 ,   3 ,   4
This moment is characterized by the aerodynamic moment coefficient, denoted as k and defined as
k = C Q ρ A R 3
which is utilized for yaw control. Furthermore, the differential thrust generated by the propellers, measured from the center of gravity (CoG), induces roll and pitch moments. These moments can be expressed as
M Q x M Q y M Q z = l l k l l k l l k l l k F 1 F 2 F 3 F 4
which quantifies the roll, pitch, and yaw moments acting on the center of mass. In these equations, k , C Q , and l denote the moment coefficient for each propeller, the dimensionless moment constant, and the distance from the center of mass to each propeller, respectively.
In fixed-wing mode, the rotation of the puller propulsion generates forward speed, as described by
F E = b ω E 2
When the aircraft is flying forward, the lift and drag forces acting on the wing can be expressed using
F L = C l ρ 2 A v 2
F D = C d ρ 2 A v 2 .
The deflection of the aileron, elevator, and rudder generates aerodynamic forces, which can be quantified as
F c s A = j   s i n δ c s A
F c s E = j   s i n δ c s E
F c s R = j   s i n δ c s R
The aerodynamic forces coefficient for the aileron, elevator, and rudder is represented by
j = C l ρ 2 A v 2
These forces induce roll, pitch, and yaw moments, which can be calculated by multiplying each force by its distance from the CoG as follows:
M E x = 2 F c s A l x
M E y = F c s E l y
M E z = F c s R l z
By incorporating the quadrotor and fixed-wing models discussed earlier, the total forces acting on each axis are expressed as
F x F y F z = 0 0 m g + 0 0 F Q + F D + F E 0 F L

2.2. Ground Effect

The thrust required for the quadrotor to hover in the ground effect and out-of-ground effect is depicted in Figure 4. In the ground effect, the quadrotor needs less thrust than the out-of-ground effect due to the downwash reflections, as illustrated in Figure 5b. In the presence of a ground effect, the lift generated will increase while the propeller angular velocity remains constant. However, when almost leaving the ground effect there will be a sudden drop in a lift. Consequently, the flight control of the VFW in the ground effect became more intricate [13]. Figure 5a demonstrates that elevating the aircraft height resulted in a more uniform flow field.
The ground effect calculation for the rotor is defined as [33]
F I G E F O G E = 1 1 r 4 h r 2 ,   m a k a   F I G E = F O G E 1 1 r 4 h r 2
and measurable until the distance to the ground h r r = 6 [34,35], where F I G E is the thrust in-ground effect, F O G E is the thrust without out-of-ground effect, r is the propeller radius, h r is the distance from the ground, and K b is coefficient lift body. In the quadrotor, the resultant ground effect will be affected by four rotors is expressed as
F I G E F O G E = 1 r 4 h r 2 r 2 h r d 2 4 h r 2 3 r 2 2 h r d 2 4 h r 2 3 3 2 r 2 h r b 2 4 h r 2 3 K b 1
After applying the ground effect to Equation (40), the force acting on the aircraft becomes
F x F y F z = 0 0 m g + 0 0 F I G E + F D + T E 0 F L

2.3. VFW Simulation Parameters

The VFW 3D model incorporates a fixed-wing airframe based on the Skysurfer and a quadrotor frame utilizing the H550. Isometric front, top, and side views of the model are depicted in Figure 6a, Figure 6b and Figure 6c, respectively, with dimensions expressed in millimeters. The Skysurfer was selected as the fixed-wing airframe due to its low stall speed, which is crucial for the VFW’s intended application. This characteristic also contributes to a shorter and simpler transition process, as documented in reference [36]. The physical and aerodynamic parameters for the model are presented in Table 1.
Proper component placement is critical in achieving an optimal center of mass, which is a crucial factor in vehicle stability. The law of leverage can be used to calculate the center of mass (CoM) based on the principle that the weight on one side of a fulcrum multiplied by its distance from the fulcrum is equal to the weight on the opposite side multiplied by its distance from the fulcrum [37]. To ensure adequate thrust, propeller selection for both fixed-wing and quadrotor configurations should be based on the thrust requirements, as detailed in the performance analysis [38,39,40].
The URDF (unified robot description format) in Algorithm 1 declares the robot parameters necessary for visual simulations of the VFW in ROS. As illustrated in Figure 7, the VFW comprises ten links, namely base_link (wing and fuselage), link_1 (front right propeller), link_2 (front left propeller), link_3 (rear left propeller), link_4 (rear right propeller), link_5 (fixed wing propeller), link_6 (right aileron), link_7 (left aileron), link_8 (elevator), and link_9 (rudder). These links are connected to base_link via joints to form the VFW structure. Each link is defined by four properties, namely inertial body properties, visual properties, collision body properties, and aerodynamic parameters. The inertial body properties include the link’s mass, inertial matrix, and center of mass. Collision body properties are defined by geometry functions with subfunctions meshes, while visual properties are the visible parts of the links that are rendered in the simulator. Finally, the aerodynamic parameters are the parameters that affect the aerodynamic forces acting on each link.
Algorithm 1: URDF Configuration.
<robot name=“quad-plane”>
<link
name=“base_link”>
<inertial> <origin … /> <mass … /> <inertia … /> </inertial>
<visual> <origin … /> <geometry> … </geometry> <material … /> </material> </visual>
<collision> <origin … /> <geometry> … </geometry> </collision>
</link>
<joint name … type … > <origin… /> <parent link …/> <child link… /> <axis … /> <limit … /> </joint>
<plugin>
The parameters in Algorithm 2 obtained the aerodynamic force based on the airfoil geometry and aerodynamic curves, as shown in Figure 4. The aerodynamic parameters of each link are presented in detail in Table 2.
Algorithm 2: Aerodynamic Plugin Parameters.
<gazebo>
<plugin name=“aero_fly” filename=“libLiftDragPlugin.so”>
<air_density> … </air_density> # Density of the fluid this model is suspended in (kg/m3)
<cla> … </cla> #The ratio of the coefficient of lift and alpha slope before stall
<cla_stall> … </cla_stall> #The ratio of coefficient of lift and alpha slope after stall
<cda> … </cda> #The ratio of the coefficient of drag and alpha slope before stall
<cda_stall> … </cda_stall> #The ratio of coefficient of drag and alpha slope after stall
<alpha_stall> … </alpha_stall> #Angle of attack at stall point (radians)
<a0> … </a0> #initial angle of attack (radians)
<area> … </area> #Surface area of the link (m2)
<upward> x y z</upward> # 3-vector representing the upward direction of motion in the link frame
<forward> x y z</forward> # 3-vector representing the forward direction of motion in the link frame
<link_name>base_link</link_name> #Name of the link affected by the group of lift/drag properties
<cp> x y z</cp> #Center of pressure in link frame. The forces due to lift and drag will be applied here
</plugin>
</gazebo>
The lift and drag forces acting on the wings, propellers, ailerons, elevator, and rudder are depicted in Figure 8a by the red lines. For generating lift, draft, and moment, the propeller must rotate, whereas the wings, ailerons, elevators, and rudder require relative airspeed to generate aerodynamic forces. The top view of the chord line wing is shown in Figure 8b and Figure 8c displays the aerodynamic parameter curves for the wing, propeller, aileron, elevator, and rudder, which correspond to the airfoil’s shape specified by the NACA number. The NACA 2412 airfoil is used for the wings and propellers, whereas the NACA 0012 airfoil is employed for the ailerons, elevators, and rudder.

3. State Space Model and Controller

Considering the dynamic equations described in Section 2.1, a state space controller is proposed to maintain the stability of VFW and track the trajectory. Figure 9 depicts the block diagram of the LQR (linear quadratic regulator) controller based on state space representation. The diagram illustrates how system state is obtained from various sensors, including those for position, orientation, linear velocity, and angular velocity, resulting in twelve state data. This state is then compared to the desired state, and the comparison results are used as input for a full state feedback controller, which generates an input for the control signal (u). The control signal is calculated by multiplying the full-state feedback controller input by a negative feedback gain K. The control signal is first interpreted by the control element as the value of the motor rotational speed and servo position for the VFW actuator before being forwarded to the process plant. It iterates continuously as long as the VFW is carrying out its mission.
The LQR controller’s ability to provide feedback on the system’s entire state, rather than just partial feedback, is crucial. This ensures that the VFW moves precisely as desired. The output of the entire system (y) is in the form of six states that describe the VFW’s position and orientation. This level of detail enables close monitoring of the VFW’s movement, ensuring that it performs its task with high accuracy and reliability.
The VFW is characterized by twelve states, namely ( X , Y , Z , V x , V y , V z , ϕ , θ , ψ , ω ϕ , ω θ , and ω θ ) , representing the position, translation velocity, orientation, and angular velocity. Utilizing the parameters listed in Table 1, the state space fixed wing model is derived using Matlab to obtain matrices A, B, C, D, Q, R, and K. Equation (44) as follows:
x ˙ = A x + B u
describes the relationship between the system’s present state, input, and future state, and for VFW, detailed as
X ˙ Y ˙ Z ˙ V x ˙ V y ˙ V z ˙ ϕ ˙ θ ˙ ψ ˙ ω ϕ ˙ ω θ ˙ ω ψ ˙ = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 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 1 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 g 0 0 0 0 0 0 0 0 0 0 g 0 0 0 0 0 0 0 0 0 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 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 X Y Z V x V y V z ϕ θ ψ ω ϕ ω θ ω θ + 0 0 0 0 0 1 / m 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 / I x 0 0 0 0 0 0 0 0 0 0 0 0 1 / I y 0 0 0 0 0 0 0 0 0 0 0 0 1 / I z u 1 u 2 u 3 u 4
x ˙ = A x + B u
The relationship among the system’s state, input, and output id defined as
y = C x + D u
and elaborated in
y 1 y 2 y 3 y 4 y 5 y 6 = 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 X Y Z V x V y V z ϕ θ ψ ω ϕ ω θ ω ψ + [ 0 ]   u 1 u 2 u 3 u 4
y = C x + D u
The control signal
u = K x
is obtained by multiplying the negative feedback gain K with the state control. A, B, C, D, x, u, y, and K denote the state matrix, input matrix, output matrix, feed-forward matrix, current state, control signal, output, and feedback gain, respectively. Feedback gained K as follows:
K = 0.0000 0.0000 1.0000 0.0000 0.0000 1.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.4515 0.0000 0.0000 1.4515 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 5.7562 0.0000 0.0000 0.0000 0.0000 6.0383 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 1.1268 0.0000 0.0000 0.0000 0.0000 1.2939 0.0000 0.0000 0.0000 0.0000 1.0646
which depends on Q and R matrix. Q matrix defines the weights on the states, and the R matrix defines the weights on the control input in the cost function. The Q and R matrices are given as identity matrices below.
Q = 1 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 1 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 1 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 1 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 1 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 1 0 0 0 0 0 0 0 0 0 0 0 0 1 R = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1
The control signals are defined as
u 2 = K 2.7   ϕ + K 2.10   ω ϕ
u 3 = K 3.8   θ + K 3.11   ω θ
u 4 = K 4.9   ψ + K 4.12   ω ψ
where u 2 is the control signal for roll, u 3 is the control signal for pitch, and u 4 is the control signal for yaw. For K a . b in Equations (51)–(53), a is the number of rows, and b is the number of columns in K. The value of u is the multiplication of the negative feedback gain K and the current state, as in Equation (46). The signal control for the quadrotor’s motor R is described as
R f r o n t _ r i g h t = u 1 + u 2 d b + u 3 d b + u 4 k
R f r o n t _ l e f t = u 1 + u 2 ( d b ) + u 3 d b + u 4 ( k )
R r e a r _ l e f t = u 1 + u 2 ( d b ) + u 3 ( d b ) + u 4 k
R r e a r _ r i g h t = u 1 + u 2 d b + u 3 ( d b ) + u 4 ( k )
Meanwhile, the control signals for fixed-wing flight control surfaces E are described as
E a i l e r o n _ l e f t = u 2 j l x
E a i l e r o n _ r i g h t = u 2 j l x
E e l e v a t o r = u 3 j l y
E r u d d e r = u 4 j l z

4. Simulation

This section presents the simulation of VFW bird take of mode transition that applied ground effect and external disturbance. Section 4.1 presents the transition trajectory of bird take-off mode, Section 4.2 presents the simulation setup, and Section 4.3 discusses the simulation result.

4.1. Bird Take-Off Mode Transition Trajectory

The transition scenario in Figure 10 imitates how a bird takes off. VFW starts the transition from the ground and gains altitude and airspeed simultaneously. The quadrotor and fixed-wing propulsion turn on together. VFW flies up and gets forward speed with a specific rate of climb until it reaches the desired altitude. After reaching the desired altitude, VFW maintains its altitude by adjusting the rate of climb value to 0 m/s until it reaches the desired airspeed. When VFW exceeds the desired airspeed, the fixed wing can fly using the lift from the wing. So the quadrotor can be slowed down and turned off. It indicates that the transition is successful.

4.2. Simulation Setup

Several key components must be set up to carry out a realistic visual 3D simulation using ROS and the Gazebo simulator. The first step is to define the physical properties of the robot, which include inertial body properties, visual properties, collision body properties, and aerodynamic parameters. These parameters are described using URDF configuration, as outlined in Algorithm 1 and Algorithm 2 in Section 2.3. Additionally, the physical parameters of the VFW and environmental parameters, such as gravity and air density, are described in detail in Table 1. The full aerodynamic parameters of each link of the VFW are described in Table 2 and implemented using a modified LiftDragPlugin. These parameters play a crucial role in determining how closely the simulation will resemble the real world.
Once the physical properties have been defined, the VFW’s sensors, actuators, and controllers must be added to the robot model using ROS packages. The robot’s position, orientation, linear velocity, and angular velocity relative to the global inertial frame are represented by twelve states (X, Y, Z, V_x, V_y, V_z, ϕ, θ, ψ, ω_ϕ, ω_θ, and ω_ψ) acquired using rosservices GetModelState as a sensor. The actuators are represented by joint_state_publisher, which enables movement between connected links.
The controller in Gazebo consists of two types: a physical controller and a dynamic controller. The physical controller uses the joint controller package to set the maximum and minimum joint angles, joint velocity limits, and actuator response. On the other hand, dynamic controllers use rostopic to implement the LQR controller. These controllers play an essential role in governing the VFW’s behavior in the simulation.

4.3. Simulation Result

In this research, a visually realistic VFW simulation has been successfully conducted using Gazebo and ROS. This simulator integrates physical parameters, aerodynamics parameters, visual parameters, sensors, and actuator control. Figure 11a shows VFW visual 3D simulation in Gazebo, while Figure 11b the inertial of each link. Each link has a center of mass, as shown in Figure 11c, where the center of pressure for aerodynamic force is applied. Figure 11d shows the joint position of each child’s link with its parent, represented by three arrows (red, green, and blue).
To evaluate the performance of the proposed system, a fully autonomous transition was tested under three flight conditions. The first test was conducted under normal conditions, without applying ground effect and external disturbance. The second test was carried out in ground effect but without external disturbance. The third test included both ground effect and external disturbance. For a successful transition, predetermined parameters were set, including an altitude of more than 5 m and airspeeds of more than 5 m/s. Additionally, the attitude stability indicators were monitored, with steady-state error roll and pitch less than 5 degrees and yaw rate less than 3 degrees/s.
Figure 12 presents a real-time visualization of the simulated transition process, beginning with the activation of both quadrotor and fixed-wing propulsions in Figure 12a. Due to the low airspeed during this phase, the quadrotor is primarily responsible for controlling the vehicle’s flight attitude, while deflections of the ailerons, elevator, and rudder have negligible effects on the vehicle’s roll, pitch, and yaw moments. As shown in Figure 12b, the vehicle aims to attain the predetermined altitude and airspeed in the middle of the transition. Figure 12c displays the end of the transition, with the vehicle cruising at its designated altitude and minimum airspeed.
Figure 13 shows the flight test result in normal conditions. In this flight test, the ground effect was not applied without disturbance, as shown in Figure 13d. Figure 13a shows the flight phase; there are preparation flight, transition, and cruising. Figure 13b,e show attitude response with a maximum roll error value is 3.87° and a maximum pitch error value is 3.66°. The maximum yaw rate is 0.3 degree/s, as shown in Figure 13g. It indicates that the controller can maintain flight stability during the transition. The transition was completed at a predetermined altitude and airspeed in 4.71 s, as shown in Figure 13a. Compared to previous research [10,12], this uninterrupted test result has better control results seen from the maximum attitude error during the transition process.
The influence of the ground effect on VFW performance can be seen in the test results presented in Figure 14a. This effect occurs up to a height of 0.762 m ( h r = 6 r [34,35]) when t = 51 to t = 210, as shown in Figure 14d. The test results indicate that aircraft stability is more challenging to maintain in ground effect than in normal conditions. The maximum roll error is 4.3 when t = 187, and the maximum error pitch is 8.17 when t = 183 due to the momentary lack of lift when leaving the ground effect. Despite these challenges, the controller was able to maintain flight stability with a settling time roll of 0.45 s and a settling time pitch of 0.60 s. It is noteworthy that the yaw rate is not disturbed and still maintains stability, with a maximum value of 0.2 degree/s, as shown in Figure 14g.
Moreover, the ground effect provides the advantage of additional lift that increases the vertical speed of the aircraft. This effect enables the predetermined altitude to be reached more quickly, resulting in a faster transition time of 4.53 s, as shown in Figure 14c. In comparison to previous research on handling ground effect [41,42,43], the proposed controller in this study demonstrated a good attitude response. Despite the challenges posed by the ground effect, the results of the test suggest that with the right controller, it is possible to maintain flight stability and take advantage of the increased lift provided by this phenomenon.
The flight test results with external disturbances are presented in Figure 15. At time t = 358, external disturbances were applied at body coordinates 0.6, −0.3, and 0 with a force of 600 N to each of the x, y, and z axes, as shown in Figure 15d. Despite being disturbed, the controller was able to maintain stability with a settling time of 0.93 s for roll and 0.90 s for pitch, while the maximum roll and pitch errors were 9.43° and 9.40°, respectively. The maximum yaw rate was 1.1 degrees/s. The results indicate that the proposed controller can quickly handle external disturbances. The transition was completed in 5.15 s, as Figure 15a shows.
In short, the simulation results showed that the proposed controller was able to maintain flight stability during the transition of the proposed VFW model and trajectory, even in the presence of ground effect and external disturbance, which caused greater instability. The ground effect provided additional lift, which increased the vertical speed and allowed for a quicker transition to the predetermined altitude. Despite external disturbances, the proposed controller could handle them quickly and maintain stability during the transition. Overall, the study emphasized the importance of considering the ground effect in transitioning from vertical take-off and landing to fixed-wing flight and the effectiveness of the proposed controller in maintaining flight stability during the transition under various flight conditions.
However, it is important to note that there may be some differences between the simulation and the real-world application of the ground effect on VTOL. These differences could arise from various factors, such as variations in atmospheric conditions, surface properties, and other environmental factors that may not be accurately represented in the simulation. So it is necessary to conduct experiments to validate the influence of ground effect on bird take-off mode transition of VFW in the real world.

5. Conclusions

This paper concentrated on modeling, controlling, and conducting a realistic 3D simulation of the VFW transition considered ground effect. The mathematical and 3D models were obtained with fully dynamic and aerodynamic properties. Ground effect in VFW also had been successfully modeled as a considerable influence in transition flight. Furthermore, a realistic 3D simulation of the bird take-off mode transition was conducted on Gazebo and ROS to validate the proposed model, trajectory, and controller performance. Various flight tests have been performed to verify the robustness of the LQR controller in ground effect and external disturbance. The results show that the proposed model and controller carry out flight missions successfully with the bird take-off mode transition in any conditions. The controller has reliability and robustness in compensating for attitude errors caused by ground effects and external disturbance.
In the future, there are several areas of research that could build on the findings of this study. One promising direction is to optimize the flight performance of the system by implementing adaptive control techniques. Adaptive control can enhance the system’s ability to operate effectively even in the face of significant uncertainties, such as those caused by ground effects and external disturbances. Additionally, it would be valuable to conduct a bird takeoff mode transition experiment to further validate the effectiveness of the proposed approach in real-world scenarios.

Author Contributions

Conceptualization, E.I., A.H. and A.D.; data curation, E.I. and A.D.; formal analysis, E.I., A.H. and A.D.; methodology, E.I. and A.D.; resources, E.I. and A.D.; mathematical and 3D model, E.I.; software and simulation, E.I.; investigation, E.I.; supervision, A.H. and A.D.; writing—original draft, E.I.; writing-review and editing, E.I., A.H. and A.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jo, D.; Kwon, Y. Analysis of VTOL UAV Propellant Technology. J. Comput. Commun. 2017, 5, 76–82. [Google Scholar] [CrossRef]
  2. Govdeli, Y.; Tran, A.T.; Kayacan, E. Multiple Modeling and Fuzzy Switching Control of Fixed-Wing VTOL Tilt-Rotor UAV. In Fuzzy Techniques: Theory and Applications IFSA/NAFIPS, 5th ed.; Kearfott, R., Batyrshin, I., Reformat, M., Ceberio, M., Kreinovich, V., Eds.; Springer: Cham, Switzerland, 2019; Volume 1000, pp. 270–284. [Google Scholar]
  3. Cobleigh, B.R. Ikhana: A NASA UAS Supporting Long Duration Earth Science Missions; Technical Report NASA TM-2007-214614; NASA Dryden Flight Research Center: Edwards, CA, USA, 2007. [Google Scholar]
  4. Logan, M.; Chu, J.; Motter, M.; Carter, D.; Ol, M.; Zeune, C. Small UAV research and evolution in long endurance electric powered vehicles. In Proceedings of the AIAA Infotech at Aerospace Conference and Exhibit, Rohnert Park, NV, USA, 7–10 May 2007. [Google Scholar]
  5. Tielin, M.; Chuanguang, Y.; Wenbiao, G.; Zihan, X.; Qinling, Z.; Xiaoou, Z. Analysis of technical characteristics of fixed-wing VTOL UAV. In Proceedings of the IEEE International Conference on Unmanned Systems (ICUS), Beijing, China, 27–29 October 2017; pp. 293–297. [Google Scholar]
  6. Nugroho, G.; Zuliardiansyah, G.; Rasyiddin, A.A. Performance Analysis of Empennage Configurations on a Surveillance and Monitoring Mission of a VTOL-Plane UAV Using a Computational Fluid Dynamics Simulation. Aerospace 2022, 9, 208. [Google Scholar] [CrossRef]
  7. Gunarathna, J.K.; dan Munasinghe, R. Development of a Quad-rotor Fixed-wing Hybrid Unmanned Aerial Vehicle. In Proceedings of the Moratuwa Engineering Research Conference (MERCon), Moratuwa, Sri Lanka, 30 May–1 June 2018; pp. 72–77. [Google Scholar]
  8. Saeed, A.S.; Younes, A.B.; Cai, C.; Cai, G. A survey of hybrid Unmanned Aerial Vehicles. Prog. Aerosp. Sci. 2018, 98, 91–105. [Google Scholar] [CrossRef]
  9. Zaludin, Z.; Gires, E. Automatic Flight Control Requirements for Transition Flight Phases When Converting Long Endurance Fixed Wing UAV to VTOL Aircraft. In Proceedings of the 2019 IEEE International Conference on Automatic Control and Intelligent Systems (I2CACIS), Selangor, Malaysia, 29 June 2019; pp. 273–278. [Google Scholar]
  10. Gu, H.; Lyu, X.; Li, Z.; Shen, S.; Zhang, F. Development and experimental verification of a hybrid vertical take-off and landing (VTOL) unmanned aerial vehicle(UAV). In Proceedings of the 2017 International Conference on Unmanned Aircraft Systems (ICUAS), Miami, FL, USA, 13–16 June 2017; pp. 160–169. [Google Scholar]
  11. Dundar, O.; Bilici, M.; Unler, T. Design and performance analyses of a fixed wing battery VTOL UAV. Int. J. Eng. Sci. Technol. 2020, 23, 1182–1193. [Google Scholar] [CrossRef]
  12. Gunarathna, J.K.; Munasinghe, R. Simultaneous Execution of Quad and Plane Flight Modes For Efficient Take-Off of Quad-Plane Unmanned Aerial Vehicles. Appl. Sci. 2021. accepted. [Google Scholar] [CrossRef]
  13. Kamal, A.M.; Ramirez, S.A. Design methodology for hybrid (VTOL + Fixed Wing) unmanned aerial vehicles. Aeronaut. Aerosp. Open Access J. 2018, 2, 165–176. [Google Scholar]
  14. Curtiss, J.H.C.; Sun, M.; Putman, W.F.; Hanker, J.E.J. Rotor aerodynamics in ground effect at low advance ratios. J. Am. Helicopter Soc. 1984, 29, 48–55. [Google Scholar] [CrossRef]
  15. Lee, T.E.; Leishman, J.G.; Ramasamy, M. Fluid dynamics of interacting blade tip vortices with a ground plane. J. Am. Helicopter Soc. 2010, 55, 2200500–2200516. [Google Scholar] [CrossRef]
  16. Tanner, P.; Overmeyer, A.; Jenkins, L.; Yao, C.S.; Bartram, S. Experimental investigation of rotorcraft outwash in ground effect. In Proceedings of the 71th Annual Forum of the American Helicopter Society, Virginia Beach, Virginia, 5–7 May 2015; pp. 1–26. [Google Scholar]
  17. Cuevas, P.S.; Heredia, G.; Ollero, A. Characterization of the Aerodynamic Ground Effect and Its Influence in Multirotor Control. Int. J. Aerosp. Eng. 2017, 2017, 1–17. [Google Scholar] [CrossRef]
  18. Bernard, D.D.C.; Riccardi, F.; Giurato, M.; Lovera, M. A dynamic analysis of ground effect for a quadrotor platform. Int. Fed. Autom. Control. 2017, 50, 10311–10316. [Google Scholar] [CrossRef]
  19. Alaez, D.; Olaz, X.; Prieto, M.; Porcellinis, P.; Villadangos, J. HIL Flight Simulator for VTOL-UAV Pilot Training Using X-Plane. Information 2023, 13, 585. [Google Scholar] [CrossRef]
  20. Deng, Z.; Guo, Z.; Wu, L.; You, Y. Trajectory Planning for Emergency Landing of VTOL Fixed-Wing Unmanned Aerial Vehicles. Mob. Inf. Syst. 2021, 2021, 6289822. [Google Scholar] [CrossRef]
  21. Wang, W.; Bai, X.; Zheng, L. PD control of VTOL aircraft trajectory tracking based on double-loop design. J. Phys. Conf. Ser. 2022, 2358, 012008. [Google Scholar] [CrossRef]
  22. Dharmawan, A.; Ashari, A.; Aprilia, A.G.; Handayani, A.M. Auto VTOL System on Quadrotor Using Madgwick Quaternion Kalman Filter and LQR. In Proceedings of the 2018 4th International Conference on Science and Technology (ICST), Yogyakarta, Indonesia, 7–8 August 2018; pp. 1–6. [Google Scholar]
  23. Pan, K.; Chen, Y.; Wang, Z.; Wu, H.; Cheng, L. Quadrotor Control based on Self-Tuning LQR. In Proceedings of the 2018 37th Chinese Control Conference (CCC), Wuhan, China, 25–27 July 2019; pp. 9974–9979. [Google Scholar]
  24. Shehzad, M.F.; Bilal, A.; Ahmad, H. Position & Attitude Control of an Aerial Robot (Quadrotor) With Intelligent PID and State feedback LQR Controller: A Comparative Approach. In Proceedings of the 2019 16th International Bhurban Conference on Applied Sciences and Technology (IBCAST), Islamabad, Pakistan, 8–12 January 2019; pp. 340–346. [Google Scholar]
  25. Kumar, K.K.S.; Arya, H.; Joshi, A. Automatic Control of Harrier Maneuver of an Agile Fixed-Wing UAV. In Proceedings of the 2019 Fifth Indian Control Conference (ICC), New Delhi, India, 9–11 January 2019; pp. 136–141. [Google Scholar]
  26. Jeevan, H.L.; Narahari, H.K.; Sriram, A.T. Development of pitch control subsystem of autopilot for a fixed wing unmanned aerial vehicle. In Proceedings of the 2018 2nd International Conference on Inventive Systems and Control (ICISC), Coimbatore, India, 9–10 January 2018; pp. 1233–1238. [Google Scholar]
  27. Dharmawan, A.; Putra, A.E.; Tresnayana, I.M.; Wicaksono, W.A. The Obstacle Avoidance System In A Fixed-Wing UAV When Flying Low Using LQR Method. In Proceedings of the 2019 International Conference on Computer Engineering, Network, and Intelligent Multimedia (CENIM), Surabaya, Indonesia, 19–20 November 2019; pp. 1–7. [Google Scholar]
  28. Nelissen, L.P.M. Development of Gazebo-Based High-Fidelity Simulation Environment for Morphing-Wing Uavs. Bachelor’s Thesis, University of Twente, Enschede, The Netherlands, December 2021. [Google Scholar]
  29. Alaez, D.; Olaz, X.; Prieto, M.; Villadangos, J.; Astrain, J.J. VTOL UAV digital twin for take-off, hovering and landing in different wind conditions. Simul. Model. Pract. Theory 2022, 123, 102703. [Google Scholar] [CrossRef]
  30. Meyer, J.; Sendobry, A.; Kohlbrecher, S.; Klingauf, U.; von Stryk, O. Comprehensive Simulation of Quadrotor UAVs Using ROS and Gazebo. In Proceedings of the Third International Simulation, Modeling, and Programming for Autonomous Robots, Tsukuba, Japan, 5–8 November 2012; pp. 400–411. [Google Scholar]
  31. Sciortino, C.; Fagiolini, A. ROS/Gazebo-Based Simulation of Quadcopter Aircrafts. In Proceedings of the 2018 IEEE 4th International Forum on Research and Technology for Society and Industry (RTSI), Palermo, Italy, 10–13 September 2018; pp. 1–6. [Google Scholar]
  32. Ale Isaac, M.S.; Luna, M.A.; Ragab, A.R.; Ale Eshagh Khoeini, M.M.; Kalra, R.; Campoy, P.; Flores Peña, P.; Molina, M. Medium-Scale UAVs: A Practical Control System Considering Aerodynamics Analysis. Drones 2022, 6, 244. [Google Scholar] [CrossRef]
  33. Danjun, L.; Yan, Z.; Zongying, S.; Geng, L. Autonomous landing of quadrotor based on ground effect modelling. In Proceedings of the 2015 34th Chinese Control Conference (CCC), Hangzhou, China, 28–30 July 2015; pp. 5647–5652. [Google Scholar]
  34. Davis, E.; Spollard, J.; Pounds, P. Passive height stability and trajectory repeatability of a quadrotor maneuvering in ground effect with regulated voltage bus. In Proceedings of the Australasian Conference on Robotics and Automation (ACRA 2015), Canberra, Australia, 2–4 December 2015. [Google Scholar]
  35. Sharf, I.; Nahon, M.; Harmat, A.; Khan, W.; Michini, M.; Speal, N.; Trentini, M.; Tsadok, T.; Wang, T. Ground effect experiments and model validation with Draganflyer X8 rotorcraft. In Proceedings of the 2014 International Conference on Unmanned Aircraft Systems (ICUAS), Orlando, FL, USA, 27–30 May 2014; pp. 1158–1166. [Google Scholar]
  36. Jo, D.; Kwon, Y. Development of Autonomous VTOL UAV for Wide Area Surveillance. World J. Eng. Technol. 2019, 7, 227–239. [Google Scholar] [CrossRef]
  37. 8083-1A; Aircraft Weight and Balance Handbook. U.S Department of Transportation, Federal Aviation Administration: Washington, DC, USA, 2007.
  38. Azman, S. Design the VTOL Aircraft for Land Surveying Purposes. Mechanical Engineering Study Program; Faculty of Mechanical Engineering, University Teknologi Malaysia: Skudai, Malaysia, 2018. [Google Scholar]
  39. Saengphet, W.; Thumthae, C. Conceptual Design of Fixed Wing-VTOL UAV for AED Transport. In Proceedings of the 7th TSME International Conference on Mechanical Engineering, Chiang Mai, Thailand, 13–16 December 2016. [Google Scholar]
  40. Anderson, J.D. Aircraft Performance and Design; Mc-Graw-Hill Education: New York, NY, USA, 1998. [Google Scholar]
  41. Wu, Z.; Wu, T.W. Influence of downwash flow on helicopter airborne missile initial trajectory. J. Proj. Rocket. Missiles Guid. 2009, 29, 202–204. [Google Scholar]
  42. Qin, Y.H.; Zhu, Q.H.; Shao, S. Aerodynamic characteristics analysis for hovering coaxial rotors in ground effect. J. Nanjing Univ. Aeronaut. Astronaut. 2015, 47, 266–274. [Google Scholar]
  43. Zheng, Y.; Yang, S.; Liu, X.; Wang, J.; Norton, T.; Chen, J.; Tan, Y. The computational fluid dynamic modeling of downwash flow field for a six-rotor UAV. Front. Agric. Sci. Eng. 2018, 5, 159–167. [Google Scholar] [CrossRef]
Figure 1. VFW transition strategies.
Figure 1. VFW transition strategies.
Drones 07 00330 g001
Figure 2. Bird take-off mode transition.
Figure 2. Bird take-off mode transition.
Drones 07 00330 g002
Figure 3. Configuration of VFW.
Figure 3. Configuration of VFW.
Drones 07 00330 g003
Figure 4. Thrust required versus height.
Figure 4. Thrust required versus height.
Drones 07 00330 g004
Figure 5. Rotor downwash flow field, (a) out-of-ground effect, (b) in-ground effect.
Figure 5. Rotor downwash flow field, (a) out-of-ground effect, (b) in-ground effect.
Drones 07 00330 g005
Figure 6. VFW dimension (a) isometric front view, (b) isometric top view, (c) isometric side view.
Figure 6. VFW dimension (a) isometric front view, (b) isometric top view, (c) isometric side view.
Drones 07 00330 g006
Figure 7. Links of the VFW.
Figure 7. Links of the VFW.
Drones 07 00330 g007
Figure 8. (a) Aerodynamics properties of airfoil, (b) top of view chord line wing, (c) coefficient lift and drag versus angle of attack.
Figure 8. (a) Aerodynamics properties of airfoil, (b) top of view chord line wing, (c) coefficient lift and drag versus angle of attack.
Drones 07 00330 g008
Figure 9. LQR controller block diagram.
Figure 9. LQR controller block diagram.
Drones 07 00330 g009
Figure 10. Flowchart bird take-off mode transition.
Figure 10. Flowchart bird take-off mode transition.
Drones 07 00330 g010
Figure 11. VTOL fixed-wing visual 3D simulation in Gazebo, (a) simulation environment, (b) inertial of each link, (c) center of mass of each link, (d) joint’s position of each link.
Figure 11. VTOL fixed-wing visual 3D simulation in Gazebo, (a) simulation environment, (b) inertial of each link, (c) center of mass of each link, (d) joint’s position of each link.
Drones 07 00330 g011
Figure 12. Transition mode in Gazebo simulator, (a) beginning of transition, (b) middle of transition, (c) end of transition.
Figure 12. Transition mode in Gazebo simulator, (a) beginning of transition, (b) middle of transition, (c) end of transition.
Drones 07 00330 g012
Figure 13. Flight test result when ground effect is not applied and without disturbance; (a) flight phase versus time; (b) roll versus time; (c) altitude versus time; (d) flight condition versus time; (e) pitch versus time; (f) airspeed versus time; (g) yaw rate versus time.
Figure 13. Flight test result when ground effect is not applied and without disturbance; (a) flight phase versus time; (b) roll versus time; (c) altitude versus time; (d) flight condition versus time; (e) pitch versus time; (f) airspeed versus time; (g) yaw rate versus time.
Drones 07 00330 g013aDrones 07 00330 g013b
Figure 14. Flight test result when ground effect is applied and without disturbance; (a) flight phase versus time; (b) roll versus time; (c) altitude versus time; (d) flight condition versus time; (e) pitch versus time; (f) airspeed versus time; (g) yaw rate versus time.
Figure 14. Flight test result when ground effect is applied and without disturbance; (a) flight phase versus time; (b) roll versus time; (c) altitude versus time; (d) flight condition versus time; (e) pitch versus time; (f) airspeed versus time; (g) yaw rate versus time.
Drones 07 00330 g014aDrones 07 00330 g014b
Figure 15. Flight test result when ground effect is applied and with disturbance; (a) flight phase versus time; (b) roll versus time; (c) altitude versus time; (d) flight condition versus time; (e) pitch versus time; (f) airspeed versus time; (g) yaw rate versus time.
Figure 15. Flight test result when ground effect is applied and with disturbance; (a) flight phase versus time; (b) roll versus time; (c) altitude versus time; (d) flight condition versus time; (e) pitch versus time; (f) airspeed versus time; (g) yaw rate versus time.
Drones 07 00330 g015
Table 1. VFW parameters.
Table 1. VFW parameters.
ParametersValue
wing span1.470 m
aircraft length0.924 m
mass1.9835 kg
airfoil wingNACA 2412
airfoil propellerNACA 2412
airfoil aileron, elevator dan rudderNACA 0012
initial AoA fixed wing
fixed wing and quadrotor propeller diameter0.254 m (10 inches)
area formed propeller quadrotor rotate5.06 × 10−2 m2
distance center of mass and motor quadrotor0.38891 m
center of massX = 0.34, Y = 0.00, Z = 0.036
gravity9.81 kg/m2
inertia about x-axis (Ix)0.023426 kg/m2
inertia about y-axis (Iy)0.044475 kg/m2
inertia about z-axis (Iz)0.066685 kg/m2
air density ( ρ )1.225 kg/m3
coefficient torque (CT)1.3602 × 10−1
coefficient thrust (CQ)5.5555 × 10−2
Table 2. Aerodynamic parameters of each link.
Table 2. Aerodynamic parameters of each link.
Aerodynamics ParameterBase_LinkLink_1, Link_2, Link_3, Link_4, Link_5Link_6, Link 7Link_8Link_9
air_density1.2251.2251.2251.2251.225
cla4.754.755.95.95.9
cla_stall−3.85−3.85−3.85−3.85−3.85
cda−0.4−0.4−0.4−0.4−0.4
cda_stall−0.92−0.92−0.92−0.92−0.92
alpha_stall0.3391.50.3390.3390.339
a00.050.40.00.00.0
area0.20.020.040.050.04
upward0 0 10 0 1 (link_1, link_2, link_3, link_4)
1 0 0 (link_5)
0 0 10 0 10 1 0
forward−1 0 00 −1 0 (link_1, link_4)
0 1 0 (link_2, link_3, link_5)
−1 0 0−1 0 0−1 0 0
cp0.340 0.000 0.0300.065 0.275 0.036 (link_1)
0.065 −0.275 0.036 (link_2)
0.615 −0.275 0.036 (link_3)
0.615 0.275 0.036 (link_4)
0 0 0 (link_5)
0.394 0.312 0.046 (link_6)
0.394 −0.312 0.046 (link_7)
0.900 0.000 0.0151.571 0.000 0.010
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

Irmawan, E.; Harjoko, A.; Dharmawan, A. Model, Control, and Realistic Visual 3D Simulation of VTOL Fixed-Wing Transition Flight Considering Ground Effect. Drones 2023, 7, 330. https://doi.org/10.3390/drones7050330

AMA Style

Irmawan E, Harjoko A, Dharmawan A. Model, Control, and Realistic Visual 3D Simulation of VTOL Fixed-Wing Transition Flight Considering Ground Effect. Drones. 2023; 7(5):330. https://doi.org/10.3390/drones7050330

Chicago/Turabian Style

Irmawan, Erwhin, Agus Harjoko, and Andi Dharmawan. 2023. "Model, Control, and Realistic Visual 3D Simulation of VTOL Fixed-Wing Transition Flight Considering Ground Effect" Drones 7, no. 5: 330. https://doi.org/10.3390/drones7050330

APA Style

Irmawan, E., Harjoko, A., & Dharmawan, A. (2023). Model, Control, and Realistic Visual 3D Simulation of VTOL Fixed-Wing Transition Flight Considering Ground Effect. Drones, 7(5), 330. https://doi.org/10.3390/drones7050330

Article Metrics

Back to TopTop