Design and Experimental Study of a Turbojet VTOL Aircraft with One-Dimensional Thrust Vectoring Nozzles

: Turbojet engines have a slow response and require mechanical assistance to perform vertical take-offs and landings. In this paper, we propose a vertical take-off and landing (VTOL) aircraft that contains one-dimensional thrust vectoring nozzle modules to solve the problem of dependence on the response of the turbojet. The turbojet engine is placed horizontally in the module, the rotational speeds are ﬁxed, and the force along the turbojet axis remains constant as the nozzle rotates from the horizontal to the vertical position. The aircraft is composed of four modules, whose position and attitude control are realized by using the thrust difference caused by the different angles of the four vector nozzles. The modules no longer rely on the response speed of the turbojet engine, and controlling the angles of four nozzles makes it easier to operate the aircraft. The horizontal placement of the turbojet can reduce the thickness of the aircraft fuselage and decrease air resistance. The process of moving the nozzle from a horizontal to a vertical position prevents the engine from blowing directly on the ground and, thus, reduces the ground effect.


Introduction
In the past, turbojet vertical take-off and landing (VTOL) aircraft were mainly used in the military, and their advantages of a large thrust-to-weight ratio and high fuel efficiency also have broad application prospects in emergency transportation and rescue and other fields. VTOL aircraft have simple requirements for take-off and landing sites, can provide precise hovering and good maneuverability, and are very flexible [1]. They can not only independently complete intelligent inspections [2], but they can also protect agricultural developments [3], support logistics and distribution [4], and complete other tasks. Additionally, VTOLs can provide mobile platforms for cross-domain operations, such as airborne mobile arms operations [5] and Leonardo [6], which have particularly high requirements for load capacity and mobility.
Common VTOL power units include electric propellers and turbojets [7]. The authors of [8] summarized the detailed design of specific drone elements and propulsion components. There are many platforms powered by electric propellers, such as ducted [9], dual-rotor, tri-rotor, quad-rotor, and vector quad-rotor propulsion systems. Their overall characteristics are ease of use, fast response speed, and high stability of the motor. The disadvantage is that the propeller is bulky when there is a certain load, and the battery is located in the fuselage, where it is a fixed structural weight. The energy density of the battery is less than that of fuel oil, and the specific energy of the lithium-ion battery is approximately 120~200 w·h/kg, which is much less than that of fuel oil, i.e., 12 kw·h/kg [10]. In contrast, the advantages of the turbojet VTOL aircraft are the large thrust-to-weight ratio and compact structure. Under the same load conditions, the volume is smaller than that of propeller-type platforms. At present, the turbojet VTOL aircraft that have been applied in the military field mainly include the British Harrier, AV-8A of the US Marine Corps, Jacques-38 and Jacques-141 developed by the former Soviet Union, and F-35B of the US [11].
The mathematical modeling of the Jetcat P200SX micro turbojet engine was reported in the literature [12]. Turbojet aircraft usually use vector thrusting rather than the differential speed of the turbojet engine to achieve vertical take-off and landing, as the speed response process of the turbojet engine is much slower than that of motors. There are three ways to achieve thrust vectoring: to tilt the engine, to use the lift engine and main engine at the same time, and to deflect the nozzle of the engine. In [13,14], a vector quad-rotor was proposed in which the engine was rotated along two axes to achieve zero-pitch-angle flight of the aircraft in a redundant manner. However, the engine weight was large, which imposed higher requirements on the tilting mechanism. Moreover, tilting the engine disrupted the streamlines of the fuselage and increased the resistance during cruising. The Russian Jacques-38 and Jacques-141 are VTOL aircraft that use a combination of lift engines and deflected main engine nozzles [15]. In this layout, the lift engine occupies a large space in the fuselage, and the engine is not used during cruise flight. Vectoring nozzles are the key to engine thrust vectoring technology, which can, in principle, be divided into mechanical regulation, plasma jet regulation, and pneumatic jet regulation. Among them, the mechanical regulation of the vectoring nozzle is the most mature and widely used. A novel large-deflection-angle thrust vectoring nozzle called ACHEON was adopted, and a complete self-position control scheme for this type of fixed-wing UAV with thrust vectoring was developed in [16]. In [17], the thrust vectoring nozzle was controlled by the steering gear, which could deflect the pitching and yawing directions. The installation positions of the connecting rods of the two channels were orthogonal, which eliminated the coupling motion between the pitching and yawing channels. In addition, flaps with an airfoil were placed in the nozzle outlet of the ducted fan, and the attitude control was generated by changing the angle of attack of the flaps in [18]. However, the nozzle rotation angle was limited, and when the engine was perpendicular to the surface of the fuselage, it increased the difficulty of the design of the streamline of the VTOL aircraft.
In terms of attitude assistance control, the F35B adopts a vertical thrust system, which can generate lift force during takeoff. It also leads the compression of air generated by the engine to the wings on both sides, and it is possible to adjust the direction and volume of the jet exhaust to control the attitude of the aircraft during take-off or landing. Additional electric ducted fans were introduced to assist in the attitude control in [19]. A Pegasus-type thrust steering engine was installed in the middle of the Harrier fuselage [20]. In the front and rear, two pairs of rotatable nozzles were located on both sides of the fuselage. They could be rotated from 0 • to 98 • [21] and were symmetrical with respect to the center of gravity of the fuselage. As the four nozzles must be linked, and the resultant force of the thrust generated by the four nozzles is always required to pass through the center of the gravity of the aircraft, the Harrier was equipped with three pairs of attitude control nozzles. According to [13], a fixed engine thrust was selected, and attitude control was performed by using a two-dimensional vector method. The above methods made the structure of the aircraft more complicated.
In order to solve the problem of the slow speed response of turbojet engines, we designed one-dimensional thrust vectoring nozzle modules that arranged the turbojet engine horizontally and rotated the nozzle from horizontal to vertical positions. In the process of downward rotation of the nozzle, the direct blowing of the turbojet engine on the ground was prevented, which reduced the ground effect. Referring to the design of the guide vane in two references [22,23], the nozzle with an arc-shaped guide vane had a more uniform airflow streamline and better thrust efficiency. The thrust distribution characteristics of one-dimensional vector nozzle modules and the thrust change process during the nozzle rotation were obtained. The thrust of the engine was fixed, and the thrust difference caused by the difference between the angles of multiple vector nozzles was used for attitude stability control, which reduced the difficulty of control.
The rest of this paper is organized as follows. The structure of the one-dimensional thrust vectoring nozzle modules is introduced in Section 2. The computational fluid dynamics (CFD) simulations and experimental tests conducted in this study are described in Section 3. The control strategy of a one-dimensional thrust vectoring nozzle turbojet aircraft and the results of the simulations are discussed in Section 4. Section 5 presents the experimental results, and conclusions are drawn in Section 6.

Hardware
As shown in Figure 1, the mechanical design of our aircraft consisted of two main parts: (i) a central rigid body that carried the fuel tanks, batteries, sensors, and control systems required for flight and (ii) four one-dimensional thrust vectoring nozzle modules. The choice of four modules represented a compromise between simplicity and the ability to generate enough thrust to compensate for gravity, with an associated mechanical complexity that could result in a bulky drone with limited flight time and payload. Each one-dimensional thrust vectoring nozzle module included a turbojet engine, a rotating pair of modules composed of nozzles and bearings, and a four-bar linkage that drove by using the steering gear. As shown in Figure 2, the nozzle was fixed to the rear of the turbojet through structural parts such as bearings, the axis of the turbojet engine coincided with the axis of the bearing, and the airflow generated by the turbojet was turned through the nozzle. The crank-rocker mechanism, which comprised crank 1 and rocker 2, could realize one-dimensional rotation with a steering-gear driving nozzle. A rotation range of no more than 180 • could be achieved by adjusting the length of the four-link rod. To avoid a singularity of the four-bar linkage mechanism, the rotation range here was selected as 90 • .
To improve the load capacity of the VTOL aircraft, an engine with a high thrust-toweight ratio was needed. The SWIWIN SW400Pro turbojet engine was selected here; it can provide a maximum thrust of 400 N. To improve the stability of the vector system, the KST HS25-30-M-5545-1 A steering gear was selected. Table 1 shows the complete list of electronic components included in this platform.

Study of Nozzle Characteristics
The presence of the nozzle changed the direction of the high-speed airflow generated by the turbojet. The airflow along the axial direction turned along the radial direction of the engine. Due to the obstruction and wear of the nozzle, its force efficiency on the airflow decreased. The force characteristics of the one-dimensional thrust vectoring nozzle module were related to the stability of the aircraft. The thrust distribution during the nozzle rotation was essential for modeling the aircraft. As shown in Figure 3, the diameter of nozzle D, the length of the straight pipe section L, the end angle α, and the presence or absence of guide vanes all affected the high-speed airflow generated by the turbojet engine. The arc-shaped guide vane was at 95 • with a thickness of 3 mm. To determine the impact of the nozzle, a CFD simulation and experimental analysis were conducted.

CFD Simulation Analysis
We focused on the viscous compressible turbulent motion of the high-speed airflow generated by the turbojet in the nozzle; thus, the gravitational effect, thermal radiation, chemical reaction, and the influence of the airflow before the nozzle of the turbojet were neglected in the computational domain. The Reynolds-averaged Navier-Stokes (RANS) formulations with the corresponding turbulence closures are the preferred methodsparticularly in industry-for the analysis of compressible turbulent jet flows. Within this framework, two-equation turbulence models are used the most. SST k-ω had a major improvement in performance over both the original k-ω and the standard k-ε model [24]. The turbulent viscosity, which was modified to account for the transport of the principal turbulent shear stress, gave it an advantage in terms of performance over both the standard k-ω and k-ε formulations. It behaved as a k-ω model in the inner region of the boundary layers and as a k-ε model in the outer region of the boundary layers, in the free-stream zones. The authors of [25] proved the capabilities of the SST k-ω turbulence model within the steady-state RANS formulation for predicting a high-subsonic, coaxial-heated jet issued from a bypass ratio of a five-turbofan engine model, which was consistent with our cases.
The structured mesh was divided by using Fluent meshing. The far-field of the nozzle outlet was cylindrical, the distance from the axis of the nozzle outlet was 40 times the nozzle radius, and the nozzle extended 20 times the distance of the radius in the circumferential direction. After encrypting the grids at the nozzle wall, the grid independence check was performed. For the same example, we divided the number of grids into 8.27 × 10 6 and 1.52 × 10 7 , respectively, and the calculated results were consistent. The calculated results are shown in Table 2. As the nozzles calculated in the simulations had different dimensions, the grid numbers differ from 7 × 10 6 to 9 × 10 6 . Referring to the settings of the methods in [25], the compressible flow simulations were carried out by using a pressure-based solver. For the RANS simulations, the pressurevelocity coupling was achieved through the Coupled scheme, which provided superior performance compared with the segregated scheme for steady-state flows. Second-order upwind finite-volume schemes were employed for the spatial discretization of the terms in the RANS equations. The specific values of the model boundary conditions were set as follows: The total pressure of the nozzle inlet was 1053 K, the total pressure was 1.446 × 10 5 Pa, the back pressure of the nozzle outlet interface was 1.0 × 10 5 Pa, the external far-field static pressure was 1.0 × 10 5 Pa, and the total temperature was fixed at 300 K.
The thrust calculation method was F = mv + A(P n − P 0 ), where m is the gas mass flow rate, v is the outlet gas velocity, A is the outlet area, P n is the outlet static pressure, and P 0 is the outside atmospheric pressure. According to the simulation results, the turbojet thrust was 362.743 N. Figure 4 shows the effects of different nozzle modules on the thrust. As shown in Figure 4a, when L gradually increased to 140 mm, F y increased to the maximum value of 317.183 N, and then the thrust decreased, with L increasing gradually. Figure 4b shows that F y gradually decreased and F z gradually increased as the angle α increased. Figure 4c shows that the thrusts in the Y and Z directions improved due to the guide vane. Figure 4d shows that the presence of the guide vane made the airflow distribution more uniform. In summary, the nozzle with an arc-shaped guide vane had the best thrust efficiency with the appropriate straight pipe length L and an angle equal to zero.

Nozzle Test Analysis
To validate the simulation results, different nozzles were printed and tested. From the simulation results, we were able to find that increasing the pipe length L improved the thrust, but also increased the aircraft volume and weight. As a comparison, we set L to 60 and 80 mm, α to 0 • and 10 • , and the presence or absence of an arc-shaped guide vane. As shown in Figure 5, we printed nozzles in four sizes, and their detailed sizes are shown in Table 3. We built an experimental platform based on a six-dimensional force balance, as shown in Figure 6. Full-throttle thrust tests were conducted in the same environment (0°C, 50 m above sea level), and the force test results are shown in Table 3.   In the comparison with the simulation results, the experiments and simulations showed the same trends when the structural parameters of the different nozzles changed. Nozzle IV, which had an arc-shaped guide vane, L of 80 mm, and α of 0 • , had a higher efficiency (87.56%) compared with that of the turbojet engine. Thus, nozzle IV was used in the aircraft.
In order to confirm how the thrust changed during the rotation of the nozzle, the nozzle was driven by a steering gear, as shown in Figure 6, and changes in the force were recorded. The experimental results are shown in Figure 7. The forces along the axial direction of the turbojet (F y ) remained at 125 N during the rotation of the nozzle. F x and F z , which were distributed in the nozzle rotation plane, were a trigonometric function of the nozzle angle. The total thrust through the nozzle remained 340 N.

Attitude Control Strategy
After obtaining the force variation characteristics of the one-dimensional thrust vectoring nozzle module, we analyzed how the position and attitude of the aircraft could be controlled by using the force and moment generated by the modules; this analysis is described in this section. We kept the turbojet at the same throttle and used the rotation of the nozzle to control the thrust because of the slow speed response of the turbojet. As shown in Section 3, although there were forces along the axial direction of the one-dimensional thrust vectoring nozzle module, the modules were placed symmetrically in the aircraft, and the forces along the Y-axis of the aircraft offset each other. Then, the thrust was divided along the XZ-plane, where F i1 (i = 1, 2, 3, 4) was along the Z-axis and F i2 (i = 1, 2, 3, 4) was along the X-axis of the body. The way of using the force distribution to realize attitude control is shown in Figure 8.
(1) Figure 8. Attitude control schematic: (a) nozzles motion to produce pitch angular acceleration; (b) nozzles motion to produce roll angular acceleration; (c) nozzles motion to produce yaw angular acceleration. Figure 8a shows how the aircraft realized pitch. When δ 1 = δ 2 > δ 3 = δ 4 , the force followed Equation (2) and generated positive torque around the Y-axis of the body, and the aircraft had positive pitch acceleration. In contrast, when δ 1 = δ 2 < δ 3 = δ 4 , it generated negative torque, and the aircraft had negative pitch acceleration. Notably, as the thrust was fixed, (2) Figure 8b shows how the aircraft realized roll. When δ 1 = δ 4 > δ 2 = δ 3 , the force followed Equation (3) and generated positive torque around the X-axis of the body; the aircraft had positive roll acceleration. In contrast, when δ 1 = δ 4 < δ 2 = δ 3 , it generated negative torque, and the aircraft had negative roll acceleration. Figure 8c shows how the aircraft realized yaw. When δ 2 = δ 4 > δ 1 = δ 3 , the force followed Equation (4) and generated positive torque around the Z-axis of the body; the aircraft had positive yaw acceleration. In contrast, when δ 2 = δ 4 < δ 1 = δ 3 , it generated negative torque, and the aircraft had negative yaw acceleration.

Modeling of the Aircraft
The following assumptions were introduced to simplify the airframe model: It was considered that the aircraft was a rigid symmetric structure and that no structural deformation occurred during flight; the airflows of the turbojets did not interfere with each other, and the thrust values of the four turbojets were identical when they had same throttle; changes in mass and moment of inertia were ignored. The mass change was caused by fuel consumption, which was uniform and continuous. As the fuel tank was designed to be in the center of the aircraft, near the center of mass, the moment of inertia changed slowly. As the aircraft had sufficient control force and torque, we used the initial mass and moment of inertia.
The inertial frame E was fixed on the ground, and its z-axis pointed upward. In contrast, the z-axis of the body frame B pointed downward, as shown in Figure 9; its origin was at the center of gravity. The symbols are listed in Table 4. R E B defined the transformation matrix that allowed the transformation of a vector from the body-fixed frame to the inertial frame, where s(·) and c(·) are the sin(·) and cos(·) of a given angle, respectively. Figure 9. The inertial and body coordinate systems. η = [φ,θ,ψ] T represents the attitude angular velocity coordinates in the Earth coordinate system, and ω B = [p, q, r] T represents the corresponding attitude angular velocity coordinates in the body coordinate system. According to the coordinate transformation matrix R E B , the relationship betweenη and ω B is described as follows, where s(·) and c(·) are the sin(·) and cos(·) of a given angle, respectively:

Symbol Meanings
attitude under the inertial system ω B = [p, q, r] T angular velocity in the body-fixed system I = diag[I xx , I yy , I zz ] moment of inertia along the body axis 1, 2, 3, 4) rotation angle of the i-th nozzle l 1 , l 2 distance between the thrust center and centroid m aircraft mass In Section 3, the forces along the axial direction of the turbojet remained unchanged during the rotation of the nozzle. As the one-dimensional thrust vectoring nozzle modules were placed symmetrically on the aircraft, the resultant force along the aircraft's Y-axis was offset to zero. The force produced by the one-dimensional thrust vectoring nozzle module was distributed in the X-and Z-directions relative to the aircraft; its resultant force was T. For the modeling of the body dynamics, the Newton-Euler formalism [26] was used. F and τ b represented the thrust and moment separately generated by the aircraft, F g represented the gravity, F d represented the wind resistance force, and τ d represented the resistance torque. The general form is the following: The forces and moments produced by the aircraft are the following: T is a positive constant matrix. The gravity and wind resistance are described as follows:

Position and Attitude Control
The input was the rotation angle of the four nozzles, and the output was the three-axis position and attitude, which comprised an underactuated system. Since the output of the model was more than the input, not all outputs could be controlled independently. The aircraft was controlled by a cascade control structure, and the outermost control loop was the position controller. This control method was called time-scale separation, which worked as long as the inner loop was much faster than the outer loop. The position loop and attitude loop were controlled by cascade PIDs.
The position error is defined as the difference between the desired and estimated position ξ err = ξ des − ξ. The position control law is given by Equation (12), where k p,p , k i,p , and k d,p are the proportional, derivative, and integral gains, respectively. The desired attitude angles are given by Equations (13) and (14). a des = k p,pξerr + k i,p ξ err dt + k d,pξerr (12) φ des = −arctan a des (1) × sinψ + a des (2) × cosψ g (13) θ des = arctan a des (1) × cosψ + a des (2) × sinψ g The attitude error is defined as the difference between the desired and estimated attitude angles η err = η des − η. The attitude control law is given by Equation (15), where k p,a , k i,a , and k d,a are the proportional, derivative, and integral gains, respectively. w des = k p,aηerr + k i,a η err dt + k d,aηerr

Simulation
Our proposed attitude control strategy is verified with a simulation in this section. MATLAB and Flight-Gear co-simulations were conducted; the dynamic model and controller of the aircraft were built in MATLAB, the route navigation was performed by Flight-Gear, and the control frequency was 100 Hz, which did not exceed the PX4 control frequency. The aircraft parameters are specified in Table 5. Trajectory tracking is the first simulation purpose to be examined when there are uncertainties in a system. The results are presented for the designed controllers to display the effectiveness of the new control form. The step signals at 0-5 m along the Z-, X-, and Y-axes in the inertial reference frame were given in Simulink, and then the aircraft flew back over the takeoff point and landed on the origin of the inertial reference. Then, the desired coordinate point at the next moment was given by the Flight-Gear software; the "Earth-North-Up"(ENU) coordinate system was adopted, and the desired trajectory and the actual trajectory are shown in Figure 10a. As shown in Figure 10b, the position tracking errors were within 0.1 m, and Figure 11 shows the velocity responses. There was an approximately 10% overshoot of velocity. The position loop was affected by the attitude loop, which was the inner loop of the position loop. Figure 12 shows that the rolling angle of the aircraft was basically zero, and there was only a slight fluctuation during the disturbance. As a change in the pitch angle corresponded to a change in the aircraft speed and a change in the yaw angles corresponded to a change in the flight path, the aircraft changed course by yawing and used pitch to control flight speed. Since the actuators needed a certain amount of time to respond, there was a phase delay of 0.2 s. The simulation results showed that the aircraft had a good position and attitude tracking ability, thus proving the effectiveness of the designed control method and controller. The change process of the nozzles' angles during the flight simulation is recorded in Figure 13. By comparing the change in the nozzle angles with the flight intentions, we found that: there existed a hover nozzle angle δ g for which the resultant force along the vertical direction was equal to gravity, and the aircraft hovered when the four nozzles' angles were equal to δ g ; the aircraft had upward acceleration when the angles of the nozzles were smaller than δ g , and the aircraft had downward acceleration when the angles of the nozzles were greater than δ g ; the aircraft had positive yaw acceleration when the angles of nozzle 1 and nozzle 3 were smaller than δ g and the angles of nozzle 2 and nozzle 4 were greater than δ g ; the aircraft had negative pitch acceleration when the angles of nozzle 1 and nozzle 2 were smaller than δ g and the angles of nozzle 3 and nozzle 4 were greater than δ g ; the aircraft had positive roll acceleration when the angles of nozzle 2 and nozzle 3 were smaller than δ g and the angles of nozzle 1 and nozzle 4 were greater than δ g . These changes were consistent with the control schematic shown in Figure 8. As mentioned above, 4 ∑ i=1 F i2 = 0 when the aircraft pitched. We recorded the resultant force along the X-axis of the aircraft that was generated by the four turbojet engines during the simulation, as shown in Figure 14. We found that the resultant force was less than 10 N. As discussed in Section 4.1, the angles of nozzle 1 and nozzle 2 were different from those of nozzle 3 and nozzle 4 when the aircraft pitched, so the resultant force along the aircraft X-axis was not zero. The angle difference between the nozzles was larger due to the larger pitch angle acceleration, which caused a larger resultant force along the aircraft's X-axis.

Experiment
In this section, trajectory tracking experiments on an outdoor aircraft are presented in order to illustrate the effectiveness of the proposed control schematic. The experimental prototype is shown in Figure 1. The position and attitude information of the aircraft were obtained by sensors integrated into pixhawk. The information collected by the sensors was filtered and fused by using EKF2 and then fed to the pose controller. The "North-East-Down" (NED) coordinate system was adopted; Figure 15 shows the desired and actual three-dimensional trajectories, and Figure 16 displays the aircraft attitude changes in sequence during the turn marked in Figure 15. During the flight, the turbojet speed was fixed at 90,000 rpm, corresponding to 30 kgf at 15°C, and 25 m above sea level. The onboard fuel capacity was 20 L, and the endurance was about five minutes. The outdoor temperature of the test environment was 20°C and the wind speed was around 0-1 m/s, 50 m above sea level. The initial mass of the aircraft was 61 kg, and its dimensions were consistent with those in the simulation.

Start point
End point Ⅰ Ⅲ Ⅳ Ⅱ Figure 15. Three-dimensional trajectories of the aircraft. The simulation proved that the aircraft could turn by using the yaw; in order to ensure the safety of the flight process and show the aircraft had a wide range of adaptability, as shown in Figure 17a, we kept the yaw angle unchanged during the course of flight and relied on tilting the fuselage to change direction. With reference to Figures 16 and 17, the turning process shown in Figure 16 corresponded to the period of 17-20 s in Figure 17, where the pitch angle and roll angle changed in coordination to complete the turning process. The direction of forward flight was changed at an obtuse angle during the turning phase; thus, the aircraft completed the turning operation with a large roll angle and a small pitch angle. Corresponding to Figure 18,ẋ changed from −2 to 2 m/s, and in Figure 17, the direction of movement changed along the X-axis during the period of 17-20 s. In Figure 17a, there are five major changes in the pitch angle and roll angle. The first and last changes corresponded to the acceleration from the start point and deceleration before reaching the endpoint, and the middle three changes corresponded to the three turns in Figure 15. Aircraft can steer without yaw, and there are more options when the terrain is limited or a load is carried.
As shown in Figure 17, the aircraft had a negative pitch angular acceleration and positive roll angular acceleration between 0 and 1 s. The aircraft headed down and rolled to the right, and it flew from the starting point to the right of the front, which was in the negative direction on the X-axis negative and the positive direction on the Y-axis in the inertial reference frame. In Figure 19, the angles of nozzle 1 and nozzle 2 increased, while the angles of nozzle 3 and nozzle 4 decreased between 0 and 1 s. To generate negative pitch angular acceleration, the angles of nozzle 1 and nozzle 2 were increased and the angles of nozzle 3 and nozzle 4 were decreased; to generate positive roll angular acceleration, the angles of nozzle 1 and nozzle 4 were increased and the angles of nozzle 2 and nozzle 3 were decreased. As the aircraft had the two angular accelerations at the same time, the nozzle angles were superimposed, so we found that the angles of nozzle 1 and nozzle 3 changed very much and the angles of nozzle 2 and nozzle 4 changed little. This is consistent with the results of the analysis in Section 5. The one-dimensional thrust vectoring nozzle modules no longer relied on the response speed of the turbojet engine, and the prototype completed the trajectory tracking experiment successfully. The tracking results in three-dimensional space are shown in Figure 20a; there was a position tracking error of approximately 0.2 m, as shown in Figure 20b. Figure 18 shows that there was a maximum velocity overshoot of 15%. There was a phase delay of 0.25 s, as shown in Figure 17. The attitude was stable during flight, and the aircraft followed the desired trajectory well. Trajectory tracking was realized through the experiment and the simulation, which proved that the proposed control schematic is effective for turbojet aircraft. Compared with the simulation, the error range was larger for the experimental results. First, limitations due to sensor accuracy and measurement and estimation biases were present. Second, the combustion process was uncontrolled, and the thrusts of the four engines were inconsistent. Third, air resistance and gust disturbance occurred during flight. Finally, the center of gravity and inertia of the aircraft changed as fuel was consumed. The values of the angle of the nozzle were estimated with the pulse-width modulation (PWM) output of the steering gear because of the presence of the parallel four-bar mechanism. The changes in the angles of the nozzles are shown in Figure 19. As time progressed, the angles of the four nozzles increased gradually as fuel consumption reduced the weight of the aircraft. The angle of nozzle 1 was obviously larger than those of the other nozzles, which was probably because the aircraft was not completely symmetrical, or because there were deviations in the initial installation angles of the steering gear.

Conclusions
In this study, we demonstrated how an aircraft with one-dimensional thrust vectoring nozzle modules can perform vertical takeoffs and landings. According to the results of simulations and experiments, the one-dimensional thrust vectoring nozzle modules could achieve 87.56% recovery of the turbojet thrust, and the thrust along the engine axis remained unchanged during the rotation of the nozzles. We controlled the position and attitude of the aircraft through the combined movements of the one-dimensional thrust vectoring nozzle modules. The dynamic modeling of the aircraft in the simulations verified the effectiveness of our strategy. We built a prototype for trajectory tracking experiments, and the experimental results showed that the prototype performed well by following the designed control method, which included a good ability to control attitude and position. In addition, the prototype no longer relied on the response speed of the turbojet engine, and it also had the ability to compensate for continuous weight change and uncertain interference.
Throughout the development process, we encountered a series of problems that will form the basis of our future study for developing this platform. Next, we will continue to optimize the aircraft on the basis of this platform to take full advantage of the large thrust-to-weight ratio of the turbojet engine and to improve the load capacity and control accuracy of the aircraft.